王凱 顏慶智 王中輝 劉作勛

摘要：

石化結(jié)構(gòu)經(jīng)常因其設(shè)備管線復(fù)雜的布置而造成結(jié)構(gòu)的剛度分布不均勻，設(shè)計經(jīng)常會過于保守。為了對基于位移的鋼框架結(jié)構(gòu)可靠度進(jìn)行研究，利用有限元軟件對煉化平臺鋼框架結(jié)構(gòu)進(jìn)行建模、計算、分析，在PDS（概率設(shè)計系統(tǒng)）模塊中運用Monte Carlo法對結(jié)構(gòu)位移的可靠度進(jìn)行分析，獲取煉化平臺鋼框架結(jié)構(gòu)的最大位移與隨機(jī)輸入變量之間關(guān)系的計算結(jié)果。該結(jié)果表明該煉化平臺鋼框架結(jié)構(gòu)的設(shè)計比較保守，造成了鋼材的浪費。為了實現(xiàn)結(jié)構(gòu)設(shè)計的經(jīng)濟(jì)性，應(yīng)根據(jù)計算得到的靈敏度圖，在滿足結(jié)構(gòu)可靠性的基礎(chǔ)上，將圖中非敏感性因子的截面尺寸進(jìn)行優(yōu)化再設(shè)計，得出安全經(jīng)濟(jì)的設(shè)計方案，為以后石油化工行業(yè)鋼框架結(jié)構(gòu)的設(shè)計提供參考。

關(guān)鍵詞：

鋼框架；結(jié)構(gòu)位移；Monte Carlo；可靠度

中圖分類號：TU 391文獻(xiàn)標(biāo)識碼：A文章編號：1001-005X（2018）01-0091-06

Abstract：

Petrochemical structures are often characterized by uneven distribution of stiffness due to the complicated arrangement of its equipment pipelines，and designs are often overly conservative.In order to analysis of the steel frame structure based on displacement，the finite element software was used to model，calculate and analyze the steel frame structure of refining platform.In the PDS（the probability design system）module，the Monte Carlo method was used to analyze the reliability of structural displacement in the PDS module to obtain the calculation results of the relationship between the maximum displacement and the random input variables of the steel frame structure of the refining platform.The results revealed that the section size of the steel frame structure was more conservative，resulting in the waste of steel.In order to realize the design of the structures economy，we should optimize and redesign the crosssectional dimensions of nonsensitive factors in the graph based on the calculated sensitivity diagram to meet the structural reliability，and draw a safe and economical design plan for the future petrochemical industry，which could provide reference for the future design of steel frame structure in petrochemical industry.

Keywords：

steel frame structure；displacement；Monte Carlo；reliability

0引言

鋼結(jié)構(gòu)在石油化工領(lǐng)域中應(yīng)用非常廣泛，石化產(chǎn)業(yè)的煉化平臺鋼框架結(jié)構(gòu)因其生產(chǎn)設(shè)備與生產(chǎn)工藝上的要求，而具有自己獨特的特點。石化設(shè)備大多是露天或敞開式布置，由于生產(chǎn)工藝的要求，管線的布置往往是縱橫交錯的，數(shù)量多而且質(zhì)量以及尺寸都很大，因而造成結(jié)構(gòu)的剛度和質(zhì)量分布不均；大多數(shù)的石化設(shè)備以及管線的布置是不均勻的，造成了結(jié)構(gòu)受力的復(fù)雜，設(shè)計方法也同鋼結(jié)構(gòu)設(shè)計規(guī)范有所差異，進(jìn)行結(jié)構(gòu)設(shè)計時，設(shè)計經(jīng)常會過于保守，而且對于鋼框架平臺的研究較少，因此，需要對煉化平臺鋼框架結(jié)構(gòu)進(jìn)行優(yōu)化分析[1]。

1蒙托卡羅法

結(jié)構(gòu)體系可靠度的計算方法大致包括：失效模式法、Monte Carlo法、響應(yīng)面法和隨機(jī)有限元法等。Monte Carlo法[2-3]不受功能函數(shù)非線性的影響，可以通過多次模擬試驗就可以得到比較精確的結(jié)果。Monte Carlo法又稱為概率統(tǒng)計模擬方法，該方法是按隨機(jī)抽樣規(guī)律方法產(chǎn)生子樣本，而后進(jìn)行子樣本的統(tǒng)計分析，分析解出問題的值。在ANSYS可靠度數(shù)值分析模塊中Monte Carlo法的抽樣方法包括直接抽樣法與拉丁超立方抽樣方法，也可以自定義抽樣類型，拉丁超立方抽樣方法[4]其基本思路是把隨機(jī)變量的采樣區(qū)域通過細(xì)分為等概率空間將樣本均勻分布在整個問題域中，并在每個段中繪制一個樣本。拉丁超立方抽樣方法[5-6]可以實現(xiàn)樣本數(shù)量的減少，而且拉丁超立方體抽樣方法對于估計隨機(jī)結(jié)構(gòu)中的平均值與方差非常有效，而且抽樣效率比單純的Monte Carlo抽樣效率略高。本文采用的是Monte Carlo法的拉丁超立方抽樣法對煉化平臺空間鋼框架進(jìn)行可靠度分析。

2基于ANSYS及Monte Carlo數(shù)值模擬

2.1工程實例

本工程項目名稱為1.60 Mt/a DCC生產(chǎn)烯烴項目，建設(shè)地點位于山東墾利。建筑工程等級為一級；耐火等級為二級。結(jié)構(gòu)類型為多層鋼框架結(jié)構(gòu)；建筑層數(shù)為地上三層；建筑總高度為22.495 m。設(shè)計使用年限為50 a?？拐痤悇e為丙類，抗震設(shè)防烈度為 7度，設(shè)計基本地震加速度值為 0.10 g，鋼框架結(jié)構(gòu)在多遇地震荷載作用下的阻尼比采用0.04。材料采用Q235鋼，這里給出結(jié)構(gòu)二層平面布置，其它結(jié)構(gòu)層除設(shè)備梁的截面不同外布置均相同，結(jié)構(gòu)平面布置如圖1所示。

該建筑結(jié)構(gòu)形式為輕型鋼框架結(jié)構(gòu)，長度63 m，其柱距為7 m，跨度4.35 m，該石油煉化平臺為多層框架結(jié)構(gòu)，其中樓面材料采用5 mm厚花紋鋼板。結(jié)構(gòu)共兩跨，柱與柱之間采用支撐連接，設(shè)備梁與框架梁間采用簡支連接。在計算地震荷載作用于結(jié)構(gòu)時，其正常運行時重力荷載代表值應(yīng)取永久荷載與可變重力荷載組合值之和，平臺活荷載為安裝檢修與正常生產(chǎn)之和取14 kN/m2，活荷載組合系數(shù)采用0.5，永久荷載的組合系數(shù)采用1.0[7]。

材料數(shù)據(jù)：

彈性模量E=2.1×1011Pa，泊松比μ=0.3，質(zhì)量密度為7 850 kg/m3。

2.3ANSYS建立模型的分析步驟：

對工程實例運用ANSYS對煉化平臺鋼框架結(jié)構(gòu)進(jìn)行建模計算[9-11]。與其它鋼框架結(jié)構(gòu)不同，直接或者間接支承石化設(shè)備的設(shè)備梁是結(jié)構(gòu)抗震的第一道防線，設(shè)備梁的抗震性能直接影響抗震設(shè)計的成敗，因此，簡化模型計算時，省略對平臺板加勁肋的建模，同時，選取了鋼框架結(jié)構(gòu)中具有代表性的三榀框架進(jìn)行建模計算，鋼框架結(jié)構(gòu)建模選用的單元類型如下：

梁柱單元采用Beam188單元；板單元選用Shell181單元；支撐桿采用Link180單元。

鋼框架結(jié)構(gòu)建模選用的單元的截面形式見表1。

鋼框架結(jié)構(gòu)建模選用單元的邊界條件為：柱底邊界條件為全約束，設(shè)備梁約束為X、Y、Z方向自由度，建立的模型如圖2所示。

ANSYS分析步驟：首先建立幾何模型、選擇單元類型、設(shè)置材料數(shù)據(jù)、對幾何模型賦予屬性。施加位移邊界條件，施加模型荷載，進(jìn)行計算，進(jìn)入ANSYS的PDS可靠度分析模塊對模型進(jìn)行模擬分析計算，模擬步驟如圖3所示。

模型計算結(jié)果如圖4所示。

在石油化工鋼結(jié)構(gòu)冷換框架設(shè)計規(guī)范中結(jié)構(gòu)在地震荷載作用下的結(jié)構(gòu)變形驗算為層位移與層高之比不宜大于1/250[12]。地震作用下結(jié)構(gòu)的位移限值為：

7 500×1/250=30 mm

由圖4結(jié)構(gòu)位移云圖可以看出結(jié)構(gòu)的最大位移在第二層，最大位移約為5.85 mm，由表 2 ANSYS模型計算結(jié)果得置信區(qū)間為[0.984 556，0.999 540]，置信水平為95%條件下的DXMAX小于8×10-3的概率值為99.599%，該位移遠(yuǎn)小于30 mm。

2.4結(jié)構(gòu)可靠度分析

結(jié)構(gòu)的可靠性分析一般是通過概率模型來估算失效概率以此來估算可靠度的，而可靠性分析方法的主要目標(biāo)就是近似的估算失效概率Pf[13]，失效概率估計：

其中，N是模擬總數(shù)，n表示系統(tǒng)失敗的模擬次數(shù)。采用Monte Carlo法對于復(fù)雜的工程問題計算是比較容易的，因此，概率設(shè)計過程中采用蒙特卡羅模擬方法進(jìn)行概率設(shè)計。在利用ANSYS進(jìn)行數(shù)值模擬時，可靠度分析部分定義的輸出變量分別為DXMAX、DZMAX（X、Z兩個方向的最大位移），對數(shù)據(jù)的抽樣方法采取拉丁超立方抽樣方法進(jìn)行仿真計算，設(shè)置仿真循環(huán)次數(shù)為300次，假定地震荷載、鋼結(jié)構(gòu)的構(gòu)件截面面積分布服從對數(shù)正態(tài)分布，彈性模量、材料密度、泊松比服從高斯正態(tài)分布[15-17]。具體輸入變量的參數(shù)見表3。

2.5數(shù)值模擬結(jié)果分析

由圖5、6隨機(jī)變量輸入的樣本均值分布圖、樣本標(biāo)準(zhǔn)差分布圖可知樣本的抽樣次數(shù)足夠，均值、標(biāo)準(zhǔn)差均收斂，可以滿足結(jié)構(gòu)可靠度計算的水平。

由圖7、8位移樣本值的直方圖可以看出在地震荷載作用下，空間框架結(jié)構(gòu)的X、Z向最大位移樣本值的概率分布服從極值Ⅰ型分布。

由圖9可以看出該煉化平臺框架在Z向的最大位為0.21 328 mm，Z向位移非常小，該石油煉化平臺鋼框架結(jié)構(gòu)為框架-支撐體系，柱間支撐限制了結(jié)構(gòu)Z向的位移，結(jié)構(gòu)的最大位移出現(xiàn)在結(jié)構(gòu)的X方向，只需分析結(jié)構(gòu)X向位移即可。

在結(jié)構(gòu)優(yōu)化方面，經(jīng)常會權(quán)衡結(jié)構(gòu)的經(jīng)濟(jì)性、適用性，以及其它影響結(jié)構(gòu)性能的因素，因此，在進(jìn)行可靠度分析時，研究各個因素對結(jié)構(gòu)失效的影響大小對于結(jié)構(gòu)的優(yōu)化具有重要意義[18]。在輸入的隨機(jī)變量中，確定性輸入變量和參數(shù)變量都稱為輸入?yún)?shù)。靈敏度是分析結(jié)構(gòu)可靠性的度量，這些輸入?yún)?shù)的變化對于簡單地評估結(jié)構(gòu)的可靠性變化以及對已定的結(jié)構(gòu)設(shè)計變更具有很重要的作用。靈敏度分析就是結(jié)合優(yōu)化程序使用迭代方法，對結(jié)構(gòu)進(jìn)行優(yōu)化再設(shè)計分析以使設(shè)計的總成本最小[19]。

圖10、11柱狀圖是根據(jù)靈敏度大小從左到右依次排列，可以看出作用在結(jié)構(gòu)二層的外荷載的靈敏度最大，對結(jié)構(gòu)X向位移影響最大，彈性模量的靈敏度小于零表明它對結(jié)構(gòu)X向位移起到負(fù)作用。在ANSYS的后處理模塊中可以對已定的輸入變量進(jìn)行靈敏度分析，然后判斷哪些參數(shù)可以實現(xiàn)結(jié)構(gòu)的優(yōu)化，以達(dá)到結(jié)構(gòu)優(yōu)化的目的[20]。該煉化平臺鋼框架結(jié)構(gòu)的可靠度是可以接受的，可以對結(jié)構(gòu)非敏感性因子進(jìn)行優(yōu)化，即對截面1-8進(jìn)行優(yōu)化設(shè)計，在保證結(jié)構(gòu)安全可靠的情況下，適當(dāng)?shù)臏p少構(gòu)件的截面面積[21]。以此來達(dá)到提高煉化平臺鋼框架結(jié)構(gòu)設(shè)計經(jīng)濟(jì)性的目的。

3結(jié)論

本文利用有限元軟件對煉化平臺鋼框架結(jié)構(gòu)進(jìn)行建模，對地震荷載作用下結(jié)構(gòu)的位移進(jìn)行可靠度分析。煉化平臺鋼框架結(jié)構(gòu)在地震荷載的作用下，水平最大位移為5.9 mm，遠(yuǎn)小于規(guī)范設(shè)計值30 mm，可以看出該框架結(jié)構(gòu)的設(shè)計比較保守，鋼材比較浪費。由靈敏度圖可知，需要將該煉化平臺鋼框架結(jié)構(gòu)的各梁柱敏感性因子低的截面尺寸進(jìn)行優(yōu)化再設(shè)計。該煉化平臺鋼框結(jié)構(gòu)可為以后尋求結(jié)構(gòu)可靠性與經(jīng)濟(jì)性的設(shè)計方法提供借鑒。

由于本文對于鋼框架結(jié)構(gòu)體系并沒有進(jìn)行優(yōu)化再設(shè)計分析，因此，需要對截面優(yōu)化再設(shè)計后的鋼框結(jié)構(gòu)做進(jìn)一步的可靠度分析和研究。

【參考文獻(xiàn)】

[1]高劍，王忠凱，潘毅，等.罕遇地震下石化鋼結(jié)構(gòu)減震的關(guān)鍵影響因素[J].土木建筑與環(huán)境工程，2016，38（1）：93-94.

[2]蔡雁蕊.基于靈敏度分析的鋼框架支撐結(jié)構(gòu)體系抗側(cè)性能研究[D].廣州：廣東工業(yè)大學(xué)，2016.

[3]戚毅婷.基于隨機(jī)有限元的框架結(jié)構(gòu)可靠性分析[D].昆明：昆明理工大學(xué)，2012.

[4]Baalousha H M.Approximation of the exponential integral（well function）using sampling methods[J].Journal of Hydrology，2015，523：278-280.

[5]Olsson A，Sandberg G，Dahlblom O.On Latin hypercube sampling for structural reliability analysis[J].Structural Safety，2003，25（1）：47-68.

[6]Olsson A，Sandberg G E.Latin hypercube sampling for stochastic finite element analysis[J].Journal of Engineering Mechanics，2002：128（1）：121-125.

[7]中國石化工程建設(shè)公司.化工、石化建（構(gòu)）筑物荷載設(shè)計規(guī)范：HG/T20674-200[S].北京：化工出版社，2005.

[8]中華人民共和國和城鄉(xiāng)建設(shè)部.建筑抗震設(shè)計規(guī)范：GB50011-2010[S].北京：中國建筑工業(yè)出版社，2010.

[9]張仙平，胡仁喜，康士廷，等.ANSYS16.0土木工程有限元分析從入門到精通[M].北京：機(jī)械工業(yè)出版社，2016.

[10]徐鶴山.ANSYS建筑鋼結(jié)構(gòu)工程實例分析[M].北京：機(jī)械工業(yè)出版社，2007.

[11]博弈創(chuàng)作室.ANSYS9.0經(jīng)典產(chǎn)品基礎(chǔ)教程與實例詳解[M].北京：中國水利水電出版社，2006.

[12]中國石化工程建設(shè)有限公司.石油化工鋼結(jié)構(gòu)冷換框架設(shè)計規(guī)范：SH/T 3077—2012[S].北京：中國石化出版社，2013.

[13]Keshtegar B，Kisi O.M5 model tree and Monte Carlo simulation for efficient structural reliability analysis[J].Applied Mathematical Modelling，2017，48：899-901.

[14]Zhang H，Shayan S，Rasmussen K J R，et al.Systembased design of planar steel frames，I：Reliability framework[J].Journal of Constructional Steel Research，2016，123：135-143.

[15]張偉.結(jié)構(gòu)可靠性理論與應(yīng)用[M].北京：科學(xué)出版社，2008.

[16]趙國藩，金偉良，貢金鑫.結(jié)構(gòu)可靠度理論[M].北京：中國建筑工業(yè)出版社，2000.

[17]王湛，何嘉年.基于蒙特卡羅法的鋼框架結(jié)構(gòu)可靠度計算方法的研究[C]//全國鋼結(jié)構(gòu)學(xué)術(shù)年會論文集.中國鋼結(jié)構(gòu)協(xié)會，2008，10：55-59.

[18]李遠(yuǎn)瑛，張德生.基于ANSYS的平面框架結(jié)構(gòu)位移可靠度分析[J].湘潭大學(xué)自然科學(xué)學(xué)報，2011，33（2）：45-49.

[19]Papadrakakis M，Lagaros N D，Plevris V.Design optimization of steel structures considering uncertainties[J].Engineering Structures，2005，27（9）：1408-1418.

[20]Ditlevsen O D，Madsen H O.Structural reliability methods[M].Chichester：John Wiley & Sons，Inc.，1996.

[21]謝鐳，于英霞，張偉.基于強(qiáng)度的鋼框架可靠性分析與計算[J].煉油技術(shù)與工程，2009，39（7）：45-49.