Capacity Spectrum-Based Retroﬁtting Method and Quick Design for Viscously Damped Structures by Utilizing the Concept of Uniform Damping Ratio

: Viscous dampers have proven to be effective in enhancing the seismic performance of existing structures. Despite this, there is still a need for rapid and simpliﬁed design methods and formulae for viscous dampers that can take into account the elastic–plastic performance of structures. This study introduces a retroﬁt design method for existing structures using viscous dampers, based on the concept of uniform damping ratio (UDR), with the aim of fully utilizing each damper. The UDR concept assumes that each damper in the structure provides the same UDR when subjected to seismic excitations of identical intensity. In this method, the ﬁrst step involves deﬁning the equivalent damping ratio (EDR) of the damper. Then, based on the capacity spectrum of the structure, the response mitigation ratio can be determined, which helps to determine the additional EDR required from the dampers. Once the UDR and additional EDR from the viscous damper have been determined, the parameters of the dampers at each story can be rapidly obtained. To demonstrate the effectiveness of this method, a six-story reinforced concrete frame was utilized as a benchmark structure. A comparison between this UDR-based approach and a traditional design approach was also conducted. The study ﬁndings reveal that the UDR concept enables the maximum utilization of energy dissipation capacity of viscous dampers installed in the structure, leading to a more effective and economical design approach.


Basic Concept of EDR for Single Viscous Damper
A viscous damper is a kind of hydraulic device that takes advantage of the viscosity of a fluid to dissipate the input energy. The ideal viscous fluid is able to dissipate energy and cannot store energy at all, which does not exist in reality. For general dampers, silicon oil is typically selected as the viscous material. However, silicon oil has great compressibility, which leads to the dynamic stiffness of dampers under vibration processes. The dynamic stiffness depends on the type of silicon oil, the frequency of excitation, and the size of the valve. For products produced by different manufacturers, stiffness varies. According to an experimental study, dynamic stiffness is related to the damping coefficient of the damper as follows [32]: where K d represents the dynamic stiffness, β represents the dynamic stiffness coefficient, and C d represents the damping coefficient. As the dampers can both store and dissipate energy under vibration processes, the Maxwell model can be used [46]. The Maxwell model consists of a spring element and a damping element, which are depicted in Figure 1. The restoring force of the Maxwell model is shown in Equation (2): where α is the damping exponent of the viscous damper, F d is the output force of the Maxwell model, d K is the relative deformation of the spring element, and . d c is the relative velocity of the damping element.
where d is the total deformation of the Maxwell model and c d is the deformation of the damping element. Therefore, the EDR of a single damper can be defined according to the concept of EDR in structural dynamics [47]. The most common method for defining the EDR is to make the energy dissipated by the damper in a vibration cycle equal to an equivalent viscous system, which is given by the following equation [45]:

Kd Cd
Spring element Damping element Maxwell model The total deformation of the Maxwell model, d, comprises the deformation of the spring element and the damping element: where d is the total deformation of the Maxwell model and d c is the deformation of the damping element. Therefore, the EDR of a single damper can be defined according to the concept of EDR in structural dynamics [47]. The most common method for defining the EDR is to make the energy dissipated by the damper in a vibration cycle equal to an equivalent viscous system, which is given by the following equation [45]: where ζ d represents the EDR of the viscous damper, E d is the energy dissipated by the damper in one cycle, and E Sd is the maximum strain energy of the damper structure in one cycle. For most kinds of viscous dampers used for resisting seismic excitation, the damping exponent is commonly in the range of 0.1~0.4. Figure 2 shows a typical hysteresis curve of a nonlinear viscous damper with a damping exponent of 0.2. For practical engineering, it would be complicated to calculate the dissipated and restored energy of the viscous damper. To simplify the calculation, a simplified parallelogram encompassing the hysteresis curve of the damper is introduced, which is represented by the dashed line in Figure 2. The energy dissipation capacity of the viscous damper can be approximately obtained by the area of the simplified parallelogram.  Therefore, the EDR of a single damper can be defined according to the concept of EDR in structural dynamics [47]. The most common method for defining the EDR is to make the energy dissipated by the damper in a vibration cycle equal to an equivalent viscous system, which is given by the following equation [45]:  Figure 2 shows a typical hysteresis curve of a nonlinear viscous damper with a damping exponent of 0.2. For practical engineering, it would be complicated to calculate the dissipated and restored energy of the viscous damper. To simplify the calculation, a simplified parallelogram encompassing the hysteresis curve of the damper is introduced, which is represented by the dashed line in Figure  2. The energy dissipation capacity of the viscous damper can be approximately obtained by the area of the simplified parallelogram. To determine if the approximation of the calculation of energy dissipation is acceptable, the area error between the ideal hysteresis curve and the simplified parallelogram is explored. First, without considering the dynamic stiffness of the damper, namely K d = ∞ and d = d c , the hysteresis curve is plotted, as shown in Figure 3. The area of the hysteresis curve can be calculated as follows [48]: where ω is the natural circular frequency of the system with dampers, d c,max is the maximum deformation of the damping element in one cycle, and Γ is the gamma function. Without considering the dynamic stiffness of damper K d , the simplified parallelogram becomes rectangular, which is also shown in Figure 3. The area of the parallelogram can be obtained as: where ϕ is the reduction coefficient of the parallelogram. It is obvious that the error is only related to the damping exponen damper. Figure 4 shows that the ratio varies with α in the range o expression of the reduction coefficient of the parallelogram ϕ is give It is obvious that the error is only related to the damping exponent α of the viscous damper. Figure 4 shows that the ratio varies with α in the range of 0 to 1. The fitting expression of the reduction coefficient of the parallelogram ϕ is given as follows: Through ϕ, the energy dissipation of the damper can be estimated, which is shown as: Buildings 2023, 13, x FOR PEER REVIEW 7 of 22    Next, when the dynamic stiffness K d is considered, d no longer equals d c . The calculation of E d is same as Equation (5), but d cmax in Equation (5) is difficult to obtain. E d,paral can be approximately obtained by: where d 0 is the deformation of the Maxwell model at zero force, which equals the maximum deformation of the damping element, namely d 0 = d c,max . In general, obtaining d c,max is complicated; however, we can obtain the maximum deformation of the Maxwell model d max easily through analysis. In this light, the relation between d 0 and d max can be explored. This relationship depends on the damping exponent α, the dynamic stiffness K d , the circular frequency ω of the structure, and the maximum deformation of the Maxwell model d max .
The first three parameters are decided by the type of viscous damper and the structural characteristics. Therefore, κ varies with d max . The JSSI Manual provides a fitting expression with high accuracy [32]. Before the expression is given, the pseudo loss of stiffness of the viscous damper is defined: Next, when the dynamic stiffness d K is considered, d no longer equals c d . The calculation of d E is same as Equation (5), but max  [32]. Before the expression is given, the pseudo loss of stiffness of the viscous damper is defined: where p d K ″ represents the loss of the stiffness of the viscous damper.
The ratio of the loss of stiffness to the dynamic stiffness is: and the ratio of ,max c d to max d is given as follows [32]: where κ represents the ratio of ,max c d to max d . The applicable condition of Equation (13) is / 5  Next, when the dynamic stiffness d K is considered, d no longer equals c d . The calculation of d E is same as Equation (5), but max c d in Equation (5) [32]. Before the expression is given, the pseudo loss of stiffness of the viscous damper is defined: where p d K ″ represents the loss of the stiffness of the viscous damper.
The ratio of the loss of stiffness to the dynamic stiffness is: and the ratio of ,max c d to max d is given as follows [32]: where κ represents the ratio of ,max c d to max d . The applicable condition of Equation (13) is / 5 ≤ , which can fully cover the common range of the parameters of the viscous dampers in practical design. Therefore, Equation (10) can be written as: (12) and the ratio of d c,max to d max is given as follows [32]: Next, when the dynamic stiffness d K is considered, d no longer equals c d . The calculation of d E is same as Equation (5), but max c d in Equation (5) [32]. Before the expression is given, the pseudo loss of stiffness of the viscous damper is defined: where p d K ″ represents the loss of the stiffness of the viscous damper.
The ratio of the loss of stiffness to the dynamic stiffness is: and the ratio of ,max c d to max d is given as follows [32]: where κ represents the ratio of ,max c d to max d . The applicable condition of Equation (13) is / 5 ≤ , which can fully cover the common range of the parameters of the viscous dampers in practical design. Therefore, Equation (10) can be written as: where κ represents the ratio of d c,max to d max . The applicable condition of Equation (13) is Next, when the dynamic stiffness d K is considered, d no longer equals c d . The calculation of d E is same as Equation (5), but max c d in Equation (5) [32]. Before the expression is given, the pseudo loss of stiffness of the viscous damper is defined: where p d K ″ represents the loss of the stiffness of the viscous damper.
The ratio of the loss of stiffness to the dynamic stiffness is: and the ratio of ,max c d to max d is given as follows [32]: where κ represents the ratio of ,max c d to max d . The applicable condition of Equation (13) is / 5 ≤ , which can fully cover the common range of the parameters of the viscous dampers in practical design. Therefore, Equation (10) can be written as: , which can fully cover the common range of the parameters of the viscous dampers in practical design. Therefore, Equation (10) can be written as: The strain energy of the viscous damper can be calculated as half of the product of the maximum damping force F d,max and the maximum deformation d max in one cycle. Consistent with the definition of Equation (4), the EDR of viscous damper ζ d can be approximately calculated as follows: It can be found that ζ d can be determined by ϕ and κ. As κ varies with d max and ϕ is decided by the given α, ζ d is the function of d max . Therefore, d max can reflect the energy dissipation capacity of the viscous damper, and the dampers will be more effective when configured at locations with large story drift. In addition, the limit of story drifts of the damped structure is determined according to the target structural performance. Hence, as the structural performance is determined, the UDR of dampers can be estimated, which will be discussed in a later section.

Basic Concept of Equivalent Damping Ratio for Viscously Damped Structure
Pushover analysis is a common method to evaluate structural seismic performance. According to performance-based seismic design, several pushover methods have been proposed [49,50]. In this study, the CSM was adopted.
As a kind of nonlinear analysis method, the nonlinear behaviors of the structural members should be taken into consideration. When the earthquake drives the structure into an inelastic range, the additional structural damping, which is mainly evoked by the hysteretic effect of the structure members, will be generated. The structural hysteretic damping can also be evaluated by structural EDR, which is a function of the structural ductility factor µ s [51]: where ζ s is the damping ratio of the primary structure and ζ 0 is the inherent damping ratio of the structure. When it affects the energy dissipation structure, the damping induced by the dampers should also be considered. The EDR of the damped structure ζ can be written as: where ζ a is the additional EDR induced by the dampers.
In the CSM, a spectral reduction factor SR is utilized to decrease the elastic spectrum to a reduced response spectrum. The performance point is the intersection point of the reduced response spectrum and the structural capacity spectrum. As the EDR of the damped structure is obtained, the reduction factor SR is given by [50]: where SR A is the spectral reduction value in a constant acceleration range of the spectrum and SR V is spectral reduction value in a constant velocity range of spectrum. The reduced response spectrum is illustrated in Figure 5; T is the natural period of the primary structure and T g is the characteristic period of the site soil.
Buildings 2023, 13, x FOR PEER REVIEW reduced response spectrum is illustrated in Figure 5; T is the natural period of mary structure and g T is the characteristic period of the site soil.

Target Response Mitigation Ratio (RMR)
Since the performance point can be obtained through CSM, the inter-story hereafter referred as drift, can be obtained at this point. Based on the concept o mance-based seismic design, the limits of the drift [ ] θ are assigned according to

Target Response Mitigation Ratio (RMR)
Since the performance point can be obtained through CSM, the inter-story drift, θ i , hereafter referred as drift, can be obtained at this point. Based on the concept of performance-based seismic design, the limits of the drift [θ] are assigned according to different performance levels.
For the structures with additional dampers, the primary structure is supposed to be elastic under frequent earthquakes and exempt from collapse under rare earthquakes. A proper limit for the drift in the design stage should be chosen with respect to the design requirements. The target RMR λ can be determined based on the drift limit [θ], which is shown as follows: The story drifts at the performance point can also be utilized to determine the location of dampers. For all kinds of dampers, it is more effective to install them on stories with large drifts, because that indicates that the story is relatively weak and needs strengthening. Furthermore, minimum deformation is necessary for a viscous damper to dissipate input energy, and with a larger deformation, the dampers can dissipate more input energy. The stories that satisfy the condition of Equation (20) are appropriate for the installation of dampers: where θ denotes the average of story drifts and λ d denotes the limit value of the ratio of story drift to average drift. The value of λ d should be set with consideration of the required story drift, the performance of the dampers, etc.

Determination of UDR
As discussed in Section 2.1, ζ d depends on d max . By obtaining the average d max for the installed dampers, the UDR of the dampers can be realized, which suggests that all dampers can offer uniform viscous damping under seismic excitations. Given that the required structural performance has been determined, the average d max should be calculated according to limit value of the drift [θ]. The equivalent SDOF system of the structure is adopted. The diagram of the SDOF system with a viscous damper is depicted in Figure 6, in which the damper is represented by a Maxwell model. The equivalent story height of the SDOF system can be calculated by Equation (21): where H e denotes the equivalent story height of the SDOF system, which is the massweighted average of story height, m i denotes the mass at level i, H i denotes the height between level i and level i-1, and N is the total number of stories. The displacement can be obtained according to the limit of the drift [θ] by Equation (22): In the SDOF system, the maximum deformation of the damper equals the displacement of the system, namely d max = u max . κ of the viscous damper can be determined on the basis of the limit value of drift [θ], the natural circular frequency of the structure ω, the dynamic stiffness coefficient β, and the damping exponent of the damper α. Then the uniform EDR of dampers ζ d , namely the UDR, can be determined through Equation (15).
where e H denotes the equivalent story height of the SDOF system, which is the mass weighted average of story height, i m denotes the mass at level i, i H denotes the heigh between level i and level i-1, and N is the total number of stories.

Distribution of Damping Force
Given that λ has been determined, the target spectral displacement dt S can be ac quired from Equation (23): where dp S is the spectral displacement at the performance point, at S is the target spec tral acceleration, ap S is the spectral acceleration at the performance point, and t ζ is th EDR of the damped structure at the target spectral displacement.

Distribution of Damping Force
Given that λ has been determined, the target spectral displacement S dt can be acquired from Equation (23): where S dp is the spectral displacement at the performance point, S at is the target spectral acceleration, S ap is the spectral acceleration at the performance point, and ζ t is the EDR of the damped structure at the target spectral displacement. Equation (23) implies that the equivalent period T eq of the structure stays the same after the installation of viscous dampers, and the target spectral displacement S dt is shown in Figure 5. The equivalent period T eq can be obtained by: Then the target additional EDR offered by dampers ζ at can be derived by: where µ st is the structural ductility factor at target spectral displacement. According to the definition of EDR, the additional EDR induced by dampers ζ at can be calculated through Equation (26): where W di is the energy consumed by the viscous dampers at level i, W si is the maximum strain energy of structure at level i, and W sdi is the maximum strain energy of dampers at level i.
By rational design and bracing of the viscous damper, the deformation of viscous damper d i at level i in the horizontal direction can be approximately equal to the inter-story displacement of the structure between level i and level i-1: where ∆ i is the inter-story displacement between level i and level i-1.
As the purpose of adding a viscous damper is to improve the original structural performance, and the structural performance is reflected by story drifts, the allocation of viscous dampers should be optimized with the story drifts. Based on this principle, the design damping forces of viscous dampers on the stories with large drifts can be increased, while the damping forces of viscous dampers on the stories with small drifts can be decreased. The damping force of a viscous damper at level i, F di , is expressed as: where f d is the damping force factor of the viscous dampers. Equation (29) can be obtained by substituting Equations (27) and (28) into Equation (26): where q i is the story shear at level i. Transforming Equation (29), f d can be solved by: As the damping force factor f d is known, the damping force F di of the viscous damper on each story can be obtained through Equation (28). According to F di , the damping coefficient C di and damping exponent α of dampers can be determined.

Uniform Damping Ratio-Oriented Design Procedure
The design flow chart is depicted in Figure 7. The detailed design procedure based on the UDR concept is shown as follows: As the damping force factor d f is known, the damping force di F of the viscous damper on each story can be obtained through Equation (28). According to di F , the damping coefficient di C and damping exponent α of dampers can be determined.

Uniform Damping Ratio-Oriented Design Procedure
The design flow chart is depicted in Figure 7. The detailed design procedure based on the UDR concept is shown as follows: According to the design demand, the parameters of the seismic response spectrum can be determined. By the adoption of CSM, in which the nonlinear behavior of members is considered, the original structural performance can be evaluated.
Corresponding to the performance point obtained by means of CSM, the structural seismic responses can be retrieved, i.e., story drift i θ , story shear i q , and inter-story dis-   According to the design demand, the parameters of the seismic response spectrum can be determined. By the adoption of CSM, in which the nonlinear behavior of members is considered, the original structural performance can be evaluated.
Corresponding to the performance point obtained by means of CSM, the structural seismic responses can be retrieved, i.e., story drift θ i , story shear q i , and inter-story displacement ∆ i . Then, the RMR λ can be calculated according to the target and original structural performance. Considering the utilization efficiency of energy dissipation of viscous dampers, the stories for damper installation can be determined by a reasonable λ d . Therefore, based on the calculated λ, the additional EDR ζ at induced by the viscous damper is established.
The UDR of the viscous damper is determined according to the required seismic performance level. The key parameter, uniform EDR ζ d of the viscous damper, namely UDR, can be estimated by the limit value of drift [θ], the natural circular frequency of the system ω, the dynamic stiffness coefficient β, and the damping exponent of the damper α. As UDR ζ d , structural responses, and supplemental EDR ζ at are known, the damping force factor f d can be solved through Equation (29). After that, the damping force F di at each story is calculated through Equation (28).
The proper viscous damper is selected in accordance with the damping force F di at each story. The damping coefficient C di and damping exponent α of dampers can be determined. The viscous dampers should be installed at the locations in which the dampers can be fully utilized and prevent torsion effects.

Model Illustration
A benchmark six-story RC frame model was adopted to illustrate the design procedure based on the UDR concept. The plan and elevation of the facility are shown in Figure 8. The characteristic period of the site was 0.4 s. According to the Chinese Code of Seismic Design [52], the normalized demand seismic response spectrum is as follows:

Model Illustration
A benchmark six-story RC frame model was adopted to illustrate the design procedure based on the UDR concept. The plan and elevation of the facility are shown in Figure  8. The characteristic period of the site was 0.4 s. According to the Chinese Code of Seismic Design [52], the normalized demand seismic response spectrum is as follows:  The story height and mass are shown in Table 2. The cross sections of the members are listed in Table 3. The members are made of concrete with a compression strength of 30 MPa according to the Code for Design of Concrete Structures [52]. The strength of the longitudinal and stirrup rebars was 345 MPa. The model of the structure was established using the software SAP2000 [53]. The beams and columns were modeled as elastic bars The story height and mass are shown in Table 2. The cross sections of the members are listed in Table 3. The members are made of concrete with a compression strength of 30 MPa according to the Code for Design of Concrete Structures [52]. The strength of the longitudinal and stirrup rebars was 345 MPa. The model of the structure was established using the software SAP2000 [53]. The beams and columns were modeled as elastic bars with two inelastic hinges at both ends. The moment hinge was applied for the beams, while the P-M-M hinge was applied for the columns. The moment-rotation relationship was obtained with the sTakeda model [54]. The thin shell element was utilized to simulate the floors and roof, and the assumption of rigid diaphragms for the floors was adopted. The distribution loads were exerted on the floors. The dead load was 5.5 kN/m 2 for every story. The live load was 3.0 kN/m 2 for stories 1-5 and 5.0 kN/m 2 for story 6. The selfweights of the floors were included in the dead load. The story mass was calculated by (1.0 × Dead load + 0.5 × Live load)/g. The inherent damping ratio of the structure was 0.05 and the natural period of the structure was 1.266 s.  Table 3. Cross section of members.

1-2 3-4 5-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 mm The stiffness of the structure in the Y direction was weak, so viscous dampers were required to be installed in this direction. For simplicity, only the seismic responses in the Y direction were analyzed.

Detailed Design Procedure
The pushover analysis was carried out on the primary structure and the force-displacement curve of the structure was recorded. After that, the force-displacement curve and the design response spectrum were transferred to an acceleration-displacement response spectrum.
The limit values of story drifts [θ] under the frequent, moderate, and rare earthquake scenarios adopted in this example are listed in Table 4. The performance point of the structure was obtained by CSM under a moderate seismic level. The responses of each story at the performance point were retrieved, which are plotted in Figure 9 and listed in Table 5. The limit value of story drifts [θ] was set to 0.33% and the value of RMR λ was 65.1%. In this case, the value of λ d was set to 0.5, which indicated that all stories required the installation of viscous dampers. Table 4. Limit value of story drifts under different seismic levels.

Seismic Levels
[ ] θ Frequent earthquake 0.125% Moderate earthquake 0.33% Rare earthquake 0.67%  According to λ, the target point (S dt , S at ) of the structure was determined. The equivalent period T eq at the target point was 1.5 s. Table 6 lists all the calculation parameters that are needed for Equation (25). Next, the deformation coefficient of the dampers was calculated. By means of modal analysis, the natural circular frequency of the structure was found to be 4.96 rad/s. In this project, the damping exponent α and the dynamic stiffness coefficient β of the dampers were set to be 0.2 and 1 (mm 0.8 s 0.2 ), respectively. Because the limit value of drift [θ], the natural circular frequency ω of the structure, the damping exponent α and the dynamic stiffness coefficient K d were known, κ could be determined through (13), which was 0.884. Therefore, the UDR of viscous damper ζ d was 0.53.
The damping force factor f d was derived by substituting ζ at and story responses q i and θ i into Equation (29). As a result, f d was determined to be 1135 kN. Therefore, the damping force of the dampers at each story was determined. The damping coefficient C di of the viscous damper at each story an be calculated according to Equation (32): (32) where v i is the relative velocity at level i. The velocity at each story was taken to be 100 mm/s based on engineering experience. Table 7 shows the design parameters of the viscous damper. On each story, two viscous dampers were added along axis X3 and X9, respectively. The diagram of configurations along X3 is depicted in Figure 10. viscous dampers 1~6 represent the viscous dampers installed on the corresponding stories.
Buildings 2023, 13, x FOR PEER REVIEW On each story, two viscous dampers were added along axis X3 and X9, resp The diagram of configurations along X3 is depicted in Figure 10. viscous damp represent the viscous dampers installed on the corresponding stories.

Seismic Performance Verification
A nonlinear time history analysis was carried out to validate the mitigation e the seismic response of the structure with the viscous dampers designed by the p method. The response of the damped structure (ST1) was compared with thos structure without dampers (ST0). According to Chinese code [52], the site chara period for the soil is 0.4 s, and the acceleration response spectrum of the input

Seismic Performance Verification
A nonlinear time history analysis was carried out to validate the mitigation effect on the seismic response of the structure with the viscous dampers designed by the proposed method. The response of the damped structure (ST1) was compared with those of the structure without dampers (ST0). According to Chinese code [52], the site characteristic period for the soil is 0.4 s, and the acceleration response spectrum of the input seismic wave can be determined accordingly. Based on the spectrum, seven earthquake waves were selected, five of which were natural seismic records. These natural seismic records were selected from the Pacific Earthquake Engineering Research Center ground motion database [55]. The spectra of the selected waves were matched with the design acceleration spectrum. Detailed information of the natural seismic records is described in Table 8. The normalized spectra of all the waves are shown in Figure 11, in which AW1 and AW2 represent the two artificial records. Peak ground acceleration (PGA) was 0.7 m/s 2 , 2 m/s 2 , and 4 m/s 2 corresponding to the frequent, moderate, and rare earthquake scenarios, respectively, according to the specifications in Chinese code [52].  The story drifts of ST1 and ST0 were compared to illustrate the mitigation effects of the designed viscous dampers. Figures 12 and 13 show a comparison of average story drifts and shears under frequent, moderate, and rare earthquakes. The story drifts were controlled within the limit values and the shears were simultaneously reduced.  The story drifts of ST1 and ST0 were compared to illustrate the mitigation effects of the designed viscous dampers. Figures 12 and 13 show a comparison of average story drifts and shears under frequent, moderate, and rare earthquakes. The story drifts were controlled within the limit values and the shears were simultaneously reduced. The story drifts of ST1 and ST0 were compared to illustrate the mitigation effects of the designed viscous dampers. Figures 12 and 13 show a comparison of average story drifts and shears under frequent, moderate, and rare earthquakes. The story drifts were controlled within the limit values and the shears were simultaneously reduced.   Viscous dampers can dissipate great input energy under different earthquake scenarios. The hysteresis curves of viscous dampers installed on the first, third, and fifth stories under AW1 for PGA = 200 gal are displayed in Figure 14, which indicates that the capacity of energy dissipation of the dampers was fully utilized. It should be noted that in practical design, the parameters of viscous dampers need not be strictly identical to the design values. Generally, for the purpose of convenience and efficiency, different viscous dampers with similar parameters are merged into one type.

Comparison with Classic Design Method and Discussion
To illustrate the advantages of the design method based on the UDR concept, it was  Table 9 shows the average EDR of dampers ζ d of ST1 under different seismic excitations. Because the story drift is well controlled within the limit value of story drift [θ], the average value of ζ d under the earthquakes with PGA = 200 gal is close to and slightly less than the anticipated UDR 0.53. In addition, it can be observed that ζ d under different earthquakes with the same PGA are close to each other, which verifies the rationality of the UDR concept in a different way.

Comparison with Classic Design Method and Discussion
To illustrate the advantages of the design method based on the UDR concept, it was compared to the classic method in the JSSI Manual [32]. The JSSI method adopts the elastic response spectrum, and the parameters of viscous dampers are determined according to lateral stiffness at each story.
The limit values of story drifts under different seismic levels for the JSSI design method were set to be the same as those for the UDR design method. The expected damping force and the damping coefficient of the viscous damper at each story obtained by the JSSI method are shown in Table 10. The seismic responses of the structure designed by the JSSI method under the same earthquake waves were compared with those of the structure designed by the UDR method. Figure 15 represents a comparison of story drifts under frequent, moderate, and rare earthquakes, respectively. From Figure 15, it can be observed that the seismic responses of the two damped structures designed by the two different methods are similar. The seismic responses of the structure designed by the JSSI method under the same earthquake waves were compared with those of the structure designed by the UDR method. Figure 15 represents a comparison of story drifts under frequent, moderate, and rare earthquakes, respectively. From Figure 15, it can be observed that the seismic responses of the two damped structures designed by the two different methods are similar.  Figure 16 shows the comparison of hysteresis curves of viscous dampers installed on the first story of the two damped structures under AW1 with PGA = 200 gal and NW3 with PGA = 400 gal. By reducing the damping force of the viscous dampers installed on the first story, the dampers of the structure designed by the UDR method obtain larger deformation, and the energy dissipation capacity of these dampers is more efficiently utilized. Moreover, the redundant damping forces of the dampers in the damped structure designed by the JSSI method can increase the forces in the members of the primary structure and lower the utilization efficiency of the energy dissipation capacity of the dampers.  Figure 16 shows the comparison of hysteresis curves of viscous dampers installed on the first story of the two damped structures under AW1 with PGA = 200 gal and NW3 with PGA = 400 gal. By reducing the damping force of the viscous dampers installed on the first story, the dampers of the structure designed by the UDR method obtain larger deformation, and the energy dissipation capacity of these dampers is more efficiently utilized. Moreover, the redundant damping forces of the dampers in the damped structure designed by the JSSI method can increase the forces in the members of the primary structure and lower the utilization efficiency of the energy dissipation capacity of the dampers. with PGA = 400 gal. By reducing the damping force of the viscous dampers installed o the first story, the dampers of the structure designed by the UDR method obtain larg deformation, and the energy dissipation capacity of these dampers is more efficiently u lized. Moreover, the redundant damping forces of the dampers in the damped structu designed by the JSSI method can increase the forces in the members of the primary stru ture and lower the utilization efficiency of the energy dissipation capacity of the damper