Numerical operational methods for time-dependent linear problems

A general and systematic discussion on the use of the operational method of Laplace transform for numerically solving complex time-dependent linear problems is presented. Application of Laplace transform with respect to time on the governing differential equations as well as the boundary and initial conditions of the problem reduces it to one independent of time, which is solved in the transform domain by any convenient numerical technique, such as the finite element method, the finite difference method or the boundary integral equation method. Finally, the time domain solution is obtained by a numerical inversion of the transformed solution. Eight existing methods of numerical inversion of the Laplace transform are systematically discussed with respect to their use, range of applicability, accuracy and computational efficiency on the basis of some framework vibration problems. Other applications of the Laplace transform method in conjunction with the finite element method or the boundary integral equation method in the areas of earthquake dynamic response of frameworks, thermaliy induced beam vibrations, forced vibrations of cylindrical shells, dynamic stress concentrations around holes in plates and viscoelastic stress analysis are also briefly described to demonstrate the generality and advantages of the method against other known methods.     

Transient Analysis of Elastic Wave Propagation in Multilayered Structures

In this article, explicit transient solutions for one-dimensional wave propagation behavior in multi-layered structures are presented. One of the objectives of this study is to develop an effective analytical method for constructing solutions in multilayered media. Numerical calculations are performed by three methods: the generalized ray method, numerical Laplace inversion method (Durbin's formula), and finite element method (FEM). The analytical result of the generalized ray solution for multilayered structures is composed of a matrix-form Bromwich expansion in the transform domain. Every term represents a group of waves, which are transmitted or reflected through the interface. The matrix representation of the solution can be used to calculate the transient response, without tracing the ray path manually. Numerical inversion of the Laplace transform by Durbin's formula is also used to construct transient responses. This numerical Laplace inversion technique has the advantage of calculating long-time transient responses for complicated multilayered structures. FEM results agree well with calculations obtained by the generalized ray method and numerical Laplace inversion.     

Two-dimensional generalized thermal shock problem of a thick piezoelectric plate of infinite extent

The theory of generalized thermoelasticity, based on the theory of Green and Lindsay with two relaxation times, is used to deal with a thermoelastic–piezoelectric coupled two-dimensional thermal shock problem of a thick piezoelectric plate of infinite extent by means of the hybrid Laplace transform-finite element method. The generalized thermoelastic–piezoelectric coupled finite element equations are formulated. By using Laplace transform the equations are solved and the solutions of the temperature, displacement and electric potential are obtained in the Laplace transform domain. Then the numerical inversion is carried out to obtain the temperature, displacement and electric potential distributions in the physical domain. The distributions are represented graphically. From the distributions, it can be found the wave type heat propagation in the piezoelectric plate. The heat wavefront moves forward with a finite speed in the piezoelectric plate with the passage of time. This indicates that the generalized heat conduction mechanism is completely different from the classic Fourier’s in essence. In generalized thermoelasticity theory heat propagates as a wave with finite velocity instead of infinite velocity in media.     

Application of 3D time domain boundary element formulation to wave propagation in poroelastic solids

The dynamic responses of fluid-saturated semi-infinite porous continua to transient excitations such as seismic waves or ground vibrations are important in the design of soil-structure systems. Biot's theory of porous media governs the wave propagation in a porous elastic solid infiltrated with fluid. The significant difference to an elastic solid is the appearance of the so-called slow compressional wave. The most powerful methodology to tackle wave propagation in a semi-infinite homogeneous poroelastic domain is the boundary element method (BEM). To model the dynamic behavior of a poroelastic material in the time domain, the time domain fundamental solution is needed. Such solution however does not exist in closed form. The recently developed ‘convolution quadrature method’, proposed by Lubich, utilizes the existing Laplace transformed fundamental solution and makes it possible to work in the time domain. Hence, applying this quadrature formula to the time dependent boundary integral equation, a time-stepping procedure is obtained based only on the Laplace domain fundamental solution and a linear multistep method. Finally, two examples show both the accuracy of the proposed time-stepping procedure and the appearance of the slow compressional wave, additionally to the other waves known from elastodynamics.     

冲击波作用下悬臂梁瞬态响应的传递函数方法

利用传递函数方法对爆炸冲击波作用下的悬臂梁进行分析，得到其瞬态响应的封闭形式的解析解。首先，针对给定的爆炸冲击波载荷和初始条件，对梁的控制方程和边界条件进行Laplace变换，然后通过引入状态向量将其改写成状态空间形式，并利用传递函数方法求得其在频域内的解析解，最后利用Crump方法并结合ε算法进行Laplace逆变换，求得悬臂梁在时域内的瞬态响应。给出数值算例，通过与有限元的比较验证了方法的正确性。     

Finite elements for field problems in cylindrical co-ordinates

This paper deals with the extension of the finite element method as applied to the solution of the Laplace or wave equation to cylindrical co-ordinate systems. The base matrices required for solving problems governed by these equations are derived for circular polar, elliptic cylinder and parabolic cylinder co-ordinates. The matrices allow problems whose boundaries are described as co-ordinate surfaces in cylinder co-ordinates to be attacked directly by the finite element method. The subject is discussed from the point of view of one interested in electromagnetic wave propagation in uniform waveguide structures.     

A new efficient numerical method and the dynamic analysis of composite laminates with piezoelectric layers

Firstly, a numerical method for the inversion of Laplace transform is developed and its accuracy is shown through examples. Then, a state-vector equation for the dynamic problems of piezoelectric plates is deduced directly from a modified mixed variational principle for piezoelectric bodies and its exact solution for the dynamic problems of simply supported rectangle piezoelectric plate is simply given. For multilayered hybrid plates, we derive the solution in terms of the propagator matrices. The techniques accounts for the compatibility of generalized displacements and generalized stresses on the interface both the elastic layers and piezoelectric layers, and the transverse shear deformation and the rotary inertia of laminate are also considered in the global algebraic equation of structure. Meanwhile, there is no restriction on the thickness and the number of layers. As an application of the numerical inversion of Laplace transform presented in this paper, typical numerical examples of the harmonic vibration and transient response are proposed and discussed. Since the highly accurate numerical results, they can serve as benchmarks to test various thick plate theories and various numerical methods, such as the finite and boundary element methods for transient response problems.     

A numerical treatment to the solution of quasiparabolic partial differential equations

This paper presents results obtained by the implementation of a hybrid Laplace transform finite element method to the solution of quasiparabolic problem. The present method removes the time derivatives from the quasiparabolic partial differential equation using the Laplace transform and then solves the associated equation with the finite element method. The numerical inverse of the Laplace transform is realized by solving linear overdetermined systems and a polynomial equation of the kth order. Test examples are used to show that the numerical solution is comparable to the exact solution of the initial-boundary value problem at the given grid points.     

DRBEM solution of exterior nonlinear wave problem using FDM and LSM time integrations

The nonlinear wave equation is solved numerically in an exterior region. For the discretization of the space derivatives dual reciprocity boundary element method (DRBEM) is applied using the fundamental solution of Laplace equation. The time derivative and the nonlinearity are treated as the nonhomogenity. The boundary integrals coming from the far boundary are eliminated using rational and exponential interpolation functions which have decay properties far away from the region of interest. The resulting system of ordinary differential equations in time are solved using finite difference method (FDM) with a relaxation parameter and least squares method (LSM). The proposed methods are examined with numerical test problems in which the behaviours of solutions are known. Although it gives almost the same accuracy with the DRBEM+FDM procedure, DRBEM+LSM solution procedure is preferred, since it is a direct method without the need of a parameter.     

Green's functions for the stress intensity factor evolution in finite cracks in orthotropic materials

The elastodynamic response of an infinite orthotropic material with finite crack under concentrated loads is examined. Solution for the stress intensity factor history around the crack tips is found. Laplace and Fourier transforms are employed to solve the equations of motion leading to a Fredholm integral equation on the Laplace transform domain. The dynamic stress intensity factor history can be computed by numerical Laplace transform inversion of the solution of the Fredholm equation. Numerical values of the dynamic stress intensity factor history for some example materials are obtained. This solution can be used as a Green's function to solve dynamic problems involving fini te cracks.     

瞬态弹性动力学问题的半解析边界元法

本文用半解析有限元法对边界积分方程作离散化处理,通过引入基本解函数和半解析半离散试函数的二次半解析过程,使三维弹性动力学问题简化为一维数值计算。文中又采用移动边界元法来模拟波在半无限介质中传播的表面积分问题,分析计算了各种瞬态波在介质内传播,绕射及地面运动问题。计算结果表明,半解析边界元法不仅计算精度高,而且工作量大大降低,具有较高的经济效益与应用价值。     

Torsional impact response of a penny-shaped crack lying on a bimaterial interface

The torsional impact response of a penny-shaped crack lying on a bimaterial interface is considered in this study. Laplace and Hankel transforms are used to reduce the problem to the solution of a pair of dual integral equations. The solution to the dual integral equations is expressed in terms of a Fredholm integral equation of the second kind with a finite integral kernel. A numerical Laplace inversion routine is used to recover the time dependence of the solution. The dynamic stress intensity factor is determined and its dependence on time and material constants is discussed.     

Boundary integral equation method for linear porous-elasticity with applications to soil consolidation

For physical phenomena governed by the Biot model of porous-elasticity, a reciprocal relation, similar to the Betti's recoprocal theorem in elasticity, is constructed in Laplace transformed space. Integrating the reciprocal relation enables one to formulate boundary integral equations. The fundamental kernels for the integral equations are solved in closed forms for the case of isotropic material. Numerical implementation of two-dimensional problems includes finite element ideas of discretization and polynomial interpolation, and numerical inversion of a Laplace transform. Practical applications of the method are found in consolidation problems in soils which contain compressible as well as incompressible pore fluids. Also, as a numerical experiment, consolidation of partially saturated soil is simulated and interesting phenomena are observed. The currently developed boundary integral equation method (BIEM) for porous-elasticity may be viewed as an efficient and accurate alternative of existing finite element and finite difference methods. For linear consolidation problems, application of BIEM is always preferred to the other numerical methods whenever possible.     

CQM-based BEM formulation for uncoupled transient quasistatic thermoelasticity analysis

A boundary element method based on the convolution quadrature method for the numerical solution of uncoupled transient thermoelasticity problems is presented. In the proposed formulation, the time-domain integral equation is numerically approximated by a quadrature formula whose weight factors are computed by means of integral expressions involving the Laplace transform of the fundamental solution. Compared with other numerical methods that operate directly in the Laplace transformed domain, the proposed formulation requires only the definition of the time-step used by the procedure of integration, and does not need special techniques of inversion from the Laplace-domain to the time-domain. Numerical examples of transient thermoelasticity problems are presented to show the versatility and accuracy of the method.     

A hybrid discontinuous in space and time Galerkin method for wave propagation problems

Medium‐frequency regime and multi‐scale wave propagation problems have been a subject of active research in computational acoustics recently. New techniques have attempted to overcome the limitations of existing discretization methods that tend to suffer from dispersion. One such technique, the discontinuous enrichment method, incorporates features of the governing partial differential equation in the approximation, in particular, the solutions of the homogeneous form of the equation. Here, based on this concept and by extension of a conventional space–time finite element method, a hybrid discontinuous Galerkin method (DGM) for the numerical solution of transient problems governed by the wave equation in two and three spatial dimensions is described. The discontinuous formulation in both space and time enables the use of solutions to the homogeneous wave equation in the approximation. In this contribution, within each finite element, the solutions in the form of polynomial waves are employed. The continuity of these polynomial waves is weakly enforced through suitably chosen Lagrange multipliers. Results for two‐dimensional and three‐dimensional problems, in both low‐frequency and medium‐frequency regimes, show that the proposed DGM outperforms the conventional space–time finite element method. Copyright © 2014 John Wiley & Sons, Ltd.     

Initial conditions contribution in a BEM formulation based on the convolution quadrature method

This work presents a two‐dimensional boundary element method (BEM) formulation for the analysis of scalar wave propagation problems. The formulation is based on the so‐called convolution quadrature method (CQM) by means of which the convolution integral, presented in time‐domain BEM formulations, is numerically substituted by a quadrature formula, whose weights are computed using the Laplace transform of the fundamental solution and a linear multistep method. This BEM formulation was initially developed for scalar wave propagation problems with null initial conditions. In order to overcome this limitation, this work presents a general procedure that enables one to take into account non‐homogeneous initial conditions, after replacing the initial conditions by equivalent pseudo‐forces. The numerical results included in this work show the accuracy of the proposed BEM formulation and its applicability to such kind of analysis. Copyright © 2006 John Wiley & Sons, Ltd.     

Thermoelastic analysis of a layered slab subjected to a moving line heat source

Abstract

This paper reports a theoretical study of the transient thermal stress distributions of a layered slab which is composed of two different materials. The layered slab is heated by a moving line heat source on its upper surface and cooled convectively on the lower surface. In order to solve the initial and boundary value problem, a general hybrid Laplace transform/finite element method is utilized. Finally, a numerical procedure, the Fourier series technique, is used to obtain the inversion of the Laplace transform. The effect of the number of mesh elements in the X‐direction is also investigated to verify the accuracy and convergence of the finite element method. In addition, a typical result is compared with the analytic solution. The numerical results of the transient temperature and thermal stress distribution of the layered slab are presented to demonstrate the effect of the physical properties.     

Boundary type finite element method for surface wave motion based on trigonometric function interpolation

There are many physical phenomena which can be handled by the Helmholtz equation. The equation explains certain phenomena of wave propagation. This paper presents a new finite element method to analyse surface wave motion. The characteristic point of this method is that the interpolation equation is chosen to satisfy the governing Helmholtz equation using trigonometric functions. This follows that the variational functional to be minimized can be formulated such that the integration is limited to the boundary of the element. The numerical solutions obtained are compared with analytical and experimental solutions. From these comparative studies, it is concluded that the present method provides a useful tool for the analysis of surface wave motion.     

A frequency-domain approach for modelling transient elastodynamics using scaled boundary finite element method

This study develops a frequency-domain method for modelling general transient linear-elastic dynamic problems using the semi-analytical scaled boundary finite element method (SBFEM). This approach first uses the newly-developed analytical Frobenius solution to the governing equilibrium equation system in the frequency domain to calculate complex frequency-response functions (CFRFs). This is followed by a fast Fourier transform (FFT) of the transient load and a subsequent inverse FFT of the CFRFs to obtain time histories of structural responses. A set of wave propagation and structural dynamics problems, subjected to various load forms such as Heaviside step load, triangular blast load and ramped wind load, are modelled using the new approach. Due to the semi-analytical nature of the SBFEM, each problem is successfully modelled using a very small number of degrees of freedom. The numerical results agree very well with the analytical solutions and the results from detailed finite element analyses.     

纤维增强塑料齿轮轮齿弯曲强度分析

本文由纤维增强塑料齿轮轮齿的细观结构建立了各向异性粘弹性力学的计算模型,采用Laplace变换和有限单元法相结合求解轮齿的粘弹性应力场。通过引入时间折算因子,改进了Schapery的数值反演过程,给出了一个可行的复合材料齿轮弯曲强度实用分析方法。