An explicit method for numerical simulation of wave equations: 3D wave motion

In this paper, an explicit method is generalized from 1D and 2D models to a 3D model for numerical simulation of wave motion, and the corresponding recursion formulas are developed for 3D irregular grids. For uniform cubic grids, the approach used to establish stable formulas with 2M-order accuracy is discussed in detail, with M being a positive integer, and is illustrated by establishing second order (M=1) recursion formulas. The theoretical results presented in this paper are demonstrated through numerical testing.     

An explicit method for numerical simulation of wave equations

In this paper, a method to develop a hierarchy of explicit recursion formulas for numerical simulation in an irregular grid for scalar wave equations is presented and its accuracy is illustrated via 2-D and 1-D models. Approaches to develop the stable formulas which are of 2M-order accuracy in both time and space with M being a positive integer for regular grids are discussed and illustrated by constructing the second order (M = 1) and the fourth order (M = 2) recursion formulas.     

求解运动方程的一种不等时间步长的显式数值积分方法

在波动有限元模拟中, 若采用传统的显式数值积分方法求解运动方程, 计算时间步长需采用计算区内满足稳定条件要求的最小时间步长. 然而, 对于大部分计算区域, 这一时间步长过小, 是不必要的. 本文提出了一种不等时间步长的显式数值积分方法, 其基本思想是不同的计算区域采用满足各自稳定条件的计算时间步长. 最后, 本文通过数值试验检验了这一方法的可行性及其对数值计算精度的影响.      

大型结构地震反应值模拟中的波动输入

解耦的时域有限元数值模拟技术在考虑土－结构相互作用的大型结构地震反应分析中的应用日趋成熟。为完善其中的波动输入技术，本提出 了一种波动输入时步数值模拟的简便方法，替换对自由场的频域计算，从而使结构反应分析可以完全通过时步数值模拟实现。同时，通过数值实现，对这一实施方案的有效性进行了检验。     

近断层强地面运动影响场显式有限元数值模拟的示意性算例

本文通过一个简单的示意性算例，对近断层强地面运动影响场显式有限元数值模拟的主要环节作了概要说明，给出了算例的主要计算结果并对其作了定性和初步的讨论分析。主要观点和结论是：本文的计算结果较好地定性反映和印证了人们对近断层强地面运动特征（如上盘效应、竖向效应、Fling Step效应和速度大脉冲现象）的现有认识水平，验证了采用的计算模型和计算方法的可行性和科学性。本文关于位移峰值场的计算结果对于估算近断层地表破裂带的位置、长度和宽度具有重要参考价值，应进一步深入研究二者之间的相关性。     

Kinematic source model for simulation of near-fault ground motion field using explicit finite element method

1 Introduction The prediction of broadband near-fault ground motion is an important aspect of seismic risk analysis, and is essential for the evaluation of urban seismic safety. Many methodologies are used to model the near-fault ground motions, such as the Green function method, ?nite element (FE) / ?nite difference (FD) method and a hybrid method (Aagaard, 1999; Oprsal and Zahradnik, 2002; Oprsal et al., 2003 a and b). In the hybrid method, the FE / FD method is used to model the low fr…     

某正倾滑断层引起的近断层强地面运动的有限元数值模拟

本文使用与文献[1]相同的数值模拟方法,预测了某城市设定矩震级为MW=6.75的活断层产生的近断层强地面运动。主要定性计算结果是:正倾滑断层(倾角δ=75°)产生的近断层强地面运动也不同程度地存在上盘效应、F ling Step效应、速度大脉冲效应和竖向效应。表明本文采用的震源模型和计算方法可以实际应用于预测城市近断层强地面运动影响场。     

有限元算法在声波方程数值模拟中的频散分析

针对有限元算法在地震波数值模拟中的数值频散问题,利用集中质量矩阵双线性插值有限元算法,推导了二维声波方程的频散函数.在此基础上采用定量分析方法,对比分析了网格纵横长度比变化时的入射方向、空间采样间隔、地震波频率以及地层速度对数值频散的影响.数值算例和模型正演结果表明：当采用集中质量矩阵双线性插值有限元算法时,为了有效地压制数值频散,在所使用震源子波的峰值频率对应的波长内,采样点数目应不少于20个;减小网格长度的纵横比可以有效地抑制入射角(波传播方向与z轴的夹角)较小的地震波的数值频散;地震波频率越高,传播速度越慢,频散越严重,尤其是当相速度与其所对应的频率比值小于2倍空间采样间隔时,不仅会出现严重的数值频散,还会出现假频现象.     

Improved numerical dissipation for explicit methods in pseudodynamic tests

There is no second-order accurate, dissipative, explicit method in the currently available step-by-step integration algorithms. Two new families of second-order accurate, dissipative, explicit methods have been successfully developed for the direct integration of equations of motion in structural dynamics. These two families of methods are numerically equivalent and possess the desired numerical dissipation which can be continuously controlled. These two families of algorithms are very useful for pseudodynamic tests since the favourable numerical damping can be used to suppress the spurious growth of high-frequency modes due to the presence of numerical and/or experimental errors in performing a pseudodynamic test. © 1997 by John Wiley & Sons, Ltd.     

波动有限元方程显式逐步积分格式稳定性分析

对常用的波动有限元方程的两种显式格式的稳定性进行了分析,利用单元的最大频率乃是系统的特征频率的上界的概念结合模态分析方法,给出了便于实际应用的稳定性必要条件。同时,利用Von Neumann方法给出了这两种格式稳定性的充分条件,并通过算例对这两种稳定性条件进行验证。     

Development of a family of unconditionally stable explicit direct integration algorithms with controllable numerical energy dissipation

The implicit dissipative generalized‐ α method is analyzed using discrete control theory. Based on this analysis, a one‐parameter family of explicit direct integration algorithms with controllable numerical energy dissipation, referred to as the explicit KR‐α method, is developed for linear and nonlinear structural dynamic numerical analysis applications. Stability, numerical dispersion, and energy dissipation characteristics of the proposed algorithms are studied. It is shown that the algorithms are unconditionally stable for linear elastic and stiffness softening‐type nonlinear systems, where the latter indicates a reduction in post yield stiffness in the force–deformation response. The amount of numerical damping is controlled by a single parameter, which provides a measure of the numerical energy dissipation at higher frequencies. Thus, for a specific value of this parameter, the resulting algorithm is shown to produce no numerical energy dissipation. Furthermore, it is shown that the influence of the numerical damping on the lower mode response is negligible. It is further shown that the numerical dispersion and energy dissipation characteristics of the proposed explicit algorithms are the same as that of the implicit generalized‐ α method. A numerical example is presented to demonstrate the potential of the proposed algorithms in reducing participation of undesired higher modes by using numerical energy dissipation to damp out these modes. Copyright © 2014 John Wiley & Sons, Ltd.     

求解振动方程的一种显式积分格式及其精度与稳定性

介绍了一种求解有阻尼体系振动方程的显式积分格式及其逐步求解的过程，并对其计算精度和稳定性进行了分析。该方法不但能同时求得体系的位移、速度和加速度反应，而且所得到的加速度反应的精度能满足工程需要。     

A semi-analytical time integration for numerical solution of Boussinesq equation

A semi-analytical time integration method is proposed for the numerical simulation of transient groundwater flow in unconfined aquifers by the nonlinear Boussinesq equation. The method is based on the analytical solution of the system of ordinary differential equations with constant coefficients. While it is unconditionally stable and more accurate than the finite difference methods, the computational cost is much more expensive than (can be more than 10 times) that of the finite difference methods for a single time step. However, by partitioning the nonlinear parameters into linear and nonlinear parts, the costly computation can be performed only once. With larger and less variable time step sizes, the total computational cost can be significantly reduced. Three examples are included to illustrate the advantages and limitations of the proposed method.     

基于精细积分法的三维弹性波数值模拟(英文)

波动方程有限差分法是地震数值模拟中的一种重要的方法,对理解和分析地震传播规律、分析地震属性和解释地震资料有着非常重要的意义。但是有限差分法由于其离散化的思想,产生了不稳定性。精细积分法在有限差分法的基础上,在时间域采用解析解的表达形式,在空间域保留任意差分格式,发展成为半解析的数值方法。本文结合并发展了以往学者的成果,推导了任意精细积分法的三维弹性波正演模拟计算公式,并对其稳定性进行了数值分析。在计算实例中,实现了精细积分法二维和三维弹性波模型的地震正演模拟,对计算结果的分析表明,精细积分法反射信号走时准确,稳定性好,弹性波场相较于声波波场,弹性波波场成分更为丰富,包含了更多波型成分(PP-和PS-反射波、透射波和绕射波),这对实际地震资料的解释和储层分析有重要的意义。实践证明,该方法可直接应用到弹性波的地质模型的数值模拟中。     

Simulation of near-fault bedrock strong ground-motion field by explicit finite element method

Based on presumed active fault and corresponding model, this paper predicted the near-fault ground motion filed of a scenario earthquake （Mw=6 3/4 ） in an active fault by the explicit finite element method in combination with the source time function with improved transmitting artificial boundary and with high-frequency vibration contained. The results indicate that the improved artificial boundary is stable in numerical computation and the predicted strong ground motion has a consistent characteristic with the observed motion.     

采用显式有限元方法求解开放系统动力响应的主要环节和问题

分析了采用显式有限元方法求解开放系统动力响应的主要环节中影响精度和稳定性的因素及其相互关系,指出数值解的精度和稳定性是这些因素共同作用的结果。归纳了合理输入数字荷载应满足的条件以及系统离散化准则。基于有限元方法推导了吸收边界的一般表达式,并据此提出吸收边界和透射边界之间存在确定的函数关系,并不是两类独立人工边界的新见解。认为内域的阻尼稳定性和人工边界的近场稳定性是大型复杂开放系统动力响应数值方法研究和应用中的关键科学问题,并从应用的角度提出了解决问题的思路和看法。     

二维波动方程参数反演的微分进化算法

利用微分进化算法进行二维波动方程参数反演的研究.微分进化算法属于全局最优化算法,具有鲁棒性强、收敛速度快、计算精度高的优点.把参数反演问题转化为非线性优化问题,通过保持在解空间不同区域中各个点的搜索,以最大的概率找到问题的全局最优解.同时还进行了数值模拟,结果令人满意.     

Combined implicit or explicit integration steps for hybrid simulation

Explicit integration procedures have been widely adapted and applied to hybrid simulations of the seismic response of structures due to their simplicity. However, these procedures are only conditionally stable and have limited recent applications of hybrid simulations to simple structural models with few degrees of freedom. A novel integration procedure is proposed herein, in which a fully implicit formulation is applied to solve the equation of motion for the hybrid model, but defaults to an explicit or noniterative formulation in steps that fail to converge. The advantages to this approach are the ensured continuity of the simulation and the reduced accumulation of errors that occur during consecutive explicit steps that may lead to instability. The implicit procedure is applied by loading the experimental substructures beyond the expected displacement for the current step, then using the displacements and forces measured through the load path in the iterative implicit scheme. This approach captures the instantaneous behaviour of experimental substructures without physically imposing iterations. Numerical and experimental simulations demonstrate the effectiveness of the proposed integration scheme for multi‐degree‐of‐freedom models, especially in utilization of longer time steps that exceed stability limits of explicit methods, prevention of excitation of higher modes, and testing of stiff systems. Copyright © 2007 John Wiley & Sons, Ltd.     

瑞利阻尼介质有限元离散模型动力分析的数值稳定性

本文针对几种有一般阻尼的动力系数数值积分的显式方法，讨论了阻尼对稳定性的影响，并建议了瑞利阻尼介质有限元离散模型中动力分析数值稳定性的实用稳定判别方法。     

复杂地表边界元-体积元波动方程数值模拟

复杂近地表引起来自深部构造的地震反射信号振幅和相位的异常变化,是影响复杂近地表地区地震资料品质的主要原因.本文采用边界元-体积元方法,通过求解含复杂地表的波动积分方程,来模拟地震波在复杂近地表构造中的传播.其中,边界元法模拟地形起伏和表层地质结构对地震波传播的影响;体积元法模拟起伏地表下非均质低降速层的影响.与其他数值模拟方法比较,其主要优点为几何上精确描述不规则地表界面,实现精确模拟自由表面对地震波的边界散射;显式应用近地表地层界面的连续边界条件,实现半解析的数值模拟;分区处理近地表复杂结构,有效模拟复杂地表下非均匀介质对地震波场的体散射.数值试验结果表明了该方法的实用性和有效性.