Next Article in Journal
Thermal, Mechanical, and Microstructural Study of PBO Fiber during Carbonization
Next Article in Special Issue
Magnetostrictive Guided Wave Technique Verification for Detection and Monitoring Defects in the Pipe Weld
Previous Article in Journal
Improvement of Asphalt-Aggregate Adhesion Using Plant Ash Byproduct
Previous Article in Special Issue
Computation of Propagating and Non-Propagating Lamb-Like Wave in a Functionally Graded Piezoelectric Spherical Curved Plate by an Orthogonal Function Technique
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Nonlinear Acoustic Modeling and Measurements during the Fatigue Process in Metals

1
State Key Laboratory of Acoustics, Institute of Acoustics, Chinese Academy of Sciences, Beijing 100190, China
2
University of Chinese Academy of Sciences, Beijing 100049, China
3
Université Polytechnique Hauts-de-France, CNRS, Univ. Lille, YNCREA, Centrale Lille, UMR 8520IEMN-DOAE, F-59313 Valenciennes CEDEX 9, France
*
Author to whom correspondence should be addressed.
Materials 2019, 12(4), 607; https://doi.org/10.3390/ma12040607
Submission received: 11 January 2019 / Revised: 9 February 2019 / Accepted: 15 February 2019 / Published: 18 February 2019
(This article belongs to the Special Issue Non-destructive Inspection)

Abstract

:
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.

1. 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.

2. Theoretical Models

2.1. 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 position of one dislocation P in the pair relative to another one O at ( 0 , 0 ) . The two dislocation lines are parallel and lie perpendicular to the “xy” plane. The unit length glide force [9,10] that acts between the two dislocations in the x direction is given as follows:
F x = G b 2 2 π ( 1 v ) x ( x 2 h 2 ) ( x 2 + h 2 ) 2
where G is the shear modulus, v is Poisson’s ratio, b is the magnitude of Burger’s vector, and h represents the spacing between the two glide planes. Propagation of an ultrasonic wave in the material produces the strains of two distinguished origins: the elastic strain ε l due to the lattice structure of the material and the strain ε d associated with the motion of dislocation dipoles.
ε = ε l + ε d
The relationship between the elastic strain ε l and stress σ is given by [9,28]:
ε l = 1 A 2 l σ 1 2 A 3 l ( A 2 l ) 3 σ 2 +
where A 2 l ,   A 3 l are the second- and third-order elastic constants, respectively, while the contribution of strain ε d to the stress σ can be obtained from the dislocation motion equation [10]. Only the relative motion of the dipoles along the direction parallel to the glide planes is considered:
m e f f 2 x t 2 = F x + b σ R
where m e f f is the effective mass of the dislocation dipoles, and R is the Schmid factor. Let ζ = x h , which is the relative displacement of dislocation dipoles to the position h. Combining Equations (1) and (4) and expanding Equation (1) in the series method with ζ , the following is obtained:
m e f f 2 ζ t 2 + G b 2 4 π h 2 ( 1 v ) ζ G b 2 8 π h 3 ( 1 v ) ζ 2 + = b σ R
Solving this equation and using the relationship ε d = Ω Λ d b ζ , where, Ω is the conversion factor and Λ d is the dislocation density, the following relationship is obtained:
ε d = 1 A 2 d σ 1 2 A 3 d ( A 2 d ) 3 σ 2 +
with:
A 2 d = G 4 π Ω R Λ d h 2 ( 1 v ) ,     A 3 d = G 4 π Ω 2 R Λ d 2 h 3 ( 1 v ) b
Substituting Equations (3) and (6a,b) into Equation (2) yields:
ε = ( 1 A 2 l + 1 A 2 d ) σ 1 2 [ A 3 d ( A 2 d ) 3 + A 3 l ( A 2 l ) 3 ] σ 2 +
The inverse expression can be written in the following form approximately:
σ = A 2 l ε 1 2 [ A 3 l + A 3 d ( A 2 l ) 3 ( A 2 d ) 3 ] ε 2 +
Substituting Equation (8) and Equation (6b) into the following wave motion equation:
ρ 2 u t 2 = σ x ,           ε = u x
we have:
2 u t 2 c 2 2 u x 2 = c 2 ( β l + β d ) 2 u x 2 u x +
with:
β l = A 3 l A 2 l ,               β d = A 3 d ( A 2 l ) 2 ( A 2 d ) 3 = 16 π 2 Ω R 2 Λ d h 3 ( 1 v ) 2 ( A 2 l ) 2 G 2 b
where u is the total displacement, β l ,     β d are the nonlinear coefficients due to the material lattices and dislocations, respectively, c = A 2 l / ρ 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:
u = i = 1 n M i u i .
where the Mach number M = V 0 / c and V 0 is the velocity amplitude of the vibration. Substituting Equation (11) into Equation (10a):
1 c 2 i = 1 n 2 t 2 ( M i u i ) = i = 1 n 2 x 2 ( M i u i ) + β i = 1 n x ( M i u i ) 2 x 2 ( M i u i ) +
with:
β = β l + β d
Factoring the above equation in terms of M i ,
M 1 :     1 c 2 2 u 1 t 2 = 2 u 1 x 2 ,
M 2 :   1 c 2 2 u 2 t 2 = 2 u 2 x 2 + β u 1 x 2 u 1 x 2 ,
M :                                                                                              
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:
u 1 = A 0 cos ( k x ω t ) ,
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),
1 c 2 2 u 2 t 2 2 u 2 x 2 = 1 2 β A 0 2 k 3 sin ( 2 k x 2 ω t )
whose solution can be found to be,
u 2 = 1 8 β A 0 2 k 2 x cos ( 2 k x 2 ω t )
In the same way, the former order solutions are used to determine the next order one:
u 3 = β 2 [ 1 32 A 0 3 k 4 x 2 cos ( k x ω t ) 3 32 A 0 3 k 3 x sin ( k x ω t ) + 1 32 A 0 3 k 4 x 2 cos ( 3 k x 3 ω t ) + 5 96 A 0 3 k 3 x sin ( 3 k x 3 ω t ) ]
u 4 = β 3 [ ( 1 96 A 0 4 k 6 x 3 + 1 12 A 0 4 k 4 x ) cos ( 2 k x 2 ω t ) 33 512 A 0 4 k 5 x 2 sin ( 2 k x 2 ω t ) + ( 1 96 A 0 4 k 6 x 3 11 512 A 0 4 k 4 x ) cos ( 4 k x 4 ω t ) + 5 128 A 0 4 k 5 x 2 sin ( 4 k x 4 ω t ) ]
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,
u = n A n sin ( n k x n ω t + arctan c n s n )
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 n th 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:
s 1 = 3 32 M 3 β 2 A 0 3 k 3 x ,         c 1 = A 0 M 1 32 M 3 β 2 A 0 3 k 4 x 2
s 2 = M 4 β 3 33 512 A 0 4 k 5 x 2 ,         c 2 = 1 8 M 2 β A 0 2 k 2 x + M 4 β 3 ( 1 96 A 0 4 k 6 x 3 + 1 12 A 0 4 k 4 x )
s 3 = 5 96 M 3 β 2 A 0 3 k 3 x ,           c 3 = 1 32 M 3 β 2 A 0 3 k 4 x 2
s 4 = 5 128 M 4 β 3 A 0 4 k 5 x 2 ,             c 4 = M 4 β 3 ( 1 96 A 0 4 k 6 x 3 11 512 A 0 4 k 4 x )
M = V 0 c = ω A 0 c
Through the above formulations, the relationship between the dislocation density and high order harmonics can be deduced.

2.2. 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.
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 c t ) , and g ( x + c t ) , h ( x c t ) 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],
ρ c y ˙ ( t ) = 2 G ( y ) 2 F ( t )
where F ( t ) is caused by the incident ultrasonic wave,
F ( t ) = ( ρ c 2 f ( x c t ) x ) x = x 0
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.
G ( y ) = K 1 y + K 2 y 2 +
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):
f ( x c t ) = n A n sin ( n k x 0 n ω t + φ n )
with:
φ n = arctan c n s n , n = 1 , 2 , 3
Substitute Equations (23) and (24a,b) into Equation (21), then,
y ˙ + 2 K 1 ρ c y 2 K 2 ρ c y 2 = 2 n n ω A n cos ( n k x 0 n ω t + φ n )
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
Then,
y 1 ˙ + y 2 ˙ + 2 K 1 ρ c ( y 1 + y 2 + ) 2 K 2 ρ c ( y 1 + y 2 + ) 2 = 2 n n ω A n cos ( n k x 0 n ω t + φ n )
The linear part y 1 is easy to obtain by solving Equation (28):
y 1 ˙ + 2 K 1 ρ c y 1 = 2 n n ω A n cos ( n k x 0 n ω t + φ n )
Then,
y 1 = 2 n n A n 4 K 1 2 / n 2 ω 2 ρ 2 c 2 + n 2 cos ( n k x 0 n ω t + φ n + ϕ n )
where ϕ n = arctan n ω ρ c 2 K 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.
y 2 ˙ + 2 K 1 ρ c y 2 = 2 K 2 ρ c y 1 2
The harmonics of four orders are considered in the result,
y = 2 A 1 cos ( k x 0 ω t + φ 1 + ϕ 1 ) 4 K 1 2 / ω 2 ρ 2 c 2 + 1 + 4 A 2 cos ( 2 k x 0 2 ω t + φ 2 + ϕ 2 ) 4 K 1 2 / ω 2 ρ 2 c 2 + 4 + 6 A 3 cos ( 3 k x 0 3 ω t + φ 3 + ϕ 3 ) 4 K 1 2 / 9 ω 2 ρ 2 c 2 + 9 + 8 A 4 cos ( 4 k x 0 4 ω t + φ 4 + ϕ 4 ) 4 K 1 2 / 16 ω 2 ρ 2 c 2 + 16 2 A 1 2 K 2 sin ( 2 k x 0 2 ω t + 2 φ 1 + 2 ϕ 1 + ϕ 2 ) ( 4 K 1 2 / ω 2 ρ 2 c 2 + 1 ) 4 K 1 2 + 4 ω 2 ρ 2 c 2 8 K 2 A 2 2 sin ( 4 k x 0 4 ω t + φ 4 + 2 ϕ 2 + ϕ 4 ) ( 4 K 1 2 / ω 2 ρ 2 c 2 + 4 ) 4 K 1 2 + 16 ω 2 ρ 2 c 2 + 4 K 2 A 2 A 1 ( 4 K 1 2 / ω 2 ρ 2 c 2 + 1 ) ( 4 K 1 2 / ω 2 ρ 2 c 2 + 4 ) ( cos ( k x 0 ω t + φ 3 φ 1 + ϕ 2 ) 4 K 1 2 + ω 2 ρ 2 c 2 + cos ( 3 k x 0 3 ω t + φ 3 + φ 1 + ϕ 1 + ϕ 2 + ϕ 3 ) 4 K 1 2 + 9 ω 2 ρ 2 c 2 ) +
The transmitted wave through the spring interface can be obtained with the following equations [16,17],
h ( x c t ) = f ( x c t ) + 1 2 y ( t x c )
Then,
h ( x c t ) = A 1 sin ( k x 0 ω t + k x + φ 1 ) + A 2 sin ( 2 k x 0 2 ω t + 2 k x + φ 2 ) + A 3 sin ( 3 k x 0 3 ω t + 3 k x + φ 3 ) + A 4 sin ( 4 k x 0 4 ω t + 4 k x + φ 4 ) + A 1 cos ( k x 0 ω t + k x + φ 1 + ϕ 1 ) 4 K 1 2 / ω 2 ρ 2 c 2 + 1 + 2 A 2 cos ( 2 k x 0 2 ω t + 2 k x + φ 2 + ϕ 2 ) 4 K 1 2 / 4 ω 2 ρ 2 c 2 + 4 + 3 A 3 cos ( 3 k x 0 3 ω t + 3 k x + φ 3 + ϕ 3 ) 4 K 1 2 / 9 ω 2 ρ 2 c 2 + 9 + 4 A 4 cos ( 4 k x 0 4 ω t + 4 k x + φ 4 + ϕ 4 ) 4 K 1 2 / 16 ω 2 ρ 2 c 2 + 16 2 A 1 2 K 2 sin ( 2 k x 0 2 ω t + 2 k x + 2 φ 1 + 2 ϕ 1 + ϕ 2 ) ( 4 K 1 2 / ω 2 ρ 2 c 2 + 1 ) 4 K 1 2 + 4 ω 2 ρ 2 c 2 4 K 2 A 2 2 sin ( 4 k x 0 4 ω t + 4 k x + 2 φ 2 + 2 ϕ 2 + ϕ 4 ) ( 4 K 1 2 / ω 2 ρ 2 c 2 + 4 ) 4 K 1 2 + 16 ω 2 ρ 2 c 2 + 2 K 2 A 2 A 1 ( 4 K 1 2 / ω 2 ρ 2 c 2 + 1 ) ( 4 K 1 2 / ω 2 ρ 2 c 2 + 4 ) ( cos ( k x 0 ω t + k x + φ 3 φ 1 + ϕ 2 ) 4 K 1 2 + ω 2 ρ 2 c 2 + cos ( 3 k x 0 3 ω t + 3 k x + φ 3 + φ 1 + ϕ 1 + ϕ 2 + ϕ 3 ) 4 K 1 2 + 9 ω 2 ρ 2 c 2 ) +
The transmitted coefficient of each order harmonic can be obtained with the above solution,
T n = 4 K 1 2 / n 2 ω 2 ρ 2 c 2 + 2 n 2 + 2 n sin ϕ n 4 K 1 2 / n 2 ω 2 ρ 2 c 2 + n 2 4 K 1 2 / n 2 ω 2 ρ 2 c 2 + n 2 ,   n = 1 , 2 , 3 , 4

2.3. 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.
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 = 50 mm ,     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.
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.

3. 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.
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.
Comparing the results in Figure 7a and Figure 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 Figure 7b and Figure 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,b), 175 V (Figure 11c,d), and 200 V (Figure 11e,f), respectively.
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.

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.

Author Contributions

Conceptualization, W.L.; validation, X.W. and W.X.; experiments’ design and sample preparation, W.L. and X.W.; data curation, W.L. and W.X.; original draft preparation, W.L.; review and editing, X.W. and W.X.; supervision, X.W. and W.X.; funding acquisition, X.W.

Funding

This research was funded by the National Natural Science Foundation of China (11674351) and China Scholarship Council (201704910691).

Acknowledgments

The authors would like to express their gratitude for the help with the fatigue tests by the Technical Institute of Physics and Chemistry, Chinese Academy of Sciences.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Klepka, A.; Staszewski, W.; Jenal, R.; Szwedo, M.; Iwaniec, J.; Uhl, T. Nonlinear acoustics for fatigue crack detection—Experimental investigations of vibro-acoustic wave modulations. Struct. Health Monit. 2012, 11, 197–211. [Google Scholar] [CrossRef]
  2. Tsyfansky, S.; Beresnevich, V. Non-linear vibration method for detection of fatigue cracks in aircraft wings. J. Sound Vib. 2000, 236, 49–60. [Google Scholar] [CrossRef]
  3. Palitsagar, S.; Das, S.; Parida, N.; Bhattacharya, D. Non-linear ultrasonic technique to assess fatigue damage in structural steel. Scr. Mater. 2006, 55, 199–202. [Google Scholar] [CrossRef]
  4. Nagy, P.B. Fatigue damage assessment by nonlinear ultrasonic materials characterization. Ultrasonics 1998, 36, 375–381. [Google Scholar] [CrossRef]
  5. Yang, Z.; Tian, Y.; Li, W.; Zhou, H.; Zhang, W. Experimental Investigation of the Acoustic Nonlinear Behavior in Granular Polymer Bonded Explosives with Progressive Fatigue Damage. Materials 2017, 10, 660. [Google Scholar] [CrossRef] [PubMed]
  6. Landau, L.D.; Lifshitz, E.M. Theory of Elasticity, 3rd ed.; Pergamon Press: New York, NY, USA, 1959; pp. 104–107. [Google Scholar]
  7. Murnaghan, F.D. Finite Deformations of an Elastic Solid; John Wiley: New York, NY, USA, 1951; pp. 60–66. [Google Scholar]
  8. Granato, A.V.; Lücke, K. The Vibrating String Model of Dislocation Damping. In Physical Acoustics; Mason, W.P., Ed.; Academic Press: New York, NY, USA, 1966; Volume 3, pp. 225–276. [Google Scholar]
  9. Cantrell, J.H.; Yost, W.T. Nonlinear ultrasonic characterization of fatigue microstructures. Int. J. Fatigue 2001, 23, 487–490. [Google Scholar] [CrossRef]
  10. Cantrell, J.H.; Yost, W.T. Acoustic harmonic generation from fatigue-induced dislocation dipoles. Philos. Mag. A 1994, 69, 315–326. [Google Scholar] [CrossRef]
  11. Mayergoyz, I.D. Hysteresis models from the mathematical and control theory points of view. J. Appl. Phys. 1985, 57, 3803–3805. [Google Scholar] [CrossRef]
  12. Gonzalez-Valadez, M.; Baltazar, A.; Dwyer-Joyce, R. Study of interfacial stiffness ratio of a rough surface in contact using a spring model. Wear 2010, 268, 373–379. [Google Scholar] [CrossRef] [Green Version]
  13. Wu, X.; Li, H.; Lv, W.; Zhang, S.; Zhang, D. Numerical analysis on acoustic nonlinear phenomena between rough contact interfaces. In Proceedings of the 2016 IEEE International Ultrasonics Symposium (IUS), Tours, France, 18–21 September 2016; pp. 1–4. [Google Scholar]
  14. An, Z.; Wang, X.; Deng, M.; Mao, J.; Li, M. A nonlinear spring model for an interface between two solids. Wave Motion 2013, 50, 295–309. [Google Scholar] [CrossRef]
  15. Takahashi, S.; Motegi, R. Measurement of third-order elastic constants and applications to loaded structural materials. SpringerPlus 2015, 4, 325. [Google Scholar] [CrossRef] [PubMed]
  16. Richardson, J.M. Harmonic generation at an unbonded interface—I. Planar interface between semi-infinite elastic media. Int. J. Eng. Sci. 1979, 17, 73–85. [Google Scholar] [CrossRef]
  17. Buck, O. Acoustic harmonic generation at unbonded interfaces and fatigue cracks. J. Acoust. Soc. Am. 1978, 33, 371–373. [Google Scholar]
  18. Cantrell, J.H.; Pan, E. Substructural organization, dislocation plasticity and harmonic generation in cyclically stressed wavy slip metals. Proc. R. Soc. Lond. A Math. Phys. Eng. Sci. 2004, 460, 757–780. [Google Scholar] [CrossRef]
  19. Cantrell, J.H.; Yost, W.T. Effect of precipitate coherency strains on acoustic harmonic generation. J. Appl. Phys. 1997, 81, 2957–2962. [Google Scholar] [CrossRef]
  20. Deng, M. Cumulative second-harmonic generation of Lamb-mode propagation in a solid plate. J. Appl. Phys. 1999, 85, 3051–3058. [Google Scholar] [CrossRef]
  21. Pei, J.; Deng, M. Assessment of accumulated fatigue damage in solid plates using nonlinear Lamb wave approach. Appl. Phys. Lett. 2007, 90, 121902. [Google Scholar]
  22. Shah, A.; Ribakov, Y. Non-linear ultrasonic evaluation of damaged concrete based on higher order harmonic generation. Mater. Des. 2009, 30, 4095–4102. [Google Scholar] [CrossRef]
  23. Yan-Xun, X.; Fu-Zhen, X.; Ming-Xi, D. Evaluation of Thermal Degradation Induced Material Damage Using Nonlinear Lamb Waves. Chin. Phys. Lett. 2010, 27, 016202. [Google Scholar] [CrossRef]
  24. Khelladi, H.; Rahmi, F. Perturbation methods for the spectral analysis of a weakly nonlinear acoustic field generated by a transient insonation. J. Acoust. Soc. Am. 2013, 19, 045079. [Google Scholar]
  25. Liu, G.; Jayathilake, P.G.; Khoo, B.C. Perturbation method for the second-order nonlinear effect of focused acoustic field around a scatterer in an ideal fluid. Ultrasonics 2014, 54, 576–585. [Google Scholar] [CrossRef] [PubMed]
  26. Momani, S.; Odibat, Z. Homotopy perturbation method for nonlinear partial differential equations of fractional order. Phys. Lett. A 2007, 365, 345–350. [Google Scholar] [CrossRef]
  27. Blackstock, D.T. Convergence of the Keck-Beyer Perturbation Solution for Plane Waves of Finite Amplitude in a Viscous Fluid. J. Acoust. Soc. Am. 1966, 39, 411–413. [Google Scholar] [CrossRef]
  28. Born, M. Dynamical Theory of Crystal Lattices. Am. J. Phys. 1955, 23, 474. [Google Scholar] [CrossRef]
Figure 1. The dislocation dipole model: (a) CT specimen (b) dislocation dipole.
Figure 1. The dislocation dipole model: (a) CT specimen (b) dislocation dipole.
Materials 12 00607 g001
Figure 2. Nonlinear spring model.
Figure 2. Nonlinear spring model.
Materials 12 00607 g002
Figure 3. Relationship between the transmitted coefficient of each order harmonic (the fundamental wave T1, the second-order harmonics T2, the third-order harmonics T3, the fourth-order harmonics T4) and spring stiffness.
Figure 3. Relationship between the transmitted coefficient of each order harmonic (the fundamental wave T1, the second-order harmonics T2, the third-order harmonics T3, the fourth-order harmonics T4) and spring stiffness.
Materials 12 00607 g003
Figure 4. Amplitude variation of (a) the fundamental wave A1, (b) the second-order harmonics A2, (c) the third-order harmonics A3, and (d) the fourth-order harmonics A4 during the fatigue process, which is expressed by dislocation density for early stage and spring stiffness for the later stage.
Figure 4. Amplitude variation of (a) the fundamental wave A1, (b) the second-order harmonics A2, (c) the third-order harmonics A3, and (d) the fourth-order harmonics A4 during the fatigue process, which is expressed by dislocation density for early stage and spring stiffness for the later stage.
Materials 12 00607 g004aMaterials 12 00607 g004b
Figure 5. (a) The experiment system. (b) The stress distribution during the fatigue test.
Figure 5. (a) The experiment system. (b) The stress distribution during the fatigue test.
Materials 12 00607 g005
Figure 6. (a) The fatigue test. (b) The fatigue sample with a crack.
Figure 6. (a) The fatigue test. (b) The fatigue sample with a crack.
Materials 12 00607 g006
Figure 7. (a) Transmitted signals received from the fatigue sample (Loading Cycle 0) in the time-domain and (b) in the frequency-domain (excitation voltage 200 V).
Figure 7. (a) Transmitted signals received from the fatigue sample (Loading Cycle 0) in the time-domain and (b) in the frequency-domain (excitation voltage 200 V).
Materials 12 00607 g007
Figure 8. (a) Transmitted signals received from the fatigue sample (Loading Cycle 5000) in the time-domain and (b) in the frequency-domain (excitation voltage 200 V).
Figure 8. (a) Transmitted signals received from the fatigue sample (Loading Cycle 5000) in the time-domain and (b) in the frequency-domain (excitation voltage 200 V).
Materials 12 00607 g008
Figure 9. (a) Transmitted signals received from the fatigue sample (Loading Cycle 20,000) in the time-domain and (b) in the frequency-domain (excitation voltage 200 V).
Figure 9. (a) Transmitted signals received from the fatigue sample (Loading Cycle 20,000) in the time-domain and (b) in the frequency-domain (excitation voltage 200 V).
Materials 12 00607 g009
Figure 10. (a) Transmitted signals received from the fatigue sample (Loading Cycle 29800) in the time-domain and (b) in the frequency-domain (excitation voltage 200 V).
Figure 10. (a) Transmitted signals received from the fatigue sample (Loading Cycle 29800) in the time-domain and (b) in the frequency-domain (excitation voltage 200 V).
Materials 12 00607 g010
Figure 11. Amplitude variation of (a) the fundamental wave A1, (b) the high order harmonics An (excitation voltage 150 V), (c) the fundamental wave A1, (d) the high order harmonics An (excitation voltage 175 V), (e) the fundamental wave A1, and (f) the high order harmonics An (excitation voltage 200 V) during the whole fatigue process.
Figure 11. Amplitude variation of (a) the fundamental wave A1, (b) the high order harmonics An (excitation voltage 150 V), (c) the fundamental wave A1, (d) the high order harmonics An (excitation voltage 175 V), (e) the fundamental wave A1, and (f) the high order harmonics An (excitation voltage 200 V) during the whole fatigue process.
Materials 12 00607 g011aMaterials 12 00607 g011b
Table 1. Computational parameters of the aluminum alloy.
Table 1. Computational parameters of the aluminum alloy.
ParametersValueUnit
A 2 l 109GPa
A 3 l −510GPa
G 28.6GPa
b0.4mm
Ω 0.33/
R0.33/
h 5nm
v 0.33/
c6250m/s
ρ 2700kg/m3
f1MHz

Share and Cite

MDPI and ACS Style

Lyu, W.; Wu, X.; Xu, W. Nonlinear Acoustic Modeling and Measurements during the Fatigue Process in Metals. Materials 2019, 12, 607. https://doi.org/10.3390/ma12040607

AMA Style

Lyu W, Wu X, Xu W. Nonlinear Acoustic Modeling and Measurements during the Fatigue Process in Metals. Materials. 2019; 12(4):607. https://doi.org/10.3390/ma12040607

Chicago/Turabian Style

Lyu, Wenhan, Xianmei Wu, and Weijiang Xu. 2019. "Nonlinear Acoustic Modeling and Measurements during the Fatigue Process in Metals" Materials 12, no. 4: 607. https://doi.org/10.3390/ma12040607

APA Style

Lyu, W., Wu, X., & Xu, W. (2019). Nonlinear Acoustic Modeling and Measurements during the Fatigue Process in Metals. Materials, 12(4), 607. https://doi.org/10.3390/ma12040607

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop