A Review of Damage, Void Evolution, and Fatigue Life Prediction Models

Degradation, damage evolution, and fatigue models in the literature for various engineering materials, mostly metals and composites, are reviewed. For empirical models established under the framework of Newtonian mechanics, Gurson–Tvergaard–Needleman (GTN) type model, Johnson-Cook (J-C) type damage model, microplasticity model, some other micro-mechanism based damage models, and models using irreversible entropy as a metric with an empirical evolution function are thoroughly discussed. For Physics-based models, the development and applications of unified mechanics theory is reviewed.


Introduction
The degradation, damage evolution, and fatigue behavior of materials are closely related to structural performance and safety. It is well-understood that engineering materials (such as metals and composites) have different micro-mechanisms, degradation processes, damage accumulation, and different failure modes dependent on many factors. For example, when the strain rate is around 10 −6 to 10 −5 s −1 , creep can be a dominant mechanism; when around 10 −4 to 10 −3 s −1 , it is defined as a quasi-static process; when above 10 3 s −1 , it is usually regarded as a high strain rate, where inertia effects, thermal effects (e.g., adiabatic shear banding) and wave propagation influences must be taken into account [1].
In the following, we classify the damage, void evolution and fatigue life prediction models in the literature into two categories: 1. Empirical models are established under a Newtonian mechanics framework (this also includes Hamiltonian and Lagrangian mechanics). Regardless of the different techniques used to characterize the damage evolution with equations and some parameters, it primarily relies on the test data curve to fit the empirical void/damage evolution function. Examples include the Gurson-Tvergaard-Needleman (GTN) model [2], Rice-Tracey model [3], Gunawardena model [4], the well-known, strain rate dependent Johnson-Cook (J-C) damage model [5], some micro-mechanism based damage models, and models using irreversible entropy as a metric with an empirical evolution function. These empirical models are popular for engineering applications due to their simplicity. Still, the identification of parameters is costly, time-consuming, and lacks any scientific basis due to a lack of physical and mathematical foundations [6]. 2. The physics-based models, on the other hand, as the name suggests, are based on the physical foundations and do not require curve fitting empirical functions to test data. They can be classified under the framework of Unified Mechanics Theory (UMT), incorporating the second law of thermodynamics directly into Newton's laws at the ab-initio level [7]. As a result, governing differential equations of any system automatically include energy loss, entropy generation, and the system's degradation in a non-empirical way. The UMT based models are pure physics-based and do not need curve fitting to any test data for the evolution of void/damage. However, they do require deriving analytical thermodynamic fundamental equations of the material without curve fitting. Thermodynamic fundamental equation and the second law of thermodynamics controls the evolution of damage along the fifth axis (Thermodynamic State Index axis), according to Boltzmann equation.
Some recent models for various engineering materials are reviewed in the following sections. They are categorized based on their approach and the type of material being investigated. In Section 2, empirical curve fitting models for metals including Gurson-Tvergaard-Needleman (GTN) type models, Johnson-Cook (J-C) type models, microplasticity models, and some other empirical models are discussed. Their equations are presented. In Section 3, some empirical curve fitting models developed for non-metallic materials are discussed. In Section 4, entropy based models are introduced and presented, including 1. Models using entropy as a damage metric with an empirical evolution function, and 2. Pure physics-based models using unified mechanics theory. Experimental verifications of physics-based models without simulations are also included. In Section 5, the effects of mechanical treatment including manufacturing techniques and surface finishing on surface integrity and fatigue performance are discussed.

Gurson-Tvergaard-Needleman (GTN) Model and Its Modified Form
GTN model [2,8] is a micromechanical empirical model based on the void nucleation, growth and coalescence in materials. Although it is widely used for the prediction of material ductile fracture behavior, it is limited to relatively high-stress conditions and gives inaccurate results when shear is involved [9][10][11]. For the strain rate effect, the incorporation of strain rate related terms on void growth and nucleation in the GTN model only allows the calculation for relatively low strain rate [6]. Therefore, it is incapable of solving high strain rate problems.
Considering these drawbacks, many researchers have adopted the GTN model and performed further modifications such as: 1. Including the damage mechanism due to shear, 2. Improving the applicability of the model under low stress triaxiality, 3. Combined with Johnson-Cook model to better describe the material ductile fracture process under high strain rate of loading. In the following sub-section, the original GTN model and the modified GTN models are introduced. Acharya and Dhar [12] used the GTN model to predict the ductile failure of pipe. The GTN model's empirical curve fitting parameters were properly chosen after the experimental results were obtained. Using the obtained parameters, they developed an elasticplastic finite element code with GTN model for void nucleation and growth. The load versus load point displacement and the load versus crack growth of a pipe under four point bending is well computed by their FEM code incorporating GTN model. The GTN model they used can be expressed in the following form.
The yield function equation is given by where , , are empirical curve fitting parameters, is the von Mises equivalent stress, is the hydrostatic stress, is the flow stress that characterizes the microscopic stress state of the matrix. * is the total effective void volume fraction proposed by Tvergaard and Needleman to account for the onset of void coalescence [8]: where the void volume fraction serves as the damage parameter, is the critical value for void coalescence and is the failure value when the material loses its loading capacity. and are obtained by empirical curve fitting to a test data. The increment of equivalent plastic strain is given by The evolution of void volume fraction is composed of the void growth rate and the void nucleation rate, which is expressed as follows where is a function of in some statistical sense. The previous GTN model, however, is governed only by the first and second invariants and the void fraction of the material. It does not take the void shearing effect which depends on the third stress invariant into account, hence will fail to capture the damage growth under intense shear loading such as micro-/macro-indentation tests. Xue [9] modified the original empirical GTN model by incorporating the void shearing effect into damage evolution function, and used this model to study the ductile fracture behavior of porous materials. In his work, a new void shearing related damage variable was introduced to properly consider the additional damage.
The damage evolution is given by where , , are empirical constants, is an empirical damage rate coefficient, is the void volume fraction, is damage associated with void shearing, accounts for the azimuthal dependence on an octahedral plane, is the equivalent strain. The same approach was used by Youbin Chen et al. [13], where they studied the damage induced by spherical indentation deformation by using a modified GTN model.
Jin et al. [14] established a modified GTN model that incorporates thermal cycling effect into the evolution of void-damage to investigate the damage mechanism of 5A06 aluminum alloy welded joint subjected to thermal cycling. It was shown from the microstructural and fractographic observation results that the void nucleation around the second phase particle is the dominant reason for performance degradation. The modified empirical evolution law of Gurson's void nucleation equation is therefore obtained by implementing the theoretical results from micromechanical analysis using a spherical unit cell model containing particle. However, only the strain-controlled nucleation mechanism is considered in the study. This model cannot be applied to mechanisms such as stresscontrolled void nucleation or fatigue controlled void nucleation.
The increment of equivalent plastic strain in the matrix due to thermal cycling is given by where is Young's modulus of the alloy matrix, the radius of the spherical unit cell model containing a second phase particle, the radius of the plastic region. A zone of reversed plastic flow will form under thermal cycling.
is the radius of the reversed plastic region, is the cyclic period, is the yield stress in the interface and the neighboring matrix material.
The modified void nucleation equation has the same form as Equation (4), however ̇ here is obtained from Equation (6).
The standard local continuum models have a high sensitivity on spatial discretization length when it is used to describe ductile damage, which leads to the high mesh sensitivity of FEM for crack problems. However, in ductile-brittle transition region damage can precede to final failure by cleavage. In that case models require very fine discretization to resolve the high stress gradients at crack tips in order to evaluate the probability of cleavage fracture at a certain load state, which conflicts with the mesh sensitivity in local ductile damage models. Linse et al. [15] used a GTN model based on non-local modification to investigate the evolution of ductile damage and stress state at the crack tip of a typical steel pressure vessel. The non-local GTN model is established by replacing the dilatational part of the plastic strain rate ̇ by its non-local spatial average ̅ in the void growth empirical damage evolution equation. The modified empirical evolution of void volume fraction is given by where ̇ is the non-local modification of void growth rate, ̇ the empirical damage evolution rate, ̇ the void nucleation rate, the volumetric plastic strain, ̅ the nonlocal volumetric plastic strain, the non-local length parameter. The results show that compared to local damage model, the proposed non-local model predicts an initial stage of stable crack tip blunting followed by a distinct point of crack initiation.
As previously mentioned, the GTN model cannot correctly capture the damage mechanism due to shear. Among all the research works that attempted to incorporate shear into GTN's empirical damage evolution equation, Xue's shear mechanism [9] and Nahshon-Hutchinson's shear mechanism [10] have received the most attention. Xu et al. [16] applied the modified GTN model based on Nahshon-Hutchinson's shear mechanism to simulate the failure behavior of the clinched joint on the Al6061 sheet, because it was observed that during the clinching process, ductile fracture is a main failure mode that belongs to shear-dominated failure. The modified GTN empirical damage evolution equation for void growth is given by where is the deviatoric stress, ( ) the non-dimensional metric of stress state, the von Mises equivalent stress, is an empirical shear damage coefficient. Using the same modified model, Gatea et al. [17] simulated the ductile fracture in Single Point Incremental Forming (SPIF) process due to void nucleation and coalescence, and results are compared with the original GTN model. It was shown that Nahshon-Hutchinson's mechanism has better accuracy on fracture than Xue's mechanism under tensile/shear loading conditions.
Malcher et al. [18] proposed a modified GTN model to improve the limitation of original GTN model on low stress triaxiality. The modified GTN model is established by introducing two new empirical damage parameters, one is related to the hydrostatic stress component and the other is related to the deviatoric stress component.
where , , are empirical constants, is a Lode angle dependence function, ̇ is the evolution of the damage, ̇ is the evolution of void nucleation and ̇ is the evolution of shear effects, represents the fraction of all second-phase particles with potential for nucleation, and are the mean strain for void nucleation and its standard deviation. The variable ̅ represents the equivalent plastic strain and ̅ is the rate of the accumulated plastic strain. is the stress triaxiality parameter, is a numerical constant that needs to be calibrated for each material by curve fitting, is the ligament size ratio defined for two or three-dimensional problems. The two damage parameters affect the evolution of internal variables and allows more accurate values at the time of crack formation. The new micro-defects nucleation mechanism improved the performance for a wide range of stress triaxiality. Furthermore, the stress triaxiality dependence in the evolution of shear damage parameter also enhanced the prediction of the fracture location under combined loading conditions. Wang et al. [19] also used a modified GTN model that considers shear induced damage to analyze the tearing failure of ultra-thin sheet-metal, including the size effect in blanking process. In comparison to Nahshon-Hutchinson's shear mechanism [10] that extends the GTN model by incorporating the third stress invariant into the damage evolution function to distinguish the shear dominated states, the model used here extends the GTN model by coupling the volumetric damage and shear damage into yield function and flow potential, respectively. This modified GTN model is given by In which is the von Mises equivalent stress, the hydrostatic stress, the flow stress that characterizes the microscopic stress state of the matrix, * the total effective void volume fraction.
is an empirical shear damage parameter, the equivalent plastic strain of material matrix, the failure strain under pure shear state, is an empirical weakening exponential, ( , * ) is an empirical weight factor.
Bambach and Imran [20] proposed a modified GTN model that considers hot forming, in which the damage is not only controlled by stress state but also the softening process-including recovery and recrystallization. The proposed empirical void nucleation function which couples softening processes at inclusions is implemented into the GTN empirical damage evolution equation. Their empirical void nucleation evolution function is given by where ̇ is the void nucleation rate, the fracture toughness, ̇ the plastic strain rate, is a normalized third stress invariant, , , , and are curve fitting parameters. The derived extended GTN model is able to describe void formation at elevated temperatures and its retardation when softening occurs.
The GTN model can also be combined with various constitutive models to predict the material behavior of interest. Chen et al. [21] investigated a single impact loading process for twinning induced plasticity (TWIP) in steels using GTN damage model combined with the dislocation density-based viscoplastic constitutive model. Results show that the induced residual stress, plastic strain, and damage fields can be predicted. The fracture analysis of shape memory alloys based on GTN model's voids nucleation/growth in martensite and austenite phase has also been studied recently [22,23]. Bahrami et al. [23] proposed a constitutive model to investigate the pseudoelastic-plastic behavior of the shape memory alloys, SMAs, up to fracture. The proposed model based on the Boyd and Lagoudas phase transformation model [24] is extended to take the plastic deformation and the fracture behavior of the SMAs into account by applying the GTN model shown in Equations (1)-(4).

Johnson-Cook (J-C) Type Models
An empirical constitutive relation developed by Johnson and Cook (J-C) [5], or the so-called J-C model, is widely used to simulate the damage evolution and predict failure in many strain rate sensitive engineering materials. The Johnson-Cook model is a plasticity model that is based on von-Mises plasticity [25] with closed-form analytical equations specifying the hardening behavior and the strain-rate dependence of the yield stress. Compared to the micromechanics-based GTN model which only allows the calculation for relatively low strain rate, Johnson-Cook model enables the calculation of the material response under high strain rate (above 10 3 s −1 ) deformation for most metals. However, this model is incapable of properly capturing the static and dynamic recovery, and reflecting the effects of load path and strain-rate history in large deformation processes [26]. In the following, some research works based on Johnson-Cook models and their modified forms are introduced.
Laser shock peening (LSP) is a surface treatment process for airfoils that is achieved by the induction of compressive stress. Improper operating during LSP process may cause formation of small subsurface cracks in the processed material. Li et al. [27] simulated the formation and predicted the sizes of cracks generated by inappropriate laser shock peening (LSP) processing in airfoil specimens by a 3D FEM model. Their model was developed based on Johnson-Cook plastic and fracture models to properly calculate the plastic and fracture behaviors at high strain rate.
The Johnson-Cook empirical plastic model is given by and the Johnson-Cook empirical damage initiation criterion is given by where is the stress; the plastic strain; the initial yield stress; ̇ the reference strain rate; the reference temperature; * the dimensionless pressure-stress ratio; ̅ ̇ the plastic strain rate; * the homologous temperature; , , , , and are empirical coefficients; , , , , and represent different empirical failure parameters obtained by curve fitting to test data.
Nam et al. [28] used the Johnson-Cook model to investigate the crack tip stress and strain fields at crack initiation of A106 Gr. B carbon steels under high strain rates. The empirical parameters for Johnson-Cook model are found by fitting the model to the tensile test results at different strain rate. Nam et al. [29] performed a ductile fracture simulation based on fracture strain energy, which is calculated based on the assumption that the fracture strain energy for dimple fracture depends on strain rate, and also depends exponentially on stress triaxiality. The high strain rate effect in the simulation is captured by Johnson-Cook model, and the damage can be calculated by summing the incremental damage ∆ , by where is the fracture strain energy, ∆ is the equivalent plastic strain energy increment, / the stress triaxiality, ̇ the equivalent plastic strain rate, ̇ the reference strain rate, and , , , are empirical curve fitting constants. This model provides the damage initiation criteria for dynamic failure.
Chen et al. [30] proposed a modified Johnson-Cook model that includes the influence of the corrosion to study the mechanical behavior of corroded high strength reinforcing steel bars under static and dynamic loading. Corrosion reduction factors that relate the tensile behaviors with corrosion degree and strain rate are obtained from tensile test results under various strain rates.
The modified model is given by where is the corrosion degree, , , and are curve fitting parameters of corroded steel bars, , , and are curve fitting parameters of the Johnson-Cook model of the uncorroded steel bars. , , are curve fitting parameters also determined from test data. It was found that the strain rate effect of the yield strength and ultimate strength decreases with the increase of corrosion rate.
Chen et al. [31] developed a coupling J-CM (modified Johnson-Cook) plastic and energy density based damage model to characterize the ductile damage behaviors of Ti-6Al-4V alloy, including plastic, failure initial strain and damage evolution, for high strain rate compression tests under wide strain rate and temperature range. The modified Johnson-Cook model is achieved by incorporating an empirical temperature dependent function to the original Johnson-Cook model, as follows = ( + ( ) ) 1 + × ln where , , , , and are curve fitting empirical parameters. ( ) is an empirical temperature dependent function to account for the temperature dependent work hardening behavior in flow curves.
is the equivalent plastic flow stress, ̇ the equivalent plastic strain rate, ̇ the reference strain rate, and the workpiece ambient and melting temperature, respectively. is a critical temperature related to microstructure evolution. In this model, the temperature dependent work hardening behavior in flow stress is characterized by the empirical temperature dependent term related to microstructural evolution.
Wang et al. [32] investigated the deformation and fracture behaviors of superalloy GH3536 over a wide range of temperatures (298 K-1073 K), strain rates (0.1 s −1 -5200 s −1 ) and also stress triaxiality (0.6-1.1). A modified Johnson-Cook model was developed to overcome the limitation of original Johnson-Cook model which cannot characterize the temperature and strain rate dependences of the deformation behavior, and the stress triaxiality, temperature and strain rate dependences of the fracture behavior of the superalloy. The modified Johnson-Cook model is given by where , , , , , , , , , and are curve fitting parameters, ̇ * = / ̇ is the dimensionless strain rate, * = ( − )/( − ) is homologous temperature. The proposed model successfully predict the impact resistance of the superalloy honeycomb which cannot be captured by the traditional J-C constitutive model and J-C fracture criterion.
Wang and Hassani [33] investigated the deformation of the spherical microparticles of pure titanium impacting a rigid substrate at an ultra-high strain rates (10 6 -10 10 s −1 ). The impact deformation was recorded in real-time. The simulation was conducted by finite element approach using two constitutive equations: Johnson-Cook and Zerilli-Armstrong [34]. By comparing the deformed geometries from experimental data and simulated results, they evaluated the capability of the two constitutive equations-originally calibrated at 10 3 -10 4 s −1 to describe deformation at ultra-high strain rates. Being mechanistically based, the Zerilli-Armstrong model was found to have a better performance than the Johnson-Cook model at higher strain rates because the extrapolation of the Johnson-Cook strain rate sensitivity to a higher strain rate leads to a significant deviation from the experimental data. The Zerilli-Armstrong model is given by where is the equivalent plastic strain, ̇ the plastic strain rate, the temperature, , , , , , , , are curve fitting parameters. Here, is the athermal component of the material yield strength. and are strain hardening constants. and are thermal softening parameters. Finally, and are strain rate sensitivity parameters. In comparison, both Johnson-Cook model and Zerilli-Armstrong model are plasticity models that incorporate the dependence of flow stress on the plastic strain, strain rate, and temperature. Despite both of them needing to calibrate empirical parameters, the Zerilli-Armstrong model is more micro-mechanism based. The Zerilli-Armstrong model is established based on the thermally activated dislocation motion of the material undergoing plastic deformation. It uses an Arrhenius type relation to describe the probability of dislocations to overcome barriers and obstacles. The Zerilli-Armstrong model considers that the activation volume of face-centered cubic (FCC) metals is dependent on strain while that of body-centered cubic (BCC) metals is independent of strain. For HCP titanium crystal tested in this study, the ZA model uses the superposition of FCC and BCC approach.
The dynamic mechanical responses of metallic alloys are widely described by the Johnson-Cook model and the Khane-Huange-Liang (KHL) [35,36] model. KHL model is a constitutive viscoplastic model that combines effect of stress, strain frequency, and temperature to determine work-hardening behavior. Zhang et al. [37] proposed a combined JC-KHL (Johnson Cook-Khane Huange Liang) model (or CJK) to model constitutive response, such as true strain rate evolution and true stress-strain data, for the necking evolution of near Ti3Al2.5 V alloy at high strain rates. By using the proposed CJK model, the decreasing strain hardening with strain rate, and logarithmic strain rate dependent flow stress, with the reference strain rate 0.001 s −1 is described. The CJK model is given by where is the true stress, the plastic strain, ̇ the reference strain rate, and the upper bound strain rate chosen arbitrarily. , , , , , and are the curve fitting constants. , , and are current, reference, and melting temperatures.
Chiyatan and Uthaisangsuk [38] investigated the strain rate effects on mechanical properties and fracture mechanism of ferritic-martensitic dual phase (DP) steels. The microstructure effects on local deformation and damage of DP steels under varying strain rates are simulated by FEM using 2D representative volume elements (RVEs). The macroscopic flow curve is obtained from the coupling of Swift-Voce hardening law [39,40] and Johnson-Cook rate-dependent model, as follows where , , n, B, Q, β in the first term of Equation (21) are the Swift and Voce curve fitting material parameters, and parameter 0 ≤ ≤ 1 is a weighting coefficient. The second term is Johnson-Cook strain rate hardening equation including the parameters C and the reference strain rate ̇ . The flow curve of phase constituents at different strain rates were characterized by dislocation-based theory and local chemical composition combined with the Johnson-Cook hardening model, as follows where and are the von Mises stress and equivalent plastic strain, respectively. The first term represents the Peierls stress [41] and the effect of alloying elements in the solid solution state. The second term, ∆ described the material strengthening by precipitation or carbon in solution. The last term demonstrated the effect of dislocation strengthening and material softening, which contains the material constant , Taylor factor [42] , shear modulus , Burger's vector for both phases. is the recovery rates and is the dislocation mean free path of each phase.
Like GTN model, Johnson-Cook model can also be combined with various constitutive models to predict the material behavior of interest. Jeunechamps and Ponthot [43] proposed an elasto-viscoplastic model based on the coupling between the Johnson-Cook model and the Perzyna viscosity model [44] to perform a thermomechanical simulation of elastic-viscoplastic materials subjected to high strain rate. Wang et al. [45] proposed model based on the combination of Johnson-Cook and GTN models to characterize the ductile fracture process of steel and predict structural damage during ship collisions and grounding accidents.

Micro-Plasticity Models
Besides the studies in macroscopic or continuum scale, many studies focus on micro plasticity in mesoscale or the dislocation motions and atomic vacancy generation/annihilation in microscale to unveil the fatigue hotspots, fatigue crack nucleation, crack path, and fatigue life prediction. Still, these micro-mechanism-based models need the calibration of material coefficients, and the damage evolution function needs to be obtained by curve fit experimental data. They are introduced in the following.
Wan et al. [46] used the high resolution electron backscatter diffraction (HR-EBSD) to investigate the microstructural residual stress distributions and important stress states for defect nucleation on a deformed copper polycrystal, then utilized computational crystal plasticity to capture the microstructural residual stress components, effective stress, hydrostatic stress and stress triaxiality and compared with experimental results. The fatigue crack nucleation equation is established based on a stored energy criterion.
The empirical microstructure-sensitive fatigue crack nucleation equation is given by This criterion is established based on the local slip activity and local storage volume Δ due to the accumulation of geometrically necessary and statistically stored dislocation . is a fraction for stored plastic energy and accumulated crystal slip used to calculate the local stored energy per cycle. This criterion for fatigue crack nucleation indicates that preferential sites for fatigue crack nucleation are local to grain boundaries, and that hard-soft grain interfaces where high GND densities develop are preferable.
Wilson et al. [47] utilized both experiments and integrated crystal plasticity eXtended Finite Element (XFEM) modelling to study the fatigue crack growth in titanium alloy, ferritic steel, nickel superalloy, and zirconium alloy (zircaloy 4) with BCC, FCC, and HCP crystallography. The crack propagation was found to be controlled by crack tip stored energy, and the crack direction by anisotropic crystallographic slip at the crack tip. The fatigue crack path tortuosities and growth rates in the materials are also captured.
The crystal plasticity model is as follows: where is the deformation gradient which can be decomposed into elastic part and plastic part . is the plastic velocity gradient, is the total number of slip systems, ̇ is the slip rate on slip system , and and are the corresponding normal slip plane and slip direction.
is the density of mobile dislocations, the frequency of attempts of dislocations to jump obstacle energy barriers, the Burgers vector, Δ the thermal activation energy, the Boltzmann constant, the temperature (295 K), and the resolved shear stress and critical resolved shear stress on slip system respectively, and Δ is the activation volume, is the hardening coefficient, ̇ the rate of accumulated plastic strain, and the time increment. They used empirical fatigue crack nucleation equation given by Equation (23). The empirical crystal plasticity model shown in Equation (24) is essentially a constitutive and medium mechanics theory based on crystalline material's dislocation slip mechanism. It is widely used to describe the anisotropic, texture evolution, the twin deformation, damage fracture, creep, recrystallization in polycrystal materials. Bandyopadhyay et al. [48] proposed that the microstructure-sensitive critical plastic strain energy density (SPSED) is the driving mechanism of fatigue crack initiation and can be a parameter for fatigue life prediction. They used the crystal plasticity finite element simulations to compute the (local) SPSED at each material point within polycrystalline aggregates of a nickel-based superalloy. The obtained SPSED is calibrated and then used to predict fatigue life at various strain ranges. The empirical fatigue life prediction model is given by where is a set of empirical parameters defining the loading conditions, such as the applied strain range, strain ratio, temperature; ( , * ), Δ ( , * ) and ( ) are obtained from the crystal plasticity finite element simulations, and is the parameter calibrated from laboratory test data.
The evolution of fatigue damage can be ascribed to multiple micro-mechanisms related to atomic vacancy accumulation and dislocation motions. A quantified value , called slip irreversibility, can be found in literature as a material parameter to incorporate those micro-mechanisms in fatigue life predication models. Mughrabi [49] proposed a microstructure-based reformulation of the Coffin-Manson fatigue life law, which includes an empirically determined slip irreversibility. Ho et al. [50] used the slip irreversibility based Coffin-Manson fatigue law to study the relation between fatigue crack initiation and cumulative slip irreversibility of nickel-base alloys. Their empirical equation is given by where cum is the cumulative slip irreversibility, defined as the ratio of microstructurally irreversible cyclic plastic strain with respect to the total cyclic plastic strain. is the interband spacing, ℎ is the mean extrusion height, is the average grain size, is the volume fraction of one slip band, is the Taylor factor [42], and are the fatigue ductility coefficient and exponent.

Other Empirical Models for Metals
Besides the above mentioned GTN model, Johnson-Cook model, micro-plasticity models, and their modified versions, many other empirical models such as energy based model, Lemaitre's damage model, and various micro-voiding models can capture the mechanical response of metals under various strain rates. They are discussed in the following section.
Dondeti et al. [51] performed a computational analysis of ductile fracture in heterogeneous porous ductile materials containing brittle inclusions based on a rate-dependent homogenization-based continuum plasticity damage (HCPD) model. The rate-dependent HCPD model has a similar structure to the GTN model for yield function and void growth but a novel empirical void nucleation criterion that can simulate the decreasing load-carrying capacity in the material to account for inclusion cracking and void growth. The proposed empirical void nucleation criterion is given by is a material parameter that is calibrated from micro-mechanically simulated volume fractions of cracked inclusions at a given strain, is a function related to the rate of evolution of the area fraction of cracked inclusions.
The Energy-based model is another common approach. For example, Darras et al. [52] observed the damage evolution of 5083 Aluminum alloy at various strain rates. They calculated the accumulated plastic strain and plastic strain energies from the true stresstrue strain curves. Then, an energy-based empirical model was used to predict the damage evolution of the 5083 Aluminum alloy at various strain rates. Abed et al. [53] performed a mechanical study on EN08 steel at quasi-static and dynamic strain rates. They quantified the micro defects density by scanning electron microscope (SEM) images and proposed it as an input for the energy-based model to predict the material's damage evolution. Their energy-based empirical damage model is given by where is the damage during deformation, the damage at fracture obtained by SEM images, the dissipated energy, the total dissipated energy, and is an empirical constant obtained by curve fitting to test data to determine the damage evolution trend throughout deformation. However, the accuracy of this SEM image-aided energy-based model proposed in Equation (28) is affected by the loading conditions. For example, it cannot accurately estimate the damage (such as in drop hammer dynamic test) because the specimen is not fractured. The SEM image obtained by cutting the specimen's section is therefore inaccurate.
Khoei et al. [54] used a rate-dependent damage-viscoplasticity model to simulate the crack growth in ductile materials under dynamic loading. The model is established by eliminating the yield surface concept and modifying the plastic multiplier to consider the damaging effect, as follows where ̇ is the plastic multiplier, and Ξ are empirical material constants, the von-Mises effective stress, the material yield strength, ̅ is the equivalent plastic strain, is an empirical damage parameter.
Shojaei et al. [55] proposed a viscoplastic constitutive theory for brittle to ductile damage in polycrystalline materials under dynamic loading. They first developed a viscoplastic model to properly obtain the material response between low and high strain rate, then they correlated the microscopic degradation process (microcracking, microvoiding) to macroscopic failure modes (ductile or brittle). A fracture mechanism-based damage model is established to capture the microcracking process when various dynamic energy densities are applied to polycrystalline materials. This model is suited for deviatoric stress dominant problems.
The proposed micro-void nucleation and growth are given by where ̇( ) is the total void volume fraction evolution rate, ̇ ( ) is the rate of void nucleation, ̇ ( ) is the rate of void growth. and are empirical curve fitting parameters, is hydrostatic threshold stress for microvoid nucleation, is stress for microvoid growth, Σ is the real damage stress, is the Boltzmann constant, is a viscosity parameter.
Chen et al. [56] proposed a reliability assessment model focusing on low-cycle fatigue and high energy impact loads. The fatigue damage with time considering impact effect was developed. A modified Coffin-Manson equation for analyzing the effect of impact loads on fatigue damage was established by data fitting, as follows Equation (31) is the modified empirical Coffin-Manson equation, in which , , , and are the empirical curve fitting constants, is the original static tensile strength, the original static fracture ductility, is the plastic strain amplitude. The impact-affected fatigue can be obtained from Equation (31) by iteration method. Carniel et al. [57] analyzed the geometrically nonlinear trusses associated with viscoelastic and viscoplastic members considering mechanical degradation, inertia, high strain rate deformations, and strain rate related effects. Lemaitre's damage model [58] was applied to capture the material's damage evolution under a high strain rate.
Lemaitre's empirical damage evolution law is given by where is the total area of intersection of a given plane with a representative volume element, ̅ the effective resisting area so that the damage variable can assume values 0 ≤ ≤ 1, , and are empirical damage evolution parameters, ̅ the damage threshold and (− ) the damage strain energy density release rate. However, it is important to point that having plastic strain alone as a criterion for damage potential as in Equation (32) violates the second law of thermodynamics. Because if the load is applied continuously below the critical strain rate(for example, one million times), there can be no damage, which is not true. Shen et al. [59] employed a damage-coupled elastoplastic constitutive model considering nonlinear kinematic hardening to evaluate the fatigue life of notched specimens with plastic deformation at the notch tip. At the notch tip, damage is mainly induced by plastic strain while for the material not at the tip the damage is induced by cyclic stress. Therefore, they established a strain-based damage model and a stress-based model, respectively, as follows to properly estimate the damage.
The plastic strain-based damage model is given by The stress-based damage model is given by where * is the maximum value of the damage equivalent stress over a loading cycle, is the elastic modulus. The parameters and are fitted from the experimental curve of plastic strain versus number of cycles to failure. and , are the amplitude of the octahedral shear stress and the mean value of the hydrostatic stress in a loading cycle, respectively. The term , is the maximum equivalent stress over a loading cycle, is the fatigue limit at the fully reversed loading condition, and is the ultimate tensile stress. The five parameters, , , , , and , are determined by using plain fatigue tests of standard specimens. In Equations (33) and (34), the strain-based model depends on the cumulative plastic strain while the stress-based model depends on the stress quantities in one loading cycle.
Tang et al. [60] modified a damage model proposed by Kachanov [61] and later on developed by Lemaitre [62] to predict the formability of high-strength steel sheets at elevated temperatures. The modified Lemaitre-based damage model's empirical material parameters are obtained from tensile test data at temperatures ranging from 550 °C to 850 °C and at various strain rates.
Lemaitre's empirical damage potential is given by whereas the modified damage potential is given by In which and are empirical materials parameters and are functions of the strain rate and temperature, is the damage strain energy density release rate, ̅ is the accumulated plastic strain, is another empirical parameter obtained by curve fitting to test data. The modified Lemaitre damage potential fixes the inaccuracy of classical Lemaitre's potential when describing highly ductile materials.
Wu et al. [63] modeled the low cycle fatigue of 1.4848 cast austenitic steel at temperatures ranging from room temperature to 1173 K and at strain rate from 2 × 10 −4 to 2 × 10 −2 s −1 based on the integrated creep-fatigue theory (ICFT). They observed that at intermediate temperatures, 673 K and 873 K, dynamic strain aging (DSA) promotes slip inhomogeneity with dislocation pile-ups, forming concentrated slip bands or dislocation walls, which are embryos of cracks. Their empirical fatigue damage evolution equation based on dynamic strain aging is given by where Δ is the amplitude of cyclic hardening, the dislocation density level below which there is no instantaneous crack nucleation, the proportional constant for dislocation-nucleated cracks, β the empirical proportional constant for creep damage, and the creep strain. Equation (37) assumes that the reduced factor D for low cycle fatigue life depends on DSA effect and creep deformation.
Phase field method [64] is a continuous interface description method which has been widely used to simulate the dynamic crack propagation and ductile fracture in metals because of the advantage of not using any explicit fracture criterion and easy trace of complex fracture surfaces. Mozaffari and Voyiadjis [65,66] showed that the established non-local gradient type damage model through the phase-field method can couple with a viscoplastic model to assess the inelastic behavior of rate-dependent material. The empirical damage evolution law incorporating the viscoplastic deformation is given by where is effective stiffness tensor, ̅ is effective strain tensor, ̅ is effective viscoplastic strain tensors in undamaged configuration.
is an empirical scalar function to map the stress state between the damaged configuration and effective undamaged configuration for the case of isotropic damage, = ( ) . is the material constant with energy dimensions which contains dissipation during the whole process of damage through elastic and viscoplastic deformations. The constant is considered in the form = to separate the effect of energy type constant ( ) and length unit, in which corresponds to the length scale due to damage. Furthermore, Badnava et al. [67] incorporated the phase-field effect in the Peric's viscoplastic model [68] to simulate the influence of the loading rate on the ductile fracture, Schreiber et al. [69] utilized the framework of phase-field modeling for the fracture to simulate fatigue crack growth.
Chu et al. [70] proposed a unified phase-field damage model to capture the transition of typical high strain rate failure modes of metals, such as dynamic brittle fracture and shear banding. The proposed method distinguished the transition of the material's failure by setting the critical energy release rate and energy density threshold for damage to vary with stress triaxiality. The failure energy excluding the plastic dissipation in the fracture process zone before damage evolution is defined to model the ductile failure more physically. A degradation function of the yield stress is introduced to provide damage softening mechanism for ductile fracture and increase the simulation accuracy of brittle fracture. The failure energy density is given by where is an empirical length scale parameter associated with the regulation of sharp discontinuities, the empirical damage parameter with = 0 defines the intact state, and = 1 defines the fully damaged state of the material.
is the equivalent critical energy release rate corresponding to the evolution of the internal discontinuous boundary, is the critical energy release rate obtained by experiments, and < is the energy release rate failure parameter before damage initiation in the fracture process zone (FPZ). The evolution equation for the phase field is given by For a homogeneous stretch problem of ideal elasto-plastic material where is empirical fraction of plastic work converted to heat, is compression part of the inherent elastic strain energy density, is the inherent plastic stored energy density, is an empirical energy density threshold before the damage initiation. and are one-dimensional strain and inherent stress, respectively.

Empirical Curve Fitting Models for Non-Metallic Materials
The models introduced in earlier sections fail to describe the damage evolution for composite materials. Because in composite materials, matrix and fiber must be treated separately in the simulation, and the damage evolution equation in different directions also differs. In the following, some recent research work focusing on the damage evolution of composites, laminates, cement, and asphalt materials are introduced.

Composite Materials
Zhang et al. [71] performed a high strain rate compression experiment on a 2D plain woven composites (2DPWC) along the thickness direction, and modeled its mechanical behavior by FEM at microstructural level. The failure morphologies of 2DPWCs were found to be different depending on the strain rate of the loading.
They proposed the following empirical criterion for damage initiation: Their empirical damage evolution law for ductile and shear failure is given by where is the stress triaxiality, ̅ ̇ the equivalent plastic strain rate, ̇ the effective displacement rate, and ̇ the maximum value of the effective displacement at the point of failure. After damage initiation following Equation (42), the material stiffness was degraded progressively according to Equation (43). The damage evolution law used in ductile and shear failure was based on the energy dissipated during the damage process, which allowed the removal of elements from the mesh.
Alemi-Ardakani et al. [72] proposed two simplified empirical models incorporating correction factors to perform the fast simulation of out-of-plane impact response of fiberreinforced polymer composites (FRPs). The empirical model considered: (a) strain rate dependency of the mechanical properties, (b) difference between tensile and flexural bending responses, (c) delamination, and (d) the geometry of fixture. The first approach is achieved by applying multiple correction factors to the quasistatic material properties, while the second approach is by applying only one single correction factor over all material properties at once. Results are validated with the simulation by Abaqus with built in Hashin's empirical damage criterion [73,74]. The Hashin's empirical damage criterion considered the progressive damage of the fiber reinforced composites. This criterion considers fiber and matrix failure in tensile and compression.
Hashin's failure criterion is given by , , ̂ are the inplane normal and shear stresses (the 1-direction denotes the fibers direction). The empirical coefficient defines the contribution of the shear stress to the fiber tensile failure initiation. In this model it is assumed that fiber or matrix fails if ultimate stress is exceeded. However, applying a stress little below the critical stress million times doesn't lead to fatigue in this model.
Chen and Morozov [75] proposed a consistency elasto-viscoplastic damage model to capture composite material's rate-dependent plastic response and progressive post-failure behavior. The yield criterion is rate-dependent in the model, and the standard Kuhn-Tucker conditions [76] govern the plastic loading and unloading. The use of this rate-dependent criterion ensures that the Kuhn-Tucker plastic loading/unloading conditions and plastic consistency condition remain valid for the proposed viscoplastic model The damage initiation and propagation criterion in their work is given by where is the loading function in the form of Hashin's failure criterion, is an empirical damage threshold corresponding to each failure mechanism that controls the size of the expanding damage surface depending on loading history, is an empirical parameter that defines the exponential softening law. The initial value of is 1; subscripts 1 and 2 represent the fiber and transverse directions of the unidirectional ply; subscripts and denote tension and compression. The damage variable represents the damage effects on the shear stiffness due to matrix fracture caused by a combined action of transverse and shear stresses. The damage initiation and propagation criterion predict the onset and further evolution of each damage mechanism, and the damage evolution laws govern the evolution of damage variables.
Park et al. [77] used a rate-dependent constitutive model based on the multi-scale approach to model the rate-dependent damage of polymeric composites. For dynamic loading, the strain-rate affects the damage behavior of composite materials and the matrix's behavior. They employed an enhanced micromechanical model that improves the in-plane shear behavior to analyze the fiber and matrix constituents' rate-dependent behaviors. The rate-dependent elastic composite damage model is introduced to account for the damage from the micro-cracks in the matrix and the matrix/fiber interface. The empirical damage evolution laws for the rate-dependent damage model are given by where is an empirical damage variable, , , , and are empirical damage constants representing the transverse damage initiation, transversely critical damage, inplane shear damage initiation, and in-plane shear critical damage, respectively.
Seman et al. [78] investigated the high strain rate damage behavior of Kenaf fiberreinforced composite materials based on multi-scale finite element modeling. In the FEM modeling scheme, a micro-scale structure model with periodic boundary conditions was used to homogenize the heterogeneous fiber/resin system into a unit cell. A mesoscale model incorporating with constituent's failure criterion was responsible for the stiffness degradation and subsequent element removal. They employed Hashin's 3D failure criterion for the fiber to predict the composite's failure initiation under a high strain rate, similar to Equations (44) and (45). For the matrix, two different dynamic failure criterion are used, as follows. The ductile criteria for damage initiation is met when: The shear criteria for damage initiation is met when: In which the equivalent plastic strain is given by where and correspond to the equivalent plastic strain at ductile damage initiation for uniaxial tensile and uniaxial compressive deformation, respectively, = 1/3 and = −1/3 are the stress triaxiality in uniaxial tensile and compressive deformation state, is an empirical parameter, = (1 − )/ , = , and correspond to the equivalent plastic strain at shear damage initiation for uniaxial tensile and compressive deformation, respectively. The parameters and correspond to the values of at = and, = respectively. and are curve fitting empirical parameters. Using these structural models together with finite element analysis, they obtained the mechanical properties and also the stress state in Kenaf composites.
Alabdullah and Ghoniem [79] investigated the non-linear mechanical behavior of SiC/SiC composite materials in nuclear applications. The neutron irradiation-induced cyclic thermal and mechanical loading leads to wide-spread and progressive micro-cracking that reduces thermal conductivity and enhances thermo-mechanical damage. An empirical damage model coupling mechanical, thermal, and irradiation damage and considering the loss of thermal conductivity is developed based on continuum thermodynamic damage mechanics for wide-spread micro-cracking.
The total damage from sources mentioned above is given by where empirical parameters * , , , , are the maximum effective thermodynamic force obtained during the loading history, initial damage threshold, thermodynamic normalizing force constant, shaping parameter, and maximum damage obtained, respectively. , , are the maximum thermodynamic force obtained during the radiation history, initial damage threshold, shaping parameter. The parameter is a material constant related to the damage associated with volumetric swelling (micro voids and/or small dislocation loops). The linear decomposition of the total damage to three linearly independent dissipation function with controlled threshold provides numerical benefit when implementing the model. The convergence of all the decomposed damage parameters are based on only the convergence of three total damage scalar variables. However, this model is based on superposition principle, which is only true for entropy. But the author's damage is based on forces.

Laminates
Morinière et al. [80] published a review of modeling of impact damage and dynamics in fiber-metal laminates(FMLs). Sitnikova et al. [81] developed a 3-D nonlinear finite element model to simulate blast failure of fiber metal laminates. In their work, an empirical damage evolution law is incorporated into the composite constitutive behavior to obtain fiber-metal laminate (FML) panels' blast response. Since the simulation of aluminum alloys is relatively simple, they focused their work on modelling the composites. The proposed formulation can be applied to composites with a plain weave architecture and possible through-thickness damage. The empirical damage evolution law for the glass fiber reinforced composites in the model is as follows: where the are the damage variables, namely, and correspond to the failure in warp and weft fiber directions, respectively, describes through the thickness composite crushing failure, and refers to the in-plane shear failure. The subscripts "t" and "c" denote tensile and compressive failure.
are effective stresses, are material strengths. Empirical coefficients and in the above equation govern the rate of growth of damage.

Cement and Asphalt Mortar
Fu et al. [82] investigated the high strain rate effects (impact loading) on the compressive strength, elastic modulus, peak strain, and specific energy absorption of cement and asphalt mortar using a split hopkinson pressure bar (SHPB) test. The results indicated that the compressive strength and specific energy absorption increased with increasing strain rate. They proposed a statistical continuous rate-dependent damage constitutive model for cement and asphalt (CA).
The empirical damage evolution equation is given by where is the total strain of CA mortar, is any strain rate under impact loading, is an empirical strain rate sensitivity index of the damage variable, is an empirical scale parameter related to the strength; is the empirical morphological parameter of the Weibull distribution.

Models Using Irreversible Entropy as a Metric with an Empirical Evolution Function
The development of this type of model can be traced back to 1998, when Basaran and Yan [83] proposed using entropy as a damage metric. And in 1999, Chandaroy and Basaran [84,85] used energy dissipated during thermal cycling in a microelectronic solder joint to model fatigue life. It was found that the total energy dissipation (cumulative entropy generation) for a specimen under test is a material property. It differs according to the material but is a constant for the same material at failure. It is also independent of loading conditions. Therefore, entropy was proposed to be a metric for damage estimation in many research works since then. However, among these studies many of them are still using the empirical evolution function that need to be calibrated for the empirical coefficients. In this section, we discussed the models using irreversible entropy as a metric with an empirical evolution function. The focus is on their entropy generation (intrinsic dissipation) equation, damage evolution equation, and entropy based application.
In recent years, using entropy as a damage metric has caught more attention. Naderi et al. [86] used the entropy generated during the specimen's fatigue life to characterize the material degradation. Results indicated that the cumulative entropy generation is constant at failure, independent of material geometry, loading type, and loading frequency. The value of this constant, or called the fatigue fracture entropy (FFE), varies with material, implying that it is a material property. For example, stainless steel and Al have different cumulative entropy generations at the fracture point (different FFE). Figures 1 and 2 show the FFE of Al 6061-T6 and SS 304 for different type of loads and frequencies. The results demonstrate the validity of the constant cumulative entropy at failure for Al and SS specimens. Because the cumulative entropy for the final fracture of Al 6061-T6 was 4 MJm −3 K −1 , regardless of the test frequency, the thickness of the specimen, and the stress state.
For SS 304 specimens, the cumulative entropy was around 60 MJm −3 K −1 , also regardless of test frequency, the thickness of the specimen, and the stress state.   [86].This picture is ©2009 The Royal Society. Naderi and Khonsari [87] investigated the fatigue damage in metals based on the irreversible thermodynamic process. Their work showed that the cyclic plastic strain energy is the primary source of entropy generation in the low cycle fatigue test, and it reaches a critical value at the onset of fracture. They determined the degradation of different specimens subjected to low cyclic bending, tension-compression, and torsional fatigue based on entropy generation. The entropy generation equation is given by where is the total entropy generation at the onset of fracture, is the cyclic plastic energy, is the absolute temperature. The cyclic plastic energy determined by Morrow's cyclic plastic energy dissipation formula is given below where is the cyclic strain hardening exponent, is the fatigue ductility coefficient, denotes the fatigue strength coefficient, is the final number of cycles when failure occurs, is fatigue strength exponent, and is fatigue ductility exponent. Damage evolution is given by is the initial damage, is the critical value of entropy generation at the time when temperature starts to rise just after the steady-state phase, and is a curve fitting parameter.
Applying the same entropy generation equation, Amiri and Khonsari [88] investigated the degradation of metal fatigue. They showed that empirical fatigue models such as Miner's rule, Coffin-Manson equation, and Paris law could be applied to calculate the irreversible thermodynamic entropy production. In work by Naderi and Khonsari [89], they again assessed the degradation of machinery components associated with fatigue. They established a real-time fatigue damage monitoring methodology based on calculating the cumulative entropy production and is independent of the geometry, loads, and frequency. The developed fatigue monitoring unit automatically shuts down the machine before the FFE based on a given factor of safety. A laboratory test using Al 6061-T6 and SS 304 is provided. Their entropy calculation scheme is given by where ̇ represents the entropy production rate ( ̇≥ 0), the heat flux, the surface temperature, the cyclic plastic energy per unit volume, is the final number of cycles when failure occurs, and are curve fitting parameters, is fatigue strength exponent and is fatigue ductility exponent, is the fatigue ductility coefficient, denotes the fatigue strength coefficient. Similar work was also done by Teng et al. [90] for normalized SAE1045 steel.
Liakat and Khonsari [91] studied the high cycle fatigue of medium carbon steel 1018 using the concept of thermodynamic entropy generation in a degradation process. The evolution of cyclic plastic strain energy and temperature gradient in the specimen is proposed to be the mechanisms for entropy generation. They observed similar results as Naderi et al. [86] that the FFE is constant at failure. Therefore the concept of tallying entropy is valid for the fatigue life prediction evolution and damage evolution of a material subjected to cyclic fatigue. Their entropy generation equation is given by Ontiveros et al. [92] examined the fatigue crack initiation of AA7075-T651 by series of experimental tests. In his work, entropy generation is proposed to be an index of fatigue crack initiation. It is calculated based on the hysteresis energy and temperature increment of the material, as the Equation (60). He later addressed that fatigue life is a function of entropy and also the likelihood of microcrack formation. This finding is supported by the scanning electron microscope analysis performed on the fractured surfaces, which can capture the fatigue striations, and optical microscope analysis, which can capture the persistent slip bands.
Guo et al. [93] introduced a model based on intrinsic dissipation to predict metallic materials' high-cycle fatigue life. The intrinsic dissipation was defined as the unrecoverable microstructure motion that induced damage in their work. It was derived by incorporating two critical stress amplitudes related to the corresponding microstructure. A fatigue life prediction model was then established by taking the intrinsic dissipation as a fatigue damage indicator, and energy dissipation threshold as failure criterion. The quantitative assessment of cumulative damage caused by unrecoverable microstructure motion is given by The high cycle fatigue life prediction model is given by where is the cycle number, is the energy dissipation threshold for fatigue failure which can be identified using the proposed intrinsic dissipation model with experimental data, is the frequency, is the stress amplitude, is the critical stress amplitude for the onset of the unrecoverable deformation mechanism, which is considered as a material fatigue limit = in the proposed model, and are experimental curve fitting parameters. H(・) denotes the Heaviside step function.
Ribeiro et al. [94] estimated the FFE for the low cycle fatigue of the Al-2024 specimen. Various estimation methods, including empirical mechanical models and temperature measurements, were discussed. The first approach employed Park and Nelson's empirical model to calculate the entropy generation based on cyclic plastic work in the fatigue process. The second approach estimated the entropy generation in the fatigue process based on experimentally determined temperature by thermography. Results again indicated a constant FFE for specimens at failure. The FFE based on temperature is given by where is the total entropy generated during the fatigue tests obtained from thermal balance. In Equation (63), the first term is the contribution from heat conduction, the second term is from heat accumulation, third term is from convection and radiation when we consider them as heat source.
is the mean temperature, is the thermal conductivity, is a parameter obtained by parabolic curve fitting to the temperature profile, is the density, is the heat capacity, ℎ is the global heat transfer coefficient, is exchange surface, is the specimen volume. Roslinda Idris et al. [95] proposed an entropy-based model to predict the fatigue crack growth rate for dual-phase steel under spectrum loading. They monitored the crack length and temperature evolution during the fatigue crack test until failure to validate the estimation. Results showed that the proposed model capture the fatigue crack growth rate of the material accurately under spectrum loading.
The number of cycles for crack length to propagate from distance to , is developed as where is the number sequence, is the difference between the numbers of cycles for crack length interval, is the interval between the crack length and represents / . Osara and Bryant [96] developed a Degradation-Entropy-Generation (DEG) methodology based on the original work by Basaran and Yan [83] and Basaran and Nie [97] for system and process characterization and failure analysis on metal low-cycle fatigue. They related desired fatigue measures (stress, strain, cycles, or time to failure) to the loads, materials, and environmental conditions (including temperature and heat) through the irreversible entropies generated by the dissipative processes that degrade the fatigued material based on degradation-entropy generation theorem. The proposed formulation was validated with a steel shaft's fatigue data under bending and torsion. Their entropy generation equations are given by where is the load entropy and is called MicroStructuroThermal (MST) entropy. Indices 1, 2, 3, …, m correspond to times , , ,…., , ∆ is the total number of cycles within sampling time increment, is the cyclic strain hardening coefficient, , , , are elastic and plastic strain, stress, absolute temperature at , respectively, is the density, is the heat capacity, is isothermal loadability, is thermal expansion coefficient.
Sun et al. [98] measured the dependence among different degradation processes by proposing a copula entropy approach, which is a combination of the copula function and information entropy theory. The copula function was employed to identify the complex dependence structure of performance features, and information entropy theory was used to quantify the degree of dependence. In their work, the multiple degradation data of a microwave electronic assembly were used to verify the proposed approach. The copula entropy is given by where ( , , … , ) is the probability density function of the copula function; = ( ) = ( ≤ ), = 1,2, … , , represents the marginal distribution function of random variables. The probability density function, ( ), of the th performance feature degradation increment is calculated by where is the time interval; is the width of the time interval; ℎ is the width of the form smooth parameter, ∆ is the increment of degradation data, and () is a kernel function, which is a standard Gaussian distribution.
Yun and Modarres [99] used entropy as a damage metric for the fatigue damage on metallic materials. They calculated the entropy from three energy dissipation sources, including mechanical, thermal, and acoustic emission (AE), all collected experimentally. These experimentally obtained data are estimated, and their corresponding entropy generations were correlated with the observed fatigue damages.
The classical thermodynamic entropy equation is given by The AE information (Shannon) entropy is given by The entropy from the definition of statistical mechanics is given by where is the thermodynamic flux due to heat conduction, is the thermodynamic flux due to diffusion, is the chemical potential, is the mechanical stress, is the plastic strain, is the chemical reaction rate, is the chemical affinity, is the coupling constant, is the thermodynamic flux due to the external field, and is the potential of the external field. ( ) is a corresponding discrete histogram for processed digitized data, (+ ) and (− ) in the fatigue damage process are interpreted as the forward and reverse work distributions over many load cycles, respectively.
Young and Subbarayan [100] proposed the maximum entropy principle as a systematic approach to derive a cumulative distribution function to model damage. They performed low cycle fatigue tests on aluminum 2024-T351 to measure entropy from the cyclic plastic work, where Ramberg-Osgood constitutive model is employed to fit the hysteresis loops to calculate the dissipation per cycle and corresponding damage. The developed damage models predict fatigue life more accurately and consistently than models such as the Weibull distribution function and the Coffin-Manson relation. The Ramberg-Osgood plasticity model for stress-strain loops are where Δ is plastic strain range, ∆ is stress range, is Ramberg-Osgood strength parameter and 1 ⁄ is Ramberg-Osgood exponent. The damage per reversal is proposed as a function of inelastic dissipation per reversal with power law fit = ( And the inelastic dissipation for a monotonic test is where f (_) denotes a functional relationship with the argument, is the true fracture stress, is the true fracture strain, 2 is the total reversals to failure, is the total inelastic dissipation (per unit volume) to failure.
Based on Boltzmann's statistical entropy equation and continuum damage mechanics, Wang and Yao [101] proposed an entropy-based failure prediction model for the creep of metallic materials. The relation between entropy generation rate and normalized creep failure time was developed empirically and validated experimentally. Then the empirical entropy-based creep strain prediction model was established, as follows where is the creep strain, is the initial value of cumulative plasticity in the microstructure of material, which represents the value at the beginning of creep damage accumulation; is the threshold value of cumulative plastic variable when creep failure occurs. The corresponding creep time when = and = are and , respectively. The parameter is related to the applied stress, temperature, and material properties.
Sosnovskiy and Sherbakov [102] proposed the main principles of the physical discipline-mechano-thermodynamics (MTD), which is based on using entropy as a bridge between mechanics and thermodynamics. A unified mechano-thermodynamical function of limiting states (critical according to damageability) of polymers and metals is determined from more than 600 experimental results. These limiting states are also known as FFE states. The generalized expressions for entropy in the MTD system consists of a liquid (gas) medium of volume and a solid of volume , given by where is the specific tribo-fatigue entropy, is the limiting density of the internal energy treated as the initial activation energy of the disintegration process, is the total effective energy of the system, is a dimensionless parameter of local energy damageability, = . is the heat flux, is the chemical potential, is the number of mols per unit mass, are experimentally found coefficients.

Physics Based Evolution Functions: Unified Mechanics Theory
The models presented in Sections 2 and 3 that predict material degradation evolution are based on empirical curve fitting of the experimental data. Most do not satisfy the 2nd law of thermodynamics. Because according to the 2nd law of thermodynamics only entropy can be a damage or degradation criteria, not stress or strain or displacement. Moreover, only entropy has additive property. In Section 4.1, there are indeed many research works that use intrinsic dissipation or entropy as a damage metric based on the first and second law of thermodynamics. Still, they only consider entropy as a threshold to failure. For example, the term "fatigue fracture entropy" that repeatedly appeared in Section 4.1 was used as an indicator at fatigue failure. They did not derive a damage evolution equation based solely on the entropy evolution, which does not require an empirical evolution function. On the other hand, the physical based model introduced in this section eliminates the need for curve fitting a damage evolution function to a test data. Instead, the analytical thermodynamic fundamental equation, the second law of thermodynamics, and Boltzmann evolution are used.
Unified mechanics theory (UMT) [7,103] is a purely physics-based approach that doesn't need any experimental data curve fitting for degradation evolution function. It is obtained by combining the universal laws of motion of Newton and the first and second laws of thermodynamics directly at the ab-initio level. The second and third laws of unified mechanics are given by, [7] The second law The third law = where is the acting force, and is the stiffness and change in length of the reacting system, is the thermodynamic state index (TSI) which accounts for the dissipation. The meaning of thermodynamic index will be explained in the next paragraph.
In unified mechanics the material is treated as a thermodynamic system. As a result, governing partial differential equations of any system automatically include energy loss, entropy generation, and degradation of the system. A damage evolution is calculated along the Thermodynamic State Index axis, TSI, given by [7] where is the critical value of TSI, Δ is the change in entropy, is the molar mass and is the gas constant. Equation (78) is the normalized form of the second law of thermodynamics as defined by Boltzmann. When a material in ground (reference) state, it is assumed to be free of any possible defects, i.e., "damage", it can be assumed that "damage" in material is equal to zero. Therefore, TSI will be = 0. However, does not have to be zero initially. In final stage, material reaches a critical state, such that entropy is maximum (system reaches thermodynamic equilibrium, material is defined as failed). At this stage, entropy production rate will become zero. Therefore, TSI will be = 1.
Thermodynamic State Index is an additional linearly independent axis complementing the Newtonian x, y, z, time coordinates. UMT maps the entropy generation rate between zero and one on the TSI axis. It predicts the lifespan of any closed system without curve fitting an empirical model to a test data just using the thermodynamic fundamental equation of the material/system. Figure 3 shows the coordinate system in UMT. It is important to emphasize that in the new coordinate system derivative of displacement with respect to entropy is not zero because TSI is a linearly independent axis. The following example can explain the UMT coordinate system: "A 5-year-old boy with terminal illness and a 100-year-old sick person would have different coordinates in the time axis. But they have the same coordinate on TSI axis. Because in Newtonian space-time coordinates without incorporating thermodynamics, the physical state cannot be defined. On the other hand, in the unified mechanics theory their physical state is represented by the TSI axis in addition to the space-time coordinates. On TSI axis at Φ = 0.999, 5-year-old boy with terminal illness and a sick 100-year-old man will have the same thermodynamic state index coordinate. It is shown that they have about the same remaining life on TSI axis. Essentially, a person's age according to calendar does not give any information about that person's thermodynamic physical state". In UMT, in addition to nodal displacements, the entropy generation rate is also necessary to relate microstructural changes in the material with spatial and temporal coordinates. In the following, studies that adopted this concept and have been experimentally and mathematically validated will be listed and discussed.
UMT was used for fatigue life prediction under thermo-mechanical loading [83,97,[103][104][105][106][107][108][109][110][111][112][113][114][115][116][117]. Life prediction simulations were validated by experiments. Noushad et al. [103] used the UMT to predict the fatigue life of Ti-6Al-4V alloys based on fundamental equation of the material. Similarly, Egner et al. [118] analyzed the low cycle fatigue behavior of P91 steel using the exponential damage evolution equation derived in UMT. They performed experimental tests to obtain different aspects of material behavior and then developed a constitutive model within the framework of thermodynamics, including irreversible processes with internal state variables. Two different fatigue damage modeling approaches are applied to describe the final stage of material cyclic softening. In the classical ductile damage model (Chaboche-Lemaitre ductile damage model [62,119]), the process of micro-cracks and micro-voids development starts when the accumulated plastic strain reaches a certain threshold, while the entropy-based model includes damage evolution from the very beginning of the loading. Egner et al. [118] concluded that UMT gives better fatigue predictions, in spite of the fact that UMT has no curve fitting parameters.
The fundamental equation of the steel under low cycle fatigue is given by: and the exponential damage evolution equation is derived as follows: where is the stress tensor, is the plastic strain tensor, is the absolute temperature, is the unit mass, is thermal conductivity, is the distributed internal heat source per unit mass, is the TSI, Δ is the change in entropy, is the molar mass and is the gas constant. UMT was also used for fatigue under electrical-thermal-mechanical loading [120][121][122][123][124][125][126][127][128][129][130][131][132][133][134]. In the literature, a model was implemented into a finite element procedure to predict nanoelectronics solder joint's time to failure under high current density. The nonlinear viscoplastic time-dependent nature of the material and current crowding effects are considered in the formulation. Simulations were validated by experiments.
The fundamental equation of the nanoelectronics solder joint is given by: where is effective vacancy diffusivity, is the equilibrium vacancy concentration in the absence of a stress field, is instantaneous atomic vacancy concentration, is normalized vacancy concentration = ⁄ , * is vacancy effective charge number, is electronic charge of an electron, is Boltzmann's constant, is absolute temperature, is metal resistivity, ⃗ is current density vector, is atomic vacancy relaxation ratio, Ω is atomic volume, = trace( ) 3 ⁄ , is the plastic strain rate tensor, * is heat of transport, the isothermal heat transmitted by the moving atom in the process of jumping a lattice site. Δ is the entropy production, is Avogadro's constant. Besides metals, UMT has also been used for fatigue life prediction in particle filled composites. These models were verified experimentally [97,[135][136][137][138][139]. In short, the models based on unified mechanics theory all used the same damage evolution equation derived from Boltzmann equation shown in Equation (78) to predict the material failure without curve fitting. However, the fundamental equation (entropy generation equation) differs according to the material and the micro-mechanisms activated. Deriving these fundamental equations is the essential key in using unified mechanics approach.

Mechanical Treatment
The models discussed in previous sections did not take the effect of manufacturing process or the surface finishing of the specimen into consideration. However, it is well understood that manufacturing processes have effect on the residual stress, surface roughness, and surface integrities of the metallic material that can be correlated to fatigue performance. In this section, some research work about the effects of mechanical treatments that improve fatigue life are introduced.

Manufacturing Process
Manufacturing processes are known to drastically impact the performance and lifetime of a component subjected to fatigue in service. Therefore, it is important to understand the effect of manufacturing processes on surface integrity. Xu et al. [140] investigated the influence of shot peening on the fatigue response of Ti-6Al-4V surfaces subject to different machining processes. They showed that the milled, abrasive waterjet (AWJ) and wire electrical discharge machined (WEDM) surfaces has better fatigue strength than the corresponding surface ground and large area electron beam melting (LAEBM) specimens. On the other hand, the fatigue strength of shot-peened specimens was in general superior to those in the as-processed state, which is attributed the induced compressive residual stresses on the test face under maximum tension during testing regardless of the base surface processing method.
Choi [141] studied the Influence of feed rate on surface integrity and fatigue performance of hard machined surfaces. The fatigue life was predicted based on the residual stress and micro-hardness distribution to investigate the influence of the feed rate on fatigue performance. They showed that a higher feed rate increases crack initiation life and crack propagation life. Higher feed rate also induces more compressive residual stresses. Models of crack initiation life, crack propagation life, and rolling contact fatigue life are given as follows: Crack initiation life model Crack propagation life model Rolling contact fatigue life model where ∆ is the range of maximum shear stress during a cycle at the initiation point, and is the frictional stress, and are the Knoop Hardness number at the bulk material and the local Knoop Hardness number, respectively, is the material constant, ∆ is the stress intensity factor range, is an empirical index, and is the half length of the crack.
Choi [142] also studied the influence of the rake angle on the surface integrity and fatigue performance of hard machined surfaces. Using the models shown in Equations (82)- (84), the crack initiation life, crack propagation life, and rolling contact fatigue life of specimens machined at various rake angles can be estimated. The results showed that 1. A higher rake angle induces more compressive residual stresses in the circumferential and radial directions, and a lower rake angle induces a more softened layer. 2. A higher rake angle increases the crack initiation life of the specimen, especially evident in lower loading conditions. Same for the crack propagation life. 3. The rolling contact fatigue tests confirmed that the rake angle influences the fatigue life significantly and that the effect is further increased if the loading is reduced.
The Electrical Discharge Machining (EDM) process is a metal fabrication process used to machine any electrically conductive material. It is widely used for machining hard metals or manufacturing parts with complex shapes, such as aeronautical industry components. Still, the inevitable thermal input caused by the EDM process that leads to metallurgical change in the material impedes the EDM application in selected aeronautical components. Therefore, the effect of EDM process on fatigue service life becomes an important consideration. Ayesta et al. [143] investigated the influence of the wire electrical discharge machining (WEDM) process on the fatigue behavior of the nickel-based Inconel ® 718. From the fatigue test results conducted under a load ratio R = 0.1, they presented the relation coefficient between the fatigue strengths of WEDM and ground specimens as = • = 1.61 . 10 ≤ ≤ 10 (85) where is the fatigue strength of the WEDM specimens, is the fatigue strength of the ground specimens, is the number of cycles. They showed that conducting the WEDM process with the new generation of machines has a detrimental effect at high fatigue cycles, reducing fatigue strength by approximately 10% in comparison with ground specimens, as calculated in Equation (85). However, at low fatigue cycles, no significant differences were observed between ground and WEDM samples.
Additive manufacturing is a novel manufacturing method that gained popularity in the early 2010s and brought increasing applications in recent years. However, the metallic parts after the additive manufacturing process show a high surface roughness and an irreproducible and inhomogeneous surface morphology that limit its application in highly stressed and cyclically loaded areas. Bagehorn et al. [144] investigated the effects of different mechanical surface finishing processes on the surface morphology and subsequent fatigue performance of additive manufactured Ti-6Al-4V specimens, including milling, abrasive blasting, vibratory grinding, and micro machining. They showed that Ti-6Al-4V is sensitive to notches, hence the built surface roughness strongly affects its fatigue performance. All surface finishing processes achieved a substantial roughness improvement and therefore also a superior fatigue behavior (increased fatigue life). However, they did not find a direct correlation between a specific roughness value and fatigue life.
The non-conventional machined processes such as abrasive water Jet (AWJ), wire electrical discharge machining (WEDM), and ultrasound vibration assisted milling (UVAM) mentioned earlier have received more attention in recent years. Suárez et al. [145] presented the differences in the surface integrity and fatigue performance of non-conventional machined surfaces to conventional milling for Alloy 718 parts. These machining processes have great influence on the residual stresses and surface roughness. They showed that both residual stresses and surface roughness are correlated with fatigue strength, and that low surface roughness typically results in high surface fatigue strength. However for roughness values in the range 2.5-5 μm Ra, it is primarily dependent on workpiece residual stress and surface microstructure, rather than roughness. Their experimental results demonstrated that among these machining processes, WEDM samples have the lowest fatigue strength while UVAM sample has the highest fatigue strength (also improved fatigue life). Besides, UVAM results showed an improvement on the surface integrity of the final workpiece while AWJ and WEDM showed poorer results.
Holmberg et al. [146] also performed surface integrity investigations for prediction of fatigue properties after different machining (conventional milling and non-conventional machining) of alloy 718. Besides on the study of topology, residual stress, surface roughness, and fatigue life for non-conventional machining, they proposed a statistical multivariate regression model based on the results of fatigue testing and surface integrity investigations to predict the fatigue life, as follows: are regression constants, 10 is the ten-point height parameter for the surface roughness, is the residual stress at 10 , is the full width at half maximum at 10 μm (10 μm below the surface). To secure a robust prediction result using Equation (86), multiple observations per method is necessary.

Surface Finishing
As fatigue cracks initiate predominantly at the free surface of a material, the condition of the surface can be assumed to be critical with regards to fatigue strength. Obiukwu et al. [147] performed a series of low cycle fatigue tests on low and medium carbon steel to investigate the effect of surface finish on fatigue life. They concluded that polished surface increases the fatigue life of machined components when compared to ground, lathe turned and filled specimens because rough surfaces form stress concentration centers that decrease endurance limit. However, the fatigue strength was impaired in low-cycle regime (<10 4 ) of low carbon steels since hardness is a more important factor.
Avilés et al. [148] quantified the improvement in the high-cycle fatigue strength of AISI 1045 normalized steel after low-plasticity ball burnishing. They showed that the fatigue strength of the ball-burnished specimens improved for 3.25 × 10 6 cycles, and the bending fatigue limit is increased by 21.25% compared to non-treated specimens. The variation of the coefficients that relate the non-treated specimens to low plasticity burnished (LPB) specimens is presented in an analytical expression. For their test specimen, the value of the LPB coefficient is found by where is the fatigue strength of the ball-burnished specimens at N cycles, is the fatigue strength at N cycles, is the bending fatigue limit of the non-treated specimen. Besides the fatigue strength, the law of relaxation of the effective mean stress, the equation of residual stress relaxation due to LPB are also shown in their work. The modified coefficient shown in Equation (87) considers following four effects caused by LPB: Enhanced surface roughness, Reduced grain size in the surface zone, Increased surface biaxial compressive residual stresses, and Increased HBN surface hardness. In addition, the effect of ball burnishing process on the roughness, hardness, and residual stress can also be found in Rodriguez et al. [149], in which a G10380 steel specimens milled and treated with a ball burnishing process assisted by vibrations were analyzed.
Laser polishing (LP) is a surface finishing method that can reduce the surface roughness of machined parts by laser beam. When performing laser polishing in the air, the pores, inclusions, and oxides in the surface layer result in a detrimental effect on fatigue strength. These unwanted surface features can be reduced by performing the laser polishing in a N2 atmosphere (LPN). Avilés et al. [150] studied the high-cycle fatigue behavior of normalized 1045 AISI steel after laser polishing in a N2 atmosphere. The fatigue test results showed that LP improves the HCF fatigue strength of the base material if its average surface roughness of the base material before LP is larger than 0.35 μm. For lower number of cycles, a small reduction of fatigue strength in LP specimen is assumed to be attributed to a combination of the microstructural changes in the heat-affected zone (HAZ), surface features of LPN specimens, and the relaxation of the moderate compressive residual stresses produced by LP. Quantitative data of the modifying factors of the fatigue limit, fatigue strength for LP and LPN due to surface roughness at different cycles, and the calculation procedures of the fatigue strength and fatigue limit for LP and LPN at different cycles and loading ratios, are shown in their work. The modifying coefficient value (ratio between the fatigue limit of untreated mirror finish specimens and treated specimens) attributed to LPN show significant improvement compared to those coefficient attributed to LP in air Avilés et al. [151] performed a comparative study on the effect of shot-peening (SP) and low-plasticity burnishing (LPB) on the high cycle fatigue of a DIN 34CrNiMo6 alloy steel. The experimental data showed that compared with the machined specimens, the fatigue limit of the SP specimens and LPB specimens increased by 39% and 52%, respectively, as follows = ≈ 1.39 = ≈ 1.52 (88) where and is the improvement index for SP and LPB specimens, respectively. , , is the rotating-bending fatigue limit of the machined specimens, SP specimens and LPB specimens, respectively. The improvement index and modification coefficient for SP and LPB specimens compared to mirror-polished specimens are also established based on curve fitting function using number of cycles as an input. Another focus in their study on the relaxation of the surface residual stresses under cyclic load. They developed a model to predict the evolution of residual stresses with time of SP and LPB specimens.
Selective laser melting (SLM) is a manufacturing technique that melts metal powder particles using heat source to form complex components layer by layer. Yu et al. [152] evaluate the effects of different surface finishing processes, hot isostatic pressing (HIP) and heat treatments on the fatigue performance of selective laser melted (SLMed) Ti-6Al-4V. They found that surface finishing processes of turning, grinding, grinding followed by sandblasting, and polishing can effectively reduce the surface roughness of as-SLMed surfaces (Rz = 68.66 μm) and significantly improve the fatigue performances. On the other hand, the HIP can reduce or eliminate the pores and modify the microstructure. Thus significantly extending the period of fatigue crack initiation, leading to an increased fatigue limit More surface finishing related research work was published recently. Rodriguez et al. [153] presented a theoretical and experimental validation using different burnishing conditions. They showed that roller burnishing technique greatly improves surface roughness, and eliminates the kinematic-driven roughness pattern of turning, thus achieving isotropic surface topography (isotropic finishing) on cylindrical components made of austempered ductile iron (ADI). A detailed comparison between ball burnishing and roller burnishing were presented in their work. Tan et al. [154] investigated the evolution of surface integrity and fatigue properties after milling, polishing, and shot peening of TC17 alloy blades. Kahlin et al. [155] performed a comparative study on the improved fatigue strength of additively manufactured Ti6Al4V by 5 different surface post processing-including shot peening, laser shock peening, centrifugal finishing, laser polishing, and linishing.
All these works identified the importance of mechanical treatment on fatigue performance.

Conclusions
This review paper discussed some widely used damage, void evolution, and fatigue life prediction models. For example, empirical curve fitting models including GTN model, Johnson-Cook model, micro-plasticity models, phase-field models; Entropy based models including models using entropy as a damage metric with an empirical evolution function, and pure physics-based (unified mechanics theory based) models.
Among all the empirical curve fitting models, the GTN model is the most popular. It is based on void nucleation, growth, and coalescence in materials. It is suitable for simulating the fracture in ductile materials. However, it is inaccurate in low stress triaxiality conditions, shear dominated mechanisms, and high strain rate loading. The modified GTN models improved some of the mentioned drawbacks by proposing various novel void growth and nucleation evolution equations, or simply adding new damage variables (shear damage). Still, the micro-mechanisms related curve fitting constants are difficult to determine, and some values even need the help of stochastic methods, which limit its industrial applications.
The Johnson-Cook model, on the other hand, is a plasticity model that is capable of describing the material response (equivalent plastic stress, failure equivalent plastic strain) under high strain rate considering the effects of strain hardening, strain-rate (viscosity) and thermal softening, separately. The calculated equivalent plastic strain is used to develop a dynamic failure model called Johnson-Cook dynamic failure model. However, the model may not provide accurate results when temperature and strain rate dependent hardening is considered. It is also incapable of capturing the static and dynamic recovery, and reflecting the effects of load path and strain-rate history in large deformation processes. The modified J-C models fixed the drawbacks by improving the functional relationships of strain effects, strain rate effects, and temperature effects in various empirical ways. Compared to the GTN type model, Johnson-Cook model's empirical constants are obtained purely from curve fitting a flow stress test data without strong mechanics-based theory. The simplicity of determining empirical constants makes it widely accepted in engineering applications.
The microplasticity model includes crystal plasticity and some other models based on microscopic dislocation motions. Crystal plasticity is an inherently multi-scale process that is established based on the dislocation slip mechanism to construct constitutive laws for analyzing responses in polycrystalline materials, such as the texture evolution, inhomogeneous deformation, damage fracture, fatigue crack nucleation, crack path, fatigue life prediction, creep, and recrystallization. It provides accurate results but require careful calibrating of the material parameters and defining crystallographic directions and orientations. The simulation is usually limited to very small scales.
Besides above mentioned models, Zerilli-Armstrong model has shown superior performance compared to the Johnson-Cook model when facing ultra-high strain rate deformation; Lemaitre damage model can also be used to simulate ductile fracture of materials; phase-field models are widely used to simulate the dynamic crack propagation and ductile fracture in metals due to the advantage of eliminating the need of explicit fracture criterion. For non-metallic composites, Hashin's failure criterion and energy based ductile criterion are usually utilized for fiber and matrix, respectively.
All of the above models provide useful tools for simulations. However, according to the 2nd law of thermodynamics only entropy can be a damage or degradation criteria, not stress or strain or displacement or force. The cumulative entropy production was found to be a material property independent of the specimen geometry, loading type, and frequencies. Because for the same material the cumulative entropy production is constant at failure, it is widely used as a damage metric. These models using entropy as a damage metric are classified as entropy based models. For the entropy based models, the most difficult part is to establish the entropy generation equation (the fundamental equation). The damage evolution equation can be either empirical or purely physics-based. The unified mechanics theory is a purely physics based approach that eliminates the need for curve fitting damage evolution function to a test data. Instead, it utilizes the analytical thermodynamic fundamental equation (entropy generation equation), the second law of thermodynamics, and Boltzmann evolution to derive the damage evolution equation without curve fitting. The entropy based models have been verified against test data extensively.
Finally, the manufacturing techniques and surface treatment are known for affecting fatigue performance. These effects are neglected in previous discussions. To give a good overview for readers, some recent research works about the effect of metal treatment on fatigue performance and surface integrity are discussed.