scholarly journals On the Calculation of Stress Intensity Factors and J-Integrals Using the Submodeling Technique

2010 ◽  
Vol 132 (4) ◽  
Author(s):  
Eduard Marenić ◽  
Ivica Skozrit ◽  
Zdenko Tonković

In the present paper, calculations of the stress intensity factor (SIF) in the linear-elastic range and the J-integral in the elastoplastic domain of cracked structural components are performed by using the shell-to-solid submodeling technique to improve both the computational efficiency and accuracy. In order to validate the submodeling technique, several numerical examples are analyzed. The influence of the choice of the submodel size on the SIF and the J-integral results is investigated. Detailed finite element solutions for elastic and fully plastic J-integral values are obtained for an axially cracked thick-walled pipe under internal pressure. These values are then combined, using the General Electric/Electric Power Research Institute method and the reference stress method, to obtain approximate values of the J-integral at all load levels up to the limit load. The newly developed analytical approximation of the reference pressure for thick-walled pipes with external axial surface cracks is applicable to a wide range of crack dimensions.

2011 ◽  
Vol 133 (6) ◽  
Author(s):  
Masayuki Kamaya

When multiple cracks approach one another, the stress intensity factor and J-integral value change due to the interaction of the stress field. Since the changes in these parameters are not always conservative in structural reliability evaluations, the interaction between multiple cracks should be taken into account. Section XI of the ASME Boiler and Pressure Vessel Code provides a flaw characterization rule for interacting multiple cracks. In Section XI, adjacent cracks are replaced with a coalesced single crack when the distance between the cracks is less than half of the crack depth. However, the criterion for the offset distance is given as an absolute value, although the magnitude of the interaction depends on the crack size. In the current study, an alternative criterion for the offset distance was examined. Linear-elastic and elastic–plastic analyses were performed for interacting semicircular and semi-elliptical surface cracks by the finite element method under a tensile or bending load. The change in the stress intensity factors and J-integral values due to the relative spacing of cracks was investigated. Based on the relationship between the magnitude of the interaction and the relative position of the cracks, the allowable ctriterion for the offset distance was discussed.


Author(s):  
G. Wilkowski ◽  
S. Kalyanam ◽  
Y. Hioe ◽  
F. W. Brust ◽  
S. Pothana ◽  
...  

Abstract Work published for the first time at the ASME PVP 2017 conference showed that when on the upper-shelf, the toughness measured directly from surface-cracked pipe tests decreased as the flaw depth increased. A similar trend existed in SENT tests. Initially it was found that this flaw depth sensitivity of the toughness occurred for a very tough material like TP304 stainless steel. The significance of that result was that even for a material where limit-load was thought to exist, as the flaw depth increased the toughness dropped appreciably, and the failure analysis mode changed from limit-load to elastic-plastic fracture. Experimentally, this made sense because it explained the observed phenomena of load-controlled leak-versus-break behavior for circumferential surface-cracked pipes (as will be shown for several pipe tests), but that LBB behavior is not predictable from circumferential flaw limit-load analysis. Furthermore, the flaw depth effect on toughness also exists for axial surface cracks and even in flat plates with surface cracks. For axial surface cracks the implication was that the long-used empirical surface-crack bulging factor from Maxey/Kiefner (incorporated in many international codes and standards) actually incorporated both the bulging factor and the toughness changes with flaw depth. Because of the change in toughness with flaw depth, when using detailed finite-element fracture analyses for the crack-driving force it is possible to have more error in the failure stress predictions if a constant toughness is assumed for all surface-flaw depths. In fact, in another paper in the ASME 2019 PVP conference it will be shown that the toughness in a wrought TP304 elbow at crack initiation of a circumferential surface crack that was 68% of the thickness was about 1/3rd of the toughness from a standard 1T CT specimen made from the same material. Those results will also be reviewed. Similar results of toughness decreasing with flaw depth in surface-cracked pipes and SENT specimens for various materials over a large range of strain-hardening behavior will show the toughness decrease trend with flaw depth is consistent. To understand these trends more theoretically, 3D FE analyses were also conducted for one initial set of TP304 SENT specimens with a wide range of a/w values (0.3 < a/w < 0.9). The initiation toughness decreased by a factor of 5 to 6 as the crack depth increased; however, the Q value coinciding to the load at the start of ductile tearing was constant for the wide range of a/W values. Q at the start of ductile tearing in the SENT (Qi) was more consistent at normalized distances from the crack tip, rσo/J that were in the range from 0.25 to 1.5 rather than just the popularly considered rσo/J = 2. Hence, by having one SENT test result with a single a/W value, the Ji value for any other a/W can then be calculated. This is consistent with the experimental trends to date, but unfortunately Ji was found to be not proportional to the Q values as is conventionally assumed by many researchers at this time.


1997 ◽  
Vol 32 (3) ◽  
pp. 229-236 ◽  
Author(s):  
X B Lin ◽  
R A Smith

Stress intensity factors for semi-elliptical surface cracks located at the centre of a semicircular edge notch in a finite thickness plate subjected to a remote tensile load are presented in a tabulated format. A wide range of geometry ratios are considered. They are all combinations of the following ratios: the ratio of crack surface half-length to plate half-thickness, c/t = 0.2, 0.4, 0.6, 0.8 and 0.95; the ratio of crack depth to surface half-length, a/c = 0.2, 0.4, 0.6, 0.8 and 1; and the ratio of notch radius to plate half-thickness, r/t = 0.5, 1, 2 and 3. Both the quarter-point displacement and J.-integral methods based on three-dimensional finite element analyses were employed for the calculation of stress intensity factors. The calculation accuracy was studied by analysing the J.-integral path independence and comparing stress intensity factor results with other solutions available in the literature.


Author(s):  
Greg Thorwald ◽  
Pedro Vargas

The reference stress for axial (longitudinal) surface cracks in cylinders is compared using equations from the 2016 API 579-1/ASME FFS-1 and BS 7910:2013 engineering standards, and by using J-integral values from elastic-plastic Finite Element Analysis of three-dimensional crack meshes to compute crack front reference stress. The cylinder axial surface crack reference stress solutions from the two standards differ, and further examination and comparison is desired. To evaluate if a crack is unstable and may cause catastrophic structural failure, the Failure Assessment Diagram method provides an evaluation using two ratios: brittle fracture and plastic collapse. The FAD vertical axis gives the Kr stress intensity to toughness ratio, and the FAD horizontal axis gives the Lr reference stress to yield strength ratio. The details of the FAD method are described in both standards, along with stress intensity and reference stress solutions for various geometries and crack shapes. Since the cylinder axial surface crack reference stress solutions from API 579 and BS 7910 differ, J-integral values are used to compute reference stress trends that provide additional insight and reveal if there is agreement with one or the other or neither standard. Computing reference stress from crack front J-integral results is described in API 579 Annex 9G Section 9G.4. A 3D crack mesh is created for each crack and cylinder size. Along the crack front the focused mesh pattern uses initially coincident groups of nodes at each crack front position. The group of nodes at each location on the crack front are initially coincident and can separate to help model the blunting at the crack front as the loading increases and local plasticity occurs. Post processing calculations use the J-integral versus load trend and the material specific Kr at Lr = 1 ratio to determine the reference stress geometry factor. The reference stress is computed at each crack front node to find the maximum crack front reference stress value for comparison to the engineering standards’ reference stress solutions. A range of surface crack sizes in thin to thick wall cylinders with internal pressure are used to examine reference stress trends. Standard pipe sizes and typical pipeline steel material is used in the analysis. The difference in reference stress solutions was found during an engineering critical assessment, so the J-integral approach was used to improve the solution to reduce conservatism and allow the component to remain in service.


2011 ◽  
Vol 52-54 ◽  
pp. 43-48 ◽  
Author(s):  
Al Emran Ismail ◽  
Ahmad Kamal Ariffin ◽  
Shahrum Abdullah ◽  
Mariyam Jameelah Ghazali ◽  
Ruslizam Daud

This paper presents a non-linear numerical investigation of surface cracks in round bars under bending moment by using ANSYS finite element analysis (FEA). Due to the symmetrical analysis, only quarter finite element (FE) model was constructed and special attention was given at the crack tip of the cracks. The surface cracks were characterized by the dimensionless crack aspect ratio, a/b = 0.6, 0.8, 1.0 and 1.2, while the dimensionless relative crack depth, a/D = 0.1, 0.2 and 0.3. The square-root singularity of stresses and strains was modeled by shifting the mid-point nodes to the quarter-point locations close to the crack tip. The proposed model was validated with the existing model before any further analysis. The elastic-plastic analysis under remotely applied bending moment was assumed to follow the Ramberg-Osgood relation with n = 5 and 10. J values were determined for all positions along the crack front and then, the limit load was predicted using the J values obtained from FEA through the reference stress method.


Author(s):  
Shinji Yoshida ◽  
Hideo Machida

This paper describes applicability of the 2 parameter assessment method using a reference stress method from the viewpoint of reliability. The applicability of the reference stress method was examined comparing both the GE-EPRI method. As a result, J-integral and limit load at the time of fracture evaluated by the reference stress method is almost equivalent to that by the GE-EPRI method. Furthermore, the partial safety factor (PSF) evaluated by reliability assessment has little difference between two methods, and the required safety factor is enveloped by the safety factor for Service Level-A and B defined in fitness for service (FFS) codes. These results show that of the reference stress method is applicable for J-integral calculation in fracture assessment.


Author(s):  
Bruno Michel ◽  
Jean-Philippe Sermage ◽  
Philippe Gilles ◽  
Bruno Barthelet ◽  
Patrick Le Delliou

The RSE-M Code [1] provides rules and requirements for in-service inspection of French Pressurized Water Reactor power plant components. Non mandatory guidance is given in the Code for analytical flaw evaluation in a wide range of situations. In Appendix 5.4 of the Code, influence coefficients are provided to calculate stress intensity factors in pipes and shells containing semi-elliptical surface defects. The J assessment method is based on the reference stress concept with two options for reference loads evaluation: “CEP elastic plastic stress” and “CLC modified limit load”. In this paper recent advances concerning J assessment under mechanical loading for a crack located in a pipe-elbow junction are presented. Reference stress evaluation with “CLC” option is developed and mechanical foundations of the equation of large scale yielding under complex loading (pressure, in-plane and out-of-plane bending) are presented.


Sign in / Sign up

Export Citation Format

Share Document