Seismic Damage Assessment of Reinforced Concrete Containment Structures
Seismic Damage Assessment of Reinforced Concrete Containment Structures
Seismic Damage Assessment of Reinforced Concrete Containment Structures
ABSTRACT
This paper presents a procedure for assessing seismic damage of concrete containment structures using the
nonlinear time-history numerical analysis. For this purpose, two kinds of damage index are introduced at finite
element and structural levels. Nonlinear finite element analysis for the containment structure applies PSC shell
elements using a layered approach leading to damage indices at finite element and structural levels, which are then
used to assess the seismic damage of the containment structure. As an example of such seismic damage assessment,
seismic damages of the containment structure of Wolsong I nuclear power plant in Korea are evaluated against 30
artificial earthquakes generated with a wide range of PGA according to US NRC regulatory guide 1.60. Structural
responses and corresponding damage index according to the level of PGA and nonlinearity are investigated. It is also
shown that the containment structure behaves elastically for earthquakes corresponding to or lower than DBE.
KEY WORDS: earthquake, damage assessment, prestressed concrete, shell, finite element, layered approach,
nonlinear, damage index, containment structure, artificial earthquake, stress, prestressed tendon.
INTRODUCTION
Appropriate and safe performance under seismic occurrence is required for containment structures in nuclear
power plants to avoid completely radioactive leakage from nuclear reactor and power supply equipment. Pertaining to
such requirement, the Korea Institute of Nuclear Safety (KINS) is operating the Earthquake Monitoring Center, which
records seismic ground motions at nuclear power plant sites, evaluates Korean earthquake characteristics and assesses
seismic damage in nuclear structures. The center has conducted a multi-year research project along with Seoul
National University to establish a seismic damage assessment system which estimates probable seismic damages of
concrete containment structures by performing inelastic time-history analysis [1]. In the assessment procedure, two
kinds of damage index are introduced at finite element and structural levels. Nonlinear finite element analysis for the
containment structure applies PSC shell elements using a layered approach leading to damage indices at finite element
and structural levels, which are then used to assess the seismic damage of the containment structure. As an example of
such seismic damage assessment, seismic damages of the containment structure of Wolsong I nuclear power plant in
Korea are evaluated against 30 artificial earthquakes generated with a wide range of PGA according to US NRC
regulatory guide 1.60. Wolsong I is a nuclear power plant in the southern part of Korea, which has been designed with
respect to the Canadian seismic design codes with 0.2g Design Basis Earthquake(DBE) and 0.1g Site Design
Earthquake(SDE).
Since seismic damages in structures are caused by inelastic deformations, the proposed seismic damage
assessment makes use of nonlinear time history analysis so as to account for material nonlinearities and earthquake
characteristics. An improved Newmark method that is, the HHT-α method, is selected as direct integration method
because of its stability and accuracy.
Most of the containment structures in Korea being made of PSC, the seismic damage assessment system performs
on PSC containment structures, adopting a special nonlinear PSC shell finite element using layered approach and
appropriate inelastic material models. The accuracy of the system to describe correctly the behaviour of PSC
containment structures is validated by means of comparative study between analytical and experimental results at
element and structural levels. Introducing appropriate inelastic material models and including also a nonlinear
prestressed tendon model, the special PSC shell element improves the reinforced concrete shell element using a four-
node quadrilateral thin flat shell finite element with 6 DOFs per node proposed by Kim et al. (2001) [2].
1
layered element formulation, the shell is divided into several paneled layers and two-dimensional constitutive models
were applied to take into account material nonlinearities. The constitutive matrix can be rotated from the local axes to
the global axes and added to the constitutive matrix for the concrete, or it can be used to define the properties as an
overlaying element, adding stiffness to the connected nodes. This technique allows for the addition of any number of
additional layers of steel with different orientations to be added to the finite element model. More detailed development
is presented in [6].
DAMAGE INDEX
Damage index is introduced for the seismic damage assessment of the containment structure in order to quantify
numerically the degree of damage. The concept of damage index can provide the means to quantify damage and relate
it to costs and other consequences such as potential risk after earthquake. Hence, damage index can play an important
role in retrofit decision-making and disaster planning in earthquake region. Two kinds of damage index are introduced
at finite element and structural levels. A simple lumped mass model is used to represent global damage level of the
containment structure as presented in [6]. The damage index at finite element level is explained in the following.
N 2 fc represents the number of loading cycles until concrete reaches the fatigue failure, ε cu is the ultimate strain of
concrete and ε cs is the principle compression strain of concrete at each analysis step.
based on the fatigue model of steel proposed Miner[9]. N 2 fr represents the number of loading cycles until steel reaches
the fatigue failure, ε tu is the ultimate strain of steel and ε ts is the principle tension strain of steel at each analysis step.
2
Table 1. Damage Index at Finite Element Level
The structure selected for seismic damage assessment is the PSC containment structure of Wolsong I nuclear
power plant(Fig. 1a). Wolsong I NPP has been designed with respect to the Canadian seismic design codes with 0.2g
Design Basis Earthquake(DBE) and 0.1g Site Design Earthquake(SDE). The structure, built in 1983, is the object of
contention concerning seismic risks due to its proximity with the Yangsan geologic fault. As shown in Fig. 1b, the top
of the structure is constituted by an upper dome and a lower dome on which dousing water has been poured. The
thickness of the perimeter wall is 1.08m and the one of the ring beam is 1.36m. The equipment opening (3.95×4.27m)
is located on the western face at 9m height from the ground.
(a) (b)
Fig. 1 (a) Wolsong I containment structure, (b) Typical vertical section of the structure
3
Response Spectrum of Artificial Acceleration Time Response Spectrum of Artificial Acceleration Time Response Spectrum of Artificial Acceleration Time
History (Horizontal, NS) History (Horizontal, EW) History (Vertical)
Acceleration (g)
Acceleration (g)
Acceleration (g)
Artificial Acceleration Time History (Horizontal, NS) Artificial Acceleration Time History (Vertical)
Artificial Acceleration Time History (Horizontal, EW)
0.30 0.30
0.30
0.20 0.20 0.20
Acceleration (g)
Acceleration (g)
Acceleration (g)
0.10 0.10 0.10
0.00 0.00 0.00
Push-over Analysis
30000
25000
20000
Reaction(KN)
15000
10000
5000
0
0 1 2 3 4 5 6 7 8
Top Disp.(cm)
D.I.s 0.1g 0.2g 0.4g 0.6g 0.8g 1.0g 1.2g 1.4g 1.6g
Mean -0.130 -0.052 0.105 0.257 0.407 0.569 0.766 1.008 1.219
St. Dev. 0.006 0.011 0.023 0.041 0.058 0.087 0.111 0.137 0.181
Max. -0.119 -0.030 0.152 0.344 0.541 0.793 1.024 1.292 1.578
Min. -0.143 -0.079 0.051 0.173 0.305 0.448 0.538 0.750 0.915
4
Table 3 Summary of maximum damage index at finite element level for 30 seismic analyses results
D.I.e 0.1g 0.2g 0.4g 0.6g 0.8g 1.0g 1.2g 1.4g 1.6g
Mean 0.112 0.115 0.133 0.162 0.227 0.311 0.406 0.506 0.564
St. Dev. 0.000 0.00147 0.00379 0.0117 0.03689 0.06216 0.05720 0.07823 0.09987
Max. 0.112 0.118 0.139 0.200 0.304 0.464 0.542 0.679 0.746
Min. 0.112 0.112 0.125 0.145 0.172 0.209 0.309 0.340 0.356
Max. Max.
1.4 0.7
Mean Mean
1.2 Min. Min.
0.6
1.0
0.5
Max(D.I.e)
0.8
0.4
D.I.s
0.6
0.3
0.4
0.2
0.2
0.0 0.1
-0.2 0.0
0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6
PGA(g) PGA(g)
(a) (b)
Fig. 4 (a) Damage index at structural level and (b) maximum damage index at finite element level according to PGA
Fig. 5 depicts the location distribution of maximum damage index at finite element level around the equipment
opening. The values in brackets stand for the number of occurrences of maximum damage of the concerned element
among the 30 seismic analyses. Among all the elements in the neighbourhood of the equipment opening, the lower
part tends to be the weakest under seismic occurrence. In addition, damage occurrence tends to move toward the
foundation as the PGA increases. Such trend may be explained by the low ratio of height to diameter of the
containment structure.
(3)
5
CONCLUSION
A seismic damage assessment procedure that makes it possible to assess the actual seismic resistance capacity and
the damage level of containment structures has been presented. Damage indices at finite element and structural levels
have been computed to provide quantitative assessment of damages that structures may suffer under seismic event. The
damage index at element level was used to detect the heaviest local damage and damage index at structural level was
used to provide global damage state of the structure.
Resulting from such assessment, it appeared that the structural damage index and the maximum damage index at
finite element level increased with the PGA. In the same way, the corresponding standard deviations increased as
nonlinearity dominance of the seismic behaviour grows also with the PGA. Lower elements located in the
neighbourhood of the equipment opening have been verified to be the weakest under seismic occurrence. Damage
occurrence tends to move toward the foundation as the PGA increases.
Recalling that appropriate and safe performance under seismic occurrence is required for containment structures in
nuclear power plants to avoid completely radioactive leakage from nuclear reactor and power supply equipment, results
obtained from such seismic damage assessment procedure are believed to supply valuable and reliable data not only for
seismic risk assessment but also for the prediction of seismic damage in nuclear containment structures.
ACKNOWLEDGEMENT
The authors wish to thank the Korea Institute of Nuclear Safety(KINS) for providing funding for this research. The
authors are also grateful to the Korea Earthquake Engineering Research Center(KEERC) at Seoul National University
granted by the Korea Science and Engineering Foundation(KOSEF).
REFERENCES
1 Seoul National University, Development of a Seismic Damage Assessment Program for Nuclear Power Plant
Structures, Korea Institute of Nuclear Safety (2002)
2 Kim, T. H., Lee, K. M., and Shin, H. M., Nonlinear Analysis of Reinforced Concrete Shells Using Layered
Elements with Drilling Degree of Freedom. In: ACI Structural Journal. V. 99, No. 4 (2002), pp.418-426.
3 Zienkiewicz, O.C. and Taylor, R.L., The Finite Element Method, Vol. 2-Solid and Fluid Mechanics, Dynamics and
Non-linearity. McGraw Hill Book Co., 4th ed. (2001)
4 Lanheng Jin, Analysis and Evaluation of a Shell Finite Element with Drilling Degree of Freedom. Masters thesis,
University of Maryland at College Park (1994)
5 Hilton, E. and Owen, D.R.J., Finite Element Software for Plates and Shells. Pineridge Press Ltd., Swansea, UK.
(1984)
6 Hyun, C. H., Koh, H. M., Cho, H. H., Lee, S. K., and Choi, K. R., Development of Seismic Damage Assessment
System for Nuclear Power Plant Structures in Korea, to be printed in Proc. 17th Structural Mechamics in Reactor
Technology, No. K564, Prague, Chech Republic.
7 Maekawa, K., Irawan, P., and Okamura, H., Three-Dimensional Constitutive Laws of Reinforced Concrete. In:
International Conference on Applied Concrete Mechanics APCOM, Seoul/1996
8 Priestley, M.J.N., Seible, F. and Calvi, G.M., Seismic Design and Retrofit of Bridges. John Wiley & Sons, Inc.
(1996)
9 Miner, M.A., Cumulative Damage in Fatigue. In: Journal of Applied Mechanics. V.67, No.3 (1945) pp.A159-A164.
10 United State Nuclear Regulatory Commition, Design Response Spectra for Seismic Design of Nuclear Power
Plants, US NRC Reg. Guide, 1.60 (1973)