李振等

摘 要：通過對罐體進行全尺寸建模模擬，對比研究充液比、密度和黏度對液體晃動的影響。結(jié)果表明：罐車橫向運動時，液體晃動水擊內(nèi)壁壓力隨著充液比的增加先增大后減小，充液比為0.4時達到最大；且水擊內(nèi)壁壓力隨著液體密度的增加而增大，隨著液體黏度的增加而減小。當防波板面積大于60%時，最大水擊內(nèi)壁壓力顯著降低，防波效果明顯。

關(guān)鍵詞：罐車；液體晃動；多相流；VOF模型

Simulation of Liquid Slosh in Lateral Movement of Tank Truck

Li Zhen，Xie Moyi，Zhu Bing，Zhao Yunyun

（College of Mechanical and Automobile Engineering， Chongqing Jiao Tong University， Chongqing 400074， China）

Abstract：Based on full-scale modeling simulation of tank， the VOF model based on multi-phase flow， contrast research liquid ratio， the influence of the density and viscosity of liquid sloshing.Results show that the tank car lateral movement， the liquid sloshing water hammer internal pressure with the increase of liquid than first increases then decreases， maximum when liquid ratio 0.4；And the inner wall of the water hammer pressure increased with the increase of liquid density increases， decreases with the increase of liquid viscosity.When washboard area greater than 60%， the floods hit wall pressure significantly reduced， the wave effect is obvious.

Key words：tank truck； liquid sloshing； multi-phase flow； volume of fliud （VOF） model

罐車部分載荷做橫向運動時，罐內(nèi)液體晃動對罐體產(chǎn)生沖擊載荷，使罐車重心偏移，易造成車輛側(cè)滑、側(cè)傾從而引發(fā)交通事故。

1劉奎和康寧[1]通過實驗與數(shù)值模擬相結(jié)合的方法，得出防波板面積大于橫截面的40%時，增加防波板面積能顯著改善罐體受力的結(jié)論；孫麗娜等人[2]對罐車緊急制動時的安全性問題進行了研究，得出罐車液體充裝率、密度和黏度對罐體水擊壓強影響。文獻[3]采用VOF方法研究了不同充液比、不同防波板面積及不同罐體半徑情況下，罐車制動時罐體內(nèi)液體晃動動力特性。

本文針對罐車轉(zhuǎn)向運動時液體的橫向晃動問題，運用VOF方法對罐體內(nèi)液體的晃動動力學特性進行數(shù)值模擬，研究不同充液比、密度、黏度和不同防波板面積對罐體水擊壓強的影響，以期為提高罐車側(cè)向穩(wěn)定性提供設(shè)計參考。

1 液體晃動水擊過程模擬與結(jié)果分析

1.1 罐體模型與邊界條件確定

本文以某罐車罐體為研究對象，進行全尺寸建模，罐體長度6000mm，端面為圓，半徑R=1000mm，罐體壁厚度均為4mm。在進行數(shù)值模擬分析時，不考慮罐體變形，設(shè)定其內(nèi)部介質(zhì)為水和空氣。采用Fluent VOF歐拉兩相流模型，初始時刻保持靜止，上方為空氣，下方為水，t=0時，罐體受到橫向激勵，加速度ɑ=0.8g，持續(xù)時間為4s。

1.2 過程參數(shù)和液體性能參數(shù)對晃動影響

1.2.1 充液比對罐壁液體晃動水擊壓強影響

充液比是影響罐內(nèi)液體晃動的重要參數(shù)之一。在本文中定義充液比k=h/2R，h為液體深度。通過改變罐內(nèi)液體的充液比，研究罐車轉(zhuǎn)向過程中罐內(nèi)液體充液比對液體晃動的影響。充液比分別取0.1，0.2，0.3，0.4，0.5，0.6，0.7，0.8，0.9，共九種工況。

充液比對液體晃動水擊內(nèi)壁壓力影響的模擬結(jié)果如圖1所示。研究結(jié)果表明，橫向運動過程中液體晃動水擊內(nèi)壁壓力先隨充液比增加而增大，充液比大于0.4后，液體晃動水擊內(nèi)壁壓力隨充液比增加而減小，充液比為0.4時，液體晃動水擊內(nèi)壁壓力達到最大。

1.2.2 罐內(nèi)介質(zhì)特性對罐壁壓強影響

液體的主要性能參數(shù)是密度和黏度。本文通過改變液體的密度和黏度，研究罐內(nèi)介質(zhì)特性對罐壁水擊壓力的影響。

液體密度分別取500kg/m3，750kg/m3，1000kg/m3，1250kg/m3，1500kg/m3五種情況。數(shù)值模擬結(jié)果如圖2所示，液體晃動水擊內(nèi)壁壓力在t=2.3s時達到最大值，且水擊內(nèi)壁壓力隨著液體密度的增大而增大。

液體黏度分別取0.0001kg/（ms），0.001kg /（ms），0.01kg /（ms），0.1kg/（ms），1kg /（ms）五種情況。數(shù)值模擬結(jié)果如圖3所示，液體晃動水擊內(nèi)壁壓力在t=2.6s時達到最大值，且水擊內(nèi)壁壓力隨著液體黏度的增大而減小。但液體黏度對水擊內(nèi)壁壓力的影響遠小于密度。

2 不同防波板面積對罐壁沖擊壓力的影響

文獻[4]研究表明布置一塊較大面積的防波板即可達到較好的防波效果，但沒有具體分析防波板面積對罐壁沖擊壓力的影響。

《液化氣體運輸車》中規(guī)定，“每個防波板的有效面積應大于罐體橫斷面積的40%，防波板的安裝位置，應使上部弓形面積小于罐體橫斷面積的20%”，參考以上對橫向防波板的尺寸規(guī)定，設(shè)置縱向防波板高度，示意圖如圖4所示，S1為縱截面上部空間面積，S2為縱向防波板面積，S3為縱截面下部空間面積。令a=S1/（S1+S2+S3）；b=S2/（S1+S2+S3），應保證a≤0.2；b≥0.4。

在計算不同面積防波板的影響時，為研究方便，令a為定值，b為變量，本文統(tǒng)一取a=0.2，b的取值范圍0～0.8。側(cè)向加速度ɑ=0.8g，充液比k=0.5，模擬結(jié)果如圖5所示。

隨著防波板面積的增加，液體晃動水擊內(nèi)壁壓力先增大后減小，當b=0.5時，水擊內(nèi)壁壓力達到最大值，這是因為液體晃動時受到防波板阻礙后回流，與罐壁沖擊的次數(shù)增多，使得液體在較為狹小的空間內(nèi)更快的沖擊罐壁，從而使得水擊內(nèi)壁壓力劇烈增大。b=0.6時，水擊內(nèi)壁壓力大幅度降低，且b=0.6～0.8之間，圖5中曲線斜率較小，b=0.6時即可有效抑制液體晃動，取得較好的防波效果。

3 結(jié)論

1） 罐車橫向運動過程中，水擊內(nèi)壁壓力先隨著充液比的增加而增大，當充液比大于0.4后，液體晃動水擊內(nèi)壁壓力隨著充液比的增加而減小，充液比為0.4時，液體晃動水擊內(nèi)壁壓力達到最大。

2） 液體晃動水擊內(nèi)壁壓力隨著液體密度的增加而增大，而隨著液體黏度的增加而減小，且液體黏度影響程度遠小于密度影響程度。

3） 布置防波板后，能夠顯著降低水擊內(nèi)壁壓力最大值，當防波板面積大于60%后，最大水擊內(nèi)壁壓力顯著降低，防波效果明顯。

參考文獻

1 劉奎，康寧.罐車制動時液體晃動的仿真分析[J].北京航空航天大學學報，2009，35 （7）：799-803.

2 孫麗娜，周國發(fā).罐式集裝箱液體晃動過程數(shù)值模擬研究[J].振動與沖擊，2012，31（22）：147-150.

3 Kang N，Liu K. Influence of baffle position on liquid sloshing during braking and turning of a tank truck[J]. Journal of Zhejiang University： Science A，2010，11（5）：317-324.

4 趙樹恩，趙靈鶴.汽車罐車橫向運動液體晃動動力學特性模擬[J].應用數(shù)學和力學，2014，35（11）：1259-1270.

作者簡介

李振（1988-），男，碩士研究生，研究方向：液罐車轉(zhuǎn)向側(cè)向穩(wěn)定性