直线电机本体建模.docx
《直线电机本体建模.docx》由会员分享,可在线阅读,更多相关《直线电机本体建模.docx(36页珍藏版)》请在冰豆网上搜索。
直线电机本体建模
基于Simulink的直线电机本体建模
电磁发射课题组
2015年10月29日
1直线感应电动机的等效电路
直线电机在结构上可看作是沿径向剖开并将圆周展为直线的旋转电机,如图1所示。
直线感应电动机的稳态特性近似计算方法基本可以沿用旋转感应电动机的等效电路[1]
定子
图1旋转电机演变为直线电机示意图
对于旋转异步电机而言,与电机绕组交链的磁通主要有两类:
类是穿过气隙的相间互感磁通;另一类是只与一相绕组交链而不穿过气隙的漏磁通,前者是主要的。
定子各相漏磁通所对应的电感称作定子漏感Lls,由于绕组的对称性,各相漏感值均相等。
同样,转子各相漏磁通则对应于转子漏感Llr。
对于每一相绕组来说,它所交链的磁通是互感磁通和漏感磁通之和,因此,定子各相自感为:
LaA二Lbb二Lee二Lms'Lls
(1)
转子各相自感为:
Laa—Lbb一Lcc一LmrLl^_LmsLlr
(2)
两相绕组之间只有互感,互感又分为两类:
1)定子三相彼此之间和转子三相彼此之间位置都是固定的,故互感为常值;
2)定子任一相与转子任一相之间的位置是变化的,互感是角位移的函数。
由于三相绕组轴线彼此在空间的相位差为-120,因此互感为:
定转子绕组间的互感由于相互间的位置的变化,为:
LAb二LbA=Lbc=Lcb=Lac=LC^=LmsCOs二'120(7)
LAc二Lca二LBa二La^Lb^LC^LmscO^~120(8)
以上是针对旋转异步电机的参数的推到过程,而对于直线电机,
文献[3]中作者给出了圆筒形直线感应电动机的等效电路,如所示:
图2圆筒形直线电机的等效电路
图2中,Rs和Xs分别代表初级绕组的电阻和漏抗;Rm代表励磁电阻;Xm代表励磁电抗;r20代表次级表面电阻;X20代表次级表面电抗;Red代表边端效应影响纵向边电功率产生的损耗折算成的等效电
阻;Rr代表在次级铜层中的折算的电阻值。
在该文献[3]中次级使用的是导电层和导磁层所构成的复合材料至于图2中的相关参数的计算过程,在该文献中都有详细的说明,不再赘述。
文献[1]中给出了计及边端效应的等效电路,如所示:
图3计及边端效应的等效电路
图3中,bo为励磁电纳(Q);ri为初级绕组电阻;xi为初级绕组漏电抗;a为次级导体电阻折算到初级的换算值,Re为边端效应消耗功率的等效电阻折算到初级的换算值。
2直线感应电机的数学模型
(1)电压方程
参看海军工程大学鲁军勇在文献[4]中给出的电压方程,即:
S-Rsids+Mds-EVYqs
0二RJdr+D屮dr—3(Ve—V^qr''
O=Rriqr+D%r+0(Ve—V^dr
式中:
Rs为通电段定子绕组电阻;尺为通电段定子绕组电阻;乂
为同步速度,V为动子实际速度;D为微分算子;1=二。
注释:
对上式进行简要的推导:
利用三相静止坐标系到两相任意旋转坐标系间的转换矩阵
C3s2r
可将三相静止坐标系下的定子电压方程转换到任意旋转坐标系
dqO
坐标系下,即:
UaRs
Ub
0
Rs
0
C3s/2r
Ua
Ub
-C3s/2r
0iA0Rs」Lic」
Rs00Rs
iB
d
+——
dt
a
屮B
严Cj
0iA
0
iB
%00
Uds_idsd'dsuqs
3s/2r3s/2r
/■ds
H」
=RsI
+
-I
Llf
■
」qs_
dtVqs一
〔-餌0一
dsI
RiC
01
Ud^RsidsDds-‘TqsUqs=RsiqsDq^■■/ds
(10)
对于转子电压方程的推导过程类似,
只是转子坐标系转换矩阵与
定子坐标系的转换矩阵不一样,即:
cost
C3s32r
、2兀
cos(4)
3
_2兀
sin(片)
3
1
2
r2兀1COS^r-——)
3
2兀、
sin(亠-〒)
3
1
、2
利用该转换矩阵将转子电压方程由三相静止坐标系转换到两相
任意旋转坐标系下,即:
Ub
0
Rr
0
_Uc
1
一
一0
0
Rr
=C3s/2r
Ua
=C3s/2r
—I
IJ
=I_Uqr
Udr
Uqr
Rr0
0Rr
'00
0ia
ii.i+dr%i(
一R1一
+—Lr
_1
dt丫qr一
--国s0
qr_
sdr
0
R^JjcJ
0
Ub|2c」
Udr—idrddr
二RridrDdr-sqr二RJqr-D*甲…
2c
dt
C/dr
3s/2rC3s/2rL1.
丿Fqr一
(12)
综上,将电压方程归结为:
(13)
考虑角速度与速度间的关系,即:
R=V—
(14)
—
~R=
将式(14)带入到式(13)中可得:
、卜\、、八
注意:
31
而式(15)中的速度M
v转角速度(机械量),折算关系是:
,是同步速度(定子磁场的速度),V是动子实际的运动速度(机械运动速度)。
(2)磁链方程
鲁军勇在文献[4]中给出的直线电机的磁链方程为:
'ds二LisLui-LmidsLmidr
屮=fL+L+L\i+Li
qs\—Isu1—m/・qs■—m■qr
dr二LmidsL.Lmidr
tqr=Lmiqs•L|rLmiqr
(3)电磁推力方程
文献[4]中给出的直线电机电磁推力方程为:
3.L
Fe二刁十driqs-q」ds(17)
r
注释:
对上式(17)进行简要的推导:
从电磁功率的角度入手,贝心
FeVr
Te
JI
因此,电磁推力与电磁转矩的关系为:
np
T«
em
巧nen
m
(18)
(19)
而我们知道对于旋转异步电机而言,其电磁转矩的表达式为:
Te-npLm(is:
・-isjr|.:
‘)(20)
结合磁链方程将式(20)中的转子电流分量消掉,则:
即"LmQ+Lr。
Lr
r皆Lmis:
Lr
将式(21)带入到式(20)中可得:
TenpLm(isqirdisdirq)
疑问:
式(17)中的系数如何理解?
?
?
?
?
?
分析中采用的都是恒功率转换矩阵,而在鲁军勇的文献中所使用的转换矩阵是恒幅值转换矩阵,下面我们验证这种猜测:
由文献[7]可知恒功率转换矩阵C3s2r和C3s2r分别为:
。
O1
cos八120cos二120-sin二-120-sin一120(23)
11
V2近一
-
COST
2I
C3s2rsin,
3
1
-2
恒幅值转换矩阵为:
cosv-120cos二120
-sin二-120-sin二120(25)
11
22一
cos日-sin日11
C3,=cos(e—120)-sin(e—120)1(26)
cosn120-sin^1201
仍然借助旋转异步电机的电磁转矩来推导电磁推力,将式(25)
和(26)带入到文献[7]中给出的电磁转矩表达式中,即:
Te=npLms[i」aipibicicsi"i』ipici」asin(,120)
(27)i」ciBiai」bsin(—120)]
利用恒幅值转换矩阵将ABC坐标系上的定、转子电流转换到dq0坐标系,即:
iAcosr-sinr1isd
L卜|cos(日1-120)-sin(日1-120)1|卜]ijlcosW+120)-sin(巧+120)1血」
iacost-sins1ird
ib|=|cos(TrT20”)-sin®T20)1jirq
_ic_cos3120-sin1201丄°
由于推导过程相当复杂,但是我们发现在文献[7]中作者指出:
在化简过程中的零轴分量完全抵消了,所以对比两种情况的转换矩阵,
可做如下的推导:
当使用恒功率转换矩阵时:
Te=npLms[i」aTbLDeSi"iJbDa120)
iAiciBiaiCibsin^-120)](28)
3丿
=nL
pms
当使用恒幅值转换矩阵时:
Te二npLms[i」aipLiJcSixi」bipiciciasin^120)
iAiciBiaicibSin(「120)](29)
二npLmsK
因此,采用恒幅值转换矩阵运算时的电磁转矩为采用恒功率转换
3_
矩阵运算时电磁转矩的1.5倍,即T^2n丄m(isqird_isdirq)。
将其带入到
式(19)中可得电磁推力为:
.T3-3-I
F—益”5也pi…-i亠)(30)
(4)运动方程
文献[4]中给出的直线电机在发射阶段的运动方程为:
Mm一=Fe-Bv?
-丄:
;:
Mmgdt
(31)
式中:
风摩系数;
M为负载质量;m为动子本体质量;Fe为电磁推力;B为"为滑动摩擦系数。
注释:
个人认为如果按照式(23)来编写状态方程时,较难列写出速度
的状态方程,因为我们知道状态方程的形式为:
炸Ax,Bu,考虑是
否能将运动方程简化为这种容易列写状态方程的形式,为此,参看文
献[5][6]中给出的运动方程的形式,即:
MmFe-斤-BvV(32)
dt
式中:
Fl—负载阻力;
V—机械运动速度;
Bv—与速度有关的阻尼系数;
将电磁推力的表达式带入到式(32)中,得:
文献⑹中作者将粘滞阻尼系数取:
Bv=°.2Nsm
3状态方程推导
状态方程是指刻画系统输入和状态关系的表达式。
状态向量所满足的向量常微分方程称为控制系统的状态方程,状态方程是控制系统数学模型的重要组成部分。
对于线性系统而言,我们知道其状态方程的形式为:
X=AtxBtu<(33)
y=CtxDtu
状态变量的选取:
直线电机作为异步电机的一种,同样具有4阶电压方程和1阶运动方程,因此其状态方程也应该是5阶的,因此必须选取5个状态变量[7]。
在旋转异步电机中可选的变量共有9个,即转速「、4个电流
变量isd、isq、ird、irq和4个磁链变量’sd、一sq、‘rd、一rq。
个人认为针对直线电机而言,将其中的转速换成速度V,另外,转子电流ird、
爲是不可测的,因此不宜作为状态变量,故只能选择定子电流匚和isq,另两个状态变量必须是转子磁链'rd、‘rq,或定子磁链'sd、‘sq。
为了推导出状态方程,需要结合电压方程(15)和磁链方程,现将两个方程重新列出:
电压方程:
磁链方程:
式(34)中:
rd、1rq•
将磁链方程代入到电压方程中,消去其中的叽、%、ird、i
Usq-RisqP(Ls「sq+Lmir^)"*"3dqs(Ls「sd十Lmird)
1.
0=尺厂(叽-Lmisd
Lr
1
0=R5q-Lmisq
ILr
)+PWrd_(3dqs_3)5q
)+PWrq+(0qs-3)Wrd
由上式的第3、4式可得:
Wrd二Rrisd
dtLr
dLm-
"rq=Rrisq
dtqLrq
c1
_RrL"rd3dqs-3%
Lr
c1
-Rr厂"rq-3dqs-3"rd
Lr
g.dtisd
dtisq
RsL2+RL:
j+i
2isd3dqsisq(i_sL
RsL2RLmi
isq
-3dqsisd
5叽oLsLrT吆
Lm3
Lm3
Lm
_Usd_
cLs
其中:
EsLrTr"rq
Usq
◊——电机漏磁系数,
严1-丄d
LsLr
Tr――转子电磁时间常数,
Tr
Lr
R
•'dqs「’1
将上述推导出的状态方程写成矩阵的形式,则:
RLr+RLmLmLm3
2
(LsLr
COdqs
_isd1disqdt%
3dqs
RsL;RrL:
(LsL;
Lm
Tr
(LsLrTr
Lm3
也Lr
_1_
一Tr
(LsLr
Lm
(EsLrTr
■iSdl
isq
"rd
艸q一
Lm
Tr
-(3qs—3
Tr
■1
ol_s
1
(s
0
0
上式与式(
33)
B
v亠
—v
M总LrM
i」i-
dr1qsqr'ds
旦组
M总
心、
成直线电机的状态方程。
从上述的状态方程中可知,状态变量为:
输入变量为:
-_T
u=Usd,Usq^i,FJ(37)
如果在推导状态方程时使用的是鲁军勇文献[4]中给出的电压方程和磁链方程,则只需对上述的状态方程做如下的修改:
4S-Function的编写
4.1S函数的原理
Simulink中的大部分模块都具有一个输入向量u、一个输出向量y和一个状态向量x,如所示:
引入S函数的引入S函数的目的是为了使Simulink有能力构作一般的仿真框图,去处理如下各种系统的仿真:
连续系统、离散系统、离散和连续混合系统等。
通常S函数的调用格式为:
其中,sfuntmpl为模型文件名,t,x,u分别为当前时间、状态向量,输入向量,而变量flag的值是仿真过程中的状态标志(用它来判断当前是初始化还是运行等),sys输出根据flag的不同而不同,x0是状
态变量的初始化,str是保留参数,
4.2S函数的m文件
K:
\3\asynchronousmotor\ays_m_3.m
function[sys,xO,str,ts]=asy_m(t,x,u,flag,J,np,Rs,Rr,Ls,Lr,Lm)
switchflag
case0
[sys,xO,str,ts]=mdllnitializeSizes;
case1
sys=mdlDerivatives(t,x,u);
case2
sys=mdlUpdata(t,x,u);
case3
sys=mdlOutputs(t,x,u);
case9
sys=mdlTerminate(t,x,u);
otherwise
end
function[sys,x0,str,ts]=mdlInitializeSizessizes=simsizes;
x0=[0,0,0,0,0];%初始化
str=[];%固有格式,预留的
ts=[00];
functionsys=mdlDerivatives(t,x,u)
Rs=6.33;Rr=32.45;Lm=0.06212;Lr=0.08;Ls=0.125;rou=1-Lm*Lm/(Ls*Lr);Tr=Lr/Rr;np=2;J=0.002;Bv=0.02;m=10;M=10;miu=0.05;g=9.8;tao=0.0616;
sys
(1)=-(Rs*Lr*Lr+Rr*Lm*Lm)/(rou*Ls*Lr*Lr)*x
(1)+u(3)*x
(2)+Lm/(rou*Ls*Lr*Tr)*x(3)+Lm*x(5)*pi/tao*x(4)/(rou*Ls*Lr)+u
(1)/(rou*Ls);
sys
(2)=-u(3)*x
(1)-(Rs*Lr*Lr+Rr*Lm*Lm)/(rou*Ls*Lr*Lr)*x
(2)-Lm*x(5)*pi/tao*x(3)/(rou*Ls*Lr)+Lm*x(4)/(rou*Ls*Lr*Tr)+u
(2)/(rou*Ls);
sys(3)=Lm*x
(1)/Tr-x(3)/Tr+x(4)*(u(3)-x(5)*pi/tao);sys(4)=x
(2)*Lm/Tr-x(3)*(u(3)-x(5)*pi/tao)-x(4)/Tr;
sys(5)=-Bv*x(5)/(m+M)+pi/tao*Lm*(x(3)*x
(2)-x(4)*x
(1))/(Lr*(m+M))-FL/(
M+m);
functionsys=mdlOutputs(~,x,~)sys=[x
(1);x
(2);x(3);x(4);x(5)];
functionsys=mdlUpdata(t,x,u)
sys=[];
functionsys=mdlTerminate(t,x,u)
sys=[];
5直线电机的矢量控制仿真
5.1仿真模型
直线电机作为异步电机的一种,其动态数学模型同样是一个高阶、非线性、强耦合的多变量系统。
虽然通过坐标变换可以使之降阶并简化,但是并没有改变其非线性、多变量的本质[7]。
因此,仍然需要采用相应的解耦控制策略来实现直线电机调速系统的高动态性能。
目
前应用最多的方案有:
1)按转子磁链定向的矢量控制系统;
2)按定子磁链控制的直接转矩控制系统;
将仿真模型分成如下几个部分分别单独介绍:
(1)转子磁链观测模型
仿真模型是按转子磁场定向的,因此转子磁链位置的精确观测是控制系统能否实现定子电流转矩分量和励磁分量成功解耦的关键。
转
子磁链位置的表达式为:
"丨F亠心12dt(38)
在文献[9]中动子磁链角为:
\'sdt(39)
式中r为动子运动的电角频率。
注释:
个人对式(38)和(39)的理解是:
'r应该是而不是电角频率而不是动子运动的实测角频率。
正因为如此,在实际建模的时候,动子磁链角米用的表达式才改写为:
(40)
其中:
np为感应电机极对数。
当采用转子磁链定向时,满足:
L「,R为动子时间常数
实际上,此时的电磁推力的表达式也相应的变为:
此时的转差角速度满足:
式(40)、(41)和(43)构成转子磁链观测模块的主要方程。
从式(43)中可知:
转子磁链处在分母的位置上,因此在电机启
动的时候,转子磁链为零,则式(43)出现奇异点发散现象,造成仿
真错误,因此,可将上式改写为:
Lmi
sq
T/rk
可将式(44)中的k取为一个很小的常数
注释:
这样的做法应该是相当于给转子磁链设定了一个初始值,如果我们不采用这种做法,而采用预励磁控制,是否也可以起到这样的作用,有待验证。
(2)励磁电流和转矩电流计算模块
励磁电流和转矩电流的计算表达式为:
1TrP.
内部模块为:
2*u
(2)*Lr*tao/(3*pi*Lm*u
(1))
Phir
v-
V
图5励磁电流和转矩电流计算模块
(3)电流调节器
由于控制方案中采用的是SVPWM调制方法,其输入量则必须为电压,因此需要将电流转矩分量和励磁分量转换为相应的电压分量,因此需要推导出其数学表达式。
在文献[9]中给出的矢量控制条件下直线感应电机的定子电压方
程的表达式为:
-Lmisq
Lrirq
irq
■m
—Lmisd
Lrird
屮d
rd
Lm.
i
Lr
—Lmisd
sq
=Uds=RsidsL
二RJdsL
Rsids
Uqs
I屮
rq
屮
二rd
disd
dT
disd
Ls
=Rids二L
二RsiqsPsq
-Lmisq
Lrirq
=Uqs
Lmisd
Lrird
二RsiqsL
二RsiqsL
=Rsiqs+L
pl:
Lm±-Lisq-「Lmirqdt
Lmdisd
dtLrdt
__LmdL
LrdtLrdt
%Ldr
dtLr
irq
L2
1LsisqFisq
LrdtLr
isq
Lm■
sqr
一Li
rdm■sd
(47)
L2
Lm
s
Lr丿
disq,
卫L
dtmdt
disqLmdisq
dtLrdt
L2!
disq
Lrdt
disqLd®
丝—•一丄sisd
dtLr
dirq
■m
■m
sq
dt
1Ljsd「Lmird
丄sisd
Lm;/
1Lr
r
因此,按照上述的推导得:
Uds
二Rsids
Uqs二Rsiqs
+©・屮
1Lr
L2
i
1sd
r
L2
ImI■
「匸I
sd
disdLmd'r..
竺m匚…•.Ii
1—s■sq
dtLrdt
disqLd屮
聖m「:
;’:
*应丨i
1fd
dtLrdt
s'sd
(48)
在式(46)中作者指出:
该式是在保持定子励磁分量isd不变的情况下,而观察式(49)可知,还需默认转子磁链恒定,才能得到式(46),但是貌似该式中没有系数匚,是不是作者推导错误。
通常采用PI调节的方式来求解电压分量,但是文献[9]中给出的式(46)可知,需要采用前馈解耦控制测量对交叉耦合电势进行解耦处理。
其内部模块为:
isd*
图6电流调节器模块
(4)SVPWM模块
(5)LIM模块
在上文中利用S函数建立了直线电机的数学模型,将其封装成模块,并建立如所示的LIM模块。
图7LIM模块
(6)
矢量控制模型
初步搭建的仿真模型为:
图8矢量控制框图
该模型是初步建立,需要进行如下的修改:
A、转子磁链闭环修改成磁链开环
因为计算转子磁链要受到电机本身参数的影响,当这些参数发生变化时,造成磁链估算不准确,但是注意到在转子磁链角即电磁推力
计算中都要用到转子磁链,即:
B、求解定子电压dq轴分量时是否采用前馈补偿
在图8中求解定子电压dq轴参量参考值时直接采用了PI调节,未使用前馈补偿,下一步需进行改进。
(7)
整体仿真模型
Subsysteml
ia
isd
ib
isq
ic
theta
图9整体系统控制框图
5.2仿真结果
0
斗
5
6
78910
图10三相电压波形
图11电磁推力波形
图12动子运动速度波形
参考文献:
[1]叶云岳.直线电机原理与应用[M].