Estimation of Cyclic Demand in Metallic Yielding Dampers Installed on Frame Structures

Featured Application: Assessment and design of hysteric dampers considering the impulsive e ﬀ ects of earthquakes. Abstract: The e ﬃ cacy of hysteretic dampers can be formulated as the number or equivalent inelastic cycles, or the ratio of normalised dissipated energy to displacement ductility. This parameter is used in the design of framed structures with supplemental dampers and it is strongly inﬂuenced by the impulsive e ﬀ ects of earthquakes and other structural parameters. This paper presents an estimate of the cyclic demand of dampers installed in reinforced concrete frames, based on nonlinear time history analyses. Statistical analyses of the results are used to highlight relevant parameters and calibrate a predictive formula and upper-bound design values. The collapse pattern of the frames seems to have no e ﬀ ect on the e ﬃ cacy of the dampers; thus, the seismological parameters that describe impulsivity should drive the design of hysteric dampers.


Introduction
Metallic yielding dampers are a cost-effective technology for passive seismic control of the vibration of structures. In buildings, these types of dampers are generally installed into new or retrofitted framed structures to ensure reliable seismic performance through inelastic energy dissipation. Unlike velocity dependent dampers such as viscous and viscoelastic or other design strategies [1,2], the force output and amount of dissipated energy does not depend on the loading rate [3], but rather on the magnitude of displacement and on accumulated displacement. Maholtra [4] observed that pulse-like ground motions have dramatic effects on structural response, such as reduced apparent flexibility of high-rise and base-isolated buildings, increased base shear and inter-storey drifts in high-rise buildings, reduced effectiveness of supplemental damping, and increased ductility demand.
Consequently, near-fault ground motions can expose aseismic devices to critical working conditions and a significant reduction in their efficacy. This issue has been studied in base isolated systems [5][6][7][8][9][10], raising concerns about the viability of base isolation and leading to specific design provisions in the near-fault to accommodate lateral deformation and optimise supplemental damping. The effects of near-fault seismicity on the seismic performance of aseismic dampers has also been addressed. Although dampers still provide seismic protection compared to bare frames in the near-fault [11], it has been shown that earthquake impulsive effects lead to the need for higher yield strength and to a decrease in the available local ductility in the dampers [12]. This results in lower energy dissipation capacity and requires specific design considerations [13][14][15][16][17]. The number of equivalent cycles of the hysteretic dampers is also relevant to fatigue issues related to heating due to the rate of energy dissipation.
Near-fault records are characterised by shorter durations and tend to produce less inelastic cycles than far-field records. In the case of hysteretic dampers, energy dissipation occurs due to metallic yielding and it is strongly related to cumulative inelastic deformation or cyclic demand in the damper. One approach to tackle cyclic demand is by considering the relationship between dissipated energy at a maximum or target displacement.
Zahran and Hall [18] investigated the response of elastic-perfectly plastic single-degree-of-freedom (SDOF) systems and noted that structures can absorb more or less energy despite experiencing the same peak displacement. Because energy dissipation must be related to cyclic or cumulative demand, they formulated the concept of equivalent number of yield cycles N. This is defined in Equation (1) and shown in Figure 1, as the ratio of the total energy, E H that is dissipated by yielding within the structure when subjected to the ground motion of an earthquake to the area under the monotonic load-displacement curve at the same maximum displacement, δ m , experienced during the vibration.
lower energy dissipation capacity and requires specific design considerations [13][14][15][16][17]. The number of equivalent cycles of the hysteretic dampers is also relevant to fatigue issues related to heating due to the rate of energy dissipation. Near-fault records are characterised by shorter durations and tend to produce less inelastic cycles than far-field records. In the case of hysteretic dampers, energy dissipation occurs due to metallic yielding and it is strongly related to cumulative inelastic deformation or cyclic demand in the damper. One approach to tackle cyclic demand is by considering the relationship between dissipated energy at a maximum or target displacement.
Zahran and Hall [18] investigated the response of elastic-perfectly plastic single-degree-offreedom (SDOF) systems and noted that structures can absorb more or less energy despite experiencing the same peak displacement. Because energy dissipation must be related to cyclic or cumulative demand, they formulated the concept of equivalent number of yield cycles N. This is defined in Equation (1) and shown in Figure 1, as the ratio of the total energy, EH that is dissipated by yielding within the structure when subjected to the ground motion of an earthquake to the area under the monotonic load-displacement curve at the same maximum displacement, δm, experienced during the vibration.
Zahran and Hall's results show that the equivalent number of yield cycles N is a useful parameter to define the severity or damage potential of the ground motion. It is influenced by (i) the characteristics of the ground motion: N increases with duration and tends to decrease for near-field records; (ii) the fundamental period: N increases for structures in the 0.2-2 Hz range; and (iii) ductility: N greatly increases as the displacement ductility increases. The number of equivalent cycles N is sometimes termed as neq in the literature. Akiyama [19][20][21] addressed the correspondence between the cumulative inelastic deformation ratio, η defined in Equation (2) (i.e., dissipated hysteretic energy normalised to the yield force Qy and displacement δy) and the inelastic deformation ratio μ defined in Equation (3) (i.e., plastic deformation normalised to δy) through numerical analyses with elastic-plastic systems with 5 degrees of freedom. Equation (4) shows that N and the ratio η/μ take the same value. The ratio η/μ is also referred to as the "efficacy" of the structural element. Zahran and Hall's results show that the equivalent number of yield cycles N is a useful parameter to define the severity or damage potential of the ground motion. It is influenced by (i) the characteristics of the ground motion: N increases with duration and tends to decrease for near-field records; (ii) the fundamental period: N increases for structures in the 0.2-2 Hz range; and (iii) ductility: N greatly increases as the displacement ductility increases. The number of equivalent cycles N is sometimes termed as n eq in the literature.
Akiyama [19][20][21] addressed the correspondence between the cumulative inelastic deformation ratio, η defined in Equation (2) (i.e., dissipated hysteretic energy normalised to the yield force Q y and displacement δ y ) and the inelastic deformation ratio µ defined in Equation (3) (i.e., plastic deformation normalised to δ y ) through numerical analyses with elastic-plastic systems with 5 degrees of freedom. Equation (4) shows that N and the ratio η/µ take the same value. The ratio η/µ is also referred to as the "efficacy" of the structural element.
Akiyama highlighted that, if known in advance, the ratio η/µ allows us to estimate the maximum expected displacement for a certain amount of dissipated energy demand. Regardless of any other parameters, he proposed an approximate design value of η/µ = 2 for bare elastic-perfectly plastic frames with no strength deterioration. This expression suggests that the efficacy of the structure (in terms of the relationship between dissipated energy and maximum displacement) can be predicted for a certain lateral carrying system under the design earthquake.
Akiyama also proposed design values [19][20][21] of the ratio η/µ in Equation (5) for flexible-stiff mixed frame multi-storey systems, consisting of (i) a flexible part or frame element able to restrain the lateral deformations within the elastic range for the maximum drift expected under the seismic loading (Figure 2a), and (ii) a stiff part or energy absorbing elements (namely, dampers) with high elastic rigidity and large deformation capacity ( Figure 2b). Figure 2c represents the behaviour of the mixed system. Because of the added complexity of the mixed multi-storey system, a parametric study was used to focus on structures with predetermined lateral stiffness and strength distributions, and changing the parameter r q = f Q m / s Q y , that is, the ratio of the maximum force in the flexible part, f Q m , by the yield strength of the stiff part, s Q y . From the results of these studies, it can be concluded from Equation (5) that as r q increases, the number of equivalent yield cycles in the dampers increases, leading to the following relationship: Appl. Sci. 2019, 9, x FOR PEER REVIEW 3 of 17 Akiyama highlighted that, if known in advance, the ratio η/μ allows us to estimate the maximum expected displacement for a certain amount of dissipated energy demand. Regardless of any other parameters, he proposed an approximate design value of η/μ = 2 for bare elastic-perfectly plastic frames with no strength deterioration. This expression suggests that the efficacy of the structure (in terms of the relationship between dissipated energy and maximum displacement) can be predicted for a certain lateral carrying system under the design earthquake.
Akiyama also proposed design values [19][20][21] of the ratio η/μ in Equation (5) for flexible-stiff mixed frame multi-storey systems, consisting of (i) a flexible part or frame element able to restrain the lateral deformations within the elastic range for the maximum drift expected under the seismic loading (Figure 2a), and (ii) a stiff part or energy absorbing elements (namely, dampers) with high elastic rigidity and large deformation capacity (Figure 2b). Figure 2c represents the behaviour of the mixed system. Because of the added complexity of the mixed multi-storey system, a parametric study was used to focus on structures with predetermined lateral stiffness and strength distributions, and changing the parameter rq = fQm/sQy, that is, the ratio of the maximum force in the flexible part, fQm, by the yield strength of the stiff part, sQy. From the results of these studies, it can be concluded from Equation (5) that as rq increases, the number of equivalent yield cycles in the dampers increases, leading to the following relationship: Malhotra [22] proposed a cyclic demand spectra that specifies the number of cycles for which the seismic load must be resisted by a structure, and noted that both the amplitude and the number of cycles should be considered in assessing the damage potential of ground motions. Malhotra's alternative definition of the number of cycles Ncy is based on dividing cumulative damage by the damage caused by a full cycle of largest amplitude umax. The Coffin-Manson law was used for damage and damage is accumulated by the Miner's rule. The number of equivalent cycles of amplitude umax Malhotra [22] proposed a cyclic demand spectra that specifies the number of cycles for which the seismic load must be resisted by a structure, and noted that both the amplitude and the number of cycles should be considered in assessing the damage potential of ground motions. Malhotra's alternative definition of the number of cycles N cy is based on dividing cumulative damage by the damage caused by a full cycle of largest amplitude u max . The Coffin-Manson law was used for damage and damage is accumulated by the Miner's rule. The number of equivalent cycles of amplitude u max that cause the same damage to the structure as the entire deformation history of n half-cycles of amplitude u i is given by Equation (6), where c is a structural damage parameter that requires experimental calibration. Cyclic demand spectra are site-specific, so they contain all uncertainties related to the ground motion and yield level, so, the only variable is the fundamental period. Exponent c must account for the rest of the structural characteristics. Based on the definition of the number of cycles in Equation (6), three cyclic spectral parameters N A , N V , N D are introduced to characterise cyclic response, depending on whether the structural response is to be governed by acceleration, velocity or displacement histories.
The assumption behind this reasoning is that cyclic response depends primarily on the stiffness of the system. It was found that N A > N V > N D so a stiff structure is likely to experience a larger number of cycles than a flexible structure. Relative stiffness was measured by the ratio of the fundamental period to the earthquake central period.
Kunnath and Chai [23] claimed that the estimation of cyclic deformation demand resulting from earthquake load is crucial if the damage of the system needs to be evaluated, and they produced a cyclic demand spectra to assess the number of full inelastic cycles, N f . The expression for N f in Equation (7) includes an energy shape factor α h that depends on both the ductility and the shape of the hysteresis loops. Structural properties such as Q y and δ y are eliminated from Equation (7) with a focus on seismic design by introducing the relationships between fundamental period T, displacement ductility µ, and the seismic response modification factor R. As a result, the number of inelastic cycles in a SDOF can be calculated with Equation (8) as a function of the hysteretic energy E H , the mass m, the displacement ductility (µ + 1) = δ max /δ y , the force reduction factor R, fundamental period T, energy shape factor α h , and spectral acceleration S a . Use of Equation (8) implies prior knowledge of several variables including ground motion parameters, structural parameters, and target ductility. For each specific structural system, a series of empirical expressions can be calibrated to relate E H , R, α h , and µ, to arrive at simple expressions with only one variable.
Manfredi [24] underlined that the equivalent number of yield cycles varies largely depending on the characteristics of the ground motion, with values varying from 1 for impulsive records and up to 40 for long-duration records. Statistical regressions of nonlinear time history analyses using two components of 122 European ground motion records revealed that the equivalent number of yield cycles can be predicted with the exponential type expression in Equation (9). It was found that (i) the elastic behaviour is well described by the ratio τ of the fundamental period, T, to the corner period within short and medium/long periods in the spectrum, T G , where the spectrum changes from the constant acceleration segment to the constant velocity segment; (ii) the plastic behaviour is described with the strength reduction factor R; and (iii) the seismological parameter with the highest correlation factor is the dimensionless parameter I d = 2g·I A /(π·PGA·PGV), where I A is the Arias intensity, PGA is peak ground acceleration, and PGV is peak ground velocity.
Further research by Manfredi et al. [25] has shown that I d correlates well to the equivalent number of yield cycles N both in the vicinity of the source and far from it. Using regression analyses with 128 near-fault and 122 far-field records, sets of parameters in Equation (9) were calibrated in order to assess the number of yield cycles experienced by a SDOF structure in the far-field and in the near-fault. They observed that N tends to increase with I d and that I d tends to increase with the distance to the earthquake source. These observations led them to conclude that the relative importance of the cyclic demand grows at higher distances from the fault, while structural response is governed by peak demand in the near-fault. Using this equation requires calibration of the regression parameters C, χ, ϕ, and κ.
It can be deduced from Equations (6)-(9) that (i) accurate predictions of the number of inelastic cycles requires calibration of parameters for each specific structural system or hysteretic behaviour, (ii) all equations recognise the influence of the yield level (either via R or µ), and (iii) impulsive near-fault ground motions has been found to significantly reduce the number of inelastic cycles. This last point is of paramount importance in the design of structures with hysteretic dampers in the near-fault. Because the response of hysteric dampers is displacement dependent, the decrease in the number of yield cycles necessarily reduces the efficacy of the dampers.
The purpose of this investigation is to highlight the relevant parameters to predict the equivalent number of yield cycles (η/µ) of hysteretic dampers installed in reinforced frames within the framework of performance-based seismic design, with explicit consideration of near-fault effects. To this end, this paper analyses the seismic response of twelve representative reinforced concrete (RC) frames equipped with hysteretic dampers in a region with medium to high seismicity in Europe. The numerical models account for three building heights (3, 6, 9 stories), the collapse pattern of the bare RC frames (with or without capacity design), and two setups of hysteretic dampers (designed for near-fault or far-field seismicity).

Materials and Methods
The methodology is based on nonlinear time history analyses of 12 numerical models of RC frames equipped with hysteretic dampers. These analyses were performed using a research code to implement the stiffness method. The solver uses the well-known Newmark-beta method. The nonlinear behaviour was accounted for using a hysteretic model for RC sections with normal cyclic degradation, and a hysteretic model for yielding dampers accounting for the Bauschinger effect. The response of each model was analysed under two seismic scenarios, far-field (FF) or near-fault (NF). The main differences among the models were: (i) building height, 3, 6, and 9 storeys, and (ii) capacity design (CD), that is, enforcing or not of capacity design rules to achieve a weak beam strong column lateral collapse pattern in the RC frame. Capacity design rules to achieve a weak-beam strong-column lateral collapse pattern consist of a ratio of ultimate moment capacity of columns to beam that is larger than 1.71 [21]. Table 1 shows a summary of the 12 models, where each model is labelled as follows, number of storeys-seismic scenario-capacity design (e.g., 3 storeys under near-fault and designed without capacity rules would be 3NF. The same, with capacity rules would be 3NFCD) Table 1. Summary of buildings models.

Storeys
Far-Field Seismicity Near-Fault Seismicity

Design of Prototype Structures
Representative frames were extracted from 3, 6, and 9 storey prototype residential buildings with three 5.1 m bays and storey heights of 3.5 m in the lower floor and 3.1 m for the rest of the storeys. Figures A1-A3 in Appendix A show the dimensions and reinforcement of the frame elements. RC frames were designed to the limit states, with concrete C25/30 and steel B500S following the Eurocode 8 conventional force-based modal response spectrum analysis, with design ground acceleration a g = 0.23 g, soil type C and a behaviour factor q = 3. The frames lateral behaviour is characterised by the parameters in Table 2, that is, the storey lateral stiffness f k, and the storey lateral yield displacement f δ y . Hysteretic dampers in the RC structures are designed so that the RC frame remains elastic under the action of the design earthquake and damage concentrates on the dampers. This is achieved by applying an energy-based method [17] that spreads damage evenly in the dampers along the height of the building. In this procedure, the design earthquake is characterised in terms of a bilinear velocity spectrum with equivalent velocity V E = 112.7 cm/s and V D = 72 cm/s, and a number of seismological parameters: the initial period of the medium period range T NH , the main period of the ground motion is T G , the dimensionless parameter I d , and the exponential regression parameters c 1 and c 2 as formulated by Manfredi et al. [25].  Table 3, which shows the lateral stiffness s k and lateral yield displacement s δ y of the dampers at each storey level.

Numerical Models and Analysis
The RC frames with dampers were modelled as two nonlinear springs connected in parallel with lumped masses at each storey level as shown in Figure 3. Beams and columns in the frame (referred to as flexible elements) are idealised as rigid bars connected by rotational springs or hinges. The nonlinear behaviour of the hinges is given by elastic-perfectly plastic envelopes that are derived from their cross-section dimensions and reinforcement. Four parameters control the cyclic behaviour of the hinges to model the effects of pinching (HS = 0.50), stiffness degradation (HC = 2.00), and strength degradation (HBE = 0.10 and HBD = 0.10) in RC [26,27]. These parameters correspond to normal degradation of reinforced concrete elements. Hysteretic dampers are referred to as stiff elements, and their nonlinear load-deformation behaviour is idealised by a trilinear envelope proportional to the elastic stiffness s k, and by two parameters (α = 0.66, β = 0.61) that account for the Bauschinger effect in steel. Damper properties are taken from available component testing [28]. An example of the hysteretic model of the dampers is depicted in Figure A4, which shows the nonlinear behaviour of a damper under incremental cyclic displacement.
Strong Motion Database [29]. One set includes 22 far-field records in Table A2, and the other set includes other 22 near-fault or pulse-like records in Table A1. The records are arranged in Appendix A by the parameter Id, which is representative of the intensity of velocity pulses in the ground motion records. These accelerograms are scaled iteratively in amplitude to two seismic hazards levels: the first hazard level corresponds to the design earthquake and is characterised by an energy input in terms of equivalent velocity VE = 112.7 cm/s. The scaling of the ground motion changes for each prototype structure as is shown in the elastic response spectra at 5% damping in Figure A5. The second hazard level corresponds to an extreme earthquake that takes the structure to the onset of strength degradation due to the capacity loss of one of the dampers.

Results
The results for the two sets of analyses are presented separately. The first set corresponds to the design earthquake. At this hazard level, seismic damage and inelastic energy dissipation is sustained by the dampers and the results correlate to the efficacy sη/sμ of the dampers with significant structural and seismological parameters. The second sets of analyses correspond to the ultimate capacity of the mixed structural system under the action of very rare earthquakes. In this case, the results are upper bound values of the efficacy of the two parts of the structural system.

Estimation of Average sη/sμ of Dampers under Design Earthquake Hazard Level
The results at the end of the design earthquake hazard level showed that the inelastic demand concentrated at the dampers while the structure remained elastic, in both the with and without capacity design cases. Therefore, just the plastic response of the dampers is of interest herein. For each nonlinear time history analysis, the following variables were collected in a dataset: the ratio τ of the fundamental period T1 to the initial period of medium Tg, the reduction factor R, the seismological parameter Id, and the average cyclic demand parameter sη/sμ along the height of the building. Regular statistical procedures were used to analyse the results and find significant relationships between the variables. A preliminary principal component analysis was performed to determine the minimum number of variables that account for the most variation in the data, and found that the first two components (Id and R) cover 92% of the variability of the data. On the other hand, a poor correlation between the parameter τ and the average value of sη/sμ was found. Figure 4 shows a scatter of these principal variables (R-1) and Id versus the average cyclic demand parameter sη/sμ for each analysis in the dataset. Data is arranged around two clusters: near-fault seismicity in red, and far-field seismicity in blue. Seismic action is represented by selecting two sets of ground motion records from the European Strong Motion Database [29]. One set includes 22 far-field records in Table A2, and the other set includes other 22 near-fault or pulse-like records in Table A1. The records are arranged in Appendix A by the parameter I d , which is representative of the intensity of velocity pulses in the ground motion records. These accelerograms are scaled iteratively in amplitude to two seismic hazards levels: the first hazard level corresponds to the design earthquake and is characterised by an energy input in terms of equivalent velocity V E = 112.7 cm/s. The scaling of the ground motion changes for each prototype structure as is shown in the elastic response spectra at 5% damping in Figure A5. The second hazard level corresponds to an extreme earthquake that takes the structure to the onset of strength degradation due to the capacity loss of one of the dampers.

Results
The results for the two sets of analyses are presented separately. The first set corresponds to the design earthquake. At this hazard level, seismic damage and inelastic energy dissipation is sustained by the dampers and the results correlate to the efficacy s η/ s µ of the dampers with significant structural and seismological parameters. The second sets of analyses correspond to the ultimate capacity of the mixed structural system under the action of very rare earthquakes. In this case, the results are upper bound values of the efficacy of the two parts of the structural system.

Estimation of Average s η/ s µ of Dampers under Design Earthquake Hazard Level
The results at the end of the design earthquake hazard level showed that the inelastic demand concentrated at the dampers while the structure remained elastic, in both the with and without capacity design cases. Therefore, just the plastic response of the dampers is of interest herein. For each nonlinear time history analysis, the following variables were collected in a dataset: the ratio τ of the fundamental period T 1 to the initial period of medium T g , the reduction factor R, the seismological parameter I d , and the average cyclic demand parameter s η/ s µ along the height of the building. Regular statistical procedures were used to analyse the results and find significant relationships between the variables. A preliminary principal component analysis was performed to determine the minimum number of variables that account for the most variation in the data, and found that the first two components (I d and R) cover 92% of the variability of the data. On the other hand, a poor correlation between the Appl. Sci. 2020, 10, 4364 8 of 16 parameter τ and the average value of s η/ s µ was found. Figure 4 shows a scatter of these principal variables (R-1) and I d versus the average cyclic demand parameter s η/ s µ for each analysis in the dataset. Data is arranged around two clusters: near-fault seismicity in red, and far-field seismicity in blue. It can be observed that far-field seismicity produces significantly higher values of sη/sμ. In Figure  4, the results from building models incorporating (or not) capacity design rules are plotted with triangles pointing downward or upwards respectively, and show very little variation one to each other. Using a non-linear least squares fitting, the average cyclic demand parameters of the dampers in a structure can be predicted with Equation (10) as a function of the principal variables, the strength reduction factor R and the seismological parameter Id. The prediction with this equation is depicted in Figure 4 by the grey surface and black dots. Equation (10) correlates moderately with a regression score R 2 = 0.59.

Behaviour of sη/sμ of Dampers and fη/fμ of Frames under the Maximum Considered Earthquake
In this study, the maximum considered earthquake that the structures were able to withstand was defined as that in which one of the dampers in the building reaches its ultimate energy dissipation capacity. The ultimate energy dissipation capacity was assessed using the damage model proposed by [30] and calibrated with the experimental results from a previous study [28]. At this earthquake level, plastic excursions are likely to occur on the reinforced concrete frame as well. For each nonlinear time history analysis under the maximum earthquake, the following variables were collected in a dataset: the parameter Id of the ground motion, the number of equivalent inelastic cycles in the dampers sη/sμ, and the average fη/fμ in the reinforced concrete frame. Figure 5 shows the number of equivalent inelastic cycles for dampers and reinforced concrete frames against the seismological parameter Id. Results for the hysteric dampers are shown with triangular markers with error bars (±1 standard deviation). In the case of frame elements, the average results are shown with triangular markers without error bars. The markers are pointing downwards for models that incorporate capacity design.
Under the maximum considered earthquake, the ratio sη/sμ in the dampers showed a tendency to increase with the parameter Id, yielding lower values for near-fault records, and significantly larger values for far-field records. The maximum average values obtained for far-field records are about twice as much as the maximum values in the near-fault (14 vs. 7). The minimum average values are also double in the far-field (4 vs. 2). The average number of equivalent inelastic cycles in the reinforced concrete frame fη/fμ remains constant in the range 1-2, suggesting that there is no It can be observed that far-field seismicity produces significantly higher values of s η/ s µ. In Figure 4, the results from building models incorporating (or not) capacity design rules are plotted with triangles pointing downward or upwards respectively, and show very little variation one to each other. Using a non-linear least squares fitting, the average cyclic demand parameters of the dampers in a structure can be predicted with Equation (10) as a function of the principal variables, the strength reduction factor R and the seismological parameter I d . The prediction with this equation is depicted in Figure 4 by the grey surface and black dots. Equation (10) correlates moderately with a regression score R 2 = 0.59. s η/ s µ = 1 + 1.18· (R − 1) 0.24 · I 0.44 d (10)

Behaviour of s η/ s µ of Dampers and f η/ f µ of Frames under the Maximum Considered Earthquake
In this study, the maximum considered earthquake that the structures were able to withstand was defined as that in which one of the dampers in the building reaches its ultimate energy dissipation capacity. The ultimate energy dissipation capacity was assessed using the damage model proposed by [30] and calibrated with the experimental results from a previous study [28]. At this earthquake level, plastic excursions are likely to occur on the reinforced concrete frame as well. For each nonlinear time history analysis under the maximum earthquake, the following variables were collected in a dataset: the parameter I d of the ground motion, the number of equivalent inelastic cycles in the dampers s η/ s µ, and the average f η/ f µ in the reinforced concrete frame. Figure 5 shows the number of equivalent inelastic cycles for dampers and reinforced concrete frames against the seismological parameter I d .
Results for the hysteric dampers are shown with triangular markers with error bars (±1 standard deviation). In the case of frame elements, the average results are shown with triangular markers without error bars. The markers are pointing downwards for models that incorporate capacity design. correlation with the seismological parameter Id. Regarding the lateral collapse pattern of the frames, models of the weak-beam strong-column type are shown with triangular markers that point down, while frames with other collapse patterns are shown with markers that point up. Given the scatter of the results, it seems that the values of both ratios, sη/sμ and fη/fμ are not influenced by capacity design rules being enforced or not.

Discussion
The number of equivalent inelastic cycles N can be formulated as the ratio η/μ between the normalised dissipated energy η and the displacement ductility μ. This ratio is crucial to the design of hysteric dampers within the framework of performance-based design, as it quantifies the capacity to dissipate seismic energy for a given displacement. The ratio η/μ is also referred to as the efficacy of the dampers. The literature shows that the number of equivalent inelastic cycles varies for each specific lateral carrying system, and that it is strongly influenced by seismological parameters, especially by the intensity and the impulsive effects of the ground motion [16,17,25]. This study focused on mixed flexible-stiff structures, made of reinforced concrete frames with supplemental hysteric dampers under two levels of seismic hazard: the design earthquake and the maximum considered earthquake.
One interesting result is that the capacity design provisions in the gravity load resisting system, i.e., reinforced concrete frames has a negligible effect on the number of equivalent inelastic cycles of the dampers for both hazard levels. Because the efficacy of the dampers is independent of the lateral collapse pattern, the frame design can be neglected when designing the dampers. This is important in the retrofitting of pre-code frames that have undefined lateral collapse mechanisms. In light of these results, retrofitting strategies that focus on increasing the lateral deformability of the frames or preventing brittle failures (such as FRP wrapping) will have a greater impact on the performance of the dampers compared to strategies that focus on strengthening or reallocating plastic hinges.
Under the design earthquake, the ratio sη/sμ can be predicted as a function of the structural parameter R (the seismic force reduction factor, which relates the intensity of the ground motion and Under the maximum considered earthquake, the ratio s η/ s µ in the dampers showed a tendency to increase with the parameter I d , yielding lower values for near-fault records, and significantly larger values for far-field records. The maximum average values obtained for far-field records are about twice as much as the maximum values in the near-fault (14 vs. 7). The minimum average values are also double in the far-field (4 vs. 2). The average number of equivalent inelastic cycles in the reinforced concrete frame f η/ f µ remains constant in the range 1-2, suggesting that there is no correlation with the seismological parameter I d . Regarding the lateral collapse pattern of the frames, models of the weak-beam strong-column type are shown with triangular markers that point down, while frames with other collapse patterns are shown with markers that point up. Given the scatter of the results, it seems that the values of both ratios, s η/ s µ and f η/ f µ are not influenced by capacity design rules being enforced or not.

Discussion
The number of equivalent inelastic cycles N can be formulated as the ratio η/µ between the normalised dissipated energy η and the displacement ductility µ. This ratio is crucial to the design of hysteric dampers within the framework of performance-based design, as it quantifies the capacity to dissipate seismic energy for a given displacement. The ratio η/µ is also referred to as the efficacy of the dampers. The literature shows that the number of equivalent inelastic cycles varies for each specific lateral carrying system, and that it is strongly influenced by seismological parameters, especially by the intensity and the impulsive effects of the ground motion [16,17,25]. This study focused on mixed flexible-stiff structures, made of reinforced concrete frames with supplemental hysteric dampers under two levels of seismic hazard: the design earthquake and the maximum considered earthquake.
One interesting result is that the capacity design provisions in the gravity load resisting system, i.e., reinforced concrete frames has a negligible effect on the number of equivalent inelastic cycles of the dampers for both hazard levels. Because the efficacy of the dampers is independent of the lateral collapse pattern, the frame design can be neglected when designing the dampers. This is important in the retrofitting of pre-code frames that have undefined lateral collapse mechanisms. In light of these results, retrofitting strategies that focus on increasing the lateral deformability of the frames or preventing brittle failures (such as FRP wrapping) will have a greater impact on the performance of the dampers compared to strategies that focus on strengthening or reallocating plastic hinges.
Under the design earthquake, the ratio s η/ s µ can be predicted as a function of the structural parameter R (the seismic force reduction factor, which relates the intensity of the ground motion and the lateral capacity of the structure), and the seismological parameter I d , which has been shown to be an indicator of the impulsive effects of the ground motion. The tendency in the results is for larger ratios of s η/ s µ as the parameters R and I d increase. This latter parameter I d accounts for most of the variability in the results, therefore it must drive the design of hysteretic dampers. Other parameters such as the ratio τ between the fundamental period and ground motion corner period show a poor correlation with s η/ s µ in this type of structure and were neglected in the regression.
The maximum considered earthquake takes the structure to the onset of strength degradation in one of the dampers along the building, and allows for some structural damage to the reinforced concrete frames. Because some dampers reach their ultimate dissipation capacity in this scenario, the results of the ratio s η/ s µ in the dampers should be taken as an upper-bound limit in the damper design. Our study confirmed that the efficacy of the dampers is significantly larger and more stable in the far-field. However, the ratio f η/ f µ in the frame remains in the range of 1~2 for all analyses, which agrees with Akiyama's [19] approximate design value of η/µ = 2 for bare elastic-plastic frames.
Interestingly, the ratio f η/ f µ is unchanged despite the impulsive effects of the ground motions, as it shows non-correlation to I d . Because the efficacy of the frames is similar for far-field and near-fault seismicity, the matter of impulsivity can be neglected in the overall design and assessment of frames equipped with dampers.