ML20134E159

From kanterella
Jump to navigation Jump to search
Forwards Brief Summary Rept on Leakage Characteristics of Nuclear Containment Hatches During Severe Accident. Draft of Complete Rept Will Be Mailed Later in Aug
ML20134E159
Person / Time
Site: Sequoyah, 05000000
Issue date: 08/02/1984
From: Greimann L
AMES LABORATORY, ENERGY & MINERAL RESOURCES RESEARCH
To: Bagchi G
Office of Nuclear Reactor Regulation
Shared Package
ML20134E148 List:
References
CON-FIN-A-4135-4, FOIA-85-457 NUDOCS 8508200054
Download: ML20134E159 (2)


Text

o, 1>..

. ta, . . s--

S5 '

_10BOROTORY Energy & Mineral Resources Research Institute Iowa State University l Ames, Iowa 50011 i August 2,1984 Mr. Goutam Bagchi, Leader Seismic Qualification Section -

Office of Nuclear Reactor Regulation U.S. Nuclear Regulatory Conmission Phillips Suilding, P-234 7920 Norfolk Avenue Bethesda, MD 20034

SUBJECT:

Summary Report for Ames Laboratory Technical Assistance to the Division of Engincaring, Nuclear Reactor Regulation--

" Leakage Characteristics of Nuclear Containment Hatches During a Severe Accident" FIN NO. A4135-4

Dear Gautam:

Enclosed, please find a brief summary report on the findings of our study of the Sequoyah equipment hatch. A draft of the complete report will be mailed later in August. .

If you have any questions, please call.

Sincerely, l(19 /LU4h& W )

Lowell Greimann Senior Engineer LG:cb Enclosure cc: Delwyn Bluhm Ames Laboratory Tom Bridges EGtG, Idaho, Inc. ,

iAnicles Hofmayer Brookhaven National Labs 8508200054 850724 PDR FDIA SHOLLY85-457 PDR n.-n: *o ~ u s w.sm .v t,..g. u s. cs , < w. ct t,; a2

9

  • ' 1 ..,

SEQUOYAH EQUIPMENT HATCH SEAL LEAKAGE L. Greimann, F. Fanous, D. Bluhm July 1984

SUMMARY

The following summary constitutes an abstract of a soon-to-be-published report by the same title in which complete details of the analysis are presented.

Background and Objective As part of the Containment Performance Working Group study, an analyti-cal effort was undertaken to investigate the leakage characteristics of the Sequoyah equipment hatch seal during severe accidents. The effects of containment deformation and hatch postbuckling on the distortions at the sealing surface were to be included. Only leakage at the seal surface is considered in this work. Leakage through a crack in the steel plate was not considered.

Equipment 'iatch Model The equipment hatch is one of several penetrations in the Sequayah containment vessel. If penetration effects are neglected, gross yield-ing of the 1/2-inch shell pressure of 50 to 60 psig.A plate near the springline three-dimensional will occur finite element modelat of a the equipment hatch was developed which includes: shell elements for the hatch, sleeve, containment plate and stiffeners; prestressed bar elements for the swing bolts; and friction elements for the seal interface. One-quarter symmetry of.the assembly was enforced.

Geometric and material nonlinearities were included. A single lobe .

imperfection was included in the hatch shell model. The imperfection amplitude was equal to the ASME tolerance allowable and the wave length was equal to the elastic buckling length. Loading which corresponds to pressure inside the containment was applied in increments. Analyses were made with two coefficients of friction for the steel at the ' seal interface: 0.3 and 0.6.

Results The maximum strain of about 1.4 percent occurred in the top of the '

sleeve, near its intersection with the containment. The maximum relative sliding and rotation at the seal interf ace was 0.9 inch and 1 1/2 degrees, respectively, for the 0.3 friction coefficient. Hatch buckling occurred in the range of 85 to 90 psi.

Conclusions The Sequoyah equipment hatch seals should not leak before very large strains develop in the 1/2-inch containment shell plate near the springline, which occurs between 50 and 60 psig. In the unlikely event of hatch buckling, postbuckling deformations would not introduce leakage.

1 E s

O AI i . . -

' NUREG/CR-3952 15-4862

- 7%

NETB Igr FLc1LEL3S!ON PURPUSES ONLY l

r i

, SEQUOYAH EQUIPMENT HATCH SEAL LEAKAGE l

Final Report \

August 1984 l

. 1 i

Prepared by l

L. Greimann, F. Fanous, D. Bluhm l Ames Laboratory

,. Iowa State University Ames, IA 50011 Prepared for Division of Engineering Office of Nuclear Reactor Regulation l ' U.S. Nuclear Regulatory Commission l

Washington, D.C. 20555

. NRC FIN NO. A4135-4 e m e M *w _

  • f _q u v W u p p v_w g j ig

NUREG/CR-3952

' 15-4862 ann mr FLa L;SCUSSION PURPJSES ONLY SEQUOYAH EQUIPMENT HATCH SEAL LEAKAGE Final Report August 1984 Prepared by L. Greimann, F. Fanous, D. Bluhm knes Laboratory Iowa State University Ames, IA 50011 l l

l l

Prepared for Division of Engineering Office of Nuclear Reactor Regulation U.S. Nuclear Regulatory Comission Washington, D.C. 20555 NRC FIN NO. A4135-4

.. t~._.

i 6e

. . _ , _ -e,- , . , , - - , , " - , , -+-- -,--

l

. DISCLAIMER This book ves prepared as an account of work sponsored by an agency of the United States Government. Neither the United States Government nor any agency thereof, nor any of their employees' makes any war canty, express or implied, or as-sumes any legal liability or responsibility for the accuracy, completeness or usefulness of any information, appar atus, pro-duct, or process disclosed, or represents that its use would not infringe pr !vately owned rights. Reference herein to any specific commercial product, process, or service by trade name , trade-mark, manufacturer, or otherwise , does not necessarily con-stitute or imply its endorsement, recommendation, or favoring by the United States Government or any agency thereof. The views and opinions of authors expressed herein do not neces-sarily stata or reflect those of the United States Government or any agency thereof. .

Printed in the United States of America Available from Natter.at Technical information Service U.S. Department of Commerce 5265 Port Royal Road Springfteld, VA 22161 l

e.

-. ., . _ _ . -,_ _ -. . - . - - - - - ~-,r.

1

~

. iii TABLE OF CONTENTS LIST OF FIGURES . . . . . . . . . . . . . . . . . . . . . . . . . iv ACKNOWLEDGMENT ......................... vi EXECUTIVE

SUMMARY

. . . . . . . . . . . . . . . . . . . . . . . . 1

1. INTRODUCTION ........................ 2 1.1 Background ....................... 2 1.2 Objective ....................... 2
2. EQUIPMENT HATCH DESCRIPTION . . ... . . . . . . . . . . . . . 4 2.1 Hatch Sel ecti on . . . . . . . . . . . . . . . . . . . . . 4 2.2 Geometry and Materials ................. 4
3. T IMPERFECTION S'UDY ..................... 6
4. THREE-0IMENSIONAL ANALYSIS ................. 7 4.1 Model . . . . . . . . . . . . . . . . . . . . . . . . . . 7
4. 2 Lo ad s an d Ex ec u t i o n . . . . . . . . . . . . . . . . . . . 8 4.3 Results . . . . . . . . . . . . . . . . . . . . . . . . . 9 4.3.1 Prebuckling Results . . . . . . . . . . . . . . . . 9 4.3.2 Postbuckling Results ............... 10 .

4.4 Summary Observations .................. 10

5.

SUMMARY

. . . . . . . . . . . . . . . . . . . . . . . . . . . 11 5.1 Summary ........................ 11 5.2 Conclusions ...................... 12 5.3 Recommendations .................... 12

6. REFERENCES ......................... 13 4

I

  • e p M*mpi

\

iv LIST OF FIGURES Figure 2.1 Sequoyah Containment - Azimuth 285* . . . . . . . . . 15 Figure 2.2 Sequoyah Containment Shell Plate Thicknesses. . . . . 16 Figure 2.3 Sequoyah Equipment Hatch Elevation ......... 17 Figure 2.4 Detail A - Sequoyah Equipment Hatch . . . . . . . . . 17 Figure 2.5 ASME Tolerance for Shells . . . . . . . . . . . . . . 18 Figure 2.6 Seal Configuration ................. 19 Figure 2.7 Idealized Stress Strain. Curve for Steel . . . . . . . 20 Figure 3.1 BOSOR5 Axisymmetric Finite Difference Model . . . . . 21 Figure 3.2 Idealized Imperfection at Crown . . . . . . . . . . . 22 Figure 3.3 Idealized Imperfection at Flange .......... 22 Figure 3.4 Influence of Imperfection Wavelength ........ 23 Figure 3.5 Influence of Imperfection Magnitude . . . . . . . . . 23 Figure 3.6 Assumed Imperfection for Three-Dimensional Analysis . 24 .

Figure 4.1 Sequoyah Containment Equipment Hatch - Outside View . 25 Figure 4.2 Sequoyah Containment Equipment Hatch - Inside View . 26 Figure 4.3 Flange / Sleeve Interface Idealization . . . . . ... . 27 Figure 4.4 Sequoyah Containment Shell Reference Points . . . . . 28 (see also Figures 4.7 and 4.8)

Figure 4.5 Rad i al Di s pl ac ement . . . . . . . . . . . . . . . . . 29 Figure 4.6 Strains at Selected Location in the Hatch Model . . . 30 Figure 4.7 Displace Shape in Vertical Symetry Plane . . . . . . 31 (82 psig, coefficient of friction = 0.3)

Figure 4.8 Displace Shape in Horizontal Symmetry Plane . . . . . 32 (82.psig, coefficient of friction = 0.3)

Figure 4.9 Flange / Sleeve Mismatch (82 psig, coefficient .... 33  !'

of friction = 0.3) .,

l Figure 4.10 Radial Displacement a': Seal Surf ace, Vertical . . . . 34 <

Symmetry Plane, F _,,

, l

, i .

  • v Figure 4.11 Radial Displacement at Seal Surface, Horizontal . . . 35 i Symmetry Plane, G Figure 4.12 Rotation at Seal Surface, Vertical Symmetry Plane . . 36 Figure 4.13 Rotation at Seal Surf ar.e, Horizontal Symmetry . . . . 37 Plane Figure 4.14 Displaced Shape during the Postbuckling process . . . 38 at Vertical Synmetry Plane (90 psig, coefficient of friction = 0.3) i h

- , , ~ _ . ~

me l i

e

.i i

vi l-ACKNOWLEDGMENT The authors would like to express their appreciation for three members  !

of the U.S. Nuclear Regulatory Commission, Mr. Goutam Bagchi. Leader, Seismic Qualification Section; Mr. Vince Noonan, Chief, Equipment Qualifications Branch; and, Mr. Robert Wright, NRC Project Manager, for their help throughout the course of this work. Mr. Charles Hofmayer of Brookhaven National Laboratory and Mr. Tom Bridges of EG&G, Idaho, t

Inc., supplied information which was very useful for the conduct of i this work. The authors also wish to acknowledge the able assistance of the project secretaries, Connie Bates and Beth Lott, for the word processor operations and the secretarial services associated with this project.

9 1

1 e o e $

  • pm , *

,i .

EXECUTIVE

SUMMARY

Nuclear containments, which will not leak during a design accident, may leak during a severe accident when the pressure increases beyond the design level. Small leaks which are visualized as occurring at local details may occur before complete vessel failure. As part of the NRC Containment Performance Working program, this study was undertaken to investigate the leakage-before-break potential of a typical equipment hatch seal. Buckling of the hatch door, large deformations and ovaling of the hatch sleeve are potential causes of mismatch at the sealing surface which can result in a leakage path.

Among the several plants being considered by the Containment Performance Working Group, the Sequoyah equipment hatch was selected.

If penetrations effects are neglected, gross yielding of the'1/2-inch shell plate near the springline of the Sequoyah containment will occur at an internal pressure of between 50 and 60 psi.

The effect of an imperfection on the buckling strength of the hatch was first studied using an axisymmetric model of the hatch, containment and sleeve. A single axisymetric imperfection lobe (cosine curve) was placed at the crown and adjacent to the seal. The buckling load for various values of the imperfection wave icngth and amplitutde was calculated using the BOSORS code. The nonlinear behavior of the material was approximated by a piecewise linear relationship which accounted for residual stress effects. The critical value of the imperfection wave length, amplitude and location was selected to be used in a three-dimensional finite element model. .

The ANSYS finite element computer code was used to perform a three-dimensional finite element analysis of the containment / sleeve / hatch assembly. One-quarter symmetry is assumed about the hatch vertical and horizontal diameters. The hatch system was modeled using flat triangular shell elements for the door, sleeve, containment, ring stiffeners and stringers. Prestressed bar elements were used to model the swing bolts and friction elements were employed to idealize the seal interf ace. Geometric and material nonlinear behavior were also included. The model was analyzed twice, using coefficients of friction of 0.6 and 0.3. In both analyses the pretension forces in the swing bolts were first applied and then the internal pressure was increased in increments as convergence was insured.

The results of the finite element analysis showed that a maximtan of 0.9 4 inch of relative sliding occurred at the seal interface at 82 psi. The maximum relative rotation of the sealing surfaces was 1-1/2 degrees.

The buckling load was predicted to occur in the range of 85 to 90 psi, but postbuckling displacements did not cause large seal motions.

In conclusion, the Sequoyah egaipment hatch should not leak before very large strains develop in the 1/2-inch containment shell plate near the -

springline, which occurs between 50 and 60 osi. In the unlikely event-of hatch buckling, postbuckling deformations would not introduce leakage. __.4_. -

. 2

1. INTRODUCTION

1.1 Background

Nuclear containments are designed to prevent leakage of radioactive material . As the pressure and temperature increase during a severe accident, which is beyond the design level, the containment may reach a point where leakage begins to occur. There are at least two models used to characterize this leakage--the threshold model and the leak-before-break model.

In the threshold model, the containment is assured to be leak-tight until certain pressure / temperature conditions are reached. At this point, a very.large rupture or burst area is postulated and large quantitites of fission products are released. Typically, the predicted threshold pressure is based upon a model of the containment shell which only permits such gross failure modes. Local discontinuities, e.g., ,

penetrations, hatches, seals, are usually omitted from such models

[1.1]. l As an alternative model, the leak-before-break model is probably more )

realistic. In this model, it is hypothesized that small leak paths will develop as the containment is being pressurized at levels below the threshold. Hence, pressure is being released at earlier stages and l the threshold pressure may never be reached. Typically, these small i leaks are visualized as occurring at local details, e.g., penetrations, l hatches, and seals which are often omitted from the threshold model.

The NRC has established a Containment Performance Working Group (CPWG), l at the request of the Severe Accident Research Plan (SARP) Senior

.I Review Group, to develop containment leakage Inodels. Besides the containment loads aspect of this problem, this group is studying containment leakage in several categories. The current hypothesis of the group is that leakage will, most likely, begin at certain penetrations or other discontinuities (leak-before-break) and not by gross rupture of the containment shell. As a test for this hypothesis, the group is currently involved in quantifying predictions of the amount of leakage which may occur in six representative containments.

1.2 Objective i

The objective of this particular work is to investigate the leak-before-break potential of a typical equipment hatch seal. The results will be incorporated into the work of the entire Containment Performance Working Group, i In the initial stages of the work, postbuckling and other large I displacements of the equipment hatch shell itself were considered to be l possible causes of leakage as local deformations distorted the sealing

. surface. As the project and comunication with others involved in _

similar work developed, another possibility became evident--ovaling of the penetration sleeve. That is, as the containment shell is so m

3 pressurized, circumferential stresses in the containment tend to distort the hole and sleeve for the equipment penetration into an elliptical shape with a horizontal major axis. This distortion causes a mismatch of the sealing surf ace with ' leakage potential . (This behavior is documented and described more fully herein.)

d

.*~

I I

1 -r i

4

- 4 l

2. EQUIPMENT HATCH DESCRIPTION 2.1 Hatch Selection There are many different equipment hatch configurations in existing l containments. Si.'.ce an initial consideration was postbuckling displacement, a hatch under compression was selected (pressure seating) as opposed to a hatch in tension (pressure unseating). Again, since
buckling was important, it was decided to select a hatch with the most likely possibility of buckling, i.e., a large r/t value. From the six plants being studied by the Containment Performance Working , Group, the Sequoyah equipment hatch was selected.

2.2 Geometry and Materials The Sequoyah equipment hatch is located on Azimuth 285* at Elevation 741' 1 1/2" of the containment (Fig. 2.1). More geometric details of the containment in the vicinity of th9 hatch are shown in Figs. 2.2, 2.3 and 2.4. Of particular interest in these figures are:

(1) The containment shell is 1/2 inch thick just below the springline. Previous studies [2.1] using axisymmetric (threshold type) models of the Sequoyah containment which neglect penetrations and other nonsymmetric discontinuities, predict shell failure in the 1/2 inch plate at a static pressure between 50 and 60 psig. (The variation in this value is the result of different approaches by the different investigations.) .

(2) There is 5/8 inch shell plate ia the imediate vicinity of the equipment hatch penetration (Fig. 2.2). The local plate thickness of 1-1/2 inch is presumably used to satisfy the ASME area replacement rule.

(3) The sleeve length between the containment wall and the sealing surfaces is only 14.5 inches at the horizontal diameter of the sleeve. This distance is quite short and, as will be seen, couples the containment and hatch shell deformations.

~

(4) The actual containment it the upper right quadrant of the hatch region is shown in Fig. 2.2, except a ring actually located at Elevation 740' 6 5/8" was moved to the hatch centerline to introduce symmetry. The containment in the remaining three quadrants is not symetrical, e.g., the 5/8 inch plate does not exist in the same locations in the other quadrants.

Shell geometric imperfections have a significant effect on buckling strength. Particularly for spherical shells under external pressure, small deviations from perfect sphericity greatly reduce the buckling strength below classical buckling theory. Figure 2.5 illustrates the w ASME tolerance criteria on shell imperfections [2.2]. For the Sequoyah hatch, the maximum deviation from a perfect circular arc, e, is m.

5 permitted to be one shell thickness (3/4 inch) over a chord length of 67 inches.

The details of the seals as furnished by TVA are shown in Fig. 2.6.

Compression set values (ASTM D 395) for typical temperature / radiation loadings are also listed.

The containment steel is A516, Grade 60 with a specified minimum yield of 32 ksi. Actual mill test yield strengths for typical plates with the various thicknesses are listed in Fig. 2.7. Residual stresses have a significant affect on the buckling strength of steel components.

When regions with residual stress yield, the material tangent modulus is reduced and the structural tangent stiffness is, thereby, reduced.

In lieu of knowing the actual initial residual stress pattern, this effect can be approximated by using an effective stress-strain curve with a proportional limit below the yield strength [2.4, 2.5]. For the present case, the proportional limit will be taken as one-half yield which approximates residual stresses of the order of one-half yield.

For each thickness, the stress-strain curve is approximated by a modified Ramberg-Osgood equation [2.2] with a strain hardening exponent n of five (Fig. 2.7). For purposes of the following analysis, this curve is approximated by a piecewise linear curve.

l l

>@ hyes 4

em

.i -

6

3. IMPERFECTION STUDY As mentioned previously, imperfections in the spherical cap hatch can have a large effect on its buckling behavior. To study this effect in an economical manner, an axisymmetric model of the hatch and sleeve was developedfor8050R5[3.1]. Since the model is axisymmetric, the cylindrical containment was approximated by a portion of a sphere (Fig.

3.1). No relative motion was permitted at the sealing surface.

A single axisymetric imperfection lobe (cosine curve) was placed at the crown (Fig. 3.2) and adjacent to the flange (Fig. 3.3). 8050R5 was executed to calculate the buckling load for various valves of the imperfection wave length, L o

, and amplitude. The resulting buckling pressures are plotted in Figs. 3.4 and 3.~5, respectively. The critical buckling length for classical elastic buckling of a perfect sphere is approximately [3.2]

(3.1)

L cr - 3.5 /iT = 47 in.

Figure 3.4 illustrates that this is approximately the wave length of the crown imperfection which minimizes the buckling strength of the Sequoyah hatch. The length of the (axisymetric) imperfection adjacent to the flange has little effect on buckling strength. As expected, Fig. 3.5 confirms that an increasing imperfection amplitude decreases buckling strength. Note again that the maximum permitted value of e for the Sequoyah hatch is 3/4 inch (see Sec. 2.3).

The " design" buckling pressure for the hatch as calculated from the ASME Code, Case N-284 [3.3] is 50 psi. (In terms of Case N-284, the capacity reduction factor is 0.146 (M = 240//(140)(0.75) = 17.9), the classical elastic buckling pressure is 343 psig and the plasticity reduction f actor is 1.0, i.e., elastic buckling.) This value is shown in Fig. 3.5. Note that a large imperfection is necessary to reduce the buckling load to that of the design equations.

For the subsequent three-dimensional analysis, a single imperfection lobe was used with a length, L o , of 47 inches and a deviation for perfect circularity, e, of 3/4 inch (Fig. 3.6). The imperfection was also assumed to be L in the circumferential direction, on the basis thatthisshapewillOesympatheticwiththeelasticbucklingmodein-which n lobes form around the circumference. (8050R5 indicated n was about 8 for this case [3.4].) The imperfection was located near the flange to maximize its effect on possible distortion at the sealing surface. In anticipation of the possible ovaling effect (Sec.1.2),

compression in the hatch will be largest in the vertical direction.

Hence, the imperfection was placed on a vertical diameter.

l 1

e 6 y T q.

l 7

4. THREE-DIMENSIONAL ANALYSIS 4.1 Model A three-dimensional finite element model of the Sequoyah equipment hatch was constructed for the ANSYS [4.1] computer program. ANSYS is a general purpose finite element computer program with extensive pre- and post-processing capability. The geometric and material nonlinear analysis capabilities of ANSYS were used in this model. The geometric and material properties presented in Sec. 2.2 and the imperfection in Fig. 3.6 were incorporated into the model. One-quarter. symmetry is assumed about the hatch vertical and horizontal diameters, as suggested in Fig. 2.2, even though the containment is only approximately
  • symmetrical.

Symmetry boundary conditions were used in the plane of the vertical and horizontal diameters of the hatch. The model is extended only 20 degrees circumferentially from the meridional plane tnrough the center of the hatch (Fig. 2.2). (This certainty represents a minimum.)

Symetry boundary conditions are also used on this boundary. The top of the model is constrained to move uniformly vertically, i.e., the top of the model remains a plane section'.

The mesh in the hatch and sleeve was layed out such that the maximum element dimension was less than 6t/2 (6.7 inch for the hatch and 9.5 inch for the sleeve). Flat triangular shell elements were used (STIF48 in ANSYS). The elements in the containment shell are noticeably larger than this guideline (frf/2 of 16 inch), but sufficiently small to -

approximate its membrane behavior. Rings and stringers were also modeled tiy shell elements. Figure 4.1 and 4.2 are plots of the finite element mesh from outside and inside the containment, respectively (915 elements, 541 nodes).

The swing bolts which are pretensioned to hold the hatch against the sealing surface (Fig. 2.4) were included in the model as uniaxial, pretensioned bar elements (STIF8). The swing bolts are rigidly connected to the middle surface of the sleeve shell.

The material model is the von Mises yield surface with the Prandtl-Reuss flow rule with isotropic strain hardening. The piecewise linear approximation to the uniaxial stress-strain curve is illustrated in Fig. 2.7.

The interface between the flange and the sleeve is modelled using interface elements (STIF 52) which have friction and opening capability l

(Fig. 4.3). Normal to the interface, the element is very stiff in compression and very flexible in tension. (Infinite and zero stiffnesses introduce numerical problems.) Tangential to the interface, the element permits relative sliding if the tangential force exceeds the normal (compressive) force multiplied by the coefficient of -

friction. Two interf ace elements are used at each nadal meridian of the sleeve--one on the inside and one on the outside of the sleeve.

8 Such an arrangement permits rotational and sliding discontinuities at the seal interface. (The results to follow illustrate this.) Two different coefficients of friction, 0.6 and 0.3, were used in the analysis of the hatch model. A coefficient of friction of 0.6 represents the clean steel on steel case [4.2].

4.2 Loads and Execution The finite element model described in the previous section was loaded to correspond to internal pressure on the containment, i.e., pressure on the inside face of the appropriate elements and a vertical load at the top of the model corresponding to the meridional membrane stress resultant in the containment. Two different analysis approaches available in the ANSYS Code were used-to analyze the hatch model. A static analysis was performed with.the 0.6 coefficient of friction case, while a " slow-dynamic" analysis was used in conjunction with the 0.3 friction coefficient. In both approaches, the prestress in the suing bolts was first applied (25 kips / bolt which is about 1/4 the yield strength) and then the pressure inside the containment was applied.

In the static analysis, the pressure was increased in increments, typically 5 psi. After each increment in pressure was applied, a sufficient number of iterations was run to permit convergence to the nonlinear solution. Convergence criteria were: (1) the change in displacement between consecutive iterations must be less than 0.1 inch and (2) the change in plastic strain divided by the yield strain be less than 0.1. Beyond 80 psi, these criteria were reduced to 0.05 and -

0.05, respectively. Through the 82 psi load level, the solution required iterations which used about 55 cpu hours on a VAX 11/780.

When the pressure was incree. sed to 84 psi, a zero-stiffness condition occurred, indicating the hatch had buckled at the imperfection. A couple unsuccessful attempts were made to follow the postbuckling behavior using the traditional static analysis.

In an attempt to track the postbuckling behavior with a physically, more realistic analysis, the " slow dynamic" analysis suggested in [4.1]

was used for the 0.3 coefficient of friction case. In principle the method is more attractive because the actual buckling process is dynamic. A high damping ratio is used to minimize the vibration

. - response. Physically, this would correspond to placing the hatch in a

! viscous fluid during pressurization. The pressure in this analysis was l increased linearly over a rise time larger than several times the structure natural period (about 1 psi per second). The integration time increnents were initially quite large (1 second) as the structure behaved linearly and the high damping produced an essentially static

, response. As material and geometric nonlinearities began to dominate

! the problem, the time steps were successively cut in half. For l example, a time step of 0.125 seconds was used at 82 psi. At this point, the pressure was increased to 85 psi as a step function and 24 -

time increments of 0.01 second were run. Since r.o buckling occurred and since the postbuckling shape is important (not ndcessarily the buckling pressure), the pressure was stepped to 90 psi and the solution -- +_-

was performed for 100 time steps. During these steps, the time -

9 increment size varied between 0.01 to 0.1 seconds was used. The authors recognize that the larger time step size would yield an 1 inaccurate result, however, it was considered acceptable since the interest was to investigate the postbuckling behavior not to reach a converged solution.

' Initially, it was hoped that the time step optimizer within ANSYS may reduce the number of steps and, hence, the amount of cpu time.

However, to 82 psi, the " slow-dynamic" required about 70 cpu hours time than the static analysis.

4.3 Results 4.3.1 Prebuckling Results ,

Typical strain and displacement results are presented in this section.

Figure 4.4 will be used to reference the location of the points of interest. Figure 4.5 is a plot of the radial displacement at the center of the hatch, A, the center of the imperfection, B, and the 5/8-inch containment shell plate, C. The displacement at B is the maximum in the model. Strains of interest are plotted in Fig. 4.6. The maximum strain in the hatch occurs at the surface in the imperfection, B. The maximum membrane strain in the containment occurs in the 5/8 inch plate near C. The maximum membrane strain in the entire model occurs in the horizontal direction at the top of the sleeve, D, just outside the containment. Most of the differences between the results for the 0.3 and 0.6 friction cases are caused by the difference in the analysis techniques. -

The distortions at the sealing surface can be visualized by examining the deformed shapes at the vertical and horizontal symmetry planes of the hatch, Fig. 4.7 and 4.8, respectively. These figures are plotted for the 0.3 coefficient of friction at 82 psig. Displacements are drawn at an exaggerated scale in these figures. The inset figures in each of these plots are drawn to scale and illustrate the translation and rotation relative to the flange and the flange grooves. These figures and Fig. 4.9 illustrate the ovaling effect, i.e., the forces in the containment shell tend to distort the sleeve into an ellipse. If the frictional forces at the sealing surface where zero, the sleeve would become elliptical and the hatch flange would remain circular with a uniform outward radial motion. This mismatch in shapes causes sliding and rotation at the sealing surface. Friction reduces this mismatch, i.e., tends to force the two surfaces to move together, but it does not eliminate the mismatch.

Figures 4.10 and 4.11 are plots of the vertical displacement in the vertical symmetry plane, F, and the horizontal displacement in the horizontal synnetry plane, G, respectively, at the sealing surface.

Down is positive in Fig. 4.10 and to the right is positive in Fig.

4.11. The difference in the displacement for the sleeve and the flange -

represent the slip at the seal surface. Notice that, as expected, the slip is larger when the coefficient of friction was reduced to 0.3.

Figures 4.9 and 4.10 present the rotations at the sealing surface in _.. __.

the vertical and horizontal symmetry planes, respectively. Positive rotations are shown in the inset. Again, the distance between the two

10 i

curves represents the relative rotation of the two surfaces. The relative rotation was not much affected by the change in the coefficient of friction. )

4.3.2 Postbuckling Results As previously mentioned, the solution was continued to track the postbuckling behavior of the hatch, especially the seal surface.

During the postbuckling process, the radial displacement in the imperfection region continued to grow rapidly compared to the rest of the hatch model. The displacement at the center of the hatch, center of imperfection and the 5/8 inch containment plate are 8.2 inches, 23 inches, and 48 inches, respectively. Despite these large displace-ments, no noticeable changes in the pre- and postbuckling relative

-rotatien and translation were found. The maximum displacement of 27.8 inches occurred along the vertical symmetry plane in the vicinity of the imperfection region. Figure 4.14 illustrates the unexaggerated defurmed shape of the model. Imposed also on Fig. 4.14 is an inset figure drawn to scale illustrating the relative translation and rotation of the flange and flange grooves. The maximum values of the relative sliding and rotation reached during the postbuckling investigation were 1.12 inches and 1.80 degrees, respectively.

4.4 Summary Observations The maximun strain in the model .at 82 psig, for both friction cases, is about 0.014 in the penetration sleeve (Fig. 4.4). Although this strain is probably well below the ultimate strain for the ductile steel -

) used in the containment, the possibility exists that localized plane strain conditions and local flaws may initiate a leakage point.

Elastic-plastic f acture mechanics principles could usefully be applied here[4.3].

The maximum relative sliding at the seal surface are:

Coefficient of Relative Relative Rotation Friction Sliding (in.) (Degrees) 0.6 0.6 1 1/2 0.3 0.9 1 1/2 With these small amount of displacement of the seals, leakage is very unlikely [4.4]. Even postbuckling displacements did not introduce  :

large displacements at the seals (Sec. 4.3.2).  ;

  1. @ p gges em

,, ,,-. - _s.__.,. . . - _ . .~4,._ , .

l.. -

l 11

5. SUMARY 5.1 Summary Nuclear containments, which will not leak during a design accident, may  ;

leak during a severe accident when the pressure increases beyond the l design level. Small leaks which are visualized as occurring at local details may occur before complete vessel failure. As part of the NRC Containment Performance Working program, this study was undertaken to investigate the leakage-before-break potential of a typical equipment hatch seal. Buckling of the hatch door, large deformations and ovaling of the hatch sleeve are potential causes of mismatch at the sealing surface which can result in a leakage path.

Among the several plants being considered by the Containment Performance Working Group, the Sequoyah equipment hatch was selected.

If penetrations effects are neglected, gross yielding of the 1/2-inch shell plate near the springline of the Sequoyah containment will occur at an internal pressure of between 50 and 60 psi.

The effect of an imperfection on the buckling strength of the hatch was first studied using an axisymmetric model of the hatch, containment and sleeve. A single axisymetric imperfection lobe (cosine curve) was placed at the crown and adjacent to the seal. The buckling load for various values of the imperfection wave length and amplitutde was calculated using the 8050R5 code. The nonlinear behavior of the material was approximated by a piecewise linear relationship which ccounted for residual stress effects. The critical value of the .

, imperfection wave length, amplitude and location was selected to be used in a three-dimensional finite element model, i The ANSYS finite element computer code was used to perform a three-dimensional finite element analysis of the containment / sleeve / hatch assembly. One-quarter symetry is assumed about the hatch vertical and horizontal diameters. The hatch system was modeled using flat triangular shell elements for the door, sleeve, containment, ring stiffeners and stringers. Prestressed bar elements were used to model the swing bolts and friction elements were employed to idealize the seal interface. Geometric and material nonlinear behavior were also included. The model was analyzed twice, using coefficients of friction of 0.6 and 0.3. In both analyses the pretension forces in the swing bolts were first applied and then the internal pressure was increased in increments as convergence was insured.

The results of the finite element analysis showed that a maxima of 0.9 inch of relative sliding occurred at the seal interface at 82 psi. The maxima relative rotation of the sealing surfaces was 1-1/2 degrees.

The buckling load was predicted to occur in the range of 85 to 90 psi, but postbuckling , displacements did not cause large seal motions.

~

l d

e. e a T *w +

12 5.2 Conclusions In conclusion, the Sequoyah equipment hatch should not leak before very large strains develop in the 1/2-inch containment shell plate near the springline, which occurs between 50 and 60 psi. In the unlikely event of hatch buckling, postbuckling deformations would not introduce leakage.

5.3 Recommendations Some of the items which were not included in this study which could be looked at in the future include other hatch imperfections, different friction coefficients and other pretensions in the swing bolts. The analysis methods should continue to be calibrated with existing experimental results, many of which have been published, e.g., shell buckling, pressure vessel experiments. The behavior and the seals with flange separation and rotation needs more study. The effect and temperature gradients on displacements at the seal surface could be important.

High strains in the sleeve, near the containment, should be investigated further. As pressure is increased, leakage will most likely begin at these high strain locations. The potential size of such a path should be studied. Several other penetrations with possibly higher strains should also be considered.

e b

w a gy f

l

13

6. REFERENCES 1.1 Greimann, L, Fanous, F., and Bluhm, D., " Containment Analysis Techniques, A State-of-the-Art Summary," Report to Sandia National Laboratory, October 1983.

2.1.Greimann, L., et. al, " Reliability Analysis of Containment Strength, Sequoyah and McGuire Ice Condenser Containments,"

NUREG/CR-1891, August 1982.

2.2 ASME Boiler and Pressure Vessel Code, Division III.

2.3 Desai, C.S. and Wu, T.H., "A General Function for Stress-Strain Curves," Numerical Methods in Geometrics, ASCE,1976.

2.4 Brockenbrough, R.L'. and Johnston, B.G., Steel Design Manual, U.S.

Steel, Pittsburgh, PA,1974, pg. 47.

2.5 Salmon, C.G. and Johnson, J.E., Steel Structures Design and Behavior, 2nd Edition, Harper and Row, New York,1980, pg. 263.

3.1 Bushnell, D., "B050RS, A Program for Buckling of Elastic-plastic Complex Shells of Revolution Including large Deflections and Creep," Structural Mechanics Laboratory, Lockheed Missile and Space Co., Palo Alto, CA,1974.

3.2 Bruch, 0.0. and Almarth, B.0., Buckling of Bars, Plates, and Shells, McGraw-Hill, New York,1975, Chap. 7. ,

3.3 ASME Boiler and Pressure Vessel Code, Case N-284, " Metal containment shell Buckling Design Methods," Supplement #2 to Nuclear Code Case Book, 1980.

3.4 Wackel, N., Jullien, J.F. and Ledermann, P., " Experimental Studies of the Instability of Cylindrical Shells with Initial Imperfections," Recent Advances in Nuclear Component Testing and Theoretical Studies in Buckling, PVP-Vol. 89, A5MF., June 17-21, M.

4.1 ANSYS Engineering Analysis System, User's Manual and Theoretical Manual, Swanson Analysis Systems, Inc., Houston, PA,1984 (Version 4.1C).

4.2 Handbook of Chemistry and Physics, 41 Edition, Chemical Rubber Publishing Co., Gleveland, 1960, pg. 2152.

4.3 Idelsohn, S.R., et.al., " Failure Internal Pressure of Spherical Steel Containments," Second Workshop on Containment Integrity, Washington, D.C., June 1984.

l ..

. 14 4.4 Barnes, B.L., " Containment Penetration Elastomer Seal Test,"

Second Workshop on Containment Integrity, Washington, D.C., June 1984.

m as @ >

me

l 15 t = IS/l6" t = 9/16" t e l/2" als' 4i-6 so9'sk" so3'g"*

799'9t in 7M O" wi-

.R 7ang-9 -

3 778'6{"

3 769'I" 8

T_

., 759'6["

.s

= 74'IE-p _

.i , 73o 33-

_ i, 721' sg"

  1. I

, 7:3' s-

- g' , roi q-

  • y . 691'2)"

soo95-Figure 2.1 Sequoyah Containment - Azimuth 285*

~

\

16 I

i i

STRINGERS AT 4' (9fx f) g ELE. 764' 31' ELE. 759' 6f ELE.753'6" s?

\ ELE.750'1F

= .

U SLEEVE i

, _ ELE. 741' 71 20' ARC Oe'AreTsicxuesscie.)

19t*x it' RINGS Figure 2.2 Sequoyah Containment Shell Plate Thicknesses ~

m.

17 l

l l

CONTAIPMENT SHELL , DETAIL A w (FIG. 2.4) l '

E I SLEEVE HATCH t R=240" i f =i I

3/4" s l

$ I .

l E.

. . u. . . .

ELEV. 741' 11/2" Figure 2.3 Sequoyah Equipment Hatch Elevation

  • I k" DIA. SOLTS CONTAINMENT -

~M _._

6*(Tee) 25'(TOP S BOTTOM) Sh IGF($10ES) 14.S"($10ES) N.T.S.

~

Figure 2.4 Detail A - Sequoyah Equipment Hatch '-

l 18 CHORD LENGTH g/'(Uh

, CIRCULAR TEMPLATE SHELL r

SURFACE CHORD LENGTH = TWICE THE ARC LENGTH ARC LENGTH = FIG. 4221.2.2 ASME CODE e= DEVIATION FROM TRUE SHAPE M9. 4221.2.1 ASME CODE r= INSILE OR OUTSIDE SHELL RADIUS Figure 2.5 A./E Tolerance for Shells w

.o 't " w a

- - _ - - , . - - - r ,-. - .

l 19 0.6 41" O.718 " .282" O.718" 0.641" H-

- SEALS 15/16" 54 PRESRAY EPDM COMPOUND E603 COMPRESSION SET 158'F FOR 22 HRS 14.1 %

5 7 x 10 RADS AT 120*F FOLLOWED 13.1%

BY 108RADS AT 250*F FOR 72 HRS

! Figure 2.6 Seal Configuration i ...,._s -

( .

20 i

f_ MODIFIED RAMBERG OSGOOO 77 e /s y e/a y = ;n=5

( l + l < /cy l" ) "

1.0 - --

,e

/

- PIECEWISE LINEAR

/

O.5 -

/

/O.2 % OFFSET

/

6 "Y -

l l

A 516 Gr. 60 ft. THICKNESS (in.) ey(KSI) 5" T 49.2 3"

y 43.7 1[ 45.5 1[ 46.1 l

~

3" 40.8

~ ~ ' ' ~ ~ ~

Figure 2.7 Idealized Stress-Strain curve for Steel

21 AXIS OF SYMMETRY HATCH SLEEVE CONTAINMENT (SPHERICAL APPROXIMATION)

Figure 3.1 B050R5 Axisymetric Finite Difference Model i

es a, e T

  • ep 6-

,e -- -- - , - , - - - - _ . - -g.-- - .-- - - - .- - . . , _ - . , -

9 4

Lo/E

%< /

\

3 s #

7%ufe 3.1 ided i**6 M9"#gect\00 gt (f0VO

  • v 4

s T -

Lo i\94f' 3.3 ideo S**d W9"#fect\0D flat 98 g 4

- - - , - m. - -- - w

23 150 i l l 1 I l

.. . . /

3 5 100 -

/. / . / __ -

o

$ IMPERFECTION AT CROWN -

w 50 -

[ --lMPERFECTION AT FLANGE I I I I I I 10 20 30 40 50 60 70 IMPERFECTK)N WAVELENGTH Lo (in.)

Figure 3.4 Influence of Imperfection Wavelength 200 IMPERFECTION AT CROWN -

IMPERFECTION AT 150 FLANGE E. .N,N.N w

m g 100 en

- N'N.N %,%

w E i ASME DESIGN 50 -

CASE N-284 I I O.5 1.0 1.5 1

DEVIATION FROM TRUE SHAPE e (in.) .. . . _ . .

Figure 3.5 Influence of Imperfection Magnitude ,,

i _.

24 O

'v' O

Lo ELEV.

741'l "

i I I Figure 3.6 Assumed Imperfection for Three-Dimensional Analysis l

l

. . . + .

Oe

)

. . .-.. -.. . _ _ - . - - . _ - _ _ _ _ _ _ _ _ _ _ _ _ _ _ . . _ - _ - .- __ -_ 1

- '- - - - + -- - -- e- - m- -u -- - _, -,_ _ _

e 4 e

\

7,6

\

l d 3

\ ,\

6 /Q kh 1, L

s4 a  %

\n

\

\

\

k

& \N /\

\

L \

\

t W ' g,g St@0I conto *## g,459# geg.Outd

  • 1\t*

$ g h-4,.

_ , _ _ , - . , , , _ , - - - - - - - --------,m- ,,,e-wym_ , , - - - - - --v-,,,ww-,n----, e-,------e-c,-,-----pg,- w,-,w w-n ' ~*~ ' ' ~ ' ' -="" ~M M ='w' " '

26 I

6  %

6 i 4 /- - - -.

/

8 .

Figure 4.2 Sequoyah Containment Equipment Hatch-Inside View

..,s_. -

j ma

l l

27 3"

F>8 n

,; Fs,s WRY SMALL MFLANGE

[ INTERFACE ELEMENT =  :

l l

SLEEVE y r

Fh Fs a u

.Fn- - - -

58 =s CONTACT SLIDING  !

l Figure 4.3 Flange / Sleeve Interface Idealization )

en e T we me

- , - -- ,-----w - - - - - _ - . - , . - - ,.

l 28 l ELE.764'31' ELE. 759' 6Y

)

ELE. 753' 6" l

/

ELE. 750' 1E

, Hatch SLEEVE h j i ELE.741' 7E 6

Figure 4.4 Sequoyah Containment Shell Reference Points (see also Figure 4.7 and 4.8) m 9 e me

29 100 , , , ,

Contai ent 5/8" .CenterofImperfection@

Plate C 80 - # '- - -

/ -

, sf, G //

A ,/

60 - Center of Hatch @ _

w /

m S COEFFICIENT m 40 -

OF FRICTION -

w z ----

O.3 O.6 .

20 -

O I I ' l O.00 1.00 2.00 300 4.00 DISPLACEMENT (IN.)

Figure 4.5 Radial Displacement nem e e 7 w.

.+

_ _ _ _ _ _ _ _ . - r

l l

30 I

100 l I l J.__ _

ImperfectiAn ' Sleeve, Top (Surface) UB (Membrane) D 80

\i N

2 / '

/

~

60 -

w / Containment 1 5/8" Plate g (Membrane)(C)

" COEFFICIENT

$ 40 -

OF FRICTION -

E -----

O.3 0.6 .

20 - -

0 I I ' I O.00 4.00 8.00 12.0 0 16.0 0 20.00 MICROSTRAIN 2 l (x 10 )

Figure 4.6 Strains at Selected Location in the hatch Model

  • O , T w.

ee

31

" Y

+

s v II Containment

=

@ F  ;

d Sleeve Hatch Figure 4.7 Displaced S'hape in Vertical Symetry Plane (82 psig, coefficient of friction = 0.3) l

. e , 7 ~w .

  • ge

- --- - w

32 I

l I

r bb

)

r Hatch Containment Sleeve j Figure 4.8 Displaced Shape in Horizontal Symetry Plane (82 psig, coefficient of friction = 0.3)

    • e ,

~

m.

33 Flange Sleeve Figure 4.9 Flange / Sleeve Mismatch (82 psig, coefficient of friction = 0.3) l i

96

l 34 100 , g i i Flange Sleeve Bo sb- 8 ,3- -

s

=

/

\

u) t /

1

~

60 -al I

/s' -

w /

m I /

3 I /

$ 4o -

,' / COEFFICIENT -

W r / oF FRICTION T /

n. / ----

o.3

/ ~

o.6 2o -

/ _

t j /

I I I I I o'

o.co o.zo o.4o o.6o o.so i.co DISP LACEMENT (IN.)

Figure 4.10 Radial Displacement at Se Surface.

Vertical Symetry Plane,

. . . ~ .

ee

- - - s-w-

l l l

\

35 I

i l

10 0 i i i i Flange Sleeve sb -

\

80 -

/

s',

/- ~

/

s e / /

& I /

S 60 w

I l /

/ ~

$ I [ COEFFICIENT m I m 40 -I /

/ OF FRICTION -

W l / ----

03

@ l j' O.6 .

20 - -

O ' I I I O.00 0.20 0.40 0.60 0.80 1.00 DISPL ACEME NT ' (IN.)

Figure 4.11 Radial Displacement at Seal Surface.

Horizontal Symetry, G ees l

I

.' i 1

~  !

36 l l

100 i i i i i Sleeve Flange 80 -

~~~ < -

\ #>

g 3

^ s*s\ I f

G a

60 -

s -

- 1 ', ,

1 LAl E \

o COEFFICIENT \ 6 l

40 - f

's OF FRICTION ]

y -----

O.3 ( ,

a O.6 g \'

20 -

g s

I s

l  %%

~~

O ' ' ' '

l -1.60 -l.20 -0.80 -0.40 0.00 0.40 0.80 i ROTATION (DEGREES) l l

Figure 4.12 Rotation at Seal Surface, Vertical Symetry Plane l

1 e O p we e-1

I

, l l

37 l

\ .

i 10 0 i i i i

%A -

~

\ _

,./

t E 60 - ik -

b Sleeve _

l a

j Flange m 40 -

t -

m i COEFFICIENT i y OF FRICTION .

o. _____ o,3 20 -

0.6

~

\

I I I O

-1.20 -0.80 -0.40 0.00 0.40 0.80 ROTATION (DEGREES)

Figure 4.13 Rotation at Seal Surface Horizontal 1 Symetry Plane e em , ya em e m e-w-, m

38 1

?

k E

5 q,

@ _j h 1.80' Figure 4.14 Displaced Shape during the Postbuckling process at Vertical Symetry Plane (90 psig, coefficient of friction = 0.3)

  • e p a h

. _ . _