第四章 连续体的振动优秀PPT.ppt

第四章 连续体的振动第一页，本课件共有52页本章只讨论理想弹性体的振动理想弹塑性体满足以下假设条件各向同性；均质分布；服从虎克定律4.1弦的振动讨论两端受到张力 拉紧的弦，弦上还受到横向干扰力的作用第四章 连续体的振动第二页，本课件共有52页设弦的密度为 （质量/单位体积）假设小变形，弦力不随挠度变化。则弦上的任意一点的位移y应为位置x与时间t的函数，即沿方向作用在微小区间的外力之和为第三页，本课件共有52页根据牛顿第二定律，弦的单元微段ds沿 方向的运动微分方程为：代入得：设代入得：C为波沿长度方向的传播速度第四页，本课件共有52页如无干扰力作用时，称为波动方程 弹性体系统作某阶主振动时，其各点也应当作同样的频率及相位运动，各点也应当同时通过静平衡位置和到达最大偏离位置，即系统具有一定的与时间无关的振型为振型函数第五页，本课件共有52页得故 4个待定常数，可由弦的边界条件及振动的两个初始条件来确定。由于两端固定，故有第六页，本课件共有52页得 则得第七页，本课件共有52页强迫振动 对于长为的两端固定受分布力作用下的弦的强迫振动，其运动微分方程为：振型函数 令则有设其解为 代入方程第八页，本课件共有52页将上式两边同乘以并从得：整理后得到：其通解为：第九页，本课件共有52页4.2杆的纵向振动设杆的横截面在振动时仍保持为平面并作整体运动。略去杆纵设杆的横截面在振动时仍保持为平面并作整体运动。略去杆纵向伸缩而引起的横截面变形。向伸缩而引起的横截面变形。取杆的纵向作为取杆的纵向作为 x 轴，各个截面的纵向位移表示为轴，各个截面的纵向位移表示为u(x,t)，如如图图所示。杆的微元杆的微元dx在自由振动中的受力图也在图在自由振动中的受力图也在图5-3中给出。中给出。设杆单位体积的质量为设杆单位体积的质量为，杆长为，杆长为 l，截面积为，截面积为A，材料的弹性模量为材料的弹性模量为E。再设任一。再设任一x截面处，纵向应变为截面处，纵向应变为(x)，纵向拉力表示为纵向拉力表示为N(x)第十页，本课件共有52页4.2杆的纵向振动第十一页，本课件共有52页由牛顿第二定律：当杆为均质，等截面时设第十二页，本课件共有52页例-1：两端固定的等值杆代入得则有 第十三页，本课件共有52页各阶主振型 各阶主振动：第十四页，本课件共有52页左端：右端：杆受到弹簧力代入由上式，根据不同的值，可解出不同的固有频率例5.2 如图所示，一长为 的等截面均质直杆，左端固定，右端联结一刚度为 的弹簧，试求其纵向自由振动的运动方程。第十五页，本课件共有52页 相当于自由端 主振型 即相当于固定端，其频率方程为第十六页，本课件共有52页4.3圆轴的扭转振动 取取圆圆轴轴的的轴轴心心线线作作为为x 轴轴，图图示轴轴任任一一 x截截面面处处的的转转角角表表示示为为(x,t)。设设轴轴长长为为l ，单单位位体体积积的的质质量量为为，圆圆截截面面对对其其中中心心的的极极惯惯量量矩矩为为Jp，材材料料的的剪剪切切弹弹性性模模量量为为 G 。轴轴的的扭扭转转应应变变为为 ，作作用用于于微微元元dx 两截面上的扭矩分别为两截面上的扭矩分别为 ，及，及 。假设轴的横截面在扭转振动中仍保持为平面作整体转动。假设轴的横截面在扭转振动中仍保持为平面作整体转动。第十七页，本课件共有52页4.3圆轴的扭转振动杆的扭转振动，抗扭刚度 截面处的扭转角由材料力学知识微元段扭矩的增量为极惯性矩：圆截面极惯性矩：密度的转动惯量）第十八页，本课件共有52页得：设 则有：:表示剪切波在杆内的传播速度它的解为式中四个待定常数 由系统的边界条件和初始条件确定。一般解为：第十九页，本课件共有52页解：边界条件：对于圆盘的运动微分方程：例-3：上端固定，下端装有转动惯量为的圆盘，圆轴的极惯性矩为JP，其剪切弹性模量为G，试分析其扭转振动的频率方程。JPIPJJ第二十页，本课件共有52页由于边界条件得：B=0简化得：得：频率方程代入得：的物理意义为轴的转动惯量与圆盘转动惯量之比。的物理意义为轴的转动惯量与圆盘转动惯量之比。设其中(a）第二十一页，本课件共有52页如近似地取如近似地取，则（，则（a a）式化简为）式化简为下表给出对应于各个不同的下表给出对应于各个不同的 值时，基本值时，基本特征值特征值1的值。的值。0.010.010.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）上式也就是略去轴的质量后所得单自由度系统的上式也就是略去轴的质量后所得单自由度系统的固有频率公式固有频率公式。第二十二页，本课件共有52页这时有这时有（c c）上式也就是将轴的转动惯量的三分之一加到圆盘上后所得单自上式也就是将轴的转动惯量的三分之一加到圆盘上后所得单自由度扭振系统的固有频率公式。只要轴的转动惯量不大于圆盘的转由度扭振系统的固有频率公式。只要轴的转动惯量不大于圆盘的转动惯量，那么计算基频的近似式（动惯量，那么计算基频的近似式（c c）在实用上已足够准确了。）在实用上已足够准确了。综上所述，弦的横振、杆的纵振与轴的扭振都导致同一形式的波动综上所述，弦的横振、杆的纵振与轴的扭振都导致同一形式的波动方程。它们的运动具有共同的规律，如表方程。它们的运动具有共同的规律，如表4-14-1。进一步的近似可取进一步的近似可取第二十三页，本课件共有52页弦的横振弦的横振杆的纵振杆的纵振轴的扭振轴的扭振物物理理参参数数弦的张力弦的张力弦的体密度弦的体密度弹性模量弹性模量截面积截面积密度密度剪切弹性模量剪切弹性模量截面极惯性矩截面极惯性矩密度密度截面的截面的位移位移横向位移横向位移纵向位移纵向位移转转角角单位长度单位长度的质量或的质量或转动惯量转动惯量截面处力截面处力（或扭矩）（或扭矩）表表4-14-1弦的横振、杆的纵振与轴的扭振对比表弦的横振、杆的纵振与轴的扭振对比表第二十四页，本课件共有52页弦的横振弦的横振杆的纵振杆的纵振轴的扭振轴的扭振波速波速运动微运动微分方程分方程通解通解边界条件边界条件两端固定两端固定两端自由两端自由一端固定一端固定一端自由一端自由固有频率固有频率振型函数振型函数第二十五页，本课件共有52页作业：p228:4-1;4-5第二十六页，本课件共有52页4.4梁的弯曲振动梁的振动都在同一平面内，即在xoy平面之内满足平截面假设梁上各点的运动只需用轴线的横向位移来表示建立运动微分方程：按牛顿第二定律，微段沿方向的运动方程为：.（1）第二十七页，本课件共有52页由材料力学可知：代入（1）式中得：则有：为梁弯曲振动的偏微分方程梁自由振动的解分离变量：代入得：第二十八页，本课件共有52页设 则有，代入设 代入 的通解为：第二十九页，本课件共有52页则：第三十页，本课件共有52页求固有频率和主振型例-5因此有：第三十一页，本课件共有52页代入B=D=0则有：第三十二页，本课件共有52页例-6：悬臂梁的情况 参见书上的例题。第三十三页，本课件共有52页4.5 主振型的正交性对于等截面的均匀梁，其截面面积，线密度刚度分布都是常数时，有对于一般情况，各阶主振型正交，表明各阶主振动之间没有耦合，各主振动之间是独立无关的，它们之间没有能量交换，都可以采用解耦分析法（振型分析法），简化为类似单自由度系统那样的微分方程来求解。作业：作业：p228:4-8；4-11第三十四页，本课件共有52页4.6薄板的横向弯曲振动基本假设：1直法线假设；2为主应力；3厚度变化不计；4板弯曲的中面不产生变形，即中面为中性面。一、建模第三十五页，本课件共有52页矩形薄板的振动微分方程惯性载荷：二、求解-分离变量法设：得：第三十六页，本课件共有52页三、边界条件1.简支边：2固定边3自由边第三十七页，本课件共有52页四、四边简支板自由振动的解设：固有频率第三十八页，本课件共有52页结构动力学分析工程应用介绍结构动力学分析工程应用介绍 随着形形色色飞行器的出现，动力固有特性分析也出现许多特殊的随着形形色色飞行器的出现，动力固有特性分析也出现许多特殊的问题。例如，火箭、导弹虽然基本上近似一个一维梁系统，但它们往往问题。例如，火箭、导弹虽然基本上近似一个一维梁系统，但它们往往存在着级间连接与数目较多的分离面存在着级间连接与数目较多的分离面(接头接头),这使得结构的刚度分布、阻这使得结构的刚度分布、阻尼分布发生变化，给固有特性的精确分析带来困难；大型液体燃料火箭尼分布发生变化，给固有特性的精确分析带来困难；大型液体燃料火箭广泛存在着液体晃动问题，它对固有特性有不可忽视的影响；人造卫星广泛存在着液体晃动问题，它对固有特性有不可忽视的影响；人造卫星经常采用自旋稳定方案，并安装有柔度很大的附件，如太阳能帆板、天经常采用自旋稳定方案，并安装有柔度很大的附件，如太阳能帆板、天线等，构成具有转动的刚体及柔性体的复合结构，使它的固有特性分析线等，构成具有转动的刚体及柔性体的复合结构，使它的固有特性分析更加复杂化。更加复杂化。本章对这些问题做一概括介绍。本章对这些问题做一概括介绍。第三十九页，本课件共有52页一、一、飞行器分离面（接头）对固有特性的影响飞行器分离面（接头）对固有特性的影响 大部分火箭、导弹类型的飞行器都具有较多的接头，它们使全弹的刚度分布大部分火箭、导弹类型的飞行器都具有较多的接头，它们使全弹的刚度分布发生局部扰动，如图所示。发生局部扰动，如图所示。第四十页，本课件共有52页一、一、飞行器分离面（接头）对固有特性的影响飞行器分离面（接头）对固有特性的影响导导弹弹名名称称接接头头数数刚刚度度损损失失K注注“响尾蛇”413近程攻击导弹SRAM631中程攻击导弹MR731远程标准导弹ER733“不死鸟”1049接头使刚度损失可达接头使刚度损失可达(3040)%,如表如表81所示所示表表 8 81 1第四十一页，本课件共有52页一、一、飞行器分离面（接头）对固有特性的影响飞行器分离面（接头）对固有特性的影响 对弹体固有特性的影响对弹体固有特性的影响1.1.使全弹的固有频率下降使全弹的固有频率下降 由于接头削弱了附近弹体的刚度，同时接头处往往存在空隙，因而它们都使全由于接头削弱了附近弹体的刚度，同时接头处往往存在空隙，因而它们都使全弹的固有频率降低。影响程度与接头的数量、类型、位置相关。表弹的固有频率降低。影响程度与接头的数量、类型、位置相关。表8-2中列举了一中列举了一些导弹的一阶频率由于接头而引起的下降情况。些导弹的一阶频率由于接头而引起的下降情况。导导弹弹名名称称接接头头数数一阶固有频率下降一阶固有频率下降“响尾蛇”47近程攻击导弹SRAM617中程攻击导弹MR717远程标准导弹ER718“不死鸟”1033表表 8-2 8-2第四十二页，本课件共有52页一、一、飞行器分离面（接头）对固有特性的影响飞行器分离面（接头）对固有特性的影响 对弹体固有特性的影响对弹体固有特性的影响2.2.使全弹的振动发生畸形使全弹的振动发生畸形 接头的存在，改变了刚度分布，必然使振型形状、节点位置发生变化。在控接头的存在，改变了刚度分布，必然使振型形状、节点位置发生变化。在控制系统的设计中节点位置是个重要参数，所以，为了精确确定振型，必须考虑接制系统的设计中节点位置是个重要参数，所以，为了精确确定振型，必须考虑接头的响应。图头的响应。图8-2所示为某弹前三阶振型受接头影响的变换情况。所示为某弹前三阶振型受接头影响的变换情况。图中图中x为弹体轴向坐标，坐标原点设在弹体头部理论顶点，为弹体轴向坐标，坐标原点设在弹体头部理论顶点，为振型幅值。为振型幅值。应当注意，同样的接头，所处的位置不同影响也不同。一般来讲，在导弹弹应当注意，同样的接头，所处的位置不同影响也不同。一般来讲，在导弹弹体中部的接头影响更加突出。体中部的接头影响更加突出。第四十三页，本课件共有52页一、一、飞行器分离面（接头）对固有特性的影响飞行器分离面（接头）对固有特性的影响 分析方法分析方法 精确地用纯分析方法考虑接头进行固有特性计算是困难的，一般都采用实验与分析精确地用纯分析方法考虑接头进行固有特性计算是困难的，一般都采用实验与分析结合的方法。全弹的计算模型可以选用一维梁式模型或三维壳体模型，接头则可分为处结合的方法。全弹的计算模型可以选用一维梁式模型或三维壳体模型，接头则可分为处理为集中弯曲弹簧或沿分离面周线分布的弹性组件。处理这类模型的理为集中弯曲弹簧或沿分离面周线分布的弹性组件。处理这类模型的关键是这些弹性件关键是这些弹性件柔度的确定柔度的确定。对于集中弯度弹簧，其柔度确定的方法有以下几种。对于集中弯度弹簧，其柔度确定的方法有以下几种。1.1.类比法类比法 参照已有的导弹接头的柔度数据，用相似类比方法，推测所设计的接头的参照已有的导弹接头的柔度数据，用相似类比方法，推测所设计的接头的柔度。柔度。2.2.经验公式法经验公式法Alley和和Leadbetfer根据大量实验的统计，归纳出以下经验公式根据大量实验的统计，归纳出以下经验公式式中式中 第第i个接头的弯曲柔度；个接头的弯曲柔度；接头处弹身直径（英寸）接头处弹身直径（英寸）柔度系数，根据接头的不同类型已制成表格供查（表柔度系数，根据接头的不同类型已制成表格供查（表8-3）第四十四页，本课件共有52页一、一、飞行器分离面（接头）对固有特性的影响飞行器分离面（接头）对固有特性的影响接头分类Ai额定值范围优10-10310-10良10-9310-10310-9中10-8310-9310-8差10-7310-8310-7表表 9-3此经验公式是以英制给出的，使用时应予注意。此经验公式是以英制给出的，使用时应予注意。3.3.实验测定法实验测定法4.4.其它方法其它方法 当已具有实体结构时，可通过静力实验或动力实验得到各个接头实际当已具有实体结构时，可通过静力实验或动力实验得到各个接头实际柔度。不过，实验中应注意消除弹性弯曲的影响。柔度。不过，实验中应注意消除弹性弯曲的影响。也可采用有限元或最佳拟合的分析方法来确定接头的柔度，不过方法复杂而精也可采用有限元或最佳拟合的分析方法来确定接头的柔度，不过方法复杂而精度并不理想。度并不理想。分析方法分析方法第四十五页，本课件共有52页一、一、飞行器分离面（接头）对固有特性的影响飞行器分离面（接头）对固有特性的影响 总之，随着飞行器对固有频率、振型、振型斜率的数量与精度要求日益提高，总之，随着飞行器对固有频率、振型、振型斜率的数量与精度要求日益提高，接头产生的影响必须予以考虑。由于接头类型较多，单纯的分析方法尚未完善，接头产生的影响必须予以考虑。由于接头类型较多，单纯的分析方法尚未完善，目前主要依靠实验来确定特性参数目前主要依靠实验来确定特性参数柔度柔度。将所得柔度参量代入系统模型，。将所得柔度参量代入系统模型，即可计得较为精确得全弹固有特性。即可计得较为精确得全弹固有特性。分析方法分析方法第四十六页，本课件共有52页 问题的特点问题的特点二、二、贮箱内液体晃动对固有特性的影响贮箱内液体晃动对固有特性的影响晃动频率及晃动激烈程度均与下列因素有关：晃动频率及晃动激烈程度均与下列因素有关：(1)容器形状；容器形状；(2)推进器性质；推进器性质；(3)阻尼隔板设置情况；阻尼隔板设置情况；(4)推进剂液面高度；推进剂液面高度；(5)加速度场的情况。加速度场的情况。在飞行过程中，随着燃料不断燃烧，推进剂液面情况不断发生变化，从而对系统固有特性的影响在飞行过程中，随着燃料不断燃烧，推进剂液面情况不断发生变化，从而对系统固有特性的影响也随之变化。这是本问题的主要特点。一般说来，为了掌握整个飞行过程的固有特性情况，就要分析也随之变化。这是本问题的主要特点。一般说来，为了掌握整个飞行过程的固有特性情况，就要分析各个不同的推进剂燃烧阶段的频率与振型。各个不同的推进剂燃烧阶段的频率与振型。另外一个特点是，一般只需考虑低阶情况，特别是一阶情况。因为经研究指出，对于圆柱另外一个特点是，一般只需考虑低阶情况，特别是一阶情况。因为经研究指出，对于圆柱壳体，二阶晃动质量仅为一阶晃动质量的壳体，二阶晃动质量仅为一阶晃动质量的3，而且在高阶情况下液体内部将产生紊乱的扰动，而且在高阶情况下液体内部将产生紊乱的扰动，使阻尼激增，故二阶以上可不予考虑。使阻尼激增，故二阶以上可不予考虑。第四十七页，本课件共有52页 推进剂晃动频率的确定推进剂晃动频率的确定二、二、贮箱内液体晃动对固有特性的影响贮箱内液体晃动对固有特性的影响 工程上常采用当量变换的方法，在对壳体壁作用的力与力矩相等、频率相当的条件下，将工程上常采用当量变换的方法，在对壳体壁作用的力与力矩相等、频率相当的条件下，将液体晃动模型等价代换为机械力学模型。一旦建立了当量机械模型，对于各种飞行器的液体晃液体晃动模型等价代换为机械力学模型。一旦建立了当量机械模型，对于各种飞行器的液体晃动问题，可根据它们的液体参量、飞行状态参量、飞行器参量很容易地确定出当量机械模型参动问题，可根据它们的液体参量、飞行状态参量、飞行器参量很容易地确定出当量机械模型参量，从而确定出晃动频率。一般采用的液体侧向晃动当量机械模型有两种。量，从而确定出晃动频率。一般采用的液体侧向晃动当量机械模型有两种。1.1.弹簧质量模型弹簧质量模型 根据壳体半径及液面高度，可按流体动力根据壳体半径及液面高度，可按流体动力学分析导得的当量公式确定贮箱中液体固定质学分析导得的当量公式确定贮箱中液体固定质量量m0、晃动质量、晃动质量m1。晃动质量的运动受到弹簧与阻。晃动质量的运动受到弹簧与阻尼器的约束，其模型如图尼器的约束，其模型如图8-10所示。所示。模型中阻尼系数模型中阻尼系数C、弹簧刚度、弹簧刚度K1都按流体动力都按流体动力学导得的当量公式确定。当阻尼较小时，阻尼的学导得的当量公式确定。当阻尼较小时，阻尼的作用可以忽略。作用可以忽略。第四十八页，本课件共有52页 推进剂晃动频率的确定推进剂晃动频率的确定二、二、贮箱内液体晃动对固有特性的影响贮箱内液体晃动对固有特性的影响2.2.自由摆模型自由摆模型 从晃动的物理现象来看，将它等价为一个当量摆是无可非议的。此模型如图从晃动的物理现象来看，将它等价为一个当量摆是无可非议的。此模型如图8-11所示。其中所示。其中 等参量也是由流体动力学分析得到的当量公式来确定。等参量也是由流体动力学分析得到的当量公式来确定。进行这种模型代换后，液体晃动问题的处理就与结构系统的其它部件的处理方式完全相同。考虑到将它进行这种模型代换后，液体晃动问题的处理就与结构系统的其它部件的处理方式完全相同。考虑到将它们并入全系统模型的方便性，在飞行器动态分析中更习惯与采用弹簧质量模型。由两种模型均取一个运动们并入全系统模型的方便性，在飞行器动态分析中更习惯与采用弹簧质量模型。由两种模型均取一个运动质量点质量点 m1 可知，模型中值考虑晃动的一阶模态。可知，模型中值考虑晃动的一阶模态。第四十九页，本课件共有52页 并入全系统动态分析的两种方式并入全系统动态分析的两种方式二、二、贮箱内液体晃动对固有特性的影响贮箱内液体晃动对固有特性的影响液体晃动的动态特征引入全系统的方式有以下两种。液体晃动的动态特征引入全系统的方式有以下两种。(1)将它从基本模型中分离出来，单独拟定模型，导出其刚度、阻尼、惯性的等特性参量，然后将它从基本模型中分离出来，单独拟定模型，导出其刚度、阻尼、惯性的等特性参量，然后以独立的广义坐标耦合到总系统中去构成总运动方程，联立求解。这样求解的优点是，可以形象地看以独立的广义坐标耦合到总系统中去构成总运动方程，联立求解。这样求解的优点是，可以形象地看到晃动的液体在整个系统动态特性中的地位与作用。到晃动的液体在整个系统动态特性中的地位与作用。(2)将晃动液体的质量特性、刚度特性直接并入将晃动液体的质量特性、刚度特性直接并入全系统模型对应位置上，如图全系统模型对应位置上，如图8-12所示。对于流体晃所示。对于流体晃动的影响，则在总系统中增加一个分支系统反映。这动的影响，则在总系统中增加一个分支系统反映。这种考虑方式简单，但较为粗糙。有点飞行器在动态分种考虑方式简单，但较为粗糙。有点飞行器在动态分析时，初始阶段采用此法处理，在最后分析阶段则改析时，初始阶段采用此法处理，在最后分析阶段则改用第一种方式处理。用第一种方式处理。“土星土星V”发射器就经历了这一过发射器就经历了这一过程。程。第五十页，本课件共有52页 并入全系统动态分析的两种方式并入全系统动态分析的两种方式二、二、贮箱内液体晃动对固有特性的影响贮箱内液体晃动对固有特性的影响 关于晃动液体对全系统动态固有特性的影响，其处理方法虽然基本上可以满足工程关于晃动液体对全系统动态固有特性的影响，其处理方法虽然基本上可以满足工程需求，但是，由于液体的晃动而使结构固有特性频带变宽，增加了发生耦合振动的可能需求，但是，由于液体的晃动而使结构固有特性频带变宽，增加了发生耦合振动的可能性；另外，它提供了一个激励源，是导致系统动态失稳的渠道之一；同时，由于燃料量性；另外，它提供了一个激励源，是导致系统动态失稳的渠道之一；同时，由于燃料量在整个飞行器中不断变化，从而使全系统固有频率与振型都成为时间的变化量；这些都在整个飞行器中不断变化，从而使全系统固有频率与振型都成为时间的变化量；这些都给动力设计带来了附加困难。为此，工程上往往采用结构措施（如设计阻尼挡板，采用给动力设计带来了附加困难。为此，工程上往往采用结构措施（如设计阻尼挡板，采用集束式贮箱）或系统化措施（燃料按程序转移，燃烧次序化等），尽量降低液体晃动所集束式贮箱）或系统化措施（燃料按程序转移，燃烧次序化等），尽量降低液体晃动所产生的动力影响。产生的动力影响。第五十一页，本课件共有52页作业：第五十二页，本课件共有52页