Hybrid absorbing boundary condition based on transmitting boundary and its application in 3D fractional viscoacoustic modeling

An accurate numerical simulation for wave equations is essential for understanding of wave propagation in the earth's interior as well as full waveform inversion and reverse time migration. However, due to computational cost and hardware capability limitations, numerical simulations are often performed within a finite domain. Thus, an adequate absorbing boundary condition (ABC) is indispensable for obtaining accurate numerical simulation results. In this study, we develop a hybrid ABC based on a transmitting boundary, which is referred to as THABC, to eliminate artificial boundary reflections in 3D second-order fractional viscoacoustic numerical simulations. Furthermore, we propose an adaptive weighted coefficient to reconcile the transmitting and viscoacoustic wavefields in THABC. Through several numerical examples, we determine that the proposed THABC approach is characterized by the following benefits. First, with the same number of absorbing layers, THABC exhibits a better ability in eliminating boundary reflection than traditional ABC schemes. Second, THABC is more effective in computation, since it only requires the wavefields at the current and last time steps to solve the transmitting formula within the absorbing layers. Benefiting from a simple but effective combination between the transmitting equation and the second-order wave equation, our scheme performs well in the 3D fractional Laplacian viscoacoustic numerical simulation.

FINITE ELEMENT-ARTIFICIAL TRANSMITTING BOUNDARY METHOD FOR WAVE SCATTERING FROM IRREGULAR CYLINDER

The finite element artificial transmitting boundary method is employed here to treat the near field scattering of a cylindrical wave from an irregular cylinder. A comparison is made between this method and the analytical one. And then examples are given to demonstrate the solution of several problems of the irregular object scattering. The method can not only produce clear physical pictures, but can efficiently handle many complicated scattering problems.

Time-Domain Analysis of Underground Station-Layered Soil Interaction Based on High-Order Doubly Asymptotic Transmitting Boundary

Based on the modified scale boundary finite element method and continued fraction solution,a high-order doubly asymptotic transmitting boundary(DATB)is derived and extended to the simulation of vector wave propagation in complex layered soils.The high-order DATB converges rapidly to the exact solution throughout the entire frequency range and its formulation is local in the time domain,possessing high accuracy and good efficiency.Combining with finite element method,a coupled model is constructed for time-domain analysis of underground station-layered soil interaction.The coupled model is divided into the near and far field by the truncated boundary,of which the near field is modelled by FEM while the far field is modelled by the high-order DATB.The coupled model is implemented in an open source finite element software,OpenSees,in which the DATB is employed as a super element.Numerical examples demonstrate that results of the coupled model are stable,accurate and efficient compared with those of the extended mesh model and the viscous-spring boundary model.Besides,it has also shown the fitness for long-time seismic response analysis of underground station-layered soil interaction.Therefore,it is believed that the coupled model could provide a new approach for seismic analysis of underground station-layered soil interaction and could be further developed for engineering.

Comparison of perfectly matched layer and multi-transmitting formula artificial boundary condition based on hybrid finite element formulation

The theory of perfectly matched layer （PML） artificial boundary condition （ABC）, which is characterized by absorption any wave motions with arbitrary frequency and arbitrarily incident angle, is introduced. The construction process of PML boundary based on elastodynamic partial differential equation （PDE） system is developed. Combining with velocity-stress hybrid finite element formulation, the applicability of PML boundary is investigated and the numerical reflection of PML boundary is estimated. The reflectivity of PML and multi-transmitting formula （MTF） boundary is then compared based on body wave and surface wave simulations. The results show that although PML boundary yields some reflection, its absorption performance is superior to MTF boundary in the numerical simulations of near-fault wave propagation, especially in comer and large angle grazing incidence situations. The PML boundary does not arise any unstable phenomenon and the stability of PML boundary is better than MTF boundary in hybrid finite element method. For a specified problem and analysis tolerance, the computational efficiency of PML boundary is only a little lower than MTF boundary.

Analysis of high-frequency local coupling instability induced by multi- transmitting formula - P-SV wave simulation in a 2D waveguide

Local coupling instability will occur when the numerical scheme of absorbing boundary condition and that of the field wave equation allow energies to spontaneously enter into the computational domain. That is, the two schemes support common wave solutions with group velocity pointed into the computation domain. The key to eliminate local coupling instability is to avoid such wave solutions. For lumped-mass finite element simulation of P-SV wave motion in a 2D waveguide, an approach for stable implementation of high order multi-transmitting formula is provided. With a uniform rectangular mesh, it is proven and validated that high-freqaency local coupling instability can be eliminated by setting the ratio of the element size equal to or greater than x/2 times the ratio of the P wave velocity to the S wave velocity. These results can be valuable for dealing instability problems induced by other absorbing boundary conditions.

Dynamic characterization of tailing dam using fully coupled dynamic analysis with different boundary conditions-a case study

The seismic response analysis of a tailing dam is studied using a fully coupled effective stress approach in conjunction with an advanced multi yield surface plastic constitutive model for tailing material.Strain controlled static and cyclic triaxial tests were carried out to obtain the constitutive model for the tailing material.The tailing materials were collected from the Rampura Agucha tailing dam(Rajasthan State,India).A 2D nonlinear finite element(FE)model was then developed using different boundary conditions from the tailing embankment constructed using the downstream and upstream method of rising using OpenSees software.In first case,the model boundary was fixed in both the X and Y directions,and in the second case,viscous dashpots were introduced for both side and horizontal boundaries.The model was validated with experimental results on tailing material.Analyses were carried out considering five different earthquake motions,which were applied at the base.Comparisons of the different boundary conditions in terms of displacement flow vectors,pore pressure and stress-strain curves during shaking are presented.From the analysis,it was observed that the viscous boundary condition replicates the actual field conditions more accurately than the fixed boundary condition.In addition,it was found that the tailing embankment constructed by the downstream and upstream method of rising is not susceptible to liquefaction and lateral spreading for earthquake motions,even for a magnitude>5.5.

Transmitted and Reflected Coefficients for Horizontal or Vertical Plate Type Breakwater

Surface or submerged horizontal or vertical plate can be considered as a new concept breakwater. This paper investigates the wave-plate interaction of this type of breakwater by use of the boundary element method. The relationships of wave transmitted and reflected among plate thickness, submergence and length are carefully studied by numerical simulation. It is shown that： （1） The transmitted coefficients of submerged horizontal plate or vertical plate will become larger with the increase of plate thickness and reduce rapidly with the decrease of plate submergence. （2） Both surface horizontal and vertical plate are efficient for intermediate and short wave elimination, but vertical plate is more effective. （3） Submerged horizontal plate can act more effectively than submerged vertical plate does. With all wave frequencies, the vertical plate almost has no wave elimination effect.

Identifying boundary between near field and far field in ground vibration caused by surface loading

The boundary between the near and far fields is generally defined as the distance from the vibration source beyond which ground vibrations are mainly dominated by Rayleigh waves. It is closely related to the type of vibration source and the soil properties. Based on the solutions of the Lamb＇s problem, the boundary at the surface between the near and far fields of ground vibration was investigated for a harmonic vertical concentrated load and an infinite line load at the surface of a visco-elastic half-space. Particularly, the variation of the boundary with the material damping was investigated for both cases. The results indicate that the material damping slightly contributes to the attenuation of vibrations in the near-source region, but significantly reduces the vibrations in the region that is at some distance away from the source. When taking the material damping into consideration, the boundary between the near and far fields tends to move towards the vibration source. Compared with the vibrations caused by a concentrated load, the vibrations induced by an infinite line load can affect a larger range of the surrounding environment, and they attenuate more slowly. This means the boundary between the near field and far field should move fitrther away from the source. Finally, the boundaries are defined in terms of R-wave length （2R） and Poisson ratio of the ground （o）. For the case of a point load, the boundary is located at the distance of （5.0-6.0）2R for v≤0.30 and at the distance of （2.0--3.0）2R for v≥0.35. For the case of an infinite line load, the boundary is located at the distance （5.5-6.5）2rt for v≤0.30 and at the distance （2.5--3.5）2R for v≥0.35.

New artificial boundary condition for saturated soil foundations

Anew artificial boundary model based on multi-directional transmitting and viscous-spring artificial boundary theories is proposed to absorb stress waves in a saturated soil foundation in dynamic analysis. Since shear waves （S-waves） are the same in a saturated soil foundation and a single-phase medium foundation, a tangential visco-elastic boundary condition for a single-phase medium foundation can also be used for saturated soil foundations. Thus, the purpose of the artificial boundary proposed in this paper is primarily to absorb two types of P-waves in a saturated soil foundation. The main idea is that the stress of the P-waves in the saturated soil foundation is decomposed into two types. The first type of stress, δra＇ is absorbed by the first artificial boundary. The second type of stress, δrb, is balanced by the stress generated by the second artificial boundary. Ultimately, both types of P-waves （fast-P-waves and slow-P-waves） are absorbed by the artificial boundary model proposed in this paper. In particular, note that the fast-P-waves and slow-P-waves are absorbed at the position of the first boundary. Thus, the artificial boundary model proposed herein can simultaneously absorb P-fast waves, P-slow waves and shear waves. Finally, a numerical example is given to examine the proposed artificial boundary model, and the results show that it is very accurate.

Normal transmitting boundary conditions

The multi-transmitting formula (MTF) governed by a single artificial speed is analytically developed into a generalized MTF governed by a few artificial speeds to improve its capacity in simultaneous simulation of several one-way waves propagating at different speeds.The generalized MTF is then discretized and further generalized using the space extrapolation to improve its accuracies in numerical simulation of transient waves at large angles of incidence.The above two successive generalizitions of MTF based on the notion of normal transmission lead to a compact formula of local non-reflecting boundary condition.The formula not only provides a general representation of the major schemes of existing local boundary conditions but can be used to generate new schemes,which combine advantages of different schemes.

Finite Element - Artificial Transmitting Boundary Method for Acoustical Field on Tapered Waveguide

In earlier approach, the 2-D acoustical field profiles on the substrate region are often calculated with BPM. In this paper, we present a new approach based on the finite element -artificial transmitting boundary method and calculate acoustical field on the substrate region.

波动数值模拟中透射边界的稳定实现

从波动能量在计算区内累积增大的观点出发 ,通过简单的一维弹性波模型 ,系统地阐明了在近场波动数值模拟中透射边界两类数值失稳———“高频振荡”和“零频飘移”的机理 :前者源于对波动数值模拟无意义的高频波动在人工边界上的放大和波在有限计算区内多次反射产生的反复放大 ;后者则源于透射边界允许零频和接近零频的分量不断进入计算区 .由此提出了稳定实现透射边界的完整方案包括两项简单措施 :第一 ,在全部计算区内按文中建议的方法注入小阻尼 ,以消除高频振荡 ;第二 ,给出一种具有明确物理意义的消除零频飘移的算子算法 .最后 ,提供了三维波源问题和散射问题的详细数值试验结果 .

利用瑞利面波特性计算动力响应的进一步探讨

本文根据地基中瑞利面波的能流分布规律，提出了利用面波的第一模态波模拟水平向传播波，而竖向传播波则采用半无限单元模拟。通过算例验证了本文方法的可靠性和精确性，并与传统的传递边界法进行了比较，最后讨论了衰减系数中对计算结果的灵敏性。

基于黏弹性边界的拱坝地震反应分析方法

无限地基的能量辐射效应是影响拱坝地震反应的一个重要因素。本文提出将一种新的黏弹性人工边界结合显式有限元的时域波动求解方法用于拱坝—地基系统的地震反应分析,地震动输入通过分别将底边界入射位移时程和侧边界自由场时程对应的解析应力时程作为边界荷载实现。该方法具有较高的计算精度,且不存在高、低频计算失稳问题。文中以小湾拱坝地震反应分析为例,比较了透射边界、黏弹性边界两种计算方法的计算结果,论证了本文建议方法的有效性。

盾构隧道联络通道地震响应规律研究

联络通道与盾构隧道的连接处是隧道抗震的薄弱环节,但目前有关联络通道抗震性能的研究开展得还很少。对于隧道与联络通道连接部位这类典型的空间结构,三维动力有限元分析是最为有效的方法,而其中的人工边界条件和地震动输入方法是实现这一方法的关键技术之一。本文以显式形式将多次透射边界条件引入三维动力分析控制方程,采用分割算法分别求解边界节点和内部节点的响应,然后以对武汉长江隧道工程为背景,利用这种算法研究了隧道与联络通道连接处的地震响应规律,并对采用柔性搭接代替刚性接头和地基加固两种措施的减振效果进行了计算分析。

土-结构相互作用通用程序的改进

本文提出了用具有工程实用价值的透射边界替换通用程序FLUSH中－致边界的方法。通过实例证明，二阶透射边界与一致边界具有相当的精度。改进后的程序WMNS－Fl比FLUSH具有更高的计算效率和实用性。

水平成层场地地震反应非线性分析

本文首先推导了李小军积分格式(中心差分与Newmark平均加速度法相结合)的增量形式,并据此离散动力平衡方程,同时,采用Pyke提出的土动力本构模型以及多次透射人工边界条件,提出了一种水平成层场地地震反应非线性分析的显式有限元方法,并据此编制了计算机程序。数值实验表明,这种方法能较好地模拟土层在强地震作用下的非线性特性。

基于高阶双渐近透射边界的大坝-库水动力相互作用直接耦合分析模型

高阶双渐近时域透射边界能够同时模拟行波和快衰波的传播,并且能够在全频范围内迅速逼近准确解,具有优良的收敛性能和计算效率.本文将动水压力波高阶双渐近透射边界直接嵌入到近场有限元方程中,建立了大坝-库水动力相互作用的直接耦合分析模型.该模型的整体控制方程保留了近场有限元方程系数矩阵对称稀疏的优势,可以方便地利用现有的通用有限元求解器求解.基于有限元开源软件框架体系OpenSees(Open System for Earthquake Engineering Simulation),编程实现了直接耦合分析模型,并将其应用于二维重力坝、三维拱坝与库水动力相互作用分析.数值算例表明,该直接耦合分析模型具有很高的精度和计算效率.

三维土-结构动力相互作用的一种时域直接分析方法

本文提出了一种分析三维土-结构动力相互作用的时域直接方法。该方法采用集中质量显式有限元和透射人工边界模拟无限域地基,通过编制的FORTRAN程序实现;采用ANSYS软件对上部结构进行建模分析,并通过FORTRAN程序对ANSYS软件的调用,实现了土与结构系统在地震作用下的整体分析。该方法为显隐式相结合的方法,地基和上部结构可采用不同的时间步距进行分析,可大大提高效率。通过两算例,验证了该方法的可行性。

流体饱和两相多孔介质动力反应分析的显式有限元法

根据Biot流体饱和两相多孔介质波动方程 ,采用解耦技术 ,提出了分析流体饱和两相多孔介质中动力反应的显式有限元法。根据所推公式 ,编制了计算程序并进行了实例计算。计算结果与解析解对比 ,两者符合非常好 ,表明本文方法是处理流体饱和两相多孔介质动力问题的一种非常有效方法。