The mixed-type integral equations for transient elastodynamic plane crack problems

In the present paper, by use of the boundary integral equation method and the techniques of Green fundamental solution and singularity analysis, the dynamic infinite plane crack problem is investigated. For the first time, the problem is reduced to solving a system of mixed-typed integral equations in Laplace transform domain. The equations consist of ordinary boundary integral equations along the outer boundary and Cauchy singular integral equations along the crack line. The equations obtained are strictly proved to be equivalent with the dual integral equations obtained by Sih in the special case of dynamic Griffith crack problem. The mixed-type integral equations can be solved by combining the numerical method of singular integral equation with the ordinary boundary element method. Further use the numerical method for Laplace transform, several typical examples are calculated and their dynamic stress intensity factors are obtained. The results show that the method proposed is successful and can be used to solve more complicated problems.     

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.     

Dynamic stress concentration studies by boundary integrals and Laplace transform

The dynamic stress field and its concentrations around holes of arbitrary shape in infinitely extended bodies under plane stress or plane strain conditions are numerically determined. The material may be linear elastic or viscoelastic, while the dynamic load consists of plane compressional waves of harmonic or general transient nature. The method consists of applying the Laplace transform with respect to time to the governing equations of motion and formulating and solving the problem numerically in the transfomed domain by the boundary integral equation method. The stress field can then be obtaind by a numerical inversion of the trasformed solution. The correspondence principle is invoked for the case of viscoelastic material behavious. The method is simplified for the case of harmonic waves where no numerical inversion is involved.     

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.     

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.     

An indirect boundary element method for three-dimensional dynamic fracture mechanics

Indirect boundary element methods (fictitious load and displacement discontinuity) have been developed for the analysis of three-dimensional elastostatic and elastodynamic fracture mechanics problems. A set of boundary integral equations for fictitious loads and displacement discontinuities have been derived. The stress intensity factors were obtained by the stress equivalent method for static loading. For dynamic loading the problem was studied in Laplace transform space where the numerical calculation procedure, for the stress intensity factor K_I(p), is the same: as that for the static problem. The Durbin inversion method for Laplace transforms was used to obtain the stress intensity factors in the time domain K_I(t). Results of this analysis are presented for a square bar, with either a rectangular or a circular crack, under static and dynamic loads.     

A finite element Laplace transform solution technique for the wave equation

A technique is described for the solution of the wave equation with time dependent boundary conditions. The finite element solution accompanied by the numerical Laplace inversion process seems to be an efficient procedure to treat such problems. The programming involved is straightforward in the sense that numerical Laplace inversion routines can be directly used as a time integration procedure after obtaining standard finite element differential equation solutions in the transformed domain. Some results are presented for one- and two- dimensional applications, such as wave propagation in longitudinal bars and wave propagation in harbours.     

Three-dimensional dynamic fracture analysis with the dual reciprocity method in Laplace domain

In this paper, the dual reciprocity boundary element method in the Laplace domain has been developed for the analysis of three-dimensional elastodynamic fracture mechanics mixed-mode problems. The boundary element method is used to calculate the unknowns of transformed boundary displacement and traction and the domain integrals in the elastodynamic equation are transformed into boundary integrals by the use of the dual reciprocity method. The transformed dynamic stress intensity factors are determined by the crack opening displacement (COD) directly in the Laplace domain. By using Durbin's inversion technique, the dynamic stress intensity factors in the time domain are obtained. Several numerical examples are presented to demonstrate the good agreement with existing solutions.     

Boundary element analysis of stress intensity factor KI in some two-dimensional dynamic thermoelastic problems

This paper presents a numerical technique for the calculation of stress intensity factor as a function of time for coupled thermoelastic problems. In this task, effect of inertia term considering coupled theory of thermoelasticity is investigated and its importance is shown.A boundary element method using Laplace transform in time-domain is developed for the analysis of fracture mechanic considering dynamic coupled thermoelasticity problems in two-dimensional finite domain. The Laplace transform method is applied to the time-domain and the resulting equations in the transformed field are discretized using boundary element method. Actual physical quantities in time-domain is obtained, using the numerical inversion of the Laplace transform method.The singular behavior of the temperature and stress fields in the vicinity of the crack tip is modeled by quarter-point elements. Thermal dynamic stress intensity factor for mode I is evaluated using J-integral method. By using J-integral method effects of inertia term and other terms such as strain energy on stress intensity factor may be investigated separately and their importance may be shown. The accuracy of the method is investigated through comparison of the results with the available data in literature.     

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.     

Dynamic stress intensity factors studied by boundary integro-differential equations

The boundary integro-differential equation method is illustrated by two numerical examples concerning the study of the dynamic stress intensity factor around a penny-shaped crack in an infinite elastic body. Harmonic and impact load on the crack surface has been considered. Applying the Laplace transform with respect to time to the governing equations of motion the problem is solved in the transformed domain by the boundary integro-differential equations. The Laplace transformed general transient problem can be used to solve the steady-state problem as a special case where no numerical inversion is involved.     

Infinite elements in the time domain using a prolate spheroidal multipole expansion

Many phenomena in acoustically loaded structural vibrations are better understood in the time domain, particularly transient radiation, shock, and problems involving non-linearities, cavitation, and bulk structural motion. In addition, the geometric complexity of structures of interest drives the analyst toward domain-discretized solution methods, such as finite elements or finite differences, and large numbers of degrees of freedom. In such methods, efficient numerical enforcement of the Sommerfeld radiation condition in the time domain becomes difficult. Although a great many methodologies for doing so have been demonstrated, there seems to exist no consensus on the optimal numerical implementation of this boundary condition in the time domain. Here, we present theoretical development of several new boundary operators for conventional finite element codes. Each proceeds from successful domain-discretized, projected field-type harmonic solutions, in contrast to boundary integral equation operators or those derived from analyses of outgoing waves. We exploit the separable prolate-spheroidal co-ordinate system, which is sufficiently general for a large variety of problems of naval interest, to obtain finite element-like operators (matrices) for the boundary points. Use of this co-ordinate system results in element matrices that can be analytically inverse transformed from the frequency to the time domain, without imposing continuity requirements on the solution above those imposed by the underlying partial differential equation. In addition, use of element-like boundary operators does not alter the banded structure of assembled system matrices. Results presented here include theoretical derivation of the infinite elements, resolution of the Fourier inversion issues, and element matrices for the boundary operators which introduce no new continuity requirements on the fluid field variable. The simplest infinite elements are verified in a coupled three-dimensional context against DAA2 and Helmholtz integral equation results. © 1998 John Wiley & Sons, Ltd.     

Free and forced vibrations of plates by boundary and interior elements

A direct boundary element method is developed for the dynamic analysis of thin elastic flexural plates of arbitrary planform and boundary conditions. The formulation employs the static fundamental solution of the problem and this creates not only boundary integrals but surface integrals as well owing to the presence of the inertia force. Thus the discretization consists of boundary as well as interior elements. Quadratic isoparametric elements and quadratic isoparametric or constant elements are employed for the boundary and interior discretization, respectively. Both free and forced vibrations are considered. The free vibration problem is reduced to a matrix eigenvalue problem with matrix coefficients independent of frequency. The forced vibration problem is solved with the aid of the Laplace transform with respect to time and this requires a numerical inversion of the transformed solution to obtain the plate dynamic response to arbitrary transient loading. The effect of external viscous or internal viscoelastic damping on the response is also studied. The proposed method is compared against the direct boundary element method in conjunction with the dynamic fundamental solution as well as the finite element method primarily by means of a number of numerical examples. These examples also serve to illustrate the use of the proposed method.     

Dynamic Stress Intensity Factor of Interfacial Finite Cracks in Orthotropic Materials

The elastodynamic response of an infinite non-homogeneous orthotropic material with an interfacial finite crack under distributed normal and shear impact 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 materials are obtained. Interfacial cracks between two different materials and between two pieces of the same material but different fiber orientation are considered. Bimaterial formulation of a crack problem is shown to converge to the mono-material formulation, derived independently, in the limiting case when both materials are the same.     

Boundary element analysis of dynamic coupled thermoelasticity problems

Some possibility of numerical analysis of coupled dynamic problems of linear elastic heat conductors on classical thermoelasticity theory by using the boundary element method is shown in this paper. The boundary integral equation formulation and its numerical implementation of the two-dimensional problem are developed in the manner by the newly derived fundamental solution for the coupled equations of elliptic type in the transformed space and the numerical inversion of Laplace transformation. The boundary element unsteady solutions of the first and second Danilovskaya problems and the Sternberg and Chakravorty problem in the half-space are demonstrated through comparison with the existing solutions.     

Dynamic stress intensity factor (Mode I) of a penny-shaped crack in an infinite poroelastic solid

This paper considers the transient stress intensity factor (Mode I) of a penny-shaped crack in an infinite poroelastic solid. The crack surfaces are impermeable. By virtue of the integral transform methods, the poroelastodynamic mixed boundary value problems is formulated as a set of dual integral equations, which, in turn, are reduced to a Fredholm integral equation of the second kind in the Laplace transform domain. Time domain solutions are obtained by inverting Laplace domain solutions using a numerical scheme. A parametric study is presented to illustrate the influence of poroelastic material parameters on the transient stress intensity. The results obtained reveal that the dynamic stress intensity factor of poroelastic medium is smaller than that of elastic medium and the poroelastic medium with a small value of the potential of diffusivity shows higher value of the dynamic stress intensity factor.     

A finite element analysis for time-dependent problems

The finite element method is applied to dynamic linear viscoelastic analysis. The arising system of differential equations is integrated by means of the Laplace transformation. The direct numerical evaluation of the inversion integral seems to be appropriate in case of vibrational problems. The method is applied to examples of transient vibrations of elastic bodies. Initial conditions and arbitrary damping matrices can be considered.     

Convolution quadrature method‐based symmetric Galerkin boundary element method for 3‐d elastodynamics

The boundary integral equations in 3‐d elastodynamics contain convolution integrals with respect to the time. They can be performed analytically or with the convolution quadrature method. The latter time‐stepping procedure's benefit is the usage of the Laplace‐transformed fundamental solution. Therefore, it is possible to apply this method also to problems where analytical time‐dependent fundamental solutions might not be known. To obtain a symmetric formulation, the second boundary integral equation has to be used which, unfortunately, requires special care in the numerical implementation since it involves hypersingular kernel functions. Therefore, a regularization for closed surfaces of the Laplace‐transformed elastodynamic kernel functions is presented which transforms the bilinear form of the hypersingular integral operator to a weakly singular one. Supplementarily, a weakly singular formulation of the Laplace‐transformed elastodynamic double layer potential is presented. This results in a time domain boundary element formulation involving at least only weakly singular integral kernels. Finally, numerical studies validate this approach with respect to different spatial and time discretizations. Further, a comparison with the wider used collocation method is presented. It is shown numerically that the presented formulation exhibits a good convergence rate and a more stable behavior compared with collocation methods. Copyright © 2008 John Wiley & Sons, Ltd.     

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.     

Impact response of a strip composite with a finite crack

The dynamic response of a central crack in a strip composite under normal impact is analyzed. The crack is oriented normally to the interfaces. Laplace and Fourier transform techniques are used to reduce the elastodynamic problem to a pair of dual integral equations. The integral equations are solved by using an integral transform technique and the result is expressed in terms of a Fredholm integral equation of the second kind. 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, the material properties and the geometrical parameters is discussed.