摘要
在显式有限元方法结合黏弹性人工边界的时域波动方法的基础上,建立了地震波垂直输入时的一种简化输入方法。将近场有限元模型沿高度方向进行分层,将地震动的入射运动转化为作用于人工边界底面及分层后每层侧面上的均布力,以实现地震动的输入。与以等效节点力的方式实现的地震动输入相比,施加均布力的方式简化了地震波输入的前处理工作,且又能保证与等效节点力方式相同的精度。自由场数值算例表明:当该方法中的分层高度与波动有限元网格离散要求的最大尺寸相等时,该方法与等效节点力方法具有相同的精度;当局部区域按网格离散要求的最大尺寸进行分层而在其它区域放大分层高度时,局部区域上的近场波动响应仍可保证具有相当高的精度。另外,某隧道结构地震响应算例的计算结果同样说明了该方法的有效性。
On the basis of the time-domain wave method combined with the explicit finite element method and the viscous-elastic artificial boundary condition,a simplified seismic wave input method for near-field wave analysis was proposed. With this approach,the near-field model was stratified in the height direction,and then the seismic motion was converted into stresses acted on the viscous-elastic artificial boundary to realize seismic wave input. Compared with the equivalent node force input method,the new method simplified the antecedent processing with an acceptable accuracy.The numerical results obtained from a free ground model indicated that when the stratified size equals the maximum size of the element mesh in dynamical analysis,the accuracy of the simplified method is the same as that of the equivalent node method; when the local parts of the near-field model were stratified with the maximum element mesn size of dynamical analysis while the other parts were stratified with a lager size,the dynamic responses of the local parts have a higher accuracy. Moreover,the dynamic response computation results of a rock tunnel demonstrated the effectiveness of the proposed simplified approach.
引文
[1]杜修力.工程波动理论与方法[M].北京:科学出版社,2009:15-21.
[2]Lysmer J,Kuhlemeyer R L.Finite dynamic model forinfinite media[J].Journal of the Engineering Mechanics Division ASCE 95(EM4),1969:869-877.
[3]廖振鹏.工程波动理论导论[M].北京:科学出版社,2002.
[4]Deeks A J,Randolph M F.Axisymmetric time-domain transmitting boundaries[J].Journal of Engineering Mechanics,ASCE,1994,120(1):25-42.
[5]Joyner W B,Chen A T F.Calculation of nonlinear ground response in earthquake[J].Bull.Seism.Soc.Amer.,1975,65(5):1315-1336.
[6]刘晶波,吕彦东.结构-地基动力相互作用问题分析的一种直接方法[J].土木工程学报,1998,31(3):55-64.LIU Jing-bo,LU Yan-dong.A direct method for analysis of dynamic soil-structure interaction[J].China Civil Engineering Journal,1998,31(3):55-64.
[7]王振宇,刘晶波.成层地基非线性波动问题人工边界与波动输入研究[J].岩石力学与工程学报,2004,23(7):1169-1173.WANG Zhen-yu,LIU Jing-bo.Study of wave motion input and artificial boundary for problem of nonlinear wave motion in layered soil[J].Chinese Journal of Rock Mechanics and Engineering,2004,23(7):1169-1173.
[8]杜修力,赵密.基于黏弹性边界的拱坝地震反应分析方法[J].水利学报,2006,37(9):1063-1069.DU Xiu-li,ZHAO Mi.Analysis method for seismic response of arch dams in time domain based on viscous-spring artificial boundary conditions[J].Shuli Xuebao,2006,37(9):1063-1069.
[9]张顺福,王海波,杨会臣,等效黏性边界单元及黏性边界波动输入方法[J].水利学报,2008,39(10):1248-1255.ZHANG Shun-fu,WANG Hai-bo,YANG Hui-chen.Equivalent viscous boundary elements and the method of wave input for viscous boundary[J].Shuli Xuebao,2008,39(10):1248-1255.
[10]张波,李术才,杨学英,等.三维黏弹性人工边界地震波动输入方法研究[J].岩土力学,2009,30(3):774-778.ZHANG Bo,LI Shu-cai,YANG Xue-ying,et al.Research on seismic wave input with three-dimensional viscoelastic artificial boundary[J].Rock and Soil Mechanics,2009,30(3):774-778.
[11]赵武胜,陈卫忠,黄胜,等.地下结构抗震分析中地震动输入方法研究[J].土木工程学报,2010,43(S):554-559.ZHAO Wu-sheng,CHEN Wei-zhong,HUANG Sheng,et al.Research on method of seismic motion input in a seismic analysis for underground structure[J].China Civil Engineering Journal,2010,43(S):554-559.
[12]刘晶波,王振宇,杜修力,等.波动问题中的三维时域粘弹性人工边界[J].工程力学,2005,22(6):46-51.LIU Jing-bo,WANG Zhen-yu,DU Xiu-li,et al.Threedimensional viscos-elastic artificial boundaries in time domain for wave motion problems[J].Engineering Mechanics,2005,22(6):46-51.