SH Wave Scattering Problems for Multiple Orthotropic Elliptical Inclusions

Abstract

A volume integral equation method (VIEM) is applied for the effective analysis of elastic wave scattering problems in unbounded solids containing general anisotropic inclusions. It should be noted that this numerical method does not require use of Green's function for anisotropic inclusions to solve this class of problems since only Green's function for the unbounded isotropic matrix is necessary for the analysis. This new method can also be applied to general two-dimensional elastodynamic problems involving arbitrary shapes and numbers of anisotropic inclusions. A detailed analysis of SH wave scattering problems is developed for an unbounded isotropic matrix containing multiple orthotropic elliptical inclusions. Numerical results are presented for the displacement fields at the interfaces of the inclusions in a broad frequency range of practical interest. Through the analysis of plane elastodynamic problems in an unbounded isotropic matrix with multiple orthotropic elliptical inclusions, it is established that this new method is very accurate and effective for solving plane elastic problems in unbounded solids containing general anisotropic inclusions of arbitrary shapes.

1. Introduction

A micrograph of the cross-section of a phosphate glass fiber/polymer composite is shown in Figure 1 [1]. Figure 1 indicates that the fibers are close to ellipses and the major axis of the elliptical fibers is not aligned in any one direction. Furthermore, the micrograph cross section of the glass fiber composite in Figure 2 [2] shows that the fibers are triangular in shape.

Figure 1:

Phosphate glass fiber/polymer composite cross section [1]. Used with permission from Journal of Materials Science: Materials in Medicine.

Figure 2:

Micrograph cross-section of a unidirectional 60 μm triangular glass fiber composite (V_f = 0.5) [2]. Used with permission from Composites Science and Technology.

A number of analytical techniques are available for solving stress analysis of isotropic inclusion problems when the geometry of the inclusions is simple (i.e., cylindrical, spherical, or ellipsoidal) and when they are well separated [3 –6]. However, these approaches cannot be applied to more general problems where the inclusions are both anisotropic and arbitrary in shape when their concentration is high. Thus, analysis of elastic wave scattering problems in heterogeneous solids often requires the use of numerical techniques based on the finite element method (FEM) or boundary element method (BIEM). Unfortunately, both methods encounter limitations in dealing with elastic wave scattering problems involving multiple anisotropic inclusions of arbitrary shapes. However, it has been demonstrated that a newly developed numerical method based on a volume integral equation formulation can overcome such difficulties in solving this class of inclusion problems [7, 8]. In contrast to the conventional boundary integral equation method (BIEM), where the infinite medium Green's functions for both the matrix material and the anisotropic inclusion material are needed, the present method does not require the latter. Since elastodynamic Green's functions for anisotropic media are extremely difficult to calculate, the present method offers a definite advantage over the boundary integral equation method. In addition, the VIEM is not sensitive to the geometry or concentration of the inclusions. Moreover, in contrast to the finite element method (FEM), where the full domain needs to be discretized, the VIEM requires discretization of the inclusions only.

Therefore, in order to investigate the influence of orthotropic elliptical inclusions on the interfacial field, a detailed analysis of the displacement field at the interface between the matrix and the central inclusion is first carried out for an unbounded isotropic matrix containing five orthotropic elliptical inclusions. The major axis of each elliptical inclusion is 0° or 90° with the x-axis and the aspect ratio is assumed to be 1.0 (circular) and 0.8. The inclusion separation distance (s) is assumed to be 3.24a and 2.51a (a: the radius of each circular inclusion or the major radius of each elliptical inclusion). It should be noted that the central inclusion is selected since the interaction effects of the multiple inclusions are likely to be highest in value there. The incident wave represents SH waves propagating parallel to the x-axis with particle motion being parallel to the z-axis.

It is demonstrated that this new method is very accurate and effective for solving plane elastodynamic problems in unbounded solids containing general anisotropic inclusions of arbitrary shapes. In forthcoming papers, in order to accurately predict the failure and damage mechanisms in real composites, the formulations developed herein will be used to calculate other quantities of practical interest in realistic models of composites containing multiple strong anisotropic heterogeneities [9].

2. Volume Integral Equation Method (VIEM)

The geometry of the general elastodynamic problem considered here is shown in Figure 3 where an unbounded isotropic elastic solid containing a number of isotropic or anisotropic inclusions of arbitrary shape are subjected to prescribed dynamic loading at infinity. The symbols $ρ^{(1)}$ and ${c_{i j k l}}^{(1)}$ denote the density and the elastic stiffness tensor of the inclusion. ρ⁽²⁾ and c_ijkl⁽²⁾ denote the density and the elastic stiffness tensor of the unbounded matrix material. The matrix is assumed to be homogeneous and isotropic so that c_ijkl⁽²⁾ is a constant isotropic tensor, while ${c_{i j k l}}^{(1)}$ can be arbitrary; that is, the inclusions may, in general, be inhomogeneous and anisotropic. The interfaces between the inclusions and the matrix are assumed to be perfectly bonded ensuring continuity of the displacement and stress vectors.

Figure 3:

Geometry of the general elastodynamic problem.

Let u_m^o(x, ω)e^−iωt denote the mth component of the displacement vector due to the incident field at x in the absence of the inclusions. Also, let u_m(x, ω)e^−iωt denote the same in the presence of the inclusions, where ω is the circular frequency of the waves. In what follows, the common time factor e^−iωt and the explicit dependence on ω for all field quantities will be suppressed.

It has been shown in Mal and Knopoff [10] that the elastodynamic displacement, u_m(x), in the composite satisfies the volume integral equation

u_{m} (x) = u_{m}^{o} (x) + \int_{R}^{} [δ ρ ω^{2} g_{i}^{m} (ξ, x) u_{i} (ξ) - δ c_{i j k l} g_{i, j}^{m} (ξ, x) u_{k, l} (ξ)] d ξ,

(1)

where the integral is over the domain occupied by the inclusions, δρ and δc_ijkl represent contrasts in the density and elastic tensor of the inclusions and the matrix (i.e., δρ = ρ⁽¹⁾ – ρ⁽²⁾ and δc_ijkl = c_ijkl⁽¹⁾ – c_ijkl⁽²⁾), and g_i^m(ξ, x)e^−iωt is elastodynamic Green's function for the unbounded homogeneous matrix material; that is, ${g_{i}}^{m} (ξ, x) e^{- i ω t}$ represents the ith component of the displacement at ξ due to unit concentrated force, e_me^−iωt, at x in the mth direction. In (1), the summation convention and comma notation have been used and the differentiations are with respect to the integration variable, ξ_i. It should also be pointed out that the integrand is nonzero only within the inclusions, since δc_ijkl = 0, outside the inclusions.

If x ∊ R (within the inclusions), then (1) is an integro-differential equation for the unknown displacement vector u(x). It can, therefore in principle, be determined through the solution of this equation. An algorithm for the solution of (1) was developed by Lee and Mal [7] by discretizing the inclusions using conventional finite elements. Once u(x) within the inclusions is determined, the displacement field in the matrix can be calculated from (1) by evaluating the integral for x ∉ R. The stress field within and outside the inclusions can also be determined in a similar manner. Extensive details of the numerical treatment of (1) for plane elastodynamic problems can be found in [7, 8]. Further explanation of the elastostatic volume integral equation method for isotropic inclusions in an isotropic matrix can also be found in Section 4.3 “Volume Integral Equation Method” by Buryachenko [11].

3. Scattering of SH Waves

3.1. Scattering of SH Waves in an Unbounded Isotropic Matrix Containing Multiple Orthotropic Circular Inclusions

We first consider multiple orthotropic circular inclusions in an unbounded isotropic matrix. The incident wave is represented as an SH wave with particle motion along the z-axis (see Figure 4). Let the coordinate axes x₁(x), x₂(y), x₃(z) be taken parallel to the symmetry axes of the orthotropic material, let ρ₁ and c₄₄, c₅₅ denote the density and elastic constants, respectively, of the orthotropic inclusion, and let ρ₂ and μ₂ denote the density and shear modulus of the isotropic matrix, respectively.

Figure 4:

SH wave interaction with five orthotropic elliptical inclusions in an unbounded isotropic matrix.

For SH waves, the volume integral equation (1) reduces to

u_{3} (x) = u_{3}^{o} (x) + \int_{R}^{} [δ ρ ω^{2} g_{3}^{3} (ξ, x) u_{3} (ξ) - δ c_{55} g_{3,1}^{3} u_{3,1} - δ c_{44} g_{3,2}^{3} u_{3,2}] d ξ_{1}^{} d ξ_{2}^{},

(2)

where u₃(x) is the antiplane displacement component (i.e., the displacement component in the z-axis), δρ = ρ₁ – ρ₂, δc₄₄ = c₄₄ – μ₂, and δc₅₅ = c₅₅ – μ₂.

In (2), g_i^m(ξ, x) is Green's function for the unbounded isotropic matrix material. The Green's function for the SH problem is given by Kitahara [12] as

g_{3}^{3} (ξ, x) = \frac{i H_{0}^{} (k_{2}^{} r)}{4 μ_{2}^{}},

(3)

where r = |x – ξ|, k₂ = ω/β₂ is the S wavenumber, β₂ is the shear wave speed in the unbounded matrix material, and H₀ is the Hankel function of the first kind of the zeroth order. Thus, the volume integral equation method does not require use of the Green's function for the orthotropic material of the inclusions. This is in contrast to the boundary integral equation method, where the infinite medium Green's functions for both the isotropic matrix and the orthotropic inclusions are involved in the formulation of the equations.

Finite element discretization of the inclusions in (2) results in a system of linear algebraic equations for the unknown nodal displacements inside the orthotropic inclusions. Once the displacement field, u₃(x), within the inclusions is determined, that outside the inclusions can be obtained from (2) by evaluating the integrals. The stress fields within and outside the orthotropic inclusions can also be readily determined. Further details of this numerical treatment can be found in Lee and Mal [7, 8].

It should be noted that the singularities in VIEM are weaker (integrable) than those in BIEM, where they are of the Cauchy type. In response, we have used the direct integration scheme introduced by Cerrolaza and Alarcon [13], Li et al. [14], and Lu and Ye [15] after suitable modifications to address these singularity concerns. A description of the modified method used in the discretization of the volume integral equation is given by Lee and Mal [7, 16].

In the SH wave case, the incident wave is given as

u_{3}^{o} = e_{}^{i k_{2}^{} x},

(4)

where k₂ is the S wavenumber in the matrix material. The elastic constants for the materials of the isotropic matrix and the orthotropic inclusions are listed in Table 1. Two different elastic constants for the orthotropic inclusions are considered: for “orthotropic inclusion #1,” the density values, c₄₄ and c₅₅ in the orthotropic inclusions, are smaller than those in the matrix, and for “orthotropic inclusion #2,” the density values, c₄₄ and c₅₅ in the orthotropic inclusions, are greater than those in the matrix. The calculations are carried out for three different normalized wavenumbers (i.e., k₂a = 3.14, 5.03, and 7.85). The ratios of the corresponding wavelengths to the radius of each circular inclusion are approximated to be between 0.8 and 2.0. This represents a useful range of intermediate frequencies for dynamic loading and ultrasonic testing.

Table 1:

Material properties of the isotropic matrix and the orthotropic inclusions.

	Isotropic matrix	Orthotropic inclusion #1	Orthotropic inclusion #2

ρ (g/cm³)	2.706	1.6	3.181
μ (GPa)	37.88	—	—
c₄₄ (GPa)	37.88	5.81	46.48
c₅₅ (GPa)	37.88	7.46	59.68

Figure 5 shows a typical discretized model used in the volume integral equation method [17]. Standard eight-node quadrilateral and six-node triangular elements were used. The number of elements in each circular inclusion used in the VIEM was 3,600 and was determined based on a convergence test. Figure 6 shows real and imaginary parts of the displacement component in the z-axis, u₃, at three separate interfaces: (a) for the single orthotropic circular inclusion (#1), (b) for the central inclusion with five orthotropic circular inclusions (#1) where s = 3.24a, and (c) for the central inclusion with five orthotropic circular inclusions (#1) where s = 2.51a, all of which use the VIEM for three different frequencies (i.e., k₂a = 3.14, 5.03, and 7.85). Figure 7 shows real and imaginary parts of the displacement component in the z-axis, u₃, at three separate interfaces: (a) for the single orthotropic circular inclusion (#2), (b) for the central inclusion with five orthotropic circular inclusions (#2) where s = 3.24a, and (c) for the central inclusion with five orthotropic circular inclusions (#2) where s = 2.51a, all of which use the VIEM for three different frequencies (i.e., k₂a = 3.14, 5.03, and 7.85).

Figure 5:

A typical discretized model in the volume integral equation method.

Figure 6:

Real and imaginary parts of the displacement component in the z-axis, u₃, at the interface for: (a) the single orthotropic circular inclusion (#1), (b) the central inclusion with five orthotropic circular inclusions (#1) where s = 3.24a, and (c) the central inclusion with five orthotropic circular inclusions (#1) where s = 2.51a. The frequencies used are k₂a = 3.14, 5.03, and 7.85 where k₂ is the S wavenumber in the unbounded isotropic matrix and a is the radius of the circular inclusion. The inclusion cited represents orthotropic inclusion #1 in Table 1.

Figure 7:

Real and imaginary parts of the displacement component in the z-axis, u₃, at the interface for: (a) the single orthotropic circular inclusion (#2), (b) the central inclusion with five orthotropic circular inclusions (#2) where s = 3.24a, and (c) the central inclusion with five orthotropic circular inclusions (#2) where s = 2.51a. The frequencies used are k_2a = 3.14, 5.03, and 7.85 where k₂ is the S wavenumber in the unbounded isotropic matrix and a is the radius of the circular inclusion. The inclusion cited represents orthotropic inclusion #2 in Table 1.

3.2. Verification of the Viem Solution

To the best of the authors' knowledge, neither an analytical solution to this problem nor a closed form solution for the two-dimensional time-harmonic elastodynamic Green's function for the orthotropic material is presently available in the literature. Thus, in order to compare the calculated results with available analytical solutions, simple packing sequences (hexagonal and square) and isotropic circular inclusions in an unbounded isotropic elastic solid are considered (see Figure 8). Let ρ₁ and μ₁ denote the density and elastic constants, respectively, of the isotropic inclusion and let ρ₂ and μ₂ denote the density and shear modulus of the isotropic matrix, respectively. The material properties used are those for typical graphite/epoxy composites and are given in Table 2.

Table 2:

Material properties of graphite/epoxy composites.

Material	Density (g/cm³)	μ (GPa)

Graphite	1.79	27.580
Epoxy	1.26	1.595

Figure 8:

SH wave interaction with multiple isotropic inclusions.

For SH waves, the volume integral equation (1) reduces to

u_{3} (x) = u_{3}^{o} (x) + \int_{R}^{} [δ ρ ω^{2} g_{3}^{3} (ξ, x) u_{3} (ξ) - δ μ g_{3,1}^{3} u_{3,1} - δ μ g_{3,2}^{3} u_{3,2}] d ξ_{1}^{} d ξ_{2}^{},

(5)

where u₃(x) is the antiplane displacement component, δρ = ρ₁ – ρ₂, δμ = μ₁ – μ₂.

In (5), g_i^m(ξ, x) is Green's function for the unbounded isotropic matrix material. Finite element discretization of the inclusions in (5) results in a system of linear algebraic equations for the unknown nodal displacements inside the isotropic inclusions. Once the displacement field, u₃(x), within the inclusions is determined, that outside the inclusions can be obtained from (5) by evaluating the integrals. The stress fields within and outside the isotropic inclusions can also be readily determined.

The incident wave is given as

u_{3}^{o} = e_{}^{i k_{2}^{} x},

(6)

where k₂ is the S wavenumber in the matrix material.

Figure 9 shows discretized models used in the VIEM. The standard eight-node quadrilateral and six-node triangular elements were used in the VIEM. The total number of elements for each inclusion used in the VIEM was 256 and was determined based on a convergence test. Tables 3 and 4 show calculated average strains in the central isotropic fiber for real and imaginary parts as determined by an analytical method using the generalized self-consistent model (GSCM) [5] and the volume integral equation method for both (i) the square packing of 9 and 25 isotropic circular inclusions and (ii) the hexagonal packing of 7 and 19 isotropic circular inclusions. The material used was standard graphite/epoxy composites with an inclusion volume fraction of 0.6 and a frequency of 10 MHz. The percentage differences in the two sets of results were less than 1% in all cases [7].

Table 3:

Calculated average strains in the central fiber for the real parts using both analytical and volume integral equation methods. The material is graphite/epoxy with a volume concentration of 0.6 and a frequency of 10 MHz.

Number of inclusions	Average strain (Re)
Number of inclusions	Analytical	VIE	Difference

7	0.09416	0.09474	0.61%
9	0.09035	0.09048	0.14%
19	0.08072	0.08089	0.21%
25	0.08047	0.08070	0.29%

Table 4:

Calculated average strains in the central fiber for the imaginary parts using both analytical and volume integral equation methods. The material is graphite/epoxy with a volume concentration of 0.6 and a frequency of 10 MHz.

Number of inclusions	Average strain (Im)
Number of inclusions	Analytical	VIE	Difference

7	−0.0270	−0.0272	0.74%
9	−0.0329	−0.0330	0.30%
19	−0.0519	−0.0523	0.77%
25	−0.0645	−0.0652	1.10%

Figure 9:

Discretized models for multiple inclusion problems used in the volume integral equation method.

3.3. Scattering of SH Waves in an Unbounded Isotropic Matrix Containing Multiple Orthotropic Elliptical Inclusions

We next consider multiple orthotropic elliptical inclusions in an unbounded isotropic matrix (see Figure 4). For elliptical cylindrical inclusions, the value of the aspect ratio for each elliptical inclusion was chosen to be 0.8 (see Figures (4) and (5)). Two different types of elliptical inclusions are considered where (1) the major axis of each inclusion coincides with the x axis (see Figure 4) and (2) the major axis of each inclusion coincides with the y axis (see Figure 5). The general features of the volume integral equation method for multiple orthotropic elliptical inclusions are similar to those for multiple orthotropic circular inclusions.

In the SH wave case, the incident wave is given as

u_{3}^{o} = e_{}^{i k_{2}^{} x},

(7)

where k₂ is the S wavenumber in the matrix material. The elastic constants for the materials of the isotropic matrix and the orthotropic inclusions are listed in Table 1. The calculations are carried out for three different normalized wavenumbers (i.e., k₂a = 3.14, 5.03, and 7.85). The ratios of the corresponding wavelengths to the major radius of each elliptical inclusion are approximated to be between 0.8 and 2.0. This is a useful range of intermediate frequencies for dynamic loading and ultrasonic testing.

When the major axis of each elliptical inclusion is parallel to the x-axis, Figure 10 shows real and imaginary parts of the displacement component in the z-axis, u₃, at three separate interfaces: (a) for the single orthotropic elliptical inclusion (#1), (b) for the central inclusion with five orthotropic elliptical inclusions (#1) where s = 3.24a, and (c) for the central inclusion with five orthotropic elliptical inclusions (#1) where s = 2.51a, all of which use the volume integral equation method for three different frequencies (i. e., k₂a = 3.14, 5.03, and 7.85). Figure 11 shows real and imaginary parts of the displacement component in the z-axis, u₃, at three separate interfaces: (a) for the single orthotropic elliptical inclusion (#2), (b) for the central inclusion with five orthotropic elliptical inclusions (#2) where s = 3.24a, and (c) for the central inclusion with five orthotropic elliptical inclusions (#2) where s = 2.51a, all of which use the VIEM for three different frequencies (i.e., k₂a = 3.14, 5.03, and 7.85).

Figure 10:

Real and imaginary parts of the displacement component in the z-axis, u₃, at the interface for: (a) the single orthotropic elliptical inclusion (#1), (b) the central inclusion with five orthotropic elliptical inclusions (#1) where s = 3.24a, and (c) the central inclusion with five orthotropic elliptical inclusions (#1) where s = 2.51a. The frequencies used are k₂a = 3.14, 5.03, and 7.85 where k₂ is the S wavenumber in the unbounded isotropic matrix and a is the major radius of the elliptical inclusion. The inclusion cited represents orthotropic inclusion #1 in Table 1.

Figure 11:

Real and imaginary parts of the displacement component in the z-axis, u₃, at the interface for: (a) the single orthotropic elliptical inclusion, (b) the central inclusion with five orthotropic elliptical inclusions (#2) where s = 3.24a, and (c) the central inclusion with five orthotropic elliptical inclusions (#2) where s = 2.51a. The frequencies used are k₂a = 3.14, 5.03, and 7.85 where k₂ is the S wavenumber in the unbounded isotropic matrix and a is the major radius of the elliptical inclusion. The inclusion cited represents orthotropic inclusion #2 in Table 1.

When the major axis of each elliptical inclusion is parallel to the y-axis, Figure 12 shows real and imaginary parts of the displacement component in the z-axis, u₃, at three separate interfaces: (a) for the single orthotropic elliptical inclusion (#1), (b) for the central inclusion with five orthotropic elliptical inclusions (#1) where s = 3.24a, and (c) for the central inclusion with five orthotropic elliptical inclusions (#1) where s = 2.51a, all of which use the VIEM for three different frequencies (i.e., k₂a = 3.14, 5.03, and 7.85). Figure 13 shows real and imaginary parts of the displacement component in the z-axis, u₃, at three separate interfaces: (a) for the single orthotropic elliptical inclusion (#2), (b) for the central inclusion with five orthotropic elliptical inclusions (#2) where s = 3.24a, and (c) for the central inclusion with five orthotropic elliptical inclusions (#2) where s = 2.51a, all of which use the VIEM for three different frequencies.

Figure 12:

Real and imaginary parts of the displacement component in the z-axis, u₃, at the interface for: (a) the single orthotropic elliptical inclusion, (b) the central inclusion with five orthotropic elliptical inclusions (#1) where s = 3.24a, and (c) the central inclusion with five orthotropic elliptical inclusions (#1) where s = 2.51a. The frequencies used are k₂a = 3.14, 5.03, and 7.85 where k₂ is the S wavenumber in the unbounded isotropic matrix and a is the major radius of the elliptical inclusion. The inclusion cited represents orthotropic inclusion #1 in Table 1.

Figure 13:

For different shapes of orthotropic elliptical inclusions, the interaction effect of the inclusions on the displacement component in the z-axis, u₃, at the interface of the central inclusion appeared to be significant for orthotropic elliptical inclusion #1. In this case, the density values, c₄₄ and c₅₅, in the inclusions were determined to be smaller than those in the matrix. However, the interaction effect of the inclusions appeared to be minimal for orthotropic elliptical inclusion #2. In that case, the density values, c₄₄ and c₅₅, in the inclusions were found to be greater than those in the matrix.

4. Concluding Remarks

The volume integral equation method was applied to the solution of SH wave scattering problems with multiple orthotropic inclusions involving variations in orientation angle and inclusion separation distances for three different frequencies in an unbounded isotropic matrix. The aspect ratio of each elliptical inclusion was given as 1.0 (circular) and 0.8. It was determined that the displacement component in the z-axis, u₃, at the interface of the central inclusion for orthotropic elliptical inclusion #1 (soft inclusion) differed significantly. However, for orthotropic elliptical inclusion #2 (hard inclusion), the interaction effect of orthotropic elliptical inclusions on the displacement component, u₃ (in the z-axis), at the interface of the central inclusion appeared to be minimal.

The main advantage of this technique over those based on finite elements is that it requires discretization of the inclusions only as opposed to having to discretize the entire domain. This approach is similar to the boundary integral equation method except for the presence of the volume integral over the inclusions instead of the surface integrals over the two sides of the interface. It should be noted that if the medium contains a small number of isotropic inclusions, this method may not be as advantageous when compared to the BIEM. However, in the presence of multiple nonsmooth inclusions, the BIEM numerical treatment becomes noticeably cumbersome. Specifically, in elastodynamic problems involving multiple anisotropic inclusions, BIEM numerical treatment becomes extremely difficult since closed-form expressions for elastodynamic Green's functions for anisotropic media are currently not available. However, the volume integral equation method does not require the use of Green's functions for anisotropic inclusions. Furthermore, since standard finite elements are used in the VIEM, it is easier and more convenient to handle multiple nonsmooth anisotropic inclusions. Therefore, the formulations developed in this paper in conjunction with other useful methods [18, 19] can be used to calculate the dynamic stress intensity factors and other quantities of practical interest in realistic models of materials containing strong heterogeneities.

Through the analysis of plane elastodynamic problems in an unbounded isotropic matrix with multiple orthotropic elliptical inclusions, it is established that this new method is very accurate and effective for solving plane elastodynamic problems in unbounded solids containing general anisotropic inclusions of arbitrary shapes. In forthcoming papers, the formulations developed in this research will be used to calculate other quantities of practical interest in realistic models of composites containing multiple strong anisotropic heterogeneities [9].

Footnotes

Appendix

For three-dimensional elastodynamic problems, u₃(x) in the volume integral equation (1) for fully anisotropic inclusions can be expressed in the form

where u₁(x), u₂(x), and u₃(x) are the three-dimensional displacement components, ω is the circular frequency of the waves, and δρ = ρ⁽¹⁾ – ρ⁽²⁾, where ρ⁽¹⁾ denotes the density of the inclusions while ρ⁽²⁾ represents that for the matrix material. In addition, ${δ c}_{α β} = c_{α β}^{(1)} - c_{α β}^{(2)} (α, β = 1, 6$ ), where c_αβ⁽¹⁾ denote the elastic stiffness constants of the inclusions while c_αβ⁽²⁾ represent those for the matrix material. In (A.1), g_i^m(ξ, x) (i, m = 1, 3) represents the three-dimensional elastodynamic Green's function for the unbounded isotropic matrix material and is given by Pao and Varatharajulu [20] as

where $R = | x - ξ | = \sqrt{{(x_{1} - ξ_{1})}_{2} + {(x_{2} - ξ_{2})}_{2} + {(x_{3} - ξ_{3})}_{2}}$ , k₁ = ω/c₁, and k₂ = ω/c₂ where c₁ is the P wave speed while c₂ is the S wave speed.

Acknowledgments

This research was supported by Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education, Science and Technology (Grant no. 2012-0007636) and Business for Cooperative R&D between Industry, Academy, and Research Institute funded by Korea Small and Medium Business Administration in 2013 (Grant no. C0096465). The authors would also like to express their sincere appreciation to Professor Jason Buschman for his excellent English proofreading and to acknowledge the support from Korea Institute of Science and Technology Information (KISTI) supercomputing center through the strategic support Program for the supercomputing application research (no. KSC-2012-C1-06).

References

Brauer

D. S.

Rüssel

Vogt

Weisser

, and Schnabelrauch

, “Degradable phosphate glass fiber reinforced polymer matrices: mechanical properties and cell response,” Journal of Materials Science: Materials in Medicine, vol. 19, no. 1, pp. 121–127, 2008.

Bond

Hucker

Weaver

Bleay

, and Haq

, “Mechanical behaviour of circular and triangular glass fibres and their composites,” Composites Science and Technology, vol. 62, no. 7–8, pp. 1051–1061, 2002.

Bose

S. K.

and Mal

A. K.

, “Longitudinal shear waves in a fiber-reinforced composite,” International Journal of Solids and Structures, vol. 9, no. 9, pp. 1075–1085, 1973.

Bose

S. K.

and Mal

A. K.

, “Elastic waves in a fiber-reinforced composite,” Journal of the Mechanics and Physics of Solids, vol. 22, no. 3, pp. 217–229, 1974.

Yang

R. B.

and Mal

A. K.

, “Phase velocity and attenuation of SH waves in a fiber reinforced composite,” in Review of Progress in Quantitative Nondestructive Evaluation, Thompson

D. O.

and Chimenti

D. E.

, Eds., vol. 12A, chapter 1, pp. 155–162, 1993.

Ying

C. F.

and Truell

, “Scattering of a plane longitudinal wave by a spherical obstacle in an isotropically elastic solid,” Journal of Applied Physics, vol. 27, no. 9, pp. 1086–1097, 1956.

Lee

and Mal

A. K.

, “A volume integral equation technique for multiple scattering problems in elastodynamics,” Applied Mathematics and Computation, vol. 67, no. 1–3, pp. 135–159, 1995.

Lee

, and Mal

, “A mixed volume and boundary integral equation technique for elastic wave field calculations in heterogeneous materials,” Wave Motion, vol. 39, no. 1, pp. 1–19, 2004.

Sumiya

Biwa

, and Haïat

, “Computational multiple scattering analysis of elastic waves in unidirectional composites,” Wave Motion, vol. 50, pp. 253–270, 2013.

10.

Mal

A. K.

and Knopoff

, “Elastic wave velocities in two-component systems,” Journal of the Institute of Mathematics and its Applications, vol. 3, no. 4, pp. 376–387, 1967.

11.

Buryachenko

V. A.

, Micromechanics of Heterogeneous Materials, Springer, New York, NY, USA, 2007.

12.

Kitahara

, Boundary Integral Equation Methods in Eigenvalue Problems of Elastodynamics and Thin Plates, Elsevier, London, UK, 1985.

13.

Cerrolaza

and Alarcon

, “A bi-cubic transformation for the numerical evaluation of the Cauchy Principal Value integrals in boundary methods,” International Journal for Numerical Methods in Engineering, vol. 28, no. 5, pp. 987–999, 1989.

14.

H.-B.

Han

G.-M.

, and Mang

H. A.

, “A new method for evaluating singular integrals in stress analysis of solids by the direct boundary element method,” International Journal for Numerical Methods in Engineering, vol. 21, no. 11, pp. 2071–2098, 1985.

15.

and Ye

T. Q.

, “Direct evaluation of singular integrals in elastoplastic analysis by the boundary element method,” International Journal for Numerical Methods in Engineering, vol. 32, no. 2, pp. 295–311, 1991.

16.

Lee

and Mal

, “A volume integral equation technique for multiple inclusion and crack interaction problems,” Journal of Applied Mechanics, vol. 64, no. 1, pp. 23–30, 1997.

17.

“PATRAN user's manual,” Tech. Rep. version 7.0, MSC/PATRAN, 1998.

18.

Chen

J. T.

Chen

C. T.

Chen

P. Y.

, and Chen

I. L.

, “A semi-analytical approach for radiation and scattering problems with circular boundaries,” Computer Methods in Applied Mechanics and Engineering, vol. 196, no. 25–28, pp. 2751–2764, 2007.

19.

Dravinski

and Sheikhhassani

, “Scattering of a plane harmonic SH wave by a rough multilayered inclusion of arbitrary shape,” Wave Motion, vol. 50, no. 4, pp. 836–851, 2013.

20.

Pao

Y.-H.

and Varatharajulu

, “Huygens' principle, radiation conditions, and integral formulas for the scattering of elastic waves,” Journal of the Acoustical Society of America, vol. 59, no. 6, pp. 1361–1369, 1976.