Next Article in Journal
Physically Based Soap Bubble Synthesis for VR
Next Article in Special Issue
Sharing Soil and Building Geophysical Data for Seismic Characterization of Cities Using CLARA WebGIS: A Case Study of Matera (Southern Italy)
Previous Article in Journal
Output-Only Modal Estimation Using Eigensystem Realization Algorithm with Nonstationary Data Correlation
Previous Article in Special Issue
Seismic Performance Evaluation of Steel Buildings with Oil Dampers Using Capacity Spectrum Method
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Influence of Maxwell Stiffness in Damage Control and Analysis of Structures with Added Viscous Dampers

Escuela Técnica Superior de Arquitectura, Universidad Politécnica de Madrid, Av. Juan de Herrera 4., 28040 Madrid, Spain
*
Author to whom correspondence should be addressed.
Appl. Sci. 2021, 11(7), 3089; https://doi.org/10.3390/app11073089
Submission received: 12 March 2021 / Revised: 25 March 2021 / Accepted: 29 March 2021 / Published: 30 March 2021

Abstract

:

Featured Application

In the absence of more accurate data, stiffness of added viscous damping systems based on brace extender properties should be substantially reduced for analysis. Factors between 0.25 and 0.50 are recommended.

Abstract

Viscous damping systems are often implemented in structures to reduce seismic damage. The stiffness of these elements is dominated by the most flexible part of the set including brace extender, auxiliary mounting elements and damping unit. Existing experimental data are used in this study to show that the actual stiffness of the set is about 25% to 50% of the value generally adopted in current engineering practice, which is based solely on the brace extender. A numerical study shows that this reduction has large implications for several variables related to damage control: residual drift ratio, storey acceleration and plastic strain energy dissipated by the frame members. Other variables, such as member forces and rotations, can experience large variations, particularly for non-linear dampers and high damping levels, especially in the top part of the building and more conspicuously for moderate earthquake intensities. In the absence of accurate data, Maxwell stiffness for analysis based on brace extender properties should be substantially reduced, with recommended factors between 0.25 and 0.50. Given the scarcity of experimental data, these results should be considered preliminary.

1. Introduction

During the last years of the past century, several seismic events (i.e., 1994, Northridge, USA; 1995, Kobe, Japan) took place that substantially shifted the focus of earthquake engineering [1,2]. Even though many of the buildings involved were designed according to modern seismic codes, non-structural damage resulted in enormous economic losses and traumatic post-event recovery, making the engineering community realize that enhanced seismic performance, beyond the mere life-safety level, was necessary [3,4]. Recent seismic events (2010, Chile; 2011, Christchurch, New Zealand) have thrust the need for even more demanding objectives, with emphasis on structural reparability, damage-protection and resilience [5,6,7]. In this regard, declassification of the military viscous damper technology (developed during the cold war) during the 1990s allowed application of these devices to vibration and earthquake protection of civil structures [8]. Addition of viscous dampers to structures is beneficial in several ways: first, seismic energy input is dissipated through viscous mechanisms rather than structural hysteresis, thus providing an effective protection to the main structural elements [9,10]; second, interstorey drift, peak storey acceleration and peak velocity are reduced, preventing damage to non-structural elements [11,12]. For these reasons, viscous dampers have been widely used to provide added seismic resilience to structures [13,14,15]. The Japanese case is particularly illustrative: The Japan Society of Seismic Isolation (JSSI) was founded in 1993, and, soon after, the 1995 Kobe earthquake boosted the application of passive protective technologies in this country [16]; by 2011, supplemental damping systems had been already validated with full-scale tests [17] and around 1000 buildings in Japan already counted with passive protection systems [18]; during the 2011 Tohoku-Oki Earthquake, 11 buildings equipped with these types of systems were instrumented with sensors; the recorded signals and post-earthquake inspection showed their excellent seismic behavior [19]. The first codified provisions for supplemental damping systems were developed in the US [20]. Revision of the European seismic code currently under work contemplates inclusion of clauses relative to these systems [21].
The behavior of these devices has been thoroughly studied [22,23,24,25,26,27,28,29]. A fractional derivative model has been proposed as the most accurate mathematical representation. For practical purposes, this model has often been simplified to a pure dashpot. Added viscous damping systems can be modeled by the serial combination of a spring and a dashpot, where the former represents the flexibility of the different elements linking the pure viscous element to the main structure, and the damper flexibility itself; this last term is often neglected. For a very rigid spring, the system can be assimilated to a pure dashpot, and a closed-form solution under harmonic load is possible; this conveniently allows for the use of simplified methods, such as response spectrum analysis (RSA) [30] implemented in ASCE/SEI 7–16 [31]. However, experimental [32,33] and analytical [34] studies have shown that inclusion of damping system flexibility is essential in achieving an accurate approximation to the dynamic structural response of systems equipped with added viscous damping. Moreover, a minimum level of stiffness is required for adequate performance of the damping system, as it has been shown that excessive flexibility leads to both substantial reduction of dissipated energy and appearance of damper forces in phase with structural restoring forces, the latter being detrimental for the main structure. Rules have been proposed to estimate the minimum amount of stiffness required to achieve a reasonable level of system efficiency [35,36]. However, they are usually based on the behavior of Single-Degree-Of-Freedom (SDOF) systems. One of the objectives of this study is to investigate their adequacy for Multi-Degree-Of-Freedom (MDOF) systems.
If the stiffness of the auxiliary components of the bracing system is included in the analysis model, it is common practice to estimate its value based only on the brace extender cross section and length, on the assumption that this approach will render a suitable approximation to the system properties [8]. The main reason for this simplification is the absence of accurate data regarding the stiffness of other auxiliary elements (such as gusset plates, cleats, bolts, clevises, frontal plates, etc.) included in the set. Moreover, the stiffness component of the damper element itself is generally unavailable in catalogs and manufacturer data, implying that this variable can be neglected by the designer. A second objective of this study is to examine the validity of this assumption, and to determine to which extent the stiffness of secondary elements—other than the brace extender and including the damper stiffness component—should be included in analysis models.
In order to achieve the study targets, this study comprises two parts. First, available experimental data are used to assess the range of influence of auxiliary elements on the overall flexibility of the added damping system. Second, a parametric study is carried out in which two frames equipped with different viscous damping systems are modeled with and without consideration of the brace elements and auxiliary system stiffness. The frame is subjected to sets of accelerograms under different intensities, and the peak values for relevant variables at different storeys are compared. The study shows that the errors incurred by considering only the brace extender stiffness (without inclusion of other auxiliary elements) are not negligible. In fact, for some variables and storeys, they are of the same order of magnitude than the errors incurred by totally neglecting the stiffness. Thus, an accurate result can only be obtained if the flexibility of all auxiliary elements (including the damper itself) are considered. Given the scarcity of experimental data, these results should be considered as preliminary. However, they suggest the need for future experimental characterization of complete damper units, including all mounting elements.

2. Materials and Methods

2.1. Mathematical Models

Figure 1 shows several mathematical models used to represent viscous damping systems. The pure dashpot (Figure 1a) is an extreme case of the Maxwell system in which the stiffness is infinite. The Maxwell system (Figure 1b) includes both linear stiffness and viscous damping and renders an accurate description of the damper and its mounting elements. The Kelvin model (Figure 1c) presents the spring and damper in parallel and has been proposed by several authors as a simple approximation to the more complex Maxwell system. The standard model (Figure 1d) includes a Maxwell system in parallel with a linear or non-linear spring and can be used to represent the viscous damper inserted within a structural system. In all cases, hysteretic loops for steady state response under harmonic excitation of amplitude uo and circular frequency ω are presented.
The constitutive equation of the pure dashpot is:
f = c u ˙ a s g n ( u ˙ )
where, f indicates force, u is displacement, c is the damping coefficient, sgn(·) is the sign function and a is the damper exponent, with a = 1 corresponding to linear dampers. For earthquake mitigation, exponents below 1 are common. The constitutive equation for the non-linear Maxwell model is:
u ˙ = f ˙ k d + f 1 / a c 1 / a s g n ( f ) ,
where, kd is the damping system stiffness. If kd takes a very small value, the first term on the right side of Equation (2) prevails. Neglecting the second term and after integration, the equation reduces to that of a pure spring with small stiffness and negligible damping. If kd takes a very large value, the same term vanishes, and the equation reduces to a pure dashpot. These extreme cases show the behavior dependence of the Maxwell system on its stiffness. For linear dampers, Equation (2) can be expressed as
f ˙ = k d u ˙ f λ ;
where, λ (=c/kd), expressed in seconds, is referred to as relaxation time. For use in earthquake analysis, the behavior of the models is generally characterized for forced harmonic excitation on the assumption that, under earthquake excitation with broad frequency content, the structure vibrates mainly with its fundamental period. For forced harmonic excitation of circular frequency ω and amplitude u0, the steady-state solution to Equation (3) is:
f = β s k d u 0 s i n ( ω t ) + β d c u ˙ 0 c o s ( ω t ) ,
β s = ω 2 λ 2 1 + ω 2 λ 2 ,
β d = 1 1 + ω 2 λ 2 .
Equation (4) can be compared to the steady-state solution for the Kelvin system, which adopts the same form with βs = βd = 1. For this reason, the Kelvin system has been proposed as an approximation to the Maxwell system. βs and βd are dependent on excitation frequency and relaxation time and represent the relative magnitudes of spring and damper response, respectively. Figure 2 presents a plot of these coefficients for two different values of relaxation time as a function of the excitation period. For very short periods the predominant behavior is that of a pure spring, whereas for very long periods pure damping dominates the response. In between these two extremes, the behavior is a blend of damper and spring, resulting in the characteristic hysteretic shape of the Maxwell system. It is clear that for MDOF structures, the performance of the added damping system will be different for the fundamental mode and higher modes where spring-like, non-dissipative behavior is to be expected.
The energy dissipated in a cycle of harmonic motion for a pure damper is [26]:
W d = π b c ω a u 0 1 + a ,
where, b is a function of a representing the relationship between dissipated energy in the non-linear damper and a linear damper with the same peak force response and amplitude. b can be found as
b = 2 2 + a π Γ 2 1 + a / 2 Γ 2 + a ;
in this expression Γ(·) represents the gamma function. For any non-linear damper, energy equivalence can be used to propose an equivalent linear damper with damping coefficient:
c 1 , e q = b c ω a 1 u 0 a 1 ,
and the relaxation time λa for the non-linear Maxwell system can be defined as c1,eq/kd. For the standard solid model presented in Figure 1d, the viscous damping ratio ξ can then be expressed using Jacobsen’s equation:
ξ = W d 4 π W s ;
where, Ws = ks·u02/2 is the strain energy at peak displacement. Combining the previous equations, and recalling that ω = 2π/T:
k s k d = π λ a ξ T .
Some authors recommend a maximum value of ks/kd in order to keep the ratio λa/T low, thus ensuring prevalence of the damper behavior. For instance, Lin and Chopra [35] suggested ks/kd ≤ 0.2; with a 30% added damping ratio, this condition is transformed approximately into λ/T ≤ 0.02, and the authors conclude that, if this condition is met, the influence of brace stiffness on response is negligible regardless of the damper exponent. Fu and Kasai [36] suggested a more restrictive limit ks/kd ≤ 0.1; for the same viscous damping ratio, this limit results in λ/T ≤ 0.01. The large difference between both limits indicates a certain ambiguity in this condition. Moreover, Equation (11) applies only to SDOF systems. Despite that, its use has been extended to MDOF systems. For this type of systems, however, it has been proposed that a Maxwell stiffness for all storeys kd,n = ks1 (first storey stiffness) is enough for a reasonable performance of the added damping system [37]. Verification of the validity of these conditions is of interest.
Figure 3 presents hysteresis loops for Maxwell models with different values of kd, for the linear (a = 1) and non-linear (a = 0.5) cases, both subjected to the same forced harmonic excitation and energy-equivalent damping coefficients as per Equation (9). The figure shows the features of behavior with limited stiffness: (i) large reduction of dissipated energy expressed as area enclosed within the loop; (ii) appearance of forces at maximum displacement, i.e., increase in storage stiffness; (iii) reduction of forces at maximum velocity, i.e., decrease in loss stiffness; and (iv) moderate reduction of maximum forces.

2.2. Estimate of Maxwell Stiffness

Viscous dampers are generally mounted using brace extenders and attached to the structure by means of strong gusset plates. Figure 4 shows a typical mounting arrangement for a single diagonal element. The stiffness and damping coefficient related to the damper axis are referred to as kb and cb. A series of auxiliary elements (gusset plates, cleats, bolts, clevises, end plates, etc.) are needed in addition to the brace extender and damper. Elastic behavior is assumed for all elements in the set. The damper element itself presents a certain stiffness kbd, which in typical application is assumed as infinite. The stiffness of other elements in the set is generally uncertain at the design stage, as many of these elements are typically provided by the damper manufacturer. For other elements, the stiffness can only be grossly approximated; for instance, the flexibility of joint elements can be estimated using the component-based approach in Eurocode 3, part 1–8 [38]. For this reason, the overall stiffness is often calculated using the brace extender cross section stiffness over the total centerline length. It is assumed that the combined stiffness of the set does not deviate significantly from this value. The model is suitable for seismic analysis not involving limit states of the device, for which a more sophisticated mathematical model has been developed [39].
The actual total damping system stiffness kb can be calculated as:
1 k b = 1 k b o + 1 k b b + 1 k b d .
Figure 5 shows the resulting value of kb considering only two components (kbb and kbo) showing that, according to Equation (12), very flexible elements in the set tend to dominate the overall stiffness, despite large stiffness values of other components.
Accurate measures of real-scale damping system stiffness are not abundant in the scientific literature. In this study, experimental values reported by Akcelyan et al. [34] (listed in Table 1) are used. These data correspond to real-scale state-of-the-art tests performed in the Tokyo Institute of Technology on damping units composed of damper, brace extender and the corresponding connection elements (clevises and brackets) in both edges. The set did not include gusset plates or frontal plate connections; addition of these elements would result in further reduction of stiffness. In the table, kbb has been calculated as Ab·E/Lb, where E is the young modulus of steel (taken as 200 Gpa), kbo has been found using Equation (10) and kbb* has been estimated as Ab·E/L, thus corresponding to the value typically adopted in practical analysis. The last column shows that the actual stiffness (kb) varies between 13% and 28% of the stiffness kbb* found with this usual assumption, a remarkable deviation. The very moderate values of the damper stiffness (kbd) and auxiliary elements (kbo) are noteworthy, as these two values dominate the total stiffness of the set kb. After testing, the devices were installed in a five-storey real-scale structure subjected to further tests in the E-Defense shaking table. The direction, storey of installation of device, and fundamental period T of the test structure are listed.
These data have been used to find equivalent values for a longer damper unit under the assumption of similar damper forces, which leads to equal damper and mounting elements. Thus, kbd and kbo are unchanged, as is the difference between L and Lb representing the total length of damper and auxiliary elements. Length L is replaced by L′ = 8006 mm, thereby finding Lb′ as L′ − (LLb). kbb, kb and kbb* are recalculated accordingly. For this case, the ratio kb/kbb* takes values between 23% and 40%, that is, a smaller influence than in the previous case but still very remarkable. It can be argued that, because the brace element is longer, the cross-section area should be increased due to buckling and robustness considerations. If that is the case, the ratio kb/kbb* decreases further, thus strengthening the conclusion that kb is actually a fraction of kbb* somewhere between, say, 10% and 50%. These limits are, of course, loosely set, based on existing data and should be taken cautiously.
A relevant question is whether this conclusion can be applied to linear dampers. Unfortunately, to these authors’ knowledge, no real-scale accurate data has been published for systems with this type of devices, probably because current seismic practice favors smaller exponents. In order to obtain a range of magnitude, the values in Table 1 have been used with certain assumptions and simplifications. First, non-linear dampers are replaced by linear dampers resulting in similar damper forces; this allows for a similar choice connection and brace extender design. The design forces are not listed in the reference, so they are estimated indirectly using peak drift velocity values v0 derived from peak drift ratios θ assuming harmonic vibration at the natural period T:
v 0 2 π θ h T ,
where h is the corresponding interstorey height. Then, from Equation (1)
c b , 1 v 0 = c b v 0 a .
The corresponding approximate values are listed in Table 2. A second important consideration is the value of kbd. This value is obtained using an assumption on the relaxation time λ of the linear damper. Constantinou and Symans [23] reported values of λ = 0.006 s for dampers with cb = 15.45 N·s/mm; Reinhorn et al. [24] reported λ = 0.014 s for dampers with cb = 201.4 N·s/mm; Seleemah and Constantinou [25] reported λ = 0.008 s for dampers with cb = 17.7 N·s/mm; these values were found under forced harmonic excitation, for reduced-scaled linear dampers with small capacity, and suggest that larger values of cb might be associated with larger values of λ. However, in the absence of data to confirm this observation, a conservative value of 0.006 s (minimum of the values reported) is adopted here. The values of kbb, kb0 and kbb* presented in Table 1 are unchanged, whereas kb is computed using Equation (12). The ratios kb/kbb* obtained (19–42%) are slightly larger than those for non-linear dampers, but still below 50%. Changing the assumed relaxation time does not fundamentally change the results: a relaxation time 10 times shorter (0.0006 s) results in a 26–57% range; a relaxation time 10 times larger (0.06 s) results in a 6–12% range. In both cases, the main point stands valid: there is an important reduction in Maxwell stiffness due to auxiliary elements and damper contribution. A comparison of results for kb/kbb* between Table 1 and Table 2 hint that similar reductions are to be expected for lower exponents. This point, however, needs experimental confirmation.

2.3. Numerical Study

A numerical study is performed in order to investigate the influence of the reduced stiffness in models including added viscous damping systems. In this study, a case is defined as a moment resisting frame (MRF) equipped with an added damping system (ADS) and subjected to earthquake excitation represented by a ground motion set (GMS) at different intensities. Peak results for significant variables obtained for the same case using different assumptions for the Maxwell system stiffness, are compared to assess the influence of this parameter in the result variability. Hereby, the study is described in detail:

2.3.1. Description of Moment Resisting Frames (MRFs)

A five-storey frame (shown in Figure 6), referred to as 05D, is analyzed in this study. In the y-direction, two of these frames provide lateral strength and stiffness to a building whose floor layout is also presented. Seismic weight for each frame is 2601 kN at a typical storey and 2440 kN at roof. The frames were designed for gravity and wind loads, but taking into account capacity criteria according to Eurocode 8: at every joint the condition ∑Mc,Rd ≥ 1.3∑Mb,Rd was imposed, where ∑Mc,Rd and ∑Mb,Rd represent the addition of column and beam resisting moments concurrent at the joint, respectively. Column sizing was adjusted iteratively to ensure a complete global mechanism in a pushover analysis under a 1st mode lateral force pattern. The dynamic properties are: fundamental period T1 = 2.010 s, modal participation factor Γ1 = 1.364, total weight Wtot = 12,847 kN, first-mode effective weight Weff,1 = 9955 kN. Estimative interstorey yield drift ratios θy and storey stiffnesses ks are found using the simplified procedure suggested by Akiyama [9] (Chapter 12), rendering θyT = [0.699; 0.985; 0.997; 1.098; 1.232]%; ks = [63.3; 30.5; 26.7; 19.1; 14.3] kN/mm. ks is used to calculate the ratio ks/kd, as indicated below.

2.3.2. Description of Ground Motion Sets (GMS)

Two different GMSs have been considered: a single representative accelerogram (referred to as EUK), and a set of accelerograms (referred to as ESB). All motions have been taken from the European Strong Motion Database [40,41]. Only horizontal components are considered.
Set EUK comprises only the X-component of the Kalamata Earthquake (000414XA in the database), a far-field accelerogram recorded at ground type B (360 m/s ≤ vs,30 < 800 m/as defined in Eurocode 8) with peak horizontal acceleration PHA = 0.24 g, peak horizontal velocity PHV = 31.5 cm/s and peak horizontal displacement PHD = 6.57 cm, epicentral distance of 11 km, focal depth of 1 km, Arias Intensity IA = 0.553 m/s and significant duration (as defined by Bommer and Mendis [42]) D5–95% = 5.13 s. The accelerogram is normalized by its PHA. The 5% elastic spectrum of this normalized accelerogram, shown in Figure A1 (Appendix A) approaches reasonably well the Eurocode 8 5% normalized elastic response spectrum throughout a wide range of periods.
Set ESB is formed by 20 non-impulsive motions, with magnitude Mw ≥ 5.5, epicentral distance from 10 km to 50 km, recorded at firm soil (ground type B). The motion selection is listed in Table A1, and the statistical properties of the motions are given in Table A2 (Appendix A). Both horizontal components from every record are consistently included in the set, adding up a total of 40 accelerograms. The motions are individually normalized by their PHA and pre-scaled so that the mean 5% elastic spectrum of the pre-scaled set approaches well the 5% Eurocode 8 elastic response spectrum over a wide range of periods. Spectra of pre-scaled accelerograms are shown in Figure A2 (Appendix A). TH response for set ESB is taken as the median ( x ^ ) of individual responses to each accelerogram in the set, assuming a lognormal distribution; correspondingly, the 16% (x16) and 84% (x84) percentiles are taken as indicative of data dispersion:
x ^ = exp μ ln x i ,
x 16 = x ^ / exp ( σ ln x i ) ,
x 84 = x ^ exp ( σ ln x i ) ,
where, xi is the individual response to the i-th accelerogram, μ(·) represents mean and σ(·) standard deviation. Time history analysis is extended 10 s after the end of the recorded signal, in order to achieve a state of negligible motion in the structural system.

2.3.3. Seismic Intensities

Intensities are defined by the value of spectral acceleration at the first mode period and 5% damping Sa (T1,5%), a common intensity measure hereby referred to as Sa1. Use of 5% damping allows for a direct relationship to code-defined spectral values with no intervention of damping correction factors and is consistent with the Ground Motion Set pre-scaling procedure adopted. The design intensity level for the frames was fixed at Sa1,D = 0.22 g, corresponding to a Eurocode 8–defined 5% elastic spectrum anchored at PHAD = 0.36 g. The intensity level was fixed so that the maximum interstorey drift ratio for GMS ESB was approximately 2%. According to Eurocode 8 [43] (section 4.4.3.2), this condition implies an interstorey drift ratio of 1% (importance class II) to 0.8% (importance classes III or IV) in the damage limitation check, which for these frames happens to be a more severe check than frame strength. This situation is usual for MRFs. The inclusion of ADSs is targeted at performance enhancement, measured as reduction of interstorey drift ratio. The design intensity level is referred to with the letter ‘D’ (from design) and it is assumed to correspond to a return period of 475 years. Besides the design level, two additional intensity levels are defined, namely, one for which the frames with added damping remain totally elastic (Sa1,E, referred to with letter ‘E’, from elastic) and a rare earthquake level (Sa1,C, referred to with letter ‘C’, from near collapse) taken as 1.5·Sa1,D, in accordance with the relationship between design and maximum considered earthquake spectra established in ASCE/SEI 7-16 [31] (section 11.4.7). These two levels are assumed to correspond roughly to return periods of less than 95 years and about 2475 years. The analysis for GMS ESB is conducted by scaling at these three levels. For GMS EUK, the analysis is performed at 20 equally spaced intensity levels; the peak intensity was Sa1,max = 0.55 g, corresponding to the Eurocode 8 5% elastic spectrum anchored at PHA = 1.08 g.

2.3.4. Description of Added Damping Systems (ADS)

Added Damping Systems (ADSs) are included in order to obtain enhanced performance, measured as a reduction of interstorey drift ratio. Different damper exponents, damping ratios, and height-wise distributions of dampers are considered. An ADS in this study is described by a label of the type “Q.2.1,” with the following meaning:
  • The first letter (“Q” in the example) indicates the height-wise distribution (Q, R, S); the “Q” distribution has constant damping coefficients at all storeys; in the “R” and “S” distribution the damping coefficients decrease linearly with height up to a fraction of the first storey damping coefficient, 2/3 for the “R” distribution and 1/3 for the “S” distribution. Only complete vertical distributions are considered.
  • The second number (“2” in the example) indicates the relative amount of damping in the first mode; for linear dampers, “1” indicates 11%, “2” indicates 23% and “3” indicates 34%; for non-linear dampers, the values are drift-dependent, but the damping coefficients keep a similar progression.
  • The third number (“1” in the example) indicates the damper exponent: “1” indicates linear (a = 1) and “5” indicates non-linear (a = 0.5). All dampers in one ADS are either linear or non-linear.
The “Q” distribution is a realistic option, as in some cases only one type of damper is used throughout the structure. The “S” distribution represents a more efficient option, as the damper coefficients are arranged loosely following the frame storey stiffness distribution. Percentages of 11%, 23%, 34% added damping ratios in the first mode were established for ADSs Q.1.1, Q.2.1 and Q.3.1 (linear dampers); values of 20% have been recommended for optimal performance [8,44]. Non-linear damping coefficients were determined to obtain approximately similar damping ratios at Sa1,C, resulting in higher damping ratios at design level Sa1,D (damping ratios for non-linear dampers are intensity-dependent). For ADSs type “R” and “S”, the coefficients are arranged so that added damping ratio in the first mode ξa,1 is the same as in the corresponding ADS type “Q”. Although optimized vertical distributions of damping coefficients have been proposed in the literature [14,27,29], this topic is out of the scope of this paper, as is the distribution of dampers among bays of the frame. Thus, the distributions in the study are chosen as simple as possible within realistic limits, and the position of dampers is limited to the central bay of the frame.
Brace sections are defined as European hot-finished square hollow sections (SHSH) with the following criteria: (i) S355 grade steel (characteristic yield strength fy = 355 MPa) conforming to Eurocode 3 [45] is adopted; (ii) a maximum reduced slenderness of 2 is allowed; (iii) only compact sections class 1 with a minimum wall thickness of 8 mm are considered. Buckling is assessed considering the whole centerline length. The properties of ADSs are listed in Table 3, where kb has been found considering the cross-sectional properties of the brace with the complete centerline length of the damping system. Typical values for kb suggested in literature are 175 to 525 kN/mm [8,39,44]. The values in Table 3 belong in this range and lie above ks,1 (63.3 kN/mm), as suggested by Chen and Chai [37]. In this table, ξa,1 is listed at Sa1,D for non-linear dampers.

2.3.5. Description of Analysis Models

The analysis is performed using OpenSees 3.2.1 [46,47]. The analysis is planar; no torsional effects are included. The frame is modeled with material and geometrical non-linearity. Columns and beams are modeled using force-based beam-column elements, with five integration points. Cross sections are discretized with fibers, with four elements per flange width and four elements per web height, resulting in a total of 12 midpoints. The model includes P-M interaction but no deterioration effects, which is deemed acceptable for the moderate ductility levels reached in the systems studied, due to the beneficial effect of added viscous damping. The material properties are defined using the Giuffré-Menegotto-Pinto hysteretic model with 0.01% isotropic strain hardening ratio. The column panel zone is modeled with rigid links and an elastic spring with equivalent tangent stiffness [48]; the panel is assumed to be reinforced with cover plates to ensure elastic behavior. To include P-Δ effects, seismic vertical loads are applied on a very flexible lean column with high axial stiffness, attached to the main structure at each storey level (minus the frame tributary gravity loads, which are entered directly in the frame). Inherent damping is modeled as 2% Rayleigh damping in the first and third modes. A program-independent convergence check was performed at the end of analysis to ensure that the Energy Balance Error (EBE, as defined by Christopoulos and Filiatrault [14]) fell under 5%, otherwise reducing the analysis time step until this condition was met.

2.3.6. Description of Brace Stiffness Model Options

Viscous dampers are modeled using the OpenSees Viscous Damper material [49], defined for every storey n by damper exponent an, damping coefficient cb,n and Maxwell stiffness kb,n (cb,n and kb,n indicate properties related to the damper axis orientation, as shown in Figure 4). For every frame and damping system, four different models are created:
  • RB, modeled with a very large stiffness (1000 kN/mm).
  • HB, in which the stiffness is calculated using the brace extender cross section and total diagonal centerline length.
  • MB, in which the stiffness is taken as 50% that of the HB model.
  • LB, in which the stiffness is taken as 25% that of the HB model.
In the first case, the Maxwell stiffness is assumed as very large, and behavior of the model is almost purely viscous; the second case corresponds to a conventional design situation in which the Maxwell stiffness is based on dependable design data; the third and fourth cases include in an approximate way the influence of the auxiliary elements and damper stiffness. As shown in previous sections, the reduction in stiffness is expected to fall somewhere between 25% and 50%; therefore, cases MB and LB correspond approximately to a lower and upper bound of the influence of the stiffness of auxiliary elements.

2.3.7. Variables Analyzed

The results for several relevant variables are examined using (as customary in earthquake engineering) their peak unsigned values throughout time history as representative. For energy variables the values at the end of motion are used. Because added damping systems are included as protective systems, emphasis is placed on those variables closely related to damage measure: interstorey drift ratio θ (or drift Δ); residual interstorey drift ratio θr (or residual drift Δr); absolute storey acceleration at; plastic strain energy at beams (Wbeam), columns (Wcolumn) or frames (Wframe); relative input energy EI [10]; energy dissipated by added damping Wdamper. The following variables are also examined, but only summarized results are provided here for brevity: beam peak rotation θb; column peak rotation θc; structural shear Vs; damper shear Vd; total shear Vt; column moment Mc; beam moment Mb; column unsigned peak axial force Ncol,min. Results for strain energy in structural elements have been obtained by numerical integration of generalized internal forces over generalized internal displacements. Other energy variables have also been similarly obtained by numerical integration [14]. Storey shear Vs excludes damper shear Vd. Total shear Vt includes damper shear Vd and storey shear Vs. Due to uncoupling Vt is not the addition of Vd and Vs for individual motions.

3. Results

3.1. Results for One Motion (GMS EUK)

Results for GMS EUK at different intensity levels are presented in Figure 7 (interstorey drift Δ), Figure 8 (residual interstorey drift Δr), Figure 9 (absolute acceleration at) and Figure 10 (energy dissipated by dampers Wdamper). The plots show the ratio between analysis results of models MB to HB (first and second columns in plots) or LB to HB (third and fourth column in plots), at every intensity level. Every storey is treated as an independent datum; for every intensity level the 50% percentile is highlighted and its value at the maximum intensity labeled. The results are binned by amount of damping at design intensity (low or high), damper type (linear or non-linear) and storey (first of roof). In the case of residual drift ratio Δr, the results are restricted to levels for which Δr at HB is not negligible. Concise results for other variables are given in Appendix B, including: storey results in Figure A3 (drift velocity ∇, damper shear Vd, storey shear Vs, total shear Vt); element results in Figure A4 (beam Moment Mb, beam rotation θb, column moment Mc, column rotation θc, column axial force Nc, storey moment Ms); and energy results in Figure A5 (relative energy input EI, plastic strain energy dissipated by frame Wframe, column Wcolumn and beam Wbeam). For Wframe, Wcolumn, Wbeam only intensity levels with a non-negligible value of the variables have been selected. Even so, very low values of plastic strain are obtained for specific storeys in model type HB and, as a result, the ordinates of the ratio LB/HB and MB/HB present large peaks; the scale of the ordinate axis has been set ignoring these values.

3.2. Results for Set of Motions (GMS ESB)

Results of the study for GMS ESB at different intensity levels are presented in Figure 11 (residual interstorey drift ratio θr), Figure 12 (interstorey drift ratio θ), Figure 13 (absolute acceleration at) and Figure 14 (plastic strain energy for frame Wframe, columns Wcolumn, and beams Wbeam). For brevity, these results are presented graphically only for ADS type “S”; results for types “Q” and “R” follow similar trends, summarized numerically in tabular form, as discussed below. At intensity level E (Sa1,E), the structure remains elastic; therefore, the values of residual drift ratio and energy dissipated by plastic strain are negligible and have been omitted. Appendix B.2 presents a summary of results for other additional variables: relative energy input EI (Figure A6), energy dissipated by damping system Wdamper (Figure A7); plots for other relevant variables (storey shear Vs; damper shear Vd; total shear Vt; beam moment Mb; column moment Mc; peak column axial force Nc; beam rotation θb; and column rotation θc) are available online as Supplementary Materials.
Table A3 in Appendix B.2 presents comprehensive values of ratio between median drift values Δ ^ obtained for models MB and LB to those obtained for model HB, at all intensity levels (E, D, C) and storeys, for all damping systems. These data are summarized in Table 4, which lists average values binned by ADS properties for the same variable (Δ), thus clarifying the influence of type of damper (linear or non-linear), amount of damping (1, 2 or 3) and damper coefficient distribution (Q, R, S); an average for storeys 1 and 5 is also given, to assess the vertical variation of the ratio. The process is repeated for each variable of interest and similar summary tables are produced for: residual drift Δr (Table 5), absolute acceleration at (Table 6) and plastic strain energy dissipated by frame Wframe (Table 7). Appendix B.2 contains additional tables for other relevant variables: relative input energy EI (Table A4), energy dissipated by added damping Wdamper (Table A5), and plastic strain energy dissipated by columns Wcolumn (Table A6) or beams Wbeam (Table A7); additional tables for other relevant variables (storey shear Vs; damper shear Vd; total shear Vt; beam moment Mb; column moment Mc; column axial force Nc; beam rotation θb; and column rotation θc) are available online as Supplementary Materials.

4. Discussion

This discussion focuses on relative change of variable results with the type of model. Values cited hereby do not refer to variable results, but to the ratio between results for two different types of model. Thus, “LB” refers to the ratio between results obtained for model type LB and those obtained for model type HB. Likewise, “MB” refers to the ratio between results for model type MB and HB. All values cited are obtained from the median and dispersion values listed on summary tables (Section 3 and Appendix B). To simplify the discussion, ratios listed in those tables under heading “MB/HB,” sub-column “ALL” are referred simply as “MB.” Ratios listed in those tables under heading “LB/HB,” sub-column “ALL,” are referred simply as “LB.” The results are discussed in detail for every relevant variable hereby.

4.1. Drift

Variation in drift with Maxwell stiffness is moderate, on average about 1.08 (type MB) to 1.25 (type LB). Focusing on type LB, the variation in drift is:
  • Larger for non-linear dampers (1.34) than for linear dampers (1.17).
  • Larger for increasing damping (1.11 for low damping to 1.37 for high damping).
  • Almost unaffected by the type of damper distribution (1.27, 1.25, 1.24 for distributions Q, R, S).
  • Much larger for the top half of the building (1.15 for first storey, 1.53 for roof).
The influence is much larger for moderate intensity (1.35 for intensity level E, 1.18 for intensity level C). Dispersion is virtually unaffected by the type of model, except for high damping at low intensity levels with non-linear dampers (reaching up to 1.26).

4.2. Residual Drift

Residual drift is largely affected by the Maxwell stiffness adopted. On average, the ratio to model type HB is about 1.17 (type MB) to 2.01 (type LB). Focusing on type LB, the variation in residual drift is:
  • Uninfluenced by the damper type; the ratio is very similar for linear dampers (2.07) and non-linear dampers (1.94).
  • Larger for increasing damping (1.38 for low damping to 2.50 for high damping).
  • Moderately influenced by the type of damper distribution (1.74; 1.88; 2.40 for distributions Q, R, S).
  • Much larger for the top half of the building (1.56 for first storey, 2.79 for roof).
The influence is much larger for moderate intensity (2.35 for intensity level D vs. 1.66 for intensity level C).
Result dispersion follows similar trends to the median: dispersion is larger for LB than for MB, increases with increasing damping, is moderately affected by the distribution type, and is much larger at the roof than in the first storey.

4.3. Total Acceleration

Total storey acceleration is remarkably affected by the Maxwell model adopted, with an average of variation 1.19 (type MB) to 1.47 (type LB). Focusing on type LB, the variation in total acceleration is:
  • Uninfluenced by the damper type (1.47 for linear dampers, 1.47 for non-linear dampers).
  • Moderately influenced by the amount of damping (1,36 for low damping, 1.50 for high damping).
  • Almost uninfluenced by the damper distribution (1.50 for Q, 1.44 for S).
  • Larger for the roof (1.65) than for the first storey (1.15).
  • Almost uninfluenced by the level of intensity (1.48 for level E vs. 1.47 for level C).
The dispersion follows similar trends, except that it increases steadily with the intensity level (1.67 at level E vs. 1.96 at level C).

4.4. Relative Energy Input

Relative energy input has been pointed out as a very stable quantity [9,50]; the results obtained in this study confirm that this variable is virtually unaffected by the type of model (1.01 for MB, 1.03 for LB). The largest variation for LB obtained is 1.15, in the first storey at intensity E. The change in dispersion is also very moderate, albeit slightly larger for LB (1.22) than for MB (1.17). The dispersion is slightly larger for smaller intensities and non-linear dampers.

4.5. Energy Dissipated by Damping System and Damper Shear

When the average throughout storeys is considered, the variations of Wdamper with Maxwell stiffness present a remarkable stability, with a slight tendency to decrease for type LB; the average values are 1.02 for MB and 0.97 for LB. These results are virtually unaffected by intensity, type of damper or type of distribution. The dispersion results show a moderate dependency on the type of damper (larger for non-linear dampers) and amount of damping (larger for higher damping).
The height-wise distribution of the variation of Wdamper is largely affected by the type of model, as shown in Figure A7 (Appendix B.2). When the Maxwell stiffness is reduced, the energy dissipated in the bottom half of the building decreases whereas the top half experiences an increase in energy dissipation. The average values for the first storey are 0.83 (MB) and 0.64 (LB). The average values for roof are 1.27 (MB) and 1.52 (LB). Similar trends are obtained for dispersion, with larger values for roof than first storey.
The observations for damper shear Vd are very similar to those for Wdamper, with a stable average throughout storeys close to 1, but with a reduction in first storey value (0.93 for MB, 0.83 for LB) and an increment at roof (1.12 for MB 1.23 for LB).

4.6. Plastic Strain Energy Dissipated by Frame, Columns and Beams

Plastic strain energy is highly influenced by the type of Maxwell model. Considering the plastic strain energy in the whole frame Wframe, the ratio increases from 1.71 (MB) to 12.13 (LB). Focusing on type LB, the variation of Wframe is:
  • Largely influenced by the type of damper (8.32 for linear dampers versus 15.95 for non-linear dampers).
  • Largely influenced by the amount of damping (2.7 for low damping to 23.69 for high damping).
  • Largely influenced by the damper distribution (7.33 for type Q, 19.43 for type S).
  • Larger for the roof (25.30) than for the first storey (1.67).
The results are much larger for beam plastic strain Wbeam (MB 1.81, LB 27.24) than for column plastic strain Wcolumn (MB 1.34, LB 4.19). For columns the variation is particularly high in the bottom of the building, due to the development of plastic hinges at the column bases. The dispersion of plastic strain energy follows similar trends.

4.7. Storey Shear and Total Shear

Storey shear Vs is clearly influenced by the Maxwell stiffness. The average variations are 1.14 (MB) and 1.40 (LB). Focusing on type LB, the trends are similar to those pointed out for drift. The variation is:
  • Larger for non-linear dampers (1.50) than for linear dampers (1.31).
  • Larger for increasing damping (1.19 for low damping to 1.57 for high damping).
  • Almost unaffected by the type of damper distribution (1.46, 1.41, 1.34 for distributions Q, R, S).
  • Much larger for the top half of the building (1.12 for first storey, 2.03 for roof).
  • Larger for moderate intensity (1.52 for intensity level E, 1.33 for intensity level C).
Dispersion is larger for model type LB, particularly for high damping, non-linear dampers, distribution type Q, top storey and low intensity levels.
For total shear, the trends are similar, but the values are more moderate; this is consistent with total shear depending on both storey and damper shear.

4.8. Member Internal Forces

Member internal forces (column moment Mc and beam moment Mb) are influenced by the Maxwell stiffness in a similar way. Considering Mc, the average variations are 1.12 (MB) and 1.34 (LB). Considering Mb, 1.09 (MB) and 1.24 (LB). Considering beam moment Mb and models type LB, the variation is:
  • Larger for non-linear dampers (1.32) than for linear dampers (1.16).
  • Larger for increasing damping (1.36 for low damping to 1.09 for high damping).
  • Almost unaffected by the type of damper distribution (1.28, 1.24, 1.20 for distributions Q, R, S).
  • Much larger for the top half of the building (1.09 for first storey, 1.58 for roof).
  • Larger for moderate intensity (1.38 for intensity level E, 1.13 for intensity level C).
Considering column moment Mc and models type LB, the variation is:
  • Larger for non-linear dampers (1.44) than for linear dampers (1.25).
  • Larger for increasing damping (1.16 for low damping to 1.49 for high damping).
  • Almost unaffected by the type of damper distribution (1.37, 1.34, 1.32 for distributions Q, R, S).
  • Much larger for the top half of the building (1.09 for first storey, 1.81 for roof).
  • Larger for moderate intensity (1.47 for intensity level E, 1.25 for intensity level C).
Column axial force shows little dependency on Maxwell stiffness, with ratios which are systematically close to 1, albeit slightly larger for the roof.

4.9. Member Rotation

Column and beam rotations are quite influenced by the Maxwell stiffness. The average variations are 1.14 (MB) and 1.45 (LB). The trends displayed are very similar to those discussed for beam and column moment.

4.10. Summary

The data presented proves that the value of Maxwell stiffness adopted in analysis exerts a remarkable influence on the analysis output. The results obtained for a single motion (GMS EUK) are consistent with those obtained for a set of motions (GMS ESB). Models type MB and LB show considerable variations in most of variable outputs when compared to model HB. In general, the variations are:
  • Quite large for the top part of the frame and moderate or negligible at the bottom part of the frame.
  • Larger for models type LB (25% Maxwell stiffness) than for models type MB (50% Maxwell stiffness).
  • Larger for moderate seismic intensities than for high seismic intensities.
  • Larger for high damping than for low damping.
  • Larger for non-linear dampers than for linear dampers.
To obtain an overview, Table 8 lists a single value for each variable calculated averaging the values for MB and LB at all intensity levels. With the purpose of providing a practical classification of the effects of Maxwell stiffness, ratios below 1.10 are considered as “small” and left unmarked; ratios above 1.30 are considered as “large” and marked with bold characters and red font, and ratios in between are considered as “moderate” and marked with italics and blue font. The variables can then be classified as follows:
  • Variables with small sensitivity to Maxwell stiffness: relative energy input, damper shear, column axial force.
  • Variables with large sensitivity to Maxwell stiffness: residual drift, absolute acceleration, plastic strain energy dissipated by frame, columns and beams.
For all other variables the sensitivity is, in general, moderate, but:
  • For non-linear dampers and/or high damping, the sensitivity of storey shear, column moment, beam and column rotation is large.
  • For the top part of the frame, the sensitivity of all variables except energy input, column axial force and damper shear is large.
  • Likewise, for the bottom part of the frame the sensitivity for all variables except residual drift, plastic strain energy and column rotation is small.
These results indicate that several fundamental variables for damage control (residual drift, absolute acceleration, plastic strain energy dissipated by structural members) are strongly influenced by the actual value of damping system stiffness.

4.11. Minimum Maxwell Stiffness

As indicated above, several authors have proposed an upper limit to the ratio ks/kd as a mean to obtain a value of stiffness that allegedly allows the Maxwell system to be assimilated to a pure viscous damper. Limiting ratios of 0.1–0.2 have been proposed [35,36]. The data of this study are used to assess the validity of this approach. Figure 15 shows the scatter plot of eight relevant variables (drift Δ, residual drift Δr, total acceleration at, energy dissipated by damper Wdamper, storey shear Vs, damper shear Vd, beam moment Mb, column moment Mc) versus ks/kd, using the analysis data for Ground Motion Set ESB at intensity level C; the y-axis shows the ratio of results for models type HB, MB or LB over results for model RB for each variable; limits for ks/kd suggested in the literature have been marked with a solid (ks/kd = 0.1) and dashed (ks/kd = 0.2) lines; the value of kd used to calculate ks/kd is the reduced value used in analysis. The plots show a poor global correlation of the ratio ks/kd with the variability with stiffness for all variables. Many cases with high ratios ks/kd present a variable value close to 1 (indicating an almost pure viscous behavior); conversely, cases with low ratios ks/kd present high values of the variables (indicating large influence of the Maxwell stiffness). The condition ks/kd < 0.1 (marked in solid line) is sufficient but not necessary: when imposed, it bounds the value of the variables, but a majority of cases with small values of the variables are discarded; additionally, if the flexibility of auxiliary elements is accounted for (cases MB, LB) this limit may not be attainable even with very stiff bracing, because the overall stiffness will be dominated by the weaker element. Thus, limiting the ks/kd ratio does not appear as a convenient condition for the practical choice of brace stiffness if the auxiliary elements are considered. From data presented above it is clear that a condition based solely on damper and storey stiffness does not capture the dependency of the variables on the Maxwell stiffness. As an example, considering drift Δ for ADS R.2.5 at design level, the ratio MB/HB is 1.04 and 1.23 for storeys 1 and 5, respectively, with ks/kd 0.44 and 0.17; it is clear that larger ks/kd ratios are acceptable at first storey than at roof. Thus, existing recommendations for damper stiffness based on SDOF behavior cannot be directly applied to MDOF systems. It is also apparent that the rule to determine kd simply as larger than the first storey stiffness (kdks,1) proposed by Chen and Chai [37] does not constitute a valid condition to determine a satisfactory Maxwell stiffness (from the point of view of damage control).

5. Conclusions

The behavior of added damping systems is largely dependent on the flexibility of elements connecting the damper with the main structure; these include brace extenders, damping units and other auxiliary elements such as gusset plates, bolts, pins, clevises, etc. The elements are connected in series, so the overall stiffness is dominated by the most flexible element in the set. Current engineering practice neglects the influence of the auxiliary elements and favors an estimate of the system stiffness based solely on the bracing element. This is partly motivated by the limited amount of available stiffness measurements performed on complete sets.
In this work, previously published experimental data have been used to estimate that the actual stiffness of the set lies between 50% and 25% of the stiffness based purely on the bracing element. The implications of this result are studied through a numerical study on a 5-storey frame, equipped with different added damping systems, and subjected to different seismic intensity levels, in which the added damping system stiffness is consecutively scaled by 1, 0.5 and 0.25, and results for reduced and unreduced stiffness are compared. The sensitivity of the variables is defined according to their variation with damping stiffness as “small” (less than 10%), “large” (more than 30%) and “moderate” (in between). These results show that:
Considering the average throughout all storeys:
  • Relative energy input, damper shear and column axial force show a small sensitivity to the reduction in damping system stiffness.
  • Residual drift, absolute acceleration and plastic strain energy dissipated by frame elements show a large sensitivity to the reduction in damping system stiffness.
  • For non-linear dampers and/or high damping ratios: storey shear, column moment, beam rotation and column rotation show a large sensitivity to reduction in damping system stiffness.
  • All other variables show a moderate sensitivity to reduction in damping system stiffness.
Considering the height-wise distribution of results:
5.
At the bottom part of the frame, only residual drift, plastic strain energy and column rotation are highly sensitive to reduction in damping system stiffness.
6.
At the top part of the frame all variables are highly sensitive to variation in damping system stiffness, except energy input, column axial force and damper shear.
The main goal of the inclusion of added viscous dampers is damage control; the study shows that three key variables related to damage (residual drift, absolute acceleration and plastic strain energy dissipated by main frame elements) are largely dependent on an accurate evaluation of the added damping system stiffness. Moreover, other relevant variables for damage control (drift, beam and column rotation, beam and column moment) also show a moderate to large dependency on this parameter. Thus, the main conclusion of this study is the need of experimental work to characterize the stiffness of added damping systems. In the absence of more accurate information at the design stage, it seems cautious to perform analyses with a reduced Maxwell stiffness based on a moderate fraction (0.25 to 0.50) of the brace extender stiffness.
As most elements in the set are currently provided by the damping unit manufacturer, it would be convenient if tests were performed on damping systems (including all auxiliary elements whose stiffness cannot be reliably estimated) and the stiffness data published as part of the damper technical information. In that way, at design stage, the engineer would choose a complete set (with clearly defined properties except for the contribution of the brace extender) instead of a damper unit (with undefined stiffness properties and unknown auxiliary elements).
The study results also show that current rules of thumb for the estimation of a minimum stiffness compatible with efficient performance of the added damping system need further refinement.
The study is subjected to the following limitations:
  • Analysis was carried out for only two damper exponents (a = 1, a = 0.5). Comparison between trends for linear and non-linear cases in the study suggests an equal or larger sensitivity to Maxwell stiffness for lower damper exponents (a < 0.5). Additional work is needed to validate this assumption.
  • The influence of horizontal distribution of dampers among bays has not been explored.
  • Only simple vertical distributions of damping coefficients have been considered.
The conclusions are based on analyses performed over a flexible 5-storey moment resisting steel frame with different added damping systems subjected to ground motions recorded in firm soil. Specific results for some variables might be different for taller frames, other structural types, soft soils or incomplete vertical damper distributions. Due to scarcity of available experimental data, the results presented must be considered as preliminary.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/app11073089/s1, Figures: Figure S1 (Storey shear Vs for frame 05D with Added Damping System type S, analyzed under Ground Motion Set ESB at intensity level D); Figure S2 (Damper shear Vd for frame 05D with Added Damping System type S, analyzed under Ground Motion Set ESB at intensity level D); Figure S3 (Total shear Vt for frame 05D with Added Damping System type S, analyzed under Ground Motion Set ESB at intensity level D); Figure S4 (Beam moment Mb for frame 05D with Added Damping System type S, analyzed under Ground Motion Set ESB at intensity level D); Figure S5 (Column moment Mc for frame 05D with Added Damping System type S, analyzed under Ground Motion Set ESB at intensity level D); Figure S6 (Column peak axial force Nc for frame 05D with Added Damping System type S, analyzed under Ground Motion Set ESB at intensity level D); Figure S7 (Beam rotation θb for frame 05D with Added Damping System type S, analyzed under Ground Motion Set ESB at intensity level D); Figure S8 (Column rotation θc for frame 05D with Added Damping System type S, analyzed under Ground Motion Set ESB at intensity level D). Tables: Table S1 (Storey shear Vs; summary of ratios between median values obtained from models MB, LB and HB under Ground Motion Set ESB binned by ADS and storey); Table S2 (Damper shear Vd; summary of ratios between median values obtained from models MB, LB and HB under Ground Motion Set ESB binned by ADS and storey); Table S3 (Total shear Vt; summary of ratios between median values obtained from models MB, LB and HB under Ground Motion Set ESB binned by ADS and storey); Table S4 (Beam moment Mb; summary of ratios between median values obtained from models MB, LB and HB under Ground Motion Set ESB binned by ADS and storey); Table S5 (Column moment Mc; summary of ratios between median values obtained from models MB, LB and HB under Ground Motion Set ESB binned by ADS and storey); Table S6 (Column peak axial force Nc; summary of ratios between median values obtained from models MB, LB and HB under Ground Motion Set ESB binned by ADS and storey); Table S7 (Beam rotation θb; summary of ratios between median values obtained from models MB, LB and HB under Ground Motion Set ESB binned by ADS and storey); Table S8 (Column rotation θc; summary of ratios between median values obtained from models MB, LB and HB under Ground Motion Set ESB binned by ADS and storey).

Author Contributions

Conceptualization, J.C., A.B.; methodology, J.C., A.B.; software, J.C.; validation, J.C., A.B.; writing—original draft preparation, J.C.; writing—review and editing, J.C., A.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data available from authors on demand.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Ground Motion Selection Information

Figure A1. Set EUK (Kalamata, 000414XA) 5% elastic response spectrum vs. Eurocode 8 5% elastic response spectrum; (a) pseudo-acceleration; (b) displacement.
Figure A1. Set EUK (Kalamata, 000414XA) 5% elastic response spectrum vs. Eurocode 8 5% elastic response spectrum; (a) pseudo-acceleration; (b) displacement.
Applsci 11 03089 g0a1
Table A1. Selection of accelerograms for Ground Motion Set ESB.
Table A1. Selection of accelerograms for Ground Motion Set ESB.
NameStationDateFaultMwEpic. Dist. (km)PHAX (g)PHAY (g)Prescaleλx (1)Prescaleλy (1)
Friuli (aftershock)Forgaria-Cornio15.09.1976thrust6170.2640.2182.8612.085
Friuli (aftershock)Forgaria-Cornio15.09.1976thrust6170.3460.3362.2301.537
MontenegroBar-SkupstinaOpstine15.04.1979thrust6.9160.3750.3630.9090.720
Montenegro (aftershock)Tivat-Aerodrom24.05.1979thrust6.2210.1660.1332.4812.198
Montenegro (aftershock)Petrovac-Hotel Oliva15.05.1979oblique5.8240.0990.0891.4511.211
Campano LucanoCalitri23.11.1980normal6.9160.1560.1760.7070.768
KalamataKalamata-Prefecture13.09.1986normal5.9100.2150.2970.8190.982
KylliniZakynthos-OTE Building16.10.1988strike slip5.9140.1510.1460.6851.367
ErzincanErzincan-Meteorologij13.03.1992strike slip6.6130.3890.5130.7651.164
TithoreaAigio-OTE Building18.11.1992normal5.9250.0380.0280.8071.932
Umbria MarcheGubbio-Piana26.09.1997normal6380.0910.0970.8380.760
PotenzaBrienza05.05.1990strike slip5.8280.0960.0802.2081.534
AnoLlosiaAthens 2 (Chalandri District)07.09.1999normal6200.1100.1612.1212.411
GrivaEdessa-Prefecture21.12.1990normal6.1360.1010.0960.8991.100
South AegeanHeraklio-Technical University23.05.1994oblique6.1450.0610.0411.2120.814
StrofadesZakynthos-OTE Building18.11.1997oblique6.6380.1310.1161.1521.542
KozaniKastoria-OTE Building13.05.1995normal6.5500.0190.0201.0241.272
AigionPatra-San Dimitrios Church15.06.1995normal6.5430.0840.0930.7130.976
Duzce 1LDEO Station No C1058 BV12.11.1999oblique7.2110.1110.0731.2571.044
FiruzabadFiroozabad20.06.1994strike slip5.9220.2500.2783.3663.481
(1) Prescale factor.
Table A2. Statistical properties of Ground Motion Set ESB.
Table A2. Statistical properties of Ground Motion Set ESB.
Vs30PHAPHVPHDEpicentral DistanceMwArias Intensity IaSignificant Duration D5–95%
m/sgmm/smmkm m/ss
min365.00.01911.53.410.05.800.012.77
max800.00.5131017.7276.050.07.203.0247.61
mean505.40.165168.039.124.96.240.4616.05
median488.00.12493.215.021.56.050.2210.41
std104.40.118200.953.512.30.410.6512.17
Note: properties obtained before scaling.
Figure A2. Ground motion set ESB: 5% elastic response spectrum vs. Eurocode 8 5% elastic response spectrum; (a) pseudo-acceleration; (b) displacement; (c) pseudo-velocity (Note: Eurocode 8 spectrum for pseudo-velocity is undefined).
Figure A2. Ground motion set ESB: 5% elastic response spectrum vs. Eurocode 8 5% elastic response spectrum; (a) pseudo-acceleration; (b) displacement; (c) pseudo-velocity (Note: Eurocode 8 spectrum for pseudo-velocity is undefined).
Applsci 11 03089 g0a2aApplsci 11 03089 g0a2b

Appendix B. Additional Results

Appendix B.1. Additional Results for Set EUK

Figure A3. Summary of results for storeys (drift velocity ∇ damper shear Vd; storey shear Vs; total Scheme 05D) under Ground Motion Set EUK at different intensity levels, binned by damping system properties; (a) ratio between results for model type MB and model type HB; (b) ratio between results for model type LB and model type HB.
Figure A3. Summary of results for storeys (drift velocity ∇ damper shear Vd; storey shear Vs; total Scheme 05D) under Ground Motion Set EUK at different intensity levels, binned by damping system properties; (a) ratio between results for model type MB and model type HB; (b) ratio between results for model type LB and model type HB.
Applsci 11 03089 g0a3
Figure A4. Summary of results for structural elements (beam moment Mb; beam rotation θb; column moment Mc; column rotation θc; column axial force Nc; storey moment Ms) for frame 05D under Ground Motion Set EUK at different intensity levels, binned by damping system properties; (a) ratio between results for model type MB and model type HB; (b) ratio between results for model type LB and model type HB.
Figure A4. Summary of results for structural elements (beam moment Mb; beam rotation θb; column moment Mc; column rotation θc; column axial force Nc; storey moment Ms) for frame 05D under Ground Motion Set EUK at different intensity levels, binned by damping system properties; (a) ratio between results for model type MB and model type HB; (b) ratio between results for model type LB and model type HB.
Applsci 11 03089 g0a4
Figure A5. Summary of results for energy variables (relative input energy EI; plastic strain energy dissipated by frame Wframe, beams Wbeam and columns Wcolumn) for frame 05D under Ground Motion Set EUK at different intensity levels, binned by damping system properties; (a) ratio between results for model type MB and model type HB; (b) ratio between results for model type LB and model type HB.
Figure A5. Summary of results for energy variables (relative input energy EI; plastic strain energy dissipated by frame Wframe, beams Wbeam and columns Wcolumn) for frame 05D under Ground Motion Set EUK at different intensity levels, binned by damping system properties; (a) ratio between results for model type MB and model type HB; (b) ratio between results for model type LB and model type HB.
Applsci 11 03089 g0a5

Appendix B.2. Additional Results for Set ESB

Table A3. Interstorey drift Δ; ratio between median values obtained from models MB and LB to median values obtained from model HB, Δ ^ M B / Δ ^ H B , Δ ^ L B / Δ ^ H B .
Table A3. Interstorey drift Δ; ratio between median values obtained from models MB and LB to median values obtained from model HB, Δ ^ M B / Δ ^ H B , Δ ^ L B / Δ ^ H B .
IntensityE (Sa1,E)D (Sa1,D)C (Sa1,C)
Storey123451234512345
MB/HB05D2.Q.1.11.011.011.011.021.021.011.011.011.021.021.011.011.011.021.02
05D2.Q.2.11.021.021.031.061.121.031.021.031.061.101.031.021.031.061.10
05D2.Q.3.11.041.041.061.121.251.041.041.061.111.230.990.991.021.091.20
05D2.R.1.11.011.011.011.021.021.011.011.011.021.021.011.011.011.021.03
05D2.R.2.11.021.021.031.061.101.021.021.031.061.101.021.021.031.061.09
05D2.R.3.11.041.041.061.111.221.041.041.061.111.201.051.041.061.111.19
05D2.S.1.11.011.011.021.021.041.021.011.021.021.041.011.011.011.021.04
05D2.S.2.11.031.041.041.071.111.041.041.041.061.101.041.031.041.071.10
05D2.S.3.11.041.051.071.111.181.051.051.071.101.171.051.051.071.111.16
05D2.Q.1.51.061.051.071.171.321.021.021.031.051.081.011.011.021.041.06
05D2.Q.2.51.091.061.131.281.541.041.031.051.131.250.960.961.011.081.18
05D2.Q.3.51.111.101.191.321.631.051.041.101.201.421.041.041.071.171.33
05D2.R.1.51.041.041.061.111.201.021.021.031.051.071.011.011.021.041.06
05D2.R.2.51.071.051.111.221.401.041.031.061.131.231.021.031.051.101.17
05D2.R.3.51.121.111.201.301.541.051.041.111.191.351.031.041.081.151.26
05D2.S.1.51.081.071.091.151.241.031.021.041.061.081.021.011.031.041.06
05D2.S.2.51.101.091.121.211.311.051.051.071.121.191.041.041.061.111.16
05D2.S.3.51.111.131.181.261.371.071.071.121.191.271.051.061.091.151.22
LB/HB05D2.Q.1.11.031.031.041.061.101.031.031.031.061.091.041.031.031.061.10
05D2.Q.2.11.091.081.101.201.401.101.071.091.191.361.041.021.051.171.34
05D2.Q.3.11.151.141.161.371.781.171.131.151.331.691.111.071.111.311.66
05D2.R.1.11.031.031.041.061.101.041.031.041.061.091.081.071.071.091.13
05D2.R.2.11.091.081.101.191.351.091.081.101.181.321.091.071.091.191.32
05D2.R.3.11.151.141.171.351.681.171.131.161.321.601.131.091.121.291.57
05D2.S.1.11.051.041.061.081.121.051.041.061.071.121.041.021.041.071.13
05D2.S.2.11.111.111.131.201.321.131.101.131.191.291.131.101.121.191.29
05D2.S.3.11.161.171.211.321.521.171.161.201.311.471.181.161.191.311.46
05D2.Q.1.51.211.191.211.431.821.081.071.071.141.211.051.041.041.101.18
05D2.Q.2.51.341.281.361.762.541.181.111.151.341.681.041.011.091.271.55
05D2.Q.3.51.421.381.481.872.711.291.201.261.532.181.111.051.141.411.96
05D2.R.1.51.151.141.161.291.511.081.061.081.141.191.041.041.051.111.18
05D2.R.2.51.301.251.291.602.101.191.131.171.331.611.121.101.131.281.49
05D2.R.3.51.451.411.491.822.481.281.211.271.501.941.201.141.211.411.75
05D2.S.1.51.261.241.251.381.591.091.071.101.141.211.051.041.061.121.20
05D2.S.2.51.361.321.321.511.771.221.171.191.311.471.121.091.121.241.39
05D2.S.3.51.471.441.461.691.981.331.271.301.481.711.231.181.221.391.60
Figure A6. Relative energy input for frame 05D with Added Damping System Type S, analyzed under Ground Motion Set ESB at intensity levels E (Sa1,E), D (Sa1,D) and C (Sa1,C).
Figure A6. Relative energy input for frame 05D with Added Damping System Type S, analyzed under Ground Motion Set ESB at intensity levels E (Sa1,E), D (Sa1,D) and C (Sa1,C).
Applsci 11 03089 g0a6
Table A4. Relative energy input EI; summary of ratios between median and dispersion values obtained from models MB, LB and HB under Ground Motion Set ESB binned by ADS and storey.
Table A4. Relative energy input EI; summary of ratios between median and dispersion values obtained from models MB, LB and HB under Ground Motion Set ESB binned by ADS and storey.
ADSStorey Median   E ^ I Dispersion   E I 84 E I 16
MB/HBLB/HBMB/HBLB/HB
TypeAmountaEDCALLEDCALLEDCALLEDCALL
ALLALL1ALL1.011.000.991.001.051.031.021.031.131.141.141.141.201.201.191.20
ALLALL0.5ALL1.021.011.011.011.051.041.021.041.221.201.171.201.291.261.191.25
ALL1ALLALL1.000.960.950.991.041.001.001.041.081.101.121.121.131.151.161.17
ALL2ALLALL0.980.980.991.001.001.011.001.021.141.111.071.131.201.171.121.19
ALL3ALLALL0.970.980.970.991.001.000.971.011.191.191.171.211.281.251.171.26
QALLALLALL1.011.000.991.001.041.041.001.031.161.161.131.151.221.191.131.18
RALLALLALL1.010.991.001.001.061.041.041.051.181.201.171.181.231.221.201.22
SALLALLALL1.031.031.011.031.041.031.021.031.181.161.171.171.291.271.241.27
ALLALLALL11.071.061.061.061.151.131.121.131.111.081.071.091.141.111.091.12
ALLALLALL51.081.051.031.051.051.031.001.031.081.081.081.081.161.171.091.14
ALLALLALLALL1.021.011.001.011.051.041.021.031.171.171.161.171.241.231.191.22
Figure A7. Energy dissipated by damping system Wdamper for frame 05D with Added Damping System Type S at levels D (Sa1,D) and C (Sa1,C); absolute values for RB, HB, MB, LB models.
Figure A7. Energy dissipated by damping system Wdamper for frame 05D with Added Damping System Type S at levels D (Sa1,D) and C (Sa1,C); absolute values for RB, HB, MB, LB models.
Applsci 11 03089 g0a7
Table A5. Energy dissipated by added damping Wdamper; summary of ratios between median and dispersion values obtained from models MB, LB and HB under Ground Motion Set ESB binned by ADS and storey.
Table A5. Energy dissipated by added damping Wdamper; summary of ratios between median and dispersion values obtained from models MB, LB and HB under Ground Motion Set ESB binned by ADS and storey.
ADSStorey Median   W ^ d a m p e r Dispersion   W d a m p e r 84 W d a m p e r 16
MB/HBLB/HBMB/HBLB/HB
TypeAmountaEDCALLEDCALLEDCALLEDCALL
ALLALL1ALL1.031.021.011.021.020.990.970.991.021.011.011.011.041.021.011.02
ALLALL0.5ALL1.011.021.021.010.950.960.950.951.171.091.071.111.251.131.081.16
ALL1ALLALL0.970.970.970.990.930.930.930.951.000.980.971.001.020.980.971.01
ALL2ALLALL0.990.990.991.010.960.960.950.981.051.011.001.051.111.041.011.08
ALL3ALLALL0.990.990.991.010.960.950.920.971.111.071.051.101.191.111.051.14
QALLALLALL1.021.021.011.020.990.980.950.971.091.051.021.051.141.061.011.07
RALLALLALL1.021.021.021.020.990.990.980.991.091.051.041.061.141.081.061.09
SALLALLALL1.011.021.021.010.960.970.960.961.091.071.061.071.151.101.061.10
ALLALLALL10.820.840.840.830.630.650.640.640.970.970.960.970.860.870.830.86
ALLALLALL51.291.271.261.271.561.511.501.521.241.171.171.191.501.391.401.43
ALLALLALLALL1.021.021.021.020.980.980.960.971.091.051.041.061.141.081.041.09
Table A6. Plastic strain energy dissipated by columns Wcolumn; summary of ratios between median values obtained from models MB, LB and HB under Ground Motion Set ESB binned by ADS and storey.
Table A6. Plastic strain energy dissipated by columns Wcolumn; summary of ratios between median values obtained from models MB, LB and HB under Ground Motion Set ESB binned by ADS and storey.
Storey Median   W ^ c o l u m n Dispersion   W c o l u m n 84 W c o l u m n 16
MB/HBLB/HBMB/HBLB/HB
TypeAmountaDCALLDCALLDCALLDCALL
ALLALL1ALL1.271.551.414.315.384.851.251.331.293.093.963.53
ALLALL0.5ALL1.421.541.483.926.215.061.281.291.293.004.023.51
ALL1ALLALL1.261.261.262.572.652.611.071.101.091.511.941.72
ALL2ALLALL1.261.391.333.534.453.991.201.221.212.493.042.77
ALL3ALLALL1.391.841.615.859.737.791.411.481.454.846.605.72
QALLALLALL1.291.341.323.283.873.571.161.131.142.212.552.38
RALLALLALL1.271.491.383.535.144.341.221.271.252.623.272.94
SALLALLALL1.471.811.645.538.396.961.421.541.484.316.155.23
ALLALLALL11.171.061.111.881.371.631.111.121.111.631.311.47
ALLALLALL51.351.331.343.674.714.191.251.171.213.213.533.37
ALLALLALLALL1.351.551.454.115.804.961.271.311.293.053.993.52
Table A7. Plastic strain energy dissipated by beams Wbeam; summary of ratios between median values obtained from models MB, LB and HB under Ground Motion Set ESB binned by ADS and storey.
Table A7. Plastic strain energy dissipated by beams Wbeam; summary of ratios between median values obtained from models MB, LB and HB under Ground Motion Set ESB binned by ADS and storey.
Storey Median   W ^ b e a m Dispersion   W b e a m 84 W b e a m 16
MB/HBLB/HBMB/HBLB/HB
TypeAmountaDCALLDCALLDCALLDCALL
ALLALL1ALL1.511.641.586.6610.128.391.321.391.355.046.275.65
ALLALL0.5ALL1.921.951.9318.0418.4618.251.741.641.6910.879.3610.12
ALL1ALLALL1.291.221.253.112.432.771.191.131.161.781.571.67
ALL2ALLALL1.611.681.648.939.379.151.471.471.475.114.814.96
ALL3ALLALL2.082.312.1923.8129.6926.751.781.791.7816.2116.3216.26
QALLALLALL1.511.521.517.0710.048.551.281.251.275.025.945.48
RALLALLALL1.551.671.619.3411.9710.651.461.501.486.246.856.54
SALLALLALL2.082.192.1420.6520.8720.761.841.781.8112.6110.6611.63
ALLALLALL11.481.151.324.921.833.381.241.071.161.961.141.55
ALLALLALL51.562.061.8117.6836.8027.241.511.981.7513.5422.6918.12
ALLALLALLALL1.711.791.7512.3514.2913.321.531.511.527.967.827.89

References

  1. Norton, J.A.; King, A.B.; Bull, D.K.; Chapman, H.E.; McVerry, G.H.; Larkin, T.J.; Spring, K.C. Northridge earthquake reconnaissance report. Bull. N. Z. Soc. Earthq. Eng. 1994, 27, 235–344. [Google Scholar] [CrossRef] [Green Version]
  2. Park, R.; Billings, I.J.; Clifton, G.C.; Cousins, J.; Filiatrault, A.; Jennings, D.N.; Jones, L.C.P.; Perrin, N.D.; Rooney, S.L.; Sinclair, J.; et al. The Hyogo-Ken Nanbu earthquake (the Great Hanshin Earthquake) of 17 January 1995. Bull. N. Z. Soc. Earthq. Eng. 1995, 28, 1–98. [Google Scholar] [CrossRef]
  3. Bertero, V.V. The need for multi-level seismic design criteria. In Proceedings of the 11th World Conference on Earthquake Engineering, Acapulco, Mexico, 23–28 June 1996. [Google Scholar]
  4. Reitherman, R.K. Five major themes in the history of earthquake engineering. In Proceedings of the 15th World Conference in Earthquake Engineering, Lisbon, Portugal, 24–28 September 2012. [Google Scholar]
  5. Buchanan, A.H.; Bull, D.; Dhakal, R.; MacRae, G.; Palermo, A.; Pampanin, S. Base Isolation and Damage-Resistant Technologies for Improved Seismic Performance of Buildings: A Report Written for the Royal Commission of Inquiry into Building Failure Caused by the Canterbury Earthquakes; Research Report 2011-02; Department of Civil and Natural Resources Engineering, University of Canterbury: Christchurch, New Zealand, 2011. [Google Scholar]
  6. Miranda, E.; Mosqueda, G.; Retamales, R.; Pekcan, G. Performance of Nonstructural Components during the 27 February 2010 Chile Earthquake. Earthq. Spectra 2012, 28, 453–471. [Google Scholar] [CrossRef]
  7. Almufti, I.; Willford, M. REDi Rating System: Resilience based Earthquake Design Initiative for the Next Generation of Buildings, version 1.0; ARUP: San Francisco, CA, USA, 2013. [Google Scholar]
  8. Berquist, M.; DePasquale, R.; Winters, C. Fluid Viscous Dampers, General Guidelines for Engineers Including a Brief History; Taylor Devices Inc.: North Tonawanda, NY, USA, 2019. [Google Scholar]
  9. Akiyama, H. Earthquake-Resistant Limit-State Design for Buildings; University of Tokyo Press: Tokyo, Japan, 1985. [Google Scholar]
  10. Uang, C.M.; Bertero, V.V. Use of Energy as a Design Criterion in Earthquake-Resistant Design; Earthquake Engineering Research Center, University of California: Berkeley, CA, USA, 1988. [Google Scholar]
  11. Sindel, Z.; Akbas, R.; Tezcan, S.S. Drift control and damage in tall buildings. Eng. Struct. 1996, 18, 957–966. [Google Scholar] [CrossRef]
  12. Pavlou, E.; Constantinou, M.C. Response of nonstructural components in structures with damping systems. J. Struct. Eng. 2006, 132, 1108–1117. [Google Scholar] [CrossRef]
  13. Soong, T.; Spencer, B. Supplemental energy dissipation: State-of-the-art and state-of-the-practice. Eng. Struct. 2002, 24, 243–259. [Google Scholar] [CrossRef]
  14. Christopoulos, C.; Filiatrault, A. Principles of Passive Supplemental Damping and Seismic Isolation; Iuss Press: Pavia, Italy, 2006. [Google Scholar]
  15. Symans, M.D.; Charney, F.A.; Whittaker, A.S.; Constantinou, M.C.; Kircher, C.A.; Johnson, M.W.; McNamara, R.J. Energy dissipation systems for seismic applications: Current practice and recent developments. J. Struct. Eng. 2008, 134, 3–21. [Google Scholar] [CrossRef] [Green Version]
  16. Nakashima, M.; Chusilp, P. A partial view of Japanese post-Kobe seismic design and construction practices. Earthq. Eng. Eng. Seism. 2003, 4, 3–13. [Google Scholar]
  17. Kasai, K.; Nakai, M.; Nakamura, Y.; Asai, H.; Suzuki, Y.; Ishii, M. Building passive control in Japan. J. Disaster Res. 2009, 4, 261–269. [Google Scholar] [CrossRef]
  18. The Japan Society of Seismic Isolation. Available online: https://www.jssi.or.jp/english/ (accessed on 9 February 2021).
  19. Kasai, K.; Mita, A.; Kitamura, H.; Matsuda, K.; Morgan, T.A.; Taylor, A.W. Performance of seismic protection technologies during the 2011 Tohoku-Oki earthquake. Earthq. Spectra 2013, 29, 265–293. [Google Scholar] [CrossRef]
  20. Federal Emergency Management Agency. FEMA 273: NEHRP Guidelines for the Seismic Rehabilitation of Buildings; FEMA: Washington, DA, USA, 1997.
  21. Labbé, P. Outlines of the revision of the Eurocode 8, part 1-generic clauses. In Proceedings of the 16th European Conference on Earthquake Engineering, Thessaloniki, Greece, 18–21 June 2018. [Google Scholar]
  22. Makris, N.; Constantinou, M.C. Viscous Dampers: Testing, Modeling and Application in Vibration and Seismic Isolation; Technical Report NCEER-90-0028; National Center for Earthquake Engineering Research: Buffalo, NY, USA, 1990. [Google Scholar]
  23. Constantinou, M.C.; Symans, M.D. Experimental and Analytical Investigation of Seismic Response of Structures with Supplemental Fluid Viscous Dampers; National Center for Earthquake Engineering Research: Buffalo, NY, USA, 1992. [Google Scholar]
  24. Reinhorn, A.M.; Li, C.; Constantinou, M.C. Experimental and Analytical Investigation of Seismic Retrofit of Structures with Supplemental Damping: Part 1-Fluid Viscous Damping Devices; National Center for Earthquake Engineering Research: Buffalo, NY, USA, 1995. [Google Scholar]
  25. Seleemah, A.A.; Constantinou, M.C. Investigation of Seismic Response of Buildings with Linear and Nonlinear Fluid Viscous Dampers; National Center for Earthquake Engineering Research: Buffalo, NY, USA, 1997. [Google Scholar]
  26. Tubaldi, E.; Gioiella, L.; Scozzese, F.; Ragni, L.; Dall’Asta, A. A design method for viscous dampers connecting adjacent structures. Front. Built Environ. 2020, 6. [Google Scholar] [CrossRef] [Green Version]
  27. Pollini, N.; Lavan, O.; Amir, O. Minimum-cost optimization of nonlinear fluid viscous dampers and their supporting members for seismic retrofitting. Earthq. Eng. Struct. Dyn. 2017, 46, 1941–1961. [Google Scholar] [CrossRef]
  28. Scozzese, F.; Dall’Asta, A.; Tubaldi, E. Seismic risk sensitivity of structures equipped with anti-seismic devices with uncertain properties. Struct. Saf. 2019, 77, 30–47. [Google Scholar] [CrossRef] [Green Version]
  29. De Domenico, D.; Ricciardi, G.; Takewaki, I. Design strategies of viscous dampers for seismic protection of building structures: A review. Soil Dyn. Earthq. Eng. 2019, 118, 144–165. [Google Scholar] [CrossRef]
  30. Ramirez, O.M.; Constantinou, M.C.; Kircher, C.A.; Whittaker, A.; Johnson, M.; Gomez, J.D.; Chrysostomou, C.Z.l. Development and Evaluation of Simplified Procedures of Analysis and Design for Structures with Passive Energy Dissi-Pation Systems; Technical Report No. MCEER-00–0010, Revision 1; Multidisciplinary Center for Earthquake Engineering Research, University of Buffalo: Buffalo, NY, USA, 2001. [Google Scholar]
  31. American Society of Civil Engineers (ASCE), A.S.O.C. Minimum Design Loads and Associated Criteria for Buildings and Other Structures; American Society of Civil Engineers (ASCE): Reston, VA, USA, 2017. [Google Scholar] [CrossRef]
  32. Kasai, K.; Ito, H.; Ooki, Y.; Hikino, T.; Kajiwara, K.; Motoyui, S.; Ishii, M. Full scale shake table tests of 5-story steel building with various dampers. In Proceedings of the 7th International Conference on Urban Earthquake Engineering (7CUEE) & 5th International Conference on Earthquake Engineering (5ICEE), Tokyo Institute of Technology, Tokyo, Japan, 3–5 March 2010; pp. 11–22. [Google Scholar]
  33. Dong, B.; Sause, R.; Ricles, J.M. Seismic response and performance of a steel MRF building with nonlinear viscous dampers under DBE and MCE. J. Struct. Eng. 2016, 142, 04016023. [Google Scholar] [CrossRef]
  34. Akcelyan, S.; Lignos, D.G.; Hikino, T.; Nakashima, M. Evaluation of simplified and state-of-the-art analysis procedures for steel frame buildings equipped with supplemental damping devices based on e-defense full-scale shake table tests. J. Struct. Eng. 2016, 142, 04016024. [Google Scholar] [CrossRef]
  35. Lin, W.-H.; Chopra, A.K. Earthquake response of elastic single-degree-of-freedom systems with nonlinear viscoelastic dampers. J. Eng. Mech. 2003, 129, 597–606. [Google Scholar] [CrossRef]
  36. Fu, Y.; Kasai, K. Comparative study of frames using viscoelastic and viscous dampers. J. Struct. Eng. 1998, 124, 513–522. [Google Scholar] [CrossRef]
  37. Chen, Y.-T.; Chai, Y.H. Effects of brace stiffness on performance of structures with supplemental Maxwell model-based brace-damper systems. Earthq. Eng. Struct. Dyn. 2010, 40, 75–92. [Google Scholar] [CrossRef]
  38. European Committee for Standardization. EN-1993-1-8: 2005. Eurocode 3: Design of Steel Structures. Part 1–8: Design of Joints; CEN: Brussels, Belgium, 2005. [Google Scholar]
  39. Miyamoto, H.K.; Gilani, A.S.J.; Wada, A.; Ariyaratana, C. Limit states and failure mechanisms of viscous dampers and the implications for large earthquakes. Earthq. Eng. Struct. Dyn. 2010, 39, 1279–1297. [Google Scholar] [CrossRef]
  40. Ambraseys, N.; Smit, P.; Douglas, J.; Margaris, B.; Sigbjörnsson, R.; Olafsson, S.; Suhadolc, P.; Costa, G. Internet site for European strong-motion data. Boll. Geofis. Appl. 2004, 45, 113–129. [Google Scholar]
  41. ESD, European Strong-Motion Database. Available online: http://www.isesd.hi.is/ESD_Local/frameset.htm (accessed on 9 February 2021).
  42. Bommer, J.J.; Mendis, R. Scaling of spectral displacement ordinates with damping ratios. Earthq. Eng. Struct. Dyn. 2005, 34, 145–165. [Google Scholar] [CrossRef]
  43. European Committee for Standardization. En-1998-1: 2004. Eurocode 8: Design of Structures for Earthquake Resistance. Part 1: General Rules, Seismic Actions and Rules for Buildings; CEN: Brussels, Belgium, 2004. [Google Scholar]
  44. Federal Emergency Management Agency. FEMA-1051: 2015 NEHRP Recommended Seismic Provisions: Design Examples; FEMA: Washington, DC, USA, 2015.
  45. European Committee for Standardization. En-1993-1-1: 2005. Eurocode 3: Design of Steel Structures. Part 1-1: General Rules and Rules for Buildings; CEN: Brussels, Belgium, 2005. [Google Scholar]
  46. McKenna, F. Opensees: A framework for earthquake engineering simulation. Comput. Sci. Eng. 2011, 13, 458–466. [Google Scholar] [CrossRef]
  47. Open System for Earthquake Engineering Simulation Framework; Version 2.1.0; University of California: Berkeley, CA, USA, 2009; Available online: http://opensees.berkeley.edu/ (accessed on 20 July 2012).
  48. Gupta, A.; Krawinkler, H. Seismic Demands for Performance Evaluation of Steel Moment Resisting Frame Structures; Report No. 132; John, A., Ed.; Blume Earthquake Engineering Centre, Dept. of Civ. Eng, Stanford University: Standford, CA, USA, 1999. [Google Scholar]
  49. Akcelyan, S.; Lignos, D.G.; Hikino, T. Adaptive numerical method algorithms for nonlinear viscous and bilinear oil damper models subjected to dynamic loading. Soil Dyn. Earthq. Eng. 2018, 113, 488–502. [Google Scholar] [CrossRef] [Green Version]
  50. Housner, G.W. Limit design of structures to resist earthquakes. In Proceedings of the 1st WCEE, Berkeley, CA, USA, 1 June 1956; pp. 5.1–5.13. [Google Scholar]
Figure 1. Mathematical models for viscous dampers and their steady-state hysteresis loops for forced harmonic excitation; (a) pure dashpot; (b) Maxwell model; (c) Kelvin model; (d) Standard solid model.
Figure 1. Mathematical models for viscous dampers and their steady-state hysteresis loops for forced harmonic excitation; (a) pure dashpot; (b) Maxwell model; (c) Kelvin model; (d) Standard solid model.
Applsci 11 03089 g001
Figure 2. Coefficients βs and βd in Equation (4) as a function of harmonic excitation period T: (a) λ = 0.02 s; (b) λ = 0.1 s.
Figure 2. Coefficients βs and βd in Equation (4) as a function of harmonic excitation period T: (a) λ = 0.02 s; (b) λ = 0.1 s.
Applsci 11 03089 g002
Figure 3. Hysteretic force-displacement loops for Maxwell systems subjected to forced harmonic vibration with different values of kd = 10,000 kN/mm, 20 kN/mm, 10 kN/mm, 5 kN/mm; (a) linear case; (b) non-linear case.
Figure 3. Hysteretic force-displacement loops for Maxwell systems subjected to forced harmonic vibration with different values of kd = 10,000 kN/mm, 20 kN/mm, 10 kN/mm, 5 kN/mm; (a) linear case; (b) non-linear case.
Applsci 11 03089 g003
Figure 4. Stiffness of damper and auxiliary mounting elements.
Figure 4. Stiffness of damper and auxiliary mounting elements.
Applsci 11 03089 g004
Figure 5. Overall stiffness kb as a function of component stiffnesses kbb and kbo.
Figure 5. Overall stiffness kb as a function of component stiffnesses kbb and kbo.
Applsci 11 03089 g005
Figure 6. Description of building and frames used in the study: (a) plan layout; (b) definition of frame.
Figure 6. Description of building and frames used in the study: (a) plan layout; (b) definition of frame.
Applsci 11 03089 g006
Figure 7. Interstorey drift Δ for frame 05D under Ground Motion Set EUK at different intensity levels, binned by damping system properties; (a) ratio between results for model type MB and model type HB; (b) ratio between results for model type LB and model type HB.
Figure 7. Interstorey drift Δ for frame 05D under Ground Motion Set EUK at different intensity levels, binned by damping system properties; (a) ratio between results for model type MB and model type HB; (b) ratio between results for model type LB and model type HB.
Applsci 11 03089 g007
Figure 8. Residual interstorey drift Δr for frame 05D under Ground Motion Set EUK at different intensity levels, binned by damping system properties; (a) ratio between results for model type MB and model type HB; (b) ratio between results for model type LB and model type HB.
Figure 8. Residual interstorey drift Δr for frame 05D under Ground Motion Set EUK at different intensity levels, binned by damping system properties; (a) ratio between results for model type MB and model type HB; (b) ratio between results for model type LB and model type HB.
Applsci 11 03089 g008
Figure 9. Absolute acceleration at for frame 05D under Ground Motion Set EUK at different intensities 10. Energy dissipated by dampers Wdamper for frame 05D under Ground Motion Set EUK at different intensity levels, binned by damping system properties; (a) ratio between results for model type MB and model type HB; (b) ratio between results for model type LB and model type HB.
Figure 9. Absolute acceleration at for frame 05D under Ground Motion Set EUK at different intensities 10. Energy dissipated by dampers Wdamper for frame 05D under Ground Motion Set EUK at different intensity levels, binned by damping system properties; (a) ratio between results for model type MB and model type HB; (b) ratio between results for model type LB and model type HB.
Applsci 11 03089 g009
Figure 10. Energy dissipated by dampers Wdamper for frame 05D under Ground Motion Set EUK at different intensity levels, binned by damping system properties; (a) ratio between results for model type MB and model type HB; (b) ratio between results for model type LB and model type HB.
Figure 10. Energy dissipated by dampers Wdamper for frame 05D under Ground Motion Set EUK at different intensity levels, binned by damping system properties; (a) ratio between results for model type MB and model type HB; (b) ratio between results for model type LB and model type HB.
Applsci 11 03089 g010
Figure 11. Residual drift ratio θr for frame 05D with Added Damping System Type S under Ground Motion Set ESB at intensity levels D (Sa1,D) and C (Sa1,C).
Figure 11. Residual drift ratio θr for frame 05D with Added Damping System Type S under Ground Motion Set ESB at intensity levels D (Sa1,D) and C (Sa1,C).
Applsci 11 03089 g011
Figure 12. Interstorey drift ratio θ for frame 05D with Added Damping System type S, under Ground Motion Set ESB at intensity levels E (Sa1,E), D (Sa1,D) and C (Sa1,C); (a) absolute values for RB, HB, MB, LB models; (b) median ( x ^ ) and dispersion (x16x84) values for RB, HB, MB, LB models relative to RB median.
Figure 12. Interstorey drift ratio θ for frame 05D with Added Damping System type S, under Ground Motion Set ESB at intensity levels E (Sa1,E), D (Sa1,D) and C (Sa1,C); (a) absolute values for RB, HB, MB, LB models; (b) median ( x ^ ) and dispersion (x16x84) values for RB, HB, MB, LB models relative to RB median.
Applsci 11 03089 g012
Figure 13. Absolute acceleration at for frame 05D with Added Damping System type S, analyzed under Ground Motion Set ESB at intensity levels E (Sa1,E), D (Sa1,D) and C (Sa1,C); (a) absolute values for RB, HB, MB, LB models; (b) median ( x ^ ) and dispersion (x16x84) values for RB, HB, MB, LB models relative to RB median.
Figure 13. Absolute acceleration at for frame 05D with Added Damping System type S, analyzed under Ground Motion Set ESB at intensity levels E (Sa1,E), D (Sa1,D) and C (Sa1,C); (a) absolute values for RB, HB, MB, LB models; (b) median ( x ^ ) and dispersion (x16x84) values for RB, HB, MB, LB models relative to RB median.
Applsci 11 03089 g013
Figure 14. Plastic strain energy for frame 05D with Added Damping System Type S, analyzed under Ground Motion Set ESB at intensity levels D (Sa1,D) and C (Sa1,C); (a) total frame plastic strain energy Wframe; (b) column plastic strain energy Wcolumn; (c) beam plastic strain energy Wbeam.
Figure 14. Plastic strain energy for frame 05D with Added Damping System Type S, analyzed under Ground Motion Set ESB at intensity levels D (Sa1,D) and C (Sa1,C); (a) total frame plastic strain energy Wframe; (b) column plastic strain energy Wcolumn; (c) beam plastic strain energy Wbeam.
Applsci 11 03089 g014
Figure 15. Scatter plots of ratios for several variables versus the parameter ks/kd, for Ground Motion Set ESB at intensity level C; (a) drift; (b) residual drift; (c) absolute acceleration; (d) energy dissipated by damper; (e) storey shear; (f) damper shear; (g) beam moment; (h) column moment.
Figure 15. Scatter plots of ratios for several variables versus the parameter ks/kd, for Ground Motion Set ESB at intensity level C; (a) drift; (b) residual drift; (c) absolute acceleration; (d) energy dissipated by damper; (e) storey shear; (f) damper shear; (g) beam moment; (h) column moment.
Applsci 11 03089 g015
Table 1. Properties of non-linear viscous damping systems (after Akcelyan et al. [34]) 1.
Table 1. Properties of non-linear viscous damping systems (after Akcelyan et al. [34]) 1.
DirectionTStoreyLLbAb 2cbakbdkbbkbokbkbb*kb/kbb*
s mmmmmm2kN·(s/mm)0.38-kN/mmkN/mmkN/mmkN/mmkN/mm-
X0.5364–3402524299121490.381197511446045313%
2394721048380980.3819379731510442524%
1470628648380980.3819358533210135628%
Y0.5754–3394721048380980.3819379731510442524%
23849154215,3231960.38438198735717979622%
14629232215,3231960.38438132035617166226%
1 Inferred data is shown in italics. 2 Cross-section area of brace extender. For other symbols refer to main text and Figure 4.
Table 2. Properties assumed for linear viscous damping systems.
Table 2. Properties assumed for linear viscous damping systems.
DirectionTStoreyθhv0fbcb,1λkbdkbbkbokbkbb*kb/kbb*
s %mmmm/skNkN·(s/mm)skN/mmkN/mmkN/mmkN/mmkN/mm-
X0.53640.5430001893591.900.0063167511448845319%
30.5630001983651.850.0063087511448745319%
20.6430002257683.410.00656879731516142538%
10.6234852538033.170.00652858533215135642%
Y0.57540.6530002137513.530.00658979731516342538%
30.7430002427893.260.00654479731515942538%
20.79300025916186.260.0061043198735723579629%
10.76348528916885.840.006974132035621866233%
Table 3. Definition of Added Damping Systems (ADSs).
Table 3. Definition of Added Damping Systems (ADSs).
MRFADSξa,1acb (kN·sa/mma), at Storeykb (kN/mm), at Storey
1234512345
05DQ.1.10.1112.332.332.332.332.33167198167167174
05DQ.2.10.2314.674.674.674.674.67222287222187167
05DQ.3.10.3417.007.007.007.007.00237237287222167
05DR.1.10.1112.872.632.392.151.91198198167174162
05DR.2.10.2315.735.264.784.303.82287287222187174
05DR.3.10.3418.607.887.176.455.73282282287222167
05DS.1.10.1113.723.102.481.861.24187187167147112
05DS.2.10.2317.446.204.963.722.48237237222198162
05DS.3.10.34111.169.307.445.583.72367282237222174
05DQ.1.50.150.530.330.330.330.330.3136167167174174
05DQ.2.50.330.560.060.060.060.060.0222287222187198
05DQ.3.50.540.589.189.189.189.189.1237237237287198
05DR.1.50.150.537.134.030.927.824.7167198167174162
05DR.2.50.330.573.467.361.255.149.0287287222187167
05DR.3.50.540.5109.099.990.981.872.7282282237287198
05DS.1.50.150.547.739.831.823.915.9187198167162112
05DS.2.50.330.594.678.863.147.331.5237237222198162
05DS.3.50.540.5140.5117.193.670.246.8367282237222167
Table 4. Interstorey drift Δ; summary of ratios between median and dispersion values obtained from models MB, LB and HB under Ground Motion Set ESB, binned by ADS and storey.
Table 4. Interstorey drift Δ; summary of ratios between median and dispersion values obtained from models MB, LB and HB under Ground Motion Set ESB, binned by ADS and storey.
ADSStorey Median   Δ ^ r Dispersion Δ84 − Δ16
MB/HBLB/HBMB/HBLB/HB
TypeAmountaEDCALLEDCALLEDCALLEDCALL
ALLALL1ALL1.051.051.051.051.181.161.161.171.011.011.011.011.051.041.051.05
ALLALL0.5ALL1.181.101.071.121.521.281.211.341.091.011.001.031.241.051.011.10
ALL1ALLALL1.030.990.991.031.151.051.041.110.970.970.970.990.990.960.991.00
ALL2ALLALL1.081.041.021.071.311.181.141.241.010.970.961.001.080.990.961.03
ALL3ALLALL1.131.081.061.121.451.321.251.371.060.990.981.041.261.081.041.15
QALLALLALL1.131.081.051.091.391.231.171.271.051.000.971.011.161.040.961.05
RALLALLALL1.111.071.061.081.331.221.191.251.041.001.011.021.141.041.061.08
SALLALLALL1.111.081.061.081.321.221.181.241.051.021.031.031.141.061.081.09
ALLALLALL11.061.031.021.041.211.151.101.151.061.051.021.041.131.111.031.09
ALLALLALL51.261.161.131.181.711.461.411.531.050.940.960.991.190.961.001.05
ALLALLALLALL1.121.071.061.081.351.221.181.251.051.011.001.021.151.041.031.07
Table 5. Residual interstorey drift Δr; summary of ratios between median and dispersion values obtained from models MB, LB and HB under Ground Motion Set ESB binned by ADS and storey.
Table 5. Residual interstorey drift Δr; summary of ratios between median and dispersion values obtained from models MB, LB and HB under Ground Motion Set ESB binned by ADS and storey.
Storey Median   Δ ^ r Dispersion   Δ r 84 Δ r 16
MB/HBLB/HBMB/HBLB/HB
TypeAmountaDCALLDCALLDCALLDCALL
ALLALL1ALL1.271.181.232.411.732.071.161.121.141.791.391.59
ALLALL0.5ALL1.141.101.122.291.591.941.161.101.131.621.241.43
ALL1ALLALL1.041.011.031.551.201.381.041.051.051.151.031.09
ALL2ALLALL1.161.141.152.301.601.951.161.081.121.521.151.34
ALL3ALLALL1.301.161.232.982.012.501.171.101.142.281.631.96
QALLALLALL1.161.091.122.021.461.741.121.041.081.561.171.36
RALLALLALL1.121.101.112.111.641.881.181.111.151.681.331.51
SALLALLALL1.331.221.282.931.882.401.181.191.191.871.441.66
ALLALLALL11.161.051.101.841.281.561.111.041.081.451.121.28
ALLALLALL51.291.281.283.202.392.791.201.161.182.021.701.86
ALLALLALLALL1.211.141.172.351.662.011.161.111.141.711.311.51
Table 6. Absolute acceleration at; summary of ratios between median and dispersion values obtained from models MB, LB and HB under Ground Motion Set ESB binned by ADS and storey.
Table 6. Absolute acceleration at; summary of ratios between median and dispersion values obtained from models MB, LB and HB under Ground Motion Set ESB binned by ADS and storey.
ADSStorey Median   a ^ t Dispersion   a t 84 a t 16
MB/HBLB/HBMB/HBLB/HB
TypeAmountaEDCALLEDCALLEDCALLEDCALL
ALLALL1ALL1.191.191.201.191.471.471.471.471.221.331.381.311.561.821.861.75
ALLALL0.5ALL1.211.181.191.191.491.451.481.471.371.441.461.421.782.052.061.96
ALL1ALLALL1.131.091.091.131.361.311.311.361.171.181.181.201.441.591.551.56
ALL2ALLALL1.181.171.181.201.471.471.491.511.281.391.441.401.692.042.101.99
ALL3ALLALL1.181.181.191.211.461.461.481.501.311.451.501.451.721.982.041.95
QALLALLALL1.211.191.201.201.491.481.511.501.301.411.441.381.681.992.051.91
RALLALLALL1.201.191.191.191.481.471.481.481.291.401.441.381.691.992.041.91
SALLALLALL1.201.181.191.191.461.431.431.441.281.351.381.341.641.821.791.75
ALLALLALL11.051.051.051.051.151.121.101.121.071.081.081.071.201.191.141.18
ALLALLALL51.281.251.251.261.691.631.631.651.181.231.331.251.481.822.001.77
ALLALLALLALL1.201.191.191.191.481.461.471.471.291.381.421.371.671.931.961.85
Table 7. Plastic strain energy dissipated by frame Wframe; summary of ratios between median and dispersion values obtained from models MB, LB and HB under GMS ESB binned by ADS and storey.
Table 7. Plastic strain energy dissipated by frame Wframe; summary of ratios between median and dispersion values obtained from models MB, LB and HB under GMS ESB binned by ADS and storey.
Storey Median   W ^ f r a m e Dispersion   W f r a m e 84 W f r a m e 16
MB/HBLB/HBMB/HBLB/HB
TypeAmountaDCALLDCALLDCALLDCALL
ALLALL1ALL1.451.591.526.5710.068.321.331.361.344.905.275.09
ALLALL0.5ALL1.851.941.9015.1316.7615.951.631.471.559.197.828.51
ALL1ALLALL1.301.211.253.012.402.701.141.121.131.741.541.64
ALL2ALLALL1.641.671.668.569.098.831.361.351.354.694.264.47
ALL3ALLALL1.862.252.0619.9327.4523.691.781.641.7114.0313.2013.61
QALLALLALL1.421.451.446.188.487.331.291.251.274.395.454.92
RALLALLALL1.531.641.587.8511.419.631.371.411.395.575.815.69
SALLALLALL2.012.222.1118.5220.3419.431.761.591.6711.188.389.78
ALLALLALL11.181.071.131.961.391.671.121.111.111.691.311.50
ALLALLALL51.612.041.8216.8133.8025.301.481.801.6411.7717.9514.86
ALLALLALLALL1.651.771.7110.8513.4112.131.481.411.457.056.546.80
Table 8. Summary of variable sensitivity to Maxwell stiffness.
Table 8. Summary of variable sensitivity to Maxwell stiffness.
TypeAmountaStoreyΔΔratEIWdamperWframeWcolumnWbeamVdVsVtMbMcNcθbθc
ALLALL1ALL1.171.591.331.021.006.923.207.541.001.271.191.161.231.021.261.29
ALLALL0.5ALL1.231.531.331.030.988.923.2710.091.011.341.191.221.301.021.341.37
ALL1ALLALL1.071.201.241.010.971.981.932.010.981.121.131.061.111.001.111.14
ALL2ALLALL1.151.551.361.010.995.242.665.400.991.261.201.151.221.011.251.28
ALL3ALLALL1.241.861.351.000.9912.874.7014.470.991.381.211.241.321.011.381.42
QALLALLALL1.181.431.351.010.994.382.455.031.001.301.201.191.241.021.261.29
RALLALLALL1.161.491.341.021.005.612.866.131.001.271.191.161.231.021.251.29
SALLALLALL1.161.841.311.030.9910.774.3011.451.001.231.191.141.221.021.271.30
ALLALLALL11.101.331.091.100.741.401.372.350.881.081.031.051.060.981.121.15
ALLALLALL51.362.041.461.041.4013.562.7714.531.171.691.461.391.541.081.481.55
ALLALLALLALL1.171.591.331.021.006.923.207.541.001.271.191.161.231.021.261.29
Note: Values above or equal 1.30 marked in red bold font. Values above or equal 1.10 and below 1.30 marked in blue italic font.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Conde, J.; Bernabeu, A. Influence of Maxwell Stiffness in Damage Control and Analysis of Structures with Added Viscous Dampers. Appl. Sci. 2021, 11, 3089. https://doi.org/10.3390/app11073089

AMA Style

Conde J, Bernabeu A. Influence of Maxwell Stiffness in Damage Control and Analysis of Structures with Added Viscous Dampers. Applied Sciences. 2021; 11(7):3089. https://doi.org/10.3390/app11073089

Chicago/Turabian Style

Conde, Jorge, and Alejandro Bernabeu. 2021. "Influence of Maxwell Stiffness in Damage Control and Analysis of Structures with Added Viscous Dampers" Applied Sciences 11, no. 7: 3089. https://doi.org/10.3390/app11073089

APA Style

Conde, J., & Bernabeu, A. (2021). Influence of Maxwell Stiffness in Damage Control and Analysis of Structures with Added Viscous Dampers. Applied Sciences, 11(7), 3089. https://doi.org/10.3390/app11073089

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop