A transmitting boundary for the dynamic finite element analysis of cross anisotropic soils

The problem of the transmitting boundary used in the dynamic finite element analysis of layered axisymmetric soil models with non-axisymmetric displacements is considered. Each layer is modelled as a homogeneous, viscoelastic cross anisotropic medium with a vertical axis of material symmetry.     

A high-order time-domain transmitting boundary for cylindrical wave propagation problems in unbounded saturated poroelastic media

Based on the u–p formulation of Biot equation with an assumption of zero permeability coefficient, a high-order transmitting boundary is derived for cylindrical elastic wave propagation in infinite saturated porous media. By this transmitting boundary the total stresses on the truncated boundaries of a numerical model, such as a finite element model, are replaced by a set of spring, dashpot and mass elements, with some additionally introduced auxiliary degrees of freedom. The transmitting boundaries are incorporated into the DIANA SWANDYNE II program and an unconditionally stable implicit time integration algorithm is adopted. Despite the assumption made in the derivation of the transmitting boundary, numerical examples show that it can provide highly accurate results for cylindrical elastic wave propagation problems in infinite saturated porous medium in case the u–p formulation is applicable. Although the direct applications of the proposed transmitting boundary to general two dimensional wave problems in infinite saturated porous media are not highly accurate, acceptable accuracy can still be achieved by placing the transmitting boundary at relatively large distance from the wave source.     

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

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

Numerical instabilities of a local transmitting boundary

A local transmitting boundary is presented in a compact form, which can be directly incorporated into finite elements. The accuracy of the boundary is studied thoroughly for a one-dimensional model in order to clarify numerical instabilities introduced by the boundary. Discretization of the model and reflection from the boundary are rigorously considered in the study, and the mechanism of the instability is then illuminated in the frequency domain by the amplification of reflection from the boundary and the multi-reflection of wave motion in a finite computational region. Typical characteristics of the instability in the time domain are illustrated by numerical results of the simple model and explained completely by the mechanism. On the basis of this understanding of the mechanism, a modified transmitting boundary is devised and its stability criterion is given for the one-dimensional model.     

一种基于多项式外推的局部透射边界位移解(外行波为非平面波情形)

一般情况下,通过人工边界向外透射的波动常为非平面波。对于近场非平面波在人工边界上的透射问题,不能采用简单的平面波透射方法。本文针对近场非平面波在人工边界上的透射特点,分析了非平面波沿人工边界法向视波速的变化规律及其近场失稳机制,推导了适合近场非平面波的、人工边界点位移解的不等步距线性外推公式。     

一种基于多项式外推的局部透射边界位移解(外行波为平面波情形)

基于外行波局部法向透射概念，从多项式外推的角度，结合一维和广义二维几何解释，推导建立了外行波为平面波时，人工边界节点位移解的计算公式。对该位移解的精度进行了时域分析，论证了该位移解收敛于数值精确解的必要条件和线性外推位移解的精度与数值精确解匹配的充要条件。在此基础上提出了一次“精确”透射的概念和相应的数学方法，实现了该位移解逼近于数值精确解的优化计算。     

Absorbing boundary conditions for dynamic analysis of fluid-saturated porous media

Based on Biot's two-phase mixture theory and the paraxial approximation, the absorbing boundary conditions in the time domain for u-w, u-U and u-p formulations are presented for the dynamic analysis of fluid-saturated porous media. These absorbing boundaries are equivalent to the viscous boundaries in the fundamental mode. The expressions for the energy ratio and reflection coefficient between the reflected and incident waves along the absorbing boundary are given. The numerical results show that the proposed absorbing boundaries can greatly suppress the spuriously reflected wave and model the far field of the system. These results also dynamic analysis of infinite fluid-saturated porous media. for the transient dynamic analysis of infinite fluid-saturated porous media.     

Local transmitting boundaries for transient elastic analysis

The aim of this paper is to investigate and develop alternative methods of analyzing problems in dynamic soil–structure-interaction (SSI). The interaction means that the amplitude of structural response is effected by additional energy dissipation through radiation and material damping in the soil. The surrounding soft soil behaves as a natural damper for a massive and stiff structure supported or embedded in it. The main focus is the major difficulty posed by such an analysis — the phenomena of waves that radiate outward from the excited structures towards infinity. In numerical calculations only a finite region of the foundation medium is analyzed and something is done to prevent the outgoing radiation waves from reflecting at the boundary region.Development of a simple and efficient finite element (FE) procedure for the solution directly in the time domain of transient SSI problems is the main concern. The central feature of the procedure is local absorbing boundaries used to render the computational domain finite. These boundaries are local in both time and space and are completely defined by a pair of symmetric stiffness and damping matrices. As the effort for implementing them is the same as for the impedance boundary condition (BC) considering the angle of incidence, standard assembly procedure can be used. Due to the local nature they also preserve the overall structure of the global equations of motion. Even though the focus is in the time domain the same equations of motions can be used to determine the solution under time-harmonic excitation directly in the frequency domain. Explicit formulae for the element matrices are included in the paper and numerical examples for transient radiation model problems to illustrate the validity and accuracy of the new procedures, are given.     

General absorbing boundary conditions for dynamic analysis of fluid-saturated porous media

General absorbing boundary conditions based on Biot's two-phase mixture theory and paraxial approximation is presented for the dynamic analysis of fluid-saturated porous media with isotropic, transverse isotropic, and anisotropic properties. For the last two cases, the equivalent Lame's constants, under conditions of uniqueness, are introduced to facilitate the analytical solutions. The numerical results show that the proposed absorbing boundary can greatly suppress spuriously reflected waves and efficiently model the far field of the system with sufficient accuracy.     

A viscous-spring transmitting boundary for cylindrical wave propagation in saturated poroelastic media

Based on the u–U formulation of Biot equation and the assumption of zero permeability coefficient, a viscous-spring transmitting boundary which is frequency independent is derived to simulate the cylindrical elastic wave propagation in unbounded saturated porous media. By this viscous-spring boundary the effective stress and pore fluid pressure on the truncated boundary of the numerical model are replaced by a set of spring, dashpot and mass elements, and its simplified form is also given. A u–U formulation FEA program is compiled and the proposed transmitting boundaries are incorporated therein. Numerical examples show that the proposed viscous-spring boundary and its simplified form can provide accurate results for cylindrical elastic wave propagation problems with low or intermediate values of permeability or frequency content. For general two dimensional wave propagation problems, spuriously reflected waves can be greatly suppressed and acceptable accuracy can still be achieved by placing the simplified boundary at relatively large distance from the wave source.     

A local time-domain transmitting boundary for simulating cylindrical elastic wave propagation in infinite media

Finite element simulation of the time-dependent wave propagation in infinite media requires enforcing the transmitting boundary to replace the truncated far-field infinite domain so as to model the effect of the wave radiation towards infinity. This paper proposed a novel local time-domain transmitting boundary for simulating the cylindrical elastic wave radiation problem. This boundary is a mechanical model consisting of the spring, dashpot and mass elements, with the auxiliary degrees of freedom introduced, which is dynamically stable and easily implemented into the commercial finite element codes. Numerical analysis of the cylindrical elastic wave radiation problem indicates that the proposed transmitting boundaries with the order N=3 for cylindrical P and SV waves and with the order N=4 for cylindrical SH wave have very high accuracy, even when the artificial boundary at wave source. The proposed transmitting boundary with order N=0 can be applied approximately to the general two-dimensional infinite elastic wave problems that contain the more complex outgoing wave fields at artificial boundary than the cylindrical waves. The plane-strain Lamb problem is analyzed with the acceptable engineering accuracy achieved. On the other hand, the proposed transmitting boundary with higher order can be a tool to localize the temporal convolution that appears in an exact time-domain transmitting boundary for the general infinite wave problems. This potential applicability is mentioned.     

土-结构动力相互作用人工边界分析及试验验证

利用有限元软件Ansys二次开发语言apdl进行二次开发,将多次透射边界(MTF边界)和弹阻边界(V-S边界)添加到软件中,从而实现对土-结构动力相互作用问题中无限域场地的模拟。以野外大比例(1/2)土-箱基-框架结构试验模型的牵引释放自由振动试验(波源问题)和爆破震动试验(散射问题)数据为基础,通过有限元数值计算结果和试验结果的对比可知,施加了人工边界的土体模型可以较好地模拟波在无限域土体中的传播,有效地减少有限元计算模拟时的计算规模。同时多次透射边界在计算精度以及散射问题中地震波输入的便捷性方面要优于弹阻边界。     

An efficient time-domain soil-structure interaction analysis based on the dynamic stiffness of an unbounded soil

An approximate method and a rigorous method are presented for the time-domain soil-structure interaction analysis, both of which use the stiffness of soil obtained numerically or experimentally in the limited frequency range. In the aseismic design of a large-scale structure, the approximate and rigorous methods are intended to be used in a preliminary analysis and a detailed study, respectively. Both methods are based on the approximation that the first few terms of the Fourier or Taylor expansions of a frequency-dependent function are used. The difference lies in the number of coefficients of each series, and also in the manner in which the coefficients can be determined. In order to demonstrate the validity of the proposed methods, a soil-structure system, whose soil has a complex profile under the foundation, is analysed. The dynamic stiffness of the soil is calculated by a 3-D hybrid approach that combines the finite element method and the boundary element method. As a result, the present methods are capable of evaluating the complexity of the soil more precisely and efficiently than conventional methods.     

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.     

On causality in dynamic response analysis by time-dependent boundary element methods

The ramifications of a particular type of causality constraint, namely so-called shadow-zones, are explored in the context of time-dependent boundary element methods. In particular, wave propagation problems in 2-D elastic soil media with a non-convex shape are analysed.     

基于高阶双渐近透射边界的重力坝-层状地基动力相互作用分析

时域高阶双渐近透射边界能够同时模拟层状介质中行波和快衰波的传播,具有很高的计算精度和计算效率.本文将高阶双渐近透射边界推广应用到多层层状地基系统弹性波传播问题的模拟,采用广义特征值分解分析该透射边界的数值稳定性,通过移谱法消除导致数值不稳定的虚假模态.将高阶双渐近透射边界以超单元的形式直接嵌入到近场有限元方程,建立了有限元-高阶双渐近透射边界时域耦合分析模型,并将其应用到重力坝-层状地基动力相互作用分析.数值算例分析结果表明,该时域耦合分析模型具有很高的精度和计算效率,适用于实际重力坝工程的地震响应分析.

     

Response of unbounded soil in scaled boundary finite‐element method

The scaled boundary finite‐element method is a powerful semi‐analytical computational procedure to calculate the dynamic stiffness of the unbounded soil at the structure–soil interface. This permits the analysis of dynamic soil–structure interaction using the substructure method. The response in the neighbouring soil can also be determined analytically. The method is extended to calculate numerically the response throughout the unbounded soil including the far field. The three‐dimensional vector‐wave equation of elasto‐dynamics is addressed. The radiation condition at infinity is satisfied exactly. By solving an eigenvalue problem, the high‐frequency limit of the dynamic stiffness is constructed to be positive definite. However, a direct determination using impedances is also possible. Solving two first‐order ordinary differential equations numerically permits the radiation condition and the boundary condition of the structure–soil interface to be satisfied sequentially, leading to the displacements in the unbounded soil. A generalization to viscoelastic material using the correspondence principle is straightforward. Alternatively, the displacements can also be calculated analytically in the far field. Good agreement of displacements along the free surface and below a prism foundation embedded in a half‐space with the results of the boundary‐element method is observed. Copyright © 2001 John Wiley & Sons, Ltd.     

透射边界高频耦合失稳机理及稳定实现——SH波动

就大型近场波动的高效数值模拟而言,稳定实现高阶人工边界是一个尚未圆满解决的问题.本文针对使用多次透射公式的SH波动集中质量有限元模拟,依据GKS定理的群速度解释,进一步阐明了人工边界与内域离散格式耦合所导致高频失稳的机理,即两者支持群速度指向内域的外行高频平面谐波,波动能量自发地从人工边界进入內域,从而导致失稳,而这类谐波是由集中质量有限元离散引入的.本文提出了消除此种耦合失稳的一种方法:通过修改有限元刚度阵来改变内域离散格式,并保证修改格式的精度不低于原有格式的精度.理论分析和数值实验表明此法能稳定实现透射边界.本文研究结果具有推广应用前景.     

ABAQUS岩土动力分析的静-动力统一人工边界研究

岩土动力学分析结果可靠性受人工边界条件及动力作用输入方式影响.在剖析ABAQUS无限元静、动力理论基础上,结合二维简单水平地基模型,研究了无限元边界生成初始地应力的计算精度及其改进方法,提出了综合考虑内、外源振动问题的改进静-动力统一人工边界.研究表明:因ABAQUS中无限元缺乏对边界节点力等效约束刚度的校核,其直接应...     

基于简谐波动模型的法向透射边界及其向任意波动的推广

基于外行简谐波动模型、多项式插值理论和具有解析表达式的外插序列,导出了N阶透射边界及其误差项与收敛的充分条件。利用外插序列的时空坐标满足同一直线方程的几何特征,完成了简谐波动N阶透射边界向任意外行波动阶透射边界的推广。依据外插序列的解析表达式,导出了后者误差项的近似表达式和收敛性的充分条件。N阶透射边界的主部与N阶多次透射公式相同,但误差项的表达式不同,精度与收敛性的分析方法与结果也不同。提出了与虚拟透射边界相关的频率放大倍数概念,认为它们是影响虚拟透射边界精度的关键因素。