第四章 多自由度系统的动力学特性
多体系统的动力学特性研究

多体系统的动力学特性研究多体系统的动力学研究是物理学中一个关键领域,涵盖了许多重要的科学和工程应用。
这些系统由许多相互作用的自由度组成,其行为具有复杂性和非线性特性。
在本文中,我们将探讨多体系统动力学研究的一些重要方面,并介绍一些常见的方法和技术。
首先,我们需要了解多体系统中的动力学行为如何受到它的微观结构和相互作用的影响。
这包括粒子间的相互作用力、碰撞、传输过程等。
在许多实际的应用中,我们经常需要研究领域特定的多体动力学模型,如分子动力学、固体力学、流体力学等。
研究多体系统的动力学特性的一个重要方面是探索系统的宏观行为和微观结构之间的关系。
这种关系通常通过建立连续力学模型来实现,例如通过偏微分方程来描述宏观行为。
通过将微观信息转化为宏观描述,我们可以更好地理解系统的非线性行为和相变现象。
在多体系统的动力学研究中,统计力学是一种非常重要的方法。
统计力学研究的是大量微观粒子组成的系统,利用概率分布函数来描述微观状态的出现概率。
统计力学可以解释系统的平衡态和非平衡态,并为系统的动力学性质提供了重要的理论基础。
基于统计力学的方法可以用来计算系统的热力学性质、输运性质和相变等。
另一个重要的多体动力学研究方法是计算模拟。
计算模拟利用计算机来模拟多体系统的运动和相互作用。
通过数值算法和计算技术,我们可以模拟和预测不同尺度下的多体系统的行为。
计算模拟方法已经被广泛应用于材料科学、生物物理学等领域,提供了对复杂系统行为的深入理解。
除了统计力学和计算模拟,实验方法也是多体系统动力学研究中不可或缺的一部分。
实验方法可以用于测量和验证理论模型的预测结果,并为理论研究提供实验数据。
通过实验观察和测量,我们可以获得关于多体系统行为的定量信息,从而更好地理解系统的动态特性。
总之,多体系统的动力学特性研究是一个宽广而充满挑战的领域。
通过深入研究多体系统的微观结构和相互作用,建立宏观描述模型,利用统计力学、计算模拟和实验方法进行研究,我们可以获得对系统行为的深入认识。
多自由度振动系统的动力学模型构建

多自由度振动系统的动力学模型构建引言:多自由度振动系统是指由多个自由度的质点组成的系统,在这样的系统中,每个自由度都可以独立地进行运动。
动力学模型的构建是研究和理解振动系统行为的基础。
本文将介绍多自由度振动系统动力学模型的构建方法及应用。
一、质点模型多自由度振动系统的最基本组成单位是质点。
质点的运动可以用坐标形式以及质点的质量、刚性等参数来描述。
对于一个有n个自由度的振动系统,可以通过将每个自由度的质点模型相连接构成整个系统。
二、约束关系与广义坐标在多自由度振动系统中,质点之间相互约束,其运动不再是自由的,而是受到约束的影响。
为了描述约束关系,引入广义坐标来表示系统各个自由度的相对运动。
广义坐标是将实际坐标通过约束条件变换得到的坐标表示。
三、拉格朗日方程与振动方程拉格朗日方程是多自由度振动系统的基本动力学方程。
通过对系统的动能和势能进行推导和求导,可以得到描述系统运动的拉格朗日方程。
对于振动系统而言,通过求解拉格朗日方程,可以得到系统的振动方程,进一步描述系统的运动行为。
四、模态分析与特征频率模态分析是研究振动系统固有特性的方法。
对于多自由度振动系统,可以通过模态分析得到系统的固有模态和特征频率。
固有模态是指系统在自由振动时,各个自由度的振动模式。
特征频率是指系统在不同固有模态下的振动频率。
五、系统的耦合与动态响应多自由度振动系统中的各个质点之间存在耦合关系,一个自由度的振动会对其他自由度的振动产生影响。
通过研究系统的耦合关系,可以得到系统的动态响应。
动态响应是指系统对外界激励的响应行为,可以通过求解振动方程得到。
六、应用案例:建筑结构振动多自由度振动系统的应用广泛,尤其在建筑结构的振动研究中起到了重要作用。
通过对建筑结构的多自由度振动系统进行建模和分析,可以评估结构的稳定性、抗震性能等。
振动模型的构建和分析可以提供设计和改进建筑结构的依据。
结论:多自由度振动系统的动力学模型构建是研究振动系统行为的关键步骤。
结构动力学4-1

(−ω 2 [M ] + [K ]){φ }sin(ωt + θ ) = {0}
因为sin(ωt + θ)为任意的,可以消去,因此,
([K ] − ω [M ]){φ } = {0}
2
上式是关于{φ}的N阶齐次线性方程组,表征了振型和自 振频率的关系 ,称为运动方程的特征方程。 由特征方程可解得自振频率ω和振型{φ}。
1
k22=1800
k23=-600
(c)
(d)
1 多自由度体系的自振振型和自振频率
算例1 结构的质量阵、刚度阵:
1.0 u 1=1 u3 600 1.5 1200 2.0 1800 (a) (b) u1 1 k11=3000 k12=-1200 k13=0 u2 k21=-1200 k31=0 k 32=-600 1 k 33=600 u2=1 u3=1
1 多自由度体系的自振振型和自振频率
算例1 运动方程的特征方程:
0⎤ ⎡ 2. 0 0 ⎢ 0 1. 5 0 ⎥ [M ] = ⎢ ⎥ ⎢0 0 1 .0 ⎥ ⎦ ⎣ 0 ⎤ ⎡ 3000 − 1200 ⎢− 1200 1800 − 600⎥ [K ] = ⎢ ⎥ ⎢ 0 − 600 600 ⎥ ⎦ ⎣
算例1 如图(a)所示三层框架结构,各楼层的质量和层间 刚度示于图中,确定结构的自振频率和振型。 结构模型及各刚度元素:
1.0 600 1.5 1200 2.0 1800 (a) (b) u1 1 k11=3000 k12=-1200 k13=0 u2 k21=-1200 u 1=1 u3 k31=0 k 32=-600 1 k 33=600 u2=1 u3=1
多自由度机械系统动力学

Chapter4多自由度机械系统动力学
机械动力学
Chapter4
多自由度机械系统动力学
2021年6月18日
机械动力学
Chapter4多自由度机械系统动力学
本章解决的主要问题及内容
解决的问题: 解决两自由度机械系统的动力学问题。采 用方法为拉格朗日方程的分析方法。
主要的内容:
一、拉格朗日方程;
工程中的非自由质点系,受到的约束大多是稳定的完整 约束(约束方程仅与质点系的位置有关)。
确定一个受完整约束的质点系的位置所需的独立坐标的数 目,称为该质点系的自由度的数目,简称为自由度数。
对一个非自由质点系,受s个完整约束,(3n-s )个独
立坐标。其自由度 为 N=3n-s 。
机械动力学
Chapter4多自由度机械系统动力学
机械动力学
Chapter4多自由度机械系统动力学
例:铅直平面内摆动的双摆。
▼确定A、B两点位置(平面问题) 需四个独立坐标 ▼系统受两个完整约束,其约束方程:
x12 y12 a2 , (x2 x1)2 ( y2 y1)2 b2
▼系统的自由度:N=2n-s=4-2=2
★两个自由度, 取广义坐标,
Qk 0 (k 1,2,, N )
机械动力学
Chapter4多自由度机械系统动力学
以广义坐标 表示的质点系的平衡条件:
Qk
n
(Xi
i 1
xi qk
Yi
yi qk
Zi
zi ) 0 qk
(k 1,2,, N)
解决质点系的平衡问题的关键是如何计算广义力
※广义力的计算
方法1:计算广义力 Qk 的步骤
N
xi
多自由度机械系统动力学

例:图示系统中,杆OA和AB以铰链相连,O端为圆柱绞, B端自由,杆重及摩擦不计,杆长 OA=l1,AB=l2,设二 杆均在铅垂面内,OA杆与铅垂线成φ1角,杆AB与铅垂 线成φ2角.今在点A和B分别作用铅垂向下的力F1和F2, 求在图示位置时的广义力。
解:1、定义法求广义力 此为具有二个自由度的双摆系统,选取φ1和φ2为广义 坐标,对应的广义虚位移为φ1和φ2,由定义得:
抖 yA Q1 = F1 + F2 抖 f1 抖 yA Q2 = F1 + F2 抖 f2 yB f1 yB f2
因
y A = l1 cos f 1 , y B = l1 cos f 1 + l2 cos f 2
求出相应的偏导数,代入广义力公式有:
Q1 = - ( F1 + F2 )l1 sin f 1 Q2 = - F2l2 sin f 2
2
cos ) l 2x
系统势能:(以弹簧原长为弹性势能零点,滑块A所在平面 为重力势能零点)
1 2 U kx m2 gl cos 2
拉格朗日函数:
L T U 1 1 2 m2l 2 (m1 m2 ) x 2 2
2
1 2 l cos kx m2 gl cos m2 x 2 L L (m1 m2 ) x m2l cos , kx x x
2、用虚功方法求Q1和Q2,可先令φ2=0,可得: 由于 d rA = d rB = l1df 1 代入上式得: Q1 = - ( F 1 +F 2 )l1 sin f 1 再令φ1=0,可得:
' d WF - F1 sin f 1d rA - F2 sin f 1d rB Q1 = = df 1 df 1
多自由度系统振动

的方法。
传递矩阵法适用于线性时不变系 统,能够处理多自由度系统的振
动问题,计算效率较高。
传递矩阵法的精度取决于系统参 数和边界条件的准确性,对于复 杂系统和非线性问题,需要采用
其他方法进行求解。
模态叠加法
模态叠加法是一种基于模态展开的数值 计算方法,通过将系统的振动表示为一 系列模态的线性组合,求解每个模态的
振动方程,得到系统的动态特性。
模态叠加法适用于线性时不变系统,能 够处理多自由度系统的振动问题,计算
精度较高。
模态叠加法需要选择合适的模态数目和 模态提取方法,对于大规模系统和复杂
未来研究方向
深入研究多自由度系统振动的 非线性特性,探索更精确的数
学模型和数值模拟方法。
针对复杂多自由度系统,研究 多因素耦合振动和多场耦合振
动的理论和方法。
发展多自由度系统振动主动控 制和智能控制技术,提高系统 振动控制精度和响应速度。
将多自由度系统振动理论应用 于实际工程领域,解决重大装 备和结构的振动问题,提高其 稳定性和安全性。
THANKS FOR WATCHING
感谢您的观看
02
它涉及到多个振动子之间的相互 作用和耦合,其动力学行为比单 自由度系统更为复杂。
研究背景和意义
随着科技的发展,多自由度系统在许多领域中得到了广泛应用,如大型机械装备、 精密仪器、高层建筑等。
由于多自由度系统在受到外部激励或内部参数变化时,会产生复杂的振动行为,这 不仅会影响系统的性能和稳定性,还可能引发安全问题。
航天器振动控制
总结词
结构动力学多自由度

pbT
~ fpa
paT ~fpb paT ~fpb T pbT ~f T pa pbT ~fpa
故 ~f 、 k 均为对称矩阵。
单元刚度矩阵
单元刚度系数表示由单位节点位移所引起的节点力。
单元刚度系数由虚位移法求得。
例如,课本P106图11-5所示简支梁中,令a端发生单位转角, 并给该处一竖向虚位移,零外力所做的功,等于内力所做的 功。
表示一个自由度发生相应单位位移而其他节点不动时在结构中所 产生的的力。
弹性特性
柔度的定义:
~ fij —在j坐标施加单位荷载而引起的i坐标的挠度。
则任意荷载组合下: vi ~fi1 p1 ~fi2 p2 ~fiN pN
用矩阵表示:
v1
vi
v N
略去阻尼矩阵和施加的荷载向量的影响: mv kv 0
假定以上多自由度体系的振动是简谐振动:
v(t) vˆ sin(t )
vˆ 表示体系的形状,不随时间变化。
v 2vˆ sin(t ) 2v 2mvˆ sin(t ) kvˆ sin(t ) 0
k 2m vˆ 0
无阻尼自由振动—振动频率分析
k 2m vˆ 0
即: k 2m 0
上式的N个根,表述体系可能存在的N个振型的频率。
1
2
3
WE va pa v1 k13
Lபைடு நூலகம்
WI v1 0 EI ( x) 1''( x) 3''( x)dx
多自由度系统的振动、响应和求解

P
f Pb(xl2x2b2), 0xa 6EIl
x
a
b
l
f Pb[l(xa)3(l2b2)xx3], axl
6EIlb
例4.1 写出图示梁的柔度矩阵,梁的抗弯刚度为EI。如果 将梁的质量按分段区间均分到区间的两个端点,写出梁的质
量矩阵,设梁单位长度的质量为 l。
0 0 y y2 1 m 3 y3
qrki qrkj qiqj 12in1
n
mijqiqj
j1
其中mij
N
mk
k1
rk rk qi qj
பைடு நூலகம்
mji
(4.4)
(4.2)、(4.4)式可写成矩阵形式
V1qTKq, 2
(4.5)
T1qTMq
(4.6)
2
其中q[q1, q2,, qn]T,K[kij]nn,M[mij]nn
矩阵K 称为刚度矩阵,它是一个对称正定或半正定矩阵;矩 阵M 称为质量矩阵,它是一个对称正定矩阵。
因此,柔度矩阵的第一列为
{f11,f21,f31}T76l83EI{9,11,7}T
类似可算出柔度矩阵的第二、第三列。柔度矩阵为
f11 [F]f21
f12 f22
ff123376l83EI191
11 16
7 11
f31 f32 f33
7 11 9
系统的动能为
T1 2(m 1y & 1 2m 2y & 2 2m 3y & 3 2)1 2{y1,y2,y3} m 01 m 02 0 0
§4.1 多自由度系统的动力学方程
我们先来考察多自由度线性系统动能和势能的数学结
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
q
T
K q [ ]
T
i
K
第i阶特征值 i { X } T [ M ]{ X } M i i
{ X } i [ K ]{ X } i
K
i
i
K
i
M i i M i i
2
展开定理
方程 主振型
X j M
T
X i T X j K X i
主振型关 于质量矩 阵和刚度 矩阵正交
两式相减 得
i
j
X M X
T i
j
0
X iT M X j 0 X iT K X j 0
X iT M X i M i X iT K X i K i
x1 (1) x 2
(1)
A1 (1) A 2
(1)
sin( ω 1 t θ 1 )
x1 (2) x 2
(2)
A1 (2) A 2
(2)
sin( ω 2 t θ 2 )
1 (1) 2
(1)
Φ
(1 )
1 (1) β
Φ
(2)
1 1 ( 2 ) (2) β 2
(2)
P (t) A 1 sin( ω 1 t θ 1 ) P (2) (2) A sin( ω t θ ) P (t) 1 2 2
k2 k2 k3 k3
x1 k3 x2 k 3 k 4 x3 0
F A e j t 0 0
一、两自由度系统的自由振动
(6)主振动的确定:
系统以某一阶固有频率按其相应的主振型作振动,称为 系统的主振动。故第一和第二阶主振动分别为
系统作主振动时,各广义坐标同时经过静平衡位置和达到最 大偏离位置,是一种有确定的频率和振型的简谐振动
一、两自由度系统的自由振动
(7)一般情况下,自由振动的通解为:
上述两种主振动的叠加,即
x1 x 2 x1 (1) x 2
(1) (1)
x1 (2) x 2
(1) (1)
一、两自由度系统的自由振动
一般情况下,系统的自由振动是由两种不同频率的主振动的 线性组合,其结果不一定是简谐振动,这是和单自由度系统 自由振动有很大区别的地方。 只有在特定的初始条件下,其自由振动才表现为单一振型 的主振动,才是一种简谐振动,如: – A 1( 2 ) 0 ,系统作第一主振动; – A (21 ) 0 ,系统作第二主振动; – 其他;
2
上式为广义特征值问题
KA ω MA
2
方程有非零解的条件:
K M 0
2
二、多自由度系统的自由振动
即: k 11 m 11 2 k 21 m 21
2
k 12 m 12
2 2
k 1n m 1n
2 2
k 22 m 22 k n2 m n2
j
1 2
T q Cq
j 1
L R 代入拉格朗日方程: d L Q dt q
j
q
j
q
F j
(i 1,2, ,n)
j
得
M q ( t ) C q ( t ) Kq( t ) Q
F
(t )
多自由度系统的运动是由一组二阶常微分方程描 述的,多采用数值法求解。 如:直接积分法、有限元法和无限元法。
设有以下解:
q j A j sin( t )
( j 1, 2 , , n )
二、多自由度系统的自由振动
矩阵形式:
q A sin( t )
其中 A [ A A A ] T 为振型向量 1 2 n 代入动力学方程得:
K
ω M A 0
T
正则化的振型矩阵
q 1 M
1
X
1
, ,
1 M
2
X
2
, ,
正则化振型矩阵的正交性
q
T
M q I [ M ]
模态刚度、广义刚度或主刚度
K i { X } i [ K ]{ X } i
T
正则化振型矩阵的正交性
二、多自由度系统的自由振动
多自由度系统,无阻尼自由振动方程式的一般形式:
M q ( t ) Kq( t ) 0
展开式为:
m 11 m n1 m 1 n q1 k 11 m nn qn k n 1 k 1n q 1 0 k nn q n 0
固有振型示意图
两质量-弹簧系统
三质量-弹簧系统
m 1 0 0 m2 0 0
0 c 1 c 2 - c2 0 x1 k 1 k 2 x1 0 - c 2 c2 c3 - c3 x2 k2 x2 m 3 3 0 - c3 c 3 c 4 x3 0 x
一、两自由度系统的自由振动
(5)将 代入下式,对应于 n 1 和 n 2 , 分别得到振幅A1和A2之间的两个确定的比值:
2
ni ( i 1, 2 )
2
2
(1 )
A A
(1 ) 2 (1 ) 1
(2)
A
(2) 2 (2)
A1
一、两自由度系统的自由振动
当系统按某一固有频率振动时,振幅比只取决于系统本身 的物理性质,而与初始条件无关 在振动过程中,系统各广义坐标的位移之相对比值可以由 该振幅比确定 该比值确定了整个振动系统的振动形态,称为主振型或固 有振型
一、两自由度系统的自由振动
无阻尼两自由度振动系统的运动微分方程式
以静平衡位置为坐标原点,设质量M1和M2的广义坐标为 x1和x2,并假设它们足够小,以保证系统在线性范围内运动; 得到以下二阶常系数线性齐次常微分方程组
一、两自由度系统的自由振动
更一般形式:
求解步骤: (1)假设简谐形式的解: 设振动时两质量块按相同频率和相位作简谐振动:
二、多自由度系统的自由振动
每个特征值 i 由小到大按序排列为:
1 2 n 1 n
多自由度系统的固有频率也是由系统本身的物理参数决定, 与起始运动状态无关;多自由度系统有多个固有频率,其 中的最低固有频率称为系统的基频;
振型矩阵 (模态矩阵) n自由度系统 ,有n个主振型{ X } i ( i = 1, 2, …, n ),振型矩阵为
q { X
}1
{ X }i
{ X }n
特征值矩阵 (频率矩阵)
n个特征值组成的对角矩阵,称为特征值矩阵或谱矩阵。
1 0 2 ω λ 0 0 0 n 0
(2)
1 A 1 sin( ω 1 t θ 1 ) (2) (2) β A 1 sin( ω 2 t θ 2 )
(1)
Φ
(1)
P (t) (2) Φ (2) ΦP P (t)
(1)
Φ 是模态矩阵,P是主坐标。
i j
i
X iT K X i X iT M X i
Ki M
i
第i阶主刚度 第i阶主质量
振型矩阵
q { X
T
}1
{ X }i
{ X }n
q
M q M
0 M
0 K
2
q
2
T
K q K
第四章
多自由度系统(MDOF)的动力学特性
用一个独立坐标描述的单自由度系统,是实际振动系 统的最简单模型; 用两个或更多有限个独立坐标描述的振动系统称作多 自由度系统; 工程上各种机械的结构物,总是由杆、梁、板、壳等 元件组成的弹性体,它们的质量与刚度都具有分布的 性质,理论上是无限自由度系统,然而在多数情况下, 无限自由度问题可以简化为有限多个自由度系统进行 研究;
2
k 2n m 2n k nn m nn
2
k n1 m n1
2
0
展开得到系统的特征方程:
2n
a 1
2
2 ( n 1 )
a n 1 a n 0
2
该特征方程存在 的n个正实根,即系统的特征值,也叫系 统的固有频率。固有频率对应的特征向量{X}i叫固有振型或 主振型。
0 0 M n
0 0 Kn
模态质量矩阵
M 1 0 M 0
K 1 0 K 0
0
0
模态刚度矩阵
模态质量、广义质量或主质量
M i {X }i M {X }i
一、两自由度系统的自由振动
(2)将上述假设解代入运动方程,得代数特征值问题:
(3)该线性齐次代数方程组,非零解的条件是系数行列式为零 :
2 2 —— n 的特征方程,是 n 的二次多项式,又称频率方程。