基于CFD的內(nèi)傾船迎浪中縱向非線性運(yùn)動(dòng)研究

劉 樂(lè) 毛筱菲

(武漢理工大學(xué)交通學(xué)院 武漢 430063)

基于CFD的內(nèi)傾船迎浪中縱向非線性運(yùn)動(dòng)研究

劉 樂(lè) 毛筱菲

(武漢理工大學(xué)交通學(xué)院 武漢 430063)

相比傳統(tǒng)船型而言，有航速內(nèi)傾船在大波陡波浪中，縱向運(yùn)動(dòng)可能會(huì)產(chǎn)生甲板上浪和非線性水動(dòng)力問(wèn)題.以DTMB 5613船模為研究對(duì)象，運(yùn)用CFD軟件OpenFOAM開(kāi)展了Fr=0,0.2及0.367，波陡H/λ=0.01,0.025及0.05的迎浪規(guī)則波下內(nèi)傾船縱搖和垂蕩運(yùn)動(dòng)的數(shù)值模擬.結(jié)果表明,高航速內(nèi)傾船在大波陡波浪中的縱向運(yùn)動(dòng)具有強(qiáng)烈的非線性，船體的內(nèi)傾設(shè)計(jì)引起的甲板上浪等非線性現(xiàn)象是導(dǎo)致船體縱向運(yùn)動(dòng)非線性的主要因素.

內(nèi)傾船；縱搖；垂蕩；波陡；非線性

0 引 言

與常規(guī)船型不同，內(nèi)傾船水線以下為外凸的曲線形式，水線以上部分采取舷側(cè)內(nèi)傾10°～20°，上層建筑側(cè)壁內(nèi)傾7°～15°設(shè)計(jì)，極大的降低了船體的雷達(dá)反射面積，表現(xiàn)出出色的隱身性能，其穿浪首設(shè)計(jì)可以減小航行中首部興波阻力，降低高海況中縱搖幅值的同時(shí)提高航行速度，但內(nèi)傾船除了具有常規(guī)船型的耐波性問(wèn)題之外，還有一些獨(dú)有的運(yùn)動(dòng)特性：首柱和舷側(cè)內(nèi)傾設(shè)計(jì)導(dǎo)致船體垂蕩運(yùn)動(dòng)劇烈，容易產(chǎn)生甲板上浪現(xiàn)象，使縱搖和垂蕩運(yùn)動(dòng)具有強(qiáng)烈的非線性，在迎浪中，這種非線性運(yùn)動(dòng)可能會(huì)導(dǎo)致內(nèi)傾船產(chǎn)生參數(shù)橫搖現(xiàn)象，威脅船舶安全；而隨浪中，船舷內(nèi)傾的設(shè)計(jì)使船舶在發(fā)生橫搖后，回復(fù)力矩不足，可能導(dǎo)致純穩(wěn)性喪失和橫甩，同樣不利于船舶安全；橫浪中，內(nèi)傾的舷側(cè)以及浮心低于重心的設(shè)計(jì)使內(nèi)傾船在發(fā)生橫傾后更難回復(fù)到平衡位置，高聳的上層建筑導(dǎo)致風(fēng)載荷較大，船舶橫向穩(wěn)定性較差.

文中運(yùn)用粘性CFD軟件OpenFOAM計(jì)算了內(nèi)傾船在迎浪規(guī)則波中的縱搖和垂蕩頻響曲線，與基于勢(shì)流理論的AQWA計(jì)算結(jié)果進(jìn)行了比較，然后探究了不同航速的內(nèi)傾船在不同波陡的迎浪規(guī)則波下的縱搖和垂蕩加速度時(shí)歷曲線，預(yù)報(bào)了內(nèi)傾船在縱浪中的非線性運(yùn)動(dòng)特性，為后續(xù)內(nèi)傾船型研究奠定基礎(chǔ).

1 數(shù)值模擬

1.1 光體計(jì)算模型

DTMB 5613是NSWCCD (Naval Surface Warfare Center, Carderock Division)提出的研究模型，其主要參數(shù)見(jiàn)表1，二維橫剖線圖及三維光體模型見(jiàn)圖1～2.

表1 DTMB 5613模型主要參數(shù)

圖1 DTMB 5613橫剖線圖

圖2 DTMB 5613光體模型

1.2 計(jì)算原理

OpenFOAM以RANS方程作為計(jì)算內(nèi)傾船在波浪中運(yùn)動(dòng)的基本方程，時(shí)均連續(xù)性方程具體形式為

(1)

動(dòng)量方程具體形式為

(2)

式中：ρ為流體密度；μ為流體粘度；p為靜水壓力；Si為質(zhì)量力；ui，uj為速度分量.

湍流模型選擇SSTk-ω模型，船體運(yùn)動(dòng)使用動(dòng)網(wǎng)格技術(shù)，自由液面使用VOF法實(shí)時(shí)捕獲.

1.3 計(jì)算域設(shè)置

3.《秀才胡同》與《東風(fēng)破》的區(qū)別，猶如曹子建之于陶淵明的古文、辛稼軒之于姜白石詞句的區(qū)別，前者辭藻華美、引經(jīng)據(jù)典，后者相較質(zhì)樸平易、清白寫(xiě)意。前者得稼軒之長(zhǎng)，善用典故，辭藻華麗動(dòng)人；后者雖用典故，但旨在白描，清白明了，得白石道人之精髓。

由于船模關(guān)于中線面對(duì)稱，為了減小計(jì)算量，取半模進(jìn)行計(jì)算，計(jì)算域大小設(shè)置為長(zhǎng)×寬×高=8.0L×2.0L×3.0L,邊界設(shè)置：①入口距離船體首垂線1.0L，邊界條件設(shè)置為速度入口；②出口距離船體尾垂線6.0L，邊界條件設(shè)置為壓力出口；③流域的上邊界設(shè)置成速度入口；④流域的左、右邊界，以及底部均設(shè)置為對(duì)稱面；⑤船體表面設(shè)置成無(wú)滑移、不可穿透壁面.

1.4 網(wǎng)格劃分

OpenFOAM采用blockMesh工具進(jìn)行背景網(wǎng)格劃分，船體附近的貼壁網(wǎng)格利用snappyHexMesh工具自動(dòng)劃分，同時(shí)為了準(zhǔn)確捕捉自由面變化，對(duì)自由面附近的網(wǎng)格進(jìn)行細(xì)化加密，波高使用30個(gè)網(wǎng)格捕捉，波長(zhǎng)使用100個(gè)網(wǎng)格捕捉，計(jì)算迭代次數(shù)5次，同時(shí)要局部細(xì)化船模貼壁網(wǎng)格尺寸，保證船體壁面y+值在30～300之間，半模網(wǎng)格數(shù)最大為170萬(wàn)左右.

2 結(jié)果及分析

2.1 線性頻響計(jì)算

文中分別計(jì)算了內(nèi)傾船在無(wú)航速Fr=0和高航速Fr=0.367時(shí)的縱搖運(yùn)動(dòng)和垂蕩運(yùn)動(dòng)傳遞函數(shù)曲線.圖3將OpenFOAM計(jì)算的無(wú)航速內(nèi)傾船的頻響曲線與勢(shì)流軟件AWQA計(jì)算結(jié)果進(jìn)行比較，由圖3可知，內(nèi)傾船在零航速時(shí)的峰值點(diǎn)在角頻率為3.917 rad/s處，即λ/L=1.27，文中波浪頻率均為此角頻率，由曲線可知，高航速內(nèi)傾船的峰值點(diǎn)在3.5 rad/s處，即λ/L=1.6，而同樣波浪條件下，內(nèi)傾船在高航速時(shí)的縱搖和垂蕩運(yùn)動(dòng)明顯比無(wú)航速時(shí)更劇烈.OpenFOAM和AQWA縱搖頻響計(jì)算結(jié)果的誤差最大為3%，垂蕩頻響的誤差最大為9%，均在可接受誤差范圍內(nèi)，說(shuō)明OpenFOAM計(jì)算內(nèi)傾船在縱浪中的運(yùn)動(dòng)較為準(zhǔn)確.

圖3 OpenFOAM和AQWA結(jié)果對(duì)比

2.2 大波陡非線性運(yùn)動(dòng)時(shí)域模擬

文中計(jì)算了角頻率為3.917 rad/s的迎浪規(guī)則波下，波陡H/λ=0.01,0.025及0.05時(shí)，F(xiàn)r=0，0.2及0.367的內(nèi)傾船的縱向運(yùn)動(dòng)時(shí)歷圖，通過(guò)對(duì)比分析各種時(shí)歷圖來(lái)探究?jī)?nèi)傾船在不同波陡的迎浪規(guī)則波中運(yùn)動(dòng)的非線性.

圖4為是零航速時(shí)，不同波陡對(duì)縱搖的影響.由圖4可知，波陡0.01時(shí)，上下幅值相等，隨著波陡增大，上下幅值差距越來(lái)越大，波陡0.05時(shí)，零線以上幅值高出零線以下幅值15%左右，這是因?yàn)榇├耸着潘w積小，而且首柱內(nèi)傾，容易發(fā)生埋首，同時(shí)尾部出水，水線面面積急劇減小，使得復(fù)原力矩大幅減小，尤其是在波陡0.05時(shí)，發(fā)生埋首和甲板上浪現(xiàn)象，致使向首部?jī)A斜的幅值更大；相反，當(dāng)模型向尾部?jī)A斜時(shí)，雖然也會(huì)出現(xiàn)首部出水情況，但水線面面積的變化相對(duì)于首部上浪時(shí)要小得多，因此，向尾部?jī)A角度較向首部?jī)A角度小，這是縱搖復(fù)原力矩的不對(duì)稱性導(dǎo)致內(nèi)傾船縱搖運(yùn)動(dòng)的非線性.

圖5為零航速時(shí)，各波陡下的重心處的垂向加速度時(shí)歷圖.波陡從0.01到0.05，零線上下幅值差距越來(lái)越大，零線以下幅值高出零線以上幅值最大23%左右，除此之外，隨著波陡的升高，上部肥短，下部尖瘦的現(xiàn)象明顯，垂向加速度的非線性主要是由兩方面原因共同引起的，一方面是船舶干舷內(nèi)傾，垂向復(fù)原力的非線性；另一方面是在縱搖過(guò)程中，首部埋首和上浪，尾部出水和抨擊等非線性現(xiàn)象，導(dǎo)致了垂向加速度的非線性，而重心處的垂向加速度是船模各點(diǎn)垂向加速度最小的地方，它的非線性也間接反映了垂蕩運(yùn)動(dòng)的非線性.

圖4 Fr=0時(shí)波陡對(duì)縱搖的影響

圖5 Fr=0時(shí)波陡對(duì)重心處垂向加速度的影響

圖6為低航速時(shí)，各波陡下的縱搖時(shí)歷圖.和零航速時(shí)一樣，波陡為0.01時(shí)，上下幅值相等.而隨著波陡的增加，上下幅值差距逐漸加大，波陡為0.05時(shí)，零線以上幅值高出零線以下幅值75%左右，一方面有航速內(nèi)傾船在航行時(shí)會(huì)有一定的初始首傾角，另一方面是因?yàn)樵谟泻剿贂r(shí)，船舶發(fā)生埋首，帶有速度的水流沖上甲板，沖擊上層建筑，加劇了船體向首部的傾斜角度；當(dāng)船模向尾部?jī)A斜時(shí)，方形的尾部改善了尾部的流場(chǎng)，水流順著船底流出，減小了尾部?jī)A斜角度.首尾不同的流場(chǎng)狀態(tài)加劇了零線上下幅值的差距，使低航速內(nèi)傾船縱搖的非線性更加明顯.

圖6 Fr=0.2時(shí)波陡對(duì)縱搖的影響

圖7是低航速時(shí)，各波陡下的重心處垂向加速度時(shí)歷圖.與零航速相比，有航速船模重心處的垂向加速度曲線明顯出現(xiàn)了上部肥大，下部尖廋的現(xiàn)象，波陡為0.05時(shí)，零線以下幅值高出零線以上幅值46%左右，同時(shí)曲線局部出現(xiàn)波動(dòng)現(xiàn)象，造成這種非線性現(xiàn)象，除了干舷內(nèi)傾致使的垂向回復(fù)力的非線性之外，首柱內(nèi)傾和方尾導(dǎo)致船體首部和尾部復(fù)雜的流場(chǎng)狀態(tài)也有影響.

圖8為高航速時(shí)，各波陡下的縱搖時(shí)歷圖.波陡從0.01開(kāi)始，上下幅值已經(jīng)出現(xiàn)差距，到波陡為0.05時(shí)，零線以上幅值高出零線以下幅值135%左右，說(shuō)明高航速時(shí)，船體縱搖的非線性越發(fā)明顯，但與低航速時(shí)相比，雖然向首部?jī)A斜的幅值和向尾部?jī)A斜幅值之間差值更大，但幅值大小卻小于低航速時(shí)，出現(xiàn)這種現(xiàn)象的原因是：高航速下，大量有航速水團(tuán)沖上甲板沖擊上層建筑，抑制了首部出水，減小了向尾部?jī)A斜幅值；同時(shí)，穿浪首的穿浪特性減小了埋首程度，使船體的濕表面沒(méi)有出現(xiàn)大幅減小，降低了向首傾斜幅值，當(dāng)然高航速內(nèi)傾船和低航速內(nèi)傾船相比，峰值點(diǎn)的位置發(fā)生了偏移，也是造成這種縱搖減緩的原因之一.

圖7 Fr=0.2時(shí)波陡對(duì)重心處垂向加速度的影響

圖8 Fr=0.367時(shí)波陡對(duì)縱搖的影響

圖9是高速時(shí)，各波陡下的重心處垂向加速度時(shí)歷圖.與低航速時(shí)相比，零線上下幅值差距更大，波陡為0.05時(shí)，零線以上幅值高出零線以下幅值100%左右，同時(shí)上部更加肥短，而下部也越發(fā)尖瘦，非線性更加明顯，這是干舷內(nèi)傾和甲板上帶速度水團(tuán)沖擊船體共同作用的結(jié)果.

圖10為高航速時(shí)，不同波陡下船體甲板最大上浪量圖.由圖10可知，在高航速時(shí)，甲板的最大上浪量隨著波陡的增大而增加，而船體在高航速下的縱搖幅值隨著波陡的增大而逐漸減小，呈現(xiàn)出與甲板最大上浪量相反的趨勢(shì)，這可能是因?yàn)閮?nèi)傾船以高航速航行時(shí)，大量的水團(tuán)流入甲板有效抑制了首部出水程度，導(dǎo)致縱向運(yùn)動(dòng)的劇烈程度有所改善.

圖9 Fr=0.367時(shí)波陡對(duì)重心處垂向加速度的影響

圖11為高航速下，波陡為0.05時(shí)，船體出現(xiàn)埋首和首部出水現(xiàn)象的運(yùn)動(dòng)圖.由圖11中可知，船體埋首時(shí)，尾部并無(wú)明顯的出水現(xiàn)象，上層建筑之前的部分完全入水；并且船體首部出水時(shí)，船體尾傾角度并不大，甲板水一部分在沖擊上層建筑，一部分出現(xiàn)回流，或順著內(nèi)傾干舷流出船體，水團(tuán)運(yùn)動(dòng)的非線性尤其明顯.

圖10 Fr=0.367時(shí)不同波陡下最大甲板上浪量

圖11 Fr=0.367,H/λ=0.05時(shí)運(yùn)動(dòng)圖

3 結(jié) 論

1) 內(nèi)傾船在大波陡波浪中的縱向運(yùn)動(dòng)存在明顯的非線性現(xiàn)象，這是船體內(nèi)傾和甲板上浪共同作用的結(jié)果.

2) 低航速時(shí)，內(nèi)傾船型的甲板上浪現(xiàn)象會(huì)加大船體的縱向運(yùn)動(dòng)，受非線性運(yùn)動(dòng)的影響，內(nèi)傾船的縱向運(yùn)動(dòng)會(huì)比傳統(tǒng)船型大.

3) 通過(guò)比較波陡1/20時(shí)，不同航速下內(nèi)傾船的縱搖運(yùn)動(dòng)幅值可以看出，船體的縱搖并沒(méi)有隨著航速的增加變得更加劇烈，說(shuō)明高航速內(nèi)傾船在大波陡波浪下依然具有較好的航行能力.

4) 粘性CFD軟件更適合模擬內(nèi)傾船的非線性運(yùn)動(dòng)，從高航速內(nèi)傾船在大波陡波浪中運(yùn)動(dòng)的動(dòng)態(tài)圖可以看出，OpenFOAM能準(zhǔn)確捕捉了內(nèi)傾船的埋首和甲板上浪等非線性現(xiàn)象.

本文對(duì)迎浪狀態(tài)，不同波陡的規(guī)則波下內(nèi)傾船的縱向運(yùn)動(dòng)展開(kāi)研究，發(fā)現(xiàn)船體艏柱、干舷的內(nèi)傾設(shè)計(jì)，使船體在迎浪中運(yùn)動(dòng)具有強(qiáng)烈的非線性，今后還將進(jìn)一步研究迎浪不規(guī)則波下內(nèi)傾船縱向非線性運(yùn)動(dòng)以及規(guī)則波下甲板上浪對(duì)內(nèi)傾船縱向運(yùn)動(dòng)的影響.

[1]HASHIMOTO H, UMEDA N, MATSUDA A. Broaching prediction of a wave-piercing tumblehome vessel with twin screws and twin rudders[J]. Journal of Marine Science and Technology,2011,16(4):448-461.

[2]MCCUE A, LEIGH S, CAMPBEL L, et al. On the parametric resonance of tumblehome hullform in a longitudinal seaway[J]. Naval Engineers Journal,2007,119(3):35-44.

[3]ALISTAIR G, EVANGELOS B. Parametric rolling of the tumblehome hull using CFD[C]. Proceedings of the 12th International Conference on the Stability of Ships and Ocean Vehicles, Glasgow, UK,2015.

[4]HASHIMOTO H. Pure loss of stability of a tumblehome hull in following seas[C]. Proc. of the 9th Int. Offshore and Polar Engineering Conference，Osaka, Japan,2009.

[5]張進(jìn)豐,顧民,魏建強(qiáng).低干舷隱身船波浪中縱向運(yùn)動(dòng)的模型試驗(yàn)及理論研究[J].船舶力學(xué),2009(2):169-176.

[6]賴國(guó)軍，曾凡明，劉金林，等.船舶軸系中推力軸承油膜對(duì)縱向振動(dòng)影響研究[J].武漢理工大學(xué)學(xué)報(bào)(交通科學(xué)與工程版)，2015(5)：1005-1009.

[7]劉寧，劉曉媛.迎浪工況下氣墊船運(yùn)動(dòng)響應(yīng)研究[J].武漢理工大學(xué)學(xué)報(bào)(交通科學(xué)與工程版)，2016，40(6)：978-982.

[8]白鑫，張寶吉，李棱，等.半潛船波浪載荷長(zhǎng)期預(yù)報(bào)方法研究與分析[J].武漢理工大學(xué)學(xué)報(bào)(交通科學(xué)與工程版)，2016，40(1)：190-194.

Study on the Longitudinal Nonlinear Motion of the Tumblehome Hull in Head Wave Based on CFD

LIU Le MAO Xiaofei

(SchoolofTransportation,WuhanUniversityofTechnology,Wuhan430063,China)

Compared with the traditional ship, in the highly steep wave, the vertical movement of speeding tumblehome hull may produce green water and nonlinear hydrodynamic problems. Taking the DTMB 5613 model as the research object, the CFD software OpenFOAM is used to carry out the numerical simulation of pitch and heave under the Froude numberFr=0, 0.2 and 0.367, wave steepnessH/λ=0.01,0.025 and 0.05. The results show that the longitudinal motion of the high speed tumblehome hull in the highly steep wave is strongly nonlinear. The non-linearity phenomena such as the wave on the deck caused by the tumblehome design of the hull are the main factors leading to the nonlinearity of the longitudinal motion.

tumblehome hull; pitch; heave; wave steepness; non-linearity

2017-06-02

U661.32

10.3963/j.issn.2095-3844.2017.04.030

劉樂(lè)(1990—)：男，碩士生，主要研究領(lǐng)域?yàn)榇八畡?dòng)力學(xué)