Next Article in Journal
A Framework for Prefabricated Component Hoisting Management Systems Based on Digital Twin Technology
Next Article in Special Issue
Vibration-Based and Near Real-Time Seismic Damage Assessment Adaptive to Building Knowledge Level
Previous Article in Journal
Study on Breaking Concrete Structures by Pulse Power Technology
Previous Article in Special Issue
Assessment of the Seismic Behavior of a Precast Reinforced Concrete Industrial Building with the Presence of Horizontal Cladding Panels
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Peak and Cumulative Response of Reinforced Concrete Frames with Steel Damper Columns under Seismic Sequences

Department of Architecture, Faculty of Creative Engineering, Chiba Institute of Technology, Chiba 275-0016, Japan
Buildings 2022, 12(3), 275; https://doi.org/10.3390/buildings12030275
Submission received: 21 January 2022 / Revised: 15 February 2022 / Accepted: 25 February 2022 / Published: 28 February 2022
(This article belongs to the Collection Structural Analysis for Earthquake-Resistant Design of Buildings)

Abstract

:
The steel damper column is an energy-dissipating member that is suitable for reinforced concrete (RC) buildings and those used for multistory housing in particular. However, the effectiveness of steel damper columns may be affected by the behavior of surrounding members, and this effect can be severe in the case of seismic sequences. This article investigates the nonlinear response of building models with an RC moment-resisting frame (MRF) with and without steel damper columns under seismic sequences. The applicability of the concept of the momentary energy input to the prediction of the peak response of RC MRFs with damper columns under seismic sequences is also investigated. The main findings of the study are summarized as follows. (1) The peak response of RC MRFs with damper columns subjected to sequential accelerations is similar to the peak response obtained considering only the mainshock, whereas the cumulative strain energy of RC MRFs accumulates more for sequential accelerations. (2) The steel damper column is effective in reducing the peak and cumulative responses of RC MRFs in the case of sequential seismic input. (3) The relation of the hysteretic dissipated energy during a half cycle of the structural response and the peak displacement of the first modal response can be properly evaluated using the simple model proposed in this study.

1. Introduction

1.1. Background

In an earthquake-prone country such as Japan, controlling seismic damage to a structure is an important issue in the seismic design of building structures. A popular and classical strategy for improving the damage control ability of the moment-resisting frame (MRF) is the so-called weak-beam strong-column concept. This strategy is widely accepted and is recommended for the seismic design of MRFs. When this strategy is adopted, most of the seismic energy is absorbed by plastic hinges set at each beam end. However, because the beams also carry gravitational loads, those MRF buildings may not continue to be usable after a huge earthquake due to the severe damage to their beams. In addition, as was the case in the 2016 Kumamoto earthquake, there may be a sequence of large foreshocks and the mainshock or a sequence of the mainshock and large aftershocks. In such cases, seismic energy accumulates at the plastic hinges and causes damage. Therefore, an MRF designed solely according to the weak-beam strong-column concept may be insufficient in the case of such seismic sequences. A dual system with sacrificial members that absorb the seismic energy prior to the beams and columns (e.g., in a damage-tolerant structure) [1] is one solution for creating structures with superior seismic performance.
The steel damper column [2] is an energy-dissipating sacrificial member. Figure 1 compares the design collapse mechanism of the traditional MRF and MRF with steel damper columns. In the traditional MRF, shown in Figure 1a, most of the seismic energy is absorbed at the plastic hinges of the beam ends and the bottom end of the first story columns. Meanwhile, for the MRF with steel damper columns shown in Figure 1b, the damper panel within a damper column absorbs seismic energy prior to the plastic hinges in the beams and columns. The energy absorbed by the plastic hinges can thus be reduced using steel damper columns.
Figure 2 shows the design concept of a reinforced concrete (RC) MRF with steel damper columns. Unlike the use of buckling-restrained braces, the use of steel damper columns provides usable space for corridors, as shown in Figure 2a. The steel damper column is thus suitable for high-rise RC housing in that it minimizes obstacles in architectural planning. Figure 2b presents an example of the RC beam–steel damper column joint. The steel beam embedded in the RC beam transfers the bending moment from the RC beam to the damper column.
The author previously investigated the development of the rational seismic design procedure of an RC MRF with steel damper columns [3,4,5]. Specifically, a displacement-controlled seismic design method was proposed for an RC MRF with steel damper columns [3,4]. With this method, the strength demand of the whole structure is determined through equivalent linearization. Although the method is simple, it does not consider the accumulated strain energy of members (both RC members and steel damper columns) or the effect of seismic sequences. Another important issue with an RC MRF with steel damper columns is that the effect of the steel damper columns depends on the behavior of the surrounding RC beams. The author previously found [5] that a proper strength balance of the steel damper columns and surrounding beams is important in maximizing the dissipation of energy into the damper columns; i.e., the beam-end section connected to the damper columns requires sufficient strength to avoid premature yielding prior to any energy dissipation.
The previous results reveal the importance of studying the nonlinear seismic behavior of an RC MRF with steel damper columns designed using the previously proposed method [3] under seismic sequences and the effectiveness of steel damper columns for the reduction of structural damage.

1.2. Brief Review of Related Studies

1.2.1. Studies on the Responses of Structures under Seismic Sequences

Many studies have investigated the response of structures under seismic sequences (e.g., [6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46,47,48,49,50]). To the best of the knowledge of the author, the first study on the nonlinear response of structures to seismic sequences was conducted by Mahin [6]. Later, Amadio et al. investigated the nonlinear response of the idealized single-degree-of-freedom (SDOF) model and two-dimensional frame model under repeated earthquake ground motions [7]. In their study, identical ground motions were applied several times as the seismic sequence. Hatzigeorgiou studied the nonlinear response of an SDOF model [8,9,10] and RC frame models [11,12,13]. In those studies, the artificial seismic sequences were created by selecting ground motions at random, neglecting the difference in the frequency characteristics of the mainshock and aftershocks (or foreshocks). Ruiz-García et al. pointed out the problem with the assumptions made in creating artificial seismic sequences [14,15,16,17,18,19,20,21,22]. Specifically, they showed that the predominant period of aftershocks is shorter than that of the mainshock because the magnitudes of the aftershocks are smaller than the magnitude of the mainshock (e.g., [14,15,17]). They thus concluded that using artificial seismic sequences with repeating identical accelerations may lead to the overestimation of the effects of seismic sequences [15]. In addition, they pointed out that the ratio of the predominant periods of aftershocks and the mainshock is an important parameter in discussing the effect of a seismic sequence [14,16,18,20]. They also conducted three-dimensional frame analyses considering seismic sequences [19,22]. Similar findings were obtained by Goda et al. [23,24,25,26], who proposed a method of generating artificial ground motion sequences by considering the difference in magnitude between the mainshock and aftershocks [23]. They also pointed out the importance of the record selection of aftershocks in incremental dynamic analysis [25]. Tesfamariam et al. studied the seismic vulnerability of RC frames with unreinforced masonry infill due to mainshock–aftershock sequences [27]. Tesfamariam and Goda proposed the seismic performance evaluation framework considering maximum and residual inter-story drift ratio of on-code conforming RC buildings [28] and also energy-based seismic evaluation method of tall RC buildings [29]. In [28], the authors concluded that although the influence of the mainshock–aftershock sequences to the maximum inter-story drift is limited, the seismic performance of non-code-conforming RC buildings are influenced by the mainshock–aftershock sequences because the residual inter-story drift of RC buildings increases due to the seismic sequences. Similarly, they concluded that the influence of major aftershocks on the damage potential is significant because the energy-based damage index of tall RC buildings increases due to the seismic sequences. In addition to these studies, Zhai et al. [30,31,32], Di Sarno et al. [33,34,35,36], Abdelnaby et al. [37,38,39,40,41], Kagermanov and Gee [42,43], Yang et al. [44,45], Yaghmaei-Sabegh et al. [46], Qiao et al. [47], Orlacchio et al. [48], Hoveidae and Radpour [49], and Pirooz et al. [50] investigated the nonlinear response of building structures under seismic sequences. Although most of those studies were analytical in nature, one experimental study used a shaking table [47].
From the author’s viewpoint, few studies have investigated the application of the pushover-based procedure using a nonlinear equivalent SDOF model for the seismic evaluation of a building considering seismic sequences; the existing studies are those of Guerrero et al. [21], Kagermanov and Gee [42,43], and Orlacchio et al. [48]. Such studies are essential in the author’s opinion. One reason is that the analysis of the nonlinear response of structures using the equivalent SDOF model can lead to a better understanding of the nonlinear behavior of structures. As an example, relations between the seismic intensity parameter and response quantities (e.g., the peak response and cumulative energy) can be clearly discussed using an equivalent SDOF model. In addition, from a practical point of view, such pushover-based procedures provide structural designers and analysts with basic information on the nonlinear behavior of the analyzed building.
Another issue to be addressed is that there have been few studies on structures with dampers, with the existing studies being those of Guerrero et al. [21], Yang et al. [45], and Hoveidae [49]. As described above, the main motivation for installing such dampers is the energy absorption of the sacrificial members prior to beams and columns. It is therefore essential to discuss the effect of dampers in terms of the cumulative strain energy in the event of seismic sequences. However, few discussions have been presented in the studies cited above.

1.2.2. Studies on the Seismic Energy Input

The quantification of structural damage to members, such as RC beams, columns, and dampers, is an important issue in rational seismic design and evaluation. Several indices of structural damage have been proposed, e.g., the Park–Ang index, which is defined as the combination of the peak deformation and cumulative energy [51]. Because the cumulative energy is directly evaluated from the seismic energy input, it is rational to consider the seismic intensity according to energy-related parameters. The total input energy [52,53] is a seismic intensity parameter related to the cumulative strain energy. Several studies have investigated total input energy spectra (e.g., [54,55,56,57,58]).
Inoue and his team proposed the maximum momentary input energy [59,60,61] as an energy-related seismic intensity parameter related to the nonlinear peak displacement. They predicted the peak displacement by equating the maximum momentary input energy and hysteretic dissipated energy in a half cycle of the structural response. The definition of the momentary input energy is described in Appendix A.
Following their work, the present author formulated the time-varying function of the energy input using a Fourier series [62]. This formulation shows that two seismic intensity parameters, namely the maximum momentary input energy and total input energy, can be evaluated based on the properties of the system and the complex Fourier coefficient of the ground motion. The concept of the maximum momentary input energy has been extended for bidirectional excitation [63] and implemented in the prediction of the peak and cumulative responses of a one-mass, two-degree-of-freedom model representing a ductile RC structure [64] and an irregular base-isolated building subjected to bidirectional horizontal ground motions [65]. The application of the momentary energy input to the RC MRF subjected to seismic sequences is thus promising.

1.3. Objectives

Against the above background, the following questions are addressed in this article.
  • What are the differences in the peak and cumulative responses of RC MRFs with and without steel damper columns between a single acceleration and sequential accelerations?
  • Is the steel damper column effective in reducing the peak and cumulative responses of an RC MRF in the event of a seismic sequence?
  • In the prediction of the peak response of the RC MRF with steel damper columns based on the momentary energy input, the relation between the hysteretic dissipated energy during a half cycle of the structural response and the peak displacement must be properly modeled. How can this relationship be modeled from the results of pushover analysis?
The present article investigates the nonlinear response of 10-story RC MRF building models with steel damper columns designed according to the previously proposed method [3] in a case study of such an RC MRF with steel damper columns subjected to seismic sequences. Ground-motion records obtained from three stations managed by the National Research Institute for Earth Science and Disaster Resilience (NIED) during the foreshock and mainshock of the 2016 Kumamoto earthquake [66] are used. The effects of the seismic sequences on the peak and cumulative responses of the RC MRFs are investigated using the results of nonlinear time-history analysis. The first modal response of RC MRFs is then calculated from the results of the nonlinear time-history analysis, and the applicability of the momentary input energy to the prediction of the peak response under a seismic sequence is discussed.
The remainder of the article is organized as follows. Section 2 presents four RC MRF building models with and without damper columns as well as the ground motion data used in the nonlinear time-history analysis. The results of nonlinear time-history analysis are presented and discussed in Section 3. This section further analyzes the peak and cumulative responses of the RC MRF and the effectiveness of the steel damper column. The first modal response is evaluated using the results of nonlinear time-history analysis and discussed in Section 4. Conclusions and future directions of study are discussed in Section 5.

2. Building and Ground Motion Data

2.1. Building Data

Figure 3 shows the simplified structural plan and elevation of the MRF building models with steel damper columns considered in this study. Two 10-story building models, namely RCDC1 and RCDC2, are designed using the simplified design method [3,4]. Specifically, RCDC1 is the building model presented in the previous study [3], and RCDC2 is the building model obtained by shortening the span of the RC beam, while the change in the total strength of the overall building model (comprising both the RC frame and damper column) is minimized [4]. The model RCDC1 represents an MRF that has been intentionally designed to be flexible, whereas the model RCDC2 represents an ordinary MRF. The unit mass per floor is assumed to be 1.2 t/m2. The story height is assumed to be 4.5 m for the first story and 3.2 m for upper stories. Details of the members are given in the Appendix B.
To investigate the effect of the steel damper column in reducing the seismic response, two other models are considered by removing all steel damper columns from RCDC1 and RCDC2. Here, the models with dampers removed from RCDC1 and RCDC2 are, respectively, referred to as BareRC1 and BareRC2. Four MRF building models are thus considered in this study.
Figure 4 shows the modeling of the MRF with steel damper columns. In this study, all RC MRFs are designed according to the weak-beam, strong-column concept. A potential hinge is set at all RC beam ends (except the beam end connected to the steel damper columns) and the bottom end of the first-story columns as shown in Figure 4a.
The building is modeled as having a planar frame, as shown in Figure 4b. All RC members are modeled as an elastic beam with a nonlinear flexural spring at both ends. The steel damper columns are modeled as an elastic column with a nonlinear damper panel at the middle of the steel damper column. The beam-column joints are assumed to have rigid behavior. For RC beam–RC column joint, the proper reinforcement is assumed to be provided to prevent premature failure. In addition, the proper dimensions (section sizes and lengths) of steel beam embedded in RC beams and proper reinforcement is assumed to be provided in RC beam–damper column joint, to prevent premature failure until the damper panels reach their ultimate stage.
Figure 5 shows the envelope of the force–deformation relationship of each member. The modeling of each member is the same as that in the previous study [3] and is summarized as follows. The envelopes are assumed to be symmetric under positive and negative loading. At the potential hinge of RC members, the crack and yielding of the section are considered as shown in Figure 5a. The yield moment M y (at point Y) is calculated according to the AIJ standard [67], whereas the crack moment M c (at point C) is assumed to be one-third of M y . The secant stiffness degradation ratio at point Y ( α y ) is calculated using the equation proposed by Sugano and Koreishi [68]. Meanwhile, at the non-hinge ends of RC members, only the cracking of the section is considered, as shown in Figure 5b, except for the beam at the foundation level. The cracking moment ( M c ) of the non-hinge beam end is assumed to be the same as that at the opposite end, and the tangent stiffness degradation ratio after cracking ( α 1 ) is assumed to be the same as the secant stiffness degradation ratio at yielding ( α y ) calculated following Sugano and Koreishi [68]. Meanwhile, at the non-hinge column end, the cracking moment ( M c ) is calculated considering the axial force attributed to the vertical load, and the tangent stiffness degradation ratio after cracking ( α 1 ) is assumed to be 0.2. The flexural behavior of a beam at the foundation level is assumed to be linearly elastic. The shear behavior is assumed to be linearly elastic of all RC members. In nonlinear static analysis, the bilinear envelope shown in Figure 5c is assumed for the damper panel. Here, Q y D L and Q y D U respectively denote the initial and upper bound yield strengths of the damper panel. The axial behavior is assumed to be linearly elastic for all vertical members.
Figure 6 shows the hysteresis rule of the nonlinear spring. In this study, the Muto model [69] with two modifications is used for the flexural spring in RC members. The first modification is the unloading of the stiffness after yielding to represent the degradation of RC members after yielding, as shown in Figure 6a. The second modification is the consideration of stiffness degradation after yielding due to cyclic loading. In this study, the model proposed by Umemura et al. [70] is implemented with the Muto model, as shown in Figure 6a. A parameter χ is introduced to represent the effect of cyclic degradation. The value of χ can be taken as zero or positive. In the case that χ is zero, there is no cyclic degradation, as shown in Figure 6a, which corresponds to the model used in the previous study (e.g., [3,5,64]). In the case that χ is positive, there is stiffness degradation due to the reloading target point shifting from point Pp to Pn, as shown in Figure 6a. Umemura et al. [70] showed that χ depends on (i) the compressive strength of the concrete, (ii) the shear reinforcement ratio, (iii) the compressive stress of the section normalized by the compressive strength of the concrete, and (iv) the shear-span-to-depth ratio of the member. In this study, the value of χ is set as given in Table 1 following Umemura et al. [70]. Note that a larger value of χ is set for the short-span beam; i.e., the stiffness degradation of the short-span beam is more severe than that of the long-span beam. Additionally, it is noted that the cyclic stiffness degradation is more severe in RCDC2 than in RCDC1.
It should be mentioned that the pinching behavior observed in RC members is not considered in this study. There are several hysteresis models that implement the pinching behavior, e.g., the model proposed by Baker and Noori [71]. Such pinching behavior may affect the response of damper columns. However, this issue would be the next phase of this study.
The model proposed by Ono and Kaneko [72] shown in Figure 6b is used to model the hysteresis behavior of the damper panel. In this model, the strain-hardening effect of the low-yield-strength steel is controlled by the two parameters α and β . Here, the two parameters are set as α = 0.022 and β = 0.013 . Note that the model has been verified by comparing the test results of shear damper panel and braces using low-yield-strength steel in the literature [72], and this model is implemented in several commercial computer programs used for the actual structural design in Japan. It should be mentioned that Vaiana et al. proposed a more generic hysteretic model for rate-independent mechanical systems [73]. They have applied their models for the modeling of rate-independent passive energy-dissipation devices [74]: the calibration of their model has been made by comparisons of the experimental tests of steel shear link tested by Nuzzo et al. [75]. Therefore, such a model would be more suitable for the hysteresis model of damper panel than the model applied herein. However, this issue is out of scope of this study.
The damping matrix is assumed to be proportional to the instantaneous stiffness matrix without a damper column. The damping ratio of the first elastic mode of the model without a damper column is assumed to be 0.03. Second-order effects, including the P-Δ effect, are neglected in this study. The soil–structure interaction (SSI) effect is neglected for the simplicity of the analysis.
Table 2 gives the natural periods of the first three modes in the initial stage for each model. As shown here, the natural period of RCDC1 is longer than that of RCDC2.
In examining the nonlinear behavior of the four models, pushover analysis (i.e., displacement-based mode-adaptive pushover analysis [76]) is carried out to obtain the relationship between the equivalent acceleration A 1 * and equivalent displacement D 1 * . Let f n R and d n be the restoring force vector and horizontal displacement vector, respectively, of the building model at each loading step n obtained in the pushover analysis. The equivalent displacement and acceleration at step n (namely D n 1 * and A n 1 * ) are determined from Equations (1) and (2), respectively, assuming that the vector d n is proportional to the first mode vector ( Γ n 1 φ n 1 ) at each loading step:
D n 1 * = Γ n 1 φ n 1 T M d n M n 1 * = d n T M d n d n T M 1 ,
A n 1 * = Γ n 1 φ n 1 T f n R M n 1 * = d n T f n R d n T M 1 ,
M n 1 * = Γ n 1 φ n 1 T M 1 = ( d n T M 1 ) 2 d n T M d n
d n = { y n 1 y n N } T ,
f n R = { f n R 1 f n R N } T ,
M = [ m 1 0 0 m N ] ,
1 = { 1 1 } T .
In Equations (1)–(6), M is the mass matrix, M n 1 * is the effective first modal mass at loading step n , and m j is the floor mass of the j th floor. Figure 7 shows the obtained A n 1 * D n 1 * relationship for each model. In the figure, the point labeled “design target” is the point assumed as the displacement limit D 1 * l i m i t (=0.2833 m), which is assumed to be 1/82.5 of the assumed equivalent height H 1 * (=23.37 m). As shown in Figure 7a,b, the displacement upon the first yielding of the RC member is larger for BareRC1 than for BareRC2. A similar observation is made in the comparison of Figure 7c,d, whereas the displacement upon the first yielding of the damper column is smaller for RCDC1 than for RCDC2.

2.2. Ground Motion Data

The present study uses records of accelerations of the foreshock (event time: 14 April 2016, 21:26 (JST), JMA magnitude 6.5) and mainshock (event time: 16 April 2016, 01:25 (JST), JMA magnitude 7.3) obtained at three stations managed by NIED, namely K-NET Kumamoto (KMM), K-NET Uto (UTO), and KIK-NET Mashiki (MAS). Table 3 lists the ground motions, and Figure 8 presents the primary and shear wave profiles for the different stations. The soil properties for each station are available from the K-NET website [66].
Figure 9 shows the recorded accelerations observed at the three stations. The present study uses the first 60 s of the as-recorded acceleration records shown in the figure for nonlinear time-history analysis.
Figure 10 shows the elastic pseudo-velocity spectra ( S p V ) of the accelerations. The viscous damping for the calculation of S p V is set at 0.05. The “design” earthquake shown in this figure is the code-specified spectrum (soil condition: type-2) of the Building Standard Law of Japan [77] defined as
S p V ( T , 0.05 ) = T 2 π S p A ( T , 0.05 ) ,
S p A ( T , 0.05 ) = { 4.8 + 45 T m / s 2 T 0.16   s 12.0 0.16   s < T 0.864   s 12.0 ( 0.864 / T ) T > 0.864   s .
Note that the design earthquake spectrum is used for the design of models RCDC1 and RCDC2. The calculated spectra of the mainshock are close to the design earthquake spectrum at KMM (Figure 10a). In contrast, the calculated spectrum of the mainshock exceeds the design earthquake spectrum at MAS (Figure 10c). Note that the spectrum of the east–west (EW) component of the foreshock at MAS also exceeds the design earthquake spectrum, which implies that the EW components of both the foreshock and mainshock at MAS are more intense than the design earthquake.
Figure 10 also shows the predominant periods of the foreshock and mainshock, denoted T g F and T g M , respectively. The predominant period is defined as the peak period of S p V , following Miranda [78] and Ruiz-García [14,17]. The figure reveals that the relation of the two predominant periods ( T g F and T g M ) depends on the site and component. In the cases of the EW-component at KMM and both components at MAS, T g M is longer than T g F , whereas T g M is shorter than T g F for the north–south (NS) component at KMM and the EW-component at UTO.
Table 4 lists the cases of ground motion considered in this study. Here, Cases F and M are, respectively, the single acceleration of only the foreshock and that of the mainshock, whereas cases FM and MF are sequential accelerations, with Case FM following the recorded order of first the foreshock (e.g., KMM0414EW) and second the mainshock (e.g., KMM0416EW) and Case MF following the opposite sequence of first the mainshock and second the foreshock. A time interval of 30 s is set between the first and second accelerations.

3. Analysis Results

3.1. Peak Response

3.1.1. Relative Displacement and Story Drift

Figure 11 compares the distribution of the peak relative displacement in the four cases for each model. The input ground motion set is the EW component at KMM (KMM-EW). The figure shows that the peak response in Case FM is larger than that in Case F, whereas the peak response in Case MF is the same as that in Case M for all models. Specifically, as shown in Figure 11a, the peak response of BareRC1 obtained in Case FM is larger than that obtained in Cases MF and M. This implies that, in Case FM, the damage to members due to the foreshock affects the response of BareRC1 during the mainshock (following the foreshock).
For the MRFs with steel damper columns, the peak response of RCDC1 (Figure 11c) in Case FM is within the design target, whereas it slightly exceeds the design target in Case MF. Similar results are obtained for RCDC2 as shown in Figure 11d. Figure 10a shows that the pseudo-velocity spectra of the foreshock and mainshock obtained from the EW components at KMM are close to those of the design earthquake. Figure 11c,d thus show that the strength demands of those two models are properly determined using the simplified design method proposed in the previous study [3].
Figure 12 compares the distribution of the peak relative displacement in the case that the input ground is set as MAS-EW. The figure shows that the peak displacement is larger than that in the case of KMM-EW shown in Figure 11. Specifically, the peak response of BareRC1 in Case FM is the same as that in Case F, whereas the peak responses of the other models in Case FM are larger than those in Case F. Additionally, the peak responses of RCDC1 and RCDC2 exceed the design target in all four cases, as shown in Figure 12c,d. This result is consistent with the observation from Figure 10c that the pseudo-velocity spectra obtained from the EW components at MAS are larger than those of the design earthquake in a wide range of the period.
Figure 13 and Figure 14 compare the peak story drift in the four cases for each model. The input ground motion sets are KMM-EW (Figure 13) and MAS-EW (Figure 14). The figures confirm the results in Figure 11 and Figure 12. Note that the largest peak story drifts of RCDC1 and RCDC2 are close to 1/75 in the case of KMM-EW. The value of 1/75 is the assumed target drift in the design of RCDC1 and RCDC2, and the responses of both models under the sequential accelerations of KMM-EW (both Cases FM and MF) are close to the assumed design limit. In contrast, the responses of both RCDC1 and RCDC2 greatly exceed the assumed design limit in the case of MAS-EW.
Figure 15 compares the peak story drifts in the cases of single and sequential accelerations to clarify the effect of sequential accelerations on the peak story drift. It is seen that there is a notable difference in the peak drift between Cases F and FM. The difference in the peak drift between Cases F and FM is more pronounced for the MRFs with steel damper columns (RCDC1 and RCDC2) than for the MRFs without dampers (BareRC1 and BareRC2). In contrast, the difference in the peak drift between Cases M and MF is negligibly small for all models. These observations imply that the peak responses of all models under sequential accelerations studied here are governed by the mainshock. This point is discussed further in Section 4.

3.1.2. Member Deformation

Next, the peak responses of the member deformation (plastic rotation of the RC beam end and shear strain of the damper panel) are investigated. Figure 16 compares the peak rotation of the plastic hinge at the beam end. MAS-EW is used as the input ground motion set because it gives the largest response among all input ground motion sets. The figure shows the plastic hinge rotation at beam ends in the inner frame (i.e., the right-side ends of beam X2X3 in Frame Y3 for BareRC1 and RCDC1 and the left-side ends of beam X3X4 in Frame Y3 for BareRC2 and RCDC2). It is seen that the peak plastic rotation in Case FM is larger than that in Case F except for the model BareRC1, whereas the plastic rotation in Case MF is the same as that in Case M. This result is consistent with Figure 14. Figure 16 also shows that the plastic rotation of beams at upper floor levels (higher than level Z7) is zero; there is no yielding of these beams.
Figure 17 compares the peak plastic hinge rotation at the beam end in the cases of single and sequential acceleration. There is a notable difference in the peak plastic rotation between Cases F and FM, whereas the difference in the peak drift between Cases M and MF is negligibly small for all models. This result is consistent with Figure 15.
Figure 18 compares the peak shear strain of the damper panel in RCDC1 and RCDC2, where the input ground motion set is MAS-EW. It is seen that the peak shear strain of the damper panel is greater for the lower stories. In addition, the peak shear strain in Case FM is greater than that in Case F, whereas the peak shear strain in Case MF is similar to that in Case M.
Figure 19 compares the peak shear strain of the damper panel in cases of single and sequential acceleration. There is a notable difference in the peak shear strain between Cases F and FM, whereas the difference in the peak drift between Cases M and MF is negligibly small for both RCDC1 and RCDC2.

3.2. Cumulative Response

This subsection discusses the cumulative responses of the four building models. First, the cumulative responses of the overall building (i.e., cumulative input energy E I , cumulative strain of the RC MRF and steel damper columns, E S f and E S d , respectively, and cumulative viscous damping energy E D ) are discussed. The cumulative strain energy of the member (plastic hinge of the RC beam and damper panel) is then discussed.

3.2.1. Cumulative Response of the Overall Building

Figure 20 shows the cumulative strain energy of the RC MRF and steel damper columns and the cumulative viscous damping energy per unit mass. Note that the total input energy E I is expressed as the sum of E S f , E S d , and E D .
The following observations are made for Figure 20.
  • The total input energy E I of the sequential accelerations is greater than that of the single acceleration: e.g., E I in Case FM is greater than that in Case F.
  • In most cases for BareRC1 and BareRC2, E I is mostly absorbed as cumulative strain energy of the RC MRF, E S f .
  • In most cases for RCDC1 and RCDC2, the total input energy E I is greater than that of BareRC1 and BareRC2. However, a large proportion of E I is absorbed as the cumulative strain energy of the steel damper columns, E S d . The relative amounts of E S f and E S d depend on the model and analysis case.

3.2.2. Cumulative Strain Energy of a Member

First, the cumulative strain energy of the plastic hinges at the RC beam end is investigated. The normalized strain energy of the k th plastic hinge, N E S f k , is defined as
N E S f k = E S f k M y k θ y k = 1 M y k θ y k 0 t d M k ( t ) θ ˙ k ( t ) d t ,
where M k ( t ) and θ k ( t ) are, respectively, the moment and rotation of the k th plastic hinge at time t , and t d is the time length of the nonlinear time-history analysis.
Figure 21 compares N E S f k in the cases of single and sequential accelerations for BareRC1 and BareRC2. In both models, N E S f k at beam ends in the inner frame (the same beam ends considered in Section 3.1.2) is investigated. The figure shows that N E H f k is notably greater in Case FM than in Case F. In addition, the increase in N E S f k from Case M to Case MF is not negligible. This implies that the effect of sequential accelerations on the cumulative strain energy is more pronounced than that on the peak response. The figure also shows that, in general, N E S f k is larger for BareRC2 than for BareRC1. This is because the yield deformation angle of beams ( θ y k ) is smaller for BareRC2 than for BareRC1.
Figure 22 compares N E S f k in the cases of single and sequential accelerations for RCDC1 and RCDC2. For RCDC1 and RCDC2, N E S f k at the long-span beam ends (the same beam ends considered in Section 3.1.2) and short-span beam ends (the right-side ends of beam X2X3 in Frame Y2 for RCDC1 and the left-side ends of beam X3X4 in Frame Y2 for RCDC2) is investigated.
The figure shows a notable increase in N E S f k due to sequential accelerations for RCDC1 and RCDC2. Specifically, for RCDC1, the increase in N E S f k in Frame Y2 is larger than that in Frame Y3. The reason is that (i) the yield deformation angle of beams ( θ y k ) in Frame Y2 is smaller than that in Frame Y3 and (ii) the cyclic stiffness degradation is more pronounced for the short-span beam (Frame Y2) than for the longer span beam (Frame Y3).
Next, the cumulative strain energy of the damper panel is investigated. The normalized strain energy of the k th damper panel, N E S d k , is defined as
N E S d k = E S d k Q y D L k γ y D L k h d 0 k = 1 Q y D L k γ y D L k 0 t d Q D k ( t ) γ ˙ D k ( t ) d t ,
where Q D k ( t ) and γ D k ( t ) are, respectively, the shear force and shear deformation angle of the k th damper panel at time t .
Figure 23 compares N E S d k in the cases of single and sequential acceleration. It is seen that the increase in N E S d k due to the sequential accelerations is pronounced in both models. In addition, the value of N E S d k for RCDC1 is greater than that for RCDC2.
It is concluded from the results presented in this subsection that the effect of the sequential accelerations on the cumulative strain energy in the member is not negligible. In contrast with the peak deformation, the cumulative strain energy accumulates in the event of sequential accelerations.

3.3. Effectiveness of Steel Damper Columns in Reducing the Seismic Response

In this subsection, the effectiveness of the steel damper column in reducing the peak and cumulative responses of the RC members is discussed by comparing results for MRFs without dampers (BareRC1, BareRC2) and those with steel damper columns (RCDC1, RCDC2).
Figure 24 compares the peak story drift for models without and with steel damper columns. As shown in Figure 24a, in general, the drift for RCDC1 is smaller than that for BareRC1. The effect of the steel damper columns in the reduction of the peak drift is obvious in the case of sequential accelerations as well as the case of single acceleration. A similar observation can be made in comparing RCDC2 and BareRC2.
Figure 25 compares the peak plastic hinge rotation at the beam ends for Frame Y3. It is seen that the peak plastic rotation of the models with steel damper columns (RCDC1, RCDC2) is smaller than that of the models without dampers (BareRC1, BareRC2). It is concluded from the comparisons that the steel damper column is effective for the reduction in the peak response of the RC MRF, as far as the models studied herein are concerned.
The cumulative responses are next compared. Figure 26 compares the cumulative strain energy of the RC MRF per unit mass. As shown in Figure 26a, in most cases, the cumulative strain energy of RCDC1 is less than that of BareRC1 for both single acceleration and sequential acceleration. Similar observations are made for Figure 26b.
Figure 27 compares the normalized cumulative strain energy at the beam ends ( N E S f k ). It is seen that N E S f k with steel damper columns (RCDC1, RCDC2) is smaller than that without dampers (BareRC1, BareRC2). It is concluded from the comparisons that the steel damper column is effective for the reduction of the cumulative response of the RC MRF, as far as the models studied herein are concerned. Note that the total input energy of MRFs with steel damper columns is greater than that of MRFs without dampers in some cases, as shown in Figure 20. Even in such cases, those steel damper columns installed in the RC MRF play an important role of absorbing the seismic energy, as shown in Figure 27.

3.4. Summary of Results

The results presented in this section are summarized as follows.
  • The peak responses of all models under sequential accelerations studied here are governed by the mainshock. The peak responses of all models in Case FM are notably greater than those in Case F. In contrast, the difference in the peak response between Case MF and Case M is limited.
  • The effect of the sequential accelerations on the cumulative strain energy is not negligible. Unlike the peak response, the cumulative strain energy accumulates in the event of sequential accelerations.
  • The steel damper column is effective in reducing the peak and cumulative responses of RC members, irrespective of whether a single acceleration or sequential accelerations is considered as the seismic input.

4. Evaluation of the First Modal Response

4.1. Method of Calculating the First Modal Response from the Results of Time-History Analysis

This section calculates the equivalent displacement D 1 * ( t ) and acceleration A 1 * ( t ) from the results of nonlinear time-history analysis. This study assumes that the building model considered here oscillates predominantly in the first mode, and the first mode vector at the peak response ( Γ p e a k 1 φ p e a k 1 ) is then assumed from the relative horizontal displacement vector ( d ( t ) ). Note that in other studies (e.g., [65]), the first mode vector is assumed by referring to pushover analysis considering the change in the first mode vector at each nonlinear stage (displacement-based mode-adaptive pushover analysis [76]). However, because the present study considers the cyclic stiffness degradation of RC members and the strain-hardening effect of damper panels, the assumption of the first mode vector at the peak response from the pushover analysis results cannot be applied. Therefore, in this study, the assumption of Γ p e a k 1 φ p e a k 1 is based on the results of nonlinear time-history analysis. The calculation method is described below.

4.1.1. Step 1: Determination of the “Peak Response Point”

The displacement at “the center of the mass” D * ( t ) is calculated as
D * ( t ) = 1 T M d ( t ) 1 T M 1 .
The time t p e a k at which the absolute value of D * ( t ) is a maximum is then found.

4.1.2. Step 2: Determination of the First Mode Vector at the Peak Response Point

The first mode vector at the peak response point ( Γ p e a k 1 φ p e a k 1 ) is assumed to be proportional to the relative horizontal displacement vector at time t p e a k ( d ( t p e a k ) ):
Γ p e a k 1 φ p e a k 1 d ( t p e a k ) T M 1 d ( t p e a k ) T M d ( t p e a k ) d ( t p e a k ) .
The effective first modal mass ( M p e a k 1 * ) is then calculated as
M p e a k 1 * = Γ p e a k 1 φ p e a k 1 T M 1 .

4.1.3. Step 3: Calculation of the Equivalent Displacement and Acceleration of the First Modal Response

The equivalent displacement D 1 * ( t ) and acceleration A 1 * ( t ) are calculated as
D 1 * ( t ) = Γ p e a k 1 φ p e a k 1 T M d ( t ) M p e a k 1 * ,
A 1 * ( t ) = Γ p e a k 1 φ p e a k 1 T f R ( t ) M p e a k 1 * .
In Equation (16), f R ( t ) is the restoring force vector at time t .

4.1.4. Step 4: Calculation of the Momentary Input Energy and Hysteretic Dissipated Energy in a Half Cycle of the First Modal Response

The momentary input energy of the first modal response per unit mass ( Δ E 1 * / M 1 * ) is calculated from the time-derivative of the equivalent displacement ( D ˙ 1 * ( t ) ) and the ground acceleration ( a g ( t ) ) as
Δ E 1 * M 1 * = t t + Δ t a g ( t ) D ˙ 1 * ( t ) d t .
In Equation (17), t and t + Δ t are, respectively, the beginning and ending times of a half cycle, following the definition of Inoue and his group [59,60,61]. The maximum momentary input energy of the first modal response per unit mass ( Δ E 1 * max / M 1 * ) is defined as the maximum value of Δ E 1 * / M 1 * over the course of the whole seismic input.
Similarly, the hysteretic dissipated energy in a half cycle per unit mass ( Δ E H 1 * / M 1 * ) is calculated as
Δ E H 1 * M 1 * = t t + Δ t A 1 * ( t ) D ˙ 1 * ( t ) d t .
The hysteresis dissipated energy per unit mass in a half cycle at the maximum momentary energy input ( Δ E H 1 * max / M 1 * ) is defined as Δ E H 1 * / M 1 * in the half cycle that Δ E 1 * max / M 1 * occurs.
For the convenience of discussion, the equivalent velocities of Δ E 1 * max / M 1 * and Δ E H 1 * max / M 1 * are defined as
V Δ E 1 * = 2 Δ E 1 * max / M 1 * , V Δ E H 1 * = 2 Δ E H 1 * max / M 1 * .

4.2. Calculation Results

Figure 28 shows the time history of the first modal response for RCDC1 subjected to Uto-NS (Case FM). It is seen that the peak response occurs during UTO0416NS (the second earthquake). The momentary energy input is a maximum in the half cycle from t = 118.60 s to t + Δ t = 119.30 s.
Figure 29 shows the time history of the first modal response for RCDC1 subjected to Uto-NS (Case MF). It is seen that the peak response occurs during UTO0416NS (the first earthquake). The momentary energy input is a maximum in the half cycle from t = 28.62 s to t + Δ t = 29.32 s.
The following observations are made from the comparisons of Case FM and Case MF.
  • The values of the maximum momentary input energy per unit mass ( Δ E 1 * max / M 1 * ) are similar in the two cases. Δ E 1 * max / M 1 * = 0.5430 m2/s2 in Case FM, whereas Δ E 1 * max / M 1 * = 0.5283 m2/s2 in Case MF.
  • The values of the peak equivalent displacement ( D 1 * max ) are similar in the two cases. D 1 * max = 0.2597 m in Case FM, whereas D 1 * max = 0.2431 m in Case MF.
  • The residual equivalent displacement after the first earthquake observed during the interval from 60 to 90 s is small in both cases.
  • The hysteresis loop during UTO0414NS acceleration is different. In Case FM (where UTO0414NS acceleration is used for the first earthquake), the hysteresis loop during UTO0414NS acceleration is not thick (i.e., there is little hysteretic energy dissipation). In contrast, the hysteresis loop during UTO0414NS acceleration is thick (i.e., there is much hysteretic energy dissipation) in Case MF (where UTO0414NS acceleration is used for the second earthquake).
Figure 28 and Figure 29 indicate that the order of ground accelerations in sequential accelerations may affect the cumulative response notably, whereas the effect on the peak response may be limited.
Figure 30 compares the peak equivalent displacement ( D 1 * max ) (a) between Cases F and M, (b) between Cases FM and MF, and (c) between the maximum of Case F, Case M and the maximum of Case FM, Case MF.
The following observations are made for Figure 30.
  • According to the single acceleration shown in Figure 30a, the value of D 1 * max obtained in Case M (mainshock) is larger than that obtained in Case F (foreshock).
  • The effect of the order of ground accelerations in sequential accelerations on D 1 * max is limited, as shown in Figure 30b. The value of D 1 * max obtained in Case FM is similar to that obtained in Case MF.
  • D 1 * max obtained as the maximum of single accelerations (Case F, Case M) is similar to that obtained as the maximum of sequential accelerations (Case FM, Case MF), as shown in Figure 30c.
To understand the third observation, the residual displacement after the first earthquake is investigated. In this study, the residual equivalent displacement after the first earthquake is defined as the absolute value of the equivalent displacement at t = t e n d 1 . Here, t e n d 1 is defined as the end time of the interval (=90 s). The ratio of the residual equivalent displacement is then defined as
r r e s D = | D 1 * ( t e n d 1 ) | / D 1 * max [ 0 , t e n d 1 ] .
In Equation (20), D 1 * max [ 0 , t e n d 1 ] is defined as the local peak equivalent displacement within the range [ 0 , t e n d 1 ] . Note that D 1 * max [ 0 , t e n d 1 ] equals D 1 * max obtained from the first earthquake. In Case FM, D 1 * max [ 0 , t e n d 1 ] equals D 1 * max obtained in Case F.
Figure 31 shows the relation between the ratio r r e s D and the peak equivalent displacement during the first earthquake ( D 1 * max [ 0 , t e n d 1 ] ). It is seen that the ratio r r e s D increases with D 1 * max [ 0 , t e n d 1 ] . However, in Case FM, the ratio r r e s D is less than 0.03, except for MAS-EW. Therefore, the effect of the residual displacement of the first earthquake on D 1 * max is small in Case FM. In addition, the peak equivalent displacement occurs during the first earthquake in Case MF. There is thus no effect of the residual displacement in Case MF, although the ratio r r e s D is greater than that in Case FM.
One reason why the ratio r r e s D is small (e.g., mostly less than 0.03 in Case FM) is that the unloading slope is degrading as shown in Figure 28 and Figure 29 and was pointed out by Ruiz-García [15].
The discussion moves next to the relation between the maximum momentary energy input and the peak displacement of the first modal response. Figure 32 shows the relation of the equivalent velocity of the maximum momentary input energy ( V Δ E 1 * ) and the peak equivalent displacement ( D 1 * max ). A clear relationship is observed between V Δ E 1 * and D 1 * max for both sequential accelerations and single accelerations; i.e., the peak displacement D 1 * max increases with V Δ E 1 * . In addition, the plots obtained for sequential accelerations and single accelerations may be expressed by the same curve. Therefore, the concept of the maximum momentary energy input may be applicable to the prediction of the peak response under sequential accelerations.
For the prediction of the peak response of the structure with viscous damping, the effect of viscous damping on the hysteretic energy in a half cycle must be evaluated. To this end, Figure 33 shows the relation between the ratio V Δ E H 1 * / V Δ E 1 * and the peak equivalent displacement D 1 * max , where V Δ E H 1 * is the equivalent velocity of the hysteretic dissipated energy in a half cycle at the maximum momentary energy input (Equations (18) and (19)). The figure shows that the range of the ratio V Δ E H 1 * / V Δ E 1 * is between 0.80 and 0.97 for BareRC1 and BareRC2 and between 0.90 and 0.98 for RCDC1 and RCDC2. Therefore, the effect of the viscous damping on the ratio V Δ E H 1 * / V Δ E 1 * is small for the building model with damper columns. This is because the contribution of the stiffness of steel damper columns is excluded when calculating the damping matrix in the nonlinear time-history analysis.

4.3. Modeling of the Hysteretic Energy Dissipation of the First Modal Response in a Half Cycle

This section models the hysteresis energy absorption of the first modal response in a half cycle using the results of pushover analysis. The relationship between the equivalent acceleration A 1 * and the equivalent displacement D 1 * obtained from the results of pushover analysis is first validated by comparing with the first modal responses calculated from the results of nonlinear time-history analysis.
Figure 34 compares the A n 1 * D n 1 * relationship obtained from the results of pushover analysis (the same as shown in Figure 7) and the plot of the peak response point ( D 1 * max , A 1 * max ) obtained from the results of time-history analysis. In Figure 34a,b, the plot peak response points fit well with the results of pushover analysis for BareRC1 and BareRC2. However, there are points plotted below the pushover analysis results, especially for BareRC2. In addition, Figure 34c,d show that the plot peak response points fit well with the pushover analysis results for RCDC1 and RCDC2, although there are points plotted above the results of pushover analysis. This is because the strain hardening effect of steel damper columns strengthens the restoring force of the overall building models for RCDC1 and RCDC2. It is concluded from Figure 34 that the A n 1 * D n 1 * relationship can be evaluated from the results of pushover analysis.
Next, the hysteresis energy absorption of the first modal response in a half cycle is modeled using the results of pushover analysis as follows. Let f n R f and f n R d be the restoring force vectors of RC frames and steel damper columns at step n of the pushover analysis, respectively. The contributions of RC frames and steel damper columns to the equivalent acceleration, A n 1 f * and A n 1 d * , respectively, are calculated as
A n 1 f * = Γ n 1 φ n 1 T f n R f M n 1 * = d n T f n R f d n T M 1 ,
A n 1 d * = Γ n 1 φ n 1 T f n R d M n 1 * = d n T f n R d d n T M 1 .
For simplicity, the A n 1 f * D n 1 * and A n 1 d * D n 1 * relationships are idealized as bilinear curves, where the “yield” point of the idealized A 1 f * D 1 * relationship is YF ( D 1 y f * , A 1 y f * ) and that of the idealized A 1 d * D 1 * relationship is YD ( D 1 y d * , A 1 y d * ) . Figure 35 shows the bilinear idealization of the A 1 f * D 1 * relationship for BareRC1 and BareRC2. In the figure, the point P1 ( D P 1 1 * , A P 1 1 * ) is the point at which the first yielding of the RC member occurs, whereas the point P2 ( D P 2 1 * , A P 2 1 * ) is the point at which the equivalent displacement D n 1 * reaches the assumed displacement limit D 1 * l i m i t (=0.2833 m). The equivalent acceleration at point YF ( A 1 y f * ) is determined as the equivalent acceleration at point P2 ( P 2 A 1 * ). The equivalent displacement at point YF ( D 1 y f * ) is then determined as
D 1 y f * = ( A 1 y f * / A P 1 1 f * ) D P 1 1 * = ( A P 2 1 f * / A P 1 1 f * ) D P 1 1 * .
Figure 36 shows the bilinear idealization of the A 1 f * D 1 * and A n 1 d * D n 1 * relationships for models RCDC1 and RCDC2. The bilinear idealization is conducted in the same manner for the two relationships. Specifically, the point P1 is taken as the point at which the first yielding of the steel damper column occurs for the bilinear idealization of the A 1 d * D 1 * relationship.
The hysteresis energy absorption of the first modal response in a half cycle per unit mass ( Δ E μ 1 * / M 1 * ) is then formulated as
Δ E μ 1 * M 1 * ( D 1 * , η ) = A 1 y f * D 1 y f * f F ( μ f , η ) + A 1 y d * D 1 y d * f D ( μ d , η ) ,
μ f = D 1 * / D 1 y f * , μ d = D 1 * / D 1 y d * .
In Equation (24), η is the ratio of displacements in the positive and negative directions. Figure 37 shows simplified models for calculation of the hysteretic dissipated energy during a half cycle of the structural response. From Figure 37, functions f F ( μ f , η ) and f D ( μ d , η ) are, respectively, calculated as
f F ( μ f , η ) = { 1 2 μ f 2 ( 1 η 2 ) : μ f < 1 μ f 1 2 { 1 + ( η μ f ) 2 } : μ f 1   and   0 η < 1 / μ f μ f η μ f : μ f 1   and   1 / μ f η 1 ,
f D ( μ d , η ) = { 1 2 μ d 2 ( 1 η 2 ) : μ d < 1 μ d 1 2 { 1 + ( η μ d ) 2 } : μ d 1   and   0 η < 1 / μ d ( 1 + η ) μ d 2 : μ d 1   and   1 / μ d η 1 .
The equivalent velocity of Δ E μ 1 * / M 1 * is defined as
V Δ E μ 1 * = 2 Δ E μ 1 * / M 1 * .
Figure 38 compares the calculated V Δ E μ 1 * D 1 * relationship obtained from the results of pushover analysis and the plot of the point ( D 1 * max , V Δ E H 1 * ) obtained from the results of time-history analysis. In the figure, the ratio of displacements in the positive and negative directions is set as η = 0, 0.5, and 1.0. It is seen that most of the plots of the point ( D 1 * max , V Δ E H 1 * ) are within the area bounded by the two V Δ E μ 1 * D 1 * curves assuming η = 0 and 1.0, for all four models. Note that for RCDC1 and RCDC2, there are a few points plotted above the V Δ E μ 1 * D 1 * curve assuming η = 0 (upper bound curve). This is because the simplified models shown in Figure 37 neglect the strain-hardening effect of the steel damper columns. The comparisons in Figure 38 confirm that the V Δ E μ 1 * D 1 * relationship of the RC MRF building models with damper columns can be properly evaluated using the simplified model presented herein, for the case of sequential accelerations.

4.4. Summary of Discussions

The discussions presented in this section are summarized as follows.
  • A method of calculating the first modal response from the results of time-history analysis is proposed. When this method is adopted, the first mode vector is assumed from the time history of the relative horizontal displacement vector.
  • The effect of the order of sequential ground accelerations on the peak equivalent displacement of the first modal response ( D 1 * max ) is limited.
  • A clear relationship is observed between the equivalent velocity of the maximum momentary input energy of the first modal response ( V Δ E 1 * ) and D 1 * max for both sequential accelerations and single accelerations; i.e., the peak displacement D 1 * max increases with V Δ E 1 * . Therefore, the concept of the maximum momentary energy input may be applicable to the prediction of the peak response under sequential accelerations.
  • A simplified model capable of evaluating the hysteretic dissipated energy of the first modal response in a half cycle for given equivalent displacement ( D 1 * ) is proposed. The simplified model can be used to evaluate the relationship between the equivalent velocity of the hysteretic dissipated energy in a half cycle ( V Δ E μ 1 * ) and D 1 * with accuracy. The simplified model is applicable to RC MRFs with and without steel damper columns subjected to sequential accelerations.

5. Conclusions

This article investigated the nonlinear response of 10-story RC MRF building models with steel damper columns designed according to a simplified method [3] as a case study of RC MRFs with steel damper columns subjected to seismic sequences. The main results and conclusions of the study as summarized as follows.
  • The peak response of RC MRFs with damper columns subjected to sequential accelerations recorded in the 2016 Kumamoto earthquake is similar to the peak response obtained considering only the mainshock. However, the cumulative strain energy of RC MRFs accumulates in the event of sequential accelerations.
  • The steel damper column is effective for reducing the peak and cumulative responses of RC MRFs in the case of sequential seismic input. The results of nonlinear time-history analysis presented in this study indicate that the installation of steel damper columns can reduce the story drift of RC MRFs and the peak plastic rotation and cumulative strain energy of RC beam ends. However, the designer and structural engineer should pay attention to the behavior of short-span beams in the presence of steel damper columns; i.e., the use of a short-span beam may result in severe damage if its energy absorption capacity is insufficient.
  • The relation of the hysteretic dissipated energy during a half cycle of the structural response and the peak displacement of the first modal response can be properly evaluated using the simple model proposed in this study. The proposed simple model can be applied for RC MRFs with and without hysteresis dampers.
It is obvious that the conclusions drawn in this article are valid only for the limited conditions of four 10-story RC MRFs and six sets of seismic sequences recorded during the foreshock and mainshock of the 2016 Kumamoto earthquake. Despite such limitations, three points can be made:
  • The first point is that the results presented in Section 3 emphasize the importance of the cumulative response of the structures in the case of seismic sequences. Unlike the peak deformation of the members, the cumulative strain energy accumulates in the event of seismic sequences. The evaluation of the cumulative response is important because the peak deformation and the cumulative strain energy are essential parameters for evaluating the structural damage to members.
  • The second point is that the method of calculating the first model response presented here is effective for the post-analysis of nonlinear time-history analysis. This calculation can help the analyst further understand the nonlinear response; e.g., the calculated first modal response can be compared directly with the response of the idealized SDOF model. Note that this calculation is applicable to the post-analysis of the nonlinear time-history analysis and also experimental results, provided the building considered oscillates in the first mode.
  • The third point is that the simplified modeling of the hysteretic dissipated energy during a half cycle of the structural response for a given peak equivalent displacement discussed in Section 4 is essential to the prediction of the peak response of RC MRFs with steel damper columns. Using the maximum momentary input energy spectrum introduced by Inoue and his coauthors [59,60,61], the peak equivalent displacement of the first modal response can be predicted.
Note that the cyclic degradation of members in RC MRFs studied herein is relatively mild. The strength degradation is not considered because there is shear reinforcement in all RC members sufficient to prevent shear failure. However, if both the stiffness and strength degradation due to cyclic loading are notable, the simplified model proposed in Section 4 may need to be revised to include the effect of cyclic degradation. In addition, the effect of cyclic loading would be more pronounced in the case of seismic sequences with long durations, e.g., the records of the 2011 Tohoku earthquake in Japan. In the RC MRF with steel damper columns studied herein, the effectiveness of steel damper columns may deteriorate if the short-span RC beams connected to damper columns are severely damaged. Such issues will be investigated in the next phase of this study.

Funding

This research received financial support from JFE Civil and Construction Corp.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

The original frame model data used in this study were provided by Riho Mukouyama, a former graduate student of Chiba Institute of Technology. Valuable comments from Chizuru Irie and Mitsuhide Yoshinaga, JFE Civil Engineering, and Construction Corp. are appreciated. The ground motions used in this study were obtained from the website of the National Research Institute for Earth Science and Disaster Resilience (NIED) (https://www.kyoshin.bosai.go.jp/, last accessed on 7 January 2022).

Conflicts of Interest

The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

Appendix A. Definition of the Momentary Input Energy

Here, the definition of the momentary input energy is described. Consider the SDOF model with mass ( m ), damping force ( f D ), and restoring force ( f R ). The equation of motion of the SDOF model can be written as
m y ¨ + f D + f R = m a g .
By multiplying both sides of Equation (A1) by y ˙ d t and integrating from 0 to t , the equation of energy balance from time 0 to t is obtained as
E V ( t ) + E D ( t ) + E H ( t ) = E I ( t ) ,
E V ( t ) = 0 t m y ¨ y ˙ d t , E D ( t ) = 0 t f D y ˙ d t , E H ( t ) = 0 t f R y ˙ d t , E I ( t ) = 0 t m a g y ˙ d t .
Here, E V ( t ) is the kinetic energy, E D ( t ) is the damping dissipated energy, E H ( t ) is the hysteretic dissipated energy, and E I ( t ) is the (relative) input energy.
Next, the momentary input energy is defined as follows. Following the work of Inoue and his coauthors [59,60,61], we consider the energy balance during a half cycle of structural response (from t to t + Δ t ). The beginning and end times of a half cycle, t and t + Δ t , respectively, are defined as times when the displacement ( y ( t ) ) is a local maximum. The equation of the energy balance during a half cycle of structural response is expressed as
Δ E V + Δ E D + Δ E H = Δ E ,
{ Δ E V = t t + Δ t m y ¨ y ˙ d t = [ 1 2 m y ˙ 2 ( t ) ] t t + Δ t , Δ E D = t t + Δ t f D y ˙ d t , Δ E H = t t + Δ t f R y ˙ d t , Δ E = t t + Δ t m a g y ˙ d t
In Equation (A4), Δ E V = 0 because the velocities at time t and t + Δ t are zero ( y ˙ ( t ) = y ˙ ( t + Δ t ) = 0 ). Therefore, Equation (A4) can be rewritten as
Δ E D + Δ E H = Δ E .
Equation (A6) indicates that the input energy during a half cycle ( Δ E ) equals to the sum of the damping dissipated energy ( Δ E D ) and hysteretic dissipated energy ( Δ E H ). Based on discussions above, we denote the input energy during a half cycle ( Δ E ) as the momentary input energy and consider as the seismic intensity parameter related to the peak response.

Appendix B. Model Properties

Here, the properties of the members of the four models analyzed in this study are described. Note that the properties of the models RCDC1 and BareRC1 are taken from previous studies [3,4] with some updates.
Table A1 presents the sections at potential hinges of RC members of RCDC1 and BareRC1. Note that the cross-sections of all RC columns have dimensions of 900 mm × 900 mm, which are the same as those of the cross-section at the bottom of the first story. The cross section of the RC beams at the foundation level (Z0) has dimensions of 800 mm × 1900 mm. The yield strength of the longitudinal reinforcement is assumed to be 1.1 × 390 = 429 N/mm2. The assumed compressive strength of concrete is 33 N/mm2 for the first and second stories, 30 N/mm2 for the third to fifth stories, and 27 N/mm2 at and above the sixth story.
Table A2 presents the selected damper columns of RCDC1. The initial normal yield stress of the steel used for the damper panel is assumed to be 205 N/mm2, whereas the normal yield stress after appreciable cyclic loading is assumed to be 300 N/mm2.
Table A1. Sections at potential hinges of RC members for RCDC1 and BareRC1.
Table A1. Sections at potential hinges of RC members for RCDC1 and BareRC1.
MemberLocationWidth (mm)Depth (mm)Longitudinal Reinforcement
Z6 to Z1050090010-D29 (Top and bottom)
Z55509008-D32 (Top and bottom)
BeamZ46009009-D32 (Top and bottom)
Z2 to Z36009008-D35 (Top and bottom)
Z18009009-D38 (Top and bottom)
Column1st Story (Bottom)90090010-D38
Table A2. Steel damper columns for RCDC1.
Table A2. Steel damper columns for RCDC1.
StoryYield StrengthPanel
Thickness
(mm)
Panel
Sectional Area
(mm2)
Column
(mm × mm × mm × mm)
QyDL
(kN)
QyDU
(kN)
1043864163700H-600 × 200 × 12 × 25
962691695290H-600 × 250 × 16 × 32
8755110596380H-700 × 300 × 16 × 28
7968141798180H-900 × 300 × 16 × 28
5 to 612511831910,580H-600 × 250 × 16 × 32
(Doubled)
1 to 415512211912,760H-700 × 300 × 16 × 32
(Doubled)
Table A3 presents the sections at potential hinges of RC members for RCDC2 and BareRC2. Note that the cross-sections of all RC columns have dimensions of 800 mm × 800 mm, which are the same as those of the cross section at the bottom of the first story. The cross section of the RC beams at the foundation level (Z0) has dimensions of 800 mm × 1900 mm. The yield strength of the longitudinal reinforcement and the assumed compressive strength of concrete are the same as for RCDC1 (BareRC1).
Table A4 presents the selected damper columns of RCDC2. The properties of the steel used for the damper panel are the same as those for RCDC1.
Table A3. Sections at potential hinges of RC members for RCDC2 and BareRC2.
Table A3. Sections at potential hinges of RC members for RCDC2 and BareRC2.
MemberLocationWidth (mm)Depth (mm)Longitudinal Reinforcement
Z6 to Z105009007-D29 (Top and bottom)
Z55509007-D29 (Top and bottom)
BeamZ46009006-D32 (Top and bottom)
Z2 to Z360090010-D25 (Top and bottom)
Z18009009-D32 (Top and bottom)
Column1st Story (Bottom)8008008-D35
Table A4. Steel damper columns for RCDC2.
Table A4. Steel damper columns for RCDC2.
StoryYield StrengthPanel
Thickness
(mm)
Panel
Sectional Area
(mm2)
Column
(mm × mm × mm × mm)
QyDL
(kN)
QyDU
(kN)
9 to 1043864163700H-600 × 200 × 12 × 25
7 to 862691695290H-600 × 250 × 16 × 32
6755110596380H-700 × 300 × 16 × 28
5862126197280H-800 × 300 × 16 × 28
2 to 4968141798180H-900 × 300 × 16 × 28
1125118411210,630H-900 × 300 × 19 × 32

References

  1. Wada, A.; Huang, Y.H.; Iwata, M. Passive damping technology for building in Japan. Prog. Struct. Eng. Mater. 2000, 8, 335–350. [Google Scholar] [CrossRef]
  2. Katayama, T.; Ito, S.; Kamura, H.; Ueki, T.; Okamoto, H. Experimental study on hysteretic damper with low yield strength steel under dynamic loading. In Proceedings of the 12th World Conference on Earthquake Engineering, Auckland, New Zealand, 30 January–4 February 2000. [Google Scholar]
  3. Mukouyama, R.; Fujii, K.; Irie, C.; Tobari, R.; Yoshinaga, M.; Miyagawa, K. Displacement-controlled Seismic Design Method of Reinforced Concrete Frame with Steel Damper Column. In Proceedings of the 17th World Conference on Earthquake Engineering, Sendai, Japan, 27 September–2 October 2021. [Google Scholar]
  4. Mukouyama, R. Displacement-controlled Seismic Design of R/C Frame with Steel Damper Column. Master’s Thesis, Chiba Institute of Technology, Chiba, Japan, March 2022. (In Japanese). [Google Scholar]
  5. Fujii, K.; Kato, M. Strength balance of steel damper columns and surrounding beams in reinforced concrete frames. In Proceedings of the 13th International Conference on Earthquake Resistant Engineering Structures, Online, 26–28 May 2021. [Google Scholar]
  6. Mahin, S.A. Effect of duration and aftershocks on inelastic design earthquakes. In Proceedings of the Seventh World Conference on Earthquake Engineering, Istanbul, Turkey, 8–13 September 1980. [Google Scholar]
  7. Amadio, C.; Fragiacomo, M.; Rajgelj, S. The effects of repeated earthquake ground motions on the non-linear response of SDOF systems. Earthq. Eng. Struct. Dyn. 2003, 32, 291–308. [Google Scholar] [CrossRef]
  8. Hatzigeorgiou, G.D.; Beskos, D.E. Inelastic displacement ratios for SDOF structures subjected to repeated earthquakes. Eng. Struct. 2009, 31, 2744–2755. [Google Scholar] [CrossRef]
  9. Hatzigeorgiou, G.D. Behavior factors for nonlinear structures subjected to multiple near-fault earthquakes. Comput. Struct. 2010, 88, 309–321. [Google Scholar] [CrossRef]
  10. Hatzigeorgiou, G.D. Ductility demand spectra for multiple near- and far-fault earthquakes. Soil Dyn. Earthq. Eng. 2010, 30, 170–183. [Google Scholar] [CrossRef]
  11. Hatzigeorgiou, G.D.; Liolios, A.A. Nonlinear behaviour of RC frames under repeated strong ground motions. Soil Dyn. Earthq. Eng. 2010, 30, 1010–1025. [Google Scholar] [CrossRef]
  12. Faisal, A.; Majid, T.A.; Hatzigeorgiou, G.D. Investigation of story ductility demands of inelastic concrete frames subjected to repeated earthquakes. Soil Dyn. Earthq. Eng. 2013, 44, 42–53. [Google Scholar] [CrossRef]
  13. Hatzivassiliou, M.; Hatzigeorgiou, G.D. Seismic sequence effects on three-dimensional reinforced concrete buildings. Soil Dyn. Earthq. Eng. 2015, 72, 77–80. [Google Scholar] [CrossRef]
  14. Ruiz-García, J.; Negrete-Manriquez, J.C. Evaluation of drift demands in existing steel frames under as-recorded far-field and near-fault mainshock–aftershock seismic sequences. Eng. Struct. 2011, 33, 621–634. [Google Scholar] [CrossRef]
  15. Ruiz-García, J. Issues on the response of existing buildings under mainshock-aftershock seismic sequences. In Proceedings of the 15th World Conference on Earthquake Engineering, Lisbon, Portugal, 24–28 September 2012. [Google Scholar]
  16. Ruiz-García, J.; Terán-Gilmore, A.; Díaz, G. Response of essential facilities under narrow-band mainshock-aftershock seismic sequences. In Proceedings of the 15th World Conference on Earthquake Engineering, Lisbon, Portugal, 24–28 September 2012. [Google Scholar]
  17. Ruiz-García, J. Mainshock-aftershock ground motion features and their influence in building’s seismic response. J. Earthq. Eng. 2012, 16, 719–737. [Google Scholar] [CrossRef]
  18. Ruiz-García, J.; Marín, M.V.; Terán-Gilmore, A. Effect of seismic sequences on the response of RC buildings located in soft soil sites. In Proceedings of the 2013 World Congress on Advances in Structural Engineering and Mechanics (ASEM13), Jeju, Korea, 8–12 September 2013. [Google Scholar]
  19. Ruiz-García, J. Three-dimensional building response under seismic sequences. In Proceedings of the 2013 World Congress on Advances in Structural Engineering and Mechanics (ASEM13), Jeju, Korea, 8–12 September 2013. [Google Scholar]
  20. Díaz-Martínez, G.; Ruiz-García, J.; Terán-Gilmore, A. Response of structures to seismic sequences corresponding to Mexican soft soils. Earthq. Struct. 2014, 7, 1241–1258. [Google Scholar] [CrossRef]
  21. Guerrero, H.; Ruíz-García, J.; Escobar, J.A.; Terán-Gilmore, A. Response to seismic sequences of short-period structures equipped with buckling-restrained braces located on the lakebed zone of Mexico City. J. Constr. Steel Res. 2017, 137, 37–51. [Google Scholar] [CrossRef]
  22. Ruiz-García, J.; Yaghmaei-Sabegh, S.; Bojórquez, E. Three-dimensional response of steel moment-resisting buildings under seismic sequences. Eng. Struct. 2018, 175, 399–414. [Google Scholar] [CrossRef]
  23. Goda, K. Nonlinear response potential of mainshock–aftershock sequences from Japanese earthquakes. Bull. Seismol. Soc. Am. 2012, 102, 2139–2156. [Google Scholar] [CrossRef]
  24. Goda, K. Peak ductility demand of mainshock-aftershock sequences for Japanese earthquakes. In Proceedings of the15th World Conference on Earthquake Engineering, Lisbon, Portugal, 24–28 September 2012. [Google Scholar]
  25. Goda, K. Record selection for aftershock incremental dynamic analysis. Earthq. Eng. Struct. Dyn. 2014, 44, 1157–1162. [Google Scholar] [CrossRef]
  26. Goda, K.; Wenzel, F.; De Risi, R. Empirical assessment of non-linear seismic demand of mainshock–aftershock ground-motion sequences for Japanese earthquakes. Front. Built Environ. 2015, 1, 6. [Google Scholar] [CrossRef]
  27. Tesfamariam, S.; Goda, K.; Mondal, G. Seismic vulnerability of reinforced concrete frame with unreinforced masonry infill due to mainshock–aftershock earthquake sequences. Earthq. Spectra 2015, 31, 1427–1449. [Google Scholar] [CrossRef]
  28. Tesfamariam, S.; Goda, K. Seismic performance evaluation framework considering maximum and residual inter-story drift ratios: Application to non-code conforming reinforced concrete buildings in Victoria, BC, Canada. Front. Built Environ. 2015, 1, 18. [Google Scholar] [CrossRef] [Green Version]
  29. Tesfamariam, S.; Goda, K. Energy-based seismic risk evaluation of tall reinforced concrete building in Vancouver, BC, Canada, under Mw9 megathrust subduction earthquakes and aftershocks. Front. Built Environ. 2017, 3, 29. [Google Scholar] [CrossRef] [Green Version]
  30. Zhai, C.H.; Wen, W.P.; Li, S.; Xie, L.L. The influences of seismic sequence on the inelastic SDOF systems. In Proceedings of the 15th World Conference on Earthquake Engineering, Lisbon, Portugal, 24–28 September 2012. [Google Scholar]
  31. Zhai, C.H.; Wen, W.P.; Chen, Z.Q.; Li, S.; Xie, L.L. Damage spectra for the mainshock–aftershock sequence-type ground motions. Soil Dyn. Earthq. Eng. 2013, 45, 1–12. [Google Scholar] [CrossRef]
  32. Zhai, C.H.; Wen, W.P.; Li, S.; Chen, Z.Q.; Chang, Z.; Xie, L.L. The damage investigation of inelastic SDOF structure under the mainshock–aftershock sequence-type ground motions. Soil Dyn. Earthq. Eng. 2014, 59, 30–41. [Google Scholar] [CrossRef]
  33. Di Sarno, L. Effects of multiple earthquakes on inelastic structural response. Eng. Struct. 2013, 56, 673–681. [Google Scholar] [CrossRef]
  34. Di Sarno, L.; Amiric, S. Period elongation of deteriorating structures under mainshock-aftershock sequences. Eng. Struct. 2019, 196, 109341. [Google Scholar] [CrossRef]
  35. Di Sarno, L.; Wu, J.R. Fragility assessment of existing low-rise steel moment-resisting frames with masonry infills under mainshock-aftershock earthquake sequences. Bull. Earthq. Eng. 2021, 19, 2483–2504. [Google Scholar] [CrossRef]
  36. Di Sarno, L.; Pugliese, F. Effects of mainshock-aftershock sequences on fragility analysis of RC buildings with ageing. Eng. Struct. 2021, 232, 111837. [Google Scholar] [CrossRef]
  37. Abdelnaby, A.E.; Elnashai, A.S. Performance of degrading reinforced concrete frame systems under the Tohoku and Christchurch Earthquake sequences. J. Earthq. Eng. 2014, 18, 1009–1036. [Google Scholar] [CrossRef]
  38. Abdelnaby, A.E.; Elnashai, A.S. Numerical modeling and analysis of RC frames subjected to multiple earthquakes. Earthq. Struct. 2015, 9, 957–981. [Google Scholar] [CrossRef]
  39. Abdelnaby, A.E. Fragility Curves for RC Frames Subjected to Tohoku Mainshock-Aftershocks Sequences. J. Earthq. Eng. 2016, 22, 902–920. [Google Scholar] [CrossRef]
  40. Hosseinpour, F.; Abdelnaby, A.E. Effect of different aspects of multiple earthquakes on the nonlinear behavior of RC structures. Soil Dyn. Earthq. Eng. 2017, 92, 706–725. [Google Scholar] [CrossRef]
  41. Oyguc, R.; Toros, C.; Abdelnaby, A.E. Seismic behavior of irregular reinforced-concrete structures under multiple earthquake excitations. Soil Dyn. Earthq. Eng. 2018, 104, 15–32. [Google Scholar] [CrossRef]
  42. Kagermanov, A.; Gee, R. Cyclic pushover method for seismic performance assessment under multiple earthquakes. In Proceedings of the 16th European Conference on Earthquake Engineering, Thessaloniki, Greece, 18–21 June 2018. [Google Scholar]
  43. Kagermanov, A.; Gee, R. Cyclic pushover method for seismic assessment under multiple earthquakes. Earthq. Spectra 2019, 35, 1541–1588. [Google Scholar] [CrossRef]
  44. Yang, F.; Wang, G.; Ding, Y. Damage demands evaluation of reinforced concrete frame structure subjected to near-fault seismic sequences. Nat. Hazards 2019, 97, 841–860. [Google Scholar] [CrossRef]
  45. Yang, F.; Wang, G.; Li, M. Evaluation of the seismic retrofitting of mainshock-damaged reinforced concrete frame structure using steel braces with soft steel dampers. Appl. Sci. 2021, 11, 841. [Google Scholar] [CrossRef]
  46. Yaghmaei-Sabegh, S.; Mahdipour-Moghanni, R. State-dependent fragility curves using real and artificial earthquake sequences. Asian J. Civ. Eng. 2019, 20, 619–625. [Google Scholar] [CrossRef]
  47. Qiao, Y.H.; Lu, D.G.; Hui Yu, X.H. Shaking table tests of a reinforced concrete frame subjected to mainshock-aftershock sequences. J. Earthq. Eng. 2020, 1–20. [Google Scholar] [CrossRef]
  48. Orlacchio, M.; Baltzopoulos, G.; Iervolino, I. State-dependent seismic fragility via pushover analysis. In Proceedings of the 17th World Conference on Earthquake Engineering, Sendai, Japan, 20 September–2 November 2021. [Google Scholar]
  49. Hoveidae, N.; Radpour, S. Performance evaluation of buckling-restrained braced frames under repeated earthquakes. Bull. Earthq. Eng. 2021, 19, 241–262. [Google Scholar] [CrossRef]
  50. Pirooz, R.M.; Habashi, S.; Massumi, A. Required time gap between mainshock and aftershock for dynamic analysis of structures. Bull. Earthq. Eng. 2021, 19, 2643–2670. [Google Scholar] [CrossRef]
  51. Park, Y.J.; Ang, A.H.S. Mechanistic seismic damage model for reinforced concrete. J. Struct. Eng. ASCE 1985, 111, 722–739. [Google Scholar] [CrossRef]
  52. Akiyama, H. Earthquake–Resistant Limit–State Design for Buildings; University of Tokyo Press: Tokyo, Japan, 1985. [Google Scholar]
  53. Uang, C.M.; Bertero, V.V. Evaluation of seismic energy in structures. Earthq. Eng. Struct. Dyn. 1990, 19, 77–90. [Google Scholar] [CrossRef]
  54. Fajfar, P.; Vidic, T. Consistent inelastic design spectra: Hysteretic and input energy. Earthq. Eng. Struct. Dyn. 1994, 23, 523–537. [Google Scholar] [CrossRef]
  55. Sucuoǧlu, H.; Nurtuǧ, A. Earthquake ground motion characteristics and seismic energy dissipation. Earthq. Eng. Struct. Dyn. 1995, 24, 1195–1213. [Google Scholar] [CrossRef]
  56. Nurtuǧ, A.; Sucuoǧlu, H. Prediction of seismic energy dissipation in SDOF systems. Earthq. Eng. Struct. Dyn. 1995, 24, 1215–1223. [Google Scholar] [CrossRef]
  57. Chai, Y.H.; Fajfar, P.; Romstad, K.M. Formulation of duration-dependent inelastic seismic design spectrum. J. Struct. Eng. 1998, 124, 913–921. [Google Scholar] [CrossRef]
  58. Cheng, Y.; Lucchini, A.; Mollaioli, F. Ground-motion prediction equations for constant-strength and constant-ductility input energy spectra. Bull. Earthq. Eng. 2020, 18, 37–55. [Google Scholar] [CrossRef]
  59. Nakamura, T.; Hori, N.; Inoue, N. Evaluation of damaging properties of ground motions and estimation of maximum displacement based on momentary input energy. J. Struct. Constr. Eng. Trans. AIJ 1998, 63, 65–72. (In Japanese) [Google Scholar] [CrossRef]
  60. Inoue, N.; Wenliuhan, H.; Kanno, H.; Hori, N.; Ogawa, J. Shaking Table Tests of Reinforced Concrete Columns Subjected to Simulated Input Motions with Different Time Durations. In Proceedings of the 12th World Conference on Earthquake Engineering, Auckland, New Zealand, 30 January–4 February 2000. [Google Scholar]
  61. Hori, N.; Inoue, N. Damaging properties of ground motion and prediction of maximum response of structures based on momentary energy input. Earthq. Eng. Struct. Dyn. 2002, 31, 1657–1679. [Google Scholar] [CrossRef]
  62. Fujii, K.; Kanno, H.; Nishida, T. Formulation of the time-varying function of momentary energy input to a SDOF system by Fourier series. J. Jpn. Assoc. Earthq. Eng. 2021, 21, 28–47. [Google Scholar] [CrossRef]
  63. Fujii, K.; Murakami, Y. Bidirectional momentary energy input to a one-mass two-DOF system. In Proceedings of the 17th World Conference on Earthquake Engineering, Sendai, Japan, 20 September–2 November 2021. [Google Scholar]
  64. Fujii, K. Bidirectional seismic energy input to an isotropic nonlinear one-mass two-degree-of-freedom system. Buildings 2021, 11, 143. [Google Scholar] [CrossRef]
  65. Fujii, K.; Masuda, T. Application of Mode-Adaptive Bidirectional Pushover Analysis to an Irregular Reinforced Concrete Building Retrofitted via Base Isolation. Appl. Sci. 2021, 11, 9829. [Google Scholar] [CrossRef]
  66. Strong-Motion Seismograph Networks (K-NET, KiK-net). Available online: https://www.kyoshin.bosai.go.jp/ (accessed on 7 January 2022).
  67. AIJ. AIJ Standard for Lateral Load-Carrying Capacity Calculation of Reinforced Concrete Structures (Draft); Architectural Institute of Japan: Tokyo, Japan, 2016. (In Japanese) [Google Scholar]
  68. Sugano, S.; Koreishi, I. An empirical evaluation of inelastic behaviour of structural elements in reinforced concrete frames subjected to lateral forces. In Proceedings of the 5th World Conference on Earthquake Engineering, Rome, Italy, 25–29 June 1974. [Google Scholar]
  69. Muto, K.; Hisada, T.; Tsugawa, T.; Bessho, S. Earthquake resistant design of a 20 story reinforced concrete buildings. In Proceedings of the 5th World Conference on Earthquake Engineering, Rome, Italy, 25–29 June 1974. [Google Scholar]
  70. Umemura, H.; Ichinose, T.; Ohashi, K.; Maekawa, J. Hysteresis Model for RC Members Considering Strength Deterioration. In Proceedings of the Japan Concrete Institute, Tsukuba, Japan, 19–21 June 2002. (In Japanese). [Google Scholar]
  71. Baber, T.T.; Noori, M.N. Random Vibration of Degrading, Pinching Systems. J. Eng. Mech. ASCE 1985, 111, 1010–1026. [Google Scholar] [CrossRef]
  72. Ono, Y.; Kaneko, H. Constitutive rules of the steel damper and source code for the analysis program. In Proceedings of the Passive Control Symposium 2001, Yokohama, Japan, 14–15 December 2001. (In Japanese). [Google Scholar]
  73. Vaiana, N.; Sessa, S.; Marmo, F.; Rosati, L. A class of uniaxial phenomenological models for simulating hysteretic phenomena in rate-independent mechanical systems and materials. Nonlinear Dyn. 2018, 93, 1647–1669. [Google Scholar] [CrossRef]
  74. Vaiana, N.; Napolitano, C.; Rosati, L. Some recent advances on the modeling of the hysteretic behavior of rate-independent passive energy dissipation devices. In Proceedings of the 8th ECCOMAS Thematic Conference on Computational Methods in Structural Dynamics and Earthquake Engineering (COMPDYN 2021), Athens, Greece, 28–30 June 2021. [Google Scholar]
  75. Nuzzo, I.; Losanno, D.; Caterino, N.; Serino, G.; Rotondo, L.M.B. Experimental and analytical characterization of steel shear links for seismic energy dissipation. Eng. Struct. 2018, 172, 405–418. [Google Scholar] [CrossRef]
  76. Fujii, K. Prediction of the largest peak nonlinear seismic response of asymmetric buildings under bi-directional excitation using pushover analyses. Bull. Earthq. Eng. 2014, 12, 909–938. [Google Scholar] [CrossRef]
  77. Building Center of Japan (BCJ). The Building Standard Law of Japan on CD-ROM; The Building Center of Japan: Tokyo, Japan, 2016. [Google Scholar]
  78. Miranda, E. Evaluation of site-dependent inelastic seismic design spectra. J. Struct. Eng. ASCE 1993, 119, 1319–1338. [Google Scholar] [CrossRef]
Figure 1. Design collapse mechanism of MRFs: (a) traditional MRF designed according to the strong-column weak-beam concept; (b) MRF with steel damper columns.
Figure 1. Design collapse mechanism of MRFs: (a) traditional MRF designed according to the strong-column weak-beam concept; (b) MRF with steel damper columns.
Buildings 12 00275 g001
Figure 2. Design concept of the RC MRF with steel damper columns: (a) conceptual structural elevation; (b) details of the RC beam–steel damper column joint.
Figure 2. Design concept of the RC MRF with steel damper columns: (a) conceptual structural elevation; (b) details of the RC beam–steel damper column joint.
Buildings 12 00275 g002
Figure 3. RC MRF building models with damper columns: (a1) structural plan of RCDC1; (a2) structural elevation of RCDC1 (frame Y2); (b1) structural plan of RCDC2; (b2) structural elevation of RCDC2 (frame Y2).
Figure 3. RC MRF building models with damper columns: (a1) structural plan of RCDC1; (a2) structural elevation of RCDC1 (frame Y2); (b1) structural plan of RCDC2; (b2) structural elevation of RCDC2 (frame Y2).
Buildings 12 00275 g003
Figure 4. Modeling of the MRF with a damper column: (a) potential hinge and non-hinge at the beam end and (b) structural model.
Figure 4. Modeling of the MRF with a damper column: (a) potential hinge and non-hinge at the beam end and (b) structural model.
Buildings 12 00275 g004
Figure 5. Envelope of the force–deformation relationship of each member: (a) flexural behavior of the potential hinge of RC members; (b) flexural behavior of non-hinge ends of RC members; (c) shear behavior of the damper panel in the steel damper column.
Figure 5. Envelope of the force–deformation relationship of each member: (a) flexural behavior of the potential hinge of RC members; (b) flexural behavior of non-hinge ends of RC members; (c) shear behavior of the damper panel in the steel damper column.
Buildings 12 00275 g005
Figure 6. Hysteresis rule of each member: (a) flexural behavior of RC members and (b) shear behavior of the damper panel in the steel damper column.
Figure 6. Hysteresis rule of each member: (a) flexural behavior of RC members and (b) shear behavior of the damper panel in the steel damper column.
Buildings 12 00275 g006
Figure 7. Relationship between the equivalent acceleration A n 1 * and equivalent displacement D n 1 * of each model: (a) BareRC1; (b) BareRC2; (c) RCDC1; (d) RCDC2.
Figure 7. Relationship between the equivalent acceleration A n 1 * and equivalent displacement D n 1 * of each model: (a) BareRC1; (b) BareRC2; (c) RCDC1; (d) RCDC2.
Buildings 12 00275 g007
Figure 8. Primary and shear wave profiles: (a) KMM; (b) UTO; (c) MAS.
Figure 8. Primary and shear wave profiles: (a) KMM; (b) UTO; (c) MAS.
Buildings 12 00275 g008
Figure 9. Ground motion time histories: (a) KMM; (b) UTO; (c) MAS.
Figure 9. Ground motion time histories: (a) KMM; (b) UTO; (c) MAS.
Buildings 12 00275 g009
Figure 10. Comparisons of the pseudo-velocity spectrum of the observed ground motions (foreshock, aftershock) and the design earthquake: (a) KMM; (b) UTO; (c) MAS.
Figure 10. Comparisons of the pseudo-velocity spectrum of the observed ground motions (foreshock, aftershock) and the design earthquake: (a) KMM; (b) UTO; (c) MAS.
Buildings 12 00275 g010
Figure 11. Distribution of the peak relative displacement (EW component of KMM-EW): (a) BareRC1; (b) BareRC2; (c) RCDC1; (d) RCDC2.
Figure 11. Distribution of the peak relative displacement (EW component of KMM-EW): (a) BareRC1; (b) BareRC2; (c) RCDC1; (d) RCDC2.
Buildings 12 00275 g011
Figure 12. Distribution of the peak relative displacement (MAS-EW): (a) BareRC1; (b) BareRC2; (c) RCDC1; (d) RCDC2.
Figure 12. Distribution of the peak relative displacement (MAS-EW): (a) BareRC1; (b) BareRC2; (c) RCDC1; (d) RCDC2.
Buildings 12 00275 g012
Figure 13. Distribution of the peak story drift (KMM-EW): (a) BareRC1; (b) BareRC2; (c) RCDC1; (d) RCDC2.
Figure 13. Distribution of the peak story drift (KMM-EW): (a) BareRC1; (b) BareRC2; (c) RCDC1; (d) RCDC2.
Buildings 12 00275 g013
Figure 14. Distribution of the peak story drift (MAS-EW): (a) BareRC1; (b) BareRC2; (c) RCDC1; (d) RCDC2.
Figure 14. Distribution of the peak story drift (MAS-EW): (a) BareRC1; (b) BareRC2; (c) RCDC1; (d) RCDC2.
Buildings 12 00275 g014
Figure 15. Comparisons of the peak story drift in the case of single and sequential acceleration: (a) BareRC1; (b) BareRC2; (c) RCDC1; (d) RCDC2.
Figure 15. Comparisons of the peak story drift in the case of single and sequential acceleration: (a) BareRC1; (b) BareRC2; (c) RCDC1; (d) RCDC2.
Buildings 12 00275 g015
Figure 16. Distribution of the peak plastic hinge rotation at the beam end (MAS-EW): (a) BareRC1; (b) BareRC2; (c) RCDC1; (d) RCDC2.
Figure 16. Distribution of the peak plastic hinge rotation at the beam end (MAS-EW): (a) BareRC1; (b) BareRC2; (c) RCDC1; (d) RCDC2.
Buildings 12 00275 g016
Figure 17. Comparisons of the peak plastic hinge rotation at the beam end in cases of single and sequential acceleration: (a) BareRC1; (b) BareRC2; (c) RCDC1; (d) RCDC2.
Figure 17. Comparisons of the peak plastic hinge rotation at the beam end in cases of single and sequential acceleration: (a) BareRC1; (b) BareRC2; (c) RCDC1; (d) RCDC2.
Buildings 12 00275 g017
Figure 18. Distribution of the peak shear strain of the damper panel (MAS-EW): (a) RCDC1; (b) RCDC2.
Figure 18. Distribution of the peak shear strain of the damper panel (MAS-EW): (a) RCDC1; (b) RCDC2.
Buildings 12 00275 g018
Figure 19. Comparisons of the peak shear strain of the damper panel in the cases of single and sequential acceleration: (a) RCDC1; (b) RCDC2.
Figure 19. Comparisons of the peak shear strain of the damper panel in the cases of single and sequential acceleration: (a) RCDC1; (b) RCDC2.
Buildings 12 00275 g019
Figure 20. Comparisons of the cumulative energy per unit mass of the overall building in the cases of single acceleration and sequential acceleration: (a1) KMM-EW; (a2) KMM-NS; (b1) UTO-EW; (b2) UTO-NS; (c1) MAS-EW; (c2) MAS-NS.
Figure 20. Comparisons of the cumulative energy per unit mass of the overall building in the cases of single acceleration and sequential acceleration: (a1) KMM-EW; (a2) KMM-NS; (b1) UTO-EW; (b2) UTO-NS; (c1) MAS-EW; (c2) MAS-NS.
Buildings 12 00275 g020
Figure 21. Comparisons of the cumulative strain energy of the plastic hinges at the beam end in the cases of single and sequential acceleration: (a) BareRC1; (b) BareRC2.
Figure 21. Comparisons of the cumulative strain energy of the plastic hinges at the beam end in the cases of single and sequential acceleration: (a) BareRC1; (b) BareRC2.
Buildings 12 00275 g021
Figure 22. Comparisons of the cumulative strain energy of the plastic hinge at the beam end for single and sequential accelerations: (a1) RCDC1 (Frame Y3); (a2) RCDC1 (Frame Y2); (b1) RCDC2 (Frame Y3); (b2) RCDC2 (Frame Y2).
Figure 22. Comparisons of the cumulative strain energy of the plastic hinge at the beam end for single and sequential accelerations: (a1) RCDC1 (Frame Y3); (a2) RCDC1 (Frame Y2); (b1) RCDC2 (Frame Y3); (b2) RCDC2 (Frame Y2).
Buildings 12 00275 g022
Figure 23. Comparisons of the cumulative strain energy of the damper panel in the cases of single and sequential acceleration: (a) RCDC1; (b) RCDC2.
Figure 23. Comparisons of the cumulative strain energy of the damper panel in the cases of single and sequential acceleration: (a) RCDC1; (b) RCDC2.
Buildings 12 00275 g023
Figure 24. Comparisons of the peak story drift for models without and with damper columns: (a) BareRC1 and RCDC1; (b) BareRC2 and RCDC2.
Figure 24. Comparisons of the peak story drift for models without and with damper columns: (a) BareRC1 and RCDC1; (b) BareRC2 and RCDC2.
Buildings 12 00275 g024
Figure 25. Comparisons of the peak plastic rotation at beam ends in models without and with damper columns: (a) BareRC1 and RCDC1; (b) BareRC2 and RCDC2.
Figure 25. Comparisons of the peak plastic rotation at beam ends in models without and with damper columns: (a) BareRC1 and RCDC1; (b) BareRC2 and RCDC2.
Buildings 12 00275 g025
Figure 26. Comparisons of the cumulative strain energy of RC MRF for models without and with damper columns: (a) BareRC1 and RCDC1; (b) BareRC2 and RCDC2.
Figure 26. Comparisons of the cumulative strain energy of RC MRF for models without and with damper columns: (a) BareRC1 and RCDC1; (b) BareRC2 and RCDC2.
Buildings 12 00275 g026
Figure 27. Comparisons of the cumulative strain energy at the plastic hinge at the beam end for models without and with damper columns: (a) BareRC1 and RCDC1; (b) BareRC2 and RCDC2.
Figure 27. Comparisons of the cumulative strain energy at the plastic hinge at the beam end for models without and with damper columns: (a) BareRC1 and RCDC1; (b) BareRC2 and RCDC2.
Buildings 12 00275 g027
Figure 28. Calculated first modal response of RCDC1 subjected to UTO-NS (Case FM): (a) time history of the ground acceleration; (b) time history of the equivalent displacement; (c) time history of the momentary energy input per unit mass; (d) hysteresis loop of the A 1 * D 1 * relationship (0–90 s); (e) hysteresis loop of the A 1 * D 1 * relationship (90–150 s); (f) hysteresis loop of the A 1 * D 1 * relationship (at the time of the maximum momentary energy input).
Figure 28. Calculated first modal response of RCDC1 subjected to UTO-NS (Case FM): (a) time history of the ground acceleration; (b) time history of the equivalent displacement; (c) time history of the momentary energy input per unit mass; (d) hysteresis loop of the A 1 * D 1 * relationship (0–90 s); (e) hysteresis loop of the A 1 * D 1 * relationship (90–150 s); (f) hysteresis loop of the A 1 * D 1 * relationship (at the time of the maximum momentary energy input).
Buildings 12 00275 g028
Figure 29. Calculated first modal response of RCDC1 subjected to UTO-NS (Case MF): (a) time history of the ground acceleration; (b) time history of the equivalent displacement; (c) time history of the momentary energy input per unit mass; (d) hysteresis loop of the A 1 * D 1 * relationship (0–90 s); (e) hysteresis loop of the A 1 * D 1 * relationship (90–150 s); (f) hysteresis loop of the A 1 * D 1 * relationship (at the time of the maximum momentary energy input).
Figure 29. Calculated first modal response of RCDC1 subjected to UTO-NS (Case MF): (a) time history of the ground acceleration; (b) time history of the equivalent displacement; (c) time history of the momentary energy input per unit mass; (d) hysteresis loop of the A 1 * D 1 * relationship (0–90 s); (e) hysteresis loop of the A 1 * D 1 * relationship (90–150 s); (f) hysteresis loop of the A 1 * D 1 * relationship (at the time of the maximum momentary energy input).
Buildings 12 00275 g029
Figure 30. Comparisons of the peak equivalent displacement: (a) single acceleration; (b) sequential acceleration; (c) maximum of sequential accelerations and maximum of single accelerations.
Figure 30. Comparisons of the peak equivalent displacement: (a) single acceleration; (b) sequential acceleration; (c) maximum of sequential accelerations and maximum of single accelerations.
Buildings 12 00275 g030
Figure 31. Relations of the ratio r r e s D and D 1 * max [ 0 , t e n d 1 ] : (a) BareRC1; (b) BareRC2; (c) RCDC1; (d) RCDC2.
Figure 31. Relations of the ratio r r e s D and D 1 * max [ 0 , t e n d 1 ] : (a) BareRC1; (b) BareRC2; (c) RCDC1; (d) RCDC2.
Buildings 12 00275 g031
Figure 32. Relations of V Δ E 1 * and D 1 * max : (a) BareRC1; (b) BareRC2; (c) RCDC1; (d) RCDC2.
Figure 32. Relations of V Δ E 1 * and D 1 * max : (a) BareRC1; (b) BareRC2; (c) RCDC1; (d) RCDC2.
Buildings 12 00275 g032
Figure 33. Relations of the ratio V Δ E H 1 * / V Δ E 1 * and D 1 * max : (a) BareRC1; (b) BareRC2; (c) RCDC1; (d) RCDC2.
Figure 33. Relations of the ratio V Δ E H 1 * / V Δ E 1 * and D 1 * max : (a) BareRC1; (b) BareRC2; (c) RCDC1; (d) RCDC2.
Buildings 12 00275 g033
Figure 34. Comparisons of the A n 1 * D n 1 * relationship obtained from the results of pushover analysis and the plot of the peak response point ( D 1 * max , A 1 * max ) obtained from the results of time-history analysis: (a) BareRC1; (b) BareRC2; (c) RCDC1; (d) RCDC2.
Figure 34. Comparisons of the A n 1 * D n 1 * relationship obtained from the results of pushover analysis and the plot of the peak response point ( D 1 * max , A 1 * max ) obtained from the results of time-history analysis: (a) BareRC1; (b) BareRC2; (c) RCDC1; (d) RCDC2.
Buildings 12 00275 g034
Figure 35. Bilinear idealization of the A 1 f * D 1 * relationship from the results of pushover analysis: (a) BareRC1; (b) BareRC2.
Figure 35. Bilinear idealization of the A 1 f * D 1 * relationship from the results of pushover analysis: (a) BareRC1; (b) BareRC2.
Buildings 12 00275 g035
Figure 36. Bilinear idealization of the A 1 f * D 1 * and A 1 d * D 1 * relationships from the results of pushover analysis: (a) RCDC1; (b) RCDC2.
Figure 36. Bilinear idealization of the A 1 f * D 1 * and A 1 d * D 1 * relationships from the results of pushover analysis: (a) RCDC1; (b) RCDC2.
Buildings 12 00275 g036
Figure 37. Simplified models for calculation of the hysteretic dissipated energy during a half cycle of the structural response: (a) contribution of RC MRFs; (b) contribution of steel damper columns.
Figure 37. Simplified models for calculation of the hysteretic dissipated energy during a half cycle of the structural response: (a) contribution of RC MRFs; (b) contribution of steel damper columns.
Buildings 12 00275 g037
Figure 38. Comparisons of the V Δ E μ 1 * D 1 * relationship obtained from the results of pushover analysis and the plot of the point ( D 1 * max , V Δ E H 1 * ) obtained from the results of time-history analysis: (a) BareRC1; (b) BareRC2; (c) RCDC1; (d) RCDC2.
Figure 38. Comparisons of the V Δ E μ 1 * D 1 * relationship obtained from the results of pushover analysis and the plot of the point ( D 1 * max , V Δ E H 1 * ) obtained from the results of time-history analysis: (a) BareRC1; (b) BareRC2; (c) RCDC1; (d) RCDC2.
Buildings 12 00275 g038
Table 1. Values of parameter χ representing the cyclic degradation of RC members.
Table 1. Values of parameter χ representing the cyclic degradation of RC members.
BareRC1, RCDC1BareRC2, RCDC2
Short-span Beam
(beams connected to damper column)
χ = 0.074
(only RCDC1)
χ = 0.081
(only RCDC2)
Long-span Beam
(other beams)
χ = 0.034 χ = 0.049
Column
(bottom-end of the first story)
χ = 0.060 χ = 0.046
Table 2. Natural periods of the first three modes in the initial stage.
Table 2. Natural periods of the first three modes in the initial stage.
BareRC1BareRC2RCDC1RCDC2
T 1 e (s)0.85470.71060.70440.6152
T 2 e (s)0.28340.24270.23440.2094
T 3 e (s)0.15770.13690.13300.1195
Table 3. Ground motions.
Table 3. Ground motions.
Station NameEvent DateDistanceGround Motion IDPGA (m/s2)
EWNS
K-NET Kumamoto (KMM)14 April 2016 6 kmKMM04143.8145.744
16 April 20165 kmKMM04166.1628.272
K-NET Uto
(UTO)
14 April 201615 kmUTO04143.0422.635
16 April 201612 kmUTO04167.7116.515
KIK-NET Mashiki (MAS)14 April 20166 kmMAS04149.2507.598
16 April 20167 kmMAS041611.5696.530
Table 4. Cases of ground motion.
Table 4. Cases of ground motion.
CaseAcceleration Sequence
Case-FForeshock (0414) only
Case-FMForeshock (0414) + Mainshock (0416)
Case-MMainshock (0416) only
Case-MFMainshock (0416) + Foreshock (0414)
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Fujii, K. Peak and Cumulative Response of Reinforced Concrete Frames with Steel Damper Columns under Seismic Sequences. Buildings 2022, 12, 275. https://doi.org/10.3390/buildings12030275

AMA Style

Fujii K. Peak and Cumulative Response of Reinforced Concrete Frames with Steel Damper Columns under Seismic Sequences. Buildings. 2022; 12(3):275. https://doi.org/10.3390/buildings12030275

Chicago/Turabian Style

Fujii, Kenji. 2022. "Peak and Cumulative Response of Reinforced Concrete Frames with Steel Damper Columns under Seismic Sequences" Buildings 12, no. 3: 275. https://doi.org/10.3390/buildings12030275

APA Style

Fujii, K. (2022). Peak and Cumulative Response of Reinforced Concrete Frames with Steel Damper Columns under Seismic Sequences. Buildings, 12(3), 275. https://doi.org/10.3390/buildings12030275

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