第二章：动力学系统的微分方程模型.docx

《第二章：动力学系统的微分方程模型.docx》由会员分享，可在线阅读，更多相关《第二章：动力学系统的微分方程模型.docx（27页珍藏版）》请在第一文库网上搜索。

1、第二章：动力学系统的微分方程模型利用计算机进行仿真时，一般情况下要给出系统的数学模型，因此有必要掌握一定的建立数学模型的方法。在动力学系统中，大多数情况下可以使用微分方程来表示系统的动态特性，也可以通过微分方程可以将原来的系统简化为状态方程或者差分方程模型等。在这一章中，重点介绍建系统动态问题的微分方程的基本理论和方法。在实际工程中，一般把系统分为两种类型，一是连续系统；其数学模型一般是高阶微分方程；另一种是离散系统，它的数学模型是差分方程。2.1 动力学系统统基本元件任何机械系统都是由机械元件组成的，在机械系统中有3种类型的基本机械元件：惯性元件、弹性元件和阻尼元件。1惯性元件：惯性元件是指

2、具有质量或转动惯量的元件，惯量可以定义为使加速度（或角加速度）产生单位变化所需要的力（或力矩）。2弹性元件：它在外力或外力偶作用下可以产生变形的元件，这种元件可以通过外力做功来储存能量。按变形性质可以分为线性元件和非线性元件，通常等效成一弹簧来表示。对于线性弹簧元件，弹簧中所受到的力与位移成正比，比例常数为弹簧刚度Z。F=Ax这里女称为弹簧刚度，Ar是弹簧相对于原长的变形量，弹性力的方向总是指向弹簧的原长位移，出了弹簧和受力之间是线性关系以外，还有所谓硬弹簧和软弹簧，它们的受力和弹簧变形之间的关系是一非线性关系。3阻尼元件：这种元件是以吸收能量以其它形式消耗能量，而不储存能量，可以形象的表示为

3、一个活塞在一个充满流体介质的油缸中运动。阻尼力通常表示为：R=cxa阻尼力的方向总是速度方向相反。当=1,为线性阻尼模型。否则为非线性阻尼模型。应注意当a等于偶数情况时，要将阻尼力表示为：R=-dr,这里的“-”表示与速度方向相反2.2 动力学建模基本定理1动力学普遍定理对于大多数力学问题，可以使用我们熟知的牛顿动力学基本定理来解决，动力学普遍定理包括动量定理、动量矩定理和动能定理，以及其他变形形式，普遍定理的特点是比较直观，针对不同的问题可以选择不同的力学定理，在一般情况下利用普遍定理可以得到大多数动力学系统的数学模型。1)动量定理与质心运动定理：设系统在任意瞬时的动量矢为作用在系统上的外力

4、矢量和为ZE,则任意瞬时的动量对时间的导数等于作用在系统中所有外力的矢量和构成了动量定理。-=Yf(2-1)dt乙通常将该式投影到直接坐标轴系、自然坐标轴系等，(更详细的情况请参阅理论力学有关知识)利用质心坐标的计算表达式，可以将动量定理转化为质心运动定理，即：MGZFi或：Z哂=EFi(2-2)其中：M是系统的总质量，应是系统的质心；?是分刚体是质心，心,是分刚体的质心。2)动量矩定理：系统在任意瞬时的动量矩对时间的导数等于作用在系统中所有外力矩的矢量和。哼*M(F)(2-3)其中，X。是系统对固定点。的动量矩，MO(F)力F对0点的矩.除了对固定点的动量矩定理外，还有对质心的动量矩定理，对

5、速度瞬心的动量矩定理和对加速度瞬心的动量矩定理。3)动能定理：动能定理的导数形式：系统在任意瞬时的动能对时间的导数等于作用在系统中所有力的功率的代数和。W=Zn(2-4)动能定理的积分形式：系统在任意两瞬时的动能的变化等于作用在系统中所有力的功的代数和。与一(=ZW2动力学普遍方程将达朗伯原理与虚位移原理相结合，得到了建立动力学模型的另一种方法。1）达朗伯原理达朗伯原理提供了研究动力学问题的一个新的方法，即借助于惯性力（Q=-ma）的概念，可用研究静力学平衡的方法来研究动力学问题，这种方法常称为动静法。即：在任意时刻，质点在主动力、约束力和惯性力的主矢作用下处于平衡；E+M+2=o（2-5）以

6、及主动力、约束力和惯性力对某点的矩矢等于零，即：府。（E）+河。）+Mo（Qi）=O通常先计算惯性力的主矢和主矩，从而得到质点系的达朗伯原理。2）虚位移原理虚位移原理本身是通过虚功的引入，提出了求解静力学问题的一种方法,它与达朗伯原理相结合得到了建立动力学模型的另一种方法。对于理想约束的完整系统，质点（质点系）在其给定位置上处于平衡的必要充分条件是作用在该质点（质点系）上的所有主动力E在其作用点的虚位移近上所做的虚功和等于零，即：M=O或（乙,*+&母+K回）=3）动力学的普遍方程受理想约束的系统，作用在质点系上的所以主动力和惯性力在各自的虚位移上所做的虚功和等于零，即：（ETnmi）师=U/

7、=1或（心-nX）西+（即-明%）双十（4-gZ）的=0/=1在具体应用这个方程的时候，可以先引入广义坐标，使得问题处理简单。解：杆绕。轴做定轴转动，水平位置为系统例2T质量为加均质的杆可以绕。轴定动,试求系统做微幅振动时的微分方程。的平衡状态，取杆绕O轴转动的角度。为坐标，可以方便的使用动量矩定理来建立动力学方程。(假定在微小转动情况下)J=f(t)a-cia+k3a)3a这里/是杆绕O轴转动的转动惯量。这是关于的二阶线性微分方程。如果不计杆的质量，则微分方程为：9ca+9ka=/.X这个方程是关于的一阶线性微分方程，I二fn称该系统模型为一阶系统。河”FR广W例2-2怂浮摆的动力学建模下图

8、所示为-Mr0小型起重机简图，机,相2是吊车和吊重的质一/量，吊绳长为/且不计质量，吊车的驱动力为F,考虑轨道的阻力为a,试以Xe为广义坐A-,叫标，建立系统的动力学控制方程。利用水平方向的质心运动定理，即：m1x+m2-(x+sin6)=F-cx(1)dt或：/n1x+n2(x+cos-sin2)=F-cx重物做平面曲线运/J,则可以直接利用牛顿定律得到切线方向的动力学方程：m2(J+xcqs)=-m2gsin(2)(1),(2)两式是耦合的非线性动力学方程。当系统被限制在。=0附近运动时，可将其在。=0处线性化处理，则可以得到系统的方程为：(w1+m2)x+cx-m210)=Fwix+n2

9、(x+1-12)=F当给定尸二Fa)时，可以建立仿真模型。请读者考虑，如果要考虑摆杆的质量，则动力学方程如何？例2-3：车辆悬架系统的动力学模型考虑图2.2所示的汽车悬架系统示意图。设计悬架缓冲系统的匕2；右,。2的目的是减小车辆在崎岖道路上行驶时产生的震动，因为道路表面的不平坦会引起悬架沿垂直方向的移动和绕某个轴的转动。图2.2悬架系统示意图图2.3架系统的受力分析示意图我们将整个系统的质量中心作为坐标的原点,因此系统在不平道路上的振动运动可以看作是质心的沿垂直方向的平移运动以及绕质心的旋转运动。车架质量为m,转动惯量为J。输入车轮的位置信息%、为表明路况信息。假设每个车轴的缓冲系统由具有阻

10、尼特性的弹簧构成。忽略轮胎的质量,每个车轮受到的外力为弹簧弹力与阻尼力之和,即fA=(G1+K)%=(CJA+幻A)FB=(。2+k2)y(S)=(C2yB+k2y)at其中：。一切yB=y-h(p-yiy和力分别表示每个弹簧距离参考位置的瞬时距离。代入上式后=(Ci4K)(y+。-%)atFB=(c2-+k2)(y-b-y2)atdy根据质心运动与相对于质心的动量矩定理得：M-=-Fa-Fr或者：my=-c1(y+a-yi)-ki(y+a-yj-c2(y-b-y2)-k2(y-b-y2)整理后得到：my+(c1+c2)y+(占+k2)y+(cc2b)-(kxa-k2b)=c1y1+c2y2+

11、kiy1+k2y2用y)和。分别表示系统质心的平移位移和沿质心的旋转角度。上式中假定在很小的角度位置条件下满足Sinoko,并且e取顺时针的旋转方向为正方向。再根据系统相对于质心的动量矩定理可得：J-j=FBbCGS(P-FAaCGS(PFbb-Faa其中J是车驾相对于质心的转动惯量，将上式整理后可得：J-r=(c2-+k2)(y-b-y2)b-(ci-+k1)(y+a-y1)ad厂atat或：J+(c2b+CIa)0+(C1-c2)y+k2b2+ka2)+(ka-k2b)y=Cqj1-c2by2+k1ayy-k2by2将系统的动力学方程写成矩阵形式：昂Bi2G1C2昂巴2B2I%1G1简写为

12、：A）：+By+C,其中：A=C1Ga,OJ网=+A,B；+A-,C尸2122J1Jy,=必+用1y2C=仁+&招一k?bkia-k2bak2b小=A-iE5+A-,F1)2当川为非奇异阵时，可以通过矢量信号我们可以得到系统的仿真模型如（图2-5)。图2.5悬架系统仿真框图以上系统中假定M、y2是系统两个相互独立的输入变量,但实际上,后轮与前车轮的位置时间相差A时间。这样,实际系统满足%=y（f-由于借助了拉斯变换，将微分方程换成了代数方程，如果要得到时域响应则需要借助拉斯反变换。根据第章的基本知识，给出基于微分方程的仿真模型，具体计算过程留给读者练习。例2-4机构运动学建模曲柄滑块机构的运动

13、学仿真建模（速度分析与建模）曲柄滑块机构如图所示：该机构只有一个自由度，首先给出机构的运动学分析模型，（1）机构的封闭的矢量方程尸=A弓（2）矢量方程的分解式COSeI+r2cos=r外sine+Sin*2=0(3)关于机构速度问题的运动学方程；-4sinxx-RSin11=rCOSe01+弓COS0202=0机构的输入运动量为,(j,输出量为2,2,r,r,写成矩阵形式12一qsinq=jJ-CoSe1可以写成显式表达式02r2SIr1%r2COSe2一OSin/-cosr2Sine2r2CoS外Simu1ink仿真模型建立在该仿真模型中，设系统的输入角速度为：=150弧度/秒，通过一次积分

14、可以得到角度怯，将这两个输入量通过一个信号混合器（以向量形式混合为一路信号），输入给MAT1ABFCN模块，通过该函数模块中的代码入，从而可以得到输出量（物，户），再进一步积分后,得到位移量物”。在MAT1ABFUNTION模块中写上函数过程文件名：Compv,其它不变,建立m脚本文件如下：(函数子程序)function(x=compv(u);x输出，(U)输入。%参数说明：r1曲柄长度，r2连杆长度%U(I)曲柄角速度；u(2)曲柄角度，u(3)连杆角度r1=15;r2=55;a=r2*sin(u(3)1;r2*cos(u(3)0;b=-u(1)*r1*sin(u(2);cos(u(2);x=inv(a)*b;将该文件名储存为compv.m,然后运行仿真模型，得如下结果。图2-10图2-11连杆的角速度与角度的变化规律滑快的速度与位移变化规律曲柄滑块机构的运动学仿真(加速度分析)加速度表达式-rx(cos(p、:+sini)-r2(cos21+sin22)=r4(一Sine0；+cos910J+g(-Sin/帮+cos22)=0机构的输入运动量为r2S1ne2r2cos2，输出量为2,2,2yr,r,f,写成矩阵形式02r2SIne2