The effect of local irregular topography on seismic ground motion

Ｔｈｅｅｆｆｅｃｔｏｆｌｏｃａｌｉｒｒｅｇｕｌａｒｔｏｐｏｇｒａｐｈｙｏｎｓｅｉｓｍｉｃ　ｇｒｏｕｎｄ　ｍｏｔｉｏｎＪＩＧ－ｈｏＬＩＵ（刘晶波）（ＤｅｐａｒｔｍｅｎｔｏｆＣｉｖｉｌＥｎｇｉｎｅｅｒｉｎｇ，ＴｓｉｎｇｈｕａＵｎｉｖｅｒｓｉｔｙＢｅｉｊｉ...     

Analysis on effect of surface fault to site ground motion using finite element method

Dynamic contact theory is applied to simulate the sliding of surface fault.Finite element method is used to analyze the effect of surface fault to site ground motions.Calculated results indicate that amplification effect is obvious in the area near surface fault,especially on the site that is in the downside fault.The results show that the effect of surface fault should be considered when important structure is constructed in the site with surface fault.     

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.     

An explicit finite element method for dy—namic analysis in three—medium coupling system and its application

In this paper,an explicit finite element method to analyze the dynamic responses of three-medium coupled systems with any terrain is developed on the basis of the numerical simulation of the continuous conditions on the bounda-ries among fluid saturated porous medium,elastic single-phase medium and ideal fluid medium.This method is a very effective one with the characteristic of high calculating speed and small memory needed because the formulae for this explicit finite element method have the characteristic of decoupling,and which does not need to solve sys-tem of linear equations.The method is applied to analyze the dynamic response of a reservoir with considering the dynamic interactions among water,dam,sediment and basement rock.The vertical displacement at the top point of the dam is calculated and some conclusions are given.     

Research on TSIP method for medium－term earthquake prediction

ＲｅｓｅａｒｃｈｏｎＴＳＩＰｍｅｔｈｏｄｆｏｒｍｅｄｉｕｍ－ｔｅｒｍｅａｒｔｈｑｕａｋｅｐｒｅｄｉｃｔｉｏｎＤｅ－ＹｕＨＵＡＮＧ（黄德瑜）；Ｙｕａｎ－ＱｉｎｇＺＨＵ（朱元清）；ＹｏｎｇＣＨＥＮ（陈颙）ａｎｄＹｉｎｇＪＩ（季颖）（ＣｅｎｔｒｅｆｏｒＡｎ...     

Hybrid Laplace transform finite element method for solving the convection–dispersion problem

It can be very time consuming to use the conventional numerical methods, such as the finite element method, to solve convection–dispersion equations, especially for solutions of large-scale, long-term solute transport in porous media. In addition, the conventional methods are subject to artificial diffusion and oscillation when used to solve convection-dominant solute transport problems. In this paper, a hybrid method of Laplace transform and finite element method is developed to solve one- and two-dimensional convection–dispersion equations. The method is semi-analytical in time through Laplace transform. Then the transformed partial differential equations are solved numerically in the Laplace domain using the finite element method. Finally the nodal concentration values are obtained through a numerical inversion of the finite element solution, using a highly accurate inversion algorithm. The proposed method eliminates time steps in the computation and allows using relatively large grid sizes, which increases computation efficiency dramatically. Numerical results of several examples show that the hybrid method is of high efficiency and accuracy, and capable of eliminating numerical diffusion and oscillation effectively.     

Mechanical features of Shanxi Datong earthquake swarm sequence and their numerical siamulation by using finite element method

ＭｅｃｈａｎｉｃａｌｆｅａｔｕｒｅｓｏｆＳｈａｎｘｉＤａｔｏｎｇｅａｒｔｈｑｕａｋｅｓｗａｒｍｓｅｑｕｅｎｃｅａｎｄｔｈｅｉｒｎｕｍｅｒｉｃａｌｓｉａｍｕｌａｔｉｏｎｂｙｕｓｉｎｇｆｉｎｉｔｅｅｌｅｍｅｎｔｍｅｔｈｏｄＺｈｉ－ＬｉＺＨＡＮＧ，ＨｕａＴ...     

Development of attenuation relation for the near fault ground motion from the characteristic earthquake

A composite source model has been used to simulate a broadband strong ground motion with an associated fault rupture process. A scenario earthquake fault model has been used to generate 1 000 earthquake events with a magni-tude of Mw8.0. The simulated results show that, for the characteristic event with a strike-slip faulting, the character istics of near fault ground motion is strongly dependent on the rupture directivity. If the distance between the sites and fault was given, the ground motion in the forward direction (Site A) is much larger than that in the backward direction (Site C) and that close to the fault (Site B). The SH waves radiated from the fault, which corresponds to the fault-normal component plays a key role in the ground motion amplification. Corresponding to the sites A, B, and C, the statistical analysis shows that the ratio of their aPG is 2.15:1.5:1 and their standard deviations are about 0.12, 0.11, and 0.13, respectively. If these results are applied in the current probabilistic seismic hazard analysis (PSHA), then, for the lower annual frequency of exceedance of peak ground acceleration, the predicted aPG from the hazard curve could reduce by 30% or more compared with the current PSHA model used in the developing of seismic hazard map in the USA. Therefore, with a consideration of near fault ground motion caused by the rupture directivity, the regression model used in the development of the regional attenuation relation should be modified accordingly.     

Study on strength reduction factors considering the effect of classification of design earthquake

Introduction The strength reduction factor is defined as the lateral yielding strength required to avoid yielding in the system when subjected to a given ground motion, to the lateral yielding strength required to maintain the displacement ductility ratio demand equal to a pre-determined target duc-tility ratio under the same ground motion. The strength reduction factors are not only the key fac-tors in determining seismic action for force-based seismic design, but also one of the key parame-t…     

Inversion for negative dislocation on elastic block boundaries in the northeastern margin of Qinghai-Xizang block and prediction for strong earthquake location

Introduction According to the negative dislocation model (Matsu′ura et al, 1986), the relative motion be-tween active blocks under contemporary crustal movement is likely to be partially blocked on the boundaries. Suppose the lower ductile zone of boundary could slip freely, while due to the fric-tional resistance, etc., the upper brittle zone would restrict such kind of relative motion, so as to give rise to stress and strain accumulation. Namely, the surface displacement in the block bound-…     

黏弹介质波动方程有限差分解的稳定性研究

稳定性问题是地震波数值模拟的一个重要问题.基于地震波传播理论,从黏弹介质本构方程出发,对矩形网格下不同黏弹模型波动方程有限差分解的稳定性进行了理论分析,导出了Kelvin-Voigt黏弹模型和Maxwell黏弹模型在任意空间差分精度下稳定性条件的表达式;给出了品质因子Q≥5时的简化式,并通过数值算例验证了理论研究结论的正确性;总结了地震波速度、频率、空间网格大小、差分系数以及品质因子与稳定性条件的关系;通过误差分析给出了近似公式的使用条件.     

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

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

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.     

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…     

An optimized finite difference method based on a polar coordinate system for regional-scale irregular topography

The finite difference method (FDM) is an important numerical approach for simulating the propagation of seismic waves, and some FDMs can be used to study the impact of the Earth’s curvature and topography over large distances. To efficiently model the effects of the Earth’s irregular topography on the propagation of seismic waves, here we optimize a previously proposed grid mesh method and develop a novel two-dimensional boundary-conforming FDM based on a curvilinear polar coordinate system. This method efficiently simulates the propagation of seismic waves in an arc-shaped model with large variations in surface topography. Our method was benchmarked against other reported methods using several global-scale models. The consistency of the results confirms the validity of our proposed optimization strategy. Furthermore, our findings indicate that the proposed optimization strategy improves computational efficiency.     

显式积分格式对局部透射边界高频失稳的抑制和消除作用

针对显式有限元-有限差分法结合局部透射边界进行波动数值模拟的情况,推导并计算了循环系数,给出了计算过程的稳定条件,并通过数值试验对该判定准则进行验证。在给出的循环系数的基础上,定量地分析了积分格式的能耗特性对由边界条件引入的高频失稳的抑制和消除作用。探讨了如何更有效地利用能耗特性保证计算稳定性,即实际计算模型中时间离散步距以及空间模型大小等参数如何取值既能确保计算稳定又能最大程度减小计算量。     

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

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