Application of Differential Entropy in Characterizing the Deformation Inhomogeneity and Life Prediction of Low-Cycle Fatigue of Metals

The relation between deformation inhomogeneity and low-cycle-fatigue failure of T2 pure copper and the nickel-based superalloy GH4169 under symmetric tension-compression cyclic strain loading is investigated by using a polycrystal representative volume element (RVE) as the material model. The anisotropic behavior of grains and the strain fields are calculated by crystal plasticity, taking the Bauschinger effect into account to track the process of strain cycles of metals, and the Shannon’s differential entropies of both distributions of the strain in the loading direction and the first principal strain are employed at the tension peak of the cycles as measuring parameters of strain inhomogeneity. Both parameters are found to increase in value with increments in the number of cycles and they have critical values for predicting the material’s fatigue failure. Compared to the fatigue test data, it is verified that both parameters measured by Shannon’s differential entropies can be used as fatigue indicating parameters (FIPs) to predict the low cycle fatigue life of metal.


Introduction
Extensive studies have been conducted in computational mechanics, materials science, and physics to investigate the fatigue lifetime problem and the deformation mechanisms of metals subjected to alternating stress and strain. These have improved our understanding of the fatigue resistance and service life of metallic engineering components and structures. However, the established and widely applied estimation methods, formulas and analyses for fatigue life prediction are based on the conclusions drawn from extensive fatigue experiments [1,2]. Among the competing approaches, Basquin's formula [3] describes the relationship between the stress state and cyclic fatigue life, and the 5] characterizes the relation between plastic strain and cyclic fatigue life. Both formulas, as well as their combination [6], rely on the lifetime measurement gained in cyclic fatigue experiments under the condition of controlled stress or strain, and the parameters are numerically fitted to match the measurements. These kinds of empirical methods, which generally depend on measured fatigue life data, are still widely used [7][8][9][10][11][12][13][14].
Considerable efforts are, therefore, underway towards establishing an analytical model which is capable of reflecting the evolution of fatigue damage in a material undergoing cyclic loading until final failure. For instance, models based on accumulated plastic strain or dissipation energy describe the degree of material damage at the macro-scale [15][16][17][18][19][20], but since they lack any interpretation

Methodology of Strain Inhomogeneity Analysis at the Grain Level for a Polycrystal under Cyclic Loading
The conventional analysis of a polycrystalline metal considers it as a continuous medium, whereas actually the material has a complex micro-structure. Hence, the deformation field is non-uniform at the grain level and becomes increasingly uneven with increasing numbers of cycles. In the present paper, the polycrystalline structure of the material is taken into account; thus, the inhomogeneous elastic-plastic deformation at the grain level and its evolution with loading cycles can be calculated. To avoid any excess of computation time, information on second phases and other micro-structural details within the grain and at the area of the grain boundaries are not included in the model.

Modeling the Material as Representative Volume Element
With reference to [53,54], a Voronoi polyhedron aggregation is employed as RVE. The shape, size and crystal orientation of the grains in this model are generated randomly. The plastic deformation of the grains is described by the slip driven by the resolved shear stresses acting at the respective slip systems. The RVE contains 27,000 8-node hexahedral elements and 29,791 nodes, as shown in Figure 1. It is subjected to symmetrical tensile-compressive loading cycles and its surfaces are assumed to remain plane during deformation.
Materials 2018, 11, x FOR PEER REVIEW 3 of 23 deformation inhomogeneity of a polycrystalline metal with loading cycles can be represented by the numerical change of entropy, and this change is related to the microstructure evolution of the metal. Thereby, we hope to further understand if it can be adopted as a FIP to predict fatigue failure of materials. To do this, it is necessary to verify whether the calculation of plastic deformation tracking the cycles up to fatigue failure can be carried out, and the entropy takes a critical value that is independent of the load amplitude within a certain range.
In the present paper, we focus on the low-cycle fatigue of polycrystalline metals under strain-controlled symmetrical tensile-compressive cyclic loading, and investigate T2 pure copper at room temperature and the nickel-based superalloy GH4169 at a temperature of 650 °C.
The outline of the paper is as follows: (1) the Shannon's differential entropies of longitudinal strain and first principal strain distributions are calculated by crystal plasticity for a polycrystalline RVE subjected to cyclic loading at different strain amplitudes. (2) Their suitability as FIPs is verified by comparing them with the standard deviation of longitudinal strain and the mean of the first principal strain, and their ability for predicting low-cycle fatigue life is discussed.

Methodology of Strain Inhomogeneity Analysis at the Grain Level for a Polycrystal under Cyclic Loading
The conventional analysis of a polycrystalline metal considers it as a continuous medium, whereas actually the material has a complex micro-structure. Hence, the deformation field is non-uniform at the grain level and becomes increasingly uneven with increasing numbers of cycles. In the present paper, the polycrystalline structure of the material is taken into account; thus, the inhomogeneous elastic-plastic deformation at the grain level and its evolution with loading cycles can be calculated. To avoid any excess of computation time, information on second phases and other micro-structural details within the grain and at the area of the grain boundaries are not included in the model.

Modeling the Material as Representative Volume Element
With reference to [53], a Voronoi polyhedron aggregation is employed as RVE. The shape, size and crystal orientation of the grains in this model are generated randomly. The plastic deformation of the grains is described by the slip driven by the resolved shear stresses acting at the respective slip systems. The RVE contains 27,000 8-node hexahedral elements and 29,791 nodes, as shown in Figure  1. It is subjected to symmetrical tensile-compressive loading cycles and its surfaces are assumed to remain plane during deformation. The macroscopic Cauchy stress tensor Σ and logarithmic strain tensor E are defined as the mean of the local Cauchy stresses σ and logarithmic strains ε over the RVE and calculated as load per unit area and displacements of the surfaces. Crystal plasticity is employed as a constitutive relation between local Cauchy stresses and logarithmic strains within each grain. The macroscopic Cauchy stress tensor Σ and logarithmic strain tensor E are defined as the mean of the local Cauchy stresses σ and logarithmic strains ε over the RVE and calculated as load per unit area and displacements of the surfaces. Crystal plasticity is employed as a constitutive relation between local Cauchy stresses and logarithmic strains within each grain.

Constitutive Equations of Crystal Plasticity
Plastic slip in the grains is described by the constitutive equations of crystal plasticity proposed by Zhang et al. [53,55], which are based on the fundamental theoretical framework established by Hill and Rice [56], Asaro and Rice [57], Peirce et al. [58] and Needleman et al. [59]. However, the latter did not consider the plastic behavior under cyclic loading.
Extending Hutchinson's [60] power relation between the shear strain rate, . γ (α) , at a slip system, α, and the resolved shear stress, τ (α) , Feng et al. [61] introduced back-stresses and nonlinear kinematic hardening, referring to the Armstrong-Frederick model [62] of the Bauschinger effect, where . γ 0 is the reference strain rate, being a constant for all slip systems. The resolved back-stress on the α-slip system is denoted by x (α) , k' is the rate sensitivity parameter and g (α) defines the domain where the material behaves elastic; its evolution is described as [63] .
h αβ (γ) denote the hardening moduli proposed by Hutchinson [64], q is a constant and h(γ) is given according to Chang and Asaro [65] as where h 0 is the initial hardening rate, τ 0 is the initial critical resolved shear stress, and τ s is the saturation value. These parameters are regarded as material constants. The evolution of back-stresses, x (α) , is introduced in [55], where a, c, e 1 , e 2 and λ are material constants. This formulation includes a strain hardening term, a dynamic recovery term, and a static recovery term. The identification of material constants in Equations (1), (3)-(5) is based on cyclic tests combined with numerical simulations. The initial unit normal vector of the slip surface of the slip system α and the unit vector of the initial direction of the system are denoted as n (α) and m (α) , respectively. Referring to Hill and Rice [56], Asaro and Rice [57], Peirce et al. [58] and Needleman et al. [59], the Schmid tensor, which establishes the relationship between the shear strain and shear stress in the slip system and its corresponding strain and stress in the Cartesian coordinate system, is given by: with m (α) * = F * ·m (α) , n (α) * = n (α) ·F * −1 ; F * is the elastic term of the deformation gradient tensor, F. The plastic deformation rate tensor can be calculated as Materials 2018, 11, 1917 5 of 22 and the resolved shear stress by Schmid's law, Assuming that the elastic deformations are small, the constitutive relation can be expressed as where . σ J is the Jaumann rate of Cauchy stress,

<4>
C is the fourth-order elasticity tensor with respect to the global coordinate axes. During the calculation, the crystal coordinate axes for each grain are rotated along with the changing configuration according to the lattice rotation, since the global coordinate system is fixed.
The incremental change of the Cauchy stress tensor can be calculated as where the increments ∆σ J , ∆ε and ∆ε p are determined by integrating the corresponding rates . σ J , D and D p . For specific numerical implementation as user-supplied subroutine UMAT in the FE code ABAQUS [66], see Zhang et al. [53].

Characterization of the Deformation Inhomogeneity by Statistical Parameters
Due to the randomness of grain shape, size and orientation the strain distribution within the RVE appears to be randomly distributed. As mentioned in the introduction, the deformation of a polycrystalline metal is non-uniform even under uniform macroscopic loading. The distribution of strains at the grain level can be described by statistical parameters like the average of strain components, k v k − ε ij 2 , and the maximum of local strains, ε max ij = max n RVE ε ij ; where n RVE is the total number of finite elements in the RVE, ε ij are the components of the local logarithmic strain tensor, and v k = ∆V k /V RVE , with ∆V k being the volume of k-th element and V RVE the total volume of the RVE. The strain distribution in the polycrystalline RVE changes with the loading cycles and so do the statistical parameters defined above. Zhang et al. [53] proposed to measure the inhomogeneity of micro-strains by the standard deviation,ε (ll) , (no summation over subscripts in brackets) of the local longitudinal strain, ε (ll) , which is the strain in the macroscopic loading direction and is also the direction of macroscopic first principal strain. Its value was found to grow with the number of cycles and can thus be related to the material fatigue failure. Further investigations demonstrated the same behavior for the average, ε I , and the maximum, ε max I , of first principal strains [54], which can hence be employed as FIPs and characterize low-cycle fatigue life by comparison with their respective critical values.
This present work attempts to quantify the inhomogeneity of the deformation by using the Shannon entropy calculated from the strain distribution. Entropy initially refers to a probability function characterizing the state of a thermodynamic system and is a measure of its disorder. Shannon and Weaver [48] transferred the concept to quantitative measurement of information and introduced the formula of the information entropy, where p i is the probability that an event will occur. The base of the logarithm can be taken as 10, e or 2, depending on the object investigated. Here, the value e is taken. In Equation (11a), the relative volume fraction of the subdivision region,v i = ∆V i /V RVE (with the strain interval ε i ≤ ε ≤ ε i+1 ) is analogized with p i . For a model such as the Voronoi aggregation in which the strain and stress are the continuous variables people can adopt the concept of entropy for continuous distributions which is also referred to as the differential entropy (or continuous entropy) and is given [47], where f x denotes the probability density function for a random variable. The differential entropy retains many of the properties of its discrete counterpart, but with some important differences that the details may refer to the literature [67]. In a RVE, the interval of the strain [ε min , ε max ] is divided into n subintervals ∆ε i , i = 1, 2 . . . , n. In addition, then the integral of the entropy described by Equation (11b) can be numerically calculated according to the following formula: where, p ε i , is the relative volume fraction, ∆V i /V RVE , of the region where ε i ≤ ε ≤ ε i+1 . For calculation convenience, the interval [ε i , ε i+1 ] in calculation is taken as a constant and ∆ε i = ε i+1 − ε i = (ε max − ε min )/n; n is the division number. It is necessary to point out that the value of the continuous entropy for strain distribution may be negative due to the term in Equation (12) being larger than 1, because generally the strain distribution range [ε min , ε max ] is very small. People used to assume entropy to be a positive value; therefore, we use 0.001ε (1000 µε) as the unit for measuring the strain interval in Equation (12); thus, the value of p ε i ∆ε i will be in the range [0, 1], and then the entropy must be positive. This treatment only changes the zero point (reference point), and does not change the entropy difference of different processes.
Thus, Equation (12) is used to describe the inhomogeneity of the distribution in the RVE. The inhomogeneity of the distributions of local longitudinal strain ε (ll) and local first principal strain ε I , respectively, are measured by the Shannon's differential entropies H ε (ll) and H ε I , where p ε (ll)i is the relative volume fraction, ∆V i /V RVE , of the region where ε (ll)i−1 ≤ ε (ll) ≤ ε (ll)i , and p ε 1i is the relative volume fraction of the region where ε Ii−1 ≤ ε I ≤ ε Ii , respectively. Equation (13) provides a measure of the deformation inhomogeneity; when the deformation is completely uniform the entropy is zero, generally it is greater than zero, and the more uneven the deformation is, the larger is H. The calculated results by using Equation (13) are convergent with the division number, n, which is verified by taking a Gaussian distribution function; when the division number is no less than 20 the calculated entropy result will be very close to the saturation value and the error can be ignored.
To verify the rationality of H ε (ll) and H ε I as the parameters characterizing the deformation inhomogeneity and as FIPs, the low-cycle fatigue predictions of the pure copper and the nickel-based superalloy GH4169 are performed by applying both FIPs. In addition, they are further compared with that using the parametersε (ll) and ε I , which have been proposed as FIPs in previous publications [53,54].

Results
For the pure copper and the nickel-based superalloy GH4169, the fatigue experiments are conducted under symmetrical cyclic strain-controlled loading, at room temperature for copper and at a temperature of 650 • C for GH4169. Both materials have the fcc lattice and they are approximated as fcc-equiaxed crystals [53,54]. The effect of grain size is ignored in the material model, and the same geometric partitioning of the RVE (see Figure 1) is used for both. Due to the presence of grain boundaries, inhomogeneous distribution of dispersed phase particles and defects, the very strong strain gradient will lead to a grain size effect. It is difficult to establish a quantitative analysis model based on these complicated factors. Therefore, we regard each grain to be an equivalent crystal with an ideal lattice, and establish a simplified polycrystalline RVE model without size effect, by which the numerical simulation of the whole cycle process can be performed and the difference in grain size for different metals can be ignored.
The chemical compositions of the pure copper T2 and the alloy GH4169 are shown in Tables 1  and 2, respectively, and the mechanical properties at room temperature for copper and at 650 • C for GH4169 are exhibited in Tables 3 and 4, respectively.  Table 3. Mechanical properties of the pure copper T2.

Young's Modulus (E) GPa
Yield Stress Based on the finite-element simulation of the RVE of the polycrystals, the parameters of the crystal-plastic model of the materials can be calibrated according to the tested hysteresis-loops. The material model parameters for the pure copper and the cast nickel-base superalloy GH4169 are displayed in Tables 5 and 6, respectively, according to references [53,54], where the detailed procedure for parameter calibration can be found. Table 5. Elastic constants and crystal plasticity parameters of the pure copper T2 [53].

Elastic Constants
Material Parameters of the Crystal Viscoplastic Model  Table 6. Elastic constants and crystal plasticity parameters of the GH4169 superalloy [54]. Applying the parameters of Tables 5 and 6, the stable hysteresis-loops of copper T2 for strain amplitudes 0.003, 0.004, 0.005 and 0.006 and of nickel-based superalloy GH4169 for strain amplitudes 0.0045, 0.006, 0.008, 0.009, 0.01 and 0.013, respectively, as calculated by the RVE simulations are shown in Figure 2. The actually measured steady hysteresis-loops are plotted for comparison. The numerically simulated results accord well with the test data.

Elastic Constants
Material Parameters of the Crystal Viscoplastic Model  Applying the parameters of Tables 5 and 6, the stable hysteresis-loops of copper T2 for strain amplitudes 0.003, 0.004, 0.005 and 0.006 and of nickel-based superalloy GH4169 for strain amplitudes 0.0045, 0.006, 0.008, 0.009, 0.01 and 0.013, respectively, as calculated by the RVE simulations are shown in Figure 2. The actually measured steady hysteresis-loops are plotted for comparison. The numerically simulated results accord well with the test data.

Strain Distribution in RVE with Increasing Numbers of Cycles
If the factors characterizing the polycrystalline structure-namely, grain orientation and anisotropy properties-are considered in the model, the deformation of the material is not uniform, even under uniform macroscopic loading. Figure 3 displays the contours of the local longitudinal strain, ) (ll ε , i.e., the normal strain along macroscopic loading direction, and the local first principal strain, I ε , at the second and 2787th tensile peak for a strain amplitude a E = 0.004. It can be seen from Figure 3a,c that the distribution of the longitudinal strains at the 2787th cycle is much more uneven than that at the second cycle, which implies that the difference of maximum and minimum values in the RVE has become much larger. The same holds for the first principal strain, see Figure  3b,d.

Strain Distribution in RVE with Increasing Numbers of Cycles
If the factors characterizing the polycrystalline structure-namely, grain orientation and anisotropy properties-are considered in the model, the deformation of the material is not uniform, even under uniform macroscopic loading. Figure 3 displays the contours of the local longitudinal strain, ε (ll) , i.e., the normal strain along macroscopic loading direction, and the local first principal strain, ε I , at the second and 2787th tensile peak for a strain amplitude E a = 0.004. It can be seen from Figure 3a,c that the distribution of the longitudinal strains at the 2787th cycle is much more uneven than that at the second cycle, which implies that the difference of maximum and minimum values in the RVE has become much larger. The same holds for the first principal strain, see Figure 3b,d.
It needs to be pointed out that the strain contours are differently scaled and the plots cannot clearly visualize the changes from the 2nd to 2787th cycle. A quantitative evaluation of the strain distribution in the RVE and its change with cyclic loading is presented in Figure 4, showing the volume fractions of both the longitudinal strain and the first principal strain at the respective tensile peaks for increasing numbers of cycles number for different strain amplitudes (E a = 0.003, 0.004, 0.005 and 0.006). The variation of strains in the RVE becomes larger and larger and the non-zero range expands with increasing number of cycles.     It needs to be pointed out that the strain contours are differently scaled and the plots cannot clearly visualize the changes from the 2nd to 2787th cycle. A quantitative evaluation of the strain distribution in the RVE and its change with cyclic loading is presented in Figure 4

Predicting Low-Cycle Fatigue Failure by Entropy
By applying Equation (12), we calculate and plot the entropy curves, H ε (ll) (E a , N) and H ε 1 (E a , N), as shown in Figure 5. The horizontal axes are linearly scaled in Figure 5a,b, and logarithmically scaled in Figure 5c , as shown in Figure 5. The horizontal axes are linearly scaled in Figure 5a,b, and logarithmically scaled in Figure 5c,d. The information entropies of longitudinal strain and the first principal strain distributions always increase with the number of cycles. The greater the strain amplitude, the faster the entropy values increase. On the logarithmically scaled x-axes, the H(N) curves tend to become straight with similar asymptotic slopes for high numbers of cycles. According to 53, based on the results of symmetric tensile-compressive cyclic experiments and respective calculations of pure copper, the statistical standard deviation of longitudinal strain over the RVE can be applied as a FIP for the fatigue life assessment of a polycrystalline material. Likewise, the results for a GH4169 nickel-based superalloy verify that the statistical standard deviation of the longitudinal strain and the statistical mean of the first principal strain can also be used as FIPs 54. The suitability of a parameter to be used as a FIP is proved based on its ability to reach a critical value that indicates the occurrences of fatigue failure. In the following, the information entropy of the longitudinal strain and the first principal strain distribution will be examined in terms of whether they are suitable as FIPs.  Table 7 column 3) and the upper and lower values with respect to According to [53], based on the results of symmetric tensile-compressive cyclic experiments and respective calculations of pure copper, the statistical standard deviation of longitudinal strain over the RVE can be applied as a FIP for the fatigue life assessment of a polycrystalline material. Likewise, the results for a GH4169 nickel-based superalloy verify that the statistical standard deviation of the longitudinal strain and the statistical mean of the first principal strain can also be used as FIPs [54]. The suitability of a parameter to be used as a FIP is proved based on its ability to reach a critical value that indicates the occurrences of fatigue failure. In the following, the information entropy of the longitudinal strain and the first principal strain distribution will be examined in terms of whether they are suitable as FIPs.

The critical entropy values resulting from
The critical entropy values resulting from H ε (ll) f ∼ = H ε (ll) E a , N f and H ε 1 f ∼ = H ε 1 E a , N f , N f is the fatigue-life cycle number at the specific strain amplitude E a , and can be determined from Figure 5. For the respective strain amplitude, the average measured critical values with respect to N f (see Table 7 column 3) and the upper and lower values with respect to corresponding N f (see Table 7 column 2), obtained from the curves, are listed in Table 7 columns 4 and 5. To prove their rationality, we need to verify whether the critical entropy value determined at any strain amplitude can be used to reasonably predict the fatigue failure at another strain amplitude. Taking a critical value in Table 7 column 4 or 5 at any strain amplitude, one can get the intersects from the curves at different strain amplitudes in Figure 5 and get the corresponding fatigue-failure cycles, that is, a fatigue life prediction ( E a ∼ N f ) based on the test at single strain amplitude. In the same way, the upper and lower prediction can also be obtained by using the corresponding upper and lower values in Table 7 column 4 or 5. All the results using the present method using the critical value determined from the tests of different strain amplitudes to predict the fatigue lives of the material are listed in Table 8. In addition, the corresponding curves are showed in Figure 6. corresponding f N (see Table 7 column 2), obtained from the curves, are listed in Table 7 columns 4 and 5. To prove their rationality, we need to verify whether the critical entropy value determined at any strain amplitude can be used to reasonably predict the fatigue failure at another strain amplitude. Taking a critical value in Table 7 column 4 or 5 at any strain amplitude, one can get the intersects from the curves at different strain amplitudes in Figure 5 and get the corresponding fatigue-failure cycles, that is, a fatigue life prediction ( a f E~N ) based on the test at single strain amplitude. In the same way, the upper and lower prediction can also be obtained by using the corresponding upper and lower values in Table 7 column 4 or 5. All the results using the present method using the critical value determined from the tests of different strain amplitudes to predict the fatigue lives of the material are listed in Table 8. In addition, the corresponding curves are showed in Figure 6.   Table 7).
The predictions based on  Table 8 and Figure 6, since no matter the average critical value at which the strain amplitude test is taken, the prediction of the fatigue lives for other strain amplitudes are in agreement with the tests. Also, the scattered feature of the tested fatigue life can be illustrated by descriptions of the upper and lower curves in Figure 6. However, the tests are not enough to describe the scattered feature of the fatigue lives; in order to improve the predictions, more tests at the specified strain amplitude are needed for determining the FIP critical value, as well as the upper and lower bounds, more accurately.
The error between the predicted and the experimental fatigue lives is shown in Figure 7, including predictions based on standard deviation and mean values of longitudinal and principal  Table 7).
The predictions based on H ε (ll) f and H ε 1 f are proved to be rational in Table 8 and Figure 6, since no matter the average critical value at which the strain amplitude test is taken, the prediction of the fatigue lives for other strain amplitudes are in agreement with the tests. Also, the scattered feature of the tested fatigue life can be illustrated by descriptions of the upper and lower curves in Figure 6. However, the tests are not enough to describe the scattered feature of the fatigue lives; in order to improve the predictions, more tests at the specified strain amplitude are needed for determining the FIP critical value, as well as the upper and lower bounds, more accurately.
The error between the predicted and the experimental fatigue lives is shown in Figure 7, including predictions based on standard deviation and mean values of longitudinal and principal strains,ε (ll) , and ε I , data treated according to [53]. The horizontal axis and the vertical axis indicate the actually measured and the predicted life, respectively. The red solid line represents the ideal prediction, and the area between the two dashed lines is the interval of a factor of 2. The error of the prediction results is regarded as acceptable if the data points fall within the double factor region. Hence, Figure 7 shows that using parameters H ε (ll) f , H ε I f ,ε (ll) f , and ε I f for the strain amplitudes can give equally reasonable predictions of fatigue lives. Therefore, it is proved possib0le for this method to predict the LCF life of a metal using the present method with its FIP determined based on the test at one strain amplitude. Hence, Figure 7 shows that using parameters ε , and f I ε for the strain amplitudes can give equally reasonable predictions of fatigue lives. Therefore, it is proved possib0le for this method to predict the LCF life of a metal using the present method with its FIP determined based on the test at one strain amplitude.

Prediction of Fatigue Failure of a Nickel-Based Superalloy by Entropy
Although they both possess FCC lattices, the nickel-base superalloy GH4169 is highly unlike pure copper T2 in terms of mechanical properties, see the hysteresis loops in Figure 2, and the microstructure of GH4169 is far more complex than that of copper T2. By using the same method for analyzing the fatigue-life, we hope to further verify the suitability of Shannon's differential entropy as FIP not only for pure materials like copper, but also for alloy materials.

Entropy Increase with Cyclic Deformation of GH4169
Corresponding to different strain amplitudes, the entropy curves,

Prediction of Fatigue Failure of a Nickel-Based Superalloy by Entropy
Although they both possess FCC lattices, the nickel-base superalloy GH4169 is highly unlike pure copper T2 in terms of mechanical properties, see the hysteresis loops in Figure 2, and the microstructure of GH4169 is far more complex than that of copper T2. By using the same method for analyzing the fatigue-life, we hope to further verify the suitability of Shannon's differential entropy as FIP not only for pure materials like copper, but also for alloy materials.

Entropy Increase with Cyclic Deformation of GH4169
Corresponding to different strain amplitudes, the entropy curves, H ε (ll) (E a , N) and H ε 1 (E a , N), are calculated dependen on the number of load cycles based on the distribution of longitudinal strain and the distribution of the first principal strain as shown in Figure 8a,b. The horizontal axis is logarithmically scaled. Since the numbers of cycles span three powers of ten, display on a linear horizontal axis is no longer feasible. Similar to the previous analysis of pure copper, the entropy curves tend to a linear growth and the slopes are nearly equal. longitudinal strain and the distribution of the first principal strain as shown in Figure 8a,b. The horizontal axis is logarithmically scaled. Since the numbers of cycles span three powers of ten, display on a linear horizontal axis is no longer feasible. Similar to the previous analysis of pure copper, the entropy curves tend to a linear growth and the slopes are nearly equal.  Figure 9a,b. Due only one test for 4 of a total of 6 strain amplitudes having been performed, the upper and lower bounds for each FIP curve cannot be obtained, and the other 2 only had two tests conducted for each. However, the scattered feature of fatigue-life data can also be observed from this figure.  The curves H ε (ll) (E a , N) and H ε 1 (E a , N) for each specified strain amplitude E a are shown in Figure 8a,b; the critical values of the FIPs H ε (ll) f and H ε 1 f can be picked from each curve by specifying the average of the specific measurement of the fatigue life, these data are listed in Table 9. Then the fatigue lives at different strain amplitudes can be obtained by using any one of these critical values, further drawing a transverse line intersecting with the curves H ε (ll) (E a , N) or H ε 1 (E a , N). In addition, abscissa of intersection points on different curve are then taken to represent the estimated fatigue-life value (shown in Table 10) with respect to the specified strain amplitude. The predicted-life curves (E a ∼ N f ) based on the FIPs H ε (ll) f and H ε 1 f are separately displayed in Figure 9a,b. Due only one test for 4 of a total of 6 strain amplitudes having been performed, the upper and lower bounds for each FIP curve cannot be obtained, and the other 2 only had two tests conducted for each. However, the scattered feature of fatigue-life data can also be observed from this figure.
The fatigue-life predictions by applying FIPs H ε (ll) f and H ε I f together withε (ll) f and , which were suggested in the previous literatures [53,54], are all verified with test fatigue-life data and are shown in Figure 10. For the verification of the FIPsε (ll) f and ε 1 f in this figure, only the predictions based on the highest and lowest critical values are given for the sake of simplicity. From this figure, the predictions are observed to provide reasonable fatigue-life estimations over the range of strain amplitudes under consideration.  H ε I f . The test data in these figures are cited from literature [9] (see Table 9).

The Influence of Model Mesh Size on the Result of Entropy Calculation
The generation of metal polycrystalline RVE models involves many aspects, such as the number of grains, the division of finite element meshes, and randomly generated grain distributions and crystal orientations of grains. Due to random generation, the deformation fields of different polycrystal RVEs obtained by calculation are usually different, and it is necessary to verify whether approximately the same statistical deformation results can be obtained and that they are not sensitive to the random generation of the RVE model. Therefore, whether the model differences influence their description and result with respect to the fatigue law is verified preliminarily by a simple example. In the following, the effect on the entropy calculation caused by the different element size of the model mesh is discussed, which is performed to verify the rationality of the method employing entropy as a FIP.
Considering the RVE model of Section 2.1, the grain number and the orientation for each grain are kept unchanged, but the model is re-divided into a finite element mesh. The new model shown in Figure 11 is the hexahedron RVE divided into 20 × 20 × 20 equal parts along the three directions of length, width, and height, having 8,000 eight-node hexahedral elements and 9,261 nodes. Compared with the original model with 30 × 30 × 30 equal parts, 27,000 elements and 29,791 nodes, the present model, which is much lower in finite element number, will decrease computation time greatly. For the strain amplitudes E a involved in Figure 8, the entropy curves H ε (ll) (E a , N) and H ε 1 (E a , N) (red solid lines) using the model with less element number and the corresponding curves (black solid lines) using the original model are plotted together in Figure 12. The curves in the figure show that the difference in the size of the mesh has little effect on the entropy calculation result, which means that the difference of FEM mesh divisions will not cause a noticeable difference in the entropy calculation results.

The Influence of Model Mesh Size on the Result of Entropy Calculation
The generation of metal polycrystalline RVE models involves many aspects, such as the number of grains, the division of finite element meshes, and randomly generated grain distributions and crystal orientations of grains. Due to random generation, the deformation fields of different polycrystal RVEs obtained by calculation are usually different, and it is necessary to verify whether approximately the same statistical deformation results can be obtained and that they are not sensitive to the random generation of the RVE model. Therefore, whether the model differences influence their description and result with respect to the fatigue law is verified preliminarily by a simple example. In the following, the effect on the entropy calculation caused by the different element size of the model mesh is discussed, which is performed to verify the rationality of the method employing entropy as a FIP.
Considering the RVE model of Section 2.1, the grain number and the orientation for each grain are kept unchanged, but the model is re-divided into a finite element mesh. The new model shown in Figure 11 is the hexahedron RVE divided into 20 × 20 × 20 equal parts along the three directions of length, width, and height, having 8,000 eight-node hexahedral elements and 9,261 nodes. Compared with the original model with 30 × 30 × 30 equal parts, 27,000 elements and 29,791 nodes, the present model, which is much lower in finite element number, will decrease computation time curves (black solid lines) using the original model are plotted together in Figure 12. The curves in the figure show that the difference in the size of the mesh has little effect on the entropy calculation result, which means that the difference of FEM mesh divisions will not cause a noticeable difference in the entropy calculation results.

Conclusions
In the present paper, the Voronoi polycrystalline RVE is employed, combined with crystal plastic analysis, to simulate the tensile-compressive symmetric strain cycle of the materials T2 pure copper and nickel-base superalloy GH4169. The entropies calculated from of the distribution of the strain in the loading direction and from the distribution of the first principal strain of the materials are employed, respectively, to characterize the inhomogeneous deformation at the tension peak of the cyclic loading and to predict the fatigue lives of materials. Based on the analysis, the following conclusions are obtained: 1. The greater the strain amplitude a E , the larger the growing rates with cycles will be for the

Conclusions
In the present paper, the Voronoi polycrystalline RVE is employed, combined with crystal plastic analysis, to simulate the tensile-compressive symmetric strain cycle of the materials T2 pure copper and nickel-base superalloy GH4169. The entropies calculated from of the distribution of the strain in the loading direction and from the distribution of the first principal strain of the materials are employed, respectively, to characterize the inhomogeneous deformation at the tension peak of the cyclic loading and to predict the fatigue lives of materials. Based on the analysis, the following conclusions are obtained: 1.
The greater the strain amplitude E a , the larger the growing rates with cycles will be for the entropies H ε (ll) (E a , N) and H ε 1 (E a , N).

2.
Applying the critical values H ε (ll) f and H ε 1 f to determine the occurrence of metal low-cycle fatigue failure, the predictions are proved rational, and they are also approaches identical to that by applying the critical valuesε (ll) f and ε 1 f of the FIPs [53,54].

3.
Even in the absence of fatigue life data, the H ε (ll) (E a , N) and H ε 1 (E a , N) can be obtained by simulation depending only on the material parameters for crystal plasticity. Once the critical values of H ε (ll) f and H ε 1 f are determined by using the fatigue tests at only single strain amplitude, the fatigue lives can be predicted for other fatigue cycle at different strain amplitudes. 4.
The difference of the statistical results of H ε (ll) f and H ε 1 f from the models with different mesh sizes is very small. This proves that the Shannon's differential entropy calculation of strain using the present method is not sensitive to the mesh division.
It is necessary to point out that the investigation in the present paper only discusses the prediction of the fatigue lives of a metal under constant strain amplitude loading, and the condition for variable amplitude is not dealt with. In Figures 5 and 8, one can observe that the values of the functions H ε (ll) (E a , N) and H ε 1 (E a , N) are larger for the larger strain amplitude loading when the same relative life is consumed for different strain amplitude loading exerted on a material. Therefore, qualitative inferences on fatigue life in variable amplitude condition can be obtained that the life may be longer for such a strain amplitude load sequence from lower to higher, and vice versa. The method may also be considered to apply to multi-axial fatigue issue, but would need to consider more factors, like the influence of tri-axial stress state, how to characterize the non-proportional deformation, etc.