Earthquake response analysis in the time domain for 2D soil–structure systems using analytical frequency‐dependent infinite elements

This paper presents a time domain method for soil–structure interaction analysis under seismic excitations. It is based on the finite element formulation incorporating analytical frequency‐dependent infinite elements for the far‐field soil region. Equivalent earthquake input forces are calculated based on the free‐field responses along the interface between the near‐ and far‐field soil regions using the fixed exterior boundary method in the frequency domain. Then, the input forces are transformed into the time domain by using inverse Fourier transform. The dynamic stiffness matrices of the far‐field soil region formulated using the analytical frequency‐dependent infinite elements in the frequency domain can be easily transformed into the corresponding matrices in the time domain. Hence, the response can be analytically computed in the time domain. A recursive procedure is proposed to compute the interaction forces along the interface and the responses of the soil–structure system in the time domain. Earthquake response analyses have been carried out on a multi‐layered half‐space and a tunnel embedded in a layered half‐space, and results are compared with those obtained by the conventional method in the frequency domain. Copyright © 2003 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.     

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

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

Diagonalization procedure for scaled boundary finite element method in modeling semi‐infinite reservoir with uniform cross‐section

To improve the ability of the scaled boundary finite element method (SBFEM) in the dynamic analysis of dam–reservoir interaction problems in the time domain, a diagonalization procedure was proposed, in which the SBFEM was used to model the reservoir with uniform cross‐section. First, SBFEM formulations in the full matrix form in the frequency and time domains were outlined to describe the semi‐infinite reservoir. No sediments and the reservoir bottom absorption were considered. Second, a generalized eigenproblem consisting of coefficient matrices of the SBFEM was constructed and analyzed to obtain corresponding eigenvalues and eigenvectors. Finally, using these eigenvalues and eigenvectors to normalize the SBFEM formulations yielded diagonal SBFEM formulations. A diagonal dynamic stiffness matrix and a diagonal dynamic mass matrix were derived. An efficient method was presented to evaluate them. In this method, no Riccati equation and Lyapunov equations needed solving and no Schur decomposition was required, which resulted in great computational costs saving. The correctness and efficiency of the diagonalization procedure were verified by numerical examples in the frequency and time domains, but the diagonalization procedure is only applicable for the SBFEM formulation whose scaling center is located at infinity. Copyright © 2009 John Wiley & Sons, Ltd.     

Precorrected FFT accelerated BEM for large‐scale transient elastodynamic analysis using frequency‐domain approach

A precorrected fast Fourier transform (pFFT) accelerated boundary element method (BEM) for large‐scale transient elastodynamic analysis is developed and described in this paper. The frequency‐domain approach is used. To overcome the ‘wrap‐around’ problem associated with the discrete Fourier transform, the exponential window method (EWM) is employed and incorporated in the frequency‐domain BEM. An improved implementation scheme of the pFFT method based on polynomial interpolation technique is developed and applied to accelerate the elastodynamic BEM. This new scheme reduces the memory required to save the convolution matrix by a factor of 8. To further improve the efficiency of the code, a newly developed linear system solver based on the induced dimension reduction method is employed. Its performance is investigated and compared with that of the well‐known GMRES. The accuracy and computational efficiency of the method are evaluated and demonstrated by three examples: a classical benchmark, a plate subject to an impact loading and a porous cube with nearly half million DOFs subject to a step traction loading. Both analytical and experimental results are employed to validate the method. It has been found that the EWM can effectively resolve the wrap‐around problem and accurate time responses for an arbitrarily chosen time period can be obtained. Copyright © 2011 John Wiley & Sons, Ltd.     

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.     

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.     

An improved continued‐fraction‐based high‐order transmitting boundary for time‐domain analyses in unbounded domains

A high‐order local transmitting boundary to model the propagation of acoustic or elastic, scalar or vector‐valued waves in unbounded domains of arbitrary geometry is proposed. It is based on an improved continued‐fraction solution of the dynamic stiffness matrix of an unbounded medium. The coefficient matrices of the continued‐fraction expansion are determined recursively from the scaled boundary finite element equation in dynamic stiffness. They are normalised using a matrix‐valued scaling factor, which is chosen such that the robustness of the numerical procedure is improved. The resulting continued‐fraction solution is suitable for systems with many DOFs. It converges over the whole frequency range with increasing order of expansion and leads to numerically more robust formulations in the frequency domain and time domain for arbitrarily high orders of approximation and large‐scale systems. Introducing auxiliary variables, the continued‐fraction solution is expressed as a system of linear equations in iω in the frequency domain. In the time domain, this corresponds to an equation of motion with symmetric, banded and frequency‐independent coefficient matrices. It can be coupled seamlessly with finite elements. Standard procedures in structural dynamics are directly applicable in the frequency and time domains. Analytical and numerical examples demonstrate the superiority of the proposed method to an existing approach and its suitability for time‐domain simulations of large‐scale systems. Copyright © 2011 John Wiley & Sons, Ltd.     

Frequency‐domain analysis of time‐integration methods for semidiscrete finite element equations—part I: Parabolic problems

Time‐integration methods for semidiscrete equations emanating from parabolic differential equations are analysed in the frequency domain. The discrete‐time transfer functions of three popular methods are derived, and subsequently the forced response characteristics of single modes are studied in the frequency domain. To enable consistent comparison of the frequency responses of different algorithms, three characteristic numbers are identified. Frequency responses and L₂‐norms of the phase and magnitude errors are compared for the three time‐integration algorithms. The examples demonstrate that frequency‐domain analysis provides substantial insight into the time‐domain properties of time‐integration algorithms. Copyright © 2001 John Wiley & Sons, Ltd.     

A boundary condition in Padé series for frequency‐domain solution of wave propagation in unbounded domains

A boundary condition satisfying the radiation condition at infinity is frequently required in the numerical simulation of wave propagation in an unbounded domain. In a frequency domain analysis using finite elements, this boundary condition can be represented by the dynamic stiffness matrix of the unbounded domain defined on its boundary. A method for determining a Padé series of the dynamic stiffness matrix is proposed in this paper. This method starts from the scaled boundary finite‐element equation, which is a system of ordinary differential equations obtained by discretizing the boundary only. The coefficients of the Padé series are obtained directly from the ordinary differential equations, which are not actually solved for the dynamic stiffness matrix. The high rate of convergence of the Padé series with increasing order is demonstrated numerically. This technique is applicable to scalar waves and elastic vector waves propagating in anisotropic unbounded domains of irregular geometry. It can be combined seamlessly with standard finite elements. Copyright © 2006 John Wiley & Sons, Ltd.     

Scaled boundary finite‐element analysis of a non‐homogeneous elastic half‐space

As a result of stresses experienced during and after the deposition phase, a soil strata of uniform material generally exhibits an increase in elastic stiffness with depth. The immediate settlement of foundations on deep soil deposits and the resultant stress state within the soil mass may be most accurately calculated if this increase in stiffness with depth is taken into account. This paper presents an axisymmetric formulation of the scaled boundary finite‐element method and incorporates non‐homogeneous elasticity into the method. The variation of Young's modulus (E) with depth (z) is assumed to take the form E=m_Ez^α, where m_E is a constant and αis the non‐homogeneity parameter. Results are presented and compared to analytical solutions for the settlement profiles of rigid and flexible circular footings on an elastic half‐space, under pure vertical load with αvarying between zero and one, and an example demonstrating the versatility and practicality of the method is also presented. Known analytical solutions are accurately represented and new insight regarding displacement fields in a non‐homogeneous elastic half‐space is gained. Copyright © 2003 John Wiley & Sons, Ltd.     

A continued‐fraction‐based high‐order transmitting boundary for wave propagation in unbounded domains of arbitrary geometry

A high‐order local transmitting boundary is developed to model the propagation of elastic waves in unbounded domains. This transmitting boundary is applicable to scalar and vector waves, to unbounded domains of arbitrary geometry and to anisotropic materials. The formulation is based on a continued‐fraction solution of the dynamic‐stiffness matrix of an unbounded domain. The coefficient matrices of the continued fraction are determined recursively from the scaled boundary finite element equation in dynamic stiffness. The solution converges rapidly over the whole frequency range as the order of the continued fraction increases. Using the continued‐fraction solution and introducing auxiliary variables, a high‐order local transmitting boundary is formulated as an equation of motion with symmetric and frequency‐independent coefficient matrices. It can be coupled seamlessly with finite elements. Standard procedures in structural dynamics are directly applicable for evaluating the response in the frequency and time domains. Analytical and numerical examples demonstrate the high rate of convergence and efficiency of this high‐order local transmitting boundary. Copyright © 2007 John Wiley & Sons, Ltd.     

Frequency‐domain analysis of time‐integration methods for semidiscrete finite element equations—part II: Hyperbolic and parabolic–hyperbolic problems

Time‐integration methods for semidiscrete finite element equations of hyperbolic and parabolic– hyperbolic types are analysed in the frequency domain. The discrete‐time transfer functions of six popular methods are derived, and subsequently the forced response characteristics of single modes are studied in the frequency domain. Three characteristic numbers are derived which eliminate the parameter dependence of the frequency responses. Frequency responses and L₂‐norms of the phase and magnitude errors are calculated, and comparisons are given of the methods. As shown; the frequency‐domain analysis explains all time‐domain properties of the methods, and gives more insight into the subject. Copyright © 2001 John Wiley & Sons, Ltd.     

An efficient and high‐order frequency‐domain approach for transient acoustic–structure interactions in three dimensions

This paper is concerned with numerical solution of the transient acoustic–structure interaction problems in three dimensions. An efficient and higher‐order method is proposed with a combination of the exponential window technique and a fast and accurate boundary integral equation solver in the frequency‐domain. The exponential window applied to the acoustic–structure system yields an artificial damping to the system, which eliminates the wrap‐around errors brought by the discrete Fourier transform. The frequency‐domain boundary integral equation approach relies on accurate evaluations of relevant singular integrals and fast computation of nonsingular integrals via the method of equivalent source representations and the fast Fourier transform. Numerical studies are presented to demonstrate the accuracy and efficiency of the method. Copyright © 2016 John Wiley & Sons, Ltd.     

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.     

A Frobenius solution to the scaled boundary finite element equations in frequency domain for bounded media

The scaled boundary finite element method (FEM) is a recently developed semi‐analytical numerical approach combining advantages of the FEM and the boundary element method. Although for elastostatics, the governing homogeneous differential equations in the radial co‐ordinate can be solved analytically without much effort, an analytical solution to the non‐homogeneous differential equations in frequency domain for elastodynamics has so far only been obtained by a rather tedious series‐expansion procedure. This paper develops a much simpler procedure to obtain such an analytical solution by increasing the number of power series in the solution until the required accuracy is achieved. The procedure is applied to an extensive study of the steady‐state frequency response of a square plate subjected to harmonic excitation. Comparison of the results with those obtained using ABAQUS shows that the new method is as accurate as a detailed finite element model in calculating steady‐state responses for a wide range of frequencies using only a fraction of the degrees of freedom required in the latter. Copyright © 2006 John Wiley & Sons, Ltd.     

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

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

A time‐integration method for the viscoelastic–viscoplastic analyses of polymers and finite element implementation

The present study introduces a time‐integration algorithm for solving a non‐linear viscoelastic–viscoplastic (VE–VP) constitutive equation of isotropic polymers. The material parameters in the constitutive models are stress dependent. The algorithm is derived based on an implicit time‐integration method (Computational Inelasticity. Springer: New York, 1998) within a general displacement‐based finite element (FE) analysis and suitable for small deformation gradient problems. Schapery's integral model is used for the VE responses, while the VP component follows the Perzyna model having an overstress function. A recursive‐iterative method (Int. J. Numer. Meth. Engng 2004; 59 :25–45) is employed and modified to solve the VE–VP constitutive equation. An iterative procedure with predictor–corrector steps is added to the recursive integration method. A residual vector is defined for the incremental total strain and the magnitude of the incremental VP strain. A consistent tangent stiffness matrix, as previously discussed in Ju (J. Eng. Mech. 1990; 116 :1764–1779) and Simo and Hughes (Computational Inelasticity. Springer: New York, 1998), is also formulated to improve convergence and avoid divergence. Available experimental data on time‐dependent and inelastic responses of high‐density polyethylene are used to verify the current numerical algorithm. The time‐integration scheme is examined in terms of its computational efficiency and accuracy. Numerical FE analyses of microstructural responses of polyethylene reinforced with elastic particle are also presented. Copyright © 2009 John Wiley & Sons, Ltd.     

Interface dynamic stiffness matrix approach for three‐dimensional transient multi‐region boundary element analysis

A novel substructuring method is developed for the coupling of boundary element and finite element subdomains in order to model three‐dimensional multi‐region elastodynamic problems in the time domain. The proposed procedure is based on the interface stiffness matrix approach for static multi‐region problems using variational principles together with the concept of Duhamel integrals. Unit impulses are applied at the boundary of each region in order to evaluate the impulse response matrices of the Duhamel (convolution) integrals. Although the method is not restricted to a special discretization technique, the regions are discretized using the boundary element method combined with the convolution quadrature method. This results in a time‐domain methodology with the advantages of performing computations in the Laplace domain, which produces very accurate and stable results as verified on test examples. In addition, the assembly of the boundary element regions and the coupling to finite elements are greatly simplified and more efficient. Finally, practical applications in the area of soil–structure interaction and tunneling problems are shown. Copyright © 2009 John Wiley & Sons, Ltd.     

Time‐decomposed parallel time‐integrators: theory and feasibility studies for fluid,structure, and fluid–structure applications

A methodology for squeezing the most out of massively parallel processors when solving partial differential evolution equations by implicit schemes is presented. Its key components include a preferred implicit time‐integrator, a decomposition of the time‐domain into time‐slices, independent time‐integrations in each time‐slice of the semi‐discrete equations, and Newton‐type iterations on a coarse time‐grid. Hence, this methodology parallelizes the time‐loop of a time‐dependent partial differential equation solver without interfering with its sequential or parallel space‐computations. It is particularly interesting for time‐dependent problems with a few degrees of freedom such as those arising in robotics and protein folding applications, where the opportunities for parallelization over the degrees of freedom are limited. Error and stability analyses of the proposed parallel methodology are performed for first‐ and second‐order hyperbolic problems. Its feasibility and impact on reducing the solution time below what is attainable by methods which address only parallelism in the space‐domain are highlighted for fluid, structure, and coupled fluid–structure model problems. Copyright © 2003 John Wiley & Sons, Ltd.