Energy-based design criterion of dissipative bracing systems for the seismic retrofit of frame structures

Direct sizing criteria represent useful tools in the design of dissipative bracing systems for the advanced seismic protection of existing frame structures, especially when incorporated dampers feature a markedly non-linear behaviour. An energy-based procedure is proposed herein to this aim, focusing attention on systems including fluid viscous devices. The procedure starts by assuming prefixed reduction factors of the most critical response parameters in current conditions, which are evaluated by means of a conventional elastic finite element analysis. Simple formulas relating the reduction factors to the equivalent viscous damping ratio of the dampers, ξeq, are proposed. These formulas allow calculating the ξeq values that guarantee the achievement of the target factors. Finally, the energy dissipation capacity of the devices is deduced from ξeq, finalizing their sizing process. A detailed description of the procedure is presented in the article, by distinguishing the cases where the prevailing structural deficiencies are represented by poor strength of the constituting members, from the cases having excessive horizontal displacements. A demonstrative application to the retrofit design of a reinforced concrete gym building is then offered to explicate the steps of the sizing criterion in practice, as well as to evaluate the enhancement of the seismic response capacities generated by the installation of the dissipative system.


Introduction
Passive energy dissipation and base isolation systems are increasingly adopted in the anti-seismic design of new frame structures as well as to retrofit existing ones.Remarkable innovations have been noticed in the last decade for both strategies.Concerning base isolation, single and multiple curved surface sliders [1][2][3][4] and roller-type isolators [5][6][7] became a competitive alternative to the widely used class of high damping rubber bearings.At the same time, new technologies, or improved versions of existing ones, have been implemented in the field of energy dissipation, capable of supplying supplemental damping and horizontal stiffness in different proportions, depending on the mechanical characteristics of dampers and their installation layout.Among the latest achievements in this field, systems incorporating updated models of metallic yielding devices like ADAS (Added Damping and Stiffness) and TADAS (Triangular Added Damping and Stiffness) elements [8][9][10][11][12], shear panels [13][14][15][16][17][18], shear panels with openings [19][20][21][22][23], or buckling-restrained braces [24][25][26][27], have been proposed, which typically provide significant contributions in terms of both properties.Viscous dampers can offer stiffening effects too, in addition to supplemental damping [28][29][30].At the same time, special types of viscous dampers, when mounted at the tip of supporting braces in parallel with the overlying beam axis [31,32], slightly increase the horizontal stiffness of the structural system, while supplying high additional damping.
Differently from the base isolation systems, whose preliminary design is straightforward, a further spreading of dissipative bracing technologies in the professional community strongly depends on the availability of intuitive procedures, especially concerning the sizing of dampers.The first methods offered in literature start from setting the desired damping ratio (i.e., the ratio of the damping coefficient to the critical damping coefficient) in the fundamental mode of vibration of the structure, in the hypothesis that the relevant effective modal mass (EMM) is a predominant portion of the total seismic mass [33][34][35].In general, the practical application of these methods consists in examining the response spectra at various damping ratios and choosing the value that allows constraining the maximum "global" response parameters (base shear, top lateral displacements, etc.) within targeted acceptable limits.When the devices are characterized by nonlinear viscous properties, the same objectives can be reached by transforming relevant damping capacity into equivalent linear viscous coefficients [28,33].These studies have provided the basis for the design procedures of buildings incorporating passive energy dissipation systems included in ASCE 41-06 [36] and  Standards.
Along the same conceptual line, some procedures based on the use of normative response spectra scaled by reduction factors corresponding to the damping capacity of the devices have been proposed more recently [38,39], where reference is made to damping ratio values no greater than 0.3.Other approaches use equivalent linear or non-linear static analyses to evaluate the design actions and reduce their effects through added damping [40].All the above-mentioned procedures are conceived for substantially regular structures.Few solutions are found for problems characterized by significant irregularities in plan and/or in elevation.Among these, a method based on properly calibrated expressions of the damping ratio derived from the results of non-linear dynamic analyses is formulated in [41,42].
Alternative approaches refer to the energy balance computed from time-history analyses [43,44], or optimization procedures [45].Among the energy-based design methods, a criterion determining the minimum damping coefficients of the devices required to assign them the capability of dissipating a prefixed fraction β of the seismic input energy E I computed on each story [31,32] or the entire structure [12] was proposed by the author and co-authors.To facilitate the choice of β, preferable ranges were provided for several different structural types and checked in relation to the assumed design targets [12,31,32].However, as the method requires a preliminary evaluation of the seismic input energy demand on the original structure, a finite element time-history analysis must be carried out first, and E I post-calculated from the results.Although an energy calculation can be performed with the help of commercial finite element programs by means of simple input instructions, professional engineers are not always familiar with this design approach and may be discouraged from using it.
In view of this, a new procedure that bypasses this initial step by directly estimating the minimum damping capacity to be assigned to the dampers is proposed in this paper, where attention is focused on the retrofit design of reinforced concrete (RC) structures.The procedure starts by assuming prefixed α reduction factors of the most critical response parameters in current conditions, which are evaluated by means of a conventional elastic finite element analysis.Simple formulas relating the α factors to the equivalent viscous damping ratio of the dampers, ξ eq , are proposed.These formulas allow calculating the ξ eq values guaranteeing the achievement of the target α factors.Finally, the energy dissipation capacity of the devices is deduced from ξ eq , finalizing their sizing process.
A detailed description of the method is presented in the next Section.A demonstrative application to the retrofit design of an RC gym building is then offered to explicate the relevant steps in practice.Finally, a performance assessment analysis of the structure in original and retrofitted conditions is carried out to evaluate the enhancement of seismic response capacities produced by the incorporation of the dissipative bracing system.

Design Procedure
The design procedure is based on the assumption that, for relatively stiff frame structures, a substantial improvement of seismic performance can be reached by incorporating a supplemental damping system with limited stiffening capacity.For more deformable structures, a supplemental stiffness contribution helps control the lateral displacements better, which prevents an overdissipation demand to the damping technology adopted.
In the following, attention is focused on the first class of frame structures, which can be identified with low-and low-to-medium-rise multistorey building structures typical of European countries, i.e., with a height below 20 m.In order to obtain for this class of structures an approximate evaluation of the maximum fundamental translational vibration period T 1 , the simple relation provided by Eurocode 8 [46] and adopted by several national Standards, among which the Italian Technical Standards [47], can be used: , where the C t coefficient is set as equal to 0.075 for steel structures and 0.085 for RC ones, and H is the height of the structure measured from the foundation, expressed in meters.For H = 20 m, T 1 = 0.71 s (steel) and T 1 = 0.8 s (RC) is obtained.By jointly considering the two values, T 1 = 0.8 s is assumed in the following as the approximate upper limit of the fundamental period for relatively stiff frame building structures.
As observed above, for these structures the retrofit design objectives can be met essentially by an added damping, with a marginal role of supplemental stiffness.This prompts to select bracing systems that incorporate pressurized fluid viscous (FV) spring-dampers as a protective strategy.Indeed, the spring component of the devices, acting in series with the steel braces, remarkably constrains the total horizontal translational stiffness of these systems.This results in a moderate contribution (normally below 10%) to the lateral stiffness of the retrofitted structures.
The α reduction factors are calibrated on the mechanical characteristics of this type of devices.For dampers adding a non-negligible supplemental stiffness, the procedure follows the same steps, but different criteria for the evaluation of α are required, which will be detailed in a further step of this research.
In view of the practical application of the design method, a distinction is made between the cases where the most significant structural deficiencies are represented by poor shear or bending moment strength of columns and/or beams-where α is intended as "stress-related" reduction factor, named α s in the following-and the cases having excessive horizontal displacements, where α is intended as "deformation-related" reduction factor, α d .Although the two types of lacks normally coexist, the former is frequently prevailing in structures with undersized member sections, albeit stiffened by structural elements (e.g., RC shear walls) or non-structural components (e.g., masonry infills interacting with the frame members), whereas the latter is likely to prevail when these elements are not present-and thus the structural system highlights high deformability to lateral loads-but the frame members have only moderate deficiencies in terms of sizes and structural details.

Structures with Poor Shear or Bending Moment Strength of Constituting Members
The basic design objective of supplemental damping-based retrofit interventions is always a significant enhancement of the seismic performance as compared to the original structural conditions.When high-damping capacity dissipaters are adopted, like pressurized FV devices, the attainment of the Immediate Occupancy (IO) seismic performance level can be planned up to the maximum considered earthquake (MCE), with reasonable costs and architectural impact.This objective corresponds to keeping the response of the structural members within the elastic field, thus replacing the plastic demand in current conditions with the dissipative action of the protective system.The design procedure is based on this target and articulated in the four steps described below.
1. Development of an elastic finite element verification analysis and relevant stress state checks in current conditions.
The shear and bending moment-related stress states in columns and beams are calculated by means of an elastic finite element (either modal superposition with response spectrum or time-history) analysis of the structure in current conditions.Then, the maximum stress states are checked by comparison with the corresponding elastic limit values for the two axes in plan of the reference Cartesian coordinate system, X and Y (Z being the vertical axis).The shear and bending moment values in the member that proves to be subjected to the highest unsafe stress conditions are named V a j,c and M a j,c , where index j refers to the axis (j = X,Y), c denotes the "critical" member(s), and a means that the values are derived from the elastic assessment analysis.
2. Evaluation of the target stress reduction factors α s .This evaluation is carried out by referring to the scheme in Figure 1a, where d, s represent the general deformation-related and stress-related parameters for the critical member(s).The objective of the retrofit design consists in reducing the seismic demand in current state calculated by the elastic analysis in where d e , s e are the elastic limit deformation-related and stress-related parameters for the critical member(s), respectively.On the contrary, in a traditional design approach, the reduction of the force-related parameter is pursued by exploiting the ductile plastic response of the structural members, symbolically represented by the dashed line K-B, by moving the maximum response point from A to K, with coordinates (d A , s B ).
Hence, the force reduction factor α s , targeted to reach an elastic response of the critical member (and thus of all remaining members) when passing from original to retrofitted conditions, is given by: The s A and s B parameters in Figure 1 are detailed below according to the specific lack of strength affecting the structural members in original conditions.The shear and bending moment-related stress states in columns and beams are calculated by means of an elastic finite element (either modal superposition with response spectrum or time-history) analysis of the structure in current conditions.Then, the maximum stress states are checked by comparison with the corresponding elastic limit values for the two axes in plan of the reference Cartesian coordinate system, X and Y (Z being the vertical axis).The shear and bending moment values in the member that proves to be subjected to the highest unsafe stress conditions are named a c j, V and a c j, M , where index j refers to the axis (j = X,Y), c denotes the "critical" member(s), and a means that the values are derived from the elastic assessment analysis.

Evaluation of the target stress reduction factors αs.
This evaluation is carried out by referring to the scheme in Figure 1a, where d, s represent the general deformation-related and stress-related parameters for the critical member(s).The objective of the retrofit design consists in reducing the seismic demand in current state calculated by the elastic analysis in step 1-denoted by point A with coordinates (dA = dmax, sA = smax), where dmax, smax are the maximum d and s computed values-below point B with coordinates (dB = de, sB = se), where de, se are the elastic limit deformation-related and stress-related parameters for the critical member(s), respectively.On the contrary, in a traditional design approach, the reduction of the force-related parameter is pursued by exploiting the ductile plastic response of the structural members, symbolically represented by the dashed line K-B, by moving the maximum response point from A to K, with coordinates (dA, sB).
Hence, the force reduction factor αs, targeted to reach an elastic response of the critical member (and thus of all remaining members) when passing from original to retrofitted conditions, is given by: The sA and sB parameters in Figure 1 are detailed below according to the specific lack of strength affecting the structural members in original conditions.2.a Lack of shear strength.
In this hypothesis, s A is the shear force in the critical member V a j,c , according with the notation introduced in step 1. Said V e j,c = s B the elastic limit shear of this member, α s -specified in this hypothesis as α V j -is evaluated as follows: 2.b Lack of bending moment strength.
Similarly to point 2.a, when the critical stress-related parameter s A is the bending moment M a j,c (associated with the concurrent axial force N c if the critical member is a column), said M e j,c = s B the corresponding elastic limit moment, the reduction factor α s -denoted with symbol α Mj in this case-is given by: 2.c Passage from member to storey (α s →α F ).
The A→B transition in the most critical member, governed by factor α V j or α Mj , implies a similar shift in the response of the frame structure storey to which the member belongs.This is visualized in the graph of Figure 1b, where the storey shear, F, and the interstorey drift, ID, are assumed as response parameters, and the member response points A and B are replaced by the corresponding storey response points A , with coordinates (F A = F max , ID A = ID max ), and B , with coordinates (F B = F e , ID B = ID e ), where the indexes "max" and "e" denote the maximum response value and the corresponding elastic limit in this case too.Based on this correlation, the storey response points A , B in Figure 1b are reached when the critical member attains points A, B in Figure 1a.Therefore, the reduction factor at the storey level, α F : coincides with α V j or α Mj , depending on the lack of strength assessed in the most critical member.
3. Correlation of α F to the dissipated energy E D and the equivalent viscous damping ratio ξ eq of the spring-dampers.
The correlation is established consistently with the design objective of reducing the storey response from point A to point B , thanks to the incorporation of the protective system.The sizing of the devices is performed first by evaluating their equivalent viscous damping ratio ξ eq , using the general expression: where E D is the energy dissipated by the set of FV spring-dampers installed on the storey containing the critical member, and E e is the strain energy of the system estimated by referring to the global response cycle of the set of devices, schematically drawn in Figure 2. Therein, d dmax is the maximum displacement, F D is the damping force component, F e is the elastic reaction force, and F ed = F e + F D is the total reaction force.In order to meet the design objective of passing from point A' to point B', FD is set as equal to the difference between FA' and FB': In doing so, the reduction of storey response is totally assigned to the dissipative action of the protective system, as planned by this retrofit strategy.When substituting ( 4) in ( 6), the following FD In order to meet the design objective of passing from point A to point B , F D is set as equal to the difference between F A and F B : In doing so, the reduction of storey response is totally assigned to the dissipative action of the protective system, as planned by this retrofit strategy.When substituting (4) in ( 6), the following F D expression is obtained: from which it follows: The F e elastic force component is set as equal to the elastic limit value of the storey shear, F B .Therefore, F ed can be alternatively expressed as a function either of F D : or of F e , by substituting ( 8) in ( 9): By referring to the response cycle in Figure 2 and relevant parameters, the following expressions of E D and E e are derived: By substituting ( 11) and ( 12) in ( 5), the following ξ eq expression is deduced: By introducing ( 8) and ( 9), ( 13) becomes: which allows directly quantifying the equivalent viscous damping ratio demanded to obtain the targeted elastic response up to the MCE, simply as a function of the reduction factor α F .Moreover, by solving (5) to express E D as a function of ξ eq and considering (10), the following relation is obtained: 4. Evaluation of the energy dissipation capacity of the dampers, E D , and selection of the devices with the nearest mechanical characteristics.
The basic design objective of reaching an elastic response of the most critical member(s), and thus of the relevant storey, implies that the maximum storey drift in retrofitted conditions, ID max , is constrained below the corresponding elastic limit ID B = ID e displayed in Figure 1b.ID max is given by the sum of d dmax and the interstorey drift contribution provided by the braced structure.The latter is normally small, because relatively stiff structures-i.e., with T 1 below about 0.8 s-are dealt with, as mentioned above, and also because the bracing system produces in any case a stiffening effect, although limited by the special installation layout of the FV spring-dampers.In view of this, in order to quickly pre-estimate the energy dissipation capacity of the devices, which is a function of d dmax -other than of F D -d dmax is set as equal to ID e at this stage.Based on this assumption, the energy dissipation capacity to be assigned to the FV devices on the considered storey is drawn from ( 11): In order to express E D as a function of structure-related terms only, ( 8) is substituted in ( 16), and gives: Relation ( 15) can be equivalently adopted to this aim: The sizing process of the spring-dampers is concluded-for each storey-by selecting from the manufacturer's catalogue [48] the device with the nearest energy dissipation capacity to the E D value estimated by (17) or (18) and a stroke approximately equal to ID e .

Structures with Excessive InterStorey Drifts
The attainment of an elastic structural response up to the MCE is assumed as the basic design objective also when the poorest response capacity is in terms of storey drifts.This is more likely verified in structures with a translational period T 1 next to the 0.8 s value approximately fixed as the upper limit for the application of the procedure.The design objective is reached by reducing the maximum interstorey drift computed in current conditions, identified by ID A = ID max in Figure 1b, to the corresponding elastic limit ID B = ID e , i.e., by scaling the drift response by a deformation-related reduction factor α d : When an elastic finite element assessment analysis is carried out on the original structure, α d is proportional to α F .For cases where an inelastic evaluation analysis is developed, assuming a typical degrading strength and degrading stiffness post-elastic behaviour-like the one qualitatively schematized by curve B -K in Figure 1b-α d significantly differs from α F .Therefore, the expressions of E D and ξ eq must be reformulated as a function of α d to allow quickly estimating both quantities also for the structures where the interstorey drift is the critical response parameter.
Starting from relation (13), the two changes introduced as compared to Section 2.1 consist in assuming d dmax = ID A − ID B = ID max − ID e , and F ed = F K = F D .The first assumption corresponds to assigning to the dampers the capacity of absorbing the post-elastic displacement demand of the structure computed in original conditions, so as to meet the design objective of limiting its response to the elastic field after retrofit.The second assumption derives from the fact that the displacement performance enhancement must be achieved essentially by means of the dissipative capacity of the FV devices, by neglecting at the sizing stage the slight stiffening effects related to their elastic spring function.Then, according to (19), the following ξ eq expression is obtained: By referring to (13), the corresponding energy dissipation capacity to be assigned to the dampers on the considered storey results as follows: E D,αd = 2πF e ξ eq,αd ID e (21) Similarly to the case of poor strength of the members, the tentative design choice of the device is carried out by referring to the spring-damper with the nearest energy dissipation capacity to E D , estimated by (20), and a stroke not less than (ID max − ID e ).

Geometrical and Structural Characteristics of the Case-Study Building
The case study examined for a first demonstrative application of the design criterion is the gym in a school in Florence, built in 1971, two external views and an internal view of which are displayed in Figure 3. Figures 4 and 5 show the structural plan and the transversal cross section, respectively.The reference Cartesian coordinate system assumed in the analyses is indicated in Figures 4 and 5, too.As highlighted in Figure 4, the plan is rectangular, with sides of 14.4 m in the transversal direction, parallel to X, and 24.25 m in the longitudinal direction, parallel to Y.The height of the roof top is equal to 9.17 m, whereas the height of the façades, measured on top of the end section of the roof beams, is equal to 8.67 m.The structure is constituted by nine identical RC frames of two columns each, numbered C1-C2 through C17-C18 in Figure 4, placed at a mutual distance of 3 m.The cross sections of the beams and columns and relevant reinforcement details, redrawn from the original structural design drawings, are illustrated in Figure 6.
The columns have a mutual rectangular section with sides of 700 mm along X and 250 mm along Y.The roof beams have a rectangular section with a base of 250 mm and a height varying from 1030 mm, at the ends, to 1540 mm, at half-span.In the longitudinal direction, the columns are connected on top by a rectangular beam, named TB in Figure 5, with dimensions of 250 mm × 400 mm; at a height of 3.17 m, by an intermediate beam, named IB, having a polygonal section with a base of 700 mm and maximum lateral side of 500 mm; and at the base, by a rectangular beam constituting the lateral edge beam of the ground floor, named LEB, with dimensions of 200 mm × 700 mm.The IB beam, which supports the curtain wall-type glazed portions of the façades, subdivides all columns in two levels along the height.The structure of the roof floor is 160 mm thick and made of 120 mm-high and 100 mm-wide partly prefabricated RC joists, parallel to Y and placed at a mutual distance of 400 mm, clay lug bricks, and a 40 mm-thick upper RC slab.The ground floor only differs for the height of the joists, equal to 160 mm, which determines a total thickness of 200 mm.The foundation consists of a 400 mm-thick continuous slab, with 1400 mm-high (slab thickness included) and 250 mm-wide transversal rib beams, which connect the column base sections in X direction and support the ground floor.
A selective investigation campaign was carried out on materials and structural members, including on-site Son-Reb, pacometric and Vickers-type micro-durometer analyses, and laboratory tests on concrete and steel bar samples.On the basis of the prescriptions of Italian Standards [47,49], as well as of professional protocols [50,51], the tests met the basic knowledge level (named LC1 in [49]) for the structural assessment analysis of public buildings in Italy.The corresponding value of the confidence factor, i.e., the additional knowledge level-related safety coefficient to be introduced in stress state checks, is equal to 1.35.The following main properties resulted from the characterization tests: mean cubic compressive strength of concrete equal to 19.6 N/mm 2 ; yield stress and limit stress of steel equal to 417 MPa and 594 MPa, respectively.The total seismic weight of the building is equal to 2960 kN.A selective investigation campaign was carried out on materials and structural members, including on-site Son-Reb, pacometric and Vickers-type micro-durometer analyses, and laboratory tests on concrete and steel bar samples.On the basis of the prescriptions of Italian Standards [47,49], as well as of professional protocols [50,51], the tests met the basic knowledge level (named LC1 in [49]) for the structural assessment analysis of public buildings in Italy.The corresponding value of the confidence factor, i.e., the additional knowledge level-related safety coefficient to be introduced in stress state checks, is equal to 1.35.The following main properties resulted from the characterization tests: mean cubic compressive strength of concrete equal to 19.6 N/mm 2 ; yield stress and limit stress of steel equal to 417 MPa and 594 MPa, respectively.The total seismic weight of the building is equal to 2960 kN.

Verification Analysis in Current Conditions (Step 1 of the Design Procedure)
The verification enquiry in current conditions, constituting Step 1 of the design procedure, is articulated in a modal analysis to calculate the vibration periods and associated modal masses, and in a time-history analysis to assess the seismic performance in terms of stress states and displacements.

Verification Analysis in Current Conditions (Step 1 of the Design Procedure)
The verification enquiry in current conditions, constituting Step 1 of the design procedure, is articulated in a modal analysis to calculate the vibration periods and associated modal masses, and in a time-history analysis to assess the seismic performance in terms of stress states and displacements.

Modal Analysis
The finite element model of the structure, a perspective view of which is displayed in Figure 7, was generated by the SAP2000NL calculus program [52] using frame type elements for all structural members.The modal analysis carried out by the model shows two first horizontal translational modes along X and Y, with vibration periods of 0.89 s (Y) and 0.35 s (X), respectively, and effective modal mass (EMM) equal to 79% along Y and 88.1% along X.The fourth and fifth mode are translational along X and Y too, with periods of 0.26 s (Y) and 0.11 s (X), and EMM equal to 20.9% (Y) and 11.8% (X), which provide summed modal masses with the corresponding first mode-related EMMs nearly equal to 100%, along both axes.The third and sixth mode are purely rotational around the vertical axis Z, with periods of 0.34 s and 0.04 s, and EMM equal to 84.4% and 12.4%, giving a summed modal mass of 96.8%.
The modal parameters quantitatively confirm a notably different translational behaviour of the structure along the two directions in plan, as a consequence of the markedly different sides of the columns along X and Y and the much higher flexural stiffness of the roof beams in comparison to the longitudinal beams.

Time-History Verification and Performance Assessment Analysis
The performance evaluation enquiry was carried out for the four reference seismic levels fixed in the Italian Standards [47], that is, Frequent Design Earthquake (FDE, with 81% probability of being exceeded over the reference time period VR); Serviceability Design Earthquake (SDE, with 50%/VR probability); Basic Design Earthquake (BDE, with 10%/VR probability); and Maximum Considered Earthquake (MCE, with 5%/VR probability).The VR period is fixed at 75 years, which is obtained by multiplying the nominal structural life VN of 50 years by a coefficient of use Cu equal to 1.5, imposed to structures whose seismic resistance is of importance in view of the consequences associated with their possible collapse, like the case-study school gym building.By referring to topographic category T1 (flat surface) and B-type soil, the resulting peak ground accelerations for the four seismic levels referred to the city of Florence are as follows: 0.065 g (FDE), 0.078 g (SDE), 0.181 g (BDE), and 0.227 g (MCE).The relevant pseudo-acceleration elastic response spectra at linear viscous damping ratio ξ = 5% are plotted in Figure 8. Time-history analyses were developed by assuming artificial ground motions as inputs, generated in families of seven by the SIMQKE-II software [53] from the spectra above.As required by the Italian Standards [47], as well as by several other international seismic Codes and Regulations [36,46], in each time-history analysis the accelerograms were assumed in groups of two simultaneous

Time-History Verification and Performance Assessment Analysis
The performance evaluation enquiry was carried out for the four reference seismic levels fixed in the Italian Standards [47], that is, Frequent Design Earthquake (FDE, with 81% probability of being exceeded over the reference time period V R ); Serviceability Design Earthquake (SDE, with 50%/V R probability); Basic Design Earthquake (BDE, with 10%/V R probability); and Maximum Considered Earthquake (MCE, with 5%/V R probability).The V R period is fixed at 75 years, which is obtained by multiplying the nominal structural life V N of 50 years by a coefficient of use C u equal to 1.5, imposed to structures whose seismic resistance is of importance in view of the consequences associated with their possible collapse, like the case-study school gym building.By referring to topographic category T1 (flat surface) and B-type soil, the resulting peak ground accelerations for the four seismic levels referred to the city of Florence are as follows: 0.065 g (FDE), 0.078 g (SDE), 0.181 g (BDE), and 0.227 g (MCE).The relevant pseudo-acceleration elastic response spectra at linear viscous damping ratio ξ = 5% are plotted in Figure 8.
1.5, imposed to structures whose seismic resistance is of importance in view of the consequences associated with their possible collapse, like the case-study school gym building.By referring to topographic category T1 (flat surface) and B-type soil, the resulting peak ground accelerations for the four seismic levels referred to the city of Florence are as follows: 0.065 g (FDE), 0.078 g (SDE), 0.181 g (BDE), and 0.227 g (MCE).The relevant pseudo-acceleration elastic response spectra at linear viscous damping ratio ξ = 5% are plotted in Figure 8. Time-history analyses were developed by assuming artificial ground motions as inputs, generated in families of seven by the SIMQKE-II software [53] from the spectra above.As required by the Italian Standards [47], as well as by several other international seismic Codes and Regulations [36,46], in each time-history analysis the accelerograms were assumed in groups of two simultaneous horizontal components, with the first one selected from the first generated family of seven motions, and the second one selected from the second family.
The results of the analyses carried out at the FDE and the SDE are evaluated in terms of interlevel drift ratio (i.e., the ratio of the interlevel drift to the interlevel height of the columns) ILDr, which is equivalent to the interstorey drift ratio in the presence of a system of continuous intermediate beams, although without a floor.The maximum ILDr values induced by the most severe among the seven groups of input motions, ILDr,max, are as follows: 0.07% (FDE), 0.09% (SDE) in X, and 0.06% Time-history analyses were developed by assuming artificial ground motions as inputs, generated in families of seven by the SIMQKE-II software [53] from the spectra above.As required by the Italian Standards [47], as well as by several other international seismic Codes and Regulations [36,46], in each time-history analysis the accelerograms were assumed in groups of two simultaneous horizontal components, with the first one selected from the first generated family of seven motions, and the second one selected from the second family.
The results of the analyses carried out at the FDE and the SDE are evaluated in terms of interlevel drift ratio (i.e., the ratio of the interlevel drift to the interlevel height of the columns) ILD r , which is equivalent to the interstorey drift ratio in the presence of a system of continuous intermediate beams, although without a floor.The maximum ILD r values induced by the most severe among the seven groups of input motions, ILD r,max , are as follows: 0.07% (FDE), 0.09% (SDE) in X, and 0.06% (FDE), 0.07% (SDE) in Y, on the first level; 0.13% (FDE), 0.16% (SDE) in X, and 0.53% (FDE), 0.64% (SDE) in Y, on the second level.The drift ratios in X are far below the 0.33% limitation adopted by [47] at the Operational (OP) performance level for frame structures interacting with drift-sensitive non-structural elements, like the masonry infills on the first level and the curtain wall-type windows on the second level, for the main façades of the building, and the infills situated on both levels, for the side façades.The ILD r,max values obtained on the second level in Y are 1.6 times (FDE) and about twice (SDE) the OP-related limit, and also greater than the drift threshold adopted by [47] for the Immediate Occupancy (IO) performance level, equal to 0.5%.
The ILD r,max values computed for the second level in Y are equal to 1.36% at the BDE, and 1.69% at the MCE, assessing moderate (BDE) to high (MCE) potential plastic demands on the columns-should an inelastic finite element analysis be carried out-and severe (BDE) to very severe (MCE) damage of infills and curtain-wall windows.Consequently, the performance level attained in terms of displacement response is Life Safety (LS), both for the BDE and the MCE.At the same time, ILD r,max is no greater than 0.18% (BDE) and 0.22% (MCE) on the first level in Y, i.e., only 13% of the second level values.This identifies a cantilever-like response of the structure along Y, with structural and non-structural damage located on the second level.As discussed in the following Sections, this suggests incorporating the dissipaters on the upper level only, in order to adequately exploit their damping capacity and limit the cost of the retrofit intervention.In X direction, ILD r,max is equal to 0.2% (BDE), 0.25% (MCE) on the first level, and 0.36% (BDE), 0.45% (MCE) on the second level.The interlevel drift profile depicts a frame-like layout along this axis, which approaches a shear-type shape on the second level, as a consequence of the high flexural stiffness of the roof beams in the X-Z vertical plan (which determines nearly a sliding-clamped constraint condition on the top section of the columns).
The BDE-and MCE-related response was assessed also in terms of stress levels.The shear-related checks are met in both directions and for both levels, up to the MCE.On the other hand, the combined axial-force and biaxial-bending-moment stress state checks are met only for the internal columns (C3 through C16, according to the numbering in Figure 4) on the first level at the BDE.The response of the corner columns (C1, C2, C17 and C18) at this level, as well as of all columns on the second level is unsafe starting from the BDE.By way of example, the M X,c -M Y,c biaxial moment interaction curves (being M X,c , M Y,c the bending moments around the X and Y axes) graphed by jointly plotting the two bending moment response histories obtained from the most demanding among the seven groups of MCE-scaled accelerograms, are plotted in Figure 9 for a corner column, namely C17.The boundary of the M X,c -M Y,c elastic interaction domain-which is a function of the time-history variation of the axial force in the columns-is also traced out in the two graphs for the value of the axial force conventionally referred to the basic combination of gravity loads, i.e., N c = 104 kN.
The response curves relevant to the first level highlight maximum M X,c -M Y,c combined values slightly exceeding the safe domain boundary at the BDE, and 1.77 times greater than the corresponding values situated on the boundary, with prevailing contribution of M Y,c , at the MCE.The curves traced out for the second level show more marked unsafe conditions at the BDE, as compared to the first-level ones, and exceed the boundary by a factor equal to 2.07 at the MCE, but with inverted role of the moments (i.e., with prevailing contribution of M X,c for the second level).
at the Operational (OP) performance level for frame structures interacting with drift-sensitive non-structural elements, like the masonry infills on the first level and the curtain wall-type windows on the second level, for the main façades of the building, and the infills situated on both levels, for the side façades.The ILDr,max values obtained on the second level in Y are 1.6 times (FDE) and about twice (SDE) the OP-related limit, and also greater than the drift threshold adopted by [47] for the Immediate Occupancy (IO) performance level, equal to 0.5%.
The ILDr,max values computed for the second level in Y are equal to 1.36% at the BDE, and 1.69% at the MCE, assessing moderate (BDE) to high (MCE) potential plastic demands on the columns-should an inelastic finite element analysis be carried out-and severe (BDE) to very severe (MCE) damage of infills and curtain-wall windows.Consequently, the performance level attained in terms of displacement response is Life Safety (LS), both for the BDE and the MCE.At the same time, ILDr,max is no greater than 0.18% (BDE) and 0.22% (MCE) on the first level in Y, i.e., only 13% of the second level values.This identifies a cantilever-like response of the structure along Y, with structural and non-structural damage located on the second level.As discussed in the following Sections, this suggests incorporating the dissipaters on the upper level only, in order to adequately exploit their damping capacity and limit the cost of the retrofit intervention.In X direction, ILDr,max is equal to 0.2% (BDE), 0.25% (MCE) on the first level, and 0.36% (BDE), 0.45% (MCE) on the second level.The interlevel drift profile depicts a frame-like layout along this axis, which approaches a shear-type shape on the second level, as a consequence of the high flexural stiffness of the roof beams in the X-Z vertical plan (which determines nearly a sliding-clamped constraint condition on the top section of the columns).The BDE-and MCE-related response was assessed also in terms of stress levels.The shear-related checks are met in both directions and for both levels, up to the MCE.On the other hand, the combined axial-force and biaxial-bending-moment stress state checks are met only for the internal columns (C3 through C16, according to the numbering in Figure 4) on the first level at the BDE.The response of the corner columns (C1, C2, C17 and C18) at this level, as well as of all columns on the second level is

Characteristics of the Protective System
Fluid viscous devices are among the most widely used types of rate-dependent passive energy dampers installed in dissipative bracing technologies worldwide.This is owed to their high damping capacities, stable mechanical properties over time, simple installation procedures, limited architectural and visual impact, competitive costs and, in the case of pressurized elements, inherent self-centering qualities [54][55][56][57].
Within this class, a special type of pressurized FV devices has been studied for several years by the author and co-authors, focusing attention on their mechanical characterization, the implementation of analytical and numerical models to simulate their dynamic response, the formulation of sizing and design criteria, and the application to several different protective technologies and structural typologies.Concerning their analytical modelling, the time-dependent F D damping and F ne non-linear elastic reaction forces corresponding to the damper and spring functions are effectively simulated by the following expressions [54,56]: where t = time variable; c = damping coefficient; sgn(•) = signum function; .
x(t) = device velocity; |•| = absolute value; γ = fractional exponent, ranging from 0.1 to 0.2; F 0 = static pre-load force; k 1 , k 2 = stiffness of the response branches situated below and beyond F 0 ; and x(t) = device displacement.For the development of the numerical analyses, the finite element model of a FV spring-damper is obtained by combining in parallel a non-linear dashpot and a non-linear spring with reaction forces given by expressions ( 22) and ( 23), respectively.Both types of elements are currently incorporated in commercial structural analysis programs, such as the SAP2000NL code used in this study.
The installation layout of the spring-dampers and the dissipative bracing system within the frame skeleton is illustrated by the drawings in Figure 10, referred to the case-study building and corresponding to the basic configuration devised for RC buildings.
Therein, a pair of interfaced devices is placed in parallel with the connecting beam axis at the tip of each couple of supporting braces.A half-stroke initial position is imposed on site to the pistons of both spring-dampers, so as to obtain symmetrical tension-compression response cycles-like the one traced out in the scheme of Figure 2-starting from a compressive-only response of the single devices.This position is obtained during the assembly operations by acting on a pair of threaded steel bars crossing the interfacing plate of each device and connected to two other bored plates, screwed into the external casing of the spring-dampers.

Application of the Design Method to the Case-Study Building
As shown in the building plan in Figure 11, the dissipative braces are placed in four alignments parallel to X (named Al.X1 through Al.X4) and four alignments parallel to Y (Al.Y1-Al.Y4).The latter are constituted by pairs of adjacent columns (C1-C3-Al.Y1, C15-C17-Al.Y2, C2-C4-Al.Y3, C16-C18-Al.Y4).Concerning the X-parallel alignments, because the beam span is about 11 m long, four additional RC columns with mutual section 250 mm × 250 mm, named AC1 through AC4 in Figure 11, are built at a distance of 3 m from the corner columns prior to mounting the bracing members.

Application of the Design Method to the Case-Study Building
As shown in the building plan in Figure 11, the dissipative braces are placed in four alignments parallel to X (named Al.X1 through Al.X4) and four alignments parallel to Y (Al.Y1-Al.Y4).The latter are constituted by pairs of adjacent columns (C1-C3-Al.Y1, C15-C17-Al.Y2, C2-C4-Al.Y3, C16-C18-Al.Y4).Concerning the X-parallel alignments, because the beam span is about 11 m long, four additional RC columns with mutual section 250 mm × 250 mm, named AC1 through AC4 in Figure 11, are built at a distance of 3 m from the corner columns prior to mounting the bracing members.

Application of the Design Method to the Case-Study Building
As shown in the building plan in Figure 11, the dissipative braces are placed in four alignments parallel to X (named Al.X1 through Al.X4) and four alignments parallel to Y (Al.Y1-Al.Y4).The latter are constituted by pairs of adjacent columns (C1-C3-Al.Y1, C15-C17-Al.Y2, C2-C4-Al.Y3, C16-C18-Al.Y4).Concerning the X-parallel alignments, because the beam span is about 11 m long, four additional RC columns with mutual section 250 mm × 250 mm, named AC1 through AC4 in Figure 11, are built at a distance of 3 m from the corner columns prior to mounting the bracing members.Based on the results of the assessment analysis reported in Section 4, the dampers are incorporated on the second level only, since the first level drifts, computed in current conditions, can produce only a marginal activation even of the smallest FV devices in standard manufacturing.Consequently, traditional non-dissipative braces are installed on the first level in the same four plus four vertical alignments, so as to provide the necessary structural continuity with the dissipative bracing system placed on the second level, as illustrated in the elevation view of Figure 10, but without adding any further supplemental damping contribution.The application of steps 2 through 4 of the design procedure is summarized below for X and Y directions.Based on the results of the assessment analysis reported in Section 4, the dampers are incorporated on the second level only, since the first level drifts, computed in current conditions, can produce only a marginal activation even of the smallest FV devices in standard manufacturing.Consequently, traditional non-dissipative braces are installed on the first level in the same four plus four vertical alignments, so as to provide the necessary structural continuity with the dissipative bracing system placed on the second level, as illustrated in the elevation view of Figure 10, but without adding any further supplemental damping contribution.The application of steps 2 through 4 of the design procedure is summarized below for X and Y directions.

X Direction-Lack of Bending Moment Strength in the Columns
Step 2. The verification analysis in current conditions highlights that the most critical response parameters in X direction are the bending moments around Y (i.e., M Y ) in the first-level columns, with the highest unsafe conditions checked in the four corner columns.By referring to the nomenclature in Section 2.1, the maximum moment M a Y,c corresponding to the peak response point in Figure 9c and associated with the concurrent axial force N c = 104 kN mentioned above, is equal to 398.7 kNm.The elastic limit moment M e Y,c of the corner columns around the Y axis is equal to 224.8 kNm.Thus, the stress reduction factor α s = α MY of the critical members in X direction results as follows: Passing from the member to the frame-structure level (equivalent to the frame structure storey for this case-study building), since all columns have the same cross section, α F ratio coincides with α MY : Step 3. On the basis of this value, the equivalent viscous damping ratio of the set of spring-dampers to be installed on the second level is evaluated by means of relation (14), obtaining: Step 4. The E D energy dissipation capacity of the spring-dampers is calculated by expression (18).The elastic limit values of the level shear F e (i.e., the sum of the elastic limit shear forces of the columns) and the first interlevel drift ILD e (replacing ID e in this case) computed in X direction, named F e,X and ILD e,1L,X , are equal to 969 kN and 22 mm, respectively.Introducing these values, as well as α F and ξ eq,αF values given by ( 25) and (26), in (18), the following E D estimate is derived: E D,αF = 2πα F F e,X ξ eq,αF ID e,1L,X = 65.6 kJ (27) Dividing E D by the number of spring-dampers placed in X, the minimum energy dissipation capacity E D,X,d to be assigned to each of the eight devices in order to reach the target performance at the MCE results as follows: E D,X,d = 8.2 kJ.The spring-damper type with the nearest nominal energy dissipation capacity, E n , to E D,X,d has the following mechanical properties, drawn from the manufacturer's catalogue [48]: E n = 9 kJ; stroke s max = ±30 mm; damping coefficient c = 9.9 kN(s/mm) γ , with γ = 0.15; F 0 = 17 kN; and k 2 = 1.74 kN/mm.

Y Direction-Lack of Bending Moment Strength in the Columns and Excessive Inter Level Drift
Step 2. The critical response parameters in Y direction are represented by the bending moments around X (M X ) in the second-level columns, with the highest unsafe conditions checked in the four corner columns too, and the second interlevel drifts.The maximum moment M a X,c , corresponding in this case to the peak response point in Figure 9d, is equal to 174.2 kNm, whereas the elastic limit moment M e X,c is equal to 84.2 kNm.Therefore, the stress reduction factor α s = α MX of the critical members in Y direction is: = 2.07 (28) and thus: The deformation-related reduction factor α d given by ( 19) is calculated for the ILD max (replacing ID max ) and ILD e (replacing ID e ) values computed on the second level in Y direction, named ILD max,2L,Y , ILD e,2L,Y in the following.ILD max,2L,Y -corresponding to the ILD r,max value of 1.69% mentioned in Section 4.2-is equal to 72.7 mm, and ILD e,2L,Y to 36.8 mm, yielding: Step 3. The equivalent viscous damping ratio is calculated in this case by referring both to α F and to α d , using expressions ( 14) and (20), respectively: Step 4. Named F e,Y the elastic limit level shear in Y direction, by applying the E D energy dissipation capacity expressions ( 18) and ( 21 where F e,Y = 638 kN.By comparing ξ eq,αF with ξ eq,αd , and E D,αF with E D,αd , it can be observed that the relevant ratios are rather different.Indeed: ξ eq,αd /ξ eq,αF = 1.89,E D,αF /E D,αd = 1.1.This is due to the fact that, consistently with the general ξ eq expression ( 5), the damping coefficient depends on E e , and thus on the elastic properties of the device, which are a function of the maximum displacement and force reached in the time-history response, in addition to the hysteretic response.On the other hand, the dissipated energy E D is only determined by the area covered by the response cycles, which identifies it as a more stable and reliable parameter for the design of the FV devices.ξ eq is only a useful synthetic measure of their limit damping capacity.
The design process is completed by referring to the largest of the energy values, E D,αF , E D,αd , i.e., E D,αF = 100.8kJ.Similarly to the X direction, the minimum energy dissipation capacity of each of the eight devices placed in Y, E D,Y,d , in order to achieve the target performance at the MCE, is obtained by dividing E D,αF by the number of spring-dampers: E D,Y,d = 12.6 kJ.The device with the nearest nominal energy dissipation capacity to E D,Y,d has the following mechanical properties: E n = 14 kJ; stroke s max = ±40 mm; damping coefficient c = 14.16 kN(s/mm) γ , with γ = 0.15; F 0 = 28 kN; and k 2 = 2.1 kN/mm.

Numerical Verification of the Retrofit Solution
A perspective view of the model including the dissipative bracing system is displayed in Figure 12.The modal analysis carried out in retrofitted conditions confirms the sequence of modes computed in current state, with differences on periods and EMMs lower than 10%, as a consequence of the small stiffening effect of this technology.The periods and EMMs of the two first horizontal translational modes along X and Y pass to 0.83 s (Y) and 0.32 s (X), and to 79.5% (Y) and 88.1% (X); the periods and EMMs of the fourth and fifth mode, translational along X and Y, to 0.25 s (Y) and 0.105 s (X), and to 20.3% (Y) and 11.7% (X).The relevant summed modal masses are nearly equal to 100% along both axes, in this case too.The third and sixth mode, rotational around Z, have periods of 0.33 s and 0.037 s, and EMM equal to 84.9% and 13.8%, giving a summed modal mass of 98.7%.The results of the time-history verification analyses in rehabilitated configuration are synthesized in Figures 13-16, all referred to the response induced by the most demanding of the seven groups of input ground motions scaled at the MCE level.The MXc-MYc interaction curves of the first and second level base sections of column C17, plotted in Figure 9c,d above for the original structure, are duplicated in Figure 13 in retrofitted conditions.The two graphs show that the dissipative action of the protective system allows confining the interaction curves within the biaxial moment safe domain, reducing the maximum MYc (Figure 13a) and MXc (Figure 13b) moments nearly by the targeted αMY and αMX factors of 1.77 and 2.07, as given by ( 24) and (29).The results of the time-history verification analyses in rehabilitated configuration are synthesized in Figures 13-16, all referred to the response induced by the most demanding of the seven groups of input ground motions scaled at the MCE level.The M Xc -M Yc interaction curves of the first and second level base sections of column C17, plotted in Figure 9c,d above for the original structure, are duplicated in Figure 13 in retrofitted conditions.The two graphs show that the dissipative action of the protective system allows confining the interaction curves within the biaxial moment safe domain, reducing the maximum M Yc (Figure 13a) and M Xc (Figure 13b) moments nearly by the targeted α MY and α MX factors of 1.77 and 2.07, as given by ( 24) and (29).
The response cycles of the pairs of spring-dampers situated along the opposite vertical alignments Al.X1, Al.X3, and Al.Y2, Al.Y4 are visualized in Figure 14.The cycles exhibit peak displacements equal to 6.3 mm (Al.X1, Al.X3) and 8.8 mm (Al.Y2, Al.Y4), far below the available stroke limits of ±30 mm (in X) and ±40 mm (in Y) mentioned above.Furthermore, the response of the devices placed in Al.X1 and Al.X3 are nearly coincident, and the same occurs for the devices situated in Al.Y2 and Al.Y4, highlighting that torsion effects in plan are virtually null.
The energy time-histories graphed in Figure 15, where E I and E D denote the total input and dissipated energies, and E I,X , E D,X , E I,Y , E D,Y the relevant portions in X and Y, show E D,X , E D,Y values of 66.9 kJ and 101.9 kJ, which differ only by 2% and 1% from the corresponding E D,αF estimates ( 27) and (33), respectively.
Furthermore, practically identical values of the E D /E I ratio are found in X and Y, namely: E D,X /E I,X = 0.855; E D,Y /E I,Y = 0.85, identifying a well-balanced energy dissipation demand in the two directions.
The roof top displacement time-histories illustrated in Figure 16 show a reduction factor on the peak values equal to about 2.3 (X) and 2.2 (Y) when passing from current (CS) to retrofitted (RS) conditions.In the most deformable direction Y, this corresponds to a drop of the second-level drift from 72.7 mm to 33.8 mm, i.e., below the corresponding elastic limit value ILD e,2L,Y = 36.8mm, as targeted in the design.The results of the time-history verification analyses in rehabilitated configuration are synthesized in Figures 13-16, all referred to the response induced by the most demanding of the seven groups of input ground motions scaled at the MCE level.The MXc-MYc interaction curves of the first and second level base sections of column C17, plotted in Figure 9c,d above for the original structure, are duplicated in Figure 13 in retrofitted conditions.The two graphs show that the dissipative action of the protective system allows confining the interaction curves within the biaxial moment safe domain, reducing the maximum MYc (Figure 13a) and MXc (Figure 13b) moments nearly by the targeted αMY and αMX factors of 1.77 and 2.07, as given by ( 24) and (29).The energy time-histories graphed in Figure 15, where EI and ED denote the total input and dissipated energies, and EI,X, ED,X, EI,Y, ED,Y the relevant portions in X and Y, show ED,X, ED,Y values of 66.9 kJ and 101.9 kJ, which differ only by 2% and 1% from the corresponding , estimates ( 27) and (33), respectively.
Furthermore, practically identical values of the ED/EI ratio are found in X and Y, namely: ED,X/EI,X = 0.855; ED,Y/EI,Y = 0.85, identifying a well-balanced energy dissipation demand in the two directions.
The roof top displacement time-histories illustrated in Figure 16 show a reduction factor on the peak values equal to about 2.3 (X) and 2.2 (Y) when passing from current (CS) to retrofitted (RS) conditions.In the most deformable direction Y, this corresponds to a drop of the second-level drift from 72.7 mm to 33.8 mm, i.e., below the corresponding elastic limit value ILDe,2L,Y = 36.8mm, as targeted in the design.

Conclusions
The energy-based design criterion formulated in this study for the seismic retrofit of frame buildings by means of dissipative bracing technologies does not require any preliminary evaluation of the input energy demand on the original structure.At the same time, the most critical response parameters in current conditions-the reduction of which within the boundary of relevant safe domains (in case of lack of strength), or below the limits preventing damage to structural and non-structural elements (in case of excessive lateral displacements)-are evaluated by a conventional elastic finite element analysis.Both aspects of the initializing step of the sizing procedure allow simplifying the design of supplemental damping-based retrofit solutions, which can be useful especially for professional engineers not familiar with seismic energy computation and the development of non-linear time-history analyses.
The criterion was detailed here for relatively stiff structures, i.e., with a fundamental translational vibration period in original conditions below 0.8 s, where the retrofit design objectives can be met by an added damping, with a marginal role of supplemental stiffness.This prompted to select the bracing systems that incorporate FV spring-dampers as protective devices, because they provide a moderate contribution to the lateral stiffness of the retrofitted structures.However, the procedure can be extended with little modifications to dissipative bracing technologies that significantly increase the translational stiffness too, like the systems including metallic dampers, which will be the subject of a further step of the study.
The demonstrative application to the considered case-study structure allowed checking the quick sizing characteristics of the design criterion, even when stress state-related and drift-related deficiencies are both found in the original structure (as occurs in Y direction of the gym building).Furthermore, the values of the equivalent damping coefficient ratio calculated as a function of the reduction factors α F , α d relevant to stress states and drifts, using formulas ( 14) and (20), respectively, resulted to be notably different.On the other hand, a slight difference was found between the corresponding energy dissipation measures, E D,αF and E D,αd .This identified E D as a more stable and reliable design parameter, as compared to ξ eq , consistently with the fact that E D is only determined by the area covered by the response cycles of the dissipaters.This was also confirmed by the fact that the E D values in X and Y computed from the time-history verification analysis were very similar to the E D,αF estimates.
As targeted in the retrofit design, the incorporation of the protective system in the gym building allows reaching an elastic and safe response of all members, as well as constraining the interlevel drifts below the Immediate Occupancy drift limit, up to the MCE, starting from a rather poor performance in original conditions.
step 1-denoted by point A with coordinates (d A = d max , s A = s max ), where d max , s max are the maximum d and s computed values-below point B with coordinates (d B = d e , s B = s e ),

Figure 1 .
Figure 1.Schematic response of the critical structural member (a) and the frame structure storey to which it belongs (b).

2. a V
Lack of shear strength.In this hypothesis, sA is the shear force in the critical member = sB the elastic limit shear of this member, αs-specified in this hypothesis as Vj α -is evaluated as follows: Lack of bending moment strength.

Figure 1 .
Figure 1.Schematic response of the critical structural member (a) and the frame structure storey to which it belongs (b).

Figure 2 .
Figure 2. Schematic response cycle of the set of fluid viscous (FV) spring-dampers installed on a storey, and parameters for evaluating ξeq.

Figure 2 .
Figure 2. Schematic response cycle of the set of fluid viscous (FV) spring-dampers installed on a storey, and parameters for evaluating ξ eq .

Figure 4 .
Figure 4. Structural plan of the building (dimensions in millimeters).

Figure 5 .
Figure 5. Cross section of the building (dimensions in millimeters).

Figure 4 .
Figure 4. Structural plan of the building (dimensions in millimeters).

Figure 5 .
Figure 5. Cross section of the building (dimensions in millimeters).

Figure 4 .
Figure 4. Structural plan of the building (dimensions in millimeters).

Figure 4 .
Figure 4. Structural plan of the building (dimensions in millimeters).

Figure 5 .
Figure 5. Cross section of the building (dimensions in millimeters).Figure 5. Cross section of the building (dimensions in millimeters).

Figure 5 .
Figure 5. Cross section of the building (dimensions in millimeters).Figure 5. Cross section of the building (dimensions in millimeters).

22 Figure 7 .
Figure 7.View and reference coordinate system of the finite element model.

Figure 7 .
Figure 7.View and reference coordinate system of the finite element model.

Figure 9 .
Figure 9. Current state (CS).MX,c-MY,c biaxial moment interaction curves at the base section of column C17 on the first level (a,c) and second level (b,d) obtained from the most demanding BDE-scaled (a,b) and MCE-scaled (c,d) group of accelerograms.

Figure 9 .
Figure 9. Current state (CS).M X,c -M Y,c biaxial moment interaction curves at the base section of column C17 on the first level (a,c) and second level (b,d) obtained from the most demanding BDE-scaled (a,b) and MCE-scaled (c,d) group of accelerograms.
Appl.Sci.2018, 8, x FOR PEER REVIEW 14 of 22 steel bars crossing the interfacing plate of each device and connected to two other bored plates, screwed into the external casing of the spring-dampers.

Figure 10 .
Figure 10.Installation details of the dissipative bracing system in the case-study building.

Figure 10 .
Figure 10.Installation details of the dissipative bracing system in the case-study building.

Figure 11 .
Figure 11.Positions of the dissipative bracing system alignments in plan.

Figure 11 .
Figure 11.Positions of the dissipative bracing system alignments in plan.

22 Figure 12 .
Figure 12.View of the finite element model incorporating the dissipative bracing system.

Figure 13 .
Figure 13.Retrofitted structure (RS).MX,c-MY,c biaxial moment interaction curves at the base section of column C17 on the first level (a) and second level (b) obtained from the most demanding MCE-scaled

Figure 12 .
Figure 12.View of the finite element model incorporating the dissipative bracing system.

Figure 12 .
Figure 12.View of the finite element model incorporating the dissipative bracing system.

Figure 13 .
Figure 13.Retrofitted structure (RS).MX,c-MY,c biaxial moment interaction curves at the base section of column C17 on the first level (a) and second level (b) obtained from the most demanding MCE-scaled group of accelerograms.

Figure 13 .
Figure 13.Retrofitted structure (RS).M X,c -M Y,c biaxial moment interaction curves at the base section of column C17 on the first level (a) and second level (b) obtained from the most demanding MCE-scaled group of accelerograms.

Figure 15 .
Figure 15.Retrofitted structure (RS).Energy time-histories obtained from the most demanding MCE-scaled group of accelerograms.

Figure 16 .
Figure 16.Rooftop displacement time-histories in X and Y direction obtained from the most demanding MCE-scaled group of accelerograms.

Figure 15 .
Figure 15.Retrofitted structure (RS).Energy time-histories obtained from the most demanding MCE-scaled group of accelerograms.

Figure 15 .
Figure 15.Retrofitted structure (RS).Energy time-histories obtained from the most demanding MCE-scaled group of accelerograms.

Figure 16 .
Figure 16.Rooftop displacement time-histories in X and Y direction obtained from the most demanding MCE-scaled group of accelerograms.

Figure 16 .
Figure 16.Rooftop displacement time-histories in X and Y direction obtained from the most demanding MCE-scaled group of accelerograms.