陳雷，張克正，封星

(1.煙臺(tái)職業(yè)學(xué)院 汽車與船舶工程系，山東 煙臺(tái) 264670；2.魯東大學(xué) 交通學(xué)院，山東 煙臺(tái) 264025； 3.大連海事大學(xué) 輪機(jī)工程學(xué)院，遼寧 大連 116026)

0 引言

液化天然氣(liquefied natural gas,LNG)動(dòng)力船舶營(yíng)運(yùn)時(shí)，非滿載燃料儲(chǔ)罐內(nèi)部自由液面會(huì)受到船舶營(yíng)運(yùn)狀態(tài)及外界作用力的影響而呈現(xiàn)非線性和不規(guī)則的復(fù)雜運(yùn)動(dòng)形態(tài)[1]，對(duì)儲(chǔ)罐壁面產(chǎn)生巨大的液體沖擊力。液體沖擊可能造成罐體破裂、液體泄漏等事故，對(duì)船員與船舶的營(yíng)運(yùn)安全帶來(lái)嚴(yán)重危害[2-5]。目前，國(guó)內(nèi)外對(duì)燃料儲(chǔ)罐內(nèi)部沖擊的研究多集中于比對(duì)數(shù)值模擬與試驗(yàn)工況下的液面振型，以此來(lái)證實(shí)數(shù)值模擬的正確性與適用性[6-9]，或者是觀察液體沖擊時(shí)，黏度、密度、加速度等參數(shù)變化對(duì)罐體內(nèi)部液體的動(dòng)壓力和自由液面波動(dòng)情況的影響，但很少分析充注率與容積變化對(duì)罐內(nèi)沖擊力的影響，或直觀揭示液體最大沖擊力與各影響因素之間的具體關(guān)系[10-15]。

本文以30 m3臥式筒形LNG燃料儲(chǔ)罐為研究對(duì)象，模擬罐內(nèi)液體縱蕩時(shí)充注率和容積對(duì)罐內(nèi)最大沖擊力的影響，并將其與慣性力進(jìn)行比對(duì)，分析最大沖擊力與慣性力之比隨參數(shù)變化的規(guī)律，進(jìn)而得到最大液體沖擊力與影響因素的關(guān)系。

1 建模

圖1 儲(chǔ)罐模型

液體縱蕩沖擊的根本原因是儲(chǔ)罐內(nèi)部自由液面導(dǎo)致的晃動(dòng)?；蝿?dòng)與儲(chǔ)罐內(nèi)部液體的充注率、容積以及液體的黏度、密度、加速度等參數(shù)存在很大關(guān)系。研究對(duì)象為容積為30 m3的C型LNG儲(chǔ)罐，長(zhǎng)8270 mm，內(nèi)徑DN=2200 mm，筒體長(zhǎng)(不包含封頭直邊)7090 mm，采用標(biāo)準(zhǔn)橢圓封頭。設(shè)定坐標(biāo)原點(diǎn)于罐左側(cè)(后端)封頭頂點(diǎn)位置，x軸正方向?yàn)榇棒疾糠较?，y軸負(fù)方向?yàn)橹亓Ψ较颍瑉軸方向?yàn)榇白笥蚁戏较?，如圖1所示。假定罐體是剛體，壁面無(wú)滑移，不考慮氣-液分界面的表面張力。

在Fluent求解器中采用瞬態(tài)模型，湍流模型采用標(biāo)準(zhǔn)k-ε模型，多相流模型采用流體體積函數(shù)，即目標(biāo)流體的體積與網(wǎng)格體積的比值(volume of fluid，VOF)模式，采用速度與壓力耦合PISO算法[16]；對(duì)流項(xiàng)離散格式采用一階迎風(fēng)，壓力的空間離散采用體積力加權(quán)(Body-Force-Weighted)方式，采用UDF程序加載定義動(dòng)量源項(xiàng)，對(duì)流體施以加速度。

根據(jù)求解器設(shè)置與流體動(dòng)力學(xué)要求，得到罐內(nèi)液體縱蕩應(yīng)遵循的基本方程[17]。

1)連續(xù)性方程

式中：ρ為流體密度，t為流動(dòng)時(shí)間，u、v、w為速度矢量在x、y、z方向上的分量。

2)動(dòng)量平衡方程

式中：t為時(shí)間，s；ρf為單位體積力作用于流體質(zhì)量上的非接觸力；p為流體流場(chǎng)的壓力；τ為黏性應(yīng)力分量，σ為應(yīng)力張量，可以表達(dá)為：

3)湍流動(dòng)能方程k與擴(kuò)散方程ε表達(dá)式分別為：

2 模型驗(yàn)證

2.1 流體網(wǎng)格驗(yàn)證

圖2 儲(chǔ)罐初始內(nèi)部靜壓強(qiáng)分布

高質(zhì)量的網(wǎng)格是Fluent計(jì)算湍流運(yùn)動(dòng)的必要條件[18]。流體空間采用六面體網(wǎng)格形式，共60 988個(gè)，為驗(yàn)證劃分網(wǎng)格的精確度，暫不對(duì)流體施加速度邊界條件。介質(zhì)的類型、黏度對(duì)介質(zhì)晃動(dòng)的影響極小，在研究液體晃動(dòng)時(shí)，國(guó)內(nèi)外學(xué)者多采用理想無(wú)黏性介質(zhì)。因此，儲(chǔ)罐內(nèi)部以水和空氣為流體介質(zhì)，替代LNG，定義水和空氣的體積比為1：1，初始化后觀察儲(chǔ)罐內(nèi)部靜壓強(qiáng)的分布(如圖2所示，淺灰色為液相，深灰色為氣相，圖中單位為kPa)，并與計(jì)算理論值進(jìn)行比對(duì)[19]。

由圖2可見：罐內(nèi)壓強(qiáng)分布層次清晰，自罐頂?shù)焦薜壮示€性分布，罐底靜壓強(qiáng)最大，為10.4 kPa。

最大理論靜壓強(qiáng)

Pmax=ρgh，

式中：ρ為液體密度,g為重力加速度,h為液體深度。

根據(jù)儲(chǔ)罐尺寸與充裝率可計(jì)算儲(chǔ)罐內(nèi)部Pmax=10.8 kPa。

計(jì)算靜壓強(qiáng)與理論靜壓強(qiáng)非常相近，故流體網(wǎng)格符合要求。

2.2 參數(shù)設(shè)置驗(yàn)證

圖3 縱向晃動(dòng)頻率的計(jì)算結(jié)果與測(cè)量結(jié)果對(duì)比

為驗(yàn)證模型參數(shù)設(shè)置的正確性，參照文獻(xiàn)[6]中的參數(shù)建立罐體模型，計(jì)算不同充注率下內(nèi)部介質(zhì)的縱向晃動(dòng)頻率，并與測(cè)量結(jié)果進(jìn)行對(duì)比，如圖3所示。

由圖3可知:在充注率為0.1～0.9時(shí)，內(nèi)部介質(zhì)縱向晃動(dòng)頻率的計(jì)算結(jié)果與試驗(yàn)測(cè)量結(jié)果較為吻合。只在充注率為0.2時(shí)兩者偏差最大，偏差為4.91%。因此，計(jì)算結(jié)果具有較高的可信度，參數(shù)設(shè)置可行。

2.3 案例計(jì)算

圖4 儲(chǔ)罐前后封頭縱向液體沖擊力時(shí)域響應(yīng)

文獻(xiàn)[10]證明以水作為液體介質(zhì)進(jìn)行研究有一定的通用性?，F(xiàn)以水與空氣作為介質(zhì)，模擬LNG動(dòng)力船舶在制動(dòng)過(guò)程中燃料儲(chǔ)罐內(nèi)部液體縱向晃動(dòng)過(guò)程，觀察儲(chǔ)罐內(nèi)部液體沖擊力的變化。文獻(xiàn)[20]規(guī)定，沿船舶運(yùn)動(dòng)方向向前和向后碰撞的最大設(shè)計(jì)加速度為2g。假設(shè)在時(shí)間t=0時(shí)，對(duì)充注率為0.50的罐內(nèi)流體沿運(yùn)動(dòng)方向施以2g的加速度，使內(nèi)部液體與儲(chǔ)罐產(chǎn)生相對(duì)速度，內(nèi)部產(chǎn)生晃動(dòng)；在t=2 s時(shí)將加速度減為0，觀察0～10 s內(nèi)罐中液體與空氣的晃動(dòng)情況，監(jiān)視整個(gè)晃動(dòng)過(guò)程中前后封頭所受液體沖擊力F隨時(shí)間的變化規(guī)律，結(jié)果如圖4所示。

圖5 t=0.60 s時(shí)儲(chǔ)罐內(nèi)部氣-液兩相體積分布

由圖4可以看出：當(dāng)船舶采取緊急制動(dòng)或碰撞后，儲(chǔ)罐內(nèi)部液體向前移動(dòng)沖擊前側(cè)封頭，在t=0.60 s時(shí)，前封頭液體沖擊力出現(xiàn)峰值，F(xiàn)max=932.179 kN，此刻氣-液兩相分布如圖5所示(圖中為氣體體積分?jǐn)?shù)，淺灰色為液相，深灰色為氣相)。在t=0.945 s時(shí)，后封頭液體沖擊力達(dá)到最大值，F(xiàn)max=307.446 kN，僅為前側(cè)封頭受力的1/3。在t=2 s，即加速度減為0的瞬間，前后封頭受力突然降低，罐內(nèi)介質(zhì)進(jìn)入自由晃動(dòng)階段，沖擊力維持在較小范圍內(nèi)變化，并逐漸趨向于0。這與船舶制動(dòng)時(shí)罐內(nèi)介質(zhì)的運(yùn)動(dòng)狀態(tài)相符，因此，模型是可信的。

3 液體沖擊影響因素分析

本文主要研究液體沖擊力對(duì)儲(chǔ)罐及船舶營(yíng)運(yùn)安全的影響，故只監(jiān)測(cè)前后封頭最大液體沖擊力Fmax，并將其與介質(zhì)慣性力F0進(jìn)行比較，找出液體沖擊力的變化規(guī)律。定義介質(zhì)慣性力

F0=ma=ρVa，

(1)

式中：a為加速度，m為液體慣性質(zhì)量，V為液體體積。

定義放大系數(shù)

(2)

定義絕對(duì)偏差率

(3)

此前，學(xué)者們針對(duì)密度、黏度、加速度等參數(shù)與儲(chǔ)罐內(nèi)部液體沖擊力的關(guān)系進(jìn)行了大量試驗(yàn)與研究，結(jié)果表明：在其他參數(shù)保持不變的情況下，最大液體沖擊力與黏度無(wú)關(guān)，與密度及加速度呈線性正比關(guān)系，液體沖擊放大系數(shù)不隨密度、黏度、加速度的變化而變化。本文主要分析儲(chǔ)罐容積與充注率對(duì)罐內(nèi)液體沖擊的影響，探討放大系數(shù)的變化規(guī)律。

3.1 儲(chǔ)罐容積

圖6 不同儲(chǔ)罐容積縱向沖擊力的時(shí)域響應(yīng)

不同型式燃料儲(chǔ)罐的區(qū)別主要是罐體長(zhǎng)度，而儲(chǔ)罐直徑相差不大。保持其他參數(shù)不變，計(jì)算不同儲(chǔ)罐容積對(duì)縱蕩過(guò)程的液體沖擊力。針對(duì)內(nèi)徑DN=2200 mm的標(biāo)準(zhǔn)橢圓封頭儲(chǔ)罐，選取長(zhǎng)度L分別為4700、8270、10 100 mm的儲(chǔ)罐作為研究對(duì)象，對(duì)其晃蕩過(guò)程進(jìn)行數(shù)值模擬，觀察不同容積儲(chǔ)罐內(nèi)部瞬時(shí)液體沖擊力的時(shí)域響應(yīng)，如圖6所示。

由圖6可以看出：儲(chǔ)罐容積不同時(shí)，瞬時(shí)最大液體沖擊力的出現(xiàn)時(shí)刻隨儲(chǔ)罐容積的增大而滯后；當(dāng)儲(chǔ)罐容積增大時(shí)，儲(chǔ)罐內(nèi)部瞬時(shí)最大液體沖擊力隨之增大。

不同工況下瞬時(shí)Fmax出現(xiàn)時(shí)間、F0及K等數(shù)據(jù)如表1所示。

表1 不同儲(chǔ)罐容積的液體沖擊參數(shù)

圖7 不同充注率下縱向沖擊力的時(shí)域響應(yīng)

從表1可以看出：若液體密度、黏度、加速度等參數(shù)保持不變，最大液體沖擊力與容積成線性正比關(guān)系。在容積發(fā)生變化時(shí)，K平均為3.16，最大偏差為2.3%，隨容積變化不大。

3.2 充注率

文獻(xiàn)[21]規(guī)定，C型LNG燃料儲(chǔ)罐的充裝極限不大于90%，因此，充注率λ≤0.90?，F(xiàn)分別模擬充注率λ為0.05、0.25、0.50、0.63、0.75、0.862、0.90等7種工況下的內(nèi)部介質(zhì)晃動(dòng)情況，觀察不同充注率下沖擊力的時(shí)域響應(yīng)，如圖7所示。

從圖7可以看出，隨著充注率的增大，最大液體沖擊力出現(xiàn)的時(shí)刻不斷提前，這是因?yàn)槌渥⒙试龃蠛螅后w與前封頭的接觸面積增大，向其運(yùn)動(dòng)進(jìn)行沖擊所需時(shí)間減少。λ接近極限充注率0.90時(shí)，最大沖擊力沒(méi)有出現(xiàn)在前封頭位置，而是出現(xiàn)在后封頭位置，相應(yīng)時(shí)間也有所延遲。不同充注率下液體Fmax出現(xiàn)時(shí)間與k等參數(shù)數(shù)據(jù)，如表2所示。

表2 不同充注率時(shí)液體的沖擊參數(shù)

圖8 Fmax、K隨λ的變化曲線

Fmax、K隨λ的變化曲線如圖8所示。

由圖8可知，K與λ呈多項(xiàng)式函數(shù)關(guān)系，隨λ的增大而減小。儲(chǔ)罐受到的Fmax呈現(xiàn)先增大后減小的趨勢(shì)：在λ較小時(shí)，液體沖擊力隨λ的增大而增大，這主要是受到液體慣性力增加的影響；在λ達(dá)到0.63左右時(shí)，F(xiàn)max達(dá)到極值后呈下降趨勢(shì)，這是因?yàn)棣溯^大時(shí)，罐內(nèi)液體不像低充裝率時(shí)在罐內(nèi)翻騰晃動(dòng)，此時(shí)，介質(zhì)涌向前封頭的時(shí)間明顯縮短，液體全部涌向前封頭，氣體被擠壓至后封頭位置，儲(chǔ)罐內(nèi)部的液體沖擊力很快達(dá)到穩(wěn)定值。

4 結(jié)果分析

4.1 擬合關(guān)系式

通過(guò)液體沖擊影響因素分析可以看出，液體沖擊放大系數(shù)與容積、黏度、密度和加速度無(wú)關(guān)，與充注率呈一定函數(shù)關(guān)系。根據(jù)放大系數(shù)隨充注率變化曲線的特點(diǎn)，擬合充注率不同時(shí)放大系數(shù)的變化曲線，采用多項(xiàng)式表示K-λ曲線為：

K=f(λ)=a1λn+a2λn-1+…+anλ+an+1，

(4)

式中：n為多項(xiàng)式的最高次冪；ai為冪次降序的多項(xiàng)式系數(shù)，i=1，2，…，n，n+1。

定義a=(a1a2…anan+1)為多項(xiàng)式系數(shù)，λ=(λn,λn-1,…λ,1)T為高度冪矢量，則可用矩陣表達(dá)為：

K=aλ。

采用Matlab軟件的指令格式f(λ)=polyfit(λ,N,n)，對(duì)已知的λ和K的離散數(shù)據(jù)及最高冪次n(n取2～5)進(jìn)行最小二乘多項(xiàng)式曲線擬合，得到系數(shù)a，得出2～5次多項(xiàng)式的擬合曲線特性方程[22-23]。表3中列出了不同最高次冪的最大偏差，從最大偏差可以看到：4次多項(xiàng)式的擬合偏差為2.74%，更高冪次的

表3 不同冪次擬合曲線偏差

擬合精度也沒(méi)有得到顯著提高，而且會(huì)使擬合曲線的光滑性變差，多項(xiàng)式出現(xiàn)擺振特性，反映更多的數(shù)據(jù)量化誤差[24-25]。

由此確定：

K=23.714 4λ4-46.584 9λ3+28.133λ2-
8.365 2λ+4.684 3。

(5)

結(jié)合式(1)～(5)，可得Fmax與F0的關(guān)系式為：

Fmax=(23.714 4λ4-46.584 9λ3+28.133λ2-8.365 2λ+4.684 3)F0。

(6)

4.2 結(jié)果驗(yàn)證

現(xiàn)數(shù)值模擬5種不同介質(zhì)在容積、加速度、充注率等不同工況下的結(jié)果，驗(yàn)證式(6)的正確性與適用范圍，對(duì)比結(jié)果如表4所示。

表4 Fmax-F0關(guān)系式驗(yàn)證比對(duì)

圖9 不同形狀燃料儲(chǔ)罐縱向沖擊力時(shí)域響應(yīng)

由表4可以得到：定容積燃料儲(chǔ)罐沿運(yùn)動(dòng)方向制動(dòng)時(shí)，內(nèi)部介質(zhì)最大沖擊力的計(jì)算公式數(shù)值模擬的計(jì)算誤差均在2%以內(nèi)，符合工程精度要求；關(guān)系式(6)在不同容積、不同介質(zhì)、不同加速度、不同充裝率等條件下都是有效的。為了驗(yàn)證式(6)的適用性，針對(duì)不同形狀的燃料儲(chǔ)罐進(jìn)行數(shù)值模擬，選擇30 m3方形儲(chǔ)罐與平封頭儲(chǔ)罐進(jìn)行驗(yàn)證，給定液體的充注率λ=0.50，加速度為2g，模擬計(jì)算結(jié)果如圖9所示。

從圖9中可以看出，儲(chǔ)罐的形狀對(duì)儲(chǔ)罐內(nèi)部最大液體沖擊力影響較大，不同形狀的儲(chǔ)罐最大沖擊力不同，與公式計(jì)算結(jié)果Fmax=938.30 kN偏差較大，故上述沖擊力計(jì)算公式并不適用于不同形狀的燃料儲(chǔ)罐。

5 結(jié)論

1)針對(duì)標(biāo)準(zhǔn)橢圓頭筒型儲(chǔ)罐的內(nèi)部液體沖擊，采用單一變量法，分析容積、充注率對(duì)最大瞬時(shí)液體沖擊力的影響。結(jié)果表明：最大瞬時(shí)液體沖擊力的出現(xiàn)時(shí)刻隨儲(chǔ)罐容積的增大而滯后，隨充注率的增大而提前；最大瞬時(shí)液體沖擊力隨容積的增大而增大，隨充注率的增大呈現(xiàn)先增大后減小的趨勢(shì)。

2)通過(guò)分析仿真結(jié)果得出最大沖擊力與介質(zhì)慣性力的關(guān)系式，表明：放大系數(shù)與充注率呈4次多項(xiàng)式關(guān)系，與容積、黏度、密度和加速度等參數(shù)無(wú)關(guān)。

3)分析擬合的Fmax-F0關(guān)系式的適用性。結(jié)果表明：針對(duì)標(biāo)準(zhǔn)橢圓頭筒型儲(chǔ)罐，關(guān)系式在不同容積、黏度、密度和加速度下都是有效的，計(jì)算誤差在2%以內(nèi)，符合工程精度要求，但不適用于其它形狀的儲(chǔ)罐。