劉智華，吳昕怡

（武漢輕工大學土木工程與建筑學院，湖北 武漢 430023）

1 研究背景及方法

中國處于地下深埋的水電站洞室群常采用地下廠房型式，而地下洞室群埋深大、跨度大，穿越的地質條件復雜，常遭遇高應力引起圍巖流變等影響構筑物的長期安全穩(wěn)定性問題[1]。巖石的流變力學特性，是指在外部條件作用下，巖石應力和應變隨時間緩慢變化的過程與現(xiàn)象[2]。圍巖流變是影響地下廠房運營期圍巖長期安全穩(wěn)定性和支護結構安全性的主要因素，分析圍巖流變對支護結構的影響，并根據(jù)圍巖和襯砌支護結構的位移變形和應力變化特征[3]，可以預測地下深埋洞室的長期安全穩(wěn)定性。

目前地下工程圍巖流變作用下的長期穩(wěn)定性研究主要采用綜合判據(jù)的評判準則，利用數(shù)值仿真的手段評價地下工程圍巖的長期安全穩(wěn)定性[4]。本文主要的方法是通過研究巖體的蠕變特性參數(shù)，確定合適的蠕變本構模型[5]，利用ABAQUS有限元軟件計算圍巖沉降位移、位移速率、應力以及支護結構的受力狀況[6]等狀態(tài)信息進行綜合定性評價。

2 流變作用原理

研究巖石流變性質，對解決巖土工程的維護設計和長期穩(wěn)定性問題有十分重要的意義。其中，蠕變現(xiàn)象是巖土工程中顯現(xiàn)最明顯、對工程穩(wěn)定性影響最大的流變現(xiàn)象[7]，是巖石流變理論研究中的極為重要的內容。

在恒定載荷作用下，只要有充分長的時間，巖土體均能產生蠕變現(xiàn)象。但在不同的恒定載荷下，變形隨時間增長的蠕變曲線卻有差異。以應變ε為縱坐標，時間t為橫坐標，做應變與時間的關系曲線圖，如圖1所示，該曲線就是典型蠕變曲線。

圖1 巖石的典型蠕變曲線

當在巖石試件上施加一個恒定載荷σ時，巖石立即產生一瞬間彈性應變εoa（oa段）。這種變形可以近似認為在t=0完成，其應變?yōu)棣舘a=σ/E。

若載荷保持恒定且持續(xù)作用，應變則隨時間緩慢地增長，此時已進入到蠕變變形階段。一般可將蠕變變形分成3個階段：減速（過渡）蠕變階段、等速（穩(wěn)定）蠕變階段、加速蠕變階段。減速（過渡）蠕變階段是蠕變過程中的第一階段，即圖1中的ab段，此階段應變速率隨時間增加而減小。等速（穩(wěn)定）蠕變階段是蠕變過程中的第二階段，即圖1中的bc段，蠕變曲線近似一傾斜直線，即應變速率保持不變，由b點一直持續(xù)到c點。加速蠕變階段是蠕變過程中的第三階段，即圖1中的cd段，這一階段中，應變速率隨時間增加而增加，應變率由c點開始迅速增加，達到d點，巖石即發(fā)生破壞。

在實際地下工程的圍巖穩(wěn)定性計算中，通常只考慮前兩個階段的蠕變，而第三階段實際上試件已經達到了破壞狀態(tài)，其性質比較復雜，通常不考慮。

3 工程地質和相關計算參數(shù)確定

3.1 工程地質條件和工程概況

某水電站是修建于四川省雅礱江干流下游河段的控制性水庫梯級電站。地下廠房洞室群主要由廠房、母線洞、主變室、尾調室及尾水連接管等組成，母線洞布置于廠房和主變室之間。地下廠區(qū)圍巖類別以Ⅲ1類大理巖為主。地下洞室建成后，受自重應力和高地應力的影響，地下洞室圍巖會隨時間發(fā)生流變效應。

3.2 相關材料物理力學計算參數(shù)

根據(jù)大理巖力學試驗結果，模擬隧洞處于運營期的圍巖和襯砌的變形情況來分析地下洞室的穩(wěn)定性，在已知Ⅲ1級大理巖圍巖下采用最不利力學參數(shù)，相關材料物理力學參數(shù)如表1所示。某水電站大理巖流變模型參數(shù)如表2所示。

表1 材料物理力學參數(shù)

表2 水電站大理巖H-K流變模型參數(shù)

3.3 ABAQUS流變模型計算參數(shù)

本文在ABAQUS的黏彈性模型中Prony級數(shù)選取根據(jù)大理巖蠕變試驗的相關試驗數(shù)據(jù)，經過計算推導和轉換，Prony級數(shù)計算取值為=0.327；因為體積模量參數(shù)與時間無關，一欄保持空白；因要計算多場情況下圍巖-支護結構的長期穩(wěn)定性，由流變模型參數(shù)計算得=5年。在考慮土體自重應力和初始地應力場的情況下，在模型上表面施加11.34 MPa的均布荷載來模擬上部山體自重的作用；對模型施加σ11=35.7 MPa，σ22=20 MPa，σ33=10 MPa，模擬初始應力場的作用。

4 數(shù)值計算模型

4.1 模型概況

結合某水電站地下廠區(qū)洞室相關資料建立母線洞的隧道模型。為了模擬施工完成后隧道的圍巖變形和支護結構內力變化情況，采用ABAQUS建立三維深埋隧道模型和支護模型，分別如圖2、圖3所示。其中地下隧洞水平兩側最大距離為9.2 m，拱頂?shù)焦暗變蓚茸畲缶嚯x為7.89 m，模型上表面距離山體表面埋深為300 m，這里深埋隧洞模型的邊界尺寸取水平兩側長為50 m、上下兩側高為50 m、徑向洞身長度為30 m的長方體。模型左右邊界約束水平方向，前后面約束縱向位移，底邊約束水平和垂直方向。所建模型中隧洞土體移除后施加Φ28、長4.5 m和Φ32、長6.0 m兩種鋼筋型號的錨桿，所施加的錨桿尾部外漏0.5 m。兩種鋼筋型號的錨桿相間布置組成一排錨桿，每排錨桿間距1.5 m。在施加錨桿支護結構后，再施加厚度為80 cm混凝土襯砌支護結構，最終組成的整個錨桿襯砌支護結構，如圖3所示。整個錨桿和襯砌支護結構的一次性施加以模擬地下洞室建成后的支護條件。

圖3 錨桿和襯砌支護結構

4.2 計算結果分析

在本隧洞計算模型中，主要考慮隧洞建成后圍巖的應力變化和位移變形情況以及觀察襯砌支護結構的應力變化情況[8-9]，來確定以及預測隧洞的長期安全性和穩(wěn)定性。

4.2.1 圍巖應力和位移分析

隧洞建成后，由于圍巖自重應力和高地應力的作用，且圍巖土體是黏彈性材料，巖土體會發(fā)生流變效應。如圖4所示，流變40年后，巖土體達到穩(wěn)定狀態(tài)時，圍巖拱頂拱底受到的較大拉應力作用，其最大拉應力達到53.68 MPa，拱腰兩側會受到最大壓應力為3.83 MPa的作用。

圖4 流變作用第40年圍巖應力云圖

如圖5所示，流變作用下，圍巖拱頂和拱底的拉應力逐漸減小，達到穩(wěn)定時，圍巖拱頂處最大拉應力由69.16 MPa減小到53.68 MPa，圍巖拱底處最大拉應力由48.13 MPa減小到32.30 MPa，圍巖拱頂處的應力變化幅度最大，因此在研究地下深埋洞室圍巖的長期穩(wěn)定性時，需重點研究圍巖拱頂處的應力變化情況。

圖5 流變作用40年圍巖應力曲線圖

由圍巖應力變化規(guī)律可知，在研究圍巖位移變化規(guī)律時，須重點研究隧洞拱頂處圍巖的位移。由圖6所示的流變作用40年圍巖拱頂拱底位移曲線圖可知，在流變作用下，圍巖拱頂位移沉降最大，最大值達到11.75 mm。圍巖初期變形速率較大，前8年平均變形速率為1.04 mm/年，第9—20年的平均變形速率為0.24 mm/年，第21—40年的平均變形速率為0.028 mm/年，第21—40年總的合位移值為0.56 mm。圍巖拱底處的沉降位移值最大為7.93 mm。

圖6 流變作用40年圍巖拱頂拱底位移曲線圖

4.2.2 支護結構分析

在ABAQUS中，壓應力表現(xiàn)為正值，拉應力為正值。在進行地下深埋洞室的長期安全穩(wěn)定性研究時，襯砌結構應力分布特征和位移變化規(guī)律是進行流變分析的重要評價指標。襯砌結構在流變作用第40年襯砌最小主應力云圖如圖7所示，作用在襯砌上的應力表現(xiàn)為壓應力，最大值為11.09 MPa，由于襯砌是混凝土結構，作用在襯砌上的壓應力遠小于混凝土的抗壓強度值。由圖8所示的流變作用第40年襯砌沉降位移云圖可知，襯砌結構的位移變形主要位于拱頂處，當達到穩(wěn)定狀態(tài)時，襯砌拱頂?shù)某两盗孔畲螅两底畲笾禐?1.99 mm，其變形量占整個襯砌結構的1.49%，對整個襯砌結構影響很小。由流變作用下第40年襯砌所受壓應力及沉降位移值可知，襯砌結構滿足流變作用下的長期安全穩(wěn)定性要求。

圖7 流變作用第40年襯砌最小主應力云圖

圖8 流變作用第40年襯砌沉降位移云圖

由圖9所示的流變作用40年錨桿應力分布云圖可知，流變作用達到穩(wěn)定狀態(tài)時，錨桿結構受到最大拉應力主要位于拱腰兩側位置，最大拉應力達到119.8 MPa。由于圍巖拱頂區(qū)域為拉應力區(qū)，也是變形最大的區(qū)域，因此對這部分區(qū)域的錨桿進行應力監(jiān)測和研究可知錨桿所受拉應力較小，有的錨桿沒有拉應力作用。錨桿所受最大拉應力始終遠小于錨桿的抗拉設計值，錨桿支護結構滿足長期穩(wěn)定性要求。

圖9 流變作用第40年錨桿應力分布云圖

綜合襯砌和錨桿的應力應變規(guī)律研究分析可知襯砌錨桿支護結構滿足流變作用下的長期穩(wěn)定性要求。

5 結束語

本文以某水電站母線洞為例，采用三參量廣義Kelvin黏彈性本構模型，推導出相關力學參數(shù)，在ABAQUS中建立黏彈性流變模型來模擬隧洞建成后的長期變形以確定地下深埋隧洞的安全穩(wěn)定性。主要結論如下。

采用黏彈性流變模型來表征隧洞建成后圍巖變形效應特征是合理的。土體的自重應力使得水電站地下深埋洞室圍巖在相當長的一段時間內還會產生緩慢變形，這是影響建成后的隧洞的安全穩(wěn)定性的主要原因。

本文是在隧洞建成后的基礎上對其進行安全穩(wěn)定性分析，所計算出的結果體現(xiàn)了隧洞在流變時間內會產生較大的位移變形，同時也體現(xiàn)了隧洞的支護結構能起到非常好的支護作用。

計算結果表明地下深埋隧洞的圍巖變形和支護結構達到相對穩(wěn)定的狀態(tài)需要相當長的時間，期間襯砌結構應力變化情況表明所用支護結構方式滿足實際工程安全需求。所采用的數(shù)值模擬方法可以運用到相類似的工程應用當中，以分析和評價工程構筑物的建設完成后的長期安全穩(wěn)定性能。