*Article* **Proposal of Damage Index Ratio for Low- to Mid-Rise Reinforced Concrete Moment-Resisting Frame with Setback Subjected to Uniaxial Seismic Loading**

**Taufiq Ilham Maulana, Badamkhand Enkhtengis and Taiki Saito \***

Department of Architecture and Civil Engineering, Toyohashi University of Technology, Aichi 441-8580, Japan; taufiq.ilham.maulana.mt@tut.jp (T.I.M.); badamkhand.enkhtengis.do@tut.jp (B.E.) **\*** Correspondence: saito.taiki.bv@tut.jp

**Abstract:** A vertical irregularity setback in reinforced concrete (RC) building affects its performance and response especially subjected to earthquake ground motions. It is necessary to understand how the seismic damage is established due to setbacks and avoid damage concentration on the irregularity section. The objective of this study is to propose a formula to estimate the damage distribution along the height of the setback building from a geometric measure of the degree of irregularity. First, previous experimental tests for two types of setback buildings, a towered and a stepped setback frames, were analyzed to verify the accuracy of the frame analysis. The results of the frame analysis considerably matched the experimental test results. Furthermore, to study the relationship between the degree of setback and the distribution of damage, a parametric study was conducted using 35 reinforced concrete setback frames, consisting of models with stepped setback type and towered setback type with different degrees of setback. The inelastic dynamic analyses of all the frames under three earthquake ground motions were conducted. The irregularity indices proposed in literature were adopted to express the degree of setback and the structural damage was expressed by the Park–Ang damage index. Using nonlinear regression analysis, formulas to estimate damage index ratio between two main structure parts (tower and base) from setback indices were proposed. Finally, the proposed formula was applied to the experimental test results to confirm its validity.

**Keywords:** setbacks; earthquake; damage index; seismic evaluation performance; damage distribution

#### **1. Introduction**

#### *1.1. Background*

In recent modern multistory building construction, the implementation of horizontal and vertical irregularities is getting popular due to different building functions and architectural needs. Vertical irregularity often regarded as a setback affects discontinuities in mass, stiffness, and structural strength distributions [1,2]. Under seismic excitations, this irregularity has a considerable effect on structure responses [3]. The existence of the setback may cause torsional response [4] and stress and damage concentration on the irregular part [5].

Studies on the seismic response of buildings with setbacks had been conducted by several researchers. Shakib and Pirizadeh [6] assessed the effect of simultaneous action of two orthogonal ground motion components on the performance of structure with different setback ratios and found that the effect of ground motion components is significant for the seismic performance of setback structures. Lin et al. [7] investigated the modal response history analysis of setback buildings and instead of using the conventional SDOF modal system, the 2DOF modal system was employed in modal response history analysis of building with specific vertical irregularities, to characterize the distinctly different modal response of tower and base structures. It was found that the 2DOF approach adequately captured the characteristic of drastic change in the peak interstory drift ratios between

**Citation:** Maulana, T.I.; Enkhtengis, B.; Saito, T. Proposal of Damage Index Ratio for Low- to Mid-Rise Reinforced Concrete Moment-Resisting Frame with Setback Subjected to Uniaxial Seismic Loading. *Appl. Sci.* **2021**, *11*, 6754. https://doi.org/10.3390/app11156754

Academic Editor: Maria Favvata

Received: 9 June 2021 Accepted: 20 July 2021 Published: 22 July 2021

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

**Copyright:** © 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https:// creativecommons.org/licenses/by/ 4.0/).

tower and base. Georgoussis and Mamou [8] proposed the equations to minimize the rotational response due to the mass eccentricity of setback buildings by changing the location of lateral load resisting elements. The parameters in the proposed equations involved the first mode frequency and the displacement ratio between the tower and the base structure. The proposed equations were also applied for minimizing the rotational response of medium-rise buildings with irregular plans under seismic excitations [9].

Habibi and Asadi [10] developed equations to determine the Park–Ang damage index [11,12] of setback buildings using three independent parameters, namely overall drift (which is the ratio between the maximum roof displacement to the total height of the frame), the natural period, and the irregularity indices of the building. The adopted irregularity indices were originally introduced by Karavasilis et al. [13], involving the geometry of setback buildings, such as the number of stories and bays and the height and the width of setback buildings. Varadharajan et al. [14,15] also predicted the Park–Ang damage index of setback structures using three primary parameters: the ratio of modal participation factor between irregular and regular buildings, the beam to column stiffness ratio, and the displacement ductility. Hait et al. [16] assessed the Park–Ang damage index of low to mid-rise buildings using parameters generated from dynamic analyses, such as the maximum interstory drift, the peak roof displacement, and the maximum joint rotation of the members. Hait et al. [17]. also predicted the Park–Ang damage index by multi-variable regression using the artificial neural network. In these studies, the Park–Ang damage index was computed through nonlinear dynamic analysis for both tower and base sections, and the damage distribution could be retrieved after. This procedure could take time and resources for the dynamic analysis.

The objective of this paper is to propose the formula to determine the damage distribution of setback buildings based on the geometric measure of the degree of irregularity without conducting dynamic analysis. First, to verify the accuracy of the analytical method, the experimental tests for two types of setback buildings are analyzed. The experimental tests are a towered structure of fifteenth-scale 9-story 3-bay RC frame with a setback at level 3 by Wood [18] and a quarter model of 6-story 2-bay setback RC frame by Shahrooz and Moehle [19]. Then, the relationship between damage distribution of the Park–Ang damage index and the irregularity indices are examined by the parametric study of 35 models with different setback types. Finally, the nonlinear regression equations are proposed to determine the building damage distribution from the geometrical irregularity indices. Then, these equations are validated using the previous experimental results.

#### *1.2. Park and Ang Damage Index*

#### 1.2.1. Definition

The damage index (DI) introduced by Park and Ang [11,12] is expressed by Equation (1). This DI value is based on the structural deformation and the hysteretic energy response due to seismic excitations.

$$DI = \frac{u\_m}{u\_u} + \beta \frac{E\_h}{F\_y u\_u} \tag{1}$$

where:

*um*: Maximum displacement response of structure element due to earthquake,

*uu*: Ultimate displacement capacity under a monotonic loading,

*Eh*: Hysteretic energy dissipated by the structural element,

*Fy*: Yield force,

β: Non-negative parameter based on repeated loading effect.

The relationship between the DI and the physical damage appearance for RC elements suggested by Park and Ang [12] is shown in Table 1. The parameters in Equation (1) are explained in the following passages.


**Table 1.** Damage Index for RC elements.

#### 1.2.2. Determination of uu

The ultimate displacement *uu* is formulated as follows.

$$
u\_{
u} = 
\mu\_{
u} 
u\_{
y} \tag{2}$$

where:

*uu*: Ultimate displacement under monotonic loading,

*μu*: Ultimate ductility factor,

*uy*: Yield displacement.

The yield displacement *uy* is defined for each member based on element geometry and material properties. Due to the limitation of the available references stating clearly about the ultimate ductility factor *μu*, in this study, for each member, the *μ<sup>u</sup>* is taken value as 15 by considering the probable failure mode occurs at the beam and column elements to be dominated as flexural failure [11].

#### 1.2.3. Determination of β

The parameter *β* by Park and Ang [12] is presented in Equations (3) and (4).

$$\beta = (-0.447 + 0.073\frac{l}{d} + 0.24n\_0 + 0.314p\_t) \times 0.7^{p\_w} \tag{3}$$

$$m\_0 = \frac{N}{f\_c b d} \tag{4}$$

where:

*l/d*: Shear span ratio (= 1.7 if *l/d* less than 1.7), *n0*: Normalized axial stress (= 0.2 if *n0* less than 0.2), *pt*: Longitudinal bars ratio, in percentage (= 0.75% if *pt* less than 0.75%), *pw*: Transversal bars ratio, *N*: Axial load, *fc*: Concrete compressive strength, *b*: Width of cross section, *d*: Effective height of cross section.

#### 1.2.4. Damage Index of Building Component

The DI value is also suggested for a part of a structure, such as the individual story and the entire structure, by considering a weighting factor that is based on the dissipated hysteretic energy (*Eh*) of each member or the component as shown below, that are proposed originally by Park and Ang [12] and developed by Belkacem et al. [20]. Figure 1 shows the difference between *DIstory* and *DIcomponent*. In the component weighting factor calculation, the dissipated energy (*Eh*) considered is only in the same story, while the calculation of story weighting factor, all energies in all stories are considered.

$$DI\_{Storage} = \sum\_{i=1}^{n\_m} \left(\lambda\_{i,component} DI\_{i,component}\right) \tag{5}$$

$$
\lambda\_{i, component} = \begin{bmatrix} E\_{h,i} \\ \frac{\sum\_{i=1}^{l\_{\text{th}}} E\_{h,i}}{\sum\_{i=1}^{l\_{\text{th}}} E\_{h,i}} \end{bmatrix}\_{component} \tag{6}
$$

$$DI\_{overall} = \sum\_{j=1}^{n\_k} \left(\lambda\_{j,story} DI\_{j,story}\right) \tag{7}$$

$$
\lambda\_{j,story} = \left[\frac{E\_{h,j}}{\sum\_{j=1}^{\mathcal{I}\_k} E\_{h,j}}\right]\_{story} \tag{8}
$$

where:

*i*: number of element member, starting at 1, *j*: number of story, starting at 1, *nk*: total number of building's story, *nm*: total element member (beams and columns) in the reviewed story, *DIi,component*: Damage index for each component such as beam and column, *DIStory*: Damage index for each story, *DIoverall*: Overall damage index, *λi,component*: Weighting factor of each component, *λj,story*: Weighting factor of each story.

**Figure 1.** Illustration about damage index component and damage index story in a frame structure.

#### *1.3. Setback Indices*

Two kinds of vertical irregularities are selected to be studied: stepped and towered. The stepped type is a setback with the upper structure part at one side of the building, while the towered type is the condition when the upper structure part is in the middle region of the building. The illustration of the stepped type is shown in Figure 2a, and the towered type is shown in Figure 2b.

To consider the setback effect, two irregularity indices are adopted. These indices were initially introduced by Mazzolani and Piluso [21], which were then developed by Karavasilis et al. [13]. The indices describe the irregularity due to the presence of setbacks through simple geometrical indices *ϕ<sup>b</sup>* and *ϕ<sup>s</sup>* given by the following Equations (9) and (10).

$$\varphi\_s = \frac{1}{n\_s - 1} \times \sum\_{i=1}^{i=n\_s - 1} \frac{L\_i}{L\_{i+1}} \tag{9}$$

$$\varphi\_b = \frac{1}{n\_b - 1} \times \sum\_{i=1}^{i=n\_b - 1} \frac{H\_i}{H\_{i+1}} \tag{10}$$

where:


**Figure 2.** Frame geometry for irregularity indices definition for: (**a**) stepped type of setbacks; (**b**) towered type of setbacks.

In the case of stepped setback type, since the height *Hi* is counted from the left side, the value of *ϕ<sup>b</sup>* will be different if the tower part is on the right side or on the left side. Therefore, in this study, *Hi* is always counted from the side that has the highest level.

#### **2. Simulation of Setback RC Buildings Tested by Shaking Table**

To verify the accuracy of the frame analysis, shaking table experimental tests of scaled setback RC buildings conducted by Wood in 1985 [18] and Shahrooz and Moehle in 1987 [19] were chosen to be numerically analyzed, and the raw data of both experimental tests were retrieved from DataCenterHub repository [22,23]. A software developed by one of the co-authors, STERA\_3D, is used to perform the nonlinear dynamic structural analysis [24]. The following subchapters elaborate about the numerical tool and the detail of specimens.

#### *2.1. Simulation Method of RC Frame Analysis*

Figure 3a,b shows the interface of the STERA\_3D for the test specimens of towered setback and stepped setback respectively. The RC beam element is modeled as a line element with two nonlinear flexural springs at both ends and one nonlinear shear spring at the middle as shown in Figure 4. The end displacement vector is obtained from Equation (11) as the sum of the displacement vector of each component. The dynamic earthquake response analysis is conducted by adopting Newmark-β numerical integration method [25].

$$\left\{ \begin{array}{c} \theta\_A\\ \theta\_B\\ \delta\_x \end{array} \right\} = \left\{ \begin{array}{c} \tau\_A\\ \tau\_B\\ \delta\_x \end{array} \right\} + \left\{ \begin{array}{c} \varphi\_A\\ \varphi\_B\\ 0 \end{array} \right\} + \left\{ \begin{array}{c} \eta\_A\\ \eta\_B\\ 0 \end{array} \right\} \tag{11}$$

where:

*θ*: the total rotation at the element joint, *δx*: is the element deformation at direction x,


**Figure 3.** STERA\_3D interface on (**a**) a towered setback and (**b**) a stepped setback building models.

**Figure 4.** Elastic, nonlinear bending, and nonlinear shear springs for elements modeled by STERA\_3D.

#### *2.2. Simulation of Experimental Test of Towered Setback Structure*

#### 2.2.1. Specimen Description

The specimen tested by Wood [18] represents a scaled building with three bays and nine stories where the first story has a slightly higher elevation, and the setback location is in between level 2 and 3. The size of the experiment test frame model was 15 times smaller than the real building. The structure consists of two 2D beam-column frames, where the mass is attached to both frames using additional structures. For the first and second story, the weight is 5.04 kN, while the rest are 1.73 kN per level. The concrete compressive strength was 42.33 MPa, and the beam and column longitudinal rebar yielding strengths are 380.59 MPa and 388.17 MPa, respectively. The elevation of the specimen is presented in Figure 5.

**Figure 5.** Experimental test specimen of towered setback structure: (**a**) front view, and (**b**) side view.

#### 2.2.2. Input Data

For specimen of Wood [18], the input motion is only one direction parallel to the 2D frames. The input motions were a set of successive earthquake acceleration records, based on scaled 1940 El Centro NS with different magnitude. In this study, the first input of scaled 1940 El Centro NS with the peak acceleration of 382.8 cm/s<sup>2</sup> was used in the numerical study. The detail of input motion properties and the earthquake time history are presented in Table 2 and Figure 6.

**Figure 6.** Time history acceleration of input motion.

**Table 2.** Detail of input motion properties.


#### 2.2.3. Results of Comparison

The maximum story drift and the maximum acceleration response of each story are shown in Figure 7a,b, and the time history of top floor displacement response is presented in Figure 8. Although there are some differences between analytical and experimental results, the trends of both displacement and acceleration responses are similar between simulation

and experiment. The demonstrated results show that STERA\_3D software has the sufficient accuracy to conduct numerical dynamic analyses of towered setback structures.

**Figure 7.** Comparison of experiment of towered setback structure with STERA\_3D response: (**a**) interstory drift; (**b**) maximum acceleration.

**Figure 8.** Roof displacement history responses of experiment of towered setback structure compared with STERA\_3D.

#### *2.3. Simulation of Experimental Test of Stepped Setback Structure*

#### 2.3.1. Specimen Description

The specimen tested by Shahrooz and Moehle [19] was designed for combined gravity and seismic effects according to the 1982 Uniform Building Code requirements to satisfy the seismic provisions of the American Concrete Institute Building Code (ACI 318-83) and the ACI-ASCE Committee 352 recommendations. The prototype is represented by the 1/4 scale model as depicted in Figure 9, in three dimensions with 2 2 bays of 1143 mm × 1905 mm floor panel size. Design concrete strength was 27.5 MPa, and all reinforcement was Grade 60 (minimum yield stress of 413 MPa). The weight of structures is 72 kN per level for first to third story and 41 kN per level for fourth to sixth story, and the inter story height is 914.4 mm.

#### 2.3.2. Input Data

The specimen of Shahrooz and Moehle [19] was subjected to two input motions sequentially: unidirectional and bidirectional. The first stage consists of three different scaled 1940 El Centro NS ground motions, while the second phase includes three unidirectional horizontal motions inputted at an angle of 45 degrees to simulate bidirectional motions. In this study, the numerical analysis only focuses on the longitudinal study and neglects the bidirectional simulation. The detail of input motion properties and the earthquake time histories are presented in Table 3 and Figure 10.

**Figure 9.** Elevation view of setback experimental test specimen: (**a**) front view, (**b**) side left view, and (**c**) side right view.

**Figure 10.** Time history acceleration of input motions.



#### 2.3.3. Results of Comparison

The relative maximum displacement response and the absolute maximum acceleration response of unidirectional experimental tests by Shahrooz and Moehle [19] and analytical tests are depicted in Figure 11a,b, and the time history of displacement response at the top floor against the last input motion is presented in Figure 12. The result shows that there are only slight differences between experiment and simulation in both acceleration and displacement response.

**Figure 11.** Comparison of experiment of stepped setback structure specimen with STERA\_3D response: (**a**) interstory drift; (**b**) maximum acceleration.

**Figure 12.** Roof displacement history responses of experiment of stepped setback structure compared with STERA\_3D.

#### **3. Parametric Study of Generated Setback Frames**

To examine the relationship between the degree of setback and the distribution of damage, a parametric study was conducted using 35 reinforced concrete setback frames, consisting of 20 models with stepped setback type and 15 models with towered setback type with different degree of setback as depicted in Figures 13 and 14. All the frames have six stores with the number of bays ranging from two to five depending on the model. Each of generated models is represented with two numbers at the above frame illustration. The first number is the difference of the number of bays between top and base structures and the second number is the number of stories in the top structure.

The typical plan of the structure is shown in Figure 15. All the frames use one type of beam section, namely B1, and one type of column section, namely C1. The concrete compressive strength is 30 MPa, while the yield strength of both longitudinal and transversal rebars are 390 MPa. The slab thickness is 15 cm with reinforcement of D13 with the spacing of 10 cm. The section details of elements are presented in Table 4. The total gravity load implemented for every floor is 12 kN/m2.

**Figure 13.** Generic 2D frame structures of stepped setback model.

**Figure 14.** Generic frame structures of towered setback model.



**Figure 15.** 2D plan of generated structure of 2 × 5 bays for 6-story building.

Nonlinear dynamic analyses were performed using STERA\_3D, and the story damage indices were populated. The list of input earthquake ground motions is shown in Table 5. Ten input ground motions are recorded acceleration waves scaled to have the maximum velocity of 50 cm/s. The ground motion data are retrieved from PEER [26] and COSMOS Virtual Data Center [27]. The intensity of these input ground motions corresponds to the level of design earthquakes for the safety limit state in the Japanese seismic standard. The acceleration response spectrum is presented in Figure 16, and the wave shapes of the earthquake input motions are depicted in Figure 17.


**Figure 16.** Acceleration response spectra (with 5% damping factor) of input earthquake motions.

**Figure 17.** Time history acceleration waves of input earthquake motions.

#### **4. Proposal of Tower-Base Damage Index Ratio**

The nonlinear dynamic analysis using STERA\_3D was conducted for all 35 frames with three input earthquake ground motions to calculate the damage index (*DI*) of each story. The damage index of each story is calculated using Equation (5). The following *DI*ratio is proposed to examine the distribution of the *DI* along the height of the frame.

$$DI\_{ratio} = \frac{\max(DI\_{tower})}{\max(DI\_{base})} \tag{12}$$

where the *DItower* and *DIbase* are taken as the average for the ten input ground motions. If the value of *DIratio* is larger than one, the damage will be concentrated at the tower structure. Conversely, if the value is less than one, the base structure will suffer more damage than the tower.

The mathematical model is developed for predicting *DIratio* from irregularity indices *ϕ<sup>s</sup>* and *ϕ<sup>b</sup>* using the nonlinear regression analysis so that the damage distribution is easily evaluated without conducting the nonlinear dynamic analysis. The relation between the calculated *DIratio* of all the frames and irregularity indices *ϕ<sup>s</sup>* and *ϕ<sup>b</sup>* are shown in Figure 18 for the stepped model and Figure 19 for the towered model type, respectively.

**Figure 18.** *DIratio* with different *ϕ<sup>s</sup>* and *ϕ<sup>b</sup>* for stepped model.

**Figure 19.** *DIratio* with different *ϕ<sup>s</sup>* and *ϕ<sup>b</sup>* for towered model.

Values of *ϕ<sup>s</sup>* and *ϕ<sup>b</sup>* equal to 1 correspond to a structure with a uniform configuration. As values of *ϕ<sup>s</sup>* and *ϕ<sup>b</sup>* increase, the degree of the setback in the frame structure becomes larger. A large value of *DIratio* indicates a relatively large damage concentration in the tower. The nonlinear regression analysis is performed to derive the quadratic polynomial equation for estimating the damage index ratios from irregularity indices,*ϕ<sup>s</sup>* and *ϕb*. The obtained equations are shown in Equation (13) for stepped type and Equation (14) for towered type. The red lines in Figures 18 and 19 indicate the prediction of the damage index ratio using Equations (13) and (14). The values of R-squared of red lines, which indicate the goodness-of-fit measure for the regression model, vary between 0.865 and 0.999.

$$DI\_{r, \text{hyperpl}} = 98.19 - 124.29q\_s - 105.06q\_b + 142.75q\_sq\_b + 38.34q\_s^2 + 2547q\_b^2 - 37.09q\_sq\_b^2 - 45.44q\_s^2q\_b + 13.22q\_s^2q\_b^2 \tag{13}$$

$$DI\_{t\text{forward}} = -3.09 + 16.89q\_s + 17.22q\_b - 35.99q\_sq\_b - 12.88q\_s^2 - 12.53q\_b^2 + 17.34q\_sq\_b^2 + 19.59q\_s^2q\_b - 6.22q\_s^2q\_b^2 \tag{14}$$

#### **5. Validation of Proposed Damage Index Ratio**

To validate the proposed method to calculate *DIratio*, the towered setback frame by Wood [18] and the stepped setback frame tested by Shahrooz and Moehle [19] are examined. The input wave is the 1940 El Centro NS record scaled to have the maximum velocity of 50 cm/s; however, the time step is adapted for each specimen. The specimen of Wood [18] is subjected to the earthquake with the time step divided by 4 by rounding up the value of <sup>√</sup>15, while the specimen of Shahrooz and Moehle [19] is subjected to the earthquake with the time step divided by 2, by considering the scale factor of test specimen. Dynamic time history analysis was performed using STERA\_3D to calculate the *DI* of each story. The calculated damage index ratios were compared with the predicted damage index ratios from the proposed equations and presented in Table 6.

**Table 6.** Numerical analysis result compared to prediction result using Regression Equation.


The accuracy of the prediction from Equation (13) is adequate, 66% for the stepped setback and the prediction accuracy from Equation (14) is relatively high, 93% for the towered setback.

#### **6. Conclusions**

To understand better the earthquake response and damage concentration of setback buildings, numerical investigations were conducted and a practical method to evaluate the damage distribution was proposed.

First, two previous experimental tests of reinforced concrete buildings comprising a towered setback and a stepped setback were modeled and analyzed by STERA\_3D and the test results were successfully simulated.

Then, to study the building damage distribution, 2D reinforced concrete frame models containing 20 stepped and 15 towered setback frames were generated to inspect the relationship between the damage index ratio and the irregularity indices under safety level earthquake input ground motions. Based on the results, two nonlinear regression equations were proposed as an alternative to the dynamic analysis procedures.

Finally, the proposed equations were applied to the previous experimental tests, and it was demonstrated that the proposed formulas could predict the damage to vertically irregular reinforced concrete moment-resisting frame with reasonable accuracy. Although the prediction formula works, it is important to be noted that the setback buildings considered in this study are limited for low- to mid-rise reinforced concrete frame buildings with uniaxial seismic loading.

**Author Contributions:** Conceptualization, T.S., B.E. and T.I.M.; methodology, B.E., T.I.M.; software, T.S.; writing—original draft preparation, T.I.M.; writing—review and editing, T.S.; supervision, T.S. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research received no external funding.

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** The data presented in this study are available on request from the corresponding author.

**Acknowledgments:** The first author appreciates AUN/SEED-Net JICA for the support through scholarship of doctoral degree.

**Conflicts of Interest:** The authors declare no conflict of interest.

#### **References**

