第12章结构动力学课件.ppt

《第12章结构动力学课件.ppt》由会员分享，可在线阅读，更多相关《第12章结构动力学课件.ppt（155页珍藏版）》请在金锄头文库上搜索。

1、第十四章 结构动力学14-1 概 述14-2 结构振动的自由度14-3 单自由度结构的自由振动14-4 单自由度结构在简谐荷载作用下的强迫振动14-6 多自由度结构的自由振动14-8 振型分解法14-9 无限自由度结构的振动14-5 单自由度结构在任意荷载作用下的强迫振动14-7 多自由度结构在简谐荷载作用下的强迫振动14-10 计算频率的近似法14-1 概 述动力荷载作用下，结构将发生振动，各种量值均随时间而变化。动力荷载作用下，结构将发生振动，各种量值均随时间而变化。一、结构动力计算的特点一、结构动力计算的特点（2 2）研究单自由度及多自由度的自由振动、强迫振动。）研究单自由度及多自由度的

2、自由振动、强迫振动。1 1、内容：、内容：（1 1）研究动力荷载作用下，结构的内力、位移等计算原理和计算方法。）研究动力荷载作用下，结构的内力、位移等计算原理和计算方法。求出它们的最大值并作为结构设计的依据。求出它们的最大值并作为结构设计的依据。2 2、静荷载和动荷载、静荷载和动荷载（1 1）静荷载：荷载的大小和方向不随时间变化（如梁板自重）。）静荷载：荷载的大小和方向不随时间变化（如梁板自重）。（2 2）动荷载：荷载的大小和方向随时间变化，需要考虑惯性力。）动荷载：荷载的大小和方向随时间变化，需要考虑惯性力。 3 3、特点、特点 （2 2）内力与荷载不能构成静平衡。必须考据惯性力。依达朗伯原

3、理，内力与荷载不能构成静平衡。必须考据惯性力。依达朗伯原理，加惯性力后，将动力问题转化为静力问题。加惯性力后，将动力问题转化为静力问题。 (1)(1)必须考虑惯性力。必须考虑惯性力。 (3)(3)分析自由振动即求自振频率、振型、阻尼参数等是求强迫振动分析自由振动即求自振频率、振型、阻尼参数等是求强迫振动动力反应的前提和准备。动力反应的前提和准备。动力荷载的种类动力荷载的种类(1) 周期荷载：随时间按一定规律变化的周期性荷载，如按正弦周期荷载：随时间按一定规律变化的周期性荷载，如按正弦 (或余弦或余弦)规律变化的称为简谐周期荷载，也称为规律变化的称为简谐周期荷载，也称为 振动荷载振动荷载。14-

4、1 概 述F(t)toF(t)=F0sin t(2) 冲击荷载：很快地把全部量值加于结构而作用时间很短即行消失的荷。冲击荷载：很快地把全部量值加于结构而作用时间很短即行消失的荷。F(t)totd(3) 突加荷载：在一瞬间施加于结构上并继续留在结构上的荷载。突加荷载：在一瞬间施加于结构上并继续留在结构上的荷载。(4) 快速移动的荷载。高速移动的列车、汽车等。快速移动的荷载。高速移动的列车、汽车等。F(t)totd14-1 概 述(5) 随机荷载：变化规律不能用确定的函数关系表示的荷载。随机荷载：变化规律不能用确定的函数关系表示的荷载。如风的脉动作用、地震等。如风的脉动作用、地震等。14-1 概

5、述14-1 概 述结构振动的形式结构振动的形式(1) 自由振动：结构受到外部因素干扰发生振动，而在振动自由振动：结构受到外部因素干扰发生振动，而在振动 过程中不再受外部干扰力作用。过程中不再受外部干扰力作用。(2) 强迫振动：在振动过程中不断受外部干扰力作用。强迫振动：在振动过程中不断受外部干扰力作用。 如图所示在跨中支承集中质量的如图所示在跨中支承集中质量的简支梁，把质点简支梁，把质点m拉离原有的弹性平拉离原有的弹性平衡位置，然后突然放松，则质点将在衡位置，然后突然放松，则质点将在原有平衡位置附近往复振动。在振动原有平衡位置附近往复振动。在振动过程中不受外来干扰，这时的振动即过程中不受外来干

6、扰，这时的振动即是是自由振动自由振动。14-2 结构振动的自由度结构振动的结构振动的自由度自由度：结构在弹性变形过程中确定全部质点位置所需的独：结构在弹性变形过程中确定全部质点位置所需的独 立参数的数目。立参数的数目。 图图a所示简支梁跨中固定一个重量较所示简支梁跨中固定一个重量较大的物体，如果梁本身的自重较小可大的物体，如果梁本身的自重较小可略去，把重物简化为一个集中质点，略去，把重物简化为一个集中质点，得到图得到图b所示的计算简图。所示的计算简图。梁在振动中的自由度梁在振动中的自由度=1单自由度结构单自由度结构具有一个自由度的结构。具有一个自由度的结构。多自由度结构多自由度结构自由度大于自

7、由度大于1的结构。的结构。14-2 结构振动的自由度图图a所示结构有三个集中质点。所示结构有三个集中质点。自由度自由度=1图图b所示简支梁上有三个集中质量。所示简支梁上有三个集中质量。自由度自由度=3图图c所示刚架有一个集中质点。所示刚架有一个集中质点。自由度自由度=2自由度的数目不完全取决于质点的数目自由度的数目不完全取决于质点的数目14-2 结构振动的自由度 图刚架上有四个集中质点，但只需要加三根链杆图刚架上有四个集中质点，但只需要加三根链杆便可限制全部质点的位置。如图便可限制全部质点的位置。如图e。自由度自由度=3图示梁，其分布质量集度为图示梁，其分布质量集度为m，可看作有无穷多，可看作

8、有无穷多个个mdx的集中质量，是无限自由度结构。的集中质量，是无限自由度结构。自由度的数目与结构是否静定或超静定无关自由度的数目与结构是否静定或超静定无关动力自由度的确定方法：加附加链杆约束质点位移，最少链杆数即为自动力自由度的确定方法：加附加链杆约束质点位移，最少链杆数即为自由度由度或或 图图a所示机器的块式基础，当机所示机器的块式基础，当机器运转时，若只考虑基础的垂直振器运转时，若只考虑基础的垂直振动，可用弹簧表示地基的弹性，用动，可用弹簧表示地基的弹性，用一个集中质量代表基础的质量。使一个集中质量代表基础的质量。使结构转化为图示的单自由度结构。结构转化为图示的单自由度结构。14-2 结构

9、振动的自由度 图图b所示的水塔，顶部水池较重，所示的水塔，顶部水池较重，塔身重量较轻，略去次要因素后，塔身重量较轻，略去次要因素后，可简化为图示的直立悬臂梁在顶端可简化为图示的直立悬臂梁在顶端支承集中质量的单自由度结构。支承集中质量的单自由度结构。实际结构针对具体问题可以进行简化实际结构针对具体问题可以进行简化14-3 单自由度结构的自由振动 图图a所示为一个简单的质点弹簧模型。取重所示为一个简单的质点弹簧模型。取重物的静力平衡位置为计算位移物的静力平衡位置为计算位移y的原点，规的原点，规定位移定位移y和质点所受的力都已和质点所受的力都已向下为正向下为正。(1) 列动力平衡方程（刚度法）列动力

10、平衡方程（刚度法）弹簧拉力弹簧拉力(恢复力恢复力) Fe=k11y惯性力惯性力 质点处于动力平衡状态质点处于动力平衡状态可得可得一、不考虑阻尼时的自由振动一、不考虑阻尼时的自由振动取振动任一时刻的质点为隔离体如图取振动任一时刻的质点为隔离体如图b。1、振动微分方程的建立、振动微分方程的建立14-3 单自由度结构的自由振动(2) 列位移方程（柔度法）如图列位移方程（柔度法）如图c。 质点质点m振动时，把惯性力振动时，把惯性力FI看作是静力荷看作是静力荷载作用在体系上，则质点处的位移为载作用在体系上，则质点处的位移为对单自由度结构有对单自由度结构有可得与可得与(1)相同的结果相同的结果或为或为(a

11、)命命上式即为单自由度结构自由振动微分方程上式即为单自由度结构自由振动微分方程则有则有(a)14-3 单自由度结构的自由振动建立振动微分方程的例：建立振动微分方程的例：14-3 单自由度结构的自由振动建立图示体系的振动微分方程：建立图示体系的振动微分方程：方程方程振动的初始条件为振动的初始条件为则有则有可得可得2、运动方程的解：、运动方程的解：为一常系数线性齐次微分方程，其通解为为一常系数线性齐次微分方程，其通解为A1和和A2为任意常数，可有初始条件来确定。为任意常数，可有初始条件来确定。式中式中y0初位移，初位移， 初速度。初速度。14-3 单自由度结构的自由振动14-3 单自由度结构的自由

12、振动结构的自由振动由两部分组成：结构的自由振动由两部分组成：一部分是初位移一部分是初位移y0引起的，为余弦规律；引起的，为余弦规律；一部分是初速度一部分是初速度 引起的，为正弦规律。如图引起的，为正弦规律。如图a、b。14-3 单自由度结构的自由振动令令则有则有式式(b)可写为可写为(c)简谐振动如图简谐振动如图ca 为振幅，表示质点的最大位移；为振幅，表示质点的最大位移； 为初相角。为初相角。周期周期工程频率工程频率角频率或频率角频率或频率讨论：结构振动主要由三个参数讨论：结构振动主要由三个参数a、和和 有关。有关。a和和与外因（初与外因（初位移、初速度）有关，位移、初速度）有关， 只与结构

13、特性有关，是结构固有特性，只与结构特性有关，是结构固有特性，决定了结构的动力特性，即两个结构只要决定了结构的动力特性，即两个结构只要相同，动力反应相同。相同，动力反应相同。14-3 单自由度结构的自由振动(d)g重力加速度；重力加速度；st重量重量mg所产生静力位移。所产生静力位移。式式(d)表明：表明：随随st的增大而减小，即把质点放在结构最大位的增大而减小，即把质点放在结构最大位 移处，则可得到最低的自振频率和最大的振动周期。移处，则可得到最低的自振频率和最大的振动周期。刚度法刚度法柔度法柔度法重力法重力法讨论：质量自重力对自振频率的影响。讨论：质量自重力对自振频率的影响。yst为重力为重

14、力mg产生的静位移，产生的静位移，y为动位移，总位移为为动位移，总位移为yst+y，动平衡方程：动平衡方程：考虑质量重力，不影响频率和考虑质量重力，不影响频率和动位移动位移。14-3 单自由度结构的自由振动解：可用柔度法计算，即先求单位力产生的位移解：可用柔度法计算，即先求单位力产生的位移11，代入公式计算，代入公式计算例例14-1 当不考虑梁的自重时，比较图中所示三种支承情况的梁当不考虑梁的自重时，比较图中所示三种支承情况的梁 的自振周期。的自振周期。自乘自乘14-3 单自由度结构的自由振动说明：随着结构刚度的增大，其自振频率说明：随着结构刚度的增大，其自振频率也相应地增高。也相应地增高。(

15、a),(b)图互乘图互乘据此有据此有11为为超静定结构位移计算，虚设状态为静定结超静定结构位移计算，虚设状态为静定结构，可取上页的图构，可取上页的图(a) ，(a),(c)图互乘图互乘上面几种情况刚上面几种情况刚度系数度系数k1111=?=?14-3 单自由度结构的自由振动例例 求下面结构的自振频率。求下面结构的自振频率。解：解：水平方向振动时，水平方向振动时，总质量为总质量为2m，故，故14-3 单自由度结构的自由振动例例 求下面结构的自振频率。求下面结构的自振频率。解：解： 此题结构为剪切型刚架，用刚度法计算此题结构为剪切型刚架，用刚度法计算频率较简单。频率较简单。14-3 单自由度结构的

16、自由振动上面杆端剪力称为杆件的侧移刚度，同层各杆侧移刚度之和称为上面杆端剪力称为杆件的侧移刚度，同层各杆侧移刚度之和称为结构的层间刚度。结构的层间刚度。杆件的侧移刚度与杆端约束有关杆件的侧移刚度与杆端约束有关例如例如14-3 单自由度结构的自由振动例例 分别用柔度法和刚度法计算图示结构的自振频率，分别用柔度法和刚度法计算图示结构的自振频率，EI=常数。常数。解：解：1、柔度法、柔度法解：解：2、刚度法、刚度法14-3 单自由度结构的自由振动求刚度系数求刚度系数附加链杆的水平单位位移引起的附加反力附加链杆的水平单位位移引起的附加反力即为刚度系数。即为刚度系数。用力矩分配法计算用力矩分配法计算AB

17、BABCCB0.60.40-12i/l007.2i/l4.8i/l2.4i/l0-4.8i/l4.8i/l2.4i/l14-3 单自由度结构的自由振动例例 图示结构杆件刚度为无穷大，试求其自振频率。图示结构杆件刚度为无穷大，试求其自振频率。(不计杆件质量不计杆件质量)解：本问题可由转动惯量的方法计算。解：本问题可由转动惯量的方法计算。J为转动惯量，为转动惯量，k11为转动刚度。为转动刚度。绕绕A点的转动惯量点的转动惯量J为：为：绕绕A点的转动刚度点的转动刚度k11为：为：14-3 单自由度结构的自由振动推导运动微分方程，将加速度项的系数简化推导运动微分方程，将加速度项的系数简化为为1，则位移项

18、的系数为自振频率的平方。，则位移项的系数为自振频率的平方。14-3 单自由度结构的自由振动例例 试求图示结构的自振频率。试求图示结构的自振频率。(不计杆件质量不计杆件质量)解：解：14-3 单自由度结构的自由振动2、考虑阻尼作用时的自由振动、考虑阻尼作用时的自由振动阻尼力的产生：外部介质的阻力，支承的摩擦等；阻尼力的产生：外部介质的阻力，支承的摩擦等； 物体内部的作用，材料分子之间的摩擦等。物体内部的作用，材料分子之间的摩擦等。 粘滞阻尼力：阻尼力与其振动的速度成正比，与速度的方向相反。粘滞阻尼力：阻尼力与其振动的速度成正比，与速度的方向相反。c 称为阻尼系数称为阻尼系数考虑阻尼力时，质点考虑

19、阻尼力时，质点m的受力图如图所示的受力图如图所示由动力平衡得由动力平衡得即即令令则有则有线性常系数齐次微分方程线性常系数齐次微分方程建立振动微分方程建立振动微分方程14-3 单自由度结构的自由振动设其解为设其解为代入微分方程，得特征方程代入微分方程，得特征方程两个根为两个根为讨论讨论(1) 大阻尼情况：大阻尼情况：r1、r2是两个负实数，方程的通解为是两个负实数，方程的通解为 是非周期函数，不产生振动，结构偏离平衡位置后将缓慢回复到原有位置。是非周期函数，不产生振动，结构偏离平衡位置后将缓慢回复到原有位置。(3) =临界阻尼情况：临界阻尼情况：r1=r2=-，方程的通解为，方程的通解为是非周期

20、函数，不发生振动。是非周期函数，不发生振动。此时阻尼比此时阻尼比=1，=m，可得临界阻尼系数，可得临界阻尼系数故有故有阻尼比为阻尼系数与临界阻尼系数之比。阻尼比为阻尼系数与临界阻尼系数之比。临界阻尼系数临界阻尼系数ccr是结构的固有特性。是结构的固有特性。14-3 单自由度结构的自由振动小小 结结自由振动分为无阻尼自由振动和有阻尼自由振动。自由振动分为无阻尼自由振动和有阻尼自由振动。无阻尼自由振动特点：无阻尼自由振动特点：结构作简谐振动，振动由初始位移和初始速度引起。结构作简谐振动，振动由初始位移和初始速度引起。振动的主要动力特性与自振频率有关，自振频率与结构的刚振动的主要动力特性与自振频率有

21、关，自振频率与结构的刚度成正比，与结构质量成反比。度成正比，与结构质量成反比。有阻尼自由振动特点：有阻尼自由振动特点：当阻尼系数小于临界阻尼系数或阻尼比小于当阻尼系数小于临界阻尼系数或阻尼比小于1时，结构振动，且时，结构振动，且为衰减的简谐振动。为衰减的简谐振动。有阻尼的自振频率小于无阻尼时的自振频率。有阻尼的自振频率小于无阻尼时的自振频率。阻尼和自振频率是反映有阻尼振动的主要特性，阻尼特性可由阻阻尼和自振频率是反映有阻尼振动的主要特性，阻尼特性可由阻尼比表示，阻尼比可由实测结构的相邻周期的振幅比值来计算。尼比表示，阻尼比可由实测结构的相邻周期的振幅比值来计算。有阻尼的振幅是随时间减少的。有阻

22、尼的振幅是随时间减少的。当阻尼系数大于等于临界阻尼系数或阻尼比大于等于当阻尼系数大于等于临界阻尼系数或阻尼比大于等于1时，结构不时，结构不振动，结构在偏离平衡位置后将缓慢回复到原有位置。振动，结构在偏离平衡位置后将缓慢回复到原有位置。自由振动是结构初始状态引起的振动，振动中无外界干扰。自由振动是结构初始状态引起的振动，振动中无外界干扰。14-4 单自由度结构在简谐荷载作用下的强迫振动强迫振动强迫振动结构在外来干扰力结构在外来干扰力F(t)作用下产生的振动。作用下产生的振动。如图所示，干扰力如图所示，干扰力F(t)直接直接作用在质点作用在质点m上上即即或或结构在使用过程中遇到振动问题主要是强迫振

23、动。结构在使用过程中遇到振动问题主要是强迫振动。 前面讨论的结构自由振动主要是讨论结构的自身的动力特性，而前面讨论的结构自由振动主要是讨论结构的自身的动力特性，而结构在强迫振动的时的效果是与结构的动力特性有关的。结构在强迫振动的时的效果是与结构的动力特性有关的。强迫振动运动微分方程的推导。强迫振动运动微分方程的推导。简谐荷载：简谐荷载：F(t)=F0sin t ，F0为荷载幅值（荷载的最为荷载幅值（荷载的最大值），大值），为荷载频率（圆频率）。为荷载频率（圆频率）。取质点为隔离体，可得取质点为隔离体，可得14-4 单自由度结构在简谐荷载作用下的强迫振动微分方程的解为齐次解与特解之和。微分方程的

24、解为齐次解与特解之和。特解特解 为满足方程的任意解。为满足方程的任意解。右端项为零的齐次方程的右端项为零的齐次方程的 y0，即为自由振动解，即为自由振动解1、无阻尼强迫振动、无阻尼强迫振动或或设设代入微分方程，得代入微分方程，得解得解得即即14-4 单自由度结构在简谐荷载作用下的强迫振动方程通解为：方程通解为：任意常数任意常数B1和和B2可由初始条件可由初始条件t=0时，时，来确定。来确定。结构振动由两部分叠加：初始条件引起的自由振动和动荷载引起的结构振动由两部分叠加：初始条件引起的自由振动和动荷载引起的纯强迫振动。前者按结构自振频率振动，后者按荷载频率振动。纯强迫振动。前者按结构自振频率振动

25、，后者按荷载频率振动。短时间内，自由振动部分会衰减而被忽略，只剩纯强迫振动部分。短时间内，自由振动部分会衰减而被忽略，只剩纯强迫振动部分。振动可分为两个阶段振动可分为两个阶段开始时各部分振动同时存在，此阶段称为过渡阶段。开始时各部分振动同时存在，此阶段称为过渡阶段。短时间后，自由振动影响可忽略，仅有强迫振动，此阶段短时间后，自由振动影响可忽略，仅有强迫振动，此阶段称为平稳阶段。称为平稳阶段。过渡阶段比较短，实际问题中平稳阶段比较重要，因此这里着过渡阶段比较短，实际问题中平稳阶段比较重要，因此这里着重讨论平稳阶段的纯强迫振动，也称为稳态强迫振动。重讨论平稳阶段的纯强迫振动，也称为稳态强迫振动。1

26、4-4 单自由度结构在简谐荷载作用下的强迫振动平稳阶段的纯强迫振动或稳态强迫振动的表达式为平稳阶段的纯强迫振动或稳态强迫振动的表达式为式中式中为已知，主要讨论振幅为已知，主要讨论振幅A。式中式中得得式中式中yst为荷载最大值产生的静位移，为荷载最大值产生的静位移， 为位移动力系数，为位移动力系数，A为最大动位移。为最大动位移。表明考虑动荷载时结构的最大动位移为荷载最大值产生的静位移的表明考虑动荷载时结构的最大动位移为荷载最大值产生的静位移的 倍。倍。位移动力系数位移动力系数相当于放大系数。相当于放大系数。14-4 单自由度结构在简谐荷载作用下的强迫振动在结构设计中，以结构反应的最大值作为设计依

27、据，结构动力反应的最在结构设计中，以结构反应的最大值作为设计依据，结构动力反应的最大值可以通过计算荷载最大值产生的反应乘以动力系数得到，这样将动大值可以通过计算荷载最大值产生的反应乘以动力系数得到，这样将动力问题转化为静力问题来解答，简化了计算。因此，动力系数的计算是力问题转化为静力问题来解答，简化了计算。因此，动力系数的计算是强迫振动计算中的一项重要内容。强迫振动计算中的一项重要内容。动力系数仅与结构的自振频率和荷载频率有关，相同荷载作用下，结构动力系数仅与结构的自振频率和荷载频率有关，相同荷载作用下，结构自振频率不同，结构的最大反应不同。自振频率不同，结构的最大反应不同。结构反应分为内力反

28、应和位移（变形）反应，动力系数也分为内力动力结构反应分为内力反应和位移（变形）反应，动力系数也分为内力动力系数和位移动力系数。系数和位移动力系数。对单自由度体系，当荷载作用于质点上时（即惯性力与荷载作用点和方对单自由度体系，当荷载作用于质点上时（即惯性力与荷载作用点和方向相同时），内力动力系数与位移动力系数相同，统称为动力系数。向相同时），内力动力系数与位移动力系数相同，统称为动力系数。对于多自由度体系和当荷载不作用于质点上时，没有统一的位移动力系对于多自由度体系和当荷载不作用于质点上时，没有统一的位移动力系数和内力动力系数。数和内力动力系数。14-4 单自由度结构在简谐荷载作用下的强迫振动当

29、当时：时：为负，动力位移与动力荷载反向。为负，动力位移与动力荷载反向。 随随/ 而变化，当干扰力频率而变化，当干扰力频率接近于结构的自振频率接近于结构的自振频率时，动力时，动力系数迅速增大；系数迅速增大； =时，理论上时，理论上无穷大，此时内力和位移都将无限大无穷大，此时内力和位移都将无限大共振共振。工程设计中应尽量避免发生共振工程设计中应尽量避免发生共振 通常取通常取/ 1.25。14-4 单自由度结构在简谐荷载作用下的强迫振动例例 计算下面结构在干扰力作用下质量的最大位移和结构的最大弯矩。计算下面结构在干扰力作用下质量的最大位移和结构的最大弯矩。已知干扰力频率为已知干扰力频率为 ，F0=1

30、0kN，l =4m。 解：解：14-4 单自由度结构在简谐荷载作用下的强迫振动例例14-2 如图发电机的重量如图发电机的重量G=35kN，梁的，梁的I=8.810-5m4，E=210GPa，发，发电机转动时离心力的垂直分力幅值电机转动时离心力的垂直分力幅值F=10kN。不考虑阻尼，试求当发。不考虑阻尼，试求当发电机转数为电机转数为n=500r/min时，梁的最大弯矩和挠度（不计梁的自重）。时，梁的最大弯矩和挠度（不计梁的自重）。解：在解：在G作用下，梁中点的最大静位移为作用下，梁中点的最大静位移为自振频率为自振频率为干扰力频率为干扰力频率为求得动力系数求得动力系数梁中点的最大弯矩梁中点的最大弯

31、矩梁中点最大挠度梁中点最大挠度14-4 单自由度结构在简谐荷载作用下的强迫振动例例 计算下面结构在干扰力作用下质量的最大位移和结构的最大弯矩。计算下面结构在干扰力作用下质量的最大位移和结构的最大弯矩。已知干扰力频率为已知干扰力频率为 ，F0=10kN，l =4m。14-4 单自由度结构在简谐荷载作用下的强迫振动 图图a所示简支梁，干扰力不作用在质点上。建立质点所示简支梁，干扰力不作用在质点上。建立质点m的振动方程。的振动方程。 F=1作用在点作用在点1时使点时使点1产生的位移为产生的位移为11，如图，如图b。 F=1作用在点作用在点2时使点时使点1产生的位移为产生的位移为12，如图，如图c。作

32、用在质点作用在质点m上的惯性力为上的惯性力为 在惯性力在惯性力FI和干扰力和干扰力F(t)共同作用下，共同作用下，任一时刻质点任一时刻质点m处的位移为处的位移为即即 可以看成作用于质量上的等效可以看成作用于质量上的等效荷载。与荷载。与F(t)产生相同的质点位移，但产生相同的质点位移，但其他位移和内力不同。其他位移和内力不同。14-4 单自由度结构在简谐荷载作用下的强迫振动此时，位移动力系数与内力动力系数不同。此时，位移动力系数与内力动力系数不同。等效质点动荷载等效质点动荷载综上所述，结构在简谐荷载综上所述，结构在简谐荷载F=F0sint作用下，作用下，速度为速度为加速度为加速度为作用于质点上的

33、惯性力为作用于质点上的惯性力为位移为位移为可以看出，干扰力、位移和惯性力是同步的，即同时达到最大值。可以看出，干扰力、位移和惯性力是同步的，即同时达到最大值。14-4 单自由度结构在简谐荷载作用下的强迫振动因此结构的最大反应可以看成是由荷载最大值和惯性力最大值共同因此结构的最大反应可以看成是由荷载最大值和惯性力最大值共同作用产生的。作用产生的。惯性力最大值为：惯性力最大值为：将将F0和和 2Am 同时作用于结构上，按静力学的方法就可求出结构的同时作用于结构上，按静力学的方法就可求出结构的最大位移和最大内力。最大位移和最大内力。14-4 单自由度结构在简谐荷载作用下的强迫振动例例 计算下面结构在

34、干扰力作用下质量的最大位移和结构的最大弯矩。计算下面结构在干扰力作用下质量的最大位移和结构的最大弯矩。已知干扰力频率为已知干扰力频率为 ，F0=10kN，l =4m。解得：解得：14-4 单自由度结构在简谐荷载作用下的强迫振动将荷载最大值和惯性力最大值同时作用到结构上，并画出弯矩图将荷载最大值和惯性力最大值同时作用到结构上，并画出弯矩图由弯矩图可知最大弯矩为由弯矩图可知最大弯矩为14-4 单自由度结构在简谐荷载作用下的强迫振动例例 计算下面结构在干扰力作用下质量的最大位移和结构的最大弯矩。计算下面结构在干扰力作用下质量的最大位移和结构的最大弯矩。已知干扰力频率为已知干扰力频率为 ，F0=10k

35、N，l =4m。14-4 单自由度结构在简谐荷载作用下的强迫振动动弯矩幅值图动弯矩幅值图最大位移发最大位移发生在哪里？生在哪里？值为多少？值为多少？14-4 单自由度结构在简谐荷载作用下的强迫振动2、有阻尼强迫振动、有阻尼强迫振动微分方程的解：齐次方程的解微分方程的解：齐次方程的解 y0，与干扰力与干扰力F(t)相应的特解相应的特解设特解为设特解为代入方程代入方程解出解出14-4 单自由度结构在简谐荷载作用下的强迫振动将将y0与特解合并，由初始条件与特解合并，由初始条件可得可得14-4 单自由度结构在简谐荷载作用下的强迫振动(1) 由初始条件决定的自由振动；由初始条件决定的自由振动；(2) 伴

36、随干扰力的作用发生的振动频率为伴随干扰力的作用发生的振动频率为，称为伴生自由振动；，称为伴生自由振动；(3) 按干扰力频率按干扰力频率振动，称为振动，称为纯强迫振动纯强迫振动或或稳态强迫振动稳态强迫振动如图。如图。前两部分振动很快衰减掉，最后只剩下纯强迫振动。前两部分振动很快衰减掉，最后只剩下纯强迫振动。过渡阶段过渡阶段振动开始的一段时间内几种振动同时存在的阶段；振动开始的一段时间内几种振动同时存在的阶段；平稳阶段平稳阶段纯强迫振动阶段。纯强迫振动阶段。由解答的表达式可知，振动由三部分组成：由解答的表达式可知，振动由三部分组成：振动可分为两个阶段振动可分为两个阶段过渡阶段比较短，实际问题中平稳

37、阶段比较重要，因此这里着过渡阶段比较短，实际问题中平稳阶段比较重要，因此这里着重讨论平稳阶段的纯强迫振动，即解答的第三部分。重讨论平稳阶段的纯强迫振动，即解答的第三部分。14-4 单自由度结构在简谐荷载作用下的强迫振动将第三项写为将第三项写为振幅振幅相位差相位差振幅振幅A可写为可写为动力系数动力系数* * 注：有阻尼的纯强迫振动不是衰减的。注：有阻尼的纯强迫振动不是衰减的。注：有阻尼的纯强迫振动不是衰减的。注：有阻尼的纯强迫振动不是衰减的。14-4 单自由度结构在简谐荷载作用下的强迫振动动力系数动力系数与与/及及的关系如的关系如图所示。图所示。14-4 单自由度结构在简谐荷载作用下的强迫振动相

38、位差相位差与与/及及的关系如的关系如图所示。图所示。14-4 单自由度结构在简谐荷载作用下的强迫振动讨论讨论(1) 时，时，很小，质量近似于不动或作振幅很微小的颤动。结构的很小，质量近似于不动或作振幅很微小的颤动。结构的Fe、FR可以忽略，位移与荷载的相位差为可以忽略，位移与荷载的相位差为180。荷载主要由惯性力平荷载主要由惯性力平衡。衡。(3) 时，时，增加很快，增加很快，受阻尼的影响很大受阻尼的影响很大 。荷载主要由阻尼。荷载主要由阻尼力平衡。当阻尼较小时，力平衡。当阻尼较小时，值很大，共振现象仍很危险。值很大，共振现象仍很危险。工程设计中一般常取工程设计中一般常取14-5 单自由度结构在

39、任意荷载作用下的强迫振动瞬时冲量：荷载瞬时冲量：荷载F(t)在极短的时间在极短的时间t0内给与内给与振动物体的冲量振动物体的冲量瞬时冲量作用下的振动问题瞬时冲量作用下的振动问题 图图a所示荷载大小为所示荷载大小为F，作用时间为，作用时间为t ，其冲量其冲量I=Ft ，即图中阴影部分的面积。，即图中阴影部分的面积。瞬时冲量作用下质点的动量增值为瞬时冲量作用下质点的动量增值为由由可得可得 当质点获得初速度后冲量即时消失，质点在这种冲击下将产生自由振当质点获得初速度后冲量即时消失，质点在这种冲击下将产生自由振动。将初始条件代入式动。将初始条件代入式(g)可得瞬时冲量可得瞬时冲量I作用下质点作用下质点

40、m的位移方程为的位移方程为14-5 单自由度结构在任意荷载作用下的强迫振动 若瞬时冲量不是在若瞬时冲量不是在t=0而是在而是在t=时加于质点上，其位移方程为时加于质点上，其位移方程为 图图b所示一般形式的干扰力所示一般形式的干扰力F(t)可认为可认为是一系列微小冲量是一系列微小冲量F()d连续作用的结果，连续作用的结果，应此有应此有(k)不考虑阻尼不考虑阻尼=0，=则有则有(m)式式(k)及式及式(m)称为称为杜哈梅积分杜哈梅积分14-5 单自由度结构在任意荷载作用下的强迫振动若在若在t=0质点原来还具有初始位移和初始速度，则质点位移为质点原来还具有初始位移和初始速度，则质点位移为若不考虑阻尼

41、则有若不考虑阻尼则有(n)14-5 单自由度结构在任意荷载作用下的强迫振动（1）突加荷载。变化规律如图）突加荷载。变化规律如图a所示。所示。设：加载前结构处于静止状态，将设：加载前结构处于静止状态，将 F()=F代入式代入式(k)求得求得其振动曲线如图其振动曲线如图b。时最大动位移时最大动位移yd为为动力系数为动力系数为不考虑阻尼不考虑阻尼14-5 单自由度结构在任意荷载作用下的强迫振动（2）短期荷载。变化规律如图所示。）短期荷载。变化规律如图所示。当当t=0时，时，有突加荷载加入并一直作用在结构上；有突加荷载加入并一直作用在结构上；当当t=t0时，时，有一个大小相等方向相反的突加荷载加入。有

42、一个大小相等方向相反的突加荷载加入。利用（利用（1）得到的突加荷载作用下的计算公式按叠加法求解：）得到的突加荷载作用下的计算公式按叠加法求解：自由振动自由振动当当t0T/2时，最大位移发生在前一阶段。时，最大位移发生在前一阶段。短期荷载的最大动力效应与突加荷载相同。短期荷载的最大动力效应与突加荷载相同。14-5 单自由度结构在任意荷载作用下的强迫振动14-6 多自由度结构的自由振动 工程实际中有很多结构是不宜简化为单自由度体工程实际中有很多结构是不宜简化为单自由度体系计算的。例如多层房屋、多跨不等高工业厂房以及系计算的。例如多层房屋、多跨不等高工业厂房以及烟囱等烟囱等, , 都必须按多自由度体

43、系来处理。都必须按多自由度体系来处理。 图示等截面烟囱，图示等截面烟囱，将其分为八段，从上到将其分为八段，从上到下将每两段的质量集中下将每两段的质量集中于其中点，将一个无限于其中点，将一个无限自由度的体系简化为四自由度的体系简化为四个自由度体系。个自由度体系。1、振动微分方程的建立、振动微分方程的建立14-6 多自由度结构的自由振动柔度法柔度法柔度矩阵柔度矩阵质量矩阵质量矩阵加速度向量加速度向量位移向量位移向量14-6 多自由度结构的自由振动刚度法刚度法弹性恢复力如何计算？弹性恢复力如何计算？Fe1=- -FR1， Fe2=- -FR2FR1和和FR2如何计算？如何计算？由叠加原理：由叠加原理

44、：由质量的平衡条件，得由质量的平衡条件，得其中其中其中其中( ( ( (质量矩阵）质量矩阵）质量矩阵）质量矩阵）( ( ( (位移向量）位移向量）位移向量）位移向量）( ( ( (加速度向量）加速度向量）加速度向量）加速度向量）或或或或( ( ( (刚度矩阵）刚度矩阵）刚度矩阵）刚度矩阵）14-6 多自由度结构的自由振动14-6 多自由度结构的自由振动由由和和可知上面两个方程是同一方程的两个表达方式可知上面两个方程是同一方程的两个表达方式与单自由度的运动方程比较，将质量、刚度、加速度和位移分别与单自由度的运动方程比较，将质量、刚度、加速度和位移分别换成质量矩阵，刚度矩阵，加速度向量和位移向量，

45、方程的形式换成质量矩阵，刚度矩阵，加速度向量和位移向量，方程的形式相似。（便于记忆）如相似。（便于记忆）如柔度法柔度法将各质点的惯性力看作是静荷载将各质点的惯性力看作是静荷载如图如图a。结构上任一质点结构上任一质点mi处的位移应为处的位移应为ii、ij为柔度系数其物理意义见图为柔度系数其物理意义见图b、c。由此，可以建立由此，可以建立n个位移方程个位移方程多自由度结构无阻尼自由多自由度结构无阻尼自由振动微分方程振动微分方程推广到推广到n个自由度的情况个自由度的情况14-6 多自由度结构的自由振动14-6 多自由度结构的自由振动写成矩阵形式为写成矩阵形式为简写为简写为为结构的柔度矩阵，是对称矩阵

46、。为结构的柔度矩阵，是对称矩阵。M为质量矩阵，在集中质点的结构中是对角矩阵；为质量矩阵，在集中质点的结构中是对角矩阵； 为加速度列向量；为加速度列向量；Y为位移列向量。为位移列向量。14-6 多自由度结构的自由振动刚度法刚度法 图图a所示无重量简支梁，所示无重量简支梁，略去梁的轴向变形和质点的转略去梁的轴向变形和质点的转动，为动，为n个自由度的结构。个自由度的结构。加入附加链杆阻止所有质点的加入附加链杆阻止所有质点的位移，如图位移，如图b。各质点的惯性力为各质点的惯性力为各链杆的反力为各链杆的反力为令各链杆发生与各质点实际位令各链杆发生与各质点实际位置相同的位移，如图置相同的位移，如图c。各链

47、杆上所需施加的力为各链杆上所需施加的力为14-6 多自由度结构的自由振动不计阻尼，各链杆上的总反力应等于零。不计阻尼，各链杆上的总反力应等于零。以质点以质点mi为例有为例有kii、kij为刚度系数其物理意义见图为刚度系数其物理意义见图d、e。可得可得i质点的动力平衡方程为质点的动力平衡方程为14-6 多自由度结构的自由振动对每个质点都列出一个动力平衡方程，于是可得对每个质点都列出一个动力平衡方程，于是可得写成矩阵形式为写成矩阵形式为多自由度结构无阻尼自由振动微分方程多自由度结构无阻尼自由振动微分方程14-6 多自由度结构的自由振动简写为简写为式中：式中：K 为刚度矩阵，是对称矩阵；为刚度矩阵，

48、是对称矩阵；由由柔度矩阵与刚度矩阵是互为逆阵。柔度矩阵与刚度矩阵是互为逆阵。记忆方法：和单自由度比较记忆方法：和单自由度比较刚度法刚度法柔度法柔度法（形式与单自由度相同）（形式与单自由度相同）（形式与单自由度相同）（形式与单自由度相同） 微分方程解答微分方程解答振动特点振动特点14-6 多自由度结构的自由振动设位移方程的特解为设位移方程的特解为代入位移方程可得代入位移方程可得振幅方程振幅方程1、按柔度法求解、按柔度法求解这是一组各质点按同一频率振动，只是振幅不同。这是一组各质点按同一频率振动，只是振幅不同。14-6 多自由度结构的自由振动写成矩阵形式写成矩阵形式式中式中振幅列向量振幅列向量I

49、为单位矩阵。为单位矩阵。要得到振幅不全为零的解答，振幅方程组的系数行列式为零。要得到振幅不全为零的解答，振幅方程组的系数行列式为零。频率方程频率方程14-6 多自由度结构的自由振动 将行列式展开将行列式展开含含 2的的n次代数方程，从而可得到次代数方程，从而可得到n个自个自振频率振频率1，2，n，将频率从小到大排列，分别称为第，将频率从小到大排列，分别称为第一，第二，一，第二， ，第，第n频率。频率。将任一将任一k代入特解得代入特解得此时各质点按同一频率此时各质点按同一频率k作同步简谐振动，各质点位移的比值为作同步简谐振动，各质点位移的比值为任何时刻结构的振动都保持同一形状。任何时刻结构的振动

50、都保持同一形状。主振动主振动多自由度结构按任一自振频率多自由度结构按任一自振频率k进行的简谐振动。进行的简谐振动。主振型主振型相应的特定振动形式，简称振型。相应的特定振动形式，简称振型。14-6 多自由度结构的自由振动将将k代回振幅方程得代回振幅方程得可写为可写为 系数行列式为零，系数行列式为零，n个方程中只有个方程中只有(n-1)个是独立的，不能确定各质个是独立的，不能确定各质点的幅值，但可确定其比值即振型。点的幅值，但可确定其比值即振型。14-6 多自由度结构的自由振动振型向量振型向量设设 ，即可求出其余各元素的值，此时振型称为标准化振型。，即可求出其余各元素的值，此时振型称为标准化振型。

51、主振动的线性组合构成振动微分方程的一般解：主振动的线性组合构成振动微分方程的一般解：各主振动分量的振幅、初相角各主振动分量的振幅、初相角由初始条件确定。由初始条件确定。自振频率、振型：与结构的质量分布和柔度系数有关；自振频率、振型：与结构的质量分布和柔度系数有关； 反映了结构本身固有的动力特性。反映了结构本身固有的动力特性。两个自由度体系的自由振动两个自由度体系的自由振动设设设设( ( ( (柔度矩阵）柔度矩阵）柔度矩阵）柔度矩阵）则则则则14-6 多自由度结构的自由振动14-6 多自由度结构的自由振动振幅方程为振幅方程为频率方程为频率方程为令令解得解得设解为设解为设解为设解为14-6 多自由

52、度结构的自由振动可得两个自振频率可得两个自振频率求第一阵型求第一阵型将将=1代入振幅方程可得代入振幅方程可得求第二阵型求第二阵型将将=2代入振幅方程可得代入振幅方程可得代入振幅方程，得振型方程代入振幅方程，得振型方程两个方程不独立，可由第一个方程解出两个主振型质点位移的比值两个方程不独立，可由第一个方程解出两个主振型质点位移的比值【例例】 试求结构的自振频率和振型试求结构的自振频率和振型. .1 1l l/4/41 1l l/2/2图图图图m m1 1= =m mm m2 2=2=2m ml l/2/2l l/2/2l l/2/2EIEI= =常数常数解解(1)(1)求柔度系数求柔度系数(2)

53、(2)求频率求频率14-6 多自由度结构的自由振动(3)(3)求振型求振型第一振型第一振型第二振型第二振型1 10.3050.3051 11.6391.63914-6 多自由度结构的自由振动【例例】 图图示示刚刚架架，在在梁梁跨跨中中D处处和和柱柱顶顶A处处有有大大小小相相等等的的集集中中质质量量m，支支座座C处为弹性支承，弹簧的刚性系数处为弹性支承，弹簧的刚性系数k=(3EI)/l3。试求自振频率和振型。试求自振频率和振型。 1. 求柔度系数求柔度系数 解：体系有两自由度，解：体系有两自由度，A处质点的水处质点的水平位移和平位移和D处质的竖向位移处质的竖向位移。 绘制绘制M1、M2图图, ,

54、由图乘及弹簧内力虚由图乘及弹簧内力虚功计算得功计算得14-6 多自由度结构的自由振动14-6 多自由度结构的自由振动2. 写出振型方程写出振型方程（a） 3. 写出频率方程，求频率写出频率方程，求频率展开式为展开式为解得解得 相应的频率为相应的频率为14-6 多自由度结构的自由振动第一主第一主振型振型第二主振型第二主振型2=2.9174. 求振型并绘出振型图求振型并绘出振型图由所得结果绘出振型由所得结果绘出振型1=27.083振型向量为振型向量为14-6 多自由度结构的自由振动例例14-3 试求图试求图a所示等截面简支梁的自振频率并确定主振型。所示等截面简支梁的自振频率并确定主振型。解：自由度

55、解：自由度=2，由图，由图b、c可得可得求得求得得到得到14-6 多自由度结构的自由振动第一阵型第一阵型第二阵型第二阵型如图如图d，振型是正对称的。，振型是正对称的。如图如图e，振型是反对称的。，振型是反对称的。结构的刚度和质量分布是对称的，结构的刚度和质量分布是对称的，则其主振型是正对称的或反对称的。则其主振型是正对称的或反对称的。取一半结构计算。取一半结构计算。14-6 多自由度结构的自由振动例例 计算图示结构的自振频率和主振型。计算图示结构的自振频率和主振型。EI=常数。常数。解：解：结构对称，振型分为对称和反对称。结构对称，振型分为对称和反对称。可取半结构计算，简化为单自由度计算可取半

56、结构计算，简化为单自由度计算对称振型对称振型反对称振型反对称振型第一振型第一振型第二振型第二振型14-6 多自由度结构的自由振动例例14-4 图图a所示刚架各杆所示刚架各杆EI都为常数，假设其质量集中于各结都为常数，假设其质量集中于各结 点处，点处，m2=1.5m1。试确定其自振频率和相应的振型。试确定其自振频率和相应的振型。解：结构是对称的，其振型为正、反对称两种。由受弯直杆的解：结构是对称的，其振型为正、反对称两种。由受弯直杆的 假定，判定不可能发生正对称形式的振动，其振型只能是假定，判定不可能发生正对称形式的振动，其振型只能是 反对称的。可取图反对称的。可取图b所示一半结构计算。所示一半

57、结构计算。超静定结构超静定结构14-6 多自由度结构的自由振动作超静定结构在作超静定结构在F1=1和和F2=1作用下的弯矩图，如图作用下的弯矩图，如图a、b。取静定的基本结构作取静定的基本结构作 图，如图图，如图c、d。计算得计算得14-6 多自由度结构的自由振动有有可得可得第一阵型第一阵型第二阵型第二阵型反对称反对称振动，振动，质点同质点同向振动向振动反对称反对称振动，振动，质点反质点反向振动向振动例题例题 试求结构的自振频率和振型试求结构的自振频率和振型. .EIEI= =常数常数m mm ml l/4/4l l/4/4l l/4/4l l/4/4m m1 13 3l l/16/161 1

58、l l/4/4图图图图1 13 3l l/16/16图图解解 (1)(1)求柔度系数求柔度系数14-6 多自由度结构的自由振动(2)求频率求频率14-6 多自由度结构的自由振动(3)(3)求振型求振型令每个振型的第一个元素为令每个振型的第一个元素为1 1，得，得1 11.411.414 41 1第三振型第三振型( (正对称正对称) )第二振型第二振型( (反对称反对称) )1 11 1第一振型第一振型( (正对称正对称) )1 11.411.414 41 114-6 多自由度结构的自由振动设解为设解为设解为设解为代入微分方程代入微分方程代入微分方程代入微分方程振型向量振型向量振型向量振型向量振

59、动时，方程必须有非零解，方程系数行列式为零：振动时，方程必须有非零解，方程系数行列式为零：振动时，方程必须有非零解，方程系数行列式为零：振动时，方程必须有非零解，方程系数行列式为零：频率方程频率方程由频率方程可解出自振频率，代回振幅方程得由频率方程可解出自振频率，代回振幅方程得振幅方程振幅方程确定相应主振型确定相应主振型2、按刚度法求解、按刚度法求解14-6 多自由度结构的自由振动或或或或14-6 多自由度结构的自由振动两个自由度两个自由度设特解为设特解为设特解为设特解为即即即即即即即即14-6 多自由度结构的自由振动代入运动方程，约去代入运动方程，约去代入运动方程，约去代入运动方程，约去si

60、n(sin( t t+ + ) )有非零解时，系数行列式等于零有非零解时，系数行列式等于零得振幅方程为得振幅方程为14-6 多自由度结构的自由振动两个自由度的结构频率方程为两个自由度的结构频率方程为展开展开解得解得两个方程不独立，可由第一个方程解出两个主振型质点位移比值为两个方程不独立，可由第一个方程解出两个主振型质点位移比值为代入振幅方程，得振型方程代入振幅方程，得振型方程14-6 多自由度结构的自由振动例例 计算图示结构的自振频率和主振型。计算图示结构的自振频率和主振型。解：解： 1、解释概念：层间侧移刚度、解释概念：层间侧移刚度层间侧移刚度：第层间侧移刚度：第i层之间发生单位侧移时水平力

61、。层之间发生单位侧移时水平力。用用ki 表示，也称为层间总剪表示，也称为层间总剪力。其值等于该层各柱的侧力。其值等于该层各柱的侧移刚度之和。移刚度之和。2、用层间刚度表示的刚度系数、用层间刚度表示的刚度系数14-6 多自由度结构的自由振动柱的侧移刚度：发生柱端单位侧移时的侧向力柱的侧移刚度：发生柱端单位侧移时的侧向力（即柱的剪力）。（即柱的剪力）。14-6 多自由度结构的自由振动两个主振型为两个主振型为14-6 多自由度结构的自由振动例例 计算图示结构的自振频率和主振型。已知横梁刚度无穷大，立柱计算图示结构的自振频率和主振型。已知横梁刚度无穷大，立柱刚度为刚度为EI=常数。常数。解：解：结构对

62、称，振型分为对称和反对称。结构对称，振型分为对称和反对称。可取半结构计算，简化为单自由度计算可取半结构计算，简化为单自由度计算请同学们自己绘出两个主振型。请同学们自己绘出两个主振型。例例14-5 图图a所示三层刚架横梁的刚度可视为无穷大，设刚架的所示三层刚架横梁的刚度可视为无穷大，设刚架的 质量集中在各层的横梁上。试确定其自振频率和主振型。质量集中在各层的横梁上。试确定其自振频率和主振型。14-6 多自由度结构的自由振动解：刚架振动时各横梁只能水平移动，自由度解：刚架振动时各横梁只能水平移动，自由度=3，结构的刚度，结构的刚度 系数如图系数如图b、c、d。14-6 多自由度结构的自由振动建立刚

63、度矩阵为建立刚度矩阵为质量矩阵为质量矩阵为14-6 多自由度结构的自由振动有有由频率方程得由频率方程得展开展开解得解得自振频率自振频率14-6 多自由度结构的自由振动确定主振型确定主振型将将k=1即即k=1=0.392代入振幅方程有代入振幅方程有设设标准化的第一振型为标准化的第一振型为同理可求得同理可求得14-6 多自由度结构的自由振动第一、二、三振型分别如图第一、二、三振型分别如图a、b、c。14-6 多自由度结构的自由振动主振型的正交性主振型的正交性 n个自由度的结构有个自由度的结构有n个自振频率及个自振频率及n个主振型，每一频率个主振型，每一频率及相应的主振型均满足振幅方程即：及相应的主

64、振型均满足振幅方程即：分别设分别设k=i，k=j，可得，可得两边左乘以两边左乘以两边左乘以两边左乘以则有则有K、M均为对称矩阵，将第二个式子两边转置有均为对称矩阵，将第二个式子两边转置有和和14-6 多自由度结构的自由振动将第一式减去第三式得将第一式减去第三式得当当ij时，时，i j，应有，应有此式表明对于质量矩阵此式表明对于质量矩阵M，不同频率的两个主振型是彼此正交的。，不同频率的两个主振型是彼此正交的。也可由上页的式子得也可由上页的式子得此式表明对于刚度矩阵此式表明对于刚度矩阵K，不同频率的两个主振型是彼此正交的。，不同频率的两个主振型是彼此正交的。 主振型的正交性是结构本身固有的特性，可

65、以用来简主振型的正交性是结构本身固有的特性，可以用来简化结构的动力计算，可用以检验所得主振型是否正确。化结构的动力计算，可用以检验所得主振型是否正确。14-6 多自由度结构的自由振动例例 验证图示结构主振型的正交性。验证图示结构主振型的正交性。解：前面已经解除了主振型，解：前面已经解除了主振型，质量矩阵质量矩阵刚度矩阵刚度矩阵表明关于表明关于M和和k都满足正交性，解答正确。都满足正交性，解答正确。14-6 多自由度结构的自由振动例例 两自由度结构的质量矩阵为两自由度结构的质量矩阵为两主振型中一个为两主振型中一个为求另一个主振型。求另一个主振型。设另一个主振型为设另一个主振型为由振型正交性由振型

66、正交性得得解得解得所求振型为所求振型为简谐荷载作用下的纯强迫振动简谐荷载作用下的纯强迫振动 图图(a)所示无重量简支梁，所示无重量简支梁，用柔度法建立振动微分方程。用柔度法建立振动微分方程。任一质点任一质点mi的位移的位移yi为为式中式中各动力荷载幅值在质点各动力荷载幅值在质点mi处引起的静力位移处引起的静力位移对对n个质点有个质点有14-7 多自由度结构在简谐荷载作用下的强迫振动运动微分方程的建立运动微分方程的建立1、柔度法、柔度法写成矩阵形式写成矩阵形式式中式中荷载幅值引起的静力位移向量荷载幅值引起的静力位移向量纯强迫振动的解答为纯强迫振动的解答为为质点为质点mi的振幅。的振幅。代入位移方

67、程可得代入位移方程可得振幅方程振幅方程14-7 多自由度结构在简谐荷载作用下的强迫振动或写为或写为式中式中I是单位矩阵，是单位矩阵，Y0是振幅向量。求解此方程即得各质点是振幅向量。求解此方程即得各质点在纯强迫振动中的振幅，从而得各质点的惯性力为在纯强迫振动中的振幅，从而得各质点的惯性力为惯性力的最大值惯性力的最大值结论：位移、惯性力、干扰力将同时达到最大值。结论：位移、惯性力、干扰力将同时达到最大值。 计算最大动力位移和内计算最大动力位移和内力时，可将惯性力、干扰力时，可将惯性力、干扰力的幅值作为静力荷载加力的幅值作为静力荷载加于结构上计算，如图于结构上计算，如图b。14-7 多自由度结构在简

68、谐荷载作用下的强迫振动将振幅方程改写为将振幅方程改写为可写为可写为最大惯性力向量最大惯性力向量当当=k (k=1,2,n)，振幅、惯性力、内力值均为无限大，振幅、惯性力、内力值均为无限大共振共振14-7 多自由度结构在简谐荷载作用下的强迫振动14-7 多自由度结构在简谐荷载作用下的强迫振动例例14-6 图图a为一等截面刚架，已知为一等截面刚架，已知m1=1kN， m2=0.5kN，F=5kN，每分钟振动，每分钟振动300次，次，l=4m， EI=5103kNm2。试作刚架的最大动力弯矩图。试作刚架的最大动力弯矩图。解：此对称刚架承受反对称荷载，可取图解：此对称刚架承受反对称荷载，可取图b所示半

69、刚架计算。所示半刚架计算。三个自由度：三个自由度：m1的水平位移的水平位移m2的水平位移的水平位移m3的竖向位移的竖向位移m1的最大惯性力的最大惯性力m2沿水平、竖向最大惯性力沿水平、竖向最大惯性力则有则有(1)14-7 多自由度结构在简谐荷载作用下的强迫振动求系数和自由项，作相应弯矩图如图求系数和自由项，作相应弯矩图如图cf。由图乘法得由图乘法得14-7 多自由度结构在简谐荷载作用下的强迫振动集中质量的数值为集中质量的数值为振动荷载的频率为振动荷载的频率为代入式代入式(1)得得解得解得由叠加法由叠加法最大动力弯矩图如图最大动力弯矩图如图g。14-7 多自由度结构在简谐荷载作用下的强迫振动 图

70、图a所示所示n个自由度个自由度的结构，当干扰力均作的结构，当干扰力均作用在质点处时，可得动用在质点处时，可得动力平衡方程为力平衡方程为写成矩阵形式写成矩阵形式若干扰力为同步简谐荷载若干扰力为同步简谐荷载式中式中F=( F1 F2 Fn )T，为荷载幅值列向量。，为荷载幅值列向量。14-7 多自由度结构在简谐荷载作用下的强迫振动在平稳阶段各质点均按频率在平稳阶段各质点均按频率作同步简谐振动。作同步简谐振动。代入动力平衡方程整理得代入动力平衡方程整理得求得各质点振幅值求得各质点振幅值各质点的惯性力为各质点的惯性力为可得可得求得惯性力幅值求得惯性力幅值 位移、惯性力、干扰力同时达到最大值，将位移、惯

71、性力、干扰力同时达到最大值，将FI、F(t)最大最大值作为静力荷载作用于结构，计算最大动力位移和内力。值作为静力荷载作用于结构，计算最大动力位移和内力。14-7 多自由度结构在简谐荷载作用下的强迫振动14-8 振型分解法多自由度结构无阻尼强迫振动微分方程为多自由度结构无阻尼强迫振动微分方程为只有集中质量的结构，只有集中质量的结构，M为对角阵，为对角阵，K不是对角阵不是对角阵方程藕联方程藕联各质点的位移向量各质点的位移向量几何坐标几何坐标坐标变换坐标变换结构标准化的主振型向量表示为结构标准化的主振型向量表示为设设位移向量按主振型分解位移向量按主振型分解展开展开14-8 振型分解法简写为简写为把几

72、何坐标把几何坐标Y变换成数目相同的另一组新坐标变换成数目相同的另一组新坐标正则坐标正则坐标主振型矩阵，几何坐标与正则坐标主振型矩阵，几何坐标与正则坐标 之间的转换矩阵之间的转换矩阵令令第第i个主振型的广义质量个主振型的广义质量广义质量矩阵，对角矩阵广义质量矩阵，对角矩阵14-8 振型分解法广义刚度矩阵，对角矩阵广义刚度矩阵，对角矩阵主对角线上的任一元素主对角线上的任一元素利用振型正交性可得利用振型正交性可得令令i=j，可得，可得或或与单自由度结构的频率公式相似与单自由度结构的频率公式相似14-8 振型分解法设设有有广义荷载向量广义荷载向量相应第相应第i个主振型的广义荷载个主振型的广义荷载振动方

73、程变换为振动方程变换为解除藕联，各自独立解除藕联，各自独立14-8 振型分解法整理得整理得与单自由度结构无阻尼强迫振动方程形式相同。与单自由度结构无阻尼强迫振动方程形式相同。初位移、初速度为零时，由杜哈梅积分求得初位移、初速度为零时，由杜哈梅积分求得n个自由度结构的计算简化为个自由度结构的计算简化为n个单自由度计算问题个单自由度计算问题振型分解法振型分解法（振型叠加法）：将位移（振型叠加法）：将位移Y分解为各主振型的叠加分解为各主振型的叠加14-8 振型分解法振型分解法计算步骤振型分解法计算步骤(1) 求自振频率和振型求自振频率和振型(2) 计算广义质量和广义荷载计算广义质量和广义荷载(3)

74、求解正则坐标的振动微分方程求解正则坐标的振动微分方程(4) 计算几何坐标计算几何坐标求出各质点位移求出各质点位移计算其他动力反应。计算其他动力反应。与单自由度问题一样求解。与单自由度问题一样求解。14-8 振型分解法例例14-7 图图a所示结构在结点所示结构在结点2处受有突加荷载作用，试求两处受有突加荷载作用，试求两 结点的位移和梁的弯矩。结点的位移和梁的弯矩。解：解：(1) 结构的自振频率和振型结构的自振频率和振型(图图b、c)(2) 广义质量广义质量14-8 振型分解法广义荷载广义荷载(3) 求正则坐标求正则坐标(4) 求位移求位移14-8 振型分解法两质点位移图形状如图两质点位移图形状如

75、图d。14-8 振型分解法(5) 求弯矩求弯矩两质点的惯性力为两质点的惯性力为由图由图e可求梁的动弯矩，如可求梁的动弯矩，如14-9 无限自由度结构的振动 图图a所示具有均布质量的单跨梁，其振动时弹性曲线上所示具有均布质量的单跨梁，其振动时弹性曲线上任一点的位移任一点的位移y是横坐标是横坐标x和时间和时间t的函数：的函数：设：梁的均布自重为设：梁的均布自重为q， 单位长度的质量单位长度的质量m=q/g， 惯性力的集度为惯性力的集度为取微段隔离体如图取微段隔离体如图b。由材料力学可得由材料力学可得14-9 无限自由度结构的振动如梁上承受均布简谐荷载如梁上承受均布简谐荷载psint，则梁的振动微分

76、方程为，则梁的振动微分方程为或或微分方程的解有两部分：相应齐次方程的一般解微分方程的解有两部分：相应齐次方程的一般解-梁的自由振动梁的自由振动 特解特解-梁的强迫振动梁的强迫振动(1) 梁的自由振动梁的自由振动微分方程为微分方程为设位移设位移y为坐标位置函数为坐标位置函数F(x)和时间函数和时间函数T(t)之积，即之积，即代入微分方程有代入微分方程有14-9 无限自由度结构的振动上式可写为上式可写为左边为变量左边为变量t的函数的函数右边为变量右边为变量x的函数的函数可设可设得得(1)(2)方程方程(1)的解为的解为令令或或频率特征值频率特征值式式(2)可写为可写为14-9 无限自由度结构的振动

77、上式通解为上式通解为位移为位移为振幅曲线为振幅曲线为A、B、C、D待定任意常数待定任意常数引入新的常量引入新的常量代入代入yx式中有式中有克雷洛夫函数克雷洛夫函数14-9 无限自由度结构的振动克雷洛夫函数有如下关系克雷洛夫函数有如下关系由这些关系可写出梁的挠度由这些关系可写出梁的挠度yx、角位移、角位移 、弯矩和剪力的公式、弯矩和剪力的公式(3)14-9 无限自由度结构的振动当当x=0时，设时，设有有可得可得全解为各特解的线性组合全解为各特解的线性组合(4)14-9 无限自由度结构的振动例例14-8 试求图试求图a所示等截面梁的自振频率和振型。所示等截面梁的自振频率和振型。解：由梁的边界条件，

78、解：由梁的边界条件，由式由式(4)可得可得系数行列式为零系数行列式为零展开展开化简为化简为14-9 无限自由度结构的振动由双曲函数和三角函数的图形可估计出由双曲函数和三角函数的图形可估计出试算法可求得前四个值为试算法可求得前四个值为相应的自振频率为相应的自振频率为可求得可求得由式由式(4)可得可得任意常数任意常数14-9 无限自由度结构的振动M0为待定值为待定值 将将k=k1，k2，分别代入分别代入yx可得出第一、第二、可得出第一、第二、 主振型主振型曲线，其形状如图曲线，其形状如图be。14-9 无限自由度结构的振动(2) 简谐均布干扰力作用下的振动简谐均布干扰力作用下的振动此时微分方程为此

79、时微分方程为设特解为设特解为代入上式有代入上式有令令可得方程可得方程14-10 计算频率的近似法(1) 能量法能量法 由能量守恒原理，结构在无阻尼自由振动时，动能由能量守恒原理，结构在无阻尼自由振动时，动能T和和应变能应变能V之和应为常数，即之和应为常数，即应有应有即即设梁的振动方程为设梁的振动方程为速度为速度为动能为动能为14-10 计算频率的近似法应变能为应变能为由由Tmax=Vmax得得如结构上还有集中质量如结构上还有集中质量mi（i=1，2，n），上式为），上式为14-10 计算频率的近似法 计算时，通常采用结构自重作用下的弹性曲线作为计算时，通常采用结构自重作用下的弹性曲线作为y(x

80、)，此时应变能可用外力功来代替，即此时应变能可用外力功来代替，即频率计算公式改写为频率计算公式改写为如求水平方向振动的频率，则重力应沿水平方向作用。如求水平方向振动的频率，则重力应沿水平方向作用。14-10 计算频率的近似法例例14-9 试用能量法求图试用能量法求图a所示等截面梁的第一自振频率。所示等截面梁的第一自振频率。解：取梁在自重解：取梁在自重q作用下的挠曲线作用下的挠曲线 作为第一振型，如图作为第一振型，如图b，即，即q=mg，因而，因而代入公式得代入公式得精确值为精确值为14-10 计算频率的近似法例例14-10 试能量法求图试能量法求图a所示刚架的最低自振频率。所示刚架的最低自振频

81、率。解：将各层重量解：将各层重量mi g作为水平力加于结构如图作为水平力加于结构如图b，此时位移作，此时位移作 为第一振型。为第一振型。14-10 计算频率的近似法一般说，一般说，n层刚架中第层刚架中第i层位移为层位移为代入公式得代入公式得精确值精确值14-10 计算频率的近似法(2) 集中质量法集中质量法例例14-11 试求图试求图a所示具有均布质量所示具有均布质量m的简支梁的自振频率。的简支梁的自振频率。解解：(1) 将梁分为两段，将每段将梁分为两段，将每段的的 质量集中于该段的两端，质量集中于该段的两端， 梁化为单自由度结构。梁化为单自由度结构。精确解精确解(2) 如求第一、第二频率，至

82、少把结构化为有两个自由度如图如求第一、第二频率，至少把结构化为有两个自由度如图b。 精确解精确解14-10 计算频率的近似法(3) 如求第一、第二、第三频率，至少把结构化为有三个如求第一、第二、第三频率，至少把结构化为有三个自自 由度如图由度如图c。 精确解精确解结论：集中质量法能给出较好的近似结果，在工程上结论：集中质量法能给出较好的近似结果，在工程上 常被采用。常被采用。(3) 用相当梁法计算桁架的最低频率用相当梁法计算桁架的最低频率14-10 计算频率的近似法相当梁相当梁：一个在某一特征点处位移与桁架位移相等的梁。：一个在某一特征点处位移与桁架位移相等的梁。当桁架变形时，任一结点当桁架变

83、形时，任一结点k的竖向位移的竖向位移kP为为相当梁同一点相当梁同一点k的竖向位移的竖向位移vkP将是其惯性矩将是其惯性矩I的函数，可写为的函数，可写为二者相等即可算出相当梁的惯性矩二者相等即可算出相当梁的惯性矩I。14-10 计算频率的近似法 简支桁架的自重为简支桁架的自重为q，则在具有相同重量的相当梁重点，则在具有相同重量的相当梁重点的竖向位移为的竖向位移为令令可按简支梁的频率公式求其最低频率可按简支梁的频率公式求其最低频率将将m=q/g和和EI代入得代入得14-10 计算频率的近似法 如图所示对称桁架，设如图所示对称桁架，设m=10kg，E=200GPa，在所示，在所示质量的重力作用下，结点质量的重力作用下，结点3的竖向位移可求得的竖向位移可求得代入频率公式可得代入频率公式可得