張興仙

[摘要]本文研究的工程旁距離大約250m處有一斷層經(jīng)過。經(jīng)歷八級(jí)地震后，斷層垂直抬升約3.5m，水平錯(cuò)動(dòng)2～3m，在工程中出現(xiàn)了次級(jí)破碎帶，為模擬該斷層錯(cuò)動(dòng)及次級(jí)破碎帶的存在對(duì)工程產(chǎn)生的影響，采用有限元軟件MIDAS/GTS對(duì)該工程及工程所在山體整體進(jìn)行數(shù)值模擬。將F1斷層的水平和豎直錯(cuò)動(dòng)以施加水平和豎向位移場的方式進(jìn)行模擬，以軟弱層的方式模擬次級(jí)破碎帶。研究次級(jí)破碎帶對(duì)工程的影響。

[關(guān)鍵詞]數(shù)值模擬 斷層 破碎帶 位移場

0概述

地下結(jié)構(gòu)周圍受巖土體介質(zhì)的約束，地震發(fā)生時(shí)隨著巖土體一起運(yùn)動(dòng)而具有良好的抗震性能，受地震的破壞與地面結(jié)構(gòu)相比較輕。但是二十世紀(jì)以來，破壞性大地震頻繁發(fā)生，均對(duì)地下結(jié)構(gòu)造成了不同程度的破壞，使得大量專家、學(xué)者開始廣泛、深入、系統(tǒng)的開展對(duì)地下工程抗震的研究。由于受周圍介質(zhì)的約束，地下結(jié)構(gòu)的地震動(dòng)反應(yīng)特點(diǎn)與地面結(jié)構(gòu)明顯不同，所以兩者的抗震分析方法也不相同。隨著地下結(jié)構(gòu)抗震理論和計(jì)算機(jī)技術(shù)的飛速發(fā)展，有限元?jiǎng)恿Ψ磻?yīng)分析的應(yīng)用逐漸變得普遍，相比靜力法而言它可以更全面更真實(shí)的再現(xiàn)地下結(jié)構(gòu)在地震荷載下的動(dòng)態(tài)特性。該方法在計(jì)算模型及參數(shù)的確定、地震波及其輸入方式的確定等方面還存在一些問題，有待于進(jìn)一步完善。本文主要是對(duì)實(shí)際工程進(jìn)行三維模擬，模型采用實(shí)體工程的尺寸進(jìn)行建模，結(jié)構(gòu)形式較復(fù)雜，單元網(wǎng)格數(shù)量龐大，計(jì)算機(jī)計(jì)算條件還不能滿足有限元?jiǎng)恿Ψ治?，所以采用擬靜力法進(jìn)行分析研究。

1工程概況

工程場地位于青藏高原東緣的龍門山構(gòu)造帶內(nèi)，具體處于龍門山構(gòu)造帶中的主干斷裂北川映秀斷裂附近。2008年5月21號(hào)汶川8.0級(jí)地震就發(fā)生在龍門山構(gòu)造帶上，主干發(fā)震斷裂為北川-映秀斷裂，并同時(shí)牽動(dòng)彭縣灌縣斷裂發(fā)生了同震地表斷裂，給龍門山及其附近的四川盆地西緣地區(qū)帶來了毀滅性的破壞，也本工程造成了嚴(yán)重的破壞。在場地附近產(chǎn)生了北東向的地表破裂，主地表破裂在場地附近寬約30m，總體呈N50-60°E方向延伸，表現(xiàn)為逆沖性質(zhì)，最大垂直位錯(cuò)約3.5m，主地表破裂帶未通過工程場地，它距離工程的最近距離約250m，但工程場地處于發(fā)震斷層的上盤，在地震時(shí)主破裂帶強(qiáng)烈變形過程中，位于斷層上盤的工程場地內(nèi)部可能產(chǎn)生了次級(jí)破裂，致使工程內(nèi)部遭受了嚴(yán)重破壞，如圖1所示。

2有限元模型的建立

2.1基本假定

工程場地處于發(fā)震斷層的上盤，在地震時(shí)主破裂帶強(qiáng)烈變形過程中，位于斷層上盤的工程場地內(nèi)部可能產(chǎn)生了次級(jí)破裂，致使工程內(nèi)部遭受了嚴(yán)重破壞，本文假設(shè)次級(jí)破碎帶寬度0.5m。

本文還對(duì)模型做了部分假定，以便于將復(fù)雜抽象的實(shí)際問題轉(zhuǎn)化為具體簡單的數(shù)值分析模型，只研究幾個(gè)特定的因素對(duì)坑道的影響。本文所做的基本假定如下：

（1）假定沿著坑道的軸線方向，周圍巖體的特性與坑道的特性并不發(fā)生改變，地質(zhì)與地形條件相同，建立簡化的三維模型。

（2）假定巖層各向同性，沿水平方向均勻分布，不考慮地下水滲流的影響。

（3）因?yàn)榭紤]的是位移場對(duì)使用的坑道的響應(yīng)，所以不考慮施工階段的影響。

（4）假設(shè)整個(gè)模型的各部分都是在線彈性范圍內(nèi)工作，所以采用靜力分析中的線彈性分析。

2.2幾何模型的建立

本文采用有限元軟件MIDSA/GTS和制圖軟件AutoCAD進(jìn)行模型的建立。坑道襯砌模型的建立基于實(shí)際工程的設(shè)計(jì)圖紙，依據(jù)圖紙上的尺寸大小建模。工程的軸線圖如圖2所示，隧道截面形狀如圖3所示。在建立幾何模型的時(shí)候，把把工程所在的山體采用圓錐體的形狀來模擬。幾何模型的建立大體分為三步，第一步分別拉伸導(dǎo)入的隧道截面形狀內(nèi)外輪廓線組為實(shí)體；第二步建立能包圍住隧道實(shí)體的“圓餅”，“圓餅”的建立如圖4所示。工程坑道埋深225m，山體坡度為35°，山體模型的建立如圖5所示，地下隧道的模型建立如圖6所示。

2.3材料參數(shù)及網(wǎng)格劃分

工程模型中的材料特性如表1所示。網(wǎng)格劃分作為有限元模型建立的關(guān)鍵環(huán)節(jié)，其網(wǎng)格形式及單元尺寸對(duì)計(jì)算精度和計(jì)算規(guī)模將產(chǎn)生直接的影響。本文中網(wǎng)格劃分也分為三個(gè)層面應(yīng)。第一層面是對(duì)坑道襯砌的網(wǎng)格劃分。根據(jù)精度要求，在橫截面方向劃一層網(wǎng)格，尺寸為200-300mm，而在縱向上采用的是漸變式的方法，即端頭較密，中間較疏。這樣做既能夠保證所要分析結(jié)果的精確性，又能夠減少網(wǎng)格的數(shù)量，提高計(jì)算的效率。見圖7為襯砌單元。第二個(gè)層面是對(duì)“圓餅”山體的網(wǎng)格劃分。主要是把與坑道接近的這一層巖石劃得較細(xì)一些，外側(cè)采用漸進(jìn)的方式，為內(nèi)側(cè)尺寸的兩到三倍。其余部分就可劃到內(nèi)側(cè)尺寸的四到六倍，見圖8所示。最后一個(gè)層面是對(duì)“圓餅”以外山體（包括斷層破碎帶及隔震層）的網(wǎng)格化分。因?yàn)檫@部分不是分析研究的主要對(duì)象，為了節(jié)約計(jì)算時(shí)間，網(wǎng)格就會(huì)劃得比較稀疏，見圖9所示。本文中整個(gè)模型劃分的單元控制在了50萬內(nèi)，其中巖體單元40多萬個(gè)、襯砌單元7萬多個(gè)。

2.4外荷載情況及邊界條件

在距1#口所在平面250m的地方模擬F1斷裂帶，在其上施加沿截面的水平位移場和豎向位移場來模擬地震作用過程中斷層的運(yùn)動(dòng)并使之作為外荷載條件。所施加的位移場的位置和方向見圖10。邊界條件主要有三組，邊界條件一為距離施加強(qiáng)制位移（外荷載）處的F1斷層580m的上盤底面固定X，Y，Z向位移，此邊界下破碎帶部分也在固定范圍內(nèi)，即破碎帶不完全，詳見圖11；邊界條件二為邊界條件一遠(yuǎn)離位移場方向平移120m，與邊界條件一相比破碎帶剛好完全釋放，詳見圖12；邊界條件三平行次級(jí)破碎帶并相距200m以外的底部固定X，Y，Z向位移，詳見圖13。

3計(jì)算結(jié)果

模型位移分析時(shí)選取監(jiān)控截面見圖14。各截面距離1#口所在截面的距離詳見表2。其中，圖14截面3、4相距5m，截面6、7相距5m，截面4、6相距10.5m。把次級(jí)破碎帶看做是一個(gè)規(guī)模較小的斷層，那么次級(jí)破碎帶左側(cè)的部分稱為上盤，右側(cè)的部分稱為下盤。

邊界條件二的各點(diǎn)的位移反應(yīng)均比邊界條件一、三的情況大，這與邊界條件固定范圍有關(guān)，邊界條件二較邊界條件一和三固定范圍較小，所以位移反應(yīng)較大。且邊界條件二為破碎帶完全釋放（破碎帶完全不在固定的邊范圍內(nèi)），是三種邊界中最合理的邊界。

結(jié)構(gòu)的整體變形與施加位移場的方向一致，距離F1斷層帶越遠(yuǎn)，結(jié)構(gòu)的位移反應(yīng)越小，這表現(xiàn)出遠(yuǎn)離斷層錯(cuò)動(dòng)帶震害逐漸減小，結(jié)構(gòu)的變形方向與斷層的錯(cuò)動(dòng)方向密切相關(guān)。在實(shí)際工程中，應(yīng)該加強(qiáng)近斷層處一定范圍的地下結(jié)構(gòu)的抗震措施，減輕地震中斷層的活動(dòng)對(duì)結(jié)構(gòu)的破壞性損壞。

與無破碎帶結(jié)構(gòu)相比，破碎帶的存在導(dǎo)致結(jié)構(gòu)的變形不再均勻變化，在破碎帶處的1#隧道的襯砌變形較大造成剪斷破壞，此處的位移明顯突變。

次級(jí)破碎帶的存在對(duì)結(jié)構(gòu)的位移影響較大，對(duì)于下盤結(jié)構(gòu)的位移明顯增大，而對(duì)于上盤結(jié)構(gòu)的位移與無破碎帶相比卻減小。破碎帶的出現(xiàn)對(duì)于遠(yuǎn)端來說吸收了大部分能量，減輕了遠(yuǎn)端的破壞。

4結(jié)論

從計(jì)算結(jié)果分析來看，首先影響地下結(jié)構(gòu)位移的主要原因是地震時(shí)斷層的活動(dòng)，然而次級(jí)破碎帶的存在改變了原來結(jié)構(gòu)的位移大小及分布，使斷層近端增大了結(jié)構(gòu)的位移反應(yīng)，而在斷層與破碎帶之外（遠(yuǎn)端）則起到耗能減震的效果，反而有減小結(jié)構(gòu)的位移反應(yīng)的趨勢(shì)。充分說明破碎帶在地震作用中，由于其剛度較小，變形較大，吸收了大部分能量，導(dǎo)致破碎帶與斷層之間的山體及結(jié)構(gòu)變形增大，而在破碎帶與斷層之外的部分結(jié)構(gòu)及山體則由于破碎帶的耗能阻斷作用而得到保護(hù)，變形反而成減小趨勢(shì)。