陳世英, 李維國(guó), 鹿曉陽

(1.中國(guó)石油大學(xué)(華東)儲(chǔ)運(yùn)與建筑工程學(xué)院,山東青島 266580; 2.山東建筑大學(xué)土木工程學(xué)院,山東濟(jì)南 250101)

網(wǎng)殼結(jié)構(gòu)截面優(yōu)化設(shè)計(jì)的思想為離散變量?jī)?yōu)化設(shè)計(jì),以桿件截面尺寸為設(shè)計(jì)變量,以結(jié)構(gòu)質(zhì)量等為目標(biāo)函數(shù),結(jié)構(gòu)的強(qiáng)度、剛度、局部(桿件)穩(wěn)定性和整體穩(wěn)定性為約束條件,尋求滿足約束條件且目標(biāo)函數(shù)最小的桿件截面取值[1-3]?，F(xiàn)行的空間網(wǎng)格結(jié)構(gòu)技術(shù)規(guī)程[4]規(guī)定單層網(wǎng)殼結(jié)構(gòu)必須進(jìn)行整體穩(wěn)定性驗(yàn)算,其彈塑性穩(wěn)定承載力須大于兩倍的均布外荷載。整體穩(wěn)定性是單層網(wǎng)殼結(jié)構(gòu)設(shè)計(jì)的關(guān)鍵問題[5-8],在優(yōu)化設(shè)計(jì)過程中必須將整體穩(wěn)定承載力作為約束條件充分考慮[9-10],否則無法保證優(yōu)化后結(jié)構(gòu)的穩(wěn)定承載力安全。單層網(wǎng)殼結(jié)構(gòu)是缺陷敏感結(jié)構(gòu),整體穩(wěn)定承載力隨缺陷的不同而變化顯著[11-14],沈世釗等[5,15]試圖尋找最不利缺陷以獲得最不利的彈塑性穩(wěn)定承載力,但準(zhǔn)確找到結(jié)構(gòu)的最不利缺陷分布幾乎不可能[16]。長(zhǎng)期以來,均采用一致缺陷法或隨機(jī)缺陷法來解決結(jié)構(gòu)的初始幾何缺陷問題[17-19]。在通常的網(wǎng)殼結(jié)構(gòu)截面優(yōu)化設(shè)計(jì)“桿件分組”方式中,會(huì)將網(wǎng)殼結(jié)構(gòu)桿件按照空間位置的對(duì)稱性分成若干組[20-21],記做空間位置分組法。當(dāng)整體穩(wěn)定性約束條件和初始幾何缺陷同時(shí)引入到截面優(yōu)化設(shè)計(jì)中,通過分析發(fā)現(xiàn)[22],基于空間位置分組法優(yōu)化后的結(jié)構(gòu)僅在設(shè)定缺陷下其穩(wěn)定承載力滿足設(shè)計(jì)要求,在其他可能出現(xiàn)的缺陷下往往不能滿足?？梢?基于空間位置分組法的優(yōu)化設(shè)計(jì),使網(wǎng)殼結(jié)構(gòu)的穩(wěn)定承載力對(duì)缺陷分布的敏感性增加,無法保證優(yōu)化后結(jié)構(gòu)的穩(wěn)定承載力安全。筆者結(jié)合網(wǎng)殼結(jié)構(gòu)的截面優(yōu)化設(shè)計(jì)原理提出基于等穩(wěn)定要求的單層球面網(wǎng)殼結(jié)構(gòu)截面優(yōu)化設(shè)計(jì)方法,在充分考慮網(wǎng)殼結(jié)構(gòu)桿件強(qiáng)度、剛度、雙重非線性整體穩(wěn)定性和隨機(jī)初始缺陷條件下,對(duì)跨度為60 m的K6網(wǎng)殼結(jié)構(gòu)按缺陷敏感區(qū)域進(jìn)行桿件分組并進(jìn)行基于序列兩級(jí)優(yōu)化算法的截面優(yōu)化設(shè)計(jì)。

1 基于等穩(wěn)定要求截面優(yōu)化設(shè)計(jì)數(shù)學(xué)模型

設(shè)計(jì)變量為

x={A,Iy,Iz,Wy,Wz,ρ}.

式中,x對(duì)應(yīng)一種型鋼的截面屬性;A為截面積,m2;Iy、Iz分別為對(duì)y、z軸的慣性矩,m4;Wy、Wz分別為對(duì)y、z軸的抗彎截面模量,m3;ρ為線密度,kg/m。

桿件分組:按網(wǎng)殼結(jié)構(gòu)的缺陷敏感區(qū)域?qū)U件進(jìn)行分組。

目標(biāo)函數(shù):同等用鋼量情況下網(wǎng)殼結(jié)構(gòu)穩(wěn)定承載力隨缺陷波動(dòng)最小,表示為

minD=D(Pcr1,Pcr2,…,Pcrl) .

(1)

式中,pcri為第i個(gè)隨機(jī)缺陷下結(jié)構(gòu)的穩(wěn)定承載力,kN/m2;l為隨機(jī)缺陷個(gè)數(shù);D為l個(gè)穩(wěn)定承載力的標(biāo)準(zhǔn)差。

(2)

式中,W為網(wǎng)殼的用鋼量;m1為桿件的分組數(shù);ρk為第k組桿件的線密度;Ak為第k組桿件的橫截面積;lk為第k組桿件的總長(zhǎng)度,m。

應(yīng)力約束[23]為

(3)

式中,N為桿件軸力,kN;My和Mz分別為繞y軸和z軸的彎矩,kN·m;Wy和Wz分別為兩主軸的截面抗彎模量,m3;γy和γz為與截面模量相應(yīng)的截面塑性發(fā)展系數(shù),均取1.15;f為桿件的屈服應(yīng)力,kPa。

彎矩作用下平面內(nèi)穩(wěn)定性約束為

(4)

其中

彎矩作用下平面外穩(wěn)定約束為

(5)

式中,φz為彎矩作用平面外的軸心受壓構(gòu)件穩(wěn)定系數(shù);φb為均勻彎曲的受彎構(gòu)件整體穩(wěn)定系數(shù),取1.0;η為截面影響系數(shù),閉口截面η=0.7;βty為等效彎矩系數(shù),取1.0。

長(zhǎng)細(xì)比約束[4]為

λ≤[λ].

(6)

式中,[λ]為桿件的許用長(zhǎng)細(xì)比,壓桿取150,拉桿取250。

最大位移約束[4]為

Δmax≤ln/400 .

(7)

整體穩(wěn)定性約束[4]為

Pcr≤2Q.

(8)

式中,Δmax為網(wǎng)殼結(jié)構(gòu)最大節(jié)點(diǎn)位移,m;ln為網(wǎng)殼結(jié)構(gòu)跨度,m;Q為網(wǎng)殼結(jié)構(gòu)均布外載荷,kN/m2;Pcr為網(wǎng)殼結(jié)構(gòu)彈塑性穩(wěn)定承載力,kN。

2 基于離散變量的序列兩級(jí)截面優(yōu)化方法

約束條件按其性質(zhì)可分為局部性與整體性約束兩類。在截面優(yōu)化過程中求解應(yīng)力、長(zhǎng)細(xì)比等局部性約束條件時(shí)采用一維搜索方法,得到滿足局部約束條件的桿件截面尺寸,即第一級(jí)優(yōu)化。整體性約束條件含有各個(gè)單元的設(shè)計(jì)變量,采用相對(duì)差商法[24]進(jìn)行整體尋優(yōu),以第一級(jí)優(yōu)化的結(jié)果作為第二級(jí)優(yōu)化的初始值,反復(fù)迭代直至整體性約束條件滿足。進(jìn)行0階修正可改進(jìn)該算法的結(jié)果精度。

采用APDL語言編寫結(jié)構(gòu)建模程序、缺陷施加程序和進(jìn)行彈塑性大變形分析。桿件模型采用beam189單元,屈服準(zhǔn)則采用von Mise準(zhǔn)則[25],材料本構(gòu)關(guān)系為理想彈塑性。用FORTRAN編寫基于離散變量的序列兩級(jí)算法截面優(yōu)化程序,通過接口程序完成FORTRAN與ANSYS的相互調(diào)用和數(shù)據(jù)交換。桿件采用熱軋無縫鋼管Q235,從型鋼表[26]中以抗彎截面系數(shù)的升序排列為原則共提取63種型鋼組成桿件截面離散集,依次編號(hào)為1#～63#,1#截面取Φ48 mm×3 mm。文中使用的桿件截面編號(hào)和對(duì)應(yīng)尺寸包括：23#/Φ108 mm×4.0 mm，30#/Φ146 mm×4.5 mm，31#/Φ152 mm×4.5 mm，32#/Φ159 mm×4.5 mm，33#/Φ168 mm×4.5 mm，34#/Φ168 mm×5.0 mm，35#/Φ180 mm×5.0 mm。

3 結(jié)構(gòu)缺陷敏感區(qū)域分析

網(wǎng)殼結(jié)構(gòu)是缺陷敏感結(jié)構(gòu),許多學(xué)者對(duì)其敏感區(qū)域進(jìn)行了分析[27-29]。從缺陷與穩(wěn)定承載力波動(dòng)關(guān)系的角度提出網(wǎng)殼結(jié)構(gòu)的穩(wěn)定承載力敏感區(qū)域,并將敏感區(qū)域與桿件分組相結(jié)合,提出了一種新的截面優(yōu)化分組方法。

算例1:K6型單層凱威特網(wǎng)殼,跨度為60 m,矢高為15 m,環(huán)向桿件圈數(shù)為10圈;當(dāng)某一個(gè)節(jié)點(diǎn)有缺陷時(shí),分析結(jié)構(gòu)穩(wěn)定承載力值的變化規(guī)律。由于網(wǎng)殼結(jié)構(gòu)的對(duì)稱性,可僅計(jì)算缺陷出現(xiàn)在圖1所示節(jié)點(diǎn)上的情況,缺陷值取跨度的1/300。單個(gè)節(jié)點(diǎn)存在缺陷時(shí)結(jié)構(gòu)的穩(wěn)定承載力變化規(guī)律見圖2。

圖1 1/6結(jié)構(gòu)節(jié)點(diǎn)編號(hào)Fig.1 Node number of one-sixth structure

圖2 單個(gè)節(jié)點(diǎn)有缺陷時(shí)結(jié)構(gòu)的穩(wěn)定承載力Fig.2 Critical-limit load under single joint imperfection

由圖2可以看出,相比于凸缺陷,節(jié)點(diǎn)發(fā)生凹缺陷時(shí),穩(wěn)定承載力更低。其中,1#～94#節(jié)點(diǎn)上存在凹缺陷時(shí),承載力低于4.5 kN/m2,94#～175#節(jié)點(diǎn)上存在凹缺陷時(shí),承載力為4.5～5.5 kN/m2。175#～223#節(jié)點(diǎn)上存在凹缺陷時(shí),承載力為5.5～6.5 kN/m2。272#～278#節(jié)點(diǎn)(支座節(jié)點(diǎn))上存在凹缺陷時(shí),承載力為6.5～7.5 kN/m2?？梢?節(jié)點(diǎn)缺陷所處位置在殼體中上部時(shí)對(duì)結(jié)構(gòu)整體穩(wěn)定性影響較大;節(jié)點(diǎn)缺陷在殼體中部偏下時(shí)對(duì)結(jié)構(gòu)整體穩(wěn)定性影響較小。支座節(jié)點(diǎn)缺陷對(duì)結(jié)構(gòu)穩(wěn)定性能影響較小。

定義相同緯度的節(jié)點(diǎn)為一組,網(wǎng)殼節(jié)點(diǎn)可分為10組,緯度最高的一組節(jié)點(diǎn)定義為第一環(huán)節(jié)點(diǎn)(圖1中,2#、3#節(jié)點(diǎn)所在的環(huán)),緯度最低一組節(jié)點(diǎn)定義為第10環(huán)節(jié)點(diǎn)(圖1中272#～282#節(jié)點(diǎn)所在的環(huán))。缺陷出現(xiàn)在各環(huán)節(jié)點(diǎn)上時(shí)結(jié)構(gòu)的穩(wěn)定承載力變化曲線見圖3。

圖3 十種環(huán)向缺陷下結(jié)構(gòu)的穩(wěn)定承載力Fig.3 Critical-limit load of structure with 10 hoop imperfections

由圖3可以看出,相比于凸缺陷節(jié)點(diǎn)存在凹缺陷時(shí),絕大多數(shù)情況下穩(wěn)定承載力更低?？紤]凹缺陷的情況,當(dāng)殼體中上部節(jié)點(diǎn)整環(huán)出現(xiàn)缺陷時(shí)(1～7環(huán)),結(jié)構(gòu)的穩(wěn)定承載力為3.0～4.0 kN/m2,當(dāng)殼體靠近支座處節(jié)點(diǎn)整環(huán)出現(xiàn)缺陷時(shí)(8～9環(huán)),結(jié)構(gòu)的穩(wěn)定承載力為4.0～5.0 kN/m2,當(dāng)支座位置的環(huán)向節(jié)點(diǎn)(10環(huán))發(fā)生缺陷時(shí),對(duì)結(jié)構(gòu)的整體穩(wěn)定性承載力影響很小。

4 基于等穩(wěn)定要求的網(wǎng)殼結(jié)構(gòu)截面優(yōu)化設(shè)計(jì)

網(wǎng)殼結(jié)構(gòu)存在缺陷敏感區(qū)域,且敏感程度不同。提出了基于敏感區(qū)域的桿件分組方法,將缺陷敏感性相近區(qū)域的桿件分為一組。

算例2: 結(jié)構(gòu)模型尺寸與算例1相同,外載荷Q為1.9 kN/m2,滿跨均勻布置。基于結(jié)構(gòu)敏感區(qū)域?qū)U件共分為3組(圖4)。縱向所有桿件為第1組;環(huán)向桿件分為兩組,1～7環(huán)桿件為第2組,8～10環(huán)桿件為第3組;基于空間位置的桿件分組見圖5,桿件共分為20組。隨機(jī)生成28種均值為0,均方差為R/3的正態(tài)分布隨機(jī)缺陷(所有節(jié)點(diǎn)上均有缺陷),最大缺陷值R取跨度的1/300。

當(dāng)所有桿件取相同桿件截面時(shí)(桿件不分組),在28種缺陷下網(wǎng)殼結(jié)構(gòu)的穩(wěn)定承載力(桿件截面為Φ168 mm×4.5 mm(33#))分別為5.002 4、4.064 7、4.568 7、4.596 1、3.962 1、4.459 3、4.461 8、4.801 9、3.820 4、4.298 9、4.475 8、4.596 9、4.216 4、4.138 9、4.387 5、3.758 1、4.693 7、4.391 4、4.787 2、4.297 9、4.393 4、4.823 0、4.795 7、4.724 1、4.398 9、4.418 5、4.849 4、4.734 1 kN/m2。當(dāng)桿件取33#截面時(shí),在第16組缺陷下穩(wěn)定承載力不滿足約束條件(3.758 1<2Q)。為提高第16組缺陷下結(jié)構(gòu)的穩(wěn)定承載力,全部桿件增大為34#截面,此時(shí)結(jié)構(gòu)的耗鋼量增大約7 t。

為對(duì)比分析未優(yōu)化結(jié)構(gòu)和優(yōu)化結(jié)構(gòu)的穩(wěn)定承載力對(duì)缺陷的敏感性,按照缺陷敏感區(qū)域分組方法(圖4)和空間位置分組方法(圖5)對(duì)網(wǎng)殼結(jié)構(gòu)進(jìn)行截面優(yōu)化設(shè)計(jì),取第16#隨機(jī)缺陷為設(shè)計(jì)缺陷。

圖4 基于敏感區(qū)域的桿件分組Fig.4 Group based on sensitive area

圖5 基于空間位置的桿件分組(標(biāo)號(hào)相同桿件同組)Fig.5 Group based on space position(bars grouped together with the same number)

采用空間位置分組方法優(yōu)化后，桿件的截面取值編號(hào)為(1～20組)：24#、 24#、29#、24#、33#、33#、33#、38#、33#、24#、31#、24#、31#、23#、30#、23#、30#、21#、30#、21#，在28種缺陷下的穩(wěn)定承載力分別為2.895 7、2.452 7、3.529 7、2.473 7、2.723 6、2.483 9、3.185 3、2.893 2、3.223 3、2.743 3、2.497 3、2.795 3、2.182 2、2.915 9、2.706 4、3.896 9、2.542 6、2.656 1、2.976 8、2.977 6、2.736 8、3.458 7、3.261、2.510 8、2.458 3、3.173 3、3.162 3、3.275 3 kN/m2。

采用等穩(wěn)定方法優(yōu)化后桿件截面取值編號(hào)為(1～3組)：33#、35#、23#，在28種缺陷下的穩(wěn)定承載力分別為4.891 2、4.949 9、5.278 1、4.694 4、4.372 3、4.805 6、4.599 5、4.705 1、4.242 6、4.966 8、4.776 8、4.427 6、4.531 2、4.491 3、4.732 9、4.107 1、5.211 4、4.488 1、5.198 2、4.431 3、4.769 7、5.154 1、4.853 3、4.949 4、4.431 1、4.688 3、4.304 3、5.125 9 kN/m2。

采用空間位置分組方法優(yōu)化得到的結(jié)構(gòu)僅在設(shè)計(jì)缺陷(16#缺陷)下滿足穩(wěn)定承載力設(shè)計(jì)要求，在其他27種隨機(jī)缺陷下，穩(wěn)定承載力均小于2Q，不滿足整體穩(wěn)定承載力設(shè)計(jì)要求。采用基于缺陷敏感區(qū)域分組方式優(yōu)化得到的等穩(wěn)定結(jié)構(gòu)模型，在28種隨機(jī)缺陷下穩(wěn)定承載力均滿足設(shè)計(jì)要求。

未經(jīng)截面優(yōu)化、空間位置分組方法優(yōu)化和等穩(wěn)定方法優(yōu)化得出的3種結(jié)構(gòu)模型，在28種缺陷下其穩(wěn)定承載力曲線見圖6，穩(wěn)定承載力與缺陷的敏感性關(guān)系見表1。在未截面優(yōu)化、等穩(wěn)定優(yōu)化、空間位置優(yōu)化和空間位置優(yōu)化[30]情況下，不同缺陷下結(jié)構(gòu)穩(wěn)定承載力標(biāo)準(zhǔn)差分別為0.32、0.31、0.40和0.46,文獻(xiàn)[30]中僅考慮了3種缺陷進(jìn)行分析，本文中考慮了28種隨機(jī)缺陷。

圖6 三種結(jié)構(gòu)模型的穩(wěn)定承載力Fig.6 Critical-limit loads of three structural models

結(jié)構(gòu)模型穩(wěn)定承載力/(kN·m-2)均值最大值最小值耗鋼量/kg 未截面優(yōu)化(33#)4.465.003.7663 193.9 空間位置分組優(yōu)化2.893.902.1851 063.9 等穩(wěn)定優(yōu)化4.725.284.1161 093.3

由圖6可見,本算例中穩(wěn)定承載力下限值為2Q=3.8 kN/m2。用空間位置分組優(yōu)化方法得到的截面優(yōu)化結(jié)果僅在16#缺陷(設(shè)計(jì)缺陷)下滿足穩(wěn)定承載力設(shè)計(jì)要求,在其他27種缺陷下均不能滿足穩(wěn)定性設(shè)計(jì)要求。未經(jīng)截面優(yōu)化的結(jié)構(gòu)模型,當(dāng)截面取值為33#截面時(shí),16#缺陷下其穩(wěn)定承載力不能滿足設(shè)計(jì)要求,9#缺陷出現(xiàn)時(shí)其承載力也是臨界滿足?；诘确€(wěn)定要求優(yōu)化得到的設(shè)計(jì)模型,在所有缺陷下均滿足穩(wěn)定承載力設(shè)計(jì)要求,且有較合理的安全儲(chǔ)備。

相比于未經(jīng)截面優(yōu)化設(shè)計(jì)的網(wǎng)殼結(jié)構(gòu),采用基于等穩(wěn)定要求設(shè)計(jì)出的結(jié)構(gòu)具有均值更大和均方差更小的穩(wěn)定承載力。采用基于空間位置分組優(yōu)化得到的結(jié)構(gòu)均值小且標(biāo)準(zhǔn)差大,即優(yōu)化后結(jié)構(gòu)的穩(wěn)定承載力對(duì)缺陷更加敏感,不能保證結(jié)構(gòu)在不同初始缺陷下的穩(wěn)定承載能力均滿足設(shè)計(jì)要求。

5 結(jié) 論

(1)將彈塑性整體穩(wěn)定承載力和隨機(jī)初始幾何缺陷引入到網(wǎng)殼結(jié)構(gòu)的截面優(yōu)化設(shè)計(jì)過程中,保證優(yōu)化設(shè)計(jì)后的網(wǎng)殼結(jié)構(gòu)具有足夠安全的穩(wěn)定承載力。

(2)在同等耗鋼量下,將較大的截面分配給缺陷敏感性大的區(qū)域,較小截面尺寸分配給缺陷敏感性小的區(qū)域,確保了網(wǎng)殼結(jié)構(gòu)在隨機(jī)出現(xiàn)的缺陷下具有魯棒性較強(qiáng)的穩(wěn)定承載力。

(3)基于等穩(wěn)定要求的截面優(yōu)化設(shè)計(jì)方法可以使結(jié)構(gòu)具有均值更大、均方差更小的穩(wěn)定承載力。通常采用的基于位置分組的截面優(yōu)化設(shè)計(jì)方法不能保證結(jié)構(gòu)在不同初始缺陷下的穩(wěn)定承載能力均滿足設(shè)計(jì)要求,是不推薦采用的優(yōu)化分組方法。

:

[1] SAKA M P, KAMESHKI E S. Optimum design of nonlinear elastic framed domes[J]. Advances in Engineering Software, 1998,29(7/8/9):519-528.

[2] TOKLU Y C, TEMUR R, BEKDAS G. Analysis of trusses by total potential optimization method coupled with harmony search[J]. Structural Engineering & Mechanics,2013,45(2):183-199.

[3] ABBASNIA R, SHAYANFAR M, KHODAM A. Reliability-based design optimization of structural systems using a hybrid genetic algorithm[J]. Structural Engineering & Mechanics, 2014,52(6):1099-1120.

[4] 中華人民共和國(guó)住房和城鄉(xiāng)建設(shè)部. 空間網(wǎng)格結(jié)構(gòu)技術(shù)規(guī)程:GJ7—2010[S]. 北京:中國(guó)建筑工業(yè)出版社, 2010.

[5] 沈世釗,陳昕.網(wǎng)殼結(jié)構(gòu)穩(wěn)定性[M].北京:科學(xué)出版社,1999:4-6.

[6] BATHE K, BOLOURCHI S. Large displacement analysis of three dimensional beam structures[J]. International Journal for Numerical Methods in Engineering, 2010,14(7):961-986.

[7] NERUBAILO B V, ZUBKOV G D, MOCHALOV M V. Stability of spirally stiffened shells under external pressure[J]. Journal of Engineering Physics and Thermophysics, 2006,9(1):202-204.

[8] 關(guān)富玲,高博青.單層網(wǎng)殼的穩(wěn)定性分析[J].工程力學(xué),1996,13(3):93-104.

GUAN Fuling, GAO Boqing. Stability analysis of single layer lattice domes[J].Engineering Mechanics, 1996,13(3):93-104.

[9] KHOT N S, VENKAYYA V B, BERKE L. Optimum structural design with stability constraints[J].International Journal for Numerical Methods in Engineering, 2010,10(5):1097-1114.

[10] SEDAGHATI R, TABARROK B. Optimum design of truss structures undergoing large deflections subject to a system stability constraint[J]. International Journal for Numerical Methods in Engineering, 2015,48(3):421-434.

[11] CEDERBAUM G, ARBOCZ J. On the reliability of imperfection-sensitive long isotropic cylindrical shells[J]. Structural Safety, 1996,18(1):1-9.

[12] CEDERBAUM G, ARBOCZ J. Reliability of shells via Koiter formulas[J]. Thin-walled Structures, 1996,24(2):173-187.

[13] AMAZIGO J C. Buckling under axial compression of long cylindrical shells with random axisymmetric imperfections[J]. Quarterly of Applied Mathematics, 1969,26(4):537.

[14] HUTCHINSON J W, MUGGERIDGE D B, TENNYSON R C. Effect of a local axisymmetric imperfection on the buckling behavior of a circular cylindrical shell under axial compression[J]. Aiaa Journal, 2012,9(9):48-52.

[15] KIYOHIRO I, KAZUO M. Critical initial imperfection of structures[J]. International Journal of Solids & Structures, 1990,26(8):865-886.

[16] 唐敢,尹凌峰,馬軍.單層網(wǎng)殼結(jié)構(gòu)穩(wěn)定性分析的改進(jìn)隨機(jī)缺陷法[J].空間結(jié)構(gòu),2004,10(4):44-47.

TANG Gan, YIN Lingfeng, MA Jun. Advanced stochastic imperfection method for stability analysis of single layer lattice domes[J].Spatial Structures, 2004,10(4):44-47.

[17] 徐長(zhǎng)航,陳國(guó)明,謝靜,等.基于支持向量機(jī)和蒙特卡洛的結(jié)構(gòu)可靠性分析方法及應(yīng)用[J].中國(guó)石油大學(xué)學(xué)報(bào)(自然科學(xué)版),2008,32(4):103-108.

XU Changhang, CHEN Guoming, XIE Jing, et al. Structural reliability analysis method based on support vector machines and Monte Carlo and its application[J].Journal of China University of Petroleum (Edition of Natural Science), 2008,32(4):103-108.

[18] 魏德敏,涂家明.單層網(wǎng)殼結(jié)構(gòu)非線性穩(wěn)定的隨機(jī)缺陷模態(tài)法研究[J].華南理工大學(xué)學(xué)報(bào)(自然科學(xué)版), 2016,44(7):83-89.

WEI Demin, TU Jiaming. A probe into nonlinear stability of single-layer reticulated shells by means of random imperfection modal method[J]. Journal of South China University of Technology(Natural Science Edition), 2016,44(7):83-89.

[19] HANSEN J S. General random imperfections in the buckling of axially loaded cylindrical shells[J]. Aiaa Journal, 2015,15(9):1250-1256.

[20] LU X Y, ZHAO X W, HUANG L L. Shape optimizing design of K6 spherical reticulated shell[J]. Advanced Materials Research, 2012,424:324-329.

[21] 王法武,唐敢.考慮整體穩(wěn)定的單層網(wǎng)殼截面優(yōu)化設(shè)計(jì)[J].空間結(jié)構(gòu),2006,12(3):31-34.

WANG Fawu, TANG Gan. Sectional optimum design of single-layer lattice shells considering structural stability[J]. Spatial Structures,2006,12(3):31-34.

[22] 陳世英,張素娟,李青.初始幾何缺陷對(duì)網(wǎng)殼截面優(yōu)化結(jié)果影響研究[J].力學(xué)與實(shí)踐,2015,37(6):708-712.

CHEN Shiying, ZHANG Sujuan, LI Qing. The influence of initial geometric imperfection on the section optimization result of reticulated shell [J]. Mechanics in Engineering, 2015,37(6):708-712.

[23] 中華人民共和國(guó)建設(shè)部.鋼結(jié)構(gòu)設(shè)計(jì)規(guī)范:GB 50017-2003[S].北京:中國(guó)計(jì)劃出版社,2003.

[24] 孫煥純.離散變量結(jié)構(gòu)優(yōu)化設(shè)計(jì)[M].大連:大連理工大學(xué)出版社,2002:79-85.

[25] 王新敏.ANSYS工程結(jié)構(gòu)數(shù)值分析[M]. 北京:人民交通出版社,2007:435-437.

[26] 尹德鈺,劉善維,錢若軍.網(wǎng)殼結(jié)構(gòu)設(shè)計(jì)[M]. 北京:中國(guó)建筑工業(yè)出版社, 1996:276-278.

[27] SCHENK C A, SCHUELLER G I. Buckling analysis of cylindrical shells with random geometric imperfections[J]. International Journal of Non-Linear Mechanics, 2003,38:1119-1132.

[28] OYESANYA M O. Influence of extra terms on asymptotic analysis of imperfection sensitivity of toroidal shell segment with random imperfection[J]. Mechanics Research Communications,2005,32(4):444-453.

[29] ROUDSARI M T, GORDINI M. Random imperfection effect on reliability of space structures with different supports[J]. Structural Engineering & Mechanics, 2015,55(3):461-472.

[30] 陳世英,李維國(guó),李青.截面優(yōu)化后網(wǎng)殼結(jié)構(gòu)幾何缺陷敏感性研究[J].空間結(jié)構(gòu),2016,22(2):16-21.

CHEN Shiying, LI Weiguo, LI Qing. Sensitivity study of geometrical imperfections of optimized reticulated shell structures[J]. Spatial Structures,2016, 22(2):16-21.