深埋圓形隧洞非飽和土-襯砌結構系統的動力特性

劉洪波， 劉　潔 ， 聞敏杰

(1. 鄭州航空工業管理學院，鄭州　450015，2. 嘉興職業技術學院，浙江　嘉興　314036)

深埋圓形隧洞非飽和土-襯砌結構系統的動力特性

劉洪波1， 劉潔1， 聞敏杰2

(1. 鄭州航空工業管理學院，鄭州450015，2. 嘉興職業技術學院，浙江嘉興314036)

采用解析方法在頻率域內研究了深埋圓形隧洞非飽和土-襯砌結構系統的動力特性。基于Bishop非飽和土有效應力原理和單相流固結理論，建立了單相流條件下非飽和土的運動方程。通過位移勢函數解耦，得到了非飽和土在簡諧荷載作用下的動力響應解答。將襯砌結構視為均勻彈性介質，基于彈性理論，得到了簡諧荷載作用下圓柱形襯砌結構的穩態形式解答。利用土體與襯砌界面的連續性條件和襯砌結構內邊界上的邊界條件，確定了表達式的待定系數。在此基礎上，考察了飽和度、吸力折減系數、襯砌厚度等參數對非飽和土響應幅值的影響。結果表明：隨著襯砌厚度的增加，系統響應幅值逐漸減小；飽和度和吸力折減系數對孔隙水壓力幅值的影響較大，而對徑向位移和環向應力幅值的影響較小。

非飽和土；單相流；圓柱形襯砌結構；動力特性；解析解

土體具有豐富且較為復雜的力學特性，尤其是地下結構受到爆炸、沖擊等內源動力荷載作用時常導致土體發生變形，結構和土體之間時有發生擾動，由此，土體的力學性能對地鐵、電力管線、污水管道和海底隧道等地下工程的安全運營具有直接的影響。近年來，國內外學者在研究具有圓形隧道周圍土體的動力響應方面取得了大量成果，可用于地下隧道工程抗爆設計、穩定性分析等。這些研究或是把周圍土視為單相介質[1-5]，或是把周圍土視為流固兩相介質[6-10]，分析土體的動力特性。

然而，我國隧道分布密集區都以軟土的形式存在，以往的單相彈性地基模型和流固耦合兩相介質模型與實際工程中常見的軟土有很大差異。這些軟土常為流、固、氣三相介質，將土體視為單相或兩相介質，不能精確地反映土體的真實特性[11]。為此，徐長節等[12-14]假定土體中氣相僅以分散氣泡形式存在于水中，將水氣混合物視為均勻流體，采用Biot兩相介質理論分別研究了具有高飽和度的準飽和彈性土和黏彈性土動力特性。但是，上述研究僅針對飽和度在95%以上的土體。為了更好地描述具有圓柱形襯砌結構非飽和土的動力行為，本文基于Bishop非飽和土有效應力原理和單相流固結理論，采用解析方法在頻率域內研究了深埋圓形隧洞非飽和土-襯砌結構系統的動力特性，對比分析了無襯砌和彈性襯砌、邊界透水與不透水、非飽和土與飽和土等情形下的動力響應差異，并考察了物性和幾何參數對非飽和土響應幅值的影響。

1　非飽和土體的控制方程及求解過程

如圖1所示，非飽和土體中有一內外半徑分別為R1和R2，厚度為d=R2-R1，圓心到中曲面的距離為a的圓柱形隧洞，其內邊界承受圓頻率為ω的徑向均布簡諧激振力q0eiωt(i2=-1)。假設非飽和土-襯砌系統為小變形，且襯砌和土體完全接觸，即滿足位移和應力連續。不計體力時，該問題可視為軸對稱平面應變問題。

圖1　非飽和土中圓形襯砌隧洞Fig.1 Circular lined tunnel in unsaturated soil

根據單相固結理論[15]，非飽和土體滿足單相流條件，即飽和度適中，忽略排氣過程，各處氣體壓力均與大氣壓力相等，僅考慮水力梯度和吸力梯度下孔隙水的流動。因此，土體的應力關系用Bishop提出的非飽和土有效應力原理表示為：

(1)

式中：σSE,σ*,σ分別表示土體的有效應力、凈應力和總應力；Sa和χ分別為吸力和吸力折減系數；pa為土體的孔隙氣壓力，當氣體壓力始終為大氣壓時，孔隙氣壓pa可取為0；{δ}={1,1,1,0,0,0}T

于是，在極坐標下土體的總應力-應變關系可表示為：

(2)

土體的運動平衡方程為：

(3)

孔隙流體運動方程為：

(4)

式中：k為液體的滲透系數，而Sr為飽和度。

非飽和土中孔隙水連續方程為：

(5)

流相相對于固相的滲流速度表示為：

(6)

根據沈珠江等[15]研究成果，可知：

(7)

結合式(5)～(7)，可得土體的孔隙水壓力滿足如下本構關系：

(8)

將式(2)代入式(3)和(8)，式(4)代入式(8)，可得：

(9)

式中：λc=λ+χαM,M=Hw/n,α=Sr(1+n),b=ρwg/k,m=ρw/nSr

對于系統作穩態振動，令p=GPeiωt,ur=R2Uηeiωt,wr=R2Wηeiωt，其中，P,Uη,Wη分別為p,ur,wr的無量綱量，并引入無量綱量和常數：

(10)

同時，利用文獻[14]的解法和無窮遠處位移為零的邊界條件，并略去eiωt項，土體的徑向位移和流體相對于土骨架的徑向位移表達式分別為

(11)

利用本構關系式(8)，可得孔隙水壓力的表達通式為

(12)

根據關系式(2)，可得土體的徑向和環向總應力分別為

(13)

(14)

m1=

(15a)

(15b)

(16)

2　非飽和土體的控制方程及求解過程

將襯砌視為均勻彈性介質，根據彈性理論，動力荷載作用下以位移表示的運動方程為[6-7]：

(17)

式(17)易解得：

(18)

根據應力-位移的本構關系，可得徑向和環向應力分別為：

(19)

(20)

滿足如下邊界條件：

(21)

② 在襯砌與土體界面處，有

(22)

③ 若外邊界不透水

wr=0r=R2

(23)

④ 若外邊界自由透水

p=0，r=R2

(24)

將式(11)～(14)和式(18)～(20)代入式(21)～(24)，易求得待定系數Ai,Di(i=1,2)的具體表達式，即可得到深埋圓形隧洞非飽和土-襯砌系統穩態振動的解析解。

3　圖形分析及討論

據文獻[13-14]取值如下：

η=1，γ=d/R2=0.005，n=0.4，

(25)

3.1對比分析

(a) 徑向位移隨頻率變化(a) 環向應力隨頻率變化(a) 孔隙水壓力隨頻率變化(a) 孔隙水壓力隨頻率變化

(b) 徑向位移隨半徑變化圖2 有無襯砌時徑向位移幅值的比較Fig.2Comparisonoftheamplitudeoftheradialdisplacementwithorwithoutlining(b) 環向應力隨半徑變化圖3 有無襯砌時環向應力幅值的比較Fig.3Comparisonoftheamplitudeofthehoopstresswithorwithoutlining(b) 孔隙水壓力隨半徑變化圖4 有無襯砌時孔隙水壓力幅值的比較Fig.4Comparisonoftheamplitudeoftheporewaterpressurewithorwithoutlining(b) 孔隙水壓力隨半徑變化圖5 透水情形下孔隙水壓力幅值Fig.5Theamplitudeoftheporewaterpressureforpermeablecondition

3.2飽和度的影響

3.3吸力折減系數的影響

3.4襯砌厚度的影響

(a) 徑向位移隨頻率變化(a) 環向應力隨頻率變化(a) 孔隙水壓力隨頻率變化(a) 徑向位移隨頻率變化

(b) 徑向位移隨半徑變化圖6 飽和度對徑向位移幅值的影響Fig.6Influenceofsaturationontheamplitudeofradialdisplacement(b) 環向應力隨半徑變化圖7 飽和度對環向應力幅值的影響Fig.7Influenceofsaturationontheamplitudeofthehoopstress(b) 孔隙水壓力隨半徑變化圖8 飽和度對孔隙水壓力幅值的影響Fig.8Influenceofsaturationontheamplitudeoftheporewaterpressure(b) 徑向位移隨半徑變化圖9 透水情形下飽和度對徑向位移幅值的影響Fig.9Influenceofsaturationontheamplitudeofradialdisplacementforpermeablecondition

(a) 徑向位移隨頻率變化(a) 環向應力隨頻率變化(a) 孔隙水壓力隨頻率變化(a) 徑向位移隨頻率變化

(b) 徑向位移隨半徑變化圖10 吸力折減系數對徑向位移幅值的影響Fig.10Influenceofreductioncoefficientofsuctionontheamplitudeoftheradialdisplacement(b) 環向應力隨半徑變化圖11 吸力折減系數對環向應力幅值的影響Fig.11Influenceofreductioncoefficientofsuctionontheamplitudeofthehoopstress(b) 孔隙水壓力隨半徑變化圖12 吸力折減系數對孔隙水壓力幅值的影響Fig.12Influenceofreductioncoefficientofsuctionontheamplitudeoftheporewaterpress(b) 徑向位移隨半徑變化圖13 襯砌厚度對徑向位移幅值的影響Fig.13Influenceofliningthicknessontheamplitudeoftheradialdisplacement

(a) 環向應力隨頻率變化(b) 環向應力隨半徑變化(a) 孔隙水壓力隨頻率變化(b) 孔隙水壓力隨半徑變化圖14 襯砌厚度對環向應力幅值的影響Fig.14Influenceofliningthicknessontheamplitudeofthehoopstress圖15 襯砌厚度對孔隙水壓力幅值的影響Fig.15Influenceofliningthicknessontheamplitudeoftheporewaterpressure

4　結　論

基于單相流固結理論和Bishop非飽和土有效應力原理，采用解析方法在頻率域內研究了深埋圓形隧洞非飽和土-襯砌結構系統的動力特性，并對比了無襯砌和彈性襯砌兩種條件下非飽和土的動力特性差異，得到以下結論：

(1) 有彈性襯砌時非飽和土的徑向位移幅值明顯小于無彈性襯砌的，且隨著襯砌厚度γ的增加，系統響應幅值逐漸減小。適當增加襯砌厚度γ可避免系統的共振。

[1] Forrestal M L, Sagartz M J. Radiated pressure in an acoustic medium produced by pulsed cylindrical and spherical shells[J]. AMSE Journal of Applied Mechanics, 1971,38：1057-1060.

[2] Zakout U, Akkas N. Transient response of a cylindrical cavity with or without a boned shell in infinite elatic medium[J]. International Journal of Engineering Science, 1997, 35(12/13): 1203-1220.

[3] Seyyed M, Hasheminejad Mojtaba K. Effect of imperfect bonding on axisymmetric elastodynamic resonpose of a lined circular tunnel in poroelastic soil due to a moving ring load[J]. International Journal of Solids and Structures, 2009, 46: 398-411.

[4] Gao M, Wang Y, Gao G Y. Dynamic response of cylindrical lined cavity in elastic medium[J]. Journal ofCentral South University, 2013, 20(10): 2849-2855.

[5] 姚署霖, 聞敏杰. 黏彈性土中襯砌隧道振動響應的解析解[J]. 工程力學, 2014, 31(3): 109-115.

YAO Shu-lin, WEN Min-jie.Analytical solutions of dynamic responses of lined tunnel in viscoelastic soil[J]. Engineering Mechanics, 2014, 31(3): 109-115.

[6] 高華喜, 聞敏杰. 內水壓力作用下黏彈性飽和土-隧洞襯砌相互作用[J]. 工程力學, 2013, 30(3): 289-296.

GAO Hua-xi, WEN Min-jie.Interaction of viscoelastic saturated soil and tunnel lining subject to inner water pressure[J]. Engineering Mechanics, 2013, 30(3): 289-296.

[7] 楊驍, 聞敏杰. 飽和分數導數型黏彈性土-深埋圓形隧洞襯砌系統的動力特性[J].工程力學,2012,29(12):248-255.

YANG Xiao, WEN Min-jie.Dynamic characteristics of saturated fractional derivative type viscoelastic soil and lining system with a deeply embedded circular tunnel[J]. Engineering Mechanics, 2012, 29(12): 248-255.

[8] 高盟, 高廣運, 李大勇. 考慮耦合質量影響的均布突加荷載作用下襯砌結構的瞬態響應[J]. 巖土工程學報, 2011, 33(6): 862-868.

GAO Meng, GAO Guang-yun, LI Da-yong.Transient response of lining structure subjected to sudden internal uniform loading considering effects of coupling mass[J]. Chinese Journal of Geotechnical Engineering, 2011, 33(6): 862-868.

[9] 高盟, 高廣運. 飽和地基任意形狀襯砌內部均布突加荷載作用下的動力響應[J]. 固體力學學報, 2011, 32(6): 655-664.

GAO Meng, GAO Guang-yun.Dynamic response of lining in saturated foundation subjected to sudden internal uniform loading[J]. Chinese Journal of Solid Mechanics,2011,32(6):655-664.

[10] 蔡袁強, 陳成振, 孫宏磊. 爆炸荷載作用下飽和土中隧道的瞬態動力響應[J]. 巖土工程學報,2011,33(3):361-367.

CAI Yuan-qiang, CHEN Cheng-zhen, SUN Hong-lei. Transient dynamic response of tunnels sujected to blast loads in saturated soil[J]. Chinese Journal of Geotechnical Engineering, 2011, 33(3): 361-367.

[11] 張智卿, 王奎華, 李強, 等. 非飽和土中端承樁縱向振動問題簡化解[J]. 工程力學, 2010, 27(5): 159-165.

ZHANG Zhi-qing, WANG Kui-hua, LI Qiang, et al. Simplified solution for vertical vibration of an end bearing pile embedded in unsturated soil[J]. Engineering Mechanics, 2010, 27(5): 159-166.

[12] 徐長節，馬曉華，蔡袁強. 彈性準飽和土中球空腔的動力響應[J]. 巖土工程學報，2005,27(8): 934-938.

XU Chang-jie, MA Xiao-hua, CAI Yuan-qiang. Dynamic response of spherical cavity in nearly saturated elastic soils[J]. Chinese Journal of Geotechnical Engineering, 2005,27(8): 934-938.

[13] 高華喜, 聞敏杰. 準飽和黏彈性土-隧道襯砌系統的動力特性[J]. 工程力學, 2013, 30(5): 90-96.

GAO Hua-xi, WEN Min-jie.Dynamic characteristics of nearly saturated viscoelastic soil and lining system with a tunnel[J]. Engineering Mechanics, 2013, 30(5): 90-96.

[14] 高華喜, 聞敏杰, 張斌. 具有圓形隧道的準飽和黏彈性土振動響應[J]. 浙江大學學報:工學版,2013,47(4):615-621.

GAO Hua-xi, WEN Min-jie, ZHANG Bin. Dynamic responses of nearly saturated viscoelastic soil with a deeply embedded circular tunnel[J]. Journal of Zhejiang University:Engineering Science, 2013, 47(4): 615-621.

[15] 沈珠江. 理論土力學[M]. 北京: 中國水利水電出版社, 2000.

Dynamic characteristics of an unsaturated soil and lining structure system with a deeply embedded circular tunnel

LIU Hong-bo1, LIU Jie1, WEN Min-jie2

(1. Zhengzhou Institute of Aeronautical Industry Management, Zhengzhou 450015, China; 2. Jiaxing Vocational Technical College, Jiaxing 314036, China)

Using the analytical method, the dynamic characteristics of an unsaturated soil and lining structure system with a deeply embedded circular tunnel were studied in frequency domain. Based on Bishop’s effective stress principle and one-phase flow consolidation theory, the motion equations of unsaturated soil subjected to one-phase flow were established. An analytical solution to the dynamic response of unsaturated soil subjected to harmonic load was obtained by means of the displacement potential function. The lining structure was treated as a homogeneous elastic medium. Based on the elastic theory, the steady-state dynamic response of a cylinder lining structure subjected to harmonic load was obtained. By utilizing the inner boundary condition of the lining structure and the continuity conditions between soil and lining, the unknown coefficients in the expressions were determined. Then, the effects of saturation, suction reduction coefficient and lining thickness on response amplitudes of unsaturated soil were investigated. It was shown that the response amplitudes of the system decrease gradually with increase in lining thickness; the saturation and suction reduction coefficient have a larger effect on the amplitudes of pore water press, while their effects on the amplitudes of radial displacement and hoop stress are smaller.

unsaturated soil; one-phase flow; cylinder lining structure; dynamic characteristics; analytical solution

10.13465/j.cnki.jvs.2016.13.007

國家自然科學基金(51508524)；河南省科技攻關項目(112102310633)；浙江省自然科學基金項目(LY12E09006)

2015-06-23修改稿收到日期：2015-09-17

劉洪波 男，博士，副教授，1975年生

聞敏杰 男，博士生，講師，1986年生

E-mail: 7279026@163.com

TU435

A