University of Huddersfield Repository On the development of creep damage constitutive equations: modified hyperbolic sine law for minimum creep strain rate and stress and creep fracture criterion based on cavity area fraction along grain boundaries On the development of creep damage constitutive equations: a modified hyperbolic sine law for minimum creep strain rate and stress and creep fracture criteria based on cavity area fraction along grain boundaries On the development of creep damage constitutive equations: a modified hyperbolic sine law for minimum creep strain rate and stress and creep fracture criteria based on cavity area fraction along grain boundaries

This paper reports (1) the latest development and application of a modified hyperbolic sine law for minimum creep strain rate and stress for both low Cr and high Cr steels, and (2) the development of a creep fracture criterion based on cavity area fraction along grain boundaries for high Cr steel. This work is part of the fundamental development of creep damage constitutive equations which were identified through a critical literature review. In the former the application of the new law results in an improved fitting; in the latter, a new creep fracture criterion based on cavity area fraction along grain boundaries was derived and quantitatively calibrated using the latest detailed cavity nucleation and growth kinetics models for high Cr steel. Furthermore, this paper revealed the trend of nucleation rate coefficient with stress, and the trend of creep life time coefficient with stress, which provide reliable and universal prediction capabilities. This paper contributes to the specific knowledge on the minimum creep strain rate and stress function, the development of a scientifically sound and novel creep rupture criterion based on the cavity area fraction along grain boundaries for high Cr steel, and the provision of creep damage/life prediction tools.


Stress breakdown
In the development of high temperature materials for long term use, the accelerated creep test (short-term) may be used although it may not be reliable in predicting long time service. During high temperature service, this microstructure may progressively degrade, leading to unexpected consequences for long-term (100,000) creep performance.
The stress breakdown phenomenon for high Cr steel received positive reports in past literature [1] and was praised by other research papers and reviews. Basically, it can be summarized as (1) there is a change of mechanisms in creep deformation, fracture, and the coupling between creep deformation and creep damage; and (2) it is understood that the nature of the cavity nucleation, growth and coalesce changes with the stress level, resulting in a change from a ductile to brittle fracture.

Constitutive modelling of creep damage
In an attempt to model the creep behaviour of high Cr steel, Yin et al. [2] originally proposed a phenomenological relationship between the creep cavity damage and creep strain, which departed from the firm mechanism based relationship proposed by Dyson [3]. The relevant equations are listed below for further use: where f u is the unaixal strain at fracture, k N has an upper limit of ≈ 1/3; Such an approach cannot be extrapolated into lower stress levels beyond the range it has been calibrated [4].
Basirat [5] adopted Yin's approach, but relaxed the definition of material constant A, by allowing there to be a function of stress; the material constant A is a function of temperature and stress. The specific variation of cavitation coefficient A for P91 steel is shown in Table 1 [5] and graphically in Figure 1 [6].

OPEN ACCESS
Due to the lack of a clear trend for A, it is difficult to use it in prediction with confidence. The current creep damage modelling for long-term service is not yet satisfactory.

Challenges involved with creep damage modelling
The challenges involved with creep damage mechanics and creep damage modelling have been explored in a literature review [7], and are summarised below: (1) Characterising and quantifying the creep cavitation and develop damage criterion for parent metal and weld, respectively; experimental work (uniaxial and multi-axial interrupted creep test) to be carried out or gathered under low stress; cavitation to be quantified, ideally using X-ray micro-tomography. new damage criterion shall be developed; (2) Quantifying the microstructural evolutions and their effects on the creep deformation; (3) Developing and/or applying the novel creep formulation suitable for a wider range of stress, and incorporating the damage criterion developed in (1); (4) Generalising uniaxial version into three dimensional one.

Aims of this paper
This paper aims to overcome the inaccuracies involved in the process of modelling minimum creep strain rate and stress, and in turn develop and apply a more scientifically sound creep fracture criterion.

The development of new minimum creep strain rate and stress function
The process for this part of work is: (1) to identify the deficiencies related to the existing functions, (2) to propose and apply a new mathematical function, and (3) to compare the predictions with experimental data in order to demonstrate and/or verify its capability. The chosen experimental data of minimum creep strain rate and stress for this part of the work are: (1) Creep Test One: low Cr steel (2·25Cr-1Mo steel), stress range of 60-180 MPa, temperature of 565 °C [8], according Xu [9]; (2) Creep Test Two: low Cr Steel (0·5Cr-0·5Mo-0·25 V steel), stress range of 70-180 MPa, temperature of 565 °C [10], according to Xu [9]; (3) Creep Test Three: high Cr steel, P91, stress range of 70-200 MPa; temperature of 600 °C [11].

The method of the development and application of creep fracture criterion based on cavity area fraction along grain boundary for high Cr steel
(1) To derive a specific quantitative equation for cavity area fraction along grain boundary based on the appropriate cavity nucleation and cavity growth models for a high Cr steel (namely P91); (2) To determine the material parameter U′ of cavity area fraction equation according to the experimental data on creep rupture time under different stress levels and temperature; (3) To apply the obtained cavity nucleation model to other high Cr steels in order to reveal the trend of cavity nucleation rate coefficient with stress.
The experimental data chosen for this part of work include: (1) The creep rupture time for P91 [11], stress range of 70-200 MPa, temperature of 600 °C, lifetime between 971 and 80736 h; temperature of 625 °C, lifetime 99-21372 h, for the work described in step (2) of the above; (2) Cavitation information (cavity number density, stress) of high Cr steel (CB8, 10%wt Cr) [12], stress range of 120-180 MPa, temperature of 600 °C lifetime between 2800 and 51406 h, for the work described in step (3) of the above.

Modified hyperbolic sine law for minimum creep strain rate and stress for low Cr steel
The conventional hyperbolic sine law [3,13] is given as: The effectiveness and/or deficiency of its application to 2·25Cr-1Mo and 0·5Cr-0·5Mo-0·25 V low Cr steel is clearly demonstrated by Figures 2 [9] and 3 [9], respectively. Hence, a modified hyperbolic sine law was proposed by the first author and applied in [9], it is given as: where q is stress exponent.

The application of the modified hyperbolic sine law to high Cr steel
Deficiencies of conventional functions that measure the relationship between minimum creep strain rate and stress are demonstrated below.
The conventional minimum creep strain rate and stress functions were proposed and compiled in Table 2 [14].
The specific data of minimum creep strain rate for P91 (9Cr-1Mo-V-Nb) for this research was taken from NIMS creep data sheet as shown in Figure 6. [11].   [10] of 2·25Cr-1Mo steel [9]. Table 2. the typical functions between minimum creep strain rate and stress [14].
Power law creep [15,16] ̇m in = A n linear +power law [17,18] ̇̇m Hyperbolic sine law [3,13] ̇m in = Asinh(B ) Figure 6. experimental data of minimum strain rate and stress at 600 °C under 70-200 Mpa for p91 steel [11].  where the dimensionless factor I( , , ) is a definite integral: The Equation (3.5) N(R, t) can be converted into the density function of cavity equivalent radius, f(R), and is given as [20]: To completely calibrate the cavity area fraction along the grain boundaries, Equation (3.7), the material constants , , , A 1 , A 2 need to be found.

Specific values of material constants α, β, γ
The physical meaning of these values has been given as: (1) With the value of α = 1, β = 2, the cavity size distribution Equation (3.5) fitted well with experiment of the equivalent diameter against number of cavities [20,21], and it is well known from the literature that the value of β = 1·95 ± 0·05 (closely to 2) is characterised by the constrained diffusional mechanism of void growth [19]; (2) γ = 1 is characterised for continuum cavity nucleation with cavity constrained growth [19].

Further confirmation on the simulation of the value of γ = 1
Determination and proofing of the value of γ is given in this section. The characteristics of the cavity size distribution Equation (3.5) were discussed by Riedel [19], and its shape with γ = 1 is similar to that from the experiment for P91. Additionally, the value of γ = 1 is consistent with (3.8) The modified hyperbolic sine law, ̇m in = Asinh(B q ) , is applied to P91 steel. The modelling result is shown in Figure 9 with the calibrated material parameters of A = 6·09867E-7, B = 2·137E-4 and q = 2. As can be observed, a positive agreement was achieved.
The comparison of the above two functions is shown in Figure 10.

Functions for cavity area fraction and the cavity size distribution
A generic theory of cavity area fraction along grain boundary was proposed by Riedel [19], and it can be summarised as: If the non-stationary growth rate of the cavity radius, Ṙ , and the nucleation rate of cavity, J * , are, respectively: where the unknown constants A 1 , A 2 , α, β and γ are material constants.
Then, the cavity size distribution function, N(R, t), represents the number of voids with radii between R and + dR in the time interval t and t + dt, is given as: Finally, the cavity area fraction, w, can be obtained, based on the above three equations [16]: Figure 9. the modelling result of modified hyperbolic sine law compared with experimental data for p91 steel.

The application of the cavity area fraction along grain boundary on high Cr steel
Inputting the above obtained values of α = 1, β = 2 and γ = 1 into Equations (3.7) and (3.8), the cavity area fraction of w is expressed as: The cavitation damage equation should be incorporated in the future development of creep damage constitutive equations for high chromium steel, due to its much sounder scientific base than the conventional phenomenological approach.
constrained diffusion growth and continuums nucleation [19]. This section confirms γ = 1 numerically by comparing predictions with experimental data reported in literature [21]. It is necessary to point out that the actual experimental data is much denser and only a limited number was taken from the graph for comparison.
Firstly, integrating Equation (3.3) yields: Secondly, using the minimum and maximum cavity diameter of 1·2 and 5·7 μm at 10,200 h, respectively, the values of A 1 , A 2 and C can be solved with = 1, β = 2 and given the value . Finally, it was found that γ is equal to 1, the theoretical modelling of f (R), and by using the above values, it agrees with the experimental data very well as displayed in Figure 11. Hence, this confirms that the value of γ is 1. Figure 10. the comparison between different function of minimum creep strain rate and applied stress for p91 steel. Figure 11. probability density function of cavity equivalent r for p91, experimental data from ref [21]. dependent upon stress and temperature. With = 1, the Equation (3.4) can be converted as: where J is the number density of cavities.
This section investigates the dependency between coefficient A 2 and stress.
Fortunately, there is another set of 3D cavity data by X-ray micro-tomography. The number density of cavities at failure for high Cr steel (CB8) under a range of stress levels (120-180 MPa) are available [12]. The relevant data was extracted by reading the published graph and is shown here: Using the above data in Table 5, the individual value of A 2 at different stress levels is calculated and shown in Table 6 and Figure 13.
It is clear that A 2 only reduces slightly with stress below 150 MPa. This is significant as the above trend could be used for extrapolation to lower stress. (3.14)

The trend of creep life time coefficient U′ with stress levels
Rupture is assumed to occur when w reaches its critical vale, w f , and it is assumed to be ∕4 according to [19]. Hence, Based on creep data sheets of creep rupture time under different stress levels and temperature on typical ASME Grade 91(9Cr-1Mo-V-Nb) steel [15] in Table 3, the value of U ′ was calculated as shown in Table 4. There is a clear trend that can be observed in Figure  12, which provides much more stable scientific findings that can be used to convince researchers, than compared to the trend displayed by the value of A in other cavitation damage Equations (1.3) as shown in Figure 1. This potentially provides a scientifically sound, novel, and rather simple life time prediction method.

The variation between cavity nucleation rate coefficient A 2 and stress
The coefficient A 2 in the cavity nucleation rate Equation (3.4) does not change with time, however it might be (3.13)

Discussion and conclusions
Key conclusions for this research include: (1) A modified sine law for minimum creep strain rate and stress was proposed and used. As can be observed it works well with both low Cr and high Cr steels. Hence it is recommended to be utilised in future creep damage constitutive equations for a wide range of stress situations; (2) A specific advanced cavity area fraction along grain boundary model related to the cavity nucleation and growth kinetics models was developed/ calibrated based on scientific reasoning and 3D X-ray tomography information for high Cr steel; (3) The application of cavity area fraction along grain boundary model w f = U � × t 2 f gave a scientifically sound and simple creep life time prediction capability. Furthermore, there is a clear trend between creep life coefficient and stress and temperature, which allows us to be much more confident in its application; (4) The application of the developed/calibrated creep cavity nucleation model gives a scientifically sound and simple cavity nucleation rate dependence of stress. More so, there is a clear trend between creep cavity nucleation coefficient and stress, thus adding to the reliability of this model.
Suggestions for further work include: (1) To optimize the values for A, B and/or q by using optimization software; (2) To couple creep deformation and novel cavitation equations in a novel hyperbolic sine law; (3) To consider w f in different conditions, which in turn can model and validate the novel creep damage constitutive equations. Figure 13. the trend of cavity nucleation rate coefficient A 2 and stress.