Abstract

The European seismic code 8 (Eurocode 8) classifies buildings as planwise regular according to four criteria which are mostly qualitative and a fifth one, which is based on parameters such as stiffness, eccentricity, and torsional radius, that can be only approximately defined for multistory buildings. Therefore, such plan-regularity criteria are in need of improvement. ASCE seismic code, according to a different criterion, considers plan (or “torsional”) irregularity in a building when the maximum story drift, at one end of the structure, exceeds more than 1.2 times the average of the story drifts at the two ends of the structure under equivalent static analysis. Nevertheless, both the ASCE approach and the threshold value of 1.2 need to be supported by adequate background studies, based also on nonlinear seismic analysis. In this paper, a numerical analysis is carried out, by studying the seismic response of an existing R/C school building taken as the reference structure. Linear static analysis is developed by progressively shifting the centre of mass, until the ratio between the maximum lateral displacement of the floor at the level is considered and the average of the horizontal displacements at extreme positions of the floor at the same level matches and even exceeds the value of 1.2. Then, nonlinear dynamic analyses are carried out to check the corresponding level of response irregularity in terms of uneven plan distribution of deformation and displacement demands and performance parameters. The above comparison leads to check the suitability of the ASCE approach and, in particular, of the threshold value of 1.2 for identifying buildings plan irregularity.

1. Introduction

In-plan irregularity is very common in existing buildings, and it is one of the most frequent sources of severe damage during earthquakes [1]. In 1985, during the earthquake of Mexico City, 42% of damaged or collapsed structures were corner buildings. Furthermore, many buildings failed in torsion due to the asymmetric layout of masonry infills [2]. The need to define adequate provisions accounting for the torsional effect due to structural asymmetry was reconfirmed by this event [3].

During the 1970s–1980s, the equivalent static analysis represented the most common method for computing seismic loads. Such a method, applied to asymmetric structures, underestimates the actions on flexible side elements, as it does not take into account the dynamic amplification of the torsional response. For this reason, researchers suggested the introduction of a design eccentricity in order to provide a torsional moment corresponding to each story of the building. In the last decades, several studies on the design eccentricity were carried out [312]. Later, spatial models and modal analysis became widely used by designers. Thus, the elastic response determination of in-plan irregular buildings turned easier.

In recent years, researchers focused on the inelastic response of in-plan irregular buildings. The response of asymmetric buildings was investigated by varying several parameters like centre of mass (CM), eccentricity [13], uneven distribution of concrete strength [14, 15], torsional stiffness, and periods of vibration [16]. Furthermore, several studies were developed focusing on the application of nonlinear static analysis to in-plan irregular buildings [1722]. Then, an important contribution that deals with torsion issue and plan regularity of building structures is provided by Anagnostopoulos et al. [23].

This paper is focused on the evaluation of the ASCE torsional provision (ASCE Standard 7-10) [24]. The response of a case study building is investigated in order to evaluate the torsional irregularity provided by the code and the corresponding uneven distribution of inelastic demand which is detected using the nonlinear dynamic analysis method.

2. Torsional Provisions

The previous generation of seismic codes was providing design eccentricities for equivalent static analysis as follows:where is the design eccentricity, is the distance between the centre of mass (CM) and the centre of rigidity (CR) measured orthogonally to the loading direction, is the building dimension perpendicular to the loading direction, is a coefficient accounting for dynamic amplification of the torsional response, and is a coefficient accounting for aleatoric position of the CM.

In Table 1, coefficients and are provided according to Eurocode 8 1993 (EC8-93) [25], National Building Code of Canada 1995 (NBCC-95) [26], and ASCE Standard 7-10 (ASCE 7-10) [24], where is the amplification factor provided by ASCE 7-10 (defined in the following) and is the additional eccentricity provided by EC8-93, equal to the smaller of the following values:where is the radius of gyration of the floor mass in-plan (square root of the ratio of the polar moment of inertia of the floor mass in-plan (a) with respect to the centre of mass of the floor to (b) the floor mass (b)) and is the square root of the ratio of the torsional stiffness to the lateral stiffness in the loading direction (“torsional radius”).

In current code provisions, such as Eurocode 8 2004 (EC8-04) [27] and National Building Code of Canada 2010 (NBCC-10) [28], in the case of torsional irregularity, only a 3D dynamic analysis is allowed. Even when adopting a spatial model and the CQC modal combination, the accidental eccentricity is considered accounting for the random position of CM, and the design eccentricity is computed by the following equation:

However, criteria for assessing in-plan irregularity are still very important; indeed, it affects the choice of the method of analysis and the definition of the behaviour factor.

Qualitatively, in-plan irregularity is due to asymmetric distributions of mass and stiffness, but quantitatively, a definition universally shared does not exist. Eurocode 8 2004 (EC8-04) provides a list of conditions to classify a building as regular in-plan:(i)With respect to the lateral stiffness and mass distribution, the building structure shall be approximately symmetrical in plan with respect to two orthogonal axes.(ii)The plan configuration shall be compact, i.e., each floor shall be delimited by a polygonal convex line. If in-plan setbacks (i.e., reentrant corners or edge recesses) exist, regularity in plan may still be considered as being satisfied provided that these setbacks do not affect the floor in-plan stiffness and that for each setback, the area between the outline of the floor and a convex polygonal line enveloping the floor does not exceed 5% of the floor area.(iii)The in-plan stiffness of the floors shall be sufficiently large in comparison with the lateral stiffness of the vertical structural elements, so that the deformation of the floor shall have a small effect on the distribution of the forces among the vertical structural elements. In this respect, the L, C, H, I, and X plan shapes should be carefully examined, notably as concerns the stiffness of the lateral branches, which should be comparable to that of the central part, in order to satisfy the rigid diaphragm condition. The application of this paragraph should be considered for the global behaviour of the building.(iv)The slenderness λ = Lmax/Lmin of the building in-plan shall not be higher than 4, where Lmax and Lmin are, respectively, the larger in-plan and smaller dimension of the building, measured in two orthogonal directions.(v)At each level and for each direction of analysis, the structural eccentricity e0 and the torsional radius shall conform to the following conditions:

The first four criteria are almost qualitative, while the fifth one, if strictly applied, is valid only for single-story buildings. For multistory buildings only, approximate definitions of the centre of stiffness and the torsional radius are possible [27]. Similarly, ASCE Standard 7-10 [24] provides a list of conditions to detect horizontal irregularity in buildings, as follows:(i)Torsional irregularity: it is defined to exist if the maximum story drift, computed including accidental torsion with , at one end of the structure transverse to an axis, is more than 1.2 times the average of the story drifts at the two ends of the structure. Torsional irregularity (equation (5)) requirements apply only to structures in which the diaphragms are rigid or semirigid:where is the average deflection determined by an elastic analysis (Figure 1) and is the amplification factor of the accidental torsional moment.

(ii)Extreme torsional irregularity: it is defined to exist if the maximum story drift, computed including accidental torsion with  = 1.0, at one end of the structure transverse to an axis, is more than 1.4 times the average of the story drifts at the two ends of the structure. Extreme torsional irregularity requirements (equation (6)) apply only to structures in which the diaphragms are rigid or semirigid:(iii)Reentrant corner irregularity: it is defined to exist where both plan projections of the structure beyond a reentrant corner are greater than 15% of the plan dimension of the structure in the given direction.(iv)Diaphragm discontinuity irregularity: it is defined to exist where there is a diaphragm with an abrupt discontinuity or variation in stiffness, including one having a cut-out or open area greater than 50% of the gross enclosed diaphragm area or a change in effective diaphragm stiffness of more than 50% from one story to the next.(v)Out-of-plane offset irregularity: it is defined to exist where there is a discontinuity in a lateral force-resistance path, such as an out-of-plane offset of at least one of the vertical elements.(vi)Nonparallel system irregularity: it is defined to exist where vertical lateral force-resisting elements are not parallel to the major orthogonal axes of the seismic force-resisting system.

The ASCE Code (ASCE 7-10) [24] provides a simpler approach for detecting torsional irregularity. Such a method involves directly the structural response of the building, and it does not require the knowledge of specific characteristics (that are mandatory in Eurocode 8), such as the torsional radius, that can be only approximately defined for multistory buildings. The criterion provides the implementation of a linear static analysis taking into account the accidental eccentricity (β = 5%). As shown previously, when the maximum story drift is more than 1.2 times the average drift, the torsional irregularity is detected. The use of this elastic index (i.e., the ratio between the maximum and the average elastic drift) appears to be very simple and effective; however, some further researches are needed. Both the ASCE approach and the threshold value of 1.2 must be supported by adequate background studies.

With this aim, in this paper, a numerical analysis is carried out, by studying the seismic response of an existing R/C school building. As provided by ASCE 7-10, static elastic analyses are performed varying, step by step, the CM position and checking the ratio between the maximum and the average lateral displacements (δmax/δavg). Lastly, nonlinear dynamic analyses are carried out in order to evaluate the corresponding level of response irregularity in terms of uneven distribution of inelastic demand.

3. Case Study

The school building under analysis is about 40 years of age, and it is situated on a flat ground in the Municipality of Prato (PO), Italy. It consists of two independent blocks: the school and the gym.

The school block has been chosen as the case study for the seismic analyses because it is made of a R/C framed structure (conversely, the gym block has a precast structure).

The school block (Figure 2) has an elongated rectangular plan, approximately 49 by 13 meters, and it is a three-story building above the ground level. The floor-to-floor height is 3.30 meters. At the ground floor, along the shorter side of the plan, there is a portico that covers the first span of the frames and that provides access to the entrance hall. After the entrance hall, a central corridor leads to the classes and offices along the two long sides of the plan. This distribution pattern is repeated, almost unchanged, at the upper floors as well. The external cladding, in precast panels, and the interior walls, which divide the classes, are made of hollow bricks and follow the frame scheme, and they also have ribbon windows outwards.

The building structure (Figure 3) is entirely made of R/C columns and beams, which form only three plane moment-resisting frames (MRF), running along the longitudinal direction. In the transversal direction, there are not MRF, excepting for the frames at the ends of the structure. The floor slabs can be considered as the rigid diaphragms.

The structural symmetry of the plan fails, in particular, at the ground floor where, due to the presence of a larger room devoted to refectory, in the central longitudinal frame, the second column is missing. Therefore, the frame beam has a span doubled and supports the column of the upper floor at midspan, and the rectangular beam cross section is 25 by 90 cm, while the others are flat beams with a cross section of 80 by 22 cm. Cross section 25 by 90 cm was also used for the corresponding beams at the two upper floors. The story floors are made of casted-in-place one-way hollow block slabs.

Instead, the rafters and the edge beams supporting the attic floor have a section of 12 by 40 cm and form a hip roof. The covering is made of insulated sandwich panels. The attic floor is not serviceable.

The columns have three different cross sections: 25 by 30 cm, 25 by 40 cm, and 25 by 55 cm. In the two longitudinal edge frames, the longer section side is arranged orthogonally to the frames, while in the central frame, the longer side follows the frame direction. In Table 2, the characteristics of the column sections at the ground floor are provided.

The concrete compressive strength was assumed to be equal to fc = 25 MPa, and the steel yield stress was assumed to be fy = 430 MPa.

A finite element structural model has been created, and numerical analyses have been carried out using the SAP2000 software [29]. In addition to the dead weight of the structural elements, the following loads have been considered: the weights of the story and roof floors (), the permanent loads (), such as the partition walls, the screed, and the tiles, and the live loads (). In the following, the seismic combination (7) for the ultimate limit states is provided according to the Italian Code (NTC 2018) [30]:where is the seismic action along the considered direction. The seismic base shear has been computed in accordance with equation (8), provided by the NTC 2018:where is the design spectral acceleration calculated according to the fundamental period of the building (as in this case, the elastic response spectrum has been adopted, and is equal to ), is the fundamental period of the building, is the effective seismic weight, is the gravity acceleration, and factor equal to 0.85 in the case of a multistory building (with the story number greater than or equal to 3).

The elastic response spectrum is used to calculate the seismic input for SLV verifications, assuming a building life of 50 years, and importance of class III (Cu = 1.5) has been obtained considering the building site (Figure 4).

In the structural model, rigid diaphragms have been considered at each building floor, so as to easily define the position centre of masses (CM).

Firstly, the model has been seismically analyzed by the equivalent static method [31], by applying the horizontal action at the CM of each story floor. The seismic force, computed according to the formulae of the Italian NTC 2018 [30], has been assumed acting in the direction of the shorter building side (y direction), where only the two edge frames exist, in order to catch the greatest displacements.

Then, nonlinear dynamic analyses have been performed using a natural earthquake accelerogram (Figure 5).

The adopted accelerogram is the record 147ya of Friuli (aftershock) earthquake happened in 1976, September 15th. The recorded magnitude (Mw) is 6, and the peak ground acceleration (PGA) is equal to 2.32 m/s2.

This accelerogram was chosen because it is spectrum compatible with the adopted elastic response spectrum. In the following, a comparison between the elastic response spectrum from the Italian NTC 2018 and the 147ya spectrum is proposed (Figure 6).

To assess the seismic performance of a building, through a nonlinear dynamic analysis, the Italian NTC 2018 requires, at least, three accelerograms. Therefore, the choice of using only one accelerogram could appear inadequate. Actually, such request aims to reduce the uncertainties of the analysis, and more in general, it aims at increasing the structural safety. Since the assessment of the structural safety of the building is out of the purpose of this paper, that deals, instead, with the influence of the in-plan irregularity on the structural inelastic response, only one accelerogram has been considered. It is evident that also the seismic action, both in terms of intensity and frequency content, has an influence on the response of in-plan irregular buildings; however, the evaluation of this influence is out of the scope of this work.

The nonlinear dynamic analyses were conducted using lumped plasticity. In particular, at both extremities of all beams and columns, potential plastic hinges were placed. The latter were developed according to the ASCE 41-13 [32], and the inelastic demand (D) is expressed in terms of the plastic rotation at the chord (θpl). The capacity (C) is given by the limit life safety (LS), as defined by ASCE 41-13. In Figure 7, the idealized full backbone curve, defined according to ASCE 41-13, is shown.

The first part of the backbone curve (A-B) represents the elastic behaviour of the section, and the activation of the plastic hinge occurs beyond the point B. The peak of the strength is reached in correspondence to the point C (note that the point C does not represent the capacity), and then the strength decreases rapidly until reaching the point D where a small residual strength is still provided by the section. The failure occurs at point E where the strength drops to zero.

Therefore, points A to E describe the nonlinear behaviour of the section, while the acronyms IO, LS, and CP indicate, respectively, the acceptance criteria for immediate occupancy (IO), life safety (LS), and collapse prevention (CP). For both life safety and collapse prevention, the acceptance criteria are differentiated for primary components (P) and secondary components (S). For primary components, LS = 0.75 C and CP = C but not greater than CP = 0.75 E, while for secondary components, LS = 0.75 E and CP = E.

Therefore, the capacity of the section (C) is represented by the plastic component of the LS deformation calculated as C = LS − B.

4. Results and Discussion

In this section, results of numerical analyses are presented and discussed. The gross member flexural rigidity was considered and the concrete elasticity modulus was assumed to be equal to 30 GPa.

In Table 3, the periods of vibration, the percentage of participating mass along the longitudinal direction of the building (x direction), the percentage of the participating mass along the transversal direction of the building (y direction), and the percentage of the mass participating to the torsional mode are listed for each mode of vibration considered.

Referring to the spectral ordinate , associated to the fundamental period of vibration along the y direction, the static forces at each floor level are calculated. Then, equivalent static force analyses are performed in the y direction, by varying the CM position along the x direction. The same procedure is repeated for the response spectrum analysis and the ratios between the maximum and the average lateral displacements (δmax/δavg), calculated with both methods at varying of the CM position, are presented in Figure 8. The effect of variation of the periods of vibration, induced by the variation of CM position, is considered negligible.

Two considerations arise from previous results. When the CM is placed in the actual position, the floor displacement includes rotational components (indeed δmax/δavg ≠ 1): this means that CM and CR are not aligned (CR is detected when δmax/δavg = 1). The two methods, the equivalent static force and the response spectrum analysis, give almost the same results in term of displacements when CM is close to CR. However, increasing the distance between CM and CR, the agreement between the results gets worse, but this difference is in accordance with several researches [11, 12] and it is widely discussed in Torsional Provisions. In the present case, according to ASCE Code (ASCE 7-10) and considering the 5% of accidental eccentricity, the building resulted in regular in-plan since the ratio between maximum and average displacement resulted lower than δmax/δavg < 1.2 (Figure 8).

Then nonlinear dynamic analyses are performed in order to assess the uneven distribution of inelastic response at varying CM position. In particular, four analyses are performed: the first considers CM coincident with CR, the second considers CM in its actual position, and the third and fourth consider CM in the positions that give δmax/δavg equal to 1.2 and 1.4, respectively, in the equivalent static force analysis (Figure 8).

In Figure 9, the ratios between inelastic demand D and capacity C (D/C), calculated by means of nonlinear dynamic analyses, of all columns at the ground floor, are shown in order to evaluate the damage levels. As previously said, in the y direction, the only frames are placed on the left and right side of the building. All the other columns, placed between these two frames, have to be considered as cantilever since such columns are not connected by beams in the y direction. This fact explains why the columns placed on the left and right side of the building are much more affected by the inelastic demand than the others.

Furthermore, the columns placed in the centerline have the longer side of the cross section oriented along the x direction, instead the columns placed on the first and third line have the longer side of the cross section oriented along the y direction. For this reason, the inelastic demand of the centerline columns is lower than that of the columns placed on the first and third lines, due to the greater shear force brought by the latter. Although the uneven distribution of stiffness, the demand to capacity ratios are symmetrically distributed, excepted for the abrupt variation shown at the extremities, such ratios are almost constant.

When the CM is placed in its actual position, the scenario is different because the symmetrical distribution of the ratios D/C is lost. The inelastic demands on the flexible side are greater than those shown on the stiff side. Moving from the stiff to the flexible side, the inelastic demand increases progressively, as shown in Figure 10.

As previously shown (Figure 8), assuming an eccentricity of 6% for the CM, the ratio δmax/δavg, obtained by performing an equivalent linear static analysis becomes equal to 1.2. In this configuration, the inelastic demand appears clearly irregular. In approximately 40% of the building (Figure 11), the ratios D/C are almost zero, while in the other part, the inelastic demand increases rapidly moving to the flexible side.

When CM is placed with 14% of eccentricity, the ratio δmax/δavg, calculated by an equivalent linear static analysis, is equal to 1.4 (Figure 8). Adopting this configuration, the inelastic response is strongly irregular. As shown in Figure 12, approximately 50% of the building provides no inelastic demand. Conversely, moving to the flexible side, the plastic demand increases rapidly providing high values of the ratio D/C.

In Figure 13, the demand to capacity ratios, calculated by means of nonlinear dynamic analyses, of the three columns on stiff and flexible side (respectively, left columns of frame Y1 and right columns of frame Y2 in Figure 3(b)) are plotted at varying δmax/δavg calculated by the equivalent linear static analysis. In this case, it is worth noting that the inelastic demand is highly dependent on δmax/δavg and that the demand decreases on the stiff side and increases on the flexible side.

In Figure 14, the demand to capacity ratios are normalized by D/C calculated in configuration CM≡CR, where δmax/δavg is equal to 1. In correspondence of δmax/δavg = 1.2, the normalized D/C on the stiff side is about equal to zero, and on the flexible side, it is close to 2.

To investigate the scattering of the damage levels over the building plan, the maximum (D/C)max and average (D/C)avg demand to capacity ratios of all columns at ground floor were computed in the four studied cases (Figures 912).

In Figure 15, the inelastic indices (D/C)max/(D/C)avg were normalized by the corresponding (D/C)max/(D/C)avg computed in the configuration in which CM≡CR (i.e., δmax/δavg = 1).

Figure 15 clearly shows that the value δmax/δavg = 1.2 may be considered as a threshold value beyond which plan scattering of the damage level D/C does not increase significantly. Of course, this result may be affected by the particular condition of the analyzed building, and therefore, it needs to be validated by means of further investigations on different building types.

5. Conclusions

This paper is focused on the evaluation of the ASCE code approach for detecting building plan irregularity. The response of an existing R/C school building, taken as reference case study, is investigated to evaluate how well torsional irregularity detected by the code corresponds to uneven distribution of inelastic demand found by nonlinear dynamic analyses. In particular, for the studied case, the threshold value of 1.2 for identifying buildings plan irregularity seems to be appropriate when plan scattering of damage levels is considered. In this specific case, the value of 1.2 seems already to characterize an extreme torsional irregularity condition, and therefore, the value 1.4 provided by the ASCE code may be too large. However, it is worth noting that this trend could also arise from the specificity of the building under investigation, and therefore, a larger number of existing buildings with different structural arrangements and plan configurations are currently under analysis to provide support to this thesis.

Data Availability

The authors state that data in the manuscript are freely available.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

The financial support of ReLUIS through the Progetto Esecutivo Convenzione DPC/ReLUIS 2017—PR2 Strutture in cemento armato—is gratefully acknowledged.