# Study on the Single Scattering of Elastic Waves by a Cylindrical Fiber with a Partially Imperfect Bonding Using the Collocation Point Method

Authors: Jun Zhang; Longhai Zeng; Chuanlin Hu; Wensheng Yan
Journal: Shock and Vibration
 (2018)

---

## Abstract

The single scattering of P- and SV-waves by a cylindrical fiber with a partially imperfect bonding to the surrounding matrix is investigated, which benefits the characterization of the behavior of elastic waves in composite materials. The imperfect interface is modelled by the spring model. To solve the corresponding single scattering problem, a collocation point (CP) method is introduced. Based on this method, influence of various aspects of the imperfect interface on the scattering of P- and SV-waves is studied. Results indicate that (i) the total scattering cross section (SCS) is almost symmetric about the axisα=π/2 with respect to the location (α) of the imperfect interface, (ii) imperfect interfaces located at α=0 and α=π highly reduce the total SCS under a P-wave incidence and imperfect interfaces located at α=π/2 reduce the total SCS most significantly under SV-incidence, and (iii) under a P-wave incidence the SCS has a high sensitivity to the bonding level of imperfect interfaces when α is small, while it becomes more sensitive to the bonding level when α is larger under SV-wave incidence.


---

## Body

## 1. Introduction

A large number of interfaces exist in composite materials and play a very important role in the performance of composite materials, for example, the transmission of the load between the matrix and fibers. To simplify the analysis and calculation, the fibers/particulates are generally assumed to be perfectly bonded to the surrounding matrix [1–4]. However, this is not always the case, which has consequently led to the increased research effort to take the bonding conditions into account [5–16]. Among them, the debonding case [9–13, 16–18], which is assumed to be a crack with noncontacting faces, has been well investigated. Specifically, Norris and Yang [11–13] investigated the single scattering of antiplane shear and P- and SV-waves by a partially debonded fiber through the expansion of the crack opening displacement (COD) or the stress at the neck in terms of Chebyshev polynomials. The results indicated that a resonance happens at low frequencies when the neck joining the fiber to the matrix is thin. Recently, Fang et al. [16] extended this method to the case of pipelines. Based on the same method, Kim [9] studied the effects of interface debonding on the attenuation and speed of antiplane shear wave propagating in fiber-reinforced composites. The corresponding static effective mechanical properties of fiber-reinforced composites were obtained from the asymptotic solutions of dynamic harmonic analysis by Liu and Kriz [10]. Besides the technique proposed by Norris and Yang [11–13], Y. S. Wang and D. Wang [17, 18] proposed a different method, where the dislocation density function was introduced and the interface conditions were transformed to standard Cauchy singular integral equations. Solving these equations gives rise to the COD at each debonding/crack. Compared with the method of Norris and Yang, the one proposed by Y. S. Wang and D. Wang has a high convergence rate especially when the size of cracks is small. Besides these theoretical investigations, numerical simulations using the boundary element methods were also conducted, for example, Biwa and Shibata [19, 20] and Sato and Shindo [21].The perfect bonding and the debonding are the two extreme cases of the bonding conditions of interfaces. Actually, there are other conditions which extensively exist between the perfect bonding and the debonding, that is, imperfect condition. For this reason, several theoretical models were proposed during the past few years [22–24]. Among them, the spring model [22] attracted the most attention. In the spring model, the stresses are assumed to be continuous across the interface, while the displacements jump. Furthermore, the displacement jump is assumed to be linearly proportional to the stress at the interface. The stress to the displacement jump ratio is defined as the stiffness of springs. Thus, the stiffness of the spring reflects the bonding level of the imperfect interface. It is clear that when the stiffness of the spring tends towards infinity, the imperfect bonding becomes the usual perfect bonding, while it approaches the debonding when the stiffness approaches zero.Using the above-introduced spring model, the effective mechanical properties of composite materials, such as the effective modulus [25, 26], thermal expansion coefficient [27], phase velocity, and attenuation [5, 8, 28], have been extensively investigated by taking the imperfect bonding into account. The spring model has also been used to study the influence of the imperfect bonding effect on the behaviors of the underground tunnels under various dynamic loadings [15, 24, 29, 30].It is worthwhile to mention that, in most of the previous studies, the whole interface between each fiber/particulate and the matrix was usually assumed to simultaneously enter the imperfect bonding condition. However, this is not real case, assuming that the interface between each fiber and the surrounding matrix deteriorates as a process would be more reasonable. In this regard, less attention has been paid except for the partial debonding case. The limited works include those of Lopez-Realpozo et al. [31] and Guinovart-Díaz et al. [32], where the effect of the partially imperfect bonding on the effective mechanical properties of fiber-reinforced composites was investigated. To the best of our knowledge, investigation on the influence of this partially imperfect bonding of interfaces on the behavior of P- and SV-waves in composites is still scarce. It is well known that the study of the single scattering is crucial to understanding the overall behavior of elastic waves in composites. Towards this end, in this work, the single scattering of P- and SV-waves by a cylindrical fiber with a partially imperfect bonding to the surrounding matrix is studied. Compared with the single scattering with a fully perfect or imperfect interface, solving the single scattering with a partially imperfect bonding is much more complicated, since the conditions at the interface become discontinuous. In this work, the imperfect bonding is characterized by the spring model and a CP method is introduced to solve the corresponding single scattering problem.In the following sections, firstly, the single scattering problem is formulated. Then, the CP method is introduced to solve the single scattering problem. After that, the influence of various aspects of the partially imperfect bonding on the single scattering is extensively studied numerically. Finally, a short conclusion is drawn.


## 2. Single Scattering of P/SV-Waves by a Fiber with a Partially Imperfect Bonding

### 2.1. Governing Equations for P- and SV-Waves in 2D Problems

The governing equations for P- and SV-waves in 2D problems with a homogeneous medium are decoupled based on the Helmholtz decomposition. In harmonic analysis, the two equations for P- and SV-waves can be expressed as(1)∇2+kl2ϕ=0∇2+ks2φ=0,where ϕ and φ denote the displacement potentials for P- and SV-waves, respectively; and ∇2 is the Laplace operator. The symbols kl and ks are the wavenumbers of the P- and SV-waves, and the subscripts l and s represent the longitudinal and transverse wave. The wavenumbers are related to material properties as(2)kl=ωλ+2μ/ρ;ks=ωμ/ρ,where ω is the circular frequency, λ is the Lame constant, μ is the shear modulus, and ρ is the mass density. In the polar coordinate system, the displacements and stresses are expressed as(3)ur=∂ϕ∂r+1r∂φ∂θuθ=1r∂ϕ∂θ-∂φ∂rσrr=λ+2μ∂2ϕ∂r2+λ1r2∂2ϕ∂θ2+1r∂ϕ∂r+2μ1r∂2φ∂r∂θ-1r2∂φ∂θσrθ=μ-2r2∂ϕ∂θ+2r∂2ϕ∂r∂θ+1r2∂2φ∂θ2-∂2φ∂r2+1r∂φ∂r,where ur and uθ represent the radial and circumferential displacements, and σrr and σrθ are the radial and tangential stresses. The formula in ((1)–(3)) can be found in numerous classic books, of which we mention the monograph by Mow and Pao [33]. They are copied here to maintain the integrity of this work.


### 2.2. Problem Formulation of the 2D Singe Scattering of P- and SV-Waves

Figure1 shows a schematic of the single scattering of a plane P- or SV-wave by a cylindrical fiber partially imperfectly bonded to the surrounding matrix. The radius of the fiber is a. The range of the imperfect interface is represented by angle β, as indicated by the red solid curve shown in Figure 1. The symmetric axis of the imperfect interface is at angle α with the horizontal axis. It is clear that when β=0, the whole fiber is fully perfectly bonded to the matrix, and the whole interface becomes a fully imperfect bonding as β=2π. In this work, the imperfect bonding is modelled using the spring model introduced before, where it is characterized using two springs distributed along the radial and circumferential directions, with stiffness Kn and stiffness Kt, respectively. The incident wave is a plane P- or SV-wave propagating in the horizontal direction. In the polar coordinate system, general solutions to (1) are the products of Bessel functions and trigonometric functions. Hence, for harmonic analysis, the incident wave, the wave scattered by the fiber, and the wave fields inside the fiber can be expressed as series of cylindrical wave functions with the omission of e-iωt:(4)ϕincr,θ=eiklmx=∑n=-∞∞inJnklmreinθϕmr,θ=∑n=-∞∞AnHn1klmreinθφmr,θ=∑n=-∞∞BnHn1ksmreinθϕfr,θ=∑n=-∞∞CnJnklfreinθφfr,θ=∑n=-∞∞DnJnksfreinθ,where (r,θ) are the polar coordinates with origins at the center of the fiber. The wavenumbers kl,sα=ω/cl,sα (α=mandf), where cl,sα is the wave speed of longitudinal (l) wave or transverse (s) wave. The superscripts m and f indicate quantities associated with the matrix and the fiber, respectively. Jn(·) and Hn(1)(·) are the Bessel and Hankel functions of the first kind; i=-1 is the imaginary unit; and An, Bn, Cn, and Dn are the unknown coefficients, which are to be determined by the conditions at the interface. In the current case, the conditions of stresses at the interface are listed as(5)σrrmr=a=σrrfr=a,σrθmr=a=σrθfr=a,θ∈0,2π.And the conditions for the radial and circumferential displacements at the interface can be expressed as(6)Knurmr=a-urfr=a=σrrfr=a,θ∈β0,otherwiseKtuθmr=a-uθfr=a=σrθfr=a,θ∈β0,otherwise.Figure 1
Scattering of a plane P- or SV-wave by a cylindrical fiber with a partially imperfect bonding to the matrix. The range of the imperfect bonding is denoted by the red solid arc.The total displacement potential in the matrix is the summation of the incident wave and the scattered wave, while the quantities in the fiber are only of the refracted one. By substituting the total displacement potentials in the matrix and the fiber expressed in (4) into (3), the displacements and stresses in the matrix and the fiber are obtained. Then, the interface conditions in (5) and (6) can be expressed in terms of the expansion coefficients An, Bn, Cn, and Dn. Firstly, based on (5), two equations are established for the expansion coefficients as(7)M31M32M33M34M41M42M43M44AnBnCnDn=b3b4,n=-∞,∞.The detailed expressions for Mij and bi(i=3,4andj=1,2,3,4) are listed in Appendix A. After a simple mathematical operation on (7), the coefficients An and Bn can be expressed in terms of Cn and Dn as(8)AnBn=β1nβ2n+α11nα12nα21nα22nCnDn,n=-∞,∞with β1(n)β2(n)=M31M32M41M42-1b3b4 and α11(n)α12(n)α21(n)α22(n)=-M31M32M41M42-1M33M43M43M44, where the superscript “−1” means the inverse operation on a square matrix. Then, the remaining conditions for the radial and circumferential displacements in (6) can be expressed in terms of coefficients Cn and Dn:(9)∑n=-∞∞cj1n,θCn+∑n=-∞∞cj2n,θDn=∑n=-∞∞hjn,θ,j=1,2,3,4.The details for cij and hi(i=1,2,3,4;j=1,2) are referred to in Appendix B. Note that, in (9), θ is within the imperfectly bonded interface as j=1 and 3 and is within the perfectly bonded interface for j=2 and 4. Since θ does not vary from 0 to 2π in any equation of (9), the orthogonality for the sine and cosine functions is not available. Therefore, analytical solutions for Cn and Dn are difficult to be obtained.


## 2.1. Governing Equations for P- and SV-Waves in 2D Problems

The governing equations for P- and SV-waves in 2D problems with a homogeneous medium are decoupled based on the Helmholtz decomposition. In harmonic analysis, the two equations for P- and SV-waves can be expressed as(1)∇2+kl2ϕ=0∇2+ks2φ=0,where ϕ and φ denote the displacement potentials for P- and SV-waves, respectively; and ∇2 is the Laplace operator. The symbols kl and ks are the wavenumbers of the P- and SV-waves, and the subscripts l and s represent the longitudinal and transverse wave. The wavenumbers are related to material properties as(2)kl=ωλ+2μ/ρ;ks=ωμ/ρ,where ω is the circular frequency, λ is the Lame constant, μ is the shear modulus, and ρ is the mass density. In the polar coordinate system, the displacements and stresses are expressed as(3)ur=∂ϕ∂r+1r∂φ∂θuθ=1r∂ϕ∂θ-∂φ∂rσrr=λ+2μ∂2ϕ∂r2+λ1r2∂2ϕ∂θ2+1r∂ϕ∂r+2μ1r∂2φ∂r∂θ-1r2∂φ∂θσrθ=μ-2r2∂ϕ∂θ+2r∂2ϕ∂r∂θ+1r2∂2φ∂θ2-∂2φ∂r2+1r∂φ∂r,where ur and uθ represent the radial and circumferential displacements, and σrr and σrθ are the radial and tangential stresses. The formula in ((1)–(3)) can be found in numerous classic books, of which we mention the monograph by Mow and Pao [33]. They are copied here to maintain the integrity of this work.


## 2.2. Problem Formulation of the 2D Singe Scattering of P- and SV-Waves

Figure1 shows a schematic of the single scattering of a plane P- or SV-wave by a cylindrical fiber partially imperfectly bonded to the surrounding matrix. The radius of the fiber is a. The range of the imperfect interface is represented by angle β, as indicated by the red solid curve shown in Figure 1. The symmetric axis of the imperfect interface is at angle α with the horizontal axis. It is clear that when β=0, the whole fiber is fully perfectly bonded to the matrix, and the whole interface becomes a fully imperfect bonding as β=2π. In this work, the imperfect bonding is modelled using the spring model introduced before, where it is characterized using two springs distributed along the radial and circumferential directions, with stiffness Kn and stiffness Kt, respectively. The incident wave is a plane P- or SV-wave propagating in the horizontal direction. In the polar coordinate system, general solutions to (1) are the products of Bessel functions and trigonometric functions. Hence, for harmonic analysis, the incident wave, the wave scattered by the fiber, and the wave fields inside the fiber can be expressed as series of cylindrical wave functions with the omission of e-iωt:(4)ϕincr,θ=eiklmx=∑n=-∞∞inJnklmreinθϕmr,θ=∑n=-∞∞AnHn1klmreinθφmr,θ=∑n=-∞∞BnHn1ksmreinθϕfr,θ=∑n=-∞∞CnJnklfreinθφfr,θ=∑n=-∞∞DnJnksfreinθ,where (r,θ) are the polar coordinates with origins at the center of the fiber. The wavenumbers kl,sα=ω/cl,sα (α=mandf), where cl,sα is the wave speed of longitudinal (l) wave or transverse (s) wave. The superscripts m and f indicate quantities associated with the matrix and the fiber, respectively. Jn(·) and Hn(1)(·) are the Bessel and Hankel functions of the first kind; i=-1 is the imaginary unit; and An, Bn, Cn, and Dn are the unknown coefficients, which are to be determined by the conditions at the interface. In the current case, the conditions of stresses at the interface are listed as(5)σrrmr=a=σrrfr=a,σrθmr=a=σrθfr=a,θ∈0,2π.And the conditions for the radial and circumferential displacements at the interface can be expressed as(6)Knurmr=a-urfr=a=σrrfr=a,θ∈β0,otherwiseKtuθmr=a-uθfr=a=σrθfr=a,θ∈β0,otherwise.Figure 1
Scattering of a plane P- or SV-wave by a cylindrical fiber with a partially imperfect bonding to the matrix. The range of the imperfect bonding is denoted by the red solid arc.The total displacement potential in the matrix is the summation of the incident wave and the scattered wave, while the quantities in the fiber are only of the refracted one. By substituting the total displacement potentials in the matrix and the fiber expressed in (4) into (3), the displacements and stresses in the matrix and the fiber are obtained. Then, the interface conditions in (5) and (6) can be expressed in terms of the expansion coefficients An, Bn, Cn, and Dn. Firstly, based on (5), two equations are established for the expansion coefficients as(7)M31M32M33M34M41M42M43M44AnBnCnDn=b3b4,n=-∞,∞.The detailed expressions for Mij and bi(i=3,4andj=1,2,3,4) are listed in Appendix A. After a simple mathematical operation on (7), the coefficients An and Bn can be expressed in terms of Cn and Dn as(8)AnBn=β1nβ2n+α11nα12nα21nα22nCnDn,n=-∞,∞with β1(n)β2(n)=M31M32M41M42-1b3b4 and α11(n)α12(n)α21(n)α22(n)=-M31M32M41M42-1M33M43M43M44, where the superscript “−1” means the inverse operation on a square matrix. Then, the remaining conditions for the radial and circumferential displacements in (6) can be expressed in terms of coefficients Cn and Dn:(9)∑n=-∞∞cj1n,θCn+∑n=-∞∞cj2n,θDn=∑n=-∞∞hjn,θ,j=1,2,3,4.The details for cij and hi(i=1,2,3,4;j=1,2) are referred to in Appendix B. Note that, in (9), θ is within the imperfectly bonded interface as j=1 and 3 and is within the perfectly bonded interface for j=2 and 4. Since θ does not vary from 0 to 2π in any equation of (9), the orthogonality for the sine and cosine functions is not available. Therefore, analytical solutions for Cn and Dn are difficult to be obtained.


## 3. Scattering Cross Section and the Far-Field Scattering Magnitude

In order to characterize the overall behavior of coherent waves in composites with randomly distributed fibers/particulates, over the past several decades, extensive theoretical models have been proposed including theT-matrix method [34–36], the self-consistent method [37–40], and the phenomenological model (also known as the independent scattering model) [1, 2, 41, 42]. It is noted that all these theoretical models necessitate the solution of the corresponding single scattering problem to obtain the parameters, such as the total SCS and the far-field scattering magnitude. Therefore, the current study is indispensable, especially taking the imperfection of interfaces into account. In the current case (2D P- and SV-waves), the total SCS can be calculated by(10)σsα,β=4klm∑n=-∞∞An2+Bn2.The far-field scattering magnitude (r→∞) is defined as(11)u⃑s=urse⃑r+uθse⃑θ=iklm2πklmreiklmr-π/4flθe⃑r+iksm2πksmreiksmr-π/4fsθe⃑θ.Here a variable with an arrow represents a vector. The superscript s denotes quantities of the scattered wave. By substituting the expression of the scattered wave into the above equation, the far-field scattering magnitudes for the scattered P- and SV-waves can be expressed as [13](12)flθ=∑n=-∞∞-inAneinθfsθ=∑n=-∞∞-inBneinθ.σs(θ), flθ, and fs(θ) are evaluators of the scattering capability of the fiber. Accurate calculation of these parameters plays a critical role in these theoretical models.


## 4. The Collocation Point Method

In this study, the CP method is adopted to solve (9). First, the infinite series in (4) are truncated with the maximum order nmax. Then, the number of the unknowns Cn and Dn in (9) becomes (2nmax+1). In order to solve these unknowns, both the first (j=1,2) and second (j=3,4) group equations in (9) are forced to be satisfied at (2nmax+1) CPs. In this work, these CPs are distributed symmetrically about the partially imperfectly and perfectly bonded interfaces, as shown in Figure 2. Additionally, these CPs are required to be distributed as uniformly as possible on the whole interface. Then, a linear system of size 2×(2nmax+1) is established for Cn and Dn. Solving this linear system yields Cn and Dn (n=-nmax,nmax); the coefficients An and Bn can then be calculated according to (8). Then the displacements and stresses at each point inside the fiber and the matrix can be calculated.Figure 2
Distribution of CPs denoted by the “×” symbol. 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.To check the convergence of this CP method, the total SCS by a partially debonded fiber (α=0 and β=2π/3) under different nmax is calculated and plotted in Figure 3. In this case, the materials for the fiber and the matrix are glass and epoxy, whose detailed properties are listed in Table 1. The results show that the calculated SCS converges to a constant as nmax increases. Based on this convergence analysis, in this study, nmax is set as 10.Table 1
Materials used in the current work [4, 13].
Material λ (GPa) μ (GPa) ρ (kg/m3) SiC 92.1 177.0 3200 Ti-alloy 103.0 44.8 5400 Glass 0.7367 1.43 2550 Epoxy 0.8895 1.28 1250Figure 3
Convergence study of the CP method using the debonding case considered in [13]. In the CP method, the stiffness Kn/Kn∗=Kt/Kt∗=10-6 represents the debonding. In this work, Kn∗ and Kt∗ are defined as Kn∗=0.5(λm+λf)/0.1a and Kt∗=0.5(μm+μf)/0.1a.Next, three examples are calculated to validate the application of the CP method. Figure4 shows the calculated SCS by a partially imperfectly bonded fiber changing with the stiffness at klma=3.0. The materials used are the same as those in Figure 3. The dash lines represent values of the debonding and the perfect bonding cases obtained in [13]. The results from this study approach the exact solutions when the stiffness tends to zero and infinity. Moreover, the calculated results under different β tend to the same value as the stiffness increases, which follows the fact that all the imperfect bonding of different widths becomes the perfect one as the stiffness is large enough. In the following two calculations, the matrix is Ti-alloy and the fiber is made up of SiC. The detailed properties of the two materials are tabulated in Table 1. Figure 5 shows the calculated far-field scattering magnitude patterns for the scattered P- and SV-waves under a plane P-wave incidence with a normalized frequency klma=1.0. The imperfect bonding is located at α=0 with a width β=π/2. As shown in Figure 5, patterns at several different bonding levels are plotted and compared with that of the fully perfect bonding case. Results show that the pattern with a partially imperfect bonding approaches that of the fully perfect bonding case as the stiffness Kn and stiffness Kt increase, which agrees with the theoretical expectation.Figure 4
The calculated SCS by a partially imperfectly bonded fiber as a function of the stiffness in the range of10-4 to 104.Figure 5
The far-field scattering magnitude patterns of the scattered (a) P- and (b) SV-waves by a fiber with a partially imperfect bonding to the matrix, under a plane P-wave incidence (dashed line:Kn/Kn∗=Kt/Kt∗=10-1; dotted line: Kn/Kn∗=Kt/Kt∗=100; dash dotted line: Kn/Kn∗=Kt/Kt∗=101; solid line: perfect bonding).
(a) (b)Figure6 shows the calculated far-field scattering magnitude patterns of the scattered P- and SV-waves by a fiber with a fully (β=2π) imperfect bonding under a P-wave incidence at a normalized frequency klma=3.0. Results show that patterns of both the scattered P- and SV-waves approach those of the cavity one as the stiffness Kn and stiffness Kt decrease, which is also in agreement with the theoretical expectation.Figure 6
Comparison of the calculated far-field scattering patterns of the scattered (a) P-wave and (b) SV-wave by a fiber with a fully imperfect bonding to the matrix under a P-wave incidence versus those of the debonding case (cavity).
(a) (b)Figures4, 5, and 6 indicate that the CP method performs well to calculate the total SCS and the far-field scattering magnitude for the single scattering problem of P- or SV-waves by a fiber with a partially imperfect bonding to the matrix. But, it should be reminded that the current CP method is not a proper technique to calculate the corresponding crack opening displacements and the stresses at the interface, especially around the tips of the imperfect interface. Moreover, the results show that the spring model with the setting of Kn/Kn∗=Kt/Kt∗=101 and Kn/Kn∗=Kt/Kt∗=10-3 properly modelled the perfect bonding and debonding situations, respectively.


## 5. Effects of the Partially Imperfect Bonding on the Single Scattering

### 5.1. Effects ofα and β

In this section, the influence ofα and β on the single scattering of P- and SV-waves is investigated. The stiffness of the imperfectly bonded interface is set as Kn/Kn∗=Kt/Kt∗=10-2, which represents a serious imperfection close to the debonding. Figure 7 shows the relative divergence of the calculated σs under different α and β from that of the fully perfect bonding, under a P-wave incidence. Based on the results shown in this figure, several conclusions are obtained. First, it is clear that σs is almost symmetric about the axis α=π/2. Second, within the range of investigated β, the partially imperfect bonding located at α=0 and π reduces σs the most. In contrast, the partially imperfect bonding located at α=π/2 has almost no effect on σs compared with the fully perfect bonding, especially as the frequency is relatively high; all figures except Figure 7(a) show this. Theoretically, this fact benefits the identification of the specific location of the imperfection. Finally, as the frequency increases to values larger than klma=2.0, the partially imperfect bonding always reduces σs compared with the fully perfect bonding, wherever the partially imperfect bonding is located.Figure 7
Relative difference (%) of the calculatedσs by a fiber with a partially imperfect bonding and with the fully perfect bonding, under the P-wave incidence at several different normalized frequencies.
(a)
k l m a = 1.0 (b)
k l m a = 1.5 (c)
k l m a = 2.0 (d)
k l m a = 2.5 (e)
k l m a = 3.0The corresponding results under the SV-wave incidence are plotted in Figure8. The frequencies of the incident wave are the same for each pair of figures in Figures 7 and 8. Similar to the results of the P-wave incidence, σs under the SV-wave incidence is also symmetric about the axis α=π/2. Moreover, for each specific β, the partially imperfect bonding located at α=π/2 reduces σs the most. Similarly, this fact is also beneficial to the identification of the location of the imperfection. In addition, contrary to the P-wave incidence, the partially imperfect bonding always reduces σs at relatively lower frequencies rather than at higher frequencies.Figure 8
Relative difference (%) of the calculatedσs by a fiber with a partially imperfect bonding and with the fully perfect bonding, under the SV-wave incidence at several different normalized frequencies.
(a)
k s m a = 2.07 (b)
k s m a = 3.11 (c)
k s m a = 4.15 (d)
k l m a = 5.18 (e)
k l m a = 6.22


### 5.2. The Stiffness Effect

In this section, influence of stiffnessKn and stiffness Kt of the imperfect interface on the single scattering is studied. Figure 9 shows the calculated σs changing with Kn and Kt at three different frequencies under the P- and SV-wave incidence, respectively. The partially imperfect bonding is located at α=0 with a width β=π/4. It is clear that σs is sensitive to Kn and insensitive to Kt under the P-wave incidence. On the contrary, under the SV-wave incidence, σs becomes sensitive to Kt rather than Kn.Figure 9
The influence of stiffnessKn and stiffness Kt on σs under a P- and SV-wave incidence, respectively. (a) and (b) show results under the P-wave incidence and (c) and (d) show results under the SV-wave incidence. In (a) and (c) Kt/Kt∗=1, and in (b) and (d) Kn/Kn∗=1.
(a) (b) (c) (d)Definitely, the sensitivity ofσs to Kn and Kt is affected by the location of the partially imperfect bonding, based on the results shown in Figures 7 and 8. In this work, this sensitivity is defined using the slope κ of the solid line shown in Figure 10, which is the linear fitting to the curve within 10-1,100. A larger absolute value of κ means σs has a higher sensitivity to the stiffness of the imperfectly bonded interface. Figure 11 shows values of κ under the P-wave incidence as the imperfect bonding is located at different positions. In this figure, β=π/4. The results indicate that σs is more sensitive to Kn than to Kt. Moreover, σs changes significantly in response to Kn when α is small, as Figure 11(a) shows. In addition, the sensitivity of σs to Kt reaches a maximum as α equals approximately π/4. Figure 11(c) shows the sensitivity of σs to both Kn and Kt.Figure 10
A schematic to define the sensitivity ofσs to Kn or Kt.Figure 11
The sensitivity ofσs to (a) Kn and (b) Kt and (c) both Kn and Kt, under the P-wave incidence. Kt/Kt∗=1 in (a), Kn/Kn∗=1 in (b), and Kn/Kn∗=Kt/Kt∗ in (c).
(a) (b) (c)The corresponding results under the SV-wave incidence are plotted in Figure12. Results show that σs is sensitive to Kn as α is larger, as Figure 12(a) shows. In contrast, σs becomes sensitive to Kt for smaller α, as shown in Figure 12(b). Figure 12(c) shows the sensitivity of σs to both Kn and Kt under the SV-wave incidence. Comparison of the results in Figures 11(c) and 12(c) indicates that when α is small, the P-wave incidence at higher frequencies is more applicable to evaluate the changing of the stiffness of the imperfectly bonded interface. In contrast, the SV-wave incidence at higher frequencies is more applicable to evaluate the changing of the stiffness of the imperfectly bonded interface as α is large.Figure 12
The sensitivity ofσs to (a) Kn and (b) Kt and (c) both Kn and Kt, under the SV-wave incidence. Kt/Kt∗=1 in (a), Kn/Kn∗=1 in (b), and Kn/Kn∗=Kt/Kt∗ in (c).
(a) (b) (c)


## 5.1. Effects ofα and β

In this section, the influence ofα and β on the single scattering of P- and SV-waves is investigated. The stiffness of the imperfectly bonded interface is set as Kn/Kn∗=Kt/Kt∗=10-2, which represents a serious imperfection close to the debonding. Figure 7 shows the relative divergence of the calculated σs under different α and β from that of the fully perfect bonding, under a P-wave incidence. Based on the results shown in this figure, several conclusions are obtained. First, it is clear that σs is almost symmetric about the axis α=π/2. Second, within the range of investigated β, the partially imperfect bonding located at α=0 and π reduces σs the most. In contrast, the partially imperfect bonding located at α=π/2 has almost no effect on σs compared with the fully perfect bonding, especially as the frequency is relatively high; all figures except Figure 7(a) show this. Theoretically, this fact benefits the identification of the specific location of the imperfection. Finally, as the frequency increases to values larger than klma=2.0, the partially imperfect bonding always reduces σs compared with the fully perfect bonding, wherever the partially imperfect bonding is located.Figure 7
Relative difference (%) of the calculatedσs by a fiber with a partially imperfect bonding and with the fully perfect bonding, under the P-wave incidence at several different normalized frequencies.
(a)
k l m a = 1.0 (b)
k l m a = 1.5 (c)
k l m a = 2.0 (d)
k l m a = 2.5 (e)
k l m a = 3.0The corresponding results under the SV-wave incidence are plotted in Figure8. The frequencies of the incident wave are the same for each pair of figures in Figures 7 and 8. Similar to the results of the P-wave incidence, σs under the SV-wave incidence is also symmetric about the axis α=π/2. Moreover, for each specific β, the partially imperfect bonding located at α=π/2 reduces σs the most. Similarly, this fact is also beneficial to the identification of the location of the imperfection. In addition, contrary to the P-wave incidence, the partially imperfect bonding always reduces σs at relatively lower frequencies rather than at higher frequencies.Figure 8
Relative difference (%) of the calculatedσs by a fiber with a partially imperfect bonding and with the fully perfect bonding, under the SV-wave incidence at several different normalized frequencies.
(a)
k s m a = 2.07 (b)
k s m a = 3.11 (c)
k s m a = 4.15 (d)
k l m a = 5.18 (e)
k l m a = 6.22


## 5.2. The Stiffness Effect

In this section, influence of stiffnessKn and stiffness Kt of the imperfect interface on the single scattering is studied. Figure 9 shows the calculated σs changing with Kn and Kt at three different frequencies under the P- and SV-wave incidence, respectively. The partially imperfect bonding is located at α=0 with a width β=π/4. It is clear that σs is sensitive to Kn and insensitive to Kt under the P-wave incidence. On the contrary, under the SV-wave incidence, σs becomes sensitive to Kt rather than Kn.Figure 9
The influence of stiffnessKn and stiffness Kt on σs under a P- and SV-wave incidence, respectively. (a) and (b) show results under the P-wave incidence and (c) and (d) show results under the SV-wave incidence. In (a) and (c) Kt/Kt∗=1, and in (b) and (d) Kn/Kn∗=1.
(a) (b) (c) (d)Definitely, the sensitivity ofσs to Kn and Kt is affected by the location of the partially imperfect bonding, based on the results shown in Figures 7 and 8. In this work, this sensitivity is defined using the slope κ of the solid line shown in Figure 10, which is the linear fitting to the curve within 10-1,100. A larger absolute value of κ means σs has a higher sensitivity to the stiffness of the imperfectly bonded interface. Figure 11 shows values of κ under the P-wave incidence as the imperfect bonding is located at different positions. In this figure, β=π/4. The results indicate that σs is more sensitive to Kn than to Kt. Moreover, σs changes significantly in response to Kn when α is small, as Figure 11(a) shows. In addition, the sensitivity of σs to Kt reaches a maximum as α equals approximately π/4. Figure 11(c) shows the sensitivity of σs to both Kn and Kt.Figure 10
A schematic to define the sensitivity ofσs to Kn or Kt.Figure 11
The sensitivity ofσs to (a) Kn and (b) Kt and (c) both Kn and Kt, under the P-wave incidence. Kt/Kt∗=1 in (a), Kn/Kn∗=1 in (b), and Kn/Kn∗=Kt/Kt∗ in (c).
(a) (b) (c)The corresponding results under the SV-wave incidence are plotted in Figure12. Results show that σs is sensitive to Kn as α is larger, as Figure 12(a) shows. In contrast, σs becomes sensitive to Kt for smaller α, as shown in Figure 12(b). Figure 12(c) shows the sensitivity of σs to both Kn and Kt under the SV-wave incidence. Comparison of the results in Figures 11(c) and 12(c) indicates that when α is small, the P-wave incidence at higher frequencies is more applicable to evaluate the changing of the stiffness of the imperfectly bonded interface. In contrast, the SV-wave incidence at higher frequencies is more applicable to evaluate the changing of the stiffness of the imperfectly bonded interface as α is large.Figure 12
The sensitivity ofσs to (a) Kn and (b) Kt and (c) both Kn and Kt, under the SV-wave incidence. Kt/Kt∗=1 in (a), Kn/Kn∗=1 in (b), and Kn/Kn∗=Kt/Kt∗ in (c).
(a) (b) (c)


## 6. Conclusions

The single scattering of P- and SV-waves by a cylindrical fiber with a partially imperfect bonding to the matrix is investigated using the CP method. This study benefits the characterization of coherent waves in composites with randomly distributed inclusions. In this study, the imperfect interface is modelled using the spring model. Three examples are calculated and the results show that this CP method performs well. Also, the parametric analysis of location, width, and bonding level of the imperfect bonding on the total SCS is extensively studied. From this study, the following conclusions are obtained:(a)
The influence of the partially imperfect bonding on the single scattering of P- and SV-waves is almost symmetric about the axisα=π/2.(b)
Under the P-wave incidence, theσs value reduced by the partially imperfect bonding at the location of α=0 and π is greater than that at any other location, compared with the fully perfect bonding. The partially imperfect bonding located at α=π/2 has almost no effect on σs. Under the SV-wave incidence, the partially imperfect bonding located at α=π/2 significantly reduces σs.(c)
Whenα is small, σs under the P-wave incidence has a high sensitivity to the level of the imperfect bonding. However, σs under the SV-wave incidence becomes more sensitive to the level of the imperfect bonding as α is large.In summary, besides its significance in the understanding of the single scattering, the present study also benefits for the application of the theoretical models to evaluate phase velocity and attenuation coefficient of coherent elastic waves in composites with randomly distributed inclusions and the nondestructive evaluation of composites using the ultrasonic techniques.



---

*Source: 7516402-2018-04-23.xml*