第二章动力学系统的微分方程模型.docx
- 文档编号:10296901
- 上传时间:2023-02-09
- 格式:DOCX
- 页数:41
- 大小:562.94KB
第二章动力学系统的微分方程模型.docx
《第二章动力学系统的微分方程模型.docx》由会员分享,可在线阅读,更多相关《第二章动力学系统的微分方程模型.docx(41页珍藏版)》请在冰豆网上搜索。
第二章动力学系统的微分方程模型
第二章:
动力学系统的微分方程模型
利用计算机进行仿真时,一般情况下要给出系统的数学模型,因此有必要掌握一定的建立数学模型的方法。
在动力学系统中,大多数情况下可以使用微分方程来表示系统的动态特性,也可以通过微分方程可以将原来的系统简化为状态方程或者差分方程模型等。
在这一章中,重点介绍建系统动态问题的微分方程的基本理论和方法。
在实际工程中,一般把系统分为两种类型,一是连续系统;其数学模型一般是高阶微分方程;另一种是离散系统,它的数学模型是差分方程。
§2.1动力学系统统基本元件
任何机械系统都是由机械元件组成的,在机械系统中有3种类型的基本机械元
件:
惯性元件、弹性元件和阻尼元件。
1惯性元件:
惯性元件是指具有质量或转动惯量的元件,惯量可以定义为使加速度
(或角加速度)产生单位变化所需要的力(或力矩)。
这种元件可以通过外
2弹性元件:
它在外力或外力偶作用下可以产生变形的元件,
力做功来储存能量。
按变形性质可以分为线性元件和非线性元件,通常等效成一弹簧来表示。
对于线性弹簧元件,弹簧中所受到的力与位移成正比,比例常数为弹簧刚度k。
FWx
这里k称为弹簧刚度,级是弹簧相对于原长的变形量,弹性力的方向总是指向弹
簧的原长位移,出了弹簧和受力之间是线性关系以外,还有所谓硬弹簧和软弹簧,它们的受力和弹簧变形之间的关系是一非线性关系。
3阻尼元件:
这种元件是以吸收能量以其它形式消耗能量,而不储存能量,可以形
象的表示为一个活塞在一个充满流体介质的油缸中运动。
阻尼力通常表示为:
Da
R=ex
阻尼力的方向总是速度方向相反。
当1,为线性阻尼模型。
否则为非线性阻
尼模型。
应注意当:
等于偶数情况时,要将阻尼力表示为:
R--ex|x4|这里的"-”表示与速度方向相反
§22动力学建模基本定理
1动力学普遍定理
对于大多数力学问题,可以使用我们熟知的牛顿动力学基本定理来解决,动力学普遍定理包括动量定理、动量矩定理和动能定理,以及其他变形形式,普遍定理的特点是比较直观,针对不同的问题可以选择不同的力学定理,在一般情况下利用普遍定理可以得到大多数动力学系统的数学模型。
1)动量定理与质心运动定理:
设系统在任意瞬时的动量矢为K,作用在系统上的外力矢量和为Fj,则
任意瞬时的动量对时间的导数等于作用在系统中所有外力的矢量和构成了动量定理。
坐八F(2-1)
dt
通常将该式投影到直接坐标轴系、自然坐标轴系等,(更详细的情况请参阅理论力学有关知识)
利用质心坐标的计算表达式,可以将动量定理转化为质心运动定理,即:
Mac="Fi或:
二.miaci="Fi(2-2)
ad是分
其中:
M是系统的总质量,ac是系统的质心;mi是分刚体是质心,
刚体的质心。
2)动量矩定理:
系统在任意瞬时的动量矩对时间的导数等于作用在系统中所有外力矩的矢量和。
其中,H0是系统对固定点O的动量矩,MO(F)力F对O点的矩.
除了对固定点的动量矩定理外,还有对质心的动量矩定理,对速度瞬心的动量矩定理和对加速度瞬心的动量矩定理。
3)动能定理:
动能定理的导数形式:
系统在任意瞬时的动能对时间的导数等于作用在系统中所有力的功率的
代数和。
匹八*N(2-4)
dt
动能定理的积分形式:
系统在任意两瞬时的动能的变化等于作用在系统中所
有力的功的代数和。
T2-Ti=為W
2动力学普遍方程
将达朗伯原理与虚位移原理相结合,得到了建立动力学模型的另一种方法。
1)达朗伯原理达朗伯原理提供了研究动力学问题的一个新的方法,即借助
于惯性力(Q-_ma)的概念,可用研究静力学平衡的方法来研究动力学问题,这种方法常称为动静法。
即:
在任意时刻,质点在主动力、约束力和惯性力的主矢作用下处于平衡;
ZFi+瓦Ni+£Qi=0(2-5)
以及主动力、约束力和惯性力对某点的矩矢等于零,即:
'M°(Fi)'M°(NJ'M°(QJ=0
通常先计算惯性力的主矢和主矩,从而得到质点系的达朗伯原理。
2)虚位移原理
虚位移原理本身是通过虚功的引入,提出了求解静力学问题的一种方法,
它与达朗伯原理相结合得到了建立动力学模型的另一种方法。
对于理想约束的完整系统,质点(质点系)在其给定位置上处于平衡的必要
充分条件是作用在该质点(质点系)上的所有主动力Fi在其作用点的虚位移T上所做的虚功和等于零,即:
或送(Fix輙+Fiy话yi+Fiz也)=0
3)动力学的普遍方程
受理想约束的系统,作用在质点系上的所以主动力和惯性力在各自的虚位移上所做的虚功和等于零,即:
n
迟(Fj-miai)^r=0
i=1
n
或'[(Fxi-mix)%(Fyj-miiyjyi亿-rnijZ)z]=0
iA
在具体应用这个方程的时候,可以先引入广义坐标,使得问题处理简单。
例2-1质量为m均质的杆可以绕O轴定动,试求系统做微幅振动时的微分方程。
解:
杆绕0轴做定轴转动,水平位置为系统
■a一3a
的平衡状态,取杆绕O轴转动的角度「为坐标,可以方便的使用动量矩定理来建
立动力学方程。
(假定在微小转动情况下)
J密=f(t)a—(c3ak3a)3a
这里J是杆绕O轴转动的转动惯量。
这是关于:
的二阶线性微分方程。
如果不计杆
的质量,则微分方程为:
9ca「9ka「二f(t)
这个方程是关于「的一阶线性微分方程,称该系统模型为一阶系统。
例2-2悬浮摆的动力学建模下图所示为
小型起重机简图,mi,m2是吊车和吊重的质量,吊绳长为I且不计质量,吊车的驱动力为F,考虑轨道的阻力为cx,试以XR为广义坐标,建立系统的动力学控制方程。
利用水平方向的质心运动定理,即:
或:
m|Xm2(xlvcoS—lsin2)=F-cx
重物做平面曲线运动,则可以直接利用牛顿定律得到切线方向的动力学方程:
m2(Pxcos^)=-m2gs2
(2)
(1),
(2)两式是耦合的非线性动力学方程。
当系统被限制在v-0附近运动时,可将其在v-0处线性化处理,则可以得
到系统的方程为:
(叶m2)xcx-m2l寸)=F叶乂m2(xl:
T妇2)=F
当给定F=F(t)时,可以建立仿真模型。
请读者考虑,如果要考虑摆杆的质量,则动力学方程如何?
例2-3:
车辆悬架系统的动力学模型
考虑图2.2所示的汽车悬架系统示意图。
设计悬架缓冲系统的k1,G;k2,C2
的目的是减小车辆在崎岖道路上行驶时产生的震动,因为道路表面的不平坦会引起悬架沿垂直方向的移动和绕某个轴的转动。
我们将整个系统的质量中心作为坐标的原点,因此系统在不平道路上的振动运
动可以看作是质心的沿垂直方向的平移运动以及绕质心的旋转运动。
车架质量为m,转动惯量为J。
输入车轮的位置信息y1、y2表明路况信息。
假设每个车轴的缓冲系统由具有阻尼特性的弹簧构成。
忽略轮胎的质量,每个
车轮受到的外力为弹簧弹力与阻尼力之和,即
d-
Fa=(Gki)yA(s)十皿ky)
dt
d-
FB=(C21,k2)yB(S)=(C2『B'k2『B)
dt
其中:
yAy-yiYby—b「—y?
yA和Yb分别表示每个弹簧距离参考位置的瞬时距离。
代入上式后
F^(Ci—ki)(ya,-yj
dt
或者:
my〉-G(ya;-yj-kjya:
-yj-C2(y-b「-y?
)-k2(y-b:
-y?
)
整理后得到:
my(c1c2)y(k|k2)y(^a-c2bp(«a-k2b):
=CiyiC2Y2kiyik2Y2
用y(t)和:
(t)分别表示系统质心的平移位移和沿质心的旋转角度。
上式中假定在很小的角度位置条件下满足sin「”「,并且「取顺时针的旋转方向
为正方向。
再根据系统相对于质心的动量矩定理可得:
d2:
J2=Fbbcos'-Faacos':
Fbb-F&adt2
其中J是车驾相对于质心的转动惯量,将上式整理后可得:
2
ddd
Jr=(C2k2)(y-b:
-y2)b-(Gki)(ya-yja
dtdtdt
或:
J「(c2bciap(ci-c2)y(k2b2kia2p(^a-k2b)y
二cay_c2by2kiay^k2by2
将系统的动力学方程写成矩阵形式:
mo_
112
BB
rl]L
+
Hl]J„yT-oJ
2
B1
cnc_
+
21
巳ErilHL
-
2c
2
2
c
2
已
2
2
E
1
n2FFr]L
F12FF
简写为:
其中:
[A]7[B]7[C]¥曰
0[BHC1c2
J&a-Qb
"ki
[F:
kia
y
W」
[A]
【A]=m
0
lc1
cia
E].|yiL[F]p]
“2.“2
db[eJi
CiaC2b||kia-k?
b
k2
-k2b
yi
k2
「k2bkiak?
b
C2
-C2b
7[A]」[B]7[A]」[C】_:
=[A]」[E]_y:
[A]」[F];
_V2
图2.5
当[A]为非奇异阵时,可以通过矢量信号我们可以得到系统的仿真模型如(图
2-5)。
Produetl
悬架系统仿真框图
;:
2,2,r,r,写成矩阵形式
r2sin
r2cos
以上系统中假定y1、y2是系统两个相互独立的输入变量,但实际上,后轮与前
车轮的位置时间相差△t=L/V时间。
这样,实际系统满足y2=%(t—fit)。
由于借助
了拉斯变换,将微分方程换成了代数方程,如果要得到时域响应则需要借助拉斯反变换。
根据第一章的基本知识,给出基于微分方程的仿真模型,具体计算过程留给读者练习。
例2-4机构运动学建模
曲柄滑块机构的运动学仿真建模(速度分析与建模)
曲柄滑块机构如图所示:
该机构只有一个自由度,首先给出机构的运动学分析模型,
(1)机构的封闭的矢量方程
r=riD
(2)矢量方程的分解式
r1cos1r2cos「2=r
r1sin1r2sin「2=0
(3)关于机构速度问题的运动学方程;
-r1sin〔〔-r2sin22=r
Acos1〔r2cos:
22二0
机构的输入运动量为>,,输出量为
1.J一5—
0.Ilr_-r1cos1
可以写成显式表达式
:
2_r2sin「2
_r|^2cos;:
2
Simulink仿真模型建立
在该仿真模型中,设系统的输入角速度为:
「1=150弧度/秒,通过一次积分可
以得到角度\,将这两个输入量通过一个信号混合器(以向量形式混合为一路信号),输入给MATLABFCN模块,通过该函数模块中的代码入,从而可以得到输
出量(;:
2,r'),再进一步积分后,得到位移量\(t),r(t)。
Compv,其它不变,
在MATLABFUNTION模块中写上函数过程文件名:
建立m脚本文件如下:
(函数子程序)
function]x]=compv(u);[x]输出,(u)输入。
%%参数说明:
r1曲柄长度,r2连杆长度
%%u
(1)曲柄角速度;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
滑快的速度与位移变化规律
-Ri(cos1sin1J「r2(cos:
2sin22)=r
rj—sincos11)r2(—sin2cos2\)=0
机构的输入运动量为S,喑,■:
1,输出量为:
:
2,2/2,r,r,r,写成矩阵形式
"r2sin®2
11
_[
-r1co^P^'2-r2co^P^'f-r1sin®^1
]r2cos®2
0_
r1sin鸣-r1cos^Q;+r2sin®2®:
P2]r2sin申21]」_口—r2cos申2弟;一口sin申retI
]r"-Ecos^0」[nsin督叫2-rrcos昭气+r2sin杓。
;一
和速度仿真一样,请读者建立机构的加速度仿真模型。
如果要对此机构的动力学仿真,可以再列写出系统的动力学方程,与运动学方程联立求解。
例2-5建立如下系统的振动微分方程,并使
用子系统封装技术。
=-k2(x2「xjc2(x2「xjK%「gm
m2X2=一k2(X2-xj-C2(X2-xjf(t)
改写上式为:
..1..
x[k2x2-k2xrc2x^c2x^k1x^c1x1]
m-i
设:
m^21kg,m2=9kg,O|=c2=2Ns/m,
kr=400N/m,
k2=600N/m,f(t)=sin(t)
利用子系统技术,我们可以建立相应的仿真模型,利用摸态分析方法可以得到系统的解析解和仿真解进行比较。
若将激励作用在左边质量块上,取f(t)二sin(5t),并分析当m2取值为多大时,质量m^j的振幅接近于零
(动力消振器原理)。
并进一步分析,当d=C2=0时,主系统的消振效果。
说明有阻尼消振效果好还是无阻尼消振效果好。
§2.3Hamilton动力学建模体系
方法,这些建模的基础理论有Lagrange第二类方程,Hamilton原理、Hamilton正则方程、APPEL方程和凯恩方程等.
1.Lagerange第二类方程
dt
其中:
L=T-V称为Lagerange函数。
这里:
T是系统的总动能,V是系统的总势能。
对于具有保守力作用和非保守力作用的混合系统,其方程为:
(2-
其中Qj是对应非保守力的广义力。
拉格朗日方程式是一组关于m个广义坐标的二阶微分方程,它有统一的格式和步骤,因此在动力学建立模型时经常采用。
2系统有耗散元件的拉格朗日方程
在工程实际问题中,如果存在有与速度有关的阻力。
例如当物体在空气、液体中运动时会受到流体介质的阻力作用。
实验表明,流动介质的阻力与相对速度有关,并且使系统的总能量不断减少。
这种阻力统称为耗散力,将这类元件统称为耗散元件。
作用于系统的耗散力一般可以表示为如下形式
-vi
Fi二—kjfi(Vi)」(i=1,2,,,,n)
Vi
其中vi表示第i个质点的速度,Fj表示第i个质点受到的耗散力,ki是阻力系数、fi(vj是与广义速度有关的函数,其中的负号表示阻尼与速度方向相反。
在系统中如果存在有耗散力时,只需将耗散力的广义力添加在拉格朗日方程的右边即可。
关于耗散广义力计算可参考下式:
根据广义力的定义
考虑到匕
v
—,kifi(Vi)--
Vi两
其中Vi•—巴
Bj
因此有
LVi
二為ki0fi(Vi)dVi
i1
:
L
这样容易得到具有耗散系统的拉格朗日方程为:
:
D
d/汀、汀:
U:
D门
或者:
()0
dt:
qj:
qj:
qj:
q」
因此对于耗散系统,只需将耗散力的广义力加进Lagerange方程的普通广义力中即可。
例如,在线性动力学系统中,一般当阻尼力是广义速度的一次式,即:
VVik2
F=-kv,则对应的耗散函数为:
D二kvdvv,对应的广义力为:
v^02
cD
Qkv。
&
例2-6一旋转摆如图所示,摆长为L,摆锤质量为m用光滑铰链连接在铅直
轴上,如果要考虑On构件的质量为M,当铅直轴以任意角速度转动时求出对应的动力学模型。
解:
当3为任意时,此时系统有两个自由度,
分别取
为:
「和二为广义坐标,其动能和势能分别
T=m(\22^2\2sin2:
)V二_mglcos:
Lagerange函数为:
L二T_V二m(l2:
2lS2sin2)mglcos:
在通常情况下,在转轴上作用有外加力偶矩M,根据Lagerange方程:
眷日今亠:
MyM
ml2「:
一;ml2'sin2:
「'川mglsin=0
以上两式仍为耦合非线性动力学方程。
1ED
D\l2,耗散力的广义力为:
Q一「
3Hamlton原理
Hamilton原理是以变分为基础的建模方法,设系统的动能为T,势能为V,
非保守力的虚元功为w,则Hamilton原理可以表示为:
t1
.(4•-w)dt=0其中:
L=T-V称为拉格朗日函数
to
Hamilton原理常用来建立连续的质量分布和连续刚度分布的系统(弹性系统)的动力学模型。
例2-7弹性系统的动力学建模
所谓的弹性系统是指具有连续的质量分布和连续刚度分布的系统,下面通过梁的横向振动来说明弹性体的建模方法。
o
1
I
A
I
/
Fy(x,t)一^^2
\x
y
1
IRwImRII
F
设梁的长度为l,截面的弯曲刚度EI为常数,单位长度质量为匸,在x截面形心处横向位移为y(x,t),忽略剪切变形,
xIrv.rv.rv.
【1I•…厂'''•-''
t[0「(Vy)--yy-(Elyy)(Ely)]dxdt
t00:
txx
根据Hamilton原理,满足时间端点的条件当:
t=t0和t二鮎时有:
y(to)」y(ti)=0
tic2”,
于是我们可以得到:
["•-y(Ely)]'ydxdt二0
to0:
x
根据-y的任意性,满足上式条件为:
_2_2_2_
'y2(EIy)=0[=(Ely")、y—Elyly'];=0
.:
t;:
x;:
x;x
第一式为梁的自由振动方程,第二式是变分问题中自然满足的边界条件。
可以
使用模态分析方法,将偏微分方程化为常微分方程,然后就可以利用前面的方法来建立数学模型。
当梁上作用有分布载荷力和分布力偶时,如下图:
贝打系统的虚功可以表示为:
如在梁上某点a处作用集中力P和b点处作用有集中力偶矩M时,这时,其
右边的广义力可以表示为:
(M、(x-b))P-;(x_a),和
ex
并注意到:
m(x—b))M(x-b)
-X
在一般情况下,一个连续系统的动态特性可以用一个高阶微分方程或微分方程组来表示;
n-1n-2
Cn_|U
du丄du
anyF77-ci廿.
dtndt
(2-1)
其中:
y是系统的输出,u表示系统的输入量,如果引进微分算子
dt
则有:
a°pnyaiPn'y…色畀pn'uqp^u…c」
n4n-1
即:
'、an4_jPjy二為Cn^P^u
j£j=0
一个动力学系统的数学模型建立起来以后,还需要对该系统响应规律进行分析,
以便揭示真正的运动规律。
或者通过建立仿真模型来揭示运动规律。
§2.4一维弹性体有限元建模
有限元的基本思想是先把结构分割成N个不同单元,分别对单元和节点编号
1,2….N。
单元划分越细,计算精度越高,但是计算工作量也越大,因此,要根据具体情况合理的划分单元数,本节将介绍一维梁单元有限元建模方法。
1
2
N
2.5.1梁单元质量矩阵与刚度矩阵
设梁单元中的第i个单元的坐标xe(局部坐标),单元长度为丨,该单元有两个节点,而每个节点有两个广义坐标,这样一个梁单元共有4个广义坐标,分别的左界面的位移qe1,qe2与转角和右截面的位移和转角qe3>qe4,有:
Wgt),
qe1一W(Xe,t)|笔出;qe2|xe=0,
cw(xe,t)
qe3-w(xe,t)|xezt,qe4|xe=t
欣e
32
设单元的位移模式为W(Xe,t)=&Xe'C2xe-C3XeC4
将单元边界条件带入上式,可得:
1
C1=13(2qe1Iqe2—2qe3Iqe4)
1
C2二F(—3qe1—3lqe23qe3Tqe4)
qe1
整理后可得:
W(Xe,t)=[:
:
1(Xe):
:
1(Xe):
:
1(Xe)1(xe)]^2
qe3
qe4
3223223
其中:
1(Xe)=(1卡),2(Xe^(X--^令),
C2r323
f"3Xe丄2XefnXe丄Xe
3(Xe)=(t2e门,4(Xe)十;^2)
设梁的单位长度质量为’,系统的动能为
11dw2144
一2°h)dx云:
严qej
其中:
mji(X^j(Xe)dXe
_156
221
54
-131〕
PI
412
131
-312
可得单兀质量矩阵为:
Me=
420
156
-221
4l2一
系统的势能为;
11(W214
U=2oEI^xe)dXe^L
1
其中:
kj=0EIl(Xe)[(Xe)dXe
功表达式为:
2.5.2总体系统动力学微分方程:
以上仅仅给出了单元系统的微分方程,通过个单元的对接条件,我们可以得到
总体坐标下的动力学微分方程,为了得到总体坐标系中的动力学方程,先引入总体节点位移向量:
{q=[qiq2…qn]T对于两单元,有:
n=2(N1)个位移分量,与单元节点位移向量,q;=01qe2qe3TedT
设局部位移向量与总体位移向量的关系为:
{q}ej=[s]{q},i=1,2…N
则系统的总动能为:
1NiN
g]T[Me佔心鈕曲
NN
其中:
Mei=[Si]T[Mei][S]
其中:
Qi=[S]T[Kei][Si]
T~
得:
M八[Si][Mei][Si]八M
ii
NN
T~
同理:
K=7[S][Kei][S]=
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 第二 章动 力学 系统 微分方程 模型