王洋洋,景月嶺,黃友保,劉丙瑞,周召虎,何李浩

(1.合肥工業(yè)大學(xué) 土木工程結(jié)構(gòu)與材料安徽省重點(diǎn)實(shí)驗(yàn)室,安徽 合肥 230009;2.合肥工業(yè)大學(xué) 土木與水利工程學(xué)院,安徽 合肥 230009;3.安徽水安建設(shè)集團(tuán)股份有限公司,安徽 合肥230601)

0 引言

地震作用下的結(jié)構(gòu)易損性分析起源于對核電站的研究,引入土木工程領(lǐng)域后得到了充分的發(fā)展,取得了一系列的研究成果[1-6]。目前來看,水利工程領(lǐng)域的地震易損性分析研究較少[7-9],且多側(cè)重于結(jié)構(gòu)類型的多樣化和構(gòu)件層次的精細(xì)化建模等方面,并且簡化了外源地震波的輸入模式,即通過在結(jié)構(gòu)底部直接輸入地震波的方式模擬地震動能量的輸入,忽略了地震波的場地效應(yīng)和土-結(jié)構(gòu)動力相互作用對結(jié)構(gòu)易損性分析結(jié)果的影響。最近Samanta[10]對比研究了高中低三種不同樓層的房屋在軟土基礎(chǔ)上的易損性,研究發(fā)現(xiàn)考慮土-結(jié)構(gòu)動力相互作用會影響建筑物的動力響應(yīng)和破損概率。黃博等[11]研究發(fā)現(xiàn),局部河谷地形和土-結(jié)構(gòu)動力相互作用對多層框架結(jié)構(gòu)的動力響應(yīng)和易損性曲線影響較大。但其考慮的是剪切波垂直入射的情形,沒有考慮近場地震動斜入射和地表非一致性運(yùn)動。

綜上所述,現(xiàn)有的結(jié)構(gòu)易損性研究還未將地震波斜入射影響考慮在內(nèi),由于地震波斜入射輸入會在地表產(chǎn)生放大效應(yīng)。為了以地表實(shí)測地震動分量為輸入目標(biāo),本文基于斜入射SV波和P波波場編制了考慮地表實(shí)測地震動水平分量的地震波入射程序。程序可以在斜入射體系下實(shí)現(xiàn)到達(dá)地表的地震動為實(shí)測地震動的目標(biāo)。接著采用增量動力學(xué)分析法(IDA),選取太平洋地震工程研究中心(PEER)數(shù)據(jù)庫中26條實(shí)測地震動的水平向分量為目標(biāo)地震動,以我國西南地區(qū)某壩后式廠房為實(shí)例,計算得到了該壩后式廠房上下游墻體易損性曲線,并對比分析了與傳統(tǒng)垂直入射模式得到的易損性曲線的差異性。

1 基本理論

1.1 基于斜入射波場的地表實(shí)測地震動模擬

一般考慮土-結(jié)構(gòu)動力相互作用時會截取結(jié)構(gòu)周圍的部分地基,并在其周圍施加一定的人工邊界以模擬半無限地基的輻射阻尼效應(yīng)。為了在大型通用有限元軟件ANSYS中實(shí)現(xiàn)地震波的斜入射輸入,本文采用黏彈性邊界[12-13]模擬半無限地基的輻射阻尼作用。地震波通過等效節(jié)點(diǎn)力的方式施加到黏彈性邊界上,由此實(shí)現(xiàn)地震波的近場入射模擬。此前人們多認(rèn)為地震波通過近場基礎(chǔ)入射到自由面是垂直傳播的,近年來國內(nèi)外學(xué)者發(fā)現(xiàn)近場地震波也有一定的比例是斜入射傳播到地表的。為了建立反映地表非一致性運(yùn)動的地震動場,考慮斜入射地震波與地表實(shí)測地震動之間的聯(lián)系,苑舉衛(wèi)等[14-15]通過固體介質(zhì)中波場疊加技術(shù),將地表地震動分量看作是斜入射SV波和P波的疊加,構(gòu)造了基于地表實(shí)測地震動的斜入射體系,采用的思想如下:地震波在傳至地表時會發(fā)生波形轉(zhuǎn)換,生成反射SV波和反射P波,如圖1和圖2所示,并且入射角和反射角滿足Snell定律。地表地震動水平分量可以看做是由垂直入射的剪切波引起的,豎向分量是由垂直入射的膨脹波引起的。遵循這個思路,將地表地震動水平分量和豎直分量分別看作是由斜入射SV波和P波疊加引起的。

圖1 P波斜入射自由場Fig.1 A free-field system under obliquely incident P-wave

圖2 SV波斜入射自由場Fig.2 A free-field system under obliquely incident SV-wave

為了編程實(shí)現(xiàn)上述思路,設(shè)半無限空間自由面中心點(diǎn)O為坐標(biāo)原點(diǎn),如圖1、圖2所示。該處水平向?qū)崪y地震動為uh(t),在截取的部分地基左下角C點(diǎn)處SV波的入射時程為f(t),入射角度為θ0,該處P波的入射時程為g(t),入射角度為θ′0。它們疊加作用得到的地表中心點(diǎn)的地震動水平分量為uh(t-ti),其中ti表示疊加波傳播到地表O點(diǎn)的時間延遲,故入射的SV波和P波豎直分量之間存在以下關(guān)系。

(1)

式中:x,y分別代表地表自由面任意點(diǎn)水平向和豎直向坐標(biāo),角度上標(biāo)中帶有“′”符號的表示入射P波與其對應(yīng)的兩個反射角。f1和f2分別代表入射SV波反射生成的P波和SV波與其的波幅比;g1和g2分別代表入射P波反射生成的P波和SV波與入射P波的波幅比。根據(jù)Snell定律,P波和SV波的入射角度確定,f1、f2、g1和g2是可以計算得到的常數(shù)。故式(1)化簡為

(2)

為了使上式在地表任意點(diǎn)處均成立,入射SV波和P波必須同時達(dá)到自由面,且其入射角之間必須滿足如下等式

(3)

則入射P波可以根據(jù)入射SV波得到

(4)

其中

(5)

注意,根據(jù)式(4)可知入射SV波和P波并不是同時從C點(diǎn)同時入射的,為使疊加波場在地表產(chǎn)生的水平向地震動與實(shí)測地震動水平分量相等,需滿足

[(cosθ0+f1sinθ1-f2sinθ2)+

fv(sinθ′0+g1sinθ′1-g2cosθ′2)]

(6)

其中:

(7)

其中:

fsh=[(cosθ0+f1sinθ1-f2cosθ2)+

fv(sinθ′0+g1sinθ′1-g2cosθ′2)]-1

(8)

根據(jù)上述簡單推導(dǎo)[14-15],入射P波可由式(4)計算得到。然后根據(jù)波的傳播理論,斜入射SV波和P波在地表共同產(chǎn)生的波場就可以得到,其水平向地震動與地表實(shí)測地震動相等且具有非一致特性。

1.2 混凝土損傷模型

為了真實(shí)描述壩后式廠房的動力損傷過程,需采用合理的混凝土損傷本構(gòu)。由于ANSYS軟件中缺乏合理評價混凝土材料在地震工況下經(jīng)歷反復(fù)拉壓荷載后的損傷本構(gòu)。本文利用APDL語言編程,將混凝土單元的損傷通過彈性模量的折減來實(shí)現(xiàn)。應(yīng)力應(yīng)變關(guān)系如圖3所示,考慮了混凝土的拉伸損傷和其在地震動期間由于拉伸損傷產(chǎn)生單元受力性能下降而導(dǎo)致的壓剪損傷,分別采用最大拉應(yīng)變準(zhǔn)則和摩爾-庫倫準(zhǔn)則來判斷拉伸損傷和壓剪損傷。

圖3 混凝土應(yīng)力應(yīng)變曲線Fig.3 Stress-strain curves of concrete

2 壩后式廠房易損性分析

2.1 案例背景

我國西南某壩后式水電站地處兩個強(qiáng)震發(fā)生帶之間,其所屬壩后式廠房上游墻高38.1 m,下游廠房墻高23.6 m。根據(jù)已有的研究[16],壩后式廠房上游的大壩在地震動期間對其動力響應(yīng)有一定影響,為了更實(shí)際的反映壩后式廠房的動力響應(yīng),本文基于ANSYS軟件建立如圖4所示的廠壩二維模型。

圖4 有限元模型Fig.4 Finite element model

基礎(chǔ)長寬分別為1 000 m和500 m,材料參數(shù)取值如表1所示。

表1 材料參數(shù)表Table 1 Material parameters

經(jīng)有限元程序彈性分析可知,圖4中整個結(jié)構(gòu)的薄弱位置處于壩后式廠房的上下游墻部位,如圖5所示,在程序中對這部分結(jié)構(gòu)賦予2.1節(jié)中的混凝土損傷模型,以此來考察壩后式廠房的易損性。

圖5 上下游墻頂、墻底位置Fig.5 Location of the top and bottom of upstream and downstream walls

2.2 實(shí)測地震波選取與處理

太平洋地震工程研究中心(PEER)中的強(qiáng)震數(shù)據(jù)庫(NGA)記錄了全球淺層地殼地震中記錄的大量地震動,由于地震波的頻譜特性對結(jié)構(gòu)響應(yīng)差別影響很大,將水電工程水工建筑物抗震設(shè)計規(guī)范[17](NB35047—2015)中標(biāo)準(zhǔn)設(shè)計反應(yīng)譜輸入NGA數(shù)據(jù)庫中,再依據(jù)震級(6.5

表2 地震動記錄Table 2 List of earthquake records

由于實(shí)測地震動長持時會引起結(jié)構(gòu)剛度和強(qiáng)度退化,根據(jù)文獻(xiàn)[9]中使用的截斷方法,使用顯著持時將實(shí)測地震波截斷,只計算地震波能量由0.5%累計至95%所需的時間段。能量累計的計算公式根據(jù)式(9)計算。

(9)

式中:AI表示Arias[18]強(qiáng)度,tt表示地震波總持時,a(t)表示地震波加速度時程。這里采用MATLAB程序包對PEER中選取的26條實(shí)測地震波進(jìn)行截斷。為了進(jìn)行IDA分析,對選取的26條近場地震動進(jìn)行調(diào)幅,每條地震波加速度峰值按0.1g的級差調(diào)成0.1g至1.0g共10組。各條地震波的反應(yīng)譜如圖6所示,圖中圖例的數(shù)字代表的是表2中地震波的編號。

圖6 選取的26條地震波反應(yīng)譜Fig.6 Response spectra of 26 seismic waves

2.3 易損性指標(biāo)選取

結(jié)構(gòu)易損性分析中易損性指標(biāo)的選取至關(guān)重要。由于本文研究的對象是單層壩后式廠房。眾多研究者[1,3,19]對單層廠房的地震易損性研究均采用最大層間位移角作為廠房破壞等級劃分的參數(shù)。本文將壩后式廠房破壞等級按層間位移角劃分為輕微破壞、中等破壞、嚴(yán)重破壞三個級別。各自對應(yīng)的層間位移角限值[3]如表3所示。

表3 壩后式廠房結(jié)構(gòu)破壞等級與層間位移角的關(guān)系Table 3 Relationship between damage grade and story drift ratio of the power house at dam toe

2.4 易損性曲線分析

已有的易損性曲線通常有兩種繪制方法。第一種是直接擬合法,第二種是概率分布法。本文采用應(yīng)用范圍較廣的直接擬合法。為了對比傳統(tǒng)垂直入射模式與本文斜入射體系下基于地表實(shí)測地震動輸入模式對壩后式廠房的易損性曲線的影響,選取同樣的地震波對同一個有限元模型進(jìn)行增量動力學(xué)分析。經(jīng)過有限元程序計算后得到的易損性曲線如圖7所示。

圖7 上下游墻易損性曲線Fig.7 Vulnerability curves of the upstream and downstream walls

從圖7可見,斜入射體系下基于地表地震動入射的壩后式廠房上下游墻的易損性曲線和傳統(tǒng)的垂直入射得到的易損性曲線有明顯不同?？傮w上來看,各個破壞等級下的易損性曲線形狀類似,均是在PGA小于0.5g時超越概率增長迅速,PGA超過0.5g時超越概率增長緩慢直至100%。壩后式廠房上下游墻體在垂直入射模式下的破壞概率比斜入射體系下高。根據(jù)中國地震局地質(zhì)研究所《金沙江龍開口水電站工程場地地震安全性評價和水庫誘發(fā)地震評價報告》成果,本文的壩后式廠房所在場地50年超越概率5%的設(shè)計加速度為0.24g,100年超越概率2%設(shè)計加速度為0.394g。在三種破壞模式下將這兩條加速度直線與相應(yīng)的易損性曲線相交。在輕微破壞等級下,垂直入射時,50年超越概率5%的設(shè)計加速度時壩后式廠房上游墻的超越概率是0.89,斜入射體系下基于地表實(shí)測地震動入射時這一值為0.63;100年超越概率2%的設(shè)計加速度下,兩種入射體系的輕微破壞概率分別是1和0.84。下游墻的輕微破壞概率比上游墻大,50年超越概率5%的設(shè)計加速度時分別為0.94和0.82;100年超越概率2%的設(shè)計加速度下輕微破壞概率分別是1和0.93,這一峰值加速度下兩種入射體系計算得到的差別最小,只有7%。在中等破壞等級時,垂直入射下,50年超越概率5%的設(shè)計加速度時壩后式廠房上游墻的超越概率是0.71,斜入射體系下基于地表實(shí)測地震動入射時這一值為0.53;100年超越概率2%的設(shè)計加速度下,兩種入射體系的輕微破壞概率分別是0.91和0.65,這一峰值加速度下兩種入射體系計算得到的差別最大,達(dá)到26%。下游墻的輕微破壞概率比上游墻大,50年超越概率5%的設(shè)計加速度時分別為0.81和0.59;100年超越概率2%的設(shè)計加速度下輕微破壞概率分別是0.96和0.73。在嚴(yán)重破壞等級下,垂直入射時,50年超越概率5%的設(shè)計加速度時壩后式廠房上游墻的超越概率是0.61,斜入射體系下基于地表實(shí)測地震動入射時這一值為0.44;100年超越概率2%的設(shè)計加速度下,兩種入射體系的輕微破壞概率分別是0.76和0.58。下游墻的輕微破壞概率比上游墻大,50年超越概率5%的設(shè)計加速度時分別為0.69和0.53;100年超越概率2%的設(shè)計加速度下輕微破壞概率分別是0.88和0.68。從地震波斜入射體系下基于地表實(shí)測地震動的壩后式廠房上下游墻的易損性曲線和傳統(tǒng)的垂直入射得到的結(jié)果對比可以看出,斜入射作用對易損性曲線的影響不可忽略,需在同類壩后式廠房的地震安全性評價中考慮地震波斜入射作用的影響。

3 結(jié)論

本文編寫的程序克服了斜入射地震波輸入難以按地表實(shí)測地震動參數(shù)為目標(biāo)來調(diào)幅進(jìn)行結(jié)構(gòu)易損性分析的困難,并在程序中嵌入了混凝土損傷模型。以我國西南地區(qū)某壩后式廠房為工程實(shí)例進(jìn)行了結(jié)構(gòu)易損性分析。經(jīng)過520組數(shù)值實(shí)驗(yàn),與傳統(tǒng)垂直入射模式下的易損性曲線做對比后得到了以下結(jié)論:

(1)根據(jù)波的傳播理論,推導(dǎo)了合成地表水平向地震動的斜入射SV波和P波公式表達(dá),使得地表水平向地震動與實(shí)測地震動相等且具有非一致特性。

(2)傳統(tǒng)的垂直入射方式得到的壩后式廠房上下游墻的易損性曲線的各級破壞超越概率比地震波斜入射體系下基于地表實(shí)測地震動模式得到的易損性曲線的超越概率大。最大時相差26%,故進(jìn)行壩后式廠房的易損性分析時需考慮地震波斜入射的影響。

(3)壩后式廠房下游墻各個等級的破壞概率比上游墻大。在地震波斜入射體系下,下游墻體各級損壞超越概率中最高比上游墻體超出19%,對同類壩后式廠房進(jìn)行抗震設(shè)計時應(yīng)當(dāng)予以關(guān)注。

(4)該方法僅適用于線彈性均質(zhì)地基。在實(shí)際工程中,地基往往是非均質(zhì)的,并且非巖質(zhì)地基在強(qiáng)震作用下具有非線性特性。因此,如何將該方法推廣到非均質(zhì)地基和非線性地基是今后研究的重點(diǎn)。