Three‐dimensional analyses of in‐plane and out‐of‐plane crack‐tip constraint characterization for fracture specimens

Three‐dimensional elastic–plastic finite element analyses have been conducted for 21 experimental specimens with different in‐plane and out‐of‐plane constraints in the literature. The distributions of five constraint parameters (namely T‐stress, Q, h, T_z and A_p) along crack fronts (specimen thickness) for the specimens were calculated. The capability and applicability of the parameters for characterizing in‐plane and out‐of‐plane crack‐tip constraints and establishing unified correlation with fracture toughness of a steel were investigated. The results show that the four constraint parameters (T‐stress, Q, h and T_z) based on crack‐tip stress fields are only sensitive to in‐plane or out‐of‐plane constraints. Therefore, the monotonic unified correlation curves with fracture toughness (toughness loci) cannot obtained by using them. The parameter A_p based on crack‐tip equivalent plastic strain is sensitive to both in‐plane and out‐of‐plane constraints, and may effectively characterize both of them. The monotonic unified correlation curves with fracture toughness can be obtained by using A_p. In structural integrity assessments, the correlation curves may be used in the failure assessment diagram (FAD) methodology for incorporating both in‐plane and out‐of‐plane constraint effects in structures for improving accuracy.     

In‐plane and out‐of‐plane constraint parameters along a three‐dimensional crack‐front stress field under creep loading

Full‐field three‐dimensional (3D) numerical analyses was performed to determine in‐plane and out‐of‐plane constraint effect on crack‐front stress fields under creep conditions of finite thickness boundary layer models and different specimen geometries. Several parameters are used to characterize constraint effects including the non‐singular T‐stresses, the local triaxiality parameter, the T_z ‐factor of the stress‐state in a 3D cracked body and the second‐order‐term amplitude factor. The constraint parameters are determined for centre‐cracked plate, three‐point bend specimen and compact tension specimen. Discrepancies in constraint parameter distribution on the line of crack extension and along crack front depending on the thickness of the specimens have been observed under different loading conditions of creeping power law hardening material for various configurations of specimens.     

Unified correlation of in‐plane and out‐of‐plane constraints with fracture toughness

In this paper, the specimens with different geometries and loading configurations were used to study the unified correlation of in‐plane and out‐of‐plane constraints with fracture toughness by using numerical simulation method. The results show that the unified constraint parameter A_p which was defined on the basis of the areas surrounded by the equivalent plastic strain isolines ahead of crack tip can characterise both in‐plane and out‐of‐plane constraints induced by different specimen geometries and loading configurations. A sole linear relation between the normalised fracture toughness J_IC/J_ref and was obtained. The J_IC/J_ref ‐ line is a unified correlation line of in‐plane and out‐of‐plane constraints with fracture toughness of a material, and the constraint dependent fracture toughness of a material can be determined from the unified correlation line. The results also demonstrate that the out‐of‐plane constraint effect is related to the in‐plane constraint effect, and there exists interaction between them. The higher in‐plane constraint strengthens the out‐of‐plane constraint effect, whereas the lower in‐plane constraint is not sensitive to the out‐of‐plane constraint effect.     

Unified characterisation of in‐plane and out‐of‐plane constraint based on crack‐tip equivalent plastic strain

Constraint can be divided into two conditions of in‐plane and out‐of‐plane, and each of them has its own parameter to characterize. However, in most cases, there exists a compound change of both in‐plane and out‐of‐plane constraint in structures, a unified measure that can reflect both of them is needed. In this paper, the finite element method (FEM) was used to calculate the equivalent plastic strain (ɛ_p) distribution ahead of crack tips for specimens with different in‐plane and out‐of‐plane constraints, and the FEM simulations based on Gurson–Tvergaard–Needleman (GTN) damage model and a small number of tests were used to obtain fracture toughness for the specimens with different constraints. Unified measure and characterisation parameter of in‐plane and out‐of‐plane constraints based on crack‐tip equivalent plastic strain has been investigated. The results show that the area A_PEEQ surrounded by the ɛ_p isoline ahead of crack tips can characterize both in‐plane and out‐of‐plane constraints. Based on the area A_PEEQ, a unified constraint characterisation parameter A_p was defined. It was found that there exists a sole linear relation between the normalised fracture toughness J_IC/J_ref and regardless of the in‐plane constraint, out‐of‐plane constraint and the selection of the ɛ_p isolines. The unified J_IC/J_ref−reference line can be used to determine constraint‐dependent fracture toughness of materials. The FEM simulations with the GTN damage model (local approach) can be used in obtaining the unified J_IC/J_ref−reference line for materials with ductile fracture.     

Effects of side‐groove depth on creep crack‐tip constraint and creep crack growth rate in C(T) specimens

The effects of side‐groove depth on creep crack‐tip constraint and creep crack growth (CCG) rate in C(T) specimens have been quantitatively studied. The results indicate that with increasing side‐groove depth, the constraint level and CCG rate increase and constraint distribution along crack front (specimen thickness) becomes more uniform. The constraint and CCG rate of thinner specimen are more sensitive to side‐groove depth. Two new creep constraint parameters (namely R^* and A_c) both can quantify constraint levels of the specimens with and without side‐grooves, and the quantitative correlations of CCG rate with constraint have been established. The mechanism of the side‐groove depth effect on the CCG rate has also been analyzed.     

In‐plane and out‐of‐plane constraint for single edge notched bending specimen and cruciform specimen under uniaxial and biaxial loading

Considering fracture constraint is an efficient way to describe stress–strain field and fracture toughness more accurately, so it is necessary to realise the relationship with in‐plane and out‐of‐plane constraint for different standard specimens. In this paper, three‐dimensional finite element method is applied to study the in‐plane and out‐of‐plane constraint for both cruciform specimen and single edge notched bending specimen made from commercial pure titanium. Crack length and in‐plane loading as the factors affecting in‐plane constraint, and thickness as the factor affecting the out‐of‐plane constraint are used to study the effect on both in‐plane and out‐of‐plane constraint in this paper. From the results, in‐plane and out‐of‐plane constraint are both related to specimen geometries and loading styles. And there exist relationships with in‐plane and out‐of‐plane constraint because of factors for different specimens. Depending on crack length, out‐of‐plane constraint increases with in‐plane constraint. While depending on transverse loading, out‐of‐plane constraint decreases with in‐plane constraint. In addition, when the in‐plane constraint of a specimen is higher, in‐plane constraint increases with out‐of‐plane constraint (thickness). When the in‐plane constraint is lower, in‐plane constraint almost remains unchanged with out‐of‐plane constraint.     

Correlation of material constraint with fracture toughness of interface regions in a dissimilar metal welded joint

In this work, the constraint parameter A_p based on crack‐tip equivalent plastic strain was calculated by finite element analyses for the cracks located at different locations in two interface regions in a dissimilar metal weld joint (DMWJ). The capabilities of the parameter A_p for characterizing material constraint and establishing correlation of material constraint with fracture toughness of the interface region cracks have been examined. The results show that the parameter A_p can characterize material constraint effect caused by material mismatch and initial crack positions in the interface regions. Based on the A_p, the correlation lines and formulae of material constraint with fracture toughness of the interface region cracks in the DMWJ can be established, and they may be used for obtaining material constraint‐dependent fracture toughness for the interface region cracks. The results in this work combining with those in the previous studies indicate that the parameter A_p may be a unified constraint parameter that can characterize both geometry constraint (including in‐plane and out‐of‐plane constraints) and material constraint, and it may be used in accurate fracture assessments of welded components with different geometry and material constraints.     

基于C(t)积分及参数Ac的P92钢高温蠕变裂纹扩展研究

对于在高温环境下工作的构件,蠕变裂纹扩展是一种主要的失效机制,而裂纹尖端的拘束水平对蠕变裂纹扩展率有很大的影响。通过数值仿真与相关试验数据对比的方法,对裂纹扩展尖端的应力应变率场表征参量C(t)积分进行了相关研究,并基于参数Ac研究了P92材料裂纹尖端的拘束水平对蠕变裂纹扩展的影响。研究结果表明,C(t)积分值随裂纹扩展急剧减小,其数值及变化与积分路径到裂纹尖端的距离相关性很强,并且与拘束水平有一定的关系;拘束水平影响蠕变裂纹扩展率,拘束越大,裂纹扩展速率越快;参数Ac可以有效表征裂纹尖端拘束水平,其在寿命预测方面的应用有待进一步研究,同时在含裂纹的高温工作构件寿命评估方面有重大的意义。     

Effects of creep properties of materials on unified creep constraint parameter Ac for cracked pipes

Extensive finite element analyses of cracked pipes with different crack sizes and orientations have been conducted to investigate effects of creep properties of materials on the unified creep constraint parameter A_c. The results show that the constraint parameter A_c is independent on Norton’s coefficient A, and it is only affected by the creep exponent n of materials. For a given crack size, with increasing n, A_c decreases and constraint level increases. The A_c of lower constraint cracks is more sensitive to n. The unified correlation equations between A_c and n have been obtained for cracked pipes with a wide range of crack sizes and constraint levels. They may be used to estimate the constraint parameter A_c at different positions along the crack fronts in cracked pipes made of materials with different n values. The two-parameter C*-A_c approach for assessing creep life of cracked pipes has also been discussed.     

Estimation of high‐temperature fracture parameters for small punch specimen with a surface crack

An outline of a newly proposed methodology for evaluating creep crack growth (CCG) parameters using cracked small‐punch (SP) specimens is explained. Three‐dimensional finite element analyses were performed to calculate the stress intensity factor along the crack front for a surface crack formed at the centre of a SP specimen. Effects of crack ratio, (a/t); crack aspect ratio, (a/c); and thickness of the specimen, (t), on the fracture parameters were studied. It was observed that the minimum variation of K‐value along the crack front can be achieved when a/c was 0.50 except the location very near the intersection of the crack and free surface. This condition is similar to the case of constant K‐values along the crack front of the conventional compact tension specimen. Thus, it can be argued that the SP specimen with a surface crack is a suitable specimen geometry for CCG testing. The proposed CCG test method was found to be practically applicable for the crack geometry of 0.10 to 0.30 of a/t with constant aspect ratio of 0.50. An estimation of the K and C_t‐parameter under the small scale creep condition was derived. Future work for further development of the suggested CCG testing is discussed.     

Numerical investigation of creep crack growth in cross‐weld CT specimens. Part II: influence of specimen size

A numerical investigation of the influence of specimen size on creep crack growth in cross‐weld CT specimens with material properties of 2.25Cr1Mo at 550 °C is performed. A three‐dimensional large strain and large displacement finite element study is carried out, where the material properties and specimen size are varied under constant load for a total of eight different configurations. The load level is chosen such that the stress intensity factor becomes 20 MPa √m regardless of specimen size. The creep crack growth rate is calculated using a creep ductility‐based damage model, in which the creep strain rate ahead of the crack tip perpendicular to the crack plane is integrated taking the degree of constraint into account. Although the constraint ahead of the crack tip is higher for the larger specimens, the results show that the creep crack growth (CCG) rate is higher for the smaller specimens than for the larger ones. This is due to much higher creep strain rates ahead of the crack tip for the smaller specimens. If, on the other hand, the CCG rate is evaluated under a constant C * condition, the creep crack growth rate is found to be higher for the larger specimens, except when the crack is located in a HAZ embedded in a material with a lower minimum creep strain rate; then, the creep crack growth rate is predicted to be higher for the smaller specimen. In view of these results, it is obvious that the size effect needs to be considered in assessments of defected welded components using results from CCG testing of cross‐weld CT specimens.     

Prediction of creep crack initiation in Cr–Mo–V steel specimens with different geometries

By using stress dependent creep ductility and strain rate model in a ductility exhaustion based damage model, the creep crack initiation (CCI) behaviour in Cr–Mo–V steel specimens with different geometries and dimensions (different constraints) over a wide range of C^* has been predicted by finite element simulations. The predicted creep crack initiation time agree well with the existing experimental data. In low and transition C^* regions, the constraint induced by specimen geometries and dimensions has obvious influence on CCI time. With increasing constraint level of specimens, the CCI time decreases due to the increase of stress triaxiality ahead of crack tip. Different CCI trends and constraint effects on CCI behaviour in a wide range of C^* result from the interaction of crack-tip stress state and stress dependent creep ductility of the steel. It is suggested that in CCI life assessments of high temperature components, the long-term CCI time data at low C^* region should be obtained and used, and the constraint effects need to be considered by using constraint dependent CCI data.     

Creep crack initiation prediction considering constraint effect for pressurized pipelines with circumferential surface cracks

Theoretical 2‐parameter C*‐Q* approaches and numerical simulation were conducted to investigate the creep crack initiation (CCI) time and the effect of constraints induced by the geometrical sizes of pipelines with circumferential surface cracks. The theoretical enhancement model of the C*‐Q* approaches under the transient creep condition, which considered the load‐independent constraint parameter Q*, was proposed to predict the CCI time around the crack front. The results revealed that the distribution regulation of Q* along the crack front for circumferential internal surface cracks and external surface cracks was similar. The maximum constraint level occurred near the deepest crack front part for cracks with small a/c (a/c < 0.4), while it occurred near the free surface for cracks with large a/c (a/c > 0.4). The constraint values at the same position (2Φ/π) increased with the increasing of the crack depth when a/c kept constant. In addition, the circumferential internal surface cracks of pipelines were proved more dangerous than the external surface cracks with the same geometrical size. Furthermore, the CCI times were decided by the peak values of constraint, or the CCI firstly occurred at the position where the constraint level was maximum. Additionally, the variation of hydrostatic stresses, triaxiality, and multiaxial strain factor considering the constraint parameter Q* was discussed. The suitability of the analytical C*‐Q* approaches was verified to predict CCI. The comparison of CCI times between the analytical approach and the BS 7910 as well as the finite element results demonstrated that the solutions under stress intensity factor–Riedel‐Rice control (initially by stress intensity factor and then by transient creep stress or Riedel‐Rice conditions)—were more accurate when internal pressure P < 18 MPa, but the solutions under Hutchinson‐Rice‐Rosengren–Riedel‐Rice control (initially by plastic Hutchinson‐Rice‐Rosengren control and then by Riedel‐Rice conditions) were more appropriate when P > 18 MPa. Finally, the accuracy of C*‐Q* 2‐parameter approaches and the suitability of finite element simulation were validated.     

Experimental and numerical analysis of in-plane and out-of-plane crack tip constraint characterization by secondary fracture parameters

In this study, several two-parameter- concepts are analyzed experimentally and numerically with respect to their capability of characterizing in-plane and out-of-plane crack tip constraint effects. Different approaches utilizing the second term T _stress of the linear-elastic crack tip stress field, a higher term A ₂ of the power-law hardening crack tip stress field, a hydrostatic correction term Q for a reference stress field or the local triaxiality parameter h are compared. Experimental results for a pressure vessel steel 22NiMoCr3-7 are investigated by means of the different approaches regarding their capability of constraint characterization for enhanced transferability. Theoretical aspects are investigated in a modified boundary layer analysis and in three-dimensional nonlinear elastic-plastic finite element analyses of the specimens. It is found that, with respect to their capability of quantifying combined in-plane and out-of-plane constraint effects, the investigated concepts differ significantly.     

Effect of residual stresses on crack behaviour in single edge bending specimens

Residual stresses due to manufacturing processes, such as welding, change the load bearing capacity of cracked components. The effects of residual stresses on crack behaviour in single edge bending specimens were investigated using Finite element analyses. Three parameters (J, Q and R) were used to study the crack behaviour. The J‐integral predicts the size scale over which large stresses and strains exist, the constraint parameter Q describes the crack‐tip constraint as a result of geometry, loading mode and crack depth and the constraint parameter R is used to describe the constraint resulting from residual stresses. To carry out a systematic investigation on the effect of residual stresses on the J‐integral and crack‐tip constraints, models under different combinations of residual stresses and external loads with different crack depths were analysed. It has been shown that the crack‐tip constraint R increased by tensile residual stresses around the crack‐tip. On the other hand, the constraint parameter R decreased and tended to zero at high external load levels.     

DAMAGE MECHANICS APPROACH TO REMOVE THE CONSTRAINT DEPENDENCE OF ELASTIC–PLASTIC FRACTURE TOUGHNESS

It is now generally agreed that the applicability of a one-parameter J-based ductile fracture approach is limited to so-called high constraint crack geometries, and that the elastic-plastic fracture toughness J_1c, is not a material constant but strongly specimen geometry constraint-dependent. In this paper, the constraint effect on elastic-plastic fracture toughness is investigated by use of a continuum damage mechanics approach. Based on a new local damage theory for ductile fracture(proposed by the author) which has a clear physical meaning and can describe both deformation and constraint effects on ductile fracture, a relationship is described between the conventional elastic-plastic fracture toughness, J_1c, and crack tip constraint, characterized by crack tip stress triaxiality T. Then, a new parameter J_dc (and associated criterion, J_d=J_dc) for ductile fracture is proposed. Experiments show that toughness variation with specimen geometry constraint changes can effectively be removed by use of the constraint correction procedure proposed in this paper, and that the new parameter J_dc is a material constant independent of specimen geometry (constraint). This parameter can serve as a new parameter to differentiate the elastic-plastic fracture toughness of engineering materials, which provides a new approach for fracture assessments of structures. It is not necessary to determine which laboratory specimen matches the structural constraint; rather, any specimen geometry can be tested to measure the size-independent fracture toughness J_dc. The potential advantage is clear and the results are very encouraging.     

Modelling the effect of creep–fatigue interaction on crack growth

Time-dependent creep–fatigue crack growth (CFCG) is a major consideration in estimating the remaining life of elevated temperature components. Fracture mechanics approaches have proven useful in providing a framework for characterizing crack growth under service conditions, and in defining safe operating conditions and selecting inspection criteria and intervals. Experimental and analytical approaches have been developed to characterize crack growth under combined creep and fatigue loading conditions using (C_t )_avg as the crack tip parameter. The analytical approaches that have been proposed to characterize CFCG are limited in their application because they do not completely account for the effect of creep–fatigue interactions in modelling crack tip deformation, and thus, accurately estimating the (C_t )_avg value. A new creep-reversal parameter, C_R , is defined in this study to quantify the extent of creep–fatigue interaction at the crack tip, and is used in an analytical scheme, suitable for components, for calculating (C_t )_avg . This approach does not rely on any simplifying assumptions regarding the extent of reinstatement of C_t , which is dependent on the amount of creep reversal due to cyclic plasticity, during the unloading part of a trapezoidal loading waveform cycle. The (C_t )_avg values calculated by this approach compare well with the experimentally obtained values for compact type (CT) specimens, thus providing an experimental verification of the approach.     

Numerical prediction of creep crack growth in different geometries using simplified multiaxial void growth model

Abstract

This paper considers the prediction of creep crack growth (CCG) in different fracture mechanics geometries using finite element (FE) analysis based on a material independent simplified multiaxial failure strain model at the crack tip. The comparison is first made by modelling C(T) specimen tests under plane stress and plane strain conditions using creep properties of a C–Mn steel at 360°C. In addition, in order to examine CCG due to different geometries, a single edge notch specimen (SENT), centre cracked tension specimen (CCT) and three-point bending (3PB) specimen have been modelled and analysed. In all cases, it is found, depending on the geometry, that for this steel at low creep temperatures the applied load develops a high reference stress/yield stress (σ_ref/σ_y) ratio, which helps reduce constraint at the crack tip. The predictions are analysed under plane stress/plane strain loading conditions identifying the effects of geometry on cracking rates and the implications for predicting long term test or component failure times exceeding where the applied σ_ref/σ_y<<1.     

General time‐dependent C(t) and J(t) estimation equations for elastic‐plastic‐creep fracture mechanics analysis

This paper presents equations for estimating the crack tip characterizing parameters C(t) and J(t), for general elastic‐plastic‐creep conditions where the power‐law creep and plasticity stress exponents differ, by modifying the plasticity correction term in published equations. The plasticity correction term in the newly proposed equations is given in terms of the initial elastic‐plastic and steady‐state creep stress fields. The predicted C(t) and J(t) results are validated by comparison with systematic elastic‐plastic‐creep FE results. Good agreement with the FE results is found.     

Transferability of specimen J–R curve to straight pipe with circumferential surface flaw

Specimen J–R curve is extensively used for structural integrity of large components. It is well known that J–R curve heavily depends on constraint level ahead of crack tip in remaining ligament. In earlier work, it was demonstrated that J–R curve from Three Point Bending (TPB) specimen is transferable to straight pipe with circumferential through wall crack. In this paper, the transferability of J–R curve is investigated from TPB specimen to pipe with circumferential surface crack. A 16 in. diameter pipe with circumferential surface crack and TPB specimen machined from same piping material (SA333Gr6 Steel) are tested. Consequently, 3D finite element analysis (FEA) has been performed on surface cracked pipe and TPB specimen. Crack‐initiation load is also predicted for surface cracked pipe by FEA and compared with experimental result. J–R curve is calculated for the pipe using experimental data, that is, load, load line displacement and crack growth. J–R curve of pipe is compared with TPB specimen and it is found that the pipe is predicting much higher J–R curve than TPB. This difference of J–R curve is investigated by evaluating stress triaxiality in remaining ligament for both cases. Stress triaxiality is quantified using triaxiality factor (h) ahead of crack tip for pipe and TPB specimen. It is found that the TPB specimen has considerably higher constraint level than pipe with surface crack, which is well supported by trend of J–R curves for specimen and pipe. A study has also been carried out to investigate the effect of internal pressure on the stress triaxiality. It is found that there is negligible difference in stress triaxiality because of internal pressure. The stress triaxiality is re‐established as a qualitative parameter to assess the transferability of J–R curve from specimen to component.