周曉松，梅志遠(yuǎn)

海軍工程大學(xué)艦船工程系，湖北武漢430033

0 引 言

板架結(jié)構(gòu)具有承載能力強(qiáng)、重量輕的特點(diǎn)，被廣泛用于船舶結(jié)構(gòu)，且隨著復(fù)合材料在艦船結(jié)構(gòu)上的應(yīng)用，也得到了越來(lái)越廣泛的關(guān)注。對(duì)復(fù)合材料板架結(jié)構(gòu)開展優(yōu)化設(shè)計(jì)方法應(yīng)用基礎(chǔ)研究，有利于科學(xué)、合理地評(píng)價(jià)艦船復(fù)合材料結(jié)構(gòu)的優(yōu)越性，具有重要的工程價(jià)值和理論意義。

復(fù)合材料板架結(jié)構(gòu)由于具有各向異性，且界面復(fù)雜、可設(shè)計(jì)變量眾多，因而優(yōu)化設(shè)計(jì)難度大，而傳統(tǒng)的數(shù)學(xué)規(guī)劃法或準(zhǔn)則法卻很難解決此類問題。近年來(lái)，國(guó)內(nèi)外學(xué)者對(duì)此開展了廣泛的研究。吳莉莉等［1］提出了加筋板穩(wěn)定性約束下的二級(jí)協(xié)同優(yōu)化算法。喬巍等［2］采用等效彎曲剛度法實(shí)現(xiàn)了鋪層順序的優(yōu)化。趙群等［3］提出了基于壓彎剛度復(fù)合材料加筋板的布局優(yōu)化設(shè)計(jì)方法。李爍等［4］應(yīng)用神經(jīng)網(wǎng)絡(luò)對(duì)復(fù)合材料帽形筋加筋板進(jìn)行了優(yōu)化。張鐵亮等［5］采用試驗(yàn)設(shè)計(jì)建立了代理模型方法。劉克龍等［6］在對(duì)機(jī)翼的氣動(dòng)結(jié)構(gòu)進(jìn)行多學(xué)科優(yōu)化設(shè)計(jì)時(shí)運(yùn)用了低自由度協(xié)同優(yōu)化方法。Blair 等［7］將快速建模方法應(yīng)用到了翼面結(jié)構(gòu)的布局優(yōu)化設(shè)計(jì)中。Nagendra 等［8］應(yīng)用改進(jìn)的遺傳算法對(duì)加筋板結(jié)構(gòu)進(jìn)行了優(yōu)化設(shè)計(jì)。

本文將采用分級(jí)遞進(jìn)優(yōu)化設(shè)計(jì)方法，對(duì)復(fù)合材料板架結(jié)構(gòu)的設(shè)計(jì)變量進(jìn)行系統(tǒng)分解，通過對(duì)子系統(tǒng)的遞進(jìn)優(yōu)化和相互協(xié)調(diào)，解決艦船工程實(shí)踐［9］中復(fù)合材料板架結(jié)構(gòu)的多變量設(shè)計(jì)問題。

1 典型板架結(jié)構(gòu)強(qiáng)度特性的仿真分析

1.1 典型板架結(jié)構(gòu)仿真分析模型的建立

本文針對(duì)某型艦上層建筑甲板室舷側(cè)外側(cè)壁夾層板架建立初始計(jì)算模型，展開強(qiáng)度特性分析，建模采用的單位制為kg，mm 和s。板架由夾層板、4 根橫筋和4 根縱筋構(gòu)成。夾層板由上、下玻璃鋼表層和PVC 泡沫芯材組成；筋材為帽形筋結(jié)構(gòu)，包括內(nèi)部芯材與玻璃鋼表層。板架尺寸參數(shù)如表1 和圖1～圖3 所示。

表1 板架尺寸參數(shù)Tab.1 Dimension parameters of stiffened panel

圖1 結(jié)構(gòu)布置圖Fig.1 Structural layout

圖2 構(gòu)件示意圖Fig.2 Schematic diagram of structural components

圖3 局部截面示意圖Fig.3 Schematic diagram of local section

板架初始計(jì)算模型的邊界條件：考慮板架最危險(xiǎn)的承載工況，板架邊界、板架與中間甲板相接處均設(shè)置固支約束。單元類型：基于材料的力學(xué)性能和結(jié)構(gòu)的受力特征，玻璃鋼采用Shell 單元，芯材采用Solid 單元進(jìn)行建模。網(wǎng)格劃分：對(duì)板架模型進(jìn)行網(wǎng)格劃分，全局種子密度為9；筋材相接區(qū)域的網(wǎng)格進(jìn)一步細(xì)化，局部種子密度為6；板材厚度劃分為6 層；網(wǎng)格數(shù)量為1 557 591。載荷：舷側(cè)外側(cè)壁板加載均布70 kPa 載荷。材料參數(shù)和板架初始計(jì)算模型如表2、表3 和圖4 所示。

表2 材料參數(shù)（玻璃鋼）Tab.2 Material parameters（glassfiber reinforced plastics）

表3 材料參數(shù)（PVC 泡沫芯材）Tab.3 Material parameters（PVC foam core）

圖4 板架初始計(jì)算模型Fig.4 Initial computational model of stiffened panel

1.2 計(jì)算結(jié)果

通過計(jì)算，得到舷側(cè)外側(cè)壁板架結(jié)構(gòu)的位移云圖和應(yīng)力云圖如圖5～圖7 所示，其中位移和應(yīng)力單位分別為mm 和kPa。

分析計(jì)算結(jié)果可知，復(fù)合材料板架結(jié)構(gòu)在均布載荷作用下的變形規(guī)律與鋼質(zhì)板架基本相同。玻璃鋼的應(yīng)力峰值為450.2 MPa，整體應(yīng)力水平較低，約為80～220 MPa；芯材的應(yīng)力峰值為2.29 MPa，整體應(yīng)力水平約為0.8～1.5 MPa。

圖5 位移云圖Fig.5 Displacement contours

圖6 玻璃鋼應(yīng)力云圖Fig.6 Stress contours of glassfiber reinforced plastics

圖7 芯材應(yīng)力云圖Fig.7 Stress contours of PVC foam core

由于板架結(jié)構(gòu)建模對(duì)連接區(qū)域進(jìn)行了簡(jiǎn)化，因而局部應(yīng)力集中問題突出，主要集中在筋材相接區(qū)域和邊界區(qū)域。在整體優(yōu)化設(shè)計(jì)階段，將僅以中心區(qū)域板格（去除邊界板格與局部應(yīng)力集中區(qū)域）的應(yīng)力水平來(lái)評(píng)價(jià)板架的承載能力，局部連接結(jié)構(gòu)與邊界約束較強(qiáng)區(qū)域峰值過高的問題暫不予考慮，可通過局部結(jié)構(gòu)細(xì)節(jié)設(shè)計(jì)和結(jié)構(gòu)形式優(yōu)化進(jìn)行重點(diǎn)分析。板架結(jié)構(gòu)質(zhì)量為486 kg，質(zhì)量較大，有必要通過優(yōu)化來(lái)減輕質(zhì)量，提高板架結(jié)構(gòu)的承載效率。

板架結(jié)構(gòu)的優(yōu)化設(shè)計(jì)變量眾多，包括夾層板芯材厚度、縱橫加強(qiáng)筋截面尺寸、復(fù)合材料鋪層厚度、復(fù)合材料鋪層角度和鋪層順序，要想一次全部?jī)?yōu)化極為困難，需要對(duì)優(yōu)化設(shè)計(jì)變量進(jìn)行分析，采用分步優(yōu)化的策略進(jìn)行。

2 分級(jí)遞進(jìn)優(yōu)化設(shè)計(jì)方法介紹

針對(duì)復(fù)合材料板架結(jié)構(gòu)優(yōu)化設(shè)計(jì)這個(gè)系統(tǒng)問題，采用了分級(jí)遞進(jìn)優(yōu)化設(shè)計(jì)方法。該方法是將一個(gè)復(fù)雜的系統(tǒng)優(yōu)化問題分解為若干較簡(jiǎn)單的子系統(tǒng)，通過對(duì)各子系統(tǒng)分別進(jìn)行優(yōu)化和級(jí)間協(xié)調(diào)來(lái)實(shí)現(xiàn)原系統(tǒng)的優(yōu)化。各個(gè)子系統(tǒng)雖然相互關(guān)聯(lián)，但又有其相對(duì)獨(dú)立性，根據(jù)對(duì)使用功能、承載特點(diǎn)和結(jié)構(gòu)形式的分析，若在獨(dú)立性強(qiáng)、關(guān)聯(lián)性弱的部位加以分解，只需很少幾次迭代，就能收斂到滿意的優(yōu)化解。

對(duì)板架結(jié)構(gòu)全部?jī)?yōu)化設(shè)計(jì)變量進(jìn)行靈敏度分析，將優(yōu)化設(shè)計(jì)變量分為2 類：第1 類為夾層板芯材厚度、加強(qiáng)筋芯材尺寸和表層玻璃鋼厚度；第2類為給定角度下玻璃鋼的鋪層數(shù)量和鋪層順序。分析每類優(yōu)化設(shè)計(jì)變量的特點(diǎn)，進(jìn)而將板架優(yōu)化工作分為宏觀尺寸優(yōu)化和鋪層方式優(yōu)化兩級(jí)進(jìn)行。該方法思路簡(jiǎn)明，實(shí)施方便，易于與常用的工程結(jié)構(gòu)分析算法和強(qiáng)度衡準(zhǔn)協(xié)調(diào)，并且還可降低所有優(yōu)化設(shè)計(jì)變量一起優(yōu)化所帶來(lái)的耦合影響，從而達(dá)到最佳優(yōu)化效果，如圖8 所示。

圖8 二級(jí)系統(tǒng)示意圖Fig.8 Schematic diagram of two-stage system

板架結(jié)構(gòu)優(yōu)化設(shè)計(jì)的數(shù)學(xué)模型可表述如下：

最小化優(yōu)化目標(biāo)函數(shù)

對(duì)于各子系統(tǒng)，狀態(tài)變數(shù)Xi與輸入、輸出變數(shù)之間存在下述等式與不等式約束：

其中：Xi為子系統(tǒng)i 的狀態(tài)變數(shù)，即板架結(jié)構(gòu)使用材料的極限強(qiáng)度值；Ui為子系統(tǒng)i 的輸入變數(shù)，即系統(tǒng)i 優(yōu)化前的設(shè)計(jì)變量值；Yi為子系統(tǒng)i 的輸出變數(shù)，即系統(tǒng)i 優(yōu)化后的設(shè)計(jì)變量值，同時(shí)又是其他子系統(tǒng)的干涉變數(shù)。

目標(biāo)函數(shù)被分離到每一個(gè)子系統(tǒng)，但由于存在相互聯(lián)系的干涉變數(shù)，因而需要采用逐次迭代的方式來(lái)實(shí)現(xiàn)。首先，將干涉變數(shù)人為地加以固定，使各子系統(tǒng)完全被分割，即給定某一組變量Zi，并使Y=Z 。然后，在第1 級(jí)優(yōu)化的基礎(chǔ)上逐次修正干涉變數(shù)，使總目標(biāo)函數(shù)取極值。因此，二級(jí)優(yōu)化結(jié)構(gòu)成為下述形式。

第1 級(jí)：

第2 級(jí)：

在求解時(shí)，首先給定干涉變數(shù)的一個(gè)初始估計(jì)值Z，在第1 級(jí)對(duì)每一個(gè)子系統(tǒng)分別進(jìn)行優(yōu)化，然后再將第1 級(jí)所得結(jié)果傳遞到第2 級(jí)，在第2 級(jí)通過目標(biāo)函數(shù)h（Z）的最優(yōu)化來(lái)修改Z，并最后確定其最優(yōu)值。

2.1 一級(jí)優(yōu)化分析

宏觀尺寸優(yōu)化階段的數(shù)學(xué)模型可具體表述為：

最小化目標(biāo)函數(shù)fi（u）

約束gi（u）gj

U≤0，（j=1，…，X）

其中：f （iu）為目標(biāo)函數(shù)；g （iu）和gjU分別為第j 個(gè)約束響應(yīng)和響應(yīng)的上限值；X 為所有約束的數(shù)目；YE 為單元數(shù)目；Yp為一級(jí)優(yōu)化設(shè)計(jì)變量的數(shù)目；uik為第k 個(gè)單元所在的第i 個(gè)一級(jí)優(yōu)化設(shè)計(jì)變量的厚度，且設(shè)計(jì)變量的厚度是可以連續(xù)變化的。在宏觀尺寸優(yōu)化階段，只考慮全局響應(yīng)和非強(qiáng)制的工藝約束。本文中，此階段考慮的全局響應(yīng)包括靜強(qiáng)度響應(yīng)和質(zhì)量響應(yīng)，將靜強(qiáng)度響應(yīng)（即板架結(jié)構(gòu)使用材料的極限強(qiáng)度值）作為約束，質(zhì)量響應(yīng)作為目標(biāo)。非強(qiáng)制的工藝約束包括一級(jí)優(yōu)化設(shè)計(jì)變量尺寸的設(shè)計(jì)范圍，即優(yōu)化后的板架構(gòu)件尺寸值不大于初始板架構(gòu)件尺寸。

2.2 二級(jí)優(yōu)化分析

鋪層方式優(yōu)化可以進(jìn)一步分為層組尺寸優(yōu)化和鋪層次序優(yōu)化2 個(gè)步驟。

宏觀尺寸優(yōu)化的結(jié)果是復(fù)合材料鋪層的總厚度。在層組尺寸優(yōu)化設(shè)計(jì)中，同一方向和布局的多層被認(rèn)為是一個(gè)層組，這可以在很大程度上減少層數(shù)，便于通過改變每個(gè)層組的厚度來(lái)模擬鋪層的增加或減少。本文采用的鋪層為垂直正交的纖維布，采用0°和45°鋪層，從而可以將復(fù)合材料層合板以兩級(jí)層組的形式進(jìn)行建模。由于層組厚度可以獨(dú)立變化，這種建模方式有利于強(qiáng)調(diào)層合板的優(yōu)化位置及層組厚度，因而可以從最終的厚度優(yōu)化中得到需要的鋪層數(shù)。

層組尺寸優(yōu)化后，鋪層布局與鋪層細(xì)節(jié)便均完成，但詳細(xì)的工藝約束可能還不滿足。因此，本階段在滿足所有設(shè)計(jì)約束的條件下，將對(duì)所有的鋪層進(jìn)行順序優(yōu)化，以滿足所有的工藝約束。

3 典型板架結(jié)構(gòu)優(yōu)化實(shí)例分析

3.1 靈敏度分析

圖9 設(shè)計(jì)變量靈敏度分析圖Fig.9 Screenshot of design variable sensitivity

對(duì)上文建立的板架結(jié)構(gòu)模型進(jìn)行優(yōu)化設(shè)計(jì)變量的靈敏度分析，其結(jié)果如圖9 所示。通過分析結(jié)果，表明在初始板架結(jié)構(gòu)模型分析中，夾層板芯材厚度對(duì)板架重量影響的靈敏度為58%；玻璃鋼鋪層厚度對(duì)板架的應(yīng)力水平影響最大，靈敏度為32%；縱向加強(qiáng)筋尺寸對(duì)板架應(yīng)力水平的影響大于橫筋的，其靈敏度分別為29%和26%；其中縱筋高度和橫筋寬度對(duì)板架應(yīng)力水平的影響較大，靈敏度分別為23%和21%。綜合分析所有優(yōu)化設(shè)計(jì)變量后可知，夾層板芯材厚度、玻璃鋼鋪層厚度與加強(qiáng)筋芯材尺寸等宏觀設(shè)計(jì)變量對(duì)優(yōu)化目標(biāo)的影響較大，而玻璃鋼鋪層角度和鋪層順序的影響則相對(duì)較小，因而可將復(fù)合材料板架結(jié)構(gòu)優(yōu)化分為宏觀尺寸優(yōu)化和鋪層方式優(yōu)化兩級(jí)進(jìn)行。

3.2 一級(jí)優(yōu)化分析

在上層建筑結(jié)構(gòu)設(shè)計(jì)中，質(zhì)量是重要的設(shè)計(jì)指標(biāo)，上層建筑結(jié)構(gòu)質(zhì)量的增加會(huì)導(dǎo)致艦船重心升高，從而影響艦船穩(wěn)定性。因此，建立了以質(zhì)量為目標(biāo)函數(shù)，以靜強(qiáng)度為約束條件，以?shī)A層板芯材厚度、加強(qiáng)筋芯材截面尺寸和玻璃鋼鋪層厚度為一級(jí)優(yōu)化設(shè)計(jì)變量的有限元優(yōu)化模型，優(yōu)化目標(biāo)為板架結(jié)構(gòu)質(zhì)量最小。通過優(yōu)化，確定夾層板芯材厚度、加強(qiáng)筋芯材截面尺寸和玻璃鋼鋪層厚度。

3.3 二級(jí)優(yōu)化分析

本階段需要綜合所有的設(shè)計(jì)響應(yīng)和工藝約束，確定規(guī)定角度下的鋪層數(shù)量和鋪層順序。

層組尺寸優(yōu)化階段的模型與宏觀尺寸優(yōu)化階段的相同，但優(yōu)化設(shè)計(jì)變量變?yōu)榱藢咏M厚度，并且增加了非強(qiáng)制的工藝約束條件。最終，確定復(fù)合材料各角度層組的厚度和鋪層數(shù)量，進(jìn)而得到總的鋪層數(shù)量。在整個(gè)層組尺寸優(yōu)化過程中，所有的約束違反均為0%。層組尺寸優(yōu)化后，將復(fù)合材料層組離散為由基本厚度為0.31 mm 鋪層組成的不等厚的復(fù)合材料層合板。

在層組尺寸優(yōu)化的基礎(chǔ)上開展鋪層次序優(yōu)化，用以確定最終的鋪層順序。本階段優(yōu)化時(shí)的所有響應(yīng)約束與層組優(yōu)化階段相同，但設(shè)置了更為詳細(xì)的工藝約束，如外表面鋪層使用45°的鋪層，以提高抗剪能力；保證鋪層數(shù)量和位置對(duì)稱分布，減小制造過程中可能產(chǎn)生的應(yīng)力；限制同一角度的鋪層數(shù)量，以減小邊緣分層現(xiàn)象的發(fā)生。優(yōu)化后的復(fù)合材料鋪層次序?yàn)椤?5/0/45/07/45/0/45】。

4 優(yōu)化前、后結(jié)果比較分析

復(fù)合材料板架結(jié)構(gòu)經(jīng)過兩級(jí)優(yōu)化后，其位移云圖、應(yīng)力云圖和質(zhì)量比較分別如圖10～圖12 及表4、表5 所示。

圖10 優(yōu)化后的位移云圖Fig.10 Displacement contours after optimization

圖11 優(yōu)化后的璃鋼應(yīng)力云圖Fig.11 Stress contours of glass fiber reinforced plastics after optimization

圖12 優(yōu)化后的芯材應(yīng)力云圖Fig.12 Stress contours of PVC foam core after optimization

表4 復(fù)合材料優(yōu)化前、后宏觀尺寸比較Tab.4 Comparision of macroscopic size of composite materials before and after optimization

表5 板架結(jié)構(gòu)優(yōu)化前、后性能參數(shù)比較Tab.5 Comparision of performance parameters of stiffened panel before and after optimization

通過優(yōu)化，夾層板芯材厚度由40 mm 減少至32 mm，減少了20%；玻璃鋼厚度由5 mm 減至4 mm，減少了20%；板架結(jié)構(gòu)總質(zhì)量從486 kg 減輕至391 kg，減輕了19.5%，減重效果明顯。其中橫筋尺寸相對(duì)于縱筋尺寸變化較小，表明在板架結(jié)構(gòu)中橫筋所起的承載作用更大。優(yōu)化后，結(jié)構(gòu)的最大位移及最大應(yīng)力與優(yōu)化前相比均增大，邊界約束區(qū)域的應(yīng)力峰值已達(dá)到材料的強(qiáng)度極限，但仍滿足板架結(jié)構(gòu)的設(shè)計(jì)約束要求。在此處，可只關(guān)注整體板架的承載能力，局部連接結(jié)構(gòu)與邊界約束較強(qiáng)區(qū)域峰值過高的問題暫不予考慮，應(yīng)通過局部結(jié)構(gòu)細(xì)節(jié)設(shè)計(jì)和結(jié)構(gòu)形式優(yōu)化予以重點(diǎn)分析。

5 結(jié) 論

1）基于Isight 優(yōu)化平臺(tái)的分級(jí)遞進(jìn)優(yōu)化設(shè)計(jì)方法能夠滿足復(fù)合材料夾層板架結(jié)構(gòu)的優(yōu)化設(shè)計(jì)要求，優(yōu)化后，夾層板表層和芯材的應(yīng)力水平趨向于約束要求值，能滿足設(shè)計(jì)要求。

2）均布載荷作用下，對(duì)復(fù)合材料板架結(jié)構(gòu)承載能力影響最大的設(shè)計(jì)變量為夾層板表層玻璃鋼厚度和芯材厚度，其次為加強(qiáng)筋芯材尺寸。在復(fù)合材料板架結(jié)構(gòu)優(yōu)化設(shè)計(jì)中，應(yīng)重點(diǎn)分析，提高優(yōu)化效率。

3）采用0°和45°層組進(jìn)行鋪層，經(jīng)宏觀尺寸優(yōu)化和鋪層方式優(yōu)化這兩個(gè)階段后，在滿足板架結(jié)構(gòu)強(qiáng)度要求的前提下，復(fù)合材料板架結(jié)構(gòu)的質(zhì)量從486 kg 下降到了391 kg，共減少95 kg，減重達(dá)19.5%，說明優(yōu)化方法效率很高，且優(yōu)化結(jié)果滿足制造工藝要求，具有很高的工程應(yīng)用價(jià)值。

4）文中板架的優(yōu)化結(jié)果僅以中心區(qū)域板格的應(yīng)力水平來(lái)評(píng)價(jià)整體板架的承載能力，是對(duì)板架模型進(jìn)行簡(jiǎn)化后的優(yōu)化結(jié)果，而對(duì)局部連接結(jié)構(gòu)與邊界約束較強(qiáng)區(qū)域峰值過高的問題則并未考慮，這可在下一步局部結(jié)構(gòu)細(xì)節(jié)設(shè)計(jì)和結(jié)構(gòu)形式優(yōu)化時(shí)進(jìn)行重點(diǎn)分析。

［1］吳莉莉，姚衛(wèi)星.復(fù)合材料加筋板結(jié)構(gòu)的二級(jí)協(xié)同優(yōu)化設(shè)計(jì)方法［J］. 南京航空航天大學(xué)學(xué)報(bào)，2011，43（5）：645-649.WU Lili，YAO Weixing. Two-level collaborative opti?mum design method for composite stiffened panel［J］. Journal of Nanjing University of Aeronautics ＆Astronautics，2011，43（5）：645-649.

［2］喬巍，姚衛(wèi)星.復(fù)合材料加筋板鋪層優(yōu)化設(shè)計(jì)的等效彎曲剛度法［J］.計(jì)算力學(xué)學(xué)報(bào)，2011，28（1）：158-162.QIAO Wei，YAO Weixing. Equivalent bending stiff?ness method for stacking sequence optimization of com?posite stiffened panel［J］. Chinese Journal of Computa?tional Mechanics，2011，28（1）：158-162.

［3］趙群，丁運(yùn)亮，金海波. 一種基于復(fù)合材料加筋板結(jié)構(gòu)效率的穩(wěn)定性優(yōu)化方法［J］.復(fù)合材料學(xué)報(bào)，2010，27（3）：169-176.ZHAO Qun，DING Yunliang，JIN Haibo . Buckling optimization method based on structure efficiency of composite stiffened panels［J］. Acta Materiae Composi?tae Sinica，2010，27（3）：169-176.

［4］李爍，徐元銘，張俊. 基于神經(jīng)網(wǎng)絡(luò)響應(yīng)面的復(fù)合材料結(jié)構(gòu)優(yōu)化設(shè)計(jì)［J］. 復(fù)合材料學(xué)報(bào)，2005，22（5）：134-140.LI Shuo，XU Yuanming，ZHANG Jun. Composite structural optimization design based on neural network response surfaces［J］. Acta Materiae Compositae Sini?ca，2005，22（5）：134-140.

［5］張鐵亮，丁運(yùn)亮.復(fù)合材料加筋壁板的結(jié)構(gòu)布局優(yōu)化設(shè)計(jì)［J］. 南京航空航天大學(xué)學(xué)報(bào)，2010，42（1）：8-12.ZHANG Tieliang，DING Yunliang. Structural layout optimization of composite stiffened panel［J］. Journal of Nanjing University of Aeronautics ＆ Astronautics，2010，42（1）：8-12.

［6］劉克龍，姚衛(wèi)星，余熊慶. 運(yùn)用低自由度協(xié)同優(yōu)化的機(jī)翼結(jié)構(gòu)氣動(dòng)多學(xué)科設(shè)計(jì)優(yōu)化［J］.航空學(xué)報(bào)，2007，28（5）：1025-1032.LIU Kelong，YAO Weixing，YU Xiongqing. Multidis?ciplinary structural-aerodynamic design optimization of wings with low degree-of-freedom collaborative opti?mization［J］. Acta Aeronautica et Astronautica Sinica，2007，28（5）：1025-1032.

［7］BLAIR M，HILL S，WEISSHAAR T A. Rapid modeling with innovative structural comcepts. AIAA-98-1755［R］. American Institute of Aeronautics and Astronau?tics，1998.

［8］NAGENDRA S，JESTIN D，GURDAL D，et al. Im?proved genetic algorithm for the design of stiffened composite panels［J］. Computers and Structures，1996，58（3）：543-555.

［9］MOURITZ A P，GELLERT E，BURCHILL P，et al.Review of advanced composite structures for naval ships and submarines［J］. Composite Structures，2001，53（1）：21-42.