Dynamic analysis of large 3-D underground structures by the bem

An up to date literature survey on the dynamics of underground structures is presented briefly. The dynamic response of large three-dimensional underground structures to external or internal dynamic forces or to seismic waves is numerically determined by the frequency domain boundary element method. This method is used to model both the structure and the soil medium, which are assumed to behave as linear elastic or viscoelastic bodies. The full-space dynamic fundamental solution is employed in the formulation and this requires a free soil surface discretization, confined to a finite portion around the area of interest, in addition to soil—structure interface and free structural surface discretizations. The dynamic disturbances can have a harmonic or a transient time variation. The transient case is treated with the aid of numerical Laplace transforms with respect to time. Various numerical examples involving lined cavities and long lined tunnels buried in the full- or the half-space subjected to harmonic or transient external forces or seismic waves are presented to illustrate the method and demonstrate its advantages.     

3-D dynamic foundation-soil-foundation interaction on layered soil

In this work the interaction between adjacent rigid, surface foundations resting on a viscoelastic layered soil medium is studied. A 3-D frequency domain BEM formulation in conjunction with infinite space fundamental solutions and the so called `successive stiffness method', initially developed for elastostatics and adapted here for the solution of elastodynamic problems, are used for the simulation of a layered soil medium. As a result, a discretization of the soil-foundation interface and the surrounding free surface as well as the soil layers' interfaces is necessary. However, it is shown in this work that reasonably accurate results can be obtained by using a substantially reduced discretization scheme involving only a small portion of the free surface surrounding the foundation and the corresponding interfaces of the soil layers. The presented numerical results demonstrate the importance of the dynamic foundation-soil-foundation interaction phenomenon which becomes even more pronounced where the supporting soil medium is made up of relatively shallow layers close to its free surface.     

Direct BEM for high-resolution global gravity field modelling

The paper presents a high-resolution global gravity field modelling by the boundary element method (BEM). A direct BEM formulation for the Laplace equation is applied to get a numerical solution of the linearized fixed gravimetric boundary-value problem. The numerical scheme uses the collocation method with linear basis functions. It involves a discretization of the complicated Earth’s surface, which is considered as a fixed boundary. Here 3D positions of collocation points are simulated from the DNSC08 mean sea surface at oceans and from the SRTM30PLUS_V5.0 global topography model added to EGM96 on lands. High-performance computations together with an elimination of the far zones’ interactions allow a very refined integration over the all Earth’s surface with a resolution up to 0.1 deg. Inaccuracy of the approximate coarse solutions used for the elimination of the far zones’ interactions leads to a long-wavelength error surface included in the obtained numerical solution. This paper introduces an iterative procedure how to reduce such long-wavelength error surface. Surface gravity disturbances as oblique derivative boundary conditions are generated from the EGM2008 geopotential model. Numerical experiments demonstrate how the iterative procedure tends to the final numerical solutions that are converging to EGM2008. Finally the input surface gravity disturbances at oceans are replaced by real data obtained from the DNSC08 altimetryderived gravity data. The ITG-GRACE03S satellite geopotential model up to degree 180 is used to eliminate far zones’ interactions. The final high-resolution global gravity field model with the resolution 0.1 deg is compared with EGM2008.     

Dynamic interaction numerical models in the time domain based on the high performance scaled boundary finite element method

Consideration of structure-foundation-soil dynamic interaction is a basic requirement in the evaluation of the seismic safety of nuclear power facilities. An efficient and accurate dynamic interaction numerical model in the time domain has become an important topic of current research. In this study, the scaled boundary finite element method （SBFEM） is improved for use as an effective numerical approach with good application prospects. This method has several advantages, including dimensionality reduction, accuracy of the radial analytical solution, and unlike other boundary element methods, it does not require a fundamental solution. This study focuses on establishing a high performance scaled boundary finite element interaction analysis model in the time domain based on the acceleration unit-impulse response matrix, in which several new solution techniques, such as a dimensionless method to solve the interaction force, are applied to improve the numerical stability of the actual soil parameters and reduce the amount of calculation. Finally, the feasibility of the time domain methods are illustrated by the response of the nuclear power structure and the accuracy of the algorithms are dynamically verified by comparison with the refinement of a large-scale viscoelastic soil model.     

Coordinate mapping of analytical contaminant transport solutions to non-uniform flow fields

Existing analytical solutions to 2D and 3D contaminant transport problems are limited by the mathematically convenient assumption of uniform flow. An approximate method is developed herein for coordinate mapping of 2D (vertically-averaged) transport solutions to non-uniform steady-state irrotational and divergence-free flow fields in single-layer aquifers. The method enables existing analytical transport solutions to be applied to aquifer systems with wells, non-uniform saturated thickness, surface water features, and (to a limited degree) heterogeneous hydraulic conductivity and recharge. This mass-conservative coordinate mapping approach is inexact in its approximation of the dispersion process but is still sufficiently accurate for many simple flow systems. The degree of model error is directly proportional to the variation of velocity magnitude within the domain. These mapped analytical solutions are compared to numerical simulation results and the coordinate mapping errors are investigated. The methods described herein may be used in the traditional capacity of analytical transport models, i.e., screening and preliminary site assessment, without sacrificing accuracy by assuming locally uniform flow conditions or applying an ad-hoc coordinate transformation. The solutions benefit from the traditional advantages of analytical methods, particularly the removal of artifacts due to spatial and temporal discretization: no time-stepping or numerical discretization is required.     

三维时间域航空电磁有理Krylov正演研究

常规的三维时间域航空电磁模拟通常采用隐式步长方法进行时间离散,需要几次矩阵分解和上百次右端源项回带,计算效率较低.为了提高正演计算效率,本文提出使用有理Krylov方法求解时间域电场扩散方程.首先使用非结构四面体网格进行空间离散,采用Nédélec矢量基函数近似四面体单元内的电场;然后基于有限元离散给出矩阵指数和矢量乘积表示的电场显式解;最后采用有理Arnoldi算法构造Krylov子空间内的正交基函数并进一步求解矩阵指数与矢量的乘积,直接得到任意时刻的电场解向量,避免步长离散过程.此外,本文还提出一种指数加权偏移参数优化方法,使得有理Arnoldi近似在瞬变衰减晚期具备更高的精度,从而降低Krylov子空间阶数并提高计算效率.通过和层状模型解析解的对比验证了有理Krylov方法的精度.针对三维异常体模型使用全局网格和局部网格剖分并和其他数值方法比较,进一步说明了有理Krylov方法的有效性.     

Non-singular time domain BEM with applications to 3D inertial soil–structure interaction

A direct time domain boundary element method is presented based on the Stokes fundamental solutions, discretized in both time and space, and an efficient time step-by-step solution that minimizes the accumulation of errors. A non-singular numerical integration procedure, in the Cauchy sense, is proposed for the generation of the associated influence matrices. This methodology is shown to be efficient for the solution of a number of computationally intensive problems in the area of soil–structure interaction. In addition, an algorithm for the direct calculation of the response of massive foundations to externally applied forces and/or obliquely incident seismic waves is introduced. The accuracy and computational efficiency of the proposed methodologies is established through a number of comparison studies.     

Dynamic soil–structure interaction analysis via coupled finite-element–boundary-element method

In this paper, a study on the transient response of an elastic structure embedded in a homogeneous, isotropic and linearly elastic half-plane is presented. Transient dynamic and seismic forces are considered in the analysis. The numerical method employed is the coupled Finite-Element–Boundary-Element technique (FE–BE). The finite element method (FEM) is used for discretization of the near field and the boundary element method (BEM) is employed to model the semi-infinite far field. These two methods are coupled through equilibrium and compatibility conditions at the soil–structure interface. Effects of non-zero initial conditions due to the pre-dynamic loads and/or self-weight of the structure are included in the transient boundary element formulation. Hence, it is possible to analyse practical cases (such as dam–foundation systems) involving initial conditions due to the pre-seismic loads such as water pressure and self-weight of the dam. As an application of the proposed formulation, a gravity dam has been analysed and the results for different foundation stiffness are presented. The results of the analysis indicate the importance of including the foundation stiffness and thus the dam–foundation interaction.     

Dynamic response of a multi-layered poroelastic medium

An exact stiffness matrix method is presented to evaluate the dynamic response of a multi-layered poroelastic medium due to time-harmonic loads and fluid sources applied in the interior of the layered medium. The system under consideration consists of N layers of different properties and thickness overlying a homogeneous half-plane or a rigid base. Fourier integral transform is used with respect to the x-co-ordinate and the formulation is presented in the frequency domain. Fourier transforms of average displacements of the solid matrix and pore pressure at layer interfaces are considered as the basic unknowns. Exact stiffness (impedance) matrices describing the relationship between generalized displacement and force vectors of a layer of finite thickness and a half-plane are derived explicitly in the Fourier-frequency space by using rigorous analytical solutions for Biot's elastodynamic theory for porous media. The global stiffness matrix and the force vector of a layered system is assembled by considering the continuity of tractions and fluid flow at layer interfaces. The numerical solution of the global equation system for discrete values of Fourier transform parameter together with the application of numerical quadrature to evaluate inverse Fourier transform integrals yield the solutions for poroelastic fields. Numerical results for displacements and stresses of a few layered systems and vertical impedance of a rigid strip bonded to layered poroelastic media are presented. The advantages of the present method when compared to existing approximate stiffness methods and other methods based on the determination of layer arbitrary coefficients are discussed.     

On rigid foundations subjected to seismic waves

In the analysis of structural foundations for seismic loads, it is customary to distinguish two types of soil-structure interaction effect: kinematic interaction (or wave passage), and inertial interaction. The former refers to the phenomenon of wave scattering, which occurs because the foundation is much stiffer than the surrounding soil and cannot accommodate to its distortions. Inertial interaction, on the other hand, is caused by feedback of kinetic energy of the structure into the soil. This paper is concerned only with the first phenomenon. The rigorous analysis of rigid, embedded foundations subjected to seismic disturbances requires, in general, substantial computational effort. Indeed, a typical analysis would normally require models with finite elements and/or boundary elements. Although such methods may be used to find an accurate solution to the problem of kinematic interaction, their use is not always warranted, given the many uncertainties involved and the multitude of assumptions that must be considered. Hence, approximate solutions are attractive for this problem. One such approximate method is the remarkably simple algorithm proposed by Iguchi.³ This paper presents first an appraisal of this method by way of a comparison with accurate numerical solutions for cylindrical foundations; next the algorithm is applied to rectangular (prismatic) foundations. It is found that Iguchi's method gives results that are adequate for engineering purposes, even if not entirely accurate.     

Analytical solutions to nonlinear conservative oscillator with fifth-order nonlinearity

This paper describes analytical and numerical methods to analyze the steady state periodic response of an oscillator with symmetric elastic and inertia nonlinearity. A new implementation of the homotopy perturbation method (HPM) and an ancient Chinese method called the max-min approach are presented to obtain an approximate solution. The major concern is to assess the accuracy of these approximate methods in predicting the system response within a certain range of system parameters by examining their ability to establish an actual (numerical) solution. Therefore, the analytical results are compared with the numerical results to illustrate the effectiveness and convenience of the proposed methods.     

求解弹性波方程的辛RKN格式

将弹性波方程变换至Hamilton体系,构造适用于弹性波模拟的高效显式二阶辛Runge-Kutta-Nystrm(RKN)格式,运用根数理论得到此格式的阶条件方程组.通过给定系数的限定条件,得到方程的对称解.为了使时间离散误差达到极小,提出数值频率与真实频率比较,通过Taylor展开,得到关于辛系数的限定方程,求解方程组得到最小频散辛RKN格式.对比分析时间演进方程的稳定性,得到使库朗数达到极大值的限定方程,求解方程组得到最稳定辛RKN格式.发现此两种格式为同一格式.新得到的辛RKN格式不依赖于空间离散方法,为了对比的需要,选取有限差分法进行空间离散.在频散、稳定性分析中,与常见辛格式对比,从理论上分析了本文提出的格式在数值频散压制、稳定性提升等方面的优势,数值实验进一步证实了理论分析的正确性.     

One-dimensional modelling of the non-linear far field in soil–structure-interaction analysis

In a bounded domain elasto-plastic wave propagation can be modelled accurately using the finite-element method. As is even the case for an elastic analysis, an unbounded domain, e.g. a semi-infinite soil or fluid, can, however, not be represented in this manner, as any spatial discretization has to be avoided. For one-dimensional wave propagation with a bi-linear elasto-plastic material law involving one stress component an analytical solution exists. The latter is used in modelling the non-linear far field of an unbounded medium using a rigorous bookkeeping procedure of the generated elastic and plastic waves propagating in both directions. The need for a non-linear model of the far field arises, as in a two-dimensional representation of soil-structure interaction the surface waves do not decay.     

3-D seismic response analysis of long lined tunnels in half-space

The dynamic response of infinitely long lined tunnels with a uniform cross-section buried into an elastic or viscoelastic half-space to body and surface harmonic seismic waves is numerically determined by a special direct boundary element method in the frequency domain. The waves have an arbitrary direction of propagation with respect to the axis of the tunnel and this renders the problem three-dimensional. However, this problem is effectively reduced to a two-dimensional one by a coordinate transformation and appropriate integration of the full-space dynamic fundamental solution along the direction of the tunnel axis. Quadratic isoparametric boundary line elements and advanced numerical integration techniques for the treatment of singular integrals produce results of high accuracy. Numerical results are presented for the case of an infinitely long lined tunnel of circular cross-section and compared against those of a full three-dimensional boundary element analysis, as well as those of other methods. Thus the proposed method is illustrated and its merits demonstrated.     

Analytical and numerical analysis of pressure drawdown in a poroelastic reservoir with complete overburden effect considered

Starting from an analytical reservoir model that incorporates full interaction with an elastic overburden, a new hybrid mathematical approach is developed by combining two numerical discretization methods. A tabular reservoir (petroleum reservoir or an aquifer) in an infinite or semi-infinite domain is viewed as a macroscopic displacement discontinuity, allowing use of the efficient displacement discontinuity mathematical method to calculate stresses and displacements that arise because of pressure changes. A 3-D finite element method using a poroelastic formulation is used to discretize the reservoir itself. By coupling the displacement discontinuity and finite element methods, a 3-D large-scale poroelastic reservoir can be simulated within an infinite or semi-infinite domain. The numerical model has been verified through comparison to known solutions, and some time-dependent pressure drawdown problems are analyzed. Results indicate that including the complete overburden (reservoir surroundings) response has a significant effect on pressure drawdown in a poroelastic reservoir during pumping, and should be incorporated in appropriate applications such as well test equations and subsidence analyses.     

基于比例边界有限元方法的拱坝子结构分区形式研究

比例边界有限元方法是一种半解析的数值计算方法,具有降维、网格灵活、严格模拟无限域和无需基本解等特点。比例边界有限元方法的基本理论是在整体坐标与局部坐标的比例边界转换基础之上建立的,相似中心的选取是否合理对分析计算具有重要的影响,导致在模拟拱坝这种不规则的空间壳体结构时,具有一定的局限性。采用子结构方法,将坝体分为若干满足相似性要求的区域可解决上述问题,以某拱坝为例给出了合理的坝体子结构分区形式,验证了子结构方法的精确性,为建立基于比例边界有限元方法的坝体-库水-地基系统的计算模型奠定了基础。     

地震波动方程的局部间断有限元方法数值模拟

近年来,随着地震波数值模拟对计算精度和效率的要求越来越高,间断有限元方法开始受到越来越多的关注.本文中,针对具有吸收边界条件的二维地震声波波动方程,作者提出了一种基于局部间断有限元方法的数值模拟算法.该算法在空间上使用局部间断有限元方法进行离散,在时间上采用了显式蛙跳格式.在这种时空离散的组合方式下,每个时间步上,此算法在空间剖分的每个单元上的求解计算是相互独立的,因而具有极高的并行性.通过数值算例,我们将该算法与连续有限元方法进行了比较.结果表明,本算法不仅具有对起伏构造的良好适应性,而且在计算效率和计算精度等方面,都具有优越性.     

多源条件下直流电阻率法有限元三维数值模拟中一种近似边界条件

在多源直流电阻率法有限元三维数值模拟中,传统混合边界条件由于刚度矩阵与源点位置相关,无法实现线性方程组的快速回代求解,目前常用齐次边界条件或无限元边界进行替代,虽然实现了快速回代求解,但同时也降低了数值模拟的精度.为了实现快速回代求解,并确保数值模拟的计算精度,本文提出了一种近似边界条件方法.其核心思想是将与场源位置相关的边界系数矩阵从刚度矩阵中分离出来,使得分离后的刚度矩阵与场源位置无关.并将边界系数矩阵与边界处一次电场向量的乘积移到线性方程组右端源项中,当场源位置发生改变时,只需要重新计算右端源项就可以实现快速回代求解.理论模型数值计算表明,在水平地形条件下,本文边界条件数值精度优于混合边界条件;在起伏地形条件下,与齐次边界条件相比,本文边界条件数值结果与混合边界条件吻合度更高.     

Time-domain two-dimensional site response analysis of non-homogeneous topographic structures by a hybrid BE/FE method

In this paper, an advanced formulation of the time-domain, two-dimensional hybrid finite element–boundary element method (FEM/BEM) is presented, and applied to carry out site response analysis of homogeneous and non-homogeneous topographic structures subjected to incident in-plane motions. Seismic responses of half-plane, horizontally layered site, alluvial valley and ridge sections subjected to incident P and SV waves are analyzed in order to demonstrate the applicability and efficiency of the presented method. The numerical results show that hybrid BE/FE methods require smaller time steps than those needed by BEM schemes. They also show that in case of surface irregularities with height to half-width ratio of up to one, the topography effect could be noticeable, if incident waves have wavelengths of less than approximately eight times the width.