李文倩,佟大威,王　振,朱曉斌,魯文妍,2

(1. 天津大學(xué)水利工程仿真與安全國家重點(diǎn)實(shí)驗(yàn)室,天津　300072;2. 南京水利科學(xué)研究院材料結(jié)構(gòu)研究所,江蘇 南京　210029)

大型地下洞室群相鄰塊體地震響應(yīng)

李文倩1,佟大威1,王振1,朱曉斌1,魯文妍1,2

(1. 天津大學(xué)水利工程仿真與安全國家重點(diǎn)實(shí)驗(yàn)室,天津300072;2. 南京水利科學(xué)研究院材料結(jié)構(gòu)研究所,江蘇 南京210029)

摘要：為了更好地模擬地下洞室中存在的相鄰塊體結(jié)構(gòu)在地震荷載作用下的運(yùn)動(dòng)與動(dòng)態(tài)響應(yīng)特性,建立了一種考慮塊體間黏結(jié)強(qiáng)度影響的地下洞室相鄰塊體動(dòng)力抗滑模型,并通過數(shù)值模型分析了該抗滑模型的有效性與優(yōu)越性。在工程實(shí)例分析中,采用經(jīng)典抗滑模型時(shí)Ⅰ類相鄰塊體中兩塊體水平X向位移在整個(gè)地震階段的均方根值分別為2.16 cm和2.02 cm,采用修正抗滑模型時(shí)兩塊體對(duì)應(yīng)響應(yīng)值分別為2.00 cm和1.94 cm;Ⅱ類相鄰塊體中兩塊體在經(jīng)典抗滑模型與修正抗滑模型條件下的對(duì)應(yīng)響應(yīng)值分別為3.68 cm和3.65 cm與3.55 cm和3.54 cm,上述結(jié)果反映了相鄰塊體在黏結(jié)特性的作用下具有更好的穩(wěn)定性,說明該模型具有較好的工程適用性與應(yīng)用價(jià)值。

關(guān)鍵詞：地下洞室;相鄰塊體;黏結(jié)特性;抗滑模型;地震響應(yīng)

大型水利工程地下洞室往往處于極端復(fù)雜的地質(zhì)條件下,洞室周邊分布的塊體中既包括單一塊體,也包括相鄰塊體,其中相鄰塊體結(jié)構(gòu)在地震荷載激勵(lì)下產(chǎn)生的動(dòng)力響應(yīng)特性是判斷洞室穩(wěn)定性的主要安全參量,也是評(píng)估水利工程結(jié)構(gòu)整體安全性的重要指標(biāo)。

針對(duì)大型地下洞室塊體結(jié)構(gòu)的動(dòng)力特性模擬,不僅涉及結(jié)構(gòu)所受地震荷載的激勵(lì)特性,同時(shí)也需要考慮來自塊體自身的動(dòng)態(tài)穩(wěn)定機(jī)理及其接觸抗滑特性的影響。最早的巖體工程中塊體與動(dòng)接觸理論由Goodman等[1-2]提出并得以廣泛應(yīng)用,石根華根據(jù)塊體穩(wěn)定分析的幾何原理與巖體的影響因素建立了滑動(dòng)體的數(shù)學(xué)模型并提出塊體穩(wěn)定分析方法。至今,考慮地下洞室塊體結(jié)構(gòu)動(dòng)力特性的方法主要包括連續(xù)數(shù)值方法[4]、非連續(xù)數(shù)值方法[5-6]和耦合數(shù)值方法[7-10],這些方法多局限于單一塊體的研究而對(duì)相鄰塊體在強(qiáng)震作用下的動(dòng)態(tài)響應(yīng)研究較少,同時(shí)現(xiàn)有理論與研究中往往只重視塊體本身的動(dòng)力性質(zhì)而忽視巖體間黏結(jié)強(qiáng)度對(duì)接觸面動(dòng)態(tài)特性的影響效應(yīng),在解決實(shí)際工程問題上具有一定的局限性。

本文在研究相鄰塊體地震響應(yīng)的基礎(chǔ)上,考慮實(shí)際工程中裂隙結(jié)構(gòu)面的不完整性與黏結(jié)性,建立了地下洞室相鄰塊體周邊接觸面動(dòng)力抗滑模型,在理論上給出了相鄰塊體間接觸面抗滑模型的基本原理與適用條件,同時(shí)分別從數(shù)值模型與實(shí)際工程角度驗(yàn)證說明了該抗滑模型在解決相鄰塊體結(jié)構(gòu)抗滑與分析接觸面動(dòng)力特性上的適用性與優(yōu)越性。建立的抗滑模型可以更好地模擬實(shí)際工程中塊體周邊接觸面上的黏結(jié)與抗滑性質(zhì),并將單一塊體動(dòng)力問題擴(kuò)展到相鄰塊體上,更好地模擬局部塊體在強(qiáng)震作用下的動(dòng)態(tài)響應(yīng)特性。

1相鄰塊體接觸面抗滑模型

1.1塊體與周邊巖體間接觸面抗滑模型

地下洞室塊體在滑動(dòng)過程中受到的抗滑力既包括切向摩擦力,也包括黏聚力[11-12]。結(jié)合廣義摩爾庫倫模型，可將塊體與基巖結(jié)構(gòu)間接觸面上的切向力描述為以下形式[13-14]:

(1)

式中:τct——塊體與巖體接觸面上切向應(yīng)力,Pa;τcn——塊體與巖體接觸面上法向應(yīng)力,Pa;kt——接觸面切向剛度,Pa/m;ut——接觸面上切向相對(duì)變形,m;f——接觸面上的摩擦系數(shù);c——塊體與巖體結(jié)構(gòu)接觸面上的黏聚力,Pa;sgn( )——符號(hào)函數(shù),給定接觸面抗滑力的符號(hào)與方向。

1.2相鄰塊體間接觸面抗滑模型

地下洞室相鄰塊體的特點(diǎn)是塊體間有一個(gè)或多個(gè)面處于接觸狀態(tài),塊體在外力作用下除了與巖體結(jié)構(gòu)有接觸作用外,塊體間的接觸面上也有相對(duì)變形和接觸作用。在相鄰雙塊體模型中,假定相鄰雙塊體間接觸面為平面,按照含有臨空面塊體的數(shù)量,分為2種情況：當(dāng)相鄰塊體間的非接觸面中每個(gè)塊體均存在地下洞室的洞室開挖面,此類型相鄰塊體在本文中簡稱“Ⅰ類相鄰塊體”, 模型簡化的典型情況見圖1;當(dāng)相鄰塊體間的非接觸面中有且只有一個(gè)塊體存在地下洞室的洞室開挖面,在本文中簡稱“Ⅱ類相鄰塊體”,模型簡化的典型情況見圖2。

圖1?、耦愊噜弶K體示意圖Fig. 1　Adjacent blocks of type Ⅰ

圖2?、蝾愊噜弶K體示意圖Fig. 2　Adjacent blocks of type Ⅱ

1.2.1Ⅰ類相鄰塊體

(2)

1.2.2Ⅱ類相鄰塊體

圖2中,同樣假設(shè)2個(gè)塊體間的接觸面為P,塊體1簡稱為“內(nèi)側(cè)塊體”,其在P上的變形為u″t1,塊體2簡稱為“外側(cè)塊體”,其在P上的變形為u″t2,接觸界面間的切向相對(duì)變形為Δu″t,法向變形為u″n,切向應(yīng)力為τ″ct,切向剛度為k″t,黏聚力為c″,相應(yīng)塊體間接觸面上的抗滑模型動(dòng)力特性可以按以下公式進(jìn)行描述:

(3)

2數(shù)值模型分析

2.1模型建立

圖3　數(shù)值模型Fig. 3　Numerical model

為研究理論模型的適用性,建立了圖3所示的數(shù)值模型。其中塊體1尺寸為20 m×10 m×10 m(長×寬×高),塊體2尺寸為60 m×30 m×10 m,下部基巖結(jié)構(gòu)尺寸為120 m×60 m×10 m。各結(jié)構(gòu)賦予相同材料屬性,密度為2 500 kg/m3,彈性模量為30 GPa,泊松比為0.20,對(duì)下部基巖結(jié)構(gòu)底面采用全約束,側(cè)面采用法向約束。在加載時(shí)首先對(duì)塊體1表面施加一個(gè)豎直方向且大小為250 kPa的壓力,平衡后再在水平X向?qū)?個(gè)塊體側(cè)面施加一個(gè)時(shí)長為8 s的動(dòng)態(tài)推力,且推力最大值為200 kPa,水平荷載歸一化時(shí)程線如圖4所示。

圖4　水平荷載歸一化時(shí)程線Fig. 4　Normalized time-history curve of horizontal loads

2.2模型結(jié)果分析

圖5給出2種抗滑模型條件下,上、下塊體的加速度響應(yīng)情況。從圖5中可以看出,采用經(jīng)典抗滑模型(經(jīng)典庫倫模型)與修正抗滑模型時(shí),下部塊體滑動(dòng)分別出現(xiàn)在總加載時(shí)間的第2.7 s和2.9 s時(shí),上部塊體的最終穩(wěn)定加速度分別為2.0 m/s2和1.33 m/s2，下部塊體的最終穩(wěn)定加速度分別為1.80 m/s2和1.25 m/s2,后者動(dòng)態(tài)響應(yīng)均小于前者計(jì)算結(jié)果,說明修正抗滑模型對(duì)于塊體結(jié)構(gòu)運(yùn)動(dòng)狀態(tài)影響的合理性。

圖6給出了2種抗滑模型條件下2塊塊體間及下部塊體與基巖間接觸面上剪應(yīng)力曲線。上部塊體相對(duì)于下部塊體滑動(dòng)后的最終恒定剪應(yīng)力分別為40.0 kPa和44.6 kPa。下部塊體與巖體間接觸面上,塊體滑動(dòng)后出現(xiàn)的最大剪應(yīng)力分別為25.9 kPa和30.6 kPa,結(jié)果體現(xiàn)了修正模型中黏聚力的重要作用。

3工 程 實(shí) 例

3.1工程概況與模型建立

工程實(shí)例位于瀾滄江下游,采用NURBS-TIN-BRep混合數(shù)據(jù)結(jié)構(gòu)與幾何建模技術(shù)建立工程中的復(fù)雜地形體與地下結(jié)構(gòu),并基于塊體理論耦合隨機(jī)結(jié)構(gòu)面網(wǎng)絡(luò)模擬技術(shù)快速建立地下洞室區(qū)域內(nèi)的曲面塊體[6]。選取2類相鄰塊體為研究對(duì)象,建立的地形、洞室及塊體模型分別見圖7和圖8,模型尺寸為837 m×446 m×(299～551)m(長×寬×高)。地下洞室?guī)r體以Drucker-Prager強(qiáng)度準(zhǔn)則作為屈服準(zhǔn)則,巖體的彈性模量為16 GPa,密度為2.71 g/cm3,泊松比為0.23,抗拉強(qiáng)度為1.6 MPa,內(nèi)摩擦角為52°,采用瑞利阻尼模擬。地震激勵(lì)采用唐山余震天津醫(yī)院地震波,輸入方向?yàn)樗絏向,地震波峰值加速度為1.46 m/s2,加載時(shí)間為19.20 s,根據(jù)文獻(xiàn)[15]對(duì)地表地震波進(jìn)行處理后,獲得基巖處的地震波加速度時(shí)程曲線。

圖5　塊體滑動(dòng)加速度曲線Fig. 5　Curves of sliding acceleration of blocks

圖6　接觸面剪應(yīng)力曲線Fig. 6　Curves of shear stress on contact surface

圖7　Ⅰ類相鄰塊體與地基位置的注釋圖Fig. 7　Position of foundation and adjacent blocks of type Ⅰ

圖8?、蝾愊噜弶K體與地基位置的注釋圖Fig. 8　Position of foundation and adjacent blocks of type Ⅱ

3.2Ⅰ類相鄰塊體響應(yīng)分析

圖9給出2種抗滑模型條件下相鄰兩塊體水平X向相對(duì)位移對(duì)比。采用經(jīng)典模型與修正模型時(shí)左側(cè)與右側(cè)塊體水平X向最大位移分別為5.42 cm和4.87 cm與5.06 cm和4.77 cm,地震階段的均方根值分別為2.16 cm和2.02 cm與2.00 cm和1.94 cm。修正模型條件下的位移在數(shù)值上相對(duì)小于經(jīng)典模型計(jì)算結(jié)果,說明在考慮塊體黏結(jié)效應(yīng)的抗滑模型在Ⅰ類相鄰兩塊體應(yīng)用中具有較好適用性。

圖9?、耦愊噜弶K體水平X向相對(duì)位移Fig. 9　Relative displacement of adjacent blocks of type Ⅰ in horizontal X direction

圖10給出了在2種模型條件下兩塊體與基巖接觸面上代表點(diǎn)的剪切應(yīng)力對(duì)比。基于修正抗滑模型計(jì)算獲得的剪切應(yīng)力在數(shù)值上大于經(jīng)典模型結(jié)果,修正模型條件下兩塊體最大剪切應(yīng)力值分別為0.80 MPa與0.58 MPa,經(jīng)典模型為0.60 MPa與0.56 MPa,平均剪切應(yīng)力分別為0.68 MPa與0.44 MPa和0.52 MPa與0.35 MPa,說明塊體與基巖間的黏聚力對(duì)結(jié)構(gòu)動(dòng)態(tài)特性起到的作用。

圖10?、耦愊噜弶K體剪切應(yīng)力Fig. 10　Shear stress of adjacent blocks of type Ⅰ

3.3Ⅱ類相鄰塊體響應(yīng)分析

圖11給出了2種抗滑模型條件下兩塊體運(yùn)動(dòng)的水平X向相對(duì)位移?？梢钥闯?修正模型條件下的位移在數(shù)值上相對(duì)小于經(jīng)典模型計(jì)算結(jié)果,外側(cè)塊體與內(nèi)側(cè)塊體在整個(gè)地震階段經(jīng)典與修正模型最大水平X向位移分別為8.42 cm和8.24 cm與8.34 cm和8.13 cm,整個(gè)地震階段水平位移均方根值分別為3.68 cm和3.65 cm與3.55 cm和3.54 cm。以上計(jì)算結(jié)果說明在地震荷載作用下,發(fā)生滑動(dòng)的塊體在修正模型的影響下發(fā)生的運(yùn)動(dòng)程度要小于在經(jīng)典模型下的激勵(lì)效果。

圖11?、蝾愊噜弶K體水平X向相對(duì)位移Fig. 11　Relative displacement of adjacent blocks of type Ⅱ in horizontal X direction

圖12給出了2種抗滑模型條件下兩塊體與基巖接觸面上代表點(diǎn)的剪切應(yīng)力對(duì)比。由圖12可知,基于修正模型計(jì)算結(jié)果在數(shù)值上要大于經(jīng)典模型結(jié)果,外側(cè)與內(nèi)側(cè)塊體在修正模型條件下最大剪切應(yīng)力值分別為0.63 MPa與0.094 MPa,經(jīng)典模型為0.49 MPa與0.099 MPa,平均剪切應(yīng)力為0.44 MPa與0.072 MPa和0.36 MPa與0.043 MPa,進(jìn)一步說明塊體與基巖間的黏聚力對(duì)結(jié)構(gòu)動(dòng)態(tài)特性起到的作用。

圖12?、蝾愊噜弶K體剪切應(yīng)力Fig. 12　Shear stress of adjacent blocks of type Ⅱ

4結(jié)語

本文提出的抗滑模型考慮塊體間的相互影響與作用機(jī)理,將單一塊體穩(wěn)定性問題擴(kuò)展到相鄰塊體上,能夠更好地模擬實(shí)際工程中塊體周邊接觸面上的摩擦與黏結(jié)性質(zhì),真實(shí)地反映局部塊體在強(qiáng)震作用下的動(dòng)態(tài)響應(yīng)特性?；跀?shù)值模型響應(yīng)結(jié)果,修正抗滑模型水平X向相對(duì)位移計(jì)算結(jié)果相對(duì)小于經(jīng)典抗滑模型計(jì)算的動(dòng)力響應(yīng),說明地下洞室相鄰塊體周邊接觸面動(dòng)力抗滑模型的可實(shí)現(xiàn)性與適用性。針對(duì)工程實(shí)例,全面分析了相鄰塊體結(jié)構(gòu)在經(jīng)典抗滑模型與修正抗滑模型條件下的地震響應(yīng)與動(dòng)態(tài)特性。研究表明，本文建立的抗滑模型能夠更好地模擬實(shí)際工程中塊體周邊接觸面上的抗滑與黏結(jié)性質(zhì),具有較好的工程應(yīng)用價(jià)值。

參考文獻(xiàn)：

[1] GOODMAN R E, SHI G H. Block theory and its application to rock engineering[M]. Englewood Cliffs:Prentice-Hall,1985.

[2] GOODMAN R E,TAYLOR R L,BREKKE T L. A model for the mechanics of jointed rock[J]. Journal of Soil Mechanics & Foundations Div,1968,94:637-660.

[3] 石根華.巖體穩(wěn)定分析的幾何方法[J].中國科學(xué),1981(4):487-495.(SHI Genhua. Geometric method of stability analysis of rockmass[J].Science China,1981(4):487-495.(in Chinese))

[4] 鄭銀河,夏露,于青春.考慮巖橋破壞的塊體穩(wěn)定性分析方法[J].巖土力學(xué),2013,34(增刊1):197-203.(ZHENG Yinhe, XIA Lu, YU Qingchun. Stability analysis method of block considering cracking of rock bridge[J].Rock and Soil Mechanics,2013,34(Sup1):197-203.(in Chinese))

[5] BEYABANAKI S A R,MIKOLA R G,HATAMI K. Three-dimensional discontinuous deformation analysis (3-D DDA)using a new contact resolution algorithm [J].Computers and Geotechnics,2008,35(3):346-356.

[6] OLIVEIRA D A F,INDRARATNA B. Comparison between models of rock discontinuity strength and deformation[J]. Journal of Geotechnical and Geoenvironmental Engineering,2010,136(6):864-874.

[7] ZHONG Denghua, LI Mingchao, LIU Jie. 3D integrated modeling approach togeo-engineering objects of hydraulic and hydroelectric projects[J].Science in China Series E (Technological Science),2007,50(3):329-342.

[8] 鐘登華,魯文妍,劉杰,等. 基于三維地質(zhì)模型的地下洞室曲面塊體分析[J]. 巖石力學(xué)與工程學(xué)報(bào),2011,30(增刊2):3696-3702.(ZHONG Denghua, LU Wenyan, LIU Jie, et al. Surface-block analysis of underground cavern based on 3D geological modeling[J]. Chinese Journal of Rock Mechanics and Engineering,2011,30(Sup2):3696-3702.(in Chinese))

[9] KAIDI S, ROUAINIA M, OUAHSINE A. Stability of breakwaters under hydrodynamic loading using a coupled DDA-FEM approach[J]. Ocean Engineering,2012,55(15):62-70.

[10] 鐘登華,魯文妍,劉杰,等. 復(fù)雜地質(zhì)條件下地下洞室曲面塊體地震響應(yīng)分析[J]. 天津大學(xué)學(xué)報(bào)(自然科學(xué)與工程技術(shù)版),2014,47(6):471-478.(ZHONG Denghua, LU Wenyan, LIU Jie,et al. Surface-block identification and seismic response analysis of underground structures under complicated geological conditions[J].Journal of Tianjin University(Science and Technology), 2014,47(6):471-478.(in Chinese))

[11] 馬永政,鄭宏,朱合華,等. DDA 法計(jì)算邊坡安全系數(shù)的黏聚力影響分析[J]. 巖土工程學(xué)報(bào),2009,31(7):1088-1093.(MA Yongzheng, ZHENG Hong, ZHU Hehua, et al. Effect of cohesion on evaluating slope stability factor of safety by DDA method[J]. Chinese Journal of Geotechnical Engineering, 2009, 31(7):1088-1093.(in Chinese))

[12] 付曉東,盛謙,張勇慧. 基于矢量和一非連續(xù)變形分析的滑坡安全系數(shù)計(jì)算方法研究[J]. 巖石力學(xué)與工程學(xué)報(bào),2014,33(增刊2):4122-4128.(FU Xiaodong, SHENG Qian, ZHANG Yonghui. Research on safety factor computation of landslide based on vector sum-discontinuous deformation analysis method[J]. Chinese Journal of Rock Mechanics and Engineering,2014, 33(Sup2):4122-4128.(in Chinese))

[13] SU Huaizhi, HU Jiang, LI Jinyou, et al. Deep stability evaluation of high-gravity dam under combining action of powerhouse and dam[J]. International Journal of Geomechanics,2013,13:257-272.

[14] JIANG Yusong, SU Chao. Mixed finite element method for contact problems of multibody[J]. Earth and Space,2010: 606-620.

[15] 魯文妍.強(qiáng)震作用下大型地下廠房洞室群災(zāi)變仿真研究[D].天津:天津大學(xué),2012.

Seismic responses of adjacent blocks in large-scale underground cavern group

LI Wenqian1, TONG Dawei1, WANG Zhen1, ZHU Xiaobin1, LU Wenyan1,2

(1.StateKeyLaboratoryofHydraulicEngineeringSimulationandSafety,TianjinUniversity,Tianjin300072,China;2.MaterialsandStructuralEngineeringDepartment,NanjingHydraulicResearchInstitute,Nanjing210029,China)

Abstract:In order to simulate the movement and dynamic response of adjacent blocks in an underground cavern group under seismic loads, a dynamic anti-sliding model, considering the bond strength between the blocks, was established. The effectiveness of the anti-sliding model in comparison to other models was analyzed using a numerical model. In the practical engineering examples, the calculated root mean square values of horizontal displacement in the X direction of two adjacent blocks of type I were 2.16 cm and 2.02 cm when the classical anti-sliding model was used, and 2.00 cm and 1.94 cm when the modified anti-sliding model was used throughout the earthquake process. The corresponding simulated results of two adjacent blocks of type II were 3.68 cm and 3.65 cm when the classical model was used and 3.55 cm and 3.54 cm when the modified model was used. Hence, it is demonstrated that adjacent blocks considering bonding characteristics show much better stability, and the proposed model has high application value.

Key words:underground cavern; adjacent blocks; bonding characteristics; anti-sliding model; seismic responses

DOI：10.3876/j.issn.1000-1980.2016.03.002

收稿日期：2015-10-02

基金項(xiàng)目：國家自然科學(xué)基金創(chuàng)新研究群體科學(xué)基金(51321065)；國家重點(diǎn)基礎(chǔ)研究發(fā)展計(jì)劃(973計(jì)劃)(2013CB035904)；天津市自然科學(xué)基金青年項(xiàng)目(13JCQNJC08900)

作者簡介：李文倩(1986—)，女，天津人，博士研究生，主要從事地下洞室模擬與評(píng)價(jià)分析研究。E-mail：bxqe@163.com 通信作者： 佟大威，講師。E-mail：tongdw@tju.edu.cn

中圖分類號(hào)：TV312

文獻(xiàn)標(biāo)志碼：A

文章編號(hào)：1000-1980(2016)03-0196-07