摘要
基于有限元分析软件ABAQUS平台,编制材料本构模型子程序UMAT,采用粘弹性人工边界条件模拟能量开放系统的边界,利用二次开发后的ABAQUS软件,模拟了隧道与围岩系统的三维地震反应,对比分析了不同入射角度下隧道结构与围岩土体的地震反应规律。分析结果表明,在地震荷载作用下,土体产生了不可恢复的塑性变形,主要表现为地表沉降;在斜入射条件下,随着入射角的增大,地表沉降显著增大,且沉降最终趋于稳定的时间有延长的趋势;均质地层中的大跨度长距离隧道结构,其横向地震反应受入射角的影响较大,表现为随着入射角的增大,结构的竖向沉降以及Mises应力显著增大;其纵向地震反应受入射角的影响不大。
Material constitutive model subroutine UMAT was programmed based on the platform of finite element analysis software ABAQUS.The viscoelastic artificial boundary was applied to simulate the boundary of the energy open system.The software ABAQUS developed secondly was used to simulate the three-dimensional seismic response of tunnel structures and surrounding system.Then the law of seismic response of tunnel structure and surrounding soil mass under different seismic incident angles was analyzed comparatively.Results show that the soil mass produces irrecoverable plastic deformation under the seismic load,which is mainly manifested as surface subsidence.The surface subsidence grows significantly and the subsidence eventually stabilizing time has a tendency to lengthen with the increase of incidence angle in conditions of oblique incidence.And the transverse seismic response of structure is greatly influenced by the incident angle for long span tunnel structures in homogeneous stratum.The larger the incidence angle,the greater the vertical subsidence and Mises stress of structure.Then the longitudinal seismic response is not much affected by the incident angle.
引文
[1]Kirzhner F,Rosenhouse G.Numerical analysis of tunnel dynamic response to earth motions[J].Tunnelling and Underground Space Technology,2000,15(3):249-258.
[2]许增会,宋宏伟,赵坚.地震对隧道围岩稳定性影响的数值模拟分析[J].中国矿业大学学报,2004,33(1):41-44.Xu Z H,Song H W,Zhao J.Numerical analysis of stability of tunnel surrounding rock under earthquake[J].Journal of China University of Mining Technology,2004,33(1):41-44.
[3]Ding J H,Jin X L,Guo Y Z,et al.Numerical simulation for large-scale seismic response analysis of immersed tunnel[J].Engineering Structures,2006,28(10):1367-1377.
[4]Anastasopoulos I,Gerolymos N,Drosos V,et al.Nonlinear response of deep immersed tunnel to strong seismic shaking[J].Journal of Geotechnical and Geoenvironmental Engineering,ASCE,2007,133(9):1067-1090.
[5]丁俊宏,金先龙,郭毅之,等.盾构隧道地震响应的三维数值模拟方法及应用[J].岩石力学与工程学报,2006,25(7):1430-1436.Ding J H,Jin X L,Guo Y Zh,et al.Three-dimensional numerical simulation method and its application to seismic response of shield tunnel[J].Chinese Journal of Rock Mechanics and Engineering,2006,25(7):1430-1436.
[6]川岛一彦.地下构造の物耐震设计[M].东京:鹿岛出版社,1994.
[7]杜修力,陈维,李亮,等.斜入射条件下地下结构时域地震反应分析初探[J].震灾防御技术,2007,2(3):290-296.Du X L,Chen W,Li L,et al.Preliminary study of time-domain seismic response for underground structures to obliquely incident seismic waves[J].Technology for Earthquake Disaster Prevention,2007,2(3):290-296.
[8]禹海涛,袁勇,顾玉亮,等.非一致激励下长距离输水隧道地震响应分析[J].水利学报,2013,44(6):718-725.Yu H T,Yuan Y,Gu Y L,et al.Effect of non-uniform excitation on seismic response of long-distance water-conveyance tunnel[J].Journal of Hydraulic Engineering,2013,44(6):718-725.
[9]Yao Y P,Hou W,Zhou A N.UH model:three-dimensional unified hardening model for overconsolidated clays[J].Geotechnique,2009,59(5):451-469.
[10]姚仰平,冯兴,黄祥,等.UH模型在有限元分析中的应用[J].岩土力学,2010,31(1):237-245.Yao Y P,Feng X,Huang X,et al.Application of UH model to finite element analysis[J].Rock and Soil Mechanics,2010,31(1):237-245.
[11]路德春.基于广义非线性强度理论的土的应力路径本构模型[D].北京:北京航空航天大学,2006.Lu D Ch.A constitutive model for soils considering complex stress paths based on the generalized nonlinear strength theory[D].Beijing:Beihang University,2006.
[12]Deeks A J,Randolph M F.Axisymmetric time-domain transmitting boundaries[J].Journal of Engineering Mechanics,1994,120(1):25-42.
[13]杜修力.局部解耦的时域波分析方法[J].世界地震工程,2000,16(3):22-26.Du X L.A partially decoupling analytical method for wave propagation problems in time domain[J].World Information on Earthquake Engineering,2000,16(3):22-26.
[14]刘晶波,王振宇,杜修力,等.波动问题中的三维时域粘弹性人工边界[J].工程力学,2005,22(6):46-51.Liu J B,Wang Zh Y,Du X L,et al.Three-dimensional visco-elastic artificial boundaries in time domain for wave motion problems[J].Engineering Mechanics,2005,22(6):46-51.
[15]杜修力,赵密,王进廷.近场波动模拟的一种应力人工边界[J].力学学报,2006,38(1):49-56.Du X L,Zhao M,Wang J T.A stress artificial boundary in FEA for near-field wave problem[J].Chinese Journal of Theoretical and Applied Mechanics,2006,38(1):49-56.
[16]王苏,路德春,杜修力.地下结构地震破坏静-动力耦合模拟研究[J].岩土力学,2012,33(11):3483-3488.Wang S,Lu D Ch,Du X L.Research on underground structure seismic damage using static-dynamic coupling simulation method[J].Rock and Soil Mechanics,2012,33(11):3483-3488.
[17]吕琦.隧道与围岩结构三维地震反应分析[D].北京:北京工业大学,2014.Lv Q.Three dimensional seismic response analysis of tunnels and surrounding rock structure[D].Beijing:Beijing University of Technology,2014.