基于多體動力學(xué)的曲軸應(yīng)力分析

孫軍 倪培相 呂世杰 胡玉平

摘要： 利用有限元法對某發(fā)動機(jī)鍛鋼曲軸進(jìn)行多體動力學(xué)計算，分析曲軸在發(fā)動機(jī)實際工作過程中的應(yīng)力分布，確定曲軸斷裂原因。分析結(jié)果表明：正常情況下曲軸的第八曲柄臂與曲柄銷過渡圓角處應(yīng)力最大，為350.0 MPa；當(dāng)?shù)诙鬏S承螺栓松動時，第二曲柄臂與曲柄銷過渡圓角處應(yīng)力最大，達(dá)977.4 MPa；第二主軸承螺栓脫落時，第二曲柄臂與曲柄銷過渡圓角的最大應(yīng)力達(dá)995.3 MPa。在主軸承螺栓松動或脫落的情況下，曲軸曲柄銷圓角的最大應(yīng)力遠(yuǎn)超曲軸材料的屈服強(qiáng)度，曲軸容易在此處產(chǎn)生裂紋并斷裂。

關(guān)鍵詞： 曲軸；斷裂；有限元；多體動力學(xué)；應(yīng)力

中圖分類號：TK423 文獻(xiàn)標(biāo)志碼：A 文章編號：1673-6397（2024）01-0058-06

引用格式： ?孫軍，倪培相，呂世杰，等.基于多體動力學(xué)的曲軸應(yīng)力分析［J］.內(nèi)燃機(jī)與動力裝置，2024，41（1）：58-63.

SUN Jun， NI Peixiang，L ?Shijie， et al.Stress analysis of crankshaft based on multi-body dynamics［J］.Internal Combustion Engine & Powerplant， 2024，41（1）：58-63.

0 引言

曲軸是發(fā)動機(jī)的重要部件之一，決定著發(fā)動機(jī)的工作壽命和可靠性。曲軸在工作過程中承受復(fù)雜的彎曲、扭轉(zhuǎn)交變載荷作用，在周期性變化的各種載荷相互作用下，容易彎曲、扭轉(zhuǎn)變形，甚至產(chǎn)生裂紋和斷裂[1-2]。曲軸形狀復(fù)雜，易出現(xiàn)應(yīng)力集中，特別在曲柄臂與曲柄銷過渡圓角、曲柄銷油孔部位的應(yīng)力集中尤為嚴(yán)重[3]。統(tǒng)計分析表明，80%左右的曲軸斷裂是由彎曲疲勞導(dǎo)致[3]，因此應(yīng)對曲軸的彎曲疲勞強(qiáng)度進(jìn)行重點研究。某6缸發(fā)動機(jī)的鍛鋼曲軸在使用過程中頻繁發(fā)生斷裂事故，但未發(fā)生軸頸拉瓦和化瓦現(xiàn)象，拆卸發(fā)動機(jī)發(fā)現(xiàn)第二主軸承螺栓也發(fā)生斷裂。該曲軸材料組織性能和加工質(zhì)量良好，發(fā)動機(jī)在使用過程中無超載運行情況。本文中采用有限元分析方法對曲軸及相關(guān)零件進(jìn)行多體動力學(xué)計算，分析曲軸在發(fā)動機(jī)實際工作過程中的應(yīng)力分布，確定曲軸斷裂原因。

1 分析方案及過程

1.1 分析方案

通過3種方案對該型號曲軸進(jìn)行多體動力學(xué)計算：1）缸體無異常；2）第二主軸承螺栓松動；3）第二主軸承螺栓脫落。通過比較各個方案下的應(yīng)力分布，確定曲軸斷裂的原因。

1.2 有限元模型的建立

利用AVL Power Unit模塊建立分析模型，為準(zhǔn)確反映曲軸實際工作狀況，采用簡易軸承壁模型代替機(jī)體模型。曲軸應(yīng)力集中最嚴(yán)重的位置為過渡圓角和油孔處，且該曲軸在過渡圓角處斷裂，因此主要計算分析過渡圓角處的應(yīng)力。

1.2.1 幾何模型

應(yīng)用Pro/E建立曲軸和軸承壁的三維幾何模型，如圖1所示。利用該模型可以準(zhǔn)確獲得多體系統(tǒng)動力學(xué)仿真分析所需要的旋轉(zhuǎn)半徑、主軸承軸向坐標(biāo)位置和旋轉(zhuǎn)質(zhì)量等模型參數(shù)。

1.2.2 有限元模型

為保證計算精確性和減少計算量，對原結(jié)構(gòu)進(jìn)行適當(dāng)簡化，忽略小的倒角和曲軸前、后端螺紋孔，保留主要影響強(qiáng)度的曲柄銷過渡圓角。分別采用C3D10M、C3D10二階四面體網(wǎng)格構(gòu)建曲軸和軸承壁的有限元模型，如圖2所示。曲軸的單元數(shù)和節(jié)點數(shù)分別為177 707、108 217個，軸承壁的單元數(shù)和節(jié)點數(shù)分別為2 555、1 296個。

1.3 多體動力學(xué)計算

曲軸軸系和機(jī)體組有限元模型的自由度數(shù)量龐大，無法直接用于動力學(xué)計算，需通過模態(tài)減縮計算減少自由度數(shù)量。采用Abaqus軟件進(jìn)行模態(tài)減縮計算，用EXCITE PU軟件的轉(zhuǎn)化模塊將模態(tài)減縮轉(zhuǎn)化為可用的質(zhì)量和剛度矩陣等文件，進(jìn)行多體動力學(xué)計算[5-6]。

在模態(tài)縮減過程中改變了材料的彈性模量和剛度矩陣，曲軸軸系和機(jī)體組振動的共振頻率和陣型等均隨之改變。但仿真的目的是模擬螺栓松動對曲軸應(yīng)力的影響，確定曲軸失效的原因，因此，忽略共振頻率和振型。

動力學(xué)模型主要由體單元、連接單元、載荷及參數(shù)設(shè)置等組成[7]。體單元包括曲軸軸系、機(jī)體組和連桿組，連接單元包括連桿大頭軸承、主軸承、止推軸承以及活塞缸套導(dǎo)向單元。曲軸連桿機(jī)構(gòu)的多體動力學(xué)模型如圖3所示。

連桿組件包括連桿桿身、大頭蓋、連桿螺栓、連桿大頭軸瓦、連桿小頭襯套。在EXCITE PU軟件中通過對連桿組基本尺寸和材料的設(shè)置，結(jié)合軟件自帶的模塊，將連桿組件簡化為集中質(zhì)量點模型。集中質(zhì)量點共有3處，分別位于連桿小頭中心處、連桿大頭中心處以及連桿組質(zhì)心處。連桿組基本尺寸設(shè)置如圖4所示。

發(fā)動機(jī)工作時，曲柄連桿機(jī)構(gòu)除了傳遞載荷外，自身慣量也是系統(tǒng)動載的主要影響因素。連桿組件大頭集中質(zhì)量點包含連桿組旋轉(zhuǎn)質(zhì)量，小頭集中質(zhì)量點包含連桿組往復(fù)質(zhì)量。由于活塞組的往復(fù)質(zhì)量對動載有一定的影響，因此將活塞組的往復(fù)質(zhì)量并入連桿組件小頭集中質(zhì)量點。連桿組件各集中質(zhì)量點采用二維梁單元連接，其截面屬性為連桿桿身最小截面處的拉、彎、扭截面模量[8]。

該發(fā)動機(jī)氣缸直徑為114 mm，行程為144.5 mm，點火順序為1—5—3—6—2—4，額定轉(zhuǎn)速為2 200 r/min，其余多體動力學(xué)模型主要參數(shù)設(shè)置如表1所示。

選用發(fā)動機(jī)額定轉(zhuǎn)速工況的示功圖數(shù)據(jù)作為載荷邊界的初始條件，EXCITE PU軟件根據(jù)各缸的點火順序自動生成每個時刻各缸的活塞載荷、各燃燒室的爆壓載荷以及曲軸輸出端的反力轉(zhuǎn)矩等。各缸的活塞載荷分別作用于連桿小頭自由度節(jié)點上，各缸燃燒室的爆壓載荷分別作用于相應(yīng)缸蓋的燃燒室自由度節(jié)點上，曲軸的反力轉(zhuǎn)矩作用在曲軸輸出端自由度節(jié)點上。

額定轉(zhuǎn)速工況的示功圖數(shù)據(jù)如圖5所示。該圖為第一缸的氣體壓力變化曲線，其他缸的氣體壓力曲線間隔相應(yīng)的點火角。

發(fā)動機(jī)工作時，活塞在燃燒上止點和換氣上止點附近敲擊缸套，采用線性彈簧單元模擬GUID類型活塞缸套導(dǎo)向單元的缸套導(dǎo)向和活塞敲擊作用。采用GUID單元中推薦的剛度和阻尼進(jìn)行計算[9]，分別為100 kN/mm、20 N · s/mm。

曲軸-機(jī)體耦合動力學(xué)計算使用的主軸承單元全部設(shè)為NONL點對面耦合的非線性彈簧單元。計算時，考慮軸承結(jié)構(gòu)的彈性變形、軸承間隙以及軸承不對中等因素，將每個軸承單元設(shè)置3層彈簧，每層彈簧的剛度不同，使總體剛度呈拋物線分布。根據(jù)經(jīng)驗，通常取軸承的三分之一圓周為傳遞力時的有效接觸部位[10]。止推軸承單元阻止曲軸工作時有較大的軸向竄動，AXBE單元通過線性或非線性彈簧和阻尼模擬止推軸承的作用，彈簧單元按載荷、間隙及等效彈簧數(shù)計算單個彈簧的剛度[11]。

2 結(jié)果分析

隨著曲軸轉(zhuǎn)角變化，最大主應(yīng)力點在不同圓角上變化。用EXCITE PU進(jìn)行多體動力學(xué)計算時取3個循環(huán)，應(yīng)力恢復(fù)時只取最后一個工作循環(huán)的應(yīng)力，按間隔曲軸轉(zhuǎn)角5°輸出最后一個工作循環(huán)的過渡圓角最大主應(yīng)力分布結(jié)果，將生成的結(jié)果文件導(dǎo)入Abaqus，分析曲軸應(yīng)力分布。

2.1 缸體無異常

缸體無異常時，計算所用的有限元模型如圖6所示。當(dāng)螺栓不松動時，軸承蓋和軸承壁間緊密接觸，可視為一個整體。經(jīng)多體動力學(xué)計算，正常情況下曲軸應(yīng)力分布如圖7所示。由圖7可知：在曲軸上未發(fā)現(xiàn)較大的應(yīng)力集中點，曲軸運行正常；曲軸轉(zhuǎn)角為30°時，第八曲柄臂與曲柄銷過渡圓角處應(yīng)力最大，最大主應(yīng)力為350.0 MPa。

2.2 第二主軸承螺栓松動

將第二主軸承處上、下部分網(wǎng)格分別存入不同的Component，下部賦予很小的彈性模量，模擬螺栓松動時的情況。第二主軸承螺栓松動的有限元模型如圖8所示。

螺栓松動時第二曲柄臂與曲柄銷過渡圓角處應(yīng)力分布如圖9所示。由圖9可知：曲軸轉(zhuǎn)角為20°時，第二曲柄臂與曲柄銷過渡圓角處應(yīng)力最大，最大主應(yīng)力為977.4 MPa。

與第二主軸承毗鄰的第三曲柄臂與曲柄銷過渡圓角處應(yīng)力也較大，最大主應(yīng)力達(dá)到644.9 MPa，如圖10所示。其它位置的曲柄銷或者主軸頸圓角應(yīng)力值較小，最大主應(yīng)力不超過400.0 MPa。

2.3 第二主軸承螺栓脫落

完全去掉第二主軸承下缸蓋，模擬下缸蓋脫落情況，有限元模型如圖11所示。

螺栓脫落時第二曲柄臂與曲柄銷過渡圓角處應(yīng)力分布如圖12所示。由圖12可知：曲軸轉(zhuǎn)角為20°時，第二曲柄臂與曲柄銷過渡圓角處應(yīng)力最大，最大主應(yīng)力為995.3 MPa，其它位置的曲柄銷、主軸頸圓角應(yīng)力比較小，最大主應(yīng)力不超過400.0 MPa。

由上述3種狀況下的多體動力學(xué)計算結(jié)果可知：第二主軸承螺栓松動或脫落的情況下，曲柄臂與曲柄銷過渡圓角的最大主應(yīng)力達(dá)到977.4 MPa，遠(yuǎn)超曲軸材料的屈服強(qiáng)度，是該曲軸在第二曲柄臂處斷裂較多的原因。

3 結(jié)論

1）該型號曲軸在正常工作過程中曲柄銷過渡圓角處最大主應(yīng)力為350.0 MPa左右，應(yīng)力較小，曲軸不發(fā)生破壞。

2）主軸承螺栓松動嚴(yán)重影響曲軸應(yīng)力分布，在螺栓松動或脫落的情況下，曲軸曲柄臂與曲柄銷過渡圓角的最大主應(yīng)力達(dá)到977.4 MPa，遠(yuǎn)大于曲軸材料的屈服強(qiáng)度，曲軸容易從此處產(chǎn)生裂紋并發(fā)生斷裂。

3）運用有限元法進(jìn)行多體動力學(xué)計算，能夠?qū)崿F(xiàn)曲軸故障再現(xiàn)，從而可以確定曲軸斷裂原因。

參考文獻(xiàn)：

[1] ?楊連生.內(nèi)燃機(jī)設(shè)計[M].北京：中國農(nóng)業(yè)機(jī)械出版社，1981.

[2] ?趙陽，呂夢國，姜全會，等. 一種調(diào)質(zhì)型汽車發(fā)動機(jī)曲軸的斷裂失效分析[J].金屬熱處理，2013，38（4）：109-111.

[3] ?HOU X， LI Y， JIANG T. Fracture failure analysis of ductile cast Iron crankshaft in a vehicle engine[J].Journal of Failure Analysis and Prevention，2011，11（1）：10-16.

[4] ?WITEK L， SIKORA M， STACHOWICZ F， et al.Stress and failure analysis of the crankshaft of diesel engine[J].Engineering Failure Analysis， 2017，82（6）：703-712.

[5] ?何聯(lián)格，蘇建強(qiáng)，周藍(lán).內(nèi)燃機(jī)曲軸結(jié)構(gòu)可靠性數(shù)值仿真計算的研究現(xiàn)狀與展望[J].重慶理工大學(xué)學(xué)報（自然科學(xué)），2020，34（12）：65-73.

[6] ?陸健.基于多體動力學(xué)的機(jī)體可靠性分析[D].濟(jì)南：山東大學(xué)，2011.

[7] ?郭磊，郝志勇，劉波，等.曲軸強(qiáng)度多體動力學(xué)與有限元子模型法仿真[J].浙江大學(xué)學(xué)報，2009，43（9）：1638-1643.

[8] ?于洪亮. 船舶柴油機(jī)曲軸軸系多體動力學(xué)仿真研究[D].大連：大連海事大學(xué)，2009.

[9] ?MOURELATOS Z P.A crankshaft system model for structural dynamic analysis of internal combustion engines[J].Computers & Structures，2001，79（20/21）：2009-2027.

[10] ?蘇鐵熊，張儒華，蔡坪，等.利用有限元法研究曲軸彎曲應(yīng)力的變化規(guī)律[J].車用發(fā)動機(jī)，1995（4）：35-40.

[11] ?TIAN Q， SUN Y L， LIU C.ElastoHydro dynamic lubricated cylindrical joints for rigid-flexible multibody dynamics[J].Computers & Structures，2013，114/115：106-120.

Stress analysis of crankshaft based on multi-body dynamics

SUN Jun1， NI Peixiang1*， L ?Shijie1， HU Yuping2*

1.TianRun Industrial Technology Co. ， Ltd. ，Weihai 264400， China;

2.School of Energy and Power Engineering， Shandong University， Jinan 250061， China

Abstract： The multi-body dynamic calculation of forged steel crankshaft of an engine is carried out using finite element method， the stress distribution of the crankshaft in the actual working process of the engine is analyzed， and the fracture reason of the crankshaft is determined. The results showed that under normal circumstances， the stress at the transition fillet between the eighth crank arm and the crank pin of the crankshaft is the largest， which is 350.0 MPa. When the second main bearing bolt loosens， the stress at the transition fillet of the second crank arm and crank pin is the largest， reaching 977.4 MPa. When the second main bearing bolt falls off， the maximum stress of the transition fillet between the second crank arm and crank pin reaches 995.3 MPa. In the case of loosening or falling off of the main bearing bolt， the maximum stress value of the crankshaft crankpin fillet far exceeds the yield strength of the crankshaft material， and the crankshaft is easy to crack and fracture from here.

Keywords： ?crankshaft; fracture; finite element; multi-body dynamics; stress

（責(zé)任編輯：郎偉鋒）