謝仲安, 姜 哲

(上海交通大學,上海200240)

基于API規(guī)范的Spar硬艙殼體結(jié)構(gòu)優(yōu)化設(shè)計

謝仲安, 姜 哲

(上海交通大學,上海200240)

根據(jù)API規(guī)范編寫設(shè)計程序,提出Spar硬艙殼體結(jié)構(gòu)的優(yōu)化模型,采用第二代非支配排序遺傳算法NSGA-II對Spar硬艙殼體結(jié)構(gòu)進行優(yōu)化設(shè)計。通過設(shè)計實例證明,該優(yōu)化設(shè)計方法的優(yōu)化結(jié)果以及收斂性都比較理想,與原始方案相比,優(yōu)化設(shè)計方案顯著降低結(jié)構(gòu)自重。文章最后提出的建議,可作為今后實際設(shè)計工作的參考。

Spar硬艙殼體結(jié)構(gòu);API規(guī)范;NSGA-II;優(yōu)化設(shè)計

0 前言

Spar平臺主要由頂部模塊、主體結(jié)構(gòu)、系泊系統(tǒng)和立管系統(tǒng)等組成。其中,主體結(jié)構(gòu)可分為硬艙和軟艙(Truss Spar還包括桁架結(jié)構(gòu)和垂蕩板)。占主體結(jié)構(gòu)絕大部分的硬艙,又可分為甲板結(jié)構(gòu)和硬艙殼體(如圖1所示),甲板結(jié)構(gòu)徑向支撐硬艙殼體,后者直接承受軸向和徑向的外界載荷。

硬艙殼體的結(jié)構(gòu)尺寸、布置型式?jīng)Q定著Spar平臺的整體結(jié)構(gòu)強度和穩(wěn)定性。各層甲板間的殼體結(jié)構(gòu)型式基本類似,可進行模塊設(shè)計。目前,國內(nèi)對Spar平臺的研究主要集中于平臺動力響應、垂蕩板和側(cè)板的設(shè)計、系泊系統(tǒng)、關(guān)鍵節(jié)點疲勞分析等方面[1],而對平臺結(jié)構(gòu)優(yōu)化設(shè)計探討較少。本文將以兩層甲板間的殼體結(jié)構(gòu)為研究對象,探討Spar平臺硬艙殼體結(jié)構(gòu)的優(yōu)化設(shè)計方法。依據(jù)API設(shè)計規(guī)范(具體指APIBULL 2U)建立結(jié)構(gòu)設(shè)計的數(shù)學模型,采用第二代非支配排序遺傳算法NSGA-II進行優(yōu)化設(shè)計,以求減輕結(jié)構(gòu)自重、提高平臺有效載荷。

圖1 Spar硬艙殼體結(jié)構(gòu)示意圖(沿直徑剖切)

1 APIBULL 2U結(jié)構(gòu)設(shè)計方法

《API圓柱殼體穩(wěn)定性設(shè)計通告》(APIBULL 2U,以下簡稱“規(guī)范”)[2],以滿足圓柱殼體結(jié)構(gòu)穩(wěn)定性要求進行結(jié)構(gòu)設(shè)計。殼體穩(wěn)定性的計算基于線性分叉理論,通過能力折減系數(shù)和塑性折減系數(shù)將幾何非線性、材料非線性及初始缺陷等因素所起的作用考慮在內(nèi)。

規(guī)范將結(jié)構(gòu)的失效模式分為以下五類(圖2所示為各種失效模式,柱狀屈曲除外):

圖2 各失效模式示意圖(不包含柱狀屈曲)

(1)局部屈曲——加強筋之間的外殼板發(fā)生屈曲,縱向筋和環(huán)向筋仍保持原有形狀;

(2)板架屈曲——兩層環(huán)向筋之間的縱向筋連著外殼板一同發(fā)生屈曲,但環(huán)向筋仍保持原有形狀;

(3)總體屈曲——環(huán)向筋、縱向筋及所附著的外殼板一同發(fā)生屈曲;

(4)加強筋屈曲——環(huán)向或縱向加強筋單獨發(fā)生屈曲;

(5)柱狀屈曲——整個圓柱殼體發(fā)生柱狀屈曲。

規(guī)范針對各種失穩(wěn)模式計算出對應的軸向和周向臨界應力,依據(jù)畸變能理論計算出在軸向和周向載荷組合作用下,計入兩者相互作用后的結(jié)構(gòu)軸向和周向臨界應力,分別記為Fφcj和Fθcj(其中:φ表示軸向;θ表示周向;c表示考慮非線性因素折減后的結(jié)果;j表示各種失效模式),計及安全系數(shù)后,分別取軸向和周向各自的最小值為結(jié)構(gòu)軸向和周向許用應力,記為Fa和Fθ;同時根據(jù)外載計算出結(jié)構(gòu)各部位軸向和周向上的實際應力,分別記為fa和fθ。設(shè)計方案必須滿足:

若不滿足(1)式,則需調(diào)整結(jié)構(gòu)尺寸直至達到要求為止。

圖3展示了APIBULL 2U的設(shè)計流程。Spar硬艙殼體結(jié)構(gòu)通常采用環(huán)、縱向聯(lián)合加強布置型式,故實際設(shè)計流程如圖3中實線所示。

圖3 APIBULL 2U設(shè)計流程(實線表示本文的設(shè)計過程)

2 Spar硬艙殼體結(jié)構(gòu)優(yōu)化模型

APIBULL 2U的設(shè)計過程比較復雜,存在眾多約束條件和大量離散型函數(shù)關(guān)系。對于這種情況,適合采用遺傳算法(Genetic A lgo rithm s,簡稱GA)進行優(yōu)化計算。GA是模擬生物在自然環(huán)境中的遺傳和進化過程而形成的一種自適應全局優(yōu)化算法[3],通過編碼、選擇、交叉、變異等操作實現(xiàn)全局概率搜索。非支配排序遺傳算法(NSGA)是GA算法的一次升級,是一種基于Pareto最優(yōu)概念的遺傳算法,由Srinivas和Deb于1995年提出[4]。它采用非支配分層方法,使優(yōu)的個體有更大的機會遺傳到下一代,并采用適應度共享策略保證非劣最優(yōu)解均勻分布,以保持群體多樣性,避免優(yōu)化結(jié)果過早地收斂到非最優(yōu)解。本文采用的優(yōu)化算法為NSGA-II,由Deb等人于2002年提出[5],它在計算的快速性和全局性方面有了進一步的改進。

2.1 設(shè)計變量

Spar平臺硬艙段殼體結(jié)構(gòu)如圖1所示。設(shè)計變量以向量形式表示為

式中:t為外殼板厚;RW為環(huán)向加強筋翼板寬度;RH為環(huán)向加強筋腹板高度;Rtf為環(huán)向加強筋翼板厚度;Rtw為環(huán)向加強筋腹板厚度;Nr為環(huán)向加強筋數(shù)量;SW為縱向加強筋翼板寬度;SH為縱向加強筋腹板高度;Stf為縱向加強筋翼板厚度;Stw為縱向加強筋腹板厚度;Ns為縱向加強筋數(shù)量。

2.2 目標函數(shù)

設(shè)計目標是:在滿足結(jié)構(gòu)穩(wěn)定性要求的前提下,使硬艙段殼體結(jié)構(gòu)的自重盡可能小。故目標函數(shù)為

式中:X為設(shè)計變量向量;G(X)表示結(jié)構(gòu)自重,為設(shè)計變量向量X的函數(shù)。由于所有構(gòu)件都采用相同材料,具有相同的密度,故

式中:ρ為材料密度;V(X)表示結(jié)構(gòu)體積,為設(shè)計變量向量X的函數(shù)。因為ρ為定值,所以可化簡得到最終的目標函數(shù):

式中:Lb為艙段長度;R為外殼半徑。

2.3 約束條件

根據(jù)APIBULL 2U的規(guī)定,設(shè)計方案必須滿足以下約束條件。

(1)軸向和周向的許用應力需分別大于軸向和周向的實際應力:

式中:FSj表示各失效模式下的安全系數(shù);其余符號意義與上節(jié)所述相同。

(2)為保證不發(fā)生加強筋屈曲,其尺寸必須滿足:

式中:E為彈性模量;Fy為屈服應力。

3 優(yōu)化實例及結(jié)果分析

本文運用Microsoft Office Excel與MA TLAB編寫結(jié)構(gòu)設(shè)計程序,并將該程序嵌入優(yōu)化軟件Isight,采用NSGA-II算法實現(xiàn)Spar硬艙殼體結(jié)構(gòu)的優(yōu)化設(shè)計。遺傳算法的運行參數(shù)包括群體大小(Population Size)、遺傳代數(shù)(Num ber of Generations)、交叉概率(Crossover Probablity)、交叉分布指數(shù)(Crossover Distribution Index)和變異分布指數(shù)(Mutation Distribution Index)。這些運行參數(shù)對遺傳算法的求解結(jié)果和求解效率都有一定影響,本文采用Isight默認值,如表1所示。

Spar硬艙段殼體結(jié)構(gòu)的主尺度參照CTSpar概念設(shè)計方案[6]所給出的數(shù)據(jù),假設(shè)兩層甲板間的距離為10 m,置于水面以下10 m～20 m處。結(jié)構(gòu)承受平臺頂部模塊的壓力和四周海水壓力,如圖4所示。根據(jù)概念設(shè)計方案[6],分攤到單個殼體上的軸向載荷P約為700 t(即6.68×106N),均勻分布在軸向截面上;上端徑向載荷pu取水面以下10 m處壓強(即0.10045 MPa),下端徑向載荷pl取水面以下20 m處壓強(即0.2009 MPa),兩端之間線性插值。

表1 遺傳算法運行參數(shù)

圖4 殼體結(jié)構(gòu)受力示意圖

為了檢驗本文優(yōu)化設(shè)計方法的收斂性,對設(shè)計變量選取三個互不相同的初始方案分別進行優(yōu)化,初始值與優(yōu)化結(jié)果匯總于表2。

表2 優(yōu)化結(jié)果匯總 單位:mm

圖5 各方案優(yōu)化前后體積比較

分析優(yōu)化計算結(jié)果,可以得到以下結(jié)論:

(1)NSGA-II算法的優(yōu)化效果顯著。對于三個方案分別進行優(yōu)化設(shè)計后,結(jié)構(gòu)體積比優(yōu)化前下降了21.23%、36.19%和36.48%,如圖5所示。從算例的優(yōu)化結(jié)果來看,運用NSGA-II進行結(jié)構(gòu)優(yōu)化設(shè)計是有效的。

(2)各方案優(yōu)化結(jié)果基本趨同。三個方案優(yōu)化后的目標函數(shù)值(體積)最大相差僅1.59%,基本達到一致。在三個設(shè)計方案中,方案一優(yōu)化后的目標函數(shù)值最小,結(jié)構(gòu)自重最輕,所以該方案的優(yōu)化結(jié)果為最佳設(shè)計方案。

(3)各方案的優(yōu)化結(jié)果呈現(xiàn)整體趨同,局部多樣的特點。各方案經(jīng)優(yōu)化后的設(shè)計變量雖存在差異,但影響結(jié)構(gòu)總體布局的設(shè)計變量(如外殼板厚t、環(huán)向筋數(shù)量Nr)優(yōu)化結(jié)果基本一致。這說明了外載條件一定時,最優(yōu)的整體結(jié)構(gòu)布局方案是唯一的。但細化到局部構(gòu)件尺寸,可能存在多種設(shè)計變量組合同時達到最優(yōu)或次優(yōu)。產(chǎn)生這種情況的原因主要有兩點:①遺傳算法是一種基于概率的隨機搜索方法,得到的結(jié)果存在一定的偶然性,即使從同一初始值開始優(yōu)化計算,每次計算結(jié)果也未必相同;②優(yōu)化結(jié)果的偶然性隨著遺傳代數(shù)的增加而逐步降低,若增加本算例的遺傳代數(shù),可得到更為收斂的結(jié)果,但同時也要為此付出更大的計算成本和時間成本。

優(yōu)化的結(jié)果是否可行需經(jīng)過強度校核。由于優(yōu)化過程基于APIBULL 2U結(jié)構(gòu)設(shè)計規(guī)范,而此規(guī)范從穩(wěn)定性的角度出發(fā)設(shè)計結(jié)構(gòu),所以優(yōu)化結(jié)果必然滿足穩(wěn)定性要求,故只需對結(jié)構(gòu)強度進行校核。運用Patran/Nastran建立有限元模型并分析計算,得到三種優(yōu)化方案的Von Mises應力云圖,如圖6所示。三個方案的最大Von Mises應力分別為40.2 MPa、45.1 MPa和40.2 MPa,皆遠低于材料屈服應力550 MPa,強度滿足要求。

圖6 三種優(yōu)化方案結(jié)構(gòu)Von Mises應力云圖

4 結(jié)論

本文以APIBULL 2U為設(shè)計指導,將硬艙殼體結(jié)構(gòu)設(shè)計程序與優(yōu)化軟件Isight相結(jié)合,運用NSGA-II遺傳算法實現(xiàn)硬艙殼體結(jié)構(gòu)的優(yōu)化設(shè)計。通過實例計算和強度校核,驗證了本文所提出的優(yōu)化設(shè)計方法的可行性,為日后實際設(shè)計提供參考。

最后,本文作以下幾點總結(jié):

(1)采用NSGA-II算法對Spar硬艙殼體結(jié)構(gòu)進行優(yōu)化設(shè)計是可行的,能夠在保證平臺結(jié)構(gòu)性能的前提下,達到降低結(jié)構(gòu)自重,提高有效載重能力的目的;

(2)要得到更好的優(yōu)化結(jié)果,需增加遺傳代數(shù),提高計算時間成本。所以在設(shè)計時需要注意權(quán)衡優(yōu)化效果和計算成本之間的矛盾,選取合適的遺傳代數(shù);

(3)實際設(shè)計工作中需要考慮更多因素的限制,如加工工藝、生產(chǎn)條件及型材通用性等。但計入這些因素后優(yōu)化方案不一定理想。故建議在設(shè)計伊始選取多個初始設(shè)計變量向量同時優(yōu)化,得到多個優(yōu)化方案,從中遴選出便于實際生產(chǎn)的方案,以此解決理論設(shè)計與實際生產(chǎn)之間的矛盾。

[1] 張帆,楊建民,李潤培.Spar平臺的發(fā)展趨勢及其關(guān)鍵技術(shù)[J].中國海洋平臺,2005,20(2):6-11.

[2] American Petroleum Institute.Bulletin on Stability Design of Cylindrical Shells[S],APIBULLETIN 2U,Second Edition,2000.

[3] 周明,孫樹棟.遺傳算法原理及應用[M].北京:國防工業(yè)出版社,1999.

[4] Srinivas N,Deb K.Multiobjective op timization using nondominated sorting genetic algo rithms[J].Evolutionary Computation,1995,2(3):221-248.

[5] Deb K,Pratap A,Agarwal S,Meyarivan T.A fast elitistmulti-objective genetic algorithm:NSGA-II(2000)[J].IEEE Transactions on Evolutionary Computation,2002,6(2):182-197.

[6] 石紅珊.深水單立柱采油平臺方案設(shè)計及穩(wěn)性研究[D].上海:上海交通大學船舶海洋與建筑工程學院,2007.

Optim ization Design of Hard Tank Shell of Spar Based on API Rules

XIE Zhong-an, JIANG Zhe

(Shanghai Jiao Tong University,Shanghai 200240,China)

The op timization method of hard tank shell of Spar is investigated in this paper.Based on API rules,the op timization formulation is proposed and design codes are written.Non-dominated Sorting Genetic A lgorithm II(NSGA-II)is adop ted in the op timization.The op timization technique converges well and the results are good.After the op timization,the structuralweight is reduced remarkably.Helpful advice isalso provided concerning practical design of similar structures.

Spar hard tank shell;API rules;NSGA-II;op timization

P752

A

1001-4500(2010)03-0036-05

2009-09-28; 修改稿收到日期:2010-02-26

上海市科學技術(shù)委員會重大基礎(chǔ)課題資助項目(05DJ14001)

謝仲安(1985-),男,碩士研究生,主要從事海洋工程結(jié)構(gòu)設(shè)計研究。