增壓器渦輪箱熱-結(jié)構(gòu)強(qiáng)度分析及試驗(yàn)驗(yàn)證

龍?chǎng)危?陳少林， 靳鵬， 劉朝峰， 劉艷東， 張濤

(湖南天雁機(jī)械有限責(zé)任公司技術(shù)中心， 湖南 衡陽(yáng) 421005)

增壓器渦輪箱熱-結(jié)構(gòu)強(qiáng)度分析及試驗(yàn)驗(yàn)證

龍?chǎng)危?陳少林， 靳鵬， 劉朝峰， 劉艷東， 張濤

(湖南天雁機(jī)械有限責(zé)任公司技術(shù)中心， 湖南 衡陽(yáng) 421005)

為了能夠直接驗(yàn)證熱-結(jié)構(gòu)耦合仿真計(jì)算的精確性，以某型增壓器渦輪箱為研究對(duì)象, 運(yùn)用CFD分析軟件與FEA分析軟件同步耦合仿真計(jì)算渦輪端溫度場(chǎng)，然后將求解得到的溫度場(chǎng)節(jié)點(diǎn)溫度賦給結(jié)構(gòu)分析模型，再求解渦輪箱模型的熱應(yīng)力。在發(fā)動(dòng)機(jī)臺(tái)架上運(yùn)行增壓器，發(fā)動(dòng)機(jī)運(yùn)行時(shí)采用與增壓器CFD仿真計(jì)算相同的邊界參數(shù)，采用紅外熱像儀測(cè)定考察區(qū)域的溫度場(chǎng)分布，通過(guò)非接觸式應(yīng)變測(cè)量系統(tǒng)測(cè)定考察區(qū)域的應(yīng)變分布，提取考察區(qū)域標(biāo)記考察點(diǎn)的溫度及應(yīng)變結(jié)果，并與仿真計(jì)算的溫度及應(yīng)變結(jié)果進(jìn)行對(duì)比分析，結(jié)果表明，采用熱-結(jié)構(gòu)耦合分析方法求解的渦輪箱溫度和應(yīng)變精度能滿足工程設(shè)計(jì)要求。

渦輪箱； 溫度場(chǎng)； 應(yīng)變； 結(jié)構(gòu)強(qiáng)度； 仿真

增壓器是一種壓縮新鮮空氣進(jìn)入氣缸提升發(fā)動(dòng)機(jī)動(dòng)力的裝置。研究渦輪箱的溫度場(chǎng)分布和應(yīng)力分布情況可以有效預(yù)測(cè)渦輪箱的熱-機(jī)械疲勞壽命，降低渦輪箱開(kāi)裂的風(fēng)險(xiǎn)，提高渦輪箱的耐久可靠性。渦輪箱失效最主要的方式是熱應(yīng)力，其次是蠕變和氧化，因此研究渦輪箱熱應(yīng)力是當(dāng)前渦輪箱設(shè)計(jì)研究的重要工作。近幾年來(lái)國(guó)內(nèi)外科技工作者致力于熱應(yīng)力方面的研究[1-4]，熱-機(jī)械應(yīng)力分析理論及熱-機(jī)械疲勞理論都得到了長(zhǎng)足的進(jìn)步，其中疲勞分析理論Sehitoglu方法已運(yùn)用到工程實(shí)踐中[5-7]?；裟犴f爾公司(Honeywell)、美國(guó)材料與試驗(yàn)協(xié)會(huì)(ASTM)及其他科研人員對(duì)增壓器渦輪箱熱-機(jī)械疲勞方面進(jìn)行了深入的研究[8-16]。

1 渦輪端分析模型的建立

增壓器渦輪端模型包括流體模型與固體模型，其中固體模型分為熱分析模型與結(jié)構(gòu)分析模型，包括渦輪箱、渦輪、放氣閥門、簡(jiǎn)化的排氣歧管及排氣管、預(yù)緊螺絲等，流體模型包括渦輪端流道、渦輪、放氣閥門等組成的封閉腔體。

1.1 渦輪端分析模型的離散

渦輪端流體模型(見(jiàn)圖1)總體上可分為旋轉(zhuǎn)域與靜止域，局部可以分為一個(gè)旋轉(zhuǎn)域和一個(gè)固定域。旋轉(zhuǎn)域命名為Rot_Fluid，旋轉(zhuǎn)域的邊界層初始值為0.01 mm，共設(shè)置10層，增長(zhǎng)率為1.2，旋轉(zhuǎn)域網(wǎng)格數(shù)為2 103 642個(gè)。靜止域流道入口段，命名為Housing_Fluid，邊界層初始值為0.02 mm，共設(shè)置10層，增長(zhǎng)率為1.3，靜止域網(wǎng)格數(shù)為836 587個(gè)。

由圖2可見(jiàn)，渦輪端的模型組件由簡(jiǎn)化的排氣歧管、渦輪箱、渦輪轉(zhuǎn)子、放氣閥門、簡(jiǎn)化的排氣管組成，該模型全部采用二階四面體單元。

圖1 CFD仿真計(jì)算模型

圖2 渦輪端有限元離散模型

圖3中標(biāo)記的1～8號(hào)點(diǎn)為數(shù)據(jù)處理使用的考察點(diǎn)，在增壓器相同位置有相應(yīng)的標(biāo)記測(cè)點(diǎn)。

圖3 渦輪箱離散模型考察點(diǎn)位置

1.2 載荷邊界條件

流體計(jì)算邊界條件包括進(jìn)出口邊界、旋轉(zhuǎn)域轉(zhuǎn)速、壁面邊界條件等。發(fā)動(dòng)機(jī)工況為3 500 r/min(0.5% FS)、全負(fù)荷，增壓器渦前溫度(渦輪箱入口溫度)為924 ℃(溫度超過(guò)600 ℃時(shí)，誤差為±1 ℃)，渦前壓力(渦輪箱入口壓力)為59.8 kPa(0.5% FS)，渦后溫度(渦輪箱出口溫度)為834 ℃，渦后壓力(渦輪箱出口壓力)為16.2 kPa(0.5% FS)，渦輪轉(zhuǎn)速為(180 000±500) r/min。分析模型采用k-omega，SST湍流模型。將初始條件設(shè)為出口邊界，初始化后進(jìn)行求解。

渦輪端的邊界則需要與實(shí)際工況相結(jié)合。增壓器安裝在發(fā)動(dòng)機(jī)上運(yùn)行時(shí)受到冷卻風(fēng)機(jī)的強(qiáng)制冷卻(冷卻區(qū)域見(jiàn)圖4)，為了得到較好的計(jì)算結(jié)果，對(duì)強(qiáng)制冷卻區(qū)域進(jìn)行了劃分，不同冷卻區(qū)域的對(duì)流傳熱系數(shù)與層膜溫度不盡相同。

圖4 標(biāo)定轉(zhuǎn)速工況示功圖

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

2.1 流體計(jì)算結(jié)果及分析

20世紀(jì)80年代以后，隨著改革開(kāi)放以及社會(huì)經(jīng)濟(jì)的發(fā)展。閩南語(yǔ)合唱作品更趨多元化，在保留部分原汁原味的“鄉(xiāng)土”色彩外也開(kāi)始融入一些情感與唱腔，使得合唱歌曲更貼近生活氣息。如:無(wú)伴奏女聲合唱領(lǐng)唱《雨來(lái)檐頭流》(泉州作曲家楊雙智)，同聲三部閩南語(yǔ)合唱《十二月果子歌》(臺(tái)灣作曲家劉晏良創(chuàng)作)等作品。至今，這些作品仍然流傳于海峽兩岸，影響甚廣。

從溫度分布來(lái)看，渦輪端氣體膨脹做功，溫度逐步降低。截取流道YZ平面(渦輪徑向)的切面作為考察對(duì)象，由圖5可見(jiàn)，流道進(jìn)口溫度高，沿著流道到渦輪輪轂方向，溫度逐漸降低；旋轉(zhuǎn)域的溫度分布呈外部高內(nèi)部低的規(guī)律。整個(gè)渦輪流道切面溫度最高為922 ℃，在流道入口處，溫度最低為543 ℃，在輪轂表面附近。截取流道ZX平面的切面(渦輪軸向)作為考察對(duì)象, 由圖6可見(jiàn)，渦輪流道入口側(cè)的溫度高于流道出口側(cè)的溫度；放氣閥門孔附近的溫度高于其他位置處的溫度，這是由于在該工況下，放氣閥門處于開(kāi)啟狀態(tài)。由圖7和圖8可見(jiàn)，壓力分布基本與溫度分布一致，在放氣閥門孔處有局部壓力突變。

圖5 流道YZ平面溫度分布

圖6 流道ZX平面溫度分布

圖7 流道YZ平面壓力分布

圖8 流道ZX平面壓力分布

2.2 固體模型溫度場(chǎng)計(jì)算結(jié)果及分析

由圖9可見(jiàn)，渦輪端組件的最高溫度為880 ℃，在渦輪箱喉口位置處；放氣閥門處溫度也比較高，但是閥門軸與軸套溫度梯度較大，這也是造成軸套變形的重要影響因素之一。從渦輪的溫度分布來(lái)看，流體進(jìn)入部位的溫度高，流體流出部位溫度低，這與渦輪實(shí)際溫度分布規(guī)律一致，說(shuō)明溫度結(jié)果具有一定正確性。

圖9 渦輪端溫度分布

渦輪箱溫度分布見(jiàn)圖10。渦輪箱最高溫度為830 ℃，在放氣閥門孔與渦輪出口中間的隔板上，這是由于放氣閥門開(kāi)啟后，高溫排氣直接從放氣閥門孔流出作用在該隔板上，造成較高的溫度。

圖10 渦輪箱溫度分布

2.3 結(jié)構(gòu)計(jì)算結(jié)果及分析

渦輪箱熱應(yīng)力分布見(jiàn)圖11。最大熱應(yīng)力在外壁面的出口法蘭下方，該處結(jié)構(gòu)有局部突變，內(nèi)壁隔板溫度高，造成該處溫度梯度高，熱應(yīng)力大。由于渦輪箱各個(gè)位置的溫度分布不均勻，單獨(dú)考察應(yīng)力不能很好地預(yù)測(cè)疲勞裂紋的位置，而通過(guò)渦輪箱累積塑性應(yīng)變值更能直觀地發(fā)現(xiàn)渦輪箱的風(fēng)險(xiǎn)位置。由圖12可見(jiàn)，喉口、隔板邊緣倒角處、流道外壁面倒角位置、隔板與渦輪箱壁面交接的外壁面累積塑性應(yīng)變均較大，這些位置均可能是渦輪箱裂紋的風(fēng)險(xiǎn)位置。

圖11 渦輪箱等效應(yīng)力分布

圖12 渦輪箱累積塑形應(yīng)變分布

總應(yīng)變包括彈性應(yīng)變、非線性應(yīng)變、熱應(yīng)變等，非線性應(yīng)變包括塑性應(yīng)變、蠕變應(yīng)變。由圖13至圖16可見(jiàn)：渦輪箱的總應(yīng)變最大值為2.94×10-2；熱應(yīng)變?yōu)?.92×10-2，在總應(yīng)變中占據(jù)較大比例；彈性應(yīng)變相對(duì)總應(yīng)變來(lái)說(shuō)數(shù)值較小，彈性應(yīng)變值為2.1×10-3；在局部出現(xiàn)塑性應(yīng)變和蠕變位置的非線性應(yīng)變較大，非線性應(yīng)變值為1.38×10-2，大多數(shù)位置的非線性應(yīng)變?yōu)?。

圖13 渦輪箱溫度總應(yīng)變

圖14 渦輪箱彈性應(yīng)變

圖15 渦輪箱非線性應(yīng)變

圖16 渦輪箱熱應(yīng)變

3 試驗(yàn)及數(shù)據(jù)對(duì)比分析

進(jìn)行了溫度和應(yīng)變的試驗(yàn)測(cè)試以驗(yàn)證仿真的準(zhǔn)確性。試驗(yàn)中分別采用Vic-3D/2D非接觸應(yīng)變測(cè)量系統(tǒng)和紅外熱像儀測(cè)量應(yīng)變和溫度。Vic-3D/2D非接觸應(yīng)變測(cè)量系統(tǒng)采用最佳化的3D數(shù)字圖像相關(guān)性運(yùn)算法則，為機(jī)械試驗(yàn)提供三維/二維空間內(nèi)全視野的形貌、位移及應(yīng)變數(shù)據(jù)測(cè)量。由于此系統(tǒng)運(yùn)用白光散斑相關(guān)性技術(shù)，因此物體的真實(shí)位移可以被測(cè)量,并且可獲得物體表面的每一點(diǎn)上的拉格朗日應(yīng)變張量。

3.1 溫度場(chǎng)試驗(yàn)及應(yīng)變測(cè)量

渦輪箱實(shí)體考察點(diǎn)的位置與圖3標(biāo)記的考察點(diǎn)位置相同。Vic-3D使用DIC(數(shù)字圖像相關(guān))技術(shù)，因此，在測(cè)量之前需在待測(cè)物表面制作斑點(diǎn)。制作好的實(shí)體模型考察點(diǎn)和散斑見(jiàn)圖17。

首先安裝紅外熱像儀，然后安裝非接觸式應(yīng)變測(cè)量相機(jī)，在捕捉影像之前只需進(jìn)行影像對(duì)焦及光圈調(diào)整、表面斑點(diǎn)處理及3D校正。調(diào)校好非接觸式應(yīng)變測(cè)量相機(jī)后，需進(jìn)行相應(yīng)的標(biāo)定工作。在增壓器未運(yùn)行情況下，用指定的溫度標(biāo)定板加熱至高于室溫，通過(guò)標(biāo)定板與待標(biāo)定物的溫度差來(lái)標(biāo)定紅外熱像儀。對(duì)于應(yīng)變測(cè)量相機(jī)的標(biāo)定，需要將拍得的形變過(guò)程影像儲(chǔ)存起來(lái)，過(guò)后重新載入軟件作形變的分析處理。如果待測(cè)物本身已具有足夠的斑點(diǎn)，則不需再做額外的斑點(diǎn)處理。

圖17 渦輪箱標(biāo)記線和散斑圖

3.2 試驗(yàn)結(jié)果及仿真計(jì)算結(jié)果對(duì)比分析

由表1可見(jiàn)，1～4,7～8測(cè)點(diǎn)溫度測(cè)試值與仿真計(jì)算值誤差均小于5%，在誤差允許范圍內(nèi)。5，6測(cè)點(diǎn)的誤差大于5%，這是由于在5，6測(cè)點(diǎn)的正前方有一臺(tái)冷卻風(fēng)機(jī)對(duì)增壓器進(jìn)行強(qiáng)制冷卻。冷卻風(fēng)機(jī)的冷卻過(guò)程是一個(gè)強(qiáng)制對(duì)流傳熱過(guò)程，本次仿真計(jì)算沒(méi)有將該冷卻域加入計(jì)算，僅采用經(jīng)驗(yàn)的對(duì)流傳熱邊界不能完全模擬該冷卻的真實(shí)過(guò)程，因此造成溫度計(jì)算結(jié)果與測(cè)試結(jié)果誤差較大。

表1 溫度計(jì)算結(jié)果與測(cè)試結(jié)果對(duì)比分析

由表2可見(jiàn)，1，2，4，7測(cè)點(diǎn)的總應(yīng)變誤差均小于5%，測(cè)試值與仿真計(jì)算值在誤差允許范圍內(nèi)。由于5、6測(cè)點(diǎn)仿真計(jì)算的溫度結(jié)果與測(cè)試結(jié)果誤差較大(見(jiàn)表1)，導(dǎo)致仿真計(jì)算的熱應(yīng)變誤差較大，而熱應(yīng)變占總應(yīng)變的較大比例，因此，5，6測(cè)點(diǎn)總應(yīng)變的仿真計(jì)算結(jié)果與測(cè)量結(jié)果的誤差較大。3，8測(cè)點(diǎn)的誤差大于5%，原因是被測(cè)渦輪箱的區(qū)域是一個(gè)球面，應(yīng)變測(cè)量相機(jī)捕捉微小變形時(shí)會(huì)造成一定的誤差，同時(shí)還有發(fā)動(dòng)機(jī)臺(tái)架振動(dòng)、高溫使散斑效果變差等因素，累積誤差均會(huì)增大。總之，造成誤差的可能原因包括：1)溫度測(cè)量值與仿真計(jì)算值有一定的誤差；2)應(yīng)變測(cè)量?jī)H能測(cè)定面應(yīng)變；3)發(fā)動(dòng)機(jī)振動(dòng)較大，測(cè)量誤差較大；4)數(shù)據(jù)處理及散斑等原因造成；5)測(cè)量邊界的誤差。

表2 應(yīng)變計(jì)算結(jié)果與測(cè)試結(jié)果對(duì)比分析

4 結(jié)論

a) 在溫度對(duì)比研究中，溫度場(chǎng)試驗(yàn)結(jié)果與仿真計(jì)算結(jié)果總體分布一致，散斑質(zhì)量差的位置溫度誤差大；

b) 在應(yīng)變對(duì)比研究中，由于熱應(yīng)變占據(jù)總應(yīng)變的比重大，因此，溫度誤差大的位置熱應(yīng)變誤差大；

c) 在試驗(yàn)過(guò)程中采用冷卻風(fēng)機(jī)對(duì)增壓器進(jìn)行強(qiáng)制冷卻，而在仿真計(jì)算中沒(méi)有計(jì)算該冷卻域，該位置溫度的測(cè)量值和計(jì)算值誤差較大。

[1]ZienkiewiczOC,MorganK.Finiteelementsandapproximation[M].Swansea:PineridgePress,1983:324-325.

[2]BatheKJ.FiniteElementproceduresinengineeringanalysis[M].Englewood:Prentice-Hall,1982:143-144.

[3]Shiang-WoeiChyuan.Finiteelementsimulationofatwin-cam16-valvecylinderstructure[J].FiniteElementsinAnalysisandDesign,2000,35(7):199-212.

[4] 廖日東，左正興,鄒文勝.對(duì)高速大功率柴油機(jī)用氣缸蓋機(jī)械負(fù)荷有限元分析邊界條件施加方式的探討[J].兵工學(xué)報(bào),2001,21(1)：5-10.

[5]XumingSu,MichaelZubeck,JohnLasecki,etal.Thermo-mechanicalfatigueanalysisofcastaluminiumenginecomponents[J].FiniteElementsinAnalysisandFatigue,2003,30(4):159-173.

[6]NeuRW,HuseyinSehitoglu.Thermo-mechanicalfati-gue,oxidation,creep:PartⅠ.DamageMechanisms[J].MetallurgicalTransactionsA,1989,20:1755-1767.

[7]NeuRW,HuseyinSehitoglu.Thermo-mechanicalfati-gue,oxidation,creep:PartⅡ.LifePrediction[J].MetallurgicalTransactionsA,1989,20:1769-1781.

[8]MichaelJVerrilli,MichaelGCastelli.Thermo-mechanicalfatiguebehaviorofmaterials:secondvolume[C].[S.l.]:ASTM,1994:27-41.

[9]AVL.Turbochargertestrigbriefintroductionintothesystem[M].Graz:AVLListGmbH,2005:1-4.

[10]JohnBHeywood.Internalcombustionenginefundamentals[M].NewYork:McGraw-Hill,1988.

[11]FaridAHDAD,CarolCAI.FiniteElementanalysisofturbinehousing[J].HoneywellTurboTechnologies,2007:7-16.

[12]MichalVajdak.Thermalstressanalysisofturbinehousing[J].HoneywellTurboTechnologies,2006:06-0854.

[13]MickaelCormerais,PascalChesse,Jean-FrancoisHetet.Turbochargerheattransfermodelingundersteadyandtransientconditions[J].Int.J.ofThermodynamics,2009,12(4):193-202.

[14]LicuDN,FindlayMJ,GartshoreIS,etal.Transientheattransfermeasurementsusingasinglewide-bandliquidcrystaltest[J].JournalofTurbomachinery,2000,122:546-552.

[15]FaridAhdad,RagupathyKannusamy,DamodharanB.AnalyticalapproachofTMFpredictionforcomplexloading[C]//AsmeTurboExpo:PowerforLand,Sea,&Air.[S.l.]:ASTM,2010:581-590.

[16]Sun-puiNg,RogerNg,WinnieYu.Bilinearapproximationofanisotropicstress-strainpropertiesofwovenfabrics[J].ResearchJournalofTextile&Apparel,2005,9(4):50-56.

[編輯： 姜曉博]

Thermo-Mechanical Strength Analysis and Experimental Validation for Turbine Housing of Turbocharger

LONG Xin, CHEN Shaolin, JIN Peng, LIU Chaofeng, LIU Yandong， ZHANG Tao

(Technology Center, Hunan Tyen Machinery Co., Ltd., Hengyang 421005, China)

In order to verify the accuracy of thermal structural coupling simulation analysis, a certain type of turbine housing of turbocharger was taken as the research object. The synchronous coupling simulation and calculation of turbine end temperature field was conducted by using CFD and FEA software, the gained joint temperature of temperature field was input to structural analysis model and the thermal stress of turbine housing model was solved. With the same boundary parameters as simulation, the turbocharger ran on the engine test bench and the temperature field and strain distributions of inspection area were measured by using the thermal infrared imager and non-contact strain measurement system respectively. The temperature and stain results were acquired and compared with simulation results. The results show that the temperature and strain accuracy of turbine housing can meet engineering design requirements.

turbine housing; temperature field; strain; structural strength; simulation

2017-01-14；

2017-04-11

龍?chǎng)?1983—)，男，工程師，主要研究方向?yàn)槠囓囉脺u輪增壓器；hntylx@163.com。

10.3969/j.issn.1001-2222.2017.04.011

TK423.5

B

1001-2222(2017)04-0053-06