黃友欽，王艷萍，傅繼陽，林俊宏，吳玖榮

(廣州大學(xué)廣東省結(jié)構(gòu)安全與健康監(jiān)測工程技術(shù)研究開發(fā)中心，廣東廣州510006)

0 引言

空間桁架結(jié)構(gòu)中構(gòu)件可靈活布置成各種空間外形來滿足建筑需要，因此在現(xiàn)代公共建筑中得到了廣泛應(yīng)用［1］.但此類結(jié)構(gòu)通常具有固有頻率低、柔性大、阻尼小等特點，是典型的風敏感結(jié)構(gòu)，風荷載常成為結(jié)構(gòu)設(shè)計的主要控制荷載［2］.因此，有必要研究這類結(jié)構(gòu)的抗風優(yōu)化設(shè)計方法，使桿件間受力更加協(xié)調(diào)以節(jié)省材料和降低造價.

然而，目前結(jié)構(gòu)抗風優(yōu)化研究基本集中于高層建筑結(jié)構(gòu)，較少涉及大跨屋蓋.Chan等以等效靜力風荷載下高層建筑的樓頂和層間位移響應(yīng)為約束條件，工程造價作為目標函數(shù)，采用最優(yōu)準則法對高層建筑的構(gòu)件尺寸進行優(yōu)化［3］.Spence進一步考慮了場地氣候概率信息，建立了基于可靠度的抗風優(yōu)化數(shù)學(xué)模型［4］.

大跨空間結(jié)構(gòu)的優(yōu)化中通常設(shè)計變量眾多，準則法通常被認為是有效方法［5］.筆者首先闡述基于滿應(yīng)力準則的抗風優(yōu)化設(shè)計理論和程序?qū)崿F(xiàn)，然后介紹用于本研究的雙層柱面網(wǎng)殼及其風荷載獲得方法，最后基于計算結(jié)果討論了截面積下限、滿應(yīng)力步長等的影響.

1 抗風優(yōu)化設(shè)計理論

空間桁架結(jié)構(gòu)的抗風優(yōu)化數(shù)學(xué)模型為:

式中:A={A1，A2，…，An}表示所有n根桿件的截面積;W表示結(jié)構(gòu)重量，ρi和li分別表示桿件i的截面積;σi和［σi］分別表示桿件i的風致應(yīng)力和容許應(yīng)力;Amin表示截面積下限值，即幾何約束下界.

滿應(yīng)力準則法是指通過調(diào)整構(gòu)件的截面尺寸，使各構(gòu)件的應(yīng)力盡量達到容許應(yīng)力，此時認為結(jié)構(gòu)承受荷載的能力被充分發(fā)揮，結(jié)構(gòu)達到最優(yōu)狀態(tài)［6］.空間桁架結(jié)構(gòu)的抗風優(yōu)化可按下面步驟進行.

(1)根據(jù)結(jié)構(gòu)的初始設(shè)計方案確定截面積初值A(chǔ)，且令k=1(k表示迭代次數(shù)).

(2)計算桿件應(yīng)力比.等效靜力風荷載下計算桿件的應(yīng)力，通過與容許應(yīng)力比較來得到應(yīng)力比，并選擇非零桿進行優(yōu)化.

(3)走射線步.由于結(jié)構(gòu)最優(yōu)解位于主約束曲面，所以為了能夠得到最優(yōu)解，應(yīng)使迭代點落在主約束曲面上，因此在兩次應(yīng)力比計算之間增加射線步使所有桿件滿足約束條件［6］.

(4)判斷是否終止迭代.為了判別設(shè)計點是否為最優(yōu)點，可在每次射線步后記錄一次結(jié)構(gòu)重量，當發(fā)現(xiàn)經(jīng)某一射線步后結(jié)構(gòu)重量大于前一次的重量時，就取前一次的設(shè)計點為最優(yōu)點.射線步后通過下式計算結(jié)構(gòu)重量.

(5)走滿應(yīng)力步.走射線步后若計算未收斂則需要將應(yīng)力比還原至原本值，即將各桿件的應(yīng)力同時改變β(k)max倍以得到下次迭代的截面積初值.

當滿應(yīng)力步的步長太大時，可能會錯過真正的最優(yōu)解，而使所求解遠離最優(yōu)解，因此可通過修改應(yīng)力比來縮短滿應(yīng)力步長，使相鄰兩射線步的迭代點之間相距更近，從而提高最優(yōu)解的精度.可通過下式來縮短滿應(yīng)力步長.式中:λ為修改系數(shù)，其值越大表示對滿應(yīng)力步長的改變越小.

(6)越界處理.當某些桿件的截面積小于幾何約束中的截面積下限值時，可由下式進行截面積修正.

進行越界處理后獲得更新的桿件截面積，再回到步驟(2)進行下一次迭代計算，直至結(jié)構(gòu)重量出現(xiàn)增加而停止迭代.

基于MATLAB數(shù)值計算平臺編制優(yōu)化計算程 序 (Wind ResistantOptimization Program，WROP)，整合了等效靜力風荷載、風致響應(yīng)和優(yōu)化計算等三個部分［7］.

2 結(jié)構(gòu)簡介及風荷載

用于分析的空間桁架結(jié)構(gòu)為正放四角錐雙層柱面網(wǎng)殼，高度和跨度分別為為40 m和103 m，長度為140 m，采用縱向邊緣落地支承［8］.結(jié)構(gòu)中桿件的布置如圖1所示，桿件總數(shù)為10 080，節(jié)點總數(shù)為2 592，初始總重量為603.5 t.桿件間的連接視為鉸接，桿件截面積的最大值為1.41e-2 m2，最小值為8.51e-4 m2.桿件的彈性模量和密度分別為206 GPa和7 850 kg/m3.初始設(shè)計的桿件截面積分布如圖2所示，將初始截面積作為優(yōu)化設(shè)計的初值A(chǔ)(1)，將其最小值8.51e-4 m2作為幾何約束下界Amin.

風洞試驗是獲得結(jié)構(gòu)上風荷載的有效方法，該雙層柱面網(wǎng)殼的剛性模型測壓風洞試驗在同濟大學(xué)TJ-2風洞中完成，每隔10°在90°～180°之間進行7個風向角的測試［8］.基于試驗得到的非定常風荷載進行平穩(wěn)激勵隨機振動分析可得到結(jié)構(gòu)上的等效靜力風荷載［9］.

考慮到風向角為150°時結(jié)構(gòu)的風效應(yīng)較為顯著，因此選取150°風向角下的等效靜力風荷載=1.6進行優(yōu)化計算.圖3給出了計算該風向等效靜力風荷載時得到的應(yīng)力峰值桿件的應(yīng)力響應(yīng)功率譜.

圖3 150°風向下桿件的應(yīng)力響應(yīng)功率譜Fig.3 Auto-power spectrum of stress at 150°

3 優(yōu)化結(jié)果與討論

以初始設(shè)計截面對應(yīng)的風致應(yīng)力最大值作為所有桿件的約束容許應(yīng)力，采用編制的計算程序?qū)﹄p層柱面網(wǎng)殼進行抗風優(yōu)化.以下給出優(yōu)化計算結(jié)果并對關(guān)鍵參數(shù)進行討論.

3.1 優(yōu)化計算結(jié)果

優(yōu)化計算76次后迭代停止，圖4給出了結(jié)構(gòu)總重和最大應(yīng)力比隨迭代次數(shù)的變化.可以看出，結(jié)構(gòu)重量約降低48%，最大降幅發(fā)生在首次迭代;最大應(yīng)力比的變化較為平穩(wěn)，由1.18逐漸穩(wěn)定至收斂值.因此，通過設(shè)置射線步保證了所有桿件的應(yīng)力均小于容許應(yīng)力，即不違反強度約束.

3.2 截面積下限的影響

圖5給出了λ=1.0且不設(shè)幾何約束下界時第2次迭代中射線步得到的桿件應(yīng)力比和截面積.由圖可見，最大應(yīng)力比為6.21，遠大于設(shè)定截面積下限時的1.18.同時，由于未設(shè)定截面積下限，滿應(yīng)力步后的最小截面積僅為3.03e-8 m2，小于規(guī)范允許值［10］.

迭代4次后結(jié)構(gòu)重量出現(xiàn)增加，因此取第4次迭代中射線步的結(jié)果作為優(yōu)化結(jié)果.表1給出了優(yōu)化過程中桿件的最大應(yīng)力比、截面積和結(jié)構(gòu)重量的變化，并與設(shè)定幾何約束下界的結(jié)果進行比較.可以看出，不設(shè)截面積下限時，隨著迭代進行，桿件的最小截面積遠小于約束下限值(僅為1.25e-12 m2)，而最大桿件截面積(4.84e-2 m2)卻大于設(shè)定截面積下限時的結(jié)果(1.12e-2 m2)，從而得到的結(jié)構(gòu)總重反而大于設(shè)定截面積下限的結(jié)果.因此，在復(fù)雜桁架結(jié)構(gòu)的抗風優(yōu)化中，雖然射線步已將迭代點拉回主約束曲面，但仍有必要設(shè)定截面積下限來使各桿件間的受力更加協(xié)調(diào)，從而充分利用材料性能.

表1 前3次迭代和停止迭代時的部分計算結(jié)果Tab.1 Results of the first 3 and final iterations

3.3 滿應(yīng)力步長的影響

λ=0.5時，第1次迭代中滿應(yīng)力步后的桿件應(yīng)力比和截面積如圖6所示.顯然，將滿應(yīng)力步長折減后應(yīng)力比的范圍變小，提高了應(yīng)力比的最小值，從而增大了桿件截面積的最小值.

圖7給出了λ=0.5時結(jié)構(gòu)總重隨迭代次數(shù)的變化.對比圖4可見，優(yōu)化得到的結(jié)構(gòu)重量與未修改滿應(yīng)力步長的結(jié)果接近，這是因為設(shè)定截面積下限時桿件的最大應(yīng)力比本身較小，故優(yōu)化結(jié)果對滿應(yīng)力步長縮短并不敏感.由于第2次迭代時最大應(yīng)力比就已十分接近1.0，因此計算收斂所需迭代次數(shù)小于未修改滿應(yīng)力步長的情況.

圖6 第1次迭代中滿應(yīng)力步后的應(yīng)力比和截面積Fig.6 Stress ratio and area after the full stress step in 1st iteration

圖7 結(jié)構(gòu)重量隨迭代次數(shù)的變化Fig.7 Variation in structural weight with iteration

4 結(jié)論

(1)優(yōu)化設(shè)計使網(wǎng)殼總重降低了約48%，所有桿件的應(yīng)力均小于初始設(shè)計下的最大應(yīng)力值，因此在滿足初始設(shè)計要求的前提下使結(jié)構(gòu)總重有效降低;

(2)通過設(shè)定設(shè)計變量下限避免了優(yōu)化得到的桿件截面積小于規(guī)范容許值，且使桿件間受力更加協(xié)調(diào)，充分利用材料性能;

(3)對于最大應(yīng)力比并未明顯超出可行范圍的情況，優(yōu)化結(jié)果對滿應(yīng)力步長并不敏感，但通過修改滿應(yīng)力步長可能使迭代次數(shù)減少.

［1］LU Chun-ling，LI Qiu-sheng，HUANG Sheng-hong，et al.Large eddy simulation of wind effects on a longspan complex roof structure［J］.Journal of Wind Engineering and Industrial Aerodynamics，2012，100(1):1-18.

［2］孫文斌，孫芳錦.大跨度屋蓋風振控制的遺傳算法［J］.鄭州大學(xué)學(xué)報:工學(xué)版，2012，33(1):40-42，50.

［3］CHAN C，HUANG Ming-feng，KWOK K.Stiffness optimization for wind-induced dynamic serviceability design of tall buildings［J］.Journal of Structural Engineering-ASCE.2009，135(8):985-997.

［4］SPENCE S，GIOFFRE M.Large scale reliabilitybased design optimization of wind excited tall buildings［J］.Probabilistic Engineering Mechanics， 2012(28):206-215.

［5］錢令希.工程結(jié)構(gòu)優(yōu)化設(shè)計［M］.北京:科學(xué)出版社，2011.

［6］朱杰江.建筑結(jié)構(gòu)優(yōu)化及應(yīng)用［M］.北京:北京大學(xué)出版社，2011.

［7］BATHE K.Finite element procedures［M］.New Jersey:Prentice-Hall Inc.，1996.

［8］顧明.浙江省華能玉環(huán)電廠干煤棚風效應(yīng)研究［R］.上海:同濟大學(xué)土木工程防災(zāi)國家重點實驗室，2006.

［9］DAVENPORT A.Gust loading factors［J］.ASCEJournal of Structural Division，1967(93):11-34.

［10］中華人民共和國建設(shè)部.JGJ 61—2003網(wǎng)殼結(jié)構(gòu)技術(shù)規(guī)程［S］.北京:中國建筑工業(yè)出版社，2003.