Document: NRC Regulatory Guide
Document ID: 2704425a-c58a-45c4-93ab-8761721c3e7a
Document Type: regulatory_guide
Title: Evaluation of Reactor Pressure Vessels with Charpy Upper-Shelf Energy Less Than 50 Ft-Lb
Source: NRC Regulatory Guide Division 1
Source URL: https://www.nrc.gov/docs/ML0037/ML003740038.pdf
Revision Date: 2023-06
Chapter: 
Section ID: RG-1.161
CFR Part: 
CFR Title: 

Content:
(3)-56*CC(4)+1 1 *CC(5))*(2.0/3 0) A4 = (-5*CC(1)+18*CC(2)-24*CC(3)+14*CC(4)-3*CC(5))*(16.0I3.0) A5 - (CC(l)-4*CC(2)+6*CC(3)-4*CC(4)+CC(5))*(32.0I3.0) SIG) = AMf.0 + A3/3.0 + A414.0 + A5/5.0 SIG2 =-A2 SIG3 = -A3 SIG4 = -A4 IGS = -A5 C CALCULATE STRESS AT CRACK TIPS DO 170 CRACK - 1. ICMAX X =Z(CRACK)/TH STRESS(CRACK,TIME) - SlGI + S102*X + SIG3*(X**2) I + SIG4*(X**3) + SIGS*(X**4) C CALCULATE INFLUENCE FUNCTIONS DO 160M= 1,5 I(M = ZZ(M1) +X*ZZ(M,2)+ (X**2)*ZZQM.43)+ (X**3)*ZZ(M.4) IC(M) =ZZC(Mj)+X*ZZC(M2)4(X**2)*ZZC(MK3)+(X**3)*ZZC(MA) 160 CONTINUE A = Z(CRACK) C EQUATION FOR THE THERMAL STRESS INTENSITY TK(CRACKTIE) - ((3.1416*A)**.5)*(SIGI 1(1) I +51G2*I(2)*X+SIG3*I(3)*X**2 2 +SIG4*IK4)*X**3+SI050I(5)*X**4) TKC(CRACK,TilvE) =(3. 141 6A)**.5)*(SIO1 *IC(1)+SIG2*IC(2) I *X+51G3*IC(3)*X**2+5104*IC(4)*X**3+5105*IC(5)*X**4) B-7 . I I I 170 CONTINUE RETURN END SUBROUTINE KICLAD C THIS SUBROUTINE CALCULATES STRESSES AND STRESS INTENSITY FACTORS C DUE TO THE PRESENCE OF "CLADDING" ON THE I.D. SURFACE OF THE VESSEL INTEGER CRACK, TIME INTEGER CONSTE, CONSTK, Q REAL 10, 11 DO 170 TIME- I, 10 C CALCULATE STRESS DISTRIBUTION THROUGH VESSEL WALL C TEMP AT CLAD/BASE METAL INTERFACE TI = 0.5*(T(2,TIME) + T(3,TIME)) C TEMPERATURE AT THE VESSEL I.D. TO - T(,TIME) C STRESS-FREE TEMPERATURE TI= SFREET C CALCULATE STRESS DISTRIBUTION DUE TO CLAD K C SIGCI = STRESS IN CLAD AT VESSEL I.D. C SIGC2 = STRESS IN CLAD AT CLAD/BASE METAL INTERFACE C SIGBI = STRESS IN BASE METAL AT CLAD/BASE METAL INTERFACE C SIGB2 - STRESS IN BASE METAL AT VESSEL O.D. DELEA - CLADE*CALPHA*(I-ARATIO)/(I-CLADNU) C CALCULATE STRESS IN CLAD (KSI) SIGCI - DELEA* (TI - TO) SIGC2 = DELEA ( -TI) C CALCULATE FORCE DEVELOPED IN CLAD FCLAD = CTH*0.5*(SIGCI + SIGC2) C CALCULATE STRESSES IN BASE METAL (KSI) RO -RAD RI = RAD +CTH R2 =RAD+TH CONST = 1.0/((R2/RI) 2.0-1.0)*(RO-RI)/RI *DELEA I *(TI-O.5*(TO+TI)) SIGBI = CONST * (I + (R2/RI)**2.0) B-8 SIGB2 - CONST * 2.0 C CALCULATE FORCE DEVELOPED IN BASE M[ETAL FBASE -