An Explicit Creep-Fatigue Model for Engineering Design Purposes

Background: Creep-fatigue phenomena are complex and difficult to model in ways that are useful from an engineering design perspective. Existing empirical-based models can be difficult to apply in practice, have poor accuracy, and lack economy. Need: There is a need to improve on the ability to predict creep-fatigue life, and do so in a way that is applicable to engineering design. Method: The present work modified the unified creep-fatigue model of Liu and Pons by introducing the parameters of temperature and cyclic time into the exponent component. The relationships between them were extracted by investigating creep behavior, and then a reference condition was introduced. Outcomes: The modified formulation was successfully validated on the materials of 63Sn37Pb solder and stainless steel 316. It was also compared against several other models. The results indicate that the explicit model presents better ability to predict fatigue life for both the creep fatigue and pure fatigue situations. Originality: The explicit model has the following beneficial attributes: Integration—it provides one formulation that covers the full range of conditions from pure fatigue, to creep fatigue, then to pure creep; Unified—it accommodates multiple temperatures, multiple cyclic times, and multiple metallic materials; Natural origin—it provides some physical basis for the structure of the formulation, in its consistency with diffusion-creep behavior, the plastic zone around the crack tip, and fatigue capacity; Economy—although two more coefficients were introduced into the explicit model, the economy is not significantly impacted; Applicability—the explicit model is applicable to engineering design for both manual engineering calculations and finite element analysis. The overall contribution is that the explicit model provides improved ability to predict fatigue life for both the creep-fatigue and pure-fatigue conditions for engineering design.


Introduction
Creep-fatigue damage is defined as the damage caused by reversed loading at elevated temperatures, hence combines the effects of fatigue and creep.This is a complex process since fatigue and creep behaviours are based on significantly different mechanisms at the microstructural level.Observationally, fatigue occurs via cracks through the grains, while creep involves the grain boundary cracking [1].The creep-fatigue phenomenon is relevant to a wide range of industries, such as aerospace, naval, nuclear and industries [2], hence cannot be ignored in engineering design.

The Design-Based Method
To provide an easier method for engineering practitioners to evaluate fatigue behaviour, a design-based method was proposed by Marin (Equation ( 1)) [3]: where S e is the endurance limit at the critical location of a machine part in the geometry and condition of use, S e is the rotary-beam test specimen endurance limit, k a is the surface condition modification factor, k b is the size modification factor, k c is the load modification factor, k d is the temperature modification factor, k e is the reliability factor, k f is the miscellaneous-effects modification factor.Engineering designers typically use this simple equation to determine the acceptable fatigue stress in a part.This modified endurance limit is based on the endurance limit at the reference condition and several multiplicative factors for surface condition, part size, type of load, operating temperature, etc.The only mechanical property included here is the endurance limit.This property can be related to ultimate tensile strength, such as the values of the endurance limit for steels are half of the ultimate tensile strength [4].Therefore, the benefits of this approach (Equation ( 1)) are ease of use since the tensile strength is readily known or easily measured.The detriments of this approach are that it only includes temperature effect when creep is active, and the fatigue evaluation is merely an approximation.Furthermore, all the modification factors have to be determined experimentally.Some degree of creep may be accommodated in the temperature factor, but the equation does not present a robust treatment for creep-fatigue.
Although the design-based method (Equation ( 1)) is simple enough for engineering practitioners, the poor accuracy is of concern.In addition, the consideration of multiple effects (such as shape, size and surface) is redundant if engineers merely aim to select materials.However, making further improvements to this formulation would not seem to be a viable way forward, since this numerical structure is only one of convenience rather than representing any deeper mechanics at the material science level.
There is a need for a more robust design method for creep-fatigue.Ideally such a method would have a formulation that directly related applied stress to life, included macrostructural rather than microstructural properties, and was economical to validate.The various attempts at addressing this problem are reviewed below.

The Conventional Empirical Methods
For mechanical design, a pre-evaluation of fatigue life (or damage) is normally applied at the initial stage of design to make a material selection or structural optimization.Normally, in the creep-fatigue situation, the total damage is numerically evaluated through the theory of damage accumulation and conventional-fatigue-based formulations.However, they present significant limitations.
Specifically, the creep-fatigue evaluation based on damage accumulation is normally conducted by the linear damage rule [1,5] or crack growth law [6], wherein the fatigue damage and creep damage are evaluated separately and then are numerically added.However, this is untrue to the physics of failure in that the fatigue and creep effects are not independent.Rather the effects compound each other.Existing methods based on summation of fatigue damage and creep damage ignore the interaction between fatigue and creep, and thus result in less reliable findings.Although the improved representations of creep and fatigue components have been proposed in the literature, such as the non-linear accumulated damage models for creep [7,8] and fatigue [9,10], the issue caused by ignoring the interacted effect of fatigue and creep is still not fundamentally solved.
In addition, the conventional formulations typically assume a power-law relationship between life and applied loading, as evident in the Basquin equation [11,12] and Coffin-Manson equation [13,14].Although this approach is simple, the coefficients need to be recalculated with changed temperature and/or frequency.Hence, this makes the design process inefficient and expensive because a large number of empirical data are required and must be re-fitted for each condition.To improve this limitation, others have attempted to introduce the variables of temperature and frequency into modified models, resulting in the Coffin-Manson-based creep-fatigue models proposed by Solomon [15], Shi [16], Jing [17] and Wong & Mai [18], and the Basquin-based creep-fatigue models developed by Kohout [19] and Mivehchi [20].However, these models may only be applied at the situations for which they were derived.They do not represent the creep-fatigue behaviour for other materials, hence the formulations Metals 2018, 8, 853 3 of 32 cannot present a unified characteristic.Furthermore, these models are determined by curve fitting, the accuracy of which is strongly determined by the number of empirical datum points.This results in poor empirical economy.Consequently, the conventional-based creep-fatigue models are severely limited in their applicability to engineering design.

Models Based on Observation of Microstructural Damage (Mechanism-Based)
The curve-fitting method, which is applied to build the Basquin-based and the Coffin-Manson-based models (Section 1.2), provides the simplest process to construct a numerical model, and thus is well-accepted in the field of mechanical engineering.However, from the perspective of material science, the fatigue models ideally should be constructed through observations of physical phenomena (such as the crack growth, diffusion creep, and void growth).This approach has resulted in the development of several mechanism-based models.These models are variously based on micromechanical cyclic void growth [21], partition of energy and micro-crack growth [22], and multistage fatigue theory [23].
These models are attractive because they relate physically measureable microstructural properties to life or total damage.Some of these models already include the ability to accommodate multiple forms of damage (including creep, fatigue, or oxidation), and represent both creep and fatigue in one formulation.However, this class of models suffers from limitations from the perspective of an engineering designer:

•
They relate to life evolution in some way, but often not in ways that are accessible to engineering design.This is a particular limitation of the damage models.

•
The mechanism-based models need to be validated.They require the measurement of microstructural parameters of damage.This information is not readily available to design engineers, certainly not at the onset of design.Also, designers do not select materials based on microstructure, but rather on mechanical properties.Furthermore, microstructural data are also not easily available during the service life of the part without resorting to destructive testing.

•
They have abstract mathematical formulations that are not easy to conceptualise, and are difficult to apply to design.

•
They typically have multiple coefficients in power law formulations, and each equation has sub-coefficients that can only be determined empirically by fitting.

•
They are not convenient for mechanical design.For example, for the material selection at the initial stage of design, it is not easy to investigate and determine the microstructural damage caused by fatigue, creep and oxidation.It is also not reasonable to assume multiple materials have the same damage.However, for the empirical-based models, the fatigue evaluation can be conveniently calculated through inputting the temperature, frequency, and applied loading.
From the perspective of mechanical engineering design, it is desirable that a creep-fatigue model should have a clear structure that is understood by engineering practitioners, include the general variables at the engineering scale (such as temperature, time, and loading), include parameters that are measureable or knowable, and be easily mathematically solved.This is not the case for the mechanism-based models.Furthermore, material standards are invariably based on assurances of mechanical properties and element composition, not on microstructural properties.Hence, while designers may be interested in microstructure, they cannot rely on in their specifications.

Extension of the Empirical Models
As mentioned in Section 1.2, the damage-accumulation-based models ignore the interaction between fatigue and creep.The microstructural interactions between creep and fatigue are beginning to be understood at a qualitative level, e.g., [24].Various mathematical expressions for this interaction are also available, with some (albeit limited) basis in microstructure or loading partitioning, e.g., [25,26].Hence a possible way to move the field forward from a materials and design perspective is to further improve the conventional Coffin-Manson-based creep-fatigue models.
A recent development in that direction has been the development of a model that includes temperature, cyclic time, applied loading, and with applicability to multiple (metallic) materials [27].This 'unified' model takes the form of a mathematical representation of plastic strain with functions including empirically determined coefficients: where ε p is the plastic strain which reflects fatigue capacity, N f is the creep-fatigue life, C 0 and β 0 are the fatigue ductility coefficient and fatigue ductility exponent respectively, which are related to fatigue capacity at the pure-fatigue condition, T is the temperature, t c is the cyclic time which presents the reciprocal of loading frequency, c 1 (σ) is the stress moderating equation which reflects the creep effect caused by the applied loading, c 2 is the constant, and σ reflects the applied loading which can be related to plastic strain through the cyclic strain-stress relation.
The equation also includes the concept of a reference condition.Here T re f is the reference temperature, which is defined as 35% of the melting temperature, t re f is the reference cyclic time which is suggested as a small value of 1 s.
The limitations presented by the existing Coffin-Manson-based models are improved by this model.The improvements are that: the structure includes the parameters of typical engineering problems, is easily mathematically solved, may be applied in multiple situations on multiple metallic materials, and covers the full range of conditions from pure fatigue to creep fatigue and then to pure creep.In particular, the model provides a more economic method for fatigue-life prediction since less empirical data are required than other empirical methods such as [15,17].In addition, the model is applicable for engineering design at the initial stage through combining with finite element analysis (FEA) [28].
Nonetheless from an engineering design perspective, the model has room for improvement.There is a need to have a representation that can predict fatigue life at a given applied loading (or can be used to evaluate the critical value of applied loading under a given life).This process of engineering calculation is applied at the early stages of engineering design, when candidate materials are being considered in relation to the functional requirements.Furthermore, it is necessary to represent the full range of fatigue, creep-fatigue, and creep conditions.From a design perspective it is essential that any model is able to be applied using the type of information available to a design engineer (which may be tentative or incomplete).

Opportunities for Modifying the Unified Model
There is something of philosophical debate between proponents of the mechanism-based models, and the empirical models.From the perspective of the mechanism-based models, design ought to be conducted by detailed examination of microstructure and the determination of multiple material parameters, some based on properties of the crystal lattice, defect sizes, oxidation factors, and curve-fitting parameters.The methods are valuable because they can relate say critical crack length to life.However, they have other limitations as described above.From the perspective of the empirical models, design ought to be conducted by performing macroscopic tests (no microstructural tests required) at various environmental conditions, and then curve-fitting to obtain coefficients for a formulation.The methods are valuable because they can be highly accurate, and they readily relate loading to life.However, they have other limitations as described above.
Both methods have strengths and weaknesses.Proponents of the different schools of thought tend not appreciate the approach taken by the other, which is strange since both rely on fitting of many coefficients, and formulations encapsulating many assumptions.In the longer term the mechanism-based models may prove to be superior, if they can eventually link the physical features of the virgin and damaged microstructure to life, using parameters that are easy to measure and available at design time.However, at present, the empirical models are superior, at least for engineering design purposes.Hence the further improvement of the existing models is still worthwhile attempting from an engineering perspective.
The development of the unified creep-fatigue model [27] was based on an assumption, which is the change rate (β 0 ) of applied loading to fatigue life is constant for different temperatures and cyclic times.Graphically, the curves of applied loading vs. fatigue life at the situations with different temperatures and cyclic times at the log-log coordinate are parallel.The model applied this assumption because the slopes of loading-life curves at the log-log coordinate change only slightly among the situations with different temperatures and cyclic times.Although the accuracy of fatigue-life prediction is acceptable [27], this assumption still suggests some opportunities for future improvements.
Firstly, the accuracy of the fatigue-life prediction could be further improved.Specifically, although the influence of temperature and cyclic time to the exponent (β 0 ) is slight, this influence may not be negligible.However, in the unified model (Equation ( 2)), this exponent is a constant, not a function of temperature and cyclic time.This implies that inclusion of this influence may improve accuracy of the fatigue-life prediction.In addition, the derivation of creep-fatigue-related coefficients was conducted by applying numerical optimization.This is a curve-fitting-based method, and thus the fitting quality strongly depends on the number of power series and coefficients.
Such methods generally benefit, as regards fitting accuracy, from provision of higher power series and more tunable coefficients.There are examples in the literature that specialize in this approach, and result in exceptionally good fits [15][16][17].However, this comes with two significant costs: (a) parameter non-identification becomes problematic in that multiple different combinations of parameters give similar results, hence the model becomes degenerate, and (b) it becomes difficult, even impossible, to link the coefficients to any meaningful parameters of physical properties or microstructure, hence the ontological power is depleted.Therefore, it is prudent to exercise restraint when expanding the terms within predictive models.It is preferable to add parameters that have some basis in physical reality.Consequently, we propose that the unified model might be improved by introducing new parameters for temperature and cyclic time into the exponent component (β 0 ) (See Sections 4.3 and 5.2).
Secondly, the description of the pure-fatigue condition could be further improved.Specifically, the unified model can be restored to the Coffin-Manson equation at the pure-fatigue condition which is described by the coefficients of C 0 and β 0 .These two coefficients are derived from the empirical data by numerical optimization.As mentioned above, the assumption may impact the accuracy of these two coefficients, and thus the quality of pure-fatigue description may be reduced.In this case, the modification for exponent component (β 0 ) may improve the accuracy of C 0 and β 0 , and then a better description for pure fatigue might be obtained (see Section 4.4).
In summary, we propose that the unified creep-fatigue model [27] could potentially be further improved through introducing the parameters of temperature and cyclic time into the exponent component.This has the potential to improve the accuracy of the fatigue-life prediction for both creep-fatigue and pure-fatigue conditions.
In the present work, we propose an explicit creep-fatigue model.

Methodology
The present work aims to further improve the unified creep-fatigue model [27].This new explicit model should present improved accuracy of the fatigue-life prediction for both the creep-fatigue and pure-fatigue conditions.We are also mindful of the need to make such models accessible for engineering design.This has not always been a strong feature of models in the literature.This requires consideration of the type of information available to designers, and an understanding of what they are trying to achieve.
To improve the model, we removed the assumption that β 0 in Equation ( 2) is constant, and then introduced the parameters of temperature and cyclic time into the exponent component.We retained from [27] the concept that the fatigue capacity is reduced due to active creep behavior, which is influenced by temperature and time [1,4].These two elements were included into the unified model (Equation ( 2)) through introducing a creep moderating function to the c component in Equation ( 2) [27].In the present work the additional change is the introduction of an additional creep moderating function (a function of temperature and cyclic time) to modify the fatigue ductility exponent (β 0 ).The numerical relationships among temperature, cyclic time and exponent component were extracted from the understanding of creep behaviour (diffusion creep).Then, to build a bridge between pure fatigue and creep fatigue, the reference condition was also introduced.By this way, the exponent component can be restored to β 0 at the pure-fatigue condition.
Creep mechanisms are normally divided into Nabarro-Herring creep, Coble creep, grain boundary sliding and dislocation creep [1,29].Nabarro-Herring creep and Coble creep show a strong dependency on temperature, where the diffusional flow of atoms occurs under conditions of relatively high temperature.Grain boundary sliding involves displacements of grains against each other.This is a particularly important mechanism for the creep failure of ceramics at high temperature because of the glassy-phase formation which provides a good sliding condition along the grain boundary.Dislocation creep presents progressive disruption through the crystal lattice, which results from both line defects and point defects, and can occur at relatively low temperature.This process is sensitive to the applied stress on the material, with a secondary dependency on temperature [1].
Based on the brief description of these four creep mechanisms, the diffusion creep (including Nabarro-Herring creep and Coble creep), which has strong temperature dependence, is used to extract the creep effect.(In the Discussion we briefly comment on the effect of ignoring these other creep mechanisms).
Then, an explicit creep-fatigue model was developed, see Section 3.1.This model was then validated on the materials of 63Sn37Pb solder and stainless steel 316 (see Sections 4.1 and 4.2).The coefficients were determined by the empirical data (including pure-creep data and creep-fatigue data) which were extracted from the literature.Ideally, the creep-fatigue data applied to obtain the coefficients and applied to validate the model should be extracted from two different literature sources.However, in the present work the empirical data are limited so, we extracted the empirical data from one source in the literature, and then the data were divided into two groups.One group was used to extract the coefficients of this model, and the other group used to validate this model.Hence if the experiments are conducted by following the experimental standard, the data at one specific condition (temperature, loading and cyclic time) should not be impacted by the location and operator.
After this, this model was compared with the unified and other models to evaluate the accuracy of the fatigue-life prediction at creep fatigue and pure fatigue (see Sections 4.3, 4.4 and 5.2), and the economy (see Section 5.3).In this process, the accuracy of life prediction for both the creep-fatigue and pure-fatigue conditions is discussed by evaluating the average error and prediction ratio (which are defined in Sections 3.2 and 3.3).In addition, the unified and integrated characteristics of the explicit model were investigated.Although the explicit model presents better the fatigue-life prediction, introducing more parameters into a numerical representation may result in poor economy for engineering designers because more empirical data may be required.Specifically, for an economical method for creep-fatigue-life prediction, the coefficients of this model should be obtained by fewer creep-fatigue experiments, because conducting creep-fatigue test is an expensive and time-consuming process.Thus, reduced empirical effect means better economy.This potential issue of economy is discussed (see Section 5.3).
Finally, the explicit model was applied to engineering design calculation (see Appendix A.1) and finite element analysis (see Appendix A.2). We provide specific directions for how the model may be used under both approaches, and the limitations thereof.
This new explicit model was developed with engineering design in mind.In particular, the general variables at the engineering scale (such as temperature, time and loading) were introduced to this explicit model, but the variables at the microstructural level (such as crack configuration, damage size, inter-void spacing, and oxidation) were not included.Although the explicit model still relies on empirical data, it is not a purely curve-fitting-based model.Specifically, the relationships between the different variables were derived from the understanding of creep and fatigue behaviours at the microstructural level, and the formulation was constructed by harmoniously integrating these relationships.This is not simply a curve-fitting-based process, thus gives an improved method for life prediction.During the process of engineering design, the coefficients are determined from the empirical data.

The Explicit Creep-Fatigue Model
We introduce the parameters of temperature and cyclic time into the exponent component.The modification process is presented in Section 3.1, and the method of determining the coefficients is presented in Section 3.2.

Development of the Explicit Creep-Fatigue Model
As mentioned in Section 1.4, the previous research applied an assumption that the fatigue ductility exponent (β 0 ) is constant at different temperatures and cyclic times.Removing this assumption gives an opportunity to further improve the unified model (Equation ( 2)).The unified model aimed to be applied for engineering design, thus the general variables at the engineering scale (temperature, time and applied loading) were included.However, at one specific temperature and cyclic time, applied loading does not influence the slope of life-loading curve, thus this parameter is not included into the exponent component and only the variables of temperature and cyclic time are included.In addition, according to the concept of fatigue capacity presented in [27], the slopes of life-loading curves gradually trend to zero with an increased creep effect (elevated temperature and prolonged cyclic time).
To resolve these issues, we introduce a creep moderating function b(T, t c ) to modify the fatigue ductility exponent, and then is further expanded as the form of '1 − x': We assume that time and temperature are not convoluted with each other, and thus the overall effect caused by temperature and time are additive.Later we show that this assumption gives sufficiently accurate outcomes.Then, function b (T, t c ) is split into a thermal component and a time component: Then, we determine the relationships among temperature, cyclic time and exponent component.This is achieved through investigating creep behaviour.Specifically, function β 0 b(T, t c ) implies the rate of fatigue-capacity decreases or increases between different temperatures and/or cyclic times.This rate can be described by diffusion-creep rate, and described by Fick's law [30]: Metals 2018, 8, 853 where J is the diffusion flux which shows that the amount of substance flowing through a unit area at a unit time (thus reflects the diffusion rate), D is the diffusion coefficient, x is the position and ϕ reflects the concentration of vacancies.The diffusion process is identified as a thermodynamic system due to the strong driving force of temperature.In this system, the transfer of atoms and the formation of vacancies are numerically evaluated by free energy at atomic level [31], and then the equilibrium atomic fraction of vacancies (N v ) is given by Equation ( 7): where ∆G f is the Gibbs free energy for formation of a vacancy, k is the Boltzmann's constant and T is the temperature.In Equation ( 6), ϕ is defined as the number of vacancies per unit volume, and thus is related to the atomic fraction by Equation ( 8): where Ω is the atomic volume.Therefore, a natural exponential relation between the diffusion flux and the temperature component can be presented: The expression of exp(1/T) can be simplified to a linear dependence when the temperature is relatively high enough, which is higher than the temperature where the creep behavior is activated (normally 0.35 of melting temperature), and usually the case when creep-fatigue is being considered in an engineering application.This provides a linear relationship, but the coefficient of the temperature (the slope of this straight line) should be determined from the empirical data.Thus, a linear relationship between diffusion rate and temperature arises: In addition, Equation (9) shows that the diffusion-creep behaviour gives a logarithmical relationship between temperature and diffusion flux.The definition of 'diffusion flux' indicates that this term measures the amount of substance flowing through a cross sectional area during a unit time.Thus, a time dependence is included in this parameter in the form of a rate function.Then, Equation ( 9) can be presented as: where D v reflects the amount of substance flowing through a unit area.This equation gives a logarithmical relation between temperature and cyclic time: Then, the linear relationship of temperature vs. exponent component and the logarithmical relationship of temperature vs. cyclic time are integrated into Equation (5).The moderating function b(T, t c ) is presented by Equation (13): where b 1 and b 2 are constant and determined by empirical data.To build a bridge between pure fatigue and creep fatigue, we introduce the thermal and cycle time reference condition into Equation (13), then this equation is modified as: Finally, the explicit creep fatigue model is given as:

The Method of Determining the Coefficients
The coefficients of the explicit model (Equation ( 15)) are determined by the empirical data, including pure-creep data and creep-fatigue data.

Selecting the Reference Condition
The creep damage is assumed to be active above the reference temperature and the reference cyclic time.The reference temperature is defined as 35% of the melting temperature [32], and the reference cyclic time is suggested as a small value, nominally 1 s.

Deriving the Coefficients of Function c(σ, T, t c )
The method to derive the coefficients of c(σ, T, t c ) proposed in [27] is extended to the present work.In this case, function c 1 (σ) and constant c 2 are presented by Equations ( 17) and ( 18): In Equation ( 17), P MH (σ) is a function which represents the relationship between the Manson-Haferd parameter and applied stress (σ).The Manson-Haferd parameter under one specific stress is numerically presented as: where T is the absolute temperature, t is the creep-rupture time, and (log t a , T a ) is the point of convergence of the log t-T lines.In particular, T a is defined as the reference time (T re f ) below which creep is dormant.Both Equations ( 17) and ( 18) are obtained by the empirical data of pure creep.Specifically, during creep-rupture tests, the temperatures (T), stresses (σ) and creep-rupture times (t) are recorded.Then, the relationships between T and log t under different stresses are plotted (Figure 1), wherein the temperature at the point of convergence is identified as the reference temperature, and the value of log t at this convergence point (log t a ) is given by the average value of the logt(T re f ) at different stresses.The value of log t then gives c 2 .
Both Equations ( 17) and ( 18) are obtained by the empirical data of pure creep.Specifically, during creep-rupture tests, the temperatures (T), stresses (σ) and creep-rupture times (t) are recorded.Then, the relationships between T and log  under different stresses are plotted (Figure 1), wherein the temperature at the point of convergence is identified as the reference temperature, and the value of log  at this convergence point (log  ) is given by the average value of the logt( ) at different stresses.The value of log  then gives  .According to Figure 1, the Manson-Haferd parameters under different stresses are given, then the relationship between the Manson-Haferd parameter and applied stress (P MH (σ)) can be obtained through curve fitting.Then, function c 1 (σ) is expanded.

Deriving the Coefficients
The remaining coefficients in the explicit model are determined by the empirical data of creep fatigue.Specifically, during the creep-fatigue tests, the temperatures (T), cyclic time (t), stresses (σ), plastic strain (ε P ) and fatigue life (N) are recorded.In particular, with the empirical data of plastic strain vs. stress, the coefficients (K and n ) of the cyclic strain-stress relation under different temperature-cyclic time conditions are obtained.In the present work, these two coefficients are applied to describe the engineering quantities-based relationship, and a power-law-based transition between strain and stress is included.They then are involved in the function c 1 (σ) for transforming stress into plastic strain (Equation ( 20)), and a moderating factor ( f m ) is introduced to compress the stress effect on creep-related damage.We did not separate the whole of applied loading (σ) into two components.This is because we cannot say one part the applied loading contributes to creep, and another part contributes to fatigue.Therefore, we defined that the whole of applied loading works for both fatigue and creep damage.
In the present work, we define f m as a stress-moderating factor which is applied to compress the cyclic stress to an equivalent constant stress.This moderating factor is related to the shape of the loading wave, and presents the average level of the cyclic loading.Illustratively, the area below the contour of the cyclic loading along the time dimension should be equal to the area below the contour of the equivalent constant loading at the same time period.This is based on an assumption that creep makes the same contribution to tensile and compressional portions.Although this assumption may be not appropriate for some materials [33,34], it simplifies the method of extracting this factor.For example, f m is defined as 0.6366 for the sinusoidal wave and as 0.5 for the triangular wave.
Then, numerical optimisation was applied to derive the coefficients of C 0 , β 0 , b 1 and b 2 by minimizing the average difference (δ a ) between the predicted fatigue life (N pre,ij ) and the experimental results (N exp,ij ) (Equation ( 21)).
Metals 2018, 8, where n is the number of data, and N ij presents the fatigue life obtained at multiple conditions of (T, t c ) j and strain amplitude i.

Evaluation of the Explicit Model
The quality of fatigue-life prediction is evaluated by the prediction ratio.Specifically, the prediction ratio (Equation ( 22)) gives the ratio of predicted creep-fatigue life to experimental creep-fatigue life: In the present work, we define that: An acceptable prediction ratio should be between 0.75 and 1.25.This range is narrower (more conservative than the range shown in other literature, wherein a factor of 2 or 1.5 is normally given [25,35,36].This also can be shown illustratively, where all data points of N pre vs. N exp under multiple temperatures and cyclic times should fall between the upper bound (+25%) and the lower bound (−25%) relative to ideal correlation (H = 1).

Validation
The explicit model is validated on the materials of 63Sn37Pb solder and stainless steel 316.The coefficients are determined by using the method proposed in Section 3.2, where the empirical data are extracted from the literature.The quality of fatigue-life prediction is evaluated by the method proposed in Section 3.3.

Deriving the Coefficients
The reference temperature for 63Sn37Pb solder was chosen as 160 K and the reference cyclic time was defined as 1 s.The creep-rupture data [37] are plotted in Figure 2, and the point of convergence (T re f , log t a ) is evaluated as (160 K, 8.232).This gives = 0.1215 (23) and the relationship between stress and the Manson-Haferd Parameter: Then, substituting into Equation ( 24), function c 1 (σ) is expressed as: and the magnitude of f m is given as 0.6366 for the sinusoidal wave.
•  •  and the magnitude of  is given as 0.6366 for the sinusoidal wave.
To evaluate the explicit creep-fatigue model, another groups of creep-fatigue data (Table 3) [16] are used to compare with predicted fatigue life which is supported by the results shown in Section 4.1.1.The prediction ratio (N pre /N exp ) under multiple temperatures and cyclic times are plotted in Figure 3, where all data points fall between the upper bound (+25%) and the lower bound (−25%).The upper bound and the lower bound present the prediction ratios are 0.75 and 1.25 respectively.This implies that the explicit creep-fatigue equation provides a high quality of fatigue-life prediction, specifically, a relatively high correlation between predicted and experimental creep-fatigue life.

Deriving the Coefficients
The reference temperature for stainless steel 316 was chosen as 585 K and the reference cyclic time was defined as 1 s.The creep-rupture data [38] are plotted in Figure 4, and the point of convergence (T ref , logt a ) is evaluated as (585 K, 10.783).This gives = 0.09274 (26) and the relationship between stress and the Manson-Haferd Parameter: Then, substituting into Equation ( 27), function c 1 (σ) is expressed as: and the magnitude of f m is given as 0.5 for the triangular wave.The creep-fatigue coefficients [39] obtained from the literature are tabulated in Table 4. Minimizing the difference between the predicted creep-fatigue life (N pre,ij ) and the experimental creep-fatigue life (N exp,ij ) yields C 0 = 7.768, β 0 = 0.571, b 1 = −0.000225and b 2 = −0.0223,and returns an average error (δ a ) (Equation ( 21)) of 0.00255.

𝑃 σ
and the magnitude of  is given as 0.5 for the triangular wave.
The creep-fatigue coefficients [39] obtained from the literature are tabulated in Table 4. Minimizing the difference between the predicted creep-fatigue life ( , ) and the experimental creep-fatigue life ( , ) yields  = 7.768,  = 0.571,  = −0.000225and  = −0.0223,and returns an average error ( ) (Equation ( 21)) of 0.00255.Consequently, the coefficients of the explicit creep-fatigue equation for stainless steel 316 are collected in Table 5: Table 5.The coefficients of the explicit formulation for stainless steel 316.

.2. Evaluation
To evaluate the explicit creep-fatigue model, another groups of creep-fatigue data (Table 6) [39] are used to compare with predicted fatigue life which is supported by the results shown in Section 4.2.1.
The prediction ratio (N pre /N exp ) under multiple temperatures and cyclic times are plotted in Figure 5, where all data points fall between the upper bound (+25%) and the lower bound (−25%).The upper bound and the lower bound present the prediction ratios are 0.75 and 1.25 respectively.This implies that the explicit creep-fatigue equation provides a high quality of fatigue-life prediction, specifically, a relatively high correlation between predicted and experimental creep-fatigue life.

Accuracy Comparison: Explicit vs. Unified Models
The present work aims to improve the accuracy of the creep-fatigue-life prediction through further modifying the unified model.Thus, the ability of life prediction by applying the explicit model should be better than applying the unified model.This is proved through comparing the explicit model with the unified model on the materials of 63Sn37Pb and stainless steel 316.
Specifically, the explicit formulation removes the assumption applied in the unified creep-fatigue model (Equation ( 2)), and then a creep moderating function was introduced into the exponent component.In this way, the explicit model has better ability to describe creep fatigue.To prove this, we applied the creep-fatigue data (Table 1 for 63Sn37Pb solder and Table 3 for stainless steel 316) to extract the coefficients of the unified model (Equation ( 2)) and the explicit model (Equation ( 15)).Then, these coefficients were applied to predict the fatigue life for the situations shown in Table 4 for 63Sn37Pb solder and Table 6 for stainless steel 316.The empirical data, and predicted life given by the unified model and the explicit model are illustrated in Figure 6 for 63Sn37Pb and Figure 7 for stainless steel 316.Table 7.The average errors and prediction ratios given by the unified model and the explicit model.Figures 7 and 8 show that life-loading curves given by the explicit model are closer to the empirical data, thus we conclude that the explicit model has better ability to predict fatigue life at the creep-fatigue condition.This is also proved by the average errors and prediction ratios (see Equations ( 21) and ( 22) in Sections 3.2.3 and 3.3 for the definitions of the average error and prediction ratio respectively) in Table 7.In particular, the value of the prediction ratio in Table 7 is represented by a range which is given by the maximum and minimum prediction ratios of the whole results.This representation is also shown in Table 8.

General Process of Validation for Other Materials
The explicit model could be further validated through involving more empirical data on more materials.The general process of validation can be summarised as follows: (1) Obtain the empirical data for one specific material.The data include pure-creep data and creep-fatigue data, which could be extracted from the literature, or collected by performing testing.In particular, the creep-fatigue data under multiple temperatures and cyclic times are divided into  Table 7 shows that the explicit model provides smaller average errors and narrower ranges of prediction ratio for both the materials of 63Sn37Pb solder and stainless steel 316.This demonstrates that the explicit model has better accuracy for quantitatively representing creep fatigue.present work, this ability was evaluated through comparing the predicted life with the empirical data on stainless steel 316 [39] (Figure 8).
Figure 8 shows that the loading-life curve formulated by the coefficients of C 0 and β 0 in the explicit model is closer than the unified model to the empirical data.This is also described by the prediction ratio.The prediction ratios for these two models is presented in Figure 9.The dotted lines (bounds) which are labeled by 1, −25%, −50% and −75% in Figure 9 represent the prediction ratios of 1, 0.75, 0.5 and 0.25 respectively.

General Process of Validation for Other Materials
The explicit model could be further validated through involving more empirical data on more materials.The general process of validation can be summarised as follows: (1) Obtain the empirical data for one specific material.The data include pure-creep data and creep-fatigue data, which could be extracted from the literature, or collected by performing testing.In particular, the creep-fatigue data under multiple temperatures and cyclic times are divided into two groups (3 to 4 sub-group data at different temperatures and cyclic times for each group).One group data (named Group1) are applied to determine the coefficients of the explicit model, and other group data (named Group2) are applied to compare with the predicted life (evaluate accuracy of the fatigue-life prediction).
(2) Determine the coefficients of the explicit model.The coefficients are determined by the method presented in Section 3.2.Figure 9 shows that the pure-fatigue prediction ratios are around 0.75 for applying the coefficients of the explicit model, but the prediction ratios are lower, between 0.5 and 0.25, for using the coefficients of the unified model.
Both Figures 8 and 9 imply that the errors between the empirical data and the predicted life given by the explicit model are smaller.Therefore, we conclude that the explicit model has better ability to describe pure fatigue than the unified model.

General Process of Validation for Other Materials
The explicit model could be further validated through involving more empirical data on more materials.The general process of validation can be summarised as follows: (1) Obtain the empirical data for one specific material.The data include pure-creep data and creep-fatigue data, which could be extracted from the literature, or collected by performing testing.In particular, the creep-fatigue data under multiple temperatures and cyclic times are divided into two groups (3 to 4 sub-group data at different temperatures and cyclic times for each group).One group data (named Group1) are applied to determine the coefficients of the explicit model, and other group data (named Group2) are applied to compare with the predicted life (evaluate accuracy of the fatigue-life prediction).
(2) Determine the coefficients of the explicit model.The coefficients are determined by the method presented in Section 3.2.
(3) Predict fatigue life.With the coefficients obtained in step 2, the predicted life under the situations presented in 'Group2' are calculated through using the explicit model.
(4) Evaluate the explicit model.The evaluation of the explicit model is conducted by the method given in Section 3.3.This process is numerically and illustratively presented by the prediction ratios.If the predicted data satisfy the range of acceptation given in Section 3.3, we can conclude that the explicit model can be applied on this material to predict creep-fatigue life.

The Characteristics of the Explicit Creep-Fatigue Model
The explicit creep-fatigue model was validated on the materials of 63Sn37Pb solder and stainless steel 316 (see Section 4).This implies that this model has ability to be applied at multiple temperatures and cyclic times, and the relationships between different variables (temperature, cyclic time, applied loading and life) in the explicit model are applicable for different materials.
In addition, at the reference (the pure-fatigue) condition (where T = T re f and t c = t re f ), the explicit creep-fatigue model can be restored to the Coffin-Manson equation.At the pure-creep condition (where c(σ, T, t c ) = 0), the explicit creep-fatigue model can be reformed as the Manson-Haferd parameter for creep.Consequently, the explicit formulation recovers both of the standard fatigue and creep formulations.
Provides one formulation that covers the full range of conditions from pure fatigue, to creep fatigue, then to pure creep.

2.
Recovers the mathematical formulation of both of the standard fatigue and creep formulations (Coffin-Manson and Manson-Haferd respectively).

3.
Accommodates multiple temperatures.Specifically, the explicit model can be applied to predict fatigue life at situations with different temperatures.

4.
Accommodates multiple cyclic times.The explicit model is applied at the cyclic loading without hold time, thus the cyclic time refers to the period of one cycle of this loading condition.This is a limitation of this model, which will be discussed in Section 5.4.

5.
Accommodates multiple materials.The explicit model was not a purely empirical-based model because the physical meaning was indirectly introduced into the explicit model.This process is quite different from the curve-fitting method.Thus, we conclude that the explicit model is potential able to be applied for multiple materials: we have demonstrated validation for 63Sn37Pb solder and stainless steel 316.Further validation on different materials is needed: this will be discussed in Section 5.4.

6.
Provides a physical basis for the structure of the formulation.The basis of the c term has been explained previously [40].Specifically, diffusion-creep behaviour gives a linear relationship of temperature vs. loading and a logarithmical relationship of temperature vs. cyclic time.Plastic zone around the crack tip gives a power-law relation between life and loading.The new b term is justified on principles of diffusion-creep rate and represents Fick's law (see Section 3.1).Both c and b terms were built on the concept of fatigue capacity, which was formulated as '1 − x'.This formulation numerically presents the negative effect of creep on fatigue.In addition, the introduction of the reference condition gives an opportunity to connect pure fatigue with creep fatigue.
Attributes 1-2 may be considered 'integrated' attributes, 3-5 'unified' attributes, and 6 a 'natural origin' attribute.Regarding integration, the models based on microstructural features, e.g., the integrated creep-fatigue theory [41,42], also offer an integrated characteristic.However, the determination of microstructural variables (such as crack, damage size and inter-void spacing) is challenging in the engineering situation.The explicit model is potentially easier to use in the engineer case.
Furthermore, the explicit model can be restored to the Coffin-Mason equation at pure fatigue and can be reformed as the Manson-Hefard parameter at pure creep.Both of these two formulations are conventional engineering models, and can be applied for engineering design without the need of observations at the microstructural level.This is a positive feature.
By natural origin we do not necessarily mean that the model has a physical basis traceable to microstructure and mechanical properties.Rather that the formulation of the model is consistent with existing representations of principles of physics (e.g., laws).We acknowledge that a full connection of all parameters in the explicit model to measureable variables of microstructure remains elusive.This limitation applies to all creep and fatigue models.
While there are other creep-fatigue formulations that also have high accuracy, they lack one or more of the features of the explicit model: they do not have the integrated characteristic; they are typically accurate only for specific cases (poor unified attribute); or they rely on the inclusion of many coefficients (typically into power series) which have no natural origin.Many of the competing models are so over-endowed with coefficients, e.g., [16,18], that they also have the risk of parameter non-identifiability.

The Ability to Describe Creep Fatigue
The unified model (Equation ( 2)) presents better ability to predict life at the creep-fatigue condition.This was proved through comparing the unified model with the existing creep-fatigue models.For example, the unified model was compared with Solomon's model [15], Jing's model [17], and Wong & Mai's model [18].Both Solomon's model and Jing's model use fixed coefficients.When they are applied to other situations, Solomon's model results in a poor average error (23.96), and Jing's model cannot give any numerical solution.Thus, they only can be used in the situations where they were derived, and cannot be extended to other situations and other materials where there are no empirical data.This is because these models determine their coefficients by numerical optimisation across all variables (including temperature, frequency, fatigue life and applied loading).Hence when changing to a different material it is necessary to recalculate all the coefficients: it is not possible to simply change only some of the coefficients.However, Wong & Mai's model has potential to be applied to multiple materials.This is because this model has seven independent coefficients which are required to be recalculated for different materials.The accuracy of these models comes at the cost of high specificity of the coefficients, and the risk of parameter non-identifiability.Also, the coefficients in the power series terms have no physical identity, but only exist to provide improve mathematical fit.
To allow a comparison with the explicit model, we re-calculated the coefficients for Solomon's model, Jing's model, and Wong & Mai's model, as follows.
Solomon's model (Equation ( 29)) is: where T is temperature in • C, f is the frequency, N f is the fatigue life, and c 1 , c 2 , c 3 , c 4 , k and β 0 are constants derived from the empirical data.Jing's model (Equation ( 31)) is: with Wong & Mai's model (Equation ( 33)) is: with where n is cyclic hardening index, σ yield is the yield stress, T is the temperature in Kelvin, f is the frequency, T re f is the reference temperature below which creep becomes dormant, f re f is the reference frequency, and C 0 , β 0 , A , c 1 , c 2 , b 1 and b 2 are constant.We recalculated all these constants, for these three models, for stainless steel 316, and plotted the results in Figure 10 to compare with the explicit model and empirical data.Both Figures 10 and 11 show that the explicit model (Equation ( 15)) has better numerical accuracy for describing creep fatigue than the models of Solomon, Jing, and Wong & Mai models.For Solomon's and Jing's models, the average error given by these two models are much higher than the explicit model.This is because the relationships between different variables in these two models were completely derived from the empirical data of one specific material, thus they cannot be extended to other materials.However, Wong & Mai's model presents better ability for life prediction than Solomon's and Jing's models (except the situation of 973 K-4%/min).This may be because Wong & Mai's model involves a material property (yield stress), and it also applied a concept that would later be referred to as fatigue capacity.However, this model has seven independent coefficients which are required to be determined by empirical data.Consequently, this leads to another issue, that of economy.

The Economy
The economy is an important factor which is considered during the process of engineering design [4].An economical method should provide a good balance between the accuracy and cost.Although the mechanism-based models may not need any creep-fatigue tests, observing and measuring microstructure is not a simple and economic process for engineering practitioners.The empirical-based models are more suited for engineering purposes, hence are the point of comparison for the explicit model.We selected Wong & Mai's equation [18] to compare with the explicit model regarding to the economy, since the Wong & Mai's equation shows better life-prediction ability than other existing models (see Section 5.2).
In the present work, empirical data under different temperatures and cyclic times was taken from the literature (Table 4).The first stage took seven groups of data and split this into a group of six and one.The six groups of empirical data were applied to derive the coefficients of Wong & Mai's equation and the explicit model.Then, these coefficients were used to predict fatigue life at the condition of the remaining data set (named 'predicted condition').The discrepancy was noted.
The second stage repeated this analysis but with three and four groups respectively, and again the discrepancy was noted.Finally, the average errors and prediction ratios obtained in these two situations were compared.Thus, it becomes possible to infer how sensitive each model is to the available quantity of data.A model with better economy would be one where the degradation in accuracy was less sensitive to the quantity of data.Both Figures 10 and 11 show that the explicit model (Equation ( 15)) has better numerical accuracy for describing creep fatigue than the models of Solomon, Jing, and Wong & Mai models.For Solomon's and Jing's models, the average error given by these two models are much higher than the explicit model.This is because the relationships between different variables in these two models were completely derived from the empirical data of one specific material, thus they cannot be extended to other materials.However, Wong & Mai's model presents better ability for life prediction than Solomon's and Jing's models (except the situation of 973 K-4%/min).This may be because Wong & Mai's model involves a material property (yield stress), and it also applied a concept that would later be referred to as fatigue capacity.However, this model has seven independent coefficients which are required to be determined by empirical data.Consequently, this leads to another issue, that of economy.

The Economy
The economy is an important factor which is considered during the process of engineering design [4].An economical method should provide a good balance between the accuracy and cost.Although the mechanism-based models may not need any creep-fatigue tests, observing and measuring microstructure is not a simple and economic process for engineering practitioners.The empirical-based models are more suited for engineering purposes, hence are the point of comparison for the explicit model.We selected Wong & Mai's equation [18] to compare with the explicit model regarding to the economy, since the Wong & Mai's equation shows better life-prediction ability than other existing models (see Section 5.2).
In the present work, empirical data under different temperatures and cyclic times was taken from the literature (Table 4).The first stage took seven groups of data and split this into a group of six and one.The six groups of empirical data were applied to derive the coefficients of Wong & Mai's equation and the explicit model.Then, these coefficients were used to predict fatigue life at the condition of the remaining data set (named 'predicted condition').The discrepancy was noted.
The second stage repeated this analysis but with three and four groups respectively, and again the discrepancy was noted.Finally, the average errors and prediction ratios obtained in these two situations were compared.Thus, it becomes possible to infer how sensitive each model is to the available quantity of data.A model with better economy would be one where the degradation in accuracy was less sensitive to the quantity of data.if the modification was conducted for the coefficient component, since this already has a several tunable coefficients.
Consequently, we conclude that the explicit creep-fatigue model presents an economical method for the fatigue-life prediction, and introducing two more coefficients into the exponent component does not significantly impact the economy.

Limitations and Implications for Future Research
Limitations that designers need to note are that strictly speaking the method has only been validated for the materials of 63Sn37Pb and stainless steel 316.However, this explicit model is potentially usable for other metallic materials.We anticipate difficulties applying this model for plastics and composites because they present totally different material characteristics and failure mechanisms.However, it is not impossible that this explicit model may be further improved and extended to other material categories as more empirical data are included.In particular, nylon is widely used in the engineering industry for load bearing parts.Thus, it may be an interesting future project to check and adapt this model to engineering nylon (such as nylon 6).
The explicit model is not yet ideal regarding natural origin as it does not include quantifiable microstructural properties.To achieve this, it would be necessary to better understand the microstructural processes of fatigue & creep-especially their interactions-and how those affect plastic strain and life.Some work is available in this area, e.g., [24], but there is still a long way to go before the values of coefficients in a creep-fatigue model can be predicted ab initio from microstructural inspection.In addition, as mentioned in Section 2, the explicit model ignores the dislocation creep and grain boundary sliding.In this case, introducing these two behaviours to reflect creep effect at high stress may be beneficial.
Another potential avenue of future research is to continue the process of extending existing models towards a more complete theory, as has been illustrated here with the redevelopment of the unified model into the explicit.During the process of development, more microstructural-level-based parameters may be included, with a corresponding inclusion of new terms into the model.We suggest that it is worthwhile designing these extensions to include other well-established phenomena, as we demonstrated in Section 3, rather than merely chasing better accuracy by adding more power terms and coefficients.
The situation of cyclic loading without hold time (dwell-fatigue) is not covered by the explicit model.In this loading condition, fatigue makes more of a contribution than creep, because the total time is too small to produce marked creep damage.However, for cyclic loading with hold time, the creep effect gradually intensifies as the hold time increases.Then more creep damage is produced than fatigue damage, and the failure finally occurs due to the creep effect.We could imagine that in the situation with a relatively short hold time, the explicit model may still present a reasonable prediction of fatigue life, but the accuracy of this prediction may become worse when the hold time is prolonged.This implies that the explicit formulation has an opportunity to be further improved to cover the situation with hold time or relatively long cyclic time.To achieve this, it would seem necessary to modify the formulation (especially, the creep component in this explicit model) to include new terms of as yet-unknown mathematical form.Conceptual works, e.g., [24], may be useful in identifying the basic form of these relationships.
At elevated temperature, the crack surface is oxidized, and then the material becomes more brittle.This results in further crack propagation.Therefore, the oxidation effect should ideally be included.The class of models based on observation of microstructure, e.g., [41,42], are superior in this regard because they can measure the voids and internal damage.
The class of models based on macroscopic empirical testing, to which the explicit model belongs, lack the microstructural parameters of crack length, oxidation, etc.At least not as primary variables, but the effects are partly accommodated through other means.In the explicit model, the coefficients are determined from the empirical data through numerical optimization, thus any oxidation effects are incorporated into the fitting process.Although the accuracy of fatigue-life prediction may be impacted, the results still show acceptable accuracy (see Sections 4.3 and 5.2).
In future work it might be possible to include oxidization in the explicit equation.Superficially this might involve simply including a power series term.However this may not be entirely successful, because our observation is that simply adding more terms and coefficients does improve accuracy, but at the cost of introducing model degeneracy.This has the further consequence of making the model more highly dependent on the specific situation, i.e., reduces the ability of the model to generalize to other materials and situations.The challenge is to include the oxidation effect, in a way that is coherent with how the effect operates physically, and to do so using parameters that are identifiable by the engineering designer.This opens an opportunity to further improve this engineering-based model.

Application to Engineering Design and Structural Mechanics
The present work aims to develop a creep-fatigue model for engineering design, thus this section is included to briefly explain how this model is applied by manual engineering calculations and finite element analysis at the engineering design process.
Fundamentally, the explicit creep-fatigue model can be used to predict fatigue life at a given applied loading, or can be used to evaluate the critical value of applied loading under a given life.This process of engineering calculation is normally applied at the initial stage of engineering design.For example, the explicit model can be applied to select a material.
In addition, the explicit model can also represent the pure fatigue condition.This is because it can be restored to the Coffin-Manson equation at the reference condition (T = T re f and t c = t re f ), which represents pure fatigue wherein the creep effect is dormant.In this case, this restored equation can be used to predict the fatigue life or critical value of applied loading at pure fatigue.The accuracy of pure-fatigue description was demonstrated in Section 5.3, which implies the coefficients of C 0 and β 0 obtained in creep fatigue can be extended to predict fatigue life at pure fatigue.
The method may be applied to manual calculation or finite element analysis, as shown in Appendix A.

Conflicts of Interest:
The authors declare no conflict of interest.The research was conducted without personal financial benefit from any funding body, and no such body influenced the execution of the work.This example is also summarised in Table A1.Table A1 shows that these two life-evaluation processes give similar results (fatigue life), with the explicit method being slightly more conservative.The explicit method is faster to implement even for a single pass through the design, and has further time advantages when there are revisions and loops in the design process.We conclude that combining the explicit model with FEA can reduce the difficulty and complexity of analysis regarding fatigue evaluation, and speed up the design process.These benefits are particularly attractive early design stages, when the design is still experimental and the loading conditions are not finalized.This example is also summarised in Table Table  Table A1 shows that these two life-evaluation processes give similar results (fatigue life), with the explicit method being slightly more conservative.The explicit method is faster to implement even for a single pass through the design, and has further time advantages when there are revisions and loops in the design process.We conclude that combining the explicit model with FEA can reduce the difficulty and complexity of analysis regarding fatigue evaluation, and speed up the design process.These benefits are particularly attractive early design stages, when the design is still experimental and the loading conditions are not finalized.

Figure 1 .
Figure 1.Plotting of T vs. log t.

Figure 6 .
Figure 6.Empirical data and predicted life for 63Sn37Pb.

Figure 7 .Table 7 .
Figure 7. Empirical data and predicted life for Stainless steel 316.

Figure 7 .
Figure 7. Empirical data and predicted life for Stainless steel 316.

Figure 7 .
Figure 7. Empirical data and predicted life for Stainless steel 316.

Figure 9 .
Figure 9. Prediction ratio for stainless steel 316 at pure fatigue.

Figure 8 .
Figure 8. Predicted life and empirical data of stainless steel 316 for pure fatigue.

4. 4 .
The Ability to Describe Pure Fatigue Both the unified creep-fatigue model and the explicit model can be restored into the Coffin-Manson equation at pure fatigue.This loading condition is numerically presented by the coefficients of C 0 and β 0 , thus the accuracy of pure-fatigue description is determined by them.In the Metals 2018, 8, 853 18 of 32

Figure 8 .
Predicted life and empirical data of stainless steel 316 for pure fatigue.

Figure 9 .
Figure 9. Prediction ratio for stainless steel 316 at pure fatigue.

Figure 9 .
Figure 9. Prediction ratio for stainless steel 316 at pure fatigue.

4 Figure 10 .
Figure 10.Empirical data for stainless steel 316, and predicted life given by the explicit model, Solomon's model, Jing's models, and Wong & Mai's model.Figure 10.Empirical data for stainless steel 316, and predicted life given by the explicit model, Solomon's model, Jing's models, and Wong & Mai's model.

Figure 10 .
Figure 10.Empirical data for stainless steel 316, and predicted life given by the explicit model, Solomon's model, Jing's models, and Wong & Mai's model.Figure 10.Empirical data for stainless steel 316, and predicted life given by the explicit model, Solomon's model, Jing's models, and Wong & Mai's model.
The average errors calculated by Equation(21), for the explicit model, Solomon's model, Jing's model, and Wong & Mai's model, for stainless steel 316, are plotted in Figure11:The average errors calculated by Equation (21), for the explicit model, Solomon's model, Jing's model, and Wong & Mai's model, for stainless steel 316, are plotted in Figure11:

Figure 11 .
Figure 11.Average error for predicted life given by the explicit model, Solomon's model, Jing's model, and Wong & Mai's model.

Figure 11 .
Figure 11.Average error for predicted life given by the explicit model, Solomon's model, Jing's model, and Wong & Mai's model.

A cross-sectional area A constant b 1 , b 2 , b 3 , b 4 and b 5 constants C 0 fatigue ductility coefficient c 1 , c 2 , c 3 , c 4 ,
c 5 and c 6 constants D diffusion coefficient D v amount of substance flowing through a

Figure A2 .
Figure A2.Fatigue life obtained at the creep-fatigue condition using conventional method.Figure A2.Fatigue life obtained at the creep-fatigue condition using conventional method.

Figure A2 .
Figure A2.Fatigue life obtained at the creep-fatigue condition using conventional method.Figure A2.Fatigue life obtained at the creep-fatigue condition using conventional method.

Figure A2 .
Figure A2.Fatigue life obtained at the creep-fatigue condition using conventional method.

Figure A3 .
Figure A3.Fatigue life obtained at the pure-fatigue condition using explicit method.

Figure A3 .
Figure A3.Fatigue life obtained at the pure-fatigue condition using explicit method.

Table 2 :Table 2 .
The coefficients of the explicit formulation for 63Sn37Pb solder.

Table 7 .
The average errors and prediction ratios given by the unified model and the explicit model.

Table 8 .
Accuracy of prediction regarding the empirical-data number.

Table A1 .
Fatigue evaluation by FEA.