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:
(CRACKTIME) = TWATER(TIME) - SUM(1) - SUM(2) - SUM(3) - SUM(4) C CONTROL FOR THE CONSTANT KAPPA OPTION IF (CONSTK.EQ. 1) GO TO 130 C TEST FOR KAPPA ACCURACY AND CONTROL OF KAPPA OPTION KTEST = 1.030 - (5.97E-7)*((T(CRACKTIvE))**2) IF ((ABS(KTEST-K)) IE. 0.0001) GO TO 130 K =KTEST GO TO 100 130 CONTINUE RETURN END B-5 ,I i , SUBROUTINE SKIT C Calculate Stress and Temperature at Crack-Tip and Thermal Stress C Intensity Factor, SKIt REAL E(5,10), CC(5), I(5), IC(5) INTEGER CRACK, TIME INTEGER Q, CONSTE, CONSTK C DETERMIINE POLYNOMIAL REPRESENTATION OF TEMPERATURE PROFILE C CONVERT CLAD TIHERMAL CONDUCTIVITY TO INCH AND MINUTE UNITS CCOND = CCOND / (12.0*60.0) COND = COND /(12.0*60.0) DO 105 TIME - 1, 10 TQI = TQ(I.TIME) TQ2 - TQ(2,TIME) TQ3 - TQ(3,TIME) TQ4 = TQ(4.TIME) TQ5 = TQ(5,TIME) Cl =TQI C2 = (-25*TQI+48*TQ2-36*TQ3+16*TQ4-3*TQS)/(3*TH) C3 = (35*TQI-104*TQ2+114*TQ3-56*TQ4+1 I *TQS)*(2.0/3.0*TH**(-2)) C4 = (-5*TQI+1 S*TQ2-24*TQ3+14*TQ4-3*TQS)*(16.03.0*TH**(-3)) C5 = (TQI-4*TQ2+6*TQ3-4*TQ4+TQ5)*(32.0/3.0*TH**(-4)) C CALCUATE TEMPRATURE AT THE CRACK TIPS DO 100 CRACK = 1, ICMAX T(CRACKTIME) = C1+C2*Z(CRACK)+C3*(Z(CRACK)**2) I -C4*(Z(CRACK)**3)+C5*(Z(CRACK)**4) 100 CONTINUE IF (CTH .LE. 0.0) GO TO 105 T(,TIME) = T(2,TIME) - (COND/CCOND)*(r(2,TIME)-T(ITME)) 105 CONTINUE IF (CONSTE .EQ. 1) GO TO 120 DO 10 TIME = 1, 10 DO 110 CRACK = 1,5 E(CRACKTIME) = 0.286+(5.400E-5 * (TQ(CRACK,TIMED))) 1 -(2.600E-8 * (TQ(CRACK,TIME))**2) 110 CONTINUE B-6 K K GO TO 140 120 DO 130 TIME = 1.10 DO 130 CRACK - 1, 5 E(CRACK,TIME) - EDATA 130 CONTINUE C DETERMIN POLYNOMIAL REPRESENTATION OF STRESS DIST 140 DO 170 TIME - 1, 10 DO 150 CRACK - 1, 5 CC(CRACK) - E(CRACKTIMAE)*TQ(CRACK.TIIE 150 CONTINUE Al =CCWI A2 - (-25*CC(1)48*CC(2)-36*CC(3)+16*CC(4)-3*CC(5)Yt3.0 A3- (35*CC(1)d104*CCC)+1 14*CC(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