Sensitivity analysis method of vehicle steering system under braking condition
-
摘要: 为了解制动工况下转向系统敏感度对前轮摆振的影响程度, 将车轮、前轴与钢板弹簧作为一个整体进行考虑, 分析了钢板弹簧和转向系统的空间状态变化情况, 提出了前轮摆动角的计算方法, 以转向系统空间状态引起的前轮摆动角为指标评价转向系统敏感度。计算结果表明: 在空载状态下, 当制动力达到35kN时, 前轮摆动角最大值为0.5°。可见, 此车转向系统敏感度较低, 评价结果符合实际, 因此, 计算方法可行。Abstract: In order to determine the effect of the sensitivity of steering system on the shimmy of front wheel under braking condition, wheel, front axle and leaf springs were considered as a whole, the changes of dimensional states of leaf spring and steering system were analyzed, a mathematic method of the shimmy angle of front wheel was proposed, and the sensitivity of steering system was evaluated with the shimmy angle resulted from the change of dimensional state of steering system. Computation result indicates that the maximum of the shimmy angle is 0.5° when braking force is 35 kN. Obviously, the sensitivity of vehicle steering system is lower, the evaluation result accords with its real state, so the method is feasible.
-
Key words:
- automobile engineering /
- steering system /
- sensitivity /
- dimensional state /
- shimmy of front wheel
-
0. 引言
在制动过程中, 车辆载荷发生纵向转移, 前悬架板簧发生弯曲变形和扭转变形, 使得车架与前轴在垂直方向和汽车前进方向的相对位移发生变化, 使转向系统产生运动, 引起前轮的定位参数发生变化, 并产生一定的转向角。如果转向角较大, 就会引起前轮摆振, 汽车转向系统空间敏感度的高低, 对前轮的摆振与整车的操纵稳定性有很大的影响。现阶段, 不仅国外学者已经对转向系统对汽车制动性能和操纵稳定性的影响进行了深入研究[1-2], 而且国内学者也对汽车转向轮摆振的原因进行了大量的理论和试验研究[3-4], 研究表明实际车辆前轮的摆振问题涉及车桥、车轮和转向系统等多个系统[5]。在转向系统中, 除了主销间隙对前轮摆振有影响之外[6], 转向系统空间状态的敏感性是影响制动情况下前轮摆振的一个重要因素, 因此, 对转向系统及其空间状态敏感度进行分析, 有助于解决由转向系统方面的问题引起的前轮摆振, 对提高车辆的制动稳定性和操纵稳定性极为重要。
1. 制动时前悬架上的作用力
根据右手定则建立空间坐标系, 原点固结在前轴中心, x轴沿汽车行驶方向, y轴垂直向上, z轴指向汽车左侧。制动时把车轮、前轴以及板簧作为一个整体进行考虑。在不考虑车身自重作用的情况下, 作用在其上的力有车轮的制动力F、前轴(加板簧) 的惯性力Fa、纵向垂直载荷转移引起的地面的附加力G以及板簧两卷耳处的支撑力Fcx、Fcy、Fdx和Fdy (图 1)。如不考虑板簧之间的摩擦力, 可将整个板簧简化为简支梁。根据受力分析可知, 作用在扭曲中心[7] (板簧主片中心点) 点b处有转矩Tb、力Fb和G
Τb=Fh1-Fah2 (1)Fb=F-Fa (2)G=ΜhBa (3)
式中: h1为扭曲中心离地高度; h2为前轴(加板簧) 质心离地高度; M为车身质量; a为制动加速度; h为车身质心离地高度; B为轴距。另外, 惯性力Fa和车轮制动力F求解式为
Fa=ma (4)F=(G0+G)φ (5)
式中: m为前轴(加板簧) 质量; φ为附着系数; G0为车轮静载荷。
当前、后车轮都有制动力时, 有
a=φg (6)
式中: g为重力加速度。
当只有前轮有制动力时, 有
a=F/Μ (7)
2. 钢板弹簧空间状态分析
由受力分析可知, 制动时板簧在垂直方向由力G作用下, 板簧发生弹性变形, 扭曲中心与车架产生相对位移, 竖向和纵向的位移量分别为
y=GC (8)x=4z3l (9)
式中: C为板簧刚度; l为板簧长度; z为点c与b之间的垂直距离。
另外, 力Fb和转矩Tb会造成板簧卷曲变形而使前轴连同主销、转向节臂都产生一定的转角, 从而使得转向系统的空间状态发生变化。根据受力分析, 可把整个板簧简化为变刚度曲梁, 曲梁圆弧中心为点o, 圆弧半径为R, 扭曲中心为点b, 两卷耳中心分别为点c、d, 见图 2。
当板簧的安装角为γ时, 板簧扭曲中心点b处在纵向和竖向受到的力分别为
Fbx=Fbcos(γ)Fby=-G+Fbsin(γ)
在曲梁上任意一点和钢板弹簧弧线圆心o之间的连线与oc之间的夹角可用θ表示。根据参考文献[8], 当0≤θ≤θcob时, 任意一点处的转矩为
当θcob < θ≤θcod时, 转矩为
点b在纵向、竖向及竖直平面上的3个位移Δx、Δy和Δθ分别为
Δx=∫θcob0Τ1(θ)EΙ(θ)∂Τ1(θ)∂FbxRdθ+∫θcodθcobΤ2(θ)EΙ(θ)∂Τ2(θ)∂FbxRdθ (12)
Δy=∫θcob0Τ1(θ)EΙ(θ)∂Τ1(θ)∂FbyRdθ+∫θcodθcobΤ2(θ)EΙ(θ)∂Τ2(θ)∂FbyRdθ (13) Δθ=∫θcob0Τ1(θ)EΙ(θ)∂Τ1(θ)∂ΤbRdθ+∫θcodθcobΤ2(θ)EΙ(θ)∂Τ2(θ)∂ΤbRdθ (14)
Ι(θ)=bLhθ312 (15)
当0≤θ≤θcob时, 有
hθ=Ηn+(n-1)Ηnθcobθ (16)
当θcob < θ≤θcod时, 有
hθ=Η-(n-1)Ηnθcob(θ-θcob) (17)
式中: E为弹簧弹性模量; I (θ) 为角度区间惯性矩; bL为板簧宽度; hθ为板簧整体厚度的函数; H为板簧总厚度; n为板簧片数。
通过以上分析可知, 制动时板簧的角位移为Δθ。由于制动时车轴会随板簧转动, 扭曲中心点b的角位移Δθ与其引起的前轴转动角β之间的关系为
β=Δθ (18)
3. 转向系统空间状态分析
3.1 转向节臂绕主销转角求解
设定在静止情况下转向机坐标为A (xA, yA, zA), 摇臂与直拉杆铰接点坐标为B (xB, yB, zB), 转向节坐标为C (xC, yC, zC), 转向节臂与主销交点坐标为D (xD, yD, zD)。由钢板弹簧空间状态分析知道, 制动时前轴和车架以及转向机产生相对位移, 转向机坐标变为A1 (xA+x, yA-y-Δy, zA), 摇臂与直拉杆铰接点坐标变为B1 (xB+x, yB-y-Δy, zB), 转向节坐标变为C2 (xC2, yC2, zC2), 转向节臂与主销交点坐标变为D2 (xD2, yD2, zD2)。设定在此状态下转向节臂绕主销转角为θ1。
在制动状态下, 转向节臂(在转向节臂的转向范围内) 只在一个角度上与前轴转动角β满足约束条件, 使得转向节臂与主销交点D1 (xD1, yD1, zD1) 和转向节点C1 (xC1, yC1, zC1) 之间的距离为直拉杆长度L, 而这个角度可以通过以下搜索法求出。
(1) 分别取θ1为∧θ和˜θ‚根据式(20)、(21) 求转向节点坐标C1 (xC1, yC1, zC1) 和D1 (xD1, yD1, zD1), 再分别求不同θ1值对应的点C1与D1之间的距离L1和L2。其中、分别为转向节臂绕主销转动的极限转角, 其区间范围被N等分, 步长为
(2) 分别取θ1为和。
(3) 计算
如果K < 0, K > t, t为精度, 取s为s/N, 回到(2);如果K > 0, 取回到(2);如K < 0, K≤t, 取θ1为, 结束搜索。
3.2 车轮摆动角求解
通过运动分析, 可知转向节点坐标C (xC, yC, zC) 与其绕主销旋转角θ1之后的坐标C2 (xC2, yC2, zC2) 有如下关系[9]
\left(x_{C_{2}}, y_{C_{2}}, z_{C_{2}}, 1\right)=\left(x_{C}, y_{C}, z_{C}, 1\right) \boldsymbol{R} (19) 式中: λ1为主销内倾角; λ2为后倾角; R、T、Rx、Ry与Rz为转动关系矩阵[9]。
前轴有转动角β时, 转向节点坐标C1 (xC1, yC1, zC1) 和转向节臂与主销交点坐标D1 (xD1, yD1, zD1) 有如下关系
\left(x_{C_{1}}, y_{C_{1}}, z_{C_{1}}, 1\right)=\left(x_{C_{2}}, y_{C_{2}}, z_{C_{2}}, 1\right) \boldsymbol{R}_{\beta} (20) \left(x_{D_{1}}, y_{D_{1}}, z_{D_{1}}, 1\right)=\left(x_{D_{2}}, y_{D_{2}}, z_{D_{2}}, 1\right) \boldsymbol{R}_{\beta} (21) 式中: Rβ为转动关系矩阵[9]。
因此, 制动时转向节臂的摆动角α为
4. 算例分析
依照相关的制动试验标准[10], 对某重型货车进行了制动工况下的转向系统敏感度分析。通过实际测量和计算, 得到该车在空载情况下的参数: 安装角γ为0°, 板簧长度l为1 800 mm, 板簧总厚度H为220 mm, 板簧宽度bL为200 mm, 扭曲中心离地高度h1为640 mm, 前轴(加板簧) 质心离地高度h2为300 mm, 板簧弹性刚度C为237 000 N·m, 车身质心离地高度h为1 000 mm, 车身质量M为9 000 kg, 前轴(加板簧) 质量m为150 kg, 直拉杆L为1 030 mm, 主销内倾角λ1为5.75°, 主销外倾角λ2为3.50°。
结合汽车试验标[7], 对该货车在空载、车载质量为10、20 t时的制动工况进行了限定, 在Matlab环境下经过计算, 得到只用前轮制动和前后轮同时制动时制动力与前轮的摆动角的关系, 见图 3, 分析结果如下。
(1) 车载质量相同时, 在同等制动力作用下, 前后轮同时制动比只用前轮制动引起的前轮转向角大, 且车载质量越大, 制动力对前轮转向角的影响越大。
(2) 空载情况下, 制动力为35 kN时前轮的最大转向角能达到0.5° (与试验车的实测数据非常接近); 非空载情况下, 制动力达到一定值时, 板簧基本被压平, 弧高会达到最小, 前轮转向角达到最大值(车载质量为10 t时达到0.25°, 车载质量为20 t时达到0.12°), 但这些角度都小于空载时能达到的角度0.5°。
(3) 制动时前轮能达到的最大转向角为0.5°, 即使制动时制动力有变化, 由转向系统空间状态引起的前轮摆振也只能在小于0.5°的范围内进行。
根据结果分析可以得出, 该车辆转向系统空间状态的敏感度较低, 不会对前轮的摆振产生较大影响。该重型货车生产厂家多次试验结果表明, 以上得到的结论正确。
5. 结语
转向系统的技术状况好坏将影响到汽车的操纵性, 在实际中汽车的转向节臂、转向节弯曲变形等都会造成汽车转向系统空间状态敏感度发生变化, 引起汽车跑偏或前轮摆振, 因此, 分析汽车转向系统空间状态敏感度对研究和解决前轮摆振极为重要。本文通过对制动工况下转向系统的空间状态进行分析, 得到了制动力与车轮摆振角之间的关系式, 通过实车算例分析表明, 应用本文提出的方法得到的结论与厂家提供的结论相符, 为更深层次地研究和解决转向系统与前轮摆振问题提供了基础。
-
[1] Hac A, Bodie M O. I mprovements in vehicle handling through integrated control of chassis systems[J]. International Journal of Vehicle Autonomous Systems, 2002, 11 (1): 83-110. [2] Pilutti T, Ulsoy G, Hrovat D, et al. Vehicle steering intervention through differential braking[J]. Transactions of the ASME, 1998, 120 (9): 314-321. [3] 李胜, 林逸. 汽车转向轮摆振研究[J]. 汽车技术, 2004, 35 (11): 16-19. https://www.cnki.com.cn/Article/CJFDTOTAL-QCJS200411004.htmLi Sheng, Lin Yi. Study on shimmy of automotive steering wheels[J]. Automotive Technology, 2004, 35 (11): 16-19. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-QCJS200411004.htm [4] 贺丽娟, 林逸. 汽车操纵稳定性与前轮摆振的非线性仿真分析[J]. 汽车工程, 2007, 29 (5): 389-392. doi: 10.3321/j.issn:1000-680X.2007.05.007He Li-juan, Lin Yi. Nonlinear simulation on vehicle handling stability and front wheel shimmy[J]. Automotive Engineering, 2007, 29 (5): 389-392. (in Chinese) doi: 10.3321/j.issn:1000-680X.2007.05.007 [5] 韩同群, 邹汉华. EQ2102越野汽车前轮摆振影响因素试验分析及改进措施[J]. 汽车技术, 2007, 38 (1): 34-37. https://www.cnki.com.cn/Article/CJFDTOTAL-QCJS200701009.htmHan Tong-qun, Zou Han-hua. The experi mental study and analysis of influence factors on front wheel shimmy and improving measure of EQ2120 off-road vehicle[J]. Automotive Technology, 2007, 38 (1): 34-37. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-QCJS200701009.htm [6] 顾款. 考虑转向系间隙的汽车前轮摆振系统研究[D]. 合肥: 合肥工业大学, 2007. [7] 郭孔辉. 汽车操纵动力学[M]. 长春: 吉林科学技术出版社, 1991. [8] 王志宏. 后钢板弹簧悬架变形对传动轴最大滑出量的影响[J]. 汽车工程, 2001, 23 (3): 167-171. https://www.cnki.com.cn/Article/CJFDTOTAL-QCGC200103006.htmWang Zhi-hong. The Effect of deformation of rear leaf-spring suspension on the variation of total length of propeller shaft[J]. Automotive Engineering, 2001, 23 (3): 167-171. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-QCGC200103006.htm [9] 孙家广. 计算机图形学[M]. 北京: 清华大学出版社, 1998. [10] GB12676—1999, 汽车制动系统结构、性能和试验方法[S]. -