文丘里管設(shè)計仿真一體化平臺的設(shè)計

閆樂樂, 李楠, 王俊莉, 姚旺, 孟慶當

1.內(nèi)燃機可靠性國家重點實驗室，山東 濰坊 261061；2. 濰柴動力股份有限公司，山東 濰坊 261061

0 引言

柴油機在經(jīng)濟性、動力性、可靠性等方面具有優(yōu)勢，是汽車、工程機械、農(nóng)機、電站和船舶的主要動力來源[1]。對柴油機NOx排放要求越來越高，廢氣再循環(huán)(exhaust gas recirculation, EGR)已成為控制柴油機NOx排放的關(guān)鍵技術(shù)，實現(xiàn)EGR閉環(huán)精確控制尤為重要[2-5]。文丘里管的測量精度較高，對發(fā)動機、整車等應(yīng)用邊界的適應(yīng)性好，已廣泛應(yīng)用于柴油機的設(shè)計和生產(chǎn)制造，利用傳感器測量文丘里管入口與喉口的壓差,可準確計算流過的廢氣流量。

為滿足壓差傳感器量程要求，文丘里管設(shè)計過程中需要對仿真設(shè)計進行多次迭代。本文中基于優(yōu)化設(shè)計軟件Isight、Creo二次開發(fā)工具 Protoolkit、計算流體動力學(xué)(computational fluid dynamics, CFD)仿真分析軟件STAR-CCM+及參數(shù)化設(shè)計原理開發(fā)文丘里管設(shè)計仿真一體化平臺，計算過程中通過不斷優(yōu)化迭代，得到最優(yōu)的文丘里管設(shè)計，實現(xiàn)文丘里管自動優(yōu)化設(shè)計，滿足國六柴油機開發(fā)要求。

1 總體方案

1.1 結(jié)構(gòu)和設(shè)計流程

文丘里管由圓筒段(也稱入口段)、收縮段、喉部和擴散段等組成。喉部直徑小于圓筒段直徑，喉口流速大于入口流速，流經(jīng)文丘里管的廢氣質(zhì)量流量[6-7]

(1)

式中：C為表征通過裝置的實際流量與理論流量之間關(guān)系的流出系數(shù)；β為喉口直徑與入口直徑之比；ε為膨脹因子；d為喉口直徑，mm；Δp為入口與喉口的靜壓差，kPa；ρ為通過的廢氣密度，kg/m3。

文丘里管結(jié)構(gòu)如圖1所示。傳統(tǒng)的文丘里管設(shè)計流程如圖2所示。

圖1 文丘里管結(jié)構(gòu)圖

圖2 傳統(tǒng)文丘里管設(shè)計過程

傳統(tǒng)文丘里管設(shè)計過程包括理論計算、三維數(shù)字模型曲面修改、網(wǎng)格劃分、CFD流場計算、采樣和優(yōu)化設(shè)計等流程[8]。大量重復(fù)性迭代工作枯燥、費時，不僅延長了設(shè)計和仿真周期， 也耗費設(shè)計人員過多精力，存在出錯的潛在風險。對此，本文中依托優(yōu)化設(shè)計軟件Isight和Creo二次開發(fā)工具Protoolkit，基于參數(shù)化設(shè)計原理[9-12]開發(fā)了文丘里管設(shè)計仿真一體化平臺，實現(xiàn)文丘里管自動最優(yōu)化設(shè)計，滿足國六柴油機開發(fā)要求。

文丘里管設(shè)計仿真一體化總體設(shè)計流程如圖3所示。根據(jù)一維熱力學(xué)仿真得到EGR管路直徑D1，新鮮進氣量、EGR率、溫度T、壓力p等參數(shù)，以壓差傳感器量程的1/2(10 kPa，且滿足精度±0.2 kPa)要求作為目標，根據(jù)式(1)反向計算文丘里管喉口直徑d。在現(xiàn)有定型的文丘里管中選取一款與計算的d最接近的文丘里管進行計算，計算應(yīng)用流量系數(shù)和膨脹系數(shù)均默認為1，判斷計算出的壓差Δp是否為8～15 kPa，若符合要求，Isight啟動CCM+單向仿真計算程序，仿真計算壓差滿足目標要求，則定型該方案進行開發(fā)，輸出仿真報告作為開發(fā)依據(jù)；若不滿足目標要求或數(shù)據(jù)庫中所有產(chǎn)品均無法滿足壓差為8～15 kPa要求，則重新輸入結(jié)構(gòu)參數(shù)，在系統(tǒng)中生成三維模型，Isight啟動CCM+進行迭代仿真計算，迭代仿真計算過程通過不斷調(diào)整d以滿足目標要求，滿足目標要求的d作為此次定型方案，輸出仿真報告作為開發(fā)依據(jù)。

圖3 文丘里管設(shè)計仿真一體化總體設(shè)計流程

1.2 輸入?yún)?shù)設(shè)計及優(yōu)化目標

設(shè)計輸入?yún)?shù)分為性能參數(shù)、結(jié)構(gòu)參數(shù)和其他相關(guān)參數(shù)3類，輸入?yún)?shù)分類如表1所示。

表1 輸入?yún)?shù)分類

實測壓差超過壓差傳感量器量程，會導(dǎo)致測量失準，實際測量壓差信號如圖4所示。

圖4 實際壓差信號

實際壓差信號是不斷波動的，由于測量壓差超過傳感器量程限值(傳感器量程為20 kPa)時曲線a按照量程限值進行計算，因此壓差傳感量器得到的壓差信號比實際壓差偏小。為避免壓差超過傳感器量程限值，在設(shè)計文丘里管時將壓差設(shè)計在量程一半左右(如曲線b所示)。在廢氣流量不變情況下，喉口直徑對壓差信號影響顯著，實際應(yīng)用中可通過調(diào)整喉口直徑將壓差信號控制在壓差傳感量器量程一半左右。

2 文丘里管設(shè)計仿真一體化的實現(xiàn)

2.1 過程集成

Isight具有強大的過程集成能力[13]，能夠?qū)⒍鄬W(xué)科代碼集成并使流程自動化，并提供實時監(jiān)控。本文中設(shè)計的系統(tǒng)在Isight基礎(chǔ)上集成Creo二次開發(fā)代碼、CFD仿真分析軟件STAR-CCM+。

Isight集成優(yōu)化界面如圖5所示，箭頭代表數(shù)據(jù)及過程傳遞路徑；Creo代碼實現(xiàn)對三維模型的數(shù)據(jù)傳遞；STAR-CCM+和腳本實現(xiàn)邊界條件以及有限分析求解，輸出壓差信號；Script通過壓差信號判斷喉口直徑變化，當Δp<9.8 kPa時，喉口直徑遞減0.1 mm，當Δp>10.2 kPa時，喉口直徑遞增0.1 mm，觸發(fā)Creo二次開發(fā)代碼及STAR-CCM+腳本計算，直至壓差信號為9.8～10.2 kPa時終止計算。

圖5 Isight集成優(yōu)化界面

2.2 Creo二次開發(fā)

由于錄制Creo宏文件容易造成Creo閃退，終止整個系統(tǒng)，且錯誤不易發(fā)現(xiàn)，本文中采用Protoolkit[14-15]開發(fā)Creo代碼取代錄制Creo宏文件，Protoolkit開發(fā)Creo代碼更加穩(wěn)健，其API函數(shù)涵蓋了Creo的大部分功能，Protoolkit能夠?qū)崿F(xiàn)深層次開發(fā)。

本系統(tǒng)采用Protoolkit異步開發(fā)模式作為獨立程序進行Creo開發(fā)，每次Isight調(diào)用STAR-CCM+計算完成后，結(jié)果不滿足要求時執(zhí)行以下程序：獲取Creo路徑并通過ProEngineerStart()函數(shù)觸發(fā)Creo啟動，獲取模板路徑并通過ProMdlRetrieve()函數(shù)加載模板模型，通過ProModelitemInit()、ProParameterInit()函數(shù)初始化當前模板模型及模型參數(shù)，通過ProParameterValueSet()函數(shù)將參數(shù)賦值給模板對象，通過ProSolidRegenerate()實現(xiàn)模板對象再生；ProIntf3DFileWrite()函數(shù)完成模板對象轉(zhuǎn)碼為.stp格式給STAR-CCM+提供輸入模型。

2.3 STAR-CCM+宏文件

STAR-CCM+宏文件主要集成模型導(dǎo)入、網(wǎng)格劃分、邊界條件輸入、輸出結(jié)果設(shè)置等功能。相比STAR-CCM+中Geometry 3D-CAD models模塊，Creo是更為專業(yè)的繪圖軟件?？紤]測壓部位形狀會隨著喉口位置發(fā)生微調(diào)整，選擇繪圖Creo進行三維建模，保證三維模型的魯棒性，進而保證迭代計算正常進行。

3 仿真實例

3.1 發(fā)動機參數(shù)

某6 L發(fā)動機設(shè)計參數(shù)如表2所示。

表2 發(fā)動機設(shè)計參數(shù)

3.2 計算結(jié)果

經(jīng)過7輪迭代優(yōu)化計算，壓差由8.49 kPa優(yōu)化至10.03 kPa，滿足9.8～10.2 kPa要求，喉口直徑由15.5 mm優(yōu)化至14.9 mm，同時輸出仿真結(jié)果，優(yōu)化喉口直徑d后文丘里管的流動壓力和速度云圖仿真結(jié)果如圖6所示。整個計算過程，無需人為干涉，無需重新修改模型、重新劃分網(wǎng)格、重新輸入邊界條件等重復(fù)性工作，全程用時15 min左右，極大提高了工作效率。

a)壓力云圖 b)速度云圖圖6 優(yōu)化喉口直徑d的文丘里管流動壓力和速度云圖仿真結(jié)果

4 結(jié)論

采用優(yōu)化設(shè)計軟件Isight、Creo二次開發(fā)工具Protoolkit和CFD仿真分析軟件STAR-CCM+，基于參數(shù)化設(shè)計原理開發(fā)了文丘里管設(shè)計仿真一體化平臺，實現(xiàn)文丘里管自動優(yōu)化設(shè)計，滿足國六柴油機開發(fā)要求。

1)設(shè)計過程實現(xiàn)自動計算及優(yōu)化，無需重新修改模型、重新劃分網(wǎng)格、重新輸入邊界條件等重復(fù)性工作，經(jīng)過多輪迭代優(yōu)化，迅速得到最優(yōu)化的計算結(jié)果，極大提高工作效率。

2)仿真迭代優(yōu)化計算過程無需人為干涉，沒有報錯，證明模型及程序穩(wěn)健，模型及程序的魯棒性滿足文丘里管自動計算優(yōu)化要求。