结构动力学读书笔记.doc

《结构动力学读书笔记.doc》由会员分享，可在线阅读，更多相关《结构动力学读书笔记.doc（10页珍藏版）》请在三一文库上搜索。

1、结构动力学读书报告 学院专业学号姓名指导老师2013 年 5月 28日摘要：本书在介绍基本概念和基础理论的同时,也介绍了结构动力学领域的若干前沿研究课题。既注重读者对基本知识的掌握,也注重读者对结构振动领域研究发展方向的掌握。主要内容包括运动方程的建立、单自由度体系、多自由度体系、无限自由度体系的动力学问题、随机振动、结构动力学的前沿研究课题。侧重介绍单自由度体系和多自由度体系,重点突出,同时也着重介绍了在抗震中的应用。1 概述1.1结构动力学的发展及其研究内容：结构动力学,作为一门课程也可称作振动力学,广泛地应用于工程领域的各个学科,诸如航天工程,航空工程,机械工程,能源工程,动力工程,交

2、通工程,土木工程,工程力学等等。作为固体力学的一门主要分支学科,结构动力学起源于经典牛顿力学,就是牛顿质点力学。质点力学的基本问题是用牛顿第二定律来建立公式的。牛顿质点力学,拉格朗日力学和哈密尔顿力学是结构动力学基本理论体系组成的三大支柱。经典动力学的理论体系早在19世纪中叶就已建立,。但和弹性力学类似,理论体系虽早已建立,但由于数学求解上的异常困难,能够用来解析求解的实际问题实在是少之又少,能够通过手算完成的也不过仅仅限于几个自由度的结构动力体系。因此,在很长一段时间内,动力学的求解思想在工程实际中并未得到很好的应用,人们依然习惯于在静力学的范畴内用静力学的方法来解决工程实际问题。随着汽车,

3、飞机等新时代交通工具的出现,后工业革命时代各种大型机械的创造发明,以及越来越多的摩天大楼的拔地而起,工程界日新月异的发展和变化对工程师们提出了越来越高的要求,传统的只考虑静力荷载的设计理念和设计方法显然已经跟不上时代的要求了。也正是从这个时候起,结构动力学作为一门学科,也开始受到工程界越来越高的重视,从而带动了结构动力学的快速发展。结构动力学这门学科在过去几十年来所经历的深刻变革,其主要原因也正是由于电子计算机的问世使得大型结构动力体系数值解的得到成为可能。由于电子计算机的超快速度的计算能力,使得在过去凭借手工根本无法求解的问题得到了解决。目前,由于广泛地应用了快速傅立叶变换FFT,促使结构动

4、力学分析发生了更加深刻地变化,而且使得结构动力学分析与结构动力试验之间的相互关系也开始得以沟通。总之,计算机革命带来了结构动力学求解方法的本质改变。作为一门课程,结构动力学的基本体系和内容主要包括以下几个部分：单自由度系统结构动力学,；多自由度系统结构动力学,；连续系统结构动力学。此外,如果系统上所施加的动力荷载是确定性的,该系统就称为确定性结构动力系统；而如果系统上所施加的动力荷载是非确定性的,该系统就称为概率性结构动力系统。1.2主要理论分析结构的质量是一连续的空间函数,因此结构的运动方程是一个含有空间坐标和时间的偏微分方程,只是对某些简单结构,这些方程才有可能直接求解。对于绝大多数实际结

5、构,在工程分析中主要采用数值方法。作法是先把结构离散化成为一个具有有限自由度的数学模型,在确定载荷后,导出模型的运动方程,然后选用合适的方法求解。1.3载荷确定 载荷有三个因素,即大小、方向和作用点。如果这些因素随时间缓慢变化,则在求解结构的响应时,可把载荷作为静载荷处理以简化计算。载荷的变化或结构的振动是否 缓慢, 只是一个相对的概念。如果载荷的变化周期在结构自由振动周期的五、六倍以上,把它当作静载荷将不会带来多少误差。若载荷的变化周期接近于结构的自由振动周期,即使载荷很小,结构也会因共振见线性振动而产生很大的响应,因而必须用结构动力学的方法加以分析。 动载荷按其随时间的变化规律可以分为：周

6、期性载荷,其特点是在多次循环中载荷相继呈现相同的时间历程,如旋转机械装置因质量不平衡而引起的离心力。周期性载荷可借助傅里叶分析分解成一系列简谐分量之和。 冲击载荷, 其特点是载荷的大小在极短的时间内有较大的变化。冲击波或爆炸是冲击载荷的典型来源。随机载荷,其时间历程不能用确定的时间函数而只能用统计信息描述。由大气湍流引起的作用在飞行器上的气动载荷和由地震波引起的作用在结构物上的载荷均属此类。对于随机载荷,需要根据大量的统计资料制定出相应的载荷时间历程载荷谱。对于前两种载荷,可以从运动方程解出位移的时间历程并进一步求出应力的时间历程。对于随机载荷,只能求出位移响应的统计信息而不能得到确定的时间历

7、程,因而须作专门分析才能求出应力响应的统计信息。 1.4体系的动力自由度为了确定一个体系在振动过程中全部质量的位置所需独立几何参数的数目,称为动力自由度或简称自由度。这些参数通常表示质量的线位移或转角,它们也就是动力计算中的基本未知量。实际结构的质量是连续分布的,是无限自由度体系。为了简化计算,常按下面的方法进行简化。1集中质量法从物理的角度提供一种减少动力自由度的简化方法。把连续分布的质量根据静力等效原则集中为几个质点。这样就把无限自由度体系简化成有限自由度体系。具体分为：不计轴向变形的均质简支梁；三层平面刚架在水平力作用下计算侧向振动和块形基础。2广义位移法具有分布质量的简支梁的振动曲线位

8、移曲线,可近似地用三角级数表示为式中,是一组给定的函数,称作位移函数或形状函数,与时间无关。是一组待定参数,称作广义坐标,随时间而变化。因此,体系在任一时刻的位置是由广义坐标来确定的。注意：这里的形状函数应满足位移边界条件,所选的函数形式可以是任意的连续函数。因此,上式可写成更一般的形式式中,是自动满足位移边界条件的函数集合中任意选取的n个函数。广义坐标法将应用于后面的振型叠加法和能量法。3有限单元法可看作广义坐标法的一种特殊应用。把体系的离散化和单元的广义坐标二者结合起来,就构成了有限单元的概念。其具体作法是：第一,将结构离散为有限个单元本例为3个单元；第二,取结点的位移参数挠度y和转角作为

9、广义坐标,本例为、和、。第三,分别给出与结点位移参数均为1时相应的形状函数称作插值函数它们确定了指定结点位移之间的形状；第四,仿照广义位移法的公式,体系的位移曲线可用4个广义坐标及其形状函数表示为：可事先给定,让其满足边界条件,这样就把无限自由度体系简化为4个自由度体系和。有限元法综合集中质量法和广义坐标法的优点：与广义坐标法相似,有限元法采用了形函数的概念,但不同于广义坐标法在全部体系上插值,而是采用了分片的插值,因此形函数的公式可以相对简单。与集中质量法相比,有限元法中的广义坐标也采用了真实的物理量,具有直接、直观的优点,这与集中质量法相同。1.5体系振动时能量的耗散与阻尼力实际结构在自由

10、振动时有衰减现象,振幅随时间逐渐减小,最后趋于静止；在强迫振动时,外荷载需对结构不断做功,才能维持振幅不变稳态振动。这都表明在振动工程中会产生能量的耗散,这种消耗能量并使振动衰减的因素,成为阻尼。在动力计算时,要先建立结构的振动方程,为了能反映振动过程中的能量耗散,在建立方程时须引入一个造成能量耗散的阻尼力。而这个力的引入提高了运动方程计算的难度。在结构动力分析时,由于粘滞阻尼力的分析比较简单,其他类型的阻尼力也可以简化为等效粘滞阻尼力来分析。因此,本书只讨论粘滞阻尼力的情形。1.6建立振动方程的方法动力问题主要是求出位移或位移参数随时间变化的反应。建立振动方程的常用方法有四种,分别介绍如下。

11、1动力平衡法此法也称达朗伯原理的直接平衡法。根据牛顿第二运动定律,任何质量m的动量的变化率等于作用在这个质量上的力式中y为动位移。若m不随时间变化,上式可写成上式中第一项为作用在质量上的力,第二项可以称为质量m的惯性力。质量所产生的惯性力,与它的加速度成正比,但方向相反。这个概念称作达朗伯原理。有第二式可以看出,在引入达朗伯原理后,与静力学中的平衡方程的表达式相识,及作用于质量上的所有里保持平衡,常称此法为动静法。本方法的优点在于物理概念清楚,形象鲜明。缺点是解决复杂问题时困难较大,且不便用它来推证某些结论。2虚功法当结构比较复杂,如所包含的各种力可以容易的用位移自由度来表示,而它们的平衡规律

12、可能不清楚或很复杂。此时,运用给予虚位移原理的虚功法来建立运动方程就较方便。按照虚位移原理,虚位移时所作的总虚功为0是与平衡条件等价的。在建立体系的方程时,先确定作用于质量上的所有力,包括惯性力；然后引入相应于每个自由度的虚位移,并使所做的总虚功等于0,从而得出振动方程。此方法的优点是适应性强,可用它推出运动的普遍规律；虚功是标量可以按照代数规则计算避免复杂的矢量计算。缺点是比较抽象。3变分法用基于哈密顿原理以变分形式表示的能量关系来建立动力平衡方程。哈密顿原理可以表达为式中,T为体系的总动能,V为体系的势能,包括应变能及任何保守外力的势能,Wnc为作用于体系上的非保守力所做的功,为在指定时间

13、区间内所取得变分。哈密顿原理表明,在任何时间区间内,动能和势能的变分加上所考虑的费保守力所做的功的变分必须等于0。应用这个原理可以直接导出任何体系的振动方程。这个方法和虚功法的区别是：此方法中,不明显使用惯性力和弹性力,而是用动能和势能的变分项来代替。4能量法基于能量守恒原理的能量法,不仅可以用来建立体系的振动方程,而且可以用来直接计算体系的自振频率。2 单自由度系统的振动振动系统可分为离散模型和连续模型两种不同的类型。离散模型具有有限个自由度,而连续模型则具有无限个自由度。系统的自由度定义为能完全描述系统运动所必须的独立的坐标个数。在离散模型中,最简单的是单自由度线性系统,它用一个二阶常系数

14、微分方程来描述。这类模型常用来作为较复杂系统的初步近似描述。另外在后续章节将会讲到,复杂系统的数学模型可通过模态分析技术转化为一组独立的二阶常微分方程,其中每一个方程都类似于单自由度系统的运动方程。因此,对单自由度系统进行详细深入的分析是十分必要的。单自由度振动系统,指用一个独立参量便可确定系统位置的振动系统。只要以它的平衡位置取为坐标原点,任一瞬时的质点坐标x线位移或q角位移就可以决定振动质点的瞬时位置。 2.1单自由度系统无阻尼自由振动 根据牛顿定律：令2=k/m,求微分方程的解,得将其合成一个简谐振动,并代入初始条件：t=0时,x=x0,2.1.1固有频率系统的圆频率和频率只与系统本身的

15、物理性质弹性和惯性有关,因此当振动系统的结构确定后,系统的振动频率就固定不变,而不管运动的初始条件如何,也和振幅的大小无关,因此成为固有圆频率和固有频率。2.1.2固有频率计算方法1公式法。根据公式计算2静变形法。根据质量块所处平衡位置的弹簧变形计算。3能量法。根据能量守恒定律,由于无阻尼,无能量损失,将的方程代入上式,系统的最大动能等于系统的最大弹性势能,计算求出。4瑞利法。考虑到系统弹簧质量的计算方法,如假设系统的静态变形曲线作为假定的振动形式,根据推倒,得出系统的固有频率为 ,式中加入的部分为弹簧等效质量不同振动系统的等效质量不同,只需先算出弹性元件的动能,根据,计算即可。2.2单自由度

16、有阻尼系统始终与物体振动方向相反,对系统做负功,不断消耗系统能量的阻力,统称为阻尼,单位为N。2.2.1分类：1干摩擦阻尼两个干燥表面相互压紧并做相对运动,阻力的大小取决于摩擦副的材料和法相压力：F=mN；2粘滞阻尼物体以中等速度在流体中运动时所产生的阻尼有润滑油的润画面之间,与速度的一次方成正比,即：,其中阻尼系数取决于物体的形状、尺寸及润滑剂介质的粘性,单位为Ns/cm。当以较大速度运动时3m/s以上,阻力与速度平方成正比3结构阻尼材料在变形过程中由内部晶体之间的摩擦所产生的阻尼,大小取决于材料性质。2.2.2固有频率1强阻尼状态,an系统不再是振动,是逐渐回复到平衡位置的非周期运动2临界

17、阻尼状态,a=n也不是振动,是逐渐回到平衡位置的非周期运动。3弱阻尼状态,an,是振幅被限制在之内的,随时间不短衰减的衰减振动。只有在这种状态下才存在频率和周期。2.2.3共振当激振力F周期力频率与固有频率相同时,振幅最大,称为共振。1产生共振的原因以简谐激振力为例,设F=P0sint,q=P0/m得解方程得由于第一项很快将衰减掉,故只研究第二项,并将其求导代入原式,经过化简整理,将两式子联立,令B0=P0/k称为静变位；l=/n称为频率比；z=a/n称为阻尼比,代入以B/B0为纵坐标,以/n为横坐标,以阻尼比z为参变量,做幅频响应曲线。l1时,振幅的大小主要取决于系统惯性。2阻尼对共振振幅的

18、影响阻尼仅在共振区对振幅有较大影响,阻尼增大,共振振幅减小。但在共振区域外对振幅影响较小。2.3产生振动的原因1自由振动。由于系统平衡受到破坏,靠弹性恢复力来维持的振动2受迫振动。在外界激振力的持续作用下,或位移干扰或偏心质量的作用下产生的3自激振动。由于系统的非震荡性能源和反馈特性,从而引起的一种稳定的周期性振动。3多自由度系统的振动3.1多自由度系统的数学模型实际上的动力学系统都是连续系统,为了能够在一个连续系统中生成一个N个自由度的多自由度模型,可以假定该连续系统的位移能够近似表示为：其中为假定的N个振型函数,即相当于把这一连续系统看作包含N个独立广义坐标的多自由度系统。然后以代入动能,

19、位能和非保守力虚功的表达式,应用拉格朗日方程推导得到多自由度系统的动力学基本方程为：这里应注意所选用的振型函数必须为一组线性独立的函数,每个所具有的导数必须等于位能中所出现的阶数,并且必须满足所有的位移边界条件。3.2多自由度系统的自由振动求解多自由度系统自由振动的动力特性,首先仍然是要求解出系统的无阻尼固有圆频率。这是一个N阶代数特征值问题,需要求解特征方程。其结果为N次的多项式公式,根为特征值,即固有频率的平方。这些频率可以从低到高排列：相应于每一个特征值将有一个特征向量,即固有模态：模态仅是确定内部的一个定数乘子,可用任何方便的方式度量。3.3多自由度系统的动力响应对线性多自由度系统,且

20、为非耦连阻尼效应时,可采用振型叠加法求解其动力响应。振型叠加法包括振型位移法和振型加速度法,核心思想是：首先解得自由振动下的频率和模态,然后将模态集成为模态矩阵,利用模态矩阵进行坐标变换,将物理坐标变换为主坐标,从而实现动力学方程的解耦。在求得系统的固有频率和模态后,将其正规化,根据模态来集成模态矩阵在振型叠加的过程中的关键步骤是要导出坐标变换：坐标称为主坐标。将上式方程代入多自由度系统的动力学基本方程,同时将方程结果乘以,得到主坐标的运动方程,即：其中： 模态质量矩阵 模态阻尼矩阵 模态刚度矩阵 模态力向量 4自振频率和振型的计算4.1能量法求自振频率瑞利能量法设有一n自由度系统,其质量矩阵

21、为M,刚度矩阵为K,系统的主振动x第一瑞利商当X取系统各阶主振型Xi时,瑞利商RI就相应地给出各阶固有频率pi的平方式中,X代表假设的主振型,那末由瑞利商只能得出相应固有频率的估值。理论上,它适用于求各阶固有频率。实际上,因为系统高阶主振型很难作出合理假设,所以该式往往只有在估算系统基频p1时才是切实可行的。从柔度法可得到另一形式的系统固有频率的估式。第二瑞利商式中的X取为系统的主振型时,RII也就给出相应的特征值 i ,i=1,2,n两种形式瑞利能量法各有所长。前者适用于刚度矩阵已知的情形,后者适用于柔度矩阵已知的情形。一般说来,前者较为简便,后者较为准确。4.1.2李兹能量法里兹法是瑞利法

22、的改进,用里兹法不仅可以计算系统的基频,还可以算出系统的前几阶频率和模态,瑞利法算出的基频的精度取决于假设的振型对第一阶主振型的近似程度,而且得到的基频总是精确值的上限,里兹法将对近似振型给出更合理的假设,从而使算出的基频值进一步下降。里兹法基于与瑞利法相同的原理,但将瑞利使用的单个假设模态改进为若干个独立的假设模态的线性组合：,代入瑞利商：,由于在系统中的真实主振型处取驻值,所以 A 的各个元素应当从下式确定：。代入：其中是 r 阶单位矩阵的第 j 列。上面 r 个方程可合成为：表示将函数分别对 A 的各个元素依次求偏导,然后排列成列向量。同理,有： 两项代入得：由于的阶数 r 一般远小于系

23、统自由度数 n,上式所示的矩阵特征值问题比原来系统的矩阵特征值问题解起来容易得多,因此里兹法实际上是一种缩减系统自由度求解固有振动的近似方法,就是自由度缩减为 r 的新系统的刚度矩阵和质量矩阵。可求出 r 个特征根,及相应的特征向量,原来系统的前 r 阶固有频率可近似取为：,相应的前 r 阶主振型近似取为：。正交性分析即：因此：得出的近似主振型式关于矩阵 M 和 K 相互正交。4.2瑞利李兹法虽然应用瑞利法可以较为精确地求解得到系统的第一振型,但在结构动力分析中,为了得到较为精确的结果,常常需要一个以上的振型。李兹把瑞利法加以推广,是计算前几个振型的最为方便的方法之一。其基本假设是用一组假设的

24、形状矢量和一组幅值来表示位移矢量,其实质是缩减自由度方法的一种。通过坐标的变换公式,把具有N个自由度的体系转化为用2s个广义坐标和相应的假定李兹基表示的2s个自由度体系,这样就使得求解变得相对容易些。用瑞利李兹法求得的近似频率,对最低几个振型一般精度较高,而对较高阶的振型精度就相对差得多。其解的精度与基矢量的选取有关。4.3子空间迭代法前述几种求解结构体系自由振动频率和振型的解法大多适用于自由度数目较小的小型动力体系。而对于实际的结构,其自由度的数目往往会达到几千个,甚至有时会达到几十万个。针对如此超大型的结构体系,不可能也没有必要求解其所有的振型和频率,而往往仅需要知道其最低的若干阶振型。因

25、此,为适应这种求解需要,产生了基于多种计算技术基础之上的子空间迭代法。子空间迭代法可以求解结构体系的前面最低的p阶频率和振型,是用来解决大型结构振动问题的发展较早的富有成效的方法之一。其基本点在于假设r个起始向量同时进行迭代以求得矩阵的前ss个特征值和特征向量,一般可取r2s。它的一个很大的优点在于,可以按照任意要求的精度求得振型坐标,而其它的坐标缩减方法内包含的近似性,使得结果的精度无法估计。5 结构动力学在抗震设计中的应用5.1序言地震时地面运动是一个复杂的时间空间过程。结构地震响应应取决于地震动特性和结构特性,特别是结构的动力特性。结构地震响应分析的水平也是随着人们对这两方面认识的逐步深

26、入而提高的。近几十年来,人们对地震动的谱成分和各类结构的动力特征有了深入认识。因此,结构的分析也随之有了相应的进展。结构地震反应分析的发展经过了静力法、反应谱法、动力法三个阶段。反应谱法根据单自由度系统的地震响应,既考虑了结构动力特性与地震动特性之间的动力关系,又保持了静力法的形式,在各国结构抗震设计规范中已被广泛采用。 现行的抗震设计方法包括反应谱法和时程分析法。5.2方法比较根据建筑结构抗震规范,对单自由度体系,给定场地条件以及结构的自振周期和阻尼比,便可以从反应谱中获得结构的最大地震响应位移、速度和加速度,进而可求出结构的地震力。对于多自由度体系,首先采用多自由度体系的反应谱理论,即先利

27、用模态分析法将多自由度体系分解为一系列广义单自由度体系,最后将各振型的最大值用一定的振型组合方法组合出结构的最大地震反应。由于反应谱方法基本正确地反映了地震动特性,并考虑了结构的动力特性,所以对于一般的结构而言,具有良好的精度,且概念明确,计算方便。静力法假设结构各部分水平加速度与地面运动水平加速度完全一样。 地震地面运动是一个非平稳随机过程,而随机振动法充分考虑了地震发生的概率特性,所以普遍认为随机振动法是一种合理的分析方法。但是,随机振动法的缺点是它的计算量庞大而且对于非线性问题可能引起较大的误差,在处理罕遇地震下的强非线性问题时有其局限性。时程分析法是确定性动力分析方法的一种,是发展较为

28、成熟、应用较多的一种方法。由于这种分析方法是在离散时间点上一步一步地求响应的数值解,所以该法可以在任一时间点上随时修改结构参数,很适合于处理参数随时间变化的非线性问题。它既可虑地震波的多维多点输入,还可以考虑结构几何非线性、物理非线性、非比例阻尼和桩土结构相互作用等的地震反应。常用的积分方法有线性加速度法。这里主要介绍比较先进的时程分析法：逐步积分数值方法特别适用于计算大型结构在地震作用下的动力响应,其无需像振型叠加法那样要预先花费很多的工作量计算频率和振型。此外,由于计算中考虑几何非线性大变形的影响,本文中采用Newmark逐步积分方法求解。注意： 1在进行时程分析过程中,利用上述方法计算结

29、构反应关键的是地震动的描述,即恰当地输入地震波。 2分析和结果存在一定的局限性,即计算结果仅仅是选择地震波的反应,若选择另外一条地震波,计算结果差别很大。 3为得到结构反应的统计结果,必须对多条地震波进行分析,工作量较大。6总结通过本课程的学习,我了解到：结构的动力计算与静力计算有很大的区别。静力计算是研究静荷载作用下的平衡问题。这时结构的质量不随时间快速移动,因而无惯性力。动力计算研究的是动荷载作用下的运动问题,这时结构的质量随时间快速运动,惯性力的作用成为必须考虑的重要问题。根据达朗贝尔原理,动力计算可以转化为静力平衡问题来处理。但是,这是一种形式上的平衡。也就是说,动力计算中,虽然形式仍是是在列平衡方程,但是这里要注意两个问题：所考虑的力系中要包括惯性力这个新力,考虑的是瞬间的平衡,荷载、内力等都是时间的函数。我们首先学习了单自由度系统自由振动和受迫振动的概念,所以在学习多自由度系统和弹性体的振动分析时,则重点学习后者的振动特点以及前者的联系和区别,这样既节省了时间,又抓住了重点。参考文献1建筑抗震设计规范GB 50011-2010北京：出版社,20102结构动力学 张子明 周星德姜冬菊编著 中国电力出版社3高层建筑混凝土结构技术规程JGJ3-2002.北京：出版社,2001.