Geometric characteristics of elements from each modeled building
1. Introduction
Along its history, Venezuela has been severely affected by destructive earthquakes [1]. Approximately 80% of the population lives in seismically active areas, where have occurred destructive earthquakes even in recent times [2]; The seismic hazard, inadequate design and construction of buildings as well as the damage occurred from previous earthquakes, demonstrate a high vulnerability in existing buildings. Then it is essential to continuously make progress and research in the field of earthquake engineering and upgrade the seismic design codes. Seismic upgrade requires the evaluation or predictions of the expected damage to structures at the time of an earthquake of a certain severity occur. From this prediction it can be defined solutions for the reduction of structural vulnerability [3].
The damage occurred in buildings after an earthquake indicates the need for reliable methodologies for the evaluation of seismic behavior of the existing buildings. According to current technical and scientific advances, seismic evaluation of reinforced concrete (RC) structures can be done by two different approaches: empirical methods and mechanical methods [4]. The current tendency of earthquake engineering in the evaluation of structural behavior is the application of simplified mechanical methods based on performance, involving the capacity spectrum [5], because there are developed refined models and detailed analysis.
This study used a mechanical method that involves non-linear analysis with deterministic and probabilistic approaches, as well as procedures of analysis based on Limits States defined by displacements [6], in order to evaluate the behavior of a low rise RC building with plan irregularity, designed according to Venezuelan codes [7]-[9] and subjected to seismic action effect. Through the use of mathematical models and computational tools, seismic behavior of the building is obtained in a suitable way. Among these tools any procedure was chosen: the quadrants method, which leads to the rapid assessment of the seismic capacity of a structure through its non-linear response [10]. Results of the research shown that the current design of this kind of structures is not safe when they are under the maximum seismic actions prescribed by codes, then it is necessary to review the design procedures in order to find more realistic designs that fulfill the goals of the performance-based design.
2. Case studied
A two story RC framed building was studied, (Figure 1a), which contains internal staircase and 220 m2 total plan area. This structure represents a common typology used for residential buildings in Venezuela, for prone seismic zones. The structure was designed and detailed for a high ductility value (response reduction factor of 6).
The building was modeled according its original design, called
|
|
|
||
OB | 20x35 | 20x35 | 20x30 | 20x30 |
RB | 20x45 | 20x35 | 30x30 | 30x30 |
DBDB | 20x40 | 20x40 | 35x35 | 30x30 |
3. Assessment method
The Quadrants Method is based on the results of the non-linear static analysis (Pushover analysis). This analysis results are plotted in a displacement vs. base shear format, this generate the capacity curve which represents the overall capacity of the whole structure against lateral forces. In order to evaluate the capacity curve two of the main structural parameters are taken into account. The first one is the design elastic shear, obtained from the elastic analysis of the structure using the elastic design spectrum. The second parameter is the threshold that defines the Repairable Limit State, obtained from [14] for RC framed buildings with similar characteristics to the studied ones. The thresholds have been computed from characteristic values of three levels of damage proposed in [15] and are showed in Table 2. Both values are used to define two axes over the capacity curve, the elastic base shear defines an horizontal axis and the damage threshold defines a vertical axis, then Capacity Curve is divided in four spaces or quadrants, see Figure 2.
|
|
|
|
|
Service | Non-structural | Frequent | 50% in 50 years | 0,2<δ<0,5 |
Damage control | Moderate structural | Occasional | 10% in 50 years | 0,5<δ<1,5 |
Collapse prevention | Severe structural | Rare | 2% in 50 years | 1,5<δ<3,0 |
The performance point is a common procedure accepted among the scientific community to evaluate the seismic performance of a structure under a specific demand. It is usually obtained from the idealized shape of the Capacity Curve as is shown in the Figure 2 [16]. The Quadrants Method also uses this parameter in order to define the roof displacement of the case studied, defined according to the N2 method [5]. If the performance point is under the axis defined by the elastic base shear (Quadrants III or IV), the design does not meet the basic objective of the seismic design because the building does not have enough lateral strength. If the performance point is on the right side of the vertical axis (Quadrant I) means that the building has adequate stiffness, otherwise (Quadrant II) it means that the stiffness is very low and the displacements can be longer than the displacements that can produce advanced structural damage, technically or economically irreparable. These lateral displacements are usually computed from the dynamic response of the structure submitted to a strong motion with a return period of 475 years, or an occurrence probability of 10% in 50 years [17].
The Quadrants Method can provide an objective criterion in order to upgrade the seismic capacity of a structure. If the performance point is on the Quadrant I, the structure has enough lateral strength and stiffness, so does not need to be reinforced. If the structure is on the Quadrant II, it is necessary to provide additional stiffness by using conventional procedures like RC or steel jacketing. If the performance point is on Quadrant III, the structure requires a more radical intervention, adding stiffness and lateral strength. In this case it is possible to combine some traditional reinforcement techniques with new ones like FRP jacketing. In this case the columns are the subject of the main intervention. Finally, if the performance point is on the Quadrant IV, the structure does not has enough lateral strength and then the reinforcement technique must be FRP jacketing.
4. Nonlinear analysis
The structures are modeled by incorporating the structural response when it incurs in the material and geometrical non-linear range, produced by high deformations caused by accidental excitations (earthquakes) [11]. The analyses were performed using ZEUS-NL software [18], which allows to model complex structures with "n" number of finite elements, thus to know the elements in the building which are most vulnerable to damage. Each building is modeled in two dimensions, spitting each frame to get a more detailed response for the seismic behavior of each frame; a 3D dynamic analysis was applied to the ER model.
The static Pushover analysis is performed once the frames have been subjected to action of gravity loads, based on the pseudo-static application of lateral forces equivalent to displacements of seismic action [5]. The pattern of lateral seismic loads consist in increasing loads with height (triangular distribution) applied in a monotonic way until the structure reaches its maximum capacity [20].
This procedure applies a solution of equilibrium equations in an incremental iterative process form. In small increments of linear loads, equilibrium is expressed as:
Where
While this procedure is applied, the strength of the structure is evaluated from it is balance internal conditions, updating at each step the tangent stiffness matrix. Unbalanced loads are applied again until it can satisfy a convergence criterion. Then, a new load increase is applied. The increases are applied until a predetermined displacement is reached or until the solution diverges.
From the capacity curve provided in this analysis, it is determined the structural ductility (μ) by the quotient between the ultimate displacement and cadence point displacement, as shown in the following expression:
Where ∆u is Ultimate displacement and ∆y is the global yield displacement. Both values are computed from the idealized capacity curve of the structure.
By the other hand, the dynamic analysis is an analysis method that can be used to estimate structural capacity under seismic loads. It provides continuous response of the structural system from elastic range until it reaches collapse. In this method the structure is subjected to one or more seismic records scaled to intensity levels that increase progressively. The maximum values of response are plotted against the intensity of seismic signal [21-22]. The procedure to perform the dynamic analysis from the seismic signal is:
To define a seismic signal compatible with the design scenario;
To define the scaled earthquake intensity a monotonic way;
To define the extent of damage or damage Limit States;
To study a seismic record for the dynamic analysis of a structural model parameterized to measure earthquake intensity;
The non-linear dynamic analyses provide a set of curves which are a graphical representation of the evolution of the drifts respect time. Results let to compute the damage lumped in specific elements of the structure, but these results are beyond the objective of this Chapter.
|
|
|
|
|
Frequent | Serviceability | 95 | 50 % | δ < 0,5 |
Rare | Reparable damage | 475 | 10 % | δ < 1,5 |
Very Rare | Collapse Prevention |
2475 | 2 % | δ < 3,0 |
For the dynamic analysis the structures were subjected to seismic action (see Table 2) defined by accelerograms built on the basis of a likely value of maximum acceleration of the soil and the hazard level associated with the location of the structure and other seismic characteristic design parameters [16]. These accelerograms called "synthetic accelerograms" are generated through the implementation of a set of earthquakes with wide frequency content, using the PACED program [17], based on the Venezuelan code’s elastic design spectrum. For the dynamic analyses of the three buildings (
Non-linear dynamic analysis was applied to all buildings in order to verify if the performance evaluated by the Quadrants Method is reliable in order to evaluate the fulfilment of the thresholds defined in the precedent section. For this purpose they has been computed three synthetic elastic design spectrum-compatible accelerograms by means of the PACED program [23]. In Figure 3 are shown the Venezuelan rigid-soil elastic design spectrum with the response spectra obtained from the synthetic accelerograms.
These three earthquakes were applied to all frames from the three buildings evaluated, in order to obtain maximum displacement that can be reached by each one. In the software used [18], it was required the implementation of dynamic loads in direction X and the assignation of a control node located in the gravity center of the roof level.
The 3D non-linear dynamic analysis is based on the procedure explained in [20]. The RB building is analyzed, defining its geometry, materials and sections, serviceability loads in Y direction in all beams-columns joints, and dynamic loads on outer nodes with directions and combinations shown in Table 4. One direction ribbed slabs were modeled as rigid diaphragms in its plane by using additional elements with no flexural capacity (Figure 4).
Once built the model, there were applied all the accelerograms with the combinations shown in Table 3, for the interstorey drifts and maximum torsional moments on supports. These combinations are based on the Venezuelan seismic code [7] and following established by [24] about the seismic response of asymmetric structural systems in the inelastic range.
|
|
1 | 100 % (X) |
2 | 100 % (Z) |
3 | 100% (X) y 30% (Z) |
4 | 100% (Z) y 30% (X) |
5. Results
From classic elastic analysis, the verification of interstorey drifts of the
To determine the values of structural ductility it was necessary to plot the idealized curve in function of the capacity curve obtained from non-linear pseudo-static analysis (pushover analysis), in order to know the point at which the structure begins to yield. Figure 5 shows an example of the normalized capacity curve with the idealized (bi-linear) curve of Frame C of
Structural Ductility |
|
|
||
EO | ER | DBDB | ||
|
5,56 | 5,52 | 4,77 | |
|
2,22 | 6,04 | 5,38 | |
|
2,17 | 4,69 | 5,25 | |
|
2,21 | 5,54 | 5,59 | |
|
2,23 | 7,07 | 6,06 | |
|
2,66 | 5,29 | 6,69 | |
|
2,20 | 4,17 | 5,92 | |
|
2,83 | 5,95 | 6,24 |
From the obtained capacity curves there were computed the Performance point (Pp) of every frame of each evaluated building. Table 5 presents the values of Pp of all the frames of evaluated buildings. Figure 6 shows the determination of the Pp of Frame C from
|
|
||
EO | ER |
|
|
A | 5,94 | 2,42 | 2,52 |
B | 13,89 | 9,47 | 7,43 |
C | 15,22 | 9,50 | 9,38 |
D | 14,01 | 9,50 | 7,57 |
E | 13,45 | 9,55 | 6,60 |
1 | 12,62 | 9,35 | 6,07 |
2 | 15,74 | 11,48 | 9,29 |
3 | 10,92 | 7,57 | 4,23 |
From dynamic analyses, there were determined global and interstorey drifts of each frame from all three models studied. Both types of drifts were calculated on the basis of the application of synthetic accelerograms with different intensities, representing the lateral forces applied to frames in order to generate their respective maximum displacements. Figures 7 to 9 show the evolution of the global (∆/H) drifts expressed as a percentage, respect to time (sec) of the frame C from
Figures 10 to 12 show the results for interstorey drifts of frame C from
|
|
||||||||
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
A | - | - | - | - | OK | OK | OK | OK | OK |
B | - | - | - | OK | OK | OK | OK | OK | OK |
C | - | - | - | OK | OK | OK | OK | OK | OK |
D | - | - | - | - | OK | OK | OK | OK | OK |
E | - | - | - | - | OK | OK | OK | OK | OK |
1 | - | - | - | - | OK | OK | OK | OK | OK |
2 | - | - | - | - | OK | OK | OK | OK | OK |
3 | - | - | - | - | OK | OK | OK | OK | OK |
Interstorey drifts in frames of
In the 2D model greater drifts were obtained, while in 3D model the drifts were reduced by the contribution of the diaphragms. Also it were determined the maximum torsional moments in each column before the implementation of R1_3 earthquake in all supports for the four combinations described in Table 8. In Figure 13 have been plotting torsional moments in function of time for the four combinations, where nodes appointed by n111 until the n513 are corresponding to supports, while Figure 14 shows the maximum torsional moment range for each column from three-dimensional analysis. It is evident that for the accelerograms used, the maximum torsional moments occurs in the extreme columns and in the columns located in the intersection of the structure. This is an important feature that confirms the negative effect of the irregularity combined with the seismic action. The torsional moments for the other seismic combinations used in this study were obtained using the same procedure.
|
||
1 | Corner column. n513 | 64225 |
2 | Corner column. n512 | 76000 |
3 | Corner column. n513 | 41000 |
4 | Corner column. n512 | 65000 |
6. Conclusions
In order to know the seismic response of the studied building it were used analytical methods considering the seismic hazard level and structural regularity criteria. The elastic analysis applied to the
From dynamic analysis there were computed the global and interstorey drifts for all three evaluated models determining the dynamic response of these structures and controlling the damage level reached in them. With the global drifts, it was evaluated the threshold of the collapse Limit State, which corresponds to the maximum value of 2.5%.
Three-dimensional dynamic analysis applied to
Interstorey drifts of
Inelastic static analysis is more reliable than linear methods in the prediction of the parameters of response of buildings, although this method has no response on the effects of higher modes of vibration. A more reliable and sophisticated method is the 2D no linear dynamic analysis, where it can be better determined the likely behavior of the building in response to the earthquake. However, the uncertainties associated with the definition of accelerograms used in these analysis and properties of coplanar structural models can be reduced with the implementation of the dynamic 3D analysis because there are considered factors associated with structural redundancy and are used more actual values in terms of rigidity of resistant structural lines.
The Quadrants Method presented in this paper is suitable for rapid and reliable evaluation of structures, with a low calculation effort. The cases studied demonstrated that the method can provide a reliable criterion to predict if any structure would have an inadequate seismic performance based on the results of static non-linear analysis. Results obtained from dynamic non-linear analysis confirmed the results obtained from the application of the Quadrants Method.
Despite the plan irregularity of the studied building, the Quadrants Method was suitably in order to predict that its lateral stiffness was not enough. Dynamics analysis confirmed this feature, then the cross sections of this building was resizing and details of the confinement were improved in order to meet the regulations of the current version of the Venezuelan seismic code. The seismic performance of the new designed building was tested with the Quadrants Method and dynamic analysis, showing that the resizing structure met all the Limit States used in this research.
Acknowledgments
Authors wish to acknowledge to the Scientific and humanistic Council of Lisandro Alvarado University for the financial support of the research in the field of non-regular structures. Also the authors would highlight the role of the International Center of Numerical Methods for Engineering in the collaborative research.References
- 1.
Grases J., Altez R. and Lugo M. Destructives Earthquakes Catalogue. Venezuela 1530/1998. Central University of Venezuela. Natural Sciiences, Physics and Mathematics Academy. Engineering School. Caracas, Venezuela, 1999. - 2.
Pérez O. and Mendoza J. Seismicity and tectonics in Venezuela and surroundings areas. Earth Physics,1998, 10, pp 87-110. - 3.
Barbat A., Mena U. and Yépez F. Probabilistic evaluation of seismic risk in urban zones. International magazine for numerical methods for Calculus and engineering projects. 1998, 14, 2, 247-268. - 4.
Calvi, G., Pinho, R., Magenes, G., Bommer, J., Restrepo, L and, Crowley, H. Development of Seismic Vulnerability Assessment Methodologies over the Past 30 Years, ISET Journal of Earthquake Technology, 2006, Paper No. 472, 43, 3, 75-104 - 5.
Fajfar P. Nonlinear analysis method for performance based seismic design. Earthquake Spectra, EERI, United States of America, 2000, 16, 3, 573-591. - 6.
Vielma, J. C., Barbat, A. H. and Oller, S. Framed structures earthquake resistant design. International Center for Numerical Methods in Engineering (CIMNE) Monograph. Earthquake Engineering Mongraphs. Barcelona, Spain, 2011. - 7.
Covenin 1756:01. Earthquake-resistant Design code. Part 1. Fondo Norma. Caracas, Venezuela, 2001. - 8.
Covenin 1753:06. Design and construction of buildings with structural concrete. Fondo Norma. Caracas, Venezuela, 2006. - 9.
Covenin 2002:88. Minimum design loads and criteria for buildings code. Fondo Norma. Caracas, Venezuela, 1988. - 10.
Vielma, J. C., Barbat, A. and Martínez, Y. The Quadrants Method: A procedure to evaluate the seismic performance of existing buildings. 15 World Conference on Earthquake Engineering. Lisbon, Portugal, 2012. - 11.
Herrera, R., Vielma, J. C., Ugel, R., and Barbat, A. Optimal design and earthquake resistant design evaluation of low rise framed rc building. Natural Science, Earthquake Special Issue. August 2012. California, USA. Doi:10.4236/ns. - 12.
Vielma, J. C., Barbat, A. H. and Oller, S. Seismic safety of RC framed buildings designed according modern codes. Journal of Civil Engineering and Architecture. 2011, 5:7, 567-575. David Publishing Company. Chicago, United States of America. - 13.
Priestley M., Calvi G. and Kowalski M. Displacement-based seismic design of structures. IUSS Press. Pavia. Italia, 2007. - 14.
Vielma, J. C., Barbat, A. H. and Oller, S. Nonlinear structural analysis. Application to evaluating the seismic safety. Nova Science Publishers. New York,2009. - 15.
Di Sarno, L. and Elnashai, A. Fundamentals of Earthquake Enginering. John Wiley and Sons. Chichester, United Kingdom, 2008. - 16.
Park, R. State-of-the-art report: ductility evaluation from laboratory and analytical testing. Proceedings 9th WCEE, IAEE, Tokyo-Kyoto, Japan, Vol VIII: 605-616, 1988. - 17.
Vielma, J. C., Barbat, A. H. and Oller, S. Seismic sizing of RC buildings according to energy-based amplification factors (In Spanish). Hormigón y acero. 2011. 63:263, 83-96. - 18.
Elnashai A., Papanikolau V. and Lee, D. ZEUS-NL, A system for Inelastic Analysis of Structures. User Manual. Mid-America Earthquake Center report no. MAE, Illinois University. Urban, Champagne, Illinois, 2011. - 19.
Papanikolaou, V. and Elnashai, A. Evaluation of conventional and adaptive pushover analysis I: methodology, Journal of Earthquake Engineering, 2005, 9, 6, 923-941 - 20.
Mwafy, A. and Elnashai, A. Static pushover versus dynamic collapse analysis of RC buildings, Engineering Structures, 2001, 23, 407–424 - 21.
Vamvatsikos D. and C. Allin Cornell (2002). Incremental dynamic analysis, Earthquake Engineering and Structural Dynamics, 2002, 31, 491–514. - 22.
Kappos, A. and Stefanidou, S. A deformation-based seismic design method for 3D R/C irregular buildings using inelastic dynamic analysis. Bulletin of Earthquake Engineering, Springer, Netherlands, 2010, 8, 4, 875-895 - 23.
UCLA-CIMNE. Compatible accelerograms with elastic design spectrums generation software. (PACED). International Center for Numerical Methods in Engineering. Universidad Centroccidental Lisandro Alvarado. Venezuela, 2009. - 24.
Fajfar, P., Marusic, and D., Perus, I. Torsional effects in the pushover-based seismic analysis of buildings. Journal of Earthquake Engineering, 2005, 9, 6, 831-854.