0
Research Papers: Materials and Fabrication

Simulation of Stress Corrosion Cracking in In-Core Monitor Housing of Nuclear Power Plant

[+] Author and Article Information
Yuichi Shintaku

Tokyo University of Science,
2641 Yamazaki, Noda,
Chiba 278-8510, Japan
e-mail: shin@me.noda.tus.ac.jp

Fuminori Iwamatsu

Hitachi Research Laboratory,
Hitachi, Ltd.,
11 Saiwai-cho 3-chome, Hitachi,
Ibaraki 317-8511Japan
e-mail: fuminori.iwamatsu.vt@hitachi.com

Kazuhiro Suga

Tokyo University of Science,
2641 Yamazaki, Noda,
Chiba 278-8510, Japan
e-mail: ksuga@rs.noda.tus.ac.jp

Yoshitaka Wada

Kinki University,
341, Kowakae,
Higashiosaka 577-8502, Japan
e-mail: wada@mech.kindai.ac.jp

Masanori Kikuchi

Tokyo University of Science,
2641 Yamazaki, Noda,
Chiba 278-8510, Japan
e-mail: kik@me.noda.tus.ac.jp

Contributed by the Pressure Vessel and Piping Division of ASME for publication in the JOURNAL OF PRESSURE VESSEL TECHNOLOGY. Manuscript received March 22, 2013; final manuscript received October 1, 2014; published online February 20, 2015. Assoc. Editor: Xian-Kui Zhu.

J. Pressure Vessel Technol 137(4), 041401 (Aug 01, 2015) (13 pages) Paper No: PVT-13-1054; doi: 10.1115/1.4028735 History: Received March 22, 2013; Revised October 01, 2014; Online February 20, 2015

In the in-core monitor (ICM) housing of a reactor pressure vessel (RPV), residual stress has been widely reported to cause stress corrosion cracking (SCC) damage in the weld heat-affected zone. For this reason, it is important to evaluate the crack growth conservatively, and with high confidence to demonstrate fitness for service. This paper presents crack growth simulations in an ICM housing, which is welded at two different angles to the RPV. One weld angle is at the bottom of the RPV, and the welding area of the ICM housing is axisymmetric. The other angle is at the curved position of the RPV, and the weld area of the ICM housing is asymmetric. In these weld areas, crack growth behavior is estimated by superposed-finite element method (S-FEM), which allows generation of a global finite model and a detailed local mesh representing the crack independently. In the axisymmetric weld area, axial, slant and circumferential surface cracks are assumed at two locations where the residual stress fields are different from each other: one is isotropic and the other is circumferential. It is shown that crack growth behaviors are different under different residual stress fields. The results of S-FEM are compared with those of the influence function method (IFM), which assumes that an elliptical crack shape exists in a plate. It is shown that the IFM result is conservative compared to that of S-FEM. Next, an axial surface crack is assumed at the uphill, downhill, and midhill asymmetric weld areas. The midhill crack growth behavior is different from the uphill and downhill behaviors. Finally, two surface cracks are simulated in the asymmetric weld area and two initial crack arrangements are assumed. These results show the differences of the crack interaction and the crack growth process.

Copyright © 2015 by ASME
Your Session has timed out. Please sign back in to continue.

References

Nuclear and Industrial Safety Agency, 2004, Minister of Economy, Trade and Industry, Subcommittee on Stress Test for Atomic Power Plant (No. 10): Documents Distributed 10-3 Stress Test on Core Shroud and Primary Loop Recirculation System (in Japanese).
JSME S NA1-2008, 2008, “Codes for Nuclear Power Generation Facilities: Rule of Fitness-for-Service for Nuclear Power Plants,” Tokyo, Japan.
The Japan Society of Mechanical Engineers, 2004, “Rules on Fitness-for-Service for Nuclear Power Plants,” Report No. JSME S NA1-2004.
Shiratori, M., Miyoshi, T., and Tanikawa, K., 1985, “Analysis of Stress Intensity Factors for Surface Cracks Subjected to Arbitrarily Distributed Surface Stresses,” Trans. JSME, Ser. A, 51(467), pp. 1828–1833 (in Japanese). [CrossRef]
Buckner, H. F., 1958, “The Propagation of Cracks and the Energy of Elastic Deformation,” Trans. ASME, 80E, pp. 1225–1230.
Matsushita, H., Omata, S., Matsuda, H., and Shiratori, M., 2004, “Fatigue Surface Crack Analysis Software SCANP,” Nippon Kaiji Kyoukai, 266, pp. 18–35 (in Japanese).
Iwamatsu, F., Miyazaki, K., and Shiratori, M., 2009, “Development of Evaluation Method of Stress Intensity Factor and Fatigue Crack Growth Behavior of Surface Crack Under Arbitrarily Stress Distribution by Using Influence Function Method,” Trans. JSME, Ser. A, 77(782), pp. 1613–1624 (in Japanese). [CrossRef]
Belytchko, T., Lu, Y. Y., and Gu, L., 1995, “Element-Free Galerkin Method,” Int. J. Numer. Methods Eng., 37(2), pp. 229–256. [CrossRef]
Belytchko, T., Lu, Y. Y., and Gu, L., 1955, “Crack Propagation by Element-Free Galerkin Methods,” Eng. Fract. Mech., 51(2), pp. 295–315. [CrossRef]
Belytschko, T., and Black, T., 1999, “Elastic Crack Growth in Finite Elements With Minimal Remeshing,” Int. J. Numer. Methods Eng., 45(5), pp. 602–620. [CrossRef]
Sakumasu, K., and Nagashima, T., 2007, “Three-Dimensional Crack Propagation Analysis Using X-FEM,” JSME Annual Meeting, pp. 21–22.
Fish, J., Markolefas, S., Guttal, R., and Nayak, P., 1994, “On Adaptive Multilevel Superposition of Finite Element Meshes for Linear Elastostatics,” Appl. Numer. Math., 14(1–3), pp. 135–164. [CrossRef]
dell'Erba, D. N., and Aliabadi, M. H., 2000, “Three-Dimensional Thermo-Mechanical Fatigue Crack Growth Using BEM,” Int. J. Fatigue, 22(4), pp. 261–273. [CrossRef]
Kikuchi, M., Wada, Y., Shimizu, Y., and Li, Y., 2011, “Crack Growth Analysis in a Weld-Heat-Affected Zone Using S-Version FEM,” Int. J. Press. Vessels Pip., 90–91, pp. 2–8. [CrossRef]
Kikuchi, M., Mattireymu, M., and Sano, H., 2010, “Fatigue Crack Growth Simulation Using S-Version FEM: 3rd Report, Fatigue of 3D. Surface Crack,” Trans. JSME, Ser. A, 75(755), pp. 918–924 (in Japanese).
Kikuchi, M., Wada, Y., Suyama, H., and Shimizu, Y., 2009, “Crack Growth Analysis in Weld Heat Affected Zone Using S-Version FEM,” Trans. JSME, Ser. A, 75(758), pp. 1381–1386 (in Japanese).
Okada, H., Higashi, M., Kikuchi, M., Fukui, Y., and Kumazawa, N., 2005, “Three-Dimensional Virtual Crack Closure-Integral Method (VCCM) With Skewed and Non-Symmetric Mesh Arrangement at the Crack Front,” Eng. Fract. Mech., 72(11), pp. 1717–1737. [CrossRef]
Erdogan, F., and Shi, G. C., 1963, “On the Crack Extension in Plates Under Plane Loading and Transverse Shear,” ASME J. Basic Eng., 85(4), pp. 519–525. [CrossRef]
Sih, G. C., 1990, Mechanics of Fracture Initiation and Propagation, Kluwer Academic, Boston.
Pook, L. P., 2000, Linear Elastic Fracture Mechanics for Engineers: Theory and Application, WIT Press.
Schöllmann, M., Kullmer, G., Fulland, M., and Richard, H. A., 2001, “A New Criterion for 3D Crack Growth Under Mixed-Mode (I+II+III) Loading,” Proceedings of the 6th International Conference on Biaxial/Multiaxial Fatigue and Fracture, Lisbon, Portugal, June 25–28, Vol. 2, pp. 589–596.
Richard, H. A., Fulland, M., and Sander, M., 2005, “Theoretical Crack Path Prediction,” Fatigue Fract. Eng. Mater. Struct., 28(1–2), pp. 3–12. [CrossRef]
Mochizuki, M., Enomoto, K., Okamoto, N., Saito, H., and Hayashi, E., 1994, “Study on Production Mechanism of Welding Residual Stress at the Juncture of a Pipe Penetrating a Thick Plate,” Q. J. Jpn. Weld. Soc., 12(4), pp. 561–567 (in Japanese). [CrossRef]
Miyazaki, K., Kanno, S., Mochizuki, M., Hayashi, M., Shiratori, M., and Yu, Q., 1999, “Analysis of Stress Intensity Factor Due to Surface Crack Propagation by Using Influence Function Method and Inherent Strain Analysis in Residual Stress Fields Caused by Welding,” Trans. JSME, Ser. A, 65(636), pp. 1709–1715 (in Japanese). [CrossRef]
Shiratori, M., Miyoshi, T., and Skai, T., 2010, “Analysis of Stress Intensity Factors for Surface Cracks Subjected to Arbitrarily Distributed Surface Stresses (4th Report, Application of Influence Coefficients for the Cracks Originating at the Notches and Welding Joints),” Trans. JSME, Ser. A, 53(755), pp. 918–924 (in Japanese).
Kikuchi, M., Wada, Y., Utsunomiya, A., and Suyama, H., 2009, “Study on Fatigue Crack Growth Criterion: 1st Report, Paris' Law of a Surface Crack Under Pure Mode I Loading,” Trans. JSME, Ser. A, 76(764), pp. 516–522 (in Japanese).

Figures

Grahic Jump Location
Fig. 2

Crack growth direction

Grahic Jump Location
Fig. 3

Geometry of Pressure Vessel Reactor Lower Head

Grahic Jump Location
Fig. 4

Geometry of ICM housing (a) axisymmetric weld area (b) antisymmetric weld area

Grahic Jump Location
Fig. 5

Global mesh representing ICM housing (a) axisymmetric weld area (b) antisymmetric weld area

Grahic Jump Location
Fig. 6

Initial crack shape

Grahic Jump Location
Fig. 7

Definition of crack shape (a) definition of crack length and depth (b) definition for eccentric angle

Grahic Jump Location
Fig. 8

Initial crack in symmetric weld area of ICM housing and distribution of residual stress (a) hoop stress σθ (b) axial stress σz

Grahic Jump Location
Fig. 9

Initial crack shape (a) axial crack (b) slant crack (c) circumferential crack

Grahic Jump Location
Fig. 10

Stress intensity factor KI of axial crack in symmetric weld area (ψ = 0 deg, z = 15 mm) (a) a/t = 0.2 (b) a/t = 0.6

Grahic Jump Location
Fig. 11

Axial crack growing process in symmetric weld area (ψ = 0 deg, z = 15 mm) and distribution of hoop stress

Grahic Jump Location
Fig. 12

Stress intensity factor KI of slant crack in symmetric weld area (z = 15 mm)

Grahic Jump Location
Fig. 13

Normalized stress intensity factor KII/KI of slant crack in symmetric weld area (ψ = 45 deg, z = 15 mm)

Grahic Jump Location
Fig. 14

Normalized stress intensity factor KIII/KI of slant crack in symmetric weld area (ψ = 45 deg, z = 15 mm)

Grahic Jump Location
Fig. 15

Slant crack growing process in symmetric weld area (ψ = 45 deg, z = 15 mm) (a) a/t = 0.2 (b) a/t = 0.3 (c) a/t = 0.6

Grahic Jump Location
Fig. 16

Stress intensity factor KI of slant crack in symmetric weld area (ψ = 45 deg, z = 25 mm)

Grahic Jump Location
Fig. 17

Slant crack growing process in symmetric weld area (ψ = 45 deg, z = 25 mm) (a) a/t = 0.2 (b) a/t = 0.4 (c) a/t = 0.6

Grahic Jump Location
Fig. 18

Circumferential crack growing process in symmetric weld area (ψ = 90 deg, z = 25 mm) (a) a/t = 0.2 (b) a/t = 0.4 (c) a/t = 0.6

Grahic Jump Location
Fig. 19

Relationship between crack speed and initial crack angle (a) comparison of crack depth (b) comparison of crack length

Grahic Jump Location
Fig. 21

Comparison of S-FEM and IFM due to difference in crack growth behavior and geometry of ICM housing (a) comparison of crack depth (b) comparison of crack length

Grahic Jump Location
Fig. 22

Initial crack in antisymmetric weld area of ICM housing and distribution of residual stress

Grahic Jump Location
Fig. 23

Stress intensity factor KI in uphill of antisymmetric weld area

Grahic Jump Location
Fig. 24

Crack growing process in uphill of antisymmetric weld area (a) a /t = 0.2 (b) a/t = 0.6

Grahic Jump Location
Fig. 25

Slant crack growing process in middle hill of antisymmetric weld area (a) a /t = 0.2 (b) a/t = 0.6

Grahic Jump Location
Fig. 26

Initial crack arrangement in symmetric weld area (a) parallel arrangement to weld line (b) vertical arrangement to weld line

Grahic Jump Location
Fig. 27

Crack growing process of two vertical cracks to weld line (a) 0.0 [year] (b) 3.3 [year] (c) 6.8 [year]

Grahic Jump Location
Fig. 28

Stress intensity factor KI of two vertical cracks to weld line (a) 0.0 [year] (b) 3.3 [year] (c) 6.8 [year]

Grahic Jump Location
Fig. 29

Crack growing process of two parallel cracks to weld line (a) 0.0 [year] (b) 6.6 [year] (c) 19.8 [year]

Grahic Jump Location
Fig. 30

Stress intensity factor KI of two parallel cracks to weld line (a) 0.0 [year] (b) 6.6 [year] (c) 19.8 [year]

Grahic Jump Location
Fig. 31

Relationship between crack speed and difference of the initial arrangement of two cracks

Tables

Errata

Discussions

Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In