Efficient time truncation in two-dimensional bem analysis of transient wave propagation problems

An efficient time truncation algorithm for transient two-dimensional wave propagation analysis by the boundary element method is proposed. First the complete (without any truncation) formulation is reviewed and subsequently time truncation concepts are discussed. Three applications are presented, where the efficiency of complete and truncated algorithms is compared. The most important conclusions inferred from the numerical analyses are: (i) time truncation errors can be controlled and made as small as required and (ii) time truncation algorithms lead to considerable savings in both CPU time and storage area.     

Assessment of errors in the newmark method in structural dynamics

In the Newmark and other approximate step-by-step methods, having introduced assumptions in order to transform the differential equations, which are characteristic of response problems, into simultanéous equations, successive solutions lead to a response-time history. In this paper numerical results and formulae are given for the errors which are generated by this procedure. These errors are oscillatory in nature and, in general, the oscillations increase in magnitude as the number of time steps increases. Recommendations for upper limits on the time step, which will provide acceptable accuracy for a wide range of system and excitation parameters, are presented.     

局部解耦的时域波分析方法

从模拟散射波的观点出发,推导了波动问题局部解耦的人工边界条件,与有限元内点的时域显式方法相结合建立了波动分析的局部解耦的直接分析方法,文中还从波动分解的概念,建立了外源问题（如地震作用）输入方法。     

Macro–micro analysis method for wave propagation in stochastic media

This paper presents a new analysis method, called macro–micro analysis method (MMAM) for numerical simulation of wave propagation in stochastic media, which could be used to predict distribution of earthquake strong motion with high accuracy and spatial resolution. This MMAM takes advantage of the bounding medium theory (BMT) and the singular perturbation expansion (SPE). BMT can resolve uncertainty of soil and crust structures by obtaining optimistic and pessimistic estimates of expected strong motion distribution. SPE leads to efficient multi‐scale analysis for reducing a huge amount of computation. The MMAM solution is given as the sum of waves of low resolution covering a whole city and waves of high resolution for each city portion. This paper presents BMT and SPE along with the formulation of MMAM for wave propagation in three‐dimensional elastic media. Application examples are presented to verify the validity of the MMAM and demonstrate potential usefulness of this approach. In a companion paper (Earthquake Engng. Struct. Dyn., this issue) application examples of earthquake strong motion prediction are also presented. Copyright © 2005 John Wiley & Sons, Ltd.     

瞬态面波法在滑坡勘察中的应用效果

通过对瞬态面波法的工作原理、野外工作方法、资料整理和成果解释的描述和对应用实例的比较与分析,认为瞬态面波法在滑坡勘察中确定滑坡体的岩土分层、软弱带、滑动面是可行和可靠的。     

基于时频分析方法提取瞬态瑞利波频散曲线

作为一种无损浅层勘探方法,瑞利波越来越多地应用于地震工程、环境工程及岩土工程测试中。其中,一项重要的工作是提取瑞利波的频散曲线。根据激振方式不同,瑞利波勘探可分为稳态法、瞬态法和地脉动台阵法。基于相位谱法提取瞬态瑞利波的频散曲线,一般将两通道瑞利波的相位差限制在2π之内,因此限制了两拾振器之间的布点距离。本文利用时频分析方法,分析了瞬态瑞利波法频散曲线的提取方法,避免了对相位差的限制。数值试验结果表明,采用时频分析方法提取到的频散曲线具有较高的精度。     

局部角度域波传播步进算法研究

本文从非均匀介质中波动方程出发,提出了基于一般标架的相空间（局部角度域）波传播的步进算法. 该方法在构造单程波的步进算法时,在选择标架或正交基等方面有更大的自由度. 我们以不随频率及深度变化且具有变尺度特性的Gabor_Daubechies紧标架为例,给出了单程波传输算子的具体形式及相应的波场步进算法;详细讨论了基于Gabor_Daubechies标架的传输算子的高频渐近展开问题,得出了在高频、小传输步长条件下传输算子的近似解析表达式,并给出使用条件. 通过模型算例,比较了精确传输算子与高频近似传输算子用于非均匀介质中波传播的结果,说明在一定条件下由两者得出的波场几乎是相同的.     

Condensed hyperelements method of non-vertical consistent boundaries for wave propagation analysis in irregular media

The study of wave propagation in finite/infinite media has many applications in geotechnical and structural earthquake engineering and has been a focus of research for the past few decades. This paper presents an analysis of 2D anti- plane problems （Love waves） and 2D in-plane problems （Rayleigh waves） in the frequency domain in media consisting of a near-field irregular and a far-field regular part. The near field part may contain structures and its boundaries with the far-field can be of any shape. In this study, the irregular boundaries of the near-field are treated as consistent boundaries, extending the concept of Lysmer＇s vertical consistent boundaries. The presented technique is called the Condensed Hyperelements Method （CHM）. In this method, the irregular boundary is limited to a vertical boundary at each end that is a consistent boundary at the far-field side. Between the two ends, the medium is discretized with hyperelements. Using static condensation, the stiffness matrix of the far-field is derived for the nodes on the irregular boundary. Examples of the application of the CHM illustrate its excellent accuracy and efficiency.     

A post-processing technique and an a posteriori error estimate for the newmark method in dynamic analysis

In this paper, we present a post-processing technique and an a posteriori error estimate for the Newmark method in structural dynamic analysis. By post-processing the Newmark solutions, we derive a simple formulation for linearly varied third-order derivatives. By comparing the Newmark solutions with the exact solutions expanded in the Taylor series, we achieve the local post-processed solutions which are of fifth-order accuracy for displacements and fourth-order accuracy for velocities in one step. Based on the post-processing technique, a posteriori local error estimates for displacements, velocities and, thus, also the total energy norm error estimate are obtained. If the Newmark solutions are corrected at each step, the post-processed solutions are of third-order accuracy in the global sense, i.e. one-order improvement for the original Newmark solutions is achieved. We also discuss a method for estimating the global time integration error. We find that, when the total energy norm is used, the sum of the local error estimates will give a reasonable estimate for the global error. We present numerical studies on a SDOF and a 2-DOF example in order to demonstrate the performance of the proposed technique.     

Time-domain analysis of wave propagation in 3-D unbounded domains by the scaled boundary finite element method

Transient wave propagation in three-dimensional unbounded domains is studied. An efficient numerical approach is proposed, which is based on using the displacement unit-impulse response matrix representing the interaction force–displacement relationship on the near field/far field interface. Spatially, an approximation is used to reduce the computational effort associated with the large size of three-dimensional problems. It is based on subdividing the fully coupled unbounded domain into multiple subdomains. The displacement unit-impulse response matrices of all subdomains are calculated separately. The error associated with this spatial decoupling can be reduced by placing the near field/far field interface further away from the domain of interest. Detailed parameter studies have been conducted using numerical examples, in order to provide guidelines for the proposed spatially local schemes, and to demonstrate the accuracy and high efficiency of the proposed method for three-dimensional soil–structure interaction problems.     

饱水介质波动分析的谐波输入初始间断的处理方法

本文以饱水两相介质的土力学模型为研究对象,在假定两相介质为弹性介质条件下,采用了显式有限元法和透射边界进行了饱和弹性半空间动力响应问题的研究。为避免谐波输入初始间断的影响,文中提出了一个处理函数,并以弹性半空间为算例,对饱水介质和单相介质分别进行了在底边界P波垂直入射时的动力响应分析,验证了该处理函数的有效性和实用性。     

Normally incident wave propagation across a joint set with the virtual wave source method

This paper presents an analytical study on normally incident wave propagation across a joint set. Combined with the displacement discontinuity method (DDM), the virtual wave source method (VWSM) is adopted to take account of multiple wave reflections among joints. After theoretical formulation, a comparison between results from VWSM and those from the method of characteristics (MC) is carried out. These results agree well with each other, and the applicability of VWSM into wave propagation study is verified. Subsequently, an interesting phenomenon, i.e., the magnitude of transmission coefficient increases with joint number on some circumstances, is pointed out and studied. By analyzing differently arriving transmitted waves and their superposition, such phenomenon is explained.     

Axisymmetric wave propagation in multilayered poroelastic grounds due to a transient acoustic point source

This paper deals with the study of axisymmetric wave propagation in various acoustic/porous stratified media coupling configurations. It presents the theoretical developments of a semi-analytical method, its validation for a limit test-case half-space ground, and an extension to a realistic multilayered seabed, when spherical waves are emitted from a transient point source in water.     

A simplified input output relation method using AR model for earthquake wave propagation analysis

This paper introduces a Simplified Input Output Relation Method (SIORM) using multiple-variable autoregression (AR) model which can be used to determine ground wave propagation properties. Using the AR model, a method is developed to establish the basis for the formulation of SIORM. The degree of accuracy of this method is evaluated. SIORM is applied to actual ground acceleration records taken at Minamisuna and Etchujima sites in Japan and its results are discussed in this paper. © 1997 John Wiley & Sons, Ltd.     

Obliquely incident wave propagation across rock joints with virtual wave source method

Due to the presence of joints, waves are greatly attenuated when propagating across rock masses. Zhu et al. (2011) (Normally incident wave propagation across a joint set with virtual wave source method. J. Appl. Geophys.73, 283–288.) studied normally incident wave propagation across a joint set with the virtual wave source method (VWSM). The introduced VWSM has merits in some aspects, especially the capability of separating differently arriving transmitted waves. However, normal wave incidence is only the special case for wave incidence with arbitrary incident angles. Obliquely incident wave propagation across a joint set is more complicated than normally incident wave propagation due to wave transformation at the joints. As a continuation of the previous paper, this work is extended to analytically study obliquely incident wave propagation across joints with VWSM. Complete theoretical reflection and transmission coefficients across single joint described by displacement discontinuity model are derived through plane wave analysis. The superposition of P wave and S wave is for the first time mathematically expressed and studied. The VWSM is verified through comparison with the propagation matrix method. Through extensive parametric studies on wave transmission across single and multiple parallel joints, it is shown that transmitted wave energy is mainly constrained in the transmitted wave of the same type as the incident wave. And with increasing joint stiffness, the transmission coefficients across single joint increases except those whose wave type is different from the incident wave. The amplitude of superposed transmitted wave for P wave incidence increases with incident angle, which is coincident with field observations. Both joint spacing and number of joints have significant effects on transmission coefficients. We find that when joint spacing is sufficiently large, the transmission coefficient is no longer a constant as the normally incident wave propagation case (Zhu et al., 2011). And when joints are very closely spaced, wave attenuation depends little on the number of joints, which is different from the conclusions from equivalent medium method.     

Pure acoustic wave propagation in transversely isotropic media by the pseudospectral method

Seismic wave propagation in transversely isotropic (TI) media is commonly described by a set of coupled partial differential equations, derived from the acoustic approximation. These equations produce pure P‐wave responses in elliptically anisotropic media but generate undesired shear‐wave components for more general TI anisotropy. Furthermore, these equations suffer from instabilities when the anisotropy parameter ε is less than δ. One solution to both problems is to use pure acoustic anisotropic wave equations, which can produce pure P‐waves without any shear‐wave contaminations in both elliptical and anelliptical TI media. In this paper, we propose a new pure acoustic transversely isotropic wave equation, which can be conveniently solved using the pseudospectral method. Like most other pure acoustic anisotropic wave equations, our equation involves complicated pseudo‐differential operators in space which are difficult to handle using the finite difference method. The advantage of our equation is that all of its model parameters are separable from the spatial differential and pseudo‐differential operators; therefore, the pseudospectral method can be directly applied. We use phase velocity analysis to show that our equation, expressed in a summation form, can be properly truncated to achieve the desired accuracy according to anisotropy strength. This flexibility allows us to save computational time by choosing the right number of summation terms for a given model. We use numerical examples to demonstrate that this new pure acoustic wave equation can produce highly accurate results, completely free from shear‐wave artefacts. This equation can be straightforwardly generalized to tilted TI media.     

岩石样品中超声波传播数值模拟及分析

为了深入研究超声波速度测试中超声波在岩石样品中的传播规律,对岩石超声波速度实验采用交错网格有限差分算法进行了数值模拟.选取了适合于超声波震源的差分精度、稳定性条件和吸收边界.通过模拟波形与实验记录波形的对比分析,结果表明:激发器和接收器直接位于岩心与探头接触面时,接收到的波形中干扰波较少,合成记录波形特征明显,易于分析;使用雷克子波作为超声波数值模拟的震源子波,计算得到实际首波初至时间为21.55μs,而波形图中首波初至时间为23.90μs,这是由于子波自身的长度使得波初至时间比实际波初至的时间晚;相对于岩心与探头接触面的反射波,岩心柱体与空气接触面的反射波对有效波形信号的影响更大.