弹簧质量阻尼系统模型

自动控制原理综合训练项目题目：关于MSD系统控制的设计目录1设计任务及要求分析................................................................................................31.1初始条件.............................................................................................................................31.2要求完成的任务.................................................................................................................31.3任务分析.............................................................................................................................42系统分析及传递函数求解........................................................................................42.1系统受力分析.....................................................................................................................42.2传递函数求解....................................................................................................................92.3系统开环传递函数的求解.................................................................................................93.用MATLAB对系统作开环频域分析.......................................................................103.1开环系统波特图...............................................................................................................103.2开环系统奈奎斯特图及稳定性判断...............................................................................124.系统开环频率特性各项指标的计算......................................................................14总结..............................................................................................................................17参考文献......................................................................................................................18弹簧-质量-阻尼器系统建模与频率特性分析1设计任务及要求分析1.1初始条件已知机械系统如图。2b1kyp2kx图1.1机械系统图1.2要求完成的任务（1）推导传递函数)(/)(sXsY，)(/)(sPsX，m（2）给定mNkmNkmsNbgm/5,/8,/6.0,2.0212，以p为输入)(tu（3）用Matlab画出开环系统的波特图和奈奎斯特图，并用奈奎斯特判据分析系统的稳定性。（4）求出开环系统的截止频率、相角裕度和幅值裕度。（5）对上述任务写出完整的课程设计说明书，说明书中必须进行原理分析，写清楚分析计算的过程及其比较分析的结果，并包含Matlab源程序或Simulink仿真模型，说明书的格式按照教务处标准书写。1.3任务分析由初始条件和要求完成的主要任务，首先对给出的机械系统进行受力分析，列出相关的微分方程，对微分方程做拉普拉斯变换，将初始条件中给定的数据代入，即可得出)(/)(sXsY，)(/)(sPsX两个传递函数。由于本系统是一个单位负反馈系统，故求出的传递函数即为开环传函。后在MATLAB中画出开环波特图和奈奎斯特图，由波特图分析系统的频率特性，并根据奈奎斯特判据判断闭环系统位于右半平面的极点数，由此可以分析出系统的稳定性。最后再计算出系统的截止频率、相角裕度和幅值裕度，并进一步分析其稳定性能。2系统分析及传递函数求解2.1系统受力分析单自由度有阻尼振系的力学模型如图2-1所示，包括弹簧、质量及阻尼器。以物体的平衡位置0为原点，建立图示坐标轴x。则物体运动微分方程为kxxcxm－＝－（2-1）式中：xc为阻尼力，负号表示阻尼力方向与速度方向相反。图2-1将上式写成标准形式，为0kxxcxm(2-2）令p2=mk,mcn2,则上式可简化为022pxnx(2-3)这就是有阻尼自由振动微分方程。它的解可取stex，其中s是待定常数。代入（2-1）式得0)2(22stepnss，要使所有时间内上式都能满足，必须0222pnss，此即微分方程的特征方程，其解为222,1pnns（2-4）于是微分方程（2-1）的通解为)(2222212121tpntpnnttstsececeececx（2-5）式中待定常数c1与c2决定与振动的初始条件。振动系统的性质决定于根式22pn是实数、零、还是虚数。对应的根s1与s2可以是不相等的负实根、相等的负实根或复根。若s1与s2为等根时，此时的阻尼系数值称之为临界阻尼系数，记为cc，即cc＝2mp。引进一个无量纲的量，称为相对阻尼系数或阻尼比。cccmpcpn/2//（2-6）当np或1,根式22pn是实数，称为过阻尼状态，当np或1，根式22pn是虚数，称为弱阻尼状态，当n＝p，即＝1，称为临界阻尼状态。现分别讨论三种状态下的运动特性。1.过阻尼状态此时1，即22pnn，（b）式中s1及s2均为负值，则tse1及tse2是两根下降的指数曲线，故（2-2）式所表示的是两条指数曲线之和，仍按指数衰减，不是振动。图3-2所示为c1c2,c10时的情况。图2-22.临界阻尼状态此时＝1，（b）式中s1＝s2＝－n＝－p，特征方程的根是重根，方程（2-1）的另一解将为te－pt，故微分方程（2-1）的通解为x＝（c1＋c2t）e－pt（2-7）式中等号右边第一项c1e－pt是一根下降的指数曲线，第二项则可应用麦克劳林级数展开成以下形式：!/!3/!2//12322/22ntptptpptcectecnntptpt（2-8）从上式看出，当时间t增长时，第二项c2te－pt也趋近于零。因此（c）式表示的运动也不是振动，也是一个逐渐回到平衡位置的非周期运动。3.弱阻尼状态此时pn,或1。利用欧拉公式tnpitnpeetnptpn2222sincos2222（2-9）可将（2-2）式改写为)sincos()(222221212222tnpDtnpDeeCeCexnttnpitnpint（2-10）或)sin(22tnpAexnt（1-11）令22nppd，则)sin(tpAexdnt（2-12）式中A与为待定常数，决定于初始条件。设t＝0时，x＝x0，0xx，则可求得000120020,)(xnxxpxtgpnxxAdd（2-13）将A与代入（2-4）式，即可求得系统对初始条件的响应，由式（2-13）可知，系统振动已不再是等幅的简谐振动，而是振幅被限制在曲线ntAe之内随时间不断衰减的衰减振动。如图3-3所示。图2-3这种衰减振动的固有圆频率、固有频率和周期分别为2221PnPPd（2-14）222222221111221122TPnPTfPnPfdd（2-15）式中P、f、T是无阻尼自由振动的固有圆频率、固有频率和周期。由上可见，阻尼对自由振动的影响有两个方面：一方面是阻尼使自由振动的周期增大、频率减小，但在一般工程问题中n都比P小得多，属于小阻尼的情况。例=n/p=0.05时，fd=0.9990f，Td=1.00125T；而在=0.20时，fd=0.98f，Td=1.02T，所以在阻尼比较小时，阻尼对系统的固有频率和周期的影响可以略去不计，即可以近似地认为有阻尼自由振动的频率和周期与无阻尼自由振动的频率和周期相等。另一方面，阻尼对于系统振动振幅的影响非常显著，阻尼使振幅随着时间不断衰减，其顺次各个振幅是：t=t1时，A1=Ae-nt1；t=t1+Td时，A2=A)(1dTtne；t=t1+2Td时，A3=A)2(1dTtne，…..。而相邻两振幅之比是个常数。即nTdjjeAA1/（2-16）式中η称为减幅系数或振幅衰减率，n称为衰减系数，n越大表示阻尼越大，振幅衰减也越快。当＝0.05时，η＝1.37，A2=A1/1.37=0.73A1，每一个周期内振幅减少27％，振幅按几何级数衰减，经过10次振动后，振幅将减小到初值的4.3％。可见，衰减是非常显著的。在工程上，通常取（2-6）式的自然对数以避免取指数的不便，即djjnTAALn)/(1（2-17）式中δ称为对数减幅或对数衰减率。将22/2npTd代入，得2221/2/2npn（2-18）当1时，δ≈2π（2-19）因为任意两个相邻的振幅之比是一个常数enTd，即eeAAAAAAAAnTdjj1433221/......///故有)/)......(/)(/(/1322111jjjAAAAAAAAje因此对数减幅δ也可表达为)1(11jAALnj（2-20）此外，根据（3-6）式，可以用实测法来求得系统的阻尼系数。因为111121jjdjjddjjAALnTmcAALnTnnTAALn故12jjdAALnTmc（2-21）所以只要实测得出衰减振动的周期Td及相邻两次振幅Aj和Aj+1，即可计算出系统的阻尼系数C。根据弹簧和阻尼器的特性可得以下关系式：Fk1(t)=k1x(t)，Fk2(t)=k2[x(t)－y(t)]，Fb2(t)=b2dy(t)/dt设不加p(t)时，质量块处于平衡状态，此时x=0，y=0，即x(0)=0，y(0)=0，根据受力平衡方程，在不计重力时，可得出以下方程：k2[x(t)-y(t)]=b2dy(t)/dt(2-22)又根据牛顿第二定律，有方程：md2x(t)/dt2=p(t)－Fk1(t)－Fk2(t)－Fb2(t)(2-23)2.2传递函数求解（1）求Y(s)/X(s)：对式(2-1)进行拉普拉斯变换，得：k2X(s)－k2Y(s)=b2*sY(s)，化简得传递函数：Y(s)/X(s)=k2/(b2s+k2)(2-24)（2）求X(s)/P(s)：对式(