CN104298835A - 内偏置式驾驶室稳定杆橡胶衬套内圆套筒壁厚的设计方法 - Google Patents

Abstract

本发明涉及内偏置式驾驶室稳定杆橡胶衬套内圆套筒壁厚的设计方法，属于驾驶室悬置技术领域。本发明可根据内偏置式驾驶室稳定杆系统的结构参数和材料特性参数，利用侧倾角刚度与稳定杆的等效线刚度、橡胶衬套的径向刚度和等效组合线刚度及内圆套筒壁厚之间的关系，建立内圆套筒壁厚的设计数学模型，并利用Matlab对其求解设计。通过实例设计及ANSYS仿真验证可知，该方法可得到准确可靠的内偏置式驾驶室稳定杆橡胶衬套内圆套筒壁厚的设计值，为稳定杆系统提供了可靠的设计方法。利用该方法可在不增加成本前提下，提高稳定杆系统的设计水平和性能，提高车辆行驶平顺性和安全性；同时，还可降低设计及试验费用，加快产品开发速度。

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)扭转橡胶衬套的载荷系数β_F的计算：

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

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

(4)内偏置式稳定杆橡胶衬套的等效组合线刚度表达式K_x(δ)建立：

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

根据橡胶套的长度L_x，外圆半径r_b，弹性模量E_x和泊松比μ_x，内圆套筒的内圆直径d_x，以内圆套筒壁厚δ为参变量，则橡胶套的内圆半径r_a可表示为橡胶衬套的径向刚度可表达为k_x(δ)，即

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

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

$y\left(\mathrm{\δ}\right)=a_1\left(\mathrm{\δ}\right)I(0,{\mathrm{\α}}_b)+a_2\left(\mathrm{\δ}\right)K(0,{\mathrm{\α}}_b)+{\mathrm{\α}}_3\left(\mathrm{\δ}\right)+\frac{1+{\mathrm{\μ}}_x}{5\mathrm{\π}{E}_x{L}_x}[\ln r_b+\frac{r_b^2}{{(d_x/2+\mathrm{\δ})}^2+r_b^2}],$

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

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

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

$b_1\left(\mathrm{\δ}\right)=[I(1,{\mathrm{\α}}_a)K(0,{\mathrm{\α}}_a)+K(1,{\mathrm{\α}}_a)I(0,{\mathrm{\α}}_a)](\frac{d_x}{2}+\mathrm{\δ})({(\frac{d_x}{2}+\mathrm{\δ})}^2+3r_b^2),$

$b_2\left(\mathrm{\δ}\right)=[I(1,{\mathrm{\α}}_b)K(0,{\mathrm{\α}}_a)+K(1,{\mathrm{\α}}_b)I(0,{\mathrm{\α}}_a)]\mathrm{rb}(r_b^2+3{(\frac{d}{2}+\mathrm{\δ})}^2),$

$b_3\left(\mathrm{\δ}\right)={\mathrm{\α}}_{\mathrm{ab}}[I(1,{\mathrm{\α}}_a)K(1,{\mathrm{\α}}_b)-K(1,{\mathrm{\α}}_a)I(1,{\mathrm{\α}}_b)][{(\frac{d}{2}+\mathrm{\δ})}^2+({(\frac{d}{2}+\mathrm{\δ})}^2+r_b^2)\ln (\frac{d}{2}+\mathrm{\δ})],$

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

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

$\mathrm{\α}=2\sqrt{15}/L_x,{\mathrm{\α}}_x=\mathrm{\α}(\frac{d}{2}+\mathrm{\δ}),{\mathrm{\α}}_b={\mathrm{\αr}}_b,{\mathrm{\α}}_{\mathrm{ab}}=\mathrm{\α}(\frac{d}{2}+\mathrm{\δ})r_b;$

②建立内偏置式稳定杆橡胶衬套的等效组合线刚度的表达式K_x(δ)：

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

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

(5)稳定杆橡胶衬套的内圆套筒壁厚δ设计数学模型的建立及设计：

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

K_TK_X(δ)-K_wsK_X(δ)-K_TK_ws＝0；

利用Matlab程序，求解上述关于δ的方程，便可得到内偏置式驾驶室稳定杆橡胶衬套内圆套筒壁厚δ的设计量，且橡胶套的内圆半径

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

I利用ANSYS有限元仿真软件，根据设计得到的橡胶衬套的内圆套筒壁厚δ和橡胶套的内圆半径r_a，及该内偏置式驾驶室稳定杆系统的其他结构参数和材料特性参数，建立相应的ANSYS仿真模型，划分网格，并在摆臂的悬置位置处施加载荷F，对稳定杆系统的变形进行ANSYS仿真，得到稳定杆系统在摆臂最外端的变形位移量f_A；

II根据步骤(5)中设计得到的内圆套筒壁厚δ和橡胶套的内圆半径r_a，及橡胶衬套的其他结构参数和材料特性参数，利用步骤(4)中的①步骤中的径向刚度计算式，计算得到所设计橡胶衬套的径向刚度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，并以橡胶衬套的内圆套筒壁厚δ为待设计参变量，建立了橡胶衬套的径向刚度表达式k_x(δ)和橡胶衬套的等效组合线刚度表达式K_x(δ)；根据内偏置式驾驶室稳定杆系统的结构参数和材料特性参数，利用稳定杆系统的侧倾线刚度，与稳定杆的等效线刚度和橡胶衬套的等效组合线刚度之间的关系，建立了驾驶室稳定杆橡胶衬套的内圆套筒壁厚δ的设计数学模型；可根据驾驶室对稳定杆系统侧倾角刚度的设计要求，稳定杆系统的结构参数和材料特性参数，利用橡胶衬套的内圆套筒壁厚δ的设计数学模型和Matlab计算程序，对内偏置式驾驶室稳定杆橡胶衬套的内圆套筒壁厚δ进行求解设计。通过设计实例及ANSYS仿真验证可知，该方法可得到准确可靠的内圆套筒壁厚δ的设计值，为驾驶室悬置及稳定杆系统的设计提供了可靠的设计方法，并且为驾驶室稳定杆系统CAD软件开发奠定了可靠的技术基础。利用该方法，可在不增加产品成本的情况下，提高驾驶室悬置及稳定杆系统的设计水平、质量和性能，满足驾驶室悬置对稳定杆系统侧倾刚度的设计要求，进一步提高车辆的行驶平顺性和安全性；同时，还降低设计及试验费用，加快产品开发速度。

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

图1是内偏置式驾驶室稳定杆橡胶衬套内圆套筒壁厚的设计流程图；

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

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

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

图5是实施例一的橡胶衬套的径向刚度k_x随内圆套筒壁厚δ的变化曲线；

图6是实施例一的稳定杆橡胶衬套的等效组合线刚度K_x随内圆套筒壁厚δ的变化曲线；

图7是实施例一的稳定杆系统侧倾角刚度随内圆套筒壁厚δ的变化曲线；

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

图9是实施例二的橡胶衬套的径向刚度k_x随内圆套筒壁厚δ的变化曲线；

图10是实施例二的稳定杆橡胶衬套的等效组合线刚度K_x随内圆套筒壁厚δ的变化曲线；

图11是实施例二的稳定杆系统侧倾角刚度随内圆套筒壁厚δ的变化曲线；

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

具体实施方案

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

实施例一：某内偏置式驾驶室稳定杆系统的结构左右对称，如图2所示，包括：摆臂1，悬置橡胶衬套2，扭转橡胶衬套3，扭管4；其中，扭管4与扭转橡胶衬套3不同轴，扭管4的内偏置量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，壁厚δ为待设计参量；橡胶套6的外圆半径r_b＝34.5mm，弹性模量E_x＝7.84MPa，泊松比μ_x＝0.47。该驾驶室稳定杆设计所要求的侧倾角刚度对该内偏置式驾驶室稳定杆系统橡胶衬套的内圆套筒壁厚δ进行设计，并对在载荷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)扭转橡胶衬套的载荷系数β_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;$

(4)内偏置式稳定杆橡胶衬套的等效组合线刚度表达式K_x(δ)建立：

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

根据橡胶套的长度L_x＝25mm，外圆半径r_b＝34.5mm，弹性模量E_x＝7.84MPa和泊松比μ_x＝0.47，内圆套筒的内圆直径d_x＝35mm，以内圆套筒壁厚δ为参变量，则橡胶套的内圆半径r_a可表示为r_a＝d_x/2+δ，橡胶衬套的径向刚度的表达式k_x(δ)，可表示为：

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

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

$y\left(\mathrm{\δ}\right)=a_1\left(\mathrm{\δ}\right)I(0,{\mathrm{\α}}_b)+a_2\left(\mathrm{\δ}\right)K(0,{\mathrm{\α}}_b)+{\mathrm{\α}}_3\left(\mathrm{\δ}\right)+\frac{1+{\mathrm{\μ}}_x}{5\mathrm{\π}{E}_x{L}_x}[\ln r_b+\frac{r_b^2}{{(d_x/2+\mathrm{\δ})}^2+r_b^2}],$

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

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

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

$b_1\left(\mathrm{\δ}\right)=[I(1,{\mathrm{\α}}_a)K(0,{\mathrm{\α}}_a)+K(1,{\mathrm{\α}}_a)I(0,{\mathrm{\α}}_a)](\frac{d_x}{2}+\mathrm{\δ})({(\frac{d_x}{2}+\mathrm{\δ})}^2+3r_b^2),$

$b_2\left(\mathrm{\δ}\right)=[I(1,{\mathrm{\α}}_b)K(0,{\mathrm{\α}}_a)+K(1,{\mathrm{\α}}_b)I(0,{\mathrm{\α}}_a)]\mathrm{rb}(r_b^2+3{(\frac{d}{2}+\mathrm{\δ})}^2),$

$b_3\left(\mathrm{\δ}\right)={\mathrm{\α}}_{\mathrm{ab}}[I(1,{\mathrm{\α}}_a)K(1,{\mathrm{\α}}_b)-K(1,{\mathrm{\α}}_a)I(1,{\mathrm{\α}}_b)][{(\frac{d}{2}+\mathrm{\δ})}^2+({(\frac{d}{2}+\mathrm{\δ})}^2+r_b^2)\ln (\frac{d}{2}+\mathrm{\δ})],$

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

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

$\mathrm{\α}=2\sqrt{15}/L_x,{\mathrm{\α}}_x=\mathrm{\α}(\frac{d}{2}+\mathrm{\δ}),{\mathrm{\α}}_b={\mathrm{\αr}}_b,{\mathrm{\α}}_{\mathrm{ab}}=\mathrm{\α}(\frac{d}{2}+\mathrm{\δ})r_b;$

其中，在外圆半径r_b及长度L_x给定情况下，橡胶衬套的径向刚度k_x随内圆套筒壁厚δ的变化曲线，如图5所示；

②建立内偏置式稳定杆橡胶衬套的等效组合线刚度的表达式K_x(δ)：

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

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

其中，稳定杆橡胶衬套的等效组合线刚度K_x随内圆套筒壁厚δ的变化曲线，如图6所示；

(5)稳定杆橡胶衬套的内圆套筒壁厚δ设计数学模型的建立及设计：

根据步骤(1)中计算得到的K_ws＝2.514×10⁵N/m，步骤(2)中计算得到的K_T＝3.90387×10⁵N/m，及步骤(4)中所建立的K_x(δ)＝0.33446k_X(δ)，建立稳定杆橡胶衬套的内圆套筒壁厚δ的设计数学模型，即

K_TK_X(δ)-K_wsK_X(δ)-K_TK_ws＝0；

利用Matlab程序，求解上述关于δ的方程，可得到该内偏置式驾驶室稳定杆橡胶衬套内圆套筒壁厚δ，即

δ＝2.0mm；

其中，橡胶套的内圆半径r_a＝d_x/2+δ＝19.5mm；稳定杆系统侧倾角刚度随内圆套筒壁厚δ的变化曲线，如图7所示；

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

I利用ANSYS有限元仿真软件，根据设计得到的橡胶衬套的内圆套筒壁厚δ＝2.0mm和内圆半径r_a＝19.5mm，及该内偏置式驾驶室稳定杆系统的其他结构参数和材料特性参数，建立ANSYS仿真模型，划分网格，并在摆臂的悬置位置C处施加载荷F＝5000N，对稳定杆系统的变形进行ANSYS仿真，所得到的变形仿真云图，如图8所示，其中，稳定杆系统在摆臂最外端A处的变形位移量f_A为

f_A＝19.811mm；

II根据步骤(5)中设计得到的内圆套筒壁厚δ＝2.0mm和橡胶套的内圆半径r_a＝19.5mm，及橡胶衬套的其他结构参数和材料特性参数，利用步骤(4)中的①步骤中的径向刚度计算式，计算得到所设计橡胶衬套的径向刚度k_x为；

k_x＝2.1113×10⁶N/m；

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

采用与实施例一相同的步骤，对该内偏置式驾驶室稳定杆橡胶衬套的内圆套筒壁厚δ进行设计，即：

(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}\text{=4.62894\×}{10}^{5}N/m;$

(3)扭转橡胶衬套的载荷系数β_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;$

(4)内偏置式稳定杆橡胶衬套的等效组合线刚度表达式K_x(δ)建立：

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

根据橡胶套的长度L_x＝40mm，外圆半径r_b＝37.5mm，E_x＝7.84MPa和泊松比μ_x＝0.47，内圆套筒的内圆直径d_x＝35mm，以内圆套筒壁厚δ为参变量，则橡胶套的内圆半径r_a可表示为r_a＝d_x/2+δ，橡胶衬套的径向刚度可表达为k_x(δ)，即

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

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

$y\left(\mathrm{\δ}\right)=a_1\left(\mathrm{\δ}\right)I(0,{\mathrm{\α}}_b)+a_2\left(\mathrm{\δ}\right)K(0,{\mathrm{\α}}_b)+{\mathrm{\α}}_3\left(\mathrm{\δ}\right)+\frac{1+{\mathrm{\μ}}_x}{5\mathrm{\π}{E}_x{L}_x}[\ln r_b+\frac{r_b^2}{{(d_x/2+\mathrm{\δ})}^2+r_b^2}],$

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

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

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

$b_1\left(\mathrm{\δ}\right)=[I(1,{\mathrm{\α}}_a)K(0,{\mathrm{\α}}_a)+K(1,{\mathrm{\α}}_a)I(0,{\mathrm{\α}}_a)](\frac{d_x}{2}+\mathrm{\δ})({(\frac{d_x}{2}+\mathrm{\δ})}^2+3r_b^2),$

$b_2\left(\mathrm{\δ}\right)=[I(1,{\mathrm{\α}}_b)K(0,{\mathrm{\α}}_a)+K(1,{\mathrm{\α}}_b)I(0,{\mathrm{\α}}_a)]\mathrm{rb}(r_b^2+3{(\frac{d}{2}+\mathrm{\δ})}^2),$

$b_3\left(\mathrm{\δ}\right)={\mathrm{\α}}_{\mathrm{ab}}[I(1,{\mathrm{\α}}_a)K(1,{\mathrm{\α}}_b)-K(1,{\mathrm{\α}}_a)I(1,{\mathrm{\α}}_b)][{(\frac{d}{2}+\mathrm{\δ})}^2+({(\frac{d}{2}+\mathrm{\δ})}^2+r_b^2)\ln (\frac{d}{2}+\mathrm{\δ})],$

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

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

$\mathrm{\α}=2\sqrt{15}/L_x,{\mathrm{\α}}_x=\mathrm{\α}(\frac{d}{2}+\mathrm{\δ}),{\mathrm{\α}}_b={\mathrm{\αr}}_b,{\mathrm{\α}}_{\mathrm{ab}}=\mathrm{\α}(\frac{d}{2}+\mathrm{\δ})r_b;$

其中，在外圆半径r_b及长度L_x给定情况下，橡胶衬套的径向刚度k_x随内圆套筒壁厚δ的变化曲线，如图9所示；

②建立内偏置式稳定杆橡胶衬套的等效组合线刚度的表达式K_x(δ)：

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

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

其中，稳定杆橡胶衬套的等效组合线刚度K_x随内圆套筒壁厚δ的变化曲线，如图10所示；

(5)稳定杆橡胶衬套的内圆套筒壁厚δ设计数学模型的建立及设计：

根据步骤(1)中计算得到的K_ws＝3.0308×10⁵N/m，步骤(2)中计算得到的K_T＝4.62894×10⁵N/m，及步骤(4)中所建立的K_x(δ)＝0.20860k_X(δ)，建立稳定杆橡胶衬套的内圆套筒壁厚δ的设计数学模型，即

K_TK_X(δ)-K_wsK_X(δ)-K_TK_ws＝0；

利用Matlab程序，求解上述关于δ的方程，可得到该内偏置式驾驶室稳定杆橡胶衬套内圆套筒壁厚δ，即

δ＝5.0mm；

其中，橡胶套的内圆半径r_a＝d_x/2+δ＝22.5mm；稳定杆系统侧倾角刚度随内圆套筒壁厚δ的变化曲线，如图11所示；

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

I利用ANSYS有限元仿真软件，根据设计得到的橡胶衬套的内圆套筒壁厚δ＝5.0mm和内圆半径r_a＝22.5mm，及该内偏置式驾驶室稳定杆系统的其他结构参数和材料特性参数，建立ANSYS仿真模型，划分网格，并在摆臂的悬置位置C处施加载荷F＝5000N，对稳定杆系统的变形进行ANSYS仿真，所得到的变形仿真云图，如图12所示，其中，稳定杆系统在摆臂最外端A处的变形位移量f_A为

f_A＝17.637mm；

II根据步骤(5)中设计得到的内圆套筒壁厚δ＝5.0mm和橡胶套的内圆半径r_a＝22.5mm，及橡胶衬套的其他结构参数和材料特性参数，利用步骤(4)中的①步骤中的径向刚度计算式，计算得到所设计橡胶衬套的径向刚度k_x为；

k_x＝4.2085×10⁶N/m；

III根据ANSYS仿真所得到的摆臂最外端A处的变形位移量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)扭转橡胶衬套的载荷系数β_F的计算：

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

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

(4)内偏置式稳定杆橡胶衬套的等效组合线刚度表达式K_x(δ)建立：

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

根据橡胶套的长度L_x，外圆半径r_b，弹性模量E_x和泊松比μ_x，内圆套筒的内圆直径d_x，以内圆套筒壁厚δ为参变量，则橡胶套的内圆半径r_a可表示为橡胶衬套的径向刚度可表达为k_x(δ)，即

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

其中， $u\left(\mathrm{\δ}\right)=[1n\frac{r_b}{(\frac{d_x}{0}+\mathrm{\δ})}-\frac{r_b^2-{(\frac{d_x}{2}+\mathrm{\δ})}^2}{(\frac{d_x}{2}+\mathrm{\δ})+r_b^2}]\frac{1+{\mathrm{\μ}}_x}{2\mathrm{\π}{E}_x}{L}_x,$

$y\left(\mathrm{\δ}\right)=a_1\left(\mathrm{\δ}\right)I(0,{\mathrm{\α}}_b)+a_2\left(\mathrm{\δ}\right)K(0,{\mathrm{\α}}_b)+a_3\left(\mathrm{\δ}\right)+\frac{1+{\mathrm{\μ}}_x}{5\mathrm{\π}{E}_x{L}_x}[1n{r}_b+\frac{r_b^2}{{(d_x/2+\mathrm{\δ})}^2+r_b^2}],$ $a_1\left(\mathrm{\δ}\right)=\frac{(1-{\mathrm{\μ}}_x)[K(1,\mathrm{\α}(\frac{d_x}{2}+\mathrm{\δ}))(\frac{d_x}{2}+\mathrm{\δ})({(\frac{d_x}{2}+\mathrm{\δ})}^2+{3r}_b^2)-K(1,{\mathrm{\α}}_b)r_b(3{(\frac{d_x}{2}+\mathrm{\δ})}^2+r_b^2)]}{5\mathrm{\π}{E}_x{L}_x{\mathrm{\α}}_{\mathrm{ab}}[I(1,{\mathrm{\α}}_a)K(1,{\mathrm{ar}}_b)-K(1,{\mathrm{\α}}_a)I(1,{\mathrm{\α}}_b)]({(\frac{d_x}{2}+\mathrm{\δ})}^2+r_b^2)},$

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

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

$b_1\left(\mathrm{\δ}\right)=[I(1,{\mathrm{\α}}_a)K(0,{\mathrm{\α}}_a)+K(1,{\mathrm{\α}}_a)I(0,{\mathrm{\α}}_a)](\frac{d_x}{2}+\mathrm{\δ})({(\frac{d_x}{2}+\mathrm{\δ})}^2+{3r}_b^2),$

$b_2\left(\mathrm{\δ}\right)=[I(1,{\mathrm{\α}}_b)K(0,{\mathrm{\α}}_a)+K(1,{\mathrm{\α}}_b)I(0,{\mathrm{\α}}_a)]r_b(r_b^2+3{(\frac{d}{2}+\mathrm{\δ})}^2),$

$b_3\left(\mathrm{\δ}\right)={\mathrm{\α}}_{\mathrm{ab}}[I(1,{\mathrm{\α}}_a)K(1,{\mathrm{\α}}_b)-K(1,{\mathrm{\α}}_a)I(1,{\mathrm{\α}}_b)][{(\frac{d}{2}+\mathrm{\δ})}^2+({(\frac{d}{2}+\mathrm{\δ})}^2+r_b^2)\ln (\frac{d}{2}+\mathrm{\δ})],$

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

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

$\mathrm{\α}=2\sqrt{15}/L_x,{\mathrm{\α}}_a=\mathrm{\α}(\frac{d}{2}+\mathrm{\δ}),{\mathrm{\α}}_b={\mathrm{\αr}}_b,{\mathrm{\α}}_{\mathrm{ab}}=\mathrm{\α}(\frac{d}{2}+\mathrm{\δ})r_b;$

②建立内偏置式稳定杆橡胶衬套的等效组合线刚度的表达式Kx(δ)：

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

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

(5)稳定杆橡胶衬套的内圆套筒壁厚δ设计数学模型的建立及设计：

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

K_TK_X(δ)-K_wsK_X(δ)-K_TK_ws＝0；

利用Matlab程序，求解上述关于δ的方程，便可得到内偏置式驾驶室稳定杆橡胶衬套内圆套筒壁厚δ的设计量，且橡胶套的内圆半径

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

I利用ANSYS有限元仿真软件，根据设计得到的橡胶衬套的内圆套筒壁厚δ和橡胶套的内圆半径r_a，及该内偏置式驾驶室稳定杆系统的其他结构参数和材料特性参数，建立相应的ANSYS仿真模型，划分网格，并在摆臂的悬置位置处施加载荷F，对稳定杆系统的变形进行ANSYS仿真，得到稳定杆系统在摆臂最外端的变形位移量f_A；

II根据步骤(5)中设计得到的内圆套筒壁厚δ和橡胶套的内圆半径r_a，及橡胶衬套的其他结构参数和材料特性参数，利用步骤(4)中的①步骤中的径向刚度计算式，计算得到所设计橡胶衬套的径向刚度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仿真验证值与设计要求值进行比较，从而对本发明所提供的内偏置式驾驶室稳定杆橡胶衬套内圆套筒壁厚的设计方法及参数设计值进行验证。