Cleavage Dynamic Propagation Analysis in a Nuclear Reactor Pressure Vessel Steel Using a High-Speed Camera

Author(s):  
Amaury Bousquet ◽  
Stéphane Marie ◽  
Philippe Bompard

Initiation stage of cracks is considered as a key issue, but more and more component integrity analyses investigate the crack propagation and arrest possibility. This study deals with physical mechanisms of cleavage crack propagation and numerical computations related to brittle fracture. Dynamic effects, involved in unstable cleavage crack propagation, have to be taken into account to properly depict brittle crack propagation, arrest and possible propagation re-initiation events. Experiments were carried out on thin CT specimens made of 16MND5 PWR vessel steel at five temperatures (−150°C, −125°C, −100°C, −75°C, −50°C). In addition to standard crack gages, an innovative experimental technique has been used to determine crack propagation. By the means of developments on the experimental protocol (improvements of isolation and airtightness of the thermal chamber, optimization of the experimental protocol to eliminate ice in the thermal chamber and in order to have a good acquisition quality), use of a high-speed framing camera was made possible to measure crack propagation on a CT mirror polished surface. This optical device, combined with this optimized experimental process, has allowed the study of straight and branching crack paths with high accuracy. The framing camera (520 000 fps up to 1 100 000 fps) has allowed to have a very accurate estimation of crack speed even up to 1000 m.s−1 and also to detect some phases of crack branching during propagation and phases of arrest-re-initiation. Numerical computations, based on X-FEM and combining a local non linear dynamic approach with a RKR type fracture stress criterion, have been performed to depict experimental crack behavior. This paper describes this innovative experimentation and the interpretation by FE calculations and SEM observations associated with quantitative 3D optical microscopy.

Author(s):  
Amaury Bousquet ◽  
Stéphane Marie ◽  
Philippe Bompard

The integrity assessment of Reactor Pressure Vessels, mainly based on crack initiation, can be completed by studying crack propagation and arrest. Whereas engineering approaches do not take into account dynamic effects, these effects are important in unstable cleavage crack propagation, arrest and possible propagation re-initiation events. This study deals with physical mechanisms of cleavage crack propagation and numerical computations related to brittle fracture in the framework of local approach to fracture. Experiments were carried out on thin CT 25 specimens made of 16MND5 PWR vessel steel at five temperatures (−150°C, −125°C, −100°C, −75°C, −50°C). Two kinds of crack path, straight or branching path, were observed. Branching cracks appear for the highest critical loadings at initiation, that increase the elastic stored energy and the effect of plasticity. The elastic-viscoplastic behavior of the ferritic steel was studied up to a strain rate of 104 s−1 and taken into account in the numerical simulations. The eXtended Finite Element Method (X-FEM) was used in CAST3M FE software to model crack propagation. Numerical computations combine a local non linear dynamic approach with a RKR type fracture stress criterion. The different physical micro-mechanisms, involved in cleavage fracture, were examined by the means of SEM fracture surface analyses at different temperatures and strain rates for the two kinds of crack path. The links of the critical fracture stress with both temperature and strain rate for straight crack path as well as analyses of branching crack phenomena were considered by the means of Scanning Electron Microscopy (SEM) fracture surface analyses, 3D quantitative optical microscopy and FE computations in order to aim at a robust physical justification of the propagation model which has already been developed at CEA in the frame of the B. Prabel PhD.


Author(s):  
Xiaoyu Yang ◽  
Stéphane Marie ◽  
Clémentine Jacquemoud

Cleavage crack propagation has been tested for three different geometries of Compact Tensile (CT) specimens: CT25, CT50 and extended CT25 (CT25 with CT50 width) (Figure 3). The experimental results show that the crack paths are straight for CT25 and CT50, but they are unstable and curved for extended CT specimens (Figure 5 to 7). Numerical computation had been performed by extended finite element method (XFEM) in CAST3M software. 2D modeling was used in order to predict the crack path. The analysis was based on a local non-linear dynamic approach with a RKR fracture stress criterion depending on temperature and strain rate. In order to simulate the curvature of the cracks path, a statistical effect was introduced in the model to take into account the spatial distribution of cleavage initiators, which is the characteristic of cleavage fracture. At each step of propagation during the modeling, the direction was randomly chosen, according to a uniform defects distribution. The numerical results show a good agreement with experience. The different crack paths were curved in extended CT25, but remained almost straight in CT25 and CT50 specimens, despite of the instability introduced in the modeling in the propagation direction. These results show that the statistics of micro-defects can induce, jointly with the geometry of specimen, a large scatter of crack propagation paths.


2018 ◽  
Vol 191 ◽  
pp. 486-503 ◽  
Author(s):  
Xiaoyu Yang ◽  
Stéphane Marie ◽  
Clémentine Jacquemoud ◽  
Philippe Bompard

Author(s):  
Xiuqin Chu ◽  
Wenting Guo ◽  
Jun Wang ◽  
Feng Wu ◽  
Yuhuan Luo ◽  
...  

1988 ◽  
Vol 110 (3) ◽  
pp. 408-413 ◽  
Author(s):  
L. J. Ghosn

Crack propagation in a rotating inner raceway of a high-speed roller bearing is analyzed using the boundary integral method. The model consists of an edge plate under plane strain condition upon which varying Hertzian stress fields are superimposed. A multidomain boundary integral equation using quadratic elements was written to determine the stress intensity factors KI and KII at the crack tip for various roller positions. The multidomain formulation allows the two faces of the crack to be modeled in two different subregions making it possible to analyze crack closure when the roller is positioned on or close to the crack line. KI and KII stress intensity factors along any direction were computed. These calculations permit determination of crack growth direction along which the average KI times the alternating KI is maximum.


Author(s):  
B. Prabel ◽  
S. Marie ◽  
A. Combescure

In the frame of analysis of the pressure thermal shock in a PWR RVP and the associated R&D activities, some developments are performed at CEA on the dynamic brittle propagation and crack arrest. This paper presents a PhD work on the modeling of the dynamic brittle crack growth. For the analyses, an important experimental work is performed on different geometries using a French RPV ferritic steel: Compact Tension specimens with different thickness, isothermal rings under compression with different positions of the initial defect to study a mixed mode configuration, and a ring submitted to thermal shock. The first part of this paper details the test conditions and main results. To propose an accurate interpretation of the crack growth, a viscous-elastic-plastic dynamic model is used. The strain rate influence is taken into account based on Cowper-Symond’s law (characterization was made from Split Hopkinson Pressure Bar tests). To model the crack propagation in the Finite Element calculation, eXtended Finite Element Method (X-FEM) is used. The implementation of these specific elements in the CEA F.E. software CAST3M is described in the second part of this paper. This numerical technique avoids re-meshing, because the crack progress is directly incorporated in the degrees of freedom of the elements crossed by the crack. The last part of this paper compares the F.E. predictions to the experimental measurements using different criteria. In particular, we focused on a RKR (Ritchie-Knott-Rice) like criterion using a critical principal stress in the front of the crack tip during the dynamic crack extension. Critical stress is found to depend on crack speed, or equivalently on strain rate. Good results are reported concerning predictive simulations.


1959 ◽  
pp. 345-364
Author(s):  
F. H. Nadig ◽  
J. L. Bohn ◽  
T. Korneff

Author(s):  
Frank Liebold ◽  
Ali A. Heravi ◽  
Oliver Mosig ◽  
Manfred Curbach ◽  
Viktor Mechtcherine ◽  
...  

The determination of crack propagation velocities can provide valuable information for a better understanding of damage processes of concrete. The spatio-temporal analysis of crack patterns developing at a speed of several hundred meters per second is a rather challenging task. In the paper, a photogrammetric procedure for the determination of crack propagation velocities in concrete specimens using high-speed camera image sequences is presented. A cascaded image sequence processing which starts with the computation of displacement vector fields for a dense pattern of points on the specimen’s surface between consecutive time steps of the image sequence chain has been developed. These surface points are triangulated into a mesh, and as representations of cracks, discontinuities in the displacement vector fields are found by a deformation analysis applied to all triangles of the mesh. Connected components of the deformed triangles are computed using region-growing techniques. Then, the crack tips are determined using principal component analysis. The tips are tracked in the image sequence and the velocities between the time stamps of the images are derived. A major advantage of this method as compared to established techniques is in the fact of its allowing for spatio-temporally resolved, full-field measurements rather than point-wise measurements and that information on crack width can be obtained simultaneously. To validate the experimentation, the authors processed image sequences of tests on four compact-tension specimens performed on a split-Hopkinson tension bar. The images were taken by a high-speed camera at a frame rate of 160,000 images per second. By applying to these datasets the image sequence processing procedure as developed, crack propagation velocities of about 800 m/s were determined with a precision in the order of 50 m/s.


Sign in / Sign up

Export Citation Format

Share Document