成层饱和介质平面波斜入射问题的一维化时域方法

李伟华; 夏佩林; 张奎; 赵成刚

doi:10.6052/0459-1879-17-393

成层饱和介质平面波斜入射问题的一维化时域方法

doi: 10.6052/0459-1879-17-393

李伟华^1,*, ,,
夏佩林²,
张奎¹,
赵成刚¹

1 北京交通大学土木建筑工程学院, 北京 100044
2 中铁工程设计咨询集团有限公司, 北京 100055

基金项目: 国家自然科学基金项目(51378058,51778386)和国家重点基础研究发展计划项目(2015CB057800)资助.

详细信息

作者简介:
^*通讯作者：李伟华,副教授,主要研究方向：土动力学和地震工程. E-mail：whli@bjtu.edu.cn

通讯作者:
李伟华

中图分类号: TU45;
计量
- 文章访问数: 1171
- HTML全文浏览量: 179
- PDF下载量: 326
- 被引次数: 0
出版历程
- 收稿日期: 2017-11-22
- 刊出日期: 2018-03-18

ONE-DIMENSIONAL TIME-DOMAIN METHOD FOR FREE FIELD IN LAYERED SATURATED POROELASTIC MEDIA BY PLANE WAVE OBLIQUE INCIDENCE

1 School of Civil Engineering, Beijing Jiaotong University, Beijing 100044, China
2 China Railway Engineering Consulting Group Limited Company, Beijing 100055, China

摘要

摘要: 地震波斜入射下自由场的输入是大型结构抗震分析中亟待解决的问题之一,尤其是成层饱和多孔介质自由场问题,由于问题的复杂性,目前研究甚少. 本文基于Biot提出的饱和多孔介质动力方程,建立了一种新的求解平面波斜入射下基岩上覆饱和多孔介质成层场地自由场分析的一维化时域计算方法. 该方法首先根据Snell定律将饱和多孔介质二维空间问题转化为一维时域问题,通过对深度方向的有限元离散,得到饱和多孔介质波动问题的一维化有限元方程,然后采用单相弹性介质精确人工边界条件模拟基岩半空间的波动辐射和输入特征,通过考虑基岩与饱和多孔介质间透水或不透水边界条件以及不同饱和多孔介质交界面边界条件,形成基岩上覆成层饱和介质系统的整体有限元方程,最后采用中心差分法与Newmark平均加速度近似格式相结合的方法对时间进行离散,得到节点的动力时程的显式表达. 典型场地的地震反应分析表明,本文方法的计算结果与传递矩阵法结合傅里叶变换的计算结果完全吻合,证明了其有效性.
- 成层饱和多孔介质 /
- 自由场 /
- 斜入射 /
- Snell定律 /
- 时域有限元
Abstract: The input of free field under oblique incidence of seismic waves is one of the urgent problems to be solved in the seismic analysis of large structures. Because of the complexity of the problem, there are few studies on the free field of layered saturated poroelastic media at present. In this paper, a 1-D time-domain finite element method is proposed to simulate the plane wave motion in layered saturated poroelastic media overlaid on bedrock subjected to the oblique incidence plane wave. The method is on the basis of Biot dynamic theory for saturated poroelastic media. Firstly, the spatially 2-D problem is transformed into a 1-D time-domain problem along the vertical direction according to Snell’s law. 1-D finite element equations for poroelastic media are established by discretization principle and finite element. Then an exact artificial boundary condition for elastic media is used to model the wave absorption and input effects of the truncated bedrock half space. The global finite element equations for the system of layered saturated poroelastic media overlaid on bedrock are developed according to the drained or undrained boundary conditions between the poroelastic medium and the bedrock. By solving the 1D equations, the displacements of nodes in any vertical line can be obtained combining the method of central differences and Average acceleration of Newmark, and the wave motions in layered poroelastic medium system are finally determined based on the characteristic of traveling wave. The method is verified and by comparing with the frequency-domain transfer matrix method with fast Fourier transform in analyzing two engineering site.
- layered saturated poroelastic media /
- free field /
- oblique incident /
- Snell law /
- time-domain finite element

HTML全文

参考文献(1)

[1]	Jin X, Liao ZP.Statistical research on S-wave incident angle.Earthquake Research in China, 1994, 8(1): 121-131
[2]	Takahiro S, Takeshi U, Ryoichi T, et al.Estimation of earthquake motion incident angle at rock site. Proc of 12th World Conference Earthquake Engineering, NZ National Society for Earthquake Engineering. Auckland, New Zealand, 2000: 956
[3]	Wolf JP, Obernhuber P.Effects of horizontally traveling waves in soil-structure interaction.Nuclear Engineering and Design, 1979, 57(2): 221-244
[4]	李山有, 廖振鹏, 周正华. 大型结构地震反应数值模拟中的波动输入. 地震工程与工程振动, 2001, 6(2): 1-5
[4]	(Li Shanyou, Liao Zhenpeng, Zhou Zhenghua.Wave motion input in numerical simulation of seismic response for large-scale structure.Earthquake Engineering and Engineering Vibration, 2001, 6(2): 1-5 (in Chinese))
[5]	林皋. 地下结构地震响应的计算模型. 力学学报, 2017, 49(3): 528-542
[5]	(Lin Gao.A computational model for seismic response analysis of underground stauctures.Chinese Journal of Theoretical and Applied Mechanics, 2017, 49(3): 528-542 (in Chinese))
[6]	章小龙, 李小军, 陈国兴等. 黏弹性人工边界等效载荷计算的改进方法. 力学学报, 2016, 48(5): 1126-1135
[6]	(Zhang Xiaolong, Li Xiaojun, Chen Guoxing, Zhou Zhenghua.An improved method of the calculation of equivalent nodal forces in viscous-elastic artificial boundary.Chinese Journal of Theoretical and Applied Mechanics, 2016, 48(5): 1126-1135 (in Chinese))
[7]	Thomson WT.Transmission of elastic waves through a stratified solid medium.Journal of Applied Physics, 1950, 21(1): 89-93
[8]	Kausel E, Roësset JM.Stiffness matrices for layered soils.Bulletin of the Seismological Society of America, 1981, 71(6): 1743-1761
[9]	尤红兵, 赵凤新, 荣棉水. 地震波斜入射时水平层状场地的非线性地震反应. 岩土工程学报, 2009, 31(2): 234-240
[9]	(You Hongbing, Zhao Fengxin, Rong Mianshui.Nonlinear seismic response of horizontal layered site due to inclined wave.Chinese Journal of Geotechnical Engineering, 2009, 31(2): 234-240 (in Chinese))
[10]	王笃国, 赵成刚. 地震波斜入射时二维成层介质自由场求解的等效线性化方法. 岩土工程学报, 2016, 38(3): 554-561
[10]	(Wang Duguo, Zhao Chenggang.Two-dimensional equivalent linear seismic analysis of free field in layered half-space due to oblique incidence.Chinese Journal of Geotechnical Engineering, 2016, 38(3): 554-561 (in Chinese))
[11]	廖河山, 陈清军. SH波倾斜入射时土层的非线性响应分析. 同济大学学报(自然科学版), 1994, 22(4): 517-522
[11]	(Liao Heshan, Chen Qingjun.Nonlinear responses of layered soils to obliquely incident SH waves. Journal of Tongji University( Natural Science), 1994, 22(4): 517-522 (in Chinese))
[12]	李山有, 王学良, 周正华. 地震波斜入射情形下水平成层半空间自由场的时域计算. 吉林大学学报(地球科学版), 2003, 33(3): 372-375
[12]	(Li Shanyou, Wang Xueliang, Zhou Zhenghua.The time-step numerical simulation of free field motion of layered half-space for inclined seismic waves.Journal of Jilin University ( Earth Science Edition), 2003, 33(3): 372-375 (in Chinese))
[13]	刘晶波, 王艳. 成层半空间出平面自由波场的一维化时域算法. 力学学报, 2006, 38(2): 219-225
[13]	(Liu Jingbo, Wang Yan.A 1-D time-domain method for 2-D wave motion in elastic layered half-space by antiplane wave oblique incidence.Chinese Journal of Theoretical and Applied Mechanics, 2006, 38(2): 219-225 (in Chinese))
[14]	刘晶波, 王艳. 成层介质中平面内自由波场的一维化时域算法. 工程力学, 2007, 24(7): 16-22
[14]	(Liu Jingbo, Wang Yan.1D time-domain method for in-plane wave motion of free field in layered media.Engineering Mechanics, 2007, 24(7): 16-22 (in Chinese))
[15]	Liu Jingbo, Wang Yan.A 1D time-domain method for in-plane wave motions in a layered half-space. Acta Mechanica Sinica, 2007, 23(6): 673-680
[16]	赵密, 杜修力, 刘晶波等. P-SV 波斜入射时成层半空间自由场的时域算法. 地震工程学报, 2013, 35(1): 84-90
[16]	(Zhao Mi, Du Xiuli, Liu Jingbo, et al.Time-domain method for free field in layered half space under P-SV Waves of oblique incidence.China Earthquake Engineering Journal, 2013, 35(1): 84-90 (in Chinese))
[17]	卓卫东, 高智能, 谷音. P-SV 波斜入射时有阻尼成层介质自由波场的一维化时域算法. 水利与建筑工程学报, 2016, 14(6): 18-24, 34
[17]	(Zhuo Weidong, Gao Zhineng, Gu Yin.A 1D time-domain method for in-plane wave motion of free field in layered media with damping under obliquely incident P-SV waves.Journal of Water Resources and Architectural Engineering, 2016, 14(6): 18-24, 34 (in Chinese))
[18]	Zhao M, Yin HQ, Du XL, et al.1D finite element artificial boundary method for layered half space site response from obliquely incident earthquake.Earthquakes and Structures, 2015, 9(1): 173-194
[19]	Biot MA.Theory of propagation of elastic wave in fluid-saturated porous soil.The Journal of the Acoustical Society of America, 1956, 28(2): 168-178
[20]	Schanz M. Poroelastodynamics: linear models, analytical solutions, and numerical methods. Applied Mechanics Reviews, 2009, 62(3): 030803-1-15
[21]	王子辉. 饱和两相与单相土互层场地中地铁车站地震反应分析. [博士论文]. 北京:北京交通大学, 2008
[21]	(Wang Zihui.Seismic response analysis of subway station in saturated phase and single phase soil. [PhD Thesis]. Beijing: Beijing Jiaotong University, 2008 (in Chinese))
[22]	赵成刚, 王进廷, 史培新等. 流体饱和两相多孔介质动力反应分析的显式有限元法. 岩土工程学报, 2001, 23(2): 178-182
[22]	(Zhao Chenggang, Wang Jinting, Shi Peixin, et al.Dynamic analysis of fluid-saturated porous media by using explicit finite element method.Chinese Journal of Geotechnical Engineering, 2001, 23(2): 178-182 (in Chinese))
[23]	Biot MA, Willis DG.The elastic coefficients of the theory of consolidation.Journal of Applied Mechanics, 1957, 15(2): 594-601
[24]	周正华, 周扣华. 有阻尼振动方程常用显式积分格式稳定性分析. 地震工程与工程振动, 2001, 21(3): 22-28
[24]	(Zhou Zhenghua, Zhou Kouhua.Stability analysis of an integral method for damped vibration equation.Earthquake Engineering and Engineering Vibration, 2001, 21(3): 22-28 (in Chinese))
[25]	王进廷, 张楚汉, 金峰. 有阻尼动力方程显式积分方法的精度研究. 工程力学, 2006, 23(3): 1-5
[25]	(Wang Jinting, Zhang Chuhan, Jin Feng.On the accuracy of several explicit integration schemes for dynamic equation with damping.Engineering Mechanics, 2006, 23(3): 1-5 (in Chinese))
[26]	李亮, 杜俢力, 李立云等. 两相介质波动问题显式有限元方法稳定性研究. 西北地震学报, 2011, 33(3): 218-222, 227
[26]	(Li Liang, Du Xiuli, Li Liyun, et al.Study on stability of explicit finite element method for wave motion of fluid-saturated porous media.Northwestern Seismological Journal, 2011, 33(3): 218-222, 227 (in Chinese))
[27]	Lin CH, Lee VW, Trifunac MD.The reflection of plane waves in a poroelastic half-space saturated with inviscid fluid.Soil Dynamics and Earthquake Engineering, 2005. 25(13): 205-223
[28]	傅淑芳, 刘宝诚. 地震学教程. 北京: 地震出版社, 1991: 73-84
[28]	(Fu Shufang, Liu Baocheng. Seismology Tutorial.Beijing: Seismic Press, 1991: 73-84 (in Chinese))
[29]	王艳. 非一致地震动场数值方法研究及在结构动力分析中的应用. [博士论文]. 北京: 清华大学, 2007
[29]	(Wang Yan.Research on the numerical method for asynchronous seismic wave motions and its application in dynamic analysis of atructures. [PhD Thesis]. Beijing: Tsinghua University, 2007 (in Chinese))
[30]	李伟华, 赵成刚. 地下水位变化对地震地面运动的影响. 地震学报, 2015, 37(3): 482-492
[30]	(Li Weihua, Zhao Chenggang.Effects of the groundwater level variation on earthquake ground motions.Acta Seismologica Sinica, 2015, 37(3): 482-492 (in Chinese))