# Ultrasonic Nonlinearity Evaluation of the Cracked Interface.

1. Introduction

Adhesion bonding technology for fiber reinforced composite has been widely applied in aerospace industry. Generally, the cohesive interface is very thin but can provide the powerful strength and resiliency. However, the capability of the cohesive interface could be greatly affected by the imperfect condition, such as debond, fatigue damage, and microcracks. Thus, there is a great need for development of nondestructive testing methods to detect the imperfect condition.

Nonlinear ultrasonic methods have the powerful ability to characterize the material nonlinearity change caused by plasticity [1], imperfect interface [2-4], microcracks [5, 6], and fatigue damage [7,8]. When a time-harmonic longitudinal wave propagates through a cracked solid or interface, it will cause the tension and compression elastic asymmetry; then the waveform will be distorted and higher-order harmonic waves are generated [9-11]. And considerable experimental evidence [12,13] has shown that ultrasonic waves do interact with microcracks in a nonlinear fashion, but those researches may fail when the crack size is only tens of microns, which is one of the initial factors leading to interface degradation.

For the researches on acoustic nonlinearity induced by microcracks [14-16], majority are concentrated either on the scattering of elastic waves by a single crack or an array of cracks [17-20] or on the propagation of elastic waves in a cracked medium [21-23]. The first study on such contact-induced acoustic nonlinearity is probably the paper by Richardson [24] who considered the contact interface between two semi-infinite half-spaces. He has analyzed one-dimensional nonlinear wave propagation in a system composed of an unbounded planer interface separating two semi-infinite linear elastic media. The nonlinearity is caused by the opening and closing of the interface. However, in Richardson's analysis, the interface stiffness varying continuously is not accounted for. Improving Richardson's theory, Biwa et al. [25] have analyzed a nonlinear interface stiffness model, where the stiffness property of the contact interface is described as a function of the nominal contact pressure. Furthermore, by assuming the interface of the adhesive as a nonlinear spring, Achenbach and Parikh [26] have investigated theoretically to obtain information on the adhesive bond strength from ultrasonic test results, and it is shown that the nonlinear adhesive bond behavior could cause the generation of higher harmonics.

In this paper, based on the researches of Richardson and Biwa, one theoretical solution for one-dimensional nonlinear wave propagation in a system composed of a multicrack-included interface separating two semi-infinite elastic media is derived, from which the acoustic nonlinearity caused by the tension and compression elastic asymmetry is clearly shown. To validate the analytical solution, numerical results from the FEM simulations are presented. Comparison between the analytical predictions and the FEM simulation results shows good agreement. Finally, we also study two factors affecting acoustic nonlinearity. One is the tension and compression elastic asymmetry of the cracked interface and another is the frequency of incident wave.

2. Solution for Harmonic Wave Incidence

We consider a system composed of a cracked interface separating two semi-infinite elastic media. In the presence of microcracks, the interface will respond to the tensile and compressive loadings differently, which is the tension and compression elastic asymmetry. According to [11], the elastic constants of the cracked medium can be considered to be dependent with crack density, the friction of the crack faces, and frequency of incident wave:

[mathematical expression not reproducible], (1)

where E' and v' are Young's modulus and Poisson's ratio of the uncracked solid, respectively. And [[bar.E].sub.[alpha]] and [[bar.v].sub.[alpha]] are, respectively, the corresponding effective Young's modulus and Poisson's ratio of the cracked solid under tension ([alpha] = t) and under compression ([alpha] = c). Meanwhile, E' = E and v' = v for plane stress, and E' = E/( 1 - [v.sup.2]) and v = v/(1 - v) for plane strain. The crack density d is defined as d = [Na.sup.2]/A, where a is the average half-length of the cracks; A is the area of the solid containing N randomly distributed and randomly oriented two-dimensional microcracks.

We have

[mathematical expression not reproducible]. (2)

In the above, [gamma] is the Euler-Mascheroni constant, [eta] = [k.sub.T]a = [omega]a/[c.sub.T] is a dimensionless wavenumber, [kappa] = [c.sub.L]/[c.sub.T] = [square root of (2(1 - v)/(1 - 2v))] is a dimensionless parameter, [mu] is the coefficient of friction, and

[mathematical expression not reproducible]. (3)

Therefore, in this paper, the conclusion of [11] will be adopted, and the interface can be regarded as an equivalent medium consisting of the random distribution cracks. Actually this assumption has been proven in [25] for investigating the void inclusion medium subjected to the ultrasonic wave loading conditions.

2.1. The General Solution. Firstly, we assume that the elasticity of the cracked interface is a function of time, and the stress in the interface is uniform. In Figure 1 we illustrate a one-dimensional system schematically to consider elastic longitudinal wave propagation along x-axis, where f(x - ct) and g(x - ct) represent the incident and transmitted wave functions, respectively. The blue area represents the cracked interface.

The two semi-infinite elastic solids are located in the regions x < [0.sup.-] and x > [0.sup.+], respectively. The equation of responses of the interface and the stress-strain relation are defined as follows:

[mathematical expression not reproducible], (4)

where [rho] is mass density, [E.sub.s](t) is the longitudinal stiffness of the cracked interface, u(x, t) is the displacement in the x-direction of an element at time t from its position x, and [sigma](x, t) is the stress. The boundary conditions in locations [0.sup.-] and [0.sup.+] are as follows:

[sigma](t) = [sigma]([0.sup.-], t) = [sigma]([0.sup.+], t), (5)

u([0.sup.+], t) - u([0.sup.-], t) = [sigma](t) x h/[E.sub.s], (6)

where h is the thickness of the cracked interface.

According to [24,25], the governing equation of (4) can be solved based on travelling wave method. Therefore, this paper will employ travelling wave method to obtain the general displacement solution when longitudinal wave propagates through the interface. Then considering the cracked interface as the tension and compression elastic asymmetry, this paper will achieve the corresponding special displacement solution.

According to the travelling wave method, as a solution to (4), the following forms are considered [24,25]:

u(x, t) = f(x - ct) - 1/2V(t + x/c), x < [0.sup.-], (7)

u(x, t) = f(x - ct) + 1/2v(t - x/c), x [greater than or equal to] [0.sup.+], (8)

V(t) = u([0.sup.+], t) - u([0.sup.-], t), (9)

where c = [square root of (E/[rho])] is the wave velocity, E is the longitudinal stiffness of the two semi-infinite elastic solids, f(x - ct) represents the incident waves, the term -(l/2)V(t + x/c) represents reflected waves, and (8) represents the transmitted waves. Equation (9) represents the equivalent displacement of the cracked interface.

According to (6) and (9), we can obtain

[mathematical expression not reproducible]. (10)

The incident wave f(x - ct) is now assumed in the following form:

f(x - ct) = [F.sub.0] cos [k(x - ct)], (11)

where k is the wave number and [F.sub.0] is the maximum amplitude.

Substitution of (11) in (10) yields

dV(t)/dt + [xi]V(t) = 2[F.sub.0][omega] sin ([omega]t), (12)

where [omega] is the angular frequency and [xi] = [xi](t) = 2c[E.sub.s](t)/hE.

Therefore, the general solution of (12) is

[mathematical expression not reproducible]. (13)

When we know the function form of [xi](t), we can obtain the special solution from (13).

2.2. The Special Solution. Assuming the different tensile and compression effective moduli of the cracked interface from [11], the function [E.sub.s](t) can be defined by

[mathematical expression not reproducible], (14)

where [E.sup.t.sub.s] and [E.sup.c.sub.s] are the tensile and compression effective moduli, respectively. And the parameter [gamma] = ([E.sup.t.sub.s] - [E.sup.c.sub.s])/E is introduced to represent the tension and compression asymmetry in the elastic moduli of the cracked interface.

According to (13) and (14), we can obtain

[mathematical expression not reproducible], (15)

where [phi] = [sin.sup.-1] (-[omega]/[square root of ([[xi].sup.2] + [[omega].sup.2]))].

Since V([t.sub.n]) = 0 and the calculation could begin from [t.sub.0] = 0 without loss of generality, the time points [t.sub.n+1] to switch the property come from the following equation:

[mathematical expression not reproducible]. (16)

The displacement becomes

[mathematical expression not reproducible], (17)

where

[mathematical expression not reproducible]. (18)

Substitution of (17) back to (7) and (8) gives the finial displacement field of the problem.

3. Analytical and FEM Results

3.1. Analytical Results. Numerical solutions using Matlab are performed to solve (16). For a demonstrative purpose, dashed line in Figure 2 shows one example of transmitted waveforms for an incident wave with the center frequency of 1 MHz, obtained by numerical solutions of (16) and (17), with the parameters [rho] = 2700 kg/[m.sup.3], E = 7 x [10.sup.10] Pa, [E.sup.t.sub.s] = 7 x [10.sup.9] Pa, [E.sup.c.sub.s] = 7 x [10.sup.10] Pa, [F.sub.0] = 1 x [10.sup.-5] m, and h = 1 x [10.sup.-4] m. The time-domain waveform with the dashed line is shown in Figure 2(a), and the locals in Figures 2(b1) and 2(b2) are to compare the difference of the displacement curve subjected to the first-half and second-half waveform in one-cycle loading conditions. It is seen that the first-half waveform in one cycle is distinctly different from the second-half waveform caused by the tension and compression elastic asymmetry. And Figure 2(c) shows the frequency spectrum using FFT method; it is clearly shown that the transmitted wave contains the term with angular frequency 2[omega], which is second harmonic generated by the cracked interface.

3.2. Comparison with FEM Results. The analytical solutions will be validated in this section by conducting numerical simulations of wave propagation through a cracked interface. The numerical simulations are performed using the finite element method (FEM) for the case of two-dimensional plane strain deformation. The commercial FEM software ABAQUS is used for this purpose.

To this end, a two-dimensional FEM model is constructed using the four-node plane strain (CPE4R) elements. To simulate plane waves, a rectangular strip of 0.09 m x 0.003 m is used in our FEM model. The model consists of ~300,000 four-node plane strain (CPE4R) elements. The user defined constitutive law is performed to present the different tension and compression elasticity of the cracked interface. And the interface is located in the middle of the sample with the thickness of 1 x [10.sup.-4] m. Other regions are the isotropic elastic material of aluminum. A plane longitudinal wave is generated from one end of the sample. The amplitude longitudinal wave with the frequency of 1 MHz is 1 x [10.sup.-5] m. The receiver is located at right end of the interface.

Shown in Figure 2(a) with the solid line is the FEM result of time-domain waveform of transmitted wave. The excellent comparison shows that the analytical solution results are identical to that of the FEM simulation. Furthermore, the two obvious amplitudes, being 9.37 x [10.sup.-6] m and 7.76 x [10.sup.-7] m, respectively, exist with the frequencies of ultrasonic wave at 1 MHz and 2 MHz in Figure 2(c).

4. Acoustic Nonlinearity Parameter

The time-domain waveform curve can be obtained by the analytical solution derived in Section 2; then the frequency spectrum using FFT method can also be solved to obtain [A.sub.1] and [A.sub.2], which are the amplitudes of the fundamental and second harmonic, respectively. Through the analysis of data from the analytical solution, we can obtain several observations. First, [A.sub.2] is linearly related to [A.sub.1]. This is significantly different from the second harmonic generated by quadratic nonlinearity in the solid, where [A.sub.2] is proportional to the square of [A.sub.1]. The fundamental cause of this relationship is the tension and compression elastic asymmetry of the cracked interface. Second, [A.sub.2] is scaled by parameter [gamma]. Therefore, according to [9], for acoustic nonlinearity parameter in cracked solids, the similar acoustic nonlinearity parameter of the cracked interface is introduced by

[beta] = 3[pi][A.sub.2]/[A.sub.1]kh. (19)

For investigating the relation between the linearity parameter and the frequencies of the ultrasonic acoustic, a series of analytic solutions have been performed under different frequencies including 0.5 MHz, 1 MHz, 1.5 MHz, 2.0 MHz, and 2.5 MHz. Figure 3 shows the acoustic nonlinearity parameter versus frequency. It is seen that although the acoustic nonlinearity parameter decreases linearly with increasing frequency, the decrease is not significant. In other words, the acoustic nonlinearity parameter has a rather weak dependence on frequency.

Figure 4 shows the acoustic nonlinearity parameter versus parameter [gamma]. It is seen that the acoustic nonlinearity parameter increases linearly with increasing [gamma], which can represent the degree of damage for the cracked interface. In other words, the degree of damage for the cracked interface can be characterized by the acoustic nonlinearity parameter defined as (19).

5. Conclusions

As a time-harmonic longitudinal wave propagates in a system composed of a cracked interface separating two semi-infinite elastic media, a second harmonic wave may be generated by the cracked interface. Assuming the longitudinal stiffness of the interface as a function of time, this paper solves the problem and obtains the analytical general solution (any interface) and the special solution (the cracked interface). It is clearly shown that the tension and compression elastic asymmetry can result in acoustic nonlinearity. Furthermore, to verify the developed theoretical solution, we carry out detailed numerical simulation of wave propagation through a cracked interface by using the FEM. The results indicate that the micromechanics model predictions agree well with the FEM simulations. Finally, because the amplitude of the second harmonic is linearly related to the amplitude of the fundamental harmonic, we introduce the acoustic nonlinearity parameter induced by the cracked interface based on the analytical solution. It is shown that acoustic nonlinearity parameter is a function of two factors, which are the tension/compression elastic asymmetry and the frequency of incident wave, respectively. This study is beneficial to develop nonlinear ultrasonic quantitative NDE techniques for assessing microcrack-induced damage of interfaces.

In future, we will design and develop the reasonable experiments to validate the solution of this paper.

http://dx.doi.org/10.1155/2016/5023127

Competing Interests

The authors declare that there are no competing interests regarding the publication of this paper.

Acknowledgments

Youxuan Zhao acknowledges the support from the Fundamental Research Funds for the Central Universities (Project no. 0903005203383). Yewang Su acknowledges the support from Chinese Academy of Sciences via the "Hundred Talent program" and support from NSFC (no. 11572323).

References

[1] A. J. Croxford, P. D. Wilcox, B. W. Drinkwater, and P. B. Nagy, "The use of non-collinear mixing for nonlinear ultrasonic detection of plasticity and fatigue," The Journal of the Acoustical Society of America, vol. 126, no. 5, pp. EL117-EL122, 2009.

[2] J. Jiao, J. Sun, N. Li, G. Song, B. Wu, and C. He, "Micro-crack detection using a collinear wave mixing technique," NDT and E International, vol. 62, pp. 122-129, 2014.

[3] S. Zhou and Y. Shui, "Nonlinear reflection of bulk acoustic waves at an interface," Journal of Applied Physics, vol. 72, no. 11, pp. 5070-5080, 1992.

[4] G. Shui, Y.-S. Wang, P. Huang, and J. Qu, "Nonlinear ultrasonic evaluation of the fatigue damage of adhesive joints," NDT and E International, vol. 70, pp. 9-15, 2015.

[5] V. V. Kazakov and A. M. Sutin, "Pulsed sounding of cracks with the use of the modulation of ultrasound by vibrations," Acoustical Physics, vol. 47, no. 3, pp. 308-312, 2001.

[6] D. T. Zeitvogel, K. H. Matlack, J.-Y. Kim, L. J. Jacobs, P. M. Singh, and J. Qu, "Characterization of stress corrosion cracking in carbon steel using nonlinear Rayleigh surface waves," NDT and E International, vol. 62, pp. 144-152, 2014.

[7] P. B. Nagy, "Fatigue damage assessment by nonlinear ultrasonic materials characterization," Ultrasonics, vol. 36, no. 1-5, pp. 375-381, 1998.

[8] G. P. M. Fierro, F. Ciampa, D. Ginzburg, E. Onder, and M. Meo, "Nonlinear ultrasound modelling and validation of fatigue damage," Journal of Sound and Vibration, vol. 343, pp. 121-130, 2015.

[9] Y. Zhao, Y. Qiu, L. J. Jacobs, and J. Qu, "A micromechanics model for the acoustic nonlinearity parameter in solids with distributed microcracks," AIP Conference Proceedings, vol. 1706, Article ID 060001, 2016.

[10] Y. Zhao, Y. Qiu, L. J. Jacobs, and J. Qu, "Frequency-dependent tensile and compressive effective moduli of elastic solids with distributed penny-shaped microcracks," Acta Mechanica, vol. 227, no. 2, pp. 399-419, 2016.

[11] Y. Zhao, Y. Qiu, L. J. Jacobs, and J. Qu, "Frequency-dependent tensile and compressive effective moduli of elastic solids with randomly distributed two-dimensional microcracks," Journal of Applied Mechanics, vol. 82, no. 8, Article ID 081006, 2015.

[12] K. Jinno, A. Sugawara, Y. Ohara, and K. Yamanaka, "Analysis on nonlinear ultrasonic images of vertical closed cracks by damped double node model," Materials Transactions, vol. 55, no. 7, pp. 1017-1023, 2014.

[13] J. Riviere, M. C. Remillieux, Y. Ohara et al., "Dynamic acousto-elasticity in a fatigue-cracked sample," Journal of Nondestructive Evaluation, vol. 33, no. 2, pp. 216-225, 2014.

[14] G. Sih and J. Loeber, "Normal compression and radial shear waves scattering at a penny-shaped crack in an elastic solid," The Journal of the Acoustical Society of America, vol. 46, no. 3, pp. 711-721, 1969.

[15] A. K. Mal, "Interaction of elastic waves with a Griffith crack," International Journal of Engineering Science, vol. 8, no. 9, pp. 763-776, 1970.

[16] A. K. Mal, "Interaction of elastic waves with a penny-shaped crack," International Journal of Engineering Science, vol. 8, no. 5, pp. 381-388, 1970.

[17] J. D. Achenbach, Ray Methods for Waves in Elastic Solids: With Applications to Scattering by Cracks, Edited by A. K. Gautesen and H. McMaken, Pitman Advanced Publishing, Boston, Mass, USA, 1982.

[18] J. D. Achenbach and A. N. Norris, "Loss of specular reflection due to nonlinear crack-face interaction," Journal of Nondestructive Evaluation, vol. 3, no. 4, pp. 229-239, 1982.

[19] Y. C. Angel and J. D. Achenbach, "Reflection and transmission of elastic waves by a periodic array of cracks," Journal of Applied Mechanics, vol. 52, no. 1, pp. 33-41, 1985.

[20] Y. C. Angel and J. D. Achenbach, "Reflection and transmission of elastic waves by a periodic array of cracks: oblique incidence," Wave Motion, vol. 7, no. 4, pp. 375-397, 1985.

[21] V. E. Nazarov and A. M. Sutin, "Nonlinear elastic constants of solids with cracks," Journal of the Acoustical Society of America, vol. 102, no. 6, pp. 3349-3354, 1997.

[22] A. S. Eriksson, A. Bostrom, and S. K. Datta, "Ultrasonic wave propagation through a cracked solid," Wave Motion, vol. 22, no. 3, pp. 297-310, 1995.

[23] D. Gross and C. Zhang, "Wave propagation in damaged solids," International Journal of Solids and Structures, vol. 29, no. 14-15, pp. 1763-1779, 1992.

[24] J. M. Richardson, "Harmonic generation at an unbonded interface-I. Planar interface between semi-infinite elastic media," International Journal of Engineering Science, vol. 17, no. 1, pp. 73-85, 1979.

[25] S. Biwa, S. Nakajima, and N. Ohno, "On the acoustic nonlinearity of solid-solid contact with pressure-dependent interface stiffness," Journal of Applied Mechanics, vol. 71, no. 4, pp. 508-515, 2004.

[26] J. D. Achenbach and O. K. Parikh, "Ultrasonic analysis of nonlinear response and strength of adhesive bonds," Journal of Adhesion Science and Technology, vol. 5, no. 8, pp. 601-618, 1991.

Yanjun Qiu, (1) Yilin Zhang, (1) Yewang Su, (2) Peng Cao, (3) and Youxuan Zhao (4)

(1) School of Civil Engineering, Southwest Jiaotong University, Chengdu 610031, China

(2) State Key Laboratory of Nonlinear Mechanics, Institute of Mechanics, Chinese Academy of Sciences, Beijing 100190, China

(3) Department of Hydraulic Engineering, Tsinghua University, Beijing 100084, China

(4) College of Aerospace Engineering, Chongqing University, Chongqing 400044, China

Correspondence should be addressed to Youxuan Zhao; youxuan.zhao@gmail.com

Received 14 December 2015; Revised 5 May 2016; Accepted 15 May 2016

Academic Editor: Roman Lewandowski

Caption: Figure 1: One-dimensional simplistic model.

Caption: Figure 2: Waveforms (a), (b1), and (b2) and frequency spectrum (c) of transmitted wave at location x = [0.sup.+].

Caption: Figure 3: Acoustic nonlinearity parameter versus frequency.

Caption: Figure 4: Acoustic nonlinearity parameter versus parameter [gamma].

Adhesion bonding technology for fiber reinforced composite has been widely applied in aerospace industry. Generally, the cohesive interface is very thin but can provide the powerful strength and resiliency. However, the capability of the cohesive interface could be greatly affected by the imperfect condition, such as debond, fatigue damage, and microcracks. Thus, there is a great need for development of nondestructive testing methods to detect the imperfect condition.

Nonlinear ultrasonic methods have the powerful ability to characterize the material nonlinearity change caused by plasticity [1], imperfect interface [2-4], microcracks [5, 6], and fatigue damage [7,8]. When a time-harmonic longitudinal wave propagates through a cracked solid or interface, it will cause the tension and compression elastic asymmetry; then the waveform will be distorted and higher-order harmonic waves are generated [9-11]. And considerable experimental evidence [12,13] has shown that ultrasonic waves do interact with microcracks in a nonlinear fashion, but those researches may fail when the crack size is only tens of microns, which is one of the initial factors leading to interface degradation.

For the researches on acoustic nonlinearity induced by microcracks [14-16], majority are concentrated either on the scattering of elastic waves by a single crack or an array of cracks [17-20] or on the propagation of elastic waves in a cracked medium [21-23]. The first study on such contact-induced acoustic nonlinearity is probably the paper by Richardson [24] who considered the contact interface between two semi-infinite half-spaces. He has analyzed one-dimensional nonlinear wave propagation in a system composed of an unbounded planer interface separating two semi-infinite linear elastic media. The nonlinearity is caused by the opening and closing of the interface. However, in Richardson's analysis, the interface stiffness varying continuously is not accounted for. Improving Richardson's theory, Biwa et al. [25] have analyzed a nonlinear interface stiffness model, where the stiffness property of the contact interface is described as a function of the nominal contact pressure. Furthermore, by assuming the interface of the adhesive as a nonlinear spring, Achenbach and Parikh [26] have investigated theoretically to obtain information on the adhesive bond strength from ultrasonic test results, and it is shown that the nonlinear adhesive bond behavior could cause the generation of higher harmonics.

In this paper, based on the researches of Richardson and Biwa, one theoretical solution for one-dimensional nonlinear wave propagation in a system composed of a multicrack-included interface separating two semi-infinite elastic media is derived, from which the acoustic nonlinearity caused by the tension and compression elastic asymmetry is clearly shown. To validate the analytical solution, numerical results from the FEM simulations are presented. Comparison between the analytical predictions and the FEM simulation results shows good agreement. Finally, we also study two factors affecting acoustic nonlinearity. One is the tension and compression elastic asymmetry of the cracked interface and another is the frequency of incident wave.

2. Solution for Harmonic Wave Incidence

We consider a system composed of a cracked interface separating two semi-infinite elastic media. In the presence of microcracks, the interface will respond to the tensile and compressive loadings differently, which is the tension and compression elastic asymmetry. According to [11], the elastic constants of the cracked medium can be considered to be dependent with crack density, the friction of the crack faces, and frequency of incident wave:

[mathematical expression not reproducible], (1)

where E' and v' are Young's modulus and Poisson's ratio of the uncracked solid, respectively. And [[bar.E].sub.[alpha]] and [[bar.v].sub.[alpha]] are, respectively, the corresponding effective Young's modulus and Poisson's ratio of the cracked solid under tension ([alpha] = t) and under compression ([alpha] = c). Meanwhile, E' = E and v' = v for plane stress, and E' = E/( 1 - [v.sup.2]) and v = v/(1 - v) for plane strain. The crack density d is defined as d = [Na.sup.2]/A, where a is the average half-length of the cracks; A is the area of the solid containing N randomly distributed and randomly oriented two-dimensional microcracks.

We have

[mathematical expression not reproducible]. (2)

In the above, [gamma] is the Euler-Mascheroni constant, [eta] = [k.sub.T]a = [omega]a/[c.sub.T] is a dimensionless wavenumber, [kappa] = [c.sub.L]/[c.sub.T] = [square root of (2(1 - v)/(1 - 2v))] is a dimensionless parameter, [mu] is the coefficient of friction, and

[mathematical expression not reproducible]. (3)

Therefore, in this paper, the conclusion of [11] will be adopted, and the interface can be regarded as an equivalent medium consisting of the random distribution cracks. Actually this assumption has been proven in [25] for investigating the void inclusion medium subjected to the ultrasonic wave loading conditions.

2.1. The General Solution. Firstly, we assume that the elasticity of the cracked interface is a function of time, and the stress in the interface is uniform. In Figure 1 we illustrate a one-dimensional system schematically to consider elastic longitudinal wave propagation along x-axis, where f(x - ct) and g(x - ct) represent the incident and transmitted wave functions, respectively. The blue area represents the cracked interface.

The two semi-infinite elastic solids are located in the regions x < [0.sup.-] and x > [0.sup.+], respectively. The equation of responses of the interface and the stress-strain relation are defined as follows:

[mathematical expression not reproducible], (4)

where [rho] is mass density, [E.sub.s](t) is the longitudinal stiffness of the cracked interface, u(x, t) is the displacement in the x-direction of an element at time t from its position x, and [sigma](x, t) is the stress. The boundary conditions in locations [0.sup.-] and [0.sup.+] are as follows:

[sigma](t) = [sigma]([0.sup.-], t) = [sigma]([0.sup.+], t), (5)

u([0.sup.+], t) - u([0.sup.-], t) = [sigma](t) x h/[E.sub.s], (6)

where h is the thickness of the cracked interface.

According to [24,25], the governing equation of (4) can be solved based on travelling wave method. Therefore, this paper will employ travelling wave method to obtain the general displacement solution when longitudinal wave propagates through the interface. Then considering the cracked interface as the tension and compression elastic asymmetry, this paper will achieve the corresponding special displacement solution.

According to the travelling wave method, as a solution to (4), the following forms are considered [24,25]:

u(x, t) = f(x - ct) - 1/2V(t + x/c), x < [0.sup.-], (7)

u(x, t) = f(x - ct) + 1/2v(t - x/c), x [greater than or equal to] [0.sup.+], (8)

V(t) = u([0.sup.+], t) - u([0.sup.-], t), (9)

where c = [square root of (E/[rho])] is the wave velocity, E is the longitudinal stiffness of the two semi-infinite elastic solids, f(x - ct) represents the incident waves, the term -(l/2)V(t + x/c) represents reflected waves, and (8) represents the transmitted waves. Equation (9) represents the equivalent displacement of the cracked interface.

According to (6) and (9), we can obtain

[mathematical expression not reproducible]. (10)

The incident wave f(x - ct) is now assumed in the following form:

f(x - ct) = [F.sub.0] cos [k(x - ct)], (11)

where k is the wave number and [F.sub.0] is the maximum amplitude.

Substitution of (11) in (10) yields

dV(t)/dt + [xi]V(t) = 2[F.sub.0][omega] sin ([omega]t), (12)

where [omega] is the angular frequency and [xi] = [xi](t) = 2c[E.sub.s](t)/hE.

Therefore, the general solution of (12) is

[mathematical expression not reproducible]. (13)

When we know the function form of [xi](t), we can obtain the special solution from (13).

2.2. The Special Solution. Assuming the different tensile and compression effective moduli of the cracked interface from [11], the function [E.sub.s](t) can be defined by

[mathematical expression not reproducible], (14)

where [E.sup.t.sub.s] and [E.sup.c.sub.s] are the tensile and compression effective moduli, respectively. And the parameter [gamma] = ([E.sup.t.sub.s] - [E.sup.c.sub.s])/E is introduced to represent the tension and compression asymmetry in the elastic moduli of the cracked interface.

According to (13) and (14), we can obtain

[mathematical expression not reproducible], (15)

where [phi] = [sin.sup.-1] (-[omega]/[square root of ([[xi].sup.2] + [[omega].sup.2]))].

Since V([t.sub.n]) = 0 and the calculation could begin from [t.sub.0] = 0 without loss of generality, the time points [t.sub.n+1] to switch the property come from the following equation:

[mathematical expression not reproducible]. (16)

The displacement becomes

[mathematical expression not reproducible], (17)

where

[mathematical expression not reproducible]. (18)

Substitution of (17) back to (7) and (8) gives the finial displacement field of the problem.

3. Analytical and FEM Results

3.1. Analytical Results. Numerical solutions using Matlab are performed to solve (16). For a demonstrative purpose, dashed line in Figure 2 shows one example of transmitted waveforms for an incident wave with the center frequency of 1 MHz, obtained by numerical solutions of (16) and (17), with the parameters [rho] = 2700 kg/[m.sup.3], E = 7 x [10.sup.10] Pa, [E.sup.t.sub.s] = 7 x [10.sup.9] Pa, [E.sup.c.sub.s] = 7 x [10.sup.10] Pa, [F.sub.0] = 1 x [10.sup.-5] m, and h = 1 x [10.sup.-4] m. The time-domain waveform with the dashed line is shown in Figure 2(a), and the locals in Figures 2(b1) and 2(b2) are to compare the difference of the displacement curve subjected to the first-half and second-half waveform in one-cycle loading conditions. It is seen that the first-half waveform in one cycle is distinctly different from the second-half waveform caused by the tension and compression elastic asymmetry. And Figure 2(c) shows the frequency spectrum using FFT method; it is clearly shown that the transmitted wave contains the term with angular frequency 2[omega], which is second harmonic generated by the cracked interface.

3.2. Comparison with FEM Results. The analytical solutions will be validated in this section by conducting numerical simulations of wave propagation through a cracked interface. The numerical simulations are performed using the finite element method (FEM) for the case of two-dimensional plane strain deformation. The commercial FEM software ABAQUS is used for this purpose.

To this end, a two-dimensional FEM model is constructed using the four-node plane strain (CPE4R) elements. To simulate plane waves, a rectangular strip of 0.09 m x 0.003 m is used in our FEM model. The model consists of ~300,000 four-node plane strain (CPE4R) elements. The user defined constitutive law is performed to present the different tension and compression elasticity of the cracked interface. And the interface is located in the middle of the sample with the thickness of 1 x [10.sup.-4] m. Other regions are the isotropic elastic material of aluminum. A plane longitudinal wave is generated from one end of the sample. The amplitude longitudinal wave with the frequency of 1 MHz is 1 x [10.sup.-5] m. The receiver is located at right end of the interface.

Shown in Figure 2(a) with the solid line is the FEM result of time-domain waveform of transmitted wave. The excellent comparison shows that the analytical solution results are identical to that of the FEM simulation. Furthermore, the two obvious amplitudes, being 9.37 x [10.sup.-6] m and 7.76 x [10.sup.-7] m, respectively, exist with the frequencies of ultrasonic wave at 1 MHz and 2 MHz in Figure 2(c).

4. Acoustic Nonlinearity Parameter

The time-domain waveform curve can be obtained by the analytical solution derived in Section 2; then the frequency spectrum using FFT method can also be solved to obtain [A.sub.1] and [A.sub.2], which are the amplitudes of the fundamental and second harmonic, respectively. Through the analysis of data from the analytical solution, we can obtain several observations. First, [A.sub.2] is linearly related to [A.sub.1]. This is significantly different from the second harmonic generated by quadratic nonlinearity in the solid, where [A.sub.2] is proportional to the square of [A.sub.1]. The fundamental cause of this relationship is the tension and compression elastic asymmetry of the cracked interface. Second, [A.sub.2] is scaled by parameter [gamma]. Therefore, according to [9], for acoustic nonlinearity parameter in cracked solids, the similar acoustic nonlinearity parameter of the cracked interface is introduced by

[beta] = 3[pi][A.sub.2]/[A.sub.1]kh. (19)

For investigating the relation between the linearity parameter and the frequencies of the ultrasonic acoustic, a series of analytic solutions have been performed under different frequencies including 0.5 MHz, 1 MHz, 1.5 MHz, 2.0 MHz, and 2.5 MHz. Figure 3 shows the acoustic nonlinearity parameter versus frequency. It is seen that although the acoustic nonlinearity parameter decreases linearly with increasing frequency, the decrease is not significant. In other words, the acoustic nonlinearity parameter has a rather weak dependence on frequency.

Figure 4 shows the acoustic nonlinearity parameter versus parameter [gamma]. It is seen that the acoustic nonlinearity parameter increases linearly with increasing [gamma], which can represent the degree of damage for the cracked interface. In other words, the degree of damage for the cracked interface can be characterized by the acoustic nonlinearity parameter defined as (19).

5. Conclusions

As a time-harmonic longitudinal wave propagates in a system composed of a cracked interface separating two semi-infinite elastic media, a second harmonic wave may be generated by the cracked interface. Assuming the longitudinal stiffness of the interface as a function of time, this paper solves the problem and obtains the analytical general solution (any interface) and the special solution (the cracked interface). It is clearly shown that the tension and compression elastic asymmetry can result in acoustic nonlinearity. Furthermore, to verify the developed theoretical solution, we carry out detailed numerical simulation of wave propagation through a cracked interface by using the FEM. The results indicate that the micromechanics model predictions agree well with the FEM simulations. Finally, because the amplitude of the second harmonic is linearly related to the amplitude of the fundamental harmonic, we introduce the acoustic nonlinearity parameter induced by the cracked interface based on the analytical solution. It is shown that acoustic nonlinearity parameter is a function of two factors, which are the tension/compression elastic asymmetry and the frequency of incident wave, respectively. This study is beneficial to develop nonlinear ultrasonic quantitative NDE techniques for assessing microcrack-induced damage of interfaces.

In future, we will design and develop the reasonable experiments to validate the solution of this paper.

http://dx.doi.org/10.1155/2016/5023127

Competing Interests

The authors declare that there are no competing interests regarding the publication of this paper.

Acknowledgments

Youxuan Zhao acknowledges the support from the Fundamental Research Funds for the Central Universities (Project no. 0903005203383). Yewang Su acknowledges the support from Chinese Academy of Sciences via the "Hundred Talent program" and support from NSFC (no. 11572323).

References

[1] A. J. Croxford, P. D. Wilcox, B. W. Drinkwater, and P. B. Nagy, "The use of non-collinear mixing for nonlinear ultrasonic detection of plasticity and fatigue," The Journal of the Acoustical Society of America, vol. 126, no. 5, pp. EL117-EL122, 2009.

[2] J. Jiao, J. Sun, N. Li, G. Song, B. Wu, and C. He, "Micro-crack detection using a collinear wave mixing technique," NDT and E International, vol. 62, pp. 122-129, 2014.

[3] S. Zhou and Y. Shui, "Nonlinear reflection of bulk acoustic waves at an interface," Journal of Applied Physics, vol. 72, no. 11, pp. 5070-5080, 1992.

[4] G. Shui, Y.-S. Wang, P. Huang, and J. Qu, "Nonlinear ultrasonic evaluation of the fatigue damage of adhesive joints," NDT and E International, vol. 70, pp. 9-15, 2015.

[5] V. V. Kazakov and A. M. Sutin, "Pulsed sounding of cracks with the use of the modulation of ultrasound by vibrations," Acoustical Physics, vol. 47, no. 3, pp. 308-312, 2001.

[6] D. T. Zeitvogel, K. H. Matlack, J.-Y. Kim, L. J. Jacobs, P. M. Singh, and J. Qu, "Characterization of stress corrosion cracking in carbon steel using nonlinear Rayleigh surface waves," NDT and E International, vol. 62, pp. 144-152, 2014.

[7] P. B. Nagy, "Fatigue damage assessment by nonlinear ultrasonic materials characterization," Ultrasonics, vol. 36, no. 1-5, pp. 375-381, 1998.

[8] G. P. M. Fierro, F. Ciampa, D. Ginzburg, E. Onder, and M. Meo, "Nonlinear ultrasound modelling and validation of fatigue damage," Journal of Sound and Vibration, vol. 343, pp. 121-130, 2015.

[9] Y. Zhao, Y. Qiu, L. J. Jacobs, and J. Qu, "A micromechanics model for the acoustic nonlinearity parameter in solids with distributed microcracks," AIP Conference Proceedings, vol. 1706, Article ID 060001, 2016.

[10] Y. Zhao, Y. Qiu, L. J. Jacobs, and J. Qu, "Frequency-dependent tensile and compressive effective moduli of elastic solids with distributed penny-shaped microcracks," Acta Mechanica, vol. 227, no. 2, pp. 399-419, 2016.

[11] Y. Zhao, Y. Qiu, L. J. Jacobs, and J. Qu, "Frequency-dependent tensile and compressive effective moduli of elastic solids with randomly distributed two-dimensional microcracks," Journal of Applied Mechanics, vol. 82, no. 8, Article ID 081006, 2015.

[12] K. Jinno, A. Sugawara, Y. Ohara, and K. Yamanaka, "Analysis on nonlinear ultrasonic images of vertical closed cracks by damped double node model," Materials Transactions, vol. 55, no. 7, pp. 1017-1023, 2014.

[13] J. Riviere, M. C. Remillieux, Y. Ohara et al., "Dynamic acousto-elasticity in a fatigue-cracked sample," Journal of Nondestructive Evaluation, vol. 33, no. 2, pp. 216-225, 2014.

[14] G. Sih and J. Loeber, "Normal compression and radial shear waves scattering at a penny-shaped crack in an elastic solid," The Journal of the Acoustical Society of America, vol. 46, no. 3, pp. 711-721, 1969.

[15] A. K. Mal, "Interaction of elastic waves with a Griffith crack," International Journal of Engineering Science, vol. 8, no. 9, pp. 763-776, 1970.

[16] A. K. Mal, "Interaction of elastic waves with a penny-shaped crack," International Journal of Engineering Science, vol. 8, no. 5, pp. 381-388, 1970.

[17] J. D. Achenbach, Ray Methods for Waves in Elastic Solids: With Applications to Scattering by Cracks, Edited by A. K. Gautesen and H. McMaken, Pitman Advanced Publishing, Boston, Mass, USA, 1982.

[18] J. D. Achenbach and A. N. Norris, "Loss of specular reflection due to nonlinear crack-face interaction," Journal of Nondestructive Evaluation, vol. 3, no. 4, pp. 229-239, 1982.

[19] Y. C. Angel and J. D. Achenbach, "Reflection and transmission of elastic waves by a periodic array of cracks," Journal of Applied Mechanics, vol. 52, no. 1, pp. 33-41, 1985.

[20] Y. C. Angel and J. D. Achenbach, "Reflection and transmission of elastic waves by a periodic array of cracks: oblique incidence," Wave Motion, vol. 7, no. 4, pp. 375-397, 1985.

[21] V. E. Nazarov and A. M. Sutin, "Nonlinear elastic constants of solids with cracks," Journal of the Acoustical Society of America, vol. 102, no. 6, pp. 3349-3354, 1997.

[22] A. S. Eriksson, A. Bostrom, and S. K. Datta, "Ultrasonic wave propagation through a cracked solid," Wave Motion, vol. 22, no. 3, pp. 297-310, 1995.

[23] D. Gross and C. Zhang, "Wave propagation in damaged solids," International Journal of Solids and Structures, vol. 29, no. 14-15, pp. 1763-1779, 1992.

[24] J. M. Richardson, "Harmonic generation at an unbonded interface-I. Planar interface between semi-infinite elastic media," International Journal of Engineering Science, vol. 17, no. 1, pp. 73-85, 1979.

[25] S. Biwa, S. Nakajima, and N. Ohno, "On the acoustic nonlinearity of solid-solid contact with pressure-dependent interface stiffness," Journal of Applied Mechanics, vol. 71, no. 4, pp. 508-515, 2004.

[26] J. D. Achenbach and O. K. Parikh, "Ultrasonic analysis of nonlinear response and strength of adhesive bonds," Journal of Adhesion Science and Technology, vol. 5, no. 8, pp. 601-618, 1991.

Yanjun Qiu, (1) Yilin Zhang, (1) Yewang Su, (2) Peng Cao, (3) and Youxuan Zhao (4)

(1) School of Civil Engineering, Southwest Jiaotong University, Chengdu 610031, China

(2) State Key Laboratory of Nonlinear Mechanics, Institute of Mechanics, Chinese Academy of Sciences, Beijing 100190, China

(3) Department of Hydraulic Engineering, Tsinghua University, Beijing 100084, China

(4) College of Aerospace Engineering, Chongqing University, Chongqing 400044, China

Correspondence should be addressed to Youxuan Zhao; youxuan.zhao@gmail.com

Received 14 December 2015; Revised 5 May 2016; Accepted 15 May 2016

Academic Editor: Roman Lewandowski

Caption: Figure 1: One-dimensional simplistic model.

Caption: Figure 2: Waveforms (a), (b1), and (b2) and frequency spectrum (c) of transmitted wave at location x = [0.sup.+].

Caption: Figure 3: Acoustic nonlinearity parameter versus frequency.

Caption: Figure 4: Acoustic nonlinearity parameter versus parameter [gamma].

Printer friendly Cite/link Email Feedback | |

Title Annotation: | Research Article |
---|---|

Author: | Qiu, Yanjun; Zhang, Yilin; Su, Yewang; Cao, Peng; Zhao, Youxuan |

Publication: | Mathematical Problems in Engineering |

Date: | Jan 1, 2016 |

Words: | 3431 |

Previous Article: | A Mathematical Method for Eliminating Spin Losses in Toroidal Traction Drives. |

Next Article: | Research on Sustainable Development of Resource-Based Cities Based on the DEA Approach: A Case Study of Jiaozuo, China. |

Topics: |