張曉輝, 謝永勝,曾超,程市,李超
1.內(nèi)燃機可靠性國家重點實驗室,山東 濰坊 261061;2.濰柴動力股份有限公司,山東 濰坊 261061
皮帶輪[1-3]是各種機械設(shè)備中的重要傳動部件,應(yīng)用領(lǐng)域廣泛,其中在汽車工業(yè)、農(nóng)業(yè)機械和輕型紡織工業(yè)等領(lǐng)域應(yīng)用最廣。鑄造法、沖壓法、液壓脹形法和焊接成形法等是皮帶輪的傳統(tǒng)制造方法。近年來,綜合了擠壓、環(huán)軋、滾壓等工藝特點的旋壓成型新技術(shù)成為主流,它是一種先進的少、無切削加工工藝,效率高、材料利用率高、制件力學(xué)性能優(yōu)良。旋壓皮帶輪以其精度高、傳動適應(yīng)性好、重量輕等優(yōu)勢,被廣泛應(yīng)用于發(fā)動機前端輪系中[4-8]。農(nóng)業(yè)機械傳動中,皮帶傳動系統(tǒng)工況錯綜復(fù)雜、作業(yè)條件惡劣,每年都有大量的張緊輪及皮帶輪部件需要更換。機器作業(yè)時,張緊輪失效會造成張緊輪安裝軸承、傳送皮帶和相關(guān)零部件的損壞。旋壓皮帶輪的可靠性直接影響農(nóng)機用發(fā)動機的運行狀況。
傳統(tǒng)有限元分析采用靜力學(xué)分析,得到位移分布和應(yīng)力分布云圖,進行強度、剛度分析,但旋壓皮帶輪的失效通常屬于疲勞失效范疇,靜力學(xué)分析很難確定失效危險點[9-13]。本文根據(jù)旋壓皮帶輪實際受力工況,綜合應(yīng)用有限元分析方法,分析故障旋壓皮帶輪疲勞失效的原因,提出改進方案,解決發(fā)動機旋壓皮帶輪故障問題,對提高發(fā)動機的可靠性具有重要意義。
某發(fā)動機旋壓皮帶輪在試驗過程中發(fā)生故障,如圖1所示。觀察旋壓皮帶輪故障件的斷裂位置及斷面失效形式,發(fā)現(xiàn)斷裂位置發(fā)生在旋壓皮帶輪安裝面的螺栓孔周圍,且失效斷面分層明顯,可以觀察到裂紋擴展痕跡。
圖1 旋壓皮帶輪故障
根據(jù)旋壓皮帶輪的工作狀態(tài),判斷其是在發(fā)動機正常工作載荷下,受周期脈動皮帶力的作用發(fā)生受迫運動,旋壓皮帶輪安裝面的螺栓孔周圍及內(nèi)部出現(xiàn)滑移或錯位,當(dāng)損傷積累達(dá)臨界值時,零部件發(fā)生疲勞破壞,從而導(dǎo)致旋壓皮帶輪斷裂或損壞[14]。
在旋壓皮帶輪結(jié)構(gòu)設(shè)計或故障診斷過程中,通常只考慮正常工況的靜態(tài)工作載荷。但在不同應(yīng)用場合下,發(fā)動機曲軸的周期旋轉(zhuǎn)會產(chǎn)生循環(huán)振動激勵,因而疲勞失效是旋壓皮帶輪部件最常見的失效形式。靜態(tài)失效與疲勞失效存在機理差異,靜力學(xué)分析無法準(zhǔn)確診斷發(fā)動機旋壓皮帶輪部件的失效,需進行疲勞失效分析研究,曲軸周期旋轉(zhuǎn)產(chǎn)生的循環(huán)振動激勵使用AVL公司的多體動力學(xué)Excite Power Unit分析軟件計算或者測試獲得時,建模計算工作量大,仿真滯后。本文采用靜力學(xué)分析,結(jié)合最大主應(yīng)力,對比材料許用疲勞強度,進行疲勞失效診斷。
發(fā)動機旋壓皮帶輪總成由風(fēng)扇連接盤、墊塊、旋壓皮帶輪和軸端法蘭通過螺栓連接組成,如圖2所示。旋壓皮帶輪為薄壁零件,為保證仿真精度,需要在旋壓皮帶輪壁厚方向保持2層以上網(wǎng)格單元,圓角處網(wǎng)格要求細(xì)化。采用四面體十節(jié)點單元(C3D10I)對旋壓皮帶輪進行實體網(wǎng)格劃分,為了提高模型的計算精度,選用 2~5 mm網(wǎng)格,保證在厚度方向上至少2層網(wǎng)格單元,同時避免計算資源不足或計算時間過長,旋壓皮帶輪網(wǎng)格模型如圖3所示;螺栓不考慮螺紋細(xì)節(jié),螺栓帽及螺桿均采用六面體一階單元網(wǎng)格(C3D8I)進行劃分,如圖4所示;風(fēng)扇連接盤、墊塊、軸端法蘭等其余部件按照各自幾何特點均采用二階四面體網(wǎng)格(C3D10I),旋壓皮帶輪總成有限元模型如圖5所示。
1—螺栓;2—風(fēng)扇連接盤;3—墊塊;4—軸端法蘭;5—旋壓皮帶輪。 a)圓角細(xì)化 b)皮帶輪
圖4 螺栓六面體網(wǎng)格 圖5 旋壓皮帶輪總成有限元模型
該失效旋壓皮帶輪材料為 SPHE,為鈑金旋壓成型的一體化結(jié)構(gòu),板材厚度為 2 mm,仿真分析所需的零件材料屬性如表1所示。
表1 零件材料屬性
根據(jù)旋壓皮帶輪總成實際裝配情況,在螺栓法蘭面、風(fēng)扇連接盤、墊塊、旋壓皮帶輪和軸端法蘭的裝配面建立接觸,采用有限滑移,建立各個零部件之間、零部件與螺栓之間的相互接觸關(guān)系,根據(jù)實際模型將容差設(shè)置為0.1,接觸屬性按照摩擦因數(shù)設(shè)置為0.12,如圖6所示。
圖6 接觸模型 圖7 Tie綁定模型
建立部件與螺栓螺紋區(qū)域的綁定約束,根據(jù)實際模型將容差設(shè)置為2.0。螺栓螺紋部分與軸端法蘭用Tie綁定模擬螺紋連接,如圖7所示。
軸端法蘭內(nèi)表面建立固定約束,約束六方向自由度,螺栓施加預(yù)緊力,如圖8所示。旋壓皮帶輪槽內(nèi)與皮帶接觸區(qū)域按照包角和皮帶力的方向建立耦合,并施加發(fā)動機最大扭矩工況對應(yīng)的最大皮帶力F1(旋壓皮帶輪與空調(diào)傳遞力)和F2(旋壓皮帶輪與發(fā)動機曲軸傳遞力)。皮帶力F1、F2如圖9、10所示。皮帶力和螺栓預(yù)緊力分別建立在不同的分析步中。
圖8 約束及螺栓預(yù)緊力 圖9 皮帶力F1 圖10 皮帶力F2
旋壓皮帶輪應(yīng)力云圖如圖11所示(圖中單位為MPa)。由圖11可知,在最大皮帶力工況下,旋壓皮帶輪受到的最大主應(yīng)力為158.4 MPa,高于材料SPHE的許用疲勞強度121.5 MPa,靜強度不滿足設(shè)計要求。
從“走出去”到“引進來”,從“賣全球”到“買全球”,40年來,中國上演了擁抱全球化和不斷擴大開放的“三步曲”。
圖11 旋壓皮帶輪應(yīng)力云圖
模型中最大應(yīng)力位于皮帶輪安裝面螺栓孔周圍,與圖1所示實際故障位置相吻合,仿真模型正確有效。
由于旋壓皮帶輪在螺栓孔周圍受到交變的拉、壓載荷作用產(chǎn)生裂紋,進而導(dǎo)致旋壓皮帶輪大面積失效。
分析旋壓皮帶輪的變形情況,發(fā)現(xiàn)墊塊與軸端法蘭直徑不同,相差3.0 mm,導(dǎo)致旋壓皮帶輪兩側(cè)接觸面受壓不均,危險點正好與軸端法蘭同直徑,如圖12所示。
a)整體變形 b)危險點
由于旋壓皮帶輪變形較大、自身剛度較弱,因此將軸端法蘭直徑增加3 mm,與墊塊的直徑相同。
對優(yōu)化后旋壓皮帶輪進行仿真,計算其在最大皮帶力載荷工況下的應(yīng)力。優(yōu)化后旋壓皮帶輪應(yīng)力云圖如圖13所示(圖中單位為MPa)。
由圖13可知,改進后旋壓皮帶輪受到的最大主應(yīng)力為118.6 MPa,低于材料SPHE的許用疲勞強度限值121.5 MPa,滿足設(shè)計要求。
運用疲勞應(yīng)力有限元分析方法,分析旋壓皮帶輪故障產(chǎn)生原因。采用靜力學(xué)分析,得到最大主應(yīng)力,并對比材料疲勞強度限值的方法,進行旋壓皮帶輪疲勞壽命預(yù)測分析,發(fā)現(xiàn)由于墊塊與軸端法蘭直徑,相差3.0 mm,導(dǎo)致旋壓皮帶輪兩側(cè)接觸面受壓不均,而危險點與軸端法蘭同直徑。提出旋壓皮帶輪的優(yōu)化方案,通過仿真與市場應(yīng)用對優(yōu)化方案進行驗證,改善效果明顯。