1、第14章 结构动力学,动力计算概述 单自由度体系的自由振动 单自由度体系的强迫振动 多自由度体系的自由振动 多自由度体系的强迫振动 频率的近似计算,知识点,教学基本要求,了解结构动力计算的特点，能够判断动力计算自由度； 掌握单体系振动微分方程的建立方法。,掌握单自由度体系在不同的动荷载作用下强迫振动 的分析方法以及动力特性。掌握阻尼对单自由度体 系动力特性的影响。,理解柔度法和刚度法建立振动微分方程的思路。掌 握两个自由度体系的频率方程和自振频率的求解， 理解主振型和主振型正交性，掌握振型分解法。,了解计算频率的几种近似法,能够正确计算单自由度体系的固有频率和周期。,14.1 动力计算概述,1
2、、动力计算的特点、目的和内容,1)特点：静力荷载与动力荷载的特点及其效应。,静力荷载是指其大小、方向和作用位置不随时间而变化的荷载。这 类荷载对结构产生的惯性力可以忽略不计，由它所引起的内力和变形都 是确定的。,动力荷载是指其大小、方向和作用位置随时间而变化的荷载。这类 荷载对结构产生的惯性力不能忽略，因动力荷载将使结构产生相当大的加 速度，由它所引起的内力和变形都是时间的函数。,与静力计算的对比：两者都是建立平衡方程，但动力计算，利用动静法，建立的是形式上的平衡方程。力系中包含了惯性力，考虑的是瞬间平衡，荷载、内力都是时间的函数。建立的平衡方程是微分方程。,2)目的和内容,计算结构的动力反应
3、：内力、位移、速度与加速度，使结构在动内力 与静内力共同作用下满足强度和变形的要求。,动力计算的内容：研究结构在动荷载作用下的动力反应的计算原理和 方法。,涉及到内外两方面的因素：(1）确定动力荷载（外部因素，即干扰力）；(2）确定结构的动力特性（内部因素，如结构的自振频率、周期、振型和阻尼等等），类似静力学中的I、S等；,计算动位移及其幅值；计算动内力及其幅值。,简谐荷载（按正余弦规律变化）,一般周期荷载,2、动力荷载分类 按起变化规律及其作用特点可分为：1）周期荷载：随时间作周期性变化。（转动电机的偏心力）,2）冲击荷载：短时内剧增或剧减。（如爆炸荷载）,tr,P,tr,P,3、动力计算中
4、体系的自由度确定体系上全部质量位置所需独立参数的个数称为体系的振动自由度。,实际结构的质量都是连续分布的，严格地说来都是无限自由度体系。 计算困难，常作简化如下：1)集中质量法把连续分布的质量集中为几个质点，将一个无限自由度的问题简化成有 限自由度问题。,3）随机荷载：(非确定性荷载) 荷载在将来任一时刻的数值无法事先确定。（如地震荷载、风荷载）,2个自由度,2个自由度,自由度与质量数不一定相等,m,mm梁,m,+m梁,I,I,2I,m,+m柱,厂房排架水平振时的计算简图,单自由度体系,水平振动时的计算体系,多自由度体系,构架式基础顶板简化成刚性块,(t),v(t),u(t),4个自由度,m1
5、,m2,m3,2个自由度,无限自由度体系,2)广义座标法：,如简支梁的变形曲线可用三角级数来表示,用几条函数曲线来描述体系的振动曲线就称它是几个自由度体系，其中, 是根据边界约束条件选取的函数，称为形状函数。,ak(t) 称广义座标，为一组待定参数，其个数即为自由度数，用此法可将无限自由度体系简化为有限自由度体系。,a1， a2， an,4、动力计算的方法,动力平衡法（达朗伯尔原理）,运动方程,设其中,P(t),I(t),平衡方程,I(t)惯性力，与加速度成正比，方向相反,改写成,虚功原理（拉格朗日方程）,哈米顿原理（变分方程）,都要用到抽象的虚位移概念,自由振动：体系在振动过程中没有动荷载的
6、作用。,自由振动产生原因：体系在初始时刻（t=0）受到外界的干扰。,研究单自由度体系的自由振动重要性在于：,1)它代表了许多实际工程问题，如水塔、单层厂房等。,2)它是分析多自由度体系的基础，包含了许多基本概念。,自由振动反映了体系的固有动力特性。,要解决的问题包括：,建立运动方程、计算自振频率、周期和阻尼.,14.2 单自由度体系的自由振动,1、运动微分方程的建立,方法：达朗伯尔原理,应用条件：微幅振动（线性微分方程）,1)刚度法：研究作用于被隔离的质量上的力，建立平衡方程。,m,yj,yd,质量m在任一时刻的位移 y(t)=yj+yd,k,力学模型,yd,m,m,W,I(t),重力 W,弹
7、性力,恒与位移反向,惯性力,(a),其中 kyj=W,上式可以简化为,或,由平衡位置计量。以位移为未知量的平衡方程式，引用了刚度系数，称 刚度法。,(a),2)柔度法：研究结构上质点的位移，建立位移协调方程。,可得与 (b) 相同的方程,刚度法常用于刚架类结构，柔度法常用于梁式结构。,2、自由振动微分方程的解,改写为,其中,它是二阶线性齐次微分方程，其一般解为：,积分常数C1，C2由初始条件确定,设 t=0 时,（d）式可以写成,由上式可知，位移是由初位移y引起的余弦运动和由初速度v引起的正弦运动的合成.,由上式可知，位移是由初位移y引起的余弦运动和由初速度v引起的正弦运动的合成，为了便于研究
8、合成运动,令,(e)式改写成,它表示合成运动仍是一个简谐运动。其中A和可由下式确定,振幅,相位角,3、结构的自振周期和频率,由式,及图可见位移方程是一个周期函数。,周期,工程频率,园频率,计算频率和周期的几种形式,其中 是沿质点振动方向的结构柔度系数，它表示在质点上沿振动方向加单位荷载使质点沿振动方向所产生的位移。 k使质点沿振动方向发生单位位移时，须在质点上沿振动方向施加的力。st=W在质点上沿振动方向施加数值为W的荷载时质点沿振动方向所产生的位移。 计算时可根据体系的具体情况，视、 k、 st 三参数中哪一个最便于计算来选用。,一些重要性质： （1）自振周期与且只与结构的质量和结构的刚度有
9、关，与外界的干扰因素无关。干扰力只影响振幅。 （2）自振周期与质量的平方根成正比，质量越大，周期越大（频率越小）；自振周期与刚度的平方根成反比，刚度越大，周期越小（频率越大）；要改变结构的自振周期，只有从改变结构的质量或刚度着手。 （3）两个外形相似的结构，如果周期相差悬殊，则动力性能相差很大。反之，两个外形看来并不相同的结构，如果其自振周期相近，则在动荷载作用下的动力性能基本一致,是结构动力特性的重要数量标志。,例1. 计算图示结构的频率和周期。,例2.计算图示结构的水平和竖向振动频率。,例3.计算图示刚架的频率和周期。,由截面平衡,例4、图示三根单跨梁，EI为常数，在梁中点有集中质量m，不
10、考虑梁的质量，试比较三者的自振频率。,解：1）求,3l/16,5l/32,l/2,据此可得：1 2 3= 1 1.512 2,结构约束越强,其刚度越大,刚度越大,其自振动频率也越大。,例5、求图示结构的自振圆频率。,解法1：求 k,=1/h,MBA=kh = MBC,解法2：求 ,例6、求图示结构的自振频率。,解：求 k,对于静定结构一般计算柔度系数方便。 如果让振动体系沿振动方向发生单位位移时，所有刚节点 都不能发生转动（如横梁刚度为刚架)计算刚度系数方便。,一端铰结的杆的侧移刚度为：,两端刚结的杆的侧移刚度为：,4、简谐自由振动的特性,由式,可得，加速度为：,在无阻尼自由振动中，位移、加速
11、度和惯性力都按正弦规律变化，且作相位相同的同步运动，即它们在同一时刻均达极值，而且惯性力的方向与位移的方向一致。,它们的幅值产生于,时，其值分别为：,既然在运动的任一瞬时质体都处于平衡状态，在幅值出现时间也一样，于是可在幅值处建立运动方程，此时方程中将不含时间t，结果把微分方程转化为代数方程了，使计算得以简化。,惯性力为：,例7. 计算图示体系的自振频率。,解：单自由度体系，以表示位移参数的幅值, 各质点上所受的力为：,建立力矩平衡方程,化简后得,5、阻尼对振动的影响,实验证明，振动中的结构，不仅产生与变形成比例的弹性内力，还产生非弹性的内力，非弹性力起阻尼作用。在不考虑阻尼的情况下所得出的某
12、些结论也反应了结构的振动规律，如：,事实上，由于非弹性力的存在，自由振动会衰减直到停止；共振时振幅也不会无限增大，而是一个有限值。非弹性力起着减小振幅的作用，使振动衰减，因此，为了进一步了解结构的振动规律，就要研究阻尼。,1)阻尼的存在,忽略阻尼的振动规律,考虑阻尼的振动规律,结构的自振频率是结构的固有特性，与外因无关。,简谐荷载作用下有可能出现共振。,自由振动的振幅永不衰减。,自由振动的振幅逐渐衰减。,共振时的振幅趋于无穷大。,共振时的振幅较大但为有限值。,2)在建筑物中产生阻尼、耗散能量的因素,(1）结构在变形过程中材料内部有摩擦，称“内摩擦”，耗散能量；,(2）建筑物基础的振动引起土壤发
13、生振动，此振动以波的形式向周围扩散，振动波在土壤中传播而耗散能量；,(3）土体内摩擦、支座上的摩擦、结点上的摩擦和空气阻尼等等。,振动的衰减和能量的耗散都通过非弹性力来考虑，由于对非弹性力的描述不同，目前主要有两种阻尼理论：,*粘滞阻尼理论非弹性力与变形速度成正比：,*滞变阻尼理论,关于阻尼，有两种定义或理解：,(1）使振动衰减的作用；,(2）使能量耗散。,3)阻尼力的确定：总与质点速度反向；大小与质点速度有如下关系：(1）与质点速度成正比（比较常用，称为粘滞阻尼）。(2）与质点速度平方成正比（如质点在流体中运动受到的阻力）。(3）与质点速度无关（如摩擦力）。,其他阻尼力也可化为等效粘滞阻尼力
14、来分析。,k,m,4)有阻尼的自由振动，动平衡方程：,( 阻尼比）,设解为：,特征方程为：,(1)1（低阻尼）情况,令,ae-t,低阻尼y- t曲线,无阻尼y- t曲线,阻尼对自振频率的影响.,当0.2，则存在0.96r/1。在工程结构问题中，0.010.1，可近似取:,称为振幅的对数递减率.,设yk和yk+n是相隔n个周期的两个振幅则:,经过一个周期后，相邻两振幅yk和yk+1的比值的对数为:,工程中常用此方法测定阻尼,阻尼对振幅的影响.振幅ae- t 随时间衰减，相邻两个振幅的比,振幅按等比级数递减.,例8、图示一单层建筑物的计算简图。屋盖系统和柱子的质量均集中在横梁处共计为m,，加一水平
15、力P=9.8kN，测得侧移A0=0.5cm，然后突然卸载使结构发生水平自由振动。在测得周期T=1.5s 及一个周期后的侧移A1=0.4cm。求结构的阻尼比和阻尼系数c。,解：,临界阻尼常数cr为=1时的阻尼常数。（振与不振的分界点）,阻尼比,反映阻尼情况的基本参数。,(3)1 强阻尼：不出现振动，实际问题不常见。,(2)=1（临界阻尼）情况,这条曲线仍具有衰减性， 但不具有波动性。,受迫振动（强迫振动）：结构在动力荷载作用下的振动。,k,弹性力ky、惯性力,和荷载P(t)之间的平衡方程为:,1、简谐荷载：,单自由度体系强迫 振动的微分方程,特解：,14.3 单自由度体系的受迫振动,最大静位移y
16、st（是把荷载幅值当作静荷载作用时结构所产生 的位移）。,特解可写为：,通解可写为：,设t=0时的初始位移和初始速度均为零，则：,过渡阶段：振动开始两种振动同时存在的阶段； 平稳阶段：后来只按荷载频率振动的阶段。（由于阻尼的存在）,按自振频率振动,按荷载频率振动,平稳阶段：,最大动位移（振幅）为：,动力系数为：,重要的特性： 当/0时,1，荷载变化得很慢，可当作静荷载处理。 当01，并且随/的增大而增大。 当/ 1时,。即当荷载频率接近于自振频率时，振幅会无限增大。称为“共振”。通常把0.75 / 1.25称为共振区。,当/ 1时,的绝对值随/ 的增大而减小。当很大时，荷载变化很快，结构来不及
17、反应。,当动荷载作用在单自由度体系的质点上时，由于体系上各 截面的内力、位移都与质点处的位移成正比，故各截面的最大动内力和最大动位移可采用统一的动力系数，只需将干扰力幅值乘以动力系数按静力方法来计算即可。,例11：已知m=300kg，EI=90105N.m2 ,k=48EI/l3 ,P=20kN,=80s-1求梁中点的位移幅值及最大动力弯矩。,解：1)求,2)求,3)求ymax， Mmax,如何求最大位移和最大弯矩？,动荷载不作用于质点时的计算,令,仍是位移动力系数,是内力动力系数吗?,运动方程,稳态解,振幅,列幅值方程求内力幅值,解:,例12求图示体系振幅、动弯矩幅值图.,动弯矩幅值图,解:
18、,例13求图示体系右端的质点振幅,o,设体系在t=0时静止，然后有瞬时冲量S作用。,2、一般荷载,一般荷载作用下的动力反应可利用瞬时冲量的动力反应来推导,1)瞬时冲量的动力反应,瞬时冲量S引起的振动可视为由初始条件引起的自由振动。 由动量定理：,2)任意荷载P(t)的动力反应,时刻的微分冲量对t瞬时(t )引起的动力反应:,初始静止状态的单自由度体系在任意荷载作用下的位移公式:,(Duhamel 积分),初始位移y0和初始速度v0不为零在任意荷载作用下的位移公式:,3)几种典型荷载的动力反应,(1）突加荷载,yst=P0=P0 /m2,质点围绕静力平衡位置作简谐振动,(2）短时荷载,阶段(0t
19、u)：与突加荷载相同。,阶段(tu)：无荷载，体系以t=u时刻的位移,和速度,为初始条件作自由振动。,或者直接由Duhamel积分作,另解：短时荷载可认为由两个突加荷载叠加而成。,当0t u,当t u,最大动反应,1）当 u T/2 最大动位移发生在阶段,2）当u T/2 最大动位移发生在阶段, =2,动力系数反应谱 (与T和u之间的关系曲线),(3）线性渐增荷载,这种荷载引起的动力反应同样可由Duhamel积分来求解:,对于这种线性渐增荷载,其动力反应与升载时间的长短有很大关系。其动力系数的反应谱如下：,动力系数反应谱,动力系数介于1与2之间。 如果升载很短，tr4T,则接近于1,即相当于静
20、荷载情况。 常取外包虚线作为设计的依据。,3、有阻尼的强迫振动,单独由v0引起的自由振动：,瞬时冲量ds=Pdt=mv0所引起的振动，可视为以v0=Pdt/m，y0=0为初始条件的自由振动：,将荷载P(t)的加载过程 看作一系列瞬时冲量：,总反应,（1）突加荷载P0,低阻尼y- t曲线,无阻尼y- t曲线,静力平衡位置,具有阻尼的体系在 突加荷载作用下， 最初所引起的最大 位移接近于静位移 yst=P0/m2的两倍， 然后逐渐衰减，最 后停留在静力平衡 位置。,（2）简谐荷载P(t)=Fsint,设特解为：y=Asin t +Bcos t代入得：,+Asin t +Bcos t ,齐次解加特解
21、得到通解：,自由振动，因阻尼作用， 逐渐衰减、消失。,纯强迫振动，平稳振动， 振幅和周期不随时间而变化。,结论：在简谐荷载作用下，无论是否计入阻尼的作用，纯强迫振动部分总是稳定的周期运动，称为平稳振动。,y=Asin t +Bcos t =yPsin(t ),振幅：yp， 最大静力位移：yst=F/k=F/m2,动力系数与频率比/和阻尼比有关,几点注意： 随增大曲线渐趋平缓，特别是在/=1附近的峰值下降的最为显著。,当接近 时， 增加很快， 对的数值影响也很大。在0.75 / 1.25(共振区)内，阻尼大大减小了受迫振动的位移，因此, 为了研究共振时的动力反映, 阻尼的影响是不容忽略。在共振区
22、之外阻尼对的影响较小，可按无阻尼计算。,max并不发生在共振/=1时，而发生在，,由y=yPsin(t ) 可见，阻尼体系的位移比荷载P=Fsin t 滞后一个相位角 ，,但因很小，可近似地认为：,当时,0体系振动得很慢，FI、R较小，动荷主要由S平衡，S与y反向，y与P基本上同步；荷载可作静荷载处理。,当时,180体系振动得很快，FI很大，S、R相对说来较小，动荷主要由FI 平衡， FI 与y同向，y与P反向；,弹性力S，惯性力FI， 阻尼力R分别为：,当=时,90,由此可见：共振时（=），S与FI刚好互相平衡，,yst,有无阻尼均如此。动荷恰与阻尼力平衡，故运动呈现稳态故不会出现内力为无穷
23、大的情况。而在无阻尼受迫振动时，因不存在阻尼力与动荷载平衡，才出现位移为无限大的现象。,k=m2=m2,例14 图示块式基础.机器与基础的质量为 ;地基竖向刚度为 ;竖向振动时的阻尼比为 机器转速为N=800r/min,其偏心质量引起的离心力为P=30kN.求竖向振动时的振幅。,解：,例15 求突加荷载作用下的位移，开始时静止，不计阻尼。,解：,动力系数为 2,m1,m2,14.4 两自由度体系的自由振动,自由振动分析的目的是确定体系的动力特性.可不计阻尼。,1.运动方程及其解,或,设方程的特解为,代入方程,得,-频率方程,振型方程,显然 Y1=Y2=0 为其解，为了求得不全为零的解，令,解频
24、率方程得 的两个根,将 频率代入振型方程,特解2,特解1,通解,2.频率与振型,m1,m2,体系按特解振动时有如下特点,1)各质点同频同步;,2)任意时刻,各质点位移的比值保持不变,几点说明：,1)按振型作自由振动时，各质点的速度的比值也为常数，且与位移比值相同。,2)发生按振型的自由振动是有条件的.,4)N自由度体系有N个频率和N个振型,频率方程,解频率方程得 ,从小到大排列,依次称作第一频率,第二频率.,第一频率称作基本频率,其它为高阶频率.,将频率代入振型方程,得N个振型,N个振型是线性无关的.,3)振型与频率是体系本身固有的属性,与外界因素无关.,5)柔度法,设解为,此时惯性力,幅值,
25、在自由振动过程中任意时刻t，质量m1、m2的位移y1(t)、y2(t)应当等于体系在当时惯性力作用下的静力位移。,当然解 Y1=Y2=0， 为了求得不全为零的解，令,令,主振型,振型可看作是体系按振型振动时， 惯性力幅值作为静荷载所引起的静位移,采用柔度矩阵时,6)求振型、频率可列幅值方程.,振型方程,频率方程,按振型振动时,例16 设图示刚架横梁刚度为无限大，层间侧移刚度分别为k1和 k2 ，试求刚架水平振动时的自振动频率和主振型。,m1,m2,k1,k2,解：（1）求频率方程中的刚度系数,k11=k1+k2,k12=k21=-k2,k22=k2,（3）求主振型,第1振型,第2振型,（2）求
26、频率,代公式,若有,（5）求主振型,（4）求频率,若有,若 n=90,则第一振型和第二振型分别为：,可见当顶端质点的质量和刚度很小时，顶端水平侧移很大。,建筑结构抗震设计中，将这种因顶端质点质量和刚度突变，而导致顶端巨大反应的现象，称为鞭梢效应。,如：屋顶消防水池、上人屋面设计的楼电梯间，女儿墙或屋顶建筑物等。,例17 质量集中在楼层上，,层间侧移刚度如图。求自振频率,解：1）求刚度系数：,k,k33=k/5,刚度矩阵K和质量矩阵M：,展开得：234222252250 解得：1=1.293， 2=6.680， 3=13.027,2）求频率：代入频率方程： K2 M0,3）求主振型：振型方程：（
27、K2 M）Y0的后两式：（令Y3i=1）,（a）,Yij为正时表示质量mi的运动方向与单位位移方向相同，为负时，表示与单位位移方向相反。,例18 试求图示梁的自振频率和主振型，梁的EI已知。,解：（1）计算频率,（2）振型,第一振型,第二振型,3、主振型及主振型的正交性,由功的互等定理：,整理得：,因 ，则存在：,两个主振型相互正交，因与质量有关，称为第一正交关系。,由功的互等定理：,上式分别乘以12、22，则得：,第一主振型惯性力在第二主振型位移上所做的功等于零；第二主振型惯性力在第一主振型位移上所做的功等于零；,某一主振型的惯性力在其它主振型位移上不做功，其能量不会转移到其它主振型上，不会
28、引起其它主振型的振动；,各个主振型能单独存在，而不相互干扰。,1、柔度法（忽略阻尼）因为在简谐荷载作用下， 荷载频率在共振区之外，阻尼 影响很小；在共振区之内时， 计不计阻尼，虽对振幅影响很大，但都能反映共振现象。,（2）动位移的解答及讨论 通解包含两部分：齐次解对应按自振频率振动的自由振动， 由于阻尼而很快消失；特解对应按荷载频率振动的简谐振 动是平稳阶段的纯强迫振动。,（1）建立振动微分方程,各简谐荷载频率相同相位相同，否则用其他方法,14.5 两自由度体系的受迫振动,n各自由度体系，存在n个可能的共振点,设纯强迫振动解答为：,代入：,（3）动内力幅值的计算,荷载、位移、惯性力同频、同相、
29、同时达到最大。位移达到最大时，内力也达到最大。求内力时可将动荷载和惯性力的幅值作为静荷载作用于结构，用静力法求出内力，即为动内力幅值。或用叠加公式求：,由Y1 ，Y2值可求得位移和惯性力。,惯性力的幅值为：,代入位移幅值方程,可得求惯性力幅值的方程 （直接求惯性力幅值）,在平稳阶段，各质点也作简谐振动：,Y1=D1/D0,Y2=D2/D0,如果荷载频率与任一个自振频率 1、 2重合，则D0=0, 当D1、D2 不全为零时，则出现共振现象,2.刚度法,例19 质量集中在楼层上m1、m2 ，层间侧移刚度为k1、k2,解：荷载幅值：P1=P，P2=0，求刚度系数：,k11=k1+k2 , k21=k
30、2 ,k22=k2 , k12=k2,当m1=m2=m，k1=k2=k,两个质点的 位移动力系 数不同。,当,趋于无穷大。 可见在两个自由度体系中，在两种情况下可能出现共振。也有例外情况,yst1,yst2=P/k,荷载幅值产生的静位移和静内力,yst1= yst2=P/k,层间剪力: Qst1= P,动荷载产生的位移幅值和内力幅值,由此可见，在多自由度体系中，没有一个统一的动力系数。,层间动剪力:,这说明在右图结构上，,适当增加m2、k2系统,可以消除m1的振动（动力吸振器原理）。,吸振器不能盲目设置，必须在干扰力使体系产生较大振动时才有必要设置。,设计吸振器时，先根据m2的许可振幅Y2，选
31、定,，再确定,例21 如图示梁中点放一电动机。重2500N，电动机使梁中点产生的静位移为1cm，转速为300r/min，产生的动荷载幅值P=1kN，问：1）应加动力吸振器吗？2）设计吸振器。(许可位移为1cm),解：1）,频率比在共振区之内应设置吸振器。,2）由,弹簧刚度系数为：,N/m,=102 kg,例22 图示简支梁EI=常数，=0.751求动位移幅值和动弯矩幅值。,解：1)求柔度系数,2)作MP图，求1P 2P,5)计算动内力,1.4119P,0.2689P,0.8740P,Qd 图,0.3530Pl,0.2180Pl,Md 图,6)比较动力系数,因此，多自由度体系没有统一的动力系数。
32、,如图示对称结构在对称荷载作用下。,与2相应的振型是,=1,当=2 ，D0=0 ，也有：,不会趋于无穷大，不发生共振， 共振区只有一个。,对称体系在对称荷载作用下时，只有当荷载频率与对称主振型的自振频率相等时才发生共振；当荷载频率与反对称主振型的自振频率相等时不会发生共振。同理可知：对称体系在反对称荷载作用下时，只有当荷载频率与反对称主振型的自振频率相等时才发生共振。,例23 已知图a刚架受简谐荷载作用，=0.6,绘出动力弯矩图Md，并求柱顶最大位移 ymax。,解：利用对称性取半边结构如图所示。,柱顶位移,只考虑稳态振动，设方程的特解,代入方程解得,，,例24 求图a所示体系的自振频率及主振
33、型。梁EI =常数。,解 ：将原结构化成正对称和反对称半结构分别计算（图b、c）。,，,当=1时，振型为正对称，则,当=2时，振型为反对称，则,1、能量法求第一频率Rayleigh法,根据能量守恒定律，当不考虑阻尼自由振动时，振动体系在任何时刻的动能T 和应变能U 之和应等于常数。 根据简谐振动的特点可知：在体系通过静力平衡位置的瞬间，速度最大（动能具有最大值），动位移为零（应变能为零）；当体系达到最大振幅的瞬间（变形能最大），速度为零（动能为零）。对这两个特定时刻，根据能量守恒定律得： Umax=Tmax,求Umax ，Tmax,求频率,如梁上还有集中质量mi，,位移幅值,.,Yi为集中质量
34、mi处的位移幅值。,14.6 近似法求自振频率,假设位移幅值函数Y(x)必须注意以下几点：,1)必须满足运动边界条件：（铰支端：Y=0；固定端：Y=0，Y=0）,尽量满足弯矩边界条件，以减小误差。剪力边界条件可不计。 2)所设位移幅值函数应与实际振型形状大致接近；如正好与第 n 主振型相似，则可求的n的准确解。但主振型通常是未知的，只能假定一近似的振型曲线，得到频率的近似值。由于假定高频率的振型困难，计算高频率误差较大。故 Rayleigh法主要用于求1的近似解。 3)相应于第一频率所设的振型曲线，应当是结构比较容易出现的变形形式。曲率小，拐点少。 4)通常可取结构在某个静荷载q（x）（如自重
35、）作用下的弹性曲线作为Y（x）的近似表达式。此时应变能可用相应荷载q（x）所作的功来代替，即,2）假设均布荷载q作用下的挠度曲线作为Y(x),例25 试求等截面简支梁的第一频率。1）假设位移形状函数为抛物线,满足边界条件且与第一振型相近,3）假设,第一振型的精确解。,精 确 解,例26 求楔形悬臂梁的自振频率。 设梁截面宽度为 1，高度为 h=h0x/l。,解：,单位长度的质量：,设位移形状函数：,满足边界条件：,Rayleigh 法所得频率的近似解总是比精确解偏高。其原因是假设了一振型曲线代替实际振型曲线，迫使梁按照这种假设的形状振动，相当于给梁加上了某种约束，增大了梁的刚度，致使频率偏高。
36、当所设振型越接近于真实，则相当于对体系施加的约束越小，求得的频率越接近于真实，即偏高量越小。,截面惯性矩：,1)假设多个近似振型,都满足前述两个条件。,2)将它们进行线性组合,（a1、a2、an是待定常数）,3)确定待定常数的准则是：获得最佳的线性组合，这样的Y（x）代入频率计算公式中得到的2 的值虽仍比精确解偏高，但对所有的a1，a2，an的可能组合，确实获得了最小的2值。,所选的a1，a2，an使,2 获得最小值的条件是,这是以a1，a2，an为未知量的n个奇次线性代数方程。令其系数行列式等于零，得到频率方程，可以解出原体系最低 n 阶频率来。阶次越低往往越准。,为了使假设的振型尽可能的接
37、近真实振型，尽可能减小假设振型对体系所附加的约束， Ritz 提出了改进方法：,例27 用RayleighRitz 法求等截面悬臂梁的最初几个频率。,解：悬臂梁的位移边界条件为：,只取第一项,代入：,代入频 率方程：,其精确解：,与精确解相比，误差为27%。,例28 用RayleighRitz法求等截面悬臂梁的最初几个频率。,解：,取两项,代入：,代入频率方程：,求得kij，mij：,求得最 初两个 频率近 似值：,（0.48%）,（58%） 说明,说明：1)由于1、2均近似于第一振型，由它们组合的第二振型自然很差，故第二频率不准。2)RayleighRitz法所得结果仍然偏高，其原因同瑞利法
38、。,2、集中质量法,在计算无限自由度体系的自振频率时，可以用若干个集中质量来代替连续分布的质量。关于质量的集中方法有多种，最简单的是静力等效的集中质量法。,该法既可求基本频率，也可求较高频率。且适用于各类结构。,集中质量的数目越多结果越精确，但工作量也就越大。,等效原则：使集中后的重力与原重力互为静力等效，即两者的合力相等。 作 法：将杆分为若干段，将每段质量集中于其质心或集中于两端。,例29 试用集中质量法求简支梁自振频率。,(0.7%）,(0.1%）,(3.1%）,(0.05%）,(4.8%）,(0.7%）,对于对称刚架，可分别用不同的集中质量方案求出对称振动和反对称振动的自振频率。,最小频率对应 着反对称振型,
