生長(zhǎng)進(jìn)化拓?fù)鋬?yōu)化算法在油船中剖面結(jié)構(gòu)優(yōu)化上的運(yùn)用

劉宏亮，夏利娟，吳嘉蒙

（1.上海交通大學(xué)海洋工程國(guó)家重點(diǎn)實(shí)驗(yàn)室，上海200240；2.中國(guó)船舶及海洋工程設(shè)計(jì)研究院，上海200011）

生長(zhǎng)進(jìn)化拓?fù)鋬?yōu)化算法在油船中剖面結(jié)構(gòu)優(yōu)化上的運(yùn)用

劉宏亮1，夏利娟1，吳嘉蒙2

（1.上海交通大學(xué)海洋工程國(guó)家重點(diǎn)實(shí)驗(yàn)室，上海200240；2.中國(guó)船舶及海洋工程設(shè)計(jì)研究院，上海200011）

單元進(jìn)化生長(zhǎng)拓?fù)鋬?yōu)化算法較之其他拓?fù)鋬?yōu)化算法有其獨(dú)特之處，文章通過(guò)引入單元生長(zhǎng)進(jìn)化算法對(duì)VLCC中剖面橫撐結(jié)構(gòu)進(jìn)行優(yōu)化設(shè)計(jì)。結(jié)合單元應(yīng)力指標(biāo)函數(shù)及單元權(quán)重系數(shù)確定優(yōu)化模型的目標(biāo)函數(shù)，通過(guò)計(jì)算單元應(yīng)力確定單元權(quán)重系數(shù)的迭代遞推關(guān)系式，結(jié)合有限元平衡方程建立優(yōu)化數(shù)學(xué)模型，通過(guò)解析計(jì)算知優(yōu)化模型目標(biāo)函數(shù)的取值主要取決于體積函數(shù)值，依此建立完整的單元生長(zhǎng)進(jìn)化優(yōu)化流程。通過(guò)與ANSYS連續(xù)體拓?fù)鋬?yōu)化計(jì)算的比較可知單元生長(zhǎng)進(jìn)化算法更加簡(jiǎn)單高效，對(duì)VLCC輕量化設(shè)計(jì)具有一定的理論參考價(jià)值。

單元生長(zhǎng)進(jìn)化算法；拓?fù)鋬?yōu)化；應(yīng)力指標(biāo)函數(shù)；VLCC中剖面優(yōu)化

0 引言

圖1 VLCC典型的橫撐布置Fig.1 Typical VLCC’s transverse brace arrangement

目前，我國(guó)設(shè)計(jì)和建造的油船與國(guó)外同類(lèi)型船相比還是存在一定差距，主要反映在空船重量偏高、油耗偏高、建造效率偏低等，其中空船重量的差異直接影響了諸如油耗、建造成本、船舶能效設(shè)計(jì)指數(shù)（EEDI）等指標(biāo)。具體地說(shuō)，目前國(guó)內(nèi)建造的32萬(wàn)噸VLCC的平均空船重量約為46 500噸，而韓國(guó)建造的同尺度VLCC的平均空船重量?jī)H為44 500噸左右。因此，有必要運(yùn)用結(jié)構(gòu)優(yōu)化設(shè)計(jì)方法，改進(jìn)現(xiàn)有的油船結(jié)構(gòu)設(shè)計(jì)，以取得更好的經(jīng)濟(jì)效益。

橫撐結(jié)構(gòu)是VLCC特有的結(jié)構(gòu)構(gòu)件，如圖1所示，其主要作用是支撐油船縱艙壁上的垂直桁結(jié)構(gòu)，可布置在中貨艙和邊貨油艙。目前，30萬(wàn)噸級(jí)VLCC結(jié)構(gòu)中1根橫撐的結(jié)構(gòu)重量約為14噸左右，而整個(gè)貨艙區(qū)的橫撐數(shù)量為25～35根，如能取消橫撐，預(yù)估將減輕約200噸左右的結(jié)構(gòu)重量，因此，有必要對(duì)橫撐的中剖面結(jié)構(gòu)進(jìn)行優(yōu)化研究。

對(duì)于單元生長(zhǎng)進(jìn)化拓?fù)鋬?yōu)化國(guó)內(nèi)外很多學(xué)者取得了很多新的研究進(jìn)展。劉遠(yuǎn)東等[1]運(yùn)用單元生長(zhǎng)進(jìn)化方法對(duì)圓形阻尼層合板的阻尼特性進(jìn)行優(yōu)化設(shè)計(jì)，以阻尼材料用量為約束條件，以提高阻尼層的利用率為目標(biāo)，找到了圓形阻尼層合板的拓?fù)鋬?yōu)化構(gòu)型，并用相關(guān)實(shí)驗(yàn)驗(yàn)證了優(yōu)化設(shè)計(jì)的準(zhǔn)確性；王仁華、Chan和Kaveh等[2-6]將智能算法形狀退火算法運(yùn)用到結(jié)構(gòu)拓?fù)鋬?yōu)化設(shè)計(jì)中，通過(guò)與其他拓?fù)鋬?yōu)化方法的比較驗(yàn)證了該方法的可行性；何林偉、Amstutz等[7-8]基于基結(jié)構(gòu)，通過(guò)比較前后兩次結(jié)構(gòu)單元的性能指標(biāo)來(lái)增加和刪除單元，提供了一種類(lèi)似單元生長(zhǎng)進(jìn)化方法的連續(xù)體結(jié)構(gòu)拓?fù)鋬?yōu)化的簡(jiǎn)單實(shí)現(xiàn)；張宇、Lotfi等[9-10]基于SIMP插值模型，給出了一種基于A(yíng)NSYS的濾波控制算法，通過(guò)具體的數(shù)值算例驗(yàn)證了方法的可行性；胡興國(guó)、Rong等[11-12]改進(jìn)了傳統(tǒng)的結(jié)構(gòu)進(jìn)化優(yōu)化算法，提出每次刪除的單元數(shù)量根據(jù)保留單元數(shù)成固定比例，提高了進(jìn)化優(yōu)化的效率；王雷等[13]運(yùn)用仿生優(yōu)化方法，增加承載力大區(qū)域的材料，減小承載小區(qū)域的材料，提出了一種連續(xù)體結(jié)構(gòu)拓?fù)鋬?yōu)化的新方法；左中杰、Huang、Xie等[14-15]運(yùn)用ANSYS的單元生死技術(shù)給出了一種基于漸進(jìn)結(jié)構(gòu)優(yōu)化方法的拱壩結(jié)構(gòu)曲線(xiàn)函數(shù)優(yōu)化。

上述研究提出了很多新思路,進(jìn)一步完善了結(jié)構(gòu)拓?fù)鋬?yōu)化理論體系。因此，本文根據(jù)漸進(jìn)式單元生長(zhǎng)算法，結(jié)合有限元以及變密度SIMP拓?fù)鋬?yōu)化理論，建立較完整的VLCC中剖面橫撐結(jié)構(gòu)拓?fù)鋬?yōu)化模型。將優(yōu)化結(jié)果與ANSYS連續(xù)體優(yōu)化準(zhǔn)則拓?fù)鋬?yōu)化結(jié)果進(jìn)行對(duì)比，驗(yàn)證了算法的有效性，并且提出了一種減輕VLCC空船重量的可行設(shè)計(jì)思路，對(duì)工程實(shí)際具有一定的理論參考價(jià)值。

1 考慮應(yīng)力指標(biāo)函數(shù)的單元生長(zhǎng)拓?fù)鋬?yōu)化數(shù)學(xué)模型

拓?fù)鋬?yōu)化算法很難提供解析的數(shù)學(xué)優(yōu)化模型，常用的變密度理論如SIMP法是將離散的單元變量通過(guò)相應(yīng)的假設(shè)轉(zhuǎn)變?yōu)閱卧牧蠈傩裕ㄈ鐔卧獎(jiǎng)偠?、材料密度）的連續(xù)變量問(wèn)題，從而能夠建立初步的數(shù)學(xué)模型，并進(jìn)行近似的解析數(shù)學(xué)優(yōu)化計(jì)算。單元生長(zhǎng)進(jìn)化算法是通過(guò)尋找結(jié)構(gòu)最優(yōu)傳力路徑，通過(guò)在迭代計(jì)算中增加或刪除單元從而獲得最優(yōu)的結(jié)構(gòu)拓?fù)錁?gòu)型。它是建立在單元有限元計(jì)算基礎(chǔ)上的結(jié)構(gòu)智能生長(zhǎng)優(yōu)化模式，因此可以從有限元分析如剛度矩陣的求解和等效平衡方程的建立出發(fā)，通過(guò)相應(yīng)的數(shù)學(xué)解析計(jì)算建立具體的拓?fù)鋬?yōu)化數(shù)學(xué)模型。

1.1 優(yōu)化模型自變量參數(shù)的確定

單元生長(zhǎng)進(jìn)化算法在迭代的初始要確立簡(jiǎn)單的結(jié)構(gòu)初始構(gòu)型以及相應(yīng)載荷邊界條件，在優(yōu)化迭代計(jì)算過(guò)程中，尋找結(jié)構(gòu)的最優(yōu)傳力路徑，通過(guò)在應(yīng)力較大的單元周?chē)a(chǎn)生新單元，刪除應(yīng)力較小的單元來(lái)獲得拓?fù)錁?gòu)型。因此，單元生長(zhǎng)進(jìn)化算法需建立前后的迭代關(guān)系：

判斷單元的生長(zhǎng)和刪除要考慮單元的應(yīng)力水平，但是不能直接以單元的應(yīng)力值作為單元生長(zhǎng)和刪除的標(biāo)準(zhǔn)，因?yàn)閱卧獞?yīng)力值是迭代計(jì)算后得到的結(jié)果，并且結(jié)構(gòu)的應(yīng)力水平會(huì)隨著結(jié)構(gòu)的拓?fù)渖L(zhǎng)發(fā)生很大變化，因此，建立單元的權(quán)重系數(shù)xi（i=0,1,2L n）并以此為標(biāo)準(zhǔn)判斷單元的增減，單元權(quán)重系數(shù)值在每一次迭代計(jì)算中都會(huì)隨單元的應(yīng)力水平而變化，設(shè)定兩個(gè)臨界值判斷迭代計(jì)算中單元的增減：

xi≤0.01；表示單元權(quán)重太小，單元承載較小，應(yīng)該刪除該單元。

0.01 ≤xi≤1；表示單元有一定承載，繼續(xù)保留該單元。

xi＞1；表示單元承載較大，單元權(quán)重大，應(yīng)該在該單元周?chē)a(chǎn)生新的單元，單元的生長(zhǎng)示意圖如圖2[16]所示。

另一方面，結(jié)合SIMP變密度函數(shù)方法，單元的材料屬性跟單元權(quán)重系數(shù)滿(mǎn)足如下迭代關(guān)系式：

圖2 單元生長(zhǎng)示意圖（根據(jù)文獻(xiàn)[16]）Fig.2 Element growth diagram based on the reference[16]

式中：ρ0，E0是初始的材料密度和彈性模量，xi為單元在第i次迭代計(jì)算中的單元權(quán)重系數(shù)值。

因此，單元的權(quán)重系數(shù)大小決定單元的材料屬性，并且有一個(gè)確定的臨界值確定單元的增減，以此為優(yōu)化模型的自變量函數(shù)能簡(jiǎn)化模型的迭代計(jì)算。

1.2 優(yōu)化模型目標(biāo)函數(shù)的確定

油船中剖面橫撐結(jié)構(gòu)優(yōu)化設(shè)計(jì)模型中，主要考慮在結(jié)構(gòu)應(yīng)力的許可范圍內(nèi)減輕結(jié)構(gòu)的重量；主要考慮兩個(gè)方面：拓?fù)浣Y(jié)構(gòu)的Von Mises應(yīng)力、拓?fù)浣Y(jié)構(gòu)的重量m。對(duì)于結(jié)構(gòu)的相當(dāng)應(yīng)力，只要結(jié)構(gòu)的應(yīng)力水平在一個(gè)合理范圍內(nèi)就可接受，并不要求應(yīng)力最小。另外，還需要考慮結(jié)構(gòu)拓?fù)錁?gòu)型總體的應(yīng)力水平，所以應(yīng)該根據(jù)拓?fù)錁?gòu)型的最大單元應(yīng)力σmax以及拓?fù)錁?gòu)型的平均單元應(yīng)力σmean來(lái)確定一個(gè)應(yīng)力水平指標(biāo)函數(shù)，結(jié)構(gòu)拓?fù)錁?gòu)型受力均勻程度可以用單元最大相當(dāng)應(yīng)力和平均單元相當(dāng)應(yīng)力的差值來(lái)衡量，所以，建立以下的應(yīng)力水平指標(biāo)函數(shù)：

式中：a、b是權(quán)重系數(shù)，且a+b=1；分別對(duì)應(yīng)單元應(yīng)力均勻程度以及最大的單元應(yīng)力。在油船中剖面橫撐結(jié)構(gòu)優(yōu)化設(shè)計(jì)模型中，a、b都取為0.5。

上述應(yīng)力水平指標(biāo)函數(shù)同時(shí)考慮了結(jié)構(gòu)的受力均勻程度和結(jié)構(gòu)最大應(yīng)力，較好的拓?fù)錁?gòu)型應(yīng)該有較優(yōu)的結(jié)構(gòu)傳力路徑且結(jié)構(gòu)的最大應(yīng)力在合理的范圍內(nèi)，要求有較小的最大單元應(yīng)力及較小的應(yīng)力差值，所以一個(gè)較好的拓?fù)錁?gòu)型應(yīng)該有較小的結(jié)構(gòu)應(yīng)力水平指標(biāo)函數(shù)。

考慮優(yōu)化目標(biāo)的另一個(gè)方面結(jié)構(gòu)重量m，原則上油船中剖面橫撐結(jié)構(gòu)優(yōu)化設(shè)計(jì)最直接的目標(biāo)是使得縱艙壁支撐結(jié)構(gòu)重量最小，但是結(jié)構(gòu)的應(yīng)力水平指標(biāo)函數(shù)不能太大，設(shè)定結(jié)構(gòu)的理想應(yīng)力指標(biāo)函數(shù)值為因此可以確定結(jié)構(gòu)拓?fù)鋬?yōu)化的目標(biāo)函數(shù)：所以最終優(yōu)化模型的目標(biāo)函數(shù)為：

結(jié)合有限元計(jì)算平衡方程，同時(shí)設(shè)定應(yīng)力指標(biāo)函數(shù)所允許的最大值，油船中剖面橫撐結(jié)構(gòu)拓?fù)鋬?yōu)化的數(shù)學(xué)模型如下所示：

2 優(yōu)化數(shù)學(xué)模型的解析分析

油船中剖面橫撐結(jié)構(gòu)的拓?fù)鋬?yōu)化模型中用到平面四邊形單元計(jì)算，根據(jù)平面四邊形單元的有限元計(jì)算理論，平面四邊形單元的剛度矩陣：

根據(jù)（6）式可知，油船中剖面橫撐結(jié)構(gòu)的單元生長(zhǎng)進(jìn)化優(yōu)化模型的目標(biāo)函數(shù)僅僅隨結(jié)構(gòu)的密度函數(shù)和體積函數(shù)而改變?？紤]到模型中的密度函數(shù)基本保持不變，所以?xún)?yōu)化目標(biāo)即為保持單元應(yīng)力水平指標(biāo)函數(shù)值在允許的范圍內(nèi)讓結(jié)構(gòu)的體積最小，也就是結(jié)構(gòu)重量最輕，也就是說(shuō)在單元的生長(zhǎng)和進(jìn)化過(guò)程中盡量刪除掉承載較小的單元。

3 具體的單元生長(zhǎng)進(jìn)化流程

單元生長(zhǎng)進(jìn)化算法的核心過(guò)程是單元隨著迭代進(jìn)程增加或刪除，單元生長(zhǎng)基于單元權(quán)重系數(shù)值的變化，所以（1）式所示的單元生長(zhǎng)遞推關(guān)系就等同于單元權(quán)重系數(shù)的遞推關(guān)系式。根據(jù)前述對(duì)優(yōu)化數(shù)學(xué)模型的解析計(jì)算分析可知，優(yōu)化目標(biāo)函數(shù)大小由結(jié)構(gòu)拓?fù)錁?gòu)型的體積函數(shù)V（）x決定，因此，通過(guò)尋找結(jié)構(gòu)最優(yōu)的傳力路徑，盡可能使單元受力均勻，減小承載較小的單元數(shù)目是得到最優(yōu)解的必然條件。根據(jù)拓?fù)浣Y(jié)構(gòu)單元的平均應(yīng)力σmean設(shè)定三個(gè)應(yīng)力區(qū)間來(lái)建立單元權(quán)重函數(shù)的遞推關(guān)系式，三個(gè)應(yīng)力區(qū)間如下：（0，λσmean），[λσmean，σmean），[σmean，σmax），系數(shù)λ是設(shè)計(jì)參數(shù)，在油船中剖面橫撐結(jié)構(gòu)的優(yōu)化模型中，λ＝0.3，因此單元權(quán)重的遞推關(guān)系式如下式所示：

根據(jù)拓?fù)錁?gòu)型的單元應(yīng)力計(jì)算值，當(dāng)單元應(yīng)力值在縮減區(qū)間（0，λσmean），單元的權(quán)重系數(shù)值減小，直到減小到0.01，刪除該單元；單元應(yīng)力值在平衡區(qū)間[λσmean，σmean），單元權(quán)重系數(shù)值保持不變，單元受力均勻；單元應(yīng)力值在生長(zhǎng)區(qū)間[σmean，σmax），單元權(quán)重系數(shù)將增加，直到增加到1，在單元周?chē)磮D2所示增加新單元，新單元的權(quán)重系數(shù)取為1，并且原單元的單元權(quán)重系數(shù)減小為1。因此，形成新的拓?fù)錁?gòu)型，并在此基礎(chǔ)上進(jìn)行下一次的迭代計(jì)算。

根據(jù)單元生長(zhǎng)計(jì)算的迭代關(guān)系可以確定迭代的終止條件，在單元的生長(zhǎng)進(jìn)化過(guò)程中伴隨著單元的增加或刪除，并且在每次迭代過(guò)程都會(huì)有單元的增加或刪除；另一方面，優(yōu)化模型的目標(biāo)函數(shù)要求模型單元數(shù)最小，所以確定單元生長(zhǎng)迭代的終止條件為：在N次迭代計(jì)算中，優(yōu)化模型的目標(biāo)函數(shù)僅在微小的范圍內(nèi)變動(dòng)。N的確定由具體的模型確定，油船中剖面橫撐結(jié)構(gòu)的優(yōu)化模型中N取為3。

至此，建立了完整的單元生長(zhǎng)進(jìn)化拓?fù)鋬?yōu)化算法，優(yōu)化流程圖如圖3所示。

圖3 優(yōu)化流程圖Fig.3 The optimization flow chart

4 VLCC中剖面橫撐結(jié)構(gòu)拓?fù)鋬?yōu)化

首先，確定拓?fù)鋬?yōu)化設(shè)計(jì)的邊界條件和載荷，VLCC縱艙壁是船體縱向主要構(gòu)件，垂向連接船體主甲板和船體內(nèi)底板與底縱桁，因此在單元生長(zhǎng)進(jìn)化拓?fù)鋬?yōu)化模型中定義結(jié)構(gòu)的邊界為簡(jiǎn)支約束。VLCC縱艙壁承載的極限狀態(tài)為一邊貨艙滿(mǎn)載，一邊空載，載荷示意圖如圖4所示。

圖4 拓?fù)鋬?yōu)化模型初始構(gòu)型Fig.4 The initial configuration of the optimization model

圖5 單元生長(zhǎng)進(jìn)化優(yōu)化結(jié)果示意圖Fig.5 The result from the evolutionary structural optimization

從優(yōu)化結(jié)果可知，在迭代初期，在初始構(gòu)型基礎(chǔ)上，根據(jù)響應(yīng)的邊界條件和載荷工況，拓?fù)錁?gòu)型產(chǎn)生大量新單元形成較完整的初始拓?fù)錁?gòu)型，單元數(shù)量增加明顯即拓?fù)錁?gòu)型的重量增加；進(jìn)入迭代中期階段，根據(jù)拓?fù)鋬?yōu)化目標(biāo)函數(shù)設(shè)定，對(duì)新生成的拓?fù)錁?gòu)型進(jìn)行優(yōu)化計(jì)算，在滿(mǎn)足優(yōu)化約束條件的基礎(chǔ)上，使優(yōu)化模型的目標(biāo)函數(shù)最小，即增加新單元的同時(shí)更注重尋找結(jié)構(gòu)最優(yōu)的傳力路徑，刪除應(yīng)力較小的單元，使結(jié)構(gòu)拓?fù)錁?gòu)型的單元數(shù)量盡可能減少；最后通過(guò)進(jìn)一步優(yōu)化計(jì)算得到最終的拓?fù)鋬?yōu)化構(gòu)型。

為了驗(yàn)證單元生長(zhǎng)拓?fù)鋬?yōu)化算法在VLCC中剖面橫撐結(jié)構(gòu)優(yōu)化設(shè)計(jì)中的有效性，將基于A(yíng)NSYS平臺(tái)的連續(xù)結(jié)構(gòu)拓?fù)鋬?yōu)化運(yùn)用到VLCC中剖面橫撐結(jié)構(gòu)的拓?fù)鋬?yōu)化過(guò)程中；ANSYS拓?fù)鋬?yōu)化算法采用優(yōu)化準(zhǔn)則法，優(yōu)化目標(biāo)函數(shù)為結(jié)構(gòu)構(gòu)型的柔度，以最小柔度即結(jié)構(gòu)最大剛度為目標(biāo)，優(yōu)化模型如下：

式中：Uc為結(jié)構(gòu)變性能。

優(yōu)化結(jié)果如圖6所示。

圖6 ANSYS拓?fù)鋬?yōu)化結(jié)果示意圖Fig.6 The result from the ANSYS topology optimization

對(duì)比圖5和圖6兩種方法的優(yōu)化結(jié)果示意圖可知，單元生長(zhǎng)進(jìn)化拓?fù)鋬?yōu)化算法所得的最終優(yōu)化拓?fù)錁?gòu)型更加清晰簡(jiǎn)單，通過(guò)尋找最優(yōu)的結(jié)構(gòu)傳力路徑，在應(yīng)力較大單元周?chē)尚聠卧?，刪除應(yīng)力較小的單元，得到的拓?fù)錁?gòu)型形狀清晰簡(jiǎn)單。基于A(yíng)NSYS平臺(tái)的連續(xù)結(jié)構(gòu)拓?fù)鋬?yōu)化能控制結(jié)構(gòu)的總體積，但是結(jié)構(gòu)的拓?fù)錁?gòu)型輔助結(jié)構(gòu)較多，構(gòu)型較復(fù)雜。因此總體上講，單元生長(zhǎng)進(jìn)化拓?fù)鋬?yōu)化算法在VLCC中剖面橫撐結(jié)構(gòu)的優(yōu)化上更加簡(jiǎn)單和有效，更加符合實(shí)際的優(yōu)化設(shè)計(jì)。

因此根據(jù)圖5最終的拓?fù)錁?gòu)型，VLCC縱艙壁支撐結(jié)構(gòu)的橫撐可以取消，只需加靠近船底的扶強(qiáng)材尺寸，具體的可行方案如圖7所示。

圖7 VLCC中剖面橫撐取消優(yōu)化設(shè)計(jì)圖Fig.7 The VLCC mid-section optimal design by excluding the transverse brace

5 結(jié)語(yǔ)

自1992年Xie和Steven提出了進(jìn)化結(jié)構(gòu)優(yōu)化算法（ESO）[16]以來(lái)，單元生長(zhǎng)進(jìn)化方法受到了廣大學(xué)者的關(guān)注，其拓?fù)鋬?yōu)化的簡(jiǎn)單高效性使其得到越來(lái)越廣泛的運(yùn)用。文章基于進(jìn)化結(jié)構(gòu)優(yōu)化算法，通過(guò)定義單元應(yīng)力指標(biāo)函數(shù)，結(jié)合結(jié)構(gòu)輕量化目標(biāo)設(shè)計(jì)，確定優(yōu)化模型的目標(biāo)函數(shù)；結(jié)合變密度SIMP材料理論模型，定義單元權(quán)重函數(shù)為優(yōu)化函數(shù)的自變量參數(shù)，并根據(jù)單元平均應(yīng)力給出單元權(quán)重函數(shù)的遞推迭代關(guān)系式；結(jié)合有限元平衡方程以及應(yīng)力指標(biāo)上限建立單元生長(zhǎng)進(jìn)化拓?fù)鋬?yōu)化算法的數(shù)學(xué)模型。通過(guò)對(duì)其解析計(jì)算分析，指出優(yōu)化目標(biāo)函數(shù)最小即要求拓?fù)錁?gòu)型體積函數(shù)最??；最后結(jié)合迭代終止條件給出具體的考慮應(yīng)力指標(biāo)的單元生長(zhǎng)拓?fù)鋬?yōu)化流程。

通過(guò)對(duì)VLCC中剖面橫撐結(jié)構(gòu)的優(yōu)化計(jì)算，對(duì)比了單元生長(zhǎng)進(jìn)化算法和基于A(yíng)NSYS平臺(tái)的連續(xù)體拓?fù)鋬?yōu)化準(zhǔn)則算法，從優(yōu)化結(jié)果對(duì)比中看出，單元生長(zhǎng)進(jìn)化算法能得到更加簡(jiǎn)單清晰的拓?fù)錁?gòu)型。本文最后根據(jù)優(yōu)化的拓?fù)錁?gòu)型給出了一種取消VLCC中剖面橫撐結(jié)構(gòu)的設(shè)計(jì)方案。

[1]劉遠(yuǎn)東,尹益輝,胡文軍,等.圓形阻尼層合板的拓?fù)鋬?yōu)化設(shè)計(jì)及其實(shí)驗(yàn)驗(yàn)證[J].振動(dòng)與沖擊,2012(22)：76-79. Liu Yuandong,Yin Yihui,et al.Topologic Optimal design of a circular laminated damping plate with its experimental verification[J].Journal of Vibration and Shock,2012(22)：76-79.

[2]王仁華,趙憲忠,等.平面桁架結(jié)構(gòu)拓?fù)鋬?yōu)化設(shè)計(jì)的改進(jìn)智能算法[J].工程力學(xué),2012(11)：205-211. Wang Renhua,Zhao Xianzhong,et al.Improved intelligent algorithm for topology optimization design of planar truss structure[J].Engineering Mechanics,2012(11)：205-211.

[3]Chan C M,Zhang L M,Ng J T M.Optimization of pile groups using hybrid genetic algorithms[J].Journal of Geotechnical and Geoenvironmental Engineering,2009,135(4)：497-505.

[4]Kaveh A,Talatahari S.Particle swarm optimizer,ant-colony strategy and harmony search scheme hybridized for optimization of truss structures[J].Computers&Structures,2009,87(5/6)：267-283.

[5]Kaveh A,Talatahari S.Optimum design of skeletal structures using imperialist competitive algorithm[J].Computers& Structures,2010,88(21/22)：1220-1229.

[6]Gomes H M.Truss optimization with dynamic constraints using a particle swarm algorithm[J].Expert Systems with Applications,2011,38(1)：957-968.

[7]何林偉,蔡國(guó)平.連續(xù)結(jié)構(gòu)拓?fù)鋬?yōu)化的一種簡(jiǎn)單實(shí)現(xiàn)[J].機(jī)械強(qiáng)度,2012,34(5)：687-691. Hen Linwei,Cai Guoping.Simple approach to structural topology optimization for continuous structures[J].Journal of Mechanical Strength,2012,34(5)：687-691.

[8]Amstutz S,Novotny A A．Topological optimization of structures subject to von Mises stress constraints[J]．Structural and Multidisciplinary Optimization,2010,41：407-420．

[9]張宇,羅陸峰.基于SIMP插值模型的棋盤(pán)格現(xiàn)象濾波控制算法研究[J].制造業(yè)自動(dòng)化,2012(10)：81-84. Zhang yu,Luo Lufeng.Research on filter control algorithm for the checkerboard phenomenon based on SIMP interpolation model[J].Journal of Manufacture and Automation,2012(10)：81-84.

[10]Lotfi K G,Sherif A M,Ashraf O N.A particle swarm-based genetic algorithm for schedul-ing in an agile environment[J]. Computers&Industrial Engineering,2008,55(3)：707.

[11]胡興國(guó),程赫明,位向飛.根據(jù)排序移除低效材料的漸進(jìn)結(jié)構(gòu)優(yōu)化方法研究[J].科學(xué)技術(shù)與工程,2012(10)：7261-7265. Hu Xingguo,Cheng Heming,Wei Xiangfei.Evolutionary structural optimization by removing inefficient material according to the rank of material’s efficiency[J].Science Technology and Engineering,2012(10)：7261-7265.

[12]Rong J H,Xie Y M,Yang X Y.An improved method for evolutionary structural optimization against buckling[J].Computers and Structures,2001,79(3)：253-263.

[13]王雷.連續(xù)結(jié)構(gòu)體的SKO優(yōu)化方法研究[J].設(shè)計(jì)與分析,2012(27)：135-137. Wanglei.The SKO optimization algorithm applying on the continuous structure[J].Design and Analysis,2012(27)：135-137.

[14]左中杰,干泉,王運(yùn)興,等.基于漸進(jìn)結(jié)構(gòu)優(yōu)化法的拱壩結(jié)構(gòu)曲線(xiàn)函數(shù)優(yōu)化[J].土木工程與管理學(xué)報(bào),2012(6)：93-98. Zuo Zhongjie,Gan chuan,et al.Optimization research of Arch Dam structure curve function based on the evolutionary structural optimization algorithm[J].Joural of Civil Engineering and Management,2012(6)：93-98.

[15]Huang X,Xie Y M．Bi-directional evolutionary topology optimization of continuum structures with one or multiple materials[J]．Computational Mechanics,2009,43(3)：393-401．

[16]Xie Y M,Steven G P.Shape and layout optimisation via an evolutionary procedure[C].In：Proceedings of International Conference on Computational Engineering Science.Hong Kong University of Science and Technology,1992：17-22.

Applying the evolution topology algorithm to the VLCC’s mid-section structure optimization

LIU Hong-liang1,XIA Li-jüan1,WU Jia-meng2
(1.State Key Laboratory of Ocean Engineering,Shanghai Jiaotong University,Shanghai 200240,China; 2.Marine Design and Research Institute of China,Shanghai 200011,China)

Evolutionary structural topology optimization shows some specific features compared with other topology optimization algorithms.Which can be applied on the optimization design of the VLCC mid-section cross-brace structure.According to the element’s stress index and the element weight coefficient,the objective function was established.Then,the iteration equation of the element’s weight coefficient was created based on the obtained element stress.Finally,combining the FEM’s balance equation the total optimization mathematical model was created.The mathematical calculation results show that the objective function depends mainly on the model’s volume function.Compared with the result of the ANSYS continuous structure topology optimization,the calculation with the evolution algorithm is more simple and efficient. This method using in this paper may provide something useful to the VLCC’s lightweight design.

evolutionary structural optimization;topology optimization;stress index; VLCC mid-section optimization;VLCC(very large crude oil carrier)

U662

A

10.3969/j.issn.1007-7294.2017.04.011

1007-7294（2017）04-0464-08

2016-10-21

國(guó)家自然科學(xué)基金資助項(xiàng)目（50909060）；海洋工程國(guó)家重點(diǎn)實(shí)驗(yàn)室青年創(chuàng)新基金資助（GKZD010059-20）

劉宏亮（1989-），男，碩士研究生，E-mail：76908358liang@sjtu.edu.cn；夏利娟（1975-），女，副教授。