Monday, December 10, 2012

与单自由度系统不同,n 自由度系统一般有n 个固有频率,因此可能出现n 次共振


自由度系统不同n 自由度

系统一
n 频率因此能出n 振。

70



4 章多自由度系统的振动

实际的物体与工程结构,其质量和弹性是连续分布的,系统具有无限多个自由度。为简化研究

和便于计算,可采用质量聚缩法或其它方法离散化,使系统简化为有限多个自由振动系统,或称为

多自由度振动系统。它的运动需要

n 个独立的坐标来描述。

4.1

变分原理与拉格朗日(Lagrange)运动方程

4.1.1
虚位移原理

在力学中遇到的第一个变分原理是虚位移原理。它是处理力学系统平衡问题的最基本原理,也

是分析力学的基础。

虚位移是指满足

固定时刻约束条件、假象、任意无限位移。可变系统,

虚位移

必须满足变形相容条件(连续条件)一个系统的虚位移个系统的广义坐标的变分。

假设

一个系统的广义坐标是q1 , q2 ,L, ql,其间存非定恒完整约束

f

k
(t; q1 , q2 ,L, ql ) = 0 (k =1, 2,L, m) (4-1)

若给

系统的位一虚位移dq那么根据定义,虚位移必定约束面上

f

k
(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 称为实位移。显然,它也是在

约束面上

的,

f

k
(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-6dq 满足条件式(4-4,也


,在这种情况下,系统的实位移与虚位移是不同的。


系统的约束条件定恒完整约束

PDF created with pdfFactory Pro trial version
www.pdffactory.com

71


f

k
(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独立的,平衡条件变为

Q

i
= 0 (i =1, 2,L, n) (4-9)

如果

将广义Q 广义坐标q 看作n 维空量,那么式(4-8)表Q q


向积,其义就:广义Q与虚位移dq 正交的。

4.1.2
达朗贝尔(DAlembert)原理

达朗贝尔提

了惯性力的概念虚位移原理的范围力学到动力学的领域达朗

贝尔

原理的普遍叙述方一个力学系统运动要在动力再加性力,它的任何

个位

置都以看作是平衡的位任何动力学问题按相力学问题来处理。

根据

原理,可以得出达朗贝尔原理的述方式:一个动力学系统的动力性力


任何虚位移上所作的虚

δ

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 (即t1t2时刻虚位移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)


TV 代入变分

δ δ
( )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



TP代入(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)



t1t2的虚位移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 1sinq = 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


mg

q
= kx (e)

将式

(e)代入(d)中,到运动方程

0

= + ÷ø

ö

çè

æ +

l q gq

k

mg


&& (f)


4-3 4-3体由簧支承图示

面内

运动。图示为平衡位质量为mIO


运动分方程。


:取体质O 离平衡位xy 转角q


广义坐标,

q


1 = 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

其中质量

矩阵

ú ú ú ú

û

ù

ê ê ê ê

ë

é

=


I

O

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



位移

列阵

q


T
= {x y q}

4.1.5
完整的非保守系统的拉格朗日运动方程

如果

系统是保守的,某种类型阻尼例如阻尼这时


瑞莱Rayleigh概念广义阻尼


因此于具有瑞莱保守系统,统一的标形式拉格朗


运动方程。


n 个自由度的系统,设广义坐标为q1 , q2 ,L, qn如果在系统j个质mj上作用有阻尼

F

Dj 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)



TP代入式(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

x
1 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-4a)所自由度系统的矩阵

首先令

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= 04-4c

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 = 0x3 = 14-4d,有

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-4a)所自由度系统的影响

首先

,在m1位力,m m 2 3 上不力,1 0 1 2 3 P = P = P = (4-5a)


4-4自由度系统

PDF created with pdfFactory Pro trial version
www.pdffactory.com

80



产生位移按定义应d d d 11 21 31 、、


m1P1,第一个弹的变



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-4a)所的系统,也可用影响立其运动分方程。


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



D
Mx

&& + 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-7afs24 -

No comments:

Post a Comment