李兆琦 劉 洋 張朋天 王小敬

(1.河北建筑工程學(xué)院土木工程學(xué)院，河北 張家口 075000；2.北旺建設(shè)集團(tuán)有限公司，河北 承德 067000；3.河北省土木工程診斷、改造與抗災(zāi)重點(diǎn)實(shí)驗(yàn)室，河北 張家口 075000)

0 引 言

隨著隧道施工技術(shù)的不斷提升，我國(guó)隧道建設(shè)得到長(zhǎng)足的發(fā)展.目前，國(guó)內(nèi)在建隧道下穿建筑物的方式多為斜穿，且以下穿現(xiàn)代建筑物為主，針對(duì)軟弱圍巖地區(qū)大斷面下穿古長(zhǎng)城隧道的施工研究較少.為達(dá)到國(guó)家文物局對(duì)于古長(zhǎng)城保護(hù)的要求，需嚴(yán)格控制圍巖變形與長(zhǎng)城結(jié)構(gòu)的沉降.本文以某公路隧道下穿古長(zhǎng)城工程項(xiàng)目為案例，采用有限元軟件FLAC3D進(jìn)行數(shù)值模擬，針對(duì)在CD法、CRD法和三臺(tái)階法三種開挖工法下測(cè)得的監(jiān)測(cè)點(diǎn)沉降值對(duì)上部長(zhǎng)城的影響進(jìn)行了對(duì)比分析，選擇出更為合理的開挖方案，為今后的此類工程提供出解決的范例.

許多學(xué)者所做的研究為大斷面隧道的開挖模擬提供了必要的依據(jù)，如張孟帥等[1]通過Ansys軟件建立動(dòng)態(tài)施工模型分析，得到CD法和CRD法應(yīng)用于大斷面隧道時(shí)產(chǎn)生不同的力學(xué)響應(yīng)；來弘鵬等[2]對(duì)西河口大斷面隧道下穿明長(zhǎng)城工程進(jìn)行研究，通過室內(nèi)模擬實(shí)驗(yàn)和現(xiàn)場(chǎng)監(jiān)控沉降確定其施工方案的可行性；劉誠(chéng)等[3]對(duì)湖北省某條隧道使用有限元軟件FLAC3D進(jìn)行數(shù)值模擬，并分析多種開挖工法對(duì)隧道穩(wěn)定性的影響；孫星亮[4]對(duì)廣乘山隧道工程采用密排超前小導(dǎo)管注漿加固，進(jìn)行有限元數(shù)值模擬，并用數(shù)值計(jì)算結(jié)果和監(jiān)測(cè)結(jié)果證明所選預(yù)加固措施的有效性；唐思聰?shù)萚5]以人和場(chǎng)隧道為例，通過對(duì)比分析隧道塑性區(qū)的分布特征及開挖產(chǎn)生的擾動(dòng)效應(yīng)、圍巖變形的時(shí)空效應(yīng)等，得出安全、經(jīng)濟(jì)的最佳施工方案.

1 工程概述

某公路隧道長(zhǎng)度約為1360m，隧址區(qū)地貌單元屬侵蝕構(gòu)造低中山，地形中高側(cè)低，地勢(shì)起伏較大，坡度在10°～45°之間，平均海拔在800m左右，高差為89.02m.隧道對(duì)應(yīng)的起終點(diǎn)樁號(hào)為YK0+660～YK2+080.隧址區(qū)圍巖條件較差，為V級(jí)圍巖，巖體穩(wěn)定性較差，且存在多處斷裂.隧道從古長(zhǎng)城下穿過，隧道拱頂上地層主要以黃土為主，局部地區(qū)發(fā)現(xiàn)有卵石夾層.

2 數(shù)值模擬

2.1 建立模型

圖1 計(jì)算模型

數(shù)值模擬采用FLAC3D軟件，根據(jù)設(shè)計(jì)方案，在MIDAS軟件中建立有限元模型并劃分網(wǎng)格，將模型通過轉(zhuǎn)換軟件導(dǎo)入FLAC3D中.建立模型如圖1所示，由圣維南原理確定模型范圍為：y方向取100m，x方向上取100m(-50m～50m)，隧道埋深為50m[6].圍巖用實(shí)體單元來模擬，襯砌等支護(hù)結(jié)構(gòu)采用彈性材料提高材料參數(shù)來實(shí)現(xiàn).上部長(zhǎng)城模型用長(zhǎng)方體單元模擬，平均高度設(shè)為6m，寬度定為10m.除上邊界設(shè)定為自由邊界之外，模型的x、y方向均有法向約束約束位移.圍巖采用M-C本構(gòu)模型，控制隧道洞徑為6.5m.本文在深埋條件下，在V級(jí)圍巖地區(qū)分別采用CD法、CRD法以及三臺(tái)階法進(jìn)行隧道開挖的三維數(shù)值模擬[7-8].

根據(jù)工程資料顯示，該研究區(qū)段內(nèi)地層主要以黃土為主，力學(xué)參數(shù)參見表1[9].

表1 模型的物理力學(xué)指標(biāo)

2.2 開挖方式

CD法開挖模擬的步驟順序如下：首先在左側(cè)導(dǎo)坑開挖后進(jìn)行初期支護(hù)，然后進(jìn)行右側(cè)導(dǎo)坑開挖再施做初期支護(hù)，之后施做二次襯砌；CD法斷面如圖2所示.CRD法的開挖斷面在CD法基礎(chǔ)上設(shè)置每步增加一道臨時(shí)仰拱，CRD法斷面如圖3所示.臺(tái)階法斷面如圖4所示.

模擬三臺(tái)階法開挖的步驟順序是：按照上、中、下三層臺(tái)階依次進(jìn)行開挖，在每層開挖完成后進(jìn)行初期支護(hù)；三臺(tái)階法斷面如圖4所示，其中，上中下臺(tái)階的高度都為3m.

圖2 CD法斷面 圖3 CRD法斷面 圖4 臺(tái)階法斷面

3 結(jié)果分析

3.1 沉降分析

3.1.1 各工法開挖得到的最終沉降位移分析

本隧道作為大斷面下穿長(zhǎng)城隧道，采用3種開挖工法：CD法、CRD法和三臺(tái)階法分別對(duì)該隧道開挖與支護(hù)進(jìn)行數(shù)值模擬，開挖后的沉降位移云圖如圖5所示.

(a)CD法

(b)CRD法

(c)三臺(tái)階法

由圖5所示的位移云圖得到，CD法、CRD法和三臺(tái)階法三種工況下最大位移值分別是32.23mm、32.22mm和91.87mm，對(duì)比結(jié)果可知，在上述幾種開挖工法下，拱頂均會(huì)發(fā)生不同程度的沉降.由沉降位移云圖所示顏色深度的變化范圍[10]，云圖顯示隧道底部存在隆起較大的情況，應(yīng)及時(shí)施做仰拱初襯，并且控制襯砌材料參數(shù)以限制拱底的隆起程度，同時(shí)也可減小周邊圍巖的變形程度，對(duì)修建在軟弱圍巖上部的古長(zhǎng)城能夠控制其沉降.

同時(shí)隧道開挖后，監(jiān)測(cè)點(diǎn)的位移值如表2所示，由于CD法是先開挖左側(cè)土體，然后開挖右側(cè)土體的順序，較其他工法來說釋放土壓力更早，故左側(cè)拱底隆起值大于拱頂沉降.CRD法相比三臺(tái)階法測(cè)得拱頂沉降量減少61.7%，比CD法拱頂沉降量減少0.8%；CRD法比三臺(tái)階法在拱肩豎向位移減少52.5%，最跨度處水平收斂減少66.6%，拱腳豎向位移減少52.9%，拱底隆起減少64.9%.各關(guān)鍵點(diǎn)監(jiān)測(cè)到的位移CRD法均要小于CD法，上述兩種工法相較三臺(tái)階法對(duì)圍巖變形約束的效果都更加突出.

表2 隧道最終位移監(jiān)測(cè)值 mm

3.2 各監(jiān)測(cè)點(diǎn)沉降值對(duì)比

3.2.1 長(zhǎng)城沉降值對(duì)比

隧道開挖時(shí)會(huì)對(duì)圍巖及上部長(zhǎng)城造成一定的影響[11-13]，為保證隧道上方長(zhǎng)城的結(jié)構(gòu)安全，需嚴(yán)格控制隧道開挖時(shí)長(zhǎng)城結(jié)構(gòu)產(chǎn)生的沉降值.在本文中，為研究3種工法開挖對(duì)上部長(zhǎng)城結(jié)構(gòu)沉降值的影響，在本數(shù)值模擬中沿x軸方向在長(zhǎng)城底部軸線上均勻布設(shè)監(jiān)測(cè)點(diǎn)，設(shè)置x坐標(biāo)范圍為-50m～50m，測(cè)點(diǎn)間距為10m.據(jù)監(jiān)測(cè)點(diǎn)記錄的z方向位移數(shù)據(jù)表明，3種工法引起長(zhǎng)城底部的沉降值如圖6所示.

圖6 三種開挖工法下長(zhǎng)城沉降值與x坐標(biāo)軸的關(guān)系

由圖6可知，3種開挖工法下長(zhǎng)城沉降值與x軸坐標(biāo)關(guān)系趨勢(shì)基本相同，在拱頂中心位置正上方x=0處，長(zhǎng)城底部中心，測(cè)得該處的沉降值最大，沉降值最大為；當(dāng)隨著x軸的坐標(biāo)由長(zhǎng)城底部中心向兩側(cè)移動(dòng)，測(cè)得的長(zhǎng)城基礎(chǔ)的沉降值由最大沉降值開始逐漸減小，且以x=0直線為對(duì)稱軸，x軸左右兩側(cè)長(zhǎng)城底部沉降值呈軸對(duì)稱式分布，沉降曲線走向基本呈“v”字型.

3.2.2 拱頂監(jiān)測(cè)點(diǎn)沉降值

與地表變形值類似，在隧道拱頂沿y軸軸線均勻布設(shè)監(jiān)測(cè)點(diǎn)，控制坐標(biāo)范圍為0m～100m，布設(shè)間距為10m，測(cè)得在3種工況下對(duì)應(yīng)的隧道拱頂沉降值與y軸坐標(biāo)的關(guān)系，如圖7所示.

圖7 三種開挖工法下拱頂沉降值與y軸坐標(biāo)的關(guān)系

由圖7所示可知，采用CRD法開挖拱頂沉降值最小，為13.93mm；其次是CD法開挖的拱頂沉降值，為14.04mm；造成拱頂沉降值最大的是臺(tái)階法，拱頂沉降值為36.37mm.

3.3 應(yīng)力分析

隧道開挖完成后隧道上初期支護(hù)上關(guān)鍵點(diǎn)的最大、最小主應(yīng)力值如表3所示(符號(hào)為拉為正,壓為負(fù)).由表3可知，3種工法的拱頂拉應(yīng)力均為最大拉應(yīng)力處，拱肩次之，拱腳處最大主應(yīng)力為負(fù)值,即拱腳全部承受壓應(yīng)力.

3種工法開挖在最大跨處的壓應(yīng)力要比其余關(guān)鍵點(diǎn)壓應(yīng)力值都大，而拱頂位置處壓應(yīng)力次之.最小壓應(yīng)力CD法和CRD法位于拱底，臺(tái)階法最小壓應(yīng)力位于拱腳.分析可知結(jié)論，CRD法在各關(guān)鍵監(jiān)測(cè)點(diǎn)的拉應(yīng)力與壓應(yīng)力都要小于CD法、臺(tái)階法.故3種開挖工法下CRD法應(yīng)為最優(yōu)考慮.

表3 3種開挖工法下關(guān)鍵監(jiān)測(cè)點(diǎn)應(yīng)力值

4 結(jié) 論

通過使用有限元軟件FLAC3D，對(duì)CD法、CRD法和三臺(tái)階法3種開挖工法下隧道下穿古長(zhǎng)城的施工過程進(jìn)行數(shù)值模擬，分析不同工況下的長(zhǎng)城沉降值以及各監(jiān)測(cè)點(diǎn)沉降值與隧道位置的關(guān)系，得到以下結(jié)論：

1.使用CRD法進(jìn)行開挖，沉降變形多發(fā)生在拱頂中心靠右的位置.左半部分開挖完成后，施做襯砌支護(hù)封閉，此時(shí)拱頂處于懸空狀態(tài)且左側(cè)圍巖位移受限，因此右側(cè)未開挖完全的圍巖就會(huì)產(chǎn)生較大位移，右側(cè)開挖累計(jì)沉降大于左側(cè)，證明CRD法的開挖工序會(huì)對(duì)圍巖變形與位移造成一定影響.故在偏壓圍巖開挖前，應(yīng)探明圍巖情況，通過模擬計(jì)算去避免不利狀況的發(fā)生.

2.長(zhǎng)城軸線與隧道拱頂中心軸線的相交處為長(zhǎng)城沉降最大值處，且以拱頂中心軸線為對(duì)稱軸，長(zhǎng)城沉降值沿軸線兩側(cè)逐漸減小，呈軸對(duì)稱式分布.

3.CD法開挖的沉降與應(yīng)力值相比CRD法稍大，三臺(tái)階法最大.CRD法比三臺(tái)階法在拱肩豎向位移減少52.5%，最跨度處水平收斂減少66.6%，拱腳豎向位移減少52.9%，拱底隆起減少64.9%.根據(jù)國(guó)家文物局于2019年12月3日回復(fù)青海文物局《國(guó)家文物局辦公室關(guān)于國(guó)道227線涉及青海長(zhǎng)城、長(zhǎng)寧遺址和寺溝遺址保護(hù)區(qū)部分工程建設(shè)方案意見的函》提出的指導(dǎo)意見，長(zhǎng)城沉降值需盡可能減小，不能超過5mm的限制.使用CRD法進(jìn)行下穿施工引起的長(zhǎng)城沉降為4.98mm，能夠滿足控制古長(zhǎng)城沉降、保證結(jié)構(gòu)穩(wěn)定性的施工要求，故CRD法為軟弱圍巖地區(qū)隧道下穿古長(zhǎng)城綜合效益最高的施工方法.