第四章 连续体的振动精选文档.ppt

《第四章 连续体的振动精选文档.ppt》由会员分享，可在线阅读，更多相关《第四章 连续体的振动精选文档.ppt（52页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、第四章 连续体的振动本讲稿第一页，共五十二页本章只讨论理想弹性体的振动理想弹塑性体满足以下假设条件各向同性；均质分布；服从虎克定律4.1弦的振动讨论两端受到张力 拉紧的弦，弦上还受到横向干扰力的作用第四章 连续体的振动本讲稿第二页，共五十二页设弦的密度为 （质量/单位体积）假设小变形，弦力不随挠度变化。则弦上的任意一点的位移y应为位置x与时间t的函数，即沿方向作用在微小区间的外力之和为本讲稿第三页，共五十二页根据牛顿第二定律，弦的单元微段ds沿 方向的运动微分方程为：代入得：设代入得：C为波沿长度方向的传播速度本讲稿第四页，共五十二页如无干扰力作用时，称为波动方程 弹性体系统作某阶主振动时，其

2、各点也应当作同样的频率及相位运动，各点也应当同时通过静平衡位置和到达最大偏离位置，即系统具有一定的与时间无关的振型为振型函数本讲稿第五页，共五十二页得故 4个待定常数，可由弦的边界条件及振动的两个初始条件来确定。由于两端固定，故有本讲稿第六页，共五十二页得 则得本讲稿第七页，共五十二页强迫振动 对于长为的两端固定受分布力作用下的弦的强迫振动，其运动微分方程为：振型函数 令则有设其解为 代入方程本讲稿第八页，共五十二页将上式两边同乘以并从得：整理后得到：其通解为：本讲稿第九页，共五十二页4.2杆的纵向振动设杆的横截面在振动时仍保持为平面并作整体运动。略去杆纵向伸缩设杆的横截面在振动时仍保持为平面

3、并作整体运动。略去杆纵向伸缩而引起的横截面变形。而引起的横截面变形。取杆的纵向作为取杆的纵向作为 x 轴，各个截面的纵向位移表示为轴，各个截面的纵向位移表示为u(x,t)，如图如图所示。杆的微元杆的微元dx在自由振动中的受力图也在图在自由振动中的受力图也在图5-3中给出。中给出。设杆单位体积的质量为设杆单位体积的质量为，杆长为，杆长为 l，截面积为，截面积为A，材料的弹性模量为材料的弹性模量为E。再设任一。再设任一x截面处，纵向应变为截面处，纵向应变为(x)，纵向拉力表示为纵向拉力表示为N(x)本讲稿第十页，共五十二页4.2杆的纵向振动本讲稿第十一页，共五十二页由牛顿第二定律：当杆为均质，等截

4、面时设本讲稿第十二页，共五十二页例-1：两端固定的等值杆代入得则有 本讲稿第十三页，共五十二页各阶主振型 各阶主振动：本讲稿第十四页，共五十二页左端：右端：杆受到弹簧力代入由上式，根据不同的值，可解出不同的固有频率例5.2 如图所示，一长为 的等截面均质直杆，左端固定，右端联结一刚度为 的弹簧，试求其纵向自由振动的运动方程。本讲稿第十五页，共五十二页 相当于自由端 主振型 即相当于固定端，其频率方程为本讲稿第十六页，共五十二页4.3圆轴的扭转振动 取取圆圆轴轴的的轴轴心心线线作作为为x 轴轴，图图示轴轴任任一一 x截截面面处处的的转转角角表表示示为为(x,t)。设设轴轴长长为为l ，单单位位体

5、体积积的的质质量量为为，圆圆截截面面对对其其中中心心的的极极惯惯量量矩矩为为Jp，材材料料的的剪剪切切弹弹性性模模量量为为 G 。轴轴的的扭扭转转应应变变为为 ，作作用用于于微微元元dx 两截面上的扭矩分别为两截面上的扭矩分别为 ，及，及 。假设轴的横截面在扭转振动中仍保持为平面作整体转动。假设轴的横截面在扭转振动中仍保持为平面作整体转动。本讲稿第十七页，共五十二页4.3圆轴的扭转振动杆的扭转振动，抗扭刚度 截面处的扭转角由材料力学知识微元段扭矩的增量为极惯性矩：圆截面极惯性矩：密度的转动惯量）本讲稿第十八页，共五十二页得：设 则有：:表示剪切波在杆内的传播速度它的解为式中四个待定常数 由系统

6、的边界条件和初始条件确定。一般解为：本讲稿第十九页，共五十二页解：边界条件：对于圆盘的运动微分方程：例-3：上端固定，下端装有转动惯量为的圆盘，圆轴的极惯性矩为JP，其剪切弹性模量为G，试分析其扭转振动的频率方程。JPIPJJ本讲稿第二十页，共五十二页由于边界条件得：B=0简化得：得：频率方程代入得：的物理意义为轴的转动惯量与圆盘转动惯量之比。的物理意义为轴的转动惯量与圆盘转动惯量之比。设其中(a）本讲稿第二十一页，共五十二页如近似地取如近似地取，则（，则（a a）式化简为）式化简为下表给出对应于各个不同的下表给出对应于各个不同的 值时，基本值时，基本特征值特征值1的值。的值。0.010.01

7、0.100.100.300.300.500.500.700.700.900.901.001.001.501.500.100.100.320.320.520.520.650.650.750.750.820.820.860.860.980.982.002.003.003.004.004.005.005.0010.010.020.020.01001001.081.081.201.201.271.271.321.321.421.421.521.521.571.57（b b）上式也就是略去轴的质量后所得单自由度系统的上式也就是略去轴的质量后所得单自由度系统的固有频率公式固有频率公式。本讲稿第二十二页，共

8、五十二页这时有这时有（c c）上式也就是将轴的转动惯量的三分之一加到圆盘上后所得单自由上式也就是将轴的转动惯量的三分之一加到圆盘上后所得单自由度扭振系统的固有频率公式。只要轴的转动惯量不大于圆盘的转动惯度扭振系统的固有频率公式。只要轴的转动惯量不大于圆盘的转动惯量，那么计算基频的近似式（量，那么计算基频的近似式（c c）在实用上已足够准确了。）在实用上已足够准确了。综上所述，弦的横振、杆的纵振与轴的扭振都导致同一形式的波综上所述，弦的横振、杆的纵振与轴的扭振都导致同一形式的波动方程。它们的运动具有共同的规律，如表动方程。它们的运动具有共同的规律，如表4-14-1。进一步的近似可取进一步的近似可

9、取本讲稿第二十三页，共五十二页弦的横振弦的横振杆的纵振杆的纵振轴的扭振轴的扭振物物理理参参数数弦的张力弦的张力弦的体密度弦的体密度弹性模量弹性模量截面积截面积密度密度剪切弹性模量剪切弹性模量截面极惯性矩截面极惯性矩密度密度截面的截面的位移位移横向位移横向位移纵向位移纵向位移转转角角单位长度单位长度的质量或的质量或转动惯量转动惯量截面处力截面处力（或扭矩）（或扭矩）表表4-14-1弦的横振、杆的纵振与轴的扭振对比表弦的横振、杆的纵振与轴的扭振对比表本讲稿第二十四页，共五十二页弦的横振弦的横振杆的纵振杆的纵振轴的扭振轴的扭振波速波速运动微运动微分方程分方程通解通解边界条件边界条件两端固定两端固定两

10、端自由两端自由一端固定一端固定一端自由一端自由固有频率固有频率振型函数振型函数本讲稿第二十五页，共五十二页作业：p228:4-1;4-5本讲稿第二十六页，共五十二页4.4梁的弯曲振动梁的振动都在同一平面内，即在xoy平面之内满足平截面假设梁上各点的运动只需用轴线的横向位移来表示建立运动微分方程：按牛顿第二定律，微段沿方向的运动方程为：.（1）本讲稿第二十七页，共五十二页由材料力学可知：代入（1）式中得：则有：为梁弯曲振动的偏微分方程梁自由振动的解分离变量：代入得：本讲稿第二十八页，共五十二页设 则有，代入设 代入 的通解为：本讲稿第二十九页，共五十二页则：本讲稿第三十页，共五十二页求固有频率和

11、主振型例-5因此有：本讲稿第三十一页，共五十二页代入B=D=0则有：本讲稿第三十二页，共五十二页例-6：悬臂梁的情况 参见书上的例题。本讲稿第三十三页，共五十二页4.5 主振型的正交性对于等截面的均匀梁，其截面面积，线密度刚度分布都是常数时，有对于一般情况，各阶主振型正交，表明各阶主振动之间没有耦合，各主振动之间是独立无关的，它们之间没有能量交换，都可以采用解耦分析法（振型分析法），简化为类似单自由度系统那样的微分方程来求解。作业：作业：p228:4-8；4-11本讲稿第三十四页，共五十二页4.6薄板的横向弯曲振动基本假设：1直法线假设；2为主应力；3厚度变化不计；4板弯曲的中面不产生变形，即

12、中面为中性面。一、建模本讲稿第三十五页，共五十二页矩形薄板的振动微分方程惯性载荷：二、求解-分离变量法设：得：本讲稿第三十六页，共五十二页三、边界条件1.简支边：2固定边3自由边本讲稿第三十七页，共五十二页四、四边简支板自由振动的解设：固有频率本讲稿第三十八页，共五十二页结构动力学分析工程应用介绍结构动力学分析工程应用介绍 随着形形色色飞行器的出现，动力固有特性分析也出现许多特殊的随着形形色色飞行器的出现，动力固有特性分析也出现许多特殊的问题。例如，火箭、导弹虽然基本上近似一个一维梁系统，但它们往往问题。例如，火箭、导弹虽然基本上近似一个一维梁系统，但它们往往存在着级间连接与数目较多的分离面存

13、在着级间连接与数目较多的分离面(接头接头),这使得结构的刚度分布、阻这使得结构的刚度分布、阻尼分布发生变化，给固有特性的精确分析带来困难；大型液体燃料火箭尼分布发生变化，给固有特性的精确分析带来困难；大型液体燃料火箭广泛存在着液体晃动问题，它对固有特性有不可忽视的影响；人造卫星广泛存在着液体晃动问题，它对固有特性有不可忽视的影响；人造卫星经常采用自旋稳定方案，并安装有柔度很大的附件，如太阳能帆板、天经常采用自旋稳定方案，并安装有柔度很大的附件，如太阳能帆板、天线等，构成具有转动的刚体及柔性体的复合结构，使它的固有特性分析线等，构成具有转动的刚体及柔性体的复合结构，使它的固有特性分析更加复杂化。

14、更加复杂化。本章对这些问题做一概括介绍。本章对这些问题做一概括介绍。本讲稿第三十九页，共五十二页一、一、飞行器分离面（接头）对固有特性的影响飞行器分离面（接头）对固有特性的影响 大部分火箭、导弹类型的飞行器都具有较多的接头，它们使全弹的刚大部分火箭、导弹类型的飞行器都具有较多的接头，它们使全弹的刚度分布发生局部扰动，如图所示。度分布发生局部扰动，如图所示。本讲稿第四十页，共五十二页一、一、飞行器分离面（接头）对固有特性的影响飞行器分离面（接头）对固有特性的影响导导弹弹名名称称接接头头数数刚刚度度损损失失K注注“响尾蛇”413近程攻击导弹SRAM631中程攻击导弹MR731远程标准导弹ER733

15、“不死鸟”1049接头使刚度损失可达接头使刚度损失可达(3040)%,如表如表81所示所示表表 8 81 1本讲稿第四十一页，共五十二页一、一、飞行器分离面（接头）对固有特性的影响飞行器分离面（接头）对固有特性的影响 对弹体固有特性的影响对弹体固有特性的影响1.1.使全弹的固有频率下降使全弹的固有频率下降 由于接头削弱了附近弹体的刚度，同时接头处往往存在空隙，因而它们都使全弹的固由于接头削弱了附近弹体的刚度，同时接头处往往存在空隙，因而它们都使全弹的固有频率降低。影响程度与接头的数量、类型、位置相关。表有频率降低。影响程度与接头的数量、类型、位置相关。表8-2中列举了一些导弹的一阶中列举了一些

16、导弹的一阶频率由于接头而引起的下降情况。频率由于接头而引起的下降情况。导导弹弹名名称称接接头头数数一阶固有频率下降一阶固有频率下降“响尾蛇”47近程攻击导弹SRAM617中程攻击导弹MR717远程标准导弹ER718“不死鸟”1033表表 8-2 8-2本讲稿第四十二页，共五十二页一、一、飞行器分离面（接头）对固有特性的影响飞行器分离面（接头）对固有特性的影响 对弹体固有特性的影响对弹体固有特性的影响2.2.使全弹的振动发生畸形使全弹的振动发生畸形 接头的存在，改变了刚度分布，必然使振型形状、节点位置发生变化。在控制系统的设计中节点位接头的存在，改变了刚度分布，必然使振型形状、节点位置发生变化。

17、在控制系统的设计中节点位置是个重要参数，所以，为了精确确定振型，必须考虑接头的响应。图置是个重要参数，所以，为了精确确定振型，必须考虑接头的响应。图8-2所示为某弹前三阶振型受接所示为某弹前三阶振型受接头影响的变换情况。头影响的变换情况。图中图中x为弹体轴向坐标，坐标原点设在弹体头部理论顶点，为弹体轴向坐标，坐标原点设在弹体头部理论顶点，为振型幅值。为振型幅值。应当注意，同样的接头，所处的位置不同影响也不同。一般来讲，在导弹弹应当注意，同样的接头，所处的位置不同影响也不同。一般来讲，在导弹弹体中部的接头影响更加突出。体中部的接头影响更加突出。本讲稿第四十三页，共五十二页一、一、飞行器分离面（接

18、头）对固有特性的影响飞行器分离面（接头）对固有特性的影响 分析方法分析方法 精确地用纯分析方法考虑接头进行固有特性计算是困难的，一般都采用实验与分析结合的精确地用纯分析方法考虑接头进行固有特性计算是困难的，一般都采用实验与分析结合的方法。全弹的计算模型可以选用一维梁式模型或三维壳体模型，接头则可分为处理为集中弯曲方法。全弹的计算模型可以选用一维梁式模型或三维壳体模型，接头则可分为处理为集中弯曲弹簧或沿分离面周线分布的弹性组件。处理这类模型的弹簧或沿分离面周线分布的弹性组件。处理这类模型的关键是这些弹性件柔度的确定关键是这些弹性件柔度的确定。对于。对于集中弯度弹簧，其柔度确定的方法有以下几种。集

19、中弯度弹簧，其柔度确定的方法有以下几种。1.1.类比法类比法 参照已有的导弹接头的柔度数据，用相似类比方法，推测所设计的接头的柔度。参照已有的导弹接头的柔度数据，用相似类比方法，推测所设计的接头的柔度。2.2.经验公式法经验公式法Alley和和Leadbetfer根据大量实验的统计，归纳出以下经验公式根据大量实验的统计，归纳出以下经验公式式中式中 第第i个接头的弯曲柔度；个接头的弯曲柔度；接头处弹身直径（英寸）接头处弹身直径（英寸）柔度系数，根据接头的不同类型已制成表格供查（表柔度系数，根据接头的不同类型已制成表格供查（表8-3）本讲稿第四十四页，共五十二页一、一、飞行器分离面（接头）对固有特

20、性的影响飞行器分离面（接头）对固有特性的影响接头分类Ai额定值范围优10-10310-10良10-9310-10310-9中10-8310-9310-8差10-7310-8310-7表表 9-3此经验公式是以英制给出的，使用时应予注意。此经验公式是以英制给出的，使用时应予注意。3.3.实验测定法实验测定法4.4.其它方法其它方法 当已具有实体结构时，可通过静力实验或动力实验得到各个接头实际柔度。不当已具有实体结构时，可通过静力实验或动力实验得到各个接头实际柔度。不过，实验中应注意消除弹性弯曲的影响。过，实验中应注意消除弹性弯曲的影响。也可采用有限元或最佳拟合的分析方法来确定接头的柔度，不过方法

21、复杂而精度并不也可采用有限元或最佳拟合的分析方法来确定接头的柔度，不过方法复杂而精度并不理想。理想。分析方法分析方法本讲稿第四十五页，共五十二页一、一、飞行器分离面（接头）对固有特性的影响飞行器分离面（接头）对固有特性的影响 总之，随着飞行器对固有频率、振型、振型斜率的数量与精度要求日益提高，接头产总之，随着飞行器对固有频率、振型、振型斜率的数量与精度要求日益提高，接头产生的影响必须予以考虑。由于接头类型较多，单纯的分析方法尚未完善，目前主要依靠实生的影响必须予以考虑。由于接头类型较多，单纯的分析方法尚未完善，目前主要依靠实验来确定特性参数验来确定特性参数柔度柔度。将所得柔度参量代入系统模型，

22、即可计得较为精。将所得柔度参量代入系统模型，即可计得较为精确得全弹固有特性。确得全弹固有特性。分析方法分析方法本讲稿第四十六页，共五十二页 问题的特点问题的特点二、二、贮箱内液体晃动对固有特性的影响贮箱内液体晃动对固有特性的影响晃动频率及晃动激烈程度均与下列因素有关：晃动频率及晃动激烈程度均与下列因素有关：(1)容器形状；容器形状；(2)推进器性质；推进器性质；(3)阻尼隔板设置情况；阻尼隔板设置情况；(4)推进剂液面高度；推进剂液面高度；(5)加速度场的情况。加速度场的情况。在飞行过程中，随着燃料不断燃烧，推进剂液面情况不断发生变化，从而对系统固有特性的影响也随之在飞行过程中，随着燃料不断燃

23、烧，推进剂液面情况不断发生变化，从而对系统固有特性的影响也随之变化。这是本问题的主要特点。一般说来，为了掌握整个飞行过程的固有特性情况，就要分析各个不同的推变化。这是本问题的主要特点。一般说来，为了掌握整个飞行过程的固有特性情况，就要分析各个不同的推进剂燃烧阶段的频率与振型。进剂燃烧阶段的频率与振型。另外一个特点是，一般只需考虑低阶情况，特别是一阶情况。因为经研究指出，对于圆柱壳体，二阶晃动另外一个特点是，一般只需考虑低阶情况，特别是一阶情况。因为经研究指出，对于圆柱壳体，二阶晃动质量仅为一阶晃动质量的质量仅为一阶晃动质量的3，而且在高阶情况下液体内部将产生紊乱的扰动，使阻尼激增，故二阶以上可

24、不，而且在高阶情况下液体内部将产生紊乱的扰动，使阻尼激增，故二阶以上可不予考虑。予考虑。本讲稿第四十七页，共五十二页 推进剂晃动频率的确定推进剂晃动频率的确定二、二、贮箱内液体晃动对固有特性的影响贮箱内液体晃动对固有特性的影响 工程上常采用当量变换的方法，在对壳体壁作用的力与力矩相等、频率相当的条件下，将液体晃工程上常采用当量变换的方法，在对壳体壁作用的力与力矩相等、频率相当的条件下，将液体晃动模型等价代换为机械力学模型。一旦建立了当量机械模型，对于各种飞行器的液体晃动问题，可根动模型等价代换为机械力学模型。一旦建立了当量机械模型，对于各种飞行器的液体晃动问题，可根据它们的液体参量、飞行状态参

25、量、飞行器参量很容易地确定出当量机械模型参量，从而确定出晃动据它们的液体参量、飞行状态参量、飞行器参量很容易地确定出当量机械模型参量，从而确定出晃动频率。一般采用的液体侧向晃动当量机械模型有两种。频率。一般采用的液体侧向晃动当量机械模型有两种。1.1.弹簧质量模型弹簧质量模型 根据壳体半径及液面高度，可按流体动力学分析根据壳体半径及液面高度，可按流体动力学分析导得的当量公式确定贮箱中液体固定质量导得的当量公式确定贮箱中液体固定质量m0、晃动质、晃动质量量m1。晃动质量的运动受到弹簧与阻尼器的约束，。晃动质量的运动受到弹簧与阻尼器的约束，其模型如图其模型如图8-10所示。所示。模型中阻尼系数模型

26、中阻尼系数C、弹簧刚度、弹簧刚度K1都按流体动力学导得的都按流体动力学导得的当量公式确定。当阻尼较小时，阻尼的作用可以忽略。当量公式确定。当阻尼较小时，阻尼的作用可以忽略。本讲稿第四十八页，共五十二页 推进剂晃动频率的确定推进剂晃动频率的确定二、二、贮箱内液体晃动对固有特性的影响贮箱内液体晃动对固有特性的影响2.2.自由摆模型自由摆模型 从晃动的物理现象来看，将它等价为一个当量摆是无可非议的。此模型如图从晃动的物理现象来看，将它等价为一个当量摆是无可非议的。此模型如图8-11所示。其中所示。其中 等参量也是由流体动力学分析得到的当量公式来确定。等参量也是由流体动力学分析得到的当量公式来确定。进

27、行这种模型代换后，液体晃动问题的处理就与结构系统的其它部件的处理方式完全相同。考虑到将它进行这种模型代换后，液体晃动问题的处理就与结构系统的其它部件的处理方式完全相同。考虑到将它们并入全系统模型的方便性，在飞行器动态分析中更习惯与采用弹簧质量模型。由两种模型均取一个运动质们并入全系统模型的方便性，在飞行器动态分析中更习惯与采用弹簧质量模型。由两种模型均取一个运动质量点量点 m1 可知，模型中值考虑晃动的一阶模态。可知，模型中值考虑晃动的一阶模态。本讲稿第四十九页，共五十二页 并入全系统动态分析的两种方式并入全系统动态分析的两种方式二、二、贮箱内液体晃动对固有特性的影响贮箱内液体晃动对固有特性的

28、影响液体晃动的动态特征引入全系统的方式有以下两种。液体晃动的动态特征引入全系统的方式有以下两种。(1)将它从基本模型中分离出来，单独拟定模型，导出其刚度、阻尼、惯性的等特性参量，将它从基本模型中分离出来，单独拟定模型，导出其刚度、阻尼、惯性的等特性参量，然后以独立的广义坐标耦合到总系统中去构成总运动方程，联立求解。这样求解的优点是，可然后以独立的广义坐标耦合到总系统中去构成总运动方程，联立求解。这样求解的优点是，可以形象地看到晃动的液体在整个系统动态特性中的地位与作用。以形象地看到晃动的液体在整个系统动态特性中的地位与作用。(2)将晃动液体的质量特性、刚度特性直接并入全将晃动液体的质量特性、刚

29、度特性直接并入全系统模型对应位置上，如图系统模型对应位置上，如图8-12所示。对于流体晃动所示。对于流体晃动的影响，则在总系统中增加一个分支系统反映。这种的影响，则在总系统中增加一个分支系统反映。这种考虑方式简单，但较为粗糙。有点飞行器在动态分析考虑方式简单，但较为粗糙。有点飞行器在动态分析时，初始阶段采用此法处理，在最后分析阶段则改用时，初始阶段采用此法处理，在最后分析阶段则改用第一种方式处理。第一种方式处理。“土星土星V”发射器就经历了这一过程。发射器就经历了这一过程。本讲稿第五十页，共五十二页 并入全系统动态分析的两种方式并入全系统动态分析的两种方式二、二、贮箱内液体晃动对固有特性的影响

30、贮箱内液体晃动对固有特性的影响 关于晃动液体对全系统动态固有特性的影响，其处理方法虽然基本上可以满足工程需求，但是，关于晃动液体对全系统动态固有特性的影响，其处理方法虽然基本上可以满足工程需求，但是，由于液体的晃动而使结构固有特性频带变宽，增加了发生耦合振动的可能性；另外，它提供了一个由于液体的晃动而使结构固有特性频带变宽，增加了发生耦合振动的可能性；另外，它提供了一个激励源，是导致系统动态失稳的渠道之一；同时，由于燃料量在整个飞行器中不断变化，从而使全激励源，是导致系统动态失稳的渠道之一；同时，由于燃料量在整个飞行器中不断变化，从而使全系统固有频率与振型都成为时间的变化量；这些都给动力设计带来了附加困难。为此，工程上往往系统固有频率与振型都成为时间的变化量；这些都给动力设计带来了附加困难。为此，工程上往往采用结构措施（如设计阻尼挡板，采用集束式贮箱）或系统化措施（燃料按程序转移，燃烧次序化采用结构措施（如设计阻尼挡板，采用集束式贮箱）或系统化措施（燃料按程序转移，燃烧次序化等），尽量降低液体晃动所产生的动力影响。等），尽量降低液体晃动所产生的动力影响。本讲稿第五十一页，共五十二页作业：本讲稿第五十二页，共五十二页