CN105183985A

CN105183985A - 货车前轮摆振系统横向减振器最优阻尼系数的设计方法

Info

Publication number: CN105183985A
Application number: CN201510560015.8A
Authority: CN
Inventors: 周长城; 提艳; 赵雷雷; 于曰伟
Original assignee: Shandong University of Technology
Current assignee: Shandong University of Technology
Priority date: 2015-09-06
Filing date: 2015-09-06
Publication date: 2015-12-23
Anticipated expiration: 2035-09-06
Also published as: CN105183985B

Abstract

本发明涉及货车前轮摆振系统横向减振器最优阻尼系数的设计方法，属于货车前轮摆振系统技术领域。本发明通过利用赫尔维茨稳定性判据，计算得到横向减振器的最小临界阻尼系数，并根据轮胎原地转向阻力、横向减振器阻尼力和液压助力转向器力之间的关系，计算得到横向减振器的最大临界阻尼系数，然后，利用黄金分割原理，设计得到横向减振器的最优阻尼系数。通过设计实例及Matlab/Simulink仿真验证可知，该方法可得到准确可靠的横向减振器的最优阻尼系数设计值，为货车前轮摆振系统横向减振器最优阻尼系数的设计提供了可靠的设计方法。利用该方法，不仅可提高货车横向减振器的设计水平和车辆的行驶安全性，还可降低设计及试验费用，缩短产品开发周期。

Description

货车前轮摆振系统横向减振器最优阻尼系数的设计方法

技术领域

本发明涉及货车前轮摆振系统，特别是货车前轮摆振系统横向减振器最优阻尼系数的设计方法。

背景技术

横向减振器能够有效抑制货车高速行驶时的前轮摆振运动，显著改善车辆的行驶平顺性、安全性及操纵稳定性。然而，由于货车前轮摆振系统属于多自由度振动系统，对其进行动力学分析计算非常复杂，据所查阅资料可知，目前国内、外对于货车前轮摆振系统横向减振器阻尼系数的设计，一直未能给出可靠的解析设计方法，大都是采用实验、Matlab/Simulink建模仿真或Adams多体动力学软件建模仿真的方法对其进行设计，尽管该方法可得到比较可靠的仿真数值，使车辆具有较好的动力性能，然而，由于该方法需要进行计算机建模优化设计和计算，很难让减振器工程设计人员掌握，此外，该方法不能够提供精确的解析计算式，不能实现解析设计，更不能满足货车前轮摆振系统横向减振器CAD软件开发的要求。随着车辆行驶速度的不断提高，人们对车辆行驶稳定性和安全性提出了更高的设计要求，因此，必须建立一种精确、可靠的货车前轮摆振系统横向减振器最优阻尼系数的设计方法，满足货车前轮摆振系统横向减振器设计的要求，提高产品设计水平和质量，提高车辆行驶平顺性和安全性；同时，降低设计及试验费用，加快产品开发速度。

发明内容

针对上述现有技术中存在的缺陷，本发明所要解决的技术问题是提供一种简便、可靠的货车前轮摆振系统横向减振器最优阻尼系数的设计方法，其设计流程图如图1所示；货车前轮摆振系统三自由度行驶振动模型的前轮摆振系统图如图2所示，货车前轮摆振系统三自由度行驶振动模型的前桥简化系统图如图3所示。

为了解决上述技术问题，本发明所提供的货车前轮摆振系统横向减振器最优阻尼系数的设计方法，其特征在于采用以下设计步骤：

(1)确定前轮摆振系统赫尔维茨稳定性判据的特征行列式：

根据转向横拉杆当量角刚度K₁，前轮转向系统转向柱当量角刚度K₃；转向横拉杆等效角阻尼C₁，转向柱等效角阻尼C₃；待设计横向减振器的阻尼系数C；车辆悬架系统垂向刚度K₂，减振器等效阻尼系数C₂；轮胎垂向刚度K_t，侧向刚度ρ，侧偏刚度k，机械拖矩β；左、右前轮绕主销的转动惯量I，车轮绕自身轴线的转动惯量I_w，前桥绕侧倾轴线的转动惯量J；转向节的臂长d₁，梯形臂的臂长d₂，横向减振器的摆臂T，前悬架两弹簧中心之间的距离B_f，主销中心到车轮中心平面的距离L，轮距B，车轮半径R；主销后倾角α；轮胎滚动阻力系数f；车辆行驶速度v；利用货车前轮摆振系统三自由度行驶振动模型，以左前轮的摆振角θ₁、右前轮的摆振角θ₂，车身的侧倾角ψ为坐标，确定前轮摆振系统赫尔维茨稳定性判据的特征行列式，即：

D (s) = |\begin{matrix} q_{1} - q_{0} s & - q_{2} s - q_{3} & {Is}^{2} + (q_{6} + {CT}^{2}) s + q_{7} & 0 & q_{8} \\ q_{1} - q_{0} s & {Is}^{2} + (q_{2} + {CT}^{2}) s + q_{3} & - q_{2} s - q_{3} & q_{8} & 0 \\ {Js}^{2} + q_{9} s + q_{12} & q_{0} & q_{0} & R & R \\ 0 & ρ v k + ρ k β s & 0 & - ρ v - k s & 0 \\ 0 & 0 & ρ v k + ρ k β s & 0 & - ρ v - k s \end{matrix}|;

其中，

q_{0} = I_{w} \frac{v}{R};

q₁＝BK_tL(α-f)/2+ραR²；

q_{2} = C_{1} d_{2}^{2}; q_{3} = K_{1} d_{2}^{2}; q_{4} = C_{3} d_{1}^{2};

q_{5} = K_{3} d_{1}^{2};

q₆＝q₂+q₄；q₇＝q₃+q₅；q₈＝αR+β；

q_{9} = C_{2} B_{f}^{2} / 2; q_{10} = K_{2} B_{f}^{2} / 2;

q₁₁＝B²K_t/2+2ρR²；q₁₂＝q₁₀+q₁₁；

(2)确定赫尔维茨稳定性判据特征行列式的特征方程式：

根据步骤(1)中所确定的特征行列式D(s)，将其展开，确定赫尔维茨稳定性判据特征行列式的特征方程式，即：

a₀s⁸+a₁s⁷+a₂s⁶+a₃s⁵+a₄s⁴+a₅s³+a₆s²+a₇s+a₈＝0；

式中，

a₀＝JI²k²；

a₁＝2JIk²C+I²k²q₉+2I²Jkρv+JIk²(q₂+q₆)；

a₂＝I²k²q₁₂+C²Jk²+2CIk²q₉+(2JIkρv+Jk²q₂+CJk²+Ik²q₉)(q₂+q₆)+JIk²(q₃+q₇)+JI²ρ²v²+(2I²kq₉+4CIJk)ρv+2JIβk²q₈ρ；

a₃＝(Ik²q₉+CJk²+2IJkρv)(q₃+q₇)+Jk(kq₃+2ρvq₂)(q₆-q₂)+Jk²q₂(q₇-q₃)+2Ik²q₀(q₀+Rβρ)+[k²q₂q₉+Ik²q₁₂+Ck²q₉+IJρ²v²+Jk²βρq₈+2kρv(CJ+Iq₉)](q₂+q₆)+Ck²(Cq₉+2Iq₁₂)+2kρ(βkq₈+Cv)(CJ+Iq₉)+2JIkq₈ρv(k+βρ)+Iρ²v²(2CJ+Iq₉)+2I²kq₁₂ρv；

a₄＝k²q₀(q₀+Rβρ)(3q₂+q₆)+kq₂(kq₂+2Jρv)(q₇-q₃)+2βk²q₈ρ(Cq₉+Jβρq₈)+[k²(q₃q₉+q₂q₁₂)+2kρv(Jq₃+q₂q₉)+Jq₂ρ²v²](q₆-q₂)+Iρ²v²(2Jkq₈+Iq₁₂)+[ρ²v²(JC+Iq₉)+Jkq₈ρv²(k+βρ)+βk²q₈q₉ρ+2kρv(Cq₉+Iq₁₂)](q₆+q₂)+Cρ²v²(CJ+2Iq₉)+k²q₁₂(2Iq₈βρ+C²)+(2q₀k²+2Rβk²ρ)(Cq₀-Iq₁)+[Ck²q₉+Ik²(q₁₂+q₃)+IJρ²v²+Jk²βρq₈+2kρv(CJ+Iq₉)](q₃+q₇)+[2IRkq₀ρv+2kq₈ρv(CJ+Iq₉)](k+βρ)+kρv(2C²q₉+4Iq₀ ²+CIq₁₂)；

a₅＝q₉ρ²v²(C²+2Ikq₈)+[(2kq₀ ²ρv+Rρvq₀k)(k+βρ)-k²q₁(Rβρ+q₀)](3q₂+q₆)+[(CJ+Iq₉)ρ²v²+(k+βρ)Jkρvq₈+k²q₉(βρq₈+q₃)+2kρv(Cq₉+Iq₁₂)](q₃+q₇)+ρv[(Cq₉+Iq₁₂)ρv+2kq₁₂(q₂+C)+kq₈q₉(k+βρ)+(Jkq₈-q₂q₉)ρv](q₂+q₆)+[(Jq₂ρv+2kq₂q₉ρv+2Jkq₃)ρv+k²q₂q₁₂](q₇-q₃)+k²β²ρ²q₈(2Rq₀+q₈q₉)+q₃[(Jρv+2kq₉)ρv+k²q₁₂](q₆-q₂)+2k²βρ(Cq₈q₁₂-CRq₁+q₈q₀ ²)+2kρv[q₈(Cq₉+Iq₁₂)+R(k+βρ)](Cq₀-Iq₁)+k²q₀(Rβρ+q₀)(q₇+3q₃)+2Jkq₈ρv²(Cρ+βkq₈)+2Iq₀ρ²v²(Rk+q₀)+2Cq₁₂ρv(Ck+Iρv)-2Ck²q₀q₁；

a₆＝ρv[(q₀+Rk)q₀ρv-2kq₀q₁-Rkq₁(Rk+βρ)](q₆+3q₂)-4Ckq₀q₁ρv+ρv[(q₃q₉+q₂q₁₂)ρv+2kq₃q₁₂](q₆-q₂)+2βk²q₈ρ²v(2Rq₀+q₈q₉)+ρv[(q₂q₉+Jq₃)ρv+2k(q₃q₉+q₂q₁₂)](q₇-q₃)-2βq₁q₈ρk²(q₀+Rβρ)+[(Cq₉+Iq₁₂+Jkq₈)ρ²v²+kq₈q₉ρv(k+βρ)+kq₁₂ρ(2Cv+βkq₈)](q₇+q₃)+ρv[(Cq₁₂+kq₈q₉)ρv+kq₈q₁₂(k+βρ)](q₆+q₂)+2ρ²v²(q₀+Rk)(Cq₀-Iq₁)+[2kq₀ ²ρv-k²q₁(q₀+Rρ)+Rkq₀ρv(kq₀+βρ)](q₇+3q₃)+2kq₈ρ²v²(Cq₉+Iq₁₂)+ρv[2kq₈(Cq₁₂+q₀ ²)-2CRkq₁](k+βρ)+(C²q₁₂+Jk²q₈ ²)ρ²v²+β²k²q₈ ²q₁₂ρ²；

a₇＝[(q₀ρv-2kq₁+Rkρv)q₀ρv-Rkq₁v(k+βρ)](3q₃+q₇)-2kq₀q₁q₈ρv(k+βρ)+[(Cq₁₂+q₃q₉+kq₈q₉)ρ²v²+kq₈q₁₂ρv(k+βρ)](q₃+q₇)+kq₈ρ²v²(kq₈q₉+2Cq₁₂)+2ρ²v²(kq₈q₀-Cq₁)(q₀+Rk)+2βk²q₈ρ²v(q₈q₁₂-2Rq₁)+kq₈q₁₂ρ²v²(q₂+q₆)-q₁ρ²v²(q₀+Rk)(3q₂+q₆)+q₁₂ρv(q₂ρv+2kq₃)(q₇-q₃)+q₃q₁₂ρ²v²(q₆-q₂)；

a₈＝kq₈q₁₂ρ²v²(q₃+q₇)-q₁ρ²v²(q₀+Rk)(3q₃+q₇)+(kq₈q₁₂-2Rkq₁-2q₀q₁)kq₈ρ²v²+q₃q₁₂ρ²v²(q₇-q₃)；

其中，

q_{0} = I_{w} \frac{v}{R};

q₁＝BK_tL(α-f)/2+ραR²；

q_{2} = C_{1} d_{2}^{2}; q_{3} = K_{1} d_{2}^{2}; q_{4} = C_{3} d_{1}^{2};

q_{5} = K_{3} d_{1}^{2};

q₆＝q₂+q₄；q₇＝q₃+q₅；q₈＝αR+β；

q_{9} = C_{2} B_{f}^{2} / 2; q_{10} = K_{2} B_{f}^{2} / 2;

q₁₁＝B²K_t/2+2ρR²；q₁₂＝q₁₀+q₁₁；

(3)横向减振器最小临界阻尼系数C_min的设计：

根据步骤(2)中所确定的特征方程式，利用赫尔维茨稳定性判据及货车前轮摆振系统稳定性的临界条件，求解关于C的行列式方程

|\begin{matrix} a_{1} & a_{3} & a_{5} & a_{7} & 0 & 0 & 0 & 0 \\ a_{0} & a_{2} & a_{4} & a_{6} & a_{8} & 0 & 0 & 0 \\ 0 & a_{1} & a_{3} & a_{5} & a_{7} & 0 & 0 & 0 \\ 0 & a_{0} & a_{2} & a_{4} & a_{6} & a_{8} & 0 & 0 \\ 0 & 0 & a_{1} & a_{3} & a_{5} & a_{7} & 0 & 0 \\ 0 & 0 & a_{0} & a_{2} & a_{4} & a_{6} & a_{8} & 0 \\ 0 & 0 & 0 & a_{1} & a_{3} & a_{5} & a_{7} & 0 \\ 0 & 0 & 0 & 0 & a_{2} & a_{4} & a_{6} & a_{8} \end{matrix}| = 0

的正实数根，便可得到横向减振器的最小临界阻尼系数C_min；

(4)横向减振器最大临界阻尼系数C_max的设计：

根据转向系统的角传动比i_w，液压助力转向器作用力F_h，转向盘角速度轮胎与底面积间的滑动摩擦系数f_s，转向轴负荷G，轮胎气压P，横向减振器的摆臂T，及转向横拉杆到主销的力臂长度d₃，利用轮胎原地转向阻力、横向减振器阻尼力和液压助力转向器力之间的关系，对横向减振器的最大临界阻尼系数C_max进行设计，即：

C_{m a x} = \frac{(3 F_{h} d_{3} P - f_{s} \sqrt{G^{3} P}) i_{w}}{6 {PT}^{2} {\overset{\cdot}{θ}}_{w}};

(5)横向减振器最优阻尼系数C的设计：

根据步骤(3)中所确定的横向减振器的最小临界阻尼系数C_min，及步骤(4)中确定的横向减振器的最大临界阻尼系数C_max，利用黄金分割原理，对横向减振器的最优阻尼系数C进行设计，即：

C＝C_min+(1-0.618)(C_max-C_min)。

本发明比现有技术具有的优点：

由于货车前轮摆振系统属于多自由度振动系统，对其进行动力学分析计算非常复杂，据所查阅资料可知，目前国内、外对于货车前轮摆振系统横向减振器阻尼系数的设计，一直未能给出可靠的解析设计方法，大都是采用实验、Matlab/Simulink建模仿真或Adams多体动力学软件建模仿真的方法对其进行设计，尽管该方法可得到比较可靠的仿真数值，使车辆具有较好的动力性能，然而，由于该方法需要进行计算机建模优化设计和计算，很难让减振器工程设计人员掌握，此外，该方法不能够提供精确的解析计算式，不能实现解析设计，更不能满足货车前轮摆振系统横向减振器CAD软件开发的要求。

本发明通过利用赫尔维茨稳定性判据，计算得到横向减振器的最小临界阻尼系数，并根据轮胎原地转向阻力、横向减振器阻尼力和液压助力转向器力之间的关系，计算得到横向减振器的最大临界阻尼系数，然后，利用黄金分割原理，设计得到横向减振器的最优阻尼系数。通过设计实例及Matlab/Simulink仿真验证可知，该方法可得到准确可靠的横向减振器的最优阻尼系数设计值，为货车前轮摆振系统横向减振器最优阻尼系数的设计提供了可靠的设计方法。利用该方法，不仅可提高横向减振器的设计水平和质量，提高货车高速行驶的稳定性和安全性；同时，还可降低设计及试验费用，缩短产品开发周期。

附图说明

为了更好地理解本发明下面结合附图做进一步的说明。

图1是货车前轮摆振系统横向减振器阻尼系数设计方法的设计流程图；

图2是货车前轮摆振系统三自由度行驶振动模型的前轮摆振系统图；

图3是货车前轮摆振系统三自由度行驶振动模型的前桥简化系统图。

具体实施方案

下面通过一实施例对本发明做进一步详细说明。

某货车的转向横拉杆当量角刚度K₁＝49.1kN.m/rad，前轮转向系统转向柱当量角刚度K₃＝30.41kN.m/rad；转向横拉杆等效角阻尼C₁＝30N.m.s/rad，转向柱等效角阻尼C₃＝70N.m.s/rad；悬架系统垂向刚度K₂＝6.06kN.s/m，减振器等效阻尼系数C₂＝1.72kN.s/m；轮胎垂向刚度K_t＝850kN/m，侧向刚度ρ＝190kN/m，侧偏刚度k＝94kN.m/rad，机械拖矩β＝0.06m；左、右前轮绕主销的转动惯量I＝20kg.m²，车轮绕自身轴线的转动惯量I_w＝12.96kg.m²，前桥绕侧倾轴线的转动惯量J＝31.88kg.m²；转向节的臂长d₁＝0.85m，梯形臂的臂长d₂＝1m，横向减振器的摆臂T＝0.15m，前悬架两弹簧中心之间的距离B_f＝0.74m，主销中心到车轮中心平面的距离L＝0.07m，轮距B＝1.608m，车轮半径R＝0.48m；主销后倾角α＝0.012rad；轮胎与底面积间的滑动摩擦系数f_s＝0.7，轮胎滚动阻力系数f＝0.015；转向系统的角传动比i_w＝16，转向盘角速度液压助力转向器作用力F_h＝3500N，转向轴负荷G＝6000N，轮胎气压P＝0.25MPa，转向横拉杆到主销的力臂长度d₃＝0.18m，横向减振器的摆臂T＝0.15m；待设计横向减振器的阻尼系数为C。该货车前轮摆振系统横向减振器阻尼系数设计所要求的车辆行驶速度v＝70km/h，对该货车前轮摆振系统横向减振器的最优阻尼系数进行设计。

本发明实施例所提供的货车前轮摆振系统横向减振器最优阻尼系数的设计方法，其设计流程图如图1所示，货车前轮摆振系统三自由度行驶振动模型的前轮摆振系统图如图2所示，货车前轮摆振系统三自由度行驶振动模型的前桥简化系统图如图3所示，具体步骤如下：

(1)确定前轮摆振系统赫尔维茨稳定性判据的特征行列式：

根据转向横拉杆当量角刚度K₁＝49.1kN.m/rad，前轮转向系统转向柱当量角刚度K₃＝30.41kN.m/rad；转向横拉杆等效角阻尼C₁＝30N.m.s/rad，转向柱等效角阻尼C₃＝70N.m.s/rad；待设计横向减振器的阻尼系数C；车辆悬架系统垂向刚度K₂＝6.06kN.s/m，减振器等效阻尼系数C₂＝1.72kN.s/m；轮胎垂向刚度K_t＝850kN/m，侧向刚度ρ＝190kN/m，侧偏刚度k＝94kN.m/rad，机械拖矩β＝0.06m；左、右前轮绕主销的转动惯量I＝20kg.m²，车轮绕自身轴线的转动惯量I_w＝12.96kg.m²，前桥绕侧倾轴线的转动惯量J＝31.88kg.m²；转向节的臂长d₁＝0.85m，梯形臂的臂长d₂＝1m，横向减振器的摆臂T＝0.15m，前悬架两弹簧中心之间的距离B_f＝0.74m，主销中心到车轮中心平面的距离L＝0.07m，轮距B＝1.608m，车轮半径R＝0.48m；主销后倾角α＝0.012rad；轮胎滚动阻力系数f＝0.015；车辆行驶速度v＝70km/h；利用货车前轮摆振系统三自由度行驶振动模型，以左前轮的摆振角θ₁、右前轮的摆振角θ₂，车身的侧倾角ψ为坐标，确定前轮摆振系统赫尔维茨稳定性判据的特征行列式，即：

D (s) = |\begin{matrix} q_{1} - q_{0} s & - q_{2} s - q_{3} & {Is}^{2} + (q_{6} + {CT}^{2}) s + q_{7} & 0 & q_{8} \\ q_{1} - q_{0} s & {Is}^{2} + (q_{2} + {CT}^{2}) s + q_{3} & - q_{2} s - q_{3} & q_{8} & 0 \\ {Js}^{2} + q_{9} s + q_{12} & q_{0} & q_{0} & R & R \\ 0 & ρ v k + ρ k β s & 0 & - ρ v - k s & 0 \\ 0 & 0 & ρ v k + ρ k β s & 0 & - ρ v - k s \end{matrix}|;

其中，

q_{0} = I_{w} \frac{v}{R};

q₁＝BK_tL(α-f)/2+ραR²；

q_{2} = C_{1} d_{2}^{2}; q_{3} = K_{1} d_{2}^{2}; q_{4} = C_{3} d_{1}^{2};

q_{5} = K_{3} d_{1}^{2};

q₆＝q₂+q₄；q₇＝q₃+q₅；q₈＝αR+β；

q_{9} = C_{2} B_{f}^{2} / 2; q_{10} = K_{2} B_{f}^{2} / 2;

q₁₁＝B²K_t/2+2ρR²；q₁₂＝q₁₀+q₁₁；

(2)确定赫尔维茨稳定性判据特征行列式的特征方程式：

a₀s⁸+a₁s⁷+a₂s⁶+a₃s⁵+a₄s⁴+a₅s³+a₆s²+a₇s+a₈＝0；

式中，

a₀＝JI²k²；

a₁＝2JIk²C+I²k²q₉+2I²Jkρv+JIk²(q₂+q₆)；

其中，

q_{0} = I_{w} \frac{v}{R};

q₁＝BK_tL(α-f)/2+ραR²；

q_{2} = C_{1} d_{2}^{2}; q_{3} = K_{1} d_{2}^{2}; q_{4} = C_{3} d_{1}^{2};

q_{5} = K_{3} d_{1}^{2};

q₆＝q₂+q₄；q₇＝q₃+q₅；q₈＝αR+β；

q_{9} = C_{2} B_{f}^{2} / 2; q_{10} = K_{2} B_{f}^{2} / 2;

q₁₁＝B²K_t/2+2ρR²；q₁₂＝q₁₀+q₁₁；

(3)横向减振器最小临界阻尼系数C_min的设计：

|\begin{matrix} a_{1} & a_{3} & a_{5} & a_{7} & 0 & 0 & 0 & 0 \\ a_{0} & a_{2} & a_{4} & a_{6} & a_{8} & 0 & 0 & 0 \\ 0 & a_{1} & a_{3} & a_{5} & a_{7} & 0 & 0 & 0 \\ 0 & a_{0} & a_{2} & a_{4} & a_{6} & a_{8} & 0 & 0 \\ 0 & 0 & a_{1} & a_{3} & a_{5} & a_{7} & 0 & 0 \\ 0 & 0 & a_{0} & a_{2} & a_{4} & a_{6} & a_{8} & 0 \\ 0 & 0 & 0 & a_{1} & a_{3} & a_{5} & a_{7} & 0 \\ 0 & 0 & 0 & 0 & a_{2} & a_{4} & a_{6} & a_{8} \end{matrix}| = 0

的正实数根，可得到横向减振器的最小临界阻尼系数C_min＝1952.622N.s/m；

(4)横向减振器最大临界阻尼系数C_max的设计：

根据转向系统的角传动比i_w＝16，液压助力转向器作用力F_h＝3500N，转向盘角速度轮胎与底面积间的滑动摩擦系数f_s＝0.7，转向轴负荷G＝6000N，轮胎气压P＝0.25MPa，横向减振器的摆臂T＝0.15m，及转向横拉杆到主销的力臂长度d₃＝0.18m，利用轮胎原地转向阻力、横向减振器阻尼力和液压助力转向器力之间的关系，对横向减振器的最大临界阻尼系数C_max进行设计，即：

C_{m a x} = \frac{(3 F_{h} d_{3} P - f_{s} \sqrt{G^{3} P}) i_{w}}{6 {PT}^{2} {\overset{\cdot}{θ}}_{w}} = 2724.5 N . s / m;

(5)横向减振器最优阻尼系数C的设计：

根据步骤(3)中所确定的横向减振器的最小临界阻尼系数C_min＝1952.622N.s/m，及步骤(4)中确定的横向减振器的最大临界阻尼系数C_max＝2724.5N.s/m，利用黄金分割原理，对横向减振器的最优阻尼系数C进行设计，即：

C＝C_min+(1-0.618)(C_max-C_min)＝2247.48N.s/m。

根据实施例所提供的车辆参数，利用Matlab/Simulink仿真软件，通过建模仿真验证可知，当横向减振器的阻尼系数小于最小临界阻尼系数C_min＝1952.622N.s/m时，左前轮摆振角、右前轮摆振角、车身侧倾角的幅值越来越大；当横向减振器的阻尼系数等于最优阻尼系数C＝2247.48N.s/m时，左前轮摆振角、右前轮摆振角、车身侧倾角的幅值分别平均减小32％、33.3％、26.2％，同时，转向系统液压助力转向器所需要的作用力小于3500N，能够保证货车的转向灵活性，表明本发明所提供的货车前轮摆振系统横向减振器最优阻尼系数的设计方法是正确的。

Claims

1.货车前轮摆振系统横向减振器最优阻尼系数的设计方法，其具体步骤如下：

(1)确定前轮摆振系统赫尔维茨稳定性判据的特征行列式：

D (s) = |\begin{matrix} q_{1} - q_{0} s & - q_{2} s - q_{3} & {Is}^{2} + (q_{6} + {CT}^{2}) s + q_{7} & 0 & q_{8} \\ q_{1} - q_{0} s & {Is}^{2} + (q_{2} + {CT}^{2}) s + q_{3} & - q_{2} s - q_{3} & q_{8} & 0 \\ {Js}^{2} + q_{9} s + q_{12} & q_{0} & q_{0} & R & R \\ 0 & ρ v k + ρ k β s & 0 & - ρ v - k s & 0 \\ 0 & 0 & ρ v k + ρ k β s & 0 & - ρ v - k s \end{matrix}|;

其中，

q_{0} = I_{w} \frac{v}{R};

q₁＝BK_tL(α-f)/2+ραR²；

q_{2} = C_{1} d_{2}^{2}; q_{3} = K_{1} d_{2}^{2}; q_{4} = C_{3} d_{1}^{2};

q_{5} = K_{3} d_{1}^{2};

q₆＝q₂+q₄；q₇＝q₃+q₅；q₈＝αR+β；

q_{9} = C_{2} B_{f}^{2} / 2; q_{10} = K_{2} B_{f}^{2} / 2;

q₁₁＝B²K_t/2+2ρR²；q₁₂＝q₁₀+q₁₁；

(2)确定赫尔维茨稳定性判据特征行列式的特征方程式：

根据步骤(1)中所确定的特征行列式D(s)，将其展开，确定赫尔维茨稳定性判据特征行列式

的特征方程式，即：

a₀s⁸+a₁s⁷+a₂s⁶+a₃s⁵+a₄s⁴+a₅s³+a₆s²+a₇s+a₈＝0；

式中，

a₀＝JI²k²；

a₁＝2JIk²C+I²k²q₉+2I²Jkρv+JIk²(q₂+q₆)；

其中，

q_{0} = I_{w} \frac{v}{R};

q₁＝BK_tL(α-f)/2+ραR²；

q_{2} = C_{1} d_{2}^{2}; q_{3} = K_{1} d_{2}^{2}; q_{4} = C_{3} d_{1}^{2};

q_{5} = K_{3} d_{1}^{2};

q₆＝q₂+q₄；q₇＝q₃+q₅；q₈＝αR+β；

q_{9} = C_{2} B_{f}^{2} / 2; q_{10} = K_{2} B_{f}^{2} / 2;

q₁₁＝B²K_t/2+2ρR²；q₁₂＝q₁₀+q₁₁；

(3)横向减振器最小临界阻尼系数C_min的设计：

根据步骤(2)中所确定的特征方程式，利用赫尔维茨稳定性判据及货车前轮摆振系统稳定性

的临界条件，求解关于C的行列式方程

|\begin{matrix} a_{1} & a_{3} & a_{5} & a_{7} & 0 & 0 & 0 & 0 \\ a_{0} & a_{2} & a_{4} & a_{6} & a_{8} & 0 & 0 & 0 \\ 0 & a_{1} & a_{3} & a_{5} & a_{7} & 0 & 0 & 0 \\ 0 & a_{0} & a_{2} & a_{4} & a_{6} & a_{8} & 0 & 0 \\ 0 & 0 & a_{1} & a_{3} & a_{5} & a_{7} & 0 & 0 \\ 0 & 0 & a_{0} & a_{2} & a_{4} & a_{6} & a_{8} & 0 \\ 0 & 0 & 0 & a_{1} & a_{3} & a_{5} & a_{7} & 0 \\ 0 & 0 & 0 & 0 & a_{2} & a_{4} & a_{6} & a_{8} \end{matrix}| = 0

(4)横向减振器最大临界阻尼系数C_max的设计：

C_{m a x} = \frac{(3 F_{h} d_{3} P - f_{s} \sqrt{G^{3} P}) i_{w}}{6 {PT}^{2} {\overset{\cdot}{θ}}_{w}};

(5)横向减振器最优阻尼系数C的设计：

C＝C_min+(1-0.618)(C_max-C_min)。