第7章基于MATLAB的交流電機仿真

1、第7章基于MATLAB的交流電機仿真7.1電力系統(tǒng)模塊集 Simulink中可以使用電力系統(tǒng)仿真模塊集（SimPowerSystems）。其功能非常強大，可以用于電路、電力電子系統(tǒng)、電機系統(tǒng)、電力傳輸?shù)冗^程的仿真，它提供了一種類似電路建模的方式進行模型繪制，在仿真前將自動將其變化成狀態(tài)方程描述的系統(tǒng)形式，然后才能在Simulink下進行仿真分析。 該模塊集下有許多子模塊集，雙擊每一個圖標都將打開下一級子模塊集。 電力系統(tǒng)模塊的使用與SIMULINK模塊的使用不同，電力系統(tǒng)模塊必須在回路中使用，因此每個模塊都用輸入端和輸出端，在回路中流動的都是電流，并且電流通過每個電氣元器件時產(chǎn)生電壓降。 SI

2、MULINK模塊的組成是信號流程，流入流出模塊的信號沒有特定的含義，其含義視仿真模型的對象而定。由電力系統(tǒng)模塊組成的電路和系統(tǒng)可以和SIMULINK模塊庫中的控制單元聯(lián)接、組合成控制系統(tǒng)，研究和觀察在不同控制方案下系統(tǒng)的穩(wěn)態(tài)和動態(tài)響應，為系統(tǒng)的設計提供依據(jù)。 7.1.1電源模塊庫(Electrical Sources)：交、直流電源，受控源，三相可編程電源，三相電源共七個電源模塊。 7.1.2電器元件模塊庫（Elements）：斷路器、子系統(tǒng)連接、分布參數(shù)傳輸線、接地、單線性變壓器、多繞組變壓器、互感線圈、等電位接點、RLC并聯(lián)電路、 RLC并聯(lián)負載、型參數(shù)傳輸線、飽和變壓器、 RLC串聯(lián)電路

3、、 RLC串聯(lián)負載、壓敏電阻、三相斷路器、三相動力負載、三相短路故障、三相濾波器、三相互感線圈、 三相RLC并聯(lián)電路、三相 RLC并聯(lián)負載、 三相型參數(shù)傳輸線、三相RLC串聯(lián)電路、 三相RLC串聯(lián)負載、三相三繞組變壓器、三相二繞組變壓器、十二端三相變壓器、之字形變壓器。 7.1.3其它電氣模塊庫（Extra Library） 7.1.4電機模塊庫（Machines）：異步電機（2）、直流電機（2）、同步電機勵磁系統(tǒng)、水輪機和調(diào)節(jié)器、電機測量單元、電力系統(tǒng)穩(wěn)定器、永磁式同步電機、同步電機簡化模型（2）、汽輪機和調(diào)節(jié)器、同步電機基本模型（2）、同步電機標準模型。 7.1.5測量儀器模塊庫（Meas

4、urements）：電流測量模塊、阻抗測量模塊、多路測量儀、三相電流電壓測量模塊、電壓測量模塊、 7.1.6狀態(tài)元件庫（Phasor Elements）：靜態(tài)無功補償器。 7.1.7電力電子元件模塊庫（Power Electronics）：詳細晶閘管模型、二極管、可關斷晶閘管、理想開關、絕緣柵雙極型晶體管、電力場效應管、三級橋、晶閘管、多功能橋式電路。7.2電路的仿真 電路中最常用的元件是電阻、電容和電感，雙擊電力系統(tǒng)模塊集中的Elements圖標，將出現(xiàn)元件的子模塊集。其中既包含各種電阻、電容和電感元件還包含變壓器元件、三相元件等。 從普通的電阻、電容和電感元件來看，有串聯(lián)的RLC分支和并聯(lián)

5、的RLC分支，以及它們的負載形式。 雙擊Series RLC Branch（串聯(lián)RLC分支）元件，將得出對話框，在對話框中適當?shù)剌斎腚娮?、電容和電感的參?shù)即可。填寫參數(shù)時應該注意其單位。 遺憾的是這里不包含單個的電阻、電容和電感元件，可以從串、并聯(lián)的分支來定義單獨的電路。具體情況見下表。 元件串聯(lián)RLC分支并聯(lián)RLC分支類型電阻數(shù)值電感數(shù)值電容數(shù)值電阻數(shù)值電感數(shù)值電容數(shù)值單個電阻R0 infRinf0 單個電感0 L infinfL 0 單個電容0 0 C infinfC 在一般電路中，除了一些元件外還需要一些連線器類模塊，雙擊電力系統(tǒng)模塊集中的Connectors打開連線器類模塊集。 例1考

6、慮如圖所示的感應電機的等效電路，輸入的交流電壓源為220V，50Hz，其它參數(shù)值為R1=0.428，L1=L2=1.926mH，R2=1.551，R3=1.803，L3=31.2mH。 步驟： 1、將所需的各電路元件復制到模型編輯窗口中。（對各元件點擊左鍵并按住拖入即可，對重復的元件可在編輯窗口中按右鍵拖動）。 2、對各元件進行賦值。雙擊元件圖標彈出參數(shù)對話框，修改參數(shù)并確定。 3、連接電路。Matlab7.0無連線器類模塊，只需按左鍵從一端頭拖至另一端頭即可。對分支按右鍵拖動。連接完成后還應該在電路的輸出端加一個電壓測量元件，并將輸出送給普通的Out1模塊。 4、仿真設置。建立好Simuli

7、nk模型后就可以啟動仿真過程，最簡單的方法是按下Simulink工具欄下的“啟動仿真按鈕”。啟動仿真過程后將以默認參數(shù)進行仿真。在Matlab命令窗口中鍵入plot(tout,yout);00.010.020.030.040.050.060.070.080.090.1-100-80-60-40-20020406080100 從仿真結(jié)果曲線看，該系統(tǒng)的輸出基本上是正弦波信號。 還可以對仿真電路進行穩(wěn)態(tài)分析。將電力系統(tǒng)模塊集中的powergui模塊復制到仿真框圖中，則雙擊該模塊就可以對電路圖中的信號進行穩(wěn)態(tài)分析，得出如下圖所示的分析結(jié)果。從得出的結(jié)果可以看出所測量的信號穩(wěn)態(tài)曲線為幅值94.124V

8、初始相位為-14.890的正弦信號，也即該信號的解析表達式在穩(wěn)態(tài)時趨近：089.14502sin124.94t 事實上，如果使用下面MATLAB語句將得出兩條幾乎重合的曲線。 y1=94.124*sin(2*pi*50*tout-14.89*pi/180);plot(tout,yout,tout,y1) 另外、用電力系統(tǒng)工具箱中提供的power2sys（）可以提取出從給定電源到輸出端子的狀態(tài)方程模型，根據(jù)此狀態(tài)方程模型就可以對整個電路進行頻域分析，如繪制Bode圖等。 power2sys（）函數(shù)的調(diào)用格式為： a，b，c，d=power2sys（模型名） 其中a，b，c，d為系統(tǒng)的狀態(tài)方程矩陣

9、，由tf（）還可以得出系統(tǒng)的傳遞函數(shù)模型。00.010.020.030.040.050.060.070.080.090.1-100-80-60-40-20020406080100 a,b,c,d=power2sys(ch7ex1)%獲得系統(tǒng)的狀態(tài)方程 a = -128.8763 -844.6462 -121.3833 -896.7868 b = 267.3783 251.8325 c = 0 1.8030 d = 0 G=ss(a,b,c,d);bode(G)%繪制系統(tǒng)的Bode圖-120-100-80-60-40-200Magnitude (dB)10-1100101102103104105

10、-90-4504590Phase (deg)Bode DiagramFrequency (rad/sec) 用戶可以自己設置需要的控制參數(shù)。 控制參數(shù)由Simulation/ Simulation parameters菜單選擇。 選擇該菜單項后將得到下圖所示的對話框，可以從中填寫相應的數(shù)據(jù)，控制仿真過程。 G1=tf(G)%獲得等效傳遞函數(shù)模型Transfer function: 454.1 s + 1.034e-010-s2 + 1026 s + 1.305e004 1、仿真區(qū)間設置。仿真起始、終了時間設置。 2、類型設置。步長選擇：定步長、變步長。 3、仿真算法選擇。定步長算法變步長算法

11、1、ode45 它是一種一步算法，對大多數(shù)仿真模型來說，首先使用ode45來解算模型是最佳的選擇，所以在SIMULINK的算法選擇中將ode45設為默認的算法。 2、ode23 它也是一種一步算法，在容許誤差和計算略帶剛性的問題方面，該算法比ode45算法要好。 3、ode113 它是可變階數(shù)的算法，在誤差要求很嚴時， ode113較ode45算法更適合。 4、ode15s 它是一種可變階數(shù)多步算法，當遇到帶剛性問題時或ode45算法不行時，可以試用這種算法。 5、ode23s 它是一種改進的二階算法，在容許誤差較大時ode23s比ode15s有效所以在解算一類帶剛性問題時用ode15s 不行

12、時，可以用ode23s算法。 6、 ode23t 一種采用自由內(nèi)插的梯形算法，如果模型有一定剛性，右要求解沒有數(shù)值衰減時，可以使用這種算法。 7、 ode23tb 在容許誤差比較大時， ode23tb和ode23t都比ode15s 要好。 8、discrete 處理離散系統(tǒng)的算法。 一般情況下，連續(xù)系統(tǒng)仿真應該選擇ode45變步長算法，對剛性問題可以選擇變步長的ode15s算法，離散系統(tǒng)一般默認地選擇定步長的discrete（no continous states）算法，而在仿真模型中含有連續(xù)環(huán)節(jié)時注意不能采用該仿真算法，而可以采用四階RungeKutta算法來求解問題。 定步長算法的步長應該

13、在Fixedstep size中填入?yún)?shù)，一般還可以選擇auto，依賴計算機自動選擇步長，而變步長下建議步長范圍使用auto選項。 在仿真過程中，用戶要根據(jù)各種類型模型的特點，各種數(shù)值積分方法的計算特點和適用范圍，才能正確的選擇恰當?shù)姆椒?，這一點是初學者難以掌握的，還沒有一種對所有模型都適用的算法。簡單的方法：當一種算法不能完成模型的計算時，選用另一種算法嘗試。若不行則要對模型和參數(shù)進行修改。 電力系統(tǒng)模塊集還提供了一些常用的電機仿真模塊，雙擊該模塊集中的Machines圖標，打開電機子模塊集，該模塊集中包含了各種電機模型，包括直流機（DC motors）、異步機（Asynchronous m

14、otors）、同步機（Synchronous motors）及各種其它形式等。7.3電機系統(tǒng)仿真限于篇幅，介紹星形連接下的電動機仿真。例2已知某電動機的參數(shù)為：Pn=5.5kw，Vn=380V，fn=50Hz，Rs=0.0217，xls=0.039，Rr=0.0329，xlr=0.0996，xm=3.6494，J=11.4kgm2，F(xiàn)=0，p=2。這里電感由電抗的形式給出，需用L=x/（2f）公式算出。電力系統(tǒng)工具箱中提供了兩個異步電動機模型，標幺值下的異步電動機模型（pu），國際單位制下的異步電動機模型（SI），在本例中考慮國際單位制的異步電動機模型，所以將該模型復制到空白的模型編輯窗口中。

15、電力系統(tǒng)工具箱中提供了兩個異步電動機模型，標幺值下的異步電動機模型（pu），國際單位制下的異步電動機模型（SI），在本例中考慮國際單位制的異步電動機模型，所以將該模型復制到空白的模型編輯窗口中。異步電動機模塊有4個輸入端子，4個輸出端子，前3個輸入端子（A，B，C）為電機的定子電壓輸入，一般可直接接三相電壓，可以有星形和三角形接法。第4輸入端一般接負載，為軸上機械轉(zhuǎn)矩，該端子可以直接接Simulink信號。模塊的前3個輸出端子（a，b，c）為轉(zhuǎn)子電壓輸出，一般短接在一起，或連接到其它的附加電路中。第4個輸出端為m端子，它返回一系列電機的內(nèi)部信號，共有21路信號組成，其構(gòu)成如下： 第1到第3路：

16、轉(zhuǎn)子電流ira, irb,irc；第4到第9路：qdn坐標系下的轉(zhuǎn)子信號，依次為q軸電流iqr，d軸電流idr，q軸磁通qr，d軸磁通dr，q軸電壓vqr，d軸電壓vdr；第10到第12路：定子電流isa, isb,isc；第13到第18路：qdn坐標系下的定子信號，依次為q軸電流iqr，d軸電流ids，q軸磁通qs，d軸磁通ds，q軸電壓vqs，d軸電壓vds；第13到第18路：電機轉(zhuǎn)速m，機械轉(zhuǎn)矩Tm，電機轉(zhuǎn)子角位移m。該路信號、應該接電機測試信號分路器（Machines Measurement Demux）模塊將各路信號分離出來，以便直接接顯示器或Simulink輸出端子進行顯示。雙擊異

17、步電動機模塊，將彈出該模塊的參數(shù)對話框，在對話框中需要輸入如下參數(shù)：繞組類型（Rotor type）列表框：分為繞線式（Wound）和鼠籠式（Squirrel-cage）兩種，后者將不顯示輸出端a，b，c，而直接在模塊內(nèi)部短接。參考坐標系（Reference frame）列表框：有靜止坐標系（Stationary），轉(zhuǎn)子的（Rotor）和同步的（Synchronous），一般常選擇靜止坐標系。額定參數(shù)額定功率Pn（單位：kw），線電壓Vn（單位：V），頻率fn（單位：Hz）；定子電阻Rs（Stator）（單位：Ohms）和漏感(Lls)（單位：H）；轉(zhuǎn)子電阻Rr（Rotor）（單位：Ohms）

18、和漏感(Llr)（單位：H）；互感(Mutual inductance)Lm（單位：H）；轉(zhuǎn)動慣量(inertia)J（單位：kgm2），摩擦系數(shù)F（單位：Nms），極對數(shù)。這些參數(shù)基本上都是電動機的銘牌參數(shù)。 測量信號分路器其默認的分路是針對簡化的同步機（Simplified synchronous）輸出的，雙擊該模塊彈出對話框，在電機類型（machine type）欄目中選擇異步（Asynchronous）機，則得出信號列表，從該列表中選出想輸出的信號。 異步電動機仿真框圖如下圖所示，電源采用了星形接法，故將這三個交流電壓源的電壓值設置成相電壓220V，A，B，C三相的相位分別填寫為0，1

19、20，240。 在該模型中有一個空閑的電流測試元件（Current Measurement），這是電力系統(tǒng)工具箱所要求的，因為它要求至少有一個測量元件。 另外還需要從輸出的信號中使用Selector元件提取所需的單路信號，這樣示波器中同時輸出4路信號，a相轉(zhuǎn)子電流iar，a相定子電流ias，轉(zhuǎn)速和轉(zhuǎn)矩T，設置仿真的終止時間為3秒，并選擇仿真算法為ode15s，相對允許誤差和絕對允許誤差均為10-7，對系統(tǒng)進行仿真，得到仿真結(jié)果。 將A、B兩相換序，即將A和B相的相位移分別改成120和0，則再進行仿真，即得到電機反轉(zhuǎn)時的仿真結(jié)果。仿真結(jié)果A、B互換相序仿真結(jié)果 例3利用simulink對方程組1

20、2212211xxxxxx 第1個方程可以認為是將 信號作為一個積分器的輸入端，這樣積分器的輸出則將成為 信號，類似地， 信號本身也可以認為是一個積分器的輸出，在積分器的輸入端信號應該為 ，在構(gòu)造該信號時還應該使用信號的乘積的處理，這樣可以按下圖所示的格式建立起描述該微分方程組的模型。除了上述的各個模塊，還需要添加結(jié)果輸出模塊，這里使用輸出端口模塊輸出結(jié)果。 tx2進行仿真。 tx2 tx112211xxx 可以看出，在系統(tǒng)框圖中，除了各個模塊及其連接之外，還給出了各個信號的文字描述。在Simulink模型中加文字描述的方法很簡單，在想加文字說明的位置雙擊鼠標，則將出現(xiàn)字符插入的標示，這時可將

21、任意的字符串寫入該位置即可。文字描述寫入模型后，則可用鼠標單擊并拖動到指定位置。 雙擊Gain模塊輸入u=1，雙擊兩個積分器模塊輸入初始條件x01=1，x02=-2。然后起動仿真命令。在MATLAB命令窗口中給出繪圖命令 plot(tout,yout) 將得出如下圖所示的時間響應曲線。 例4雙輸入雙輸出系統(tǒng)的狀態(tài)方程表示為012345678910-4-3-2-10123x1(t)x2(t)xyuxx20201000,2022426475. 025. 075. 125. 1125. 15 . 025. 025. 025. 125. 425. 25 . 025. 1525. 2 且輸入信號分別為s

22、in（t）cos（t），試構(gòu)造該系統(tǒng)的仿真框圖并仿真。仿真前在MATLAB命令窗口中輸入系統(tǒng)的狀態(tài)方程模型參數(shù) A=2.25 -5 1.25 -0.5;2.25 -4.25 -1.25 -0.25;0.25 -0.5 -1.25 -1;1.25 -1.75 -0.25 -0.75; B=4 6;2 4;2 2;0 2; C=0 0 0 1;0 2 0 2; D=zeros(2,2); 對整個系統(tǒng)進行仿真，則可以得出tout和yout兩個變量，其中yout的前2列為系統(tǒng)的輸出信號，后4列為系統(tǒng)的狀態(tài)變量。 用 plot(tout,yout(:,1:2)%系統(tǒng)的輸出曲線012345678910-2

23、0-15-10-505101520輸 出 信 號 y1(t)y2(t) plot(tout,yout(:,3:6)%系統(tǒng)的狀態(tài)曲線012345678910-10-8-6-4-20246810 x1x2x3x4例5利用Simulink自身的狀態(tài)方程模塊，對上例進行仿真構(gòu)建增廣矩陣 C=C;eye(4); D=D;zeros(4,2);選擇器參數(shù)設置Simulink模型 例6六極雙反饋發(fā)電機，額定功率7.5kw，額定頻率60Hz，額定電壓220V，rs=0.2943，rr =0.1442，定子漏電感Ls=0.00134H，轉(zhuǎn)子漏電感Lr=0.00056H，激磁電感Lm=0.03517H；轉(zhuǎn)動慣量J

24、=50kgm2。對任意速度旋轉(zhuǎn)的坐標系統(tǒng)dq0，在三相對稱條件下，雙饋發(fā)電機電壓方程及磁鏈方程有如下形式： dtdiRudtdiRudtdiRudtdiRuqrdrrqrrqrdrqrrdrrdrqsdsqssqsdsqsdssdsqrdrqsdssmsmmrmrqrdrqsdsLLLLLLLLLLLLLLLLiiii*00000000式中 srmLLLL2用矩陣方程表示為： UBIADCI 式中 A、B、C系數(shù)矩陣；D微分算符；I電流向量；U電壓向量。000000000000rrArrssRRRRB000000000000LLLLLLLLLLLLLLLLCsmsmmrmr00000000TqrdrqsdsiiiiI TqrdrqsdsuuuuU 消去磁鏈得： FEIDI式中 CBACCE1CUF 雙反饋發(fā)電機運動方程：JTiiiiMnnDmecrqrdsdrqsppr/23設雙反饋發(fā)電機端電壓為：32sin32sinsintUutUutUuescsesbsesas32sin32sinsintsUutsUutsUuercrerbrerar經(jīng)派克變換后，得：rqrdrsqsdsUuuUuu00 利用計算機對雙反饋發(fā)電機三相短路進行仿真分析時，