永磁同步電機控制原理

1、永磁同步電機控制原理-yyc1永磁同步電機原理1.1永磁同步電機基本結構永磁同步電機轉子是永磁鐵，定子是三柑繞組。 假設1）定子繞組為Y形接法，三相繞組對稱分布。2）忽略定子繞組的齒槽對氣隙磁場分布的影響。3）忽略鐵心損失4）忽略電動機參數(shù)（繞組電阻和繞組電感等）的變化。 三相兩極永磁同步電機的結構：a)b)a）定子繞組分布圖b）定子繞組示意圖1.2永磁同步電機動態(tài)數(shù)學模型ABC三項繞組的磁鏈可以拆分為2部分，定子繞組自感互感產生的磁鏈和轉子永磁體產生的 磁鏈。aa Lab Lac中4mba Lbb Lbc+屮BmCA “B Lgic定子繞組對稱分布S = LS-Lm cos(20)Lbb =

2、 G + Lm cos(2(8 27r/3)Lee = G + Lm cos(2(0 + 2tt/3)Ls為三相繞組自感平均值Lm為三相繞組自感二次諧波幅值B為轉子電角度電角阪=N%N為極對數(shù)叫為機械角度定子繞組互感Lca = Lac = 一M$ Lm cos(2 (% + 石 + 2w/3)注意，這堅1忖定義的是負號,跟有的書上不一樣。轉子永磁體在三相繞組產生的磁鏈與電流無關，僅與電角度冇關屮4m = Wm COS(0)屮Bm =屮mCOS(& -2兀/3)屮Cm=屮m COS(& + 27T/3)永磁同步電機模型的表達關鍵在于磁鏈數(shù)據(jù)，現(xiàn)在磁鏈數(shù)據(jù)表達為4個參數(shù)Lm Ms i|/m 轉矩方

3、程化簡表達式:Matlab 代碼:syns thetasms L_AA L_AB L_AC L_BA L_BB L_BC L_CA L_CB L_CC Ls Lm MsL-AA = Ls +Ln*cos(2*theta);L_BB Ls 4Ln*cos(2*(theta-2*pi/3);L_CC = Ls *Ln*cos(2*(theta*2*pi/3):L_AB = -Ms - Lu*cos(2*(theta + pl/6);L.BC = -Ms - Ln*cos(2*(theta + pi/6 - 2*pi/3):L_CA = -Ms - Lii*cos (2* (theta pi/6 土

4、 2*pi/3):L_BA 二 L-AB;L.CB - L_BC：L_AC = L_CA：syns L.ABC 1_A1_Ci_ABC = Ci_A： i.B： i_C：L_ABC = L_AA L_AB L_AC； L_BA L_BB L_BC： L_CA L_CB L_CC；syns phi_Aa phi_Bm phi_Cx phi_Mphi.Am = phi_M*cos(theta):phl.Bm = phl_Mcoslthea-2*pl/3):phi.Cm = phi_M*cos(.theta*2*pi/3);syns phi.ABCnphl_A3Cn = phi_Am： phi_Bm

5、： phiCmJ:syras Te NTe = N/2*(i_ABC. *J*diff(L.ABC theta)*i_ABC * (i.ABC. )*diff(phi-ABCm, theta): slnpllfy(Te|)1.3坐標變換由前面推導發(fā)現(xiàn)在ABC三相坐標下磁鏈往往都相互耦合，與三相電流，角度都相關，表達式 很復雜，因此想到坐標變換。a) 3s坐標系b) 2s坐標系c) 2r坐標系電動機內的氣隙磁場是進行電磁能量傳遞的媒介，定轉子之間的能量傳遞是靠氣隙磁場 來進行的，不同類型的繞組進行變換時，需要保持它們的總磁動勢不變。嚟=“3 ( +滋片喑=比仏+ 4力）/嚟表示3 s坐標繞組產生

6、的磁動勢在2s坐標中的衣達 珂：表示2s坐標繞組產生的磁動勢在2s坐標中的表達 於:茨示2r坐標繞組產生的磁動勢在2s坐標中的表達為了讓推導出來的三相電流與兩相電流幅值相等，取護=2/31 - 2V3-2_12V3_2同理得到:2宀2$ =cosOsin 3-shi 0cos 6 .C2sf2r =COS0一 sin 8sin 8cosJ1.4 dq坐標系的動態(tài)數(shù)學方程執(zhí)行坐標變換ld1=2s2r3s2slALSKl=22山325Matlab代碼syits thetaC.2s_2r = Leos(thetaj sin (theta): -sin(theta) cos(theta):C_cs_2

7、s = 2/3*L -1/2 -1/2： 0, sqrt(3)Z2 -sqrt /2;C.2s_3s = ：1 0;-l/2 sqrt(3) /2;l/2 -sqrt(3)/2;C_2r_2s = Icos I thetG -sin(theta; sin(theta) cos(theta); i_ABC C_2s_3s*C_2r_2s木i_d;i_Q】：u.ABC = C.2s.3s*C.2r_2s*u_d;u_q;phi.ABC = C_2s_3s*C_2r_2s* lphi_d; phi.q |;i_A = 1_ABC(1);u.A = u_ABC(l);phi_A = ohi_ABC(l

8、):得到：uA = ud cos 0 uq sin 0“ =id cos 3 iq sin 8= id cosO 一$ sin&帶入A相電壓方程:dA得到:iia cos 6 uq sin 0整理得到:才 d + R j COS 01 Uq ” $ CD書q RH I Sill 8 = 0對任惡角度&都滿足，所以:ud = cod + Rlid+ 3書 q + R liq接下來就是求dQ的表達式，并且把它和前面4個參數(shù)0 Lm Ms屮m關聯(lián)起來Matlab代碼syms tlvtasyms L AA L AB L AC L BA L BB L BC L CA L CB L CC Ls Ln Ms

9、L.AA = Ls 45*cos(2theta)：L_BB = Ls LmcosCZ*(theta-2*pi/3):L-CC = Ls Hm*cos(2*(theta*2*pi/3):L_AB = -Ms - Lm*cos(2*(theta * pl/6)：L.BC =-Ms -Lm*cos(2*theta*pi/6-2*pi/3):L_CA -Ms -Lm*cos(2*(theta4pi/62*pi/3):L_BA =L_AB：L_CB =L_BC：L_AC =L_CA：syms L_ABC i_A i_B i_Csyms i_ABC i_d 1_qi_ABC = 1_A: i_B: i_C

10、：L.ABC = L_AA L_AB L_AC： L_BA L_B3 L_BC； L_CA L_CB L_CC： syms phi_Am phl.En phiCm phl.M phl_Ax 二 phl_M*oos (theta), phl_Bx = phlM*cos(theta-2*pl/3); phl_Cx 二 phl_M*cos(theta*2*pl/3);syms phl_ABCm phl.ABC phl.ABCo - phl.An, phl_Bn, pril_Cxl, phL.ABC 二 L-ABC*L-A3OphL.ABCm;C_2s_2r - cos(thete) sin(thet

11、aj, sin(theta) cos (theta)：,C_3s_2s = 2/3*l -1/2 -1/2; 0, sqrt(3)/2 -sqrt(3)/2j;C_2s_3s 三1 0-1/2 sqrt(3)/2; T/2 -sqrt (3)/21;C_2r_2s = Icos Ltheta) sin(theta/; sin(theta) cos(thwta)2, phi.dq = C_2s_2r*C_3s_2s*phi_ABC;i-ABC = C_2s-3s*C_2r_2s*i_d,i_q：;phi.dq = subs(phi_dq, (i.A, i.B i-C. (i.ABC(l), i_

12、ABC(2), i.ABC (3);| simplify (phidq)phi_M +(l_q*(2*Ls -3*Lm + 2*Ms) /2最終得到dq軸磁鏈表達式Ed &已經(jīng)與角度無關了，JI互相解耦了Ld Lq與前面4個參數(shù)的關系如下3Lq = Ls + Ms - -LmLo = Ls - 2Ms最后，得出dq軸下轉矩方程Te = 1.5N陽q - 忖d) = 15N(%4 + 仏 - Ljdij1.5永磁同步電機的穩(wěn)態(tài)運行永磁同步電機穩(wěn)態(tài)運行相量圖。圖中，R1和XI是電樞繞組上的阻抗和漏抗。Xad和Xaq是dq軸電樞反應電抗。 Xd和Xq足dq軸電樞同步電抗。Xd = X+Xad Xq

13、=X“Xaq繞組不通電時轉子產生的感應電動勢叫空我反電動勢E0 矩角定義為E0落后于電樞電壓矢就U的角度.0可以看到U由幾部分組成EO jiXl jIXaq IR1可以得到【dXd + IqR, = U COS 9 Eq-IdRi + lqXq = UsmO得到：rXq(7cose - Eo) -RosineIdR/+XqXd,_ Ri(/cos0 - EQ)-XdU sin 61lR/+XqXd stemPlnsical s stemReferenceGcncratucd rtfIPouer IonxertcrCurrentControllerVelocity(ontrollerPositi

14、on Sensor Signal ProcessingProteclion and auxilian* fiincdonsCurrent2.2永磁同步電機控制策略永磁同步電機控制，策略的不同主要體現(xiàn)在具體Id lq怎么控制，大體上有3類。FOC：定向磁場控制id = 0MTPA：最大化Tq/I弱磁査表：Id和lq通過査表來得到2.2.1 FOC2.2.2 MTPA2.2.3弱磁査表3永磁同步電機建模3.1永磁電機模型主要有2類模型，基于dq軸表達磁鏈和基于ABC坐詠表達磁鏈3.1.1基于dq軸坐標表達磁鏈的永磁同步電機模型1) Matlab Simscape Power SysteOm 工具箱

15、中的 Machines / Permanent Magnet Rotor 模塊LibraryMachi nes / Perma nent Mag net Rotor磁鏈參數(shù)農達3 Block Parameters: Permanent Magnet Synchronous MotorXPennaiieiil Magnet Synchronous MolorThis block represents a penianent Magnet synchronous bio tor with sinusoidal flux distribuiio n.Righl-click on the block

16、and select Sinscape block choices to access variant impleraentations of this block.SettingsIk i nInit ial ConditionsNumber of pole pairs:Pcraancnt magnet flux linkage paia:neterization:PerfMnent wagnet flux linkage:S .fn! | ii meter /it ion:Srator d-axis inductance, Ld:Stator q axis inductance, La：S

17、:ator zero-sequenee induct anre. 1,0:Srator resistance per pliase, Rs:E3|ihqSpecify Ld, Lqt and LOHvgIIV*LO11v|RsOhir7Specify flux linkageUK | CancelHelpApp ly2) Matlab Powertrain Blockset 匸具箱中的 Interior PMSM block 模塊Interior PMSMImplenGnt eynarrwes & an 應nor pormanem magno synchronous motorexpand a

18、l mpkeLibrory: Propulsion / Dectnc Motor?磁鏈參數(shù)表達O Elock Parameters: Interior PMSM SS SpdXInterior PMSM (mask) (link)Model the dynamics of a three-phase interior pern)anent magnet synchronous motor (BISM) vrith sinusoidal back electroti ve force.Block(Ht ionsPort Configuration: SpeedParameters ON Canc

19、el | Heh)它是通過引mulink而不是simscape來實現(xiàn)的.因此可以看到內部實現(xiàn).相電壓數(shù)據(jù)進來進行Clarke變換和Park變換.得到dq電壓.在PMSM Equivalent Circuit子 系統(tǒng)中依據(jù)dq電壓微分方程積分得出dq電流通過轉矩方程計算轉矩，通過逆Park變換和逆Clarke變換計算出相電流。3.1.2基于ABC坐標表達磁鏈的永磁同步電機模型Matlab Simsacpe Electronics 1】具箱中的 The FEM-Parameterized PMSM block 模塊FEM-Parameterized PMSMPorma nont rnqgnot s

20、ynchro nous motor dofinod in terms of magn otic flux lin kageLibraryRotational Actuators它的磁鏈數(shù)據(jù)是通過矩陣來表達的，數(shù)據(jù)可以來源于冇限元仿真數(shù)據(jù)。磁鏈數(shù)據(jù)的表達形式有兒種2D數(shù)據(jù)模型，固定互感可 Block Parameters: FEM-Parameterized PMSMXFEM Piinimctcrizcd PMSM、block lupleoen Meclrical and aechaniral i*1* h i i u i m , wi. iwhI nugfi (motor ln limogji

21、ctic flux linkage depends nonlinearly on staler currents and rotur axiglc. Rigid click on the block and select SinsrapcdM) block choices to icccss Viiriant iBplrnh-ntiit ions.SettingsElectricalIron losses hanitdlParameterization:Current vector, i:Assuhc constant autual inductance - tabulate ith pha-

22、 rent and rotor angleRotor angle vector, theta:Flux liiikuxc partial drriviit ivc wrt currcnl. dPhi(ivthota)/di:Flux linluigc partial derivative rt angle. dPhi(i.theta)/dtheta:Gilculatc torque Matrix?:Inierpolat ion iwthod:|icc Ih xVec thdf luxdiU.it ri xTli. ASlblor rcsisiiuicc per phase Rs:CancelA

23、pplySuior Buluul induetuiicu.OK2D數(shù)據(jù)模型，正弦反電動勢可 Blo|x*fids (kjiiI incarly on stator currents ssesParameieri2ation:Direct-axis current vector id:Qutdruiurc-uxis currcni vector iq:Ld Btitrix, I i. iq):Lq rntrix, llat ion nwthnd: Slttor resisluiice per pliase. Rs:UcctuihicdlAssune sinusoidal hack ewf Tab

24、ulate v th d and q axis currentsHCancelHelpApply4D數(shù)據(jù)模型% Block Parameters: FEM-Parameterized PMSMXFEM-Paranetcr i zed PMSMThis block implements the electrical and nechanical characteristics of a permanent naRnet synchronous inotor for which magnetic flux linkage depends nonlinearl on stator currents

25、and rotor angle Right-click on the block and select Si inscape(TM) block choices to access variant implementationsSettingsElectricalIron LossesMechanicalA-phasecurremvector.iA： -200. 0. 0. 0, 200. 0B phasecurrentvector,iB: |-200.0. 0.0, 200.0：C phasecurrentvector,iC:卜200000, 2000Rotor angle vector i

26、hna:|xVec t.hA-phasc flux linkage partial derivative wrt iA. dFhiA(iA. iBJC, theta)/diA:A phase flux linkage partial derivative wrt dPhi.4(iA. iBJC, theta)/dill:A-pliase flux linkdKe partial derivative wrt iC. dPhiA(iAt iB, iC, theta)/diC:A-phase flux Iinkagc partial derivative wrt angle.dPhiA(iA. i

27、B. iC. theta)/ dtheta:zeros (3. 3, 3t 4)zeros(3, 3, 3, 4)zeros (3. 3. 3. 1)Torque matrix, T(iA, iB, iC, theta):0.3, 0; 0, 0F 0, 0; 0, - 0.3t 0.3, 0 X*wInterpolation method:Stator resistance per phase, Rs:LinearOKCancelApply3D數(shù)據(jù)模型，使用Park變換加 Block Parameters: FEM Parameterized PMSMXFElbParamcteri zed

28、PUSHThis block iraplemcnts the electrical and mechanical characteristics of a Dcrmanent macnct syncluofious motor for which magnetic flux linkage depends nonl ineai ly on stator currents and rotor anglc. Right-click on the block and select Simscapc(TM)block choices to access variant ipleaienlations

29、SettingsEl cctricalron LossesMeehani caIDirect axis current vector. iD:Quadrature-axis current vector. iQ:Rotor angle vector, theta:A-phase flux linkage part ial derivative irt iA, dPhiAfiD. iQ. theta)/diA:A-phase flux linkage part ial derivative wrt iB, dPhiACiD. iQ. rhcta)/diB:A-phase flux linkage part ial derivative vrt iC, dPhiA(iD. iQ. theta)/diC:A-phase flux linkage part ial derivative iTt angle, dPhiA(iD, iQ. theta)/ tithe ta:Torque iwtrix, T(iD, iQ. theta):lnterp(jlation ix