弹簧质量阻尼系统模型自动控制原理综合训练项目题目:关于MSD系统控制的设计目录1设计任务及要求分析 (5)1.1初始条件 (5)1.2要求完成的任务 (6)1.3任务分析 (7)2系统分析及传递函数求解 (7)2.1系统受力分析 (7)2.2传递函数求解 (13)2.3系统开环传递函数的求解 (13)3.用MATLAB寸系统作开环频域分析 (14)3.1开环系统波特图 (14)3.2开环系统奈奎斯特图及稳定性判断.164.系统开环频率特性各项指标的计算 (19)总结 (21)参考文献 (22)弹簧-质量-阻尼器系统建模与频率特性分析1设计任务及要求分析1.1初始条件已知机械系统如图。
//////////--- NN\/V-------b2NWk iIP k2x图 1.1 机械系统图1.2要求完成的任务(1)推导传递函数 Y(s)/X(s) , X(s)/P(s),(2)给定m=0.2g,b2 =0.6N-s/mK =8N/m,k2 =5N/m,以p 为输入u(t)(3)用Matlab画出开环系统的波特图和奈奎斯特图,并用奈奎斯特判据分析系统的稳定性。
(4)求出开环系统的截止频率、相角裕度和幅值裕度。
(5)对上述任务写出完整的课程设计说明书,说明书中必须进行原理分析,写清楚分析计算的过程及其比较分析的结果,并包含Matlab源程序或Simulink仿真模型,说明书的格式按照教务处标准书写1.3任务分析由初始条件和要求完成的主要任务,首先对给出的机械系统进行受力分析,列出相关的微分方程,对微分方程做拉普拉斯变换,将初始条件中给定的数据代入,即可得出 Y(s)/X(s),X(s)/P(s)两个传递函数。
由于本系统是一个单位负反馈系统,故求出的传递函数即为开环传函。
后在MATLAB中画出开环波特图和奈奎斯特图,由波特图分析系统的频率特性,并根据奈奎斯特判据判断闭环系统位于右半平面的极点数,由此可以分析出系统的稳定性。
最后再计算出系统的截止频率、相角裕度和幅值裕度,并进一步分析其稳定性能。
2系统分析及传递函数求解2.1系统受力分析单自由度有阻尼振系的力学模型如图2-1所示,包括弹簧、质量及阻尼器。
以物体的平衡位置0为原点,建立图示坐标轴X。
则物体运动微分方程为mx = —ex —kx (2-1)式中:-CX为阻尼力,负号表示阻尼力方向与速度方向相反。
将上式写成标准形式,为令p2= —, 2n = C,贝U上式可简化为m m2x 2nx p 0 (2-3)这就是有阻尼自由振动微分方程。
它的解可取x = e st,其中s是待定常数。
代入(2-1)式得(s2• 2ns p2)e st = 0,要使所有时间内上式都能满足,必须s2 2ns • p2 =0,此即微分方程的特征方程,其解为3,2 二-n 一. n2- p2(2-4) 于是微分方程(2-1 )的通解为. . . 2 2 2 2.x = C i e° +C2e s2=e』(C i e'n " +c2n ) ( 2-5)式中待定常数C1与C2决定与振动的初始条件。
振动系统的性质决定于根式■ n2二p2是实数、零、还是虚数。
对应的根s与S2可以是不相等的负实根、相等的负实根或复根。
若S1与S2为等根时,此时的阻尼系数值称之为临界阻尼系数,记为c c,即e c= 2mp。
引进一个无量纲的量,称为相对阻尼系数或阻尼比。
=n/ p = c/2mp = c/c C(2-6)当n>p或'>1,根式i n2- p2是实数,称为过阻尼状态,当n<p或<1,根式:.n2 -p2是虚数,称为弱阻尼状态,当n = p,即•二1,称为临界阻尼状态。
现分别讨论三种状态下的运动特性。
1■过阻尼状态此时【>1,即Jn2-p2<n,(b)式中S1及S2均为负值,则e sit及e s2t是两根下降的指数曲线,故(2-2)式所表示的是两条指数曲线之和,仍按指数衰减,不是振动。
图3-2所示为c i >C 2,c i <0时的情况。
图2-22■临界阻尼状态此时 =1, (b )式中s i = sz =-n =— p ,特征方程的根是重根,方程(2-1) 的另一解将为te —叭故微分方程(2-1)的通解为x =( c i + C 2t ) e pt(2-7)式中等号右边第一项c i e -pt 是一根下降的指数曲线,第二项则可应用麦克劳林级 数展开成以下形式:从上式看出,当时间t 增长时,第二项c 2te — pt也趋近于零。
因此(c )式表示的 运动也不是振动,也是一个逐渐回到平衡位置的非周期运动。
3■弱阻尼状态此时p>n,或<1 o 利用欧拉公式二cos., p 2 - n 2t 士i sin .. p 2 - n 2t (2-9)可将(2-2)式改写为x = Ae" sin(.. p 2 _ n 2t)令 P d — p 2 - n 2,则C 2te 4 C 2~pTTte________________ C2 _________________1/t p p 2t/2! p 3t 2 /3L ...... - p n t n/n! (2-8) ■ nt二 e(Gei p2-n 2t. 2 2.■ C 2e^ pe Jnt(D i cos . p 2「n 2t D 2 sin . p 2「n 2t)(2-10)(1-11)式中A 与「为待定常数,决定于初始条件。
设 t = 0时,x = x o, x = x o ,则可求 得A 品(30心赵」厘乩(2-13)V P dXo + nx o将A 与「代入(2-4)式,即可求得系统对初始条件的响应,由式(2-13)可知, 系统振动已不再是等幅的简谐振动,而是振幅被限制在曲线_Ae 』之内随时间不 断衰减的衰减振动。
如图3-3所示。
图2-3这种衰减振动的固有圆频率、固有频率和周期分别为P d = 、P 2— n 2二 PJ — 2( 2-14)(2-15)2 2 二 _ 2_ ____ 二 2 二 1 二 丁 1—P 21 n 2「P 、一1- 2「1- 2式中P、f、T是无阻尼自由振动的固有圆频率、固有频率和周期由上可见,阻尼对自由振动的影响有两个方面:一方面是阻尼使自由振动的周期增大、频率减小,但在一般工程问题中n都比P小得多,属于小阻尼的情况。
例=n/p=0.05 时,f d=0.9990f, T d=1.00125T;而在=0.20 时,f d=0.98f, T d=1.02T,所以在阻尼比较小时,阻尼对系统的固有频率和周期的影响可以略去不计,即可以近似地认为有阻尼自由振动的频率和周期与无阻尼自由振动的频率和周期相等。
另一方面,阻尼对于系统振动振幅的影响非常显著,阻尼使振幅随着时间不断衰减,其顺次各个振幅是:t=t i时,A i=Ae-nt i ;t=t i+T d时,A2=A e』(tlTd) ; t=t i+2T d时,A3=A e』(t「2Td),…..。
而相邻两振幅之比是个常数。
即(2-i6)式中n称为减幅系数或振幅衰减率,n称为衰减系数,n越大表示阻尼越大,振幅衰减也越快。
当'=0.05时,n= i・37, A2=A i/i.37=0.73A i,每一个周期内振幅减少27%,振幅按几何级数衰减,经过i0次振动后,振幅将减小到初值的4.3%。
可见,衰减是非常显著的。
在工程上,通常取(2-6)式的自然对数以避免取指数的不便,即-=Ln(A j / A j = nT d 式中S称为对数减幅或对数衰减率。
将T d =2 / ... p2-n2代入,得当<<1时,(2-17):=2”p2 -n2=2^/..1- 2(2-18)A / A ? = A 2 / A 3 = A 3 / A 4 = .......... = A j / A j 1 = e = e故有A/A H =(A/A 2)(A 2/A 3)……(A j /A® =e 诸因此对数减幅S 也可表达为A A (j 4)此外,根据(3-6)式,可以用实测法来求得系统的阻尼系数。
因为c 旦Ln 旦T d A j *所以只要实测得出衰减振动的周期 T d 及相邻两次振幅A j 和A j+i ,即可计算出系 统的阻尼系数C 。
根据弹簧和阻尼器的特性可得以下关系式:F ki (t)=k i x(t),F k2(t)=k 2[x(t) — y(t)],F b2(t)=b 2dy(t)/dt设不加p(t)时,质量块处于平衡状态,此时 x=0, y=0,即x(0)=0, y(0)=0, 根据受力平衡方程,在不计重力时,可得出以下方程:k 2[x(t)-y(t)]=b 2dy(t)/dt(2-22)(2-20)A jLnnT d -■A j *n 二丄 LnT dc 2m二丄Ln A j TdAj 1(2-21)又根据牛顿第二定律,有方程:2 2md2x(t)/dt2=p(t) —F ki(t) —F k2(t) —F b2(t) (2-23)2.2传递函数求解(1)求Y(s)/X(s):对式(2-1)进行拉普拉斯变换,得:k2X(s)—k2Y(s)=b2*sY(s),化简得传递函数:Y(s)/X(s)=k 2/(b2s+k2) (2-24) (2)求X(s)/P(s):一2对式(2-2)进行拉普拉斯变换,得:ms X(s)=P(s)—k i X(s) —2k2[X(s) —Y(s)],并将式(2-3)代入可解得传递函数:3 2X(s)/P(s)=(b2s+k2)/[mb2s3+mk2s2+b2(k i+2k2)s+k i k2] (2-25)已知条件为:给定m 二0.2g, b2 = 0.6N *s/m, k^ 8N /m,k2 = 5N /m,设 p(t)是输入 u(t)的阶跃力。
将所给参数代入传递函数式(2-3)和式(2-4)中,可求得具体的传递函数如下:Y(s)/X(s)=5/(0.6s+5) (2-26) X(s)/P(s)=(0.6s+5)/ (1.2*10A-4s3+10A-3s2+10.8s+4Q (2-27) 2.3系统开环传递函数的求解(1)对于Y(s)/X(s):由微分方程Y(s)/X(s)=5/(0.6s+5)可画出单位负反馈系统方框结构图如下:2(s) X (s)5/(0.6s+5) —故开环传递函数为:G (S) =5/(0.6s+5)(2)对于X(s)/P(s):由微分方程ms2X(s)=P(s) —k1X(s) —2k2[X(s) —Y(s)]及Y(s)/X(s)=k2/(b2s+k2)可画出系统方框结构图如下:P(s) X(s) 故开环传递G(s)=3.用MATLAB寸系统作开环频域分析3.1开环系统波特图(1)对于Y(s)/X(s):(画波特图时采用的MATLAB语句如下:>> num=[5];den=([0.6,5]);>> margin(num,den)特性运行结果如图3-1%画系统的开环对数幅频、相频10 10"Frequency (rad/B)图3-1 Y(s)/X(s)的开环波特图(2)对于 X(s)/P(s):G(s)=画波特图时采用的MATLAB 语句如下:>> num=[0.6,5];den=([]);>> margin(num,den)特性运行结果如图3-2所示:GB MC JiagremGm« iif ,Pm - -1 BO deg fat 0 rad/s}101J10%画系统的开环对数幅频、相频3od« DiagnfnGrr - Inf dB (at hf radJs) , Pn- 1&7de-g 但t3唱ra岀s)(S4pjesgulxio T101w' iMiaFrequency (ra^s)图3-2 X(s)/P(s)的开环波特图3.2开环系统奈奎斯特图及稳定性判断(1)对于Y(s)/X(s)画奈奎斯特图时MATLAB语句如下:>> num=[5];>> den=[0.6,5];>> nyquist(num,den)运行结果如图3-3所示:图 3-3Y (s ) /X (s )开环奈奎斯特图开环传函(),由于系统开环传递函数不存在右半平面的 极点,故p=0,从 变到 时,系统的开环幅相曲线不能包围(-1,j0)点周数N=0,则系统位于右半平面的闭环极点数为:Z=P-2N=0,故系统是稳定的。