楊思昭，王憲杰,2

(1 云南大學(xué)建筑與規(guī)劃學(xué)院，昆明 650500；2 浙江大學(xué)建筑工程學(xué)院，杭州 310058)

0 概述

地震易損性分析能夠預(yù)測在不同強度地震作用下，建筑結(jié)構(gòu)達(dá)到或超越不同性能水準(zhǔn)的概率[1-2]。但傳統(tǒng)的易損性分析主要是靜態(tài)的，未考慮材料性能衰退對結(jié)構(gòu)抗震性能的影響。事實上，結(jié)構(gòu)在全壽命周期內(nèi)，隨著服役齡期的增加，由于受力構(gòu)件有效截面面積的減小和材料力學(xué)性能的退化，結(jié)構(gòu)在同一地面運動加速度作用下發(fā)生破壞的概率增加[3-4]。因此，在進行地震易損性評估的研究中，考慮材料時變特性的影響更能真實反映不同服役期建筑結(jié)構(gòu)的地震易損性。

但由于耐久性研究存在時間跨度長、影響因素復(fù)雜、不確定性大等問題，目前尚未有統(tǒng)一、可靠的理論模型來刻畫材料、構(gòu)件和結(jié)構(gòu)層面的長周期、隨機衰變特性。與此同時，在實際工程的檢測、監(jiān)測和檢修等作業(yè)中，已取得了服役結(jié)構(gòu)的諸多性能參數(shù)。但實際工程中的實測數(shù)據(jù)(特別是有損檢測)由于樣本數(shù)量少而缺乏足夠的代表性，僅依靠小樣本檢測數(shù)據(jù)評估建筑結(jié)構(gòu)的實際工作性能缺乏足夠的完備性[5]。結(jié)合現(xiàn)有的實測數(shù)據(jù)，引入信息更新理論對初始設(shè)計模型不斷修正和完善，可實現(xiàn)對具體結(jié)構(gòu)性能退化的合理預(yù)測。針對實時信息更新與系統(tǒng)吸收，經(jīng)典貝葉斯方法可以有效實現(xiàn)對先驗分布的相關(guān)參數(shù)進行重新估計和檢驗，得到系統(tǒng)后驗分布數(shù)據(jù)[6-7]。因此，采用貝葉斯方法進行結(jié)構(gòu)耐久性分析既能利用多年來耐久性基礎(chǔ)研究對結(jié)構(gòu)內(nèi)在規(guī)律的理性認(rèn)識，又能考慮結(jié)構(gòu)的實際情況[8]，從而準(zhǔn)確刻化結(jié)構(gòu)實時動態(tài)性能。

退化結(jié)構(gòu)的時變結(jié)構(gòu)抗力是非平穩(wěn)隨機過程，要求必須以隨機過程理論為基礎(chǔ)，研究隨機荷載效應(yīng)超過隨機抗力的概率[9]。即通過理論分析研究結(jié)構(gòu)的時變易損性時，同時需要考慮結(jié)構(gòu)自身參數(shù)的不確定性和外界激勵的隨機性。此外，在傳統(tǒng)的地震易損性研究中，為了便于問題的求解，通常需要基于以下幾點假定進行分析[2-4,10]：結(jié)構(gòu)地震需求參數(shù)和給定地震動水平作用之間服從冪指數(shù)回歸關(guān)系；在給定地震動水平作用下，結(jié)構(gòu)地震需求參數(shù)服從對數(shù)正態(tài)分布；結(jié)構(gòu)的概率抗震能力服從對數(shù)正態(tài)分布。顯然，基于此求解得到的地震易損性結(jié)果只是經(jīng)驗估計值，不能反映結(jié)構(gòu)響應(yīng)的真實概率分布。李杰[11]提出的概率密度演化理論不僅為復(fù)合隨機問題的求解提供了理論依據(jù)，而且可以分析得到滿足真實分布的結(jié)構(gòu)響應(yīng)概率密度分布函數(shù)，為精細(xì)化評估既有結(jié)構(gòu)的地震易損性提供了新思路。

本文通過總結(jié)建筑材料耐久性的相關(guān)研究成果，建立材料力學(xué)性能退化模型；引入不同服役齡期材料力學(xué)性能的實測數(shù)據(jù)，并采用貝葉斯方法進行數(shù)據(jù)更新?；诟碌牟牧闲阅芡嘶Ｐ团c概率密度演化理論對某7層帶內(nèi)廊式對稱三跨RC框架結(jié)構(gòu)的地震易損性進行分析，得到RC框架結(jié)構(gòu)抗震性能隨服役齡期的一般變化規(guī)律。

1 材料性能時變特性分析與貝葉斯方法

1.1 隨機參數(shù)更新貝葉斯方法

已有研究[12-13]表明，一般大氣環(huán)境中既有建筑物混凝土的抗壓強度服從正態(tài)分布，其均值和標(biāo)準(zhǔn)差是服役時間t的函數(shù)。采用共軛先驗分布確定均值μ的后驗分布。當(dāng)方差σ2已知時，均值μ的先驗分布可視為服從正態(tài)分布N(θ,τ2)，即：

(1)

式中：θ，τ為超參數(shù)，可根據(jù)先驗預(yù)測公式進行確定。

先驗分布反映對研究對象內(nèi)在規(guī)律的理性認(rèn)識，不可避免的存在主觀判斷，尚需考慮其實際情況。設(shè)t時刻材料物理性能參數(shù)的實測樣本值{x1,x2,x3,…,xn}取自總體樣本N(μ,σ2)，其中方差σ2已知。則此樣本x的似然函數(shù)可表示為：

(2)

進而，由貝葉斯方法可將結(jié)構(gòu)物理參數(shù)μ的后驗分布表示為：

(3)

將式(1)，(2)帶入式(3)，經(jīng)過一定的數(shù)學(xué)推導(dǎo)可得μ的后驗分布為：

(4)

(5)

圖1給出了采用貝葉斯方法對材料性能參數(shù)進行信息更新的一般流程。

圖1 材料性能參數(shù)貝葉斯動態(tài)信息更新過程

1.2 混凝土經(jīng)時抗壓強度預(yù)測模型

一般而言，混凝土經(jīng)時抗壓強度的均值在初期隨時間增大而增大，但增長速率逐漸減緩，在后期則隨時間增大而降低；標(biāo)準(zhǔn)差隨時間增大而增大，且不受均值影響。牛荻濤等[12]通過分析長期暴露試驗和實測的服役混凝土強度，建立了一般大氣環(huán)境中混凝土強度均值μf(t)和標(biāo)準(zhǔn)差σf(t)的經(jīng)時變換模型：

μf(t)=1.452 9exp[-0.024 6(lnt-

1.715 4)2]·μ′0

(6)

σf(t)=(0.030 5t+1.236 8)·σ′0

(7)

式中μ′0和σ′0分別為混凝土28d強度均值和標(biāo)準(zhǔn)差。

本文依據(jù)《混凝土結(jié)構(gòu)設(shè)計規(guī)范》(GB 50010—2010)給出的混凝土抗壓強度均值和標(biāo)準(zhǔn)值之間的關(guān)系，將混凝土強度均值μf(t)的經(jīng)時變化模型轉(zhuǎn)換為混凝土強度標(biāo)準(zhǔn)值fcu(t)的經(jīng)時變化模型：

fcu(t)=μf(t)-1.645σf(t)

(8)

1.3 鋼筋材料性能時變模型

相較于混凝土經(jīng)時抗壓強度的研究，鋼筋經(jīng)時銹蝕深度已有較為成熟的研究成果[14]，且已在規(guī)范中進行體現(xiàn)[15]。處于一般大氣環(huán)境中的RC框架結(jié)構(gòu)，鋼筋銹蝕開始時間ti可按下式進行計算：

(9)

式中：C為混凝土保護層厚度；k為碳化系數(shù)；x0為碳化殘量。具體參數(shù)取值詳見規(guī)范[15]。

鋼筋銹蝕是一個電化學(xué)腐蝕過程，其銹蝕速率與鋼筋表面的含氧量有關(guān)，而含氧量取決于混凝土質(zhì)量、保護層厚度和環(huán)境條件。在溫度為20℃，相對濕度為75%的典型環(huán)境中，鋼筋初始銹蝕電流密度icoor(1)為[14]：

(10)

式中w/c為水灰比，當(dāng)已知混凝土抗壓強度標(biāo)準(zhǔn)值fcu(t)時，水灰比可按下式計算：

(11)

進而，鋼筋開始銹蝕之后的某時刻tc，鋼筋銹蝕電流密度icoor(tc)為：

(12)

通常情況下，當(dāng)銹蝕電流密度icoor為1μA/cm2時，相應(yīng)的鋼筋銹蝕深度速率為11.6μm/年。

2 經(jīng)時材料本構(gòu)關(guān)系定義

2.1 經(jīng)時混凝土的本構(gòu)關(guān)系

相較于大量已有的時不變或素混凝土本構(gòu)關(guān)系模型，考慮銹蝕箍筋約束影響建立的經(jīng)時混凝土本構(gòu)關(guān)系更符合實際服役情況。在考慮了體積配箍率、箍筋屈服強度等影響因素建立的Kent-Scott-Park模型基礎(chǔ)上，通過引入峰值應(yīng)力和應(yīng)變軟化修正系數(shù)[16]建立了銹蝕矩形箍筋約束混凝土本構(gòu)關(guān)系模型，如式(13)所示：

(13)

式中：σc和εc分別為銹蝕箍筋約束混凝土的應(yīng)力和應(yīng)變；εcc為銹蝕箍筋約束混凝土的峰值點應(yīng)變，εcc=0.002c2K；εcu為銹蝕箍筋約束混凝土的極限應(yīng)變；c1，c2和c3分別為非銹蝕箍筋約束混凝土的峰值應(yīng)力、峰值點應(yīng)變和應(yīng)變軟化調(diào)整系數(shù)，按式(14)確定；k1，k2分別為考慮箍筋銹蝕影響的峰值應(yīng)力和應(yīng)變軟化修正系數(shù)，按式(15)確定；K為考慮箍筋約束影響的混凝土強度增強系數(shù)，按式(16)確定；Zm為應(yīng)變軟化段斜率，按式(17)確定；fc(t)為混凝土的軸心抗壓強度。

c1=0.157 0λt+0.963 4

(14-1)

c2=-1.556 8λt+1.232 0

(14-2)

c3=17.744 0λt+0.974 2

(14-3)

k1=(0.291 4λt-0.179 1)ω+1.0

(15-1)

k2=exp[(-1.585 2λt+0.960 8)ω]

(15-2)

(16)

(17)

式中：λt為箍筋特征值；ω為銹脹裂縫寬度；ρs為體積配箍率；fyh為箍筋屈服強度；h″為從箍筋外邊緣算起的核心混凝土寬度；sh為箍筋間距。

2.2 銹蝕鋼筋的本構(gòu)關(guān)系

實際工程中的銹蝕鋼筋，當(dāng)其銹蝕率ηs<80%時，本構(gòu)關(guān)系可按式(18)進行確定[17]：

(18)

式中：σs和εs分別為銹蝕鋼筋的應(yīng)力和應(yīng)變；Es0為未銹蝕鋼筋的彈性模量；fyc，fuc分別為銹蝕鋼筋的屈服強度和極限強度，按式(19)確定；εsyc，εshc和εsuc分別為銹蝕鋼筋的屈服應(yīng)變、強化應(yīng)變和極限應(yīng)變，按式(20)確定。

(19-1)

(19-2)

(20-1)

(20-2)

εsuc=e-2.501ηsεsu0

(20-3)

式中：fy0，fu0分別為未銹蝕鋼筋的屈服強度和極限強度；εsh0，εsu0分別為未銹蝕鋼筋的強化應(yīng)變和極限應(yīng)變；ηs,cr為鋼筋屈服平臺消失時的臨界銹蝕率，帶肋鋼筋取20%，光圓鋼筋取10%。

3 基于PDEM的結(jié)構(gòu)時變易損性分析

時變易損性表征既有建筑結(jié)構(gòu)在不同服役時期遭遇可能強度地震作用時，結(jié)構(gòu)達(dá)到或超越各級性能水準(zhǔn)的概率?？杀硎緸椋?/p>

FR(a,t)=Pf[PL|SA=a,T=t]

(21)

式中：FR(a,t)為結(jié)構(gòu)的時變易損性；Pf為失效概率；PL為結(jié)構(gòu)的性能水平；SA為地震動強度指標(biāo)，SA=a；T為結(jié)構(gòu)服役期，T=t。

顯然，結(jié)構(gòu)響應(yīng)θ(SA,T)達(dá)到或超過目標(biāo)性能水準(zhǔn)θc的概率可表示為：

Pf[SA,T]=Pr[θc≤θ(SA,T)]

(22)

通過以上分析可知，劃分合理可信的目標(biāo)性能水準(zhǔn)和較為精細(xì)化求解結(jié)構(gòu)在不同地震動強度和不同服役年限條件下的結(jié)構(gòu)響應(yīng)是時變易損性分析的基礎(chǔ)，以下分別就這兩方面進行討論。

3.1 概率抗震能力分析

《建筑抗震設(shè)計規(guī)范》(GB 50011—2010)蘊含豐富的基于性能抗震設(shè)計思想，其目標(biāo)性能的選定依托于對性能水平的合理劃分。針對RC框架結(jié)構(gòu)，可采用結(jié)構(gòu)層間位移角來定義結(jié)構(gòu)的性能水平，將RC框架結(jié)構(gòu)不同破壞程度對應(yīng)的最大層間位移角限值列于表1。在此指出，為了便于問題的考慮，本文將各目標(biāo)破壞狀態(tài)限值視為確定量進行考慮。

破壞等級與最大層間位移角限值的關(guān)系 表1

3.2 基于PDEM的概率地震需求分析

傳統(tǒng)的易損性分析方法通常假定概率地震需求參數(shù)服從對數(shù)正態(tài)分布，是基于經(jīng)驗的一種近似分析方法。區(qū)別于傳統(tǒng)方法，基于概率守恒原理提出的概率密度演化理論[11,18]可以準(zhǔn)確求取結(jié)構(gòu)動力響應(yīng)的概率密度函數(shù)及其演化過程，為結(jié)構(gòu)的概率地震需求分析提供了新思路。

假設(shè)結(jié)構(gòu)反應(yīng)需求參數(shù)Z(t)為所考察的狀態(tài)量，對于某一物理解答存在、唯一且連續(xù)依賴于初始條件的概率保守系統(tǒng)(Z(t),Ψ)，其中Ψ=(Ψ1,Ψ2,…,Ψs)為隨機參數(shù)向量，s為隨機變量的個數(shù)，其聯(lián)合概率密度函數(shù)為PZΨ(z,ψ,t)，根據(jù)概率守恒原理的隨機事件描述，則有：

(23)

式中：D/Dt(·)為全導(dǎo)數(shù)；Ωt×ΩΨ為初始狀態(tài)空間在t時刻對應(yīng)的解區(qū)域，對其經(jīng)過一系列的數(shù)學(xué)推導(dǎo)，即求得廣義概率密度演化方程。

(24)

設(shè)初始條件為：

pZΨ(z,ψ,t0)=δ(z-z0)pΨ(ψ)

(25)

式中：pΨ(ψ)為Ψ的聯(lián)合概率密度函數(shù)；δ(·)為Dirac函數(shù)；z0為確定性初始值。

進而，采用TVD格式的有限差分法可求得Z(t)的時變概率密度函數(shù)：

(26)

式中ΩΨ為Ψ的分布空間。

顯然，所求的時變概率密度函數(shù)pZ(z,t)即反映了給定地震動強度與結(jié)構(gòu)響應(yīng)之間的概率分布關(guān)系。

3.3 概率密度演化方程求解的TVD格式

為方便討論，取式(24)中的任一偏微分方程，記為：

(27)

針對形如式(27)的一維偏微分方程，采用有限差分法求解是行之有效的思路之一[19],且目前已發(fā)展了單邊差分格式、Lax-Wendroff格式和TVD格式三種求解格式。其中，單邊差分格式精度較低，耗散效應(yīng)明顯；Lax-Wendroff格式不能保證結(jié)果的非負(fù)性，色散效應(yīng)突出；基于上述兩種基本格式，并施加適當(dāng)形式的通量限制器可求得具有總變差不增的TVD格式為：

(28)

(29)

式中：φ0(r)=max[0,min(2r,1),min(r,2)]；

3.4 基本求解流程

目前，國內(nèi)外已發(fā)展了多種結(jié)構(gòu)易損性分析方法[20]。為了較為精細(xì)地考慮既有建筑結(jié)構(gòu)的抗震性能，本文給出了采用貝葉斯方法進行材料參數(shù)更新，并基于概率密度演化理論研究信息更新結(jié)構(gòu)的時變易損性分析方法，具體步驟包括如下。

(1)材料性能參數(shù)的貝葉斯更新：1)選取合理的材料物理性能經(jīng)時變化模型作為先驗預(yù)測信息；2)實測不同服役齡期建筑材料的實時物理性能；3)采用貝葉斯方法對先驗?zāi)Ｐ瓦M行信息更新，獲得更符合實際情況的材料性能衰退模型。

(2)隨機非線性力學(xué)模型的建立：1)建立不同服役齡期條件下結(jié)構(gòu)的非線性力學(xué)模型；2)選取一系列滿足場地條件的不同震級、不同斷層距的地震動時程記錄；3)分析結(jié)構(gòu)參數(shù)、外部激勵存在的隨機性，確定隨機變量的個數(shù)，進而采用數(shù)論選點法進行離散代表點的選取和賦得概率的計算。

(3)時變易損性分析：1)采用TVD格式的有限差分法求解廣義概率密度演化方程，獲得不同服役齡期、不同地震動強度作用下結(jié)構(gòu)需求參數(shù)的概率密度分布函數(shù)；2)定義結(jié)構(gòu)的破壞等級及相應(yīng)的量化指標(biāo)；3)計算不同服役齡期、不同強度地震作用下結(jié)構(gòu)響應(yīng)超過不同破壞等級限值的條件概率，進而繪制以地震動參數(shù)為變量的時變易損性曲線(面)。

4 算例分析

4.1 RC框架結(jié)構(gòu)設(shè)計

本文以7層帶內(nèi)廊式對稱三跨鋼筋混凝土框架結(jié)構(gòu)辦公樓為例。結(jié)構(gòu)設(shè)計使用年限為50年，建筑抗震設(shè)防類別為丙類；抗震設(shè)防烈度為8度(0.2g)，Ⅱ類場地，設(shè)計地震分組第二組；地面粗糙度類別C類，基本風(fēng)壓0.3kN/m2。

框架結(jié)構(gòu)平面布置如圖2所示。結(jié)構(gòu)底層層高3.9m,其余層層高3.3m。柱截面尺寸為0.55m×0.55m，梁截面尺寸為0.50m×0.20m，板厚為100mm。梁板柱混凝土強度等級均為C30。梁柱主筋采用HRB400級鋼筋，梁柱箍筋、板采用HRB335級鋼筋。梁柱混凝土保護層厚度20mm。填充墻為190mm厚的混凝土空心砌塊，活荷載標(biāo)準(zhǔn)值按《建筑結(jié)構(gòu)荷載規(guī)范》(GB 50009—2012)規(guī)定取值。

圖2 框架結(jié)構(gòu)平面布置圖

采用YJK軟件進行荷載組合、內(nèi)力分析和截面配筋等初步設(shè)計工作。其中，結(jié)構(gòu)的柱配筋見圖3。

圖3 結(jié)構(gòu)的柱配筋示意圖

4.2 混凝土經(jīng)時強度更新的貝葉斯方法

文獻[13]采用回彈法和鉆芯法實測一般大氣環(huán)境中不同服役年限民用建筑物的混凝土抗壓強度。由于鉆芯法更加準(zhǔn)確，本文選取基于此檢測得到的服役混凝土相對強度樣本，如表2所示。當(dāng)已知混凝土的立方體抗壓強度標(biāo)準(zhǔn)值時，乘以相應(yīng)的相對強度即得到混凝土?xí)r變強度標(biāo)準(zhǔn)值。

實測服役混凝土相對強度樣本值 表2

將式(6)，(7)視為先驗預(yù)測模型，并通過式(8)建立混凝土經(jīng)時抗壓強度平均值和標(biāo)準(zhǔn)值的換算關(guān)系。圖4為先驗混凝土經(jīng)時相對強度標(biāo)準(zhǔn)值的預(yù)測模型和實測樣本值的對比。顯然，采用牛荻濤模型[12]預(yù)測的混凝土經(jīng)時強度標(biāo)準(zhǔn)值高于實測樣本值，即該預(yù)測模型不能很好地預(yù)測混凝土經(jīng)時強度的真實情況。若直接采用該模型進行結(jié)構(gòu)非線性響應(yīng)分析，預(yù)測結(jié)果將比真實情況更大，不利于對實際服役結(jié)構(gòu)進行檢修決策和優(yōu)化加固設(shè)計等作業(yè)。

圖4 先驗混凝土經(jīng)時強度預(yù)測值與試驗值對比

根據(jù)貝葉斯方法，將表2中的實測數(shù)據(jù)樣本值視為似然函數(shù)，并用其將先驗預(yù)測模型進行多次修正，得到信息更新后混凝土經(jīng)時強度的預(yù)測值，如圖5所示。分析可知：貝葉斯方法可有效實現(xiàn)先驗分布和實測數(shù)據(jù)的綜合。一方面，更新后的模型曲線與先驗預(yù)測曲線的發(fā)展趨勢存在相似的變化規(guī)律，即貝葉斯方法繼承了對先驗分布理性規(guī)律的認(rèn)識；另一方面，貝葉斯模型曲線與先驗分布曲線存在一定的偏移，經(jīng)過若干次貝葉斯更新即可使預(yù)測結(jié)果改善，使更新值更加接近實測樣本值，即貝葉斯方法能充分考慮結(jié)構(gòu)的實際服役情況。

將混凝土的經(jīng)時抗壓強度fcu(t)視為滿足正態(tài)分布的隨機變量，為了充分考慮材料的實際服役情況，采用貝葉斯方法對其均值進行參數(shù)更新(圖5)；而將鋼筋的時變銹蝕深度視為確定性變量。故本文考察的7層RC框架結(jié)構(gòu)中，將每一層混凝土的時變強度視為一個隨機變量，共包含7個隨機變量。表3列出了服役齡期條件下材料參數(shù)的基本時變信息。鋼筋和混凝土采用第2節(jié)定義的經(jīng)時材料本構(gòu)關(guān)系模型。

圖5 混凝土經(jīng)時相對強度均值的貝葉斯更新

不同服役齡期下的材料性能 表3

4.3 地震波的選取

地震波的合理選取直接影響結(jié)構(gòu)易損性的分析效果。本文基于PEER地震動數(shù)據(jù)庫，取場地剪切波速Vs,30為260～510m/s，對應(yīng)中國的Ⅱ類場地[21]。選取的20條地震動記錄均勻分布于以下5個地震動條帶中[22]：1)SMSR(小震級、短距離)：5.5≤M≤6.5，15km≤R≤30km；2)SMLR(小震級、長距離)：5.5≤M≤6.5，30km≤R≤50km；3)LMSR(大震級、短距離)：6.5≤M≤7.5，15km≤R≤30km；4)LMLR(大震級、長距離)：6.5≤M≤7.5，30km≤R≤50km；5)NEAR(近場地震)：6.5≤M≤7.5，0km≤R≤15km。其中M為震級，R為斷層距，SA1為結(jié)構(gòu)第一自振周期對應(yīng)的加速度反應(yīng)譜值。所選的地震波信息見表4。

地震動記錄 表4

4.4 時變易損性分析

基于最大層間位移角的概率密度分布函數(shù)求得結(jié)構(gòu)在不同服役年限、不同地震動水平、不同性能等級作用下的失效概率，采用分段三次Hermite函數(shù)擬合得到服役結(jié)構(gòu)在不同服役齡期條件下地震易損性曲線如圖6所示。為了更加直觀地反映服役期對結(jié)構(gòu)失效概率的影響，進而擬合得到三維時變易損性曲面如圖7所示。

圖7 時變易損性曲面

從圖6可知，隨著結(jié)構(gòu)破壞程度的增加，結(jié)構(gòu)易損性曲線在不同服役齡期之間的差異越大。如：當(dāng)選用結(jié)構(gòu)最大層間位移角為1/550時(完好)，易損性曲線在不同服役齡期之間差異很??；而選用結(jié)構(gòu)最大層間位移角為1/60時(不嚴(yán)重破壞)，易損性曲線在不同服役齡期之間差異性較大。即當(dāng)處于不同服役齡期的既有建筑結(jié)構(gòu)遭遇“小震”作用時，其易損性曲線不會存在過大的偏差；而當(dāng)其遭遇“大震”作用時，材料的實時力學(xué)性能很大程度上決定了結(jié)構(gòu)的安全性。另一方面，結(jié)合表3可知，結(jié)構(gòu)服役齡期分別為14年和20年時，鋼筋均未發(fā)生銹蝕，而混凝土強度均值處于相當(dāng)水平；而結(jié)構(gòu)服役齡期分別為1年和50年時，50年混凝土的相對強度略大于1年時的相對強度，且1年時鋼筋未發(fā)生銹蝕，而50年時鋼筋則發(fā)生了程度不大的銹蝕。故而，在本文考察的6個服役齡期內(nèi)，結(jié)構(gòu)服役1年和50年、14年和20年生成的地震易損性曲線是大致一致的，即結(jié)構(gòu)抗震性能相似。

圖6 不同服役齡期條件的地震易損性曲線

綜合分析圖6和圖7的結(jié)構(gòu)時變易損性變化規(guī)律，可以發(fā)現(xiàn)：結(jié)構(gòu)在服役初期，由于混凝土徐變、硬化等因素的影響，其強度處于增強的過程，在同一地面運動加速度作用下，其抗震性能隨著服役齡期的增加而提高。結(jié)構(gòu)服役齡期達(dá)10年左右，混凝土強度發(fā)展到較為穩(wěn)定的階段，并在今后10～20年之間處于平穩(wěn)過程，該時期結(jié)構(gòu)的抗震性能最好。結(jié)構(gòu)服役齡期達(dá)25～30年時，混凝土強度開始降低，且鋼筋開始銹蝕，其抗震性能開始呈現(xiàn)衰減的趨勢。結(jié)構(gòu)服役齡期為50～60年時(超過結(jié)構(gòu)設(shè)計使用年限)，混凝土強度呈現(xiàn)較為明顯的降低趨勢，且鋼筋銹蝕仍持續(xù)發(fā)展，結(jié)構(gòu)的抗震性能亦有一定程度的降低。

此外，可以注意到，傳統(tǒng)的易損性分析方法通常假定結(jié)構(gòu)地震需求參數(shù)服從對數(shù)正態(tài)分布，所得的分析結(jié)果只是經(jīng)驗近似值。而采用概率密度演化理論可直接分析求得地震需求參數(shù)的真實分布模式，其計算結(jié)果相較于傳統(tǒng)方法應(yīng)更加可靠。限于篇幅，圖8僅給出了結(jié)構(gòu)在服役齡期1年時分別采用本文選用方法與傳統(tǒng)分析方法所得的地震易損性曲線。顯然，當(dāng)?shù)卣鹱饔幂^小時，兩種分析方法的計算結(jié)果差異不大；而當(dāng)?shù)卣鹱饔幂^大時，選用本文方法所得的失效概率略大于傳統(tǒng)方法。即當(dāng)?shù)卣鹱饔幂^小時，可選用傳統(tǒng)方法進行結(jié)構(gòu)抗震性能分析，而當(dāng)?shù)卣鹱饔幂^大時，傳統(tǒng)分析方法所求得的計算結(jié)果偏于不安全。

圖8 服役齡期1年時地震易損性計算結(jié)果對比

5 結(jié)論

本文考慮材料性能的時變規(guī)律對結(jié)構(gòu)抗震性能的影響，采用實測數(shù)據(jù)進行材料信息更新，對一考慮隨機參數(shù)的7層RC框架結(jié)構(gòu)進行時變易損性分析，結(jié)論如下：

(1)貝葉斯方法可有效實現(xiàn)先驗分布和實測數(shù)據(jù)的綜合，經(jīng)過更新的貝葉斯數(shù)據(jù)既可保持已有理性認(rèn)識的變化規(guī)律，又能更加接近實際服役情況。

(2)RC框架結(jié)構(gòu)的抗震性能與材料的經(jīng)時性能有關(guān)，前期抗震性能隨著服役齡期的增加而提高，之后一段時期內(nèi)保持在較為穩(wěn)定的階段，當(dāng)服役齡期達(dá)到30年前后，由于材料性能的退化，結(jié)構(gòu)的抗震性能隨服役齡期的增加呈現(xiàn)降低的趨勢。

(3)結(jié)構(gòu)的時變易損性曲線在不同服役期是不一致的，且與劃分的性能水平有關(guān)。隨著結(jié)構(gòu)破壞等級的增加，結(jié)構(gòu)的易損性曲線在不同服役年限之間存在的差異越大。

(4)當(dāng)?shù)卣鹱饔幂^小時，采用本文方法與傳統(tǒng)方法求得的結(jié)構(gòu)失效概率差異不大；而當(dāng)?shù)卣鹱饔幂^大時，本文方法所得的結(jié)構(gòu)失效概率略大于傳統(tǒng)方法。

利用本文所提方法，可以求得滿足實際概率分布的服役結(jié)構(gòu)時變易損性，研究成果對既有建筑結(jié)構(gòu)的檢修決策和優(yōu)化加固提供參考，具有良好的工程實用性。