等效线化方法分析亚音速壁板非线性极限环颤振.docVIP

等效线化方法分析亚音速壁板非线性极限环颤振 　　摘要： 研究了受集中质量与非线性运动约束联合作用下的二维亚音速壁板的极限环颤振问题。采用Galerkin方法将非线性壁板运动方程离散为常微分方程组。分析了集中质量大小及其位置对壁板系统失稳特性的影响；采用等效线化方法研究了系统的分叉特性及极限环颤振稳定性。结果表明：系统会产生颤振失稳，质量块的大小及其位置对颤振临界速度有着重要的影响；系统会经历超临界的Hopf分叉而处于稳定的极限环运动；等效线化方法可在一定范围内较为精确地对极限环稳定性及其幅值进行判定。 　　关键词： 壁板； 亚音速流； 等效线化； Hope分叉； 极限环颤振 　　中图分类号： U270.1+1 文献标志码： A 文章编号： 1004-4523（2015）05-0748-06 　　引 言 　　近些年，随着高速铁路运行速度的不断提升，列车的气动弹性问题也越来越突显，并已经成为高速列车中亟待解决的关键基础问题之一[12]。由于高速列车采用流线型的设计，因此列车车体中存在着大量的蒙皮等壁板结构。这些壁板结构在列车高速运行时会产生明显的振动，如武广客运专线试验中，当列车运行速度达到350 km/h时，列车车身蒙皮和车窗的振动非常显著，并会产生很大的辐射噪声。因此有必要对这类特殊壁板的气动弹性问题进行相关研究。 　　就现有高速列车的运行速度而言，其基本上属于低亚音速范围（马赫数约为0.3）。针对亚音速气流中壁板的气动稳定性问题，文献[2]曾指出两端简支的壁板会在马赫数为0.125时发生颤振失稳，并得到了风洞吹风试验的证实。而文[36]却指出，两端固定支撑的壁板在亚音速气流中并不会发生颤振失稳，而仅会出现发散失稳；而一端固支一端自由的壁板却会出现颤振失稳。除了壁板的气动失稳问题，亚音速壁板结构在失稳后的复杂非线性运动特性也是学者们关注的焦点。现有的研究较为广泛地考虑壁板大变形而产生的几何非线性对系统失稳特性的影响。文[79]均针对该非线性作用下的壁板的稳定性及极限环响应进行了研究。事实上，由于生产、安装过程中的误差，壁板结构常常还会受到其他结构非线性因素的作用。这些非线性因素会约束壁板结构的位移并导致壁板呈现出复杂动力学特性。Li等[1011]考虑支撑松动产生的接触非线性因素，研究了亚音速壁板的极限环颤振及混沌响应。相比于几何非线性而言，针对壁板在位移约束下的研究还比较欠缺。另外，上述研究主要是以理论模型分析为主，均未涉及到实际的风洞模型。事实上，对于实际风洞实验中的壁板而言，不可避免地需要对壁板结构施加某些必要的集中质量，以满足模型设计要求及数据测试要求，例如在壁板关键位置安装有一定质量的传感器及某些必要的实验挂件等。这些额外的重量对壁板的动力学行为，尤其是对非线性极限环运动有何影响，也是需要关注的一个重要问题。 　　因此，本文综合考虑非线性约束及集中质量两个因素的影响，对壁板的稳定性及极限环响应进行分析。文中采用Galerkin方法对非线性亚音速黏弹性壁板的运动方程进行离散；采用等效线化方法研究壁板的非线性运动特性。着重考察集中质量对系统稳定性及极限环运动的影响，并采用数值方法进行积分验证。 　　1 壁板颤振模型 　　2 颤振边界分析 　　首先，考察集中质量及运动约束刚度对系统颤振稳定性的影响。为计算系统的临界颤振速度，将式（10）写作状态空间形式=AX+（X） 　　采用数值方法计算式（11）在零平衡点的雅克比矩阵的特征值，并依据特征值得特性对系统的颤振临界速度进行判定。 　　图2给出了当附加质量块位于端部，即em=1时，不同端部支撑刚度k1对应的临界颤振速度λf与集中质量的变化关系。由图2可知，随着集中质量的增加，颤振临界速度呈现下降趋势，系统的颤振临界速度随着刚度的增加而增加。有趣的是，当集中质量较小时，刚度对颤振临界速度影响较为明显；而当集中质量增加至0.5左右时，不同刚度对应的系统颤振临界速度将趋于相同的值。 　　图3给出了当无量纲集中质量=0.1，不同刚度k1时，系统颤振临界临界速度与集中质量的位置ξm之间的变化关系。图3中的点划线对应的是无集中质量时系统的颤振临界速度。由图3可知：不同刚度时系统的颤振临界速度随位置的变化呈现相似性；系统的颤振临界速度随位置的变化呈现非线性关系。特别值得指出的是：不同刚度对应的颤振边界与点划线的交点有着相同的横坐标，即ξmA=0.558 7，ξmA=0.845 2 　　（12） 当ξmAξmB时，集中质量却起到了相反的作用；而在ξm=ξmA及ξm=ξmB两点处，集中质量对系统的颤振稳定性没有影响。 　　图4给出了=0.1，k1=5.0，λ=49.29时，不同集中质量位置ξm对应的系统线性的响应相图。从图可知：即使在相同的动压下，由于不同集中