Damage Evolution Due to Extremely Low-Cycle Fatigue for Inconel 718 Alloy

: This paper evaluates the damage evolution process under extremely low-cycle fatigue (ELCF). The study explores the damage behavior under di ﬀ erent stress states. The inﬂuence of the multiaxial state of stress on the metal’s life is determined. Two di ﬀ erent stress states were examined: (a) axisymmetric and (b) plane-strain. The study is based on the modiﬁed Mohr–Coulomb (MMC) ductile fracture criterion that was extended to cover the ELCF regime in a previous research study. Four distinctive geometries are designed to study the e ﬀ ect of di ﬀ erent stress states on ELCF life and damage evolution. The damage model is calibrated for life prediction to agree with the ELCF experimental results. The investigation of the damage evolution behavior is dependent on equivalent plastic strain, stress triaxiality, Lode angle, and cyclic loading e ﬀ ect. The damage evolution is extracted from Abaqus ﬁnite element simulations and plotted versus the equivalent plastic strain. The damage accumulation shows nonlinear evolution behavior under cyclic loading conditions. SEM images were taken to further study the microscopic failure mechanisms of ELCF.


Introduction
Extreme dynamic and fatigue loadings are very common in engineering product designs. For example, during the shutdown and startup operations, engines and machine components are highly susceptible to this type of loading. It commonly causes engines and machine components to fail due to very high strain cyclic loading, referred to as extremely low-cycle fatigue (ELCF). The dynamic cyclic loadings are due to direct fluctuating forces exerted among the machine parts or due to massive vibration within the machine parts. ELCF is also seen in many engineering structures that are exposed to aerodynamic forces. Therefore, predicting, understanding, and evaluating the damage evolution caused by ELCF with a suitable mathematical model is crucial and challenging. The parameters in such a model should consider the complex geometries of the machine parts, the applied loading conditions, and the material properties.
An extended damage accumulation model for ELCF life prediction was developed with the aid of experimental data in Reference [1]. The original modified Mohr-Coulomb (MMC) model was introduced in Reference [2]. Years later, the MMC model was extended to predict ELCF in Reference [1]. Fatigue damage refers to the process of deterioration for engineering parts and components. The damage model parameters are also strongly related to the applied loads [3].

Plasticity Model with a Combined Hardening Rule
The plasticity behavior of fully reversed push-pull cyclic loading conditions for Inconel 718 involves a combined material hardening rule. The kinematic hardening model used in this research was proposed by Chaboche in References [17,18] that captures the Bauschinger effect. The Bauschinger effect is caused due to the microstructure deformation by the irreversible dislocation movement along the crystallographic planes that changes the material's stress-strain characteristics and displaces the yield surface. The model consists of a number of nonlinear back-stresses as in Equation (1). The Chaboche model is a superposition of three Armstrong and Frederick kinematic hardening rules [19]. The three nonlinear back-stress components (α = α 1 + α 2 + α 3 ) were calibrated by simulating the stable Inconel 718 hysteresis loop of the force-displacement curve of the smooth round bar specimen (denoted as "R0") as in Figure 1. The parameters in Equations (1) and (2) are as follows: α i is the backstress tensor, C i and B i are material-dependent constants, and . ε pl is the accumulated equivalent plastic strain rate as described in Equation (2). The isotropic hardening rule used was the Vocé isotropic strain hardening function in Equation (4). It is a nonlinear isotropic model that describes a yield surface expansion.
. ε pl = ( The parameters in Equation (6) (A, N, c , c  , c , and c ) are material parameters that are calibrated based on a scripted MATLAB code to match the simulation crack initiation to the experimental one. The first function c D in Equation (7) has the parameter that controls the damage evolution behavior as concave or convex (see Figure 3). The second function 1 c D μ incorporates the effect of the different tensor direction between the current stress and the back-stress. This function is activated when reverse loading (non-proportional) exists. The material parameters in both functions (c , c , β , β , and k) is calibrated under reverse loading tests. The damage model was initially set to predict ductile fracture damage (in Reference [2]) and was extended later in Reference [1] to predict ELCF. The damage model is based on the parameters Ɛ , ,  , and μ defined as equivalent plastic strain, stress triaxiality, Lode angle, and effect of cycling loading during ELCF, respectively. The state parameters change as the loading process continues and they can be tracked by FEA (using an Abaqus/Explicit user-defined material subroutine). The importance of each parameter in understanding the damage evolution, due to ELCF, is discussed in the following sections. All damage material parameters in Equations (6) and (7) are set in Tables 2  and 3.    Table 3. Material of non-proportional parameters used in Equation (7). −6.01 3.0 0.049 0.000001 − 1.0

Equivalent Plastic Strain
The first state variable of the four variables in the extended MMC damage model is strain. It is observable, internal, and it is split into elastic strain Ɛ and plastic strain Ɛ . The experimental observations of the plastic strain show that the accumulated plastic strain mainly governs the material damage in the form of plane slip along the crystallographic planes and dislocation movements. The accumulated equivalent plastic strain used in this model is defined in agreement with the Von-Mises criterion, Ɛ pl 2 3 Ɛ : Ɛ . Studies show that there is a strong relation between the accumulated plastic deformation, which is represented by the equivalent plastic strain, and the ductile fracture    Table 3. Material of non-proportional parameters used in Equation (7). −6.01 3.0 0.049 0.000001 − 1.0

Equivalent Plastic Strain
The first state variable of the four variables in the extended MMC damage model is strain. It is observable, internal, and it is split into elastic strain Ɛ and plastic strain Ɛ . The experimental observations of the plastic strain show that the accumulated plastic strain mainly governs the material damage in the form of plane slip along the crystallographic planes and dislocation movements. The accumulated equivalent plastic strain used in this model is defined in agreement with the Von-Mises criterion, Ɛ pl 2 3 Ɛ : Ɛ . Studies show that there is a strong relation between the accumulated plastic deformation, which is represented by the equivalent plastic strain, and the ductile fracture  Table 2. Material parameters used in the modified Mohr-Coulomb (MMC) criterion in Equation (6).  Table 3. Material of non-proportional parameters used in Equation (7).

Equivalent Plastic Strain
The first state variable of the four variables in the extended MMC damage model is strain. It is observable, internal, and it is split into elastic strain ε el and plastic strain ε pl . The experimental observations of the plastic strain show that the accumulated plastic strain mainly governs the material damage in the form of plane slip along the crystallographic planes and dislocation movements. The accumulated equivalent plastic strain used in this model is defined in agreement with the . ε pl ij : . ε pl ij . Studies show that there is a strong relation between the accumulated plastic deformation, which is represented by the equivalent plastic strain, and the ductile fracture crack formation [27]. The current study investigates the relationship between the accumulated equivalent plastic strain and the damage accumulation during the cyclic loading process.

Stress Triaxiality Parameter
The dimensionless hydrostatic pressure parameter, termed stress triaxiality (η = σ H /σ), is the hydrostatic stress (σ H = (σ 1 + σ 2 + σ 3 )/3) over the Von-Mises equivalent stress (σ = 3J 2 ). The initial stress triaxiality can be estimated using Bridgman's analysis for an axisymmetric stress state as in Reference [28] and by the analytical solutions provided in Reference [29] for the plane-strain stress state. High-stress triaxiality zones (η > 1.0) can be seen in highly constrained regions such as the center of the notched specimen, or at crack tips. On the other hand, low-stress triaxiality zones (η < 0) are often observed where the principal shear stress is considerably higher relative to the hydrostatic stress (i.e., at the surfaces and protruding corners) [27]. Stress triaxiality has been extensively applied in many ductile fracture models [30][31][32][33][34]. Studies show that material ductility depends evidently on the stress triaxiality. The equivalent strain to fracture for numerous ductile materials, over different ranges of strain rates, is also strongly dependent on the stress triaxiality [23,35]. However, quantitative studies of stress triaxiality effects on fatigue life and damage accumulation have not been sufficiently investigated [1,36], yet this area merits extensive consideration. This study focuses on how the damage accumulation due to ELCF is affected by stress triaxiality as part of the fatigue life-prediction model (Equations (6) and (7)).

Lode Angle Parameter
The Lode angle θ, which is related to the third deviatoric stress invariant, ranges between 0 < θ < π 3 . The Lode angle can be normalized by the following equation θ = 1 − (6θ/π). Hence, the normalized Lode angle ranges between −1 < θ < 1. It is used to describe the stress state along with the stress triaxiality. For example, uniaxial tension corresponds to θ = 1 and η = 1 3 whereas pure shear corresponds to θ = 0 and η = 0, and uniaxial compression gives θ = −1 and η = −1 3 . Any state of θ other than the above states is a representation of a combined multiaxial stress state (Figure 2d).

Cyclic Loading Parameter Effect
During cyclic (push-pull) loading, the material experiences a shift in its yield surface. The yield surface moves slightly due to kinematic hardening during each cyclic loading. The center of the yield surface becomes the back-stress coordinate and not the origin. A mathematical model is set to capture the effect of the change in the non-proportional loading direction by measuring the current stress and the back-stress tensors, which is represented by a scalar parameter χ in Equation (8).
The range of χ is [0, 2] where χ = 0 means that there is no difference in the direction of the current stress and the back-stress tensors, whereas χ = 2 represents a 180 • difference in direction ( Figure 4). Another scalar parameter µ in Equation (9) is the accumulated differences captured by χ during the cyclic loading process. On the other hand, the parameter µ is incorporated in the damage accumulation model Equation (7). To illustrate this, one can see that the parameter χ changes from zero to two in the beginning of a reverse loading, and it reduces to zero after a certain amount of continuing loading. The parameter µ keeps increasing if there are changes in the loading path as in Figure 5. The change in loading condition during reverse loading from tension to compression is represented by (i) and (iii) whereas (ii) represents the change from compression to tension in Figure 5.

Damage Evolution Parameter
The damage evolution parameter c in Equation (7) is considered as the critical parameter in simulating the damage evolution under ELCF. The different curves in Table 4 is for the monotonic loading case. Previously, the material damage evolution due to ELCF has not been studied intensively in terms of knowing the development behavior of the damage accumulation during the cycling loading. The importance of the damage parameter expresses how to simulate the evolution of the material damage and whether it accumulates linearly or nonlinearly. Later in this research, the investigation of the ELCF damage is explained based on this parameter to comprehend the multiaxial ELCF damage evolution thoroughly.

Experiments
The specimens' geometries were designed in a way to confirm that the fatigue occurs at a selected stress state. The calculations to achieve a desired stress states were based on Reference [23] and the computer simulations confirmed the attained stress state by the following two parameters: stress triaxiality and Lode angle (Table 4). For instance, the uniaxial tension stress state had  = 0.333 and  = 1 by calculations while the plane strain stress state had  = 0.577 and  = 0. The metal is Inconel 718 (a nickel-based super alloy) and is made in USA (Metalmen Inc., Long Island City, NY, USA) with the chemical composition in wt% shown in Table 5. This material is widely known in the aerospace and defense industry. It is heavily used in most hot environment applications due to its excellent high strength and fatigue properties at high temperatures. ELCF test

Damage Evolution Parameter
The damage evolution parameter c in Equation (7) is considered as the critical parameter in simulating the damage evolution under ELCF. The different curves in Table 4 is for the monotonic loading case. Previously, the material damage evolution due to ELCF has not been studied intensively in terms of knowing the development behavior of the damage accumulation during the cycling loading. The importance of the damage parameter expresses how to simulate the evolution of the material damage and whether it accumulates linearly or nonlinearly. Later in this research, the investigation of the ELCF damage is explained based on this parameter to comprehend the multiaxial ELCF damage evolution thoroughly.

Experiments
The specimens' geometries were designed in a way to confirm that the fatigue occurs at a selected stress state. The calculations to achieve a desired stress states were based on Reference [23] and the computer simulations confirmed the attained stress state by the following two parameters: stress triaxiality and Lode angle (Table 4). For instance, the uniaxial tension stress state had  = 0.333 and  = 1 by calculations while the plane strain stress state had  = 0.577 and  = 0. The metal is Inconel 718 (a nickel-based super alloy) and is made in USA (Metalmen Inc., Long Island City, NY, USA) with the chemical composition in wt% shown in Table 5. This material is widely known in the aerospace and defense industry. It is heavily used in most hot environment applications due to its excellent high strength and fatigue properties at high temperatures. ELCF test Figure 5. An example of the evolution of the two non-proportional parameters χ and µ during cyclic loading. This is an experimental example of the "R1" specimen under extremely low-cycle fatigue (ELCF) with N f = 4.

Damage Evolution Parameter
The damage evolution parameter c g in Equation (7) is considered as the critical parameter in simulating the damage evolution under ELCF. The different curves in Table 4 is for the monotonic loading case. Previously, the material damage evolution due to ELCF has not been studied intensively in terms of knowing the development behavior of the damage accumulation during the cycling loading. The importance of the damage parameter expresses how to simulate the evolution of the material damage and whether it accumulates linearly or nonlinearly. Later in this research, the investigation of the ELCF damage is explained based on this parameter to comprehend the multiaxial ELCF damage evolution thoroughly.

Experiments
The specimens' geometries were designed in a way to confirm that the fatigue occurs at a selected stress state. The calculations to achieve a desired stress states were based on Reference [23] and the computer simulations confirmed the attained stress state by the following two parameters: stress triaxiality and Lode angle (Table 4). For instance, the uniaxial tension stress state had η = 0.333 and θ = 1 by calculations while the plane strain stress state had η = 0.577 and θ = 0. The metal is Inconel 718 (a nickel-based super alloy) and is made in USA (Metalmen Inc., Long Island City, NY, USA) with the chemical composition in wt% shown in Table 5. This material is widely known in the aerospace and defense industry. It is heavily used in most hot environment applications due to its excellent high strength and fatigue properties at high temperatures. ELCF test series of eight specimens under different high strain amplitudes were executed until specimen full fracture. The cyclic loading hysteresis loops and the strain-life curves are shown in Figures 6 and 7, respectively. The tests strain rate was 0.003/s and the elongation range related to reference gauge length are summarized in Table 6.
The servohydraulic testing machine (model 370.10) has a ±100 kN load cell range manufactured by MTS systems corporations ® in Eden Prairie, MN, USA. The force and displacement tests data were concurrently recorded by the loading machine and an optical measurement system, respectively. The optical measurement system measured and recorded the strains during the cyclic loading using a Digital Imaging Correlation (DIC) of type VIC-2D version 5 software made by correlated Solutions Inc ® in Irmo, SC, USA.  series of eight specimens under different high strain amplitudes were executed until specimen full fracture. The cyclic loading hysteresis loops and the strain-life curves are shown in Figures 6 and 7, respectively. The tests strain rate was 0.003/s and the elongation range related to reference gauge length are summarized in Table 6. The servohydraulic testing machine (model 370.10) has a ±100 kN load cell range manufactured by MTS systems corporations ® in Eden Prairie, MN, USA. The force and displacement tests data were concurrently recorded by the loading machine and an optical measurement system, respectively. The optical measurement system measured and recorded the strains during the cyclic loading using a Digital Imaging Correlation (DIC) of type VIC-2D version 5 software made by correlated Solutions Inc ® in Irmo, SC, USA.

Test Specimen Geometry
Four distinctive shapes were designed in order to study the effect of different stress states on ELCF life and damage evolution. The four shapes were a smooth round bar, small external notch round bar "R0", large external notch round bar "R1", and flat plane-strain bar "PE". A 3D demonstration and detailed dimensions in mm are expressed in Figure 8a,b and Table 7. All specimens were fabricated from one material block to ensure property and microstructure similarity.   Table 7. Reproduced from Reference [1], with permission from Elsevier, 2019.

Loading Conditions
The testing machine's upper grip held the top part of the specimen and remained fixed while the lower grip pulled/pushed the specimen. The fatigue tests were conducted at room temperature. The four shapes of the specimens (Table 7) were designed in such a way to attain different stress states, three of which (R0, R1, and R2) are classified as axisymmetric stress states and one (PE) is

Test Specimen Geometry
Four distinctive shapes were designed in order to study the effect of different stress states on ELCF life and damage evolution. The four shapes were a smooth round bar, small external notch round bar "R0", large external notch round bar "R1", and flat plane-strain bar "PE". A 3D demonstration and detailed dimensions in mm are expressed in Figure 8a,b and Table 7. All specimens were fabricated from one material block to ensure property and microstructure similarity.

Test Specimen Geometry
Four distinctive shapes were designed in order to study the effect of different stress states on ELCF life and damage evolution. The four shapes were a smooth round bar, small external notch round bar "R0", large external notch round bar "R1", and flat plane-strain bar "PE". A 3D demonstration and detailed dimensions in mm are expressed in Figure 8a,b and Table 7. All specimens were fabricated from one material block to ensure property and microstructure similarity.   Table 7. Reproduced from Reference [1], with permission from Elsevier, 2019.

Loading Conditions
The testing machine's upper grip held the top part of the specimen and remained fixed while the lower grip pulled/pushed the specimen. The fatigue tests were conducted at room temperature. The four shapes of the specimens (Table 7) were designed in such a way to attain different stress states, three of which (R0, R1, and R2) are classified as axisymmetric stress states and one (PE) is  Table 7. Reproduced from Reference [1], with permission from Elsevier, 2019.

Loading Conditions
The testing machine's upper grip held the top part of the specimen and remained fixed while the lower grip pulled/pushed the specimen. The fatigue tests were conducted at room temperature. The four shapes of the specimens (Table 7) were designed in such a way to attain different stress states, three of which (R0, R1, and R2) are classified as axisymmetric stress states and one (PE) is classified as a plane-strain stress state. The axisymmetric stress state had two identical principal stresses (σ 2 = σ 3 ) and the Lode angle parameter equals unity (θ = 1). The smooth round bar, R0, (η = 1/3) was used for calibrating the plasticity model (Reference [1]) and the two parameters (A, N) in Equation (6). The sharp notch "R1" and large notch specimens "R2" (η = 0.74 for "R1" and η = 0.51 for "R2") were used to study the effect of the stress triaxiality on the ELCF life, which is clearly demonstrated in the strain-life curve (Figure 7). The plane-strain stress state corresponds to θ = 0, and it adds additional important details to be further investigated. By comparing the plane-strain stress state with the axisymmetric stress state, the value of θ in specimens is different. Consequently, the effect of the Lode angle (or the normalized third stress invariant) on ELCF life can be studied.

Damage Evolution Evaluation
In this section, the material damage model for ELCF life prediction is evaluated covering a wide spectrum of fatigue loading conditions (i.e., variable and constant strain amplitude loading, different loading stress states, wide range of different strain ratio, and a wide array of alternating stresses). This damage model shows a suitable explanation of fatigue damage under all the aforementioned complexities in loading conditions and stress states. A quantitative study is described here to expose the complicated underlying ELCF mechanisms throughout the loading process.

ELCF Damage Evolution
Investigating multiaxial damage evolution due to ELCF can be performed by inspecting three distinguished features on the specimen's fracture surface: the crack initiation site, the crack growth behavior, and the final fractured surface. Crack initiation likely occurs in highly stressed regions (i.e., microvoids or discontinuities due to a crystallographic plane slip). Crack growth studies the crack propagation rate and crack type: transgranular or intergranular. The specimens' fractured surface reveals the fracture mode and the grains deformation using SEM. This section will focus on the microstructure damage and damage evolution prior to crack initiation. The ELCF process in the very early loading cycles creates microstructure changes in the bulk of the material and rapidly form strain localization because of the high strain amplitude. The mechanism of the microstructure deformation is caused by the irreversible dislocation movement along the crystallographic planes. As the dislocation pattern in many bundles (veins) becomes localized and the number of persistent slip bands (PSB) increases, the strain localization occurs and forms plastic strain accordingly. A consensus of opinion among researchers is that the PSB are where most cracks initiate in isotropic materials under reverse loading [37]. This explanation of the irreversible damage mechanism is from a microscopic perspective. However, in numerical simulations, the irreversible damage mechanism is quantitatively represented as shown by the accumulated damage D in Equations (6) and (7). The tendency of the damage evolution behavior during the cycling loading is ambiguous and yet to be discovered. It is not known whether the damage under ELCF accumulates linearly or nonlinearly. In other words, is the damage accumulation rate steady during ELCF or not?
To study the damage accumulation, the parameter c g in Equation (7) is the focal parameter that will examine the damage evolution behavior under ELCF. For this reason, the damage model (Equations (6) and (7)) was set in a MATLAB code and graphed versus ε pl with a wide range of different c g values, −8 < c g < 8 (Figure 3). The various values of c g makes the damage evolution behave differently where D may behave linearly or nonlinearly. In other words, the damage D evolution is predicted in FEA by the parameter c g . Thus, the damage D may have multi-evolution trends: linear when c g is approaching to zero (i.e., c g = 0.0001), nonlinear concave curve when c g < 0, and nonlinear convex curve when c g > 0 (all different evolution trends are depicted in Figure 3). It is worth noting that a nonzero (c g 0) should be set to avoid mathematical singularity in the model.
On the other hand, damage evolution due to cycling loading was seen to behave differently. In Reference [1], a good prediction of ELCF life, using Equation (7), was obtained when c g < 0.
The fatigue life prediction shows good agreement with the experimental results. Hence, the damage accumulation of the specimen's life under monotonic loading becomes nonlinear ( Figure 9) and is in good agreement with the experiment results. Accordingly, the new unified model introduced in Reference [1] predicts both ductile fracture and ELCF. accumulation of the specimen's life under monotonic loading becomes nonlinear ( Figure 9) and is in good agreement with the experiment results. Accordingly, the new unified model introduced in Reference [1] predicts both ductile fracture and ELCF. Figure 9. Nonlinear damage evolution for monotonic loading for four specimens at different stress states when introducing the parameter Cg in Equation (7).
The damage evolution during the high strain cyclic loading for the two "R1" specimens, in the FEA, was seen to increase rapidly in the first couple of cycles until the accumulated equivalent plastic strain reaches around 0.5. Afterwards, the damage evolution starts to accumulate slowly as the cycling process continues. This slow evolution in damage becomes even slower as the specimen reaches its life limit. In Figure 10a, the specimen fails due to fatigue at four high-strain reversed cycles (reminder: the failure occurs when damage accumulation reaches unity). During the first reversed cycle, the damage evolves from zero to approximately 0.67 as the equivalent plastic strain reaches 0.12. In the second cycle, the damage accumulation evolves to about 0.90 as the equivalent plastic strain reaches 0.45. During cycles 3 and 4, the damage accumulation slightly increases to 0.97 and 1, respectively. Similarly, Figure 10b simulates a similar damage evolution behavior for nine high-strain reversed cycles. The end of the curve line represents the fatigue full failure (D = 1). In Figure 10a,b, the damage evolution during the tension loading in the first cycle is represented by (i). Similarly, damage evolution during the tension loading in the second cycle is represented by (iii), whereas the damage evolution during the compression loading in the first cycle is represented by (ii). It is seen that the damage evolution during compression loading is very small (≈0.0001) and almost negligible. The same applies for the remaining number of reverse loading cycles.
For the three "R2" specimens ( Figure 11a-c), the damage evolution was seen to increase in a similar process as in the "R1" specimens. In Figure 11a, the number of high-strain reversed cycles to failure was ten cycles. The number of reverse loading cycles to failure in Figure 11b,c are 41 and 51 Figure 9. Nonlinear damage evolution for monotonic loading for four specimens at different stress states when introducing the parameter C g in Equation (7).
The damage evolution during the high strain cyclic loading for the two "R1" specimens, in the FEA, was seen to increase rapidly in the first couple of cycles until the accumulated equivalent plastic strain reaches around 0.5. Afterwards, the damage evolution starts to accumulate slowly as the cycling process continues. This slow evolution in damage becomes even slower as the specimen reaches its life limit. In Figure 10a, the specimen fails due to fatigue at four high-strain reversed cycles (reminder: the failure occurs when damage accumulation reaches unity). During the first reversed cycle, the damage evolves from zero to approximately 0.67 as the equivalent plastic strain reaches 0.12. In the second cycle, the damage accumulation evolves to about 0.90 as the equivalent plastic strain reaches 0.45. During cycles 3 and 4, the damage accumulation slightly increases to 0.97 and 1, respectively. Similarly, Figure 10b simulates a similar damage evolution behavior for nine high-strain reversed cycles. The end of the curve line represents the fatigue full failure (D = 1). accumulation of the specimen's life under monotonic loading becomes nonlinear ( Figure 9) and is in good agreement with the experiment results. Accordingly, the new unified model introduced in Reference [1] predicts both ductile fracture and ELCF. Figure 9. Nonlinear damage evolution for monotonic loading for four specimens at different stress states when introducing the parameter Cg in Equation (7).
The damage evolution during the high strain cyclic loading for the two "R1" specimens, in the FEA, was seen to increase rapidly in the first couple of cycles until the accumulated equivalent plastic strain reaches around 0.5. Afterwards, the damage evolution starts to accumulate slowly as the cycling process continues. This slow evolution in damage becomes even slower as the specimen reaches its life limit. In Figure 10a, the specimen fails due to fatigue at four high-strain reversed cycles (reminder: the failure occurs when damage accumulation reaches unity). During the first reversed cycle, the damage evolves from zero to approximately 0.67 as the equivalent plastic strain reaches 0.12. In the second cycle, the damage accumulation evolves to about 0.90 as the equivalent plastic strain reaches 0.45. During cycles 3 and 4, the damage accumulation slightly increases to 0.97 and 1, respectively. Similarly, Figure 10b simulates a similar damage evolution behavior for nine high-strain reversed cycles. The end of the curve line represents the fatigue full failure (D = 1). In Figure 10a,b, the damage evolution during the tension loading in the first cycle is represented by (i). Similarly, damage evolution during the tension loading in the second cycle is represented by (iii), whereas the damage evolution during the compression loading in the first cycle is represented by (ii). It is seen that the damage evolution during compression loading is very small (≈0.0001) and almost negligible. The same applies for the remaining number of reverse loading cycles.
For the three "R2" specimens ( Figure 11a-c), the damage evolution was seen to increase in a similar process as in the "R1" specimens. In Figure 11a, the number of high-strain reversed cycles to failure was ten cycles. The number of reverse loading cycles to failure in Figure 11b,c are 41 and 51 In Figure 10a,b, the damage evolution during the tension loading in the first cycle is represented by (i). Similarly, damage evolution during the tension loading in the second cycle is represented by (iii), whereas the damage evolution during the compression loading in the first cycle is represented by (ii). It is seen that the damage evolution during compression loading is very small (≈0.0001) and almost negligible. The same applies for the remaining number of reverse loading cycles.
For the three "R2" specimens ( Figure 11a-c), the damage evolution was seen to increase in a similar process as in the "R1" specimens. In Figure 11a, the number of high-strain reversed cycles to failure was ten cycles. The number of reverse loading cycles to failure in Figure 11b,c are 41 and 51 cycles, respectively. The damage accumulation for the specimens in Figure 11b,c rapidly increase in the first couple of cycles then decrease to a point where the increment becomes approximately 0.01 in the final counts of cycles. Lastly, the three PE specimens (Figure 12a-c) depict a similar damage accumulation behavior of R2 specimens. The number of reverse loading cycles to failure in Figure 12a-c are 10, 21, and 43 cycles, respectively.
In summary, the damage evolution is shown to behave nonlinearly. This result verifies that the damage evolution in ELCF is very high in the first reverse cycles and decreases as the cycling loading continues until failure. Thus, the damage evolution rate during ELCF is nonlinear. Accordingly, we expect an excessive deterioration in the material's microstructure during the early cycles of ELCF and less deterioration in the remaining cycles until failure. In between the early number of cycles and the crack initiation, a transition phase in the damage evolution is seen in the damage curve, where the evolution changes from high and rapid to low and slow. Figure 13 shows the specimens' deformation in FEA during the cyclic loading and the damage accumulation within the specimen' center. The highest damage accumulation (in red) during ELCF is where the voids coalesce and therefore a crack initiate.
cycles, respectively. The damage accumulation for the specimens in Figure 11b,c rapidly increase in the first couple of cycles then decrease to a point where the increment becomes approximately 0.01 in the final counts of cycles. Lastly, the three PE specimens (Figure 12a-c) depict a similar damage accumulation behavior of R2 specimens. The number of reverse loading cycles to failure in Figure  12a-c are 10, 21, and 43 cycles, respectively.
In summary, the damage evolution is shown to behave nonlinearly. This result verifies that the damage evolution in ELCF is very high in the first reverse cycles and decreases as the cycling loading continues until failure. Thus, the damage evolution rate during ELCF is nonlinear. Accordingly, we expect an excessive deterioration in the material's microstructure during the early cycles of ELCF and less deterioration in the remaining cycles until failure. In between the early number of cycles and the crack initiation, a transition phase in the damage evolution is seen in the damage curve, where the evolution changes from high and rapid to low and slow. Figure 13 shows the specimens' deformation in FEA during the cyclic loading and the damage accumulation within the specimen' center. The highest damage accumulation (in red) during ELCF is where the voids coalesce and therefore a crack initiate.  cycles, respectively. The damage accumulation for the specimens in Figure 11b,c rapidly increase in the first couple of cycles then decrease to a point where the increment becomes approximately 0.01 in the final counts of cycles. Lastly, the three PE specimens (Figure 12a-c) depict a similar damage accumulation behavior of R2 specimens. The number of reverse loading cycles to failure in Figure  12a-c are 10, 21, and 43 cycles, respectively. In summary, the damage evolution is shown to behave nonlinearly. This result verifies that the damage evolution in ELCF is very high in the first reverse cycles and decreases as the cycling loading continues until failure. Thus, the damage evolution rate during ELCF is nonlinear. Accordingly, we expect an excessive deterioration in the material's microstructure during the early cycles of ELCF and less deterioration in the remaining cycles until failure. In between the early number of cycles and the crack initiation, a transition phase in the damage evolution is seen in the damage curve, where the evolution changes from high and rapid to low and slow. Figure 13 shows the specimens' deformation in FEA during the cyclic loading and the damage accumulation within the specimen' center. The highest damage accumulation (in red) during ELCF is where the voids coalesce and therefore a crack initiate.

Fatigue Crack Initiation and Propagation
Studying fatigue, over the years, essentially involves analyzing the loss of a metal's strength and predicting a metal's life. This approach is essential in studying metal fatigue to help comprehend the crack initiation and propagation that occurs during the cyclic loading process. In high-cycle fatigue (HCF) (Nf > 10 4 ) and low-cycle fatigue (LCF) (10 4 > Nf > 10 2 ) regimes, cracks initiate and progressively propagate until failure. In these two regimes, the fatigue crack mechanism starts once a metal's microstructure deforms due to repetitive application of stress. After that, microcracks initiate at the surface and then coalesce to form larger cracks. The final metal failure occurs as the cracks propagate inwards [38].
On the other hand, the phenomena of crack initiation and propagation (aforementioned) for HCF and LCF does not uniquely take a role in the ELCF regime. Researchers claim that the ductile fracture failure mechanism is part of the ELCF failure as well as the fatigue failure mechanism. However, the ductile fracture mechanism takes the principal and major role in ELCF failure [39].
Experimental observation reveals strong resemblances between the fracture surfaces due to ELCF and the ductile fracture. Studies presumed this resemblance as a result of the high strain amplitudes applied that lead to significant ductility exhaustion within the metal. Moreover, fracture surfaces due to ELCF show no signs of beach marks (clamshell marks). Nevertheless, fracture surfaces as a result of ELCF display a cup-cone fracture mode for the axisymmetric specimens "R2" and "R1" and slant fracture mode for the plane-strain specimens "PE" (Figure 14), similar to fractured specimens due to ductile fracture in [23]. This observation supports the notion that cracks due to ELCF tend to initiate in the gauge center of the specimen and propagate outwards towards the surface. In Figures 15-17, SEM pictures of the fracture surface profile show microvoids growth and coalescence due to the very high plastic strain during ELCF for "R1", "R2", and "PE" specimens. In the same vein, no signs of striations nor fatigue-fracture transition on the fracture surfaces are observed under high magnification.

Fatigue Crack Initiation and Propagation
Studying fatigue, over the years, essentially involves analyzing the loss of a metal's strength and predicting a metal's life. This approach is essential in studying metal fatigue to help comprehend the crack initiation and propagation that occurs during the cyclic loading process. In high-cycle fatigue (HCF) (N f > 10 4 ) and low-cycle fatigue (LCF) (10 4 > N f > 10 2 ) regimes, cracks initiate and progressively propagate until failure. In these two regimes, the fatigue crack mechanism starts once a metal's microstructure deforms due to repetitive application of stress. After that, microcracks initiate at the surface and then coalesce to form larger cracks. The final metal failure occurs as the cracks propagate inwards [38].
On the other hand, the phenomena of crack initiation and propagation (aforementioned) for HCF and LCF does not uniquely take a role in the ELCF regime. Researchers claim that the ductile fracture failure mechanism is part of the ELCF failure as well as the fatigue failure mechanism. However, the ductile fracture mechanism takes the principal and major role in ELCF failure [39].
Experimental observation reveals strong resemblances between the fracture surfaces due to ELCF and the ductile fracture. Studies presumed this resemblance as a result of the high strain amplitudes applied that lead to significant ductility exhaustion within the metal. Moreover, fracture surfaces due to ELCF show no signs of beach marks (clamshell marks). Nevertheless, fracture surfaces as a result of ELCF display a cup-cone fracture mode for the axisymmetric specimens "R2" and "R1" and slant fracture mode for the plane-strain specimens "PE" (Figure 14), similar to fractured specimens due to ductile fracture in [23]. This observation supports the notion that cracks due to ELCF tend to initiate in the gauge center of the specimen and propagate outwards towards the surface. In Figures 15-17, SEM pictures of the fracture surface profile show microvoids growth and coalescence due to the very high plastic strain during ELCF for "R1", "R2", and "PE" specimens. In the same vein, no signs of striations nor fatigue-fracture transition on the fracture surfaces are observed under high magnification.

Conclusions
Predicting and evaluating the damage evolution caused by ELCF with mathematical models is crucial and challenging. An extended damage accumulation model for ELCF life prediction was developed with the aid of experimental data in Reference [1]. Four distinctive geometries were designed to study the effect of different stress states on ELCF life and damage evolution. The four specimen's geometries were designed in such a way to attain different stress states, three of which are classified as axisymmetric stress states and one is classified as a plane-strain stress state. Investigating multiaxial damage evolution due to ELCF can be performed by inspecting three distinguished features on the specimen's fracture surface: the crack initiation site, the crack growth surface area, and the final fractured surface. The irreversible damage mechanism is quantitatively represented by the accumulated damage D in Equations (6) and (7). The damage model is calibrated for life prediction to agree with the ELCF experimental results. After that, the damage evolution is extracted from Abaqus and plotted versus the accumulated Ɛ . The results are as follows: 1. The results of the damage evolution plot show that the ELCF damage behaves in a nonlinear manner. 2. The ELCF damage evolution for all stress states (Figures 10-12) was seen to increase rapidly in the first couple of cycles until the accumulated equivalent plastic strain reaches a certain value then it starts to accumulate slowly as the cycling process continues. This slow evolution in damage becomes even slower as the specimen reaches its life limit (D = 1). 3. There is excessive deterioration in the metal microstructure during the early cycles and less deterioration in the remaining cycles until failure. 4. In all stress states, the damage evolution during compression loading is very small (≈0.0001) and almost negligible. 5. This FEA simulation of the damage accumulation models the crack initiation location in the specimens' center similar to what have been seen in the experimental results. 6. The fracture surface shows that the ductile fracture failure mechanism is part of the ELCF failure as well as the fatigue failure mechanism. Moreover, it is observed that the ductile fracture mechanism takes the dominant role in very high strain cyclic loading. 7. SEM pictures demonstrate microvoids nucleation, growth, and coalescence, with no signs of striations due to ELCF.
Finally, it is recommended for future-related research to increase the number of test specimens to replicate the test results in order to study the variability associated with the observations and conclusions. It is also recommended to use transmission electron microscopy (TEM) to perform extra examinations for the fractured surfaces to add additional insight into the dislocation configurations

Conclusions
Predicting and evaluating the damage evolution caused by ELCF with mathematical models is crucial and challenging. An extended damage accumulation model for ELCF life prediction was developed with the aid of experimental data in Reference [1]. Four distinctive geometries were designed to study the effect of different stress states on ELCF life and damage evolution. The four specimen's geometries were designed in such a way to attain different stress states, three of which are classified as axisymmetric stress states and one is classified as a plane-strain stress state. Investigating multiaxial damage evolution due to ELCF can be performed by inspecting three distinguished features on the specimen's fracture surface: the crack initiation site, the crack growth surface area, and the final fractured surface. The irreversible damage mechanism is quantitatively represented by the accumulated damage D in Equations (6) and (7). The damage model is calibrated for life prediction to agree with the ELCF experimental results. After that, the damage evolution is extracted from Abaqus and plotted versus the accumulated ε pl . The results are as follows: 1.
The results of the damage evolution plot show that the ELCF damage behaves in a nonlinear manner.

2.
The ELCF damage evolution for all stress states (Figures 10-12) was seen to increase rapidly in the first couple of cycles until the accumulated equivalent plastic strain reaches a certain value then it starts to accumulate slowly as the cycling process continues. This slow evolution in damage becomes even slower as the specimen reaches its life limit (D = 1).

3.
There is excessive deterioration in the metal microstructure during the early cycles and less deterioration in the remaining cycles until failure.

4.
In all stress states, the damage evolution during compression loading is very small (≈0.0001) and almost negligible.

5.
This FEA simulation of the damage accumulation models the crack initiation location in the specimens' center similar to what have been seen in the experimental results. 6.
The fracture surface shows that the ductile fracture failure mechanism is part of the ELCF failure as well as the fatigue failure mechanism. Moreover, it is observed that the ductile fracture mechanism takes the dominant role in very high strain cyclic loading. 7.
SEM pictures demonstrate microvoids nucleation, growth, and coalescence, with no signs of striations due to ELCF.
Finally, it is recommended for future-related research to increase the number of test specimens to replicate the test results in order to study the variability associated with the observations and conclusions. It is also recommended to use transmission electron microscopy (TEM) to perform extra examinations for the fractured surfaces to add additional insight into the dislocation configurations in ELCF specimens. In addition, in situ X-ray tomography during the cycling loading tests can further assist in monitoring and observing the fracture mode and failure mechanism [40][41][42].  -159-829-1441). The authors, therefore, gratefully acknowledge the DSR technical and financial support.

Conflicts of Interest:
The authors declare no conflict of interest.