Chemical composition (wt %) of P91 and 316 Stainless Steel
1. Introduction
Components in power plant, chemical plant, manufacturing processes, aeroengines, etc. may operate at temperatures which are high enough for creep to occur [1]. Such components may contain cracks or must be assumed to contain cracks as part of design life or remaining life analyses which are required [2]. In order to perform these analyses a number of approaches have been used, based on, for example, a fracture mechanics approach [3], or a continuum damage mechanics approach [4, 5, 6]. This paper is related to the use of the damage mechanics approach. In particular the methods used to obtain the material constants in the multiaxial form of the creep damage and creep strain equations are described. Most of the constants are obtained by fitting to uniaxial creep data; this is a wellestablished method [7]. However, in this paper, the determination of the multiaxial stress state parameter, α [8], is based on results from compact tension (CT) tests; this approach is novel and results in properties which are particularly suited for predicting creep crack growth in components, where the crack growth is defined by a damage parameter, ω. When this damage parameter reaches a critical value (0.99 chosen for the presented work) the material is regarded as ‘completely damaged’ and hence a void or crack growth is assumed to be present. A previously used technique for obtaining the multiaxial stress state parameter, based on the notch strengthening which usually occurs in Bridgman notch [9] creep rupture tests, relative to corresponding uniaxial tests, does not closely represent the stress states and constraint which occur at crack tips. The validity of the method proposed has been established by comparing finite element predictions of creep crack growth in thumbnail cracked specimens with experimental data [7] using the material constants obtained from uniaxial creep and CT creep test results.
The material chosen for the investigation is a 316 stainless steel and a P91 steel because of the ready availability of uniaxial creep, uniaxial creep rupture, compact tension creep crack growth and thumbnail creep crack growth data at temperatures of 600°C and 650°C, respectively. The particular form of damage equation chosen for the investigation is that proposed by Liu and Murakami [6]. By comparison with the more commonly used Kachanov damage equations [4], it was found that the Liu and Murakami equations do not cause the time steps in the finite element analyses to become impractically small [10] and unlike the Kachanov equations, they produce results which are relatively insensitive to element size near the crack tip. These aspects are covered further in the paper.
2. Experimental testing
Two materials have been used for the experimental testing presented, namely P91 steel and 316 stainless steel. The modified 9Cr (P91) steel was initially developed in the US in the early 1980s and was introduced to UK power plants in the early 1990s, to replace some of the components made from low alloy ferritic steels, as its high creep strength allows the use of thinner walled components, which will be less prone to thermal fatigue cracking. The 316 stainless steel is also a creep resistant steel, which is widely used in power plants at high temperature. Table I shows the chemical composition of the P91 steel and 316 stainless steel. All tests for P91 were performed at 650°C [11] and all tests for 316 stainless steel were performed at 600°C [10].















8.74    0.98  0.36  0.022  0.08  0.21      0.11  0.12  0.048  Balance 

16.8  11.8  2.15  1.42  0.5  0.49  0.08  0.07  0.03  0.02  0.02    Balance 
Three main specimen types have been used in order to obtain the experimental data shown in this paper, namely, uniaxial specimens, compact tension (CT) crack growth and thumbnail crack growth creep specimens, as shown in Figure 1, Figure 2 and Figure 3, respectively. Testing was also carried out using sidegrooved CT specimens (see Figure 4) for P91. Tunnelling behaviour was observed for the plain specimens and relatively uniform creep crack growth fronts were observed for the side grooved specimens.
3. Liu and Murakami creep damage model
3.1. Definition of the model
The governing equations for the Liu and Murakami creep damage model are shown by equations (1), (2) and (3).
where A, n, B, q_{2} and χ are material constants.
3.2. Determination of the material constants
The required material constants shown in equations (1) and (2), i.e.
3.2.1. Uniaxial constants
From equation (1) the relationship between the minimum strain rate and stress can be given by [10]:
Therefore, using experimental uniaxial creep data to plot
Therefore, plotting
In order to obtain
An example of this plot using uniaxial creep data for 316 stainless steel, at 600°C, is shown by Figure 7.
3.2.2. Multiaxiality parameter, α
Equation (3) is used for the rupture stress,
It can therefore be seen that the failure time is dependent on this multiaxial constant,
3.3. Material constants
The material constants obtained for P91 and 316 stainless steels are given in Table 2.








1.09×10^{20}  8.462  2.95×10^{16}  6.789  3.2  0.313 

1.47×10^{29}  10.147  2.73×10^{30}  10.949  6.35  0.478 
3.4. Predictive capability of the model
3.4.1. P91 at 650˚C
A typical threedimensional FE mesh and 0.99 damage (crack) zone for the CT specimen geometry is shown in (plain CT specimen), where due to two planes of symmetry in a CT specimen, only one quarter of the specimen has been modelled, with the appropriate boundary conditions applied [11]. Testing and modelling has been carried out for P91 at 650˚C for both plain and sidegrooved CT specimens (see Figure 2 and Figure 4), with the model constants being calculated as shown in section 3.2. Figure 10 shows an example of a tested CT specimen of each type and shows the difference in the corresponding characteristic crack front shapes. Examples of FE creep crack growth modelling of P91 CT specimens using the damage mechanics approach are illustrated in Figure 10 and Figure 11. Figure 10 shows the damage contours, at times close to fracture, of a plain specimen and a sidegrooved specimen. It can be seen from Figure 10 that the tunnelling effect observed in the plain specimens and the essentially uniform crack growth observed in the sidegrooved specimens, as shown in Figure 10, have been reasonably accurately reproduced by the FE analyses. In addition and more importantly, as can be seen in Figure 11, the FE damage modelling has reasonably accurately predicted the creep crack growth behaviour for the P91 CT specimens, when compared with the corresponding experimental results.
3.4.2. 316 stainless steel at 600˚C
The comparisons of the experimental and FE creep crack growths for three plain CT specimens, each subjected to a different test load, are shown in Figure 12, from which it can be seen that the crack front shapes, as well as the extents of creep crack growth, were accurately predicted.
As the multiaxial constant, α, was determined using the CT crack growth data, it is to some extent not surprising that the FE crack growth predictions correspond well to this experimental data, with all of the other material constants having been determined using data from uniaxial creep data. However, similar simulations have been performed for thumbnail crack geometries using the same constants as for the CT specimen and can therefore be considered as ‘pure prediction’. Figure 13 shows an example of the 3dimensional mesh (and 0.99 damage (crack) zone) used for the thumbnail crack growth simulations. As with the CT specimens, due to two planes of symmetry in a thumbnail crack specimen, only one quarter of the specimen has been modelled, with the appropriate boundary conditions applied [10].
The comparisons of the experimental and FE creep crack growths for five thumbnail specimens, each subjected to a different test load, are shown in Figure 14, from which it can be seen that similarly to the CT predictions, the crack front shapes, as well as the extents of creep crack growth were accurately predicted.
4. Discussion and future work
A comprehensive procedure for the determination of the material constants for the Liu and Murakami creep damage model, based on experimental data, has been described. Particular attention has been given to ensuring a constant of multiaxiality value (α) which is highly appropriate to crack tip conditions. These constants have been applied, for a 316 stainless steel at 600˚C and a P91 steel at 650˚C, to a user subroutine for the Liu and Murakami model which has been used in conjunction with Finite Element package ABAQUS, in order to provide numerical predictions for creep crack growth in both compact tension specimen and thumbnail specimen geometries. Comparisons of the model predictions to corresponding experimental data for multiple specimen geometries, in terms of both crack growth and final crack length/profile, show extremely close correlation.
Also shown is the effect that sidegrooves have on the crack profile in a CT specimen and the ability of the Liu and Murakami creep damage model to predict this more uniform crack profile observed in sidegrooved CT specimens.
Nomenclature
Acknowledgments
The authors would like to thank Dennis Cooper, Brian Webster and Shane Maskill for their assistance with the experimental work.
References
 1.
R. K. Penny and D. L. Marriott, “Design for Creep”, McGrawHill , Liverpool, 1971.  2.
T. H. Hyde, W. Sun and A. A. Becker, "Creep crack growth in welds: A damage mechanics approach to predicting initiation and growth of circumferential cracks in CrMoV weldments", Int. J. Pres. Ves. & Piping , 2001; 78, 765771.  3.
B. Dogan and B. Ptrovski, "Creep crack growth of high temperature weldment." International Journal of Pressure Vessel and Piping , 2001, 78, 795805.  4.
L. M. Kachanov, "The time to failure under creep condition", Izv. Akad. Nauk., SSSR. Tekh. Nauk , 1958, 8, 2631.  5.
Y. N. Robotnov, "Creep Problems of Structural Members”, NorthHolland , 1969.  6.
Y. Lui and S. Murakami, "Damage localization of conventional creep damage models and proposition of a new model for creep damage analysis", JSME International Journal , 1998, 41, 5765.  7.
T. H. Hyde, "Creep crack growth in 316 stainless steel at 600°C", High Temperature Technology , 1988, 6(2), 5161.  8.
D. R. Hayhurst and C. J. Morrison, "Development of continuum damage in the creep rupture of notched bars", Phil. Trans. R. Soc. Lond. (A) , 1984, 311, 103129.  9.
G. A. Webster, S. R. Holdsworth, M. S. Loveday, I. J. Perrin, K. Nikbin, H. Purper, R. P. Skelton, and M. W. Spindler, "A code of practice for conducting notched bar creep rupture tests and for interpreting the data", J. Fatigue and Fatigue of Eng. Materials and Struct , 2004, 24, 319342.  10.
C. J. Hyde, T. H. Hyde, W. Sun and A. A. Becker, “Damage mechanics based predictions of creep crack growth in 316 stainless steel”, Engineering Fracture Mechanics , 2010, 77(12), 23852402.  11.
T. H. Hyde, M. Saber and W. Sun, “Testing and modelling of creep crack growth in compact tension specimens from a P91 weld at 650°C,” Engng. Frac. Mech , 2010, 77(15), 29462957.  12.
T. H. Hyde and W. Sun, “Determining high temperature properties of weld materials,” JSME Int. J. of Solids Mechanics & Material Eng, Series A , 2000, 43(4), 408414.