雷 剛，陳華良，賀艷輝，王 靜

(重慶理工大學(xué) a.汽車零部件先進(jìn)制造技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室;b.車輛工程學(xué)院, 重慶 400054)

?

提取懸架連接點(diǎn)載荷的等效模型提載法

雷 剛a，陳華良b，賀艷輝b，王 靜b

(重慶理工大學(xué) a.汽車零部件先進(jìn)制造技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室;b.車輛工程學(xué)院, 重慶 400054)

提出了一種基于有限元分析軟件Abaqus的懸架連接點(diǎn)載荷提取方法，即等效模型提載法。以四連桿懸架為例，建立了基于 Abaqus的懸架完整模型及等效模型，分析了建立等效模型的依據(jù)，建立了基于多體動(dòng)力學(xué)分析軟件Adams的懸架模型。分別提取出各模型在典型工況下的硬點(diǎn)載荷進(jìn)行對比，發(fā)現(xiàn)等效模型與完整模型的載荷值十分接近，驗(yàn)證了懸架等效模型提載法的準(zhǔn)確性，并簡要討論了各模型之間載荷值差異的原因。該懸架等效模型提載法簡單快速，對縮短懸架的開發(fā)周期具有重要意義。

四連桿懸架；硬點(diǎn)；載荷提?。坏刃Ｐ?/p>

懸架是車架或車身與車輪之間所有傳力連接部件的總稱，它將路面作用于車輪上的垂直反力、縱向反力和側(cè)向反力以及這些反力所形成的力矩傳遞到車架或車身上，以保證汽車正常行駛[1]。為保證汽車行駛的安全性和操縱穩(wěn)定性，懸架系統(tǒng)各部件應(yīng)具有足夠的強(qiáng)度、剛度和使用壽命，所以在汽車設(shè)計(jì)初期，對懸架系統(tǒng)中的部件進(jìn)行強(qiáng)度、剛度以及耐久性分析是不可或缺的[2-4]。在相應(yīng)的分析中，懸架系統(tǒng)硬點(diǎn)處載荷值是必不可少的輸入條件。然而，由于懸架系統(tǒng)中部件較多且裝配復(fù)雜，通過試驗(yàn)的方法獲取懸架系統(tǒng)硬點(diǎn)處載荷顯得異常困難。目前提取懸架系統(tǒng)硬點(diǎn)載荷多在多體動(dòng)力學(xué)軟件Adams中完成，但該過程比較繁瑣，延長了產(chǎn)品的開發(fā)周期。

本文提出了一種基于有限元分析軟件Abaqus的懸架等效模型提載方法。以四連桿懸架為例，建立了基于 Abaqus的懸架完整模型及等效模型，以及基于多體動(dòng)力學(xué)分析軟件Adams的懸架模型。分別提取出各模型在典型工況下的硬點(diǎn)載荷進(jìn)行對比，發(fā)現(xiàn)等效模型與完整模型的載荷值十分接近。該方法簡單快速，可縮短懸架的開發(fā)周期。

四連桿式懸架是技術(shù)比較先進(jìn)的懸架形式，可以從設(shè)計(jì)上保證車輛良好的直線行駛性能，并最大程度減輕載荷的影響，在中高端汽車上有較為廣泛的應(yīng)用[5]。四連桿懸架的典型結(jié)構(gòu)如圖1所示。

1 四連桿懸架有限元模型的建立

1.1 網(wǎng)格劃分與連接搭建

將四連桿后懸架系統(tǒng)的數(shù)學(xué)模型以Step 格式導(dǎo)人HyperWorks軟件Hypermesh模塊中進(jìn)行網(wǎng)格劃分。后軸節(jié)為鑄造件，采用四面體單元進(jìn)行網(wǎng)格劃分。由于后上控制臂厚度為9 mm，因此采用六面體單元進(jìn)行網(wǎng)格劃分。穩(wěn)定桿采用六面體單元與四面體單元相結(jié)合的方法進(jìn)行網(wǎng)格劃分。穩(wěn)定桿連桿采用六面體單元和殼單元進(jìn)行網(wǎng)格劃分。其他部件均采用殼單元進(jìn)行網(wǎng)格劃分，平均尺寸為5，最小尺寸為2。經(jīng)網(wǎng)格劃分后共得到139 704個(gè)單元和91 559個(gè)節(jié)點(diǎn)。

圖1 四連桿后懸架的典型結(jié)構(gòu)

在Hypermesh Abaqus模板下進(jìn)行前處理。焊縫用四邊形殼單元表示，并賦予焊縫屬性，這樣能較好地表現(xiàn)出焊縫及附近的應(yīng)力值。螺栓連接用rigid剛性單元表示，襯套用bushing單元表示。bushing單元實(shí)際是六向彈簧，考慮到襯套剛度對懸架連接點(diǎn)載荷的影響很大，本文所有襯套6個(gè)方向上的剛度皆為非線性。后螺旋彈簧用AXIAL單元表示，減震器用cylindrical單元表示，減震器彈簧用AXIAL+ALIGN單元表示[6]。輪胎接地點(diǎn)與輪心采用剛性連接，輪心與后軸節(jié)采用剛性連接。完整的四連桿懸架有限元模型如圖2所示。

1.2 分析步與約束條件

建立兩個(gè)分析步。在分析步1，后懸架縱臂與車身連接處全約束，螺旋彈簧與車身連接處全約束，副車架與車身連接處全約束，減震器與車身連接處全約束，輪胎接地點(diǎn)全約束。在分析步2，釋放輪胎接地點(diǎn)處的約束并在輪胎接地點(diǎn)處施加工況載荷。輸出了應(yīng)力、應(yīng)變、位移、連接點(diǎn)載荷等數(shù)據(jù)[7]。轉(zhuǎn)彎制動(dòng)工況下的應(yīng)力如圖3所示，工況載荷見表1。

圖2 四連桿后懸架完整有限元模型

圖3 轉(zhuǎn)彎制動(dòng)工況下的應(yīng)力 表1 工況載荷

工況方向輪胎接地點(diǎn)載荷 左 右垂直沖擊工況Fz/N1897018970前行制動(dòng)Fz/N27542754Fx/N27542754倒車制動(dòng)Fz/N83498349Fx/N-8349-8349轉(zhuǎn)彎工況(右轉(zhuǎn))Fz/N10089985Fy/N10089985Tz/(N·mm)29556428846轉(zhuǎn)彎制動(dòng)工況(右轉(zhuǎn))Fz/N56431456Fx/N39211026Fy/N39211026Tz/(N·mm)11896531259

2 四連桿懸架等效模型的建立

影響四連桿后懸架系統(tǒng)連接點(diǎn)載荷的主要因素有硬點(diǎn)坐標(biāo)、彈性元件(襯套和彈簧)的剛度參數(shù)和各桿件的剛度。在硬點(diǎn)坐標(biāo)和彈性元件剛度不變的情況下，影響四連桿后懸架系統(tǒng)連接點(diǎn)載荷的主要因素是各桿件的剛度。在較小的位移下襯套剛度遠(yuǎn)小于各桿件的剛度，而在典型工況下襯套位移比較小，故在典型工況下各桿件的剛度對懸架系統(tǒng)連接點(diǎn)載荷的影響很小，本文基于此建立懸架等效模型。

保持硬點(diǎn)坐標(biāo)和彈性元件的剛度參數(shù)與完整模型一致，懸架的上控制臂、下控制臂、拖曳臂、拉桿均用剛性梁單元(beam)代替。取消后軸節(jié)，直接用剛性梁單元連接輪心至后軸節(jié)上的各硬點(diǎn)。取消副車架，將約束直接加在懸架各桿與副車架的連接點(diǎn)上。取消后拖曳臂與車身的安裝支架，將約束直接加在后拖曳臂與安裝支架的連接點(diǎn)上。取消減震器，減震器用cylindrical單元表示，減震器彈簧用AXIAL+ALIGN單元表示，并根據(jù)具體情況增加約束條件[8]以保證其運(yùn)動(dòng)形式與完整模型一致。保留穩(wěn)定桿，這是因?yàn)榉€(wěn)定桿在含有使汽車有側(cè)傾趨勢的工況下對連接點(diǎn)載荷影響較大，且難以用其他單元模擬。分析步、約束條件、工況載荷與完整懸架模型保持一致。建立的等效模型如圖4所示，由圖中可以看出：除了穩(wěn)定桿以外懸架的其他部件皆不需要實(shí)際模型。也就是說，除了穩(wěn)定桿以外的其他部件皆不需要?jiǎng)澐志W(wǎng)格，從而可以節(jié)省大量時(shí)間。

圖4 四連桿后懸架等效模型

3 四連桿懸架多體動(dòng)力學(xué)模型的建立

本文四連桿懸架的多體動(dòng)力學(xué)模型是在Adams/car模塊中建立的。Adams/car中的約束種類可分為鉸鏈約束(joint)和橡膠襯套約束(bushing)，通過定義不同的鉸鏈連接方式，所建立的模型可以進(jìn)行運(yùn)動(dòng)學(xué)分析或彈性運(yùn)動(dòng)學(xué)分析[9]。在Adams/car模塊中建立懸架模型同樣需要知道各硬點(diǎn)坐標(biāo)，其硬點(diǎn)坐標(biāo)與前兩個(gè)模型保持一致?？紤]到拖曳臂為薄長板(長約400 mm，厚為4.5 mm)，且在制動(dòng)工況下長度方向受到較大載荷，此處將拖曳臂建為柔性體。模型中除了穩(wěn)定桿和拖曳臂為柔體外其余桿件皆為剛體，襯套參數(shù)(包括材料參數(shù)和安裝角度)、工況載荷與前面兩個(gè)模型保持一致。具體的建模過程比較繁瑣，故在此不再敘述。建立的四連桿后懸架多體動(dòng)力學(xué)模型如圖5所示。

圖5 四連桿后懸架多體動(dòng)力學(xué)模型

4 各模型提載結(jié)果的對比

4.1 參考點(diǎn)的選擇

四連桿懸架左右兩邊共有20個(gè)連接點(diǎn)，但在靜態(tài)或準(zhǔn)靜態(tài)分析時(shí)單邊4個(gè)桿兩端連接點(diǎn)載荷大小相等、方向相反，為減小工作量，在此只提取一端的載荷。一般來說，四連桿懸架的拖曳臂與后軸節(jié)都有2個(gè)及以上的連接點(diǎn)，本模型有3個(gè)，但在Adams多體動(dòng)力學(xué)模型中一般簡化成1個(gè)連接點(diǎn)。為便于提取結(jié)果的比較，此處將拖曳臂與車身的連接點(diǎn)作為參考點(diǎn)，即左右兩邊各選取5個(gè)點(diǎn)作為參考點(diǎn)，共10個(gè)參考點(diǎn)，每個(gè)參考點(diǎn)提取3個(gè)方向上的力和3個(gè)方向的力矩。參考點(diǎn)位置如圖6所示。

圖6 參考點(diǎn)位置

圖6中：參考點(diǎn)1為后軸節(jié)與后上控制臂的襯套連接點(diǎn)；參考點(diǎn)2為后軸節(jié)與減震器的襯套連接點(diǎn)；參考點(diǎn)3為后軸節(jié)與后下控制臂的襯套連接點(diǎn)；參考點(diǎn)4為后軸節(jié)與拉桿的襯套連接點(diǎn)；參考點(diǎn)5為后拖曳臂與車身的襯套連接點(diǎn)。

4.2 分析工況

本文選取5個(gè)典型的懸架分析工況來進(jìn)行載荷提取，分別為垂直沖擊工況、前行制動(dòng)工況、倒車制動(dòng)工況、轉(zhuǎn)彎工況和轉(zhuǎn)彎制動(dòng)工況。工況載荷如表1所示。本中Fx、Fy、Fz分別表示x軸、y軸、z軸上的力，Tx、Ty、Tz分別表示繞x、y、z軸的力矩。

4.3 結(jié)果對比

每個(gè)工況的每個(gè)模型都有60個(gè)載荷(10個(gè)參考點(diǎn)，每個(gè)參考點(diǎn)提取3個(gè)力和3個(gè)力矩)。若只看相同工況下模型之間載荷的最大差值很難反映出相同工況下模型之間載荷的總體差異情況。本文用距離函數(shù)來表示兩個(gè)模型之間載荷的差異程度。距離函數(shù)有很多種，歐式距離又稱歐幾里得度量(euclidean metric)，是一種廣泛使用的距離函數(shù)[10]。歐式距離是指在n維空間中兩個(gè)點(diǎn)之間的真實(shí)距離或者向量的自然長度(即該點(diǎn)到原點(diǎn)的距離)，在二維和三維空間中的歐氏距離就是兩點(diǎn)之間的實(shí)際距離。本文采用歐式距離來評價(jià)兩個(gè)模型之間載荷的總體差異程度。兩n維向量a(x11,x12,…,x1n)與b(x21,x22,…,x2n)間的歐氏距離計(jì)算公式為

(1)

其中：d12是兩n維向量a、b的歐式距離；x1k表示第1個(gè)向量的第k個(gè)分量；x2k表示第2個(gè)向量的第k個(gè)分量。

從式(1)可以看出：歐式距離對差值較大的分量比差值較小的分量敏感，這種特性正是本文所需要的。

此次分析提取了各模型各工況下的載荷，但限于篇幅只給出受力情況比較復(fù)雜的轉(zhuǎn)彎制動(dòng)工況下的載荷值。表2～4中的L1,L2,…,L5表示懸架左邊的參考點(diǎn)，R1,R2,…,R5表示懸架右邊的參考點(diǎn)，L1表示懸架左邊的參考點(diǎn)1，以此類推。

表2 完整模型轉(zhuǎn)彎制動(dòng)工況下載荷值

表3 等效模型轉(zhuǎn)彎制動(dòng)工況下載荷值

表4 多體動(dòng)力學(xué)模型轉(zhuǎn)彎制動(dòng)工況下載荷值

表5中：Dfb表示完整懸架模型載荷值與等效懸架模型載荷值的歐式距離；Dfa表示完整懸架模型載荷值與懸架多體動(dòng)力學(xué)模型載荷值歐式距離；Dba表示等效懸架模型載荷值與懸架多體動(dòng)力學(xué)模型載荷值的歐式距離。

表5 各工況下模型之間的歐氏距離

從表2～4中的載荷數(shù)據(jù)來看，3個(gè)模型在轉(zhuǎn)彎制動(dòng)工況下各參考點(diǎn)載荷的差異很小。表2、3的最大差值出現(xiàn)在(L3，Tx)，差值絕對值為303.1，最大差值差異程度為1.28%(相對于完整模型)；表2、4的最大差值出現(xiàn)在(L3，Tx)，差值絕對值為897，最大差值差異程度為3.95%(相對于完整模型)；表3、4的最大差值出現(xiàn)在(L3，Tx)，差值絕對值為1 200，最大差值差異程度為5.28%(相對于多體動(dòng)力學(xué)模型)。從表5中的數(shù)據(jù)來看，在各工況下Dfb均要比Dfa、Dba小得多，這可能是因?yàn)橛?jì)算軟件不同導(dǎo)致。Dfa、Dba在相同工況下的值都比較接近，這從某種程度上反映了完整懸架模型和等效懸架模型的載荷值差異很小。結(jié)合表1和表5可以看出：工況載荷越大，模型之間的歐式距離也越大。這可能是因?yàn)橥暾麘壹苣Ｐ透鞑考际侨嵝缘?，而等效懸架模型和懸架多體動(dòng)力學(xué)模型大部分部件是剛性的，工況載荷越大，剛性和柔性之間的差異程度也越大?？偟膩碚f，在各工況下完整懸架模型和等效懸架模型的載荷值是十分接近的，用等效懸架模型來提取接觸點(diǎn)處載荷的方法是可行的。

5 結(jié)束語

本文提出了一種基于有限元分析軟件Abaqus的懸架連接點(diǎn)載荷提取方法，即懸架等效模型提載法。以四連桿懸架為例，建立了基于 Abaqus的懸架完整模型及等效模型，以及基于多體動(dòng)力學(xué)分析軟件Adams的懸架模型。分別提取出各模型在典型工況下的硬點(diǎn)載荷進(jìn)行對比，發(fā)現(xiàn)等效模型與完整模型的載荷值十分接近，但差值隨著工況載荷的增大而增大，分析得出是因?yàn)榈刃Ｐ蛯⒋蟛糠謼U件等效成了剛性體導(dǎo)致的。在今后的工作中考慮多用柔性單元來代替懸架各桿件以改進(jìn)懸架等效模型，但目前還沒有找到可行的柔性替代方法。

懸架等效模型提載法簡單快速，提取的載荷可為懸架系統(tǒng)各部件以及車架(副車架)的強(qiáng)度和疲勞分析等提供輸入條件。采用此方法可以減輕汽車設(shè)計(jì)人員的工作強(qiáng)度，縮短懸架系統(tǒng)的開發(fā)周期。

[1] 雷剛，張澤俊.麥弗遜前懸架系統(tǒng)硬點(diǎn)載荷提取技術(shù)[J].重慶理工大學(xué)學(xué)報(bào)(自然科學(xué))，2016，30(2)：17-23.

[2] 上官文斌，蔣翠翠，潘孝勇．汽車懸架控制臂的拓?fù)鋬?yōu)化與性能計(jì)算[J].汽車工程，2008，30(8):709-712.

[3] MIGUEL A E，HUANG M，TYAN T．Impact Testing of Lower Control Arm for Crashworthiness Simulation Simu-lation[C]//SAE Paper 2005-01-0352．USA:[s．n．]，2005.

[4] HUANG M，MIGUEL A E，TYAN T．Crashworthiness Simulation of Lower Control Arm Impact Tests[C]//SAE Paper 2005-01-0361.USA:[s.n.],2005.

[5] 鄭松林，顧晗.一種新的四連桿懸架剛?cè)狁詈夏Ｐ蚚J].上海汽車，2011(2):42-45.

[6] 王鈺棟，金磊，洪清泉，等.HyperMesh & Hyper View應(yīng)用技巧與高級(jí)實(shí)例[M].北京：機(jī)械工業(yè)出版社,2012.

[7] CHEN Song，LEI Gang，LIU Ying.Strength Analysis and Optimization of a Torsion Beam Rear Suspension[J].Hydromechatronics Engineering,2013,41(18):45-49.

[8] 莊茁，由小川，廖劍暉， 等.基于ABAQUS的有限元分析和應(yīng)用[M].北京：清華大學(xué)出版社，2009.

[9] 海辰光.雙連桿后懸架模板的開發(fā)及仿真分析[D].長春:吉林大學(xué),2007.

[10]孔超，張化，劉麗.基于興趣區(qū)域特征融合的半監(jiān)督圖像檢索算法[J].山東大學(xué)學(xué)報(bào)，2014,44(3):22-28.

(責(zé)任編輯 楊黎麗)

Extracting Suspension Connecting Point Load with Equivalent Model Method

LEI Ganga, CHEN Hua-liangb, HE Yan-huib, WANG Jingb

(a.The Key Laboratory of Automobile Parts & Test Technique, Ministry of Education;b.College of Vehicle Engineering, Chongqing University of Technology, Chongqing 400054, China)

This paper presented a method for the load extraction of the suspension connecting point based on the finite element analysis software Abaqus, the equivalent model method. Take the four-link rods type suspension as an example, a complete model and an equivalent model of suspension based on Abaqus were established, and the basis of establishing the equivalent model was analyzed, and a suspension model based on multi-body dynamics analysis software Adams was established. The hard point load was compared with each model in the typical working condition.And it finds that the load value of equivalent model is very close to the full model. This verified the accuracy of the equivalent model method. The causes of load difference between the models are discussed briefly as also. The equivalent model method is simple and fast, and it has important significance to shorten the development period of the suspension.

four-link rods type suspension；hard point；load extraction；equivalent model

2016-12-25

重慶市基礎(chǔ)與前沿研究計(jì)劃資助項(xiàng)目(cstc2013jcyjA60005)

雷剛(1967—)，男,博士,教授,主要從事計(jì)算力學(xué)、CAD/CAE研究，E-mail:ganglei4786@126.com；陳華良(1992—)，男，碩士研究生，主要從事 CAD/CAE 研究，E-mail:554036674@qq.com。

雷剛，陳華良，賀艷輝，等.提取懸架連接點(diǎn)載荷的等效模型提載法[J].重慶理工大學(xué)學(xué)報(bào)(自然科學(xué))，2017(4):21-26.

format：LEI Gang, CHEN Hua-liang, HE Yan-hui,et al.Extracting Suspension Connecting Point Load with Equivalent Model Method[J].Journal of Chongqing University of Technology(Natural Science)，2017(4):21-26.

10.3969/j.issn.1674-8425(z).2017.04.004

U463.33

A

1674-8425(2017)04-0021-06