某型通機消聲器護罩溫度場分析及其實驗研究

譚禮斌，黃 燦，劉小強，余千英，唐 琳，何 丹，冷小麗

（隆鑫通用動力股份有限公司技術中心基礎研究所CFD研究室，重慶400039）

0 引言

通機（General Purpose Engine），是指車用及特殊用途以外的發動機。通機主要由進排氣系統、冷卻系統、箱體、缸頭、缸頭蓋、箱體蓋、空濾器、消聲器及護罩等組成。其中，排氣消聲器在工作時具有較高的溫度，消聲器護罩用于防護排氣消聲器的同時，可以防止人與高溫消聲器直接接觸而造成皮膚灼傷。通機消聲器護罩一般由成型鋼板沖壓制成，相關法規對護罩表面溫度有嚴格的要求。為了滿足法規和外觀設計要求，護罩通常需多輪修模甚至重新開模才能合格，嚴重影響了開發進度，浪費人力物力。隨著計算機仿真技術的發展，計算機輔助工程（CAE）分析在工業領域得到廣泛的應用。基于計算機虛擬仿真開發（SimCenter）的集成化平臺，對工業各領域產品進行虛擬設計與開發，已逐漸成為行業的趨勢[1-2]。目前，基于計算流體動力學（Computational Fluid Dynamics）的仿真分析在航空航天、汽車、通用機械等行業運用廣泛。因此，本文以某型通機為研究對象，采用通用流體分析軟件STAR-CCM+為虛擬仿真平臺，利用數值模擬方法預測某型通機消聲器護罩表面溫度的分布情況，獲得某型通機消聲器護罩表面溫度分布，與實驗實測溫度值對比分析，驗證數值模擬分析方法的可靠性。本文的研究成果對后續通機消聲器護罩的優化設計、提升產品性能具有重要的指導意義。

1 物理模型

某型通機整機三維模型采用Unigraphics NX 8.0（Siemens PLM Software，Italy）建模獲得（圖 1a）。流體仿真分析軟件 STAR-CCM+（Siemens PLM Software，Italy）具有較高的集成度，是幾何前處理、網格劃分、計算求解、后處理功能一體化的集成軟件，且具有自帶的多面體網格技術和ThinMesher（薄壁層網格）技術。在流固耦合模擬中，采用薄壁層網格可以自動探測固體薄壁件，劃分為高質量、均勻分布的類似于棱柱一樣的網格，同時能大大控制網格數量，保證能量方程計算能夠得到很好的收斂。本文采用STARCCM+11.06為分析平臺，基于整機流場特性分析的網格參數控制策略研究結論對整機模型進行網格參數設置，流體計算域采用多面體網格和邊界層網格（邊界層6層）技術進行網格劃分，固體計算域采用多面體網格和薄壁層網格（薄壁層5層）技術進行網格劃分[3]。網格生成后的整機計算域網格模型圖如圖1（b）所示，網格數量約為 600 萬。圖 1（c）為消聲器護罩表面13個溫度測量點。常用的溫度測量方法有紅外熱成像儀、點溫計、溫度傳感器等[4]，本文選用點溫計進行消聲器護罩表面溫度的測量。

圖1 整機模型示意圖

2 數學模型

流體流動滿足三大守恒定律：質量守恒定律、動量守恒定律和能量守恒定律。在流體流動處于湍流狀態時，整個體系還要遵循湍流運輸方程。以上這些守恒定律的數學描述，統稱為控制方程。文中選用STAR-CCM+中提供的Realizable k-ε湍流模型進行數值計算。

湍流控制方程為三維不可壓縮雷諾時均Navier-Stokes方程：

（1）質量守恒方程（連續性方程）[5，6]

該方程表述為在同一時間內某一封閉流體控制體表面的體積流量為零，即流入的體積流量等于流出的體積流量。不可壓縮流體下，密度ρ為常數的質量守恒方程如式（1）所示。

（2）動量守恒方程（Navier-Stokes方程）[6]

該方程可表述為某一流動控制體，其動量的積累速率等于作用于控制體上的外力的總和。

（3）能量守恒方程[7-8]

溫度場分析需要考慮熱對流、對傳導及輻射傳熱等三種熱交換形式。含有熱交換的流動系統必須滿足能量守恒定律，即對某一流動控制體中流體所做功與熱交換傳遞給流體的熱量等于流體能量的增值。以溫度T為變量的能量守恒方程如式（3）所示。

式中：ui、uj是平均速度分量，m/s；xi、xj為坐標分量，m；p是流體微元體上的壓力，Pa；μeff是湍流有效黏性系數，Pa·s；T 是溫度，K；λ 為流體換熱系數，W/m2.K；Cp為流體比熱容，J/（kg·K）；ST是流體內熱源和由粘性作用引起流體機械能轉變為熱能，J.

（4）k-ε湍流模型方程[9]

k-ε湍流模型方程是基于湍流動能和湍流耗散率的半經驗公式，其中k方程（湍動能方程）為精確方程，而ε方程（湍流耗散率方程）則是由經驗公式經過推導形成的。k-ε湍流模型方程如式（4）所示。

式中：Gk為速度梯度產生的湍動能項；Gb為浮力產生的湍動能項；YM表為脈動擴張項；C1ε，C2ε，C3ε為經驗常數；σk，σε分別為與湍動能k和耗散率相對應的Prandtl數；Sk和Sε為用戶自定義的源項。

（5）輻射傳熱模型

輻射傳熱模型選擇為Surface to Surface Radiation，類型為Gray Thermal Radiation.表面A1和表面A2之間的輻射換熱能量的計算公式為[10]：

式中：φb12為輻射能量，W；X12、X21為角系數；A1、A2為表面積，m2；Eb1、Eb2為輻射能力，W/m2.固體計算域模型中任意兩個表面，若其表面面積和溫度已知，若角系數一旦確定，即可求出輻射換熱能量。因此問題的關鍵為求解并確定角系數。STAR-CCM+中可以通過View Factor Calculator計算角系數，通過固體熱邊界條件（溫度和換熱系數）的輸入，可以計算出受熱輻射部件的表面溫度，從而可以獲得固體部件的輻射傳熱溫度場[10-11]。

3 邊界條件

本文邊界設置均采用實測邊界進行設置。圖2所示為邊界加載區域示意圖（進出口邊界及熱邊界加載區域）。消聲器入口采用Mass Flow inlet邊界，流量為1.016 g/s，溫度為280℃；空濾器采用Mass Flow Inlet邊界，流量為0.975 g/s.虛擬計算域按照通機實驗室實際大小建立，入口設置為滯止入口（Stagnation Inlet），出口設置為壓力出口（Pressure Outlet）；實驗測試環境溫度為28℃；流體屬性選擇為理想氣體（Ideal Gas）。為了進行通機消聲器護罩溫度場的計算，本文對缸套、進排氣道、缸頭燃燒室、缸體、箱體內表面進行了相應的熱邊界設置，其中缸體、箱體內表面溫度設置為80℃，對流傳熱系數設置為100 W/m2·K.缸套劃分為 Cylinder-up、Cylinder-middle、Cylinder-down三個區域，對應的加載溫度分別為200℃、225℃、250℃，對流換熱系數都設置為300 W/m2·K；缸頭進排氣道（Intake Pipe、Exhaust Pipe）、燃燒室面四個區（Sparkplug、Zone1、Zone2、Zone3） 按空載運行狀態下熱流量進行加載，分別為43 760 W/m2，14 850 W/m2，169 000 W/m2，112 500 W/m2，72 560 W/m2，25 430 W/m2.

圖2 邊界加載區域示意圖

本文在通機消聲器護罩溫度場分析中，在流體計算域側考慮了輻射傳熱，因此需要設置固體部件材料的輻射發射率。本文將缸體、箱體、消聲器部件等處理為固體，缸體、箱體材料為ADC12，發射率為0.5；消聲器部件材料為ST14，發射率為0.5；其余非金屬部件（風扇罩、空濾器、油箱等部件）材料為PP，發射率為0.9[12]。

4 模擬結果分析

4.1 模型驗證

圖3為整機結構布置圖及發動機主體（箱體、缸體及火花塞）表面溫度分布云圖。圖中火花塞截面處的溫度仿真值為177℃，實驗測試溫度為168℃，兩者間的差異較小，誤差為5.4%，表明本文的熱邊界加載區域的劃分方法及其邊界的取值是合理的，可以用于后續消聲器護罩表面溫度場的仿真預測。

圖3 整機結構布置圖及發動機主體表面溫度分布云圖

4.2 輻射傳熱計算結果

圖4 、圖5分別為消聲器護罩的溫度分布云圖、消聲器護罩的輻射能量圖。從圖4可以看出，消聲器護罩總體平均溫度約為60℃，消聲器護罩表面背離通機風扇側的溫度比正對通機風扇側的溫度略高，原因是較多的冷卻風能夠流經正對風扇側的護罩表面，而流過背離風扇側的護罩表面的冷卻風較少，未能得到較好的冷卻。同理，圖5可以看出，背離通機風扇側的遠離風扇側的輻射能量大，溫度較高。

圖4 消聲器護罩溫度分布云圖

圖5 消聲器護罩輻射能量圖

4.3 實驗驗證

圖6 為消聲器護罩表面13個測點的溫度實測值與仿真值的對比曲線圖。從圖中可以看出，該消聲器護罩表面溫度仿真結果與實測結果變化趨勢基本一致，總體平均誤差約9%.表明本文構建的消聲器護罩溫度場計算方法是可行的，可以有效地預測消聲器護罩表面溫度分布，為通機消聲器護罩設計與開發的前期提供相應的分析支撐。造成測點溫度仿真值與實測值存在誤差的原因可能是仿真溫度提取點與實驗測點選取間可能存在偏差和模擬假定固定不變的材料物理屬性（實際材料物理屬性是隨溫度變化的）。

圖6 消聲器護罩各測點溫度仿真值與實測值的對比分析曲線

5 結論

（1）采用STAR-CCM+流體分析軟件對某型通機消聲器護罩進行了溫度場分析，與點溫計測取的溫度值進行了對比分析，13個測點的溫度仿真結果與溫度實測結果的整體變化趨勢基本一致，總體誤差較小。通過模擬與實驗的對比驗證分析，表明本文構建的通機消聲器護罩溫度場計算方法具有一定的可行性，可以有效地預測消聲器護罩表面的溫度分布情況。

（2）消聲器護罩表面溫度分布合理，護罩表面總體平均溫度約為60℃，消聲器護罩背離通機風扇側的輻射能量大，溫度較高，護罩正對通機風扇側的輻射能量小，溫度較低。消聲器護罩溫度仿真值與實驗值之間的總體平均誤差約為9%，產生誤差的原因主要是溫度測點的提取差異、幾何模型的簡化及材料物理屬性不定的假設。

（3）采用計算流體力學的數值模擬方法可以有效預測消聲器護罩表面溫度分布，評估溫度值的合理性，可為消聲器護罩的設計與開發提供分析支撐，減小實驗次數，降低開發成本，且開發成功的消聲器護罩能夠滿足相應的溫度法規和外觀設計的要求。