Research Article  Open Access
Scattering Attenuation of SHWaves in FiberReinforced Composites with Partial Imperfectly Bonded Interfaces
Abstract
The influence of partial imperfectly bonded interfaces between each fiber and the matrix on the scattering attenuation of coherent SHwaves in fiberreinforced composites is investigated. The imperfection of interfaces is modelled using the spring model, in which the level of imperfection is characterized by a parameter , which is called the stiffness of the imperfect interfaces. First, the single scattering of SHwaves by a cylindrical fiber with such a partial imperfectly bonded interface to the matrix is formulated and subsequently solved using the collocation point (CP) method. Later on, based on the analysis of the corresponding single scattering problem, effects of several parameters (i.e., frequency of the incident wave, level of imperfection of interfaces, and width of the imperfect interfaces) on the farfield scattering magnitude, scattering cross section, and scattering attenuation coefficients of coherent SHwaves are shown graphically. The potential application of the current results to nondestructive evaluation of interfaces in composites is also discussed.
1. Introduction
Various interfaces exist in the nature, such as interfaces between different constituents of composite materials and those between distinct phases in a polycrystalline material. For composite materials, the interfaces play very important roles, such as transferring loads from the matrix to fibers. Accordingly, evaluation of the interfaces of composite materials has been the focus of many researchers [1–6].
Currently, the ultrasonic technique has become a major method for the nondestructive evaluation of interfaces, due to the fact that the behavior (i.e., attenuation and wave speed) of elastic waves in composite materials is highly related to the interfaces [1, 2, 7]. It is well known that, before using this ultrasonic technique to evaluate composite materials, the relationship between characteristics of elastic waves in composites and their properties should be established. It is in this spirit that, to date, numerous theoretical models [8–20] for the evaluation of attenuation and wave speed of coherent elastic waves in composite materials have been proposed. It is worthwhile to note that, among these models, there is a simple one used for cases with dilute concentration of inclusions, which is called the independent scattering model [10, 17–19, 21]. The model indicates that the scattering attenuation coefficient of coherent elastic waves in composites can be calculated as , where is the scattering cross section [22] and is the number of fibers/particulates inside a unit area/volume. Several experiments [1, 9] have shown that this independent scattering model has a good performance.
To date, in most of published relevant works, the fibers and the surrounding matrix are usually assumed to be perfectly bonded. However, this is not always true in practice. Therefore, during the past few years, several theoretical models have been developed to model the imperfect situations of interfaces. These models include, but are not limited to, the spring model [23] and the springmass model [24]. Specifically, in the spring model, the stress is assumed to be continuous across the interfaces, while allowing the displacement has a jump across the interfaces. This jump is assumed to be proportional to the stresses at the interfaces and has the relation , where denotes the jump of the displacements across the interface and is the parameter to characterize the imperfection level of the interfaces. It is clear that, when , the spring model becomes the familiar perfect bonding and it is degenerated to the debonding case as . Other intermediate values of represent the “imperfect bonding” between the “perfect bonding” and the “debonding.” The validity of this spring model has been verified in several works [1, 7, 25]. Compared with the spring model, the additional inertia effect of the interphases has been considered in the springmass model. Detailed information on these models has been provided in [26]. This spring model can be also used to model the imperfect interfaces in other physical fields, that is, electric and magnetic [27].
In the past few years, this imperfect bonding effect on the effective properties and the behavior of elastic waves in composite materials/structures have attracted attention of several researchers [1, 28–31]. Most of these works have accounted for the imperfect bonding effect through the spring model. Specifically, for example, the influence of the imperfect bonding between inclusions and the matrix on the effective elastic properties [28, 31, 32] and thermal [33] properties of composites has been investigated. For the behavior of elastic waves in composites, Rokhlin et al. [1] have investigated the imperfect bonding effect on the scattering of SH, P, and SVwaves in fiberreinforced composites. Mal and Bose [30] have studied the phase velocity and attenuation of coherent P and SVwaves in particulate composites taking the imperfect bonding effect into account. It should be pointed out that, in most of these published works, the whole interface between each fiber/particulate and the matrix is usually assumed to enter the imperfect bonding situation simultaneously. However, this is not always the case. Assuming that the interface between each fiber and the matrix deteriorates as a process would be more reasonable. In this regard, LopezRealpozo et al. [31] have investigated the nonuniform imperfect bonding effect on the effective shear stiffness of a periodic fibrous composite. The corresponding work on the effective thermal conductivity has been conducted in [33]. Several other authors [3, 6, 34–38] have investigated the influence of partial debonding of interfaces on elastic wave propagation in composite materials, which is an extreme case of the partial imperfect bonding. To the best of our knowledge, the investigation of the influence of the partial imperfect bonding of interfaces on the behavior of coherent elastic waves in composites is still scarce.
In this work, we aimed to investigate the influence of the partial imperfect bonding of interfaces on the scattering attenuation of coherent SHwaves in fiberreinforced composite materials. The imperfection of interfaces is modelled using the spring model. The scattering attenuation is calculated using the previously introduced independent scattering model. The CP method is proposed to solve the corresponding single scattering problem with a partial imperfectly bonded interface. Compared with the existing techniques, the derivation process and implementation of this CP method are simple and it has a comparable accuracy. Based on this CP method, influence of various aspects of the imperfect interfaces between fibers and the matrix on the scattering attenuation of coherent SHwaves in composites is studied. Such investigation benefits the nondestructive evaluation of the interfaces in composites using the ultrasonic technique, as the results indicate that the working frequency of the scattering attenuation based ultrasonic technique to evaluate the interfaces should be carefully selected.
In the following sections, first, the single scattering of SHwaves by a cylindrical fiber with a partial imperfect bonding to the matrix is described in detail. The CP method is then introduced to solve such a single scattering problem. Based on the analysis of the single scattering problem, extensive parameter analysis on the farfield scattering magnitude, scattering cross section, and scattering attenuation is then conducted. At the end, a short conclusion is presented.
2. Single Scattering of SHWaves by a Fiber with a Partial Imperfect Bonding to the Matrix
2.1. Formulation of the Problem
Figure 1 shows a schematic of the single scattering of a plane SHwave by a cylindrical fiber with radius . As shown in this figure, only a portion of the fiber is perfectly bonded to the matrix while the remaining is not. The symmetric axis of the imperfect bonding is at an angle with the horizontal axis. The range of the imperfect bonding is represented by angle , as indicated by the red solid curve in Figure 1. It is clear that, when , the whole fiber is purely perfectly bonded to the matrix and the whole interface becomes a pure imperfect bonding as . The incident wave is a plane SHwave with a unit magnitude () propagating in the positive horizontal direction. For harmonic analysis, the term is omitted for simplicity. In the polar coordinates, the incident wave, the wave scattered by the fiber, and the wave refracted into the fiber can be expressed as series [22] as follows: where is the polar coordinates with origin at the center of the fiber; and are the wave numbers in the matrix and the fiber, respectively; and are the Bessel and Hankel function of the first kind; is the imaginary unit; and are the unknown coefficients, which are to be determined by the conditions of displacements and stresses at the interface . In this case, the conditions at the interface can be listed as follows:where and are the axial shear stresses at the interface; the superscript “+” means values in the matrix and “−” denotes those in the fiber; variable is the stiffness of the interface adopted in the spring model. Substituting of the expression of displacements in (1) and the corresponding expression for the axial shear stress into (2), the conditions at the interface described in (2) can be further expressed in detail as follows:with . Here, and are the shear modulus of the matrix and the fiber. The symbol Prime means the firstorder derivative. Unlike the pure perfect and imperfect bonding cases, the orthogonality of trigonometric functions has not been available to (3). Therefore, analytical solutions for and cannot be obtained as the pure perfect and imperfect bonding cases. In this work, (3) has been solved using the CP method [39].
2.2. The CP Method to Solve (3)
Previous works [39, 40] have shown that several terms of the lowest order are sufficient to obtain convergent results for the infinite series expressed in (1), for cases in which the fiber is purely perfectly bonded to the matrix. In this work, those infinite series in (1) are truncated with the maximum order . Then the number of unknown in (3) equals . In order to solve these unknowns, (3) is forced to be satisfied at CPs. In this work, these CPs are distributed symmetrically about the partial imperfectly and perfectly bonded interfaces, respectively, as seen in Figure 2. Additionally, these CPs are required to be distributed as uniform as possible along the whole interface. Then a linear system of size is established for aswith and as the collocation point is located within the imperfect bonding; otherwise, and . The angle is illustrated in Figure 2. Solving (4) produces all and , . Then displacement and axial shear stress at each point in both the matrix and the fiber can be obtained.
Additionally, the two parameters, scattering cross section () [22], and the directional farfield scattering magnitude [20], required by those theoretical models (including the independent scattering model) to evaluate attenuation and phase velocity of coherent waves in composites, can be also calculated. is defined as , where the numerator is the averaged (over one period) energy flux of the scattered wave across a closed circle with a radius enclosing the fiber. It is expressed asHere, an asterisk indicates the complex conjugate of the corresponding quantity. Substituting the displacements of the scattered wave in (1) and the corresponding axial shear stresses into (5), In the derivation of the above equation, the orthogonality of trigonometric functions has been used. As , and . Therefore, the above equation becomesas . The denominator is the corresponding energy flux of the incident wave across a unit length, which has an expression . Then, , which is a function of and . is an evaluation of the energy scattered by the fiber compared with the incident wave, which reflects the scattering capability of the fiber.
The directional farfield scattering magnitude is defined asas . Substituting the displacements of the scattered wave in (1), equals , which also depends on and . evaluates the scattering magnitude of the fiber in direction . A larger means the magnitude of the scattered wave in the direction is stronger.
Next, the process of this CP method is validated using several examples. Firstly, the single scattering by a Sic fiber embedded in the Aluminum with a pure imperfectly bonded interface () is calculated. The stiffness of the imperfectly bonded interface is set to be , where . The radius of the Sic fiber is chosen to be 71 m except if stated otherwise. Material properties of the fiber and the matrix are listed in Table 1. The comparison of the calculated between the exact solutions (see the Appendix) and the current CP method at different frequencies is shown in Figure 3. For the CP method, two different s have been considered. It is readily seen from this figure that results under these two different s agree well with each other, only a little bit of difference exists at high frequencies, and both numerical results agree well with the exact solutions. Secondly, two single scattering problems with partial imperfectly bonded interfaces are calculated. The results are presented in Figure 4. It is clear that the results converge to the exact solutions as increases. Additionally, comparison of the results in Figures 3 and 4 indicates that more terms in the infinite series in (1) are required for cases with partial imperfectly bonded interfaces than those with a pure imperfect and perfect bonding. The main reason may be caused by the discontinuity of displacements at the two junctions between the partial imperfect bonding and the remaining perfect bonding. In this work, is set to be 10 in all the following calculations.

(a)
(b)
Lastly, the single scattering by a fiber with a partial debonding from the matrix is calculated. In order to model the debonding (crack), the stiffness is set to be a small value in the current CP method. The normalized frequency is , and the crack is located at and . The crack opening displacement (COD) is calculated and plotted in Figure 5(a) and the farfield scattering magnitude pattern is illustrated in Figure 5(b). For a comparison, the results obtained from the Chebyshev expansion method [36–38] are also illustrated in Figure 5. The results obtained from the current CP method are found to be in good agreement with those of the Chebyshev expansion method. All these examples show that the current CP method has a good performance.
(a)
(b)
3. Parameter Study and Analysis
In this section, the effects of the normalized frequency of the incident wave, the stiffness and width, , of the imperfectly bonded interface on the farfield scattering magnitude, scattering cross section, and scattering attenuation coefficient of coherent SHwaves in fiberreinforced composites are investigated.
3.1. Influence of on the FarField Scattering Magnitude and
In this subsection, the influence of on the farfield scattering magnitude and is investigated. The material properties of the fiber and matrix are the same as those used before. The normalized frequency of the incident wave is . The partial imperfectly bonded interface has a stiffness and is located at . The farfield scattering magnitudes under several different are shown in Figure 6(a). These results show that the forward scattering (along the direction of the incident wave) decreases and the backward scattering (opposite to the direction of the incident wave) increases with the increment of . changing with is plotted in Figure 6(b). It is clear that, in general, increment of reduces the scattering capacity of the fiber.
(a)
(b)
3.2. Influence of Stiffness on
The changing with stiffness under several is plotted in Figure 7. In this figure, the partial imperfectly bonded interface is also located at . As evident from this figure, s under different approach the same value as increases to a large value. The reason is that all the imperfect bonding will become the perfect one as increases. under each approaches a constant as decreases, since the imperfect bonding will become the debonding one with the decrement of . Moreover, decreases as the imperfectly bonded interface deteriorates. With the increment of , becomes more sensitive to the stiffness , which indicates that the ultrasonic technique is more suitable to evaluate large imperfectly bonded interfaces.
3.3. Influence of and on the Scattering Attenuation of Coherent SHWaves in FiberReinforced Composites
In this work, we assume that the orientation of the partial imperfectly bonded interface is randomly distributed, and the width of the imperfect interface is the same for all the fibers. According to the independent scattering model, the attenuation coefficient of coherent SHwaves in such a composite can be calculated aswhere is the statistical average of over angle , with the definition as follows:where is the probability density function of . In this work, we assume that has a uniform distribution over . In the calculations, s at 100 different uniformly selected within are calculated and averaged to obtain . The contour plots of the calculated changing with and under three different s are shown in Figure 8. The results indicate that and the scattering attenuation coefficient become more sensitive to , as the normalized frequency of the incident wave is within about . The sensitivity of and to is low for both much lower and much higher frequencies of the incident wave. Accordingly, the frequency of waves should be carefully selected for the using of the ultrasonic technique to nondestructively evaluate the interfaces in composites.
(a)
(b)
(c)
3.4. Influence of and on the Scattering Attenuation
In this subsection, the influence of and the stiffness on the scattering attenuation of coherent SHwaves in fiberreinforced composites is investigated. As before, the orientation of the partial imperfectly bonded interface is randomly distributed. The contour plots of the calculated changing with and at several different are shown in Figure 9. The results show that and the scattering attenuation are not sensitive to the initial deterioration of interfaces. As the interfaces deteriorate further, the sensitivity of and scattering attenuation to increases. This fact benefits the nondestructive evaluation of interfaces in composites using the scattering attenuation based ultrasonic technique. Similar to the conclusion reached in Section 3.3, the sensitivity of and scattering attenuation to is affected by the frequency of the incident wave. In this case, and the scattering attenuation become more sensitive to as the normalized frequency of the incident wave is within about.
(a)
(b)
(c)
(d)
(e)
(f)
4. Conclusions
The influence of the partial imperfectly bonded interfaces between each fiber and the matrix on the farfield scattering magnitude, scattering cross section, and scattering attenuation of coherent SHwaves in fiberreinforced composites is investigated. The imperfection of interfaces is modelled using the spring model, and the CP method is proposed to solve the single scattering problem with such a partial imperfectly bonded interface. Based on the analysis of the single scattering problems, the effects of frequency of the incident wave and parameters of the partial imperfectly bonded interfaces on the farfield scattering magnitude, scattering cross section, and the scattering attenuation of coherent SHwaves in fiberreinforced composites are investigated numerically. The main contribution of this work includes the following:(i)The CP method is proposed to solve the single scattering of SHwaves with a partial imperfectly bonded interface between the fiber and the matrix. Compared with other existing techniques, the derivation process and implementation of the current CP method are simple and it has a comparable accuracy. The disadvantage is that it cannot give accurate solutions at points near the two tips of the imperfect interface. In general, the current CP method can be taken as an auxiliary one for the existing methods.(ii)Extensive parameter study shows that the sensitivity of scattering attenuation of coherent SHwaves in fiberreinforced composites to the imperfection/deterioration of interfaces is significantly affected by the frequency of the incident wave. Therefore, waves should be carefully selected for the use of the scattering attenuation based ultrasonic technique to nondestructively evaluate the interfaces in composites.(iii)Initial deterioration of the interfaces is relatively difficult to be detected using the scattering attenuation based ultrasonic technique.
Although this work only considers the SHwaves, the analysis framework can be extended to more complicated cases, such as the P and SVwaves, without significant difficulties.
Appendix
Exact Solutions of the Single Scattering Problem with a Pure Imperfectly Bonded Interface
Figure 10 shows the schematic of the single scattering of a plane SHwave with a unit magnitude by a fiber with a pure imperfectly bonded interface to the matrix. is the scattered wave and is the refracted wave into the fiber. The incident, scattered, and refracted waves can be expressed in series as those in (1).
In this case, the two boundary conditions of displacements and stresses are Since varies from 0 to in the above equation, the orthogonality of trigonometric functions is available. Then substitution of the displacements in (1) and the corresponding expressions for the stresses into (A.1) yields the analytical solutions for and :
Disclosure
Part of this work has been presented orally at “The 2nd International Conference on Computational Engineering and Science for Safety and Environmental Problems,” 2017.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This work was supported by the Chinese National Natural Science Fund (Grants 11502036, 11372104, and 11632004) and the Natural Science Fund of the city of Chongqing (Grant cstc2015jcyjA0577). The Key Program for International Science and Technology Cooperation Projects of the Ministry of Science and Technology of China (no. 2016YFE0125900) and the Guangxi Key Laboratory of Manufacturing Systems and Advanced Manufacturing Technology (no. 1638012014k) also provided partial financial support.
References
 S. I. Rokhlin, W. Huang, and Y. C. Chu, “Ultrasonic scattering and velocity methods for characterization of fibrematrix interphases,” Ultrasonics, vol. 33, no. 5, pp. 351–364, 1995. View at: Publisher Site  Google Scholar
 R.B. Yang and A. K. Mal, “Elastic waves in a composite containing inhomogeneous fibers,” International Journal of Engineering Science, vol. 34, no. 1, pp. 67–79, 1996. View at: Publisher Site  Google Scholar
 Y.S. Wang and D. Wang, “Scattering of elastic waves by a rigid cylindrical inclusion partially debonded from its surrounding matrix  I. SH case,” International Journal of Solids and Structures, vol. 33, no. 19, pp. 2789–2815, 1996. View at: Publisher Site  Google Scholar
 X. Ji, X. Luo, Y. Q. Yang, J. H. Lou, B. Huang, and X. R. Wang, “Research progress of nondestructive testing for continuous fiberreinforced metalmatrix composites,” Rare Metal Materials and Engineering, vol. 42, pp. 401–405, 2013. View at: Google Scholar
 L. L. Shaw, P. Karpur, and T. E. Matikas, “Fracture strength and damage progression of the fiber/matrix interfaces in titaniumbased MMCs with different interfacial layers,” Composites Part B: Engineering, vol. 29, no. 3, pp. 331–339, 1998. View at: Publisher Site  Google Scholar
 Y.S. Wang and D. Wang, “Scattering of elastic waves by a rigid cylindrical inclusion partially debonded from its surrounding matrix  II. P and SV cases,” International Journal of Solids and Structures, vol. 33, no. 19, pp. 2817–2840, 1996. View at: Publisher Site  Google Scholar
 T. H. Ju and S. K. Datta, “Dynamics of a laminated composite plate with interface layers,” Journal of Nondestructive Evaluation, vol. 11, no. 34, pp. 227–235, 1992. View at: Publisher Site  Google Scholar
 L. L. Foldy, “The multiple scattering of waves. I. General theory of isotropic scattering by randomly distributed scatterers,” Physical Review A: Atomic, Molecular and Optical Physics, vol. 67, no. 34, pp. 107–119, 1945. View at: Publisher Site  Google Scholar
 S. Biwa, S. Idekoba, and N. Ohno, “Wave attenuation in particulate polymer composites: Independent scattering/absorption analysis and comparison to measurements,” Mechanics of Materials, vol. 34, no. 10, pp. 671–682, 2002. View at: Publisher Site  Google Scholar
 J. Zhang and W. Ye, “An explicit formula for the coherent SH waves' attenuation coefficient in random porous materials with low porosities,” Ultrasonics, vol. 62, pp. 27–34, 2015. View at: Publisher Site  Google Scholar
 S. K. Kanaun and V. M. Levin, “Effective medium method in the problem of axial elastic shear wave propagation through fiber composites,” International Journal of Solids and Structures, vol. 40, no. 18, pp. 4859–4878, 2003. View at: Publisher Site  Google Scholar
 P. G. Bussink, P. L. Iske, J. Oortwijn, and G. L. Verbist, “Selfconsistent analysis of elastic wave propagation in twodimensional matrixinclusion composites,” Journal of the Mechanics and Physics of Solids, vol. 43, no. 10, pp. 1673–1690, 1995. View at: Publisher Site  Google Scholar  MathSciNet
 J.Y. Kim, J.G. Ih, and B.H. Lee, “Dispersion of elastic waves in random particulate composites,” The Journal of the Acoustical Society of America, vol. 97, no. 3, pp. 1380–1388, 1995. View at: Publisher Site  Google Scholar
 S. K. Bose and A. K. Mal, “Axial shear waves in a medium with randomly distributed cylinders,” The Journal of the Acoustical Society of America, vol. 55, no. 3, pp. 519–523, 1974. View at: Publisher Site  Google Scholar
 V. K. Varadan, Y. Ma, and V. V. Varadan, “A multiple scattering theory for elastic wave propagation in discrete random media,” The Journal of the Acoustical Society of America, vol. 77, no. 2, pp. 375–385, 1985. View at: Publisher Site  Google Scholar
 V. K. Varadan, V. V. Varadan, and Y. . Pao, “Multiple scattering of elastic waves by cylinders of arbitrary cross section. I. {SH} waves,” The Journal of the Acoustical Society of America, vol. 63, no. 5, pp. 1310–1319, 1978. View at: Publisher Site  Google Scholar  MathSciNet
 A. I. Beltzer, “The effective dynamic response of random composites and polycrystalsa survey of the causal approach,” Wave Motion, vol. 11, no. 3, pp. 211–229, 1989. View at: Publisher Site  Google Scholar
 A. I. Beltzer and N. Brauner, “The dynamic response of random composites by a causal differential method,” Mechanics of Materials, vol. 6, no. 4, pp. 337–345, 1987. View at: Publisher Site  Google Scholar
 S. Biwa, “Independent scattering and wave attenuation in viscoelastic composites,” Mechanics of Materials, vol. 33, no. 11, pp. 635–647, 2001. View at: Publisher Site  Google Scholar
 A. N. Norris and J.M. Conoir, “Multiple scattering by cylinders immersed in fluid: High order approximations for the effective wavenumbers,” The Journal of the Acoustical Society of America, vol. 129, no. 1, pp. 104–113, 2011. View at: Publisher Site  Google Scholar
 N. Brauner and A. I. Beltzer, “Waveobstacle interaction in a lossy medium: energy perturbations and negative extinction,” Ultrasonics, vol. 26, no. 6, pp. 328–334, 1988. View at: Publisher Site  Google Scholar
 C. C. Mow and Y. H. Pao, The diffraction of elastic waves and dynamic stress concentrations, 1971.
 N. M. Newmark, “Tests and analysis of composite beams with incomplete interaction,” in Proceedings of the Society for Experimental Stress Analysis, vol. 9, pp. 75–92, 1951. View at: Google Scholar
 J.M. Baik and R. B. Thompson, “Ultrasonic scattering from imperfect interfaces: A quasistatic model,” Journal of Nondestructive Evaluation, vol. 4, no. 34, pp. 177–196, 1984. View at: Publisher Site  Google Scholar
 W. Huang, S. Brisuda, and S. I. Rokhlin, “Ultrasonic wave scattering from fibermatrix interphases,” The Journal of the Acoustical Society of America, vol. 97, no. 2, pp. 807–817, 1995. View at: Publisher Site  Google Scholar
 P. A. Martin, “Boundary integral equations for the scattering of elastic waves by elastic inclusions with thin interface layers,” Journal of Nondestructive Evaluation, vol. 11, no. 34, pp. 167–174, 1992. View at: Publisher Site  Google Scholar
 J. A. Otero, R. RodriguezRamos, G. Monsivais, C. Stern, R. Martínez, and R. Dario, “Interfacial waves between two magnetoelectroelastic halfspaces with magnetoelectromechanical imperfect interface,” Philosophical Magazine Letters, vol. 94, no. 10, pp. 629–638, 2014. View at: Publisher Site  Google Scholar
 H. BritoSantana, R. de Medeiros, R. RodriguezRamos, and V. Tita, “Different interface models for calculating the effective properties in piezoelectric composite materials with imperfect fiber–matrix adhesion,” Composite Structures, vol. 151, pp. 70–80, 2016. View at: Publisher Site  Google Scholar
 J. Aboudi, “Damage in compositesModeling of imperfect bonding,” Composites Science and Technology, vol. 28, no. 2, pp. 103–128, 1987. View at: Publisher Site  Google Scholar
 A. K. Mal and S. K. Bose, “Dynamic elastic moduli of a suspension of imperfectly bonded spheres,” in Proceedings of the, vol. 76, pp. 587–600. View at: Google Scholar
 J. C. LopezRealpozo, R. RodriguezRamos, R. GuinovartDiaz et al., “Effective elastic shear stiffness of a periodic fibrous composite with nonuniform imperfect contact between the matrix and the fibers,” International Journal of Solids and Structures, vol. 51, no. 6, pp. 1253–1262, 2014. View at: Publisher Site  Google Scholar
 M. Wurkner, H. Berger, and U. Gabbert, “Numerical investigations of effective properties of fiber reinforced composites with parallelogram arrangements and imperfect interface,” Composite Structures, vol. 116, no. 1, pp. 388–394, 2014. View at: Publisher Site  Google Scholar
 E. Vilchevskaya and I. Sevostianov, “Overall thermal conductivity of a fiber reinforced composite with partially debonded inhomogeneities,” International Journal of Engineering Science, vol. 98, pp. 99–109, 2016. View at: Publisher Site  Google Scholar
 S. Biwa and T. Shibata, “Elastic and ultrasonic properties of a unidirectional composite with partially debonded fibres: Numerical analysis for longitudinal shear modes,” Composites Science and Technology, vol. 60, no. 1, pp. 83–93, 2000. View at: Publisher Site  Google Scholar
 H. Sato and Y. Shindo, “Influence of microstructure on scattering of plane elastic waves by a distribution of partially debonded elliptical inclusions,” Mechanics of Materials, vol. 34, no. 7, pp. 401–409, 2002. View at: Publisher Site  Google Scholar
 Y. Yang and A. Norris, “Longitudinal wave scattering from a partially bonded fiber,” Wave Motion, vol. 15, no. 1, pp. 43–59, 1992. View at: Publisher Site  Google Scholar
 J.Y. Kim, “Attenuation and speed of antiplane shear wave in fiberreinforced composites with random interfacial cracks,” International Journal of Solids and Structures, vol. 38, no. 4041, pp. 7121–7137, 2001. View at: Publisher Site  Google Scholar
 W. Liu and R. D. Kriz, “Shear waves in fiberreinforced composites with interfacial cracks,” International Journal of Solids and Structures, vol. 35, no. 13, pp. 1425–1449, 1998. View at: Publisher Site  Google Scholar
 S. Biwa, S. Yamamoto, F. Kobayashi, and N. Ohno, “Computational multiple scattering analysis for shear wave propagation in unidirectional composites,” International Journal of Solids and Structures, vol. 41, no. 2, pp. 435–457, 2004. View at: Publisher Site  Google Scholar
 J. Zhang, W. Ye, and T. X. Yu, “Numerical simulation of effective phase velocity and attenuation of shear elastic wave propagation in unidirectional composite materials,” Ultrasonics, vol. 53, no. 6, pp. 1200–1211, 2013. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2018 Jun Zhang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.