# Scattering Attenuation of SH-Waves in Fiber-Reinforced Composites with Partial Imperfectly Bonded Interfaces

Authors: Jun Zhang; Longhai Zeng; Lian Liu
Journal: Shock and Vibration
 (2018)

---

## Abstract

The influence of partial imperfectly bonded interfaces between each fiber and the matrix on the scattering attenuation of coherent SH-waves in fiber-reinforced composites is investigated. The imperfection of interfaces is modelled using the spring model, in which the level of imperfection is characterized by a parameterK, which is called the stiffness of the imperfect interfaces. First, the single scattering of SH-waves 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 far-field scattering magnitude, scattering cross section, and scattering attenuation coefficients of coherent SH-waves are shown graphically. The potential application of the current results to nondestructive evaluation of interfaces in composites is also discussed.


---

## Body

## 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 1/2σsn0, where σs is the scattering cross section [22] and n0 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 spring-mass 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 σr=a=Ku, where u denotes the jump of the displacements across the interface and K is the parameter to characterize the imperfection level of the interfaces. It is clear that, when K→∞, the spring model becomes the familiar perfect bonding and it is degenerated to the debonding case as K→0. Other intermediate values of K 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 spring-mass 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 SV-waves in fiber-reinforced composites. Mal and Bose [30] have studied the phase velocity and attenuation of coherent P- and SV-waves 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 toenter 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, Lopez-Realpozo 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 SH-waves in fiber-reinforced 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 SH-waves 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 SH-waves 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 far-field scattering magnitude, scattering cross section, and scattering attenuation is then conducted. At the end, a short conclusion is presented.


## 2. Single Scattering of SH-Waves by a Fiber with a Partial Imperfect Bonding to the Matrix

### 2.1. Formulation of the Problem

Figure1 shows a schematic of the single scattering of a plane SH-wave by a cylindrical fiber with radius a. 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 β=0, the whole fiber is purely perfectly bonded to the matrix and the whole interface becomes a pure imperfect bonding as β=2π. The incident wave is a plane SH-wave with a unit magnitude (uincx=eikmx-ωt) propagating in the positive horizontal direction. For harmonic analysis, the term e-iωt 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: (1)uincr,θ=∑n=-∞∞inJnkmreinθ,r≥a,usr,θ=∑n=-∞∞AnHn1kmreinθ,r≥a,ufr,θ=∑n=-∞∞BnJnkfreinθ,r≤a,where r,θ is the polar coordinates with origin at the center of the fiber; km and kf are the wave numbers in the matrix and the fiber, respectively; Jn(·) and Hn(1)(·) are the Bessel and Hankel function of the first kind; i=-1 is the imaginary unit; An and Bn are the unknown coefficients, which are to be determined by the conditions of displacements and stresses at the interface r=a. In this case, the conditions at the interface can be listed as follows:(2)σrz+r=a=σrz-r=a,θ∈02πu+r=a-u-r=a=1Kσrz+r=a,θ∈Imperfect  Bonding0,θ∈Perfect  Bonding,where σrz+r=a and σrz-r=a are the axial shear stresses at the interface; the superscript “+” means values in the matrix and “−” denotes those in the fiber; variable K 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:(3)∑n=-∞∞inJnkma+AnHn1kma-γinJn′kma+AnHn1′kmaJn′kfaJnkfaeinθ=∑n=-∞∞inJn′kma+AnHn1′kmaμmkmKeinθ,θ∈Imperfect  Bonding0,θ∈Perfect  Bondingwith γ=μmkm/μfkf. Here, μm and μf are the shear modulus of the matrix and the fiber. The symbol Prime means the first-order derivative. Unlike the pure perfect and imperfect bonding cases, the orthogonality of trigonometric functions has not been available to (3). Therefore, analytical solutions for An and Bnn=-∞,…,∞ cannot be obtained as the pure perfect and imperfect bonding cases. In this work, (3) has been solved using the CP method [39].Figure 1
Scattering of a plane SH-wave by a cylindrical fiber with a partial imperfect bonding to the matrix. The range of the imperfect bonding is denoted as the red solid arc.


### 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 nmax. Then the number of unknown An in (3) equals 2nmax+1. In order to solve these unknowns, (3) is forced to be satisfied at 2nmax+1 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 2nmax+1 is established for An,n=-nmax,…,nmax as(4)ClnAn=bll,n=-nmax,…,nmaxwith Cln=Hn1kma-γHn1′kma/Jn′kfaJnkfa-μmkm/KHn1′kmaeinθl and bl=inμmkm/KJn′kma-Jnkma+γJn′kma/Jn′kfaJnkfaeinθl as the collocation point is located within the imperfect bonding; otherwise, Cln=Hn1kma-γHn1′kma/Jn′kfaJnkfaeinθl and bl=-inJnkmaeinθl. The angle θl is illustrated in Figure 2. Solving (4) produces all An and Bn, n=-nmax,…,nmax. Then displacement and axial shear stress at each point in both the matrix and the fiber can be obtained.Figure 2
Distribution of CPs denoted by “×” symbol. It is noted that at least one CP should be distributed on the imperfectly bonded interface as long as β is not equal to 0. In this figure, =π/4β=π/2, and nmax = 6.Additionally, the two parameters, scattering cross section (σs) [22], and the directional far-field scattering magnitude fθ [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. σs is defined as σs=AveE˙/Avee˙, where the numerator is the averaged (over one period) energy flux of the scattered wave across a closed circle with a radius R enclosing the fiber. It is expressed as(5)AveE˙=-iω4∫02πσrzs×us∗-σrzs∗×usr=RRdθ.Here, 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), (6)AveE˙=-iπω2μmkmR×∑n=-nmaxnmaxAn2Hn1′kmRHn1∗kmR-Hn1′∗kmRHn1kmR.In the derivation of the above equation, the orthogonality of trigonometric functions has been used. As kmR→∞, Hn1kmR→2/πkmReikmR-n+1/2π/2 and Hn1′kmR→iHn1kmR. Therefore, the above equation becomes(7)AveE˙=2ωμm∑n=-nmaxnmaxAn2as kmR→∞. The denominator is the corresponding energy flux of the incident wave across a unit length, which has an expression Avee˙=ω/2μmkm. Then, σsα,β=4/km∑n=-nmaxn=nmaxAn2, which is a function of α and β. σs 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 far-field scattering magnitude is defined as(8)usr,θ≈fθ2πkmreikmr-π/4as kmr→∞. Substituting the displacements of the scattered wave in (1), fθ equals ∑n=-nmaxn=nmaxAn-ineinθ, which also depends on α and β. fθ evaluates the scattering magnitude of the fiber in direction θ. A larger f(θ) 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 (β=2π) is calculated. The stiffness of the imperfectly bonded interface is set to be K=0.01Kp, where Kp=μm+μf/0.1a. The radius of the Sic fiber is chosen to be a= 71 μm except if stated otherwise. Material properties of the fiber and the matrix are listed in Table 1. The comparison of the calculated σs 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 nmaxs have been considered. It is readily seen from this figure that results under these two different nmaxs 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 nmax 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, nmax is set to be 10 in all the following calculations.Table 1
Material properties of Sic and Al.
Sic
ρ
f (kg/m3)
μ
f (GPa)
3181
188.1
Al
ρ
m (kg/m3)
μ
m (GPa)
2706
26.7Figure 3
Comparison ofσs between the current CP method and exact solutions for the pure imperfect bonding case.Figure 4
Convergence analysis of the CP method through two examples with a partial imperfectly bonded interface: (a)α=π, β=π and (b) α=π/4, β=π/2.
(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 stiffnessK is set to be a small value 10-6Kp in the current CP method. The normalized frequency is kma=3.0, and the crack is located at α=0 and β=π. The crack opening displacement (COD) is calculated and plotted in Figure 5(a) and the far-field 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.Figure 5
The calculated COD (a) and polar plot of the far-field scattering magnitude pattern (b) for the single scattering by a fiber with a partial debonding from the matrix.
(a)
(b)


## 2.1. Formulation of the Problem

Figure1 shows a schematic of the single scattering of a plane SH-wave by a cylindrical fiber with radius a. 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 β=0, the whole fiber is purely perfectly bonded to the matrix and the whole interface becomes a pure imperfect bonding as β=2π. The incident wave is a plane SH-wave with a unit magnitude (uincx=eikmx-ωt) propagating in the positive horizontal direction. For harmonic analysis, the term e-iωt 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: (1)uincr,θ=∑n=-∞∞inJnkmreinθ,r≥a,usr,θ=∑n=-∞∞AnHn1kmreinθ,r≥a,ufr,θ=∑n=-∞∞BnJnkfreinθ,r≤a,where r,θ is the polar coordinates with origin at the center of the fiber; km and kf are the wave numbers in the matrix and the fiber, respectively; Jn(·) and Hn(1)(·) are the Bessel and Hankel function of the first kind; i=-1 is the imaginary unit; An and Bn are the unknown coefficients, which are to be determined by the conditions of displacements and stresses at the interface r=a. In this case, the conditions at the interface can be listed as follows:(2)σrz+r=a=σrz-r=a,θ∈02πu+r=a-u-r=a=1Kσrz+r=a,θ∈Imperfect  Bonding0,θ∈Perfect  Bonding,where σrz+r=a and σrz-r=a are the axial shear stresses at the interface; the superscript “+” means values in the matrix and “−” denotes those in the fiber; variable K 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:(3)∑n=-∞∞inJnkma+AnHn1kma-γinJn′kma+AnHn1′kmaJn′kfaJnkfaeinθ=∑n=-∞∞inJn′kma+AnHn1′kmaμmkmKeinθ,θ∈Imperfect  Bonding0,θ∈Perfect  Bondingwith γ=μmkm/μfkf. Here, μm and μf are the shear modulus of the matrix and the fiber. The symbol Prime means the first-order derivative. Unlike the pure perfect and imperfect bonding cases, the orthogonality of trigonometric functions has not been available to (3). Therefore, analytical solutions for An and Bnn=-∞,…,∞ cannot be obtained as the pure perfect and imperfect bonding cases. In this work, (3) has been solved using the CP method [39].Figure 1
Scattering of a plane SH-wave by a cylindrical fiber with a partial imperfect bonding to the matrix. The range of the imperfect bonding is denoted as the red solid arc.


## 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 nmax. Then the number of unknown An in (3) equals 2nmax+1. In order to solve these unknowns, (3) is forced to be satisfied at 2nmax+1 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 2nmax+1 is established for An,n=-nmax,…,nmax as(4)ClnAn=bll,n=-nmax,…,nmaxwith Cln=Hn1kma-γHn1′kma/Jn′kfaJnkfa-μmkm/KHn1′kmaeinθl and bl=inμmkm/KJn′kma-Jnkma+γJn′kma/Jn′kfaJnkfaeinθl as the collocation point is located within the imperfect bonding; otherwise, Cln=Hn1kma-γHn1′kma/Jn′kfaJnkfaeinθl and bl=-inJnkmaeinθl. The angle θl is illustrated in Figure 2. Solving (4) produces all An and Bn, n=-nmax,…,nmax. Then displacement and axial shear stress at each point in both the matrix and the fiber can be obtained.Figure 2
Distribution of CPs denoted by “×” symbol. It is noted that at least one CP should be distributed on the imperfectly bonded interface as long as β is not equal to 0. In this figure, =π/4β=π/2, and nmax = 6.Additionally, the two parameters, scattering cross section (σs) [22], and the directional far-field scattering magnitude fθ [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. σs is defined as σs=AveE˙/Avee˙, where the numerator is the averaged (over one period) energy flux of the scattered wave across a closed circle with a radius R enclosing the fiber. It is expressed as(5)AveE˙=-iω4∫02πσrzs×us∗-σrzs∗×usr=RRdθ.Here, 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), (6)AveE˙=-iπω2μmkmR×∑n=-nmaxnmaxAn2Hn1′kmRHn1∗kmR-Hn1′∗kmRHn1kmR.In the derivation of the above equation, the orthogonality of trigonometric functions has been used. As kmR→∞, Hn1kmR→2/πkmReikmR-n+1/2π/2 and Hn1′kmR→iHn1kmR. Therefore, the above equation becomes(7)AveE˙=2ωμm∑n=-nmaxnmaxAn2as kmR→∞. The denominator is the corresponding energy flux of the incident wave across a unit length, which has an expression Avee˙=ω/2μmkm. Then, σsα,β=4/km∑n=-nmaxn=nmaxAn2, which is a function of α and β. σs 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 far-field scattering magnitude is defined as(8)usr,θ≈fθ2πkmreikmr-π/4as kmr→∞. Substituting the displacements of the scattered wave in (1), fθ equals ∑n=-nmaxn=nmaxAn-ineinθ, which also depends on α and β. fθ evaluates the scattering magnitude of the fiber in direction θ. A larger f(θ) 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 (β=2π) is calculated. The stiffness of the imperfectly bonded interface is set to be K=0.01Kp, where Kp=μm+μf/0.1a. The radius of the Sic fiber is chosen to be a= 71 μm except if stated otherwise. Material properties of the fiber and the matrix are listed in Table 1. The comparison of the calculated σs 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 nmaxs have been considered. It is readily seen from this figure that results under these two different nmaxs 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 nmax 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, nmax is set to be 10 in all the following calculations.Table 1
Material properties of Sic and Al.
Sic
ρ
f (kg/m3)
μ
f (GPa)
3181
188.1
Al
ρ
m (kg/m3)
μ
m (GPa)
2706
26.7Figure 3
Comparison ofσs between the current CP method and exact solutions for the pure imperfect bonding case.Figure 4
Convergence analysis of the CP method through two examples with a partial imperfectly bonded interface: (a)α=π, β=π and (b) α=π/4, β=π/2.
(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 stiffnessK is set to be a small value 10-6Kp in the current CP method. The normalized frequency is kma=3.0, and the crack is located at α=0 and β=π. The crack opening displacement (COD) is calculated and plotted in Figure 5(a) and the far-field 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.Figure 5
The calculated COD (a) and polar plot of the far-field scattering magnitude pattern (b) for the single scattering by a fiber with a partial debonding from the matrix.
(a)
(b)


## 3. Parameter Study and Analysis

In this section, the effects of the normalized frequencykma of the incident wave, the stiffness K and width, β, of the imperfectly bonded interface on the far-field scattering magnitude, scattering cross section, and scattering attenuation coefficient of coherent SH-waves in fiber-reinforced composites are investigated.


### 3.1. Influence ofβ on the Far-Field Scattering Magnitude and σs

In this subsection, the influence ofβ on the far-field scattering magnitude and σs 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 kma=3.0. The partial imperfectly bonded interface has a stiffness K=0.01Kp and is located at α=0.0. The far-field scattering magnitudes under several different βsπ/6,π/3,π/2andπ 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 β. σs changing with β is plotted in Figure 6(b). It is clear that, in general, increment of β reduces the scattering capacity of the fiber.Figure 6
Influence ofβ on (a) the far-field scattering magnitude and (b) σs at kma=3.0 and α=0.0.
(a)
(b)


### 3.2. Influence of StiffnessK on σs

Theσs changing with stiffness K under several βs is plotted in Figure 7. In this figure, the partial imperfectly bonded interface is also located at α=0.0. As evident from this figure, σss under different βs approach the same value as K increases to a large value. The reason is that all the imperfect bonding will become the perfect one as K increases. σs under each β approaches a constant as K decreases, since the imperfect bonding will become the debonding one with the decrement of K. Moreover, σs decreases as the imperfectly bonded interface deteriorates. With the increment of β, σs becomes more sensitive to the stiffness K, which indicates that the ultrasonic technique is more suitable to evaluate large imperfectly bonded interfaces.Figure 7
Influence of stiffnessK on the σs at kma=3.0 under several βs.


### 3.3. Influence ofkma and β on the Scattering Attenuation of Coherent SH-Waves in Fiber-Reinforced 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 SH-waves in such a composite can be calculated as(9)αs=12σsn0,where σs is the statistical average of σs over angle α, with the definition as follows:(10)σs=∫02πpασsα,βdα,where p(α) is the probability density function of α. In this work, we assume that α has a uniform distribution over 02π. In the calculations, σss at 100 different αs uniformly selected within 02π are calculated and averaged to obtain σs. The contour plots of the calculated σs changing with kma and β under three different Ks are shown in Figure 8. The results indicate that σs and the scattering attenuation coefficient αs become more sensitive to β, as the normalized frequency of the incident wave is within about kma∈3.04.0. The sensitivity of σs and αs 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.Figure 8
Contour plots ofσs changing with kma and β at (a) K/Kp=0.1, (b) K/Kp=0.01, and (c) K/Kp=0.001.
(a)
(b)
(c)


### 3.4. Influence ofkma and K on the Scattering Attenuation

In this subsection, the influence ofkma and the stiffness K on the scattering attenuation of coherent SH-waves in fiber-reinforced composites is investigated. As before, the orientation of the partial imperfectly bonded interface is randomly distributed. The contour plots of the calculated σs changing with kma and K at several different βs are shown in Figure 9. The results show that σs and the scattering attenuation are not sensitive to the initial deterioration of interfaces. As the interfaces deteriorate further, the sensitivity of σs and scattering attenuation to K 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 σs and scattering attenuation to K is affected by the frequency of the incident wave. In this case, σs and the scattering attenuation become more sensitive to K as the normalized frequency of the incident wave is within aboutkma∈3.04.0.Figure 9
Contour plots ofσs changing with kma and K at several βs: (a) β=π/6, (b) β=π/3, (c) β=π/2, (d) β=π, (e) β=3/2π, and (f) β=2π.
(a)
(b)
(c)
(d)
(e)
(f)


## 3.1. Influence ofβ on the Far-Field Scattering Magnitude and σs

In this subsection, the influence ofβ on the far-field scattering magnitude and σs 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 kma=3.0. The partial imperfectly bonded interface has a stiffness K=0.01Kp and is located at α=0.0. The far-field scattering magnitudes under several different βsπ/6,π/3,π/2andπ 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 β. σs changing with β is plotted in Figure 6(b). It is clear that, in general, increment of β reduces the scattering capacity of the fiber.Figure 6
Influence ofβ on (a) the far-field scattering magnitude and (b) σs at kma=3.0 and α=0.0.
(a)
(b)


## 3.2. Influence of StiffnessK on σs

Theσs changing with stiffness K under several βs is plotted in Figure 7. In this figure, the partial imperfectly bonded interface is also located at α=0.0. As evident from this figure, σss under different βs approach the same value as K increases to a large value. The reason is that all the imperfect bonding will become the perfect one as K increases. σs under each β approaches a constant as K decreases, since the imperfect bonding will become the debonding one with the decrement of K. Moreover, σs decreases as the imperfectly bonded interface deteriorates. With the increment of β, σs becomes more sensitive to the stiffness K, which indicates that the ultrasonic technique is more suitable to evaluate large imperfectly bonded interfaces.Figure 7
Influence of stiffnessK on the σs at kma=3.0 under several βs.


## 3.3. Influence ofkma and β on the Scattering Attenuation of Coherent SH-Waves in Fiber-Reinforced 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 SH-waves in such a composite can be calculated as(9)αs=12σsn0,where σs is the statistical average of σs over angle α, with the definition as follows:(10)σs=∫02πpασsα,βdα,where p(α) is the probability density function of α. In this work, we assume that α has a uniform distribution over 02π. In the calculations, σss at 100 different αs uniformly selected within 02π are calculated and averaged to obtain σs. The contour plots of the calculated σs changing with kma and β under three different Ks are shown in Figure 8. The results indicate that σs and the scattering attenuation coefficient αs become more sensitive to β, as the normalized frequency of the incident wave is within about kma∈3.04.0. The sensitivity of σs and αs 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.Figure 8
Contour plots ofσs changing with kma and β at (a) K/Kp=0.1, (b) K/Kp=0.01, and (c) K/Kp=0.001.
(a)
(b)
(c)


## 3.4. Influence ofkma and K on the Scattering Attenuation

In this subsection, the influence ofkma and the stiffness K on the scattering attenuation of coherent SH-waves in fiber-reinforced composites is investigated. As before, the orientation of the partial imperfectly bonded interface is randomly distributed. The contour plots of the calculated σs changing with kma and K at several different βs are shown in Figure 9. The results show that σs and the scattering attenuation are not sensitive to the initial deterioration of interfaces. As the interfaces deteriorate further, the sensitivity of σs and scattering attenuation to K 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 σs and scattering attenuation to K is affected by the frequency of the incident wave. In this case, σs and the scattering attenuation become more sensitive to K as the normalized frequency of the incident wave is within aboutkma∈3.04.0.Figure 9
Contour plots ofσs changing with kma and K at several βs: (a) β=π/6, (b) β=π/3, (c) β=π/2, (d) β=π, (e) β=3/2π, and (f) β=2π.
(a)
(b)
(c)
(d)
(e)
(f)


## 4. Conclusions

The influence of the partial imperfectly bonded interfaces between each fiber and the matrix on the far-field scattering magnitude, scattering cross section, and scattering attenuation of coherent SH-waves in fiber-reinforced 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 far-field scattering magnitude, scattering cross section, and the scattering attenuation of coherent SH-waves in fiber-reinforced 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 SH-waves 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 SH-waves in fiber-reinforced 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 SH-waves, the analysis framework can be extended to more complicated cases, such as the P- and SV-waves, without significant difficulties.



---

*Source: 2341797-2018-02-20.xml*