Three-dimensional interpretation of cleavage fracture tests of cladded specimens with local approach to cleavage fracture

Electricité de France has conducted during these last years an experimental and numerical research programme in order to evaluate fracture mechanics analyses used in nuclear reactor pressure vessels integrity assessment, regarding the risk of brittle fracture. Two cladded specimens made of ferritic steel A508 Cl3 with stainless steel cladding, and containing shallow subclad flaws, have been tested in four point bending at very low temperature to obtain cleavage failure. The crack instability was obtained in base metal by cleavage fracture, without crack arrest. The tests have been interpreted by local approach to cleavage fracture (Beremin model) using three-dimensional finite element computations. After the elastic–plastic computation of stress intensity factor K_J along the crack front, the probability of cleavage failure of each specimen is evaluated using m, σ_u Beremin model parameters identified on the same material. The failure of two specimens is conservatively predicted by both analyses. The elastic–plastic stress intensity factor K_J in base metal is always greater than base metal fracture toughness K_1c. The calculated probabilities of cleavage failure are in agreement with experimental results. The sensitivity of Beremin model to numerical aspects is finally exposed.     

Effect of reverse cyclic loading on the fracture resistance curve in C(T) specimen

Fracture resistance (J–R) curves, which are used for elastic–plastic fracture mechanics analyses, are known to be dependent on the cyclic loading history. The objective of this paper is to investigate the effect of reverse cyclic loading on the J–R curves in C(T) specimens. The effect of two parameters was observed on the J–R curves during the reverse cyclic loading. One was the minimum-to-maximum load ratio (R) and the other was the incremental plastic displacement (δ_cycle/δ_i), which is related to the amount of crack growth that occurs in a cycle. Fracture resistance tests on C(T) specimens with varying the load ratio and the incremental plastic displacement were performed, and the test results showed that the J–R curves were decreased with decreasing the load ratio and decreasing the incremental plastic displacement. Direct current potential drop (DCPD) method was used for the detection of crack initiation and crack growth in typical laboratory J–R tests. The values of crack initiation J-integral (J_I) and crack initiation displacement (δ_i) were also obtained by using the DCPD method.     

Application of a local approach to ductile-brittle transition in a low-alloyed steel

Both tensile tests on notched specimens and fracture mechanics experiments on axisymmetrically cracked specimens were performed on one heat of A508 steel (AFNOR: 16MND5). Tensile tests on notched geometries were made to determine the characteristic parameters used in a statistical analysis of cleavage fracture proposed previously [1]. Tests on cracked specimens were carried out between −80°C and −20°C to investigate the critical values of stable crack growth, Δa_c, occurring before unstable cleavage fracture. At a given temperature a large scatter in the values of critical crack growth, Δa_c, was observed.To interpret these results a model derived previously for cleavage crack initiation [1] is used. In this model the Weibull stress is calculated by the finite element method for three different initial crack lengths covering stable growth increments observed experimentally. It is shown that this model accounts reasonably well for the observed effects.     

Large stable crack growth in fracture mechanics specimens

According to the J concept, information is reported about the crack resistance behaviour up to 8 mm crack growth of side-grooved CT-25 as well as CCT-25 specimens made from German standard steel StE 460. Numerical simulations controlled by J_R curves make the calculation of J from the stresses and strains of specimen models during large crack growth feasible. These data allow a comparison to standards and rules describing the evaluation of J from experiments. Using stress, strain and displacement fields from a plane-strain finite-element analysis, the extended J concept is discussed concerning larger ductile crack growth. Additionally, the distribution of other fracture mechanics parameters such as the crack tip opening displacement (CTOD) and the crack tip opening angle (CTOA) are presented for larger crack growth.     

The application of constraint based fracture mechanics to engineering structures

Elastic-plastic crack tip fields can be characterised by two parameters, J and T/Q which describes crack tip constraint. This forms the foundation of a constraint based fracture mechanics in which toughness is expressed as a function of constraint in the form of a J–(T/Q) locus. The enhanced toughness associated with shallow cracks and defects can be used in a systematic manner by constraint matching, and implemented through a simple modification to a failure assessment diagram. This methodology enables the conservatism associated with the use of deeply cracked fracture toughness measurements to be relaxed. In this work these methods have been applied to a tubular welded joint.     

2D crack growth simulation with an energetic approach

It is well known that the tearing resistance curve J–Δa is not a material property. A recent approach, based on an energetic critical parameter to model ductile tearing propagation, is used to model 3D effects. The approach considered in this work aims to estimate the dissipated energy in the fracture process during ductile tearing represented by an intrinsic parameter G_fr. A fracture criterion, which accounts for the crack extension length, is defined and lies on a critical energy release rate, noted G_c, which is compared to G_fr. Previously, this parameter was obtained from a numerical local energy release rate, which handicaps the application field of the approach: a fine mesh for the whole propagation area was needed and the criterion allowed only to model 1D propagation. A new manner to estimate G_c is then proposed in this article, based on the J plastic part variation, which allows to model 2D propagation by defining a local criterion. This new calculation method is validated on a CT specimen made in Tu52b ferritic steel, by comparing the results obtained from the two methods of G_c calculation. Then, the 2D crack growth case is studied, by modelling the propagation in a ring, loaded in compression. It is shown that a 3D effect, such as tunnel effect, could be successfully represented with this approach.     

Ductile–brittle transition of polycrystalline iron and iron–chromium alloys

Fracture toughness of polycrystalline Fe, Fe–3%Cr and Fe–9%Cr was measured by four-point bending of pre-cracked specimens at temperatures between 77 K and 150 K and strain rates between 4.46 × 10⁻⁴ and 2.23 × 10⁻² s⁻¹. For all materials, fracture behaviour changed with increasing temperature from brittle to ductile at a distinct brittle–ductile transition temperature (T_c), which increased with increasing strain rate. At low strain rates, an Arrhenius relation was found between T_c and strain rate in each material. At high strain rates, T_c was at slightly higher values than those expected from extrapolation of the Arrhenius relation from lower strain rates. This shift of T_c was associated with twinning near the crack tip. For each material, use of an Arrhenius relation for tests at strain rates at which specimens showed twinning gave the same activation energy as for the low strain rate tests. The values of activation energy for the brittle–ductile transition of polycrystalline Fe, Fe–3%Cr and Fe–9%Cr were found to be 0.21, 0.15 and 0.10 eV, respectively, indicating that the activation energy for dislocation glide decreases with increasing chromium concentration in iron.     

Evaluation of SCC crack behavior in zirconium and zircaloy-2 using nonlinear fracture mechanics parameters

Applicability of nonlinear fracture mechanics parameters, i.e. J-integral, crack tip opening displacement (CTOD), and crack tip opening angle (CTOA), to evaluation of stress corrosion crack (SCC) propagation rate was investigated using fully annealed zirconium plates and Zircaloy-2 tubing, both of which produce SCC with comparatively large plastic strain in an iodine environment at high temperatures.Tensile SCC tests were carried out at 300°C for center-notched zirconium plates and internal gas pressurization SCC tests at 350°C, for Zircaloy-2 tubing, to measure the SCC crack propagation rate. The J-integral around semi-elliptical SCC cracks produced in Zircaloy-2 tubing was calculated by a three-dimensional finite element method (FEM) code.The test results revealed that the SCC crack propagation rate dc/dt could be expressed as a function of the J-integral, which is the most frequently used parameter in nonlinear fracture mechanics, by the equation dc/dt = C · Jⁿ, where C and n were experimental constants.Among the other parameters, CTOD and CTOA, the latter appeared to be useful for assessing the crack propagation rate, because it had a tendency to hold a constant value at various crack depths.     

Study on the effect of the crack length on the JIC value

The ductile fracture tests are carried out using compact type (CT) specimen, three point bend specimen and center cracked tension (CCT) specimens made of A533B steel and aluminum alloy with different crack lengths. The effect of the crack tip constraint on the microscopic fracture behavior is studied by the scanning electron microscope (SEM). It is shown that the apparent J_IC value increases due to the decrease of the crack tip constraint. It is pointed out that the increase of the apparent J_IC value is partly due to the error of the conventional equation to estimate the J value. Based on the FEM analyses, these apparent J_IC values are corrected and are compared with the valid J_IC values. The good co-relation between apparent J_IC value and the Q factor, proposed by O'Dowd and Shih, is shown for all the specimens. The FEM analyses based on the Gurson's constitutive equation is also carried out. The effect of the constraint on the crack tip field is discussed.     

Evaluation of hydrogen assisted cracking: the meaning and significance of the fracture mechanics approach

Analysis of hydrogen assisted cracking is performed with emphasis on the validity of the linear elastic fracture mechanics approach. Its soundness is ensured by the uniqueness of the crack growth kinetics curve ‘crack growth rate v versus stress intensity factor K' for a material–environment system. The ability of this relationship to match the similitude of crack tip events is revised for the whole chain of process components: ‘hydrogen supply to metal—entry into and accumulation in the process zone—degradation and crack advance'. Some of the factors able to violate the uniqueness of the ‘material's' v(K)-curve are discussed: (i) the role of far field in hydrogen transportation to the fracture zone; (ii) near tip environmental interactions and process zone hydrogenation dependent on the cracking process history; (iii) mechanical nonlinear effects caused by crack tip plasticity which may be significant even under small-scale yielding. Suggestions for a refinement of the customary fracture mechanics approach are given.     

Short circumferential through-wall cracked pipes subjected to stretch-bending load

The methods for assessment of elastic–plastic fracture behaviour of cracked components include the net section plastic collapse concept, the J-integral approach, and the two-parameter R-6 failure assessment diagram, Revision 3. These failure assessment methods are usually used to obtain fracture behaviour prediction with a reasonable degree of accuracy without carrying out complicated full-length numerical fracture analysis. In the current work, fracture experiments on stainless steel pipes with short circumferential through-wall cracks under stretch-bending load were conducted. Stretch-bending load refers to the loading situation where axial load is generated that is proportional or related to the applied bending load. The J-integral values derived from the experimental load-point load–displacement data under stretch-bending and pure bending conditions are compared to investigate the effect of axial load on the J–resistance curves. The results show clear dependence of crack resistance force on axial load for short circumferential cracks. Crack resistance force decreased noticeably for increased stretch-bending loading compared to pure bending loading.     

The growth of fatigue cracks in the creep range from the point of view of ductile fracture

This paper deals with experimental and numerical investigations of crack growth in 1% CrMoV steels in the creep temperature range. For the load and displacement controlled fatigue tests with predominantly plastic deformation amplitudes, centre cracked panels and CT-specimens of different thicknesses were used.For describing the crack growth per load cycle the crack tip displacement δ_t and the J-integral were applied, whereby the question arises whether the various fracture concepts are applicable unrestrictedly or where they meet their limits of validity.Detailed theoretical-numerical calculations using the finite element method were expected to yield information about the path independence of the J-integral. Special importance was therefore attached to investigate the J-integral, both cyclical and with reference to the instantaneous state of deformation, by comparing the different J-integral values (ΔJ, J) with the experimental method of Dowling/Begley.     

The use of acoustic emission to determine characteristic dynamic strength and toughness properties of steel

The dynamic J integral at crack initiation (J_id) and dynamic yield stress (σ_yd) are useful parameters to characterize elastic-plastic material behaviour under rapid loading rates. The critical step for evaluating J_id and σ_yd under the condition of the three point impact bending test is the detection of the crack initiation and of the yield point in the impact load–deflection curve, respectively. This paper presents an acoustic emission (AE) based method to determine the ductile crack initiation and additionally the beginning of yield. The experimental techniques used to evaluate σ_yd and J_id include both instrumented pendulum impact tests with the AE transducer within the striker (tup) and medium rate three point bend (TPB) tests with additional AE transducers on the specimen surface. Results obtained from the tests indicate that the AE method is capable of detecting general yielding and the onset of ductile crack growth (initiation). Different types of pulse shaped AE signals can be observed. They were connected with characteristic features during the loading process.     

Observations on the effect of post-weld heat treatment on J-resistance curves of sa-106b seamless piping welds

Ontario Hydro has developed a leak-before-break (LBB) methodology for application to large diameter piping (21, 22 and 24 inch) Schedule 100 SA106B heat transport (HT) piping as a design alternative to pipe whip restraints and in recognition of the questionable benefits of providing such devices. Ontario Hydro's LBB approach uses elastic-plastic fracture mechanics (EPFM).In order to assess the stability of HT piping in the presence of hypothetical flaws, the value of the material J-integral associated with crack extension (J— R curve) must be known. In a material test program J-resistance curves were determined from various pipe heats and four different welding procedures that were developed by Ontario Hydro for nuclear Class 1 piping. The test program was designed to investigate and quantify the effect of various factors such as test temperature, crack plane orientation and welding effects which have an influence on fracture properties. An acceptable lower bound J-resistance curve for the piping steels and welds were obtained by machining maximum thickness specimens from the pipes and weldments and by testing side-grooved compact tension specimens. This paper addresses the effect of test temperature and post-weld heat treatment on the J-resistance curves from the welds.The fracture toughness of all the welds at 250°C was lower than that at 20°C. Welds that were post-weld heat treated showed high crack initiation toughness, J_lc, rising J-resistance curves and stable and ductible crack extension. Non post-weld heat treated welds, while remaining tough and ductile, showed comparatively lower J_Ic, and J-resistance curves at 250°C. This drop in toughness is possibly due to a dynamic strain aging mechanism evidenced by serrated load-displacement curves. The fracture toughness of non post-weld heat treated welds increased significantly after a comparable post-weld heat treatment.The test procedure was validated by comparing three test results against independent tests conducted by Materials Engineering Associates (MEA) of Lanham, Maryland. The J_Ic and J-resistance curves obtained by Ontario Hydro and MEA were comparable.     

Effect of loading rate upon fracture behavior of ferritic steel under large scaled yielding

Dynamic loading to ferromagnetic materials and large scaled yielding result in peak or valley and non-linear curve, respectively, on the Direct Current Potential Drop (DCPD) versus Crack Opening Displacement (COD) plots, which make it difficult to determine the crack initiation point. In this work high intensity of current up to 100 A was applied to the specimens of SA106Gr.C ferritic steel and the crack growth behavior was directly monitored by a high speed camera to obtain the crack initiation point. The effects of loading rate up to 1200 mm min⁻¹ upon the fracture resistance were explored. As the results, it has been shown that, although no substantial difference was seen in the load–COD plots, the crack initiation and then J_i and J–R curve were quite sensitive to the loading rate. That is, under the loading rate of 300 mm min⁻¹ the material showed the worst fracture resistance than under static loading and even under the higher loading rates of 600 and 1200 mm min⁻¹. Also applying the high speed camera and high current source have been proved to be an effective way to find out the accurate crack initiation point and to compensate the pulse of DCPD due to the ferromagnetic effect.     

Failure assessment of RPV nozzle under loss of coolant accident

In case of a postulated loss of coolant accident (LOCA) of a reactor pressure vessel (RPV), the nozzle region experiences higher stresses and lower temperatures than the remaining part of the RPV. Thus, the nozzle is to be considered in the RPV safety assessment. For a LOCA event, three-dimensional elastic–plastic finite element calculations of stresses and strains in the intact RPV were performed. Using the substructure technique, fracture mechanics analyses were then carried out for several postulated cracks in the nozzle corner and in the circumferential weld below the nozzle. For different crack geometries and locations, the J-integral and the stress intensity factor were calculated as functions of the crack tip temperature. Based on the K_IC-reference curve and the J_R curve, both brittle and ductile instability of the postulated cracks were excluded. In order to reduce the expenses of three-dimensional finite element analyses for various crack geometries, an analytical procedure for calculating stress intensity factors of subclad cracks in cylindrical components was extended for cracks in the nozzle corner.     

Fracture initiation and propagation in a cylindrical shell containing an initial surface flaw

The paper deals with the problem of fracture initiation, propagation, and arrest in a pressurized cylindrical vessel which contains an initial surface flaw. It is assumed that the flaw has the most unfavorable geometry and orientation, namely, it is a relatively long part-through axial crack.First we consider the problem of a crack which is sufficiently ‘shallow’ so that the plastic deformations are confined to the neighborhood of the crack border and part of the net section near the inner wall is still elastic. The plasticity-corrected stress intensity factor obtained from this analysis is the controlling load factor in failure considerations related to fatigue crack propagation, stress corrosion cracking, and static fracture (with the use of fracture toughness, COD, or a K_R curvetype failure criterion).The problem of relatively deep crack with fully-yielded net ligament is then considered. Plastic deformations are also assumed to spread around the crack ends through the entire wall thickness. A perfectly plastic strip model (with an eight order shell theory) is used to calculate the plastic zone size and the crack opening displacement along the crack border. Previous studies indicate that for the analysis of the type of stable and subsequent unstable crack propagation problems under consideration, the crack opening displacement δ is a more suitable load factor than the stress intensity factor K, or the crack extension force G. Thus, in this paper a ‘crack opening stretch’ type material characterization will be used.After the rupture of the net ligament under the crack, the axial crack propagation is accompanied by the depressurization of the vessel caused by leakage. From this point on the fracture problem is coupled with the related fluid mechanics or gas dynamics problem where the primary unknowns are the pressure and the crack length as functions of time. In the present study it is assumed that the volume of the vessel is finite and the crack propagation is quasi-static (this assumption, which is necessary to keep the problem within manageable proportions, is justified by the relatively low crack velocities, i.e. v_c < 0.25 c₂, c₂ being the shear wave velocity).     

The Weibull stress parameters calibration upon the toughness scaling model between cracks having different constraint

An extensive investigation has been carried out on the sensitivity parameters determination describing the fracture behaviour of the body with a crack with respect to the character change of the true stress–strain curve with the dominant region of Lueders deformation. This paper presents the consideration on the change judgement of the J-integral and the constraint as the base parameters of two-parameter fracture mechanics. The Weibull stress model for cleavage fracture originally proposed by Beremin group requires calibration of two micromechanics parameters (m, σ_u). The Weibull stress σ_w seems to be a parameter for the prediction of cleavage failure of cracked bodies and the study is focused on the assessment of the effects of constraint loss on the cleavage fracture toughness (J_c). To quantify the effects of the constraint variation on the cleavage fracture toughness the form of the toughness scaling model based on the Weibull stress σ_w is investigated. Local material parameters have been calculated from Beremin approach and the calibration is performed on various approaches. Methods are based on the weakest link assumption and the incremental fracture probability, which depends not only on the maximum principal stress, but also on the equivalent plastic strain. The fracture resistance has been assessed using data from static tests of three point bend specimens.     

The effect of irradiation on the toughness of pressurized water reactor vessel steels under different service conditions

The different toughness tests performed on two pressure vessel steels with very different upper shelves served to make a number of observations concerning the shifts in the transition temperature due to the effect of irradiation, as well as changes in toughness with temperature in the ductile region.With respect to shifts in the transition temperature, the following was observed: the shift obtained with precracked charpy test specimens was narrower than with the others; the shift obtained with charpy V impact tests was substantially equal to that obtained with CT test specimens — wider in the case of steel A, but slightly narrower in that of steel H.With respect to toughness values in the ductile region: the toughness values obtained using precracked charpy test specimens are significantly higher than those obtained with CT test specimens for static tests; 25and 12.5 mm thick CT test specimens display comparable variations in J_1C and dJ/da, but with wide scattering; the effect of irradiation, if any, is of the same order of magnitude as the scattering of the results — however, a test temperature effect is observed; the variation in toughness with temperature is determined more easily by considering a J value corresponding to a stable crack propagation of 1 mm, so that ; this variation of J_Δal with temperature is substantially the same for both steels, or about −30% at 70 or 80°C, and −50% at 290°C.