王 毅，楊 靖?，鄧幫林，張云飛，孫 承，劉凱敏

(1.湖南大學汽車車身先進設計制造國家重點實驗室，湖南 長沙 410082；2.湖南大學 先進動力總成技術研究中心，湖南 長沙 410082）

試驗與仿真相結(jié)合的發(fā)動機活塞熱負荷分析*

王 毅1，2，楊 靖1，2?，鄧幫林1，2，張云飛1，2，孫 承1，2，劉凱敏1，2

(1.湖南大學汽車車身先進設計制造國家重點實驗室，湖南 長沙 410082；2.湖南大學 先進動力總成技術研究中心，湖南 長沙 410082）

借助硬度塞溫度測試、數(shù)值分析手段，結(jié)合對改進前活塞的溫度場和熱應力的計算，對發(fā)動機改進后活塞的熱負荷狀況進行評估.通過計算活塞傳熱邊界條件，利用有限元分析軟件，對活塞進行溫度場計算，然后利用試驗實測數(shù)據(jù)對其進行標定，在計算誤差小于5%的情況下，得到改進前與改進后活塞在A工況和B工況下的溫度場分布，最后進行了對應工況的熱應力計算.結(jié)果表明在相同工況下發(fā)動機改進后活塞最高溫度均比改進前活塞低，最大熱應力比改進前活塞小.在改進前活塞滿足熱負荷要求的前提下，改進后活塞滿足熱負荷要求.

活塞；硬度塞法；溫度場；熱應力

活塞是發(fā)動機的主要受熱零件，工作時，處于高溫、高壓、高負荷的惡劣環(huán)境，經(jīng)受周期性交變機械負荷和熱負荷，容易發(fā)生故障［1］.熱負荷是造成活塞失效的主要因素之一，隨著發(fā)動機強化指標的不斷提高，它的影響作用也越來越突出［2-5］.為準確計算活塞在各個工況下的應力大小，全面考察活塞的強度，對活塞進行疲勞分析、敲缸分析以及潤滑磨損分析，確定活塞的溫度場分布是進行上述分析的前提［6］.因此，活塞溫度場計算的精度顯得尤為重要.目前，活塞溫度場測試試驗主要是硬度塞法，硬度塞安裝方便，同時能夠保證一定的精度.謝琰等人研究了發(fā)動機活塞的溫度場［4］，但研究工況比較單一.夏飛等人利用有限元方法研究了活塞溫度場及熱應力［5］，但缺乏試驗的驗證，其研究結(jié)果缺乏可信度.由于發(fā)動機生產(chǎn)工藝發(fā)生變化，需要對發(fā)動機活塞進行再設計以滿足現(xiàn)有的工藝要求.活塞熱負荷是否滿足要求是活塞再設計成功與否的重要評價指標.本文通過發(fā)動機多次試驗表明，改進前活塞滿足熱負荷要求.計算活塞傳熱邊界條件，利用ABAQUS有限元分析軟件，對活塞進行溫度場計算，然后利用試驗實測數(shù)據(jù)對其進行標定，在計算誤差小于5%的狀況下，得到改進前與改進后活塞在A工況和B工況下的溫度場分布，最后進行了對應工況的熱應力計算，評估改進后活塞熱負荷是否滿足工程應用要求.

1 活塞溫度測試

1.1 硬度塞回火曲線的標定

硬度塞法測試溫度的原理是利用某種合金經(jīng)過淬火后產(chǎn)生永久性硬度變化，在不同的溫度下進行回火，其表面硬度也將隨之變化的現(xiàn)象.對硬度塞材料而言，其回火溫度和硬度變化最好呈直線關系或近似線性關系，同時要求材料在同一回火溫度下硬度值穩(wěn)定.因此，本次試驗中選用材料為GCr15滾珠軸承鋼.

硬度塞的淬火處理是在氣體保護電爐中進行，硬度測量采用維氏硬度計.淬火后抽取10只硬度塞進行硬度測量，在每只硬度塞上取3個測量點.測量數(shù)據(jù)表明，淬火后的硬度塞硬度均勻性很好，硬度值偏差在± 5 HV內(nèi)，滿足硬度塞材料要求，如圖1所示.

圖1 硬度塞淬火后檢測值Fig.1 Values of plug hardness after quenching

硬度塞達到所需硬度要求后，進行硬度塞回火試驗，整理回火數(shù)據(jù)，并利用最小二乘法擬合得到硬度 -溫度曲線，如圖2所示，誤差分析見表1.

圖2 硬度 -回火溫度曲線Fig.2 The curve of temperature under different HVs

表1 最小二乘法誤差分析Tab.1 The error analysis by least square method

1.2 活塞溫度測點布置

活塞主要幾何參數(shù)如表2所示.在滿足活塞潤滑性能的前提下，改進前后活塞結(jié)構對比如圖3所示，改進后活塞油環(huán)槽高度減小，第3氣環(huán)槽高度加大，頂部厚度增加.為了準確預測活塞溫度場分布，本次試驗選用20個測點.根據(jù)活塞結(jié)構，布置硬度塞測點，改進前后活塞測點位置相同，如圖4所示.

表2 活塞主要幾何參數(shù)Tab.2 The main parameterof_piston

1.3 發(fā)動機試驗

本次試驗分4組進行，即改進前后活塞分別在最大扭矩工況(工況A）和最大功率工況(工況B）下進行試驗，試驗條件見表3.試驗時最大限度地滿足外界條件的一致性，保證發(fā)動機前后試驗均在相同狀態(tài)下運行.

圖3 活塞改進前后對比Fig.3 The comparison of design between original and modified

圖4 活塞硬度塞測點布置Fig.4 The position of hardness plug on surface of piston

表3 試驗工況及試驗條件Tab.3 The operating condition specification

1.4 活塞測點溫度及數(shù)據(jù)分析

改進前后活塞在兩種不同工況下的測點溫度分別如圖5和圖6所示.從圖中可以看出，在A工況下，改進后活塞大部分溫度測點比改進前相對應測點低，B工況亦如此.同時，A工況下活塞測點溫度較之B工況下測點溫度波動大，這與發(fā)動機運行狀況相關，溫度波動大表明A工況運行工況不穩(wěn)定.

圖5 A工況下改進前后活塞溫度對比Fig.5 The piston temperature comparison between original and the modified under condition A

圖6 B工況下改進前后活塞溫度對比Fig.6 The piston temperature comparison between original and the modified under condition B

2 活塞熱負荷計算

2.1 活塞有限元模型

活塞有限元模型如圖7所示，本次計算為了保證計算精度，采用了二階四面體，節(jié)點及單元信息見表4.

圖7 活塞改進前與改進后有限元模型Fig.7 FEM model of original and modified piston

表4 活塞網(wǎng)格信息Tab.4 Mesh information of piston

2.2 活塞材料屬性

活塞材料為鍛鋁2A80，其密度為2.77 g/cm3，彈性模量為72 GPa，泊松比為0.33，抗拉強度為375 MPa；其他隨溫度變化屬性如表5所示.

表5 活塞材料屬性Tab.5 Material properties of piston

2.3 活塞傳熱邊界條件的確定

2.3.1 活塞頂部傳熱邊界條件

發(fā)動機工作過程中，燃料燃燒釋放的熱能轉(zhuǎn)化為機械能，該過程中活塞頂部直接與高溫燃氣接觸，熱量通過活塞頂面?zhèn)鞯交钊^部，然后通過冷卻油腔和活塞環(huán)將熱量傳給其他冷卻介質(zhì)［7］.根據(jù)周期瞬態(tài)溫度波 動 理論［5，8］，活 塞 頂面的溫度沿 活 塞 頂法線方向迅速衰減，而這個溫度的波動只發(fā)生在活塞頂面1～2 mm的表層，不對活塞的溫度場產(chǎn)生較大的影響，所以在某一特定工況下，一般將活塞溫度場近似為穩(wěn)定的溫度場.

活塞頂面的傳熱邊界條件包括傳熱系數(shù)和燃氣溫度.此次試驗測試了發(fā)動機的氣缸壓力，已知燃氣壓力及氣缸容積后，根據(jù)式(1）計算氣缸瞬時溫度，B工況氣缸壓力曲線如圖8所示.瞬時放熱系數(shù)采用艾歇伯格(Eickelberg）公式［9］計 算 ，即

圖8 發(fā)動機P-Ψ圖Fig.8 P-Ψmap of engine

式中：um為 活 塞 平 均 速 度；Pg為氣體 瞬 時 壓 力(MPa）；Tg為氣體瞬時溫度(K）.

計算氣缸內(nèi)瞬時壓力和溫度后，根據(jù)式(3）和式(4）計算每個循環(huán)的平均對流放熱系數(shù)和平均溫度［］.

氣缸工作容積隨曲軸轉(zhuǎn)角的變化關系為［8］：

式中：Vs為氣缸總?cè)莘e；εc為發(fā)動機壓縮比；λs為發(fā)動機曲柄連桿比.

為了更精確地分析活塞頂面?zhèn)鳠徇吔?，對活塞頂部進行分區(qū).根據(jù)活塞試驗所得式(6），(7）來進行活塞頂面?zhèn)鳠嵯禂?shù)分區(qū)［6］，根據(jù)本文研究發(fā)動機結(jié)構參數(shù)，計算中N=38.8 mm是式(6）和式(7）的分界點.

2.3.2 活塞內(nèi)腔傳熱邊界條件

活塞內(nèi)腔換熱系數(shù)跟活塞內(nèi)腔溫度與曲軸箱內(nèi)部氣流的流動狀況相關.活塞內(nèi)腔表面的換熱系數(shù)為：

式中：T1，T2和Toil分別為活塞頂部溫度、活塞內(nèi)腔底部溫度和曲軸箱氣體溫度；k，σ分別為活塞的導熱系數(shù)、活塞頂厚度；αoil為活塞內(nèi)腔表面與曲軸箱氣體的換熱系數(shù).

2.3.3 活塞側(cè)面?zhèn)鳠徇吔鐥l件

活塞火力岸、環(huán)槽和裙部換熱系數(shù)比較難確定，目前一般采用經(jīng)驗公式來確定.傳熱模型如圖9所示，影響這些區(qū)域的傳熱系數(shù)的因素主要有：氣膜、油膜、活塞環(huán)、缸套的厚度，還有油膜、活塞環(huán)、缸套的換熱系數(shù)和缸套與冷卻水之間的換熱系數(shù)，式(9）～(14）為活塞側(cè)面區(qū)域的換熱系數(shù)求解公式.

圖9 活塞側(cè)面?zhèn)鳠崮Ｐ虵ig.9 Thermal transfer model of side of piston

火力岸換熱系數(shù)為：

第1氣環(huán)槽下緣換熱系數(shù)為：

第2氣環(huán)槽下緣換熱系數(shù)為：

油環(huán)上下緣表面換熱系數(shù)為：

環(huán)槽內(nèi)側(cè)表面換熱系數(shù)為：

活塞裙部表面換熱系數(shù)為：

式中：λ1為燃氣的導熱系數(shù)；λ2為缸套的導熱系數(shù)；λ3為活塞環(huán)的導熱系數(shù)；λ0為機油的導熱系數(shù)；αw為缸套和水腔之間的換熱系數(shù)；Δ2為活塞裙部與缸套的間隙；b為缸套厚度；l0為油環(huán)厚度；d為環(huán)傳熱中心間距；n1，n2分別為油環(huán)與第3氣環(huán)的環(huán)外間隙.活塞兩個工況下各個區(qū)域的對流換熱邊界條件見表6.

表6 2種工況活塞傳熱邊界條件Tab.6 Two kinds of piston boundaries of thermal transfer

3 計算結(jié)果及分析

3.1 溫度場計算結(jié)果及分析

不同工況下改進前與改進后活塞與試驗對比如圖10～圖13所示，對比結(jié)果表明計算值與試驗值吻合較好，誤差均小于5%，有效驗證了此次計算的精度達到工程應用要求.

改進前與改進后活塞在2種不同工況下溫度場計算結(jié)果如圖14～圖17所示.從結(jié)果可以看出，兩種狀態(tài)下，改進前最高溫度為298.8℃，改進后活塞最高溫度為278.5℃，改進后最高溫度比改進前有所降低，改進后活塞頂部區(qū)域溫度比改進前低.改進后油環(huán)槽溫度與活塞裙部區(qū)域溫度改進后比改進前活塞略高.這主要是因為改進后活塞油槽結(jié)構所引起的，油環(huán)槽高度降低，活塞頭部熱量更好的向活塞裙部傳遞，活塞裙部的熱量則通過第三氣環(huán)槽的活塞環(huán)散熱，減少了熱量的累積，同時減小了溫度梯度.

圖10 A狀態(tài)下改進前活塞溫度場計算值與試驗值比較Fig.10 The original piston temperature comparison between simulation and experiment under condition A

圖11 A狀態(tài)下改進后活塞溫度場計算值與試驗值比較Fig.11 The modified piston temperature comparison between simulation and experiment under condition A

圖12 B狀態(tài)下改進前活塞溫度計算值與試驗值比較Fig.12 The original piston temperature comparison between simulation and experiment under condition B

圖13 B狀態(tài)下改進后活塞溫度計算值與試驗值比較Fig.13 The modified piston temperature comparison between simulation and experiment under condition B

圖14 A狀態(tài)下改進前活塞溫度場Fig.14 Temperature field of original piston under condition A

圖15 A狀態(tài)下改進后活塞溫度場Fig.15 Temperature field of modified piston under condition A

圖16 B狀態(tài)下改進前活塞溫度場Fig.16 Temperature field of original piston under condition B

圖17 B狀態(tài)下改進后活塞溫度場Fig.17 Temperature field of modified piston under condition B

活塞內(nèi)腔頂部路徑溫度分布如圖18和圖19所示，改進后活塞第3氣環(huán)槽加寬使得活塞散熱更好，因此，A狀態(tài)下，改進前活塞溫度梯度比改進后大，表明沿此路徑，改進前活塞熱應力比改進后活塞大. B狀態(tài)下所表現(xiàn)出來的此種現(xiàn)象更加明顯.

3.2 活塞熱應力分析

根據(jù)溫度場分布，結(jié)合有限元分析軟件，計算出活塞的熱應力，其結(jié)果對評價活塞熱強度具有一定的指導意義.活塞在不同工況下熱應力計算結(jié)果如圖20～圖23所示，結(jié)果表明各處的熱應力均在50 MPa以下，滿足材料強度要求；其最大熱應力主要集中在活塞內(nèi)腔頂部、活塞油環(huán)槽下邊緣、活塞減重孔邊緣處.改進后的活塞熱應力均有明顯下降，尤其是活塞頂部的內(nèi)表面處.

圖18 A狀態(tài)下活塞熱應力分析Fig.18 Thermal stress analysis under condition A

圖19 B狀態(tài)下活塞熱應力分析Fig.19 Thermal stress analysis under condition B

圖20 A狀態(tài)下改進前活塞熱應力分布Fig.20 Thermal stress distribution of original piston under condition A

圖21 A狀態(tài)下改進后活塞熱應力分布Fig.21 Thermal stress distribution of modified piston under condition A

圖22 B狀態(tài)下改進前活塞熱應力分布Fig.22 Thermal stress distribution of original piston under condition B

圖23 B狀態(tài)下改進后活塞熱應力分布Fig.23 Thermal stress distribution of modified piston under condition B

4 結(jié) 論

1）在相同工況下，改進后活塞最高溫度比改進前最高溫度低，改進后活塞頂部區(qū)域溫度比改進前低，油環(huán)槽溫度與活塞裙部區(qū)域溫度改進后比改進前略高.

2）在狀態(tài)A和狀態(tài)B下，改進后活塞熱應力最大值比改進前活塞熱應力均有減小.

3）通過試驗與仿真相結(jié)合分析活塞溫度場及熱應力分布狀況是解決活塞熱負荷問題的有效途徑，結(jié)構改進后溫度分布和熱應力均得到有效改善.

［1］ 楊進華.微型HCCI自由活塞發(fā)動機工作過程實驗及數(shù)值模擬研究［D］.鎮(zhèn)江：江蘇大學能源與動力工程學院，2012. YANG Jin-hua.Study on experiments and numerical simulation of working process of micro HCCI free-piston engine［D］. Zhenjiang：School of Energy and Power Engineering，Jiangsu University，2012.(In Chinese）

［2］ CERIT M，AYHAN V，PARLAK A.Thermal analysis of a partially ceramic coated piston：effect on cold start HC emission in a spark ignition engine［J］.Applied Thermal Engineering，2010，31(2011）：336-341.

［3］ FORMOSA F.Coupled thermodynamic-dynamic semi-analytical model of free piston stirling engines［J］.Energy Conversion and Management，2010，52(2011）：2098-2109.

［4］ 謝琰，席明智，劉曉麗.柴油機縮口四角 ω燃燒室活塞溫度場試驗研究［J］.內(nèi)燃機，2010，2(1）：38-42. XIE Yan，XI Ming-zhi，LIU Xiao-li.Experimental research to neckingωcombustion chamber-piston temperature field［J］.Internal Combustion Engine，2010，2(1）：38-42.(In Chinese）

［5］ 夏飛，蘇鐵熊，馮耀楠，等.基于質(zhì)量限制的全鋼活塞溫度場的三維數(shù)值模擬［J］.拖拉機與農(nóng)用運輸車，2010，37(1）：42-45. XIA Fei，SU Tie-xiong，F(xiàn)ENG Yao-nan，et al.3D numerical simulation for temperature field of steel piston with quality restrictions［J］.Tractor&Farm Transporter，2010，37(1）：42-45.(In Chinese）

［6］ 張桂昌.基于熱機耦合的柴油機活塞系統(tǒng)敲擊噪聲與潤滑研究及優(yōu)化設計［D］.天津：天津大學機械工程學院，2012. ZHANG Gui-chang.Research on impact noise and lubrication of diesel engine piston assembly considering thermal-mechanical coupling loads and optimization design［D］.Tianjin：School of Mechanical Engineering，Tianjing University，2012.(In Chinese）

［7］ 仲杰.活塞噴油振蕩冷卻的穩(wěn)、瞬態(tài)模擬計算及活塞溫度場分析［D］.濟南：山東大學能源與動力工程學院，2012. ZHONG Jie.Steady and transient numerical simulation of oilinject shaking cooling and temperature distribution in piston［D］.Jinan：School of Energy and Power Engineering，Shangdong University，2012.(In Chinese）

［8］ 周龍保.內(nèi)燃機學［M］.北京：機械工業(yè)出版社，2005：43-44. ZHOU Long-bao.Fundamentals of engine［M］.Beijing：Machinery Industry Press，2005：43-44.(In Chinese）

［9］ 王國強.船舶柴油機活塞和缸套溫度場有限元分析［D］.大連：大連海事大學輪機工程學院，2012. WANG Guo-qiang.FEA for the thermal field of the piston and cylinder liner in marine diesel engine［D］.Dalian：College of Marine，Dalian Maritime University，2012.(In Chinese）

Thermal Load Analysis of Engine Piston Combined Experiment with Simulation

WANG Yi1，2，YANG Jing1，2?，DENG Bang-lin1，2，ZHANG Yun-fei1，2，SUN Cheng1，2，LIU Kai-min1，2

（1.State Key Laboratory of Advanced Design and Manufacturing for Vehicle Body，Hunan Univ，Changsha，Hunan 410082，China；2.Research Center for Advanced Powertrain Technology，Hunan Univ，Changsha，Hunan 410082，China）

In order to investigate the thermal load distribution of modified piston，it is necessary to analyze temperature field and thermal stress combined simulation with experiments.The FE software was used to simulate the temperature field of the piston after the boundary conditions of thermal transfer were calculated，and then，the results were verified according to experiment data，and the error of calculation was less than 5%.Finally，temperature field and thermal stress were derived under different operating conditions including conditions A and B.The results of investigation show that the maximum temperature and maximum value of thermal stress of modified piston decrease compared with the original piston.The modified piston meets the demand of thermal load based on the original piston meeting the demand of thermal load.

pistons；plug of hardness method；temperature field；thermal stress

U464.134

A

1674-2974（2014）08-0023-07

2013- 10- 25

國家高技術研究發(fā)展計劃(863計劃）資助項目(2008AA11A114）；湖南省研究生科研創(chuàng)新基金資助項目(CX20108125）

王 毅(1986-），男，湖北利川人，湖南大學博士研究生

?通訊聯(lián)系人，E-mail：yangjing10@vip.com