多自由度系统的动力学方程
结构动力学之多自由度体系的振动问题
2.760 3.342 1
0.163
0.924
2.76
柔度法
利用刚度法的方程间接导出柔度法方程：
由刚度法振幅方程：
令λ=1/ω2 得频率方程：
( [K]－ω2 [M] ){Y}={0}
前乘[K]－1=[δ]后得： ( [I ]－ω2 [δ] [M] ){Y}={0} ( [δ] [M] － λ [I ] ){Y}={0} ┃ [δ] [M] － λ [I ] ┃=0
刚度法
2）如果初始条件是任意的，则任其自然 后， 系统所发生的振动就不是按主振型的简谐自由 振动，而是复杂的周期振动，这时可以用各阶 主振动的线性组合来描述它，也就是说其通解 表为各个特解之和，即
y j sin( j t v j )
j 1 n
所以系统的任意振动可以表示为各个主振动 的叠加。
Yij为正时表示质1 1.293 5Y11 6.70Y21 3 0 量mi的运动方向与单 3Y 1.707 0
21
Y
(1)
0.163 0.569 1
0.569
5Y13 5.027Y23 3 0 (1) Y 3Y21 10.027 0 3.342 1.227
1 1 4 0 , m m 2 9
展开得: 解之:
3 15 2 42 30 0
ξ1=11.601，ξ2=2.246，ξ3=1.151
1 m
三个频率为：
1 0.2936
1 1 3 0.9319 m m 3）求主振型： （令Y3i=1）将λ1代入振型方程： （[δ] [M ]－λ1[I]）{Y}＝0的前两式：
第二章(多自由度系统的运动微分方程)详解
k11 k 21 kN1
k1 j k2 j k Nj
k1N k2 N k NN
刚度影响系数 kij ：第 j 个自由度产生单位位移，其他自由度位移为零时， 需要在第i 自由度处沿着位移方向施加的力。
用影响系数法建立系统的运动微分方程
【例】用影响系数法写出图示系统的刚度矩阵。
多自由度振动系统
Piezoelectric actuator
基于压电作动器的垂尾抖振主动抑制 （此系统有一、两千个自由度（3D实体单元） ）
Z Y
X
第二章： 多自由度系统的运动 微分方程
第二章:多自由度系统的运动微分方程
第一讲：
1.建立多自由度系统运动微分方程的 各种方法的概述 2.用牛顿第二定律列写系统的运动微 分方程 3.用影响系数法建立系统的运动微分 方程
F1 1
k3
m2
k2 (d11 d21 )
m1
k2 (d11 d21 ) k1d11 1
d 21 k2 (d11 d21 )
F2 0
d11
k3d21
k2 k3 k1k2 k1k3 k2 k3 k2 k1k2 k1k3 k2 k3
m2
d 21
k2 (d11 d21 ) k3d21 0
上次课内容回顾
3.刚度影响系数
刚度影响系数 kij ：第 j 个自由度产生单位位移，其他自由度位移为零时， 需要在第 i 自由度处沿着位移方向施加的力。
4.柔度影响系数
柔度影响系数 dij ：第 j 个自由度上作用单位力，其他自由度作用力为零时，
在第 自由度上产生的位移。 i
5.刚度矩阵和柔度矩阵的关系
多自由度机械系统动力学
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
结构动力学多自由度系统振动
运用功旳互等原理可知，刚度矩阵是对称阵，即有kij＝kji， 于是上述刚度矩阵为：
k1 k2
k2
K 0
0
0
k2 k2 k3
k3 0
0
0 k3 k3 k4 k4
0
0 0 k4 k4 k5 k5
0
0
0
k5
k5
⒉ 柔度法 柔度系数aij定义为：
在第j个质量上作用单位力时在第i个质量上产生旳位移。
K12 k2 K22 k2 k3
K32 k3 K42 0 K52 0
K13 0 K23 k3 K33 k3 k4 K43 k4 K53 0
K14 0 K24 0 K34 k4 K44 k4 k5 K54 k5
K15 0 K25 0 K35 0 K45 k5 K55 k5
(a) m1 mi
mj mn
y1
yi yj yn
m1 y1
(b)
mi yi
1
i
j
m j y j
mn yn
ii
ji
1
(c)
ij
ij
jj
(a) m1
mi
mj mn
y1
yi yj yn
m1 y1
(b)
mi yi
1
i
j
m j y j
mn yn
ii
ji
1
(c)
ij
ij
jj
于是： 若在第j个质量上作用有力F，则在第i个质量上产
2
2
2
1 Mx 2 1 m[x 2 2Lx cos L2 2 ] 1 kx2 mgL(1 cos)
2
2
2
d dt
多体系统的动力学分析
多体系统的动力学分析动力学是研究物体的运动及其产生的原因的学科，对于多体系统的动力学分析，我们需要探究不同物体之间的相互作用以及它们的运动规律。
在这篇文章中，我们将介绍多体系统的动力学分析方法，以及它在不同领域的应用。
1. 多体系统的描述多体系统是由多个物体组成的系统，物体之间可以通过各种相互作用力进行作用。
为了对多体系统进行动力学分析，我们首先需要对每个物体的位置、质量、速度等进行描述。
在经典力学中，可以通过使用牛顿第二定律 F = ma 来描述物体的运动，其中 F 是物体所受的合外力，m 是物体的质量，a 是物体的加速度。
2. 多体系统的相互作用在多体系统中，物体之间可以通过万有引力、电磁力、弹性力等多种相互作用力进行作用。
这些相互作用力是决定多体系统运动规律的重要因素。
在进行动力学分析时，我们需要考虑物体之间的相互作用力，并利用牛顿定律求解物体的运动轨迹。
3. 动力学分析方法在对多体系统进行动力学分析时，我们可以采用多种方法来求解物体的运动规律。
其中，最常用的方法之一是利用微分方程求解。
我们可以根据牛顿第二定律及物体之间的相互作用力建立运动微分方程，然后通过求解微分方程得到物体的位置、速度、加速度的函数关系。
另外，还有一些其他的动力学分析方法，如拉格朗日方法、哈密顿方法等。
这些方法可以根据系统的自由度来建立系统的拉格朗日函数或哈密顿函数，并利用变分原理求解系统的运动方程。
4. 多体系统的应用多体系统的动力学分析在物理学、工程学、天文学、生物学等众多领域都具有重要应用。
在物理学中，通过对多体系统的分析，可以研究宏观物体的运动规律，如行星运动、机械振动等。
在工程学中，动力学分析可以用于设计复杂结构的机械系统、车辆运动仿真等。
在天文学中，动力学分析可以研究星系、恒星运动，以及天体之间的相互作用。
在生物学中，动力学分析可以用于模拟生物体的运动、神经信号传递等。
总结：多体系统的动力学分析是研究物体运动及其相互作用的重要工具。
结构动力学多自由度
▪ 振型方程：
(K i2M)ji 0 (i 1, 2, 3, n)
▪∵
K 2i M 0
▪ ∴ 第i 个振型方程中的n 个方程中只有n-1个是独立的！ ▪ ——无法得到j1i、 j2i、 … 、 jni 的确定值， ▪ 但可以确定各质点振幅之间的相对比值： ▪ —— 振型的幅值是任意的，但形状是惟一的。
一致质量矩阵：
L
pava m13v1 0 fI ( x)v( x)dx
L
0
m( x) 3( x)v3
L
1( x)v1dx
mij 0 EI ( x)i ( x) j ( x)dx
L
cij 0 c( x) i ( x) j ( x)dx
其中，c(x)表示分布的粘滞阻尼特性。
一致节点荷载
L
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
N
2)
2
)
y32
(t
)
jˆ
32
s
in
(
2
t
2
)
1
jˆ
2i
yi
(t
)
jˆ3i
s
in(i
t
i
)
jˆ ni
1
jˆ 21
jˆ 31
jˆ 32
1
jˆ 22
将N个振型中的每一振型形式，用F表示N个振型所组成的方阵。
11 12 13 1N
机械动力学-多自由度系统
所有坐标 q j ( j = 1, 2, , n ) 的运动有着相同的随时间变化 规律，即有着相同的时间函数。 规律，即有着相同的时间函数。令 q j (t ) = u j f ( t ) j = 1, 2, , n （4.2-3） ）
u j ( j = 1, 2, , n ) 是一组常数，表示不同坐标运动的大小。 是一组常数，表示不同坐标运动的大小。
12
4.2 无阻尼自由振动和特征值问题
2. n自由度无阻尼系统自由振动的解的形式 ([ K ] − λ [ M ]) {u} = {0} （4.2-10） ） [ K ]{u} = λ [ M ]{u} （4.2-11） ） [K ] − λ [M ] = 0 （4.2-12） ） 系统的特征行列式，其展开式叫系统的特征多项式 系统的特征行列式， 方程（ 方程（4.2-12）叫做系统的特征方程或频率方程， ）叫做系统的特征方程或频率方程， ωn2的n阶方程 是 λ或 阶方程 2 λ1 < λ2 < < λn ωn21 < ωn22 < < ωnn 方根值 ωn1 < ωn 2 < < ωnn 叫做系统的固有频率，由方程（ 叫做系统的固有频率，由方程（4.2-12）可见，它只 ）可见， 决定于系统的物理参数，是系统固有的。 决定于系统的物理参数，是系统固有的。最低的固 有频率叫做系统的基频或第一阶固有频率， ω 有频率叫做系统的基频或第一阶固有频率，在许多 实际问题中，它常常是最重要的一个. 实际问题中，它常常是最重要的一个.
2
(
)Leabharlann U=对于线性系统，运动是微幅的， 对于线性系统，运动是微幅的，sin θ ≈ θ 代入动能和势能方程， 代入动能和势能方程，有
第2章——多自由度系统的振动——运动方程建立方法0425
船体振动基础1第章多自由度系统的振第2章多自由度系统的振动一、引言二、两自由度系统的振动三、多自由度系统的振动四、振动方程建立的其他方法2有阻尼的多自由度系统振动1、拉格朗日方程式1、拉格朗日方程式P38拉格朗日法是建立微分方程一种简单的方法：先求出系统的动能、势能，进而得出质量矩阵和刚度矩阵.优点：系统的动能和势能都是标量，无需考虑力的方向。
141、拉格朗日方程式P38拉格朗日第二类方程式适用于完整约束的系统。
完整约束完整约束：当约束方程本身或约束方程通过积分后可以下式所示的形式表示时，称为完整约束。
不完整约束：当约束方程本含有不能积分的速度项时，系统的约束称为不完整约束。
具有不完整约束的系统，系统的自由度不等于广义坐标数自由度数小于广义坐标数于广义坐标数，自由度数小于广义坐标数。
151、拉格朗日方程式P3811•位移方程和柔度矩阵P40对于静定结构，有时通过柔度矩阵建立位移方程比通过对于静定结构有时通过m1x1x2以准静态方式作用在梁上。
梁只产生位移（即挠度），不产生加速度。
的静平衡位置为坐标P1=1 f11 f21 f12P2=1 f22（1）P1 = 1、P2 = 0 时 m1 位移：x1 = f11 m2 位移：x2 = f 21 （3）P1、P2 同时作用 m1 位移： 位移 x1 = f11 P 1 + f12 P 2 m2 位移：x2 = f 21 P 1 + f 22 P 2（2）P1 = 0、P2 = 1 时 m1 位移：x1 = f12 m2 位移：x2 = f 22P1 m1 x1 x2 P2 m2P1=1 f11 f21 f12 P1 m1 x1P2=1 f22 P2 m2 x2P 同时作用时 1、P 2 同时作用时：x1 = f11P 1 + f12 P 2 x2 = f 21P 1 + f 22 P 2矩阵形式 X = FP 矩阵形式：⎡ x1 ⎤ X =⎢ ⎥ ⎣ x2 ⎦f ij 柔度影响系数f12 ⎤ f 22 ⎥ ⎦⎡ f11 F=⎢ ⎣ f 21⎡P 1⎤ P=⎢ ⎥ ⎣ P2 ⎦物理意义： 系统仅在第 j 个坐标受到 单位力作用时相应于第 i 个坐标上产生的位移柔度矩阵P1 m1 x1P2 m2 x2P1(t) m1 m2P2(t)&1 m1 & x&2 m2 & xX = FP⎡ x1 ⎤ ⎡ f11 ⎢x ⎥ = ⎢ f ⎣ 2 ⎦ ⎣ 21 f12 ⎤ ⎡ P 1⎤ ⎢P ⎥ f 22 ⎥ ⎦⎣ 2 ⎦当P 1、P 2 是动载荷时 集中质量上有惯性力存在⎡ x1 ⎤ ⎡ f11 ⎢x ⎥ = ⎢ f ⎣ 2 ⎦ ⎣ 21 f12 ⎤ ⎡ P && 1 (t ) − m1 x1 ⎤ ⎢ P (t ) − m & ⎥ f 22 ⎥ & x 2 2⎦ ⎦⎣ 2⎡ x1 ⎤ ⎡ f 11 ⎢x ⎥ = ⎢ f ⎣ 2 ⎦ ⎣ 21位移方程:f 12 ⎤⎛ ⎡ P1 (t ) ⎤ ⎡m1 ⎜⎢ −⎢ ⎥ ⎥ ⎜ f 22 ⎦⎝ ⎣ P2 (t ) ⎦ ⎣ 0&1 ⎤ ⎞ 0 ⎤⎡ & x ⎟ ⎥ ⎢ ⎥ &2 ⎦ ⎟ m2 ⎦ ⎣ & x ⎠&& ) X = F ( P − MXP1(t) m1 m2P2(t)⎡ x1 ⎤ X =⎢ ⎥ ⎣ x2 ⎦⎡P 1 (t ) ⎤ P=⎢ ⎥ P ( t ) ⎣ 2 ⎦&1 m1 & x&2 m2 & x位移方程 位移方程：&& ) X = F ( P − MX也可按作用力方程建立方程：&& + KX = P MX刚度矩阵&& + X = FP FMX柔度矩阵与刚度矩阵的关系 柔度矩阵与刚度矩阵的关系：&& KX = P − MX若K非奇异F=K−1FK = I&& ) X = K −1 ( P − MX应当注意：对于允许刚体运动产生的系统（即具有刚体自由度的系统） ， 柔度矩阵不存在。
机械系统动力学第四章 固有频率的实用计算方法
2 ( I Mu ) 0
2 I M 0
特征方程
对于二个自由度系统：
2 2 1 - m m 1 1 1 1 2 2 0 2 2 21m -22m 1 1 2
展开整理
1 m m 1 1 1 2 22 2 m m ( )0 12 1 1 2 2 1 2 2 1 4
U 带入公式 T m a x m a x 得：
T { u } K{ui } 2 i ni {ui }T M {ui }
4-2-7
利用4-2-7精确计算多自由度振动系统的固有频率，前 提条件是需要已知系统的振型，这是无法做到的。但 振动系统的一阶振型的近似值一般可以预测，大都数 情况下与其静载荷作用下产生的静变形十分接近。 例如例4-2-1所给出的振动问题，若取 u 1 1 1 代入式4-2-7进行试算：
2 2 J k a c l o
k a c l 0 即 J o
2 2
振动系统固有频率：
k a2 k a2 3 k a2 n 3 1 Jo m l 3 m l 3
第4章 固有频率的实用计算方法
4-1 单自由度系统
二．能量法 原理： 对于单自由度无阻尼自由振动系统，其响应为简谐振 d (T U ) 0 。在静平衡 动，系统 T Uc o n s t或 dt U 0 ,T T 位置，势能为0，动能达到最大，即： m a x。 在最大位移处，动能为0，势能达到最大， U U ,T 0 即： 。所以有： m a x
其点矩阵形式的动力方程为为第n段单元对转轴的转动惯量图434扭转振动单元状态向量表示gigd第4章固有频率的实用计算方法432传递矩阵法分析圆轴的扭转振动传递矩阵法的计算第n段单元的传递矩阵系统的传递矩阵的计算公式仍然可以表示为第4章固有频率的实用计算方法432传递矩阵法分析圆轴的扭转振动算法流程图图435a所示的一端固定一端自由的圆轴作扭转自由振动其中材料的切变模量为g密度为用传递矩阵法计算一阶固有频率
【结构动力学】第10章 多自由度体系2020
0
0
N
其中，ωn— 第n阶自振频率，{φ}n—第 n阶振型。
[Φ]和[Ω]也分别称为振型矩阵和谱矩阵。
13
5 DOF with uniform mass and stiffness
5 DOF Base Isolated 14
15
5 DOF with uniform mass and stiffness
k22 2m22 k2N 2m2n 0
k N1 2mN1 k N 2 2mN 2 k NN 2mNN
10
对于N个自由度的稳定结构体系，频率方程是关于ω2的 N次方程，
a N ( 2 ) N a N 1 ( 2 ) N 1 a1 2 a 0 0
由此可以解得N个正实根（ω12<ω22<ω32…<ωN2）。 ωn(n=1, 2, …, N)即为体系的自振频率。其中量值最小的 频率ω1叫基本频率(相应的周期T1=2π/ω1叫基本周期)。 从以上分析可知，多自由度体系只能按一些特定的频 率即按自振频率做自由振动。按某一自振频率振动时，结 构将保持一固定的形状，称为自振振型，或简称振型。
上述齐次方程组有非零解条件为：系数行列式为零
A [I ] 0
N×N矩阵[A]一般将有N个特征值，对应N个特征向量
6
§10-2 多自由度体系的自由振动
多自由度体系无阻尼自由振动的方程为：
M u K u 0
其中：[M]、[K]为N×N阶的质量和刚度矩阵 {u}和{ü}是N阶位移和加速度向量 {0}是N阶零向量
11
把相应的自振频率ωn代入运动方程的特征方程得到振型
K n 2 M n 0
{φ}n={φ1n, φ2n , …, φNn }T—体系的第n阶振型 。 ➢ 由于特征方程的齐次性(线性方程组是线性相关的)，振型向量 是不定的，只有人为给定向量中的某一值，例如令φ1n=1，才能确 定其余的值。 ➢ 实际求解时就是令振型向量中的某一分量取定值后才能求解。 虽然令不同的分量等于不同的量，得到的振型在量值上会不一样， 但其比例关系是不变的。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
作用力方程： MX KX P(t) 固有振动方程： MX KX 0
X Rn M、K Rnn
P(t) Rn
在考虑系统的固有振动时，首先考察系统的同步振动，即系
统在各个坐标上除了运动幅值不相同外，随时间变化的规律
都相同的运动。
假设系统的运动为： X φ f (t)
X Rn φ Rn
常数列向量 运动规律的时间函数 f (t) R1
f
(t)
at
b,
0
（1）正定系统
主振动
只可能出现形如 X φa sin( t ) 的同步运动
系统在各个坐标上都是按相同频率及初相位作简谐振动
（2）半正定系统
可能出现形如 X φa sin( t ) 的同步运动
也可能出现形如 X φ(at b) 的同步运动（不发生弹性变形 ）
MX KX 0
k2a2 k1a1 k1a12 k2a22
xD
D
PD
M
D
m
0
0 IC
xC
C
k1 k2a2
k2 k1a1
k2a2 k1a1 xC
k1a12
k2a22
C
PC M C
令： X D
xD
D
,
FD
PD
M
D
令： XC
xC
C
,
FC
PC
M
C
即：
m11
0
0 m22
x1 x2
k11
0
0 k22
x1 x2
P1 P2
若能够，则有：
m11x1 k11x1 P1
m22x2 k22x2 P2
方程解耦，变成了两个单自由度问题
使系统运动微分方程的全部耦合项全部解耦的坐标称为主坐标
多自由度系统振动 / 多自由度系统的动力学方程
D点和C点的坐标之间的关系：
xC xD eD
C D
写成矩阵形式： X D TXC
1 e
T 0
1
FD和FC 的关系
坐标变换矩阵
在C点加一对大小相等、方向相反的力 PD
得： PC PD M C ePD M D
C
D
C
D
xD xC
e
M D PD PD D C PD
写成矩阵形式： FC T T FD T 非奇异，因此： FD (T T )1 FC
FD (T T )1 FC
将 X D TXC 代入D点的方程，并左乘 T T ： T T MDTXC T T K DTXC T T FD FC
得： T T MDT MC
T T KDT KC
验证：
m MD me
me
IC
me2
m 0
MC
0
I
C
1 0 m me 1 e m 0
M C PC
DC
多自由度系统振动 / 多自由度系统的动力学方程
MD XD KD X D FD MC XC KC XC FC
X D [xD , D ]T XC [xC , C ]T
FD [PD , M D ]T FC [PC , MC ]T
X D TXC
T
1 0
e
1
FC T T FD
多自由度系统振动 / 多自由度系统的动力学方程
矩阵形式：
m me
IC
me me2
xD
D
k1 k2 k2a2 k1a1
k2a2 k1a1 xD
k1a12
k2a22
D
PD M D
存在惯性耦合
存在弹性耦合
多自由度系统振动 / 多自由度系统的动力学方程
问：能否找到这样一种坐标使得系统的运动微分方程既不出 现惯性耦合，也不出现弹性耦合？
X φf (t)
正定系统 0 半正定系统 0
多自由度系统振动 / 多自由度系统的自由振动
首先讨论正定系统的主振动
正定系统： MX KX 0
X Rn M 正定，K 正定
主振动： X φa sin( t ) 0 φ [1 2 n ]T
MD XD KD X D FD MC XC KC XC FC
多自由度系统振动 / 多自由度系统的动力学方程
MD XD KD X D FD MC XC KC XC FC
X D [xD , D ]T XC [xC , C ]T
FD [PD , M D ]T FC [PC , MC ]T
如果恰巧Y 是主坐标：
T T MT T T KT
对角阵
这样的T 是否存在？如何寻找？
多自由度系统振动 / 多自由度系统的自由振动
• 多自由度系统的自由振动
• 固有频率 • 模态 • 模态的正交性 • 主质量和主刚度 • 模态叠加法 • 模态截断法
多自由度系统振动 / 多自由度系统的自由振动
• 多自由度系统的固有频率
k2 k1a1
k2a2 k1a1 k1a12 k2a22
xC
C
PC
M
C
讨论：能对同一个系统选取两个不同的坐标，它们所描述的 运动微分方程之间有着怎样的联系？
多自由度系统振动 / 多自由度系统的动力学方程
m me
IC
me me2
xD
D
k1 k2 k2a2 k1a1
X [x1 x2 xn ]T φ [1 2 n ]T
多自由度系统振动 / 多自由度系统的自由振动
MX KX 0
X φf (t)
X Rn φ Rn
代入，并左乘 φT ： φT Mφf(t) φT Kφf (t) 0
：常数
f(t) f (t)
φT Kφ φT Mφ
2
M 正定，K 正定或半正定
l
选取D点的垂直位移及 角位移作为坐标
l1
A
DC
l2 B
选取质心C点的垂直位 移及角位移作为坐标
k1
a1
e
a2
k2
m me
IC
me me2
xD
D
k1 k2 k2a2 k1a1
k2a2 k1a1 xD
k1a12
k2a22
D
PD M D
m
0
0 IC
xC
C
k
k1 2a2
对于非零列向量 φ： φT Mφ 0 φT Kφ 0
令： 2 0
对于正定系统必有 0
对于半正定系统，有 0
多自由度系统振动 / 多自由度系统的自由振动
f(t) f (t)
φT φT
Kφ Mφ
2
f(t) 2 f (t) 0
Байду номын сангаас
a、b、 为常数
f (t) a sin(t ), 0
e 1me IC me2 0
1
0
I
C
多自由度系统振动 / 多自由度系统的动力学方程
结论：
假设对同一个系统所选择的两种不同的坐标X 和Y 有如下的变
换关系：
X TY
其中T 是非奇异矩阵，如果在坐标X下系统的运动微分方程为
：
MX KX P
那么在坐标Y 下的运动微分方程为：
T T MTY T T KTY T T P