机械振动大作业——简支梁的各情况分析

机械振动大作业
姓名：徐强
学号：SX1302106
专业：航空宇航推进理论与工程
能源与动力学院
2013年12月
简支梁的振动特性分析
题目:针对简支梁、分别用单、双、三、十个自由度以及连续体模型，计算其固有频率、固有振型。

单、双、三自由度模型要求理论解；十自由度模型要求使用李兹法、霍尔茨法、矩阵迭代法、雅可比法、子空间迭代法求解基频；连续体要求推导理论解，并通过有限元软件进行数值计算。

解答：
一、 单自由度简支梁的振动特性
如图1，正方形截面（取5mm ×5mm ）的简支梁，跨长为l =1m ，质量m 沿杆长均匀分布，将其简化为单自由度模型，忽略阻尼，则运动微分方程为0=+•
•kx x m ,固有频率ωn =
eq
eq m k ，其中k 为等效刚度，
eq m 为等效质量。

因此，求出上述两项即可知单自由度简支梁的固有
频率。

根据材料力学的结果，由于横向载荷F 作用在简支梁中间位置而
引起的变形为）（2
24348EI F -)(x l x x y -=（2
0l x ≤≤）, 48EI F -3max l y =为最大挠
度，则： eq k =δF
=
348EI
l
梁本身的最大动能为：
）（224348EI F -
)(x l x
x y -==）（223
max
43x l l x y -
T max =2×dx x y l m l 2
20)(21⎭
⎬⎫
⎩⎨⎧•⎰=2max 351721•y m ）
（
如果用eq m 表示简支梁的质量等效到中间位置时的大小，它的最大动能可表示为：
T max =2max
21
•y m eq
所以质量为m 的简支梁，等效到中间位置的全部质量为： m m eq 35
17=
故单自由度简支梁横向振动的固有频率为：
ωn =
eq
eq m k =
3
171680ml
EI
m
k
图1 简支梁的单自由度模型
二、 双自由度简支梁的振动特性
如图2，将简支梁简化为双自由度模型，仍假设在简支梁中间位置作用载荷，根据对称性，等效质量相等,因此只要求出在3/l 处的等效质量即可。

在6/l 至2/l 之间积分，利用最大动能进行质量等效，略去小量得：
m m eq 258≈
所以，质量矩阵为：
⎥⎦⎤⎢⎣⎡=→
1001258m m
双自由度简支梁的柔度矩阵：
在b=3/2l 处作用单位力，挠曲线方程为：）（222
6EI b -
)(b x l l
x x y --=则3/l 处的变形为：δ712=a ，同理可求：δ721=a ，δ82211==a a ，其中EI
l 4863
=δ。

所以，柔度矩阵为：
⎥
⎦⎤⎢⎣⎡=→
8778δa
动力矩阵：
⎥⎦⎤⎢⎣⎡=
→
877825
8δm D
令特征行列式为零，得到频率方程为：
=-=∆→
→
D I λ
其中，2
1
ωλ=
，将上式整理得：
116158177812=+-=---=
∆a a a
a a a
-
其中，25
82582
δωλδm m a ==。

解上述方程的根为：
15
1
1=
a ,δωm 245
1= 12=a ,δ
ωm 825
2=
由式→
→
→
→=-0)(）（i i X D I λ，2,1=i
其中⎭
⎬⎫
⎩⎨⎧=→
)(2)(1i i i X X X
）（，分别将1ω、2ω代入上式，得 第一、二阶主振型分别为：
⎭⎬⎫⎩⎨⎧=→
11)1(1
1X X
）（， ⎭
⎬⎫⎩⎨⎧=→1-1)1(12X X ）（
图2 简支梁的双自由度模型
三、 三自由度简支梁的振动特性
如图3，将简支梁简化为三自由度模型，按照双自由度类似的等效思想，可得等效质量：
m m m 41m 231≈
≈=
因此，质量矩阵为：
⎥⎥⎥⎦
⎤⎢⎢⎢⎣⎡=→
1000100014m m 由机械振动中文教材例6.6可知，系统的柔度矩阵为：
⎥⎥⎥⎦⎤
⎢⎢⎢⎣⎡=→
91171116117119δa
其中，EI
l 7683=δ。

动力矩阵：
⎥⎥⎥⎦
⎤
⎢⎢⎢⎣⎡=→
911711161171194δm D
令特征行列式为零，得到频率方程为： 0=-=∆→
→D I λ
其中，2
1
ω
λ=
，将上式整理得：
091117111611171191=---------=∆a
a a a a a a
a
a
其中，442
δωλδm m a ==。

利用Matlab 软件,求解上述方程的根为：
0317
.01=a ,δωm a 1
14=
5.02=a ,δωm a 2
24=
254.23=a ,δ
ωm a 3
34=
由式→
→
→
→=-0)(）（i i X
D I λ，3,2,1=i
其中⎪⎭
⎪⎬⎫⎪⎩⎪⎨⎧=→
)(3)(2)
(1i i i i X X X X ）（，分别将1ω、2ω、3ω代入上式，得 第一、二、三阶主振型分别为：
⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧=→
121)1(11X X
）（， ⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧=→1-01)2(12X X ）（, ⎪⎭
⎪⎬⎫⎪⎩⎪⎨⎧=→
12-1)1(13X X ）（
图3 简支梁的三自由度模型
四、 十自由度简支梁的数值方法
将简支梁简化为十自由度模型（如图4）。

图4 简支梁的十自由度模型
通过在一点施加单位力，计算其余点的挠度，可得柔度矩阵：。

为挠度变形矩阵，如表，其中1,63→
→
→
==y EI
l y a δδ
0.0137 0.0240 0.0306 0.0339 0.0344 0.0324 0.0284 0.0227 0.0158 0.0081 0.0240 0.0443 0.0579 0.0650 0.0664 0.0628 0.0552 0.0443 0.0309 0.0158 0.0306 0.0579 0.0787 0.0904 0.0934 0.0891 0.0787 0.0633 0.0443 0.0227 0.0339 0.0650 0.0904 0.1071 0.1131 0.1093 0.0973 0.0787 0.0552 0.0284 0.0344 0.0664 0.0934 0.1131 0.1229 0.1212 0.1093 0.0891 0.0628 0.0324 0.0324 0.0628 0.0891 0.1093 0.1212 0.1229 0.1131 0.0934 0.0664 0.0344 0.0284 0.0552 0.0787 0.0973 0.1093 0.1131 0.1071 0.0904 0.0650 0.0339 0.0227 0.0443 0.0633 0.0787 0.0891 0.0934 0.0904 0.0787 0.0579 0.0306 0.0158 0.0309 0.0443 0.0552 0.0628 0.0664 0.0650 0.0579 0.0443 0.0240 0.0081 0.0158 0.0227 0.0284 0.0324 0.0344 0.0339 0.0306 0.0240 0.0137
表1 十自由度挠度变形矩阵→
y
十自由度简支梁为十个集中质量的振动模型，每个质量都近似等于
m 11
1
，因此，质量矩阵为： ⎥⎥⎥
⎥
⎥⎥
⎥
⎦
⎤⎢⎢⎢⎢⎢⎢⎢⎣⎡=
→4484476O 10
10000000010000111m m
动力矩阵为：
→
→
=y
m D 11
δ
下面，用如下几种方法计算十自由度简支梁的固有频率与振型。

1、邓克莱法
利用邓克莱法求基频（比准确值小）：
n
nn m a m a m a +++≈Λ2221112
1
1
ω
因此，将柔度矩阵主对角线上各元素相加并乘以
m 11
1
，可求得：
3
19015ml EI
m δ≈≈
ω
2、瑞利法
（1）瑞利第一商
柔度矩阵求逆得刚度矩阵：
→-→→
==z y k δ
δ3131010，其中，→
z 矩阵见表2。

2.0433
-1.9003 0.7778 -0.2649 0.1647 0.0356 -0.2817 0.2834 -0.0446 -0.0926 -1.9003 2.9228 -2.4675 1.4375 -0.6862 0.0299 0.5193 -0.6226 0.3193 -0.0446 0.7778 -2.4675 3.8387 -3.3468 1.6331 -0.1417 -0.6684 0.8588 -0.6226 0.2834 -0.2649 1.4375 -3.3468 4.2339 -2.9123 0.778 0.4309 -0.6684 0.5193 -0.2817 0.1647 -0.6862 1.6331 -2.9123 3.4193 -2.2932 0.778 -0.1417 0.0299 0.0356 0.0356
0.0299 -0.1417 0.778 -2.2932 3.4193 -2.9123 1.6331 -0.6862 0.1647
-0.2817 0.5193 -0.6684 0.4309 0.778 -2.9123 4.2339 -3.3468 1.4375 -0.2649 0.2834 -0.6226 0.8588 -0.6684 -0.1417 1.6331 -3.3468 3.8387 -2.4675 0.7778 -0.0446 0.3193 -0.6226 0.5193 0.0299 -0.6862 1.4375 -2.4675 2.9228 -1.9003 -0.0926 -0.0446 0.2834 -0.2817 0.0356 0.1647 -0.2649
0.7778 -1.9003 2.0433
表2 矩阵→
z 各元素
假设力作用在简支梁中间位置而得到各点的静变形，可以表示为：
[]T
A 11.933 2.7333.3313.6633.663 3.3312.7331.9331κ=→
其中，EI
l 483
=κ。

因此，可以假设振型：
[]T
A 11.933 2.7333.3313.6633.663 3.3312.7331.9331=→
则由瑞利第一商公式：→
→→→→→=
A
M A A K A A R T
T
I )(，可得：
δm
δm m A R I 46.1646.164965.111)(1≈=⨯=
ωδ， （2）瑞利第二商
同样假设力作用在简支梁中间位置，由瑞利第二商公式：
→
→
→→→→→→∆=
A
M M A A
M A A R T
T
Ⅱ)( 可得: δm
.δm m A R Ⅱ24
1624.164762.111)(2≈=⨯=ωδ， 瑞利法中，→
M 代表质量矩阵，→
K 代表刚度矩阵，→
∆代表柔度矩阵，
→
A 为模态向量。

3、李兹法
将十自由度简支梁缩减为三自由度，假设振型为：
[]T
11.9 2.73.33.73.7 3.32.71.911=→
ψ []T 11.8 2.510.20.2- 1-2.5-1.8-1-2=→ψ
[]T 1-2- 1-121 01-2-1-3=→
ψ
则可求出：
由式 →
→
*
→
*
→*
=0-2
A M K ）（ω
，得：
0017
.01=a ,1415.02=a ,2959.03=a
其中，⎥⎥⎥
⎦
⎤⎢⎢⎢⎣⎡⨯=⎥⎥⎥⎦
⎤⎢⎢⎢⎣⎡=***→
23222133211011ωωωδm a a a a ，因此可得： δ
δωδδωδωm m a m m a m m δa 9
.32541011,7.15561011,7.181011333232131
=
⨯⨯==⨯⨯==⨯⨯=***
以及：
⎪⎪⎪⎭
⎫ ⎝⎛=⎪⎪⎪⎭⎫ ⎝⎛=⎪⎪⎪⎭⎫ ⎝⎛=→*→*→
*8372.05468.0-0122.0,6271.07788.00135.0,0349.00099.09993.0-)3()2()
1(A A A
所以系统的前三阶主振型的近似为：
4、矩阵迭代法
单位力作用在简支梁中间位置得到各点的挠度变形，将首项化一，得：
[]
T
11.933 2.7333.3313.6633.663 3.3312.7331.9331κμ=
⨯==→→→→
*11
m
M M T
ψψ
72.96 0 -0.6 0 23.06 1.2 -0.6 1.2 18
⨯==→→→→
*
δψψ3
10K K T
0.1279 0 0.1354
0 4.3927 -1.2322 0.1354 -1.2322 4.2842
其中，δκ8
1
483==EI l 。

因此可假设振型：
[]T
A 11.933 2.7333.3313.6633.663 3.3312.7331.93310=→
利用矩阵迭代求第一阶固有频率和主振型：
[]T m δ
A D 0.69661.3381 1.87072.25232.45032.4503 2.25231.87071.33816966.0110=
→
→
[]
T
m 11.9209 2.68553.2333
3.51753.5175 3.23332.68551.9209
1 0633.0⨯=δ
[]T
A 11.9209 2.68553.23333.51753.5175 3.23332.68551.920911=→
[]T m δ
A D 0.67781.3018 1.81962.19052.38282.3828 2.19051.81961.30186778.0111=
→
→
[]
T
m 11.9206 2.68463.2318
3.51553.5155 3.23182.68461.9206
1 0616.0⨯=δ
[]T
A 11.9206 2.68463.23183.51553.5155 3.23182.68461.920612=→
[]T m δ
A D 0.67751.3013 1.81892.18952.38182.3818 2.18951.81891.30136775.0112=
→
→
[]
T
m 11.9207 2.68473.2318
3.51553.5155 3.23182.68471.9207
1 0616.0⨯=δ
[]
T
A 11.9207
2.6847
3.2318
3.5155
3.5155 3.23182.68471.9207
13=→
由上，仅3次矩阵迭代后，→
3A 与→
2A 基本相等，因此可以认为系统的第一阶主振型为：
[]T
A 11.9207 2.68473.2318
3.51553.5155 3.23182.68471.9207
11=→）
（ 第
一阶固有频率为：
δ
δωm m a 23.160616.01131≈==
5、雅可比法
根据雅可比法原理，依次找出上三角非对角线上（考虑对称性）的最大元素，利用公式jj
ii ij d d d -=
22tan θ得到θ值，代入旋转矩阵，可得：
=••=→
→→451R R R Λ
0.349 -0.056 -0.417 0.215 0.396 -0.394 -0.217 0.415 0.058 -0.347 0.396 0.323 -0.134 -0.420 -0.212 0.212 0.420 0.134 -0.324 -0.395 0.398 -0.369 -0.032 0.334 -0.185 -0.215 0.340 0.002 -0.436 0.455 0.321 0.424 0.230 -0.122 -0.387 -0.387 -0.122 0.231 0.423 0.322 0.161 -0.324 0.419 -0.346 0.087 0.244 -0.460 0.454 -0.282 0.102 0.232 0.388 0.422 0.322 0.120 -0.120 -0.322 -0.422 -0.388 -0.232 0.104 -0.231 0.354 -0.442 0.466 -0.424 0.351 -0.263 0.155 -0.053 0.120 0.231 0.322 0.388 0.422 0.422 0.388 0.322 0.231 0.120 0.439 0.101 -0.383 -0.219 0.304 0.337 -0.236 -0.391 0.127 0.415 -0.404 0.453 -0.158 -0.184 0.324 -0.242 0.044 0.217 -0.448 0.399
则：
=••=→
→→→R D R D
T
)
45(
⨯11
δm 0.0005 0 0 0 0 0 0 0 0 0
0 0.0027 0 0 0 0 0 0 0 0 0 0 0.0003 0 0 0 0 0 0 0 0 0 0 0.0084 0 0 0 0 0 0 0 0 0 0 0.0001 0 0 0 0 0 0 0 0 0 0 0.0424 0 0 0 0 0 0 0 0 0 0 0.0001 0 0 0 0 0 0 0 0 0 0 0.6775 0 0 0 0 0 0 0 0 0 0 0.0011 0 0 0 0 0 0 0 0 0 0 0.0002
由此得到系统的固有频率：
δ
δδδδδδδδδωm m m m m m m m m m 110000,110000,55000,67.36666,22000
,10000
,07.4074,52.1309,43.259,24.16101=
Λ
对应各特征值的特征向量即为振型，即→
R 的列向量为各阶振型。

6、子空间迭代法
取假设振型：
⎥
⎥
⎥
⎥
⎥
⎥
⎥
⎥
⎥
⎥
⎥
⎥
⎥
⎥
⎦
⎤
⎢
⎢
⎢
⎢
⎢
⎢
⎢
⎢
⎢
⎢
⎢
⎢
⎢
⎢
⎣
⎡
=
→
1.5
1-
1
1
0.5
2-
1.8
1.9
1.5
-
1-
2.5
2.7
1-
1
1
3.3
0.5
-
2
0.2
3.7
2
1
0.2
-
3.7
1
1-
3.3
1
2.5
-
2.7
2-
2-
1.8
-
1.9
1-
1-
1-
1
A
由动力矩阵迭代得到：
⨯
=
∆
→
→
→
11
δ
m
A
M
0.6941-0.0460.0341-0.0176
1.3333-0.07680.0706-0.0287
1.8642-0.08340.108-0.0312
2.2447-0.06290.138-0.0304
2.4421-0.02340.1537-0.0318
2.44210.02340.1506-0.0361
2.24470.06290.1282-0.0388
1.86420.08340.0932-0.0347
1.33330.07680.0568-0.0237
0.69410.0460.0258-0.011
将各列分别归一化得：
T
⎥
⎥
⎥
⎥
⎦
⎤⎢
⎢
⎢
⎢
⎣
⎡
=
→
2822
.0
6108
.0
8943
.0
1
9317
.0
8209
.0
7835
.0
8041
.0
7397
.0
4536
.0
1679
.0
3696
.0
6064
.0
8341
.0
9798
.0
1
8979
.0
7027
.0
4593
.0
2219
.0
5517
.0
9213
.0
1
7548
.0
2802
.0
2802
.0-
7548
.0-
1-
9213
.0-
5517
.0-
2842
.0
5460
.0
7633
.0
9191
.0
1
1
9191
.0
7633
.0
5460
.0
2842
.0
Ⅰ
ψ
求得→*
Ⅰ
M和
→
*
Ⅰ
K分别为：
⎥
⎥
⎥
⎥
⎦
⎤
⎢
⎢
⎢
⎢
⎣
⎡
-
-
=
=
→
→
→
*
8041
.5
0899
.5
0721
.0
6327
.5
0899
.5
7483
.4
2625
.0
1344
.5
0721
.0
2625
.0
5995
.5
6327
.5
1344
.5
6129
.5
11
Ⅰ
m
M
M T
Ⅰ
Ⅰ
ψ
ψ
⎥⎥⎥⎥⎦
⎤⎢⎢⎢⎢⎣⎡--=
=→→→*041.0001.00007.00083.0001.00119.00063.00076.00007.00063.01324.000083.00076.000083
.0103Ⅰδψψ
ⅠⅠ
K K T 再由李兹法得特征值问题为：
→
→*→
*→
*=0-2
Ⅰ
Ⅰ
Ⅰ
）（A M K ω
解出：
00147
.01=a ,0236.02=a ,1198.03=a ,3835.04=a
其中，⎥⎥⎥
⎦
⎤⎢⎢⎢⎣⎡⨯=⎥⎥⎥⎦
⎤⎢⎢⎢⎣⎡=***→
23222133211011ωωωδm a a a a ，相应的主振型为： ⎪⎪⎪⎪⎪⎭
⎫ ⎝⎛-=⎪⎪⎪⎪⎪⎭⎫ ⎝⎛---=⎪⎪⎪⎪⎪⎭⎫ ⎝⎛=⎪⎪⎪⎪⎪⎭⎫ ⎝⎛=→*→*→*→*0479.00854.09872.01262.03275.051.00195.07952.0,0268.07504.0-0355.0-6595.0,0001.0-0044.00002.01)4()3()2()1(Ⅰ
ⅠⅠⅠ
，A A A A
所以：
=→*→
Ⅰ
ⅠA ψ
0.285
0.0527 -0.0249 -0.5397 0.5477 0.0679 -0.0244 -0.9033 0.7661 0.0332 0.0048 -0.985 0.9229 -0.0198 0.0311 -0.7463 1.0042 -0.0589 0.0222 -0.2785 1.0042 -0.0607 -0.0147 0.279 0.9229 -0.0196 -0.0367 0.7477 0.7661 0.037 -0.0146 0.9855 0.5477 0.0665 0.0277 0.901 0.285
0.0495
0.0368 0.5365
各列分别归一化后，得：
=→
ⅠA
0.283808 0.776141
-0.67663 -0.54764 0.545409 1
-0.66304 -0.91659 0.762896 0.488954 0.130435 -0.99949 0.91904
-0.29161 0.845109 -0.75728 1 -0.86745 0.603261 -0.2826 1 -0.89396 -0.39946 0.283105 0.91904 -0.28866 -0.99728 0.758701
0.762896
0.544919
-0.39674
1
0.545409 0.979381 0.752717 0.914257 0.283808 0.729013 1 0.544394
重复上述过程进行第二次迭代。

由：
⨯=
∆→
→→11
δ
m A M Ⅰ 0.1926
0.0065 -0.0021 -0.0232 0.3699 0.0085 -0.0017 -0.0388 0.5171 0.0048 0.0008 -0.0422 0.6225 -0.002 0.0025 -0.0322 0.6771 -0.007 0.0015 -0.012 0.6771 -0.007 -0.001 0.0121 0.6225 -0.0019 -0.0024 0.0322 0.5171 0.0049 -0.0007 0.0422 0.3699 0.0086 0.0022 0.0388 0.1926
0.0064
0.0026
0.0232
归一化后得：
=→
Ⅱψ
0.284448 0.755814 -0.80769 -0.54976 0.5463 0.988372 -0.65385 -0.91943
0.763698 0.55814 0.307692 -1 0.919362
-0.23256 0.961538 -0.76303 1 -0.81395 0.576923 -0.28436 1
-0.81395 -0.38462 0.28673 0.919362 -0.22093 -0.92308 0.763033
0.763698 0.569767
-0.26923 1
0.5463 1
0.846154
0.919431 0.284448
0.744186
1
0.549763
则有：
⨯=
=→
→
→
*11m M M T ⅡⅡ
Ⅱψψ
5.6156 0.3295 0.4168 0.0024 0.3295 5.166 0.1758 0.0229 0.4168 0.1758 5.2204 0.0837 0.0024 0.0229 0.0837 5.6227
⨯
=
=→
→
→*δ
ψψ
3
10ⅡⅡ
Ⅱ
K K T
0.0083 0.0009 0.0006 0 0.0009 0.6107 0.0157 0.0004 0.0006
0.0157 1.9186 0.0023 0
0.0004 0.0023 0.1328
由：
→
→*→
*→*=0-2
Ⅱ
Ⅱ
Ⅱ
）（A M K ω 解出： 00147.01=a ,0236.02=a ,1186.03=a ,375.04=a
其中，⎥⎥⎥
⎦
⎤⎢⎢⎢⎣⎡⨯=⎥⎥⎥⎦
⎤⎢⎢⎢⎣⎡=***→
23222133211011ωωωδm a a a a ，相应的主振型为： ⎪⎪⎪⎪⎪⎭
⎫ ⎝⎛=⎪⎪⎪⎪⎪⎭⎫ ⎝⎛--=⎪⎪⎪⎪⎪⎭⎫ ⎝⎛=⎪⎪⎪⎪⎪⎭⎫ ⎝⎛=→*→*→*→*0145.09968.0-031.00722.00043.00017.09983.00581.0,10002.0-0003.00004.0-,000006.0-1)4()3()2()1(ⅡⅡⅡⅡ
，A A A A
=→*→
Ⅱ
ⅡA ψ
0.284
-0.739 0.8411 -0.5495 0.5457 -0.9578 0.7085 -0.9193 0.7633 -0.5176 -0.2488 -1.0002 0.9195 0.2806 -0.9104 -0.7637 1.0005 0.8685 -0.5323 -0.2851 1.0005 0.8726 0.4345 0.2862 0.9195 0.2789 0.9907 0.7628 0.7633 -0.5197 0.3556 0.9999 0.5457 -0.9641 -0.7597 0.9193 0.284
-0.7258
-0.9453 0.5496
各列分别归一化后，得：
=→
ⅡA
0.283858 0.766518 0.848996 0.54939 0.545427 0.993465 0.715151 0.919116
0.762919 0.536874 -0.25114 1
0.91904
-0.29105 -0.91895 0.763547 1 -0.90084 -0.5373 0.285043 1 -0.90509 0.438579
-0.28614 0.91904 -0.28929 1
-0.76265 0.762919 0.539052
0.358938 -0.9997 0.545427 1
-0.76683 -0.91912 0.283858
0.752826
-0.95417 -0.54949
进行第三次迭代：
⨯=
∆→
→→11
δm A M Ⅱ 0.1926 0.0065 0.0024 0.0233 0.37 0.0086 0.002 0.0389 0.5171 0.0049 -0.0008 0.0424 0.6225 -0.002 -0.0026 0.0323 0.6772 -0.0072 -0.0014 0.0121 0.6772 -0.0072 0.0013 -0.0121 0.6225 -0.002 0.0026 -0.0323 0.5171 0.0049 0.0009 -0.0424 0.37
0.0086
-0.002
-0.0389
0.1926 0.0065 -0.0025 -0.0233
归一化得：
=→
Ⅲψ
0.284448 0.755814 -0.80769 -0.54976 0.5463 0.988372 -0.65385 -0.91943 0.763698 0.55814 0.307692 -1 0.919362 -0.23256 0.961538 -0.76303
1 -0.81395 0.576923 -0.28436 1 -0.81395 -0.3846
2 0.28673
0.919362 -0.22093 -0.92308 0.763033 0.763698 0.569767 -0.26923 1 0.5463 1 0.846154 0.919431 0.284448 0.744186 1 0.549763 ⨯=
=→
→
→*11m M M
T
ⅢⅢⅢ
ψψ
5.6149 0.2908 -0.02 0 0.2908 5.3018 0.025 0 -0.02 0.025 5.7145 -0.0267
0 0 -0.0267 5.6109 ⨯
=
=→
→
→*δ
ψψ3
10ⅢⅢⅢ
K K T
0.0083 0.0008 0 0 0.0008 0.6267 0.0031 0
0 0.0031 2.1067 -0.0008 0 0 -0.0008 0.1325
由：
→
→*→*→
*
=0-2
Ⅲ
Ⅲ
Ⅲ
）（A M K ω 解出：00147.01=a ,0236.02=a ,1185.03=a ,3687.04=a
其中，⎥⎥⎥
⎦
⎤⎢⎢⎢⎣⎡⨯=⎥⎥⎥⎦
⎤
⎢⎢⎢⎣⎡=***→
23222133211011ωωωδm a a a a ，相应的主振型为： ⎪⎪⎪⎪⎪⎭
⎫ ⎝⎛=⎪⎪⎪⎪⎪⎭⎫ ⎝⎛-=⎪⎪⎪⎪⎪⎭⎫ ⎝⎛=⎪⎪⎪⎪⎪⎭⎫ ⎝⎛=→*→*→*→
*10001.0000047.010049.00038.0,009987.00511.0-,000006.0-1)4()3()2()1(ⅢⅢⅢⅢ
，A A A A
由于→*Ⅲ
A 近似于单位矩阵，所以有：
→
→*→
→
==ⅢⅢ
ⅢⅢψψA A
本次迭代的结果与第一次迭代的结果相差不大，所以可以结束迭代。

故系统的前二阶固有频率及相应的主振型的近似分别为：
δ
ωδωδωδωm m m m 7
.40055
.1303,6.259,17
.164321====
，
==→→→→→
)()
4()
3()2()1(A A
A A A Ⅲ 0.284448 0.755814 -0.80769 -0.54976 0.5463 0.988372 -0.65385 -0.91943 0.763698 0.55814 0.307692 -1 0.919362
-0.23256 0.961538 -0.76303 1 -0.81395 0.576923 -0.28436 1
-0.81395 -0.38462 0.28673 0.919362 -0.22093 -0.92308 0.763033
0.763698 0.569767
-0.26923 1 0.5463 1 0.846154
0.919431 0.284448 0.744186 1 0.549763
五、 连续分布简支梁的振动特性
（1）理论解
（2）数值解。