一种直线感应电机的参数计算方法与流程

1.本发明涉及机电领域，尤其是一种直线感应电机的参数计算方法。


背景技术：

2.随着直线驱动技术的不断发展，直线电机已经越来越多地应用在轨道交通、电磁快速驱动、航空航天、垂向牵引等领域中。直线感应电机在初级绕组中通入电流，产生的磁动势通过气隙使次级导体板上形成感应涡流，从而直接获得直线推力。直线感应电机继承了异步电机结构简单、制造方便、运行可靠和成本较低等优点，成为了直线驱动领域的主要装置。
3.然而，直线感应电机也存在着参数获取难、调速性能不高等问题，难以在宽速度范围，如0～马赫级的宽速域下精确运行。直线感应电机由其自身结构引起的磁路开断等问题，会带来纵向和横向端部效应等现象，同时还存在铁心饱和等非线性特性，这些都将直接影响电机的电磁参数，尤其是应用于持续加速的直线感应电机，其励磁电感和次级电阻需要在不同速度下进行修正。
4.对于处在设计阶段的直线电机，常用类似于旋转电机的计算公式通过电机结构尺寸的设计参数获得相应的电磁参数。但是，除了计算方法比较复杂，直线电机在工程中的各种特殊结构也使得通过计算公式难以获取准确的结果。此外，直线电机的端部效应也不易通过补偿系数应用到工程实际中。对于直线感应电机实物，可以采用类似旋转电机的堵动和空载试验来获得电机的电磁参数，但是对直线电机进行正式试验十分困难。
5.文献“何晋伟,史黎明.一种基于静态特性的直线感应电机参数辨识方法[j].电工电能新技术,2009,28(04):50-53 70.”中描述了一种基于变频器驱动的直线感应电机等效电路参数辨识方法，利用不同频率时的电机推力、变频器输出电压和电流等静态堵转特性得到电机参数。但该方法是建立在静态参数的基础上，适用的速度范围较小，也没有考虑铁心饱和对电磁参数的影响。
[0006]
文献“徐伟,孙广生,李耀华,金能强.一种新型的直线异步电机参数测量方案[j].电工技术学报,2007(06):54-58.”中描述了一种通过开路试验的方法测定单边直线感应电机的电磁参数的方案，并考虑了铁损和附加损耗。但是，该方法需要建立直线感应电机测试平台，通入实际的电压电流，将直线推力转换为转速进行测量，对于工作在高速、大电流和大推力的直线感应电机是无法实现的。
[0007]
中国专利“李健,贺昆.一种直线感应电机参数离线测量方法[p].湖北：cn111257749a, 2020-06-09.”公布了一种通过实验的方式离线测量直线感应电机参数的方法。该方法对直线感应电机进行额定频率电流的短路实验和空载试验，确定被测电机的短路和空载特性曲线，从而获得电机参数。但是，该方法需要将被测电机的次级拖动至所需的速度，不适用于大型电机和高速电机。
[0008]
中国专利“吕刚,马云双,孙守光.用于通用变频器的直线感应牵引电机参数自动识别方法[p].山东：cn102193065a,2011-09-21.”公布了一种利用变频器应用矢量控制算
法测量直线感应电机参数的方法。该方法通过直流实验确定电机的矢量模型，再通过空载实验测定电机的参数。但是该方法未考虑直线感应电机的非线性特性，也没有考虑次级对电机电磁特性的影响，不适用于需要精确控制的直线感应电机。
[0009]
总结上述现有技术可知，电机试验虽然能获得电机的电气特性曲线，获得范围内的电机参数，但是建立直线电机试验平台十分困难，也无法兼顾不同速度、不同电流、不同结构和不同型式的直线感应电机，不能适合于高速大推力的工程场合。


技术实现要素：

[0010]
为了简洁准确的获得直线感应电机的电磁参数，本发明提出一种直线感应电机的参数计算方法，将直线感应电机的端部效应、铁心饱和以及交流损耗等影响定量体现，特别适合在短时间应用于宽速域、长距离、大功率直线驱动的直线感应电机，适用于不同结构型式的直线感应电机。
[0011]
本发明的技术方案为一种直线感应电机的参数计算方法，包括如下步骤：
[0012]
(1)建立直线感应电机的三维瞬态电磁场有限元模型，将直线电机初级供电单元的铁心轭部尺寸、齿槽尺寸、极距极数、绕组形式、绕组电阻率、初级和次级之间的气隙尺寸、次级反应板的结构型式和尺寸参数输入到三维瞬态电磁场有限元模型中；
[0013]
(2)将铁心的磁化曲线、次级反应板的电阻率、初级绕组电流电磁参数输入到三维瞬态电磁场有限元模型中；
[0014]
(3)设次级反应板的运动速度为0，此时电流的角频率ω1即为转差角频率ωs，设置电机初级绕组激励为有效值相同、频率不同的电流，对直线感应电机的有限元模型进行计算，获得不同频率的电流对应的电机推力f，得到一组推力f和转差角频率ωs一一对应的计算结果；
[0015]
(4)根据如下公式，建立以转差角频率ωs为自变量、电机推力f为因变量、次级电阻 r'2和励磁电感l
m0
为参数的非线性拟合模型，将步骤(3)中的有限元模型计算得到的f-ωs输入到非线性拟合模型中，采用levenberg-marquardt算法，设置迭代停止判据为简化的卡方检定值(reduced chi-sqr)，求取以转差角频率ωs为自变量、电机推力f为因变量的拟合曲线以及对应的次级电阻r'2和励磁电感l
m0
的拟合值；
[0016][0017]
式中，m1为初级绕组的相数，i1为初级相电流的有效值，τ为极距，ωs为转差角频率，r'2为次级电阻，l
m0
为励磁电感；
[0018]
(5)以步骤(4)中获得的次级电阻r'2和励磁电感l
m0
的拟合值作为参数的初始值，将步骤(3)中的计算结果分n段输入步骤(4)中的非线性拟合模型中，采用levenberg-marquardt 算法，设置迭代停止判据为简化的卡方检定值(reduced chi-sqr)，求取以转差角频率ωs为自变量、电机推力f为因变量的n段拟合曲线，以及对应的n组次级电阻r'2和励磁电感的拟合值l
m0
；
[0019]
(6)改变初级绕组电流设定值，重复步骤(3)、(4)、(5)，获得不同电流下被测直线感应电机的次级电阻r'2和励磁电感l
m0
；
[0020]
(7)建立直线感应电机的无次级三维瞬态电磁场有限元模型，电机参数与步骤(1)
和 (2)中的模型参数相同；
[0021]
(8)对直线感应电机注入不同频率的电流i1，通过步骤(7)中的无次级三维瞬态电磁场有限元模型计算感应电势，获得电机的初级电阻r1和初级漏感值l1，通过步骤(1)中的模型的得到的感应电势e和注入的不同频率的电流i1，修正有次级情况下的电机的初级漏感值；
[0022]
(9)通过步骤(1)至步骤(8)，获得直线感应电机在不同电流、不同速度下的电磁参数。
[0023]
进一步的，步骤(1)中的计算模型可以包括电机底座、紧固件等影响电磁性能的结构件。
[0024]
进一步的，步骤(1)和(2)中的结构尺寸和电磁参数可以作为变量进行并行计算。
[0025]
进一步的，步骤(4)中的电流频率设定值应集中在对应最大推力值的电流频率值附近。
[0026]
进一步的，步骤(5)中的n为2时，可以将对应最大推力值的电流频率值作为分界点。
[0027]
进一步的，步骤(5)中，可以将参数的初始值替换为分段后获得的次级电阻和励磁电感的拟合值，不断增加n的值，重复步骤(5)获得更准确的结果。
[0028]
进一步的，步骤(8)中的模型的电流频率设定值可以在电机运行工况的全速度范围内取对应的电流频率值。
[0029]
通过本发明所构思的以上技术方案与现有技术相比，能够取得下列有益效果：
[0030]
(1)本发明所提供的直线感应电机的参数计算方法，采用电磁场有限元计算的方式，将电机的铁心齿槽、绕组、气隙、次级反应板的结构型式、尺寸以及铁心饱和特性曲线、次级反应板电阻率等进行参数化，对不同电流给定值和频率进行并行计算，可以获得更高的计算效率，更高的准确性，更广泛的适用性，特别适合于宽速域、长距离、大功率直线驱动的直线感应电机。
[0031]
(2)本发明所提供的直线感应电机的参数计算方法，建立了直线感应电机的三维瞬态电磁场有限元计算模型，将直线电机的纵向和横向端部效应纳入计算，还将电机底座、紧固件等影响电磁性能的结构件包含在计算模型中，可以获得更精确的电机参数。
[0032]
(3)本发明所提供的直线感应电机的参数计算方法，利用直线感应电机的推力计算公式，建立以电流频率为自变量、电机推力为因变量、次级电阻和励磁电感为参数的非线性拟合模型，通过多次分段非线性拟合计算获得次级电阻和励磁电感值，可以获得更准确的电机参数。
[0033]
(4)本发明所提供的直线感应电机的参数计算方法，在不同电流给定值和频率下进行参数计算，将电机在铁心饱和特性曲线上不同运行点的工况以及交流损耗的影响体现在计算结果中，可以获得更准确的电机参数。
[0034]
(5)本发明所提供的直线感应电机的参数计算方法，从工程实际出发，将相同电流输入情况下的有次级模型和无次级模型对比计算，对不同磁场激励下的初级漏感和励磁电感值进行修正，可以获得更准确的电机参数。
附图说明
[0035]
图1为采用本发明的方法计算电机参数的直线感应电机实施例一主要结构模型截面图；
[0036]
图2为实施例一初级铁心的饱和特性曲线；
[0037]
图3为直线感应电机等效电路；
[0038]
图4为实施例一的堵转推力频率特性初步拟合结果与仿真结果对比图；
[0039]
图5为实施例一的堵转推力频率特性分段拟合结果与仿真结果对比图；
[0040]
图6为实施例一的堵转推力电流特性的预计算结果和仿真结果对比图；
[0041]
图7为实施例一的初级电阻的计算结果；
[0042]
其中，1初级绕组、2次级反应板、3初级铁心、4电机固定板。
具体实施方式
[0043]
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整的描述，显然，所描述的实施例仅为本发明的一部分实施例，而不是全部的实施例，基于本发明中的实施例，本领域的普通技术人员在不付出创造性劳动的前提下所获得的所有其他实施例，都属于本发明的保护范围。
[0044]
根据本发明的实施例，提出一种直线感应电机的参数计算方法，步骤如下：
[0045]
(1)图1为实施例一的主要电磁结构的纵截面图，包括初级铁心3、初级绕组1、次级反应板2以及电机固定板4。对实施例一建立三维瞬态电磁场有限元计算模型，将初级铁心 3、初级绕组1、气隙、次级反应板2的结构型式和尺寸作为设计参数输入到计算模型中，将电机底座、紧固件等影响电磁性能的结构件也输入到模型中。计算模型的求解域中涉及到的矢量有：磁场强度h、磁感应强度b、电场强度e、电场密度j、导电媒质在磁场中的运动速度v。模型求解域中电磁场各量的关系为：
[0046][0047]
式中，σ为电导率、μ为磁导率。由洛伦兹方程可得，电磁力密度f为：
[0048]
f＝j
×b[0049]
(2)将实施例一的初级铁心的磁化曲线(图3)、次级反应板的电阻率、初级绕组电流等电磁参数输入到模型中。
[0050]
(3)将实施例一的次级反应板的运动速度在有限元计算模型中设定为0，将初级绕组电流的频率设置为1hz、3hz、5hz、7hz等多个不同的值，进行有限元计算，获得不同电流频率的电压电流，初次计算后可知最大推力点的电流频率在1hz至3hz之间，故进行电流频率为1.5hz、2hz、2.5hz等的补充计算。根据步骤(1)中的公式，通过离散化求取偏微分方程的数值解，求解过程简要归纳为以下几个主要步骤：
[0051]
①
建立以磁位a为变量的偏微分方程组，在给定的区域和边界条件下，求出能量泛函达到极小值时的磁位a；
[0052]
②
把求解区域用三角单元或四边形单元划分为有限的网格，然后在单元中构造一个插值函数；
[0053]
③
将能量泛函的极值问题转化为能量函数的极值问题，建立相应的线性方程组；
[0054]
④
利用第一类边界条件对方程组进行修正，利用高斯消去法或三角分解法，求解上述建立的线性方程组，得到所有节点的磁位，进一步求取其他需要的磁通密度、绕组电感、感应电动势和电磁力。
[0055]
(4)实施例一的直线感应电机等效电路如图4所示，图4中初级电阻r1、初级漏感l1、励磁电感l
m0
、次级电阻r'2为待计算的参数，转差s为计算时需要设置的参数。直线感应电机推力f计算公式为：
[0056][0057]
式中，m1为初级绕组的相数，i1为初级相电流的有效值，τ为极距，ωs为转差角频率。
[0058]
利用直线感应电机的推力计算公式，建立以电流频率(转差角频率ωs)为自变量、电机推力f为因变量、次级电阻r'2和励磁电感l
m0
为参数的非线性拟合模型，输入步骤(3)中的计算结果，采用levenberg-marquardt算法，设置迭代停止判据为简化的卡方检定值 (reduced chi-sqr)，获得以ωs为自变量、f为因变量的拟合曲线以及对应的r'2和l
m0
的拟合值，图5为实施例一的第一次拟合结果与仿真结果的对比图。
[0059]
(5)以步骤(4)中获得的r'2和l
m0
的拟合值作为参数的初始值，将步骤(3)中的计算结果分n段输入步骤(4)中的非线性拟合模型中，采用levenberg-marquardt算法，设置迭代停止判据为简化的卡方检定值(reduced chi-sqr)，获得以ωs为自变量、f为因变量的n段拟合曲线以及对应的n组次级电阻和励磁电感的拟合值，图6为分2段(n＝2)的拟合结果和仿真结果的对比图。
[0060]
(6)图7为按照直线感应电机的计算公式计算的不同i1下的f和仿真结果对比图，图 7说明不同电流下的r'2和l
m0
是变化的，改变i1的设定值，重复步骤(3)、(4)、(5)，获得不同电流下实施例一的r'2和l
m0
。
[0061]
(7)建立实施例一的无次级的三维瞬态电磁场有限元计算模型，电机的设计参数和电磁参数与步骤(1)和(2)中的模型参数相同。
[0062]
(8)对实施例一的有限元计算模型进行电流频率为1hz、10hz、50hz、100hz、200hz 等不同值的有限元计算，通过步骤(7)中的模型的感应电势和电流的计算结果，获得电机的初级电阻r1和初级漏感l1的值，通过步骤(1)中的模型的感应电势和电流的计算结果，修正有次级情况下的电机的l1的值。
[0063]
通过步骤(1)至步骤(8)，可以获得直线感应电机在不同电流、不同速度下的电磁参数，包含端部效应、铁心饱和及交流损耗的影响。图7为计算的各同步角频率下初级电阻r1的值。
[0064]
直线感应电机是短初级、长次级形式，或是长初级、短次级形式。
[0065]
本发明描述的直线感应电机的参数计算方法，适用于短初级、长次级形式直线感应电机，适用于长初级、短次级形式直线感应电机，改变自变量参数后，同样适用于直线同步电机，以及其他形式的直线电机。
[0066]
尽管上面对本发明说明性的具体实施方式进行了描述，以便于本技术领域的技术人员理解本发明，且应该清楚，本发明不限于具体实施方式的范围，对本技术领域的普通技术人员来讲，只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内，这些变化是显而易见的，一切利用本发明构思的发明创造均在保护之列。