ML19114A202
ML19114A202 | |
Person / Time | |
---|---|
Site: | Peach Bottom |
Issue date: | 04/08/2019 |
From: | Exelon Generation Co |
To: | Office of Nuclear Reactor Regulation |
Shared Package | |
ML19114A265 | List:
|
References | |
Download: ML19114A202 (7) | |
Text
PBAPS UFSAR EXHIBIT V - SEISMIC ANALYSIS OF REACTOR VESSEL V.1 ABSTRACT AND RESULTS The maximum seismic shears and moments in the Peach Bottom reactor pressure vessel and vessel internals due to the maximum credible earthquake (MCE) have been determined. The maximum credible earthquake response spectra used in this analysis is shown in Figure C.3.2. The results for those components that develop significant seismic loads are given in the following table:
SIGNIFICANT SEISMIC SHEARS AND MOMENTS MCE RPV Support Skirt Moment 6,320 kip-ft RPV Support Skirt Shear 205 kips Fuel Assy. Moment 0.89 kip-ft/assy CRD Housing Moment 2.97 kip-ft each Total Fuel Shear on Top Guide 133 kips Shroud Support Moment 6,219 kip-ft Shroud Support Shear 287 kips RPV Stabilizer Force 267 kips V.2 ANALYSIS V.2.1 Mathematical Model The Peach Bottom system was modeled with 40 discrete lumped masses connected by massless springs. Relatively lightweight components such as the jet pumps, steam dryer, in-core housing, in-core guide tube, and spargers were not modeled separately; only their mass effects being accounted for. This simplification is made to limit the complexity of the analysis and is valid since the dynamic effects of these components on the rest of the system is slight.
Should the seismic response of these components be of interest, the response of the reactor vessel at their points of attachment to the vessel would then be used as input forcing functions.
Because of the solidity of the site foundation, the reactor building was modeled as fixed at its base.
V.2.2 Mass and Elastic Properties To develop equivalent lumped masses, the system components are visualized as consisting of vertical segments. The lumped mass, located in general at the center of mass of each segment, represents the sum of: (1) the segment mass, (2) components in the segment not modeled separately (such as the dryer mass being lumped into a vessel mass), (3) the enclosed water mass effects discussed in detail in paragraph V.2.3.
APPENDIX K K.V-1 REV. 21, APRIL 2007
PBAPS UFSAR The system flexibility matrix is determined by using the elastic properties of cross-sectional moment of inertia, shear area, length, Young's Modulus, and Poisson's ratio for each component segment between lumped masses to derive individual segment matrices which are then combined as discussed in paragraph V.2.6.
Shearing deformations are accounted for in the derivation of the segment flexibility coefficients.
V.2.3 Hydrodynamic Mass In order to account for the effects of the core cooling water on the dynamic characteristics of the system, hydrodynamic mass (or apparent mass) terms have been included in the system mass matrix.
The hydrodynamic mass manifests itself as a dynamic coupling between the real masses and therefore appears in the mass matrix at diagonal as well as off-diagonal locations. These hydrodynamic masses are added algebraically to the lumped metal masses to form the total mass matrix.
V.2.4 Natural Frequencies and Mode Shapes The overall system flexibility matrix is derived by using the matrix force method which is described in detail in reference 1.
Given the flexibility matrix, f, the stiffness matrix, K, of the system can be found by inversion. From the mass matrix, (non-diagonal in this case), and the stiffness matrix, the undamped natural frequencies and normalized mode shapes are computed.
V.2.5 Equivalent Modal Damping Damping coefficients for individual component elements are given as follows:
MCE Coefficient Component (% Critical Damping)
Reactor Vessel 2.0 Vessel Support Skirt 2.0 Shroud 2.0 Shroud Head and Separators 2.0 Fuel 7.0 Control Rod Guide Tubes 2.0 CRD Housings 3.5 Shield and Pedestal 5.0 Reactor Building 5.0 Modal damping values are then computed for each mode on the basis of modal-displacement and mass.
APPENDIX K K.V-2 REV. 21, APRIL 2007
PBAPS UFSAR V.2.6 Earthquake Response The equation of motion of the system in matrix form is as follows:
M X + Y + C X + K X = O (V.1) where:
M = mass matrix, n x n (this includes the hydro-dynamic mass)
X = column vector of displacement relative to ground (n x 1)
C = damping matrix (n x n)
K = stiffness matrix (n x n)
Y = column vector of ground accelerations (n x 1)
( ) = second derivative with respect to time
(*) = first derivative with respect to time Upon removing the ground acceleration vector to the right side of Equation (V.1), the equation reduces to the classical form:
M X + CX + K X = M Y (V.2)
To uncouple equation (V.2) set X = q, = modal matrix (V.3)
Equation (V.2) then becomes M q + C q + K q = M Y (V.4)
Pre-multiplying by the transpose of and using the orthogonality conditions:
T M q + T C q + T K q = T M Y (V.5) q + T C q + T Kq = T M Y (V.6)
This procedure for uncoupling the equation of motion by using the modal matrix of the undamped system requires that damping in the system be small. It is further assumed that the damping matrix C is such that TC is a diagonal matrix. The elements of this diagonal matrix are the modal damping values. With these assumptions, Equation (V.6) may be written as:
APPENDIX K K.V-3 REV. 21, APRIL 2007
PBAPS UFSAR qi + 2ii qi + i 2qi = S Ug (+) (V.7) i = 1, 2, -----N where:
th i = damping ratio for the i mode expressed as percent of critical damping th i = i natural angular frequencies of the system th TMI Si = modal participation factor for the i mode = i u Ug = earthquake ground acceleration time history T= th i transpose of the i mode shape Iu = (n x 1) column vector whose elements are all unity The system of one degree of freedom equation represented by Equation (V.7) subjected to the initial conditions qi (0) = 0 qi (0) = 0 determines the modal response qi (+). Using Equation (V.7), the maximum seismic modal response (displacement, acceleration, or load response) can be determined for each natural frequency of interest from the response spectrum curve. Response spectrum curves are essentially plots of the maximum responses of single-degrees-of-freedom systems described by Equation (V.7) with Si =
1.0 against various natural periods or frequencies. Having found the maximum modal responses, qi, i = l, ...m, the maximum physical th displacement for the i mode is given by:
X = Si q (V.8) i i i where:
APPENDIX K K.V-4 REV. 21, APRIL 2007
PBAPS UFSAR X i = X 1i i = 1i X 2i 2i O
O X ni Ni As shown in reference 2, the maximum physical response for each mass can be taken to be the square root of the sums of the squares of each of the maximum responses for each mode; i.e.,
1/2 m _
2 (xi ) max X ij j = 1 i = 1,. . . , N th Similarly, the maximum load response for the i mode if found from:
L = S X i i where:
L 1i Li =
L 2i Where n is the number of rows in the S matrix O
L ni S = Stress or load matrix APPENDIX K K.V-5 REV. 21, APRIL 2007
PBAPS UFSAR The maximum load response is again taken to be the square root of the sums of the squares of each of the maximum responses for each mode; i.e.,
m _ 2 1/ 2 (Li ) max = L ij j = 1 i = 1, . . . n APPENDIX K K.V-6 REV. 21, APRIL 2007
PBAPS UFSAR K.V SEISMIC ANALYSIS OF REACTOR VESSEL REFERENCES
- 1. Leckie, F. A. and Pestel, E. C., Matrix Methods in Elasto-Mechanics, McGraw-Hill, 1963.
- 2. Clough, R. W. "Earthquake Analysis by Response Spectrum Superposition," Bull, Seismological Society of America, Vol 52, No. 3, July, 1962.
APPENDIX K K.V-7 REV. 21, APRIL 2007