# Estimation of Additional Equivalent Damping Ratio of the Damped Structure Based on Energy Dissipation

Authors: Xiuyan Hu; Qingjun Chen; Dagen Weng; Ruifu Zhang; Xiaosong Ren
Journal: Advances in Civil Engineering
 (2019)

---

## Abstract

In the design of damped structures, the additional equivalent damping ratio (EDR) is an important factor in the evaluation of the energy dissipation effect. However, previous additional EDR estimation methods are complicated and not easy to be applied in practical engineering. Therefore, in this study, a method based on energy dissipation is developed to simplify the estimation of the additional EDR. First, an energy governing equation is established to calculate the structural energy dissipation. By means of dynamic analysis, the ratio of the energy consumed by dampers to that consumed by structural inherent damping is obtained under external excitation. Because the energy dissipation capacity of the installed dampers is reflected by the additional EDR, the abovementioned ratio can be used to estimate the additional EDR of the damped structure. Energy dissipation varies with time, which indicates that the ratio is related to the duration of ground motion. Hence, the energy dissipation during the most intensive period in the entire seismic motion duration is used to calculate the additional EDR. Accordingly, the procedure of the proposed method is presented. The feasibility of this method is verified by using a single-degree-of-freedom system. Then, a benchmark structure with dampers is adopted to illustrate the usefulness of this method in practical engineering applications. In conclusion, the proposed method is not only explicit in the theoretical concept and convenient in application but also reflects the time-varying characteristic of additional EDR, which possesses the value in practical engineering.


---

## Body

## 1. Introduction

Energy dissipation technology, which has been used in civil engineering for decades, has proven to be an effective tool to resist earthquakes [1, 2]. Through the balanced installation of dampers in a structure, a large amount of seismic energy can be consumed, thereby protecting the structure from the effects of severe earthquakes [3, 4]. As a commonly used form of seismic technology, energy dissipation devices have been widely studied and extensively developed. Various types of dampers based on different energy-consuming mechanics have been proposed and used in engineering practice worldwide [4–8].In many research studies and in engineering practice, it has been found that the mitigation effects of the seismic response resulting from the installation of dampers can be evaluated according to the additional equivalent damping ratio (EDR) [9, 10]. Moreover, researches show that the energy dissipation capacity of dampers can be reflected by the additional EDR [11, 12]. The additional EDR can connect the mitigation effects of the seismic responses to the energy dissipation capacity of dampers, which may simplify the procedure of retrofit design of structures. In this manner, the common methods used for designing damped structures are related to the required damping ratio [13–17]. Based on the required damping ratio, the parameters of the dampers are determined. On the basis of the abovementioned reasons, it is obvious that the additional EDR obtained from the installed dampers is an important factor in the design procedure of damped structures. The EDR is a basic concept in structural dynamics, and there are some classical methods to estimate the EDR of the structure, such as the free decay oscillation method and the half-power bandwidth method. The half-power bandwidth method is verified to be inaccurate to estimate the EDR of the structure in some cases, which needs some corrections [18, 19]. The most common method is the strain energy method, which is expressed as follows [20]:(1)ξeq=14π⋅EDESo,where ξeq is the EDR, ED is the energy dissipated by the system in one cycle, and ESo is the maximum strain energy in one cycle. The strain energy method is employed in many design codes [21, 22]. However, the method is accurate only when the system is in resonance. In addition, the method has some deficiencies in practical applications, such as complexity in calculating energy dissipation.Because a large variety of dampers are available, the estimation of the additional EDR is not easy. Nevertheless, because of the importance of the additional EDR, many studies have been dedicated to the estimation of this parameter. Lee et al. [23] utilized the modal energy-formed Lyapunov function to obtain a closed-form solution for the additional EDR. Two-dimensional shear buildings with various kinds of dampers were introduced to verify the accuracy of the solution. Occhiuzzi [24] proposed a method based on the state-space representation of a dynamical system to calculate the modal damping ratio. This method was adopted to estimate the additional EDR of structures with dampers that were configured by some existing design methods. The results showed that these damped structures yielded a similar additional EDR, which was approximately 0.2. Diotallevi et al. [25] defined a damping index and presented a method to directly estimate the EDR using the index. Iterations were not needed in the method, which was verified using a single-degree-of-freedom (SDOF) system. Park [26] investigated the damping of an inelastic SDOF system with an added damper. Damping correction factors and regression equations were utilized to estimate the additional EDR. Guerrero et al. [27] proposed an experimental approach to measure the damping of the structures with buckling-restrained braces (BRBs). The approach was verified by the shaking table tests, and the results showed that the energy dissipation of BRBs must be considered in the seismic design of the structure. Papagiannopoulos and Beskos [28] developed a modal damping identification model to estimate the modal damping ratios of the structures. This method was utilized by Katsimpini et al. [29] to obtain the modal damping ratios of the structure with the seesaw system, and the results indicated that the seesaw system could provide significant damping capacity. Li et al. [30] determined the damping in an SDOF system with a Maxwell damper. By adopting the properties of the structure in free vibration, they obtained derivations for the first- and second-order equivalent damping.The abovementioned researches indicate that the additional EDR can be accurately calculated by many methods, but these methods are generally complicated and not convenient for practical application. For practical application, the easiness of the method in understanding and application is as important as the accuracy. Therefore, it is necessary to put forward an approach that is easy to understand, convenient to apply, and accurate to estimate the additional EDR of damped structures.In this study, based on the energy dissipation through structural inherent damping and through the use of dampers, a simple method is developed to estimate the additional EDR induced by the dampers. Compared to the strain energy method, the proposed method is more accurate and the concept is explicit and intuitive, which is also convenient for practitioners to apply. First, the energy governing equation of the system is established. According to the definition of the EDR in structural dynamics [20], the relationship between energy dissipation and the EDR is established, which indicates that the additional EDR reflects the energy dissipation capacity of dampers. In this manner, the ratio of the energy dissipated by structural inherent damping to that dissipated by dampers is utilized to estimate the additional EDR. After that, the determination of the calculation time for the energy dissipation is discussed. The result of energy dissipation can be obtained by the structural energy governing equation. Within the range of significant duration of the ground motion, energy dissipation can be calculated by the integrals of the damping forces of the primary structure and the dampers with respect to the corresponding displacements, respectively, of which the calculation can be completed programmatically or by the structural analysis program. Next, the procedure of the proposed method is presented. In addition, the feasibility of this method is verified by using an SDOF system with an added viscous damper. A six-story concrete benchmark model with metallic yielding dampers (MYDs) is introduced to illustrate the usefulness of this method in practical engineering applications. Finally, the accuracy of results obtained by the proposed method is shown by conducting a comparison study.


## 2. Method Based on Energy Dissipation

### 2.1. Establishment of Energy Governing Equation

The simplified model of an SDOF system with an added viscous damper is shown in Figure1, where the viscous damper is simulated by the Maxwell model. According to structural dynamics, the differential equation of movement subjected to a horizontal earthquake can be written as follows [31]:(2)mx¨+cx˙+fsx,x˙+fdx,x˙=−mx¨g,where x, x˙, and x¨ are the relative displacement, velocity, and acceleration of the mass, respectively; x¨g is the ground acceleration; m is the mass; k is the stiffness of the system; c is the inherent damping coefficient of the system; fsx,x˙ is the elastic or inelastic resisting force of the system; and fdx,x˙ is the damping force of the damper. In the Maxwell model, ka is the dynamic stiffness of the damper and ca is the damping coefficient of the damper.Figure 1
SDOF system with the damper.By integratingx in equation (2) in the range of the earthquake duration [0, t], the energy response equation can be obtained:(3)m∫0tx¨dx+c∫0tx˙dx+∫0tfsdx+∫0tfddx=−m∫0tx¨gdx.Equation (3) can be rewritten as equation (4) in terms of energy:(4)Ek+Ec+Eh+Ee+Ed=Ein,where Ek is the kinetic energy, Ec is the energy dissipated by the inherent damping of the primary structure, and Eh and Ee are the hysteretic energy and elastic strain energy of the primary structure, respectively. Here, Eh+Ee=∫0tfsx˙dt; when the primary structure is in an elastic state, Eh=0. Ed is the energy dissipated by the dampers, and Ein is the energy input by ground motion.


### 2.2. Calculation of Additional EDR

In general, the energy dissipation capacity of the primary structure is indicated in terms of the EDR. Similarly, the energy dissipation capacity of additional devices can be indicated in terms of the additional EDRξeq. According to the classic definition method for the additional EDR, the energy dissipated by the dampers is equalized with that dissipated by an equivalent viscous system during one vibration cycle [20]:(5)Ed=∫02π/ωfdx˙dt=πx02ωceq=2πξeqkx02ωωn,where ω is the circular frequency of excitation, x0 is the maximum displacement of the structure in one cycle, ceq is the equivalent damping coefficient of the dampers in one cycle (ceq=2kξeq/ωn), and ωn is the natural circular frequency of the structure.This definition can be used to rewrite the equation for the additional EDR provided by the dampers during0,t as follows:(6)Edt=∫0tfdx˙dt=∫0tceqx˙2dt=2ξ¯akωn∫0tx˙2dt,where ξ¯a denotes the additional EDR provided by the dampers in the seismic duration 0,t and t denotes the duration of the seismic record.Meanwhile, the energy dissipated by the inherent damping of the primary structure can be expressed as follows:(7)Ect=∫0tcx˙2dt=2ξ0kωn∫0tx˙2dt,where ξ0 denotes the inherent damping ratio of the primary structure.From equations (6) and (7), the additional EDR of the dampers in the duration 0,t can be obtained:(8)ξ¯at=ξ0EdEc=ξ0∫0tfdx˙dt∫0tcx˙2dt.It is obvious that the additional EDRξ¯a can be obtained from the inherent damping ratio and the ratio of the energy dissipated by the dampers to that dissipated by the inherent damping during 0,t. It should be noted that the method focuses on the amount of energy dissipated and the manner to dissipate energy is not involved, which indicates that the method can be applied to various kinds of dampers. Equation (8) provides the expression for calculating the energy dissipation. When the damping force of the dampers and the inherent damping coefficient are known, the energy dissipation can be calculated by obtaining the integrals of the energy in equation (8).By assuming that the additional EDR provided by the dampers is uniform throughout all the stories, i.e.,ξ¯ai=ξ¯a, where ξ¯ai is the additional EDR of the i-th story, and the vibration of the structure is dominated by the 1st mode, equation (8) can be applied in the multi-degree-of-freedom (MDOF) system. The energy dissipated by the dampers at i-th story Edi is(9)Edi=∑j=1m∫0tfdjx˙idt=2ξ¯aikiω1∫0tx˙i2dt=2ξ¯akiω1∫0tx˙i2dt,where fdj denotes the damping force of the j-th damper at i-th story, m denotes the number of dampers installed at i-th story, x˙i denotes the relative velocity between i-th and (i-1)-th story, and ki denotes the stiffness of i-th story. And the energy dissipated by the inherent damping at the i-th story Eci is(10)Eci=∫0tcix˙i2dt=2ξ0kiω1∫0tx˙i2dt,where ci is the damping coefficient of the i-th story. The ratio of energy dissipated by the dampers and the inherent damping of all the stories can be(11)EdEc=∑i=1n1Edi∑i=1nEci=∑i=1n1∑j=1m∫0tfdj,ix˙idt∑i=1n∫0tcix˙i2dt=2ξ¯a/ω1∑i=1nki∫0tx˙i2dt2ξ0/ω1∑i=1nki∫0tx˙i2dt=ξ¯aξ0,where fdj,i is the j-th damper installed at the i-th story, n1 is the number of stories installed with dampers, and n is the number of stories.The additional EDR of the structureξ¯a can be obtained by rewriting equation (11) as follows:(12)ξ¯at=ξ0∑i=1n1∑j=1m∫0tfdj,ix˙idt∑i=1n∫0tcix˙i2dt.Generally, as the dampers are not configured uniformly throughout all the stories,ξ¯at denotes the average additional EDR of the entire structure in the duration 0,t.


### 2.3. Determination of Calculation Time

In equation (8), the additional EDR of the dampers ξ¯a is a function of the seismic duration. An earthquake is a random vibration process. The seismic response varies with time, and the additional EDR depends on the seismic response [31]. Therefore, the additional EDR ξ¯a varies with time during the entire duration 0,t. To reflect the energy dissipation capacity of the installed dampers during the most intensive period in the entire seismic motion duration, the calculation time should be determined. Generally, the dampers can ensure the maximum utilization of the energy dissipation capacity under large deformation, and large structural responses can be induced by intensive ground motion. Therefore, the time range covering the most intensive seismic motion can be adopted as the calculation time for the energy dissipation.Significant duration parameters are generally utilized to describe the duration characteristic of the ground motion recorders [32]. The significant duration is defined as the time intervals in which a specified amount of energy is dissipated. The integral of the square of the ground acceleration is called Arias intensity (IA), which is(13)IA=π2g∫0Ta2tdt,where g is the gravity acceleration, at is ground acceleration time history, and T is the entire duration of the ground motion record. Da5−75 is one of the common types of the signification duration [33], which is the time interval between 5% and 75% of IA and is employed in this study.Here, the calculation time is defined ast1,t2, where t1 is the moment arriving at 5% of IA and t2 is the moment arriving at 75% of IA. The additional EDR in the duration can be derived as follows:(14)ξ¯a,peak=ξ0∫t1t2fdx˙dt∫t1t2cx˙2dt=ξ0∫0t2fdx˙dt−∫0t1fdx˙dt∫0t2cx˙2dt−∫0t1cx˙2dt,where ξ¯a,peak is the additional EDR of the dampers in the calculation time t1,t2.Similarly, equation (12) can be modified to obtain the additional EDR for the MDOF system in t1,t2:(15)ξ¯a,peak=ξ0∑i=1n1∑j=1m∫t1t2fdj,ix˙idt∑i=1n∫t1t2cix˙i2dt=ξ0∑i=1n1∑j=1m∫0t2fdj,ix˙dt−∫0t1fdj,ix˙dt∑i=1n∫0t2cix˙2dt−∫0t1cix˙2dt.


### 2.4. Procedure

The procedure of the calculation method based on energy dissipation is as follows:(1)
Time history analysis: because the energy dissipation should be calculated during a certain period, this method is applicable only for time history analysis.(2)
Calculation of the energy dissipation: the energy dissipation of natural damping and the additional dampers can be calculated by obtaining the integrals of the energy in equation (8) for SDOF systems and in equation (12) for MDOF systems.(3)
Calculation of the additional EDR during the entire seismic motion duration: the time history of the additional EDR provided by the dampers can be obtained by equation (8) for SDOF systems and equation (12) for MDOF systems during the entire seismic motion duration 0,t.(4)
Determination of the durationt1,t2: the duration covering the intensive ground motion should be determined. t1 is the moment arriving at 5% of IA and t2 is the moment arriving at 75% of IA.(5)
Calculation of the additional EDR in the durationt1,t2: the energy dissipation of natural damping and the additional dampers during the duration t1,t2 can be calculated by obtaining the integrals in equation (14) for SDOF systems and equation (15) for MDOF systems. Then, the additional EDR induced by the dampers can be obtained using equation (14) and equation (15) for SDOF and MDOF systems, respectively.


## 2.1. Establishment of Energy Governing Equation

The simplified model of an SDOF system with an added viscous damper is shown in Figure1, where the viscous damper is simulated by the Maxwell model. According to structural dynamics, the differential equation of movement subjected to a horizontal earthquake can be written as follows [31]:(2)mx¨+cx˙+fsx,x˙+fdx,x˙=−mx¨g,where x, x˙, and x¨ are the relative displacement, velocity, and acceleration of the mass, respectively; x¨g is the ground acceleration; m is the mass; k is the stiffness of the system; c is the inherent damping coefficient of the system; fsx,x˙ is the elastic or inelastic resisting force of the system; and fdx,x˙ is the damping force of the damper. In the Maxwell model, ka is the dynamic stiffness of the damper and ca is the damping coefficient of the damper.Figure 1
SDOF system with the damper.By integratingx in equation (2) in the range of the earthquake duration [0, t], the energy response equation can be obtained:(3)m∫0tx¨dx+c∫0tx˙dx+∫0tfsdx+∫0tfddx=−m∫0tx¨gdx.Equation (3) can be rewritten as equation (4) in terms of energy:(4)Ek+Ec+Eh+Ee+Ed=Ein,where Ek is the kinetic energy, Ec is the energy dissipated by the inherent damping of the primary structure, and Eh and Ee are the hysteretic energy and elastic strain energy of the primary structure, respectively. Here, Eh+Ee=∫0tfsx˙dt; when the primary structure is in an elastic state, Eh=0. Ed is the energy dissipated by the dampers, and Ein is the energy input by ground motion.


## 2.2. Calculation of Additional EDR

In general, the energy dissipation capacity of the primary structure is indicated in terms of the EDR. Similarly, the energy dissipation capacity of additional devices can be indicated in terms of the additional EDRξeq. According to the classic definition method for the additional EDR, the energy dissipated by the dampers is equalized with that dissipated by an equivalent viscous system during one vibration cycle [20]:(5)Ed=∫02π/ωfdx˙dt=πx02ωceq=2πξeqkx02ωωn,where ω is the circular frequency of excitation, x0 is the maximum displacement of the structure in one cycle, ceq is the equivalent damping coefficient of the dampers in one cycle (ceq=2kξeq/ωn), and ωn is the natural circular frequency of the structure.This definition can be used to rewrite the equation for the additional EDR provided by the dampers during0,t as follows:(6)Edt=∫0tfdx˙dt=∫0tceqx˙2dt=2ξ¯akωn∫0tx˙2dt,where ξ¯a denotes the additional EDR provided by the dampers in the seismic duration 0,t and t denotes the duration of the seismic record.Meanwhile, the energy dissipated by the inherent damping of the primary structure can be expressed as follows:(7)Ect=∫0tcx˙2dt=2ξ0kωn∫0tx˙2dt,where ξ0 denotes the inherent damping ratio of the primary structure.From equations (6) and (7), the additional EDR of the dampers in the duration 0,t can be obtained:(8)ξ¯at=ξ0EdEc=ξ0∫0tfdx˙dt∫0tcx˙2dt.It is obvious that the additional EDRξ¯a can be obtained from the inherent damping ratio and the ratio of the energy dissipated by the dampers to that dissipated by the inherent damping during 0,t. It should be noted that the method focuses on the amount of energy dissipated and the manner to dissipate energy is not involved, which indicates that the method can be applied to various kinds of dampers. Equation (8) provides the expression for calculating the energy dissipation. When the damping force of the dampers and the inherent damping coefficient are known, the energy dissipation can be calculated by obtaining the integrals of the energy in equation (8).By assuming that the additional EDR provided by the dampers is uniform throughout all the stories, i.e.,ξ¯ai=ξ¯a, where ξ¯ai is the additional EDR of the i-th story, and the vibration of the structure is dominated by the 1st mode, equation (8) can be applied in the multi-degree-of-freedom (MDOF) system. The energy dissipated by the dampers at i-th story Edi is(9)Edi=∑j=1m∫0tfdjx˙idt=2ξ¯aikiω1∫0tx˙i2dt=2ξ¯akiω1∫0tx˙i2dt,where fdj denotes the damping force of the j-th damper at i-th story, m denotes the number of dampers installed at i-th story, x˙i denotes the relative velocity between i-th and (i-1)-th story, and ki denotes the stiffness of i-th story. And the energy dissipated by the inherent damping at the i-th story Eci is(10)Eci=∫0tcix˙i2dt=2ξ0kiω1∫0tx˙i2dt,where ci is the damping coefficient of the i-th story. The ratio of energy dissipated by the dampers and the inherent damping of all the stories can be(11)EdEc=∑i=1n1Edi∑i=1nEci=∑i=1n1∑j=1m∫0tfdj,ix˙idt∑i=1n∫0tcix˙i2dt=2ξ¯a/ω1∑i=1nki∫0tx˙i2dt2ξ0/ω1∑i=1nki∫0tx˙i2dt=ξ¯aξ0,where fdj,i is the j-th damper installed at the i-th story, n1 is the number of stories installed with dampers, and n is the number of stories.The additional EDR of the structureξ¯a can be obtained by rewriting equation (11) as follows:(12)ξ¯at=ξ0∑i=1n1∑j=1m∫0tfdj,ix˙idt∑i=1n∫0tcix˙i2dt.Generally, as the dampers are not configured uniformly throughout all the stories,ξ¯at denotes the average additional EDR of the entire structure in the duration 0,t.


## 2.3. Determination of Calculation Time

In equation (8), the additional EDR of the dampers ξ¯a is a function of the seismic duration. An earthquake is a random vibration process. The seismic response varies with time, and the additional EDR depends on the seismic response [31]. Therefore, the additional EDR ξ¯a varies with time during the entire duration 0,t. To reflect the energy dissipation capacity of the installed dampers during the most intensive period in the entire seismic motion duration, the calculation time should be determined. Generally, the dampers can ensure the maximum utilization of the energy dissipation capacity under large deformation, and large structural responses can be induced by intensive ground motion. Therefore, the time range covering the most intensive seismic motion can be adopted as the calculation time for the energy dissipation.Significant duration parameters are generally utilized to describe the duration characteristic of the ground motion recorders [32]. The significant duration is defined as the time intervals in which a specified amount of energy is dissipated. The integral of the square of the ground acceleration is called Arias intensity (IA), which is(13)IA=π2g∫0Ta2tdt,where g is the gravity acceleration, at is ground acceleration time history, and T is the entire duration of the ground motion record. Da5−75 is one of the common types of the signification duration [33], which is the time interval between 5% and 75% of IA and is employed in this study.Here, the calculation time is defined ast1,t2, where t1 is the moment arriving at 5% of IA and t2 is the moment arriving at 75% of IA. The additional EDR in the duration can be derived as follows:(14)ξ¯a,peak=ξ0∫t1t2fdx˙dt∫t1t2cx˙2dt=ξ0∫0t2fdx˙dt−∫0t1fdx˙dt∫0t2cx˙2dt−∫0t1cx˙2dt,where ξ¯a,peak is the additional EDR of the dampers in the calculation time t1,t2.Similarly, equation (12) can be modified to obtain the additional EDR for the MDOF system in t1,t2:(15)ξ¯a,peak=ξ0∑i=1n1∑j=1m∫t1t2fdj,ix˙idt∑i=1n∫t1t2cix˙i2dt=ξ0∑i=1n1∑j=1m∫0t2fdj,ix˙dt−∫0t1fdj,ix˙dt∑i=1n∫0t2cix˙2dt−∫0t1cix˙2dt.


## 2.4. Procedure

The procedure of the calculation method based on energy dissipation is as follows:(1)
Time history analysis: because the energy dissipation should be calculated during a certain period, this method is applicable only for time history analysis.(2)
Calculation of the energy dissipation: the energy dissipation of natural damping and the additional dampers can be calculated by obtaining the integrals of the energy in equation (8) for SDOF systems and in equation (12) for MDOF systems.(3)
Calculation of the additional EDR during the entire seismic motion duration: the time history of the additional EDR provided by the dampers can be obtained by equation (8) for SDOF systems and equation (12) for MDOF systems during the entire seismic motion duration 0,t.(4)
Determination of the durationt1,t2: the duration covering the intensive ground motion should be determined. t1 is the moment arriving at 5% of IA and t2 is the moment arriving at 75% of IA.(5)
Calculation of the additional EDR in the durationt1,t2: the energy dissipation of natural damping and the additional dampers during the duration t1,t2 can be calculated by obtaining the integrals in equation (14) for SDOF systems and equation (15) for MDOF systems. Then, the additional EDR induced by the dampers can be obtained using equation (14) and equation (15) for SDOF and MDOF systems, respectively.


## 3. Feasibility of the Proposed Method

To verify the feasibility of the proposed method, three viscously damped SDOF systems with different natural periods are introduced. In addition, a comparison study is conducted between this method and the classic strain energy method.


### 3.1. SDOF Systems

The natural periods of the three SDOF systems are 0.5 s, 1.0 s, and 2.0 s, respectively. Their masses are all equal to 204t, and the inherent damping ratio is 0.05. The stiffnesses of the three systems are 32214.4 kN/m, 8048.6 kN/m, and 2013.4 kN/m. A linear viscous damper is added to the three SDOF systems, respectively. The Maxwell model is employed to simulate the mechanical behavior of the damper. The damping coefficient of the damper is set as 250 kN/(m/s). The diagram of the SDOF systems is the same as that shown in Figure 1. The theoretical values of the additional EDRs induced by the added damper are 0.049, 0.097, and 0.195 for the three SDOF systems, respectively, which can be calculated using the differential equation of motion of the SDOF system.Three sinusoidal waves with different periods, an artificial ground motion AWX, and a natural ground motion NRX are selected for the time history analysis. The periods of the three sine waves are 0.5, 1, and 2 s. And the PGAs of all the excitations are 0.035 g. The normalized spectra of the AWX and NRX records are plotted in Figure2. The information of NRX is shown in Table 1, in which RSN is the record sequence number in the PEER ground motion database.Figure 2
Normalized spectra of AWX and NRX.Table 1
Information of NRX.
Name
RSN
Year
Magnitude
Resource
Station
t
1 (s)
t
2 (s)
NRX
760
1989
6.93
Loma Prieta
Foster City—Menhaden Court
5.79
13.055


### 3.2. Numerical Analysis

Because of the adoption of a linear viscous damper, equation (8) becomes(16)ξ¯at=ξ0∫0tfdx˙dt∫0tcx˙2dt=ξ0∫0tcax˙2dt∫0tcx˙2dt=ξ0cac=ca2mωn,which indicates that when the SDOF system is determined, the variation in the additional EDR ξ¯at is related to only the damping coefficient of the added damper ca.The time histories of the cumulative energy dissipated by the damper and inherent damping can be obtained from dynamic analysis. The additional EDR can be obtained using equation (8). The time histories of the cumulative energy of the SDOF system with the period equaling to 1.0 s are shown in Figure 3. The time histories of the additional EDR fluctuate in the first few seconds when the system is not in steady motion. After that, the values of the additional EDR tend to be stable. The values for all the ground motions are the same in this system, equaling 0.097.Figure 3
Cumulative energy dissipated by inherent damping and viscous damper and time histories of additional EDR of the SDOF system with the period equal to 1.0 s (the lines of inherent damping and the damper denote the cumulative energy dissipated by inherent damping and damper, respectively). (a) Sine wave,T=0.5s; (b) sine wave, T=1.0s; (c) sine wave, T=2.0s; (d) AWX; (e) NRX.
(a)
(b)
(c)
(d)
(e)The values of additional EDR of the three SDOF systems are listed in Table2. The frequencies of the excitations do not affect the calculation results of the additional EDR. The results show the feasibility of the proposed method.Table 2
Additional EDRs obtained by the strain energy method and proposed method for different cases.
Case
Strain energy method
Proposed method based on energy dissipation
T
n
=
0.5
s
T
n
=
1.0
s
T
n
=
2.0
s
T
n
=
0.5
s
T
n
=
1.0
s
T
n
=
2.0
s
T
=
0.5
s
0.048
0.193
0.774
0.048
0.096
0.192
T
=
1.0
s
0.024
0.097
0.389
0.048
0.097
0.194
T
=
2.0
s
0.012
0.049
0.195
0.049
0.097
0.195
AWX
0.046
0.084
0.240
0.048
0.097
0.194
NRX
0.037
0.087
0.244
0.048
0.097
0.195
Note. Tn stands for the natural period of the SDOF systems, and T stands for the period of excitation period of the sine waves.The additional EDRs are also calculated by the strain energy method for comparison, of which the results are listed in Table2. The additional EDRs are calculated using equation (1). Given that the added damper is a linear viscous damper, the energy dissipation can be calculated by the formula for the area of an ellipse:(17)ED=π⋅Fdmax⋅Δdmax,where Fdmax is the maximum output force of the damper in a cycle and Δdmax is the maximum deformation of the damper in a cycle.In the strain energy method, because the periods of the sine waves are known, the additional EDR in each vibration cycle can be obtained for the three cases involving the sine waves. Figure4 shows the variation in the additional EDR from the calculation results of the strain energy method for the SDOF system with the period equaling to 1.0 s. The hysteresis curve of the damper in this system under the sine wave with T=0.5s is shown in Figure 4(d); the curve indicates that the energy dissipation capacity of the damper is fully utilized.Figure 4
Time histories of additional EDRs obtained by the strain energy method under different sine waves of the SDOF system with the period equal to 1.0 s and hysteresis curve of the viscous damper for the system under the sine wave withT=0.5s. (a) Sine wave, T=0.5s; (b) sine wave, T=1.0s; (c) sine wave, T=2.0s; (d) hysteresis curve of damper.
(a)
(b)
(c)
(d)It can be observed from Table2 that only when the system is in resonance, the result of the additional EDR in this case is consistent with the result obtained from structural dynamics [20]. In the other cases, because the ratio of the excitation frequency to the structural natural frequency causes errors in the strain energy method, the results of the additional EDR are not an accurate reflection of the energy dissipation capacity of the damper.The comparison study between the proposed method and strain energy method shows that the proposed method not only is an easier approach to obtain the additional EDR but also provides calculation results that are more stable. This method can avoid the influence of the excitation frequency on the calculation result, which is a major problem of the strain energy method.


## 3.1. SDOF Systems

The natural periods of the three SDOF systems are 0.5 s, 1.0 s, and 2.0 s, respectively. Their masses are all equal to 204t, and the inherent damping ratio is 0.05. The stiffnesses of the three systems are 32214.4 kN/m, 8048.6 kN/m, and 2013.4 kN/m. A linear viscous damper is added to the three SDOF systems, respectively. The Maxwell model is employed to simulate the mechanical behavior of the damper. The damping coefficient of the damper is set as 250 kN/(m/s). The diagram of the SDOF systems is the same as that shown in Figure 1. The theoretical values of the additional EDRs induced by the added damper are 0.049, 0.097, and 0.195 for the three SDOF systems, respectively, which can be calculated using the differential equation of motion of the SDOF system.Three sinusoidal waves with different periods, an artificial ground motion AWX, and a natural ground motion NRX are selected for the time history analysis. The periods of the three sine waves are 0.5, 1, and 2 s. And the PGAs of all the excitations are 0.035 g. The normalized spectra of the AWX and NRX records are plotted in Figure2. The information of NRX is shown in Table 1, in which RSN is the record sequence number in the PEER ground motion database.Figure 2
Normalized spectra of AWX and NRX.Table 1
Information of NRX.
Name
RSN
Year
Magnitude
Resource
Station
t
1 (s)
t
2 (s)
NRX
760
1989
6.93
Loma Prieta
Foster City—Menhaden Court
5.79
13.055


## 3.2. Numerical Analysis

Because of the adoption of a linear viscous damper, equation (8) becomes(16)ξ¯at=ξ0∫0tfdx˙dt∫0tcx˙2dt=ξ0∫0tcax˙2dt∫0tcx˙2dt=ξ0cac=ca2mωn,which indicates that when the SDOF system is determined, the variation in the additional EDR ξ¯at is related to only the damping coefficient of the added damper ca.The time histories of the cumulative energy dissipated by the damper and inherent damping can be obtained from dynamic analysis. The additional EDR can be obtained using equation (8). The time histories of the cumulative energy of the SDOF system with the period equaling to 1.0 s are shown in Figure 3. The time histories of the additional EDR fluctuate in the first few seconds when the system is not in steady motion. After that, the values of the additional EDR tend to be stable. The values for all the ground motions are the same in this system, equaling 0.097.Figure 3
Cumulative energy dissipated by inherent damping and viscous damper and time histories of additional EDR of the SDOF system with the period equal to 1.0 s (the lines of inherent damping and the damper denote the cumulative energy dissipated by inherent damping and damper, respectively). (a) Sine wave,T=0.5s; (b) sine wave, T=1.0s; (c) sine wave, T=2.0s; (d) AWX; (e) NRX.
(a)
(b)
(c)
(d)
(e)The values of additional EDR of the three SDOF systems are listed in Table2. The frequencies of the excitations do not affect the calculation results of the additional EDR. The results show the feasibility of the proposed method.Table 2
Additional EDRs obtained by the strain energy method and proposed method for different cases.
Case
Strain energy method
Proposed method based on energy dissipation
T
n
=
0.5
s
T
n
=
1.0
s
T
n
=
2.0
s
T
n
=
0.5
s
T
n
=
1.0
s
T
n
=
2.0
s
T
=
0.5
s
0.048
0.193
0.774
0.048
0.096
0.192
T
=
1.0
s
0.024
0.097
0.389
0.048
0.097
0.194
T
=
2.0
s
0.012
0.049
0.195
0.049
0.097
0.195
AWX
0.046
0.084
0.240
0.048
0.097
0.194
NRX
0.037
0.087
0.244
0.048
0.097
0.195
Note. Tn stands for the natural period of the SDOF systems, and T stands for the period of excitation period of the sine waves.The additional EDRs are also calculated by the strain energy method for comparison, of which the results are listed in Table2. The additional EDRs are calculated using equation (1). Given that the added damper is a linear viscous damper, the energy dissipation can be calculated by the formula for the area of an ellipse:(17)ED=π⋅Fdmax⋅Δdmax,where Fdmax is the maximum output force of the damper in a cycle and Δdmax is the maximum deformation of the damper in a cycle.In the strain energy method, because the periods of the sine waves are known, the additional EDR in each vibration cycle can be obtained for the three cases involving the sine waves. Figure4 shows the variation in the additional EDR from the calculation results of the strain energy method for the SDOF system with the period equaling to 1.0 s. The hysteresis curve of the damper in this system under the sine wave with T=0.5s is shown in Figure 4(d); the curve indicates that the energy dissipation capacity of the damper is fully utilized.Figure 4
Time histories of additional EDRs obtained by the strain energy method under different sine waves of the SDOF system with the period equal to 1.0 s and hysteresis curve of the viscous damper for the system under the sine wave withT=0.5s. (a) Sine wave, T=0.5s; (b) sine wave, T=1.0s; (c) sine wave, T=2.0s; (d) hysteresis curve of damper.
(a)
(b)
(c)
(d)It can be observed from Table2 that only when the system is in resonance, the result of the additional EDR in this case is consistent with the result obtained from structural dynamics [20]. In the other cases, because the ratio of the excitation frequency to the structural natural frequency causes errors in the strain energy method, the results of the additional EDR are not an accurate reflection of the energy dissipation capacity of the damper.The comparison study between the proposed method and strain energy method shows that the proposed method not only is an easier approach to obtain the additional EDR but also provides calculation results that are more stable. This method can avoid the influence of the excitation frequency on the calculation result, which is a major problem of the strain energy method.


## 4. Example Illustration

### 4.1. Basic Information

A benchmark reinforcement concrete model is used to illustrate the usefulness of this method in practical engineering applications [34]. According to the Chinese code for the seismic design of buildings [22], the design PGA for dynamic analysis is 0.2 g and the characteristic period of the site soil is 0.4 s. The structural plan and elevation are depicted in Figure 5, and the dimensions of the cross sections of the members are listed in Table 3. The model of the structure is created using the software SAP2000 [35]. The beams and columns are modeled as elastic bars with plastic hinges at the two ends. The mechanical model of the hinges is the ideal rigid-plastic model, and the types of plastic hinges are moment M3 for the beams and interacting P-M2-M3 for the columns. The material of the members is reinforced concrete, and the design strength of the concrete is 30 MPa according to the code for the design of concrete structures [36]. The thin shell element is employed to simulate the floor slabs. The thickness of all the floors is set to be 0 mm so as to achieve optimal load transmission. The floors are subjected to distributed loads. The distributed dead loads are 5.5kN/m2, including the self-weight of the floor, for all the stories. The live loads are 3.0kN/m2 for stories 1–5 and 5.0kN/m2 for story 6. The story heights and loads are listed in Table 4. The Rayleigh damping is adopted in the structural model for dynamic analysis. The damping ratio is set to be 0.05 for the reinforcement concrete structure for the first and the fourth modes.Figure 5
Structural plan and elevation (MYDs are installed in theY direction) (unit: m): (a) Structural plan; (b) structural elevation.
(a)
(b)Table 3
Dimensions of member cross sections.
Story
1 and 2
3 and 4
5 and 6
Beam
350
mm
×
600
mm
350
mm
×
600
mm
350
mm
×
600
mm
Column
700
mm
×
700
mm
600
mm
×
600
mm
500
mm
×
500
mmTable 4
Story heights and loads.
Story
Height (m)
Dead load (kN/m2)
Live load (kN/m2)
1
4.6
5.5
3.0
2
4.2
5.5
3.0
3
3.6
5.5
3.0
4
3.6
5.5
3.0
5
3.6
5.5
3.0
6
4.6
5.5
5.0Two MYDs are installed at each story along the axes X3 and X7. The damper parameters are listed in Table5. The damper installation diagram is also shown in Figure 5.Table 5
Parameters of dampers at each story.
Story
Initial stiffnessKI (kN/m)
Yielding deformationΔdy (mm)
Postyielding stiffness ratioαd
1
1800
1
0.02
2
1800
1
0.02
3
1500
1
0.02
4
1500
1
0.02
5
1200
1
0.02
6
600
1
0.02Next, nonlinear dynamic analysis is performed. To comply with the specifications of the Chinese code [22], seven earthquake ground motions are selected, of which five are natural ground motions. The natural ground motions are selected from the ground motion database of the Pacific Earthquake Engineering Research Center; the spectra of these ground motions are matched with the design acceleration spectrum [37]. The artificial ground motions are generated using a trigonometric method [38] and adjusted to match the design response spectrum [39, 40]. For each earthquake, t1 and t2 are determined according to the definitions in Section 2.3. The details of the ground motions are listed in Table 6, and the normalized spectra of these ground motions are shown in Figure 6. The PGAs of the ground motions are set to 0.2 g.Table 6
Details of input ground motion records.
Name
RSN
Year
Magnitude
Resource
Station
t
1 (s)
t
2 (s)
NW1
1768
1999
7.13
Hector Mine
Barstow
9.56
19.52
NW2
1838
1999
7.13
Hector Mine
Whitewater Trout Farm
6.36
16.38
NW3
4489
1999
6.3
L’Aquila, Italy
Montereale
9.34
23.44
NW4
6896
1999
7
Darfield, New Zealand
DORC
7.28
14.08
NW5
167
1999
6.53
Imperial Valley-06
Compuertas
15.24
31.36
AW1
—
—
—
—
5.9
20.58
AW2
—
—
—
—
6.42
22.4Figure 6
Spectra of normalized input ground motion records.


### 4.2. Additional EDR Calculation

According to the proposed method, the additional EDR for the entire earthquake duration0,t and t1,t2 is obtained using equations (12) and (15). In these two equations, the energy dissipated by inherent damping and the added dampers can be retrieved from the results of dynamic analysis in SAP2000. As a result, the time histories of ξa and the energy dissipated by natural damping and the MYDs can be easily obtained, which are shown in Figure 7.Figure 7
Cumulative energy dissipated by inherent damping and MYDs and time histories of additional EDR (the lines of inherent damping and the damper denote the cumulative energy dissipated by inherent damping and damper, respectively): (a) NW1; (b) NW2; (c) NW3; (d) NW4; (e) NW5; (f) AW1; (g) AW2.
(a)
(b)
(c)
(d)
(e)
(f)
(g)For multistory structures, unlike the SDOF system, the time histories of additional EDR does not keep in a stable value, which fluctuate with the time throughout the duration of entire seismic motions. It can be observed that, aftert1, the inherent damping and the additional dampers started to dissipate the input energy and additional EDR increases sharply in the same time. And the energy dissipation of the dampers almost stops, and the fluctuation of the additional EDR tends to be gentle after the time t2. It indicates that the dampers are fully utilized within the significant duration of ground motions.Because of the difference of the frequency component and the duration of the ground motions, the responses of the structure, which are influenced by the higher modes, will be different. Moreover, except for the linear viscous dampers, the additional EDR induced by dampers are closely related to the amplitude of structural response [31], which results in the fluctuation of the time histories of additional EDR and the slight difference in the values of additional EDR obtained from different excitations.Next, the strain energy method is adopted to calculate the additional EDR induced by the MYDs using equation (1). In the equation, the dissipated energy ED and maximum strain energy ESo in one cycle are calculated as follows:(18)ED=4⋅∑j=1MKIj⋅Δdyj2⋅1−αdjμdj−1,(19)ESo=12∑i=1NQi⋅Δi+12∑jMKIj⋅Δdyj2⋅μdj⋅1+μdj⋅αdj−αdj,where KIj denotes the initial stiffness of the j-th damper, Δdyj denotes the yielding deformation of the j-th damper, αdj denotes the postyielding stiffness ratio of the j-th damper, and μdj denotes the ductility factor of the j-th damper. Here, μdj=Δdj/Δdyj, where Δdj denotes the deformation of the j-th damper, Qi denotes the shear force at level i, Δi denotes the interstory displacement between level i and level i-1, M denotes the number of dampers, and N denotes the number of stories. Figure 8 shows the hysteresis curves of the MYDs at levels 1, 3, and 5 under the seismic motion AW1. The curves indicate that the energy dissipation capacity of the MYDs is fully utilized. The values obtained from equations (12), (15), and (1) are listed in Table 7.Figure 8
Hysteresis curves of MYDs under AW1: (a) damper at story 1; (b) damper at story 3; (c) damper at story 5.
(a)
(b)
(c)Table 7
Additional EDRs obtained by different approaches.
Case
ξ
a,end
ξ
a,peak
ξ
a,strain
NW1
0.085
0.118
0.213
NW2
0.055
0.068
0.185
NW3
0.053
0.074
0.204
NW4
0.053
0.075
0.170
NW5
0.061
0.078
0.184
AW1
0.079
0.083
0.210
AW2
0.059
0.073
0.201
Average
0.063
0.081
0.195
Note:ξa,end is the final value of ξ¯a, ξa,peak is the additional EDR induced by the dampers during the duration t1,t2, and ξa,strain is the value obtained from the strain energy method.


### 4.3. Comparative Analysis

The responses of the structures with the additional EDR obtained from different methods (ST-end, ST-peak, and ST-strain) are compared with the response of the original structure (ST-original) to determine whether the results of the additional EDR can reflect the energy dissipation capacity of the MYDs. ST-end denotes the structure withξa,end, ST-peak denotes the structure with ξa,peak, and ST-strain denotes the structure with ξa,strain. In these structures, the energy dissipation capacity of the MYDs is represented by the additional EDR obtained from different approaches and the Rayleigh damping are also adopted in these structures, of which the values are 0.05+ξa for the first and fourth modes. The MYDs are modeled as elastic bars with equivalent stiffness, which can be calculated as follows [21]:(20)kj,eff=Fdjmax−FdjminΔdjmax−Δdjmin,where kj,eff denotes the equivalent stiffness of the j-th damper; Fdjmax and Fdjmin denote the maximum and minimum damping forces of the j-th damper, respectively, under a certain earthquake; and Δjmax and Δjmin denote the maximum and minimum deformations of the j-th damper, respectively, under a certain earthquake. Here, kj,eff is the average of the equivalent stiffnesses under the seven seismic ground motions.Next, the dynamic analysis is carried out. The story shears and story drifts of the three structures are compared with that of the original structure, as shown in Figures9 and 10, respectively.Figure 9
Story shears.Figure 10
Story drifts.The above comparison clearly indicates that all the responses are close. The response of ST-peak is closest to that of ST-original. The response of ST-end is larger than that of ST-original, whereas the response of ST-strain is smaller than that of ST-original. Because the response of ST-end is larger than that of ST-original, the value ofζa,end is small. Generally, during most of the entire earthquake duration, the seismic responses of the structures are considerably smaller than the maximum responses of the structures. Therefore, the MYDs are not always utilizing their maximum capacity during the entire seismic duration. The energy dissipation capacity of the dampers can be underestimated within the entire earthquake duration. Therefore, the additional EDR is estimated during the most intensive period in the entire seismic motion duration.The responses of ST-strain are smaller than those of ST-original, which means thatξa,strain is too large. Currently, the strain energy method is the most common approach used to estimate the additional EDR. However, for practical applications, the method not only is complicated in terms of calculation but also overestimates the additional EDR in some cases. In the example considered in this study, the error in ξa,strain may be attributed to the following causes:(1)
The seismic motion is a stochastic process, which can be represented by a combination of a series of sine waves with a wide range of frequencies. As discussed in Section3.2, the frequency of the excitation will cause obvious errors in the results of the additional EDR obtained by the strain energy method.(2)
The energy dissipation is obtained according to the maximum deformations of the dampers based on equation (18). However, the maximum deformations of the dampers in different levels usually do not occur at the same time during the vibration. The sum of the energy dissipation values of the dampers based on the maximum deformations during the entire earthquake will overestimate the mitigation effect of the MYDs.(3)
The strain energy method is derived using an elastic SDOF system. When the primary structure is driven into the inelastic range, the method can cause obvious errors in the calculation of the strain energy and dissipated energy; hence, the calculation equations for these energies should be modified.


## 4.1. Basic Information

A benchmark reinforcement concrete model is used to illustrate the usefulness of this method in practical engineering applications [34]. According to the Chinese code for the seismic design of buildings [22], the design PGA for dynamic analysis is 0.2 g and the characteristic period of the site soil is 0.4 s. The structural plan and elevation are depicted in Figure 5, and the dimensions of the cross sections of the members are listed in Table 3. The model of the structure is created using the software SAP2000 [35]. The beams and columns are modeled as elastic bars with plastic hinges at the two ends. The mechanical model of the hinges is the ideal rigid-plastic model, and the types of plastic hinges are moment M3 for the beams and interacting P-M2-M3 for the columns. The material of the members is reinforced concrete, and the design strength of the concrete is 30 MPa according to the code for the design of concrete structures [36]. The thin shell element is employed to simulate the floor slabs. The thickness of all the floors is set to be 0 mm so as to achieve optimal load transmission. The floors are subjected to distributed loads. The distributed dead loads are 5.5kN/m2, including the self-weight of the floor, for all the stories. The live loads are 3.0kN/m2 for stories 1–5 and 5.0kN/m2 for story 6. The story heights and loads are listed in Table 4. The Rayleigh damping is adopted in the structural model for dynamic analysis. The damping ratio is set to be 0.05 for the reinforcement concrete structure for the first and the fourth modes.Figure 5
Structural plan and elevation (MYDs are installed in theY direction) (unit: m): (a) Structural plan; (b) structural elevation.
(a)
(b)Table 3
Dimensions of member cross sections.
Story
1 and 2
3 and 4
5 and 6
Beam
350
mm
×
600
mm
350
mm
×
600
mm
350
mm
×
600
mm
Column
700
mm
×
700
mm
600
mm
×
600
mm
500
mm
×
500
mmTable 4
Story heights and loads.
Story
Height (m)
Dead load (kN/m2)
Live load (kN/m2)
1
4.6
5.5
3.0
2
4.2
5.5
3.0
3
3.6
5.5
3.0
4
3.6
5.5
3.0
5
3.6
5.5
3.0
6
4.6
5.5
5.0Two MYDs are installed at each story along the axes X3 and X7. The damper parameters are listed in Table5. The damper installation diagram is also shown in Figure 5.Table 5
Parameters of dampers at each story.
Story
Initial stiffnessKI (kN/m)
Yielding deformationΔdy (mm)
Postyielding stiffness ratioαd
1
1800
1
0.02
2
1800
1
0.02
3
1500
1
0.02
4
1500
1
0.02
5
1200
1
0.02
6
600
1
0.02Next, nonlinear dynamic analysis is performed. To comply with the specifications of the Chinese code [22], seven earthquake ground motions are selected, of which five are natural ground motions. The natural ground motions are selected from the ground motion database of the Pacific Earthquake Engineering Research Center; the spectra of these ground motions are matched with the design acceleration spectrum [37]. The artificial ground motions are generated using a trigonometric method [38] and adjusted to match the design response spectrum [39, 40]. For each earthquake, t1 and t2 are determined according to the definitions in Section 2.3. The details of the ground motions are listed in Table 6, and the normalized spectra of these ground motions are shown in Figure 6. The PGAs of the ground motions are set to 0.2 g.Table 6
Details of input ground motion records.
Name
RSN
Year
Magnitude
Resource
Station
t
1 (s)
t
2 (s)
NW1
1768
1999
7.13
Hector Mine
Barstow
9.56
19.52
NW2
1838
1999
7.13
Hector Mine
Whitewater Trout Farm
6.36
16.38
NW3
4489
1999
6.3
L’Aquila, Italy
Montereale
9.34
23.44
NW4
6896
1999
7
Darfield, New Zealand
DORC
7.28
14.08
NW5
167
1999
6.53
Imperial Valley-06
Compuertas
15.24
31.36
AW1
—
—
—
—
5.9
20.58
AW2
—
—
—
—
6.42
22.4Figure 6
Spectra of normalized input ground motion records.


## 4.2. Additional EDR Calculation

According to the proposed method, the additional EDR for the entire earthquake duration0,t and t1,t2 is obtained using equations (12) and (15). In these two equations, the energy dissipated by inherent damping and the added dampers can be retrieved from the results of dynamic analysis in SAP2000. As a result, the time histories of ξa and the energy dissipated by natural damping and the MYDs can be easily obtained, which are shown in Figure 7.Figure 7
Cumulative energy dissipated by inherent damping and MYDs and time histories of additional EDR (the lines of inherent damping and the damper denote the cumulative energy dissipated by inherent damping and damper, respectively): (a) NW1; (b) NW2; (c) NW3; (d) NW4; (e) NW5; (f) AW1; (g) AW2.
(a)
(b)
(c)
(d)
(e)
(f)
(g)For multistory structures, unlike the SDOF system, the time histories of additional EDR does not keep in a stable value, which fluctuate with the time throughout the duration of entire seismic motions. It can be observed that, aftert1, the inherent damping and the additional dampers started to dissipate the input energy and additional EDR increases sharply in the same time. And the energy dissipation of the dampers almost stops, and the fluctuation of the additional EDR tends to be gentle after the time t2. It indicates that the dampers are fully utilized within the significant duration of ground motions.Because of the difference of the frequency component and the duration of the ground motions, the responses of the structure, which are influenced by the higher modes, will be different. Moreover, except for the linear viscous dampers, the additional EDR induced by dampers are closely related to the amplitude of structural response [31], which results in the fluctuation of the time histories of additional EDR and the slight difference in the values of additional EDR obtained from different excitations.Next, the strain energy method is adopted to calculate the additional EDR induced by the MYDs using equation (1). In the equation, the dissipated energy ED and maximum strain energy ESo in one cycle are calculated as follows:(18)ED=4⋅∑j=1MKIj⋅Δdyj2⋅1−αdjμdj−1,(19)ESo=12∑i=1NQi⋅Δi+12∑jMKIj⋅Δdyj2⋅μdj⋅1+μdj⋅αdj−αdj,where KIj denotes the initial stiffness of the j-th damper, Δdyj denotes the yielding deformation of the j-th damper, αdj denotes the postyielding stiffness ratio of the j-th damper, and μdj denotes the ductility factor of the j-th damper. Here, μdj=Δdj/Δdyj, where Δdj denotes the deformation of the j-th damper, Qi denotes the shear force at level i, Δi denotes the interstory displacement between level i and level i-1, M denotes the number of dampers, and N denotes the number of stories. Figure 8 shows the hysteresis curves of the MYDs at levels 1, 3, and 5 under the seismic motion AW1. The curves indicate that the energy dissipation capacity of the MYDs is fully utilized. The values obtained from equations (12), (15), and (1) are listed in Table 7.Figure 8
Hysteresis curves of MYDs under AW1: (a) damper at story 1; (b) damper at story 3; (c) damper at story 5.
(a)
(b)
(c)Table 7
Additional EDRs obtained by different approaches.
Case
ξ
a,end
ξ
a,peak
ξ
a,strain
NW1
0.085
0.118
0.213
NW2
0.055
0.068
0.185
NW3
0.053
0.074
0.204
NW4
0.053
0.075
0.170
NW5
0.061
0.078
0.184
AW1
0.079
0.083
0.210
AW2
0.059
0.073
0.201
Average
0.063
0.081
0.195
Note:ξa,end is the final value of ξ¯a, ξa,peak is the additional EDR induced by the dampers during the duration t1,t2, and ξa,strain is the value obtained from the strain energy method.


## 4.3. Comparative Analysis

The responses of the structures with the additional EDR obtained from different methods (ST-end, ST-peak, and ST-strain) are compared with the response of the original structure (ST-original) to determine whether the results of the additional EDR can reflect the energy dissipation capacity of the MYDs. ST-end denotes the structure withξa,end, ST-peak denotes the structure with ξa,peak, and ST-strain denotes the structure with ξa,strain. In these structures, the energy dissipation capacity of the MYDs is represented by the additional EDR obtained from different approaches and the Rayleigh damping are also adopted in these structures, of which the values are 0.05+ξa for the first and fourth modes. The MYDs are modeled as elastic bars with equivalent stiffness, which can be calculated as follows [21]:(20)kj,eff=Fdjmax−FdjminΔdjmax−Δdjmin,where kj,eff denotes the equivalent stiffness of the j-th damper; Fdjmax and Fdjmin denote the maximum and minimum damping forces of the j-th damper, respectively, under a certain earthquake; and Δjmax and Δjmin denote the maximum and minimum deformations of the j-th damper, respectively, under a certain earthquake. Here, kj,eff is the average of the equivalent stiffnesses under the seven seismic ground motions.Next, the dynamic analysis is carried out. The story shears and story drifts of the three structures are compared with that of the original structure, as shown in Figures9 and 10, respectively.Figure 9
Story shears.Figure 10
Story drifts.The above comparison clearly indicates that all the responses are close. The response of ST-peak is closest to that of ST-original. The response of ST-end is larger than that of ST-original, whereas the response of ST-strain is smaller than that of ST-original. Because the response of ST-end is larger than that of ST-original, the value ofζa,end is small. Generally, during most of the entire earthquake duration, the seismic responses of the structures are considerably smaller than the maximum responses of the structures. Therefore, the MYDs are not always utilizing their maximum capacity during the entire seismic duration. The energy dissipation capacity of the dampers can be underestimated within the entire earthquake duration. Therefore, the additional EDR is estimated during the most intensive period in the entire seismic motion duration.The responses of ST-strain are smaller than those of ST-original, which means thatξa,strain is too large. Currently, the strain energy method is the most common approach used to estimate the additional EDR. However, for practical applications, the method not only is complicated in terms of calculation but also overestimates the additional EDR in some cases. In the example considered in this study, the error in ξa,strain may be attributed to the following causes:(1)
The seismic motion is a stochastic process, which can be represented by a combination of a series of sine waves with a wide range of frequencies. As discussed in Section3.2, the frequency of the excitation will cause obvious errors in the results of the additional EDR obtained by the strain energy method.(2)
The energy dissipation is obtained according to the maximum deformations of the dampers based on equation (18). However, the maximum deformations of the dampers in different levels usually do not occur at the same time during the vibration. The sum of the energy dissipation values of the dampers based on the maximum deformations during the entire earthquake will overestimate the mitigation effect of the MYDs.(3)
The strain energy method is derived using an elastic SDOF system. When the primary structure is driven into the inelastic range, the method can cause obvious errors in the calculation of the strain energy and dissipated energy; hence, the calculation equations for these energies should be modified.


## 5. Conclusion

In this paper, a simple method is proposed for the fast estimation of the additional EDR of damped structures. The method is derived from the definition of the EDR in structural dynamics. Through time history analysis, the ratio of the energy consumption of inherent damping to that of the added dampers can be obtained under external excitation. The additional EDR can be estimated using the inherent damping ratio and energy consumption ratio. Based on the study results, the following conclusions can be drawn:(1)
The concept of the method is explicit, intuitive, and easy to understand for structural practitioners. By using the inherent damping ratio and the ratio of the energy dissipated by dampers to that dissipated by inherent damping, the additional EDR can be accurately estimated.(2)
Unlike the strain energy method, the proposed method can eliminate the influence of the excitation frequency on the result of the additional EDR, which means that the additional EDR obtained by this method is a better reflection of the energy dissipation capacity of dampers.(3)
The time-varying characteristic of the additional EDR is revealed in this study. Except for the linear viscous dampers, the additional EDR induced by the dampers in the structure is not constant and fluctuates in the duration of ground motions. To improve the accuracy, the estimation of additional EDR should be calculated within the significant duration of ground motion.(4)
The proposed method focuses on the energy dissipation of inherent damping and the dampers and hence is not related to the state of the structure. Therefore, this method is applicable to not only elastic structures but also structures in the inelastic state.



---

*Source: 8052413-2019-08-15.xml*