弹簧-质量-阻尼模型

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

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

3、=6Ns/cm。由图2.1，根据牛顿第二定律,，建立系统的动力学模型如下：对有： （2-1）对有： （2-2）3 建立状态空间表达式令,则原式可化为：化简得： （2-3） （2-4）整理得： （2-5）代入数据得： 则系统的状态空间表达式为4 化为对角标准型当系统矩阵A有n个不相等的特征根时，相应的有n个不相等的特征向量，所以有矩阵A的特征矩阵根据矩阵论线性变换得： 可以使用matlab进行对角标准型的运算，matlab作为一种数学运算工具，很大程度的方便了了我们的计算，对于这个弹簧-质量-阻尼系统是一个四阶的状态空间表达式，所以可以用matlab简化计算。（1） 求特征值与特征向量A=0 0

4、 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)求得结果：P = 0.0007 - 0.0402i 0.0007 + 0.0402i 0.0401 - 0.0698i 0.0401 + 0.0698i -0.0171 + 0.0157i -0.0171 - 0.0157i 0.0176 - 0.0792i 0.0176 + 0.0792i 0.8650 0.8650 0.6682 + 0.2084i 0.6682 - 0.2084i -0.3442 - 0.3621i

5、-0.3442 + 0.3621i 0.7050 0.7050 J = 0.3667 +21.5183i 0 0 0 0 0.3667 -21.5183i 0 0 0 0 1.8833 + 8.4864i 0 0 0 0 1.8833 - 8.4864i（2） P矩阵求逆PN=inv(P)求得结果：PN = 3.4167 + 9.7803i -2.1017 - 9.2399i 0.3466 - 0.2323i -0.4703 - 0.1054i 3.4167 - 9.7803i -2.1017 + 9.2399i 0.3466 + 0.2323i -0.4703 + 0.1054i -3.35

6、54 + 3.4224i 3.7199 + 3.2032i 0.2886 - 0.0353i 0.5337 - 0.2409i -3.3554 - 3.4224i 3.7199 - 3.2032i 0.2886 + 0.0353i 0.5337 + 0.2409i（3） 带入公式 解得对角标准型为：5求状态空间表达式的解（1） 求状态转移矩阵其中，T为特征向量状态转移矩阵为：5 可控性与可观性不同于经典控制理论，能控性和能观性，是一个具有实际意义的概念，经典控制理论中用传递函数描述系统的输入-输出特性，输出量即被控量，只要系统是因果系统并且稳定，输出量便可以受控，且输出量总是可以被测量的，因而

7、不需要能控能观性的提出。但是现代控制理论是建立在状态空间表达式描述系统的基础上的，状态方程描述输入u（t）引起状态x（t）的变化过程，输出方程描述有状态变化引起的输出y（t）的变化。能控能观便是定性的描述输入u（t）对状态x(t)的控制能力，输出y（t）对状态x（t）的反应能力，他们分别回答了“输入能否控制状态的变化”-可控性“状态的变化能否有输出反映出来”-可观性另外在工程上常用状态变量作为反馈信息，可是状态x（t）的值通常是难测的，往往需要从测量到的y（t）中估计出状态，如果输出y（t）不能完全反映出系统的状态x（t），那么就无法实现对状态的估计。能控性定义：当系统用状态方程描述时，给定系

8、统的任意初始状态，可以找到允许的输入量，在有限的时间内使系统的所有状态达到任一终止状态，则称系统是完全可控的。有状态方程x(t)=Ax(t)+Bu(t) 其解为：如果有限的时间内0 t t1内通过输入量u(t)的作用把系统的所有状态引向状态x(t1)设x(t1)=0 ，则应有：即在给定x(0-)和A、B的条件下求可以使x(t)=x(t1)的u(t)。换言之：上述方程有解则系统能控。根据凯莱哈米尔顿定理， e-At、 eAt可写成有限级数：如果方程有解，等式右边左侧矩阵应满秩=n，此时系统是可控的。求可控性：n=4 满秩所以系统是可控的可观性定义：当系统用状态方程描述时，给定控制后，如果系统的每

9、一个初始状态x(0-)都可以在有限的时间内通过系统的输出y(t)唯一确定，则称系统完全可观。若只能确定部分初始状态，则称系统部分可观。有状态方程x(t)=Ax(t)+Bu(t) y(t)=Cx(t) 其解为由于在讨论能观性问题时，输入是给定的，上式右侧第二项是确知的，设u(t)=0。y(t)=CeAtx(0-)。根据凯莱哈米尔顿定理， e-At 、eAt可写成有限级数：如果方程有解，等式右侧中间侧矩阵应满秩。其中，秩=n(系统的阶数)求可观性：n=4满秩所以系统是可观的6 求系统的输入输出传递函数对于两输入两输出的系统求得的传递函数是一个二阶的传递函数阵,其中包含四个传递函数Transfer

10、function from input 1 to output. s2 + 4.5 s + 200 #1: - s4 - 4.5 s3 + 541.5 s2 - 1800 s + 3.5e004 3 s + 150 #2: - s4 - 4.5 s3 + 541.5 s2 - 1800 s + 3.5e004 Transfer function from input 2 to output. 3 s + 150 #1: - s4 - 4.5 s3 + 541.5 s2 - 1800 s + 3.5e004 0.5 s2 - 4.5 s + 200 #2: - s4 - 4.5 s3 + 541

11、.5 s2 - 1800 s + 3.5e004矩阵函数阵：7 分析开环稳定性稳定性定义是系统在受到小的外界扰动后，系统状态方程解的收敛性。系统正常工作要求是系统在受到外界扰动后，虽然其原平衡状态被打破，但在扰动消失后，仍然能恢复到原来的平衡状态，或者趋于另一平衡状态继续工作，且线性系统稳定性与输入作用无关。研究系统的稳定性对于研究系统能否正常工作具有很重要的意义，稳定性是自动控制系统正常工作的必要条件，是系统的重要特征。我们不仅要分析一个系统是否稳定 还要解决的问题便是怎样使一个系统稳定。经典控制理论稳定性判别方法有很多，例如代数判据，niquist判据，根轨迹判据等。而现代控制理论经常用李

12、雅普诺夫第二法求稳定性。（1）利用特征根的方法根据上述结果求得的特征根为 0.3667 +21.5183i 0.3667 +21.5183i 1.8833 + 8.4864i 1.8833 - 8.4864i ，四个特征值全部都在坐标轴的右半平面，所以系统是不稳定的。（2）利用利亚普诺夫第二法求解其中：将其换成矩阵形式可以看出A不是正定的，所以系统不稳定。8 利用状态反馈将系统的闭环极点配置到合适的值状态反馈是将系统每一个状态变量乘以相应的反馈系数馈送到输入端与参考输入相加，其和作为受控系统的控制输入。原受控对象为经过状态反馈后得到的闭环系统为闭环系统期望极点选取原则为以下几点：1）n维控制系

13、统有n个期望极点；2）期望极点是物理上可实现的，为实数或共轭复数对；3）期望极点的位置的选取，需考虑它们对系统品质的影响（离虚轴的位置），及与零点分布状况的关系。4）离虚轴距离较近的主导极点收敛慢，对系统性能影响最大，远极点收敛快，对系统只有极小的影响。闭环极点 0.3667 +21.5183i 0.3667 -21.5183i 1.8833 + 8.4864i 1.8833 - 8.4864i配置状态反馈后，系统应稳定，所以期望极点应在虚轴左侧，所以期望闭环极点 -1.1001+64.5549i -1.1001-64.5549i -5.6499+25.4952 -5.6499-25.4952

14、得到极点配置矩阵K = 1.0e+003 * 1.2556 -0.0375 0.0157 -0.0332 0.9718 2.8969 0.0839 0.004验证极点配置结果是正确的：ans = -1.1001 -64.5549i -1.1001 +64.5549i -5.6499 -25.4952i -5.6499 +25.4952i求得开环传递函数阶跃响应曲线（没有经过状态反馈的）：没有上升时间经过状态反馈的传递函数：状态空间表达式为Matlab解得闭环传递函数： s2 + 6.758 s + 1648 #1: - s4 + 13.5 s3 + 4875 s2 + 4.86e004 s +

15、 2.843e006 -38.95 s - 335.9 #2: - s4 + 13.5 s3 + 4875 s2 + 4.86e004 s + 2.843e006 Transfer function from input 2 to output. 19.59 s + 168.7 #1: - s4 + 13.5 s3 + 4875 s2 + 4.86e004 s + 2.843e006 0.5 s2 + 3.371 s + 827.8 #2: - s4 + 13.5 s3 + 4875 s2 + 4.86e004 s + 2.843e006反馈后的阶跃响应：阶跃响应上升时间是0.034s，配置后

16、系统最终稳定。9 设计全维状态观测器不是所有的系统状态物理上都能够直接测量得到。需要从系统的可量测参量，如输入u和输出y来估计系统状态 。状态观测器基于可直接量测的输出变量y和控制变量u来估计状态变量，是一个物理可实现的模拟动力学系统。如果系统状态是完全能观测的，那么根据输出y的测量，可以唯一的确定系统的初始状态，所以只要满足一定的条件，便可以从可测量y和u中把x间接重构出来。全维渐进状态结构图：系统的原极点为0.3667 +21.5183i 0.3667 +21.5183i 1.8833 + 8.4864i 1.8833 - 8.4864i 期望极点应是原极点的2-5倍，并不是越快越好。期望

17、极点为 -1.4668 +86.0732i -1.4668 -86.0732i -7.5332 +33.9456i -7.5332 -33.9456i相应的全维观测器是： 10 带观测器的输出状态空间表达式： D=0 分别得到输出和观测状态的传递函数：Transfer function from input y1 to output. 17.51 s3 + 5899 s2 + 1.23e005 s + 7.911e006 y1_e: - s4 + 18 s3 + 8664 s2 + 1.152e005 s + 8.96e006 -52.62 s3 + 359.1 s2 - 6844 s + 4

18、.369e005 y2_e: - s4 + 18 s3 + 8664 s2 + 1.152e005 s + 8.96e006 17.51 s3 + 5899 s2 + 1.23e005 s + 7.911e006 x1_e: - s4 + 18 s3 + 8664 s2 + 1.152e005 s + 8.96e006 -52.62 s3 + 359.1 s2 - 6844 s + 4.369e005 x2_e: - s4 + 18 s3 + 8664 s2 + 1.152e005 s + 8.96e006 3200 s3 + 9.295e004 s2 + 7.698e006 s + 3.96

19、7e006 x3_e: - s4 + 18 s3 + 8664 s2 + 1.152e005 s + 8.96e006 122.3 s3 + 1.405e005 s2 - 8571 s + 5.737e007 x4_e: - s4 + 18 s3 + 8664 s2 + 1.152e005 s + 8.96e006 Transfer function from input y2 to output. 51.13 s3 + 1379 s2 + 2.945e004 s + 8.583e005 y1_e: - s4 + 18 s3 + 8664 s2 + 1.152e005 s + 8.96e006

20、 4.985 s3 + 5001 s2 + 1.152e005 s + 8.303e006 y2_e: - s4 + 18 s3 + 8664 s2 + 1.152e005 s + 8.96e006 51.13 s3 + 1379 s2 + 2.945e004 s + 8.583e005 x1_e: - s4 + 18 s3 + 8664 s2 + 1.152e005 s + 8.96e006 4.985 s3 + 5001 s2 + 1.152e005 s + 8.303e006 x2_e: - s4 + 18 s3 + 8664 s2 + 1.152e005 s + 8.96e006 16

21、09 s3 - 1.337e005 s2 + 1.373e006 s - 1.853e007 x3_e: - s4 + 18 s3 + 8664 s2 + 1.152e005 s + 8.96e006 2246 s3 + 2.437e004 s2 + 6.754e006 s - 4.844e007 x4_e: - s4 + 18 s3 + 8664 s2 + 1.152e005 s + 8.96e006 Input groups: Name Channels Measurement 1,2 Output groups: Name Channels OutputEstimate 1,2 StateEstimate 3,4,5,6 分析其稳定性根据阶跃响应图来看系统是稳定的，并且上升时间为0.0015s附录：