第十三章 振动分析的矩阵迭代法
线性振动理论和振动近似解法简略史 正文
线性振动理论和振动近似解法简略史摘要：读史使人明智，本文意在对线性振动理论和工程振动近似解法的发展做简要明了的阐述，其中线性振动理论史以科学家对具体模型的解答为路线，依次阐述：单摆、弦线、梁、膜、板壳、三维弹性体理论、激励响应和强迫振动理论。
工程近似解法以时间为顺序依次阐述各近似解法，依次简要阐述：邓克莱法、逐步近似法、阵型叠加法，传递矩阵法、瑞立法、里茨法、有限元法。
部分近似解法做了较为详细的解释。
关键词线性振动近似解法简略史1线性振动理论1.1单摆单摆是最早引起人们注意的振动之一，真正对单摆的研究要追溯到16世纪，早在1581年，伽利略发现了摆的等时性，之后科学家对单摆的研究主要就是计算摆的周期，当然也包括伽利略本人。
伽利略在1638年用落体公式推得摆动周期正比于摆长与重力加速度比的平方根,还从能量的角度讨论摆的周期，但始终没得到正确的比例系数。
结束摆周期的计算是在17世纪中后叶，惠更斯利用几何方法，得到摆振动周期的正确公式。
1678年牛顿在其划时代的《自然哲学的数学原理》中建立运动变化与受力的关系，使振动问题的动力学研究成为可能,假设了介质阻力与速度及速度平方成正比,形成阻尼概念的雏形，在1728年欧拉考察了摆在有阻尼介质中的运动建立并求解了相应的二阶常微分方程，至此单摆在无阻尼和有阻尼的条件下的周期计算基本结束，后期对摆的研究主要集中在摆的大幅振动和其具有的非线性特征。
1638年伽利略摆动周期正比于摆长与］重力加速度比的平方根，即：T二1LJ(1673年伽利略利用几何方法得到单摆振动周期的正确公式\准确解：T二4j(l/g)*K(sin(a/2))广1728年欧拉建立并求解了摆在有阻尼、介质中的运动相应的二阶常微分方程I)图一单摆周期的发现及求解简略图1.2弦线在振动力学研究兴起之前,有两个典型的振动问题引起注意,一个是单摆摆动，另一个就是弦线振动。
弦线振动是无穷多自由度连续系统的振动,单摆摆动是单自由度离散系统的振动,振幅不大时都可认为是线性的。
结构动力学振动分析的矩阵迭代法
（13-28）
在试探向量中消除第一振型分量的方便方法是应用滤型 矩阵 S1
~
其中
（0） 2
（0） 2
1 T (0 ) (0 ) - φ1 φ 1 m 2 S1 2 M1
（13-29）
1 T S1 Ι - φ 1φ 1 m M1
（13-30）
在这种情况下，式（13-5）可以写成
1 （1） （0） D ~ ~2 2 2 2
克拉夫书P205 例题E13-1 通过计算图E11-1的三层建筑框架的第一振型和频率来说 明矩阵迭代法。虽然用例题E11-1中导得的刚度矩阵求逆可以很容易 求得该结构的柔度矩阵，但是为了说明柔度矩阵的求法，这里对每 一个自由度相继施加单位荷载进行推导。根据定义，由这些单位荷 载所产生的位移表示柔度影响系数。
（0） 1
(0 )
φY
（0） 1 1
φ Y
（0） 2 2
φ Y
（0） 3 3
（13-14）
第一振动频率的振动形状所对应的惯性力为
(0 ) 2 f(0) 12m1 1 m(0)
2 2 2 ( / ) 记 1 ，展开得 n 1 n
（13-15）
（ 1） 1
D
（0） 1
（13-7）
（ 1） 来进行规格 1
该向量除以向量中最大的元素max（ 化，得到改进的迭代向量：
）
（1） 1
max （ ）ຫໍສະໝຸດ （1） 1 （1） 1（13-8）
设 k 为向量中任一自由度，频率近似值为：
2 1
(1)
(0 ) k1 (1) k1
机械振动教学大纲
《机械振动》教学大纲一、课程基本信息二、课程目的和任务《机械振动》是理论与应用力学等力学类本科专业必修的专业课程，同时也是机械、土建等工程学科本科和研究生培养的一门专业基础课程。
《机械振动》是一门系统地研究自然界和工程技术领域中振动现象的产生机理、运动规律、描述和控制方法的科学。
本课程教学应立足于加强学生的振动力学基础理论素养和相关基本技能培养，并着眼于拓宽学生的相关工程背景，提高科学建模能力，为今后学生能够创造性的从事相关理论研究或工程技术实践奠定必要的基础。
三、本课程与其它课程的关系本课程学习所需的主要选修课程为微分方程、矩阵理论、概率与统计、理论力学、材料力学等一系列数学、力学基础课程。
本课程教学应紧密结合相关的实验力学教学共同完成。
通过本课程的学习，为学生完成相关毕业设计课题奠定必备的基础。
四、教学内容、重点、教学进度、学时分配第一章绪论（2学时）1、主要内容机械振动的概念、振动理论研究体系、振动系统分类、简谐振动以及振动发展历史概述（选）2、本章重点机械振动的概念，振动理论研究体系，简谐振动3、本章难点振动系统分类4、教学要求从工程实践方面介绍广泛存在的振动现象，概括其特点和共同性，由此给出机械振动的科学概念。
指出振动理论的研究体系，分类的方法及振动力学的发展历史与现状，特别是指出振动力学在工程中的应用前景和应用价值；介绍相关参考书，提示学生在今后的学习中，从全书观点逐步理解分类的系统性。
第二章单自由度系统的自由振动（10学时）1、主要内容单自由度系统的无阻尼自由振动、等效质量与等效刚度、等效黏性阻尼和有阻尼自由振动。
2、本章重点建立振动微分方程、固有频率和振型、阻尼比、幅频和相频曲线与共振。
3、本章难点建立微分方程、固有频率、振幅减缩率和阻尼比。
4、教学要求介绍单自由度振动系统的工程实际背景，给出描述这一自然现象的力学模型，通过牛顿法和拉氏法建立数学模型及其简化理由和适用条件。
给出固有频率、阻尼特性及它们在自由振动中的物理意义，着重讲解幅频特性、相频特性曲线的物理意义及其在工程设计、控制中的重要作用。
第十三章 振动分析的矩阵迭代法
(13-27)
M1
不包含第一振型的试探形状为
0 0 20 v v 1Y 2 1
(13-28)
§13.4 高阶振型分析 在试探向量中消除第一振型分量的方便方法是 应用淘汰矩阵S1
2 v 0 v 0 2 1 11T mv 0 S1v 0 2 2 M1
ˆ ˆ v=n 2Dvn
v1 =Dv1
1 0
(13-5)
(13-5a)
§13.2 基本振型分析—Stodola法
高等结构动力学
1
2 1
1 =v1 v
1 1
(13-6*)
考虑任一点k的位移
1
12
1 vk0 =vk1 1
(13-7*)
2 1
0 v
k1 1 vk1
该方程最后一个等号成立,第一振型的系数远大于其它的系数.1 3
2s
(13-23)
最终结果可视为
v1
s
1Y1(0) 1 (0) max(1Y1 )
(13-24)
§13.3 收敛性的证明
高等结构动力学
(13-32)
D2 =DS1
(13-33)
此时可用下式计算频率
2 2 (1)
(1) T (v 2 ) mv(0) 2 (1) (1) v 2 mv 2
(13-34)
v 2 D 2 v(0) 2
§13.4 高阶振型分析 第三和更高振型的分析
高等结构动力学
净化了的第三振型的形状为
(0) 3 v (0) Y (0) Y (0) v 3 1 1 2 2
N
k
i
传递矩阵法
传递矩阵法适用于计算链状结构的固有频率和主振型
多个圆盘的扭振，连续梁，气轮机和发电机的转轴系统
特征：可简化为无质量的梁上带有若干个集中质量的横向振动 系统 特点：将链状结构划分为一系列单元，每对相邻单元之间的传 递矩阵的阶数等于单元的运动微分方程的阶数，因此传递矩阵 法对全系统的计算分解为阶数很低的各个单元的计算，然后加 以综合，从而大大减少计算工作量。 （1）轴盘扭转振动系统 （2）梁的横向弯曲振动系统
第 i 个梁段左端与第i-1梁段右端状态变量的传递关系：
1 li li /(2 Ei I i ) li /(6 Ei I i ) y y 2 0 1 li /( Ei I i ) li /(2 Ei I i ) 0 0 M M 1 li F 0 1 s i 0 0 Fs i 1
Z0R H0P Z0L , Z1LA H1f Z0R , Z1RA H1P Z1LA , Z2R H2Z1RA , Z3R H3Z2R
• 这时需要考虑分支系统对齿轮A的影响，重新推导。
• 假定齿轮A、B的转动惯量可以忽略不计，其传动比为 n，由于是外啮合，则其转角关系为： 1B n1A R • 扭矩关系为： M1R A nM1B • 分支系统的传递关系为： Z R H Z R
M 1 A
R
1 2 2 n I4 2 1 I4 k4
0 L 1 M 1 A
其中
1 2 2 n I4 P H1 2 I 4 1 k4
0 1
L
I1
1 R R 1 1 1 k 98 2 M I M M 2 2 2 2 I 1 1 1 9.8 1 9.8 / 98 k 1 R R 1 1 R 1 k 98 2 M M 2 2 3 2 I 1 I M 2 19.6 1 19.6 / 196 2 k
振动分析矩阵迭代法ppt课件
精选编辑ppt
10
此课件下载可自行编辑修改，供参考！ 感谢您的支持，我们努力做得更好！
精选编辑ppt
11
(s)
max(v1 )
v(4) 11
v(5) 11
210.77
12
(1)
(v1
)T
mv1(0)
(1)
(v1
)T
(1)
mv1
(22.5, 24.5,18)(1,1,1)T 3600 (22.5, 24.5,18)(22.5,16.5,9)
218
精选编辑ppt
6
二、高阶阵型分析
第二振型以及更高振型的分析可以通过引入淘汰矩阵的方法 实现更高一阶收敛，依次就能求得相关振型。
• 计算自由振动的频率常常是利用柔度法的公 式进行迭代。 由用公迭式代法(m反复12进)X行0迭，代可，以就得可到以：得m到X 体1系2 X的 最低频率和相应的阵型。
精选编辑ppt
2
迭代收敛性的证明：
Stodola法和Holzer法 Stodola法是先假定初始振型并不断迭代调整
获得一个实际振型的适当近似形式为止， 然后从运动方程确定振动频率。 Holzer法，恰好与之相反，先假定振动频率 并迭代调整直到满足边界条件为止，频率 和振型就在满足边界条件的过程中确定了。 Stodola法（迭代分析法）
(o) 2
T 1
m
1Y
(o 1
)
T 1
m
1Y
2
(o
)
由阵型的正交性可知：
Y1(o )
T 1
m
(o) 2
M1
消去
~
( 2
o
)
(o) 2
Y 1( o ) 1
线性振动近似计算方法 振动力学课件
1
i 1
i 1
以外，第二阶以上的固有频率 2 ,
远,小n 于 可近似1 地
1 12
忽略，导出基频近似公式。
导出基频近似公式
1
12
n
1
2
=
n
i1 i i1
fii mi
此公式算出的基频必小于实际基频，成为实际基频的下限。
例题1： 用邓克利法估计系统的基频下限。
解：
1 0 0
M m 0 1 0
a j
aT Ka
2
aT a j
Ka ,
a j
aT Ma
2
aT a j
Ma
其中
aT a j
Ka
eTj
j为1r,阶,n单 位阵地第阶，得到的r个方程综合为
K 2M a 0 又归结为本征值问题。
注意：
与原系统的本征值问题比较，矩阵的阶数r小于原系统的 阶数n。所以说，里茨法实质上起着使坐标缩并的作用。 缩并后的本征值问题计算与原系统类似，可导出r个固有频 率和r个模态，此缩并后的模态同样具有正交性。由于满足 瑞利商 的驻值条件，用里茨法计算模态比瑞利法更为合理， 但毕竟不是真实的模态，所得出的固有频率仍高于真实值。
必须满足i 方程。
Di ii
其中 i为第 阶i 固有频率平方的倒数。
任意选定系统的一个假设模态 ，它一般不是真实模
态，但总能表示为真实模态的线性组合式
n
ai i
i 1
左乘D矩
阵D
n i 1
ai D i
n
iai i
i 1
1
a1
1
n i2
ai
i 1
i
再左乘一次D矩阵，
用于振动分析的有限元方法概要
对于基本的一维元素进行有限元分析，能得到质量矩 阵与刚度矩阵和所需的力矢量，对于二维三维，元素矩阵 会转换成相关的更高维的空间。使用一致的和集中质量矩 阵的有限元方程并结合边界条件能为复杂系统提供解释。
最后,使用序得到在轴向载荷下的指定节点 位移,固有振动频率和特征值分析。
本章目的
用于振动分析的有限元方法
指导老师：陈益
报告人：成志斌 韩宗彪
何瑜 宁鹏
内容
有限元介绍 单个元素的运动方程 单个元素的质量矩阵、刚度矩阵、力矢量及其转化 整个系统的运动方程 整个系统的边界条件的加载及质量矩阵
MATLAB实例及总结
有限元法简介
有限元法是一种可用于精确地(但近似)解决许多复杂的 振动问题的数值方法。
把上式写成矩阵形式： （13）
所以等效节点力可以写成：
t k ut f t mu
（14）
如下图所示，一个梁单元也是有两个节点，但是有四个自由度，每个节 点处，有两种位移形式，一个是线位移，即挠度，一种是角位移。
二 梁单元
图中， f1 t , f 3 t 是力，
二 单元矩阵的坐标变换
局部坐标系：以各个单元本身的轴线为基准所设立的坐标系。便于计算节点位移。 缺点：如果整个系统里各个单元取向各异，各个节点位移方向不一致，如下图。如 何使汇交于一个节点的各个杆件的节点位移真正相等? 解决方法：进行坐标变换
结构动力学7
{ψ } [ K ]{ψ } ρ (ψ ) = {ψ }T [ M ]{ψ }
T
若假设振型{ψ}接近结构的基本振型，则Rayleigh熵为：
{ψ }T [ K ]{ψ } ρ (ψ ) = ≈ ω12 {ψ }T [ M ]{ψ }
7.1 Rayleigh法 ——近似的证明
假设振型可表示为结构固有振型的线性组合
{ } = ∑{φ}iYi = [φ ]{Y } ψ
i =1
N
设振型为正交归一化振型，则Rayleigh熵可表示为
{Y }T [φ ]T [K ][φ ]{Y } = N Y 2ω 2 / N Y 2 ρ (ψ ) = T T ∑i i ∑i {Y } [φ ] [M ][φ ]{Y } i =1 i =1
7.2 Rayleigh-Ritz法 虽然用Rayleigh法能获得较为满意的结构基频的 近似解，但在动力分析中，为得到足够精确的 结果，常常需要使用一阶以上的振型和频率。 Rayleigh法的Ritz扩展可以求得结构前若干阶固 有频率的近似值，同时还可以获得相应阶数的 振型。 Rayleigh-Ritz法首先通过假设一组振型，要求其 Rayleigh熵取极值，从而获得一低阶的特征方程 组，由此低阶方程组可以获得体系的一组自振 频率和自振振型。
7.1 Rayleigh法
结构最大动能：
T=
1 2 2 Z ω {ψ }T [ M ]{ψ } cos2 ωt 2
E=
1 2 Z {ψ }T [ K ]{ψ }sin 2 ωt 2
Tmax
：
Emax
1 2 = Z {ψ }T [ K ]{ψ } 2
结构动力学
清华大学土木工程系 刘晶波 2005年秋
结构动力学
刚度矩阵
(2)
2-2）自相关函数 当X(t2)=Y(t2)时,上式结果称相关函数,记作RX(t1,t2).
6.1.3 随机过程的数字特征
3）方差函数
E{[ X (t1 ) m X (t1 )][Y (t 2 ) mY (t 2 )]} C XY (t1 , t 2 ) (3)
5.2 迭代法求频率、振型
5.2.4 求高阶振型和频率 从 {A}0=aiin{X}i (b) 可以看到，如果a1=0（初始向量不包含第一振型时） 像上小节证明，迭代将收敛于第二振型。如果a1=0、 a2=0，则将收敛于第三振型，其余可类推。 教材上利用振型正交性具体推导了滤去低阶振型的 过滤矩阵求法，限于学时，这里不再细说，仅给出过 滤矩阵和动力矩阵的一般公式 过滤矩阵 [Sj]=[Sj-1]-{Aj}{Aj}T[M] (14) 动力矩阵 [Dj+1]=[D][Sj] (15) 利用改造后的动力矩阵，进行迭代即可得到高阶频 率和振型。Vibra程序中包含这种解法。
6.1.4 随机过程按统计特征分类
如果随机过程X(t)的均值为常数，且相关函数只与 时间差=t2-t1有关，则称此过程为弱平稳随机过程。 若随机过程X(t)的任意阶分布函数当时间参数平移 时都保持不变，则称此过程为强平稳随机过程。 一般只讨论弱平稳随机过程。
6.1.5 随机过程的时域特征
补充讲义里还介绍了五点特征，其主要的特征可用 下图示意。(P.3~4)
5.1 能量法求基本频率
5.1.1 单自由度求频率 单自由度无阻尼自由振动解答为Asin(t+)， 当t+=n时位移等于零，因此势能为零，速度为 A，动能为m(A)2/2。 当t+=(n+1)/2时，速度为零、动能等于零，位 移为A，势能为kA2/2。 由无阻尼、能量守恒可得 Tmax=m(A)2/2=kA2/2=EP,max 设=1时最大动能为Tmax，由此即可得 =(EP,max/Tmax)1/2 (1)
基于MATLAB的振动模态分析
摘要振动系统是研究机械振动的运动学和动力学，研究单自由系统的振动有着实际意义，因为工程上有许多问题通过简化，用单自由度系统的振动理论就能得到满意的结果。
模态是振动系统的一种固有振动特性，模态一般包含频率、振型、阻尼。
振动系统问题是个比较虚拟的问题，比较抽象的理论分析，对于问题的分析可以实体化建立数学模型，通过MATLAB可以转化成为图像。
单自由度频率、阻尼、振型的分析，我们可以建立数学模型，最后通过利用MATLAB编程实现数据图形；多自由度主要研究矩阵的迭代求解，我们在分析抽象的理论的同时根据MATLAB编程实现数据的迭代最后可以得到所要的数据，使我们的计算更加简便。
利用MATLAB编程并验证程序的正确性。
通过程序的运行，能快速获得多自由度振动系统的固有频率以及主振型，为设计人员提供了防止系统共振的理论依据，也为初步分析各构件的振动情况以及解耦分析系统响应奠定了基础。
关键词：振动系统；单自由度；MATLAB；多自由度AbstractVibration system is to study the kinematics and dynamics of mechanical vibration, the vibration of a single free system has practical significance, because there are many engineering problems by simplifying, using the vibration theory of a single degree of freedom system can be satisfied with the results.Vibration system problems is a relatively virtual problems, more abstract and theoretical analysis, problem analysis for a mathematical model can be materialized by MATLAB can be converted into images. Single degree of freedom frequency, damping, mode shape analysis, we can create mathematical models, the final program data through the use of MATLAB graphics; many degrees of freedom main matrix iterative solution, our analysis based on abstract theory, while MATLAB programming The last iteration of data can be the desired data, so our calculations easierUsing MATLAB programming and verify the correctness of the program.Through the process of operation, can quickly obtain multiple degrees of freedom vibration system and the main vibration mode natural frequency for the design to prevent resonance provide the theoretical basis for the preliminary analysis of the vibration of each component, and laid the decoupling of system response basis.Key words:vibrating system; Single Degree of Freedom ;MATLAB; multiple degree offreedom辽宁工程技术大学毕业设计（论文）1 绪论1.1问题的提出机械振动是一门既古老又年轻的科学，随着人类科学技术的不断进步振动理论得到不断的发展和完善。
环状分块覆盖ACLD圆柱壳振动分析的新矩阵方法
V ( 3) x
M ( 3) x
0
g g ( 3)
( 3)
12
13
0
g ( 3) 15
上某一参考点( x 0 , 0) 的径向位移与速度[ 5]
U ( ) = - ( k p + kd i ) w ( x 0, 0 )
( 3)
式中 kp 和 kd 分别为位移和速度增益系数。
极化方向为径向的横观各向同性压电层的本构 关系为[ 12]
( 3)
1 ∀3
0
x
( 3)
!( 3) x
=
E3 1 - ∀23
如何选择控制电压的分布形覆盖层对振动响应的抑制效果文中将用具体算例来研究电压分布函数采用pd控制方案施加控制电压受控于基壳上某一参考点irchhoff等效剪力为13环状分块覆盖acl圆柱壳振动分析的新矩阵方法277结合式11沿周向展开成傅立叶级数并将各物理量进行无量钢处理后参照文献10思路和方法经推导和整理可得压电层的一阶常微分矩阵方程5u5v5u5v125w6862105q5u5v叶级数的展开系数其具体表达式为5u5v12中的系数11意义和表达式参见文献105x5u112acld圆柱壳段的整合一阶常微分矩阵方程15可退化为光壳的控制方程参见文献它也是一个一阶常微分矩阵方程10中关于约束层一般弹性材料的一阶常微分矩阵方程用式12替代通过类似推导可得acl向量y中元素称为acl211未知边界状态变量的确定合式17的解由两部分组成
振动力学考题集
1、四个振动系统中，自由度为无限大的是（）。
A. 单摆；B. 质量-弹簧；C. 匀质弹性杆；D. 无质量弹性梁；2、两个分别为c1、c2的阻尼原件,并连后其等效阻尼是（）。
A. c1+c2；B. c1c2/(c1+c2)；C. c1-c2；D. c2-c1；3、（）的振动系统存在为0的固有频率。
A. 有未约束自由度；B. 自由度大于0；C. 自由度大于1；D. 自由度无限多；4、多自由度振动系统中，质量矩阵元素的量纲应该是（）。
A. 相同的，且都是质量；B. 相同的，且都是转动惯量；C. 相同的，且都是密度；D. 可以是不同的；5、等幅简谐激励的单自由度弹簧-小阻尼-质量振动系统，激励频率（）固有频率时，稳态位移响应幅值最大。
A. 等于；B. 稍大于；C. 稍小于；D. 为0；6、自由度为n的振动系统，且没有重合的固有频率，其固有频率的数目（A ）。
A. 为n；B. 为1；C. 大于n；D. 小于n；7、无阻尼振动系统两个不同的振型u(r)和u(s)，u(r)T Mu(s)的值一定（）。
A. 大于0；B. 等于0；C. 小于0；D. 不能确定；8、无阻尼振动系统的某振型u(r)，u(r)T Ku(r)的值一定（）。
A. 大于0；B. 等于0；C. 小于0；D. 不能确定；9、如果简谐激励力作用在无约束振动系统的某集中质量上，当激励频率为无限大时，该集中质量的稳态位移响应一定（）。
A. 大于0；B. 等于0；C. 为无穷大；D. 为一常数值；10、相邻固有频率之间的间隔呈近似无限等差数列的振动系统是（）。
A. 杆的纵向振动；B. 弦的横向振动；C. 一般无限多自由度系统；D. 梁的横向振动；11、两个刚度分别为k1、k2串连的弹簧,其等效刚度是（）。
A. k1+k2；B. k1k2/(k1+k2)；C. k 1-k 2；D. k 2-k 1；12、 无阻尼振动系统两个不同的振型u (r )和u (s )，u (r )T Ku (s )的值一定（ ）。
振动力学课程设计题目
振动力学课程设计题目采用MATLAB 对所选的问题进行数值计算和作图，采用高于MATLAB7.4（2007）版本所编写的程序需转换为文本（.txt ）文件, 早于MATLAB7.4（2007）版本所编写的程序可直接采用M 文件传送至QQ ：296637844。
题目如下，其中1，2，3题为必做题，4-38选二题(第一轮:一班01号为第4题, 一班02号为第5题…一班28号为第25题, 二班01号为第26题,…二班17号为第38题, 二班18号为第4题,…二班27号为第13题;第二轮:一班01号为第14题…)。
文件名采用自己的姓名。
考核时间暂定于12月30日。
题目： 1. 编写MA TLAB 程序，根据书本公式（3.1-10）、（3.1-10）作出单自由度系统强迫振动的幅频特性曲线、相频特性曲线。
0.1,0.2,0.3,0.5,0.7,1.0,1.2ς=。
2. 根据书本图4.5-3，分析有阻尼动力减振器的特性。
包括在不同的质量比，频率比，阻尼比条件下结构的响应。
3. 对于图2所示体系，用矩阵迭代法计算其固有频率及振型。
1231,2m m m ===，1230c c c ===，1231,5,8k k k ===，1230,0,0F F F ===，1231,1,1ωωω===。
4. 采用中心差分法计算单自由度体系10105sin(/2)x cx x t ++= ，当c=3和c=20，000,0x x== 前10s 内的位移，作出其时间位移曲线图。
5. 采用Houbolt 法计算单自由度体系10105sin(/2)x cx x t ++= ，当c=3和c=20，000,0x x== 前10s 内的位移，作出其时间位移曲线图。
6. 采用Wilson-θ法计算单自由度体系10105sin(/2)x cxx t ++= ，当c=3和c=20，000,0x x== 前10s 内的位移，作出其时间位移曲线图。
《振动力学》8线性振动的近似计算方法
ω2 = 1.3213 k / m ω3 = 2.0286 k / m
代入邓克利法公式：
1
2 1
ω
≈
1
ω
2 1
+
1
ω
2 2
+L+
1
ω
2 n
ω1 = 0.3535 k / m
9
线性振动的近似计算方法 / 瑞利法
• 瑞利法
基于能量原理的一种近似方法 可用于计算系统的基频 算出的近似值为实际基频的上限
配合邓克利法算出的基频下限，可以估计实际基频的大致范围。
代入，得：
2 2 R (ϕ ) ≈ ω + ∑ (ω 2 − ω ) ε j k j 2 k j =1 n
ε j << 1
因此，若 ϕ 与 φ( k ) 的差异为一阶小量，则瑞利商与 ωk2的差别 为二阶小量。 2 2 对于基频的特殊情况，令k＝1，则由于 ω j − ω1 > 0( j = 2 ~ n) 瑞利商在基频处取极小值， 利用瑞利商估计系统的基频所得的结果必为实际基频的上限 ϕ 愈接近系统的真实模态，算出的固有频率愈准确。 15
因此，瑞利商的极小值为 ω12
2 ω 同理可证明，瑞利商的极大值为 n
12
线性振动的近似计算方法 / 瑞利法
T n a T ΦN K Φ N a a T Λa 2 R (ϕ ) = T T = T = ∑ a2 ω j j a Φ N MΦ N a a Ia j =1
2 a ∑ j
j =1
n
2 ω12 ≤ R (ϕ ) ≤ ω n
Mφ(i )
= ωi2
(1) N
( 2) N
(n) N
n
《结构动力学》教学日志知识资料b
第
24
次
总结复习
知识点串讲
年月日
学生考核成绩记录
序号
项目
出勤
作业
学号
姓名
/
/
/
/
/
/
/
/
/
/
成
绩
/
/
/
/
/
/
/
/
成
绩
1
5
杨金银
2
2
甄一帆
3
4
周叙霖
4
1
史宝红
5
2
李明聪
6
3
桑胜涛
7
4
崔亚歌
8
5
贾世宁
9
6
连娜
10
7
周文丽
11
8
熊治凯
12
9
薛涛
13
0
周翱翔
14
1
赵锦涛
15
2
田里
16
3
孙可锋
17
4
王浩
教研室主任主管教学院（部）长
年月日年月日
教学计划内容
授课实施记录
课内
课外作业、实验
第
1
次
第1章绪论和概述
1.1结构动力分析主要目的
1.2荷载的分类(持时和来源)
1.3动力问题的基本特性
重点：结构动力分析意义及基本概念。
难点：动力问题与静力问题区别与联系。
寻找1-2本国外结构动力学相关的教材，供学习参考。
（自愿上交）
年月日
第
2
次
第1章绪论和概述
1.4离散化主意
1.5运动方程的建立
振动力学考题集[]
![振动力学考题集[]](https://img.taocdn.com/s3/m/2c2ccd19680203d8cf2f245e.png)
1、四个振动系统中,自由度为无限大的就是( )。
A、单摆;B、质量-弹簧;C、匀质弹性杆;D、无质量弹性梁;2、两个分别为c1、c2的阻尼原件,并连后其等效阻尼就是( )。
A、c1+c2;B、c1c2/(c1+c2);C、c1-c2;D、c2-c1;3、( )的振动系统存在为0的固有频率。
A、有未约束自由度;B、自由度大于0;C、自由度大于1;D、自由度无限多;4、多自由度振动系统中,质量矩阵元素的量纲应该就是( )。
A、相同的,且都就是质量;B、相同的,且都就是转动惯量;C、相同的,且都就是密度;D、可以就是不同的;5、等幅简谐激励的单自由度弹簧-小阻尼-质量振动系统,激励频率( )固有频率时,稳态位移响应幅值最大。
A、等于;B、稍大于;C、稍小于 ;D、为0;6、自由度为n的振动系统,且没有重合的固有频率,其固有频率的数目(A )。
A、为n;B、为1;C、大于n;D、小于n;7、无阻尼振动系统两个不同的振型u(r)与u(s),u(r)T Mu(s)的值一定( )。
A、大于0;B、等于0;C、小于0;D、不能确定;8、无阻尼振动系统的某振型u(r),u(r)T Ku(r)的值一定( )。
A、大于0;B、等于0;C、小于0;D、不能确定;9、如果简谐激励力作用在无约束振动系统的某集中质量上,当激励频率为无限大时,该集中质量的稳态位移响应一定( )。
A、大于0;B、等于0;C、为无穷大;D、为一常数值;10、相邻固有频率之间的间隔呈近似无限等差数列的振动系统就是( )。
A、杆的纵向振动;B、弦的横向振动;C、一般无限多自由度系统;D、梁的横向振动;11、两个刚度分别为k1、k2串连的弹簧,其等效刚度就是( )。
A、k1+k2;B、k1k2/(k1+k2);C、k1-k2;D、k2-k1;12、 无阻尼振动系统两个不同的振型u (r )与u (s ),u (r )T Ku (s )的值一定( )。
A 、 大于0;B 、 等于0;C 、 小于0;D 、 不能确定;13、 无阻尼振动系统的某振型u (r ),u (r )T Mu (r )的值一定( )。
机械振动5多自由度振动7矩阵迭代法.
采用矩阵迭代法，基频：
k 1 0.37308 m
第一阶模态：
u(1) w4 ＝ (0.4626 0.8608 1)T
采用常规方法，基频：
第一阶模态：
u (1) 0.4626 0 . 8608 1
13
1 0.3730 k / m
基频误差为0.2%。模态在精确到四位小数时，完全一致。
2019年2月21日 《振动力学》
从迭代过程看出，获得模态（收敛）的速度取决于 (r / 1 ) k
(r 2,, n) 趋于零的速度。 主要体现在两个方面：
一是，λ1比λ2大多少，相差越大，收敛越快，迭代次数越少；
二是，假设模态选取的准确性，w越接近于第一阶模态u1 ，
收敛速度越快，迭代次数越少。
归一化：
w4 ＝ (0.4626 0.8608 1)T
1 1 2 0.4626 3.3234 第五次迭代： m m Dw 4 ＝ 1 2 4 0.8608 6.1842 k k 1 2 5 1 7.1842
归一化： w5 ＝ (0.4626 0.8608 1)T w4
n n
上式再左乘一次D矩阵： n j ( 1 ) ( j) 2 Du D( Dw) D w 1 C1 Du C j 1 j 2
2019年2月21日 《振动力学》 3
w C1u(1) C2u(2) ＋Cn u(n)
n j ( 1 ) ( j) 2 Du D( Dw) D w 1 C1 Du C j 1 j 2 2 2 n n j ( j) j ( j) 2 (1) (1) u 1 C11u C j u 1 C1u C j 1 1 j 2 j 2 k次左乘D矩阵（相当于k次迭代）后：
迭代法的收敛条件
的情形.
收敛的必要条件，而非充要条件， 因为Gauss-Seidel 迭代即为
1
定理3.6虽然给出了判别迭代法收敛的充要条件， 但实际使用是很不方便 。因为求逆矩阵和特征值的 难度并不亚于用直接方法求解线性方程组。推论１与 推论２使用起来方便得多， 但它们分别给出收敛的 充分条件与必要条件，许多情形下，不能起作用.
( A) A
3
解线性方程组的迭代法
定理3.4 设A为n阶方阵， 则对任意正数 种矩阵范数
使得
,
存在一
A ( A)
证明参看［１］ . 对任意n 阶方阵 A, 一般不存在矩阵范数 使得
,
( A) A . 但若
A为对称矩阵，则
2
( A) A
下面的结论对建立迭代法的收敛性条件非常重要。
k
A lim k
0
所以
k
lim Ak 0
6
解线性方程组的迭代法
3.5.2
迭代法的收敛条件
(0) 定理３.６ 对任意初始向量 x 和右端项 g ,
由迭代格式
x( k 1) Mx( k ) g
产生的向量序列 x ( k ) 收敛的充要条件是
(k 0,1,2,)
19
解线性方程组的迭代法
(2) 若
n阶方阵 A (aij ) 满足
n
aii aij
j 1 j i
(i 1,2,, n)
(3 25)
且至少有一个i 值，使上式中不等号严格成立，则 称为A弱对角占优阵。 定义3.5 如果矩阵Ａ不能通过行的互换和相应列的互 换成为形式
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
k
i
k , xk 1
1 0
§13.3 收敛性的证明 由(11-39)得
高等结构动力学
2 n
n fmn
(13-25*)
则由(13-17)可写为
v1 n Yn Y
1 1 1 n=1 N
（13-26*）
作下一次迭代循环得到第二次循环产生的挠度
2 1
v v
1
1 1 1
T
mv1
0
T
mv11
(13-11)
§13.2 基本振型分析—Stodola法
高等结构动力学
重复上述过程s次，能求得较近似的解,即s次循环之后
v1
s
1
2 1
v1
s 1
1
2 1
1
(13-12)
2 1
max(v1 ) 1 (s) s max(v1 ) max(v1 )
收敛性的证明（第2种方法）
i
1
i K 1Mi
即:
K M
迭代格式
Kxk 1 Mxk
x1 i
i i
x 0 1
Kx2 Mx1
x2 K Mx1 K Mii i K Mi i
1 1 1 i 1 i 1 i 1 N N N
（13-9）
§13.2 基本振型分析—Stodola法
高等结构动力学
真正的第一振型频率介于上式求得的最大值和最小值 之间:
vk0 1 1 v k1 vk0 2 1 1 v 1 min k1
max
(13-10)
取平均值求频率的近似值
(13-43)
§13.4 高阶振型分析
高等结构动力学
最高振型的分析
由(13-1)式得
ˆ ˆ v=Ev
2
(13-44) (13-45)
E m-1k D1
若代入最高振型的试探形状
1 2 N vN
0 Ev N
（13-46）
§13.4 高阶振型分析 第N振型频率的近似值为
2 N
1 T (0) 2 mv3 M2
(13-36a) (13-36b)
代入(13-35)得
(0) 3 v(0) 1 v 3
或
M1
(0) T mv3 1 1
1 (0) 22T mv3 M2
(0) 3 I 1 T mv(0) 1 T mv(0) v(0) v 1 1 3 M2 2 2 3 3 M1
2s 2s
(13-23)
§13.4 高阶振型分析
高等结构动力学
§13.4 高阶振型分析
第二振型分析 展开(13-21)得
0 0 0 v 1 Y1 2 Y2 3 Y3 1 2 3
D
i 1 n
n
2 n
Yn
0
1 2
2
(13-17)
n Dn
2 n
(13-18)
§13.3 收敛性的证明 则由(13-17)可写为
高等结构动力学
1 v1 n Yn （13-19） n=1 n 1 然后，用最大的基准元素 max(v1 )去除 v11 ，使之规格化， 从而得到改进的第一次迭代循环的形状 ，因此
§13.1 引言
高等结构动力学
§13.1 引言
振型位移叠加法提供了一种计算结构动力反应的有效方 法，即无阻尼振型用于对结构运动方程组解耦。 问题——如何获得结构的无阻尼振型？ Stodola法以迭代为基础，先假定初始振型并迭代调整至 实际振型的适当近似,再由运动方程确定振动频率.
§13.2 基本振型分析—Stodola法
高等结构动力学
vkN
1 kN
0 v
T 1 mvN 0 mvN
（13-47a）
2 N
v
vN
1
1 N
T
（13-47b）
其中
vN EvN
1 0
证明最高振型收敛与最低振型收敛的区别是
N 1 N 2 N 3 1 N N N
1
1
i
x3 K 1Mx2
1 i i i 1
N
2
i
§13.3 收敛性的证明
高等结构动力学
1 xk 1 K Mxk i i 1 i
1
N
k
i
1 1
k
1 i i 1 i
(13-16)
§13.3 收敛性的证明 由这些惯性力产生的挠度是
v1
1
高等结构动力学
0 0 2 k -1f I k -1m 11 Y1
2 0 1 22 Y2 2
2
或
v1
1
(13-13)
s-1
§13.3 收敛性的证明
高等结构动力学
§13.3 收敛性的证明
最初假定的形状用正规坐标表示为
v1 Y 1Y 2Y2 3Y3 1
0 0 0 0 0
(13-14)
第一振型频率的振动形状所对应的惯性力为
f I
0
该方程最后一个等号成立,第一振型的系数远大于其它的系数.
1 1 2
2s
1 3
2s
(13-23)
最终结果可视为
v1
s
1Y1(0) 1 (0) max(1Y1 )
(13-24)
§13.3 收敛性的证明
高等结构动力学
(13-40)
§13.4 高阶振型分析 求第四振型计算第三淘汰矩阵
S3 S2
(0) v
4
高等结构动力学
1 T 33 m M3
(13-41) (13-42)
S3 v(0) 4
D4 DS3
相应的动力矩阵 依次类推
Sn Sn-1
1 T nn m Mn
Dn+1 DSn
s s 1 v1 Y n
2s
0 Yn
0 1Y1
(13-28*)
该方程最后一个等号成立,第一振型的系数远大于其它的系数.
1 1 1 2 3
§13.4 高阶振型分析
高等结构动力学
由第一淘汰矩阵减去第二振型的项可得第二淘汰矩阵
1 T S2 S1 22 m M2
(13-38)
(0) (0) v3 S2 v3
其中淘汰矩阵的运算表示为 第三振型的Stodola关系式为
1
(13-39)
2 3
v
(1) 3
(0) 3 =DS v(0) D v(0) Dv 2 3 3 3
高等结构动力学
(13-29)
其中
S1 I 1 11T m M1
(13-30)
在这种情形中,(13-5)能写成
1 v (1) 2 2 2
(0) Dv 2
(13-31)
§13.4 高阶振型分析 代入(13-31)得
高等结构动力学
2
其中
1
v(1) DS1v(0) D2v(0) 2 2 2 2
(13-35)
利用正交特性
T (0) 1 mv3 T (0) 2 mv3
(0) 0 1T mv3 M1Y1(0) (0) 0 2T mv3 M 2Y2(0)
§13.4 高阶振型分析 则
Y1(0)
(0) Y2
高等结构动力学
1 T (0) 1 mv3 M1
ˆ ˆ v=n 2Dvn
v1 =Dv1
1 0
(13-5)
(13-5a)
§13.2 基本振型分析—Stodola法
高等结构动力学
1
2 1
1 =v1 v
1 1
(13-6*)
考虑任一点k的位移
1
12
1 vk0 =vk1 1
(13-7*)
2 1
0 v
k1 1 vk1
(13-27)
M1
不包含第一振型的试探形状为
0 0 20 v v 1Y 2 1
(13-28)
§13.4 高阶振型分析 在试探向量中消除第一振型分量的方便方法是 应用淘汰矩阵S1
2 v 0 v 0 2 1 11T mv 0 S1v 0 2 2 M1
高等结构动力学
பைடு நூலகம்
§13.2 基本振型分析—Stodola法
fI n
ˆ =n 2 mvn
n
（13-1） （13-2）
ˆ vn=k -1f I
或者用式（13-1）则为
ˆ ˆ v=n 2 k-1mvn
可记作
(13-3)
D=k -1m
(13-4)
§13.2 基本振型分析—Stodola法
高等结构动力学
先假定试探形状,它尽可能接近第一振型的形状,而振 幅是任意的,即：
N 4
0
v1
2
（13-21）
s次循环后
v1
s s v1 1 (0) (0) 1 2 S Y 2Y2 ( ) ... (s) 1 1 s 2 max(v1 ) max(v1 )