鄒 炎， 景立平， 崔 杰， 李永強

(1.中國地震局工程力學(xué)研究所，黑龍江 哈爾濱 150080； 2.廣州大學(xué)土木工程學(xué)院，廣東 廣州 510006)

飽和砂土中箱型隧道的地震反應(yīng)分析①

鄒 炎1， 景立平1， 崔 杰2， 李永強1

(1.中國地震局工程力學(xué)研究所，黑龍江 哈爾濱 150080； 2.廣州大學(xué)土木工程學(xué)院，廣東 廣州 510006)

飽和砂土地基在地震作用下存在液化的潛在危險，液化引起的地基失效可能會導(dǎo)致地下結(jié)構(gòu)的嚴(yán)重震害。以O(shè)pensees作為計算平臺，對飽和砂土中帶中柱箱型隧道的地震反應(yīng)進(jìn)行輸入不同幅值地震動時的動力數(shù)值計算，研究場地和結(jié)構(gòu)的加速度反應(yīng)及其頻譜特性、場地的永久變形、隧道的震后位移以及隧道的內(nèi)力分布。計算結(jié)果表明，飽和砂土中箱型隧道的地震附加內(nèi)力仍受周圍土體的相對位移控制，此外在震后隧道可能會產(chǎn)生側(cè)移和上浮的永久位移，并且可能存在殘余內(nèi)力。

砂土液化； 地下結(jié)構(gòu)； 隧道； 數(shù)值模擬； Opensees

0 引言

近江或近海區(qū)域的上覆土層往往由含水率很高的沖積砂層組成，這種近飽和的砂土地基在地震作用下存在液化的潛在危險，液化引起的地基失效可能會導(dǎo)致地下結(jié)構(gòu)的嚴(yán)重震害。1989年美國加利福尼亞州的Loma Prieta地震中，地基液化引起部分隧道的上浮，造成隧道的破壞[1]；1995年日本阪神地震中，神戶市海岸一帶的可液化填筑砂礫層存在液化現(xiàn)象，部分地鐵區(qū)間隧道和地下管道的震害可能與此有關(guān)[2]。因此，對可液化土層中地下結(jié)構(gòu)的力學(xué)性能和震害機理進(jìn)行深入研究是十分必要的。

以前在實際工程中更多的依賴場址初勘階段的現(xiàn)場液化判別，通過剪切波速試驗和標(biāo)準(zhǔn)貫入試驗等現(xiàn)場試驗進(jìn)行經(jīng)驗性簡化判別[3]。然而在水底隧道等地下結(jié)構(gòu)工程中很難進(jìn)行現(xiàn)場液化判別試驗，因此需要通過其他途徑進(jìn)行場地液化影響的研究。隨著砂土液化本構(gòu)的發(fā)展和完善，近年來數(shù)值模擬已經(jīng)成為研究可液化場地地震反應(yīng)的一種重要手段。劉華北和宋二祥等[4]在動力兩相體非線性有限元軟件Dyna-Swandyne-II中采用Pastor-Zienkiewicz III廣義塑性模型模擬可液化土的動力特性，研究了地下結(jié)構(gòu)在地震液化作用下的響應(yīng)。莊海洋等[5]在在ABAQUS商用軟件的計算平臺上，建立了可液化土地基-地鐵車站結(jié)構(gòu)非線性靜動力耦合相互作用分析模型，對埋有雙層三跨島式大型地鐵車站的周圍場地進(jìn)行了地震反應(yīng)規(guī)律數(shù)值模擬分析。

本文以O(shè)pensees作為計算平臺,對飽和砂土中帶中柱箱型隧道的地震反應(yīng)進(jìn)行輸入不同幅值地震動時的動力數(shù)值計算，研究場地和結(jié)構(gòu)的加速度反應(yīng)及其頻譜特性、場地的永久變形、隧道的變形和位移、以及隧道的內(nèi)力分布，分析可液化土層中箱型隧道的地震反應(yīng)規(guī)律和震害機理，對可液化土層中地下結(jié)構(gòu)地震反應(yīng)的研究具有重要意義。

1 有限元模型

1.1 砂土液化本構(gòu)模型

Opensees材料庫中的Pressure Dependent Multi Yield Material基于YANG Zhao-hui提出的砂土液化本構(gòu)模型[6]，可以反映飽和砂土在循環(huán)荷載作用下的變形行為和孔壓增長特性。根據(jù)Prevost提出的J2型迭套屈服面模型[7]，YANG Zhao-hui定義飽和砂土的屈服面為：

(1)

圖1 主應(yīng)力空間和偏應(yīng)力平面的屈服面示意圖Fig.1 Conical yield surfaces in principal stress space and deviatoric plane

圖2 循環(huán)加載過程的階段劃分Fig.2 Stages of the cyclic loading process

圖3 計算模型及觀測點示意圖Fig.3 The calculation model and observation points

1.2 計算模型的建立

在Opensees中建立如圖3所示的二維平面應(yīng)變模型，計算范圍水平方向選取40 m，豎直方向埋深0～1 m之間是人工填土，1～20 m是可液化的飽和砂土，20 m更深為基巖。區(qū)間隧道為尺寸10 m×6 m的帶中柱箱型樣式，頂面埋深為5 m，中柱截面為0.4 m×0.6 m的矩形，因為在垂直于計算平面方模型厚度取1 m，頂板、側(cè)墻和底板的截面為0.6 m×1 m的矩形。人工填土和飽和砂土的的本構(gòu)分別采用Opensees中的PIMY材料和PDMY材料，相關(guān)參數(shù)[10]分別見表1和表2；隧道材料為C30混凝土，采用彈性本構(gòu)，參數(shù)見表3。取圖中的兩個豎直觀測面A和B，分別位于土層中和結(jié)構(gòu)所在位置，每個觀測面各有6個觀測點，其中B3、B4、B5位于結(jié)構(gòu)上。土層采用網(wǎng)格尺寸為1 m×1 m的quad單元劃分，隧道采用尺寸為1 m的elastic Beam Column單元劃分，側(cè)邊界采用自由場邊界，底邊界采用加速度邊界，直接施加地震動的加速度時程。輸入地震波采用0.05 g、0.1 g、0.2 g和0.4 g的El-Centro地震波，其中0.1 g的加速度時程和加速度β譜如圖4所示。

表1 人工填土本構(gòu)參數(shù)

圖4 輸入0.1g的El-Centro波Fig.4 Input El-centro wave of 0.1 g

2 計算結(jié)果分析

2.1 初始靜力分析

在進(jìn)行動力分析之前，首先要對整個計算模型施加重力荷載，并進(jìn)行初始靜力分析，最終得到隧道各關(guān)鍵部位的初始內(nèi)力如表4。

表2 飽和砂土本構(gòu)參數(shù)

表3 隧道材料參數(shù)

表4 隧道各部位的初始內(nèi)力

2.2 超靜孔壓分析

在輸入0.05 g的El-Centro地震波時，A3、A4、A5和A6觀測點的超靜孔壓比時程如圖5所示，其他幅值下的超靜孔壓比時程不再列出。綜合比較超靜孔壓比的值見表5，可知：在0.1 g時，A3點和A5點率先達(dá)到完全液化；在0.2 g時，A4點也接近液化；而到0.4 g時A4點和A6點達(dá)到完全液化。說明箱型隧道頂板和底板所在深度的A3點和A5點受到隧道運動的影響較大，更容易達(dá)到完全液化，而隧道中心所在深度的A4點反而不那么容易完全液化。

圖5 0.05 g時超靜孔壓比時程圖Fig.5 Hyperstatic pore pressure ratio time history at 0.05 g

表5 不同輸入地震動幅值下各觀測點的超靜孔壓比

2.3 加速度放大系數(shù)分析

各幅值輸入波下A、B兩個觀測面的加速度放大系數(shù)沿深度的分布如圖6所示，其中兩個橫線之間是箱型隧道所在高度?？芍?，在小震(0.05 g、0.1 g)下，從基巖到地表加速度峰值逐漸增大；在大震下(0.2 g、0.4 g)，從基巖到地表加速度峰值從埋深15 m處開始逐漸減小，直至埋深4 m處才開始增大。這是因為在小震下飽和砂土還沒有液化，場地對地震動具有放大作用，導(dǎo)致了加速度峰值的增大；在大震下飽和砂土產(chǎn)生液化，埋深15 m至4 m處為液化區(qū)，該區(qū)域內(nèi)的砂土完全液化后喪失了剪切模量，損失了傳播剪切波的能力，導(dǎo)致地震波能量大量損耗，加速度峰值反而減小。另一方面，隨著輸入地震動峰值的逐漸增大，箱型隧道中柱的加速度峰值分布與相同埋深位置A觀測面砂土的加速度峰值分布的差別也逐漸增大，這是因為飽和砂土液化后，作為地基土其剪切模量大大衰減，剛度也大大降低，出現(xiàn)地基失效，整體對隧道的約束作用大大減弱，隧道隨周圍土體運動的同步性也大大減小，而呈現(xiàn)更強的隨機性。

圖6 不同幅值地震動時觀測面A和B的加速度放大系數(shù)分布Fig.6 Acceleration amplification factors distribution of surface A and B under earthquake motion with different amplitudes

2.4 地面反應(yīng)譜分析

對不同輸入地震動峰值下地面觀測點B1的加速度反應(yīng)β譜進(jìn)行比較(圖7)?？芍狠斎?.05 g、0.1 g和0.2 g的El-Centro地震波時，地面觀測點B1的加速度反應(yīng)β譜基本相同；而輸入0.4 g的El-Centro地震波時，地面觀測點B1的加速度反應(yīng)β譜發(fā)生了較大變化，存在兩個很明顯的峰值，第一個頻率大致在0.1 Hz處。第二個大致在0.6 Hz處。這是因為埋深4～15 m的液化區(qū)內(nèi)的飽和砂土已經(jīng)完全液化，剪切模量大大減小，而埋深15～20 m的非液化區(qū)內(nèi)的飽和砂土沒有液化，剪切模量變化不大，因此地基土出現(xiàn)了明顯的分層現(xiàn)象，并且液化區(qū)內(nèi)飽和砂土剪切模量的減小反映0.6 Hz處反應(yīng)譜峰值的衰減。

圖7 不同幅值地震動時地面觀測點B1的加速度反應(yīng)β譜Fig.7 Acceleration response spectrum of point B1 under earthquake motion with different amplitudes

圖8 不同幅值地震動時地面的震后豎向位移分布Fig.8 Post-earthquake vertical displacements distribution of ground under earthquake motion with different amplitudes

2.5 地面震后豎向位移分析

不同輸入地震動峰值下計算區(qū)域內(nèi)水平地面所有點的震后豎向位移分布如圖8所示?？芍簩τ陲柡蜕巴量梢夯瘓龅?，地震后隧道正上方地面會上浮，而隧道兩側(cè)上方地面會下沉，并且隨著輸入地震動峰值的增大和液化區(qū)的擴展，以及飽和砂土液化程度的增大，地面的上浮量和下沉量也會增大。

2.6 隧道震后位移分析

對箱型隧道上各點的震后位移(包括水平和豎向)進(jìn)行監(jiān)測，得到震后隧道的震后位移如圖9所示，其中實線為隧道的初始位置，虛線為隧道的震后位置?？芍涸谛≌鹣?0.05 g和0.1 g)，隧道的震后位置與初始位置基本重合，隧道的震后位移很??；在中大震下(0.2 g和0.4 g)，隧道在震后出現(xiàn)了明顯的側(cè)移和上浮。其中隧道的中柱頂端即觀測點B3的震后位移與輸入地震動峰值的關(guān)系如圖10所示。可知：隨著輸入地震動峰值的增大，隧道的側(cè)移量和上浮量也會增大，并且側(cè)移量明顯增加的更快。

圖9 不同幅值地震動時隧道的震后位移Fig.9 Post-earthquake displacement of tunnel under earthquake motion with different amplitudes

圖10 B3點的震后位移vs地震動幅值Fig.10 Post-earthquake displacement vs amplitude of earthquake motion at point B3

2.7 隧道動內(nèi)力峰值分析

表6列出不同輸入地震動幅值下隧道各關(guān)鍵構(gòu)件連接部位的最大動內(nèi)力幅值。可知隧道的動軸力在底板側(cè)端最大，動剪力在側(cè)墻墻底最大，動彎矩在底板側(cè)端和側(cè)墻墻底最大。隧道的動內(nèi)力峰值與隧道頂?shù)装遄畲笙鄬ξ灰频年P(guān)系如圖11所示?？芍S著隧道最大相對位移的增大，隧道的動內(nèi)力峰值也逐漸增大，并且兩者之間呈近似線性關(guān)系，這說明即使對于可液化的飽和砂土場地，場地中的地下結(jié)構(gòu)仍然滿足內(nèi)力受變形尤其是周圍土體相對位移控制的規(guī)律。

圖11 隧道的動內(nèi)力峰值vs.頂?shù)装遄畲笙鄬ξ灰艶ig.11 Peak dynamic internal force of tunnel vs maximum relative displacement of roof and floor

2.8 隧道殘余動內(nèi)力分析

輸入0.2 g El-Centro波時隧道左墻墻頂處的動剪力時程如圖12所示，可以看到該處在震后仍然存在較大的殘余動剪力，這說明可液化場地中的地下結(jié)構(gòu)在震后可能存在較大的殘余動內(nèi)力。列出不同輸入地震動幅值下隧道各關(guān)鍵構(gòu)件連接部位的殘余動內(nèi)力于表7?？芍湫退淼罋堄鄤觾?nèi)力最大的位置與其動內(nèi)力峰值最大的位置基本相同，但隧道中柱的殘余彎矩與殘余剪力基本為零。隧道的殘余動內(nèi)力與輸入地震動幅值的關(guān)系如圖13所示?？芍簭?.05 g到0.1 g，隧道的殘余動內(nèi)力大幅增加；而從0.1 g到0.2 g再到0.4 g，則增幅變緩。這是因為隧道的殘余動內(nèi)力是由隧道周圍的飽和砂土在液化過程中超靜孔壓的增加引起的，從0.05 g到0.1 g，隧道周圍的飽和砂土從基本無液化到接近完全液化，其超靜孔壓增長較大，因此隧道的殘余動內(nèi)力也迅速增大；0.1 g以后隧道周圍的飽和砂土已經(jīng)基本液化，其超靜孔壓幾乎不再增長，因此殘余動內(nèi)力也幾乎不再增大。

表6 不同輸入地震動幅值下隧道各關(guān)鍵連接部位的動內(nèi)力峰值

表7 不同輸入地震動幅值下隧道各關(guān)鍵連接部位的殘余動內(nèi)力

圖12 0.2 g時隧道左墻墻頂處的動剪力時程Fig.12 Dynamic shear force of left wall top at 0.2 g

圖13 隧道的殘余動內(nèi)力vs地震動幅值Fig.13 Remnant dynamic internal forces vs amplitudes

3 結(jié)論

本文以O(shè)pensees作為計算平臺,對飽和砂土中帶中柱箱型隧道的地震反應(yīng)進(jìn)行了輸入不同幅值地震動時的動力數(shù)值計算，得到如下結(jié)論：

(1) 隧道頂板和底板所在深度的飽和砂土比隧道中心所在深度的飽和砂土更容易液化；

(2) 飽和砂土在小震下對地震波具有放大作用，在大震下完全液化后反而具有隔震效應(yīng)；

(3) 震后隧道正上方地面會上浮，而隧道兩側(cè)上方地面會下沉；

(4) 隧道在震后可能會出現(xiàn)永久的側(cè)移和上?。?/p>

(5) 飽和砂土中的箱型隧道的地震附加內(nèi)力仍然受周圍土體的相對位移控制；

(6) 震后隧道可能存在超靜孔隙水壓力引起的殘余內(nèi)力。

References)

[1] Schmidt B，Hashash Y， Stimac T.US Immersed Tube Retrofit[J].Tunnels Tunneling International Magazine，1998，30(11)：22-24.

[2] 劉惠珊.1995年阪神大地震的液化特點[J].工程抗震，2001(1)：21-26.LIU Hui-shan.Some Features of Liquefaction During the 1995 Great Hanshin-awaji Eartgqyake[J].Earthquake Resistant Engineering，2001(1)：21-26.(in Chinese)

[3] 中華人民共和國國家標(biāo)準(zhǔn)編寫組.GB50011-2001，建筑抗震設(shè)計規(guī)范[S].北京：中國建筑工業(yè)出版社，2001.The National Standards Complication Group of People’s Republic of China.GB5001 1-2001，Code for Seismic Design of Building[S].Beijing：China Architecture & Building Press，2001.

[4] 劉華北，宋二祥.可液化土中地鐵結(jié)構(gòu)的地震響應(yīng)[J].巖土力學(xué)，2005，26(3)：381-386，391.LIU Hua-bei，SONG Er-xiang.Earthquake Induced Liquefaction Response of Subway Structure in Liquefiable Soil[J].Rock and SoilL Mechanics，2005，26(3)：381-386，391.(in Chinese)

[5] 莊海洋，龍慧，陳國興，等.可液化地基中地鐵車站周圍場地地震反應(yīng)分析[J].巖土工程學(xué)報，2012，34(1)：83-88.ZHUANG Hai-yang，LONG Hui，CHEN Guo-xing，et al.Seismic Responses of Surrounding Site of Subway Station in Liquefiable Foundation[J].Chinese Journal of Geotechnical Engineering，2012，34(1)：83-88.(in Chinese)

[6] Yang Z，Elgamal A，Parra E.Computational Model for Cyclic Mobility and Associated Shear Deformation[J].Journal of Geotechnical and Geoenvironmental Engineering，2003，129(12)：1119-1127.

[7] Prevost J H.A Simple Plasticity Theory for Frictional Cohesionless Soils[J].International Journal of Soil Dynamics and Earthquake Engineering，1985，4(1)：9-17.

[8] Ishihare K，Tatsuoka Fumio，Yasuda S.Undrained Deformation and Liquefaction of Sand Under Cyclic Stresses[J].Soils and Foundations，1975，15(1)：29-44.

[9] Mroz Z.On the Description of Anisotropic Workhardening[J].Journal of the Mechanics and Physics of Solids，1967，15(3)：163-175.

[10] Silvia Mazzoni，F(xiàn)rank McKenna，Michael H Scott，et al.OpenSees Command Language Manual[EB/OL].http://opensees.Berkeley.Edu/OpenSees/manuals/usermanual/，2007-07-01.

Seismic Response Analysis of Box Tunnels in Saturated Sand

ZOU Yan1， JING Li-ping1， CUI Jie2， LI Yong-qiang1

(1.InstituteofEngineeringMechanics，ChinaEarthquakeAdministration，Harbin，Heilongjiang150080，China；2.SchoolofCivilEngineering，GuangzhouUniversity，Guangzhou，Guangdong510006，China)

Foundations in saturated sand have the potential to liquefy under the effects of earthquakes，and foundation failure caused by liquefaction can lead to a serious damage of underground structures.Considering Opensees as the computing platform，a dynamic numerical simulation was performed to calculate the seismic response of a box tunnel with mid-columns in saturated sand when seismic waves of different amplitudes were input.Acceleration responses and the spectral properties of ground and structure，as well as the permanent deformation of the ground，aftershock displacement of the tunnel，and internal force distribution of the tunnel，were studied.The results show that the additional internal seismic forces in a box tunnel in saturated sand are controlled by the relative displacement of the surrounding soil.The permanent displacement of the tunnel such as side shift and floating probably appear after the main shock；in addition，remnant internal forces also probably occur.

liquefaction of sand soil； underground structure； tunnel； numerical simulation； Opensees

2014-08-20

國家自然科學(xué)基金項目(51408566)；中央級公益性研究所基本科研業(yè)務(wù)費專項(2014B03)；國家自然科學(xué)基金項目(51438004)；博士點基金(20124410110004)；973課題(2011CB013606)

鄒 炎(1987- )，男，博士研究生，主要從事地下結(jié)構(gòu)抗震方面的研究.E-mail：latro@qq.com

TU43

A

1000-0844(2015)02-0329-07

10.3969/j.issn.1000-0844.2015.02.0329