張永鋒，趙徐林，王海艷，王慎平

[浦林成山（青島）工業(yè)研究設(shè)計(jì)有限公司，山東 青島 266042]

輪胎的磨耗是一個(gè)十分復(fù)雜的過(guò)程，輪胎的材料、結(jié)構(gòu)、動(dòng)態(tài)力學(xué)性能以及汽車懸架系統(tǒng)、路面條件、駕駛?cè)藛T的駕駛習(xí)慣都會(huì)造成輪胎出現(xiàn)不均勻磨耗，影響輪胎的使用壽命[1-2]。此外，輪胎的不均勻磨耗還會(huì)影響輪胎的通過(guò)噪聲、駕乘舒適性以及牽引和制動(dòng)性等諸多性能，因此輪胎磨耗始終是國(guó)內(nèi)外研究的一個(gè)熱點(diǎn)問(wèn)題[3]。

輪胎磨耗的研究主要有室外路試、室內(nèi)轉(zhuǎn)鼓測(cè)試以及數(shù)值分析方法[4-5]。輪胎磨耗路試周期較長(zhǎng)，因此數(shù)值分析方法近些年被廣泛采用。J.C.CHO等[6]建立了帶有復(fù)雜花紋的輪胎仿真模型，采用等參數(shù)單元的方法分析了胎面負(fù)荷的分布特點(diǎn)，并將負(fù)荷作為仿真輸入條件，建立了室內(nèi)磨耗仿真分析方法，仿真與實(shí)測(cè)結(jié)果一致性較好。W.L.SANG等[7]對(duì)具有簡(jiǎn)單花紋的半鋼子午線輪胎進(jìn)行初始磨耗測(cè)試，同時(shí)建立了輪胎磨耗仿真模型，并將仿真與實(shí)測(cè)結(jié)果進(jìn)行了對(duì)比。李釗[8]建立了輪胎與路面的統(tǒng)一摩擦模型，采用網(wǎng)格更新的方法對(duì)輪胎的磨耗過(guò)程進(jìn)行分析，得到了較為理想的分析結(jié)果。許順凱等[9]采用幾何更新的方法，分析了輪胎在自由滾動(dòng)過(guò)程中的磨耗現(xiàn)象，并與路試輪胎進(jìn)行對(duì)比，得到了相對(duì)近似的磨耗結(jié)果。

以上研究對(duì)象均為半鋼子午線輪胎，本工作以全鋼載重子午線輪胎為研究對(duì)象，參照國(guó)內(nèi)某測(cè)試機(jī)構(gòu)的室內(nèi)轉(zhuǎn)鼓磨耗試驗(yàn)方法，建立與之對(duì)應(yīng)的室內(nèi)磨耗仿真方法，驗(yàn)證仿真方法的可靠性。

1 輪胎磨耗發(fā)生機(jī)理及室內(nèi)測(cè)試方法

1.1 發(fā)生機(jī)理

輪胎主要有自由滾動(dòng)、驅(qū)動(dòng)、制動(dòng)和轉(zhuǎn)彎等工作狀態(tài)，輪胎與地面接觸過(guò)程中，胎面與路面之間會(huì)產(chǎn)生縱向和橫向微觀滑移，進(jìn)而胎面膠會(huì)產(chǎn)生剪切應(yīng)力。研究表明單位面積的摩擦能與輪胎的磨耗具有很大的相關(guān)性，因此本工作用該指標(biāo)評(píng)價(jià)輪胎磨耗性能[10]。

單位面積的摩擦能在微觀上主要分為縱向和橫向，表達(dá)式如下：

式中，Ixy為輪胎接地點(diǎn)總摩擦能密度，Ix為縱向摩擦能密度，Iy為橫向摩擦能密度，τx為縱向剪切應(yīng)力，μx為縱向滑移距離，τy為橫向剪切應(yīng)力，μy為橫向滑移距離。

1.2 室內(nèi)測(cè)試方法

本仿真方法參照國(guó)內(nèi)某輪胎測(cè)試機(jī)構(gòu)的室內(nèi)轉(zhuǎn)鼓磨耗測(cè)試方法，輪胎有自由滾動(dòng)、驅(qū)動(dòng)、制動(dòng)和受橫向力作用等諸多工況，測(cè)試過(guò)程可以選擇不同材質(zhì)的路面，在恒定的滾動(dòng)速度下進(jìn)行一定時(shí)間的磨耗測(cè)試，試驗(yàn)完成后測(cè)量輪胎花紋溝深度。室內(nèi)磨耗測(cè)試如圖1所示，F(xiàn)x為輪胎旋轉(zhuǎn)過(guò)程中受到的驅(qū)動(dòng)力和制動(dòng)力，F(xiàn)y為輪胎在旋轉(zhuǎn)過(guò)程受到的橫向力，F(xiàn)z為垂直負(fù)荷，ω為轉(zhuǎn)鼓旋轉(zhuǎn)角速度。

圖1 室內(nèi)磨耗測(cè)試示意

2 有限元仿真模型的建立

2.1 材料模型

選用有限元分析軟件Abaqus進(jìn)行輪胎的充氣、加載和滾動(dòng)分析[11]。根據(jù)輪胎材料及各部件的力學(xué)特性，選用Yeoh本構(gòu)模型對(duì)橡膠材料進(jìn)行描述，骨架材料如胎體、帶束層結(jié)構(gòu)選用Rebar單元定義。二維充氣模型中，橡膠材料的單元類型分別采用CGAX3H和CGAX4H，骨架材料的單元類型采用SFMGAX1。三維加載分析中，橡膠材料分別采用雜交單元C3D6H和C3D8H。為了節(jié)約計(jì)算成本，有限元模型只考慮輪胎縱溝，簡(jiǎn)化了花紋橫溝和細(xì)小鋼片。

輪輞和轉(zhuǎn)鼓定義為解析剛體，輪胎與輪輞、輪胎與轉(zhuǎn)鼓表面分別采用接觸對(duì)的接觸方式，鑒于橡膠材料與路面之間摩擦模型的復(fù)雜性，本研究采用庫(kù)倫摩擦，摩擦因數(shù)選擇0.85，輪胎與轉(zhuǎn)鼓模型見圖2。

圖2 室內(nèi)磨耗仿真模型

2.2 仿真分析條件

以12R22.5規(guī)格輪胎為研究對(duì)象，參照實(shí)際的室內(nèi)磨耗測(cè)試方法，結(jié)合建立有限元模型，分別對(duì)輪胎進(jìn)行自由滾動(dòng)、驅(qū)動(dòng)、制動(dòng)以及帶有橫向力的數(shù)值分析，總測(cè)試?yán)锍虨? 440 km，詳細(xì)工況如表1所示。

表1 12R22.5規(guī)格輪胎磨耗測(cè)試工況

3 仿真結(jié)果分析

3.1 胎面滑移量的確定

根據(jù)建立的模型和工況條件，進(jìn)行輪胎室內(nèi)磨耗數(shù)值求解。需要指出，Abaqus求解器輸出的滑移是隨時(shí)間積分的結(jié)果，但輪胎磨耗的機(jī)理是胎面與不光滑接觸面之間發(fā)生瞬時(shí)微觀滑移導(dǎo)致，故需要將胎面節(jié)點(diǎn)滑移量進(jìn)行微分處理。處理方式見圖3，將胎面接地區(qū)域所有節(jié)點(diǎn)分配在沿胎面寬度和周向的二維坐標(biāo)系中，節(jié)點(diǎn)區(qū)域要完全覆蓋接地印痕區(qū)域，分別以某測(cè)試階段i和j（j＞i）時(shí)刻胎面周向均勻分布的4個(gè)節(jié)點(diǎn)為例，i時(shí)刻4個(gè)節(jié)點(diǎn)的滑移總量分別為Si1，Si2，Si3，Si4，j時(shí)刻4個(gè)節(jié)點(diǎn)的滑移總量分別為Sj1，Sj2，Sj3，Sj4，每個(gè)節(jié)點(diǎn)滑移量（Δ）的推導(dǎo)過(guò)程如下：

圖3 胎面滑移量處理

根據(jù)以上推導(dǎo)過(guò)程，計(jì)算出胎面所有節(jié)點(diǎn)的縱向和橫向瞬態(tài)滑移。

3.2 各工況下胎面摩擦能仿真結(jié)果

根據(jù)以上滑移距離計(jì)算結(jié)果，結(jié)合胎面節(jié)點(diǎn)的剪切應(yīng)力數(shù)據(jù)，分別計(jì)算出各工況下胎面摩擦能分布情況。

3.2.1 自由滾動(dòng)狀態(tài)

輪胎自由滾動(dòng)狀態(tài)下胎面摩擦能分布曲線見圖4。由圖4可以看出，輪胎花紋溝兩側(cè)的摩擦能較大。輪胎在徑向力作用下，胎面受力產(chǎn)生壓縮效應(yīng)，進(jìn)而導(dǎo)致縱向花紋塊邊緣產(chǎn)生較大的橫向剪切應(yīng)力，在周期性旋轉(zhuǎn)的作用下，花紋溝邊緣出現(xiàn)較嚴(yán)重的磨耗現(xiàn)象。

圖4 輪胎自由滾動(dòng)狀態(tài)下胎面摩擦能分布

3.2.2 驅(qū)動(dòng)和制動(dòng)狀態(tài)

輪胎驅(qū)動(dòng)和制動(dòng)狀態(tài)下胎面摩擦能分布曲線見圖5。由圖5可見，驅(qū)動(dòng)和制動(dòng)狀態(tài)下胎面摩擦能明顯高于自由滾動(dòng)狀態(tài)，這主要是由于輪胎在驅(qū)動(dòng)和制動(dòng)狀態(tài)下，胎面除受橫向剪切力外，還受到較大的縱向剪切力，這顯著加劇了磨耗現(xiàn)象，由此也可以說(shuō)明減少輪胎驅(qū)動(dòng)和制動(dòng)工作狀態(tài)，可大幅減小輪胎磨耗量。

圖5 輪胎驅(qū)動(dòng)和制動(dòng)狀態(tài)下胎面摩擦能分布

3.2.3 橫向力作用

輪胎在正常使用過(guò)程中除了自由滾動(dòng)、驅(qū)動(dòng)和制動(dòng)狀態(tài)外，還受到較多橫向力的作用。輪胎橫向力作用下胎面摩擦能分布見圖6。由圖6可見，橫向力作用下花紋塊中間部位產(chǎn)生相對(duì)較大的摩擦能，導(dǎo)致花紋塊中間部位的磨耗量相對(duì)較大。

圖6 輪胎橫向力作用下胎面摩擦能分布

3.3 磨耗仿真與測(cè)試結(jié)果對(duì)比

根據(jù)測(cè)試設(shè)定的工況，將每個(gè)工況的摩擦能加權(quán)求和，得到的全工況下胎面摩擦能分布見圖7。輪胎磨耗測(cè)試后花紋溝深度的變化見圖8。實(shí)際測(cè)試的花紋溝1—4的磨耗量分別為0.365，0.247，0.378和0.418 mm，與仿真結(jié)果的變化趨勢(shì)一致。由于測(cè)試?yán)锍梯^短，測(cè)試過(guò)程只能看作輪胎初始磨耗階段，胎面摩擦能可以看作與輪胎初始磨耗量正相關(guān)。

圖 7 胎面總摩擦能分布

圖8 花紋溝磨耗情況

4 結(jié)語(yǔ)

輪胎室內(nèi)磨耗測(cè)試具有測(cè)試周期短和成本低等特點(diǎn)。本工作建立了多工況組合的輪胎室內(nèi)磨耗仿真方法，其計(jì)算結(jié)果與磨耗試驗(yàn)結(jié)果具有較好的一致性。說(shuō)明在輪胎新產(chǎn)品開發(fā)階段，采用有限元方法預(yù)測(cè)輪胎初始磨耗性能以及優(yōu)化胎面輪廓是相對(duì)可靠的。