ML16049A081
ML16049A081 | |
Person / Time | |
---|---|
Site: | Monticello |
Issue date: | 01/27/2016 |
From: | Haroon R, Lomas G, Pieper J AREVA |
To: | Division of Spent Fuel Management, Office of Nuclear Reactor Regulation |
Shared Package | |
ML16049A083 | List: |
References | |
L-MT-16-003, TAC L25058 11042-0205, Rev. 2 | |
Download: ML16049A081 (88) | |
Text
L-MT-1 6-0 03Enclosure 2ENCLOSURE 2AREVA, INC.CALCULATION 11042-0205, Revision 2TITLE:61 BTH ITCP and OTCP Closure Weld Flaw Evaluation 87 pages follow Form 3.2-1Calculation Cover Sheet ...Revision No.. 2A RE VA TIP 3.2 (Revision
- 10) Pageli of 87OCR NO (if applicable):
1104-023Revii~n PROECT 61 BTH Type 1 DS~s for114-2 eiin0POETNM:Monticello Nuclear Generating PlantPROJECT NO: 110942 CLIENT: Xcel EnergyCALCULATION
'tiTLE:61IBTH ITCP and OTCP Closure Weld Flaw Evaluation SUMMARY DESCRIPTION:
- 1) Calculation SummaryThis calculation qualifies Monticello DSC-16, a 61BTIH Type 1 DSCI for all design basis loads inconsideration of observed flaws in the Inner Top Cover Plate (ITCP) and Outer Top Cover Plate (OTOP)closure welds.2) Storage Media Location-Coldstor
-lareva_tn/i 1042/11042-0205.000
-Coldstor
-/areva~tn/i11042/110942-0205-002 If original issue, ,is licensing review per TIP 3.5 required?
NIAYes E] No [] (explain below) Licensing Review No.:SoftWare utilized (subject to test requirements of TIP 3.3): Software Version:
Software LogANSYS14.0Revision R-3 1Calculation is complete
..'Date:f " " ' "1120/2016 Originator Name .and Signature:
Jeff Pieper ________Cralcul'ation has been checked for consistency, Completeness, and correctness Date:Checker Name and Signature:
Gabriel LomasCalculation is approved for use HAROON Raheel 6.01 .27 Date:Project Engineer Name and Signature:
13:40:07
-05'00' 1/2712016 A Calculation No. 11042-0205 Aevi_ _ _ __.AR EVA Calculation Page 2of87REVISION SUMMARYAffected AffectedRev. Description Pages Data0 Initial Issue All All1-10, 13,Revised per DCR 11042-022 Revision
- 0. Made 14, 17, 18, Removed1 editorial clarifications, updated information from 21, 22, 24- etaeurevised Reference calculations, removed 34, 36, 45, data.extraneous sensitivity analyses.
46, 48-50,59-771-4, 6-10,13-17, 19-Revied er CR 1042023Revsio 0.22, 24, 27-Revied perlDsRi11042-023 Reviys sion 0pedi. 29, 31, 33 Added files in2 Addedeatic-pastdiscuso anals caictions appendixA 35, 44, 60, Appendix AAddiionl dscusionandclaifictios aded 62-63, 68-throughout.
78,Appendix A(79-87)
A Calculation No. 11042-0205 SRevision No. _ __AR EVA Calculation Page 3of87TABLE OF CONTENTSPage1.0 PURPOSE.......................................................................................................
72.0 ASSUMPTIONS.................................................................................................
83.0 DESIGN INPUT/DATA
.........................................................................................
93.1 DSC Geometry
...........................................................................................
93.2 Flaw Details and Geometry............................................................................
103.2.1 Outer Top Cover Plate.........................................................................
103.2.2 Inner Top Cover Plate ............................................................
.............
113.3 Material Properties......................................................................................
133.4 Design Criteria
..........................................................................................
134.0 METHODOLOGY..............................................................................................
144.1 Analysis Method and Acceptance Criteria............................................................
144.2 Load Cases..............................................................................................
174.3 FEA Model Details......................................................................................
224.3.1 Axisymmetric Case #1 .........................................................................
244.3.2 Axisymmetric Case #2.........................................................................
244.3.3 Axisymmetric Case #0.........................................................................
244.3.4 Half-Symmetry (3D) Case #1 .................................................................
254.3.5 Half-Symmetry (3D) Case #0 .................................................................
274.4 Limit Load Solution Details.............................................................................
2
75.0 REFERENCES
................................................................................................
286.0 ANALYSIS
.....................................................................................................
296.1 Axisymmetric Analyses for Internal Pressure
........................................................
296.1.1 Axisymmetric Case #1 -Initial Mesh Model.................................................
296.1.2 Axisymmetric Case #1 -Refined Mesh Models.............................................
296.1.3 Axisymmetric Case #2.........................................................................
316.1.4 Axisymmetric Case #0.........................................................................
316.2 Half Symmetry Analyses for Internal Pressure (Benchmark Cases)...............................
326.3 Half Symmetry Analyses for Side Drop Loading .....................................................
336.3.1 Half-Symmetry Case #1.......................................................................
336.3.2 Half-Symmetry Case #0.......................................................................
346.4 Evaluation of the 25g Corner Drop....................................................................
347.0 DISCUSSION AND CONCLUSIONS
........................................................................
358.0 LISTING OF COMPUTER FILES............................................................................
369.0 TABLES AND FIGURES......................................................................................
3810.0 APPENDIX A -ADDITIONAL RESULTS PLOTS AND ELASTIC-PLASTIC ANALYSES..............
79 A Calculation No. 11042-0205 Revision No. 2AR EVA Calculation Page 4of87LIST OF TABLESPageTable 1 -Summary of Design Basis Load Combinations for the 61 BTH DSC [Ref. 5.8]......................
38Table 2 -Internal Pressure in the 61 BTH Type I OSO ..........................................................
41Table 3 -Maximum Temperatures in the 61 BTH Type I DSC Shell............................................
41Table 4 -Properties of SA-240 Type 304. [Ref. 5.11] ............................................................
42Table 5 -Properties of SA-36. [Ref. 5.11] .........................................................................
43Table 6 -Summary of Load Cases and Results ...................................................................
44Table A-I -Summary of Elastic Plastic Analysis Results ........................................................
82 A Calculation No. 11042-0205 A Revision No. 2AR EVA Calculation Page 5of87LIST OF FIGURESPageFigure 1 -Sketch of the 618BTH DSC Top End and Transfer Cask from Reference 5.1 ......................
45Figure 2 -Details of the 618TH Top End Component Interfaces
.... ..........................................
46Figure 3 -ITCP and OTCP Closure Weld Details from Reference 5.5..........................................
47Figure 4 -DSC Top End Detailed Dimensions....................................................................
48Figure 5 -OTCP Flaws -Raw Data from Reference 5.1.........................................................
49Figure 6 -OTCP Flaws -Main Flaw Group Reduced and Bounded............................................
49Figure 7 -OTCP Flaws -Bounding Set #1 for ANSYS Collapse Analysis.....................................
50Figure 8 -OTCP Flaws -Bounding Set #2 for ANSYS Collapse Analysis
.....................................
50Figure 9 -ITCP Flaws -Raw Data from Reference 5.1..........................................................
51Figure 10 -ITCP Flaws -Bounding Flaw Set for ANSYS Collapse Analysis
..................................
51Figure 11 -Overview of the Axisymmetric Model .................................................................
52Figure 12- Mesh Details Near the Lid Regions of the Axisymmetric Model ...................................
52Figure 13 -Mesh Details at the Welds for Axisymmetric Case #1...............................................
53Figure 14 -Flaw Locations for Axisymmetric Case #1 ...........................................................
53Figure 15 -Refined Mesh (Weld Region) for Axisymmetric Case #1 ...........................................
54Figure 16 -Refined Mesh (Weld and Lid Interior Region) for Axisymmetric Case #1 .........................
54Figure 17- Mesh Details at the Welds for Axisymmetric Case #2...............................................
55Figure 18 -Flaw Locations for Axisymmetric Case #2 ...........................................................
55Figure 19 -Overview of the Half-Symmetry Model ...............................................................
56Figure 20 -Detail Views and Mesh Plots of the Half Symmetry Model.........................................
57Figure 21 -Isometric Views of Half-Symmetry Model ............................................................
58Figure 22 -Isometric Views of Half-Symmetry Model (Refined Circumferential Mesh) .......................
59Figure 23 -Results for Axisymmetric Case #1 -Initial Mesh -Service Level A/B ............................
60Figure 24 -Deflection History of the Center of the OTCP for the Axisymmetric Case #1 Initial Mesh......61 Figure 25 -Results for Axisymmetric Case #1 -Refined Mesh in Weld Region -Service Level A/B......62 Figure 26 -Results for Axisymmetric Case #1 -Refined Mesh in Weld and Lid Interior Region -ServiceLevel A/B.............................................................................................
63Figure 27 -Deflection History of the Center of the OTCP for the Axisymmetric Case #1 Refined Mesh....64 Figure 28 -Comparison of Maximum Displacement Histories for Axisymmetric Model Sensitivity Studies.
65Figure 29 -Comparison of Maximum Displacement Histories for Axisymmetric Model with Lid ContactDefined using Nodal DOF Couples vs. Contact Elements.......................................
66Figure 30 -Comparison of Maximum Displacement Histories for Axisymmetric Model With and WithoutPressure Loading Applied to the ITCP Weld Root Flaw Faces .................................
67Figure 31 -Results for Axisymmetric Case #2 -Refined Mesh in Weld and Lid Interior Region -ServiceLevel A/B.............................................................................................
68Figure 32 -Results for Axisymmetric Case #0 -Refined Mesh in Weld and Lid Interior Region -ServiceLevel A/B.............................................................................................
69Figure 33 -Comparison of Maximum Center-of-Lid Displacement Histories for the Various Flaw Models ..70Figure 34 -Results for Half-Symmetry Case #1 Internal Pressure Loading Benchmark Analysis
-ServiceLevelA/B
.............................................................................................
71Figure 35 -Benchmark of the Half Symmetry model with the Axisymmetric Analysis.........................
72Figure 36 -Equivalent Stress and Plastic Strain Plots from the Half-Symmetry
- 1 Side Drop Analysis
....73Figure 37 -Additional Results Plots from the Half-Symmetry
- 1 Side Drop Analysis
.........................
74Figure 38 -Equivalent Stress and Plastic Strain Plots from the Half-Symmetry
- 1 Side Drop Analysis withOff-Normal Internal Pressure.......................................................................
75Figure 39 -Equivalent Stress and Plastic Strain Plots from the Half-Symmetry
- 1 Side Drop Analysis withRefined Circumferential Mesh .....................................................................
76 A Calculation No. 11042-0205 Revision No. 2AR EVA Calculation Page 6of87Figure 40 -Equivalent Stress and Plastic Strain Plots from the Half-Symmetry
- 0 (No Flaws) Side DropAnalysis
..............................................................................................
77Figure 41 -Comparison of Maximum Displacement Histories for the Various Half-Symmetry Analyses....78 Figure A-I -Ramberg-Osgood Derived Stress Strain Curve for SA-240 Type 304 at 500 °F................
83Figure A-2 -Ramberg-Osgood Stress Strain Curves for SA-240 Type 304 from ANSYS Model at VariousTemperatures
.......................................................................................
84Figure A-3 -Service Level A Internal Pressure
-Equivalent Plastic Strain at 32 psi..........................
85Figure A-4 -Service Level D Internal Pressure
-Equivalent Plastic Strain at 65 psi..........................
86Figure A-5 -Service Level D Side Drop -Equivalent Plastic Strain at 75g.....................................
87 A Calculation No. 11042-0205
.Revision No. 2AR EVA Calculation Page 7 of 871.0 PURPOSEThe purpose of this calculation is to evaluate DSC-16 at the Monticello Nuclear Generating Plant (MNGP)per ASME Section III criteria in consideration of flaws observed in the Inner and Outer Top Cover Plate(ITCP and OTCP) closure welds. The flaws are documented in the Reference 5.1 Phased Array Ultrasonic Testing (PAUT) inspection report. The canister is a 61BTH Type 1 design. The ASME Section IIISubsection NB Code limits on primary stress are evaluated using the limit load analysis criteria prescribed inthe Code [Ref. 5.7].
A Calculation No. 11042-0205
.... ,Revision No. 2AR EVA Calculation Page 8of572.0 ASSUMPTIONS
- 1. The ITCP weld to the siphon and vent block and the welds of the siphon and vent port cover platesare inaccessible for PAUT inspection.
Approximately 11" are obscured due to the location of thesiphon and vent block. Whereas the main circumferential lid-to-shell welds are made with anautomated welding machine, some manual welding was performed around the siphon/vent block andports. As discussed in Section 3.4, a strength reduction factor of 0.8 is considered for both the ITCPand OTCP welds. This factor accounts for the siphon and vent block welds and uncertainties in theUT technique.
Note that the bounding flaws evaluated in this analysis are treated as full circumferential flaws. Inother words, it is not assumed that the siphon and vent block is free of flaws, but rather contains thesame bounding flaws as the examined welds. The geometry of the siphon and vent block is notassumed in this analysis.
It is assumed that the stresses in the circular configuration bound thestresses that would be computed for a configuration that explicitly includes the siphon and vent block.2. The longitudinal seams in the canister shell caused attenuation in the PAUT energy beam atlocations 24.3" to 24.8" and 129.5" to 130" [Ref. 5.11 that can potentially diminish the effectiveness ofthe examination in these half inch areas. These regions are considered limited examination zones. Itis assumed that the flaws observed outside of these regions are representative, and that no larger ormore bounding flaws exist in the regions behind the canister seam welds. The use of the 0.8 weldstrength reduction factor discussed above in Assumption 1 accounts for any uncertainty in this region.3. [Not used]4. The flaws are considered to be planar cracks lying on circumferential planes, parallel with thelongitudinal axis of the cask. (I.e. the crack tips are pointed in the axial directions of the cask). Thisis a conservative flaw orientation since the welds primarily resist normal stresses in the plane of thelids due to plate bending caused by DSC internal pressure.
Also, during the side drop loading,normal stresses in the plane of the lids resist the ovalizing mode of shell deformation.
This flaw orientation is also conservative for through-thickness shear stresses in the lid welds since itmaximizes the reduction in available shear area. (A flaw of equal length, but placed at an angle,would result in less reduction of the weld throat thickness).
- 5. Many of the flaws identified in the Reference 5.1 PAUT examination report lie in very similar locations within the weld cross section.
As discussed in detail in Section 3.2, flaws that lie in similar radial andaxial positions within the weld are considered bounded by a representative "group flaw." Thelocations and sizes of the "group flaws" are chosen conservatively to ensure they are bounding of theindividual flaws.6. The analysis is based on the nominal dimensions of the components as shown in the designdrawings
[References 5.3 and 5.4] including the as-fabricated radial gap between the outer diameterof the lids and the inner diameter of the DSC shell. Although weld shrinkage will close this gap duringclosure operations, the resulting compressive load path between the lids and shell is conservatively ignored.
Further discussion is provided in Section 4.3.
A Calculation No. 11042-0205
.... .Revision No. 2A.R EVA Calculation Page 9of873.0 DESIGN INPUT/DATA 3.1 DSC GeometryThe 61BTH Type 1 DSC geometry is detailed in the Reference 5.3 and 5.4 drawings.
The Reference 5.5drawing shows the details for the final ITCP and OTOP closure field welds. Sketches from Reference 5.1and details from References 5.3 and 5.4 are shown in Figure 1 through Figure 4.The material for all structural components (DSC Shell, OTCP, and ITCP) is SA-240 Type 304 stainless steel.The shield plug material is SA-36 carbon steel.The DSC shell is 0.5" nominal thickness.
The ITCP is 0.75" nominal thickness.
Per the Reference 5.5 drawing, it is welded to the DSC shell andvent/siphon block with a 3/16" groove weld. However, the ITCP lid groove (weld prep) is 0.25" minimum, andit was confirmed that the weld is also 0.25" [Ref. 5.1].The OTCP is 1.25" nominal thickness.
It is welded to the 050 shell with a 1/2" groove weld.The ITOP and OTCP closure welds (with the exception of the ITCP welds around the vent/siphon block andthe welds of the vent and siphon port cover plates) are made using the GTAW process with an automated welder. This is a non-flux type of weld. The vent/siphon block and the vent and siphon port cover platewelds are performed
- manually, also using a non-flux process.
A Calculation No. 11042-0205 Revision No. 2AR EVA Calculation Page 10of873.2 Flaw Details and GeometryVarious sets of bounding flaws are chosen for the detailed analyses based on the flaw dimensions inReference 5.1 and the discussion below. Note that flaws are identified in this calculation using the numerical flaw listings in the Reference 5.1 inspection report.3.2.1 Outer Top Cover Plate3.2.1.1 Case 1Figure 5 shows all of the OTCP weld flaws from Reference 5.1 plotted on an outline of the DSC geometry.
Figure 6 shows a similar plot but with the main cluster of flaws bounded by a box, and showing arepresentative "group flaw" for this region. The longest flaw within the group region is 31 .7" long and thetallest flaw is 0.14" high. Therefore, the bounding flaw for this region is taken as a full circumferential flaw,0.14" in height.Note that all flaws in the group region were reviewed to ensure that no two flaws in close circumferential proximity, considered as being joined, would produce a taller flaw. For example, OTCP Flaw #9 and OTOPFlaw #10 are within 0.17" of each other in the circumferential direction, but their combined height is only0.47-0.38=0.09".
Therefore these flaws, considered
- combined, are bounded by the 0.14" high group flaw.The radial and axial positions of the bounding flaw were chosen to be at the center of the group region. Thisradial position is within the critical failure plane of the weld (i.e. a plane containing the minimum weld throatthickness of 0.5").Figure 6 also shows additional information about the flaws outside of the group region. OTCP Flaw #2 isintermittent around the entire circumference of the DSC. Therefore this flaw, at 0.12" in height, is considered a full circumferential flaw. Since OTCP Flaw #14 is in close proximity to Flaw #2, it is conservatively considered joined to OTOP Flaw #2, and the combined flaw height is considered to be present around theentire circumference.
The combined flaw height is determined based on the geometry to be 0.195".As seen in Figure 6, OTCP Flaw #20 is remote from the group region and from OTOP Flaw #s 2 and 14.OTCP Flaw #20 is only 0.32" in length, and only 0.07" in height. This flaw is separated from OTCP Flaw #19by 0.36" in the circumferential direction and by 0.19" in the axial direction.
It is separated from OTOP Flaw#21 by 1.66" in the circumferential direction and by 0.23" in the axial direction.
Since extension of the flawsunder the postulated loading is negligible (since only one cycle of the critical loads is applied) this flaw willnot join with the adjacent flaws. Additionally, since OTOP Flaw #20 is much smaller than the critical surfaceflaw size of 0.29" from Reference 5.17, it is not considered explicitly in the FEA analyses and is considered bounded by the other modeled flaws which are very conservative.
Similarly, OTCP Flaw #3 is remote from all flaws with the exception of OTCP Flaw #2. However, OTCP Flaw#3 is very small, only 0.18" long and 0.09" tall. Inspection of the PAUT plots (see Page 22 of Reference 5.1)also shows that OTCP Flaw #2, which is considered as fully continuous in this analysis, is actually veryintermittent at the circumferential position of OTCP Flaw #3. Furthermore, OTCP Flaw #3 is much smallerthan the critical subsurface flaw size of 0.29" from Reference 5.17. Therefore, it is not considered explicitly in the FEA analyses and is considered bounded by the other modeled flaws which are very conservative.
Figure 7 shows the first bounding flaw set considered for the OTCP in the ANSYS collapse analyses.
A Calculation No. 11042-0205 S ,Revision No. _ __AIREVA Calculation Page li of 873.2.1.2 Case 2The discussion above and the flaw locations shown in Figure 5 through Figure 7 are based primarily on thetabulated flaw data from Reference 5.1. Since OTCP Flaw #2 is intermittent around the circumference of theweld, a closer inspection of the PAUT scan images is performed, and an additional flaw set for the OTCP iscreated.
In this additional case, the location of OTCP Flaw #2 is based on the PAUT scan image of the flawat the circumferential position of OTCP Flaw #14, which is the only additional flaw that could be considered to interact with OTCP Flaw #2. Based on the PAUT scan images, the flaws are located as seen in Figure 8.In this case the height of both Flaw #2 and Flaw #7 are estimated based on the PAUT scan images and areconservatively larger than the flaw heights tabulated in Reference 5.1.3.2.2 Inner Top Cover PlateFigure 9 shows all of the ITCP weld flaws from Reference 5.1 plotted on an outline of the DSC geometry.
Allbut two of the flaws are clustered in the region of the weld root at the inner surface of the DSC shell. Figure10 shows the bounding flaw set considered for the ITCP in the ANSYS collapse analyses.
Both therepresentative group flaw and ITCP Flaw #7 are considered to be full circumferential flaws. ITCP Flaw #11is remote from all other flaws (in the circumferential direction) and is therefore considered bounded by therepresentative group flaw. The representative group flaw for the ITCP is conservatively placed at the tensionside of the weld when resisting internal pressure.
All of the ITCP flaws documented in Reference 5.1 were reviewed to ensure that no two (or more) flaws,which are in close proximity to each other, could be considered as combined and therefore creating a morecritical flaw. The following cases are considered in particular:
- ITCP Flaw #2 and Flaw #3 are within 0.12" from each other in the circumferential
- position, but theirmaximum combined height (1.58-1.49
= 0.09") is bounded by the group flaw height of 0.09".*ITCP Flaw #5 and Flaw #8 partially overlaps with Flaw #6 in the circumferential direction and wouldhave a combined height of 0.15". However, Flaw #5 (0.15" in length) and Flaw #8 (0.14" in length)are extremely small. Due to their overlap in the circumferential direction, their combined length wouldbe only 0.16", and therefore would not affect the global or local stability of the weld. This very shortregion with a potential 0.15,, high flaw is bounded by the full-circumferential representation of themodeled flaws.*ITCP Flaw #10 is within 0.04" of Flaw #12 in the circumferential direction.
The individual flaws are0.05" tall and 0.04" tall, respectively, and 0.49" long and 0.18" long, respectively.
They are alsoseparated in the axial direction by 0.09". Postulating a flaw from the bottom of Flaw #12 to the top ofFlaw #10 would imply a height of 0.18". However, the combined-height region would be over a veryshort length and would not affect the global or local stability of the weld. Therefore this postulated combined flaw is considered bounded by the full-circumferential representation of the modeled flaws.It is noted that based on Figure 9 and Figure 10, ITCP Flaw #7 appears to be in the base metal of the innertop cover plate. It is likely that the flaw is actually at the fusion / heat affected region between the weld metal A Calculation No. 11042-0205
- , Revision No. 2A REVA Calculation Page 12of87and the base metal. The ANSYS models used in this calculation place the flaw at 0.81" inward from theouter surface of the DSC shell whereas the tabulated data in Reference 5.1 places the flaw at 0.80" from theouter surface.
The 0.01" discrepancy is considered negligible.
The exact location of the flaw is notconsidered critical in light of the significant margin that is available (See Section 7.0) and the generally veryconservative idealization of the flaws (i~e. full circumferential).
A Calculation No. 11042-0205
,,, Revision No. 2A REVA Calculation Page 13of873.3 Material Properties The material properties for the DSC structure are taken from Reference 5.11. The properties of the twomaterials of construction, SA-240 Type 304 stainless steel and SA-36 carbon steel, are provided in Table 4and Table 5, respectively.
The weld metal is considered to be composed of the same properties as the basemetal, as the welds are made with the non-flux GTAW method [Reference 5.14] using bare metal ER308(stainless) filler material.
The tensile strength of the ER308 electrode (80 ksi at room temperature
[Ref.5.15]) is slightly greater than the type 304 base metal (75 ksi at room temperature
[Ref. 5.16]). The yieldstress value of the weld metal is assumed to be equal to or greater than the base metal. Therefore, thetreatment of the weld metal as being identical to the base metal is appropriate for the Section III limit loadanalyses performed in this calculation.
Temperatures used for material properties are discussed in Section 4.2 and are shown in Table 3.Poisson's ratio for all modeled parts is taken as 0.29.Weight density for SA-240 Type 304 is taken as 0.285 lb/in3.Weight density for SA-36 is taken as 0.284 lb/in3.3.4 Design CriteriaAll of the applicable design basis loading conditions are considered in accordance with the requirements ofASME Section III Subsection NB [Ref. 5.7]. Section 4.1 details the methods used to perform the code [Ref.5.7] qualifications.
Section 4.2 details the selection of the bounding load cases.The mockup used in the PAUT process development contained weld manufacturing flaws intentionally distributed in locations that would be expected with the weld process used for the OSO lid closure welds.Approximately 30% of those flaws were placed at the weld root and 27% were placed near the weld toe todemonstrate that they could be reliably detected in the presence of typical geometric responses from thoseregions.
The flaws include incomplete root penetration, lack of fusion, and tungsten inclusions.
AREVAdocument 54-PQ-114-001
[Ref. 5.19], Section 8.0, provides images of the UT responses for these flaws anddemonstrates that the PAUT process can effectively detect these flaws. Furthermore, the qualification performed on the blind mockup provides objective evidence that detection of flaws in these regions of theweld is not a problem.
The blind mockup used for qualification contained a similar percentage/number anddistribution of flaws as the development mockup. Although the flaw information for the blind mockup cannotbe disclosed in order to preserve the security of the mockup for future qualifications, EPRI and NRC3personnel present at the demonstration have reviewed that information.
In addition, uncertainties in thePAUT examination are accounted for by using a 0.8 reduction factor on the limit load. This factor, which is inagreement with ISG-1 5 [Ref. 5.20], conservatively accounts for any additional limitations in the efficacy of thePAUT examinations and also accounts for the inaccessible area around the vent and siphon block as well asthe geometric reflectors at the root and near the toe of the weld.
A Calculation No. 11042-0205 Revision No. 2AR EVA Calculation Page 14 of 874.0 METHODOLOGY 4.1 Analysis Method and Acceptance CriteriaThe 61BTH DSC including the ITCP and OTCP welds are designed and analyzed per ASME Section IIISubsection NB (the Code) [Ref. 5.7] in the Reference 5.2 calculation.
The presence of the ITCP and OTCPweld flaws will cause high local stresses and complex stress fields that will render an elastic analysis (suchas those performed in Reference 5.2) very difficult.
Therefore, the flaws are explicitly included in the finiteelement models as "design features",
and the applicable ASME code [Ref. 5.7] stress limits are evaluated asdescribed below.Primary Stress LimitsIn order to satisfy the primary stress limits of Reference 5.7 paragraphs NB-3221 .1, NB-3221 .2, and NB-3221.3, a Limit Analysis will be performed per Paragraph NB-3228.1.
The acceptance criterion is that thespecified loadings not exceed two-thirds of the lower bound collapse load, as determined using an ideallyplastic (non-strain hardening) material model, with the yield stress set at a value of 1 .5*Smn. This criterion isused for evaluation of the Service Level A and B load cases discussed in Section 4.2.Note that Service Level C acceptance criteria are generally 20% greater than Service Level A criteria, perParagraph NB-3224 of Reference 5.7. This information is used in the discussion in Section 4.2 to eliminate some non-critical load cases.For the Service Level D loadings (accident level internal pressure and side drop), the rules of ASME SectionIll Appendix F Paragraph F-1 341 .3 [Ref. 5.9] are used, which indicate that the loads "shall not exceed 90%of the limit analysis collapse load using a yield stress which is the lesser of 2.3Sm and 0.7Su." This criterion is used for evaluation of the Service Level D load cases discussed in Section 4.2.An additional increase factor of 1/0.8=1.25 is applied to the required limit load collapse pressure in order toaccount for the weld strength reduction factor of 0.8 to account for UT sensitivity and inaccessible weldregions discussed in Section 3.4. Typically, the weld strength reduction factor is applied to the weldallowable stress during qualification.
In the case of limit-load
- analysis, reduction of the material yield stressis applicable.
The reduction in yield stress would have a direct, 1:1 correlation to the calculated lower boundcollapse pressure due to the perfectly-plastic (i.e. non-strain hardening) material model. In this analysis, rather than decrease the material yield stress the required calculated collapse pressure is increased by thefactor of 1.25.Note that the Service Level D criterion is essentially 2.1 times greater than the Service Level A/B criterion, ascalculated below. This information is used in the discussion in Section 4.2 to eliminate some non-critical load cases.At a temperature of 500 °F, the limit load yield stress for SA-240 Type 304 for Service LevelsA/B and D are 26.3 ksi and 40.3 ksi, respectively.
The code [Ref. 5.7] required factors against the lower bound collapse load as determined bythe limit load analyses for Service Levels A/B and D are 1.5 and 1.11, respectively.
(40.3X1.E)
The ratio of the acceptance criteria is therefore:
(2.,.)= 2.1.(i.e. the Service level A/B criteria are 2.1 times as severe)
~Calculation No. 11042-0205
,. .:,Revision No. 2AREVA Calculation Page 15of87Limit Load Analysis Backaqround ASME Section III Subsection NB provides only a basic description of the Limit Load analysis technique.
Amore thorough description is provided in ASME Section VIII Division 2 Paragraph 5.2.3 [Ref. 5.18]:Limit-load analysis addresses the failure modes of ductile rupture and the onsetof gross plastic deformation (plastic collapse)...
Limit-Load analysis provides an altemnative to elastic analysis and stresslinearization and the satisfaction of prim ary stress limits...
Displacements and strains indicated by a limit analysis have no physicalmeaning.Limit load analysis is based on the theory of limit analysis that defines a lowerbound to the limit load of a structure as the solution of a numerical model withthe following properties:
- 1. The material model is elastic-perfectly plastic with a specified yield strength.
- 2. The strain-displacement relations are those of small displacement theory.The limit load is the load that causes overall structural instability.
This point isindicated by the inability to achieve an equilibrium solution for a small increasein load (i.e. the solution will not converge).
Separately, in order to address questions on the potential for material rupture due to potentially high plasticstrains, supplemental elastic-plastic analyses are performed in Section 10.0 (Appendix A).Primary Plus Secondary Stress LimitsThe Code [Ref. 5.71 also prescribes limits on primary plus secondary stresses for Service Levels A and B[Ref. 5.7 Paragraph NB-3222.2].
Secondary stresses may be developed in the DSC due to differential thermal expansion of the interconnected parts and thermal gradients within the structure.
The code stresslimit for primary plus secondary stress (calculated on an elastic basis) is 3Sin. However, as shown in Ref. 5.7Figure NB-3222-1, rules for exceeding the 3Sin limit are provided in Paragraph NB-3228.5, which states that"the 3Sm limit .., may be exceeded provided that the requirements of (a) through (f) below are met."Requirement (a) states that "the range of primary plus secondary membrane plus bending stress intensity, excluding thermal bending stresses, shall be < 3Srn." This provision is related to the potential for "plasticstrain concentrations" occurring in "localized areas of the structure",
and the potential for theseconcentrations to affect the "fatigue
- behavior, ratcheting
- behavior, or buckling behavior of the structure"
[Ref.5.7 Paragraph NB-3228.1].
Requirements (b) through (d) are also limitations related to fatigue and thermalstress ratchet.
As detailed in Section 10.5 of Reference 5.2, the DSC is exempt from fatigue analysisrequirements since all of the criteria in NB-3222.4 of Reference 5.7 are satisfied.
Similarly, since the DSCthermal loads are not cyclic in nature (other than small daily and seasonal fluctuations),
thermal stressratchet is not a concern.
Therefore, the 3Sm limit as it relates to fatigue is not applicable.
A Calculation No. 11042-0205
~Revision No. 2AR EVA Calculation Page 16 of 87Requirement (e) requires that the component temperature be less than 800 °F for austenitic stainless steel.The maximum DSC shell temperature (entire shell including the lid region) is 611 °F (See Table 3).Therefore this requirement is satisfied.
Requirement (f) states that the material must have a specified yield stress to ultimate stress ratio of less than0.8. For the 61 BTH DSC which used SA-240 Type 304 steel, the ratio is 30/75 = 0.4. Therefore thisrequirement is satisfied.
Based on the discussion above (primarily the fact that cyclic conditions are not a design factor for the DSC),there is no need to consider limits on primary plus secondary stresses.
Therefore, thermal stresses are notincluded in this analysis.
Special Stress LimitsIn addition to the primary and primary plus secondary stress limits the Code [Ref. 5.7] also imposes SpecialStress Limits as detailed in paragraph NB-3227.
The applicable special stress limits are discussed below inrelation to the DSC top end cover plate welds.Bearing Loads: There are no significant bearing loads affecting the ITCP and OTCP closure welds duringService Level A, B, or C loading.
During the Service Level D side drop event, bearing stress exists at thecontact surface between the DSC and Transfer Cask. However, as noted in ASME Section III Appendix F[Ref. 5.9] paragraph F-1331.3, bearing stress need only be evaluated for pinned and/or bolted joints.Therefore this special stress limit is not applicable to this evaluation.
Pure Shear: Although the ITCP and OTCP closure welds are loaded in shear by internal pressure loading,the stress state is not pure shear due to the additional bending stresses.
Paragraph NB-3227.2 ofReference 5.7 clarifies that this stress limit is applicable to "for example, keys, shear rings, screw threads."
Therefore this special stress limit is not applicable to this evaluation.
Progressive Distortion of Nonintegral Connections:
The ITCP and OTCP closure welds are integral andtherefore not nonintegral connections.
Furthermore, there are no sources of significant cyclic loading thatwould cause progressive distortion of the DSC. Therefore this special stress limit is not applicable to thisevaluation.
Triaxial Stress: The purpose of the code [Ref. 5.7] limit on triaxial stress is to provide protection againstfailure due to uniform triaxial tension [Ref. 5.13 Chapter 4.5]. Internal pressure in the DSC and bending ofthe cover plates may cause tension in the weld in the radial and circumferential directions, but there is nosource for tension in the axial direction.
Therefore failure due to hydrostatic tension in the weld metal is notcredible.
Therefore this special stress limit is not applicable to this evaluation.
Fracture and Flaw Extension Although linear-type flaws have been identified in the structure, the critical failure mode of the welds is plasticcollapse.
Under one-time
- loading, elastic and plastic crack extension are not a concern for the very toughtype 304 stainless steel materials of the DSC shell, OTCP, and ITCP. This conclusion is supported byASME Section Xl Article C-4000 "Determination of Failure Model" [Ref. 5.10] which states that for austenitic wrought material and non-flux welds, "plastic collapse is the controlling failure mode." Note that the 61 BTHType I DSC OTCP and ITCP closure welds are made with the GTAW method [Reference 5.14] which is anon-flux type of weld.
A Calculation No. 11042-0205 A _ _ _ __io N.A R EVA Calculation Page 17of87Additionally, there is no source for fatigue flaw extension.
The only cyclic loads on the 050 are minor dailyand seasonal temperature fluctuations.
Therefore, cyclic fatigue growth of the flaws in not a crediblephenomenon.
Combined with the discussions above, the limit load analysis of the DSC top cover plates and closure weldsis sufficient to satisfy all of the applicable stress criteria of the Code [Ref. 5.7].Residual StressResidual stress due to welding is a secondary stress and therefore is not considered in the limit loadanalyses performed in this calculation, as the Section III Code [Ref. 5.7] does not require it in the limit loadanalysis.
4.2 Load CasesTable 1 lists the design basis load combinations for the 61 BTH DSC. This calculation is concerned with allload cases beginning with the inner top cover plate weld, identified as Load FL-6 in Table 1.The loading conditions of interest in this evaluation are internal and external pressure and inertial loads dueto handling,
- transfer, seismic, and accidental drop conditions.
As discussed in Section 4.1, secondary (thermal) loading is not considered.
Note that the discussions below, and the analyses performed in this calculation, are based on theconservative design values for internal pressure
- loading, rather than the actual calculated values of internalpressure.
Table 2 summarizes the conservative design values as well as the actual calculated values.Temperatures used for the material properties for each Service Level condition are listed in Table 3 anddiscussed further in the paragraphs below.Service Level AThe bounding Service Level A load combination for the DSC top end cover plates and welds is load caseTR-5 which combines the hot ambient condition with internal pressure and lg axial inertial loading.
Theother directions of inertial loading are not considered critical since their effects are not directly additive to theinternal pressure
- loading, and furthermore they are bounded by the 75g side drop load discussed furtherbelow.The lg axial load will cause the DSC payload weight (fuel, basket, holddown ring, shield plug) to bearagainst the ITCP. The total maximum payload weight is 75,811 lbs conservatively including the weights ofthe ITCP and OTCP [Ref. 5.2 Section 10.2]. The equivalent uniform pressure applied to the top-endcomponents is therefore:
Pfuellg = 75,811 =2. sP .et l =x (66.25 tn)2=2. sWhere the inner diameter of the DSC shell is 66.25 inches.
A Calculation No. 11042-0205 Revision No. 2A R EVA. Calculation Page 18of87Therefore, the bounding Service Level A case is a uniform 10 psi internal pressure (for a Type 1 DSC) plusan additional 22.0 psi acting on the shield plug in the outward axial direction of the 050 Shell.Conservatively, this analysis considers the combined 10+22=32 psi load as a uniform internal pressure in theDSC Shell. This is very conservative since the fuel pressure load which is applied to the inner surface of theshield plug would in reality be distributed to the perimeter of the ITOP as a line load by the significant stiffness of the 7-inch thick shield plug. In other words, the approach used in this calculation maximizes thebending loads on the cover plates and therefore maximizes the loading on the closure welds.Note that the cases with external pressure loading are discussed below.
A Calculation No. 11042-0205 Revision No. _ __--AR EVA Calculation Page 19of87Service Level BThe bounding Service Level B load combination for the DSC top end cover plates and welds is thecombination of the hot ambient condition with the off-normal internal pressure of 20 psi (LD-6). All of theother Service Level B conditions, such as ram push/pull loads, do not affect the top end components.
Therefore, the bounding Service Level B case is a uniform 20 psi internal pressure.
Since the pressureloading is smaller (20 psi for SL B versus 32 psi for SL A as described above), the temperature used for SL A(500 °F) bounds the maximum SL B temperatures (416 °F), and since the same limit load acceptance criterion is used for Service Levels A and B, this case is bounded by Service Level A.Service Level CThe bounding Service Level C load combination for the DSC top end cover plates and welds is one thatbounds HSM-4 and HSM-8 which combines the hot ambient condition, normal internal pressure (20 psi), andseismic loading.
- However, the seismic loads are bounded by the handling loads [Ref. 5.2 Section 7.8]discussed above for Service Level A. In addition, the acceptance criteria for Level C limit load analysis isgreater than Service Levels A and B. Therefore, all Service Level C conditions are bounded by the ServiceLevel A case described above.Note that the other Service Level C cases (such as LD-7 and UL-7) are for accident condition DSC rampush/pull loads. These loads do not affect the DSC top end components.
Therefore they are not applicable to this analysis.
Note that cases with external pressure loading are discussed below.
A Calculation No. 11042-0205 Revision No. 2AR EVA Calculation Page 20of87Service Level DThree load combinations are found to be critical for Service Level D loading of the DSC top end components, namely:* accident level internal pressure* corner drop* side dropThe first load combination is bounded by HSM-5 or HSM-6 which consist of 65 psi internal pressure due toHSM blocked vent thermal conditions.
This load is not combined with any other load that affects the top-endcomponents.
Therefore, the first bounding Service Level D load case considered in this analysis is 65 psiinternal pressure.
Note that in this condition the maximum DSC shell temperature is 611 °F and 625 0F isconservatively used in this analysis (See Table 3).The other Service Level 0 conditions consist of the drop events and accident-level seismic loading.
Theaccident seismic loads are bounded by the handling loads [Ref. 5.2 Section 7.8] discussed above for ServiceLevel A. The end-drop load is not a credible event [see footnote 12 to Table 1] but was used in the originalcalculation
[Ref. 5.2] to bound the corner drop event. However, that analysis produced negligible load in thetop cover plate welds due to the idealized boundary conditions.
As a result of an RAI by the NRC, the cornerdrop is considered using an alternate idealization that maximizes the load in the top cover plate welds. Inthis case, the 25-g corner drop load has an axial component that may be considered to load the top endcover plates with the inertia of the fuel, shield plug, hold-down ring, ITOP and OTCP. This case is evaluated in Section 6.4.The 75g side drop load TR-1 0 is considered a critical load case and is evaluated in detail. Note that this loadcase represents 75x more load than the Service Level A lg inertial loads. As discussed in Section 4.1, theService Level D acceptance criterion is only 2.1 times less stringent than the Service Level A/B criterion.
Therefore, evaluation of the 75g side drop case using the Service Level D criterion is bounding of theService Level A transverse inertial loading.
(Also, as discussed in Section 4.3, the boundary conditions usedfor the 75g side drop analysis are conservative and representative of the boundary conditions encountered for the Service Level A inertial loads and seismic loading.)
The 75g side drop case also includes the off-normal internal pressure of 20 psi, as shown in Table 1.Note that the side drop event TR-10 occurs during transfer operations which result in a maximum DSC shelltemperature of 500 °F as shown in Table 3. The higher Service Level D temperature of 625 0F discussed above occurs only during DSC storage in the HSM, and therefore is not combined with the side drop loading.
A Calculation No. 11042-0205 Revision No. _ __AR EVA Calculation Page 21lof 87External Pressure LoadingExternal pressure is present on the DSC in load cases DD-2 (vacuum drying, Service Level A) and HSM-9/10 (flood load, Service Level C). (The load cases with hydrostatic external pressure are due to the caskannulus being filled with water while the cask and DSC are in the vertical position.
In this case the pressureload varies from zero at the top of the DSC to a maximum value at the bottom of the cask. Since theexternal pressure near the cover plates is essentially zero, these cases are not critical and are notconsidered further in this calculation.)
In load case DD-2, the external pressure is 14.7 psi (full vacuum).
This pressure is bounded by the ServiceLevel B off-normal pressure (20 psi) and therefore primary stresses in the cover plates and welds arebounded by the internal pressure load cases. Stability concerns of the DSC shell are not affected by thepresence of weld flaws since they are at the end of the cask, remote from the locations at which bucklingwould occur. Additionally, the external pressure is resisted directly by the shield plug and the shield plugsupport ring, rather than by the OTCP and ITCP welds. Therefore external pressure load case DD-2 is notcritical and is not considered further in this analysis.
In load case HSM-9/1 0, the flood load is due to a 50-foot static head of water, which is equivalent to 22 psiexternal pressure
[Ref. 5.2 Section 7.9]. This pressure is bounded by the 32 psi internal pressureconsidered for Service Level A discussed above. Therefore the flood load case HSM-9/1 0 is bounded by theother internal pressure load cases.SummaryThe bounding load cases considered for the limit load collapse analyses are therefore:
(See Table 3 for temperature references)
(See Section 4.1 for explanation of the 1.5 and 1.11 factors for Service levels A/B and D, respectively, andalso for the 0.8 factor which accounts for limitations in the weld examination and inaccessible weld regions,as discussed in Section 2.0, Assumption No. 1.Service Level A/B: 32 psi Uniform Internal
- Pressure, Properties at 500 0F(Accounts for internal pressure
+ inertial load of DSC contents onto Lid)Limit load collapse pressure required to satisfy criteria:
1.5*32/0.8
= 60 psiService Level D-1: 65 psi Uniform Internal
- Pressure, Properties at 625 °FLimit load collapse pressure required to satisfy criteria:
1.11*65/0.8
= 90.2 psiService Level D-2: 75g Side Drop Acceleration plus 20 psi Uniform Internal
- Pressure, Properties at 500°F.Limit load collapse acceleration required to satisfy criteria:
1.11*75/0.8
= 104 g
~Calculation No. 11042-0205
- Revision No. 2AR EVA Calculation Page 22of874.3 FEA Model DetailsSeveral finite element models of the top half of the 61 BTH DSC are constructed in ANSYS based on theReference 5.3, 5.4, 5.5 drawings.
The models fall into two basic categories:
axisymmetric (2D) and half-sym metric (3D).The axisymmetric models use ANSYS plane element type PLANE182, a 4-node axisymmetric plane elementwith non-linear capabilities.
Each node has 2 degrees of freedom (translation in the X (radial) and Y (axial)directions).
The default element options are used in the analysis.
Sensitivity studies were performed toensure that there were no adverse effects on the results due to the potential shear locking of the elements.
(Sensitivity runs used KEYOPTION 1=3 to invoke the simplified enhanced strain formulation to relieve shearlocking.)
Additional discussion of the sensitivity analyses is provided in Section 6.0.Contact between the ITCP and OTCP is simulated using nodal coupling in the Y (axial) direction.
(SeeSection 6.1.2 for a sensitivity study using contact elements at this interface.)
No contact is defined between the opposing faces of the weld flaws. In other words, whereas compressive loading normal to the plane of the flaw may in reality be transmitted via compression through the crack facesurfaces, this load path is ignored.
This is conservative, and considered necessary since it is difficult (orimpossible) to deduce from the PAUT data what separation may exist between the two faces of the flaws.Also, no contact is considered between the DSC shell inner diameter and the ITCP and OTCP outerdiameters.
As seen in Figure 4, the fabricated dimensions of the lids and shell result in small radial gapsbetween the outer diameter of the lids and the inner surface of the shell. During the welding process, thesegaps close, but since a small remaining gap cannot be ruled out, this analysis conservatively assumes thatthe as-fabricated gap exists, as shown in Figure 4. Even if the lids deflect in the analysis such that the gapswould close, the resulting contact/compressive load path is conservatively neglected.
This is conservative since it forces all loads in the lid to travel through the weld, rather than through compression between the lidsand shell.Figure 11 and Figure 12 show images of the axisymmetric model. Loading and boundary conditions arediscussed in the following sections.
A Calculation No. 11042-0205
~Revision No. 2A REVA Calculation Page 23 of 87The 3D, half-symmetric model uses ANSYS solid element type SOLID1 85, an 8-node brick (or 6-node prism)element with non-linear capabilities.
Each node has 3 degrees of freedom (translation in the X, Y, and Zdirections).
The default element options are used in the analysis.
Sensitivity studies were performed toensure that the mesh was adequate.
Additional discussion of the sensitivity analyses is provided in Section6.0.Contact in the half-symmetry model is defined using ANSYS element types CONT173 and TARGE170.
Contact is defined between the following interfaces:
- OTCPto ITCP* ITCP to Shield Plug* Shield Plug outer diameter to DSC Shell* Shield Plug bottom surface to Support Ring* Support Ring to DSC ShellThe default contact parameters are used, although the contact stiffness is reduced in some cases to aid inconvergence.
Due to the large contact areas and since the contact areas are generally remote from thecritical stress regions, the contact stiffness is not considered a critical parameter.
The default contactparameters include:
[Reference 5.6]* Penetration tolerance factor: Default value = 0.1. This parameter controls the acceptable level ofpenetration of the contact node into the target surface, based on the depth of the element underlying the target element.* Pinball region scale factor: Default Value = 1.0. This parameter controls the extents of the regionaround each contact node that is checked for contact with target segments.
The default volume is asphere of radius 4*depth of the underlying element.* KeyOption 2: Contact algorithm:
Default = Augmented Lagrangian.
The contact method is aniterative penalty method where the contact pressure is augmented during the equilibrium iterations sothat the final penetration is within the acceptable tolerance.
- KeyOption 4: Location of contact detection point: Default = On Gauss Point. Other options includeusing the nodal points, normal to either the contact surface or the target surface.
The default optionis suggested for general cases.Other features and controls of the CONTA1 73 elements are related to advanced features (bonded contact,cohesion, etc.) and initial penetration and gap controls which are not utilized in this analysis.
Figure 19 through Figure 21 show images of the half-symmetry model. Loading and boundary conditions arediscussed in the following sections.
Table 6 shows a summary of the ANSYS models and analyses which are performed.
Further details on thevarious ANSYS models are provided below.
A Calculation No. 11042-0205
~Revision No. 2AR EVA Calculation Page 24of874.3.1 Axisymmetric Case #1The first case considered is a combination of OTCP Flaw Set #1 and the ITCP bounding flaw set discussed in Sections 3.2.1.1 and 3.2.2, respectively.
The mesh and flaw details for this case, called Axisymmetric Case #1, are shown in Figure 13 and Figure 14.The mesh shown in these figures was created based on a basic goal of having at least 4 elements acrossthe thickness of the net sections of the weld, as reduced by the flaws. In order to investigate the effects ofmesh density, a refined mesh was created for this case, as shown in Figure 15. Since the sensitivity modelshown in Figure 15 only refined the weld region an additional model was created as shown in Figure 16 toensure a sufficient mesh in the lid interior region.This model, and all of the other axisymmetric models discussed below, are used for analysis of uniforminternal pressure loading.
The model is constrained in the radial direction at the axis of symmetry and in theaxial direction at the bottom cut of the DSC shell near the mid-length of the cask (remote from the top endcomponents of interest.)
The pressure loading is applied to the internal pressure boundary (bottom surfaceof ITCP, surface of ITCP weld to Shell, and Shell inner surface).
(See Section 6.1.2 for a sensitivity analysiswhere internal pressure is included on the ITCP weld root flaw internal surfaces.)
4.3.2 Axisymmetric Case #2The second case considered is a combination of OTCP Flaw Set #2 and the ITCP bounding flaw setdiscussed in Sections 3.2.1.2 and 3.2.2, respectively.
The mesh and flaw details for this case, calledAxisymmetric Case #2, are shown in Figure 17 and Figure 18 for the refined mesh. Based on the results ofthe Axisymmetric Case #1 (See Section 6.1.2), the initial mesh level described above for Case #1 would besufficient.
- However, since the run times remained reasonable, a refined mesh model (weld and lid interiorregions) was generated and is used for Case #2.4.3.3 Axisymmetric Case #0In order to study the effect of the flaws, a 3rd case is considered in which the flaws are removed and the as-designed collapse load is determined.
Only the refined mesh model (weld and lid interior regions) isconsidered.
The mesh is identical to Figure 16 but the coincident nodes along the crack faces are merged.
A Calculation No. 11042-0205
.... .Revision No. 2A R EVA Calculation Page 25 of 874.3.4 Half-Symmetry (3D) Case #1The 3D model is based on the Axisymmetric Case #1. (Analysis results showed that there was negligible difference in the results from Axisymmetric Case #1 and Case #2. The total projected cross-sectional areaof the flaws in Case #1 is greater than Case #2. Therefore, Case #1 is considered critical for the side droploading).
The same flaw pattern is modeled, but the initial mesh is slightly less refined in order to obtain reasonable run times. Mesh sensitivity studies are described below. The half-symmetry model is used for internalpressure loading (as a benchmark case to study the effects of mesh refinement) and also for side-drop loading.The shield plug support ring is connected to the DSC shell at the two corners using nodal DOE couples torepresent the fillet welds used to join the two parts.In order to improve the numerical stability of the ANSYS model, soft springs (COMBIN 14) elements are usedto connect the shield plug to the support ring. The springs have a stiffness of 1 lb/in. The low stiffness combined with the very small relative deflections between these parts results in negligible internal force inthe springs.
The forces in the springs at the final converged solution are reviewed to confirm that the springforces are small.In all load cases, symmetry conditions are applied to the cut face of the model. Axial constraints are appliedat the bottom cut of the DSC shell near the mid-length of the cask (remote from the top end components ofinterest.)
For the internal pressure load case, the model was further reduced to a 90-degree model andsymmetry constraints were placed on both cut faces of the model.The purpose of this calculation is to evaluate the effects of the closure weld flaws and qualify the welds andany other components affected by the welds. All other aspects of the DSC (such as the shell remote fromthe welds) are not in the scope of this calculation.
The modeling approach (loads and boundary conditions) for the side drop event are considered in light of this purpose and are described in the following paragraphs.
For the side drop cases, the OD of the canister shell is constrained in the vertical (drop) direction for a smallsector (approximately 1.5,, inches or 2.8 degrees) of assumed contact.
In reality the DSC is supported insidethe Transfer Cask (TC) during this event. Therefore the true boundary condition would either be a line ofcontact along a TC rail (which is 3" wide) or a line of contact at areas remote from the rails. As deformations
- increase, the area of contact would also increase.
As discussed below in Section 4.4, deflections are over-estimated in a limit load analysis.
Therefore, the area of contact with the TC rail or inner surface is assumedto be constant.
This conservatively neglects the increase in contact area that would occur during the dropdeformations.
Additionally, this boundary condition is representative of the DSC storage condition inside theHSM, where the DSC rests on the 3-inch wide steel rails.As discussed in the Reference 5.2 calculation, the DSC payload (basket and fuel) are located approximately 21.5 inches away from the ITCP and are therefore considered to have no effect on the DSC lid components.
The effect of the basket and fuel loading on the DSC shell is considered in the basket design-basis calculation for side-drop loading.
The basket hold-down ring is.a grid-type structure that does not represent significant weight and is of sufficient strength and stiffness to be self-supporting during the side drop and notsignificantly affect the DSC shell and adjacent regions.
Therefore, as in the Reference 5.2 calculation, theDSC payload is not considered as affecting the top-end components and the weight is applied as a pressurealong a strip of elements at the impact region, beginning approximately 23" below the ITCP. Since the loads A Calculation No. 11042-0205
- ....
No. 2-AR EVA Calculation Page 26 of 87are essentially applied directly over the supported (impact) region of the DSC shell, they have no appreciable effect on the shell deformations.
Images of the Half-Symmetry model are shown in Figure 19 to Figure 21.In order to study the adequacy of the mesh for the half-symmetry model, an internal pressure load case wasperformed and compared to the results of the axisym metric case refined mesh. This study confirms theadequacy of the mesh in the cross-section of the 3D) model. In order to evaluate the mesh in thecircumferential direction, a model was created with a refined mesh in the regions of the model showing largeplastic strains (the impact region) and locations where tensile stress is expected in the weld (at the 90-degree location where the lid resists ovalization of the DSC shell). This model is shown Figure 22.
A Calculation No. 11042-0205
~Revision No. 2AR EVA ~~CalculationPae2o8 4.3.5 Half-Symmetry (3D) Case #0In order to study the effect of the flaws, an additional case is considered in which the flaws are removed fromthe model and the as-designed side drop limit load capacity is determined.
4.4 Limit Load Solution DetailsAs discussed in Section 4.1, this calculation is based on predicting the lower-bound collapse loads of theDSC based on limit load analysis.
All materials are modeled as elastic-perfectly plastic1, with yield stressvalues based on the limit load analysis requirements of the ASME code [Ref. 5.7]. Table 3 lists thetemperatures used for each load case, and the values of the material properties are shown in Table 4 andTable 5.The prescribed loads are applied to the model, and then are increased linearly until the solution fails toconverge.
The analyses use small deflection theory (NLGEOM,OFF).
This is conservative since deflections areunrealistically high in a limit load analysis due to the lower-bound non-strain-hardening material properties that are used. If large deflections were to be considered, the beneficial effects of OTCP and ITCPmembrane action and of increased contact areas would be over-estimated, resulting in non-conservative effects.
This was verified with a sensitivity study using NLGEOM,ON, which resulted in much higher collapsepressures.
This confirmed that using NLGEOM,OFF is appropriate, and conservative.
In addition, Paragraph 5.2.3.1 of Reference 5.18 states that small displacement theory is to be used in a limitload analysis.
I "Elastic-perfectly plastic is standard mechanics of materials term that describes an idealized material that behaves in alinear-elastic manner up to the yield point, and thereafter is perfectly-plastic, i.e. non-strain hardening.
A Calculation No. 11042-0205
- Revision No. 2AR EVA Calculation Page 28of8
75.0 REFERENCES
5.1. AREVA Document No. 180-9236022-000.
NDE Services Final Report. "Monticello, DSC-16,Phased Array UT Examination Results of the Inner and Outer Top Cover Lid Welds."Revision 0.5.2. AREVA (Transnuclear)
Calculation No. NUH61BTH-0200 Revision
- 0. "NUHOMS-61BTH Type 1 Dry Shielded Canister Shell Assembly Structural Analysis."
5.3. ARE VA (Transnuclear)
Drawing No. NUH61 BTH-3000 Revision
- 8. "NUHOMS 61 BTH Type 1DSC Main Assembly."
5.4. AREVA (Transnuclear)
Drawing No. NUH61 BTH-3001 Revision
- 4. "NUHOMS 61 BTH Type 1DSC Shell Assembly."
5.5. AREVA (Transnuclear)
Drawing No. NUH61 BTH-4008 Revision
5.6. ANSYS Version 14.0. ANSYS Inc. (Including the ANSYS Mechanical APDL Documentation).
5.7. ASME Boiler and Pressure Vessel Code,Section III Subsection NB. 1998 Edition withAddenda through 2000.5.8. ARE VA (Transnuclear)
Document Number NUH-003 Revision
- 14. "Updated Final SafetyAnalysis Report for the Standardized NUHOMS Horizontal Modular Storage System forIrradiated Nuclear Fuel."5.9. ASME Boiler and Pressure Vessel Code,Section III Appendices.
1998 Edition with Addendathrough 2000.5.10. ASME Boiler and Pressure Vessel Code,Section XI. Rules for Inservice Inspection of NuclearPower Plant Components.
1998 Edition with Addenda through 2000.5.11. AREVA (Transnuclear)
Document No. NUH61BTH1-0101 Revision
- 2. "Design CriteriaSpecification for the NUHOMS-61BTH Transportable Storage Canister."
5.12. AREVA Calculation No. 11042-0204 Revision
- 3. "Allowable Flaw Size Evaluation in the InnerTop Cover Plate Closure Weld for DSC #16"5.13. Chattopadhyay, Somnath.
"Pressure Vessels Design and Practice."
CRC Press. 2004.5.14. TriVis Incorporated Welding Procedure Specification No. SS-8-M-TN Revision 10.5.15. ASME Boiler and Pressure Vessel Code,Section II, Part C. "Specifications for Welding Rods,Electrodes, and Filler Metals."
1998 Edition with Addenda through 2000.5.16. ASME Boiler and Pressure Vessel Code,Section II, Part D. "Properties."
1998 Edition withAddenda through 2000.5.17. AREVA (Transnuclear)
Calculation No. NUH61BTH-0253 Revision
- 0. "NUHOMS 6IBTHType I DSC Shell Assembly Outer Top Cover Plate Critical Flaw Size of Weld."5.18. ASME Boiler and Pressure Vessel Code,Section VIII Division
- 2. 2010.5.19. AREVA Document No. 54-PQ-1 14-001 Revision
- 0. "Phased Array Ultrasonic Examination ofDry Storage Canister Lid Welds."5.20. NRC Spent Fuel Project Office Interim Staff Guidance
-15. Materials Evaluation.
(ISG-15).
1110/2001.
A Calculation No. 11042-0205
.Revision No. 2AR EVA Calculation Page 29 of 876.0 ANALYSISTable 6 shows a summary of the results of all of the analyses performed for this calculation and includes acomparison of the results with the acceptance criteria.
Each case is discussed in more detail below.6.1 Axisymmetric Analyses for Internal Pressure6.1.1 Axisymmetric Case #1 -Initial Mesh ModelTwo analyses are performed with the Axisymmetric Case #1 initial-mesh model described in Section 4.3.1:one case using the Service Level A/B material properties and one case using the Service Level D materialproperties.
The collapse pressures were determined to be 95.9 psi for Service Level A/B and 136.6 psi forService Level 0. Figure 23 shows various plots of the plastic strain in the initial-mesh model for ServiceLevel A/B at various locations and levels of loading.
These strain plots are also representative of thebehavior of the Service Level D analysis.
Figure 24 shows the deflection history at the center of the lid, andindicates the expected plastic instability that occurs as the limit load is approached.
Note that both thestrains and displacements presented in these figures have no physical meaning and the displacement plotsshow only the loading (pressure) at which the solution fails to converge.
Since the initial mesh contains several element divisions at each critical cross-section, it is not expected thatelement shear locking (due to the default fully-integrated elements) will be significant.
To confirm this, a testcase was done using the Service Level A/B model but with the Simplified Enhanced Strain elementformulation (KEYOP 1=3). The collapse pressure was found to be 96.1 psi, which is essentially identical tothe initial results.6.1.2 Axisymmetric Case #1 -Refined Mesh ModelsAdditional analyses are performed using the Service Level A/B material properties with the refined meshmodels described in Section 4.3.1. Figure 25 and Figure 26 show the plastic strain results for the refinedmesh at the weld region and the refined mesh at the weld and lid interior
- regions, respectively.
The collapsepressures were found to be 94.8 psi and 93.8 psi, respectively, for these models. The OTCP deflection histories are shown in Figure 27. Note that both the strains and displacements presented in these figureshave no physical meaning and the displacement plots show only the loading (pressure) at which the solutionfails to converge.
Figure 28 shows a comparison of the maximum displacement history curves for the various Axisymmetric Case #1 models, done as part of the mesh sensitivity study. As seen in the figure, the results match verywell. The results of the refined mesh models deviate at most (95.9-93.8)/93.8
= 2.2% from the initial meshresults.
This is very close agreement particularly due to the non-linear nature of the analysis.
Therefore, theinitial mesh is considered sufficient.
However since the analysis run times for the axisymmetric cases arereasonable even for the refined mesh model, the remaining axisymmetric cases use a refined mesh.The Axisymmetric Case #1 with refined weld and lids for Service Level D criteria reported a collapsepressure of 132.6 psi.Note that the nodal coupling in the axial direction between the ITCP and OTCP is a valid method to modelthe contact between the plates since the internal pressure loading ensures that the ITCP lid will bear againstthe OTCP, and since the nodes that are coupled remain coincident throughout the analysis, with only very A Calculation No. 11042-0205
, Revision No. 2AR EVA Calculation Page 30of87minor differences in radial position occurring at the later load steps. In order to confirm the behavior of thenodal coupling, the Axisymmetric Case #1 model with refined welds and lids was modified to include contactbetween the ITCP and OTCP. The model replaces the nodal coupling with CONTA171 and TARGE169elements, using the default element parameters.
Figure 29 shows a comparison between the model usingDOE couples and the model using contact elements.
As seen in the figure, the results are very similar, withthe DOF-couple-model showing slightly more conservative results.
Therefore, the nodal coupling isacceptable and is used in all other axisymmetric models.Note that in all of the FEA models, the internal pressure loading was not applied to the faces of the ITCPweld root flaw that is exposed to the internal region of the cask. Pressure loading on this crack face isnegligible since the flaw is only 0.09" high, and in reality the ITOP flaws are generally very short (i.e. not full-circumferential flaws). In order to support this conclusion, a sensitivity analysis is performed where thepressure loading is applied to the ITCP weld root crack faces. The results, shown in Figure 30, confirm thatpressure loading on the faces of this flaw are negligible.
A Calculation No. 11042-0205 Revision No. 2AR EVA Calculation Page 3l of 876.1.3 Axisymmetric Case #2Two analyses are performed with the Axisymmetric Case #2 refined-mesh model described in Section 4.3.2:one case using the Service Level A/B material properties and one case using the Service Level D materialproperties.
The collapse pressures were determined to be 93.7 psi for Service Level A/B and 132.9 psi forService Level 0. Figure 31 shows various plots of the plastic strain for Service Level A/B at various locations and levels of loading.
These strain plots are also representative of the behavior of the Service Level Danalysis.
Note that both the strains and displacements presented in these figures have no physical meaningand the displacement plots show only the loading (pressure) at which the solution fails to converge.
6.1.4 Axisymmetric Case #0One analysis is performed with the Axisymmetric Case #0 refined-mesh model described in Section 4.3.3using the Service Level A/B material properties.
The collapse pressures were determined to be 94.5 psi forService Level A/B. Figure 32 shows various plots of the plastic strain at various locations and levels ofloading.
Note that both the strains and displacements presented in these figures have no physical meaningand the displacement plots show only the loading (pressure) at which the solution fails to converge.
Figure 33 shows a comparison of the maximum center-of-lid displacement history for all three axisymmetric cases. As seen in the figure, there is essential no difference between Axisymmetric Case #0, Case #1 andCase #2. The Case #1 and Case #2 analyses show slightly larger deflections early in the analysis due to theslightly reduced rotational fixity of the welds. However, the final collapse pressure are within (94.5-93.7)/93.7=0.9%
of each other. This supports a supposition that the observed flaws have negligible impacton the governing failure mode of the top end closure plates and welds. Again, displacements from the limitload analyses have no physical
- meaning, other than to show the onset of non-convergence of the FE model.
A Calculation No. 11042-0205
...Revision No. _ __A R EVA Calculation Page 32 of876.2 Half Symmetry Analyses for Internal Pressure (Benchmark Cases)The model described in Section 4.3.4 is used for an internal pressure collapse analysis in order tobenchmark the model against the axisymmetric cases. The collapse pressure was calculated to beapproximately 97 psi. (The run was terminated at 95 psi and the final collapse pressure was estimated toavoid excessive computer run time). Figure 34 shows various plots of the plastic strain at various locations and levels of loading.
A comparison of the half-symmetry case to the refined-mesh axisymmetric case isshown in Figure 35. As seen in the figure, the half-symmetry case closely matches the behavior of therefined mesh axisymmetric model although the results indicate a slightly greater collapse pressure.
Therefore, the half-symmetry model is considered sufficiently accurate for this analysis.
As shown by theresults, and as discussed in Section 7.0, there is significant safety margin available such that further meshrefinement of the half-symmetry model is not warranted.
- However, the effects of circumferential meshdensity for the half-symmetry model can be seen in Section 6.3.1.
A Calculation No. 11042-0205
.... ,Revision No. 2AREVA Calculation Page 33 of 876.3 Half Symmetry Analyses for Side Drop Loading6.3.1 Half-Symmetry Case #1The model described in Section 4.3.4 is used to perform two side-drop limit load analysis.
One caseincludes side-drop acceleration loading only, while the second case includes the DSC off-normal internalpressure of 20 psi. For this later case, the 20 psi internal pressure is applied simultaneously with a 75gacceleration, and then both the pressure and the acceleration are increased linearly until the collapse g-loadis obtained.
(For example, for collapse occurring at 181g, the internal pressure at collapse is 20*181/75
=48.3 psi.)Note that the side drop loading is combined with the design-basis off-normal internal pressure of 20 psi, asopposed to the internal pressure value of 32 psi used for the SL A/B cases which was the sum of the 10 psinormal pressure and an additional 22 psi to account for inertial handling/seismic loads. See Section 4.2The collapse g-load for side-drop-only loading was found to be approximately 181g. The collapse g-loadwhen internal pressure loading was included was found to be greater than 181g. This later run terminated at181g, but based on the collapse behavior (see Figure 41) it is expected that smaller time steps would allowthe solution to continue to larger loads.Various images of the stress and strain in the side drop analyses are shown in Figure 36 to Figure 38. Notethat both the strains and displacements presented in these figures have no physical meaning and thedisplacement plots show only the loading (pressure) at which the solution fails to converge.
The Half-Symmetry Case #1 model with refined mesh in the circumferential direction was used to evaluatethe side drop load case (without internal pressure).
This analysis was performed up until a load of 185 g's,at which time the analysis was terminated manually to avoid large file sizes and excessive run time. As seenin Figure 41, this model showed a greater resistance to the side-drop
- loading, and would eventually result incollapse g-loads in excess of 185 g if smaller timesteps and longer run times were provided.
Images of thestress and strain from this analysis are shown in Figure 39. Note that both the strains and displacements presented in these figures have no physical meaning and the displacement plots show only the loading(pressure) at which the solution fails to converge.
This analysis confirms that the mesh used in the other half-symmetry cases is adequate, and conservative.
A Calculation No. 11042-0205 Revision No. 2AR EVA Calculation Page 34of876.3.2 Half-Symmetry Case #0One side drop analysis is performed with the Half-Symmetry Case #0 model (no flaws) described in Section4.3.5. Based on the results discussed above, only the case without internal pressure loading wasconsidered.
This analysis resulted in a collapse load of 189g. Stress and strain plots from this analysis areshown in Figure 40. As noted previously, these strain plots have no physical meaning.
As shown in Figure41, the collapse behavior was nearly identical to the case with weld flaws, indicating that the flaws hadnegligible effect on the results.6.4 Evaluation of the 25g Corner DropReference 5.2 Section 10.2 evaluated the OTOP weld to resist a 25g inertial load on the entire DSC contentsand neglecting the strength of the ITCP weld. Furthermore, a conservative stress was assumed in the welddue to internal pressure.
The Reference 5.2 calculation is revised below to account for the strength of bothwelds and include a reduction in the weld thickness due to the observed flaws. The total weld thickness istaken as the combined weld throats from the ITCP and OTCP minus the height of the flaws present in thewelds. (See Reference 5.2 Section 10.2 for the basis of the following values and calculations.)
Note that the allowable weld stress noted below includes a joint efficiency factor of 0.7 as described inReference 5.2 Section 6.2. This reduction factor conservatively bounds the reduction factor of 0.8 discussed in Section 3.4. Therefore, no further reduction factor is applied, and the calculation below is conservative.
WTor = 75,81l1bs (total weight of fuel + basket + lids and shield plug)Wp= 68,943 lbs (load due to axial component of pressure)
WTOT2Sg = 25 x 75,811 + 68,943 = 1,964,218 lbsWTOT2EO -1,964,218 947L-b (ethof weld is20.3"W2=- Lweld 208.131 in UOIL3 1tweld + -(0.23 + 0.11) = 0,3475 in (, See Note)wT~g = --- 9,.3 = 27,157 psi (weld shear stress due to 25g corner drop) 0.347'5-T2o0s = 4,120 psi (weld stress radial component due to 20 psi internal pressure)
TTOT -= T2Sg_ +¢ T2opsi -= 27,157 + 4,120 =31,277 psiTALLOW = 32,400 psiTTOT _31,277TALW- = 0.97<1 (0K)T~tw 32,400*Note: the reduction of the weld to account for the flaws is based on the maximum flaw heights in any oneplane through each of the welds. This is taken as 0.23" for the OTCP weld and 0.11" for the IT-CP weld.)Therefore, the top end closure welds, with the observed flaws, are OK for the Service Level 0 corner dropevent.
A Calculation No. 11042-0205 Revision No. 2AR EV/A Calculation Page 35 of 877.0 DISCUSSION AND CONCLUSIONS The lower bound collapse pressure for Service Level A/B criteria was found to be 93.7 psi which is greaterthan the required pressure of 1 .5x32/0.8=60 psi (where 1.5 is the code-required
[Ref. 5.7] factor on the 32psi design pressure loading and 0.8 is the weld strength reduction factor- see Section 4.2). Therefore theService Level A/B criteria is satisfied.
The lower bound collapse pressure for Service Level D criteria was found to be 132.6 psi which is greaterthan the required pressure of 1.1 1x65/0.8=90.2 psi (where 1.11 is the code-required
[Ref. 5.7] factor on the65 psi design pressure loading and 0.8 is the weld strength reduction factor -see Section 4.2). Therefore the Service Level D criteria for internal pressure is satisfied.
As noted in Section 6.1.4 and as shown in Figure 33, there is essentially no difference in the collapsepressure and extremely little difference in the overall collapse behavior of the DSC subjected to internalpressure loading with and without flaws in the weld.The lower bound collapse acceleration for side drop (Service Level D) loading was found to be 181g which isgreater than the required load of 1.1 1x75/0.8=104g.
Therefore the Service Level D criteria for side droploading is satisfied.
As noted in Section 6.3.2 and as shown in Figure 41, there is essentially no difference in the collapse loadand behavior between the as-designed DSC and the DSC with closure weld flaws.The Reference 5.12 and 5.17 calculations document the ITOP and OTOP closure weld critical flaw sizes,respectively, based on the maximum radial stresses in the welds. The guidance and safety factors ofReference 5.10 are used in the critical flaw size analysis.
The critical flaw sizes are determined to be 0.19and 0.29 inches for surface and subsurface flaws, respectively, in the OTCP weld and 0.15 inches forsurface and subsurface flaws in the ITOP weld. The largest single OTOP flaw size documented in Reference 5.1 is 0.14 inches. As discussed in Section 3.2 a very conservative maximum combined flaw height of 0.195inches is assumed in this analysis.
The largest single ITCP flaw size documented in Reference 5.1 is 0.11inches. Therefore, the observed flaws actually are smaller than the critical flaw size limits and therefore it isnot surprising that the flaws are shown to have little effect on the capacity of the structure.
This analysisshows that the quantity and close proximity of some of the flaws also has no significant adverse effects onthe structural capacity of the OSO.Even though all observed flaws in the ITICP and OTCP welds are included in the analysis models usingconservative representations, an additional weld strength reduction factor of 0.8 is considered by increasing the limit load acceptance criteria by a factor of 1/0.8=1.25 times. The 0.8 factor, which is the samemagnitude reduction factor as in ISG-1 5 [Ref. 5.20], conservatively accounts for any additional limitations inthe efficacy of the PAUT examinations and also accounts for the inaccessible area around the vent andsiphon block as well as the geometric reflectors at the root and near the toe of the weld.Therefore it is concluded that Monticello DSC-16, remains in compliance with the ASME Section IIISubsection NB [Ref. 5.7] stress limits with the presence of the ITCP and OTCP closure weld flaws asdocumented in Reference 5.1.
A Calculation No. 11042-0205 A _ _ _ __n N.A REVA Calculation Page 36 of878.0 LISTING OF COMPUTER FILESAnalyses performed on Computer HEA-0213A using ANSYS Version 14.0 [Ref. 5.6].File Date & Time listing is as displayed by the Windows 7 Operating System -Differences may occur duelocal time zone and daylight savings settings.
Analysis Case File Name Date & TimeA~xisymmetric 1 61BTHWeldFlaw 1FAX_2_DETACH.db 4/7/2015 10:45AMInitial Mesh 618THWeldFlaw 1F AX_2_DETACH.rst 4/7/2015 10:45 AMInternal Pressure 61BTH_WeldFlaw 1F AX_2_DETACH.mntr 4/7/2015 10:45 AMSL A/B SOLUTION_AXISYMMIPLimitLoad.INP 4/7/2015 10:20 AMAxisymmetric 1 6IBTHWeldFlaw1F AX 2.DETACH.db 4/7/215 11:59 AMRefined Weld Mesh 61 BTH_WeldFlaw1 1FAX_2_DETACH.rst 4/7/215 11:58 AMInternal Pressure 61BTH_WeldFlaw 1 FAX 2_DETACH.mntr 4/7/215 11:59 AMSL A/B SOLUTION_AXISYMMIPLimitLoad.INP 4/7/2015 11:55 AMAxisymmetric 1 61BTHWeld~law IF AX_2_DETACH.db 4/21/2015 9:04 AMRefined Weld and Lid Mesh 61BTHWeld~law IF AX_2_DETACH.rst 4/21/2015 9:03 AMInternal Pressure 61BTH_Weld~lawIFAX_2_DETACH.mntr 4/21/2015 9:04 AMSL A/B SOLUTION_AXISYMMIPLimitLoad.INP 4/7/2015 11:55 AMAxisymmetric 1 618TH_Weld~law 1FAX_2_DETACH.db 4/20/2015 10:43 AMInitial Mesh 618THWeld~law 1FAX_2_DETACH.rst 4/20/2015 11:09 AMInternal Pressure 618THWeld~lawIFAX_2_DETACH.rnntr 4/20/2015 11:09 AMSL D SOLUTION_AXISYMMIPLimitLoad_SLD.INP 4/7/2015 11:20AMAxisymmetric 1 618TH_WeldFlaw 1FAX_2_DETACH.db 4/30/2015 8:12 AMRefined Weld and Lid Mesh 618THWeldFlaw iF AX_2_DETACH.rst 4/30/2015 8:12 AMInternal Pressure 618TH_WeldFlaw iF AX 2_DETACH.mntr 4/30/2015 8:12 AMSL D SOLUTION_AXISYMM IP LimitLoad_SLD.INP 4/7/2015 12:02 PMAxisymmetric 2 618TH_WeldFlaw_2GAX_2.db 4/21/2015 3:06 PMRefined Weld and Lid Mesh 618THWeldFlaw_2GAX_2.rst 4/21/2015 2:58 PMInternal Pressure 618TH_WeldFlaw_2GAX_2.mntr 4/21/2015 3:06 PMSL A/B SOLUTION_AXISYMMIPLimitLoad.INP 4/7/2015 11:55 AMAxisymmetric 2 618TH_WeldFlaw_2GAX_2.db 4/21/2015 3:10 PMRefined Weld and Lid Mesh 618TH_Weld~law_2GAX_2.rst 4/21/2015 3:10 PMInternal Pressure 61BTH_Weld~law_2GAX_2.mntr 4/21/2015 3:10 PMSL D SOLUTION_AXISYMMIPLimitLoad_SLD.INP P 4/7/2015 12:02 PMAxisymmetric 0 618THWeldFlaw 1FAX_2_DETACH.db 4/21/2015 10:39 AMRefined Weld and Lid Mesh 618THWeldFlawIFAX_2_DETACH.rst 4/21/2015 10:31 AMInternal Pressure 618TH_WeldFlaw IFAX_2_DETACH.mntr 4/21/2015 10:39 AMSL A/B SOLUTION_AXISYMMIPLimitLoad.INP 4/15/2015 11:07 AMIntilAxsmmtic1618TH_WeldFlawIFAX_2_DETACH.db 4/17/2015 5:40 PMIiilMesh with Keyoption 1=3 618TH_Weld~lawIFAX_2_DETACH.rst 4/17/2015 5:40 PMInternal Pressure 618THWeldFlaw iF AX 2 DETACH.mntr 4/17/2015 5:40 PMSLA/B SOLUTION_AXISYMM_IP_LimitLoad.INP 4/16/2015 12:27 PM A Calculation No. 11042-0205 Revision No. 2AR EVA Calculation Page 37 of87Analysis Case File Name Date & TimeAxisymmetric 1 61BTHWeld~law_1FAX_2_DETACH.db 5/19/2015 8:45 AMRefined Weld and Lid Mesh 6IBTHWeldFlaw IF AX 2 DETACH.rst 5/19/2015 8:02 AMInternal
- Pressure, SL N/BITCP/OTCP couples replaced with 61BTH_WeldFlaw_1FAX_2_DETACH.mntr 5/19/2015 8:45 AMContact SOLUTION_AXISYMM_
IPLimitLoad.INP 5/18/2015 5:03 PMAxisymmetric 1 61BTH_WeldFlaw 1F AX_2_DETACH.db 5/18/2015 2:42 PMRefined Weld and Lid Mesh 61BTH_WeldFlaw iF AX_2_DETACH.rst 5/18/2015 1:37 PMInternal
- Pressure, SL NBWith Pressure on ITCP Weld Root 61BTH_WeldFlaw 1F AX 2 DETACH.mntr 5/18/2015 2:42 PMFlaw Surfaces SOLUTION_AXISYMMIPLimitLoad.INP 5/18/2015 1:25 PMHalf Symmetry 1 61 BTH_WeldFlaw_1 GC.db 4/29/2015 2:10 PMInitial Mesh 61BTH_WeldFlawI GC.rst 4/29/2015 4:52 PMInternal Pressure 61BTH_WeldFlaw_I1GC.mntr 4/29/2015 4:52 PMSL NB SOLUTIONHALFSYMLimitLoad.INP 4/29/2015 2:10 PMHalf Symmetry 1 61BTHWeldFlawI GC.db 4/30/2015 8:20 AMInitial Mesh 61BTHWeldFlaw_
1GC.rst 4/30/2015 3:35 PMSide Drop 61BTHWeldFlaw_
1GC.mntr 4/30/2015 3:35 PMSL D SOLUTIONHALESYMSD.INP 4/30/2015 8:21 AMHalf Symmetry 1 61 BTH_WeldFlaw_I GC.db 5/1/2015 6:58 PMInitial Mesh 6IBTH_WeldFlawl 1GC.rst 5/1/2015 4:29 PMSide Drop + Internal Pressure 61BTHWeldFlaw_1GC.mntr 5/1/2015 4:13 PMSL D SOLUTION_HALFSYMSD.INP 4/30/2015 10:26 PMHalf Symmetry 1 61BTHWeldFlaw_I GD_Refined.db 5/6/2015 1:54 PMRefined Circumferential Mesh 61 BTH_WeldFlaw_1 GDRefined.rst 5/6/2015 11:48 AMSide Drop 61 BTH_WeldFlaw 1 GD_Refined.mntr 5/6/2015 11:47 AMSL D SOLUTIONHALFSYMSD.INP 5/5/2015 9:01 PMHalf Symmetry 0 61BTH_WeldFlaw_
1GC.db 5/2/2015 6:53 AMInitial Mesh 61 BTH_WeldFlaw_I GC.rst 5/2/2015 6:53 AMSide Drop 61BTH_Weld Flawl IGC.mntr 5/2/2015 3:37 AMSL D SOLUTION_HALESYMSD.INP 4/30/2015 8:21 AM AAR EVACalculation No. 11042-0205 Calculation Revision No. 2Page 38 of 879.0 TABLES AND FIGURESTable 1 -Summary of Design Basis Load Combinations for the 61 BTH DSC [Ref. 5.8]LdCaeHorizontalfDW Vertical DW Internal
.ExternaL Theanal Lit~ing Ot0e Sea-sireLdCaeDSC
.Fuel DSC Fud!
Presasur Condition
- Loads .oad 'LevelNon-Operational Load Cases0-l~ Fab. Leak Ta .tie_ .....- 14.7piog 70"F -.I55lripax~al Test.FabI" e tTne ....... -1512ps i '~ -* :l 70'F -155 !alp ax ial TestNO-S DSC Upsielatin
-... ..- 701E x -ANO-4DSC Vertica!Lift
--x -.. ... 7O5* -A:FL-i DSC$'Css Filling -. -:Cask --Hydrostatic 12(0F.Caak x x AFL-2 DSC/Cask Fillin .. -Cask -~
Hydrosttic Carl: x x A'FL-S
- -. -Cask --, 'Hyd,,.tatic Hyrotaic IOF'0I Cask -A-FL.-S Xftexs o Decon --- Cask- x Hydrostatic Hydrostatic 12OAP Cask --AFL-6
--Cask x' -Hvdro~tiati i2O0FCaal
--.A"FL-7? Fuel Deck .Semisiroadinr
--. -Cask x. "H-drostatic Hydrosta.cL 2OA l":.Cask
--Note 10 CDanizagrningtoad C.'ises,Hydrostatac DlD-1 DSC Blowdosas
--Cask x lO!15 pdig Hydrosfatic 12OAF Ca-,sk -. A÷D2V1UtD4ng- &~a~ I~~ --ADD-3 HeliurnBackfill -Cask x: 12pss Hyd-ostatic 12(O"F Cask -ADD1-4 Final HelismsBackftlU
--Cask x" ,. pai'" .Hydr'ostatic 120"1= Cask -- -- A"DD2-50OssleCover Plate Weld -~ -Cask x .3. Hydrzositic i2OAF Cask --- ATransfer Traier Loa,'ing2L-1 VerticalXfer to Traier -'Cask x !lO5psaig OA(F Cak --AT2-2 Vertical.eftoTrailer
--: -'Cask: ," l0,bps Ii-- .12O-AFCa-k
-_ .- ATL-S Laydoros Cask X " --
-O"FCasI --- ATI_-4 LnSidown Cask X --10/15!p -: I20F Cask" --ALo)ad Case .Horizental DW' "Vertacl.DW
,Internl
,External Thermal LIafdang Other Servite......___________
DSC- Fuel *DSC F~uel Pressm'e Condanano LoadA -Loads !Tr.ansfer Tofl~rom ISFScITR-4 AxialLo~ad-Cold Cask X --10/lvS~ia
-OAF 1 A~xial .-. ATR-2 Trnaue-se Load- Cold Cask X -- -- l0!15 paig -OAF 'ig Transverse
-A*.TR-3 ,Verticai Lo -Cold -_ask X --10,15 psi. "- OAF i! gVertical
-"ATR-4Oblique Lotd-'Cold Cask X -- -I0,q5 vsig -- O'F %n£Axial
-A4**TR-5 AxialLout-H~ot
.Cask .X ---J -l0.!l~psi
-lO' "!0F lg-Axial
--- A:TR.-6 Transvrse Io ad -Hot Cask: X --si -. 101rF .lg Traa. -- ATR-7 Vertical Load Cask X -10/1 paig -- g Verical -A-iR-8 Oblique Load- Hot Cask X. -L0/15 pag -. 10OAF Z-gAxia1
-~ A+-% g-TraiisI'P.9 25g
.Note 1:14 Note 1.14 2Op-si, -10OAFA' -5g D2________Cofiet DropTR-10 75g Side Drotp!A Note 1 -20 psig -- -1O0"PA 75g 1)TR-1i :Top orBottonsEnd DropaA":
Note: 1. 12 20 psig -1C0XIP' ________
- gEndDrop D
A Calculation No. 11042-0205 Revision No. 2AR EVA Calculation Page 39of87Table 1 (Continued)
-Summary of Design Basis Load Combinations for the 61 BTH DSC [Ref. 5.8]HS.MLOADING
-Cokd -Cask X --I0u15 prig --÷80 ip -- AiLD-2 Nonsul Lo adina -Hot .-Cask X [- -101l5 ssik "- O0'oFCask
+S0Kip --I.D-3 ,Cask X -" -l0/!5 prig ,117=F. ÷80 Kip -ALD-4 Off-Norman L~oadisa-'Co1d C~ask X -- -20prig -0=FCas * +.SO ip FE ,BL.D-50ff-Nomnal osdin -Hot Cads X --2Op.sig -"l0O'f Kip FE BLD-4 Cask X 20- psig -. ÷40 Kip IFE Bw/shded'HSM STOR-AGE fIlorizontiI 3W' "Vetrtical 13W Internal External
,Thermal Other Service'_________________
DSC FPuel lTiel Pressuret' 5 Pr~essuretO)
- Conditiion Loads -Loads L~evel'HSM-I Off-Normsal HSM X .. .. 15psie" --:40'FHSM
-, -" BHSMl-2 NormalStorage HSM X -l5psiF -- 0'.FHSM -: -- A:HSM-30ff-Xormsl HSM X -, -!Spsig -I17'FHSM
-- -B*.HSM-4 ff-Normal Temp. -Fsiiled Furl ,HSM X --20 psig -
-EF :C,HSM-5 Blocked Venf Storage HSY, X. --65/120 pria -117' F --- f-HSM-6 B.V. +.Failed Fuel Stora~e HSM X .-65SI20ptig
-- *
-FE 13_HSMBV~x" HS.M-7Es-rdsquake -Cold -HSM X -10115 ps.F -0O FHS*M -EQ: !Y:HSM-8 Ear hquakaLoadi.ma-Hot HSM X -- -10/I -- 100'FHSM
-- EQ C;D#HSM-9 FloodtLoad (50' H-O)--Co~d HSM X --!0/lS~psig 22 psig 0= FHSM -Flood* CHSM-10 Flood ibad (50' H2O)-Hot ,HSM. X --10;15pss 22,psig 100FHSM -: FloodO CH$* 'ODIGHorizontalfllW VereicalflW Internal External Thermal Other Se-wera" -" DSC Fuel DSC 'Fuel Pre~sure('q Condiion
,Loads- .Lotds .o leelUL-1 Normal Urdoadiag-Cold HSM X --- 10/15 p-ig -0"F HSM -+60 Ksp A,UL-2 NormallUnloading-Hot HSM X --- I0f15 psig -100 F HSM +~60 Ki'p -A-"- HSM X -- -- 10{1 prig -117* F -60 Kip -- AUL-4 Off-Nonnal U~iioading--Cold HSM. X --" 20 poig -0° F HSM +60 FE BU'L-5 Off-'ormal Ltiloadi-Hot X --20 psig -100' F HSM +60Kip FEUL-6 HSM X, --- 20psi.¢ -l!7'F' +60Kip FE B.UL-7 Off. Norm.
11 HSM X --- 20 psig -wf/slsade
-*10 K.p FE CUS~gAccideit Unloodlng
-FEFAot,'
HSM X,-.- 65I120wp,!ig
'-
-S0 Hip FE DRF-1 DSC Refiood --Ca.s: .X 20 psia(ns.x)
Hyvdrostatic 120o'F Cask --13 A Calculation No. 11042-0205 Revision No. 2AR EVA Calculation Page 40of87Table 1 (Concluded)
-Summary of Design Basis Load Combinations for the 61 BTH DSC [Ref. 5.8](Notes for the preceding portions of Table 1)1. 25g and 75g drop acceleration includes gravity effects.
Therefore, it is not necessary to add an additional 1.Og load.2. For Level D events, only maximum temperature case is considered.
(TIhermal stresses are not limited for level D events and maximumtemperatures give mininium allowables).
- 3. Flood load is an external pressure equivalent to 50 feet of water.4. BV = HSM vents are blocked.5. At temperature over lOO°F a sunshade is required over the Transfer Cask Temperatures for these cases are enveloped by the l00OF (withoutsunshade) case.6. As described in Section T.4, this pressure assumes release of the fuel cover gas and 30% of the fission gas. Since unloading requires theUSM door to be removed, the pressure and temperatures are based on the normal (unblocked vent) condition.
Pressure is applied to theconfinement boundary.
- 7. As described in Section T.4, this pressu~re ass~umes release of the fuel cover gas ,and 30% of the fission gas. Although unloading requires theUSM door to be removed, thle pressure and temperatures are based on the blocked vent condition.
Pressutre is applied to the shell, innerbottom and inner and outer top cover plates.S. Not used.9. Unless noted otherwise, pressure is applied to the confinement boundary.
10 psig and 65 psig are applicable to Type 1 DSC, while 15 psig,and 120 psig are applicable to Type 2 DSC.10. Fuel deck seismic loads are assumed enveloped by handling loads.11. Load Cases UL-7 and UL-S envelop loading cases where the stresses due to insertion loading of SO kips are added to stresses due to internalpressure (in reality, the insertion force is opposed by internal pressure).
- 12. The 60g top end drop and bottom end drop ,are not credible events, therefore these drop ,analyses are not required.
- However, consideration of60g end drop and 75g side drop conservatively envelops the effect of 25g corner drop.13. Conservatively based on normal operating pressure times 1.5 to cover future 1OCFR Part 71 requirements.
- 14. A 25g corner drop analysis (300 fr-om horizontal) of 61BTH DSC without support from the TC is to be documented.
- 15. Service Level C is fo the standard seismic event and Service D is for the high seismic event.
A Calculation No. 11042-0205 Revision No. 2AR EVA Calculation Page 41lof 87Table 2 -Internal Pressure in the 61 BTH Type I DSCMaximumDesign PressureMaximumused in Ref. 5.2Calculated adTiDesign Condition Pressure Calculation Reference
[psi] [psi]Normal 7.3 10 Ref. 5.8 Table T.4-16Off-Normal 10.9 20 Ref. 5.8 Table T.4-20Accident 56.1 65 Ref. 5.8 Table T.4-24Table 3 -Maximum Temperatures in the 61 BTH Type 1 DSC ShellMaximum DesignTemperature Design Calculated used in ThisCondition Temperature Calculation Reference
[OF] [OF]Storage 374 500Normal Ref. 5.8 Table T.4-13Transfer 439 500Storage 399 500Off-Normal Ref. 5.8 Table T.4-18Transfer 416 500Storage 611 625Accident
________________________
Ref. 5.8 Table T.4-22Transfer 467 500 A Calculation No. 11042-0205 Revision No. 2AR EVA Calculation Page 42of87Table 4 -Properties of SA-240 Type 304. [Ref. 5.11]Sm SuE Aloal ~ liae Yield Stress for SL A/B Yield Stress for SL DTemp Modulus of Strwabes Yil Strs Utenilae Limit Load Analysis Limit Load Analysis[°F] Elasticity Inest ki tegh(Note
- 1) (Note 2)[s][ksi]
[ksi] [ksi] [ksi]70 28,300 20.0 30.0 75.0 30.0 46.0100 28,138 20.0 30.0 75.0 30.0 46.0200 27,600 20.0 25.0 71.0 30.0 46.0300 27,000 20.0 22.4 66.2 30.0 46.0400 26,500 18.7 20.7 64.0 28.1 43.0500 25,800 17.5 19.4 63.4 26.3 40.3600 25,300 16.4 18.4 63.4 24.6 37.7625 25,175 16.3 18.2 63.4 24.5 37.5700 24,800 16.0 17.6 63.4 24.0 36.8(1) The yield strength to be used in a Limit Analysis for Service Level A and B Loading is 1.5"Sm, per Paragraph NB-3228.1 of Reference 5.7.(2) The yield strength to be used in a Limit Analysis for Service Level D Loading is the lesser of 2.3*Sm and 0.7*Su, perParagraph F-1341.3 of Reference 5.9.
A Calculation No. 11042-0205 Revision No. 2AR EVA Calculation Page 43of87Table 5 -Properties of SA-36. [Ref. 5.111E SmYield Stress for SL A/B Yield Stress for SL DTemp Modulus of Allowable SV Ultimate Limit Load Analysis Limit Load Analysis[°F] Elasticity Srs YilStes Tnle(Note
- 1) (Note 2)[ksi] Intensity
[ksi] Strength
[ksi] [ksi][ksi] [ksi]70 29,500 19.3 36.0 58.0 29.0 40.6100 29,338 19.3 36.0 58.0 29.0 40.6200 28,800 19.3 33.0 58.0 29.0 40.6300 28,300 19.3 31.8 58.0 29.0 40.6400 27,700 19.3 30.8 58.0 29.0 40.6500 27,300 19.3 29.3 58.0 29.0 40.6600 26,700 17.7 27.6 58.0 26.6 40.6625(3) 26,400 17.6 27.2 58.0 26.4 40.4700 25,500 17.3 25.8 58.0 26.0 39.8(1) The yield strength to be used in a Limit Analysis for Service Level A and B Loading is 1.5*Sm, per Paragraph NB-3228.1 of Reference 5.7.(2) The yield strength to be used in a Limit Analysis for Service Level D Loading is the lesser of 2.3*Sin and 0.7*Su, perParagraph F-1341.3 of Reference 5.9.(3) All values are interpolated from the 600 0F and 700 °F values.
A Calculation No. 11042-0205 Revision No. 2AR EVA Calculation Page 44of87Table 6 -Summary of Load Cases and ResultsRequired!
Required Pressure toTep. Aalss esg Pesue Satisfy Limit Load Limit Load Code LimitName Mesh Level Loading Criteria Collapse Pressure Load Criteria[e] Ciei Nt )(NoteS5)
[psi] Satisfied?
[psi][piInitial Inena S00 SL A/B 32 60.0 95.9 YesPressureInternalRefined Welds 500 SL A/B 32 60.0 94.8 Yes.u PressureReieaedndLd ItraEPressure 500 SL A/B 326009.Ye Internal 625 SL D 65 90.2 136.6 YesPressureRefined Welds and Lid Inenl 625 SL D 65 90.2 132.6 YesPressureInternalRefined Welds and Lid 500 SL A/B 32 60.0 93.7 YesN PressureRefined Welds and Lid Internal 625 SL D 65 90.2 132.9 YesPressureRefned eld andLid InternaleiedWls n i Pressure 500 SL A/B 32 60.0 94.5 YeaE :Initial Internal5= (NtS)Pesr 00 SL A/B 32 60.0 97 YesRequired G-Load toTep nlss Required!
Satisfy Limit Load Collapse Code LimitName Mesh Level Loading Criterinayi Design G-Load Criteria G-Load Load Criteria[~][g] (Note 5) [g] Satisfied?
Initial Side Drop 500 SL D 75 104.0 181 -Yes~RefinedCircumferential Mesh Side Drop 500 SL D 75 104.0 185 YesE~(Note 3)"z Initial Side Drop with 50 S e(Note 3) Off-Normal I P 7 0.8~ ~Initial Side Drop 500 SL D 75 104.0 189 YesNotes:1) [Not Used]2) The 97 psi colipase load is estimated!/
extrapolated from the final obtained solution at 95 psi. Excessive run times make more precise results impractical.
- 3) The reported collapse load is conservative
-based on the collapse behavior it is expected that smaller analysis time steps would yield larger collapse loads. Thisiwas deemed impractical due to the long rntime and the large magnavailable.
- 4) The Service Level A/B required pressure of 32 psi is based on the design internal pressure of 10 psi plus an equivalent internal pressure of 22 psi which accounts forother Service Level A loads. See Section 4.25) Includes both the ASME code required factor and the 0.8 weld strength reduction factor to account for examination limitations.
A Calculation No. 11042-0205 Revision No. 2A R E VA Calculation Page 45 of 87 ]Figure 1 -Sketch of the 61 BTH DSC Top End and Transfer Cask from Reference 5.1
A Calculation No. 11042-0205 Revision No. 2AR EVA Calculation Page 47of87Vent/Siphon
~frATDSCFigure 3 -ITCP and OTCP Closure Weld Details from Reference 5.5(Note: Section B-B above is rotated 90 degrees clockwise from Figure 2.)
A Calculation No. 11042-0205 Revision No. 2AR EVA Calculation Page 49of8710+16 v-2+7315+31- 6 F+18+28+29/t //8+11+24+25+262 OTCPWeld Metal(ITCP Not Shown)Figure 5 -OTOP Flaws -Raw Data from Reference 5.1(Numbers indicate the Flaw Number from Reference 5.1)REPRESENTATIVE GROUP FLAW20 (0.32L Lx .07 H, NEAR 19 IND Weld MetalEXTENTS OF, GROUP FLAW REGION3 (0.18 L x .09"H, REMOTE FROM ALL BUT #2)14(./x.7 EOEFO L #2)
CIRCUMFERENTIAL, INTERMITTENT)
OTCP(ITCP Not Shown)Figure 6 -OTOP Flaws -Main Flaw Group Reduced and Bounded A Calculation No. 11042-0205
., Revision No. 2AR EVA Calculation Page 50of87-2+140,195"H, FULL CIRCUMFERENTIAL OTCP(ITCP Not Shown)Weld MetalFigure 7 -OTCP Flaws -Bounding Set #1 for ANSYS Collapse Analysis0.14"I-,Weld Metal-14 ,0.07 'H, FULL CIRCUMFERENTIAL OTCP(ITOP Not Shown)Figure 8 -OTCP Flaws -Bounding Set #2 for ANSYS Collapse Analysis A Calculation No. 11042-0205 Revision No. 2A R EVA Calculation Page 51 of 87/OTCPWeld Metal </--ENVELOPE OF ALL BUT 7 AND 11/MAX INDIV. HEIGHT IS 0.09'MAX INDIV. LENGTH IS 2.09"/ 11 (7.17"L x .09"H, REMOTE FROM ALL OTHERS,-- 7(10.34"Lx 0.11"H, NEAR FLAWS 1 THRU 9)ITCPFigure 9 -ITCP Flaws -Raw Data from Reference 5.1OTCPREPRESENTATIVE GROUP FLAW0.09"H, FULL CIRCUMFERENTIAL Weld MetalITCPFigure 10 -ITCP Flaws -Bounding Flaw Set for ANSYS Collapse Analysis A Calculation No. 11042-0205 Revision No. 2A R EVA Calculation Page 52 of 87ANSYS 14.0MAP. 25 2015~10:39:49 ELEMENTSE1ACET=ITYPE MUMSV -1DIST..53.
912XE -16.9125YE -146.68Z-BUFEERFigure 11 -Overview of the Axisymmetric ModelANSYS 14.0MAR 25 201510:37:29ELEMENTSPowerGiaphics E!ACET-1TYPE MUMZV -1 4 9924'XE -31.9932,YE -190.433Z-BIrEEER iiIFigure 12 -Mesh Details Near the Lid Regions of the Axisymmetric Model(Small differences in the mesh exist amongst the sub-models)
A Calculation No. 11042-0205 Revision No. 2AR EVA Calculation Page 53of87~MAR 25 201510:44 :56ELEMENTSPowerGraphics EFACET-1TYPE MUMZV -1*D1ST-1.27014
- XF "32. 6238Y7F -194.594Z -BUFPERFigure 13 -Mesh Details at the Welds for Axisymmetric Case #1ANYS 14.0MAR 25 201510:45:35ELEMENTSPowe EFACET.ITYPE MUM2V -1*DIST-] .210144238*YP 2 PUFFERPEDGEFigure 14 -Flaw Locations for Axisymmetric Case #1 A Calculation No. 11042-0205 Revision No. 2AR EVA Calculation Page 54of87ANSYS 14./APR21 S.111:.44' 5P er~raphlcs EFACET-iTYPE NtJMZV =1*DIST-3.
78456*xF =30. 6287*YF -193, 906Z-BUJFFER Figure 15 -Refined Mesh (Weld Region) for Axisymmetric Case #1AY =194.043 Figure 16 -Refined Mesh (Weld and Lid Interior Region) for Axisymmetric Case #1 A Calculation No. 11042-0205 Revision No. _ __AR EVA Calculation Page 55of87 .......A9% 21 2Q1514=52=50/
0D2PT2. 5 .=30. 7029*YF 2,-3thFkER S.Figure 17 -Mesh Details at the Welds for Axisymmetric Case #2ANSYS 14.0MR25 2015I0: 38:01PowerGraphics EFACET-ITYPE NUMZV -i*DIST-1.
23285*XF -32.7711*YF -194.*652 Z-BUJEEREDGEFigure 18 -Flaw Locations for Axisymmetric Case #2 A Calculation No. 11042-0205 Revision No. 2AR EVA Calculation Page 56of87ANSYS 14.0APR 21. 201514:38:39ELEMENT SPove rGraphi c sEFACET=ITYPE NUN~YV =-iDIST-53.
812i YE =16.8125ZE =146.88iZ-BUFFEP.
EDGEFigure 19 -Overview of the Half-Symmetry Model A Calculation No. 11042-0205 Revision No. 2A R EVA Calculation Page 57of87t41 =5 1 1 M AY 1a43 1 04=3EPACET I A? .114t ic0 4
- 013.(c))WLid Region Solid Viewwt lw iil (b) Wlid Region MeshFigure 20 -Detail Views and Mesh Plots of the Half Symmetry Model A Calculation No. 11042-0205 Revision No. _ __A R E VA Calculation Page 58 of 87ANSYS 14.0MAY 4 201510:04:39ELEMENTSPowe rtraphi csEFACET-1TYPE NUNXV --.373471YV --. 679348ZV -. 631669*D1ST-35.7666
- XF -12.2625*YF -29.0871*ZB' -167.331A-ZS-41 .1572Z-BUJFFER ANSYS 14.0MAY 4 201510:04:54ELEMENTSPowe rGraphics EFACET=ITYPE NUNXV --.373471YV =-. 679348ZV -.631669*DIST-4.
86713i*X =44.6517*YF -26.3642*Z1 =167.034A-ZS-41 .1572Z-BUJFFER Figure 21 -Isometric Views of Half-Symmetry Model A Calculation No. 11042-0205 Revision No. 2A RE VA Calculation Page 59 of 87AflSYS 14.0NAY 4 201510: 22:15ELEM4ENTS Powe rGraphlc3 EYACET-ITYPE NUNXV --. 290468Yv --.652092Zv -. 700299*YF -30. 5999A-ZS=33.8 517Z-BUFFEJR ANSY$ 14.0NAY 4 201510 22=3TYPE NUNXiv -.290460Yv=.652092 2'. .700299' DIT-4.Th42b
- XF 43.55'TY ;0952ZF 165. 967A-ZS-33.,851 712-BUFFERFigure 22 -Isometric Views of Half-Symmetry Model (Refined Circumferential Mesh)
A Calculation No. 11042-0205 Revision No. 2AREVA Calculation Page 60of8714,0 514373 14.5APR 21 2015~ APR 22 251510:41 12 1042:4152 NJOSAL 50O5710O.
NODAL SOLTIJTON sos -t SITB -TOM,.- U T01ME=65EIMPLEQ2 IA E SPLEQV IAVG)SowerGraph es PeweerGIraphi csEFA'57-.
£5A05T-SAVF.SS.Nst S3M .29 1142 -.84 7490141
- 002675 1542 -.047214.297E.3 Il.O05 246m .020492.892E03 .015736.OO196 .020294.001783 .021423* 002376 % .0419 6m .002 75 m .047214(a) Equivalent Plastic Strain in Weld Region [in/in] (b) Equivalent Plastic Strain in Weld Region [in/in]at 20 psi Internal Pressure at 65 psi Internal Pressure~gTS14.OAllOT 14.0APK TI 2405 APR 22 2025ASi121 0:41:541 0:42:53IIDA ODAL SOLUT1ION
( SS" =18SUB -19S IME..95 .9249 5E 5.94K~PLSQV (AVISEPLQ; A~I 5- rls aphr cr Sewe toaphicoEFACE?=I E ACN'I" 1DX 52.0242 2942 ..51.026759221{ -1.36914 0242 -1.25092,2J 1.1230.22794197 1.4971 ) .555962S 2.24612 .8339432.6204 7 .9729142.99403 1 I.112921: 36918 1.25091(c) Equivalent Plastic Strain in Weld Region [in/in] (d) EQV Plastic Strain in the Cover Plates at 95.9at 95.9 psi Internal Pressureps Figure 23 -Results for Axisymmetrdc Case #1 -Initial Mesh -Service Level A/B(Note that the magnitude of the strains and deflections has no true physical meaning due to the nature of limit loadanalysis)
A Calculation No. 11042-0205 Revision No. 2AR EVA Calculation Page 61 of 874A:4'4AN3iS .1!ZV *13r -.5Z 4U S4 J4.~44 *~1;~~ 1.1254,4 44.744 ,4,.~:4,4.44; iaiii 47.54.4 41.141 44.137TCCrno; f,,i',nar.
Ipti>(a) Service Level A/B Material Properties 4..:4:4:74'CDIS'T-.T5 YF 1.54, !UFFER4 Z~275~ 4,5.71.54 41414 7772574, U4.44~444 47.477 44.:745 44,4417 127.111r tsu. P$1J(b) Service Level D Material Properties Figure 24 -Deflection History of the Center of the OTCP for the Axisymmetrdc Case #1 Initial Mesh(Maximum deflection occurs at the center point of the lids, in the outward axial direction)
(Note that the magnitude of the deflections has no true physical meaning due to the nature of limit load analysis)
A Calculation No. 11042-0205 Revision No. 2A REVA Calculation Page 62of87E[PLE\ tAWS EFPLEQV d .WA-Al l I .0921454AJVAE 4.ET- EII'A -14.2)A .404.171 I IIA .512r~44.
311 -I 11T1414.74114.103 4A~73 1APEV A(c) quialet Plsti Stainin Wld egin [i/in (d EQ PlaticStrin i th Coer latesat 4A.8at 94. psi Iterna Pressre psFigure 25 -Results for Axisym metric Case #1 -Refined Mesh in Weld Region -Service Level A/B(Note that the magnitude of the strains and deflections has no true physical meaning due to the nature of limit loadanalysis)
I A Calculation No. 11042-0205 Revision No. 2AR EVA Calculation Page 63of87FAm.2YS I4. 11A-i SIUB -7AVP, S4LM~ti 4(a) Equivalent Plastic Strain in Weld Region [in/in] (b) Equivalent Plastic Strain in Weld Region [in/in]at 20 psi Internal Pressure at 65 psi Internal PressureI~ *I SrE8E .i41AC flL.V A'05"xl(c) Equivalent Plastic Strain in Weld Region [in/in] (d) EQV Plastic Strain in the Cover Plates atat 93.7 psi Internal Pressure 93.7 psiLevel A/B(Note (c) and (d) are plotted one timestep before the collapse pressure)
(Note that the magnitude of the strains and deflections has no true physical meaning due to the nature of limit loadanalysis)
A Calculation No. 11042-0205
- , Revision No. 2AR EVA Calculation Page 64of8772.117., 1.2 I.* I --.....l. .J! 2157DI*Th -19;XFv I.(a) Service Level A/B Material Properties Refined Mesh at Weld Region Only71 I I .7 40 7,71 1171 O.2T.n7. 9 I; II RI:'l' 19.104 14')POST2.?YFv -,5-VI 11,5~Refined Mesh at the Weld and Lid Interior RegionsFigure 27 -Deflection History of the Center of the OTCP for the Axisymmetric Case #1 Refined Mesh(Maximum deflection occurs at the center point of the lids, in the outward axial direction)
(Note that the magnitude of the deflections has no true physical meaning due to the nature of limit load analysis)
A Calculation No. 11042-0205 Aevi_ _ _ __.AR EVA Calculation Page 65of873 0 -.. ..... .In itial M e sh ... .... .... ... .... .... .... ..... .. ....... l .. ................Initial Mesh, Enhanced ormul atios !I2.50 -.... ..-Refined M esh (W elds) ... .... .... ... ..... ... ... .. ..... ..-ii-- -Refined Mesh (Welds end Lids) =E /15 .-. t."!050 i --.... .... ... -- --..... .........
.... ...... .... ..0.00 --.. ..............
-... .........
..... ...-- ... ------.... ........
.. .. .. ...............
.... --...........
... -.. ...... .... .... .... .020 40 60 80 100 120Internal Pressure
[psi]Figure 28 -Comparison of Maximum Displacement Histories for Axisymmetric Model Sensitivity Studies(Maximum deflection occurs at the center point of the lids, in the outward axial direction)
(Note that the magnitude of the deflections has no true physical meaning due to the nature of limit load analysis)
(Service Level A/B materiaJ Properties)
A Calculation No. 11042-0205
- Revision No. 2AR EVA Calculation Page 66of871.80 i --/ 116 -I -I I~00 -0 2 0 0 5 0 0 8/9 0Internal Pressure
[psi]Figure 29 -Comparison of Maximum Displacement Histories for Axisymmetric Model with Lid ContactDefined using Nodal DOF Couples vs. Contact Elements(Maximum deflection occurs at the center point of the lids, in the outward axial direction)
(Note that the magnitude of the deflections has no true physical meaning due to the nature of limit load analysis)
(Service Level A/B material properties)
A Calculation No. 11042-0205 Revision No. 2A REVA Calculation Page 67 of 87S 4.00 , .. .. .3 .00 * ..... Re fined Mesh (W'elds and Lids) : .....i i RefinedvMesh (Welds andtLids) w/ Pressure onlTCP Crack Faces' i ,-i100i 0 10 20 30 40 ,50 60 70 80 90 1001 Internal Pressure (psi)Figure 30 -Comparison of Maximum Displacement Histories for Axisymmetric Model With and WithoutPressure Loading Applied to the ITCP Weld Root Flaw Faces(Maximum deflection occurs at the center point of the lids, in the outward axial direction)
(Note that the magnitude of the deflections has no true physical meaning due to the nature of limit load analysis)
(Service Level ANB material properties)
A Calculation No. 11042-0205 Revision No. _ __AR EVA Calculation Page 68of87mU.IE,:V R,+f+)LW IA,'.0 -I'.03.0.+g+,
47m +O +3 m .0*(a) Equivalent Plastic Strain in Weld Region [in/in] (b) Equivalent Plastic Strain in Weld Region [in/in]at 20 psi Internal Pressure at 65 psi Internal PressureI Is ; o, + ,o+1U07I :f#C .-+' :+,++++,1ALOOU~VZ E' -!+101(c) Euivaent lastc Stain n Wel Regon [n/in (d)EQV lasti Strin i theCove.Plaes aat 936 ps Intenal ressue 936 psLevel A/B(Note (c) and (d) are plotted one timestep before the collapse pressure)
(Note that the magnitude of the strains and deflections has no true physical meaning due to the nature of limit loadanalysis)
A Calculation No. 11042-0205 Revision No. 2AR EVA Calculation Page 69of87AP.2 0- Apa .1 I. : I OLU U'EEACET' 3a .at2Ps nenlPesr at 64psPIntrna2Presurt t-+ SNX P,.,13¶1I I Ii .t 040!s J(a) Equivalent Plastic Strain in Weld Region [in/in] (b) EquValn Plastic Strain in Whel Roegionte aininat 240 psi Internal Pressure at65pi .nera presurFigure 32 -Results for Axisymmetric Case #0 -Refined Mesh in Weld and Lid Interior Region -ServiceLevel A/B(Note (c) and (d) are plotted one timestep before the collapse pressure)
(Note that the magnitude of the strains and deflections has no true physical meaning due to the nature of limit loadanalysis)
I A Calculation No. 11042-0205 Revision No. 2AR EVA Calculation Page 70of874.003.503.002,500S2.00E1.501.000.50-AAxisynmretric Case #1......................
Axisymnmetric Case #2---- -- Axisymmnetric Case #0 (No Flaws):'Ii IIII/i#I//I /I/0.0010 20 30 40 50Internal Pressure
[psi!60 70 80 90 100IFigure 33 -Comparison of Maximum Center-of-Lid Displacement Histories for the Various Flaw Models(Service Level NB material properties)
(Note that the magnitude of the defiections has no true physical meaning due to the nature of limit load analysis)
A Calculation No. 11042-0205 Revision No. 2A REVA Calculation Page 71lof 87?W(A 4 ZO I M(AY4 2015gF&CET- TI.4 I0 (a) Equivalent Plastic Strain in Weld Region [in/in] ()Euvln lsi tani edRgo i/nat 20 psi Internal Pressureat6psInealPsur (d) QV Pasti Stain n th CovrPates at 9 psFigure 34 -Results for Half-Symmetry Case #1 Internal Pressure Loading Benchmark Analysis
-ServiceLevel NB(Note that the magnitude of the strains and deflections has no true physical meaning due to the nature of limit loadanalysis)
I A Calculation No. 11042-0205 Revision No. 2AR EVA Calculation Page 72of874.(003.503.00'E 2.500IIE 1.50t1.000.50-0.00 1050Internal Pressure
[osij100Figure 35 -Benchmark of the Half Symmetry model with the Axisymmetric Analysis(Service Level NB material properties)
(Note that the magnitude of the deflections has no true physical meaning due to the nature of limit load analysis)
A Calculation No. 11042-0205 Revision No. 2AR EVA Calculation Page 73of87Atl4r'3+3 14, KAY I 402501 11 5t,2 ,1 4 I0+TIKE-(a qivln (o iss trs ps]a 7g () qialn latcStani Wel Regio [in/in3455 93at4475g.
M.4Y ++ :O£,0 1 ll-59 a8"4. .15234(c31iaet vnMss Srs pi at' 18g () qiaen3lstcSrini el4ein i/nE 54at4181g.
Figure 36 -Equivalent Stress and Plastic Strain Plots from the Half-Symmetry
- 1 Side Drop Analysis(Note that the magnitude of the strains and deflections has no true physical meaning due to the nature of limit loadanalysis)
I A Calculation No. 11042-0205 Revision No. 2AR EVA Calculation Page 74of8714.0MAY 4 201510: 55"29N4OfAL SOLOTIOMSTEP-4SU]9 TIMSE 5USUK iAV')RSYS1-0KFACET-1AVP.E3-Mat DM4X -. 340414i SMX -.340414i m -:40414(a) Deformed Shape Plot -Axial View -Exaggerated Scale14.0MAY 4 0= 5V:19~NODAL TIME-1S0, 545USUM lAVG)\ ~Powz'Gt aphi csEFACET-I DMX -404140.0Y624.075047m.113471.226942(b) Deformed Shape Plot of DSC Shell -Axial View -Exaggerated ScaleFigure 37 -Additional Results Plots from the Half-Symmetry
- 1 Side Drop Analysis(Note that the magnitude of the deflections has no true physical meaning due to the nature of limit load analysis)
A Calculation No. 11042-0205 Revision No. 2AR EVA Calculation Page 75of87AIOOY 14.014A1 4 14 15', ¶5314 4 SOLI.1 212TI1975 EPPL2 [AVOI' []~2 1142 97 .01 V2124...I0-1.01* 4 .0 O'D& L I 34.1.1 419(c) quialen (vn Mses)Strss psi]at 81g (b) Equivalent Plastic Strain in Weld Region fin/in]at 18114.Figure 38 -Equivalent Stress and Plastic Strain Plots from the Half-Symmetry
- 1 Side Drop Analysis withOff-Normal Internal Pressure(Note that the magnitude of the strains and deflections has no true physical meaning due to the nature of limit loadanalysis)
I A Calculation No. 11042-0205 Revision No. 2AR EVA Calculation Page 76of87O1514021 4.013t147 01 13434,14 IOI04Q'. A2,07AVE..e 0220,9 -7,io. 3319.02 4441 !,3046. 944402,(a) Equivalent (von Mises) Stress [psi] at 75g (b) Equivalent Plastic Strain in Weld Region [in/in]at 75g.433 ' 14.0 3 149 20MA 3 215 : 53 3OLJT7ON3 L 5?LU 71 TP4 1' , 420m)AVPF .41e3m0403
-.132032343 1 ..45347408 51325 447 1 04.4-1441atU185g.
Figure 39 -Equivalent Stress and Plastic Strain Plots from the Half-Symmetry
- 1 Side Drop Analysis withRefined Circumferential Mesh(Note that the magnitude of the strains and deflections has no true physical meaning due to the nature of limit loadanalysis)
Calculation No. 11042-0205 A _____Revision No. 2AR EVA Calculation Page 77of8714171 14 ..14.14 14 444. 4 2~7J51441 4 211I ~74 ~ ?44XAL SOtU11744 1J74A1. I. 1754?r1141-1114~ -44 141~ s.2za~e1754 5 A 1.444414.75 4114L 1444 .47 '154 1 144. ..A4734~144.-'s 474S442 2~4~ .'7~42743 4 114it*-'724~4.5
-0155 157~l.i-1752~ 34119~.4
-147242-L(b) Equivalent Plastic Strain in Weld Region (in/in](a) Equivalent (von Mises) Stress [ps~ at 75g at 75g.N4221 14.27411211 ~ 14741 4 201114,21 4 2274I15XAL 111011SIXPL 4)L1TIIN Sti5~.45715 ~'l751515457V IAVI471{72~ ~flA1~5't'l 741514.4isi 4.41757-44,75 0117 ~.O11014441 4441 ~,1.41744 1441 -1142.410 11541 '-512) i744)771042.319 4441122122342 40.11-72114 3 2:474271-14414,',
-~-144201.4 4,4444,~. ~ 54~49ti-I(d) Equivalent Plastic Strain in Weld Region (in/in](C) Equivalent (von Mises) Stress [psi] at 189g at 189g.rigure 'iu -c-qulvaieni ano riastic r',ots irom me ra-wr-ymmeiry W.u Uro r-iawsj oiue uropAnalysis(Note that the magnitude of the strains and deflections has no true physical meaning due to the nature of limit loadanalysis)
A Calculation No. 11042-0205 Revision No. 2AR EVA Calculation Page 78of870.800.700.60"2E 0.50.20.40S0.300.200.100.000 20 40 60 80 100 120 140 160 180 200Acceleration
[g]Figure 41 -Comparison of Maximum Displacement Histories for the Various Half-Symmetry Analyses(Service Level D material properties)
(Note that the magnitude of the strains and deflections has no true physical meaning due to the nature of limit loadanalysis)
A Calculation No. 11042-0205 Revision No. 2AR EVA Calculation Page 79of8710.0 APPENDIX A -ADDITIONAL RESULTS PLOTS AND ELASTIC-PLASTIC ANALYSESPurposeThe purpose of this appendix is to document elastic-plastic analyses of DSC-16. The models listed beloware used as a basis for the analyses.
Axisym metric 1 with Refined Welds and Lid MeshHalf Symmetry 1 with Initial MeshThese models produced the bounding results using the limit load methodology.
The models are updated toinclude the elastic-plastic material properties described below. In addition, these new runs consider theeffects of large-deformations (NLGEOM,ON).
The intent of these analyses is to provide a more realistic prediction of the actual material strains that would occur under the design basis loading, as opposed to theover-estimated strains and deformations which result from the limit-load analysis methodology.
Material Properties The elastic-plastic behavior of SA-240 Type 304 stainless steel is idealized using Ramberg-Osgood stress-strain curve constants calculated using the equations in Appendix B of Reference A12.The constants arecalculated using the ASME code [Ref.5.16]
specified minimum yield and ultimate strength values at theapplicable temperatures.
In order to incorporate the curves into the ANSYS analysis, the initial slope of thecurves must match the defined elastic modulus.
Therefore the first data point in the curves is defined at the(strain,stress) data point (Sy/E,Sy).
The material behavior is based on true stress and true strain, since theANSYS analysis accounts for changes in geometry (e.g. necking).
The following equations from Reference Al were used to develop the curves:E = true strainEo = true strain at yield= true stress6o= true stress at yield1n=ln(1 + eu)ln(1 + eu) c'u(1 + eu) IIau(1 + eu)eu= engineering strain at the ultimate tensile strength2 The euations in Reference Al to develop the full-range true stress-strain curve are based on curve fits of tensile testdata. The resulting curve is not indicative of a specific failure type or analysis approach.
Rather, it is a method todevelop a full-range Stress-strain curv of a material using a limited set of data (iLe. minimum specified yield andultimate strengths.)
A Calculation No. 11042-0205 Revision No. 2AR EVA Calculation Page 80 of 87637= engineering yield stressau= engineering ultimate strengthThe value of eu is taken as 0.35, which is assumed to be etot-0.05, where etot is taken as the minimumspecified elongation of the material (40%), per Reference A2.The relationship between true and engineering stress and strain is per the following equations:
6true = aeng(eeng
+ 1)Etrue = ln(eeng + 1)Figure A-I shows both the true and engineering stress strain curves based on the Ramberg-Osgood equations.
The curves at various temperatures as coded into the ANSYS analysis are shown in Figure A-2.The SA-36 shield plugs use a bi-linear stress strain curve with a tangent modulus of 1% of the initial elasticmodulus.
This results in a less stiff representation of the shield plug, which will result in conservatively greater strains in the DSC.Load CasesAnalyses are performed for the following load cases:1. Internal pressure loading (32 psi) for Service Level A/B.2. Internal pressure loading (65 psi) for Service Level D.3. Side drop Loading (75g) for Service Level D.As discussed in Section 4.2, these three load cases bound all of the design loading conditions for the DSCOTCP and ITCP welds.Results and Conclusion Plots of the equivalent plastic strain for the three analyses are shown in Figure A-3 through A-5. The resultsare summarized in Table A-I. As shown by the results, the strain levels remain well below the minimumspecified elongation limits of Type 304 steel and Type 308 weld electrodes
[Ref. A2 and A3]. Therefore, material rupture will not occur at the design conditions.
A Calculation No. 11042-0205 Revision No. _ __AR EVA Calculation Page 81lof 87References Al. EPRI NP-5531.
Evaluation of High-Energy Pipe Rupture Experiments.
January 1988.A2. ASME Section II Part A. Ferrous Material Specifications.
1998 Edition with Addenda through 2000.A3. ASME Section II Part C. Specifications for Welding Rods, Electrodes, and Filler Metals. 1998 Editionwith Addenda through 2000.Cornputer FilesAnalyses performed on Computer HEA-0213A using ANSYS Version 14.0 [Ref. 5.6]File date & time listing is as displayed by the Windows 7 Operating System -Differences may occur due tolocal time zone and daylight savings settings.
Analysis Case File Name Date & TimeElastic-Plastic 61BTH_WeldFlaw_1F A)( 2 DETACH.db 11/29/2015.8:41 AMAxisymmetric 1Refined Lids and Welds 61BTH_WeldFlaw 1FAX_2_DETACH.rst 11/29/2015 8:18 AMInternal PressureSLNAB SOLUTIONAXISYMMIP_500F.INP 11/27/2015 4:03 PMElastic-Plastic 61BTH_WeldFlaw 1FAX_2_DETACH.db 11/27/2015 3:36 PMAxisymmetric 1Refined Lids and Welds 61BTH_WeldFlaw 1FAX_2_DETACH.rst 11/27/2015 3:33 PMInternal PressureSL D SOLUTION_AXISYMMIP_625F.INP 11/19/2015 11:46 AMElastic-Plastic 61 BTH_Weld Flawl 1GC.db 11/29/2015 8:16 AMHalf-Symmetry 1Initial Mesh 61 BTH_WeldFlaw 1GC.rst 11/27/2015 6:26 PMSide DropSL D SOLUTION_HALESYM_SD.INP 11/20/2015 10:08 AMStress-Strain Curve Development Stress-Strain.xls 11/30/2015 10:59 AM A Calculation No. 11042-0205 Revision No. 2AR EVA Calculation Page 82of87Table A-i -Summary of Elastic-Plastic Analysis Results.ValueAnalysis Case Result[in/in]Internal Pressure Equivalent Plastic StrainService Level A at 32 psi Internal 0.0183Axisymmetric Pressure (1.83%)(Note 1) (Note 1)Internal Pressure Equivalent Plastic Strain 009Service Level 0 at 65 psi Internal (5.97%)Axisymmetric PressureSide DropService Level D Equivalent Plastic Strain 0.0609Half-Symmetry at 75g Acceleration (6.09%)Note 1: The 32 psi internal pressure is bounding for Service Levels A and B and includes designinternal pressure of 10 psi plus an additional 22 psi to account for inertial loading of the DSCcontents onto the lid. See Section 4.2 for details.
A Calculation No. 11042-0205 Revision No. _ __AR EVA Calculation Page 83of8788 type 304100.00______
_____ ______ ___ __90.00 ____ _80.00 ______ _____70.00 ___-___60.00 ______40.00 _____ ___ __30.00 ---o-- Engineering
--_20.0010.00 ___ __0,00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40Straln, lnhlnFigure A-I -Ramberg-Osgood Derived Stress Strain Curve for SA-240 Type 304 at 500 °F.
A Calculation No. 11042-0205 Revision No. 2AR EVA Calculation Page 84of87ANSYS 14.91(PIN Thbi-e Eor:
1 NOV 27 201515: 25: 00zeuo 500.002V -1to, 01T1'-. 75Irs 5520.04 .1 .20 .30AZ40V1 14.01(IN TeJ021 Fot Naterial 1 NOV 27 201515:26i05Table flatST1= 500.O00zV *12000 DI S-.,75ZF .5l~o Z-BOtFEr1
.4 -8 2.2 1.6 2.2 .0 11.4q 1.0Figure A-2 -Ramberg-Osgood Stress Strain Curves for SA-240 Type 304 from ANSYS Model at VariousTemperatures.
Upper image shows full range of curves (strain axis from 0 to 40%), lower image shows detail of the initialyield point (strain axis from 0 to 0.2%)
A Calculation No. 11042-0205 Revision No. 2AR EVA Calculation Page 85of87ANSYS 14.0NOV 29 201508:41:45NODAL SOLUTIONSTEP=3SUB =iTIME=32EPPLEQV (AVG)Powe rGr aphi csEFACET=IAVRES=Mat DMX =. 307147SMX =.018306.002034i.004068.012204mm.016272
.018306Figure A-3 -Service Level A Internal Pressure
-Equivalent Plastic Strain at 32 psi *Note*Note: The 32 psi internal pressure is bounding for Service Levels A and B and includes design internal pressure of 10psi plus an additional 22 psi to account for inertial loading of the DSC contents onto the lid. See Section 4.2 for details.
A Calculation No. 11042-0205 Revision No. 2AR EVA Calculation Page 86of87ANSYS 14.0NOV 27 201515:36:02NODAL SOLUTIONSTEP=3SUB =3TIME=65EPFLEQV (AVG)Powe rGraphics EFACET=IAVKE S Ma tDMX ~=.818411 SMX =.0597290.039891.053093Figure A-4 -Service Level D Internal Pressure
-Equivalent Plastic Strain at 65 psi A Calculation No. 11042-0205 Revision No. 2A R EVA Calculation Page 87of87ANSYS 14.0NOV 29 201508:15:52NODAL SOLUTIONSTE P-2SUB -5TIME-75EPPLEQV (AVG)PowerGraphicsAVRES MatDMX =. 080518SMX =. 060912o.0406708
-.046091ANSYS 14.0NOV 29 201508:16:20NODAL SOLUTIONSTEP=2SUB -5TIME=75EPPLEQV (AVG)Powe rGraphics EFACET=1AVR.ES-lat DM( -. 033851SNX =.0548530006095.036569.042664.054853Figure A-5 -Service Level D Side Drop -Equivalent Plastic Strain at 75g.Upper image shows all DSC components, lower image is without shell to allow view of the weld surface.
Thepeak strain of 6.09% occurred on the surface of the shell. Therefore, when the shell was removed for thelower image, the peak strain reported reduced to 5.49%.
L-MT-1 6-0 03Enclosure 2ENCLOSURE 2AREVA, INC.CALCULATION 11042-0205, Revision 2TITLE:61 BTH ITCP and OTCP Closure Weld Flaw Evaluation 87 pages follow Form 3.2-1Calculation Cover Sheet ...Revision No.. 2A RE VA TIP 3.2 (Revision
- 10) Pageli of 87OCR NO (if applicable):
1104-023Revii~n PROECT 61 BTH Type 1 DS~s for114-2 eiin0POETNM:Monticello Nuclear Generating PlantPROJECT NO: 110942 CLIENT: Xcel EnergyCALCULATION
'tiTLE:61IBTH ITCP and OTCP Closure Weld Flaw Evaluation SUMMARY DESCRIPTION:
- 1) Calculation SummaryThis calculation qualifies Monticello DSC-16, a 61BTIH Type 1 DSCI for all design basis loads inconsideration of observed flaws in the Inner Top Cover Plate (ITCP) and Outer Top Cover Plate (OTOP)closure welds.2) Storage Media Location-Coldstor
-lareva_tn/i 1042/11042-0205.000
-Coldstor
-/areva~tn/i11042/110942-0205-002 If original issue, ,is licensing review per TIP 3.5 required?
NIAYes E] No [] (explain below) Licensing Review No.:SoftWare utilized (subject to test requirements of TIP 3.3): Software Version:
Software LogANSYS14.0Revision R-3 1Calculation is complete
..'Date:f " " ' "1120/2016 Originator Name .and Signature:
Jeff Pieper ________Cralcul'ation has been checked for consistency, Completeness, and correctness Date:Checker Name and Signature:
Gabriel LomasCalculation is approved for use HAROON Raheel 6.01 .27 Date:Project Engineer Name and Signature:
13:40:07
-05'00' 1/2712016 A Calculation No. 11042-0205 Aevi_ _ _ __.AR EVA Calculation Page 2of87REVISION SUMMARYAffected AffectedRev. Description Pages Data0 Initial Issue All All1-10, 13,Revised per DCR 11042-022 Revision
- 0. Made 14, 17, 18, Removed1 editorial clarifications, updated information from 21, 22, 24- etaeurevised Reference calculations, removed 34, 36, 45, data.extraneous sensitivity analyses.
46, 48-50,59-771-4, 6-10,13-17, 19-Revied er CR 1042023Revsio 0.22, 24, 27-Revied perlDsRi11042-023 Reviys sion 0pedi. 29, 31, 33 Added files in2 Addedeatic-pastdiscuso anals caictions appendixA 35, 44, 60, Appendix AAddiionl dscusionandclaifictios aded 62-63, 68-throughout.
78,Appendix A(79-87)
A Calculation No. 11042-0205 SRevision No. _ __AR EVA Calculation Page 3of87TABLE OF CONTENTSPage1.0 PURPOSE.......................................................................................................
72.0 ASSUMPTIONS.................................................................................................
83.0 DESIGN INPUT/DATA
.........................................................................................
93.1 DSC Geometry
...........................................................................................
93.2 Flaw Details and Geometry............................................................................
103.2.1 Outer Top Cover Plate.........................................................................
103.2.2 Inner Top Cover Plate ............................................................
.............
113.3 Material Properties......................................................................................
133.4 Design Criteria
..........................................................................................
134.0 METHODOLOGY..............................................................................................
144.1 Analysis Method and Acceptance Criteria............................................................
144.2 Load Cases..............................................................................................
174.3 FEA Model Details......................................................................................
224.3.1 Axisymmetric Case #1 .........................................................................
244.3.2 Axisymmetric Case #2.........................................................................
244.3.3 Axisymmetric Case #0.........................................................................
244.3.4 Half-Symmetry (3D) Case #1 .................................................................
254.3.5 Half-Symmetry (3D) Case #0 .................................................................
274.4 Limit Load Solution Details.............................................................................
2
75.0 REFERENCES
................................................................................................
286.0 ANALYSIS
.....................................................................................................
296.1 Axisymmetric Analyses for Internal Pressure
........................................................
296.1.1 Axisymmetric Case #1 -Initial Mesh Model.................................................
296.1.2 Axisymmetric Case #1 -Refined Mesh Models.............................................
296.1.3 Axisymmetric Case #2.........................................................................
316.1.4 Axisymmetric Case #0.........................................................................
316.2 Half Symmetry Analyses for Internal Pressure (Benchmark Cases)...............................
326.3 Half Symmetry Analyses for Side Drop Loading .....................................................
336.3.1 Half-Symmetry Case #1.......................................................................
336.3.2 Half-Symmetry Case #0.......................................................................
346.4 Evaluation of the 25g Corner Drop....................................................................
347.0 DISCUSSION AND CONCLUSIONS
........................................................................
358.0 LISTING OF COMPUTER FILES............................................................................
369.0 TABLES AND FIGURES......................................................................................
3810.0 APPENDIX A -ADDITIONAL RESULTS PLOTS AND ELASTIC-PLASTIC ANALYSES..............
79 A Calculation No. 11042-0205 Revision No. 2AR EVA Calculation Page 4of87LIST OF TABLESPageTable 1 -Summary of Design Basis Load Combinations for the 61 BTH DSC [Ref. 5.8]......................
38Table 2 -Internal Pressure in the 61 BTH Type I OSO ..........................................................
41Table 3 -Maximum Temperatures in the 61 BTH Type I DSC Shell............................................
41Table 4 -Properties of SA-240 Type 304. [Ref. 5.11] ............................................................
42Table 5 -Properties of SA-36. [Ref. 5.11] .........................................................................
43Table 6 -Summary of Load Cases and Results ...................................................................
44Table A-I -Summary of Elastic Plastic Analysis Results ........................................................
82 A Calculation No. 11042-0205 A Revision No. 2AR EVA Calculation Page 5of87LIST OF FIGURESPageFigure 1 -Sketch of the 618BTH DSC Top End and Transfer Cask from Reference 5.1 ......................
45Figure 2 -Details of the 618TH Top End Component Interfaces
.... ..........................................
46Figure 3 -ITCP and OTCP Closure Weld Details from Reference 5.5..........................................
47Figure 4 -DSC Top End Detailed Dimensions....................................................................
48Figure 5 -OTCP Flaws -Raw Data from Reference 5.1.........................................................
49Figure 6 -OTCP Flaws -Main Flaw Group Reduced and Bounded............................................
49Figure 7 -OTCP Flaws -Bounding Set #1 for ANSYS Collapse Analysis.....................................
50Figure 8 -OTCP Flaws -Bounding Set #2 for ANSYS Collapse Analysis
.....................................
50Figure 9 -ITCP Flaws -Raw Data from Reference 5.1..........................................................
51Figure 10 -ITCP Flaws -Bounding Flaw Set for ANSYS Collapse Analysis
..................................
51Figure 11 -Overview of the Axisymmetric Model .................................................................
52Figure 12- Mesh Details Near the Lid Regions of the Axisymmetric Model ...................................
52Figure 13 -Mesh Details at the Welds for Axisymmetric Case #1...............................................
53Figure 14 -Flaw Locations for Axisymmetric Case #1 ...........................................................
53Figure 15 -Refined Mesh (Weld Region) for Axisymmetric Case #1 ...........................................
54Figure 16 -Refined Mesh (Weld and Lid Interior Region) for Axisymmetric Case #1 .........................
54Figure 17- Mesh Details at the Welds for Axisymmetric Case #2...............................................
55Figure 18 -Flaw Locations for Axisymmetric Case #2 ...........................................................
55Figure 19 -Overview of the Half-Symmetry Model ...............................................................
56Figure 20 -Detail Views and Mesh Plots of the Half Symmetry Model.........................................
57Figure 21 -Isometric Views of Half-Symmetry Model ............................................................
58Figure 22 -Isometric Views of Half-Symmetry Model (Refined Circumferential Mesh) .......................
59Figure 23 -Results for Axisymmetric Case #1 -Initial Mesh -Service Level A/B ............................
60Figure 24 -Deflection History of the Center of the OTCP for the Axisymmetric Case #1 Initial Mesh......61 Figure 25 -Results for Axisymmetric Case #1 -Refined Mesh in Weld Region -Service Level A/B......62 Figure 26 -Results for Axisymmetric Case #1 -Refined Mesh in Weld and Lid Interior Region -ServiceLevel A/B.............................................................................................
63Figure 27 -Deflection History of the Center of the OTCP for the Axisymmetric Case #1 Refined Mesh....64 Figure 28 -Comparison of Maximum Displacement Histories for Axisymmetric Model Sensitivity Studies.
65Figure 29 -Comparison of Maximum Displacement Histories for Axisymmetric Model with Lid ContactDefined using Nodal DOF Couples vs. Contact Elements.......................................
66Figure 30 -Comparison of Maximum Displacement Histories for Axisymmetric Model With and WithoutPressure Loading Applied to the ITCP Weld Root Flaw Faces .................................
67Figure 31 -Results for Axisymmetric Case #2 -Refined Mesh in Weld and Lid Interior Region -ServiceLevel A/B.............................................................................................
68Figure 32 -Results for Axisymmetric Case #0 -Refined Mesh in Weld and Lid Interior Region -ServiceLevel A/B.............................................................................................
69Figure 33 -Comparison of Maximum Center-of-Lid Displacement Histories for the Various Flaw Models ..70Figure 34 -Results for Half-Symmetry Case #1 Internal Pressure Loading Benchmark Analysis
-ServiceLevelA/B
.............................................................................................
71Figure 35 -Benchmark of the Half Symmetry model with the Axisymmetric Analysis.........................
72Figure 36 -Equivalent Stress and Plastic Strain Plots from the Half-Symmetry
- 1 Side Drop Analysis
....73Figure 37 -Additional Results Plots from the Half-Symmetry
- 1 Side Drop Analysis
.........................
74Figure 38 -Equivalent Stress and Plastic Strain Plots from the Half-Symmetry
- 1 Side Drop Analysis withOff-Normal Internal Pressure.......................................................................
75Figure 39 -Equivalent Stress and Plastic Strain Plots from the Half-Symmetry
- 1 Side Drop Analysis withRefined Circumferential Mesh .....................................................................
76 A Calculation No. 11042-0205 Revision No. 2AR EVA Calculation Page 6of87Figure 40 -Equivalent Stress and Plastic Strain Plots from the Half-Symmetry
- 0 (No Flaws) Side DropAnalysis
..............................................................................................
77Figure 41 -Comparison of Maximum Displacement Histories for the Various Half-Symmetry Analyses....78 Figure A-I -Ramberg-Osgood Derived Stress Strain Curve for SA-240 Type 304 at 500 °F................
83Figure A-2 -Ramberg-Osgood Stress Strain Curves for SA-240 Type 304 from ANSYS Model at VariousTemperatures
.......................................................................................
84Figure A-3 -Service Level A Internal Pressure
-Equivalent Plastic Strain at 32 psi..........................
85Figure A-4 -Service Level D Internal Pressure
-Equivalent Plastic Strain at 65 psi..........................
86Figure A-5 -Service Level D Side Drop -Equivalent Plastic Strain at 75g.....................................
87 A Calculation No. 11042-0205
.Revision No. 2AR EVA Calculation Page 7 of 871.0 PURPOSEThe purpose of this calculation is to evaluate DSC-16 at the Monticello Nuclear Generating Plant (MNGP)per ASME Section III criteria in consideration of flaws observed in the Inner and Outer Top Cover Plate(ITCP and OTCP) closure welds. The flaws are documented in the Reference 5.1 Phased Array Ultrasonic Testing (PAUT) inspection report. The canister is a 61BTH Type 1 design. The ASME Section IIISubsection NB Code limits on primary stress are evaluated using the limit load analysis criteria prescribed inthe Code [Ref. 5.7].
A Calculation No. 11042-0205
.... ,Revision No. 2AR EVA Calculation Page 8of572.0 ASSUMPTIONS
- 1. The ITCP weld to the siphon and vent block and the welds of the siphon and vent port cover platesare inaccessible for PAUT inspection.
Approximately 11" are obscured due to the location of thesiphon and vent block. Whereas the main circumferential lid-to-shell welds are made with anautomated welding machine, some manual welding was performed around the siphon/vent block andports. As discussed in Section 3.4, a strength reduction factor of 0.8 is considered for both the ITCPand OTCP welds. This factor accounts for the siphon and vent block welds and uncertainties in theUT technique.
Note that the bounding flaws evaluated in this analysis are treated as full circumferential flaws. Inother words, it is not assumed that the siphon and vent block is free of flaws, but rather contains thesame bounding flaws as the examined welds. The geometry of the siphon and vent block is notassumed in this analysis.
It is assumed that the stresses in the circular configuration bound thestresses that would be computed for a configuration that explicitly includes the siphon and vent block.2. The longitudinal seams in the canister shell caused attenuation in the PAUT energy beam atlocations 24.3" to 24.8" and 129.5" to 130" [Ref. 5.11 that can potentially diminish the effectiveness ofthe examination in these half inch areas. These regions are considered limited examination zones. Itis assumed that the flaws observed outside of these regions are representative, and that no larger ormore bounding flaws exist in the regions behind the canister seam welds. The use of the 0.8 weldstrength reduction factor discussed above in Assumption 1 accounts for any uncertainty in this region.3. [Not used]4. The flaws are considered to be planar cracks lying on circumferential planes, parallel with thelongitudinal axis of the cask. (I.e. the crack tips are pointed in the axial directions of the cask). Thisis a conservative flaw orientation since the welds primarily resist normal stresses in the plane of thelids due to plate bending caused by DSC internal pressure.
Also, during the side drop loading,normal stresses in the plane of the lids resist the ovalizing mode of shell deformation.
This flaw orientation is also conservative for through-thickness shear stresses in the lid welds since itmaximizes the reduction in available shear area. (A flaw of equal length, but placed at an angle,would result in less reduction of the weld throat thickness).
- 5. Many of the flaws identified in the Reference 5.1 PAUT examination report lie in very similar locations within the weld cross section.
As discussed in detail in Section 3.2, flaws that lie in similar radial andaxial positions within the weld are considered bounded by a representative "group flaw." Thelocations and sizes of the "group flaws" are chosen conservatively to ensure they are bounding of theindividual flaws.6. The analysis is based on the nominal dimensions of the components as shown in the designdrawings
[References 5.3 and 5.4] including the as-fabricated radial gap between the outer diameterof the lids and the inner diameter of the DSC shell. Although weld shrinkage will close this gap duringclosure operations, the resulting compressive load path between the lids and shell is conservatively ignored.
Further discussion is provided in Section 4.3.
A Calculation No. 11042-0205
.... .Revision No. 2A.R EVA Calculation Page 9of873.0 DESIGN INPUT/DATA 3.1 DSC GeometryThe 61BTH Type 1 DSC geometry is detailed in the Reference 5.3 and 5.4 drawings.
The Reference 5.5drawing shows the details for the final ITCP and OTOP closure field welds. Sketches from Reference 5.1and details from References 5.3 and 5.4 are shown in Figure 1 through Figure 4.The material for all structural components (DSC Shell, OTCP, and ITCP) is SA-240 Type 304 stainless steel.The shield plug material is SA-36 carbon steel.The DSC shell is 0.5" nominal thickness.
The ITCP is 0.75" nominal thickness.
Per the Reference 5.5 drawing, it is welded to the DSC shell andvent/siphon block with a 3/16" groove weld. However, the ITCP lid groove (weld prep) is 0.25" minimum, andit was confirmed that the weld is also 0.25" [Ref. 5.1].The OTCP is 1.25" nominal thickness.
It is welded to the 050 shell with a 1/2" groove weld.The ITOP and OTCP closure welds (with the exception of the ITCP welds around the vent/siphon block andthe welds of the vent and siphon port cover plates) are made using the GTAW process with an automated welder. This is a non-flux type of weld. The vent/siphon block and the vent and siphon port cover platewelds are performed
- manually, also using a non-flux process.
A Calculation No. 11042-0205 Revision No. 2AR EVA Calculation Page 10of873.2 Flaw Details and GeometryVarious sets of bounding flaws are chosen for the detailed analyses based on the flaw dimensions inReference 5.1 and the discussion below. Note that flaws are identified in this calculation using the numerical flaw listings in the Reference 5.1 inspection report.3.2.1 Outer Top Cover Plate3.2.1.1 Case 1Figure 5 shows all of the OTCP weld flaws from Reference 5.1 plotted on an outline of the DSC geometry.
Figure 6 shows a similar plot but with the main cluster of flaws bounded by a box, and showing arepresentative "group flaw" for this region. The longest flaw within the group region is 31 .7" long and thetallest flaw is 0.14" high. Therefore, the bounding flaw for this region is taken as a full circumferential flaw,0.14" in height.Note that all flaws in the group region were reviewed to ensure that no two flaws in close circumferential proximity, considered as being joined, would produce a taller flaw. For example, OTCP Flaw #9 and OTOPFlaw #10 are within 0.17" of each other in the circumferential direction, but their combined height is only0.47-0.38=0.09".
Therefore these flaws, considered
- combined, are bounded by the 0.14" high group flaw.The radial and axial positions of the bounding flaw were chosen to be at the center of the group region. Thisradial position is within the critical failure plane of the weld (i.e. a plane containing the minimum weld throatthickness of 0.5").Figure 6 also shows additional information about the flaws outside of the group region. OTCP Flaw #2 isintermittent around the entire circumference of the DSC. Therefore this flaw, at 0.12" in height, is considered a full circumferential flaw. Since OTCP Flaw #14 is in close proximity to Flaw #2, it is conservatively considered joined to OTOP Flaw #2, and the combined flaw height is considered to be present around theentire circumference.
The combined flaw height is determined based on the geometry to be 0.195".As seen in Figure 6, OTCP Flaw #20 is remote from the group region and from OTOP Flaw #s 2 and 14.OTCP Flaw #20 is only 0.32" in length, and only 0.07" in height. This flaw is separated from OTCP Flaw #19by 0.36" in the circumferential direction and by 0.19" in the axial direction.
It is separated from OTOP Flaw#21 by 1.66" in the circumferential direction and by 0.23" in the axial direction.
Since extension of the flawsunder the postulated loading is negligible (since only one cycle of the critical loads is applied) this flaw willnot join with the adjacent flaws. Additionally, since OTOP Flaw #20 is much smaller than the critical surfaceflaw size of 0.29" from Reference 5.17, it is not considered explicitly in the FEA analyses and is considered bounded by the other modeled flaws which are very conservative.
Similarly, OTCP Flaw #3 is remote from all flaws with the exception of OTCP Flaw #2. However, OTCP Flaw#3 is very small, only 0.18" long and 0.09" tall. Inspection of the PAUT plots (see Page 22 of Reference 5.1)also shows that OTCP Flaw #2, which is considered as fully continuous in this analysis, is actually veryintermittent at the circumferential position of OTCP Flaw #3. Furthermore, OTCP Flaw #3 is much smallerthan the critical subsurface flaw size of 0.29" from Reference 5.17. Therefore, it is not considered explicitly in the FEA analyses and is considered bounded by the other modeled flaws which are very conservative.
Figure 7 shows the first bounding flaw set considered for the OTCP in the ANSYS collapse analyses.
A Calculation No. 11042-0205 S ,Revision No. _ __AIREVA Calculation Page li of 873.2.1.2 Case 2The discussion above and the flaw locations shown in Figure 5 through Figure 7 are based primarily on thetabulated flaw data from Reference 5.1. Since OTCP Flaw #2 is intermittent around the circumference of theweld, a closer inspection of the PAUT scan images is performed, and an additional flaw set for the OTCP iscreated.
In this additional case, the location of OTCP Flaw #2 is based on the PAUT scan image of the flawat the circumferential position of OTCP Flaw #14, which is the only additional flaw that could be considered to interact with OTCP Flaw #2. Based on the PAUT scan images, the flaws are located as seen in Figure 8.In this case the height of both Flaw #2 and Flaw #7 are estimated based on the PAUT scan images and areconservatively larger than the flaw heights tabulated in Reference 5.1.3.2.2 Inner Top Cover PlateFigure 9 shows all of the ITCP weld flaws from Reference 5.1 plotted on an outline of the DSC geometry.
Allbut two of the flaws are clustered in the region of the weld root at the inner surface of the DSC shell. Figure10 shows the bounding flaw set considered for the ITCP in the ANSYS collapse analyses.
Both therepresentative group flaw and ITCP Flaw #7 are considered to be full circumferential flaws. ITCP Flaw #11is remote from all other flaws (in the circumferential direction) and is therefore considered bounded by therepresentative group flaw. The representative group flaw for the ITCP is conservatively placed at the tensionside of the weld when resisting internal pressure.
All of the ITCP flaws documented in Reference 5.1 were reviewed to ensure that no two (or more) flaws,which are in close proximity to each other, could be considered as combined and therefore creating a morecritical flaw. The following cases are considered in particular:
- ITCP Flaw #2 and Flaw #3 are within 0.12" from each other in the circumferential
- position, but theirmaximum combined height (1.58-1.49
= 0.09") is bounded by the group flaw height of 0.09".*ITCP Flaw #5 and Flaw #8 partially overlaps with Flaw #6 in the circumferential direction and wouldhave a combined height of 0.15". However, Flaw #5 (0.15" in length) and Flaw #8 (0.14" in length)are extremely small. Due to their overlap in the circumferential direction, their combined length wouldbe only 0.16", and therefore would not affect the global or local stability of the weld. This very shortregion with a potential 0.15,, high flaw is bounded by the full-circumferential representation of themodeled flaws.*ITCP Flaw #10 is within 0.04" of Flaw #12 in the circumferential direction.
The individual flaws are0.05" tall and 0.04" tall, respectively, and 0.49" long and 0.18" long, respectively.
They are alsoseparated in the axial direction by 0.09". Postulating a flaw from the bottom of Flaw #12 to the top ofFlaw #10 would imply a height of 0.18". However, the combined-height region would be over a veryshort length and would not affect the global or local stability of the weld. Therefore this postulated combined flaw is considered bounded by the full-circumferential representation of the modeled flaws.It is noted that based on Figure 9 and Figure 10, ITCP Flaw #7 appears to be in the base metal of the innertop cover plate. It is likely that the flaw is actually at the fusion / heat affected region between the weld metal A Calculation No. 11042-0205
- , Revision No. 2A REVA Calculation Page 12of87and the base metal. The ANSYS models used in this calculation place the flaw at 0.81" inward from theouter surface of the DSC shell whereas the tabulated data in Reference 5.1 places the flaw at 0.80" from theouter surface.
The 0.01" discrepancy is considered negligible.
The exact location of the flaw is notconsidered critical in light of the significant margin that is available (See Section 7.0) and the generally veryconservative idealization of the flaws (i~e. full circumferential).
A Calculation No. 11042-0205
,,, Revision No. 2A REVA Calculation Page 13of873.3 Material Properties The material properties for the DSC structure are taken from Reference 5.11. The properties of the twomaterials of construction, SA-240 Type 304 stainless steel and SA-36 carbon steel, are provided in Table 4and Table 5, respectively.
The weld metal is considered to be composed of the same properties as the basemetal, as the welds are made with the non-flux GTAW method [Reference 5.14] using bare metal ER308(stainless) filler material.
The tensile strength of the ER308 electrode (80 ksi at room temperature
[Ref.5.15]) is slightly greater than the type 304 base metal (75 ksi at room temperature
[Ref. 5.16]). The yieldstress value of the weld metal is assumed to be equal to or greater than the base metal. Therefore, thetreatment of the weld metal as being identical to the base metal is appropriate for the Section III limit loadanalyses performed in this calculation.
Temperatures used for material properties are discussed in Section 4.2 and are shown in Table 3.Poisson's ratio for all modeled parts is taken as 0.29.Weight density for SA-240 Type 304 is taken as 0.285 lb/in3.Weight density for SA-36 is taken as 0.284 lb/in3.3.4 Design CriteriaAll of the applicable design basis loading conditions are considered in accordance with the requirements ofASME Section III Subsection NB [Ref. 5.7]. Section 4.1 details the methods used to perform the code [Ref.5.7] qualifications.
Section 4.2 details the selection of the bounding load cases.The mockup used in the PAUT process development contained weld manufacturing flaws intentionally distributed in locations that would be expected with the weld process used for the OSO lid closure welds.Approximately 30% of those flaws were placed at the weld root and 27% were placed near the weld toe todemonstrate that they could be reliably detected in the presence of typical geometric responses from thoseregions.
The flaws include incomplete root penetration, lack of fusion, and tungsten inclusions.
AREVAdocument 54-PQ-114-001
[Ref. 5.19], Section 8.0, provides images of the UT responses for these flaws anddemonstrates that the PAUT process can effectively detect these flaws. Furthermore, the qualification performed on the blind mockup provides objective evidence that detection of flaws in these regions of theweld is not a problem.
The blind mockup used for qualification contained a similar percentage/number anddistribution of flaws as the development mockup. Although the flaw information for the blind mockup cannotbe disclosed in order to preserve the security of the mockup for future qualifications, EPRI and NRC3personnel present at the demonstration have reviewed that information.
In addition, uncertainties in thePAUT examination are accounted for by using a 0.8 reduction factor on the limit load. This factor, which is inagreement with ISG-1 5 [Ref. 5.20], conservatively accounts for any additional limitations in the efficacy of thePAUT examinations and also accounts for the inaccessible area around the vent and siphon block as well asthe geometric reflectors at the root and near the toe of the weld.
A Calculation No. 11042-0205 Revision No. 2AR EVA Calculation Page 14 of 874.0 METHODOLOGY 4.1 Analysis Method and Acceptance CriteriaThe 61BTH DSC including the ITCP and OTCP welds are designed and analyzed per ASME Section IIISubsection NB (the Code) [Ref. 5.7] in the Reference 5.2 calculation.
The presence of the ITCP and OTCPweld flaws will cause high local stresses and complex stress fields that will render an elastic analysis (suchas those performed in Reference 5.2) very difficult.
Therefore, the flaws are explicitly included in the finiteelement models as "design features",
and the applicable ASME code [Ref. 5.7] stress limits are evaluated asdescribed below.Primary Stress LimitsIn order to satisfy the primary stress limits of Reference 5.7 paragraphs NB-3221 .1, NB-3221 .2, and NB-3221.3, a Limit Analysis will be performed per Paragraph NB-3228.1.
The acceptance criterion is that thespecified loadings not exceed two-thirds of the lower bound collapse load, as determined using an ideallyplastic (non-strain hardening) material model, with the yield stress set at a value of 1 .5*Smn. This criterion isused for evaluation of the Service Level A and B load cases discussed in Section 4.2.Note that Service Level C acceptance criteria are generally 20% greater than Service Level A criteria, perParagraph NB-3224 of Reference 5.7. This information is used in the discussion in Section 4.2 to eliminate some non-critical load cases.For the Service Level D loadings (accident level internal pressure and side drop), the rules of ASME SectionIll Appendix F Paragraph F-1 341 .3 [Ref. 5.9] are used, which indicate that the loads "shall not exceed 90%of the limit analysis collapse load using a yield stress which is the lesser of 2.3Sm and 0.7Su." This criterion is used for evaluation of the Service Level D load cases discussed in Section 4.2.An additional increase factor of 1/0.8=1.25 is applied to the required limit load collapse pressure in order toaccount for the weld strength reduction factor of 0.8 to account for UT sensitivity and inaccessible weldregions discussed in Section 3.4. Typically, the weld strength reduction factor is applied to the weldallowable stress during qualification.
In the case of limit-load
- analysis, reduction of the material yield stressis applicable.
The reduction in yield stress would have a direct, 1:1 correlation to the calculated lower boundcollapse pressure due to the perfectly-plastic (i.e. non-strain hardening) material model. In this analysis, rather than decrease the material yield stress the required calculated collapse pressure is increased by thefactor of 1.25.Note that the Service Level D criterion is essentially 2.1 times greater than the Service Level A/B criterion, ascalculated below. This information is used in the discussion in Section 4.2 to eliminate some non-critical load cases.At a temperature of 500 °F, the limit load yield stress for SA-240 Type 304 for Service LevelsA/B and D are 26.3 ksi and 40.3 ksi, respectively.
The code [Ref. 5.7] required factors against the lower bound collapse load as determined bythe limit load analyses for Service Levels A/B and D are 1.5 and 1.11, respectively.
(40.3X1.E)
The ratio of the acceptance criteria is therefore:
(2.,.)= 2.1.(i.e. the Service level A/B criteria are 2.1 times as severe)
~Calculation No. 11042-0205
,. .:,Revision No. 2AREVA Calculation Page 15of87Limit Load Analysis Backaqround ASME Section III Subsection NB provides only a basic description of the Limit Load analysis technique.
Amore thorough description is provided in ASME Section VIII Division 2 Paragraph 5.2.3 [Ref. 5.18]:Limit-load analysis addresses the failure modes of ductile rupture and the onsetof gross plastic deformation (plastic collapse)...
Limit-Load analysis provides an altemnative to elastic analysis and stresslinearization and the satisfaction of prim ary stress limits...
Displacements and strains indicated by a limit analysis have no physicalmeaning.Limit load analysis is based on the theory of limit analysis that defines a lowerbound to the limit load of a structure as the solution of a numerical model withthe following properties:
- 1. The material model is elastic-perfectly plastic with a specified yield strength.
- 2. The strain-displacement relations are those of small displacement theory.The limit load is the load that causes overall structural instability.
This point isindicated by the inability to achieve an equilibrium solution for a small increasein load (i.e. the solution will not converge).
Separately, in order to address questions on the potential for material rupture due to potentially high plasticstrains, supplemental elastic-plastic analyses are performed in Section 10.0 (Appendix A).Primary Plus Secondary Stress LimitsThe Code [Ref. 5.71 also prescribes limits on primary plus secondary stresses for Service Levels A and B[Ref. 5.7 Paragraph NB-3222.2].
Secondary stresses may be developed in the DSC due to differential thermal expansion of the interconnected parts and thermal gradients within the structure.
The code stresslimit for primary plus secondary stress (calculated on an elastic basis) is 3Sin. However, as shown in Ref. 5.7Figure NB-3222-1, rules for exceeding the 3Sin limit are provided in Paragraph NB-3228.5, which states that"the 3Sm limit .., may be exceeded provided that the requirements of (a) through (f) below are met."Requirement (a) states that "the range of primary plus secondary membrane plus bending stress intensity, excluding thermal bending stresses, shall be < 3Srn." This provision is related to the potential for "plasticstrain concentrations" occurring in "localized areas of the structure",
and the potential for theseconcentrations to affect the "fatigue
- behavior, ratcheting
- behavior, or buckling behavior of the structure"
[Ref.5.7 Paragraph NB-3228.1].
Requirements (b) through (d) are also limitations related to fatigue and thermalstress ratchet.
As detailed in Section 10.5 of Reference 5.2, the DSC is exempt from fatigue analysisrequirements since all of the criteria in NB-3222.4 of Reference 5.7 are satisfied.
Similarly, since the DSCthermal loads are not cyclic in nature (other than small daily and seasonal fluctuations),
thermal stressratchet is not a concern.
Therefore, the 3Sm limit as it relates to fatigue is not applicable.
A Calculation No. 11042-0205
~Revision No. 2AR EVA Calculation Page 16 of 87Requirement (e) requires that the component temperature be less than 800 °F for austenitic stainless steel.The maximum DSC shell temperature (entire shell including the lid region) is 611 °F (See Table 3).Therefore this requirement is satisfied.
Requirement (f) states that the material must have a specified yield stress to ultimate stress ratio of less than0.8. For the 61 BTH DSC which used SA-240 Type 304 steel, the ratio is 30/75 = 0.4. Therefore thisrequirement is satisfied.
Based on the discussion above (primarily the fact that cyclic conditions are not a design factor for the DSC),there is no need to consider limits on primary plus secondary stresses.
Therefore, thermal stresses are notincluded in this analysis.
Special Stress LimitsIn addition to the primary and primary plus secondary stress limits the Code [Ref. 5.7] also imposes SpecialStress Limits as detailed in paragraph NB-3227.
The applicable special stress limits are discussed below inrelation to the DSC top end cover plate welds.Bearing Loads: There are no significant bearing loads affecting the ITCP and OTCP closure welds duringService Level A, B, or C loading.
During the Service Level D side drop event, bearing stress exists at thecontact surface between the DSC and Transfer Cask. However, as noted in ASME Section III Appendix F[Ref. 5.9] paragraph F-1331.3, bearing stress need only be evaluated for pinned and/or bolted joints.Therefore this special stress limit is not applicable to this evaluation.
Pure Shear: Although the ITCP and OTCP closure welds are loaded in shear by internal pressure loading,the stress state is not pure shear due to the additional bending stresses.
Paragraph NB-3227.2 ofReference 5.7 clarifies that this stress limit is applicable to "for example, keys, shear rings, screw threads."
Therefore this special stress limit is not applicable to this evaluation.
Progressive Distortion of Nonintegral Connections:
The ITCP and OTCP closure welds are integral andtherefore not nonintegral connections.
Furthermore, there are no sources of significant cyclic loading thatwould cause progressive distortion of the DSC. Therefore this special stress limit is not applicable to thisevaluation.
Triaxial Stress: The purpose of the code [Ref. 5.7] limit on triaxial stress is to provide protection againstfailure due to uniform triaxial tension [Ref. 5.13 Chapter 4.5]. Internal pressure in the DSC and bending ofthe cover plates may cause tension in the weld in the radial and circumferential directions, but there is nosource for tension in the axial direction.
Therefore failure due to hydrostatic tension in the weld metal is notcredible.
Therefore this special stress limit is not applicable to this evaluation.
Fracture and Flaw Extension Although linear-type flaws have been identified in the structure, the critical failure mode of the welds is plasticcollapse.
Under one-time
- loading, elastic and plastic crack extension are not a concern for the very toughtype 304 stainless steel materials of the DSC shell, OTCP, and ITCP. This conclusion is supported byASME Section Xl Article C-4000 "Determination of Failure Model" [Ref. 5.10] which states that for austenitic wrought material and non-flux welds, "plastic collapse is the controlling failure mode." Note that the 61 BTHType I DSC OTCP and ITCP closure welds are made with the GTAW method [Reference 5.14] which is anon-flux type of weld.
A Calculation No. 11042-0205 A _ _ _ __io N.A R EVA Calculation Page 17of87Additionally, there is no source for fatigue flaw extension.
The only cyclic loads on the 050 are minor dailyand seasonal temperature fluctuations.
Therefore, cyclic fatigue growth of the flaws in not a crediblephenomenon.
Combined with the discussions above, the limit load analysis of the DSC top cover plates and closure weldsis sufficient to satisfy all of the applicable stress criteria of the Code [Ref. 5.7].Residual StressResidual stress due to welding is a secondary stress and therefore is not considered in the limit loadanalyses performed in this calculation, as the Section III Code [Ref. 5.7] does not require it in the limit loadanalysis.
4.2 Load CasesTable 1 lists the design basis load combinations for the 61 BTH DSC. This calculation is concerned with allload cases beginning with the inner top cover plate weld, identified as Load FL-6 in Table 1.The loading conditions of interest in this evaluation are internal and external pressure and inertial loads dueto handling,
- transfer, seismic, and accidental drop conditions.
As discussed in Section 4.1, secondary (thermal) loading is not considered.
Note that the discussions below, and the analyses performed in this calculation, are based on theconservative design values for internal pressure
- loading, rather than the actual calculated values of internalpressure.
Table 2 summarizes the conservative design values as well as the actual calculated values.Temperatures used for the material properties for each Service Level condition are listed in Table 3 anddiscussed further in the paragraphs below.Service Level AThe bounding Service Level A load combination for the DSC top end cover plates and welds is load caseTR-5 which combines the hot ambient condition with internal pressure and lg axial inertial loading.
Theother directions of inertial loading are not considered critical since their effects are not directly additive to theinternal pressure
- loading, and furthermore they are bounded by the 75g side drop load discussed furtherbelow.The lg axial load will cause the DSC payload weight (fuel, basket, holddown ring, shield plug) to bearagainst the ITCP. The total maximum payload weight is 75,811 lbs conservatively including the weights ofthe ITCP and OTCP [Ref. 5.2 Section 10.2]. The equivalent uniform pressure applied to the top-endcomponents is therefore:
Pfuellg = 75,811 =2. sP .et l =x (66.25 tn)2=2. sWhere the inner diameter of the DSC shell is 66.25 inches.
A Calculation No. 11042-0205 Revision No. 2A R EVA. Calculation Page 18of87Therefore, the bounding Service Level A case is a uniform 10 psi internal pressure (for a Type 1 DSC) plusan additional 22.0 psi acting on the shield plug in the outward axial direction of the 050 Shell.Conservatively, this analysis considers the combined 10+22=32 psi load as a uniform internal pressure in theDSC Shell. This is very conservative since the fuel pressure load which is applied to the inner surface of theshield plug would in reality be distributed to the perimeter of the ITOP as a line load by the significant stiffness of the 7-inch thick shield plug. In other words, the approach used in this calculation maximizes thebending loads on the cover plates and therefore maximizes the loading on the closure welds.Note that the cases with external pressure loading are discussed below.
A Calculation No. 11042-0205 Revision No. _ __--AR EVA Calculation Page 19of87Service Level BThe bounding Service Level B load combination for the DSC top end cover plates and welds is thecombination of the hot ambient condition with the off-normal internal pressure of 20 psi (LD-6). All of theother Service Level B conditions, such as ram push/pull loads, do not affect the top end components.
Therefore, the bounding Service Level B case is a uniform 20 psi internal pressure.
Since the pressureloading is smaller (20 psi for SL B versus 32 psi for SL A as described above), the temperature used for SL A(500 °F) bounds the maximum SL B temperatures (416 °F), and since the same limit load acceptance criterion is used for Service Levels A and B, this case is bounded by Service Level A.Service Level CThe bounding Service Level C load combination for the DSC top end cover plates and welds is one thatbounds HSM-4 and HSM-8 which combines the hot ambient condition, normal internal pressure (20 psi), andseismic loading.
- However, the seismic loads are bounded by the handling loads [Ref. 5.2 Section 7.8]discussed above for Service Level A. In addition, the acceptance criteria for Level C limit load analysis isgreater than Service Levels A and B. Therefore, all Service Level C conditions are bounded by the ServiceLevel A case described above.Note that the other Service Level C cases (such as LD-7 and UL-7) are for accident condition DSC rampush/pull loads. These loads do not affect the DSC top end components.
Therefore they are not applicable to this analysis.
Note that cases with external pressure loading are discussed below.
A Calculation No. 11042-0205 Revision No. 2AR EVA Calculation Page 20of87Service Level DThree load combinations are found to be critical for Service Level D loading of the DSC top end components, namely:* accident level internal pressure* corner drop* side dropThe first load combination is bounded by HSM-5 or HSM-6 which consist of 65 psi internal pressure due toHSM blocked vent thermal conditions.
This load is not combined with any other load that affects the top-endcomponents.
Therefore, the first bounding Service Level D load case considered in this analysis is 65 psiinternal pressure.
Note that in this condition the maximum DSC shell temperature is 611 °F and 625 0F isconservatively used in this analysis (See Table 3).The other Service Level 0 conditions consist of the drop events and accident-level seismic loading.
Theaccident seismic loads are bounded by the handling loads [Ref. 5.2 Section 7.8] discussed above for ServiceLevel A. The end-drop load is not a credible event [see footnote 12 to Table 1] but was used in the originalcalculation
[Ref. 5.2] to bound the corner drop event. However, that analysis produced negligible load in thetop cover plate welds due to the idealized boundary conditions.
As a result of an RAI by the NRC, the cornerdrop is considered using an alternate idealization that maximizes the load in the top cover plate welds. Inthis case, the 25-g corner drop load has an axial component that may be considered to load the top endcover plates with the inertia of the fuel, shield plug, hold-down ring, ITOP and OTCP. This case is evaluated in Section 6.4.The 75g side drop load TR-1 0 is considered a critical load case and is evaluated in detail. Note that this loadcase represents 75x more load than the Service Level A lg inertial loads. As discussed in Section 4.1, theService Level D acceptance criterion is only 2.1 times less stringent than the Service Level A/B criterion.
Therefore, evaluation of the 75g side drop case using the Service Level D criterion is bounding of theService Level A transverse inertial loading.
(Also, as discussed in Section 4.3, the boundary conditions usedfor the 75g side drop analysis are conservative and representative of the boundary conditions encountered for the Service Level A inertial loads and seismic loading.)
The 75g side drop case also includes the off-normal internal pressure of 20 psi, as shown in Table 1.Note that the side drop event TR-10 occurs during transfer operations which result in a maximum DSC shelltemperature of 500 °F as shown in Table 3. The higher Service Level D temperature of 625 0F discussed above occurs only during DSC storage in the HSM, and therefore is not combined with the side drop loading.
A Calculation No. 11042-0205 Revision No. _ __AR EVA Calculation Page 21lof 87External Pressure LoadingExternal pressure is present on the DSC in load cases DD-2 (vacuum drying, Service Level A) and HSM-9/10 (flood load, Service Level C). (The load cases with hydrostatic external pressure are due to the caskannulus being filled with water while the cask and DSC are in the vertical position.
In this case the pressureload varies from zero at the top of the DSC to a maximum value at the bottom of the cask. Since theexternal pressure near the cover plates is essentially zero, these cases are not critical and are notconsidered further in this calculation.)
In load case DD-2, the external pressure is 14.7 psi (full vacuum).
This pressure is bounded by the ServiceLevel B off-normal pressure (20 psi) and therefore primary stresses in the cover plates and welds arebounded by the internal pressure load cases. Stability concerns of the DSC shell are not affected by thepresence of weld flaws since they are at the end of the cask, remote from the locations at which bucklingwould occur. Additionally, the external pressure is resisted directly by the shield plug and the shield plugsupport ring, rather than by the OTCP and ITCP welds. Therefore external pressure load case DD-2 is notcritical and is not considered further in this analysis.
In load case HSM-9/1 0, the flood load is due to a 50-foot static head of water, which is equivalent to 22 psiexternal pressure
[Ref. 5.2 Section 7.9]. This pressure is bounded by the 32 psi internal pressureconsidered for Service Level A discussed above. Therefore the flood load case HSM-9/1 0 is bounded by theother internal pressure load cases.SummaryThe bounding load cases considered for the limit load collapse analyses are therefore:
(See Table 3 for temperature references)
(See Section 4.1 for explanation of the 1.5 and 1.11 factors for Service levels A/B and D, respectively, andalso for the 0.8 factor which accounts for limitations in the weld examination and inaccessible weld regions,as discussed in Section 2.0, Assumption No. 1.Service Level A/B: 32 psi Uniform Internal
- Pressure, Properties at 500 0F(Accounts for internal pressure
+ inertial load of DSC contents onto Lid)Limit load collapse pressure required to satisfy criteria:
1.5*32/0.8
= 60 psiService Level D-1: 65 psi Uniform Internal
- Pressure, Properties at 625 °FLimit load collapse pressure required to satisfy criteria:
1.11*65/0.8
= 90.2 psiService Level D-2: 75g Side Drop Acceleration plus 20 psi Uniform Internal
- Pressure, Properties at 500°F.Limit load collapse acceleration required to satisfy criteria:
1.11*75/0.8
= 104 g
~Calculation No. 11042-0205
- Revision No. 2AR EVA Calculation Page 22of874.3 FEA Model DetailsSeveral finite element models of the top half of the 61 BTH DSC are constructed in ANSYS based on theReference 5.3, 5.4, 5.5 drawings.
The models fall into two basic categories:
axisymmetric (2D) and half-sym metric (3D).The axisymmetric models use ANSYS plane element type PLANE182, a 4-node axisymmetric plane elementwith non-linear capabilities.
Each node has 2 degrees of freedom (translation in the X (radial) and Y (axial)directions).
The default element options are used in the analysis.
Sensitivity studies were performed toensure that there were no adverse effects on the results due to the potential shear locking of the elements.
(Sensitivity runs used KEYOPTION 1=3 to invoke the simplified enhanced strain formulation to relieve shearlocking.)
Additional discussion of the sensitivity analyses is provided in Section 6.0.Contact between the ITCP and OTCP is simulated using nodal coupling in the Y (axial) direction.
(SeeSection 6.1.2 for a sensitivity study using contact elements at this interface.)
No contact is defined between the opposing faces of the weld flaws. In other words, whereas compressive loading normal to the plane of the flaw may in reality be transmitted via compression through the crack facesurfaces, this load path is ignored.
This is conservative, and considered necessary since it is difficult (orimpossible) to deduce from the PAUT data what separation may exist between the two faces of the flaws.Also, no contact is considered between the DSC shell inner diameter and the ITCP and OTCP outerdiameters.
As seen in Figure 4, the fabricated dimensions of the lids and shell result in small radial gapsbetween the outer diameter of the lids and the inner surface of the shell. During the welding process, thesegaps close, but since a small remaining gap cannot be ruled out, this analysis conservatively assumes thatthe as-fabricated gap exists, as shown in Figure 4. Even if the lids deflect in the analysis such that the gapswould close, the resulting contact/compressive load path is conservatively neglected.
This is conservative since it forces all loads in the lid to travel through the weld, rather than through compression between the lidsand shell.Figure 11 and Figure 12 show images of the axisymmetric model. Loading and boundary conditions arediscussed in the following sections.
A Calculation No. 11042-0205
~Revision No. 2A REVA Calculation Page 23 of 87The 3D, half-symmetric model uses ANSYS solid element type SOLID1 85, an 8-node brick (or 6-node prism)element with non-linear capabilities.
Each node has 3 degrees of freedom (translation in the X, Y, and Zdirections).
The default element options are used in the analysis.
Sensitivity studies were performed toensure that the mesh was adequate.
Additional discussion of the sensitivity analyses is provided in Section6.0.Contact in the half-symmetry model is defined using ANSYS element types CONT173 and TARGE170.
Contact is defined between the following interfaces:
- OTCPto ITCP* ITCP to Shield Plug* Shield Plug outer diameter to DSC Shell* Shield Plug bottom surface to Support Ring* Support Ring to DSC ShellThe default contact parameters are used, although the contact stiffness is reduced in some cases to aid inconvergence.
Due to the large contact areas and since the contact areas are generally remote from thecritical stress regions, the contact stiffness is not considered a critical parameter.
The default contactparameters include:
[Reference 5.6]* Penetration tolerance factor: Default value = 0.1. This parameter controls the acceptable level ofpenetration of the contact node into the target surface, based on the depth of the element underlying the target element.* Pinball region scale factor: Default Value = 1.0. This parameter controls the extents of the regionaround each contact node that is checked for contact with target segments.
The default volume is asphere of radius 4*depth of the underlying element.* KeyOption 2: Contact algorithm:
Default = Augmented Lagrangian.
The contact method is aniterative penalty method where the contact pressure is augmented during the equilibrium iterations sothat the final penetration is within the acceptable tolerance.
- KeyOption 4: Location of contact detection point: Default = On Gauss Point. Other options includeusing the nodal points, normal to either the contact surface or the target surface.
The default optionis suggested for general cases.Other features and controls of the CONTA1 73 elements are related to advanced features (bonded contact,cohesion, etc.) and initial penetration and gap controls which are not utilized in this analysis.
Figure 19 through Figure 21 show images of the half-symmetry model. Loading and boundary conditions arediscussed in the following sections.
Table 6 shows a summary of the ANSYS models and analyses which are performed.
Further details on thevarious ANSYS models are provided below.
A Calculation No. 11042-0205
~Revision No. 2AR EVA Calculation Page 24of874.3.1 Axisymmetric Case #1The first case considered is a combination of OTCP Flaw Set #1 and the ITCP bounding flaw set discussed in Sections 3.2.1.1 and 3.2.2, respectively.
The mesh and flaw details for this case, called Axisymmetric Case #1, are shown in Figure 13 and Figure 14.The mesh shown in these figures was created based on a basic goal of having at least 4 elements acrossthe thickness of the net sections of the weld, as reduced by the flaws. In order to investigate the effects ofmesh density, a refined mesh was created for this case, as shown in Figure 15. Since the sensitivity modelshown in Figure 15 only refined the weld region an additional model was created as shown in Figure 16 toensure a sufficient mesh in the lid interior region.This model, and all of the other axisymmetric models discussed below, are used for analysis of uniforminternal pressure loading.
The model is constrained in the radial direction at the axis of symmetry and in theaxial direction at the bottom cut of the DSC shell near the mid-length of the cask (remote from the top endcomponents of interest.)
The pressure loading is applied to the internal pressure boundary (bottom surfaceof ITCP, surface of ITCP weld to Shell, and Shell inner surface).
(See Section 6.1.2 for a sensitivity analysiswhere internal pressure is included on the ITCP weld root flaw internal surfaces.)
4.3.2 Axisymmetric Case #2The second case considered is a combination of OTCP Flaw Set #2 and the ITCP bounding flaw setdiscussed in Sections 3.2.1.2 and 3.2.2, respectively.
The mesh and flaw details for this case, calledAxisymmetric Case #2, are shown in Figure 17 and Figure 18 for the refined mesh. Based on the results ofthe Axisymmetric Case #1 (See Section 6.1.2), the initial mesh level described above for Case #1 would besufficient.
- However, since the run times remained reasonable, a refined mesh model (weld and lid interiorregions) was generated and is used for Case #2.4.3.3 Axisymmetric Case #0In order to study the effect of the flaws, a 3rd case is considered in which the flaws are removed and the as-designed collapse load is determined.
Only the refined mesh model (weld and lid interior regions) isconsidered.
The mesh is identical to Figure 16 but the coincident nodes along the crack faces are merged.
A Calculation No. 11042-0205
.... .Revision No. 2A R EVA Calculation Page 25 of 874.3.4 Half-Symmetry (3D) Case #1The 3D model is based on the Axisymmetric Case #1. (Analysis results showed that there was negligible difference in the results from Axisymmetric Case #1 and Case #2. The total projected cross-sectional areaof the flaws in Case #1 is greater than Case #2. Therefore, Case #1 is considered critical for the side droploading).
The same flaw pattern is modeled, but the initial mesh is slightly less refined in order to obtain reasonable run times. Mesh sensitivity studies are described below. The half-symmetry model is used for internalpressure loading (as a benchmark case to study the effects of mesh refinement) and also for side-drop loading.The shield plug support ring is connected to the DSC shell at the two corners using nodal DOE couples torepresent the fillet welds used to join the two parts.In order to improve the numerical stability of the ANSYS model, soft springs (COMBIN 14) elements are usedto connect the shield plug to the support ring. The springs have a stiffness of 1 lb/in. The low stiffness combined with the very small relative deflections between these parts results in negligible internal force inthe springs.
The forces in the springs at the final converged solution are reviewed to confirm that the springforces are small.In all load cases, symmetry conditions are applied to the cut face of the model. Axial constraints are appliedat the bottom cut of the DSC shell near the mid-length of the cask (remote from the top end components ofinterest.)
For the internal pressure load case, the model was further reduced to a 90-degree model andsymmetry constraints were placed on both cut faces of the model.The purpose of this calculation is to evaluate the effects of the closure weld flaws and qualify the welds andany other components affected by the welds. All other aspects of the DSC (such as the shell remote fromthe welds) are not in the scope of this calculation.
The modeling approach (loads and boundary conditions) for the side drop event are considered in light of this purpose and are described in the following paragraphs.
For the side drop cases, the OD of the canister shell is constrained in the vertical (drop) direction for a smallsector (approximately 1.5,, inches or 2.8 degrees) of assumed contact.
In reality the DSC is supported insidethe Transfer Cask (TC) during this event. Therefore the true boundary condition would either be a line ofcontact along a TC rail (which is 3" wide) or a line of contact at areas remote from the rails. As deformations
- increase, the area of contact would also increase.
As discussed below in Section 4.4, deflections are over-estimated in a limit load analysis.
Therefore, the area of contact with the TC rail or inner surface is assumedto be constant.
This conservatively neglects the increase in contact area that would occur during the dropdeformations.
Additionally, this boundary condition is representative of the DSC storage condition inside theHSM, where the DSC rests on the 3-inch wide steel rails.As discussed in the Reference 5.2 calculation, the DSC payload (basket and fuel) are located approximately 21.5 inches away from the ITCP and are therefore considered to have no effect on the DSC lid components.
The effect of the basket and fuel loading on the DSC shell is considered in the basket design-basis calculation for side-drop loading.
The basket hold-down ring is.a grid-type structure that does not represent significant weight and is of sufficient strength and stiffness to be self-supporting during the side drop and notsignificantly affect the DSC shell and adjacent regions.
Therefore, as in the Reference 5.2 calculation, theDSC payload is not considered as affecting the top-end components and the weight is applied as a pressurealong a strip of elements at the impact region, beginning approximately 23" below the ITCP. Since the loads A Calculation No. 11042-0205
- ....
No. 2-AR EVA Calculation Page 26 of 87are essentially applied directly over the supported (impact) region of the DSC shell, they have no appreciable effect on the shell deformations.
Images of the Half-Symmetry model are shown in Figure 19 to Figure 21.In order to study the adequacy of the mesh for the half-symmetry model, an internal pressure load case wasperformed and compared to the results of the axisym metric case refined mesh. This study confirms theadequacy of the mesh in the cross-section of the 3D) model. In order to evaluate the mesh in thecircumferential direction, a model was created with a refined mesh in the regions of the model showing largeplastic strains (the impact region) and locations where tensile stress is expected in the weld (at the 90-degree location where the lid resists ovalization of the DSC shell). This model is shown Figure 22.
A Calculation No. 11042-0205
~Revision No. 2AR EVA ~~CalculationPae2o8 4.3.5 Half-Symmetry (3D) Case #0In order to study the effect of the flaws, an additional case is considered in which the flaws are removed fromthe model and the as-designed side drop limit load capacity is determined.
4.4 Limit Load Solution DetailsAs discussed in Section 4.1, this calculation is based on predicting the lower-bound collapse loads of theDSC based on limit load analysis.
All materials are modeled as elastic-perfectly plastic1, with yield stressvalues based on the limit load analysis requirements of the ASME code [Ref. 5.7]. Table 3 lists thetemperatures used for each load case, and the values of the material properties are shown in Table 4 andTable 5.The prescribed loads are applied to the model, and then are increased linearly until the solution fails toconverge.
The analyses use small deflection theory (NLGEOM,OFF).
This is conservative since deflections areunrealistically high in a limit load analysis due to the lower-bound non-strain-hardening material properties that are used. If large deflections were to be considered, the beneficial effects of OTCP and ITCPmembrane action and of increased contact areas would be over-estimated, resulting in non-conservative effects.
This was verified with a sensitivity study using NLGEOM,ON, which resulted in much higher collapsepressures.
This confirmed that using NLGEOM,OFF is appropriate, and conservative.
In addition, Paragraph 5.2.3.1 of Reference 5.18 states that small displacement theory is to be used in a limitload analysis.
I "Elastic-perfectly plastic is standard mechanics of materials term that describes an idealized material that behaves in alinear-elastic manner up to the yield point, and thereafter is perfectly-plastic, i.e. non-strain hardening.
A Calculation No. 11042-0205
- Revision No. 2AR EVA Calculation Page 28of8
75.0 REFERENCES
5.1. AREVA Document No. 180-9236022-000.
NDE Services Final Report. "Monticello, DSC-16,Phased Array UT Examination Results of the Inner and Outer Top Cover Lid Welds."Revision 0.5.2. AREVA (Transnuclear)
Calculation No. NUH61BTH-0200 Revision
- 0. "NUHOMS-61BTH Type 1 Dry Shielded Canister Shell Assembly Structural Analysis."
5.3. ARE VA (Transnuclear)
Drawing No. NUH61 BTH-3000 Revision
- 8. "NUHOMS 61 BTH Type 1DSC Main Assembly."
5.4. AREVA (Transnuclear)
Drawing No. NUH61 BTH-3001 Revision
- 4. "NUHOMS 61 BTH Type 1DSC Shell Assembly."
5.5. AREVA (Transnuclear)
Drawing No. NUH61 BTH-4008 Revision
5.6. ANSYS Version 14.0. ANSYS Inc. (Including the ANSYS Mechanical APDL Documentation).
5.7. ASME Boiler and Pressure Vessel Code,Section III Subsection NB. 1998 Edition withAddenda through 2000.5.8. ARE VA (Transnuclear)
Document Number NUH-003 Revision
- 14. "Updated Final SafetyAnalysis Report for the Standardized NUHOMS Horizontal Modular Storage System forIrradiated Nuclear Fuel."5.9. ASME Boiler and Pressure Vessel Code,Section III Appendices.
1998 Edition with Addendathrough 2000.5.10. ASME Boiler and Pressure Vessel Code,Section XI. Rules for Inservice Inspection of NuclearPower Plant Components.
1998 Edition with Addenda through 2000.5.11. AREVA (Transnuclear)
Document No. NUH61BTH1-0101 Revision
- 2. "Design CriteriaSpecification for the NUHOMS-61BTH Transportable Storage Canister."
5.12. AREVA Calculation No. 11042-0204 Revision
- 3. "Allowable Flaw Size Evaluation in the InnerTop Cover Plate Closure Weld for DSC #16"5.13. Chattopadhyay, Somnath.
"Pressure Vessels Design and Practice."
CRC Press. 2004.5.14. TriVis Incorporated Welding Procedure Specification No. SS-8-M-TN Revision 10.5.15. ASME Boiler and Pressure Vessel Code,Section II, Part C. "Specifications for Welding Rods,Electrodes, and Filler Metals."
1998 Edition with Addenda through 2000.5.16. ASME Boiler and Pressure Vessel Code,Section II, Part D. "Properties."
1998 Edition withAddenda through 2000.5.17. AREVA (Transnuclear)
Calculation No. NUH61BTH-0253 Revision
- 0. "NUHOMS 6IBTHType I DSC Shell Assembly Outer Top Cover Plate Critical Flaw Size of Weld."5.18. ASME Boiler and Pressure Vessel Code,Section VIII Division
- 2. 2010.5.19. AREVA Document No. 54-PQ-1 14-001 Revision
- 0. "Phased Array Ultrasonic Examination ofDry Storage Canister Lid Welds."5.20. NRC Spent Fuel Project Office Interim Staff Guidance
-15. Materials Evaluation.
(ISG-15).
1110/2001.
A Calculation No. 11042-0205
.Revision No. 2AR EVA Calculation Page 29 of 876.0 ANALYSISTable 6 shows a summary of the results of all of the analyses performed for this calculation and includes acomparison of the results with the acceptance criteria.
Each case is discussed in more detail below.6.1 Axisymmetric Analyses for Internal Pressure6.1.1 Axisymmetric Case #1 -Initial Mesh ModelTwo analyses are performed with the Axisymmetric Case #1 initial-mesh model described in Section 4.3.1:one case using the Service Level A/B material properties and one case using the Service Level D materialproperties.
The collapse pressures were determined to be 95.9 psi for Service Level A/B and 136.6 psi forService Level 0. Figure 23 shows various plots of the plastic strain in the initial-mesh model for ServiceLevel A/B at various locations and levels of loading.
These strain plots are also representative of thebehavior of the Service Level D analysis.
Figure 24 shows the deflection history at the center of the lid, andindicates the expected plastic instability that occurs as the limit load is approached.
Note that both thestrains and displacements presented in these figures have no physical meaning and the displacement plotsshow only the loading (pressure) at which the solution fails to converge.
Since the initial mesh contains several element divisions at each critical cross-section, it is not expected thatelement shear locking (due to the default fully-integrated elements) will be significant.
To confirm this, a testcase was done using the Service Level A/B model but with the Simplified Enhanced Strain elementformulation (KEYOP 1=3). The collapse pressure was found to be 96.1 psi, which is essentially identical tothe initial results.6.1.2 Axisymmetric Case #1 -Refined Mesh ModelsAdditional analyses are performed using the Service Level A/B material properties with the refined meshmodels described in Section 4.3.1. Figure 25 and Figure 26 show the plastic strain results for the refinedmesh at the weld region and the refined mesh at the weld and lid interior
- regions, respectively.
The collapsepressures were found to be 94.8 psi and 93.8 psi, respectively, for these models. The OTCP deflection histories are shown in Figure 27. Note that both the strains and displacements presented in these figureshave no physical meaning and the displacement plots show only the loading (pressure) at which the solutionfails to converge.
Figure 28 shows a comparison of the maximum displacement history curves for the various Axisymmetric Case #1 models, done as part of the mesh sensitivity study. As seen in the figure, the results match verywell. The results of the refined mesh models deviate at most (95.9-93.8)/93.8
= 2.2% from the initial meshresults.
This is very close agreement particularly due to the non-linear nature of the analysis.
Therefore, theinitial mesh is considered sufficient.
However since the analysis run times for the axisymmetric cases arereasonable even for the refined mesh model, the remaining axisymmetric cases use a refined mesh.The Axisymmetric Case #1 with refined weld and lids for Service Level D criteria reported a collapsepressure of 132.6 psi.Note that the nodal coupling in the axial direction between the ITCP and OTCP is a valid method to modelthe contact between the plates since the internal pressure loading ensures that the ITCP lid will bear againstthe OTCP, and since the nodes that are coupled remain coincident throughout the analysis, with only very A Calculation No. 11042-0205
, Revision No. 2AR EVA Calculation Page 30of87minor differences in radial position occurring at the later load steps. In order to confirm the behavior of thenodal coupling, the Axisymmetric Case #1 model with refined welds and lids was modified to include contactbetween the ITCP and OTCP. The model replaces the nodal coupling with CONTA171 and TARGE169elements, using the default element parameters.
Figure 29 shows a comparison between the model usingDOE couples and the model using contact elements.
As seen in the figure, the results are very similar, withthe DOF-couple-model showing slightly more conservative results.
Therefore, the nodal coupling isacceptable and is used in all other axisymmetric models.Note that in all of the FEA models, the internal pressure loading was not applied to the faces of the ITCPweld root flaw that is exposed to the internal region of the cask. Pressure loading on this crack face isnegligible since the flaw is only 0.09" high, and in reality the ITOP flaws are generally very short (i.e. not full-circumferential flaws). In order to support this conclusion, a sensitivity analysis is performed where thepressure loading is applied to the ITCP weld root crack faces. The results, shown in Figure 30, confirm thatpressure loading on the faces of this flaw are negligible.
A Calculation No. 11042-0205 Revision No. 2AR EVA Calculation Page 3l of 876.1.3 Axisymmetric Case #2Two analyses are performed with the Axisymmetric Case #2 refined-mesh model described in Section 4.3.2:one case using the Service Level A/B material properties and one case using the Service Level D materialproperties.
The collapse pressures were determined to be 93.7 psi for Service Level A/B and 132.9 psi forService Level 0. Figure 31 shows various plots of the plastic strain for Service Level A/B at various locations and levels of loading.
These strain plots are also representative of the behavior of the Service Level Danalysis.
Note that both the strains and displacements presented in these figures have no physical meaningand the displacement plots show only the loading (pressure) at which the solution fails to converge.
6.1.4 Axisymmetric Case #0One analysis is performed with the Axisymmetric Case #0 refined-mesh model described in Section 4.3.3using the Service Level A/B material properties.
The collapse pressures were determined to be 94.5 psi forService Level A/B. Figure 32 shows various plots of the plastic strain at various locations and levels ofloading.
Note that both the strains and displacements presented in these figures have no physical meaningand the displacement plots show only the loading (pressure) at which the solution fails to converge.
Figure 33 shows a comparison of the maximum center-of-lid displacement history for all three axisymmetric cases. As seen in the figure, there is essential no difference between Axisymmetric Case #0, Case #1 andCase #2. The Case #1 and Case #2 analyses show slightly larger deflections early in the analysis due to theslightly reduced rotational fixity of the welds. However, the final collapse pressure are within (94.5-93.7)/93.7=0.9%
of each other. This supports a supposition that the observed flaws have negligible impacton the governing failure mode of the top end closure plates and welds. Again, displacements from the limitload analyses have no physical
- meaning, other than to show the onset of non-convergence of the FE model.
A Calculation No. 11042-0205
...Revision No. _ __A R EVA Calculation Page 32 of876.2 Half Symmetry Analyses for Internal Pressure (Benchmark Cases)The model described in Section 4.3.4 is used for an internal pressure collapse analysis in order tobenchmark the model against the axisymmetric cases. The collapse pressure was calculated to beapproximately 97 psi. (The run was terminated at 95 psi and the final collapse pressure was estimated toavoid excessive computer run time). Figure 34 shows various plots of the plastic strain at various locations and levels of loading.
A comparison of the half-symmetry case to the refined-mesh axisymmetric case isshown in Figure 35. As seen in the figure, the half-symmetry case closely matches the behavior of therefined mesh axisymmetric model although the results indicate a slightly greater collapse pressure.
Therefore, the half-symmetry model is considered sufficiently accurate for this analysis.
As shown by theresults, and as discussed in Section 7.0, there is significant safety margin available such that further meshrefinement of the half-symmetry model is not warranted.
- However, the effects of circumferential meshdensity for the half-symmetry model can be seen in Section 6.3.1.
A Calculation No. 11042-0205
.... ,Revision No. 2AREVA Calculation Page 33 of 876.3 Half Symmetry Analyses for Side Drop Loading6.3.1 Half-Symmetry Case #1The model described in Section 4.3.4 is used to perform two side-drop limit load analysis.
One caseincludes side-drop acceleration loading only, while the second case includes the DSC off-normal internalpressure of 20 psi. For this later case, the 20 psi internal pressure is applied simultaneously with a 75gacceleration, and then both the pressure and the acceleration are increased linearly until the collapse g-loadis obtained.
(For example, for collapse occurring at 181g, the internal pressure at collapse is 20*181/75
=48.3 psi.)Note that the side drop loading is combined with the design-basis off-normal internal pressure of 20 psi, asopposed to the internal pressure value of 32 psi used for the SL A/B cases which was the sum of the 10 psinormal pressure and an additional 22 psi to account for inertial handling/seismic loads. See Section 4.2The collapse g-load for side-drop-only loading was found to be approximately 181g. The collapse g-loadwhen internal pressure loading was included was found to be greater than 181g. This later run terminated at181g, but based on the collapse behavior (see Figure 41) it is expected that smaller time steps would allowthe solution to continue to larger loads.Various images of the stress and strain in the side drop analyses are shown in Figure 36 to Figure 38. Notethat both the strains and displacements presented in these figures have no physical meaning and thedisplacement plots show only the loading (pressure) at which the solution fails to converge.
The Half-Symmetry Case #1 model with refined mesh in the circumferential direction was used to evaluatethe side drop load case (without internal pressure).
This analysis was performed up until a load of 185 g's,at which time the analysis was terminated manually to avoid large file sizes and excessive run time. As seenin Figure 41, this model showed a greater resistance to the side-drop
- loading, and would eventually result incollapse g-loads in excess of 185 g if smaller timesteps and longer run times were provided.
Images of thestress and strain from this analysis are shown in Figure 39. Note that both the strains and displacements presented in these figures have no physical meaning and the displacement plots show only the loading(pressure) at which the solution fails to converge.
This analysis confirms that the mesh used in the other half-symmetry cases is adequate, and conservative.
A Calculation No. 11042-0205 Revision No. 2AR EVA Calculation Page 34of876.3.2 Half-Symmetry Case #0One side drop analysis is performed with the Half-Symmetry Case #0 model (no flaws) described in Section4.3.5. Based on the results discussed above, only the case without internal pressure loading wasconsidered.
This analysis resulted in a collapse load of 189g. Stress and strain plots from this analysis areshown in Figure 40. As noted previously, these strain plots have no physical meaning.
As shown in Figure41, the collapse behavior was nearly identical to the case with weld flaws, indicating that the flaws hadnegligible effect on the results.6.4 Evaluation of the 25g Corner DropReference 5.2 Section 10.2 evaluated the OTOP weld to resist a 25g inertial load on the entire DSC contentsand neglecting the strength of the ITCP weld. Furthermore, a conservative stress was assumed in the welddue to internal pressure.
The Reference 5.2 calculation is revised below to account for the strength of bothwelds and include a reduction in the weld thickness due to the observed flaws. The total weld thickness istaken as the combined weld throats from the ITCP and OTCP minus the height of the flaws present in thewelds. (See Reference 5.2 Section 10.2 for the basis of the following values and calculations.)
Note that the allowable weld stress noted below includes a joint efficiency factor of 0.7 as described inReference 5.2 Section 6.2. This reduction factor conservatively bounds the reduction factor of 0.8 discussed in Section 3.4. Therefore, no further reduction factor is applied, and the calculation below is conservative.
WTor = 75,81l1bs (total weight of fuel + basket + lids and shield plug)Wp= 68,943 lbs (load due to axial component of pressure)
WTOT2Sg = 25 x 75,811 + 68,943 = 1,964,218 lbsWTOT2EO -1,964,218 947L-b (ethof weld is20.3"W2=- Lweld 208.131 in UOIL3 1tweld + -(0.23 + 0.11) = 0,3475 in (, See Note)wT~g = --- 9,.3 = 27,157 psi (weld shear stress due to 25g corner drop) 0.347'5-T2o0s = 4,120 psi (weld stress radial component due to 20 psi internal pressure)
TTOT -= T2Sg_ +¢ T2opsi -= 27,157 + 4,120 =31,277 psiTALLOW = 32,400 psiTTOT _31,277TALW- = 0.97<1 (0K)T~tw 32,400*Note: the reduction of the weld to account for the flaws is based on the maximum flaw heights in any oneplane through each of the welds. This is taken as 0.23" for the OTCP weld and 0.11" for the IT-CP weld.)Therefore, the top end closure welds, with the observed flaws, are OK for the Service Level 0 corner dropevent.
A Calculation No. 11042-0205 Revision No. 2AR EV/A Calculation Page 35 of 877.0 DISCUSSION AND CONCLUSIONS The lower bound collapse pressure for Service Level A/B criteria was found to be 93.7 psi which is greaterthan the required pressure of 1 .5x32/0.8=60 psi (where 1.5 is the code-required
[Ref. 5.7] factor on the 32psi design pressure loading and 0.8 is the weld strength reduction factor- see Section 4.2). Therefore theService Level A/B criteria is satisfied.
The lower bound collapse pressure for Service Level D criteria was found to be 132.6 psi which is greaterthan the required pressure of 1.1 1x65/0.8=90.2 psi (where 1.11 is the code-required
[Ref. 5.7] factor on the65 psi design pressure loading and 0.8 is the weld strength reduction factor -see Section 4.2). Therefore the Service Level D criteria for internal pressure is satisfied.
As noted in Section 6.1.4 and as shown in Figure 33, there is essentially no difference in the collapsepressure and extremely little difference in the overall collapse behavior of the DSC subjected to internalpressure loading with and without flaws in the weld.The lower bound collapse acceleration for side drop (Service Level D) loading was found to be 181g which isgreater than the required load of 1.1 1x75/0.8=104g.
Therefore the Service Level D criteria for side droploading is satisfied.
As noted in Section 6.3.2 and as shown in Figure 41, there is essentially no difference in the collapse loadand behavior between the as-designed DSC and the DSC with closure weld flaws.The Reference 5.12 and 5.17 calculations document the ITOP and OTOP closure weld critical flaw sizes,respectively, based on the maximum radial stresses in the welds. The guidance and safety factors ofReference 5.10 are used in the critical flaw size analysis.
The critical flaw sizes are determined to be 0.19and 0.29 inches for surface and subsurface flaws, respectively, in the OTCP weld and 0.15 inches forsurface and subsurface flaws in the ITOP weld. The largest single OTOP flaw size documented in Reference 5.1 is 0.14 inches. As discussed in Section 3.2 a very conservative maximum combined flaw height of 0.195inches is assumed in this analysis.
The largest single ITCP flaw size documented in Reference 5.1 is 0.11inches. Therefore, the observed flaws actually are smaller than the critical flaw size limits and therefore it isnot surprising that the flaws are shown to have little effect on the capacity of the structure.
This analysisshows that the quantity and close proximity of some of the flaws also has no significant adverse effects onthe structural capacity of the OSO.Even though all observed flaws in the ITICP and OTCP welds are included in the analysis models usingconservative representations, an additional weld strength reduction factor of 0.8 is considered by increasing the limit load acceptance criteria by a factor of 1/0.8=1.25 times. The 0.8 factor, which is the samemagnitude reduction factor as in ISG-1 5 [Ref. 5.20], conservatively accounts for any additional limitations inthe efficacy of the PAUT examinations and also accounts for the inaccessible area around the vent andsiphon block as well as the geometric reflectors at the root and near the toe of the weld.Therefore it is concluded that Monticello DSC-16, remains in compliance with the ASME Section IIISubsection NB [Ref. 5.7] stress limits with the presence of the ITCP and OTCP closure weld flaws asdocumented in Reference 5.1.
A Calculation No. 11042-0205 A _ _ _ __n N.A REVA Calculation Page 36 of878.0 LISTING OF COMPUTER FILESAnalyses performed on Computer HEA-0213A using ANSYS Version 14.0 [Ref. 5.6].File Date & Time listing is as displayed by the Windows 7 Operating System -Differences may occur duelocal time zone and daylight savings settings.
Analysis Case File Name Date & TimeA~xisymmetric 1 61BTHWeldFlaw 1FAX_2_DETACH.db 4/7/2015 10:45AMInitial Mesh 618THWeldFlaw 1F AX_2_DETACH.rst 4/7/2015 10:45 AMInternal Pressure 61BTH_WeldFlaw 1F AX_2_DETACH.mntr 4/7/2015 10:45 AMSL A/B SOLUTION_AXISYMMIPLimitLoad.INP 4/7/2015 10:20 AMAxisymmetric 1 6IBTHWeldFlaw1F AX 2.DETACH.db 4/7/215 11:59 AMRefined Weld Mesh 61 BTH_WeldFlaw1 1FAX_2_DETACH.rst 4/7/215 11:58 AMInternal Pressure 61BTH_WeldFlaw 1 FAX 2_DETACH.mntr 4/7/215 11:59 AMSL A/B SOLUTION_AXISYMMIPLimitLoad.INP 4/7/2015 11:55 AMAxisymmetric 1 61BTHWeld~law IF AX_2_DETACH.db 4/21/2015 9:04 AMRefined Weld and Lid Mesh 61BTHWeld~law IF AX_2_DETACH.rst 4/21/2015 9:03 AMInternal Pressure 61BTH_Weld~lawIFAX_2_DETACH.mntr 4/21/2015 9:04 AMSL A/B SOLUTION_AXISYMMIPLimitLoad.INP 4/7/2015 11:55 AMAxisymmetric 1 618TH_Weld~law 1FAX_2_DETACH.db 4/20/2015 10:43 AMInitial Mesh 618THWeld~law 1FAX_2_DETACH.rst 4/20/2015 11:09 AMInternal Pressure 618THWeld~lawIFAX_2_DETACH.rnntr 4/20/2015 11:09 AMSL D SOLUTION_AXISYMMIPLimitLoad_SLD.INP 4/7/2015 11:20AMAxisymmetric 1 618TH_WeldFlaw 1FAX_2_DETACH.db 4/30/2015 8:12 AMRefined Weld and Lid Mesh 618THWeldFlaw iF AX_2_DETACH.rst 4/30/2015 8:12 AMInternal Pressure 618TH_WeldFlaw iF AX 2_DETACH.mntr 4/30/2015 8:12 AMSL D SOLUTION_AXISYMM IP LimitLoad_SLD.INP 4/7/2015 12:02 PMAxisymmetric 2 618TH_WeldFlaw_2GAX_2.db 4/21/2015 3:06 PMRefined Weld and Lid Mesh 618THWeldFlaw_2GAX_2.rst 4/21/2015 2:58 PMInternal Pressure 618TH_WeldFlaw_2GAX_2.mntr 4/21/2015 3:06 PMSL A/B SOLUTION_AXISYMMIPLimitLoad.INP 4/7/2015 11:55 AMAxisymmetric 2 618TH_WeldFlaw_2GAX_2.db 4/21/2015 3:10 PMRefined Weld and Lid Mesh 618TH_Weld~law_2GAX_2.rst 4/21/2015 3:10 PMInternal Pressure 61BTH_Weld~law_2GAX_2.mntr 4/21/2015 3:10 PMSL D SOLUTION_AXISYMMIPLimitLoad_SLD.INP P 4/7/2015 12:02 PMAxisymmetric 0 618THWeldFlaw 1FAX_2_DETACH.db 4/21/2015 10:39 AMRefined Weld and Lid Mesh 618THWeldFlawIFAX_2_DETACH.rst 4/21/2015 10:31 AMInternal Pressure 618TH_WeldFlaw IFAX_2_DETACH.mntr 4/21/2015 10:39 AMSL A/B SOLUTION_AXISYMMIPLimitLoad.INP 4/15/2015 11:07 AMIntilAxsmmtic1618TH_WeldFlawIFAX_2_DETACH.db 4/17/2015 5:40 PMIiilMesh with Keyoption 1=3 618TH_Weld~lawIFAX_2_DETACH.rst 4/17/2015 5:40 PMInternal Pressure 618THWeldFlaw iF AX 2 DETACH.mntr 4/17/2015 5:40 PMSLA/B SOLUTION_AXISYMM_IP_LimitLoad.INP 4/16/2015 12:27 PM A Calculation No. 11042-0205 Revision No. 2AR EVA Calculation Page 37 of87Analysis Case File Name Date & TimeAxisymmetric 1 61BTHWeld~law_1FAX_2_DETACH.db 5/19/2015 8:45 AMRefined Weld and Lid Mesh 6IBTHWeldFlaw IF AX 2 DETACH.rst 5/19/2015 8:02 AMInternal
- Pressure, SL N/BITCP/OTCP couples replaced with 61BTH_WeldFlaw_1FAX_2_DETACH.mntr 5/19/2015 8:45 AMContact SOLUTION_AXISYMM_
IPLimitLoad.INP 5/18/2015 5:03 PMAxisymmetric 1 61BTH_WeldFlaw 1F AX_2_DETACH.db 5/18/2015 2:42 PMRefined Weld and Lid Mesh 61BTH_WeldFlaw iF AX_2_DETACH.rst 5/18/2015 1:37 PMInternal
- Pressure, SL NBWith Pressure on ITCP Weld Root 61BTH_WeldFlaw 1F AX 2 DETACH.mntr 5/18/2015 2:42 PMFlaw Surfaces SOLUTION_AXISYMMIPLimitLoad.INP 5/18/2015 1:25 PMHalf Symmetry 1 61 BTH_WeldFlaw_1 GC.db 4/29/2015 2:10 PMInitial Mesh 61BTH_WeldFlawI GC.rst 4/29/2015 4:52 PMInternal Pressure 61BTH_WeldFlaw_I1GC.mntr 4/29/2015 4:52 PMSL NB SOLUTIONHALFSYMLimitLoad.INP 4/29/2015 2:10 PMHalf Symmetry 1 61BTHWeldFlawI GC.db 4/30/2015 8:20 AMInitial Mesh 61BTHWeldFlaw_
1GC.rst 4/30/2015 3:35 PMSide Drop 61BTHWeldFlaw_
1GC.mntr 4/30/2015 3:35 PMSL D SOLUTIONHALESYMSD.INP 4/30/2015 8:21 AMHalf Symmetry 1 61 BTH_WeldFlaw_I GC.db 5/1/2015 6:58 PMInitial Mesh 6IBTH_WeldFlawl 1GC.rst 5/1/2015 4:29 PMSide Drop + Internal Pressure 61BTHWeldFlaw_1GC.mntr 5/1/2015 4:13 PMSL D SOLUTION_HALFSYMSD.INP 4/30/2015 10:26 PMHalf Symmetry 1 61BTHWeldFlaw_I GD_Refined.db 5/6/2015 1:54 PMRefined Circumferential Mesh 61 BTH_WeldFlaw_1 GDRefined.rst 5/6/2015 11:48 AMSide Drop 61 BTH_WeldFlaw 1 GD_Refined.mntr 5/6/2015 11:47 AMSL D SOLUTIONHALFSYMSD.INP 5/5/2015 9:01 PMHalf Symmetry 0 61BTH_WeldFlaw_
1GC.db 5/2/2015 6:53 AMInitial Mesh 61 BTH_WeldFlaw_I GC.rst 5/2/2015 6:53 AMSide Drop 61BTH_Weld Flawl IGC.mntr 5/2/2015 3:37 AMSL D SOLUTION_HALESYMSD.INP 4/30/2015 8:21 AM AAR EVACalculation No. 11042-0205 Calculation Revision No. 2Page 38 of 879.0 TABLES AND FIGURESTable 1 -Summary of Design Basis Load Combinations for the 61 BTH DSC [Ref. 5.8]LdCaeHorizontalfDW Vertical DW Internal
.ExternaL Theanal Lit~ing Ot0e Sea-sireLdCaeDSC
.Fuel DSC Fud!
Presasur Condition
- Loads .oad 'LevelNon-Operational Load Cases0-l~ Fab. Leak Ta .tie_ .....- 14.7piog 70"F -.I55lripax~al Test.FabI" e tTne ....... -1512ps i '~ -* :l 70'F -155 !alp ax ial TestNO-S DSC Upsielatin
-... ..- 701E x -ANO-4DSC Vertica!Lift
--x -.. ... 7O5* -A:FL-i DSC$'Css Filling -. -:Cask --Hydrostatic 12(0F.Caak x x AFL-2 DSC/Cask Fillin .. -Cask -~
Hydrosttic Carl: x x A'FL-S
- -. -Cask --, 'Hyd,,.tatic Hyrotaic IOF'0I Cask -A-FL.-S Xftexs o Decon --- Cask- x Hydrostatic Hydrostatic 12OAP Cask --AFL-6
--Cask x' -Hvdro~tiati i2O0FCaal
--.A"FL-7? Fuel Deck .Semisiroadinr
--. -Cask x. "H-drostatic Hydrosta.cL 2OA l":.Cask
--Note 10 CDanizagrningtoad C.'ises,Hydrostatac DlD-1 DSC Blowdosas
--Cask x lO!15 pdig Hydrosfatic 12OAF Ca-,sk -. A÷D2V1UtD4ng- &~a~ I~~ --ADD-3 HeliurnBackfill -Cask x: 12pss Hyd-ostatic 12(O"F Cask -ADD1-4 Final HelismsBackftlU
--Cask x" ,. pai'" .Hydr'ostatic 120"1= Cask -- -- A"DD2-50OssleCover Plate Weld -~ -Cask x .3. Hydrzositic i2OAF Cask --- ATransfer Traier Loa,'ing2L-1 VerticalXfer to Traier -'Cask x !lO5psaig OA(F Cak --AT2-2 Vertical.eftoTrailer
--: -'Cask: ," l0,bps Ii-- .12O-AFCa-k
-_ .- ATL-S Laydoros Cask X " --
-O"FCasI --- ATI_-4 LnSidown Cask X --10/15!p -: I20F Cask" --ALo)ad Case .Horizental DW' "Vertacl.DW
,Internl
,External Thermal LIafdang Other Servite......___________
DSC- Fuel *DSC F~uel Pressm'e Condanano LoadA -Loads !Tr.ansfer Tofl~rom ISFScITR-4 AxialLo~ad-Cold Cask X --10/lvS~ia
-OAF 1 A~xial .-. ATR-2 Trnaue-se Load- Cold Cask X -- -- l0!15 paig -OAF 'ig Transverse
-A*.TR-3 ,Verticai Lo -Cold -_ask X --10,15 psi. "- OAF i! gVertical
-"ATR-4Oblique Lotd-'Cold Cask X -- -I0,q5 vsig -- O'F %n£Axial
-A4**TR-5 AxialLout-H~ot
.Cask .X ---J -l0.!l~psi
-lO' "!0F lg-Axial
--- A:TR.-6 Transvrse Io ad -Hot Cask: X --si -. 101rF .lg Traa. -- ATR-7 Vertical Load Cask X -10/1 paig -- g Verical -A-iR-8 Oblique Load- Hot Cask X. -L0/15 pag -. 10OAF Z-gAxia1
-~ A+-% g-TraiisI'P.9 25g
.Note 1:14 Note 1.14 2Op-si, -10OAFA' -5g D2________Cofiet DropTR-10 75g Side Drotp!A Note 1 -20 psig -- -1O0"PA 75g 1)TR-1i :Top orBottonsEnd DropaA":
Note: 1. 12 20 psig -1C0XIP' ________
- gEndDrop D
A Calculation No. 11042-0205 Revision No. 2AR EVA Calculation Page 39of87Table 1 (Continued)
-Summary of Design Basis Load Combinations for the 61 BTH DSC [Ref. 5.8]HS.MLOADING
-Cokd -Cask X --I0u15 prig --÷80 ip -- AiLD-2 Nonsul Lo adina -Hot .-Cask X [- -101l5 ssik "- O0'oFCask
+S0Kip --I.D-3 ,Cask X -" -l0/!5 prig ,117=F. ÷80 Kip -ALD-4 Off-Norman L~oadisa-'Co1d C~ask X -- -20prig -0=FCas * +.SO ip FE ,BL.D-50ff-Nomnal osdin -Hot Cads X --2Op.sig -"l0O'f Kip FE BLD-4 Cask X 20- psig -. ÷40 Kip IFE Bw/shded'HSM STOR-AGE fIlorizontiI 3W' "Vetrtical 13W Internal External
,Thermal Other Service'_________________
DSC FPuel lTiel Pressuret' 5 Pr~essuretO)
- Conditiion Loads -Loads L~evel'HSM-I Off-Normsal HSM X .. .. 15psie" --:40'FHSM
-, -" BHSMl-2 NormalStorage HSM X -l5psiF -- 0'.FHSM -: -- A:HSM-30ff-Xormsl HSM X -, -!Spsig -I17'FHSM
-- -B*.HSM-4 ff-Normal Temp. -Fsiiled Furl ,HSM X --20 psig -
-EF :C,HSM-5 Blocked Venf Storage HSY, X. --65/120 pria -117' F --- f-HSM-6 B.V. +.Failed Fuel Stora~e HSM X .-65SI20ptig
-- *
-FE 13_HSMBV~x" HS.M-7Es-rdsquake -Cold -HSM X -10115 ps.F -0O FHS*M -EQ: !Y:HSM-8 Ear hquakaLoadi.ma-Hot HSM X -- -10/I -- 100'FHSM
-- EQ C;D#HSM-9 FloodtLoad (50' H-O)--Co~d HSM X --!0/lS~psig 22 psig 0= FHSM -Flood* CHSM-10 Flood ibad (50' H2O)-Hot ,HSM. X --10;15pss 22,psig 100FHSM -: FloodO CH$* 'ODIGHorizontalfllW VereicalflW Internal External Thermal Other Se-wera" -" DSC Fuel DSC 'Fuel Pre~sure('q Condiion
,Loads- .Lotds .o leelUL-1 Normal Urdoadiag-Cold HSM X --- 10/15 p-ig -0"F HSM -+60 Ksp A,UL-2 NormallUnloading-Hot HSM X --- I0f15 psig -100 F HSM +~60 Ki'p -A-"- HSM X -- -- 10{1 prig -117* F -60 Kip -- AUL-4 Off-Nonnal U~iioading--Cold HSM. X --" 20 poig -0° F HSM +60 FE BU'L-5 Off-'ormal Ltiloadi-Hot X --20 psig -100' F HSM +60Kip FEUL-6 HSM X, --- 20psi.¢ -l!7'F' +60Kip FE B.UL-7 Off. Norm.
11 HSM X --- 20 psig -wf/slsade
-*10 K.p FE CUS~gAccideit Unloodlng
-FEFAot,'
HSM X,-.- 65I120wp,!ig
'-
-S0 Hip FE DRF-1 DSC Refiood --Ca.s: .X 20 psia(ns.x)
Hyvdrostatic 120o'F Cask --13 A Calculation No. 11042-0205 Revision No. 2AR EVA Calculation Page 40of87Table 1 (Concluded)
-Summary of Design Basis Load Combinations for the 61 BTH DSC [Ref. 5.8](Notes for the preceding portions of Table 1)1. 25g and 75g drop acceleration includes gravity effects.
Therefore, it is not necessary to add an additional 1.Og load.2. For Level D events, only maximum temperature case is considered.
(TIhermal stresses are not limited for level D events and maximumtemperatures give mininium allowables).
- 3. Flood load is an external pressure equivalent to 50 feet of water.4. BV = HSM vents are blocked.5. At temperature over lOO°F a sunshade is required over the Transfer Cask Temperatures for these cases are enveloped by the l00OF (withoutsunshade) case.6. As described in Section T.4, this pressure assumes release of the fuel cover gas and 30% of the fission gas. Since unloading requires theUSM door to be removed, the pressure and temperatures are based on the normal (unblocked vent) condition.
Pressure is applied to theconfinement boundary.
- 7. As described in Section T.4, this pressu~re ass~umes release of the fuel cover gas ,and 30% of the fission gas. Although unloading requires theUSM door to be removed, thle pressure and temperatures are based on the blocked vent condition.
Pressutre is applied to the shell, innerbottom and inner and outer top cover plates.S. Not used.9. Unless noted otherwise, pressure is applied to the confinement boundary.
10 psig and 65 psig are applicable to Type 1 DSC, while 15 psig,and 120 psig are applicable to Type 2 DSC.10. Fuel deck seismic loads are assumed enveloped by handling loads.11. Load Cases UL-7 and UL-S envelop loading cases where the stresses due to insertion loading of SO kips are added to stresses due to internalpressure (in reality, the insertion force is opposed by internal pressure).
- 12. The 60g top end drop and bottom end drop ,are not credible events, therefore these drop ,analyses are not required.
- However, consideration of60g end drop and 75g side drop conservatively envelops the effect of 25g corner drop.13. Conservatively based on normal operating pressure times 1.5 to cover future 1OCFR Part 71 requirements.
- 14. A 25g corner drop analysis (300 fr-om horizontal) of 61BTH DSC without support from the TC is to be documented.
- 15. Service Level C is fo the standard seismic event and Service D is for the high seismic event.
A Calculation No. 11042-0205 Revision No. 2AR EVA Calculation Page 41lof 87Table 2 -Internal Pressure in the 61 BTH Type I DSCMaximumDesign PressureMaximumused in Ref. 5.2Calculated adTiDesign Condition Pressure Calculation Reference
[psi] [psi]Normal 7.3 10 Ref. 5.8 Table T.4-16Off-Normal 10.9 20 Ref. 5.8 Table T.4-20Accident 56.1 65 Ref. 5.8 Table T.4-24Table 3 -Maximum Temperatures in the 61 BTH Type 1 DSC ShellMaximum DesignTemperature Design Calculated used in ThisCondition Temperature Calculation Reference
[OF] [OF]Storage 374 500Normal Ref. 5.8 Table T.4-13Transfer 439 500Storage 399 500Off-Normal Ref. 5.8 Table T.4-18Transfer 416 500Storage 611 625Accident
________________________
Ref. 5.8 Table T.4-22Transfer 467 500 A Calculation No. 11042-0205 Revision No. 2AR EVA Calculation Page 42of87Table 4 -Properties of SA-240 Type 304. [Ref. 5.11]Sm SuE Aloal ~ liae Yield Stress for SL A/B Yield Stress for SL DTemp Modulus of Strwabes Yil Strs Utenilae Limit Load Analysis Limit Load Analysis[°F] Elasticity Inest ki tegh(Note
- 1) (Note 2)[s][ksi]
[ksi] [ksi] [ksi]70 28,300 20.0 30.0 75.0 30.0 46.0100 28,138 20.0 30.0 75.0 30.0 46.0200 27,600 20.0 25.0 71.0 30.0 46.0300 27,000 20.0 22.4 66.2 30.0 46.0400 26,500 18.7 20.7 64.0 28.1 43.0500 25,800 17.5 19.4 63.4 26.3 40.3600 25,300 16.4 18.4 63.4 24.6 37.7625 25,175 16.3 18.2 63.4 24.5 37.5700 24,800 16.0 17.6 63.4 24.0 36.8(1) The yield strength to be used in a Limit Analysis for Service Level A and B Loading is 1.5"Sm, per Paragraph NB-3228.1 of Reference 5.7.(2) The yield strength to be used in a Limit Analysis for Service Level D Loading is the lesser of 2.3*Sm and 0.7*Su, perParagraph F-1341.3 of Reference 5.9.
A Calculation No. 11042-0205 Revision No. 2AR EVA Calculation Page 43of87Table 5 -Properties of SA-36. [Ref. 5.111E SmYield Stress for SL A/B Yield Stress for SL DTemp Modulus of Allowable SV Ultimate Limit Load Analysis Limit Load Analysis[°F] Elasticity Srs YilStes Tnle(Note
- 1) (Note 2)[ksi] Intensity
[ksi] Strength
[ksi] [ksi][ksi] [ksi]70 29,500 19.3 36.0 58.0 29.0 40.6100 29,338 19.3 36.0 58.0 29.0 40.6200 28,800 19.3 33.0 58.0 29.0 40.6300 28,300 19.3 31.8 58.0 29.0 40.6400 27,700 19.3 30.8 58.0 29.0 40.6500 27,300 19.3 29.3 58.0 29.0 40.6600 26,700 17.7 27.6 58.0 26.6 40.6625(3) 26,400 17.6 27.2 58.0 26.4 40.4700 25,500 17.3 25.8 58.0 26.0 39.8(1) The yield strength to be used in a Limit Analysis for Service Level A and B Loading is 1.5*Sm, per Paragraph NB-3228.1 of Reference 5.7.(2) The yield strength to be used in a Limit Analysis for Service Level D Loading is the lesser of 2.3*Sin and 0.7*Su, perParagraph F-1341.3 of Reference 5.9.(3) All values are interpolated from the 600 0F and 700 °F values.
A Calculation No. 11042-0205 Revision No. 2AR EVA Calculation Page 44of87Table 6 -Summary of Load Cases and ResultsRequired!
Required Pressure toTep. Aalss esg Pesue Satisfy Limit Load Limit Load Code LimitName Mesh Level Loading Criteria Collapse Pressure Load Criteria[e] Ciei Nt )(NoteS5)
[psi] Satisfied?
[psi][piInitial Inena S00 SL A/B 32 60.0 95.9 YesPressureInternalRefined Welds 500 SL A/B 32 60.0 94.8 Yes.u PressureReieaedndLd ItraEPressure 500 SL A/B 326009.Ye Internal 625 SL D 65 90.2 136.6 YesPressureRefined Welds and Lid Inenl 625 SL D 65 90.2 132.6 YesPressureInternalRefined Welds and Lid 500 SL A/B 32 60.0 93.7 YesN PressureRefined Welds and Lid Internal 625 SL D 65 90.2 132.9 YesPressureRefned eld andLid InternaleiedWls n i Pressure 500 SL A/B 32 60.0 94.5 YeaE :Initial Internal5= (NtS)Pesr 00 SL A/B 32 60.0 97 YesRequired G-Load toTep nlss Required!
Satisfy Limit Load Collapse Code LimitName Mesh Level Loading Criterinayi Design G-Load Criteria G-Load Load Criteria[~][g] (Note 5) [g] Satisfied?
Initial Side Drop 500 SL D 75 104.0 181 -Yes~RefinedCircumferential Mesh Side Drop 500 SL D 75 104.0 185 YesE~(Note 3)"z Initial Side Drop with 50 S e(Note 3) Off-Normal I P 7 0.8~ ~Initial Side Drop 500 SL D 75 104.0 189 YesNotes:1) [Not Used]2) The 97 psi colipase load is estimated!/
extrapolated from the final obtained solution at 95 psi. Excessive run times make more precise results impractical.
- 3) The reported collapse load is conservative
-based on the collapse behavior it is expected that smaller analysis time steps would yield larger collapse loads. Thisiwas deemed impractical due to the long rntime and the large magnavailable.
- 4) The Service Level A/B required pressure of 32 psi is based on the design internal pressure of 10 psi plus an equivalent internal pressure of 22 psi which accounts forother Service Level A loads. See Section 4.25) Includes both the ASME code required factor and the 0.8 weld strength reduction factor to account for examination limitations.
A Calculation No. 11042-0205 Revision No. 2A R E VA Calculation Page 45 of 87 ]Figure 1 -Sketch of the 61 BTH DSC Top End and Transfer Cask from Reference 5.1
A Calculation No. 11042-0205 Revision No. 2AR EVA Calculation Page 47of87Vent/Siphon
~frATDSCFigure 3 -ITCP and OTCP Closure Weld Details from Reference 5.5(Note: Section B-B above is rotated 90 degrees clockwise from Figure 2.)
A Calculation No. 11042-0205 Revision No. 2AR EVA Calculation Page 49of8710+16 v-2+7315+31- 6 F+18+28+29/t //8+11+24+25+262 OTCPWeld Metal(ITCP Not Shown)Figure 5 -OTOP Flaws -Raw Data from Reference 5.1(Numbers indicate the Flaw Number from Reference 5.1)REPRESENTATIVE GROUP FLAW20 (0.32L Lx .07 H, NEAR 19 IND Weld MetalEXTENTS OF, GROUP FLAW REGION3 (0.18 L x .09"H, REMOTE FROM ALL BUT #2)14(./x.7 EOEFO L #2)
CIRCUMFERENTIAL, INTERMITTENT)
OTCP(ITCP Not Shown)Figure 6 -OTOP Flaws -Main Flaw Group Reduced and Bounded A Calculation No. 11042-0205
., Revision No. 2AR EVA Calculation Page 50of87-2+140,195"H, FULL CIRCUMFERENTIAL OTCP(ITCP Not Shown)Weld MetalFigure 7 -OTCP Flaws -Bounding Set #1 for ANSYS Collapse Analysis0.14"I-,Weld Metal-14 ,0.07 'H, FULL CIRCUMFERENTIAL OTCP(ITOP Not Shown)Figure 8 -OTCP Flaws -Bounding Set #2 for ANSYS Collapse Analysis A Calculation No. 11042-0205 Revision No. 2A R EVA Calculation Page 51 of 87/OTCPWeld Metal </--ENVELOPE OF ALL BUT 7 AND 11/MAX INDIV. HEIGHT IS 0.09'MAX INDIV. LENGTH IS 2.09"/ 11 (7.17"L x .09"H, REMOTE FROM ALL OTHERS,-- 7(10.34"Lx 0.11"H, NEAR FLAWS 1 THRU 9)ITCPFigure 9 -ITCP Flaws -Raw Data from Reference 5.1OTCPREPRESENTATIVE GROUP FLAW0.09"H, FULL CIRCUMFERENTIAL Weld MetalITCPFigure 10 -ITCP Flaws -Bounding Flaw Set for ANSYS Collapse Analysis A Calculation No. 11042-0205 Revision No. 2A R EVA Calculation Page 52 of 87ANSYS 14.0MAP. 25 2015~10:39:49 ELEMENTSE1ACET=ITYPE MUMSV -1DIST..53.
912XE -16.9125YE -146.68Z-BUFEERFigure 11 -Overview of the Axisymmetric ModelANSYS 14.0MAR 25 201510:37:29ELEMENTSPowerGiaphics E!ACET-1TYPE MUMZV -1 4 9924'XE -31.9932,YE -190.433Z-BIrEEER iiIFigure 12 -Mesh Details Near the Lid Regions of the Axisymmetric Model(Small differences in the mesh exist amongst the sub-models)
A Calculation No. 11042-0205 Revision No. 2AR EVA Calculation Page 53of87~MAR 25 201510:44 :56ELEMENTSPowerGraphics EFACET-1TYPE MUMZV -1*D1ST-1.27014
- XF "32. 6238Y7F -194.594Z -BUFPERFigure 13 -Mesh Details at the Welds for Axisymmetric Case #1ANYS 14.0MAR 25 201510:45:35ELEMENTSPowe EFACET.ITYPE MUM2V -1*DIST-] .210144238*YP 2 PUFFERPEDGEFigure 14 -Flaw Locations for Axisymmetric Case #1 A Calculation No. 11042-0205 Revision No. 2AR EVA Calculation Page 54of87ANSYS 14./APR21 S.111:.44' 5P er~raphlcs EFACET-iTYPE NtJMZV =1*DIST-3.
78456*xF =30. 6287*YF -193, 906Z-BUJFFER Figure 15 -Refined Mesh (Weld Region) for Axisymmetric Case #1AY =194.043 Figure 16 -Refined Mesh (Weld and Lid Interior Region) for Axisymmetric Case #1 A Calculation No. 11042-0205 Revision No. _ __AR EVA Calculation Page 55of87 .......A9% 21 2Q1514=52=50/
0D2PT2. 5 .=30. 7029*YF 2,-3thFkER S.Figure 17 -Mesh Details at the Welds for Axisymmetric Case #2ANSYS 14.0MR25 2015I0: 38:01PowerGraphics EFACET-ITYPE NUMZV -i*DIST-1.
23285*XF -32.7711*YF -194.*652 Z-BUJEEREDGEFigure 18 -Flaw Locations for Axisymmetric Case #2 A Calculation No. 11042-0205 Revision No. 2AR EVA Calculation Page 56of87ANSYS 14.0APR 21. 201514:38:39ELEMENT SPove rGraphi c sEFACET=ITYPE NUN~YV =-iDIST-53.
812i YE =16.8125ZE =146.88iZ-BUFFEP.
EDGEFigure 19 -Overview of the Half-Symmetry Model A Calculation No. 11042-0205 Revision No. 2A R EVA Calculation Page 57of87t41 =5 1 1 M AY 1a43 1 04=3EPACET I A? .114t ic0 4
- 013.(c))WLid Region Solid Viewwt lw iil (b) Wlid Region MeshFigure 20 -Detail Views and Mesh Plots of the Half Symmetry Model A Calculation No. 11042-0205 Revision No. _ __A R E VA Calculation Page 58 of 87ANSYS 14.0MAY 4 201510:04:39ELEMENTSPowe rtraphi csEFACET-1TYPE NUNXV --.373471YV --. 679348ZV -. 631669*D1ST-35.7666
- XF -12.2625*YF -29.0871*ZB' -167.331A-ZS-41 .1572Z-BUJFFER ANSYS 14.0MAY 4 201510:04:54ELEMENTSPowe rGraphics EFACET=ITYPE NUNXV --.373471YV =-. 679348ZV -.631669*DIST-4.
86713i*X =44.6517*YF -26.3642*Z1 =167.034A-ZS-41 .1572Z-BUJFFER Figure 21 -Isometric Views of Half-Symmetry Model A Calculation No. 11042-0205 Revision No. 2A RE VA Calculation Page 59 of 87AflSYS 14.0NAY 4 201510: 22:15ELEM4ENTS Powe rGraphlc3 EYACET-ITYPE NUNXV --. 290468Yv --.652092Zv -. 700299*YF -30. 5999A-ZS=33.8 517Z-BUFFEJR ANSY$ 14.0NAY 4 201510 22=3TYPE NUNXiv -.290460Yv=.652092 2'. .700299' DIT-4.Th42b
- XF 43.55'TY ;0952ZF 165. 967A-ZS-33.,851 712-BUFFERFigure 22 -Isometric Views of Half-Symmetry Model (Refined Circumferential Mesh)
A Calculation No. 11042-0205 Revision No. 2AREVA Calculation Page 60of8714,0 514373 14.5APR 21 2015~ APR 22 251510:41 12 1042:4152 NJOSAL 50O5710O.
NODAL SOLTIJTON sos -t SITB -TOM,.- U T01ME=65EIMPLEQ2 IA E SPLEQV IAVG)SowerGraph es PeweerGIraphi csEFA'57-.
£5A05T-SAVF.SS.Nst S3M .29 1142 -.84 7490141
- 002675 1542 -.047214.297E.3 Il.O05 246m .020492.892E03 .015736.OO196 .020294.001783 .021423* 002376 % .0419 6m .002 75 m .047214(a) Equivalent Plastic Strain in Weld Region [in/in] (b) Equivalent Plastic Strain in Weld Region [in/in]at 20 psi Internal Pressure at 65 psi Internal Pressure~gTS14.OAllOT 14.0APK TI 2405 APR 22 2025ASi121 0:41:541 0:42:53IIDA ODAL SOLUT1ION
( SS" =18SUB -19S IME..95 .9249 5E 5.94K~PLSQV (AVISEPLQ; A~I 5- rls aphr cr Sewe toaphicoEFACE?=I E ACN'I" 1DX 52.0242 2942 ..51.026759221{ -1.36914 0242 -1.25092,2J 1.1230.22794197 1.4971 ) .555962S 2.24612 .8339432.6204 7 .9729142.99403 1 I.112921: 36918 1.25091(c) Equivalent Plastic Strain in Weld Region [in/in] (d) EQV Plastic Strain in the Cover Plates at 95.9at 95.9 psi Internal Pressureps Figure 23 -Results for Axisymmetrdc Case #1 -Initial Mesh -Service Level A/B(Note that the magnitude of the strains and deflections has no true physical meaning due to the nature of limit loadanalysis)
A Calculation No. 11042-0205 Revision No. 2AR EVA Calculation Page 61 of 874A:4'4AN3iS .1!ZV *13r -.5Z 4U S4 J4.~44 *~1;~~ 1.1254,4 44.744 ,4,.~:4,4.44; iaiii 47.54.4 41.141 44.137TCCrno; f,,i',nar.
Ipti>(a) Service Level A/B Material Properties 4..:4:4:74'CDIS'T-.T5 YF 1.54, !UFFER4 Z~275~ 4,5.71.54 41414 7772574, U4.44~444 47.477 44.:745 44,4417 127.111r tsu. P$1J(b) Service Level D Material Properties Figure 24 -Deflection History of the Center of the OTCP for the Axisymmetrdc Case #1 Initial Mesh(Maximum deflection occurs at the center point of the lids, in the outward axial direction)
(Note that the magnitude of the deflections has no true physical meaning due to the nature of limit load analysis)
A Calculation No. 11042-0205 Revision No. 2A REVA Calculation Page 62of87E[PLE\ tAWS EFPLEQV d .WA-Al l I .0921454AJVAE 4.ET- EII'A -14.2)A .404.171 I IIA .512r~44.
311 -I 11T1414.74114.103 4A~73 1APEV A(c) quialet Plsti Stainin Wld egin [i/in (d EQ PlaticStrin i th Coer latesat 4A.8at 94. psi Iterna Pressre psFigure 25 -Results for Axisym metric Case #1 -Refined Mesh in Weld Region -Service Level A/B(Note that the magnitude of the strains and deflections has no true physical meaning due to the nature of limit loadanalysis)
I A Calculation No. 11042-0205 Revision No. 2AR EVA Calculation Page 63of87FAm.2YS I4. 11A-i SIUB -7AVP, S4LM~ti 4(a) Equivalent Plastic Strain in Weld Region [in/in] (b) Equivalent Plastic Strain in Weld Region [in/in]at 20 psi Internal Pressure at 65 psi Internal PressureI~ *I SrE8E .i41AC flL.V A'05"xl(c) Equivalent Plastic Strain in Weld Region [in/in] (d) EQV Plastic Strain in the Cover Plates atat 93.7 psi Internal Pressure 93.7 psiLevel A/B(Note (c) and (d) are plotted one timestep before the collapse pressure)
(Note that the magnitude of the strains and deflections has no true physical meaning due to the nature of limit loadanalysis)
A Calculation No. 11042-0205
- , Revision No. 2AR EVA Calculation Page 64of8772.117., 1.2 I.* I --.....l. .J! 2157DI*Th -19;XFv I.(a) Service Level A/B Material Properties Refined Mesh at Weld Region Only71 I I .7 40 7,71 1171 O.2T.n7. 9 I; II RI:'l' 19.104 14')POST2.?YFv -,5-VI 11,5~Refined Mesh at the Weld and Lid Interior RegionsFigure 27 -Deflection History of the Center of the OTCP for the Axisymmetric Case #1 Refined Mesh(Maximum deflection occurs at the center point of the lids, in the outward axial direction)
(Note that the magnitude of the deflections has no true physical meaning due to the nature of limit load analysis)
A Calculation No. 11042-0205 Aevi_ _ _ __.AR EVA Calculation Page 65of873 0 -.. ..... .In itial M e sh ... .... .... ... .... .... .... ..... .. ....... l .. ................Initial Mesh, Enhanced ormul atios !I2.50 -.... ..-Refined M esh (W elds) ... .... .... ... ..... ... ... .. ..... ..-ii-- -Refined Mesh (Welds end Lids) =E /15 .-. t."!050 i --.... .... ... -- --..... .........
.... ...... .... ..0.00 --.. ..............
-... .........
..... ...-- ... ------.... ........
.. .. .. ...............
.... --...........
... -.. ...... .... .... .... .020 40 60 80 100 120Internal Pressure
[psi]Figure 28 -Comparison of Maximum Displacement Histories for Axisymmetric Model Sensitivity Studies(Maximum deflection occurs at the center point of the lids, in the outward axial direction)
(Note that the magnitude of the deflections has no true physical meaning due to the nature of limit load analysis)
(Service Level A/B materiaJ Properties)
A Calculation No. 11042-0205
- Revision No. 2AR EVA Calculation Page 66of871.80 i --/ 116 -I -I I~00 -0 2 0 0 5 0 0 8/9 0Internal Pressure
[psi]Figure 29 -Comparison of Maximum Displacement Histories for Axisymmetric Model with Lid ContactDefined using Nodal DOF Couples vs. Contact Elements(Maximum deflection occurs at the center point of the lids, in the outward axial direction)
(Note that the magnitude of the deflections has no true physical meaning due to the nature of limit load analysis)
(Service Level A/B material properties)
A Calculation No. 11042-0205 Revision No. 2A REVA Calculation Page 67 of 87S 4.00 , .. .. .3 .00 * ..... Re fined Mesh (W'elds and Lids) : .....i i RefinedvMesh (Welds andtLids) w/ Pressure onlTCP Crack Faces' i ,-i100i 0 10 20 30 40 ,50 60 70 80 90 1001 Internal Pressure (psi)Figure 30 -Comparison of Maximum Displacement Histories for Axisymmetric Model With and WithoutPressure Loading Applied to the ITCP Weld Root Flaw Faces(Maximum deflection occurs at the center point of the lids, in the outward axial direction)
(Note that the magnitude of the deflections has no true physical meaning due to the nature of limit load analysis)
(Service Level ANB material properties)
A Calculation No. 11042-0205 Revision No. _ __AR EVA Calculation Page 68of87mU.IE,:V R,+f+)LW IA,'.0 -I'.03.0.+g+,
47m +O +3 m .0*(a) Equivalent Plastic Strain in Weld Region [in/in] (b) Equivalent Plastic Strain in Weld Region [in/in]at 20 psi Internal Pressure at 65 psi Internal PressureI Is ; o, + ,o+1U07I :f#C .-+' :+,++++,1ALOOU~VZ E' -!+101(c) Euivaent lastc Stain n Wel Regon [n/in (d)EQV lasti Strin i theCove.Plaes aat 936 ps Intenal ressue 936 psLevel A/B(Note (c) and (d) are plotted one timestep before the collapse pressure)
(Note that the magnitude of the strains and deflections has no true physical meaning due to the nature of limit loadanalysis)
A Calculation No. 11042-0205 Revision No. 2AR EVA Calculation Page 69of87AP.2 0- Apa .1 I. : I OLU U'EEACET' 3a .at2Ps nenlPesr at 64psPIntrna2Presurt t-+ SNX P,.,13¶1I I Ii .t 040!s J(a) Equivalent Plastic Strain in Weld Region [in/in] (b) EquValn Plastic Strain in Whel Roegionte aininat 240 psi Internal Pressure at65pi .nera presurFigure 32 -Results for Axisymmetric Case #0 -Refined Mesh in Weld and Lid Interior Region -ServiceLevel A/B(Note (c) and (d) are plotted one timestep before the collapse pressure)
(Note that the magnitude of the strains and deflections has no true physical meaning due to the nature of limit loadanalysis)
I A Calculation No. 11042-0205 Revision No. 2AR EVA Calculation Page 70of874.003.503.002,500S2.00E1.501.000.50-AAxisynmretric Case #1......................
Axisymnmetric Case #2---- -- Axisymmnetric Case #0 (No Flaws):'Ii IIII/i#I//I /I/0.0010 20 30 40 50Internal Pressure
[psi!60 70 80 90 100IFigure 33 -Comparison of Maximum Center-of-Lid Displacement Histories for the Various Flaw Models(Service Level NB material properties)
(Note that the magnitude of the defiections has no true physical meaning due to the nature of limit load analysis)
A Calculation No. 11042-0205 Revision No. 2A REVA Calculation Page 71lof 87?W(A 4 ZO I M(AY4 2015gF&CET- TI.4 I0 (a) Equivalent Plastic Strain in Weld Region [in/in] ()Euvln lsi tani edRgo i/nat 20 psi Internal Pressureat6psInealPsur (d) QV Pasti Stain n th CovrPates at 9 psFigure 34 -Results for Half-Symmetry Case #1 Internal Pressure Loading Benchmark Analysis
-ServiceLevel NB(Note that the magnitude of the strains and deflections has no true physical meaning due to the nature of limit loadanalysis)
I A Calculation No. 11042-0205 Revision No. 2AR EVA Calculation Page 72of874.(003.503.00'E 2.500IIE 1.50t1.000.50-0.00 1050Internal Pressure
[osij100Figure 35 -Benchmark of the Half Symmetry model with the Axisymmetric Analysis(Service Level NB material properties)
(Note that the magnitude of the deflections has no true physical meaning due to the nature of limit load analysis)
A Calculation No. 11042-0205 Revision No. 2AR EVA Calculation Page 73of87Atl4r'3+3 14, KAY I 402501 11 5t,2 ,1 4 I0+TIKE-(a qivln (o iss trs ps]a 7g () qialn latcStani Wel Regio [in/in3455 93at4475g.
M.4Y ++ :O£,0 1 ll-59 a8"4. .15234(c31iaet vnMss Srs pi at' 18g () qiaen3lstcSrini el4ein i/nE 54at4181g.
Figure 36 -Equivalent Stress and Plastic Strain Plots from the Half-Symmetry
- 1 Side Drop Analysis(Note that the magnitude of the strains and deflections has no true physical meaning due to the nature of limit loadanalysis)
I A Calculation No. 11042-0205 Revision No. 2AR EVA Calculation Page 74of8714.0MAY 4 201510: 55"29N4OfAL SOLOTIOMSTEP-4SU]9 TIMSE 5USUK iAV')RSYS1-0KFACET-1AVP.E3-Mat DM4X -. 340414i SMX -.340414i m -:40414(a) Deformed Shape Plot -Axial View -Exaggerated Scale14.0MAY 4 0= 5V:19~NODAL TIME-1S0, 545USUM lAVG)\ ~Powz'Gt aphi csEFACET-I DMX -404140.0Y624.075047m.113471.226942(b) Deformed Shape Plot of DSC Shell -Axial View -Exaggerated ScaleFigure 37 -Additional Results Plots from the Half-Symmetry
- 1 Side Drop Analysis(Note that the magnitude of the deflections has no true physical meaning due to the nature of limit load analysis)
A Calculation No. 11042-0205 Revision No. 2AR EVA Calculation Page 75of87AIOOY 14.014A1 4 14 15', ¶5314 4 SOLI.1 212TI1975 EPPL2 [AVOI' []~2 1142 97 .01 V2124...I0-1.01* 4 .0 O'D& L I 34.1.1 419(c) quialen (vn Mses)Strss psi]at 81g (b) Equivalent Plastic Strain in Weld Region fin/in]at 18114.Figure 38 -Equivalent Stress and Plastic Strain Plots from the Half-Symmetry
- 1 Side Drop Analysis withOff-Normal Internal Pressure(Note that the magnitude of the strains and deflections has no true physical meaning due to the nature of limit loadanalysis)
I A Calculation No. 11042-0205 Revision No. 2AR EVA Calculation Page 76of87O1514021 4.013t147 01 13434,14 IOI04Q'. A2,07AVE..e 0220,9 -7,io. 3319.02 4441 !,3046. 944402,(a) Equivalent (von Mises) Stress [psi] at 75g (b) Equivalent Plastic Strain in Weld Region [in/in]at 75g.433 ' 14.0 3 149 20MA 3 215 : 53 3OLJT7ON3 L 5?LU 71 TP4 1' , 420m)AVPF .41e3m0403
-.132032343 1 ..45347408 51325 447 1 04.4-1441atU185g.
Figure 39 -Equivalent Stress and Plastic Strain Plots from the Half-Symmetry
- 1 Side Drop Analysis withRefined Circumferential Mesh(Note that the magnitude of the strains and deflections has no true physical meaning due to the nature of limit loadanalysis)
Calculation No. 11042-0205 A _____Revision No. 2AR EVA Calculation Page 77of8714171 14 ..14.14 14 444. 4 2~7J51441 4 211I ~74 ~ ?44XAL SOtU11744 1J74A1. I. 1754?r1141-1114~ -44 141~ s.2za~e1754 5 A 1.444414.75 4114L 1444 .47 '154 1 144. ..A4734~144.-'s 474S442 2~4~ .'7~42743 4 114it*-'724~4.5
-0155 157~l.i-1752~ 34119~.4
-147242-L(b) Equivalent Plastic Strain in Weld Region (in/in](a) Equivalent (von Mises) Stress [ps~ at 75g at 75g.N4221 14.27411211 ~ 14741 4 201114,21 4 2274I15XAL 111011SIXPL 4)L1TIIN Sti5~.45715 ~'l751515457V IAVI471{72~ ~flA1~5't'l 741514.4isi 4.41757-44,75 0117 ~.O11014441 4441 ~,1.41744 1441 -1142.410 11541 '-512) i744)771042.319 4441122122342 40.11-72114 3 2:474271-14414,',
-~-144201.4 4,4444,~. ~ 54~49ti-I(d) Equivalent Plastic Strain in Weld Region (in/in](C) Equivalent (von Mises) Stress [psi] at 189g at 189g.rigure 'iu -c-qulvaieni ano riastic r',ots irom me ra-wr-ymmeiry W.u Uro r-iawsj oiue uropAnalysis(Note that the magnitude of the strains and deflections has no true physical meaning due to the nature of limit loadanalysis)
A Calculation No. 11042-0205 Revision No. 2AR EVA Calculation Page 78of870.800.700.60"2E 0.50.20.40S0.300.200.100.000 20 40 60 80 100 120 140 160 180 200Acceleration
[g]Figure 41 -Comparison of Maximum Displacement Histories for the Various Half-Symmetry Analyses(Service Level D material properties)
(Note that the magnitude of the strains and deflections has no true physical meaning due to the nature of limit loadanalysis)
A Calculation No. 11042-0205 Revision No. 2AR EVA Calculation Page 79of8710.0 APPENDIX A -ADDITIONAL RESULTS PLOTS AND ELASTIC-PLASTIC ANALYSESPurposeThe purpose of this appendix is to document elastic-plastic analyses of DSC-16. The models listed beloware used as a basis for the analyses.
Axisym metric 1 with Refined Welds and Lid MeshHalf Symmetry 1 with Initial MeshThese models produced the bounding results using the limit load methodology.
The models are updated toinclude the elastic-plastic material properties described below. In addition, these new runs consider theeffects of large-deformations (NLGEOM,ON).
The intent of these analyses is to provide a more realistic prediction of the actual material strains that would occur under the design basis loading, as opposed to theover-estimated strains and deformations which result from the limit-load analysis methodology.
Material Properties The elastic-plastic behavior of SA-240 Type 304 stainless steel is idealized using Ramberg-Osgood stress-strain curve constants calculated using the equations in Appendix B of Reference A12.The constants arecalculated using the ASME code [Ref.5.16]
specified minimum yield and ultimate strength values at theapplicable temperatures.
In order to incorporate the curves into the ANSYS analysis, the initial slope of thecurves must match the defined elastic modulus.
Therefore the first data point in the curves is defined at the(strain,stress) data point (Sy/E,Sy).
The material behavior is based on true stress and true strain, since theANSYS analysis accounts for changes in geometry (e.g. necking).
The following equations from Reference Al were used to develop the curves:E = true strainEo = true strain at yield= true stress6o= true stress at yield1n=ln(1 + eu)ln(1 + eu) c'u(1 + eu) IIau(1 + eu)eu= engineering strain at the ultimate tensile strength2 The euations in Reference Al to develop the full-range true stress-strain curve are based on curve fits of tensile testdata. The resulting curve is not indicative of a specific failure type or analysis approach.
Rather, it is a method todevelop a full-range Stress-strain curv of a material using a limited set of data (iLe. minimum specified yield andultimate strengths.)
A Calculation No. 11042-0205 Revision No. 2AR EVA Calculation Page 80 of 87637= engineering yield stressau= engineering ultimate strengthThe value of eu is taken as 0.35, which is assumed to be etot-0.05, where etot is taken as the minimumspecified elongation of the material (40%), per Reference A2.The relationship between true and engineering stress and strain is per the following equations:
6true = aeng(eeng
+ 1)Etrue = ln(eeng + 1)Figure A-I shows both the true and engineering stress strain curves based on the Ramberg-Osgood equations.
The curves at various temperatures as coded into the ANSYS analysis are shown in Figure A-2.The SA-36 shield plugs use a bi-linear stress strain curve with a tangent modulus of 1% of the initial elasticmodulus.
This results in a less stiff representation of the shield plug, which will result in conservatively greater strains in the DSC.Load CasesAnalyses are performed for the following load cases:1. Internal pressure loading (32 psi) for Service Level A/B.2. Internal pressure loading (65 psi) for Service Level D.3. Side drop Loading (75g) for Service Level D.As discussed in Section 4.2, these three load cases bound all of the design loading conditions for the DSCOTCP and ITCP welds.Results and Conclusion Plots of the equivalent plastic strain for the three analyses are shown in Figure A-3 through A-5. The resultsare summarized in Table A-I. As shown by the results, the strain levels remain well below the minimumspecified elongation limits of Type 304 steel and Type 308 weld electrodes
[Ref. A2 and A3]. Therefore, material rupture will not occur at the design conditions.
A Calculation No. 11042-0205 Revision No. _ __AR EVA Calculation Page 81lof 87References Al. EPRI NP-5531.
Evaluation of High-Energy Pipe Rupture Experiments.
January 1988.A2. ASME Section II Part A. Ferrous Material Specifications.
1998 Edition with Addenda through 2000.A3. ASME Section II Part C. Specifications for Welding Rods, Electrodes, and Filler Metals. 1998 Editionwith Addenda through 2000.Cornputer FilesAnalyses performed on Computer HEA-0213A using ANSYS Version 14.0 [Ref. 5.6]File date & time listing is as displayed by the Windows 7 Operating System -Differences may occur due tolocal time zone and daylight savings settings.
Analysis Case File Name Date & TimeElastic-Plastic 61BTH_WeldFlaw_1F A)( 2 DETACH.db 11/29/2015.8:41 AMAxisymmetric 1Refined Lids and Welds 61BTH_WeldFlaw 1FAX_2_DETACH.rst 11/29/2015 8:18 AMInternal PressureSLNAB SOLUTIONAXISYMMIP_500F.INP 11/27/2015 4:03 PMElastic-Plastic 61BTH_WeldFlaw 1FAX_2_DETACH.db 11/27/2015 3:36 PMAxisymmetric 1Refined Lids and Welds 61BTH_WeldFlaw 1FAX_2_DETACH.rst 11/27/2015 3:33 PMInternal PressureSL D SOLUTION_AXISYMMIP_625F.INP 11/19/2015 11:46 AMElastic-Plastic 61 BTH_Weld Flawl 1GC.db 11/29/2015 8:16 AMHalf-Symmetry 1Initial Mesh 61 BTH_WeldFlaw 1GC.rst 11/27/2015 6:26 PMSide DropSL D SOLUTION_HALESYM_SD.INP 11/20/2015 10:08 AMStress-Strain Curve Development Stress-Strain.xls 11/30/2015 10:59 AM A Calculation No. 11042-0205 Revision No. 2AR EVA Calculation Page 82of87Table A-i -Summary of Elastic-Plastic Analysis Results.ValueAnalysis Case Result[in/in]Internal Pressure Equivalent Plastic StrainService Level A at 32 psi Internal 0.0183Axisymmetric Pressure (1.83%)(Note 1) (Note 1)Internal Pressure Equivalent Plastic Strain 009Service Level 0 at 65 psi Internal (5.97%)Axisymmetric PressureSide DropService Level D Equivalent Plastic Strain 0.0609Half-Symmetry at 75g Acceleration (6.09%)Note 1: The 32 psi internal pressure is bounding for Service Levels A and B and includes designinternal pressure of 10 psi plus an additional 22 psi to account for inertial loading of the DSCcontents onto the lid. See Section 4.2 for details.
A Calculation No. 11042-0205 Revision No. _ __AR EVA Calculation Page 83of8788 type 304100.00______
_____ ______ ___ __90.00 ____ _80.00 ______ _____70.00 ___-___60.00 ______40.00 _____ ___ __30.00 ---o-- Engineering
--_20.0010.00 ___ __0,00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40Straln, lnhlnFigure A-I -Ramberg-Osgood Derived Stress Strain Curve for SA-240 Type 304 at 500 °F.
A Calculation No. 11042-0205 Revision No. 2AR EVA Calculation Page 84of87ANSYS 14.91(PIN Thbi-e Eor:
1 NOV 27 201515: 25: 00zeuo 500.002V -1to, 01T1'-. 75Irs 5520.04 .1 .20 .30AZ40V1 14.01(IN TeJ021 Fot Naterial 1 NOV 27 201515:26i05Table flatST1= 500.O00zV *12000 DI S-.,75ZF .5l~o Z-BOtFEr1
.4 -8 2.2 1.6 2.2 .0 11.4q 1.0Figure A-2 -Ramberg-Osgood Stress Strain Curves for SA-240 Type 304 from ANSYS Model at VariousTemperatures.
Upper image shows full range of curves (strain axis from 0 to 40%), lower image shows detail of the initialyield point (strain axis from 0 to 0.2%)
A Calculation No. 11042-0205 Revision No. 2AR EVA Calculation Page 85of87ANSYS 14.0NOV 29 201508:41:45NODAL SOLUTIONSTEP=3SUB =iTIME=32EPPLEQV (AVG)Powe rGr aphi csEFACET=IAVRES=Mat DMX =. 307147SMX =.018306.002034i.004068.012204mm.016272
.018306Figure A-3 -Service Level A Internal Pressure
-Equivalent Plastic Strain at 32 psi *Note*Note: The 32 psi internal pressure is bounding for Service Levels A and B and includes design internal pressure of 10psi plus an additional 22 psi to account for inertial loading of the DSC contents onto the lid. See Section 4.2 for details.
A Calculation No. 11042-0205 Revision No. 2AR EVA Calculation Page 86of87ANSYS 14.0NOV 27 201515:36:02NODAL SOLUTIONSTEP=3SUB =3TIME=65EPFLEQV (AVG)Powe rGraphics EFACET=IAVKE S Ma tDMX ~=.818411 SMX =.0597290.039891.053093Figure A-4 -Service Level D Internal Pressure
-Equivalent Plastic Strain at 65 psi A Calculation No. 11042-0205 Revision No. 2A R EVA Calculation Page 87of87ANSYS 14.0NOV 29 201508:15:52NODAL SOLUTIONSTE P-2SUB -5TIME-75EPPLEQV (AVG)PowerGraphicsAVRES MatDMX =. 080518SMX =. 060912o.0406708
-.046091ANSYS 14.0NOV 29 201508:16:20NODAL SOLUTIONSTEP=2SUB -5TIME=75EPPLEQV (AVG)Powe rGraphics EFACET=1AVR.ES-lat DM( -. 033851SNX =.0548530006095.036569.042664.054853Figure A-5 -Service Level D Side Drop -Equivalent Plastic Strain at 75g.Upper image shows all DSC components, lower image is without shell to allow view of the weld surface.
Thepeak strain of 6.09% occurred on the surface of the shell. Therefore, when the shell was removed for thelower image, the peak strain reported reduced to 5.49%.