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.     

A high‐order approach for modelling transient wave propagation problems using the scaled boundary finite element method

A high‐order time‐domain approach for wave propagation in bounded and unbounded domains is proposed. It is based on the scaled boundary FEM, which excels in modelling unbounded domains and singularities. The dynamic stiffness matrices of bounded and unbounded domains are expressed as continued‐fraction expansions, which leads to accurate results with only about three terms per wavelength. An improved continued‐fraction approach for bounded domains is proposed, which yields numerically more robust time‐domain formulations. The coefficient matrices of the corresponding continued‐fraction expansion are determined recursively. The resulting solution is suitable for systems with many DOFs as it converges over the whole frequency range, even for high orders of expansion. A scheme for coupling the proposed improved high‐order time‐domain formulation for bounded domains with a high‐order transmitting boundary suggested previously is also proposed. In the time‐domain, the coupled model corresponds to equations of motion with symmetric, banded and frequency‐independent coefficient matrices, which can be solved efficiently using standard time‐integration schemes. Numerical examples for modal and time‐domain analysis are presented to demonstrate the increased robustness, efficiency and accuracy of the proposed method. Copyright © 2013 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.     

High‐order doubly asymptotic open boundaries for scalar wave equation

High‐order doubly asymptotic open boundaries are developed for transient analyses of scalar waves propagating in a semi‐infinite layer with a constant depth and a circular cavity in a full‐plane. The open boundaries are derived in the frequency domain as doubly asymptotic continued fraction solutions of the dynamic stiffness of the unbounded domains. Each term of the continued fraction is a linear function of the excitation frequency. The constants of the continued fraction solutions are determined recursively. The continued fraction solution is expressed in the time domain as ordinary differential equations, which can be solved by standard time‐stepping schemes. No parameters other than the orders of the low‐ and high‐frequency expansions need to be selected by users. Numerical experiments demonstrate that evanescent waves and long‐time (low‐frequency) responses are simulated accurately. In comparison with singly asymptotic open boundaries, significant gain in accuracy is achieved at no additional computational cost. Copyright © 2009 John Wiley & Sons, Ltd.     

The scaled boundary finite element method in structural dynamics

The scaled boundary finite element method is extended to solve problems of structural dynamics. The dynamic stiffness matrix of a bounded (finite) domain is obtained as a continued fraction solution for the scaled boundary finite element equation. The inertial effect at high frequencies is modeled by high‐order terms of the continued fraction without introducing an internal mesh. By using this solution and introducing auxiliary variables, the equation of motion of the bounded domain is expressed in high‐order static stiffness and mass matrices. Standard procedures in structural dynamics can be applied to perform modal analyses and transient response analyses directly in the time domain. Numerical examples for modal and direct time‐domain analyses are presented. Rapid convergence is observed as the order of continued fraction increases. A guideline for selecting the order of continued fraction is proposed and validated. High computational efficiency is demonstrated for problems with stress singularity. Copyright © 2008 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.     

求解无限域动力刚度矩阵的双渐近算法

采用连分式算法可以有效地求解无限域动力刚度表示的比例边界有限元方程, 它具有收敛范围广、收敛速度快等优点. 该文在高频渐近连分式算法的基础上考虑了低频渐近, 发展了一种针对矢量波动方程的双渐近算法. 随着展开阶数的增加, 双渐近算法可以在全频域范围内快速逼近准确解. 引入了系数矩阵?X(i)来增强连分式算法的数值稳定性. 通过在高频极限、低频极限时满足动力刚度表示的比例边界有限元方程, 建立了递推关系以求得动力刚度矩阵. 通过二维半无限楔形体、三维均质弹性半空间数值算例表明, 双渐近算法比单渐近算法更稳定、优越.     

离散高阶Higdon-like透射边界

基于比例边界有限元法(SBFEM)半离散思想和Higdon透射微分算子提出了一种用于模拟二维层状介质标量波传播的高效离散高阶Higdon-like透射边界。对无限介质边界进行迦辽金有限元离散后,描述标量波的偏微分方程转换为局部坐标系下半离散矩阵方程组;然后使用高阶Higdon透射算子和辅助变量,在时域内得到了一个阶数不超过2阶的离散高阶透射边界。透射边界是由一组常微分方程构成,可以采用通常的时步积分方法求解,它在截断边界上非局部,在时间域局部。算例表明:该文提出的透射边界的计算精度可以随着辅助变量的增加而提高,但计算量却呈线性化增加,因而计算效率较全局方法有了显著提高。另外,由于该文的边界条件是直接建立在离散节点上的,所以它很方便与近场有限单元法耦合。     

An optimal high‐order non‐reflecting finite element scheme for wave scattering problems

A new finite element scheme is proposed for the numerical solution of time‐harmonic wave scattering problems in unbounded domains. The infinite domain in truncated via an artificial boundary ?? which encloses a finite computational domain Ω. On ?? a local high‐order non‐reflecting boundary condition (NRBC) is applied which is constructed to be optimal in a certain sense. This NRBC is implemented in a special way, by using auxiliary variables along the boundary ??, so that it involves no high‐order derivatives regardless of its order. The order of the scheme is simply an input parameter, and it may be arbitrarily high. This leads to a symmetric finite element formulation where standard C⁰ finite elements are used in Ω. The performance of the method is demonstrated via numerical examples, and it is compared to other NRBC‐based schemes. The method is shown to be highly accurate and stable, and to lead to a well‐conditioned matrix problem. Copyright © 2002 John Wiley & Sons, Ltd.     

An efficient complex‐frequency shifted‐perfectly matched layer for second‐order acoustic wave equation

In the context of simulations of wave propagations in unbounded domain, absorbing boundary conditions are often used to truncate the simulation domain to a finite space. Perfectly matched layer (PML) has proven to be an excellent absorbing boundary conditions. However, as this technique was primarily designed for the first‐order equation system, it cannot be applied to the second‐order equation system directly. In this paper, based on a complex‐coordinate stretching technique, we developed a novel, efficient auxiliary‐differential equation form of the complex‐frequency shifted‐PML for the second‐order equation system. This facilitates the use of complex‐frequency shifted‐PML in acoustic simulations based upon wave equations of second‐order form. Compared with previous state‐of‐the‐art methods, the proposed one has the advantage of simpler implementation. It is an unsplit‐field scheme that can be extended to higher‐order discretization schemes conveniently. Numerical results from both homogeneous and heterogeneous computational domains are provided to illustrate the validity of the method. Copyright © 2013 John Wiley & Sons, Ltd.     

Accurate H-shaped absorbing boundary condition in frequency domain for scalar wave propagation in layered half-space

An accurate absorbing boundary condition (ABC) is developed in frequency domain for finite element analysis of scalar wave propagation in unbounded layered half-space. The proposed ABC is H-shaped line that consists of two parts: a new ABC at horizontal bottom boundary of finite domain to replace semiinfinite strip below horizontal boundary and between two vertical boundaries, and a general consistent ABC at vertical lateral boundary to replace semiinfinite layered half-space outside vertical boundary. The key point for constructing the ABC is that a new continued fraction (CF) is presented to expand dynamic stiffness of underlying half-space, and the CF-based stress-displacement relationship is then transformed into an auxiliary variable system with square of horizontal wavenumber. The ABC has only one undetermined real parameter that is the CF-order independent of frequency and incidence angle of propagating outgoing waves. The parameter can be chosen relatively small value to achieve an accurate ABC. Moreover, the ABC can couple seamlessly with finite element method of finite domain. The finite domain can be chosen very small size due to high accuracy of the ABC. Numerical examples are finally given to demonstrate the effectiveness of the ABC.     

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.     

High‐order Higdon‐like boundary conditions for exterior transient wave problems

Recently developed non‐reflecting boundary conditions are applied for exterior time‐dependent wave problems in unbounded domains. The linear time‐dependent wave equation, with or without a dispersive term, is considered in an infinite domain. The infinite domain is truncated via an artificial boundary ??, and a high‐order non‐reflecting boundary condition (NRBC) is imposed on ??. Then the problem is solved numerically in the finite domain bounded by ??. The new boundary scheme is based on a reformulation of the sequence of NRBCs proposed by Higdon. We consider here two reformulations: one that involves high‐order derivatives with a special discretization scheme, and another that does not involve any high derivatives beyond second order. The latter formulation is made possible by introducing special auxiliary variables on ??. In both formulations the new NRBCs can easily be used up to any desired order. They can be incorporated in a finite element or a finite difference scheme; in the present paper the latter is used. In contrast to previous papers using similar formulations, here the method is applied to a fully exterior two‐dimensional problem, with a rectangular boundary. Numerical examples in infinite domains are used to demonstrate the performance and advantages of the new method. In the auxiliary‐variable formulation long‐time corner instability is observed, that requires special treatment of the corners (not addressed in this paper). No such difficulties arise in the high‐derivative formulation. Published in 2005 by 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.     

Local transmitting boundaries for the generalized interpolation material point method

The boundary‐value problems of mechanics can be solved using the material point method with explicit solver formulations. In explicit formulations, even quasi‐static problems are solved as if dynamic, which means that waves are reflected at computational boundaries, generating spurious oscillations in the solution to the boundary‐value problem. Such oscillations can be reduced to a level such that they are barely noticeable with the use of transmitting boundaries. Current implementations of transmitting boundaries in the material point method are limited to the standard viscous boundary. The absence of any stiffness component in the standard viscous boundary may lead to an undesirable finite rigid‐body motion over time. This motion can be minimized through the adoption of the transmitting cone boundary that approximates the stiffness of the unbounded domain. This paper lays out the implementation of the transmitting cone boundary for the generalized interpolation material point method. The cone boundary reflection‐canceling tractions can be applied to either the edges or the centroids of material points; this paper discusses the implications of both approaches.     

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.     

A wideband fast multipole accelerated boundary integral equation method for time‐harmonic elastodynamics in two dimensions

This article presents a wideband fast multipole method (FMM) to accelerate the boundary integral equation method for two‐dimensional elastodynamics in frequency domain. The present wideband FMM is established by coupling the low‐frequency FMM and the high‐frequency FMM that are formulated on the ingenious decomposition of the elastodynamic fundamental solution developed by Nishimura's group. For each of the two FMMs, we estimated the approximation parameters, that is, the expansion order for the low‐frequency FMM and the quadrature order for the high‐frequency FMM according to the requested accuracy, considering the coexistence of the derivatives of the Helmholtz kernels for the longitudinal and transcendental waves in the Burton–Muller type boundary integral equation of interest. In the numerical tests, the error resulting from the fast multipole approximation was monotonically decreased as the requested accuracy level was raised. Also, the computational complexity of the present fast boundary integral equation method agreed with the theory, that is, Nlog N, where N is the number of boundary elements in a series of scattering problems. The present fast boundary integral equation method is promising for simulations of the elastic systems with subwavelength structures. As an example, the wave propagation along a waveguide fabricated in a finite‐size phononic crystal was demonstrated. Copyright © 2012 John Wiley & Sons, Ltd.     

A semi‐analytical curved element for linear elasticity based on the scaled boundary finite element method

This work introduces a semi‐analytical formulation for the simulation and modeling of curved structures based on the scaled boundary finite element method (SBFEM). This approach adapts the fundamental idea of the SBFEM concept to scale a boundary to describe a geometry. Until now, scaling in SBFEM has exclusively been performed along a straight coordinate that enlarges, shrinks, or shifts a given boundary. In this novel approach, scaling is based on a polar or cylindrical coordinate system such that a boundary is shifted along a curved scaling direction. The derived formulations are used to compute the static and dynamic stiffness matrices of homogeneous curved structures. The resulting elements can be coupled to general SBFEM or FEM domains. For elastodynamic problems, computations are performed in the frequency domain. Results of this work are validated using the global matrix method and standard finite element analysis. Copyright © 2016 John Wiley & Sons, Ltd.     

High‐order global–regional model interaction: Extension of Carpenter's scheme

A two‐dimensional global–regional model interaction problem for linear time‐dependent waves is considered. The setup, which is sometimes called ‘one‐way nesting,’ arises in numerical weather prediction as well as in other fields concerning waves in very large domains. It involves the interaction of a coarse global model and a fine limited‐area (regional) model through an ‘open boundary.’ The multiscale nature of this general problem is described. The Carpenter scheme, originally proposed in a note by K. M. Carpenter in 1982 for this type of problem, is then revisited, in the context of the linear scalar wave equation. The original Carpenter scheme is based on the Sommerfeld radiation operator and thus is associated with low‐order accuracy. By replacing the Sommerfeld operator with the high‐order Hagstrom–Warburton absorbing operator, a modified Carpenter open‐boundary condition emerges, which possesses high‐order accuracy. This boundary condition is incorporated in a computational scheme, which uses finite element discretization in space and Newmark time‐stepping. Error analysis and numerical tests for wave guides demonstrate the performance of the modified scheme for combinations of incoming and outgoing waves. Copyright © 2008 John Wiley & Sons, Ltd.