多自由度系统近似计算方法

在线性多自由度系统振动中，振动问题归结为刚度矩阵和质量矩阵的广义特征值问题，缺点：当系统自由度较大时，求解计算工作量非常大。

本章介绍邓克利法，瑞利法，里茨法，传递矩阵法等计算方法，可作为实用的工程计算方法对系统的振动特性作近似计算。

1、邓克利法由邓克利（Dunkerley ）在实验确定多圆盘的横向振动固有频率时提出的，便于作为系统基频的计算公式 。

自由振动作用力方程：0KX XM =+ n R ∈X 左乘柔度矩阵F = K -1，位移方程：0X X FM =+ 定义D=FM 为系统的动力矩阵：0X XD =+ 作用力方程的特征值问题：φφM K 2ω= 位移方程的特征值问题：φφλ=D 特征值：22221n ωωω<<< ，n λλλ>>> 21 关系：2/1i i ωλ=位移方程的最大特征根：211/1ωλ=，对应着系统的第一阶固有频率。

位移方程的特征方程：0=-I D λ展开：0)()1(1111=++++---n n n nna a a λλλD tr d d d a nn -=+++-=)(22111例：022211211=--λλd d d d0)]()([)1(21122211221122=-++--d d d d d d λλ当 M 为对角阵时：)(FM D tr tr =∑==ni iii m f 1特征方程又可写为：0)())((21=---n λλλλλλ有：∑=-=ni i a 11λtrD -=∑=-=ni i ii m f 1∑∑===ni iii ni im f 11λ∑∑===ni i ii ni im f 1121ω如果只保留第 i 个质量，所得的单自由度系统的固有频率为：iii ii i m f m k 12==ω例：两自由度系统柔度矩阵：⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡+=2111111111k k k k kF （1）只保留 m 1 时1111k f =，1121m k =ω（2）只保留 m 2 时122122111k k k f =+=，21222m k =ω将2i ω代入：22221121111nni iωωωω+++=∑=对于梁结构系统，第二阶及第二阶以上的固有频率通常远大于基频，因此左端可只保留基频项，有：22221211111nωωωω+++≈例：三自由度系统⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡----+⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡0002223101220010*********x x x k x x x m 采用常规方法，固有频率：m k /3730.01=ω，m k /3213.12=ω，m k /0286.23=ω邓克利法当 m1 单独存在时：m k /21=ω 当 m2 单独存在时：k k k k k k 21212112=+=，m k /1222=ω当 m3 单独存在时：kk k k k 251111321123=++=，52123k k =，mk 523=ω代入邓克利法公式：22221211111nωωωω+++≈，mk /3535.01=ω2、瑞利法瑞利法是基于能量原理的一种近似方法，可用于计算系统的基频，算出的近似值为实际基频的上限，配合邓克利法算出的基频下限，可以估计实际基频的大致范围。

n 自由度保守系统：0KX XM =+ n R ∈X 主振动 ：)sin(ϕω+=t φX 动能与势能：XM X TT 21=，KX X TV 21=最大值：φφM TT 2max 21ω=，φφK TV 21max =maxmax V T =得2)(ω==φφφφφM K TTR 成为瑞利商。

对于第 i 阶模态：2)()()()()()(i i Ti i Ti i R ω==φφφφφM K当φ为一般向量时（不是实际模态），总能展开为 n 个正则模态的线性组合：)()2(2)1(1n NN NNa a a φφφ＋ ++==ϕφ∑==nj j Nja 1)(φaΦN =其中],,,[)()2()1(n N N N N φφφ =Φ，T n a a a ],,,[21 =a 代入瑞利商：aM ΦΦa a K ΦΦa N TN TN TN TR =)(ϕIaa Λa a TT =∑∑===nj jn j jjaa12122ω可以证明：21ω和 2n ω分别为瑞利商的极小值和极大值，即：221)(n R ωϕω≤≤ 分析：若将瑞利商右端分子内的所有j ω换为1ω，由于1ω是最低阶固有频率，因此：21121212)(ωωϕ=≥∑∑==nj j nj jaaR由瑞利商公式知，当)1(φ=ϕ确为第一阶模态时，有：21)(ωϕ=R 。

因此，瑞利商的极小值为21ω，同理可证明，瑞利商的极大值为2n ω 如果ϕ接近第 k 阶真实模态)(k φ=ϕ，比起 a k ，其它系数很小kj n j a a k j j ≠==,,,2,1 ，ε ，1<<j ε代入，得：∑=-+≈nj j k jkR 12222)()(εωωωϕ 例如 k ＝1，n j a a j j ,,21 ==，ε222212222222121)(nnn a a a a a a R ++++++=ωωωϕ21221222122122221222121a a a a a a n nn εεωεωεω++++++=约去a 1()222222222211nnn R εεωεωεωϕ++++++=分子上加减1项())1/()()(222212212321222122123212222222221n n n n n R εεωεωεωεωεωεωεωεωεωϕ+++⎪⎪⎭⎫⎝⎛+++-+++++++= 22222122222211)()1(nni i inεεωωεεεω+++-++++=∑= 22222122211)(nni i i εεωωεω+++-+=∑=∑=-+≈ni i i 2212221)(ωωεω因此，若ϕ与)(k φ的差异为一阶小量，则瑞利商与2k ω的差别为二阶小量。

对于基频的特殊情况，令k ＝1，则由于)~2(0212n j j=>-ωω瑞利商在基频处取极小值。

利用瑞利商估计系统的基频所得的结果必为实际基频的上限，ϕ愈接近系统的真实模态，算出的固有频率愈准确。

例题：接上例mk /3730.01=ω，采用邓克利法，基频：m k /3535.01=ω取在2m 质量上施加力P 所产生的“静变形曲线”作为近似的第一阶主振型，即：T]5.2,2,1[=ϕ代入瑞利商公式：mk R 142857.0)(=ϕ，mk 3780.01=ω与精确值相比，相对误差1.34%。

3、里茨法里兹法是瑞利法的改进，用里兹法不仅可以计算系统的基频，还可以算出系统的前几阶频率和模态，瑞利法算出的基频的精度取决于假设的振型对第一阶主振型的近似程度，而且得到的基频总是精确值的上限，里兹法将对近似振型给出更合理的假设，从而使算出的基频值进一步下降。

里兹法基于与瑞利法相同的原理，但将瑞利使用的单个假设模态改进为若干个独立的假设模态的线性组合：ΠA ηηηη==+⋯++=∑=rj j jr r aa a a 1)()()2(2)1(1ϕ，1)(⨯∈n i R ηrn r R⨯∈=],,,[)()2()1(ηηηΠ ，121],,,[⨯∈=r T r R a a a A代入瑞利商：)()(ΠA R R =ϕAM ΠΠA A K ΠΠA TTTT =AM A A K A T T=2ω=r r TR⨯∈=K ΠΠK ，r r T R ⨯∈=M ΠΠM由于()ϕR 在系统中的真实主振型处取驻值，所以 A 的各个元素应当从下式确定：)2,1(,0r j a R j⋯⋯==∂∂。

代入：),,2,1(,0)()(2r j a a TjTj⋯==∂∂-∂∂A M A A K A ω),,2,1(,2)(2)()()(r j a a a a TjTjjTTjTj ⋯==∂∂=∂∂+∂∂=∂∂A K e A K A A K A A K A A K A其中j e 是 r 阶单位矩阵的第 j 列。

上面 r 个方程可合成为：AK A K A A2)(=∂∂T表示将函数分别对 A 的各个元素依次求偏导，然后排列成列向量A∂∂。

同理，有：AM A M A A2)(=∂∂T两项代入得：0A M K =-)(2ω由于M K 、的阶数 r 一般远小于系统自由度数 n ，上式所示的矩阵特征值问题比原来系统的矩阵特征值问题解起来容易得多，因此里兹法实际上是一种缩减系统自由度求解固有振动的近似方法，M K 、就是自由度缩减为 r 的新系统的刚度矩阵和质量矩阵。

可求出 r 个特征根22221,,r ωωω⋯，，及相应的特征向量)()2()1(,,,r A A A ⋯，原来系统的前 r 阶固有频率可近似取为：),,2,1(,22r j j j ⋯==ωω，相应的前 r阶主振型近似取为：),,2,1(,)()(r j j j ⋯==ΠA ϕ。

正交性分析ji j Ti j T Ti j Ti ≠===当 0)()()()()()(AM AAM ΠΠAM ϕϕ即：00)()()()(==j T i j T i A K A A M A ，因此：得出的近似主振型式关于矩阵 M 和 K 相互正交。

例题：接前例采用常规方法，固有频率：mk /3730.01=ω，m k /3213.12=ω，m k /0286.23=ω采用邓克利法，基频：m k /3535.01=ω，采用瑞利法，基频：m k /3780.01=ω。

将假设的振型取为⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡132211][21－＝＝）（）（ηηΠ 缩减后的新系统的刚度矩阵和质量矩阵：⎥⎦⎤⎢⎣⎡--=k kk k 20444TK ΠΠΚ＝，⎥⎦⎤⎢⎣⎡--==m mm m TT 723M ΠΠM 特征值问题：⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡-+-+--0072044234αααα，km 2ωα=139853.01=α，860147.22=α，[]T1,927547.4)1(=A，[]T 1,018449.0)2(-=A固有频率：mk mk 373969.0111===αωω，mk mk 691197.1222===αωω主振型：[]T1860147.0430073.01)1()1(βϕ==ΠA[]T1860148.1930074.02)2()2(--==βϕΠA主振型精确值：[]T18608.04626.0)1(=φ，[]T17458.09339.2)2(--=φ里兹法得到的基频精度比用瑞利法的高，但第二阶固有频率的精度还欠佳。