结构力学专题十三(多自由度体系的动力计算)
合集下载
结构力学 结构的动力计算
输出 （动力反应）
第二类问题：参数（或系统）的识别 结构 （系统）
输入 （动荷载）
输出 （动力反应）
第三类问题：荷载识别 输入 （动荷载） 结构 （系统）
输出 （动力反应）
第四类问题：控制问题
输入 （动荷载）
结构 （系统）
控制系统 （装置、能量）
输出 （动力反应）
2．结构动力计算的目的 研究结构在动荷载作用下的反应规 律，找出动荷载作用下结构的最大动内 力和最大动位移，为结构的动力可靠性 设计提供依据。
第13章
结构的动力计算
§13－1 动力计算的特点和动力自由度
一．动荷载及其分类
动荷载是指其大小、方向和作用位置随 时间变化的荷载．由于荷载随时间变化较快 ，所产生的惯性力不容忽视。因此，考虑惯 性力的影响是结构动力学的最主要特征。 静荷载只与作用位置有关，而动荷载 是坐标和时间的函数。
动荷载按其随时间的变化规律进行分类：
质量 m 在 t 时刻的位移y(t)是由此时作 用在质量上的惯性力产生的，位移方程为：
y(t ) [m(t )] y
整理，
m(t ) y 1
y (t ) 0
(a) (b)
单自由度体系： k
1
式（13－1）或（a）称为单自由度体系 自由振动运动方程（微分方程）
二.自由振动运动方程的解
由式（13－4）
y (t ) A sin(t ) A sin(t 2 ) A sin[ (t 2 ) ] y (t 2 )
y(t)是周期函数
T 2
－自振周期（固有周期） －自振频率（固有频率）
2 T
1. 结构自振周期 T 和自振频率 的各种等 价计算公式
结构动力计算多自由PPT教案
Mi
sin
(it i )
0 (i 1, 2)
第i 阶振型产生的惯性力在第j 阶振型的位移上所做的
虚功为零，也即由某振型产生的惯性力在非自身振型上不
做功
第11页/共67页
8. 5 两自由度体系的振动分析
3）振型正交性的利用
（1）可用振型的正交性来检验所求得的振型是否正 确。 例题 8-17 试检验例题 8-15所求得振型的正确性。
11
22
43 243EI
12
21
73 486EI
1
(11
12 )m
15 486
m3 EI
2
( 22
21)m
1 486
m3 EI
第8页/共67页
8. 5 两自由度体系的振动分析
1
1/ 1 5.69
EI m3
1
1
1c1
2
1/ 2 22
EI m3
2
1 1c2
第一主振 型 (正对称)
，
A2
1 B0
k21F01
k11 2m1
F02
0.1
F0h3 EI
（3）计算惯性力幅值
m1 2Y1
16
EI h3
0.075
Fh3 EI
1.2F
m2 2Y2
16
EI h3
0.1
Fh3 EI
1.6F
（4）计算内力：将荷载幅值和惯性力幅值作用在结构上， 按静力进行计算
第21页/共67页
8. 5 两自由度体系的振动分析
k21 k12 k2
(2)求位移幅值 由已知条件知:
k11 k1 k2 k22 k2
k11 2m1
结构力学中体系的计算自由度
结构力学中体系的计算自由度写在前面：①【自由度】：体系运动时所具有的独立运动方式数目，也就是体系运动时可以独立变化的几何参考数目，或者说确定体系位置所需的独立坐标数目。
如：一个点在平面中的自由度为2（两个平动）；一个刚片在平面中的自由度数为3（两个平动和一个转动）。
②【约束】：限制运动的装置称为联系（或约束），体系的自由度可因加入联系而减少，能减少一个自由度的装置称为一个联系（或约束）。
常用的联系有链杆和铰。
一根链杆为一个联系，一个单铰为两个联系，也就是相当于两根链杆的作用。
联结n 个刚片的复铰可以当做（n-1）个单铰，将减少2（n-1）个约束；联结n 个刚片的复刚结点可以当做（n-1）个单刚结点，将减少3（n-1）个约束。
在体系中加入一个联系，而并不能减少体系的自由度，这样的联系称为多余联系（联系）。
使体系成为几何不变而必须的约束称为必要联系（约束）。
方法一：结点法以结点为对象，以链杆为约束计算式：bj W -=2其中，W —计算自由度j —自由结点数（与地面紧连的不算）b —链杆数（包括支座链杆）方法二：刚片法以自由的刚片为对象，以结点和链杆为约束计算式：)23(m 3r b g W ++-=其中，W —计算自由度m —自由刚片数g —单刚结点的数目（复铰）b —单铰结点的数目（复铰）r —链杆的数目结点法例子：体系计算自由度020210=-⨯=W 说明：体系中除了与大地紧连的结点外一共有10个，链杆有20根。
刚片法例子：体系计算自由度[]116263317=⨯+⨯-⨯=W 说明：将刚片拆分成17根，刚结点4，铰结点14。
NOTE ：①刚片法中的单铰结点数和单刚结点数以及链杆数指刚片间的；②大地也是一大刚片，但大地作为参考系没有自由度，算自由刚片的时候不能算上，但是若是将支座链杆也看做是刚片的话，则需要考虑支座链杆与大地之间的单铰接点或单刚结点；③可以用计算自由度的方法分析体系的几何构造:0>W ,表明体系中缺少足够的约束，为几何常变体系；0=W ,表明体系具有成为几何不变所需要的最少联系数目。
1005多自由度体系自由振动(力学)
FI 2 (t ) FS 2 (t ) FE 2 (t ) 0
其中
FIi mi i y
FSi kij y j
j 1
2
( i 1,2)
m1 1 k11 y1 k12 y2 FE1 (t ) y
m2 2 k 21 y1 k 22 y2 FE2 (t ) y
主振型的位移幅值等于 主振型惯性力幅值作用下产 生的静力位移。
（2）振型方程
( 11 m1 2 ) A1 12 m2 A2 0 1 21 m1 A1 ( 22 m2 2 ) A2 0 1
A1=A2= 0 ？
（3）频率方程
D
11m1
y 11 12 m1 0 1 y1 Δ1P (t ) 0 m y Δ (t ) 2 y2 21 22 2 2P
m1 0 0 1 k11 y k m2 y 2 21 k12 y1 FE1 (t ) y F ( t ) k 22 2 E 2
1
2
12m2 22m2
1
21m1
0
令
1
2
2
2 (11m1 22 m2 ) (11 22 m1m2 12 21m1m2 ) 0
1 1 ( 11m1 22 m2 ) ( 11m1 22 m2 ) 2 4( 11 22 12 21 )m1m2 2 2
y2 (t ) m1 1 (t ) 21 m2 2 (t ) 22 y y
设解为 y1 (t ) A1 sin(t )
其中
FIi mi i y
FSi kij y j
j 1
2
( i 1,2)
m1 1 k11 y1 k12 y2 FE1 (t ) y
m2 2 k 21 y1 k 22 y2 FE2 (t ) y
主振型的位移幅值等于 主振型惯性力幅值作用下产 生的静力位移。
（2）振型方程
( 11 m1 2 ) A1 12 m2 A2 0 1 21 m1 A1 ( 22 m2 2 ) A2 0 1
A1=A2= 0 ？
（3）频率方程
D
11m1
y 11 12 m1 0 1 y1 Δ1P (t ) 0 m y Δ (t ) 2 y2 21 22 2 2P
m1 0 0 1 k11 y k m2 y 2 21 k12 y1 FE1 (t ) y F ( t ) k 22 2 E 2
1
2
12m2 22m2
1
21m1
0
令
1
2
2
2 (11m1 22 m2 ) (11 22 m1m2 12 21m1m2 ) 0
1 1 ( 11m1 22 m2 ) ( 11m1 22 m2 ) 2 4( 11 22 12 21 )m1m2 2 2
y2 (t ) m1 1 (t ) 21 m2 2 (t ) 22 y y
设解为 y1 (t ) A1 sin(t )
考研结构力学必看精华总结结构的动力计算
杜哈梅积分(Duhamel)
零初始条件下，单自由度体系在任意荷载下的动位移公式
若 y0 0 v0 0
则
y
y0
cos t
v0
sin t
1
m
t 0
FP
(
)
sin
t
d
第26页/共77页
（1）突加荷载
y
FP 0
m 2
(1
cos t )
yst (1 cost)
质点围绕静力平衡位置作简谐振动，动 力系数为
1， 产生共振。 但振幅不会一下增加到很大。
1
动力系数的绝对值随频率比增大而减小。
第22页/共77页
例10-3 已知：跨度l=4m，惯性矩 I=7480cm4，截面系数W=534cm3 ，弹性模 量E=2.1×105MPa。电动机重量G=35kN，转速n=500r/min，离心力FP=10kN， 竖向分力FPsint。试求梁动力系数和最大正应力。
第34页/共77页
阻尼对自振特性的影响
r 1 2
阻尼对振幅的影响
★影响小,可以忽略
ln yk ln y tk
yk1
y tk T
e tk ln etk T
ln eT
T
★振幅的对数衰减率
★阻尼越大,衰减速度越快
1 ln yk 或 2 yk1
1 ln yk 2 n ykn
2004年8月
第8页/共77页
§10-2单自由度体系的自由振动 1 振动方程的建立
刚度法 体系在惯性力作用下处于动态平衡。
myt kyt 0
柔度法 质体的动位移等于质体在惯性力作用下的静位移。
y t my t my t
【哈工大 结构动力学】SD 第10章 多自由度体系2020
➢ 所谓振型就是结构不同点(自由度)变化时的比例关系。
11
以上分析方法就是代数方程中的特征值分析，自振频率相应 于特征值，而振型即是特征向量。
得到体系的N 个自振频率和振型后，可以把振型和自振频率
分别写成矩阵的形式，
1 2 N
1 0 0
0
2
0
或
1,2,3
T
n
jqj(t)
j
2.0jjqj(t )
j
2 j
qj(t
)
T j
M
Mj
I
ug(t )
振型分解法仅需知道各振型阻尼比 ξ，不需要知道阻尼矩阵[C]
定义振型参与系数γj
j
jT M I
Mj
jT M I jT M j
基本性质
[] 1
两边同时除以振型参与系数γj ，得到：（j＝1，2，…N）
得到三个根 ：
B1 0.3515, B2 1.6066, B3 3.5420
利用关系式
Bn n2 600
可得结构的三个自振频率：
12 210.88
1 14.522
22 963.96 2 31.048 (rad / s)
32 2125.20
3 46.100
19
算例10-1
求振型 ： (K n 2 M ) n 0
7
将位移向量{u}和加速度向量{ü}代入无阻尼自由振动方程：
u 2 sin( t ) u s in ( t )
M u K u 0
( 2 M K ) sin( t ) 0
因为sin(ωt+θ)为任意的，可以消去，因此，
(K 2 M ) 0
上式是关于{φ}的N阶齐次线性方程组，表征了振型和自振
11
以上分析方法就是代数方程中的特征值分析，自振频率相应 于特征值，而振型即是特征向量。
得到体系的N 个自振频率和振型后，可以把振型和自振频率
分别写成矩阵的形式，
1 2 N
1 0 0
0
2
0
或
1,2,3
T
n
jqj(t)
j
2.0jjqj(t )
j
2 j
qj(t
)
T j
M
Mj
I
ug(t )
振型分解法仅需知道各振型阻尼比 ξ，不需要知道阻尼矩阵[C]
定义振型参与系数γj
j
jT M I
Mj
jT M I jT M j
基本性质
[] 1
两边同时除以振型参与系数γj ，得到：（j＝1，2，…N）
得到三个根 ：
B1 0.3515, B2 1.6066, B3 3.5420
利用关系式
Bn n2 600
可得结构的三个自振频率：
12 210.88
1 14.522
22 963.96 2 31.048 (rad / s)
32 2125.20
3 46.100
19
算例10-1
求振型 ： (K n 2 M ) n 0
7
将位移向量{u}和加速度向量{ü}代入无阻尼自由振动方程：
u 2 sin( t ) u s in ( t )
M u K u 0
( 2 M K ) sin( t ) 0
因为sin(ωt+θ)为任意的，可以消去，因此，
(K 2 M ) 0
上式是关于{φ}的N阶齐次线性方程组，表征了振型和自振
计算结构动力学 多自由度体系的振动
tgi=2i/i(1-i2)
(36)
将式(34)代回
{u}=ii(t){A}i ， 得
{u(t)}=[iisin(it-i)/i2]{A}i
(37)
无阻尼情况自然可以当作有阻尼情况的特例，在上
述结果中令i=0得到。
4.3 多自由度的受迫振动
4.3.3 简谐荷载受迫振动反应分析步骤
左乘{A}jT、后一式左乘{A}iT，再将两式相减，由于质 量、刚度的对称性，可得
由此可得
(i2-j2){A}jT[M]{A}i=0
(11)
{A}jT[M]{A}i=0
(12)
上式乘j2，考虑到j2[M]{A}j物理意义是第j振型对应
的惯性力幅值，因此式（12）表明第j振型对应的惯性
力在第i振型位移上不做功，反之亦然。
(e)
方程两边同时左乘{A}jT，根据正交性则有
Mj*ÿj(t)+Kj*yi(t)=0
(20)
从式(20)可得(根据单自由度自由振动结果)
yi(t)=aisin(it+ci)
(f)
代回多自由度所假设的解，即可得
{u(t)}=aisin(it+ci){A}i
(21)
5）式(21)中的待定常数ai、ci可由初始条件确定。如何
22求无阻尼自由振动的振型求无阻尼自由振动的振型aaii频率频率ii33用阻尼比用阻尼比1122和频率和频率1122求瑞利阻尼的求瑞利阻尼的00和和44求求ii振型振型参与系数振型振型参与系数iiaaiittppaaiittmmaa55求求ii振型阻尼比振型阻尼比12120066求求ii振型动力系数振型动力系数iiii222244ii22ii22121277求求ii振型相位角振型相位角iiarctg2arctg2iiii2288求求ii振型广义位移振型广义位移iittiisinsiniittiiii2299将各振型广义位移代回将各振型广义位移代回uuiittaaii则得最终则得最终结果结果uuttiisinsiniittiiii2237374444441441基本原理基本原理对动力问题设单元位移场仍表示成对动力问题设单元位移场仍表示成ddnnddee只是现在只是现在ddddxtddee设杆单元的密度为设杆单元的密度为将微段惯性力将微段惯性力aaaaddxx作为作为体积力则这一单元荷载的总虚功为体积力则这一单元荷载的总虚功为dxdx3838引入单元一致质量矩阵引入单元一致质量矩阵mmeedx39394444由式3939代入形函数并积分对质量均匀分布的平代入形函数并积分对质量均匀分布的平面弯曲单元其单元一致质量矩阵面弯曲单元其单元一致质量矩阵mmee13221561354221354221564204040作业
(结构动力学)多自由度体系运动方程
fsi ki1u1 ki2u2 kiNuN
系数kij称为刚度影响系数，简称刚度系数，物理意义是：
kij—由第j自由度的单位位移所引起的第i自由度的力
即j自由度给定一个单位位移， 而其余自由度都不动时， 所需要的力（反力）。
6.1 直接平衡法
弹性恢复力 fsi ki1u1 ki2u2 kiNuN
只要能用广义坐标给出体系总动能T和位能V的表 达式，以及确定相应于每一广义坐标的非保守 力Qi，就可以直接由Lagrange运动方程建立结构 体系的运动控制方程。
下 面 通 过 算 例 来 介 绍 如 何 应 用 Lagrange 方 程 ， 从 算例中可以看到，用Lagrange运动方程建立的运 动方程不限于线性。
自由度方法也可以得到相当好的近似解。但对于复杂
的结构体系或作用的外荷载变化复杂时，用等效的单
自由度方法得到的解可能会导致相当大的误差。这时
就必须直接采用多自由度体系分析方法解决问题，即 必须采用更多自由度来描述体系的运动状态。
第六章 多自由度体系的运动方程
建立单自由度体系运动方程的方法均可以用来建立多自 由度体系的运动方程，例如：牛顿第二定律；直接平 衡法(d’ Alember)；虚位移原理；Hamilton方程；运动 的Lagrange方程，都可用于多自由度体系。但基于矩 阵位移法的直接平衡方程和基于变分原理的Lagrange 方法应用更广泛一些。前者对于多自由度体系直接应 用动平衡的概念以矩阵的形式建立体系的运动方程， 概念直观，易于通过各个结构单元矩阵(刚度矩阵、质 量矩阵、阻尼矩阵)建立整个结构体系的相应矩阵，进 而建立体系的运动方程，便于计算机编程，在结构动 力分析的有限元程序中基本上都基于直接平衡法。而 对于一些特殊的问题，例如，大变形（位移）问题， 采用Lagrange方法可能更有效。本章将主要介绍这两 种方法。
系数kij称为刚度影响系数，简称刚度系数，物理意义是：
kij—由第j自由度的单位位移所引起的第i自由度的力
即j自由度给定一个单位位移， 而其余自由度都不动时， 所需要的力（反力）。
6.1 直接平衡法
弹性恢复力 fsi ki1u1 ki2u2 kiNuN
只要能用广义坐标给出体系总动能T和位能V的表 达式，以及确定相应于每一广义坐标的非保守 力Qi，就可以直接由Lagrange运动方程建立结构 体系的运动控制方程。
下 面 通 过 算 例 来 介 绍 如 何 应 用 Lagrange 方 程 ， 从 算例中可以看到，用Lagrange运动方程建立的运 动方程不限于线性。
自由度方法也可以得到相当好的近似解。但对于复杂
的结构体系或作用的外荷载变化复杂时，用等效的单
自由度方法得到的解可能会导致相当大的误差。这时
就必须直接采用多自由度体系分析方法解决问题，即 必须采用更多自由度来描述体系的运动状态。
第六章 多自由度体系的运动方程
建立单自由度体系运动方程的方法均可以用来建立多自 由度体系的运动方程，例如：牛顿第二定律；直接平 衡法(d’ Alember)；虚位移原理；Hamilton方程；运动 的Lagrange方程，都可用于多自由度体系。但基于矩 阵位移法的直接平衡方程和基于变分原理的Lagrange 方法应用更广泛一些。前者对于多自由度体系直接应 用动平衡的概念以矩阵的形式建立体系的运动方程， 概念直观，易于通过各个结构单元矩阵(刚度矩阵、质 量矩阵、阻尼矩阵)建立整个结构体系的相应矩阵，进 而建立体系的运动方程，便于计算机编程，在结构动 力分析的有限元程序中基本上都基于直接平衡法。而 对于一些特殊的问题，例如，大变形（位移）问题， 采用Lagrange方法可能更有效。本章将主要介绍这两 种方法。
结构力学 第十三章 结构动力学
4)
1) 平面上的一个质点
y1
W=1
y2
y1 W=2
5)
2) W=2
W=2
*弹性支座不减少动力自由度 6)
y2 y1
3)
W=2
计轴变时 W=2
***自由度数与质点个数无关，
不计轴变时 W=1
但不大于质点个数的2倍。 7)
EI
**为减少动力自由度，梁与刚架不 计轴向变形。
W=1
8) 平面上的一个刚体
FP(t )
FP
t
t
简谐荷载（按正余弦规律变化）
一般周期荷载
2）冲击荷载：短时内剧增或剧减。（如爆炸荷载）
FP
FP(t )
FP
FP
tr
t
tr
t
3）随机荷载：(非确定性荷载) 荷载在将来任一时刻的数值无法事先确定。 （如地震荷载、风荷载）
三、动力计算中体系的自由度
确定体系上全部质量位置所需独立参数的个数称为体系的振动自由度。
.. 静平衡位置 m
y(t) FI (t) my(t)
.......... .(c)
1
k
FI(t)
my(t) y 0
可得与 (b) 相同的方程
刚度法常用于层间模型，柔度法常用于质点模型。
刚度法实质： 从静力平衡角度建立运动 微分方程,思路类同于位移 法方程的建立。
柔度法实质： 从变形协调角度建立运动 微分方程,思路类同于力法 方程的建立。
经整理后,得:
1.在质量上沿位移正向加惯性力；
2.求发生位移y所需之力； 3.令该力等于体系外力和惯性力。
my(t)
24EI l3
y(t)
Fp(t)
结构动力计算教学课件PPT_OK
k12
0
k21
(k22 2m2 )
特征方程 频率方程
(k11 2m1)(k22 2m2) k12k21 0
4
(k11 2m1)(k22 2m2) k12k21 0
2
2 1,2
1 2
k11 m1
k 22 m2
1
2
k11 m1
k 22 m2
k11k22 k12k21 m1m2
最小圆频率称为第一(基本)圆频率： 第二圆频率-------
K1 F
n1 n2 nn
FMYY 0
K Fn自M由度Y体 系作K自由Y振动 的K 0 IM运动Y方程（K柔Y度法）0
将特解带入方
程整理后：
FM
1 2
IX
0
M Y
KY 0
FM
1 2
I
0
频率方程
19
FM
1
2 j
I j
0
j(1) 1
规准化主振型方程
一般的：
n个主振型向量彼此线性无关，
( j 1,2,, n)
n个自由 度体系的
依上式可求得与ωj 相对应 主振型，我们可唯一地确 振型方程
定主振型的形状，但不能唯一地确定它的振幅。
N自由度体系有n个主振型，若体系为对称形式，则这些主振型
分为对称及反对称形式两类。
17
主振型的规准化：
为了使主振型的振幅也具有确定值，需另外补充条件， 由此得到的主振型叫规准化主振型。
则系数行列式为零：
K 2 M 0
n个自由度体系 的频率方程
n个频率（按数值大小从小到大排列）： ω1，ω2，---，ωn
令：Xj 表示与频率ωj相对应的主振型向量：
结构力学 体系的计算自由度
O13 行吗？
瞬变体系
它可 变吗？
2
有
几
个 单
3
铰？
1
讨论
2
3 1
体系W
等于多少？ 可变吗？
W=0,体系 是否一定
几何不变呢？
W=3 ×9-(2×12+3)=0
除去约束后，体系的自由度将增 加，这类约束称为必要约束。
因为除去图中任 意一根杆，体系 都将有一个自由 度，所以图中所 有的杆都是必要 的约束。
除去约束后，体系的自由度并不 改变，这类约束称为多余约束。
在m=2的情况下，刚片间没有铰结点，h=0
W＝3×2－（3×3＋7）＝－10
解法一： 所有结点都是铰结点，j＝16
包括支座在内共有连杆31根
W＝2×16－31＝1
解法二： 图示三角形视为刚片，m＝8 刚片间单铰h＝8，刚结点没有，g＝0 包括支座在内共有连杆7根
W＝3×8－(2×8＋7)＝1
例1：计算图示体系的自由度
瞬 变 体 系
常变体系
小结
几何不变体系 可作为结构
体系
几何可变体系 不可作结构
无多余联系
静定结构
有多余联系
超静定结构
常变
瞬变
分析示例 加、减二元体 无多几何不变
瞬变体系 去支座后再分析
加、减 二元体
无多几何不变
找虚铰 无多几何不变
找 刚Ⅰ 片 、 O23 找 虚 铰
无多几何不变 O12
Ⅱ Ⅲ
在m=11的情况下，刚片间没有铰结点，h=0
W＝3×11－（3×12＋7） ＝－10
解法二：
将ABCDEGHI、FGHIJ看
作刚片，m＝2
A
瞬变体系
它可 变吗？
2
有
几
个 单
3
铰？
1
讨论
2
3 1
体系W
等于多少？ 可变吗？
W=0,体系 是否一定
几何不变呢？
W=3 ×9-(2×12+3)=0
除去约束后，体系的自由度将增 加，这类约束称为必要约束。
因为除去图中任 意一根杆，体系 都将有一个自由 度，所以图中所 有的杆都是必要 的约束。
除去约束后，体系的自由度并不 改变，这类约束称为多余约束。
在m=2的情况下，刚片间没有铰结点，h=0
W＝3×2－（3×3＋7）＝－10
解法一： 所有结点都是铰结点，j＝16
包括支座在内共有连杆31根
W＝2×16－31＝1
解法二： 图示三角形视为刚片，m＝8 刚片间单铰h＝8，刚结点没有，g＝0 包括支座在内共有连杆7根
W＝3×8－(2×8＋7)＝1
例1：计算图示体系的自由度
瞬 变 体 系
常变体系
小结
几何不变体系 可作为结构
体系
几何可变体系 不可作结构
无多余联系
静定结构
有多余联系
超静定结构
常变
瞬变
分析示例 加、减二元体 无多几何不变
瞬变体系 去支座后再分析
加、减 二元体
无多几何不变
找虚铰 无多几何不变
找 刚Ⅰ 片 、 O23 找 虚 铰
无多几何不变 O12
Ⅱ Ⅲ
在m=11的情况下，刚片间没有铰结点，h=0
W＝3×11－（3×12＋7） ＝－10
解法二：
将ABCDEGHI、FGHIJ看
作刚片，m＝2
A
【结构力学】第十三章 结构动力计算(1)汇总
11
1 EI
lh 2
2h 3
lh 2 3EI
w 1 3EI m11 mlh2
15
例6
k11
k11
解：求 k
3EI
3EI
m
k l3
1
k
EI
k11 k l 3
l
w k11 3EI l3 k
m
m
•对于静定结构一般计算柔度系数方便。
•如果让振动体系沿振动方向发生单位位移时，所有刚节点
都不能发生转动（如横梁刚度为无穷大的刚架)计算刚度系数方
2
1
m
l/2
11
1 EI
(
1 2
l 2
l 2
2 3
l 2
1 2
l 2
l
l 2
2 3
)
l3 8EI
w
1
m11
8EI ml3
14
例5 1
k Am
h
I=∞ EI
B
θ
1
1
h
解法1：求 k θ=1/h
MBA=kh
=
MBC
3
EI l
3EI lh
C
k 3EI lh2
w
k11 m
3EI mh2l
解法2：求 δ
2n
6023.14500
6052.3
1 S
52.3/57.4=0.91
2）求动力系数β
1
1
5.88
1 2 w 2 152.32 5379.742 1.35
max
Ql 4W
Pl 4W
(Q P)l
4W
175.6MPa
149.2
结构动力学多自由度
k = k − kG
弹性特性
刚度的定义： 刚度的定义：
表示一个自由度发生相应单位位移而其他节点不动时在结构中所 产生的的力。 产生的的力。
弹性特性
柔度的定义： 柔度的定义：
~ ~ ~ 则任意荷载组合下： 则任意荷载组合下： v i = f i 1 p1 + f i 2 p2 + L + f iN p N ~ ~ ~ 用矩阵表示： 用矩阵表示： v1 f 11 L f 11 L f 11 p1 M M M M M ~ f 11 L ~11 L ~11 pi vi = f f M M M M M v N ~ L ~ L ~ p N f11 f11 f 11 ~ v = fp
多自由度体系的振动分析
运动方程
多自由度体系的动力平衡方程： 多自由度体系的动力平衡方程：
f I + f D + f S = p(t )
即：
&& & m v + cv + kv = p(t )
&& & m v + cv + kv − kG v = p(t )
考虑几何刚度： 考虑几何刚度： 或
&& & m v + cv + k v = p(t )
当结构的全部有限自由度的刚度系数均以求得后， 当结构的全部有限自由度的刚度系数均以求得后，只要适当 地叠加单元的刚度系数，就能得到整个结构的刚度， 地叠加单元的刚度系数，就能得到整个结构的刚度，这就叫 做直接刚度法。 做直接刚度法。 结构的任何一个刚度系数，都能通过与这些节点相连的单元， 结构的任何一个刚度系数，都能通过与这些节点相连的单元， 所对应的刚度系数叠加求得。 所对应的刚度系数叠加求得。
结构动力学多自由度
pbT
~ fpa
paT ~fpb paT ~fpb T pbT ~f T pa pbT ~fpa
故 ~f 、 k 均为对称矩阵。
单元刚度矩阵
单元刚度系数表示由单位节点位移所引起的节点力。
单元刚度系数由虚位移法求得。
例如，课本P106图11-5所示简支梁中，令a端发生单位转角， 并给该处一竖向虚位移，零外力所做的功，等于内力所做的 功。
表示一个自由度发生相应单位位移而其他节点不动时在结构中所 产生的的力。
弹性特性
柔度的定义：
~ fij —在j坐标施加单位荷载而引起的i坐标的挠度。
则任意荷载组合下： vi ~fi1 p1 ~fi2 p2 ~fiN pN
用矩阵表示：
v1
vi
v N
略去阻尼矩阵和施加的荷载向量的影响： mv kv 0
假定以上多自由度体系的振动是简谐振动：
v(t) vˆ sin(t )
vˆ 表示体系的形状，不随时间变化。
v 2vˆ sin(t ) 2v 2mvˆ sin(t ) kvˆ sin(t ) 0
k 2m vˆ 0
无阻尼自由振动—振动频率分析
k 2m vˆ 0
即： k 2m 0
上式的N个根，表述体系可能存在的N个振型的频率。
1
2
3
WE va pa v1 k13
Lபைடு நூலகம்
WI v1 0 EI ( x) 1''( x) 3''( x)dx
结构动力学之多自由度体系的振动问题ppt课件
1 536EI
448 (1 536)2
m1m2l 6 (EI )2
0
解得
21
23l3 (m1 m2 2 1 536EI
)
529(m1 m2 )2l6 41 5362 (EI )2
448m1m2l 6 1 5362 (EI )2
从而得第一和第二阶自振频率
1
1
1
2
1
2
为了确定第一阶振型，可将1代入平衡方程。
其展开式是关于λ的n次代数方程,先求出λi再求 出频率ωi
柔度法
(11m1 )
12m2
...
21m1 ( 22m2 ) ...
...
...
...
1n mn 2nmn 0
...
n1m1
n2m2 ... ( nnmn )
将λi代入 ( [δ] [M] － λi [I ] ){Y(i)}={0} 可求出n个主振型。
多个自由度体系的自由振动
结构在受迫振动时的动力响应与结构的动力特性 密切相关；另外，当用振型叠加法计算任意干扰力 作用下结构的动力响应时，往往要用到自由振动的 频率(frequency)和振型(mode)。
为此，要需要首先分析自由振动。
自振频率和振型的计算
m1
m2
mi
mn
y1(t) y2(t)
yi(t)
刚度法
其中最小的频率1 称为最低自振频率，或称
基本频率。 通常将上述每一个频率所对应的振动都称为
主振动，对应于每一个主振动的形状称为主振 型。
1）如果各质体的初速度为零，而初位移和某 一振型成比例，然后任其自然，则系统就按 这个振型作简谐自由振动，此解答就相应于 该振动的一组特解；
448 (1 536)2
m1m2l 6 (EI )2
0
解得
21
23l3 (m1 m2 2 1 536EI
)
529(m1 m2 )2l6 41 5362 (EI )2
448m1m2l 6 1 5362 (EI )2
从而得第一和第二阶自振频率
1
1
1
2
1
2
为了确定第一阶振型，可将1代入平衡方程。
其展开式是关于λ的n次代数方程,先求出λi再求 出频率ωi
柔度法
(11m1 )
12m2
...
21m1 ( 22m2 ) ...
...
...
...
1n mn 2nmn 0
...
n1m1
n2m2 ... ( nnmn )
将λi代入 ( [δ] [M] － λi [I ] ){Y(i)}={0} 可求出n个主振型。
多个自由度体系的自由振动
结构在受迫振动时的动力响应与结构的动力特性 密切相关；另外，当用振型叠加法计算任意干扰力 作用下结构的动力响应时，往往要用到自由振动的 频率(frequency)和振型(mode)。
为此，要需要首先分析自由振动。
自振频率和振型的计算
m1
m2
mi
mn
y1(t) y2(t)
yi(t)
刚度法
其中最小的频率1 称为最低自振频率，或称
基本频率。 通常将上述每一个频率所对应的振动都称为
主振动，对应于每一个主振动的形状称为主振 型。
1）如果各质体的初速度为零，而初位移和某 一振型成比例，然后任其自然，则系统就按 这个振型作简谐自由振动，此解答就相应于 该振动的一组特解；
结构力学课件 动力计算
M AC 3i A 3i 3i(1 A ) M AB i i i 3 A 3 3 ( A 1) 2 2 2
l
l
M 1 2
M AB 6i M AB , FQBA 1 2 2l 12i k k11 r11l m 1 2
弹簧反力 FRC 5m 2a,
2k FRC 5m 2 a 2 , 弹簧变形 2a 5m k k
结束放映 上一页 下一页
郑州大学土木工程学院 樊友景
例3：求图示体系的 自振频率。 ⑶间接应用公式
A
l/2
B m B θ
EI＝∞ l/2
C k C
l/2
D m D
A
设梁绕A点转角为θ, 2 2 2 l 3l 5ml •则梁绕A点转动的惯性 JA m m 2 2 2 为其对A点的转动惯量. A k B C D 11 •与其相应的刚度为梁 l θ=1 绕A点转动刚度k11 。 kl k11 kl 2 •用转动刚度k11代替刚 度系数k，转动惯量JA k11 k kl 2 2k 2 代替质量m代入公式 m J 5ml / 2 5m
1 192 EI m 5ml 3 192 90 105 1 134.16 s 5 300 43
结束放映
EI
F sin t
m
1
l/2
l/2
k
δ
k 1/2
上一页
下一页
郑州大学土木工程学院 樊友景
例5:已知图示梁中质量 m=300kg，EI=90×105N· m2， l=4m，k=48EI/l3，F=20kN，θ=80s-1。
结束放映
上一页
下一页
l
l
l
M 1 2
M AB 6i M AB , FQBA 1 2 2l 12i k k11 r11l m 1 2
弹簧反力 FRC 5m 2a,
2k FRC 5m 2 a 2 , 弹簧变形 2a 5m k k
结束放映 上一页 下一页
郑州大学土木工程学院 樊友景
例3：求图示体系的 自振频率。 ⑶间接应用公式
A
l/2
B m B θ
EI＝∞ l/2
C k C
l/2
D m D
A
设梁绕A点转角为θ, 2 2 2 l 3l 5ml •则梁绕A点转动的惯性 JA m m 2 2 2 为其对A点的转动惯量. A k B C D 11 •与其相应的刚度为梁 l θ=1 绕A点转动刚度k11 。 kl k11 kl 2 •用转动刚度k11代替刚 度系数k，转动惯量JA k11 k kl 2 2k 2 代替质量m代入公式 m J 5ml / 2 5m
1 192 EI m 5ml 3 192 90 105 1 134.16 s 5 300 43
结束放映
EI
F sin t
m
1
l/2
l/2
k
δ
k 1/2
上一页
下一页
郑州大学土木工程学院 樊友景
例5:已知图示梁中质量 m=300kg，EI=90×105N· m2， l=4m，k=48EI/l3，F=20kN，θ=80s-1。
结束放映
上一页
下一页
l
方案介绍-结构力学自由度计算13
y
y
x
y
o
o
x
x
A
1①
B
2 ②3 E
C
D
解: m 3, n 2, r 4
w 3m (2n r)
3 3 (2 2 4)
1
A 1
4 B
2D
1
C
E2
J
M
2
F
H
G
2
2
F
G
H
2
A
3
2
1
2
3
D
C2
B
1 E
A J
B 2
H 1C 2 G1源自DE2F 1
K
A
J
B2
三杆交于一点
F D B
A
C
E
刚片1
三杆平行不等长
A
C
B
三铰共线
常变体系——发生大位移的体系。
刚片2
B
D
F
A
C
E
刚片1
A
K
L
N H
B M
G
C
D E F
That's why I want to stay here. 那就是我想待在这里的原因。(从句作表 语) (四)宾语 表示谓语动词动作对象的成分叫宾语,即宾语是动作的承受者。英语 中,及物动词(或相当于及物动词的动词短语)、介词需带宾语。可以充 当宾语的有名词、代词、动词不定式、动名词、名词化的形容词以及 从句(宾语从句)等。 I am reading a book. 我在看书。(名词作动词的宾语) I'm going to Beijing with my father.我计划和我父亲去北京。(名词作介 词的宾语) Yesterday, Tom's mother looked after him at home.昨天汤姆的妈妈在家
结构力学专题
2
5 k 3 m
2
P0 sin t P0 sin t
k m
P0 sin t
P0 sin t
P0 sin t
1
1
P0 K
?
?
?
一班课后练习
例4：用振型叠加法求图示结构的稳态振幅。
k1 3103 kN / m ; k2 2 103 kN / m; m1 m2 m 10200 kg;
1
讨论： （a）在简谐荷载作用下，各质点仍做简谐振动。
(b) 当 0 时，
同单自由度体系 同单自由度体系 发生共振。
(c) 当 时，
(d ) 当 i 时， i 1,2,n
２、动内力计算——动静法。（同单自由度体系）
例1：求图示结构的动位移幅值和动内力幅值。
k1 k，k 2 2k， m1 m，m2 2m；
?
y1st ?
1 ?
y2 st ?
2 ?
例5：对下图示体系，试证明： 当 y1 0 y2 0 ，y 1 0 y 2 0 时， 体系只按第二主振型振动。 提示： m1=m
l /3
m2=m
l /3
l /3
y(t ) 1(t )1 2 (t )2
m1
m2
l /3
y1 EI y 2
l /3
l /3
P
I1
I1 0.2936P I 2 0.2689P
A1
A2
I2
A1 0.025 Pl3 / EI A2 0.023 Pl3 / EI
1.2936 P
0.2689P
0.3173Pl
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
FP1
m1
l
EI
l
FP 2
m2
l
二、任意荷载作用*
运动方程： M y(t) Ky(t) FP (t) (a)
1、主振型矩阵
1 2 n
2、广义质量、广义刚度
} M * T M 对角阵
K* T K
3、正则坐标
y(t) (t)
(b)
M y(t) Ky(t) FP(t) (a)
４、振型迭加法分析强迫振动
例1：求图示结构的动位移幅值和动内力幅值。
k1 k，k2 2k，
m1
m1 m，m2 2m；
P0 sin t
EI1
k1 m2
h
已知：
2
k m
EI1
k2
h
A
P0 k
1 0
1
1
I
F
0P0
P0
P0
P0 k
动位移幅值图
动荷载图（虚拟）
例2：求图示结构的动位移幅值和动内力幅值。
已知：
i
(t
)
i
(0)
cos
it
i (0) i
sin
it
(i 1, 2)
l
0E.I041
P0 L3 EI
sinP0 stin
m
t
EI
从以上例题的计算中可看出，一般情况下 １l 〉２ 〉l〉ｎ
故在振型迭加法中,一般是前几阶振型起主要作用。
思考：用振型叠加法求例1所示结构的位移幅值。
2
k m
2
1 3
k m
2 5 k 3m
2
k m
P0 sin t
P0 sin t
P0 sin t
2
4 EI ml3
A
1
1
P0 L3 8EI
1
L3 P0 8EI
m
l EI
EI
l
I
F
PP00/
/ 2
2
0.5P0
1.5P0
P0 sin t
m
l
1
动位移幅值图
动荷载图（虚拟）
例3：求图示体系的稳态振幅、动弯矩幅值图。
已知：m1 m2 m, 3.415
EI ml 3
I1 0.2936 P I2 0.2689 P A1 0.025 Pl 3 / EI
2 (t) 0.534 10 3 (1 cos2t)
FP 8kN m2
EI1
k2 m1
EI1
k1
yy12((tt))
1 2
1
(t
)
1 1/
2
2
(t
)
y1 (2.67 3.20 cos1t 0.53cos 2t) 103 y2 (6.67 6.40 cos1t 0.267 cos2t) 103
P0 sin t
P0 sin t
1
1
?
?
?
P0
K
例4：用振型叠加法求图示结构的稳态振幅。
k1 3103 kN / m; k2 2103 kN / m; m1 m2 m 10200 kg;
1 9.899 1/ s;2 24.25 1/ s
1
1 2
2
1 1/
2
1(t) 3.20 10 3 (cos1t)
将 y(t) (t) y(t) (t)
｝ 代入（ａ）式，得：
M * (t) K*(t) FP (t)*
(c)
展开任一式得：
i (t) i2i (t)
FPi* (t )
M
* i
形式同单自由度体系强迫振动
y(t) 1(t)1 2 (t)2
例1：用振型叠加法求图示结构的动力反映。
m1
已求出 1T 1 0.5,
1、位移幅值求法
y(t) A sin t
A K 2M 1F A I 2 M 1 F
讨论：
（a）在简谐荷载作用下，各质点仍做简谐振动。
(b) 当 0时，
同单自由度体系
(c) 当 时，
同单自由度体系
(d) 当 i 时， i 1,2,n 发生共振。
２、动内力计算——动静法。（同单自由度体系）
P sin t
m1
m2
l / 3 y1 lE/I3y2 l / 3
A2 0.023 Pl3 / EI
P
A1
I1
A2
I2
作业： 10—23
补充题：求图示结构的动位移幅值和动内力幅值。
其中:
FP1(t) P0 sin t
2
55.25EI ml3
FP2 (t) 2P0 sin t
m1 m2 m
2 T 1 1 P0 sin t
EI1
k1 m2
h
1(t)
2 3
P0 K
sin
t,
2 (t )
1 3
P0 k
sin
t
EI1
k2
h
例2：用振型叠加法求图示结构动位移。
已求出 1T 1 0.278, 2T 1 3.61m
1(t)
0.084
P0 L3 EI
sin
t,
2
(t
)
第十章 结构动力计算
§10-3 多自由度体系的动力计算
§10.3.3 多自由度体系的强迫振动
一、简谐荷载作用 ——直接解法
运动方程： M y(t) Ky(t) FP(t)
(a)
或 M y(t)y(t) FP(t)
(b)
FP (t) F sin t FT F1, F2,, Fn
一、简谐荷载作用——直接解法
? A1
A2
y1st ?
1 ?
y2st ? 2 ?
例5：对下图示体系，试证明：
当 y1 0 y2 0，y1 0 y2 0 时，
体系只按第二主振型振动。
提示：
m1=m m2=m
y(t)1(t)12(t)2 l /3 l /3 l /3
用正则坐标表示的运动方程为：
i (t) i2i (t) 0 (i 1, 2)