EI、Scopus 收录
中文核心期刊

成层半空间出平面自由波场的一维化时域算法

A 1-D time-domain method for 2-d wave motion in elastic layered half-space by antiplane wave oblique incidence

  • 摘要: 提出了一种计算出平面SH波斜入射时弹性水平成层半空间中自由波场时域计算的一维化有限元方法. 在进行有限元网格划分时,竖向单元取满足有限元模拟精度的任意尺寸,水平向网格尺寸由时间离散步长和水平视波速确定,并自动进行虚拟网格划分. 基底设置人工边界,并将波动输入转化为等效荷载施加在边界节点上. 然后将集中质量有限元法和中心差分法相结合建立节点运动方程,并将水平方向相邻节点的运动用该节点相邻时刻的运动表示,从而将求解节点运动的二维方程组转化为一维方程组. 求解此方程组,即得到自由场中竖向一列节点的运动. 最后根据行波传播的特点,可方便地确定全部自由波场. 理论分析和数值算例表明,该方法具有较高的精度和良好的稳定性.

     

    Abstract: A 1-D finite element method in time domain is developed inthis paper, which can be used to calculate the wave motion of free field inelastic layered half-space by antiplane SH wave oblique incidence. When thelayered half-space is discretized, the vertical element size is determinedconforming to the simulation accuracy; the horizontal element sizeis determined automatically by the horizontal apparent wave velocity and thediscrete time step in the step-by-step calculation, and then the elements aredivided virtually. Artificial boundary is constructed on the bottom ofthe computational area and the input wave motion is transformed into an equivalentload, which is applied on the nodes of the boundary. Then, the finite elementmethod with lumped mass and the central difference method are combined to establish the wave motion equations in 2-D finite element model.Since the displacement of any node in the finite element model can berepresented by that of the adjacent node in the horizontal direction, the 2-Dwave motion equations can be transformed into 1-D equations. By solving the1-D equations, the displacement of nodes in one vertical line can beobtained. Finally, the wave motion in elastic half-space is obtained basedon the propagation characteristics of traveling wave. Both the theoreticalanalysis and the numerical results demonstrate that the proposed methodfeatures high accuracy and good stability.

     

/

返回文章
返回