彈簧質(zhì)量阻尼系統(tǒng)模型.doc

自動(dòng)控制原理綜合訓(xùn)練項(xiàng)目題目：關(guān)于MSD系統(tǒng)控制的設(shè)計(jì) 目 錄 1設(shè)計(jì)任務(wù)及要求分析31.1初始條件31.2要求完成的任務(wù)31.3任務(wù)分析42系統(tǒng)分析及傳遞函數(shù)求解42.1系統(tǒng)受力分析42.2 傳遞函數(shù)求解92.3系統(tǒng)開環(huán)傳遞函數(shù)的求解93.用MATLAB對(duì)系統(tǒng)作開環(huán)頻域分析103.1開環(huán)系統(tǒng)波特圖103.2開環(huán)系統(tǒng)奈奎斯特圖及穩(wěn)定性判斷124.系統(tǒng)開環(huán)頻率特性各項(xiàng)指標(biāo)的計(jì)算14總結(jié)17參考文獻(xiàn)18彈簧-質(zhì)量-阻尼器系統(tǒng)建模與頻率特性分析1設(shè)計(jì)任務(wù)及要求分析1.1初始條件 已知機(jī)械系統(tǒng)如圖。 yp m x 圖1.1 機(jī)械系統(tǒng)圖 1.2要求完成的任務(wù)（1） 推導(dǎo)傳遞函數(shù)，（2） 給定，以p為輸入（3） 用Matlab畫出開環(huán)系統(tǒng)的波特圖和奈奎斯特圖，并用奈奎斯特判據(jù)分析系統(tǒng)的穩(wěn)定性。（4） 求出開環(huán)系統(tǒng)的截止頻率、相角裕度和幅值裕度。（5） 對(duì)上述任務(wù)寫出完整的課程設(shè)計(jì)說明書，說明書中必須進(jìn)行原理分析，寫清楚分析計(jì)算的過程及其比較分析的結(jié)果，并包含Matlab源程序或Simulink仿真模型，說明書的格式按照教務(wù)處標(biāo)準(zhǔn)書寫。1.3任務(wù)分析由初始條件和要求完成的主要任務(wù)，首先對(duì)給出的機(jī)械系統(tǒng)進(jìn)行受力分析，列出相關(guān)的微分方程，對(duì)微分方程做拉普拉斯變換，將初始條件中給定的數(shù)據(jù)代入，即可得出，兩個(gè)傳遞函數(shù)。由于本系統(tǒng)是一個(gè)單位負(fù)反饋系統(tǒng)，故求出的傳遞函數(shù)即為開環(huán)傳函。后在MATLAB中畫出開環(huán)波特圖和奈奎斯特圖，由波特圖分析系統(tǒng)的頻率特性，并根據(jù)奈奎斯特判據(jù)判斷閉環(huán)系統(tǒng)位于右半平面的極點(diǎn)數(shù)，由此可以分析出系統(tǒng)的穩(wěn)定性。最后再計(jì)算出系統(tǒng)的截止頻率、相角裕度和幅值裕度，并進(jìn)一步分析其穩(wěn)定性能。2系統(tǒng)分析及傳遞函數(shù)求解2.1系統(tǒng)受力分析 單自由度有阻尼振系的力學(xué)模型如圖2-1所示，包括彈簧、質(zhì)量及阻尼器。以物體的平衡位置0為原點(diǎn)，建立圖示坐標(biāo)軸x。則物體運(yùn)動(dòng)微分方程為 （2-1）式中 ： 為阻尼力，負(fù)號(hào)表示阻尼力方向與速度方向相反。圖2-1將上式寫成標(biāo)準(zhǔn)形式，為 (2-2）令p2=, , 則上式可簡(jiǎn)化為 (2-3)這就是有阻尼自由振動(dòng)微分方程。它的解可取，其中s是待定常數(shù)。代入（2-1）式得 ，要使所有時(shí)間內(nèi)上式都能滿足，必須，此即微分方程的特征方程，其解為 （2-4）于是微分方程（2-1）的通解為 （2-5）式中待定常數(shù)c1與c2決定與振動(dòng)的初始條件。振動(dòng)系統(tǒng)的性質(zhì)決定于根式是實(shí)數(shù)、零、還是虛數(shù)。對(duì)應(yīng)的根s1與s2可以是不相等的負(fù)實(shí)根、相等的負(fù)實(shí)根或復(fù)根。若s1與s2為等根時(shí)，此時(shí)的阻尼系數(shù)值稱之為臨界阻尼系數(shù)，記為cc，即cc2mp。引進(jìn)一個(gè)無量綱的量，稱為相對(duì)阻尼系數(shù)或阻尼比。 （2-6） 當(dāng)np或1,根式是實(shí)數(shù)，稱為過阻尼狀態(tài)，當(dāng)np或1，即c2,c1n,或1。利用歐拉公式（2-9）可將（2-2）式改寫為（2-10）或 （1-11） 令，則 （2-12） 式中A與為待定常數(shù)，決定于初始條件。設(shè)t0時(shí)，xx0，則可求得 （2-13）將A與代入（2-4）式，即可求得系統(tǒng)對(duì)初始條件的響應(yīng)，由式（2-13）可知，系統(tǒng)振動(dòng)已不再是等幅的簡(jiǎn)諧振動(dòng)，而是振幅被限制在曲線之內(nèi)隨時(shí)間不斷衰減的衰減振動(dòng)。如圖3-3所示。 圖2-3這種衰減振動(dòng)的固有圓頻率、固有頻率和周期分別為（2-14）（2-15）式中P、f、T是無阻尼自由振動(dòng)的固有圓頻率、固有頻率和周期。由上可見，阻尼對(duì)自由振動(dòng)的影響有兩個(gè)方面：一方面是阻尼使自由振動(dòng)的周期增大、頻率減小，但在一般工程問題中n都比P小得多，屬于小阻尼的情況。例=n/p=0.05時(shí)，fd=0.9990f，Td=1.00125T；而在=0.20時(shí)，fd=0.98f，Td=1.02T，所以在阻尼比較小時(shí)，阻尼對(duì)系統(tǒng)的固有頻率和周期的影響可以略去不計(jì)，即可以近似地認(rèn)為有阻尼自由振動(dòng)的頻率和周期與無阻尼自由振動(dòng)的頻率和周期相等。另一方面，阻尼對(duì)于系統(tǒng)振動(dòng)振幅的影響非常顯著，阻尼使振幅隨著時(shí)間不斷衰減，其順次各個(gè)振幅是：t=t1時(shí)，A1=Ae-nt1；t=t1+Td時(shí)，A2=A；t=t1+2Td時(shí)，A3=A，.。而相鄰兩振幅之比是個(gè)常數(shù)。即 （2-16）式中稱為減幅系數(shù)或振幅衰減率，n稱為衰減系數(shù)，n越大表示阻尼越大，振幅衰減也越快。當(dāng)0.05時(shí)，1.37，A2=A1/1.37=0.73A1，每一個(gè)周期內(nèi)振幅減少27，振幅按幾何級(jí)數(shù)衰減，經(jīng)過10次振動(dòng)后，振幅將減小到初值的4.3?？梢?，衰減是非常顯著的。在工程上，通常?。?-6）式的自然對(duì)數(shù)以避免取指數(shù)的不便，即 （2-17）式中稱為對(duì)數(shù)減幅或?qū)?shù)衰減率。將代入，得 （2-18）當(dāng) num=5;den=(0.6,5); margin(num,den) %畫系統(tǒng)的開環(huán)對(duì)數(shù)幅頻、相頻特性運(yùn)行結(jié)果如圖3-1 圖3-1 Y(s)/X(s)的開環(huán)波特圖（2）對(duì)于X(s)/P(s):G(s)= (0.6s+5)/(1.2*10-4s3+10-3s2+10.8s+40) 畫波特圖時(shí)采用的MATLAB語句如下: num=0.6,5;den=(1.2*10-4,10-3,10.8,40); margin(num,den) %畫系統(tǒng)的開環(huán)對(duì)數(shù)幅頻、相頻特性運(yùn)行結(jié)果如圖3-2所示: 圖3-2 X(s)/P(s)的開環(huán)波特圖3.2開環(huán)系統(tǒng)奈奎斯特圖及穩(wěn)定性判斷（1）對(duì)于Y(s)/X(s)畫奈奎斯特圖時(shí)MATLAB語句如下： num=5; den=0.6,5; nyquist(num,den)運(yùn)行結(jié)果如圖3-3所示：圖3-3 Y（s）/X（s）開環(huán)奈奎斯特圖開環(huán)傳函G（S）=5/(0.6s+5),由于系統(tǒng)開環(huán)傳遞函數(shù)不存在右半平面的極點(diǎn)，故P=0，從0變到+時(shí)，系統(tǒng)的開環(huán)幅相曲線不能包圍（-1，j0）點(diǎn)周數(shù)N=0，則系統(tǒng)位于右半平面的閉環(huán)極點(diǎn)數(shù)為:Z=P-2N=0,故系統(tǒng)是穩(wěn)定的。（2）對(duì)于X(s)/P(s)畫奈奎斯特圖時(shí)MATLAB語句如下： num=0.6,5; den=1.2*10-4,10-3,10.8,40; nyquist(num,den)運(yùn)行結(jié)果如圖3-4所示：圖3-4 X（s）/P（s）開環(huán)奈奎斯特圖開環(huán)傳函G(s)= (0.6s+5)/(1.2*10-4s3+10-3s2+10.8s+40) ，由于系統(tǒng)開環(huán)傳遞函數(shù)不存在右半平面的極點(diǎn)，故P=0，從0變到+時(shí)，系統(tǒng)的開環(huán)幅相曲線不能包圍（-1，j0）點(diǎn)周數(shù)N=0，則系統(tǒng)位于右半平面的閉環(huán)極點(diǎn)數(shù)為:Z=P-2N=0,故系統(tǒng)是穩(wěn)定的。4.系統(tǒng)開環(huán)頻率特性各項(xiàng)指標(biāo)的計(jì)算（1）對(duì)于Y(s)/X(s)：G（S）=5/(0.6s+5)計(jì)算各項(xiàng)頻率指標(biāo)時(shí)采用的MATLAB語句如下: num=5;den=(0.6,5); margin(num,den); gm,pm,wcg,wcp=margin(num,den)計(jì)算幅值裕度gm(0)、相位裕度pm（h0）、穿越頻率wcg(r0)、截止頻率wcp(c0)。運(yùn)行結(jié)果gm = Infpm = 180wcg = NaNwcp = 0 由結(jié)果可知該系統(tǒng)幅值裕度為無窮，截止頻率為0，相位裕度為180是正值，故系統(tǒng)穩(wěn)定。（2）對(duì)于X(s)/P(s) :G(s)= (0.6s+5)/(1.2*10-4s3+10-3s2+10.8s+40) 計(jì)算各項(xiàng)頻率指標(biāo)時(shí)采用的MATLAB語句如下: num=0.6,5;den=(1.2*10-4,10-3,10.8,40) margin(num,den); gm,pm,wcg,wcp=margin(num,den)計(jì)算幅值裕度gm(0)、相位裕度pm（h0）、穿越頻率wcg(r0)、截止頻率wcp(c0)。運(yùn)行結(jié)果gm = Infpm = 15.6933wcg = Infwcp = 307.8588 由結(jié)果可知該系統(tǒng)幅值裕度為無窮，截止頻率為308rad/s，相位裕度為15.7是正值，故系統(tǒng)穩(wěn)定。總結(jié)本次課設(shè)是對(duì)一個(gè)彈簧-質(zhì)量-阻尼器系統(tǒng)建模并進(jìn)行頻率特性分析。首先根據(jù)這個(gè)實(shí)際的機(jī)械系統(tǒng)的受力分析得出它的受力微分方程，對(duì)其進(jìn)行拉普拉斯變換，可以得出傳遞函數(shù)。在求開環(huán)傳遞函數(shù)的過程中我遇到了一些困難，在老師的指點(diǎn)和同學(xué)的幫助下我發(fā)現(xiàn)自己其實(shí)把問題想得過于復(fù)雜了，原來這是一個(gè)單位負(fù)反饋的穩(wěn)定系統(tǒng)，求出的傳遞函數(shù)即為開環(huán)傳函。接下來便是MATLAB的應(yīng)用，利用MATLAB可以輕松地對(duì)系統(tǒng)做出頻率特性分析，畫出Bode圖和奈奎斯特圖，并通過奈奎斯特判據(jù)判斷閉環(huán)系統(tǒng)的穩(wěn)定性。也可以利用MATLAB語句直接求出各項(xiàng)頻率特性指標(biāo)，從而可以進(jìn)一步對(duì)系統(tǒng)做出分析，完成既定目標(biāo)。通過本次課設(shè)，加強(qiáng)了我對(duì)MATLAB程序的應(yīng)用能力，這是一款功能強(qiáng)大而又實(shí)用性很強(qiáng)的程序，對(duì)于我們專業(yè)的學(xué)習(xí)有著很強(qiáng)的幫助性；另一方面也加強(qiáng)了我對(duì)課本理論知識(shí)的理解，通過MATLAB的分析也印證了平時(shí)自己學(xué)習(xí)理論知識(shí)時(shí)所用分析方法的正確性。最后通過本次課設(shè)也提高了我個(gè)人獨(dú)立思考、查閱資料和解決問題的能力，使我受益匪淺。參考文獻(xiàn)1 王子才. 控制系統(tǒng)設(shè)計(jì)手冊(cè). 北京：國(guó)防工業(yè)出版社，19932 王樹青，樂嘉謙. 自動(dòng)化與儀表工程師手冊(cè). 北京：化學(xué)工業(yè)出版社