風機基礎TLP平臺在波浪中運動性能的CFD數值分析

安筱婷，趙偉文，萬德成

(上海交通大學 船舶海洋與建筑工程學院 海洋工程國家重點實驗室 船海計算水動力學研究中心(CMHL)，上海 200240)

海上風能作為一種清潔的可再生資源，在陸上資源日益枯竭的今天越來越受到關注，海上風機是人類開發海洋風能資源的載體。隨著風能開發向深水發展，支撐風機的載體平臺越來越受到關注。深水風機的優勢在于深水環境中風速更大，風能更強，且不影響航運，但隨著水深的增加，作用在載體平臺上的水動力增加，淺水中常用的固定式平臺由于造價高、質量大、建造安裝難度大、運動性能不適合惡劣海況等原因不再適用。浮式支撐平臺可重復使用，造價低，運動性能更好，因而在深水環境中應用廣泛。

浮式平臺中，張力腿平臺(tension leg platform，簡稱TLP)是目前被廣泛應用的一種平臺，使用剩余浮力平衡張緊式系泊系統的預張力，具有半順應半剛性的特點，可抵抗波浪引起的運動，以此為風力發電提供一個相對穩定的環境，減少風力發電機功率的變化對電網的影響。同時TLP自身的直立浮筒結構也為平臺提供了良好的運動性能。

目前已有多座TLP實際投入生產使用，此外，研究學者也針對作為風機基礎平臺的TLP的特性需求，提出了許多設計方案，并針對這些設計方案進行了大量計算與驗證。Ren等[1]對一種TLP—海上風機系統進行了1/60比例模型試驗,以考察耦合風浪效應，同時還提出了TLP系泊系統的優化方案。Bachynski等[2]考慮100～200 m水深中作業的TLP實際情況，提出了該水深水平中TLP的設計方案。Adam等[3]提出了20～350 m水深中，將TLP作為風機基礎平臺的設計——GICON?設計。Gao等[4]基于邊界元方法,對風機基礎TLP 的運動響應和波浪力機制進行了分析。Bilgili等[5]使用全耦合模擬對于單立柱TLP風機的設計參數進行了參數研究。Han等[6]提出了一種半潛式TLP平臺(submerged tension leg platform，簡稱STLP)，此平臺在作業狀態下水線面面積較小，因此具有更佳的動力學性能。

與船舶不同，浮式平臺需要長期保持在特定海域進行作業，因此需要系泊系統提供回復力來抑制外界載荷引起的水平運動。另一方面，風機支撐平臺也不允許發生過大幅度的運動響應，這樣才能保證風機發電正常輸出，電網穩定。故要求浮式平臺配備相應的系泊系統，以擁有位置保持能力[7]。綜上所述，浮式平臺工作時，需要響應系泊系統配合。

計算機科學技術的不斷發展與數值計算方法的不斷優化，使得船舶與海洋工程水動力學計算領域中應用計算流體動力學(computational fluid dynamics，簡稱CFD)方法進行計算成為可能。使用自主開發的基于開源程序庫OpenFOAM開發的naoe-FOAM-SJTU求解器，該求解器可以求解船舶與海洋工程水動力學問題[8]。為了求解系泊結構物的運動響應，劉遠傳[9]以此求解器的6自由度運動模塊為基礎，聯系了浮式結構物運動與系泊系統的受力，進一步開發出了系泊系統計算模塊(naoeFOAM-ms)，完善了原求解器的功能。

由于風機平臺系統在實際投入使用時海洋環境條件非常復雜，研究不同海況下風機基礎平臺的運動響應特性對于指導風機系統的設計、保障人身與財產安全、保證風機系統正常工作有著重要意義。模擬中考慮了不規則波的高階成分，相比于目前已有的工作更貼近海洋環境載荷實際情況。采用naoe-FOAM-SJTU求解器分析了一座STLP的水動力性能，計算了該平臺的固有運動性能，并以此為基礎分別討論處于中等水深作業狀態下以及極端海況生存狀態下的運動響應情況。首先計算了STLP的固有周期，并與已有結果進行對比，驗證了求解器的有效性。在此基礎上分別研究了STLP在中等海況下和極端海況下的運動響應、平臺受到彎矩及系泊張力，研究了在考慮非線性波浪載荷的極端海況下與一般作業海況下STLP的運動響應情況，計算了不同工況中的風機基礎所受彎矩及錨鏈受力情況，并詳細展示了流場、速度場信息，分析了高階波浪成分、不同海況等條件對于STLP運動性能的影響。研究結果表明，naoe-FOAM-SJTU求解器可以有效地模擬TLP平臺在波浪中的運動響應問題，還可以詳細展示流場信息。

1 數學模型

1.1 控制方程及離散方法

關于OpenFOAM的求解原理和具體求解步驟可參見文獻[10]和文獻[11]，這里作簡單概述。文中的數值模擬基于不可壓黏性流體，其控制方程為:

·U=0

(1)

(2)

式中：U代表流場速度，Ug代表網格節點的運動速度；pd=p-ρg·x代表流場動壓力，x為網格節點位置；g、ρ和μ分別為重力加速度、流場的密度和動力黏性系數；fσ為源項，表示自由面上的表面張力。

由于模擬中無來流，雷諾數小，湍流的影響可以忽略，因此采用層流模型以提高計算效率，忽略了結構物表面的邊界層；采用VOF(volume of fluid)方法捕捉自由面[11-12]，定義每個網格中的體積分數，通常選擇ɑ=0.5等值面作為流場中的自由面。計算中，采用有限體積法(finite volume method，簡稱FVM)進行空間離散，使用隱式Euler格式進行時間離散，使用二階TVD格式進行對流項離散，使用二階中心差分格式進行擴散項離散。在進行離散化后，應用PISO算法來迭代求解和解耦壓力及速度。PISO算法由 Issa提出[13]，OpenFOAM中的 PISO算法采用的是由Rhie提出的同位網格方法[14]。最后，可以獲得每個時間步的流場變化及最終模擬結果。

1.2 動網格技術

為了模擬平臺的運動,需要不斷變換計算網格以適應結構物運動對流場的影響。文中采用了naoe-FOAM-SJTU中的動網格技術。當平臺發生位移時，網格的形狀由移動節點改變，同時保持每個網格節點與相鄰網格節點之間的拓撲關系不變。通過求解拉普拉斯方程可以得到網格節點的位置：

·(κXg)=0

(3)

其中，Xg是節點的位移，κ是節點變形系數：

(4)

其中，r是網格中心與移動邊界的距離。

1.3 造波與系泊系統模擬

影響海洋工程結構物運動最顯著的環境因素之一就是波浪。文中的造波通過開源代碼waves2foam[15]實現。waves2foam代碼為了產生實際模擬中所需的具體波浪，改變了速度入口邊界條件。文中模擬的各種波浪情況均為一階不規則波(Jonswap譜)。在實踐中，線性波浪理論被用來模擬不規則波，例如沿x軸正向傳播的長峰不規則波的波面升高[7]：

(5)

式中：Aj為第j個單元波的波幅,ωj為第j個單元波的圓頻率，εj和kj分別為相應的隨機相位角及波數。波幅與波譜關系為：

(6)

Jonswap波譜公式為：

(7)

式中：H1/3和T1分別為波浪的有效波高和平均波浪周期，峰值增強因子記為γ。

使用彈簧模擬方法來模擬STLP的垂向系泊系統，具體實現利用了naoe-FOAM-SJTU中的系泊系統求解模塊[8]。彈簧模擬方法把STLP的垂向系泊系統等效為只可拉伸不可壓縮的一根彈簧，且不考慮彈簧自身的重力以及流場力，彈簧模擬方法目前被廣泛應用于TLP類平臺運動響應計算中。

彈簧模擬方法的基本原理為胡克定律：

F=kΔl

(8)

其中，k是彈簧的彈性系數，在naoeFOAM-ms的字典文件中給定；F為張力結果，即系泊系統供給浮式結構物的回復力；彈簧的伸長量記為Δl，程序中即為系泊纜實際長度(系泊點與錨泊點的距離)與系泊纜初始長度之差。在系泊系統求解模塊中，t-1時刻的系泊力是t時刻的邊界條件，參與6自由度運動方程求解，得到本時刻的結構物位置，在錨泊系統中通過系泊點位置(與浮式結構物的位置有關)即可獲得下一時間步的系泊力，作為外力參與6自由度運算。

2 計算模型設定

文中的計算模型參考Han等[6]的工作，風機(NREL-5MV)[16-17]、風機支撐平臺與8條系泊纜共同組成風機發電系統，如圖1所示。STLP作為風機的水下支撐平臺，由中心立柱、4個垂直立柱、水平浮筒和斜撐構成，見圖2。計算域布置如圖3所示，垂直和水平浮筒完全淹沒在平均海平面(MSL)下，吃水為20 m，計算域深度為100 m。模擬中坐標軸原點o設置在STLP的水線面中心，三條坐標軸遵循右手螺旋法則，z軸垂直向上。

圖1 風機系統示意Fig. 1 Diagram of turbine system

圖2 STLP示意Fig. 2 Diagram of STLP

圖3 計算域布置示意Fig. 3 Diagram of calculation domain

限制平臺運動的錨鏈系統主要參數如表1所示，每2根錨鏈分為一組，每組錨鏈系泊在平臺4個垂直浮筒下緣外側，STLP平臺的主要參數如表2所示，為了計算系統整體參數，NREL-5MW風機的主要參數也一并給出。需要說明的是由于風機平臺模型構造上的差異，風機系統整體關于質心的橫搖、縱搖轉動慣量與文獻設定的不相同，其他參數與文獻[6]中的設置保持一致。

此算例中以作用于浮式風機系統重心上的定常外力來模擬風載荷。模擬極端海況下平臺運動實際上是模擬風機的生存條件，因此此時下風機處于關閉狀態。為了簡化模擬，一般海況下同樣簡化為關閉狀態，只考慮定常風和波浪條件對于平臺運動響應的影響。

表1 系泊系統主要參數Tab. 1 Main parameters of mooring system

表2 平臺與NREL-5MW風機的主要參數Tab. 2 Main parameters of platform and NREL-5MW turbine

圖4 計算網格Fig. 4 Calculation grids

為了保證數值模擬中可以準確生成波浪并穩定傳播，捕捉流場內部變化、流場自由面變化等，需要對流場內部及結構物表面進行網格生成，并對部分區域進行網格加密。使用snappyHexMesh工具[18]生成網格，可以捕捉出物面邊界條件，并加密背景網格，從而生成質量較高的網格。具體流程為：首先提取風機基礎平臺模型的特征邊，并對自由面范圍、平臺表面進行網格加密，以便更精準地捕捉。計算網格如圖4所示，STLP平臺表面與附件的網格如圖5、6所示。模擬中網格量約為260×104，為滿足CFL條件，計算時間步長為0.005 s。

圖5 水下平臺表面網格Fig. 5 Grids on the surface of submerged platform

3 計算結果與分析

3.1 STLP自由衰減驗證

通過對STLP單自由度自由衰減特性的模擬來驗證naoe-FOAM-SJTU求解器的準確性和可靠性。平臺的外形對稱，故其縱蕩周期等于橫蕩周期，橫搖周期等于縱搖周期，模擬中只計算縱蕩與縱搖周期。在計算STLP自由衰減特性中，風機葉片固定且無外部風浪載荷，在單自由度的條件下將STLP以某一設定初始速度值釋放以觀察平臺的固有周期。

圖7展示了不同自由度下平臺自由衰減運動隨時間變化的情況以及其振動能量的頻域分布，STLP運動的固有周期可根據該時歷曲線得到。CFD模擬(naoe-FOAM-SJTU)的固有周期結果與勢流理論計算(BModes、FAST)的固有周期結果如表3所示，其中BModes為有限元模態分析結果，FAST為Han等[6]使用FAST(NREL設計的風機載荷分析軟件)進行的全耦合時域模擬結果，以BModes為參考的模擬誤差由下式計算：

(9)

圖7 STLP各自由度自由衰減運動時歷曲線及能量譜結果Fig. 7 Time curves and energy spectrum results of free decay motion of the STLP

表3 STLP運動的固有周期Tab. 3 The natural period of motion of the STLP

由表3可知，使用naoe-FOAM-SJTU得到的模擬結果與已有文獻結果吻合良好，誤差均在3%以內(除縱搖(橫搖)運動外)。其中縱搖(橫搖)的自由衰減表現與文獻結果偏差較大，這是因為本模擬中的風機模型與文獻不同，因此對于縱搖(橫搖)旋轉軸的轉動慣量與文獻[6]中不同。該STLP的縱蕩(橫蕩)的自有周期大于25 s，垂蕩、縱搖(橫搖)的自有周期小于3.5 s，完全符合浮式平臺經驗設計標準建議[19]。若平臺的固有周期接近波浪頻率，可能發生波浪載荷引起共振的危險。比較可知，相比于典型波浪頻率，該STLP的縱搖(橫搖)、垂蕩運動的固有頻率偏高，橫蕩(縱蕩)、艏搖運動的固有頻率偏低，所以共振危險較低，作為風機發電基礎平臺可以達到穩定輸出的目的。模擬結果與文獻結果吻合良好，證明了naoe-FOAM-SJTU求解器模擬此類問題的可靠性。

3.2 波浪場中STLP水動力運動響應結果

研究不同波浪場中平臺的運動特性，初始時刻將平臺自靜止釋放所有自由度運動，研究平臺在波浪中運動的特性，繪制不同自由度運動、彎矩、系泊張力的時歷曲線，并計算其振動能量的頻域分布。最后將CFD模擬結果與已有結果對比。STLP結構與外界環境載荷都是對稱的，因此文中只討論縱蕩、縱搖及垂蕩運動響應。

3.2.1 極端海況下STLP水動力運動響應結果

此海況是平臺的生存海況，風機處于關閉狀態。將平臺置于極端海況下，波浪譜采用Jonswap譜，有效波高14.4 m，譜峰周期13.3 s，平均風速49 m/s，將各自由度的運動響應結果進行統計，結果如表4所示。

表4 極端海況下STLP運動響應統計值Tab. 4 STLP motion response statistics under extreme sea conditions

由于波浪及定常風力的聯合作用，在6個自由度運動中，縱蕩運動最顯著，運動幅值最大，波動也最為劇烈。由于TLP平臺的自身性質，平面外運動(橫搖、縱搖及垂蕩)的幅值很小；由于系泊系統、環境條件的對稱性，橫蕩和艏搖運動的幅值也很小。將此海況下CFD結果與全耦合時域分析結果(考慮二階波浪)的運動響應的平均值進行對比，如圖8所示。

圖8 極端海況下STLP的運動及系泊受力統計特性Fig. 8 Statistical characteristics of STLP motion and mooring tension under extreme sea conditions

由圖8分析可知，使用CFD方法模擬極端海況下STLP的運動，相比于全耦合時域分析，運動幅度及系泊纜張力都偏大。這是更高階波浪載荷作用在風機系統的影響所致。

將各系泊錨鏈的系泊力值進行統計，如表5所示，可見錨鏈3～6受到的平均張力最大，這是由于這4條錨鏈位于迎浪面。迎浪、背浪的錨鏈張力差距平均差值可達50 kN。

表5 系泊系統受力平均值Tab. 5 Average tension of the mooring system

將平臺受到的波浪載荷、運動響應、所受彎矩和系泊系統受力的時歷曲線進行快速傅里葉變換(FFT)處理，得到圖9至圖14。

圖9為STLP在極端海況下的波浪載荷時歷曲線與能量譜結果，可見CFD模擬考慮了波浪的各階成分。圖10為STLP在極端海況下的縱蕩運動時歷曲線與能量譜結果，此時在運動中波浪頻率響應占主要成分，縱蕩共振響應也較為明顯。波浪頻率響應成分遠大于縱蕩共振響應，雖然此時縱蕩運動的固有頻率(0.0336 Hz)接近波浪頻率范圍。這一特性與Han等[6]的全耦合時域分析結果相似。與全耦合時域分析結果對比可知，在CFD模擬中，高階波浪成分更多，能量更分散，因此波浪頻率成分峰值相比于全耦合模擬結果偏小，縱蕩共振響應也更分散。這是因為圖10(d)中最多只考慮到波浪二階載荷，而圖10(c)中波浪的各階成分均考慮在內。對比可知，縱蕩運動受到高階波浪荷載的影響不敏感。

圖9 極端海況下平臺波浪載荷隨時間變化曲線與能量譜結果Fig. 9 Time curves and energy spectrum results of wave load under extreme sea conditions

圖10 極端海況下平臺縱蕩運動隨時間變化曲線與能量譜結果Fig. 10 Time curves and energy spectrum results of surge motion of the STLP under extreme sea conditions

在縱蕩運動隨時間變化的曲線中，發現STLP平臺的縱蕩平衡位置不在0位置處，有明顯漂移，這是定常風載荷的作用影響。CFD模擬中將湍流風場等效為定常風場，結果與全耦合模擬中的結果相差不大，可知在極端海況下，風機關停狀態下風場對于STLP的縱蕩運動影響不大。

圖11為STLP在極端海況下的垂蕩運動時歷曲線與能量譜結果，可見在極端海況下，平臺的垂蕩運動響應中主要成分是波浪頻率響應和縱蕩共振響應，其中波浪頻率為主導。因為平臺初始為直立狀態，垂蕩會與縱蕩在運動中發生耦合現象，因此能量譜中會出現縱蕩共振成分。由于TLP系泊系統的固有特性，垂蕩運動的平衡位置低于水平面。

極端海況、湍流風場中的STLP垂蕩運動能量譜結果中(圖11(d))明顯可見風頻率響應，但在CFD模擬中由于將風場等效為定常場，所以無此成分，說明極端海況下風場作用對于處于關停狀態的風機STLP系統的垂蕩運動影響較大。

圖11 極端海況下平臺垂蕩運動隨時間變化曲線與能量譜結果Fig. 11 Time curves and energy spectrum results of heave motion of the STLP under extreme sea conditions

圖12為STLP在極端海況下的縱搖運動隨時間變化曲線及能量譜結果。在此情況下，能量譜中占主導成分的分別是波浪頻率響應和縱搖共振響應成分。在CFD模擬中，除一階、二階波浪成分外，還可見更高階的波浪成分，而全耦合時域模擬只有二階波浪成分；由CFD結果可知，平臺運動的縱搖共振響應成分由于波浪中的高階成分的存在而變大，且響應的頻率范圍也更大。由上述分析可知，在此情況下波浪的高階載荷對STLP的縱搖運動影響較大。與縱蕩結果類似，風場對于極端海況下風機關閉狀態下的STLP的縱搖運動影響不明顯。

圖12 極端海況下平臺縱搖運動隨時間變化曲線與能量譜結果Fig. 12 Time curves and energy spectrum results of pitch motion of the STLP under extreme sea conditions

由以上各圖可知，平臺在極端海況產生運動時，縱蕩運動幅值較大，縱搖、垂蕩運動幅值較小，平面內近似柔性，平面外近似于剛性。這種運動特性的產生是TLP自身的系泊特性決定的：TLP的張力腿平衡了大于重力的浮力，因此系泊纜始終是受拉狀態，限制了平面外運動的幅值；STLP的主體是一個直立的浮筒，在吃水方向尺寸要遠大于在水平面內的尺寸，所以受到的垂直方向波浪力要小于水平方向上的波浪力，因而在平面內具有柔性的特點[20]。STLP的這種特點使其可以很好地抵抗波浪引起的運動，并為風力發電提供一個相對穩定的環境。

圖13是極端海況下平臺受到的y向彎矩隨時間變化曲線與能量譜結果，其中y向彎矩是針對風機基座而言的。由圖可知在彎矩受力能量分布中，占主要成分的是波浪頻率響應和平臺的縱搖共振響應。將CFD結果與時域分析結果對比，發現CFD模擬中彎矩FFT結果的縱搖共振成分有顯著升高，對應的頻率范圍也更廣。出現這種現象是因為平臺y向彎矩主要受縱搖運動影響，而CFD模擬中全面地考慮了波浪更高階成分，高階波浪載荷引起了更大的縱搖運動，使得風機的重力對于y向彎矩的貢獻更大。

圖13 極端海況下平臺y向彎矩隨時間變化曲線與能量譜結果Fig. 13 Time curves and energy spectrum results of y bending moment of the STLP under extreme sea conditions

錨鏈4位于迎浪面，圖14為STLP在極端海況下錨鏈4的張力時歷曲線及其能量譜結果。錨鏈4張力的能量譜中，占據主導地位的是波浪頻率響應，且可見平臺的縱搖、垂蕩共振響應成分。而對比發現，全耦合結果中縱搖成分幾乎不可見，說明高階波浪載荷將引起系泊張力上更大的共振響應。

圖14 極端海況下錨鏈4張力隨時間變化曲線與能量譜結果Fig. 14 Time curves and energy spectrum results of #4 mooring tension of the STLP under extreme sea conditions

以上結果顯示出：極端海況下STLP運動響應主要是波浪頻率成分主導，也有共振響應成分存在，且高階波浪載荷對于縱搖共振響應的增加效果較為顯著，高階波浪載荷將引起y向彎矩與系泊張力上更大的共振響應；風場對于極端海況下處于關停狀態的風機垂蕩運動響應影響較大；CFD結果與部分勢流理論結果相比，可以明顯體現出波浪中更高階成分對于平臺運動性能的影響，還可以可視化地顯示出流場詳細情況，如圖15和圖16所示。

圖15 極端海況下瞬時流場Fig. 15 Instantaneous flow field under extreme sea conditions

圖16 STLP在極端海況下不同時刻的流場狀態Fig. 16 The state of the flow field at different times under extreme sea conditions

圖15為極端海況下平臺附近的瞬時流場，圖15(a)為波谷經過平臺，圖15(b)為波峰經過第一排斜撐，圖15(c)為波峰經過第二排斜撐。由圖可以看出由于平臺斜撐的存在，波浪在傳播到斜撐時將被阻塞，產生局部次生波。當波峰到達斜撐時，在斜撐上產生了爬高現象。波浪繼續傳播，被斜撐阻塞的波浪離開斜撐，繞過斜撐繼續傳播，與入射波疊加，從而將導致波浪變形。

圖16為極端海況下不同時刻的流場狀態。可以看出，對于文中使用的二維不規則波，在經過一定時間的演化及與平臺的相互作用后，流場變得非常復雜，整體波面受影響很大，已具有了三維特性。

3.2.2 中等海況下STLP水動力運動相應結果

將STLP置于中等海況下，波浪有效波高3 m，譜峰周期10 s，平均風速11.4 m/s，STLP運動響應統計值見表6。文中模擬只考慮風機基礎平臺的運動性能，暫不考慮與風場的耦合作用，因此關閉風機，使用作用在風機系統重心處的定常風來等效湍流風場。

表6 中等海況下STLP運動響應統計值Tab. 6 STLP motion response statistics under moderate sea conditions

與表4相比，由于環境較為緩和，中等海況統計值各項數據都比極端海況下的小很多。主要運動仍為縱蕩運動。將各系泊纜所受張力的均值進行統計，如表7所示。

表7 系泊系統受力平均值Tab. 7 Average tension of the mooring system

由表7可知中等海況下各系泊纜繩的張力分布非常均勻。說明此時系泊系統可以充分發揮能力，比較不易發生斷裂情況。

將平臺受到的波浪載荷、運動響應、所受彎矩和系泊系統受力的時歷曲線進行FFT處理，模擬結果如圖17～22所示。

圖17為STLP在中等海況下的波浪載荷時歷曲線與能量譜結果，可見CFD模擬考慮了波浪的各階成分。

圖17 中等海況下平臺受到的波浪載荷隨時間變化曲線與能量譜結果(naoe-FOAM-SJTU)Fig. 17 Time curves and energy spectrum results of wave load under moderate sea conditions (naoe-FOAM-SJTU)

圖18為STLP在中等海況下的縱蕩運動時歷曲線與能量譜結果，結果與極端海況下結果類似，主要成分為波浪頻率響應，縱蕩頻率響應也較為明顯。相比于圖10中在極端海況下的結果，由于有效波高的降低，波浪頻率成分減少，因此中等海況下縱蕩共振響應更顯著。

圖18 中等海況下平臺縱蕩運動隨時間變化曲線與能量譜結果Fig. 18 Time curves and energy spectrum results of surge motion of the STLP under moderate sea conditions

圖19為 STLP在中等海況下的垂蕩運動時歷曲線與能量譜結果。發現在中等海況下，垂蕩運動主要成分是波浪頻率響應和垂蕩共振響應，且可見縱蕩共振響應，其中波浪頻率響應最為明顯。與圖11對比可知，由于有效波高的減少，波浪頻率成分減少，垂蕩運動中的垂蕩共振成分更加明顯。與全耦合時域模擬結果對比可知，高階波浪在中等海況下對垂蕩共振運動的影響明顯，使得垂蕩共振運動對應的頻率范圍更大。

圖19 中等海況下平臺的垂蕩運動隨時間變化曲線與能量譜結果Fig. 19 Time curves and energy spectrum results of heave motion of the STLP under moderate sea conditions

圖20為STLP在中等海況下的縱搖運動時歷曲線及能量譜結果，與極端海況結果類似，縱搖運動的主要成分是波浪頻率響應及縱搖共振響應，但與極端海況相比縱搖成分更顯著。風機關停使得縱搖運動中幾乎沒有縱蕩共振成分和垂蕩共振成分，運動相對簡單。

圖20 中等海況下平臺的縱搖運動隨時間變化曲線與能量譜結果Fig. 20 Time curves and energy spectrum results of pitch motion of the STLP under moderate sea conditions

圖21是STLP在中等海況下y向彎矩時歷曲線與能量譜結果，其中主要成分為縱搖共振成分，同時可見波浪頻率成分。相比于極端海況，由于縱搖成分的增加，中等海況下y向彎矩中縱搖成分顯著增加。

圖22為STLP在中等海況下錨鏈4的張力時歷曲線及其FFT處理后的能量譜，其中能量譜主要成分為波浪頻率響應、縱搖及垂蕩運動頻率成分，垂蕩成分范圍較大，這是因為系泊張力主要是由垂蕩運動和縱搖運動決定的。

圖21 中等海況下平臺受到的y向彎矩隨時間變化與能量譜結果Fig. 21 Time curves and energy spectrum results of y bending moment of the STLP under moderate sea conditions

圖22 中等海況下平臺的4號錨鏈張力隨時間變化曲線與能量譜結果Fig. 22 Time curves and energy spectrum results of #4 mooring tension under moderate sea conditions

以上結果顯示出風機關停狀態下處于中等海況的STLP運動響應主要是波浪頻率成分主導，也有共振響應成分存在，由于有效波高減小，相比于極端海況，各自由度運動的共振響應更加明顯。高階波浪在中等海況下對垂蕩共振運動的影響明顯，使得垂蕩共振運動對應的頻率范圍更大。

由于不考慮風機作業的影響，CFD模擬的平臺各自由度運動結果仍是波浪頻率占主要成分。由于風場作用在葉片上的力沿縱蕩方向，因此全耦合時域模擬結果中縱蕩成分明顯，對比說明在中等海況下作業中風機的影響決定了STLP的各自由度運動響應，考慮風機的作用影響則風機基礎平臺的運動更加復雜，對風機和基礎平臺的強度、抗疲勞性等性能的要求更高。

STLP在中等海況下的平臺附近瞬時流場如圖23所示，圖23(a)為波峰經過第一排斜撐，圖23(b)為波峰經過第二排斜撐；中等海況下不同時刻的流場如圖24所示。

圖23 STLP在中等海況下的平臺附近瞬時流場Fig. 23 Instantaneous flow field near platform under moderate sea conditions

圖24 STLP在中等海況下不同時刻的流場狀態Fig. 24 The state of the flow field at different times under moderate sea conditions

由圖23、24可知由于平臺的存在，入射的二維不規則波浪在平臺附近的傳播受到了影響。平臺斜撐具有阻塞作用，波浪在到達第一排斜撐時出現波面抬升，繞過斜撐后匯合，在第二排斜撐的阻塞作用下再次爬高，導致波峰線出現明顯變形，不再為直線。

為了更清晰地分析波面，研究平臺附近的流場及速度場。圖25為T=375 s時平臺附近流場及速度場，此時將平臺隱去。可見由于平臺斜撐的阻塞作用，波浪在斜撐附近輻射、繞射，形成了局部的次生波。分析速度場可知，平臺斜撐的存在對于流體速度影響極大。

圖25 T=375 s時平臺附近流場及速度場Fig. 25 Flow field and velocity field near the platform at T=375 s

4 結 語

對一座STLP的運動特性進行了CFD數值模擬。首先對STLP的自由衰減運動進行了CFD模擬，計算出STLP的固有周期，并將結果與全耦合時域分析結果對比，發現二者基本吻合，驗證了利用naoe-FOAM-SJTU求解器處理此類平臺在波浪環境下運動問題的準確性和可靠性；在此基礎上分別模擬了兩種不同海況下平臺的各自由度運動響應、風機對平臺彎矩及系泊系統受力的時歷特性及能量譜特性，并對結果進行了分析，得出以下結論：

1) 文中研究的TLP在波浪中的各自由度運動、風機對于底座的y向彎矩及系泊系統張力的主要成分是波浪頻率響應及當前運動共振響應，除此之外，也耦合了其他自由度運動。

2) 極端海況下，平臺的縱蕩受波浪高階成分的影響較小，縱搖比較容易受到高階成分影響；中等海況下高階波浪載荷對垂蕩運動的影響明顯。

3) 極端海況下風機處于關閉狀態，風場對于STLP垂蕩運動影響較大，對縱蕩、縱搖影響很小；中等海況下風機的工作很大程度上決定了STLP的各自由度運動響應，考慮到風機的作用影響則風機基礎平臺的運動更加復雜，對風機和基礎平臺的強度、抗疲勞性等性能的要求更高。

4) 中等海況由于有效波高較小，因此相比于極端海況條件，波浪頻率成分偏小，共振運動成分更顯著。

5) CFD模擬結果相比部分勢流理論結果更加精細，可以捕捉到波浪的高階成分，并展示出流場可視化圖像。由于結構物的存在產生了波浪的輻射、繞射現象，極端海況下流場非線性更強。

6) 文中研究的風機基礎平臺在平面內的運動呈現出柔性特點，在平面外的運動近似于剛性，具有良好的運動性能，可有效避免平臺產生波浪頻率上的共振，作為風機發電基礎平臺可以達到穩定輸出的目的。