張淑芬，靳　麒，王文勝，魏豪杰，侯中華

(河南科技大學(xué) a.土木工程學(xué)院；b.藝術(shù)與設(shè)計學(xué)院，河南 洛陽 471023)

金屬履帶拖拉機(jī)支重輪接觸應(yīng)力分析

張淑芬a，靳麒b，王文勝a，魏豪杰a，侯中華a

(河南科技大學(xué) a.土木工程學(xué)院；b.藝術(shù)與設(shè)計學(xué)院，河南 洛陽 471023)

摘要：針對金屬履帶拖拉機(jī)支重輪斷裂、卷邊時的接觸應(yīng)力進(jìn)行分析計算。首先，對支重輪在工作過程中的靜止、重心偏離、牽引力傾斜與重心偏離3種狀態(tài)進(jìn)行了受力分析，并使用接觸理論進(jìn)行接觸應(yīng)力的理論計算；然后，利用ANSYS軟件建立支重輪結(jié)構(gòu)的三維有限元模型，對3種工況下的接觸應(yīng)力進(jìn)行有限元數(shù)值計算。計算結(jié)果表明：整機(jī)在牽引力傾斜及重心偏離工作狀態(tài)下支重輪的接觸應(yīng)力為最大，尤其對于C1302機(jī)型；最大接觸應(yīng)力在輪緣處。接觸應(yīng)力的理論計算與數(shù)值模擬分析結(jié)果非常接近，吻合較好，其誤差不超過5%。

關(guān)鍵詞：金屬履帶拖拉機(jī)；支重輪；有限元分析；接觸應(yīng)力

基金項目：河南省科技廳科技攻關(guān)基金項目(092102210250)

作者簡介：張淑芬(1953-),女,河南新鄉(xiāng)人,教授,研究方向為機(jī)械結(jié)構(gòu)強(qiáng)度分析.

收稿日期：2014-07-16

文章編號：1672-6871(2015)02-0018-04

中圖分類號：S219.2;TB121

文獻(xiàn)標(biāo)志碼：志碼：A

0引言

中國一拖集團(tuán)有限公司生產(chǎn)的東方紅金屬履帶拖拉機(jī)具有牽引附著性能好、越野性能強(qiáng)、穩(wěn)定性好等突出優(yōu)點(diǎn)，尤其在潮濕、黏重土壤上具有更好的實(shí)用性能[1]。支重輪是金屬履帶拖拉機(jī)行走系統(tǒng)中的重要部件，擔(dān)負(fù)著整車的承重任務(wù)，它與履帶板直接接觸并有相對滾動，露天作業(yè)，是數(shù)量較大的易損部件。近年來，金屬履帶拖拉機(jī)不斷改型，由東方紅75、802、1202到1302，牽引功率不斷增大，但是行走裝置卻基本上無大改動。支重輪的尺寸和技術(shù)要求一直未變，仍為每臺車8對，每側(cè)4對，每輪外徑350 mm，寬45 mm，厚20 mm，材料為RZG310-570，輪緣表面硬度不低于415HB；履帶板材料為ZGMn13，硬度為156HB~229HB[2]。因此，由支重輪-履帶板組成的系統(tǒng)在使用過程中，支重輪的斷裂與表面磨損較為嚴(yán)重，一般400 h即發(fā)生嚴(yán)重擠壓卷邊或斷裂，需更換新輪。

目前，有不少學(xué)者對支重輪接觸問題進(jìn)行了分析探討。文獻(xiàn)[3-5]對履帶拖拉機(jī)支重輪的塑變磨損進(jìn)行了深入研究。文獻(xiàn)[6]對基于振動能耗的間隔阻尼層合支重輪采取了參數(shù)優(yōu)化。文獻(xiàn)[7-8]對大型履帶行走裝置支重輪和履帶板接觸分析建立了數(shù)學(xué)模型。文獻(xiàn)[9-10]利用有限元法分析了履帶板和履帶行走裝置支重輪的接觸特性。

本文針對支重輪與履帶板在整機(jī)處于靜止、重心偏離、牽引力傾斜與重心偏離3種工作狀態(tài)下的接觸應(yīng)力，利用接觸理論[11]進(jìn)行了理論計算；并使用ANSYS軟件[12]對支重輪進(jìn)行數(shù)值模擬分析。比較理論計算結(jié)果與數(shù)值模擬分析結(jié)果的一致性及誤差，其目的是為支重輪結(jié)構(gòu)的改進(jìn)設(shè)計提供參考，從而提升支重輪與支重臺車機(jī)構(gòu)的工作性能。

1模型簡化

因為在工作狀態(tài)下支重輪與履帶板直接接觸，依據(jù)設(shè)計圖紙和制造加工實(shí)物，支重輪表面沿輪厚方向并非弧面，而履帶板上與支重輪接觸的工作面則為平面，輪與板一旦接觸，履帶板相對于地面處于靜止[13]。故可將其簡化為圓柱面與固定平面的接觸，接觸區(qū)域近似為一狹長的矩形，假設(shè)為線接觸，如圖1所示。設(shè)圓柱輪的半徑為R，由文獻(xiàn)[11]可得出式(1)~式(4)。

(1)

(2)

圖1　圓柱面與平面接觸

式中：b為接觸面的寬度，mm；R為輪半徑，mm；p為單位長度上的荷載，N/mm；Q為接觸應(yīng)力，MPa；μ1、μ2為滾動摩擦副的泊松比；E1、E2為滾動摩擦副材料的彈性模量，GPa。

當(dāng)E1=E2=E，μ1=μ2= 0.3時，

(3)

(4)

2支重輪受力分析

對東方紅75、802、1202、1302這4種機(jī)型的支重輪受力情況在理論上進(jìn)行了分析計算[4-5]。各種

表1　東方紅系列拖拉機(jī)的質(zhì)量及牽引力

機(jī)型的質(zhì)量與牽引力如表1所示。

支重輪在工作過程中常見的受力狀態(tài)簡化為3種：靜止?fàn)顟B(tài)、重心偏離狀態(tài)、牽引力傾斜及重心偏離狀態(tài)。

2.1　靜止?fàn)顟B(tài)下支重輪受力分析

支重輪之間為等距，在靜止?fàn)顟B(tài)下的受力情況如圖2所示。G為拖拉機(jī)整機(jī)質(zhì)量。

圖2　靜止?fàn)顟B(tài)受力簡圖

在靜止?fàn)顟B(tài)不考慮重心偏移時支重輪受力

FN=G/8=0.125G。

(5)

2.2　重心偏離狀態(tài)下支重輪受力分析

在重心偏離狀態(tài)下，機(jī)車重心偏離的距離，經(jīng)計算約為300 mm。這時支重輪的受力不再均勻，而是呈線性遞增[1]，受力如圖3所示。

根據(jù)平面平行力系平衡條件，在鉛垂方向投影，即∑Fy=0，可得方程

FN1+FN2+FN3+FN4=G。

(6)

對C點(diǎn)取矩，即∑MC=0，可得方程

3×540.7×FN1+2×540.7×FN2+540.7×FN3=G×(811-300)。

(7)

圖3　重心偏離狀態(tài)受力簡圖

又因為4個力之間為線性關(guān)系，則可得到如下方程：

2FN3-FN2=FN4；

(8)

2FN2-FN1=FN3，

(9)

故由以上4個方程可以得出：

FN4=0.416G。

(10)

由于對稱的原因，則輪子承受的最大荷載為：

(11)

2.3　牽引力傾斜及重心偏離狀態(tài)下支重輪的受力分析

在牽引力傾斜及重心偏離狀態(tài)下，不僅是重力產(chǎn)生了偏離，而且牽引力F也向下傾斜，傾角一般為15°。因此，在垂直方向產(chǎn)生了分力，導(dǎo)致支重輪受力加大，進(jìn)而使接觸應(yīng)力加大，其受力分析如圖4所示。

根據(jù)平面任意力系平衡條件，在鉛垂方向投影和對D點(diǎn)取矩，即∑Fy=0，∑MD=0，可得如下方程：

圖4　牽引力傾角和重心偏離狀態(tài)受力簡圖

FN1+FN2+FN3+FN4=G+F×sin15°;

(12)

3×540.7FN1+2×540.7FN2+540.7FN3+

F×1 815sin15°=G×(811-300)。

(13)

又因為4個力為線性關(guān)系，則可得到方程

2FN3-FN2=FN4；

(14)

2FN2-FN1=FN3。

(15)

由以上4個方程可以得出：

FN4=0.416G+0.442F。

(16)

因為對稱的緣故，則支重輪承受的最大荷載為：

(17)

根據(jù)式(5)、式(11)、式(17)和表1及式(1)、式(2)可得到4種機(jī)型支重輪在3種工況下共計12種情況時，所受荷載、接觸應(yīng)力以及支重輪與履帶板的接觸寬度b，分別見表2、表3和表4。

表2　支重輪在3種狀態(tài)下所受荷載　kN

表3　 支重輪在3種狀態(tài)下的

表4　 支重輪與履帶板在3種狀態(tài)下的

以上對支重輪在3種狀態(tài)下所受載荷和接觸應(yīng)力進(jìn)行了分析計算，其中，未考慮支重輪輪面的拔模斜度和輪在工作過程中所受的振動沖擊等因素的影響。因拖拉機(jī)在工作中支重輪輪面與履帶板接觸的有效長度是變化的，故所受的接觸應(yīng)力也是變化的。

3支重輪有限元模型建立及求解

在生成支重輪結(jié)構(gòu)有限元網(wǎng)格時，考慮了結(jié)構(gòu)應(yīng)力計算對網(wǎng)格密度較高的要求，采用ANSYS軟件中Solid185實(shí)體單元，取單元長度為8 mm，支重輪結(jié)構(gòu)材料為RZG310-570，即材料為精密鑄鋼，其屈服極限σs= 310 MPa，強(qiáng)度極限σb=570 MPa；材料密度ρ= 7 800 kg/m3，彈性模量E=210 GPa，泊松比μ= 0.3。由此建立支重輪有限元模型。

3.1　支重輪邊界條件和荷載施加

模型建立完成后進(jìn)行邊界條件和荷載的施加。邊界條件施加在輪軸孔內(nèi)部以此來固定住支重輪，然后將荷載施加在支重輪輪面上，用接觸理論中公式(3)計算的接觸寬度b，并假設(shè)輪面上的接觸應(yīng)力是均勻分布。邊界條件和荷載施加如圖5所示。

圖5　支重輪的邊界條件和荷載施加

3.2　各種機(jī)型支重輪的求解結(jié)果

對4種機(jī)型在3種狀態(tài)、共計12種情況下的接觸應(yīng)力進(jìn)行有限元分析，得出支重輪的應(yīng)力云圖?？紤]篇幅，這里只給出東方紅1302機(jī)型支重輪在第3種工況下，即牽引力傾角及重心偏移狀態(tài)下的應(yīng)力云圖，如圖6所示。

從圖6可以看出：最大接觸應(yīng)力集中在支重輪的輪面與輪輻支撐處，因此可以得出支重輪發(fā)生斷裂的部位與高應(yīng)力區(qū)是一致的。

輪面的接觸應(yīng)力分布具有邊緣部分大、中間部分小的特點(diǎn)，由于應(yīng)力分布的不均勻，使較低硬度和強(qiáng)度的支重輪很容易產(chǎn)生局部比較大的變形，出現(xiàn)卷邊現(xiàn)象。

4種機(jī)型支重輪接觸應(yīng)力有限元分析結(jié)果如表5所示。對4種機(jī)型在3種狀態(tài)下、12種情況時的接觸應(yīng)力理論計算與有限元分析結(jié)果進(jìn)行比對，見表6。

圖6　牽引力傾角及重心偏移狀態(tài)下的應(yīng)力云圖(單位：MPa)

機(jī)型靜止重心偏離牽引力傾斜及重心偏離7512422728580213724529612021452633151302147275355

表6　4種機(jī)型在3種狀態(tài)下、12種情況時支重輪接觸應(yīng)力理論計算與有限元分析結(jié)果

4結(jié)論

本文針對東方紅金屬履帶拖拉機(jī)支重輪結(jié)構(gòu)，利用接觸理論和ANSYS軟件完成了支重輪接觸應(yīng)力計算與分析，確定了支重輪結(jié)構(gòu)上最大接觸應(yīng)力在輪緣處。從理論計算和ANSYS分析的結(jié)果來看，3種狀態(tài)下的最大接觸應(yīng)力值變化不大。理論計算值與有限元分析值誤差在5%以內(nèi)。支重輪的接觸應(yīng)力在靜止?fàn)顟B(tài)下是最小的，而在牽引力傾斜及重心偏離狀態(tài)下則為最大。

參考文獻(xiàn)：

[1]機(jī)械電子工業(yè)部洛陽拖拉機(jī)研究所.拖拉機(jī)設(shè)計手冊[M].北京:機(jī)械工業(yè)出版社,1994:225-226.

[2]吉林工業(yè)大學(xué).拖拉機(jī)底盤結(jié)構(gòu)設(shè)計圖冊[M].北京:機(jī)械工業(yè)出版社,1994:354-357.

[3]謝敬佩,王文炎,王愛琴,等.履帶拖拉機(jī)支重輪的塑變磨損研究[J].農(nóng)業(yè)機(jī)械學(xué)報,2003,34(15):150-153.

[4]張殿武,謝敬佩,張益華,等.東方紅履帶拖拉機(jī)支重輪形變機(jī)制[J].拖拉機(jī)和農(nóng)用運(yùn)輸車,2001(5):19-22.

[5]張益華,王愛琴,陳殿云,等.支重輪履帶板接觸應(yīng)力的有限元分析[J].拖拉機(jī)和農(nóng)用運(yùn)輸車,2004(1):34-36.

[6]燕碧娟,孫大剛,宋勇,等.基于振動能耗的間隔阻尼層合支重輪參數(shù)優(yōu)化[J].農(nóng)業(yè)工程學(xué)報,2011,27(1):175-179.

[7]閆振華,王國強(qiáng),姚宗偉,等.大型履帶行走裝置支重輪和履帶板接觸分析[J].農(nóng)業(yè)工程學(xué)報,2012,28(12):51-56.

[8]張翼飛,邢星.大型履帶支重輪接觸應(yīng)力分析[J].中國新技術(shù)產(chǎn)品,2014(2):132.

[9]蘇寶程.基于有限元分析的履帶板和履帶行走裝置支重輪接觸特性分析[J].煤礦機(jī)械,2014,35(8):107-109.

[10]范京紅,張曙彩,姚雪俠,等.東方紅-C1202履帶拖拉機(jī)[J].拖拉機(jī)和農(nóng)用運(yùn)輸車,2010(3):87-88.

[11]白明華,劉洪彬.工程彈性力學(xué)基礎(chǔ)[M].北京:機(jī)械工業(yè)出版社,1996:105-106.

[12]胡仁喜,徐東升,李亞東,等.ANSYS13.0機(jī)械與接觸結(jié)構(gòu)有限元分析[M].北京:機(jī)械工業(yè)出版社,2011:60-129.

[13]曹付義,劉洋,周志立.履帶車輛軟地面穩(wěn)態(tài)轉(zhuǎn)向驅(qū)動力計算模型[J].河南科技大學(xué)學(xué)報:自然科學(xué)版,2014,35(2):29-32.