Analysis of Residual Stresses and Assessment of Postulated Cracks in a Core Shroud Weld

Author(s):  
Igor Varfolomeyev ◽  
Dieter Siegele ◽  
Dieter Beukelmann

In order to assess postulated cracks in weldments of a BWR core shroud residual stresses are calculated by simulating the welding process. In the numerical analysis, weld metal deposition and the sequence of weld passes follow the manufacture protocol. The calculations are performed using the finite element program ABAQUS and a material model with kinematic nonlinear hardening. Calculations of the crack driving parameter, the stress intensity factor, are carried out for postulated circumferential cracks using a numerical procedure, as well as by applying a weight function solution specially developed for cracks in a thin-walled cylinder. The results give rise to a discussion on the validity of linear elastic fracture mechanics for assessing defects in weldments. Additionally, for a complete circumferential crack the trend in the stress intensity factor is studied when the crack depth approaches the full wall thickness.

2005 ◽  
Vol 297-300 ◽  
pp. 743-749
Author(s):  
Min Koo Han ◽  
Mamidala Ramulu

Fatigue crack propagation life of weld toe crack through residual stress field was estimated using Elber's crack closure concept. Propagation of weld toe crack is heavily influenced by residual stresses caused by the welding process, so it is essential to take into account the effect of residual stresses on the propagation life of a weld toe crack. Fatigue cracks at transverse and longitudinal weld toe was studied, these two cases represent the typical weld joints in ship structures. Numerical and experimental studies are performed for both cases. Residual stresses near the welding area were estimated through a nonlinear thermo-elasto-plastic finite element method and the residual stress intensity factor with Glinka's weight function method. Effective stress intensity factor was calculated using the Newman-Forman-de Koning-Henriksen equation, which is based on the Dugdale strip yield model in estimating the crack closure level, U, at different stress ratios. Calculated crack propagation life coincided well with experimental results.


1989 ◽  
Vol 111 (1) ◽  
pp. 61-67 ◽  
Author(s):  
F. Erdogan

The problem of slow crack growth under residual stresses and externally applied loads in plates is considered. Even though the technique developed to treat the problem is quite general, in the solution given it is assumed that the plate contains a surface crack and the residual stresses are compressive near and at the surfaces and tensile in the interior. The crack would start growing subcritically when the stress intensity factor exceeds a threshold value. Initially the crack faces near the plate surface would remain closed. A crack-contact problem would, therefore, have to be solved to calculate the stress intensity factor. Depending on the relative magnitudes of the residual and applied stresses and the threshold and critical stress intensity factors, the subcritically growing crack would either be arrested or become unstable. The problem is solved and examples showing the time to crack arrest or failure are discussed.


Author(s):  
Zheng Liu ◽  
Xu Chen ◽  
Xin Wang

In the present paper, three-dimensional clamped SENT specimens, which is one of the most widely used low-constraint and less-conservative specimen, are analyzed by using a crack compliance analysis approach and extensive finite element analysis. Considering the test standard (BS8571) recommended specimen sizes, the daylight to width ratio, H/W, is 10.0, the relative crack depth, a/W, is varied by 0.2, 0.3, 0.4, 0.5 or 0.6 and the relative plate thickness, B/W, is chosen by 1.0, 2.0 or 4.0, respectively. Complete solutions of fracture mechanics parameters, including stress intensity factor (K), in-plane T-stress (T11) and out-of-plane T-stress (T33) are calculated, and the results obtained from above two methods have a good agreement. Moreover, the combination of the effects of a/W and B/W on the stress intensity factor K, T11 and T33 stress are thus illustrated.


Materials ◽  
2019 ◽  
Vol 12 (7) ◽  
pp. 1042 ◽  
Author(s):  
Gabriel Coêlho ◽  
Antonio Silva ◽  
Marco Santos ◽  
Antonio Lima ◽  
Neilor Santos

The purpose of this research is to compare both British standard BS 7910 (2013) and American standard API 579/ASME FFS-1 (2016) stress intensity factor (SIF) solutions by considering a series of semielliptical surface cracks located in the external surface of a pressurized hollow cylinder in the axial direction. Finite element analysis was used as a comparison basis for both standards’ SIF results. The solution from the British standard provided consistent results compared to Finite Element (FE) results for crack depth not much higher than half the thickness in the deepest and surface-breaking points. Above those limits, the British standard’s solutions diverged quite a lot from the American standard, whose results followed FE values for every crack depth/thickness ratio tested with a maximum percentage difference of 1.83%.


Author(s):  
Patrick Le Delliou ◽  
Stéphane Chapuliot

Analytical evaluation procedures for determining the acceptability of flaw detected during in-service inspection of nuclear power plant components are provided in Appendix 5.4 of the French RSE-M Code. Linear elastic fracture mechanics based evaluation procedures require calculation of the stress intensity factor (SIF). In Appendix 5.4 of the RSE-M Code, influence coefficients needed to compute the SIF are provided for a wide range of surface axial or circumferential flaws in cylinders, the through-wall stress field being represented by a cubic equation. On the other hand, Appendix C of API 579-1 FFS procedure provides also a very complete set of influence coefficients. The paper presents the comparison of the influence coefficients from both documents, focused on axial ID semi-elliptical surface flaws in cylinders. The cylinder and crack geometries are represented by three ratios: Ri/t, a/t, and a/c, where Ri, t, a, and c are respectively the inner radius, the wall thickness, the crack depth and one-half of the crack length. The solutions for the coefficients G0 and G1 at the deepest point and at the surface point are investigated. At the deepest point, the agreement between the solutions is good, the relative difference being lower than 2 %, except for the plate (Ri/t = ∞) at a/c = 0.125 and 0.0625 and a/t = 0.8 (around 5 %). At the surface point, the agreement between both solutions is not so good. At this point, the relative differences depend strongly on the a/c ratio, being larger for elongated cracks (with low a/c ratios). However, it must be recalled that the absolute values of the coefficients are low at the surface point for elongated cracks, and that for these cracks the critical point regarding the stress intensity factor is the deepest point.


Author(s):  
Igor Orynyak ◽  
Andrii Oryniak

The consideration of a geometrical nonlinearity is a common practice for the thin-walled structures. The relevance here are two well-known cases treated in ASME codes. First one is accounting for reduction of the pipe bends flexibility due to the inner pressure. The second one is the retarded increasing (and subsequent saturation) of additional local bending stress with increasing of inner pressure in a pipe with initial cross section form distortion. In both cases the rerounding effect and decreasing of local flexibilities take place. The crack can be treated as the concentrated flexibility and it is quite natural to expect that the stress intensity factor should grow nonlinearly with applied load. Two cases of SIF calculation for 1-D long axial surface crack in a pipe loaded by inner pressure are considered here: a) cross section has an ideal circular form: b) the form has a small distortion and crack is located in the place of maximal additional bending stresses. The theoretical analysis is based on: a) the well known crack compliance method [1] and b) analytical linearized solution obtained for deformation of the curved beam in case of action of fixed circumferential stress due to pressure written in the form convenient for transfer matrix method application. It was shown that for moderately deep crack (crack depth to the wall thickness ratio is 0.5 and bigger) and typical dimensions of pipes used for oil and gas transportation (radius to thickness ratio is 25–40) and loading which can reach up to 200 to 300 MPa, the effect investigated can be quite noticeable and can lead to 5–15 percent reduction of calculated SIF as compared with linear calculation. The analytical results are supported by nonlinear FEM calculation.


Author(s):  
Mayumi Ochi ◽  
Kiminobu Hojo ◽  
Itaru Muroya ◽  
Kazuo Ogawa

Alloy 600 weld joints have potential for primary water stress corrosion cracks (PWSCC). At the present time it has been understood that PWSCC generates and propagates in the Alloy 600 base metal and the Alloy 600 weld metal and there has been no observation of cracking the stainless and the low alloy steel. For the life time evaluation of the pipes or components the crack extension analysis is required. To perform the axial crack extension analysis the stress intensity database or estimation equation corresponding to the extension crack shape is needed. From the PWSCC extension nature mentioned above, stress intensity factors of the conventional handbooks are not suitable because most of them assume a semi-elliptical crack and the maximum aspect ratio crack depth/crack half length is one (The evaluation in this paper had been performed before API 579-1/ASME FFS was published). Normally, with the advance of crack extension in the thickness direction at the weld joint, the crack aspect ratio exceeds one and the K-value of the conventional handbook can not be applied. Even if those equations are applied, the result would be overestimated. In this paper, considering characteristics of PWSCC’s extension behavior in the welding material, the axial crack was modeled in the FE model as a rectangular shape and the stress intensity factors at the deepest point were calculated with change of crack depth. From the database of the stress intensity factors, the simplified equation of stress intensity factor with parameter of radius/thickness and thickness/weld width was proposed.


Author(s):  
V Sura ◽  
S Mahadevan

Shattered rim cracking, propagation of a subsurface crack parallel to the tread surface, is one of the dominant railroad wheel failure types observed in North America. This crack initiation and propagation life depends on several factors, such as wheel rim thickness, wheel load, residual stresses in the rim, and the size and location of material defects in the rim. This article investigates the effect of the above-mentioned parameters on shattered rim cracking, using finite element analysis and fracture mechanics. This cracking is modelled using a three-dimensional, multiresolution, elastic–plastic finite element model of a railroad wheel. Material defects are modelled as mathematically sharp cracks. Rolling contact loading is simulated by applying the wheel load on the tread surface over a Hertzian contact area. The equivalent stress intensity factor ranges at the subsurface crack tips are estimated using uni-modal stress intensity factors obtained from the finite element analysis and a mixed-mode crack growth model. The residual stress and wheel wear effects are also included in modelling shattered rim cracking. The analysis results show that the sensitive depth below the tread surface for shattered rim cracking ranges from 19.05 to 22.23 mm, which is in good agreement with field observations. The relationship of the equivalent stress intensity factor (Δ K eq) at the crack tip to the load magnitude is observed to be approximately linear. The analysis results show that the equivalent stress intensity factor (Δ K eq) at the crack tip depends significantly on the residual stress state in the wheel. Consideration of as-manufactured residual stresses decreases the Δ K eq at the crack tip by about 40 per cent compared to that of no residual stress state, whereas consideration of service-induced residual stresses increases the Δ K eq at the crack tip by about 50 per cent compared to that of as-manufactured residual stress state. In summary, the methodology developed in this article can help to predict whether a shattered rim crack will propagate for a given set of parameters, such as load magnitude, rim thickness, crack size, crack location, and residual stress state.


1964 ◽  
Vol 86 (4) ◽  
pp. 681-684 ◽  
Author(s):  
A. S. Kobayashi ◽  
R. D. Cherepy ◽  
W. C. Kinsel

The advantages of the complex variable method are combined with the numerical procedure of collocation for estimating the stress intensity factors in finite, cracked plates subjected to in-plane loadings. In this approach, the complex stress functions for an infinite plate problem are modified to meet the boundary conditions for a finite plate with identical crack configuration. This procedure produces a system of linear equations which can be programmed readily on high-speed computers. The procedure is used to find the elastic stress intensity factor at the crack tip in a centrally notched plate in uniaxial tension. The resulting values are nearly identical to the stress intensity values determined analytically by the theory of elasticity. This numerical procedure should be useful for designers and analysts working in the fields of fracture mechanics and fail-safe concepts.


Author(s):  
Adam Toft ◽  
David Beardsmore ◽  
Colin Madew ◽  
Huego Teng ◽  
Mark Jackson

Within the UK nuclear industry the assessment of fracture in pressurised components is often carried out using procedures to calculate the margin of safety between a lower-bound fracture toughness and the crack driving force. Determination of the crack driving force usually requires the calculation of elastic stress intensity factor solutions for primary loads and secondary loads arising from weld residual stresses and/or thermal stresses. Within established UK assessment procedures weight function solutions are available which allow the stress intensity factors to be calculated from the through-wall opening-mode stress distribution in an uncracked component. These weight-function solutions are generally based on models where either no boundary condition is applied, or where one is applied at a distance either side of the crack plane that is very long compared with the crack size and wall thickness. Such solutions do not take into account any reduction in the stress field that might occur as the distance from the crack faces increases. Weld residual stress fields may often be expected to reduce in this manner. A separate, earlier study has shown that the stress intensity factor for a cracked plate loaded in displacement control decreases substantially as the loading plane is moved closer to the crack plane. It would therefore be expected that a similar reduction in stress intensity factor would be obtained for a residual stress analysis when displacement boundary conditions are imposed at a distance relatively close to the crack plane. This paper describes an investigation of the differences, particularly in terms of a reduction in calculated stress intensity factor, which may arise from application of displacement controlled stress intensity factor solutions, as compared with load controlled solutions, when considering weld residual stresses. Consideration is also given as to how new displacement controlled stress intensity factor solutions could be developed by modification of existing load controlled solutions.


Sign in / Sign up

Export Citation Format

Share Document