http://course.tju.edu.cn/gchzhd/dianzijc/pdf/ch4.pdf
70
第 4 章多自由度系统的振动
实际的物体与工程结构,其质量和弹性是连续分布的,系统具有无限多个自由度。为简化研究
和便于计算,可采用质量聚缩法或其它方法离散化,使系统简化为有限多个自由振动系统,或称为
多自由度振动系统。它的运动需要n 个独立的坐标来描述。
4.1 变分原理与拉格朗日(Lagrange)运动方程
4.1.1 虚位移原理
在力学中遇到的第一个变分原理是虚位移原理。它是处理力学系统平衡问题的最基本原理,也
是分析力学的基础。
虚位移是指满足固定在某一时刻的约束条件的、假象的、任意的、无限小位移。对可变形系统,
虚位移必须满足变形相容条件(连续条件)。即一个系统的虚位移就是这个系统的广义坐标的变分。
假设一个系统的广义坐标是(q1 , q2 ,L, ql),其间存在非定恒的完整约束
fk (t; q1 , q2 ,L, ql ) = 0 (k =1, 2,L, m) (4-1)
若给系统的位形一虚位移dq,那么,根据定义,虚位移必定在约束面上,即
fk (t, q1 + δq1 , q2 + δq2 ,L, ql + δql ) = 0 (k =1, 2,L, m) (4-2)
将式(4-2)按泰勒级数展开,得到
δ δ δ ( 1, 2, , )
( , δ , δ , , δ ) ( ; , , , )
2
2
1
1
1 1 2 2 1 2
q k m
q
q
q
q
q
t q q q q q q t q q q
k k k
k k
L L
L L
+ =
¶
¶
+ +
¶
+
¶
+
+ + + =
l 高次项
l
l l l
f f f
f f
(4-3)
略去高次项后,得到虚位移应满足的条件为
δ δ δ 2 δ ( 1, 2, , )
2
1
1
q k m
q
q
q
q
q
k k k
k L = L
¶
¶
+ +
¶
¶
+
¶
¶
= l
l
f f f
f (4-4)
而系统的位形在dt时间内由q 运动到q +dq时,无限小的位移dq 称为实位移。显然,它也是在
约束面上的,即
fk (t + dt; q1 + dq1, q2 + dq2 ,L, ql + dql ) = 0 (k =1, 2,L, m) (4-5)
展开式(4-5),略去高次项后,得到实位移应满足的条件为
d d d 2 d d 0 ( 1, 2, , )
2
1
1
t k m
t
q
q
q
q
q
q
k k k k
k L = = L
¶
¶
+
¶
¶
+ +
¶
¶
+
¶
¶
=
f f f f
f l
l
(4-6)
从式(4-4)与式(4-6)中可以看出,满足式(4-6)的dq 不可能满足条件式(4-4),也就是
说,在这种情况下,系统的实位移与虚位移是不同的。
若系统的约束条件是定恒的完整约束,即
PDF created with pdfFactory Pro trial version www.pdffactory.com
71
fk (q1, , q2 ,L, ql ) = 0 (k =1, 2,L, m) (4-7)
则dfk 与ddfk 没有差别,真实的无限小位移属于虚位移。因此,对于自由质点系,以及只具有定恒
的完整的约束系统,真实的无限小位移可取作虚位移。
虚位移原理可表述为:力学系统在某一定位形时,平衡的必要与充分条件是:在此位形上所有
主动力在任何虚位移上所作的虚功之和等于零。即
δ δ 0
1
= =å=
i i
n
i
A Q q (4-8)
若广义坐标(q1 , q2 ,L, qn)是彼此独立的,则平衡条件变为
Qi = 0 (i =1, 2,L, n) (4-9)
如果将广义力Q 与广义坐标q 看作是n 维空间的矢量,那么,式(4-8)表达的就是Q 与q 的
无向积为零,其几何含义就是:广义力矢量Q与虚位移矢量dq 是正交的。
4.1.2 达朗贝尔(D’Alembert)原理
达朗贝尔提出了惯性力的概念,把虚位移原理的应用范围从静力学扩展到动力学的领域。达朗
贝尔原理的普遍叙述方式是:当一个力学系统运动时,只要在主动力上再加上惯性力,它的任何一
个位置都可以看作是平衡的位置。这样就可以把任何动力学问题按相当的静力学问题来处理。
根据虚功原理,可以得出达朗贝尔原理的另一种叙述方式:一个动力学系统的主动力及惯性力
在任何虚位移上所作的虚功之和等于零,即
δA + δAin = 0 (4-10)
其中dAin是惯性力所作的虚功。
当然,可以把虚功原理看作是达朗贝尔原理的一个特例。这样,达朗贝尔原理就是力学的最基
本的变分原理。
4.1.3 哈密顿原理(Hamliton)原理
从 18 世纪开始,很多数学家就致力于寻求一个能推导出牛顿力学定律的统一的力学原理。直
到19 世纪,哈密顿建立了最小作用原理- 哈密顿原理,科学家们的梦想才得以实现。
哈密顿原理表明:一个具有完整约束的力学系统的运动必定使积分作用量
I T A t
t
t
( )d 2
1
= ò + (4-11)
取驻值,即
δ δ ( )d 0 2
1
I = ò T + A t =
t
t
(4-12)
其中,I称为哈密顿作用量,T是系统的动能,A 是主动力所作的功。如果主动力有势,那么式(4-12)
可以写成
PDF created with pdfFactory Pro trial version www.pdffactory.com
72
δ δ ( )d 0 2
1
I = ò T -V t = t
t
(4-13)
式(4-12)与式(4-13)可以解释为:完整的力学系统从状态“1”到状态“2”的各种可能运动中,
唯有真实运动使哈密顿作用量取驻值。
为了证明式(4-12),哈密顿将式(4-10)变换为
δ d δ d 0 2
1
2
1
ò A t + ò Ain t =
t
t
t
t
(4-14)
其中
R t
t
m R R t m R R m R
t
A t j j j
N
j
t
t
t
t
j j j
N
j
j j j
N
j
t
in t
t
t
δ d
d
( )δ d δ d
d
δ d d
1 1 1
2
1
2
1
2
1
2
1
ò ò å & å & ò å &
= = =
+ úû
ù
êë
é
= - = -
只要取δ δ 0
2 = 1 = 1 = Rj t t R t=t (即t1与t2时刻虚位移dRj为零),则有
R t m R R t T t
t
A t m R t
j j j t
N
j
t
j j j t
N
j
t
in t
t
t
δ d δ d δ d
d
δ d d 2
1
2
1
2
1
2
1 1 1
ò = ò å = ò å = ò
= =
& & & (4-15)
其中, 2
1 2
1
j j
N
j
R m T & å=
= 是系统的动能。
将式(4-15)代入式(4-14)中,得到
δ ( )d 0 2
1
ò T + A t =
t
t
(4-16)
这就证明了泛函驻值形式的变分原理——哈密顿原理。
若系统的主动力有势,则式(4-16)可写成
δ ( )d δ d 0 2
1
2
1
ò T -V t = ò L t = t
t
t
t
(4-17 )
其中,L = T -V 称为拉格朗日函数。
若系统的主动力一部分有势,而另一部分没有势,则式(4-16)可写成
δ d δ ( )d 0 2
1
2
1 1
= - + ÷ø
ö
çè
æ ò å ò
=
Q q t T V t t
i i t
N
k
t
t
(4-18)
其中,Qi ( i = 1, 2, …, n)是与没有势的那些主动力有关的广义力。
上述哈密顿原理是对离散系统导出的,只要将连续系统的动能T 与势能V 代入式(4-17),它
对连续系统照样适用。
例 4-1 图4-1所示系统中,半径为r的均匀圆盘在槽内作不滑动的滚动。已知圆盘质量为m,
槽的半径为R。试用哈密顿原理建立系统的运动方程。
解:
若选择q 为广义坐标,则系统微幅振动时的能量为
2 2
2
[( ) ]
2
1 q& j&
T m R r I A
1
= - + (a)
其中,j&为圆盘的角速度,IA = mr2/2 是圆盘对质心的转动惯量。圆盘作不滑动的滚动时,存在有
PDF created with pdfFactory Pro trial version www.pdffactory.com
73
j&r =q&(R - r) (b)
由此,得到
j& q&
r
R - r
= (c)
将式(c)代入式(a),得到
2
2
2
4
3 q& ÷ø
ö
çè
æ -
=
r
T mr R r
(d)
而系统的位能
( ) 2
2
V = mg(R - r)(1- cosq ) » 1 mg R - r q (e)
将T与V 代入变分式
δ δ ( )d 0 2
1
I = ò T -V t = t
t
中,得到
( ) δ d 0
) δ d
2
( ) δ - 3
2
3
δ ( ) δ d
2
3
( ) d
2
1
4
δ 3
2
1
2
1
2
1
2
1
2
1
2 2
2
2
2 2
2
2
- - =
= - -
ú úû
ù
ê êë
é
- - ÷ø
ö
çè
æ -
=
ú ú
û
ù
ê ê
ë
é
- - ÷ø
ö
çè
æ -
ò
ò
ò
ò
mg R r t
m R r m(R r t
mg R r t
r
mr R r
mg R r t
r
mr R r
t
t
t
t
t
t
t
t
t
t
q q
q q q q
q q q q
q q
& &
& &
&
(f)
由于,t = t1 = t2时,哈密顿原理要求dq = 0,所以,式(f)满足时,必有
( ) ( ) 0
2
3 m R - r 2q&& + mg R - r q = (g)
式(g)就是系统微幅振动时的运动方程。
4.1.4 完整的保守系统的拉格朗日运动方程
n 个自由度的系统,在一般情况下,动能可能是时间t、广义坐标qi以及广义速度q&i的函数。即
T = T(t; q1 , q2 ,L, ql ; q&1 , q& 2 ,L, q&n )
而势函数只是广义坐标qi的函数,即
( , , , ) 1 2 n V = V q q L q
图4-1圆盘微幅振动
PDF created with pdfFactory Pro trial version www.pdffactory.com
74
将 T与P代入式(4-18)中,进行变分运算,得到:
δ d 0
d
d
( δ )d δ
( δ )d δ δ δ d
1
1 1
1 1
2
1
2
1
2
1
2
1
2
1
= úû
ù
êë
é
÷ ÷ø
ö
ç çè
æ
¶
¶
+
¶
¶
-
¶
¶
-
÷ ÷ø
ö
ç çè
æ
¶
¶
= +
÷ ÷ø
ö
ç çè
æ
¶
¶
-
¶
¶
+
¶
¶
+
ò å
ò å å
ò å ò å
=
= =
= =
q t
q
V
q
T
q
T
t
q
q
Q q t T
q t
q
q V
q
q T
q
Q q t T
i
i i i
n
i
t
t
t
t
i
i
n
i
i i
n
i
t
t
i
i
i
i
i
i
n
i
t
i i t
n
i
t
t
&
&
&
&
(4-19)
由于,t = t1 = t2时,哈密顿原理要求dqi = 0,因此,式(4-19)变成
δ d 0
d
( δ )d d
1 1
2
1
2
1
= úû
ù
êë
é
÷ ÷ø
ö
ç çè
æ
¶
¶
+
¶
¶
-
¶
¶
ò å - ò å
= =
q t
q
V
q
T
q
T
t
Q q t i
i i i
n
i
t
i i t
n
i
t
t &
(4-20)
在t1与t2区间的虚位移dqi是任意的,而且dqi彼此独立的。因此,由式(20)得到著名的拉格朗日
方程
( 1, 2, , )
d
d Q i n
q
V
q
T
q
T
t i
i i i
L
&
= =
¶
¶
+
¶
¶
- ÷
÷ø
ö
ç çè
æ
¶
¶
(4-21)
拉格朗日方程提供了解决有限自由度完整系统运动的一个普遍的简单而又统一的方法。
例 4-2 图4-2所示系统,摆的支点在水平方向受到弹性约束,其总刚度为k,摆的质量为m,
摆长为l。试用拉格朗日方程求出系统的运动方程。
解:
(1)选择x 及q 为广义坐标。
(2)动能及势能
动能: 2 [( sin ) ]2
2
[ ( cos ) ] 1
2
T = 1 m x& + l q q& + m l q q& (a)
势能: (1 cos )
2
V = 1 kx2 + mgl - q (b)
(3)广义外力为零
本例题中广义外力。
(4)运动方程
将式(a)与式(b)代入式(4-21),得到
cos sin sin 0
cos sin 0
2
2 2
+ - + =
+ - + =
q q q q q
qq qq
ml x ml mlx mg
mx ml ml kx
&& && & &
&& && &
(c)
这就是摆的运动方程。当微幅振动时,可取cosq ≈1,sinq = 0,并可略去高阶项,则式(c)可简化为
0
0
+ + =
+ + =
q q
q
mx ml mg
mx ml kx
&& &&
&& &&
(d)
两式相减得到
图4-2摆振系统
PDF created with pdfFactory Pro trial version www.pdffactory.com
75
mgq = kx (e)
将式(e)代入式(d)中,得到运动方程
0 = + ÷ø
ö
çè
æ + l q gq
k
mg && (f)
例4-3 图4-3所示的刚体由四根拉伸弹簧支承,被限制在图示平
面内运动。图示位置为平衡位置。且质量为m,转动惯量IO。试导出微
幅运动微分方程。
解:取刚体质心O 点偏离平衡位置的x、y 和刚体绕质心的转角q
为广义坐标,即
q1 = x, q2 = y, q3 =q
并且四根弹簧端点的坐标分别为
, 0, 0
, ,
4 4 1 2 3 4
1 1 2 2 3 3
= - = = = =
= + = - = +
y y a y y x x
x x a x x a y y a
q
q q q
系统的动能为
2 2 2
2
( ) 1
2
1 & & q& T = m x + y + IO
系统的势能为
2
4 4
2
3 3
2
2 2
2
1 1 ( )
2
( ) 1
2
( ) 1
2
( ) 1
2
V = 1 k x + a q + k x - a q + k y + a q + k y - a q
计算拉格朗日方程中各项导数如下
3 3 3 4 4 4
1 1 1 2 2 2
3 3 4 4
1 1 2 2
( ) ( )
( ) ( )
; 0
d
d
( ) ( )
; 0
d
d
( ) ( )
; 0
d
d
k y a a k y a a
V k x a a k x a a
T I T
t
k y a k y a
y
V
y
my T
y
T
t
k x a k x a
x
V
x
mx T
x
T
t
O
q q
q q
q
q
q
q
q q
q q
+ + - -
= + - -
¶
¶
=
¶
¶
= ÷ø
ö
çè
æ
¶
¶
= + + -
¶
¶
=
¶
¶
= ÷
÷ø
ö
ç çè
æ
¶
¶
= + + -
¶
¶
=
¶
¶
= ÷ø
ö
çè
æ
¶
¶
&&
&
&&
&
&&
&
代入拉格朗日方程,得系统运动微分方程为
图4-3刚体微幅运动
PDF created with pdfFactory Pro trial version www.pdffactory.com
76
( ) 0
( ) ( )
( ) ( ) 0
( ) ( ) 0
2
4 4
2
3 3
2
2 2
2
1 1
1 1 2 2 3 3 4 4
3 4 3 3 4 4
1 2 1 1 2 2
+ + + + =
+ - + -
+ + + - =
+ + + - =
q
q
q
q
k a k a k a k a
I k a k a x k a k a y
my k k y k a k a
mx k k x k a k a
O
&&
&&
&&
引入记号,写为矩阵形式
mq&& + kq = 0
其中质量矩阵为
ú ú ú ú
û
ù
ê ê ê ê
ë
é
=
IO
m
m
0 0
0 0
0 0
m
刚度矩阵为
ú ú ú ú
û
ù
ê ê ê ê
ë
é
- - + + +
+ -
+ -
=
2
4 4
2
3 3
2
2 2
2
1 1 2 2 3 3 4 4 1 1
3 4 3 3 4 4
1 2 1 1 2 2
0
0
k a k a k a k a k a k a k a k a
k k k a k a
k k k a k a
k
位移列阵为
qT = {x y q}
4.1.5 完整的非保守系统的拉格朗日运动方程
如果系统是非保守的,且存在某种类型的耗能力(阻尼力),例如存在粘性阻尼。这时引入所
谓的瑞莱(Rayleigh)散逸函数的概念后,象有势(位)的广义力能够从位能导出一样,阻尼力能够
从散逸函数导出。因此,对于具有瑞莱散逸函数的非保守系统,仍可求得统一的标准形式的拉格朗
日运动方程。
对于n 个自由度的系统,设广义坐标为q1 , q2 ,L, qn。如果在系统上第j个质点mj上作用有阻尼
FDj c j Rj = - & (4-22)
其中,Rj是质点mj的位置向径。
如果质点mj有虚位移dRj时,阻尼力FDj所作到虚功为
i
i
j
j j
n
i
i
i
j
Dj
n
i
Dj Dj j q
q
R
q c R
q
R
δA F δR F δ δ
1 1 ¶
¶
= -
¶
¶
= =å å
= =
& (4-23)
由于i
i
j
n
i
i
i
j
n
i
j q
q
R
t
q
q
R
R& &
¶
¶
=
¶
¶
=å å
=1 =1 d
d
(4-24)
PDF created with pdfFactory Pro trial version www.pdffactory.com
77
i
j
i
j
q
R
q
R
¶
¶
=
¶
¶
&
&
(4-25)
所以,式(4-23)可以写成
j j i
i
n
i
i
i
j
j j
n
i
Dj q c R q
q
A c R δ
2
δ δ 1 2
1 1
÷ø
ö
çè
æ
¶
¶
= -
¶
¶
= -å å
= =
&
& &
&
&
q q
(4-26)
对于具有N个质点的系统,阻尼力所作的虚功是
i
i
n
i
j j i
i
n
i
N
j
Dj
N
j
D q A A c R δq D δ
2
δ δ 1
1
2
1 1 1 q q&
&
& ¶
¶
- = ÷ø
ö
çè
æ
¶
¶
= -å = -å å å
= = = =
(4-27)
其中
2
1 2
1
j j
N
j
c D R& å=
= & & & & q& TBq&
1 1 1 1 1 2
1
2
1
2
1 = = ÷
÷ø
ö
ç çè
æ
¶
¶
¶
¶
= å å å å å
= = = = =
kl k l
n
l
n
k
k l
l
j
k
j
i
n
l
n
k
N
j
q q b q q
q
R
q
R
c (4-28)
称为瑞莱散逸函数。其中
÷ ÷
ø
ö
ç ç
è
æ
¶
¶
¶
¶
=å=
k l
j
j
N
j
kl q
R
q
R
b c
1
,
÷ ÷ ÷ ÷ ÷ ÷ ÷ ÷
ø
ö
ç ç ç ç ç ç ç ç
è
æ
=
n n nn
n
n
b b b
b b b
b b b
L
L L L L
L
L
1 2
21 22 2
11 12 1
B
是对称正定方程。
具有散逸函数的非保守系统的虚功表达式是:
δAin + δA + δAD = 0
对于具有完整约束的非保守的系统,哈密顿变分式变为
δ ( δ )d δ ( )d δ d 0 2
1
2
1
2
1 1
= ò å + ò - + ò =
=
I Q q t T V t A t D
t
t
t
i i t
n
i
t
t
(4-29)
将T与P代入式(4-30)中,进行变分运算,并利用式(4-27),则可得到
δ d δ d 0
d
d
1 1
2
1
2
1
=
¶
¶
- úû
ù
êë
é
¶
¶
-
¶
¶
+ ÷ ÷ø
ö
ç çè
æ
¶
¶
ò å - ò å
= =
q t D q t
q
V
q
T T
t
Q i
i
n
i
t
i t
i i i
i
n
i
t
t q& q&
或
δ d 0
d
d
1
2
1
= úû
ù
êë
é
¶
¶
-
¶
¶
-
¶
¶
+ ÷
÷ø
ö
ç çè
æ
¶
¶
ò å -
=
q t
q
D
q
V
q
T T
t
Q i
i i i i
i
n
i
t
t q&
由于虚位移δqi (i =1, 2,L, n)是任意的,且互相是独立的。因此,有
( 1, 2, , )
d
d Q i n
q
D
q
V
q
T
q
T
t i
i i i i
L
& &
= =
¶
¶
+
¶
¶
+
¶
¶
- ÷
÷ø
ö
ç çè
æ
¶
¶
(4-30)
这就是具有瑞莱散逸函数的非保守系统的拉格朗日运动方程的标准形式。
PDF created with pdfFactory Pro trial version www.pdffactory.com
78
4.2 多自由度系统自由振动的运动微分方程
一般情况下,n 个自由度无阻尼系统的自由振动的运动微分方程具有以下形式
m x m x m x k x k x k x
m x m x m x k x k x k x
m x m x m x k x k x k x
n n n n
n n n n
n n nn n n n nn n
11 1 12 2 1 11 1 12 2 1
21 1 22 2 2 21 1 22 2 2
1 1 2 2 1 1 2 2
0
0
0
&& && &&
&& && &&
&& && &&
+ + × + + + + + =
+ + + + + + + =
+ + + + + + + =
ì
í
ïï
î
ïï
L L
L L
LL
L L
(4-31)
若用矩阵表示,则可写成
Mx&& + Kx = 0 (4-32)
式中
( ) ( )T
n
T
n x1 x2 L x && &x&1 &x&2 L &x& x = ,x =
分别是系统的坐标矢量和加速度矢量。
M =
é
ë
êêêêê
ù
û
úúúúú
m m m
m m m
m m m
n
n
n n n n
11 12 1
21 22 2
1 2
L
L
L L L L
L
, K =
é
ë
êêêêêê
ù
û
úúúúúú
k k k
k k k
k k k
n
n
n n n n
11 12 1
21 22 2
1 2
L
L
L L L L
L
分别称为系统的质量矩阵和刚度矩阵。
4.2.1刚度影响系数作用力方程
方程(4-31)中各项均为力的量纲,因此,称之为作用力方程。其中刚度矩阵中的元素称刚度
影响系数(在单自由度系统中,简称弹性常数)。它表示系统单位变形所需的作用力。具体地说,如果
使第j 个质量沿其坐标方向产生单位位移,沿其它质量的坐标方向施加作用力而使它们保持不动,
则沿第i 个质量坐标方向施加的力,定义为刚度影响系数kij;在第j 个质量坐标方向上施加的力称
刚度影响系数k jj 。由刚度影响系数的物理意义,可直接写出刚度矩阵,从而建立作用力方程,这种
方法称为影响系数法。
现分析求出图4-4(a)所示的三自由度系统的刚度矩阵。
首先令m1有单位位移x1 = 1、而x x 2 3 、保持不动,即x x 2 3 = = 0。在此条件下系统保持平
衡,按定义需加于三物块的力为k k k 11 21 31 、、(图中箭杆上画有斜线的力)。画出各物块的受力图4-4
(b)根据平衡条件,有
0 11 1 2 21 2 31 k = k + k ,k = -k ,k =
同理,令x x x 1 2 3 = 0, = 1, = 0,画出受力图4-4(c),则有
k k k k k k k 12 2 22 2 3 32 3 = - , = + , = -
PDF created with pdfFactory Pro trial version www.pdffactory.com
79
最后令x1 = x2 = 0,x3 = 1,画出受力图4-4(d),有
k k k k k 13 23 3 33 3 = 0, = - , =
因此刚度矩阵为
K =
+ -
- + -
-
é
ë
êêê
ù
û
úúú
k k k
k k k k
k k
1 2 2
2 1 3 3
3 3
0
0
上式表明,k k ij ji = 。因此,刚度矩阵一般是对称的。实际上任何多自由度线性系统都具有这个性
质。即
K = KT (4-33)
4.2.2柔度影响系数位移方程
在单自由度的弹簧——质量系统中,若弹簧常数是k,则
1
k
就是物块上作用单位力时弹簧的变
形,也称柔度影响,用d 表示。n 自由度系统的柔度矩阵D 为n 阶方阵,其元素d ij称为柔度影响系
数,表示单位力产生的位移。具体地说,仅在第j 个质量的坐标方向上受到单位力作用时相应于在
第i个质量的坐标方向上产生的位移,即定义为d ij 。
现分析求出图4-4(a)所示的三自由度系统的柔度影响系数。
首先,在m1上施加单位力,而m m 2 3 、上不加力,即令1 0 1 2 3 P = ,P = P = (图4-5(a)),这
图4-4三自由度系统受力图
PDF created with pdfFactory Pro trial version www.pdffactory.com
80
时三物块所产生的静位移按定义应分别是d d d 11 21 31 、、。
当m1受到P1作用后,第一个弹簧的变形
为
1
1 k
,第二和第三个弹簧的变形为零。
所以三物块的位移都是
1
1 k
,即
d d d 11
1
21
1
31
1
= 1 = 1 = 1
k k k
, ,
同理,如令P P P 2 1 3 = 1, = = 0。这时
第一和第二弹簧均受单位拉力,其变形分别为
1 1
1 2 k k
, ,而第三个弹簧不受力,故其变形为零。因
此有
d d d 12
1
22
1 2
32
1 2
= 1 = 1 + 1 = 1 + 1
k k k k k
, ,
再令P P P 1 2 3 = = 0 , = 1,可得到
d d d 13
1
23
1 2
33
1 2 3
= 1 = 1 + 1 = 1 + 1 + 1
k k k k k k
, ,
因此,该系统的柔度矩阵为
D =
é
ë
êêê
ù
û
úúú
= + +
+ + +
é
ë
êêêêêêê
ù
û
úúúúúúú
d d d
d d d
d d d
11 12 13
21 22 23
31 32 33
1 1 1
1 1 2 1 2
1 1 2 1 2 3
1 1 1
1 1 1 1 1
1 1 1 1 1 1
k k k
k k k k k
k k k k k k
上式表明,d d i j ji = 。因此,柔度矩阵一般也是对称的。实际上任何多自由度线性系统都具有这个
性质。即
D = DT (4-34)
对于图4-4(a)所示的系统,也可用柔度影响系数来建立其运动微分方程。
如设x x x 1 2 3 , , 分别表示质量m m m 1 2 3 , , 的位移。系统运动时,质量m m m 1 2 3 , , 的惯性力使
弹簧产生变形,对于线性弹性体应用叠加原理可得到
图4-5三自由度系统柔度影响系数
PDF created with pdfFactory Pro trial version www.pdffactory.com
81
x m x m x m x
x m x m x m x
x m x m x m x
1 1 1 11 2 2 12 3 3 13
2 1 1 21 2 2 22 3 3 23
3 1 1 31 2 2 32 3 3 33
= - + - + -
= - + - + -
= - + - + -
( && ) ( && ) ( && )
( && ) ( && ) ( && )
( && ) ( && ) ( && )
d d d
d d d
d d d
写成矩阵形式
x
x
x
m
m
m
x
x
x
1
2
3
11 12 13
21 22 23
31 32 33
1
2
3
1
2
3
0 0
0 0
0 0
æ
è
ççç
ö
ø
÷÷÷
=
é
ë
êêê
ù
û
úúú
é
ë
êêê
ù
û
úúú
-
-
-
æ
è
ççç
ö
ø
÷÷÷
d d d
d d d
d d d
&&
&&
&&
或
x = -DMx&& (4-35a)
和
DMx&& + x = 0 (4-35b)
上式对任一个n 自由度系统都是正确的。式(4-35)称位移方程。它是振动微分方程的另一种形式。
为了比较作用力方程与位移方程,将式(4-32)改写为
Kx = -Mx&&
如果K 是非奇异的,即K 的逆矩阵K -1存在,对上式两端左乘K -1,得
x = K -1 (-Mx&&) (4-36)
比较式(4-36)与式(4-35)得出
D = K -1 (4-37)
上式即柔度矩阵与刚度矩阵之间的关系。即当刚
度矩阵是非奇异时,刚度矩阵K 与柔度矩阵D
互为逆矩阵;当刚度矩阵是奇异时,不存在逆矩
阵即无柔度矩阵,此时系统的平衡位置有无限多
或者说它有刚体运动。例如图4-6所示的系统,
该系统具有刚体运动,因此柔度矩阵不存在。
例 4-4 试写出图4-7所示刚体AB 的刚度矩阵并建立系统的运动微分方程。
解:刚体AB 在图面内的位置可以由其质心C的坐标yC (以水平位置O 为坐标原点,且水平运
动不计)和绕C的转角q 确定。
图 4-7(a)
No comments:
Post a Comment