含方形凹陷半无限非均匀介质波动问题FDM模拟  

FDM Simulation of Wave Motion in a Semi-infinite Inhomogeneous Medium with a Rectangular Depression

作　　者：杨在林[1] 孙铖[1] 王耀[1] 李志东[1] 

机构地区：[1]哈尔滨工程大学航天与建筑工程学院,黑龙江哈尔滨150001

出　　处：《地震工程学报》2015年第2期553-558,共6页China Earthquake Engineering Journal

基　　金：2015年地震行业科研专项经费项目(201508026-02);黑龙江省自然科学基金(A201310);黑龙江省博士后科研启动金(LBH-Q13040)

摘　　要：采用规则网格有限差分方法对二维平面弹性波动方程进行差分离散,得到相应的弹性波动方程的有限差分方程,再将弹性波动方程的差分格式与吸收边界、自由边界的离散形式结合形成弹性波动方程有限差分方程解决问题的主体,将其应用于含方形凹陷半无限非均匀介质的模型中进行数值模拟,得到此离散化模型中不同时刻不同节点的位移值。针对具体算例,运用上述方法结合科学计算软件MATLAB和结果后处理软件DIFEM ISOLINE PLOTER得到不同时刻的水平方向位移等值线图与接收器测量点处的合成位移记录,讨论非均匀介质、吸收边界、方形凹陷等对波动特性的影响。Using the finite difference method,the two-dimensional elastic wave equation is processed to obtain the finite difference equation.Then,the differential elastic wave equation format is combined with the absorbing boundary and the discrete form of the free boundary.The elastic wave equation finite difference equation is formed to solve the problem.Using the finite difference equation,the displacement on different nodes can be obtained at different times for the rectangleshaped depression of a semi-infinite inhomogeneous medium.Using the finite difference equation of the two-dimensional elastic wave equation combined with the scientific computation software,MATLAB,and the post-processor software,DIFEM ISOLINE PLOTER,the displacement isoline of the rectangle-shaped depression of the semi-infinite inhomogeneous medium in the xdirection can be obtained for different times.The effect of the inhomogeneous media,absorbing boundary,and the rectangle-shaped depression on the wave characteristics is analyzed numerically.

关 键 词：非均匀介质 波动 有限差分方法 凹陷 Higdon吸收边界 

分 类 号：O343.7[理学—固体力学]