劉 陽，雷基林，鄧晰文，牛小強(qiáng)，文 均，溫志高

(1.昆明理工大學(xué) 云南省內(nèi)燃機(jī)重點(diǎn)實(shí)驗(yàn)室，昆明 650500；2. 成都銀河動(dòng)力有限公司，成都 610505)

0 概述

柴油機(jī)不斷向高強(qiáng)化、輕量化、緊湊化、低排放等方向發(fā)展，單位時(shí)間內(nèi)放熱量大幅度增加，缸內(nèi)溫度和壓力急劇升高，缸內(nèi)環(huán)境逐漸惡化[1-4]?；钊鳛椴裼蜋C(jī)缸內(nèi)的核心零部件，所承受的熱沖擊越來越高，嚴(yán)重影響柴油機(jī)活塞工作的可靠性與耐久性。將熱障涂層(thermal barrier coatings， TBC)應(yīng)用到柴油機(jī)活塞上，已成為解決柴油機(jī)活塞熱負(fù)荷較高問題的一種途徑。因此開展基于熱障涂層的柴油機(jī)活塞熱分析研究，對降低活塞熱負(fù)荷，延長活塞疲勞壽命及優(yōu)化活塞運(yùn)動(dòng)學(xué)、動(dòng)力學(xué)性能具有重要意義。

國內(nèi)外學(xué)者針對TBC活塞開展了眾多研究。文獻(xiàn)[5]中采用等離子噴涂法在鋁活塞頂面噴涂0.33 mm厚度的氧化鋯陶瓷層，研究發(fā)現(xiàn)活塞第一環(huán)槽處平均溫度下降12 ℃。文獻(xiàn)[6]中研究表明噴涂熱障陶瓷的鋁合金活塞和鋼制活塞的頭部溫度分別下降48%和35%。文獻(xiàn)[7]中通過道路試驗(yàn)研究裝有 TBC活塞的發(fā)動(dòng)機(jī)性能，表明發(fā)動(dòng)機(jī)性能得到了提升。文獻(xiàn)[8]中通過單缸汽油機(jī)試驗(yàn)研究發(fā)現(xiàn)，TBC活塞最多能減少汽油機(jī)冷起動(dòng)43.2%的HC排放。文獻(xiàn)[9]中通過試驗(yàn)研究表明，陶瓷涂層能夠抑制NOx的形成，減少有害氣體的排放。綜上，TBC活塞的研究主要集中在熱障涂層對活塞基體溫度場、發(fā)動(dòng)機(jī)熱效率、燃油經(jīng)濟(jì)性及排放性能的影響等方面，但對TBC活塞系統(tǒng)的熱負(fù)荷和熱變形研究較少。

為此，本文中以某非道路高壓共軌柴油機(jī)鋁合金活塞作為研究對象，基于臺架活塞溫度場測試結(jié)果建立TBC活塞有限元仿真分析模型，研究熱障涂層對活塞基體溫度場、熱負(fù)荷及熱變形的分布特點(diǎn)和變化規(guī)律，為提升活塞性能指標(biāo)提供了參考依據(jù)。

1 活塞溫度場測試

1.1 試驗(yàn)機(jī)型

研究對象為一臺某非道路高壓共軌柴油機(jī)，其中鋁合金活塞的冷卻方式為內(nèi)冷油腔振蕩，燃燒室形狀為縮口ω形。研究機(jī)型的主要性能參數(shù)見表1。

表1 發(fā)動(dòng)機(jī)主要參數(shù)

1.2 試驗(yàn)方案

試驗(yàn)選取TT-K-30熱電偶作為試驗(yàn)熱電偶和引線式傳輸系統(tǒng)對標(biāo)定功率工況下柴油機(jī)活塞的溫度場進(jìn)行測試。該型號熱電偶可測量的溫度范圍為 0 ℃～1 250 ℃，熱電偶容差為0.4%。在活塞溫度場試驗(yàn)之前，預(yù)先對所使用的熱電偶傳感器進(jìn)行0 ℃～400 ℃ 的標(biāo)定?？紤]試驗(yàn)測試過程較短，引線過程中并未引入輔助機(jī)構(gòu)，導(dǎo)線經(jīng)油底殼和油尺孔導(dǎo)出到發(fā)動(dòng)機(jī)外部HIOKI公司的LR8401-21溫度數(shù)據(jù)記錄儀。每支活塞布置4個(gè)測點(diǎn)，如圖1所示，分別位于活塞頂面(測點(diǎn)a、d)、燃燒室中心(測點(diǎn)b)和燃燒室底部(測點(diǎn)c)，測點(diǎn)均位于距離活塞頂面2 mm位置。

圖1 活塞溫度測點(diǎn)分布示意圖

試驗(yàn)過程中，將鋁合金活塞安裝在第1缸，采用等離子噴涂技術(shù)在其頂面制成0.15 mm金屬黏結(jié)層(NiCrAl)+0.35 mm陶瓷層(MgZrO3)的TBC活塞安裝在第3缸，進(jìn)行標(biāo)定功率工況下的溫度場測試。

2 建立活塞有限元仿真模型

2.1 活塞有限元模型的建立

利用HYPERMESH軟件對模型進(jìn)行有限元網(wǎng)格劃分，采用ABAQUS軟件搭建活塞有限元仿真模型，模型包括活塞、鑲?cè)?、活塞銷、半連桿、黏結(jié)層、陶瓷層，如圖2所示。金屬黏結(jié)層作為活塞基體和陶瓷層之間的過渡層，提高陶瓷層與活塞基體之間的結(jié)合強(qiáng)度并緩解兩者之間由于熱膨脹不匹配而產(chǎn)生的應(yīng)力集中，陶瓷層起到隔熱降低活塞基體溫度的作用[10-13]。

圖2 活塞有限元模型及放大區(qū)域

在熱障涂層活塞有限元分析過程中，對活塞頂面進(jìn)行網(wǎng)格細(xì)化，將活塞頂面上的網(wǎng)格節(jié)點(diǎn)映射到熱障涂層，并采用節(jié)點(diǎn)耦合法將熱障涂層與活塞基體建立成一個(gè)整體，對熱障涂層與活塞基體分別賦予各自的材料屬性。

2.2 材料屬性

活塞、鑲?cè)?、活塞銷、連桿的具體熱物性特征參數(shù)如表2所示，各涂層的具體熱物性特征參數(shù)如表3所示。

表2 活塞組件材料屬性

表3 熱障涂層物性參數(shù)

2.3 邊界條件

準(zhǔn)確的熱邊界條件是研究活塞溫度場、熱負(fù)荷、熱變形的基礎(chǔ)，也是決定模型計(jì)算精度的關(guān)鍵因素。采用第三類邊界條件，由經(jīng)驗(yàn)公式計(jì)算對比測試計(jì)算結(jié)果進(jìn)行反復(fù)修正，從而獲得準(zhǔn)確的熱邊界條件。

在穩(wěn)定工況下，活塞整體熱負(fù)荷只有在活塞頂面附近區(qū)域的熱邊界隨缸內(nèi)燃?xì)鉅顟B(tài)的變化而變化，因此僅對活塞頂面、火力岸和第一環(huán)槽位置加載瞬態(tài)變化的熱邊界，其余區(qū)域均加載固定熱邊界。

基于該款柴油機(jī)試驗(yàn)測試數(shù)據(jù)建立一維熱力學(xué)仿真模型，獲得在標(biāo)定功率工況下缸內(nèi)燃?xì)獾臏囟群蛯α鱾鳠嵯禂?shù)，考慮柴油機(jī)缸內(nèi)氣流流動(dòng)變化和缸內(nèi)結(jié)構(gòu)等因素，計(jì)算得到活塞火力面位置的對流傳熱系數(shù)隨著徑向距離變化的無量綱關(guān)系曲線[14]，其中無量綱HTC為活塞頂面沿徑向在不同位置的傳熱系數(shù)與缸內(nèi)最大對流傳熱系數(shù)的比值，無量綱R為頂面沿徑向在不同位置與活塞半徑的比值，如圖3所示。

圖3 對流換熱系數(shù)與徑向距離的無量綱曲線

從而得到活塞火力面徑向任意位置的傳熱系數(shù)αr為：

(1)

式中，αgm為缸內(nèi)燃?xì)馄骄鶎α鱾鳠嵯禂?shù)，W/(m2·K)；N為從活塞中心軸線到最大放熱系數(shù)位置(即喉口位置)的距離，喉口位置距離活塞中心軸線27.5 mm，因此N=27.5 mm，如圖4所示。將活塞頂面區(qū)域分成12個(gè)子區(qū)域，取每個(gè)區(qū)域的中心線到活塞中心的距離作為該區(qū)域?qū)α鱾鳠嵯禂?shù)計(jì)算中的r。

圖4 活塞頂面尺寸

2.4 模型驗(yàn)證

提取鋁合金活塞和0.15 mm金屬黏結(jié)層 +0.35 mm 陶瓷層的TBC活塞仿真計(jì)算的4個(gè)仿真值與試驗(yàn)值進(jìn)行對比分析，見圖5、圖6。仿真計(jì)算溫度與試驗(yàn)測試溫度的偏差均小于1.5%，說明仿真結(jié)果具有較高的準(zhǔn)確度，可以用來進(jìn)行下一步的方案分析研究。

圖5 鋁合金活塞試驗(yàn)溫度與仿真溫度對比

圖6 TBC活塞試驗(yàn)溫度與仿真溫度對比

2.5 方案設(shè)計(jì)

考慮到金屬黏結(jié)層的厚度低于0.10 mm就難以連接活塞頂面與陶瓷層，此次研究以0.15 mm為固定黏結(jié)層厚度，以0.05 mm為一個(gè)厚度梯度增加陶瓷層的厚度。具體方案設(shè)計(jì)如表4所示。

表4 熱障涂層設(shè)計(jì)方案

3 結(jié)果分析

3.1 熱障涂層對溫度場的影響與分析

仿真計(jì)算得到不同方案的溫度場分布，如圖7所示。可發(fā)現(xiàn)各方案活塞基體的最高溫度均出現(xiàn)在活塞喉口處，A方案、B方案、C方案、D方案、E方案的最高溫度分別為320.9 ℃、315.7 ℃、310.7 ℃、306.0 ℃、301.5 ℃，較O方案最高溫度354.0 ℃下降33.1 ℃～52.5 ℃，最大降幅為14.83%，表明熱障涂層能有效阻隔燃燒室內(nèi)的熱量向活塞頭部傳遞，能夠顯著降低活塞頭部的溫度。活塞第二環(huán)岸以下的區(qū)域遠(yuǎn)離熱障涂層，受熱障涂層作用較小，且該區(qū)域傳熱系數(shù)沒有發(fā)生變化，因此溫度差異比較小。綜上，熱障涂層能夠有效降低活塞頭部的溫度，離熱障涂層越遠(yuǎn)的位置受到熱障涂層影響越小。

圖7 各方案活塞溫度場云圖

各方案活塞頂面溫度沿著活塞徑向規(guī)律性分布，為了更好地展示活塞溫度徑向分布規(guī)律，在活塞基體頂面垂直于銷孔方向取路徑P1～P5，對比各方案活塞頂面徑向濕度，如圖8所示?？砂l(fā)現(xiàn)，活塞燃燒室中心凸起區(qū)域的溫度較高，隨著徑向距離的增加和活塞幾何結(jié)構(gòu)的影響，燃燒室溫度逐漸降低，在燃燒室底面達(dá)到最低溫度；隨后活塞溫度隨著徑向距離的增加而升高，并在活塞喉口位置P3處達(dá)到最高溫度；之后隨著徑向距離的增加活塞溫度逐漸降低，但在活塞幾何結(jié)構(gòu)突變位置P4處，活塞溫度產(chǎn)生激變，小幅度上升后持續(xù)降低直到活塞邊緣。由圖8可以看出，各方案活塞頂面徑向溫度分布規(guī)律基本相同，TBC活塞基體陶瓷層厚度每增加0.05 mm，溫度降低3.7 ℃～5.3 ℃。

圖8 各方案活塞頂面徑向溫度分布規(guī)律

3.2 熱障涂層對熱應(yīng)力的影響與分析

為詳細(xì)分析熱障涂層對活塞基體熱應(yīng)力的影響，選取活塞上13個(gè)特征點(diǎn)：頂面(點(diǎn)1、2、7)，喉口(點(diǎn)3)，底圈(點(diǎn)4)，燃燒室中心(點(diǎn)5)，燃燒室過度區(qū)(點(diǎn)6)，第一環(huán)槽下沿(點(diǎn)8)，內(nèi)冷油腔(點(diǎn)9、10、11、12)和內(nèi)腔(點(diǎn)13)，見圖9。將各方案活塞在特征點(diǎn)處的熱應(yīng)力進(jìn)行對比，如圖10所示。

圖9 活塞熱應(yīng)力關(guān)鍵位置編號

圖10 各方案活塞關(guān)鍵位置熱應(yīng)力對比

由圖10可知，在活塞頂面邊緣(點(diǎn)1、7)，喉口(點(diǎn)3)和燃燒室中心(點(diǎn)5)處TBC活塞金屬黏結(jié)層與活塞頂面緊密結(jié)合，且金屬黏結(jié)層的熱膨脹系數(shù)遠(yuǎn)小于活塞的熱膨脹系數(shù)，限制了活塞基體在其結(jié)構(gòu)邊緣與突變處的膨脹變形，陶瓷層厚度增加 0.05 mm 后，產(chǎn)生高于鋁合金活塞0.89 MPa～3.16 MPa 的熱應(yīng)力，表明涂層厚度不能無限制增加。在活塞其他關(guān)鍵位置，由于熱障涂層的隔熱作用降低了活塞的溫度梯度，使得TBC活塞的熱應(yīng)力較鋁合金活塞有顯著的降低，并隨瓷層厚度的增加而整體下降，其中E方案較鋁合金活塞在特征點(diǎn)11處的58.72 MPa有最大降幅(降幅為42.44%)。

3.3 熱障涂層對熱變形的影響與分析

活塞的熱變形直接影響活塞往復(fù)運(yùn)動(dòng)過程中與缸套的配合間隙,是影響活塞運(yùn)動(dòng)學(xué)與動(dòng)力學(xué)性能的一個(gè)重要參數(shù)，由此繪制活塞頂端徑向熱變形的極坐標(biāo)圖，其中TS、ATS表示活塞主、次推力側(cè)，見圖11。由圖11可知，活塞銷偏置0.5 mm導(dǎo)致各方案活塞在沿銷軸方向的徑向熱變形大于垂直于銷軸方向的徑向熱變形；TBC活塞陶瓷層與活塞基體頂面緊密結(jié)合在一起，限制了頂端的徑向熱變形，其中E方案較鋁合金活塞的徑向熱變形減少量最大，為0.077 mm；TBC活塞陶瓷層厚度每增加0.05 mm，徑向熱變形量降低0.004 mm～0.006 mm。A、B、C、D、E、O方案的最小徑向熱變形均出現(xiàn)在90°方向，分別為0.275 mm、0.270 mm、0.266 mm、0.261 mm、0.257 mm、0.298 mm。O方案最大徑向變形量集中在290°～330°，均為0.338 mm，而A、B、C、D、E方案中TBC活塞頂端最大徑向變形量集中在192°～197°，變形量分別為0.290 mm、0.284 mm、0.278 mm、0.273 mm、0.269 mm。

圖11 各方案活塞頂端徑向熱變形

提取0°、30°、60°、180°方向上的特征線，研究活塞軸向方向的徑向熱變形，活塞底端為高度0處，見圖12。由圖可知，各方案活塞頂部徑向熱變形最大，向下逐漸減小, 底部徑向熱變形最小，其中在活塞結(jié)構(gòu)突變的環(huán)槽位置徑向熱變形大幅降低；TBC活塞陶瓷層厚度每增加0.05 mm，各角度徑向熱變形量由底部到頂部降低0～0.006 mm。0°方向，A、B、C、D、E方案最大徑向熱變形分別為0.282 mm、0.276 mm、0.271 mm、0.266 mm、0.261 mm，較鋁合金活塞(0.330 mm)最多減少20.9%；30°方向，A、B、C、D、E方案最大徑向熱變形分別為0.278 mm、0.273 mm、0.268 mm、0.263 mm、0.258 mm，較鋁合金活塞(0.317 mm)最多減少18.6%；60°方向，A、B、C、D、E方案最大徑向熱變形分別為0.274 mm、0.269 mm、0.265 mm、0.260 mm、0.256 mm，較鋁合金活塞(0.304 mm)最多減少15.8%，其中在面窗區(qū)域，由于活塞結(jié)構(gòu)的變化導(dǎo)致活塞徑向熱變形出現(xiàn)小幅度波動(dòng)；180°方向，A、B、C、D、E方案最大徑向熱變形分別為0.288 mm、0.282 mm、0.277 mm、0.272 mm、0.267 mm，較鋁合金活塞(0.325 mm)最多減少17.8%。

圖12 活塞軸向高度下的不同徑向熱變形

綜上所述，活塞徑向熱變形變化幅度的差異導(dǎo)致了活塞熱變形不均勻，可進(jìn)一步為優(yōu)化活塞運(yùn)動(dòng)學(xué)、動(dòng)力學(xué)性能及改進(jìn)配缸間隙提供改進(jìn)方向。

4 結(jié)論

(1) 熱障涂層的作用使得TBC活塞最高溫度較鋁合金活塞最多降低14.83%，且陶瓷層厚度每增加0.05 mm，溫度降低3.7℃～5.3℃。

(2) 金屬黏結(jié)層限制了活塞基體在其結(jié)構(gòu)邊緣與突變處的膨脹變形，產(chǎn)生比鋁合金活塞高0.89 MPa～3.16 MPa的熱應(yīng)力。在活塞其他特征點(diǎn)，陶瓷層的隔熱作用降低了活塞的溫度梯度，使得TBC活塞的熱應(yīng)力較鋁合金活塞有顯著的降低，最大降幅42.44%。

(3) 陶瓷層與活塞基體頂面緊密結(jié)合在一起，限制了頂端的徑向熱變形，TBC活塞較鋁合金活塞的徑向熱變形有最多0.077 mm的降低，陶瓷層厚度每增加0.05 mm，徑向熱變形量降低0.004 mm～0.006 mm?；钊谘劁N軸方向的徑向熱變形大于垂直于銷軸方向的徑向熱變形，變化幅度的差異導(dǎo)致了活塞熱變形不均勻，可進(jìn)一步為優(yōu)化活塞運(yùn)動(dòng)學(xué)、動(dòng)力學(xué)性能及改進(jìn)配缸間隙提供改進(jìn)方向。