黃 磊，劉利琴，唐友剛

（天津大學(xué)水利工程仿真與安全國家重點(diǎn)實(shí)驗(yàn)室，天津 300072）

二維矩形月池內(nèi)流體自振特性研究

黃 磊，劉利琴，唐友剛

（天津大學(xué)水利工程仿真與安全國家重點(diǎn)實(shí)驗(yàn)室，天津 300072）

研究二維矩形月池流體的自振特性?；诰€性勢流理論，采用Galerkin展開，將問題轉(zhuǎn)化為月池自振頻率特征值問題，求解自振頻率及速度勢函數(shù)的系數(shù)矩陣，獲得二維矩形月池固有頻率及振型的半解析解，分析月池幾何參數(shù)對月池流體自振特性影響。研究表明，由于內(nèi)部與外流連通，月池流體呈垂向活塞（piston）振動(dòng)及晃蕩（sloshing）運(yùn)動(dòng)兩種形式。垂向活塞振動(dòng)模態(tài)頻率變化主要受液面高度影響，且隨液面高度增加而減??；晃蕩模態(tài)頻率主要受月池寬度影響，隨月池寬度增加而減小?；钊駝?dòng)模態(tài)振型在液面較淺時(shí)出現(xiàn)波峰，隨液面高度增加而消失。利用CFD方法計(jì)算月池運(yùn)動(dòng)，結(jié)果與半解析解中振型曲線擬合良好。

月池；半解析解；自振特性；勢流理論

鉆井船、FPSO及Spar平臺(tái)等海洋結(jié)構(gòu)物常據(jù)實(shí)際需要布置垂直方向貫穿結(jié)構(gòu)主體的月池結(jié)構(gòu)，立管及其它鉆井設(shè)施穿過月池通向海底，以避免遭受過大波流載荷激勵(lì)［1－3］。由于底部開口，月池內(nèi)流體運(yùn)動(dòng)存在兩類固有振蕩形式，即流體沿深度方向“活塞（piston）”振動(dòng)及流體液面左右“晃蕩（sloshing）”運(yùn)動(dòng)［4－6］。針對月池內(nèi)流體的振動(dòng)特性及復(fù)雜海洋環(huán)境中月池流體與結(jié)構(gòu)的耦合運(yùn)動(dòng)已開展相關(guān)研究。如Aalbers［7］將月池內(nèi)水柱的活塞振蕩處理為彈簧－質(zhì)量系統(tǒng)，研究波浪作用的船體發(fā)生垂向運(yùn)動(dòng)時(shí)月池內(nèi)流體活塞振動(dòng)，分析月池內(nèi)水柱活塞振蕩特性對船舶垂蕩運(yùn)動(dòng)影響。結(jié)果表明，月池能有效降低船舶垂蕩運(yùn)動(dòng)，改善其運(yùn)動(dòng)特性。Molin［8］假設(shè)無限水深靜止船舶的矩形月池，用簡化的準(zhǔn)解析法研究月池中流體活塞振蕩、晃蕩運(yùn)動(dòng)的固有振動(dòng)特性，給出液艙內(nèi)流體活塞振動(dòng)、晃蕩的一階固有頻率近似表達(dá)式，并模擬月池內(nèi)二、三維自由液面形狀。Faltinsen等［9］對有限水深條件下二維矩形月池建立線性勢流方程，通過求解域分解，利用伽遼金法研究船舶小幅輻射運(yùn)動(dòng)對月池活塞振蕩運(yùn)動(dòng)的解析解。Kristiansen等［10－12］利用CFD建立二維月池模型，研究非線性自由液面、邊界阻尼及流動(dòng)分離對月池活塞振蕩影響，結(jié)果表明流動(dòng)分離是造成線性勢流理論結(jié)果高于實(shí)驗(yàn)測量數(shù)據(jù)的主要原因。姚熊亮等［13－14］通過實(shí)驗(yàn)研究圓形月池在波浪條件下的流激振蕩特性。周華偉等［15］據(jù)勢流理論研究具有月池的圓柱結(jié)構(gòu)輻射與繞射問題。Zhang等［16］通過分解求解域方法研究兩層流中的月池運(yùn)動(dòng)問題。

本文在文獻(xiàn)［8］基礎(chǔ)上，基于線性勢流理論研究二維矩形月池的活塞振動(dòng)、晃蕩運(yùn)動(dòng)特性，分析月池幾何參數(shù)對月池流體自振特性影響，利用CFD軟件模擬月池運(yùn)動(dòng)對比解析解振型曲線，研究成果可為相關(guān)計(jì)算分析及結(jié)構(gòu)設(shè)計(jì)提供參考依據(jù)。

1 數(shù)學(xué)模型建立

月池長度遠(yuǎn)大于寬度時(shí)，長度方向?qū)υ鲁刈哉裉匦杂绊戄^小。因此可將月池簡化為二維形式?？紤]的二維矩形月池寬度b，液面高度h，月池底部完全開放。建立坐標(biāo)系，原點(diǎn)o位于月池底部左下角，見圖1。

圖1 二維矩形月池坐標(biāo)系Fig．1 Coordinate of a two dimensional rectanglemoonpool

設(shè)船舶具有無限長寬且靜止不動(dòng)，月池內(nèi)流體為無粘、無旋、不可壓縮的理想流體，據(jù)線性勢流理論，得月池內(nèi)流體控制方程為

2 數(shù)值計(jì)算結(jié)果

2.1 系數(shù)矩陣數(shù)值計(jì)算

本文利用復(fù)化辛浦生算法［17］，數(shù)值計(jì)算式（13）中二重積分得系數(shù)tmn，并將計(jì)算結(jié)果與采用Molin計(jì)算方法［8］所得結(jié)果進(jìn)行對比，見圖2。由圖2看出，兩種計(jì)算方法所得結(jié)果吻合較好。

2.2 固有振動(dòng)頻率

通過Galerkin展開將計(jì)算月池的固有振動(dòng)特性問題轉(zhuǎn)化為計(jì)算式（18）的特征值問題。截取適當(dāng)階次，即可得各階頻率近似解。其中0階代表月池的piston模態(tài)，奇數(shù)階代表反對稱sloshing模態(tài)，其余偶數(shù)階代表對稱sloshing模態(tài)。取月池寬b＝20 m、液面高h(yuǎn)分別為1 m、10 m時(shí)，不同截?cái)嚯A次N對應(yīng)的piston模態(tài)及前三階sloshing模態(tài)頻率見表1、表2。由二表看出，各階頻率收斂速度均較快。

表1 b＝20 m，h＝1 m不同截?cái)嚯A次下各階模態(tài)頻率Tab.1 b＝20 m，h＝1m，Piston and sloshingmode frequencies for d ifferent truncation order

表2 b＝20 m，h＝10m不同截?cái)嚯A次下各階模態(tài)頻率Tab.2 b＝20 m，h＝10 m，Piston and sloshing mode frequencies for different truncation order

月池寬度b分別為10 m、20 m時(shí)piston模態(tài)頻率及1、2階sloshing模態(tài)頻率隨月池液面高度h的變化曲線見圖3～圖6。

圖2 積分結(jié)果對比Fig．2 Integration value comparison

圖3 自振頻率隨h變化（b＝10 m）Fig．3 Natural frequency change versus h（b＝10 m）

圖4 自振頻率比值隨h變化（b＝10 m）Fig．4 Ratio of natural frequency change versus h（b＝10 m）

圖3為各階自振頻率隨液面高度h變化曲線。由圖3看出，Piston模態(tài)頻率為最低階模態(tài)。此階模態(tài)振動(dòng)時(shí)月池內(nèi)流體整體上下運(yùn)動(dòng)，隨液面高度增加流體質(zhì)量持續(xù)增加，自振頻率隨之降低。而各階sloshing模態(tài)主要由自由液面左右晃動(dòng)激勵(lì)時(shí)產(chǎn)生，液面高度變化不影響sloshing模態(tài)，各階sloshing曲線為水平平行；但液面高度較低時(shí)自由液面位置離月池底部較近，sloshing頻率會(huì)隨h增高而降低，由于月池底部開放內(nèi)外相連，流體流動(dòng)對自由液面運(yùn)動(dòng)產(chǎn)生耦合作用。由此知，月池水平運(yùn)動(dòng)時(shí)底部開放對月池內(nèi)流體影響范圍較小，離月池底超5 m以上的流體已基本不受底部流體流動(dòng)影響。

圖4為對應(yīng)不同液面高度h的各階自振頻率與h＝0．001 m時(shí)自振頻率比值變化曲線。由圖4可清晰看出各階頻率變化率。h較小時(shí)頻率下降速度相同，此時(shí)底部流體流動(dòng)明顯主導(dǎo)各階模態(tài)變化；而隨h升高，僅piston模態(tài)頻率持續(xù)下降，其它各階sloshing模態(tài)頻率不變化。

圖5、圖6為月池寬度20 m時(shí)兩組頻率隨h變化曲線。與圖3、圖4對比看出，月池寬度增加未改變自振頻率隨h的變化規(guī)律。由此可見，液面高度h為決定pinston模態(tài)頻率的關(guān)鍵參數(shù)，而對sloshing模態(tài)頻率影響較小。由于piston模態(tài)頻率隨h持續(xù)降低，h較高時(shí)piston模態(tài)自振頻率可能接近平臺(tái)自振頻率及波浪入射頻率，應(yīng)特別關(guān)注。圖7為液面高度1 m時(shí)月池各階自振頻率隨月池寬度變化曲線。由圖7看出，月池寬度從1 m增長至15 m時(shí)各階頻率均逐漸減小，超過15 m后頻率變化趨于平穩(wěn)。對比各曲線知，sloshing模態(tài)頻率隨寬度變化率高于piston模態(tài)。圖8為各高度頻率與液面1 m高頻率比值變化曲線。由圖8看出，各階sloshing模態(tài)頻率隨月池寬度變化率相同，均高于piston模態(tài)頻率。

圖5 自振頻率隨h變化（b＝20 m）Fig．5 Natural frequency change h（b＝20 m）

圖6 自振頻率比值隨h變化（b＝20m）Fig．6 Ratio of natural frequency change versus versus h（b＝20 m）

圖7 自振頻率隨b變化（h＝1m）Fig．7 Natural frequency change versus b（h＝1m）

圖10 自振頻率比值隨b變化（h＝10 m）Fig．10 Ratio of natural frequency change versus b（h＝10m）

圖8 自振頻率比值隨b變化（h＝1 m）Fig．8 Ratio of natural frequency change versus b（h＝1m）

圖9 自振頻率隨b變化（h＝10 m）Fig．9 Natural frequency change versus b（h＝10m）

圖7、圖8中月池液面較淺，寬度高于液面高度，該形式月池在鉆井船、FPSO等結(jié)構(gòu)中較常見。而對深水SPAR平臺(tái)月池，其液面較高與月池長度相當(dāng)，甚至高于月池長度。圖9～圖12為液面高度10 m、20 m時(shí)自振頻率隨月池寬度變化曲線。由圖9看出，b＜10 m時(shí)sloshing模態(tài)頻率下降較明顯；b＞10 m后隨寬度繼續(xù)增加sloshing模態(tài)頻率變化逐漸平緩。piston模態(tài)曲線則始終變化較小，曲線形狀較平緩。對比圖11可知，液面高度20 m時(shí)自振頻率隨月池寬度變化與圖9中液面高度10 m曲線具有完全相同規(guī)律。圖10為不同月池寬度時(shí)各階自振頻率與b＝1 m時(shí)比值變化曲線。由圖11看出，各階sloshing頻率隨寬度增加變化規(guī)律完全相同，且變化速度由快變慢。對比圖12看出，液面高度增加未影響此變化規(guī)律。

由圖9～圖12可知，月池寬度對sloshing模態(tài)頻率影響非常明顯，此為sloshing模態(tài)由自由液面運(yùn)動(dòng)所致。月池寬度變化對自由液面形狀產(chǎn)生直接影響，因此自振模態(tài)頻率發(fā)生明顯變化。液面較低時(shí)自由液面離底部較近，其運(yùn)動(dòng)會(huì)受由底面流體流動(dòng)影響，但范圍有限。而對垂直方向振動(dòng)的piston模態(tài)頻率，基本不受月池寬度變化影響。

2.3 自由液面

將由式（18）所得各階特征值對應(yīng)的特征向量矩陣B代入式（19）得矩陣A。將矩陣A、B代入式（9）及自由液面方程

式中：η為波面升高。

由式（20）可得各階自振頻率對應(yīng)振型。

圖13～圖15為20 m寬月池不同液面高度時(shí)piston模態(tài)及1～3階sloshing模態(tài)對應(yīng)振型。由三圖看出，piston模態(tài)呈活塞運(yùn)動(dòng)特性，波面整體升高偏離平衡位置；sloshing模態(tài)奇數(shù)階為反對稱振型，偶數(shù)階為對稱振型，波面曲線在平衡位置具有三角函數(shù)特征。

圖11 自振頻率隨b變化（h＝20m）Fig．11 Natural frequency change versus b（h＝20 m）

圖12 自振頻率比值隨b變化（h＝20m）Fig．12 Ratio of natural frequency change versus b（h＝20 m）

圖13 b＝20m，h＝0．02 m時(shí)月池各階模態(tài)振型Fig．13 Free surface shape for different orders when b＝20 m，h＝0．02 m

圖14 b＝20 m，h＝2 m時(shí)月池各階模態(tài)振型Fig．14 Free surface shape for different orderswhen b＝20 m，h＝2m

圖15 b＝h＝20 m時(shí)月池各階模態(tài)振型隨h變化Fig．15 Free surface shape for different orderswhen b＝h＝20m

圖16 b＝20 m時(shí)piston模態(tài)振型Fig．16 Free surface shape of piston mode for different h when b＝20 m

對比圖13～圖15看出，不同液面高度振型具有明顯差異。h＝0．02 m時(shí)液面較低，sloshing模態(tài)各階振型在接近月池兩側(cè)壁處出現(xiàn)波峰，piston模態(tài)在月池寬度方向出現(xiàn)波峰。當(dāng)液面高度達(dá)到2 m時(shí)波峰逐漸減小，低階頻率波峰仍存在，而高階頻率此現(xiàn)象已消失。該現(xiàn)象因底部流體流動(dòng)所致。piston模態(tài)液面較低，液柱在垂向振動(dòng)時(shí)，液面向上或向下會(huì)出現(xiàn)波峰，而高液面時(shí)，液面完全呈平面，在平衡位置上下振動(dòng)。sloshing模態(tài)隨液面高度增加底面流動(dòng)對其影響衰減較迅速。底面流動(dòng)對自由液面形狀影響隨h減小速度高于對頻率影響減小速度。

圖17 b＝20 m時(shí)不同h一階sloshing模態(tài)振型Fig．17 Free surface shape for 1st orders sloshingmode for different h

圖18 b＝20 m時(shí)不同h二階sloshing模態(tài)振型Fig．18 Free surface shape for 2nd orders sloshingmode for different h

圖19 CFD模型Fig．19 CFD model

寬20 m、不同液面高度時(shí)piston模態(tài)振型曲線見圖16。該圖進(jìn)一步表明，液面越低piston模態(tài)在月池寬度方向波峰越大。液面超過12 m時(shí)波峰基本消失，月池內(nèi)流體垂向振動(dòng)類似于剛體運(yùn)動(dòng)。月池寬20 m時(shí)不同液面高度對應(yīng)二階sloshing模態(tài)振型曲線見圖17、圖18。兩圖表明，sloshing模態(tài)一、二階振型接近月池兩側(cè)壁處出現(xiàn)波峰，隨液面升高波峰逐漸降低；液面高于2 m時(shí)該波峰基本消失。

2.4 數(shù)值仿真對比

利用數(shù)值仿真方法計(jì)算月池自振并對比與分析，以驗(yàn)證結(jié)果的有效性。計(jì)算所用CFD軟件FLUENT為基于非結(jié)構(gòu)化網(wǎng)格的有限體積法，利用動(dòng)網(wǎng)格技術(shù)模擬邊界運(yùn)動(dòng)。在其前處理模塊Gambit中建立計(jì)算模型、劃分網(wǎng)格。月池模型為二維模型，包括月池內(nèi)流體、月池外流體、船體，月池內(nèi)外流體相聯(lián)通。月池及兩側(cè)船體結(jié)構(gòu)均寬20 m，計(jì)算區(qū)域總寬160 m。用非結(jié)構(gòu)性三角形單元進(jìn)行網(wǎng)格劃分，對月池內(nèi)流體網(wǎng)格局部細(xì)化，以便對月池內(nèi)流體運(yùn)動(dòng)準(zhǔn)確刻畫。計(jì)算模型見圖19。

據(jù)sloshing模態(tài)一階自振頻率計(jì)算一階周期T1＝2π／ω1，對液面高2 m、10 m時(shí)模型施加周期T1的水平簡諧激勵(lì)。月池內(nèi)流體運(yùn)動(dòng)響應(yīng)穩(wěn)定后取體積分?jǐn)?shù)0．5的自由液面曲線，與一階振型進(jìn)行對比。CFD計(jì)算液面高2 m曲線與解析解液面高1 m、2 m的一階sloshing模態(tài)振型曲線對比見圖20。由圖20看出，解析解與CFD具有相同特征?？拷鼉蓚?cè)壁處出現(xiàn)峰值。CFD模擬液面高2 m時(shí)結(jié)果與解析解中液面高1 m的曲線穩(wěn)合程度更高。在用勢流理論求解時(shí)設(shè)流體為無粘理想流體，而在CFD方法模擬中，流體粘性對側(cè)壁處流體也會(huì)產(chǎn)生影響，造成側(cè)壁處波峰較明顯。但粘性對計(jì)算結(jié)果影響相對整體振型尺度較小，可忽略。液面高10 m時(shí)CFD計(jì)算曲線與解析解一階sloshing模態(tài)振型曲線對比見圖21。由圖21看出，兩條曲線擬合較好。

圖20 Sloshing模態(tài)解析解與h＝2 m時(shí)CFD自由液面計(jì)算結(jié)果對比Fig．20 Free surface shape of analytical solution and CFD solution when h＝2 m in sloshingmode

圖22 Piston模態(tài)解析解與h＝10 m時(shí)CFD自由液面計(jì)算結(jié)果對比Fig．22 Free surface shape of analytical solution and CFD solution in pistonmode

圖21 Sloshing模態(tài)解析解與h＝10m時(shí)CFD自由液面計(jì)算結(jié)果對比Fig．21 Free surface shape of analytical solution and CFD solution when h＝2m in sloshingmode

據(jù)piston模態(tài)自振頻率計(jì)算垂蕩周期T0＝2π／ω0，對液面高2 m、10 m時(shí)模型施加周期T0的垂蕩方向簡諧激勵(lì)。月池內(nèi)流體運(yùn)動(dòng)響應(yīng)穩(wěn)定后取體積分?jǐn)?shù)0．5的自由液面曲線與piston模態(tài)振型進(jìn)行對比。CFD計(jì)算及解析解自由液面對比曲線見圖22。由圖22看出，CFD結(jié)果與解析解結(jié)果擬合較好，CFD計(jì)算中自由液面在高2 m時(shí)出現(xiàn)明顯波峰，隨液面升高其波峰減小趨于平面，與解析解呈現(xiàn)相同變化趨勢。對比圖20～圖22各曲線可知，利用勢流理論研究底部開放月池的振動(dòng)特性，求解效率高，結(jié)果精度能滿足要求。

3 結(jié) 論

本文利用線性勢流理論建立月池計(jì)算模型，研究其自振頻率及自由液面。對底板完全開放的月池，自振特性顯現(xiàn)垂向piston模態(tài)及水平方向sloshing模態(tài)。結(jié)論如下：

（1）月池流體有piston及sloshing兩種運(yùn)動(dòng)模態(tài)，其中piston模態(tài)頻率低于sloshing模態(tài)各階固有頻率，sloshing模態(tài)振型的奇數(shù)階為反對稱振型，偶數(shù)階為對稱振型。

（2）月池液面高度變化對piston模態(tài)頻率影響較大：液面增度piston模態(tài)頻率降低，液面足夠高時(shí)piston模態(tài)頻率可能接近結(jié)構(gòu)的自振頻率及波浪入射頻率，引起結(jié)構(gòu)大的振動(dòng)；月池寬度變化對sloshing模態(tài)頻率影響較大，b＜10時(shí)sloshing模態(tài)頻率下降較明顯。

（3）液面高度h較小時(shí)sloshing模態(tài)振型在接近月池兩側(cè)壁處出現(xiàn)波峰，piston模態(tài)在月池寬度方向出現(xiàn)波峰。液面增度波峰逐漸降低并消失，此時(shí)月池內(nèi)流體的垂向振動(dòng)為剛體運(yùn)動(dòng)。

（4）CFD結(jié)果與半解析解piston模態(tài)及sloshing模態(tài)振型變化規(guī)律相同，顯示半解析解計(jì)算精度良好。

［1］董艷秋．深海采油平臺(tái)波浪載荷及響應(yīng)［M］．天津：天津大學(xué)出版社，2005．

［2］Drobyshevski Y．Hydrodynamic coefficients of a twodimensional，truncated rectangular floating structure in shallow water［J］．Ocean Engineering，2004，31：305－341．

［3］Wei Y，Yang J，Chen G，et al．Experimental study on the hydrodynamic performance of FDPSO and SRV［J］．Ships and Offshore Structures，2012，7（4）：357－369．

［4］Faltinsen O M，Timokha A．Sloshing［M］．London：Cambridge University Press，2009．

［5］Miles JW．On the eigenvalue problem for fluid sloshing in a half-space［J］．Journal of Applied Mathematics and Physics，1972，23（6）：861－869．

［6］Gupta H，Blevins R，Banon H．Effect of moonpool hydrodynamics on spar heave［C］．ASME，2008．

［7］Aalbers A B．The water motions in a moonpool［J］．Ocean Engineering，1984，11（6）：557－579．

［8］Molin B．On the piston and sloshingmodes inmoonpools［J］．Journal of Fluid Mechanics，2001，430（1）：27－50．

［9］Faltinsen O M，Rognebakke O F，Timokha A N．Twodimensional resonant piston-like sloshing in a moonpool［J］．Journal of Fluid Mechanics，2007，575：359－397．

［10］Kristiansen T，F(xiàn)altinsen O M．Application of a vortex tracking method to the piston-like behaviour in a semientrained vertical gap［J］．Applied Ocean Research，2008，30（1）：1－16．

［11］Kristiansen T，F(xiàn)altinsen O M．A two-dimensional numerical and experimental study of resonant coupled ship and pistonmodemotion［J］．Applied Ocean Research，2010，32（2）：158－176．

［12］Kristiansen T，F(xiàn)altinsen O M．Gap resonance analyzed by a new domain-decomposition method combining potential and viscous flow Draft［J］．Applied Ocean Research，2012，34：198－208．

［13］姚熊亮，康莊．圓形月池流激振蕩實(shí)驗(yàn)研究［J］．力學(xué)學(xué)報(bào)，2007，39（3）：333－342．

YAO Xiong-liang，KANG Zhuang．Experimental research on flow induced oscillation of circle moonpool［J］．Chinese Journal of Theoretical and App lied Mechanics，2007，39（3）：333－342．

［14］康莊，姚熊亮．浪流條件下月池內(nèi)流體活塞振蕩特性實(shí)驗(yàn)研究［J］．哈爾濱工程大學(xué)學(xué)報(bào)，2008，29（3）：209－216．

KANG Zhuang，YAO Xiong-liang．Experimental research on piston oscillating behavior of fluid in a moonpool with waveflow［J］．Journal of Harbin Engineering University，2008，29（3）：209－216．

［15］Zhou Hua-wei，Zhang Hong-sheng．Radiation and diffraction analysis of a cylindrical body with amoon pool［J］．Journal of Hydrodynamics，Ser．B，2013，25（2）：196－204．

［16］Zhang Xin-shu，Bandyk Piotr．On two-dimensionalmoonpool resonance for twin bodies in a two-layer fluid［J］．Applied Ocean Research，2013，40：1－13．

［17］翟瑞彩，謝偉松．?dāng)?shù)值分析［M］．天津：天津大學(xué)出版社，2000．

Natural vibration characteristics of fluid in a two dimensional rectangular moonpool

HUANG Lei，LIU Li-qin，TANG You-gang
（State Key Laboratory of Hydraulic Engineering Simulation and Safety，Tianjin University，Tianjin 300072，China）

Based on the linear potential theory，the equations of fluid motions in a two dimensional rectangular moonpoolwere converted into an eigenvalue problem by Galerkin method．Via solving the coefficientmatrix of potential function，the semi-analytical solution of natural frequencies and shapes of fluid was obtained and the influences of geometric parameters of themoonpool on its natural vibration characteristics were discussed．The results show that there exist two natural vibration modes of fluid in themoonpool，namelt，the piston mode and sloshing mode．The piston mode frequency in vertical direction decreases as the liquid height increases and ismainly determined by the liquid height，whereas sloshingmode frequencies decreases as themoonpool width increases and mainly depend on themoonpoolwidth．In the piston mode shape crestappears as the liquid depth is small，and vanishes as the liquid depth increases．The CFD method was used to verify the semi-analytical solution and the free surface curves in the two results are in good coincidence．

moonpool；semi-analytical solution；natural vibration characteristics；potential theory

P751

：A

10．13465／j．cnki．jvs．2014．22．025

國家自然科學(xué)基金資助項(xiàng)目（51179125）

2013－05－15 修改稿收到日期：2013－11－21

黃磊男，博士生，1981年生

劉利琴女，博士，副教授，1977年生