王國(guó)慶 戴汝泉 張竹林 王金波

山東交通學(xué)院 山東濟(jì)南 250000

1 前言

液體運(yùn)輸車(chē)在行駛過(guò)程中，由于轉(zhuǎn)彎、制動(dòng)等運(yùn)動(dòng)狀態(tài)的改變，罐內(nèi)液體受到慣性力的作用，出現(xiàn)罐體和充裝介質(zhì)的流固耦合效應(yīng)。在這種流固耦合效應(yīng)的作用下，罐內(nèi)介質(zhì)的晃動(dòng)對(duì)罐體產(chǎn)生沖擊力，嚴(yán)重時(shí)會(huì)造成罐體的破壞[1]。當(dāng)前罐車(chē)制造企業(yè)采用不同型式的防波板，將罐內(nèi)流體區(qū)域分為多個(gè)小艙室，罐體內(nèi)部流體的沖擊由防波板分擔(dān)，從而降低罐體損壞的可能性。

本文對(duì)安裝不同型式防波板的罐體進(jìn)行流固耦合分析，研究了當(dāng)前罐車(chē)企業(yè)所采用的不同型式防波板的防晃效果。

2 計(jì)算模型和計(jì)算方法

2.1 罐體模型和網(wǎng)格劃分

本文以某型號(hào)液體運(yùn)輸車(chē)的罐體為研究對(duì)象，按一定比例縮放建立的罐體模型如圖1所示，簡(jiǎn)化的罐體總長(zhǎng)為650 mm，罐體直徑為218 mm。圖1所示的罐體由前后封頭、3個(gè)防波板和筒體組成。其中防波板將罐體分割為4個(gè)區(qū)域，各區(qū)域之間通過(guò)5個(gè)通液孔連接。圖1所示的模型是安裝了曲面防波板的罐體，其他兩種模型分別為不安裝防波板和安裝直面防波板的罐體，三種罐體的前后封頭和筒體尺寸完全一樣，兩種不同型式防波板通液孔開(kāi)孔位置和尺寸大小完全一致，在罐體安裝位置保持一致。計(jì)算時(shí)，坐標(biāo)系原點(diǎn)位于后封頭端面中心點(diǎn)位置，選定罐體前進(jìn)方向?yàn)閄軸的正方向且平行于地面，垂直于地面且向上的方向?yàn)閅軸正方向，平行于地面且向右的方向?yàn)閆軸正方向。

圖1 罐體模型

圖2 流體域網(wǎng)格

網(wǎng)格劃分包括流體域網(wǎng)格和固體域網(wǎng)格的劃分，流體域網(wǎng)格采用CFD流體單元，固體域網(wǎng)格采用結(jié)構(gòu)單元。無(wú)防波板罐體填充的流體定義為1號(hào)流體域，安裝曲面防波板和直面防波板的罐體填充的流體分別定義為2號(hào)流體域和3號(hào)流體域。1號(hào)流體域由于無(wú)防波板分割因而形狀規(guī)整，故采用sweep方法生成六面體網(wǎng)格。2號(hào)和3號(hào)流體域形狀不規(guī)整，通液孔區(qū)域采用六面體網(wǎng)格，其他區(qū)域以四面體網(wǎng)格為主。固體域包括前后封頭、防波板和筒體，固體域進(jìn)行抽中面處理，生成殼單元網(wǎng)格。流體域生成的網(wǎng)格如圖2所示，圖2所生成的流體域網(wǎng)格，從上至下依次為1號(hào)、2號(hào)和3號(hào)流體域網(wǎng)格。

2.2 流固耦合分析

本文選用的罐體材質(zhì)為鋁合金，罐內(nèi)填充的流體介質(zhì)為水，采用的罐體模型在緊急制動(dòng)工況下固體域受到流體域沖擊變形很小，故忽略固體域計(jì)算結(jié)果對(duì)流體域運(yùn)動(dòng)的影響，只把CFD分析計(jì)算的結(jié)果通過(guò)耦合交界面?zhèn)鬟f給固體域，從而得到固體域罐體變形圖。

2.3 氣液兩相流計(jì)算模型

多相流采用VOF模型[2]，VOF模型通過(guò)求解一套運(yùn)動(dòng)方程組和處理穿過(guò)區(qū)域的每一流體的體積分?jǐn)?shù)可以精確模擬氣液分界面的變化情況和液體晃動(dòng)過(guò)程。VOF模型求解的一套運(yùn)動(dòng)方程組，包括質(zhì)量守恒方程和動(dòng)量方程。

式中，t為時(shí)間，ρ為流體密度→， 為流體速度矢量，τf為剪切u力張力，ff為體積力矢量。

2.4 求解器參數(shù)設(shè)置

湍流模型采用標(biāo)準(zhǔn)k-epsilon模型，k-epsilon模型是基于湍動(dòng)能k及湍流耗散率epsilon的一種數(shù)值模型，具有可靠、收斂性好和內(nèi)存需求低的特點(diǎn)，流體域基于壓力求解，鑒于氣液兩相交界面是隨著運(yùn)算時(shí)間變化的，因此采用瞬態(tài)求解。速度和壓力的耦合處理采用PISO算法，壓力修正方程離散格式采用Body Force Weighted格式。

2.5 邊界條件

流固耦合的邊界條件設(shè)置分為固體域和流體域兩部分，固體域的邊界條件為：給定X方向的制動(dòng)減速度5.08 m/s2和Y軸負(fù)方向的重力加速度9.81 m/s2，約束位置為罐體與車(chē)架連接的腹板位置[3]。流體域的邊界條件為：給定X方向的制動(dòng)減速度5.08m/s2和Y軸負(fù)方向的重力加速度9.81 m/s2[4]，給定X方向的初速度為50 km/h,不考慮傳熱影響，設(shè)置壓力參考點(diǎn)，定義的耦合面為流體區(qū)域的外表面。邊界條件的制動(dòng)減速度5.08 m/s2，是根據(jù)行駛初速度50 km/h和最小制動(dòng)距離19 m的規(guī)定求出的。本文的主要目的是研究防波板的防晃效果，故流體域和固體域的瞬態(tài)仿真時(shí)間選取液體晃動(dòng)比較劇烈的時(shí)間段0 ～1.2 s，非定常時(shí)間步長(zhǎng)為0.001 s，流體域時(shí)間步收斂殘差數(shù)量級(jí)達(dá)到10-4，自動(dòng)保存時(shí)間步長(zhǎng)為100。

3 防波板的防晃效果分析

3.1 氣液兩相圖

氣液兩相分布圖能夠清晰地看出整個(gè)仿真過(guò)程中流體運(yùn)動(dòng)的狀況，氣液兩相圖選取的平面為流體域縱向中間軸平面，初始?xì)庖航唤缑嫫叫杏赬Z平面。三種流體域氣液兩相圖如圖3所示，由上至下分別為1號(hào)、2號(hào)和3號(hào)流體域氣液兩相圖，由左至右分別對(duì)應(yīng)的時(shí)間點(diǎn)為0 s、0.4 s、0.8 s和1.2 s。

圖3 氣液兩相圖

以1號(hào)流體域氣液兩相圖為參考對(duì)象，2號(hào)和3號(hào)流體域由于防波板的作用，液體運(yùn)動(dòng)明顯受到防波板的阻礙。1號(hào)流體域的液體在0.8 s時(shí)絕大部分位于罐體的前半部區(qū)域，2號(hào)和3號(hào)流體域的液體在0.8 s時(shí)仍有一部分液體位于罐體的后半部區(qū)域。1號(hào)流體域在1.2 s內(nèi)液體碰到前封頭出現(xiàn)了液體回流，2號(hào)和3號(hào)流體域的液體在1.2 s內(nèi)仍未出現(xiàn)明顯的液體回流。

3.2 流體最大壓力值曲線

根據(jù)不同時(shí)刻的流體壓力，繪制流體最大壓力值曲線如圖4所示，圖4除去了0.1 s、0.2 s、0.3 s的流體最大壓力值，因?yàn)檫@些時(shí)間點(diǎn)的流體壓力最大值不符合仿真要求，從圖4可以看出1號(hào)流體域流體最大壓力值高于2號(hào)和3號(hào)流體域，2號(hào)流體域流體最大壓力值在大多數(shù)時(shí)間點(diǎn)略高于3號(hào)流體域。在整個(gè)仿真時(shí)間歷程中，1號(hào)、2號(hào)和3號(hào)流體域流體壓力峰值分別出現(xiàn)在時(shí)間點(diǎn)1 s、1.2 s和1.1 s。

圖4 流體最大壓力值曲線

3.3 罐體結(jié)構(gòu)分析

罐體在0～1.2 s內(nèi)所受到的載荷包括自身的重力、流體的沖擊力和制動(dòng)力，罐體結(jié)構(gòu)約束的位置為車(chē)架腹板和罐體接觸區(qū)域。三種罐體結(jié)構(gòu)分析選取的時(shí)間點(diǎn)分別為1 s、1.2 s和1.1 s，通過(guò)對(duì)比三種罐體變形值，得出不同型式防波板的防晃效

圖5 無(wú)防波板罐體變形圖

圖5為1 s時(shí)無(wú)防波板的罐體變形圖，罐體最大變形值出現(xiàn)在罐體中部頂端區(qū)域，變形區(qū)域較為集中，最大變形值為6.3829×10-7m。圖6為1.2 s時(shí)曲面防波板的罐體變形圖，罐體最大變形值出現(xiàn)在罐體前封頭底部區(qū)域，最大變形值為7.7657×10-8m，防波板1最大變形值出現(xiàn)在中間通液孔附近區(qū)域，最大變形值為5.4825×10-8m。圖7為1.1 s時(shí)直面防波板的罐體變形圖，罐體最大變形值出現(xiàn)在罐體前封頭區(qū)域，最大變形值為7.1459×10-8m，防波板1最大變形值出現(xiàn)在中間通液孔區(qū)域，最大變形值為2.1627×10-7m。

圖6 曲面防波板罐體變形圖

圖7 直面防波板罐體變形圖

4 結(jié)語(yǔ)

本文以液體運(yùn)輸車(chē)罐體為研究對(duì)象，以不安裝防波板的罐體作對(duì)比，對(duì)安裝曲面防波板和直面防波板的罐體在制動(dòng)工況下罐內(nèi)液體晃動(dòng)進(jìn)行流固耦合分析，得到了流體域氣液兩相圖和流體最大壓力值曲線，對(duì)三種罐體結(jié)構(gòu)分析選取的時(shí)間點(diǎn)分別為1 s、1.2 s和1.1 s，對(duì)比三種罐體變形值，進(jìn)而分析出不同型式防波板的防晃效果。從中獲得的結(jié)論如下：

a. 無(wú)防波板的罐體在1.2 s仿真過(guò)程中，液體已完成對(duì)罐體壁面的沖擊出現(xiàn)了回流現(xiàn)象，罐體最大變形區(qū)域較為集中且分布在罐體中部頂端區(qū)域。

b. 對(duì)比三種罐體的變形值可以看出，曲面防波板和直面防波板的罐體最大變形值均小于無(wú)防波板的罐體，兩種型式的防波板起到了明顯的防晃作用。直面防波板比曲面防波板能更有效地減緩液體對(duì)罐體的沖擊力，但直面防波板最大變形值也大于曲面防波板。

由于計(jì)算機(jī)資源的限制[5]，本文未采用實(shí)車(chē)罐體模型進(jìn)行流固耦合分析，關(guān)于罐體在其他充液比下的液體晃動(dòng)仍需進(jìn)一步研究，進(jìn)而更全面地考察不同型式防波板的防晃效果。

[1] 張韜,柯龍燕,白光建.罐車(chē)制動(dòng)工況液體晃動(dòng)的流固耦合分析[J].專用汽車(chē),2015(08):88-91.

[2]李清.貯箱液體晃動(dòng)及穩(wěn)定性分析[D],南京航空航天大學(xué),2017.

[3] 王次安.基于流固耦合的液體罐裝車(chē)瞬態(tài)應(yīng)力分析及優(yōu)化設(shè)計(jì)[D].合肥工業(yè)大學(xué),2013.

[4] Martin Mo?ilan,Milan ?mindák,Peter Pechá?,Peter Weis.CFD Simulation of Hydraulic Tank[J]. Procedia Engineering,2017,192:46-49

[5] 劉小民,王星,許運(yùn)賓.運(yùn)動(dòng)罐內(nèi)液體晃動(dòng)的雙向流固耦合數(shù)值分析[J].西安交通大學(xué)學(xué)報(bào),2012.46(5):120-126.