Nonlinear Acoustic Modeling and Measurements during the Fatigue Process in Metals

The nonlinear spring model combined with dislocation dipole theory was applied to describe the acoustic nonlinearity during the fatigue process in metals. The spring stiffness changes with fatigue degree. For the early stage, spring stiffness approaches infinity, and the heavier nonlinearity mainly results from the increase of dislocation density. Further fatigue leads to the occurrence of micro-cracks, during which spring stiffness begins to decrease. Abundant micro-crack sprouting accelerates the crack’s expansion, and spring stiffness drops quickly, which causes the obvious decline in the transmitted harmonic amplitudes. Solutions obtained from the nonlinear wave equation with dislocation terms were added into the spring model. Varying spring stiffness was chosen for simulating the fatigue process. Then, nonlinear harmonic variation during this process was observed, which was classified into three stages: (I) the early dislocation fatigue stage; (II) the micro-crack sprouting stage; (III) the crack expansion stage. Nonlinear acoustic measurements were carried out on an aluminum alloy specimen during its fatigue process until cracks could be seen clearly. Harmonic variations in experiments can also be classified into the same three stages as the numerical results, which provides a theoretical and experimental reference for fatigue evaluation in metals using the nonlinear acoustic method.


Introduction
Metal materials are widely used in practical engineering applications, but long-term in-service work makes them prone to failure due to fatigue damages. The linear acoustic method is invalid for detecting fatigue damages of materials, while the finite amplitude sound wave can produce obvious waveform distortion when it propagates through a medium with fatigue damages, which has been proven to be an effective method for evaluating fatigue damages in metal materials [1][2][3][4][5].
Nonlinear acoustic theory in solids established by researchers such as Landau, Murnaghan, and Goldberg provides the basis for nonlinear ultrasonic testing [6,7]. Based on this, theoretical models were proposed to explain the mechanism of the nonlinear acoustic effect caused by micro-damages. For early stage fatigue damages are usually described by the dislocation string [8] or dislocation dipole models [9,10], while these are not suitable for heavier damages owing to the sprouting of micro-cracks. The PM (proposed by Preisach and Mayergoyz) model [11] and contact interface model [12][13][14][15] were proposed to explain the nonlinear acoustic mechanism caused by micro-cracks or imperfect interfaces, through which nonlinear effects such as harmonics, difference frequency waves, resonance frequency drift, and slow dynamics can be simulated. Most of the research works have been based on the above theoretical models. For example, reflection and transmission theories at an open or closed interface were established to study the contact nonlinearity of acoustics based on the contact interface model [16,17]. Some experiments [9,10,18,19] were carried out to observe the nonlinear properties of material delamination and microstructure changes caused by fatigue damages, which established the connection between micro-damages and acoustical nonlinearity. Nonlinear Lamb waves have been commonly applied to evaluate the fatigue damages or detect the bonding quality of solid boards [20,21]. Furthermore, some researchers [22] proved the effectiveness of the nonlinear acoustic method in estimating the strength degradation caused by initial damage of concrete materials. High temperature also leads to micro-damages, which can be measured through the ultrasonic nonlinear parameters to evaluate the microstructure changes during the high-temperature degradation process [23]. The spring model is usually used to study the characteristics of nonlinear acoustic waves at the solid rough contact interface to evaluate the interface boundary conditions [14].
Nonlinear acoustics are used in many research works for the non-destructive evaluation of metal materials, but the above research works are mainly focused on one of the stages during the fatigue process, such as early stage fatigue evaluation or crack detection. Actually, fatigue is a continuously changing process, so modeling the whole process is essential for its evaluation. In our work, the nonlinear spring model combined with dislocation dipoles theory was established for modeling the continuous fatigue process from the early dislocation stage to the crack expansion stage. Solutions obtained from the nonlinear wave equation with dislocations [24][25][26][27] were added into the nonlinear spring model. For the early stage, spring stiffness approaches infinity. Fatigue damage is mainly caused by dislocations in this stage. For the late stage damages, varying stiffness was chosen for modeling the different fatigue degree with the sprouting of micro-cracks. Once the crack begins to expand, spring stiffness declines quickly, and the amplitude of each order harmonic drops obviously. The nonlinear acoustic experimental system was built, and measurements were carried out during the whole fatigue process of the sample, which certified the numerical models.

Nonlinearity Due to Lattice Dislocations at Early Stage Fatigue
Early stage fatigue damage of metals mainly originates from dislocations. As the degree of fatigue further deepens, the dislocation density increases progressively, and the nonlinear acoustic effects becomes more obvious. The dislocation can be described by the dipole model [10], as shown in Figure 1. open or closed interface were established to study the contact nonlinearity of acoustics based on the contact interface model [16,17]. Some experiments [9,10,18,19] were carried out to observe the nonlinear properties of material delamination and microstructure changes caused by fatigue damages, which established the connection between micro-damages and acoustical nonlinearity. Nonlinear Lamb waves have been commonly applied to evaluate the fatigue damages or detect the bonding quality of solid boards [20,21]. Furthermore, some researchers [22] proved the effectiveness of the nonlinear acoustic method in estimating the strength degradation caused by initial damage of concrete materials. High temperature also leads to micro-damages, which can be measured through the ultrasonic nonlinear parameters to evaluate the microstructure changes during the hightemperature degradation process [23]. The spring model is usually used to study the characteristics of nonlinear acoustic waves at the solid rough contact interface to evaluate the interface boundary conditions [14]. Nonlinear acoustics are used in many research works for the non-destructive evaluation of metal materials, but the above research works are mainly focused on one of the stages during the fatigue process, such as early stage fatigue evaluation or crack detection. Actually, fatigue is a continuously changing process, so modeling the whole process is essential for its evaluation. In our work, the nonlinear spring model combined with dislocation dipoles theory was established for modeling the continuous fatigue process from the early dislocation stage to the crack expansion stage. Solutions obtained from the nonlinear wave equation with dislocations [24][25][26][27] were added into the nonlinear spring model. For the early stage, spring stiffness approaches infinity. Fatigue damage is mainly caused by dislocations in this stage. For the late stage damages, varying stiffness was chosen for modeling the different fatigue degree with the sprouting of micro-cracks. Once the crack begins to expand, spring stiffness declines quickly, and the amplitude of each order harmonic drops obviously. The nonlinear acoustic experimental system was built, and measurements were carried out during the whole fatigue process of the sample, which certified the numerical models.

Nonlinearity Due to Lattice Dislocations at Early Stage Fatigue
Early stage fatigue damage of metals mainly originates from dislocations. As the degree of fatigue further deepens, the dislocation density increases progressively, and the nonlinear acoustic effects becomes more obvious. The dislocation can be described by the dipole model [10], as shown in Figure 1.  A compact tension (CT) specimen (Figure 1a) was chosen for modeling. Stress during the fatigue loading mainly concentrates on the initial grooving. At the early stage of the loading process, the appearance of dislocations is considered as pairs of dipoles, as sketched in Figure 1b: (x, h) is the Substituting Equation (8) and Equation (6b) into the following wave motion equation: we have: with: where u is the total displacement, β l , β d are the nonlinear coefficients due to the material lattices and dislocations, respectively, c = A l 2 /ρ is the wave velocity, and ρ is the mass density of the material. Therefore, Equation (10a,b) is the nonlinear acoustic wave equation in materials with dislocations, from which the relationship between dislocation density and nonlinear harmonics can be established. The perturbation method is usually employed to solve this equation.
Assuming the following solution form that is expanded with the Mach number: where the Mach number M = V 0 /c and V 0 is the velocity amplitude of the vibration. Substituting Equation (11) into Equation (10a): with: Factoring the above equation in terms of M i , M ··· : · · · · · · (13c) from which the fundamental and higher harmonics of the propagating wave can be resolved by the perturbation method. The fundamental solution of Equation (13a) in the positive x direction has the form: where A 0 is the wave amplitude at the initial position (usually a driving source), k is the wave number, and ω is the angular frequency. For the second harmonic, substituting Equation (14) into Equation (13b), Materials 2019, 12, 607 5 of 15 whose solution can be found to be, In the same way, the former order solutions are used to determine the next order one: and so on, to obtain u 5 , u 6 , · · · · · · , u n . Substituting all the solutions into Equation (11) and combining the terms having the same harmonic frequency, we have, where A n = s n 2 + c n 2 and arctan c n s n (for n = 1, 2, 3, 4, · · ·) are the amplitude and phase of the nth harmonic wave component, respectively. A n , s n , c n are functions of propagation distance x and wave amplitude A 0 at the initial position. The first four orders are discussed in our work and are given as follows: Through the above formulations, the relationship between the dislocation density and high order harmonics can be deduced.

Nonlinearity Due to Micro-Cracks at Late Stage Fatigue
While the fatigue of the metals is gradually strengthened, a further increase of the dislocation density leads to the occurrence of micro-cracks, which means a weak interface is produced in the material. The expansion of the micro-cracks will change the nonlinear behavior of the transmission wave, in both its transmitted energy and harmonics distribution. The acoustic nonlinearity mechanism in this case differs from that of the dislocation model. It can be established with the following nonlinear spring model, as shown in Figure 2. density leads to the occurrence of micro-cracks, which means a weak interface is produced in the material. The expansion of the micro-cracks will change the nonlinear behavior of the transmission wave, in both its transmitted energy and harmonics distribution. The acoustic nonlinearity mechanism in this case differs from that of the dislocation model. It can be established with the following nonlinear spring model, as shown in Figure 2. A and B are two planes from either side of a crack interface. They are connected by the interfacial force, behaving like springs. Supposing that the displacement of an incident acoustic wave at plane A is − , and + , ℎ − represents that of the reflected and transmitted wave, respectively. P is used to simulate the static force of the closed crack, and G is the joint dynamic force caused by the motion of the interface. The midpoint of the interface is taken as the original point " ", where , represent the upper side A and the lower side B of the interface. , , , are the displacement of sides A and B, respectively. The relationship between the interface spacing = , − , and , is given as the following [16,17], Nonlinear springs A and B are two planes from either side of a crack interface. They are connected by the interfacial force, behaving like springs. Supposing that the displacement of an incident acoustic wave at plane A is f (x − ct), and g(x + ct), h(x − ct) represents that of the reflected and transmitted wave, respectively. P is used to simulate the static force of the closed crack, and G is the joint dynamic force caused by the motion of the interface. The midpoint of the interface is taken as the original point "x 0 ", where x 0− , x 0+ represent the upper side A and the lower side B of the interface. u(x 0− , t), u(x 0+ , t) are the displacement of sides A and B, respectively. The relationship between the interface spacing y(t) = u(x 0+ , t) − u(x 0− , t) and G(y), F(t) is given as the following [16,17], where F(t) is caused by the incident ultrasonic wave, From Equation (21), the variation of the crack pitch y(t) is related to both the joint dynamic force and ultrasonic force, as expressed in Equation (22). G(y) is usually nonlinear, whose relationship with interface spacing is given by the nonlinear spring model.
where K 1 , K 2 are the linear and nonlinear stiffness of the spring, respectively, which are related to the interface conditions. Considering the former dislocation model, the incident wave was already nonlinear because of the dislocation damage during the early stage before it propagated through the micro-cracks' area, so it can be obtained from the former results at position x 0 in the following form according to Equation (19): with: Substitute Equations (23) and (24a,b) into Equation (21), then, The solution of the above equation is composed of two parts: one linear part y 1 and a series of nonlinear parts y 2 , . . . y = y 1 + y 2 + · · · , y 1 y 2 · · · (26) The linear part y 1 is easy to obtain by solving Equation (28): . Then, where φ n = arctan nωρc 2K 1 , and y 2 can be resolved similarly from the following equation, while the small terms y 1 y 2 , y 2 2 , . . . are omitted. .
The harmonics of four orders are considered in the result, The transmitted wave through the spring interface can be obtained with the following equations [16,17], The transmitted coefficient of each order harmonic can be obtained with the above solution,

Computation Results
The transmission coefficients are related to the stiffness K 1 , as shown in Figure 3, the second-order terms are small enough to be ignored compared with the linear parts. The computation parameters of the aluminum alloy are chosen as an example, and the values of relevant parameters are referenced from the work of John H. Cantrell [9], which are given in Table 1. Table 1. Computational parameters of the aluminum alloy.

Parameters
Value Unit  It is obvious that the transmitted coefficient of each order harmonic approaches one, while → ∞, usually ≥ 1 × 10 / , which can be considered as full transmission, as shown in Region I: the early dislocation fatigue stage in Figure 3. Therefore, during this stage, acoustic nonlinearity mainly results from dislocations and wave energy transmitted completely without the effects of the weak interface, as described in the dislocation dipole model. Once the micro-cracks sprout, from = 1 × 10 / , spring stiffness begins to decline, and the energy of the transmitted wave drops gradually, as shown in Region II: micro-crack sprouting stage. A large amount of micro-crack sprouting leads to the rapid expansion of the crack length, which reduces the transmitted wave energy clearly; it is the last stage of the fatigue process in Region III: crack expansion stage. During the micro-crack sprouting stage, nonlinearity increased, owing to the increasing dislocation density being very weak compared to the effects of spring the stiffness decline with continuous fatigue loading. Relative results are given in Figure 4. Computation parameters are also given in Table 1. The initial amplitude of each order harmonic is chosen from the above dislocation model results at = 50mm, driving source amplitude = 1 . Assuming that micro-cracks sprout while the dislocation density is about 3 × 10 , which is the early fatigue stage, the harmonics variations are given in Figure 4a,c. After this stage, the spring stiffness in the computation model was taken from ∞ → 0, which means 1/ from 0 → ∞, as shown in Figure 4b,d; the dislocation density was assumed as constant during this process. It is obvious that the transmitted coefficient of each order harmonic approaches one, while K 1 → ∞ , usually K 1 ≥ 1 × 10 15 N/m 3 , which can be considered as full transmission, as shown in Region I: the early dislocation fatigue stage in Figure 3. Therefore, during this stage, acoustic nonlinearity mainly results from dislocations and wave energy transmitted completely without the effects of the weak interface, as described in the dislocation dipole model. Once the micro-cracks sprout, from K 1 = 1 × 10 15 N/m 3 , spring stiffness begins to decline, and the energy of the transmitted wave drops gradually, as shown in Region II: micro-crack sprouting stage. A large amount of micro-crack sprouting leads to the rapid expansion of the crack length, which reduces the transmitted wave energy clearly; it is the last stage of the fatigue process in Region III: crack expansion stage. During the micro-crack sprouting stage, nonlinearity increased, owing to the increasing dislocation density being very weak compared to the effects of spring the stiffness decline with continuous fatigue loading. Relative results are given in Figure 4. Computation parameters are also given in Table 1. The initial amplitude of each order harmonic is chosen from the above dislocation model results at x 0 = 50mm, driving source amplitude A 0 = 1µm. Assuming that micro-cracks sprout while the dislocation density is about 3 × 10 15 m −3 , which is the early fatigue stage, the harmonics variations are given in Figure 4a,c. After this stage, the spring stiffness in the computation model was taken from ∞ → 0 , which means 1/K 1 from 0 → ∞ , as shown in Figure 4b,d; the dislocation density was assumed as constant during this process. initial amplitude of each order harmonic is chosen from the above dislocation model results at = 50mm, driving source amplitude = 1 . Assuming that micro-cracks sprout while the dislocation density is about 3 × 10 , which is the early fatigue stage, the harmonics variations are given in Figure 4a,c. After this stage, the spring stiffness in the computation model was taken from ∞ → 0, which means 1/ from 0 → ∞, as shown in Figure 4b,d; the dislocation density was assumed as constant during this process.  The above results show that acoustic nonlinearity is very weak for materials without damages, as given in Figure 4c, while the dislocation density is zero. However, with the increase of dislocation density, which means the deepening of fatigue damage, acoustic nonlinearity becomes much heavier. The fundamental wave amplitude is going to decrease, while harmonic amplitudes begin to rise due to the energy distribution law that energy transfers from the lower order signals to the higher orders, as shown in Region I in Figure 4a,c. Then, acoustic nonlinearity indicates the fatigue degree in materials during the early fatigue stage. Micro-cracks sprout with further fatigue loading, and the value of spring stiffness begins to drop, whose effects on each order harmonic are much heavier than that of the increasing dislocation density, so the harmonic amplitudes are going to decrease, as shown in Region II in Figure 4b,d. The energy of each order transmitted wave declines very quickly with the decreasing stiffness, owing to the rapid expansion of cracks in Region III.

Experiment Results and Discussions
Based on the above analysis for the theoretical models during the fatigue process, the following experiments were designed and carried out. The experimental setup employed is represented in Figure 5a. It was composed of: a signal generator (Tektronix AFG1022, Tektronix, Shanghai, China), a power amplifier (Model 2100L RF Power Amplifier, ENI, Rochester, NY, USA), two transducers (central frequency of 1.1 MHz for the emitting one and 2.3 MHz for the receiving one, respectively), and a digital oscilloscope (Agilent Technologies MSO7032A, Agilent Technologies Inc., Santa Clara, CA, USA).
The receiving transducer sensitivity at 1.1 MHz was about −11.4 dB, for a signal of 2.2 MHz (the second-order signal), 3.3 MHz (the third-order one), and 4.4 MHz (the fourth-order one), and the sensitivities were −2 dB, −6 dB, and −12 dB, respectively. The gains of the receiving transducers were compensated in the following high-order harmonic amplitude measurements.  The above results show that acoustic nonlinearity is very weak for materials without damages, as given in Figure 4c, while the dislocation density is zero. However, with the increase of dislocation density, which means the deepening of fatigue damage, acoustic nonlinearity becomes much heavier. The fundamental wave amplitude is going to decrease, while harmonic amplitudes begin to rise due to the energy distribution law that energy transfers from the lower order signals to the higher orders, as shown in Region I in Figure 4a,c. Then, acoustic nonlinearity indicates the fatigue degree in materials during the early fatigue stage. Micro-cracks sprout with further fatigue loading, and the value of spring stiffness begins to drop, whose effects on each order harmonic are much heavier than that of the increasing dislocation density, so the harmonic amplitudes are going to decrease, as shown in Region II in Figure 4b,d. The energy of each order transmitted wave declines very quickly with the decreasing stiffness, owing to the rapid expansion of cracks in Region III.

Experiment Results and Discussion
Based on the above analysis for the theoretical models during the fatigue process, the following experiments were designed and carried out. The experimental setup employed is represented in Figure 5a. It was composed of: a signal generator (Tektronix AFG1022, Tektronix, Shanghai, China), a power amplifier (Model 2100L RF Power Amplifier, ENI, Rochester, NY, USA), two transducers (central frequency of 1.1 MHz for the emitting one and 2.3 MHz for the receiving one, respectively), and a digital oscilloscope (Agilent Technologies MSO7032A, Agilent Technologies Inc., Santa Clara, CA, USA).
The receiving transducer sensitivity at 1.1 MHz was about −11.4 dB, for a signal of 2.2 MHz (the second-order signal), 3.3 MHz (the third-order one), and 4.4 MHz (the fourth-order one), and the sensitivities were −2 dB, −6 dB, and −12 dB, respectively. The gains of the receiving transducers were compensated in the following high-order harmonic amplitude measurements.
Based on the above analysis for the theoretical models during the fatigue process, the following experiments were designed and carried out. The experimental setup employed is represented in Figure 5a. It was composed of: a signal generator (Tektronix AFG1022, Tektronix, Shanghai, China), a power amplifier (Model 2100L RF Power Amplifier, ENI, Rochester, NY, USA), two transducers (central frequency of 1.1 MHz for the emitting one and 2.3 MHz for the receiving one, respectively), and a digital oscilloscope (Agilent Technologies MSO7032A, Agilent Technologies Inc., Santa Clara, CA, USA).
The receiving transducer sensitivity at 1.1 MHz was about −11.4 dB, for a signal of 2.2 MHz (the second-order signal), 3.3 MHz (the third-order one), and 4.4 MHz (the fourth-order one), and the sensitivities were −2 dB, −6 dB, and −12 dB, respectively. The gains of the receiving transducers were compensated in the following high-order harmonic amplitude measurements.  Before the preparation of the fatigue specimen, the fatigue model of the aluminum alloy specimen was established using finite element software (COMSOL Multiphysics 5.3, COMSOL Inc., Stockholm, Sweden). As in the fatigue test, the same size model was established with the same material. The lower surface boundary was fixed, and the upper surface boundary was loaded with the maximum stress during the fatigue test. The total stress distribution of this specimen in the fatigue process was analyzed, and the results are given in Figure 5b. It was shown that the stress in the test piece was concentrated along the initial grooving. If the fatigue degree increased constantly, tiny cracks would start to sprout and gradually extend forward. According to the calculation results of the model, the fatigue test was carried out on the sample by selecting appropriate parameters, the process of which is shown in Figure 6a. Before the preparation of the fatigue specimen, the fatigue model of the aluminum alloy specimen was established using finite element software (COMSOL Multiphysics 5.3, COMSOL Inc., Stockholm, Sweden). As in the fatigue test, the same size model was established with the same material. The lower surface boundary was fixed, and the upper surface boundary was loaded with the maximum stress during the fatigue test. The total stress distribution of this specimen in the fatigue process was analyzed, and the results are given in Figure 5b. It was shown that the stress in the test piece was concentrated along the initial grooving. If the fatigue degree increased constantly, tiny cracks would start to sprout and gradually extend forward. According to the calculation results of the model, the fatigue test was carried out on the sample by selecting appropriate parameters, the process of which is shown in Figure 6a. The fatigue test was carried out at room temperature using the stress-loaded control mode, during which the average fatigue loading force was 5 kN and the loading frequency was 5 Hz. After a certain period of loading, nonlinear acoustic measurements were carried out and then redetected after another period of loading, until micro-cracks sprouted in the sample. The loading cycle was 500, 1500, 3000, 5000, 10,000, 15,000, 20,000, 25,000, and 29,800, respectively. Before the fatigue loading, acoustic measurement was performed to get the result of the non-damaged sample as a comparison, which is given in Figure 7. While the loading cycle was about 5000 during the fatigue test, acoustic nonlinearity increased very clearly in the received signals, as shown in Figure 8. Micro-cracks were observed from the end of the preset grooving after 20,000 cycles of loading, the results of which are given in Figure 9. Once micro-cracks sprouted, continuous loading accelerated their expansion. Fatigue loading was performed until the crack could be seen clearly in the sample, as shown in Figure  6b, whose loading cycle was about 29,800, and the measurement results are given in Figure 10. The fatigue test was carried out at room temperature using the stress-loaded control mode, during which the average fatigue loading force was 5 kN and the loading frequency was 5 Hz. After a certain period of loading, nonlinear acoustic measurements were carried out and then redetected after another period of loading, until micro-cracks sprouted in the sample. The loading cycle was 500, 1500, 3000, 5000, 10,000, 15,000, 20,000, 25,000, and 29,800, respectively. Before the fatigue loading, acoustic measurement was performed to get the result of the non-damaged sample as a comparison, which is given in Figure 7. While the loading cycle was about 5000 during the fatigue test, acoustic nonlinearity increased very clearly in the received signals, as shown in Figure 8. Micro-cracks were observed from the end of the preset grooving after 20,000 cycles of loading, the results of which are given in Figure 9. Once micro-cracks sprouted, continuous loading accelerated their expansion. Fatigue loading was performed until the crack could be seen clearly in the sample, as shown in Figure 6b, whose loading cycle was about 29,800, and the measurement results are given in Figure 10.
during which the average fatigue loading force was 5 kN and the loading frequency was 5 Hz. After a certain period of loading, nonlinear acoustic measurements were carried out and then redetected after another period of loading, until micro-cracks sprouted in the sample. The loading cycle was 500, 1500, 3000, 5000, 10,000, 15,000, 20,000, 25,000, and 29,800, respectively. Before the fatigue loading, acoustic measurement was performed to get the result of the non-damaged sample as a comparison, which is given in Figure 7. While the loading cycle was about 5000 during the fatigue test, acoustic nonlinearity increased very clearly in the received signals, as shown in Figure 8. Micro-cracks were observed from the end of the preset grooving after 20,000 cycles of loading, the results of which are given in Figure 9. Once micro-cracks sprouted, continuous loading accelerated their expansion. Fatigue loading was performed until the crack could be seen clearly in the sample, as shown in Figure  6b, whose loading cycle was about 29,800, and the measurement results are given in Figure 10. Comparing the results in Figures 7a and 8a, the signal received from the sample after the fatigue test (Loading Cycle 5000) deformed much more heavily than that of the non-fatigue situation (Loading Cycle 0), which is much more obvious in the spectrum results in Figures 7b and 8b. After fatigue loading, which means the increase of dislocation density, as explained in the dislocation dipole model, harmonics caused by the damage increased also; therefore, the spectrum was much richer in harmonic components, and each-order harmonic could be seen more clearly than the result of the non-damage situation. Comparing the results in Figures 7a and 8a, the signal received from the sample after the fatigue test (Loading Cycle 5000) deformed much more heavily than that of the non-fatigue situation (Loading Cycle 0), which is much more obvious in the spectrum results in Figures 7b and 8b. After fatigue loading, which means the increase of dislocation density, as explained in the dislocation dipole model, harmonics caused by the damage increased also; therefore, the spectrum was much richer in harmonic components, and each-order harmonic could be seen more clearly than the result of the non-damage situation.
The difference is not obvious in the time-domain result in Figure 9a, except that the transmitted energy is a little lower compared with the result in Figure 8a. However, in the spectrum results in Figure 9b, not only the harmonic amplitudes are a little lower, but another frequency component whose frequency is between the fundamental wave and second-order harmonics appears, compared with the results in Figure 8b, which may be caused by the wave scattering of the crack tip or the spatial variation in the stiffness of the springs. These results are much more apparent in Figure 10. As micro-cracks expanded very quickly with continuous loading, the energy of the transmitted wave in this stage was much lower than that in the former situations, as given in the time-domain result, which was consistent with the spring model result. In addition, more spectrum components appeared in the spectrum result, as marked with rectangles in Figure 10b. Because the crack surface was rough and curved, this may have led to more scattering waves and generated more other intermediate frequency components in the result. It is very complicated to describe in theoretical models and will be studied in our next work. Comparing the results in Figures 7a and 8a, the signal received from the sample after the fatigue test (Loading Cycle 5000) deformed much more heavily than that of the non-fatigue situation (Loading Cycle 0), which is much more obvious in the spectrum results in Figures 7b and 8b. After fatigue loading, which means the increase of dislocation density, as explained in the dislocation dipole model, harmonics caused by the damage increased also; therefore, the spectrum was much richer in harmonic components, and each-order harmonic could be seen more clearly than the result of the non-damage situation. The difference is not obvious in the time-domain result in Figure 9a, except that the transmitted energy is a little lower compared with the result in Figure 8a. However, in the spectrum results in Figure 9b, not only the harmonic amplitudes are a little lower, but another frequency component whose frequency is between the fundamental wave and second-order harmonics appears, compared with the results in Figure 8b, which may be caused by the wave scattering of the crack tip or the spatial variation in the stiffness of the springs. These results are much more apparent in Figure 10. As micro-cracks expanded very quickly with continuous loading, the energy of the transmitted wave in this stage was much lower than that in the former situations, as given in the time-domain result, which was consistent with the spring model result. In addition, more spectrum components appeared in the spectrum result, as marked with rectangles in Figure 10b. Because the crack surface was rough and curved, this may have led to more scattering waves and generated more other intermediate frequency components in the result. It is very complicated to describe in theoretical models and will be studied in our next work. The above results are some typical measurements for explaining the acoustic nonlinearity variation during different fatigue stages. Detailed amplitude variation of each order harmonic with fatigue degree, which is described by the fatigue loading cycle, is shown in Figure 11. During each measurement, excitation voltages were taken as 150 V (Figure 11a The above results are some typical measurements for explaining the acoustic nonlinearity variation during different fatigue stages. Detailed amplitude variation of each order harmonic with fatigue degree, which is described by the fatigue loading cycle, is shown in Figure 11. During each measurement, excitation voltages were taken as 150 V (Figure 11a,b), 175 V (Figure 11c,d), and 200 V (Figure 11e,f), respectively.
The above results are some typical measurements for explaining the acoustic nonlinearity variation during different fatigue stages. Detailed amplitude variation of each order harmonic with fatigue degree, which is described by the fatigue loading cycle, is shown in Figure 11. During each measurement, excitation voltages were taken as 150 V (Figure 11a,b), 175 V (Figure 11c Generally, the experiment results can also be divided into three stages: (I) the early dislocation fatigue stage, (II) the micro-crack sprouting stage, and (III) the crack expansion stage, which is consistent with the theoretical model results, as shown in the above figures. While the loading cycle was about less than 10,000, harmonic components increased obviously with continuous loading, which can be treated as the early stage of fatigue. In this stage, acoustic nonlinearity mainly came from dislocations, and it can be used to assess the early stage fatigue degree, as explained in the theoretical model of Region I. Further loading led to micro-cracks sprouting inside of the metals, and the occurrence of a weak interface reduced the transmitted wave energy. Therefore, the increasing trend of the harmonics amplitude became slow even with higher dislocation density, which was the transition stage. The boundary line of this stage was between 10,000 and 15,000 cycles, approximately, which corresponds to the middle part in Figure 4. While the loading cycle was about 20,000, cracks were seen clearly at the end of the preset grooving, which meant the beginning of the crack expansion stage. The amplitude of each order harmonic declined quickly during this stage, the results of which Generally, the experiment results can also be divided into three stages: (I) the early dislocation fatigue stage, (II) the micro-crack sprouting stage, and (III) the crack expansion stage, which is consistent with the theoretical model results, as shown in the above figures. While the loading cycle was about less than 10,000, harmonic components increased obviously with continuous loading, which can be treated as the early stage of fatigue. In this stage, acoustic nonlinearity mainly came from dislocations, and it can be used to assess the early stage fatigue degree, as explained in the theoretical model of Region I. Further loading led to micro-cracks sprouting inside of the metals, and the occurrence of a weak interface reduced the transmitted wave energy. Therefore, the increasing trend of the harmonics amplitude became slow even with higher dislocation density, which was the transition stage. The boundary line of this stage was between 10,000 and 15,000 cycles, approximately, which corresponds to the middle part in Figure 4. While the loading cycle was about 20,000, cracks were seen clearly at the end of the preset grooving, which meant the beginning of the crack expansion stage. The amplitude of each order harmonic declined quickly during this stage, the results of which are consistent with the theoretical results of Region III in Figure 4.

Conclusions
Acoustic nonlinearity during the fatigue process in metals can be divided into three stages: (I) the early dislocation fatigue stage, (II) the micro-crack sprouting stage, and (III) the crack expansion stage, which can be described by nonlinear spring models combined with dislocation theory. Varying spring stiffness is chosen to connect different stages in the fatigue process. For the early dislocation fatigue stage, spring stiffness approaches infinity, and acoustic nonlinearity mainly results from dislocations, the results of which can be used for early fatigue evaluation. Once micro-cracks sprout in Region II, spring stiffness will drop, which reduces the transmitted wave energy and makes the harmonics fall. The decreasing trends become steep quickly, owing to the crack's expansion, and additional spectrum components are generated in the results with crack occurrence in Region III. These may be the effects of the scattering wave on rough crack interface or the spatial variation in the stiffness of the springs. This phenomenon is complicated, and studying its nonlinear acoustic mechanism is our next work.