Table of Contents Author Guidelines Submit a Manuscript
Mathematical Problems in Engineering
Volume 2016, Article ID 5023127, 7 pages
http://dx.doi.org/10.1155/2016/5023127
Research Article

Ultrasonic Nonlinearity Evaluation of the Cracked Interface

1School of Civil Engineering, Southwest Jiaotong University, Chengdu 610031, China
2State Key Laboratory of Nonlinear Mechanics, Institute of Mechanics, Chinese Academy of Sciences, Beijing 100190, China
3Department of Hydraulic Engineering, Tsinghua University, Beijing 100084, China
4College of Aerospace Engineering, Chongqing University, Chongqing 400044, China

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

Academic Editor: Roman Lewandowski

Copyright © 2016 Yanjun Qiu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

This paper derives a novel analytical solution for acoustic nonlinearity evaluation of the cracked interface. When microcracks exist at the interface, the tensile and compressive effective moduli of the cracked interface are considered to be different. It is clearly shown that the tension and compression elastic asymmetry can result in acoustic nonlinearity. In addition, numerical simulations using the finite element method are conducted to validate the theoretical solution. It is shown that numerical results agree well with the analytical solution. Finally, two factors affecting acoustic nonlinearity are studied based on the analytical solution. One is the tension and compression elastic asymmetry and another is the frequency of incident wave. Different from acoustic nonlinearity parameter of the general second harmonics, it is found that acoustic nonlinearity parameter is a function of two factors.

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 [24], 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 [911]. 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 [1416], majority are concentrated either on the scattering of elastic waves by a single crack or an array of cracks [1720] or on the propagation of elastic waves in a cracked medium [2123]. 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:where and are Young’s modulus and Poisson’s ratio of the uncracked solid, respectively. And and are, respectively, the corresponding effective Young’s modulus and Poisson’s ratio of the cracked solid under tension () and under compression (). Meanwhile, and for plane stress, and and for plane strain. The crack density is defined as , where is the average half-length of the cracks; is the area of the solid containing randomly distributed and randomly oriented two-dimensional microcracks.

We haveIn the above, is the Euler-Mascheroni constant, is a dimensionless wavenumber, is a dimensionless parameter, is the coefficient of friction, and

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 -axis, where and represent the incident and transmitted wave functions, respectively. The blue area represents the cracked interface.

Figure 1: One-dimensional simplistic model.

The two semi-infinite elastic solids are located in the regions and , respectively. The equation of responses of the interface and the stress-strain relation are defined as follows:where is mass density, is the longitudinal stiffness of the cracked interface, is the displacement in the -direction of an element at time from its position , and is the stress. The boundary conditions in locations and are as follows:where 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]:where is the wave velocity, is the longitudinal stiffness of the two semi-infinite elastic solids, represents the incident waves, the term 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

The incident wave is now assumed in the following form: where is the wave number and is the maximum amplitude.

Substitution of (11) in (10) yieldswhere is the angular frequency and .

Therefore, the general solution of (12) is

When we know the function form of , 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 can be defined bywhere and are the tensile and compression effective moduli, respectively. And the parameter is introduced to represent the tension and compression asymmetry in the elastic moduli of the cracked interface.

According to (13) and (14), we can obtainwhere .

Since and the calculation could begin from without loss of generality, the time points to switch the property come from the following equation:

The displacement becomeswhere

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  kg/m3,  Pa,  Pa,  Pa,  m, and  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 , which is second harmonic generated by the cracked interface.

Figure 2: Waveforms (a), (b1), and (b2) and frequency spectrum (c) of transmitted wave at location .
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 × 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  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  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  m and  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 and , 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, is linearly related to . This is significantly different from the second harmonic generated by quadratic nonlinearity in the solid, where is proportional to the square of . The fundamental cause of this relationship is the tension and compression elastic asymmetry of the cracked interface. Second, is scaled by parameter . Therefore, according to [9], for acoustic nonlinearity parameter in cracked solids, the similar acoustic nonlinearity parameter of the cracked interface is introduced by

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 3: Acoustic nonlinearity parameter versus frequency.

Figure 4 shows the acoustic nonlinearity parameter versus parameter . It is seen that the acoustic nonlinearity parameter increases linearly with increasing , 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).

Figure 4: Acoustic nonlinearity parameter versus parameter .

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.

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. View at Publisher · View at Google Scholar · View at Scopus
  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. View at Publisher · View at Google Scholar · View at Scopus
  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. View at Publisher · View at Google Scholar · View at Scopus
  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. View at Publisher · View at Google Scholar · View at Scopus
  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. View at Publisher · View at Google Scholar · View at Scopus
  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. View at Publisher · View at Google Scholar · View at Scopus
  7. P. B. Nagy, “Fatigue damage assessment by nonlinear ultrasonic materials characterization,” Ultrasonics, vol. 36, no. 1-5, pp. 375–381, 1998. View at Publisher · View at Google Scholar · View at Scopus
  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. View at Publisher · View at Google Scholar · View at Scopus
  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. View at Google Scholar
  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. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus
  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. View at Publisher · View at Google Scholar · View at Scopus
  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. View at Publisher · View at Google Scholar · View at Scopus
  13. J. Rivière, 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. View at Publisher · View at Google Scholar · View at Scopus
  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. View at Publisher · View at Google Scholar · View at Scopus
  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. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus
  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. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus
  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. View at Publisher · View at Google Scholar · View at Scopus
  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. View at Publisher · View at Google Scholar · View at Scopus
  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. View at Publisher · View at Google Scholar · View at Scopus
  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. View at Publisher · View at Google Scholar · View at Scopus
  22. A. S. Eriksson, A. Boström, and S. K. Datta, “Ultrasonic wave propagation through a cracked solid,” Wave Motion, vol. 22, no. 3, pp. 297–310, 1995. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus
  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. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus
  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. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus
  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. View at Publisher · View at Google Scholar · View at Scopus
  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. View at Publisher · View at Google Scholar · View at Scopus