基于神經(jīng)網(wǎng)絡(luò)和遺傳算法的系泊優(yōu)化設(shè)計(jì)

許小穎，周盼，王寬

1文華學(xué)院機(jī)械與電氣工程學(xué)部，湖北武漢430074

2中國(guó)船級(jí)社武漢分社，湖北武漢430074

基于神經(jīng)網(wǎng)絡(luò)和遺傳算法的系泊優(yōu)化設(shè)計(jì)

許小穎1，周盼1，王寬2

1文華學(xué)院機(jī)械與電氣工程學(xué)部，湖北武漢430074

2中國(guó)船級(jí)社武漢分社，湖北武漢430074

［目的］為使船舶維持作業(yè)位置的穩(wěn)定，需采用系泊系統(tǒng)來(lái)減小浮體的平動(dòng)。［方法］以一艘南海作業(yè)的鋪管船為例，通過(guò)優(yōu)化錨泊線的布置方式來(lái)最大限度地減小系泊狀態(tài)下錨鏈的平動(dòng)位移量，保障船舶作業(yè)安全。首先，以拋錨距離和錨鏈方位角作為正交試驗(yàn)的因素，得到不同的系泊布置方案，利用Moses軟件計(jì)算不同布置方案在不同浪向下的時(shí)域運(yùn)動(dòng)位移和錨鏈?zhǔn)芰η闆r。然后，將結(jié)果作為樣本，對(duì)BP神經(jīng)網(wǎng)絡(luò)進(jìn)行訓(xùn)練，實(shí)現(xiàn)BP網(wǎng)絡(luò)對(duì)Moses時(shí)域計(jì)算的仿真。最后，將拋錨距離和錨鏈方位角作為優(yōu)化變量，取不同浪向下的加權(quán)平動(dòng)位移為優(yōu)化目標(biāo)，并以BP神經(jīng)網(wǎng)絡(luò)的泛化能力代替Moses的時(shí)域計(jì)算，采用遺傳算法進(jìn)行優(yōu)化求解。［結(jié)果］結(jié)果表明，該鋪管船各個(gè)浪向下的平動(dòng)位移均有了顯著的減小，優(yōu)化效果明顯，［結(jié)論］可為海上浮式結(jié)構(gòu)物的系泊布置提供參考。

系泊優(yōu)化；BP神經(jīng)網(wǎng)絡(luò)；遺傳算法；Moses；時(shí)域分析

Abstract：［Objectives］In order to maintain the stability of the position of a ship,a mooring system is required to reduce the translational motion of floating structures.［Methods］Taking a pipe-laying vessel in the South China Sea as an example,it is possible to minimize the translational displacement of the anchor chain in the mooring state by optimizing the arrangement of the anchor line to ensure the safe operation of the ship.First,we can obtain several different layouts through orthogonal testing after selecting the azimuth and distance of the anchor chain as the test factors.We then calculate the different movements and force in time domain value at different wave direction angles for each layout using Moses.With the calculation results as samples,the BP neural network method achieves time domain simulation in Moses.After choosing the azimuth and distance of the anchor chain as the optimization variables,and with each wave-weighted translational displacement probability as the optimization objective,we find that the generalization capability of the BP neural network method can replace the time domain calculation of Moses.［Results］Using a genetic algorithm optimization solution,movement is significantly reduced at different wave direction angles.［Conclusions］This conclusion can provide a reference for the mooring arrangements of floating structures.

Key words：mooring optimization；BP neural network；genetic algorithm；Moses；time domain analysis

0 引 言

眾所周知，船舶或浮動(dòng)平臺(tái)要在海上保持一個(gè)固定位置，必須采用系泊系統(tǒng)。系泊結(jié)構(gòu)運(yùn)動(dòng)位移過(guò)大會(huì)影響其正常作業(yè)，而錨泊線布置不合理則容易致其斷裂。因此，減小系泊狀態(tài)下的平動(dòng)位移，合理布置錨泊線就顯得尤為重要。余龍和譚家華［1］提出多目標(biāo)函數(shù)的多成分錨泊線優(yōu)化模型，通過(guò)遺傳算法并結(jié)合圖譜設(shè)計(jì)的方式，完成了錨泊系統(tǒng)設(shè)計(jì)實(shí)例。馬小劍等［2］對(duì)風(fēng)、流作用下碼頭的系泊船舶，以及其錨鏈張力和平移運(yùn)動(dòng)量進(jìn)行了研究。Ali和Inoue［3］則假定每一根系泊纜完全相同，采用遺傳算法，以橫蕩和縱蕩的合運(yùn)動(dòng)最小值為目標(biāo)函數(shù)，對(duì)錨鏈與船艏的相對(duì)夾角及錨鏈長(zhǎng)度進(jìn)行了優(yōu)化。嚴(yán)傳續(xù)等［4］分析了某鋪管船錨泊定位系統(tǒng)的結(jié)構(gòu)及其受力情況，通過(guò)研究錨鏈懸鏈狀態(tài)的特性，提出錨泊定位系統(tǒng)的二次規(guī)劃優(yōu)化算法，并用面向?qū)ο蟮腣C++工具開(kāi)發(fā)了錨泊定位系統(tǒng)軟件。上述學(xué)者的研究均取得了顯著成果，但也存在著一些問(wèn)題：首先，對(duì)于系泊系統(tǒng)的研究大多還停留在對(duì)簡(jiǎn)單模型的仿真模擬階段，未涉及實(shí)船的系泊優(yōu)化問(wèn)題；其次，多數(shù)研究還處于針對(duì)現(xiàn)有布置的錨鏈?zhǔn)芰εc平動(dòng)位移的計(jì)算上，并不能做到對(duì)不同錨泊線布局情況下結(jié)果的預(yù)測(cè)；再次，基于水動(dòng)力軟件模擬計(jì)算的結(jié)果雖然較合理，但耗時(shí)巨大，并不適用于錨泊線多種布局方式的推廣研究。

針對(duì)上述問(wèn)題，本文將從工程實(shí)際出發(fā)，以一艘在南海工作的鋪管船為研究對(duì)象，將Moses軟件的計(jì)算結(jié)果作為樣本，與目前流行的BP神經(jīng)網(wǎng)絡(luò)和遺傳算法進(jìn)行結(jié)合，考慮環(huán)境因素的影響，以錨鏈的方位角和錨位置作為優(yōu)化變量，對(duì)其布置進(jìn)行優(yōu)化設(shè)計(jì)。主要研究?jī)?nèi)容包括：1）根據(jù)鋪管船的圖紙資料，應(yīng)用Moses軟件建模進(jìn)行頻域分析，得到其水動(dòng)力系數(shù)，然后通過(guò)卷積積分轉(zhuǎn)化為時(shí)域下的計(jì)算結(jié)果，考慮風(fēng)、流的作用，計(jì)及船體與錨鏈的耦合影響，采用基于Newmark積分的時(shí)域分析方法計(jì)算常規(guī)系泊布置下的船體運(yùn)動(dòng)及錨鏈?zhǔn)芰Α?）依據(jù)正交試驗(yàn)原理，錨鏈的方位角和拋錨距離各取4個(gè)水平，得到32組系泊布置方案，通過(guò)編制批處理文件，得到32組不同布置下的浮體運(yùn)動(dòng)與纜繩受力值，并將其作為BP神經(jīng)網(wǎng)絡(luò)訓(xùn)練樣本進(jìn)行仿真模擬。3）以拋錨距離和錨鏈方位角為優(yōu)化變量、以不同浪向下浮體平動(dòng)位移的最小值為目標(biāo)函數(shù)，采用遺傳算法進(jìn)行優(yōu)化求解，以使該位移值顯著減小，優(yōu)化效果明顯。研究所得結(jié)果用于為海上浮式結(jié)構(gòu)物的系泊布置提供參考。

1 神經(jīng)網(wǎng)絡(luò)和遺傳算法

1.1 神經(jīng)網(wǎng)絡(luò)

神經(jīng)網(wǎng)絡(luò)具有強(qiáng)大的逼近能力和推廣作用［5］，其通過(guò)協(xié)調(diào)多種輸入關(guān)系來(lái)處理復(fù)雜的非線性問(wèn)題，一般由輸入層、隱含層和輸出層組成。隱含層的數(shù)目可自定，本研究選用3層網(wǎng)絡(luò)［6］。各層之間通過(guò)權(quán)值相連，每層均含有多個(gè)神經(jīng)元，每個(gè)神經(jīng)元自帶閾值，這樣輸入和輸出之間即為映射關(guān)系，通過(guò)權(quán)（閾）值來(lái)表現(xiàn)，如圖1所示。通常采用sigmoid型作用函數(shù)，其輸出范圍為0～1之間的小數(shù)，采用正切sigmoid（tansig）為隱層傳遞函數(shù)，采用線性purelin為輸出層傳遞函數(shù)。訓(xùn)練成功的三層網(wǎng)絡(luò)模型可用圖2和式（1）來(lái)表達(dá)。圖2中a1為中間層變量；a2為輸出變量；w1為輸入層到中間層的權(quán)值；w2為中間層到輸出層的權(quán)值；b1為中間層的閾值；b2為輸出層的閾值；n1為中間層神經(jīng)元個(gè)數(shù)；n2為輸出層神經(jīng)元個(gè)數(shù)；p為輸入變量。

式中：Oi為輸出變量；ξn為輸入變量；Wkn，Wjk，Wij為各層神經(jīng)元的權(quán)值；b1k，b2j，b3i為各層神經(jīng)元的閾值；i，j，k分別為神經(jīng)網(wǎng)絡(luò)訓(xùn)練過(guò)程中輸入層、中間隱層和輸出層的神經(jīng)元個(gè)數(shù)。

圖1 BP神經(jīng)網(wǎng)絡(luò)Fig.1 Structure of BP neural network

圖2 BP神經(jīng)網(wǎng)絡(luò)的計(jì)算過(guò)程Fig.2 Calculational process of BP neural network

網(wǎng)絡(luò)的誤差信號(hào)［7］指的是實(shí)際輸出與期望輸出之間的差值，由輸出端逐層向前傳播。在傳播過(guò)程中，網(wǎng)絡(luò)的權(quán)值隨誤差反饋進(jìn)行調(diào)節(jié)，經(jīng)過(guò)修正后制定其相應(yīng)的誤差范圍。本文采用的誤差信號(hào)為

樣本的輸入和輸出采用歸一化處理，如式（3）所示。

式中：xi為輸入（出）數(shù)據(jù)；xmin為數(shù)據(jù)變化的最小值；xmax為數(shù)據(jù)變化的最大值。

全局誤差［8］選用均方誤差表示：

式中：m為輸出節(jié)點(diǎn)數(shù)；p為訓(xùn)練樣本數(shù)；為期望輸出；yij為實(shí)際輸出。

MATLAB軟件自帶的神經(jīng)網(wǎng)絡(luò)工具箱包含多種訓(xùn)練算法［9］，如增加動(dòng)量法、自適應(yīng)學(xué)習(xí)率法、彈性BP算法、L-M算法、動(dòng)量及自適應(yīng)學(xué)習(xí)率法等。經(jīng)過(guò)試算，本文從訓(xùn)練時(shí)間和誤差兩方面考慮，選用了有動(dòng)量梯度下降的BP算法訓(xùn)練函數(shù)，通過(guò)改進(jìn)參數(shù)設(shè)置，有效避免了其收斂速度過(guò)慢、局部極值而得不到全局最優(yōu)等缺點(diǎn)。

1.2 遺傳算法

該方法的操作對(duì)象為一群二進(jìn)制串，利用初始種群，選擇滿足適應(yīng)度函數(shù)的個(gè)體，通過(guò)雜交和變異產(chǎn)生下代種群，直至滿足期望條件，其過(guò)程如圖3所示。

圖3 遺傳算法運(yùn)算過(guò)程圖Fig.3 Calculational process of genetic algorithm

MATLAB軟件的遺傳算法工具箱［10］包括各類子函數(shù)，如種群初始化函數(shù)、適應(yīng)度函數(shù)、選擇函數(shù)、交叉算子函數(shù)和變異算子函數(shù)。

對(duì)于種群初始化函數(shù)，選用crtrp函數(shù)，因其二進(jìn)制編（解）碼以及交叉、變異的實(shí)現(xiàn)較為方便。

式中：numind為種群個(gè)體數(shù)；lenind為個(gè)體長(zhǎng)度。

適應(yīng)度函數(shù)主要用于評(píng)價(jià)個(gè)體優(yōu)劣，適應(yīng)度值越高，個(gè)體遺傳下一代的概率就越高。因ranking函數(shù)是基于排序的適應(yīng)度分配，支持非線性評(píng)估，故選此函數(shù)，如式（6）所示。

式中，ObjV值越大，分配的適應(yīng)度值越小。

本文選用高級(jí)入口函數(shù)select對(duì)個(gè)體進(jìn)行優(yōu)勝劣汰，以方便多種群接口。

式中：SEL_F為選擇函數(shù)名；Chrom代表種群；FitnV包含種群Chrom中個(gè)體的適應(yīng)度值；GGAP為可選參數(shù)，指出代溝，本文取為0.9。

僅有選擇操作，并不能產(chǎn)生新的個(gè)體，故需要有交叉操作。本文選用高級(jí)入口函數(shù)recombine，因其對(duì)所有交叉操作可提供多子群支持入口。

式中：REC_F為包含了低級(jí)重組函數(shù)名的字符串；RecOpt為交叉概率的可選數(shù)值，本文選用單點(diǎn)交叉，概率值取為0.7。

選擇和交叉運(yùn)算有可能會(huì)造成某些信息丟失，而變異算子函數(shù)主要用于保證算法的有效性，故采用mut函數(shù)作為變異算子。

2 算例分析

2.1 研究背景

本文以國(guó)內(nèi)某公司的一艘鋪管船為例，計(jì)算系泊下時(shí)域運(yùn)動(dòng)和錨鏈的受力。工作錨鏈采用鍍鋅鋼絲繩，作業(yè)地點(diǎn)為南海某海區(qū)，船型及錨鏈參數(shù)、環(huán)境條件如表1所示，錨鏈布置情況如表2所示。其中，方位角表示錨鏈與船體坐標(biāo)系中X軸的夾角，指向Y軸為正。

表1 船型、錨鏈和環(huán)境主要參數(shù)Table 1 Main parameters of ship，anchor chain and environment

初始系泊分析應(yīng)用Moses軟件中的時(shí)域模塊［11］模擬不同位置處系泊的剛度變化，可進(jìn)行橫搖阻尼修正，具備一定的可信度。根據(jù)船體型線圖，建立船體各站的形狀，形成型表面。該船的型線較為規(guī)則，網(wǎng)格最大邊長(zhǎng)設(shè)為12 m，經(jīng)程序自動(dòng)劃分后共形成678個(gè)面元。水動(dòng)力模型如圖4所示，整體系泊外觀圖如圖5所示。

表2 錨鏈布置Table 2 Layout of anchor chain

圖4 水動(dòng)力分析模型Fig.4 Model of hydrodynamic analysis

圖5 系泊外觀圖Fig.5 Layout of mooring

通過(guò)Moses軟件的頻域分析模塊，可得到水動(dòng)力計(jì)算結(jié)果：附加質(zhì)量和附加質(zhì)量慣性矩系數(shù)、阻尼系數(shù)、幅值響應(yīng)算子（Response Amplitude Operator，RAO）。其中，附加質(zhì)量（慣性矩）系數(shù)為附加質(zhì)量（慣性矩）與排水量的比值。頻域計(jì)算所選的海浪譜為JONSWAP譜［8］，波浪周期范圍為4～40 s，浪向角為0°，即順浪。

對(duì)頻率響應(yīng)函數(shù)進(jìn)行傅里葉變換可得到脈沖響應(yīng)函數(shù)［12］，海況則利用JONSWAP譜形成的偽隨機(jī)波來(lái)模擬，進(jìn)而通過(guò)卷積積分得到時(shí)域運(yùn)動(dòng)響應(yīng)。計(jì)算時(shí)長(zhǎng)設(shè)為4 000 s，時(shí)間步長(zhǎng)設(shè)為0.2 s，不計(jì)前500 s的計(jì)算結(jié)果。因系泊錨鏈主要影響船體（重心位置）平動(dòng)，故僅列出了不同預(yù)張力下船體發(fā)生橫蕩和縱蕩運(yùn)動(dòng)時(shí)的時(shí)域計(jì)算結(jié)果，如圖6和表3所示，錨鏈?zhǔn)芰η闆r如表4所示?？梢?jiàn)在3種情況下，錨鏈的受力均滿足強(qiáng)度標(biāo)準(zhǔn)，Moses軟件計(jì)算結(jié)果可信。在實(shí)際工程作業(yè)中，由于錨鏈長(zhǎng)度的調(diào)節(jié)可通過(guò)改變預(yù)張力來(lái)實(shí)現(xiàn)，故表中的30%指的是錨鏈預(yù)張力取為錨機(jī)工作能力的30%（30 t）［13］，其余類推。

圖6 不同預(yù)張力下縱蕩運(yùn)動(dòng)Fig.6 Surge motion under different pre-tension

表3 不同預(yù)張力下的時(shí)域運(yùn)動(dòng)統(tǒng)計(jì)Table 3 Time-domain motion under different pre-tension

表4 錨鏈最大受力值Table 4 Maximum tension of anchor chain

2.2 網(wǎng)絡(luò)計(jì)算模型

本文以錨鏈方位角和拋錨距離作為正交試驗(yàn)的因素，總個(gè)數(shù)為16。在工程實(shí)際中，因錨鏈多沿船長(zhǎng)方向?qū)ΨQ布置，故選取船右舷一側(cè)的錨鏈即可，因素個(gè)數(shù)可減為8，極大地減小了工作量。對(duì)于每個(gè)因素，均取4個(gè)水平進(jìn)行正交試驗(yàn)設(shè)計(jì)。根據(jù)工程經(jīng)驗(yàn)，分別選取拋錨距離為850，1 000，1 500和2 000 m。為使各錨鏈能按逆時(shí)針?lè)较蛞来闻帕卸唤诲e(cuò)，將錨鏈S2的方位角設(shè)定為95°，105°，120°和135°；錨鏈S4的方位角設(shè)為5°，15°，30°和45°；錨鏈S1和S3的方位角則根據(jù)與S2和S4的夾角來(lái)確定，其中夾角的4個(gè)水平均為10°，20°，30°，和40°。采用正交表L32（49）來(lái)安排正交試驗(yàn)，具體如表5所示。

表5 0°浪向下的時(shí)域計(jì)算結(jié)果Table 5 Time domain calculation results in 0°wave direction

鑒于實(shí)驗(yàn)組數(shù)較多，且一組實(shí)驗(yàn)計(jì)算約需1 h，故編寫(xiě)了批處理文件，可計(jì)算得到時(shí)域下鋪管船在0°，45°，90°，135°和180°浪向角下的平動(dòng)位移以及錨鏈?zhǔn)芰Φ挠?jì)算結(jié)果。受篇幅所限，這里僅給出了0°浪向下的結(jié)果（表5），其余角度下的不再贅述。

根據(jù)經(jīng)驗(yàn)，將縱蕩和橫蕩運(yùn)動(dòng)仿真的輸出節(jié)點(diǎn)設(shè)為2，錨鏈?zhǔn)芰Ψ抡娴妮敵龉?jié)點(diǎn)設(shè)為8，兩者的輸入節(jié)點(diǎn)均設(shè)為8；隱含節(jié)點(diǎn)數(shù)通過(guò)循環(huán)計(jì)算設(shè)置，從5～20中選擇全局誤差最小的數(shù)值作為最終節(jié)點(diǎn)數(shù)。將表5中的1～30組正交試驗(yàn)取為訓(xùn)練樣本，31和32組正交試驗(yàn)取為測(cè)試樣本，驗(yàn)證BP神經(jīng)網(wǎng)絡(luò)的泛化能力。

受篇幅所限，這里僅給出0°浪向下浮體的運(yùn)動(dòng)情況。由表5中橫蕩運(yùn)動(dòng)位移多接近0 m，最大值僅為0.04 m，故僅對(duì)縱蕩運(yùn)動(dòng)進(jìn)行仿真，網(wǎng)絡(luò)訓(xùn)練過(guò)程如圖7所示，錨鏈?zhǔn)芰Φ腂P網(wǎng)絡(luò)訓(xùn)練過(guò)程如圖8所示。前者經(jīng)過(guò)940次訓(xùn)練，得到隱含節(jié)點(diǎn)個(gè)數(shù)為14；后者經(jīng)過(guò)2 376次訓(xùn)練，得到隱含節(jié)點(diǎn)個(gè)數(shù)為14，兩者的誤差均小于設(shè)定的目標(biāo)全局誤差0.01。用第31和32組正交試驗(yàn)對(duì)訓(xùn)練好的神經(jīng)網(wǎng)絡(luò)進(jìn)行測(cè)試，結(jié)果如表6所示。

圖7 0°浪向下縱蕩運(yùn)動(dòng)仿真的網(wǎng)絡(luò)訓(xùn)練過(guò)程Fig.7 Training procedure of surge motion in 0°wave direction

圖8 0°浪向下錨鏈?zhǔn)芰Φ木W(wǎng)絡(luò)訓(xùn)練過(guò)程Fig.8 Training procedure of anchor chain tension in 0°wave direction

表6 0°浪向下網(wǎng)絡(luò)的相對(duì)誤差Table 6 Relative errors of simulation training in 0°wave direction

由測(cè)試結(jié)果可知，對(duì)5個(gè)浪向下的時(shí)域運(yùn)動(dòng)和錨鏈?zhǔn)芰Φ姆抡嫱耆梢詽M足工程需要，可見(jiàn)依靠BP神經(jīng)網(wǎng)絡(luò)的泛化能力代替Moses進(jìn)行時(shí)域計(jì)算可行。

根據(jù)規(guī)范對(duì)系泊系統(tǒng)的規(guī)定，結(jié)合考慮該船在實(shí)際作業(yè)海域各個(gè)浪向的發(fā)生概率（表7），采用遺傳算法進(jìn)行了優(yōu)化設(shè)計(jì)。以拋錨距離和錨鏈方位角為設(shè)計(jì)變量，取浮體縱蕩和橫蕩運(yùn)動(dòng)位移最小為優(yōu)化目標(biāo)，取上述BP網(wǎng)絡(luò)得到的仿真訓(xùn)練值并以每根錨鏈的受力要求作為約束條件。其數(shù)學(xué)表達(dá)如下：

表7 不同浪向的發(fā)生概率Table 7 Occurrence probability of different wave directions

式中：α為浪向發(fā)生概率；i為浪向個(gè)數(shù)，i=5；Δx，Δy分別為縱蕩、橫蕩運(yùn)動(dòng)的位移；X為設(shè)計(jì)變量向量；m為變量個(gè)數(shù)，m=8；Xmin，Xmax分別為變量邊界的下限和上限值，拋錨距離取850～ 2 000 m；fk為每根錨鏈的受力；[f]為錨鏈破斷力除以安全系數(shù)后的值；k為錨鏈根數(shù)，k=8。

運(yùn)用MATLAB軟件的遺傳算法工具箱進(jìn)行計(jì)算。因初始種群的隨機(jī)性，使得每次所得到的最優(yōu)解均不同，故設(shè)定循環(huán)計(jì)算500次，取其中的最優(yōu)解，如表8所示。按優(yōu)化后的方案布置錨鏈，再次應(yīng)用Moses軟件進(jìn)行時(shí)域計(jì)算，并將優(yōu)化前、后的結(jié)果進(jìn)行對(duì)比，如圖9所示。

表8 優(yōu)化后的錨鏈布置Table 8 Layout of anchor chain after optimization

圖9 優(yōu)化前、后的最大平動(dòng)位移對(duì)比Fig.9 Comparisons of maximum displacement before and after optimization

2.3 小 結(jié)

1）由Moses軟件計(jì)算結(jié)果可知，提高錨鏈的預(yù)張力，浮體的縱蕩位移會(huì)逐漸減小，故可通過(guò)增大預(yù)張力來(lái)改善系泊浮體結(jié)構(gòu)的運(yùn)動(dòng)。

2）由BP神經(jīng)網(wǎng)絡(luò)對(duì)5個(gè)浪向下的時(shí)域運(yùn)動(dòng)和錨鏈?zhǔn)芰Φ姆抡?，發(fā)現(xiàn)誤差均不超過(guò)15%，可滿足工程需要。因此，用神經(jīng)網(wǎng)絡(luò)的泛化能力代替Moses軟件的時(shí)域計(jì)算被認(rèn)為是可行的。

3）錨鏈經(jīng)遺傳算法優(yōu)化布置后，除在0°浪向下結(jié)果有所增大外，其余浪向下的平動(dòng)位移均不同程度地有所減小。其中，90°浪向，也就是橫浪下的運(yùn)動(dòng)改變最為顯著。

4）本文建立的基于BP神經(jīng)網(wǎng)絡(luò)和遺傳算法的系泊系統(tǒng)優(yōu)化方法可以為海上浮式結(jié)構(gòu)物的系泊布置提供參考。

3 結(jié) 語(yǔ)

本文從工程實(shí)際出發(fā)，以一艘在南海工作的鋪管船為研究對(duì)象，考慮環(huán)境因素的影響，以錨鏈的方位角和錨位置作為優(yōu)化變量，基于BP神經(jīng)網(wǎng)絡(luò)和遺傳算法對(duì)其布置進(jìn)行了優(yōu)化設(shè)計(jì)。在此過(guò)程中，還與Moses軟件的計(jì)算結(jié)果進(jìn)行了對(duì)比，發(fā)現(xiàn)BP神經(jīng)網(wǎng)絡(luò)在滿足精度要求的基礎(chǔ)上更節(jié)省時(shí)間，無(wú)論是對(duì)樣本的模擬，還是其優(yōu)越的預(yù)測(cè)能力均得到了充分的體現(xiàn)。限于既定的硬件條件以及筆者知識(shí)儲(chǔ)備的欠缺，本文的結(jié)果仍具有一定的局限性，還有待于學(xué)者們進(jìn)行進(jìn)一步的探討。

［1］余龍，譚家華.基于準(zhǔn)靜定方法的多成分錨泊線優(yōu)化［J］.海洋工程，2005，23（1）：69-73.YU L，TAN J H.Research on optimum multi-component mooring lines based on catenary equation［J］.The Ocean Engineering，2005，23（1）：69-73（in Chinese）.

［2］馬小劍，孫昭晨，張志明，等.風(fēng)流作用下碼頭系泊船舶纜繩張力及運(yùn)動(dòng)量研究［J］.水道港口，2010，31（3）：165-169.MA X J，SUN Z C，ZHANG Z M，et al.A numerical study on mooring line forces and motions of a moored ship under combined action of wind and current［J］.Journal of Waterway and Harbor，2010，31（3）：165-169（in Chinese）.

［3］ALI M T，INOUE Y.On hydrodynamic interaction between two rectangular barges floating side-by-side in regular waves［C］//Proceedings of the International Conference on Mechanical Engineering.Dhaka，Bangladesh：ICME，2005.

［4］嚴(yán)傳續(xù)，錢宏，項(xiàng)軍毅，等.鋪管船錨泊定位系統(tǒng)優(yōu)化設(shè)計(jì)研究［J］.中國(guó)造船，2010，51（1）：84-94.YAN C X，QIAN H，XIANG J Y，et al.Study on design optimization of mooring positioning system for pipelaying vessel［J］.Shipbuilding of China，2010，51（1）：84-94（in Chinese）.

［5］蘇曉宇，金鴻章，胡曉東.自動(dòng)錨泊定位系統(tǒng)的張力混合優(yōu)化算法［J］.電機(jī)與控制學(xué)報(bào)，2015，19（4）：100-105.SU X Y，JIN H Z，HU X D.Tension hybrid optimization strategy of automatic mooring positioning system［J］.Electric Machines and Control，2015，19（4）：100-105（in Chinese）.

［6］毛建平.基于神經(jīng)網(wǎng)絡(luò)的橋梁結(jié)構(gòu)靜力有限元模型修正［D］.長(zhǎng)春：吉林大學(xué)，2011.

［7］朱文龍.基于遺傳算法的BP神經(jīng)網(wǎng)絡(luò)在多目標(biāo)優(yōu)化中的應(yīng)用研究［D］.哈爾濱：哈爾濱理工大學(xué)，2009.

［8］李揚(yáng).基于遺傳算法的BP神經(jīng)網(wǎng)絡(luò)在網(wǎng)絡(luò)數(shù)據(jù)抓取中的研究［D］.青島：中國(guó)海洋大學(xué)，2010.

［9］許小穎.基于BP神經(jīng)網(wǎng)絡(luò)的三體船布局優(yōu)化研究［D］.天津：天津大學(xué)，2013.

［10］張偉.系泊狀態(tài)下柔性連接多浮體耦合運(yùn)動(dòng)研究［D］.天津：天津大學(xué)，2008.

［11］吳廣懷，沈慶，陳徐均，等.浮體間距對(duì)多浮體系統(tǒng)水動(dòng)力系數(shù)的影響［J］.海洋工程，2003，21（4）：29-34.WU G H，SHEN Q，CHEN X J，et al.Influence of the distance between floating bodies on hydrodynamic coefficients of floating multi-body system［J］.The Ocean Engineering，2003，21（4）：29-34（in Chinese）.

［12］WANG J F，LI J D，CAI X G，et al.Hydrodynamic interactions between two bodies in waves in 3D time domain［J］.Journal of Marine Science and Application，2005，4（1）：15-20.

［13］OHKUSU M.On the heaving motion of two circular cylinders on the surface of a fluid［R］.Reports of Research Institute for Applied Mechanics，1969，17（58）：167-185.

Mooring optimization design based on neural network and genetic algorithm

XU Xiaoying1，ZHOU Pan1，WANG Kuan2
1 Department of Mechanical and Electrical Engineering，Wenhua College，Wuhan 430074，China
2 Wuhan Branch of China Classification Society，Wuhan 430074，China

U674.38+1

A

10.3969/j.issn.1673-3185.2017.05.012

2017-01-06< class="emphasis_bold">網(wǎng)絡(luò)出版時(shí)間：

時(shí)間：2017-9-26 10:38

文華學(xué)院青年基金資助項(xiàng)目（J02e0540211）

許小穎，女，1987年生，碩士，講師。研究方向：船舶水動(dòng)力學(xué)。E-mail：xuxiaoying52616@163.com

周盼（通信作者），女，1990年生，碩士，助教。研究方向：結(jié)構(gòu)強(qiáng)度分析與設(shè)計(jì)

http://kns.cnki.net/kcms/detail/42.1755.TJ.20170926.1038.014.html期刊網(wǎng)址：www.ship-research.com

許小穎，周盼，王寬.基于神經(jīng)網(wǎng)絡(luò)和遺傳算法的系泊優(yōu)化設(shè)計(jì)［J］.中國(guó)艦船研究，2017，12（5）：97-103.

XU X Y，ZHOU P，WANG K.Mooring optimization design based on neural network and genetic algorithm［J］.Chinese Journal of Ship Research，2017，12（5）：97-103.