宋振余 毛安琪 易山山
(1.廣東省南粵交通揭惠高速公路管理中心 揭陽(yáng) 515325; 2.中國(guó)建筑科學(xué)研究院 北京 100013; 3.華中科技大學(xué)土木工程與力學(xué)學(xué)院 武漢 430074)
弱質(zhì)巖體在大自然的作用下發(fā)育著大量的裂隙、孔隙和節(jié)理面,由于這些裂隙的存在,為地下水活動(dòng)提供了相應(yīng)的空間,同時(shí)對(duì)巖體的完整性造成了一定的影響[1-2]。地下水通過(guò)滲流效應(yīng)改變了巖體的天然結(jié)構(gòu),從而影響到應(yīng)力場(chǎng)的重分布,導(dǎo)致隧道圍巖產(chǎn)生變形。反過(guò)來(lái)說(shuō),巖體中應(yīng)力場(chǎng)和位移場(chǎng)也會(huì)導(dǎo)致巖體裂隙的改變,進(jìn)而影響到其滲透性能。這種相互作用反復(fù)耦合直到達(dá)到某種平衡狀態(tài),稱(chēng)之為流-固耦合效應(yīng)[3-4]。地下水在下滲過(guò)程中,易引起隧道施工現(xiàn)場(chǎng)塌方和大變形等事權(quán)和病害,對(duì)隧道施工和運(yùn)營(yíng)造成不利影響,因此,對(duì)弱質(zhì)圍巖隧道在地下水滲流效應(yīng)下的開(kāi)挖變形分析是亟待解決的問(wèn)題。目前,對(duì)于隧道圍巖應(yīng)力場(chǎng)和滲流場(chǎng)耦合作用下的變形分析較少,開(kāi)挖過(guò)程中的發(fā)展歷程也很少研究[5-6]。
山嶺隧道一般采用新奧法施工,主要靠隧道施工監(jiān)控量測(cè)及時(shí)掌握圍巖變形的情況以判斷圍巖穩(wěn)定性[7]。對(duì)于高水位軟弱圍巖隧道,在高地應(yīng)力、高水壓圍巖應(yīng)力場(chǎng)與滲流場(chǎng)耦合作用下,隧道開(kāi)挖圍巖變形值很大且變形速度很快。針對(duì)高水壓斷層破碎帶和節(jié)理發(fā)育密集帶的止水注漿及穿越軟弱、不良地質(zhì)段的加固注漿對(duì)隧道的施工安全、施工進(jìn)度及運(yùn)營(yíng)后的止排水等都有顯著效果[8-9],因此,結(jié)合工程實(shí)際開(kāi)展高水位下破碎帶圍巖開(kāi)挖變形規(guī)律的研究非常有必要。
廣東省揭陽(yáng)至惠來(lái)高速公路小北山一號(hào)隧道位于丘陵區(qū)。隧道按分離式布設(shè),左線(xiàn)隧道里程ZK14+390-ZK17+390,最大埋深280 m。右線(xiàn)隧道里程K14+380-K17+388,最大埋深270 m。
隧址區(qū)基底主要為燕山期花崗巖,局部見(jiàn)輝綠巖巖脈,覆蓋層由黏土、全~強(qiáng)風(fēng)巖組成,基巖由中~微風(fēng)化巖組成。地下水以大氣降水和山谷匯水下滲補(bǔ)給為主,排泄方式則以蒸發(fā)和側(cè)向徑流排泄為主。隧道在ZK16+450-ZK16+600(K16+400-K16+580)段存在F3區(qū)域大斷裂,該斷層處圍巖破碎,裂隙發(fā)育,工程地質(zhì)復(fù)雜,屬于IV,V級(jí)圍巖,且其上方約150 m處存在龍?zhí)斗逅畮?kù),隧道開(kāi)挖過(guò)程中改變了天然地下水的補(bǔ)徑排條件,隧道成為新的局部排泄基準(zhǔn),易發(fā)生塌方、涌水、水資源流失、生態(tài)環(huán)境破壞等不良情況。
小北山一號(hào)隧道軟巖施工現(xiàn)場(chǎng)情況見(jiàn)圖1。
圖1 隧道軟巖施工現(xiàn)場(chǎng)
根據(jù)設(shè)計(jì)要求,依據(jù)小北山隧道工程的實(shí)際情況,確定了監(jiān)測(cè)、監(jiān)控項(xiàng)目的主要范圍,并在實(shí)際的監(jiān)測(cè)過(guò)程中予以了相應(yīng)的調(diào)整和改進(jìn),本文只列出其中的必測(cè)項(xiàng)目,見(jiàn)表1。
表1 施工監(jiān)測(cè)必測(cè)項(xiàng)目
以上監(jiān)測(cè)項(xiàng)目中,最主要的是拱頂下沉及水平收斂量測(cè),也是本文分析中需要的監(jiān)測(cè)數(shù)據(jù),因此,僅對(duì)拱頂下沉及水平收斂監(jiān)測(cè)方法進(jìn)行介紹。
1) 拱頂下沉量測(cè)。拱頂下沉值主要用于確認(rèn)圍巖穩(wěn)定、預(yù)報(bào)拱頂崩塌。通過(guò)使用水準(zhǔn)儀及銦鋼尺測(cè)量觀測(cè)點(diǎn)與基準(zhǔn)點(diǎn)的相對(duì)高差變化量得出拱頂下沉量和下沉速度,拱頂沉降量測(cè)與周邊位移量測(cè)布置在同一斷面,測(cè)點(diǎn)布置如圖2所示。斷面布置:V級(jí)圍巖1個(gè)斷面/10 m,VI級(jí)圍巖1個(gè)斷面/20 m,III級(jí)圍巖1個(gè)面/30 m,II級(jí)圍巖1個(gè)斷面/50 m。量測(cè)頻率如表2所示。
表2 現(xiàn)場(chǎng)監(jiān)測(cè)數(shù)據(jù)采集頻率表
2) 水平收斂量測(cè)。量測(cè)獲取隧道周邊位移量,以此判斷隧道結(jié)構(gòu)的穩(wěn)定性,確定二次襯砌施作的合理時(shí)機(jī)。用收斂計(jì)量測(cè)凈空收斂位移,2次測(cè)量之差即為該周邊2點(diǎn)在該時(shí)間間隔內(nèi)收斂值,根據(jù)收斂值判定隧道圍巖及支護(hù)結(jié)構(gòu)的變形情況。特殊地段可采用非接觸量測(cè)方法協(xié)助進(jìn)行量測(cè),即在隧道周邊洞壁粘貼反射膜片,利用全站儀進(jìn)行周邊位移觀測(cè)。斷面布置為V級(jí)圍巖1個(gè)斷面/10 m,VI級(jí)圍巖1個(gè)斷面/20 m,III級(jí)圍巖1個(gè)斷面/30 m,II級(jí)圍巖1個(gè)斷面/50 m。根據(jù)斷面開(kāi)挖方式,采取相應(yīng)的測(cè)點(diǎn)布置,如圖2所示。量測(cè)頻率如表2所示。
圖2 水平收斂、拱頂下沉測(cè)線(xiàn)布置圖
現(xiàn)場(chǎng)開(kāi)挖時(shí),斷裂破碎帶發(fā)生在K16+442-K16+457段,開(kāi)挖進(jìn)度為2 m/d。為防止圍巖發(fā)生坍塌,對(duì)此段進(jìn)行注漿加固并加強(qiáng)對(duì)此段的施工監(jiān)控,此段的監(jiān)控?cái)嗝娴牟荚O(shè)距離由10 m減為5 m。斷面1~5與現(xiàn)場(chǎng)監(jiān)測(cè)時(shí)的斷面K16+440,K16+445,K16+450,K16+455,K16+460基本對(duì)應(yīng),隧道監(jiān)測(cè)斷面位置如圖3所示。
圖3 監(jiān)測(cè)斷面位置圖
K16+440,K16+445,K16+450,K16+455,K16+460斷面的拱頂下沉和水平收斂監(jiān)測(cè)結(jié)果如圖4所示。
圖4所示的監(jiān)控結(jié)果表明:
1) 破碎帶圍巖的開(kāi)挖變形規(guī)律(拱頂下沉及水平收斂)與普通等級(jí)圍巖基本一致,均表現(xiàn)在掌子面附近時(shí)變形速度較快,遠(yuǎn)離掌子面后變形速度越來(lái)越小。
圖4 各斷面的監(jiān)測(cè)位移場(chǎng)曲線(xiàn)圖
2) 水平收斂測(cè)線(xiàn)a比測(cè)線(xiàn)b收斂值大。因?yàn)闇y(cè)線(xiàn)a在上臺(tái)階開(kāi)挖后布設(shè),測(cè)線(xiàn)b在下臺(tái)階開(kāi)挖后布設(shè),在上臺(tái)階開(kāi)挖后,位移已得到一定程度的釋放。
3) 拱頂下沉及水平收斂位移最大值均發(fā)生在K16+450處,且沿其兩側(cè)方向斷面的位移值逐漸減小。
運(yùn)用midas GTS NX軟件,采用參數(shù)弱化法對(duì)小北山一號(hào)隧道K16+400-K16+520軟弱破碎帶段進(jìn)行數(shù)值模擬,破碎帶處于V級(jí)圍巖,巖石破碎帶斷層寬度約12 m,傾角約為50°,隧道埋深取150 m,每次進(jìn)尺的長(zhǎng)度2 m,模擬開(kāi)挖30步,模型的計(jì)算范圍為120 m×120 m×210 m,整個(gè)數(shù)學(xué)模型共劃分為97 679個(gè)實(shí)體單元,4 480 個(gè)板單元,17 469個(gè)節(jié)點(diǎn)。據(jù)地質(zhì)勘探資料確定巖體的物理力學(xué)參數(shù),采用等效連續(xù)介質(zhì)模型,D-P彈塑性屈服準(zhǔn)則。巖體采用實(shí)體單元模擬,支護(hù)結(jié)構(gòu)用板單元模擬。初噴在圍巖開(kāi)挖后跟進(jìn)施作,二襯只在開(kāi)挖前20 m施作,只作為開(kāi)挖前支護(hù)條件,不隨開(kāi)挖的進(jìn)行而跟進(jìn)施作,開(kāi)挖網(wǎng)格組分布如圖5所示。
圖5 開(kāi)挖網(wǎng)格組分布圖
隧道埋深較大,隧道在垂直方向的初始地應(yīng)力取巖體的自重,水平方向的初始地應(yīng)力按側(cè)壓力系數(shù)λ=1考慮。假定隧道開(kāi)挖前,圍巖處于飽和狀態(tài)。開(kāi)挖后,地下水在隧道開(kāi)挖區(qū)域的邊界上為自由透水邊界,滲流水壓力為零。采用超前帷幕全斷面注漿加固,注漿加固效果通過(guò)提升注漿區(qū)圍巖的物理力學(xué)參數(shù)模擬,計(jì)算參數(shù)見(jiàn)表3。
表3 計(jì)算參數(shù)表
在隧道施工中,二襯的施作發(fā)生在圍巖初支作用下變形穩(wěn)定后,另外,初噴混凝土往往伴隨著圍巖的變形開(kāi)裂,其實(shí)際堵水作用十分有限,因此,在討論注漿加固對(duì)圍巖開(kāi)挖變形的影響時(shí),往往不考慮二襯及初支的堵水作用,表中設(shè)初支及二襯的滲透系數(shù)為無(wú)窮大。
不考慮開(kāi)挖滲流影響,單獨(dú)討論破碎帶對(duì)圍巖開(kāi)挖變形的影響。由于隧道與破碎帶斜交,隧道研究斷面位置如圖6所示,因此,可分5種不同的情況分析隧道開(kāi)挖破碎帶對(duì)圍巖變形的影響:
①隧道斷面拱底存在巖石破碎帶(斷面1)②隧道斷面下臺(tái)階為巖石破碎帶(斷面2);③整個(gè)隧道斷面都在巖石破碎帶內(nèi)(斷面3)④隧道斷面上臺(tái)階為巖石破碎帶(斷面4);⑤隧道斷面拱頂存在巖石破碎帶(斷面5)。
運(yùn)用有限元模擬后,得到在不考慮滲流情況下,含破碎帶圍巖開(kāi)挖后的總位移、豎向位移,如圖6、圖7所示。
圖6 破碎帶開(kāi)挖后總位移云圖
圖7 破碎帶開(kāi)挖后豎向位移云圖
圖6、圖7所示計(jì)算結(jié)果表明,斷層破碎帶處圍巖產(chǎn)生的位移量明顯大于兩側(cè)圍巖,這是因?yàn)閿鄬悠扑閹巼鷰r與兩側(cè)巖體在物理力學(xué)特性上差異顯著,所以,應(yīng)特別關(guān)注此段的施工過(guò)程,適當(dāng)提高支護(hù)參數(shù)。
提取含破碎帶圍巖斷面1~5并對(duì)比完整V級(jí)圍巖的下拱頂下沉、拱底隆起及水平收斂的計(jì)算結(jié)果并對(duì)比完整V級(jí)圍巖計(jì)算結(jié)果,如圖8所示。
由圖8可見(jiàn):
1) 含破碎帶時(shí),各斷面在開(kāi)挖時(shí)的拱頂下沉、拱底隆起和水平收斂值均比完整圍巖時(shí)大且變形速度更快。
圖8 斷面1~5位移場(chǎng)曲線(xiàn)對(duì)比圖
2) 斜向破碎帶對(duì)圍巖拱頂下沉、拱底隆起和水平收斂變形影響不同。
①對(duì)拱頂位移的影響:各斷面拱頂下沉值從小到大依次為1,2,3,5,4,破碎帶對(duì)圍巖拱頂位移的影響在斷面2開(kāi)始逐漸顯現(xiàn),在斷面4達(dá)到最大下沉值,且斷面5的拱頂下沉值也較大,可見(jiàn)破碎帶對(duì)圍巖拱頂位移的影響在斷面5后仍會(huì)延續(xù)一段距離。
②對(duì)拱底位移的影響:各斷面拱底隆起值由小到大依次為5,4,3,2,1,圍巖拱底位移在斷面1處就達(dá)到最大隆起值,破碎帶對(duì)圍巖拱底位移的影響在斷面1前已有較大一段距離;在斷面4,5的拱底隆起值已基本接近完整圍巖,可見(jiàn)破碎帶在此處對(duì)圍巖拱頂位移的影響已基本結(jié)束。
③對(duì)水平收斂的影響:各斷面水平收斂值小到大依次是5,1,4,2,3,圍巖水平收斂在斷面3處達(dá)到最大下沉值,沿著左、右量測(cè)逐漸減小,斷面5處水平收斂值已基本接近完整圍巖,斷面1處水平收斂值較完整圍巖大,說(shuō)明破碎帶對(duì)斷面1前圍巖水平收斂仍有一段影響距離,對(duì)斷面5以后基本無(wú)影響。
由此看來(lái),拱頂下沉、拱底隆起和水平收斂值并不是所有隧道斷面都在破碎帶時(shí)達(dá)到最大值。拱頂下沉最大值發(fā)生在斷面4,5主要是由于此處開(kāi)挖部分為V級(jí)圍巖,而隧道拱頂為破碎圍巖,由于V級(jí)圍巖重度、彈性模量等均比破碎帶大,因此開(kāi)挖回彈量比全破碎帶大,拱頂下沉值更大。同理,拱底最大值發(fā)生在斷面1主要是由于此處開(kāi)挖部分為V級(jí)圍巖,而此處隧道底部為破碎圍巖,因此開(kāi)挖回彈量比全破碎帶更大,拱底隆起更大。水平收斂值在斷面3處最大,此處整個(gè)隧道處在破碎帶中,是由于此處開(kāi)挖后隧道周?chē)鶠槠扑檐泿r,對(duì)于水平收斂的限制大大減小,因此,此處水平收斂值最大。
綜上所述,當(dāng)圍巖中存在破碎帶,尤其是斜向破碎帶時(shí),隧道開(kāi)挖時(shí)不僅要注意破碎帶區(qū)的圍巖變形,也要注意破碎帶前后一定范圍內(nèi)的圍巖變形,在采取注漿加固時(shí)應(yīng)適當(dāng)擴(kuò)大注漿范圍。
4.選擇合適的教學(xué)技術(shù)。利用多媒體課件和網(wǎng)絡(luò)信息資源,做到“互聯(lián)網(wǎng)+課堂教學(xué)”,可以讓課堂教學(xué)全校同步,做到“班班通”,甚至“校校通”。
考慮開(kāi)挖滲流影響,采用超前帷幕全斷面注漿加固,注漿加固范圍為5 m,注漿圈滲透系數(shù)取2.0×10-5cm/s,其他參數(shù)同表3。
3.3.1滲流場(chǎng)分析
對(duì)破碎帶圍巖注漿前后分別進(jìn)行滲流分析,得到注漿加固前、后圍巖孔隙水壓力分布圖,如圖9所示。
圖9 注漿前后破碎帶圍巖孔隙水壓力云圖
由圖9可見(jiàn):
1) 注漿加固前,由于破碎軟巖滲透系數(shù)大于其兩側(cè)圍巖滲透系數(shù),破碎軟巖處水壓力小于兩側(cè)圍巖,因此破碎帶處孔隙水壓力比兩側(cè)圍巖孔隙水壓力小。
2) 注漿加固后,整個(gè)隧道周邊圍巖滲透系數(shù)減小,注漿圈承擔(dān)了一部分水壓力,使得整個(gè)隧道周邊孔隙水壓力增大,減小開(kāi)挖滲流的影響,同時(shí)整個(gè)隧道周邊孔隙水壓力分布更加均勻,幾乎可以忽略破碎帶的影響。
進(jìn)行應(yīng)力-滲流耦合分析,得到破碎帶圍巖在滲流作用下注漿加固前后位移云圖,如圖10所示。
圖10 加固前后破碎帶圍巖開(kāi)挖總位移云圖
提取破碎帶處不考慮滲流情況、完整圍巖150 m水頭差、破碎帶150 m水頭差,以及破碎帶注漿后的拱頂下沉、拱底隆起及水平收斂的計(jì)算結(jié)果,如圖11所示。
由圖11可見(jiàn):
1) 相較于不考慮滲流情況,在考慮滲流時(shí),高水位下破碎帶與其兩側(cè)圍巖變形差異更大,尤其是破碎帶拱底隆起值高達(dá)300 mm,此時(shí),拱底基本已產(chǎn)生隆起破壞,這也是加固前含破碎帶圍巖開(kāi)挖總位移、豎向位移云圖色彩范圍聚集在破碎帶拱底處的原因,如圖11a)所示,因此必須對(duì)圍巖進(jìn)行注漿加固處理。
圖11 不同情況位移場(chǎng)曲線(xiàn)對(duì)比圖
2) 破碎帶在注漿加固后其拱頂下沉、拱底隆起及水平收斂值均大幅減小,且最終變形值均比完整圍巖還小,變形過(guò)程較完整圍巖更加平緩。這主要是由于注漿加固后,隧道周邊圍巖滲透系數(shù)減小,整個(gè)隧道周邊孔隙水壓力分布更加均勻。另外,注漿圈承擔(dān)了一部分水壓力,減小開(kāi)挖滲流的影響,加上注漿圈材料參數(shù)的提高,幾乎可以忽略破碎帶的影響,使得圍巖變形值大幅減小。
提取注漿加固后的斷面1~5開(kāi)挖后的圍巖變形數(shù)值模擬結(jié)果,并對(duì)比K16+440,K16+445,K16+450,K16+0455,K16+460的拱頂下沉和水平收斂的監(jiān)測(cè)結(jié)果,如圖12所示。
圖12 實(shí)測(cè)結(jié)果與模擬結(jié)果對(duì)比圖
由圖12可見(jiàn):
1) 各斷面開(kāi)挖后的實(shí)測(cè)變形與數(shù)值模擬變形規(guī)律基本一致,但拱頂下沉、水平收斂實(shí)測(cè)結(jié)果均比數(shù)值模擬結(jié)果偏大,且不同位置偏大程度不同。
經(jīng)分析認(rèn)為,由于數(shù)值模擬中圍巖物理力學(xué)參數(shù)取值的誤差,加上現(xiàn)場(chǎng)施工中,支護(hù)施作的實(shí)際效果不可能達(dá)到理想設(shè)計(jì)效果,破碎帶區(qū)的注漿加固后的圍巖參數(shù)應(yīng)小于其兩側(cè)圍巖參數(shù),而不是數(shù)值模擬中取用的相同值,這是造成實(shí)測(cè)結(jié)果偏大的原因,因此實(shí)測(cè)拱頂下沉及水平收斂最大值均發(fā)生在K16+450處。不同位置偏大程度不同是由于斜向破碎帶在不同斷面破碎圍巖分布不同,例如,K16+450斷面實(shí)測(cè)拱頂下沉值約為模擬值的1.8倍,而K16+460斷面實(shí)測(cè)拱頂下沉值僅為模擬值的1.2倍,這是由于K16+450斷面完全處在破碎帶內(nèi),K16+460斷面只有一小部分處在破碎帶內(nèi)。
2) 水平收斂測(cè)線(xiàn)a比測(cè)線(xiàn)b收斂值大,模擬值約介于兩者之間。
這是因?yàn)闇y(cè)線(xiàn)a在上臺(tái)階開(kāi)挖后布設(shè),測(cè)線(xiàn)b經(jīng)下臺(tái)階開(kāi)挖后布設(shè),測(cè)線(xiàn)b經(jīng)上臺(tái)階開(kāi)挖,位移已得到一定程度的釋放,在數(shù)值模擬中,水平測(cè)線(xiàn)取的是上下臺(tái)階交界線(xiàn),因此收斂值介于測(cè)線(xiàn)a、b之間。
3) 實(shí)測(cè)圍巖變形收斂距離大于數(shù)值模擬變形收斂距離,且此后變形速度維持在一定水平,表現(xiàn)出圍巖流變效應(yīng)。
以K16+445斷面為例,數(shù)值模擬拱頂下沉變形基本收斂(變形速度≤0.1 mm/d)時(shí),斷面距掌子面2.0R,且此后變形速度越來(lái)越小趨近于零,但是實(shí)測(cè)拱頂下沉變形基本收斂時(shí),斷面距掌子面約為2.8R,且此后的3.2R范圍內(nèi)監(jiān)測(cè)變形速度一直維持在0.05~0.08 mm/d之間。
隧道開(kāi)挖過(guò)程中,洞身應(yīng)力得到釋放,圍巖應(yīng)力產(chǎn)生重分布。理論上,如果在這種重分布完成后其圍巖應(yīng)力不超過(guò)圍巖的支護(hù)條件和強(qiáng)度,則隧洞的開(kāi)挖將會(huì)保持穩(wěn)定和安全,但實(shí)際情況并非如此。在實(shí)際的開(kāi)挖過(guò)程中,不論掌子面的開(kāi)挖距離如何,甚至掌子面的空間效應(yīng)作用較大時(shí),隨著時(shí)間的推移,洞室仍有失穩(wěn)甚至坍塌的可能出現(xiàn)。這表明,隧道圍巖的變形,是隨著時(shí)間的變形而不斷繼續(xù)演變和發(fā)展的,而更多的工程實(shí)際表明,流變是普遍且長(zhǎng)期存在的。
1) 當(dāng)圍巖中存在破碎帶,尤其是斜向破碎帶時(shí),隧道開(kāi)挖不僅影響破碎帶區(qū)的圍巖變形,對(duì)破碎帶前后的圍巖變形也有較大影響,在采取注漿加固時(shí)應(yīng)適當(dāng)擴(kuò)大注漿范圍。在滲流影響下,高水位下破碎帶處圍巖開(kāi)挖變形急劇增大,對(duì)比注漿前,破碎帶注漿后圍巖變形得到明顯改善。
2) 破碎帶注漿加固后,實(shí)際監(jiān)測(cè)變形值比數(shù)值模擬變形值大很多,要區(qū)分破碎帶與其兩側(cè)圍巖的加固效果,在數(shù)值模擬時(shí)要弱化破碎帶加固后的圍巖參數(shù);同時(shí)要注重圍巖的流變效應(yīng),不能忽視距掌子面較遠(yuǎn)斷面的監(jiān)控工作,要密切關(guān)注圍巖的動(dòng)態(tài)變形,充分發(fā)揮施工監(jiān)控的作用。
[1] 張有天.巖體水力學(xué)與工程[M].北京:中國(guó)水利水電出版社,2005.
[2] 仵彥卿,張倬元.巖體水力學(xué)導(dǎo)論[M].成都:西南交通大學(xué)出版,1995.
[3] 楊天鴻,李連崇,唐春安.巖石破壞過(guò)程中滲流-損傷關(guān)系的認(rèn)識(shí):兼對(duì)唐紅俠博士提問(wèn)的答復(fù)[J].巖石力學(xué)與工程學(xué)報(bào),2004,(24):4254-4257.
[4] 朱珍德,郭海慶.裂隙巖體水力學(xué)基礎(chǔ)[M].北京:科學(xué)出版社,2007.
[5] 紀(jì)佑軍,劉建軍,程林松.考慮流-固耦合的隧道開(kāi)挖數(shù)值模擬[J].巖土力學(xué),2011,32(4):1229-1233.
[6] 吉小明,王宇會(huì),陽(yáng)志元.隧道開(kāi)挖問(wèn)題中的流固耦合模型及數(shù)值模擬[J].巖土力學(xué),2007,28(增刊1):379-384.
[7] 李小青.隧道施工技術(shù)[M].北京:中國(guó)建筑工業(yè)出版社,2010.
[8] 劉英棨,張謝東,李佳莉,等.基于反分析法的淺埋隧道圍巖應(yīng)力釋放率研究[J].武漢理工大學(xué)學(xué)報(bào)(交通科學(xué)與工程版),2015,39(6):1259-1262,1268.
[9] 任文峰.高水壓隧道應(yīng)力場(chǎng)-位移場(chǎng)-滲流場(chǎng)耦合理論及注漿防水研究[D].長(zhǎng)沙:中南大學(xué),2013.