弹簧-高质量-阻尼模型

1、弹簧-质量-阻尼系统1研究背景及意义弹簧-质量-阻尼系统是一种比较普遍的机械振动系统，研究这种系统对于我们的生活与科技也是具有意义的，生活中也随处可见这种系统,例如汽车缓冲器就是一种可以耗减运动能量的装置，是保证驾驶员行车安全的必备装置,再者在建筑抗震加固措施中引入阻尼器，改变结构的自振特性，增加结构阻尼，吸收地震能量,降低地震作用对建筑物的影响。因此研究弹簧-质量-阻尼结构是很具有现实意义。2弹簧-质量-阻尼模型的建立数学模型是定量地描述系统的动态特性，揭示系统的结构、参数与动态特性之间关系的数学表达式。其中，微分方程是基本的数学模型，不论是机械的、液压的、电气的或热力学的系统等都可以用微分

2、方程来描述。微分方程的解就是系统在输入作用下的输出响应。所以，建立数学模型是研究系统、预测其动态响应的前提。通常情况下，列写机械振动系统的微分方程都是应用力学中的牛顿定律、质量守恒定律等。弹簧-质量-阻尼系统是最常见的机械振动系统。机械系统如图2.1所示,图2.1弹簧-质量-阻尼系统简图其中mi，m2表示小车的质量，Ci表示缓冲器的粘滞摩擦系数，«表示弹簧的弹性系数，Fi（t）表示小车所受的外力，是系统的输入即Uj（t）= Fi （t） ,Xj（t）表示小车的位移，是系统的输出，即Yi（t）= Xi（t），i=1,2。设缓冲器的摩擦力与活塞的速度成正比，其中 m1 =1kg ， m2

3、 =2kg ， k1 = k3=100N/cm， k2 =300N/cm ， c1 = c3 =3N ? s/cm ，C2=6N ? s/cm。由图2.1，根据牛顿第二定律”建立系统的动力学模型如下：对mi有：IJ >_< 卜(2-1 )对m2有：I , f _>_< L ' '(2-2 )f Cj + Cn + k + 1c n kpSC? = Fj(t)i円屉 + (c2 + g3)x2 一 c2kj + (k2 +k3)x3 - kax!=吟(t)3建立状态空间表达式U2(t) k2Xi (ki k2)X2 (I312)X4l2X3(2-4)令X3

4、 &X4 X2, UiFi,U2 F2,则原式可化为：mi& m2X4(li(I212) X313) X4I2X412X3(ki(k2k2)Xik3)X2k2X2k2XiUi(t)U2(t)化简得：& ui(t)k2X2(kik2)Xi(liI2)X3I2X4（2-3）mim2整理得:00i000000iXi00X(ki k2)k2(liI2)I2X2iui0 imimimmi%miU2k2(k3 k2)I2(I3 I2)X40im2m2m2m2m2Xi1 0 0 0X20 1 0 0X3X4001000000100x4003009x61u01520034.500.51

5、000yx0100(2-5 )叶 1,m22,Kk3100,k2300l1l33,l26001000代入数据得：A40030091502003则系统的状态空间表达式为0001001000BC61001004.500.54化为对角标准型当系统矩阵A有n个不相等的特征根i(i1,2,3.)时，相应的有n个不相等的特征向量mi m2 m3 mkmi(i1,2,3.)，所以有矩阵 A的特征矩阵 Mo"V ”必jg妇心打0A.根据矩阵论A= AM =加蚣心,為)C'-CMD'=D线性变换得：t M 1 z Tx x Mzi+AfZ + BfU可以使用matlab进行对角标准型的

6、运算，matlab作为一种数学运算工具，很大程度的方便了了我们的计算，对于这个弹簧-质量-阻尼系统是一个四阶的状态空间表达式，所以可以用matlab简化计算。(1 )求特征值与特征向量A=0 0 1 0;0 0 0 1;-400 300 9 6;150 -200 3 -4.5B=0 0;0 0;1 0;0 0.5C=1 0 0 0;0 1 0 0P,J=eig(A)求得结果：0.0007-0.0402i0.0007 + 0.0402i0.0401-0.0698i0.0401 + 0.0698i-0.0171 + 0.0157i-0.0171 - 0.0157i0.0176 -0.0792i0.

7、0176 + 0.0792i0.86500.86500.6682 + 0.2084i0.6682 - 0.2084i-0.3442 -0.3621i-0.3442 + 0.3621i0.70500.70500.3667 +21.5183i00000.3667 -21.5183i00001.8833 + 8.4864i00001.8833 - 8.4864i(2)P矩阵求逆PN=i nv(P)求得结果:PN =3.4167 + 9.7803i-2.1017 - 9.2399i0.3466 - 0.2323i-0.4703 - 0.1054i3.4167 - 9.7803i-2.1017 + 9.

8、2399i0.3466 + 0.2323i-0.4703 + 0.1054i-3.3554 + 3.4224i3.7199 + 3.2032i0.2886 - 0.0353i0.5337 - 0.2409i-3.3554 - 3.4224i3.7199 - 3.2032i0.2886 + 0.0353i0.5337 + 0.2409i(3 )带入公式B PNB C CP0.0007- 0.0402i-0.0171 + 0.0157i0.0007+0.0402i-0.0171-0.0157i0.0401-0.0698i0.0176-0.0792i0.0401 +0.0698iu0.0176+0.

9、0792i解得对角标准型为:0.3667 21.5183i0000.3466-0.2323i-0.2352- 0.0527i- 00.3667 21.5183i00 -0.3466+ 0.2323i-0.2352 +0.0527ixxu001.8833 8.4864i00.2886-0.0353i0.2669-0.1205i0001.8833-8.4864i0.2886+ 0.0353i0.2669+ 0.1205i求状态空间表达式的解(1 )求状态转移矩阵T 1ATAte T其中，T为特征向量0.0007 - 0.0402i0.00070.0402i0.0401 -0.0698i0.0401

10、0.0698i0.366721.5183iteAt -0.01710.0157ie0.8650-0.3442 -0.3621i-0.0171 - 0.0157i0.0176 -0.0792i0.01760.0792i00.3667 -21.5183ite0.86500.66820.2084i0.6682 - 0.2084i-0.34420.3621i0.70500.705001.88338.4864ite001.8833 -8.4864ite*T 13.4167 + 9.7803i-2.1017 -9.2399i0.3466 -0.2323i-0.4703-0.1054i3.4167 -9.7

11、803i-3.3554 + 3.4224i-3.3554 -3.4224i-2.1017+ 9.2399i3.7199+ 3.2032i3.7199-3.2032i0.3466 + 0.2323i0.2886 -0.0353i0.2886 + 0.0353i-0.4703+ 0.1054i0.5337 - 0.2409i0.5337 + 0.2409i状态转移矩阵为：-5.5977-0.6477 + 0.0000i0.19990.4493At-4.4097 -0.0000i-0.5817 + 0.0000i0.2247 + 0.0000i0.5509 + 0.0000ie-12.5772 +

12、0.0000i-29.8835 -0.0000i-2.4502-1.4700-0.0000i-7.2316 + 0.0000i-42.7799-0.0000i-0.7350-0.0000i-1.7127 -0.0000i5可控性与可观性不同于经典控制理论， 能控性和能观性，是一个具有实际意义的概念，经典控制理论中用传递函数描述系统的输入-输出特性，输出量即被控量，只要系统是因果系统并且稳定，输出量便可以受控，且输出量总是可以被测量的， 因而不需要能控能观性的提出。 但是现代控制 理论是建立在状态空间表达式描述系统的基础上的，状态方程描述输入u (t)引起状态x(t)的变化过程，输出方程描述有状

13、态变化引起的输出y (t)的变化。能控能观便是定性的描述输入u (t)对状态x(t)的控制能力，输出y (t)对状态x( t)的反应能力，他们分 别回答了“输入能否控制状态的变化”-可控性“状态的变化能否有输出反映出来” 可观性x(t) eAtx(0 )0eA(t )Bu( )d另外在工程上常用状态变量作为反馈信息，可是状态x( t)的值通常是难测的，往往需要从测量到的y (t)中估计出状态，如果输出y (t)不能完全反映出系统的状态x( t)，那么就无法实现对状态的估计。能控性定义：当系统用状态方程描述时，给定系统的任意初始状态，可以找到允许的输入量，在有限的时间内使系统的所有状态达到任一终

14、止状态，则称系统是完全可控的。有状态方程 x' (t)=Ax(t)+Bu(t)其解为：tl Atl n 1in 1 i tlx(0 ) e A Bu( )dc( )AiBu( )dAiB ci ( )u( )d如果有限的0时间内0 < t < tl 0内通过输入量u(t)的作用把系统的所有状态引向状态 x(t1)设 人Xn 1 i令 i(tj= 0 Ci( )u( )dx(0 ) AiB i(t1)x(t1)=0，则应有：i 00(t)写成矩阵形式x(0 )= B AB An 1B 1n-1(t)即在给定x(0-)和A、B的条件下求可以使 x(t)=x(t1)的u(t)。换

15、言之：上述方程有解则系统 能控。根据凯莱哈米尔顿定理，e-At、 eAt可写成有限级数:x(t1) eAt1x(0 ):eA(t1)Bu( )d 0eAtx(0 )01eA(t Bu( )d 0x(0 )t1Ae0Bu( )d如果方程有解，等式右边左侧矩阵应满秩=n，此时系统是可控的。求可控性：001093Q B AB A* AB0000.532.251093301163.500.532.25163.580.875n=4满秩所以系统是可控的可观性定义：当系统用状态方程描述时，给定控制后，如果系统的每一个初始状态x(O-)x(t) eAtx(0 ) eA(t )Bu( )d都可以在有限的时间内通

16、过系统的输出y(t)唯一确定，则称系统完全可观。若只能确定部分初始状态，则称系统部分可观。(有状态方程X(° 0 (CAxG+BuBu( y(d=Cx(t)其解为由于在讨论能观性问题时，输入是给定的，上式右侧第二项是确知的，设u(t)=0y(t)=CeAtx(0-)。根据凯莱哈米尔顿定理，e-At、eAt可写成有限级数:如果方程有解，等式右侧中间侧矩阵应满秩。1000CAt 0eC0 (t)I1c:1(t)AQ2(t)A2n 1n 1Cn1(t)ACi(t)A0010i 0QoCA0001求可观性:其中，秩=n(系统的阶数)CA* An 1y(t) Cc400l c3(t?A E(t

17、6A2Cni(t)An1x(O )Ci(t)CA x(0 )150200 3 6i 0n=4满秩所以系统是可观的C写成矩系统的输入输出传递函数Cm(t) CA X(0 )对于两输入两输出的系统求得的传递函数是一个二阶的传递函数阵1，其中包含四个传递函数Ygig 丫i（s"U 2（s）Y2（sUi（s） Y2（s）/U2（s）Tran sfer function from in put 1 to output.sA2 + 4.5 s + 200#1:sA4 - 4.5 sA3 + 541.5 sA2 - 1800 s + 3.5e0043 s + 150#2:sA4 - 4.5 sA3

18、 + 541.5 sA2 - 1800 s + 3.5e004Tran sfer function from in put 2 to output.3 s + 150#1:sA4 - 4.5 sA3 + 541.5 sA2 - 1800 s + 3.5e0040.5 sA2 - 4.5 s + 200#2:sA4 - 4.5 sA3 + 541.5 sA2 - 1800 s + 3.5e004矩阵函数阵:3s 150sA2 4.5 s 200sA4-4.5sA3 541.5sA2-1800s 3.5e0043s 150sA4_4.5sA3 541.5 sA2-1800s 3.5e0040.5s

19、A2-4.5s 20sA4-4.5sA3 541.5sA2-1800s 3.5e004sA4-4.5sA3541.5sA2-1800s 3.5e0047分析开环稳定性稳定性定义是系统在受到小的外界扰动后，系统状态方程解的收敛性。系统正常工作要求是系统在受到外界扰动后，虽然其原平衡状态被打破，但在扰动消失后，仍然能恢复到原来的平衡状态，或者趋于另一平衡状态继续工作，且线性系统稳定性与输入作用无关。研究系统的稳定性对于研究系统能否正常工作具有很重要的意义，稳定性是自动控制系统正常工作的必要条件，是系统的重要特征。我们不仅要分析一个系统是否稳定还要解决的问题便是怎样使一个系统稳定。经典控制理论稳定性

20、判别方法有很多，例如代数判据，niquist判据，根轨迹判据等。而现代控制理论经常用李雅普诺夫第二法求稳定性。(1 )利用特征根的方法根据上述结果求得的特征根为0.3667 +21.5183i0.3667 +21.5183i1.8833 + 8.4864i1.8833 - 8.4864i，四个特征值全部都在坐标轴的右半平面，所以系统是不稳定的。(2)利用利亚普诺夫第二法求解2 2 2 2v(x) °.5(X1X2X3X4)V(X)X1X1 X2X2 X3X3 X4X42 2300x2X3 9X3 9X3X4 150 XiX4 4.5X4v(x)其中:0 0399o (A,B,C)0

21、0300将其换成矩阵形式 A0 090 0 0150199可以看出04.5k (A BK,B,C)A不是正定的，所以系统不稳定。X1X3X2X4X3400 X1300 X29X36X4X4叫1叫3X34.5X4199 X2 X4 399 X1 X38利用状态反馈将系统的闭环极点配置到合适的值状态反馈是将系统每一个状态变量乘以相应的反馈系数馈送到输入端与参考输入相加，其和作为受控系统的控制输入。原受控对象为经过状态反馈后得到的闭环系统为闭环系统期望极点选取原则为以下几点：1）n维控制系统有n个期望极点；2）期望极点是物理上可实现的，为实数或共轭复数对；3） 期望极点的位置的选取，需考虑它们对系统

22、品质的影响（离虚轴的位置），及与零点分 布状况的关系。4）离虚轴距离较近的主导极点收敛慢，对系统性能影响最大，远极点收敛快，对系统只有 极小的影响。闭环极点 0.3667 +21.5183i0.3667 -21.5183i1.8833 + 8.4864i1.8833 - 8.4864i配置状态反馈后，系统应稳定，所以期望极点应在虚轴左侧，所以期望闭环极点-1.1001+64.5549i-1.1001-64.5549i-5.6499+25.4952-5.6499-25.4952得到极点配置矩阵K =1.0e+003 *-1.1001 -64.5549i-1.1001 +64.5549i-5.64

23、99 -25.4952i-5.64991.2556-0.03750.0157-0.03320.97182.89690.08390.004验证极点配置结果是正确的：ans =+25.4952i求得开环传递函数阶跃响应曲线(没有经过状态反馈的)Stsp Responses.-.x 104 E： From： ln(1Emo BLCM¥°bRtseTime. N；ATime (b&c)没有上升时间经过状态反馈的传递函数：Gk(s) Csl (A BK) 1B状态空间表达式为k(A BK),B,CMatlab解得闭环传递函数：sA2 + 6.758 s + 1648#1:sA

24、4 + 13.5 sA3 + 4875 sA2 + 4.86e004 s + 2.843e006-38.95 s - 335.9#2:sA4 + 13.5 sA3 + 4875 sA2 + 4.86e004 s + 2.843e006Tran sfer function from in put 2 to output.19.59 s + 168.7#1:sA4 + 13.5 sA3 + 4875 sA2 + 4.86e004 s + 2.843e0060.5 sA2 + 3.371 s + 827.8#2:sA4 + 13.5 sA3 + 4875 sA2 + 4.86e004 s + 2.8

25、43e006反馈后的阶跃响应:阶跃响应上升时间是 0.034s，配置后系统最终稳定。9设计全维状态观测器不是所有的系统状态物理上都能够直接测量得到。需要从系统的可量测参量，如输入u和输出y来估计系统状态 。状态观测器基于可直接量测的输出变量 y和控制变量u来估计状 态变量，是一个物理可实现的模拟动力学系统。 如果系统状态是完全能观测的， 那么根据输 出y的测量，可以唯一的确定系统的初始状态， 所以只要满足一定的条件， 便可以从可测量y和u中把x间接重构出来。全维渐进状态结构图：1.8833系统的原极点为 0.3667 +21.5183i0.3667 +21.5183i1.8833 + 8.48

26、64i-8.4864i期望极点应是原极点的2-5倍，-33.9456iA并不是越快2igA KeC期望极点为-1.4668 +86.0732iTC-7.+33.9456i-7.533217.551.152.65.03199.91608.8122.32246.1相应的全维观测器是o=aIz0000Iz009t7乙乙 &乙即00Iz0609 Iz 00 充0a000Iz986t7 3939-00Iz0£1/旳 G"000Iz9 t7-£9m- 992乙V :(H9660£k009£-10986 t7 -S9Z901£lz lz9-

27、lz9Zk|/9t7乙乙&乙即90 09t7£01 9如乙LIZA88091，666比0 L968'80£l，666工nX0 99390 0Iz099391/旳9Zlz0 00Iz9lz99Zlz|/9t7乙乙&乙即90 000l/9t7乙乙&勿9t7£003091，A88091，666比0 L008'8091，666氏9600£00t7inX099390 000 099Z9Iz0001/旳9Zlz0 000 1/旳9Zlz0Iz00|/9t7乙乙&勿9 00|/9t7乙乙&勿9t7£00

28、3091，r 8'8091，666氏0Iz00 L 08'8091，666氏9600£00t7AnX099390000 0 L099Z9Iz0001/旳9Zlz001/旳9Zlz0Iz00Al ng x(g-| v) X分别得到输出和观测状态的传递函数Tran sfer function from in put "y1" to output.17.51 sA3 + 5899 sA2 + 1.23e005 s + 7.911e006yi_e：sA4 + 18 sA3 + 8664 sA2 + 1.152e005 s + 8.96e006-52.62

29、sA3 + 359.1 sA2 - 6844 s + 4.369e005y2_e:sA4 + 18 sA3 + 8664 sA2 + 1.152e005 s + 8.96e00617.51 sA3 + 5899 sA2 + 1.23e005 s + 7.911e006x1_e:sA4 + 18 sA3 + 8664 sA2 + 1.152e005 s + 8.96e006-52.62 sA3 + 359.1 sA2 - 6844 s + 4.369e005x2_e:sA4 + 18 sA3 + 8664 sA2 + 1.152e005 s + 8.96e0063200 sA3 + 9.295e

30、004 sA2 + 7.698e006 s + 3.967e006x3_e:sA4 + 18 sA3 + 8664 sA2 + 1.152e005 s + 8.96e006122.3 sA3 + 1.405e005 sA2 - 8571 s + 5.737e007x4_e:sA4 + 18 sA3 + 8664 sA2 + 1.152e005 s + 8.96e006Tran sfer function from in put "y2" to output.51.13 sA3 + 1379 sA2 + 2.945e004 s + 8.583e005y1_e：sA4 + 18

31、 sA3 + 8664 sA2 + 1.152e005 s + 8.96e0064.985 sA3 + 5001 sA2 + 1.152e005 s + 8.303e006y2_e:sA4 + 18 sA3 + 8664 sA2 + 1.152e005 s + 8.96e006x1_e:51.13 sA3 + 1379 sA2 + 2.945e004 s + 8.583e005x2_e:sA4 + 18 sA3 + 8664 sA2 + 1.152e005 s + 8.96e0064.985 sA3 + 5001 sA2 + 1.152e005 s + 8.303e006x3_e:sA4 + 18 sA3 + 8664 sA2 + 1.152e005 s + 8.96e0061609 sA3 - 1.337e005 sA2 + 1