數(shù)值模擬計算在淺埋暗挖通道設(shè)計中的應(yīng)用

楊克軍

（溫州設(shè)計集團(tuán)有限公司，浙江 溫州 325000）

0 引言

近年來隨著城市軌道交通的建設(shè)，采用淺埋暗挖法施工的通道日益增多。國內(nèi)外目前普遍采用的通道結(jié)構(gòu)設(shè)計模型可歸納為四類：以工程類比為依據(jù)的經(jīng)驗法；以測試為依據(jù)的實用法；作用-反作用模型稱為結(jié)構(gòu)力學(xué)法；連續(xù)介質(zhì)模型稱為連續(xù)介質(zhì)力學(xué)法[1]。而隨著有限元軟件的日益普及，數(shù)值模擬計算在通道結(jié)構(gòu)計算中的應(yīng)用也日益成熟。

1 通道設(shè)計概況

某軟土地區(qū)地下通道橫穿城市道路，道路寬度為30 m，通道呈一字形布置，通道開挖寬度7.6 m，通道頂覆土厚度為4.0 m。通道采用復(fù)合式襯砌，初期支護(hù)采用C20網(wǎng)噴混凝土+格柵鋼架+連接筋聯(lián)合支護(hù)，二次襯砌采用C30、S8現(xiàn)澆鋼筋混凝土。通道施工采用108大管棚超前支護(hù)、TSS輔助二重管無收縮漿全斷面預(yù)加固；分上下兩層共六塊，進(jìn)行開挖支護(hù)并設(shè)臨時仰拱和支撐，確保開挖安全；開挖臺階長度控制在3～5 m，二襯每次施工長度控制在6～9 m。

2 開挖方式數(shù)值模擬計算

2.1 有限元模型

采用ABAQUS通用計算軟件，土體采用DP模型，支撐、初襯、二襯采用彈性模型，建立有限元模型。

（1）支撐采用梁單元模擬，初襯、二襯采用殼單元模擬，土體采用8節(jié)點實體單元模擬。

（2）TSS注漿管對周邊土體的加固采用等效材料法近似處理。

（3）超前支護(hù)管棚作用提高管棚作用范圍內(nèi)的圍巖材料參數(shù)來近似模擬。

（4）地面車輛、人流及施工荷載按照20 kN/m2來考慮。

（5）開挖面通過施加荷載來平衡。

根據(jù)經(jīng)驗，計算范圍將開挖的影響區(qū)域包括進(jìn)去，取開挖洞室的3～5倍洞徑建立模型。具體如圖1～圖3所示。

圖1 有限元網(wǎng)格模型

圖2 斷面尺寸（單位：cm）

圖3 外圍TSS注漿加固

計算模型縱向取100 m，橫向按照最長通道31 m來考慮，深度為50 m。根據(jù)上述模擬，按照實際施工過程進(jìn)行模擬計算，以下分別按照三種開挖方式計算并分析地表沉降。

2.2 通道開挖數(shù)值模擬分析

通道開挖方式按照圖4的施工順序進(jìn)行模擬計算，即先開挖兩側(cè)，后開挖中間土體，如圖4所示。

圖4 分布開挖順序

計算結(jié)果如圖5～圖11所示。

圖5 開挖第一步豎向變形云圖

圖6 開挖第二步豎向變形云圖

圖7 開挖第三步豎向變形云圖

圖8 開挖第四步豎向變形云圖

圖9 開挖第五步豎向變形云圖

圖10 開挖第六步豎向變形云圖

圖11 地表沉降曲線

從圖中可以看出，第一、二步下，先開挖左側(cè)，左側(cè)地面沉降明顯大于右側(cè)，第三、四步下右側(cè)沉降增大，兩側(cè)基本對稱；每個開挖步下沉降都在增大，地面最大沉降最終為40 m m。第四步下地面沉降基本達(dá)到最大，第五、六步下沉降基本不再增大。前面三步開挖卸載引起的地面差異沉降值明顯大于后面兩步的。

2.3 計算結(jié)論

（1）考慮注漿加固效果，開挖引起的地表最大沉降為4.2 cm。

（2）距開挖面10 m以下的土層基本穩(wěn)定。

（3）分步開挖引起的地表沉降變形為一遞增過程，每步開挖必須嚴(yán)格控制施工工藝，尤其是第一至第三步開挖。

3 拆撐穩(wěn)定性數(shù)值模擬分析

通道開挖完成后，將進(jìn)行二次襯砌澆筑，此時需將臨時支撐逐步拆除。臨時支撐的拆除打破了原有地層、結(jié)構(gòu)體系的平衡，造成應(yīng)力的重分布。對拆撐過程進(jìn)行模擬分析計算的目的在于了解這一過程初支結(jié)構(gòu)內(nèi)力的變化，驗證拆撐順序的合理性，在保證初支結(jié)構(gòu)穩(wěn)定的前提下，為二次襯砌的施作創(chuàng)造條件。

按照上述模型進(jìn)行拆撐計算分析，支撐及初襯網(wǎng)格如圖12所示，當(dāng)通道開挖完成初襯支護(hù)做好后，通道周圍土體出現(xiàn)塑性區(qū)位置如圖13所示，可見在對通道周圍土體加固時應(yīng)重點放在這些部位。

圖12 支撐網(wǎng)格圖

圖13 開挖后塑性區(qū)云圖

3.1 拆撐過程內(nèi)力分析

首先取一品鋼架和支撐進(jìn)行二維分析，先拆除下導(dǎo)支撐，后拆除上部支撐，圍巖壓力按照等效荷載來考慮，如圖14～圖19所示。

圖14 未拆撐前N

圖15 未拆撐前M

圖16 拆除下導(dǎo)支撐時N

圖17 拆除下導(dǎo)支撐時M

圖18 拆除全部支撐時N

圖19 拆除全部支撐時M

從圖中可以看出，在拆撐過程中，初襯軸力在增大，下部拱軸力減??；上部拱及上部拐角處彎矩在增大，所以設(shè)計時應(yīng)考慮拆撐給結(jié)構(gòu)造成的不利影響。

3.2 拆撐過程三維分析

主要用于模擬施工過程拆撐縱向?qū)νǖ澜Y(jié)構(gòu)的影響，該次計算分析三種方案，一次拆撐長度分別為6 m、10 m、16 m。建立的模型如圖12所示，計算結(jié)果如圖20～圖31所示。

圖20 拆除6 m范圍內(nèi)下導(dǎo)支撐截面軸力云圖

圖21 拆除6 m范圍內(nèi)下導(dǎo)支撐截面彎矩云圖

圖22 拆除6 m范圍內(nèi)全部支撐截面軸力云圖

圖23 拆除6 m范圍內(nèi)全部支撐截面彎矩云圖

圖24 拆除10m范圍內(nèi)下導(dǎo)支撐截面軸力云圖

圖25 拆除10 m范圍內(nèi)下導(dǎo)支撐截面彎矩云圖

圖26 拆除10 m范圍內(nèi)全部支撐截面軸力云圖

圖27 拆除10 m范圍內(nèi)全部支撐截面彎矩云圖

圖28 拆除16 m范圍內(nèi)下導(dǎo)支撐截面軸力云圖

圖29 拆除16 m范圍內(nèi)下導(dǎo)支撐截面彎矩云圖

圖30 拆除16 m范圍內(nèi)全部支撐截面軸力云圖

圖31 拆除16 m范圍內(nèi)全部支撐截面彎矩云圖

3.3 計算結(jié)論

（1）拆除下導(dǎo)支撐對結(jié)構(gòu)體系的影響最大，拱頂彎矩增大120%，拆除時慎重考慮相應(yīng)的保護(hù)性措施，確保施作安全。

（2）拆除完畢后，結(jié)構(gòu)初支最大彎矩出現(xiàn)在拱部與邊墻連接處，在結(jié)構(gòu)配筋時應(yīng)考慮這一因素。

（3）拆撐長度為16m時，拱部中間位置最大應(yīng)力增大35%，使初支應(yīng)力達(dá)到1.3 M Pa。

（4）拆撐時引起結(jié)構(gòu)變形值，6 m為5.5 m m，9 m為8.5 m m，12 m 為11.6 m m。

綜上所述，拆撐長度可定為6～9m，同時在拆撐過程中應(yīng)采取提前主動換撐、加強(qiáng)通道結(jié)構(gòu)變形與應(yīng)力的監(jiān)測。

4 結(jié)語

通過上述案例分析，并通過項目實踐檢驗，連續(xù)介質(zhì)力學(xué)法即數(shù)值模擬計算完全能模擬常規(guī)地質(zhì)條件下暗挖通道施工過程中的不同施工方式、工序帶來的襯砌受力變化情況、路面沉降變化及相關(guān)建筑物、構(gòu)筑物受影響狀況。根據(jù)其模擬結(jié)果進(jìn)行分析計算，能更加精確地進(jìn)行設(shè)計和對施工過程事故進(jìn)行防范。