第3章控制系統(tǒng)建模-ppt課件

1、第3章 控制系統(tǒng)建模第第3章章 控制系統(tǒng)建?？刂葡到y(tǒng)建模 3.1 簡(jiǎn)單機(jī)械系統(tǒng)的建模簡(jiǎn)單機(jī)械系統(tǒng)的建模3.2 簡(jiǎn)單流體系統(tǒng)的建模簡(jiǎn)單流體系統(tǒng)的建模3.3 永磁體直流電機(jī)的建模永磁體直流電機(jī)的建模3.4 光源跟蹤伺服系統(tǒng)的建模與系統(tǒng)分析光源跟蹤伺服系統(tǒng)的建模與系統(tǒng)分析第3章 控制系統(tǒng)建模3.1 簡(jiǎn)單機(jī)械系統(tǒng)的建模簡(jiǎn)單機(jī)械系統(tǒng)的建模 3.1.1 彈簧振動(dòng)系統(tǒng)的建模 思索圖2.1所示的簡(jiǎn)單機(jī)械系統(tǒng)。 選擇垂直向下的方向?yàn)檎较?根據(jù)系統(tǒng)力平衡關(guān)系可以得到 mg-ks0=0 (3.1) 假設(shè)系統(tǒng)遭到正方向的外力,那么力平衡方程變成 202( )( )( )dFmy tmgksky tf tdt(3.2

2、)第3章 控制系統(tǒng)建模 圖3.1 彈簧振動(dòng)系統(tǒng)的表示圖第3章 控制系統(tǒng)建模 其中,y(t)是間隔平衡點(diǎn)的偏移間隔。以上是非阻尼條件下的系統(tǒng)方程。如今,假設(shè)系統(tǒng)浸入到一種粘性物質(zhì)中,那么系統(tǒng)將遭到與其瞬時(shí)速度方向相反的阻尼力的作用。當(dāng)系統(tǒng)以較慢速度運(yùn)動(dòng)時(shí),系統(tǒng)遭到的阻尼力與其運(yùn)動(dòng)的速度成正比,而方向相反。 假設(shè)這時(shí)的阻尼系數(shù)為常數(shù)c,整個(gè)系統(tǒng)的平衡方程為20222( )( )( )( )( )( )( )( )ddFmy tmgksky tcy tf tdtdtddmy tcy tky tf tdtdt(3.3) 第3章 控制系統(tǒng)建模 3.1.2 摩托車緩沖系統(tǒng)的建模 思索圖3.2所示的摩托車表

3、示圖。設(shè)計(jì)摩托車緩沖系統(tǒng)的目的是減小車輛在崎嶇道路上行駛時(shí)產(chǎn)生的震動(dòng)。道路外表的不平坦會(huì)引起摩托車沿垂直方向的挪動(dòng)和沿某個(gè)軸的轉(zhuǎn)動(dòng)。忽略輪胎的質(zhì)量,這樣整個(gè)系統(tǒng)由車架和駕駛員組成。 第3章 控制系統(tǒng)建模 圖3.2 摩托車系統(tǒng)表示圖 第3章 控制系統(tǒng)建模 摩托車緩沖系統(tǒng)的力平衡表示圖如圖3.3所示。 我們將整個(gè)系統(tǒng)的質(zhì)量中心作為坐標(biāo)的原點(diǎn),因此系統(tǒng)在不平道路上的振動(dòng)運(yùn)動(dòng)可以看作是質(zhì)心的沿垂直方向的平移運(yùn)動(dòng)以及沿質(zhì)心的旋轉(zhuǎn)運(yùn)動(dòng)。摩托車架以及駕駛員可以整個(gè)視作質(zhì)量為M,轉(zhuǎn)動(dòng)慣量為J的剛體。輸入車輪的位置信息Ya、Yb闡明路況信息。假設(shè)每個(gè)車軸的緩沖系統(tǒng)由具有阻尼特性的彈簧構(gòu)成。因此,每個(gè)車輪遭到的外

4、力為彈簧彈力與阻尼力之和,即()( )()( )()( )()( )aaaaaaabbbbbbbdFaky tc skysdtdFaky tc sky sdt (3.4)第3章 控制系統(tǒng)建模圖3.3 摩托車緩沖系統(tǒng)的力平衡表示圖第3章 控制系統(tǒng)建模 ya和yb分別表示每個(gè)彈簧間隔參考位置的瞬時(shí)間隔。 用Y(t)和(t)分別表示系統(tǒng)質(zhì)心的平移位移和沿質(zhì)心的旋轉(zhuǎn)角度。對(duì)于單個(gè)彈簧有()()aaabbbyYLYyYLY(3.5) 上式中假定在很小的角度位置條件下滿足sin=,并且取逆時(shí)針的旋轉(zhuǎn)方向?yàn)檎较?，如圖3.4所示。第3章 控制系統(tǒng)建模圖3.4 摩托車緩沖系統(tǒng)垂直位置與旋轉(zhuǎn)角度的幾何分析第3章

5、 控制系統(tǒng)建模將式3.5代入式3.4中,得到 Fa=(cas+ka)Ya-(Y-La) Fb=(cbs+kb)Yb-(Y+Lb) (3.6) 或者定義Za=cas+ka,Zb=cbs+kb,得到 Fa=ZaYa-(Y-La) Fb=ZbYb-(Y+Lb) (3.7)第3章 控制系統(tǒng)建模 或者 Ms2Y=ZaYa-(Y-La)+ZbYb-(Y+Lb) 整理后得到 (Ms2+Za+Zb)Y-(ZaLa-ZbLb)=ZaYa+ZbYb (3.9)最后根據(jù)牛頓第二運(yùn)動(dòng)定律,有22abdMYFFdt第3章 控制系統(tǒng)建模 上式給出了摩托車緩沖系統(tǒng)的力平衡方程,同時(shí)假定車架和駕駛員在初始位置沒有垂直方向上的

6、速度Y0=0,dY/dt|0=0。 假設(shè)對(duì)上述系統(tǒng)建立關(guān)于質(zhì)心位置的力矩平衡方程,可以得到另一個(gè)系統(tǒng)方程，即 22coscos(0,cos1)bbaabbaadJF LF LF LF Ldt(3.10) 或者 Js2=ZbLbYb-(Y+Lb)-ZaLa Ya-(Y-La)第3章 控制系統(tǒng)建模 (Js2+ZaL2a+ZbL2b)-(ZaLa-ZbLb)Y =-ZaLaYa+ZbLbYb (3.11) 再次假定初始條件為零(0=0，d/dt|0=0)，最后將力和力矩平衡方程寫成矩陣方式 整理得2222()()()()aababaabbaabbbaabbaabbYZZYMsZZZ LZ LZ LZ

7、 LYZ LZ LJsZ LZ L (3.12) 第3章 控制系統(tǒng)建模 寫成簡(jiǎn)化方式 1112111221222122abYAABBYAABBY (3.13) 用Aij和Bij可以表示Y和為111212112122212211ababYB YB YAAB YB YAA(3.14) 最終系統(tǒng)模型可以用如圖3.5所示的框圖表示。第3章 控制系統(tǒng)建模圖3.5 摩托車緩沖系統(tǒng)的方框圖 第3章 控制系統(tǒng)建模 以上系統(tǒng)中假定Ya和Yb是系統(tǒng)兩個(gè)相互獨(dú)立的輸入變量,但實(shí)踐上,后輪與前輪的位置信號(hào)相差t=L/V時(shí)間。這樣,實(shí)踐系統(tǒng)滿足Yb(t)=Ya(t-t)。 假設(shè)定義系統(tǒng)形狀分別為Y、dy/dt和d/dt

8、,還可以計(jì)算出系統(tǒng)的形狀方程描畫。另外一種得到整個(gè)系統(tǒng)傳送函數(shù)的方法是經(jīng)過模型方框圖進(jìn)展計(jì)算。然后，在此根底上可以對(duì)該系統(tǒng)進(jìn)展時(shí)域和頻域的仿真,詳細(xì)計(jì)算過程留給讀者練習(xí)。第3章 控制系統(tǒng)建模3.2 簡(jiǎn)單流體系統(tǒng)的建模簡(jiǎn)單流體系統(tǒng)的建模 3.2.1 單個(gè)蓄水槽的動(dòng)態(tài)模型 思索圖3.6所示的單個(gè)蓄水槽模型,其槽底的液體流出速度是由槽內(nèi)的液壓決議的。 各部分的含義為： A蓄水槽的外表區(qū)域； V水槽的容積；Ae水槽出口處的連通部分； P1槽底的液壓。第3章 控制系統(tǒng)建模 圖3.6 單個(gè)蓄水槽模型 第3章 控制系統(tǒng)建模 液體的輸出壓強(qiáng)為Pa,輸出液體的速率作為系統(tǒng)的輸入。系統(tǒng)的形狀變量包括槽內(nèi)液體的高度

9、,其系統(tǒng)輸出為液體流出的速率We。根據(jù)系統(tǒng)的物質(zhì)平衡,可得到iedmdt(3.15) 假設(shè)蓄水槽的周圍壁是垂直的,槽內(nèi)液體的質(zhì)量是液體的密度乘以液體的體積，有ddddmAhAhAhdtdtdtdtmAh(3.16) (3.17) 第3章 控制系統(tǒng)建模輸出液體的質(zhì)量可以寫成輸出速率的函數(shù)eeeA(3.18) 根據(jù)出口處的能量平衡w=w1=w2,可以得到2212121212()()()()02uug zzPP(3.19) 假定整個(gè)系統(tǒng)不存在能量或物質(zhì)的滯留,并且忽略內(nèi)部能量的變化u1=u2,z1=z2,那么根據(jù)能量守恒原理得到11221222()()2eaaPPPghPgh(3.20)第3章 控制

10、系統(tǒng)建模 顯然該系統(tǒng)的形狀方程是一階非線性的,槽內(nèi)液體質(zhì)量的瞬時(shí)變化等于輸入的液體速率減去輸出的液體速率,槽內(nèi)的液體質(zhì)量和輸出液體的速率都視作t時(shí)辰液體高度的函數(shù)。 下面來定義液體的阻力作用,它的作用與電路中的電阻作用類似：綜合上面的方程得出 2eidAhAghdt (3.21)( )2eeeehhR hAAgh(3.22) 第3章 控制系統(tǒng)建模 與此類似,還可以定義該系統(tǒng)的電容效應(yīng),它反映了蓄水槽存儲(chǔ)液體容量的變化。AhCAh(3.23) 將式3.22和式3.23代入式3.21,得到1( )( )( )dCx tx tu tdtR (3.24) 或者 11( )( )( )dx tx tu

11、tdtRCC (3.25) 第3章 控制系統(tǒng)建模 該非線性系統(tǒng)也可以線性化并寫成規(guī)范的形狀方程方式。定義系統(tǒng)的參考位置 h(t)=h0(t)+h(t) (3.26) wi(t)=wi0(t)+wi(t) (3.27)系統(tǒng)方程3.21可以寫成00000( , )(, )(, )(, )iiihihdhf h tdtAddhhf h tf h thdtdthAAdhf h thdthA(3.28) 第3章 控制系統(tǒng)建模參考位置處的系統(tǒng)Jacobian矩陣為00120( )21( )(2)22eeehhAf hghAAAgf hghhAAgh (3.29) (3.30) 最后將式3.30代入式3.2

12、8，得到 01211( )( )( )eidAghhdtAAghdx tx tu tdtRCC (3.31) (3.32) 第3章 控制系統(tǒng)建模圖3.7 單個(gè)蓄水槽的時(shí)域特性 第3章 控制系統(tǒng)建模圖3.8 系統(tǒng)線性化模型的Bode圖 第3章 控制系統(tǒng)建模圖3.9 線性化模型的其它頻域特性 (Nyquist曲線和Nichols圖)第3章 控制系統(tǒng)建模其相關(guān)的程序包括： TANKCS1.M：仿真程序TANK1L.M：系統(tǒng)線性化模型描畫TANK1NL.M：非線性系統(tǒng)模型描畫程序代碼如下：% TANK1L.M 單個(gè)蓄水槽的線性化模型 functionxdot=tank1l(t,x) globalc1l

13、c2l xdot=c1l*x+c2l;第3章 控制系統(tǒng)建模% TANKCS1.M 單個(gè)蓄水槽動(dòng)態(tài)特性的仿真程序 clearall, closeall, nfig=0;% 定義ODE求解器中運(yùn)用的全局變量 globalc1nlc1lc2nlc2l% 定義相關(guān)參數(shù) too=0; tff=25; %仿真計(jì)算的啟動(dòng)和終了時(shí)間(min) rho=62.4; %水的密度(lbm/ft3) g=1.1592e5; %重力加速度常數(shù)(ft/min2) wio=1900; %參考位置的入口處液體流速(lbm/min) ho=10; %參考位置處的液位高度(ft) a=5; %槽底的外表積(ft2) ae=.02

14、; %出水管的外表積(ft2)第3章 控制系統(tǒng)建模% deltainput dw1=0.1*wio; %10%的變化 dw2=0.5*wio; %50%的變化% 非線性系數(shù) c1nl=-ae*sqrt(2*g)/a; % 線性模型系數(shù) r=sqrt(2*g*ho)/(rho*g*ae); c=rho*a; c1l=-1/(r*c); % 非線性仿真 - 脈沖&階躍呼應(yīng)1 y1o=ho+dw1/c; c2nl=wio/c; 第3章 控制系統(tǒng)建模t1i,y1i=ode23(tank1nl,too,tff,y1o); y1o=ho; c2nl=(wio+dw1)/c; t1s,y1s=ode

15、23(tank1nl,too,tff,y1o);% 線性仿真 - 脈沖&階躍呼應(yīng)1 y1o=dw1/c; c2l=0; t2i,y2i=ode23(tank1l,too,tff,y1o); y2i=ho+y2i; y1o=0; c2l=dw1/c; t2s,y2s=ode23(tank1l,too,tff,y1o); y2s=ho+y2s;第3章 控制系統(tǒng)建模% 繪制結(jié)果1 nfig=nfig+1; figure(nfig) subplot(2,2,1),plot(t1i,y1i,r,t2i,y2i,g:),grid title(ImpulseResponse(10%); xlabel

16、(Time(min),ylabel(Height(ft); legend(NL,L); subplot(2,2,2),plot(t1s,y1s,r,t2s,y2s,g:),grid title(StepResponse(10%); xlabel(Time(min),ylabel(Height(ft); legend(NL,L); 第3章 控制系統(tǒng)建模% 非線性仿真 - 脈沖&階躍呼應(yīng)2 y1o=ho+dw2/c; c2nl=wio/c; t1i,y1i=ode23(tank1nl,too,tff,y1o); y1o=ho; c2nl=(wio+dw2)/c; t1s,y1s=ode23

17、(tank1nl,too,tff,y1o);% 線性仿真 - 脈沖&階躍呼應(yīng)2 y1o=dw2/c; c2l=0; t2i,y2i=ode23(tank1l,too,tff,y1o); y2i=ho+y2i; y1o=0; c2l=dw2/c; t2s,y2s=ode23(tank1l,too,tff,y1o); y2s=ho+y2s;第3章 控制系統(tǒng)建模% 繪制結(jié)果2 subplot(2,2,3),plot(t1i,y1i,r,t2i,y2i,g:),grid title(ImpulseResponse(50%); xlabel(Time(min),ylabel(Height(ft)

18、; legend(NL,L); subplot(2,2,4),plot(t1s,y1s,r,t2s,y2s,g:),grid title(StepResponse(50%); xlabel(Time(min),ylabel(Height(ft); legend(NL,L);第3章 控制系統(tǒng)建模 % 線性化模型的頻域仿真 b=r; a=r*c1; sys1=tf(b,a); w=logspace(-2,2,200); %仿真的頻域范圍(rad/min) nfig=nfig+1; figure(nfig) bode(sys1,w),xlabel(Frequency(rads/min) nfig=n

19、fig+1; figure(nfig) subplot(2,1,1),nichols(sys1,w) subplot(2,1,2),nyquist(sys1,w) Re,Im=nyquist(sys1,w); Resys1(:,1)=Re(1,1,:); Imsys1(:,1)=Im(1,1,:);第3章 控制系統(tǒng)建模 3.2.2 雙蓄水槽系統(tǒng)的模型建立 在實(shí)踐系統(tǒng)中經(jīng)常存在相互連通的多個(gè)蓄水槽系統(tǒng),這使得系統(tǒng)呈現(xiàn)出更高階動(dòng)態(tài)系統(tǒng)的特性。下面以雙連通的蓄水槽系統(tǒng)如圖3.10所示為對(duì)象,研討二階系統(tǒng)的動(dòng)態(tài)行為。第3章 控制系統(tǒng)建模圖3.10 雙蓄水槽系統(tǒng)構(gòu)造表示圖 第3章 控制系統(tǒng)建模 按照前一

20、節(jié)對(duì)單個(gè)蓄水槽的分析,可以確定1號(hào)和2號(hào)蓄水槽的液體輸出流量分別為 1221222eehhhRR(3.33) 從而得到兩個(gè)蓄水槽的平衡方程為 121111222212idhhChdtRdhhhChdtRR(3.34) 第3章 控制系統(tǒng)建模 其中Ci(i=1,2)的定義與單個(gè)蓄水槽系一致樣。 將上面的結(jié)果寫成規(guī)范的形狀方程方式,并且將液體的輸出流速作為系統(tǒng)輸出,那么111111122121221222221111111()010iehRCRChdChhdtRCRRChhyhRR(3.35) 第3章 控制系統(tǒng)建模 上述方程中的h1、h2和wi實(shí)踐上是相對(duì)于參考位置的偏移量即h1=h1等。 為了創(chuàng)建

21、系統(tǒng)的Simulink方框圖,只需針對(duì)形狀方程中每個(gè)SISO子系統(tǒng)繪制它的框圖如圖3.11所示,然后將它們銜接起來即可。最后得到的仿真框圖如圖3.12所示。第3章 控制系統(tǒng)建模圖3.11 各個(gè)子系統(tǒng)的框圖 第3章 控制系統(tǒng)建模圖3.12 整個(gè)系統(tǒng)的框圖構(gòu)造 第3章 控制系統(tǒng)建模 整個(gè)系統(tǒng)的傳送函數(shù)可以經(jīng)過矩陣運(yùn)算或相關(guān)的模塊框圖計(jì)算得到12112222221111211111111( )( )()(111( )()()Y sHsWYHHYRRC sR C s RH sUWUHHC sC sR利用H2=R2Y簡(jiǎn)化上式,得到 12111( )1URYRH sRC s第3章 控制系統(tǒng)建模 將上式代入

22、式3.36,得22122222211( )11RYURRY sYYR C s R C sR C sR 整理后得到整個(gè)系統(tǒng)的傳送函數(shù)為22212221111( ) ( )()1Y sU sR C C sR CR CRC s第3章 控制系統(tǒng)建模 圖3.13和3.14為該系統(tǒng)的脈沖和階躍呼應(yīng)以及頻域呼應(yīng)的部分仿真結(jié)果,相關(guān)程序文件為TANKCS2.M。程序?qū)﹄p蓄水槽系統(tǒng)的兩種不同情況進(jìn)展了仿真。一種是對(duì)系統(tǒng)參考位置的動(dòng)態(tài)行為進(jìn)展了仿真,第二種情況那么將1號(hào)槽的容積添加了4倍,這將導(dǎo)致系統(tǒng)時(shí)間常數(shù)的添加,讀者可以從圖3.13中清楚地了解到這一點(diǎn)。第3章 控制系統(tǒng)建模圖3.13 雙蓄水槽系統(tǒng)的時(shí)域呼應(yīng)第

23、3章 控制系統(tǒng)建模圖3.14 雙蓄水槽系統(tǒng)的Bode圖 第3章 控制系統(tǒng)建模以下是相關(guān)的程序代碼： % TANKCS2.M 雙蓄水槽系統(tǒng)的動(dòng)態(tài)仿真 clearall, closeall, nfig=0; too=0;tff=60;t=too:1:tff;t=t; %仿真的時(shí)間范圍(min) w=logspace(-3,2,100); %頻域仿真的頻率范圍(rad/min)第3章 控制系統(tǒng)建模% Example1 wio=1500; %參考位置處液體入口的液體流量(lbm/min) dw=150; %輸入流量10%的變化(lbm/min) r1=.01; %1號(hào)槽出口處的阻力常數(shù)(ft.min/

24、lbm) r2=.01; %2號(hào)槽出口處的阻力常數(shù)(ft.min/lbm) c1=250; %1號(hào)槽出口處的電容效應(yīng)常數(shù)(lbm/ft) c2=250; %2號(hào)槽出口處的電容效應(yīng)常數(shù)(lbm/ft) h1=30; %參考位置時(shí)1號(hào)槽的液位高度(ft) h2=15; %參考位置時(shí)2號(hào)槽的液位高度(ft)第3章 控制系統(tǒng)建模% 創(chuàng)建形狀方程描畫 a11=-1/(r1*c1); a12=-a11; a21=1/(r1*c2); a22=-(1/r1+1/r2)/c2; b1=1/c1; b2=0; cc1=0; cc2=1/r2; d=0; am1=a11a12;a21a22; bv1=b1b2;

25、cv1=cc1cc2; sys1=ss(am1,dw*bv1,cv1,d); % Example2 c1=1000; %1號(hào)槽出口處的電容效應(yīng)常數(shù)(lbm/ft) c2=250; %2號(hào)槽出口處的電容效應(yīng)常數(shù)(lbm/ft)第3章 控制系統(tǒng)建模% 創(chuàng)建形狀方程描畫 a11=-1/(r1*c1); a12=-a11; a21=1/(r1*c2); a22=-(1/r1+1/r2)/c2; b1=1/c1; b2=0; cc1=0; cc2=1/r2; d=0; am2=a11a12;a21a22; bv2=b1b2; cv2=cc1cc2; sys2=ss(am2,dw*bv2,cv2,d);

26、第3章 控制系統(tǒng)建模% 對(duì)上述兩個(gè)例子的脈沖和階躍呼應(yīng)進(jìn)展仿真 yi1,t,xi1=impulse(sys1,t); ys1,t,xs1=step(sys1,t); yi2,t,xi2=impulse(sys2,t); ys2,t,xs2=step(sys2,t);% 時(shí)域仿真曲線 nfig=1; figure(nfig) subplot(2,2,1),plot(t,yi1,r),grid title(ImpulseResponse(V1=Vr); xlabel(Time(min),ylabel(DelFlow(lb/m); subplot(2,2,2),plot(t,ys1,r),grid第

27、3章 控制系統(tǒng)建模title(StepResponse(V1=Vr); xlabel(Time(min),ylabel(DelFlow(lb/m); subplot(2,2,3),plot(t,yi2,g),grid title(ImpulseResponse(V1=4*Vr); xlabel(Time(min),ylabel(DelFlow(lb/m); subplot(2,2,4),plot(t,ys2,g),grid title(StepResponse(V1=4*Vr); xlabel(Time(min),ylabel(DelFlow(lb/m);第3章 控制系統(tǒng)建模% 上述兩個(gè)例子的

28、Bode圖 nfig=nfig+1; figure(nfig), bode(sys1,w); xlabel(Frequency(rad/min) nfig=nfig+1; figure(nfig), bode(sys2,w); xlabel(Frequency(rad/min)第3章 控制系統(tǒng)建模3.3 永磁體直流電機(jī)的建模永磁體直流電機(jī)的建模 3.3.1 物理描畫 永磁體直流DC電機(jī)是動(dòng)態(tài)系統(tǒng)中經(jīng)常用到的設(shè)備。這一節(jié)將引見如何對(duì)規(guī)范的DC電機(jī)建模。在了解電機(jī)的任務(wù)原理之后,運(yùn)用形狀方程和Laplace變換對(duì)電機(jī)進(jìn)展進(jìn)一步的分析,并將它運(yùn)用到前一節(jié)的光源跟蹤 伺 服 系 統(tǒng) 中 。 取 定 一

29、 組 電 機(jī) 參 數(shù) 后 , 在MATLAB/Simulink中分析該系統(tǒng)的根本時(shí)域特性。第3章 控制系統(tǒng)建模 永磁體直流DC電機(jī)是一個(gè)經(jīng)過磁性耦合將電能轉(zhuǎn)化成機(jī)械能的安裝。其最根本的組成包括兩部分：轉(zhuǎn)子或電樞和定子。電樞將在電機(jī)的定子框架中旋轉(zhuǎn)。圖3.15是該電機(jī)的根本構(gòu)造表示圖。 電機(jī)定子是由永磁體構(gòu)成的,用來產(chǎn)生磁場(chǎng)。轉(zhuǎn)子是由繞在鐵芯上的電磁線圈構(gòu)成的。轉(zhuǎn)子產(chǎn)生的電磁場(chǎng)與定子產(chǎn)生的永磁場(chǎng)相互作用,從而使轉(zhuǎn)子發(fā)生旋轉(zhuǎn)。電機(jī)的換向器是由安裝在轉(zhuǎn)子末端軸上的兩片半圓型金屬銅片構(gòu)成的,每個(gè)轉(zhuǎn)子繞組的末端都與其中的一片銅片相連。靜止的電刷緊貼換向器,為轉(zhuǎn)子繞組提供直流電流。第3章 控制系統(tǒng)建模圖3

30、.15 永磁體直流DC電機(jī)的根本構(gòu)造表示圖 第3章 控制系統(tǒng)建模 3.3.2 數(shù)學(xué)模型 圖3.16是DC電機(jī)的等效電路圖。其中,La和Ra為轉(zhuǎn)子繞組的等效電感和等效電阻,Vc是轉(zhuǎn)子繞組旋轉(zhuǎn)時(shí)切割永磁場(chǎng)的磁力線產(chǎn)生的內(nèi)部電壓常稱為電動(dòng)勢(shì)。 根據(jù)電路回路的電壓平衡關(guān)系,得到 0aRaLacRaaaLaaacvaVVVVVi RdVLidtVk同時(shí) 第3章 控制系統(tǒng)建模圖3.16 DC電機(jī)的等效電路圖 第3章 控制系統(tǒng)建模 其中，kv是由永磁體的磁通密度、轉(zhuǎn)子繞組的數(shù)目以及鐵芯的物理性質(zhì)決議的速度常數(shù)。a是轉(zhuǎn)子旋轉(zhuǎn)的角速度。 將上面的式子代入方程3.38中,得到0aaaaavadVi RLikdt再

31、根據(jù)電機(jī)的力矩平衡關(guān)系可以得到 0eLTTTT(3.39) (3.40) 第3章 控制系統(tǒng)建模 其中,Te是電機(jī)的電磁力矩,T是驅(qū)動(dòng)轉(zhuǎn)子加速度運(yùn)動(dòng)的力矩,T是轉(zhuǎn)子速度產(chǎn)生的力矩,而TL那么是電機(jī)的負(fù)載力矩。 我們知道,電機(jī)的電磁力矩與電流的大小成正比：et aTk i 其中,kt是由永磁體的磁通密度、轉(zhuǎn)子繞組的數(shù)目以及鐵芯的物理性質(zhì)決議的力矩常數(shù)。T可以寫成adTJdt第3章 控制系統(tǒng)建模 J是轉(zhuǎn)子和電機(jī)負(fù)載的轉(zhuǎn)動(dòng)慣量。與轉(zhuǎn)動(dòng)速度相關(guān)聯(lián)的力矩為 ATB B是整個(gè)機(jī)械旋轉(zhuǎn)系統(tǒng)的阻尼常數(shù)。將上面得到的關(guān)系代入式3.40,得0t aaaLdk iJBTdt(3.41) 第3章 控制系統(tǒng)建模 聯(lián)立式3

32、.39和3.41,得到DC電機(jī)的完好描畫：avaaaaaaaLaaadRkViidtLLLdkBTidtJJJ (3.42) (3.43) 第3章 控制系統(tǒng)建模 寫成形狀空間方式 1201010000100avaaaaaaataaaLRkBiiLLLddtkBJJJiVyyT(3.44) (3.45) 第3章 控制系統(tǒng)建模 下面計(jì)算該系統(tǒng)的傳送函數(shù)。對(duì)式3.42和3.43進(jìn)展Laplace變換,得到 1( )(0)( )( )( )1( )(0)( )( )( )avaaaaaaaataaaaLRksIsIIssV sLLLkBssIssT sJJJ 假設(shè)思索穩(wěn)定形狀周圍的干擾,并且假設(shè)電機(jī)的

33、初始條件為零，那么一切的變量都是指間隔參考形狀的偏移量,上述方程變成 ( )( )( )( )( )( )vaaaaataLaksV sIsL sRk IsT ssJsB(3.48) (3.49) 第3章 控制系統(tǒng)建模 以上方程可以很容易寫成模塊框圖方式。由此得到的永磁體DC電機(jī)的框圖如圖3.17所示。 假設(shè)電機(jī)的負(fù)載為常數(shù),并且只需求輸出電機(jī)轉(zhuǎn)動(dòng)的角速度,那么圖3.17可以進(jìn)一步簡(jiǎn)化成圖3.18所示的框圖,并由此得到整個(gè)電機(jī)的傳送函數(shù)為2( )( )( )()taaaaatvaaaksL JG sR JL BR Bk kV sssL JL J(3.50)第3章 控制系統(tǒng)建模圖3.17 永磁體

34、DC電機(jī)的模塊框圖 第3章 控制系統(tǒng)建模圖3.18 永磁體DC電機(jī)的簡(jiǎn)化框圖 第3章 控制系統(tǒng)建模 3.3.3 仿真分析 按照?qǐng)D3.17和圖3.18可以在Simulink中建立永磁體DC電機(jī)的仿真框圖如圖3.19所示,直接察看系統(tǒng)的時(shí)域呼應(yīng)。 下面分別采用兩種方法來進(jìn)展DC電機(jī)的時(shí)域仿真。一種是直接運(yùn)用式3.50所得到的整個(gè)電機(jī)的傳送函數(shù)。另一種方法那么是采用所建立的Simulink仿真框圖。第3章 控制系統(tǒng)建模圖3.19 永磁體直流電機(jī)的Simulink仿真框圖 (MOTORSL.MDL)第3章 控制系統(tǒng)建模相關(guān)的程序代碼如下：% MOTORTST.M 永磁體直流電機(jī)的時(shí)域動(dòng)態(tài)特性仿真程序

35、clearall, closeall% 各種參數(shù)定義 Ra=1.75; La=2.83e-3; kv=0.093; kt=0.0924; Jeq=30.e-6; Beq=5.0e-3;% CASE1:利用整個(gè)電機(jī)的傳送函數(shù) num=kt/(La*Jeq); d e n = 1 ( R a * J e q + L a * B e q ) / ( L a * J e q ) (Ra*Beq+kt*kv)/(La*Jeq); sys1=tf(num,den);第3章 控制系統(tǒng)建模% CASE2:利用Simulink仿真框圖 A,B,C,D=linmod(motorsl); sys2=ss(A,B,C

36、,D);% 確定兩種情況下電機(jī)的階躍呼應(yīng) t=linspace(0,0.025,51); y1=step(sys1,t); y2=step(sys2,t);% 繪制相關(guān)結(jié)果 plot(t,y1,ro,t,y2,b-),grid xlabel(time(sec),ylabel(angularvelocity(rad/sec) title(MotorResponsetoStepChangeinAppliedVoltage) text(.0152,.75,Case1TransferFunction-points) text(.0152,.45,Case2SimulinkModel-solidline

37、)第3章 控制系統(tǒng)建模圖3.20 永磁體DC電機(jī)的階躍呼應(yīng)曲線第3章 控制系統(tǒng)建模3.4 光源跟蹤伺服系統(tǒng)的建模與系統(tǒng)分析光源跟蹤伺服系統(tǒng)的建模與系統(tǒng)分析 3.4.1 物理模型 太陽跟蹤系統(tǒng)是太陽能系統(tǒng)不可短少的組成部分,它可以有效地提高太陽光線的利用效率。 下面將研討一種光源跟蹤侍服系統(tǒng)模擬太陽電池帆板的挪動(dòng)。光源跟蹤侍服系統(tǒng)是由永磁體直流電機(jī)、光線檢測(cè)電路和電機(jī)驅(qū)動(dòng)放大器等幾部分組成的。該系統(tǒng)的模型框圖如圖3.21所示。 第3章 控制系統(tǒng)建模圖3.21中各部分的含義為：K比例增益系數(shù)；E誤差信號(hào)；Va電機(jī)驅(qū)動(dòng)電壓；Gm電機(jī)傳送函數(shù)；W電機(jī)的轉(zhuǎn)動(dòng)角速度；A電池帆板的角度位置；L光源的角度方向

38、；C0系統(tǒng)開環(huán)傳送函數(shù)；GC系統(tǒng)閉環(huán)傳送函數(shù)。第3章 控制系統(tǒng)建模圖3.21 光源跟蹤伺服系統(tǒng)的仿真框圖 第3章 控制系統(tǒng)建模 3.4.2 數(shù)學(xué)模型 下面我們來建立光源跟蹤伺服系統(tǒng)的數(shù)學(xué)模型。整個(gè)系統(tǒng)的傳送函數(shù)是根據(jù)檢測(cè)電路、功率放大器和驅(qū)動(dòng)電機(jī)的傳送函數(shù)綜合得來的。 1) 圖像檢測(cè)電路和功率放大器 圖像檢測(cè)電路和功率放大器可以整個(gè)視作增益環(huán)節(jié)K,這里K的單位為V/rad。 2) 驅(qū)動(dòng)電機(jī) 在該伺服系統(tǒng)中,驅(qū)動(dòng)電機(jī)采用永磁體直流電機(jī),它的傳送函數(shù)和系統(tǒng)特性在前一節(jié)曾經(jīng)詳細(xì)討論過,這里不再贅述。 第3章 控制系統(tǒng)建模 根據(jù)以上討論可以很容易地創(chuàng)建該伺服系統(tǒng)的頻域方框圖如圖3.22所示。圖中顯示,

39、為了實(shí)現(xiàn)太陽電池自動(dòng)跟蹤光源變化,采用了比例反響的控制方法。圖3.22 光源跟蹤伺服系統(tǒng)的控制框圖 第3章 控制系統(tǒng)建模 根據(jù)圖3.22可以得到整個(gè)系統(tǒng)的傳送函數(shù)32( )()()taaaatvtaaaKkL JG sR JL BR Bk kKksssL JL JL J(3.51)第3章 控制系統(tǒng)建模 3.4.3 仿真分析 從式3.52表示的系統(tǒng)傳送函數(shù)中可以看出,存在三個(gè)要素會(huì)影響系統(tǒng)的動(dòng)態(tài)行為和控制系統(tǒng)的穩(wěn)定性。這三個(gè)要素是： (1) 圖像檢測(cè)電路和放大器的等效增益。 (2) 電機(jī)負(fù)載的轉(zhuǎn)動(dòng)慣量J。 (3) 機(jī)械系統(tǒng)的阻尼系數(shù)B。 下面的程序LTSERVO.M可以讓用戶自在改動(dòng)其中的參數(shù),

40、并且對(duì)相應(yīng)的階躍呼應(yīng)和頻域特性進(jìn)展仿真。第3章 控制系統(tǒng)建模% LTSERVO.M 比例反響控制條件下光源跟蹤伺服系統(tǒng)的仿真和參數(shù)敏感性分析 clearall, closeall% 驅(qū)動(dòng)電機(jī)參數(shù) Ra=1.75; La=2.83e-3; kv=0.093; kt=0.0924;% 機(jī)械系統(tǒng)參數(shù) Jeq=30.0e-6; Beq=5.0e-3;第3章 控制系統(tǒng)建模% 利用根軌跡方法尋覓缺省的參考增益K num=kt/(La*Jeq); d2=Ra/La+Beq/Jeq; d3=(Ra*Beq+kt*kv)/(La*Jeq); den=1d2d30; kk1=linspace(0,12,7); k

41、k2=linspace(12,18,19); kk3=linspace(18,30,7); kk4=linspace(30,200,15); kk=kk1kk2kk3kk4; figure(1) sys=tf(num,den); rlocus(sys,kk),sgrid axis(-50050-500500)第3章 控制系統(tǒng)建模 cont=input(Selectrootlocationforgaindetermination?. (y/n)n:,s); ifisempty(cont); cont=n; end whilecont=y K,p=rlocfind(sys); disp(GainK

42、=); K syscl=feedback(K*sys,1); damp(syscl) cont=input(Selectanotherroot?(y/n):, s); ifisempty(cont); cont=n; end end第3章 控制系統(tǒng)建模%-% 敏感度分析(根據(jù)需求反復(fù)多次)%- repeat=y; whilerepeat=y; Jeq=30.0e-6; Beq=5.0e-3; K=25.0; n=menu(Selectvariabletoanalyze,GAINParameter, . DAMPING Term, . INERTIAofSystem);第3章 控制系統(tǒng)建模 %

43、創(chuàng)建不同增益條件下的三個(gè)傳送函數(shù) ifn=1; opt=SensitivitytoGain fori=1:3; K=input(Inputgaintypical=25V/radian:); d2=Ra/La+Beq/Jeq; d3=(Ra*Beq+kt*kv)/(La*Jeq); d4=(K*kt)/(La*Jeq); NUM(i,:)=(K*kt)/(La*Jeq); DEN(i,:)=1d2d3d4; end;第3章 控制系統(tǒng)建模 % 創(chuàng)建不同機(jī)械阻尼條件下的三個(gè)傳送函數(shù) elseifn=2; opt=SensitivitytoDamping fori=1:3; Beq=input(Inp

44、utdampingfactortypical=5.0e-3Nms:); d2=Ra/La+Beq/Jeq; d3=(Ra*Beq+kt*kv)/(La*Jeq); d4=(K*kt)/(La*Jeq); NUM(i,:)=(K*kt)/(La*Jeq); DEN(i,:)=1d2d3d4; end;第3章 控制系統(tǒng)建模% 在不同轉(zhuǎn)動(dòng)慣量條件下創(chuàng)建三個(gè)傳送函數(shù) elseifn=3; opt=SensitivitytoInertia fori=1:3; Jeq=input(Inputinertiatypical=30.0e-6kgm2:); d2=Ra/La+Beq/Jeq; d3=(Ra*Beq

45、+kt*kv)/(La*Jeq); d4=(K*kt)/(La*Jeq); NUM(i,:)=(K*kt)/(La*Jeq); DEN(i,:)=1d2d3d4; end;第3章 控制系統(tǒng)建模 end; sys1=tf(NUM(1,:),DEN(1,:); sys2=tf(NUM(2,:),DEN(2,:); sys3=tf(NUM(3,:),DEN(3,:);% 繪制Bode圖和階躍呼應(yīng)曲線 w1=input(Selectbeginningfreq.exponent(rad/sec)1:); ifisempty(w1); w1=1; end; w2=input(Selectendfreq.e

46、xponent(rad/sec)3:); ifisempty(w2); w2=3; end; w=logspace(w1,w2,100); figure(2)第3章 控制系統(tǒng)建模 bode(sys1,sys2,sys3,w) % 閉環(huán)系統(tǒng)的Bode圖 title(FrequencyResponse-, opt); gtext(low),gtext(mid),gtext(high); gtext(low),gtext(mid),gtext(high);% 閉環(huán)系統(tǒng)的階躍呼應(yīng) t2=input(Specifyendtime(sec)0.08:); ifisempty(t2); t2=0.08; end; t=linspace(0,t2,101); xs(:,1)=step(sys1,t);第3章 控制系統(tǒng)建模 xs(:,2)=step(sys2,t); xs(:,3)=step(sys3,t); figure(3) plot(t,xs(:,1),r-,t,xs(:,2),b-,t,xs(:,3),g:),grid xlabel(time(sec),ylabel(position(radians); title(StepResponse-, opt)