Engineering, Technology & Applied Science Research Vol. 9, No. 1, 2019, 3796-3802 3796 www.etasr.com Xu & Yamanari: Performance of Steel Frame with Linkage System under Earthquake Excitation Performance of Steel Frame with Linkage System under Earthquake Excitation Guang Xu Department of Advanced Industrial Science Kumamoto University Kumamoto, Japan xuguang1991@hotmail.com Minoru Yamanari Faculty of Advanced Science and Technology Kumamoto University Kumamoto, Japan yamanari@kumamoto-u.ac.jp Abstract—Authors have already proposed a linkage system with rotational friction dampers that can be expected to dissipate seismic energy rather than the main frame. In this paper, the relationship between the rotational angle of the linkage system and inter-story drift angle of a structure was investigated. Consequently, a three-bay twelve-story steel frame with linkage system was analyzed by a numerical simulation program based on a theory of non-linear dynamic analysis. The distribution of slipping moment of friction damper along with the height of the main structure was discussed. Meanwhile, optimal slipping moment range was suggested according to analysis results. Keywords-steel frame; linkage system; friction damper I. INTRODUCTION There could be a big increase in the number of tremendous earthquakes around the world in the following years [1]. So far, it is difficult to predict where and when these earthquakes will occur. Therefore, in order to protect people and properties, a number of methods to decrease structure vibration under earthquake excitation have been proposed. These methods/control systems could be divided into four patterns: passive, active, semi-active and hybrid [2]. A passive control system works without requiring an external power source, which is necessary for active and semi-active control system. Generally, passive control system consists of one or more equipment designed to modify the structural properties like ductility and stiffness, causing a reduction in the structure seismic behavior. Friction damper is a typical representative of passive control systems, which is used to dissipate the seismic input energy via solid friction. This type of device was pioneered in 1982 [3]. It is similar to automobile brakes, which use the mechanism of solid friction to dissipate seismic energy reducing the motion of the structure. Commonly, the friction damper is made of a series of special steel plates that are bolted by a high strength screw. Suitable friction pads can be clamped between these steel plates, in order to create frictional performance. Because the clamping force is controlled by a screw, it is very convenient to change the frictional force to satisfy diffident purposes. However, a lower limit value of frictional force is also necessary, as it is not expected for friction damper to slip during a severe wind or earthquake. When the force or moment developing in the friction damper reaches a predetermined value, it will start slipping and keeping a steady status. That is a remarkable feature differing from other dampers, liking viscous damper, oil damper and TMD. According to the connection with the main structure, friction dampers can be divided into four types: Pall friction damper, Sumitomo friction damper, gap-screw connection damper, and rotational friction damper [4]. Pall friction damper has two friction joints on the upper and down arm-brace, respectively. Four cross braces are used to connect it with the main structure, which are pin connected at both ends. It has been widely used to strengthen the seismic behavior of many structures around the world [2, 5-10]. It is easily found that the Pall friction damper has only translational movement at the friction joints when it starts slipping [11]. Differing from it, in 2000 the DAMPTECH company presented a rotational friction damper in which the movements are rotational [12]. A common point of all dumpers is that they will dissipate seismic energy to protect structures from damage or collapse. In order to extend the dissipation energy to a maximum value, some researches about optimal friction force or moment have been done [2, 12- 15]. However, the optimal value of the friction damper in each story is the same with the one of the first story, but, seismic shear force gradually decreases along with the height of the structure. It means that if a same optimal value was used for every damper of each story, it will not get the best seismic performance of the structure because the top dampers will not slip or slip just a little bit. For this reason, the distribution of the optimal value along the height of the structure will be discussed in this paper. On the other hand, the capacity of energy dissipation within a friction damper almost depends on slip displacement and slip force. For the same slip force, if the friction damper slipped on a maximum displacement at each step, the amount of dissipated energy will be increased. Hence, the objectives of this paper are: 1) To propose a linkage system, which will enlarge the slip displacement of friction damper. 2) To discuss the distribution of slipping moment in dampers along the structure height. 3) To find the optimal value for which the best seismic performance of the structure occurs. II. LINKAGE SYSTEM As mentioned above, the capacity of energy dissipation in the friction damper almost depends on the slip displacement Corresponding author: G. Xu Engineering, Technology & Applied Science Research Vol. 9, No. 1, 2019, 3796-3802 3797 www.etasr.com Xu & Yamanari: Performance of Steel Frame with Linkage System under Earthquake Excitation and slip force. The proposed mechanism of displacement amplification is shown in Figure 1. It can be regarded as a combination of Pall friction damper and rotational friction damper. The friction joints, which are placed in the center of the side of the quadrangle in Pall friction damper, are removed to its corner. The movement in the linkage system is rotational not translational as in the Pall friction damper. When a tension or compression acts on the corner of the quadrangle, it will be shaped as a parallelogram, and the relative rotational angle can be obtained from the four corner joints. That is the slip displacement of the rotational friction damper. The connections of the other elements to the main structure are the same as in the Pall friction damper. In the linkage system, cross brace is pinned to the main structure and the other end of it is directly connected to the damper. Because the arm braces rotate on the contrary direction, there will be no moment generated in the cross braces like in the ones in Pall friction damper. Therefore, comparing to Pall friction damper, it has a simplified construction in which the number of connection bolts, slip bolts and curve slots has been reduced. Furthermore, instead of translational displacement, rotational displacement is used for energy dissipation. rotational friction damper arm brace cross brace pin connection force tension force bending moment H L R a b Fig. 1. The conception of linkage system a b c A BO φ1 φ2 A’ B’O’ c+δ δα αφ1-β φ2-α β Fig. 2. Geometrical relationshiop III. GEOMETRICAL RELATIONSHIP WITH MAIN STRUCTURE To identify the mechanism of displacement amplification of the linkage system, the relationship between the rotational angle in the linkage system and the drift angle in the main structure needs to be formulated. In this section, it is assumed that deformation in elements will not be considered. The computation schematic for the half of linkage system is shown in Figure 2, where a and b are the span and height of the linkage system, c is the length of the diagonal line, φ1 and φ2 are the inner angles of triangle AOB, α and β are the rotational angles, δ is the increment of diagonal line. Using cosine law, the relationship of a and δ may be written as: 2 2 2 1 ( ) cos( ) 2 ( ) a c b a c           (1) Using addition theorem ,(1) could be written as: 2 2 2 2 ( )cos sin ( )a c bd e a c        (2) where, d and e are expressed as: d = cosϕ 1 = a2 + c2 − b2 2ac (3) e = sinϕ 1 = 1− d2 (4) Taylor series expansion is used to simplify the rotational angle α: 2 2 2 21 2 2 ( )( )a c bd e d a c         (5) According to (5), α expressed by δ could be earned as: 2 1             ( ) ( )( )ace ae acd ad xce ecd d acd ad (6) in which, x1 is expressed as: 2 2 2 2 1 2 2 2( )x a b c acd c ad         (7) Using the same method, the other rotational angle β could be expressed as: 2 2 ( ) ( )( )bcg bg bcf bf xcg gcf f bcf bf        (8) where, x2 is: 2 2 2 2 2 2 2 2( )x a b c bcf c bf          (9) Then, the relative rotational angle of damper could be written as:     (10) Using (1)-(10), the relationship between the damper relative rotational angle and the increment of quadrangle diagonal line could be obtained. In addition, the relationship between the drift angle of main structure and the increment of diagonal line of quadrangle is necessary. It is assumed that main structure deforms as shown in Figure 3, where H is the height of main structure, L is the length of the span, and R is the drift angle. Therefore, the increment of diagonal line of quadrangle could be expressed as: 2 2 2 2( )L HR H L       (11) Using (10)-(11), the relationship between the damper’s relative rotational angle and the main structure’s drift angle can be obtained. To examine the accuracy of equations, numerical analysis is necessary. Program OpenSees based on the finite element method is used to simulate the numerical example. The model shown in Figure 1 is used for the analysis, where the lengths of arm braces a and b are 1.4m and 0.7m respectively, and height and span of structure are 3.5m and 7m respectively. The beam is assumed infinitely rigid and the columns are fixed at their bases. The steel columns are wide-flange sections with the moment of inertia of 3.4×10-5m4. The sectional area of the brace is 4×10-3m2 and the slipping moment of damper is set as Engineering, Technology & Applied Science Research Vol. 9, No. 1, 2019, 3796-3802 3798 www.etasr.com Xu & Yamanari: Performance of Steel Frame with Linkage System under Earthquake Excitation 20Nm. The dead load is assumed as 392×103N and the critical damper ratio is 5%. Newmark-β (β=0.25) method is adopted for the numerical integration. The comparison between the approximate value obtained by the equations and the geometrical exact value obtained by OpenSees is shown in Figure 4. The results indicated that the accuracy of equations is very high in normal range (i.e. R≤0.06). R H L HR Fig. 3. Structure deformation 0 0.06 0.12 0 0.4 0.8 1.2 1.6 exact approximateD ri ft , R ( m m ) Relative rotational angle, φ (rad) Fig. 4. Accuracy of equations IV. FINITE ELEMENT ANALYSIS A. Finite Element Modeling For discussing the distribution of friction force in the damper along the height of main structure, two models are created by SAP2000, as shown in Figure 5. Their difference will be introduced below. Both are the models of a 12-story two-dimensional structure with 7m span lengths and a height of 3.5m, with 3 openings. Supports are assumed to be fixed and the soil structure interaction effects are neglected. Dead and live loads are concentrated mass at the beam-to-column node. The sections of the main structure, shown in Table I, are made from SN490 (based on JIS standard). The plastic behavior of elements in main structure is considered by simple hinges based on TABLES 9-6 IN ASCE 41-13 [16], which are assigned to both ends of the elements. The cross braces and arm members should be considered strong enough to remain elastic during earthquakes. These elements are selected from circular cross section cables and beam element, respectively. Meanwhile, the ratio of initial stiffness of the linkage system to one of the main structures is set as 1.0. In order to modify friction damper in linkage system, plastic-Wen element is used in SAP2000. The friction behavior is considered in R3 direction and the other directions are fixed. A yielding exponent of 10 indicating a sharp transition from linear to plastic phase and a post yield stiffness ratio of 0.0001 indicating a rectangular hysteresis loop are used to simulate the friction behavior of dampers [14]. The details of the damper model are shown in Figure 6. In addition, P-Δ effects are included to account for the geometrical nonlinearity effects and out-of-plane distortion of all elements is not considered. 4 2 m 21m cross brace rotational friction damper pin connection arm brace L=7m H = 3 .5 ma=1.4m b = 0 .7 m :mass=20,000 (Kg) Fig. 5. Finite element model TABLE I. THE SECTION USED IN MAIN STRUCTURE Story No. Column Beam 10-12 500✕500✕19 450✕250✕12✕28 7-9 500✕500✕22 500✕300✕16✕25 4-6 500✕500✕25 550✕300✕16✕28 1-3 500✕500✕28 600✕250✕16✕32 Ms -Ms 0 φ k 0.0001k Bar element Beam element B ea m e le m en t Plastic-Wen element Fig. 6. Damper model and hysteretic curve B. Analysis Conditions and Methods In this analysis, the damping of the model is considered as a Rayleigh damping. First-order and second-order damping coefficients are both 2%. Newmark-β method with β=0.25, α=0.5 is used to analyze the model. The output time step size is set at 0.001s and seismic duration is 40s. Four observed seismic waves, shown in Table II, are used in nonlinear time history analysis. In order to observe the seismic behavior of the structure equipped with a linkage system under a tremendous earthquake, the peak ground velocity (PGV) of each wave is scaled to 1.0m/s. The input direction of seismic waves is horizontal. For determining the slipping moment of friction damper, slipping moment ratio As is defined as: 1 2 1 si sis si si M MA M M  (12) where Msi is the slipping moment of friction damper in the i-th story, Msi1 is a moment, which is generated in friction damper of the i-th story when a group of forces (base shear force=0.2g) according to Ai distribution [17] were imposed on the structure and friction damper did not slip. In the same way, when the forces equating to the held horizontal yield strength of structure Engineering, Technology & Applied Science Research Vol. 9, No. 1, 2019, 3796-3802 3799 www.etasr.com Xu & Yamanari: Performance of Steel Frame with Linkage System under Earthquake Excitation [17] were imposed on the structure, the generated moment in i- th story is determined as Msi2. In this analysis, Msi1 is set as the lower limit value of slipping moment in friction dampers and Msi2 is set as the upper limit value. Hence, the friction dampers will not slip under a slight earthquake and start to slip before the ultimate limit state of main structure. By changing As from 0.0 to 1.0, Msi basing on Ai distribution could be obtained by (12). From [18] we can see that a 5000KN capacity of damper has been designed and experimented successfully. By comparison with Msi in this paper, the slipping moment is much smaller than that and can be considered to be achieved by controlling the axial force of bolts. However, it is worth noting that the above-mentioned distribution of slipping moment is only for Model-A. For Model-B, the slipping moment of damper in each story is the same as that of the damper in the 12th story so that the damper of each story can easily slip under earthquake. C. Model Analysis Main and damped frame periods are shown in Table III. The periods of the damped frame were obtained by rigid joints instead of damper connection (i.e., plastic-wen element). TABLE II. SEISMIC WAVES Earthquake PGV (m/s) PGA (m/s2) Time (s) Elcentro NS 1.0 10.21 40 Tohoku EW 7.34 Hachinohe EW 5.10 Kumamoto EW 8.02 TABLE III. STRUCTURE PERIODS Periods Damped structure (s) Main structure (s) 1st 1.06 1.56 2nd 0.35 0.55 D. Comparison of the Two Models In order to observe the seismic behavior of two models that used different distribution of slipping moment, drift, bending moment, shear force and energy dissipation of friction dampers are used as reference indexes for a fixed value of As (As=0.0). Drifts of the two models under the four seismic waves are shown in Figure 7. It can be observed that each story drifts in Model-A are smaller than the ones of each story in Model-B under any seismic wave. However, according to various seismic waves, the drift reduction ratio is different. In the case of El centro NS, the maximum and minimum drift reduction ratios are 37% (5th story) and 4% (12th story), respectively. In another case of Kumamoto EW, maximum and minimum drift reduction ratios dwindle to 17% (1st story) and 2% (5th story), respectively. Despite all this, it still can be considered that for drift reduction, the distribution of slipping moment in Model-A is much more efficient than that in Model-B. A comparison of the bending moment of the beam element in the two models is shown in Figure 8. The beam bending moment of Model-A is smaller than that of Model-B. It means the importance of the distribution of slipping moment. In the four cases, the bending moment reduction of the sub-structure is relatively more obvious than that of the upper structure. (a) 0 0.005 0.01 0.015 Drift (rad) 2 4 6 8 10 12 S to ry N o . Model-A Model-B (a) (b) 0 0.005 0.01 0.015 Drift (rad) 2 4 6 8 10 12 S to ry N o . (b) (c) 0 0.005 0.01 0.015 Drift (rad) 2 4 6 8 10 12 S to ry N o . (c) (d) 0 0.005 0.01 0.015 Drift (rad) 2 4 6 8 10 12 S to ry N o . (d) Fig. 7. Drift of each story in the two models: (a) ElcentroNS, (b) HachinoheEW, (c) KumamotoEW, (d) TohokuEW. Because shear force of the sub-structure is larger than that of the upper structure, the bending moment of friction damper installed in the sub-structure is easy to reach the predetermined value, and then the dampers can dissipate more seismic energy. It also means that energy dissipation of friction damper installed in the sub-structure is more than that in the upper structure. Although in the case of Kumamoto EW, a part of bending moment in Model-A is almost the same with that in Model-B, as a whole, the distribution of slipping moment in Model-A is more efficient. The shear force of each story in the two models is shown in Figure 9. The reduction of shear force in Model-A exhibits the superiority of Model-A. The maximum reduction ratios in four cases are 37% (Elcentro NS, 5th story), 22% (Hachinohe EW, 2nd story), 18% (Kumamoto EW, 10th story), 26% (Tohoku EW, 9th story), respectively. To compare the capacity of energy dissipation of friction damper in the two models, the percentage of total input energy dissipated by the Engineering, Technology & Applied Science Research Vol. 9, No. 1, 2019, 3796-3802 3800 www.etasr.com Xu & Yamanari: Performance of Steel Frame with Linkage System under Earthquake Excitation dampers is shown in Figure 10 As we see the dissipated energy by the dampers in both models has exceeded the 50% of the total input energy. This means that the friction dampers have a large potential in energy dissipation. The dissipated energy in Model-A is more than that in Model-B. In the Model-A, the maximum value is about 67% in the case of Hachinohe EW. From Figures 7-10 can be concluded that the distribution of the slipping moment in dampers along the height of Model-A is more efficient. (a) 0 500 1000 1500 2000 2500 Bending moment (KNm) 2 4 6 8 10 12 S to ry N o . Model-A Model-B (a) (b) 0 500 1000 1500 2000 2500 Bending moment (kKNm) 2 4 6 8 10 12 S to ry N o. (b) (c) 0 500 1000 1500 2000 2500 Bending moment (KNm) 2 4 6 8 10 12 S to ry N o. (c) (d) 0 500 1000 1500 2000 2500 Bending moment (KNm) 2 4 6 8 10 12 S to ry N o. (d) Fig. 8. Bending moment of beam in each story: (a) ElcentroNS, (b) HachinoheEW, (c) KumamotoEW, (d) TohokuEW. E. Optimal Value of Slipping Moment The optimal value of the slipping moment in Model-A will be discussed in this section. The normalized response values of the structure should be studied as a function of the slipping moment ratio (As). The normalized drift is described by (13) in which RiAs=j is the drift in i-th story when As=j and R i As=0 is the drift in i-th story when As=0. 0 0 0 1 step 0 1 1 12 step 1, , , . , ,s s iAi jratio jA RR j iR          (13) (a) 0 1000 2000 3000 4000 Shear force (KN) 2 4 6 8 10 12 S to ry N o. Model-A Model-B (a) (b) 0 1000 2000 3000 4000 Shear force (KN) 2 4 6 8 10 12 S to ry N o. (b) (c) 0 1000 2000 3000 4000 Shear force (KN) 2 4 6 8 10 12 S to ry N o. (c) (d) 0 1000 2000 3000 4000 Shear force (KN) 2 4 6 8 10 12 S to ry N o. (d) Fig. 9. Shear force of each story: (a) ElcentroNS, (b) HachinoheEW, (c) KumamotoEW, (d) TohokuEW. Elce ntro NS Hac hino heE W Kum amo toEW Toh oku EW 50% 55% 60% 65% 70% P er ce n ta ge ( % ) Model-A Model-B Fig. 10. Percentage of total input energy dissipated by damper Similarly, the normalized bending moment of beam is presented by: Engineering, Technology & Applied Science Research Vol. 9, No. 1, 2019, 3796-3802 3801 www.etasr.com Xu & Yamanari: Performance of Steel Frame with Linkage System under Earthquake Excitation 0 0 1 step 0 1 1 2 step 1, , , . , ,s s iA ji jratio iA BMBM j iBM          (14) From Figure 11, the drift of each story expresses a different variation trend with the increasing of As. In the four cases, the drift of the sub-structure (white marking) increased with the increasing of As. However, in the cases of Elcentro NS and Kumamoto EW, the drift of the upper structure (black marking) nearly kept constant. In the case of Hachinohe EW, it kept decreasing. In the case of Tohoku EW, it firstly decreased and then increased for As within the range from 0.0 to 0.6. (a) 0 0.2 0.4 0.6 0.8 1 Slipping moment ratio 0.8 0.9 1 1.1 1.2 N o rm al iz ed d ri ft (a) (b) 0 0.2 0.4 0.6 0.8 1 Slipping moment ratio 0.8 0.9 1 1.1 1.2 N o rm al iz ed d ri ft (b) (c) 0 0.2 0.4 0.6 0.8 1 Slipping moment ratio 0.8 0.9 1 1.1 1.2 N o rm al iz ed d ri ft (c) (d) 0 0.2 0.4 0.6 0.8 1 Slipping moment ratio 0.8 0.9 1 1.1 1.2 N o rm al iz ed d ri ft (d) Fig. 11. Normalized drift of each story in Model-A: (a) ElcentroNS, (b) HachinoheEW, (c) KumamotoEW, (d) TohokuEW. The normalized beam bending moment shown in Figure 12 has a similar trend with the drift. However, in the case of Tohoku EW, the minimum bending moment of the upper structure can be obtained when As=0.4. Comparing the two Figures, it can be observed that for the dampers installed in the sub-structure, the smaller the As the better. For those in the upper structure, the range of As from 0.0 to 0.4 is suggested. Generally, the percentage of the total input energy dissipated by the dampers is an important index to determine an optimal value of slipping moment (Figure 13). As shown in Figure 13, more than 60% of total input energy has been dissipated by friction dampers, which is a very large value and obviously not suitable. The peak of the dissipated energy (black marking) is obtained for As within the range from 0.0 to 0.4. Sometimes, the dissipated energy does not necessarily indicate the best performance of the structure, as it is accompanied by increase in the total seismic energy [12]. (a) 0 0.2 0.4 0.6 0.8 1 Slipping moment ratio 0.9 0.95 1 1.05 1.1 1.15 1.2 N o rm al iz ed b en d in g m o m en t (a) (b) 0 0.2 0.4 0.6 0.8 1 Slipping moment ratio 0.9 0.95 1 1.05 1.1 1.15 1.2 N o rm al iz ed b en d in g m o m en t (b) (c) 0 0.2 0.4 0.6 0.8 1 Slipping moment ratio 0.9 0.95 1 1.05 1.1 1.15 1.2 N o rm al iz ed b en d in g m o m en t (c) (d) 0 0.2 0.4 0.6 0.8 1 Slipping moment ratio 0.9 0.95 1 1.05 1.1 1.15 1.2 N o rm al iz ed b en d in g m o m en t (d) Fig. 12. Normalized bending moment of each story in Model-A: (a) ElcentroNS, (b) HachinoheEW, (c) KumamotoEW, (d) TohokuEW. Engineering, Technology & Applied Science Research Vol. 9, No. 1, 2019, 3796-3802 3802 www.etasr.com Xu & Yamanari: Performance of Steel Frame with Linkage System under Earthquake Excitation Fig. 13. The effect of slipping moment ratio on energy dissipation V. CONCLUSIONS Basing on the feature of friction damper, (i.e., displacement dependence) a concept of linkage system that could be considered as a combination of the Pall friction damper and the rotational friction damper has been presented. It can improve the capacity of energy dissipation of friction damper by increasing the rotational angle. The relationship between rotational angle and the main structure drift has been derived as a formula by trigonometric functions and the high accuracy of the formula has been verified by numerical analysis. Two models having different distribution of slipping moment along the structure height have been analyzed and compared. Analysis results demonstrated that a properly distribution of the slipping moment is necessary and the Model-A is more efficient than Model-B. In addition, the optimal value of slipping moment, (i.e. As) in Model-A has been discussed. Time history analysis results showed that for the friction damper installed in the sub-structure, the smaller the value of As the better and for the upper ones, a range of As from 0.0 to 0.4 is suggested. Further research is needed to verify the above conclusions. In future studies, stiffness of linkage system, intensity of seismic waves, distribution of linkage system, etc. should be factored in and the construction of linkage system should be discussed and tested by a series of experimental and numerical analyses. REFERENCES [1] United States Geological Survey, USGS Real-time Earthquake , available at https://www.usgs.gov/products/data-and-tools/real-time- data/earthquakes [2] S. M. Zahraei, A. Moradi, M. Moradi, “Using pall friction dampers for seismic retrofit of a 4-story steel building in Iran”, Topics in Dynamics of Civil Structures, Vol. 4, pp. 101-107, 2013 [3] A. S. Pall, C. Marsh, “Response of Friction Damped Braced Frames”, Journal of Structural Engineering, Vol. 108, No. 9, pp. 1313-1323, 1982 [4] V. Ruiz, K. Walsh, M. M. Abdullah, Investigating The Energy Dissipating Properties of Passive Friction Devices, Florida State University, &Tokyo Institute of Technology, Tokyo, Japan, 2005 [5] R. Chandra, M. Masand, S. K. Nandi, C. P. Tripathi, R. Pall, A. Pall, “Friction-dampers for seismic control of La Gardenia towers south city, Gurgaon, India”, 12th World Conference on Earthquake Engineering, New Zealand, January 30-February 4, 2000 [6] T. H. Hale, R. Pall, “Seismic upgrade of the Freeport water reservoir, Sacramento, California”, 12th World Conference on Earthquake Engineering, New Zealand, January 30-February 4, 2000 [7] A. S. Pall, R. Pall, “Friction-damper used for seismic control of new and existing buildings in Canada”, ATC-17-1 Seminar on Seismic Isolation, Passive Energy Dissipation and Active Control, San Francisco, USA, March 11-12, 1992 [8] B. Wu, H. Li, L. Lin, M. Shan, “Seismic retrofit of a city hall in Northeast China with frictional energy dissipators”, Journal of Building Structures, Vol. 19, No. 5, pp. 28-36, 1998 (in Chinese) [9] A. Pall, S. Vezina, P. Proulx, R. Pall, “Friction‐Dampers for Seismic Control of Canadian Space Agency Headquarters”, Earthquake Spectra, Vol. 9, No. 3, pp. 547-557, 1993 [10] T. H. Hale, R. T. Pall, “Seismic upgrade of the freeport water reservoir, Sacramento, California”, 12th World Conference on Earthquake Engineering, New Zealand, January 30-February 4, 2000 [11] D. Shao, A. Pall, B. Soli, “Friction Dampers for Seismic Upgrade of A 14-Story Patient Tower With A 36-Foot Tall Soft-Story”, 8th US National Conference on Earthquake Engineering, San Francisco, USA, April 18-22, 2006 [12] H. I. Mualla, B. Borislav, “Performance of steel frames with a new friction damper device under earthquake excitation”, Engineering Structure, Vol. 24, pp. 365-371, 2002 [13] K. Jafarzadeh, M. L. Yaghin, R Sabetahd, “Evaluation pf Pall Friction Damper Performance in Near-Fault Earthquakes by Using of Nonlinear History Analysis”, World Applied Sciences Journal, Vol. 20, No. 2, pp. 264-270, 2012 [14] G. S. Adithya, H. Narendra, “Performance evaluation of friction dampers under seismic loads”, International Journal of Research in Engineering and Technology, Vol. 5, No. 14, pp. 10-15, 2016 [15] P. Hosseini Sarjou, N. Shabakhty, “Effect of the Improved Pall Friction Damper on the Seismic Response of Steel Frames”, Engineering Technology & Applied Science Research, Vol. 7, No. 4, pp. 1833-1837, 2017 [16] American Society of Civil Engineers, Seismic Evaluation and Upgrade of Existing Buildings, 2014 [17] Architectural Steel Structure - Theory and Design, Kajima Institute Publishing, 2007 (in Japanese) [18] I. H. Mualla, E. D. Jakupsson, L. O. Nielsen, “Structure behavior of 5000KN damper”, 14th European Conference on Earthquake Engineering, Republic of Macedonia, August 30-September 3, 2010