梁恩強(qiáng)， 劉 昆， 包 杰， 王自力

(江蘇科技大學(xué)， 江蘇 鎮(zhèn)江 212003)

板-梁耦合技術(shù)在自升式海洋平臺(tái)碰撞分析中的應(yīng)用研究

梁恩強(qiáng)， 劉 昆， 包 杰， 王自力

(江蘇科技大學(xué)， 江蘇 鎮(zhèn)江 212003)

該文基于非線性有限元軟件ABAQUS，針對(duì)自升式海洋平臺(tái)結(jié)構(gòu)特點(diǎn)，提出了“板-梁”耦合技術(shù)，用以簡(jiǎn)化建立自升式海洋平臺(tái)的結(jié)構(gòu)模型。以平臺(tái)典型的單K和雙K節(jié)點(diǎn)為研究對(duì)象，基于“板-梁”耦合模型及板殼模型分析了其碰撞性能，同時(shí)比較驗(yàn)證了“板-梁”耦合技術(shù)的可靠性。研究結(jié)果表明，“板-梁”耦合模型在碰撞沖擊載荷下的結(jié)構(gòu)損傷變形、碰撞力和能量吸收等動(dòng)態(tài)響應(yīng)與板殼模型結(jié)果一致性較好。該文研究成果可以為大型平臺(tái)結(jié)構(gòu)整體碰撞性能分析提供技術(shù)支撐。

自升式海洋平臺(tái)；典型節(jié)點(diǎn)；“板-梁”耦合；碰撞性能；數(shù)值仿真

0 引言

隨著世界各國(guó)對(duì)海洋資源開(kāi)發(fā)的不斷重視，海洋平臺(tái)的數(shù)量持續(xù)增加，同時(shí)，海上運(yùn)輸業(yè)的高速發(fā)展使得船舶與海洋平臺(tái)碰撞的事故越發(fā)頻繁。以挪威大陸架為例，1982年～2010年間，共發(fā)生了115起船舶與海洋平臺(tái)的碰撞事故[1]，造成了嚴(yán)重的環(huán)境污染，甚至人員傷亡。這些事故說(shuō)明，在海洋平臺(tái)的設(shè)計(jì)研發(fā)過(guò)程中應(yīng)該考慮船舶撞擊等意外載荷的影響。同時(shí)，為了降低碰撞事故對(duì)人員及海洋環(huán)境的危害，開(kāi)展海洋平臺(tái)碰撞性能的研究，揭示平臺(tái)結(jié)構(gòu)在碰撞過(guò)程中的損傷變形機(jī)理，對(duì)更好的開(kāi)展平臺(tái)耐撞結(jié)構(gòu)設(shè)計(jì)有重要的現(xiàn)實(shí)意義。

隨著非線性有限元技術(shù)的發(fā)展和計(jì)算機(jī)硬件水平的提高，數(shù)值仿真已經(jīng)成為國(guó)內(nèi)外學(xué)者研究船舶與海洋工程碰撞問(wèn)題的主要方法。楊亮[2]基于有限元軟件ANSYS建立了整體全殼導(dǎo)管架平臺(tái)，對(duì)商船與固定導(dǎo)管架平臺(tái)的碰撞過(guò)程進(jìn)行仿真；孫勇敢[3]基于有限元軟件ANSYS建立了全殼半潛平臺(tái)模型，分析了平臺(tái)受到大型商船撞擊的動(dòng)態(tài)響應(yīng)；Martin Storheim等[4]利用有限元軟件LS.DYNA開(kāi)展了船舶與完整的全殼小規(guī)模平臺(tái)的碰撞數(shù)值仿真，研究平臺(tái)的耐撞性能；Joao Travanca等[5]基于有限元軟件Patran建立了多種小規(guī)模全殼平臺(tái)，預(yù)測(cè)了平臺(tái)受到典型補(bǔ)給船高速撞擊的動(dòng)態(tài)響應(yīng)。在研究小規(guī)模平臺(tái)碰撞問(wèn)題時(shí)，建立完整的全殼平臺(tái)模型是可行的，但是在研究大型平臺(tái)的碰撞響應(yīng)時(shí)，由于其結(jié)構(gòu)龐大且復(fù)雜，建立完整的全殼模型，不僅給建模帶來(lái)困難，而且在計(jì)算精度上也無(wú)優(yōu)勢(shì)，因此，平臺(tái)碰撞模型的合理簡(jiǎn)化成為高效、真實(shí)開(kāi)展數(shù)值仿真的關(guān)鍵。該文以一自升式海洋平臺(tái)為研究對(duì)象，針對(duì)其結(jié)構(gòu)特點(diǎn)提出“板-梁”耦合建模技術(shù)，分析其遭遇3 000 t補(bǔ)給船撞擊下的動(dòng)態(tài)響應(yīng)，通過(guò)與全殼模型計(jì)算結(jié)果的比較，驗(yàn)證“板-梁”耦合技術(shù)的可靠性。

1 “板-梁”耦合技術(shù)

“板-梁”耦合技術(shù)是在有限元建模過(guò)程中，將板殼單元與梁?jiǎn)卧ㄟ^(guò)多點(diǎn)耦合約束連接在一起，相互作用，相互影響。由于自升式海洋平臺(tái)結(jié)構(gòu)龐大，進(jìn)行全殼模型的建立并開(kāi)展碰撞仿真將耗費(fèi)大量計(jì)算資源，且由于碰撞具有明顯的局部性特點(diǎn)[6]，也沒(méi)有必要建立全殼模型，合理地簡(jiǎn)化模型可以高效、準(zhǔn)確地開(kāi)展碰撞仿真分析。該文針對(duì)船舶-平臺(tái)碰撞問(wèn)題，根據(jù)自升式平臺(tái)結(jié)構(gòu)特點(diǎn)，提出了“板-梁”耦合的平臺(tái)建模技術(shù)，即僅在碰撞區(qū)域內(nèi)的弦桿和撐桿結(jié)構(gòu)采用板殼單元，碰撞區(qū)域以外的結(jié)構(gòu)采用梁?jiǎn)卧獊?lái)模擬，梁元與板殼元通過(guò)多點(diǎn)耦合約束連接在一起。

2 碰撞有限元模型

2.1 自升式海洋平臺(tái)簡(jiǎn)介

計(jì)算選用某三樁腿桁架式自升式鉆井平臺(tái)，如圖1所示，其最大工作水深120 m，最大鉆井深度1 200 m，樁腿高167 m，弦桿中心距離13.1 m，樁腿節(jié)距8.5 m。

2.2 幾何模型

以自升式海洋平臺(tái)樁腿碰撞區(qū)域典型節(jié)點(diǎn)為對(duì)象，進(jìn)行數(shù)值仿真研究。樁腿節(jié)點(diǎn)構(gòu)件尺寸見(jiàn)表1，圖2為碰撞區(qū)整體的“板-梁”耦合結(jié)構(gòu)，圖3為平臺(tái)典型節(jié)點(diǎn)形式及幾何尺寸。

樁腿構(gòu)件外徑(mm)厚度(mm)斜撐245.019.0水平撐273.025.5弦管500.080.0

圖3 典型節(jié)點(diǎn)示意圖及尺寸參數(shù)

2.3 有限元模型

根據(jù)圖3所確定的兩種平臺(tái)典型節(jié)點(diǎn)，分別采用“板-梁”耦合技術(shù)和全殼單元建立有限元模型，如圖4、圖5所示。其中，“板-梁”耦合模型弦管在中部碰撞區(qū)域5 m(弦桿總長(zhǎng)8.5 m)范圍內(nèi)采用殼單元，斜撐在與弦管連接的一半長(zhǎng)度(3.8 m)范圍采用殼單元，其余區(qū)域采用梁?jiǎn)卧M，殼單元與梁?jiǎn)卧苯油ㄟ^(guò)建立耦合點(diǎn)進(jìn)行連接。管節(jié)點(diǎn)端部(弦桿和斜撐)剛性固定，撞擊船選取某3 000 t補(bǔ)給船。為了簡(jiǎn)化問(wèn)題，將其視為一直徑為2 m的剛性球，以2 m/s的初速度水平撞擊。管節(jié)點(diǎn)采用彈塑性材料，密度為7 850 kg/m3，彈性模量為2.1×1011N/m2，泊松比為0.3，屈服應(yīng)力為2.35×108N/m3，最大塑性失效應(yīng)變?nèi)?.3[5]。由于平臺(tái)用低碳鋼的塑性性能對(duì)應(yīng)變率是高度敏感的[8]，因此在材料模型中選用考慮應(yīng)變率敏感性影響的Cowper-Symonds本構(gòu)方程，其參數(shù)和分別取40.4和5[9,10]。

圖4 單K節(jié)點(diǎn)碰撞有限元模型

圖5 雙K節(jié)點(diǎn)碰撞有限元模型

3 結(jié)果對(duì)比分析

將以上模型提交計(jì)算，得到自升式海洋平臺(tái)節(jié)點(diǎn)的動(dòng)態(tài)響應(yīng)過(guò)程，從損傷變形、能量吸收和碰撞力三個(gè)方面分析平臺(tái)的碰撞性能并通過(guò)對(duì)比分析，驗(yàn)證“板-梁”耦合技術(shù)的合理性。

3.1 損傷變形

圖6、圖7分別給出了自升式海洋平臺(tái)典型節(jié)點(diǎn)在兩種不同建模方案下的損傷變形情況。從圖中可以看出：兩典型節(jié)點(diǎn)在不同建模方式下，結(jié)構(gòu)損傷變形情況基本一致，弦桿和斜撐均發(fā)生了較為明顯的彎曲變形，在碰撞接觸區(qū)域發(fā)生較大的塑性變形，且最大應(yīng)力和等效塑性應(yīng)變均出現(xiàn)在斜撐與弦管連接的位置。單K節(jié)點(diǎn)弦桿和斜撐的彎曲變形比雙K節(jié)點(diǎn)嚴(yán)重，這是由于雙K節(jié)點(diǎn)斜撐數(shù)量多，使得其抵抗變形的能力較強(qiáng)。通過(guò)分析比較可以發(fā)現(xiàn)，采用“板-梁”耦合技術(shù)建立的平臺(tái)典型節(jié)點(diǎn)簡(jiǎn)化模型與全殼模型在碰撞載荷作用下的損傷變形情況具有較好的一致性。

圖6 單K節(jié)點(diǎn)損傷變形情況

圖7 雙K節(jié)點(diǎn)損傷變形情況

3.2 能量吸收

圖8、圖9給出兩種不同建模方式下單K型及雙K型平臺(tái)節(jié)點(diǎn)能量轉(zhuǎn)化及吸收情況。從圖中可以看出，不同建模方式下兩節(jié)點(diǎn)對(duì)應(yīng)的能量轉(zhuǎn)化曲線趨勢(shì)基本一致，撞擊動(dòng)能大部分轉(zhuǎn)化為節(jié)點(diǎn)結(jié)構(gòu)的塑性變形能，且弦桿為主要的吸能構(gòu)件。表2為不同建模方式計(jì)算得到的吸能結(jié)果，比較兩節(jié)點(diǎn)的吸能情況可以發(fā)現(xiàn)，單K節(jié)點(diǎn)總的吸能略高于雙K節(jié)點(diǎn)，這是由其損失變形程度所決定的，但是，雙K節(jié)點(diǎn)斜撐吸能所占的比重高于單K節(jié)點(diǎn)。此外，“板-梁”耦合模型總的吸能略高于其所對(duì)應(yīng)的全殼模型，但是最多也僅高出1%左右，這仍然充分說(shuō)明了采用“板-梁”耦合技術(shù)建立的平臺(tái)典型節(jié)點(diǎn)簡(jiǎn)化模型與全殼模型在碰撞載荷作用下的能量吸收情況基本一致。

圖8 單K節(jié)點(diǎn)能量吸收情況

圖9 雙K節(jié)點(diǎn)能量吸收情況

模型結(jié)構(gòu)構(gòu)件吸能比例總吸能單K節(jié)點(diǎn)模型板殼模型“板梁”耦合模型弦桿4.426MJ78.1%斜撐1.238MJ21.9%弦桿4.625MJ80.8%斜撐1.100MJ19.2%5.664MJ5.725MJ雙K節(jié)點(diǎn)模型板殼模型“板梁”耦合模型弦桿3.975MJ70.5%斜撐1.663MJ29.5%弦桿3.950MJ69.4%斜撐1.742MJ30.6%5.638MJ5.692MJ

3.3 碰撞力

碰撞力是碰撞過(guò)程中結(jié)構(gòu)發(fā)生大位移、大變形的原因，通過(guò)碰撞力時(shí)程曲線可以清楚了解碰撞進(jìn)行的過(guò)程。圖10、圖11分別給出了單K節(jié)點(diǎn)和雙K節(jié)點(diǎn)的碰撞力-撞深曲線。從圖中可以看出， 碰撞力曲線具有較強(qiáng)的非線性特點(diǎn)，雙K節(jié)點(diǎn)的碰撞力震蕩更為明顯。另外，雙K節(jié)點(diǎn)的第一個(gè)碰撞力峰值較單K節(jié)點(diǎn)要高且出現(xiàn)時(shí)間略晚，這主要是由于其有四根斜撐支撐，其相應(yīng)承載能力及加載時(shí)間較兩根斜撐的單K節(jié)點(diǎn)要高。在第一個(gè)峰值之后，碰撞力出現(xiàn)一定卸載后又呈現(xiàn)上升趨勢(shì)，直至碰撞結(jié)束。通過(guò)觀察對(duì)應(yīng)時(shí)刻的損傷變形情況可以發(fā)現(xiàn)，在結(jié)構(gòu)屈服小幅卸載之后，碰撞接觸區(qū)域逐漸增加，響應(yīng)結(jié)構(gòu)也逐漸開(kāi)始由彈性變形向塑性變形過(guò)渡，隨著接觸面積的逐漸增大，碰撞力也隨之增加。在碰撞結(jié)束時(shí)，撞擊動(dòng)能基本耗盡，碰撞力迅速下降，被撞結(jié)構(gòu)凹陷變形的線性彈性部分逐漸恢復(fù)。此外，分別對(duì)應(yīng)比較兩種不同建模方式下的碰撞力曲線可以發(fā)現(xiàn)，不同建模方式對(duì)于兩種結(jié)構(gòu)的碰撞力影響不大，“板-梁”耦合模型的碰撞力稍滯后于全殼模型。

圖10 單K節(jié)點(diǎn)碰撞力-撞深曲線 圖11 雙K節(jié)點(diǎn)碰撞力-撞深曲線

4 結(jié)論

(1) “板-梁”耦合模型在碰撞載荷下的損傷變形、碰撞力和能量吸收等均與板殼模型結(jié)果具有較高的一致性，且減少了建模及計(jì)算工作量，適用于海洋平臺(tái)結(jié)構(gòu)的碰撞仿真分析。

(2) 兩典型平臺(tái)節(jié)點(diǎn)結(jié)構(gòu)在碰撞載荷作用下，弦管變形以彎曲為主而斜撐以失穩(wěn)為主，雙K節(jié)點(diǎn)抵抗變形的能力較單K節(jié)點(diǎn)強(qiáng)，但兩節(jié)總的塑形變能相當(dāng)，且弦桿吸收了較多的能量。

[ 1 ] Arne K. Collisions between platforms and ships in Norway in the period 2001-2010[C]. OMAE,2010.

[ 2 ] 楊亮. 海洋導(dǎo)管架平臺(tái)的耐撞特性研究[D]. 大連：大連理工大學(xué),2007.

[ 3 ] 孫勇敢. 海洋平臺(tái)碰撞性能研究[D]. 大連：大連理工大學(xué),2011.

[ 4 ] Martin S. Design of offshore structures against accidental ship collisions[J]. Marine Structures, 2014,37(6):135-172.

[ 5 ] Joao T, HONG Hao. Dynamics of steel offshore platforms under ship impact[J]. Applied Ocean Research,2014,47(5):352-372.

[ 6 ] 胡永利,林一. 半潛式平臺(tái)遭遇碰撞的結(jié)構(gòu)響應(yīng)分析[J]. 船舶與海洋工程,2011,28(1):46-54.

[ 7 ] Liu Kun, Wang Zili, Tang Wenyong, 等. Experimental and numerical analysis of laterally impacted stiffened plates considering the effect of strain rate[J]. Ocean Engineering，2015，99(3)：44-54.

[ 8 ] 王自力,顧永寧. 應(yīng)變率敏感性對(duì)船體結(jié)構(gòu)碰撞性能的影響[J]. 上海交通大學(xué)報(bào),2000,34(12):1704-1707.

[ 9 ] 王自力,蔣志勇,顧永寧. 船舶碰撞數(shù)值仿真的附加質(zhì)量模型[J]. 爆炸與沖擊,2002,22(4):321-326.

[10] 陶亮. 船舶舷側(cè)結(jié)構(gòu)碰撞性能研究[D]. 大連：大連理工大學(xué),2005.

[11] 胡永利,林一等. 半潛式平臺(tái)遭遇碰撞的結(jié)構(gòu)響應(yīng)分析[J]. 船舶與海洋工程,2011,28(1):46-54.

[12] 劉昆,張延昌,王自力. 船首形狀對(duì)舷側(cè)結(jié)構(gòu)碰撞性能影響研究[J]. 船舶工程,2010,32(2):12-14.

The Application of Plate-beam Coupling Technology in Offshore Jack-up Platform Collision Analysis

LIANG En-qiang, LIU Kun, BAO Jie, WANG Zi-li

(Jiangsu University of Science and Technology, Jiangsu Zhenjiang 212003, China)

The "plate-beam" coupling technology is proposed based on the nonlinear finite element software ABAQUS for building the structural model of Jack-up platform. The collision performance between the shell models and the "plate-beam" coupling models are compared， and the results indicate that the damage deformation, impact force and energy absorption of the "plate-beam" coupling model are in good agreement with the shell model. The results can provide technical support for the analysis of large-scale platform collisions.

jack-up offshore platform; typical joint; plate-beam coupling technology; collision performance; numerical simulation

2015-03-19

國(guó)家自然科學(xué)基金(51379093)。

梁恩強(qiáng)(1989-)，男，碩士研究生。

1001-4500(2015)05-0028-07

P75

A