（海軍航空大學(xué)航空基礎(chǔ)學(xué)院 煙臺(tái) 264000）

1 引言

為了保證海洋浮標(biāo)平臺(tái)設(shè)備的可靠性和數(shù)據(jù)測(cè)量準(zhǔn)確性，海洋浮標(biāo)平臺(tái)在海上漂浮需要一個(gè)穩(wěn)定的狀態(tài)，如何在保持安全可靠的情況下兼顧經(jīng)濟(jì)性，是當(dāng)前浮標(biāo)系泊需要解決的問題，系泊系統(tǒng)的設(shè)計(jì)具有非常重要的意義。目前，常見的方法都是根據(jù)遺傳算法或神經(jīng)網(wǎng)絡(luò)方法來計(jì)算布線位置。主要通過優(yōu)化錨泊線的布置方式來最大限度地減小系泊狀態(tài)下錨鏈的平動(dòng)位移量［1］。以錨鏈拋錨距離和方位角作為變量，通過正交試驗(yàn)選取了最優(yōu)的系泊布置方案［2～3］。這些研究主要是針對(duì)船舶系泊纜索位置的優(yōu)化，本文以纜索的物理性質(zhì)的優(yōu)化為研究的重點(diǎn)，對(duì)具有對(duì)稱的特點(diǎn)浮標(biāo)平臺(tái)進(jìn)行系泊方案設(shè)計(jì)與優(yōu)化。

2 纜索優(yōu)化設(shè)計(jì)

通過采用有限元分析對(duì)纜索截面參數(shù)以及長(zhǎng)度來進(jìn)行優(yōu)化設(shè)計(jì)［4～6］，分析的過程中插入了程序語言對(duì)不同纜索參數(shù)的選用進(jìn)行了自動(dòng)化處理，從而進(jìn)行批量的水動(dòng)力計(jì)算和后處理。在該段處理中，將待選的各種截面參數(shù)化，程序可以直接調(diào)用對(duì)應(yīng)截面參數(shù)進(jìn)行計(jì)算。

2.1 實(shí)驗(yàn)變量設(shè)計(jì)

對(duì)纜索優(yōu)化的目的是盡可能減小纜索總重量，且保證浮臺(tái)的搖蕩幅值和纜索張力峰值也盡可能?。?～8］。優(yōu)化分析的輸入變量有纜索總長(zhǎng)度、單位長(zhǎng)度纜索重量、軸向剛度以及等效直徑［1］。纜索軸向剛度和等效直徑是與單位長(zhǎng)度纜索重量緊密相關(guān)的，即相同的材料，纜索越粗其重量與剛度就越大。根據(jù)纜索截面積與其重量與剛度的正比例關(guān)系，可在纜索重量、軸向剛度和等效直徑三者中選擇纜索重量作為設(shè)計(jì)變量，另外兩個(gè)變量作為其表達(dá)式關(guān)系。那么實(shí)際的輸入變量就是纜索總長(zhǎng)度與纜索重量?jī)蓚€(gè)不相關(guān)的變量?，F(xiàn)有的方案下，纜索的總長(zhǎng)為55m，纜索重量為150kg/m。鑒于纜索的總長(zhǎng)較短，如繼續(xù)減短，則有可能增大浮臺(tái)的初始拉力，因此纜索的總長(zhǎng)取更大的值，而纜索的重量為經(jīng)濟(jì)性考慮則取更小的值，各輸入變量的設(shè)計(jì)區(qū)間見表1。

表1 各輸入變量的設(shè)計(jì)區(qū)間

輸出變量指定為浮臺(tái)垂蕩運(yùn)動(dòng)的谷值與峰值，浮臺(tái)搖動(dòng)角度的谷值與峰值以及纜索張力的峰值。

2.2 實(shí)驗(yàn)設(shè)計(jì)

表2 實(shí)驗(yàn)設(shè)計(jì)中的輸入變量與輸出變量

本文通過設(shè)計(jì)實(shí)驗(yàn)［9］建立響應(yīng)面模型。因?yàn)閷?shí)際的設(shè)計(jì)變量只有兩個(gè)變量：纜索重量與纜索長(zhǎng)度，那么在實(shí)驗(yàn)設(shè)計(jì)中可將纜索重量變化8次，纜索長(zhǎng)度變化6次，來進(jìn)行參數(shù)化計(jì)算。纜索重量值分別為150、135、120、105、90、75、60、45，量綱kg/m；纜索長(zhǎng)度分別為55、56、57、58、59、60，量綱為m；一共進(jìn)行6×8=48次計(jì)算，由于篇幅所限表格中只列出部分計(jì)算結(jié)果。

2.3 響應(yīng)面分析

實(shí)驗(yàn)設(shè)計(jì)的48次離散樣本采用ANSYS的響應(yīng)面程序［10～11］，對(duì)離散的樣本進(jìn)行連續(xù)插值計(jì)算，生成了1000組數(shù)據(jù)，隨后將各個(gè)數(shù)據(jù)生成三維坐標(biāo)系的統(tǒng)計(jì)圖，形成了響應(yīng)面，通過對(duì)響應(yīng)面的分析便能夠得到各參數(shù)之間的關(guān)系。通過分析響應(yīng)面［12］，以此達(dá)到優(yōu)化的目的。

圖1 浮臺(tái)垂蕩幅值的三維響應(yīng)面

圖2 浮臺(tái)搖動(dòng)幅值的三維響應(yīng)面

圖3 纜索張力峰值的三維響應(yīng)面

其中：X軸為單位長(zhǎng)度纜索質(zhì)量；Y軸為纜索長(zhǎng)度；Z軸分別為浮臺(tái)的垂蕩幅值、浮臺(tái)的搖動(dòng)幅值和纜索的張力峰值。由圖1可得，浮臺(tái)的垂蕩振幅與繩索的長(zhǎng)度無關(guān)，而與繩索的質(zhì)量有關(guān)，質(zhì)量越小垂蕩振幅越小。

由圖2可得，浮臺(tái)的搖動(dòng)幅值與繩索的質(zhì)量關(guān)系較小，而與繩索的長(zhǎng)度關(guān)系較大，繩索較大或較小時(shí)搖動(dòng)幅值較小。

由圖3可得，纜索張力峰值與纜繩的長(zhǎng)度關(guān)系較小，與繩索的重量關(guān)系較大，繩索的重量越大，纜繩的張力峰值越大。

從變量的敏感度來看，浮臺(tái)的垂蕩運(yùn)動(dòng)和纜索的張力與纜索的重量相關(guān)性較大；而浮臺(tái)的搖動(dòng)角度與纜索長(zhǎng)度相關(guān)性較大。

2.4 目標(biāo)驅(qū)動(dòng)優(yōu)化

指定各變量的優(yōu)化目標(biāo)見表3。

表3 輸出變量與優(yōu)化目標(biāo)

圖4 浮臺(tái)垂蕩幅值的可行樣本分布圖

圖5 浮臺(tái)搖動(dòng)幅值的可行樣本分布圖

通過仿真實(shí)驗(yàn)在響應(yīng)面得到的大量連續(xù)樣本中進(jìn)行篩選，圖4～圖6分別顯示浮臺(tái)的垂蕩運(yùn)動(dòng)幅值、浮臺(tái)的搖動(dòng)角度幅值和纜索的張力峰值在滿足優(yōu)化目標(biāo)條件下的可行樣本分布圖。系統(tǒng)根據(jù)垂蕩幅值、搖動(dòng)角度幅值以及纜索張力峰值的優(yōu)化目標(biāo)對(duì)這些可行的樣本進(jìn)行最優(yōu)解的篩選。

圖6 纜索張力峰值的可行樣本分布圖

根據(jù)自動(dòng)篩選得到的最優(yōu)樣本見表4。從結(jié)果可以看出，通過優(yōu)化可以減輕纜索18%的重量。

表4 初始方案與最優(yōu)方案的變量對(duì)比

3 五級(jí)海況下最優(yōu)樣本的核算

以表4中的纜索重量和纜索長(zhǎng)度反饋給AQWA軟件，進(jìn)行水動(dòng)力時(shí)域分析。核算了在五級(jí)海況海況兩種海洋環(huán)境下的結(jié)果。圖7和圖8分別為X、Y、Z、Rx、Ry、Rz方向浮臺(tái)自由度的時(shí)間歷程曲線。

相比初始設(shè)計(jì)的纜索總重量8250kg，此時(shí)的纜索總重量為6768.65kg，減重18.0%。相比初始設(shè)計(jì)的纜索張力峰值122.12KN，此時(shí)的纜索峰值為104.80KN，降低14.2%。如表5所示。

表5 五級(jí)海況下初始方案與最優(yōu)樣本的結(jié)果對(duì)比

圖7 浮臺(tái)X、Y、Z自由度時(shí)間歷程曲線

圖8 浮臺(tái)Rx、Ry、Rz自由度時(shí)間歷程曲線

圖9 纜索的張力時(shí)間歷程曲線

4 結(jié)語

采用AQWA水動(dòng)力分析結(jié)合ANSYS Design Exploration優(yōu)化分析模塊對(duì)海洋通信浮臺(tái)的系泊系統(tǒng)進(jìn)行優(yōu)化分析。得到如下結(jié)論：

1）通過目標(biāo)驅(qū)動(dòng)的優(yōu)化設(shè)計(jì)，得到了最優(yōu)的纜索重量和纜索長(zhǎng)度，纜索共減重18.0%。

2）經(jīng)過對(duì)最優(yōu)設(shè)計(jì)參數(shù)分別在五級(jí)海況和百年一遇海況下的核算，得到浮臺(tái)的垂蕩運(yùn)動(dòng)和搖動(dòng)運(yùn)動(dòng)幅值在優(yōu)化前后變化較小，在五級(jí)海況下纜索張力峰值降低了14.2%，在百年一遇海況下纜索張力峰值降低了5.3%。

仿真計(jì)算結(jié)果證明提出的優(yōu)化方案和設(shè)計(jì)，能夠?yàn)槔|索在滿足的強(qiáng)度要求下充分考慮輕量化設(shè)計(jì)，使系泊方案具有良好的經(jīng)濟(jì)性，為海洋浮標(biāo)平臺(tái)的系泊制造設(shè)計(jì)提供了一定的解決方案。