Time‐domain soil–structure interaction analysis in two‐dimensional medium based on analytical frequency‐dependent infinite elements

A direct method for soil–structure interaction analysis in two‐dimensional medium is presented in time domain, which is based on the transformation of the analytical frequency‐dependent dynamic stiffness matrix. The present dynamic stiffness matrix for the far‐field region is constructed by assembling stiffness matrices of the analytical frequency‐dependent dynamic infinite elements, so that the equation of motion can be analytically transformed into the time‐domain equation. An efficient procedure is devised to evaluate the dynamic responses in time domain. Verification of the present formulation is carried out by comparing the compliances for a strip foundation on a homogeneous and layered half‐spaces with those obtained by other methods. Numerical analyses are also carried out for the transient responses of an elastic block and tunnel in a homogeneous and a layered half‐space. The comparisons with those by other approaches indicate that the proposed time‐domain method for soil–structure interaction analysis gives good solutions. Copyright © 2000 John Wiley & Sons, Ltd.     

A modified scaled boundary finite element method for three‐dimensional dynamic soil‐structure interaction in layered soil

This paper is devoted to the analysis of elastodynamic problems in 3D‐layered systems which are unbounded in the horizontal direction. For this purpose, a finite element model of the near field is coupled to a scaled boundary finite element model (SBFEM) of the far field. The SBFEM is originally based on describing the geometry of a half‐space or full‐space domain by scaling the geometry of the near field / far field interface using a radial coordinate. A modified form of the SBFEM for waves in a 2D layer is also available. None of these existing formulations can be used to describe a 3D‐layered medium. In this paper, a modified SBFEM for the analysis of 3D‐layered continua is derived. Based on the use of a scaling line instead of a scaling centre, a suitable scaled boundary transformation is proposed. The derivation of the corresponding scaled boundary finite element (SBFE) equations in displacement and stiffness is presented in detail. The latter is a nonlinear differential equation with respect to the radial coordinate, which has to be solved numerically for each excitation frequency considered in the analysis. Various numerical examples demonstrate the accuracy of the new method and its correct implementation. These include rigid circular and square foundations embedded in or resting on the surface of layered homogeneous or inhomogeneous 3D soil deposits over rigid bedrock. Hysteretic damping is assumed in some cases. The dynamic stiffness coefficients calculated using the proposed method are compared with analytical solutions or existing highly accurate numerical results. Copyright © 2011 John Wiley & Sons, Ltd.     

A 2.5D finite/infinite element approach for modelling visco‐elastic bodies subjected to moving loads

The objective of this study is to propose a 2.5D finite/infinite element procedure for dealing with the ground vibrations induced by moving loads. Besides the two in‐plane degrees of freedom (DOFs) per node conventionally used for plane strain elements, an extra DOF is introduced to account for the out‐of‐plane wave transmission. The profile of the half‐space is divided into a near field and a semi‐infinite far field. The near field containing loads and irregular structures is simulated by the finite elements, while the far field covering the soils extending to infinity by the infinite elements with due account taken of the radiation effects for moving loads. Enhanced by the automated mesh expansion procedure proposed previously by the writers, the far field impedances for all the lower frequencies are generated repetitively from the mesh created for the highest frequency considered. Finally, the accuracy of the proposed method is verified through comparison with a number of analytical solutions. Copyright © 2001 John Wiley & Sons, Ltd.     

Seismic response to a prescribed seismogram of a body embedded in a multilayered half space

While the actual problem is composed of an active fault surface, a soil site and a body embedded at that site; the proposed method provides an alternative smaller linear problem by replacing the propagating rupture on the fault surface by a fictitious focal point and a seismograph station in the vicinity of the given soil site. The Green's function for each of three fundamental problems of isotropic elastic and viscoelastic spaces undergoing harmonic vibration is derived. Infinite elements are adopted in the far field, and finite elements in the near field. The three fundamental problem solutions are used as the shape functions of infinite element nodal lines. The three concentrated orthogonal force components at the focal point are determined in such a way that the Fourier transforms of the three orthogonal acceleration components measured at a seismograph station are checked. For seismic analysis of a finite embedded body, consider the differential between the actual system and the seismic free field, which is the embedding half space without any embedment and being excited by the fictitious focal point forces. All along the analysis has been carried out in the frequency domain. An appropriate inverse Fourier transform algorithm will properly yield all results as time functions.     

Rayleigh wave correction for the BEM analysis of two‐dimensional elastodynamic problems in a half‐space

A simple, elegant approach is proposed to correct the error introduced by the truncation of the infinite boundary in the BEM modelling of two‐dimensional wave propagation problems in elastic half‐spaces. The proposed method exploits the knowledge of the far‐field asymptotic behaviour of the solution to adequately correct the BEM displacement system matrix for the truncated problem to account for the contribution of the omitted part of the boundary. The reciprocal theorem of elastodynamics is used for a convenient computation of this contribution involving the same boundary integrals that form the original BEM system. The method is formulated for a two‐dimensional homogeneous, isotropic, linearly elastic half‐space and its implementation in a frequency domain boundary element scheme is discussed in some detail. The formulation is then validated for a free Rayleigh pulse travelling on a half‐space and successfully tested for a benchmark problem with a known approximation to the analytical solution. Copyright © 2004 John Wiley & Sons, Ltd.     

Dynamic impedance of pile and caisson foundations

Abstract

The dynamic impedance of a pile or caisson foundation embedded in a homogeneous visco elastic half‐space was evaluated using the hybrid method. The soil‐foundation system is partitioned into a near‐field and a far‐field by choosing a cylindrical interface passing through the soil region very close to the foundation. The near‐field is modeled by the finite element method while the far‐field is characterized by a frequency‐dependent impedance matrix through the continuum approach. From the results presented, it is shown that the proposed method is very effective and can be widely used for parameter studies in engineering applications.     

SV波超临界角斜入射时层状地基地震动输入在ABAQUS中的实现

有限元数值模拟中,超临界角斜入射SV波作用下层状地基的地震动输入是一个亟待解决的问题。该文采用层状地基频域内精确的动力刚度矩阵(即频域刚度矩阵法)推导了SV波任意角度斜入射下的地震动输入等效节点力计算公式,通过ABAQUS有限元软件模拟SV波斜入射下均匀半空间、层状地基的地震波场,探讨了基于频域刚度矩阵法的层状地基任意角度斜入射地震动输入应用于ABAQUS有限元软件的有效性和准确性。结果表明,采用频域刚度矩阵法可以在ABAQUS中实现层状地基任意角度斜入射SV波地震动输入,且方法具有很高的计算精度,尤其对于SV波超临界角入射情况,有限元数值模拟很好地再现了均匀半空间地表质点的椭圆型运动轨迹和层状地基的行波特点。在此基础上进一步将频域刚度矩阵法与等效线性化方法相结合,解决了考虑土体非线性的层状地基在任意角度平面波入射条件下的地震动输入问题。     

多层介质硬币形(Penny-Shaped)交界裂纹的弹性波散射

本文研究多层介质硬币形交界裂纹的弹性波散射.文中采用Hankel积分变换,得到了含有硬币形交界裂纹多层介质模型的散射波传递矩阵,并将散射问题为转化求解矩阵形式的对偶积分方程.作为特例,文中给出了单一弹性层与半空间的硬币形交界裂纹的弹性波散射远场模式,并计算了几组不同弹性常数组合情形下的远场模式的幅频特性曲线,其结果表明有共振峰存在.     

Transient analysis of three-dimensional dynamic interaction between multilayered soil and rigid foundation

The study of dynamic soil-structure interaction is significant to civil engineering applications, such as machine foundation vibration, traffic-induced vibration, and seismic dynamic response. The scaled boundary finite element method (SBFEM) is a semi-analytical algorithm, which is used to solve the dynamic response of a three-dimensional infinite soil. It can automatically satisfy the radiation boundary condition at infinity. Based on the dynamic stiffness matrix equation obtained by the modified SBFEM, a continued fraction algorithm is proposed to solve the dynamic stiffness matrix of layered soil in the frequency-domain. Then, the SBFEM was coupled with the finite element method (FEM) at the interface to solve the dynamic stiffness matrices of the rigid surface/buried foundation. Finally, the mixed-variable algorithm was used to solve the three-dimensional transient dynamic response of the foundation in the time domain. Numerical examples were performed to verify the accuracy of the proposed algorithm in solving the dynamic stiffness matrix of the infinite domain in the frequency domain and the dynamic transient displacement response of the foundation in the time domain. Compared with the previous numerical integration technique, the dynamic stiffness matrix in the frequency domain calculated by using the proposed algorithm has higher accuracy and higher efficiency.     

层状饱和半空间中沉积谷地对斜入射平面P1波的三维散射

在作者给出层状饱和场地三维精确动力刚度矩阵和层状饱和半空间中移动荷载动力格林函数基础上,采用间接边界元方法在频域内求解了层状流体饱和场地中沉积谷地对斜入射平面P1波的三维散射问题。该方法的特点在于虚拟移动均布荷载和斜线孔隙水压可以直接施加在沉积与层状饱和半空间交界面而不存在奇异性。该文通过与已有结果的比较验证了方法的正确性,并以均匀饱和半空间和弹性基岩上单一饱和土层中沉积谷地为例进行了数值计算分析。研究表明,沉积谷地对平面P1波的三维散射与二维散射之间存在本质差别,入射角度、孔隙率、饱和土层刚度和饱和土层厚度等参数对沉积谷地附近动力响应有着显著影响。     

Improper integrals for infinite elements in vibration of half spaces

In vibration of elastic and viscoelastic multilayered half spaces, the far field of each half space is discretized into infinite elements. Far-field displacement functions of three fundamental problems are used as the shape functions of appropriate infinite element nodal lines. Each fundamental problem solution is obtained as a linear combination of discrete waves. As a result, the mass and stiffness matrices of infinite elements involve improper integrals, which are evaluated by three schemes proposed herein. Different schemes are good for different ranges of the complex argument of the integrals.     

Explicit finite element artificial boundary scheme for transient scalar waves in two‐dimensional unbounded waveguide

To simulate the transient scalar wave propagation in a two‐dimensional unbounded waveguide, an explicit finite element artificial boundary scheme is proposed, which couples the standard dynamic finite element method for complex near field and a high‐order accurate artificial boundary condition (ABC) for simple far field. An exact dynamic‐stiffness ABC that is global in space and time is constructed. A temporal localization method is developed, which consists of the rational function approximation in the frequency domain and the auxiliary variable realization into time domain. This method is applied to the dynamic‐stiffness ABC to result in a high‐order accurate ABC that is local in time but global in space. By discretizing the high‐order accurate ABC along artificial boundary and coupling the result with the standard lumped‐mass finite element equation of near field, a coupled dynamic equation is obtained, which is a symmetric system of purely second‐order ordinary differential equations in time with the diagonal mass and non‐diagonal damping matrices. A new explicit time integration algorithm in structural dynamics is used to solve this equation. Numerical examples are given to demonstrate the effectiveness of the proposed scheme. Copyright © 2011 John Wiley & Sons, Ltd.     

The fractal finite element method for unbounded problems

The fractal finite element method, previously developed for stress intensity factor calculation for crack problems in fracture mechanics, is extended to analyse some unbounded problems in half space. The concepts of geometrical similarity and two‐level finite element mesh are applied to generate an infinite number of self‐similar layers in the far field with a similarity ratio bigger than one; that is, one layer is bigger than the next in size but of the same shape. Only conventional finite elements are used and no new elements are generated. The global interpolating functions in the form of a truncated infinite series are employed to transform the infinite number of nodal variables to a small number of unknown coefficients associated with the global interpolating functions. Taking the advantage of geometrical similarity, transformation for one layer is enough because the relevant entries of the transformed matrix after assembling all layers are infinite geometric series of the similarity ratio and can be summed analytically. Accurate nodal displacements are obtained as shown in the numerical examples. Copyright © 2004 John Wiley & Sons, Ltd.     

Simulation of transient heat conduction using one‐dimensional mapped infinite elements

Many engineering problems exist in physical domains that can be said to be infinitely large. A common problem in the simulation of these unbounded domains is that a balance must be met between a practically sized mesh and the accuracy of the solution. In transient applications, developing an appropriate mesh size becomes increasingly difficult as time marches forward. The concept of the infinite element was introduced and implemented for elliptic and for parabolic problems using exponential decay functions. This paper presents a different methodology for modeling transient heat conduction using a simplified mesh consisting of only two‐node, one‐dimensional infinite elements for diffusion into an unbounded domain and is shown to be applicable for multi‐dimensional problems. A brief review of infinite elements applied to static and transient problems is presented. A transient infinite element is presented in which the element length is time‐dependent such that it provides the optimal solution at each time step. The element is validated against the exact solution for constant surface heat flux into an infinite half‐space and then applied to the problem of heat loss in thermal reservoirs. The methodology presented accurately models these phenomena and presents an alternative methodology for modeling heat loss in thermal reservoirs. Copyright © 2010 John Wiley & Sons, Ltd.     

Earthquake‐induced hydrodynamic pressures on a 3D rigid dam–reservoir system using DRBEM and a radiation matrix

The dual reciprocity method is applied to determine the hydrodynamic pressure distribution in a three‐dimensional dam–reservoir system subjected to earthquake excitation. The reservoir domain is idealized as a finite region of irregular geometry adjacent to the dam followed by an infinite domain of uniform cross‐section in the upstream direction. The reservoir hydrodynamic pressure response is governed by the Helmholtz equation subject to free surface, dam–reservoir interface, absorbing bottom/banks and radiation boundary conditions. A three‐dimensional (3D) dual reciprocity model is developed to determine the hydrodynamic pressure in the finite reservoir domain. A radiation matrix is developed and introduced to relate the hydrodynamic pressure and its normal derivative on the interface between the finite and infinite domains. The three‐dimensional radiation model used is developed by applying a two‐dimensional dual reciprocity formulation along the interface of the finite and infinite reservoirs together with a continuum solution in the upstream direction of the infinite domain. The model is compared for the hydrodynamic response of a three‐dimensional rectangular reservoir and found to be in excellent agreement with results obtained from a model based on the analytical formulation existing in the literature. Copyright © 2003 John Wiley & Sons, Ltd.     

Boundary elements for half‐space problems via fundamental solutions: A three‐dimensional analysis

An efficient solution technique is proposed for the three‐dimensional boundary element modelling of half‐space problems. The proposed technique uses alternative fundamental solutions of the half‐space (Mindlin's solutions for isotropic case) and full‐space (Kelvin's solutions) problems. Three‐dimensional infinite boundary elements are frequently employed when the stresses at the internal points are required to be evaluated. In contrast to the published works, the strongly singular line integrals are avoided in the proposed solution technique, while the discretization of infinite elements is independent of the finite boundary elements. This algorithm also leads to a better numerical accuracy while the computational time is reduced. Illustrative numerical examples for typical isotropic and transversely isotropichalf‐space problems demonstrate the potential applications of the proposed formulations. Incidentally, the results of the illustrative examples also provide a parametric study for the imperfect contact problem. Copyright © 2001 John Wiley & Sons, Ltd.     

A super‐element for crack analysis in the time domain

A super‐element for the dynamic analysis of two‐dimensional crack problems is developed based on the scaled boundary finite‐element method. The boundary of the super‐element containing a crack tip is discretized with line elements. The governing partial differential equations formulated in the scaled boundary co‐ordinates are transformed to ordinary differential equations in the frequency domain by applying the Galerkin's weighted residual technique. The displacements in the radial direction from the crack tip to a point on the boundary are solved analytically without any a priori assumption. The scaled boundary finite‐element formulation leads to symmetric static stiffness and mass matrices. The super‐element can be coupled seamlessly with standard finite elements. The transient response is evaluated directly in the time domain using a standard time‐integration scheme. The stress field, including the singularity around the crack tip, is expressed semi‐analytically. The stress intensity factors are evaluated without directly addressing singular functions, as the limit in their definitions is performed analytically. Copyright © 2004 John Wiley & Sons, Ltd.     

Transient infinite elements for contaminant transport problems

A time-dependent infinite element which can be used to simulate contaminant transport problems in infinite media is presented in this paper. Since this transient infinite element is constructed in a global co-ordinate system instead of a local one, a closed-form solution for the property matrices of the element has been derived from an advection–diffusion/dispersion problem in a homogeneous, anisotropic infinite medium. The numerical results from the present transient infinite element have excellent agreements with the corresponding analytical solutions. Compared with the previous infinite elements, the present infinite element has the following special characteristics: (1) both space and time variables were explicitly considered in the formulation; (2) its property matrices were expressed in a closed form; (3) it can be used to represent the far field of a mass/contaminant transport problem in a homogeneous, anisotropic infinite medium; (4) it was constructed in a global co-ordinate system. Therefore, it is highly recommended that the transient infinite element be used for the numerical simulation of contaminant transport problems in infinite media.     

Dynamic stability of non-linear shells of revolution under consideration of the fluid–soil–structure interaction

The dynamic behaviour of liquid-filled shells of revolution is investigated considering the soil–structure interaction and the fluid–structure interaction, respectively. In the circumferential direction the loads and variables are approximated by Fourier series. The shell is modelled through shell ring elements including non-linear behaviour, coupled with isoparametric continuum ring elements and special infinite elements for the soil and isoparametric pressure ring elements for the fluid. Transient loadings like earthquake excitation and the non-linearities of the shell and the soil require an analysis in the time domain. To reduce the size of the problem, linear parts of the system are condensed by the substructure technique. The soil region is divided into two parts, a near field permitting non-linearities like plastification or uplifting of the shell, and a far field for the treatment of radiation of energy. The boundary conditions for the shell footing have a strong influence on the distribution of the axial membrane forces and, hence, on the stability limit, which is mostly governed by plastic collapse and caused by the dynamically activated pressure acting on the tank wall. It is shown how the soil properties influence the dynamic stability of the shell under harmonic excitation and under realistic earthquake motion.     

高速列车荷载作用下高架桥和地基振动分析

提出了一种准解析方法来分析高速列车运行荷载作用下高架桥的动力响应。着重考虑了桥梁和地基的动力相互作用。通过动力子结构方法把研究对象分成两部分，一个是列车荷载作用下三维高架桥的有限元振动模型，另一个是基于傅里叶级数展开的轴对称群桩基础与周围分层地基动力相互作用模型，两者通过桩基承台节点处的连续条件进行结合。用薄层单元构建了应力波的透射边界条件来模拟远场地基对近场有限元区域的影响作用。用阻抗函数来表示群桩基础对上部桥梁结构的支撑作用。通过数值计算考察了软弱地基上新干线高架桥在高速列车荷载作用下的振动特性，分析了列车轮轴荷重，运行速度和群桩基础等因素对高架桥振动的影响；本分析模型和方法在计算上具有很高的效率。同时根据计算结果与现场实测的对比说明了本方法的可靠性。