Next Article in Journal
Experimental Comparison between Spherical and Refractive Optics in a Concentrating Photovoltaic System
Next Article in Special Issue
Uncertainty and Sensitivity Analysis of the In-Vessel Hydrogen Generation for Gen-III PWR and Phebus FPT-1 with MELCOR 2.2
Previous Article in Journal
Risks and Safety of CO2 Transport via Pipeline: A Review of Risk Analysis and Modeling Approaches for Accidental Releases
Previous Article in Special Issue
The EC MUSA Project on Management and Uncertainty of Severe Accidents: Main Pillars and Status
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Influence of Vertical Equivalent Damping Ratio on Seismic Isolation Effectiveness of Nuclear Reactor Building

1
Faculty of Infrastructure Engineering, Dalian University of Technology, Dalian 116024, China
2
Nuclear and Radiation Safety Center, Ministry of Ecology and Environment, Beijing 100082, China
*
Author to whom correspondence should be addressed.
Energies 2021, 14(15), 4602; https://doi.org/10.3390/en14154602
Submission received: 17 June 2021 / Revised: 10 July 2021 / Accepted: 13 July 2021 / Published: 29 July 2021
(This article belongs to the Special Issue Design and Safety Issues of Nuclear Plants and Installations)

Abstract

:
This paper aimed at evaluating the influence of different vertical equivalent damping ratios of a 3-dimensional combined isolation bearing (3D-CIB) as regards seismic response and isolation effectiveness. A comparative study of the seismic response in terms of acceleration floor response spectra (FRS), peak acceleration, displacement response of the nuclear reactor building, and dynamic response of the 3D-CIB was carried out. The results showed that: (1) the horizontal FRS is slightly influenced by the vertical equivalent damping ratio of 3D-CIB, whereas the increase of the vertical equivalent damping ratio has a significant effect on reducing the vertical FRS; (2) the peak vertical acceleration increased with the decrease in the vertical equivalent damping ratios of 3D-CIB and the difference of peak accelerations calculated by the damping ratio of 20 and 25% is within 10%; (3) the increase of the vertical equivalent damping ratio is capable of reducing the horizontal displacement and the rocking effect of the superstructure, and effectively controlling the vertical displacement amplitude; and (4) the vertical equivalent damping ratio of 3D-CIB has a slight effect on its axial force. Consequently, it is demonstrated that the increase of the vertical equivalent damping ratio is advantageous for isolation effectiveness. From the view of displacement control, it is suggested that the 3D-CIB with the vertical an equivalent damping ratio of 15~20% is appropriate and acceptable.

1. Introduction

Earthquakes are one of the most potentially destructive natural disasters, causing great loss of property and life. In the traditional seismic design of structures, seismic resistance is enhanced by increasing the strength of the structural components. In the past, several kinds of new passive and active techniques have been developed to reduce the disaster of earthquakes, with base isolation a popular seismic protection technology whereby some bearings with low stiffness are interposed between the structure and the ground [1,2,3,4]. In recent years, innovative methods have been developed to increase the efficiency and safety of structures against earthquakes. Since the seismic base isolation system can decrease the total energy transmitted from the ground to the superstructure, it has been studied extensively. Different types of isolation bearings have been developed as effective strategies to protect infrastructure from seismic hazards [5,6,7,8,9]. For instance, steel-rubber isolators with rubber cores [10,11], natural rubber bearings equipped with high toughness steel ring dampers [12], natural rubber bearings equipped with steel and shape memory alloys dampers [13], new rectangular rubber isolators with square rubber cores [14], and laminate rubber bearings (LRB) [15]; the above isolators are used for horizontal isolation.
Moreover, earthquakes are one of the most severe natural disasters that threaten the safety of nuclear power plants (NPPs). Although isolation technology has been widely deployed for industrial and civil buildings, bridges, and certain classes of infrastructure, the application of base isolation technology in NPPs is limited. Among all NPPs currently in commercial operation worldwide, only Koeberg in South Africa and Cruas in France have adopted base isolation technology, in which neoprene rubber bearings, including flat sliders in some installations, have been used as the horizontal base isolation [16,17]. The main reasons for the limited use of seismic isolation in NPPs, its development, and future research needs have been presented and identified [18,19]. Based on many new theories and techniques, the base isolation of NPPs in the horizontal direction has been widely studied. These conventional base isolation systems clearly show the benefits of seismic isolation in the horizontal direction. However, the vertical ground motions transmitted through the isolation system were either unchanged or even magnified [17,20,21,22].
It is known that the vertical seismic design of NPPs is as necessary and important as the horizontal directions, especially for the high intensity zone. If 3D seismic isolation is achieved, it would substantially enhance the economy and safety of NPPs. Hence, it is beneficial that a few isolation systems are invented to provide vertical isolation for NPPs. Various attempts have been made to provide enhanced protection against the vertical component of ground motion by using a complete three-dimensional (3D) seismic isolation system, or adding a vertical isolation system in combination with a horizontal isolation component or the vertical isolation of the main components associated to the horizontal base-isolation structure. Several vertical and 3D isolation systems which have potential application in modern nuclear facilities have been examined [23]. Devices such as the rolling seal-type air spring [24,25], hydraulic isolation systems [26], and coned disk springs [27,28,29,30] have been developed to provide isolation in the vertical direction. In order to achieve the effectiveness of 3D seismic isolation, they are proposed to be combined with a horizontal isolation component.
Although 3D seismic base isolation has not yet been applied to the engineering practice of NPPs, its development has become an attractive area of research. Since coned disk springs have high vertical bearing capacity and also low vertical stiffness, this method was selected to decouple the superstructure from the vertical component of ground motion. In this study, a new 3D combined isolation bearing (3D-CIB), composed of LRB coupled with the combined disk spring bearing (CDSB), is proposed. LRB is used as horizontal isolation component, with CDSB used as the vertical isolation component at the bottom. In this way, the combination of LRB and CDSB can achieve 3D seismic isolation. In contrast to the traditional bearings of coned disk springs with only one disk spring column [27,28,29,30], CDSB is composed of one main column and several auxiliary columns of coned disk springs. The advantage of CDSB is that the constitution can be flexibly adjusted according to the requirements of the bearing capacity and vertical stiffness. It should be noted that CDSB itself has a high damping ratio, i.e., 0.2–0.3, due to the large number of coned disk springs stacked [31]. Since the equivalent damping ratio of CDSB is an important parameter affecting isolation performance, it is necessary to conduct a performance test to measure the vertical equivalent damping ratio for the practical application of CDSB. To investigate the influence of different vertical equivalent damping ratios of 3D-CIB on the seismic response and isolation effectiveness, four vertical equivalent damping ratios were assumed in this paper. Subsequently, a comparative study of the seismic response in terms of acceleration floor response spectra (FRS), peak acceleration, displacement response of the nuclear reactor building, and the dynamic response of 3D-CIB was carried out.

2. Design of 3D Combined Isolation Bearing (3D-CIB)

2.1. Design Theory and Parameters of CDSB

A typical plan and section view of the CDSB is illustrated in Figure 1. The main column of disk spring is located in the center, and the auxiliary columns are uniformly distributed on a large circle with a certain radius. It is worth noting that the number of auxiliary columns m can be flexibly adjusted according to the bearing capacity and vertical stiffness requirements of the CDSB. The main and auxiliary columns are all equipped with central guide tubes, and the bottom of each central guide tube is rigidly connected to the bottom connection plate. The central guide tube is installed in the stack in order to prevent excessive side slip. The upper part of each main and auxiliary column is connected to a guide sleeve, which is integrated with the top connection plate. These components including the bottom connection plate, central guide tubes, guide sleeves, and top connection plate together ensure that the horizontal stiffness of the CDSB is sufficient.
The compound combination adopted for CDSB is achieved through stacking the coned disk springs of the same shape in parallel, and then by piling up this set more and more in several series. The main column was composed by stacking nc coned disk springs in parallel, and then ic sets in series, which made up nc × ic disk springs in total. Similarly, each auxiliary column was composed of the stack of ns coned disk springs in parallel, and then is sets in series. The characteristic parameters of a coned disk spring with bearing surface, including the outer diameter (D), inner diameter (d), thickness of the bearing surface (t), and inner coned height (h), are shown in Figure 2. Accordingly, for the main and auxiliary coned disk springs, the above parameters are denoted as Dc, dc, tc, hc and Ds, ds, ts, hs, respectively.
The principle of the coned disk spring was first presented by Almen and Laszlo [32] and then recommended for the existing national standard for coned disk springs [33]. Considering the main disk spring as an example, when the vertical load Fc acts on the inner circumference of the coned disk spring, the deformation fc occurs:
F C = 4 E 1 μ 2 t C 4 K 1 D C 2 K 4 2 f C t C K 4 2 h C t C f C t C h C t C f C 2 t C + 1
where E is the elasticity modulus, µ is Poisson’s ratio, and K1 and K4 are calculation coefficients.
When the main disk spring is flattened, that is, fc = hc, Equation (1) yields the critical load of a coned disk spring FCc, as follows:
F C c = F f C = h C = 4 E 1 μ 2 t C 3 h C K 1 D C 2 K 4 2
The critical load of a single auxiliary disk spring FSc can be obtained in a similar manner. Finally, the critical load of the CDSB can be expressed as:
F C D S B c = F V C + m F V S = n C F C c + m n S F S c
It is assumed that the weight of the superstructure borne by each CDSB is defined as F. Then, according to the deformation characteristics of the compound combination, the preloading initial deformation can be obtained as fCDSB = icfc = isfs.
According to the stiffness definition of the coned disk spring, under the initial deformation of the main disk spring fc, its vertical stiffness KC can be derived:
K C = d F C d f C = 4 E 1 μ 2 t C 3 K 1 D C 2 K 4 2 K 4 2 h C t C 2 3 h C t C f C t C + 3 2 f C t C 2 + 1
Further, the vertical stiffness of the main and auxiliary columns of the CDSB can be inferred as follows:
K V C = n C K C / i C K V S = n S K S / i S
Finally, the total vertical stiffness of the CDSB is obtained as follows:
K V   CDSB = K V C + m K V S
Owing to the large weight of the nuclear reactor building, a total of 245 bearings were set, and the average gravity load borne by each bearing was 2843 kN. Since many factors such as strength, critical load, and deformation should be considered in the design of the coned disk spring, after the trial calculation the geometrical parameters of the coned disk spring with the bearing surface were finally confirmed, as shown in Table 1.
According to Equation (2), the critical loads of the single main and auxiliary disk spring are FCc = 623.1 kN and FSc = 155.8 kN, respectively. For the compound combination of the CDSB in this paper, m = 6, nc = ns = 4, is = 2ic = 28. According to Equation (3), the critical load of the CDSB is FCDSBc = 6231.3 kN.
Based on the preloading initial displacement of the main and auxiliary disk spring, it can be obtained from Equation (4) that the vertical stiffness values of the single main and auxiliary coned disk spring are KC = 7.238 × 107 N/m and KS = 1/2KC, respectively. Finally, according to Equation (5) and Equation (6), the vertical stiffness of the CDSB can be obtained KVCDSB = 5.170 × 107 N/m.

2.2. Design Theory of 3D-CIB

3D-CIB consists of a CDSB to provide vertical isolation, which is supported by a LRB to provide horizontal isolation. The CDSB and LRB are stacked in series using high-strength bolts, as illustrated in Figure 3. The LRB has small horizontal stiffness and high vertical stiffness, whereas the CDSB exhibits high horizontal stiffness, and its vertical stiffness is significantly lower than that of the LRB. Therefore, it can be inferred that the horizontal stiffness ( K H 3 D CIB ) of the 3D-CIB only considers the horizontal stiffness of the LRB ( K H L R B ), and the vertical stiffness ( K V 3 D CIB ) is calculated according to a series combination of LRB and CDSB. In summary, the horizontal and vertical stiffness of the 3D-CIB is as follows:
K H 3 D CIB = K H   LRB K V 3 D CIB = K V   LRB K V   CDSB / K V   LRB + K V   CDSB
In this study, the LRB with an additional lead core to enhance the energy dissipation capacity was used, with the vertical stiffness K V L R B = 3.392 × 10 9   N / m , initial stiffness before yield K u = 1.32 × 10 7   N / m , stiffness after yield K d = 1.2 × 10 6   N / m , yield load Q d = 113   kN , equivalent damping ratio   ζ e q = 0.25 , and ultimate shear displacement D u l t = 40   c m . According to Equation (7), the vertical stiffness values of the 3D-CIB can be obtained K V 3 D CIB = 5.092 × 10 7   N / m .

3. Basic Data for Numerical Analysis

3.1. 3D Base-Isolated Nuclear Reactor Building

In this study, the nuclear reactor building includes three parts: raft foundation, containment, and internal structure. The outer diameter, height, and thickness of the containment were 51.2, 69.0, and 1.15 m, respectively. The commercial finite element software ANSYS was employed to establish the 3D model, as shown in Figure 4. The raft foundation, with a thickness of 3.5 m (EL. 5.2–8.7 m), was meshed by the element type of Solid 185. Containment was meshed by the element type of Shell 181, and the element sizes of the cylinder and dome were approximately 1.6 m × 1.7 m and 1.6 m × 1.2 m, respectively. To simplify the model and reduce the calculation time, a lumped mass stick model of the internal structure was developed by the element type of Beam 188 and Mass 21. Generally, the properties of the lumped mass stick model were back-calculated from an analysis of the real 3D structure, and can represent the dynamic characteristics of the real 3D structure. The total weight of the nuclear reactor building was 70,966 tons. A total of 245 isolation bearings were uniformly distributed beneath the raft foundation, as shown in Figure 5. The element type of COMBIN40, in which both the bilinear strengthening model and the effects of viscous damping are considered, was adopted to simulate the horizontal behavior of 3D-CIB. COMBIN40 includes two spring constants (K1 and K2), damping coefficient ( C ), mass constant ( M ), gap size ( G A P ), and limiting sliding force ( F S L I D E ). The relationship among initial stiffness Ku, K1, and K2 is that Ku = K1 + K2. Mass constant ( M ) and gap size ( G A P ) are set to zero for isolators. The element type of COMBIN14, whose real constants include the spring constant ( K ) and damping coefficient ( C V ), was adopted to simulate the vertical characteristics of 3D-CIB. The geometry and mechanical models of COMBIN14 and COMBIN40 are shown in Figure 6. Thus, one element of COMBIN14 and two elements of COMBIN40 were used to simulate an isolation bearing in which one node of each element was shared with the bottom node of the raft foundation, and the other node was a fixed constraint (i.e., restraining degrees of freedom (DoF) of translation along the three orthogonal directions) for the boundary condition.

3.2. Input Ground Motion

In this paper, the three translational components of the artificial acceleration time histories, plotted in Figure 7, are the input ground motions. The corresponding response spectra for a damping ratio of 5%, plotted in Figure 8, are the design response spectra of the Regulatory Guide 1.60 [34]. The peak ground acceleration in both the horizontal and vertical directions is 2.943 m/s2, i.e., 0.3 g (g = 9.81 m/s2), the total duration and time step is 28 s and 0.01 s, respectively. All three components were statistically independent of each other. Consequently, according to the Regulatory Guide 1.92 [35], the seismic response of the structure was directly obtained from the algebraic summation.

4. Results and Discussion

Firstly, a modal analysis of both a non-isolated and 3D base-isolated nuclear reactor building was performed. Subsequently, to investigate the influence of different vertical equivalent damping ratios of 3D-CIB on the seismic response and isolation effectiveness, four vertical equivalent damping ratios of 10, 15, 20, and 25% for 3D-CIB were assumed in this paper. Finally, the seismic isolation effectiveness of 3D-CIB with different vertical equivalent damping ratios was compared and discussed.

4.1. Modal Analysis

Firstly, the main frequencies and mode shapes that reflect the dynamic characteristics of the non-isolated and 3D base-isolated nuclear reactor building can be obtained. The main frequencies are provided in Table 2. In this table, “ U X , U Y , U Z ” denote the DoFs of translation in the, X, Y, and Z directions, respectively. “ R O T Z ” denotes the DoF of torsion around the Z-axis. “FIXED” and “3D-CIB” denote the results calculated for the non-isolated and 3D base-isolated structure with 3D-CIB, respectively. As can be seen in Table 2, the main frequencies of the 3D isolated structure for the DoFs of both translation and torsion are significantly reduced compared to that of the non-isolated structure. The comparison of these frequencies confirms that the design objective is reached.

4.2. Comparison of Floor Response Spectra (FRS)

The seismic demands for safety-related secondary mechanical, electrical, and piping systems in NPPs are generally using FRS. A comparative study of the FRS was focused on the representative nodes, including containment nodes located at the ring crane (NODE 12) and dome top (NODE 13), secondary top node (NODE 3), and top node (NODE 4) of the internal structure, as shown in Figure 4. It is noticed that the horizontal FRS for NODE 3 and NODE 4 of the internal structure is extracted only in the X direction, owing to the structural symmetry in the horizontal plane. For brevity, the 5% damped FRS for the horizontal and vertical directions are plotted in Figure 9 and Figure 10, respectively. The legends of “FIXED” and “25%” denote the results calculated for the non-isolated structure and 3D base-isolated structure of 3D-CIB with a vertical equivalent damping ratio of 25%, respectively. Moreover, a comparison of the spectral acceleration peak value ( S m a x ), reduction rate of S m a x ( δ ), and zero period acceleration (ZPA) of the horizontal and vertical FRS are summarized in Table 3 and Table 4, respectively.

4.2.1. Analysis of FRS in the Horizontal Direction

It can be observed from Figure 9 that there are no apparent differences among the spectral shapes and values of the 3D isolated structures with different vertical equivalent damping ratios, and the first peak states all occurred near the horizontal isolation frequency. Additionally, the spectral values in high-frequency range (i.e., above 1.0–2.0 Hz) are significantly reduced, compared to those of the non-isolated structures. Overall, the spectral peak values and ZPA slightly increased with a decrease of vertical equivalent damping ratio. Especially for containment nodes, it can be seen from Table 3 that the spectral peak values calculated by 15, 20, and 25% of vertical equivalent damping ratios are reduced by approximately 5, 10, and 15% respectively, compared with the corresponding results of the damping ratio 10%. For internal structural nodes, the spectral peak values and ZPA calculated by different equivalent damping ratios are almost the same. In summary, it can be concluded that the horizontal FRS was slightly influenced by the vertical equivalent damping ratio of 3D-CIB.

4.2.2. Analysis of the FRS in the Vertical Direction

As shown in Figure 10, the spectral shapes calculated by different vertical equivalent damping ratios are similar, and all of them have peak values at 2.2 Hz. Furthermore, the peak values and the spectral values to the right of spectral peak increased obviously with a decrease in the vertical equivalent damping ratio. As summarized in Table 4, the spectral peak values calculated by the vertical equivalent damping ratios of 15, 20, and 25% are approximately 20, 35, and 40% lower than the corresponding results of the damping ratio of 10%. Hence, it can be concluded that the increase of the vertical equivalent damping ratio has a significant effect on reducing the vertical FRS.
In addition, the spectral values near the peak frequency and to the left of the spectral peak are significantly magnified in comparison to those of the non-isolated structure. In contrast, the spectral values of the frequency range greater than approximately 3.5 Hz are obviously reduced.

4.3. Comparison of Peak Acceleration Response

In this section, the acceleration response of the nodes at different elevations of the containment and internal structure, as shown in Figure 4, were extracted. The absolute peak accelerations versus the height of containment and internal structure are plotted in Figure 11 and Figure 12, respectively. In general, it appeared that both the horizontal and vertical peak accelerations of the 3D isolated structure did not vary much with height, whereas they were obviously amplified along the height of the non-isolated structure.
As shown in Figure 11 and Figure 12, the peak values of the horizontal acceleration response calculated by different vertical equivalent damping ratios of 3D-CIB are equivalent in the middle and lower parts of the containment and internal structures, while the peak values slightly increased with the decrease of the equivalent damping ratio in the upper structure. Take the dome top node for example, compared with the vertical equivalent damping ratio of 10%, the calculated values of 15, 20, and 25% are reduced by 5.9, 10.5, and 13.8%, respectively. For the internal structure, the peak value of the Y-direction is little affected by the vertical equivalent damping ratio from the bottom to the top.
In the vertical direction, it is evident that the peak vertical acceleration increased with a decrease in the vertical equivalent damping ratios of 3D-CIB. For the top nodes of containment and internal structure, compared to the non-isolated structure, the substantial reduction ratios calculated by the vertical equivalent damping ratios of 10, 15, 20, and 25% were 63.5, 68.2, 71.0, and 72.6%, and 22.2, 28.9, 35.7, and 41.2%, respectively. Moreover, compared with the peak acceleration of the vertical equivalent damping ratio of 25%, the corresponding results of 10, 15, and 20% are increased by 33.1, 16.1, and 5.9%, and 39.2, 20.9, and 9.4%, respectively. For the sub top node of the containment, compared to the non-isolated structure, the substantial reduction ratios calculated by the vertical equivalent damping ratios of 10, 15, 20, and 25% were 21.0, 31.4, 33.8, and 36.0%, respectively. Compared with the vertical equivalent damping ratio of 25%, the calculated values of 10, 15, and 20% are increased by 24.5, 7.3, and 3.5%, respectively. From the above comparison, it can be inferred that the difference of peak accelerations calculated by the damping ratio of 20 and 25% is within 10%.

4.4. Comparison of Relative Displacement Response

The displacement response of the 3D isolated and non-isolated structures relative to the ground is presented in this section. Note that all the displacement responses were obtained by considering the action of a gravity superimposed earthquake. The peak displacements of the nodes at different elevations of the containment and internal structures are plotted in Figure 13 and Figure 14, respectively. The legends “25%+” and “25%−” denote the displacement response along the coordinate axis in the positive and negative directions, which is calculated by 3D isolated structure of 3D-CIB with a vertical equivalent damping ratio of 25%.
It is observed that the horizontal peak displacement along the positive direction of the middle and upper part of containment and internal structure, calculated by different vertical equivalent damping ratios of isolation bearings, is slightly different, increasing slightly with the decrease of damping ratio. However, the horizontal peak displacement of containment and internal structure along the negative direction increases with a decrease in the vertical equivalent damping ratio. For the dome top node, the peak horizontal displacement of the dome top in the negative direction of the X-axis for 10, 15, 20, and 25%, as shown in Figure 13a, was 14.15, 13.40, 12.72, and 12.10 cm, respectively. Compared with the vertical equivalent damping ratio of 25%, the peak values calculated by 20, 15, and 10% are increased by 5.1, 10.8, and 14.5%, respectively. In addition, it is clear that the peak horizontal displacement of the 3D isolated reactor building is undergoing some amplification along the height due to rocking effect. For the containment, the rocking ratio (i.e., the horizontal displacement of the dome top relative to the raft foundation divided by the height of dome top) for 10, 15, 20, and 25%, was 1/837, 1/891, 1/943, and 1/997, respectively. Thus, compared with the vertical equivalent damping ratio of 25%, the rocking ratios calculated by 10, 15, and 20% are increased by 19.1, 11.9, and 5.7%, respectively.
As shown in Figure 13c and Figure 14c, for the 3D isolated structure, the peak vertical displacement along the positive direction of the Z-axis is less than zero. Under earthquake action, the vertical vibration displacement calculated by different vertical equivalent damping ratios increased with the decrease of damping ratio. The vertical displacement vibration amplitudes of nodes at the top edge and center of the raft foundation as shown in Figure 5, calculated by vertical equivalent damping ratios of 10, 15, 20, and 25%, are 8.77, 8.00, 7.41, and 6.97 cm, and 4.47, 3.94, 3.53, and 3.19 cm, respectively. Consequently, compared with the vertical equivalent damping ratio of 25%, the displacement amplitude calculated by 10, 15, and 20% are increased by 25.7, 14.7, and 6.3%, and 39.9, 23.3, and 10.7%, respectively.

4.5. Dynamic Response of 3D-CIB

4.5.1. Hysteretic Response of 3D-CIB

As there are 245 bearings in the model, it is unfeasible to illustrate the response results of all the isolators. Therefore, the response of the center 3D-CIB is discussed in this section. For the 3D-CIB with different vertical equivalent damping ratios, its shape of hysteretic curves of shear and axial force-displacement response is similar, so the hysteretic shear and axial force-displacement response for the central 3D-CIB with damping ratio of 20% is shown in Figure 15. It can be seen that the hysteretic shape in the horizontal direction is different from the idealized hysteresis behavior of the bilinear mechanical model of LRB, due to the coupling effect of LRB and CDSB. For the axial force-displacement response, the hysteretic behavior performed ellipse, because the vertical stiffness of 3D-CIB is linear.

4.5.2. Vertical Displacement Response of 3D-CIB

A total of 17 bearings were arranged along the Y-axis beneath the raft foundation, as shown in Figure 5. The peak vertical displacements of 3D-CIB with different vertical equivalent damping ratios are shown in Figure 16. It is evident that the peak vertical displacements are all less than zero, which means that they are all in a compressed state under the action of the gravity superimposed earthquake. The peak vertical displacement along the negative direction and the displacement vibration amplitude of bearings with different vertical damping ratio increased with the decrease of equivalent damping ratios. Furthermore, the vertical displacement amplitude of the 3D-CIB near the edge was larger than that near the center, and the smaller the vertical damping ratio was, the more obvious the amplification effect. It was demonstrated that the rocking effect of the raft foundation occurred during the earthquake.
Figure 17 presents a comparison of the vertical displacement time histories of the center and edge bearings, as shown in Figure 5. Under the action of gravity, the initial displacements of the 3D-CIB at the center and edge are −5.56 and −5.65 cm, respectively. Then, under the gravity superposition earthquake action, the peak vertical displacements of the bearings with vertical equivalent damping ratios of 25, 20, 15, and 10% are −7.13, −7.37, −7.66, and −8.02 cm, and −9.26, −9.51, −9.81, and −10.14 cm, respectively. Compared with the peak vertical displacement of the bearing with a damping ratio of 25%, the corresponding results of 10, 15, and 20% are increased by 12.4, 7.3, and 3.3%, and 9.5, 5.9, and 2.6%, respectively. Further, the deformation of the single main coned disk spring of the center and edge 3D-CIB with vertical equivalent damping ratios of 25, 20, 15, and 10% was 5.10, 5.27, 5.47, and 5.73 mm, and 6.61, 6.79, 7.01, and 7.25 mm, corresponding to 0.58 h C , 0.60 h C , 0.62 h C , 0.65 h C , and 0.75 h C , 0.77 h C , 0.80 h C , and 0.82 h C , respectively. It is worth noting that the above deformation of the single main coned disk located at the center and edge was for all scenarios all less than h C .

4.5.3. Comparison of the Axial Force of 3D-CIB

Under the action of the gravity superimposed earthquake, the peak axial forces of the bearings along the Y-axis are extracted as shown in Figure 18, in which the legend “3D-CIB-G” denotes the initial axial force of the bearings under the gravity action of the superstructure. It should be noted that the axial force of the 3D-CIB is positive in compression and negative in tension. As plotted in Figure 18, the initial axial forces of the bearings are uniform. Additionally, it is clear that the axial forces of bearings are all greater than zero, which means that they are all in compression, without the potential uplift. When the earthquake vibrated downward, the peak axial force of bearings decreased with the increase of the vertical equivalent damping ratio. In addition, the vertical force of the isolation bearing near the edge is obviously larger than that at the center. Taking the two side edge bearings along the positive and negative direction as an example, the peak axial force of bearings with damping ratios of 25, 20, 15, and 10% is 4.953, 5.021, 5.063, and 5.200 × 106 N, and 4.838, 4.921, 5.040, and 5.189 × 106 N, and all less than the critical load F C D S B = 6231.3   kN . Moreover, compared with the damping ratio of 25%, the peak axial force of bearings of 10, 15, and 20% are reduced by 5.0, 2.2, and 1.4%, and 7.2, 2.4, and 1.7%, respectively. To sum up, the above results illustrated that the vertical equivalent damping ratio of 3D-CIB had a slight effect on axial force.

5. Conclusions

This study aimed at evaluating the influence of different vertical equivalent damping ratios of 3D-CIB on the seismic response and isolation effectiveness, with four vertical equivalent damping ratios assumed. A comparison analysis of seismic response in terms of the FRS, peak acceleration, displacement response of a nuclear reactor building, and dynamic response of a 3D-CIB was performed. It was found that: (1) the horizontal FRS is slightly influenced by the vertical equivalent damping ratio of the 3D-CIB, whereas the increase of the vertical equivalent damping ratio has a significant effect on reducing the vertical FRS; (2) the peak vertical acceleration increased with the decrease in the vertical equivalent damping ratios of the 3D-CIB, and the difference of peak accelerations calculated by the damping ratios of 20 and 25% is within 10%; (3) the increase of the vertical equivalent damping ratio is capable of reducing the horizontal displacement and the rocking effect of the superstructure, and effectively controlling the vertical displacement amplitude. When the vertical equivalent damping ratio was 15 and 20%, the horizontal displacement and the rocking ratio of the superstructure was 5~12% higher than the corresponding results of 25%, and the vertical compression displacement of the 3D-CIB is increased by within 8%; (4) the vertical equivalent damping ratio of the 3D-CIB has a slight effect on its axial force. Consequently, it was demonstrated that the increase of vertical equivalent damping ratio is advantageous for isolation effectiveness. From the point of view of displacement control, it is suggested that a 3D-CIB with the vertical equivalent damping ratio of 15~20% is appropriate and acceptable. Since rocking motion is inevitable for the 3D base-isolated structure, and excessive rocking motion is undesired and should be suppressed for NPPs, the study of a rocking suppression device will be carried out in the future.

Author Contributions

Conceptualization, X.Z. and G.L; methodology, X.Z. and J.L.; software, X.Z.; validation, X.Z.; writing—original draft preparation, X.Z. and J.L.; writing—review and editing, G.L. and R.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Science and Technology Major Project of advanced technology research on the structural health inspection and evaluation of nuclear power plants (Grant No. 2018ZX06002008). This financial support is gratefully acknowledged.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Shrimali, M.; Jangid, R.; Shrimali, M.; Jangid, R. Earthquake response of isolated elevated liquid storage steel tanks. J. Constr. Steel Res. 2003, 59, 1267–1288. [Google Scholar] [CrossRef]
  2. Petti, L.; Giannattasio, G.; De Iuliis, M.; Palazzo, B. Small scale experimental testing to verify the effectiveness of the base isolation and tuned mass dampers combined control strategy. Smart Struct. Syst. 2010, 6, 57–72. [Google Scholar] [CrossRef] [Green Version]
  3. Castaldo, P.; Palazzo, B.; Alfano, G.; Palumbo, M.F. Seismic reliability-based ductility demand for hardening and softening structures isolated by friction pendulum bearings. Struct. Control. Health Monit. 2018, 25, e2256. [Google Scholar] [CrossRef]
  4. Kelly, J. Analysis of the run-in effect in fiber-reinforced isolators under vertical load. J. Mech. Mater. Struct. 2008, 3, 1383–1401. [Google Scholar] [CrossRef] [Green Version]
  5. De Domenico, D.; Ricciardi, G. Earthquake-resilient design of base isolated buildings with TMD at basement:Application to a case study. Soil Dyn. Earthq. Eng. 2018, 113, 503–521. [Google Scholar] [CrossRef]
  6. De Domenico, D.; Gandelli, E.; Quaglini, V. Effective base isolation combining low-friction curved surface sliders and hysteretic gap dampers. Soil Dyn. Earthq. Eng. 2020, 130, 105989. [Google Scholar] [CrossRef]
  7. Markou, A.; Stefanou, G.; Manolis, G.D. Stochastic response of structures with hybrid base isolation systems. Eng. Struct. 2018, 172, 629–643. [Google Scholar] [CrossRef]
  8. De Luca, A.; Guidi, L.G. Base isolation issues in Italy: Integrated architectural and structural designs. Soil Dyn. Earthq. Eng. 2020, 130, 105912. [Google Scholar] [CrossRef]
  9. De Luca, A.; Guidi, L.G. State of art in the worldwide evolution of base isolation design. Soil Dyn. Earthq. Eng. 2019, 125, 105722. [Google Scholar] [CrossRef]
  10. Rahnavard, R.; Thomas, R.J. Numerical evaluation of steel-rubber isolator with single and multiple rubber cores. Eng. Struct. 2019, 198, 109532. [Google Scholar] [CrossRef]
  11. Rahnavard, R.; Craveiro, H.D.; Napolitano, R. Static and dynamic stability analysis of a steel-rubber isolator with rubber cores. Structures 2020, 26, 441–455. [Google Scholar] [CrossRef]
  12. Sheikhi, J.; Fathi, M.; Rahnavard, R. Natural rubber bearing incorporated with high toughness steel ring dampers. Structures 2020, 24, 107–123. [Google Scholar] [CrossRef]
  13. Sheikhi, J.; Fathi, M.; Rahnavard, R.; Napolitano, R. Numerical analysis of natural rubber bearing equipped with steel and shape memory alloys dampers. Structures 2021, 32, 1839–1855. [Google Scholar] [CrossRef]
  14. Altabani, D.; Hejazi, F.; Rashid, R.; Aziz, F. Development of new rectangular rubber isolators for a tunnel-form structure sub-jected to seismic excitations. Structures 2021, 32, 1522–1542. [Google Scholar] [CrossRef]
  15. Radkia, S.; Rahnavard, R.; Tuwair, H.; Gandomkar, F. Investigating the effects of seismic isolators on steel asymmetric structures considering soil-structure interaction. Structures 2020, 27, 1029–1040. [Google Scholar] [CrossRef]
  16. Coladant, C. Seismic isolation of nuclear power plants-EDF’s philosophy. Nucl. Eng. Des. 1991, 127, 243–251. [Google Scholar] [CrossRef]
  17. Frano, R.L. Benefits of Seismic Isolation for Nuclear Structures Subjected to Severe Earthquake. Sci. Technol. Nucl. Install. 2018, 8017394. [Google Scholar] [CrossRef] [Green Version]
  18. Whittaker, A.S.; Kumar, M. Seismic Isolation of Nuclear Power Plant. Nucl. Eng. Technol. 2014, 46, 569–580. [Google Scholar] [CrossRef] [Green Version]
  19. Whittaker, A.S.; Sollogoub, P.; Kim, M.K. Seismic isolation of nuclear power plants: Past, present and future. Nucl. Eng. Des. 2018, 338, 290–299. [Google Scholar] [CrossRef]
  20. Frano, R.L.; Forasassi, G. Isolation systems influence in the seismic loading propagation analysis applied to an innovative near term reactor. Nucl. Eng. Des. 2010, 240, 3539–3549. [Google Scholar] [CrossRef]
  21. Zhao, C.; Chen, J. Numerical simulation and investigation of the base isolated NPPC building under three directional seismic loading. Nucl. Eng. Des. 2013, 265, 484–496. [Google Scholar] [CrossRef]
  22. Wang, D.; Zhuang, C.; Zhang, Y. Seismic response characteristics of base-isolated AP1000 nuclear shield building subjected to beyond-design basis earthquake shaking. Nucl. Eng. Technol. 2018, 50, 170–181. [Google Scholar] [CrossRef]
  23. Zhou, Z.; Wong, J.; Mahin, S. Potentiality of Using Vertical and Three-Dimensional Isolation Systems in Nuclear Structures. Nucl. Eng. Technol. 2016, 48, 1237–1251. [Google Scholar] [CrossRef] [Green Version]
  24. Suhara, J.; Matsumoto, R.; Oguri, S.; Okada, Y.; Inoue, K.; Takahashi, K. 3-D seismic isolation device with rolling seal type air spring. In Proceedings of the 18th International Conference on Structural Mechanics in Reactor Technology (SMiRT18), Beijing, China, 7–12 August 2005; pp. 3381–3391. [Google Scholar]
  25. Suhara, J.; Tamura, T.; Okada, Y.; Umeki, K. Development of Three Dimensional Seismic Isolation Device With Laminated Rubber Bearing and Rolling Seal Type Air Spring. In Proceedings of the Fatigue, Fracture and Damage Analysis; ASME International, Vancouver, BC, Canada, 5–9 August 2002; pp. 43–48. [Google Scholar]
  26. Kashiwazaki, A.; Shimada, T.; Fujiwaka, T.; Moro, S. Study on 3-dimensional base isolation system applying to new type power plant reactor. In Proceedings of the Transaction of the 17th International Conference on Structural Mechanics in Reactor Technology SMiRT17, Prague, Czech Republic, 17–22 August 2003; pp. 2295–2302. [Google Scholar]
  27. Kitamura, S.; Okamura, S.; Takahashi, K. Experimental Study on Vertical Component Seismic Isolation System With Coned Disk Spring. In Proceedings of the Volume 7: Operations, Applications, and Components; ASME International, Denver, CO, USA, 17–21 July 2005; pp. 175–182. [Google Scholar]
  28. Kitamura, S.; Morishita, M.; Moro, S. Study on Vertical Component Seismic Isolation System With Coned Disk Spring. Aging Manag. Licens. Renew. 2004, 2, 21–28. [Google Scholar] [CrossRef]
  29. Morishita, M.; Kitamura, S.; Moro, S.; Kamishima, Y.; Somaki, T. Study on 3-dimensional seismic isolation system for next generation nuclear power plant-vertical component isolation system with coned disk spring. In Proceedings of the 13th World Conference on Earthquake Engineering (ASME), Vancouver, BC, Canada, 1–6 August 2004. Paper No. 620. [Google Scholar]
  30. Wang, W.; Wang, X. Tests, model, and applications for coned-disc-spring vertical isolation bearings. Bull. Earthq. Eng. 2019, 18, 357–398. [Google Scholar] [CrossRef]
  31. Zhao, Y.; Su, J.; Lu, M.; Yu, Y. Mechanical performance test of combined 3-dimensional isolation bearing for 3-dimensional base-isolated buildings. In Proceedings of the Tenth Pacific Conference on Earthquake Engineering. Building an Earthquake-Resilient Pacific Earthquake, Sydney, Australia, 6–8 November 2015. [Google Scholar]
  32. Almen, J.O.; Laszlo, A. The Uniform-Section Disc Spring, Transactions of American Society of Mechanical Engineers. 1936. Available online: https://www.shotpeener.com/library/pdf/1935000.pdf (accessed on 8 May 2021).
  33. China National Standardization Management Committee. Disc Spring (GB/T1972-2005), 2005. Available online: https://www.chinesestandard.net/PDF/English.aspx/GBT1972-2005 (accessed on 8 May 2021).
  34. U.S. Nuclear Regulatory Commission. RG 1.60, Rev.1. Design Response Spectra for Seismic Design of Nuclear Power Plants. Washington, DC, 1973. Available online: https://www.nrc.gov/docs/ML0037/ML003740207.pdf (accessed on 8 May 2021).
  35. U.S. Nuclear Regulatory Commission. RG 1.92, Rev.2. Combining Modal Responses and Spatial Components in Seismic Response Analysis. Washington, DC, 2006. Available online: https://www.nrc.gov/docs/ML0532/ML053250475.pdf (accessed on 8 May 2021).
Figure 1. Typical plan and section view of the combined disk spring bearing (CDSB). (a) Plan view; (b) Section view.
Figure 1. Typical plan and section view of the combined disk spring bearing (CDSB). (a) Plan view; (b) Section view.
Energies 14 04602 g001
Figure 2. Structure size of single disk spring with bearing surface.
Figure 2. Structure size of single disk spring with bearing surface.
Energies 14 04602 g002
Figure 3. Section view of the 3D combined isolation bearing (3D-CIB).
Figure 3. Section view of the 3D combined isolation bearing (3D-CIB).
Energies 14 04602 g003
Figure 4. Finite element model of a nuclear reactor building.
Figure 4. Finite element model of a nuclear reactor building.
Energies 14 04602 g004
Figure 5. Layout of discrete 3D-CIBs beneath the raft foundation.
Figure 5. Layout of discrete 3D-CIBs beneath the raft foundation.
Energies 14 04602 g005
Figure 6. The geometry and mechanical model of the spring-damper element. (a) COMBIN14; (b) COMBIN40.
Figure 6. The geometry and mechanical model of the spring-damper element. (a) COMBIN14; (b) COMBIN40.
Energies 14 04602 g006
Figure 7. The acceleration time-histories of seismic motion. (a) In the horizontal direction X; (b) In the horizontal direction Y; (c) In the vertical direction Z.
Figure 7. The acceleration time-histories of seismic motion. (a) In the horizontal direction X; (b) In the horizontal direction Y; (c) In the vertical direction Z.
Energies 14 04602 g007
Figure 8. The acceleration response spectra of the input ground motion.
Figure 8. The acceleration response spectra of the input ground motion.
Energies 14 04602 g008
Figure 9. Comparison of the FRS of containment and internal structures in the horizontal direction. (a) NODE 12_X; (b) NODE 13_X; (c) NODE 12_Y; (d) NODE 13_Y; (e) NODE 3_X; (f) NODE 4_X.
Figure 9. Comparison of the FRS of containment and internal structures in the horizontal direction. (a) NODE 12_X; (b) NODE 13_X; (c) NODE 12_Y; (d) NODE 13_Y; (e) NODE 3_X; (f) NODE 4_X.
Energies 14 04602 g009aEnergies 14 04602 g009b
Figure 10. Comparison of the FRS of containment and internal structures in the vertical direction. (a) NODE 12; (b) NODE 13; (c) NODE 3; (d) NODE 4.
Figure 10. Comparison of the FRS of containment and internal structures in the vertical direction. (a) NODE 12; (b) NODE 13; (c) NODE 3; (d) NODE 4.
Energies 14 04602 g010
Figure 11. Comparison of the peak acceleration response of the containment. (a) In the horizontal direction X; (b) In the horizontal direction Y; (c) In the vertical direction Z.
Figure 11. Comparison of the peak acceleration response of the containment. (a) In the horizontal direction X; (b) In the horizontal direction Y; (c) In the vertical direction Z.
Energies 14 04602 g011
Figure 12. Comparison of the peak acceleration response of the internal structure. (a) In the horizontal direction X; (b) In the horizontal direction Y; (c) In the vertical direction Z.
Figure 12. Comparison of the peak acceleration response of the internal structure. (a) In the horizontal direction X; (b) In the horizontal direction Y; (c) In the vertical direction Z.
Energies 14 04602 g012
Figure 13. Comparison of the peak displacement response of the containment. (a) In the horizontal direction X; (b) In the horizontal direction Y; (c) In the vertical direction Z.
Figure 13. Comparison of the peak displacement response of the containment. (a) In the horizontal direction X; (b) In the horizontal direction Y; (c) In the vertical direction Z.
Energies 14 04602 g013aEnergies 14 04602 g013b
Figure 14. Comparison of the peak displacement response of the internal structure (a) In the horizontal direction X; (b) In the horizontal direction Y; (c) In the vertical direction Z.
Figure 14. Comparison of the peak displacement response of the internal structure (a) In the horizontal direction X; (b) In the horizontal direction Y; (c) In the vertical direction Z.
Energies 14 04602 g014aEnergies 14 04602 g014b
Figure 15. Hysteretic shear and axial force-displacement response of central 3D-CIB. (a) In the horizontal direction X; (b) In the vertical direction Z.
Figure 15. Hysteretic shear and axial force-displacement response of central 3D-CIB. (a) In the horizontal direction X; (b) In the vertical direction Z.
Energies 14 04602 g015
Figure 16. Comparison of the peak vertical displacement response of the 3D-CIBs.
Figure 16. Comparison of the peak vertical displacement response of the 3D-CIBs.
Energies 14 04602 g016
Figure 17. Comparison of the vertical relative displacement response of the center and edge 3D-CIBs. (a) At center; (b) At edge.
Figure 17. Comparison of the vertical relative displacement response of the center and edge 3D-CIBs. (a) At center; (b) At edge.
Energies 14 04602 g017
Figure 18. Comparison of the peak axial force of the 3D-CIBs along Y-axis.
Figure 18. Comparison of the peak axial force of the 3D-CIBs along Y-axis.
Energies 14 04602 g018
Table 1. Geometrical parameters of the coned disk spring.
Table 1. Geometrical parameters of the coned disk spring.
TypeOuter
Diameter (D)
Inner
Diameter (d)
Thickness of Bearing Surface (t)Inner Cone Height (h)
Main disk spring320106.717.68.8
Auxiliary disk spring16053.38.84.4
(Unit: mm).
Table 2. Comparison of the natural frequencies of vibration.
Table 2. Comparison of the natural frequencies of vibration.
ORDERMODE3D-CIBFIXED
1 U X 0.7054.541
2 U Y 0.7104.540
3 R O T Z 1.0869.455
4 U Z 2.08413.750
(Unit: Hz).
Table 3. Comparison of peak values and ZPA of the FRS in the horizontal direction.
Table 3. Comparison of peak values and ZPA of the FRS in the horizontal direction.
TYPENODE12_XNODE12_YNODE13_XNODE13_Y
S m a x   ( g ) δ Z P A   ( g ) S m a x   ( g ) δ Z P A   ( g ) S m a x   ( g ) δ Z P A   ( g ) S m a x   ( g ) δ Z P A   ( g )
25%0.77612.4%0.2320.85514.1%0.2640.97713.2%0.2941.10215.0%0.349
20%0.8098.7%0.2450.89610.0%0.2721.0229.2%0.3051.16010.6%0.361
15%0.8454.6%0.2590.9435.2%0.2801.0714.9%0.3201.2255.6%0.374
10%0.886-0.2750.995-0.2901.126-0.3401.297-0.389
FIXED4.624-0.8714.928-0.9635.896-1.1715.955-1.179
Table 4. Comparison of peak values and ZPA of the FRS in the vertical direction.
Table 4. Comparison of peak values and ZPA of the FRS in the vertical direction.
TYPENODE3NODE4NODE12NODE13
S m a x   ( g ) δ Z P A   ( g ) S m a x   ( g ) δ Z P A   ( g ) S m a x   ( g ) δ Z P A   ( g ) S m a x   ( g ) δ Z P A   ( g )
25%1.21942.2%0.3301.22042.2%0.3311.46040.3%0.4101.28942.5%0.359
20%1.38734.2%0.3621.38934.2%0.3631.66931.8%0.4251.47934.1%0.382
15%1.66421.1%0.3991.66621.1%0.4001.97019.5%0.4411.77221.0%0.420
10%2.109-0.4602.111-0.4602.448-0.5092.243-0.479
FIXED2.444-0.5522.539-0.5702.426-0.6476.776-1.329
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhu, X.; Li, J.; Lin, G.; Pan, R. Influence of Vertical Equivalent Damping Ratio on Seismic Isolation Effectiveness of Nuclear Reactor Building. Energies 2021, 14, 4602. https://doi.org/10.3390/en14154602

AMA Style

Zhu X, Li J, Lin G, Pan R. Influence of Vertical Equivalent Damping Ratio on Seismic Isolation Effectiveness of Nuclear Reactor Building. Energies. 2021; 14(15):4602. https://doi.org/10.3390/en14154602

Chicago/Turabian Style

Zhu, Xiuyun, Jianbo Li, Gao Lin, and Rong Pan. 2021. "Influence of Vertical Equivalent Damping Ratio on Seismic Isolation Effectiveness of Nuclear Reactor Building" Energies 14, no. 15: 4602. https://doi.org/10.3390/en14154602

APA Style

Zhu, X., Li, J., Lin, G., & Pan, R. (2021). Influence of Vertical Equivalent Damping Ratio on Seismic Isolation Effectiveness of Nuclear Reactor Building. Energies, 14(15), 4602. https://doi.org/10.3390/en14154602

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop