单辆双轴汽车作用下车桥耦合系统的振动理论分析及其数值模拟.docVIP

单辆双轴汽车作用下车桥耦合系统的振动理论分析及其数值模拟 　　摘要：通过对汽车振动模型简化成平面两自由度振动体系的基础上，对车桥耦合作用的机理进行了力学分析，并结合MATLAB强大的数值计算功能对车桥耦合振动问题进行数值模拟。 　　关键词：耦合 平整度 MATLAB语言 　　中图分类号：U46文献标识码：A文章编号：1006-4311（2011）14-0044-02 　　Vibration Theory Analysis of Bridge Coupling System and Its Numerical Simulation under the Action of a Single Vehicle Shaft Vehicle 　　Yin Yushi；Wu Bin；Yin Jingkui 　　（①Liaoning Provincial College of Communication，Shenyang 110012，China；②Benxi Hengda Real Estate Development Co.，Ltd.，Benxi 117200，China） 　　Abstract: Through simplifying the vehicle vibration model to planar two-freedom vibration systems, the mechanical analysis of the role of bridge coupling mechanism are made, and the numerical simulation of bridge vibration problems with the powerful numerical calculation function of MATLAB. 　　Key words: coupling；surface roughness；MATLAB 　　0引言 　　移动荷载作用下的车桥耦合振动问题是桥梁系统中十分重要的研究课题，在过去的100年里，对于移动荷载问题，人们已经做了大量的工作[1-3]。但当考虑移动荷载质量影响时，系统的运动方程将是相互耦合的时变微分方程组，只能用数值方法进行求解，给问题的求解与研究带来很大的困难．本文利用MATLAB[4]强大的数值计算功能在对车桥耦合作用理解的基础上，对车桥耦合振动问题进行数值求解。 　　1模型的建立 　　为了便于分析汽车对梁桥产生的竖向弯曲强迫振动，可以仅研究汽车车体和前后桥的竖向振动，把车体假设为刚体，将双轴汽车振动模型简化成平面两自由度振动体系，如图1所示。 　　M1=M■（1） 　　M2=M■（2） 　　式中：M为汽车车体 （包括载重物）的质量；ρ■为Ｍ绕过重心C点水平横桥轴的回转半径，ρ■■=■，l■为汽车的前后轴距，l■=l■+l■。l■，l■分别为汽车重心C至前、后轴的距离。由于大多数汽车的实际e＝0.8~1.2，接近于1。因此，上述分解汽车质量的方法是可行的。 　　为了简化分析，假设车辆在桥梁上行使时始终与桥梁接触，既不考虑车辆的跳离现象，车辆对梁桥的作用力可表示为： 　　F■（t）=M■g+C■Z■■（t）-y■■（t）-r■■■■（t）+k■Z■（t）-y■（t）-r■（t）（3） 　　F■（t）=M■g+C■Z■■（t）-y■■（t）-r■■■■（t）+k■Z■（t）-y■（t）-r■（t）（4） 　　式中：Z■，Z■分别为汽车前、后桥弹簧上质量M的绝对位移；C■和C■是汽车的粘性阻尼系数。 　　由于实际的路面总是不平整的，所以在车桥耦合系统振动分析中有必要考虑路面的不平整度。路面平整度函数r（x）一般为随机过程，但为简化理论分析，方便计算，假设为确定性过程，如图2所示为路面平整度的断面图。 　　与行程有关的不平度函数为： 　　r（x）=h（1-cos■）（5） 　　其中，h为平整度的幅值，l为车辆轨迹的长度。l■为波长。 　　2车桥耦合作用的机理分析及数值模拟 　　2.1 车桥作用力学分析[5]在以常速v行驶时，行程x和时间t之间的关系简化为：x=vt；由此可得到“行程圆频率”Ω和“时间圆频率”之间的关系：ω=VΩ=2π■（6） 　　当波长l为常值时，时间圆频率随车速提高而提高。该车桥模型的受力分析为： 　　前悬挂质量弹簧系统的振动微分方程为： 　　M■Z■■+C■（Z■■-y■■-r■■）+k■（Z■-y■-r■）=0（7） 　　后悬挂质量弹簧系统的振动微分方程为： 　　M■Z■■+C■（Z■■-y■■-r■■）+k■（Z■-y■-r■）=0（8） 　　式中，Z■，Z■分别为汽车前、后桥弹簧上质量M■，M■的绝对