基于FLAC3D的火山灰地質(zhì)隧洞開挖支護(hù)數(shù)值模擬

劉亮亮

(華北水利水電大學(xué) 河南 鄭州 450045)

基于FLAC3D的火山灰地質(zhì)隧洞開挖支護(hù)數(shù)值模擬

劉亮亮

(華北水利水電大學(xué) 河南 鄭州 450045)

該文從一工程實(shí)例出發(fā)，利用在隧洞數(shù)值模擬領(lǐng)域具有領(lǐng)先優(yōu)勢(shì)和良好口碑的FLAC3D軟件進(jìn)行數(shù)值模擬，對(duì)不同工況下的隧洞開挖進(jìn)行模擬計(jì)算，研究了火山灰地質(zhì)隧洞在開挖過(guò)程中的應(yīng)力應(yīng)變變化情況以及隧洞支護(hù)的效果，分析結(jié)果為隧洞施工控制提供了理論支持；同時(shí)通過(guò)控制變量法，得出剪切強(qiáng)度參數(shù)控制著數(shù)值模擬計(jì)算能否收斂，也控制著塑性區(qū)的擴(kuò)展深度，在計(jì)算可以收斂的前提下，洞室位移主要由彈性模量控制。

隧洞；火山灰；開挖支護(hù)；FLAC3D；數(shù)值模擬；控制變量法

工程中，我們要了解輸水隧洞支護(hù)的效果和改進(jìn)支護(hù)的措施，就需要對(duì)開挖支護(hù)過(guò)程和完成支護(hù)后的穩(wěn)定性進(jìn)行研究。開挖過(guò)程中引起的圍巖荷載的變化更要重視。FLAC3D模擬的系統(tǒng)是靜態(tài)的，但它使用動(dòng)態(tài)方程，不存在模擬物理上不穩(wěn)定性的數(shù)值障礙，給我們的研究帶來(lái)了很大的便利[1]。該文用FLAC3D建立二維平面應(yīng)變模型，對(duì)隧洞開挖和支護(hù)的各個(gè)工況進(jìn)行模擬計(jì)算分析。

一、基本概況

某輸水隧洞地層為第四系火山崩塌堆積碎片、火山灰為主間夾河流沖洪積層，組成為火山灰、火山碎屑、礫砂層、崩積火山巖塊，火山碎屑堆積層厚超過(guò)130 m，含礫砂層屬于疏松、稍密狀。根據(jù)地質(zhì)條件，采用鉆爆法開挖，支護(hù)采用錨桿和鋼筋網(wǎng)噴混凝土。研究段輸水隧洞總長(zhǎng)525 m，隧洞東、西走向，縱坡為0.513%，埋深300～700 m，隧洞為無(wú)壓洞設(shè)計(jì)，成洞洞徑9.61 m，開挖洞徑9.81～10.21 m。

二、基本模型

按平面應(yīng)變問(wèn)題解決，取計(jì)算坐標(biāo)系為oxy，x軸為水平方向，指向右側(cè)為正；y軸為垂直方向，指向上方為正。計(jì)算范圍為-240 m≤x≤240 m，-40 m≤y≤240 m。共劃分6318個(gè)網(wǎng)格，在洞室附近劃分了1m×1m的精細(xì)網(wǎng)格，其他區(qū)域的網(wǎng)格為4m×4m。模型兩側(cè)水平移動(dòng)即x方向方位移限制為零，模型的底部垂直和水平移動(dòng)即x和y方向位移都限制為零。

三、開挖支護(hù)方式及計(jì)算工況

(一)開挖支護(hù)方式

隧洞開挖分兩步，先開挖上部區(qū)域，再開挖下部區(qū)域。根據(jù)需要，本文從開挖左上區(qū)域和全部區(qū)域并分別進(jìn)行一次支護(hù)進(jìn)行模擬研究。隧洞斷面支護(hù)斷面圖見圖1。

圖1 支護(hù)橫斷面圖

(二)計(jì)算工況

為了研究隧洞開挖支護(hù)的過(guò)程及計(jì)算參數(shù)對(duì)洞室位移的影響，對(duì)左上區(qū)域開挖和全部開挖從以下幾個(gè)工況著手分析：①原參數(shù)開挖無(wú)支護(hù)②原參數(shù)開挖有支護(hù)③增強(qiáng)剪切強(qiáng)度參數(shù)開挖有支護(hù)④增強(qiáng)變形和剪切強(qiáng)度參數(shù)開挖有支護(hù)。對(duì)所有工況進(jìn)行模擬開挖并分析其位移場(chǎng)和塑性區(qū)。

四、模擬計(jì)算及分析結(jié)果

(一)初始地應(yīng)力場(chǎng)

由計(jì)算分析得出開挖區(qū)域的垂直地應(yīng)力在4～5 MPa之間，水平地應(yīng)力在2～2.5 MPa之間。地表沉降達(dá)到了20 m。水平方向位移幾乎為0。通過(guò)分析塑性區(qū)，初始狀態(tài)下沒(méi)有發(fā)生屈服。

(二)原參數(shù)開挖

表1 計(jì)算采用的巖體力學(xué)參數(shù)

1.無(wú)支護(hù)工況

無(wú)支護(hù)時(shí)，開挖左上區(qū)域后的垂直方向最大位移分別為15 m(頂拱下降)和5 m(底板隆起)，見圖2。洞室上部塑性區(qū)深度達(dá)到了10 m，洞室下部塑性區(qū)深度5 m，見圖3。繼續(xù)開挖計(jì)算不能收斂，塑性區(qū)范圍無(wú)限擴(kuò)展，說(shuō)明該工況下的洞室不能成形，隨開挖隨塌。

圖2 無(wú)支護(hù)時(shí)垂直方向位移場(chǎng)

圖3 無(wú)支護(hù)時(shí)塑性區(qū)

2.有支護(hù)工況

有支護(hù)時(shí)，開挖左上區(qū)域后的垂直方向最大位移分別為2 m(頂拱下降)和5 m(底板隆起)，見圖4。洞室上部塑性區(qū)深度達(dá)到了5 m，洞室下部塑性區(qū)深度5 m，見圖5。繼續(xù)開挖計(jì)算不能收斂，塑性區(qū)范圍無(wú)限擴(kuò)展，這說(shuō)明該工況下的洞室不能成形，隨開挖隨塌。

圖4 有支護(hù)時(shí)垂直方向位移場(chǎng)

圖5 有支護(hù)時(shí)塑性區(qū)

綜合比較發(fā)現(xiàn)，有支護(hù)時(shí)，位移量值和塑性區(qū)深度都減小很多。表明了支護(hù)作用比較明顯。

(三)增強(qiáng)剪切強(qiáng)度參數(shù)開挖

表2 計(jì)算采用的巖體力學(xué)參數(shù)

有支護(hù)時(shí)，開挖左上區(qū)域后的垂直方向最大位移分別為0.8 m(頂拱下降)和2 m(底板隆起)，見圖6。塑性區(qū)范圍大幅度減小，頂部只有零星塑性區(qū)(不超過(guò)1 m深)，其他區(qū)域塑性區(qū)深度不超過(guò)5 m，見圖7。

圖6 有支護(hù)時(shí)垂直方向位移場(chǎng)

圖7 有支護(hù)時(shí)塑性區(qū)

繼續(xù)開挖計(jì)算可以收斂，但是位移比較大，垂直方向最大位移分別為3.6 m(頂拱下降)和2.8 m(底板隆起)，見圖8。開挖完畢后，塑性區(qū)深度一般為5 m，見圖9。

圖8 開挖完畢后垂直方向位移

圖9 開挖完畢后塑性區(qū)

與原參數(shù)相比，增強(qiáng)剪切強(qiáng)度參數(shù)后，計(jì)算可以收斂，并且位移大幅度減小。

(四)增強(qiáng)變形和剪切強(qiáng)度參數(shù)開挖

表3 計(jì)算采用的巖體力學(xué)參數(shù)

開挖左上區(qū)域有支護(hù)時(shí)，開挖左上區(qū)域后的垂直方向最大位移分別為8 cm(頂拱下降)和16 cm(底板隆起)。塑性區(qū)范圍與只增強(qiáng)剪切強(qiáng)度參數(shù)時(shí)的結(jié)果一致，頂部只有零星塑性區(qū)(不超過(guò)1 m深)，其他區(qū)域塑性區(qū)深度不超過(guò)5 m。這也說(shuō)明增強(qiáng)彈性模量的作用，大幅度減小了位移，但是不影響塑性區(qū)的范圍。

開挖完畢后繼續(xù)開挖計(jì)算可以收斂，但是位移比較大，垂直方向最大位移分別為32cm(頂拱下降)和24cm(底板隆起)。塑性區(qū)范圍與只增強(qiáng)剪切強(qiáng)度參數(shù)時(shí)的結(jié)果一致，開挖完畢后，塑性區(qū)深度一般為5m。這也說(shuō)明增強(qiáng)彈性模量的作用，大幅度減小了位移，但是不影響塑性區(qū)的范圍。

(五)綜合比較

表4 開挖左上區(qū)域后位移對(duì)比表

表5 全部開挖后位移對(duì)比表

五、結(jié)論

a)在原參數(shù)的情況下，無(wú)論是否有支護(hù)，繼續(xù)開挖都不能收斂，塑性區(qū)無(wú)限擴(kuò)展。這說(shuō)明了開挖過(guò)程中洞室很難成形，隨開隨塌。

b)在增強(qiáng)剪切強(qiáng)度參數(shù)以后，計(jì)算可以收斂，但是位移依然偏大。這說(shuō)明剪切強(qiáng)度參數(shù)控制著計(jì)算能否收斂，也控制著塑性區(qū)的擴(kuò)展深度。

c)在同時(shí)增強(qiáng)變形參數(shù)和剪切強(qiáng)度參數(shù)后，計(jì)算不僅可以收斂，而且位移大幅度減小，但是塑性區(qū)深度沒(méi)有改變。這說(shuō)明在計(jì)算可以收斂的前提下，洞室位移主要由彈性模量控制。

d)由于地應(yīng)力在4～5 MPa之間，地應(yīng)力量值并不大，所以支護(hù)的作用效果會(huì)比較明顯。

e)綜合以上分析，認(rèn)為有必要對(duì)開挖影響區(qū)域的巖體進(jìn)行注漿改性。

[1]張宇娜．某水電站引水隧洞開挖支護(hù)數(shù)值模擬研究[J]．南水北調(diào)與水利科技，2011，(02)：101-104+109．

[2]劉英，于立宏．Mohr-Coulomb屈服準(zhǔn)則在巖土工程中的應(yīng)用[J]．世界地質(zhì)，2010，(04)：633-639．

[3]邵建國(guó)．FLAC3D在隧道施工中的數(shù)值模擬分析[J]．甘肅科技，2014，(23)：130-133．

[4]尤帆帆，邱曉莉．基于FLAC3D的隧道分臺(tái)階開挖數(shù)值模擬與分析[J]．山西建筑，2014，(28)：189-190．

[5]金長(zhǎng)文，韓鵬．CCS水電站火山灰地層大斷面隧洞開挖施工技術(shù)[J]．云南水力發(fā)電，2014，(05)：66+124．

[6]段錫健．隧洞掘進(jìn)與支護(hù)過(guò)程的數(shù)值模擬[A]．中國(guó)礦業(yè)科技文匯，2013：4

劉亮亮，男，碩士生，研究方向：水工結(jié)構(gòu)模型試驗(yàn)與數(shù)值分析。