第10章 结构动力学
习 题
10-5 试确定图示各体系的动力自由度,忽略弹性杆自身的质量。 (a) (b)
EI 1=∞
EI
m
y
?
分布质量的刚度为无穷大,由广义坐标法可知,体系仅有两个振动自由度y ,?。 (c)
(d)
在集中质量处施加刚性链杆以限制质量运动体系。有四个自由度。
10-8 图示结构横梁具有无限刚性和均布质量m ,B 处有一弹性支座(刚度系数为k ),C 处有一阻尼器(阻尼系数为c ),梁上受三角形分布动力荷载作用,试用不同的方法建立体系的运动方程。
解:1)刚度法
该体系仅有一个自由度。
可设A 截面转角a 为坐标顺时针为正,此时作用于分布质量m 上的惯性力呈三角形分布。其
端部集度为..
ml a 。
取A 点隔离体,A 结点力矩为: (3)
121233
I M ml a l l mal =???=
由动力荷载引起的力矩为:
()()2121
233
t t q l l q l ??= 由弹性恢复力所引起的弯矩为:.21
33
la k l c al ?
?+ 根据A 结点力矩平衡条件0I p s M M M ++=可得:
(
)3 (322)
1393
t q l ka m al l c al ++=
整理得:()
.
..
33t q ka c a m a l l l
+
+= 2)力法
.
c
α
解:取AC 杆转角为坐标,设在平衡位置附近发生虚位移α。根据几何关系,虚
功方程为:() (20111)
0333
l t q l l k l l l c m x xdx ααααααα-?-?-?=?
则同样有:()
.
..
33t q ka c a m a l
l l
+
+=。 10-9 图示结构AD 和DF 杆具有无限刚性和均布质量m ,A 处转动弹簧铰的刚度系数为k θ,C 、E 处弹簧的刚度系数为k ,B 处阻尼器的阻尼系数为
c ,试建立体系自由振动时的运动方程。
解:
取DF 隔离体,
0F
M
=∑:
..2220.23
223
24
a R a mx dx ka R ma ka αα
αα
?=+?=+?
取AE 隔离体:
0A M =∑
..
.
32220
430a
k mx dx ca ka Ra θαααα++++=?
将R 代入,整理得: ..
32
251504
R ma ka k θ
ααα=+
+= 10-10 试建立图示各体系的运动方程。 (a)
解:(1)以支座B 处转角作为坐标,绘出梁的位移和受力图如下所示。图中惯性力为三角形分布,方向与运动方向相反。
(t )
..α
(2)画出p M 和1M 图(在B 点处作用一附加约束)
()324
t l M α-()
t p
M
3EI l
1
M
(3)列出刚度法方程
113EI
k l
=
,()..3124p t m R l M α=-,1110p k R α+=
代入1p R 、11k 的值,整理得:()..43
2472t M EI
m l l αα+=
(b) 解:
11
=
21P =2
l
1M 图 2M 图
试用柔度法解题
此体系自由度为1 。设质量集中处的竖向位移y 为坐标。 y 是由动力荷载()p t F 和惯性力矩I M 共同引起的。
11112()p t y M F δα=+
l
l 2
m (t )
l 2 l
2
由图乘法:
3
21112233l l l EI EI
δ=?=,312/252622248l l l l l l EI EI δ??=??+?=
??? 惯性力矩为..
m y l -,()33..5348p t l l y m yl F EI EI
??=?-+ ???
经整理得,体系运动方程为:()..
3
35
16p t EI m y y F l +
=。
10-11 试求图示各结构的自振频率,忽略杆件自身的质量。 (a)
解:
2
1M 图
图乘得:3
1111225222223236a a a f a a a a EI EI
??=???
??+??
?=
??? ω=
=(b)
解:此体系为静定结构,内力容易求得。
在集中质量处施加垂直力P ,使质量发生竖向单位位移,可得弹簧处位移为2
3
。
由此根据弯矩平衡可求得4
9
P
k =
。 ω=
= (c)
解:可以将两个简支梁视为两个并联的弹簧。
上简支梁柔度系数为
()3
32486l l EI EI =
下简支梁柔度系数为3
96l EI
于是两者并联的柔度系数为33
1696
102l EI EI EI l
δ==+并,ω==
(d)
解:在原结构上质量运动方向加上一根水平支杆后,施加单位水平位移后画得弯矩
图如下。
水平支杆中力为
33013EI l
,即113
30
13EI
k l
=。,ω=
(e)忽略水平位移
l 2 l 2 l
2
l 2
2a
a
a
解:
1
M图
22
11
245527
2213
362
a a a
f a
EA EA EA
????
=??+??+?=
? ?
????
ω==
(f)
解:
3
32
3
32
1
M图2
M图M图
3
1312331323162130.014974
3223323221933219364
l
l l l l l l l
EI EI
δ??
=???+????+??=
?
??
ω=
10-15 设已测得某单自由度结构在振动10周后振幅由1.188mm减小至0.060mm,试求该结构
的阻尼比ξ。解:0475
.0
06
.0
188
.1
ln
20
1
ln
2
1
=
=
≈
+
π
π
ξ
n
k
k
y
y
n
10-16 设有阻尼比ξ=0.2的单自由度结构受简谐荷载F P(t)=F tθ
sin作用,且有ω
θ75
.0
=。
若阻尼比降低至ξ=0.02,试问要使动位移幅值不变,简谐荷载的幅值应调整到多大?
解:
2
2
2
2
2
2
2
4
1
1
ω
θ
ξ
ω
θ
ω
+
??
?
?
?
?
-
?
=
m
F
A
已知ξ从0.2降低至0.02. ω
θ75
.0
=,t
F
Fθ
sin
1
=,A不变。
1
2
2
2
2
2
2
1827
.0
16
9
02
.0
4
16
9
1
16
9
2.0
4
16
9
1
F
F
F
F
=
?
?
?
+
?
?
?
?
?
-
?
?
+
?
?
?
?
?
-
=
F简谐荷载的幅值应调整到0.827F。
10-19 试求图示梁在简谐荷载作用下作无阻尼强迫振动时质量处以及动力荷载作用点的动位
移幅值,并绘制最大动力弯矩图。设
3
6
ml
EI
=
θ。
(a)
解:由力法可知,单位荷载作用在B点引起
3
3
l
EI
位移。
tθ
sin
l
l
2
l
2
ω==
θ=
()3222
1sin sin 31t F Fl y t t EI m θθθωω=
?=--
即幅值为3
3Fl EI
当幅值最大时,弯矩也最大。
Fl
max M 图
(b)
解:
1M 图 2M 图
(1)求结构运动方程
如所示弯矩图,图乘后,333
112212215,,24348l l l f f f f EI EI EI
====
()..11121112..3
sin sin 245sin 2I t C y f F f F t f m y f F t
EI F y y t
m ml θθθ??
=+=-+ ???
+=
其中2*3
245,2EI P F ml
ω==,稳态解:()*
2
22
331
sin 1512 =sin 124145 =sin 36t C
P y t
m Fl
t EI Fl t
EI
θωθω
θθ=?-?- 所示结构的运动方程为()35=sin 36t C Fl y t EI θ,C 点最大动位移幅值为
3
536Fl EI
(2)求B 点的动位移反应
()()..
21222122sin sin I t B t B y f F f P t f m y f P t θθ??
=+=-+ ???
()*
2
22
1
sin 1t B
P y t m θωθω
=?-,()*
..
2
2
22
1
sin 1
t B
P y t m θθωθω
=-?
-
2l 2
l t θ sin
()()3
2*212222232322
2322222
3
5=sin 361sin 1551 =sin 48231251 =1sin 33217132 =
3t C t B Fl y t
EI y f P Pf t
l l
P P t EI EI Pl t EI Pl EI θθθωθωθθωθωθθωθωθ???
?
?? ?
?? ?=??+?? ?
- ??????
?
??
?????+????-?????? ? ???+ ?- ?
??-22233
sin 11214 =sin 31283121 =sin 288t Pl t
EI Pl t
EI
ωθθωθθ??
? ? ?
- ?
?
???
B 点的动位移幅值为3
121288Pl EI
(3)绘制最大动力弯矩图
22
1M 图 2M 图 ()33max 2212135122812883696A Pl EI Pl EI M Pl EI EI l l =?+?= ()
3max 21213121288192
2C Pl EI M Pl EI l =?=
121
192Pl 281
96
Pl
最大动力弯矩图
10-20 试求图示集中质量体系在均布简谐荷载作用下弹簧支座的最大动反力。设杆件为无限刚性,弹簧的刚度系数为k 。
解:
α
若()t q 为静力荷载,弹簧中反力为
ql 8
9
。 已知图示体系为静定结构,具有一个自由度。设为B 点处顺时针方向转角α为坐标。建立动力方程:
?=?+?+l xdx q l l k l m l l m l 230..
..
2
332322αααα
ααα
q k m l q l k l m 8
9
89..22
22
..
=+?=+αααααα
2
211
ω
θμ-=
则弹簧支座的最大动反力为
l 8
9112
2
?-
ωθ。
2
l 2
l l
10-21 设图a 所示排架在横梁处受图b 所示水平脉冲荷载作用,试求各柱所受的最大动剪力。已知EI =6×106N ·m 2,t 1=0.1s ,F P0=8×104N 。 (a)
解:求排架自振频率,横梁无限刚性,则各排架水平侧移相同。 可将排架柱视为三个并联的弹簧。 边柱刚度柔数331
3h EI k k =
= 中柱326h EI k =,3
12h EI
k =并 s rad N
m m N m k /645.01080006106122332
6=?????==ω
s T 73.92==
ω
π
,
3
.971
73.91.01=
=T t 数值很小 所以认为当()t P F 作用结束时,结构位移很小,弹性力忽略不计,于是根据动量守恒原理可得:
s
m v v Ft v m t t t /1051
.010821
108213141511-?=????=??=?
再根据势能守恒得:
()
m
y y ky mv st st
t 0077.0103
12110510821
2121262
352max 21=????=?????=- N k y F st Q 128310610077.06=??=?=中中,N F F Q 中Q 边6422
1
==
10-22 设图a 所示排架横梁为无限刚性,并有图b 所示水平短时动力荷载作用,试求横梁的动位移。
(a) 解:在三角形冲击荷载作用下单自由度体系的质点位移反应可分两个阶段考虑。
第一阶段(10t t
≤≤)
: ()()()()???
??
???????-??? ??=?
???
??-??? ?????
? ??=?
??
?????? ??-=???
??????
??-=
-=-=??1111112
0100022sin 2sin 21
sin 1 sin 1 sin sin 1t t t T t T y t t T t t T y t t t y t t t m F dZ
Z t t Z m F dZ Z t F m y s s s P t P t
Z P t ππππ
ωωω
ωωωωωω
求T 的过程。
2
6EI h 2
6h 2
h
1M 图
3
1124h
EI
k =,311
24mh EI m k ==
ω,EI
mh T 24223
π
ωπ
==
第二阶段(1t t >)
因为不受外力作用,所以横梁以1t 时刻的位移和速度为初始值做自由振动。
(b)
10-23 设题10-22图a 所示刚架m =4000kg ,h =4m ,刚架作水平自由振动时因阻尼引起振幅的对数递减率γ=0.10。若要求振幅在10秒内衰减到最大振幅的5%,试求刚架柱子的弯曲刚度EI 至少为何值。
解:(1)求周期数。
301
.005
.0ln 05.000=-=
?=-n e y y n Y (2)求k :k
m n t n
π
2= ()()m N t m n k n
/10223.142110100.43014159.32232
3
2
2
2?=????==?π
两柱并联
2631079.312
2m N EI k h
EI ??=?=?
10-24 设某单自由度体系在简谐荷载F P (t )= F t θsin 作用下作有阻尼强迫振动,试问简谐荷载频率θ分别为何值时,体系的位移响应、速度响应和加速度响应达到最大?
解:在简谐荷载F P (t )= F t θsin 作用下,稳态位移响应可表示为()
()αθ-=t A y t sin
其中:?
????
?
?
????
??
??
??
??-==+???? ??-?=-22122
2222212tan 411ωθωθξαμωθξωθωst y m F A
(1)使动位移最大,即使μ最大,从而得出22
22
2241ωθξωθ+???? ??-最小。 设()22
22
2241ωθξωθθ+???? ??-=f ,()222
2222414ωθξωθωθθ+???
? ??--='f 使()0='θ
f ,则221ξωθ-= (2)())cos(αθθ-='t A y t 设()222
22
22
221
411
41ωξ
ωθθωθξωθθ
θ+??
? ??-=
+???
? ??-=
g
如果使速度响应最大,则()θg 最大,设()222
1141ωξωθθθ+??
?
??-=g ,显然要求
()θ1g 最小。使:()011
12221
=??? ??--??? ??-='ωθ
ωθθθg 得ωθ=。 (3)())sin(2
αθθ--=''t A y t
F P (t )
t
F P0
t 1
O
()2
222
2222
22
222
411
1
41θωξωθ
ωθξωθθθ+??? ??-=
+???
? ??-=
h
令()2222
221411
θωξωθ
θ+??? ??-=h 显然要求()θ1h 最小。
则()0211
2
2
2
1
=--
='ω
ξθ
θh 解的:2
21ξ
ω
θ-=
10-26 试用柔度法求下列集中质量体系的自振频率和主振型。 (a) 解:
2
l
2l
1M 图 2M 图
(1)EI
l f l l l l l l EIf 42322222123222123
1111=
??????+?????=
EI
l f l l l l l l EIf 1252223222123
2222
=
???+?????=,
02112==f f
(2)振型方程
???
????=???? ??-?+?=?+???
?
??-01212500014223
12123A m EI l A A A m EI l ωω 令2
312ωλ
ml EI
=
,频率方程为:0-10 00 3=-=λ
λD ()()3
32331212312 095.110123
,100103ml EI
ml EI ml EI ml EI ==
==?==?=
--?ωωλλλλ
(3)振型图如下
第二振型
(b)
解:
体系具有两个自由度。先求柔度系数,做出单位弯矩图,由图乘法可得:
()
EI l l l l l l l EI 32132221322113
11+=
??? ????+??=δ
EI l l l l EI 6223222113
1221=
??? ?????==δδ EI
l l l l EI 622322212123
22=
??? ??????=δ
得振型方程:
()
062132123
123=+??
?
? ??-+mA EI l A m EI l ω 01626222313=???? ??-+A m EI
l mA EI l ω,令λω=?3231ml EI
λ
λ-0.707 0.7070.707 414.2-=
D ,由频率方程D=0
解得:331
576.24535.03ml EI
ml EI =?=
ω,3
32060
.16675.23ml EI
ml EI =?=
ω
1773.2707.0414.211121-=--=λA A ,1
358
.0707.0414.221222=--=λA A
(c)
解:
/2
l
1M 图 2M 图
(1)
EI l f 3311=,EI l f 1213322=,EI
l f f 1253
2112=
=
(2)振型方程
???????=???? ??-?+????? ?
?=????
? ??+????
??-0121213125012513223
1323123A m EI l A m EI l A m EI l A m EI l ωω 令2
312ωλ
ml EI
=
,频率方程为:0-13 55 4=-=λ
λD 3
32331212602
.2773.112 888
.0227.1512773
.1,227.150255217ml EI
ml EI ml EI ml EI ==
==?==?=-+-?ωωλλλλ
l l l l
(3)当227.151==λλ时,设7227.010
8
112111=-=
?=λA A
当773.12==λλ
时,设6227.010
8
122212-=-=
?=λA A
绘出振型图如下:
第一振型 第二振型
(d)
解:
12
1M 图 2M 图
EI
a
k k EI a 3
213114811/21/212161=??? ??++=
δ
EI a a k k 32121124812//21/21=
??
?
??-==δδ
EI
a k k EI a 3
21322
4811/21/212161=??? ??++=δ
频率方程为:
111122
2
2112222
1
01
m f m f m f m δωω-
=-
取3
121
,3
m ma m ma ==
代入整理得: 22444003a a λλ-
+
=其中3248
EI a m λω= 1211.045, 3.625a a λλ==
1ω=
2ω=
=振型方程为:
()11111222221112222210
10m A m A a f m A f m A δδωω???
-+?= ?????
?
????+?-= ?????
将()1,11,2i i A i ωω===代入(a )式中的第一个方程中,得:
44
111
2
1212
3
122
1
0.23010.2292
0.1351
483ma ma m EI EI A m a a ma EI δωδ--=
==? 4
111
2
22223
122
1
3.6251122.125481483
ma m EI A m a a ma EI δωδ---=
=
=?
绘出振型图如下:
a
a a
第一振型 第二振型
(e)
解:
1M 图
2M 图
3M 图
(1)3112l f EI
=,3
222l f EI =,31221336l f f f EI ===
(2)振型方程
3312323
312323
123210026100621000
6l l m A m A A EI
EI l l m A m A A EI EI l A A m A EI ωωω???
??-+?+?=?
? ??????
??????
?+?-+?=?
? ???
????????+?+?-= ?????
令32
6EI
ml λω=,频率方程为:3 1 0
1 3- 000 0 2-D λλλ
-==
1231234,2λλλωωω?=
==?=
== {}{}{}1231101 1 0001A A A ?????? ? ? ?
==-= ? ? ? ? ? ???????
振型图如下:
a
a
第一振型 第二振型
第三振型
(f)
解:
1M 图 2M 图 3M 图
333333
1122332112233231131895144,,,,,33633a a a a a a
EI EI EI EI EI EI
δδδδδδδδδ=========(2)振型方程为:
333123233312323
331232154403635811440633414914033a a a m A m A m A EI
EI EI a a a m A m A m A EI EI EI a a a m A m A m A EI
EI EI ωωω???
????-+?+??=?
? ? ????????????????
?+?-+??=?
? ? ???
???????????????+??+?-= ?
? ????????? 令32
6
EI
ml λω=
,频率方程为:2 5 32
5 16- 112
08 28 216-D λλλ
-== 123123231.8, 1.936,0.2317
λλλωωω?===?==={}{}{}1231103.469 1.390 0.6876.6400.2190.052A A A ?????? ? ? ?
===- ? ? ? ? ? ?-??????
10-27
试用刚度法求下列集中质量体系的自振频率和主振型。 (a) 解:
a a a
4m
EI=常数 l l
11
k 21
k
12
k 22
k
1M 图 2M 图
11211222222
242448,,EI EI EI
k k k k l l l =
===,23
24 -240 -24 48-2y m l y y EI ω-??== ???
127.029,40.971y y ==
,12
ωω=={}{}12
11,0.7070.707A A ????
== ? ?-????
振型图如下:
第一振型
第二振型
(b) 解:
11
k
1F 图 2F 图
(
112221214
222EA EA EA k k l l l EA k k l
==+=
=
==
振型方程:
21
22
1200
m A A A m A ωω??-=??????
?+-=?????
令24m l
EA
ωλ=
,频率方程为:0D λ
λ
=
=
12124,4 λλωω?
==+?=
=()()1
211,11A A ????== ?
?-????
(c)
解:
k= EI l 3 l
l
1M 图 2M 图
作出附加连杆移动单位位移的弯矩图
112334i EI k k l l =
+=,12213EI k k l ==-,22
23
34i EI
k k l l =+= 列出频率方程:2
11112
22122
2 k 0
k m D k k m ωω-=
=-
解得:21322335EI ml EI ml ωω?=????=??结构自振频率分别为:12ωω?=????
=??
求第一振型:令111A =得211A =
求第二振型:令121A =得221A =- 结构的振型向量形式为:
()()1
211,11A A ????== ? ?-????
振型图如下:
1
1
1
1
第一振型 第二振型 (d)
解:
22
k
1M 图 2M 图
12210k k ==,112152i k l =
,222
8i
k l = 列振型方程:()()()1
2150
*160
y A y A -=???-=??
其中322ml y EI ω=
列频率方程并求解:
()()15 00151600 16-y
y D y y -=
=?--=
1215,16y y ==
12
ωω==求振型
将11115,1y A ==代入方程组(*)中得:210A =,即()
110A
??
= ??? 将22216,1y A ==代入方程组(*)中得:220A =,即()
201A ??= ???
振型图如下:
第一振型 第二振型
10-28 试说明在应用多自由度体系强迫振动的振幅方程(10-66)和(10-71)时,对动力荷载的性质、特点和作用位置分别有何要求?
10-29 试说明为什么可以将惯性力幅值与简谐荷载幅值同时作用在体系上,按静力学方法计算体系的动内力幅值。
10-30 试求图示结构B 点的最大竖向动位移(max)yB ?,并绘制最大动力弯矩图。设均布简谐荷载频率3ma EI =
θ,B 点处弹性支座的刚度系数3
a EI k =,忽略阻尼的影响。
解:
1M 图 P M 图
画1,p M M 图
3
111121115222322212a a a f a EI k EI
??=?????+??=
???
422211111121111224342234244p
a a
qa a qa qa a qa qa EI
k EI ?????=???+?+????+??= ???????
列出方程得:
334
150124a a qa I EI EI EI ??-+
= ???
解得:137
I qa =
()
333
max 3111372428yB a a qa qa qa EI EI EI
?=??+?=
根据公式11p M M I M =+画出最大动力弯矩图。
M 图
10-31 图示结构在B 点处有水平简谐荷载t t F θsin kN 1)(P ?=作用,试求集中质量处的最大水
平位移和竖向位移,并绘制最大动力弯矩图。设3
ml EI =
θ,忽略阻尼的影响。
a
q s θ t
解:作出12M M 、图
2
1M 图 2M 图 P M 图
11112132
222222233EI EI EI δ=
?????+???=
,21121142222EI EI δδ==????= 221128
222233EI EI
δ=
?????=
,11142222p F F EI EI ?=????= 21128222233p F F EI EI ?=
?????=,代入惯性力幅值方程:3123
12324403488033ml F
I I EI mEI EI EI ml F
I I EI EI mEI EI ??-++= ???
??+-+= ???
解得:12185
,1717I KN I KN =-
=-,1212220.941,0.261I I A mm A mm m m θθ
==-==- 将以上求得最大惯性力1I 、2I 和动力荷载,同时作用于结构,可得最大动力弯矩图:
12
17
KN m ?36
17
KN m ?
M 图
10-32 图示刚架各横梁为无限刚性,试求横梁处的位移幅值和柱端弯矩幅值。已知m =100t ,l =5m ,EI =5×105kN ·m 2;简谐荷载幅值F =30 kN ,每分钟振动240次;忽略阻尼的影响。
解:
31
k 21
k
32
k 12
33
k 23
k 13
层间刚度设为k ,3
24EI
k l =
112212212332332,,k k k k k k k k k k ======-=
22240
86060
n ππθπ?=
== F=30KN l=5m 动位移幅值方程为:
21233212333321333482420
2448241.524240
EI EI m A A l l EI EI EI A m A A F l
l l EI EI A m A l
l θθθ???
--= ???
?????
-+--=? ???????
-+-=? ????
将具体数值代入,解得:
1230.1353,0.0926,0.2710A mm A mm A mm =-=-=-
底柱柱端弯矩幅值:5
31133
121251050.13531016.236225
EI l M A KN m l ??=??=-???=? 中柱柱端弯矩幅值
()()53
2213126100.13530.092610 5.12425EI l M A A KN m l
?=-??=-??=?
l
l
l
m 1=2m
m 2=1.5m m 3=m
t
F θ sin
顶柱柱端弯矩幅值:53
333126100.27101032.5225EI l M A KN m l
?=??=??=?
10-33 试求图示结构两质量处的最大竖向动位移,并绘制最大动力弯矩图。设m 1=m 2=m ,
3
2ml EI =θ
解:该结构有两个自由度,使用刚度法。
11122122
333
103,,77EI EI EI
k k k k k l l l =
==-=-= 11k 的求解过程:
3
8
16
l
3
1
1161525632372112163161183896l l l l l l l EI EI
δ??=????+????=
??? 1131
967EI
k l δ==
111333
9610377EI EI EI
k k k l l l =+=
+=
22k 的求解过程:
2
l
左构件3
2
11121222326l l l l EI EI
δ??=????= ???
22
3
1
6EI
k l δ=
=
22233367EI EI EI
k k k l l l
=+=
+= 将上述刚度系数,质量值及荷载幅值代入位移幅值方程,并计θ= 123
331233310340774EI EI EI
A A l l l EI EI EI A A F l
l l ???--= ????
??
???-+-= ????? 解得:33
120.032
,0.344Fl Fl A A EI EI
== 最大动力弯矩图
求解过程:
对于AB 杆件,相当于在中点作用一集中力
1196
0.0320.4397
AB F A k F F =?=?
= 对于CD 杆件,相当于在中点作用一集中力
220.3446 2.064CD F A k F F =?=?=
10-34 试说明用振型分解法求解多自由度体系动力响应的基本思想,这一方法是利用了振动体系的何种特性?
10-35 试用振型分解法计算题10-32。 解:
刚度矩阵2 - 0- 2 -0 - k k k k k k k k ????=?????? 质量矩阵2 0 00 1.5 00 0 m M m m ??
??=??????
其中6153
249610,110EI
k N m m Kg l
-=
=??=? 由刚度矩阵和质量矩阵可得:
-0.3115 0.5774 0.2639-0.5278 0 -0.6230-0.6230 -0.5774 0.5278 A ??
??=??????
111
12312.11,30.98,45.75s s s ωωω---===
()()
11
510.3115 2 0 00.31150.52780 1.5 00.52781100.62300 0 10.6230T
T
m A MA m m kg --??
????
??????==--==???
??????????--??
????
()()2252110T m A MA m kg ===? ()()3353110T m A MA m kg ===?
()()()
110.311500.5278sin 15.83sin 0.62300T
T
P t P t F A F F t t KN θθ-????
????==-=-??
??????-??
??
()(
)()
220.577400sin 00.57740T
T
P t P t F A F F t θ??
??
????===??
??????-????
()(
)()
330.263900.6230sin 18.695sin 0.52780T
T
P t P t F A F F t t KN θθ??
??
????==-=-??
????????
??
则()1t y 应满足方程
()1..
21111
P t F y y m ω+=
其稳态响应为:
()()
(
)
3
12
5215.8310sin 0.3264sin 11012.118t y t t mm θθπ-?=
=?-
同理:()20t y =
()()
(
)
3
32
5
2
18.6910sin 0.1279sin 11045.758t y t t mm θθπ-?=
=-?-
()()()()()()112233-0.3115 0.05774 0.26390.32640.1354-0.5278 0 -0.6230sin 0.0926-0.623 0.05774 0.52780.12790.2708t t t t t t y y y A y t y y θ????-??????????????????∴===-?????????????--?
????????????sin t mm θ?????? 显然最大位移
1max 2max 3max 0.13540.09260.2708y mm y mm y mm
=-=-=- 10-36 试用振型分解法计算题10-31结构作有阻尼强迫振动时,质量处的最大位移响应。已知阻尼比ξ1=ξ2=0.10。
解:
刚度矩阵1
32 40.2143 -0.321438-0.3214 0.85714 3k EI EI -??
????==???????????
? 质量矩阵0 11 0M m ??
=????
得: 1 -0.41420.4142 1A ??
=?
???
12ωω==
()()()()
()()()()()()
1211122212 1.17161.171610sin 0.4142sin 0.414210.41420sin sin 11T T
T
T P t P t T
T
P t P t m A MA m m A
MA m
F A F F t F t F A
F F t F t θθθθ====????=?=?= ? ?????-????=?=?= ? ?????
正则坐标()1t y 应满足方程:
()...
121111111
2P t F y y y m ξωω++=
其稳态响应为:()()111sin t y A t θα=-
3
10.8133A mm =
=
()1
111
12212tan tan 0.45870.43011θξωαθω--?? ?
?==-=- ?- ???
同理可得:()()222sin t y A t θα=-
3
20.1092A mm ==
2
1122222
2tan tan 0.08130.08111θξωαθω--?? ?
?=== ?- ???
于是
()()()()120.8133sin 0.43010.1092sin 0.0811t t y t mm y t mm
θθ=+=-
()()()()()()()()120.8133sin 0.43011 -0.41420.4142 10.1092sin 0.08110.8133sin 0.43010.0452sin 0.0811 0.3369sin 0.43010.1092sin 0.0811t t y t y t t t t t θθθθθθ??+??
?? ?=???? ?-????????
+--??
=??
++-????
()()()
()()()110.8133sin 0.43010.0452sin 0.0811 0.8133sin cos 0.4301sin(0.4301)cos 0.0452sin cos 0.0811sin(0.0811)cos 0.6942sin 0.3428cos 0.7742sin t y t t t t t t t t t b mm
θθθθθθθθθ=+--=?+?-?+?????????
=+=+ ()1max 0.7742t y mm ∴=(竖直方向)
()()()
()()()220.3369sin 0.43010.1092sin 0.0811 0.3369sin cos 0.4301sin(0.4301)cos 0.1092sin cos 0.0811sin(0.0811)cos 0.4150sin 0.1316cos 0.4354sin t y t t t t t t t t t b mm
θθθθθθθθθ=+--=?+?-?+?????????=+=+ ()2max 0.4354t y mm ∴=(水平方向)
10-38 试用基于能量原理的近似法求图示梁的基本频率。 (a)
(b)
l
EI
m 2
l 2
l