浮標系泊結構的參數優化匹配方法研究

方子帆，余紅昌，鄭小偉，何雪輝，葛旭甫，徐 浩，謝旭光

(1. 三峽大學 機械與動力學院，湖北 宜昌 443002；2. 三峽大學 新能源微電網湖北省協同創新中心，湖北 宜昌 443002；3. 三峽大學 水電機械設備設計與維護湖北省重點實驗室，湖北 宜昌 443002)

浮標系泊結構的參數優化匹配方法研究

方子帆1,2,3，余紅昌1，鄭小偉1，何雪輝1，葛旭甫1，徐 浩1，謝旭光1

(1. 三峽大學 機械與動力學院，湖北 宜昌 443002；2. 三峽大學 新能源微電網湖北省協同創新中心，湖北 宜昌 443002；3. 三峽大學 水電機械設備設計與維護湖北省重點實驗室，湖北 宜昌 443002)

針對淺海浮標系泊結構姿態控制問題，進行系統參數優化匹配計算方法研究。首先進行系泊結構力學分析，建立系統的力和力矩平衡方程，推導出系泊系統狀態的描述方程；采用系泊系統起錨角和檢測設備艙的傾斜角度為約束條件，取錨鏈型號和重物球質量作為設計變量，建立錨鏈的垂直高度相對誤差最小為目標的系統參數優化匹配的數學模型；用最小二乘法求解非線性方程組，進行分層變量優化計算。所提出的系泊系統結構參數優化計算方法可為浮標系泊系統結構設計提供理論依據。

浮標系泊結構；參數優化；優化模型；分層變量優化計算

0 引 言

在海況復雜多變情況下，進行浮標系泊結構姿態控制時面臨起錨角控制問題和水聲通訊設備艙與豎直方向的夾角限制問題，不然會造成錨被拖行和水聲通訊設備工作效果變差，需選取合適系泊設備，避免上述情況發生。其中確定錨鏈的型號、長度和重物球的質量成為關鍵，同時這些參數值可以對設備艙的傾斜角度和浮標吃水深度及游動區域姿態參數進行優化匹配，從而使系泊結構達到最佳控制姿態。

潘生兵等[1]直接通過船舶舾裝數計算查表來選取系泊設備。陳宇泉等[2]不僅通過舾裝數的計算，還利用環境載荷計算系泊力來選取系泊索，并論證2種不同計算方法的有效性。胡靈斌等[3]通過錨抓持力與外力平衡進行懸鏈方程求解來確定錨及錨索的尺寸。王丹等[4]推導懸鏈方程，結合工程船實況運用懸鏈線方程計算錨鏈參數, 并確定錨泊設備選型。MEA Yassir等[5]用懸鏈線方程隱式迭代分析一個多分量懸鏈線系泊系統。李曉平等[6]應用改進有限段模型對含彈性元件的小型系泊浮標系統進行動力學仿真，得出張力幅值頻率響應函數，為系泊系統設計提供指導意義。王磊等[7]基于單點系泊系統非線性動力學模型對浮標系統靜態姿態進行計算，方子帆等[8]建立低張力纜索有限元模型為系泊繩的運動分析與設計提供理論依據。但都沒有進行系泊結構參數優化匹配計算分析。

以起錨角和設備艙的傾斜角度為約束條件，取錨鏈型號和重物球質量作為設計變量，建立以錨鏈垂直高度相對誤差最小為目標的參數優化匹配數學模型。結合實例進行浮標吃水深度和游動區域及設備艙的傾斜角度多變量優化計算，得出系泊結構姿態控制的優化結果和調節規律，為系泊結構設計提供理論依據。

1 浮標系泊系統特性與分析

1.1 浮標系泊系統示意圖及其符號約定

浮標系泊系統由圓柱體浮標、鋼管、檢測設備艙、重物球、電焊錨鏈和特制的抗拖移錨組成，錨鏈選用無檔普通鏈環，常用型號及其參數見符號列表，鋼管共4節。根據實際工況要求起錨角α不能超過16°和裝有水聲通訊設備的設備艙與豎直方向的夾角β不能超過5°，確定合適錨鏈的型號、長度和重物球的質量能起到調節的作用。簡圖如圖1所示，符號約定如表1所示。

1.2 系泊系統力學分析

整個系泊系統在一定風載和水載共同作用下會達到一個平衡狀態，故可以簡化為靜力學分析[9-10]問題，出于力學分析的考慮，需對系泊系統的工作環境及其結構進行簡化和假設，模型假設如下：

1）假設海流為平面流，且無垂直方向分量，故簡化為二維問題[7]；

2）由于靜力學分析可不考慮波浪對整個系泊系統的影響；

表 1 浮標系泊系統符號約定Tab. 1 Sign convention of buoy mooring system

3）由于錨和錨鏈的幾何形狀尺寸未知，故假設在水載作用下阻力為0；

4）水下的錨、錨鏈、設備艙、鋼管的浮力不計。

通過靜力學分析，先整體后局部，依次建立X、Y軸力平衡和力矩平衡方程式組。

式中：G總=G浮標+G錨+G0，G0=G錨鏈+G+G艙+4×G鋼管；為風在水上浮標上的載荷，N；F水載=為水流力合力，N；S為物體在風速或水速法平面的投影面積，m2；N為支持力，N；Ff為錨所受靜摩擦力，N。

式中：N1為海床對錨鏈拖地部分的支持力，N；N =N1+N0，錨鏈若沒有著地N1=0。

式中：Fs0為設備艙所受水流力，N；T1=T1′。

式中：Fs1為鋼管所受水流力，N；T2=T2′。第3、第2、第1根鋼管受力分析同上。

式中：Fs5為浮標吃水部分所受水流力，N。

在已知水深H0下，根據確定上述參數值可以求得錨鏈的實際垂直高度h1。

1.3 系泊系統懸鏈方程

根據文獻[4]可知一般懸鏈方程的推導過程，任取懸鏈線上長度為ΔL的微元，其受力如圖8所示。

兩端為張力T和T+ΔT，其作用線與x軸成θ角，自重為qΔL，張力T+ΔT在x，y軸上的投影為連續函數，用泰勒公式將2個投影分別展開，取1階微量即可，可得在x軸上的投影為：

在y軸的投影：

圖 8 懸鏈線微元受力分析圖

Fig. 8 Stress analysis of catenary element

微元受力平衡方程式：

由式（19）、式（20）可得微分方程：

式中：a=R/q，R為錨鏈受到的水平分力，N；q為錨鏈單位長度在水中的重量，N/m。

邊界條件：

可得懸鏈的一般方程[11]：

式中：L為錨鏈的長度；x、y分別為錨鏈的水平和垂直分量。該懸鏈方程的垂直高度h2=y，將h2作為擬合值。

2 浮標系泊結構參數優化模型

2.1 設計變量

在浮標系泊結構姿態控制過程中，錨鏈型號和重物球的質量調節至關重要，可以使設備艙的傾斜角度和浮標吃水深度及游動區域參數進一步優化，且錨鏈型號和重物球一旦確定，整個系泊結構姿態在一定海況下也隨之確定，設計變量為錨鏈型號和重物球的質量。

2.2 約束條件

為了使錨鏈不被拖行，導致系泊系統移動，要求起錨角α不超過16°，則約束條件為：

為了不影響艙內水聲通訊設備的工作效果，設備艙的傾斜角度β不超過5°，則約束條件為：

2.3 目標函數

懸鏈方程求解的擬合值h2，與實際錨鏈垂直高度值h1對比，取相對誤差δ=|h1-h2|/h1×100%越小，與實際吻合程度越高。故以相對誤差δ最小為目標，即

其中，為了使實際情況相符，相對誤差控制為δ≤0.1%。

2.4 系泊結構參數優化匹配數學模型

綜上所述，單目標優化數學模型為

當相對誤差δ為0.1%時，若錨鏈出現拖地現象（h1＜ h2且a=0°），應另行求解。

3 模型求解

3.1 求解原理

在復雜海況下，風速、水速、水深可變因素較多造成分析困難，利用控制變量法分別加以分析。對于設備艙的傾斜角度和浮標吃水深度及游動區域多參數優化問題，需進行單變量優化計算分析來構成分層多變量優化計算分析方法，再通過最小二乘法fsolve函數逼近求出模型中非線性方程組的解和利用窮舉法區間搜索得出最優解。

3.2 求解過程

由于海況參數條件多為范圍或最值，在滿足功能要求下，應考慮極限條件下的系泊結構參數計算，求解流程如圖9所示。

由流程圖可得，當水深H0一定時，通過調節重物球質量，可求得設備艙的傾斜角度β在0～5°的變化下所對應的系泊系統狀態，若起錨角α超過16°，即調節重物球質量失效，需選用重型型號錨鏈進行調節。另外對于錨鏈拖地情況即h1＜h2，應將錨鏈實際高度h1作為懸鏈方程的垂直高度y重新代入原力與力矩平衡方程計算求解系泊結構參數。

4 工程應用

4.1 作業海況與系泊結構

由于潮汐因素的影響，布放海域的實測水深介于16～20 m之間。布放點的海水速度最大可達到1.5 m/s、風速最大可達到36 m/s。綜合考慮風力、水流力和水深情況下的系泊系統設計，分析不同情況下的設備艙、鋼管的傾斜角度、錨鏈形狀、浮標的吃水深度和游動區域。浮標系泊結構參數：圓柱體浮標的底面直徑2 m、高2 m，其質量1 000 kg；錨的質量為6 00 kg；鋼管每節長1 m、直徑50 mm，質量10 kg；裝有水聲通訊設備的設備艙長1 m、外徑30 cm，總質量100 kg。

4.2 重物球質量調節計算結果

設定水深H0=16 m、20 m兩個極限值，依此可得一個范圍解，海水最大速度為vsmax=1.5 m/s，最大風速vfmax=36 m/s，錨鏈的型號為II，其長度L=22.05 m情況下，系泊系統在設備艙的傾斜角度β變化下的設計情況如表2和圖10所示。

表 2 在H0=16 m下系泊系統狀態參數值Tab. 2 Mooring system state parameter values at H0=16 m

表 3 在H0=20 m下系泊系統狀態參數值Tab. 3 Mooring system state parameter values at H0=20 m

當水深深度H0一定時，設備艙的傾斜角度β越小（設備的工作效果越好），所加載重物球的質量m1越大，導致浮標下降，其吃水深度h增大，但風載減小，故起錨角α和游動區域半徑R減小。

4.3 錨鏈型號優選計算結果

在海水深度H0=20 m下，當設備艙的傾斜角度β≥2°時，調節重物球質量已經不能滿足α≤16°，應選用重型型號的錨鏈，優選III型號錨鏈結果如表3和圖11所示。

當水深H0為20 m時，起錨角α和設備艙的傾斜角度β會超過允許范圍，可以通過增大重物球質量或選取重型型號的錨鏈進行調節，且選取重型型號的錨鏈能夠調節重物球調節β值失效的那一部分，另通過上述H0=20 m，β=1°和2°情況下兩者調節系泊系統狀態所得α、H、R優化參數值的比較可知選取重型型號的錨鏈更好，其使浮標的吃水深度H和游動區域R相對更小，起錨角α也更小，從而使系泊系統更趨合理。

5 結 語

針對淺海浮標系泊結構參數優化匹配問題，以起錨角和設備艙的傾斜角度為約束條件，取錨鏈型號和重物球質量作為設計變量，建立錨鏈垂直高度相對誤差最小為目標的參數優化匹配數學模型，在不同水深情況下進行了浮標吃水深度和游動區域及設備艙的傾斜角度多變量優化計算，結論如下：

1）通過分層變量優化計算，當水深和錨鏈型號一定時，在設備艙不同的傾斜角度β下，可得其與α和R成正比，與m1和h成反比。

2）以系泊結構α，β，h，R為優化參數，通過調節重物球質量和錨鏈型號進行對比分析，可知選取更重錨鏈型號比增大重物球更合理，反之若從經濟上考慮增大重物球更劃算。

3）當系泊結構參數β或α超過允許范圍時，可優先通過增大重物球質量調節，若其調節失效，則選取重型型號的錨鏈來配合調節。

[1]潘生兵, 程慧勇, 李亮, 等. 船舶舾裝數計算及系泊設備參數確定[J]. 船舶標準化工程師, 2014, 3: 7-11 PAN Sheng-bing, CHENG Hui-yong, LI Liang, et al.Calculation of ship outfitting number and parameter determination of mooring equipment[J]. Ship Standardization Engineer, 2014, 3: 7-11

[2]陳宇泉. 海船系泊系統的設計與安裝研究[D]. 哈爾濱: 哈爾濱工程大學, 2013.

[3]胡靈斌, 唐軍. 懸鏈線方程的求解及其應用[J]. 船舶, 2004, 1:17-20.

[4]王丹, 劉家新. 一般狀態下懸鏈線方程的應用[J]. 船海工程,2007, 3: 26-28.WANG Dan, LIU Jia-xin. Application of catenary equation in general state[J]. Ship and Ocean Engineering, 2007, 3: 26-28.

[5]YASSIR M E A, KURIAN V J, HARAHAP I S H, et al.Parametric study on multi-component catenary mooring lines for offshore floating structures[J]. 2010.

[6]李曉平. 多體系統動力學建模方法及在水下纜索中的應用研究[D]. 天津: 天津大學, 2004.

[7]王磊. 單點系泊系統的動力學研究[D]. 青島: 中國海洋大學,2012.

[8]方子帆, 賀青松, 向兵飛, 等. 低張力纜索有限元模型及其應用[J]. 工程力學, 2013, 3: 445-450.

[9]金德賢. 李潤培. 半潛式鉆井平臺錨泊系統靜力分析[J]. 上海交通大學學報. 1981(2): 58-68.

[10]SENJANOVIC I, MRAVAK Z. Static analysis of catenary in mooring system[J]. Zbornik Radova Fsb, 1994.

[11]SUCH M, JIMENEZ-OCTAVIO J R, CARNICERO A, et al.An approach based on the catenary equation to deal with static analysis of three dimensional cable structures[J]. Engineering Structures, 2009, 31(31): 2162-2170.

Research on parameter optimization matching method of buoy mooring structure

FANG Zi-fan1,2,3, YU Hong-chang1, ZHENG Xiao-wei1, HE Xue-hui1, GE Xu-fu1, XU Hao1, XIE Xu-guang1
(1. College of Mechanical amp; Power Engineering, China Three Gorges University, Yichang 443002, China; 2. New Energy and Micro Grid Collaborative Innovation Center in Hubei Province, China Three Gorges University, Yichang 443002, China;3. Hubei Key Laboratory of Hydroelectric Machinery Design and Maintenance,China Three Gorges University Yichang 443002, China)

In order to solve the problem of attitude control of shallow sea buoy mooring system, the method of system parameter matching calculation is studied. First carry on the analysis of mooring structure mechanics, and establish the system of the force and moment equilibrium equations describing equations of mooring system state; the anchor angle of the mooring system and the tilt angle of detection equipment cabin as constraint conditions, take the quality of anchor type and heavy ball as the design variables, parameters optimization mathematical model, establish the system of vertical height the relative error of the target chain; using the least square method for solving nonlinear equations, stratified variable optimization. The proposed method can provide theoretical basis for the design of buoy mooring system.

the buoy mooring system；parameters optimization；optimization model；stratified variable optimization

U662

A

1672-7649(2017)11-0069-06

10.3404/j.issn.1672-7649.2017.11.013

2017-01-15；

2017-04-27

方子帆(1963-)，男，博士，教授，研究方向為機械系統動力學與控制。