Transient elastic waves propagating in a multi-layered medium subjected to in-plane dynamic loadings. I. Theory (original) (raw)
Related papers
Dynamic responses of a layered medium subjected to anti-plane loadings
International Journal of Solids and Structures, 2001
The dynamic response of a layered medium subjected to anti-plane loadings is investigated in this study. An eective analytical method is used to present transient full ®eld solutions in the layered medium. The boundary value problem is solved using Laplace transform with respect to time and Fourier transform with respect to space. The transient response is then implemented analytically by Cagniard's method that obviates the contour integration diculties and whose terms represent successively re¯ected and transmitted waves. The static solution is obtained by the extreme case of transient response while the time is suciently long. The objective in this study is to investigate the transition phenomena from transient response to static value. The numerical examples of a layered half-space are illustrated and discussed in detail. It is found that the characteristic time of the transition not only depends on the distance between source and receiver but also is strongly related to the wave velocities of layer and half-space. Furthermore, very interesting transition behavior is observed when the layer has faster wave velocity than the half-space.
Wave propagation in inhomogeneous layered media: solution of forward and inverse problems
Acta Mechanica, 2004
Wave propagation in anisotropic inhomogeneous layered media due to high frequency impact loading is studied using a new Spectral Layer Element (SLE). The element can model functionally graded materials (FGM), where the material property variation is assumed to follow an exponential function. The element is exact for a single parameter model which describes both moduli and density variation. This novel element is formulated using the method of partial wave technique (PWT) in conjunction with linear algebraic methodology. The matrix structure of finite element (FE) formulation is retained, which substantially simplifies the modeling of a multi-layered structure. The developed SLE has an exact dynamic stiffness matrix as it uses the exact solution of the governing elastodynamic equation in the frequency domain as its interpolation function. The mass distribution is modeled exactly, and, as a result, the element gives the exact frequency response of each layer. Hence, one element may be as large as one complete layer which results in a system size being very small compared to conventional FE systems. The Fast-Fourier Transform (FFT) and Fourier series are used for the inversion to the time/space domain. The formulated element is further used to study the stress distribution in multi-layered media. As a natural application, Lamb wave propagation in an inhomogeneous plate is studied and the time domain description is obtained. Further, the advantage of the spectral formulation in the solution of inverse problems, namely the force identification and system identification is investigated. Constrained nonlinear optimization technique is used for the material property identification, whereas the transfer function approach is taken for the impact force identification.
Journal of Applied Mechanics, 2001
The paper develops a universal method for SH-wave propagation in a multilayered medium with an arbitrary number of interface cracks. The method makes use of the transfer matrix and Fourier integral transform techniques to cast the mixed boundary value problem to a set of Cauchy singular integral equations of the first type which can be solved numerically. The paper calculates the dynamic stress intensity factors for some simple but typical examples.
Transient Analysis of Crack in Composite Layered Medium Subjected to Dynamic Loadings
AIAA Journal, 1997
A semi-in® nite crack in a layered medium subjected to antiplane dynamic loading is investigated. In analyzing this problem, the fact that the re¯ections and diffractions of stress waves by the interface boundary and by the crack will generate an in® nite number of waves must be taken into account. A useful fundamental solution is proposed, and the solution is determined by superposition of the fundamental solution in the Laplace transform domain. The proposed fundamental problem is the application of exponentially distributed traction (in the Laplace transform domain) on the crack faces. Cagniard's method for Laplace inversion is used to obtain the transient solution in time domain. The ® nal formulations for the stress intensity factor are expressed explicitly and the dynamic effect of each wave is presented in a closed form. The results are valid for an in® nite number of waves that are scattered from the crack tip and re¯ected by the interface. Numerical calculations of dynamic stress intensity factors are evaluated and discussed in detail.
Two-dimensional transient dynamic response of orthotropic layered media
Meccanica, 2006
In this study, two-dimensional transient dynamic response of orthotropic plane layered media is investigated. The plane multilayered media consist of N different generally orthotropic, homogeneous and linearly elastic layers with different ply angles. In the generally orthotropic layer, representing a ply reinforced by unidirectional fibers with an arbitrary orientation angle, the principal material directions do not coincide with body coordinate axes. The solution is obtained by employing a numerical technique which combines the use of Fourier transform with the method of characteristics. The numerical results are displayed in curves denoting the variations of stress and displacement components with time at different locations. These curves clearly reveal, in wave profiles, the scattering effects caused by the reflections and refractions of waves at the boundaries and at the interfaces of the layers, and also the effects of anisotropy caused by fiber orientation angle. The Preliminary version presented at curves properly predict the sharp variations in the response at the neighborhood of the wave fronts, which shows the power of the numerical technique employed in the study. By suitably adjusting the elastic constants, the results for multilayered media with transversely isotropic layers, or layers with cubic symmetry, or isotropic layers can easily be obtained from the general formulation. Furthermore, solutions for some special cases, including Lamb's problem for an elastic half-space, are obtained and compared with the available solutions in the literature and very good agreement is found.
Transient analysis of three-dimensional wave propagation using the boundary element method
International Journal for Numerical Methods in Engineering, 1987
In the past the time domain solution of the wave equation has been limited to simplified problems. This was due to the limitations of analytical methods and the capacity of computers to manipulate and store 'large' blocks of spatial information. With the advent of 'super computers' the ability to solve such problems has significantly increased. This paper outlines a method for transient analysis of wave propagation in arbitrary domains using a boundary element method. The technique presented will allow the definition ofa domain, the input of impedance conditions on the domain's surface, the specification of inputs on the surface, and the specification of initial conditions within the domain. It will produce a complete solution of the wave equation inside the domain. The techniques are demonstrated using a program with a boundary element formulation of Kirchhoff s equation. The elements used are triangular and compatible. TNTRODUCTTON Numerical solution of time domain wave propagation studies is possible using a variety of techniques, including the generalized finite difference and finite element methods. However, there are certain advantages to using boundary integral methods since the complete spatial domain of interest need not be discretized. Furthermore, since a Green's function solution is inherent in the boundary integral formulation, some of the interpolation errors which are characteristic of both finite difference and finite element methods are avoided. The original boundary integral time domain wave propagation formulations were developed by Cruse and Rizzo' and Cruse2 to study stress wave propagation in arbitrary domains. The technique solved the integral equations in the Laplace domain and then inverse transformed the solution to obtain the time domain solution. Reference 3 is a more recent work by Rizzo et al., on the three-dimensional treatment of elastic waves in frequency domain. Manolis and Beskos4 numerically solved the formulation of Cruse and Rizzo for two-dimensional dynamic stress concentration studies. This transformation technique requircs a large number of boundary integral evaluations in the Laplace domain in order that the time domain solution be accurate and have a large bandwidth. Thus, the Laplace transform-based techniques are relatively expensive, especially for three-dimensional problems. Direct time domain solutions of scalar wave propagation have been developed by Mansur and
Comparison of the Matrix Techniques for Modeling Stress Waves in Multilayered Media
2015
This paper describes the development of a general purpose model for wave prop-agation in Cartesian system. When isotropic materials are used, the model can account for elastic and visco-elastic isotropic materials, single or multi-layered struc-tures, and free or leaky systems. Model determines what resonances can exist in order to satisfy the boundary conditions and the bulk wave propagation character-istics in each of the layers. These resonances control how ultrasonic waves will be guided in the system and what properties each of these waves will have. Key words: guided waves, dispersive curves 1