0 引言
發動機缸內的氣流運動對發動機的燃燒、排放物生成關系極大，是影響發動機性能的重要因素之一。要獲得良好的燃燒過程，最有效的措施是組織空氣運動，而由于氣道結構形狀決定著吸氣渦流質量的好壞，所以氣道對于柴油機性能影響非常顯著，因此設計出符合要求的氣道一直是內燃機設計的一項重要內容。
在氣道設計初期，試驗方法周期長、成本高，而CFD計算可以快速較準確獲得流量系數和缸內流動的詳細數據，為設計與優化提供指導。單相氣體的穩態流動計算較為成熟，收斂性好且可以得到較準確的解，但涉及到與實驗相匹配的各個升程下的穩態計算case多，前后處理量大、過程繁瑣，由于STAR-CCM+的網格生成功能適應性較強，且氣道的組成和幾何構型相對比較簡單，而多面體網格的適應性強，故采用STAR-CCM+進行的穩態氣道模擬分析的設置過程可以模式化，通過程序執行，可大大降低時間成本，提高工作效率。
1 STAR-CCM+的二次開發功能與Java語言簡介
由于STAR-CCM+采用Java語言進行封裝，各種功能、方法、對象的調用均可通過繼承或引用CCM+程序包來實現，所以采用Java語言進行二次開發工作。
1.1 Java語言簡介
Java是一種可以撰寫跨平臺應用軟件的面向對象的程序設計語言，是由Sun Mircrosystem公司于1995年5月推出的Java程序設計語言和Java平臺(即JaSE、JavaEE、JavaME)的總和。Java技術具有卓越的通用性、高效性、平臺移植性與安全性，廣泛應用于個人PC、數據中心、超級計算機、移動電話和互聯網等領域。
Java平臺由Java虛擬機(JVM)和Java應用編程接口(API)構成。Java應用編程接口為此提供了一個獨立于操作系統的標準接口，可分為基本部分和擴展部分。在硬件或者操作系統平臺上安裝Java平臺之后，Java應用程序就可運行。Java平臺已經嵌入了幾乎所有操作系統。這樣Java程序可以只編譯一次，就在各種系統中運行。但是由于其每次執行編譯后的字節碼需要消耗一定的時間，這同時也在一定程度上降低了Java程序的運行效率。
Java語言是面向對象的，提供類、接口和集成等原語，為了簡單起見，只支持類之間的單繼承，但支持接口之間的多繼承，并支持類與接口之間的實現機制。Java語言全面支持動態綁定，是一個純的面向對象程序設計語言。
其擁有的各項優良特性使得Java應用具有非常好的健壯性和可靠性，這將減少應用系統的維護費用。Java對對象技術的全面支持和Java平臺內嵌的API能縮短應用系統的開發時間并降低成本�；谄錇檐浖_發過程提供的有關技術與豐富的類庫，STAR-CCM+采用Java進行封裝，故其二次開發也采用Java語言。
1.2 STAR-CCM+二次開發方式簡介
封裝好的CCM+程序具有圖形操作簡單，界面友好，易學易用等特點，但帶來的后果是程序開放性不強，不支持代碼操作不便于進行高級修改。為了給用戶自定義程序提供一個入口，CCM+采用star.common包中的StarMacro類來實現這一功能，使得用戶可以通過編寫StarMacro的子類來實現在程序當中運行自己編輯的宏命令。
編輯宏文件存在兩種方式，一是通過用戶界面宏操作面板上的“錄制”按鈕自動生成需要的代碼，一種是在對CCM+的類庫以及各個常用類的方法有所了解之后自行引入CCM+的基本類庫在集成編譯環境或者自由編輯器當中自行編寫開發。在宏文件的編輯完成之后，實現自定義程序的運行可以通過兩個方式，一種與之前類似，是通過用戶界面宏操作面板上的“播放”按鈕運行程序，另一種是直接通過操作系統指令（Windows下的cmd命令或者Linux系統在shell下操作）在打開CCM+時指定宏文件執行。
這種做法雖然為用戶自定義程序的實現成為可能，但開發團隊之外的用戶想深入計算模型的修改目前看來仍十分艱難。并且，對于初學者而言，雖然在短時間內掌握錄制和播放宏文件的方法并不困難，但如果要實現大量相似但并不完全相同模型（特別是幾何模型有差異）的設置、網格劃分、計算、后處理過程，通過簡單的錄制和播放宏或難以理想的實現，仍需逐一在宏文件當中尋找并修改參數，帶來仍然不小的工作量和出錯的概率。
所以針對不需要對計算模型進行深入修改，但工作過程高度相似化流程化的計算工作，通過合理安排宏文件的運行組織結構，并將參數提取成變量從方法當中抽取出來，將方法重構就可以實現計算過程的集成。這樣便可大大減少相似工作的重復工作量和出錯幾率，還可以采用Java語言為其編寫圖形化的界面交互程序，使得不了解CCM+類庫方法調用的用戶可以方便地修改某一類型案例的可變參數進行初步的計算分析。
2 進氣道穩態模擬計算的流程
進氣道流動特性一般指流通能力(或阻力特性)和產生渦流的能力。目前，評價氣道流通能力的試驗方法仍采用在穩流試驗臺上測得的流量系數法。而評價進氣道產生渦流能力的方法除在穩流試驗臺上測量渦流比外，還有其它一些方法。在穩流試驗臺上評價進氣道進氣渦流的方法有兩種：一是風速儀法；二是氣流校直測矩儀法。氣道內的流動狀態用壓差來度量，氣體流量用標準流量計測量。為了便于對不同形狀和尺寸的氣道流動性能進行比較，可以用無量綱流量系數評價不同氣門升程下氣道的流通能力，用無量綱渦流比來評價不同氣門升程下氣道形成缸內渦流的能力。世界上幾家著名的內燃機咨詢公司（如Ricardo、AVL等）分別定義了各自的評價參數。
針對穩態實驗進行的模擬計算，操作過程可以分為前處理、計算以及后處理這樣三個過程。其中，前處理主要包括導入幾何，劃分計算區域，生成網格，選擇并設置模型、求解器參數等步驟，后處理則涉及截面的截取，云圖、矢量圖的生成以及各種變量積分的計算和評價指標的計算。
在前處理當中，根據流體區域表面所屬的零件的不同，本文將氣道幾何分為穩壓箱、氣道、座圈、氣門、氣缸幾個區域，分別為其賦予標準名稱，之后將氣門從坐標原點沿氣缸的對稱軸向下移動當前升程對應的位移，并對各個零部件的外表面進行布林運算，生成流動的幾何區域。
圖1 計算分析工作流程
生成的流體區域按照幾何表面區分邊界，之后設置面網格模型與基準尺寸，將氣道與座圈邊界進行局部加密，設置并重新劃分面網格完畢后選取體網格模型繼續生成多面體網格，在完成網格劃分之后選取物理模型以及進行相應的求解器參數設置。最后，根據穩態氣道實驗的需要進行后處理，例如按要求選取截面查看云圖或者矢量圖，統一云圖的上下界便于其比對，計算該升程下的流量系數以及無因次渦流量，最終根據得到的結果進行分析并撰寫報告。
3 用戶界面的編寫與程序的特點
針對以上工作流程可以將程序劃分為不同的模塊分別進行編寫，在實現了相應功能之后，為了方便參數的修改，通常會把需要修改的參數以變量的形式提取出來，如果要使用戶可以更直觀的修改這些變量，可以編寫用戶輸入界面將其值寫入相應的屬性文件當中，然后在子類或者宏程序當中對屬性文件進行讀取，進而采用用戶指定的值來執行各個命令和方法。
根據功能區域的不同，本文的用戶面板由四個選項卡組成，分別對工作及保存路徑、氣門升程、計算參數進行設置，“基本設置”面板當中發動機及環境信息的設置與流量系數與渦流的計算有關，后處理面板實現的是具體結果的實時調用，如圖所示。本程序允許異地調用和保存幾何文件及結果文件，待處理的幾何文件不必保存在工作目錄下，計算得到的結果文件(.sim后綴文件)和各截面云圖矢量圖也可以保存在用戶指定的目錄下。
通過“基本設置”面板當中的“檢查”按鈕和后處理面板當中的“顯示當前”按鈕，可以顯示指定幾何保存路徑下各幾何零件的組裝效果以及當前選擇的結果截圖。在后處理面板上執行了“處理結果并保存”按鈕之后，可以通過“流量系數”“無因次渦流量”兩個按鈕來查看所選case的相應結果，所有已計算升程下的流量系數和無因次渦量結果均保存在工作目錄下的“...\Cf.properties”以及“...\swirlFile.properties”兩個文件當中，用戶可以根據需要使用。
圖2 幾何模型檢查窗口
程序對于運行環境有一定要求，首先用戶的PC上需要安裝STAR-CCM+8.02和任意版本的STAR-View+，由于CCM+不同版本的程序包之間存在一定差異，所以如果安裝的版本不是8.02版本，程序有可能會發生一些錯誤。由于程序是采用Java語言編寫，輸出的程序包又并未添加自動檢測并安裝JRE(Java Runtime Environment)的功能，所以程序必須在安裝了JRE的環境下運行。
圖3 程序的運行界面
4 具體案例分析
采用本程序，本文對某發動機氣道進行了模擬計算分析。
4.1 模型設置及網格情況
流體域的幾何模型如圖所示。氣道前端是穩壓箱，氣缸長度取1.75倍的缸徑。采用多面體網格方式進行劃分，基本參照尺寸為2mm，氣門座以及氣門等有較小結構處按照程序當中根據經驗設定的比例采取了網格加密措施，網格總數為582145，網格質量全部合格，良好質量網格占98.848%。
圖4 計算網格
程序根據經驗選定的各壁面條件以及計算模型等無法通過用戶面板修改，故仍舊按照統一設置，這樣也有利于模型的統一。進出口壓差也按照文本框的默認值4100Pa進行設置，溫度在基本設置面板設置為290.00K。
4.2 計算結果及分析
圖為升程分別為4mm和7mm時的壓力場。進氣行程中進氣門打開，活塞下行，在氣缸內形成負壓。氣體進入氣道后，由于氣道截面逐漸減小，氣道內的阻力逐漸增加，順著氣流方向，氣道內壓力逐漸減小，氣道出口即氣門座附近壓力最小，而在該處氣道內阻力最大，壓力損失也大。
圖為4mm升程和7mm升程時氣道和氣缸內的速度云圖。喉口處氣體被壓縮，速度急劇增加，而氣體進入氣缸后，由于截面積增大，流速逐漸減小。
圖5 不同升程下的壓力分布
圖6 不同氣門升程下的速度矢量分布
圖7 不同氣門升程下的湍動能分布
在小升程時，氣門背側沒有流速突降區域，大升程下存在這一區域，這一回流區域是邊界層分離所造成的。通過對比不同氣門升程下的速度分布可以發現，當氣門升程較大時，缸內的氣流速度較大，但喉口處流速低于小氣門升程時，反之亦然。并且由于氣門中心與氣缸中心在水平方向的相對位置不一致，導致不同側的進氣流速不同，這種不平衡也應當會影響到缸內渦流的形成，決定渦流的形狀與漩渦中心的位置。
從圖7不同升程下的湍能強度對比中可以發現，隨著進氣量的增加，無論是在氣道內還是氣缸內，湍動能都在增加。圖8給出了該氣道各氣門升程下的流量系數曲線，從曲線走勢來看，隨著氣門開度的增加，流量增加的幅度基本呈線性關系增加。
圖8 各個氣門升程下的流量系數
5 結論
為了分析示例當中的氣道，本文共計算了六個升程的穩態流動情況，產生了240個圖片和可視化文件，而其中很大量用于對比的圖片和文件最終雖然不會被報告采用卻仍需要制作以備篩選，程序的使用節約了大量工時，并減少了設置偏差導致計算條件不統一而返工的情況。
6 展望
在前處理方面，目前程序只能處理兩氣門氣道變十氣門升程的情況，如果將四氣門的宏程序進行方法的統一和變量的提取，通過增加或替換操作類的方式可以大大增強程序的適用性；計算過程方面，改進新增Simulation類的實例化方式（目前是通過文件地址實例化）可以拓展并行計算和計算過程的控制性，這與增加可操作氣門升程的功能也是相關的；后處理方面，通過Open Office程序包的引入和報告模板的建立，可以進一步實現更標準更終端的結果輸出，目前階段的程序對于用戶提取結果的靈活性實現得也十分有限；考慮到同一單位不同產品的對比分析在統一計算模型的前提下進行更為合理，計算模型的選取在今后的開發中仍不建議對用戶開放。
核心關注：拓步ERP系統平臺是覆蓋了眾多的業務領域、行業應用，蘊涵了豐富的ERP管理思想，集成了ERP軟件業務管理理念，功能涉及供應鏈、成本、制造、CRM、HR等眾多業務領域的管理，全面涵蓋了企業關注ERP管理系統的核心領域，是眾多中小企業信息化建設首選的ERP管理軟件信賴品牌。
轉載請注明出處：拓步ERP資訊網http://m.lukmueng.com/
相關文章
