Research Papers

An Improved Technique for Elastodynamic Green's Function Computation for Transversely Isotropic Solids

[+] Author and Article Information
Samaneh Fooladi

Department of Mechanical Engineering,
Penn State University,
The Behrend College,
Erie, PA 16510;
Department of Aerospace and Mechanical Engineering,
University of Arizona,
Tucson, AZ 85721
email: sxf567@psu.edu

Tribikram Kundu

Department of Aerospace and Mechanical Engineering;
Department of Civil and
Architectural Engineering and Mechanics,
University of Arizona,
Tucson, AZ 85721
e-mail: tkundu@email.arizona.edu

Manuscript received April 3, 2019; final manuscript received April 22, 2019; published online May 21, 2019. Special Editor: Wieslaw Ostachowicz.

ASME J Nondestructive Evaluation 2(2), 021005 (May 21, 2019) (7 pages) Paper No: NDE-19-1017; doi: 10.1115/1.4043605 History: Received April 03, 2019; Accepted April 22, 2019

Elastodynamic Green's function for anisotropic solids is required for wave propagation modeling in composites. Such modeling is needed for the interpretation of experimental results generated by ultrasonic excitation or mechanical vibration-based nondestructive evaluation tests of composite structures. For isotropic materials, the elastodynamic Green’s function can be obtained analytically. However, for anisotropic solids, numerical integration is required for the elastodynamic Green's function computation. It can be expressed as a summation of two integrals—a singular integral and a nonsingular (or regular) integral. The regular integral over the surface of a unit hemisphere needs to be evaluated numerically and is responsible for the majority of the computational time for the elastodynamic Green's function calculation. In this paper, it is shown that for transversely isotropic solids, which form a major portion of anisotropic materials, the integration domain of the regular part of the elastodynamic time-harmonic Green's function can be reduced from a hemisphere to a quarter-sphere. The analysis is performed in the frequency domain by considering time-harmonic Green's function. This improvement is then applied to a numerical example where it is shown that it nearly halves the computational time. This reduction in computational effort is important for a boundary element method and a distributed point source method whose computational efficiencies heavily depend on Green's function computational time.

Copyright © 2019 by ASME
Your Session has timed out. Please sign back in to continue.


Schmerr, L., 2016, Fundamentals of Ultrasonic Nondestructive Evaluation, A Modeling Approach (Springer Series in Measurement Science and Technology), Springer International Publishing, Switzerland.
Manolis, G. D., and Dineva, P. S., 2015, “Elastic Waves in Continuous and Discontinuous Geological Media by Boundary Integral Equation Methods: A Review,” Soil Dyn. Earthq. Eng., 70, pp. 11–29. [CrossRef]
Placko, D., and Kundu, T., 2007, DPSM for Modeling Engineering Problems, John Wiley & Sons, Hoboken, NJ.
Kundu, T., Placko, D., Rahani, E. K., Yanagita, T., and Dao, C. M., 2010, “Ultrasonic Field Modeling: A Comparison Between Analytical, Semi-Analytical and Numerical Techniques,” IEEE Trans. Ultrason. Ferroelectr. Freq. Control, 57(12), pp. 2795–2807. [CrossRef] [PubMed]
Thomson, W., 1848, “Note on the Integration of the Equations of Equilibrium of an Elastic Solid,” Cambridge Dublin Math. J., 1, pp. 76–99.
Fredholm, I., 1900, “Sur les équations de L’équilibre D’un Corps Solide élastique,” Acta Math. Stockholm, 23, pp. 1–42. [CrossRef]
Leibfried, G., 1953, “Versetzungen in Anisotropem Material,” Z. Phys. 135(1), pp. 23–43. [CrossRef]
Barnett, D. M., 1972, “The Precise Evaluation of Derivatives of the Anisotropic Elastic Green’s Functions,” Phys. Status Solidi (b), 49(2), pp. 741–748. [CrossRef]
Wang, C. Y., and Achenbach, J. D., 1994, “Elastodynamic Fundamental Solutions for Anisotropic Solids,” Geophys. J. Int., 118(2), pp. 384–392. [CrossRef]
Wang, C. Y., and Achenbach, J. D., 1995, “Three-Dimensional Time-Harmonic Elastodynamic Green’s Functions for Anisotropic Solids,” Proc. R. Soc. A, 449(1937), pp. 441–458. [CrossRef]
Tewary, V. K., 1995, “Computationally Efficient Representation for Elastostatic and Elastodynamic Green’s Functions for Anisotropic Solids,” Phys. Rev. B Condens. Matter, 51(22), pp. 15695–15702 (Third Series). [CrossRef] [PubMed]
Furukawa, A., Saitoh, T., and Hirose, S., 2014, “Convolution Quadrature Time-Domain Boundary Element Method for 2-D and 3-D Elastodynamic Analyses in General Anisotropic Elastic Solids,” Eng. Anal. Bound. Elem., 39, pp. 64–74. [CrossRef]
Igumnov, L. A., and Markov, I. P., 2016, “Boundary-Element Modeling of Three-Dimensional Anisotropic Viscoelastic Solids,” Advanced Materials, I. A. Parinov, S. H. Chang, and V. Y. Topolov, eds. (Springer Proceedings in Physics 175), Springer International Publishing, Basel, Switzerland, pp. 517–526 (Chap. 36).
Fooladi, S., and Kundu, T., 2017, “Ultrasonic Field Modeling in Anisotropic Materials by Distributed Point Source Method,” Ultrasonics, 78, pp. 115–124. [CrossRef] [PubMed]
Fooladi, S., and Kundu, T., 2018, “Distributed Point Source Modeling of the Scattering of Elastic Waves by a Circular Cavity in an Anisotropic Half-Space,” Ultrasonics, 94, pp. 264–280. [CrossRef] [PubMed]
Fooladi, S., 2016, “Numerical Implementation of Elastodynamic Greens Function for Anisotropic Media,” Master’s thesis, University of Arizona, Tucson, AZ.


Grahic Jump Location
Fig. 1

Integration over the unit hemisphere for computing the regular part of the dynamic Green's function: (a) 3D view and (b) top view

Grahic Jump Location
Fig. 2

Two unit vectors n(1) and n(2) on the unit sphere with the same azimuthal angle

Grahic Jump Location
Fig. 3

The rotation of the coordinate system from (x1, x2, x3) to (x1′, x2′, x3′)

Grahic Jump Location
Fig. 4

The components GR(1, 1) and GR(1, 2) from the regular part of Green's function

Grahic Jump Location
Fig. 5

The derivative of the components GR(1, 1) and GR(1, 2) from the regular part of Green's function with respect to coordinate x1

Grahic Jump Location
Fig. 6

Ratio of time of computation for quarter-sphere to that of hemisphere



Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In