王林　朱天軍　胡偉　栗圣博

摘要：流固耦合的作用是影響危化品液罐車的運輸安全性、可靠性的一個重要領(lǐng)域。罐體內(nèi)液體的晃動所產(chǎn)生的力對車輛運動學(xué)和車輛穩(wěn)定性有著很大的影響：液罐車在轉(zhuǎn)向過程中，液體的橫向動力學(xué)會影響到車輛的行駛的穩(wěn)定性。而縱向動力學(xué)會造成制動力的損失和加速方向控制力不穩(wěn)定。在本文中，論文借助理論分析、等效模型和仿真研究等方法，對縱向運動學(xué)進行了深入研究，通過數(shù)學(xué)模型的建立，從原理數(shù)值上論證了新型防波板的設(shè)計合理性。對設(shè)計的不同防波板進行建模與仿真，通過定點分析以及對比，得到了令我們滿意的仿真結(jié)果，驗證了防波板模型的合理性和新型防波板在?；芬汗捃囍型茝V的可行性。本文為設(shè)計新型防波板以及研究不同防波板對液罐車行駛穩(wěn)定性影響提供了一定的理論依據(jù)。

關(guān)鍵詞：危化品油罐車;穩(wěn)定性;模型

1.1油罐車動力學(xué)模型

由于低控性的限制，重型車輛在規(guī)避機動時以特殊的方式發(fā)生傾翻現(xiàn)象。事實上，運載剛性貨物的卡車的重心太高，與SUV或乘用車相比，這導(dǎo)致了縱向和橫向的巨大載荷轉(zhuǎn)移。準(zhǔn)靜態(tài)側(cè)翻閥值（SRT）是指在穩(wěn)態(tài)轉(zhuǎn)彎時，內(nèi)曲線側(cè)輪的升力引起的橫向加速度的量度，用重力加速度的分?jǐn)?shù)表示。已發(fā)現(xiàn)準(zhǔn)靜側(cè)翻閥值適合于預(yù)測重型車輛的側(cè)傾穩(wěn)定性特性。忽略輪胎和懸掛的彈性變形，對靜態(tài)側(cè)翻閾值的一個粗略估計是由靜態(tài)穩(wěn)定因子（SSF）比給出的：。通過靜態(tài)穩(wěn)定因子（SSF）方程（1）中是軌道寬度，是重心的高度;明確地指出了重心高度對罐車穩(wěn)定性的影響;事實上，增加的結(jié)果就是使得靜態(tài)穩(wěn)定因子比減少。在許多研究中，已顯示側(cè)翻碰撞頻率與側(cè)翻穩(wěn)定極限之間的關(guān)系。為降低與翻車事故有關(guān)的公路安全風(fēng)險，已提出了最低0.35 g的商業(yè)運營標(biāo)準(zhǔn)。

1.2 ?罐體幾何模型的建立

本文研究的對象是基于勒洛三角原理，在傳統(tǒng)的方形和圓形油罐車的橫截面基礎(chǔ)上改善而來的，而設(shè)計出的上圓下方的方圓形油罐體幾何截面模型。具有儲量大、質(zhì)量重心低的特點。

?液化氣體運輸車?中規(guī)定：“罐體內(nèi)應(yīng)設(shè)置防波板，每個防波板間距不得超過1.75m，每個防波擋流板的有效面積不得低于罐體橫斷截面的40%”。根據(jù)我國的相關(guān)設(shè)計規(guī)范，本文設(shè)計的罐體事以我國通用標(biāo)準(zhǔn)設(shè)計而來的。罐體的幾何尺寸：罐體全長為6880mm，高度為2030，寬度為2600，各擋流板間隔為1720mm，為了簡化模型，本設(shè)計壁厚和擋流板厚度均為10mm。罐體、封頭、擋流板的材料均采用強度較高的>400MPa的低合金高強度結(jié)構(gòu)鋼。

1.3 ?等效晃動幾何模型的建立

因為液體晃動本身是一個復(fù)雜且無法準(zhǔn)確定義的運動狀態(tài)，所以我們需要簡化模型，并建立油罐車晃動的等效力學(xué)準(zhǔn)靜態(tài)模型，通過準(zhǔn)靜態(tài)方法分析不同類型的擋流板上的壓力、應(yīng)力等，來研究擋流板對油罐車穩(wěn)定性的影響。

在圖1中，是側(cè)面的加速;是垂直輪胎負荷，= 1，2;是車輛的重量，是相對于軌道的重心的橫向運動。

準(zhǔn)靜態(tài)模型是用來預(yù)測運動容器內(nèi)液體液面穩(wěn)態(tài)位置的常用方法。在這種方法中，可以對不同激勵下的液體，進行質(zhì)心位置的標(biāo)定。結(jié)果表明，隨意晃動的液面可以用直線或面來代替。根據(jù)滾轉(zhuǎn)平面的橫向加速度或傾斜平面的縱向加速度，可以計算出自由液面的位置。然后，在車輛轉(zhuǎn)向時，分析穩(wěn)態(tài)貨物在動力學(xué)性能下的載荷轉(zhuǎn)移情況。在穩(wěn)定轉(zhuǎn)向時，假定不可壓縮的油罐車均勻流體在非慣性參考系中動量守恒：

（1）

;（2）

其中，是流體密度和壓力;和分別是重力和油罐的自轉(zhuǎn)速度，是非慣性參照系的旋轉(zhuǎn)軸線和原點，即儲罐的中底點之間的距離。

這組方程為小滾角提供了一個拋物面自由面。然而，對于半徑較大的道路集裝箱相比于集裝箱的寬度而言來說，術(shù)語cos可以忽略不計。這與羅伯特等人提出的罐車液面自由度的評定可以忽略這種近似誤差結(jié)論相爭議。在方程（1.1）和（1.2）中引入橫向加速度：

（3）

（4）

通過自由表面上的p=0，可以得到自由表面下面的方程：

（5）

穩(wěn)態(tài)轉(zhuǎn)動下二維曲面模型導(dǎo)出自由液面方程：（6）

假定小滾轉(zhuǎn)角的自由表面梯度可以表示為：（7）

因此，滾平面內(nèi)的自由面是一個與它的初始位置相對應(yīng)的直線。類似的方程可以得到二維平面間距恒定縱向加速：（8）

其中P是表示在俯仰平面液體自由表面的坡度，是表示縱向加速度而P是罐車的傾斜角?；谧杂汕娴奶荻取⒁后w填充水平和罐的幾何形狀，可以計算液體體積在縱向、橫向和垂直方向上的平移。曲面平面模型中的液體重心坐標(biāo)可以得到如下：

（9）

Ω是流體取值范圍，對于傾斜平面模型，可以計算液體質(zhì)量中心坐標(biāo)：

（10）

準(zhǔn)靜態(tài)晃動模型是探索油罐車液體運動狀態(tài)下穩(wěn)態(tài)負荷變化的一種有效的方法。

2.1 ?三維模型的建立

因為罐體和擋流板的三維模型較復(fù)雜，所以本文所需要的三維模型是在CATIA軟件進行建模的。因為本文研究的是不同擋流板對油罐車穩(wěn)定性的影響，所以罐體簡化模型的尺寸和截面都是統(tǒng)一的。如圖所示：

2.2 ?網(wǎng)格劃分及網(wǎng)格質(zhì)量評價

網(wǎng)格劃分是流體動力學(xué)仿真分析最重要的步驟之一，網(wǎng)格劃分的密度和質(zhì)量直接影響到了流體仿真結(jié)果的精確度和準(zhǔn)確性。在網(wǎng)格形狀參數(shù)中，我們選擇六棱柱和四邊形相結(jié)合的網(wǎng)格形狀，在網(wǎng)格關(guān)聯(lián)質(zhì)量我們選擇高等，高級尺寸功能我們選擇按照距離和曲率進行劃分網(wǎng)格。因為本文我們研究的對象主要是防波板對油罐車穩(wěn)定性的影響，關(guān)于膨脹層的設(shè)置，我們設(shè)置的對象是在三個防波板上，膨脹層數(shù)為5層，過渡方式選擇平滑過渡。生成網(wǎng)格后檢查網(wǎng)格質(zhì)量，網(wǎng)格質(zhì)量合適回到Workbench中更新網(wǎng)格。

3 ?CFD（計算流體動力學(xué)）仿真和分析

在Workbench搭建的仿真流程，設(shè)置劃分好網(wǎng)格就可以正式進入fluent板塊里面進行相關(guān)步驟的設(shè)置。因為我們研究的是液體晃動瞬態(tài)的，所以我們的模型選擇湍流，本文也將已經(jīng)到液體體積的變化，要求設(shè)置為VOF。液體材料我們選擇煤油，主相為空氣，次相為我們選擇的煤油。

因為我們還將研究不同填充率對擋流板對油罐車穩(wěn)定性的影響，因為本文模型較正常油罐車模型復(fù)雜些，所以本文通過設(shè)置液面在罐體里的高度來簡化煤油的填充率。為了簡化罐體的受力境況，在罐體的底部看做是地面并施加一個與之平行的縱向力，設(shè)置受力20秒，前12秒加速到10m/s，后8秒緊急制動到0m/s。因為與液體相作用的是內(nèi)表面，所以我們先抑制外壁再激活流固耦合界面。其它條件均是理想情況，而且我們本次仿真不考慮液體的粘度對晃動的影響。該模型采用三種不同的橫向擋流板，在理想的縱向加速度和重力加速度場的組合下，進行仿真。

在進行有限體元仿真分析中，施加縱向加速度下的瞬時晃動力。研究表明，不同的填充率，對擋流板的影響也有所不同，本文我們主要比較的是方便觀察的俯仰力矩，可以發(fā)現(xiàn)多孔對俯仰力矩的影響最為明顯，而半圓形擋流板和圓形擋流板對俯仰力矩的抑制作用相差無幾。仿真實驗是在準(zhǔn)靜態(tài)理想化的直線制動良好路面上進行制動操作的策略下進行的，在填充率較低時，擋流板的作用遠遠不如填充率較高時來的明顯，而且其中多孔擋流板的作用尤其突出。

本文也對液體速度和加速度進行了仿真，勻速行駛速度較低，設(shè)置仿真時間為勻速行駛1s后，進行制動3.5s的過程。我們進行了在全制動的的操作下，研究最后一塊防波板周圍流體的速度以及加速度。制動4s后，液體表面已經(jīng)停止晃動，通過展示的三種不同擋流板排列的模擬過程中自由流體速度的的演變，可以看出，擋流板的存在顯著地抑制了流體運動，但由于我們設(shè)置的是低速，所以各個擋流板對液體晃動的抑制的差異不是很明顯。

隨后我們又對各個不同防波板對縱向力、垂直力還有俯仰力矩進行了綜合對比分析發(fā)現(xiàn)對流體的俯仰力矩還有縱向力影響最明顯的依然是多孔的擋流板。因為是直線行駛的全制動，所以各擋流板對流體的垂直力影響差異不是很明顯。

在主要對縱向全制動工況仿真分析了之后，本文為了研究防波板對橫向動力學(xué)的影響，對液罐車在規(guī)避變道操作時的制動轉(zhuǎn)向工況進行了仿真研究。為了方便觀察研究，我們設(shè)置填充率為50%，初始車輛速度為70km/h，規(guī)避變道操作完成時為48km/h。通過在準(zhǔn)靜態(tài)模型中對CFD液罐車的橫向加速度響應(yīng)進行比較，若僅觀察橫向的液體加速度響應(yīng)，防波板對制動轉(zhuǎn)向工況中橫向的流體晃動影響很有限，影響程度小、時間短。所以我們可以知道防波板的主要抑制縱向的液體晃動。

4 ?結(jié)論

通過合理簡化并進行構(gòu)造等效模型的方法，先進行準(zhǔn)靜態(tài)分析，在理想狀態(tài)下對罐體就行全制動的液體激勵仿真，通過在ANSYS中標(biāo)記擋流板周邊某一點作為等效力學(xué)的研究對象，對比其的速度、加速度以及罐內(nèi)流體的俯仰力矩，來進行分析擋流板對液體晃動的影響，從而得出擋流板對罐車穩(wěn)定性的重要作用。仿真結(jié)果很理想，本文通過三種不同的防波板仿真對比，從而得出所設(shè)計的多孔擋流板對流體晃動可以有效地降低過渡峰和平均其俯仰力矩，而且其質(zhì)量也是最輕，結(jié)構(gòu)得到優(yōu)化。在規(guī)避變道工況中知道了防波板對橫向加速度的影響極小。對增大側(cè)翻閥值從而減少側(cè)翻事故提供了理論依據(jù)。為設(shè)計的防波板模型的合理性以及新型防波板在?；芬汗捃囍型茝V論證了可行性。本文對于研究油罐車的側(cè)傾及其穩(wěn)定性有著一定的指導(dǎo)意義。

參考文獻

[1]王次安，王宏大，張蕾.基于流固耦合的罐裝車結(jié)構(gòu)優(yōu)化設(shè)計[J].專用汽車，2014（10）

[2]田小平，白莉.液罐車動態(tài)結(jié)構(gòu)強度有限元分析[J].現(xiàn)代制造技術(shù)與裝備，2017（08）