第四章 连续体的振动

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