SciELO - Scientific Electronic Library Online

 
 issue96Potential infiltration and morphometry in the Arroyo Grande basin, Sucre ColombiaAttraction and repulsion forces in melt-textured and sintered YBCO-superconductors: a comparative study author indexsubject indexarticles search
Home Pagealphabetic serial listing  

Services on Demand

Journal

Article

Indicators

Related links

  • On index processCited by Google
  • Have no similar articlesSimilars in SciELO
  • On index processSimilars in Google

Share


Revista Facultad de Ingeniería Universidad de Antioquia

Print version ISSN 0120-6230On-line version ISSN 2422-2844

Rev.fac.ing.univ. Antioquia  no.96 Medellín July/Sept. 2020

https://doi.org/10.17533/udea.redin.20191045 

Artículo original

Effect of the viscous damping on the seismic response of Low-rise RC frame building

Efecto de la amortización viscosa sobre la respuesta sísmica en edificio de hormigón armado bajo

Fatima-Zohra Baba-Hamed1  * 

Luc Davenne2 

1 Department of Civil Engineering, University of Sciences and Technology Mohamed Boudiaf of Oran. BP 1505 E1 M’Naouer. C. P. 31130. Oran, Algeria.

2 Laboratoire Energetique Mecanique Electromagnetisme, University of Paris Ouest. 50 rue de sèvres. C. P. 92410. Ville d’Avray, France.


ABSTRACT

The equivalent viscous damping is a key parameter in the prediction of the maximum nonlinear response. Damping constitutes a major source of uncertainty in dynamic analysis. This paper studies the effect of using viscous damping, on the reduction of the seismic responses of reinforced concrete RC frame buildings modeled as three-dimensional multi degree of freedom (MDOF) systems, and the use of nonlinear time history analysis as a method of visualized behavior of buildings in the elastic and inelastic range. This study focuses on the implications of the available modeling options on analysis. This article illustrates the effect of using the initial or tangent stiffness in Rayleigh damping in analysis of structures. Correspondingly, this work is also concerned with the estimation of Rayleigh, mass-proportional or stiffness-proportional damping on engineering demand parameters (EDPs). As a result of a series of considerations, a damping modeling solution for nonlinear time history analysis (NLTHA) was carried out to compute the damage index. The application example is a building designed according to reinforced concrete code BAEL 91 and Algerian seismic code RPA 99/Version 2003 under seven earthquake excitations. The simulations demonstrated the accuracy and effectiveness of the proposed method to account for all of the above effects.

Keywords: RC frame building; damping approach; nonlinear time history analysis; engineering demand parameters

RESUMEN

El amortiguamiento viscoso equivalente es un parámetro clave en la predicción de la respuesta máxima no lineal. La amortiguación constituye una fuente importante de incertidumbre en el análisis dinámico. Este artículo estudia el efecto del uso de la amortiguación viscosa, en la reducción de las respuestas sísmicas de los edificios de estructuras de concreto reforzado con hormigón reforzado modeladas como sistemas tridimensionales de múltiples grados de libertad, y el uso del análisis de historia de tiempo no lineal como un método de comportamiento visualizado de Construcciones en la gama elástica e inelástica. Este estudio se centra en las implicaciones de las opciones de modelado disponibles en el análisis. Además, este artículo ilustra el efecto del uso de la rigidez inicial o tangente en la amortiguación de Rayleigh en el análisis de estructuras. Este trabajo también se ocupa de la estimación de Rayleigh, la amortiguación proporcional a la masa o la rigidez proporcional en los parámetros de demanda de ingeniería. Como resultado de una serie de consideraciones, se llevó a cabo una solución de modelado de amortiguamiento para el análisis de historial de tiempo no lineal para calcular el índice de daño. El ejemplo de aplicación es un edificio diseñado de acuerdo con el código de hormigón armado BAEL91 y el código sísmico argelino RPA99/2003 bajo siete excitaciones de terremoto. Las simulaciones demostraron la precisión y la eficacia del método propuesto.

Palabras clave: Construcción de cuadros RC; enfoque de amortiguación; análisis de historia de tiempo no lineal; parámetros de demanda de ingeniería

1. Introduction

The local failures generated by earthquakes, in turn induce vibrations of the soil. These waves cause the vibration of buildings located in the environment of the epicenter of the earthquake [1]. Several methods have been proposed to evaluate the seismic structural performance in the development of performance-based seismic engineering. The fundamental question is to determine the seismic demand and the capacity of collapse proportional to the Grounds accelaration caused by earthquakes.

Different approaches for assessing structural collapse capacity with the aim to preserve life safety differ from the simplest approach, based on a simple single-degree-of-freedom (SDOF) response model to complex nonlinear dynamic analyses done in a structural model, which is analysed for ground motion records [2].

The capacity spectrum method is highly used. The seismic action is defined by means of the elastic response spectra, the fragility of the building by means of the capacity curve. The latter is calculated from an incremental nonlinear static analysis, commonly known as "Pushover Analysis" [3-7].

With the growth in computational power, new complex analysis tools have been developed to calculate nonlinear response of structures to ground motions. The procedure is accessible for preparation of analysis models, which can capture nonlinear behavior of most components.

Nonlinear time history analysis (NLTHA) tools are becoming regularly available to engineers. This method evaluates structural seismic performance by applying a series of ground motion acceleration induced to the structure. The inelastic earthquake response of a concentrated plasticity model of a building can be greatly affected by how damping is defined [8]. The modeling of damping is a difficult task in dynamic nonlinear analysis [9] caused by the diversity of mechanisms that may be present at once in the damping of the structure [10]. In the inelastic time history analyses of structures in seismic motion, part of the seismic energy that is imparted to the structure is absorbed by the inelastic structural model, and Rayleigh damping is commonly used in practice as an additional energy dissipation source [11].

Energy dissipation and ductility are part of diverse parameters that influence viscous damping ratio [12,13]. The viscous damping matrix is fundamental in nonlinear time history analysis to calculate the inelastic deformation. The effect of viscous damping on the response of the structure is greater than the hysteresis damping connected to the inelastic deformation in the case of moderate earthquakes. Therefore, a correct evaluation of viscous damping and its effects in structural dynamics and earthquake engineering is essential [14].

Early research on damping modeling problems in inelastic structures dates back to the work of [15-17].

The variation of the non-linear response indicators through the average ductility, hysteretic-to-input energy ratio and the average number of yield incursions, as a function of the additional damping model used in the case of inelastic frame structures in seismic loading have been studied by [18].

A capped viscous damping formulation to overcome some of the problems pointed out was proposed by [19] where the use of Rayleigh damping can lead to damping forces that are unrealistically large compared with the restoring force on practical situations.

Charney [17] investigated the effects of global stiffness changes on the seismic response of a five-story structure when Rayleigh damping is used. [17] also investigated the effects of local stiffness changes at the seismic structural time-history response when initial stiffness-based Rayleigh damping is used [11]. Each structural element with an equivalent combination of an elastic element with initial stiffness-proportional damping and yielding springs at the two ends without stiffness-proportional damping was investigated and modeled by [20].

Erduran [21] used tangent stiffness-based Rayleigh damping models to analyze their effect in the story drift ratios, floor accelerations, and damping forces in three-story and nine-story steel moment-resisting frame buildings for three seismic hazard levels [11].

The aims of this article are: (1) to assess the effects of the Rayleigh damping models on the results of nonlinear time history analysis and (2) to demonstrate the impact of the damping model on the assessed engineering demand parameters and conclude by the proposal of a modeling solution appropriate for RC buildings to the RPA 2003 code.

To that aim, the article is organized in five sections; Section 1 introduces the problematic studied in this paper. Section 2 focuses on the presentation of damping in nonlinear time history analysis. Section 3 is devoted to the presentation of the finite element model used in this study. In this section, the viscous damping will also be applied to the case study using different types of approaches to model the damping matrix C.

Section 4 aims to present the main results of the sensitivity analysis. Finally, we end the article with a conclusion based on the observations in the previous section.

2. Damping in time-history analyses

In earthquake engineering, damping is usually present as a ratio or fraction of critical damping, called the damping ratio ξ, which is a property of the system material and independent of its mass and stiffness [12]. The equivalent viscous damping is the sum of the elastic and the hysteric damping [22]. Damping values depend on the construction materials, vibration amplitude, fundamental period and mode shapes, type of connections and the building configuration [23].

Many studies [24-32] have been interesting to draw the exact definition and estimation of equivalent viscous damping.

Experimental investigations have shown that the vibrations of real structures usually lie between the viscous and frictional responses. However, the viscous assumption is convenient to use analytically and is sufficiently accurate for most purposes [33]. The damping of the structure is assumed to be viscous and frequency dependant for the sake of convenience in the analysis. The most popular method is to solve the equation of motion using the modal analysis; in this case, damping values are directly assigned to the modes [34]. Damping ratios can be calculated using the Caughey series [35] which Rayleigh damping [36] is a special case.

Rayleigh damping can be used in seismic simulations either to account for energy dissipation mechanisms that are external to the structure or for energy absorption mechanisms that are internal to the structure [37].

In Rayleigh damping, the damping matrix is a linear combination of mass-proportional and stiffness proportional damping terms as shown in Equation (1):

Where ( and ( are real scalars with 1/sec and sec units respectively. Modes of classical damped systems preserve the simplicity of the real normal modes.

These coefficients have been evaluated from two specified damping ratios ((i and (j) and natural frequencies (ωi and ωj) corresponding to two “references” modes (i and j) of vibration, respectively as shown in Equation (2).

Another important aspect is which characteristics of the structure, initial or inelastic, are used in the computation of the damping matrix. The damping matrix can be constructed by using the initial stiffness matrix Ko, which is constant throughout the analysis, or the tangent stiffness Kt, which is updated at each time-step following the stiffness changes. The values of the constants α and β are based on initial characteristics, which constitutes a computational convenience on the ground that the evolution of the frequency spectrum is not known. The choice of initial damping model between a constant damping matrix and tangent-stiffness proportional damping matrix could be significant, particularly for short-period structures [38]. With initial stiffness elastic damping, the damping coefficient is constant throughout the analysis, even in the inelastic range of response, and is based on the initial elastic stiffness. With tangent-stiffness damping, the damping coefficient is proportional to the instantaneous value of the stiffness and it is updated whenever the stiffness changes [39].

3. Finite element modeling

3.1 Case study

The case study is a 3-story RC frame building designed according to reinforced concrete code BAEL 91 and Algerian seismic code RPA 99/Version 2003 with the following parameters: zone of moderate seismicity, zone Z IIa. A seismic behavior factor of R = 5 was taken into account for reinforced concrete frames without masonry infill.

The Algerian paraseismic code RPA99/version2003 [40] was used for the pre-dimensioning of cross sections of columns and beams. Then, it was used to consider the seismic action as an accidental action in the sense of the limit state calculation philosophy.

Finally, this code was used to determine in beams and columns:

  • The minimum and maximum percentage of longitudinal reinforcement

  • The minimum length of recovery

  • The minimum transverse reinforcement

The BAEL91code was used to calculate longitudinal and transverse reinforcement in beams and columns.

The building has a rectangular plan configuration as shown in Figure 1. The dimensions of the building are 16.6 m long and 8.8 m wide with a total height of 9.18 m in elevation and all the floors are the same height of 3.06 m.

Figure 1 Plan view of the structure 

Concrete characteristic compressive strength equal to fck = 25 N/mm2 and steel characteristic yielding strength equal to fyk = 400 N/mm2 are adopted. The member cross-section sizes and steel bars are given in Table 1.

Table 1 Dimensions and reinforcement details of the RC beams and columns 

3.2 Numerical Modeling of structural element

To perform time history analysis, frames are modeled in SeismoStruct-V-6 [41] finite element software that is able to perform nonlinear dynamic analysis.

A uniaxial model that follows the constitutive relationship proposed by [42] and the cyclic rules proposed by [43] representing the concrete. The confinement effects provided by the lateral transverse reinforcement are incorporated through the rules proposed by [42] whereby constant confining pressure is assumed throughout the entire stress-strain range. The model adopted for the steel reinforcement is based on the Menegotto-Pinto [44] model coupled with the isotropic hardening rules proposed by [45].

The modeled structure is shown in Figure 2. The beams were divided longitudinally into 4 elements and each beam and column element was divided transversely into 200 by 200 fibers elements.

Figure 2 Finite element model of the 3 stories building 

The capacity curve of the buildings obtained from nonlinear static (pushover) analysis under an invariant first mode load pattern is represented in Figure 3.

Figure 3 The capacity curve of the typical model buildings derived from pushover analysis 

3.3 Ground motion selection earthquake

Ground motions observed at a definite site usually depend on the type of the earthquake (fault mechanism, the way seismic wave propagation) [46]

A suite of ground motion records needs to be selected to perform nonlinear time history analysis [Figure 4]. Selected ground motions are shown in Table 2.

Table 2 Earthquakes used for non-linear dynamic analyses 

Two sets of ground motions that are representative of different intensity measures are used in the response history analysis. The first set represents a moderate seismic zone (zone IIa) and the second a high seismic zone (zone III). The ground motion set selected is scaled to represent ground motions with a RPA99/V.2003 elastics designs spectrum's zone IIa and zone III.

Figure 4 Acceleration response spectrums for the selected earthquakes together with the RPA99/version2003 elastics designs spectrums 

The SeismoMatch software [47] was used in the ground motion scaling. It has the ability to scale ground motion acceleration to match to a given target response spectrum. That consists of three steps. First, load all ground motions to the SeismoMatch software as a single ground motion or multiple ground motions. Second, define a target spectrum in the software or load user defined target spectrum of the software. Third, enter interval that needs to match such as 0.2T1 to 1.5T1 and enter a scale factor if available.

Two sets of ground motions satisfy the Eurocode 8 [48] provisions:

The mean of zero period spectral response acceleration values should not be smaller than the value of S.ag (S is the soil factor, ag is the design ground acceleration); in the range of periods between 0.2T1 and 2T1, where T1 is the fundamental period of the structure, in the direction where the accelerogram is applied, no value of the mean 5% damping elastic response spectrum, calculated from all time-histories should be less than 90% of the corresponding value of the 5% damping elastic design spectrum.

3.4 Modeling of viscous damping

The application of viscous damping was only to the building superstructure with the use of many different of computing the damping matrix C that covers a wide range of options available in many existing FE programs [Table 3]. The Table 3 also shows the relevant frequencies obtained from modal analyses.

Table 3 Several modeling approaches for viscous damping 

Approach 1 used the fixed-base frequency of the superstructure to calculate the damping coefficient β. This coefficient remained constant throughout the analysis. Consequently, the damping matrix C remained also constant. In the approach 2, the modification of stiffness allows updating the damping matrix C at each step of the analysis. Approach 3 used the fixed-base frequency of the superstructure in order to compute the damping coefficient α. This coefficient remained constant throughout the analysis. Therefore, the damping matrix C also remained constant.

On approach 4, Rayleigh damping was used where the initial stiffness matrix was used to calculate the damping matrix C by Equation 1.

First and sixth mode frequencies are used to calculate the damping coefficients. This is explained by the fact that the first six modes were observed to contribute to the response significantly after a modal analysis. In approach 5, the same approach is used with a tangent stiffness matrix to calculate the damping matrix.

At last, reduced frequencies as used in Models 6 and 8, in according with the recommendations of [17], who allows us to suppose that the definition of Rayleigh damping based on reduced modal frequencies to account for yielding in these modes. Results are then more reasonable concerning the response evaluation for buildings that are estimated to deform outside their elastic limit.

The 5% Rayleigh damping is used; first and sixth frequency modes are reduced (wi = 0.707w1 and wj = 0.707w6).

Figure 5 shows the damping ratio versus frequency curves for different approaches with the damping coefficients α and β computed using ζ=5%.

Figure 5 Variation of damping ratio with frequency computed for ζ=5% 

4. Discussion and result

The selected Engineering demand parameters include maximum floor displacements, maximum story shear forces and maximum inter-story drift ratios. The ratios of damping forces to mass forces are identified to recognize the damping models that may have as potential results, irrationally high damping forces.

4.1 Damping forces

Figure 6 illustrates the mean values of peak damping force/mass force ratios. This result was obtained by the use of different damping approaches for two seismic intensity measures.

The Ratio of damping force to mass force was observed to identify the damping models which may have as potential results, unreasonably high damping forces [21].

Figure 6 Damping force to base shear force ratio for several damping models for two seismic intensity measures 

Figure 6 shows that the damping force/mass force ratio growth with an increase in the seismic intensity for all the damping models considered for the building. Moreover, for two seismic intensity measures, mass proportional damping leads to the highest damping force/base shear ratios. While, stiffness proportional damping leads to the lowest ratios. The growth in the damping force/mass force ratio with an increase in the seismic intensity is more marked for the mass-proportional damping model than the stiffness-proportional model.

Rayleigh damping models with two approaches result in average damping force/base shear ratios between mass-proportional and stiffness-proportional damping models. When using reduced frequencies to calculate the coefficients α and β as suggested by [17], damping force/base shear ratios are inferior compared with their equivalents where natural frequencies are used.

The damping force/mass force is less pronounced when the initial stiffness matrix is used. This reduction is caused by the invariability of the damping throughout the analysis.

4.2 Floor response spectra

Figure 7 shows the mean floor acceleration response spectra at the roof level of the building using different damping approaches for two seismic intensity measures.

Figure 7 Floor response spectra at the roof level of building calculated from response history analysis for (a) 0.188g, (b) 0.313g 

Before the fundamental period of the building which is equal to 0.52 s, we find that the spectral acceleration at the roof level for the mass proportional damping is the most pronounced for the two seismic intensities. As well as the spectral floor acceleration at the roof level for the initial stiffness proportional damping is the least pronounced for the two seismic intensities. In addition, the floor spectral accelerations for tangential Rayleigh damping are the closest to the spectral accelerations of the RPA spectra corresponding to the fundamental period of the building.

4.3 Floor displacements

Figure 8 presents the mean values of maximum floor displacements obtained using different damping approaches and pushover analysis for two seismic intensity measures.

Figure 8 Maximum floor displacements calculated from response history analysis for (a) 0.188g, (b) 0.313g 

For the RC building studied, very similar maximum floor displacements are obtained for all damping approaches.

The seismic performance is strongly marked by the first mode with very modest contribution from the higher modes. We obtained globally the same damping matrix for different approaches.

In pushover analysis, the maximum displacements in the second and third levels are superior relative to the nonlinear time history analysis. Particularly in seismic intensity measurements of 0.313 g. In nonlinear time history analysis, seismic ground motion is reduced by the use of a damping matrix.

Different approaches yielded structural responses whose amplitude differences in maximum floor displacements were quantified by calculating ratios for two measurements of seismic intensities. [Table 4].

Maximum floor displacements ratios are calculated as:

% = (Maximum floor displacements - Mean maximum floor displacements) / Mean maximum floor displacements

Table 4 Maximum floor displacement ratios calculated from response history analysis for (a) 0.188g, (b) 0.313g 

4.4 Story shear forces

Figure 9 presents the mean values of maximum story shear forces obtained using different damping approaches and pushover analysis for two seismic intensity measures. The same remarks are observed concerning floor displacements.

Figure 9 Maximum story shear forces calculated from response history analysis for (a) 0.188g, (b) 0.313g 

The differences in magnitude of maximum story shear forces between structural responses obtained with different approaches have been quantified by computing ratios for two seismic intensity measures [Table 5].

Maximum story shear force ratios are calculated as:

% = (Maximum story shear forces - Mean maximum story shear forces) / Mean maximum story shear forces

Table 5 Maximum story shear force ratios calculated from response history analysis for (a) 0.188g, (b) 0.313g 

4.5 Maximum inter-story drift ratio

The inter-story drift performance of a multistory building is an important index of structural and non-structural damage of the building under a range of levels of earthquake motion. The inter-story drift performance has become a principal design consideration in performance design. The system performance levels of a multistory building are evaluated on the basis of the inter-story drift values along the height of the building under different levels of earthquake motion.

The Inter-Story Drift Ratio is evaluated by the Equation (3):

Where

  • IDR is the inter-story drift

  • di is the displacement of i story

  • di-1 is the displacement of i-1 story

  • hi is the height of i story

Figure 10 presents the mean values of the Inter-story drift ratios obtained using different damping approaches and pushover analysis for two seismic intensity measures.

Figure 10 Inter-story drift ratios obtained from response history analysis for (a) 0.188g, (b) 0.313g 

The Inter-story drift ratios are more important in the second floors in all approaches; because of the importance of the difference of displacement between the first and the second floor.

Initial and tangent Rayleigh damping models with two approaches result in average Inter-Story Drift Ratio, ratios among mass-proportional and initial and tangent stiffness-proportional damping models.

To reduce Rayleigh damping approach the Inter-Story Drift Ratio is more important compared to Rayleigh damping approach. The reduction of frequencies generates the amplification on the periods, causing attenuation of the damping matrix.

The differences in magnitude of Inter-story drift between structural responses obtained with different approaches have been quantified by computing ratios for two seismic intensity measures [Table 6].

Maximum Inter-story drift ratios are calculated as:

% = (Inter-story drift - Mean inter-story drift) / Mean inter-story drift

Table 6 Inter-story drift ratios calculated, from response history analysis for (a) 0.188g, (b) 0.313g 

5. Conclusions

The investigation of the imports of modeling viscous damping nonlinear time history analysis of RC buildings by different approaches is evaluated using the response of a three-story building.

Lowest damping force/mass force ratios are given by the initial stiffness matrix. Coefficients α and β are constant and independent, with the time. The mass forces and the damping forces were limited because they are proportional to the initial stiffness matrix. Furthermore, mass-proportional damping implies highest damping force/base shear force ratios for different seismic intensity measures for RC frame buildings.The damping force/base shear force ratios resulting from the use of mass-proportional damping model intensificates with a growth in the seismic intensity unlike damping force/base shear force ratios obtained using stiffness-proportional damping which remains constant.

Initial damping reduces the engineering demand parameters. It is caused by the invariability of the damping throughout the analysis.

The Rayleigh damping approach gives smaller maximum engineering demand parameters ratios. [Table 4, Table 5 and Table 6]. It places the different results in the analysis of performance between mass-proportional and stiffness-proportional damping models.

The reduced Rayleigh damping leads to more important results in the analysis of performance compared to Rayleigh damping approach. The reduction of frequencies generates the amplification on the periods, causing attenuation of the damping matrix.

The interpretation of results recapitulated in this article recommends that neither mass-proportional nor stiffness proportional damping are appropriate for inelastic nonlinear time history analysis of RC buildings.

The damping model used in the nonlinear dynamic do not have a significant influence the Inter-Story Drift ratio in the RC buildings according to RPA2003 [40] which represents the low-rise buildings. With the models examined, tangent Rayleigh damping appears to be the ideal alternative option for the RC frame buildings.

7. Acknowledgements

We thank the Laboratory Energétique Mécanique Electromagnétisme (LEME) for their welcome.

References

[1] J. Mazars, S. Grange, and C. Desprez, “Seismic risk: Structural response of constructions,” European Journal of Environmental and Civil Engineering, vol. 15, no. sup1, 2011. [Online]. Available: https://doi.org/10.1080/19648189.2011.9695309Links ]

[2] R. Villaverde, “Methods to assess the seismic collapse capacity of building structures: State of the art,” Journal of Structural Engineering, vol. 133, no. 1, january 2007. [Online]. Available: https://doi.org/10.1061/(ASCE)0733-9445(2007)133:1(57)Links ]

[3] T. Rossetto and A. Elnashai, “Derivation of vulnerability functions for european-type rc structures based on observational data,” Engineering Structures, vol. 25, no. 10, august 2003. [Online]. Available: https://doi.org/10.1016/S0141-0296(03)00060-9Links ]

[4] A. Chopra and R. Goe, “Direct displacement-based design: Use of inelastic vs. elastic design spectra,” Earthquake Spectra, vol. 17, no. 1, february 2001. [Online]. Available: https://doi.org/10.1193/1.1586166Links ]

[5] A. Chopra andR. Goe , “A modal pushover analysis procedure for estimating seismic demands for buildings,” Earthquake Engng Struct. Dyn., vol. 31, no. 3, march 2002. [Online]. Available: https://doi.org/10.1002/eqe.144Links ]

[6] P. Fajfar, “A nonlinear analysis method for performance-based seismic design,” Earthquake Spectra , vol. 16, no. 3, august 2000. [Online]. Available: https://doi.org/10.1193/1.1586128Links ]

[7] FEMA-273. (1997, Oct.) NEHRP guidelines for the seismic rehabilitation of buildings. Federal Emergency Management Agency. Washigton D.C., EE. UU. [Online]. Available: https://www.scinc.co.jp/nanken/pdf/fema273.pdfLinks ]

[8] A. Chopra and F. McKenna, “Modeling viscous damping in nonlinear response history analysis of buildings for earthquake excitation,” Earthquake Engng Struct. Dyn. , vol. 45, no. 2, february 2016. [Online]. Available: https://doi.org/10.1002/eqe.2622Links ]

[9] D. Pant, A. Wijeyewickrema, and M. ElGawady, “Appropriate viscous damping for nonlinear time-history analysis of base-isolated reinforced concrete buildings,” Earthquake Engineering & Structural Dynamics, vol. 42, no. 4, december 2013. [Online]. Available: https://doi.org/10.1002/eqe.2328Links ]

[10] S. Chang, “Nonlinear performance of classical damping,” Earthquake Engineering and Engineering Vibration, vol. 12, no. 2, june 2013. [Online]. Available: https://doi.org/10.1007/s11803013-0171-3Links ]

[11] P. Jehel, P. Léger, and A. Ibrahimbegovic, “Initial versus tangent stiffness-based rayleigh damping in inelastic time history seismic analyses,” Earthquake Engineering and Engineering Vibration , vol. 43, march 2014. [Online]. Available: https://doi.org/10.1002/eqe.2357Links ]

[12] A. Chopra , Dynamics of structures: theory and applications to earthquake engineering. New Jersey, EE.UU.: Prentice Hall, 1995. [ Links ]

[13] T. Heitz, C. Giry, B. Richard, and F. Ragueneau, “How are the equivalent damping ratios modified by nonlinear engineering demand parameters?” presented at 6th International Conference on Computational Methods in Structural Dynamics and Earthquake Engineering (COMPDYN), Rhodes Island, GR, 2017. [ Links ]

[14] D. Pan, G. Chen, and Z. Wang, “Suboptimal rayleigh damping coefficients in seismic analysis of viscously-damped structures,” Earthquake Engineering and Engineering Vibration , vol. 13, no. 4, december 2014. [Online]. Available: https://doi.org/10.1007/s11803-014-0270-9Links ]

[15] D. Chrisp, “Damping models for inelastic structures,” M.S. thesis, University of Canterbury, Christchurch, NZ, 1980. [ Links ]

[16] P. Shing and S. Mahin, “Elimination of spurious higher-mode response in pseudodynamic tests,” Earthquake Engineering Structural Dynamics, vol. 15, no. 4, may 1987. [Online]. Available: https://doi.org/10.1002/eqe.4290150403Links ]

[17] F. Charney, “Consequences of using rayleigh damping in inelastic response history analysis,” presented at Congreso Chileno de Sismología e Ingeniería Antisísmica IX Jornadas, Concepción, CL, 2005. [ Links ]

[18] P. Léger and S. Dussault, “Seismic-energy dissipation in MDOF structures,” Journal of Structural Engineering -asce, vol. 118, no. 5, may 1 1992. [Online]. Available: https://doi.org/10.1061/(ASCE)0733-9445(1992)118:5(1251)Links ]

[19] J. Hall, “Problems encountered from the use (or misuse) of Rayleigh damping,” Earthquake Engineering & Structural Dynamics, vol. 35, no. 5, april 2006. [Online]. Available: https://doi.org/10.1002/eqe.541Links ]

[20] F. Zareian and R. Medina, “A practical method for proper modeling of structural damping in inelastic plane structural systems,” Computers & Structures, vol. 88, no. 1-2, january 2010. [Online]. Available: https://doi.org/10.1016/j.compstruc.2009.08.001Links ]

[21] E. Erduran, “Evaluation of rayleigh damping and its influence on engineering demand parameter estimates,” Earthquake Engineering & Structural Dynamics , vol. 41, no. 14, November 2012. [Online].Available: https://doi.org/10.1002/eqe.2164Links ]

[22] S. Salawdeh and J. Goggins, “Direct displacement based seismic design for single storey steel concentrically braced frames,” Earthquakes and Structures, vol. 10, no. 5, february 2016. [Online]. Available: https://doi.org/10.12989/eas.2016.10.5.1125Links ]

[23] L. Sarno, A. Elnashai , and G. Manfredi, “Assessment of RC columns subjected to horizontal and vertical ground motions recorded during the 2009 L’Aquila (Italy) earthquake,” Engineering Structures , vol. 33, no. 5, may 2011. [Online]. Available: https://doi.org/10.1016/j.engstruct.2011.01.023Links ]

[24] P. Jennings, “Equivalent viscous damping for yielding structures,” Journal of the Engineering Mechanics Division, vol. 94, no. 1, pp. 103-116, 1968. [ Links ]

[25] P. Gulkan and M. Sozen, “Inelastic responses of reinforced concrete structure to earthquake motions,” Journal of the American Concrete Institute, vol. 71, no. 12, pp. 604-610, Dec. 1974. [ Links ]

[26] W. Iwan and N. Gates, “The effective period and damping of a class of hysteretic structures,” Earthquake Engineering & Structural Dynamics , vol. 7, no. 3, may/june 1979. [Online]. Available: https://doi.org/10.1002/eqe.4290070302Links ]

[27] E. Miranda and J. Ruiz, “Evaluation of approximate methods to estimate maximum inelastic displacement demands,” Earthquake Engineering & Structural Dynamics , vol. 31, no. 3, march 2002. [Online]. Available: https://doi.org/10.1002/eqe.143Links ]

[28] W. Kwan and S. Billington, “Influence of hysteretic behavior on equivalent period and damping of structural systems,” Journal of Structural Engineering , vol. 129, no. 5, may 2003. [Online]. Available: https://doi.org/10.1061/(ASCE)0733-9445(2003)129:5(576)Links ]

[29] M. Priestley and D. Grant, “Viscous damping in seismic design and analysis,” Journal of Earthquake Engineering, vol. 9, no. sup2, january 2005. [Online]. Available: https://doi.org/10.1142/S1363246905002365Links ]

[30] H. Dwairi, M. Kowalsky, and J. Nau, “Equivalent damping in support of direct displacement-based design,” Journal of Earthquake Engineering , vol. 11, no. 4, 2007. [Online]. Available: https://doi.org/10.1080/13632460601033884Links ]

[31] T. Heitz , C. Giry , B. Richard , andF. Ragueneau , “Identification of an equivalent viscous damping function depending on engineering demand parameters,” Engineering Structures , vol. 188, no. 4, june 1 2019. [Online]. Available: https://doi.org/10.1016/j.engstruct.2019.03.058Links ]

[32] M. Abbasi and M. Moustafa, “Effect of damping modeling and characteristics on seismic vulnerability assessment of multi-frame bridges,” Journal of Earthquake Engineering , april 2019. [Online]. Available: https://doi.org/10.1080/13632469.2019.1592791Links ]

[33] J. Smith, Vibration of Structures: Applications in civil engineering design. New York, EE.UU.: Broadview Press, 1988. [ Links ]

[34] A. Alipour andF. Zareian , “Study rayleigh damping in structures; uncertainties and treatments,” presented at 14th World Conference on Earthquake Engineering, Beijing, CN, 2008. [ Links ]

[35] T. Caughey and M. O’kelly, “Classical normal modes in damped linear dynamic systems,” Journal of Applied Mechanics, vol. 32, no. 3, 1965. [Online]. Available: https://doi.org/10.1115/1.3627262Links ]

[36] L. Rayleigh and J. Strutt, The Theory Of Sound. New York, EE.UU.:Dover Publications, 1896. [ Links ]

[37] P. Jehel , “A critical look into rayleigh damping forces for seismic performance assessment of inelastic structures,” Engineering Structures , vol. 78, november 1 2014. [Online]. Available: https://doi.org/10.1016/j.engstruct.2014.08.003Links ]

[38] E. Smyrou, M. Priestley , and A. Carr, “Modelling of elastic damping in nonlinear time-history analyses of cantilever RC walls,” Bulletin of Earthquake Engineering, vol. 9, no. 5, october 2012. [Online]. Available: https://doi.org/10.1007/s10518-011-9286-yLinks ]

[39] M. Priestley , D. Grant , and C. Blandon, “Direct displacement-based seismic design,” presented at NZSEE Conference, New Zealand, 2005. [ Links ]

[40] Algeria. Regles Parasismiques Algeriennes RPA99/Version 2003, Algerian Earthquake Resistant Regulations «RP A 99»/Version 2003, Ministry of Housing and Urbanism, 2003. [ Links ]

[41] (2015) Seismosoft. SeismoStruct v7.0.6. Seismosoft-Earthquake Engineering Software Solutions. Accessed Feb. 2015. [Online]. Available: Available: https://seismosoft.com/Links ]

[42] J. Mander, M. Priestley , and R. Park, “Theoretical stress-strain model for confined concrete,” Journal of Structural Engineering , vol. 114, no. 8, september 1988. [Online]. Available: https://doi.org/10.1061/(ASCE)0733-9445(1988)114:8(1804)Links ]

[43] J. Martínez andA. Elnashai , “Confined concrete model under cyclic load,” Materials and Structures, vol. 30, no. 3, pp. 139-147, Apr. 1997. [ Links ]

[44] M. Menegotto and P. Pinto, “Method of analysis for cyclically loaded rc frames including changes in geometry and non-elastic behaviour of elements under combined normal force and bending,” presented at IABSE reports of the working commissions, Zurich, SUI, 1973. [ Links ]

[45] F. Filippou, V. Bertero, and E. Popov, “Effects of bond deterioration on hysteretic behavior of reinforced concrete joints,” Earthquake Engineering Research Center, University of California, Los Angeles, CA, Tech. Rep. UCB/EERC-83/19, Aug. 1983. [ Links ]

[46] S. Mitrović, J. Ožbolt, and V. Travaš, “Three-dimensional finite element formulation for nonlinear dynamic analysis of seismic site and structure response,” Journal European Journal of Environmental and Civil Engineering, vol. 19, no. 7, 2015. [Online]. Available: https://doi.org/10.1080/19648189.2014.973534Links ]

[47] (2016) Seismomatch- a program for spectral matching of earthquake records. Seismosoft-Earthquake Engineering Software Solutions. Accessed Mar. 2016. [Online]. Available: Available: https://bit.ly/2knOV07Links ]

[48] Eurocode 8: Design of structures for earthquake resistance-Part 1 : General rules, seismic actions and rules for buildings, European Standard Norme Europeenne Europaische Norm, en 1998-1, 2004. [ Links ]

Received: July 17, 2019; Accepted: October 21, 2019

* Corresponding author: Fátima Zohra Baba Hamed e-mail: fatimazohra.babahamed@univ-usto.dz

6. Declaration of competing interest

None declared under financial, profesional and personal competing interests.

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License