錐體海洋平臺結構冰荷載的離散單元分析

李 曄

（中國船級社 海洋工程技術中心，天津300457）

錐體海洋平臺結構冰荷載的離散單元分析

李 曄

（中國船級社 海洋工程技術中心，天津300457）

在冰區油氣開發中，錐體結構可以有效降低冰力，避免強烈的冰激振動，是目前渤海油氣平臺的主要結構形式。為研究海冰與錐體結構的相互作用過程，文章建立了適用于模擬海冰破碎特性的離散單元模型。該模型將海冰離散為若干個具有粘接-破碎功能的顆粒單元，并通過海冰彎曲試驗確定了單元間的粘接強度；然后對海冰與錐體結構的作用過程進行了數值計算，獲得了相應的動冰荷載及冰振響應；在此基礎上討論了不同錐角影響下冰荷載及結構振動響應的變化規律。結果表明，水平方向冰荷載及結構冰振響應隨錐角的增加明顯增加，而豎直方向冰荷載則顯著降低。該離散單元模型還可進一步應用于不同類型抗冰結構的冰荷載分析，有助于解決冰區結構物的抗冰結構設計和冰致疲勞分析。

海冰；錐體結構；冰荷載；冰激振動；離散單元模型

0 引 言

海冰與海洋結構相互作用時產生的冰荷載主要受海冰性質與結構形式的影響。在冰區油氣開發的導管架式海洋平臺結構設計中，人們最初采用的是直立結構形式，但直立結構存在冰力偏大和冰激振動嚴重等問題[1-3]。由于海冰是一種典型的脆性材料，其彎曲強度要明顯低于壓縮強度。為此，人們提出了在直立結構的水面處安裝破冰錐體以降低冰荷載。海冰與斜面或錐體結構相互作用過程的研究一直是人們所關注的問題[4-6]。特別是近年來，渤海錐體結構的廣泛采用，錐體結構冰荷載及相應結構冰振響應引起了人們更加廣泛的關注[7-10]。

近年來，人們通過理論分析、試驗測試和數值方法研究了海冰與錐體結構的相互作用過程[10-12]。在海冰與錐體結構相互作用過程中，海冰以彎曲破壞為主，其中海冰彎曲強度是決定冰力的重要參數；此外，錐體結構冰荷載還與冰厚和冰速等海冰參數、錐角和錐徑等結構形式等有著密切的關系。由于海冰與錐體結構作用時的破壞過程非常復雜，影響因素較多，其冰荷載特征還需要進一步深入研究。

海冰荷載的數值模擬大多采用有限元方法[11,13]。該方法能夠很好地分析海冰與海洋結構相互作用中冰內應力分布以及海冰破碎特性，但尚不能合理地模擬海冰連續破碎過程中的動冰力。在有限元方法的基礎上，人們發展了一種粘接單元模型用以模擬海冰的破碎特征，并由此確定相應的冰荷載[14]。最近，基于離散元模型的海冰與結構物相互作用研究取得了很大的進展[15-18]。它能夠合理地模擬海冰與結構作用中的破碎過程，并獲得相應的冰力時程。離散元模型在處理海冰與錐體結構作用時，將冰排離散為若干個具有粘結—破碎功能的顆粒單元，通過顆粒材料性質及單元間粘接強度模擬海冰的破碎模式[16]。海洋平臺可簡化為具有等效的質量—剛度—阻尼的計算模型。由于顆粒單元模型具有物理意義明確、計算模型簡便、接觸判斷快速等優點，可較好地模擬海冰與海洋結構相互作用中的破碎特性，并由此確定冰荷載的動力特性以及海冰作用下的結構振動。

為此，本文采用粘接顆粒離散單元模型構造海冰材料，通過對海冰彎曲試驗的離散元分析確定單元間的粘結強度；在此基礎上對海冰與錐體之間的相互作用進行數值分析，確定相應的冰荷載及結構冰振響應；最后對錐角影響下的冰荷載和冰振響應進行討論。

1 海冰的離散單元模型及彎曲強度數值模擬

海冰在與錐體結構的相互作用過程中主要呈現彎曲破碎特性。這里首先建立海冰的離散單元模型，并通過海冰的彎曲試驗進行模型的檢驗和計算參數的確定。

1.1 海冰離散單元模型

在采用離散單元模型計算海冰的動力作用及破碎過程時，將海冰離散為若干具有一定質量和大小的顆粒單元，單元間具有相應的粘接作用。顆粒之間可視為在法向上由彈簧與阻尼器串聯連接，在切向上由彈簧、阻尼器和滑動摩擦器連接。在顆粒相互作用過程中，一般考慮單元間因相對速度和彈性變形而引起的粘彈性作用力，并采用Mohr-Coulomb摩擦定律計算剪切力。

海冰彎曲強度的試驗表明，海冰彎曲強度是海冰鹵水體積（溫度、鹽度）、加載速率等因素的函數[20-21]。由于海冰的彎曲強度主要由顆粒間的粘接強度決定，而粘結強度主要受鹵水體積的影響，可將海冰單元間的粘接強度用單元間最大粘接強度σmbax表示，即

式中：β（ vb）為鹵水體積vb影響下的海冰強度折減系數。海冰單元間的最大粘接強度σmbax可通過海冰彎曲強度的敏度分析進行確定?？紤]海冰的壓縮和彎曲強度與鹵水體積有相似的對應關系[21]，則有

式中：vb可設為海冰溫度和鹽度的函數[22]，即

式中：T為海冰溫度（℃），S為海冰鹽度（‰）。

1.2 海冰彎曲強度的離散單元模擬

采用以上海冰離散元模型對海冰三點彎曲強度試驗進行數值模擬，其主要計算參數列于表1中。海冰單元按一定排開方式組合后粘接為海冰彎曲試樣，如圖1（a）所示。試樣下側的兩個支撐點固定，其上側加載點以勻速0.01 m/s向下運動。計算得到的海冰彎曲斷裂過程及所對應的冰力時程曲線分別如圖1和圖2所示。圖中顆粒顏色代表其平均接觸力，由此可發現海冰受彎時中性層的位置及其上下兩側的對稱受力情況。從圖2中可以發現，應力隨著時間大體呈線性增長，并存在一定的波動性。當t=0.033 s時，應力達到最大值，此時海冰試樣斷裂，其對應的海冰彎曲強度σf為1.57 MPa。

表1 海冰彎曲強度的離散單元模擬主要計算參數Tab.1 Main computational para meters in numerical simulation of sea ice flexural strength with DEM

圖1 離散元數值模擬的海冰試樣彎曲破碎過程Fig.1 Bending failure process of sea ice sample simulated with DEM

圖2 離散元模擬的海冰彎曲過程中的應力變化時程Fig.2 The stress versus time in sea ice bending process simulated with DEM

圖3 不同粘接強度下離散元計算的海冰彎曲強度Fig.3 Sea ice flexural strengths simulated with DEM under various bonding strengths

圖4 不同鹵水體積下計算的海冰彎曲強度Fig.4 Sea ice flexural strengths simulated under various brine volumes

2 錐體結構冰荷載的離散元模擬

采用離散單元模型對海冰與錐體結構的相互作用過程進行數值分析，由此確定冰力時程及平臺結構的冰振響應，并在此基礎上進一步討論錐角對冰力的影響。

2.1 海冰與錐體結構的離散單元模擬

在海冰與錐體結構的相互作用過程中，海冰主要以彎曲破碎為主，顆粒間發生拉伸和剪切破壞。遼東灣現場監測的海冰與錐體結構的作用狀況如圖5所示。下面采用離散元模型對海冰與錐體結構作用時的破碎過程、冰荷載時程和冰激結構振動響應進行數值計算，主要計算參數列于表2中。

圖5 遼東灣JZ20-2油氣田海冰與錐體相互作用Fig.5 Conical offshore structure in the JZ20-2 oil field of Liaodong Bay and the interaction between ice cover and cone

采用離散元方法模擬的海冰與錐體結構相互作用過程如圖6所示，圖中顆粒顏色代表其平均作用力，也反應了冰內應力場分布的大致規律。從海冰破碎后的冰塊尺度來看，采用平行粘接模型可模擬出海冰與錐體結構作用時的破碎長度。計算得到的x向樁腿冰力時程如圖6（a）所示，其最大和平均冰力分別為62.68 kN和18.41 kN。該冰力時程呈現出明顯的周期性特征，這與海冰和錐體結構相互作

用時的彎曲破碎長度和冰速密切相關[7,10]。計算得到的y向冰力如圖6（b）所示，其最大和平均值分別為18.51 kN和0.67 kN。這表明海冰對錐體側面也會產生一定的沖擊力，但其作用力是基本對稱的。z方向的冰力時程如圖6（c）所示，其最大值和均值分別為49.09 kN和16.20 kN，其峰值也顯示出一定的周期性。錐體結構的冰振加速度時程如圖6（d）所示，其中最大加速度為9.25 gal（1.0 gal=1.0 cm/s2）。以上結果表明，該離散元模型可以較好地模擬海冰與錐體結構的相互作用過程，確定動冰荷載及錐體結構的冰振響應。

表2 海冰與錐體結構作用離散元模擬的主要計算參數Tab.2 Main computational parameters of DEM simulation of dynamic ice load on conical structures

圖6 離散元模擬的海冰與錐體結構作用過程Fig.6 The interaction between sea ice and the conical structure simulated with DEM

圖7 離散元模擬的錐體海洋結構冰荷載及冰速方向冰振加速度Fig.7 Dynamic ice loads and the ice-induced acceleration of conical offshore structure simulated with DEM

2.2 錐角對錐體結構冰荷載的影響

在海冰與錐體結構相互作用中，海冰以彎曲破壞為主，并由此產生相應的動冰力。海冰發生彎曲破壞的模式及其相應的破碎長度在很大程度上受錐體角度的影響。錐體角度越大，海冰受到的錐體作用力的水平分量也就越大，豎向分量則越小。為進一步分析錐角對冰荷載的影響，這里將錐角分別設為65.8°、55.0°、46.5°和39.0°，其他計算參數采用表1和表2中的數值。當海冰與錐體的作用時間為6.5 s時，計算得到的不同錐角下海冰破碎現象以及冰力時程分別如圖8和圖9所示?？梢钥闯?，海冰破碎的尺寸在不同錐角下有明顯差異。x方向的冰力時程更清晰地表明，錐角越大，冰力周期也就越小，即海冰的破碎尺寸也就越小。

圖8 不同錐角下海冰與錐體的相互作用Fig.8 Interaction between sea ice and the conical structure with different conical angles

圖9 離散元模擬的不同錐角下的水平和豎向冰力Fig.9 Ice loads in horizontal and vertical directions simulated with DEM under various conical angles

圖10 不同錐角下的冰力及結構冰振加速度Fig.10 Mean ice loads and accelerations of ice-induced structure vibration simulated under various conical angles

在不同錐角下，錐體結構的水平和豎向冰力均值、最大值如圖10（a）和（b）所示?？梢园l現水平方向冰力的均值和最大值隨錐角的增加而增大，而豎向冰力則均呈減小趨勢。這主要是由于在錐角由小到大變化過程中，海冰破碎由彎曲模式向擠壓模式轉變，并由此導致冰力水平分量的增加和豎向分量的降低。不同錐角下平臺結構的冰振加速度均值和最大值如圖10（c）所示，其變化規律與水平方向冰力是一致的。這主要是因為平臺結構的冰振響應主要受控于水平方向冰力。

3 結 語

為研究海冰與錐體海洋平臺結構的作用過程，本文采用離散單元模型對海冰的彎曲強度、錐體結構動冰力及冰振響應進行了數值計算。通過海冰彎曲強度的離散元數值模擬，分析了鹵水體積對海冰彎曲強度的影響，并由此確定了海冰單元間的粘結強度，驗證了離散元模擬的可行性；在此基礎上對海冰與錐體結構的作用過程進行了離散元計算，得到了動冰力及結構冰振響應；通過對不同錐角下結構冰力和冰振響應的對比，發現錐角越大水平冰力和冰振響應就越大，而豎向冰力則相應減少。

在海冰與錐體結構相互作用中，海冰的斷裂長度和正倒錐結合處的海冰破碎模式是兩個重要研究內容。這需要進一步對海冰離散元模型進行改進并確定合理的計算參數以更加全面地揭示錐體結構的冰力特性。此外，對海洋平臺結構冰荷載及冰振響應的離散元模擬還需要通過室內模型試驗和現場監測結果進行驗證，以提高其計算精度和可靠性。該海冰離散元模型對深入研究海冰與海洋結構的作用特性具有很好的應用前景，可對不同類型海洋平臺的冰荷載、冰激結構振動機理和多樁腿結構的冰力掩蔽效應進行分析。

致謝：感謝大連理工大學工業裝備結構分析國家重點實驗室季順迎教授的有益討論。

[1]Karna T,Qu Y,Kuhnlein W,Yue Q,Bi X.A spectral model of ice forces due to ice crushing[J].Journal of Offshore Mechanics&Arctic Engineering,2007,129:138-145.

[2]Huang Y,Shi Q,Song A.Model test study of the interaction between ice and compliant vertical narrow structure[J].Cold Regions Science and Technology,2007,49:151-160.

[3]Sodhi D.Crushing failure during ice-structure interaction[J].Engineering Fracture Mechanics,2001,68:1889-1921.

[4]季順迎,王安良,車嘯飛,岳前進.錐體導管架海洋平臺冰激結構振動響應分析[J].海洋工程,2011,29(2):32-39.Ji S Y,Wang A L,Che X F,et al.Analysis of ice induced structure vibration of offshore jacket platform with ice breaking cone[J].Ocean Engineering,2011,29(2):32-39.

[5]Huang Y,Li X.Dynamic ice loads on conical structures[J].Theoretical&Applied Mechanics Letters,2011,1:022007.

[6]Timco G.Ice loads on conical-shaped structure[J].Cold Regions Science and Technology,2009,55:1-2.

[7]岳前進,畢祥軍,于 曉,等.錐體結構的冰激振動與冰力函數[J].土木工程學報,2001,36(2):16-19.Yue Q J,Bi X J,Yu X,et al.Ice-induced vibration and ice force function of conical structure[J].China Civil Engineering Journal,2001,36(2):16-19.

[8]史慶增,黃 焱,宋 安,等.錐體冰力的實驗研究[J].海洋工程,2004,22(1):88-92.Shi Q Z,Huang Y,Song A,et al.Experimental study of ice force on conical structures[J].Ocean Engineering,2004,22(1):88-92.

[9]李 峰,岳前進.不同破壞模式下錐體結構冰力定性分析[J].大連理工大學學報,2005,45(6):785-789.Li F,Yue Q J.Qualitative analysis of ice forces on conical structures under different failure modes[J].Journal of Dalian University of Technology,2005,45(6):785-789.

[10]Xü N,Yue Q,Bi X,et al.Experimental study of dynamic conical ice force[J].Cold Regions Science and Technology,2015,120:21-29.

[11]王 剛,武文華,岳前進.錐體接觸寬度對冰排彎曲破壞模式影響的有限元分析[J].工程力學,2008,25(1):235-240.Wang G,Wu W H,Yue Q J.FEM analysis on ice-bending failure mode with width effect of ice-cone interaction[J].Engineering Mechanics,2008,25(1):235-240

[12]Huang Yan.Model test study of the nonsimultaneous failure of ice before wide conical structures[J].Cold Regions Science and Technology,2010,63:87-96.

[13]Konuk I,Gurtner A,Yu S.Study of dynamic ice and cylindrical structure interaction by the cohesive element method[C]//Proceedings of the 20th International Conference on Port and Ocean Engineering under Arctic Conditions.Lulea,Sweden,2009.

[14]Gagnon R E.A numerical model of ice crushing using a foam analogue[J].Cold Regions Science and Technology,2011,65:335-350.

[15]Lau M,Lawrence K P,Rothenburg L.Discrete element analysis of ice loads on ships and structures[J].Ships and Offshore Structures,2011,6(3):211-221.

[16]狄少丞,季順迎.海冰與自升式海洋平臺相互作用GPU離散元模擬[J].力學學報,2014,46(4):561-571.Di S C,Ji S Y.GPU-based discrete element modelling of interaction between sea ice and jack-up platform structure[J].Chinese Journal of Theoretical and Applied Mechanics,2014,46(4):561-571.

[17]Poloj?rvi A,Tuhkuri J,Pustogvar A.DEM simulations of direct shear box experiments of ice rubble:Force chains and peak loads[J].Cold Regions Science&Technology,2015,116:12-23.

[18]季順迎,李春花,劉 煜.海冰離散元模型的研究回顧及展望[J].極地研究,2012,24(4):315-329.Ji S Y,Li C H,Liu Y.A review of advances in sea ice discrete element models[J].Chinese Journal of Polar Research,2012,24(4):315-329.

[19]Ji S,Shen H H.Internal parameters and regime map for soft polydispersed granular materials[J].Journal of Rheology,2008,52(1):87-103.

[20]Timco G W,O’Brien S.Flexural strength equation for sea ice[J].Cold Regions Science and Technology,1994,22:285-298.

[21]王安良,許 寧,畢祥軍,等.鹵水體積和應力速率影響下海冰強度的統一表征[J].海洋學報,2016,38(9):126-133.Wang A L,Xü N,Bi X J,et al.Unified representation of sea ice strengths under influences of brine volume and stress rate[J].Acta Oceanolo Gica Sinica,2016,38(9):126-133.

[22]Frankenstein G,Garner R.Equations for determining the brine volume sea ice from-0.5℃ to-22.0℃[J].Journal of Glaciology,1967,6(48):943-944.

Analysis of ice load on conical structure with discrete element model

LI Ye

(Offshore Technology Center of China Classification Society,Tianjin 300457,China)

In oil/gas exploitations of ice covered cold regions,the conical structure can reduce ice force and avoid intense structure vibration.It is applied most widely in the Bohai Sea recently.In order to investigate the interaction between ice cover and conical structure,a discrete element method is established to model the sea ice breaking characteristics.The ice cover is dispersed into a series of bonded spherical elements.The bonding strength among ice particles is determined with the flexural test of sea ice.The interaction between the ice cover and the conical structure is simulated with the discrete element method to determine the dynamic ice load and structure vibration.Moreover,the influence of conical angle is discussed based on the simulated results.The ice load and the ice-induced structure response in horizontal direction increase obviously with the increase of cone angle,while the ice load in vertical direction decreases significantly.This discrete element method can also be applied to simulate ice loads of other type structures,and to be aided for the design of ice-resistance structure and the analysis of ice induced structure fatigue.

sea ice;conical structure;ice load;ice-induced vibration;discrete element method

P731.15 P751

A

10.3969/j.issn.1007-7294.2017.10.009

1007-7294（2017）10-1254-09

2017-03-27

國家重點研發計劃（2106YCF1401505）和國家自然科學基金項目（41576179，51639004）

李 曄（1975-），男，博士，高級工程師，E-mail:liye@ccs.org.cn。