水下樁腿長(zhǎng)度變化對(duì)浮式結(jié)構(gòu)固有頻率的影響研究

武海浪，陳徐均，黃亞新，劉俊誼

(解放軍理工大學(xué) 野戰(zhàn)工程學(xué)院，江蘇 南京 210007)

水下樁腿長(zhǎng)度變化對(duì)浮式結(jié)構(gòu)固有頻率的影響研究

武海浪，陳徐均，黃亞新，劉俊誼

(解放軍理工大學(xué) 野戰(zhàn)工程學(xué)院，江蘇 南京 210007)

為了分析帶支腿浮式結(jié)構(gòu)水下樁腿處于不同的長(zhǎng)度時(shí)振動(dòng)固有頻率的變化情況，采用有限元方法和三維線性水彈性理論對(duì)一理想的帶支腿浮式結(jié)構(gòu)進(jìn)行研究，分工況計(jì)算了3種彈性模態(tài)的干、濕固有頻率。結(jié)果表明對(duì)于同一階彈性模態(tài)來(lái)說(shuō)，結(jié)構(gòu)主體平臺(tái)振型是決定水動(dòng)力系數(shù)大小的主要因素，樁腿水下長(zhǎng)度是次要因素，而振型又受到樁腿長(zhǎng)度的制約；樁腿水下長(zhǎng)度的改變對(duì)帶支腿浮式結(jié)構(gòu)的干模態(tài)固有頻率以及濕固有頻率都有較明顯地調(diào)制作用。對(duì)于干模態(tài)固有頻率來(lái)說(shuō)，樁腿關(guān)于船身越趨于對(duì)稱分布時(shí)結(jié)構(gòu)的共振頻率越高；對(duì)于濕固有頻率來(lái)說(shuō)，水下樁腿長(zhǎng)度地增加有迫使共振頻率趨向低頻方向移動(dòng)趨勢(shì)。

海洋工程；水彈性力學(xué)；帶支腿浮式結(jié)構(gòu)；固有頻率；模態(tài)分析；樁腿長(zhǎng)度

Abstract:To analyze the dry and wet natural frequencies of the floating structure with legs,the hydrodynamic coefficients and the stiffness coefficients of the floating structure with legs are calculated by utilizing the Finite Element Method and the 3 D linear hydroelasticity theory.The results show the main factor that determines the hydrodynamic coefficients of the identical elastic mode is the mode shape of the main platform,and the length of the legs is the secondary factor,however,the mode shape of the main platform is restricted by the legs.The variation of the leg's length underwater can prominently modulate both the dry natural frequencies and the wetted natural frequencies of the floating structure with legs.In the dry condition,with the legs turning to the symmetrical distribution about the main body,the natural frequency of the structure tends to move to the high frequency domain.In the wet condition,the modulation effects of the legs on the natural frequencies are manifested by adjusting the added mass.With the legs stretching out from the bottom of the main body,the wetted frequencies have the tendency to move to the low frequency domain.

Keywords:ocean engineering; hydroelasticity; floating structure with legs; natural frequency; model analysis; length of the legs

帶支腿浮式結(jié)構(gòu)是指在浮式平臺(tái)的基礎(chǔ)上設(shè)計(jì)加裝一定數(shù)量的可升降支腿而形成的新型工程結(jié)構(gòu)。該工程結(jié)構(gòu)通過(guò)升降裝置的動(dòng)作，平臺(tái)主體或樁腿可垂直升降。如自升式石油鉆井平臺(tái)和自升(自航)式風(fēng)電安裝船就是典型的帶支腿浮式結(jié)構(gòu)。國(guó)內(nèi)外學(xué)者對(duì)這種工程結(jié)構(gòu)有較多的研究[1-10]，這些研究主要是針對(duì)該結(jié)構(gòu)處于支承狀態(tài)時(shí)所展開的。

帶支腿浮式結(jié)構(gòu)在放/提樁過(guò)程中，樁腿既未觸及海底也未完全收起，結(jié)構(gòu)始終處于漂浮狀態(tài)，樁腿的存在將導(dǎo)致浮式結(jié)構(gòu)的水彈性性質(zhì)與普通浮式結(jié)構(gòu)不同。一方面，樁腿的收放會(huì)導(dǎo)致結(jié)構(gòu)質(zhì)量分布發(fā)生變化，這必然使結(jié)構(gòu)同一種干模態(tài)的振型和固有頻率也會(huì)發(fā)生變化；另一方面，水下樁腿的伸長(zhǎng)或者縮短會(huì)使得船體的重心位置、轉(zhuǎn)動(dòng)慣量，濕面積等船體自身物理屬性隨時(shí)發(fā)生改變，這樣流場(chǎng)中速度勢(shì)的定解問(wèn)題中的物面條件會(huì)發(fā)生改變導(dǎo)致其解會(huì)發(fā)生改變。最終會(huì)引起三維線性水彈性方程中的每一個(gè)參數(shù)發(fā)生相應(yīng)地調(diào)整和變化。在這個(gè)過(guò)程中，帶支腿浮式結(jié)構(gòu)的濕固有頻率必然會(huì)發(fā)生相應(yīng)地調(diào)整。

固有頻率是浮式結(jié)構(gòu)重要的屬性，一直受到關(guān)注。殷玉梅等[11]推導(dǎo)了船舶上層建筑整體縱向振動(dòng)固有頻率預(yù)報(bào)公式，將上層建筑縱向振動(dòng)固有頻率視為由上層建筑根部剛性固定在船體上的剪彎振動(dòng)固有頻率和上層建筑根部彈性固定在主船體上的剛體回轉(zhuǎn)振動(dòng)固有頻率兩部分串聯(lián)合成。為了使船舶在固有頻率變化時(shí)也有好的減搖效果，金鴻章等[12]研究了雙水艙減搖系統(tǒng)。劉文璽等[13]對(duì)系泊浮體的固有頻率、各種工況下大幅非線性運(yùn)動(dòng)、纜索張力等進(jìn)行了計(jì)算，并進(jìn)行了對(duì)比研究。王顯正等[14]用一維梁有限元方法和三維有限元方法計(jì)算3艘實(shí)船總振動(dòng)固有頻率，對(duì)計(jì)算結(jié)果進(jìn)行統(tǒng)計(jì)分析，提出對(duì)一維梁有限元方法計(jì)算結(jié)果地修正，該方法能有效改進(jìn)一維梁有限元計(jì)算方法，可快速準(zhǔn)確地預(yù)報(bào)船舶總振動(dòng)固有頻率。

樁腿在水下長(zhǎng)度的改變而引起浮體濕固有頻率變化的研究尚不多見[15-17]，這里采用有限元方法和線性三維水彈性理論對(duì)一理想化的帶支腿浮式結(jié)構(gòu)進(jìn)行分析，計(jì)算了該浮式結(jié)構(gòu)水下樁腿處于不同的長(zhǎng)度時(shí)二節(jié)點(diǎn)垂彎振動(dòng)、三節(jié)點(diǎn)垂彎振動(dòng)、一節(jié)點(diǎn)扭轉(zhuǎn)振動(dòng)時(shí)固有頻率，重點(diǎn)分析了引起浮式結(jié)構(gòu)固有頻率變化的機(jī)理和原因。

1 帶支腿浮式結(jié)構(gòu)固有頻率分析原理

1.1干固有頻率計(jì)算

用有限元方法，任意形式結(jié)構(gòu)的運(yùn)動(dòng)方程可表示為[18]：

式中：M、C、K分別為浮體總質(zhì)量矩陣、總阻尼矩陣和總剛度矩陣；U為節(jié)點(diǎn)位移向量；P為結(jié)構(gòu)分布力。

忽略外部荷載以及阻尼得到系統(tǒng)的自由振動(dòng)方程：

假設(shè)式(2)中的自由振動(dòng)為頻率為ω的簡(jiǎn)諧振動(dòng)，這樣便得到：

每一個(gè)特征值都有一個(gè)特征向量與之相對(duì)應(yīng)Dr，代表第r階模態(tài)：

1.2濕固有頻率計(jì)算

假定船體結(jié)構(gòu)周圍的流體為理想、不可壓縮流體，波浪振幅值為小量，可得自由狀態(tài)下浮體的線性水彈性方程[19]：

若令廣義輻射力有如下表達(dá)式：

式中：n是物體濕表面上的單位法向量；ur為該點(diǎn)第r階位移振型；φk(x,y,z)表示輻射波速度勢(shì)的空間分量，它表示當(dāng)浮體在平衡位置附近以入射波圓頻率ω，以第k階干模態(tài)振型作單位主坐標(biāo)幅值振蕩時(shí)所誘導(dǎo)的流場(chǎng)周圍的流體運(yùn)動(dòng)速度勢(shì)；pk表示結(jié)構(gòu)動(dòng)力學(xué)控制方程中的主坐標(biāo)響應(yīng)復(fù)數(shù)幅值；ρ為流體密度；i是虛數(shù)單位。這樣，廣義輻射力可以表示為：

式中：Trk表示浮體以頻率ω、第k階模態(tài)作單位幅值的簡(jiǎn)諧運(yùn)動(dòng)時(shí)，浮體所受到的第r階輻射力。它可以表示為：

其中，

式中：Re，Im分別表示取復(fù)數(shù)的實(shí)部和虛部，系數(shù)Ark與加速度同相位，Brk與速度同相位，所以它們分別稱為附加質(zhì)量和附加阻尼系數(shù)。

恢復(fù)力系數(shù)的表達(dá)式：

式中:Crk表示廣義恢復(fù)力系數(shù)。當(dāng)r≤6以及k≤6時(shí)，上式即退化為剛體運(yùn)動(dòng)的恢復(fù)力系數(shù)矩陣，其余情況則為彈性變形模態(tài)和剛體運(yùn)動(dòng)的耦合交叉項(xiàng)。

當(dāng)帶支腿浮式結(jié)構(gòu)處于波浪中，對(duì)于其垂蕩、橫搖、縱搖以及每一階彈性模態(tài)，當(dāng)其作給定頻率ω的運(yùn)動(dòng)時(shí)都對(duì)應(yīng)有其本身的濕固有頻率ωwr，該頻率可由如下特征方程求解：

其中，arr和crr為結(jié)構(gòu)廣義質(zhì)量矩陣和廣義剛度矩陣對(duì)角元素；Arr(ω)為浮體以頻率ω作振蕩運(yùn)動(dòng)的附加質(zhì)量矩陣對(duì)角元素；Crr為靜水恢復(fù)力矩陣對(duì)角元素。

2 帶支腿浮式結(jié)構(gòu)特性分析

2.1有限元建模及工況劃分

圖1 梁?jiǎn)卧獛е雀∈浇Y(jié)構(gòu)Fig.1 The floating structure with legs formed by beam

將船身簡(jiǎn)化為水平三維梁是水彈性分析的常用方法[20-26]。文中將主體平臺(tái)簡(jiǎn)化為水平梁的組合，同時(shí)，將船身上攜帶的樁腿簡(jiǎn)化為垂直三維梁。為了將樁腿合理的加裝到主體平臺(tái)上，具體方法及相關(guān)參數(shù)：選取某型帶支腿浮式結(jié)構(gòu)為研究對(duì)象，結(jié)構(gòu)長(zhǎng)度約為130 m，寬38 m，型深8 m。6個(gè)升降腿分布于全船兩側(cè)，每個(gè)樁腿高68 m，截面為4 m×4 m方形。如圖1，將船身簡(jiǎn)化為水平梁組合，采用有限元分析軟件建模。將代表主體平臺(tái)的梁?jiǎn)卧胖糜谒矫嬷?，若不考慮樁腿的存在，主體平臺(tái)重心始終位于靜水面上，平衡坐標(biāo)系的原點(diǎn)位置。

按樁腿長(zhǎng)度對(duì)所建立的模型分為5種工況，如表1所示。將工況5設(shè)定為工況1的對(duì)稱形式，其上部樁腿的長(zhǎng)度恰好為工況1上部樁腿的長(zhǎng)度，而工況5下部樁腿的長(zhǎng)度為工況1上部樁腿的長(zhǎng)度，這樣重心的位置就與工況1的重心位置關(guān)于靜水面相對(duì)稱。而轉(zhuǎn)動(dòng)慣量和工況1的情況完全一致。工況2和工況4也是一組對(duì)稱工況，工況3獨(dú)立于工況1、5和工況2、4自成一類。

表1 帶支腿浮式結(jié)構(gòu)工況設(shè)置Tab.1 The set of the cases

算例只是一個(gè)理想的帶支腿浮式結(jié)構(gòu)，實(shí)際帶支腿浮式結(jié)構(gòu)的配重分布未必能滿足上述要求。但本文著重于帶支腿浮式結(jié)構(gòu)濕固有頻率變化的機(jī)理解釋，故選擇了這樣的一種建模方法來(lái)進(jìn)行分析。按照上面兩兩一組的劃分工況好處在于，當(dāng)進(jìn)行干模態(tài)分析的時(shí)候，由于工況相對(duì)稱，真空中帶支腿浮式結(jié)構(gòu)的相同的模態(tài)，其振型和固有頻率必定一致。當(dāng)對(duì)兩種結(jié)果相同的干模態(tài)結(jié)果放到波浪中進(jìn)行水彈性分析時(shí)，由于同一組中的兩種工況的水下樁腿長(zhǎng)度是不一樣的，可能會(huì)導(dǎo)致結(jié)果的差異。對(duì)比這種差異性結(jié)果就可以有針對(duì)性的分析樁腿對(duì)帶支腿浮式結(jié)構(gòu)濕固有頻率的影響。增強(qiáng)了結(jié)果的可對(duì)比性，有利于對(duì)水彈性理論中濕固有頻率計(jì)算機(jī)制的認(rèn)識(shí)。

2.2模態(tài)分析

主要考察樁腿的存在對(duì)主體平臺(tái)振型和頻率的影響，暫時(shí)不考慮樁腿本身振動(dòng)而產(chǎn)生的彈性模態(tài)，故認(rèn)為樁腿的剛度遠(yuǎn)遠(yuǎn)大于船身剛度。分別提取了帶支腿浮式結(jié)構(gòu)三種振動(dòng)：二節(jié)點(diǎn)垂直彎曲；三節(jié)點(diǎn)垂直彎曲；一節(jié)點(diǎn)扭轉(zhuǎn)振動(dòng)。如表2所示，分別給出了樁腿處于5種工況下的振型和固有頻率。

表2 帶支腿浮式結(jié)構(gòu)干模態(tài)振型及固有頻率Tab.2 The mode shape and natural frequencies of the floating structure with legs

在一個(gè)模態(tài)振型下，當(dāng)樁腿向下伸長(zhǎng)，由工況1轉(zhuǎn)為工況5，浮式結(jié)構(gòu)干模態(tài)固有頻率有先增大后減小的趨勢(shì)。如在二節(jié)點(diǎn)垂彎時(shí)，工況1的固有頻率為1.459 41 Hz，工況2時(shí)為1.529 52 Hz，到工況3時(shí)固有頻率最大，為1.554 20 Hz，工況4時(shí)為1.529 52 Hz，工況5為1.459 41 Hz。三節(jié)點(diǎn)垂彎和一節(jié)點(diǎn)扭轉(zhuǎn)振動(dòng)也有類似的規(guī)律，這表明樁腿的存在對(duì)帶支腿浮式結(jié)構(gòu)的垂彎振動(dòng)固有頻率確實(shí)存在著影響，且樁腿越是關(guān)于主體平臺(tái)對(duì)稱分布，上述三種振動(dòng)的固有頻率就會(huì)越大，這反應(yīng)了樁腿對(duì)帶支腿浮式結(jié)構(gòu)干模態(tài)固有頻率地調(diào)制作用。

工況1和工況5，工況2和工況4的垂彎模態(tài)振型是相互對(duì)稱的，固有頻率相一致，振型的幅度也是一致的。這樣將工況1、工況2和工況3的垂彎歸一化振型提出來(lái)，放在一張圖中進(jìn)行對(duì)比，如圖2所示，每一幅圖中都有三種工況的帶支腿浮式結(jié)構(gòu)的振型。

對(duì)于二節(jié)點(diǎn)垂彎，如圖2(a)，對(duì)比三種工況主體平臺(tái)振型可知，工況3主體平臺(tái)的垂彎程度明顯比工況2和工況1的振幅更大，工況1的垂彎振幅為這三種振型中最小的工況?？梢姸?jié)點(diǎn)垂彎振動(dòng)中，樁腿越是趨于對(duì)稱，不僅固有頻率越大，而且主體平臺(tái)的歸一化振幅也越大。同樣的規(guī)律也存在于一節(jié)點(diǎn)扭轉(zhuǎn)振動(dòng)和三節(jié)點(diǎn)垂彎振動(dòng)當(dāng)中，如圖2(b)和圖2(c)所示。

圖2 模態(tài)振型對(duì)比圖Fig.2 Comparison of the mode shapes

3 水動(dòng)力建模及結(jié)果分析

3.1濕面網(wǎng)格劃分

已有的三維線性水彈性計(jì)算程序采用了單對(duì)稱復(fù)合勢(shì)的方法對(duì)水動(dòng)力參數(shù)進(jìn)行求解。如果所研究的浮式結(jié)構(gòu)是關(guān)于x或者y軸的對(duì)稱結(jié)構(gòu)，那么只需要定義一半濕面單元即可。

對(duì)帶支腿浮式結(jié)構(gòu)5種工況下的濕面進(jìn)行了網(wǎng)格劃分，如圖3所示。將靜水面作為濕面元?jiǎng)澐值钠鹚泓c(diǎn)，由于樁腿的升降而帶來(lái)的水線面的改變?cè)诖吮缓雎浴Ｒ虼?，平衡坐?biāo)系中零刻度線以下的面積均被認(rèn)為是浮式結(jié)構(gòu)的濕面。結(jié)構(gòu)主體平臺(tái)的型深為8 m，定吃水為4 m，為型深的一半，與上部分面積形成對(duì)稱。

5種工況中，帶支腿浮式結(jié)構(gòu)的主體船身上濕面單元的劃分是相同的，共有面元540個(gè)，各種工況之間，濕面的主要差距體現(xiàn)在樁腿上。工況1有648個(gè)面元。以此類推，工況2時(shí)結(jié)構(gòu)上共劃分了744個(gè)面元；工況3為840個(gè)面元；工況4為936個(gè)面元；工況5為1 032個(gè)面元。圖3給出了工況1和工況5的面元?jiǎng)澐质疽鈭D。

圖3 各種工況濕面元?jiǎng)澐諪ig.3 The set of the wetted panels

3.2水動(dòng)力系數(shù)分析

利用三維線性水彈性理論，通過(guò)數(shù)值方法，可求得帶支腿浮式全梁式結(jié)構(gòu)在波浪中運(yùn)動(dòng)的附加質(zhì)量和附加阻尼，即水動(dòng)力系數(shù)。這里重點(diǎn)討論二節(jié)點(diǎn)垂彎、一節(jié)點(diǎn)扭振和三節(jié)點(diǎn)垂彎這幾個(gè)模態(tài)的濕固有頻率，因此，圖4至圖6給出了這三個(gè)運(yùn)動(dòng)模態(tài)的水動(dòng)力系數(shù)對(duì)角線元素隨無(wú)因次頻率的關(guān)系曲線。圖中橫坐標(biāo)與縱坐標(biāo)所代表的物量均進(jìn)行了無(wú)量綱化處理。A為附加質(zhì)量，B為附加阻尼，ρ為流體密度，L為結(jié)構(gòu)長(zhǎng)度，Vol為浮式結(jié)構(gòu)排水體積，g為重力加速度，ω為圓頻率。

對(duì)于二節(jié)點(diǎn)垂彎來(lái)說(shuō)，如圖4(a)，工況3所代表的曲線高于其余幾種工況，而其余工況的水動(dòng)力系數(shù)也沒(méi)有按照水下樁腿的長(zhǎng)短進(jìn)行排列。這和剛體模態(tài)的特性表現(xiàn)出了差異[15]。剛體模態(tài)的分析中，通常是工況5所獲得的附加質(zhì)量最大，附加質(zhì)量與水下樁腿長(zhǎng)度有很好的正相關(guān)性。而此時(shí)，樁腿最長(zhǎng)的工況5所獲得的附加質(zhì)量還略低于水下樁腿長(zhǎng)度遠(yuǎn)低于它的工況4。由圖4(b)也可知，工況3的附加阻尼系數(shù)峰值是所有工況中最大。出現(xiàn)以上現(xiàn)象的原因在于，盡管工況3在五種工況中可積分的濕表面不是最大的，但經(jīng)過(guò)對(duì)比可知，樁腿處于工況3時(shí)，帶支腿浮式結(jié)構(gòu)主體平臺(tái)的歸一化振型的振幅最大，樁腿的濕表面積和主體平臺(tái)濕表面相比較是一個(gè)小量。因此，水動(dòng)力系數(shù)計(jì)算公式中，主體平臺(tái)的模態(tài)振型這一項(xiàng)就成為了影響水動(dòng)力系數(shù)大小的主要因素，而樁腿長(zhǎng)度變化導(dǎo)致的濕表面稱為次要因素，因此二節(jié)點(diǎn)垂彎模態(tài)中工況3的附加質(zhì)量和附加阻尼系數(shù)最大。而工況1和工況5主體平臺(tái)歸一化振型的振幅最小，因此，在二節(jié)點(diǎn)垂彎水動(dòng)力系數(shù)的計(jì)算中，這一組值相對(duì)來(lái)說(shuō)偏低。

圖4 二節(jié)點(diǎn)垂彎水動(dòng)力系數(shù)Fig.4 Hydrodynamic coefficients of 2-node VB

圖5 一節(jié)點(diǎn)扭轉(zhuǎn)振動(dòng)水動(dòng)力系數(shù)Fig.5 Hydrodynamic coefficients of 1-node tortion

圖6 三節(jié)點(diǎn)垂彎水動(dòng)力系數(shù)Fig.6 Hydrodynamic coefficients of 3-node VB

同時(shí)，盡管工況1和工況5主體平臺(tái)的模態(tài)振型相一致，但對(duì)比工況1和工況5的附加質(zhì)量發(fā)現(xiàn)，前者要低于后者。原因在于，工況1和工況5主體平臺(tái)的模態(tài)振型相同后，積分面積就成了水動(dòng)力系數(shù)積分當(dāng)中的主要因素，工況1水下的樁腿短于工況5，附加質(zhì)量系數(shù)自然略低于工況5的情況。因此當(dāng)振型一致時(shí)，樁腿水下長(zhǎng)短又成了影響水動(dòng)力系數(shù)的主要因素。工況2和工況4的附加質(zhì)量關(guān)系同樣可以用上述原理進(jìn)行解釋。一節(jié)點(diǎn)扭轉(zhuǎn)振動(dòng)和三節(jié)點(diǎn)垂彎振動(dòng)也有相類似的規(guī)律。

從上述分析可以看出，工況不同時(shí)，主體平臺(tái)的模態(tài)振型是主要因素，當(dāng)主體平臺(tái)振型一致時(shí)，樁腿在水下的長(zhǎng)度上升為影響水動(dòng)力系數(shù)的主要因素。這體現(xiàn)出對(duì)稱分析法在帶支腿浮式結(jié)構(gòu)分析中的優(yōu)勢(shì)，即在對(duì)稱與非對(duì)稱，相等與不相等中分清帶支腿浮式結(jié)構(gòu)水彈性分析中的主要因素和次要因素。按同樣的方法對(duì)下面的模態(tài)進(jìn)行分析。

3.3濕固有頻率分析

前面已經(jīng)分析了帶支腿浮式結(jié)構(gòu)各個(gè)工況的附加質(zhì)量、靜水恢復(fù)力和剛度矩陣，按式(6)可以計(jì)算出結(jié)構(gòu)不同工況下各階模態(tài)的濕固有頻率，表3給出了濕固有頻率的計(jì)算結(jié)果。

表3 各工況濕固有頻率 Tab.3 The wet resonant frequencies of the cases (rad/s)

對(duì)于帶支腿浮式結(jié)構(gòu)的對(duì)于彈性體模態(tài)，同一種振動(dòng)模態(tài)，樁腿長(zhǎng)度不同，濕固有頻率也不一樣。這三種振動(dòng)的濕固有頻率隨水下樁腿長(zhǎng)度的增加均呈現(xiàn)出減小的趨勢(shì)，如表3所示，當(dāng)帶支腿浮式結(jié)構(gòu)處于工況1時(shí)，其二節(jié)點(diǎn)垂彎濕固有頻率為5.507 0 rad/s；工況2時(shí)為5.406 9 rad/s；工況3時(shí)為5.349 5 rad/s；工況4時(shí)為5.047 8 rad/s；工況5時(shí)為4.597 6 rad/s。

這些數(shù)據(jù)都表明，和帶支腿浮式結(jié)構(gòu)干模態(tài)固有頻率一樣，樁腿同樣對(duì)帶支腿浮式結(jié)構(gòu)的濕固有頻率具有明顯地調(diào)制作用，但這種調(diào)制作用的表現(xiàn)形式和調(diào)制的機(jī)理卻各不相同。對(duì)于干模態(tài)固有頻率來(lái)說(shuō)，樁腿關(guān)于船身越趨于對(duì)稱分布時(shí)結(jié)構(gòu)的振動(dòng)固有頻率越高，樁腿的調(diào)整使得結(jié)構(gòu)的質(zhì)量分布發(fā)生了變化，使得固有頻率發(fā)生變化；對(duì)于濕固有頻率來(lái)說(shuō)，水下樁腿長(zhǎng)度地增加引起結(jié)構(gòu)附加質(zhì)量的非線性調(diào)整。從綜合效果看，有迫使共振頻率趨向低頻方向移動(dòng)趨勢(shì)。

4 結(jié) 語(yǔ)

采用有限元法和線性三維水彈性理論對(duì)一理想化的帶支腿浮式結(jié)構(gòu)進(jìn)行分析，計(jì)算了該浮式結(jié)構(gòu)水下樁腿處于不同長(zhǎng)度時(shí)各種模態(tài)的固有頻率，重點(diǎn)分析了引起結(jié)構(gòu)固有頻率變化的機(jī)理和原因。

前人的研究結(jié)果認(rèn)為對(duì)于帶支腿浮式結(jié)構(gòu)剛體模態(tài)來(lái)說(shuō)，水動(dòng)力系數(shù)大小主要取決于樁腿水下長(zhǎng)度，其水下樁腿長(zhǎng)度越大，水動(dòng)力相對(duì)系數(shù)越大。本文研究表明，對(duì)于彈性模態(tài)來(lái)說(shuō)，主體平臺(tái)振幅、水下樁腿長(zhǎng)度和水動(dòng)力系數(shù)三者之間成為一種耦合關(guān)系。結(jié)構(gòu)主體平臺(tái)振幅是決定水動(dòng)力系數(shù)大小的主要因素，樁腿處于對(duì)稱分布時(shí)，主體平臺(tái)的振幅最大，水動(dòng)力系數(shù)相對(duì)最大，同一種振型下，工況3的附加質(zhì)量可以是工況1的2～3倍。樁腿長(zhǎng)度成為了次要因素，只有當(dāng)振幅相同時(shí)，水下樁腿長(zhǎng)度越長(zhǎng)，水動(dòng)力系數(shù)才越大，此時(shí)，工況5的附加質(zhì)量可約為工況1的2倍，但同時(shí)振幅的大小又受到樁腿長(zhǎng)度的制約。

對(duì)于干模態(tài)固有頻率來(lái)說(shuō)，樁腿關(guān)于船身越趨于對(duì)稱分布，結(jié)構(gòu)的共振頻率向高頻移動(dòng)。對(duì)于濕固有頻率來(lái)說(shuō)，水下樁腿長(zhǎng)度的增加引起結(jié)構(gòu)附加質(zhì)量的非線性調(diào)整，從綜合效果看，則有迫使共振頻率趨向低頻方向移動(dòng)趨勢(shì)。由于樁腿的存在，帶支撐腿浮式結(jié)構(gòu)形成了具有一定寬度的共振頻域區(qū)間。二節(jié)點(diǎn)垂彎振動(dòng)的共振頻域區(qū)間相對(duì)較窄，而一節(jié)點(diǎn)扭轉(zhuǎn)振動(dòng)和三節(jié)點(diǎn)垂彎的共振頻域區(qū)間則相對(duì)較寬。

[1] HOWARTH M,DIER A,JONES W,et al.Jack-up response to wave-in-deck loads during extreme storms [J].Marine Structures,2004,17(3-4):275-290.

[2] TAN X M,LI J,LU C.Structural behaviour prediction for jack-up units during jacking operations [J].Computers & Structures,2003,81(24-25):2409-2416.

[3] CASSIDY M J,TAYLOR P H,EATOCK T R,et al.Evaluation of long-term extreme response statistics of jack-up platforms [J].Ocean Engineering,2002,29(13):1603-1631.

[4] JENSEN J J,CAPUL J.Extreme response predictions for jack-up units in second order stochastic waves by FORM [J].Probabilistic Engineering Mechanics,2006,21(4):330-337.

[5] BIENEN B,QIU G,TIM P.CPT correlation developed from numerical analysis to predict jack-up foundation penetration into sand overlying clay [J].Ocean Engineering,2015,(108):216-226.

[6] KOHAN O,BIENEN B,GUADIN C,et al.The effect of water jetting on spudcan extraction from deep embedment in soft clay [J].Ocean Engineering,2015,(97):90-99.

[7] 張浦陽(yáng),于曉洋,丁紅巖.海上自升式鉆井平臺(tái)插樁階段樁靴承載力計(jì)算[J].石油勘探與開發(fā),2011,38(5):613-619.(ZHANG P Y,YU X Y,DING H Y.Spudcan bearing capacity of the offshore jack-up drilling platform during the preloading process [J].Petroleum Exploration and Development,2011,38(5):613-619.(in Chinese))

[8] 鄭喜耀,周松望,姚首龍.自升式鉆井平臺(tái)插樁分析的幾個(gè)關(guān)鍵問(wèn)題 [J].海洋通報(bào),2014,33(1):84-89.(ZHENG X Y,ZHOU S W,YAO S L.Key issue for jack-up rig leg penetration analysis [J].Marine Science Bulletin,2014,33(1):84-89.(in Chinese))

[9] ZHANG A X,DUAN M L,LI H M.Prediction methods of spudcan penetration for jack-up units [J].China Ocean Engineering,2012,26(4):591-602.

[10] LI H T,YANG Q X,LI Y.Study on general performance of jack-up under elevated condition [J].China Ocean Engineering,2009,23(3):577-584.

[11] 殷玉梅,趙德有.船舶上層建筑整體縱向振動(dòng)固有頻率預(yù)報(bào)方法研究[J].船舶力學(xué),2011,15( 5):538-544.(YIN Y M,ZHAO D Y.Study on predicting method of natural frequency for superstructure's overall longitudinal vibration [J].Journal of Ship Mechanics,2011,15( 5):538-544.(in Chinese))

[12] 金鴻章,張宏瀚,賁成華.船舶雙水艙減搖系統(tǒng)設(shè)計(jì)與仿真研究[J].哈爾濱工程大學(xué)學(xué)報(bào),2008,29(3):242-250.(JIN H Z,ZHANG H H,BEN C H.An improve dual anti-rolling system for ships [J].Journal of Harbin Engineering University,2008,29(3):242-250.(in Chinese))

[13] 劉文璽,周其斗,張緯康,等.浮體與系泊系統(tǒng)的耦合動(dòng)力分析[J].船舶力學(xué),2014,18(8):940-957.(LIU W X,ZHOU Q D,ZHANG W K,et al.Dynamic analysis of the coupled floating body/mooring system [J].Journal of Ship Mechanics,2014,18(8):940-957.(in Chinese))

[14] 王顯正,趙德有,孫超,等.船舶總振動(dòng)固有頻率使用算法[J].中國(guó)艦船研究,2007,2(1):56-58.(WANG X Z,ZHAO D Y,SUN C,et al.Improved algorithm for the natural frequencies of ship vibration [J].Chinese Journal of Ship Research,2007,2(1):56-58.(in Chinese))

[15] 陳徐均,武海浪,黃亞新.水下樁腿長(zhǎng)度對(duì)帶支腿浮式結(jié)構(gòu)水動(dòng)力系數(shù)的影響分析[J].水動(dòng)力學(xué)研究與進(jìn)展,2014,29(1):83-92.(CHEN X J,WU H L,HUANG Y X.The effect of the underwater leg length on hydrodynamic coefficients of the floating structure with legs [J].Chinese Journal of Hydrodynamics,2014,29(1):83-92.(in Chinese))

[16] 黃亞新,武海浪,陳徐均.帶支腿浮式結(jié)構(gòu)水動(dòng)力建模及波浪力分析[J].哈爾濱工程大學(xué)學(xué)報(bào),2014,35(9):1067-1075.(HUANG Y X,WU H L,CHEN X J.The analysis of the modeling and wave excitation force on the floating structure with legs [J].Journal of Harbin Engineering University,2014,35(9):1067-1075.(in Chinese))

[17] WU H L,CHEN X J,HUANG Y X,et al.Influence of the legs underwater on the hydrodynamic response of the multi-leg floating structure[J].Ships and Offshore Structures，2014,9(6):578-595.

[18] 邱吉寶,向樹紅,張正平.計(jì)算結(jié)構(gòu)動(dòng)力學(xué)[M].合肥：中國(guó)科學(xué)技術(shù)大學(xué)出版社,2009:137-139.(QIU J B,XIANG S H,ZHANG Z P.Computational structure dynamics [M].Hefei:Press of USTC,2009:137-139.(in Chinese))

[19] WU Y S.Hydroelasticity of floating boddies [D].Ph.D.Thesis,Brunel Univ.,U.K.,1984.

[20] DANIELE D,ELENA C.Slamming clustering on fast ships:From impact dynamics to global response analysis [J].Ocean Engineering,2013(62):110-122.

[21] CHEN X J,WU Y S,CUI W C,et al.Nonlinear Hydroelastic analysis of a moored floating body [J].Ocean Engineering,2003(30):965-1003.

[22] 陳占陽(yáng)，任慧龍，李輝.水彈性理論與分段模型試驗(yàn)在船體振動(dòng)響應(yīng)分析中的應(yīng)用 [J].振動(dòng)與沖擊,2012,31 (24):119-124.(CHEN Z Y,REN H L,LI H.Application of hydroelasticity theory and segmented model test in hull vibration response analysis [J].Journal of Vibration and Shock,2012,31 (24):119-124.(in Chinese))

[23] CHEN Z Y,REN H L,LIU H,et al.Experimental and numerical analysis of bow slamming and whipping in different sea states [J].Journal of Ship Mechanics,2012,16(3):246-252.

[24] SENJANOVIC I,MALENICA S,TOMASEVIC S.Investigation of ship hydroelasticity [J].Ocean Engineering,2008,36(5-6):523-535.

[25] SENJANOVIC I,MALENICA S,TOMASEVIC S.Hydroelasticity of large container ships [J].Marine Structures,2009,22(2):287-314.

[26] SENJANOVIC I,TOMIC M,TOMASEVIC S.An explicit formulation for restoring stiffness and its performance in ship hydroelasticity [J].Ocean Engineering,2008,35(13):1322-1338.

The modulation effects of the legs on the natural frequencies of the floating structure

WU Hailang,CHEN Xujun,HUANG Yaxin,LIU Junyi

(College of Field Engineering,PLA University of Science and Technology,Nanjing 210007,China)

P751; U661.1

A

10.16483/j.issn.1005-9865.2016.03.006

1005-9865(2016)03-0047-08

2015-10-16

國(guó)家自然科學(xué)基金(51379213)；國(guó)家工信部高性能船舶科研項(xiàng)目(工信部聯(lián)裝[2012]533)；國(guó)家科技支撐計(jì)劃課題(2014bab16b05)

武海浪(1983-)，男，博士研究生，主要從事浮式工程結(jié)構(gòu)流固耦合動(dòng)力研究。E-mail：whlwzgrape@163.com

陳徐均。E-mail:chenxujun213@sina.com