The Simulation Analysis and Test Research on the Braking Vibration Process of the Tractor-trailer Combinations
-
摘要: 对全挂车辆在平坦路面上的制动过程进行了动力学分析, 确定了全挂车辆制动时振动的激励因素, 应用Wilson-θ逐步积分法并通过MATLAB语言编程进行了全挂车辆制动过程的仿真计算, 仿真分析结果与试验结果相吻合Abstract: The dynamic analysis for the braking process of the tractor trailer combinations running on the even road is performed.The stimulation factors during the braking vibration process are settled.The corresponding simulation analysis using the Wilson θ step by step integral method by means of MATLAB language is developed.The theoretical results are in acceptable agreement with the experimental data.
-
Key words:
- tractor trailer combinations /
- simulation analysis /
- braking tests
-
全挂车辆的制动过程较普通车辆有其特殊性, 因为在具体结构上整个全挂车辆的制动车轴并不在一个刚性车架上, 亦即地面制动力最终作用在2个刚体车架上, 故当主车与挂车的制动分配方案不同时(如主车相对于挂车的制动滞后时间不同), 以及因车辆的振动而导致作用于各制动车轴上的制动力不协调时, 都将加剧全挂车辆制动过程的振动。因此, 研究主、挂车之间的这种振动将对全挂车辆制动过程的改善具有十分重要的意义。
车辆制动时的外力是由地面和空气提供的。在车速较低、空气阻力较小的情况下, 则车辆制动的外力主要来自于地面制动力。车辆在制动过程中的振动主要是由于地面制动力的波动和变化引起的, 因而地面制动力是车辆制动时振动的主要激励因素。
在假定车辆制动系统性能稳定及制动路面较平坦的情况下, 地面制动力主要随制动过程中地面附着系数及作用于车轮上的动反力的变化而变化, 在全挂车辆的整个制动过程中, 这两个方面是交互影响的[1]。而且由于地面附着系数及动反力的变化也并不是一个平稳的随机过程, 因而受此二者影响全挂车辆的制动振动过程将是一个非平稳随机过程。
1. 全挂车辆制动过程的振动动力学分析
1.1 制动过程振动动力学微分方程的建立
假定全挂车辆在无不平度起伏的沥青路面上进行制动, 对牵引车、挂车及挂车悬架分别取隔离体进行动力学分析。考虑到试验用牵引车为无悬挂车辆, 把其简化成如图 1所示的三自由度1/2整车综合振动模型。设ktv1、ktv2、ctv1、ctv2分别为牵引车前后轮的刚度和阻尼, 则地面对车轮的动反力及制动时的制动力可分别表示为
z11=ctv1(-˙z11)+ktv1(-z11) (1)Fz12=ctv2(-˙z12)+ktv2(-z12) (2)Fx11=Fz11f (3)Fx12=(Fz12+G2)μ1(t) (4) Fz11=ctv1(-˙z11)+ktv1(-z11) (1)Fz12=ctv2(-˙z12)+ktv2(-z12) (2)Fx11=Fz11f (3)Fx12=(Fz12+G2)μ1(t) (4)
式中: f为地面滚动阻力系数; G2为平稳状况下分布于牵引车后轴上的静载荷; μ1 (t) 为牵引车后轴车轮处的地面附着系数。
滚动阻力及地面制动力所形成的力矩分别为
11=Fx11R1 (5)Μ12=Fx12R2 (6) Μ11=Fx11R1 (5)Μ12=Fx12R2 (6)
式中: R1、R2分别为牵引车前后轮半径。
对于挂车, 考虑到其为有悬挂车辆, 把其简化成如图 2所示的五自由度1/2整车综合振动模型。设ktv3、ktv4、ctv3、ctv4分别为挂车前后轮的刚度和阻尼, R3、R4分别为挂车前后轮半径, G4为平稳状况下分布于挂车后轴上的静载荷, μ2 (t) 为挂车后轴车轮处的地面附着系数, 同时考虑挂车仅有后轴装有气压式制动器, 则按与牵引车同样的分析方法有
z21=ctv3(-˙z21)+ktv3(-z21) (7)Fz22=ctv4(-˙z22)+ktv4(-z22) (8)Fx21=Fz21f (9)Fx22=(Fz22+G4)μ2(t) (10)Fzs1=cs1(˙z21-˙zs1)+ks1(z21-zs1) (11)Fzs2=cs2(˙z22-˙zs2)+ks2(z22-zs2) (12)Fxs1=Fx21-m¨x21 (13)Fxs2=Fx22-m¨x22 (14)Μ21=Fx21R3 (15)Μ22=Fx22R4 (16) Fz21=ctv3(-˙z21)+ktv3(-z21) (7)Fz22=ctv4(-˙z22)+ktv4(-z22) (8)Fx21=Fz21f (9)Fx22=(Fz22+G4)μ2(t) (10)Fzs1=cs1(˙z21-˙zs1)+ks1(z21-zs1) (11)Fzs2=cs2(˙z22-˙zs2)+ks2(z22-zs2) (12)Fxs1=Fx21-m¨x21 (13)Fxs2=Fx22-m¨x22 (14)Μ21=Fx21R3 (15)Μ22=Fx22R4 (16)
对于牵引车和挂车, 分别应用Newton-Euler矢量法建立其动力学微分方程如下
1¨z1=Fz11+Fz12 (17)Μ1¨x1=Fc0+Fk0+Fx11+Fx21 (18)Ιφ1¨φ1=Fz11l1-Fx11h1-Fz12l2-Fx12h2- (Fc0+Fk)h0-Μ11-Μ22 (19)Μ2¨z2=Fzs1+Fzs2 (20)Μ2¨x2=Fxs1+Fxs2-Fc0-Fk0 (21)Ιφ2¨φ2=Fzs1l3-Fxs1h3-Fzs2l4-Fxs2h4+ (Fc0+Fk)h3-Μ21-Μ22 (22)m¨z21=Fz21-Fzs1 (23)m¨z22=Fz22-Fzs2 (24) Μ1¨z1=Fz11+Fz12 (17)Μ1¨x1=Fc0+Fk0+Fx11+Fx21 (18)Ιφ1¨φ1=Fz11l1-Fx11h1-Fz12l2-Fx12h2- (Fc0+Fk)h0-Μ11-Μ22 (19)Μ2¨z2=Fzs1+Fzs2 (20)Μ2¨x2=Fxs1+Fxs2-Fc0-Fk0 (21)Ιφ2¨φ2=Fzs1l3-Fxs1h3-Fzs2l4-Fxs2h4+ (Fc0+Fk)h3-Μ21-Μ22 (22)m¨z21=Fz21-Fzs1 (23)m¨z22=Fz22-Fzs2 (24)
将式(1) ~ (16) 代入式(17) ~ (24), 并考虑有关牵引车与挂车的几何尺寸关系及其相应坐标变换关系, 则整理后并写成矩阵形式为
Μ0]{⋅⋅X}+[C]t{˙X}+[Κ]t{X}={Q}t (25) [Μ0]{⋅⋅X}+[C]t{˙X}+[Κ]t{X}={Q}t (25)
1.2 制动时全挂车辆系统振动的激励因素分析与描述
车辆在硬路面上制动时, 其地面附着系数μ与车轮的滑移率有关(如图 3所示), 在制动的开始阶段, 地面附着系数随车轮滑移率的增大而增大, 当轮胎的滑移率达到20%左右时, 地面附着系数达到最大值μp, 称其为峰值附着系数。随后随着滑移率的再增加, 地面附着系数有所下降, 滑移率达到100%, 即车轮完全抱死拖滑时的地面附着系数称为滑动附着系数μs。这种变化波动规律决定了地面作用于车轮上的制动力的大小及其变化波动规律, 同时也决定了车辆制动减速度的大小及其变化规律。这可以从车辆制动时汽车制动减速度与制动时间的关系曲线中明显地看出来[2]。
由此可见, 车辆制动时地面与车轮之间的附着系数确实是影响地面制动力的主要因素, 如式(4)、式(10) 所示。随着制动过程中附着系数μ1 (t)、μ2 (t) 的变化, 地面制动力也随之改变, 引起车辆的振动, 并导致制动轮上垂直动反力Fz12、Fz22的变化, 随后地面制动力受路面附着系数与垂直动反力的共同影响而不断波动, 导致车辆更强烈的制动振动。
所以, 在由方程式(25) 所表达的全挂车辆制动振动系统中, 可以把地面附着系数作为振动激励因素。文献[3]根据图 3所示附着系数与滑移率的变化关系, 通过分段回归得到了全挂车辆在沥青路面上制动时附着系数随滑移率的变化曲线, 尽管滑移率也是制动时间的函数, 但在应用上还是不便, 需寻求制动附着系数随制动时间而变化的关系表达式。
对于主车, 假定其单独制动, 考虑式(4) 中地面制动力的线性主部, 即由G2引起的地面制动力部分Fxb1, 则有
xb1=μ1(t)G2=l1μ1(t)l1+l2Μ1g (26) Fxb1=μ1(t)G2=l1μ1(t)l1+l2Μ1g (26)
则主车的制动减速度为
1=l1l1+l2μ1(t)g (27) j1=l1l1+l2μ1(t)g (27)
由此可见, 可用车辆的制动减速度近似地表达地面制动附着系数。设沥青路面的滑动附着系数为μs, 则制动减速度的简化理论模型可描述为(如图 4所示)
1={00≤t < t1t-t1t2-t1jmt1≤t < t2jmt≥t2 (28) j1={00≤t<t1t-t1t2-t1jmt1≤t<t2jmt≥t2 (28)
比较式(27)、式(28), 则主车制动时地面附着系数随时间的变化规律为
对于挂车, 设挂车先于主车制动的时间差为t0, 则同理据式(29) 挂车制动时地面附着系数的数学模型为
2(t)={00≤t < t1-t0t-t1+t0t2-t1μst1-t0≤t < t2-t0μst≥t2-t0 (30) μ2(t)={00≤t<t1-t0t-t1+t0t2-t1μst1-t0≤t<t2-t0μst≥t2-t0 (30)
2. 全挂车辆制动过程的仿真分析与试验研究
2.1 时变方程的时城Wilson-θ逐步积分法
方程式(25) 所示的矩阵形式的动力学方程除质量阵外其阻尼阵、刚度阵及激励力列阵皆变为变系数阵, 取时间步长Δt, 按分段线性化原理, 在每一时间步长Δt内将方程式(25) 视为线性时不变方程, 为此方程式(25) 改为增量形式
Μ0]{Δ⋅⋅X}+[C]t{Δ˙X}+[Κ]t{ΔX}={ΔQ}t (31) [Μ0]{Δ⋅⋅X}+[C]t{Δ˙X}+[Κ]t{ΔX}={ΔQ}t (31)
对式(31) 应用Wilson-θ逐步积分法, 在时域内进行响应分析, 即可求得全挂车辆在制动过程中对应各广义坐标的位移、速度、加速度等时域响应值, 进而可以对全挂车辆的制动过程进行时域响应分析。
2.2 全挂车辆制动过程的时域响应仿真分析与试验
按照上述方法, 取θ=1.5[4], 对式(31) 应用Wilson-θ法进行逐步积分, 并采用MATLAB语言编程进行了全挂车辆制动过程的仿真计算, 分析了全挂车辆在沥青路面上制动时牵引力的变化波动状况及其制动稳定性状况。为了验证上述全挂车辆制动振动过程动力学模型及时域仿真分析的正确性, 在沥青路面上进行了全挂车辆的制动试验, 根据试验测定全挂车辆的制动器起作用时间为0.7 s, 主车制动滞后时间为0.5 s, 试验条件、信号测量、实验方法见参考文献[1]。
图 5、图 6所示为全挂车辆在两种不同制动初速度条件下制动时其牵引架牵引力变化波动状况的时域仿真分析与试验结果对比, 图中牵引力为正时表示牵引架受拉力, 为负时表示牵引架受压力。
(1) 在各种制动初速度水平下制动时, 牵引架牵引力的仿真分析结果与试验结果相吻合, 这说明前述全挂车辆的制动模型及其制动过程动力学分析是正确的和可靠的;
(2) 无论何种初速度下制动, 牵引架牵引力的大致变化趋势都是先增大, 约经过0.5 s的时间后开始减小, 这与全挂车辆各制动车轴的制动顺序相吻合, 因为主车的制动滞后时间约为0.5 s, 当开始制动时, 挂车车轮先受到地面制动力的作用而减速, 而主车由于制动滞后还在按原速度运动;
(3) 随着制动初速度的增加, 牵引架上牵引力的波动幅值范围增大, 如当制动初速度为v=4.12 m/s时, 仿真结果和试验结果分别为16.9359 kN和18.685 kN; 制动初速度为v=6.43 m/s时, 二者分别为26.027 kN和28.1512 kN。由图示结果还可以看出, 随着制动初速度的增加, 牵引架牵引力的负峰值也增大, 由此可见, 全挂车辆的行驶速度越高, 制动时造成的挂车对主车的冲撞力越大, 对全挂车辆的行驶安全性影响也越大。
3. 结语
根据全挂车辆制动过程的特殊性, 建立了全挂车辆制动时的动力学分析模型, 进行了全挂车辆制动过程的仿真分析, 并进行了相应的试验研究, 正确描述了全挂车辆制动时的全车振动状况及牵引架牵引力的变化波动状况, 这为全挂车辆制动过程的进一步改善奠定了基础。
-
[1] 高举成. 挂接车辆的柔性牵引及其平顺性研究[D]. 长春: 吉林大学, 2000. [2] 余志生. 汽车理论(第二版)[M]. 北京: 机械工业出版社, 1989. [3] 左曙光. 路面不平度对拖拉机-挂车机组制动性能的影响——建模与仿真分析[D]. 长春: 吉林大学, 1996. [4] 陈塑寰. 结构振动分析的数值方法[M]. 长春: 吉林科学技术出版社, 1996. -