秦武 李春歸 潘兵兵 李駿 劉霏霏

收稿日期：2023?10?19

基金項目：國家自然科學基金（51806066）；博士后科學基金項目（2022M711352）；江西省自然科學基金資助項目（20232BAB214047）；江西省教育廳科技項目（GJJ210630）

文章編號：1005-0523（2024）03-0090-11

摘要：【目的】通過對橡膠襯套各向剛度的優(yōu)化設(shè)計以指導(dǎo)產(chǎn)品設(shè)計，以期大幅縮短橡膠襯套的研發(fā)周期。橡膠襯套作為新能源汽車、特種車輛、動車組等載運工具的重要連接件，對整車操縱穩(wěn)定性起著至關(guān)重要的作用?！痉椒ā繛榱吮苊夤こ虒嵺`中反復(fù)調(diào)試的繁瑣工作，提出了一種考慮材料硬度與縮徑量的橡膠襯套各向剛度（徑向剛度、軸向剛度、偏擺剛度和扭轉(zhuǎn)剛度）優(yōu)化設(shè)計方法。首先，采用Ogden模型作為橡膠襯套的本構(gòu)模型，通過試驗設(shè)計與有限元分析相結(jié)合的方法，計算不同內(nèi)芯外徑、橡膠高度、縮徑量和材料硬度下襯套各向剛度；然后，建立各向剛度的二階響應(yīng)面模型，利用拉丁超立方抽樣和方差分析，檢驗響應(yīng)面模型精度；最后，根據(jù)遺傳算法，建立含各向剛度為設(shè)計變量的多目標優(yōu)化模型，并對優(yōu)化結(jié)果進行試驗驗證。【結(jié)果】研究表明：徑向剛度、軸向剛度、偏擺剛度和扭轉(zhuǎn)剛度的試驗結(jié)果與目標值的相對誤差分別為7.72%，9.06%，-6.33%和9.16%，均在±10%以內(nèi)，滿足工程應(yīng)用的要求?！窘Y(jié)論】試驗結(jié)果驗證了所建立的各向剛度二階響應(yīng)面模型的有效性，以及所提出優(yōu)化設(shè)計方法的可行性，為橡膠襯套產(chǎn)品設(shè)計提供指導(dǎo)，可極大地縮短橡膠襯套的研發(fā)周期。

關(guān)鍵詞：橡膠襯套；響應(yīng)面方法；遺傳算法；多目標優(yōu)化

中圖分類號：U467.3 文獻標志碼：A

本文引用格式：秦武，李春歸，潘兵兵，等. 考慮材料硬度與縮徑量的橡膠襯套各向剛度優(yōu)化設(shè)計[J]. 華東交通大學學報，2024，41（3）：90-100.

Optimal Design of Isotropic Stiffness for Rubber Bushing by Considering Material Hardness and Diameter

Reduction for Conveyance

Qin Wu1，2， Li Chungui1， Pan Bingbing2， Li Jun1， Liu Feifei1

（1. School of Mechatronics & Vehicle Engineering， East China Jiaotong University， Nanchang 330013， China；

2. Jianxin Zhao's Technology Co. Ltd.， Ningbo 315615， China）

Abstract： 【Objective】As an important connecting part of new energy vehicles， special vehicles， locomotives and other means of transportation， rubber bushing plays a crucial role in the stability of the whole vehicle handling. In order to avoid the tedious work of repeated debugging in engineering practice， 【Method】an optimal design method of rubber bushing all-directional stiffness （radial， axial， yaw and torsional stiffness） taking into account the material hardness and diameter reduction is proposed. Firstly， the Ogden model is used as the constitutive model of rubber bushing. Second， by combining experimental design with finite element analysis， the isotropic stiffness of rubber bushing with different inner core outer diameter， rubber height， diameter reduction and material hardness is calculated. The second-order response surface model of isotropic stiffness is established and the accuracy of the model is verified by Latin hypercube sampling and variance analysis. Finally， a multi-objective optimization method of isotropic stiffness of rubber bushing is proposed using genetic algorithm to testify the optimization result. 【Result】 The result shows that the relative errors between the measured radial， axial， yaw and torsional stiffness and objective stiffness are 7.72%， 9.06%， -6.33% and 9.16% respectively， which are all within ±10% and meet the requirements of engineering application. 【Conclusion】 The validity of the established second-order response surface model of isotropic stiffness and the feasibility of the proposed optimization design method are verified to provide guidance for the product design of rubber bushing， which can greatly shorten the research and development cycle of rubber bushing.

Key words： rubber bushing; response surface method; genetic algorithm; multi-objective optimization.

Citation format： QIN W， LI C G， PAN B B， et al. Optimal design of isotropic stiffness for rubber bushing by considering material hardness and diameter reduction for conveyance[J]. Journal of East China Jiaotong University， 2024， 41（3）： 90-100.

【研究意義】橡膠襯套作為一種重要的連接元件，主要安裝在懸架、副車架、轉(zhuǎn)向架、穩(wěn)定桿等位置，被廣泛應(yīng)用于新能源汽車、特種車輛、動車組等載運工具和軌道交通領(lǐng)域中[1]。橡膠襯套材料硬度與結(jié)構(gòu)決定著襯套的剛度。若剛度設(shè)計過大，對載運工具振動與噪聲控制不利；若剛度設(shè)計過小，載運工具或其子系統(tǒng)行程太大，會影響整車操縱穩(wěn)定性。因此，優(yōu)化設(shè)計襯套剛度對整車性能的提升起著至關(guān)重要作用。

【研究進展】在橡膠襯套剛度影響因素研究中，縮徑量影響著襯套的各向剛度，在工程上，縮徑量可定義為：橡膠襯套外套直徑的壓縮變化量。方紅和陳寶等[2-3]研究表明，可以合理控制縮徑量尺寸來調(diào)整橡膠襯套剛度，當縮徑值達到一定程度后對橡膠襯套靜、動剛度值的影響變大。曾家興和陳寶等[4-5]研究發(fā)現(xiàn)，不同材料硬度對橡膠彈簧超彈本構(gòu)模型參數(shù)和彈簧垂向剛度有顯著影響，且硬度越大，靜剛度變化越明顯。

在橡膠襯套剛度優(yōu)化算法研究中，目前有基于各種優(yōu)化理論的算法。較為成熟的有差分進化算法[6]、自適應(yīng)模擬退火算法[7]、新型混沌磷蝦群優(yōu)化算法[8]和遺傳算法[9]等，被廣泛應(yīng)用于橡膠襯套剛度優(yōu)化中。目前，國內(nèi)外對于橡膠襯套剛度優(yōu)化開展了研究工作。

Dai等[10]采用混合神經(jīng)網(wǎng)絡(luò)模型開發(fā)了橡膠襯套的詳細模型，并識別其關(guān)鍵參數(shù)。Jung等[11]基于襯套形狀設(shè)計變量與有限元分析預(yù)測的剛度值之間的關(guān)系，在人工神經(jīng)網(wǎng)絡(luò)模型中訓練了多層感知器和卷積神經(jīng)網(wǎng)絡(luò)模型，準確預(yù)測橡膠襯套剛度?？驴〉萚12]通過有限元法計算推力桿球鉸橡膠襯套的應(yīng)變分布特征和剛度參數(shù)，建立推力桿球鉸的多目標優(yōu)化模型，采用遺傳算法得到最優(yōu)球鉸結(jié)構(gòu)參數(shù)的優(yōu)化方案。郭一鳴等[13]對懸置中非橡膠填充結(jié)構(gòu)進行優(yōu)化設(shè)計，通過謝菲爾德遺傳算法對結(jié)構(gòu)參數(shù)進行優(yōu)化，以滿足動力總成懸置系統(tǒng)匹配過程中懸置的剛度要求。

【創(chuàng)新特色】綜上所述，已有相關(guān)研究通過優(yōu)化算法，對橡膠襯套結(jié)構(gòu)參數(shù)進行優(yōu)化，設(shè)計橡膠襯套的剛度，改善整車操縱穩(wěn)定性。但是，目前基于橡膠襯套剛度優(yōu)化研究方面還存在以下問題：鮮有綜合考慮襯套各向剛度優(yōu)化設(shè)計方面的研究，絕大多數(shù)研究只考慮襯套徑向剛度和軸向剛度，而忽略了偏擺剛度和扭轉(zhuǎn)剛度，可能會導(dǎo)致車輛在俯仰及側(cè)傾方向上行程過大，影響車輛操縱穩(wěn)定性；絕大多數(shù)研究從襯套結(jié)構(gòu)角度上只針對材料硬度或縮徑量進行單一變量分析，并未考慮二者對橡膠襯套各向剛度的綜合影響。

【關(guān)鍵問題】因此，本文以某車輛懸架中橡膠襯套為研究對象，提出了一種考慮材料硬度與縮徑量的橡膠襯套各向剛度優(yōu)化設(shè)計方法。采用試驗設(shè)計與有限元分析相結(jié)合的方法，并通過試驗驗證有限元模型的準確性。建立各向剛度與參數(shù)變量的二階響應(yīng)面模型，以拉丁超立方抽樣和方差分析對響應(yīng)面模型進行精度檢驗。利用遺傳算法實現(xiàn)多目標優(yōu)化獲取最優(yōu)參數(shù)，并將優(yōu)化結(jié)果進行試驗驗證。試驗結(jié)果表明：各向剛度試驗結(jié)果與目標剛度的相對誤差分別為7.72%，9.06%，-6.33%和9.16%，均在±10%以內(nèi)，滿足橡膠襯套剛度匹配要求。

1 橡膠襯套有限元建模及試驗驗證

1.1 橡膠襯套有限元模型建立

橡膠襯套有限元建模過程如圖1所示。其中橡膠襯套主要由金屬部分的內(nèi)芯、外套和橡膠部分組成，如圖1（a）所示。所研究的橡膠襯套結(jié)構(gòu)在徑向和軸向上采用了對稱設(shè)計，兩個徑向剛度和兩個偏擺剛度在相同載荷作用下結(jié)果相同，即所述的橡膠襯套4個方向剛度可以代表6個方向剛度。故本文采用“各向剛度”表述。

選取對橡膠襯套各向剛度影響較大的因素作為設(shè)計變量，考慮內(nèi)芯外徑[a]、橡膠高度[b]、縮徑量[c]和材料硬度[d]作為參數(shù)變量[14]，如圖1（b）所示。

橡膠襯套在承載過程中主要是橡膠部分發(fā)生變形，金屬部分的變形很小，可以忽略不計。為了簡化橡膠襯套模型，在剛度計算時只對橡膠部分進行有限元分析。

根據(jù)圖1（a）所示的橡膠襯套結(jié)構(gòu)，先利用Hypermesh軟件對橡膠襯套幾何模型進行網(wǎng)格劃分，再將網(wǎng)格導(dǎo)入Abaqus軟件中，建立如圖1（c）和圖1（d）所示的有限元模型。圖1（c）為邊界條件，圖1（d）為載荷設(shè)置，加載點1為內(nèi)表面與形心處參考點的剛性耦合點，加載面2為襯套施加載荷的外表面。

不同的材料硬度影響橡膠襯套超彈性本構(gòu)模型參數(shù)，使得有限元仿真計算的剛度值相差較大。通常采用應(yīng)變能密度函數(shù)來準確描述橡膠材料的應(yīng)力-應(yīng)變關(guān)系，采用Ogden模型作為橡膠襯套的本構(gòu)模型，來描述橡膠材料的超彈性變形特性[15]。

Ogden模型的應(yīng)變能密度函數(shù)形式為

[U=i=1N2μiαi2（λ1αi+λ2αi+λ3αi-3）+i=1N1DiJ-12i] （1）

式中：[U]為應(yīng)變能密度；[μi]和[αi]為試驗測定的材料參數(shù)；[Di]為不可壓縮參數(shù)；[λi]為伸長比；N為Ogden模型的階數(shù)，通常可以選擇為3；[J]為變形體積與未變形體積之比。

與其他本構(gòu)模型比較，Ogden本構(gòu)模型能更好的適用于橡膠材料的非線性、大變形特性，具有較好的擬合精度。因3階Ogden模型參數(shù)較少，計算和實現(xiàn)相對簡單，工程上常用3階Ogden模型用于橡膠襯套本構(gòu)模型的擬合。基于單軸拉伸、雙軸拉伸和平面剪切試驗，擬合得到材料參數(shù)如表1所示。

根據(jù)橡膠襯套的實際工況載荷，在加載點1處，施加1 mm的位移或0.1°的轉(zhuǎn)角載荷；在加載面2處，加載縮徑量大小，并同時約束6個方向的自由度；最后，提取對應(yīng)的力（力矩）與位移（轉(zhuǎn)角）的變化曲線，仿真得到襯套各向剛度。

1.2 橡膠襯套有限元模型試驗驗證

采用MTS831彈性體試驗臺（力值：±25 kN，位移：±30 mm，頻率：0～400 Hz，溫度：-100～+500 ℃，精度：0.2%）測試襯套的徑向剛度和軸向剛度，通過MTS零部件試驗臺（扭矩：±565 N·m，角度：±45°，頻率：0～100 Hz，溫度：-70～+170 ℃，精度：0.2%）測試襯套的偏擺剛度和扭轉(zhuǎn)剛度，以驗證有限元模型的準確性，試驗設(shè)備如圖2所示。

橡膠襯套通過夾具固定在試驗臺架上，激勵裝置緩慢對橡膠襯套施加載荷，通過傳感器連續(xù)記錄施加在襯套上的載荷和橡膠襯套變形量。預(yù)加載1次，取第2次加載過程的力（力矩）和位移（轉(zhuǎn)角）曲線，計算出對應(yīng)方向的剛度。

對橡膠襯套模型進行有限元仿真分析，并將橡膠襯套進行試驗測定，仿真結(jié)果和試驗結(jié)果對比如圖3所示。

通過比較仿真曲線和試驗曲線，并對數(shù)據(jù)進行擬合處理，初始襯套徑向剛度[Kx]、軸向剛度[Kz]、偏擺剛度[Knx]和扭轉(zhuǎn)剛度[Knz]的有限元仿真與試驗數(shù)據(jù)對比如表2所示。

由表2可以看出，試驗結(jié)果稍大于仿真結(jié)果，相比于有限元仿真結(jié)果，徑向剛度、軸向剛度、偏擺剛度和扭轉(zhuǎn)剛度的相對誤差分別為5.12%，4.69%、

-7.27%和3.60%，均在±10%以內(nèi)，在后續(xù)的各向剛度優(yōu)化處理中，可采用橡膠襯套有限元模型進行分析；產(chǎn)生誤差的主要原因是在仿真分析時簡化了橡膠襯套模型，忽略了金屬部件（內(nèi)芯和外套）的變形。

2 橡膠襯套各向剛度響應(yīng)面建模及分析

2.1 響應(yīng)面模型建立

考慮響應(yīng)變量和自變量的線性和非線性關(guān)系以及自變量間的交互作用，通過引入響應(yīng)面方法，擬合各向剛度值與參數(shù)變量的關(guān)系?；谇捌谘芯?，相比于一階響應(yīng)面模型，二階響應(yīng)面模型擬合效果會更好，故本文采用二階響應(yīng)面模型。所建立的二階響應(yīng)面模型為

[Y=β0+i=1nβixi+i=1nβiixi2+?ij（i

式中：[Y]為響應(yīng)變量；[β0]和[βi]為多項式系數(shù)；[xi]和[xj]為設(shè)計變量；[n]為設(shè)計變量的個數(shù)。

考慮縮徑工藝的要求，橡膠襯套的縮徑量一般應(yīng)在5%以下，初始橡膠襯套的橡膠緩沖厚度為8 mm，基于上下波動原則取縮徑量分別為0.2，0.3 mm和0.4 mm。不同材料硬度對橡膠襯套的本構(gòu)模型和剛度影響較大，當材料硬度大于68 SHA時，隔振器的靜、動剛度值隨硬度的變化速率超過了50%，橡膠隔振器已失去隔振作用[16]。對于新能源汽車、特種車輛等載運工具，工程上常使用50，55 SHA和60 SHA的車用橡膠襯套。

通過試驗設(shè)計研究因素對橡膠襯套徑向剛度、軸向剛度、偏擺剛度和扭轉(zhuǎn)剛度的影響程度，按照上下波動原則分別對4個因素均勻取3個水平，初始值及區(qū)間取值如表3所示。

由表3可知，試驗設(shè)計為4因素3水平，選取試驗設(shè)計中心點數(shù)為3，確定了27組試驗[17]，試驗設(shè)計結(jié)果如表4所示。

根據(jù)式（2），建立襯套各向剛度的二階響應(yīng)面模型為

[ki=γ1+γ2a+γ3b+γ4c+γ5d+γ6ab+?????? γ7ac+γ8ad+γ9bc+γ10bd+γ11cd+?????? γ12a2+γ13b2+γ14c2+γ15d2] （3）

式中：[ki]為各向剛度的數(shù)值（[i=x，y，nx，nz]），不含單位；[γj]為系數(shù)（[j=1，2，…，15]）。上式應(yīng)為無量綱計算公式，即等式左右兩邊無單位。

將試驗設(shè)計中的27個樣本點代入式（3）中，通過最小二乘法擬合各向剛度二次多項式響應(yīng)面近似函數(shù)系數(shù)，式（3）中系數(shù)見表5。

2.2 響應(yīng)面模型精度分析

采用拉丁超立方抽樣方法，隨機選取20個樣本點對響應(yīng)面近似模型的擬合程度進行驗證[18]。通過樣本點的響應(yīng)值和預(yù)測值之間的相關(guān)程度來評價響應(yīng)面模型的精度，統(tǒng)計學中常用決定系數(shù)[R2]來表征響應(yīng)面模型的逼近程度。決定系數(shù)的計算公式為

[R2=1-QCQZ] （4）

[QC=i=1n（yi-ypi）2] （5）

[QZ=i=1n（yi-y）2] （6）

[y=1ni=1nyi] （7）

式中：[QC]為殘差平方和；[QZ]為偏差平方和；[n]為樣本點數(shù)目；[yi]和[ypi]分別為第[i]個樣本點的響應(yīng)值和預(yù)測值；[y]為所有樣本點的響應(yīng)值的均值。當[R2]越接近于1，表示響應(yīng)面擬合效果越好。表6為襯套各個方向的決定系數(shù)。

從表6可見，各向剛度的響應(yīng)面模型的決定系數(shù)均大于0.9且接近于1，響應(yīng)面模型擬合程度較高，可用于多目標優(yōu)化。

對表4數(shù)據(jù)進行因素分析，得到各個變量對橡膠襯套各向剛度的影響程度和貢獻率，分析結(jié)果見表7。表7中的[F]值為統(tǒng)計檢定值，[F]值越大，說明變量對剛度結(jié)果影響越大

[F=R2/（k-1）（1-R2）/（n-k）] （8）

式中：[F]為各向剛度統(tǒng)計檢定值；[n]為全部觀測值個數(shù)；[k]為因素水平個數(shù)。

分析表7的數(shù)據(jù)可知：內(nèi)芯外徑[a]對徑向剛度的影響最顯著，貢獻率為93.21%；內(nèi)芯外徑[a]和材料硬度[d]對軸向剛度的影響較大，貢獻率分別為67.93%和27.49%；內(nèi)芯外徑[a]對偏擺剛度的影響最顯著，貢獻率為83.79%；內(nèi)芯外徑[a]和材料硬度[d]對扭轉(zhuǎn)剛度的影響較大，貢獻率分別為80.43%和18.30%。

為了分析材料硬度和縮徑量及其交互作用對橡膠襯套各向剛度的影響，由響應(yīng)面模型繪制的部分響應(yīng)曲線如圖4所示。

由圖4分析可知，材料硬度與縮徑量對橡膠襯套徑向剛度、軸向剛度、偏擺剛度和扭轉(zhuǎn)剛度的交互作用不明顯，且在一定范圍內(nèi)，橡膠襯套徑向剛度、軸向剛度、偏擺剛度和扭轉(zhuǎn)剛度隨著材料硬度和縮徑量的增大而增大。

3 橡膠襯套各向剛度優(yōu)化設(shè)計

3.1 目標函數(shù)設(shè)計

基于建立的響應(yīng)面模型，以內(nèi)芯外徑a、橡膠高度b、縮徑量c和材料硬度d為設(shè)計變量，以橡膠襯套的徑向剛度kx、軸向剛度kz、偏擺剛度knx和扭轉(zhuǎn)剛度knz為目標，進行多目標優(yōu)化設(shè)計。

目標剛度kx取值為8 000 N/mm，kz取值為320 N/mm，knx取值為6 N·m/（°），knz取值為1.2 N·m/（°）。利用均值化法對目標函數(shù)進行無量綱化，建立目標函數(shù)如式（9）。目標函數(shù)值越小，各向剛度與目標剛度越接近

[minf=i=14Fiai?] （9）

[F1=kx-8 000F2=kz-320F3=knx-6F4=knz-1.2] （10）

[ai=j=1mkxj-kx0m] （11）

式中：[Fi]為各向剛度與目標剛度的差值；[kxj]為各向剛度；[kx0]為目標剛度；[m]為試驗設(shè)計的數(shù)目，m=27；[ai]為無量綱化后的分配系數(shù)，即試驗設(shè)計中各向剛度與目標剛度差值的絕對值求和后的平均值，具體取值見表8。

約束條件設(shè)置為：[21 mm≤a≤25 mm;25 mm≤][b≤27 mm;0.2 mm≤c≤0.4 mm;50 SHA≤d≤60 SHA。]

3.2 優(yōu)化結(jié)果分析

NSGA-Ⅱ遺傳算法引入擁擠度計算和精英策略，避免優(yōu)化過程中產(chǎn)生局部收斂，具有較高的計算效率和魯棒性，被廣泛應(yīng)用于多目標優(yōu)化問題[19]。本文采用MATLAB軟件運行NSGA-Ⅱ遺傳算法，并對橡膠襯套各向剛度進行多目標優(yōu)化。設(shè)置初始化種群數(shù)為50、最大迭代次數(shù)為200、交叉比例為0.8、變異概率為0.05，步長分別設(shè)置為0.01，0.01，0.01，5。遺傳迭代運算后，得到目標函數(shù)值及最優(yōu)參數(shù)，優(yōu)化結(jié)果如表9所示。

通過對比表9中的非劣解集，第9組數(shù)據(jù)的目標函數(shù)值最小，對第9組數(shù)據(jù)中的參數(shù)取整，取整后的最優(yōu)參數(shù)為內(nèi)芯外徑25 mm、橡膠高度26 mm、縮徑量0.4 mm，材料硬度50 SHA。根據(jù)取整后的最優(yōu)解參數(shù)設(shè)計橡膠襯套成品并進行剛度測試試驗，對比各向剛度和目標剛度，如圖5所示。

由圖5可見，橡膠襯套各向剛度的試驗值與目標剛度接近。通過二階響應(yīng)面模型計算出的結(jié)果、有限元仿真的分析結(jié)果以及試驗結(jié)果分別與目標值進行對比，見表10。

分析表10中的數(shù)據(jù)可知，根據(jù)優(yōu)化結(jié)果對橡膠襯套參數(shù)取整后分析，代入響應(yīng)面模型擬合后的軸向剛度kx，kz分別為8 049.40，319.03 N/mm，扭轉(zhuǎn)剛度knx，knz分別為6.29，1.34 N·m/（°），與目標值的誤差分別為0.62%，-0.3%，4.83%和11.67%；有限元仿真后的軸向剛度kx，kz分別為8 184.00，333.56 N/mm，扭轉(zhuǎn)剛度knx，knz分別為6.11，1.27 N·m/（°），與目標值的誤差分別為2.3%，4.24%，1.83%和5.83%；試驗測定后的軸向剛度kx，kz分別為8 617.96，348.98 N/mm，扭轉(zhuǎn)剛度knx，knz分別為5.62，1.31 N·m/（°），與目標值的誤差分別為7.72%，9.06%，-6.33%和9.16%，剛度值誤差均在±10%以內(nèi)，符合工程應(yīng)用的要求。

4 結(jié)論

本文在橡膠襯套各向剛度優(yōu)化設(shè)計中考慮了材料硬度和縮徑量的影響，基于試驗設(shè)計與有限元分析相結(jié)合的方法，建立橡膠襯套各向剛度的二階響應(yīng)面模型，采用NSGA-Ⅱ遺傳算法對各向剛度進行多目標優(yōu)化研究，得出以下結(jié)論：

1） 建立關(guān)于內(nèi)芯外徑、橡膠高度、縮徑量和材料硬度隨襯套各向剛度變化的二階響應(yīng)面模型，并驗證了模型的有效性，可為載運工具橡膠襯套各向剛度優(yōu)化設(shè)計提供必要條件。

2） 縮徑量對橡膠襯套軸向剛度和偏擺剛度影響較大，貢獻率分別為42.9%和50.5%，通過增大縮徑量，可有效提高軸向剛度，提高產(chǎn)品的耐久性；材料硬度對襯套軸向剛度和扭轉(zhuǎn)剛度影響較大，貢獻率分別為45.6%和30.4%，通過增大材料硬度，可有效提高扭轉(zhuǎn)剛度，提高整車操縱穩(wěn)定性。

3） 采用NSGA-Ⅱ遺傳算法對橡膠襯套各向剛度進行多目標優(yōu)化，并將優(yōu)化結(jié)果進行試驗驗證，各向剛度試驗結(jié)果與目標值的相對誤差分別為7.72%，9.06%，-6.33%和9.16%，均在±10%以內(nèi)，滿足工程應(yīng)用的要求，可有效的避免工程實踐中反復(fù)調(diào)試的繁瑣工作，從而極大地縮短了產(chǎn)品研發(fā)周期。

參考文獻：

[1]?? LI L， SUN B HUA H. Analysis of radial stiffness of rubber bush used in dynamic vibration absorber[J]. Journal of Southeast University， 2019， 35（3）： 281-287.

[2]?? 方紅， 李晨芳， 辛振祥. 縮徑尺寸對橡膠襯套剛度耐久性的影響[J]. 中國機械， 2020： 183-184.

FANG H， LI C F， XIN Z X. Effects of reduction size on the stiffness and durability of rubber bushings[J]. China Machinery， 2020： 183-184.

[3]?? 陳寶， 代正莉， 付江華， 等. 汽車懸架襯套靜動特性分析中縮徑量影響研究[J]. 噪聲與振動控制， 2019， 39（1）： 126-130.

CHEN B， DAI Z L， FU J H， et al. Study on the influence of shrinkage in the static-dynamic characteristic analysis of automotive suspension bushings[J]. Noise and Vibration Control， 2019， 39（1）： 126-130.

[4]?? 曾家興， 丁智平， 林勝， 等. 膠料硬度對超彈本構(gòu)模型參數(shù)及橡膠彈簧剛度的影響[J]. 湖南工業(yè)大學學報， 2018，32（3）： 92-96.

ZENG J X， DING Z P， LIN S， et al. Effect of rubber hardness on parameters of superelastic intrinsic model and rubber spring stiffness[J]. Journal of Hunan University of Technology， 2018， 32（3）： 92-96.

[5]?? 陳寶， 張鑫， 代正莉， 等. 硬度對橡膠力學特性與懸架襯套剛度的影響[J]. 噪聲與振動控制， 2020， 40（6）： 222-227.

CHEN B， ZHANG X， DAI Z L， et al. Effect of hardness on mechanical properties of rubber and suspension bushing stiffness[J]. Noise and Vibration Control， 2020， 40（6）：222-227.

[6]?? KAYA N， ERKEK M Y， GUVEN C. Hyperelastic modelling and shape optimisation of vehicle rubber bushings[J]. International Journal of Vehicle Design， 2016， 71（1/2/3/4）： 212-225.

[7]?? DONG Y， YAO X， XU X. Cross section shape optimization design of fabric rubber seal[J]. Composite Structures， 2021， 256： 113047.

[8]?? BILAL H， OZTURK F. Rubber bushing optimization by using a novel chaotic krill herd optimization algorithm.[J]. Soft Comput， 2021， 25， 14333-14355.

[9]?? 丁智平， 曾家興， 林勝， 等. 橡膠彈簧結(jié)構(gòu)強度分析及參數(shù)優(yōu)化[J]. 振動與沖擊， 2019， 38（12）： 246-251.

DING Z P， ZENG J X， LIN S， et al. Structural strength analysis and parameter optimization of rubber springs[J]. Vibration and Shock， 2019， 38（12）：246-251.

[10] DAI L， CHI M， XU C， et al. A hybrid neural network model based modelling methodology for the rubber bushing[J]. Vehicle System Dynamics， 2022， 60（9）： 2941-2962.

[11] JUNG Y W， KIM H K. Prediction of nonlinear stiffness of automotive bushings by artificial neural network models trained by data from finite element analysis[J]. Automot.Technol， 2020， 21： 1539-1551.

[12] 柯俊， 祖洪飛， 史文庫. 基于有限元法及遺傳算法的推力桿球鉸多目標優(yōu)化方法[J]. 汽車工程， 2020， 42（2）： 178-183.

KE J， ZU H F， SHI W K. Multi-objective optimization method of thrust rod ball hinge based on finite element method and genetic algorithm[J]. Automotive Engineering， 2020， 42（2）： 178-183.

[13] 郭一鳴， 吳鐘愷， 上官文斌. 襯套型橡膠懸置結(jié)構(gòu)形狀的優(yōu)化方法[J]. 振動與沖擊， 2022， 41（5）： 55-59.

GUO Y M， WU Z K， SHANGGUAN W B. Optimization method for the shape of bushing-type rubber suspension structure[J]. Vibration and Shock， 2022， 41（5）： 55-59.

[14] 鄧雄志， 劉濤， 蘇澤博. 轎車底盤橡膠襯套半經(jīng)驗設(shè)計方法研究[J]. 汽車技術(shù)， 2018（7）： 54-58.

DENG X Z， LIU T， SU Z B. Research on semi-empirical design method of rubber bushing in car chassis[J]. Automotive Technology， 2018（7）： 54-58.

[15] 史文庫， 劉國政， 陳志勇. 橡膠的可壓縮性在推力桿球鉸有限元分析中的應(yīng)用[J]. 西安交通大學學報， 2017， 51（9）： 63-68.

SHI W K， LIU G Z， CHEN Z Y. Application of rubber compressibility in the finite element analysis of thrust rod ball hinge[J]. Journal of Xian Jiaotong University， 2017， 51（9）： 63-68.

[16] 劉山尖， 孫偉星， 董興建. 基于有限元模型的橡膠減振元件性能研究及失效分析[J]. 機電設(shè)備， 2014， 31（5）： 41-44.

LIU S J， SUN W X， DONG X J. Performance study and failure analysis of rubber vibration damping element based on finite element model[J]. Electromechanical Equipment， 2014， 31（5）： 41-44.

[17] KOSARAC A， MLADJENOVIC C， ZELJKOVIC M， et al. Neural-network-based approaches for optimization of machining parameters using small dataset[J]. Materials， 2022， 15（3）： 700.

[18] 陳靜， 徐森， 劉震， 等. 基于碰撞安全性的鋁合金吸能盒輕量優(yōu)化[J]. 汽車工程， 2021， 43（2）： 241-247.

CHEN J， XU S， LIU Z， et al. Lightweight optimization of aluminium alloy energy absorption box based on crash safety[J]. Automotive Engineering， 2021， 43（2）： 241-247.

[19] ZHANG Z， LIAO C， CHAI H， et al. Multi-objective optimization of controllable configurations for bistable laminates using NSGA-Ⅱ[J]. Composite Structures， 2021， 266： 113764.

第一作者：秦武（1990—），男，講師，博士，研究方向為橡膠隔振件性能分析及動力總成懸置系統(tǒng)設(shè)計。E-mail： qw@ ecjtu.edu.cn

通信作者：李春歸（1999—），男，碩士研究生，研究方向為橡膠襯套性能分析。E-mail： 3047855129@qq.com