渦輪增壓器水冷軸承體冷卻性能仿真研究*

龔金科，章 滔，胡遼平,，張愛明，周 崢，余明果

(1.湖南大學(xué)，汽車車身先進(jìn)設(shè)計(jì)制造國家重點(diǎn)實(shí)驗(yàn)室，長沙 410082； 2.湖南天雁機(jī)械有限責(zé)任公司，衡陽 421005)

前言

渦輪增壓器水冷軸承體是增壓器結(jié)構(gòu)中的重要機(jī)件，其各部分的溫度分布不均勻而產(chǎn)生很大的熱應(yīng)力[1-2]。其冷卻性能直接影響著渦輪增壓器的可靠性，在很大程度上決定了渦輪增壓器的壽命[3-4]。

為降低渦輪增壓器全浮動軸承和密封環(huán)附近的溫度，以提高渦輪增壓器的可靠性，采取了許多優(yōu)化措施[5-8]。但為防止熱倒流現(xiàn)象的發(fā)生普遍采用機(jī)油和冷卻水同時冷卻的方式。雖然冷卻水對水冷軸承體的冷卻性能影響很大，卻缺乏對正常工況下渦輪增壓器水冷軸承體冷卻性能的相應(yīng)研究[9]。

本文中采用計(jì)算流體力學(xué)軟件和有限元軟件聯(lián)合進(jìn)行流固耦合仿真的方法進(jìn)行數(shù)值模擬，實(shí)現(xiàn)了數(shù)據(jù)在不同軟件之間的無縫傳遞，提出了一種渦輪增壓器水冷軸承體溫度場的計(jì)算方法。以某型號渦輪增壓器為例，用該方法對兩種不同冷卻水流動方式進(jìn)行了流固耦合數(shù)值模擬，得到了軸承體溫度分布情況，并與試驗(yàn)進(jìn)行了對比驗(yàn)證。

1 渦輪增壓器水冷軸承體傳熱基本原理

渦輪增壓器水冷軸承體傳熱的數(shù)學(xué)模型包括冷卻水和機(jī)油的流動與傳熱模型、軸承體固體導(dǎo)熱模型和冷卻介質(zhì)與軸承體復(fù)雜結(jié)合面的耦合傳熱模型。

1.1 冷卻水和機(jī)油的流動與傳熱

在冷卻水和機(jī)油的流動與傳熱計(jì)算中，將其視為穩(wěn)態(tài)的湍流運(yùn)動，其流動和傳熱過程都遵從質(zhì)量守恒、動量守恒和能量守恒定律。

質(zhì)量守恒方程：

div(ρU)=0

(1)

式中：ρ為流體密度；U為流體速度矢量。

動量守恒方程：

(2)

式中：p為流體壓力；μ為動力黏度；Su、Sv、Sw為動量守恒方程的廣義源項(xiàng)。

能量守恒方程：

(3)

式中：cp為流體的比熱容；T為流體溫度；K為流體的傳熱系數(shù)；ST是黏性耗散項(xiàng)。

反映湍流脈動量對流場影響的湍流動能方程和湍流應(yīng)力方程可通過k-ε方程得到，其形式為

(4)

(5)

式中：k為湍動能；μt為湍動黏度；Gk為由于平均速度梯度引起的湍動能k的產(chǎn)生項(xiàng)；Sk為湍動能源項(xiàng)；σk為湍動能k對應(yīng)Prandtl數(shù)；ε為湍動耗散率；Sε為湍動耗散源項(xiàng)；σε為湍動能耗散率ε對應(yīng)Prandtl數(shù)；C1ε和C2ε為經(jīng)驗(yàn)常數(shù)，σk=1.0,σε=1.3,C1ε=1.44,C2ε=1.92。

1.2 渦輪增壓器水冷軸承體機(jī)體傳熱

渦輪增壓器水冷軸承體機(jī)體內(nèi)傳熱過程的能量守恒方程為

(6)

1.3 流固耦合邊界的共軛傳熱

由于在渦輪增壓器軸承體中采用水冷能很好地解決熱倒流問題，目前普遍采用機(jī)油和冷卻水同時冷卻的冷卻方式[9-10]。但對于渦輪增壓器水冷軸承體冷卻腔內(nèi)流體與軸承體之間的換熱問題，在流體和固體壁面相互影響下難以預(yù)先確定熱邊界條件[11-13]。

流固耦合傳熱是一個典型的弱耦合問題。它只在邊界上存在熱量交換，其邊界上的溫度和換熱系數(shù)都應(yīng)看成是計(jì)算結(jié)果的一部分，而不是已知條件，關(guān)鍵在于解決流體與固體壁面之間的熱量傳遞問題[14-15]。

在渦輪增壓器流固耦合傳熱邊界上有：

qw|solid=qw|fluid

(7)

當(dāng)黏性流體在貼近壁面附近流動且流速很小時，其相對運(yùn)動可忽略不計(jì)[16]。在渦輪增壓器水冷軸承體壁面流體層處由傅里葉熱定律可知：

qw|fluid=-λgradt

(8)

式中：gradt為貼近壁面法線方向上流體溫度梯度。而對流傳熱的牛頓冷卻公式為

qw|fluid=h(tw-tf)

(9)

式中：h為對流傳熱表面?zhèn)鳠嵯禂?shù)；tw和tf分別為交界面和附近冷卻液的溫度。

由式(8)和式(9)可以得到對流傳熱表面換熱系數(shù)與流體溫度場的關(guān)系式為

(10)

為減少邊界條件和更好體現(xiàn)實(shí)際工作情況，采用數(shù)值法進(jìn)行研究。若只運(yùn)用單一綜合性有限元軟件對正常工況下渦輪增壓器水冷軸承體冷卻性能進(jìn)行數(shù)值模擬，很難在減少計(jì)算量的同時保證較高計(jì)算精度[14]。因此對軸承體固體和冷卻水腔流體分別建模，采用有限元軟件和計(jì)算流體力學(xué)軟件分別仿真，再通過編譯文件實(shí)現(xiàn)數(shù)據(jù)在其之間的傳遞，較好地體現(xiàn)了有限元軟件和計(jì)算流體力學(xué)軟件在固體、流體領(lǐng)域仿真的優(yōu)勢。

結(jié)合Fire和Abaqus軟件特點(diǎn)，按圖1所示流程實(shí)現(xiàn)渦輪增壓器水冷軸承體冷卻性能數(shù)值仿真模型的流固耦合計(jì)算。將Fire計(jì)算所得的流體邊界溫度T與對流傳熱系數(shù)h映射到有限元網(wǎng)格上，利用Msc patran耦合面網(wǎng)格和軸承體內(nèi)表面的溫度自由度，并在Abaqus中將其轉(zhuǎn)換成熱交換邊界條件。再將Abaqus計(jì)算所得的固體壁面溫度映射到面網(wǎng)格節(jié)點(diǎn)上作為流體計(jì)算的壁面條件，在Fire中進(jìn)行計(jì)算，按此步驟重復(fù)迭代耦合直至溫度收斂。

2 渦輪增壓器水冷軸承體冷卻性能仿真

2.1 渦輪增壓器水冷軸承體冷卻性能仿真模型

渦輪增壓器水冷軸承體在實(shí)際工作中，同時傳遞熱量給機(jī)油、冷卻水、渦輪軸、渦輪箱和壓氣機(jī)等，傳熱情況十分復(fù)雜。在進(jìn)行數(shù)值模擬時對軸承體外表面倒角和細(xì)小結(jié)構(gòu)進(jìn)行簡化，忽略對于軸承體傳熱影響很小的部件。

在三維軟件中按照軸承體實(shí)際尺寸建立幾何模型并導(dǎo)入Hypermesh中，劃分機(jī)油腔、冷卻水腔面網(wǎng)格和渦輪增壓器軸承體三維網(wǎng)格。為保證內(nèi)部細(xì)小尺寸結(jié)構(gòu)逼近程度，全部采用四面體劃分軸承體網(wǎng)格。同時為更精確模擬壁面附近的流動和傳熱特性，將機(jī)油腔、冷卻水腔表面網(wǎng)格導(dǎo)入到Fire軟件中，利用其前處理模塊fame對進(jìn)、出口和表面網(wǎng)格以及邊界層網(wǎng)格進(jìn)行細(xì)化,采用四面體、六面體混合劃分網(wǎng)格的方法劃分流體網(wǎng)格。

建立的軸承體三維網(wǎng)格模型如圖2所示，冷卻水腔進(jìn)出水口按照順時針順序分別標(biāo)記為1、2、3、4。機(jī)油的CFD網(wǎng)格模型如圖3所示。正常工況下為避免空腔現(xiàn)象的發(fā)生，冷卻時往往只采用下進(jìn)上出的冷卻水流動方式。為分析水冷對于軸承體冷卻性能的影響，對一邊進(jìn)同邊出和一邊進(jìn)另一邊出兩種冷卻方式，分甲、乙兩種方案分別建立水腔模型，如表1所示。

表1 不同冷卻方式的比較

2.2 邊界條件

本文中研究的是正常工況下渦輪增壓器水冷軸承體穩(wěn)定運(yùn)行時的冷卻情況，其傳熱邊界條件通過相應(yīng)試驗(yàn)測取。當(dāng)渦輪增壓器軸轉(zhuǎn)速為110 000r/min時，控制渦輪進(jìn)氣溫度為950℃，進(jìn)口廢氣流量為0.12kg/s，冷卻水進(jìn)口溫度為75℃，機(jī)油進(jìn)口溫度為100℃。渦輪增壓器水冷軸承體內(nèi)冷卻水和機(jī)油的流動認(rèn)為是三維不可壓縮流動，入口采用壓力邊界，冷卻水進(jìn)水壓力2.5MPa，機(jī)油進(jìn)口壓力為0.5MPa，出口采用自由出口邊界，壁面采用無滑移壁面邊界條件。

2.3 計(jì)算結(jié)果分析

2.3.1 迭代耦合收斂性分析

在軸承體機(jī)油腔和冷卻水腔表面各取一點(diǎn)C和D，并在軸承體其余位置任取兩點(diǎn)A、B，迭代9次后，其收斂情況如圖4所示?？梢钥闯?，甲方案從第7次迭代耦合開始，乙方案從第5次迭代耦合開始，溫度值開始收斂穩(wěn)定。

2.3.2 冷卻水腔表面換熱系數(shù)分布

圖5給出了采用甲、乙兩種方案在Fire中計(jì)算時渦輪增壓器水冷軸承體冷卻水腔表面的換熱系數(shù)。兩者換熱系數(shù)最高區(qū)域皆位于水腔進(jìn)出口附近，達(dá)到4 700W/(m2·K)以上。而在水腔靠近渦輪一端，乙方案換熱系數(shù)在1 269.1～2 314.8W/(m2·K)之間，略低于甲方案換熱系數(shù)。但由于甲方案進(jìn)水口與出水口距離過近造成換熱系數(shù)分布不均勻，在冷卻水腔表面的中間部位存在許多低換熱系數(shù)區(qū)域。

2.3.3 軸承體溫度分布

圖6為甲、乙兩種方案水冷軸承體軸向方向豎直平面上的溫度分布。大量熱負(fù)荷集中在軸承體與渦輪箱接觸處，導(dǎo)致溫度高達(dá)600℃，經(jīng)過冷卻水和機(jī)油起的冷卻后溫度逐漸降低，在軸承體與壓氣機(jī)接觸處溫度最終降至110℃。

甲方案在水冷軸承體的全浮動軸承部位軸溫度從靠近渦輪側(cè)的525℃迅速下降到300℃，雖然渦輪軸上的冷卻效果遠(yuǎn)優(yōu)于乙方案，但卻造成局部熱應(yīng)力過高。而正常工況時渦輪軸處在持續(xù)高速運(yùn)轉(zhuǎn)中，在復(fù)合載荷的作用下，甲方案更容易導(dǎo)致軸的變形。

圖7為采用甲、乙兩種方案水冷軸承體水腔進(jìn)出口附近溫度分布情況。圖中溫度最高區(qū)域位于全浮動軸承部位，這是因?yàn)殡m然熱量最主要是從軸承體渦輪端傳向壓氣機(jī)端，但軸在高速運(yùn)轉(zhuǎn)中產(chǎn)生的部分熱量傳遞至軸承體，引起局部溫度過高。

受不同冷卻方式影響，與乙方案相比，甲方案4個進(jìn)出水口位置溫度分布極為不均勻，特別是沒有冷卻水通過的進(jìn)出水口附近溫度達(dá)187.5℃以上。

圖8為甲、乙兩種方案水冷軸承體軸向方向水平面上的溫度分布情況。結(jié)合圖6～圖8可以看出：從渦輪增壓器水冷軸承體的渦輪端至壓氣機(jī)端冷卻水和機(jī)油起到了很好的冷卻作用，顯著減少了熱量的傳遞，特別是在水腔和油腔附近軸承體溫度下降極為明顯；但在冷卻水腔附近，乙方案溫度下降要普遍快于甲方案，且靠近壓氣機(jī)的軸承體部分溫度更低，分布也更均勻，說明乙方案的水腔能更好地將來自渦輪端的熱量隔開，因而冷卻性能更優(yōu)。

3 試驗(yàn)驗(yàn)證

3.1 渦輪增壓器水冷軸承體溫度測點(diǎn)布置

在對渦輪增壓器水冷軸承體的溫度分布進(jìn)行數(shù)值模擬時，傳熱邊界條件的確定直接影響仿真結(jié)果的精確度。正常工況下，尾氣熱量通過渦輪箱傳遞給軸承體，一部分被冷卻介質(zhì)帶走，另一部分傳導(dǎo)至壓氣機(jī)。與此同時，軸的高速旋轉(zhuǎn)所產(chǎn)生的一部分熱量通過機(jī)油傳遞至軸承體。為了獲得準(zhǔn)確的邊界條件，在軸承體渦輪端、軸承體壓氣機(jī)端和軸承體浮動軸承部位布置測溫點(diǎn)，采用WRTK-112工業(yè)鎧裝熱電偶進(jìn)行溫度測試。

在渦輪箱與水冷軸承體連接的法蘭內(nèi)布置溫度測點(diǎn)，將熱電偶球頭與渦輪端法蘭壁面之間的距離控制在0.5～1.0mm內(nèi)。同時在水冷軸承體全浮動軸承臨近的壁面開測量孔，孔徑為1mm，測量點(diǎn)底部距離浮動軸承和密封環(huán)壁面為0.5mm。

為驗(yàn)證仿真結(jié)果的準(zhǔn)確性，在軸承體不同位置取測量點(diǎn)，最終確定的軸承體測溫點(diǎn)分布如圖9所示。其中，測溫點(diǎn)1～5用來獲得準(zhǔn)確的邊界條件，6～11是為了驗(yàn)證仿真的準(zhǔn)確性。

3.2 渦輪增壓器水冷軸承體試驗(yàn)結(jié)果分析

測試在穩(wěn)態(tài)運(yùn)行工況進(jìn)行。采用外吹狀態(tài)，控制增壓器轉(zhuǎn)速110 000r/min，待增壓器穩(wěn)定運(yùn)行20min后測取邊界條件和驗(yàn)證條件。

表2為6個驗(yàn)證測點(diǎn)溫度的仿真結(jié)果與試驗(yàn)結(jié)果的對比。由表可見，最大相對誤差僅為3.1%。上述結(jié)果表明：采用計(jì)算流體力學(xué)軟件和有限元軟件相結(jié)合進(jìn)行流固耦合仿真的方法對渦輪增壓器水冷軸承體冷卻性能數(shù)值模擬是切實(shí)可行的方法。

表2 測點(diǎn)計(jì)算結(jié)果與實(shí)測結(jié)果的比較

4 結(jié)論

(1) 通過對正常工況下渦輪增壓器水冷軸承體冷卻機(jī)理的研究，采用計(jì)算流體力學(xué)軟件和有限元軟件聯(lián)合進(jìn)行流固耦合仿真的方法進(jìn)行數(shù)值模擬，提出了一種渦輪增壓器溫度場的計(jì)算方法。該方法為分析渦輪增壓器水冷軸承體各參數(shù)與熱傳導(dǎo)的關(guān)系和揭示溫度場分布規(guī)律提供了一種有效的數(shù)值模擬方法，為冷卻系統(tǒng)的優(yōu)化提供了依據(jù)。

(2) 對不同冷卻方式下渦輪增壓器水冷軸承體溫度分布的模擬結(jié)果表明，在同樣的進(jìn)口條件下，采用冷卻水一邊進(jìn)另一邊出的冷卻方式，溫度分布更為均勻，低溫區(qū)域分布更廣，且軸承受熱應(yīng)力小，冷卻性能優(yōu)于使用一邊進(jìn)同邊出的冷卻水流通方式。

(3) 通過與試驗(yàn)結(jié)果進(jìn)行比較表明，模擬計(jì)算結(jié)果與試驗(yàn)值較為接近，符合渦輪增壓器水冷軸承體正常工作時溫度分布，證明該仿真方法合理，計(jì)算結(jié)果可信。

[1] 周虹偉.渦輪增壓器主要部件結(jié)構(gòu)特性及改進(jìn)研究[J].中國鐵道科學(xué),2004,25(2):72-77.

[2] 李志江.柴油機(jī)廢氣渦輪增壓器的常見故障與使用保養(yǎng)[J].工程機(jī)械與維修,2008(4):198-199.

[3] 于啟洋,楊忠義,張琳.渦輪增壓器的使用與維護(hù)[J].河北水利水電技術(shù),2004(5):39-40.

[4] 龔金科.汽車排放污染及控制[M].北京:人民交通出版社,2005.

[5] 黃智昊.渦輪增壓器發(fā)展前景分析[J].商用汽車,2004(12):88-89.

[6] 胡友安,李曉東.渦輪增壓器蝸殼熱裂紋的試驗(yàn)研究[J].河海大學(xué)學(xué)報,2008,36(6):846-849.

[7] 張俊紅,李志剛,王鐵寧.車用渦輪增壓技術(shù)的發(fā)展回顧、現(xiàn)狀及展望[J].小型內(nèi)燃機(jī)與摩托車,2007,36(1):66-69.

[8] 楊世友,顧宏中,郭中朝.柴油機(jī)渦輪增壓系統(tǒng)研究現(xiàn)狀與進(jìn)展[J].柴油機(jī),2001(4):1-5.

[9] 朱大鑫.渦輪增壓與渦輪增壓器[M].北京:機(jī)械工業(yè)出版社,1992.

[10] 劉海亮,馬永青,韓志引.汽車渦輪增壓發(fā)動機(jī)工作原理及其優(yōu)勢研究[J].價值工程,2011,30(8):57-58.

[11] 吳新濤,何洪,徐思友,等.渦輪增壓器軸承體傳熱計(jì)算方法研究[J].車用發(fā)動機(jī),2010(1):11-14.

[12] 胡友安,李曉東,陳圖鈞.渦殼溫度場和熱應(yīng)力的有限元分析[J].機(jī)械強(qiáng)度,2007,29(1):130-134.

[13] 徐思友,吳新濤,閆瑞乾,等.增壓器軸承和密封環(huán)溫度試驗(yàn)研究[J].車用發(fā)動機(jī),2010(2):35-37.

[14] 傅松,胡玉平,李新才,等.柴油機(jī)缸蓋水腔流動與沸騰傳熱的流固耦合數(shù)值模擬[J].農(nóng)業(yè)機(jī)械學(xué)報,2010,41(4):26-30.

[15] 孔祥謙.熱應(yīng)力有限單元法分析[M].上海:上海交通大學(xué)出版社,1999.

[16] 賈力,方肇洪.高等傳熱學(xué)[M].北京:高等教育出版社,2008.