ML16127A493: Difference between revisions
StriderTol (talk | contribs) (Created page by program invented by StriderTol) |
StriderTol (talk | contribs) (Created page by program invented by StriderTol) |
||
(5 intermediate revisions by the same user not shown) | |||
Line 14: | Line 14: | ||
| document type = Updated Final Safety Analysis Report (UFSAR) | | document type = Updated Final Safety Analysis Report (UFSAR) | ||
| page count = 28 | | page count = 28 | ||
}} | }} | ||
=Text= | |||
{{#Wiki_filter:TMI-1 UFSAR APPENDIX 5D 5D-i REV. 18, APRIL 2006 | |||
APPENDIX 5D REPORT ON STRESS ANALYSIS OF BUTTRESS | |||
TMI-1 UFSAR APPENDIX 5D 5D-ii REV. 18, APRIL 2006 TABLE OF CONTENTS | |||
SECTION TITLE 1 DESIGN DETAILS | |||
===1.1 TENDON=== | |||
ANCHORAGE HARDWARE | |||
===1.2 ANCHORAGE=== | |||
ZONE REINFORCEMENT | |||
2 TENDON LOADS AND ALLOWABLE STRESSES | |||
===2.1 TENDON=== | |||
LOADS | |||
===2.2 ALLOWABLE=== | |||
STRESSES | |||
3 ANCHOR BLOCK STRESSES | |||
===3.1 BEARING=== | |||
STRESSES | |||
===3.2 EQUIVALENT=== | |||
FLUID PRESSURE METHOD | |||
===3.3 TRANSVERSE=== | |||
TENSION-PRISM METHOD | |||
====3.3.1 REINFORCEMENT==== | |||
FOR TRANSVERSE FORCES | |||
===3.4 UNBALANCED=== | |||
FORCES | |||
====3.4.1 VERTICAL==== | |||
UNBALANCED FORCES | |||
====3.4.2 HORIZONTAL==== | |||
UNBALANCED FORCES | |||
====3.4.3 EQUILIBRIUM==== | |||
CONDITION IN THE ANCHOR BLOCK | |||
4 CHECK ANALYSIS BY FINITE-ELEMENT METHOD | |||
==4.1 INTRODUCTION== | |||
===4.2 PLANE=== | |||
STRESS FINITE-ELEMENT COMPUTER PROGRAM | |||
===4.3 GEOMETRY=== | |||
OF BUTTRESSES WITH TENDONS | |||
===4.4 THERMAL=== | |||
LOADS AND TEMPERATURE PROFILES 4.5 NON THERMAL LOADS | |||
===4.6 MATERIAL=== | |||
PROPERTIES | |||
===4.7 LOADING=== | |||
CASES AND COMBINATIONS | |||
===4.8 STRESS=== | |||
ANALYSIS | |||
===4.9 ANALYSIS=== | |||
OF RESULTS | |||
4.10 CONCLUSIONS | |||
5 REFERENCES | |||
TMI-1 UFSAR APPENDIX 5D 5D-iii REV. 18, APRIL 2006 LIST OF TABLES | |||
TABLE TITLE 5D-1 Buttress Anchor Zone Transverse Stresses in Radial and Vertical Planes | |||
5D-2 Buttress Anchor Zone Spalling | |||
TMI-1 UFSAR APPENDIX 5D 5D-iv REV. 18, APRIL 2006 LIST OF FIGURES | |||
FIGURE TITLE 5D-1 DELETED | |||
5D-2 DELETED | |||
5D-3 DELETED | |||
5D-4 DELETED | |||
5D-5 DELETED | |||
5D-6 DELETED | |||
5D-7 Typical Hoop Tendon Arrangement | |||
5D-8 Buttress-Hoop Tendon Bearing Stresses | |||
5D-9 Buttress-Radial Bursting Stresses | |||
5D-10 Buttress-Vertical Bursting Stresses | |||
5D-11 Hoop Tendon Typical Anchor Force Distribution | |||
5D-12 Hoop Prestress Idealized Vertical Stress Flow in Buttress | |||
5D-13 Idealized Model for Spalling | |||
5D-14 A Triangular Plane Stress Finite-Element | |||
5D-15 Reactor Building With Buttresses | |||
5D-16 Adiabatic Sections | |||
5D-17 Model of Gridwork for Thermal Analysis | |||
5D-18 Temperature Profiles Through Wall and Buttress for Various Thermal Load Conditions | |||
5D-19 Iso-Thermal Curves for Winter Normal Operating Temperature Condition | |||
5D-20 Iso-Thermal Curves for Winter A ccident Temperature Condition at Time = 10,000 Secs. | |||
5D-21 The Hoop Tendon Loads, Internal Pressure and the Equivalent Liner Pressure | |||
5D-22 Finite-Element Grid Work for Buttress Analysis TMI-1 UFSAR LIST OF FIGURES | |||
FIGURE TITLE APPENDIX 5D 5D-v REV. 18, APRIL 2006 5D-23 Minimum Principal Stresses for Maximum Prestress Forces and Winter Normal Operating Condition (Case A) (PSI) | |||
5D-24 Maximum Principal Stresses for Maximum Prestress Forces and Winter Normal Operating Condition (Case A) (PSI) | |||
5D-25 Minimum Principal Stresses for Maximum Prestress Forces and Winter Accident Condition (Case B) (PSI) | |||
5D-26 Maximum Principal Stresses for Maximum Prestress Forces and Winter Accident Condition (Case B) (PSI) | |||
5D-27 Minimum Principal Stresses for Minim um Prestress Forces and Winter Normal Operating Condition (Case C) (PSI) | |||
5D-28 Maximum Principal Stress for Minimum Prestress Condition (Case C) (PSI)Forces and Winter Normal Operating | |||
5D-29 Minimum Principal Stresses for Minim um Prestress Condition (Case D) (PSI) | |||
Forces and Winter Accident | |||
5D-30 Maximum Principal Stresses for Minimum Prestress Forces and Winter Accident Condition (Case D) (PSI) | |||
TMI-1 UFSAR APPENDIX 5D 5D-1 REV. 18, APRIL 2006 1 DESIGN DETAILS | |||
===1.1 TENDON=== | |||
ANCHORAGE HARDWARE | |||
The prestressing system employs the BBRV anchor age system utilizing 169-1/4" diameter wires conforming to ASTM A421-65, Type BA. The dimensions and capacities of the tendons used | |||
hereafter for the anchorage zone analysis are more fully described in Appendix 5B, Part D. | |||
===1.2 ANCHORAGE=== | |||
ZONE REINFORCEMENT | |||
The mild steel reinforcement in the tendon anchorage zones of the buttress and the base mat | |||
are shown on GAI Drawings 421004, 421005, 421006, 421031,421032, and 421033. This | |||
reinforcement was conservatively designed by the methods enumerated hereafter in Section 3 of this Appendix and is judged to be sufficient to accommodate all predicted stress states | |||
without jeopardizing the ability to satisfactorily place the concrete. | |||
2 TENDON LOADS AND ALLOWABLE STRESSES | |||
===2.1 TENDON=== | |||
LOADS | |||
The minimum guaranteed ultimate tensile strength of the tendons is approximately 2,000 kips. | |||
Before the tendons are anchored (locked off) at 70 percent of the ultimate tendon force, which is | |||
equal to 1,400 kips, they may be temporarily ov erstressed to 1,600 kips in order to reduce friction losses. The anchor zones are designed for anchor forces of 1,400 kips per tendon. The | |||
anchor zones are also checked to resist safely temporary anchor force of 1,600 kips per tendon. | |||
===2.2 ALLOWABLE=== | |||
STRESSES | |||
The allowable stresses for concrete and steel will conform to those given in Section 5.2.2 and | |||
5.2.3 of the FSAR. | |||
Due to the specialized nature of anchor zone design the following conditions and stress criteria | |||
will be used in the design. These stress criteria include adequate conservatism in their | |||
application. | |||
First, verify that the allowable concrete compression and tensile stresses are not exceeded due | |||
to bearing and bursting stresses respectively. | |||
Second, check the tensile stresses in reinforcing steel when no credit is taken for concrete | |||
tensile capacity in the investigated anchor zones. | |||
Under these conditions, the allowable stresses will be shown as follows: | |||
TMI-1 UFSAR APPENDIX 5D 5D-2 REV. 18, APRIL 2006 CONCRETE | |||
: 1. Allowable bearing stress | |||
according to ACI/318-63, 1 Paragraph 2605(C) | |||
: 2. Allowable compression stress in the structure beyond the anchorage zone c c f 0.45 f= | |||
: 3. Allowable Tension Stress | |||
REINFORCING STEEL | |||
Allowable tension for | |||
A615-Grade 40, f s = 20 ksi. | |||
3 ANCHOR BLOCK STRESSES The investigation of the buttress anchorage zones for the hoop tendons is described in the following sections. The analytical work as presented, is based on the end anchor zone stresses | |||
and stress flow pattern studies documented by Y. Guyon (Reference 2) and F. Leonhardt (Reference 3). Their work has been substant iated through extensive strain measurements and stress flow pattern established by studies utilizing photoelastic techniques. | |||
The design presented herein is quite conservative; nevertheless, it has been adopted as the | |||
basis for determination of required reinforcing steel in the end anchor zones to overcome tensile | |||
splitting stresses and to provide for spalling stresses near the concrete surface. | |||
The existance of a three-dimensional stress field including high shear, thermal and shrinkage | |||
tensile stresses is influenced by hairline cracking due to creep and thermal stresses. In order to | |||
accommodate the most unfavorable stress combination in the anchorage zones, the following | |||
design goals were established: | |||
: 1. Maximum concrete bursting tensile stresses have been kept below maximum allowable tensile stress of conservatively c f 3using Y. Guyon's symmetrical prism method. This would indicate that no special spiral and bursting reinforcement is required. The 3 b b ci cp/A A' f 6.0 f=ksi. 5.0 f ci=ksi. 5.0 f c=c f 3ksi. 5.0 f c= | |||
TMI-1 UFSAR APPENDIX 5D 5D-3 REV. 18, APRIL 2006 combined stresses during normal operation and test loading also indicate that the concrete will remain uncracked in the anchorage zone. | |||
: 2. However, due to the existance of complex thermal stresses and possible microcracking, the ability of concrete to carry tensile stresses has been discarded, and mild steel | |||
reinforcing has been provided for: | |||
: a. the equivalent fluid pressure | |||
: b. transverse tensile stresses | |||
: c. unbalanced vertical forces | |||
: d. unbalanced horizontal forces | |||
: e. equilibrium condition | |||
The combined stresses during the design pressure for the Reactor Building (p = 55 psig) indicate that the concrete will remain uncracked in the anchorage zone. The combined | |||
stresses during the factored accident loading (refer to 5.2.3.2, cases a. and b. indicate | |||
that the concrete in the anchorage zone may develop hairline cracks. | |||
: 3. Existing theoretical solutions for anchorage zone design do not agree with more recent experimental results (References 17 and 18); maximum transverse tensile forces were | |||
made up to 5.8 times larger than determined by Y. Guyon. | |||
: 4. Possible failure mechanisms due to spalling and bursting stresses were evaluated to verify static equilibrium and acceptable stre ss limits under such hypothetical conditions. | |||
: 5. Bursting stresses in a large body are much smaller than in a member bounded by planes (Reference 19). Therefore, while it was found that two-dimensional theories underestimate | |||
the bursting stresses (Reference 17), analytical results based on isolated prisms or better, experimental values such as those developed by J. Zielinski and R. E. Rowe can be used for massive structures (Reference 19). | |||
: 6. The design of the anchorage zones provides a margin of safety similar to the entire reactor building. Furthermore, the reinforcement is designed to control cracking due to | |||
spalling, bursting, or shear of the concrete under sustained load. | |||
===3.1 BEARING=== | |||
STRESSES | |||
Figures 5D-7 through 5D-10 show the typical hoop tendon arrangement and buttress anchor | |||
zone. | |||
Bearing area: | |||
2 2 2 b in. 481 4 11-24 A== Concrete anchor surface: | |||
2 2 2 b in. 994 4 11-33 A= Allowable bearing stress: | |||
()ksi 3.82 1/3 481 994 5.0 0.6 f cp== Use of 28 day concrete strength is conservative since initial tensioning will occur significantly | |||
later in time. | |||
TMI-1 UFSAR APPENDIX 5D 5D-4 REV. 18, APRIL 2006 Bearing pressure for anchor force: | |||
: a. ksi 3.40 481 1633.5 f kips, 1633.5 P b===<ksi 3.82 f cp= b. ksi 3.88 481 1866.8 f kips, 1866.8 P b=== The depth below the bearing plate, at which concrete stresses are identical to those allowed beyond the anchorage zone area, is computed on the conservative assumption that the load | |||
spreads at a 30 o angle from the plane normal to the bearing plate surface. | |||
Obtain the required concrete area to satisfy the above condition for: | |||
P = 1633.5 kips, A g = 1633.5 = 727 in 2 2.25 | |||
If it is conservatively assumed that the anchor force spreads at an angle of 30 o , the depth (h) below the bearing plate at which the concrete stress will have been reduced to the allowable | |||
compression stress as shown above, can now be computed. See Figure 5D-8b. +=4 7-)30 2htan (20.5 A 2 2 g in. 4.5 h 38.5, -3 2h 20.5 622 2=+= Special provisions must be made to confine the concrete within the 4.5 in. below the bearing plate. From then on, spiral reinforcement m eeting the requirement for tensile splitting stresses will be investigated as described in part 3.3. | |||
===3.2 EQUIVALENT=== | |||
FLUID PRESSURE METHOD | |||
Directly below the bearing plate the necessary reinforcement for confining the concrete is | |||
computed for an assumed equivalent fluid pressure "p". The equivalent fluid pressure is equal | |||
to the concrete compression stress multiplied by the Poissons Ratio, which is assumed to be | |||
0.2. | |||
The concrete compression stress 2 in. below the bearing plate,which is about at the center of | |||
the critical 4.5 in. depth, will be used to determine the equivalent fluid pressure. The initial | |||
1,400 kips anchor force is considered. See Figure 5D-8b. | |||
[]2 2 2 g in 481.5))4 ((7- )30 (tan 4 5.20 (A=+= ksi 2.91 481.5 1400 f c== The equivalent fluid pressure will be: | |||
p = 2.91 (0.2) = 0.58 ksi TMI-1 UFSAR APPENDIX 5D 5D-5 REV. 18, APRIL 2006 The reinforcement is arranged as shown on GAI Drawing 421031 . The #5 spiral has a | |||
diameter of 20 inches with a 3 1/2 in. pitch, and a total length of 36 inches. Horizontal ties | |||
2-#10 and vertical reinforcement 1-#10 plus 1-#8 are also available to resist this equivalent fluid pressure. | |||
The required steel area for the 20 inch diameter spiral was computed for an allowable steel | |||
stress of f s = 20 ksi. | |||
Required steel area: | |||
A s = 0.58(4.5)(20) | |||
= 2.61 in 2 (steel area req'd. over a 4.5 in length below the bearing plate) 20 | |||
Steel provided: | |||
: a. Radially, Spiral 2-#5 = 0.62 in 2 Ties 2-#10 = 2.54 3.16 in 2 > 2.61 in 2 | |||
: b. Vertically, Spiral 2-#5 = 0.62 Vert. 1-#10 = 1.27 Vert. 1-#8 0.79 2.68 in 2 > 2.61 in 2 | |||
===3.3 TRANSVERSE=== | |||
TENSION-PRISM METHOD | |||
Transverse stresses are produced when the anchor force spreads out into the surrounding | |||
concrete. Directly behind the bearing plate the stress lines are parallel to the tendon axis. The | |||
angle change of the stress-flow-lines, where they bend outwards, produces in these regions | |||
compression in the concrete which is of no concern. Further away, the stress-flow-lines tend to | |||
bend back and tensile transverse stresses are created as a result of this angle change. Figure | |||
5D-9 and 5D-10 show the idealized compression stress-flow-lines in plane and the approximate location of transverse stresses. | |||
Further, it can be seen that the angle changes will be small, and that the stress intensity along | |||
the compression stress-flow-lines reduces as they spread out into the concrete mass. From | |||
these considerations it can be concluded that the resulting tensile transverse stresses will be | |||
small and within the tensile capacity of the concrete. Deformed vertical, circumferential, and | |||
radial reinforcement is provided throughout the buttr ess in addition to the spiral reinforcement at each bearing plate. The reinforcement will tie the concrete mass together and ensures that | |||
thermal cracking of the concrete mass, if it should occur, is of no significance as far as structural | |||
integrity of the Reactor Building is concerned. | |||
Figure 5D-10 shows that the critical condition, with respect to transverse tensile stresses in the | |||
vertical plane, exists when only individual tendons ar e stressed. When all tendons are stressed, the stress-flow-lines are more or less parallel. Figure 5D-9 and 5D-10 show, superimposed | |||
over the anticipated stress-flow-lines, "equivalent beams" for which the transverse stresses can be computed. Computation for transverse stresses in radial planes are performed in tabulated | |||
form in Table 5D-1. Three cases are analyzed as follows: | |||
TMI-1 UFSAR APPENDIX 5D 5D-6 REV. 18, APRIL 2006 Case 1 Finds the maximum transverse radial stress in symmetrical equivalent beam. Variation la | |||
investigates anticipated conditions during normal operation. Variations 1b and 1c are made for | |||
extremely shallow and extremely deep "equivalent beams" to establish how variations of beam depths affect the magnitude and location of the transverse stresses. | |||
Case 2 Investigates the effect of eccentric loads on equivalent beams. Guyon's "Symmetrical Prism | |||
Method" is used. According to this method the depth of the eccentric beam does not affect the | |||
magnitude of radial transverse stresses; only the dimension of the prism drawn symmetrical to the tendon center line determines their magnitude. Possible variations in eccentric beam | |||
depths are therefore of no concern. | |||
All stresses were computed on the conservative assumption that the load spreads at 30 o in the circumferential direction. The stress computations confirm the previously made prediction that transverse stresses are within the acceptable concrete stress range. | |||
Case 3 Finds the maximum transverse vertical stresses in symmetrical equivalent beam. Variations 3a | |||
and 3b investigate anticipated conditions during normal operation. Vertical transverse forces | |||
counteract each other under operating conditions largely through vertical compression forces in | |||
the buttress concrete. See Figure 5D-10b. Variation 3c investigates anticipated conditions if | |||
only one anchor force acts. See Figure 5D-10a. | |||
====3.3.1 REINFORCEMENT==== | |||
FOR TRANSVERSE FORCES | |||
The allowable concrete tensile stresses in the radial and vertical planes are not exceeded and | |||
no special transverse reinforcement is required for design loads. However, reinforcement is | |||
available in the form of buttress reinforcement in the hoop and meridional direction consisting of a radial tie-reinforcement anchored into the shell and spiral reinforcement in the maximum | |||
bursting zone. | |||
The buttress is reinforced sufficiently to resist all transverse tensile forces Zy (design) acting | |||
within the limits of the buttress without relying on concrete tensile resistance. | |||
Spiral: 5/8" ?, 3 1/2" pitch, d = 20", length = 36" Maximum tensile splitting stress, Table 5D-1, f by =0.178 ksi | |||
T = 3.5 (0.178) (20) = 6.24 K , f s = 6.24 = 20.1 ksi 2 0.31 | |||
The spiral is designed for Zy (design) = 420 kips | |||
Maximum transverse tensile force as computed is Z y = 311 kips. | |||
TMI-1 UFSAR APPENDIX 5D 5D-7 REV. 18, APRIL 2006 Steel provided: | |||
Spiral 22 - #5 = 6.82 in 2 Ties 10 - #10 = 12.70 19.52 in 2 | |||
f s = 420 = 21.5 ksi 19.52 | |||
Temperature changes in the concrete during the life of the plant or at accident condition will | |||
produce a maximum of less than 1 ksi stress in the bursting zone reinforcement. | |||
===3.4 UNBALANCED=== | |||
FORCES | |||
Figure 5D-11 shows the tendon anchor forces for six buttresses in the Reactor Building. | |||
The solid line horizontal force vector in Figure 5D-11 represents the tangential components of | |||
the anchor forces. The anchor forces are applied on opposite faces of the buttresses and | |||
counteract each other. However, since the vert ical spacing of the anchors varies, it is not possible to balance all the forces. The unbalanced horizontal forces must be transferred from | |||
the buttress into the Reactor Building wall. | |||
The two dashed force vectors, Figure 5D-11, adjacent and parallel to each anchor force vector | |||
are obtained by proportioning the anchor force according to its position with respect to two | |||
opposing anchors similar to determining the react ions of a simply supported beam with an unsymmetrical concentrated load. The dashed horizontal force vector on the buttress center | |||
line indicates the unbalanced force between two opposing anchor forces. How these horizontal | |||
forces are resisted is discussed later. | |||
Since the horizontal anchor forces on opposing sides are offset from each other, the forces | |||
must change direction in order to balance each other. Transverse forces in vertical direction are | |||
a result of this angle change. The vertical dashed force vectors shown in Figure 5D-11 | |||
represent the unbalanced vertical forces behind each anchor zone. | |||
Vertical forces in Figure 5D-11 were computed for the simplified assumption that the anchor | |||
forces travel in straight lines from anchor plate to anchor plate after first moving 18 in. | |||
horizontally into the concrete and gradually bend over into the inclined position. See Figure 5D-12. | |||
The check analysis employing the Finite-Element Me thod did reveal a limited region of surface tensile stresses which presumable would produce cracking. However after concrete cracking | |||
has occurred with the extent limited by s upplemental reinforcing a new state of equilibrium will exist. The following section discussing vertical and horizontal unbalanced forces verifies that a | |||
mechanism does exist to maintain equilibrium. | |||
====3.4.1 VERTICAL==== | |||
UNBALANCED FORCES | |||
During initial prestressing and during tendon inspection the balance of opposing tendon anchors | |||
might not exist and the total vertical transverse stress must be resisted. The maximum concrete | |||
transverse TMI-1 UFSAR APPENDIX 5D 5D-8 REV. 18, APRIL 2006 tensile stress for this condition was computed in Table 5-2, Case 3c, f by = 115 psi. Figure 5D-10a shows the anticipated typical vertical stress-flow-lines for such conditions. The figure | |||
also shows the "equivalent beam" for which transverse stresses were determined according to | |||
Guyon's Symmetrical Prism Method. Sufficient reinforcing steel in form of a spiral has been | |||
provided for Case 3c without relying on concrete tensile resistance. In addition, the vertical | |||
reinforcement in the buttress, consisting of 7-#10, As = 8.8 in.,(Reference 2) can resist a | |||
transverse splitting force of Z = 8.8 (0.95) (40) = 335 kips which is slightly more than the | |||
expected force Zy = 311 kips shown in Table 5D-1. The reinforcement is placed mainly on the buttress surfaces in order to provide access for concrete placement. There is no need to | |||
distribute this reinforcement throughout the cross-se ction of the buttress or to place it where the vertical transverse forces actually act. Wherever the reinforcement is placed, it will not prevent the concrete from cracking and will only tie t he concrete together preventing the development of uncontrolled cracks. Consequently, surface reinforcem ent is sufficient. It restrains the buttress in the vertical direction and no failure mechanism caused by vertical transverse stresses can | |||
develop. | |||
The highest vertical in-plane shear stresses at the buttress-cylinder intersection occur directly | |||
beneath the ring-girder and at the cylinder-base s lab intersection haunch area. An amplification factor of 1.64 was applied to the shear stresses in both areas to account for the increased | |||
stiffness at the buttress. Sufficient reinf orcing was provided between the buttress and cylinder to take the shear stresses. | |||
====3.4.2 HORIZONTAL==== | |||
UNBALANCED FORCES | |||
Unbalanced horizontal forces shown in Figure 5D-11 must be transferred into the wall. | |||
During initial prestressing and tendon inspection the counteracting effect of opposing anchor | |||
forces might be lost at any any anchor zone position. All anchor zones will therefore be reinforced as the analysis requires for such conditions. Figure 5D-9a shows the anticipated | |||
typical horizontal stress-flow-lines for such conditions. The figure also shows the "equivalent beams" for which transverse stresses were determined according to Guyon's "Symmetrical Prism Method." | |||
The analysis for horizontal transverse forces is tabulated in Table 5D-1. The critical loading | |||
condition exists when a tendon is stressed without being counterbalanced by opposing anchor | |||
forces. The maximum concrete transverse tens ile stress for this condition is represented by Case 2a, f by = 178 psi. Sufficient reinforcing steel in form of a spiral has been provided for Case 2a without relying on concrete tensile resistance. In addition, the radial tie reinforcement | |||
in the buttress, consisting of 10-10#, A s = 12.7 in., (Reference 2) can resist a transverse splitting force of Z = 12.7 (0.95) (40) = 482 kips which is larger than the expected force Zy = 311 kips | |||
shown in Table 5D-1. The tie reinforcement restrains the buttress in radial direction and no | |||
failure mechanism caused by radial transverse stresses can develop. | |||
The highest shear stresses for the case of no opposing prestress force in the buttress occurs at | |||
the buttress-cylinder discontinuity on the side of the buttress opposite the existing prestress force. Reinforcing of #10 at 12 in. (See Figure 5D-4) are provided to take the resulting tensile | |||
stresses. | |||
TMI-1 UFSAR APPENDIX 5D 5D-9 REV. 18, APRIL 2006 | |||
====3.4.3 EQUILIBRIUM==== | |||
CONDITION IN THE ANCHOR BLOCK | |||
The reinforcement is checked based on the assumption of an existing spalling crack in the | |||
vertical direction. The forces and the extent of cracking in the anchor block were based on a | |||
method developed by Gergely, Soyen, and Siess of arresting cracks due to transverse stresses (Reference 4). The idealized mathematical model used for this approach, see Figure 5D-13, is | |||
quite conservative. The computation is shown in Table 5D-2. It should be noted that zero | |||
moment exists at B' when c = 27 in., and that the governing moment, Mm, prevails at c = 42 in., that is, far beyond the face of the wall. | |||
Concerning the visco-elastic behavior of concrete, it can be stated that the resulting | |||
redistribution of stresses developed by Gergely will not significantly affect reinforcing requirements (Reference 5). The super-position of bursting stresses and normal stresses from vertical tendons results in a tri-axial state of stress which is less critical than during construction | |||
stages. No additional investigations for the final state of stress in the anchorage zones, therefore, appear necessary. The various sour ces of errors in a numerical tri-axial analysis (uncertainties of material properties, element definition, etc.) may lead to less accurate values | |||
than were obtained by two-dimensional methods. | |||
The shear across the cylinder at the buttress-cylinder intersection for the case of unbalanced | |||
prestressed force has been checked. Sufficient reinforcing had been provided to account for | |||
the shear stresses. For accident temperature and pressure conditions a crack across the entire | |||
cylinder at the buttress cylinder discontinuity has been assumed. Sufficient steel reinforcing is provided to take any tension where the compressive stresses due to prestress have been | |||
overcome. | |||
4 CHECK ANALYSIS BY FINITE-ELEMENT METHOD | |||
==4.1 INTRODUCTION== | |||
The finite-element method has been employed to check the design of tendon anchorage zones | |||
in the Reactor Building. The purpose of this analysis was not as a design basis but simply as a | |||
check to demonstrate the conservatism of the method described here before. This section | |||
presents the method of analysis, load conditions, and results for the buttress and adjacent wall, and in conclusion to correlate these results with the design of the buttress reinforcement which | |||
has been presented in Section 3. | |||
The dominate loads are: | |||
: a. Prestress forces | |||
: b. Thermal loads | |||
: c. Pressure forces due to loss of coolant accident (LOCA) | |||
Although the stress analysis of the buttress and adjacent wall is a three dimensional problem, it can be reduced to a two dimensional problem, as the meridional stress of the cylindrical portion of the vessel is usually constant. The finite-element method of plane stress analysis (References 6 and 7) which can account for complex loadings, arbitrary geometric configurations, and arbitrary boundary conditions, assumes a homogeneous material (i.e., no TMI-1 UFSAR APPENDIX 5D 5D-10 REV. 18, APRIL 2006 cracks) and is an ideal tool for dealing with such a problem. Therefore, a plane stress finite-element computer program (Referenc e 7) was used in the buttress analysis. | |||
For the plane stress finite-element program t he thermal loads were obtained from another computer program. The results from this progr am, for the steady state winter normal operation and the transient state during winter start-up, w inter shut-down, and winter loss of coolant accident were used as input data in the finite- element program. | |||
As described later, a system of 768 direct st iffness simultaneous equations required solution. | |||
To save computer storage, the Gauss-Seidel iter ative method (Reference 8), which requires the storage of only the non-zero elements of the sti ffness matrix, was employed. In order to expedite convergence, an over-relaxation factor (R eference 9) was also used in Reference 7; this factor usually ranges from 1.0 to 2.0. | |||
Cracking of the concrete for the winter LOCA load condition is predicted. Areas of high tensile | |||
stresses have been defined and reinforcement capable of resisting these provided. | |||
===4.2 PLANE=== | |||
STRESS FINITE-ELEMENT COMPUTER PROGRAM | |||
Due to the fact that the meridional stresses of the cylindrical portion of a containment vessel are | |||
essentially uniform and constant, the horizontal cross section of the cylinder may be considered | |||
as a section subjected to plane stresses only. The finite-element method is a versatile method | |||
of analysis and is capable of providing states of stress in the buttress for complex loadings, arbitrary geometric configurations, and arbitrary boundary conditions. Tendon holes are not | |||
included in the buttress geometry. | |||
In the plane stress finite-element analysis, the horizontal cross sectional plane of the cylindrical portion of the vessel is replaced by a finite number of elements interconnected at a finite number of nodal points. Continuity between element s of the system is maintained by requiring that within each element, "lines initially straight remain straight in their displaced positions." | |||
This requirement is satisfied if the strains within each element are assumed to be constant. | |||
Therefore, the stresses within each element are also constant. Distributed loads are replaced | |||
by equivalent nodal forces which act at the three nodes of the element. Based on these | |||
assumptions, it is possible to derive the sti ffness matrix of a typical element, which is an expression for the nodal forces resulting from unit nodal displacements: | |||
{p} = [k] {?} (1) | |||
where {p} is the vector of six nodal forces which corresponds to the vector {?} of six nodal | |||
displacements (Figure 5D-14), and (k) is the stiffness matrix of the element. | |||
The nodal force vector {p} can be expressed as: | |||
{p} = {PX i , Px j , Px k , Py i , Py j , Pyk}) (2) and the nodal displacements vector {?} can be expressed as, | |||
{?} = {u i , u j , u k , v i , v j , v k} (3) which are shown in Figure 5D-14. | |||
TMI-1 UFSAR APPENDIX 5D 5D-11 REV. 18, APRIL 2006 Formulation of Element Stiffness Matrix | |||
Corresponding to the above assumptions, the displacements of a point at (x,y) are written as | |||
linear functions of the coordinates: | |||
u = u i + C 1 x + C 2 y v = v i + C 3 x + C 4 y (4) | |||
The stiffness matrix [k] can be obtained by the following steps: | |||
: a. Substitute the nodal coordinates into Equation(4) to obtain the displacement functions in terms of the nodal displacements. | |||
: b. Substitute results of (a) into the linear strain-displacement equations to obtain the strains in terms of nodal displacements. | |||
: c. Substitute the results of (b) into the stress-strain equation to obtain the stress in terms of nodal displacements. | |||
: d. Find the nodal resultant forces which are statically equivalent to the stresses found in (c) so that the relationships between the nodal forces and the nodal | |||
displacements can be written in a matrix form as shown in Equation (1). | |||
The explicit formulations of the matrices which compose the element stiffness matrix can be found in Reference 7. | |||
Solution of the Total Stiffness Equations After the stiffness matrices are formulated for each of the elements, they are assembled to form a total system stiffness matrix. If the boundary conditions and the loading conditions are properly specified, the tot al stiffness simultaneous equations are readily solvable to give the nodal displacements. Then, the nodal displacements are used to generate the stresses in each element and at each node. The directions and magnitudes of the principal | |||
stresses can also be computed by using the Mohr's circle technique. | |||
Iterative Method and Optimum Over-Relaxation Factor The Gauss-Seidel iterative procedure (Referenc e 8) is applied to solve the stiffness simultaneous equations. The method requires the storage of only the non-zero elements of the | |||
stiffness matrix during the iterative process. | |||
Since the stiffness matrix [k] is positive definite,the iterative method will always converge (Reference 10). The covergence rate is, however, often rather slow. An over-relaxation factor (Reference 9) can be applied which greatly incr eases the convergence rate. The determination of the proper over-relaxation factor is diffi cult and is based mostly on experience and trial and error. | |||
===4.3 GEOMETRY=== | |||
OF BUTTRESSES WITH TENDONS | |||
The pertinent geometry of the Reactor Building as it affects the buttress analysis is shown in Figure 5D-15. In Figure 5D-15a, a sectional view illustrates the vessel with six buttresses. | |||
Figure 5D-15b, shows the arrangement of the hoop tendons on a stretch out of the cylindrical TMI-1 UFSAR APPENDIX 5D 5D-12 REV. 18, APRIL 2006 wall. Figure 5D-15c, shows the dimensions of one buttress and the arrangement of the hoop and vertical tendons. | |||
===4.4 THERMAL=== | |||
LOADS AND TEMPERATURE PROFILES | |||
If the temperature profile in a buttress is known, the difference between the temperature of any point and the reference temperature can be determined. The temperature difference, with the | |||
modulus of elasticity and thermal expansion coeffi cient of the concrete, provides the equivalent constant stress for each finite-element. The constant thermal stress of each element can be | |||
converted to equivalent nodal resultant forces which can be considered as external forces | |||
applied at the nodal points of a non-heated element. | |||
Temperature Profiles The temperature profiles and iso-thermal curves were calculated by a computer program developed by Brunnel (Reference 11) for the s olution of steady state and transient state problems in one, two, or three dimensional Cart esian or cylindrical systems. The governing differential equations were solved by the method of finite differences, using standard relaxation techniques. | |||
The convergence criterion: | |||
T j (step i) - T j (step i-1) < | |||
0.00001 (5) | |||
Y j (step i) | |||
where T j is the temperature at node j, was applied at every node of the difference lattice. To ensure stability of the transient solutions, the time increment was selected according to the | |||
method presented in Reference 11. | |||
The buttress and wall temperatures were determined for the following conditions: | |||
: a. Normal operation in winter | |||
: b. Start-up in winter | |||
: c. Shut-down in winter | |||
: d. Loss of coolant accident during winter | |||
The winter season was chosen because it produces the greatest thermal gradients in the shell | |||
of the Reactor Building. | |||
The mean monthly temperature was chosen due to t he fact that day- night fluctuations had little effect on the interior portion of the concrete. According to Carslow and Jeager, (Reference 12) | |||
the amplitude of the temperature diminishes as: | |||
wP con C p con exp [X | |||
] (6) 2 k con Thus a distance of 12.8 inches from the outer wall of the Reactor Building the effect of the | |||
oscillation is only 0.1 what it is on the outer face, and the middle and inner portion of the vessel | |||
will be unaffected by day-night temperature oscillations. | |||
TMI-1 UFSAR APPENDIX 5D 5D-13 REV. 18, APRIL 2006 For the start-up condition, it was assumed that the inside air temperature of the Reactor Building rose from 60 F to 110 F in one hour, and for shut-down the reverse was used. The | |||
inside film coefficient was taken as 1.490 Btu/hr/sq ft/F for cases a, b, and c above. This was | |||
modified film coefficient, which took into account the steel liner and the still air between the liner | |||
and the concrete: | |||
h m = 1 (7) 1 X steel X air h i + | |||
+ k steel k air Equation 7 neglects the circumferential conduction and heat capacity in the steel liner and air | |||
gap. However, due to the relatively small thi ckness of the liner and the air gap when compared with the thickness of the concrete wall, these effects may be neglected. The preliminary | |||
analyses verified this assumption. The following values, obtained by theoretical analysis (Reference 13), were used in the thermal analysis: | |||
inside film coefficient = 2.0 Btu/hr/sq ft/F | |||
specific heat, C p con = 0.210 Btu/lb/F | |||
thermal conductivity, k con = 10.0 Btu in./hr/sq ft/F | |||
k steel = 300.0 Btu in./hr/sq ft/F | |||
k air = 1.836 Btu in./hr/sq ft/F | |||
The density of the concrete, P con was assumed to be 140 lbs/cu ft. | |||
The above values agree closely with the results from experimental research carried out in Japan (Reference 14). | |||
C p con = 0.200 Btu/lb/F | |||
k con = 11.28 Btu in./hr/sq ft/F | |||
The density of the concrete P con , was assumed to be 146 lb/cu ft. | |||
The thickness of the steel and air gap was: | |||
X air = 1/32 inch | |||
X steel = 3/8 inch | |||
Using Reference 15, the local film coefficient for air at 20 F and a free stream velocity of 20 mph | |||
flowing past a concrete cylinder was calculated to be: | |||
h = 0.56 Btu/hr/ft 2/F As a safe approximation the outside film coefficient was chosen as: | |||
TMI-1 UFSAR APPENDIX 5D 5D-14 REV. 18, APRIL 2006 h? = 1.0 Btu/hr/sq ft/F | |||
The buttress model was replaced by a mesh of 61 nodes, with adiabatic boundaries at the | |||
centerline of the buttress and at the radial mid-section of the wall between two neighboring | |||
buttresses (Figures 5D-16 and 5D-17). Due to symmetry, there is no heat transferred across | |||
these boundaries. For Case d the liner temperature during LOCA was used as the inner | |||
boundary temperature in conjunction with a film coefficient of 5.880 Btu/hr/sq ft/F, which | |||
simulated the conduction through the air gap: | |||
k air h m = (8) X air The transient state temperature analyses were performed at various time steps. | |||
The temperature profiles for Cases a, b, c, and d are shown in Figure 5D-18. Because the most | |||
severe temperature gradients occur for case a and d, only these two cases were considered in | |||
the stress analyses. Iso-thermal curves fo r cases a and d are shown in Figures 5D-19 and 5D-20. | |||
The design pressure inside the Reactor Building during LOCA is equal to 55 psig and the | |||
maximum pressure was initially calculated to occur at approximately 200 seconds after initiation | |||
of the accident. The maximum liner temperature of 280F was initially calculated to occur at | |||
approximately 650 seconds after the initiation of the accident. However, the maximum | |||
temperature gradient, which causes the most severe thermal stresses, was calculated to occur | |||
at 10,000 seconds after initiation of the accident. | |||
Consequently, the determination of the most severe stress state, due to thermal and pressure loads, requires that a number of stress analyses be conducted at various elapsed times after | |||
commencement of the accident. A conservative estimate of the stresses was obtained, by | |||
superimposing directly the maximum values for accident temperature and pressure. | |||
Consequently although subsequent calculations alter t he time at which minimum pressure and liner temperature occur, these changes do not affect this analysis. | |||
4.5 NON THERMAL LOADS | |||
In addition to the thermal loads, there are other loads which act on the structure: | |||
: a. Hoop Tendon Anchor Pressure | |||
The tendon force was assumed to be distributed uniformly up the face of the buttress. The resulting average force (=p a) was applied to each of six nodes. | |||
The location of the nodes approximately corresponded to the position of the | |||
tendon. | |||
: b. Hoop Tendon Pressure in a Radial Direction | |||
Radial pressures, due to the hoop tendons in the wall, were calculated as average stresses in pounds/in. | |||
The two radial pressures applied (refer to Figure 5D-21): | |||
TMI-1 UFSAR APPENDIX 5D 5D-15 REV. 18, APRIL 2006 | |||
: 1. From point A to the tangent point T. Over this length the tendon spacing causing radial pressure (i.e., continuous past the buttress) is 33 in., giving | |||
a pressure of p t , and | |||
: 2. For the length from point T to point B (midway between buttresses) the average tendon spacing causing radial pressure is 16 1/2 in., giving a | |||
pressure of 2p | |||
: t. | |||
: c. Internal Pressure Due to LOCA | |||
Internal pressure acting on the inside face of the Reactor Building Wall was considered as a combination of two separate effects: | |||
: 1. Pressure due to LOCA. The maximum pressure due to LOCA was taken at 82.5 psig. (1.5 P) | |||
: 2. Pressure due to expansion of the liner caused by the increase in temperature during the LOCA condition (110F to 280F) based upon a limiting yield strength of 1.2 f | |||
: y. | |||
The above two loads were assumed to act simultaneously during the LOCA condition, giving an equivalent internal pressure of p | |||
: i. | |||
The magnitude and combination of the above non-thermal loads is noted in section 4.7 | |||
===4.6 MATERIAL=== | |||
PROPERTIES | |||
Prestressing Steel Wire | |||
The type of prestressing wire to be used in the construction is a stress relieved wire. | |||
Dimensions, properties, etc., of the wire as used in the analysis were: | |||
: a. Size of tendon = 169 wires | |||
: b. Diameter of wires = 1/4 in. | |||
: c. Minimum guaranteed ultimate tensile stress = 240,000 psi | |||
: d. Tendon force @ 70 percent GUTS = 1400 kips | |||
: e. Modulus of elasticity = 29 x 10 6 psi f. Coefficient of friction = 0.16 | |||
: g. Coefficient of wobble friction = 0.00030 TMI-1 UFSAR APPENDIX 5D 5D-16 REV. 18, APRIL 2006 Concrete The assumed properties of the concrete used in the analysis were: | |||
: a. Compressive strength @ 28 days = 5000 psi | |||
: b. Modulus of elasticity = 4 x 10 6 psi | |||
: c. Poissons ratio = 0.20 | |||
Liner The liner is fabricated from 3/8 inch thick steel plate manufactured to ASTM A-283 Grade C. | |||
The minimum yield stress for this material is 30,000 psi. In determining the equivalent liner pressure LOCA, the yield stress was assumed to be 36,000 psi. | |||
===4.7 LOADING=== | |||
CASES AND COMBINATIONS | |||
The buttress and wall were analyzed for four loading cases. The first two cases considered the | |||
buttress and wall being subjected to the maximum prestress forces, i.e., prestress forces | |||
allowing only those losses that would occur between prestressing and start-up (beginning of | |||
plant life.) The latter two cases considered the buttress and wall being subjected to the | |||
minimum prestress forces, i.e., prestress forc es after 40 year losses (end of plant life). | |||
The load cases were: | |||
: a. Winter normal operation with maximum prestress force. | |||
: b. Winter LOCA with maximum prestress force. | |||
: c. Winter normal operation with minimum prestress force. | |||
: d. Winter LOCA with minimum prestress force. | |||
These cases were considered to be more critical than the summer normal operation/LOCA | |||
conditions. The summer condition results in a smaller temperature differential in the wall which | |||
in turn results in maximum compressive and tensile stresses less than the winter condition. | |||
Hence, the summer condition was not considered in the analyses. | |||
TMI-1 UFSAR APPENDIX 5D 5D-17 REV. 18, APRIL 2006 The numerical value of the applied loads was: | |||
LOAD CASE a b c d | |||
t LT 110 150 110 150 | |||
t LP 110 280 110 280 | |||
t e 20 20 20 20 | |||
p a 6.25 6.73 5.52 6.03 | |||
pt A-T 46.07 49.75 39.25 42.2 | |||
2ptT-B 92.13 99.50 78.5 84.4 | |||
p i 0.0 96.1 0.0 96.1 | |||
where: | |||
t LT = Temperature (F) of liner for thermal gradient in wall. | |||
t LP = Temperature (F) of liner for equivalent pressure of liner due to expansion during LOCA. | |||
t e = External air temperature (F) | |||
p a = Prestress force at bearing plate (kips/inch) | |||
pt A-T = Radial prestress force from point A to tangent point T (lbs/inch) | |||
2p T-B = Radial prestress force tangent point T to point B (lbs/inch) | |||
p i = Equivalent internal pressure (LOCA pressure + liner expansion pressure lbs/inch) | |||
Figure 5D-21 shows the application of these loads on the buttress and wall. | |||
Load p a p t , and 2p t for Case a incorporate losses due to elastic shortening of the concrete, friction and creep, shrinkage, and steel relaxation losses that occur between prestressing and | |||
start-up. Case b prestress forces reflect increases from Case a, due to elastic elongation of the | |||
structure during LOCA. Thus, Case a and b reflect the maximum prestress force that will exist | |||
in the structure, at the earliest time in the life of the plant, that a LOCA could occur. | |||
Loads p a , p t , and 2 pt for Cases c and d are similar to Cases a and be respectively, except that the losses due to creep, shrinkage, and steel relaxation are those that are assumed to occur | |||
over 40 years to end of plant life. Case d pres tress forces reflect an increase from Case c due to elastic elongation of the structure during LOCA. | |||
TMI-1 UFSAR APPENDIX 5D 5D-18 REV. 18, APRIL 2006 | |||
===4.8 STRESS=== | |||
ANALYSIS | |||
Due to the symmetrical nature of the structur e only one-twelfth of the cylinder, as shown in Figure 5D-16, need be analyzed. Preliminary analy ses showed that beyond a certain point in the wall, the stress state, for all loading cases, was uniform. To conserve computer storage and | |||
give a more accurate solution, it was decided to reduce the size of the model to that shown in | |||
Figure 5D-22. Nodal displacements of both boundaries were specified. The displacements | |||
were obtained form a computer program as developed by Kalnins (Reference 16) using a shell of revolution theory. | |||
There were a total of 672 elements and 384 nodal points, so that 768 simultaneous equations | |||
required solution. For the computations, the Gauss-Seidel iterative solution procedure was used. | |||
===4.9 ANALYSIS=== | |||
OF RESULTS | |||
Figures 5D-23 through 5D-30 show iso-stress curves for the maximum and minimum principal stresses for the four loading conditions outlined in Section 4.7. It should be noted that these | |||
curves represent the algebraic maximum and minim um values of the principal stresses. The do not show the orientation of the stress trajectories. | |||
Inspection of the iso-stress curves shows tha t, as expected, high compressive stresses exist directly beneath the bearing plate for all load conditions. Under winter normal conditions, a | |||
region of tensile bursting stress occurs in the buttress at some distance removed from the | |||
bearing plate; during LOCA, however, these bursting stresses merge with those resulting from the pressure and temperature loads. | |||
The variation of temperature through both the wall and buttress leads to substantial hoop stress | |||
gradients. It should be noted that in the curves for LOCA, the direction of the stresses in the | |||
wall changes in the vicinity of the centerline. | |||
Of primary interest, however, is the high localized t ensile stress that occurs at the buttress-wall re-entrant corner during LOCA, which results fr om both the thermal gradient as well as the local bending arising from the variation of resistance of the buttress and wall to internal pressure. | |||
4.10 CONCLUSIONS | |||
The results indicate that for the most severe loading conditions, the maximum stresses in the | |||
buttress and adjacent wall are within safe limits, except for the localized area at the buttress | |||
wall corner where tensile stresses greater than 6 ?f | |||
? occur. Reinforcement to stabilize cracking in this local region are provided. All other results would indicate that the design method | |||
employed was conservative. | |||
TMI-1 UFSAR APPENDIX 5D 5D-19 REV. 18, APRIL 2006 5 REFERENCES | |||
: 1. Building Code Requirements for Reinforced Concrete , ACI 318-63, American Concrete Institute, Detroit, Michigan, 1963. | |||
: 2. Y. Guyon, Prestressed Concrete , Contractors Record Ltd. London and John Wiley and Sons Inc., New York, 1960. | |||
: 3. F. Leonhardt, Prestressed Concrete Design and Construction , Wilhelm Ernest and Sohn, Berlin-Munich, 1964. | |||
: 4. P. Gergely, M. A. Sozen and C. P. Siess, The Effect of Reinforcement on Anchorage Zone Cracks in Prestressed Concrete Members , Civil Engineering Studies, Structural Research Series No. 271, University of Illinois, Urbana, Illinois, July 1963. | |||
: 5. Anchorage Systems in Prestressed Concrete Pressure Vessels; Anchorage Zone Problems , Oak Ridge National Laboratory, ORNL-TM-2378, January, 1969. | |||
: 6. M. J. Turner, R. W. Clough, H. C. Martin, and L. J. Topp, "Stiffness and Deflection Analyses of Complex Structures," Jl. of Aero. Sci., V. 23, No. 9, Sept 1965, p. 805. | |||
: 7. E. L. Wilson, Finite-Element Analysis of Two Dimensional Structures , Ph. D. Thesis, Univ. of Cal. at Berkeley, June 1963. | |||
: 8. M. G. Salvadori and M. L. Baron, Numerical Methods in Engineering , Prentice-Hall, Englewood Cliffs, N. J., 1961. | |||
: 9. F. G. Lehan, "Simultaneous Equations S olved by Over-Relaxization, "Proc. ASCE, 2nd Conf. on Electronic Computation , Pittsburgh, Pa., USA, Sept. 1960. | |||
: 10. S. H. Crandall, Engineering Analysis , McGraw-Hill Book Co., N. Y., USA, 1956. | |||
: 11. R. D. Brunnel, Computer Program Heating and Heating 2 , Atomic International Division, Report No. A1-64-MEMO-177, North American Aviation, Inc., Detroit, Michigan, USA, 1964. | |||
: 12. H. S. Craslaw and J. C. Jaeger, Conduction of Heat in Solids , 2nd Ed., Oxford, London, 1959. | |||
: 13. Johns-Manville Research and Engineering Center Report E-455-T-266 filed in Robert Emmett Ginna Final Facility Description and Safety Analysis Report, Section 5, Appendix 5B. | |||
: 14. J. Kosaka, Experimental Research of Unsteady-State Heat Transfer Concerned with Nuclear Energy Facilities (Preliminary Report), Division of Nuclear Energy Facilities, Taisei Construction Ltd., 1969. | |||
: 15. F. Kreith, Principles of Heat Transfer , 2nd Ed. International, Scranton, Pa., 1965. | |||
: 16. A. Kalnins, Computer Programs , Crystal River Unit No., 3, Florida Power Corporation, PSAR, Pages 5C-1 and 5C-2. | |||
: 17. J. Zielinski and R. E. Rowe, An Investigation of the Stress Distribution in the Anchorage Zones of Post-Tensioned Concrete Members , Cement and Concrete Association, Research Report NO. 9, London, September 1960. | |||
TMI-1 UFSAR APPENDIX 5D 5D-20 REV. 18, APRIL 2006 | |||
: 18. J. Zielinski and R. E. Rowe, The Stress Distribution Associated with Groups of Anchorages in Post-Tensioned Concrete Members , Cement and Concrete Association, Research Report No. 13, London, October 1962. | |||
: 19. P. Gergely, Anchorage Systems in Prestressed Concrete Pressure Vessels; Anchorage Zone Problems , U.S. Atomic Energy Commission, ORNL-TM-2378, January, 1969. | |||
TMI-1 UFSAR APPENDIX 5D 5D-21 REV. 18, APRIL 2006 TABLE 5D-1 BUTTRESS ANCHOR ZONE TRANSVERSE STRESSES IN RADIAL AND VERTICAL PLANES Case Radial Plane | |||
1a 1b 1c 2a | |||
2a'(in) 20.5 20.5 20.5 20.5 | |||
2a (in) 52 26 78 32 | |||
2a'/2a 0.40 0.79 0.26 0.64 | |||
fy[1] 0.27 0.10 0.34 0.15 | |||
z (in) [2] 0.80 0.90 0.70 0.9 | |||
2b' (in) 20.5 20.5 20.5 20.5 | |||
2b=2b'+(2az)tg30 o 44.5 33.5 51.5 37.0 A=2a2b (in | |||
: 2) 2315 870 4150 1182 P (kips) [3] 1400 1400 1400 1400 | |||
f bz = P/A (ksi) 0.605 1.610 0.338 1.182 f by = f bz fy (ksi) 0.164 0.161 0.115 0.178 [4] | |||
Zy = 0.3P (1-a'/a) 252 88 311 151 (kips) [5] | |||
Zy (design) = 0.3P 420 420 420 420 (kips) [6] | |||
az (in) 20.8 11.7 27.3 14.4 | |||
[1] Guyon, "Prestressed Concrete", Figure 9. | |||
2 [2] Guyon, "Prestressed Concrete", Figure 9. | |||
2 [3] 1600 Kips and 2000 kips anchor forces is not critical. | |||
[4] Changes in tensile splitting stresses with varying "2a" dimension. | |||
[5] Zy = total transverse force Leonhardt "Prestressed Concrete" | |||
: p. 271. | |||
3 [6] For design purpose "(a'/a)" is conservatively assumed to be equal | |||
to zero. | |||
TMI-1 UFSAR APPENDIX 5D 5D-22 REV. 18, APRIL 2006 TABLE 5D-1 (continued) | |||
BUTTRESS ANCHOR ZONE TRANSVERSE STRESSES IN RADIAL AND VERTICAL PLANES | |||
Case Ultimate Load- Vertical Planes Radial Plane 2a 3a 3b 3c | |||
2a' (in) 20.5 20.5 20.5 20.5 | |||
2a (in) 32 33 26 78 | |||
2a'/2a 0.64 | |||
fy [1] 0.15 | |||
z (in) [2] 0.9 | |||
2b' (in) 20.5 | |||
2b=2b' + (2az)tg30 0 37.0 See Case 2a, 1b, and 1c A=2a2b (in | |||
: 2) 1182 P (kips) [3] 2000 | |||
f bz = P/A (ksi) 1.690 f by = f bz fy (ksi) 0.254 [4] | |||
Zy=0.3P (1-a'/a) 216 (kips) [5] | |||
Zy (design)=0.3P 600 (kips) [6] | |||
az (in) 14.4 | |||
[1] Guyon, "Prestressed Concrete", Figure 9. | |||
2 [2] Guyon, "Prestressed Concrete", Figure 9. | |||
2 [3] 1600 Kips and 2000 kips anchor forces is not critical. | |||
[4] Changes in tensile splitting stresses with varying "2a" dimension. | |||
[5] Zy = total transverse force Leonhardt "Prestressed Concrete" p. 271. | |||
3 [6] For design purpose "(a'/a)" is conservatively assumed to be equal to zero. | |||
TMI-1 UFSAR APPENDIX 5D 5D-23 REV. 18, APRIL 2006 TABLE 5D-2 BUTTRESS ANCHOR ZONE SPALLING (Refer to Figure 5D-13) | |||
M = P {c - e - (c/h) 2 [2h - 3e - c + 2ec | |||
]} h c= varies e = 15.5 in. h = 69 in. | |||
M = P {c - 15.5 - (c/69) 2 [2 x 69 - 3 x 15.5 - c ( 1 - 2x15.5 | |||
)]} 69 | |||
C (in.) c - 15.5 (c/69) 2 -c( 1 - 2 x 15.5) M/P M m (in. - K) 69 | |||
5 -10.5 0.00525 - 2.75 -11.0 10 - 5.5 0.02100 - 5.50 - 7.3 15.5 0 0.05040 - 8.50 - 4.4 20 + 4.5 0.08400 -11.00 - 2.3 25 + 9.5 0.1310 -13.80 -0.7 zero at 27 in. 30 +14.5 0.1890 -16.50 + 0.3 42 +26.5 0.370 -23.10 + 1.2 1680 50 +34.5 0.525 -27.50 + 0.9 (in. - kips) 69 +53.5 1.00 -38.00 0 | |||
F t = M m Z = 3 in h = 69 in. | |||
h-z | |||
f s allowable = 10 5 W A s A s = 1.27 sq in. (#10 bars) | |||
W = 0.005 | |||
f s = 10 5 0.005/1.27 = 6.28 ksi F T = 1680 | |||
= 25.5 K 69 - 3 | |||
A s = F T/f s = 25.5/6.28 = 4.05 sq in. | |||
A s provided: #10 @ 12 in. + 1 #10 Loop = 6.03 sq in. | |||
TMI-1 UFSAR APPENDIX 5D 5D-i REV. 18, APRIL 2006 | |||
APPENDIX 5D REPORT ON STRESS ANALYSIS OF BUTTRESS | |||
TMI-1 UFSAR APPENDIX 5D 5D-ii REV. 18, APRIL 2006 TABLE OF CONTENTS | |||
SECTION TITLE 1 DESIGN DETAILS | |||
===1.1 TENDON=== | |||
ANCHORAGE HARDWARE | |||
===1.2 ANCHORAGE=== | |||
ZONE REINFORCEMENT | |||
2 TENDON LOADS AND ALLOWABLE STRESSES | |||
===2.1 TENDON=== | |||
LOADS | |||
===2.2 ALLOWABLE=== | |||
STRESSES | |||
3 ANCHOR BLOCK STRESSES | |||
===3.1 BEARING=== | |||
STRESSES | |||
===3.2 EQUIVALENT=== | |||
FLUID PRESSURE METHOD | |||
===3.3 TRANSVERSE=== | |||
TENSION-PRISM METHOD | |||
====3.3.1 REINFORCEMENT==== | |||
FOR TRANSVERSE FORCES | |||
===3.4 UNBALANCED=== | |||
FORCES | |||
====3.4.1 VERTICAL==== | |||
UNBALANCED FORCES | |||
====3.4.2 HORIZONTAL==== | |||
UNBALANCED FORCES | |||
====3.4.3 EQUILIBRIUM==== | |||
CONDITION IN THE ANCHOR BLOCK | |||
4 CHECK ANALYSIS BY FINITE-ELEMENT METHOD | |||
==4.1 INTRODUCTION== | |||
===4.2 PLANE=== | |||
STRESS FINITE-ELEMENT COMPUTER PROGRAM | |||
===4.3 GEOMETRY=== | |||
OF BUTTRESSES WITH TENDONS | |||
===4.4 THERMAL=== | |||
LOADS AND TEMPERATURE PROFILES 4.5 NON THERMAL LOADS | |||
===4.6 MATERIAL=== | |||
PROPERTIES | |||
===4.7 LOADING=== | |||
CASES AND COMBINATIONS | |||
===4.8 STRESS=== | |||
ANALYSIS | |||
===4.9 ANALYSIS=== | |||
OF RESULTS | |||
4.10 CONCLUSIONS | |||
5 REFERENCES | |||
TMI-1 UFSAR APPENDIX 5D 5D-iii REV. 18, APRIL 2006 LIST OF TABLES | |||
TABLE TITLE 5D-1 Buttress Anchor Zone Transverse Stresses in Radial and Vertical Planes | |||
5D-2 Buttress Anchor Zone Spalling | |||
TMI-1 UFSAR APPENDIX 5D 5D-iv REV. 18, APRIL 2006 LIST OF FIGURES | |||
FIGURE TITLE 5D-1 DELETED | |||
5D-2 DELETED | |||
5D-3 DELETED | |||
5D-4 DELETED | |||
5D-5 DELETED | |||
5D-6 DELETED | |||
5D-7 Typical Hoop Tendon Arrangement | |||
5D-8 Buttress-Hoop Tendon Bearing Stresses | |||
5D-9 Buttress-Radial Bursting Stresses | |||
5D-10 Buttress-Vertical Bursting Stresses | |||
5D-11 Hoop Tendon Typical Anchor Force Distribution | |||
5D-12 Hoop Prestress Idealized Vertical Stress Flow in Buttress | |||
5D-13 Idealized Model for Spalling | |||
5D-14 A Triangular Plane Stress Finite-Element | |||
5D-15 Reactor Building With Buttresses | |||
5D-16 Adiabatic Sections | |||
5D-17 Model of Gridwork for Thermal Analysis | |||
5D-18 Temperature Profiles Through Wall and Buttress for Various Thermal Load Conditions | |||
5D-19 Iso-Thermal Curves for Winter Normal Operating Temperature Condition | |||
5D-20 Iso-Thermal Curves for Winter A ccident Temperature Condition at Time = 10,000 Secs. | |||
5D-21 The Hoop Tendon Loads, Internal Pressure and the Equivalent Liner Pressure | |||
5D-22 Finite-Element Grid Work for Buttress Analysis TMI-1 UFSAR LIST OF FIGURES | |||
FIGURE TITLE APPENDIX 5D 5D-v REV. 18, APRIL 2006 5D-23 Minimum Principal Stresses for Maximum Prestress Forces and Winter Normal Operating Condition (Case A) (PSI) | |||
5D-24 Maximum Principal Stresses for Maximum Prestress Forces and Winter Normal Operating Condition (Case A) (PSI) | |||
5D-25 Minimum Principal Stresses for Maximum Prestress Forces and Winter Accident Condition (Case B) (PSI) | |||
5D-26 Maximum Principal Stresses for Maximum Prestress Forces and Winter Accident Condition (Case B) (PSI) | |||
5D-27 Minimum Principal Stresses for Minim um Prestress Forces and Winter Normal Operating Condition (Case C) (PSI) | |||
5D-28 Maximum Principal Stress for Minimum Prestress Condition (Case C) (PSI)Forces and Winter Normal Operating | |||
5D-29 Minimum Principal Stresses for Minim um Prestress Condition (Case D) (PSI) | |||
Forces and Winter Accident | |||
5D-30 Maximum Principal Stresses for Minimum Prestress Forces and Winter Accident Condition (Case D) (PSI) | |||
TMI-1 UFSAR APPENDIX 5D 5D-1 REV. 18, APRIL 2006 1 DESIGN DETAILS | |||
===1.1 TENDON=== | |||
ANCHORAGE HARDWARE | |||
The prestressing system employs the BBRV anchor age system utilizing 169-1/4" diameter wires conforming to ASTM A421-65, Type BA. The dimensions and capacities of the tendons used | |||
hereafter for the anchorage zone analysis are more fully described in Appendix 5B, Part D. | |||
===1.2 ANCHORAGE=== | |||
ZONE REINFORCEMENT | |||
The mild steel reinforcement in the tendon anchorage zones of the buttress and the base mat | |||
are shown on GAI Drawings 421004, 421005, 421006, 421031,421032, and 421033. This | |||
reinforcement was conservatively designed by the methods enumerated hereafter in Section 3 of this Appendix and is judged to be sufficient to accommodate all predicted stress states | |||
without jeopardizing the ability to satisfactorily place the concrete. | |||
2 TENDON LOADS AND ALLOWABLE STRESSES | |||
===2.1 TENDON=== | |||
LOADS | |||
The minimum guaranteed ultimate tensile strength of the tendons is approximately 2,000 kips. | |||
Before the tendons are anchored (locked off) at 70 percent of the ultimate tendon force, which is | |||
equal to 1,400 kips, they may be temporarily ov erstressed to 1,600 kips in order to reduce friction losses. The anchor zones are designed for anchor forces of 1,400 kips per tendon. The | |||
anchor zones are also checked to resist safely temporary anchor force of 1,600 kips per tendon. | |||
===2.2 ALLOWABLE=== | |||
STRESSES | |||
The allowable stresses for concrete and steel will conform to those given in Section 5.2.2 and | |||
5.2.3 of the FSAR. | |||
Due to the specialized nature of anchor zone design the following conditions and stress criteria | |||
will be used in the design. These stress criteria include adequate conservatism in their | |||
application. | |||
First, verify that the allowable concrete compression and tensile stresses are not exceeded due | |||
to bearing and bursting stresses respectively. | |||
Second, check the tensile stresses in reinforcing steel when no credit is taken for concrete | |||
tensile capacity in the investigated anchor zones. | |||
Under these conditions, the allowable stresses will be shown as follows: | |||
TMI-1 UFSAR APPENDIX 5D 5D-2 REV. 18, APRIL 2006 CONCRETE | |||
: 1. Allowable bearing stress | |||
according to ACI/318-63, 1 Paragraph 2605(C) | |||
: 2. Allowable compression stress in the structure beyond the anchorage zone c c f 0.45 f= | |||
: 3. Allowable Tension Stress | |||
REINFORCING STEEL | |||
Allowable tension for | |||
A615-Grade 40, f s = 20 ksi. | |||
3 ANCHOR BLOCK STRESSES The investigation of the buttress anchorage zones for the hoop tendons is described in the following sections. The analytical work as presented, is based on the end anchor zone stresses | |||
and stress flow pattern studies documented by Y. Guyon (Reference 2) and F. Leonhardt (Reference 3). Their work has been substant iated through extensive strain measurements and stress flow pattern established by studies utilizing photoelastic techniques. | |||
The design presented herein is quite conservative; nevertheless, it has been adopted as the | |||
basis for determination of required reinforcing steel in the end anchor zones to overcome tensile | |||
splitting stresses and to provide for spalling stresses near the concrete surface. | |||
The existance of a three-dimensional stress field including high shear, thermal and shrinkage | |||
tensile stresses is influenced by hairline cracking due to creep and thermal stresses. In order to | |||
accommodate the most unfavorable stress combination in the anchorage zones, the following | |||
design goals were established: | |||
: 1. Maximum concrete bursting tensile stresses have been kept below maximum allowable tensile stress of conservatively c f 3using Y. Guyon's symmetrical prism method. This would indicate that no special spiral and bursting reinforcement is required. The 3 b b ci cp/A A' f 6.0 f=ksi. 5.0 f ci=ksi. 5.0 f c=c f 3ksi. 5.0 f c= | |||
TMI-1 UFSAR APPENDIX 5D 5D-3 REV. 18, APRIL 2006 combined stresses during normal operation and test loading also indicate that the concrete will remain uncracked in the anchorage zone. | |||
: 2. However, due to the existance of complex thermal stresses and possible microcracking, the ability of concrete to carry tensile stresses has been discarded, and mild steel | |||
reinforcing has been provided for: | |||
: a. the equivalent fluid pressure | |||
: b. transverse tensile stresses | |||
: c. unbalanced vertical forces | |||
: d. unbalanced horizontal forces | |||
: e. equilibrium condition | |||
The combined stresses during the design pressure for the Reactor Building (p = 55 psig) indicate that the concrete will remain uncracked in the anchorage zone. The combined | |||
stresses during the factored accident loading (refer to 5.2.3.2, cases a. and b. indicate | |||
that the concrete in the anchorage zone may develop hairline cracks. | |||
: 3. Existing theoretical solutions for anchorage zone design do not agree with more recent experimental results (References 17 and 18); maximum transverse tensile forces were | |||
made up to 5.8 times larger than determined by Y. Guyon. | |||
: 4. Possible failure mechanisms due to spalling and bursting stresses were evaluated to verify static equilibrium and acceptable stre ss limits under such hypothetical conditions. | |||
: 5. Bursting stresses in a large body are much smaller than in a member bounded by planes (Reference 19). Therefore, while it was found that two-dimensional theories underestimate | |||
the bursting stresses (Reference 17), analytical results based on isolated prisms or better, experimental values such as those developed by J. Zielinski and R. E. Rowe can be used for massive structures (Reference 19). | |||
: 6. The design of the anchorage zones provides a margin of safety similar to the entire reactor building. Furthermore, the reinforcement is designed to control cracking due to | |||
spalling, bursting, or shear of the concrete under sustained load. | |||
===3.1 BEARING=== | |||
STRESSES | |||
Figures 5D-7 through 5D-10 show the typical hoop tendon arrangement and buttress anchor | |||
zone. | |||
Bearing area: | |||
2 2 2 b in. 481 4 11-24 A== Concrete anchor surface: | |||
2 2 2 b in. 994 4 11-33 A= Allowable bearing stress: | |||
()ksi 3.82 1/3 481 994 5.0 0.6 f cp== Use of 28 day concrete strength is conservative since initial tensioning will occur significantly | |||
later in time. | |||
TMI-1 UFSAR APPENDIX 5D 5D-4 REV. 18, APRIL 2006 Bearing pressure for anchor force: | |||
: a. ksi 3.40 481 1633.5 f kips, 1633.5 P b===<ksi 3.82 f cp= b. ksi 3.88 481 1866.8 f kips, 1866.8 P b=== The depth below the bearing plate, at which concrete stresses are identical to those allowed beyond the anchorage zone area, is computed on the conservative assumption that the load | |||
spreads at a 30 o angle from the plane normal to the bearing plate surface. | |||
Obtain the required concrete area to satisfy the above condition for: | |||
P = 1633.5 kips, A g = 1633.5 = 727 in 2 2.25 | |||
If it is conservatively assumed that the anchor force spreads at an angle of 30 o , the depth (h) below the bearing plate at which the concrete stress will have been reduced to the allowable | |||
compression stress as shown above, can now be computed. See Figure 5D-8b. +=4 7-)30 2htan (20.5 A 2 2 g in. 4.5 h 38.5, -3 2h 20.5 622 2=+= Special provisions must be made to confine the concrete within the 4.5 in. below the bearing plate. From then on, spiral reinforcement m eeting the requirement for tensile splitting stresses will be investigated as described in part 3.3. | |||
===3.2 EQUIVALENT=== | |||
FLUID PRESSURE METHOD | |||
Directly below the bearing plate the necessary reinforcement for confining the concrete is | |||
computed for an assumed equivalent fluid pressure "p". The equivalent fluid pressure is equal | |||
to the concrete compression stress multiplied by the Poissons Ratio, which is assumed to be | |||
0.2. | |||
The concrete compression stress 2 in. below the bearing plate,which is about at the center of | |||
the critical 4.5 in. depth, will be used to determine the equivalent fluid pressure. The initial | |||
1,400 kips anchor force is considered. See Figure 5D-8b. | |||
[]2 2 2 g in 481.5))4 ((7- )30 (tan 4 5.20 (A=+= ksi 2.91 481.5 1400 f c== The equivalent fluid pressure will be: | |||
p = 2.91 (0.2) = 0.58 ksi TMI-1 UFSAR APPENDIX 5D 5D-5 REV. 18, APRIL 2006 The reinforcement is arranged as shown on GAI Drawing 421031 . The #5 spiral has a | |||
diameter of 20 inches with a 3 1/2 in. pitch, and a total length of 36 inches. Horizontal ties | |||
2-#10 and vertical reinforcement 1-#10 plus 1-#8 are also available to resist this equivalent fluid pressure. | |||
The required steel area for the 20 inch diameter spiral was computed for an allowable steel | |||
stress of f s = 20 ksi. | |||
Required steel area: | |||
A s = 0.58(4.5)(20) | |||
= 2.61 in 2 (steel area req'd. over a 4.5 in length below the bearing plate) 20 | |||
Steel provided: | |||
: a. Radially, Spiral 2-#5 = 0.62 in 2 Ties 2-#10 = 2.54 3.16 in 2 > 2.61 in 2 | |||
: b. Vertically, Spiral 2-#5 = 0.62 Vert. 1-#10 = 1.27 Vert. 1-#8 0.79 2.68 in 2 > 2.61 in 2 | |||
===3.3 TRANSVERSE=== | |||
TENSION-PRISM METHOD | |||
Transverse stresses are produced when the anchor force spreads out into the surrounding | |||
concrete. Directly behind the bearing plate the stress lines are parallel to the tendon axis. The | |||
angle change of the stress-flow-lines, where they bend outwards, produces in these regions | |||
compression in the concrete which is of no concern. Further away, the stress-flow-lines tend to | |||
bend back and tensile transverse stresses are created as a result of this angle change. Figure | |||
5D-9 and 5D-10 show the idealized compression stress-flow-lines in plane and the approximate location of transverse stresses. | |||
Further, it can be seen that the angle changes will be small, and that the stress intensity along | |||
the compression stress-flow-lines reduces as they spread out into the concrete mass. From | |||
these considerations it can be concluded that the resulting tensile transverse stresses will be | |||
small and within the tensile capacity of the concrete. Deformed vertical, circumferential, and | |||
radial reinforcement is provided throughout the buttr ess in addition to the spiral reinforcement at each bearing plate. The reinforcement will tie the concrete mass together and ensures that | |||
thermal cracking of the concrete mass, if it should occur, is of no significance as far as structural | |||
integrity of the Reactor Building is concerned. | |||
Figure 5D-10 shows that the critical condition, with respect to transverse tensile stresses in the | |||
vertical plane, exists when only individual tendons ar e stressed. When all tendons are stressed, the stress-flow-lines are more or less parallel. Figure 5D-9 and 5D-10 show, superimposed | |||
over the anticipated stress-flow-lines, "equivalent beams" for which the transverse stresses can be computed. Computation for transverse stresses in radial planes are performed in tabulated | |||
form in Table 5D-1. Three cases are analyzed as follows: | |||
TMI-1 UFSAR APPENDIX 5D 5D-6 REV. 18, APRIL 2006 Case 1 Finds the maximum transverse radial stress in symmetrical equivalent beam. Variation la | |||
investigates anticipated conditions during normal operation. Variations 1b and 1c are made for | |||
extremely shallow and extremely deep "equivalent beams" to establish how variations of beam depths affect the magnitude and location of the transverse stresses. | |||
Case 2 Investigates the effect of eccentric loads on equivalent beams. Guyon's "Symmetrical Prism | |||
Method" is used. According to this method the depth of the eccentric beam does not affect the | |||
magnitude of radial transverse stresses; only the dimension of the prism drawn symmetrical to the tendon center line determines their magnitude. Possible variations in eccentric beam | |||
depths are therefore of no concern. | |||
All stresses were computed on the conservative assumption that the load spreads at 30 o in the circumferential direction. The stress computations confirm the previously made prediction that transverse stresses are within the acceptable concrete stress range. | |||
Case 3 Finds the maximum transverse vertical stresses in symmetrical equivalent beam. Variations 3a | |||
and 3b investigate anticipated conditions during normal operation. Vertical transverse forces | |||
counteract each other under operating conditions largely through vertical compression forces in | |||
the buttress concrete. See Figure 5D-10b. Variation 3c investigates anticipated conditions if | |||
only one anchor force acts. See Figure 5D-10a. | |||
====3.3.1 REINFORCEMENT==== | |||
FOR TRANSVERSE FORCES | |||
The allowable concrete tensile stresses in the radial and vertical planes are not exceeded and | |||
no special transverse reinforcement is required for design loads. However, reinforcement is | |||
available in the form of buttress reinforcement in the hoop and meridional direction consisting of a radial tie-reinforcement anchored into the shell and spiral reinforcement in the maximum | |||
bursting zone. | |||
The buttress is reinforced sufficiently to resist all transverse tensile forces Zy (design) acting | |||
within the limits of the buttress without relying on concrete tensile resistance. | |||
Spiral: 5/8" ?, 3 1/2" pitch, d = 20", length = 36" Maximum tensile splitting stress, Table 5D-1, f by =0.178 ksi | |||
T = 3.5 (0.178) (20) = 6.24 K , f s = 6.24 = 20.1 ksi 2 0.31 | |||
The spiral is designed for Zy (design) = 420 kips | |||
Maximum transverse tensile force as computed is Z y = 311 kips. | |||
TMI-1 UFSAR APPENDIX 5D 5D-7 REV. 18, APRIL 2006 Steel provided: | |||
Spiral 22 - #5 = 6.82 in 2 Ties 10 - #10 = 12.70 19.52 in 2 | |||
f s = 420 = 21.5 ksi 19.52 | |||
Temperature changes in the concrete during the life of the plant or at accident condition will | |||
produce a maximum of less than 1 ksi stress in the bursting zone reinforcement. | |||
===3.4 UNBALANCED=== | |||
FORCES | |||
Figure 5D-11 shows the tendon anchor forces for six buttresses in the Reactor Building. | |||
The solid line horizontal force vector in Figure 5D-11 represents the tangential components of | |||
the anchor forces. The anchor forces are applied on opposite faces of the buttresses and | |||
counteract each other. However, since the vert ical spacing of the anchors varies, it is not possible to balance all the forces. The unbalanced horizontal forces must be transferred from | |||
the buttress into the Reactor Building wall. | |||
The two dashed force vectors, Figure 5D-11, adjacent and parallel to each anchor force vector | |||
are obtained by proportioning the anchor force according to its position with respect to two | |||
opposing anchors similar to determining the react ions of a simply supported beam with an unsymmetrical concentrated load. The dashed horizontal force vector on the buttress center | |||
line indicates the unbalanced force between two opposing anchor forces. How these horizontal | |||
forces are resisted is discussed later. | |||
Since the horizontal anchor forces on opposing sides are offset from each other, the forces | |||
must change direction in order to balance each other. Transverse forces in vertical direction are | |||
a result of this angle change. The vertical dashed force vectors shown in Figure 5D-11 | |||
represent the unbalanced vertical forces behind each anchor zone. | |||
Vertical forces in Figure 5D-11 were computed for the simplified assumption that the anchor | |||
forces travel in straight lines from anchor plate to anchor plate after first moving 18 in. | |||
horizontally into the concrete and gradually bend over into the inclined position. See Figure 5D-12. | |||
The check analysis employing the Finite-Element Me thod did reveal a limited region of surface tensile stresses which presumable would produce cracking. However after concrete cracking | |||
has occurred with the extent limited by s upplemental reinforcing a new state of equilibrium will exist. The following section discussing vertical and horizontal unbalanced forces verifies that a | |||
mechanism does exist to maintain equilibrium. | |||
====3.4.1 VERTICAL==== | |||
UNBALANCED FORCES | |||
During initial prestressing and during tendon inspection the balance of opposing tendon anchors | |||
might not exist and the total vertical transverse stress must be resisted. The maximum concrete | |||
transverse TMI-1 UFSAR APPENDIX 5D 5D-8 REV. 18, APRIL 2006 tensile stress for this condition was computed in Table 5-2, Case 3c, f by = 115 psi. Figure 5D-10a shows the anticipated typical vertical stress-flow-lines for such conditions. The figure | |||
also shows the "equivalent beam" for which transverse stresses were determined according to | |||
Guyon's Symmetrical Prism Method. Sufficient reinforcing steel in form of a spiral has been | |||
provided for Case 3c without relying on concrete tensile resistance. In addition, the vertical | |||
reinforcement in the buttress, consisting of 7-#10, As = 8.8 in.,(Reference 2) can resist a | |||
transverse splitting force of Z = 8.8 (0.95) (40) = 335 kips which is slightly more than the | |||
expected force Zy = 311 kips shown in Table 5D-1. The reinforcement is placed mainly on the buttress surfaces in order to provide access for concrete placement. There is no need to | |||
distribute this reinforcement throughout the cross-se ction of the buttress or to place it where the vertical transverse forces actually act. Wherever the reinforcement is placed, it will not prevent the concrete from cracking and will only tie t he concrete together preventing the development of uncontrolled cracks. Consequently, surface reinforcem ent is sufficient. It restrains the buttress in the vertical direction and no failure mechanism caused by vertical transverse stresses can | |||
develop. | |||
The highest vertical in-plane shear stresses at the buttress-cylinder intersection occur directly | |||
beneath the ring-girder and at the cylinder-base s lab intersection haunch area. An amplification factor of 1.64 was applied to the shear stresses in both areas to account for the increased | |||
stiffness at the buttress. Sufficient reinf orcing was provided between the buttress and cylinder to take the shear stresses. | |||
====3.4.2 HORIZONTAL==== | |||
UNBALANCED FORCES | |||
Unbalanced horizontal forces shown in Figure 5D-11 must be transferred into the wall. | |||
During initial prestressing and tendon inspection the counteracting effect of opposing anchor | |||
forces might be lost at any any anchor zone position. All anchor zones will therefore be reinforced as the analysis requires for such conditions. Figure 5D-9a shows the anticipated | |||
typical horizontal stress-flow-lines for such conditions. The figure also shows the "equivalent beams" for which transverse stresses were determined according to Guyon's "Symmetrical Prism Method." | |||
The analysis for horizontal transverse forces is tabulated in Table 5D-1. The critical loading | |||
condition exists when a tendon is stressed without being counterbalanced by opposing anchor | |||
forces. The maximum concrete transverse tens ile stress for this condition is represented by Case 2a, f by = 178 psi. Sufficient reinforcing steel in form of a spiral has been provided for Case 2a without relying on concrete tensile resistance. In addition, the radial tie reinforcement | |||
in the buttress, consisting of 10-10#, A s = 12.7 in., (Reference 2) can resist a transverse splitting force of Z = 12.7 (0.95) (40) = 482 kips which is larger than the expected force Zy = 311 kips | |||
shown in Table 5D-1. The tie reinforcement restrains the buttress in radial direction and no | |||
failure mechanism caused by radial transverse stresses can develop. | |||
The highest shear stresses for the case of no opposing prestress force in the buttress occurs at | |||
the buttress-cylinder discontinuity on the side of the buttress opposite the existing prestress force. Reinforcing of #10 at 12 in. (See Figure 5D-4) are provided to take the resulting tensile | |||
stresses. | |||
TMI-1 UFSAR APPENDIX 5D 5D-9 REV. 18, APRIL 2006 | |||
====3.4.3 EQUILIBRIUM==== | |||
CONDITION IN THE ANCHOR BLOCK | |||
The reinforcement is checked based on the assumption of an existing spalling crack in the | |||
vertical direction. The forces and the extent of cracking in the anchor block were based on a | |||
method developed by Gergely, Soyen, and Siess of arresting cracks due to transverse stresses (Reference 4). The idealized mathematical model used for this approach, see Figure 5D-13, is | |||
quite conservative. The computation is shown in Table 5D-2. It should be noted that zero | |||
moment exists at B' when c = 27 in., and that the governing moment, Mm, prevails at c = 42 in., that is, far beyond the face of the wall. | |||
Concerning the visco-elastic behavior of concrete, it can be stated that the resulting | |||
redistribution of stresses developed by Gergely will not significantly affect reinforcing requirements (Reference 5). The super-position of bursting stresses and normal stresses from vertical tendons results in a tri-axial state of stress which is less critical than during construction | |||
stages. No additional investigations for the final state of stress in the anchorage zones, therefore, appear necessary. The various sour ces of errors in a numerical tri-axial analysis (uncertainties of material properties, element definition, etc.) may lead to less accurate values | |||
than were obtained by two-dimensional methods. | |||
The shear across the cylinder at the buttress-cylinder intersection for the case of unbalanced | |||
prestressed force has been checked. Sufficient reinforcing had been provided to account for | |||
the shear stresses. For accident temperature and pressure conditions a crack across the entire | |||
cylinder at the buttress cylinder discontinuity has been assumed. Sufficient steel reinforcing is provided to take any tension where the compressive stresses due to prestress have been | |||
overcome. | |||
4 CHECK ANALYSIS BY FINITE-ELEMENT METHOD | |||
==4.1 INTRODUCTION== | |||
The finite-element method has been employed to check the design of tendon anchorage zones | |||
in the Reactor Building. The purpose of this analysis was not as a design basis but simply as a | |||
check to demonstrate the conservatism of the method described here before. This section | |||
presents the method of analysis, load conditions, and results for the buttress and adjacent wall, and in conclusion to correlate these results with the design of the buttress reinforcement which | |||
has been presented in Section 3. | |||
The dominate loads are: | |||
: a. Prestress forces | |||
: b. Thermal loads | |||
: c. Pressure forces due to loss of coolant accident (LOCA) | |||
Although the stress analysis of the buttress and adjacent wall is a three dimensional problem, it can be reduced to a two dimensional problem, as the meridional stress of the cylindrical portion of the vessel is usually constant. The finite-element method of plane stress analysis (References 6 and 7) which can account for complex loadings, arbitrary geometric configurations, and arbitrary boundary conditions, assumes a homogeneous material (i.e., no TMI-1 UFSAR APPENDIX 5D 5D-10 REV. 18, APRIL 2006 cracks) and is an ideal tool for dealing with such a problem. Therefore, a plane stress finite-element computer program (Referenc e 7) was used in the buttress analysis. | |||
For the plane stress finite-element program t he thermal loads were obtained from another computer program. The results from this progr am, for the steady state winter normal operation and the transient state during winter start-up, w inter shut-down, and winter loss of coolant accident were used as input data in the finite- element program. | |||
As described later, a system of 768 direct st iffness simultaneous equations required solution. | |||
To save computer storage, the Gauss-Seidel iter ative method (Reference 8), which requires the storage of only the non-zero elements of the sti ffness matrix, was employed. In order to expedite convergence, an over-relaxation factor (R eference 9) was also used in Reference 7; this factor usually ranges from 1.0 to 2.0. | |||
Cracking of the concrete for the winter LOCA load condition is predicted. Areas of high tensile | |||
stresses have been defined and reinforcement capable of resisting these provided. | |||
===4.2 PLANE=== | |||
STRESS FINITE-ELEMENT COMPUTER PROGRAM | |||
Due to the fact that the meridional stresses of the cylindrical portion of a containment vessel are | |||
essentially uniform and constant, the horizontal cross section of the cylinder may be considered | |||
as a section subjected to plane stresses only. The finite-element method is a versatile method | |||
of analysis and is capable of providing states of stress in the buttress for complex loadings, arbitrary geometric configurations, and arbitrary boundary conditions. Tendon holes are not | |||
included in the buttress geometry. | |||
In the plane stress finite-element analysis, the horizontal cross sectional plane of the cylindrical portion of the vessel is replaced by a finite number of elements interconnected at a finite number of nodal points. Continuity between element s of the system is maintained by requiring that within each element, "lines initially straight remain straight in their displaced positions." | |||
This requirement is satisfied if the strains within each element are assumed to be constant. | |||
Therefore, the stresses within each element are also constant. Distributed loads are replaced | |||
by equivalent nodal forces which act at the three nodes of the element. Based on these | |||
assumptions, it is possible to derive the sti ffness matrix of a typical element, which is an expression for the nodal forces resulting from unit nodal displacements: | |||
{p} = [k] {?} (1) | |||
where {p} is the vector of six nodal forces which corresponds to the vector {?} of six nodal | |||
displacements (Figure 5D-14), and (k) is the stiffness matrix of the element. | |||
The nodal force vector {p} can be expressed as: | |||
{p} = {PX i , Px j , Px k , Py i , Py j , Pyk}) (2) and the nodal displacements vector {?} can be expressed as, | |||
{?} = {u i , u j , u k , v i , v j , v k} (3) which are shown in Figure 5D-14. | |||
TMI-1 UFSAR APPENDIX 5D 5D-11 REV. 18, APRIL 2006 Formulation of Element Stiffness Matrix | |||
Corresponding to the above assumptions, the displacements of a point at (x,y) are written as | |||
linear functions of the coordinates: | |||
u = u i + C 1 x + C 2 y v = v i + C 3 x + C 4 y (4) | |||
The stiffness matrix [k] can be obtained by the following steps: | |||
: a. Substitute the nodal coordinates into Equation(4) to obtain the displacement functions in terms of the nodal displacements. | |||
: b. Substitute results of (a) into the linear strain-displacement equations to obtain the strains in terms of nodal displacements. | |||
: c. Substitute the results of (b) into the stress-strain equation to obtain the stress in terms of nodal displacements. | |||
: d. Find the nodal resultant forces which are statically equivalent to the stresses found in (c) so that the relationships between the nodal forces and the nodal | |||
displacements can be written in a matrix form as shown in Equation (1). | |||
The explicit formulations of the matrices which compose the element stiffness matrix can be found in Reference 7. | |||
Solution of the Total Stiffness Equations After the stiffness matrices are formulated for each of the elements, they are assembled to form a total system stiffness matrix. If the boundary conditions and the loading conditions are properly specified, the tot al stiffness simultaneous equations are readily solvable to give the nodal displacements. Then, the nodal displacements are used to generate the stresses in each element and at each node. The directions and magnitudes of the principal | |||
stresses can also be computed by using the Mohr's circle technique. | |||
Iterative Method and Optimum Over-Relaxation Factor The Gauss-Seidel iterative procedure (Referenc e 8) is applied to solve the stiffness simultaneous equations. The method requires the storage of only the non-zero elements of the | |||
stiffness matrix during the iterative process. | |||
Since the stiffness matrix [k] is positive definite,the iterative method will always converge (Reference 10). The covergence rate is, however, often rather slow. An over-relaxation factor (Reference 9) can be applied which greatly incr eases the convergence rate. The determination of the proper over-relaxation factor is diffi cult and is based mostly on experience and trial and error. | |||
===4.3 GEOMETRY=== | |||
OF BUTTRESSES WITH TENDONS | |||
The pertinent geometry of the Reactor Building as it affects the buttress analysis is shown in Figure 5D-15. In Figure 5D-15a, a sectional view illustrates the vessel with six buttresses. | |||
Figure 5D-15b, shows the arrangement of the hoop tendons on a stretch out of the cylindrical TMI-1 UFSAR APPENDIX 5D 5D-12 REV. 18, APRIL 2006 wall. Figure 5D-15c, shows the dimensions of one buttress and the arrangement of the hoop and vertical tendons. | |||
===4.4 THERMAL=== | |||
LOADS AND TEMPERATURE PROFILES | |||
If the temperature profile in a buttress is known, the difference between the temperature of any point and the reference temperature can be determined. The temperature difference, with the | |||
modulus of elasticity and thermal expansion coeffi cient of the concrete, provides the equivalent constant stress for each finite-element. The constant thermal stress of each element can be | |||
converted to equivalent nodal resultant forces which can be considered as external forces | |||
applied at the nodal points of a non-heated element. | |||
Temperature Profiles The temperature profiles and iso-thermal curves were calculated by a computer program developed by Brunnel (Reference 11) for the s olution of steady state and transient state problems in one, two, or three dimensional Cart esian or cylindrical systems. The governing differential equations were solved by the method of finite differences, using standard relaxation techniques. | |||
The convergence criterion: | |||
T j (step i) - T j (step i-1) < | |||
0.00001 (5) | |||
Y j (step i) | |||
where T j is the temperature at node j, was applied at every node of the difference lattice. To ensure stability of the transient solutions, the time increment was selected according to the | |||
method presented in Reference 11. | |||
The buttress and wall temperatures were determined for the following conditions: | |||
: a. Normal operation in winter | |||
: b. Start-up in winter | |||
: c. Shut-down in winter | |||
: d. Loss of coolant accident during winter | |||
The winter season was chosen because it produces the greatest thermal gradients in the shell | |||
of the Reactor Building. | |||
The mean monthly temperature was chosen due to t he fact that day- night fluctuations had little effect on the interior portion of the concrete. According to Carslow and Jeager, (Reference 12) | |||
the amplitude of the temperature diminishes as: | |||
wP con C p con exp [X | |||
] (6) 2 k con Thus a distance of 12.8 inches from the outer wall of the Reactor Building the effect of the | |||
oscillation is only 0.1 what it is on the outer face, and the middle and inner portion of the vessel | |||
will be unaffected by day-night temperature oscillations. | |||
TMI-1 UFSAR APPENDIX 5D 5D-13 REV. 18, APRIL 2006 For the start-up condition, it was assumed that the inside air temperature of the Reactor Building rose from 60 F to 110 F in one hour, and for shut-down the reverse was used. The | |||
inside film coefficient was taken as 1.490 Btu/hr/sq ft/F for cases a, b, and c above. This was | |||
modified film coefficient, which took into account the steel liner and the still air between the liner | |||
and the concrete: | |||
h m = 1 (7) 1 X steel X air h i + | |||
+ k steel k air Equation 7 neglects the circumferential conduction and heat capacity in the steel liner and air | |||
gap. However, due to the relatively small thi ckness of the liner and the air gap when compared with the thickness of the concrete wall, these effects may be neglected. The preliminary | |||
analyses verified this assumption. The following values, obtained by theoretical analysis (Reference 13), were used in the thermal analysis: | |||
inside film coefficient = 2.0 Btu/hr/sq ft/F | |||
specific heat, C p con = 0.210 Btu/lb/F | |||
thermal conductivity, k con = 10.0 Btu in./hr/sq ft/F | |||
k steel = 300.0 Btu in./hr/sq ft/F | |||
k air = 1.836 Btu in./hr/sq ft/F | |||
The density of the concrete, P con was assumed to be 140 lbs/cu ft. | |||
The above values agree closely with the results from experimental research carried out in Japan (Reference 14). | |||
C p con = 0.200 Btu/lb/F | |||
k con = 11.28 Btu in./hr/sq ft/F | |||
The density of the concrete P con , was assumed to be 146 lb/cu ft. | |||
The thickness of the steel and air gap was: | |||
X air = 1/32 inch | |||
X steel = 3/8 inch | |||
Using Reference 15, the local film coefficient for air at 20 F and a free stream velocity of 20 mph | |||
flowing past a concrete cylinder was calculated to be: | |||
h = 0.56 Btu/hr/ft 2/F As a safe approximation the outside film coefficient was chosen as: | |||
TMI-1 UFSAR APPENDIX 5D 5D-14 REV. 18, APRIL 2006 h? = 1.0 Btu/hr/sq ft/F | |||
The buttress model was replaced by a mesh of 61 nodes, with adiabatic boundaries at the | |||
centerline of the buttress and at the radial mid-section of the wall between two neighboring | |||
buttresses (Figures 5D-16 and 5D-17). Due to symmetry, there is no heat transferred across | |||
these boundaries. For Case d the liner temperature during LOCA was used as the inner | |||
boundary temperature in conjunction with a film coefficient of 5.880 Btu/hr/sq ft/F, which | |||
simulated the conduction through the air gap: | |||
k air h m = (8) X air The transient state temperature analyses were performed at various time steps. | |||
The temperature profiles for Cases a, b, c, and d are shown in Figure 5D-18. Because the most | |||
severe temperature gradients occur for case a and d, only these two cases were considered in | |||
the stress analyses. Iso-thermal curves fo r cases a and d are shown in Figures 5D-19 and 5D-20. | |||
The design pressure inside the Reactor Building during LOCA is equal to 55 psig and the | |||
maximum pressure was initially calculated to occur at approximately 200 seconds after initiation | |||
of the accident. The maximum liner temperature of 280F was initially calculated to occur at | |||
approximately 650 seconds after the initiation of the accident. However, the maximum | |||
temperature gradient, which causes the most severe thermal stresses, was calculated to occur | |||
at 10,000 seconds after initiation of the accident. | |||
Consequently, the determination of the most severe stress state, due to thermal and pressure loads, requires that a number of stress analyses be conducted at various elapsed times after | |||
commencement of the accident. A conservative estimate of the stresses was obtained, by | |||
superimposing directly the maximum values for accident temperature and pressure. | |||
Consequently although subsequent calculations alter t he time at which minimum pressure and liner temperature occur, these changes do not affect this analysis. | |||
4.5 NON THERMAL LOADS | |||
In addition to the thermal loads, there are other loads which act on the structure: | |||
: a. Hoop Tendon Anchor Pressure | |||
The tendon force was assumed to be distributed uniformly up the face of the buttress. The resulting average force (=p a) was applied to each of six nodes. | |||
The location of the nodes approximately corresponded to the position of the | |||
tendon. | |||
: b. Hoop Tendon Pressure in a Radial Direction | |||
Radial pressures, due to the hoop tendons in the wall, were calculated as average stresses in pounds/in. | |||
The two radial pressures applied (refer to Figure 5D-21): | |||
TMI-1 UFSAR APPENDIX 5D 5D-15 REV. 18, APRIL 2006 | |||
: 1. From point A to the tangent point T. Over this length the tendon spacing causing radial pressure (i.e., continuous past the buttress) is 33 in., giving | |||
a pressure of p t , and | |||
: 2. For the length from point T to point B (midway between buttresses) the average tendon spacing causing radial pressure is 16 1/2 in., giving a | |||
pressure of 2p | |||
: t. | |||
: c. Internal Pressure Due to LOCA | |||
Internal pressure acting on the inside face of the Reactor Building Wall was considered as a combination of two separate effects: | |||
: 1. Pressure due to LOCA. The maximum pressure due to LOCA was taken at 82.5 psig. (1.5 P) | |||
: 2. Pressure due to expansion of the liner caused by the increase in temperature during the LOCA condition (110F to 280F) based upon a limiting yield strength of 1.2 f | |||
: y. | |||
The above two loads were assumed to act simultaneously during the LOCA condition, giving an equivalent internal pressure of p | |||
: i. | |||
The magnitude and combination of the above non-thermal loads is noted in section 4.7 | |||
===4.6 MATERIAL=== | |||
PROPERTIES | |||
Prestressing Steel Wire | |||
The type of prestressing wire to be used in the construction is a stress relieved wire. | |||
Dimensions, properties, etc., of the wire as used in the analysis were: | |||
: a. Size of tendon = 169 wires | |||
: b. Diameter of wires = 1/4 in. | |||
: c. Minimum guaranteed ultimate tensile stress = 240,000 psi | |||
: d. Tendon force @ 70 percent GUTS = 1400 kips | |||
: e. Modulus of elasticity = 29 x 10 6 psi f. Coefficient of friction = 0.16 | |||
: g. Coefficient of wobble friction = 0.00030 TMI-1 UFSAR APPENDIX 5D 5D-16 REV. 18, APRIL 2006 Concrete The assumed properties of the concrete used in the analysis were: | |||
: a. Compressive strength @ 28 days = 5000 psi | |||
: b. Modulus of elasticity = 4 x 10 6 psi | |||
: c. Poissons ratio = 0.20 | |||
Liner The liner is fabricated from 3/8 inch thick steel plate manufactured to ASTM A-283 Grade C. | |||
The minimum yield stress for this material is 30,000 psi. In determining the equivalent liner pressure LOCA, the yield stress was assumed to be 36,000 psi. | |||
===4.7 LOADING=== | |||
CASES AND COMBINATIONS | |||
The buttress and wall were analyzed for four loading cases. The first two cases considered the | |||
buttress and wall being subjected to the maximum prestress forces, i.e., prestress forces | |||
allowing only those losses that would occur between prestressing and start-up (beginning of | |||
plant life.) The latter two cases considered the buttress and wall being subjected to the | |||
minimum prestress forces, i.e., prestress forc es after 40 year losses (end of plant life). | |||
The load cases were: | |||
: a. Winter normal operation with maximum prestress force. | |||
: b. Winter LOCA with maximum prestress force. | |||
: c. Winter normal operation with minimum prestress force. | |||
: d. Winter LOCA with minimum prestress force. | |||
These cases were considered to be more critical than the summer normal operation/LOCA | |||
conditions. The summer condition results in a smaller temperature differential in the wall which | |||
in turn results in maximum compressive and tensile stresses less than the winter condition. | |||
Hence, the summer condition was not considered in the analyses. | |||
TMI-1 UFSAR APPENDIX 5D 5D-17 REV. 18, APRIL 2006 The numerical value of the applied loads was: | |||
LOAD CASE a b c d | |||
t LT 110 150 110 150 | |||
t LP 110 280 110 280 | |||
t e 20 20 20 20 | |||
p a 6.25 6.73 5.52 6.03 | |||
pt A-T 46.07 49.75 39.25 42.2 | |||
2ptT-B 92.13 99.50 78.5 84.4 | |||
p i 0.0 96.1 0.0 96.1 | |||
where: | |||
t LT = Temperature (F) of liner for thermal gradient in wall. | |||
t LP = Temperature (F) of liner for equivalent pressure of liner due to expansion during LOCA. | |||
t e = External air temperature (F) | |||
p a = Prestress force at bearing plate (kips/inch) | |||
pt A-T = Radial prestress force from point A to tangent point T (lbs/inch) | |||
2p T-B = Radial prestress force tangent point T to point B (lbs/inch) | |||
p i = Equivalent internal pressure (LOCA pressure + liner expansion pressure lbs/inch) | |||
Figure 5D-21 shows the application of these loads on the buttress and wall. | |||
Load p a p t , and 2p t for Case a incorporate losses due to elastic shortening of the concrete, friction and creep, shrinkage, and steel relaxation losses that occur between prestressing and | |||
start-up. Case b prestress forces reflect increases from Case a, due to elastic elongation of the | |||
structure during LOCA. Thus, Case a and b reflect the maximum prestress force that will exist | |||
in the structure, at the earliest time in the life of the plant, that a LOCA could occur. | |||
Loads p a , p t , and 2 pt for Cases c and d are similar to Cases a and be respectively, except that the losses due to creep, shrinkage, and steel relaxation are those that are assumed to occur | |||
over 40 years to end of plant life. Case d pres tress forces reflect an increase from Case c due to elastic elongation of the structure during LOCA. | |||
TMI-1 UFSAR APPENDIX 5D 5D-18 REV. 18, APRIL 2006 | |||
===4.8 STRESS=== | |||
ANALYSIS | |||
Due to the symmetrical nature of the structur e only one-twelfth of the cylinder, as shown in Figure 5D-16, need be analyzed. Preliminary analy ses showed that beyond a certain point in the wall, the stress state, for all loading cases, was uniform. To conserve computer storage and | |||
give a more accurate solution, it was decided to reduce the size of the model to that shown in | |||
Figure 5D-22. Nodal displacements of both boundaries were specified. The displacements | |||
were obtained form a computer program as developed by Kalnins (Reference 16) using a shell of revolution theory. | |||
There were a total of 672 elements and 384 nodal points, so that 768 simultaneous equations | |||
required solution. For the computations, the Gauss-Seidel iterative solution procedure was used. | |||
===4.9 ANALYSIS=== | |||
OF RESULTS | |||
Figures 5D-23 through 5D-30 show iso-stress curves for the maximum and minimum principal stresses for the four loading conditions outlined in Section 4.7. It should be noted that these | |||
curves represent the algebraic maximum and minim um values of the principal stresses. The do not show the orientation of the stress trajectories. | |||
Inspection of the iso-stress curves shows tha t, as expected, high compressive stresses exist directly beneath the bearing plate for all load conditions. Under winter normal conditions, a | |||
region of tensile bursting stress occurs in the buttress at some distance removed from the | |||
bearing plate; during LOCA, however, these bursting stresses merge with those resulting from the pressure and temperature loads. | |||
The variation of temperature through both the wall and buttress leads to substantial hoop stress | |||
gradients. It should be noted that in the curves for LOCA, the direction of the stresses in the | |||
wall changes in the vicinity of the centerline. | |||
Of primary interest, however, is the high localized t ensile stress that occurs at the buttress-wall re-entrant corner during LOCA, which results fr om both the thermal gradient as well as the local bending arising from the variation of resistance of the buttress and wall to internal pressure. | |||
4.10 CONCLUSIONS | |||
The results indicate that for the most severe loading conditions, the maximum stresses in the | |||
buttress and adjacent wall are within safe limits, except for the localized area at the buttress | |||
wall corner where tensile stresses greater than 6 ?f | |||
? occur. Reinforcement to stabilize cracking in this local region are provided. All other results would indicate that the design method | |||
employed was conservative. | |||
TMI-1 UFSAR APPENDIX 5D 5D-19 REV. 18, APRIL 2006 5 REFERENCES | |||
: 1. Building Code Requirements for Reinforced Concrete , ACI 318-63, American Concrete Institute, Detroit, Michigan, 1963. | |||
: 2. Y. Guyon, Prestressed Concrete , Contractors Record Ltd. London and John Wiley and Sons Inc., New York, 1960. | |||
: 3. F. Leonhardt, Prestressed Concrete Design and Construction , Wilhelm Ernest and Sohn, Berlin-Munich, 1964. | |||
: 4. P. Gergely, M. A. Sozen and C. P. Siess, The Effect of Reinforcement on Anchorage Zone Cracks in Prestressed Concrete Members , Civil Engineering Studies, Structural Research Series No. 271, University of Illinois, Urbana, Illinois, July 1963. | |||
: 5. Anchorage Systems in Prestressed Concrete Pressure Vessels; Anchorage Zone Problems , Oak Ridge National Laboratory, ORNL-TM-2378, January, 1969. | |||
: 6. M. J. Turner, R. W. Clough, H. C. Martin, and L. J. Topp, "Stiffness and Deflection Analyses of Complex Structures," Jl. of Aero. Sci., V. 23, No. 9, Sept 1965, p. 805. | |||
: 7. E. L. Wilson, Finite-Element Analysis of Two Dimensional Structures , Ph. D. Thesis, Univ. of Cal. at Berkeley, June 1963. | |||
: 8. M. G. Salvadori and M. L. Baron, Numerical Methods in Engineering , Prentice-Hall, Englewood Cliffs, N. J., 1961. | |||
: 9. F. G. Lehan, "Simultaneous Equations S olved by Over-Relaxization, "Proc. ASCE, 2nd Conf. on Electronic Computation , Pittsburgh, Pa., USA, Sept. 1960. | |||
: 10. S. H. Crandall, Engineering Analysis , McGraw-Hill Book Co., N. Y., USA, 1956. | |||
: 11. R. D. Brunnel, Computer Program Heating and Heating 2 , Atomic International Division, Report No. A1-64-MEMO-177, North American Aviation, Inc., Detroit, Michigan, USA, 1964. | |||
: 12. H. S. Craslaw and J. C. Jaeger, Conduction of Heat in Solids , 2nd Ed., Oxford, London, 1959. | |||
: 13. Johns-Manville Research and Engineering Center Report E-455-T-266 filed in Robert Emmett Ginna Final Facility Description and Safety Analysis Report, Section 5, Appendix 5B. | |||
: 14. J. Kosaka, Experimental Research of Unsteady-State Heat Transfer Concerned with Nuclear Energy Facilities (Preliminary Report), Division of Nuclear Energy Facilities, Taisei Construction Ltd., 1969. | |||
: 15. F. Kreith, Principles of Heat Transfer , 2nd Ed. International, Scranton, Pa., 1965. | |||
: 16. A. Kalnins, Computer Programs , Crystal River Unit No., 3, Florida Power Corporation, PSAR, Pages 5C-1 and 5C-2. | |||
: 17. J. Zielinski and R. E. Rowe, An Investigation of the Stress Distribution in the Anchorage Zones of Post-Tensioned Concrete Members , Cement and Concrete Association, Research Report NO. 9, London, September 1960. | |||
TMI-1 UFSAR APPENDIX 5D 5D-20 REV. 18, APRIL 2006 | |||
: 18. J. Zielinski and R. E. Rowe, The Stress Distribution Associated with Groups of Anchorages in Post-Tensioned Concrete Members , Cement and Concrete Association, Research Report No. 13, London, October 1962. | |||
: 19. P. Gergely, Anchorage Systems in Prestressed Concrete Pressure Vessels; Anchorage Zone Problems , U.S. Atomic Energy Commission, ORNL-TM-2378, January, 1969. | |||
TMI-1 UFSAR APPENDIX 5D 5D-21 REV. 18, APRIL 2006 TABLE 5D-1 BUTTRESS ANCHOR ZONE TRANSVERSE STRESSES IN RADIAL AND VERTICAL PLANES Case Radial Plane | |||
1a 1b 1c 2a | |||
2a'(in) 20.5 20.5 20.5 20.5 | |||
2a (in) 52 26 78 32 | |||
2a'/2a 0.40 0.79 0.26 0.64 | |||
fy[1] 0.27 0.10 0.34 0.15 | |||
z (in) [2] 0.80 0.90 0.70 0.9 | |||
2b' (in) 20.5 20.5 20.5 20.5 | |||
2b=2b'+(2az)tg30 o 44.5 33.5 51.5 37.0 A=2a2b (in | |||
: 2) 2315 870 4150 1182 P (kips) [3] 1400 1400 1400 1400 | |||
f bz = P/A (ksi) 0.605 1.610 0.338 1.182 f by = f bz fy (ksi) 0.164 0.161 0.115 0.178 [4] | |||
Zy = 0.3P (1-a'/a) 252 88 311 151 (kips) [5] | |||
Zy (design) = 0.3P 420 420 420 420 (kips) [6] | |||
az (in) 20.8 11.7 27.3 14.4 | |||
[1] Guyon, "Prestressed Concrete", Figure 9. | |||
2 [2] Guyon, "Prestressed Concrete", Figure 9. | |||
2 [3] 1600 Kips and 2000 kips anchor forces is not critical. | |||
[4] Changes in tensile splitting stresses with varying "2a" dimension. | |||
[5] Zy = total transverse force Leonhardt "Prestressed Concrete" | |||
: p. 271. | |||
3 [6] For design purpose "(a'/a)" is conservatively assumed to be equal | |||
to zero. | |||
TMI-1 UFSAR APPENDIX 5D 5D-22 REV. 18, APRIL 2006 TABLE 5D-1 (continued) | |||
BUTTRESS ANCHOR ZONE TRANSVERSE STRESSES IN RADIAL AND VERTICAL PLANES | |||
Case Ultimate Load- Vertical Planes Radial Plane 2a 3a 3b 3c | |||
2a' (in) 20.5 20.5 20.5 20.5 | |||
2a (in) 32 33 26 78 | |||
2a'/2a 0.64 | |||
fy [1] 0.15 | |||
z (in) [2] 0.9 | |||
2b' (in) 20.5 | |||
2b=2b' + (2az)tg30 0 37.0 See Case 2a, 1b, and 1c A=2a2b (in | |||
: 2) 1182 P (kips) [3] 2000 | |||
f bz = P/A (ksi) 1.690 f by = f bz fy (ksi) 0.254 [4] | |||
Zy=0.3P (1-a'/a) 216 (kips) [5] | |||
Zy (design)=0.3P 600 (kips) [6] | |||
az (in) 14.4 | |||
[1] Guyon, "Prestressed Concrete", Figure 9. | |||
2 [2] Guyon, "Prestressed Concrete", Figure 9. | |||
2 [3] 1600 Kips and 2000 kips anchor forces is not critical. | |||
[4] Changes in tensile splitting stresses with varying "2a" dimension. | |||
[5] Zy = total transverse force Leonhardt "Prestressed Concrete" p. 271. | |||
3 [6] For design purpose "(a'/a)" is conservatively assumed to be equal to zero. | |||
TMI-1 UFSAR APPENDIX 5D 5D-23 REV. 18, APRIL 2006 TABLE 5D-2 BUTTRESS ANCHOR ZONE SPALLING (Refer to Figure 5D-13) | |||
M = P {c - e - (c/h) 2 [2h - 3e - c + 2ec | |||
]} h c= varies e = 15.5 in. h = 69 in. | |||
M = P {c - 15.5 - (c/69) 2 [2 x 69 - 3 x 15.5 - c ( 1 - 2x15.5 | |||
)]} 69 | |||
C (in.) c - 15.5 (c/69) 2 -c( 1 - 2 x 15.5) M/P M m (in. - K) 69 | |||
5 -10.5 0.00525 - 2.75 -11.0 10 - 5.5 0.02100 - 5.50 - 7.3 15.5 0 0.05040 - 8.50 - 4.4 20 + 4.5 0.08400 -11.00 - 2.3 25 + 9.5 0.1310 -13.80 -0.7 zero at 27 in. 30 +14.5 0.1890 -16.50 + 0.3 42 +26.5 0.370 -23.10 + 1.2 1680 50 +34.5 0.525 -27.50 + 0.9 (in. - kips) 69 +53.5 1.00 -38.00 0 | |||
F t = M m Z = 3 in h = 69 in. | |||
h-z | |||
f s allowable = 10 5 W A s A s = 1.27 sq in. (#10 bars) | |||
W = 0.005 | |||
f s = 10 5 0.005/1.27 = 6.28 ksi F T = 1680 | |||
= 25.5 K 69 - 3 | |||
A s = F T/f s = 25.5/6.28 = 4.05 sq in. | |||
A s provided: #10 @ 12 in. + 1 #10 Loop = 6.03 sq in.}} |
Latest revision as of 13:50, 9 October 2018
ML16127A493 | |
Person / Time | |
---|---|
Site: | Three Mile Island |
Issue date: | 04/11/2016 |
From: | Exelon Generation Co |
To: | Office of Nuclear Reactor Regulation |
Shared Package | |
ML16127A520 | List:
|
References | |
TMI-16-035 | |
Download: ML16127A493 (28) | |
Text
TMI-1 UFSAR APPENDIX 5D 5D-i REV. 18, APRIL 2006
APPENDIX 5D REPORT ON STRESS ANALYSIS OF BUTTRESS
TMI-1 UFSAR APPENDIX 5D 5D-ii REV. 18, APRIL 2006 TABLE OF CONTENTS
SECTION TITLE 1 DESIGN DETAILS
1.1 TENDON
ANCHORAGE HARDWARE
1.2 ANCHORAGE
ZONE REINFORCEMENT
2 TENDON LOADS AND ALLOWABLE STRESSES
2.1 TENDON
LOADS
2.2 ALLOWABLE
STRESSES
3 ANCHOR BLOCK STRESSES
3.1 BEARING
STRESSES
3.2 EQUIVALENT
FLUID PRESSURE METHOD
3.3 TRANSVERSE
TENSION-PRISM METHOD
3.3.1 REINFORCEMENT
FOR TRANSVERSE FORCES
3.4 UNBALANCED
FORCES
3.4.1 VERTICAL
UNBALANCED FORCES
3.4.2 HORIZONTAL
UNBALANCED FORCES
3.4.3 EQUILIBRIUM
CONDITION IN THE ANCHOR BLOCK
4 CHECK ANALYSIS BY FINITE-ELEMENT METHOD
4.1 INTRODUCTION
4.2 PLANE
STRESS FINITE-ELEMENT COMPUTER PROGRAM
4.3 GEOMETRY
OF BUTTRESSES WITH TENDONS
4.4 THERMAL
LOADS AND TEMPERATURE PROFILES 4.5 NON THERMAL LOADS
4.6 MATERIAL
PROPERTIES
4.7 LOADING
CASES AND COMBINATIONS
4.8 STRESS
ANALYSIS
4.9 ANALYSIS
OF RESULTS
4.10 CONCLUSIONS
5 REFERENCES
TMI-1 UFSAR APPENDIX 5D 5D-iii REV. 18, APRIL 2006 LIST OF TABLES
TABLE TITLE 5D-1 Buttress Anchor Zone Transverse Stresses in Radial and Vertical Planes
5D-2 Buttress Anchor Zone Spalling
TMI-1 UFSAR APPENDIX 5D 5D-iv REV. 18, APRIL 2006 LIST OF FIGURES
FIGURE TITLE 5D-1 DELETED
5D-2 DELETED
5D-3 DELETED
5D-4 DELETED
5D-5 DELETED
5D-6 DELETED
5D-7 Typical Hoop Tendon Arrangement
5D-8 Buttress-Hoop Tendon Bearing Stresses
5D-9 Buttress-Radial Bursting Stresses
5D-10 Buttress-Vertical Bursting Stresses
5D-11 Hoop Tendon Typical Anchor Force Distribution
5D-12 Hoop Prestress Idealized Vertical Stress Flow in Buttress
5D-13 Idealized Model for Spalling
5D-14 A Triangular Plane Stress Finite-Element
5D-15 Reactor Building With Buttresses
5D-16 Adiabatic Sections
5D-17 Model of Gridwork for Thermal Analysis
5D-18 Temperature Profiles Through Wall and Buttress for Various Thermal Load Conditions
5D-19 Iso-Thermal Curves for Winter Normal Operating Temperature Condition
5D-20 Iso-Thermal Curves for Winter A ccident Temperature Condition at Time = 10,000 Secs.
5D-21 The Hoop Tendon Loads, Internal Pressure and the Equivalent Liner Pressure
5D-22 Finite-Element Grid Work for Buttress Analysis TMI-1 UFSAR LIST OF FIGURES
FIGURE TITLE APPENDIX 5D 5D-v REV. 18, APRIL 2006 5D-23 Minimum Principal Stresses for Maximum Prestress Forces and Winter Normal Operating Condition (Case A) (PSI)
5D-24 Maximum Principal Stresses for Maximum Prestress Forces and Winter Normal Operating Condition (Case A) (PSI)
5D-25 Minimum Principal Stresses for Maximum Prestress Forces and Winter Accident Condition (Case B) (PSI)
5D-26 Maximum Principal Stresses for Maximum Prestress Forces and Winter Accident Condition (Case B) (PSI)
5D-27 Minimum Principal Stresses for Minim um Prestress Forces and Winter Normal Operating Condition (Case C) (PSI)
5D-28 Maximum Principal Stress for Minimum Prestress Condition (Case C) (PSI)Forces and Winter Normal Operating
5D-29 Minimum Principal Stresses for Minim um Prestress Condition (Case D) (PSI)
Forces and Winter Accident
5D-30 Maximum Principal Stresses for Minimum Prestress Forces and Winter Accident Condition (Case D) (PSI)
TMI-1 UFSAR APPENDIX 5D 5D-1 REV. 18, APRIL 2006 1 DESIGN DETAILS
1.1 TENDON
ANCHORAGE HARDWARE
The prestressing system employs the BBRV anchor age system utilizing 169-1/4" diameter wires conforming to ASTM A421-65, Type BA. The dimensions and capacities of the tendons used
hereafter for the anchorage zone analysis are more fully described in Appendix 5B, Part D.
1.2 ANCHORAGE
ZONE REINFORCEMENT
The mild steel reinforcement in the tendon anchorage zones of the buttress and the base mat
are shown on GAI Drawings 421004, 421005, 421006, 421031,421032, and 421033. This
reinforcement was conservatively designed by the methods enumerated hereafter in Section 3 of this Appendix and is judged to be sufficient to accommodate all predicted stress states
without jeopardizing the ability to satisfactorily place the concrete.
2 TENDON LOADS AND ALLOWABLE STRESSES
2.1 TENDON
LOADS
The minimum guaranteed ultimate tensile strength of the tendons is approximately 2,000 kips.
Before the tendons are anchored (locked off) at 70 percent of the ultimate tendon force, which is
equal to 1,400 kips, they may be temporarily ov erstressed to 1,600 kips in order to reduce friction losses. The anchor zones are designed for anchor forces of 1,400 kips per tendon. The
anchor zones are also checked to resist safely temporary anchor force of 1,600 kips per tendon.
2.2 ALLOWABLE
STRESSES
The allowable stresses for concrete and steel will conform to those given in Section 5.2.2 and
5.2.3 of the FSAR.
Due to the specialized nature of anchor zone design the following conditions and stress criteria
will be used in the design. These stress criteria include adequate conservatism in their
application.
First, verify that the allowable concrete compression and tensile stresses are not exceeded due
to bearing and bursting stresses respectively.
Second, check the tensile stresses in reinforcing steel when no credit is taken for concrete
tensile capacity in the investigated anchor zones.
Under these conditions, the allowable stresses will be shown as follows:
TMI-1 UFSAR APPENDIX 5D 5D-2 REV. 18, APRIL 2006 CONCRETE
- 1. Allowable bearing stress
according to ACI/318-63, 1 Paragraph 2605(C)
- 2. Allowable compression stress in the structure beyond the anchorage zone c c f 0.45 f=
- 3. Allowable Tension Stress
REINFORCING STEEL
Allowable tension for
A615-Grade 40, f s = 20 ksi.
3 ANCHOR BLOCK STRESSES The investigation of the buttress anchorage zones for the hoop tendons is described in the following sections. The analytical work as presented, is based on the end anchor zone stresses
and stress flow pattern studies documented by Y. Guyon (Reference 2) and F. Leonhardt (Reference 3). Their work has been substant iated through extensive strain measurements and stress flow pattern established by studies utilizing photoelastic techniques.
The design presented herein is quite conservative; nevertheless, it has been adopted as the
basis for determination of required reinforcing steel in the end anchor zones to overcome tensile
splitting stresses and to provide for spalling stresses near the concrete surface.
The existance of a three-dimensional stress field including high shear, thermal and shrinkage
tensile stresses is influenced by hairline cracking due to creep and thermal stresses. In order to
accommodate the most unfavorable stress combination in the anchorage zones, the following
design goals were established:
- 1. Maximum concrete bursting tensile stresses have been kept below maximum allowable tensile stress of conservatively c f 3using Y. Guyon's symmetrical prism method. This would indicate that no special spiral and bursting reinforcement is required. The 3 b b ci cp/A A' f 6.0 f=ksi. 5.0 f ci=ksi. 5.0 f c=c f 3ksi. 5.0 f c=
TMI-1 UFSAR APPENDIX 5D 5D-3 REV. 18, APRIL 2006 combined stresses during normal operation and test loading also indicate that the concrete will remain uncracked in the anchorage zone.
- 2. However, due to the existance of complex thermal stresses and possible microcracking, the ability of concrete to carry tensile stresses has been discarded, and mild steel
reinforcing has been provided for:
- a. the equivalent fluid pressure
- b. transverse tensile stresses
- c. unbalanced vertical forces
- d. unbalanced horizontal forces
- e. equilibrium condition
The combined stresses during the design pressure for the Reactor Building (p = 55 psig) indicate that the concrete will remain uncracked in the anchorage zone. The combined
stresses during the factored accident loading (refer to 5.2.3.2, cases a. and b. indicate
that the concrete in the anchorage zone may develop hairline cracks.
- 3. Existing theoretical solutions for anchorage zone design do not agree with more recent experimental results (References 17 and 18); maximum transverse tensile forces were
made up to 5.8 times larger than determined by Y. Guyon.
- 4. Possible failure mechanisms due to spalling and bursting stresses were evaluated to verify static equilibrium and acceptable stre ss limits under such hypothetical conditions.
- 5. Bursting stresses in a large body are much smaller than in a member bounded by planes (Reference 19). Therefore, while it was found that two-dimensional theories underestimate
the bursting stresses (Reference 17), analytical results based on isolated prisms or better, experimental values such as those developed by J. Zielinski and R. E. Rowe can be used for massive structures (Reference 19).
- 6. The design of the anchorage zones provides a margin of safety similar to the entire reactor building. Furthermore, the reinforcement is designed to control cracking due to
spalling, bursting, or shear of the concrete under sustained load.
3.1 BEARING
STRESSES
Figures 5D-7 through 5D-10 show the typical hoop tendon arrangement and buttress anchor
zone.
Bearing area:
2 2 2 b in. 481 4 11-24 A== Concrete anchor surface:
2 2 2 b in. 994 4 11-33 A= Allowable bearing stress:
()ksi 3.82 1/3 481 994 5.0 0.6 f cp== Use of 28 day concrete strength is conservative since initial tensioning will occur significantly
later in time.
TMI-1 UFSAR APPENDIX 5D 5D-4 REV. 18, APRIL 2006 Bearing pressure for anchor force:
- a. ksi 3.40 481 1633.5 f kips, 1633.5 P b===<ksi 3.82 f cp= b. ksi 3.88 481 1866.8 f kips, 1866.8 P b=== The depth below the bearing plate, at which concrete stresses are identical to those allowed beyond the anchorage zone area, is computed on the conservative assumption that the load
spreads at a 30 o angle from the plane normal to the bearing plate surface.
Obtain the required concrete area to satisfy the above condition for:
P = 1633.5 kips, A g = 1633.5 = 727 in 2 2.25
If it is conservatively assumed that the anchor force spreads at an angle of 30 o , the depth (h) below the bearing plate at which the concrete stress will have been reduced to the allowable
compression stress as shown above, can now be computed. See Figure 5D-8b. +=4 7-)30 2htan (20.5 A 2 2 g in. 4.5 h 38.5, -3 2h 20.5 622 2=+= Special provisions must be made to confine the concrete within the 4.5 in. below the bearing plate. From then on, spiral reinforcement m eeting the requirement for tensile splitting stresses will be investigated as described in part 3.3.
3.2 EQUIVALENT
FLUID PRESSURE METHOD
Directly below the bearing plate the necessary reinforcement for confining the concrete is
computed for an assumed equivalent fluid pressure "p". The equivalent fluid pressure is equal
to the concrete compression stress multiplied by the Poissons Ratio, which is assumed to be
0.2.
The concrete compression stress 2 in. below the bearing plate,which is about at the center of
the critical 4.5 in. depth, will be used to determine the equivalent fluid pressure. The initial
1,400 kips anchor force is considered. See Figure 5D-8b.
[]2 2 2 g in 481.5))4 ((7- )30 (tan 4 5.20 (A=+= ksi 2.91 481.5 1400 f c== The equivalent fluid pressure will be:
p = 2.91 (0.2) = 0.58 ksi TMI-1 UFSAR APPENDIX 5D 5D-5 REV. 18, APRIL 2006 The reinforcement is arranged as shown on GAI Drawing 421031 . The #5 spiral has a
diameter of 20 inches with a 3 1/2 in. pitch, and a total length of 36 inches. Horizontal ties
2-#10 and vertical reinforcement 1-#10 plus 1-#8 are also available to resist this equivalent fluid pressure.
The required steel area for the 20 inch diameter spiral was computed for an allowable steel
stress of f s = 20 ksi.
Required steel area:
A s = 0.58(4.5)(20)
= 2.61 in 2 (steel area req'd. over a 4.5 in length below the bearing plate) 20
Steel provided:
- a. Radially, Spiral 2-#5 = 0.62 in 2 Ties 2-#10 = 2.54 3.16 in 2 > 2.61 in 2
- b. Vertically, Spiral 2-#5 = 0.62 Vert. 1-#10 = 1.27 Vert. 1-#8 0.79 2.68 in 2 > 2.61 in 2
3.3 TRANSVERSE
TENSION-PRISM METHOD
Transverse stresses are produced when the anchor force spreads out into the surrounding
concrete. Directly behind the bearing plate the stress lines are parallel to the tendon axis. The
angle change of the stress-flow-lines, where they bend outwards, produces in these regions
compression in the concrete which is of no concern. Further away, the stress-flow-lines tend to
bend back and tensile transverse stresses are created as a result of this angle change. Figure
5D-9 and 5D-10 show the idealized compression stress-flow-lines in plane and the approximate location of transverse stresses.
Further, it can be seen that the angle changes will be small, and that the stress intensity along
the compression stress-flow-lines reduces as they spread out into the concrete mass. From
these considerations it can be concluded that the resulting tensile transverse stresses will be
small and within the tensile capacity of the concrete. Deformed vertical, circumferential, and
radial reinforcement is provided throughout the buttr ess in addition to the spiral reinforcement at each bearing plate. The reinforcement will tie the concrete mass together and ensures that
thermal cracking of the concrete mass, if it should occur, is of no significance as far as structural
integrity of the Reactor Building is concerned.
Figure 5D-10 shows that the critical condition, with respect to transverse tensile stresses in the
vertical plane, exists when only individual tendons ar e stressed. When all tendons are stressed, the stress-flow-lines are more or less parallel. Figure 5D-9 and 5D-10 show, superimposed
over the anticipated stress-flow-lines, "equivalent beams" for which the transverse stresses can be computed. Computation for transverse stresses in radial planes are performed in tabulated
form in Table 5D-1. Three cases are analyzed as follows:
TMI-1 UFSAR APPENDIX 5D 5D-6 REV. 18, APRIL 2006 Case 1 Finds the maximum transverse radial stress in symmetrical equivalent beam. Variation la
investigates anticipated conditions during normal operation. Variations 1b and 1c are made for
extremely shallow and extremely deep "equivalent beams" to establish how variations of beam depths affect the magnitude and location of the transverse stresses.
Case 2 Investigates the effect of eccentric loads on equivalent beams. Guyon's "Symmetrical Prism
Method" is used. According to this method the depth of the eccentric beam does not affect the
magnitude of radial transverse stresses; only the dimension of the prism drawn symmetrical to the tendon center line determines their magnitude. Possible variations in eccentric beam
depths are therefore of no concern.
All stresses were computed on the conservative assumption that the load spreads at 30 o in the circumferential direction. The stress computations confirm the previously made prediction that transverse stresses are within the acceptable concrete stress range.
Case 3 Finds the maximum transverse vertical stresses in symmetrical equivalent beam. Variations 3a
and 3b investigate anticipated conditions during normal operation. Vertical transverse forces
counteract each other under operating conditions largely through vertical compression forces in
the buttress concrete. See Figure 5D-10b. Variation 3c investigates anticipated conditions if
only one anchor force acts. See Figure 5D-10a.
3.3.1 REINFORCEMENT
FOR TRANSVERSE FORCES
The allowable concrete tensile stresses in the radial and vertical planes are not exceeded and
no special transverse reinforcement is required for design loads. However, reinforcement is
available in the form of buttress reinforcement in the hoop and meridional direction consisting of a radial tie-reinforcement anchored into the shell and spiral reinforcement in the maximum
bursting zone.
The buttress is reinforced sufficiently to resist all transverse tensile forces Zy (design) acting
within the limits of the buttress without relying on concrete tensile resistance.
Spiral: 5/8" ?, 3 1/2" pitch, d = 20", length = 36" Maximum tensile splitting stress, Table 5D-1, f by =0.178 ksi
T = 3.5 (0.178) (20) = 6.24 K , f s = 6.24 = 20.1 ksi 2 0.31
The spiral is designed for Zy (design) = 420 kips
Maximum transverse tensile force as computed is Z y = 311 kips.
TMI-1 UFSAR APPENDIX 5D 5D-7 REV. 18, APRIL 2006 Steel provided:
Spiral 22 - #5 = 6.82 in 2 Ties 10 - #10 = 12.70 19.52 in 2
f s = 420 = 21.5 ksi 19.52
Temperature changes in the concrete during the life of the plant or at accident condition will
produce a maximum of less than 1 ksi stress in the bursting zone reinforcement.
3.4 UNBALANCED
FORCES
Figure 5D-11 shows the tendon anchor forces for six buttresses in the Reactor Building.
The solid line horizontal force vector in Figure 5D-11 represents the tangential components of
the anchor forces. The anchor forces are applied on opposite faces of the buttresses and
counteract each other. However, since the vert ical spacing of the anchors varies, it is not possible to balance all the forces. The unbalanced horizontal forces must be transferred from
the buttress into the Reactor Building wall.
The two dashed force vectors, Figure 5D-11, adjacent and parallel to each anchor force vector
are obtained by proportioning the anchor force according to its position with respect to two
opposing anchors similar to determining the react ions of a simply supported beam with an unsymmetrical concentrated load. The dashed horizontal force vector on the buttress center
line indicates the unbalanced force between two opposing anchor forces. How these horizontal
forces are resisted is discussed later.
Since the horizontal anchor forces on opposing sides are offset from each other, the forces
must change direction in order to balance each other. Transverse forces in vertical direction are
a result of this angle change. The vertical dashed force vectors shown in Figure 5D-11
represent the unbalanced vertical forces behind each anchor zone.
Vertical forces in Figure 5D-11 were computed for the simplified assumption that the anchor
forces travel in straight lines from anchor plate to anchor plate after first moving 18 in.
horizontally into the concrete and gradually bend over into the inclined position. See Figure 5D-12.
The check analysis employing the Finite-Element Me thod did reveal a limited region of surface tensile stresses which presumable would produce cracking. However after concrete cracking
has occurred with the extent limited by s upplemental reinforcing a new state of equilibrium will exist. The following section discussing vertical and horizontal unbalanced forces verifies that a
mechanism does exist to maintain equilibrium.
3.4.1 VERTICAL
UNBALANCED FORCES
During initial prestressing and during tendon inspection the balance of opposing tendon anchors
might not exist and the total vertical transverse stress must be resisted. The maximum concrete
transverse TMI-1 UFSAR APPENDIX 5D 5D-8 REV. 18, APRIL 2006 tensile stress for this condition was computed in Table 5-2, Case 3c, f by = 115 psi. Figure 5D-10a shows the anticipated typical vertical stress-flow-lines for such conditions. The figure
also shows the "equivalent beam" for which transverse stresses were determined according to
Guyon's Symmetrical Prism Method. Sufficient reinforcing steel in form of a spiral has been
provided for Case 3c without relying on concrete tensile resistance. In addition, the vertical
reinforcement in the buttress, consisting of 7-#10, As = 8.8 in.,(Reference 2) can resist a
transverse splitting force of Z = 8.8 (0.95) (40) = 335 kips which is slightly more than the
expected force Zy = 311 kips shown in Table 5D-1. The reinforcement is placed mainly on the buttress surfaces in order to provide access for concrete placement. There is no need to
distribute this reinforcement throughout the cross-se ction of the buttress or to place it where the vertical transverse forces actually act. Wherever the reinforcement is placed, it will not prevent the concrete from cracking and will only tie t he concrete together preventing the development of uncontrolled cracks. Consequently, surface reinforcem ent is sufficient. It restrains the buttress in the vertical direction and no failure mechanism caused by vertical transverse stresses can
develop.
The highest vertical in-plane shear stresses at the buttress-cylinder intersection occur directly
beneath the ring-girder and at the cylinder-base s lab intersection haunch area. An amplification factor of 1.64 was applied to the shear stresses in both areas to account for the increased
stiffness at the buttress. Sufficient reinf orcing was provided between the buttress and cylinder to take the shear stresses.
3.4.2 HORIZONTAL
UNBALANCED FORCES
Unbalanced horizontal forces shown in Figure 5D-11 must be transferred into the wall.
During initial prestressing and tendon inspection the counteracting effect of opposing anchor
forces might be lost at any any anchor zone position. All anchor zones will therefore be reinforced as the analysis requires for such conditions. Figure 5D-9a shows the anticipated
typical horizontal stress-flow-lines for such conditions. The figure also shows the "equivalent beams" for which transverse stresses were determined according to Guyon's "Symmetrical Prism Method."
The analysis for horizontal transverse forces is tabulated in Table 5D-1. The critical loading
condition exists when a tendon is stressed without being counterbalanced by opposing anchor
forces. The maximum concrete transverse tens ile stress for this condition is represented by Case 2a, f by = 178 psi. Sufficient reinforcing steel in form of a spiral has been provided for Case 2a without relying on concrete tensile resistance. In addition, the radial tie reinforcement
in the buttress, consisting of 10-10#, A s = 12.7 in., (Reference 2) can resist a transverse splitting force of Z = 12.7 (0.95) (40) = 482 kips which is larger than the expected force Zy = 311 kips
shown in Table 5D-1. The tie reinforcement restrains the buttress in radial direction and no
failure mechanism caused by radial transverse stresses can develop.
The highest shear stresses for the case of no opposing prestress force in the buttress occurs at
the buttress-cylinder discontinuity on the side of the buttress opposite the existing prestress force. Reinforcing of #10 at 12 in. (See Figure 5D-4) are provided to take the resulting tensile
stresses.
TMI-1 UFSAR APPENDIX 5D 5D-9 REV. 18, APRIL 2006
3.4.3 EQUILIBRIUM
CONDITION IN THE ANCHOR BLOCK
The reinforcement is checked based on the assumption of an existing spalling crack in the
vertical direction. The forces and the extent of cracking in the anchor block were based on a
method developed by Gergely, Soyen, and Siess of arresting cracks due to transverse stresses (Reference 4). The idealized mathematical model used for this approach, see Figure 5D-13, is
quite conservative. The computation is shown in Table 5D-2. It should be noted that zero
moment exists at B' when c = 27 in., and that the governing moment, Mm, prevails at c = 42 in., that is, far beyond the face of the wall.
Concerning the visco-elastic behavior of concrete, it can be stated that the resulting
redistribution of stresses developed by Gergely will not significantly affect reinforcing requirements (Reference 5). The super-position of bursting stresses and normal stresses from vertical tendons results in a tri-axial state of stress which is less critical than during construction
stages. No additional investigations for the final state of stress in the anchorage zones, therefore, appear necessary. The various sour ces of errors in a numerical tri-axial analysis (uncertainties of material properties, element definition, etc.) may lead to less accurate values
than were obtained by two-dimensional methods.
The shear across the cylinder at the buttress-cylinder intersection for the case of unbalanced
prestressed force has been checked. Sufficient reinforcing had been provided to account for
the shear stresses. For accident temperature and pressure conditions a crack across the entire
cylinder at the buttress cylinder discontinuity has been assumed. Sufficient steel reinforcing is provided to take any tension where the compressive stresses due to prestress have been
overcome.
4 CHECK ANALYSIS BY FINITE-ELEMENT METHOD
4.1 INTRODUCTION
The finite-element method has been employed to check the design of tendon anchorage zones
in the Reactor Building. The purpose of this analysis was not as a design basis but simply as a
check to demonstrate the conservatism of the method described here before. This section
presents the method of analysis, load conditions, and results for the buttress and adjacent wall, and in conclusion to correlate these results with the design of the buttress reinforcement which
has been presented in Section 3.
The dominate loads are:
- a. Prestress forces
- b. Thermal loads
- c. Pressure forces due to loss of coolant accident (LOCA)
Although the stress analysis of the buttress and adjacent wall is a three dimensional problem, it can be reduced to a two dimensional problem, as the meridional stress of the cylindrical portion of the vessel is usually constant. The finite-element method of plane stress analysis (References 6 and 7) which can account for complex loadings, arbitrary geometric configurations, and arbitrary boundary conditions, assumes a homogeneous material (i.e., no TMI-1 UFSAR APPENDIX 5D 5D-10 REV. 18, APRIL 2006 cracks) and is an ideal tool for dealing with such a problem. Therefore, a plane stress finite-element computer program (Referenc e 7) was used in the buttress analysis.
For the plane stress finite-element program t he thermal loads were obtained from another computer program. The results from this progr am, for the steady state winter normal operation and the transient state during winter start-up, w inter shut-down, and winter loss of coolant accident were used as input data in the finite- element program.
As described later, a system of 768 direct st iffness simultaneous equations required solution.
To save computer storage, the Gauss-Seidel iter ative method (Reference 8), which requires the storage of only the non-zero elements of the sti ffness matrix, was employed. In order to expedite convergence, an over-relaxation factor (R eference 9) was also used in Reference 7; this factor usually ranges from 1.0 to 2.0.
Cracking of the concrete for the winter LOCA load condition is predicted. Areas of high tensile
stresses have been defined and reinforcement capable of resisting these provided.
4.2 PLANE
STRESS FINITE-ELEMENT COMPUTER PROGRAM
Due to the fact that the meridional stresses of the cylindrical portion of a containment vessel are
essentially uniform and constant, the horizontal cross section of the cylinder may be considered
as a section subjected to plane stresses only. The finite-element method is a versatile method
of analysis and is capable of providing states of stress in the buttress for complex loadings, arbitrary geometric configurations, and arbitrary boundary conditions. Tendon holes are not
included in the buttress geometry.
In the plane stress finite-element analysis, the horizontal cross sectional plane of the cylindrical portion of the vessel is replaced by a finite number of elements interconnected at a finite number of nodal points. Continuity between element s of the system is maintained by requiring that within each element, "lines initially straight remain straight in their displaced positions."
This requirement is satisfied if the strains within each element are assumed to be constant.
Therefore, the stresses within each element are also constant. Distributed loads are replaced
by equivalent nodal forces which act at the three nodes of the element. Based on these
assumptions, it is possible to derive the sti ffness matrix of a typical element, which is an expression for the nodal forces resulting from unit nodal displacements:
{p} = [k] {?} (1)
where {p} is the vector of six nodal forces which corresponds to the vector {?} of six nodal
displacements (Figure 5D-14), and (k) is the stiffness matrix of the element.
The nodal force vector {p} can be expressed as:
{p} = {PX i , Px j , Px k , Py i , Py j , Pyk}) (2) and the nodal displacements vector {?} can be expressed as,
{?} = {u i , u j , u k , v i , v j , v k} (3) which are shown in Figure 5D-14.
TMI-1 UFSAR APPENDIX 5D 5D-11 REV. 18, APRIL 2006 Formulation of Element Stiffness Matrix
Corresponding to the above assumptions, the displacements of a point at (x,y) are written as
linear functions of the coordinates:
u = u i + C 1 x + C 2 y v = v i + C 3 x + C 4 y (4)
The stiffness matrix [k] can be obtained by the following steps:
- a. Substitute the nodal coordinates into Equation(4) to obtain the displacement functions in terms of the nodal displacements.
- b. Substitute results of (a) into the linear strain-displacement equations to obtain the strains in terms of nodal displacements.
- c. Substitute the results of (b) into the stress-strain equation to obtain the stress in terms of nodal displacements.
- d. Find the nodal resultant forces which are statically equivalent to the stresses found in (c) so that the relationships between the nodal forces and the nodal
displacements can be written in a matrix form as shown in Equation (1).
The explicit formulations of the matrices which compose the element stiffness matrix can be found in Reference 7.
Solution of the Total Stiffness Equations After the stiffness matrices are formulated for each of the elements, they are assembled to form a total system stiffness matrix. If the boundary conditions and the loading conditions are properly specified, the tot al stiffness simultaneous equations are readily solvable to give the nodal displacements. Then, the nodal displacements are used to generate the stresses in each element and at each node. The directions and magnitudes of the principal
stresses can also be computed by using the Mohr's circle technique.
Iterative Method and Optimum Over-Relaxation Factor The Gauss-Seidel iterative procedure (Referenc e 8) is applied to solve the stiffness simultaneous equations. The method requires the storage of only the non-zero elements of the
stiffness matrix during the iterative process.
Since the stiffness matrix [k] is positive definite,the iterative method will always converge (Reference 10). The covergence rate is, however, often rather slow. An over-relaxation factor (Reference 9) can be applied which greatly incr eases the convergence rate. The determination of the proper over-relaxation factor is diffi cult and is based mostly on experience and trial and error.
4.3 GEOMETRY
OF BUTTRESSES WITH TENDONS
The pertinent geometry of the Reactor Building as it affects the buttress analysis is shown in Figure 5D-15. In Figure 5D-15a, a sectional view illustrates the vessel with six buttresses.
Figure 5D-15b, shows the arrangement of the hoop tendons on a stretch out of the cylindrical TMI-1 UFSAR APPENDIX 5D 5D-12 REV. 18, APRIL 2006 wall. Figure 5D-15c, shows the dimensions of one buttress and the arrangement of the hoop and vertical tendons.
4.4 THERMAL
LOADS AND TEMPERATURE PROFILES
If the temperature profile in a buttress is known, the difference between the temperature of any point and the reference temperature can be determined. The temperature difference, with the
modulus of elasticity and thermal expansion coeffi cient of the concrete, provides the equivalent constant stress for each finite-element. The constant thermal stress of each element can be
converted to equivalent nodal resultant forces which can be considered as external forces
applied at the nodal points of a non-heated element.
Temperature Profiles The temperature profiles and iso-thermal curves were calculated by a computer program developed by Brunnel (Reference 11) for the s olution of steady state and transient state problems in one, two, or three dimensional Cart esian or cylindrical systems. The governing differential equations were solved by the method of finite differences, using standard relaxation techniques.
The convergence criterion:
T j (step i) - T j (step i-1) <
0.00001 (5)
Y j (step i)
where T j is the temperature at node j, was applied at every node of the difference lattice. To ensure stability of the transient solutions, the time increment was selected according to the
method presented in Reference 11.
The buttress and wall temperatures were determined for the following conditions:
- a. Normal operation in winter
- b. Start-up in winter
- c. Shut-down in winter
- d. Loss of coolant accident during winter
The winter season was chosen because it produces the greatest thermal gradients in the shell
of the Reactor Building.
The mean monthly temperature was chosen due to t he fact that day- night fluctuations had little effect on the interior portion of the concrete. According to Carslow and Jeager, (Reference 12)
the amplitude of the temperature diminishes as:
wP con C p con exp [X
] (6) 2 k con Thus a distance of 12.8 inches from the outer wall of the Reactor Building the effect of the
oscillation is only 0.1 what it is on the outer face, and the middle and inner portion of the vessel
will be unaffected by day-night temperature oscillations.
TMI-1 UFSAR APPENDIX 5D 5D-13 REV. 18, APRIL 2006 For the start-up condition, it was assumed that the inside air temperature of the Reactor Building rose from 60 F to 110 F in one hour, and for shut-down the reverse was used. The
inside film coefficient was taken as 1.490 Btu/hr/sq ft/F for cases a, b, and c above. This was
modified film coefficient, which took into account the steel liner and the still air between the liner
and the concrete:
h m = 1 (7) 1 X steel X air h i +
+ k steel k air Equation 7 neglects the circumferential conduction and heat capacity in the steel liner and air
gap. However, due to the relatively small thi ckness of the liner and the air gap when compared with the thickness of the concrete wall, these effects may be neglected. The preliminary
analyses verified this assumption. The following values, obtained by theoretical analysis (Reference 13), were used in the thermal analysis:
inside film coefficient = 2.0 Btu/hr/sq ft/F
specific heat, C p con = 0.210 Btu/lb/F
thermal conductivity, k con = 10.0 Btu in./hr/sq ft/F
k steel = 300.0 Btu in./hr/sq ft/F
k air = 1.836 Btu in./hr/sq ft/F
The density of the concrete, P con was assumed to be 140 lbs/cu ft.
The above values agree closely with the results from experimental research carried out in Japan (Reference 14).
C p con = 0.200 Btu/lb/F
k con = 11.28 Btu in./hr/sq ft/F
The density of the concrete P con , was assumed to be 146 lb/cu ft.
The thickness of the steel and air gap was:
X air = 1/32 inch
X steel = 3/8 inch
Using Reference 15, the local film coefficient for air at 20 F and a free stream velocity of 20 mph
flowing past a concrete cylinder was calculated to be:
h = 0.56 Btu/hr/ft 2/F As a safe approximation the outside film coefficient was chosen as:
TMI-1 UFSAR APPENDIX 5D 5D-14 REV. 18, APRIL 2006 h? = 1.0 Btu/hr/sq ft/F
The buttress model was replaced by a mesh of 61 nodes, with adiabatic boundaries at the
centerline of the buttress and at the radial mid-section of the wall between two neighboring
buttresses (Figures 5D-16 and 5D-17). Due to symmetry, there is no heat transferred across
these boundaries. For Case d the liner temperature during LOCA was used as the inner
boundary temperature in conjunction with a film coefficient of 5.880 Btu/hr/sq ft/F, which
simulated the conduction through the air gap:
k air h m = (8) X air The transient state temperature analyses were performed at various time steps.
The temperature profiles for Cases a, b, c, and d are shown in Figure 5D-18. Because the most
severe temperature gradients occur for case a and d, only these two cases were considered in
the stress analyses. Iso-thermal curves fo r cases a and d are shown in Figures 5D-19 and 5D-20.
The design pressure inside the Reactor Building during LOCA is equal to 55 psig and the
maximum pressure was initially calculated to occur at approximately 200 seconds after initiation
of the accident. The maximum liner temperature of 280F was initially calculated to occur at
approximately 650 seconds after the initiation of the accident. However, the maximum
temperature gradient, which causes the most severe thermal stresses, was calculated to occur
at 10,000 seconds after initiation of the accident.
Consequently, the determination of the most severe stress state, due to thermal and pressure loads, requires that a number of stress analyses be conducted at various elapsed times after
commencement of the accident. A conservative estimate of the stresses was obtained, by
superimposing directly the maximum values for accident temperature and pressure.
Consequently although subsequent calculations alter t he time at which minimum pressure and liner temperature occur, these changes do not affect this analysis.
4.5 NON THERMAL LOADS
In addition to the thermal loads, there are other loads which act on the structure:
- a. Hoop Tendon Anchor Pressure
The tendon force was assumed to be distributed uniformly up the face of the buttress. The resulting average force (=p a) was applied to each of six nodes.
The location of the nodes approximately corresponded to the position of the
tendon.
- b. Hoop Tendon Pressure in a Radial Direction
Radial pressures, due to the hoop tendons in the wall, were calculated as average stresses in pounds/in.
The two radial pressures applied (refer to Figure 5D-21):
TMI-1 UFSAR APPENDIX 5D 5D-15 REV. 18, APRIL 2006
- 1. From point A to the tangent point T. Over this length the tendon spacing causing radial pressure (i.e., continuous past the buttress) is 33 in., giving
a pressure of p t , and
- 2. For the length from point T to point B (midway between buttresses) the average tendon spacing causing radial pressure is 16 1/2 in., giving a
pressure of 2p
- t.
- c. Internal Pressure Due to LOCA
Internal pressure acting on the inside face of the Reactor Building Wall was considered as a combination of two separate effects:
- 2. Pressure due to expansion of the liner caused by the increase in temperature during the LOCA condition (110F to 280F) based upon a limiting yield strength of 1.2 f
- y.
The above two loads were assumed to act simultaneously during the LOCA condition, giving an equivalent internal pressure of p
- i.
The magnitude and combination of the above non-thermal loads is noted in section 4.7
4.6 MATERIAL
PROPERTIES
Prestressing Steel Wire
The type of prestressing wire to be used in the construction is a stress relieved wire.
Dimensions, properties, etc., of the wire as used in the analysis were:
- a. Size of tendon = 169 wires
- b. Diameter of wires = 1/4 in.
- c. Minimum guaranteed ultimate tensile stress = 240,000 psi
- d. Tendon force @ 70 percent GUTS = 1400 kips
- e. Modulus of elasticity = 29 x 10 6 psi f. Coefficient of friction = 0.16
- g. Coefficient of wobble friction = 0.00030 TMI-1 UFSAR APPENDIX 5D 5D-16 REV. 18, APRIL 2006 Concrete The assumed properties of the concrete used in the analysis were:
- a. Compressive strength @ 28 days = 5000 psi
- b. Modulus of elasticity = 4 x 10 6 psi
- c. Poissons ratio = 0.20
Liner The liner is fabricated from 3/8 inch thick steel plate manufactured to ASTM A-283 Grade C.
The minimum yield stress for this material is 30,000 psi. In determining the equivalent liner pressure LOCA, the yield stress was assumed to be 36,000 psi.
4.7 LOADING
CASES AND COMBINATIONS
The buttress and wall were analyzed for four loading cases. The first two cases considered the
buttress and wall being subjected to the maximum prestress forces, i.e., prestress forces
allowing only those losses that would occur between prestressing and start-up (beginning of
plant life.) The latter two cases considered the buttress and wall being subjected to the
minimum prestress forces, i.e., prestress forc es after 40 year losses (end of plant life).
The load cases were:
- a. Winter normal operation with maximum prestress force.
- b. Winter LOCA with maximum prestress force.
- c. Winter normal operation with minimum prestress force.
- d. Winter LOCA with minimum prestress force.
These cases were considered to be more critical than the summer normal operation/LOCA
conditions. The summer condition results in a smaller temperature differential in the wall which
in turn results in maximum compressive and tensile stresses less than the winter condition.
Hence, the summer condition was not considered in the analyses.
TMI-1 UFSAR APPENDIX 5D 5D-17 REV. 18, APRIL 2006 The numerical value of the applied loads was:
LOAD CASE a b c d
t LT 110 150 110 150
t LP 110 280 110 280
t e 20 20 20 20
p a 6.25 6.73 5.52 6.03
pt A-T 46.07 49.75 39.25 42.2
2ptT-B 92.13 99.50 78.5 84.4
p i 0.0 96.1 0.0 96.1
where:
t LT = Temperature (F) of liner for thermal gradient in wall.
t LP = Temperature (F) of liner for equivalent pressure of liner due to expansion during LOCA.
t e = External air temperature (F)
p a = Prestress force at bearing plate (kips/inch)
pt A-T = Radial prestress force from point A to tangent point T (lbs/inch)
2p T-B = Radial prestress force tangent point T to point B (lbs/inch)
p i = Equivalent internal pressure (LOCA pressure + liner expansion pressure lbs/inch)
Figure 5D-21 shows the application of these loads on the buttress and wall.
Load p a p t , and 2p t for Case a incorporate losses due to elastic shortening of the concrete, friction and creep, shrinkage, and steel relaxation losses that occur between prestressing and
start-up. Case b prestress forces reflect increases from Case a, due to elastic elongation of the
structure during LOCA. Thus, Case a and b reflect the maximum prestress force that will exist
in the structure, at the earliest time in the life of the plant, that a LOCA could occur.
Loads p a , p t , and 2 pt for Cases c and d are similar to Cases a and be respectively, except that the losses due to creep, shrinkage, and steel relaxation are those that are assumed to occur
over 40 years to end of plant life. Case d pres tress forces reflect an increase from Case c due to elastic elongation of the structure during LOCA.
TMI-1 UFSAR APPENDIX 5D 5D-18 REV. 18, APRIL 2006
4.8 STRESS
ANALYSIS
Due to the symmetrical nature of the structur e only one-twelfth of the cylinder, as shown in Figure 5D-16, need be analyzed. Preliminary analy ses showed that beyond a certain point in the wall, the stress state, for all loading cases, was uniform. To conserve computer storage and
give a more accurate solution, it was decided to reduce the size of the model to that shown in
Figure 5D-22. Nodal displacements of both boundaries were specified. The displacements
were obtained form a computer program as developed by Kalnins (Reference 16) using a shell of revolution theory.
There were a total of 672 elements and 384 nodal points, so that 768 simultaneous equations
required solution. For the computations, the Gauss-Seidel iterative solution procedure was used.
4.9 ANALYSIS
OF RESULTS
Figures 5D-23 through 5D-30 show iso-stress curves for the maximum and minimum principal stresses for the four loading conditions outlined in Section 4.7. It should be noted that these
curves represent the algebraic maximum and minim um values of the principal stresses. The do not show the orientation of the stress trajectories.
Inspection of the iso-stress curves shows tha t, as expected, high compressive stresses exist directly beneath the bearing plate for all load conditions. Under winter normal conditions, a
region of tensile bursting stress occurs in the buttress at some distance removed from the
bearing plate; during LOCA, however, these bursting stresses merge with those resulting from the pressure and temperature loads.
The variation of temperature through both the wall and buttress leads to substantial hoop stress
gradients. It should be noted that in the curves for LOCA, the direction of the stresses in the
wall changes in the vicinity of the centerline.
Of primary interest, however, is the high localized t ensile stress that occurs at the buttress-wall re-entrant corner during LOCA, which results fr om both the thermal gradient as well as the local bending arising from the variation of resistance of the buttress and wall to internal pressure.
4.10 CONCLUSIONS
The results indicate that for the most severe loading conditions, the maximum stresses in the
buttress and adjacent wall are within safe limits, except for the localized area at the buttress
wall corner where tensile stresses greater than 6 ?f
? occur. Reinforcement to stabilize cracking in this local region are provided. All other results would indicate that the design method
employed was conservative.
TMI-1 UFSAR APPENDIX 5D 5D-19 REV. 18, APRIL 2006 5 REFERENCES
- 1. Building Code Requirements for Reinforced Concrete , ACI 318-63, American Concrete Institute, Detroit, Michigan, 1963.
- 2. Y. Guyon, Prestressed Concrete , Contractors Record Ltd. London and John Wiley and Sons Inc., New York, 1960.
- 3. F. Leonhardt, Prestressed Concrete Design and Construction , Wilhelm Ernest and Sohn, Berlin-Munich, 1964.
- 4. P. Gergely, M. A. Sozen and C. P. Siess, The Effect of Reinforcement on Anchorage Zone Cracks in Prestressed Concrete Members , Civil Engineering Studies, Structural Research Series No. 271, University of Illinois, Urbana, Illinois, July 1963.
- 5. Anchorage Systems in Prestressed Concrete Pressure Vessels; Anchorage Zone Problems , Oak Ridge National Laboratory, ORNL-TM-2378, January, 1969.
- 6. M. J. Turner, R. W. Clough, H. C. Martin, and L. J. Topp, "Stiffness and Deflection Analyses of Complex Structures," Jl. of Aero. Sci., V. 23, No. 9, Sept 1965, p. 805.
- 7. E. L. Wilson, Finite-Element Analysis of Two Dimensional Structures , Ph. D. Thesis, Univ. of Cal. at Berkeley, June 1963.
- 8. M. G. Salvadori and M. L. Baron, Numerical Methods in Engineering , Prentice-Hall, Englewood Cliffs, N. J., 1961.
- 9. F. G. Lehan, "Simultaneous Equations S olved by Over-Relaxization, "Proc. ASCE, 2nd Conf. on Electronic Computation , Pittsburgh, Pa., USA, Sept. 1960.
- 10. S. H. Crandall, Engineering Analysis , McGraw-Hill Book Co., N. Y., USA, 1956.
- 11. R. D. Brunnel, Computer Program Heating and Heating 2 , Atomic International Division, Report No. A1-64-MEMO-177, North American Aviation, Inc., Detroit, Michigan, USA, 1964.
- 12. H. S. Craslaw and J. C. Jaeger, Conduction of Heat in Solids , 2nd Ed., Oxford, London, 1959.
- 13. Johns-Manville Research and Engineering Center Report E-455-T-266 filed in Robert Emmett Ginna Final Facility Description and Safety Analysis Report, Section 5, Appendix 5B.
- 14. J. Kosaka, Experimental Research of Unsteady-State Heat Transfer Concerned with Nuclear Energy Facilities (Preliminary Report), Division of Nuclear Energy Facilities, Taisei Construction Ltd., 1969.
- 15. F. Kreith, Principles of Heat Transfer , 2nd Ed. International, Scranton, Pa., 1965.
- 16. A. Kalnins, Computer Programs , Crystal River Unit No., 3, Florida Power Corporation, PSAR, Pages 5C-1 and 5C-2.
- 17. J. Zielinski and R. E. Rowe, An Investigation of the Stress Distribution in the Anchorage Zones of Post-Tensioned Concrete Members , Cement and Concrete Association, Research Report NO. 9, London, September 1960.
TMI-1 UFSAR APPENDIX 5D 5D-20 REV. 18, APRIL 2006
- 18. J. Zielinski and R. E. Rowe, The Stress Distribution Associated with Groups of Anchorages in Post-Tensioned Concrete Members , Cement and Concrete Association, Research Report No. 13, London, October 1962.
- 19. P. Gergely, Anchorage Systems in Prestressed Concrete Pressure Vessels; Anchorage Zone Problems , U.S. Atomic Energy Commission, ORNL-TM-2378, January, 1969.
TMI-1 UFSAR APPENDIX 5D 5D-21 REV. 18, APRIL 2006 TABLE 5D-1 BUTTRESS ANCHOR ZONE TRANSVERSE STRESSES IN RADIAL AND VERTICAL PLANES Case Radial Plane
1a 1b 1c 2a
2a'(in) 20.5 20.5 20.5 20.5
2a (in) 52 26 78 32
2a'/2a 0.40 0.79 0.26 0.64
fy[1] 0.27 0.10 0.34 0.15
z (in) [2] 0.80 0.90 0.70 0.9
2b' (in) 20.5 20.5 20.5 20.5
2b=2b'+(2az)tg30 o 44.5 33.5 51.5 37.0 A=2a2b (in
- 2) 2315 870 4150 1182 P (kips) [3] 1400 1400 1400 1400
f bz = P/A (ksi) 0.605 1.610 0.338 1.182 f by = f bz fy (ksi) 0.164 0.161 0.115 0.178 [4]
Zy = 0.3P (1-a'/a) 252 88 311 151 (kips) [5]
Zy (design) = 0.3P 420 420 420 420 (kips) [6]
az (in) 20.8 11.7 27.3 14.4
[1] Guyon, "Prestressed Concrete", Figure 9.
2 [2] Guyon, "Prestressed Concrete", Figure 9.
2 [3] 1600 Kips and 2000 kips anchor forces is not critical.
[4] Changes in tensile splitting stresses with varying "2a" dimension.
[5] Zy = total transverse force Leonhardt "Prestressed Concrete"
- p. 271.
3 [6] For design purpose "(a'/a)" is conservatively assumed to be equal
to zero.
TMI-1 UFSAR APPENDIX 5D 5D-22 REV. 18, APRIL 2006 TABLE 5D-1 (continued)
BUTTRESS ANCHOR ZONE TRANSVERSE STRESSES IN RADIAL AND VERTICAL PLANES
Case Ultimate Load- Vertical Planes Radial Plane 2a 3a 3b 3c
2a' (in) 20.5 20.5 20.5 20.5
2a (in) 32 33 26 78
2a'/2a 0.64
fy [1] 0.15
z (in) [2] 0.9
2b' (in) 20.5
2b=2b' + (2az)tg30 0 37.0 See Case 2a, 1b, and 1c A=2a2b (in
- 2) 1182 P (kips) [3] 2000
f bz = P/A (ksi) 1.690 f by = f bz fy (ksi) 0.254 [4]
Zy=0.3P (1-a'/a) 216 (kips) [5]
Zy (design)=0.3P 600 (kips) [6]
az (in) 14.4
[1] Guyon, "Prestressed Concrete", Figure 9.
2 [2] Guyon, "Prestressed Concrete", Figure 9.
2 [3] 1600 Kips and 2000 kips anchor forces is not critical.
[4] Changes in tensile splitting stresses with varying "2a" dimension.
[5] Zy = total transverse force Leonhardt "Prestressed Concrete" p. 271.
3 [6] For design purpose "(a'/a)" is conservatively assumed to be equal to zero.
TMI-1 UFSAR APPENDIX 5D 5D-23 REV. 18, APRIL 2006 TABLE 5D-2 BUTTRESS ANCHOR ZONE SPALLING (Refer to Figure 5D-13)
M = P {c - e - (c/h) 2 [2h - 3e - c + 2ec
]} h c= varies e = 15.5 in. h = 69 in.
M = P {c - 15.5 - (c/69) 2 [2 x 69 - 3 x 15.5 - c ( 1 - 2x15.5
)]} 69
C (in.) c - 15.5 (c/69) 2 -c( 1 - 2 x 15.5) M/P M m (in. - K) 69
5 -10.5 0.00525 - 2.75 -11.0 10 - 5.5 0.02100 - 5.50 - 7.3 15.5 0 0.05040 - 8.50 - 4.4 20 + 4.5 0.08400 -11.00 - 2.3 25 + 9.5 0.1310 -13.80 -0.7 zero at 27 in. 30 +14.5 0.1890 -16.50 + 0.3 42 +26.5 0.370 -23.10 + 1.2 1680 50 +34.5 0.525 -27.50 + 0.9 (in. - kips) 69 +53.5 1.00 -38.00 0
F t = M m Z = 3 in h = 69 in.
h-z
f s allowable = 10 5 W A s A s = 1.27 sq in. (#10 bars)
W = 0.005
f s = 10 5 0.005/1.27 = 6.28 ksi F T = 1680
= 25.5 K 69 - 3
A s = F T/f s = 25.5/6.28 = 4.05 sq in.
A s provided: #10 @ 12 in. + 1 #10 Loop = 6.03 sq in.
TMI-1 UFSAR APPENDIX 5D 5D-i REV. 18, APRIL 2006
APPENDIX 5D REPORT ON STRESS ANALYSIS OF BUTTRESS
TMI-1 UFSAR APPENDIX 5D 5D-ii REV. 18, APRIL 2006 TABLE OF CONTENTS
SECTION TITLE 1 DESIGN DETAILS
1.1 TENDON
ANCHORAGE HARDWARE
1.2 ANCHORAGE
ZONE REINFORCEMENT
2 TENDON LOADS AND ALLOWABLE STRESSES
2.1 TENDON
LOADS
2.2 ALLOWABLE
STRESSES
3 ANCHOR BLOCK STRESSES
3.1 BEARING
STRESSES
3.2 EQUIVALENT
FLUID PRESSURE METHOD
3.3 TRANSVERSE
TENSION-PRISM METHOD
3.3.1 REINFORCEMENT
FOR TRANSVERSE FORCES
3.4 UNBALANCED
FORCES
3.4.1 VERTICAL
UNBALANCED FORCES
3.4.2 HORIZONTAL
UNBALANCED FORCES
3.4.3 EQUILIBRIUM
CONDITION IN THE ANCHOR BLOCK
4 CHECK ANALYSIS BY FINITE-ELEMENT METHOD
4.1 INTRODUCTION
4.2 PLANE
STRESS FINITE-ELEMENT COMPUTER PROGRAM
4.3 GEOMETRY
OF BUTTRESSES WITH TENDONS
4.4 THERMAL
LOADS AND TEMPERATURE PROFILES 4.5 NON THERMAL LOADS
4.6 MATERIAL
PROPERTIES
4.7 LOADING
CASES AND COMBINATIONS
4.8 STRESS
ANALYSIS
4.9 ANALYSIS
OF RESULTS
4.10 CONCLUSIONS
5 REFERENCES
TMI-1 UFSAR APPENDIX 5D 5D-iii REV. 18, APRIL 2006 LIST OF TABLES
TABLE TITLE 5D-1 Buttress Anchor Zone Transverse Stresses in Radial and Vertical Planes
5D-2 Buttress Anchor Zone Spalling
TMI-1 UFSAR APPENDIX 5D 5D-iv REV. 18, APRIL 2006 LIST OF FIGURES
FIGURE TITLE 5D-1 DELETED
5D-2 DELETED
5D-3 DELETED
5D-4 DELETED
5D-5 DELETED
5D-6 DELETED
5D-7 Typical Hoop Tendon Arrangement
5D-8 Buttress-Hoop Tendon Bearing Stresses
5D-9 Buttress-Radial Bursting Stresses
5D-10 Buttress-Vertical Bursting Stresses
5D-11 Hoop Tendon Typical Anchor Force Distribution
5D-12 Hoop Prestress Idealized Vertical Stress Flow in Buttress
5D-13 Idealized Model for Spalling
5D-14 A Triangular Plane Stress Finite-Element
5D-15 Reactor Building With Buttresses
5D-16 Adiabatic Sections
5D-17 Model of Gridwork for Thermal Analysis
5D-18 Temperature Profiles Through Wall and Buttress for Various Thermal Load Conditions
5D-19 Iso-Thermal Curves for Winter Normal Operating Temperature Condition
5D-20 Iso-Thermal Curves for Winter A ccident Temperature Condition at Time = 10,000 Secs.
5D-21 The Hoop Tendon Loads, Internal Pressure and the Equivalent Liner Pressure
5D-22 Finite-Element Grid Work for Buttress Analysis TMI-1 UFSAR LIST OF FIGURES
FIGURE TITLE APPENDIX 5D 5D-v REV. 18, APRIL 2006 5D-23 Minimum Principal Stresses for Maximum Prestress Forces and Winter Normal Operating Condition (Case A) (PSI)
5D-24 Maximum Principal Stresses for Maximum Prestress Forces and Winter Normal Operating Condition (Case A) (PSI)
5D-25 Minimum Principal Stresses for Maximum Prestress Forces and Winter Accident Condition (Case B) (PSI)
5D-26 Maximum Principal Stresses for Maximum Prestress Forces and Winter Accident Condition (Case B) (PSI)
5D-27 Minimum Principal Stresses for Minim um Prestress Forces and Winter Normal Operating Condition (Case C) (PSI)
5D-28 Maximum Principal Stress for Minimum Prestress Condition (Case C) (PSI)Forces and Winter Normal Operating
5D-29 Minimum Principal Stresses for Minim um Prestress Condition (Case D) (PSI)
Forces and Winter Accident
5D-30 Maximum Principal Stresses for Minimum Prestress Forces and Winter Accident Condition (Case D) (PSI)
TMI-1 UFSAR APPENDIX 5D 5D-1 REV. 18, APRIL 2006 1 DESIGN DETAILS
1.1 TENDON
ANCHORAGE HARDWARE
The prestressing system employs the BBRV anchor age system utilizing 169-1/4" diameter wires conforming to ASTM A421-65, Type BA. The dimensions and capacities of the tendons used
hereafter for the anchorage zone analysis are more fully described in Appendix 5B, Part D.
1.2 ANCHORAGE
ZONE REINFORCEMENT
The mild steel reinforcement in the tendon anchorage zones of the buttress and the base mat
are shown on GAI Drawings 421004, 421005, 421006, 421031,421032, and 421033. This
reinforcement was conservatively designed by the methods enumerated hereafter in Section 3 of this Appendix and is judged to be sufficient to accommodate all predicted stress states
without jeopardizing the ability to satisfactorily place the concrete.
2 TENDON LOADS AND ALLOWABLE STRESSES
2.1 TENDON
LOADS
The minimum guaranteed ultimate tensile strength of the tendons is approximately 2,000 kips.
Before the tendons are anchored (locked off) at 70 percent of the ultimate tendon force, which is
equal to 1,400 kips, they may be temporarily ov erstressed to 1,600 kips in order to reduce friction losses. The anchor zones are designed for anchor forces of 1,400 kips per tendon. The
anchor zones are also checked to resist safely temporary anchor force of 1,600 kips per tendon.
2.2 ALLOWABLE
STRESSES
The allowable stresses for concrete and steel will conform to those given in Section 5.2.2 and
5.2.3 of the FSAR.
Due to the specialized nature of anchor zone design the following conditions and stress criteria
will be used in the design. These stress criteria include adequate conservatism in their
application.
First, verify that the allowable concrete compression and tensile stresses are not exceeded due
to bearing and bursting stresses respectively.
Second, check the tensile stresses in reinforcing steel when no credit is taken for concrete
tensile capacity in the investigated anchor zones.
Under these conditions, the allowable stresses will be shown as follows:
TMI-1 UFSAR APPENDIX 5D 5D-2 REV. 18, APRIL 2006 CONCRETE
- 1. Allowable bearing stress
according to ACI/318-63, 1 Paragraph 2605(C)
- 2. Allowable compression stress in the structure beyond the anchorage zone c c f 0.45 f=
- 3. Allowable Tension Stress
REINFORCING STEEL
Allowable tension for
A615-Grade 40, f s = 20 ksi.
3 ANCHOR BLOCK STRESSES The investigation of the buttress anchorage zones for the hoop tendons is described in the following sections. The analytical work as presented, is based on the end anchor zone stresses
and stress flow pattern studies documented by Y. Guyon (Reference 2) and F. Leonhardt (Reference 3). Their work has been substant iated through extensive strain measurements and stress flow pattern established by studies utilizing photoelastic techniques.
The design presented herein is quite conservative; nevertheless, it has been adopted as the
basis for determination of required reinforcing steel in the end anchor zones to overcome tensile
splitting stresses and to provide for spalling stresses near the concrete surface.
The existance of a three-dimensional stress field including high shear, thermal and shrinkage
tensile stresses is influenced by hairline cracking due to creep and thermal stresses. In order to
accommodate the most unfavorable stress combination in the anchorage zones, the following
design goals were established:
- 1. Maximum concrete bursting tensile stresses have been kept below maximum allowable tensile stress of conservatively c f 3using Y. Guyon's symmetrical prism method. This would indicate that no special spiral and bursting reinforcement is required. The 3 b b ci cp/A A' f 6.0 f=ksi. 5.0 f ci=ksi. 5.0 f c=c f 3ksi. 5.0 f c=
TMI-1 UFSAR APPENDIX 5D 5D-3 REV. 18, APRIL 2006 combined stresses during normal operation and test loading also indicate that the concrete will remain uncracked in the anchorage zone.
- 2. However, due to the existance of complex thermal stresses and possible microcracking, the ability of concrete to carry tensile stresses has been discarded, and mild steel
reinforcing has been provided for:
- a. the equivalent fluid pressure
- b. transverse tensile stresses
- c. unbalanced vertical forces
- d. unbalanced horizontal forces
- e. equilibrium condition
The combined stresses during the design pressure for the Reactor Building (p = 55 psig) indicate that the concrete will remain uncracked in the anchorage zone. The combined
stresses during the factored accident loading (refer to 5.2.3.2, cases a. and b. indicate
that the concrete in the anchorage zone may develop hairline cracks.
- 3. Existing theoretical solutions for anchorage zone design do not agree with more recent experimental results (References 17 and 18); maximum transverse tensile forces were
made up to 5.8 times larger than determined by Y. Guyon.
- 4. Possible failure mechanisms due to spalling and bursting stresses were evaluated to verify static equilibrium and acceptable stre ss limits under such hypothetical conditions.
- 5. Bursting stresses in a large body are much smaller than in a member bounded by planes (Reference 19). Therefore, while it was found that two-dimensional theories underestimate
the bursting stresses (Reference 17), analytical results based on isolated prisms or better, experimental values such as those developed by J. Zielinski and R. E. Rowe can be used for massive structures (Reference 19).
- 6. The design of the anchorage zones provides a margin of safety similar to the entire reactor building. Furthermore, the reinforcement is designed to control cracking due to
spalling, bursting, or shear of the concrete under sustained load.
3.1 BEARING
STRESSES
Figures 5D-7 through 5D-10 show the typical hoop tendon arrangement and buttress anchor
zone.
Bearing area:
2 2 2 b in. 481 4 11-24 A== Concrete anchor surface:
2 2 2 b in. 994 4 11-33 A= Allowable bearing stress:
()ksi 3.82 1/3 481 994 5.0 0.6 f cp== Use of 28 day concrete strength is conservative since initial tensioning will occur significantly
later in time.
TMI-1 UFSAR APPENDIX 5D 5D-4 REV. 18, APRIL 2006 Bearing pressure for anchor force:
- a. ksi 3.40 481 1633.5 f kips, 1633.5 P b===<ksi 3.82 f cp= b. ksi 3.88 481 1866.8 f kips, 1866.8 P b=== The depth below the bearing plate, at which concrete stresses are identical to those allowed beyond the anchorage zone area, is computed on the conservative assumption that the load
spreads at a 30 o angle from the plane normal to the bearing plate surface.
Obtain the required concrete area to satisfy the above condition for:
P = 1633.5 kips, A g = 1633.5 = 727 in 2 2.25
If it is conservatively assumed that the anchor force spreads at an angle of 30 o , the depth (h) below the bearing plate at which the concrete stress will have been reduced to the allowable
compression stress as shown above, can now be computed. See Figure 5D-8b. +=4 7-)30 2htan (20.5 A 2 2 g in. 4.5 h 38.5, -3 2h 20.5 622 2=+= Special provisions must be made to confine the concrete within the 4.5 in. below the bearing plate. From then on, spiral reinforcement m eeting the requirement for tensile splitting stresses will be investigated as described in part 3.3.
3.2 EQUIVALENT
FLUID PRESSURE METHOD
Directly below the bearing plate the necessary reinforcement for confining the concrete is
computed for an assumed equivalent fluid pressure "p". The equivalent fluid pressure is equal
to the concrete compression stress multiplied by the Poissons Ratio, which is assumed to be
0.2.
The concrete compression stress 2 in. below the bearing plate,which is about at the center of
the critical 4.5 in. depth, will be used to determine the equivalent fluid pressure. The initial
1,400 kips anchor force is considered. See Figure 5D-8b.
[]2 2 2 g in 481.5))4 ((7- )30 (tan 4 5.20 (A=+= ksi 2.91 481.5 1400 f c== The equivalent fluid pressure will be:
p = 2.91 (0.2) = 0.58 ksi TMI-1 UFSAR APPENDIX 5D 5D-5 REV. 18, APRIL 2006 The reinforcement is arranged as shown on GAI Drawing 421031 . The #5 spiral has a
diameter of 20 inches with a 3 1/2 in. pitch, and a total length of 36 inches. Horizontal ties
2-#10 and vertical reinforcement 1-#10 plus 1-#8 are also available to resist this equivalent fluid pressure.
The required steel area for the 20 inch diameter spiral was computed for an allowable steel
stress of f s = 20 ksi.
Required steel area:
A s = 0.58(4.5)(20)
= 2.61 in 2 (steel area req'd. over a 4.5 in length below the bearing plate) 20
Steel provided:
- a. Radially, Spiral 2-#5 = 0.62 in 2 Ties 2-#10 = 2.54 3.16 in 2 > 2.61 in 2
- b. Vertically, Spiral 2-#5 = 0.62 Vert. 1-#10 = 1.27 Vert. 1-#8 0.79 2.68 in 2 > 2.61 in 2
3.3 TRANSVERSE
TENSION-PRISM METHOD
Transverse stresses are produced when the anchor force spreads out into the surrounding
concrete. Directly behind the bearing plate the stress lines are parallel to the tendon axis. The
angle change of the stress-flow-lines, where they bend outwards, produces in these regions
compression in the concrete which is of no concern. Further away, the stress-flow-lines tend to
bend back and tensile transverse stresses are created as a result of this angle change. Figure
5D-9 and 5D-10 show the idealized compression stress-flow-lines in plane and the approximate location of transverse stresses.
Further, it can be seen that the angle changes will be small, and that the stress intensity along
the compression stress-flow-lines reduces as they spread out into the concrete mass. From
these considerations it can be concluded that the resulting tensile transverse stresses will be
small and within the tensile capacity of the concrete. Deformed vertical, circumferential, and
radial reinforcement is provided throughout the buttr ess in addition to the spiral reinforcement at each bearing plate. The reinforcement will tie the concrete mass together and ensures that
thermal cracking of the concrete mass, if it should occur, is of no significance as far as structural
integrity of the Reactor Building is concerned.
Figure 5D-10 shows that the critical condition, with respect to transverse tensile stresses in the
vertical plane, exists when only individual tendons ar e stressed. When all tendons are stressed, the stress-flow-lines are more or less parallel. Figure 5D-9 and 5D-10 show, superimposed
over the anticipated stress-flow-lines, "equivalent beams" for which the transverse stresses can be computed. Computation for transverse stresses in radial planes are performed in tabulated
form in Table 5D-1. Three cases are analyzed as follows:
TMI-1 UFSAR APPENDIX 5D 5D-6 REV. 18, APRIL 2006 Case 1 Finds the maximum transverse radial stress in symmetrical equivalent beam. Variation la
investigates anticipated conditions during normal operation. Variations 1b and 1c are made for
extremely shallow and extremely deep "equivalent beams" to establish how variations of beam depths affect the magnitude and location of the transverse stresses.
Case 2 Investigates the effect of eccentric loads on equivalent beams. Guyon's "Symmetrical Prism
Method" is used. According to this method the depth of the eccentric beam does not affect the
magnitude of radial transverse stresses; only the dimension of the prism drawn symmetrical to the tendon center line determines their magnitude. Possible variations in eccentric beam
depths are therefore of no concern.
All stresses were computed on the conservative assumption that the load spreads at 30 o in the circumferential direction. The stress computations confirm the previously made prediction that transverse stresses are within the acceptable concrete stress range.
Case 3 Finds the maximum transverse vertical stresses in symmetrical equivalent beam. Variations 3a
and 3b investigate anticipated conditions during normal operation. Vertical transverse forces
counteract each other under operating conditions largely through vertical compression forces in
the buttress concrete. See Figure 5D-10b. Variation 3c investigates anticipated conditions if
only one anchor force acts. See Figure 5D-10a.
3.3.1 REINFORCEMENT
FOR TRANSVERSE FORCES
The allowable concrete tensile stresses in the radial and vertical planes are not exceeded and
no special transverse reinforcement is required for design loads. However, reinforcement is
available in the form of buttress reinforcement in the hoop and meridional direction consisting of a radial tie-reinforcement anchored into the shell and spiral reinforcement in the maximum
bursting zone.
The buttress is reinforced sufficiently to resist all transverse tensile forces Zy (design) acting
within the limits of the buttress without relying on concrete tensile resistance.
Spiral: 5/8" ?, 3 1/2" pitch, d = 20", length = 36" Maximum tensile splitting stress, Table 5D-1, f by =0.178 ksi
T = 3.5 (0.178) (20) = 6.24 K , f s = 6.24 = 20.1 ksi 2 0.31
The spiral is designed for Zy (design) = 420 kips
Maximum transverse tensile force as computed is Z y = 311 kips.
TMI-1 UFSAR APPENDIX 5D 5D-7 REV. 18, APRIL 2006 Steel provided:
Spiral 22 - #5 = 6.82 in 2 Ties 10 - #10 = 12.70 19.52 in 2
f s = 420 = 21.5 ksi 19.52
Temperature changes in the concrete during the life of the plant or at accident condition will
produce a maximum of less than 1 ksi stress in the bursting zone reinforcement.
3.4 UNBALANCED
FORCES
Figure 5D-11 shows the tendon anchor forces for six buttresses in the Reactor Building.
The solid line horizontal force vector in Figure 5D-11 represents the tangential components of
the anchor forces. The anchor forces are applied on opposite faces of the buttresses and
counteract each other. However, since the vert ical spacing of the anchors varies, it is not possible to balance all the forces. The unbalanced horizontal forces must be transferred from
the buttress into the Reactor Building wall.
The two dashed force vectors, Figure 5D-11, adjacent and parallel to each anchor force vector
are obtained by proportioning the anchor force according to its position with respect to two
opposing anchors similar to determining the react ions of a simply supported beam with an unsymmetrical concentrated load. The dashed horizontal force vector on the buttress center
line indicates the unbalanced force between two opposing anchor forces. How these horizontal
forces are resisted is discussed later.
Since the horizontal anchor forces on opposing sides are offset from each other, the forces
must change direction in order to balance each other. Transverse forces in vertical direction are
a result of this angle change. The vertical dashed force vectors shown in Figure 5D-11
represent the unbalanced vertical forces behind each anchor zone.
Vertical forces in Figure 5D-11 were computed for the simplified assumption that the anchor
forces travel in straight lines from anchor plate to anchor plate after first moving 18 in.
horizontally into the concrete and gradually bend over into the inclined position. See Figure 5D-12.
The check analysis employing the Finite-Element Me thod did reveal a limited region of surface tensile stresses which presumable would produce cracking. However after concrete cracking
has occurred with the extent limited by s upplemental reinforcing a new state of equilibrium will exist. The following section discussing vertical and horizontal unbalanced forces verifies that a
mechanism does exist to maintain equilibrium.
3.4.1 VERTICAL
UNBALANCED FORCES
During initial prestressing and during tendon inspection the balance of opposing tendon anchors
might not exist and the total vertical transverse stress must be resisted. The maximum concrete
transverse TMI-1 UFSAR APPENDIX 5D 5D-8 REV. 18, APRIL 2006 tensile stress for this condition was computed in Table 5-2, Case 3c, f by = 115 psi. Figure 5D-10a shows the anticipated typical vertical stress-flow-lines for such conditions. The figure
also shows the "equivalent beam" for which transverse stresses were determined according to
Guyon's Symmetrical Prism Method. Sufficient reinforcing steel in form of a spiral has been
provided for Case 3c without relying on concrete tensile resistance. In addition, the vertical
reinforcement in the buttress, consisting of 7-#10, As = 8.8 in.,(Reference 2) can resist a
transverse splitting force of Z = 8.8 (0.95) (40) = 335 kips which is slightly more than the
expected force Zy = 311 kips shown in Table 5D-1. The reinforcement is placed mainly on the buttress surfaces in order to provide access for concrete placement. There is no need to
distribute this reinforcement throughout the cross-se ction of the buttress or to place it where the vertical transverse forces actually act. Wherever the reinforcement is placed, it will not prevent the concrete from cracking and will only tie t he concrete together preventing the development of uncontrolled cracks. Consequently, surface reinforcem ent is sufficient. It restrains the buttress in the vertical direction and no failure mechanism caused by vertical transverse stresses can
develop.
The highest vertical in-plane shear stresses at the buttress-cylinder intersection occur directly
beneath the ring-girder and at the cylinder-base s lab intersection haunch area. An amplification factor of 1.64 was applied to the shear stresses in both areas to account for the increased
stiffness at the buttress. Sufficient reinf orcing was provided between the buttress and cylinder to take the shear stresses.
3.4.2 HORIZONTAL
UNBALANCED FORCES
Unbalanced horizontal forces shown in Figure 5D-11 must be transferred into the wall.
During initial prestressing and tendon inspection the counteracting effect of opposing anchor
forces might be lost at any any anchor zone position. All anchor zones will therefore be reinforced as the analysis requires for such conditions. Figure 5D-9a shows the anticipated
typical horizontal stress-flow-lines for such conditions. The figure also shows the "equivalent beams" for which transverse stresses were determined according to Guyon's "Symmetrical Prism Method."
The analysis for horizontal transverse forces is tabulated in Table 5D-1. The critical loading
condition exists when a tendon is stressed without being counterbalanced by opposing anchor
forces. The maximum concrete transverse tens ile stress for this condition is represented by Case 2a, f by = 178 psi. Sufficient reinforcing steel in form of a spiral has been provided for Case 2a without relying on concrete tensile resistance. In addition, the radial tie reinforcement
in the buttress, consisting of 10-10#, A s = 12.7 in., (Reference 2) can resist a transverse splitting force of Z = 12.7 (0.95) (40) = 482 kips which is larger than the expected force Zy = 311 kips
shown in Table 5D-1. The tie reinforcement restrains the buttress in radial direction and no
failure mechanism caused by radial transverse stresses can develop.
The highest shear stresses for the case of no opposing prestress force in the buttress occurs at
the buttress-cylinder discontinuity on the side of the buttress opposite the existing prestress force. Reinforcing of #10 at 12 in. (See Figure 5D-4) are provided to take the resulting tensile
stresses.
TMI-1 UFSAR APPENDIX 5D 5D-9 REV. 18, APRIL 2006
3.4.3 EQUILIBRIUM
CONDITION IN THE ANCHOR BLOCK
The reinforcement is checked based on the assumption of an existing spalling crack in the
vertical direction. The forces and the extent of cracking in the anchor block were based on a
method developed by Gergely, Soyen, and Siess of arresting cracks due to transverse stresses (Reference 4). The idealized mathematical model used for this approach, see Figure 5D-13, is
quite conservative. The computation is shown in Table 5D-2. It should be noted that zero
moment exists at B' when c = 27 in., and that the governing moment, Mm, prevails at c = 42 in., that is, far beyond the face of the wall.
Concerning the visco-elastic behavior of concrete, it can be stated that the resulting
redistribution of stresses developed by Gergely will not significantly affect reinforcing requirements (Reference 5). The super-position of bursting stresses and normal stresses from vertical tendons results in a tri-axial state of stress which is less critical than during construction
stages. No additional investigations for the final state of stress in the anchorage zones, therefore, appear necessary. The various sour ces of errors in a numerical tri-axial analysis (uncertainties of material properties, element definition, etc.) may lead to less accurate values
than were obtained by two-dimensional methods.
The shear across the cylinder at the buttress-cylinder intersection for the case of unbalanced
prestressed force has been checked. Sufficient reinforcing had been provided to account for
the shear stresses. For accident temperature and pressure conditions a crack across the entire
cylinder at the buttress cylinder discontinuity has been assumed. Sufficient steel reinforcing is provided to take any tension where the compressive stresses due to prestress have been
overcome.
4 CHECK ANALYSIS BY FINITE-ELEMENT METHOD
4.1 INTRODUCTION
The finite-element method has been employed to check the design of tendon anchorage zones
in the Reactor Building. The purpose of this analysis was not as a design basis but simply as a
check to demonstrate the conservatism of the method described here before. This section
presents the method of analysis, load conditions, and results for the buttress and adjacent wall, and in conclusion to correlate these results with the design of the buttress reinforcement which
has been presented in Section 3.
The dominate loads are:
- a. Prestress forces
- b. Thermal loads
- c. Pressure forces due to loss of coolant accident (LOCA)
Although the stress analysis of the buttress and adjacent wall is a three dimensional problem, it can be reduced to a two dimensional problem, as the meridional stress of the cylindrical portion of the vessel is usually constant. The finite-element method of plane stress analysis (References 6 and 7) which can account for complex loadings, arbitrary geometric configurations, and arbitrary boundary conditions, assumes a homogeneous material (i.e., no TMI-1 UFSAR APPENDIX 5D 5D-10 REV. 18, APRIL 2006 cracks) and is an ideal tool for dealing with such a problem. Therefore, a plane stress finite-element computer program (Referenc e 7) was used in the buttress analysis.
For the plane stress finite-element program t he thermal loads were obtained from another computer program. The results from this progr am, for the steady state winter normal operation and the transient state during winter start-up, w inter shut-down, and winter loss of coolant accident were used as input data in the finite- element program.
As described later, a system of 768 direct st iffness simultaneous equations required solution.
To save computer storage, the Gauss-Seidel iter ative method (Reference 8), which requires the storage of only the non-zero elements of the sti ffness matrix, was employed. In order to expedite convergence, an over-relaxation factor (R eference 9) was also used in Reference 7; this factor usually ranges from 1.0 to 2.0.
Cracking of the concrete for the winter LOCA load condition is predicted. Areas of high tensile
stresses have been defined and reinforcement capable of resisting these provided.
4.2 PLANE
STRESS FINITE-ELEMENT COMPUTER PROGRAM
Due to the fact that the meridional stresses of the cylindrical portion of a containment vessel are
essentially uniform and constant, the horizontal cross section of the cylinder may be considered
as a section subjected to plane stresses only. The finite-element method is a versatile method
of analysis and is capable of providing states of stress in the buttress for complex loadings, arbitrary geometric configurations, and arbitrary boundary conditions. Tendon holes are not
included in the buttress geometry.
In the plane stress finite-element analysis, the horizontal cross sectional plane of the cylindrical portion of the vessel is replaced by a finite number of elements interconnected at a finite number of nodal points. Continuity between element s of the system is maintained by requiring that within each element, "lines initially straight remain straight in their displaced positions."
This requirement is satisfied if the strains within each element are assumed to be constant.
Therefore, the stresses within each element are also constant. Distributed loads are replaced
by equivalent nodal forces which act at the three nodes of the element. Based on these
assumptions, it is possible to derive the sti ffness matrix of a typical element, which is an expression for the nodal forces resulting from unit nodal displacements:
{p} = [k] {?} (1)
where {p} is the vector of six nodal forces which corresponds to the vector {?} of six nodal
displacements (Figure 5D-14), and (k) is the stiffness matrix of the element.
The nodal force vector {p} can be expressed as:
{p} = {PX i , Px j , Px k , Py i , Py j , Pyk}) (2) and the nodal displacements vector {?} can be expressed as,
{?} = {u i , u j , u k , v i , v j , v k} (3) which are shown in Figure 5D-14.
TMI-1 UFSAR APPENDIX 5D 5D-11 REV. 18, APRIL 2006 Formulation of Element Stiffness Matrix
Corresponding to the above assumptions, the displacements of a point at (x,y) are written as
linear functions of the coordinates:
u = u i + C 1 x + C 2 y v = v i + C 3 x + C 4 y (4)
The stiffness matrix [k] can be obtained by the following steps:
- a. Substitute the nodal coordinates into Equation(4) to obtain the displacement functions in terms of the nodal displacements.
- b. Substitute results of (a) into the linear strain-displacement equations to obtain the strains in terms of nodal displacements.
- c. Substitute the results of (b) into the stress-strain equation to obtain the stress in terms of nodal displacements.
- d. Find the nodal resultant forces which are statically equivalent to the stresses found in (c) so that the relationships between the nodal forces and the nodal
displacements can be written in a matrix form as shown in Equation (1).
The explicit formulations of the matrices which compose the element stiffness matrix can be found in Reference 7.
Solution of the Total Stiffness Equations After the stiffness matrices are formulated for each of the elements, they are assembled to form a total system stiffness matrix. If the boundary conditions and the loading conditions are properly specified, the tot al stiffness simultaneous equations are readily solvable to give the nodal displacements. Then, the nodal displacements are used to generate the stresses in each element and at each node. The directions and magnitudes of the principal
stresses can also be computed by using the Mohr's circle technique.
Iterative Method and Optimum Over-Relaxation Factor The Gauss-Seidel iterative procedure (Referenc e 8) is applied to solve the stiffness simultaneous equations. The method requires the storage of only the non-zero elements of the
stiffness matrix during the iterative process.
Since the stiffness matrix [k] is positive definite,the iterative method will always converge (Reference 10). The covergence rate is, however, often rather slow. An over-relaxation factor (Reference 9) can be applied which greatly incr eases the convergence rate. The determination of the proper over-relaxation factor is diffi cult and is based mostly on experience and trial and error.
4.3 GEOMETRY
OF BUTTRESSES WITH TENDONS
The pertinent geometry of the Reactor Building as it affects the buttress analysis is shown in Figure 5D-15. In Figure 5D-15a, a sectional view illustrates the vessel with six buttresses.
Figure 5D-15b, shows the arrangement of the hoop tendons on a stretch out of the cylindrical TMI-1 UFSAR APPENDIX 5D 5D-12 REV. 18, APRIL 2006 wall. Figure 5D-15c, shows the dimensions of one buttress and the arrangement of the hoop and vertical tendons.
4.4 THERMAL
LOADS AND TEMPERATURE PROFILES
If the temperature profile in a buttress is known, the difference between the temperature of any point and the reference temperature can be determined. The temperature difference, with the
modulus of elasticity and thermal expansion coeffi cient of the concrete, provides the equivalent constant stress for each finite-element. The constant thermal stress of each element can be
converted to equivalent nodal resultant forces which can be considered as external forces
applied at the nodal points of a non-heated element.
Temperature Profiles The temperature profiles and iso-thermal curves were calculated by a computer program developed by Brunnel (Reference 11) for the s olution of steady state and transient state problems in one, two, or three dimensional Cart esian or cylindrical systems. The governing differential equations were solved by the method of finite differences, using standard relaxation techniques.
The convergence criterion:
T j (step i) - T j (step i-1) <
0.00001 (5)
Y j (step i)
where T j is the temperature at node j, was applied at every node of the difference lattice. To ensure stability of the transient solutions, the time increment was selected according to the
method presented in Reference 11.
The buttress and wall temperatures were determined for the following conditions:
- a. Normal operation in winter
- b. Start-up in winter
- c. Shut-down in winter
- d. Loss of coolant accident during winter
The winter season was chosen because it produces the greatest thermal gradients in the shell
of the Reactor Building.
The mean monthly temperature was chosen due to t he fact that day- night fluctuations had little effect on the interior portion of the concrete. According to Carslow and Jeager, (Reference 12)
the amplitude of the temperature diminishes as:
wP con C p con exp [X
] (6) 2 k con Thus a distance of 12.8 inches from the outer wall of the Reactor Building the effect of the
oscillation is only 0.1 what it is on the outer face, and the middle and inner portion of the vessel
will be unaffected by day-night temperature oscillations.
TMI-1 UFSAR APPENDIX 5D 5D-13 REV. 18, APRIL 2006 For the start-up condition, it was assumed that the inside air temperature of the Reactor Building rose from 60 F to 110 F in one hour, and for shut-down the reverse was used. The
inside film coefficient was taken as 1.490 Btu/hr/sq ft/F for cases a, b, and c above. This was
modified film coefficient, which took into account the steel liner and the still air between the liner
and the concrete:
h m = 1 (7) 1 X steel X air h i +
+ k steel k air Equation 7 neglects the circumferential conduction and heat capacity in the steel liner and air
gap. However, due to the relatively small thi ckness of the liner and the air gap when compared with the thickness of the concrete wall, these effects may be neglected. The preliminary
analyses verified this assumption. The following values, obtained by theoretical analysis (Reference 13), were used in the thermal analysis:
inside film coefficient = 2.0 Btu/hr/sq ft/F
specific heat, C p con = 0.210 Btu/lb/F
thermal conductivity, k con = 10.0 Btu in./hr/sq ft/F
k steel = 300.0 Btu in./hr/sq ft/F
k air = 1.836 Btu in./hr/sq ft/F
The density of the concrete, P con was assumed to be 140 lbs/cu ft.
The above values agree closely with the results from experimental research carried out in Japan (Reference 14).
C p con = 0.200 Btu/lb/F
k con = 11.28 Btu in./hr/sq ft/F
The density of the concrete P con , was assumed to be 146 lb/cu ft.
The thickness of the steel and air gap was:
X air = 1/32 inch
X steel = 3/8 inch
Using Reference 15, the local film coefficient for air at 20 F and a free stream velocity of 20 mph
flowing past a concrete cylinder was calculated to be:
h = 0.56 Btu/hr/ft 2/F As a safe approximation the outside film coefficient was chosen as:
TMI-1 UFSAR APPENDIX 5D 5D-14 REV. 18, APRIL 2006 h? = 1.0 Btu/hr/sq ft/F
The buttress model was replaced by a mesh of 61 nodes, with adiabatic boundaries at the
centerline of the buttress and at the radial mid-section of the wall between two neighboring
buttresses (Figures 5D-16 and 5D-17). Due to symmetry, there is no heat transferred across
these boundaries. For Case d the liner temperature during LOCA was used as the inner
boundary temperature in conjunction with a film coefficient of 5.880 Btu/hr/sq ft/F, which
simulated the conduction through the air gap:
k air h m = (8) X air The transient state temperature analyses were performed at various time steps.
The temperature profiles for Cases a, b, c, and d are shown in Figure 5D-18. Because the most
severe temperature gradients occur for case a and d, only these two cases were considered in
the stress analyses. Iso-thermal curves fo r cases a and d are shown in Figures 5D-19 and 5D-20.
The design pressure inside the Reactor Building during LOCA is equal to 55 psig and the
maximum pressure was initially calculated to occur at approximately 200 seconds after initiation
of the accident. The maximum liner temperature of 280F was initially calculated to occur at
approximately 650 seconds after the initiation of the accident. However, the maximum
temperature gradient, which causes the most severe thermal stresses, was calculated to occur
at 10,000 seconds after initiation of the accident.
Consequently, the determination of the most severe stress state, due to thermal and pressure loads, requires that a number of stress analyses be conducted at various elapsed times after
commencement of the accident. A conservative estimate of the stresses was obtained, by
superimposing directly the maximum values for accident temperature and pressure.
Consequently although subsequent calculations alter t he time at which minimum pressure and liner temperature occur, these changes do not affect this analysis.
4.5 NON THERMAL LOADS
In addition to the thermal loads, there are other loads which act on the structure:
- a. Hoop Tendon Anchor Pressure
The tendon force was assumed to be distributed uniformly up the face of the buttress. The resulting average force (=p a) was applied to each of six nodes.
The location of the nodes approximately corresponded to the position of the
tendon.
- b. Hoop Tendon Pressure in a Radial Direction
Radial pressures, due to the hoop tendons in the wall, were calculated as average stresses in pounds/in.
The two radial pressures applied (refer to Figure 5D-21):
TMI-1 UFSAR APPENDIX 5D 5D-15 REV. 18, APRIL 2006
- 1. From point A to the tangent point T. Over this length the tendon spacing causing radial pressure (i.e., continuous past the buttress) is 33 in., giving
a pressure of p t , and
- 2. For the length from point T to point B (midway between buttresses) the average tendon spacing causing radial pressure is 16 1/2 in., giving a
pressure of 2p
- t.
- c. Internal Pressure Due to LOCA
Internal pressure acting on the inside face of the Reactor Building Wall was considered as a combination of two separate effects:
- 2. Pressure due to expansion of the liner caused by the increase in temperature during the LOCA condition (110F to 280F) based upon a limiting yield strength of 1.2 f
- y.
The above two loads were assumed to act simultaneously during the LOCA condition, giving an equivalent internal pressure of p
- i.
The magnitude and combination of the above non-thermal loads is noted in section 4.7
4.6 MATERIAL
PROPERTIES
Prestressing Steel Wire
The type of prestressing wire to be used in the construction is a stress relieved wire.
Dimensions, properties, etc., of the wire as used in the analysis were:
- a. Size of tendon = 169 wires
- b. Diameter of wires = 1/4 in.
- c. Minimum guaranteed ultimate tensile stress = 240,000 psi
- d. Tendon force @ 70 percent GUTS = 1400 kips
- e. Modulus of elasticity = 29 x 10 6 psi f. Coefficient of friction = 0.16
- g. Coefficient of wobble friction = 0.00030 TMI-1 UFSAR APPENDIX 5D 5D-16 REV. 18, APRIL 2006 Concrete The assumed properties of the concrete used in the analysis were:
- a. Compressive strength @ 28 days = 5000 psi
- b. Modulus of elasticity = 4 x 10 6 psi
- c. Poissons ratio = 0.20
Liner The liner is fabricated from 3/8 inch thick steel plate manufactured to ASTM A-283 Grade C.
The minimum yield stress for this material is 30,000 psi. In determining the equivalent liner pressure LOCA, the yield stress was assumed to be 36,000 psi.
4.7 LOADING
CASES AND COMBINATIONS
The buttress and wall were analyzed for four loading cases. The first two cases considered the
buttress and wall being subjected to the maximum prestress forces, i.e., prestress forces
allowing only those losses that would occur between prestressing and start-up (beginning of
plant life.) The latter two cases considered the buttress and wall being subjected to the
minimum prestress forces, i.e., prestress forc es after 40 year losses (end of plant life).
The load cases were:
- a. Winter normal operation with maximum prestress force.
- b. Winter LOCA with maximum prestress force.
- c. Winter normal operation with minimum prestress force.
- d. Winter LOCA with minimum prestress force.
These cases were considered to be more critical than the summer normal operation/LOCA
conditions. The summer condition results in a smaller temperature differential in the wall which
in turn results in maximum compressive and tensile stresses less than the winter condition.
Hence, the summer condition was not considered in the analyses.
TMI-1 UFSAR APPENDIX 5D 5D-17 REV. 18, APRIL 2006 The numerical value of the applied loads was:
LOAD CASE a b c d
t LT 110 150 110 150
t LP 110 280 110 280
t e 20 20 20 20
p a 6.25 6.73 5.52 6.03
pt A-T 46.07 49.75 39.25 42.2
2ptT-B 92.13 99.50 78.5 84.4
p i 0.0 96.1 0.0 96.1
where:
t LT = Temperature (F) of liner for thermal gradient in wall.
t LP = Temperature (F) of liner for equivalent pressure of liner due to expansion during LOCA.
t e = External air temperature (F)
p a = Prestress force at bearing plate (kips/inch)
pt A-T = Radial prestress force from point A to tangent point T (lbs/inch)
2p T-B = Radial prestress force tangent point T to point B (lbs/inch)
p i = Equivalent internal pressure (LOCA pressure + liner expansion pressure lbs/inch)
Figure 5D-21 shows the application of these loads on the buttress and wall.
Load p a p t , and 2p t for Case a incorporate losses due to elastic shortening of the concrete, friction and creep, shrinkage, and steel relaxation losses that occur between prestressing and
start-up. Case b prestress forces reflect increases from Case a, due to elastic elongation of the
structure during LOCA. Thus, Case a and b reflect the maximum prestress force that will exist
in the structure, at the earliest time in the life of the plant, that a LOCA could occur.
Loads p a , p t , and 2 pt for Cases c and d are similar to Cases a and be respectively, except that the losses due to creep, shrinkage, and steel relaxation are those that are assumed to occur
over 40 years to end of plant life. Case d pres tress forces reflect an increase from Case c due to elastic elongation of the structure during LOCA.
TMI-1 UFSAR APPENDIX 5D 5D-18 REV. 18, APRIL 2006
4.8 STRESS
ANALYSIS
Due to the symmetrical nature of the structur e only one-twelfth of the cylinder, as shown in Figure 5D-16, need be analyzed. Preliminary analy ses showed that beyond a certain point in the wall, the stress state, for all loading cases, was uniform. To conserve computer storage and
give a more accurate solution, it was decided to reduce the size of the model to that shown in
Figure 5D-22. Nodal displacements of both boundaries were specified. The displacements
were obtained form a computer program as developed by Kalnins (Reference 16) using a shell of revolution theory.
There were a total of 672 elements and 384 nodal points, so that 768 simultaneous equations
required solution. For the computations, the Gauss-Seidel iterative solution procedure was used.
4.9 ANALYSIS
OF RESULTS
Figures 5D-23 through 5D-30 show iso-stress curves for the maximum and minimum principal stresses for the four loading conditions outlined in Section 4.7. It should be noted that these
curves represent the algebraic maximum and minim um values of the principal stresses. The do not show the orientation of the stress trajectories.
Inspection of the iso-stress curves shows tha t, as expected, high compressive stresses exist directly beneath the bearing plate for all load conditions. Under winter normal conditions, a
region of tensile bursting stress occurs in the buttress at some distance removed from the
bearing plate; during LOCA, however, these bursting stresses merge with those resulting from the pressure and temperature loads.
The variation of temperature through both the wall and buttress leads to substantial hoop stress
gradients. It should be noted that in the curves for LOCA, the direction of the stresses in the
wall changes in the vicinity of the centerline.
Of primary interest, however, is the high localized t ensile stress that occurs at the buttress-wall re-entrant corner during LOCA, which results fr om both the thermal gradient as well as the local bending arising from the variation of resistance of the buttress and wall to internal pressure.
4.10 CONCLUSIONS
The results indicate that for the most severe loading conditions, the maximum stresses in the
buttress and adjacent wall are within safe limits, except for the localized area at the buttress
wall corner where tensile stresses greater than 6 ?f
? occur. Reinforcement to stabilize cracking in this local region are provided. All other results would indicate that the design method
employed was conservative.
TMI-1 UFSAR APPENDIX 5D 5D-19 REV. 18, APRIL 2006 5 REFERENCES
- 1. Building Code Requirements for Reinforced Concrete , ACI 318-63, American Concrete Institute, Detroit, Michigan, 1963.
- 2. Y. Guyon, Prestressed Concrete , Contractors Record Ltd. London and John Wiley and Sons Inc., New York, 1960.
- 3. F. Leonhardt, Prestressed Concrete Design and Construction , Wilhelm Ernest and Sohn, Berlin-Munich, 1964.
- 4. P. Gergely, M. A. Sozen and C. P. Siess, The Effect of Reinforcement on Anchorage Zone Cracks in Prestressed Concrete Members , Civil Engineering Studies, Structural Research Series No. 271, University of Illinois, Urbana, Illinois, July 1963.
- 5. Anchorage Systems in Prestressed Concrete Pressure Vessels; Anchorage Zone Problems , Oak Ridge National Laboratory, ORNL-TM-2378, January, 1969.
- 6. M. J. Turner, R. W. Clough, H. C. Martin, and L. J. Topp, "Stiffness and Deflection Analyses of Complex Structures," Jl. of Aero. Sci., V. 23, No. 9, Sept 1965, p. 805.
- 7. E. L. Wilson, Finite-Element Analysis of Two Dimensional Structures , Ph. D. Thesis, Univ. of Cal. at Berkeley, June 1963.
- 8. M. G. Salvadori and M. L. Baron, Numerical Methods in Engineering , Prentice-Hall, Englewood Cliffs, N. J., 1961.
- 9. F. G. Lehan, "Simultaneous Equations S olved by Over-Relaxization, "Proc. ASCE, 2nd Conf. on Electronic Computation , Pittsburgh, Pa., USA, Sept. 1960.
- 10. S. H. Crandall, Engineering Analysis , McGraw-Hill Book Co., N. Y., USA, 1956.
- 11. R. D. Brunnel, Computer Program Heating and Heating 2 , Atomic International Division, Report No. A1-64-MEMO-177, North American Aviation, Inc., Detroit, Michigan, USA, 1964.
- 12. H. S. Craslaw and J. C. Jaeger, Conduction of Heat in Solids , 2nd Ed., Oxford, London, 1959.
- 13. Johns-Manville Research and Engineering Center Report E-455-T-266 filed in Robert Emmett Ginna Final Facility Description and Safety Analysis Report, Section 5, Appendix 5B.
- 14. J. Kosaka, Experimental Research of Unsteady-State Heat Transfer Concerned with Nuclear Energy Facilities (Preliminary Report), Division of Nuclear Energy Facilities, Taisei Construction Ltd., 1969.
- 15. F. Kreith, Principles of Heat Transfer , 2nd Ed. International, Scranton, Pa., 1965.
- 16. A. Kalnins, Computer Programs , Crystal River Unit No., 3, Florida Power Corporation, PSAR, Pages 5C-1 and 5C-2.
- 17. J. Zielinski and R. E. Rowe, An Investigation of the Stress Distribution in the Anchorage Zones of Post-Tensioned Concrete Members , Cement and Concrete Association, Research Report NO. 9, London, September 1960.
TMI-1 UFSAR APPENDIX 5D 5D-20 REV. 18, APRIL 2006
- 18. J. Zielinski and R. E. Rowe, The Stress Distribution Associated with Groups of Anchorages in Post-Tensioned Concrete Members , Cement and Concrete Association, Research Report No. 13, London, October 1962.
- 19. P. Gergely, Anchorage Systems in Prestressed Concrete Pressure Vessels; Anchorage Zone Problems , U.S. Atomic Energy Commission, ORNL-TM-2378, January, 1969.
TMI-1 UFSAR APPENDIX 5D 5D-21 REV. 18, APRIL 2006 TABLE 5D-1 BUTTRESS ANCHOR ZONE TRANSVERSE STRESSES IN RADIAL AND VERTICAL PLANES Case Radial Plane
1a 1b 1c 2a
2a'(in) 20.5 20.5 20.5 20.5
2a (in) 52 26 78 32
2a'/2a 0.40 0.79 0.26 0.64
fy[1] 0.27 0.10 0.34 0.15
z (in) [2] 0.80 0.90 0.70 0.9
2b' (in) 20.5 20.5 20.5 20.5
2b=2b'+(2az)tg30 o 44.5 33.5 51.5 37.0 A=2a2b (in
- 2) 2315 870 4150 1182 P (kips) [3] 1400 1400 1400 1400
f bz = P/A (ksi) 0.605 1.610 0.338 1.182 f by = f bz fy (ksi) 0.164 0.161 0.115 0.178 [4]
Zy = 0.3P (1-a'/a) 252 88 311 151 (kips) [5]
Zy (design) = 0.3P 420 420 420 420 (kips) [6]
az (in) 20.8 11.7 27.3 14.4
[1] Guyon, "Prestressed Concrete", Figure 9.
2 [2] Guyon, "Prestressed Concrete", Figure 9.
2 [3] 1600 Kips and 2000 kips anchor forces is not critical.
[4] Changes in tensile splitting stresses with varying "2a" dimension.
[5] Zy = total transverse force Leonhardt "Prestressed Concrete"
- p. 271.
3 [6] For design purpose "(a'/a)" is conservatively assumed to be equal
to zero.
TMI-1 UFSAR APPENDIX 5D 5D-22 REV. 18, APRIL 2006 TABLE 5D-1 (continued)
BUTTRESS ANCHOR ZONE TRANSVERSE STRESSES IN RADIAL AND VERTICAL PLANES
Case Ultimate Load- Vertical Planes Radial Plane 2a 3a 3b 3c
2a' (in) 20.5 20.5 20.5 20.5
2a (in) 32 33 26 78
2a'/2a 0.64
fy [1] 0.15
z (in) [2] 0.9
2b' (in) 20.5
2b=2b' + (2az)tg30 0 37.0 See Case 2a, 1b, and 1c A=2a2b (in
- 2) 1182 P (kips) [3] 2000
f bz = P/A (ksi) 1.690 f by = f bz fy (ksi) 0.254 [4]
Zy=0.3P (1-a'/a) 216 (kips) [5]
Zy (design)=0.3P 600 (kips) [6]
az (in) 14.4
[1] Guyon, "Prestressed Concrete", Figure 9.
2 [2] Guyon, "Prestressed Concrete", Figure 9.
2 [3] 1600 Kips and 2000 kips anchor forces is not critical.
[4] Changes in tensile splitting stresses with varying "2a" dimension.
[5] Zy = total transverse force Leonhardt "Prestressed Concrete" p. 271.
3 [6] For design purpose "(a'/a)" is conservatively assumed to be equal to zero.
TMI-1 UFSAR APPENDIX 5D 5D-23 REV. 18, APRIL 2006 TABLE 5D-2 BUTTRESS ANCHOR ZONE SPALLING (Refer to Figure 5D-13)
M = P {c - e - (c/h) 2 [2h - 3e - c + 2ec
]} h c= varies e = 15.5 in. h = 69 in.
M = P {c - 15.5 - (c/69) 2 [2 x 69 - 3 x 15.5 - c ( 1 - 2x15.5
)]} 69
C (in.) c - 15.5 (c/69) 2 -c( 1 - 2 x 15.5) M/P M m (in. - K) 69
5 -10.5 0.00525 - 2.75 -11.0 10 - 5.5 0.02100 - 5.50 - 7.3 15.5 0 0.05040 - 8.50 - 4.4 20 + 4.5 0.08400 -11.00 - 2.3 25 + 9.5 0.1310 -13.80 -0.7 zero at 27 in. 30 +14.5 0.1890 -16.50 + 0.3 42 +26.5 0.370 -23.10 + 1.2 1680 50 +34.5 0.525 -27.50 + 0.9 (in. - kips) 69 +53.5 1.00 -38.00 0
F t = M m Z = 3 in h = 69 in.
h-z
f s allowable = 10 5 W A s A s = 1.27 sq in. (#10 bars)
W = 0.005
f s = 10 5 0.005/1.27 = 6.28 ksi F T = 1680
= 25.5 K 69 - 3
A s = F T/f s = 25.5/6.28 = 4.05 sq in.
A s provided: #10 @ 12 in. + 1 #10 Loop = 6.03 sq in.