?；芬汗捃嚪啦ò鍖?duì)其行駛穩(wěn)定性的影響

王林　朱天軍　胡偉　栗圣博

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

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

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

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

1.2 ?罐體幾何模型的建立

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

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

1.3 ?等效晃動(dòng)幾何模型的建立

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

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

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

（1）

;（2）

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

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

（3）

（4）

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

（5）

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

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

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

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

（9）

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

（10）

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

2.1 ?三維模型的建立

因?yàn)楣摅w和擋流板的三維模型較復(fù)雜，所以本文所需要的三維模型是在CATIA軟件進(jìn)行建模的。因?yàn)楸疚难芯康氖遣煌瑩趿靼鍖?duì)油罐車穩(wěn)定性的影響，所以罐體簡(jiǎn)化模型的尺寸和截面都是統(tǒng)一的。如圖所示：

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

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

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

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

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

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

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

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

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

4 ?結(jié)論

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

參考文獻(xiàn)

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

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