單層球面網(wǎng)殼結(jié)構(gòu)隨機(jī)地震響應(yīng)與可靠度分析

徐軍　袁子豪

摘? ?要：為評(píng)價(jià)單層球面網(wǎng)殼結(jié)構(gòu)抗震性能，需充分考慮地震激勵(lì)的隨機(jī)性以及非平穩(wěn)性. 結(jié)合概率密度演化方法，可得到單層球面網(wǎng)殼結(jié)構(gòu)動(dòng)力響應(yīng)的概率信息及可靠度.首先，采用譜表示-隨機(jī)函數(shù)方法，生成符合不同設(shè)防要求的全非平穩(wěn)地震動(dòng)，其平均反應(yīng)譜與規(guī)范反應(yīng)譜擬合效果理想，通過概率密度演化方法考察了單層球面網(wǎng)殼結(jié)構(gòu)最大位移的概率密度演化全過程;進(jìn)一步地，分別從宏觀和微觀層次上選取2個(gè)失效評(píng)價(jià)指標(biāo)，通過引入等效極值分布思想，可獲得單層球面網(wǎng)殼結(jié)構(gòu)抗震可靠度分析所關(guān)心的動(dòng)力響應(yīng)極值分布;最后通過參數(shù)分析，評(píng)價(jià)了不同參數(shù)對(duì)單層球面網(wǎng)殼結(jié)構(gòu)抗震可靠度的影響. 與傳統(tǒng)方法相比，本文方法從隨機(jī)性的角度評(píng)價(jià)網(wǎng)殼結(jié)構(gòu)的抗震性能更加合理，并且兼顧了效率和精度.

關(guān)鍵詞：?jiǎn)螌忧蛎婢W(wǎng)殼;地震動(dòng);概率密度演化;隨機(jī)動(dòng)力響應(yīng);動(dòng)力可靠度

中圖分類號(hào)：TU393.3? ? ? ? ? ? ? ? ? ? ? ? ? ? ?文獻(xiàn)標(biāo)志碼：A

Stochastic Seismic Response and Reliability Analysis of

Single-layer Spherical Reticulated Shells

XU Jun，YUAN Zihao

（College of Civil Engineering，Hunan University，Changsha 410082，China）

Abstract： To evaluate the seismic performance of single-layer spherical reticulated shell， the randomness and non-stationarity of seismic excitations need to be sufficiently taken into account. In conjunction with the probability density evolution method （PDEM）， the probabilistic information of dynamic response and reliability can be readily obtained. First， the random function based spectral representation method is employed to model the fully non-stationary seismic ground motions， which meets the requirements of different fortifications. Moreover， the average response spectrum accords well with the target response spectrum. Then， the instantaneous probability density evolution process of the shells maximum displacement is investigated from the perspective of PDEM. Further， two indices， which characterize the failure of the shell structure， are selected from macro and micro levels， respectively. By introducing the idea of equivalent extreme-value event， one can obtain the extreme value distribution of response for seismic reliability analysis. Finally， through carrying out the parametric analyses， the effect of different parameters on the seismic reliability of the shell can be evaluated. Compared with the traditional methods， the proposed method can evaluate the seismic performance of single-layer spherical reticulated shells from the perspective of randomness， which is much more reasonable. Besides， the tradeoff of efficiency and accuracy can be also ensured.

Key words： single-layer spherical reticulated shell;ground motion;probability density evolution method;stochastic dynamic response;dynamic reliability

網(wǎng)殼結(jié)構(gòu)因造型優(yōu)美、受力性能良好而被廣泛應(yīng)用于諸多重大工程.隨著設(shè)計(jì)理論和建造技術(shù)的發(fā)展，其跨度在不斷增大，大跨度網(wǎng)殼在強(qiáng)震作用下的力學(xué)性能受到眾多學(xué)者關(guān)注[1-8]. 其中，相當(dāng)一部分研究集中于揭示單層球面網(wǎng)殼在地震荷載作用下的失效機(jī)理，并已取得較為系統(tǒng)的研究成果：針對(duì)動(dòng)力失穩(wěn)，文獻(xiàn)[3-5]從微觀到宏觀層次上提出了應(yīng)力變化率法、特征位移響應(yīng)法以及特征能量法等;針對(duì)強(qiáng)度失效，文獻(xiàn)[6]提出了以結(jié)構(gòu)的特征響應(yīng)為判定指標(biāo)（如最大節(jié)點(diǎn)位移和塑性桿比例等）;文獻(xiàn)[7]提出了基于能量和最大變形的雙控準(zhǔn)則;文獻(xiàn)[8]基于大量統(tǒng)計(jì)數(shù)據(jù)建立了損傷模型，通過損傷因子判定結(jié)構(gòu)失效與否. 綜上所述，地震作用下網(wǎng)殼結(jié)構(gòu)失效機(jī)理的研究已基本形成完整的理論框架，如何充分利用已有成果，成為今后的任務(wù)之一.

另一方面，目前針對(duì)復(fù)雜結(jié)構(gòu)的抗震分析，往往只選擇幾條實(shí)際地震記錄及人造地震波進(jìn)行確定性動(dòng)力分析.然而，地震動(dòng)具有顯著的隨機(jī)性，這使得工程結(jié)構(gòu)的響應(yīng)性態(tài)難以把控，特別當(dāng)?shù)卣痣S機(jī)性與網(wǎng)殼結(jié)構(gòu)非線性相耦合時(shí)，問題將變得更為復(fù)雜. 在這種情形下，如何準(zhǔn)確把握結(jié)構(gòu)地震響應(yīng)精細(xì)的概率信息顯得尤為重要，換言之，我們須考察單層球面網(wǎng)殼結(jié)構(gòu)在隨機(jī)地震激勵(lì)下的響應(yīng). 于1950年左右發(fā)展起來的隨機(jī)振動(dòng)方法較為充分地考慮地震發(fā)生的概率統(tǒng)計(jì)特性，然而其計(jì)算量巨大，很難應(yīng)用于工程實(shí)踐中.從計(jì)算力學(xué)角度，林家浩教授等對(duì)結(jié)構(gòu)的平穩(wěn)和非平穩(wěn)隨機(jī)響應(yīng)發(fā)展了高效精確的算法-虛擬激勵(lì)法[9].基于此方法，文獻(xiàn)[10-11]針對(duì)大跨空間結(jié)構(gòu)逐步推導(dǎo)了多維、多點(diǎn)非平穩(wěn)地震動(dòng)虛擬激勵(lì)法的計(jì)算公式. 盡管如此，虛擬激勵(lì)法本質(zhì)上并不適用于強(qiáng)非線性結(jié)構(gòu). 21世紀(jì)初，李杰和陳建兵從物理隨機(jī)系統(tǒng)的基本思想出發(fā)，提出廣義概率密度演化方程，為復(fù)雜工程結(jié)構(gòu)非線性隨機(jī)動(dòng)力響應(yīng)及可靠度分析提供了新的途徑[12].

此外，為進(jìn)行隨機(jī)地震響應(yīng)分析，如何建立隨機(jī)地震動(dòng)模型是關(guān)鍵.其中，譜表示方法由于其理論完善、易于實(shí)現(xiàn)而在工程實(shí)踐中被廣為應(yīng)用.在譜表示方法基礎(chǔ)之上，劉章軍等[13]僅用一或兩個(gè)基本隨機(jī)變量構(gòu)造正交隨機(jī)函數(shù)來模擬隨機(jī)地震動(dòng)，從而達(dá)到降維目的，大大降低計(jì)算工作量.結(jié)合概率密度演化方法，可實(shí)現(xiàn)復(fù)雜工程結(jié)構(gòu)的隨機(jī)地震響應(yīng)和可靠度分析.本文從隨機(jī)動(dòng)力響應(yīng)來考察單層球面網(wǎng)殼的抗震性能，采用譜表示-隨機(jī)函數(shù)方法模擬全非平穩(wěn)隨機(jī)地震動(dòng)，結(jié)合概率密度演化方法計(jì)算單層球面網(wǎng)殼隨機(jī)地震響應(yīng)和動(dòng)力可靠度.

1? ?非平穩(wěn)隨機(jī)地震動(dòng)模型

1.1? ?譜表示-隨機(jī)函數(shù)方法

（1）

這里{Xk，Yk}為一組標(biāo)準(zhǔn)正交隨機(jī)變量，需滿足：

基于隨機(jī)函數(shù)的思想，假設(shè)一組標(biāo)準(zhǔn)正交隨機(jī)變量{Xn，Yn}[13]：

（3）

式中：cas（·）=cos（·）+sin（·），為Hartley正交基函數(shù);？專1和？專2表示區(qū)間[0，2π]上相互獨(dú)立的基本隨機(jī)變量. {Xn，Yn}經(jīng)過某種特殊的隨機(jī)映射方式可得到{Xk，Yk}. 通過這一處理，隨機(jī)過程的隨機(jī)度從∞→2N→2，從而極大地降低問題的分析難度.

式（1）中，Sa（t，ω）為演變功率譜函數(shù)，可表示為：

其中S（ω）為平穩(wěn)過程功率譜密度函數(shù)，本文采用Clough和Penzien提出的雙過濾白噪聲模型[14]：

（5）

式中：ωg和ξg分別表示場(chǎng)地土的卓越圓頻率和阻尼比，ωg = 2π/Tg，Tg和ξg可通過查詢《建筑抗震設(shè)計(jì)規(guī)范》（GB50011—2010）[15]得到;ωf和ξf分別為第2過濾層的卓越圓頻率和阻尼比，文獻(xiàn)[16]建議取ωf = 0.1ωg，ξf = ξg;S0為譜強(qiáng)度因子，可表示為：

這里amax為地震動(dòng)峰值加速度平均值，取200 cm/s2;γ為峰值因子，為符合規(guī)范要求，取值參考文獻(xiàn)[16]. 本文采用文獻(xiàn)[16]建議的時(shí)-頻調(diào)制函數(shù)：

式中：

這里b = a + 0.001，c = 0.005，a可參考文獻(xiàn)[16].

1.2? ?地震動(dòng)樣本及反應(yīng)譜擬合

考慮抗震設(shè)防烈度為8度，場(chǎng)地類別為Ⅲ類，設(shè)計(jì)地震分組為第2組，結(jié)構(gòu)阻尼比為0.02，其它參數(shù)按文獻(xiàn)[16]選取. 主要步驟如下：

①按照華羅庚-王元的數(shù)論方法[17]在[0.2π）×[0.2π）對(duì)基本隨機(jī)變量？專1和？專2取點(diǎn)，本文取點(diǎn)總數(shù)為610，并基于Voronoi域[12]對(duì)樣本點(diǎn)進(jìn)行概率剖分，計(jì)算賦得概率;

②按式（3）生成標(biāo)準(zhǔn)正交函數(shù){Xn，Yn}樣本，然后利用Matlab的rand（‘state，0）和randperm（N）映射得到{Xk，Yk}樣本;

③應(yīng)用式（1）得到地震動(dòng)時(shí)程代表性樣本.

取截?cái)囝l率ωu = 240 rad/s，N = 1 601，Δω = 0.15 rad/s，地震動(dòng)持時(shí)為30 s，時(shí)間間隔為0.01 s. 圖1為演變功率譜密度圖. 由圖1可知：能量主要集中于50 rad/s以內(nèi)，峰值大致出現(xiàn)在6～7 s左右. 圖2為兩條地震動(dòng)樣本，樣本在時(shí)域和頻域上具有典型的非平穩(wěn)特征.

圖3給出本文方法與經(jīng)典譜表示方法（610條隨機(jī)樣本）所計(jì)算得到的時(shí)程樣本均值、標(biāo)準(zhǔn)差與目標(biāo)值的對(duì)比圖. 經(jīng)計(jì)算，采用本文方法和經(jīng)典譜表示方法所得的均值誤差的最大范數(shù)分別為0.003 6和9.224 8，標(biāo)準(zhǔn)差相對(duì)誤差二范數(shù)分別為1.68%和2.98%.可見：與經(jīng)典譜表示方法相比，本文方法得到的均值及標(biāo)準(zhǔn)差與目標(biāo)值擬合更好.另一方面，目前各類結(jié)構(gòu)抗震設(shè)計(jì)規(guī)范均以反應(yīng)譜作為地震動(dòng)輸入的依據(jù)，通常平均反應(yīng)譜與規(guī)范反應(yīng)譜有一定的擬合精度要求. 圖4（a）為610條地震動(dòng)代表性樣本平均反應(yīng)譜與規(guī)范反應(yīng)譜的對(duì)比圖，可見二者之間差異較大，在長(zhǎng)周期部分差異更為明顯. 對(duì)于高層建筑或大跨結(jié)構(gòu)的抗震設(shè)計(jì)，長(zhǎng)周期地震動(dòng)分量起著十分重要的作用，因此，長(zhǎng)周期部分反應(yīng)譜的擬合需有較高精度[18]. 文獻(xiàn)[19]提出一種簡(jiǎn)單的修正方法，效果顯著. 圖4（b）為采用本文方法修正后平均反應(yīng)譜與目標(biāo)反應(yīng)譜的對(duì)比圖，可知：經(jīng)過修正得到的平均反應(yīng)譜與目標(biāo)反應(yīng)譜擬合效果十分理想，其中平均誤差和最大誤差分別為1.7%和8.4%.

2? ?單層球面網(wǎng)殼隨機(jī)地震響應(yīng)分析

2.1? ?結(jié)構(gòu)模型描述

本文采用K8型單層球面網(wǎng)殼為計(jì)算模型，其平面圖和立面圖如圖5所示. 相關(guān)參數(shù)設(shè)置如下：跨度L = 60 m，矢跨比f/L = 1/3，底部支承處為三向固定鉸接，上部桿件間均為剛接;所有桿件截面尺寸設(shè)置為？椎114 mm × 3 mm;材料選彈塑性鋼材，其屈服強(qiáng)度為235 MPa，彈性模量取E = 2.06 × 105 MPa，切線模量Et = 0.02E，考慮材料非線性，采用雙線性隨動(dòng)強(qiáng)化本構(gòu)模型（BKIN）;網(wǎng)殼結(jié)構(gòu)表面作用有均布面荷載，其大小為P = 1.8 kN/m2，計(jì)算時(shí)按靜力等效原則將殼體表面荷載轉(zhuǎn)換為節(jié)點(diǎn)集中質(zhì)量凝聚在各質(zhì)量節(jié)點(diǎn)處，用點(diǎn)單元MASS21模擬. 此外，本文采用ANSYS軟件計(jì)算，選用PIPE20單元，桿件分為3段. 計(jì)算過程中采用Rayleigh阻尼，阻尼比ξ = 0.02. 為節(jié)約計(jì)算成本，地震動(dòng)持時(shí)一般截取地震動(dòng)時(shí)程中較大的前20 s部分，時(shí)間增量步長(zhǎng)取0.02 s，施加時(shí)為三向輸入，峰值加速度之比amaxX ∶ amaxY ∶ amaxZ = 1 ∶ 0.85 ∶ 0.65.

2.2? ?概率密度演化方法

在隨機(jī)激勵(lì)作用下，多自由度體系的運(yùn)動(dòng)方程可表示為：

式中：M和C分別為質(zhì)量和阻尼矩陣;f（X）為恢復(fù)力向量;，和X分別為相對(duì)加速度、速度和位移響應(yīng)向量;I為單位向量，g（？專，t）代表上述非平穩(wěn)隨機(jī)地震動(dòng)激勵(lì).地震動(dòng)的隨機(jī)性源于基本隨機(jī)變量？專 = {？專1，？專2}.

一般來講，結(jié)構(gòu)系統(tǒng)中的任意物理量，如位移、內(nèi)力、變形等都是存在且唯一依賴于？專的. 不妨設(shè)感興趣的物理量為：

相應(yīng)地，該物理量的速度形式可表示為：

基于概率守恒原理，進(jìn)一步推導(dǎo)可得[12]：

其中，pZ？專（z，θ，t）為t時(shí)刻增廣系統(tǒng)（z，？專）的聯(lián)合概率密度函數(shù). 式（12）的初始條件可表示為：

式中：δ為狄拉克函數(shù);z0為確定性的初始值;p？專（θ）為輸入隨機(jī)變量聯(lián)合概率密度函數(shù).結(jié)合運(yùn)動(dòng)方程（9），數(shù)值求解偏微分方程（12）可得到Z（t）的概率密度演化過程[12]，即：

其中，數(shù)值方法采用TVD格式和單邊差分格式的線性組合形式，組合系數(shù)為0.95.

基于上述方法，本文對(duì)單層球面網(wǎng)殼進(jìn)行610次確定性動(dòng)力時(shí)程分析，通過考察網(wǎng)殼結(jié)構(gòu)的特征響應(yīng)，如最大節(jié)點(diǎn)位移，來反映結(jié)構(gòu)在隨機(jī)地震作用下力學(xué)性能的概率信息.同時(shí)，基于經(jīng)典譜表示方法，按照蒙特卡洛模擬方法（MCS）生成10 000條地震動(dòng)，輸入到結(jié)構(gòu)進(jìn)行動(dòng)力時(shí)程分析，將計(jì)算得到的網(wǎng)殼結(jié)構(gòu)最大節(jié)點(diǎn)位移均值和標(biāo)準(zhǔn)差與概率密度演化方法（PDEM）計(jì)算結(jié)果進(jìn)行對(duì)比，如圖6所示. 可見：在4 s以前MCS結(jié)果與PDEM結(jié)果擬合較好，

4 s以后二者出現(xiàn)一定偏差，不過PDEM結(jié)果始終在MCS結(jié)果周圍波動(dòng)，且誤差在可接受范圍之內(nèi).

圖7給出單層球面網(wǎng)殼結(jié)構(gòu)在[8，11]s范圍內(nèi)最大節(jié)點(diǎn)位移的概率密度演化過程. 由圖可知，網(wǎng)殼節(jié)點(diǎn)最大位移響應(yīng)的概率分布完全不同于通常假定的正態(tài)分布，具有非規(guī)則性，演化曲面就像連綿不斷的山峰;等概率密度線就像流動(dòng)的水流，反映了概率在狀態(tài)空間內(nèi)的演變.

3? ?單層球面網(wǎng)殼動(dòng)力可靠度分析

基于概率密度演化方法，可以得到結(jié)構(gòu)特征響應(yīng)的概率密度演變信息.然而，為精細(xì)化評(píng)價(jià)結(jié)構(gòu)的可靠性，需給出結(jié)構(gòu)可靠性的概率度量，亦即給出可靠度分析.本文采用文獻(xiàn)[20]所提出的等價(jià)極值事件的動(dòng)力可靠度分析方法，結(jié)合相應(yīng)的強(qiáng)度失效評(píng)價(jià)指標(biāo)，計(jì)算單層球面網(wǎng)殼結(jié)構(gòu)的動(dòng)力可靠度.

3.1? ?基于極值分布的動(dòng)力可靠度分析

對(duì)于球面網(wǎng)殼其中某個(gè)節(jié)點(diǎn)/單元，其動(dòng)力可靠度可表示為：

中：Xi表示結(jié)構(gòu)特征響應(yīng);b為結(jié)構(gòu)失效時(shí)對(duì)應(yīng)的閾值;T代表地震動(dòng)持時(shí);i表示網(wǎng)殼結(jié)構(gòu)的節(jié)點(diǎn)/單元號(hào);m表示節(jié)點(diǎn)/單元數(shù). 對(duì)于首次超越事件的動(dòng)力可靠度分析，式（15）也可表達(dá)為：

殼結(jié)構(gòu)的整體可靠度可表示為：

max（T） = ?{{Xi（t）}}，則有：

中p max（x，T）是max（T）的概率密度函數(shù)，亦稱為等價(jià)極值分布.一旦確定網(wǎng)殼結(jié)構(gòu)響應(yīng)的等價(jià)極值分布，即可方便地計(jì)算出網(wǎng)殼結(jié)構(gòu)的整體可靠度.

max（T）為僅依賴于基本隨機(jī)向量？專和時(shí)間段[0，T]的隨機(jī)變量，進(jìn)一步有：

構(gòu)造一個(gè)時(shí)間虛擬過程：

需滿足下述條件：

本文選取的虛擬過程形式為[20]：

其中，ω = 2.5π，τ = 1.

類似地，有廣義概率密度演化方程

其初始條件為：

則有[12]：

那么，等價(jià)極值分布則可表示為：

3.2? ?數(shù)值算例

網(wǎng)殼結(jié)構(gòu)強(qiáng)度失效問題一直以來都是人們研究的熱點(diǎn)，正確選取評(píng)價(jià)指標(biāo)尤為重要. 本文按照文獻(xiàn)[6]的建議，選取最大節(jié)點(diǎn)位移和塑性桿比例分別從局部和整體、宏觀和微觀的角度來刻畫網(wǎng)殼結(jié)構(gòu)的力學(xué)性能.同時(shí)，通過選取不同參數(shù)，分別從結(jié)構(gòu)本身、地震動(dòng)兩個(gè)角度來研究不同參數(shù)對(duì)單層球面網(wǎng)殼結(jié)構(gòu)動(dòng)力可靠度的影響. 擬選取的參數(shù)見表1.

為說明本文方法的合理性，在太平洋地震工程研究中心數(shù)據(jù)庫(kù)選取符合規(guī)范反應(yīng)譜（抗震設(shè)防烈度為8度，場(chǎng)地類別為Ⅲ類，設(shè)計(jì)地震分組為第2組）要求的40條天然地震動(dòng)，對(duì)編號(hào)為N60180303的網(wǎng)殼進(jìn)行時(shí)程分析，將結(jié)果與隨機(jī)地震動(dòng)作用下的等價(jià)極值分布函數(shù)（EEVD）曲線進(jìn)行對(duì)比，如圖8所示. 天然地震動(dòng)計(jì)算得到的最大節(jié)點(diǎn)位移最小值和最大值分別為0.091 2 m和0.282 8 m;塑性桿比例最小值和最大值分別為0.152 8和0.506 6. 結(jié)合圖8可知：隨機(jī)地震動(dòng)計(jì)算所得的極值分布可將天然地震動(dòng)計(jì)算結(jié)果涵蓋在其分布范圍內(nèi)，說明本文所采用的隨機(jī)地震動(dòng)在很大程度上可真實(shí)地反映實(shí)際地震作用.

圖9為N60180302、N60180303和N60180304參數(shù)下的EEVD曲線. 由圖可知，無論是從最大節(jié)點(diǎn)位移還是從塑性桿比例的角度，各類場(chǎng)地呈現(xiàn)的結(jié)果一致，Ⅳ類場(chǎng)地EEVD曲線分布更偏右. 圖10給出對(duì)應(yīng)的累積分布函數(shù)（CDF）曲線. 由圖可直觀看出：在保證其它條件一致的情形下，Ⅱ類場(chǎng)地單層球面網(wǎng)殼結(jié)構(gòu)可靠性更好、Ⅲ類場(chǎng)地次之、Ⅳ場(chǎng)地類最差. 一方面，根據(jù)《建筑抗震設(shè)計(jì)規(guī)范》（GB50011-2010）[15]，從Ⅱ類場(chǎng)地到Ⅳ類場(chǎng)地，地基土逐漸偏軟，場(chǎng)地土卓越圓頻率依次為15.708 rad/s、11.424 rad/s和8.378 rad/s. 結(jié)構(gòu)一階、二階圓頻率為9.43 rad/s，正好處于卓越頻率范圍之內(nèi)且與Ⅳ類場(chǎng)地更為接近;另一方面，根據(jù)功率譜分布，從Ⅱ到Ⅳ類場(chǎng)地能量越來越集中，且結(jié)構(gòu)的前兩階主頻率正好處于本文地震動(dòng)功率譜模型的峰值附近，因而對(duì)單層球面網(wǎng)殼結(jié)構(gòu)的響應(yīng)有較大影響.

限于篇幅，針對(duì)矢跨比、初始缺陷，本文僅給出一種情形，其它情形皆以可靠度形式列于表2中. 圖11為N60180303和N60180503兩種參數(shù)下的累積分布函數(shù)圖. 由圖可知，在其它參數(shù)一致的情形下，矢跨比越小，可靠度越高. 圖12為考慮初始幾何缺陷（L/300）的情形，采用一致缺陷模態(tài)法施加初始缺陷，明顯地，缺陷的存在會(huì)顯著降低結(jié)構(gòu)的動(dòng)力可靠度.

基于上述分析結(jié)果，按照文獻(xiàn)[6]提出的強(qiáng)度失效指標(biāo)，分別選取最大節(jié)點(diǎn)位移L/100和塑性桿比例42%為閾值，計(jì)算不同情形下結(jié)構(gòu)動(dòng)力可靠度，如表2所示. 對(duì)比分析可見：采用不同閾值得到的可靠度幾無明顯統(tǒng)一性. 前者從局部最大位移的角度，后者則從整體塑性發(fā)展程度的角度來考察可靠度. 如何建立統(tǒng)一的強(qiáng)度失效指標(biāo)需進(jìn)一步研究.

值得指出的是，本文中可靠度分析與傳統(tǒng)普遍采用的地震易損性分析有著本質(zhì)區(qū)別. 傳統(tǒng)地震易損性分析主要步驟為：選定某一性能水準(zhǔn)，運(yùn)用增量動(dòng)力分析（IDA）方法，對(duì)計(jì)算結(jié)果進(jìn)行曲線擬合得到在不同地震強(qiáng)度下結(jié)構(gòu)失效概率，進(jìn)而得到易損性曲線;本文中可靠度分析刻畫的是：選定某一地震動(dòng)強(qiáng)度，通過隨機(jī)動(dòng)力響應(yīng)分析，運(yùn)用概率密度演化方法得到結(jié)構(gòu)在不同性能水準(zhǔn)下的可靠概率，從而得到概率累積分布（CDF）曲線. 此外，地震易損性分析一般按照不同場(chǎng)地及設(shè)防要求選取數(shù)十條實(shí)際地震記錄進(jìn)行分析，然而選取的地震波可能很難同時(shí)在場(chǎng)地特性、平均反應(yīng)譜等方面很好地滿足設(shè)防要求，并且數(shù)十條波的計(jì)算結(jié)果本質(zhì)上屬于小樣本范疇，計(jì)算結(jié)果精度亦不能保證;而本文采用概率密度演化方法可得到所關(guān)心結(jié)構(gòu)響應(yīng)精細(xì)的概率信息，且與蒙特卡洛模擬結(jié)果擬合較好，可認(rèn)為本文方法較地震易損性分析更能精確刻畫失效概率.

4? ?結(jié)? ?論

本文基于譜表示-隨機(jī)函數(shù)思想建立非平穩(wěn)隨機(jī)地震動(dòng)模型，結(jié)合概率密度演化方法，提出一種新的單層球面網(wǎng)殼結(jié)構(gòu)抗震分析方法，主要結(jié)論如下：

1）引入隨機(jī)函數(shù)，將傳統(tǒng)譜表示方法中的基本隨機(jī)變量數(shù)大大降低，減小問題分析難度;采用一定的修正方法，得到符合規(guī)范要求的隨機(jī)地震動(dòng)，其平均反應(yīng)譜、一階統(tǒng)計(jì)值及二階統(tǒng)計(jì)值均與目標(biāo)值擬合較好.

2）基于概率密度演化方法，對(duì)單層球面網(wǎng)殼結(jié)構(gòu)進(jìn)行隨機(jī)地震響應(yīng)分析. 結(jié)果表明，結(jié)構(gòu)最大位移響應(yīng)均值和標(biāo)準(zhǔn)差與10 000次蒙特卡洛模擬結(jié)果擬合較好，且其概率密度演化過程具有非規(guī)則性，不同于常用的正態(tài)分布假定.

3）基于極值分布，對(duì)單層球面網(wǎng)殼結(jié)構(gòu)進(jìn)行了動(dòng)力可靠度分析.結(jié)果表明：a）場(chǎng)地類別、矢跨比和缺陷對(duì)單層球面網(wǎng)殼結(jié)構(gòu)的動(dòng)力可靠度均有顯著影響，為工程實(shí)踐提供重要的技術(shù)依據(jù);b）基于最大節(jié)點(diǎn)位移和塑性桿比例得到的動(dòng)力可靠度沒有明顯的統(tǒng)一性，如何建立統(tǒng)一的強(qiáng)度失效指標(biāo)是今后單層球面網(wǎng)殼結(jié)構(gòu)抗震可靠度研究的重點(diǎn)之一.

參考文獻(xiàn)

[1]? ?何放龍，馬自克. 強(qiáng)震作用下雙層球面網(wǎng)殼結(jié)構(gòu)非線性動(dòng)力響應(yīng)分析[J]. 湖南大學(xué)學(xué)報(bào)（自然科學(xué)版），2007，34（10）：1—5.

HE F L，MA Z K. Nonlinear analysis of double-layer reticulated spherical dome under strong earthquake motions[J]. Journal of Hunan University（Natural Sciences），2007，34（10）：1—5. （In Chinese）

[2]? ?李永梅，胡琨，張微敬. 考慮損傷累積效應(yīng)的單層球面網(wǎng)殼動(dòng)力穩(wěn)定[J]. 湖南大學(xué)學(xué)報(bào)（自然科學(xué)版），2014，41（6）：16—21.

LI Y M，HU K，ZHANG W J. Dynamic stability of single-layer latticed shells considering damage accumulation effect[J]. Journal of Hunan University（Natural Sciences），2014，41（6）：16—21. （In Chinese）

[3]? ?杜文風(fēng)，高博青，董石麟，等. 一種判定桿系結(jié)構(gòu)動(dòng)力穩(wěn)定的新方法——應(yīng)力變化率法[J]. 浙江大學(xué)學(xué)報(bào)（工學(xué)版），2006，40（3）：506—510.

DU W F，GAO B Q，DONG S L，et al. New judgment method on dynamic instability of truss structures：Stress rate method[J].Journal of Zhejiang University（Engineering Sciences），2006，40（3）：506—510. （In Chinese）

[4]? ?郭海山，沈世釗. 單層網(wǎng)殼結(jié)構(gòu)動(dòng)力穩(wěn)定性分析方法[J]. 建筑結(jié)構(gòu)學(xué)報(bào)，2003，24（3）：1—9.

GUO H S，SHEN S Z. Analysis method of dynamic stability of single-layer reticulated domes[J]. Journal of Building Structures，2003，24（3）：1—9. （In Chinese）

[5]? ?XU J，LI J. An energetic criterion for dynamic instability of struc-tures under arbitrary excitations[J]. International Journal of Structural Stability & Dynamics，2015，15（2）：1—32.

[6]? ?沈世釗，支旭東. 球面網(wǎng)殼結(jié)構(gòu)在強(qiáng)震下的失效機(jī)理[J]. 土木工程學(xué)報(bào)，2005，38（1）：11—20.

SHEN S Z，ZHI X D. Failure mechanism of recticular shells subjected to dynamic actions[J]. China Civil Engineering Journal，2005，38（1）：11—20. （In Chinese）

[7]? ?杜文風(fēng)，高博青，董石麟. 單層球面網(wǎng)殼結(jié)構(gòu)動(dòng)力強(qiáng)度破壞的雙控準(zhǔn)則[J]. 浙江大學(xué)學(xué)報(bào)（工學(xué)版），2007，41（11）：1916—1920.

DU W F，GAO B Q，DONG S L. Double-control criterion of dynamical strength failure for single layer latticed shells[J]. Journal of Zhejiang University（Engineering Sciences），2007，41（11）：1916—1920. （In Chinese）

[8]? ?ZHI X D，F(xiàn)AN F，SHEN S Z. Failure mechanisms of single-layer reticulated domes subjected to earthquakes[J]. Journal of the International Association for Shell & Spatial Structures，2007，48（153）：29—44.

[9]? ?林家浩. 隨機(jī)地震響應(yīng)的確定性算法[J].? 地震工程與工程振動(dòng)，1985（1）：91—96.

LIN J H. A deterministic method for the computation of stochastic earthquake response[J]. Earthquake Engineering & Engineering Vibration，1985（1）：91—96. （In Chinese）

[10]? 薛素鐸，曹資，王雪生，等. 多維地震作用下網(wǎng)殼結(jié)構(gòu)的隨機(jī)分析方法[J].? 空間結(jié)構(gòu)，2002，8（1）：44—51.

XUE S D，CAO Z，WANG X S，et al. Random analysis method for lattice shells under multiple earthquake excitations[J]. Spatial Structures，2002，8（1）：44—51. （In Chinese）

[11] 王雪生. 網(wǎng)殼結(jié)構(gòu)多維多點(diǎn)非平穩(wěn)隨機(jī)地震響應(yīng)分析研究[D]. 北京：北京工業(yè)大學(xué)建筑工程學(xué)院，2002：15—48.

WANG X S. Study on multi-dimensional and multi-support non-stationary random seismic response of reticulated shell structure [D]. Beijing：College of Architecture and Civil Engineering，Beijing University of Technology，2002：15—48. （In Chinese）

[12]? LI J，CHEN J B. Stochastic dynamics of structures [M]. Singapore：John Wiley & Sons，2010：213—283.

[13]? LIU Z J，LIU W，PENG Y. Random function based spectral represe-ntation of stationary and non-stationary stochastic processes [J]. Pro-babilistic Engineering Mechanics，2016，45：115—126.

[14]? CLOUGH R W，PENZIEN J，GRIFFIN D S. Dynamics of structures [M]. New Jersey：Prentice Hall，1990：213—225.

[15]? GB50011—2010? ?建筑抗震設(shè)計(jì)規(guī)范[S]. 北京：中國(guó)建筑工業(yè)出版社，2010：16—33.

GB50011—2010? ?Code for seismic design of buildings[S].Beijing：China Architecture & Building Press，2010：16—33. （In Chinese）

[16]? 劉章軍，劉增輝，劉威. 全非平穩(wěn)地震動(dòng)過程的概率模型及反應(yīng)譜擬合[J]. 振動(dòng)與沖擊，2017，36（2）：32—38.

LIU Z J，LIU Z H，LIU W. Probability model of fully non-stationary ground motion with the target response spectrum compatible[J]. Journal of Vibration & Shock，2017，36（2）：32—38. （In Chinese）

[17]? LI J，CHEN J B. The number theoretical method in response analysis of nonlinear stochastic structures[J]. Computational Mechanics，2007，39（6）：693—708.

[18]? 袁一凡，田啟文. 工程地震學(xué)[M]. 北京：地震出版社，2012：204—208.

YUAN Y F，TIAN Q W. Engineering seismology[M]. Beijing：Seismological Press，2012：204—208. （In Chinese）

[19]? 楊慶山，田玉基. 地震地面運(yùn)動(dòng)及其人工合成[M]. 北京：科學(xué)出版社，2014：69—72.

YANG Q S，TIAN Y J. Earthquake ground motion and artificial synthesis [M]. Beijing：Science Press，2014：69—72. （In Chinese）

[20]? LI J，CHEN J B，F(xiàn)AN W L. The equivalent extreme-value event and evaluation of the structural system reliability[J]. Structural Safety，2007，29（2）：112—131.