Next Article in Journal
Water Resilience by Centipedegrass Green Roof: A Case Study
Next Article in Special Issue
Seismic Behaviour of EC8-Compliant Moment Resisting and Concentrically Braced Frames
Previous Article in Journal
A BIM-based PSS Approach for the Management of Maintenance Operations of Building Equipment
Previous Article in Special Issue
Proposed Simplified Approach for the Seismic Analysis of Multi-Storey Moment Resisting Framed Buildings Incorporating Friction Sliders
 
 
Correction published on 15 November 2019, see Buildings 2019, 9(11), 234.
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Pushover-Based Seismic Capacity Evaluation of Uto City Hall Damaged by the 2016 Kumamoto Earthquake

Department of Architecture, Faculty of Creative Engineering, Chiba Institute of Technology, Chiba 275-0016, Japan
Buildings 2019, 9(6), 140; https://doi.org/10.3390/buildings9060140
Submission received: 18 May 2019 / Revised: 31 May 2019 / Accepted: 3 June 2019 / Published: 5 June 2019

Abstract

:
The seismic capacity of the main Uto City Hall building, which was severely damaged by the 2016 Kumamoto Earthquake, was evaluated by the proposed pushover-based procedure. In this procedure, the seismic capacity index of the building is defined as the maximum scaling factor of the seismic input, for which the local responses do not exceed their limit values. From the pushover analysis result, the displacement limit of the equivalent single-degree-of-freedom model was determined. Then, the seismic capacity index was evaluated using an equivalent linearization technique. The evaluated index was re-evaluated by considering the bidirectional excitation. The pushover analysis result revealed that the torsional response is significant in the nonlinear behavior of this building. The evaluated seismic capacity implied that some structural damages, including the yielding of the beam-column joint, may have occurred during the first earthquake on 14 April 2016.

1. Introduction

1.1. Background

Buildings with highly irregular plans and/or elevations tend to be vulnerable to earthquakes. In the 2016 Kumamoto Earthquake, the main building of Uto City Hall was severely damaged [1]. This building is a five-story reinforced concrete (RC) irregular building, which was constructed in 1965. Given that the buildings around the city hall did not incur severe damages, it can be said that the main building had low seismic performance owing to structural irregularities and old construction methods based on an inadequate seismic design code. The author had carried out the preliminary seismic evaluation of this building based on the structural drawing [2]. The preliminary evaluation results had revealed that the seismic capacity of this building was insufficient to survive severe earthquakes. However, the evaluated results could not explain the damage of the building observed in upper stories. One of the reasons is that the effect of torsion is not considered in the simplified method applied for the preliminary evaluation.
Nowadays, the available simplified nonlinear analysis procedures combine the nonlinear static (pushover) analysis of a multi-degree-of-freedom (MDOF) model with the response spectrum analysis of an equivalent single-degree-of-freedom (SDOF) model [3,4,5]. These procedures have been widely implemented in seismic design codes and seismic design schemes [6,7,8,9], and work well for buildings that oscillate predominantly with a single mode. Although these simplified procedures are less accurate than the more sophisticated nonlinear time-history analysis, they provide us with basic information about the building under investigation. Thereby, they assist designers and analysts in understanding the nonlinear seismic behavior of buildings. In recent years, the author and other researchers have worked to extend these simplified procedures and improve the seismic performance estimates of buildings with plan and/or elevation irregularities [10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29].
These extended simplified procedures have been verified by comparison to nonlinear dynamic analysis results. However, it is challenging to verify them using existing irregular buildings that were damaged during an actual seismic event. Therefore, verifying these simplified procedures by considering the main Uto City Hall building is an important benchmark of their ability.

1.2. Brief Review of Related Studies

The author has determined that there are four possible approaches toward predicting the peak response of an asymmetric building with consideration given to the torsional effects. The first one is the extended N2 method [10,11,12,13,14]; the second one is a modal pushover analysis [15,16,17,18,19,20,21,22]; the third one consists of combining two pushover analyses and the envelope of the results, as proposed by Bosco et al. [23,24]; the fourth one consists of combining the analyses of two independently equivalent SDOF models (representing the first and second modes) with the envelope of four pushover analyses (including the effect of bidirectional excitation) [25,26,27,28,29].
The first approach, namely, the extended N2 method, is an extended version of the simplified procedure proposed by Fajfar and Fischinger [4]. This approach involves the estimation of the peak response for each frame in the pushover analysis results multiplied by a correction factor, which is defined using the linear dynamic analysis and pushover analysis results [10,11,12,13].
The second approach, namely, modal pushover analysis (MPA), was proposed by Chopra and Goel for regular buildings and considers the effect of higher modes. Subsequently, this approach was extended to asymmetric buildings [15]. Then, it was extended further by Reyes and Chopra [16,17], who considered the effect of bidirectional excitation. Manoukas et al. [18,19] proposed the concept of an equivalent SDOF model, which differed from that proposed by Chopra. Manoukas et al. considered the effect of bidirectional excitation in the formulation of the equivalent SDOF model. Belejo and Bento [20] applied the improved pushover analysis (IMPA) (a modified version of MPA that considers mode shape changes in the inelastic regime) to buildings with three and nine stories. Another interesting variant of MPA based on energy-based pushover analysis has been proposed by Soleimani et al. [21].
As proposed by Bosco et al. [23], the third approach estimates the peak response of the frames on the stiff side and the peak response of the frames on the flexible side by enveloping the results of two pushover analyses. Bosco et al. [24] investigated the applicability of this procedure to multistory building models with the same geometry on each floor. In this procedure, corrective eccentricity is a key parameter in the pushover analyses.
The fourth approach has been proposed by the author [25,26,27,28,29]. In this paper, this simplified procedure is called mode-adaptive bidirectional pushover analysis (MABPA), and it involves the evaluation of the first and second mode peak responses independently from the equivalent SDOF models formulated with consideration given to the principal direction of the first modal response. The prediction of the peak response at each frame is based on a set of pushover analyses that consider the combination of the two modal responses. Although the latest version of MABPA [26] can successfully predict the peak response of asymmetric buildings, it requires two critical assumptions. The first assumption is that the building oscillates predominantly in a single mode and in each set of orthogonal directions. The second assumption is that the principal directions of the first and second modal responses are almost orthogonal. The limitation of the current version of MABPA has been discussed in a previous report [29].
These studies mentioned above were mainly focused on the prediction of the peak response (mainly displacement) under given ground motions (mostly by the response spectrum). Some of these procedures may be used for the prediction of the seismic intensity, or scaling factor, corresponding to a certain (given) displacement limit. Dolšek and Fajfar [14] proposed a simplified probabilistic seismic performance assessment for plan-asymmetric buildings. In their study, the ground motion intensity at the failure point was evaluated using the extended N2 method, so-called the incremental N2 method. Also, as an application of MPA, a scaling procedure for the nonlinear dynamic analysis of asymmetric buildings with consideration given to the bidirectional excitation was proposed by Reyes et al. [22]. In their procedure, the ground acceleration was scaled such that the peak roof displacement was equal to the target displacement. If the given seismic intensity is too large for the building considered, its peak response would be too far from its limit value: sometimes its peak response may not be obtained because the analysis would be unstable. In such case, the evaluation of the scaling factor of the given seismic input, wherein the local responses do not exceed their limit value, would be helpful for the understanding of the seismic behavior of the building considered.

1.3. Objectives

In this study, the seismic capacity of the main Uto City Hall building was evaluated according to MABPA, which has been previously proposed by the author. The seismic capacity index, which is defined in this paper as the seismic scaling factor of the given seismic input, is evaluated considering bidirectional excitation. The verification of MABPA is made by comparing the observed damages of the main Uto City Hall building and the analysis results, especially concerning the effect of torsion to the structural damage.
The rest of the paper is organized as follows. In Section 2, the seismic capacity index discussed in this paper is defined. Subsequently, a description of the seismic capacity index is provided with consideration given to the unidirectional and bidirectional ground motion. In Section 3, basic information is provided with regard to the main Uto City Hall building, its model, and ground motions. In Section 4, the evaluation procedure of the seismic capacity index is verified by using nonlinear dynamic analyses. Finally, in Section 5, the seismic capacity of this building is evaluated, considering the variation of modeling about strength degradation behavior of walls and basement spring. Then the evaluated seismic capacity curves are compared to the response spectrum of the first and second earthquakes, occurred on the 14 and 16 April, respectively, and discuss the performance of this building against two earthquakes.

2. Seismic Capacity Evaluation Procedure

2.1. Definition of Seismic Capacity Index

A set of orthogonal axes (ξ-ζ) in the X-Y plane is considered, where the ξ- and ζ-axis are the major and minor axes, respectively, of the horizontal ground motion components. Another set of orthogonal axes, namely U-V in the X-Y plane, is considered, with the U- axis being the principal axis of the first modal response [25]. It is assumed that the spectra of the two horizontal components are identical, as expressed by Equation (1).
S A ξ p ( T , h ) = S A ζ p ( T , h ) = S A U p ( T , h ) = S A V p ( T , h ) = S A p ( T , h ) .
In Equation (1) pSA(T, h) is the pseudo acceleration spectrum, and T and h are the natural period and damping ratio, respectively, of the elastic SDOF oscillators. Note that this assumption, namely, the identical-component assumption, is the same as discussed by López and Torres [30] with regard to the elastic spectrum analysis. In addition, López et al. [31] reported that the ratio of the spectra of the horizontal minor and major components varies between 0.63 and 0.81. Therefore, the identical-component assumption is expected to provide conservative prediction results.
The seismic capacity index (CI) is defined in Figure 1. The building model is assumed to be converted to the equivalent SDOF model representing the first mode response. The equivalent acceleration A1U*-D1U* relationship, which is shown in this figure, will be termed hereafter as the capacity curve. The displacement limit (D1U*limit) is defined as the equivalent displacement, wherein one of the local responses, such as the story or column drift and the plastic hinge rotation angle of the beam, reaches the limit value. The pseudo-acceleration spectrum-displacement spectrum (pSA-SD) relationship of the design (or code-specified) ground motion is considered, and is also shown in this figure. The factored pSA-SD relationship considering the equivalent damping heq will be termed hereafter as the demand curve. The intersection point of the capacity and the demand curve represents the predicted peak response. The seismic capacity index is defined as the maximum scaling factor, wherein the local responses do not exceed their limit value. In this figure, the capacity index CI is equal to the scaling factor producing the peak equivalent displacement equal to D1U*limit.
In this paper, the seismic capacity index that considers the unidirectional ground motion acting in the principal direction of the first modal response is denoted as CI, uni, whereas that considering the bidirectional ground motion acting simultaneously is denoted as CI, bi. As described below, CI, uni is evaluated only on the basis of the first modal response, whereas CI, bi is evaluated by using MABPA and by considering the first and second modal responses.
The basic idea of the seismic capacity index CI is identical: to the ratio (%NBS) defined in the New Zealand (NZ) seismic assessment guideline for existing buildings [32,33,34]; to the ratio of the ultimate capacity of a building as a whole or of an individual member; to the ultimate limit state shaking demand for a similar new building on the same site. The index CI, uni is identical to the ratio (%NBS) obtained by the nonlinear static pushover analysis by considering the single mode, or obtained by the simple lateral mechanism analysis (SLaMA), the simplified nonlinear static analysis [33,34]. Veccho et al. show an application of the SLaMA method to an existing RC building which has severely damaged during the 2011 Christchurch earthquake [35]. The advantage of the proposed procedure is in considering the bidirectional excitation, which may be more significant for a building with an irregular structural plan.
Note that the presented procedure can be applicable to a building which oscillates predominantly in the first mode under U-directional (unidirectional) excitation, as discussed in the previous study [29]. This is because in this procedure the behavior of the building can be represented by the equivalent SDOF model of the first mode, even though the effect of the second mode would be included in the evaluation of CI, bi. Therefore, the effective first modal mass ratio of the building should be checked for the applicability of this procedure.

2.2. Evaluation of the Seismic Capacity Index Considering the Unidirectional Ground Motion

2.2.1. Step 1: Pushover Analysis of Building Model (First Mode)

The nonlinear equivalent acceleration-equivalent displacement (A1U*-D1U*) relationship of the equivalent SDOF model was determined based on the pushover analysis of an N-story irregular building model considering the change of the first mode’s shape at each nonlinear stage.
In this study, the DB-MAP analysis was adopted [26]. The fundamental assumptions of the current version are as follows:
  • The envelope curve of the force-deformation relationship for each nonlinear spring of all members is symmetric in the positive and negative loading direction.
  • The equivalent stiffness of each nonlinear spring can be defined by its secant stiffness at the peak deformation that was previously derived in the calculation.
  • The first mode shape at each loading step ( Γ 1 U n φ 1 n ) can be determined according to the equivalent stiffness.
  • The displacement vector at each loading step (nd) imposed on the model is similar to the first mode shape obtained in assumptions 2 and 3.
  • In the case where unloading occurs, the equivalent stiffness obtained in assumption 2 is used as the unloading stiffness for all nonlinear springs.
Note that assumption 5 was adopted from previous studies [26,27,28,29] to apply this DB-MAP to a building model with non-ductile members.
Let nfR be the restoring force vector of the building model at each loading step n obtained from the pushover analysis. The equivalent displacement and acceleration at step n (namely, nD1U* and nA1U*, respectively) are determined from Equations (2) and (3), respectively, assuming that nd is proportional to the first mode vector Γ 1 U n φ 1 n at each loading step:
D 1 U n * = Γ 1 U n φ 1 n T M n d M 1 U n * = j ( m j   x j n 2 + m j   y j n 2 + I j   θ j n 2 ) ( j m j   x j n ) 2 + ( j m j   y j n ) 2 ,
A 1 U n * = Γ 1 U n φ 1 n T f R n M 1 U n * = j ( f R X j n   x j n + f R Y j n   y j n + f M Z j n   θ j n ) ( j m j   x j n ) 2 + ( j m j   y j n ) 2 ,
M 1 U n * = Γ 1 U n 2 φ 1 n T M φ 1 n = ( j m j   x j n ) 2 + ( j m j   y j n ) 2 j ( m j   x j n 2 + m j   y j n 2 + I j   θ j n 2 ) ,
Γ 1 U n = φ 1 n T M α U n φ 1 n T M φ 1 n = ( j m j   x j n ) 2 + ( j m j   y j n ) 2 j ( m j   x j n 2 + m j   y j n 2 + I j   θ j n 2 ) ,
M = [ M 0 0 0 0 M 0 0 0 0 I 0 ] , M 0 = [ m 1 0 0 m N ] , I 0 = [ I 1 0 0 I N ] ,
{ d n = { x 1 n x N n y 1 n y N n θ 1 n θ N n } T f R n = { f R X 1 n f R X N n f R Y 1 n f R Y N n f M Z 1 n f M Z N n } T ,
φ 1 n = { ϕ X 11 n ϕ X N 1 n ϕ Y 11 n ϕ Y N 1 n ϕ Θ 11 n ϕ Θ N 1 n } T ,
α U n = { cos ψ 1 n cos ψ 1 n sin ψ 1 n sin ψ 1 n 0 0 } T ,
cos ψ 1 n = j m j   ϕ X j 1 n ( j m j   ϕ X j 1 n ) 2 + ( j m j   ϕ Y j 1 n ) 2 = j m j   x j n ( j m j   x j n ) 2 + ( j m j   y j n ) 2 ,
sin ψ 1 n = j m j   ϕ Y j 1 n ( j m j   ϕ X j 1 n ) 2 + ( j m j   ϕ Y j 1 n ) 2 = j m j   y j n ( j m j   x j n ) 2 + ( j m j   y j n ) 2 .
In Equations (2)–(11), mj and Ij are the mass and moment of inertia, respectively, of the jth floor; nM1U* is the equivalent (effective) first modal mass with respect to the U-axis at each loading step; and nψ1 is the angle of incidence of the principal axis (U-axis) at each loading step. Figure 2 shows the comparisons of parameters between an N-story irregular building and the equivalent SDOF model.
Pushover analysis was carried out until one of the local responses (e.g., the drift angle of the columns and the plastic hinge rotation angle of the beam) reached its limit value. The displacement limit (D1U*limit) was obtained as the equivalent displacement corresponding to the loading step, where one of the local responses exceeded its limit value.

2.2.2. Step 2: Calculation of the Scaling Factor Corresponding to Each Pushover Step

The scaling factor corresponding to each pushover step n (nλ1) is calculated from the given design (or code-specified) response spectrum and the capacity curve obtained in step 1 using the equivalent linearization technique [36]. The equivalent period nT1eq and equivalent damping ratio nh1eq of the equivalent SDOF model at each loading step n are determined as follows:
T 1 e q n = 2 π / ω 1 e q n = 2 π D 1 U n * / A 1 U n * ,
h 1 e q n = k ( h e q k n   W e k n ) / k W e k n .
In Equation (13), nheqk and nWek are the equivalent damping ratio and potential energy, respectively, of the kth nonlinear spring. The equivalent damping ratio of the flexural spring is calculated from Equation (14) with consideration given to the initial damping ratio h0, and the kth ductility of the flexural spring (nμk).
h e q k n = { h 0 : μ k n < 1 h 0 + 0.25 ( 1 1 / μ k n ) : μ k n > 1 .
The equivalent damping ratio of the other spring is set to h0 throughout the analysis. In this study, h0 was assumed to be 0.05.
The scaling factor corresponding to each pushover step n (nλ1) is calculated as follows:
λ 1 n = D 1 U n * F ( h 1 e q n ) S D ( T 1 e q n , 0.05 ) = 1 + 10 h 1 e q n 1.5 D 1 U n * S D ( T 1 e q n , 0.05 ) .

2.2.3. Step 3: Evaluation of the Seismic Capacity Index

The seismic capacity index that considers the unidirectional ground motion (CI, uni) was taken from the result obtained by step 2. Additionally, CI, uni is the maximum of nλ1 in the range from D1U* = 0 to D1U*limit, as follows:
C I , u n i = max { λ 1 n ( D 1 U n * ) ; 0 D 1 U n * D 1 U * l i m i t } .

2.3. Evaluation of the Seismic Capacity Index Considering the Bidirectional Ground Motion

2.3.1. Steps 1 to 3: Evaluation of the Seismic Capacity Index Considering the Unidirectional Ground Motion

The seismic capacity index that considers the unidirectional ground motion (CI, uni) was evaluated according to the steps described in Section 2.2.

2.3.2. Step 4: Prediction of the Peak Response of the Equivalent SDOF Model (first mode)

The peak response of the equivalent SDOF model representing the first mode (the peak equivalent displacement under unidirectional excitation D1U*uni, and the equivalent acceleration A1U*uni corresponding to D1U*uni) for scaling, given the ground motion by multiplying CI, uni, was obtained using the equivalent linearization technique. Note that D1U*uni may be smaller than D1U*limit. In this case, CI, uni is not a value at D1U*limit.

2.3.3. Step 5: Pushover Analysis of the Building Model (second mode)

From the results of steps 1 and 4, the first mode vector and the principal axis angle of the first modal response corresponding to D1U*uni, namely, Γ 1 U i e φ 1 i e and ψ1ie, respectively, are obtained. The second mode vector ( Γ 2 V i e φ 2 i e ) is determined from Equations (17) and (18) in terms of Γ 1 U i e φ 1 i e and the second mode vector at the elastic range ( Γ 2 V e φ 2 e ) with consideration given to the orthogonality of the mode vectors, as follows:
Γ 2 V i e = φ 2 i e T M α V i e φ 2 i e T M φ 2 i e ,
φ 2 i e = φ 2 e φ 2 i e T M φ 1 i e φ 1 i e T M φ 1 i e φ 1 i e ,
α V i e = { sin ψ 1 i e sin ψ 1 i e cos ψ 1 i e cos ψ 1 i e 0 0 } T .
Next, an additional pushover analysis for an MDOF model was carried out to obtain the force-displacement relationship representing the second mode response by applying the invariant force distribution p2, which is determined as follows:
p 2 = M ( Γ 2 V i e φ 2 i e ) .
The equivalent displacement nD2V* and the equivalent acceleration nA2V* of the equivalent SDOF model representing the second modal response at each loading step n are determined as follows:
D 2 V n * = Γ 2 V i e φ 2 i e M n T d / M 2 V i e * ,
A 2 V n * = Γ 2 V i e φ 2 i e T f R n / M 2 V i e * ,
M 2 V i e * = Γ 2 V i e 2 φ 2 i e T M φ 2 i e .
In Equation (23), M2Vie* is the equivalent second modal mass with respect to the V-axis, and is determined in terms of Γ 2 V i e φ 2 i e .

2.3.4. Step 6: Peak Response Prediction for the Equivalent SDOF Model (second mode)

The peak response of the equivalent SDOF model representing the second mode (D2V*uni, which is the peak equivalent displacement under unidirectional excitation, and A2V*uni, which is the equivalent acceleration corresponding to D2V*uni) for scaling, given the ground motion by multiplying CI, uni, was obtained using the equivalent linearization technique.
In the current version of MABPA [26,27,28,29], the principal axes of the first and second modes are assumed to be almost orthogonal. Moreover, the principal axis of the second modal response is almost identical to the V-axis. As has been discussed in the literature [28,29], this assumption is valid in the case where the torsional indices of the first and second modes, namely, Rρ1 and Rρ2, respectively (defined in Equation (24)) are smaller than 1.
R ρ i = j I j ϕ Θ j i 2 / j m j ( ϕ X j i 2 + ϕ Y j i 2 ) .
However, this assumption may not be valid if Rρ1 or Rρ2 is larger than 1, as discussed in Fujii [29]. Therefore, in this study, the response spectrum used to predict the peak response (D2V*uni and A2V*uni) was modified by dividing |sinΔψ12|, where Δψ12 (= ψ1ψ2) is the angle between the principal axes of the first and second modal responses. Details regarding this modification can be found in Appendix A.

2.3.5. Step 7: Pushover Analysis of the Building Model Considering the Bidirectional Seismic Input

In this step, the equivalent displacement limit of the first and the second mode considering the bidirectional seismic inputs D1U*bi and D2V*bi, respectively, are evaluated according to four pushover analyses considering the bidirectional seismic input.
First, from the results obtained in steps 4 and 6, the combined force distributions, namely pU+, pU, pV+, and pV, are determined as follows:
{ p U ± = M ( Γ 1 U i e φ 1 i e A 1 U * u n i ± 0.5 Γ 2 V i e φ 2 i e A 2 V * u n i ) p V ± = M ( ± 0.5 Γ 1 U i e φ 1 i e A 1 U * u n i + Γ 2 V i e φ 2 i e A 2 V * u n i ) .
Next, pushover analysis (termed as pushover 1U and 2U, respectively) is carried out using pU+ or pU until a) the equivalent displacement nDU* calculated from Equation (26) reaches D1U*uni, as obtained in step 4, or b) one of the local responses reaches its limit value, as follows:
D U n * = Γ 1 U i e φ 1 i e T M n d / M 1 U i e * ,
M 1 U i e * = Γ 1 U i e 2 φ 1 i e T M φ 1 i e .
In Equation (27), M1Uie* is the equivalent first modal mass with respect to the U-axis, and is determined in terms of Γ 1 U i e φ 1 i e . The displacement limit of the first mode considering the bidirectional ground motion (D1U*bi) is the minimum value of a) the equivalent displacement nDU*, when one of the local responses reaches its limit value in pushovers 1U and 2U, and b) the peak equivalent displacement D1U*uni obtained in step 4.
Similarly, pushover analysis (termed as pushover 1V and 2V, respectively) is performed using pV+ or pV until a) the equivalent displacement nDV* calculated from Equation (28) reaches D2V*uni, as obtained in step 6, and b) one of the local responses reaches its limit value:
D V n * = Γ 2 V i e φ 2 i e T M n d / M 2 V i e * .
The displacement limit of the second mode that considers the bidirectional ground motion (D2V*bi) is the minimum value of a) the equivalent displacement nDV*, when one of the local responses reaches its limit value in pushovers 1V and 2V, and b) the peak equivalent displacement D2V*uni obtained in step 6.

2.3.6. Step 8: Evaluation of the Seismic Capacity Index Considering the Bidirectional Ground Motion

From the results obtained in steps 2 and 7, the scaling factor obtained from the first modal response considering the bidirectional seismic input λ1bi is obtained as the maximum nλ1 in the range from D1U* = 0 to D1U*bi, as follows:
λ 1 b i = max { λ 1 n ( D 1 U n * ) ; 0 D 1 U n * D 1 U * b i } .
Similarly, the scaling factor obtained from the second modal response considering the bidirectional seismic input λ2bi is obtained as the maximum nλ2 in the range from D2V* = 0 to D2V*bi, as follows:
λ 2 b i = max { λ 2 n ( D 2 V n * ) ; 0 D 2 V n * D 2 V * b i } ,
λ 2 n = 1 | sin Δ ψ 12 | 1 + 10 h 2 e q n 1.5 D 2 V n * S D ( T 2 e q n , 0.05 ) .
The seismic capacity index considering the bidirectional ground motion (CI, bi) is obtained as the minimum of λ1bi and λ2bi, and does not exceed CI, uni.

3. Building and Ground Motion Data

3.1. Basic Information

Figure 3 shows the structural plan of the main Uto City Hall building. The structure of this building can be divided into two structural blocks, namely, the office block and the stair block. The two blocks are connected only by a concrete slab with a thickness of 110 mm. Additionally, as shown in the figure, all of the structural walls are concentrated to the stair block, whereas, in the office block the concrete columns are the only vertical members resisting the lateral loads. Moreover, it should be noted that in the office block, not all frames are oriented in the X- or Y-directions. Specifically, frames A1-A3 lie on the axis rotated 45° counterclockwise from the X-axis (denoted as A-axis), and frames B1-B3 are orthogonal to frames A1 (denoted as B-axis).
Figure 4 shows the simplified structural elevation of frames Y4, B1, and B2. As shown in Figure 4b, the beam for the second-floor slab is doubled in frame B1. The reason for this is that, for levels Z1’ and Z1, the floor level is different in the second floor, as shown in Figure 4.
Figure 5 shows the Uto City Hall after the 2016 Kumamoto Earthquake. As can be seen in the figure, the structural damage observed in frame B1 was more severe than that incurred by frame A3. That is, frame B1 partially collapsed in the fourth story.
Figure 6 shows the failure of the beam-column joint observed in frame B1. As is shown in this figure, the part above the column A2B1 in the third story had slipped down due to the joint failure. Details regarding the structural damages and other information can be found in the literature [1]. Some of the other photos of the structural damages can be found in Appendix B.
Figure 7 shows the sections of column A1B1, A2B1, and A2B2. In this building, the deformed bar was used for the longitudinal reinforcement of the beams and the columns, while the round bar was used for the shear reinforcement.
As shown in this figure, the sectional area of the column is drastically reduced from lower stories to upper stories. This is very common in the reinforced concrete buildings constructed before 1981 because at that time the design seismic force in upper stories is smaller than that in current seismic design code of Japan.
In the previous study [2], the author has made the preliminary seismic evaluation of this building according to the simplified procedure by Shiga in the 1970s, which is based the wall-area index and the average shear stress in walls and columns [37]. The result of the preliminary seismic evaluation is summarized as follows.
Figure 8 shows the evaluation results shown in reference [2]. In this figure, AWi (unit: mm2) and ACi (unit: mm2) are the sum of the sectional area of walls and columns in ith story, respectively, Ai is the coefficient of the design seismic story shear force used in the current seismic design code of Japan, Afj (unit: m2) is the area of jth floor, and τavei (unit: N/mm2) is the average shear stress in walls and columns in ith story, assuming the weight per unit floor area of the building is 10kN/m2 and base shear coefficient is 1.0. The zone A is the area corresponds to the most of buildings were heavily damaged while the zone C is the area corresponds to the most of buildings were not damaged or only slightly damaged, in 1968 Tokachi-oki Earthquake [37]. Note that in each case, the evaluation is made for the X- and Y-directions, as shown in Figure 3. In Figure 8, the border curve between zones A and B shown in red is determined by assuming the ultimate shear stress in columns equals 1.2N/mm2 while that in walls equals 3.3 N/mm2. Besides, the horizontal border line between zones B and C is determined as the line the average shear stress τavei equals to 1.2 N/mm2, and the vertical border line is determined as the line the wall-area index AWi / (Ai ΣAfj) equals to 30 × 102 mm2/m2. The curves shown in black is the contour of equal column-area index ACi / (Ai ΣAfj).
Figure 8a shows the result of the case in which the building is assumed to behave as a unite building. In this figure, the plots of Y-direction in the first and second story are within zone A, while the plots of upper stories are in zone B or C. Therefore, it may be concluded that the seismic capacity of this building is insufficient to survive strong earthquakes, if this building behaved as a unite single building. However, this result cannot explain the fact that most of the damage in this building is in the upper stories; the most critical of which is the second story in the Y-direction, where the average shear stress is the largest.
Figure 8b shows the results of the case in which the stair and office blocks are assumed to behave independently. In this figure, the plots of office block in all stories are within zone A, while the plots of stair block in all stories are within zone C. Therefore, it may be concluded that the seismic capacity of the office block is insufficient while that of stair block is sufficient, if these two blocks of this building behaved as two independent buildings. However, in the damage observation of this building [1], it is found that the walls in the fifth story of frame Y5 (stair block) are severely damaged. The results shown in Figure 8b cannot explain this damage. Therefore, the assumption that the stair and office blocks behave independently has seemed invalid, even though the floor slab at the border of two blocks in the fourth and fifth floors were severely damaged [1].
From the preliminary evaluation, it may be concluded that the seismic capacity of this building is insufficient to survive strong earthquakes. However, neither results can explain the damage of this building observed. The reasons why this simplified evaluation method fails to explain the observed damage are (i) the lateral force distribution coefficient, Ai, is smaller in upper stories because the Ai coefficient cannot reflect the drastic reduction of the sectional area in the upper stories, and (ii) the effect of torsion is not considered in this simplified method. Other reasons would be (iii) the failure of the beams and beam-column joints are not considered in this simplified method.
In the following sections, the discussions are mainly focused on the difference in the damage observed in frame B1 and A3. Because this building had an irregular plan, the effect of torsion would be significant. The author thinks the difference of the damage in these frames may be caused due to the torsional response. Therefore, the three-dimensional nonlinear static and dynamic analyses using spatial frame model were carried out in which the effect of torsion can be considered.

3.2. Analysis Model

The main Uto City Hall building is modelled as a three-dimensional spatial frame, wherein the floor diaphragms are assumed to be rigid in their own planes without out-of-plane stiffness. Moreover, in this analysis model, the slabs connecting the office and stair blocks were considered as rigid. The compressive strength of the concrete was assumed 24 N/mm2, while the yield strength of the deformed bar and the round bar were assumed 343 N/mm2 and 294 N/mm2, respectively.
Figure 9 shows the structural model of frames B1 and Y5.
As shown in Figure 9a, the double beam in the second-floor level, which is shown in Figure 4b (beams in levels Z1 and Z1’), is modeled as a single beam. The properties of the second-floor beams were determined by considering both beam sections. A one-component model with one non-linear flexural spring at each end and one linear shear spring at the middle of the line element was used for all beams and columns. The bidirectional interaction of the columns’ bending moment is not considered because the rigid zone length of some columns is different in each orthogonal direction (e.g., the second story column A2B1(X2AY2) in frames A1 and B2, and in frames X2A and Y2). The axes for determining the bending moment and rotation were taken as the A- and B-axis for the columns belonging to the frames oriented in the A- and B-axis, whereas the rest were taken as the X- and Y-axis. To determine the flexibility of the springs, we assume an antisymmetric curvature distribution. For the beams (except for the beam within the concrete wall), a rigid zone length is assumed as half the depth of the intersected column minus one-fourth of the depth of the considered beam. For the beams within the concrete wall, the rigid zone length was not considered.
Additionally, the stiffness of the beams within the wall is 100 times higher than that calculated by considering a rectangular beam. For the columns, the rigid zone length is assumed as half the depth of the intersected beam minus one-fourth of the depth of the considered column. The shear behavior of the beams and columns is assumed to be linear elastic. The concrete walls were modeled as two diagonal braces, as shown in Figure 9b, assuming that the shear behavior is predominant. The axial behavior of all vertical members is assumed to be linear elastic.
Figure 10 shows the force-deformation relationship of the nonlinear flexural spring.
The envelopes are assumed to be symmetric in the positive and negative loading directions. The yield moment My of each member is calculated according to the AIJ standard [38], whereas the crack moment Mc is assumed to be one-third of My. For simplicity, only the axial force attributed to the vertical load was considered in the calculation of each column’s yield moment. Moreover, as has been described in a previous report, various beam-column joint failures were observed [1]. Additionally, the yielding moment of the corresponding beam ends was reduced according to the AIJ standard [38] to consider the effect of the beam-column joint yield at the ultimate strength. The range of the reduction factor was 0.837 to 1.000. In Figure 10a, the secant stiffness degradation ratio (αy) of the flexural spring at the yield point is calculated according to the equation of Sugano’s equation and Koreishi [39]. The tangent stiffness degradation ratio (α2) of the flexural spring beyond the yielding point is assumed to be 0.001 for all beams and columns.
The Muto hysteresis model [40] with one modification was used to model the flexural springs, as shown in Figure 10b. Specifically, the unloading stiffness after yielding was proportionally decreased to μ−0.5 (μ is the ductility of the flexural spring) to represent the degradation of the unloading stiffness after the yielding of the RC members, in the same way as in Otani’s model [41].
Figure 11 shows the simplified concrete wall model. The properties of the diagonal braces were determined to represent the shear force-deformation behavior. The property of the vertical elements (columns at each border and within the walls) was determined to represent the bending and axial stiffness of the concrete wall, and their behavior is assumed to be linear elastic. Figure 11b,c show the force-deformation relationship of the diagonal braces. The envelope shown in Figure 11b is assumed to be symmetric in the positive and negative loading directions. The stiffness and strength of the diagonal braces were determined with consideration given to the existing opening width. The ultimate shear strength of the concrete wall (QsuW) is calculated according to the AIJ standard [38] and reduced by considering the ratio of the sum of the opening width divided by the wall length (lW). To determine the axial deformation corresponding to the point of maximum shear strength, we assume the shear strain γsu to be 0.004. Similarly, to determine the shear strength’s point of loss, we assume the shear strain γu to be 0.010. The origin-oriented model (Figure 11c) was used to model the non-linear behavior of the diagonal braces.
The properties of the vertical and rotational springs at each basement are calculated on the bases of the axial stiffness of the piles according to the method presented in FEMA 273 [7]. The vertical and rotational springs are assumed to be linear elastic.
Second-order effects (e.g., the P-Δ effect) were not considered. The damping matrix is assumed to be proportional to the tangent stiffness matrix, with 5% of the first mode’s critical damping. Note that, when its tangent stiffness is negative, non-zero dummy positive value (10−9 times smaller than the initial stiffness) is used for the calculation of the damping matrix to avoid the negative damping.
The building model developed in this study may not sufficiently represent the behavior of the actual building during the 2016 Kumamoto Earthquake. In particular, the assumption that the behavior of all nonlinear springs, including the diagonal braces representing the behavior of the concrete wall, can be determined as strictly symmetric in the positive and negative loading directions, including the axial deformation of the vertical members (all of them are assumed to have linear elastic behavior) is the most critical consideration. Owing to the nature of the concrete’s behavior, the symmetric behavior of the members in the positive and negative loading directions may not be true. However, in this study, this modeling scheme was applied to maintain the stability of the DB-MAP analysis as much as possible after the strength of the wall degraded. For the same reason, the vertical and rotational springs at the basement are assumed to have a linear elastic behavior. Influence of the strength degradation behavior of the concrete walls and stiffness of the vertical and rotational springs at the basement to the whole behavior of the building would be discussed in Section 5.
Figure 12 shows the natural modes of the building model in the elastic range. Here, Tie is the ith natural period in the elastic range (i = 1–3), mie* is the equivalent (effective) modal mass ratio of the ith mode with respect to its principal direction in the elastic range, ψie is the incidence angle of the principal direction of the ith modal response in the elastic range, and Rρie is the torsional index of the ith mode in the elastic range.
As shown in the figure, the first mode is almost purely translational (Rρ1e = 0.008 << 1), the second mode is predominantly torsional (Rρ2e = 1.438 > 1), and the third mode is predominantly translational (Rρ3e = 0.670 < 1). Therefore, this building cannot be classified as torsionally stiff (TS). However, the equivalent first modal mass ratio is larger than 0.5 (m1e* = 0.759 > 0.5). Therefore, this building may oscillate predominantly in the first mode from the U-directional (unidirectional) excitation. Additionally, the angles between the principal directions of the first two modes are close to 90° such that Δψ12 = 86.7°− (−3.6°) = 90.3°.

3.3. Ground Motion Data

In this study, the seismic excitation was bidirectional in the X-Y plane, and ten sets of artificial ground motions were generated. Additionally, the response spectra of the major and minor components were assumed to be identical. The target elastic spectra of the major and minor components with 5% critical damping (pS(T, 0.05) and pS(T, 0.05), respectively) were determined using the current Building Standard Law of Japan with regard to extremely rare earthquake events and with consideration given to the regional seismic event factor in Uto City (Z = 0.8). Moreover, the type-2 soil (normal) is calculated using Equation (32), where T represents the natural period of the SDOF model, as follows:
S A ξ p ( T , 0.05 ) = S A ζ p ( T , 0.05 ) = { 0.8 ( 4.8 + 45 T ) m / s 2 : T 0.16   s 9.60 m / s 2 : 0.16   s < T 0.864   s 9.60 ( 0.864 / T ) m / s 2 : T > 0.864   s .
The phase angle is given by the uniform random value and the Jenning-type envelope function e(t) recommended by the Building Center of Japan [42,43], as follows:
e ( t ) = { ( t / 5 ) : 0   s < t 5   s 1 : 5   s < t 35   s exp { 0.027 ( t 35 ) } : 35   s < t 120   s .
Figure 13 shows the elastic response spectra of the artificial ground motions with 5% critical damping. Note that the artificial ground motions considered in this study were generated independently, i.e., there is no correlation between each component. The correlation coefficients of all ten sets are close to zero, and the envelope functions of the two components are identical. Therefore, the two components can be considered to be independent of each other.
Figure 14 shows an example of accelerations and orbit of the artificial ground motion used in the time-history analysis shown in Section 4.

4. Validation of the Seismic Capacity Evaluation Procedure

This section presents the seismic capacity evaluation results, obtained by pushover analysis, and their validation. In this study, the displacement limit was determined based on the story drift at each column (vector value) in frames A1 to A3 (and also B1 to B3). The limit value of the story drift was assumed as 1/75 of its height. The reasons for this are that (i) no severe shear failure were observed in all columns, and (ii) the yielding of the beam-column joint, which was observed in frames A3 and B1, affects the collapse mechanism significantly when a larger drift occurs.

4.1. Pushover Analysis Results

Figure 15 shows the non-linear properties of the equivalent SDOF model representing the first mode response obtained from the pushover analysis results. As discussed above, the pushover analysis terminated when the story drift angle at column A1B1 in the second story exceeded 1/75. From Figure 15a, the displacement limit was obtained as D1U*limit = 8.420 × 10−2 m. The equivalent period T1eq at D1U*limit was 1.191 s, as shown in Figure 15b. The equivalent damping h1eq at D1U*limit was 6.115 × 10−2, as shown in Figure 15c.
Figure 16 shows the variations in the first mode parameters (the effective modal mass ratio m1U*, the incidence angle of the principal axis for the first modal response ψ1, and the torsional index of the first mode Rρ1) based on the results of the pushover analysis. As shown in Figure 16a, the effective modal mass ratio (m1U*) decreased as the equivalent displacement (D1U*) increased from m1U* = 0.759, in the elastic range, to m1U* = 0.5850 at D1U* = D1U*limit. The incidence angle of the principal axis of the first mode (ψ1) changed from 86.7° in the elastic range to 75.1° at D1U* = D1U*limit, as shown in Figure 16b. The torsional index of the first mode Rρ1 increased significantly from 0.008 to 0.516 at D1U* = D1U*limit, as shown in Figure 16c.
Figure 17 shows the roof displacement responses based on the results of the pushover analysis.
As shown in Figure 17a, the roof displacement at frame B1 was the largest of the three frames parallel to the A-axis throughout the response. Similarly, the displacement at frame A1 was the largest of the three frames parallel to the B-axis, as shown in Figure 17b. Therefore, in comparison with the responses of the three corner columns (A1B1, A3B1, and A3B3), the displacement of column A1B1 was the largest (Figure 17c).
Figure 18 shows the shape of the first natural mode at D1U* = D1U*limit. In this figure, significant torsion can be observed in the first mode. Here, frame B1 is the most critical (largest deformation observed), while frame A3 is less critical.

4.2. Evaluation of the Capacity Index

First, the capacity index that considers the unidirectional ground motion (CI, uni) was evaluated. Figure 19 shows the capacity index evaluation result. Figure 19a shows that the evaluated peak response for each scaling factor λ. From this figure, the relation between the scaling factor corresponding to each pushover step n (nλ1) and equivalent displacement nD1U* can be constructed as shown in Figure 19b. Therefore, the estimated value was the value at D1U* = D1U*limit = 8.420 × 10−2 m and CI, uni = 0.361.
Next, the capacity index that considers the bidirectional ground motion (CI, bi) was evaluated.
Figure 20 shows the peak response prediction for the first and second modal responses of the scaling factor λ = CI, uni = 0.361. Because the angle between the principal axes of the first and second modal responses (Δψ12) was not close to 90°, the demand curve for estimating the second mode was factored by 1/|sinΔψ12|. Thus, the calculated value of Δψ12 is 116.2° (ψ1 = 75.1° and ψ2 = −41.1°). Thereby, the value of | sinΔψ12| is 0.8973. Hence, the demand curve was factored by 1/0.8973 = 1.115. As shown in this figure, the estimated peak responses are D1U* = D1U*limit = 8.420 × 10−2 m and D2V*uni = 6.780 × 10−2 m, respectively.
Figure 21 shows the evaluation of the capacity index that considers the bidirectional ground motion. Figure 21a shows the plots of the roof displacement response relationships at column A1B1, as obtained from the results of steps 1 (first mode) and 7 (pushovers 1U and 2U using the invariant force vector PU+ and PU, respectively). As can be seen, the smallest equivalent displacement, where the story drift exceeded 1/75 is 7.54 × 10−2 m, and was obtained from the result of pushover 1U. Therefore, the displacement limit of the first mode considering the bidirectional ground motion D1U*bi is 7.54 × 10−2 m. From the results of pushovers 1V and 2V (not presented in this paper), the displacement limit of the second mode considering the bidirectional ground motion (D2V*bi) is equal to D2V*uni = 6.780 × 10−2 m, because the story drifts did not exceed 1/75 in the results of pushovers 1V and 2V.
Figure 21b shows the variation of the scale factor λ1 corresponding to each D1U*. The evaluated λ1bi corresponding to D1U*bi = 7.55 × 10−2 m is equal to 0.330. In the second mode, because D2V*bi is equal to D2V*uni, the evaluated λ2bi is equal to CI, uni = 0.361. Therefore, the evaluated capacity index that considers the bidirectional excitation is the minimum of λ1bi = 0.330 and λ2bi = 0.361. Thus, the evaluated CI, bi is equal to 0.330.

4.3. Validation of the Evaluation Procedure by Nonlinear Time-history Analysis

4.3.1. Analysis Cases

The validity of the procedure was evaluated as follows. Nonlinear time-history analyses were carried out using the 10 sets of artificial ground motions presented in Section 3.3. The peak story drifts at columns A3B3, A1B1, and A3B1 were compared with their drift limit (R = 1/75). In the first case, the scale factor of the ground motions (λ) was taken as the capacity index considering the unidirectional ground motion, which was evaluated above (λ = CI, uni = 0.361). In the second case, λ was taken as the capacity index that considers the bidirectional ground motion (λ = CI, bi = 0.330). Note that, in both cases, the incidence angle of the major component with respect to the X-axis (ψ) was taken as 0°, 45°, 90°, and 135°. Therefore, 10 × 4 = 40 analyses were carried out for each case.
In this study, the peak story drifts at columns A3B3, A1B1, and A3B1 were predicted by MABPA [26] with the following modification: the response spectrum used to predict the peak response of the second mode (D2V*uni and A2V*uni) was modified by dividing |sinΔψ12|, as discussed in Section 2.3.4.

4.3.2. Analysis Results

Figure 22 and Figure 23 show the peak drift of columns A1B1, A3B1, and A3B3, as obtained by the time-history analyses for the cases of λ = CI, uni = 0.361 and λ = CI, bi = 0.330, respectively. In these figures, the predicted peak drift by MABPA is also shown.
When λ = CI, uni = 0.361, the mean value of the peak drift at all three columns notably exceeded its drift limit (R = 1/75), as shown in Figure 22. In addition, the predicted peak drift by MABPA exceeds its drift limit (R = 1/75) at column A1B1 and A3B1. While in the case of λ = CI, bi = 0.330, the mean peak drift at column A1B1 and A3B1 was closer to the drift limit than the case of λ = CI, uni = 0.361, as shown in Figure 23. In addition, the predicted peak drift by MABPA is very close to R = 1/75 at column A1B1.
Therefore, the accuracy of the evaluated seismic capacity index for this building was satisfactory. Moreover, index λ = CI, bi was better than index λ = CI, uni because in this building, the effect of bidirectional excitation was significant. In addition, note that the accuracy of the predicted peak drift by MABPA for this building was satisfactory.

5. Comparisons of the Evaluated Seismic Capacity of the Main Uto City Hall Building and the Response Spectrum of Recorded Earthquakes

In this section, the seismic capacity of the main Uto City Hall building is evaluated, considering the variations of modeling of strength degradation behavior and basement spring. Then the evaluated seismic capacity curves are compared to the response spectrum of the first and second earthquakes. For the evaluation of the seismic capacity of each model, the pushover-based method described in Section 2 is applied, considering the code-specific spectrum shown in Section 3.3.

5.1. Analysis Cases

The parameters considered in this section were (i) the force-deformation of the diagonal braces for concrete wall, and (ii) the stiffness of the vertical and rotational springs at each basement, as shown in Figure 24. We assume the shear strain γu to be 0.01, 0.02, and 0.04, as shown in Figure 24a. In each case, the stiffness of springs at each basement was taken as 100%, and 50% of the calculated values described in Section 3.2, as shown in Figure 24b. Table 1 shows the list of 6 structural models for the parametric study. Note that Model RuW1-100 is the model analyzed for the validation of the evaluation procedure in the previous section.

5.2. Analysis Results

5.2.1. Capacity Curves and Seismic Capacity Indices

Table 2 summarizes the evaluation results of each model. As is shown in this table, the range of evaluated CI, bi is from 0.324 to 0.378. This table also shows that the displacement limit is reduced by about 10% to 15% by considering the bidirectional excitation, while the reduction of the capacity index is about 7% to 11%.
Figure 25 shows the capacity curve of each model. In this figure, the points of displacement limit considering unidirectional (the point at D1U*uni) and bidirectional (the point at D1U*bi) excitations are shown. As is shown in this figure, the difference of the capacity curve of Models RuW1-100 and RuW2-100 (and also Models RuW1-050 and RuW2-050) is significant, while the difference of Models RuW2-100 and RuW4-100 (and also Models RuW2-050 and RuW4-050) is limited.
Figure 26 compares the evaluation results of the capacity index that considers the bidirectional ground motion in each model. As is shown in Figure 26a, the difference between Models RuW1-100 and RuW2-100 is noticeable, while the difference between Models RuW2-100 and RuW4-100 is very small. The same trend can be found in the case of which the stiffness of basement spring is reduced by 50% (Figure 26b).

5.2.2. Calculated Peak Story Drift and the Structural Damages

Next, the peak story drift and the structural damages of each model in the case the scaling factor λ = CI, bi is calculated according to MABPA. Note that the response spectrum used in the calculation of the peak response of the second modal response is divided by |sinΔψ12|, as is discussed in Section 2.3.4.
Figure 27 and Figure 28 shows the calculated peak story drift of columns A1B1, A3B1, and A3B3. The peak story drift at columns A1B1 in the second and third story are very close to the limit value (R = 1/75), as expected in all models. This is consistent with the preliminary seismic evaluation results shown in Figure 8: the most critical story according to the preliminary evaluation results is the second story. The difference in each model is noticeable in column A3B3. As shown in Figure 27c, the peak drift of Model RuW1-100 is larger than that of Models RuW2-100 and RuW4-100 in the second story, while in the other stories the peak drift of Model RuW1-100 is smaller than those in other two models. In the case where the stiffness of the basement spring is reduced by 50%, shown in Figure 28c, the peak story drift of Model-RuW1-050 is very close to Models RuW2-050 and RuW4-050 in the second story, while in the other stories the peak drift of Model RuW1-050 is smaller than those in other two models.
Figure 29 shows the distribution of the yielding hinges in each model at frames B1 and A3. As shown in this figure, the damage of frame B1 is more significant than A3, in all models: there are no yielding hinges in frame A3 except Models RuW1-100 and RuW4-050, while there are yielding hinges at the end of the beams in levels Z2 and Z3 and the columns in the second and third stories in frame B1. Especially, in case of Models RuW4-100, RuW2-050, and RuW4-050 (Figure 29c,e,f), the flexural yielding occurs at top of column A2B1 in the third story, and the right-side end of the beam A1-A2 in level Z3 and the left-side end of the beam A2-A3 in level Z3, shown in the dotted rectangle. Note that the yielding strength of the beam ends, where the yielding hinges are made, is reduced because of the yielding of beam-column joint. Therefore, in case of these three models, the yielding of beam-column joint and the column end occur simultaneously at this “+”-shaped beam-column joint.
As shown in Figure 6, the failure of the “+”-shaped beam-column joint at frame B1 was observed above column A2B1 in the third story, which is consistent to the analysis results of Models RuW4-100, RuW2-050, and RuW4-050. Therefore, it can be concluded the analysis results of the three models may explain the part of the observed damage of the main Uto City Hall building: the reason why the damage in frame B1 was more severe than that in frame A3 is that, due to the effect of torsion, the response of frame B1 is more critical. Considering that these three models may represent a part of the real behavior of this building, the evaluated seismic capacity index of this building is approximately 0.36 to 0.38, as shown in Table 2 (shown in the highlighted).

5.2.3. Comparisons with the Evaluated Seismic Capacity and the Response Spectrum of Recorded Earthquakes

Next, the capacity curves of six models are compared to the response spectrum of ground motions recorded during the first (14 April 2016) and second (16 April 2016) earthquakes. The ground motion records used in this analysis were recorded at K-NET Uto station, and details regarding these records can be found in the Appendix C. As described in Section 2, the evaluation of the capacity curve considering bidirectional excitation is based on the assumption that the spectra of the two horizontal components are identical. Of course, this assumption is not true in the real (recorded) ground motion. In this study, the medium of the set of geometrical means obtained using all possible rotations between 0 and 90 degrees (GMRotD50) defined by Boore et al. [44] with the damping ratio equals 0.05 are used as the demand spectrum. The capacity curve of each model is factored by 1/F(nh1eq) for the consideration of equivalent damping, as follows:
1 F ( h 1 e q n ) = 1 + 10 h 1 e q n 1.5 .
Figure 30 and Figure 31 show the comparisons of the factored capacity curve and the demand curve of the first and the second earthquake, respectively.
In the case of the first earthquake shown in Figure 30, the capacity curves of all models have no intersection point to the demand curve. Therefore, the response of all building models exceeds the displacement limit D1U*bi. This implies that the peak drift may exceed 1/75, and the flexural yielding of members shown in Figure 29 may occur during the first earthquake.
Next, in the case of the second earthquake shown in Figure 31, the capacity curves of all models have no intersection point to the demand curve. From the comparisons of two figures, it is seen that the response of all building models during the second earthquake is significantly larger than that during the first earthquake.

5.3. Discussions

Since the damages to this building were first reported on the morning of 16 April 2016, after the second earthquake, the most damage to the building had occurred during the second earthquake. However, the results shown in Figure 30 implies that the main Uto City Hall building may suffer some level of damage during the first earthquake, which occurred on 14 April 2016. Based on the results shown in Figure 30 and Figure 31, the author thinks that damage during the first earthquake, especially the simultaneous yielding of the beam-column joint and the top of column A2B1 in the third story, may be the trigger of the partial collapse in the fourth story occurred during the 16 April earthquake. The simultaneous yielding of the beam-column joint and the top of column A2B1 in the third story may lead to the unstable behavior of the joint; this may cause the whipping behavior, in which the upper part of building oscillates significantly like the head of a whip, which is when, the author thinks, that the severe damage to the fourth and fifth stories occurred.

6. Conclusions

In this paper, the seismic capacity of the main Uto City Hall building, which incurred serious damages during the 2016 Kumamoto Earthquake, was evaluated using the procedure based on the MABPA previously proposed by the author. The main conclusions drawn by this study are as follows.
  • A pushover-based procedure to evaluate the seismic capacity index considering bidirectional excitation is presented and its accuracy is verified by the nonlinear time-history analysis. The results show that the accuracy of the evaluated seismic capacity index for this building was satisfactory. Moreover, index CI, bi was better than index CI, uni because in this building, the effect of bidirectional excitation was significant.
  • The calculated damages of structural models presented in this paper may explain the part of the observed damage of the main Uto City Hall building during the 2016 Kumamoto Earthquake. The reason why the damage in frame B1 was more severe than that in frame A3 can be explained that, due to the effect of torsion, the response of frame B1 is more critical. Moreover, the severe failure of the beam-column joint at frame B1 observed above column A2B1 in the third story is consistent with the analysis results of the three analysis models presented in this paper.
  • The seismic capacity index of the main Uto City Hall building evaluated by the proposed procedure is approximately 0.36 to 0.38. Therefore, this building could withstand only approximately one-third of the design ground motion specified in the current seismic code of Japan.
  • During the first earthquake, which occurred on the 14 April, the main Uto City Hall building may suffer some structural damages, including the yielding of the beam-column joint. According to the comparison of the capacity curve of this building and the demand curve of the first earthquake, the capacity curve has no intersection point to the demand curve.
So far, unfortunately, the ability of the analysis model in this study is only to reveal the possible damage of this building during the first earthquake. Therefore, the author thinks the further improvement will be needed to explain the whole behavior during multiple earthquakes as follows: (i) the modeling of the beam-column joint yielding behavior which can reproduce the unstable behavior of the joint after failure, and (ii) the upgrade of the pushover analysis technique to overcome the limitation of symmetric behavior of all nonlinear springs while keeping numerical stability, applicable for the case when the strength degradation is significant.
It should also be noted that the discussions in Section 5.2.3 rely on the assumption that the GMRotD50 spectrum can represent the bidirectional ground motion in real earthquakes. As the epicenter of the 14 and 16 April earthquakes are not so far from the Uto City Hall (epicentral distance: 15 km and 12 km, respectively [45]), the effect of directivity may not be negligible. Further work will be prepared for the discussion of this issue.

Author Contributions

The all contributions according to this article had been made by the first (single) author.

Funding

This research received no external funding.

Acknowledgments

The author wishes to thank the Uto City Hall officials who provided the original structural designs and other material related to the main buildings. For the analysis in Section 5, the records monitored by National Research Institute for Earth Science and Disaster Prevention were used.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Modification Factor for Predicting the Peak Response of The Equivalent SDOF Model Representing the Second Mode

Figure A1 shows the two equivalent SDOF models representing the first and second modal responses, respectively. Let Equation (A1) be the equation of motion for an N-story irregular building model.
Figure A1. Equivalent SDOF model representing each mode: (a) first mode and (b) second mode.
Figure A1. Equivalent SDOF model representing each mode: (a) first mode and (b) second mode.
Buildings 09 00140 g0a1
M d ¨ ( t ) + C d ˙ ( t ) + f R ( t ) = M { α ξ a g ξ ( t ) + α ς a g ς ( t ) } ,
α ξ = { cos ψ cos ψ sin ψ sin ψ 0 0 } T ,
α ς = { sin ψ sin ψ cos ψ cos ψ 0 0 } T .
According to previous studies [26,29], the equivalent SDOF model representing the first and second modal response can be expressed as Equations (A4) and (A5), respectively.
D ¨ 1 U * ( t ) + C 1 U * M 1 U * D ˙ 1 U * ( t ) + A 1 U * ( t ) = { φ 1 T M α ξ φ 1 T M α U a g ξ ( t ) + φ 1 T M α ς φ 1 T M α U a g ς ( t ) } ,
α ς = { sin ψ sin ψ cos ψ cos ψ 0 0 } T .
According to the equivalent SDOF model for the first mode, Equation (A4) can be easily rewritten as follows. Additionally, by considering the relationship of angles ψ and ψ1 (Equation (A6)), Equations (A7) and (A4) can be rewritten as Equation (A8).
φ 1 T M α ξ φ 1 T M α U = cos ( ψ ψ 1 ) = cos Δ ψ , φ 1 T M α ς φ 1 T M α U = sin Δ ψ ,
a g U ( t ) = a g ξ ( t ) cos Δ ψ + a g ς ( t ) sin Δ ψ ,
D ¨ 1 U * ( t ) + C 1 U * M 1 U * D ˙ 1 U * ( t ) + A 1 U * ( t ) = a g U ( t ) .
Equation (A8) is the motion equation representing the first modal response, as has been discussed in previous reports [26,29]. Therefore, no modifications are made according to the first mode. One issue that should be discussed is that the principal axis of the second modal response is not orthogonal to the U-axis. Let Δψ12 be the angle between the principal axes of the first and second modal responses (Δψ12 = ψ1ψ2). Then, the following relationship can be obtained:
φ 2 T M α ξ φ 2 T M α V = cos ( ψ ψ 2 ) sin ( ψ 1 ψ 2 ) = ( cos Δ ψ 12 sin Δ ψ 12 ) cos Δ ψ sin Δ ψ ,
φ 2 T M α ς φ 2 T M α V = sin ( ψ ψ 2 ) sin ( ψ 1 ψ 2 ) = cos Δ ψ + ( cos Δ ψ 12 sin Δ ψ 12 ) sin Δ ψ .
Therefore, by considering Equations (A11) and (A12), Equation (A5) can be rewritten as Equation (A13):
{ φ 2 T M α ξ φ 2 T M α V a g ξ ( t ) + φ 2 T M α ς φ 2 T M α V a g ς ( t ) } = { a g V ( t ) + cos Δ ψ 12 sin Δ ψ 12 a g U ( t ) } ,
a g V ( t ) = a g ξ ( t ) sin Δ ψ + a g ς ( t ) cos Δ ψ ,
D ¨ 2 V * ( t ) + C 2 V * M 2 V * D ˙ 2 V * ( t ) + A 2 V * ( t ) = { a g V ( t ) + cos Δ ψ 12 sin Δ ψ 12 a g U ( t ) } .
Equation (A13) is identical to the motion equation of the equivalent SDOF model representing the second modal response, which has been proposed by previous studies [26,29], if the principal axis of the second modal response is orthogonal to the U-axis (cosΔψ12 = 0). This equation also implies that the contribution of agU(t) should be considered to estimate the peak response of the second mode, if the principal axis of the second modal response is not orthogonal to the U-axis. In this study, this contribution was considered by factoring the spectral acceleration-displacement curve as follows. Additionally, the factored pseudo acceleration spectrum (pSAV’(T)) is defined by Equation (A14).
p S A V ( t ) = { p S A V ( T ) } 2 + { cos Δ ψ 12 sin Δ ψ 12 S A U p ( T ) } 2 .
Because the spectra of the two horizontal components were assumed to be identical (as expressed by Equation (1)) Equation (A14) can be rewritten as follows:
p S A V ( t ) = S A U p ( T ) 1 + ( cos Δ ψ 12 sin Δ ψ 12 ) 2 = 1 | sin Δ ψ 12 | S A p ( T ) .
Therefore, in this study, the response spectrum used to predict the peak response (D2V*uni and A2V*uni) was modified by dividing |sinΔψ12|, as described in Section 2.3.4 and Section 2.3.6.

Appendix B. Building Damage

Figure A2 shows photographs of the main Uto City Hall building taken after the 16 April earthquake. Figure A3, Figure A4 and Figure A5 show the simplified elevations and photographs of each elevation.
Figure A2. Views of the main Uto City Hall building. Photographs were taken from (a) the south, (b) the west, (c) the north, and (d) the northeast.
Figure A2. Views of the main Uto City Hall building. Photographs were taken from (a) the south, (b) the west, (c) the north, and (d) the northeast.
Buildings 09 00140 g0a2
Figure A3. Simplified elevation and photographs of the south elevation. Photographs show damage to (a) a fifth-story column, (b) a beam–column joint at the top of the fourth story, and (c) a third story column.
Figure A3. Simplified elevation and photographs of the south elevation. Photographs show damage to (a) a fifth-story column, (b) a beam–column joint at the top of the fourth story, and (c) a third story column.
Buildings 09 00140 g0a3
Figure A4. Simplified elevation and photographs of the west elevation. Photographs show (a) damage to the corner column on the fifth story, (b) partial collapse of the fourth story, (c) damage to the corner column on the fourth story.
Figure A4. Simplified elevation and photographs of the west elevation. Photographs show (a) damage to the corner column on the fifth story, (b) partial collapse of the fourth story, (c) damage to the corner column on the fourth story.
Buildings 09 00140 g0a4
Figure A5. Simplified elevation and photographs of the north elevation. Photographs show damage to (a) the top of corner columns and beams on the roof floor, (b) a beam on the fifth floor, and (c) the bottom of a fourth-story column and a fourth-floor beam.
Figure A5. Simplified elevation and photographs of the north elevation. Photographs show damage to (a) the top of corner columns and beams on the roof floor, (b) a beam on the fifth floor, and (c) the bottom of a fourth-story column and a fourth-floor beam.
Buildings 09 00140 g0a5

Appendix C. Record Observed in the 14th and 16th Earthquake at K-NET Uto Station

Figure A6 and Figure A7 show the acceleration of two horizontal components and orbit of the record observed on the 14 and 16 April, respectively, of the April earthquake at the K-NET Uto Station. The major and minor axes in the orbit are determined according to Penzien and Watabe [46]. Figure A8 shows the pseudo-acceleration response spectrum of the two recorded ground motions observed at the K-NET Uto Station.
Figure A6. Acceleration of EW- and NS-components and orbit of the record observed in the 14 April earthquake at the K-NET Uto Station.
Figure A6. Acceleration of EW- and NS-components and orbit of the record observed in the 14 April earthquake at the K-NET Uto Station.
Buildings 09 00140 g0a6
Figure A7. Acceleration of EW- and NS-components and orbit of record observed in the 16 April earthquake at the K-NET Uto Station.
Figure A7. Acceleration of EW- and NS-components and orbit of record observed in the 16 April earthquake at the K-NET Uto Station.
Buildings 09 00140 g0a7
Figure A8. Pseudo acceleration response spectrum of two recorded ground motion at the K-NET Uto Station: (a) 14 April earthquake, and (b) 16 April earthquake.
Figure A8. Pseudo acceleration response spectrum of two recorded ground motion at the K-NET Uto Station: (a) 14 April earthquake, and (b) 16 April earthquake.
Buildings 09 00140 g0a8

References

  1. Fujii, K.; Yoshida, S.; Nishimura, T.; Furuta, T. Observations of damage to Uto City Hall suffered in the 2016 Kumamoto Earthquake. In Proceedings of the 8th European Workshop on the Seismic Behaviour of Irregular and Complex Structures, Bucharest, Romania, 19–20 October 2017. [Google Scholar]
  2. Fujii, K. Preliminary evaluation of seismic capacity and torsional irregularity of Uto City Hall damaged in 2016 Kumamoto Earthquake. In Proceedings of the 8th European Workshop on the Seismic Behaviour of Irregular and Complex Structures, Bucharest, Romania, 19–20 October 2017. [Google Scholar]
  3. Saiidi, M.; Sozen, M.A. Simple nonlinear seismic analysis of R/C structures. J. Struct. Div. 1981, 107, 937–952. [Google Scholar]
  4. Fajfar, P.; Fischinger, M. N2-A method for non-linear seismic analysis of regular buildings. In Proceedings of the 9th World Conference on Earthquake Engineering, Tokyo-Kyoto, Japan, 2–9 August 1988. [Google Scholar]
  5. Fajfar, P. A nonlinear analysis method for performance-based seismic design. Earthq. Spectra 2000, 16, 573–592. [Google Scholar] [CrossRef]
  6. ATC-40. Seismic Evaluation and Retrofit of Concrete Buildings; Applied Technology Council: Redwood City, CA, USA, 1996; Volume 1. [Google Scholar]
  7. FEMA. NEHRP Guidelines for the Seismic Rehabilitation of Buildings; FEMA Publication 273; Federal Emergency Management Agency: Washington, DC, USA, 1997.
  8. CEN. Eurocode 8—Design of Structures for Earthquake Resistance. Part 1: General Rules, Seismic Actions and Rules for Buildings; European standard EN 1998-1; European Committee for Standardization: Brussels, Belgium, 2004. [Google Scholar]
  9. ASCE. Seismic Evaluation and Retrofit of Existing Buildings; ASCE standard, ASCE/SEI41-17; American Society of Civil Engineering: Reston, VA, USA, 2017. [Google Scholar]
  10. Peruŝ, I.; Fajfar, P. On the inelastic torsional response of single-storey structures under bi-axial excitation. Earthq. Eng. Struct. Dyn. 2005, 34, 931–941. [Google Scholar] [CrossRef]
  11. Fajfar, P.; Maruŝić, D.; Peruŝ, I. Torsional effects in the pushover-based seismic analysis of buildings. J. Earthq. Eng. 2005, 9, 831–854. [Google Scholar] [CrossRef]
  12. Kreslin, M.; Fajfar, P. Seismic evaluation of an existing complex RC building. Bull. Earthq. Eng. 2010, 8, 363–385. [Google Scholar] [CrossRef]
  13. Kreslin, M.; Fajfar, P. The extended N2 method considering higher mode effects in both plan and elevation. Bull. Earthq. Eng. 2012, 10, 695–715. [Google Scholar] [CrossRef]
  14. Dolšek, M.; Fajfar, P. Simplified probabilistic seismic performance assessment of plan-asymmetric buildings. Earthq. Eng. Struct. Dyn. 2007, 36, 2021–2041. [Google Scholar] [CrossRef]
  15. Chopra, A.K.; Goel, R.K. A modal pushover analysis procedure to estimate seismic demands for unsymmetric-plan buildings. Earthq. Eng. Struct. Dyn. 2004, 33, 903–927. [Google Scholar] [CrossRef] [Green Version]
  16. Reyes, J.C.; Chopra, A.K. Three-dimensional modal pushover analysis of buildings subjected to two components of ground motion, including its evaluation for tall buildings. Earthq. Eng. Struct. Dyn. 2011, 40, 789–806. [Google Scholar] [CrossRef]
  17. Reyes, J.C.; Chopra, A.K. Evaluation of three-dimensional modal pushover analysis for unsymmetric-plan buildings subjected to two components of ground motion. Earthq. Eng. Struct. Dyn. 2011, 40, 1475–1494. [Google Scholar] [CrossRef]
  18. Manoukas, G.; Athanatopoulou, A.; Avramidis, I. Multimode pushover analysis for asymmetric buildings under biaxial seismic excitation based on a new concept of the equivalent single degree of freedom system. Soil Dyn. Earthq. Eng. 2012, 38, 88–96. [Google Scholar] [CrossRef]
  19. Manoukas, G.; Avramidis, I. Evaluation of a multimode pushover procedure for asymmetric in plan buildings under biaxial seismic excitation. Bull. Earthq. Eng. 2014, 12, 2607–2632. [Google Scholar] [CrossRef]
  20. Belejo, A.; Bento, R. Improved modal pushover analysis in seismic assessment of asymmetric plan buildings under the influence of one and two horizontal components of ground motions. Soil Dyn. Earthq. Eng. 2016, 87, 1–15. [Google Scholar] [CrossRef]
  21. Soleimani, S.; Aziminejad, A.; Moghadam, A.S. Extending the concept of energy-based pushover analysis to assess seismic demands of asymmetric-plan buildings. Soil Dyn. Earthq. Eng. 2017, 93, 29–41. [Google Scholar] [CrossRef]
  22. Reyes, J.C.; Riaño, A.C.; Kalkan, E.; Arango, C.M. Extending modal pushover-based scaling procedure for nonlinear response history analysis of multi-story unsymmetric-plan buildings. Eng. Struct. 2015, 88, 125–137. [Google Scholar] [CrossRef]
  23. Bosco, M.; Ghersi, A.; Marino, E.M. Corrective eccentricities for assessment by nonlinear static method of 3D structures subjected to bidirectional ground motion. Earthq. Eng. Struct. Dyn. 2012, 41, 1751–1773. [Google Scholar] [CrossRef]
  24. Bosco, M.; Ferrara, G.A.F.; Ghersi, A.; Marino, E.M.; Rossi, P.P. Predicting displacement demand of multi-storey asymmetric buildings by nonlinear static analysis and corrective eccentricities. Eng. Struct. 2015, 99, 373–387. [Google Scholar] [CrossRef]
  25. Fujii, K. Nonlinear static procedure for multi-story asymmetric frame buildings considering bi-directional excitation. J. Earthq. Eng. 2011, 15, 245–273. [Google Scholar] [CrossRef]
  26. 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]
  27. Fujii, K. Application of the pushover-based procedure to predict the largest peak response of asymmetric buildings with buckling-restrained braces. In Proceedings of the 5th International Conference on Computational Methods in Structural Dynamics and Earthquake Engineering (COMPDYN 2015), Crete Island, Greece, 25–27 May 2015. [Google Scholar]
  28. Fujii, K. Assessment of pushover-based method to a building with bidirectional setback. Earthq. Struct. 2016, 11, 421–443. [Google Scholar] [CrossRef]
  29. Fujii, K. Prediction of the peak seismic response of asymmetric buildings under bidirectional horizontal ground motion using equivalent SDOF model. Jpn. Archit. Rev. 2018, 1, 29–43. [Google Scholar] [CrossRef] [Green Version]
  30. López, A.; Torres, R. The critical angle of seismic incidence and the maximum structural response. Earthq. Eng. Struct. Dyn. 1997, 26, 881–894. [Google Scholar] [CrossRef]
  31. López, A.; Hernández, J.J.; Bonilla, R.; Fernández, A. Response spectra for multicomponent structural analysis. Earthq. Spectra 2006, 22, 85–113. [Google Scholar] [CrossRef]
  32. NZSEE/MBIE. The Seismic Assessment of Existing Buildings. Technical Guidelines for Engineering Assessments. Part A: Assessment Objectives and Principles. New Zealand Society for Earthquake Engineering: Wellington, 2017. Available online: www.EQ-Assess.org.nz (accessed on 26 February 2019).
  33. NZSEE/MBIE. The Seismic Assessment of Existing Buildings. Technical Guidelines for Engineering Assessments. Part B: Initial Seismic Assessment. New Zealand Society for Earthquake Engineering: Wellington, 2017. Available online: www.EQ-Assess.org.nz (accessed on 26 February 2019).
  34. NZSEE/MBIE. The Seismic Assessment of Existing Buildings. Technical Guidelines for Engineering Assessments. Part C: Detailed Seismic Assessment. New Zealand Society for Earthquake Engineering: Wellington, 2017. Available online: www.EQ-Assess.org.nz (accessed on 26 February 2019).
  35. Vecchio, C.D.; Gentile, R.; Ludovico, M.D.; Uva, G.; Pampanin, S. Implementation and Validation of the Simple Lateral Mechanism Analysis (SLaMA) for the Seismic Performance Assessment of a Damaged Case Study Building. J. Earthq. Eng. 2018, 1–32. Available online: www.tandfonline.com (accessed on 26 February 2019). [CrossRef]
  36. Otani, S. New seismic design provision in Japan. In Proceedings of the Second U.S.-Japan Workshop on Performance-Based Earthquake Engineering Methodology for Reinforced Concrete Structures, PEER Report. Sapporo, Hokkaido, Japan, 3–14 October 2000. [Google Scholar]
  37. Shiga, T. Earthquake damage and the amount of walls in reinforced concrete buildings. In Proceedings of the 6th World Conference on Earthquake Engineering, New Delhi, India, 9–14 January 1977. [Google Scholar]
  38. 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]
  39. 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; Volume 1, pp. 841–844. [Google Scholar]
  40. 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; Volume 2, pp. 1960–1969. [Google Scholar]
  41. Otani, S. Hysteresis models of reinforced concrete for earthquake response analysis. J. Fac. Eng. 1981, 36, 125–156. [Google Scholar]
  42. BRI; BCJ. Guideline for the Generation of Earthquake Motions for Seismic Design of Buildings (Draft); The Building Center of Japan, 1993. (In Japanese). Available online: www.bcj.or.jp/download/wave/ (accessed on 27 July 2010).
  43. Otani, S. Japanese seismic design of high-rise reinforced concrete buildings—An example of performance-based design code and state of practices. In Proceedings of the 13th World Conference on Earthquake Engineering, Vancouver, Canada, 1–6 August 2004. Paper No. 5010. [Google Scholar]
  44. Boore, D.M.; Watson-Lamprey, J.; Abrahamson, N.A. Orientation-Independent Measures of Ground Motion. Bull. Seismol. Soc. Am. 2007, 96, 1502–1511. [Google Scholar] [CrossRef]
  45. Strong-Motion Seismograph Network (K-NET, KIK-NET) Home Page. Available online: www.kyoshin.bosai.go.jp/kyoshin/ (accessed on 4 May 2016).
  46. Penzien, J.; Watabe, M. Characteristics of 3-Dimensional Earthquake Ground Motions. Earthq. Eng. Struct. Dyn. 1975, 3, 365–373. [Google Scholar] [CrossRef]
Figure 1. Definition of capacity index CI.
Figure 1. Definition of capacity index CI.
Buildings 09 00140 g001
Figure 2. Comparisons of parameters between an N-story irregular building and the equivalent SDOF model: (a) an N-story irregular building deformed similar to the first mode shape and (b) equivalent SDOF model representing the first modal response.
Figure 2. Comparisons of parameters between an N-story irregular building and the equivalent SDOF model: (a) an N-story irregular building deformed similar to the first mode shape and (b) equivalent SDOF model representing the first modal response.
Buildings 09 00140 g002
Figure 3. Structural plan of the main Uto City Hall building: (a) level Z0, (b) level Z1 to Z4, and (c) level Z5.
Figure 3. Structural plan of the main Uto City Hall building: (a) level Z0, (b) level Z1 to Z4, and (c) level Z5.
Buildings 09 00140 g003
Figure 4. Simplified structural elevation of Uto City Hall: (a) frame Y4, (b) frame B1, and (c) frame B2.
Figure 4. Simplified structural elevation of Uto City Hall: (a) frame Y4, (b) frame B1, and (c) frame B2.
Buildings 09 00140 g004
Figure 5. View of Uto City Hall after the 2016 Kumamoto Earthquake: (a) south (frame A3) and (b) southwest (frames A3 and B1).
Figure 5. View of Uto City Hall after the 2016 Kumamoto Earthquake: (a) south (frame A3) and (b) southwest (frames A3 and B1).
Buildings 09 00140 g005
Figure 6. Failure of the beam-column joint above the column A2B1 in the third story.
Figure 6. Failure of the beam-column joint above the column A2B1 in the third story.
Buildings 09 00140 g006
Figure 7. Sections of Column A1B1, A2B1, and A2B2 [2].
Figure 7. Sections of Column A1B1, A2B1, and A2B2 [2].
Buildings 09 00140 g007
Figure 8. Preliminary seismic capacity evaluation results of the main building of the Uto City Hall [2]: (a) the case in which the building is assumed to behave as a unite building, and (b) the case in which the stair and office blocks are assumed to behave independently.
Figure 8. Preliminary seismic capacity evaluation results of the main building of the Uto City Hall [2]: (a) the case in which the building is assumed to behave as a unite building, and (b) the case in which the stair and office blocks are assumed to behave independently.
Buildings 09 00140 g008
Figure 9. Model of each frame: (a) frame B1 and (b) frame Y5.
Figure 9. Model of each frame: (a) frame B1 and (b) frame Y5.
Buildings 09 00140 g009
Figure 10. Force-deformation relationship of nonlinear flexural spring: (a) envelope and (b) hysteresis rule.
Figure 10. Force-deformation relationship of nonlinear flexural spring: (a) envelope and (b) hysteresis rule.
Buildings 09 00140 g010
Figure 11. Simplified model of the concrete wall: (a) model, (b) envelope of the force-deformation relationship of the nonlinear brace (concrete wall), and (c) hysteresis rule.
Figure 11. Simplified model of the concrete wall: (a) model, (b) envelope of the force-deformation relationship of the nonlinear brace (concrete wall), and (c) hysteresis rule.
Buildings 09 00140 g011
Figure 12. The shape of the first three natural modes of the building model in the elastic range.
Figure 12. The shape of the first three natural modes of the building model in the elastic range.
Buildings 09 00140 g012
Figure 13. Response spectra of the elastic acceleration.
Figure 13. Response spectra of the elastic acceleration.
Buildings 09 00140 g013
Figure 14. Example of accelerations and orbit of the artificial ground motion.
Figure 14. Example of accelerations and orbit of the artificial ground motion.
Buildings 09 00140 g014
Figure 15. Nonlinear properties of the equivalent SDOF model representing the first mode response obtained from the pushover analysis: (a) A1U*-D1U* relationship, (b) T1eq-D1U* relationship, and (c) h1eq-D1U* relationship.
Figure 15. Nonlinear properties of the equivalent SDOF model representing the first mode response obtained from the pushover analysis: (a) A1U*-D1U* relationship, (b) T1eq-D1U* relationship, and (c) h1eq-D1U* relationship.
Buildings 09 00140 g015
Figure 16. Variations in the first mode parameters based on the pushover analysis: (a) effective modal mass ratio (m1U*), (b) incidence angle of the principal axis of the first modal response (ψ1), and (c) torsional index of the first mode (Rρ1).
Figure 16. Variations in the first mode parameters based on the pushover analysis: (a) effective modal mass ratio (m1U*), (b) incidence angle of the principal axis of the first modal response (ψ1), and (c) torsional index of the first mode (Rρ1).
Buildings 09 00140 g016
Figure 17. Roof displacement responses based on the pushover analysis: (a) frames B1 to B3, (b) frames A1 to A3, and (c) columns A1B1, A3B1, and A3B3.
Figure 17. Roof displacement responses based on the pushover analysis: (a) frames B1 to B3, (b) frames A1 to A3, and (c) columns A1B1, A3B1, and A3B3.
Buildings 09 00140 g017
Figure 18. The shape of the first natural mode at D1U* = 8.420 × 10−2 m.
Figure 18. The shape of the first natural mode at D1U* = 8.420 × 10−2 m.
Buildings 09 00140 g018
Figure 19. Evaluation of the capacity index considering the unidirectional ground motion CI, uni: (a) the evaluation of the peak response for each scaling factor λ, (b) variation of the evaluated scaling factor λ1 corresponding to the equivalent displacement D1U*.
Figure 19. Evaluation of the capacity index considering the unidirectional ground motion CI, uni: (a) the evaluation of the peak response for each scaling factor λ, (b) variation of the evaluated scaling factor λ1 corresponding to the equivalent displacement D1U*.
Buildings 09 00140 g019
Figure 20. Prediction of the peak response for the first and second modal responses (λ = CI, uni = 0.361).
Figure 20. Prediction of the peak response for the first and second modal responses (λ = CI, uni = 0.361).
Buildings 09 00140 g020
Figure 21. Evaluation of the capacity index considering the bidirectional ground motion CI, bi: (a) roof displacement responses at column A1B1 and (b) variation of the scale factor λ.
Figure 21. Evaluation of the capacity index considering the bidirectional ground motion CI, bi: (a) roof displacement responses at column A1B1 and (b) variation of the scale factor λ.
Buildings 09 00140 g021
Figure 22. Comparison of the peak drift for columns (λ = CI, uni = 0.361): (a) column A1B1, (b) column A3B1, and (c) column A3B3.
Figure 22. Comparison of the peak drift for columns (λ = CI, uni = 0.361): (a) column A1B1, (b) column A3B1, and (c) column A3B3.
Buildings 09 00140 g022
Figure 23. Comparison of the peak drift for columns (λ = CI, bi = 0.330): (a) column A1B1, (b) column A3B1, and (c) column A3B3.
Figure 23. Comparison of the peak drift for columns (λ = CI, bi = 0.330): (a) column A1B1, (b) column A3B1, and (c) column A3B3.
Buildings 09 00140 g023
Figure 24. Parameters considered in modeling: (a) the force-deformation of the diagonal braces for concrete wall, and (b) the stiffness of the vertical and rotational springs at each basement.
Figure 24. Parameters considered in modeling: (a) the force-deformation of the diagonal braces for concrete wall, and (b) the stiffness of the vertical and rotational springs at each basement.
Buildings 09 00140 g024
Figure 25. Comparison of the capacity curve of each model: (a) Models RuW1-100, RuW2-100, and RuW4-100, (b) Models RuW1-050, RuW2-050, and RuW4-050.
Figure 25. Comparison of the capacity curve of each model: (a) Models RuW1-100, RuW2-100, and RuW4-100, (b) Models RuW1-050, RuW2-050, and RuW4-050.
Buildings 09 00140 g025
Figure 26. Comparison of the evaluated results of the capacity index considering the bidirectional ground motion: (a) Models RuW1-100, RuW2-100, and RuW4-100, (b) Models RuW1-050, RuW2-050, and RuW4-050.
Figure 26. Comparison of the evaluated results of the capacity index considering the bidirectional ground motion: (a) Models RuW1-100, RuW2-100, and RuW4-100, (b) Models RuW1-050, RuW2-050, and RuW4-050.
Buildings 09 00140 g026
Figure 27. Comparison of the peak drift for columns (λ = CI, bi) for Models RuW1-100, RuW2-100, and RuW4-100: (a) column A1B1, (b) column A3B1, and (c) column A3B3.
Figure 27. Comparison of the peak drift for columns (λ = CI, bi) for Models RuW1-100, RuW2-100, and RuW4-100: (a) column A1B1, (b) column A3B1, and (c) column A3B3.
Buildings 09 00140 g027
Figure 28. Comparison of the peak drift for columns (λ = CI, bi) for Models RuW1-050, RuW2-050, and RuW4-050: (a) column A1B1, (b) column A3B1, and (c) column A3B3.
Figure 28. Comparison of the peak drift for columns (λ = CI, bi) for Models RuW1-050, RuW2-050, and RuW4-050: (a) column A1B1, (b) column A3B1, and (c) column A3B3.
Buildings 09 00140 g028
Figure 29. Distribution of the yielding hinges at frame B1 and A3 in each model (λ = CI, bi): (a) Model-RuW1-100, (b) Model-RuW2-100, (c) Model-RuW4-100, (d) Model-RuW1-050, (e) Model-RuW2-050, and (f) Model-RuW4-050.
Figure 29. Distribution of the yielding hinges at frame B1 and A3 in each model (λ = CI, bi): (a) Model-RuW1-100, (b) Model-RuW2-100, (c) Model-RuW4-100, (d) Model-RuW1-050, (e) Model-RuW2-050, and (f) Model-RuW4-050.
Buildings 09 00140 g029
Figure 30. Comparison of the factored capacity curve and the demand curve of the first earthquake (14 April 2016): (a) Models RuW1-100, RuW2-100, and RuW4-100, and (b) Models RuW1-050, RuW2-050, and RuW4-050.
Figure 30. Comparison of the factored capacity curve and the demand curve of the first earthquake (14 April 2016): (a) Models RuW1-100, RuW2-100, and RuW4-100, and (b) Models RuW1-050, RuW2-050, and RuW4-050.
Buildings 09 00140 g030
Figure 31. Comparison of the factored capacity curve and the demand curve of the second earthquake (16 April 2016): (a) Models RuW1-100, RuW2-100, and RuW4-100, and (b) Models RuW1-050, RuW2-050, and RuW4-050.
Figure 31. Comparison of the factored capacity curve and the demand curve of the second earthquake (16 April 2016): (a) Models RuW1-100, RuW2-100, and RuW4-100, and (b) Models RuW1-050, RuW2-050, and RuW4-050.
Buildings 09 00140 g031
Table 1. List of structural models considered in the parametric study.
Table 1. List of structural models considered in the parametric study.
Model IDShear Strain γuStiffness Ratio of Spring at Each Basement
Model-RuW1-1000.01100%
Model-RuW2-1000.02
Model-RuW4-1000.04
Model-RuW1-0500.0150%
Model-RuW2-0500.02
Model-RuW4-0500.04
Table 2. Summary of the evaluation results of each model.
Table 2. Summary of the evaluation results of each model.
Model IDDisplacement LimitCapacity Index
D1U*uni
(×10−2 m)
D1U*bi
(×10−2 m)
D1U*bi/D1U*uniCI, uniCI, biCI, bi/CI, uni
RuW1-1008.4207.5470.8960.3610.3300.914
RuW2-1008.7607.8800.9000.3930.3650.929
RuW4-1009.0258.0450.8910.4120.3780.917
RuW1-0508.6167.4800.8680.3560.3240.910
RuW2-0509.1757.9000.8610.3970.3610.909
RuW4-0509.5608.1000.8470.4190.3730.890

Share and Cite

MDPI and ACS Style

Fujii, K. Pushover-Based Seismic Capacity Evaluation of Uto City Hall Damaged by the 2016 Kumamoto Earthquake. Buildings 2019, 9, 140. https://doi.org/10.3390/buildings9060140

AMA Style

Fujii K. Pushover-Based Seismic Capacity Evaluation of Uto City Hall Damaged by the 2016 Kumamoto Earthquake. Buildings. 2019; 9(6):140. https://doi.org/10.3390/buildings9060140

Chicago/Turabian Style

Fujii, Kenji. 2019. "Pushover-Based Seismic Capacity Evaluation of Uto City Hall Damaged by the 2016 Kumamoto Earthquake" Buildings 9, no. 6: 140. https://doi.org/10.3390/buildings9060140

APA Style

Fujii, K. (2019). Pushover-Based Seismic Capacity Evaluation of Uto City Hall Damaged by the 2016 Kumamoto Earthquake. Buildings, 9(6), 140. https://doi.org/10.3390/buildings9060140

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