自由曲面形狀和拓?fù)渎?lián)合優(yōu)化研究

馬 騰，趙興忠，2，高博青，吳 慧

（1.浙江大學(xué)空間結(jié)構(gòu)研究中心，浙江杭州310058；2.中國中元國際工程有限公司，北京100089；3.浙江財(cái)經(jīng)大學(xué)城鄉(xiāng)規(guī)劃與管理學(xué)院，浙江杭州310018）

自由曲面形狀和拓?fù)渎?lián)合優(yōu)化研究

馬 騰1，趙興忠1，2，高博青1，吳 慧3

（1.浙江大學(xué)空間結(jié)構(gòu)研究中心，浙江杭州310058；2.中國中元國際工程有限公司，北京100089；3.浙江財(cái)經(jīng)大學(xué)城鄉(xiāng)規(guī)劃與管理學(xué)院，浙江杭州310018）

為了得到形狀與拓?fù)溥@兩大因素共同影響下的自由曲面形狀和拓?fù)涞淖顑?yōu)解，以結(jié)構(gòu)魯棒性為優(yōu)化目標(biāo)，將自由曲面的形狀和拓?fù)渫瑫r(shí)作為優(yōu)化的變量，采用差分演化算法進(jìn)行優(yōu)化求解，實(shí)現(xiàn)了全局最優(yōu)的曲面形狀和網(wǎng)格拓?fù)?算例表明，與形狀和拓?fù)鋬?yōu)二者分開優(yōu)化相比，優(yōu)化得到的該結(jié)構(gòu)魯棒性更高，結(jié)構(gòu)力學(xué)性能更優(yōu)，求得的曲面形狀和網(wǎng)格拓?fù)淠軌蛴行У氐挚顾鶎?yīng)的荷載，體現(xiàn)結(jié)構(gòu)的主要傳力路徑.

自由曲面結(jié)構(gòu)；形狀和拓?fù)洌霍敯粜?；?lián)合優(yōu)化

對于自由曲面結(jié)構(gòu)，影響其力學(xué)性能的有2個(gè)主要因素，分別是曲面的形狀和拓?fù)?初始的建筑設(shè)計(jì)方案，在帶來視覺沖擊的同時(shí)，往往難以兼顧結(jié)構(gòu)的受力性能.為了達(dá)到建筑形態(tài)與力學(xué)性能的協(xié)調(diào)一致，眾多學(xué)者從優(yōu)化角度，探討曲面的合理形狀、桿件的合理拓?fù)?

在自由曲面結(jié)構(gòu)的合理形狀創(chuàng)建方面，崔昌禹等[1]通過桿件節(jié)點(diǎn)坐標(biāo)的移動，逐步減少結(jié)構(gòu)應(yīng)變能，獲得受力性能良好的自由曲面網(wǎng)殼結(jié)構(gòu)形式；武岳等[2]提出曲線變換型自由曲面的形態(tài)創(chuàng)建方法，通過母線沿準(zhǔn)線平移、縮放、旋轉(zhuǎn)，創(chuàng)建復(fù)雜的網(wǎng)格形式，由于曲線的B樣條控制點(diǎn)相較曲面的離散點(diǎn)數(shù)量大幅降低，提高了優(yōu)化效率；Espath等[3]研究以NURBS的控制點(diǎn)坐標(biāo)、權(quán)因子、節(jié)點(diǎn)矢量為優(yōu)化變量的曲面優(yōu)化方法；李娜等[4]根據(jù)仿生原理，進(jìn)行仿扇貝形體的網(wǎng)殼結(jié)構(gòu)創(chuàng)建，并進(jìn)行形態(tài)優(yōu)化.以上研究均以NURBS曲面形狀構(gòu)建作為理論基礎(chǔ)，已成為找形技術(shù)研究和工程應(yīng)用領(lǐng)域的主流方法[5].

自由曲面結(jié)構(gòu)的合理拓?fù)洳贾檬枪こ淘O(shè)計(jì)中的一大難點(diǎn)，一些工程實(shí)例的網(wǎng)格生成仍依賴人工操作和經(jīng)驗(yàn)，如陽光谷工程的網(wǎng)格劃分[6].為了找到結(jié)構(gòu)的合理拓?fù)?，Rasmussen等[7]通過定義節(jié)點(diǎn)與節(jié)點(diǎn)之間是否有桿件連接的方法，探討了析架拓?fù)鋬?yōu)化的離散型方法；連續(xù)體拓?fù)鋬?yōu)化[8]是尋找結(jié)構(gòu)的力學(xué)傳導(dǎo)路徑的優(yōu)化方法，經(jīng)典的方法如均勻化方法、進(jìn)化結(jié)構(gòu)法、變密度法（SIMP）等，均適用于曲面結(jié)構(gòu)的拓?fù)鋬?yōu)化設(shè)計(jì)；SIMP方法[9]是目前廣泛應(yīng)用的連續(xù)體拓?fù)鋬?yōu)化方法，它將一般變密度法的0-1離散問題松弛為[0，l]之間的連續(xù)優(yōu)化問題，降低了求解難度.

自由曲面結(jié)構(gòu)同時(shí)受到形狀與拓?fù)溥@2類變量的影響，而目前的研究多將兩者分開進(jìn)行.尋找最優(yōu)的結(jié)構(gòu)形態(tài)，應(yīng)同時(shí)將形狀和拓?fù)渥兞孔鳛閮?yōu)化變量進(jìn)行優(yōu)化，但目前研究者甚少.Wang等[10]運(yùn)用差分調(diào)整NURBS曲面高度進(jìn)行形狀優(yōu)化，然后在最優(yōu)形狀上運(yùn)用OC法進(jìn)行拓?fù)鋬?yōu)化，因此該方法是在最優(yōu)形狀基礎(chǔ)上的拓?fù)鋬?yōu)化，沒有考慮形狀和拓?fù)渥兞康呐汉闲?Ansola等[11]進(jìn)行自由曲面形狀與拓?fù)渚C合優(yōu)化研究探索，在迭代中，先進(jìn)行一次微小的Coons曲面形狀優(yōu)化；然后在優(yōu)化形狀的基礎(chǔ)上，進(jìn)行一步微小的拓?fù)鋬?yōu)化，每一次優(yōu)化均是在前一步優(yōu)化的基礎(chǔ)上進(jìn)行，如此反復(fù)循環(huán)，雖然算法能夠收斂，但是否達(dá)到全局最優(yōu)有待檢驗(yàn).

本文采用H∞結(jié)構(gòu)魯棒性評價(jià)指標(biāo)[12-13]，該指標(biāo)能夠較好地反映曲面的形狀和拓?fù)浜侠硇?以該指標(biāo)為優(yōu)化目標(biāo)，取NURBS控制點(diǎn)坐標(biāo)作為形狀控制變量，取SIMP相對密度作為拓?fù)鋬?yōu)化變量，運(yùn)用DE算法求解自由曲面的形狀和拓?fù)渎?lián)合優(yōu)化問題.算例表明，DE算法能夠適應(yīng)不同變量同時(shí)存在下的優(yōu)化問題，采用聯(lián)合優(yōu)化設(shè)計(jì)的自由曲面網(wǎng)殼，相比依次進(jìn)行形狀優(yōu)化和拓?fù)鋬?yōu)化，力學(xué)性能更加優(yōu)越.

1 自由曲面形狀拓?fù)渎?lián)合優(yōu)化原理

1.1 NURBS曲面創(chuàng)建原理及SIMP連續(xù)體拓?fù)鋬?yōu)化

NURBS是目前自由曲面創(chuàng)建的主流方法，NURBS參數(shù)化的特點(diǎn)十分適用于曲面的創(chuàng)建與修改，曲面創(chuàng)建過程高效穩(wěn)定，具體可以參見文獻(xiàn)[14].

離散型拓?fù)鋬?yōu)化本質(zhì)是一種0-1離散變量的組合優(yōu)化問題，0代表材料無，1代表材料保留.直接求解0-1優(yōu)化模型有諸多不便，而采用SIMP模型的思想，將材料有無變化變?yōu)榫鶆蛐宰兓蔀檫B續(xù)優(yōu)化問題.

SIMP拓?fù)鋬?yōu)化[15]將單元的相對密度ρ作為設(shè)計(jì)變量，以連續(xù)變量的密度函數(shù)形式顯式地表達(dá)單元相對密度與材料彈性模量之間的對應(yīng)關(guān)系.引入懲罰因子對中間密度進(jìn)行懲罰，使中間密度向0/1兩端聚集，連續(xù)變量的拓?fù)鋬?yōu)化模型逼近0-1離散變量的優(yōu)化模型.對該模型進(jìn)行有限元求解時(shí)，中間變量對應(yīng)一個(gè)很小的彈性模量，對剛度矩陣的影響變得很小.SIMP采用的彈性模量與相對密度的插值關(guān)系如下：

式中：Ee為第e個(gè)單元的插值彈性模量，E0為實(shí)際材料的彈性模量；0＜ρe≤1為材料相對密度；p為懲罰參數(shù)，一般取p≥3，本文取3；n為單元數(shù)目.1.2 聯(lián)合優(yōu)化數(shù)學(xué)模型

進(jìn)行自由曲面形狀拓?fù)渎?lián)合優(yōu)化的數(shù)學(xué)模型如下.

優(yōu)化目標(biāo)是使結(jié)構(gòu)魯棒性指標(biāo)IR最小.設(shè)計(jì)變量ρ=[ρ1，ρ2，…ρn]為單元的材料相對密度向量，Z= [z1，z2，…，zm]為NURBS曲面控制點(diǎn)坐標(biāo)向量，fr為允許保留的材料體積比，Ve和V0分別為第e個(gè)單元的體積和總體積限值.

式（2）的優(yōu)化目標(biāo)IR為基于H∞魯棒控制理論建立的結(jié)構(gòu)魯棒性評價(jià)指標(biāo).將結(jié)構(gòu)視為一個(gè)輸入輸出系統(tǒng)，則結(jié)構(gòu)系統(tǒng)傳遞函數(shù)的H∞范數(shù)即IR.IR可以衡量結(jié)構(gòu)系統(tǒng)輸出信號與輸入信號是否相稱，反映結(jié)構(gòu)對各種干擾的敏感性，體現(xiàn)結(jié)構(gòu)的健壯程度，即魯棒性，關(guān)于H∞結(jié)構(gòu)魯棒性理論的詳細(xì)研究參考文獻(xiàn)[16]的研究.

式（2）的求解采用差分演化算法（differential evolution）[17].該方法是一種有效的全局最優(yōu)解的搜索算法.對于式（2）的不等式約束，本文采用引入懲罰函數(shù)逐步逼近方法.

對于自由曲面結(jié)構(gòu)，本文采用SIMP法得出以相對剛度為連續(xù)變量的拓?fù)鋬?yōu)化模型，并按式（2）的優(yōu)化目標(biāo)，采用DE算法搜索得到最優(yōu)解.

2 算例分析

2.1 聯(lián)合優(yōu)化與單獨(dú)形狀、拓?fù)鋬?yōu)化對比

通過對比聯(lián)合優(yōu)化結(jié)果與傳統(tǒng)的單獨(dú)形狀、拓?fù)鋬?yōu)化結(jié)果，驗(yàn)證聯(lián)合優(yōu)化的合理性.

本例的結(jié)構(gòu)模型為如圖1所示的一單位正方形平面，施加均布荷載q=11k N/m2，求解結(jié)構(gòu)上凸時(shí)的最優(yōu)形狀和拓?fù)?曲面控制點(diǎn)分布、網(wǎng)格數(shù)量20×20分別如圖2、3所示，fr=0.5.基于結(jié)構(gòu)對稱性，控制點(diǎn)取2、3、6、7、8、11、12、13共8個(gè)，其余控制坐標(biāo)由對稱取得，在1/4結(jié)構(gòu)單元中引入SIMP變密度，剩余單元SIMP變密度由對稱取得.

圖1 單位正方形結(jié)構(gòu)模型Fig.1 Unit square structure model

圖2 曲面控制點(diǎn)分布Fig.2 Distribution of curve control points

圖3 有限元網(wǎng)格劃分Fig.3 Meshing of finite element

本例共進(jìn)行兩類結(jié)構(gòu)優(yōu)化：1）單獨(dú)形狀、拓?fù)鋬?yōu)化，即先對結(jié)構(gòu)進(jìn)行單獨(dú)的形狀優(yōu)化，取所有的單元相對密度ρi=1，在此優(yōu)化基礎(chǔ)上進(jìn)行單獨(dú)的拓?fù)鋬?yōu)化，如圖4所示為最終的優(yōu)化結(jié)果；2）本文提出的形狀和拓?fù)涞穆?lián)合優(yōu)化，兩類變量同時(shí)進(jìn)行優(yōu)化求解，約束條件同第一類.優(yōu)化結(jié)果如圖5所示.

優(yōu)化結(jié)果的形狀對比：以0.3 m等高線作為對比，采用單獨(dú)優(yōu)化的曲面中部突起范圍更廣、更加平坦；兩者最終矢高均達(dá)到0.5 m.

優(yōu)化結(jié)果的拓?fù)鋵Ρ龋簝烧叩淖罱K拓?fù)錁?gòu)型較接近，略有不同，但聯(lián)合優(yōu)化結(jié)果較清晰；拓?fù)錁?gòu)型顯示在均布荷載作用下，結(jié)構(gòu)主要通過對角線方向受力；支座部位顏色較深，表明此處對整個(gè)結(jié)構(gòu)的魯棒性影響較大，應(yīng)重點(diǎn)加強(qiáng).

聯(lián)合優(yōu)化與單獨(dú)優(yōu)化的數(shù)據(jù)，列于表1中.表中，IR為魯棒性指標(biāo)，U為應(yīng)變能，D為最大位移.

本文采用的優(yōu)化目標(biāo)為結(jié)構(gòu)魯棒性指標(biāo)，單獨(dú)優(yōu)化時(shí)先進(jìn)行形狀優(yōu)化，然后在形狀優(yōu)化基礎(chǔ)上進(jìn)行拓?fù)鋬?yōu)化，得到的魯棒性指標(biāo)為5.33×10-3，而聯(lián)合優(yōu)化的魯棒性指標(biāo)為4.26×10-3，明顯優(yōu)于前述方法，證明了聯(lián)合優(yōu)化對尋找最優(yōu)結(jié)構(gòu)的有效性.結(jié)構(gòu)應(yīng)變能和最大位移均出現(xiàn)了下降，表明在結(jié)構(gòu)魯棒性提高的同時(shí)，結(jié)構(gòu)的剛度變大.

圖5 聯(lián)合優(yōu)化結(jié)果Fig.5 Results of combined optimization

表1 聯(lián)合優(yōu)化與單獨(dú)拓?fù)鋬?yōu)化結(jié)果Tab.1 Results of combined optimization and topology optimization

先進(jìn)行形狀優(yōu)化再進(jìn)行拓?fù)鋬?yōu)化的方法得到的結(jié)構(gòu)，指標(biāo)劣于聯(lián)合優(yōu)化的原因如下：在單獨(dú)進(jìn)行拓?fù)鋬?yōu)化時(shí)，結(jié)構(gòu)形狀已經(jīng)確定，故在進(jìn)行拓?fù)鋬?yōu)化時(shí)，形狀變量無法改變，不能實(shí)現(xiàn)形狀、拓?fù)渥兞康呐汉嫌绊懀顑?yōu)拓?fù)涫艿较葲Q形狀的限制.之所以形狀和拓?fù)渥兞恳?lián)合優(yōu)化，是因?yàn)閮烧邔?yōu)化指標(biāo)的影響不是線性獨(dú)立的，而是非線性藕合的；在尋求最優(yōu)結(jié)構(gòu)時(shí)，應(yīng)采用聯(lián)合優(yōu)化的方法，這樣才能求解出最優(yōu)的形狀和拓?fù)?

根據(jù)均布荷載作用下聯(lián)合優(yōu)化得到的構(gòu)型，可以設(shè)計(jì)如圖6所示的曲面（更準(zhǔn)確地理解為結(jié)構(gòu)的最優(yōu)傳力路徑）.可以看出，結(jié)構(gòu)總是試圖以最短路徑傳遞荷載，即一個(gè)性狀良好的結(jié)構(gòu)，傳力方式總是簡單直接的，這與結(jié)構(gòu)設(shè)計(jì)的認(rèn)知是一致的；同時(shí)，傳力路徑呈拱形，且傳遞路徑可以分為對角拱和邊界拱，將均布荷載收集傳遞至支座.

圖6 全跨荷載結(jié)構(gòu)主要傳力路徑Fig.6 Main force travelling path of structure under allcross load

2.2 曲面形狀和拓?fù)渎?lián)合優(yōu)化算例分析

2.2.1 半跨荷載作用下角點(diǎn)支承的自由曲面聯(lián)合優(yōu)化 為了比較荷載形式對優(yōu)化結(jié)果的影響，將圖1改為施加左半邊的均布半跨荷載，由于半跨荷載具有1/2對稱性，控制點(diǎn)取圖4中的2、3、6、7、8、11、12、13、16、17、18、22、23共13個(gè)形狀變量，在圖5下半部共10×20=200個(gè)單元中引入SIMP變密度模型，其余條件參考例2.1.DE參數(shù)如下：種群NP=500，優(yōu)化步Iter=400.

圖7 優(yōu)化拓?fù)渫敢晥DFig.7 Optimized topology perspective

最優(yōu)的結(jié)構(gòu)形狀和材料布局如圖7、8所示.在半跨荷載施加的位置，結(jié)構(gòu)上凸，最高點(diǎn)達(dá)到0.461，而無荷載的半跨，結(jié)構(gòu)較低緩，這是對半跨荷載最有效的抵抗形狀；變密度值的分布集中于左半跨，除了支座部位，右半跨材料較稀疏，單元顏色較深的地方用粗線繪制，單元顏色較淺的地方用細(xì)線繪制.可以發(fā)現(xiàn)，從荷載中心附近一點(diǎn)發(fā)出，與支座相連的4條粗線和3條邊界線，是結(jié)構(gòu)的主要傳力路徑，是材料分布的集中區(qū)域.

圖8 半跨荷載結(jié)構(gòu)主要傳力路徑Fig.8 Main force travelling path of structure under semi-cross load

2.2.2 集中荷載作用下四邊中點(diǎn)支承的自由曲面聯(lián)合優(yōu)化 將平面中心位置作用5k N的集中荷載，正方形平面四周邊中點(diǎn)支承，由結(jié)構(gòu)1/4對稱性，取圖4中1、2、6、7、8、12、13共7個(gè)控制點(diǎn)作為形狀變量，1/4結(jié)構(gòu)單元的SIMP值作為拓?fù)渥兞?，共有拓?fù)渥兞?00個(gè)，其余條件參考例2.1.DE算法參數(shù)如下：NP=250，Iter=300.

優(yōu)化完成后，曲面形狀如圖9所示，該曲面的邊界保持不變，中心處形成突起，最高點(diǎn)達(dá)到0.37；曲面的拓?fù)錁?gòu)型如圖10所示.可以發(fā)現(xiàn)，荷載傳力路徑呈十字型，材料分布集中于荷載作用處以及4個(gè)支座位置.

圖9 優(yōu)化后曲面形狀Fig.9 Surface after optimization

3 結(jié) 論

（1）聯(lián)合優(yōu)化與傳統(tǒng)單獨(dú)優(yōu)化的對比結(jié)果表明，形狀和拓?fù)涔餐绊懽杂汕娼Y(jié)構(gòu)的力學(xué)性能，兩者呈非線性藕合關(guān)系.優(yōu)化數(shù)值證明了聯(lián)合優(yōu)化后的結(jié)構(gòu)是數(shù)學(xué)上的全局最優(yōu)解，而單獨(dú)優(yōu)化只是特定條件下的局部最優(yōu)解.

圖10 優(yōu)化后拓?fù)銯ig.10 Topology after optimization

（2）本文的聯(lián)合優(yōu)化結(jié)果為自由曲面網(wǎng)格結(jié)構(gòu)的找形和桿件布置提供了依據(jù).盡管拓?fù)錁?gòu)型不能直接用于實(shí)際網(wǎng)格劃分，但為網(wǎng)格布置提供了理論依據(jù)與參考，可以指導(dǎo)自由曲面的網(wǎng)格劃分，明確桿件布置的方向，以滿足力流傳導(dǎo)的需要.

（3）本文所采用的H∞結(jié)構(gòu)魯棒性指標(biāo)是描述結(jié)構(gòu)力學(xué)性能的綜合指標(biāo)，反映結(jié)構(gòu)抵抗破壞的能力.本文的算例表明，在結(jié)構(gòu)魯棒性提高的同時(shí)，結(jié)構(gòu)的應(yīng)變能下降，即結(jié)構(gòu)的剛度提高.

（

）：

[1]崔昌禹，王有寶，姜寶石，等.自由曲面單層網(wǎng)殼結(jié)構(gòu)形態(tài)創(chuàng)構(gòu)方法研究[J].土木工程學(xué)報(bào)，2013，46（4）：57-63.

CUI Chang-yu，WANG You-bao，JIANG Bao-shi，et al.Study on the structural morphogenesis technique for single-layer reticulated shells of free-curved surface[J].China Civil Engineering Journal，2013，46（4）：57-63.

[2]武岳，李欣，曹正罡.曲線變換型自由曲面結(jié)構(gòu)的形態(tài)創(chuàng)建[J].建筑結(jié)構(gòu)學(xué)報(bào)，2012，33（5）：23-30.

WU Yue，LI Xin，CAO Zheng-gang.Computational morphogenesis method of curve-generated freeform structures[J].Journal of Building Structures，2012，33（5）：23-30.

[3]ESPATH L F R，LINN R V，AWRUCH A M.Shape optimization of shell structures based on NURBS description using automatic differentiation[J].International Journal for Numerical Methods in Engineering，2011，88（7）：613-636.

[4]李娜，岑培超，羅堯治.仿扇貝形體的網(wǎng)殼結(jié)構(gòu)建模及造型優(yōu)化方法[J].科技通報(bào)，2008，24（2）：224-230.

LI Na，CEN Pei-chao，LUO Yao-zhi.Modeling methodof the scallop-shaped lattice shell and optimization with strain energy method[J].Bulletin of Science and Technology，2008，24（2）：224-230.

[5]張崢，丁潔民.大跨度結(jié)構(gòu)的復(fù)雜曲面找形技術(shù)及工程應(yīng)用[J].建筑結(jié)構(gòu)，2006，36（3）：5-8.

ZHANG Zheng，DING Jie-min.Form-finding technique and application of complex curved surfaces in large-span structures[J].Journal of Building Structures，2006，36（3）：5-8.

[6]李承銘，盧旦.自由曲面單層網(wǎng)格的智能布局設(shè)計(jì)研究[J].土木工程學(xué)報(bào)，2011，44（3）：1-7.

LI Cheng-ming，LU Dan.Study of intelligent layout design of single-layer lattice shell of free form surface[J].China Civil Engineering Journal，2011，44（3）：1-7.

[7]RASMUSSEN M H，STOLPE M.Global optimization of discrete truss topology design problems using a parallel cut-and-branch method[J].Computers and Structures，2008，86（13）：1527-1538.

[8]羅震，陳立平，黃玉盈，等.連續(xù)體結(jié)構(gòu)的拓?fù)鋬?yōu)化設(shè)計(jì)[J].力學(xué)進(jìn)展，2004，34（4）：463-476.

LUO Zhen，CHEN Li-ping，HUANG Yu-ying，et al.Topological optimization design for continuum structures [J].Advances in Mechanics，2004，34（4）：463-476.

[9]ROZVANY G.A critical review of established methods of structural topology optimization[J].Structural and Multidisciplinary Optimization，2009，37（3）：217-237.

[10]WANG Lei，ZHANG Qi-lin，YANG Bin.Combined shape and topology optimization of free form shells[J].Advanced Materials Research，2011，163：2356-2360.

[11]ANSOLA R，CANALES J，TáRRAGO J A，et al.An integrated approach for shape and topology optimization of shell structures[J].Computers and Structures，2002，80（5）：449-458.

[12]DOYLE J C，GLOVER K，KHARGONEKAR P P，et al.State-space solutions to standard H2and H∞control problems[J].IEEE Transactions on Automatic Control，1989，34（8）：831-847.

[13]梅生偉，申鐵龍，劉康志.現(xiàn)代魯棒控制理論與應(yīng)用[M].北京：清華大學(xué)出版社，2008：61-63.

[14]PIEGL L，TILLER W.非均勻有理B樣條[M].趙罡，穆國旺，王拉柱，譯.2版.北京：清華大學(xué)出版社，2010：86-101.

[15]昌俊康，段寶巖.連續(xù)體結(jié)構(gòu)拓?fù)鋬?yōu)化的一種改進(jìn)變密度法及其應(yīng)用[J].計(jì)算力學(xué)學(xué)報(bào)，2009，26（2）：188-192.

CHANG Jun-kang，DUAN Bao-yan.An improved variable density method and application for topology optimization of continuum structures[J].Chinese Journal of Computational Mechanics，2009，26（2）：188-192.

[16]張成，吳慧，高博青，等.基于H∞理論的結(jié)構(gòu)魯棒性分析[J].建筑結(jié)構(gòu)學(xué)報(bào)，2012，33（5）：87-92.

ZHANG Cheng，WU Hui，GAO Bo-qing，et al.H∞-based structural robustness analysis[J].Journal of Building Structures，2012，33（5）：87-92.

[17]STORN R，PRICE K.Differential evolution：a simple and efficient heuristic for global optimization over continuous spaces[J].Journal of Global Optimization，1997，11（4）：341-359.

Combined shape and topology optimization of free-form structure

MA Teng1，ZHAO Xing-zhong1，2，GAO Bo-qing1，WU Hui3

（1.Space Structures Research Center，Zhejiang University，Hangzhou 310058，China；2.China IPPR International Engineering Corporation，Beijing100089，China；3.College of Urban/Rural Planning and Management，Zhejiang University of Finance and Economics，Hangzhou 310018，China）

Shape and topology of free-form grid structures were simultaneously considered in order to reach the global optimum shapes and topologies.Optimum structural robustness was set to be the objective while parameters to define the structural shape and topology were taken as design variables.A differential evolution algorithm was used to solve the combined optimization problem.Global optimum shapes and topologies of free-form structures were obtained.The obtained results were compared with those achieved by optimizing solely shape or topology.Results showed that more robust structures with better mechanical properties were obtained，which implied the superiority of the method.The optimized structures were more efficient in carrying the load，and the resulted shape and topology clearly revealed the main force transferring paths.

free-form structure；shape and topology；robustness；combined optimization

TU 393

A

1008-973X（2015）10-1946-06

2014-08-22.浙江大學(xué)學(xué)報(bào)（工學(xué)版）網(wǎng)址：www.journals.zju.edu.cn/eng

國家自然科學(xué)基金面上項(xiàng)目（51378457）；高等學(xué)校博士學(xué)科點(diǎn)專項(xiàng)科研基金資助項(xiàng)目（20120101110023）；浙江省自然科學(xué)基金資助項(xiàng)目（Y15E08004）.

馬騰（1988—），男，碩士，從事大跨度空間結(jié)構(gòu)研究.E-mail：mt001mt@163.com

高博青，男，教授.E-mail：bqgao@zju.edu.cn