Abstract
The sequential ground motion effects on the dynamic responses of reinforced concrete containment (RCC) buildings with typical isolators are studied in this paper. Although the base isolation technique is developed to guarantee the security and integrity of RCC buildings under single earthquakes, seismic behavior of base-isolated RCC buildings under sequential ground motions is deficient. Hence, an ensemble of as-recorded sequential ground motions is employed to study the effect of including aftershocks on the seismic evaluation of base-isolated RCC buildings. The results indicate that base isolation can significantly attenuate the earthquake shaking of the RCC building under not only single earthquakes but also seismic sequences. It is also found that the adverse aftershock effect on the RCC can be reduced due to the base isolation applied to the RCC. More importantly, the study indicates that disregarding aftershocks can induce significant underestimation of the isolator displacement for base-isolated RCC buildings.
1. Introduction
Historical earthquakes show that mainshock-aftershock sequence-type ground motions are formed when frequent aftershocks occur after a strong mainshock [1]. Under the United States Geological Survey (USGS), a Mw8.6 earthquake occurred in Indonesia on April 11, 2012. After that, a great number of strong aftershocks happened and the largest magnitude attained 8.2. After the Tohoku-Oki earthquake with magnitude Mw9.0 occurring on March 11, 2011, in Japan, a large number of aftershocks struck the same area and the greatest one reached magnitude Mw7.9. Due to the potentially devastating effect of aftershocks, a great number of researchers attempted to investigate the effects of aftershocks on different structural systems, including the single-degree-of-freedom systems [2–10], reinforced concrete structures [11–18], steel structures [19–22], and wood-frame buildings [23–25]. As the last defense against nuclear leakage, the reinforced concrete containment (RCC) building subject to sequential ground motions is researched by [26]. It reveals that sequential ground motions can have adverse effects on the dynamic responses of the RCC including acceleration, displacements, and accumulated damage.
To deal with the variability of earthquakes, different types of isolation systems have been proposed and adopted for residential structures through the past twenty years. The isolation system has been applied to buildings [27, 28], bridges [29], and other types of structures, particularly in seismic areas with high intensity, such as Japan and California, USA. The isolators have also been employed to protect primary structures and nonstructural components for nuclear power plant (NPP). Using properly designed isolation systems and detailed three-dimensional numerical models, it is found by some researchers that seismic isolation can significantly reduce the seismic responses of NPP structures and nonstructural components [30–34]. To demonstrate the advantage of isolators for NPPs, other researchers adopted simple lumped-mass stick models [35–37]. However, these researches were mainly carried out by considering only mainshocks. Under mainshock-aftershock sequences, the effectiveness of base isolation for NPP structures and facilities may be decreased [38]. Moreover, in order to better guarantee the security and integrity of NPP structures and facilities, it is necessary and important to investigate the dynamic response of NPPs with base isolation and quantify the adverse aftershock effect and damaging aftershock sequences.
This paper is to study the effect of sequential ground motions on the seismic behavior of the RCC building with lead rubber bearing (LRB) isolators. To achieve this aim, finite element models for fixed-base and base-isolated RCC buildings are established employing the software ABAQUS 6.10 [39]. A group of as-recorded sequential ground motions is selected from the Pacific Earthquake Engineering Research (PEER) strong motion database [40]. The seismic analysis results including acceleration and displacements are compared between the RCC with and without isolators. Besides, the damaging aftershock sequence is also investigated and given in this manuscript.
2. Structural Models
2.1. The Fixed-Base RCC Building
Figure 1 shows a sectional view of the analyzed RCC building. The inside and outside diameters of the cylinder are 37 m and 38.8 m, respectively. The height of the cylinder and the total height of the containment building are 49.38 m and 61.4 m, respectively. The wall thickness of the dome is 0.8 m. The elastic modulus, Poisson’s ratio, density, tensile strength, and compressive strength of the containment concrete are assumed to be 36000 MPa, 0.2, 2400 kg/m3, 2.85 MPa, and 40 MPa, respectively. The yield strength of the steel is 350 MPa. The finite element model for the RCC building shown in Figure 2 is set up using ABAQUS 6.10. Since the damaged plasticity model [41, 42] reflects the characteristics of cracking and crushing of concrete and is widely used in concrete structures, it is adopted here as the constitutive relation for concrete. In this model, the uniaxial strength functions are composed of plastic deformation and degradation of stiffness. Therefore, the tensile damage () and the compressive damage () are addressed by normal strain at any point as shown in Figure 3. A bilinear restoring force is adopted for the rebars in the arranged direction, with the second gradient 1/100 of the elastic stiffness supplied by a nuclear power group. The soil-structure interaction is not taken into account because the RCC building is located at a rock soil site. More detailed information about the RCC building can be found in [26].



(a)

(b)
2.2. The RCC Building with Isolators
Because the lead rubber bearing (LRB) system has all desirable features for base isolation and is widely used in engineering applications, it is adopted as isolation systems in this study. LRB is formed by force-fitting the lead plug into a preformed hole in the low damping elastomeric bearing as shown in Figure 4. The steel plates force the lead plug to deform in shear. Performance of LRB is maintained during repeated strong earthquakes with proper durability and reliability. LRB not only supplies the needed damping but also provides higher elastic stiffness. The seismic behavior for LRB is shown in Figure 5, in which is the elastic stiffness, the postyield stiffness, the equivalent horizontal stiffness, and the yield capacity. The hysteresis loop of the isolator is simulated according to the method recommended by [43]. The procedure for designing LRB isolators is conducted according to [44]. Firstly, according to Regulatory Guide 1.60 in Rev. 1, 1993 [45], the safe-shutdown-earthquake event (peak ground acceleration is equal to 0.3 g), which is the maximum earthquake potential for which the RCC building is designed to sustain and remain functional, is determined, as shown in Figure 6. Due to the reduction demand for the design force (75%~85%) required by a nuclear power company, the isolator frequency is assumed to be 0.49 Hz and the damping ratio of the isolator 15%. It can be also found from Figure 6 that the value corresponding to the frequency and damping ratio is 0.136 g and the reduction of design force reaches 84.5%. Then, the spectral displacement can be computed with the chosen isolation frequency and the damping ratio, as is shown inwhere and are displacement and acceleration spectra, respectively, while is the chosen isolation frequency. Therefore, the calculated is 140.7 mm.

(a)

(b)


According to the base displacement represented by , the stiffness can be determined for each isolator. Next, a total of isolators are applied usingwhere and represent a total stiffness of isolators and a total mass of the structure equivalent to 3.71 × 107 Kg. is the stiffness of each isolator and is the selected isolation frequency. In this study, the LRB isolator having the design horizontal displacement at 100% shear strain of 140 mm and the allowable lateral displacement of 560 mm is selected. The elastic stiffness (), the postyield stiffness () at 100% strain, and the yield capacity () are 7.60 kN/mm, 1.17 kN/mm, and 66 kN, respectively. Then, a total of isolators equal to 251 are obtained and applied. The maximum load () having 2100 kN in the vertical direction for each isolator also satisfies the need of 1448 kN. Since ABAQUS 6.10 software provides the connector element which can be used to model the nonlinear behavior and the damping of isolators, the isolator is simulated using the element in this study. Specifically, the behavior of the isolator for the horizontal direction is modeled using the linear elastic-plastic behavior, as is shown in Figure 7(a). For the vertical direction, the elastic behavior is employed, as is shown in Figure 7(b). The finite element model for the base-isolated RCC building appears in Figure 8. The arrangement diagram of the isolators is displayed in Figure 9.

(a)

(b)


3. Ground Motions
In order to study the effect of sequential ground motions on the RCC with and without isolators, appropriate sequential ground motions must first be defined. Currently, a sequential ground motion can be reproduced mainly using three approaches: the repeated approach [2], the randomized approach [9], and as-recorded sequential ground motions. In this study, as-recorded sequential ground motions including the mainshock and at least one aftershock are employed for carrying out nonlinear time history analyses since those seismic sequences represent the real correlation between the mainshock and the aftershock and would not overestimate the structural response significantly in comparison with artificial seismic sequences [22]. In this study, seismic sequences are selected according to the following criteria: (a) magnitudes of most mainshocks are close to or greater than 6.0 and most major aftershocks are equal to or greater than 5.0; (b) accelerograms are recorded by stations which are located in free field; (c) peak ground acceleration of major aftershocks and mainshocks are equal to or greater than 0.1 g; (d) accelerograms are recorded on rock or stiff soil. Under these criteria, twenty as-recorded mainshock-aftershock seismic sequences, which have been recorded during a short period of time, by the same station, are downloaded from the strong motion database of PEER [40]. Table 1 presents the detailed information of selected mainshock-aftershock seismic sequences. Each of the aforementioned sequential motions has been separated into a single ground motion where a 100s gap is applied between two consecutive seismic events. This gap has zero acceleration ordinates and is sufficient to allow the structure to cease moving after the previous event. In order to represent the seismic hazard and maintain the same correlation between the mainshock and the aftershock, all the sequential ground motions are scaled to make peak ground acceleration (PGA) of mainshocks equal the safe-shutdown-earthquake level of 0.3 g, and the relevant aftershocks are multiplied with the same amplitude factor like the mainshocks. It should be noted that the relationship between the mainshock and the aftershock is beyond the scope of this research and is not considered in this paper.
4. Results and Discussion
4.1. Results for Modal Analysis
The results from the nonlinear dynamic analyses can be analyzed with the aid of modal analysis results. Hence, free vibration analysis is firstly conducted for RCC buildings with and without isolators and the first 20 natural frequencies are shown in Table 2. It shows that for the fixed-base RCC the effective mass participations in and direction are in the order of 77.1% for the first mode and second mode and 20.1% for the 19th mode and 20th mode while they are in the order 100% for the first and second modes for the base-isolated RCC. It means that the dynamic behavior for both types of RCC is mainly controlled by the fundamental mode. Furthermore, it is expected that the isolation system would lead to greater decreases in seismic demands of the RCC on the basis of the target spectra [45] because it reduces the fundamental frequency of the RCC from 4.92 Hz to 0.49 Hz. It is also found from Figure 10 that the fundamental mode for the fixed-base RCC is strongly dependent on lateral movements of the structure, whereas it is associated with pure translation movements of the isolator for the base-isolated RCC.

(a)

(b)
4.2. Acceleration Response
Nonstructural components like facilities located on the structure can be categorized into two parts: acceleration-sensitive and displacement-sensitive. For acceleration-sensitive facilities, their seismic performance must be justified using structural acceleration responses. This part quantifies the benefit of base isolation on acceleration-sensitive facilities under seismic sequence excitations. Maximum top acceleration of the RCCs with and without isolators under mainshocks and seismic sequences appears in Table 3. Table 3 also presents the relative intensity of aftershock ground motions () and the decreasing amplitude ratio (DAR) [34]. The relative intensity of aftershock ground motions is defined as the ratio of PGA for the aftershock to that for the mainshock (). The decreasing amplitude ratio (DAR) is employed to evaluate the effectiveness of the isolation system in the RCC. The expression of DAR is defined as follows:where and are structural responses without and with isolation, respectively.
It can be seen in Table 3 that, in comparison with the fixed-base RCC, average maximum absolute acceleration of the isolated RCC reduces by up to 85.7% under mainshocks and 86.3% under multiple ground motions. This demonstrates that base isolation can result in decreased acceleration demands for structures under both single and sequential earthquakes. It is also found from Table 3 that the average is approximate to 1, indicating that the average intensity of aftershocks is equal to that of mainshocks. Indeed, one of the most concern problems for the structural and nuclear engineers is the seismic performance of the RCC building under an aftershock with the same or even higher intensity than the corresponding mainshock. In Table 3, the average maximum acceleration for the fixed-base RCC increases from 13.25 m/s2 to 14.93 m/s2, rising by 12.7% due to the multiplicity of earthquakes. This increase cannot be neglected for seismic design of facilities. In contrast, the increase of average maximum acceleration for the base-isolated RCC is from 1.89 m/s2 under mainshocks to 2.04 m/s2 under seismic sequences, only rising by 7.9% which is much lower in comparison with the fixed-base RCC. It should also be noted that Coalinga seismic sequence (number (12)) with the equal to 1.81 has induced the largest acceleration increase (11.54 m/s2) for the fixed-base RCC, while the increase for the base-isolated RCC has only 1.93 m/s2. Obtained results indicate that the base isolation can be applied to the RCC to reduce the seismic shaking of the RCC and to facilitate the standardization of facility designs for seismic sequences.
The main characteristics of damaging aftershock sequences for both base-isolated and fixed-base RCCs are further investigated in this part. The damaging aftershock sequences are labeled as those resulting in larger responses (i.e., acceleration or displacement responses) of the RCC due to aftershocks. Figure 11 presents the ratio of the maximum top acceleration for the aftershock (AS) to that for the mainshock (MS) as function of . The figure shows that, for the fixed-base RCC, 9 of the 20 records lead to the ratio of maximum top acceleration more than 1, indicating that 45% of the cases are damaging aftershock sequences. Compared with the fixed-base RCC, the number of records resulting in greater aftershock acceleration responses is 5 for the base-isolated RCC, only accounting for 25% of the cases. It is also noted that when the is more than 0.89, the aftershocks would result in greater acceleration responses for the fixed-base RCC. By contrast, the base-isolated RCC would not be exposed to greater dynamic forces due to aftershocks unless the is more than 1.22. It can be concluded that the adverse effect of aftershock sequences would be decreased due to the base isolation applied to the RCC.

The ratio of the mean response spectra for the mainshock-aftershock (MSAS) to those for the mainshock is displayed in Figure 12. Figure 12(a) shows that, for the fixed-base RCC, the ratios of damaging aftershock sequences are greater between the vibration period of 0 s and 0.52 s while being lower between the vibration period of 0.52 s and 6 s in comparison with those of nondamaging sequences. The ratio for damaging aftershock sequences at the first period of the structure (0.2 s) is 1.48, reaching 1.31 times as large as that of nondamaging sequences. It reveals that the damaging aftershock sequence for the fixed-base RCC is more associated with greater aftershock response spectral values around the structural period (0.2 s) of the RCC. Similar results can be obtained for the base-isolated RCC, as shown in Figure 12(b). Additionally, it is noteworthy from Figure 12(b) that the damaging aftershock sequence for the base-isolated RCC has no direct relevance with the aftershock response spectral values around the isolation frequency (0.49 Hz).

(a)

(b)
4.3. Displacement Response
The displacement demands for both superstructures and isolators are analytically studied in this section. The displacement demands are determined as maximum relative displacements.
4.3.1. Superstructure
Table 4 shows the maximum relative displacements of the superstructure between the top point and the base under single and multiple ground motions. As can be seen from Table 4, compared with the fixed-base RCC, average maximum relative displacements for the isolated RCC reduce by up to 84.9% and 86.6%, respectively, under mainshocks and sequential ground motions. Though the average is close to 1, the average maximum relative displacement for the fixed-base RCC rises by 21.4%, growing from 0.0210 m under mainshocks to 0.0255 m under seismic sequences, while no difference of the average maximum relative displacement for the base-isolated RCC can be found between mainshocks and multiple ground motions. It should be mentioned that Coalinga seismic sequence (number (12)) induces the largest displacement increase (0.015 m) for the fixed-base RCC in comparison with the corresponding mainshock while the increase is only 0.003 m for the base-isolated RCC. Obtained results indicate that isolators can not only reduce the structural deformation significantly but also prevent strong aftershocks from causing more permanent deformation to the structure.
In order to further research the main characteristics of damaging aftershock sequences corresponding to maximum relative displacement for both base-isolated and fixed-base RCCs, the ratio of the maximum relative displacement for the aftershock to that for the mainshock only is displayed in Figure 13 as function of . In Figure 13, for the fixed-base RCC, 10 records result in greater aftershock displacement responses, making up 50% of the cases. For the base-isolated RCC, the percentage of damaging aftershock sequences decreases to 20%. Additionally, the RCC may suffer additional displacements from aftershocks when is more than 0.94 for the fixed-base RCC while being 1.22 for the RCC with isolators. It is demonstrated that the base isolation plays an important part in reducing the adverse aftershock effect on the RCC.

Figure 14 shows the ratio of the mean response spectra for the mainshock-aftershock to those for the mainshock. It can be seen from Figure 14 that damaging aftershock sequences, which refer to greater maximum relative displacements due to aftershocks in comparison with mainshocks, are more dependent on greater aftershock response spectral values around the structural period (0.2 s) of the RCC. The similar conclusion is observed for the maximum top acceleration response as mentioned in Section 4.2. Further observing Figure 14(b), it is found that while the ratio at 0.2 s for damaging aftershock sequences is 1.71 times greater than that for nondamaging sequences, the ratios are much lower in the medium- and long-period region (between the vibration period of 1.35 s and 6 s), with the lowest values being equal to 1 around the isolation period. It is indicated that the spectral values of the aftershock around the isolation period are lower than those of the corresponding mainshock and therefore the responses controlled by the isolation mode during aftershocks may be lower than mainshocks. As a consequence, the damaging aftershock effect dominated by the structural mode can be more compromised as a result of lower responses contributed by the isolation mode during aftershocks, which further demonstrates the advantage of applying the isolation system to the RCC for seismic sequences.

(a)

(b)
4.3.2. Isolators
It should be noted that, for the base-isolated RCC, the decrease of relative displacements of the superstructure is at the expense of the increase of isolator displacements. Maximum displacements of isolators under mainshocks and seismic sequences appear in Table 5. From Table 5, the average maximum isolator displacement under mainshocks is 0.042 m. By contrast, the value for seismic sequences increases to 0.057 m, rising by 35.7% in comparison with mainshocks. It reveals that while the average has only 1.033, indicating that the average intensity of the aftershocks is just that of the mainshocks, mainshock-aftershock effects must not be neglected for the design of isolators as larger deformations of isolators may increase the failure risk of pipelines between the isolated and fixed-base structures.
It is important to study the main characteristics of damaging aftershock sequences for the isolator displacement response. The ratio of the maximum isolator displacement for the mainshock-aftershock to that for the mainshock only is displayed in Figure 15 as function of . In Figure 15, it can be found that greater isolator displacements due to aftershocks are caused by 9 of the 20 records for the value between 0.667 and 1.845. The highest ratio of 2.46 occurs at the value equal to 0.934. The increase of isolator displacements due to aftershocks is disproportionate to the relative intensity of aftershock ground motions (), which means that larger isolator displacements are associated with more factors such as different frequency and spectral shape between mainshocks and seismic sequences.

Figure 16 presents the ratio of the mean response spectra for the mainshock-aftershock to those for the mainshock. Clearly, the ratios for damaging aftershock sequences are greater than those for nondamaging sequences over the whole vibration period, particularly in the medium- and long-period region (between 1 s and 6 s). It indicates that the medium- and long-period structures, like the base-isolated RCC, may suffer greater additional global displacements (or isolator displacements) from the aftershocks with the larger spectral values in the medium- and long-period region in comparison with the mainshocks.

5. Conclusions
In this paper, in order to research the effect of sequential ground motions on the base-isolated RCC building, a suitable finite element model for the base-isolated RCC is set up. For comparison, a fixed-base RCC is also built. A contrast of the analysis results including acceleration and displacements indicates that base isolation technique is an effective and efficient way to reduce the seismic demand of the RCC subject to sequential ground motions. The summarized findings are as follows:(a)Obtained acceleration results show that the average maximum top acceleration of the base-isolated RCC can reduce by 85.7% under mainshocks and 86.3% under sequential ground motions. It is indicated that base isolation system can be applied to the RCC to reduce the seismic shaking of the RCC and facilitate the standardization of facility designs for seismic sequences. The results also show that the damaging aftershock effect on acceleration responses of the RCC would be reduced if the base isolation is applied to the RCC.(b)In comparison with the fixed-base RCC, average maximum top displacements for the base-isolated RCC reduce by 84.9% under mainshocks and 86.6% under seismic sequences. By setting up the isolation system, no difference for average maximum top displacements of the superstructure can be found between mainshocks and seismic sequences. The results reveal that isolators can not only reduce the structural deformation significantly under earthquake excitations but also prevent strong aftershocks from causing additional deformation to the structure. Additionally, base isolation can reduce the damaging aftershock effect on displacement responses of the RCC significantly.(c)It is observed that the average maximum isolator displacement for seismic sequences rises by 35.7% in comparison with the mainshocks, which indicates that mainshock-aftershock effects must not be neglected for the design of isolators as larger deformations of isolators may increase the failure risk of pipelines between the isolated and fixed-base structures. The medium- and long-period structures, like the base-isolated RCC, may suffer greater additional global displacement (or isolator displacement) from the aftershocks with larger spectral values in the medium- and long-period region in comparison with the mainshocks.(d)Ratio of response spectral statistics reveals that damaging aftershock sequences for both the fixed-base RCC and the superstructure of the RCC with isolators are more associated with greater aftershock response spectral values around the structural period (0.2 s) of the RCC. Notice that damaging aftershock sequences for the superstructure of the base-isolated RCC have no direct relevance with the aftershock spectral values around the isolation period (2.04 s).
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This study received support from the China Scholarship Council.