考慮水盆效應隧道道床上浮力學特征及計算分析

摘 要：圍繞富水隧道道床上浮底鼓工程問題，本文探索將水盆效應引入隧道無砟道床上浮底鼓病害的研究中。首先，基于應變能和能量理論，結合結構層間接觸力學條件，開展隧道道床應力和變形量推導，以驗證水盆效應與道床上浮關聯的合理性。其次，構建道床變形撓曲線方程和應力和最大變形量公式，得到基底0.5 m水頭高度可以引起粘結失效道床上浮底鼓，并通過設定上浮臨界水頭值劃分上浮三階段：1） 初期入滲匯水階段；2） 局部積水階段；3） 水盆上浮階段。最后，建立不同排水狀態下道床橫向及縱向上浮量計算公式，形成基底水流量差、壓力速率分別與上浮速率及上浮量的關系曲線，形成隧道道床上浮的最小水頭高度的計算方法，以助力道床上浮控制的研究和治理措施的改進。

關鍵詞：富水隧道；道床上?。凰栊涣W特征

中圖分類號：TU43

文獻標志碼：A

我國各地的許多隧道均存在滲漏水情況。隧道滲漏水會導致道床混凝土劣化［1-2］、道床上?。?］，道床上浮是隧道工程中較為常見的一種病害［4-5］，道床上浮會嚴重危害結構穩定和車輛安全及平順性。國外已經有許多因膨脹圍巖引起的仰拱底鼓病害［6］，通過理論公式推導可以揭示道床脫空上浮的本質原理［7-8］，然而目前研究主要集中在道仰拱底鼓引起的道床上浮及隧道仰拱底鼓量的預測。仰拱底鼓的原因有許多種，ANAGNOSTOU［9］采用彈塑性理論模擬膨脹圍巖，分析了膨脹圍巖區隧道仰拱變形及其底鼓機制。BUTSCHER等［10］分別從微觀及宏觀角度研究了隧道開挖導致周邊水文條件變化，地下水流入隧道周邊膨脹性圍巖導致仰拱隆起。鐘祖良、劉新榮等［11］采用連續介質力學原理，推導了軟巖隧道仰拱底鼓量的理論公式，可為類似隧道底鼓量的預測提供依據。師亞龍、陳禮偉等［12］運用彈塑性力學原理分別推導了各種因素引起的底鼓量，并得出了底鼓量預測的組合公式。WITTKE GATTERMANN［13］采用彈塑性力學原理，考慮材料各向異性，并引入時間變量，建立了膨脹性圍巖本構模型。秦巴列維奇［14］基于松散土體壓模試驗分析，采用土體極限平衡理論，開展了底鼓巷道土壓力及其機理的計算分析。

國內外的研究主要集中在仰拱底鼓引起道床上浮，尚未見有水盆效應引起道床上浮相關的研究。上述理論研究均針對的是隧道底鼓現象，且均為外部因素如圍巖壓力、圍巖膨脹性等方面引起隧道底鼓。然而，結合隧道工程上浮底鼓問題研究，道床板粘接強度完整，需要較高臨界水頭引起道床與仰拱填充之間的上浮底鼓。但實際工程上浮水頭高度較小，有必要開展隧道道床上浮力學問題的計算分析的研究。

1 水盆效應及計算基本假定

建筑工程中首先提出的一個“水盆效應”名詞，根據李衛平［15］的研究，水盆效應是指基坑回填土的分層夯實不達標，使得地表水、地下水以及施工用水等四周水匯入坑內和地下室底板下面。在基坑四周為不透水土層和坑內水滯留情況下，當底板滯留水的上浮浮力大于地下室抗浮力時，則地下室發生上浮現象。水盆效應示意圖見圖1（a）。

在水盆效應及二襯背后水壓力共同作用下，無砟道床可以在較小的水頭高度下發生上浮，本文創新將 “水盆效應”引入隧道道床上浮分析中，如圖2所示。

本文依托貴陽市軌道交通二號線無砟道床發生上浮的案例，通過查閱相關資料，從理論分析及數值模擬分析等方面，深入研究水壓力引起的無砟道床上浮現象（如圖1（b）所示）。因此，將水盆效應引入隧道無砟道床上浮病害的研究中，運用能量理論對無砟道床的應力和變形量理論解進行推導，對水盆效應如何影響道床上浮進行探究，助力道床上浮機理的研究以及相應治理措施的改進。本文采用能量法計算道床上浮量。能量法是求解任何彈性結構的基本方法。在外力作用下，利用功能原理求解結構指定點的位移稱為能量法，它也是有限單元法的基礎。能量法有如下幾個特點：① 求解簡便、迅速；② 不受材料和形狀的限制；③ 適用于各種彈性結構、非彈性結構、塑性結構等；④ 可用來求解各種超靜定問題?；谏鲜隼碚摲椒?，道床上浮求解的基本假設如下：

1）連續性假設

隧道道床混凝土是連續的，不考慮內部孔隙影響，可以用一個連續函數表示道床混凝土的應力、應變、位移等相關力學參數。

2）均勻性假設

道床混凝土是均勻的，即混凝土內部任取一微小單元，其力學性能均與道床整體的力學性能相同。

3）各向同性假設

道床混凝土是各向同性的，即任何方向上具有相同的力學性能。

4）小變形假設

道床混凝土的變形量遠小于道床的幾何尺寸。

2 考慮水盆效應作用下道床應力和變形量推導

2.1 隧道道床上浮三階段劃分

1）應變能和功能原理

應變能：由于外部荷載引起結構變形而內部儲存的能量稱為應變能Vε。

功能原理：在彈性體中，外力所做的功全部轉化為彈性體內的應變能，可以認為在彈性體變形過程中，彈性體內積蓄的應變能等于外力所做的功，即

Vε=W（1）

根據功能原理可以計算出彈性體發生基本變形的應變能?；诏B加原理，結構組合變形（軸向拉壓+圓軸扭轉+平面彎曲）應變能表示為

卡氏定理的主要內容是：線彈性范圍內，結構應變能對某個力的偏導數，等于結構沿此力方向上的位移（轉角），即

因此，對于3種基本變形的位移，有

軸向拉壓：

圓軸扭轉：

平面彎曲：

組合變形：

2）臨界水頭公式推導

沿縱向取1 m道床板進行受力分析，此時道床相當于兩端固定的梁，設道床寬度為l，道床厚度為d，如圖3所示。

其受力簡圖如圖4所示。

道床板受到自重G、道床與仰拱之間粘接力fz以及浮力f浮影響，當f浮大于道床板自重G與粘接力fz之和時，道床板即會發生上?。ㄉ蠐希Ｉ细l件為

f浮≥G+fz（8）

γh≥ρgd+σ（9）

式中：道床自重G=ρgV=ρgdl，粘接力fz=σs=σl，浮力f浮=γhs=γhl。

取道床板厚d=200 mm；混凝土密度ρ=2 500 kg/m3；重力加速度g=9.8 m/s2，水的重度γw=9.8 kN/m3。粘接強度σ=0，則

如果道床與仰拱填充間粘接強度完全失效，則道床僅在0.5 m水頭高度下即可發生上浮。

將20 cm板厚臨界水頭高度h和粘接強度σ的關系繪制成函數圖像如圖5。

從圖6中可以看出道床板厚對臨界水頭的影響遠遠小于道床與仰拱填充之間的粘接強度。若粘接強度失效，則在很小的水頭高度對應的水壓力作用下，道床即可發生上浮。

根據上述隧道底部結構上浮與層間粘接強度的發生發展規律分析，隧道底部結構上浮過程可以分為三個階段（見圖7）： 1） 初次上浮階段，即入滲匯水；2） 第二次上浮階段，即基底積水階段；3） 第三次上浮階段，即脫離上浮階段。由隧道基底水貫通入仰拱內，填充到結構層底部引起的結構上浮的階段稱為隧道初次上?。挥捎诨姿M入仰拱內后形成滯留水引起的結構局部變形稱為第二次上浮底鼓階段；隧道第三次上浮是由于結構層間粘結失效，形成具有水盆效應的滯留水層。

2.2 道床變形撓曲線方程推導

由圖4的受力分析圖可知：將道床簡化為梁單元，得到道床發生變形后的軸線，稱為道床撓曲線。

撓曲線方程為

ω=f（x）（10）

在小變形情況下的道床轉角和撓度有如下關系：

梁在只受到彎矩作用下的曲率與彎矩關系為

由于剪力對道床變形影響很小，可以忽略不計，而曲率和彎矩均為x的函數，即

在數學上曲率與方程有如下關系：

在小變形情況下ω=l；1+ω′2≈1，此時

因在土木工程領域一般以撓度向下為正，所以有

EIω″≈－M（x）（16）

將上式積一次分可得轉角方程

再次積分可得撓曲線方程

結合兩端固定梁邊界條件：x=0和x=l時θ=0；ω=0，荷載為均布荷載q，帶入式（18）可得撓曲線方程

2.3 應力和最大變形量公式推導

采用能量法計算道床上浮量。由于道床為對稱結構受到對稱荷載的作用，可取一半結構進行分析，右端簡化為定向支撐，解除右側約束，代之以支座反力，并將浮力、粘結力、自重轉化為一個線荷載q，如圖8所示。

設右端為原點，則彎矩函數為

M（x）對MB求偏導，得

根據卡氏定理，結合邊界條件，θB=0，可得

為求得最大位移，在右端虛加集中力Ff，如圖9所示。

則彎矩函數變為

M（x）對Ff求偏導，得

根據卡式定理，得

3 考慮不同排水階段下道床上浮量的推導

根據高新強［16］推導的深埋隧道限排情況下襯砌背后水壓力公式，若隧道為淺埋隧道且為全包型防水，襯砌完全不透水，襯砌背后不設排水系統，則襯砌背后的水壓力即為靜水壓力，不進行折減。襯砌背后水壓力為

3.1 不排水情況討論

假設滲入到道床與仰拱填充面間的水壓力與襯砌外水壓力相同，若隧道為全封堵情況或排水管完全堵塞，則相當于滲透系數K1→0，根據式（28）可得水壓力為

即水壓力與初始水頭對應的水壓力相等。

1）橫向脫空情況

2）縱向脫空情況

3.2 堵水限排情況討論

假設滲入到道床與仰拱填充面間的水壓力與襯砌外水壓力相同，若隧道為堵水限排情況，排水管部分堵塞，則相當于滲透系數K1減小，根據式（28）可知水壓力會適當增大。設增大后的水壓力為ξ1γwh。

1）橫向脫空情況

2）縱向脫空情況

3.3 以排為主情況討論

假設滲入到道床與仰拱填充面間的水壓力與襯砌外水壓力相同，若隧道為以排為主的情況，排水管效率增大，則相當于滲透系數K1增大。

根據式（28）可知水壓力會大大減小。設減小后的水壓力為ξ2γwh。

1）橫向脫空情況

2）縱向脫空情況

當ξ2→0時，ξ2γwh→0，由此可知在道床上增設排水孔可大大減小道床與仰拱填充之間的水壓力。

3.4 不連通情況討論

假設地下水進入道床與仰拱填充之間后無法排出，由于水的不可壓縮性，地下水會將道床頂起。顯然道床的上浮量與滲入道床與仰拱填充之間的水量大小有一定的關系，下面將以橫向脫空為例推導滲水量與道床上浮量之間的關系（見圖10）。

1）橫向脫空情況

已知道床變形后的撓曲線方程為式（19），在數學上，積分表示函數與x軸圍成的面積，如圖11所示，將式（19）積分1次可得道床上浮后的橫截面積增加量A：

因橫向脫空情況下l=B，經計算可得

此時可得道床上浮后單位長度脫空體積增加量V：

V=A×1=A（38）

因為道床脫空體積完全被水充滿，所以可得單位長度滲水量Vw：

又因為道床上浮量最大值表達式為

因此，可以得出單位長度滲水量與橫向脫空上浮量最大值之間的關系：

若流入與流出仰拱填充面的水量不相等，有以下幾種情況：

對式（41）進行t（時間）求導，可得

假設過水斷面面積的流入與流出相等，

將式（41）帶入式（44），得

根據公式可得到如下曲線，如圖12、13所示。上浮速率與道床與仰拱填充面水壓力變化速率成正比，上浮速率會隨著流量差的增大而增大，并且在流量差越大的情況下，上浮量的變化也越大，這種變化一開始是緩慢的，隨著流量差的變化而逐漸變快。

2）縱向脫空情況

對于縱向脫空情況，設總滲水量為V′w ：

式中：l為縱向脫空長度。

因此可以得出總滲水量與縱向脫空上浮量最大值之間的關系：

若流入與流出仰拱填充面的水量不相等，則有

根據公式可得到如下曲線，如圖14、15所示。上浮速率與道床與仰拱填充面水壓力變化速率成正比，上浮速率會隨著流量差的增大而增大，并且在流量差越大的情況下，上浮量的變化也越大，這種變化一開始是緩慢的，隨著流量差的變化而逐漸變快。變化傾向與橫向脫空情況一致，但是因為縱向脫空長度的關系，在數值上發生了明顯的變化。

4 結論

從隧道道床上浮的典型背景工程問題出發，將隧底基巖等效為均勻滲流體系，引入應變能和功能原理，結合結構層面接觸力學條件，開展考慮水盆效應隧道道床上浮力學特征及計算分析研究。

1）建立板厚臨界水頭高度和粘接強度的關系函數，形成隧道道床上浮臨界水頭計算公式，得到在0.5 m水頭高度可以引起粘結失效道床的上浮，并建立了道床變形撓曲線方程和應力和最大變形量公式。

2）根據隧道基底不同的排水狀態，將承受對稱荷載的道床進行力學簡化分析可知，在粘接強度失效條件下，道床上浮僅需要較小水頭高度。并建立了流量差與上浮速率、上浮量關系，和壓力速率與上浮速率、上浮量關系曲線。隨著道床與仰拱填充面水壓力的增高，上浮量逐漸增大。

3）形成了道床上浮三階段的底鼓演化過程?？紤]粘結強度，通過理論計算引起道床上浮的最小水頭高度，并結合相關力學知識對水盆效應引起道床橫向及縱向上浮量進行估算，驗證了水盆效應引入道床上浮研究中的合理性。考慮水盆效應推導隧道道床上浮底鼓公式，可以較好地表征地下工程上浮問題。

參考文獻：

[1]許敏娟， 陳瑩瑩， 劉浩. 基于全景圖像的隧道滲漏水病害識別方法[J]. 鐵道學報， 2023， 45（8）： 184-192.

[2] NIEDOSTATKIEWICZ M. Degradation of the concrete railway track bed located in the vicinity of the loading wharf［J］. Polish Maritime Research， 2018， 25（4）： 149-153.

［3］ 郭志州， 欒忠杰. 香爐山隧道漏水原因分析及改造設計［J］. 北方交通， 2015（7）： 97-99.

［4］ 華科. 成貴鐵路隧道滲漏水缺陷整治技術［J］. 高速鐵路技術， 2021， 12（4）： 55-59.

［5］ 何昌國， 黃棋， 趙萬強. 貴廣鐵路巖溶隧道水害原因分析及整治措施［J］. 高速鐵路技術， 2018（增刊1）： 67-73.

［6］ CHARPENTIE D， MOSSER-RUCK R， CATHEL M. Oxidation of mudstone in a tunnel （Tournemire， France）： consequences for the mineralogy and crystal chemistry of clay minerals［J］. Clay Minerals， 2004， 39（2）： 135-149.

［7］ SHARIFZADEH M， TARIFARD A， MORIDI M A. Time-dependent behavior of tunnel lining in weak rock mass based on displacement back analysis method［J］. Tunnelling and Underground Space Technology Incorporating Trenchless Technology Research， 2013， 38： 348-356.

［8］ MOYA I， AGREDA E， MORALES E， et al. A review of expansive phenomena in Wagenburg North Tunnel［J］. Revista D La Academia Colombiana D Ciencias Exactas Fisicas Y Naturales， 2010， 33： 455-468.

［9］ ANAGNOSTOU G . A model for swelling rock in tunnelling［J］. Rock Mechanics amp; Rock Engineering， 1993， 26（4）： 307-331.

［10］BUTSCHER C ， MUTSCHLER T， BLUM P. Swelling of clay-sulfate rocks： a review of processes and consols ［J］. Rock Mechanics and Rock Engineering， 2016， 49（4）： 1533-1549.

［11］鐘祖良， 劉新榮， 王道良， 等. 桃樹埡隧道底鼓發生機理與防治技術研究［J］. 巖土工程學報， 2012， 34（3）： 471-476.

［12］師亞龍， 陳禮偉， 裴濤濤， 等. 軟巖隧道底鼓機理及底鼓量計算方法［J］. 鐵道建筑， 2014（5）： 54-56.

［13］WITTKE-GATTERMANN P. Computation of strains and pressures for tunnels in swelling rocks［J］. Tunneling amp; Underground Space Technology， 2004， 19（4）： 422-423.

[14]秦巴列維奇. 礦井支護[M]. 北京： 煤炭工業出版社， 1953： 98-696.

［15］李衛平. 地下水浮力計算和抗浮基礎設計［J］. 江蘇建筑， 2013（4）： 78-83.

［16］高新強. 高水壓山嶺隧道襯砌水壓力分布規律研究［D］. 成都：西南交通大學， 2005.

Calculation and Analysis of Track Bed Floating Characteristics

of Karst Tunnel Considering Water Basin Effect

Abstract：

Focusing on the problem of floating floor heave on the track bed of water-rich tunnel， this paper explores the introduction of water basin effect into the study of floating floor heave disease on the ballastless track bed of tunnels. Firstly， based on the strain energy and energy theory， combined with the contact mechanical conditions between structural layers， the stress and deformation of tunnel track bed are deduced to verify the rationality of the correlation between water basin effect and track bed floating. Secondly， the deformation deflection curve equation of track bed and the formula of stress and maximum deformation are constructed. It is obtained that the floating floor heave of track bed can cause bonding failure at the height of 0.5m water head， and the floating stage is divided into three stages by setting the critical water head value of floating： 1） initial infiltration catchment stage; 2） local ponding stage; 3） basin floating stage. Finally， the calculation formulas of horizontal and vertical floating amount of track bed under different drainage conditions are established， the curves of relationship between base water flow difference， pressure rate， floating rate and floating amount are formed， and the calculation method of minimum head height of tunnel track bed floating is constructed， so as to help the research of track bed floating control and the improvement of treatment measures.

Key words：

water-rich tunnel; ballast bed floating; basin effect; mechanical characteristics