Abstract

Second harmonic generation (SHG) is one of the common techniques in the nonlinear ultrasonic test. The contact-type defects play an important role in material damage, which are hard to be detected. The traditional nonlinear parameter β used to evaluate the micro damage in material is derived from the classical stress-strain relation, which is more suitable for the anharmonicity of crystal rather than the contact-type defects. Recently, the theoretical model based on the bilinear stiffness law was derived, and the validity and applicability need to be further studied. For this purpose, by the numerical method, the contact interface in mortar is characterized based on the damage indicator γ. The relation between the excitation voltage and γ is obtained. Moreover, the effects of the crack length and orientation on the damage indicator γ are also obtained. The experimental method is also used to characterize the contact interface in mortar. Combining with the existing work, the results obtained in this article are discussed, and further conclusions can be drawn. The conclusions in this article provide potential of quantitative detection of the contact interface and quality evaluation of bonding layers in materials.

1. Introduction

The generation and development of the micro cracks lead to stress concentration, playing an important role in material deterioration and failure [1]. For decades, nonlinear ultrasonic methods have gained significant attention, and they offer enormous potential for detecting and evaluating the change in the microstructure of the materials [26]. A lot of studies have proved that traditional linear ultrasonic methods are infeasible for the micro damage detection in a material [68]. However, in this stage, the mechanism and regularity of the contact acoustic nonlinear behavior caused by interface are still insufficient, which should be studied more comprehensively.

Experimental studies on the higher harmonic generation (HHG) technique have been extensively conducted [9, 10], and rich results have been obtained for different damage types (crystal dislocation, precipitation, micro cracks, etc.) and materials (metals, rock, concrete, etc.) [1114]. Initially, most studies focused on the cubic nonlinearity, which belongs to the field of classical nonlinear ultrasound. For two decades, it has been demonstrated that the nonlinearity generated by contact-type defects at the mesoscale can be orders of magnitude higher than that of the microscale [15, 16]. This provides the possibility of damage detection for the contact interfaces. In subsequent studies, the HHG technique for contact interface detection has been extensively studied [1719], and the nonlinear parameter β is generally used to characterize the damage. However, there are still many barriers before applying this technique in practical test. For example, there are a lot of theoretical studies on the SHG at the contact-type defects, while the validation by the experiment is still insufficient. In addition, especially for geomaterial and concrete, the contact-type defects are extensively and randomly distributed, while the unified baseline of the damage state of the materials is hard to determine [2022]. Thus, the results cannot be compared with each other to obtain further quantitative conclusions [2325].

Theoretically, three crack-wave interaction modelling approaches (i.e., bilinear stiffness, hysteresis, and rough surface contact) are commonly used to explain the HHG phenomenon, which have been comprehensively reviewed by Broda et al. [26]. The three models have their own advantages and disadvantages. The assumption of the bilinear stiffness is hard to be satisfied, unless the level of excitation voltage is high enough to drive the cyclic open-closed behavior of the surfaces. The rough surface contact is too complex when dealing with solids with distributed microcracks. Hysteresis is a phenomenological model, and the physical structure of the crack is not taken into consideration. Recently, the scattering effect accompanied with the nonlinear ultrasonic phenomenon generating at the contact surface has attracted much attention [2730]. Based on the rough surface contact, the nonlinear ultrasonic effects generated at the interface is studied, and the results have been verified through experiment [31]. In addition, a novel approach is proposed to analyze the scattered field at a plane interface based on the rough surface model. The corresponding solution is obtained by decomposing the scattered field into two components. The novel approach is simpler for the less unknowns [30]. Furthermore, the explicit analytical solution of the scattered field is derived, improving the previous theoretical studies obtained by numerical method [32]. In addition, Zhao et al. derived the nonclassical nonlinear parameter based on a bilinear stiffness [33, 34]. However, the proposed model is still insufficient for lack of experimental validation. Besides, β0 is a global indicator, which provides no information on the local damage in materials.

In actual test, acoustic nonlinearity is commonly introduced by the electronic system as well as the coupling between the transductors and specimens. Thus, it is difficult to distinguish the material nonlinearity from the obtained nonlinearity in experiment. In general, the numerical method can be considered as an effective method for the nonlinear ultrasonic study [3539]. By finite element simulations, the motion of interface obeys the contact laws, and the nonlinear effects of ultrasonic waves can be obtained. By FEM analysis, the nonlinear ultrasonic technique is beneficial in detecting a closed crack with a different orientation [40]. A finite difference time-domain technique is proposed to study the wave scattering caused by cracks with interacting faces, and the accuracy and stability of the scheme in one-dimensional (1D) and two-dimensional (2D) spaces are verified [41]. However, numerical simulations are difficult to consider all factors having effects on the results. Thus, considerable difference usually exists between the numerical results and the actual test, and the experiment method is necessary to be a supplement.

The objective of this study is to detect the contact interface in materials based on the bilinear stiffness assumption. For this purpose, the damage indicator derived from the bilinear stiffness laws is verified numerically. The effects of excitation voltage on the indicator γ are investigated. Furthermore, the effects of the interface length and orientation on the second harmonic generation are also studied. The experimental study is conducted in mortar specimens, and the regularities are coincident with those of numerical results.

The paper is structured as follows. In Section 2, the theoretical background of the nonclassical damage indicator γ is briefly introduced. Thereafter, the finite element modelling is described in Section 3, and the numerical results are given and discussed. In Section 4, the nonlinear ultrasonic test is carried out in mortar specimens and the results are presented and discussed. Finally, the conclusions of the study are presented.

2. SHG at Contact Surface

The explicit expression of SHG can be derived based on the bilinear stiffness assumption. Considering a longitudinal plane wave propagation in an isotropic and a linear elastic solid with Young’s modulus E0, the corresponding wave equation is given by

The corresponding bilinear stiffness expression can be expressed aswhere is the initial static contact strain, is the Heaviside step function, and coefficient is the stiffness weakened by the crack interface. It is assumed that for the tensile strain. is given by the following equation:

Substituting the bilinear stiffness expression into the 1D wave equation, the inhomogeneous wave equation can be obtained by

Assuming that a sine wave passes through the crack at a normal direction, the corresponding boundary-initial value problem is calculated using the following equation:

By solving the above equation, the amplitude of second harmonic A2 can be expressed aswhere is the angular frequency, is the propagation distance, and is the longitudinal velocity.

The nonclassical damage indicator γ based on the bilinear stiffness law can be written as

It should be noted that equations (6) and (7) are established on the assumption of at the interface, and the tensile part of the cyclic load cannot pass through the surface in this case.

In addition, several conclusions can be drawn according to the expression of the damage indicator γ, which are as follows.(1)The damage indicator γ is linear with the wave propagation distance x from the nonlinear source, which is the same as the classical nonlinear parameter β. The expression of β is demonstrated in the following equation:where k is the wavenumber of the incident waves.(2)The damage indicator γ is defined as the ratio of the second harmonic amplitude A2 over the fundamental A1, and the validity and reliability should be further examined.

3. Finite Element Model

3.1. Model Description

By defining the bilinear stiffness law at the surface, the CAN generated at the surface can be observed. In addition, the regularities of the crack length as well as the orientation on the damage indicator γ can also be obtained. The commercial FEM software ABAQUS is used for these purposes.

Herein, a homogeneous and isotropic solid of 4 × 4 mm2 in a 2D numerical simulation is considered. A contact interface of length L0 is introduced at the center of the domain as shown in Figure 1. The cyclic load excited at the top boundary of the model is a 7-cycle Hann-windowed longitudinal pulse with the frequency of 500 kHz. The node located at the middle of the bottom boundary is selected to collect the y-direction displacements. The material properties are given in Table 1. The mesh of the model has a size of 0.03 mm. This ensures that the considered wavelength is discretized by enough elements. The element type is a 4-node bilinear plane strain quadrilateral reduced integration hourglass control (CPE4R). The mesh quality has been checked, and no poor mesh exists in the model. The time increment is 1.6 ns, which satisfies the Courant-Friedrichs-Lewy stability condition, ( is the smallest element dimension), and Nyquist-Shannon sampling theorem, ( is the sampling frequency and is the maximum frequency). The total time period is 3 × 10−5 s. Fast Fourier transform (FFT) is performed to obtain the amplitudes of the fundamental A1 and second harmonic A2. To eliminate the nonlinearity coming from the large displacement and rotation of the elements, the geometrical nonlinearity is set to be off, and both the linear and the quadratic bulk viscosity parameters are zero. Therefore, it can be concluded that the detected nonlinearity is solely generated by the contact behavior of the interface surfaces.

3.2. Contact Laws

The bilinear stress-strain law is assumed to be concerned with the stiffness asymmetry at the crack surfaces, which can be stated as

It is assumed that is positive for the tensile strain. The weakened elastic modulus is dependent on the damage degree in material. When the interface is well bonded (i.e., no damage exists in material), equals 0 in the cyclic compression-tensile load. When the interface is ruptured, will increase to , and the tensile components of the cyclic load are unable to pass through the surface. In this numerical model, the contact laws applied at the crack can be stated as

The above contact laws are coincident with the theoretical background. The prestress is assumed to be zero in the model. In addition, it should be noted that the second harmonic generation depends on the asymmetrical normal stiffness. For consistency, the surface is considered flat and frictionless, and the vertical displacement is received and analyzed. In addition, the contact type of the interface is “hard contact.” This indicates that the surfaces cannot interpenetrate into each other, and the crack surface allows separation after the contact.

The time and frequency domain spectra of the intact and cracked models are shown in Figure 2. The amplitude of the excitation is 80 nm in this case, and the length of the crack is 1 mm. The elastic waves are normally incident at the interface. It can be concluded that, for the cracked sample, the +y direction displacement of the bottom element of the surface is obviously small. This indicates that the tensile strain is hard to pass through the unbonded interface. It can be observed from the frequency spectrum that significant second harmonic appears at the surface of the crack compared with the intact sample. As mentioned above, the finite element model has no geometric and material nonlinearity, and the acoustic nonlinearity should come from the interaction between the surfaces of the crack. The above results further prove that the contact laws are effective, which can be further applied in the numerical study.

3.3. Numerical Results
3.3.1. Effects of the Excitation Levels on γ

As introduced in the theoretical background, the damage indicator γ is defined as the ratio of the second harmonic amplitude A2 to the fundamental A1, and the damage indicator γ is a constant with increase of the excitation levels. The FEM analysis will be conducted to verify the analytical expression.

The amplitude of the incident waves varies from 10 to 150 nm with a step of 10 nm. In addition, the crack has a length of 1 mm and is normal to the incident waves. The damage indicator γ for various excitation levels is presented in Figure 3. It can be noted that when the amplitude of incident wave is weak, corresponding γ is small. With the amplitude increasing, γ increases sharply. Then, γ remains stable in case of large enough amplitude applied. It is known that the CAN effects will occur once the excitation level becomes greater than the threshold. The strong cyclic load near the surface will lead the nonlinear effects. Considering that the trend of γ is in accordance with the CAN model, the correctness and validity of γ can be examined by the FEM model.

3.3.2. Effects of the Interface Length on the Damage Indicator γ

The length of contact interface is characterized based on the damage indicator γ. The crack length L0 is arranged from 0.4 to 1 mm with an interval of 0.1 mm. The incident waves are normally excited at the top surface of the model. The damage indicator γ is determined by the linear fitting of A1 and A2. The amplitude of the incident waves varies from 60 to 100 nm with a step of 10 nm. The damage indicator γ for various crack lengths is shown in Figure 4. As can be seen in Table 2, the coefficient of determination R2 of each fit line is nearly 1, which indicates that the excitation level is independent of the damage indicator γ. In addition, Figure 5 shows that the crack length has a significant effect on the damage indicator γ. A larger damage indicator γ will be obtained with increase of the interface length. A significant positive correlation between the damage indicator γ and the crack length L0 can be obtained.

It is worth mentioning that the crack length is also characterized by other damage indicators based on the nonlinear ultrasonic method [34], as shown in Figure 6. It can be concluded that the nonlinear effects are stronger when the crack becomes longer. Thus, the FEM results provide potential of quantitively characterizing the contact interface in materials.

3.3.3. Effects of the Interface Angle on the Damage Indicator γ

To understand the regularity of the crack orientation on the damage indicator γ, seven crack orientations (α = 0°, 15°, 30°, 45°, 60°, 75°, and 90°) are considered. The crack length is 1 mm in this section. Besides, the amplitude of the incident waves varies from 60 to 100 nm with a step of 10 nm. Similarly, the damage indicator γ is determined by linear fitting of A1 and A2. With increase of the excitation level, the distribution of indicator γ can be observed in Figure 7. Figure 8 shows the damage indicator γ versus the crack orientation. It can be concluded that the damage indicator γ decreases as the crack angle increases. According to the bilinear stiffness assumption, the second harmonic is generated by the clapping effect of the surface. Thus, the significant contact nonlinearity should appear in case of the waves being normally incident at the surface.

In other studies, the similar regularity of the crack orientation on the second harmonic generation has been obtained [35]. As shown in Figure 9, in the case of P wave, normalized A1 and A2 are affected by the crack angle. Fundamental A1 increases with the interface angle, while the opposite trend can be observed on the second harmonic A2. Thus, the regularity of the damage indicator γ decreases with increase of the interface angle, which is consistent with the numerical results in this paper.

4. Experimental Procedure

4.1. Sample Preparation

The size of the mortar specimens is 100 mm × 100 mm × 100 mm. Three groups of cubic mortar specimens were prepared in this experiment. Group A contains 3 intact mortar specimens (I1∼I3), which is used to study the relation between the excitation voltage and the damage indicator γ. Group B contains 6 specimens (L1∼L6). The length of the interface is 1∼6 cm, with an interval of 1 cm. The interface angle is 0° in this group. Group C contains 7 specimens (A1∼A7). The crack orientation is 0∼90° with a step of 15°. The interface length is 4 cm. The introduction procedure of the contact interface will be given later. The mixture proportions and physical properties of the mortar specimens are given in Table 3. The Portland cement (P.O. 32.5) was used in this experiment.

The procedure for the contact interface generation is as follows. After pouring of mortar, a thin steel sheet with machine oil on the surface vertically penetrated into the mortar mixture. All the steel sheets are 12 cm in length and 1 mm in thickness. The insertion depth is 100 mm. The width of the steel sheets is 1∼6 cm. Then, the steel sheet was removed out of the specimen carefully, and the specimens were not moved before being demolded. It should be noted that the removal of the steel sheets was gradually conducted once sheets arrived at the bottom of the sample, and the whole process of interface introduction took about one minute. Consequently, an artificial contact interface was left in the specimen due to the fluidity of the mixture. The specimens were demolded after 24 h followed by curing in a standard chamber (95% relative humidity and 20°C) until 28 days of age. The mortar specimens are shown in Figure 10.

To confirm that the artificial interface is contact, the pulse velocity of each specimen is tested, and the related results are shown in Figure 11. It is concluded that the wave velocity remains steady in all specimens, which indicates that the contact interface is at micro scale and cannot be detected by the traditional ultrasonic method. Moreover, the specimen A7 was cut in two halves after the test, as shown in Figure 12. The trace along the height of specimen can be observed. It is noted that the damage degree gradually weakened along the height of specimen because of the self-weight of the mixture. In conclusion, the interface within the specimen can be considered contact, and the damage is essentially at micro scale.

It should be pointed out that equation (8) is established on the assumption of the surface rupture, while the interface within the mortar can be considered the stiffness reduction, and the corresponding analytical expression can be given as [42]where γt and γc represent the tension and compression asymmetry in the elastic modulus, respectively. The same damage indicator can be obtained in equations (8) and (11), which is defined as the ratio of A2 to A1.

4.2. Experimental Setup

The experimental setup for the measurement of fundamental A1 and second harmonic A2 is shown in Figure 13. According to the previous study [43], for mortar specimens, strong attenuation will happen when the center frequency is above 500 kHz. Therefore, a tone-burst signal of 10 cycles at 200 kHz is generated by a function generator (Rigol 1022U). The pulse repetition rate is 10 ms. Two contact-type transducers (PXR 50) are used at the transmitting and receiving ends. Figure 14 shows that PXR 50 has a good response to fundamental A1 and second harmonic A2 considering that the sensitivity at f0 and 2f0 is above 60 dB. The receiving time domain signal after modulating is recorded by a digital oscilloscope (Tektronix MDO 3104), and the sampling rate is 50 ms/s. The quality and repeatability of the measurements are improved by averaging the signals with 256 acquisitions. In addition, a thin layer of Vaseline was used as couplant between the transducers and the specimen surface. The surface of the specimens was wiped carefully to ensure no sand or dust is attached on the test point. The specimen and transducers were tightly wrapped during the test. To ensure that the transmitter and the receiver are aligned, the test points were marked before the fixation of transducers. The FFT transform was performed on the time spectrum to obtain fundamental A1 and second harmonic A2. Figure 15 shows an example of the time and frequency spectra measured from a mortar specimen with interface, respectively. Fundamental A1 and second harmonic A2 can be observed from the frequency spectrum, and the damage indicator γ can be calculated.

4.3. Experimental Results and Discussion
4.3.1. Effects of the Excitation Levels on γ

The results of the intact mortar specimens are shown in Figure 16. The value of γ is independent of the excitation voltage. In addition, the damage indicator γ in specimens I1∼I3 can be determined by the slope of the fit line, being 0.050, 0.043, and 0.052, respectively. The damage indicator γ in specimens I1∼I3 can be regarded as the intrinsic nonlinearity of the mortar, which can be used to distinguish the nonlinearity coming from the contact interface in group B and group C.

The value of γ is independent of the excitation voltage, which has been observed by analyzing the leaky Rayleigh wave passing through the crack surface [44]. As shown in Figure 17, the ratio of second harmonic A2 to fundamental A1 is almost a constant when σ is 37.1 kPa, which indicates that the damage indicator γ is independent of the excitation voltage. It should be noted that when a large σ is applied, a minor gap exists between the two blocks. In view of the damage type being the same in the two studies, the results are consistent with each other even if the materials used in the studies are different. Besides, the mechanism of second harmonic generation is not explained theoretically in Vergara et al.’s study.

4.3.2. Effects of the Interface Length on the Damage Indicator γ

The distribution of damage indicator γ in different interface lengths is shown in Figure 18. γ0 and δ represent the means and standard deviations of the damage indicator γ in intact mortar specimens. γ0 ± δ is used to distinguish which damage indicators of the mortar specimens with interface behave different from the intact specimens. It can be obtained that the damage indicator γ is obviously larger than those of intact specimens for the 2–6 cm interface length. For the 1 cm interface, although the damage indicator γ is in the band, it is still larger than the means of γ in intact mortar specimens. Thus, it can be concluded that the damage indicator γ is effective to characterize the contact interface in materials. On the other hand, positive correlation between the nonlinear parameter γ and the interface length L can also be revealed, which is coincident with the numerical results.

4.3.3. Effects of the Interface Angle on the Damage Indicator γ

Figure 19 shows the distribution of the damage indicator γ in the mortar specimens with different interface orientations. The meanings of γ0 and δ are the same as in the previous section. In general, the damage indicator γ is significantly affected by the interface angle. The decreasing trend can be observed obviously with the increase of the interface angle, which is consistent with the numerical results. When the interface degree is close to 90°, the difference between the intact and the interface becomes smaller. In actual applications, it is reasonable to conduct the test at different direction of the material to avoid missing the interface paralleling the propagation direction of ultrasound.

5. Conclusions

In this paper, the SHG technique is studied based on the bilinear stiffness assumption. The main conclusions drawn in this paper can be stated as follows:(1)By the numerical and experimental method, the regularity of the excitation levels on the damage indicator γ is attained. Both numerical and experimental results show that γ is a constant when the amplitude is large. In addition, γ is effective in characterizing the contact-type defects in mortar.(2)The regularity of the interface length on the indicator γ is also studied. The damage indicator γ becomes larger when the elastic waves are normally incident at the surface. It should be noted that the value of γ is also affected by the attenuation at the surface, and the absorption and scattering effects should be considered to obtain the quantitative results.(3)The damage indicator γ is effective in characterizing the interface orientation. γ decreases with the increase of the interface angle. Moreover, the results are discussed with the existing work and can be further verified and analyzed.

Further in-depth research should be carried out on the experimental studies of the CAN generation in materials, especially the validation for the existing theoretical results. In addition, quantitative characterization should be further studied, which will be of benefit to applying this technique in practical test.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This research was funded by Chongqing Research Program of Basic Research and Frontier Technology, Grant no. cstc2017jcyjBX0066; Key Laboratory of Hydraulic and Waterway Engineering of the Ministry of Education, Grant no. SLK2017B05; Innovative Research Program of Chongqing Jiaotong University, Grant no. 2019B0101; Jiangxi Youth Science Fund Project, Grant no. 20171BAB216042; Advantage Technology Innovation Team of Jiangxi Province (Grant no. 20171BCB24012); and Science and Technology Projects of Department of Transportation of Jiangxi Province, Grant no. 2018Q0028.