罐車制動工況液體晃動的流固耦合分析.doc

罐車制動工況液體晃動的流固耦合分析 張韜柯龍燕白光建 中航（上海）汽車技術(shù)有限公司上海200092 摘要：采用VOF模型對罐車制動過程中內(nèi)部液體的自由晃動過程進(jìn)行了模擬，計(jì)算得出罐體各個(gè)部分的應(yīng)力分布。計(jì)算結(jié)果表明，制動時(shí)前封頭和前四個(gè)防波板受力逐漸上升，在達(dá)到峰值后迅速下降；后封頭和第五個(gè)防渡板受力則呈逐漸減小的趨勢。隨著與后封頭距離的增大，防波板承受最大應(yīng)力也相應(yīng)增大，單個(gè)防波板承受最大應(yīng)力位置處于底部孔邊緣處。 關(guān)鍵詞：罐車兩相流液體晃動流固耦合 ：U469.6+1.02：A：1004-0226(xx)08-0088-04 l前言 罐車在行駛過程中由于加速、減速或轉(zhuǎn)彎等運(yùn)動狀態(tài)的改變，會產(chǎn)生慣性力，進(jìn)而引起液體的晃動，出現(xiàn)罐車罐體與介質(zhì)的流固耦臺效應(yīng)。在這種流固耦臺系統(tǒng)的作用下，一方面可能會導(dǎo)致罐車的動力不穩(wěn)定，另一方面罐體中介質(zhì)的晃動會對罐體產(chǎn)生顯著的壓力，嚴(yán)重時(shí)會造成罐體結(jié)構(gòu)破壞。罐體內(nèi)一般都布置有防波板或阻浪板，將罐體內(nèi)的液體分成多個(gè)小艙室，在流固耦合作用下，整個(gè)罐體內(nèi)液體的沖擊力將由各防波板分擔(dān)，因此需對防波板的受力過程進(jìn)行分析，以便通過在罐體內(nèi)部設(shè)置合理的防晃裝置降低罐車防波板損壞的可能性。 本文使用流固耦合分析軟件對罐車在緊急制動工況下罐內(nèi)的油料晃動過程進(jìn)行了分析計(jì)算，研究了罐體各部分受力隨時(shí)間的變化趨勢以及防波板的應(yīng)力分布情況。 2計(jì)算模型和計(jì)算方法 2.1控制方程組 本文對于流體部分采用VOF多相流模型計(jì)算罐內(nèi)氣液兩相非定常湍流流動，控制方程組包括連續(xù)性方程和動量守恒方程。 2.2罐體模型和網(wǎng)格劃分 以某型號油罐車為研究對象，其罐體模型如圖1所示。罐體由前后封頭、5個(gè)防波板和罐體組成。其中防波板將罐體分割為六個(gè)區(qū)域，區(qū)域之間通過兩個(gè)開孔相連。罐體總長約10.2m，寬約23m，高約1.6m。 計(jì)算時(shí)，坐標(biāo)系固定于罐體上，原點(diǎn)位于罐體前封頭端面底部中心位置，磚由以罐體前進(jìn)方向?yàn)檎较蚯移叫杏诘孛?，y軸垂直于地面，向上為正方向，z軸平行于地面，向右為正方向。模型網(wǎng)格劃分采用混合網(wǎng)格結(jié)構(gòu)，使用移動坐標(biāo)平面將罐體切分為6個(gè)小罐體和5個(gè)包含防波板的部分，小罐體可采用六面體掃掠網(wǎng)格，其余部分采用四面體網(wǎng)格，網(wǎng)格過渡部分保證平滑且節(jié)點(diǎn)一一對應(yīng)，如圖2所示。這樣可以大大減少網(wǎng)格總數(shù)，提高計(jì)算效率。 2.3邊界條件 罐體及防波板表面設(shè)置為無滑移壁面邊界條件。多相流采用VOF模型，可以精確模擬氣液分界面的變化情況以模擬液體晃動過程。初始?xì)庖航唤缑嫫叫杏趚z平面，液體和氣體相對罐體速度為0，氣體壓強(qiáng)為1.013xlO5Pa。湍流模型采用標(biāo)準(zhǔn)k-a模型和標(biāo)準(zhǔn)壁面函數(shù)模型。 罐車以80km/h的初始速度正向行駛，在O.ls時(shí)施加0.5g的反向加速度模擬制動工況。罐車內(nèi)介質(zhì)為油，密度900kg/m3，粘性0.048kgms，充裝率80%。實(shí)體材料均為Q235鋼，彈性模量2.06xlOlPa，泊松比為0.250制動速度曲線如圖3所示。 2.4求解器參數(shù)設(shè)置 速度和壓強(qiáng)的耦合處理采用PISO方法，壓力修正方程離散格式采用BodyForceWeighted格式，對流項(xiàng)離散格式采用二階迎風(fēng)格式。非定常計(jì)算中時(shí)間步長為0.002s。 3制動工況計(jì)算結(jié)果及分析 圖4、5為前后封頭和各防波板壓力隨時(shí)間的變化曲線。從各曲線圖可以看出，罐車制動后在0.52s時(shí)液體各部分壓力達(dá)到最大值，而后開始逐漸減小。 在制動之前，罐體各部分由于液體重力的存在而受到的壓力基本相同。在制動初期，前封頭和前4個(gè)防波板所受壓力呈現(xiàn)震蕩式上升。這是由于初期液體晃動趨勢較小，在經(jīng)過四次震蕩后逐漸向峰值接近。后封頭和第五個(gè)防波板受力曲線相似，均呈現(xiàn)逐漸減小的趨勢，同時(shí)越靠近后封頭的防波板受力越小，這是由于液體正向運(yùn)動的慣性導(dǎo)致罐體前后部的受力趨勢不同。 罐體受力為峰值時(shí)，液體晃動最為劇烈，隨后受力減小。這是由于液體受到罐體一側(cè)壁面的阻擋開始做折返運(yùn)動，涌向另一側(cè)。在整個(gè)晃動過程中前封頭受到的峰值壓力最大。 圖6、7為初始時(shí)刻、0.52s即壓力達(dá)到峰值時(shí)罐體內(nèi)液相分布圖。圖8-11為上述兩個(gè)時(shí)刻罐體和防波板壓力分布圖。壓力分布情況與前述分析結(jié)果一致。罐內(nèi)液體受晃動影響集中在每個(gè)腔內(nèi)的前部，因此每個(gè)腔內(nèi)前部壓力大于后部壓力。0.52s時(shí)液體處在向前運(yùn)動的動量最大處，因此對罐體的沖擊最大，而后由于罐體的阻擋向反方向運(yùn)動。 圖1214為0.52s時(shí)，罐體和防波板所受應(yīng)力云圖和位移云圖。從圖中可以看出初始時(shí)刻罐體和防波板受到的力較小，0.52s時(shí)，壓力達(dá)到最大值后，應(yīng)力和位移也達(dá)到最大值，防波板最大應(yīng)力約545MPa，受力點(diǎn)在前部第一、二個(gè)防波板的固定點(diǎn)折彎角位置。 圖15、16為模型封頭、防波板結(jié)構(gòu)修改后的示意圖和應(yīng)力計(jì)算云圖。主要將底部圓孔附近面修改為規(guī)則的平整截面、增加折邊寬度并加入三角支撐結(jié)構(gòu)。從計(jì)算結(jié)果中可知模型最大應(yīng)力減小為17419MPa，小于材料屈服極限要求。 4結(jié)語 使用有限體積法和FLUENT軟件，對充有空氣與油的臥式柱形罐車內(nèi)部介質(zhì)晃動進(jìn)行了模擬，確定了VOF方法用于自由液面的定義，標(biāo)準(zhǔn)ke湍流模型作為湍流計(jì)算模型，并確定了流固耦合計(jì)算模型，為進(jìn)一步研究探索了方法和方向。 罐車制動后在0.52s時(shí)液體各部分壓力達(dá)到最大值。前封頭和前四個(gè)防波板所受壓力趨勢相同，第五個(gè)防波板和后封頭受壓趨勢相同。前封頭在整個(gè)晃動過程中所受峰值壓力最大，而后各部分峰值壓力值逐漸減小。 防波板所受最大應(yīng)力約為545MPa，受力點(diǎn)在前部第一、二個(gè)防波板的固定點(diǎn)折彎角的位置。因此可考慮在此處適當(dāng)改變結(jié)構(gòu)以減小應(yīng)力集中，以提高罐車的可靠性和安全性。經(jīng)過結(jié)構(gòu)改進(jìn)后應(yīng)力降至174.19