CN104281759A - 内偏置非同轴式驾驶室稳定杆橡胶套长度的设计方法 - Google Patents

Abstract

本发明涉及内偏置非同轴式驾驶室稳定杆橡胶套长度的设计方法，属于驾驶室悬置技术领域。本发明根据稳定杆系统的侧倾角刚度要求值，与稳定杆的等效线刚度和橡胶衬套的载荷系数、径向刚度、等效组合线刚度之间及橡胶套长度之间的关系，建立了内偏置非同轴式稳定杆橡胶套长度设计数学模型，并利用Matlab对其进行求解设计，可得到准确可靠的橡胶套长度设计值，为驾驶室稳定杆系统提供了可靠的设计方法。利用该方法，可在不增加产品成本的前提下，仅通过橡胶衬套长度的调整设计，提高稳定杆系统的设计水平、质量和性能，满足驾驶室悬置对稳定杆侧倾角刚度的设计要求，提高车辆的行驶平顺性和安全性；同时，还可降低设计及试验费用。

Description

内偏置非同轴式驾驶室稳定杆橡胶套长度的设计方法

技术领域

本发明涉及车辆驾驶室悬置，特别是内偏置非同轴式驾驶室稳定杆橡胶套长度的设计方法。

背景技术

内偏置非同轴式驾驶室稳定杆系统是由摆臂、扭管和橡胶衬套所组成的，但却是一个由刚体、弹性体及柔性体三者组成的耦合体，其分析计算非常复杂。其中，橡胶衬套由内圆套筒、橡胶套和外圆套筒组成，而橡胶套的长度对驾驶室稳定杆系统的侧倾角刚度具有重要影响。在驾驶室实际设计中，经常采用在稳定杆结构确定情况下，通过对橡胶套长度的设计或调整，使驾驶室稳定杆系统达到侧倾角刚度的设计要求。然而，由于受橡胶衬套变形解析计算、扭管的扭转变形和弯曲变形的相互耦合，及扭转橡胶衬套载荷增加量等关键问题的制约，对于内偏置非同轴式驾驶室稳定杆橡胶套长度的调整设计，一直未能给出可靠的解析设计方法。目前，国内外对于驾驶室稳定杆系统的设计，大都是利用ANSYS仿真软件，通过实体建模对给定结构的驾驶室稳定杆系统的特性进行仿真验证，尽管该方法可得到比较可靠的仿真数值，然而，由于ANSYS仿真分析只能对给定参数的稳定杆特性进行仿真验证，无法提供精确的解析设计式，不能实现解析设计，更不能满足驾驶室稳定杆系统CAD软件开发的要求。因此，必须建立一种精确、可靠的内偏置非同轴式驾驶室稳定杆橡胶套长度的设计方法，满足驾驶室悬置及稳定杆系统侧倾角刚度的设计要求，提高产品设计水平和质量，提高车辆行驶平顺性和安全性；同时，降低设计及试验费用，加快产品开发速度。

发明内容

针对上述现有技术中存在的缺陷，本发明所要解决的技术问题是提供一种简便、可靠的内偏置非同轴式驾驶室稳定杆橡胶套长度的设计方法，其设计流程图如图1所示；内偏置非同轴式驾驶室稳定杆系统的结构示意图，如图2所示；稳定杆橡胶衬套的结构示意图如图3所示；稳定杆系统变形及摆臂位移的几何关系图如图4所示。

为解决上述技术问题，本发明所提供的内偏置非同轴式驾驶室稳定杆橡胶套长度的设计方法，其特征在于采用以下设计步骤：

(1)驾驶室稳定杆系统侧倾线刚度K_ws设计要求值的计算：

根据稳定杆的侧倾角刚度设计要求值悬置安装距离L_c，对该驾驶室稳定杆系统的侧倾线刚度K_ws的设计要求值进行计算，即

(2)内偏置非同轴式稳定杆的等效线刚度K_T的计算：

根据扭管的长度L_w，内偏置量T，内径d，外径D，弹性模量E和泊松比μ，及摆臂长度l₁，对稳定杆在驾驶室悬置安装位置处的等效线刚度K_T进行计算，即

$K_T=\frac{\mathrm{\πE}(D^4-d^4)}{32(1+\mathrm{\μ}){(l_1-T)}^2{L}_w};$

(3)内偏置非同轴式稳定杆橡胶衬套的等效组合线刚度表达式K_x(L_x)建立：

①建立橡胶衬套径向刚度表达式k_x(L_x)

根据橡胶套的内圆半径r_a，外圆半径r_b，弹性模量E_x和泊松比μ_x，以橡胶套长度L_x为待设计参变量，建立稳定杆橡胶衬套的径向刚度表达式k_x(L_x)，即

$k_x\left(L_x\right)=\frac{1}{u\left(L_x\right)+y\left(L_x\right)};$

其中， $u\left(L_x\right)=(\ln \frac{r_b}{r_a}-\frac{r_b^2-r_a^2}{r_a^2+r_b^2})\frac{1+{\mathrm{\μ}}_x}{2\mathrm{\π}{E}_x}{L}_x,$

$y\left(L_x\right)=a_{1}I(0,\mathrm{\α}{r}_b)+a_{2}K(0,\mathrm{\α}{r}_b)+a_3+\frac{1+{\mathrm{\μ}}_x}{5\mathrm{\π}{E}_x{L}_x}(\ln r_b+\frac{r_b^2}{r_a^2+r_b^2}),$

$a_1=\frac{(1+{\mathrm{\μ}}_x)[K(1,\mathrm{\α}{r}_a)r_a(r_a^2+3r_b^2)-K(1,\mathrm{\α}{r}_b)r_b({3r}_a^2+r_b^2)]}{5\mathrm{\π}{E}_x{L}_x\mathrm{\α}{r}_a{r}_b[I(1,\mathrm{\α}{r}_a)K(1,\mathrm{\α}{r}_b)-K(1,\mathrm{\α}{r}_a)I(1,\mathrm{\α}{r}_b)](r_a^2+r_b^2)},$

$a_2=\frac{({\mathrm{\μ}}_x+1)[I(1,\mathrm{\α}{r}_a)r_a(r_a^2+3r_b^2)-I(1,\mathrm{\α}{r}_b)r_b({3r}_a^2+r_b^2)]}{5\mathrm{\π}{E}_x{L}_x\mathrm{\α}{r}_a{r}_b[I(1,\mathrm{\α}{r}_a)K(1,\mathrm{\α}{r}_b)-K(1,\mathrm{\α}{r}_a)I(1,\mathrm{\α}{r}_b)](r_a^2+r_b^2)},$

$a_3=-\frac{(1+{\mathrm{\μ}}_x)(b_1-b_2+b_3)}{5\mathrm{\π}{E}_x{L}_x\mathrm{\α}{r}_a{r}_b[I(1,\mathrm{\α}{r}_a)K(1,\mathrm{\α}{r}_b)-K(1,\mathrm{\α}{r}_a)I(1,\mathrm{\α}{r}_b)](r_a^2+r_b^2)};$

$b_1=[I(1,\mathrm{\α}{r}_a)K(0,\mathrm{\α}{r}_a)+K(1,\mathrm{\α}{r}_a)I(0,\mathrm{\α}{r}_a)]r_a(r_a^2+3r_b^2),$

$b_2=[I(1,\mathrm{\α}{r}_b)K(0,\mathrm{\α}{r}_a)+K(1,\mathrm{\α}{r}_b)I(0,\mathrm{\α}{r}_a)]r_b(r_b^2+3r_a^2),$

$b_3=\mathrm{\α}{r}_a{r}_b[I(1,\mathrm{\α}{r}_a)K(1,\mathrm{\α}{r}_b)-K(1,\mathrm{\α}{r}_a)I(1,\mathrm{\α}{r}_b)][r_a^2+(r_a^2+r_b^2)\ln r_a],$

$\mathrm{\α}=2\sqrt{15}/L_x,$

Bessel修正函数I(0,αr_b)，K(0,αr_b)，I(1,αr_b)，K(1,αr_b)，

I(1,αr_a)，K(1,αr_a)，I(0,αr_a)，K(0,αr_a)；

②扭转橡胶衬套的载荷系数β_F的计算：

根据扭管的长度L_W，泊松比μ，内偏置量T，及摆臂长度l₁，对扭转橡胶衬套的载荷系数β_F进行计算，即

${\mathrm{\β}}_F=\frac{24(1+\mathrm{\μ})(l_1-T)T}{L_W^2};$

③内偏置非同轴式稳定杆橡胶衬套的等效组合线刚度表达式K_x(L_x)的建立：

根据摆臂长度l₁，扭管的内偏置量T，①步骤中所建立的橡胶衬套的径向刚度表达式k_x(L_x)，及②步骤中计算得到的扭转橡胶衬套的载荷系数β_F，建立稳定杆橡胶衬套的等效组合线刚度表达式K_x(L_x)，即

$K_X\left(L_x\right)=\frac{k_X\left(L_x\right)T}{{\mathrm{\β}}_F{l}_1+(1+{\mathrm{\β}}_F)T};$

(4)内偏置非同轴式驾驶室稳定杆橡胶套长度L_x设计数学模型的建立及其求解设计：

根据步骤(1)中计算得到的驾驶室稳定杆系统的侧倾线刚度的设计要求值K_ws，步骤(2)中计算得到的稳定杆的等效线刚度K_T，及步骤(3)中所建立的稳定杆橡胶衬套的等效组合线刚度表达式K_x(L_x)，建立内偏置非同轴式驾驶室稳定杆橡胶套长度L_x的设计数学模型，即

K_TK_X(L_x)-K_wsK_X(L_x)-K_wsK_T＝0；

利用Matlab程序，求解上述关于L_x的方程，便可得到内偏置非同轴式驾驶室稳定杆橡胶套长度L_x的设计量；

(5)内偏置非同轴式驾驶室稳定杆系统侧倾角刚度的ANSYS仿真验证：

I利用ANSYS有限元仿真软件，根据橡胶套长度L_x设计值及该驾驶室稳定杆系统的其他结构参数和材料特性参数，建立ANSYS仿真模型，划分网格，并在摆臂悬置位置施加载荷F，对稳定杆系统的变形进行ANSYS仿真，得到稳定杆系统在摆臂最外端的变形位移量f_A；

II根据橡胶套长度L_x设计值及其他结构参数和材料特性参数，利用步骤(3)的①步骤中的径向刚度计算式，计算求得所设计橡胶衬套的径向刚度k_x；

III根据ANSYS仿真所得到的摆臂最外端的变形位移量f_A，摆臂长度l₁，摆臂的悬置位置到最外端的距离Δl₁，稳定杆悬置距离L_c，在摆臂的悬置位置处所施加的载荷F，及II步骤中计算得到的橡胶衬套径向刚度k_x，利用稳定杆系统变形及摆臂位移的几何关系，对所设计的内偏置非同轴式驾驶室稳定杆系统侧倾角刚度的ANSYS仿真验证值，进行计算，即

$f_C=\frac{l_1{f}_A}{l_1+\mathrm{\Δ}{l}_1};$

$f_{\mathrm{ws}}=f_C+\frac{F}{k_x};$

$K_{\mathrm{ws}}=\frac{F}{f_{\mathrm{ws}}};$

将该非同轴式驾驶室稳定杆系统侧倾角刚度的ANSYS的仿真验证值与设计要求值进行比较，从而对本发明所提供的内偏置非同轴式驾驶室稳定杆橡胶套长度的设计方法及参数设计值进行验证。

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

由于受橡胶衬套变形解析计算、扭管的扭转变形和弯曲变形的相互耦合，及扭转橡胶衬套载荷增加量等关键问题的制约，对于内偏置非同轴式驾驶室稳定杆橡胶套长度的调整设计，一直未能给出可靠的解析设计方法。目前，国内外对于驾驶室稳定杆系统的设计，大都是利用ANSYS仿真软件，通过实体建模对给定结构的驾驶室稳定杆系统的特性进行仿真验证，尽管该方法可得到比较可靠的仿真数值，然而，由于ANSYS仿真分析只能对给定参数的稳定杆特性进行仿真验证，无法提供精确的解析设计式，不能满足驾驶室稳定杆系统CAD软件开发的要求，也不能满足实际设计生产过程中对通过橡胶套长度对稳定杆系统及侧倾角刚度进行调整设计的要求。

本发明根据内偏置非同轴式驾驶室稳定杆系统的结构及扭管内偏置量T、扭管的弯曲变形与扭转变形及载荷之间的关系，得到了稳定杆的等效线刚度K_T和扭转橡胶衬套载荷系数β_F，并以橡胶套长度L_x为待设计参变量，建立了橡胶衬套的径向刚度表达式k_x(L_x)和等效组合线刚度表达式K_x(L_x)；根据内偏置非同轴式驾驶室稳定杆系统的结构参数和材料特性参数，利用稳定杆系统侧倾角刚度，与稳定杆等效线刚度，扭转橡胶衬套载荷系数，橡胶套长度，橡胶衬套的径向刚度及等效组合线刚度之间的关系，建立了驾驶室稳定杆橡胶套长度设计数学模型；可根据驾驶室对稳定杆系统侧倾角刚度的设计要求，稳定杆系统的结构参数和材料特性参数，利用橡胶套长度L_x设计数学模型和Matlab程序，对内偏置非同轴式驾驶室稳定杆系统的橡胶套长度L_x进行求解设计。通过设计实例及ANSYS仿真验证可知，该方法可得到准确可靠的橡胶套长度L_x设计值，为驾驶室悬置及稳定杆系统的设计提供了可靠的设计方法，并且为驾驶室稳定杆系统CAD软件开发奠定了可靠的技术基础。利用该方法，不仅可提高驾驶室悬置及稳定杆系统的设计水平和质量，满足驾驶室悬置对稳定杆侧倾角刚度的设计要求，进一步提高车辆的行驶平顺性和安全性；同时，还可降低设计及试验费用，加快产品开发速度。

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

图1是内偏置非同轴式驾驶室稳定杆橡胶套长度的设计流程图；

图2是内偏置非同轴式驾驶室稳定杆系统的结构示意图；

图3是橡胶衬套的结构示意图；

图4是内偏置非同轴式稳定杆系统变形及摆臂位移的几何关系图；

图5是实施例一的橡胶衬套径向刚度k_x随橡胶套长度L_x的变化曲线；

图6是实施例一的稳定杆橡胶衬套的等效组合线刚度K_x随长度L_x刚的变化曲线；

图7是实施例一的稳定杆系统的侧倾角刚度随橡胶套长度L_x的变化曲线；

图8是实施例一的内偏置非同轴式驾驶室稳定杆系统的变形仿真云图；

图9是实施例二的橡胶衬套径向刚度k_x随橡胶套长度L_x的变化曲线；

图10是实施例二的稳定杆橡胶衬套的等效组合线刚度K_x随长度L_x刚的变化曲线；

图11是实施例二的稳定杆系统的侧倾角刚度随橡胶套长度L_x的变化曲线；

图12是实施例二的内偏置非同轴式驾驶室稳定杆系统的变形仿真云图。

具体实施方案

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

实施例一：某内偏置非同轴式驾驶室稳定杆系统的结构左右对称，如图2所示，包括：摆臂1，悬置橡胶衬套2，扭转橡胶衬套3，扭管4；其中，扭管4与扭转橡胶衬套3不同轴，扭管的内偏置量T＝30mm；左右两个摆臂1之间的距离L_c＝1550mm，即稳定杆的悬置距离；悬置橡胶衬套2与扭转橡胶衬套3之间的距离，即摆臂长度l₁＝380mm,摆臂悬置位置C到最外端A的距离Δl₁＝47.5mm；扭管4的长度L_w＝1500mm，内径d＝35mm，外径D＝50mm，弹性模量E＝200GPa，泊松比μ＝0.3；左右四个橡胶衬套的结构和材料特性完全相同，如图3所示，包括：内圆套筒5，橡胶套6，外圆套筒7，其中，内圆套筒5的内圆直径d_x＝35mm，壁厚δ＝2mm；橡胶套6的内圆半径r_a＝19.5mm，外圆半径r_b＝34.5mm，弹性模量E_x＝7.84MPa，泊松比μ_x＝0.47，长度L_x为待设计参量。该驾驶室稳定杆设计所要求的侧倾角刚度对该内偏置非同轴式驾驶室稳定杆系统的橡胶套长度L_x进行设计，并对在载荷F＝5000N下的侧倾角刚度进行ANSYS仿真验证。

本发明实例所提供的内偏置非同轴式驾驶室稳定杆橡胶套长度的设计方法，其设计流程如图1所示，具体步骤如下：

(1)驾驶室稳定杆系统侧倾线刚度K_ws设计要求值的计算：

根据稳定杆系统侧倾角刚度的设计要求值悬置距离L_c＝1550mm，对驾驶室稳定杆系统的侧倾线刚度K_ws的设计要求值进行计算，即

(2)内偏置非同轴式稳定杆的等效线刚度K_T的计算：

根据扭管的长度L_w＝1500mm，内偏置量T＝30mm，内径d＝35mm，外径D＝50mm，弹性模量E＝200GPa和泊松比μ＝0.3，及摆臂长度l₁＝380mm，对稳定杆在驾驶室悬置安装位置处的等效线刚度K_T进行计算，即

$K_T=\frac{\mathrm{\πE}(D^4-d^4)}{32(1+\mathrm{\μ}){(l_1-T)}^2{L}_W}=3.90387\×{10}^{5}N/m;$

(3)内偏置非同轴式稳定杆橡胶衬套的等效组合线刚度表达式K_x(L_x)建立：

①建立橡胶衬套径向刚度表达式k_x(L_x)

根据橡胶套的内圆半径r_a＝19.5mm，外圆半径r_b＝34.5mm，及弹性模量E_x＝7.84MPa和泊松比μ_x＝0.47，以橡胶套长度L_x为待设计参变量，建立稳定杆橡胶衬套的径向刚度表达式k_x(L_x)，即

$k_x\left(L_x\right)=\frac{1}{u\left(L_x\right)+y\left(L_x\right)};$

其中， $u\left(L_x\right)=(\ln \frac{r_b}{r_a}-\frac{r_b^2-r_a^2}{r_a^2+r_b^2})\frac{1+{\mathrm{\μ}}_x}{2\mathrm{\π}{E}_x}{L}_x,$

$y\left(L_x\right)=a_{1}I(0,\mathrm{\α}{r}_b)+a_{2}K(0,\mathrm{\α}{r}_b)+a_3+\frac{1+{\mathrm{\μ}}_x}{5\mathrm{\π}{E}_x{L}_x}(\ln r_b+\frac{r_b^2}{r_a^2+r_b^2}),$

$a_1=\frac{(1+{\mathrm{\μ}}_x)[K(1,\mathrm{\α}{r}_a)r_a(r_a^2+3r_b^2)-K(1,\mathrm{\α}{r}_b)r_b({3r}_a^2+r_b^2)]}{5\mathrm{\π}{E}_x{L}_x\mathrm{\α}{r}_a{r}_b[I(1,\mathrm{\α}{r}_a)K(1,\mathrm{\α}{r}_b)-K(1,\mathrm{\α}{r}_a)I(1,\mathrm{\α}{r}_b)](r_a^2+r_b^2)},$

$a_2=\frac{({\mathrm{\μ}}_x+1)[I(1,\mathrm{\α}{r}_a)r_a(r_a^2+3r_b^2)-I(1,\mathrm{\α}{r}_b)r_b({3r}_a^2+r_b^2)]}{5\mathrm{\π}{E}_x{L}_x\mathrm{\α}{r}_a{r}_b[I(1,\mathrm{\α}{r}_a)K(1,\mathrm{\α}{r}_b)-K(1,\mathrm{\α}{r}_a)I(1,\mathrm{\α}{r}_b)](r_a^2+r_b^2)},$

$a_3=-\frac{(1+{\mathrm{\μ}}_x)(b_1-b_2+b_3)}{5\mathrm{\π}{E}_x{L}_x\mathrm{\α}{r}_a{r}_b[I(1,\mathrm{\α}{r}_a)K(1,\mathrm{\α}{r}_b)-K(1,\mathrm{\α}{r}_a)I(1,\mathrm{\α}{r}_b)](r_a^2+r_b^2)};$

$b_1=[I(1,\mathrm{\α}{r}_a)K(0,\mathrm{\α}{r}_a)+K(1,\mathrm{\α}{r}_a)I(0,\mathrm{\α}{r}_a)]r_a(r_a^2+3r_b^2),$

$b_2=[I(1,\mathrm{\α}{r}_b)K(0,\mathrm{\α}{r}_a)+K(1,\mathrm{\α}{r}_b)I(0,\mathrm{\α}{r}_a)]r_b(r_b^2+3r_a^2),$

$b_3=\mathrm{\α}{r}_a{r}_b[I(1,\mathrm{\α}{r}_a)K(1,\mathrm{\α}{r}_b)-K(1,\mathrm{\α}{r}_a)I(1,\mathrm{\α}{r}_b)][r_a^2+(r_a^2+r_b^2)\ln r_a],$

$\mathrm{\α}=2\sqrt{15}/L_x,$

Bessel修正函数I(0,αr_b)，K(0,αr_b)，I(1,αr_b)，K(1,αr_b)，

I(1,αr_a)，K(1,αr_a)，I(0,αr_a)，K(0,αr_a)；

其中，橡胶衬套的径向刚度k_x随橡胶套长度L_x的变化曲线，如图5所示；

②扭转橡胶衬套的载荷系数β_F的计算

根据扭管的长度L_W＝1500mm，泊松比μ＝0.3，内偏置量T＝30mm，及摆臂长度l₁＝380mm，对扭转橡胶衬套的载荷系数β_F进行计算，即

${\mathrm{\β}}_F=\frac{24(1+\mathrm{\μ})(l_1-T)T}{L_W^2}=0.1456;$

③内偏置非同轴式稳定杆橡胶衬套的等效组合线刚度表达式K_x(L_x)的建立：

根据摆臂长度l₁＝380mm，扭管的内偏置量T＝30mm，①步骤中所建立的橡胶衬套径向刚度表达式k_x(L_x)，及②步骤中计算得到的β_F＝0.1456，以橡胶套长度L_x为待设计参变量，建立稳定杆橡胶衬套的等效组合线刚度表达式K_x(L_x)，即

$K_X\left(L_x\right)=\frac{k_X\left(L_x\right)T}{{\mathrm{\β}}_F{l}_1+(1+{\mathrm{\β}}_F)T};$

其中，稳定杆橡胶衬套的等效组合线刚度K_x随橡胶套长度L_x的变化曲线，如图6所示；

(4)内偏置非同轴式驾驶室稳定杆橡胶套长度L_x设计数学模型的建立及其求解设计：

根据步骤(1)中计算得到的K_ws＝2.514×10⁵N/m，步骤(2)中计算得到的K_T＝3.90387×10⁵N/m，及步骤(3)中所建立的表达式K_x(L_x)，建立内偏置非同轴式驾驶室稳定杆橡胶套长度L_x的设计数学模型，即

K_TK_X(L_x)-K_X(L_x)+K_wsK_w＝0；

利用Matlab程序，求解上述关于L_x的方程，可得到该非同轴式驾驶室稳定杆橡胶套长度L_x的设计值，即

L_x＝25mm；

其中，该稳定杆系统的侧倾角刚度随橡胶套长度L_x的变化曲线，如图7所示；

(5)内偏置非同轴式驾驶室稳定杆系统侧倾角刚度的ANSYS仿真验证：

I利用ANSYS有限元仿真软件，根据实施例中的参数及设计所得到的橡胶套长度L_x＝25mm，建立ANSYS仿真模型，划分网格，并在摆臂的悬置位置C处施加载荷F＝5000N，对稳定杆系统的变形进行ANSYS仿真，所得到的变形仿真云图，如图8所示，其中，稳定杆系统在摆臂最外端A处的变形位移量f_A为

f_A＝19.811mm；

II根据设计所的的橡胶衬套的长度L_x＝25mm及其他结构参数和材料特性参数，利用步骤(3)的①步骤中的径向刚度计算式，求得该橡胶衬套的径向刚度k_x为；

k_x＝2.1113×10⁶N/m；

III根据ANSYS仿真所得到的摆臂最外端的最大变形位移量f_A＝19.811mm，摆臂长度l₁＝380mm，摆臂的悬置位置C到最外端A的距离Δl₁＝47.5mm，稳定杆的悬置距离L_c＝1500mm，在摆臂的悬置位置C处所施加的载荷F＝5000N，及II步骤中计算得到的k_x＝2.1113×10⁶N/m，利用稳定杆系统变形及摆臂位移的几何关系，如图4所示，对该内偏置非同轴式驾驶室稳定杆系统侧倾角刚度的ANSYS仿真验证值，进行计算，即

$f_C=\frac{l_1{f}_A}{l_1+\mathrm{\Δ}{l}_1}=17.443\mathrm{mm};$

$f_{\mathrm{ws}}=f_C+\frac{F}{k_x}=19.812\mathrm{mm};$

$K_{\mathrm{ws}}=\frac{F}{f_{\mathrm{ws}}}=2.52374\×{10}^{5}N/m;$

可知，该驾驶室稳定杆的侧倾角刚度的ANSYS的仿真验证值与设计要求值相吻合，相对偏差仅为0.386％；表明该发明所提供的内偏置非同轴式驾驶室稳定杆橡胶套长度的设计方法是正确的，参数设计值是准确可靠的。

实施例二：某内偏置非同轴式驾驶室稳定杆系统的结构形式与实施例一的相同，如图2所示，其中，扭管4与扭转橡胶衬套3不同轴，扭管的内偏置量T＝30mm；左右两个摆臂1之间的距离L_c＝1400mm，即稳定杆的悬置距离；悬置橡胶衬套2与扭转橡胶衬套3之间的距离，即摆臂长度l₁＝350mm；摆臂的悬置位置C到最外端A的距离Δl₁＝52.5mm；扭管4的长度L_w＝1000mm，内径d＝42mm，外径D＝50mm；左右四个橡胶衬套的结构都完全相同，如图3所示，其中，内圆套筒5的内圆直径d_x＝35mm，壁厚δ＝5mm；橡胶套6的内圆半径r_a＝22.5mm，外圆半径r_b＝37.5mm，长度L_x为待设计参量。稳定杆的材料特性及橡胶衬套的材料特性，与实施例一的相同，即扭管的弹性模量E＝200GPa，泊松比μ＝0.3；橡胶套的弹性模量E_x＝7.84MPa，泊松比μ_x＝0.47。该驾驶室稳定杆设计所要求的侧倾角刚度对该内偏置非同轴式驾驶室稳定杆系统的橡胶套长度L_x进行设计，并对在载荷F＝5000N下的侧倾角刚度进行ANSYS仿真验证。

采用与实施例一相同的步骤，对该内偏置非同轴式驾驶室稳定杆系统的橡胶套长度L_x进行设计，即：

(1)驾驶室稳定杆系统侧倾线刚度K_ws设计要求值的计算：

根据稳定杆系统侧倾角刚度的设计要求值悬置距离L_c＝1400mm，对该驾驶室稳定杆系统侧倾线刚度K_ws的设计要求值进行计算，即

(2)内偏置非同轴式稳定杆的等效线刚度K_T的计算：

根据扭管的长度L_w＝1000mm，内偏置量T＝30mm，内径d＝42mm，外径D＝50mm，弹性模量E＝200GPa和泊松比μ＝0.3，及摆臂长度l₁＝350mm，对稳定杆在驾驶室悬置安装位置处的等效线刚度K_T进行计算，即

$K_T=\frac{\mathrm{\πE}(D^4-d^4)}{32(1+\mathrm{\μ}){(l_1-T)}^2{L}_W}=4.62894\×{10}^{5}N/m;$

(3)内偏置非同轴式稳定杆橡胶衬套的等效组合线刚度表达式K_x(L_x)建立：

①建立橡胶衬套径向刚度表达式k_x(L_x)：

根据橡胶套的内圆半径r_a＝22.5mm，外圆半径r_b＝37.5mm，弹性模量E_x＝7.84MPa，泊松比μ_x＝0.47，以橡胶套长度L_x为待设计参变量，建立稳定杆橡胶衬套的径向刚度表达式k_x(L_x)，即

$k_x\left(L_x\right)=\frac{1}{u\left(L_x\right)+y\left(L_x\right)};$

其中， $u\left(L_x\right)=(\ln \frac{r_b}{r_a}-\frac{r_b^2-r_a^2}{r_a^2+r_b^2})\frac{1+{\mathrm{\μ}}_x}{2\mathrm{\π}{E}_x}{L}_x,$

$y\left(L_x\right)=a_{1}I(0,\mathrm{\α}{r}_b)+a_{2}K(0,\mathrm{\α}{r}_b)+a_3+\frac{1+{\mathrm{\μ}}_x}{5\mathrm{\π}{E}_x{L}_x}(\ln r_b+\frac{r_b^2}{r_a^2+r_b^2}),$

$a_1=\frac{(1+{\mathrm{\μ}}_x)[K(1,\mathrm{\α}{r}_a)r_a(r_a^2+3r_b^2)-K(1,\mathrm{\α}{r}_b)r_b({3r}_a^2+r_b^2)]}{5\mathrm{\π}{E}_x{L}_x\mathrm{\α}{r}_a{r}_b[I(1,\mathrm{\α}{r}_a)K(1,\mathrm{\α}{r}_b)-K(1,\mathrm{\α}{r}_a)I(1,\mathrm{\α}{r}_b)](r_a^2+r_b^2)},$

$a_2=\frac{({\mathrm{\μ}}_x+1)[I(1,\mathrm{\α}{r}_a)r_a(r_a^2+3r_b^2)-I(1,\mathrm{\α}{r}_b)r_b({3r}_a^2+r_b^2)]}{5\mathrm{\π}{E}_x{L}_x\mathrm{\α}{r}_a{r}_b[I(1,\mathrm{\α}{r}_a)K(1,\mathrm{\α}{r}_b)-K(1,\mathrm{\α}{r}_a)I(1,\mathrm{\α}{r}_b)](r_a^2+r_b^2)},$

$a_3=-\frac{(1+{\mathrm{\μ}}_x)(b_1-b_2+b_3)}{5\mathrm{\π}{E}_x{L}_x\mathrm{\α}{r}_a{r}_b[I(1,\mathrm{\α}{r}_a)K(1,\mathrm{\α}{r}_b)-K(1,\mathrm{\α}{r}_a)I(1,\mathrm{\α}{r}_b)](r_a^2+r_b^2)};$

$b_1=[I(1,\mathrm{\α}{r}_a)K(0,\mathrm{\α}{r}_a)+K(1,\mathrm{\α}{r}_a)I(0,\mathrm{\α}{r}_a)]r_a(r_a^2+3r_b^2),$

$b_2=[I(1,\mathrm{\α}{r}_b)K(0,\mathrm{\α}{r}_a)+K(1,\mathrm{\α}{r}_b)I(0,\mathrm{\α}{r}_a)]r_b(r_b^2+3r_a^2),$

$b_3=\mathrm{\α}{r}_a{r}_b[I(1,\mathrm{\α}{r}_a)K(1,\mathrm{\α}{r}_b)-K(1,\mathrm{\α}{r}_a)I(1,\mathrm{\α}{r}_b)][r_a^2+(r_a^2+r_b^2)\ln r_a],$

$\mathrm{\α}=2\sqrt{15}/L_x,$

Bessel修正函数I(0,αr_b)，K(0,αr_b)，I(1,αr_b)，K(1,αr_b)，

I(1,αr_a)，K(1,αr_a)，I(0,αr_a)，K(0,αr_a)；

其中，橡胶衬套的径向刚度k_x随橡胶套长度L_x的变化曲线，如图9所示；

②扭转橡胶衬套的载荷系数β_F的计算：

根据扭管的长度L_W＝1000mm，泊松比μ＝0.3，内偏置量T＝30mm，及摆臂长度l₁＝350mm，对扭转橡胶衬套的载荷系数β_F进行计算，即

${\mathrm{\β}}_F=\frac{24(1+\mathrm{\μ})(l_1-T)T}{L_W^2}=0.29952;$

③内偏置非同轴式稳定杆橡胶衬套的等效组合线刚度表达式K_x(L_x)的建立：

根据摆臂长度l₁＝350mm，扭管的内偏置量T＝30mm，①步骤中所建立的橡胶衬套的橡胶衬套径向刚度表达式k_x(L_x)，及②步骤中计算得到的β_F＝0.29952，以橡胶套长度L_x为待设计参变量，建立稳定杆橡胶衬套的等效组合线刚度表达式K_x(L_x)，即

$K_X\left(L_x\right)=\frac{k_X\left(L_x\right)T}{{\mathrm{\β}}_F{l}_1+(1+{\mathrm{\β}}_F)T};$

其中，稳定杆橡胶衬套的等效组合线刚度K_x随橡胶套长度L_x的变化曲线，如图10所示；

(4)内偏置非同轴式驾驶室稳定杆橡胶套长度L_x设计数学模型的建立及其求解设计：

根据步骤(1)中计算得到的K_ws＝3.0308×10⁵N/m，步骤(2)中计算得到的K_T＝4.62894×10⁵N/m，及步骤(3)中所建立的稳定杆橡胶衬套的等效组合线刚度表达式K_x(L_x)，建立内偏置非同轴式驾驶室稳定杆橡胶套长度L_x的设计数学模型，即

K_TK_X(L_x)-K_X(L_x)+K_wsK_w＝0；

利用Matlab程序，求解上述关于L_x的方程，可得到该非同轴式驾驶室稳定杆橡胶套长度L_x的设计值，即

L_x＝40mm；

其中，该稳定杆系统的侧倾角刚度随橡胶套长度L_x的变化曲线，如图11所示；

(5)内偏置非同轴式驾驶室稳定杆系统侧倾角刚度的ANSYS仿真验证：

I利用ANSYS有限元仿真软件，根据实施例中的参数及设计所得到的橡胶套长度L_x＝40mm，建立ANSYS仿真模型，划分网格，并在摆臂悬置位置施加载荷F＝5000N，对稳定杆系统的变形进行ANSYS仿真，所得到的变形仿真云图，如图12所示，其中，稳定杆系统在摆臂最外端的最大变形位移量f_A为

f_A＝17.637mm；

II根据橡胶衬套的长度设计值L_x＝40mm及其他结构参数和材料特性参数，利用步骤(3)的①步骤中的径向刚度计算式，计算求得所设计橡胶衬套的径向刚度k_x为

k_x＝4.2085×10⁶N/m；

III根据ANSYS仿真所得到的摆臂最外端的变形位移量f_A＝17.637mm，摆臂长度l₁＝350mm，摆臂的悬置位置C到最外端A的距离Δl₁＝52.5mm，稳定杆的悬置距离L_c＝1400mm，在摆臂悬置位置C处所施加的载荷F＝5000N，及II步骤中计算得到的k_x＝4.2085×10⁶N/m，利用稳定杆系统变形及摆臂位移的几何关系，如图4所示，对该内偏置非同轴式驾驶室稳定杆系统侧倾角刚度的ANSYS仿真验证值，进行计算，即

$f_C=\frac{l_1{f}_A}{l_1+\mathrm{\Δ}{l}_1}=15.15\mathrm{mm};$

$f_{\mathrm{ws}}=f_C+\frac{F}{k_x}=16.5246\mathrm{mm};$

$K_{\mathrm{ws}}=\frac{F}{f_{\mathrm{ws}}}=3.02579\×{10}^{5}N/m;$

可知，该驾驶室稳定杆的侧倾角刚度的ANSYS的仿真验证值与设计要求值相吻合，相对偏差仅为0.166％；表明该发明所提供的内偏置非同轴式驾驶室稳定杆橡胶套长度的设计方法是正确的，参数设计值是准确可靠的。

Claims (1)

1.内偏置非同轴式驾驶室稳定杆橡胶套长度的设计方法，其具体设计步骤如下：

(1)驾驶室稳定杆系统侧倾线刚度K_ws设计要求值的计算：

根据稳定杆的侧倾角刚度设计要求值悬置安装距离L_c，对该驾驶室稳定杆系统的侧倾线刚度K_ws的设计要求值进行计算，即

(2)内偏置非同轴式稳定杆的等效线刚度K_T的计算：

根据扭管的长度L_w，内偏置量T，内径d，外径D，弹性模量E和泊松比μ，及摆臂长度l₁，对稳定杆在驾驶室悬置安装位置处的等效线刚度K_T进行计算，即

$K_T=\frac{\mathrm{\πE}(D^4-d^4)}{32(1+\mathrm{\μ}){(l_1-T)}^2{L}_W};$

(3)内偏置非同轴式稳定杆橡胶衬套的等效组合线刚度表达式K_x(L_x)建立：

①建立橡胶衬套径向刚度表达式k_x(L_x)

根据橡胶套的内圆半径r_a，外圆半径r_b，弹性模量E_x和泊松比μ_x，以橡胶套长度L_x为待设计参变量，建立稳定杆橡胶衬套的径向刚度表达式k_x(L_x)，即

$k_x\left(L_x\right)=\frac{1}{u\left(L_x\right)+y\left(L_x\right)};$

其中， $u\left(L_x\right)=(\ln \frac{r_b}{r_a}-\frac{r_b^2-r_a^2}{r_a^2+r_b^2})\frac{1+{\mathrm{\μ}}_x}{2\mathrm{\π}{E}_x}{L}_x,$

$y\left(L_x\right)=a_{1}I(0,\mathrm{\α}{r}_b)+a_{2}K(0,\mathrm{\α}{r}_b)+a_3+\frac{1+{\mathrm{\μ}}_x}{5\mathrm{\π}{E}_x{L}_x}(\ln r_b+\frac{r_b^2}{r_a^2+r_b^2}),$

$a_1=\frac{(1+{\mathrm{\μ}}_x)[K(1,\mathrm{\α}{r}_a)r_a(r_a^2+3r_b^2)-K(1,\mathrm{\α}{r}_b)r_b(3r_a^2+r_b^2)]}{5\mathrm{\π}{E}_x{L}_x\mathrm{\α}{r}_a{r}_b[I(1,\mathrm{\α}{r}_a)K(1,\mathrm{\α}{r}_b)-K(1,\mathrm{\α}{r}_a)I(1,\mathrm{\α}{r}_b)](r_a^2+r_b^2)},$

$a_2=\frac{({\mathrm{\μ}}_x+1)[I(1,\mathrm{\α}{r}_a)r_a(r_a^2+3r_b^2)-I(1,\mathrm{\α}{r}_b)r_b(3r_a^2+r_b^2)]}{5\mathrm{\π}{E}_x{L}_x\mathrm{\α}{r}_a{r}_b[I(1,\mathrm{\α}{r}_a)K(1,\mathrm{\α}{r}_b)-K(1,\mathrm{\α}{r}_a)I(1,\mathrm{\α}{r}_b)](r_a^2+r_b^2)},$

$a_3=-\frac{(1+{\mathrm{\μ}}_x)(b_1-b_2+b_3)}{5\mathrm{\π}{E}_x{L}_x\mathrm{\α}{r}_a{r}_b[I(1,\mathrm{\α}{r}_a)K(1,\mathrm{\α}{r}_b)-K(1,\mathrm{\α}{r}_a)I(1,\mathrm{\α}{r}_b)](r_a^2+r_b^2)};$

$b_1=[I(1,\mathrm{\α}{r}_a)K(0,\mathrm{\α}{r}_a)+K(1,\mathrm{\α}{r}_a)I(0,\mathrm{\α}{r}_a)]r_a(r_a^2+3r_b^2),$

$b_2=[I(1,\mathrm{\α}{r}_b)K(0,\mathrm{\α}{r}_a)+K(1,\mathrm{\α}{r}_b)I(0,\mathrm{\α}{r}_a)]r_b(r_b^2+3r_a^2),$

$b_3=\mathrm{\α}{r}_a{r}_b[I(1,\mathrm{\α}{r}_a)K(1,\mathrm{\α}{r}_b)-K(1,\mathrm{\α}{r}_a)I(1,\mathrm{\α}{r}_b)][r_a^2+(r_a^2+r_b^2)\ln r_a],$

$\mathrm{\α}=2\sqrt{15}/L_x,$

Bessel修正函数I(0,αr_b)，K(0,αr_b)，I(1,αr_b)，K(1,αr_b)，

I(1,αr_a)，K(1,αr_a)，I(0,αr_a)，K(0,αr_a)；

②扭转橡胶衬套的载荷系数β_F的计算：

根据扭管的长度L_W，泊松比μ，内偏置量T，及摆臂长度l₁，对扭转橡胶衬套的载荷系数β_F进行计算，即

${\mathrm{\β}}_F=\frac{24(1+\mathrm{\μ})(l_1-T)T}{L_W^2};$

③内偏置非同轴式稳定杆橡胶衬套的等效组合线刚度表达式K_x(L_x)的建立：

根据摆臂长度l₁，扭管的内偏置量T，①步骤中所建立的橡胶衬套的径向刚度表达式k_x(L_x)，及②步骤中计算得到的扭转橡胶衬套的载荷系数β_F，建立稳定杆橡胶衬套的等效组合线刚度表达式K_x(L_x)，即

$K_X\left(L_x\right)=\frac{k_X\left(L_x\right)T}{{\mathrm{\β}}_F{l}_1+(1+{\mathrm{\β}}_F)T};$

(4)内偏置非同轴式驾驶室稳定杆橡胶套长度L_x设计数学模型的建立及其求解设计：

根据步骤(1)中计算得到的驾驶室稳定杆系统的侧倾线刚度的设计要求值K_ws，步骤(2)中计算得到的稳定杆的等效线刚度K_T，及步骤(3)中所建立的稳定杆橡胶衬套的等效组合线刚度表达式K_x(L_x)，建立内偏置非同轴式驾驶室稳定杆橡胶套长度L_x的设计数学模型，即

K_TK_X(L_x)-K_wsK_X(L_x)-K_wsK_T＝0；

利用Matlab程序，求解上述关于L_x的方程，便可得到内偏置非同轴式驾驶室稳定杆橡胶套长度L_x的设计量；

(5)内偏置非同轴式驾驶室稳定杆系统侧倾角刚度的ANSYS仿真验证：

I 利用ANSYS有限元仿真软件，根据橡胶套长度L_x设计值及该驾驶室稳定杆系统的其他结构参数和材料特性参数，建立ANSYS仿真模型，划分网格，并在摆臂悬置位置施加载荷F，对稳定杆系统的变形进行ANSYS仿真，得到稳定杆系统在摆臂最外端的变形位移量f_A；

II 根据橡胶套长度L_x设计值及其他结构参数和材料特性参数，利用步骤(3)的①步骤中的径向刚度计算式，计算求得所设计橡胶衬套的径向刚度k_x；

III 根据ANSYS仿真所得到的摆臂最外端的变形位移量f_A，摆臂长度l₁，摆臂的悬置位置到最外端的距离Δl₁，稳定杆悬置距离L_c，在摆臂的悬置位置处所施加的载荷F，及II步骤中计算得到的橡胶衬套径向刚度k_x，利用稳定杆系统变形及摆臂位移的几何关系，对所设计的内偏置非同轴式驾驶室稳定杆系统侧倾角刚度的ANSYS仿真验证值，进行计算，即

$f_C=\frac{l_1{f}_A}{l_1+\mathrm{\Δ}{l}_1};$

$f_{\mathrm{ws}}=f_C+\frac{F}{k_x};$

$K_{\mathrm{ws}}=\frac{F}{f_{\mathrm{ws}}};$

将该非同轴式驾驶室稳定杆系统侧倾角刚度的ANSYS的仿真验证值与设计要求值进行比较，从而对本发明所提供的内偏置非同轴式驾驶室稳定杆橡胶套长度的设计方法及参数设计值进行验证。