基于蒙特卡洛法的油气悬架路面不自由度时间函数仿真.docxVIP

基于蒙特卡洛法的油气悬架路面不自由度时间函数仿真 自20世纪60年代末以来，karkopp开发的气挤出器以来，车辆的良好特性一直受到许多技术研究人员的关注。油气悬架的结构最先应用在德国和日本的重型车辆上,以后逐步推广应用到军用特种车辆、工程机械等车辆上。 悬架直接影响车辆的平顺性和安全性,传统车辆悬架刚度和阻尼均不可变,仅可适用于某些路面。而油气悬架具有良好的非线性弹性特性和阻尼特性,较好地满足了不同路面的使用要求,提高了车辆的平均行驶车速和坏路的通过性能。 车辆油气悬架系统是典型非线性系统,传统的基于叠加规律的频域分析方法不能适用,应采用时域分析的方法求得振动响应的解。路面对车辆的激励通常是不确定和满足统计规律,求解困难。应用MATLAB/SIMULINK进行仿真计算可以获得较好的效果。 1 油液的压缩作用 单气室油气悬架一般由悬架缸、主活塞、蓄能器和活塞杆组成。某单气室油气悬架结构如图1所示,主活塞与悬架缸构成了主油室Ⅰ,中空的活塞杆内腔为副油室Ⅱ,把主油室和蓄能器以浮动活塞分割,蓄能器的气囊和浮动活塞组成气室Ⅳ,各腔的压力、体积、面积如图1所示。 当悬架处在压缩行程(活塞杆向上运动)时,I腔的油液受到压缩向两个方向运动,一部分油液经过单向阀2和阻尼孔1流入Ⅱ腔;另一部分油液进入蓄能器,使气室容积减小,氮气压力升高。在这一过程中由于单向阀2和阻尼孔1同时使中心腔与Ⅱ腔连通,油液流经单向阀2和阻尼孔1的流速较低,产生的油液阻尼力比较小,因此主要由蓄能器内的气体受到压缩产生弹性作用来抑制活塞杆的向上运动。 伸张行程时,活塞杆向下运动,I腔的体积增大、压力减小,Ⅱ腔体积减小、压力增大,Ⅱ腔和气室Ⅳ中的部分油液流向I腔,此时因油液单向阀关闭,产生的阻尼力较大,用于衰减振动。 分析物理模型时作如下简化:油液的质量与悬架油缸组件的质量均忽略不计;各构件没有弹性变形;各密封环节工作可靠、系统没有内泄漏和外泄漏;悬架油缸润滑良好、忽略库仑摩擦的影响;油液是连续不可压缩的、不考虑油液流动过程中气穴现象的影响、不考虑管路的沿程损失和局部损失。可得 F=p1A1-p2A2(1) 静平衡位置时,各腔室的压力相同,设为ps,即 p1=p2=p4=ps(2) 在活塞被压缩和伸张时,主油室内油液体积变化等于蓄能器内油液体积变化,即A1x=Q14(3) 式中,Q14为主油室流入蓄能器的流量,流入为正,流出为负。 由孔口淹没出流理论得,主油室和副油室的压力差为 p1?p2=ρ2(Q12CdA)2sgn(x˙)(4)p1-p2=ρ2(Q12CdA)2sgn(x˙)(4) p2=p4(5) 式中,A为阻尼孔和单向阀的面积之和,在压缩行程为阻尼孔和单向阀面积之和,在伸张行程为阻尼孔的面积;Cd为节流口流量系数,取0.62。 蓄能器内气体体积变化满足 p4Vr44r=psVr4s4sr(6) 式中,V4为蓄能器气囊内气体的体积;V4s为静平衡位置蓄能器气囊内气体体积。 当汽车振动缓慢时,气体的状态变化可以近似看作等温过程,r=1.0。当汽车振动剧烈时,气体的状态变化近似于绝热过程,r=1.4,一般情况r取1.3。 式(1)～式(6)即为单气室油气分离式油气悬架的非线性数学模型。 2 双质量振动系统 当汽车的悬挂质量分配系数接近1时,可以认为车辆前、后轴上方的簧载质量在垂直方向独立振动,因此可以把汽车简化为二自由度的双质量振动系统。 2.1 油气悬架振动方程 设车身质量为m2,车桥质量为m1,轮胎刚度为kt,悬架刚度为k,悬架阻尼系数为c,x0,x1,x2分别为路面、车桥、车身的垂向位移,方向朝上,可以建立图2所示模型。 可得振动方程如下 m2x??2+c(x˙2?x˙1)+k(x2?x1)=0(7)m2x??2+c(x˙2-x˙1)+k(x2-x1)=0(7) m1x??1?c(x˙2?x˙1)?k(x2?x1)+kt(x1?x0)=0(8)m1x??1-c(x˙2-x˙1)-k(x2-x1)+kt(x1-x0)=0(8) 式(7)～式(8)中油气悬架的刚度k和阻尼c是非线性的,可以根据公式(1)～式(6)求解。 x0是路面垂直位移的数值,沿行驶方向变化,一般可以表示为路程的函数,考虑车速时,一般表示为时间的函数。 2.2 路面不清晰度的时域模型构造 非线性悬架系统的引入,使研究需在时域中进行,所需的路面不平度激励也应是时域信号。时域路面激励信号的获得方法主要有两种,即通过试验测量和使用时域模型仿真。路面不平度的常用表示方法有幂函数和有理函数两种,均满足理论计算要求。 路面不平度常用幂函数表示为 Gq(n)=Gq(n0)(nn0)?w(9)Gq(n)=Gq(n0)(nn0)-w(9) 式中,Gq(n0)为路面不平度系数,根据路面等级的不同数值不同;w为频率指数,