P波斜入射角对沉管隧道地震响应的影响1

陆耀波 崔 杰 渠建新 袁 杰



P波斜入射角对沉管隧道地震响应的影响1

陆耀波 崔 杰 渠建新 袁 杰

（广州大学，土木工程学院，广州 510006）

为了研究P波斜入射对沉管隧道地震响应的影响，以港珠澳大桥沉管隧道为工程背景，考虑上覆海水与海床、沉管隧道之间耦合作用，采用粘弹性边界和等效力的地震荷载输入方式，利用ADINA软件建立三维有限元模型进行地震响应分析。分析入射角为0°、20°、40°、50°、60°时P波对沉管隧道环向应力峰值（正应力峰值、剪应力峰值）和位移峰值的影响，结果表明：入射角为40°时，沉管隧道应力峰值最大；入射角为0°—40°时，隧道的应力峰值逐渐增大，入射角为40°—60°时，隧道的应力峰值逐渐减小；隧道截面4个转角处及隔墙与顶板、底板的连接处为隧道剪应力峰值最大处；隧道截面左侧剪应力峰值远大于右侧；隧道顶板正应力峰值最大，顶板的正应力峰值大约为底板的2倍；隧道截面左侧位移峰值远大于隧道截面右侧。

沉管隧道 P波斜入射 地震响应分析 流固耦合

引言

地铁、隧道等地下结构破坏会引起重大的安全事故及经济损失，因此人们对隧道等地下结构的安全性十分重视，针对地下结构动力分析的研究已有不少成果（韩大建等，1999；Uenishi等，2000；Choi等，2002；Huo等，2003；袁勇等，2016）。目前，研究沉管隧道纵向地震响应的手段有：①大型振动台试验和离心机试验（袁勇等，2016）。利用振动台等试验方法分析沉管隧道的地震响应具有较高的可信度，但一般用来模拟水平和竖向振动。②简化力学模型方法：梁-弹簧模型、质点-弹簧模型（或等效质点-弹簧模型），前者把隧道作为支撑在弹性地基上的连续梁，后者把隧道作为多质点系统，接头统一用弹簧模拟，用弹簧和阻尼来模拟土-结构之间的相互作用（韩大建等，1999）。这些简化方法基于众多假设条件，所以计算精确度有所欠缺。③采用ADINA、ABUQUS、ANSYS等软件进行仿真模拟。利用大型软件进行水域沉管隧道地震响应分析已得到广泛的应用，陈向红（2013）利用有限元软件分析竖向地震作用下，动水压力对浅埋海底隧道内力的影响；陈贵红（2002）分析了在竖向地震作用下，隧道的埋深、海床的地质条件等对沉管隧道的影响；李鹏（2013）利用ANSYS软件进行沉管隧道纵向地震响应分析，但并未考虑海水的作用；周晓洁等（2017）利用ABAQUS软件建立二维模型分析SV波斜入射沉管隧道地震响应。

本文利用ADINA软件，考虑地基的人工边界和海水-海床-隧道间的相互作用，建立三维沉管隧道模型，分析P波的入射角度对沉管隧道动力响应的影响。

1 流固耦合分析原理

本文所用流体为理想流体，其连续方程和Naiver-Stokes控制方程可表达为：

根据流体的边界条件和流体与海床的耦合边界条件，利用泛函对流体的连续方程进行有限元离散，可以得到流体的振动控制方程（陈贵红，2002）：

考虑海床与流体的耦合边界条件，同样对海床土与隧道进行有限元离散，可得其振动控制方程（陈贵红，2002）：

因此，存在外部荷载时，流-固耦合的振动方程可以改写为：

式中R为固态受到除流体压力之外的荷载（N）。

目前，求解公式（5）的方法有迭代法（分离法）和直接耦合法（同步求解法），具体的求解步骤可参考文献（岳戈，2010）。

2 粘弹性人工边界

粘弹性人工边界由弹簧和阻尼器组成，由于忽略质量，并把阻尼器的一端固定，原来的阻尼和刚度系数的精度有所降低，刘晶波等（2006）对粘弹性人工边界的参数进行了修正并给出了其取值范围。

切向边界的弹簧系数和阻尼系数分别为：

法向边界为：

3 等效荷载输入

在近场波动分析中，由于使用人工边界来模拟无限域地基，当外部荷载从远场无限域进入近场有限域时，在人工边界处，需要考虑荷载的输入问题，即外源输入问题。目前处理外源输入问题主要有波场分解法和等效边界力法。人工边界处的总波场可分为内行场（下标r表示）和外行场（下标s表示），人工边界处位移和作用力可以分解为（赵密，2009）：

式中：uBr和uBs分别表示内行场位移和外行场位移（m）；FBr和FBs分别是保证边界位移等于内行场所需要提供的抵抗人工边界条件的荷载和保证边界位移等于内行场所需要提供的抵抗有限域的荷载（N）；其中FBs可以表示为：

下面推导斜入射的FBr公式。本文只考虑平面内的斜入射，推导也只适用于平面斜入射。假设P波以角从左侧进入有限域。

图1 P波斜入射示意图

利用文献（赵密，2009）中的方法推导出在边界处所施加的等效荷载。

左边界：

前边界上的等效荷载表达为：

后边界上的等效荷载可表达为：

底边界上的等效荷载可表达为：

图2 输入位移和地面位移的对比

4 沉管隧道三维计算模型

4.1 工程背景

本文建立三维分析模型，分析P波与隧道横截面形成的斜入射角度对沉管隧道的影响参考港珠澳大桥沉管隧道的尺寸，高12m，宽38m，采用3个节段进行分析，每个节段长22m。计算条件为：横向长度118m，高度78m，隧道上覆土高度6m，水深10m。隧道和土体的参数见表1。接头的力学曲线（刘鹏等，2014）见图3、4，输入时程为0.1g的El Centro波。

沉管隧道三维计算模型以三维八节点六面体划分，共74532个单元，59256个节点。土体和隧道采用实体单元，海水采用势流体单元。土体底部和侧面均为粘弹性边界，海水与土体接触边界为流固耦合边界（Fluid-Structure），海水侧面采用无限域边界（Fluid infinite region），海水表面采用自由表面（Free surface）。

4.2 结果及分析

从结果（图8—14）可以得知：入射角为40°时，隧道应力峰值（包括剪应力峰值、正应力峰值）最大；入射角为0°—40°时，隧道的应力峰值逐渐增大，入射角为40°—60°时，隧道的应力峰值逐渐减小。

下面分别讨论顶板、底板和侧向板以及隔墙的应力峰值。

根据剪应力峰值图（图8）发现，隧道截面的左侧（入射侧）剪应力峰值远大于隧道截面的右侧，隧道顶板（底板）和隔墙交接处以及截面4个角处及附近剪应力峰值相对较大。正应力峰值图（图9）显示隧道顶板正应力峰值最大，其它地方正应力峰值相对较小，顶板的正应力峰值大约为底板的2倍。根据位移峰值图（图10）可知，左侧板的位移峰值比底板和顶板大。

图5 沉管隧道有限元模型

图6 中间管节截面示意图

图7 El Centro波加速度曲线

图8 中间管节截面剪应力峰值图

根据图11—14可知，左侧隔墙比右侧隔墙的应力峰值大，其中左侧为地震波斜入射一侧；隔墙中间剪应力峰值比两端大；隔墙上端正应力峰值比下端大。

图9 中间管节截面正应力峰值图

图10 中间管节截面竖向位移峰值图

图11 中间管节隔墙截面剪应力峰值图

图12 中间管节隔墙截面正应力峰值图

图13 中间管节隔墙截面水平位移峰值图

图14 中间管节隔墙截面竖向位移峰值图

5 结论

本文采用粘弹性边界和地震波以等效荷载输入的方法，利用ADINA软件建立三维模型，考虑海水和隧道、海床的流固耦合，分析P波斜入射对隧道截面内力的影响，得到的结论可为沉管隧道设计提供参考。

（1）P波斜入射角为0°—40°时，隧道的应力峰值逐渐增大，入射角为40°—60°时，隧道的应力峰值逐渐减小。

（2）隧道截面4个转角处及隔墙与顶板、底板的连接处的剪应力最大，所以隧道4角附近及隔墙与顶板、底板的连接处为薄弱位置。

（3）隧道顶板正应力最大，约为底板的2倍；隧道截面左侧位移远大于隧道截面右侧。

陈贵红，2002．沉管隧道抗震数值分析．成都：西南交通大学．

陈向红，2013．大型水下隧道与附属竖井的地震响应研究．北京：北京交通大学．

韩大建，周阿兴，1999．沉管隧道地震响应分析的等效质点系模型探讨．华南理工大学学报（自然科学版），27（11）：108—114．

李鹏，2013．饱和地基中隧道纵向地震反应的数值分析．北京：清华大学．

廖振鹏，2002．工程波动理论导论．2版．北京：科学出版社．

刘晶波，谷音，杜义欣，2006．一致粘弹性人工边界及粘弹性边界单元．岩土工程学报，28（9）：1070—1075．

刘鹏，丁文其，金跃郎等，2014．沉管隧道接头三维非线性刚度力学模型．同济大学学报（自然科学版），42（2）：232—237．

袁勇，禹海涛，燕晓等，2016．超长沉管隧道多点振动台试验模拟与分析．中国公路学报，29（12）：157—165．

岳戈，2010．ADINA流体与流固耦合功能的高级应用．北京：人民交通出版社．

赵密，2009．近场波动有限元模拟的应力型时域人工边界条件及其应用．北京：北京工业大学．

周晓洁，张越宇，何颖等，2017．地震SV波斜入射下沉管隧道的地震响应分析．地震工程学报，39（4）：600—608．

Choi J. S., Lee J. S., Kim J. M., 2002. Nonlinear earthquake response analysis of 2-D underground structures with soil-structure interaction including separation and sliding at interface. In: Proceedings of the 15th ASCE Engineering Mechanics Conference. New York: Columbia, 1—8.

Huo H. B., Bobet A., 2003. Seismic design of cut and cover rectangular tunnels-evaluation of observed behavior of Dakai station during Kobe earthquake. In: Proceedings of 1st World Forum of Chinese Scholars in Geotechnical Engineering. Shanghai: Tongji University Press, 456—466.

Uenishi K., Sakurai S., 2000. Characteristic of the vertical seismic waves associated with the 1995 Hyogo-ken Nanbu (Kobe), Japan earthquake estimated from the failure of the Daikai Underground Station. Earthquake Engineering and Structural Dynamics, 29(6): 813—821.

Analysis on Seismic Responses of Immersed Tunnel Under Inclined P Waves

Lu Yaobo, Cui Jie, Qu Jianxin and Yuan Jie

(School of Civil Engineering, Guangzhou University, Guangzhou 510006, China)

To study the seismic response of immersed tunnel under inclined P waves, we carried out a series of simulation analysis. In this paper, we presented 3D analytical models of immersed tunnel which is based on the Hong Kong-Zhuhai-Macao Bridge immersed tunnel taking account of viscous-spring artificial boundary and solid-fluid interaction modeled by ADINA, and analyzed the seismic response of immersed tunnel under P waves with angle of incidence 0°, 20°, 40°, 50°, 60°. It is shown that peak value of stress of tunnel under P waves with angle of incidence 40° reached maximum. The peak value of stress of tunnel had shown a magnifying trend when P-wave angle of incidence 0° increased to 40° and a decreasing trend when P-wave angle of incidence 40° increased to 60°. The larger peak stress point on the immersed tunnel takes place in the four corners of the tunnel and the connection between the partition wall and the bottom plate and the roof. The peak value of shear stress of left side of tunnel cross section is larger than the right side of the tunnel section. The largest peak value of normal stress of tunnel appears on roof, and the peak value of normal stress of roof is about 2 times of on the floor. The peak value of displacement of left side of tunnel cross section is larger than the right side of the tunnel section.

Immersed tunnel; Inclined P waves; Earthquake response; Tunnel-soil-fluid interaction

陆耀波，崔杰，渠建新，袁杰，2018．P波斜入射角对沉管隧道地震响应的影响．震灾防御技术，13（3）：562—570．

10.11899/zzfy20180307

国家重点研发计划（2016YFC0800205），国家自然科学基金面上项目（51578167），国家自然科学基金重点项目（51438004）

2017-12-01

陆耀波，男，生于1987年。在读博士研究生。主要从事地下结构抗震研究。E-mail：luyaobogz@163.com