机械振动学最新版本

1、 振动概念（vibration）物体经过它的静平衡位置所做的往复运动。或者说某一物理量在其平衡位置或平衡值附近来回的变动。 振动首先是一种运动。比如：地壳的运动、交流电、电磁波、潮水的涨落等。2 机械振动的研究对象和分类机械振动的研究对象和分类2.1 2.1 研究对象研究对象“振动系统振动系统”第一章第一章 绪绪 论论 系统的定义：系统的定义： 由若干个元素构成的有机组合，个元素间存在着相互作用、互相影响的关系。 机械系统的定义：机械系统的定义： 由若干个机械元件组成的系统。具体的讲，是由运动副连接的一些构件所组成的能完成一定运动的机械装置。第一章第一章 绪绪 论论2.2 2.2 机械系统研究

2、内容机械系统研究内容 系统（系统（S） 输入（X） 输出（Y）激励响应响应第一章第一章 绪绪 论论系统的研究内容包括三个方面：系统的研究内容包括三个方面：已知系统的输入（X）和系统（S）,求输出(Y)系统的动力响应分析，或叫动态分析。已知系统的输入（X）和输出(Y)，求系统（S）系统设计；系统识别或系统辨识。已知系统的系统（S）和输出(Y)，求输入（X）环境预测。 自由振动：给图中质量块一个激励，给一个初始位移后，质量块就开始振下去。 强迫振动：用一个电机作元件，给系统一个持续激励，系统会在电机的强制激励下振动。 自激振动：扬声器的鸣叫声。3 3 机械振动的分类机械振动的分类3.1 3.1 按

3、输入分按输入分mk第一章第一章 绪绪 论论 简谐振动：符合正弦（预选）规律的振动。 周期振动：x（t）x（t+kT）， 瞬态振动：风铃随风而动；地震 随机振动：不能用当前的现象预测未来，但是符合统计学规律，可以用统计的方法来研究。如，烟的运动；红旗的飘动。3.2 3.2 按输出分按输出分第一章第一章 绪绪 论论 自由度：用来描述一个物体确定运动的独立坐标。 单自由度系统： 多自由度系统： 可以是两个、三个甚至是n个自由度系统，n个独立坐标，n维空间。 连续系统：用偏微分方程描述3.3 3.3 按自由度划分按自由度划分 ),.,(vdxtHH可用微分方程描述第一章第一章 绪绪 论论 线性振动 非

4、线性振动:二阶常系数线性齐次）(0kxxm 3.4 3.4 按微分方程分按微分方程分单摆振动方程）(0sinxkxm 第一章第一章 绪绪 论论4 4 主要参考文献主要参考文献 书书+期刊期刊书：张策、张维平、邵韧平、闻邦春、书：张策、张维平、邵韧平、闻邦春、 李有堂、张义民等李有堂、张义民等期刊：期刊：噪声与振动噪声与振动 （sound and vibration）第一章第一章 绪绪 论论2.1 一些基本概念、无阻尼单自由度振动系统一些基本概念、无阻尼单自由度振动系统 2.3 有线性阻尼有线性阻尼自由振动自由振动2.4 简谐激励力作用下的强迫振动简谐激励力作用下的强迫振动 2.8 隔振原理隔振

5、原理2.5 周期激励下的响应周期激励下的响应2.6 任意激励下的响应任意激励下的响应2.7 简谐力的功和等效阻尼简谐力的功和等效阻尼2.2 固有频率的计算固有频率的计算 当物体沿当物体沿x x轴作直线运动时，惯性的大小可用质量来表轴作直线运动时，惯性的大小可用质量来表示。根据牛顿第二定律，作用在物体上的外力示。根据牛顿第二定律，作用在物体上的外力F F，物体由此，物体由此产生的加速度和物体质量产生的加速度和物体质量m m之间有下述关系：之间有下述关系：) 1-(122dtxdmF 构成机械振动系统的基本元素有惯性、恢复性和阻尼。构成机械振动系统的基本元素有惯性、恢复性和阻尼。惯性就是能使物体当

6、前运动持续下去的性质。恢复性就是能惯性就是能使物体当前运动持续下去的性质。恢复性就是能使物体位置恢复到平衡状态的性质。阻尼就是阻碍物体运动使物体位置恢复到平衡状态的性质。阻尼就是阻碍物体运动的性质。从能量的角度看，惯性是保持动能的元素，恢复性的性质。从能量的角度看，惯性是保持动能的元素，恢复性是贮存势能的元素，阻尼是使能量散逸的元素。是贮存势能的元素，阻尼是使能量散逸的元素。构成机械振动系统的基本元素构成机械振动系统的基本元素质量的单位为质量的单位为kgkg。阻尼力阻尼力Fd反映阻尼的强弱，通常是速度反映阻尼的强弱，通常是速度x的函数，阻尼力的函数，阻尼力可表示为可表示为 这种阻尼称为粘性阻尼

7、。比例常数这种阻尼称为粘性阻尼。比例常数c称为粘性阻尼系称为粘性阻尼系数，单位数，单位N.s/m。 )31 ( xcFd 典型恢复性元件是弹簧，弹簧产生的恢复力是该元件位移典型恢复性元件是弹簧，弹簧产生的恢复力是该元件位移的函数，即的函数，即Fs=Fs（x）。当。当Fs（x）是线性函数时，有：是线性函数时，有： Fs=kx (1-2) k称为弹簧常数或弹簧的刚度系数。单位为称为弹簧常数或弹簧的刚度系数。单位为N/m。质量、弹簧和阻尼器质量、弹簧和阻尼器是构成机械振动系统是构成机械振动系统物理模型物理模型的的三个基本元件。三个基本元件。 自由度与广义坐标自由度与广义坐标 自由度数自由度数: 完全

8、确定系统运动所需的独立坐标数目称为自由度数。完全确定系统运动所需的独立坐标数目称为自由度数。 刚体在空间有刚体在空间有6个自由度：三个方向的移动和绕三个方向的转动，个自由度：三个方向的移动和绕三个方向的转动，如飞机、轮船；如飞机、轮船； 质点在空间有质点在空间有3个自由度：三个方向的移动，如高尔夫球；个自由度：三个方向的移动，如高尔夫球； 质点在平面有质点在平面有2个自由度：两个方向的移动，加上约束则成为单个自由度：两个方向的移动，加上约束则成为单自由度。自由度。质量元件质量元件 无弹性、不耗能的刚体，储存动能的元件无弹性、不耗能的刚体，储存动能的元件 xmFm 平动：平动：力、质量和加速度的

9、单位分别力、质量和加速度的单位分别为为N、kg和和m / s 2。 JTm转动：转动：力矩、转动惯量和角加速度的力矩、转动惯量和角加速度的单位分别为单位分别为Nm、kg m 2和和rad / s 2 2.1 离散系统的组成离散系统的组成2.1 离散系统的组成离散系统的组成弹性元件弹性元件 无质量、不耗能，储存势能的元件无质量、不耗能，储存势能的元件 xkFs平动：平动：力、刚度和位移的单位分别为力、刚度和位移的单位分别为N、N / m和和m 。tskT 转动：转动：力矩、扭转刚度和角位移的单力矩、扭转刚度和角位移的单位分别为位分别为Nm、 Nm / rad和和rad 阻尼元件阻尼元件 无质量、

10、无弹性、线性耗能元件无质量、无弹性、线性耗能元件 xcFd平动：平动：力、阻尼系数和速度的单位分力、阻尼系数和速度的单位分别为别为N、N s/ m和和m/s。tdcT 转动：转动：力矩、扭转阻尼系数和角速度力矩、扭转阻尼系数和角速度的单位分别为的单位分别为Nm、 Nms / rad和和rad/s 2.1 离散系统的组成离散系统的组成等效弹簧刚度等效弹簧刚度 斜向布置的弹簧斜向布置的弹簧 2ecos/kxFkxx串联弹簧串联弹簧 并联弹簧并联弹簧 niikk1eniikk1e11niicc1eniicc1e11并联系统并联系统串联系统串联系统等效阻尼系数等效阻尼系数 传动系统的等效刚度传动系统的

11、等效刚度 21 te1 t/ikk传动系统的等效阻尼传动系统的等效阻尼 ct1e= ct1 / i 221e1/iJJ等效质量等效质量 传动系统的等效惯量传动系统的等效惯量 单自由度系统的类型tQkxxrxmtQkxxmkxxrxmkxxmsinsin0000 单自由度无阻尼自由振动单自由度无阻尼自由振动单自由度有粘性阻尼的自由振动单自由度有粘性阻尼的自由振动单自由度无阻尼受迫振动单自由度无阻尼受迫振动单自由度有粘性阻尼的受迫振动单自由度有粘性阻尼的受迫振动机机 械械 振振 动动 学学例：如右图，舍振动体的例：如右图，舍振动体的质量为质量为m m，它所受的重，它所受的重力为力为W W，弹簧刚度

12、为，弹簧刚度为k,k,弹簧挂上质量块的静伸弹簧挂上质量块的静伸成量为成量为j j，此时系统，此时系统处于静平衡状态，平衡处于静平衡状态，平衡位置为位置为0-00-0，求给系统，求给系统一个初始扰动后系统的一个初始扰动后系统的振动方程。振动方程。模型的建立模型的建立机机 械械 振振 动动 学学单自由度无阻尼自由振动系统单自由度无阻尼自由振动系统无阻尼自由振动：无阻尼自由振动： 振动系统受到初始扰动后，不再受到外振动系统受到初始扰动后，不再受到外力作用，也不受阻尼的影响所作的振动。力作用，也不受阻尼的影响所作的振动。静平衡振动系统产生弹性恢复力弹力重力静平衡破坏初始扰动机机 械械 振振 动动 学学

13、单自由度无阻尼自由振动系统单自由度无阻尼自由振动系统解：解：取静平衡位置为坐取静平衡位置为坐标原点，以标原点，以X轴为系统轴为系统的坐标轴，向下为正的坐标轴，向下为正方向建立坐标系。方向建立坐标系。 以以x x表示质量块的受扰表示质量块的受扰后的位移，当质量块后的位移，当质量块离开平衡位置时，在离开平衡位置时，在质量块上作用的力有质量块上作用的力有：XTWmgkxkTjW重力弹性恢复力x 由于受力不平衡，质量块产生加速度由于受力不平衡，质量块产生加速度机机 械械 振振 动动 学学单自由度无阻尼自由振动系统单自由度无阻尼自由振动系统根据牛顿第二定律建立振动微分方程：xmxkwj )(0,022x

14、xmkkxxmnn ：则上式可写成令即叫做系统的固有频率2n二阶齐次常系数微分方程，stex 机机 械械 振振 动动 学学单自由度无阻尼自由振动系统单自由度无阻尼自由振动系统扭转振动问题扭转振动问题例1-2： 右图所示，垂直轴的下端右图所示，垂直轴的下端固定一个水平圆盘。已知固定一个水平圆盘。已知轴长为轴长为l ,l ,直径为直径为d,d,剪切剪切弹性模量为弹性模量为G,G,圆盘的转动圆盘的转动惯量为惯量为I I，在盘上施加初，在盘上施加初始扰动后（如力偶），系始扰动后（如力偶），系统做自由扭转振动。若不统做自由扭转振动。若不计阻尼影响，振动将永远计阻尼影响，振动将永远持续下去。求系统的振动持

15、续下去。求系统的振动方程。方程。机机 械械 振振 动动 学学单自由度无阻尼自由振动系统单自由度无阻尼自由振动系统由材料力学知：扭转刚度为：324Gdk)/(J0JIsradkkkn系统固有角频率令即扭转振动微分方程为：建立如图所示坐标系， 0)1(212nnsHzIkf 代入微分方程得到：将系统振动的固有频率：机机 械械 振振 动动 学学单自由度无阻尼自由振动系统单自由度无阻尼自由振动系统典型的单自由度自由振动单摆例1-3：如左图所示，求t时刻刚体的角度是多少？机机 械械 振振 动动 学学单自由度无阻尼自由振动系统单自由度无阻尼自由振动系统解：以静平衡位置为原点，以以静平衡位置为原点，以角增角

16、增加的方向为正方向建立坐标系。加的方向为正方向建立坐标系。隔离物体，进行受力分析。隔离物体，进行受力分析。 使用牛顿定律建立振动模型：使用牛顿定律建立振动模型：力矩形式：力矩形式：0sin0sinsinmglJmglJJmgl 作为摆动时，即力形式：力形式：？机机 械械 振振 动动 学学单自由度无阻尼自由振动系统单自由度无阻尼自由振动系统txtxtxxbxbbbxbbxtnnnnnnnnnncossin)(0sin0cos0cos0sin0000201210210时有：代入初始条件：1-2 无阻尼单自由度系统的自由振动规律为：高等数学知方程的通解度和初位移均为零。此初始条件亦即：初速预先给定初

17、始条件：考虑00002,0 xxxxxxttn )sin(cossin)(21tAtbtbtxnnn.;A2112221为频率三要素：nbbtgbb机机 械械 振振 动动 学学单自由度无阻尼自由振动系统单自由度无阻尼自由振动系统结结 论论单自由度无阻尼自由振动系统的方程是一样的单自由度无阻尼自由振动系统的方程是一样的,规律是相同的，具有以下特点：规律是相同的，具有以下特点： 1.1.单自由度无阻尼振动是简谐的。单自由度无阻尼振动是简谐的。 2.2.振幅决定于初始条件：振幅决定于初始条件： )(;222120220bbAxxAn图中系统，用手把图中系统，用手把m m移到移到X X0 0位置，初始

18、位移的大小决位置，初始位移的大小决定于定于m m的振幅，如果放手的同时，给的振幅，如果放手的同时，给m m一个右向的初一个右向的初速度，可以通过上式计算出其最大振幅。速度，可以通过上式计算出其最大振幅。机机 械械 振振 动动 学学单自由度无阻尼自由振动系统单自由度无阻尼自由振动系统固有频率与初始条件无关。系统一定，固有频固有频率与初始条件无关。系统一定，固有频率一定。率一定。 fTfmknn1;2;的特点，座钟。应用：利用“等时性”思考：钟表的钟摆的摆角大是准确还是小准确？思考：钟表的钟摆的摆角大是准确还是小准确？结结 论论机机 械械 振振 动动 学学单自由度无阻尼自由振动系统单自由度无阻尼自

19、由振动系统 在振动研究中，计算振动系统的固有频率有很重要的意义 ，除用定义法（牛顿法）外，通常还有以下几种常用的方法，即静变形法、能量法和瑞利法，现分别加以介绍。2.2 计算系统固有频率的其它方法计算系统固有频率的其它方法1、静变形法（、静变形法（Static Deformation Method） Wkj 当单振子处于静平衡状态时，弹簧的弹性力与振动质量的重力互相平衡，即存在关系式： 由上式可得：jjmgWk故系统的固有频率为：) 1 (2121jngmkf由此可见，只要知道质量块处的弹性静变形，就可以计算出系统的固有频率。在有些实际问题中，不能直接给出系统的弹簧刚度时，利用此法计算固有频率

20、比较方便。例例1 1 设一悬臂梁长度为，抗弯刚度为，自由端有一集中质量。设一悬臂梁长度为，抗弯刚度为，自由端有一集中质量。梁本身重量忽略不计。试求这一系统的固有频率（见下图）。梁本身重量忽略不计。试求这一系统的固有频率（见下图）。自由端有集中质量的悬臂梁 解：悬臂梁在自由端由集中力mg所引起的静挠度为： EJmglj33) 1 (3213mlEJfn当不易用计算方法求出静挠度时，也可用实测方法得到静挠度，然后按（1）式计算系统固有频率。2.2 计算系统固有频率的其它方法计算系统固有频率的其它方法2 2、能量法（、能量法（Energy MethodEnergy Method）在无阻尼自由振动系统

21、中，由于没有能量的损失，所以振幅在无阻尼自由振动系统中，由于没有能量的损失，所以振幅始终保持为一常数，即在振动过程中振幅始终不衰减。我们始终保持为一常数，即在振动过程中振幅始终不衰减。我们将这样的系统称为将这样的系统称为保守系统保守系统。在保守系统中，根据机械能守恒定律，在整个振动过程的在保守系统中，根据机械能守恒定律，在整个振动过程的任一瞬时机械能应保持不变。任一瞬时机械能应保持不变。式中：式中：T T系统中运动质量所具有的动能；系统中运动质量所具有的动能；U U系统由于弹性变形而储存的弹性势能，或由于重力作功系统由于弹性变形而储存的弹性势能，或由于重力作功 而产生的重力势能。而产生的重力势

22、能。0UTdtd即： T+U=常数 或2.2 计算系统固有频率的其它方法计算系统固有频率的其它方法xkxmgxkxdxmgx0221222121kxkxmgxmgxU221xmT 对于单自由度无阻尼自由振动系统来说，系统的动能为：1. 重力势能：当质量块m低于静平衡位置时，重力势能为-mgx。2. 弹性势能：当质量块m运动至离静平衡位置距离+x时，弹簧的弹性力对质量块所作的功即为系统此时的弹性势能。如下图所示，系统的弹性势能为：故系统的势能为故系统的势能为：)2()(212122常数Ekxxm 所以：所以：系统的势能则由以下两部分组成： 22xmmgx单自由度振动系统的弹性势能这就是单自由度无

23、阻尼自由振动系统的能量方程。这就是单自由度无阻尼自由振动系统的能量方程。这一方程说明，这一方程说明，无阻尼自由振动系统的能量关系是无阻尼自由振动系统的能量关系是振动质体的能量与弹性势能的相互转化过程振动质体的能量与弹性势能的相互转化过程，而无而无能量的消耗。能量的消耗。但在振动系统中存在阻尼时，则在振但在振动系统中存在阻尼时，则在振动质体的动能与弹性势能的互相转化过程中，有一动质体的动能与弹性势能的互相转化过程中，有一部分能量将为克服阻力而不断地转化为热能，故系部分能量将为克服阻力而不断地转化为热能，故系统的振幅将逐渐减小，直至完全消失。统的振幅将逐渐减小，直至完全消失。 2.2 计算系统固有

24、频率的其它方法计算系统固有频率的其它方法tAxnsinEtkAtAmnnn22222sin21cos21EAmTn22max210T25232，时，、或、当nnntttEkAU2max21maxmaxUT)3(2121222kAAmnmkn若将无阻尼自由振动的时间历程若将无阻尼自由振动的时间历程 代入系统的能量方程（代入系统的能量方程（2）式可得：）式可得： 这说明系统的最大动能或最大势能均等于系统的总能量，且动能与势能的这说明系统的最大动能或最大势能均等于系统的总能量，且动能与势能的最大值相等，即：最大值相等，即：0U20时，、或、当nnttt根据上式即可算出系统的固有频率：根据上式即可算出

25、系统的固有频率：对弹簧质量系统（单振子）对弹簧质量系统（单振子）用上述能量法意义不大。但用上述能量法意义不大。但是复杂的单自由度系统用能是复杂的单自由度系统用能量法计算固有频率比较方便。量法计算固有频率比较方便。 2.2 计算系统固有频率的其它方法计算系统固有频率的其它方法例1：一根矩形截面梁，上面承受质量为m 的物体（如图所示）。若忽略梁的质量，试用能量法求该系统的固有频率。承受质量的矩形截面梁解：梁的刚度可用静变形法求出：解：梁的刚度可用静变形法求出：jmgk而梁的静扰度可根据材料力学公而梁的静扰度可根据材料力学公式计算：式计算：EJlbmgaj322223baEJlk 故代入（代入（3

26、3）式即可求出该系统的固有圆频率：）式即可求出该系统的固有圆频率：223bmaEJln2.2 计算系统固有频率的其它方法计算系统固有频率的其它方法例例2:2:下图所示为测量低频振幅用的传感器的一个元件下图所示为测量低频振幅用的传感器的一个元件无定无定向摆。已知向摆。已知a=3.54cma=3.54cm，mg=0.856Nmg=0.856N，k=0.3N/cmk=0.3N/cm。且整个系统对。且整个系统对转动轴转动轴o o的转动惯量。试求系统的固有频率。的转动惯量。试求系统的固有频率。无定向摆 解：解：取摇杆偏离平衡位置的角位移取摇杆偏离平衡位置的角位移 为广为广义坐标，并设义坐标，并设则则 对

27、简谐振动来说，摇杆正经过平衡位置时的速度对简谐振动来说，摇杆正经过平衡位置时的速度最大，故此时系统动能最大，而势能为零。即：最大，故此时系统动能最大，而势能为零。即：tAnsin)cos(tAnnnAAmaxmax2202max0max2121nAIIT 当摇杆摆到最大角位移处时，速度为零，故此时系当摇杆摆到最大角位移处时，速度为零，故此时系统动能为零，而势能最大，它包括以下两个部分：统动能为零，而势能最大，它包括以下两个部分： 1 1）弹簧变形后储存的弹性势能：）弹簧变形后储存的弹性势能： 2) 2) 质量块质量块m m的重心下降的重心下降 后的重力势能：后的重力势能：AkakaU22max

28、2max121222maxmaxmax22121cos1mglAmglmglmgU2.2 计算系统固有频率的其它方法计算系统固有频率的其它方法解：取摇杆偏离平衡位置的角位移 为广义坐标，并设则 故 对简谐振动来说，摇杆正经过平衡位置时的速度最大，故此时系统动能最大，而势能为零。即： 当摇杆摆到最大角位移处时，速度为零，故此时系统动能为零，而势能最大，它包括以下两个部分： 1）弹簧变形后储存的弹性势能：2) 质量块m的重心下降 后的重力势能：tAnsin)cos(tAnnAmaxnAmax2202max0max2121nAIITAkakaU22max2max121222maxmaxmax2212

29、1cos1mglAmglmglmgU2.2 计算系统固有频率的其它方法计算系统固有频率的其它方法022ImglkanHzImglkafn77. 0106 .174856. 054. 33 . 02212212202maxmaxUT因为 2222202121mglAAkaAIn故得2.2 计算系统固有频率的其它方法计算系统固有频率的其它方法前面介绍的几种计算系统固有频率的方法，都是将系统中弹簧的质量忽略不计。但是在有些系统中，弹簧本身的质量在系统总质量中占有一定的比例，此时若再忽略弹簧的质量，就将会使得计算出来的系统固有频率偏高。瑞利法则将弹簧质量对系统振动频率的影响考虑了进去，从而能得到相当准

30、确的固有频率值。3.瑞利法（瑞利法（Rayleigh Method）应用瑞利法时，必须先假定一个系统的振动形式。而且所假定的振动形式越接近实际的振动形式，则计算出来的固有频率的近似值就越接近准确值。实践证明，以系统的静态变形曲线作为假定的振动形式，则所求得的固有频率的近似值与准确值相比较，一般来说误差是很小的。现以最简单的弹簧质量系统为例来说明瑞利法的应用。在下图的系统中，若弹簧的质量与质量块的质量相比是很小的，则系统的振动形式就不会显著地受到弹簧质量的影响。在这种情况下，假设弹簧在振动过程中的变形（各截面的瞬时位移）与弹簧在受轴向静载荷作用下的变形相同是足够精确的。 2.2 计算系统固有频率

31、的其它方法计算系统固有频率的其它方法2.2 计算系统固有频率的其它方法计算系统固有频率的其它方法 弹簧质量系统lxx解:假设弹簧上距固定端距离为 处的位移为：x式中：l处于平衡位置时弹簧的长度； x 弹簧在联结质量块一端的位移。 令令表示弹簧单位长度的质量，则表示弹簧单位长度的质量，则弹簧微段弹簧微段dd的质量为的质量为d.d.其最其最大动能则为大动能则为: : dlx2max21lx x 弹簧在弹簧在处的微段处的微段dd的速度应为的速度应为: :当质量块在某一瞬时的速度为当质量块在某一瞬时的速度为 时，时，所以弹簧的全部动能为：所以弹簧的全部动能为：32212max20maxlxdlxTls

32、32212max20maxlxdlxTls) 1 (3232212max2max2maxmaxlmxlxxmT)2(22maxmaxkxU2.2 计算系统固有频率的其它方法计算系统固有频率的其它方法显然，系统的全部动能应该是质量块显然，系统的全部动能应该是质量块m m的最大动能与弹簧的最大的最大动能与弹簧的最大动能之和，即动能之和，即系统的最大势能仍与无质量弹簧的情况相同，即：系统的最大势能仍与无质量弹簧的情况相同，即：所以弹簧的全部动能为：所以弹簧的全部动能为：由动能和势能相等原理得：由动能和势能相等原理得：2322max2maxkxlmx对简谐振动来说，上式即成为：对简谐振动来说，上式即成

33、为：232222kAlmAn由此可以得出系统固有频率的计算公式为：由此可以得出系统固有频率的计算公式为：3lmkn结论：结论：为了考虑弹簧质量对系统固有频率的影响，只需要将为了考虑弹簧质量对系统固有频率的影响，只需要将1/3的弹簧质量当作一个集中质量加到质量块上去即可。的弹簧质量当作一个集中质量加到质量块上去即可。 一般将上式中的 称为“弹簧的等效质量”“effective mass of spring”,以ms表示。但是不同的振动系统，其弹簧的等效质量不同，需具体加以计算。 因为 所以 因此只要先算出系统弹性元件的动能，即可根据上式计算出系统弹性元件的等效质量。根据系统中的弹簧质量与质量块质

34、量相比很小，从而在振动过程中弹簧各截面的瞬时位移按线性变化这一假设而得出的。但是，即使弹簧的质量较大，用原式计算系统固有频率也具有足够的精确度。例如,当 时，固有频率的计算误差约为0.5；当 时，计算误差约为0.8；当 时，计算误差约为3。3l221xmTss22xTmssml5 . 0ml ml22.2 计算系统固有频率的其它方法计算系统固有频率的其它方法例如图所示的等截面简支梁上有一集中质量m，若将梁本身的重量W考虑在内，计算此系统的固有频率。图承受集中质量的等截面梁 2.2 计算系统固有频率的其它方法计算系统固有频率的其它方法解：假设梁在振动时挠度曲线与梁在图示载荷作用下的静挠度曲线一致

35、。梁上物体左侧距A点为处的静挠度为：梁上物体右侧距B点为处的静挠度为：在物体m处梁的静挠度为：设物体m在振动状态下的最大速度为 ，则在物体左右两侧梁的所有点的最大速度 、 与振动位移y1、y2之间存在以下关系： 216blaEJlmgby226albEJlmgayEJlbmgaym322my 1y 2y nmmyyyy11nmmyyyy22计算系统固有频率的其它方法计算系统固有频率的其它方法所以梁的左右两部分的最大速度为：因而梁的左右两部分的最大动能为：式中：w梁的单位长度的质量；22112blabayyyyymmm22222albabyyyyymmmdblabagywTams20224221

36、422222221581052332balbablgwaymgwaym22 bmsdalbbagywT02242222422222221028122aalbabaalgwbymgwbym22 2.2 计算系统固有频率的其它方法计算系统固有频率的其它方法22222158105233balbabl22222102812aalbabaal梁的全部动能为：根据上式可算出梁的等效质量为：所以系统的固有圆频率为：式中： ，为梁的刚度。 2212msssygwbwaTTTgwbwamssnmmk223babwawmggEJl 223baEJlk 2.2 计算系统固有频率的其它方法计算系统固有频率的其它方法4

37、632. 02wlEJgn4637. 02wlEJln2.2 计算系统固有频率的其它方法计算系统固有频率的其它方法从上式可以看出当忽略梁的质量时所计算出的系统固有频率比从上式可以看出当忽略梁的质量时所计算出的系统固有频率比用瑞利法计算出的数值要小，因而误差较大。应用瑞利法也可用瑞利法计算出的数值要小，因而误差较大。应用瑞利法也可求得无载荷的固有频率的相当准确的数值。由于无载荷的变形求得无载荷的固有频率的相当准确的数值。由于无载荷的变形曲线是对称的，所以首先需将载荷移到梁的中间，然后再令载曲线是对称的，所以首先需将载荷移到梁的中间，然后再令载荷为零（荷为零（m m0 0），即可求出无载荷梁的固有

38、圆频率为：），即可求出无载荷梁的固有圆频率为：而这一固有圆频率的精确值为：而这一固有圆频率的精确值为：可见，近似值与理论精确值之差小于可见，近似值与理论精确值之差小于1 1。内容参考2.3。2.3 等效质量和等效刚度等效质量和等效刚度振动微分方程振动微分方程 振动微分方程振动微分方程 )(tFkxxcxm 方程的解方程的解 )()()(21txtxtx其其中中， tx1为为相相应应齐齐次次方方程程的的解解 瞬瞬态态响响应应 tx2为为方方程程的的特特解解 稳稳态态响响应应或或零零初初始始条条件件的的解解 单自由度线性阻尼自由振动系统单自由度线性阻尼自由振动系统自由振动自由振动 振动微分方程振动

39、微分方程设设 0 xkxcxm tsAtxe)(02kscsm特征方程特征方程 222, 1nnns有有临界阻尼系数临界阻尼系数 kmc2c阻尼比或阻尼因子阻尼比或阻尼因子 kmccc2c定义定义12nn2, 1smkmcnnnn2;2;令阻尼比或阻尼因子令阻尼比或阻尼因子 单自由度线性阻尼自由振动系统单自由度线性阻尼自由振动系统讨论讨论 (1)系统无阻尼即,0,0n方程的解方程的解 12nn2, 1s特征值特征值系统对初始扰动的响应系统对初始扰动的响应）（tRtxncos)(2n020)/(xxR0tanarc0tanarc0n000n00 xxxxxx11100，特征值取决于单自由度线性阻

40、尼自由振动系统单自由度线性阻尼自由振动系统讨论讨论 (2)12nn2, 1s特征值特征值系统对初始扰动的响应系统对初始扰动的响应方程的解方程的解 个不等的虚数根。为2,102, 1Snn)(cose)(rntRtxt2d0n0202221xxxRDDR0tanarc0tanarc00d0n000d0n0 xxxxxxxx则令,22rnn)sincos(sincosr2r1rrr2, 1rtDtDextiteinsnti221122nrnrrT单自由度线性阻尼自由振动系统单自由度线性阻尼自由振动系统讨论讨论 (3)方程的解方程的解 12nn2, 1s特征值特征值系统对初始扰动的响应系统对初始扰动

41、的响应nssnn21,1nttCCtxe )(21)(e)(000txnxxtxtn00 xx）（00 xx）（初始条件：初始条件：单自由度线性阻尼自由振动系统单自由度线性阻尼自由振动系统讨论讨论 (4)12nn2, 1s特征值特征值系统对初始扰动的响应系统对初始扰动的响应tstsxsxsxxsstx21e)(e)(1)(01020021tstsCCtx21ee)(2100 xx）（00 xx）（1方程的解方程的解 单自由度线性阻尼自由振动系统单自由度线性阻尼自由振动系统振动特性振动特性 无阻尼无阻尼 0 0： 简谐运动简谐运动弱阻尼弱阻尼 0 1： 衰减运动衰减运动单自由度线性阻尼自由振动系

42、统单自由度线性阻尼自由振动系统小阻尼小阻尼振动对数衰减率振动对数衰减率 2112lnnnxx2224单自由度线性阻尼自由振动系统单自由度线性阻尼自由振动系统简谐激励简谐激励稳态响应稳态响应(粘性阻尼粘性阻尼)tFxkxcxmsin0 M M 2. 5 简谐激励力作用下的强迫振动简谐激励力作用下的强迫振动 求解过程求解过程 运动方程的解可以用它对应的齐次方程的通解运动方程的解可以用它对应的齐次方程的通解 和方程（和方程（2）的特解）的特解 来表示来表示 )2(sin2,2) 1 (sinsin20n00 tfxnxmFfmkmcntmFmkxmcxtFxkxcxmn则，令可化为2x)()(21t

43、xtxx1x)sin()(1tAetxrnt 在小阻尼情况下，在小阻尼情况下， 是个衰减振动，只在开始振动是个衰减振动，只在开始振动后的某一段时间内有意义。研究受迫振动中持续等幅振动时可忽略之。后的某一段时间内有意义。研究受迫振动中持续等幅振动时可忽略之。 表示系统的受迫振动，称为系统的表示系统的受迫振动，称为系统的稳态解稳态解，设，设2x)sin()(2tBtx将将 代入到方程（代入到方程（2）中可解出）中可解出B与与2x22222222arctan;4)(nnnnfB2222012arctan;)2()1 (1kFBnn;n令2. 5 简谐激励力作用下的强迫振动简谐激励力作用下的强迫振动

44、求解过程求解过程 进一步讨论：进一步讨论： kFBs0sBB令：令： 则：则： 222)2()1 (1sBB2222012arctan;)2()1 (1kFB2. 5 简谐激励力作用下的强迫振动简谐激励力作用下的强迫振动 解的讨论解的讨论二、讨论二、讨论: 图给出了以图给出了以为横坐标，为横坐标，为纵坐标，在不同阻尼比为纵坐标，在不同阻尼比下的一组曲线簇。不难理解，在简谐激振力作用下，线性下的一组曲线簇。不难理解，在简谐激振力作用下，线性系统的受迫振动也是简谐振动，振动的频率等于激励力的系统的受迫振动也是简谐振动，振动的频率等于激励力的频率，受迫振动的振幅取决于系统本身的物理特性、激励频率，受

45、迫振动的振幅取决于系统本身的物理特性、激励力的大小及频率值，但与初始条件无关。力的大小及频率值，但与初始条件无关。 受迫振动的振幅与频率比及阻尼比有关受迫振动的振幅与频率比及阻尼比有关 (1) 当频率比当频率比0.2时，即激振频率时，即激振频率远小于系统的固有频远小于系统的固有频率率n时，无论阻尼的大小如何，时，无论阻尼的大小如何，1，称为准静态区。即，称为准静态区。即振幅近似等于激励力幅作用下的静变形。故在低频区振幅振幅近似等于激励力幅作用下的静变形。故在低频区振幅主要由弹簧刚度控制。主要由弹簧刚度控制。2. 5 简谐激励力作用下的强迫振动简谐激励力作用下的强迫振动 解的讨论解的讨论（2）频

46、率比很大频率比很大(5) , 0，激振频率，激振频率远大于系统的固有远大于系统的固有频率频率n ，因激励力方向改变太快，振动物体由于惯性来不，因激励力方向改变太快，振动物体由于惯性来不及跟随，几乎停着不动。故在及跟随，几乎停着不动。故在高频区受迫振动的振幅主要取高频区受迫振动的振幅主要取决于系统的惯性，决于系统的惯性，称为惯性区，这一特性正是隔振和惯性传称为惯性区，这一特性正是隔振和惯性传感器的理论依据。感器的理论依据。（3）当频率比当频率比 =1，激振频率接近系统的固有频率，这时阻尼值越小，激振频率接近系统的固有频率，这时阻尼值越小， 则越大。当阻尼为零时，振动为无限大。习惯上把幅值则越大。

47、当阻尼为零时，振动为无限大。习惯上把幅值 的频率的频率区间称为共振区。区间称为共振区。 将（将（6）对求导，并令）对求导，并令d/d=0 ，可解得，可解得 处有最大幅值，处有最大幅值，把把 称为共振频率。称为共振频率。 2221221n2. 5 简谐激励力作用下的强迫振动简谐激励力作用下的强迫振动 解的讨论解的讨论 相位相位 与频率比的关系曲线表明与频率比的关系曲线表明 =1时，振动位时，振动位移总是滞后激振力移总是滞后激振力/2 ，频率比，频率比 1；当；当 =-/2 -，共振点前后相位差，共振点前后相位差恰好为恰好为。 2. 5 简谐激励力作用下的强迫振动简谐激励力作用下的强迫振动 简谐激

48、励简谐激励全响应全响应(粘性阻尼粘性阻尼)tFxkxcxmsin0 222012arctansin)2()1 (1)cos(Re)(tkFttxrtn2. 5 简谐激励力作用下的强迫振动简谐激励力作用下的强迫振动 简谐激励简谐激励全响应全响应(无无阻尼阻尼)：)1 (sin0tFxkxm 设其特解为：设其特解为：tBtxsin)(22nfB代入到上式得：代入到上式得：)2(sinsincos2221tftCtCxnnn方程方程（1）的通解解为：的通解解为： 0000 xxxx设初始条件为：设初始条件为：220201/nnfxCxC代入到方程（代入到方程（2）中得：）中得：)3()sin(sin

49、sincos2200ttftxtxxnnnnnn则：则：)4()sin(sin)sin(220ttftAxnnnn即：即：初始条件产生的自由振动初始条件产生的自由振动简谐激励力产生的受迫振动简谐激励力产生的受迫振动伴随受迫振动产生的自由振动伴随受迫振动产生的自由振动2. 5 简谐激励力作用下的强迫振动简谐激励力作用下的强迫振动 ttkFtxnn20sinsin1)(n 000 xx若初始条件为：若初始条件为：则：则： 2. 5 简谐激励力作用下的强迫振动简谐激励力作用下的强迫振动 2. 5 简谐激励力作用下的强迫振动简谐激励力作用下的强迫振动 简谐激励简谐激励全响应全响应(无无阻尼阻尼)tFx

50、kxmsin0 000 xxttkFtxnn20sinsin1)(nnn2.6 简谐力的功和等效阻尼简谐力的功和等效阻尼简谐力的功简谐力的功简谐力简谐力tFtFsin)(0dtxtFFdxdWsin0Q=)sin()(tBtx振动系统的稳态解为振动系统的稳态解为则激振力在微小位移则激振力在微小位移dxdx上所作的微元功应为：上所作的微元功应为：在一个周期内（在一个周期内（t=0t=02/w2/w）所作的功，也就是）所作的功，也就是F(t)F(t)输入输入系统的能量，即为系统的能量，即为dtxtFW20)(dtttBF)cos(sin200tdttBF)cos(sin200sin0BF可见，简谐

51、激振力在一个周期内所作功的大小，不仅可见，简谐激振力在一个周期内所作功的大小，不仅决定于激振力幅决定于激振力幅F F0 0 及振幅及振幅 B B 的大小，还决定于两者的大小，还决定于两者之间的相位角之间的相位角 。 当当00即外力超前位移时，作正功；即外力超前位移时，作正功；当当00即外力落后于位移时，作负功；即外力落后于位移时，作负功；而当而当 =0=0或或 =时，即外力在一个周期内作功之和时，即外力在一个周期内作功之和等于零。等于零。激振力在一个周期内所作的功激振力在一个周期内所作的功W W ，可以看成是激振，可以看成是激振力的两个分量作功的和，即与位移同相的分量力的两个分量作功的和，即与

52、位移同相的分量F = F F = F coscos和与速度同相的分量和与速度同相的分量F = F sinF = F sin所作功之和。所作功之和。sin0BFW 2.6 简谐力的功和等效阻尼简谐力的功和等效阻尼0)cos()sin(cos)cos()sin(cos)sin(200200201tdttBFdttBtFdtxtFWF)sin(cos0tF与位移相同的力：与位移相同的力：在一个周期内所作的功为：在一个周期内所作的功为：2.6 简谐力的功和等效阻尼简谐力的功和等效阻尼简谐力简谐力tFtFsin)(0激振力在一个周期内所作的功为分量作功之和，即为激振力在一个周期内所作的功为分量作功之和，

53、即为 W =W +W = F Bsin因此，激振力在一个周期内所作的功，就是其超前位移因此，激振力在一个周期内所作的功，就是其超前位移 /2 的分量所作的功。的分量所作的功。sin)(cossin)cos(02200202 BFtdtBFdtxtFWF与速度同向的力与速度同向的力F sincos（t-）在一个周期内所作的功为：在一个周期内所作的功为：2.6 简谐力的功和等效阻尼简谐力的功和等效阻尼简谐力简谐力tFtFsin)(02.6 简谐力的功和等效阻尼简谐力的功和等效阻尼时，粘性阻尼力时，粘性阻尼力 在一个振动周期中所做的功：在一个振动周期中所做的功：xcFc2202220)(sinXcd

54、ttcXdtxFDTcc在单自由度受迫振动方程中，阻尼力被设为在单自由度受迫振动方程中，阻尼力被设为 。 实际物理模实际物理模型与振动位移一阶导数成正比的是纯液体摩擦阻尼，称为型与振动位移一阶导数成正比的是纯液体摩擦阻尼，称为粘性粘性阻尼阻尼。这种阻尼是。这种阻尼是线性线性的，数学上易于处理，故常把非线性阻的，数学上易于处理，故常把非线性阻尼用等效粘性阻尼来代替。尼用等效粘性阻尼来代替。等效原则：一个振动周期中，两种阻尼耗散的能量相等。等效原则：一个振动周期中，两种阻尼耗散的能量相等。xc等效阻尼力等效阻尼力 在一个振动周期中所作的功：在一个振动周期中所作的功： 所以有：所以有： xce2Xc

55、DDec) 1 (2XDce)cos(tXx当受迫振动的位移响应为：当受迫振动的位移响应为： 干摩擦阻尼：干摩擦阻尼：干摩擦阻尼力干摩擦阻尼力F可视为一个常力，在整个受迫可视为一个常力，在整个受迫振动中力的幅值不变，方向始终与运动方向相反。振动中力的幅值不变，方向始终与运动方向相反。当质量从平衡位置移动到最大偏离位置当质量从平衡位置移动到最大偏离位置X，即在周期内，摩，即在周期内，摩擦力做功为擦力做功为 FX，故一个整周期内做功，故一个整周期内做功 代入代入（1）式，得到式，得到干摩擦的等效阻尼干摩擦的等效阻尼:XFXFXce4422.6 简谐力的功和等效阻尼简谐力的功和等效阻尼结构阻尼：结构

56、阻尼：由材料形变过程中的内摩擦产生。由材料形变过程中的内摩擦产生。材料在加载材料在加载- -卸载过程中，会形成应力卸载过程中，会形成应力- -应变迟滞曲线，它应变迟滞曲线，它包容的面积就是内摩擦所消耗的能量，它近似地与振幅平包容的面积就是内摩擦所消耗的能量，它近似地与振幅平方成正比。即：方成正比。即：其中其中 是与频率无关的比例系数，随材料不同而变。因是与频率无关的比例系数，随材料不同而变。因此，结构等效阻尼：此，结构等效阻尼：2XaDrrreaXXaXDc222ra2.6 简谐力的功和等效阻尼简谐力的功和等效阻尼周期激励周期激励稳态响应稳态响应(粘性阻尼粘性阻尼) 122211021sinc

57、os2nnnnnnnktnbtnakatx212arctannnnn1mknkmc221 1110sincos2nnntnbtnaatFtF tFxkxcxm 2. 7 非简谐周期激励力的响应分析非简谐周期激励力的响应分析 例：例：质量弹簧系统受到周期方波激励质量弹簧系统受到周期方波激励 求系统稳态响应。求系统稳态响应。 TtTFTtFtF2,20,)(00)(tF0F0F0T2/Tt610 1 . 00 2. 7 非简谐周期激励力的响应分析非简谐周期激励力的响应分析 解：解：0弹簧质量系统固有频率弹簧质量系统固有频率 610激励力的基频激励力的基频 ：T26011110)sincos(2)(

58、nnntnbtnaatFTnTnTtdtntFTbtdtntFTadttFTa110sin)(2cos)(2)(2因因 a0 一周期内总面积为一周期内总面积为0 =0T, 02Ttn1cos2T区间区间 内，内， 关于关于 为反为反对称，对称， 而而 关于关于 对称对称)(tF=0 11sinnntnb)(tF0F0F0T2/Tt2. 7 非简谐周期激励力的响应分析非简谐周期激励力的响应分析 11sin)(nntnbtF TntdtntFTb1sin)(22, 0T区间区间 内内)(tF4T关于关于为对称为对称 tn1sin4T而而n取偶数时，取偶数时，关于关于反对称反对称 ,2TT区间区间

59、内内)(tF43T关于关于为对称为对称 tn1sin43T而而n取偶数时，取偶数时，关于关于反对称反对称 0nb6 , 4 , 2n因此因此)(tF0F0F0T2/Tt12T2. 7 非简谐周期激励力的响应分析非简谐周期激励力的响应分析 11sin)(nntnbtF TntdtntFTb1sin)(2当当 n 取奇数时取奇数时 0nb6 , 4 , 2nnFtdtnFTtdtntFTbTTn04010014sin8sin)(2 5 , 3 , 1n于是，周期性激励于是，周期性激励 F(t) 可写为：可写为： 5 , 3 , 11011sin14sin)(nnntnnFtnbtF)5sin513

60、sin31(sin41110 tttF2. 7 非简谐周期激励力的响应分析非简谐周期激励力的响应分析 )5sin513sin31(sin4)(1110 tttFtF5 , 3 , 110)sin(4nnntnkFx则有：则有：2222)2()1 (1nssnnn 22112snsntgn 01s其中：其中： 当不计阻尼时：当不计阻尼时： 5 , 3 , 11220sin)1 (14ntnsnnkFx系统运动方程系统运动方程 ：)(tFkxxcxm 1110)sin()cos(2)(nnnnnnntnkbtnkakatx2. 7 非简谐周期激励力的响应分析非简谐周期激励力的响应分析 小结：小结：