單 萍，鄭忠才，李 飛，李 達

（山東建筑大學(xué)，濟南 250101）

0 引言

2016年，我國頒布的《節(jié)能與新能源汽車技術(shù)路線圖》指出，到2020年，乘用車新車平均油耗要求達到 5.0L/100km；到2025年，乘用車新車平均油耗達到4.0 L/100km；到2030年，乘用車新車平均油耗達到3.2 L/100km[2]。而歐洲鋁協(xié)研究數(shù)據(jù)表明：汽車質(zhì)量每降低100kg，每百公里可節(jié)約0.6L燃油[3]并且每千克鋁合金的使用可以減少13～20kg溫室氣體的排放[4]，針對如今日益嚴(yán)格的排放法規(guī)與日益惡劣的環(huán)境問題，實現(xiàn)整車輕量化勢在必行。

在一臺完整的車上，有著許許多多的零部件，而輪轂是汽車上必不可少的零件之一，它承受的不止僅僅是靜態(tài)時車輛本身的自重，而且在車輛運行中需要承受交變載荷，同時在啟動、制動、加速、減速過程中，車輪需承受扭轉(zhuǎn)力矩[5]，在眾多材料的輪轂中，鋁合金輪轂因為其自身所特有的優(yōu)勢：質(zhì)量輕、散熱性好等優(yōu)點[6]，一直被人們所關(guān)注。本文針對鋁制輪轂進行了拓撲優(yōu)化設(shè)計，以達到輕量化目的，并以此結(jié)果為基礎(chǔ)，對其在三維軟件中進行修改建模。同時為了保障優(yōu)化后的強度問題，在優(yōu)化修改后再一次進行了彎曲疲勞測試，并觀察其最大應(yīng)力和最大位移的變化情況。

1 鋁合金輪轂建模

1.1 鋁合金輪轂三維實體造型

本文以較為常用的14×5.5J鋁合金輪轂做為優(yōu)化對象，本次優(yōu)化采用的三維實體建模軟件為Solidworks軟件，為了在后續(xù)的分析中縮短計算時間與儲存空間，在此選擇對螺栓進行簡化處理，其三維實體造型如圖1所示。

圖1 輪轂三維實體模型

1.2 鋁合金輪轂有限元實體模型

依據(jù)已知實體尺寸，在SolidWorks中建立實體模型，并將其保存成有限元可讀取的x_t格式，之后將其導(dǎo)入進有限元軟件中，進行后續(xù)分析。由于主要研究對象為輪轂，為了減少運算時間，決定對輪轂和加載軸實現(xiàn)不同的網(wǎng)格劃分，并在連接處進行細化處理，對輪轂選用網(wǎng)格尺寸7mm的四面體單元，對加載軸選用自由網(wǎng)格進行劃分，共計總的節(jié)點個數(shù)為153486個，總單元數(shù)為94008個，本次所劃的85%網(wǎng)格單元質(zhì)量都在0.5以上，網(wǎng)格質(zhì)量情況良好，具備進行后續(xù)優(yōu)化設(shè)計的前提。劃分網(wǎng)格之后的模型如圖2所示。

圖2 輪轂有限元模型

1.3 鋁合金輪轂材料參數(shù)

表1 14×5.5J鋁合金輪轂材料參數(shù)

2 彎曲疲勞試驗

2.1 彎曲疲勞試驗工況分析

在進行彎曲疲勞實驗時，車輪可自由選擇測試工況，其一為車輪需要在固定不變的彎矩下進行旋轉(zhuǎn)，其二為車輪需保持靜止不動，其本身需要承受一個旋轉(zhuǎn)彎曲力矩[7,8]，具體實驗裝置為圖3。

圖3 彎曲疲勞實驗裝置

2.2 載荷計算及約束條件

1）載荷計算

在載荷分析時，選擇以彎曲疲勞試驗為模擬依據(jù)，因為其可以較好地反應(yīng)出輪轂在實際工況中的受力情況。試驗彎矩估算依據(jù)為公式(1)：

上面式(1)中：

M為加載的彎矩；

μ為輪胎與路面之間的摩擦系數(shù)設(shè)為0.7；

R為靜載半徑；

d為車輪的內(nèi)偏距或外偏距，取絕對值；

F為車輪的最大額定載荷設(shè)定為1.6；

S為實驗強化系數(shù)。

結(jié)合該車輪的實際參數(shù)，計算出所需加載彎矩為2613.3N.m，根據(jù)式(2)計算出需施加在加載軸上的力。

本次試驗中選取的加載軸的長度為0.5m,據(jù)此計算出F為5226.6N，為了模擬輪轂在實際中的應(yīng)用工況，將F分解為X方向與Y方向，使FX按照正弦規(guī)律變換，使Fy按照余弦規(guī)律進行變換，并同時進行加載，保證了合力的大小不發(fā)生變化，只變化方向從而可以實現(xiàn)加載旋轉(zhuǎn)載荷模擬實際工況。

其中t為時間。其變化規(guī)律可繪制如圖3所示。

圖4 彎曲疲勞模擬所施加載荷

本文將所施加載荷進行360度連續(xù)旋轉(zhuǎn)離散為在40個方向上進行分別的加載，即使力實現(xiàn)間隔9度進行多次加載，共實現(xiàn)加載40次，將式(3)中的t分別取整數(shù)為1至40，由于本次主要完成的為靜力學(xué)分析，所以無須考慮慣性力的影響。

2）約束條件

在分析時為了簡化，忽略了螺栓及其預(yù)緊力，并定義了加載軸與輪轂之間的接觸面，根據(jù)《GBT5334-2005乘用車車輪性能要求和試驗方法》：在進行彎曲疲勞試驗中車輪需要被固定在實驗臺上，以模擬車輪在行駛過程中的實際工況[9,10]，所以在模擬疲勞試驗時，約束條件為對輪轂的下邊緣6個自由度進行全約束設(shè)置。

3 鋁合金輪轂有限元模擬分析結(jié)果

車輪處于旋轉(zhuǎn)彎曲載荷工作狀態(tài)時，其所呈現(xiàn)的是非常復(fù)雜的應(yīng)力狀態(tài)，在充分研究各軸的綜合應(yīng)力的情況下，選擇采用形狀改變比能理論，即以Von Mises屈服準(zhǔn)則作為此次分析的依據(jù)，發(fā)現(xiàn)應(yīng)力在空間分布上顯示為呈中心對稱趨勢；在時間上，呈現(xiàn)出周期性變化具體如圖4所示。由于呈現(xiàn)周期性變化，在此我們僅研究半個加載周期[11]，可以看出應(yīng)力的最大值出現(xiàn)在第10載荷步，其大小為156Mpa。

圖5 應(yīng)力變化情況

基于此，為了保證安全性，選取最大應(yīng)力值載荷步進行分析，提取該載荷步下的等效應(yīng)力云圖5和等效位移云圖6進行觀察。

圖6 輪轂結(jié)構(gòu)等效應(yīng)力云圖

圖7 輪轂結(jié)構(gòu)等效位移云圖

可從上圖看出最危險的地方在節(jié)點1579處即螺栓與加載軸連接處附近，符合實際應(yīng)用中的工況，在此處應(yīng)力值達到156MPa，但是仍然遠遠小于材料強度極限，其變形量為0.18015mm，其強度與變形量仍有較大改進空間，具有輕量化的前提條件。

4 輪轂輕量化設(shè)計

4.1 拓撲優(yōu)化設(shè)計

拓撲優(yōu)化設(shè)計主要用于現(xiàn)有設(shè)計中的重量縮減部分，其優(yōu)化流程為：建模、劃分網(wǎng)格、定義載荷及約束，設(shè)置重量優(yōu)化目標(biāo)，最后進行求解。本次的減重目標(biāo)設(shè)定為原始結(jié)構(gòu)的50%，在進行優(yōu)化設(shè)計后，可以看出模型顯示出不同的顏色，其中紅色部分是建議刪減部分，從得到的圖7拓撲優(yōu)化圖中可以看出轂邊緣部分可以進行著重進行減重設(shè)計，基于拓撲優(yōu)化給出的減重范圍，將模型的輪轂厚度由初始的15.55mm縮減至12.92mm，在原輪轂其他結(jié)構(gòu)不發(fā)生變化的前提下，在SolidWorks中進行三維模型的修改。

圖8 拓撲密度云圖

4.2 強度校核

優(yōu)化設(shè)計之后輪轂質(zhì)量由有原來的8.7087kg減少至6.8169kg，重量縮減21.72%，但是是否滿足強度要求仍需進行理論分析，將修改后的模型重新導(dǎo)入有限元中，在與未優(yōu)化之前的約束與加載條件完全相同的情況下，再次進行彎曲疲勞試驗，進行求解，得到圖8改進后輪轂結(jié)構(gòu)等效應(yīng)力云與圖9改進后輪轂結(jié)構(gòu)等效位移云圖。

圖9 改進后輪轂結(jié)構(gòu)等效應(yīng)力云圖

圖10 改進后輪轂等效位移云圖

由圖8、圖9中可以分析出，在該輪轂強度方面最大等效應(yīng)力為144MPa對比未改進之前的最大等效應(yīng)力156MPa，降低了7%，最大等效位移0.13994mm對比未改進之前的最大等效位移0.18015mm，變形程度有了顯著地減輕。最大應(yīng)力與最大變形所在位置仍然處于螺栓與加載軸連接處附近，符合實際中應(yīng)用工況。該次減重優(yōu)化設(shè)計無論從從靜態(tài)結(jié)構(gòu)分析還是輕量化方面，都達到了預(yù)期的設(shè)計目標(biāo)。

5 結(jié)論

本文針對對規(guī)格為14×5.5J的鋁合金輪轂進行輕量化設(shè)計，整體實現(xiàn)減重21.72%。

進行數(shù)次優(yōu)化之后，最大應(yīng)力點與最大變形位置仍出現(xiàn)在螺栓與加載軸連接處附近，在進行實際設(shè)計時，需對此處危險點進行材料或者結(jié)構(gòu)的強化，以使輪轂更為安全可靠。