Pencil method in elastodynamics: application to ultrasonic field computation

The principles of pencil elastodynamics and, in more detail, some selected applications of pencil techniques to elastodynamics are described. It is shown how a systematic use of a matrix representation for the wave front curvature and for its transformations simplifies the handling of arbitrary pencils and, consequently, the field computations. Pencil matrix representations for the propagation into homogeneous solids made of isotropic or anisotropic media are derived. The use of matrix representations for pencil reflections on, or refractions through, arbitrarily curved interfaces, together with matrix representations for propagation into homogeneous media, allow us to derive an overall matrix formulation for elastodynamic propagation into complex heterogeneous structures. Combined with the classical Rayleigh integral to account for transducer diffraction effects, the proposed theory is applied to the prediction of ultrasonic fields radiated into complex structures by arbitrary transducers. Examples of interest for application to ultrasonic non-destructive testing are given.     

Nonlinear surface acoustic waves: realization of solitary pulses and fracture

A laser-based technique for the contact-free generation and detection of strongly nonlinear surface acoustic wave (SAW) pulses with amplitudes limited by the materials strength has been developed. The effects of nonlinear propagation of short elastic surface pulses with finite strength in isotropic solids, such as fused quartz, anisotropic solids, such as silicon, and dispersive media were investigated. Solitary surface wave propagation was observed in layered structures for normal and anomalous dispersion. In addition, a SAW-based method for evaluating the critical fracture stress of anisotropic brittle solids, such as single crystal silicon, is introduced.     

Computing the far field scattered or radiated by objects inside layered fluid media using approximate Green's functions

A numerically efficient technique is presented for computing the field radiated or scattered from three-dimensional objects embedded within layered acoustic media. The distance between the receivers and the object of interest is supposed to be large compared to the acoustic wavelength. The method requires the pressure and normal particle displacement on the surface of the object or on an arbitrary circumscribing surface, as an input, together with a knowledge of the layered medium Green's functions. The numerical integration of the full wave number spectral representation of the Green's functions is avoided by employing approximate formulas which are available in terms of elementary functions. The pressure and normal particle displacement on the surface of the object of interest, on the other hand, may be known by analytical or numerical means or from experiments. No restrictions are placed on the location of the object, which may lie above, below, or across the interface between the fluid media. The proposed technique is verified through numerical examples, for which the near field pressure and the particle displacement are computed via a finite-element method. The results are compared to validated reference models, which are based on the full wave number spectral integral Green's function.     

Analytical close-form solutions to the elastic fields of solids with dislocations and surface stress

The concept of eigenstrain is adopted to derive a general analytical framework to solve the elastic field for 3D anisotropic solids with general defects by considering the surface stress. The formulation shows the elastic constants and geometrical features of the surface play an important role in determining the elastic fields of the solid. As an application, the analytical close-form solutions to the stress fields of an infinite isotropic circular nanowire are obtained. The stress fields are compared with the classical solutions and those of complex variable method. The stress fields from this work demonstrate the impact from the surface stress when the size of the nanowire shrinks but becomes negligible in macroscopic scale. Compared with the power series solutions of complex variable method, the analytical solutions in this work provide a better platform and they are more flexible in various applications. More importantly, the proposed analytical framework profoundly improves the studies of general 3D anisotropic materials with surface effects.     

Analysis of stress field in a single anisotropic diamond/square shaped inclusion using the volume integral equation method and the numerical equivalent inclusion method

A volume integral equation method (VIEM) is used to study elastostatic problems in an unbounded elastic solid containing a single diamond/square shaped inclusion subject to uniform tensile stress at infinity. The inclusion is assumed to be a long parallel diamond/square cylinder composed of isotropic or anisotropic elastic materials and perfectly bonded to the isotropic matrix. The solid is assumed to be under plane strain on the plane normal to the cylinder. A detailed analysis of the stress field at the interface between the isotropic matrix and the single isotropic/orthotropic diamond/square shaped inclusion is carried out. The effects of a single isotropic/orthotropic diamond/square shaped inclusion on the stress field at the interface between the matrix and the inclusion are investigated in detail. The accuracy of the volume integral equation method for the interfacial stress field is validated and compared by the numerical equivalent inclusion method (NEIM) and the finite element method (FEM) using ADINA. Through detailed analysis of plane elastostatic problems using the parallel volume integral equation method (PVIEM) in an unbounded isotropic matrix with multiple isotropic diamond shaped inclusions under uniform remote tensile loading, it is demonstrated that the volume integral equation method can also be applied to solve general two- and three-dimensional elastostatic problems involving multiple isotropic/anisotropic inclusions whose shape and number are arbitrary.     

The computation of dispersion relations for axisymmetric waveguides using the Scaled Boundary Finite Element Method

This paper addresses the computation of dispersion curves and mode shapes of elastic guided waves in axisymmetric waveguides. The approach is based on a Scaled Boundary Finite Element formulation, that has previously been presented for plate structures and general three-dimensional waveguides with complex cross-section. The formulation leads to a Hamiltonian eigenvalue problem for the computation of wavenumbers and displacement amplitudes, that can be solved very efficiently. In the axisymmetric representation, only the radial direction in a cylindrical coordinate system has to be discretized, while the circumferential direction as well as the direction of propagation are described analytically. It is demonstrated, how the computational costs can drastically be reduced by employing spectral elements of extremely high order. Additionally, an alternative formulation is presented, that leads to real coefficient matrices. It is discussed, how these two approaches affect the computational efficiency, depending on the elasticity matrix. In the case of solid cylinders, the singularity of the governing equations that occurs in the center of the cross-section is avoided by changing the quadrature scheme. Numerical examples show the applicability of the approach to homogeneous as well as layered structures with isotropic or anisotropic material behavior.     

Finite element mapping for spring network representations of the mechanics of solids

We present a general finite element mapping procedure for defining spring network representations of solid mechanics. The procedure is rigorous and equally suitable for setting regular and unstructured spring network models of generally anisotropic solids. We use the procedure to define close-packed triangular and simple cubic lattice spring models of isotropic 2D and 3D elastic media, respectively. We extend the study to heterogeneous solids and show that the mapped spring network approach constitutes an appealing route for incorporating subelement level constitutive equations.     

Free vibration of antisymmetric,angle-ply laminated plates including transverse shear deformation by the finite element method

A finite element formulation of the equations governing the laminated anisotropic plate theory of Yang, Norris and Stavsky, is presented. The theory is a generalization of Mindlin's theory for isotropic plates to laminated anisotropic plates and includes shear deformation and rotary inertia effects. Finite element solutions are presented for rectangular plates of antisymmetric angle-ply laminates whose material properties are typical of a highly anisotropic composite material. Two sets of material properties that are typical of high modulus fiber-reinforced composites are used to show the parametric effects of plate aspect ratio, length-to-thickness ratio, number of layers and lamination angle. The numerical results are compared with the closed form results of Bert and Chen. As a special case, numerical results are presented for thick isotropic plates, and are compared with those for 3-D linear elasticity theory and Mindlin's thick plate theory.     

DPSM technique for ultrasonic field modelling near fluid-solid interface

Distributed point source method (DPSM) is gradually gaining popularity in the field of non-destructive evaluation (NDE). DPSM is a semi-analytical technique that can be used to calculate the ultrasonic fields produced by transducers of finite dimension placed in homogeneous or non-homogeneous media. This technique has been already used to model ultrasonic fields in homogeneous and multi-layered fluid structures. In this paper the method is extended to model the ultrasonic fields generated in both fluid and solid media near a fluid-solid interface when the transducer is placed in the fluid half-space near the interface. Most results in this paper are generated by the newly developed DPSM technique that requires matrix inversion. This technique is identified as the matrix inversion based DPSM technique. Some of these results are compared with the results produced by the Rayleigh-Sommerfield integral based DPSM technique. Theory behind both matrix inversion based and Rayleigh-Sommerfield integral based DPSM techniques is presented in this paper. The matrix inversion based DPSM technique is found to be very efficient for computing the ultrasonic field in non-homogeneous materials. One objective of this study is to model ultrasonic fields in both solids and fluids generated by the leaky Rayleigh wave when finite size transducers are inclined at Rayleigh critical angles. This phenomenon has been correctly modelled by the technique. It should be mentioned here that techniques based on paraxial assumptions fail to model the critical reflection phenomenon. Other advantages of the DPSM technique compared to the currently available techniques for transducer radiation modelling are discussed in the paper under Introduction.     

Deviatoric stress: a nuisance or a gold mine?

Both synchrotron radiation and deviatoric stress were once considered to be nuisances. Now synchrotron radiation is one of the most important tools available to scientists of all disciplines and deviatoric stress is one of the most useful aspects of x-ray diffraction at extreme conditions. Samples in high-pressure devices are under true hydrostatic pressure only when surrounded by a fluid, thus limiting true hydrostatic pressure studies at ambient temperatures to pressures below about 11?GPa. Elevated temperature is able to extend this limit but has rarely been used for this purpose. Instead, noble gases have been used as pressure media as their solids are especially soft. Deviatoric stress and resultant anisotropic elastic strain in solid samples and solid media have led to many subtle errors in determinations of elastic properties and crystal structures, especially in the days before it was realized that they could be measured and were potentially a valuable source of information. In recent years, measuring anisotropic elastic strain by x-ray diffraction has provided new insights into materials strength, elastic properties, crystal structures, mechanisms of phase transitions, slip systems, lattice preferred orientation, and, of course, ways to make corrections when deviatoric stress is indeed a nuisance.     

Air-coupled measurement of plane wave, ultrasonic plate transmission for characterising anisotropic, viscoelastic materials

Electrostatic, air-coupled, ultrasonic transducers are used to generate and detect plane waves in viscoelastic, isotropic or anisotropic solid plates. The through-transmitted field is measured and compared to numerical predictions. An inversion scheme is then applied for identifying the values of the complex Cij which are representative of the viscoelasticity properties of the materials. The issue of this work is a contact-free, ultrasonic technique for material characterisation.     

Dynamic equations for a fully anisotropic elastic plate

A hierarchy of dynamic plate equations is derived for a fully anisotropic elastic plate. Using power series expansions in the thickness coordinate for the displacement components, recursion relations are obtained among the expansion functions. Adopting these in the boundary conditions on the plate surfaces and along the edges, a set of dynamic equations with pertinent edge boundary conditions are derived on implicit form. These can be truncated to any order and are believed to be asymptotically correct. For the special case of an orthotropic plate, explicit plate equations are presented and compared analytically and numerically to other approximate theories given in the literature. These results show that the present theory capture the plate behavior accurately concerning dispersion curves, eigenfrequencies as well as stress and displacement distributions.     

Asymmetric and axisymmetric vibrations of finite transversely isotropic circular cylinders

This paper presents an approach for obtaining the exact frequency equations of axisymmetric and asymmetric free vibrations of transversely isotropic circular cylinders. The solution method is based on the three dimensional theory of linear elasticity and uses potential functions. Using this approach, the frequency spectra and vibration mode shapes are plotted for a number of transversely isotropic cylinders. The proposed approach introduces a number of merits compared to earlier approximate and exact solution methods. First, unlike numerically complicated series methods that provide approximate solutions, the proposed approach is exact. Second, combination of scalar functions employed for representing the displacement field is consistent with the physics of the problem. One scalar potential function has been considered for each component of the wave field inside the elastic cylinder. As a result, the solution is systematically divided into coupled and decoupled equations. In addition, by using this approach, there is no need to guess the final of the solution a priori. These merits make the proposed approach suitable for other vibration problems of anisotropic materials.     

Three-dimensional hybrid model for predicting air-coupled generation of guided waves in composite material plates

For contact-less, non-destructive testing (NDT) purposes using air-coupled ultrasonic transducers, it is often required to numerically simulate the propagation of ultrasonic waves in solid media, and their coupling through air with specific transducers. At that point, one could simulate the propagation in the air and then in the solid component, using a Finite Element (FE) model. However, when three-dimensional (3D) modeling becomes necessary, such a solution reveals to be extremely demanding in terms of number of degrees of freedom and computational time. In this paper, to avoid such difficulties, the propagation in air from an ultrasonic transmitter to a tested solid plate is modeled in 3D using a closed-form solution. The knowledge of the transducer characteristics (diameter, frequency bandwidth, efficiency in Pa/V) allows the spatial distribution and actual pressure (in Pa) of the acoustic field produced in the air to be predicted, for a given input voltage. This pressure field is applied in turn as a boundary condition in a 3D FE model, to predict the plate response (displacement and stress guided beams) for a given distance between the transmitter and the plate, and for a given angle of orientation of the transmitter with respect to the plate. The FE model is so restricted to modeling of the solid structure only, thus reducing very significantly the number of degrees of freedom and computational time. The material constituting the plate is considered to be an anisotropic and viscoelastic medium. To validate the whole modeling process, an air-coupled ultrasonic transducer is used and oriented at a specific angle chosen for generating one specific Lamb mode guided along a composite plate sample, and a laser probe measures the normal velocity at different locations on the surface of the plate. In the field of NDT, it is generally suitable to excite a pure Lamb mode in order to ease the interpretation of received signals that would represent waves scattered by defects. After a validation step, the numerical model is then used to investigate the effect of the material anisotropy on the purity of the incident guided mode.     

Wave based analysis of the Green's function for a layered cylindrical shell

Cylindrical shells composed of concentric layers may be designed to affect the way that elastic waves are generated and propagated, particularly when some layers are anisotropic. To aid the design process, the present work develops a wave based analysis of the Green's function for a layered cylindrical shell in which the response is given as a sum of waves propagating in the axial coordinate. The analysis assumes linear Hookean materials for each layer. It uses finite element discretizations in the radial coordinate and Fourier series expansions in the circumferential coordinate, leading to linear equations in the axial wavenumber domain that relate shell displacements and forces. Inversion to the axial domain is accomplished via a state-space formulation that is evaluated using residue integration. The resulting expression for the Green's function for each circumferential harmonic is a summation over the natural waves of the shell. The finite element discretization in the radial direction allows the approach to be used for arbitrarily thick shells. The approach is benchmarked to results from an isotropic shell and numerical examples are given for a shell composed of a fiber-reinforced material. The numerical examples illustrate the effect of fiber orientation on the Green's function.     

General approach for finding the characteristics of chiral layered media

In this paper, a systematic approach is presented for finding the resonant frequency, damping factor and EM field distribution of a layered planar structure consisting of chiral media and having arbitrary conductor sheets at the interfaces. The analysis is based on the spectral dyadic Green's function followed by a Ritz-Galerkin's process, using matrix analysis in the Fourier transform spectral domain. The great value of our formulation lies in that it provides methodology for the solutions to chiral layered media problems.     

Inclusions and inhomogeneities in electroelastic media with hexagonal symmetry

For a long time, the absence of explicit Green's functions (fundamental solutions) for electroelastic media has hindered progress in the modelling of the properties of piezoelectric materials. Michelitsch's recently derived explicit electroelastic Green's function for the infinite medium with hexagonal symmetry (transversely isotropic medium) [4] is used here to obtain compact closed-form expressions for the electroelastic analogue of the Eshelby tensor for spheroidal inclusions. This represents a key quantity for the material properties of piezoelectric solids and analysis of the related electroelastic fields in inclusions. For the limiting case of continuous fibers our results coincide with Levin's expressions [8]. The derived method is useful for an extension to non-spheroidal inclusions or inhomogeneities having an axis of rotational symmetry parallel to the hexagonal c-axis. Received 14 September 1999     

An Analytical Method for Calculating P-SV Acoustical Field Excited by a Piezoelectric Strip with Finite Width and Thickness

A piezoelectric strip with finite width and thickness is placed on top of an isotropic elastic half-space. Acoustical field can be excited when a voltage is across the piezoelectric strip. An analytical method is presented to calculate the acoustical field by the dynamics characteristics of the piezoelectric strip. Considering the piezoelectric strip as an anisotropic material of the 6 mm-type crystal system, we study the two-dimensional P-SV acoustical fields inside the piezoelectric strip and the isotropic half-space. The displacement and stress distributions are analysed thoroughly. The effects of the width and thickness of the piezoelectric strip and other parameters on the acoustical field are also analysed.     

Two-dimensional isotropic scattering in a finite thick cylindrical medium exposed to a laser beam

Graphical and tabular results are presented for the back-scattered intensity from a finite two-dimensional cylindrical medium exposed to a Gaussian beam of radiation. Also, results for the source function and flux at the boundaries are presented. The influence of optical thickness and albedo are most pronounced at large optical radii. The semi-infinite results can be used to approximate the finite case for small optical radii. Ranges for single, double, and multiple scattering are discussed. For locations far from the incident beam, the results can be expressed in terms of universal functions independent of beam size. A method is presented for extending the isotropic results to the anisotropic case where the phase function is made up of a spike superimposed on an otherwise isotropic phase function.     

Vibration of moderately thick square orthotropic stepped thickness plates

Numerous studies that address the vibration of stepped thickness plates are reported in the literature. Predominately, classical plate theory has been used to formulate studies for both isotropic and anisotropic stepped plates. Mindlin plate theory has been employed to obtain results for thick isotropic stepped thickness plates. Exact solutions, Rayleigh-Ritz, differential quadrature and finite element methods have been employed to compute results for frequency of vibration. Results for frequency of vibration for thick orthotropic stepped thickness plates are presented here using orthorhombic material properties of aragonite. The finite element method has been used to compute frequencies and determine mode shapes for simply supported and clamped square Mindlin plates.