張 碧，呂曉春

(1.陜西鐵路工程職業(yè)技術(shù)學(xué)院，陜西渭南714000；2.華北水利水電大學(xué)資源與環(huán)境學(xué)院，河南鄭州450011)

?

水工隧道動力深、淺埋劃分限界研究及試驗驗證

張 碧1，呂曉春2

(1.陜西鐵路工程職業(yè)技術(shù)學(xué)院，陜西渭南714000；2.華北水利水電大學(xué)資源與環(huán)境學(xué)院，河南鄭州450011)

隧道動力深、淺埋劃分限界屬于隧道結(jié)構(gòu)的固有特性，受圍巖類別、隧道結(jié)構(gòu)形狀、隧道跨度等因素的影響，不受外界因素(如地震波種類等)的影響。計算分析和振動臺模型試驗驗證表明，隧道動力深、淺埋劃分限界值受圍巖類別的影響較大；矩形隧道對其敏感性較大，圓形隧道的敏感性較??；隧道跨度對其有一定的影響。

隧道；動力劃分限界；振動臺試驗；地震波

0 引 言

截至目前，全球發(fā)生過許多次重大地震。在這些地震中，隧道結(jié)構(gòu)都發(fā)生了不同程度的破壞。震害統(tǒng)計發(fā)現(xiàn)，隧道埋深大于50 m 時破壞程度明顯減小，在300 m以下沒有嚴(yán)重的破壞，受破壞的地下洞室有35%的埋深小于50 m[1]。在地震作用下，隧道埋深對隧道的受力和變形有很大影響，對深、淺埋隧道需要采用不同的計算方法。目前，公路隧道抗震設(shè)計中或者沒有區(qū)分深、淺埋，或者仍然采用靜力分析中深、淺埋劃分方法[2]，不能反應(yīng)出地震波作用的動力實質(zhì)。因此，進行隧道動力深、淺埋劃分限界研究是必要的。

隧道埋深是影響地震動力響應(yīng)的一個重要因素。舒亞俐等[3]對不同埋深構(gòu)筑物的動力響應(yīng)進行研究認(rèn)為，構(gòu)筑物的振動頻率隨埋深的增加而增大，總體位移隨埋深的增加而減少，地震慣性力隨埋深增加而減??；鄭永來等[4]對國內(nèi)外城市地下結(jié)構(gòu)的震害情況進行統(tǒng)計分析后認(rèn)為，地下管道的破壞隨埋深的增加而減小，地下結(jié)構(gòu)埋深越大，由地震面波導(dǎo)致的能量越小，震害應(yīng)減輕；于翔[5- 6]分析了1995年阪神地震中地下結(jié)構(gòu)的破壞情況發(fā)現(xiàn)，上層破壞比下層破壞嚴(yán)重，地下結(jié)構(gòu)上部土層越厚，破壞越輕；李海波等[7]研究了地震荷載作用下地下巖體洞室位移特征的影響因素認(rèn)為，地下洞室位移響應(yīng)隨埋深的增加而減小，洞室位移響應(yīng)存在臨界埋深；李加林等[8]運用FLAC就珠三角軟土地層隧道埋深對地鐵運行誘發(fā)的相關(guān)振動影響規(guī)律進行分析表明，埋深越淺則動力響應(yīng)越強，埋深越大則動力響應(yīng)越弱；鄭穎人等[9]通過模型試驗與數(shù)值分析方法研究了靜力條件下隧道從淺埋到深埋的破壞過程，根據(jù)壓力拱(普氏壓力拱) 確定深、淺埋的分界線。

以上研究表明，隧道埋深對隧道地震響應(yīng)有很大的影響，但依靠靜力條件下的劃分標(biāo)準(zhǔn)確定地震作用下隧道的動力深、淺埋是不準(zhǔn)確的，有關(guān)隧道動力深、淺埋劃分方面的研究還未見相關(guān)報道。為此，本文通過數(shù)值計算和振動臺模型試驗，劃分隧道動力深、淺埋。

1 數(shù)值計算原理與方法

計算中將隧道作為平面應(yīng)變問題處理，隧道襯砌和圍巖均采用塊體(Plane42 單元，二維實體結(jié)構(gòu)模型用作平面應(yīng)變單元)，圍巖的破壞準(zhǔn)則采用D-P準(zhǔn)則，計算模型的尺寸高h×寬b=300 m×500 m，模型左、右兩側(cè)使用粘彈性邊界，使用彈簧單元combin14模擬。計算模型見圖1。具體的彈性和阻尼參數(shù)計算如下

(1)

CBT=ρcS

CBN=ρcP

式中，KBT、KBN分別為切向和法向彈簧的彈性系數(shù)；CBT、CBN分別為切向和法向的阻尼系數(shù)；ρ為圍巖的密度；R為邊界點到隧道中心的距離；cS、cP分別為地震橫波波和縱波波速；kT、kN為調(diào)整系數(shù)，在0.5～1.2的范圍內(nèi)取值，通過數(shù)值試算確定[10]。

圖1 計算模型

為了對比分析，隧道采用矩形、馬蹄形和圓形隧道結(jié)構(gòu)。動力阻尼C采用瑞利阻尼，即

C=αM+βK

(2)

式中，M為系統(tǒng)質(zhì)量矩陣；K為剛度矩陣；α為與質(zhì)量成比例的阻尼系數(shù)；β為與剛度成比例的阻尼系數(shù)。

對于多自由度系統(tǒng)，任意系統(tǒng)角頻率ωi的臨界阻尼比ξi是角頻率為ωi的第i階模態(tài)臨界阻尼的一部分，可以從下式得到

(3)

模擬計算隧道結(jié)構(gòu)在El_Centro地震波(前8.0s)和2008年汶川地震波作用下的地震響應(yīng)，兩者的頻率和峰值分布不同，對比分析可以反應(yīng)出隧道結(jié)構(gòu)的波動特性。加速度峰值均取0.3g。相關(guān)輸入地震波時程見圖2。

圖2 輸入地震波時程

2 數(shù)值計算結(jié)果分析

2.1 不同圍巖類別

本節(jié)計算時，選取隧道結(jié)構(gòu)為馬蹄形，輸入El_Centro地震波，提取了襯砌外側(cè)拱頂處的剪切應(yīng)力τxy。不同圍巖類別襯砌應(yīng)力隨隧道埋深變化見圖3。圖中，B為隧道跨度。

圖3 不同圍巖類別襯砌應(yīng)力隨隧道埋深變化

從圖3可知，在地震波作用下，隧道襯砌的剪應(yīng)力在隨埋深變化中存在拐點，拐點處的隧道埋深可作為隧道動力深、淺埋限界；隧道動力深、淺埋限界受圍巖類別影響較大；隧道跨度對隧道動力深、淺埋限界有一定的影響(圖3d中，埋深150 m處數(shù)值有突變，應(yīng)是計算誤差所致，可忽略)。

2.2 不同隧道結(jié)構(gòu)形狀和跨度對比

為了分析影響隧道動力深、淺埋限界的影響因素，本節(jié)考慮了隧道結(jié)構(gòu)形狀和隧道跨度2種因素進行了對比分析。不同隧道結(jié)構(gòu)形狀襯砌應(yīng)力隨隧道埋深變化見圖4。

圖4 襯砌應(yīng)力隨隧道埋深變化

由圖4可知，隧道結(jié)構(gòu)形狀對隧道動力深、淺埋限界的影響不大；在不同的隧道結(jié)構(gòu)形狀條件下，隧道跨度對隧道襯砌動應(yīng)力和隧道動力深、淺埋限界的敏感性不同。其中，圓形隧道的敏感性較小，矩形隧道的敏感性最大；隧道跨度對隧道動力深、淺埋限界有一定的影響，但對不同結(jié)構(gòu)形狀隧道的影響不同，對矩形隧道影響最大，對圓形隧道影響最小。

2.3 不同地震波對比

由上述分析可知，矩形隧道襯砌動應(yīng)力對隧道動力深、淺埋限界的敏感性最強。本節(jié)計算分析了矩形隧道在El_Centro波和汶川波2種不同地震波(同為0.3g)作用下隧道動力深、淺埋限界的變化(見圖5)。

圖5 不同地震波作用下襯砌應(yīng)力隨隧道埋深變化

從圖5可知，地震波種類的不同對隧道動力深、淺埋限界的影響不明顯；隧道動力深、淺埋限界應(yīng)該是隧道自身的特性，與所受的地震波種類關(guān)系不大。

3 振動臺模型試驗設(shè)計

本試驗是為了驗證地震作用下隧道結(jié)構(gòu)的響應(yīng)規(guī)律，屬于動力試驗，應(yīng)該用到動力相似準(zhǔn)則。動力學(xué)運動方程可用矩陣方程表示如下

[M][a]+[C][v]+[K][u]=[F]

(4)

式中，[M]為質(zhì)量矩陣；[a]為加速度矩陣；[C]為阻尼矩陣；[v]為速度矩陣；[K]為剛度矩陣；[u]為位移矩陣；[F]為荷載矩陣。

本試驗為對比驗證性試驗，綜合考慮振動臺承載能力、臺面尺寸、具體隧道尺寸和制作工藝等多種因素，對隧道原型進行適當(dāng)簡化，試驗用長50 cm、寬50 cm、高80 cm的模型箱作為隧道洞口段模型，幾何相似比為1∶40；隧道襯砌采用了圓形襯砌，外徑15 cm，襯砌厚度1 cm。試驗?zāi)Ｐ鸵妶D6。圖中虛線所示的襯砌為在不同深度處埋設(shè)襯砌，以模擬不同隧道埋深情況。

圖6 振動臺試驗?zāi)Ｐ?單位：cm)

振動臺提供水平方向的地震波輸入，輸入的地震波為El_Centro地震波(前8.0 s)和2008年汶川地震波，地震波的加速度峰值均取0.1g、0.2g、0.3g和0.4g。

用高強石膏模擬襯砌，用粉煤灰、石英砂、石膏、鋸末和機油的混合物模擬圍巖。石膏和機油主要等效圍巖的粘聚力c值[11]，粉煤灰和石英砂作為主要的骨料，通過調(diào)節(jié)其密度使顆粒級配滿足試驗要求。在確定襯砌材料時，根據(jù)水膏比為0.8∶1、0.9∶1、1∶1、1.1∶1和1.2∶1等5種情況分別做成標(biāo)準(zhǔn)圓柱體試件，然后放在萬能試驗機上進行分級加載壓縮試驗，確定其彈性模量在0.5～2.9 GPa。實際隧道一般采用C25混凝土作為二襯材料，其彈性模量為17.5～30 GPa。確定襯砌材料的彈性模量相似比為1∶20，據(jù)此選用水膏比為1.05∶1。將粉煤灰、石英砂、石膏、鋸末和機油等材料按照不同的配合比進行混合[12]，參照實際圍巖材料的物理性質(zhì)指標(biāo)，進行正交試驗來確定試驗材料的相似比。

將隧道襯砌外側(cè)上(Y1)、下(Y2)、左(Y3)和右(Y4)4個部位(見圖7)的剪切應(yīng)力測試結(jié)果統(tǒng)計繪圖(見圖8)。振動臺試驗結(jié)果與前述的數(shù)值計算結(jié)果對比分析可知，隧道振動臺模型試驗進一步驗證了隧道動力深、淺埋限界的存在；振動臺模型試驗的結(jié)果與數(shù)值計算結(jié)果基本一致(限界值均在45～50 m)；隧道襯砌不同部位(上、下、左、右)處的動力深、淺埋限界值基本一致。

圖7 測試點布置

圖8 振動臺試驗結(jié)果

4 結(jié) 語

通過計算分析和振動臺模型試驗，證明了隧道動力深、淺埋限界的存在，并分析了其影響因素。得到如下結(jié)論：

(1)計算分析發(fā)現(xiàn)隧道動力深、淺埋限界是確實存在的，且與隧道的靜力深、淺埋限界值有所區(qū)別。

(2)隧道動力深、淺埋限界與圍巖類別、隧道結(jié)構(gòu)形狀、隧道跨度等因素有關(guān)，屬于隧道結(jié)構(gòu)的內(nèi)在特性，與所受地震波的種類無關(guān)。

(3)隧道動力深、淺埋限界值受圍巖類別的影響較大；隧道結(jié)構(gòu)形狀影響因素中，圓形隧道的敏感性較小，矩形隧道的敏感性最大；隧道跨度對動力深、淺埋限界值有一定的影響。

(4)隧道振動臺模型試驗驗證了隧道動力深、淺埋限界的存在，并且得到了與數(shù)值計算一致的結(jié)論。

[1]SHARMA S， JUDD W R. 地震對地下洞室的破壞[J]. 雷謙榮， 譯. 地下空間， 1992， 12(4)： 335- 344．

[2]JTG/T D70—2010 公路隧道設(shè)計細則[S].

[3]舒亞俐. 構(gòu)筑物不同埋深時地震動反分析探討[J]. 特種結(jié)構(gòu)， 1997， 14(4)： 49- 55.

[4]鄭永來， 楊林德. 地下結(jié)構(gòu)震害與抗震對策[J]. 巖石力學(xué)與工程學(xué)報. 2001， 20(1)： 70- 73．

[5]于翔. 地鐵建設(shè)中應(yīng)充分考慮抗地震作用—阪神地震破壞的啟示[J]. 鐵道建筑技術(shù)， 2000(6)： 32- 35．

[6]于翔. 地下建筑結(jié)構(gòu)應(yīng)充分考慮抗震問題—1995年阪神地震破壞的啟示[J]. 工程抗震， 2002(4)： 17- 20.

[7]李海波， 馬行東， 李俊如， 等. 地震荷載作用下地下巖體洞室位移特征的影響因素分析[J]. 巖土工程學(xué)報， 2006， 28(3)： 358- 362.

[8]李加林， 鄧飛皇. 隧道埋深對地鐵運行誘發(fā)振動的影響分析[J]. 科學(xué)技術(shù)與工程， 2007， 7(15)： 3819- 3822．[9]鄭穎人， 徐浩， 王成， 等. 隧洞破壞機理及深、淺埋分界標(biāo)準(zhǔn)[J]. 浙江大學(xué)學(xué)報： 工學(xué)版， 2010， 44(10)： 1851- 1856．

[10]杜修力， 趙密， 王進廷. 近場波動模擬的人工應(yīng)力邊界條件[J]. 力學(xué)學(xué)報， 2006， 38(1)： 49- 56．

[11]王明年， 林國進， 于麗， 等. 隧道抗震與減震[M]. 北京： 科學(xué)出版社， 2012.

[12]蔣樹屏， 蔣華， 王曉雯， 等. 高烈度地震區(qū)公路隧道洞口段大型振動臺模型試驗方案設(shè)計[J]. 公路， 2009(10)： 245- 249．

(責(zé)任編輯 楊 健)

Study on the Dynamic Bound of Shallow Buried and Deep Buried Hydraulic Tunnel and Its Experimental Verification

ZHANG Bi1, Lü Xiaochun2

(1. Shaanxi Railway Institute, Weinan 714000, Shaanxi, China; 2. School of Resources & Environment,North China University of Water Resources and Electric Power, Zhengzhou 450011, Henan, China)

The dynamic bound of shallow buried and deep buried hydraulic tunnel is the intrinsic characteristics of tunnel structure, which is affected by surrounding rock types, tunnel structure shapes, tunnel span and others, and not affected by external factors (such as the type of seismic wave). The calculation and shaking table test show that, the dynamic bound values of shallow buried and deep buried tunnel are strongly influenced by rock type, the rectangular tunnel has big sensitivity and in contrary the circular tunnel has small one, and the tunnel span has a certain influence on dynamic bound values．

tunnel; dynamic bound; shaking table test; seismic wave

2016- 03- 03

國家自然科學(xué)基金資助項目(81408222)

張碧(1979—)，男，山西沁縣人，副教授，碩士，主要從事道路與鐵道工程專業(yè)的教學(xué)和研究工作．

TV312

A

0559- 9342(2016)06- 0026- 04