冰載荷仿真及其在海上風機一體化分析中的應用

王凱， 劉天輝， 萬遠琛， 施偉

(1.中山大學 海洋工程與技術學院，廣東 珠海 519082； 2.南方海洋科學與工程廣東省實驗室(珠海)，廣東 珠海 519082； 3.大連理工大學 海岸和近海工程國家重點實驗室， 遼寧 大連 116024)

面對日益嚴峻的環(huán)境污染、溫室氣體排放和氣候變化等問題,許多國家都在調(diào)整能源發(fā)展戰(zhàn)略,煤電投資熱度持續(xù)降溫,加快可再生能源的開發(fā)利用已成為目前世界各國的普遍共識。新世紀以來，風能發(fā)電在世界范圍內(nèi)迅速發(fā)展，并在可再生能源發(fā)電領域中占據(jù)重要地位。風能發(fā)電可以按照其安裝的位置分為陸上風電和海上風電。相對陸上風電, 海上風電在風資源裕度、風電場位置、風力機裝機容量等方面具有優(yōu)勢, 在風力機制造安裝復雜度、全壽命周期成本、運行維護難度、對電網(wǎng)調(diào)峰裕度的要求等方面存在劣勢。綜合考慮后海上風電的利大于弊，具有很大開發(fā)價值，海上風電是世界眾多國家推動能源低碳轉(zhuǎn)型的核心路徑。

海上風力發(fā)電在高緯度寒區(qū)海域，例如歐洲波羅的海、亞洲渤海灣、美洲五大湖等區(qū)域有很大的發(fā)展?jié)摿?，對于寒區(qū)海上風機而言海冰載荷是重要載荷之一。風能探測的報告顯示我國渤海及黃海北部海域風能密度大，然而冬季會有海冰影響風電設備。大量工程項目表明, 海冰會對海洋結(jié)構產(chǎn)生巨大的危害。1969年渤海發(fā)生特大冰封, 巨大的海冰直接將海上的“海二井”石油平臺推倒, 造成了建國以來最大的由于海冰導致的石油平臺事故。放眼全球, 阿拉斯加庫克灣的采油平臺[1]、北歐的Bothnia灣燈塔[2]等冰區(qū)海洋結(jié)構都不同程度地遭到過海冰的破壞。因此可以說, 海冰是關乎海洋結(jié)構安全的全球性問題。

1 我國冰載荷的研究現(xiàn)狀

目前我國關于冰載荷的研究發(fā)展尚不成熟，而且大多集中在海洋平臺領域，海冰是高緯度海域海洋結(jié)構物設計和安全運營的重要影響因素，研究表明海冰載荷往往是寒區(qū)海域影響海洋工程結(jié)構物穩(wěn)定的控制載荷[3]。早在1966年，我國就開始在渤海海域建設抗冰的海上平臺，至今已成功建成了數(shù)十座具有抗冰能力的海上平臺。

冰載荷的研究只在少數(shù)的幾個重視寒區(qū)海洋能源開發(fā)的國家開展，包括美國、北歐、加拿大、俄羅斯、日本、中國等，相對于波浪載荷、風載荷、地震載荷等問題的研究，冰載荷問題的工程需要相比之下沒那么廣泛，而且各個國家海域的冰情、采用的結(jié)構形式等均有著較大的差異，研究重點也根據(jù)實際海域有所不同，造成彼此研究的相對獨立。我國的冰載荷研究主要圍繞渤海展開，目前的研究主要分布在以下幾個方面：

1)渤海的冰情調(diào)查，包括歷史冰情記錄[4-5]，冰情預報方法[6]統(tǒng)計分析劃分海冰區(qū)域以及給定工程設計參數(shù)[7];

2)海冰的物理力學性能[8](包括溫度、鹽度、密度等物理參數(shù)和壓縮強度、彎曲強度等力學參數(shù));

3)冰力及其檢測方法[9];

4)冰與結(jié)構的動力相互作用和結(jié)構的冰振響應[10];

5)結(jié)構的冰激疲勞分析以及抗冰設計[11]。

2 冰載荷數(shù)值仿真方法

2.1 數(shù)值仿真概況

針對冰載荷研究主要方式有現(xiàn)場實測、原型試驗、理論分析和數(shù)值模擬。對冰載荷的研究主要針對靜冰載荷和動冰載荷。由于海洋環(huán)境復雜，影響海冰運動的因素較多，因此如何建立合理有效的數(shù)值方法來模擬海冰動力學過程一直都是研究的重點，雖然海冰與結(jié)構物相互作用的過程至今沒有公認的數(shù)學物理模型，不過研究者們還是基于模型實驗和現(xiàn)場觀測為基礎，用數(shù)值模擬方法展現(xiàn)冰和結(jié)構相互作用的過程，特別是對海冰動態(tài)破壞的過程。

在分析冰與一般的結(jié)構物作用時，通常所研究的冰尺度較小。對于較小的尺度，若假設其為均勻介質(zhì)就可以很簡單地運用目前最為常用的數(shù)值分析方法有限元法[17-21]。除此以外還有離散元法[22-23]、非連續(xù)變形分析方法[15-16]等其他數(shù)值分析方法在處理不同問題時有著各自優(yōu)勢。

2.2 有限元法

有限元分析是解決工程和數(shù)學模型時最為廣泛使用的方法。將大型的系統(tǒng)細分為小而簡單的元，在滿足工程要求的前提下設置元素，將求解域劃分為離散域，再根據(jù)離散域各自的方程進行求解，就可以用有限數(shù)量的未知量代數(shù)方程組去求解偏微分方程組。對于不同數(shù)學模型的問題有限元的求解方法是大致相同的，其核心思想均為“數(shù)值近似”和“離散化”，求解問題的基本步驟[12]通常是根據(jù)問題確定好計算域，然后將計算域離散化，再確定狀態(tài)變量的方程組和設定邊界條件開始進行計算，計算得出結(jié)果后再進行分析。處理步驟大概可分成3個階段，前處理、求解和后處理流程如圖1所示。前處理就是根據(jù)研究的問題建立有限元模型并且根據(jù)計算需求完成單元網(wǎng)格劃分；后處理則是在得到計算結(jié)果以后根據(jù)計算的設定條件分析結(jié)果，以方便其他研究者了解。有限元方法是基于連續(xù)介質(zhì)理論用來描述連續(xù)結(jié)構體運動應用最廣泛的數(shù)值計算方法，適合用來分析海洋結(jié)構在冰力作用下動力響應。

2.3 離散元法

離散元法是專門用來解決不連續(xù)介質(zhì)問題的數(shù)值模擬方法最初往往應用于巖石力學領域，把介質(zhì)離散為元，相鄰的元之間存在某種或幾種作用力，元的運動滿足牛頓第二定律。近年來，基于離散元方法模擬海冰與船體和海洋平臺結(jié)構間的相互作用研究取得了較大進展。冰離散元方法的計算單元可采用球體、圓盤、多面體和擴展多面體等不同形態(tài)。采用離散單元模型計算更加適合模擬離散顆粒組合體在準靜態(tài)或是動態(tài)條件下的變形和破碎過程，而海冰是屬于非連續(xù)介質(zhì)可以離散為具有一定質(zhì)量的顆粒單元，顆粒之間允許位移重疊或分離，通過這些單元可以靈活構建不同類型的海冰，包括平整冰、浮冰、冰脊和堆積冰等。

圖1 有限元分析示意Fig.1 Finite element analysis diagram

2.4 其他數(shù)值方法

除了前述2種主流數(shù)值分析方法之外，還有許多數(shù)值方法用于冰載荷的研究，如差異元素法[13]、三維顆粒流法[13]、幾何網(wǎng)格法[14]等。石根華[15]于 1988 年提出了一種新型數(shù)值分析方法將數(shù)學拓撲學和工程實踐相結(jié)合的非連續(xù)變形分析(discontinuous deformation analysis，DDA)。該方法最初是為了解決巖體大變形和大位移問題時介質(zhì)不連續(xù)性導致計算誤差大的問題而提出的，平行于有限元法，充分考慮結(jié)構介質(zhì)的不連續(xù)性，將結(jié)構面切割而成的塊體作為分析單元，將動力學與靜力學統(tǒng)一起來，用最小勢能原理把塊體之間的接觸問題和塊體本身的變形問題統(tǒng)一到矩陣中求解，后來被學者運用到冰與結(jié)構的相互作用模擬中。張運良等[16]就基于DDA算法模擬冰與結(jié)構相互作用的過程，將冰破碎的過程可視化，計算了冰撞擊結(jié)構后每個時間步的結(jié)構動態(tài)響應。模擬出的冰破碎過程與原型觀測到的過程較為吻合，與有限元和邊界元方法相比，DDA方法對冰的破碎模擬更為接近實際。不過目前整體而言該方法在應用于冰載荷研究還相對較少。

3 數(shù)值仿真方法在抗冰研究中的應用

在研究冰與海洋相互作用時所關注的問題往往是作用在結(jié)構上的最大靜冰力和動冰力形式、海冰引起結(jié)構振動的機理以及結(jié)構在冰載荷作用下的失效形式，按照結(jié)構與冰作用部分的形狀可以分為直立結(jié)構基礎及非直立結(jié)構基礎。現(xiàn)在數(shù)值仿真在冰載荷方面的應用往往是通過與模型試驗數(shù)據(jù)和現(xiàn)場測試數(shù)據(jù)對比以驗證數(shù)值仿真方法的適用性，積累冰載荷領域數(shù)值仿真的經(jīng)驗，并在此基礎上進一步模擬仿真不同壞境、不同結(jié)構時冰與結(jié)構相互作用過程，為未來的抗冰結(jié)構設計奠定基礎。

3.1 在直立結(jié)構中的應用

冰與直立結(jié)構作用時，可能出現(xiàn)劈裂、彎曲、屈曲和擠壓等多種破壞形式，其中以擠壓破壞為主。岳前進等[17-18]、K?rn?等[19]、李洪升等[20]對冰與直立結(jié)構的作用過程及結(jié)構的振動進行了分析，發(fā)現(xiàn)冰與直立結(jié)構作用主要發(fā)生屈曲破壞和擠壓破壞，并且不同速度的海冰與結(jié)構作用擠壓破壞時會形成如圖2示的3種不同形式的冰激振動：準靜態(tài)振動、穩(wěn)態(tài)振動、隨機振動，海冰的擠壓破碎也因為冰速的不同而表現(xiàn)為3種不同的破碎模式。當冰速很慢時，出現(xiàn)準靜態(tài)冰力，對應冰的破碎模式為準靜態(tài)韌性破壞；隨著冰速的增加會出現(xiàn)導致結(jié)構穩(wěn)態(tài)振動的冰力，此時的破壞模式為韌脆轉(zhuǎn)換破壞；隨著冰速的繼續(xù)增加會產(chǎn)生導致結(jié)構隨機振動的冰力，對應破碎模式為脆性破壞。

(a)極限應力準則 (b)極限動量準則 (c)極限力準則圖2 直立結(jié)構最大靜冰力的3種極限準則Fig.2 Three limit rules for calculating peak static ice force on vertical structures

在靜冰力領域，目前已有相應的工程規(guī)范，楊耀鵬等[21]以新建的海上平臺為例，通過有限元軟件分析其抗冰性能，發(fā)現(xiàn)結(jié)構設計靜力的安全儲備比較大，不過冰振對結(jié)構造成的潛在影響還有待進一步考究，可能會造成結(jié)構的疲勞失效。

在動冰力的研究領域，季順迎等[22]通過離散元法模擬計算了海冰與直立海洋平臺相互作用的過程,模擬海冰對不同樁徑的直立結(jié)構的動冰力過程，確定海冰的破碎規(guī)律以及冰載荷的特性，發(fā)現(xiàn)平均冰力與樁徑呈正比。賈賓等[23]在前者研究的基礎上，采用基于近場動力學理論的數(shù)值算法,數(shù)值模擬的主要參數(shù)參考了前者的模型，模擬了海冰與柱狀結(jié)構的作用過程，并考慮了相互作用時產(chǎn)生的結(jié)構振動，探究了樁徑、冰速等參數(shù)對極值靜冰力和結(jié)構振動位移的影響，冰力部分模擬結(jié)果與前者相近，可見這種新型的數(shù)值算法也是具有合理性的。

冰激振動是影響海洋結(jié)構安全以及導致其損傷和疲勞壽命的重要因素，也是冰載荷研究領域的熱點。張大勇等[24]在分析冰激振動下直立結(jié)構的疲勞壽命時，借助有限元數(shù)值分析方法，主要考慮了穩(wěn)態(tài)振動和隨機振動造成的疲勞損傷影響，選取安全壽命設計方法[25]進行冰激疲勞估計，明確冰與直立結(jié)構相互作用過程，提出一套冰激直立抗冰平臺的疲勞壽命計算流程。

在分析冰載荷問題時往往把海洋平臺視為剛體結(jié)構，而在分析海洋結(jié)構的動力特性時則把冰載荷簡化成冰力函數(shù)處理，隨著冰激振動研究的深入，需要進一步研究海冰與海洋結(jié)構的動力耦合作用，不少研究者提出了建立離散元-有限元的耦合模型來分析冰激結(jié)構振動的內(nèi)在機理。王帥霖[26]采用基于區(qū)域分解方法的離散元-有限元耦合模型，對單樁直立結(jié)構的冰激振動進行了仿真模擬，由梁-殼組合單元構成海洋結(jié)構的有限元模型，在耦合界面上發(fā)展了相應的接觸算法和耦合參數(shù)傳遞算法，考慮了冰速、厚度變化的影響，在冰況相近的前提下，仿真計算出冰載荷、冰激振動加速度和應力分布與現(xiàn)場實測數(shù)據(jù)趨勢相同，初步說明了該耦合模型在分析冰與結(jié)構相互作用過程的適用性。

3.2 在非直立結(jié)構中的應用

近年來，錐體結(jié)構在海洋工程領域不斷被采用，不少國家為冰區(qū)的海洋平臺加裝了圖3所示的錐體結(jié)構，國內(nèi)外的研究人員針對冰排與固定傾斜結(jié)構作用時受到的冰載荷做了大量的理論、實驗和數(shù)值研究，通過海冰的物理力學特性發(fā)現(xiàn)，冰的彎曲強度遠小于擠壓強度。錐體取代直立結(jié)構被公認為降低冰力的解決辦法。

圖3 抗冰單樁基礎[24]Fig.3 Ice resistant single pile foundation[24]

對于錐體結(jié)構而言，由于冰排的彎曲強度遠小于擠壓強度，而在根據(jù)現(xiàn)在設計標準，即使是直立結(jié)構也有很大的靜力安全儲備，而錐體則不言而喻，如圖4所示，對于錐體而言，雖然冰排破壞的力大幅減小了，不過會產(chǎn)生冰與斜面的摩擦力所帶來的額外載荷，因此研究椎體結(jié)構靜冰力時不僅需要考慮錐體本身的物理參數(shù)還需要考慮冰爬坡力所帶來的額外載荷。

圖4 海冰與錐體碰撞示意[42]Fig.4 The sketch of collision between ice and cone[42]

Sand等[27]和Derradji-aouat[28]，利用有限元軟件ANSYS計算了冰排與斜坡結(jié)構作用時結(jié)構受到的冰力，將冰排理想化成作用在彈塑性基礎上的線彈性板，將得到的結(jié)果與前人基于實驗數(shù)據(jù)的理論公式進行比對。龔榆峰等[29]在各項同性海冰失效準則的基礎上結(jié)合Abaqus有限元軟件，開發(fā)用于模擬海冰失效的新型單元并且結(jié)合了Ralston[30]給出的塑性極限分析解進行比較,直觀地給出冰損傷過程及對應結(jié)構受力的變化過程。Bjrnar等[27]使用非線性有限元方法模擬冰對固定圓錐體的作用，并分析錐角及摩擦力對冰荷載的影響。利用一種特殊的接觸算法模擬了冰原與錐狀結(jié)構之間的相互作用，使跟蹤冰原與結(jié)構之間逐漸發(fā)展的接觸成為可能。在相互作用過程中，由于冰原形態(tài)的變化，浮力也隨之變化。通過引入一個連續(xù)的非線性地基模型來跟蹤這一過程，該模型包括了浮力和冰的比重的影響。用2種不同的本構模型來近似冰的力學行為。將材料非線性和冰與結(jié)構摩擦的影響也納入有限元計算當中，精確追蹤冰蓋與斜坡結(jié)構之間的接觸，包括摩擦，計算部分或全部被淹沒的冰蓋上的浮力。探究了冰與結(jié)構之間錐角和摩擦系數(shù)的變化對結(jié)構的影響，并將數(shù)值計算的結(jié)果與基于塑性極限分析的分析結(jié)果進行了比較，發(fā)現(xiàn)兩者非常接近，可見該有限元模型的可靠性。

錐體結(jié)構對減小動冰力以及減小冰激振動方面相比于直立結(jié)構表現(xiàn)出色，前文已經(jīng)提及，冰激振動是目前寒區(qū)海洋工程抗冰結(jié)構設計的面臨的重大課題，不過目前仍沒有給出計算冰激結(jié)構振動的冰力函數(shù)，對海冰引起結(jié)構振動的力學機制的解釋還沒統(tǒng)一，研究者們也在不斷探究適合計算海冰動力學的數(shù)值方法。王剛等[31]利用有限元軟件，對不同接觸寬度的冰與錐體相互作用過程進行了分析，對徑向開裂型和環(huán)向開裂型2種冰排彎曲斷裂形式進行數(shù)值仿真計算，研究了接觸寬度與冰排彎曲斷裂形式的關系。李海等[32]綜合考慮了計算的效率和精確性，發(fā)展了一種將拉格朗日坐標與歐拉坐標相耦合的方法，將有限差分法和質(zhì)點流體力學方法相結(jié)合，考慮了海冰的堆積以及渦旋風場的作用，對遼東灣海冰動力過程進行了72 h的數(shù)值模擬計算，并且與衛(wèi)星的遙感圖像進行比較，發(fā)現(xiàn)模擬的海冰的厚度、密集度分布以及冰速與衛(wèi)星遙感圖像有良好的一致性，可以較好地模擬出渤海海冰的動力學過程。

3.3 數(shù)值仿真結(jié)果的驗證

目前對錐體海洋結(jié)構與海冰的相互作用主要通過現(xiàn)場實測、室內(nèi)試驗以及數(shù)值模擬展開研究，由于冰-錐相互作用過程的復雜性，很難通過現(xiàn)場實測和試驗分析不同參數(shù)對該過程的影響，而數(shù)值模擬方法可以準確修改參數(shù)，再現(xiàn)出實測情況的冰-錐相互作用過程，并在確定數(shù)值仿真方法可行性的基礎上進一步改變參數(shù)模擬各種情況，從而達到冰激振動預警分析、疲勞分析以及動力優(yōu)化設計的目的。以上數(shù)值仿真工作大多基于已有的現(xiàn)場實測或室內(nèi)試驗數(shù)據(jù)仿真研究以確定數(shù)值仿真工作的合理性。關湃等[33]基于模型試驗結(jié)果，借助有限元軟件對冰錐結(jié)構進行優(yōu)化設計，模擬不同冰況下結(jié)構的響應，不斷改變結(jié)構的各種形狀參數(shù)以達到結(jié)構優(yōu)化的目的。也有不少研究者專門設計試驗以驗證數(shù)值分析的結(jié)果是否合理，汪震宇[34]用有限元方法分析斜坡與冰的相互過程，并設計了相應的模型試驗進行誤差分析以及模型的校正。

研究冰載荷的各種研究方法都是相互驗證，不斷改進的過程，而隨著數(shù)值分析方法不斷發(fā)展成熟，而且能夠精確地改變各種參數(shù)以適應需求，對理解冰與結(jié)構相互作用機制具有成本低效率高的優(yōu)勢。

4 在海上風機一體化數(shù)值分析的應用

4.1 海上風機一體化分析需求

海上風機也作為一種海洋結(jié)構，在冰載荷方面的研究與前述的直立結(jié)構和錐面結(jié)構有諸多相似之處，在這便不再贅述了。海上風機相比于其他海洋結(jié)構，最大的特點在于風機為典型的柔性結(jié)構，海冰會引起柔性結(jié)構明顯的動力放大，冰振響應十分顯著，對結(jié)構及上部設備造成一定影響[35]。張大勇等[36]基于有限元分析法研究動冰載荷對結(jié)構的影響，選取單立柱式三樁形式的風電基礎并與相同結(jié)構形式的導管架平臺進行動力特性對比，文章中歸納總結(jié)了冰載荷作用下風電基礎的各種失效模式，包括靜冰力作用下的結(jié)構失效和動冰力作用下的疲勞失效、振動失效和土壤失效，發(fā)現(xiàn)現(xiàn)有的風電基礎設計雖能滿足極值冰載荷的要求，但在動冰力作用下對結(jié)構造成的影響還有待考量。通過有限元法仿真模擬可以得到驗證風電基礎結(jié)構主要是塔架在載荷下的響應，但并沒有模擬出動冰力振動下對內(nèi)部設備如電機、限速安全機構等的影響，但可以通過結(jié)構的冰振位移、速度等指標來間接地考慮動冰力下結(jié)構振動對內(nèi)部設備的影響程度，為寒區(qū)風電基礎的抗冰設計的改良、相關的實驗設計及安全保障提供合理依據(jù)。

4.2 常用海上風機一體化分析軟件及應用

目前用于海上風機一體化的軟件主要為FAST軟件，作為風機耦合的軟件，不僅能夠建立基礎結(jié)構的模型還能搭建旋轉(zhuǎn)葉片-機艙-塔架組件的詳細模型，更加貼近實際的模擬仿真風機在各種載荷作用下的響應。FAST軟件[37]是美國可再生能源實驗室(National Renewable Energy Laboratory，NREL)開發(fā)，用于計算氣動彈性及結(jié)構動力學仿真的開源軟件。該軟件的仿真流程圖如圖5所示，F(xiàn)AST作為一款開源的軟件，為拓展研究、增加功能提供了基礎，部分研究者在軟件基礎上進行二次開發(fā)以適應不同的研究需求。

現(xiàn)代風力機作為最大的旋轉(zhuǎn)機械，扇葉長度甚至可以超過最大的飛機翼展，風載荷對風機的影響是其他海洋結(jié)構無法比擬的，所以分析風和海冰的聯(lián)合載荷與風機作用的過程十分重要。許子非等[38]考慮了海冰與湍流風載荷的聯(lián)合作用，利用FAST軟件，以近海單樁風力機為研究對象，數(shù)值計算在不同厚度及海冰移速的海冰碰撞條件下塔架的動力學響應特性，并對比了有無抗冰錐時響應的不同，發(fā)現(xiàn)抗冰結(jié)構能夠明顯減小結(jié)構的疲勞載荷。

除了FAST以外，目前用上風機一體化數(shù)值分析主流的軟件還有HAWC2[39](horizontal axis wind turbine simulation code 2nd generation)和Bladed，都是被一些風力機廠商所認證的仿真軟件，作為風機仿真軟件應用于很多研究項目和工程當中，能夠模擬許多現(xiàn)有的風力渦輪機。

HAWC2能夠綜合考慮空氣動力學載荷、水動力載荷和土力學載荷，在時域內(nèi)使用線性和非線性的方法結(jié)合以計算結(jié)構動力響應，基于多體公式的非線性氣動彈性模型，可以處理復雜的結(jié)構。Shi等[40]利用了HAWC2代碼，建立了半經(jīng)驗的全耦合冰載荷模型，研究冰載荷對安裝了抗冰結(jié)構的海上風電機的影響，同時考慮了風載荷和冰載荷的基礎下進行了耦合分析并且與非耦合分析的結(jié)果進行了比較，討論了耦合分析的重要性，并且研究了風速、平均冰厚、漂移速度不同時，冰的物理性質(zhì)和破冰錐幾何形狀的對海上風機響應的影響。

圖5 FAST軟件仿真流程圖[38]Fig.5 Software simulation flow chart[38]

除了冰與結(jié)構的相互作用之外，海上風力機仍需考慮結(jié)冰對風力機的影響，結(jié)冰會引起額外的載荷和振動，不同葉片上冰載荷的不平衡也會在一定程度上增加機組部件的疲勞。易賢等[41]采用數(shù)值計算的方法進行結(jié)冰計算，模擬結(jié)冰形成的動態(tài)過程，分析結(jié)冰的主要部位以及厚度，還有對風機氣動特性的影響，為設計結(jié)冰探測方法和除冰方案奠定基礎。

風力機作為高聳結(jié)構，相比于其他海洋結(jié)構柔性特性更為明顯，也意味著風力機受冰激振動的影響更為顯著，而目前對于冰激振動的研究仍處于探索階段，離徹底揭示冰激振動機理和抗冰振結(jié)構設計的完善還有很長的路要走，隨著數(shù)值仿真技術的成熟，數(shù)值計算方法相比于原型觀測和室內(nèi)試驗研究的進程更短，花費小效率高，能規(guī)避很多不必要的風險。

5 結(jié)束語

目前冰載荷研究的靜冰力領域發(fā)展相對成熟，各個研究地區(qū)均有明確的規(guī)范，不過仍未形成統(tǒng)一計算公式或是數(shù)學模型，靜冰力分析理論還有待進一步發(fā)展。新型的海洋結(jié)構或是抗冰結(jié)構對于靜冰力均有較大的安全儲備，應當探索降低靜冰力工程設計值以提高經(jīng)濟效益的可能性；在動冰力和冰激振動領域，仍需探索可以用于工程的動冰力函數(shù)，進一步優(yōu)化抗冰結(jié)構的減震性能，形成相關的工程規(guī)范，并且發(fā)展數(shù)值分析方法以適應動冰力和冰激振動數(shù)值計算需求。

理論分析和數(shù)值模型的發(fā)展必然離不開工程實例與試驗，目前我國僅天津大學有冰池實驗室以及哈爾濱工程大學建設的室外冰水池實驗室，但與德國漢堡水池相比仍有著顯著的差距，有待進一步改善實驗條件，開展更多冰力實驗。

當然，冰載荷領域的研究要想應用到海上風機還需要結(jié)合其他領域知識，諸如空氣動力學、水動力學、土力學等領域的研究，并結(jié)合現(xiàn)代的數(shù)值仿真技術，不斷發(fā)展一體化數(shù)值分析技術以貼近實際。

相信隨著數(shù)值模擬技術與計算機的蓬勃發(fā)展，以及目前各國對冰載荷愈加重視，冰載荷的數(shù)值仿真方法也會越來越成熟，同時也會有更多的文獻方便學習。我國渤海海域能源豐富，不論是傳統(tǒng)能源還是新能源的開采都避不開樁式結(jié)構，而渤海易受海冰影響，因此我國在抗冰方面的步伐從未停止，特別是近年來各地興起了海上風電的熱潮。通將過數(shù)值方法與模型試驗、現(xiàn)場觀測等手段結(jié)合不斷發(fā)展冰力研究，并在工程實踐中不斷積累相關經(jīng)驗，相信在不久之后便會形成一套完善海上風機的冰力規(guī)范。