Next Article in Journal
A Study on Minimizing Welding Deformation of Joints for the Sealing of Emission After-Treatment Structure
Next Article in Special Issue
Experimental Investigation on Bending Behavior of Existing RC Beam Retrofitted with SMA-ECC Composites Materials
Previous Article in Journal
A Novel Additive Manufacturing Method of Cellulose Gel
Previous Article in Special Issue
Seismic Response of Resilient Bridges with SMA-Based Rocking ECC-Reinforced Piers
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Seismic Response Mitigation of a Television Transmission Tower by Shape Memory Alloy Dampers

1
Key Laboratory of Roadway Bridge and Structural Engineering, Wuhan University of Technology, Wuhan 430070, China
2
School of Civil Engineering, Hefei University of Technology, Hefei 230009, China
*
Author to whom correspondence should be addressed.
Submission received: 19 August 2021 / Revised: 7 November 2021 / Accepted: 9 November 2021 / Published: 18 November 2021
(This article belongs to the Special Issue Shape Memory Alloys for Civil Engineering)

Abstract

:
High-rise television transmission towers are of low damping and may vibrate excessively when subjected to strong earthquakes. Various dynamic absorbers and dampers are proposed to protect television transmission towers from excessive vibrations and damages. Up to now, the seismic damage reduction in television towers, using SMA dampers under seismic excitations, has not been conducted. To this end, the response reduction in a flexible television tower, disturbed by earthquakes using SMA dampers, is conducted in this study. A two-dimensional dynamic model is developed for dynamic computation at first. The mathematical model of an SMA damper is proposed, and the equations of motion of the tower, without and with, are established, respectively. The structural dynamic responses are examined in the time and the frequency domain, respectively. The effects of damper stiffness, service temperature, hysteresis loops, and earthquake intensity on control efficacy are investigated in detail. In addition, the power spectrum density curves, of dynamic responses and the energy responses, are compared to provide deep insights into the developed control approach. The control performance of SMA dampers is compared with that of widely-used friction dampers. The analytical observations indicate that SMA dampers with optimal parameters can substantially reduce the vibrations of TV transmission towers under seismic excitations.

1. Introduction

To be important high-rise infrastructures, television (TV) transmission towers are widely used across the world for TV transmission and sightseeing. These TV transmission towers, with low damping, may vibrate excessively in severe environments, inducing possible damage and even failure [1,2]. For instance, a real steel TV transmission tower, with a height of 136m, in China collapsed under strong winds in 2012 [3]. The collapse of transmission towers subjected to strong earthquakes in China is also reported [4]. The simplest approach to reduce vibration of TV transmission towers is to increase inherent structural ductility, although it is often uneconomical. Thus, various control devices, such as dynamic absorbers and energy-dissipating dampers, are developed to protect TV transmission towers from excessive vibrations and damages [5,6,7]. The dynamic absorbers are firstly applied in the vibration control of high-rise towers. The wind-induced responses of high-rise towers are studied by Hitchcock et al. [8] and Balendra et al. [9], using tuned liquid dampers and active tuned liquid column dampers. Similar studies are carried out using tuned mass dampers (TMD), such as Yang et al. [10] for a steel tower, Wu and Yang [11] for the Nanjing TV tower, Ghorbani-Tanha et al. [12] for the Milad Tower, and Lu et al. [13] for Shanghai World Financial Center Tower. An obvious defect of dynamic absorbers for vibration control of TV transmission towers is the large space requirement, which often conflicts with the structural function in operation.
As an alternative approach, friction dampers and fluid dampers have been utilized in the structural vibration control of truss towers for many years [14,15,16,17]. Chen et al. [18] conducted the seismic response mitigation of TV transmission towers by using friction dampers. Similar control approaches have been developed by fluid dampers due to their simple configuration and low cost. For instance, Zhang and Li [19] investigated the responses, of a flexible tower under earthquakes, by fluid viscous dampers. Recently, smart materials are widely used in the field of vibration engineering to fabricate various smart dampers for structural vibration control. Xu et al. [20] carried out the semi-active control of a wind-excited tower by piezoelectric friction dampers, associated with a local feedback control algorithm. Chen et al. [21] conducted the semi-active control of a flexible tower subjected to winds using magnetorheological dampers.
A typical smart material, shape memory alloy (SMA) has been recently gaining great popularity due to the properties of high strength, super-elasticity, energy dissipation, better fatigue, and corrosion resistance [22,23,24,25,26]. SMA can exhibit excellent super-elasticity, and this feature is explored to fabricate SMA energy dissipation systems, including SMA bracings, SMA dampers, and connection elements, for vibration control. Zheng et al. [27] proposed a self-adaptive Ni-Ti SMA-cable-based frictional sliding bearing to mitigate pounding and unseating damage of bridges subjected to earthquakes. Liu et al. [28] investigated the seismic response control of a steel frame using base isolation systems with SMA springs.
However, the seismic response suppression, of TV transmission towers with SMA dampers, has not yet been studied. In this regard, the seismic behavior of TV transmission towers with SMA dampers is examined in this study. The two-dimensional (2D) lumped mass model is developed based on the three-dimensional (3D) finite element (FE) model for structural dynamic computation. The model of an SMA damper is proposed, and the equations of motion of the tower, without and with control, are established, respectively. Three schemes are studied and compared to determine the rational location of SMA dampers. The effects of damper stiffness, service temperature, hysteresis loops, and earthquake intensity on control efficacy are examined through a detailed parametric study. In addition, the power spectrum density (PSD) curves of dynamic responses are examined to explore the vibration properties in the frequency domain. The control performance of SMA dampers is compared with that of widely-used friction dampers. The energy responses of the uncontrolled and controlled towers are also compared to provide deep insights into the developed control approach. The made observations indicate that SMA dampers with optimal parameters are beneficial to the vibration control of TV transmission towers under seismic excitations.

2. Seismic Responses of a TV Transmission Tower

The 3D FE model of the TV transmission tower is firstly constructed with the aid of the commercial package ABAQUS, as displayed in Figure 1a. The configuration of the TV transmission tower is quite complicated, and the number of the dynamic degree-of-freedoms (DOFs) of the 3D FE model is quite large. Thus, the dynamic response computation is time-consuming, and the displacement-force relationship of SMA dampers is nonlinear. The dynamic analysis of the tower-damper system needs a lot of iterations at each time step. Thus, the response computation and parametric study on control performance are quite tedious and time-consuming for a complicated 3D FE model [18,20]. To this end, 2D lump models, with satisfactory accuracy, are commonly utilized to simplify dynamic computation and parametric study [18,20], as shown in Figure 1b. Thirty nodal floors are determined, and each floor is assumed to have only one dynamic DOF. The stiffness matrix is determined based on the flexibility matrix of the original 3D FE model.
During the construction of the 2D lumped dynamic model of the TV transmission tower, the structural masses are concentrated at limited nodal floors, and the mass matrix of the reduced 2D dynamic model is established using a lumped mass matrix. The stiffness matrix of the 2D lumped model can be determined based on the flexibility matrix of the 3D FE model of the tower. Apply the same horizontal unit force at the thirty node floors one by one and compute the flexibility coefficient of each floor. Then, form the flexibility matrix of the tower and inverse the flexibility matrix to obtain the stiffness matrix of the 2D model.

3. Model of SMA Damper

The model of SMA material can be described based on thermodynamic theory, and the commonly accepted constitutive model is:
σ σ 0 = E ( ξ ) ε E ( ξ 0 ) ε 0 + Ω ( ξ ) ξ S Ω ( ξ 0 ) ξ S 0 + Θ ( T T 0 )
where σ and ε are the stress and strain of the SMA wire, respectively; σ0 and ε0 are the stress and strain of the SMA wire at the initial state, respectively; ξ is the volume percentage of the martensite; ξ0 is the martensite content of the SMA wire at initial state; T0 is the temperature of the SMA wire at initial state; ξS0 is the initial volume content of the martensite; ξS is the volume content of the martensite; E is the Young’s modulus of SMA material; Ω is the phase transition tensor; Θ is the thermo-elastic tensor of SMA material; T is the temperature at a reference state.
The Young’s modulus of an SMA wire, related to the martensite content E(ξ), is given by
E ( ξ ) = ξ E M + ( 1 ξ ) E A
where EA and EM are Young’s moduli at the austenite and martensite phases, respectively.
The inherited properties of SMA wires, such as self-centering, better fatigue, and corrosion resistance, are useful in the design of energy-dissipating dampers. An SMA damper with re-centering capacity is designed according to the super-elasticity and high damping characteristics of the SMA wire, as shown in Figure 2. The SMA damper consists of the outer tube, inner tube, SMA wires, and two circular plates. The SMA wire is designed always in tension when the damper works. The SMA wire can dissipate vibrant energy, during its reciprocal movement, to suppress the structural vibration.
The hysteretic model of the Ni-Ti SMA material is displayed in Figure 3. In the figure, Mf and Ms denote the martensite finish temperature and martensite starts temperature, respectively. Af and As denote the austenite finish temperature and austenite start temperature, respectively; σ M s and ε M s represent the critical stress and strain of the Ni-Ti SMA material at Ms, respectively; σ M f and ε M f are the critical stress and strain at Mf, respectively; σ A s and ε A s are the critical stress and strain at As, respectively; σ A f and ε A f are the critical stress and strain at Af, respectively; εL is the maximum residual strain.
When the working temperature of the SMA material is given, the mechanical property of the SMA material will be known [29,30,31]. As illustrated in Figure 3, the relationships of strain and stress, on different paths associated with Ni-Ti SMA material, can be expressed as follows:
(1)
Paths O-A and E-O are elastic (full austenite):
σ ( t ) = E A ε ( t )
The control force of the SMA damper u(t) is:
u ( t ) = E A A l s ( x d ( t ) l s )
where x d ( t ) is the length of the damper after deformation; A is the cross-sectional area of SMA wires; ls is the initial length of the SMA wires.
(2)
Path A-B (forward transformation):
σ ( t ) = σ M s + σ M f σ M s ε M f ε M s ( ε ( t ) ε M s )
u ( t ) = σ ( t ) A
(3)
Path B-C is elastic (full martensite):
σ ( t ) = σ M f + E M ( ε ( t ) ε M f )
u ( t ) = E M A l s ( x d ( t ) l s )
(4)
Path B-D is elastic (full martensite)
σ ( t ) = σ M f + E M ( ε ( t ) ε M f )
u ( t ) = σ M f A + E M [ ( x d ( t ) l s ) ε M f ] A
(5)
Path D-E (reverse transformation):
σ ( t ) = σ M s + σ M s σ A f ε A s ε A f ( ε ( t ) ε A s )
u ( t ) = σ ( t ) A
in which σ(t) and ε(t) are the stress and the strain of SMA wires, respectively.

4. Dynamic Analysis of the Controlled Tower

The equations of motion of the controlled tower with SMA dampers is given by:
M x ¨ ( t ) + C x ˙ ( t ) + K x ( t ) = M I g x ¨ g ( t ) + H u ( t )
u ( t ) = [ u 1 u 2 u m ]
where x ( t ) , x ˙ ( t ) , and x ¨ ( t ) are the displacement, velocity, and acceleration responses of the TV transmission tower, respectively; M, and K are the mass and stiffness matrices of the tower, respectively; C is the Rayleigh damping matrix; Ig is the influence vector for ground acceleration; x ¨ g ( t ) is the seismic acceleration; H is the position matrix of control force; u(t) is the control force of SMA dampers; m is the number of SMA dampers.
The incremental equations of the TV transmission tower are given by:
M Δ x ¨ ( t ) + C Δ x ˙ ( t ) + K Δ x ( t ) = M I g Δ x ¨ g ( t ) + H Δ u ( t )
where Δ x ( t ) , Δ x ˙ ( t ) , and Δ x ¨ ( t ) are the displacement increment, velocity increment, and acceleration increment, respectively; Δ x ¨ g ( t ) is the increment of the ground acceleration.
The above incremental equation is solved by the Newmark-β method with a time step of Δt. It is noted that the increment vector of SMA dampers Δu(t) can only be determined under the known Δx(t) and Δ x ˙ ( t ) in the numerical computation. However, the displacement increment Δx(t) and velocity increment Δ x ˙ ( t ) is commonly determined based on Δu(t). Thus, a dynamic iteration process is developed to ensure computational convergence. The iteration process proceeds until Δx(t) satisfies the following equation:
a b s [ ( Δ x ( t ) ( i + 1 ) Δ x ( t ) ( i ) ) / Δ x ( t ) ( i ) ] < κ
where κ is a prescribed error limitation related to an actual TV transmission tower.
Energy response is an important index to describe the vibration of TV transmission towers. The energy representation of the example tower is formed by integrating Equation (15). The absolute energy equation of the uncontrolled tower is given by:
E K E + E D E + E S E = E E
Similarly, the energy equation for the tower with SMA dampers is:
E K E + E D E + E S E + E P E = E E
where EE is the total input energy from earthquakes; EKE and ESE are kinetic energy and strain energy; EDE and EPE are the energy dissipated by structural damping and SMA dampers, respectively.

5. Case Study

5.1. Analytical Parameters

A real TV tower, with a height of 340 m, is used to examine the feasibility of the proposed control method, as displayed in Figure 1. Two turrets are developed for sightseeing. The 3D model of the TV tower is established by ABAQUS, and then, a 2D dynamic model is constructed by a developed computer program. The validity of the 2D model is examined by comparing structural dynamic characteristics, as listed in Table 1. The fundamental frequency of the 3D and 2D models are 0.129 Hz and 0.128 Hz, respectively, with an error of 0.621%. It is also found that the average error of the first eight natural frequencies is 3.641%. Thus, the dynamic characteristics of the 2D dynamic model are close to that of the original 3D model, and the validity of the 2D model can be guaranteed. To evaluate the effectiveness of the proposed control method, three far-field and two near-field historical records are selected [32,33,34,35]: (1) Taft earthquake on 21 July 1952; (2) El Centro earthquake on 18 May 1940; (3) Kobe earthquake on 17 January 1995; (4) Hachinohe earthquake on 16 May 1968; (5) Northridge earthquake on 17 January 1994, respectively. The original peak ground accelerations (PGAs) of the five seismic records are 1.559, 3.417, 8.178, 2.250, and 8.2676 m/s2, as shown in Figure 4. All the earthquakes are scaled to have the same PGA of 4.0 m/s2 for comparison. The time step and damping ratios in the dynamic analysis are 0.02 s and 0.01, respectively.
Obvious whiplash effects of the mast are observed based on analytical results. Owing to the small geometric size of the mast, sixty SMA dampers are evenly divided into four groups and installed around the two turrets. Thirty SMA dampers (group nos. 1 and 2) are installed at the low turret, and the other two groups of SMA dampers are at the upper turret. An SMA damper with an axial brace connects spherical nodes in the axial direction of a member. The physical parameters of SMA dampers are listed in Table 2. The Young’s modulus of the SMA damper brace is 2.12 × 108 kN/m, and the cross-sectional area of the SMA damper is 50 cm2.

5.2. Control Performance Comparison

The vibration reduction factor (VRF) is used in this study to evaluate the control efficacy of SMA dampers:
VRF = R n c R c o R n c
where Rnc and Rco are the peak response without and with SMA dampers, respectively.
The curves in Figure 5 indicate that the peak responses are remarkably mitigated using SMA dampers. It is found that the displacement VRFs of the tower top, upper turret, and low turret are 42.739%, 9.420%, and 37.856%, respectively. Similarly, the counterparts of the acceleration are 28.937%, 9.533%, and 26.405%, respectively. In addition, the displacement and acceleration VRFs are slightly better than those of velocity and responses. As a result, satisfactory control can be achieved using SMA dampers, particularly for the whiplash effects of the mast. The dynamic responses of the tower top and the upper turret, with and without SMA dampers under the Taft earthquake, are displayed in Figure 6. Owing to the additional damping produced by SMA dampers, the vibrations of the controlled tower damped out more quickly than that of the original tower. The application of SMA dampers can mitigate the seismic responses of the TV transmission tower to a great extent. In particular, the control efficacy on whiplash effects of the mast is better than that of the other parts of the tower.
Three damper installation schemes are considered to compare the effects of damper position on control efficacy, as listed in Table 3. Scheme 1 is mentioned above, in which sixty dampers are around the two turrets. For scheme 2, all the SMA dampers are installed around the low turret (Mass No. 4). For scheme 3, all the SMA are installed around the upper turret.
Figure 7 shows the structural peak responses of different control schemes. For scheme 2, with all the dampers around the low turret, the control efficacy of displacement responses is much better than that of velocity and acceleration responses. Similar observations can be made for scheme 3, whose control performance is slightly inferior to that of scheme 2. The overall control efficacy of scheme 1 is better than that of schemes 2 and 3. Thus, the optimal positions for SMA dampers are both two turrets instead of a single turret. It is also observed that the control efficacy of displacement of all the schemes is slightly better than that of velocity and acceleration responses. Therefore, scheme 1 is adopted for the parametric study in the following sections.
The control efficacy of SMA dampers under the Taft earthquake is also analyzed and demonstrated in the frequency domain. The power spectral density (PSD) curves of dynamic responses of the tower top and upper turret are displayed in Figure 8. The magnitudes of the PSD curves, for all the three types of responses with control, are much smaller compared with those without control. In particular, the PSD VRFs of the first vibration mode are larger than those of other higher modes.

6. Parametric Study on Control Performance

6.1. Stiffness Effects of SMA Damper

The stiffness coefficient of an SMA damper is given by:
SC = k d S M A k 0 S M A
where k d S M A is the stiffness value in computation; k 0 S M A is the initial stiffness of an SMA damper.
The stiffness effects of SMA dampers on the peak responses under the Taft earthquake are investigated, and the results are displayed in Figure 9. The peak displacement and velocity responses of the mast are gradually reduced with an increasing SC until approximately 1.0. A further increase in stiffness of SMA dampers does not generate further significant response reduction when SC is larger than 1.0. The peak acceleration of the mast reduces at first and then increases with the increasing SC, which means there exist optimal SC values for structural peak responses. It is not beneficial to set a very large SC value to avoid unnecessary cost waste. Similar results can be achieved from the dynamic responses of the tower body and turrets. Thus, an optimal SC value for all the responses does not exist, and optimal SC values for three kinds of peak responses, at different places, are different to some extent. The optimum SC value for the example tower is selected as 1.0 considering the overall control efficacy of the tower.

6.2. Properties of Hysteresis Loops

The hysteresis loop is an important index to describe the energy-dissipating capacity of an SMA damper under earthquakes. The variation of a hysteresis loop with damper stiffness is displayed in Figure 10. A small SC value (SC = 0.6) reflects the small enclosed area of the hysteresis loop. Thus, the SMA damper is easy to extend while the energy-dissipating capacity is poor. Increasing the damper stiffness (SC = 1.0 or 2.0) can improve the shape of the hysteresis loops and the energy-absorbing capacity. After the SC value reaches 2.0, a further increase (SC = 3.0) in damper stiffness may increase the peak control force but reduce the enclosed area of the hysteresis loop. Therefore, an optimal SC value of an SMA damper can also be determined based on the configuration of hysteresis loops.

6.3. Effect of Service Temperature

The variation of damper force and deformation with service temperature is displayed in Figure 11. The peak forces of damper No. 40 are 163.71, 193.82, 224.87, and 284.08 kN, respectively, at 0, 10, 20, and 40 degrees Celsius. With the increase in service temperature, the peak damper force gradually increases, while the damper deformation decreases to some extent. Similar observations can be made for the other SMA dampers. The effects of service temperature on damper performance can also be found from the variation of damper hysteresis loops. The increasing service temperature may induce the increment in damper force and reduction in damper deformation. The enclosed areas of damper hysteresis loops reduce to some extent. Thus, a very large service temperature changes the damper states, and the damper may behave as a steel brace does and lose its energy-dissipating capacity.
The effects of service temperature on maximum responses of the TV transmission tower are studied, as shown in Figure 12. The responses of the tower top, mast, tower body, and two turrets are computed, respectively, under the common temperature range between 0 to 40 degrees Celsius. The displacement responses of the TV tower are almost stable with increasing service temperature, while the velocity and acceleration responses slightly increase. It is noted that a relatively large service temperature may affect the shape of damper hysteresis loops, as displayed in Figure 11, while the control performance is stable. This is because the peak forces of SMA dampers increase automatically with increasing temperature.

6.4. Effects of Earthquake Intensity

The effects of PGA on control performance are investigated in this section. Figure 13 displays the variation of peak control force of SMA dampers with PGA. The deformation of SMA dampers almost increases linearly with the increasing PGA under constant damper stiffness. The peak control force of the SMA damper gradually increases with the increasing earthquake intensity. This means that the performance of SMA dampers presents certain robustness under external excitations.
As displayed in Figure 14, the displacement VRF of the tower top gradually increases, with the increment in PGA, while the velocity and acceleration VRFs almost keep constant. All three types of VRFs of the mast are stable with increasing earthquake intensity. Similar observations can be made for the upper turret and low turret. Thus, it is noted that the overall VRFs of the TV transmission tower are stable when subjected to different seismic excitations. The control performance by SMA dampers is robust and versatile for the earthquake-disturbed tower.

6.5. Efficacy Subjected Other Earthquakes

The seismic responses under the other three earthquakes are also computed and compared with the same PGA of 4.0 m/s2. The response VRFs of the example tower under different earthquakes are listed in Table 4. Three far-field and two near-field historical records are selected to evaluate the performance of the proposed control method based on SMA dampers [32,33,34,35]. The SMA dampers can reduce structural seismic responses while the control efficacy, under different seismic inputs, is different. In addition, the overall VRFs of the tower top and mast are much better than those of the two turrets. Thus, the whiplash effects of the example tower can be substantially reduced. In reality, the different frequency components of the example earthquakes in the frequency domain induce the difference in structural dynamic responses. Therefore, control performances under different earthquakes cannot maintain optimality at all times.

6.6. Comparison between SMA Damper and Friction Damper

The dynamic responses of the tower with the commonly-used friction damper are computed and compared with those by SMA dampers. Friction dampers have been used to improve the seismic resistant performance of various structures during the past three decades. The Coulomb law is commonly adopted in the mathematical model of friction dampers. It presumes that the frictional force is independent of the velocity and the coefficient μ is the proportional constant considered in the mathematical relation between the frictional force and the clamping force N. The frictional force fd equals the clamping force N multiplies the coefficient of friction μ. The clamping force of a friction damper is set in advance, and its frictional force is constant during the slip process. If the force of the friction damper fd is no larger than the design frictional force, then the damper is sticking and the damper behaves as if it is a common brace. If the damper force fd is larger than the design frictional force, the damper begins to slip and dissipate vibrant energy. During the structural vibration under an earthquake, if the relative displacement increment of the two ends of the friction element is large than zero, the control force is tensile. Otherwise, the control force is compressive. As displayed in Figure 10b and Figure 11a, the peak force of an SMA damper is 150 kN. Thus, the peak control forces of friction dampers are set as the same for performance comparison. The number and location of friction dampers are also the same as those of SMA dampers.
The hysteresis loops of friction dampers, under different damper forces, are displayed in Figure 15. For a friction damper, its damper force cannot be changed during the structural vibration. If the damper slipping force is too large (250 kN), the friction damper is hard to slip and behaves similarly to a pure brace. If the damper slipping force is too small (50 kN), the damper is very easy to slip while the enclosed area is small, which means a limited energy-dissipating capacity of the friction damper. Thus, the optimal slipping force of a friction damper depends on the intensity of the external excitations.
The hysteresis loops of friction dampers, under different earthquake intensities, are displayed in Figure 16. For a small seismic intensity (PGA = 1.0 or 2.0 m/s2), the damper force of 150 kN is large enough, the damper is hard to slip, and it behaves as a brace. For an appropriate seismic intensity (PGA = 4.0 m/s2), the hysteresis loop is plump, and the control performance is satisfactory. For a larger seismic intensity (PGA = 5.0 m/s2), the damper can slip more easily, the enclosed area of the hysteresis loop is limited, and the energy-dissipating capacity and control performance is limited. The reason is that the slipping force of a friction damper is set in advance and cannot change with the varying excitation intensity. Therefore, the control robustness of friction dampers is unsatisfactory.
The performance of SMA dampers is quite different in comparison with that of friction dampers. The variation of control force with seismic intensity is investigated and shown in Figure 17. It is seen that the control force of the SMA dampers can be modulated automatically with the earthquake-induced structural and damper responses. The control forces of SMA dampers quickly increase with the increasing seismic intensity. Thus, the SMA damper can move more easily under any seismic intensity. In addition, the shapes of hysteresis loops of SMA are similar, even if the seismic intensity increases. However, the performance of friction dampers is quite different in comparison with that of SMA dampers, as displayed in Figure 16 and Figure 17. The control forces of the friction dampers are set at 150 kN in advance and cannot be changed during the whole vibration. For a small seismic intensity (PGA = 1.0 or 2.0 m/s2), the friction dampers cannot slip, and they behave as a brace. Thus, the damper force is smaller than 150 kN. The energy-dissipating capacity of SMA dampers is stable and superior to that of friction dampers.
The structural VRFs, under different seismic intensities, are computed and compared, as displayed in Table 5. It is observed that the VRFs of SMA dampers are almost stable for different seismic intensities. The SMA damper-based control approach is robust and versatile for the TV transmission tower. However, the control efficacy of friction dampers with constant slipping force (150 kN) may dramatically reduce with the varying seismic intensity. If the damper force is a certain value and cannot be automatically changed, a friction damper is easy in the sticking stage and behaves much as a common brace does under a small seismic intensity. Thus, the friction damper is unable to fully slip, which causes unsatisfactory control performance. However, an SMA damper with proper parameters can automatically adjust its control force in line with the change of external loading and thereby ensure robust control performance.

7. Energy Responses of the Tower

7.1. Comparison of Energy Responses of the Tower

The vibrant intensity of the TV transmission tower can also be depicted from the viewpoint of energy. The energy properties of the example tower, subjected to various earthquakes, are displayed in Figure 18. For the uncontrolled tower under the Taft earthquake, the total inputted energy from the seismic excitations is absorbed by structural damping. Large kinetic and strain energy is seen because of the excessive vibration of the tower. Similar observations can be made from the other three ground motions. The magnitude of total inputted energy to the tower, under various earthquakes, is quite different even though their PGAs are the same. Thus, the frequency components of seismic excitations have substantial influences on structural dynamic responses.
Figure 19 displays the energy responses of the tower with SMA dampers in comparison with those without control. The application of SMA dampers can substantially reduce the vibrant energy of the TV transmission tower to a great extent. To compare the energy curves in Figure 18 and Figure 19, one can find that the inputted energy of the controlled tower is smaller in comparison with the counterpart of the controlled tower. Owing to the aid of SMA dampers, the energy absorbed by structural damping is much less than that of the uncontrolled tower. The energy responses under Northridge earthquake have the similar properties.

7.2. Variation of Energy Response with Damper Stiffness

Curves in Figure 20 indicate that the total inputted energy from earthquakes EE varies with the damper stiffness. For a small damper stiffness (SC = 0.6), the energy-dissipating capacity of the SMA dampers is poor, and the tower energy EE is even slightly larger than that of the uncontrolled tower. The tower energy EE gradually reduces with the increasing SC until approximately 1.0. A further increase in damper stiffness (SC = 2.0 or 3.0) may induce a large value of EE. However, the energy dissipated by SMA dampers EP is quite different. As displayed in Figure 20, EP gradually increases with the increasing damper stiffness until SC reaches about 1.0. If stiffness coefficient SC reaches above 2.0, EP reduces to some extent. The energy dissipated by SMA dampers EP is limited if the damper stiffness is relatively small (SC = 0.6) or large (SC = 2.0 or 3.0). Besides, similar observations can be made from the energy curves of structural damping ED. With the increase in damper stiffness, the energy ED gradually reduces because more structural energy EE is dissipated by SMA dampers. Thus, optimal damper stiffness can also be determined based on structural energy curves.

8. Concluding Remarks

The validity of the proposed control approach for TV transmission towers by SMA dampers is examined in detail. The PSD properties and energy responses of the tower are compared with those of the uncontrolled tower. The results indicate that SMA dampers, with optimal parameters, can be used in the control of TV transmission towers disturbed by earthquakes.
Dynamic responses of the TV tower are remarkably suppressed by SMA dampers particularly for the whiplash effects of the top mast. The overall control efficacy of scheme 1 is better than that of the other two schemes, which means that the optimal positions for SMA dampers are both two turrets. The hysteresis loops of SMA dampers are quite different because the damper stiffness can directly affect the damper’s working status. The optimal damper stiffness can also be determined based on the configuration of hysteresis loops. The vibrant energy of the controlled tower is substantially reduced, and the inputted energy of the controlled tower is also mitigated to a great extent. The application of SMA dampers can remarkably improve the energy-dissipating mechanism of the TV transmission tower. Energy evaluation is an effective approach to determine the optimal damper parameters in vibration control.
It is not beneficial to set a very large SC value for SMA dampers to avoid unnecessary cost waste. Thus, an optimal SC value for all the responses does not exist, and optimal SC values for three kinds of peak responses at different places are different to some extent. Thus, in the real application of SMA dampers, an optimal damper stiffness should be determined based on the structural model of TV transmission tower through detailed parametric study, which may be time-consuming and complicated. Similarly, an optimal service temperature is required based on the structural model of the TV transmission tower through a detailed and complicated parametric study.

Author Contributions

Conceptualization, B.C.; methodology, B.C. and J.W.; software, J.W.; validation, B.C. and L.Z.; formal analysis, J.W. and X.S.; investigation, B.C. and J.W.; resources, J.W. and X.S.; data curation, J.W.; writing—original draft preparation, B.C. and J.W.; writing—review and editing, B.C. and L.Z.; visualization, B.C. and L.Z.; supervision, B.C. and L.Z.; project administration, B.C.; funding acquisition, B.C. All authors have read and agreed to the published version of the manuscript.

Funding

The authors are grateful for the financial support from the National Natural Science Foundation of China (51678463, 51978549).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data sharing is not applicable to this article.

Acknowledgments

The authors sincerely appreciate the selfless support provided by Wuhan University of Technology during the progress of this study.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Amiri, M.M.; Yahyai, M. Estimation of damping ratio of TV towers based on ambient vibration monitoring. Struct. Des. Tall Spec. Build. 2013, 22, 862–875. [Google Scholar] [CrossRef]
  2. Chen, B.; Weng, S.; Zhi, L.H.; Li, D.M. Response control of a large transmission-tower line system under seismic excitations by using friction dampers. Adv. Struct. Eng. 2017, 20, 1155–1173. [Google Scholar] [CrossRef]
  3. Li, T.; Fan, J.; Lv, Y.X. Collapse analysis of a steel TV tower under wind load. Steel Constr. 2017, 32, 80–83. [Google Scholar]
  4. Liu, R.S.; Liu, J.L.; Yan, D.Q. Seismic damage investigation and analysis of electric power system in Lushan MS 7.0 earthquake. J. Nat. Disasters 2013, 22, 85–92. [Google Scholar]
  5. Housner, G.W.; Bergman, L.A.; Caughey, T.K.; Chassiakos, A.G.; Claus, R.O.; Masri, S.F.; Skelton, R.E.; Soong, T.T.; Spencer, B.F., Jr.; Yao, T.P. Structural control: Past, present, and future. J. Eng. Mech. ASCE 1997, 123, 897–971. [Google Scholar] [CrossRef]
  6. Xu, Y.L.; Qu, W.L.; Chen, B. Active/robust moment controllers for seismic response control of a large span building on top of ship lift towers. J. Sound Vib. 2003, 261, 277–296. [Google Scholar] [CrossRef]
  7. Chen, B.; Xia, Y. Advanced technologies in disaster prevention and mitigation. Adv. Struct. Eng. 2017, 20, 1141–1142. [Google Scholar] [CrossRef] [Green Version]
  8. Hitchcock, P.A.; Glanville, M.J.; Kwok, K.C.S.; Watkins, R.D. Damping properties and wind-induced response of a steel frame tower fitted with liquid column vibration absorbers. J. Wind Eng. Ind. Aerodyn. 1999, 83, 183–196. [Google Scholar] [CrossRef]
  9. Balendra, T.; Wang, C.M.; Yan, N. Control of wind-excited towers by active tuned liquid column damper. Eng. Struct. 2001, 23, 1054–1067. [Google Scholar] [CrossRef]
  10. Yang, N.; Wang, C.M.; Balendra, T. Composite mass dampers for vibration control of wind-excited towers. J. Sound Vib. 1998, 213, 301–316. [Google Scholar] [CrossRef]
  11. Wu, J.C.; Yang, J.N. LQG control of lateral-torsional motion of Nanjing TV transmission tower. Earthq. Eng. Struct. Dyn. 2000, 29, 1111–1130. [Google Scholar] [CrossRef]
  12. Ghorbani-Tanha, A.K.; Noorzad, A.; Rahimian, M. Mitigation of wind-induced motion of Milad Tower by tuned mass damper. Struct. Des. Tall Spec. Build. 2009, 18, 371–385. [Google Scholar] [CrossRef]
  13. Lu, X.L.; Li, P.Z.; Guo, X.Q.; Shi, W.X.; Liu, J. Vibration control using ATMD and site measurements on the Shanghai World Financial Center Tower. Struct. Des. Tall Spec. Build. 2014, 23, 105–123. [Google Scholar] [CrossRef]
  14. Pall, A.S.; Marsh, C. Response of friction damped braced frames. J. Struct. Div. ASCE 1982, 108, 1313–1323. [Google Scholar] [CrossRef]
  15. Grigorian, C.E.; Yang, T.S.; Popov, E.P. Slotted bolted connection energy dissipators. Earthq. Spectra 1993, 9, 491–504. [Google Scholar] [CrossRef]
  16. Chen, B.; Xiao, X.; Li, P.Y.; Zhong, W.L. Performance evaluation on transmission-tower line system with passive friction dampers subjected to wind excitations. Shock. Vib. 2015, 2015, 1–13. [Google Scholar] [CrossRef] [Green Version]
  17. Quintana, H.C.; Petkovski, M. Optimum performance of structural control with friction dampers. Eng. Struct. 2018, 172, 154–162. [Google Scholar] [CrossRef]
  18. Chen, B.; Yang, D.; Feng, K.; Ouyang, Y.Q. Response control of a high-rise television tower under seismic excitations by friction dampers. Int. J. Struct. Stab. Dyn. 2018, 18, 1850140. [Google Scholar] [CrossRef]
  19. Zhang, Z.Q.; Li, A.Q. Seismic response vibration control analyses of Hefei TV Tower with fluid viscous damper. In Proceedings of the International Conference on Health Monitoring of Structure, Materials and Environment, Nanjing, China, 16 October 2007. [Google Scholar]
  20. Xu, Y.L.; Qu, W.L.; Chen, Z.H. Control of wind-excited truss tower using semiactive friction damper. J. Struct. Eng. ASCE 2001, 127, 861–868. [Google Scholar] [CrossRef]
  21. Chen, B.; Zheng, J.; Qu, W.L. Control of wind-induced response of transmission tower-line system by using magnetorheological dampers. Int. J. Struct. Stab. Dyn. 2009, 9, 661–685. [Google Scholar] [CrossRef]
  22. Song, G.B.; Ma, N.; Li, H.N. Applications of shape memory alloys in civil structures. Eng. Struct. 2006, 28, 1266–1274. [Google Scholar] [CrossRef]
  23. Song, G.B.; Ma, N. Robust control of a shape memory alloy wire actuated flap. Smart Mater. Struct. 2007, 16, 51–57. [Google Scholar] [CrossRef]
  24. Ozbulut, O.; Hurlebaus, S.; Desroches, R. Seismic response control using shape memory alloys: A review. J. Intell. Mater. Syst. Struct. 2011, 22, 1531–1549. [Google Scholar] [CrossRef]
  25. Fang, C.; Wang, W.; He, C.; Chen, Y. Self-centring behaviour of steel and steel-concrete composite connections equipped with NiTi SMA bolts. Eng. Struct. 2017, 150, 390–408. [Google Scholar] [CrossRef]
  26. Fang, C.; Yam, M.C.; Lam, A.C.; Xie, L. Cyclic performance of extended end-plate connections equipped with shape memory alloy bolts. J. Constr. Steel Res. 2014, 94, 122–136. [Google Scholar] [CrossRef]
  27. Zheng, Y.; Dong, Y.; Chen, B.; Ghazanfar, A.A. Seismic damage mitigation of bridges with self-adaptive SMA-cable-based bearings. Smart Struct. Syst. 2019, 24, 127–139. [Google Scholar]
  28. Liu, Y.P.; Wang, H.Y.; Qiu, C.X.; Zhao, X.N. Seismic behavior of superelastic shape memory alloy spring in base isolation system of multi-story steel frame. Materials 2019, 12, 997. [Google Scholar] [CrossRef] [Green Version]
  29. Brinson, L.C. One-dimensional constitutive behavior of shape memory alloys: Thermomechanical derivation with non-constant material functions and redefined martensite internal variable. J. Intell. Mater. Syst. Struct. 1993, 4, 229–242. [Google Scholar] [CrossRef]
  30. Zheng, Y.; Dong, Y.; Li, Y. Resilience and life-cycle performance of smart bridges with Shape Memory Alloy (SMA)-cable-based bearings. Constr. Build. Mater. 2018, 158, 389–400. [Google Scholar] [CrossRef]
  31. Qiu, C.; Fang, C.; Liang, D.; Du, X.; Yam, M.C.H. Behavior and application of self-centering dampers equipped with buckling-restrained SMA bars. Smart Mater. Struct. 2020, 29, 035009. [Google Scholar] [CrossRef]
  32. Spencer, B.F., Jr.; Dyke, S.J.; Deoskar, H.S. Benchmark problems in structural control. I: Active mass driver system. Earthq. Eng. Struct. Dyn. 1998, 27, 1127–1139. [Google Scholar] [CrossRef]
  33. Spencer, B.F., Jr.; Dyke, S.J.; Deoskar, H.S. Benchmark problems in structural control. II. Active tendon system. Earthq. Eng. Struct. Dyn. 1998, 27, 1141–1147. [Google Scholar] [CrossRef]
  34. Ohtori, Y.; Christenson, R.E.; Spencer, B.F., Jr.; Dyke, S.J. Benchmark control problems for seismically excited nonlinear buildings. J. Eng. Mech. ASCE 2004, 130, 366–385. [Google Scholar] [CrossRef] [Green Version]
  35. Chen, C.C.; Chen, G.D. Nonlinear control of a 20-story steel building with active piezoelectric friction dampers. Struct. Eng. Mech. 2002, 14, 21–38. [Google Scholar] [CrossRef]
Figure 1. Analytical model of a TV transmission tower. (a) 3D FE model; (b) 2D dynamic model.
Figure 1. Analytical model of a TV transmission tower. (a) 3D FE model; (b) 2D dynamic model.
Materials 14 06987 g001
Figure 2. Schematic configuration of SMA damper.
Figure 2. Schematic configuration of SMA damper.
Materials 14 06987 g002
Figure 3. Flag-shaped hysteretic model of an SMA wire.
Figure 3. Flag-shaped hysteretic model of an SMA wire.
Materials 14 06987 g003
Figure 4. Time histories of four historical seismic records. (a) Taft; (b) El Centro; (c) Kobe; (d) Hachinohe; (e) Northridge.
Figure 4. Time histories of four historical seismic records. (a) Taft; (b) El Centro; (c) Kobe; (d) Hachinohe; (e) Northridge.
Materials 14 06987 g004
Figure 5. Maximum responses of the tower without/with control. (a) Displacement; (b) Velocity; (c) Acceleration.
Figure 5. Maximum responses of the tower without/with control. (a) Displacement; (b) Velocity; (c) Acceleration.
Materials 14 06987 g005
Figure 6. Time histories of dynamic responses of the television tower. (a) Displacement No. 30; (b) Displacement No. 22; (c) Velocity No. 30; (d) Velocity No. 22; (e) Acceleration No. 30; (f) Acceleration No. 22.
Figure 6. Time histories of dynamic responses of the television tower. (a) Displacement No. 30; (b) Displacement No. 22; (c) Velocity No. 30; (d) Velocity No. 22; (e) Acceleration No. 30; (f) Acceleration No. 22.
Materials 14 06987 g006
Figure 7. Comparison of different control schemes. (a) Displacement; (b) Velocity; (c) Acceleration.
Figure 7. Comparison of different control schemes. (a) Displacement; (b) Velocity; (c) Acceleration.
Materials 14 06987 g007
Figure 8. Comparison of PSD curves of response time histories. (a) Displacement PSD No. 30; (b) Displacement PSD No. 22; (c) Velocity PSD No. 30; (d) Velocity PSD No. 22; (e) Acceleration PSD No. 30; (f) Acceleration PSD No. 22.
Figure 8. Comparison of PSD curves of response time histories. (a) Displacement PSD No. 30; (b) Displacement PSD No. 22; (c) Velocity PSD No. 30; (d) Velocity PSD No. 22; (e) Acceleration PSD No. 30; (f) Acceleration PSD No. 22.
Materials 14 06987 g008aMaterials 14 06987 g008b
Figure 9. Effects of damper stiffness on peak responses. (a) Displacement at mast; (b) Displacement at tower; (c) Velocity at mast; (d) Velocity at tower; (e) Acceleration at mast; (f) Acceleration at tower.
Figure 9. Effects of damper stiffness on peak responses. (a) Displacement at mast; (b) Displacement at tower; (c) Velocity at mast; (d) Velocity at tower; (e) Acceleration at mast; (f) Acceleration at tower.
Materials 14 06987 g009
Figure 10. Variation of hysteresis loop with damper stiffness. (a) SC = 0.6; (b) SC = 1.0; (c) SC = 2.0; (d) SC = 3.0.
Figure 10. Variation of hysteresis loop with damper stiffness. (a) SC = 0.6; (b) SC = 1.0; (c) SC = 2.0; (d) SC = 3.0.
Materials 14 06987 g010
Figure 11. Variation of damper force and deformation with service temperature. (a) Peak force; (b) Peak deformation; (c) T = 0 °C; (d) T = 10 °C; (e) T = 20 °C; (f) T = 40 °C.
Figure 11. Variation of damper force and deformation with service temperature. (a) Peak force; (b) Peak deformation; (c) T = 0 °C; (d) T = 10 °C; (e) T = 20 °C; (f) T = 40 °C.
Materials 14 06987 g011
Figure 12. Effects of service temperature on maximum responses of the tower. (a) Displacement; (b) Velocity; (c) Acceleration.
Figure 12. Effects of service temperature on maximum responses of the tower. (a) Displacement; (b) Velocity; (c) Acceleration.
Materials 14 06987 g012
Figure 13. Variation of peak control force of SMA dampers with PGA. (a) Peak force; (b) Peak deformation.
Figure 13. Variation of peak control force of SMA dampers with PGA. (a) Peak force; (b) Peak deformation.
Materials 14 06987 g013
Figure 14. Variation of VRFs with seismic PGA. (a) VRF No. 30; (b) VRF No. 28; (c) VRF No. 22; (d) VRF No. 4.
Figure 14. Variation of VRFs with seismic PGA. (a) VRF No. 30; (b) VRF No. 28; (c) VRF No. 22; (d) VRF No. 4.
Materials 14 06987 g014
Figure 15. Hysteresis loops of friction dampers under different damper forces. (a) fd = 50 kN; (b) fd = 150 kN; (c) fd = 250 kN.
Figure 15. Hysteresis loops of friction dampers under different damper forces. (a) fd = 50 kN; (b) fd = 150 kN; (c) fd = 250 kN.
Materials 14 06987 g015
Figure 16. Hysteresis loops of friction dampers under different earthquake intensities. (a) PGA = 1.0 m/s2; (b) PGA = 2.0 m/s2; (c) PGA = 3.0 m/s2; (d) PGA = 4.0 m/s2.
Figure 16. Hysteresis loops of friction dampers under different earthquake intensities. (a) PGA = 1.0 m/s2; (b) PGA = 2.0 m/s2; (c) PGA = 3.0 m/s2; (d) PGA = 4.0 m/s2.
Materials 14 06987 g016
Figure 17. Variation of control force with seismic intensity.
Figure 17. Variation of control force with seismic intensity.
Materials 14 06987 g017
Figure 18. Energy responses of the tower without friction dampers. (a) Taft; (b) El Centro; (c) Kobe; (d) Hachinohe.
Figure 18. Energy responses of the tower without friction dampers. (a) Taft; (b) El Centro; (c) Kobe; (d) Hachinohe.
Materials 14 06987 g018
Figure 19. Energy responses of the tower with SMA dampers. (a) Taft; (b) El Centro; (c) Kobe; (d) Hachinohe.
Figure 19. Energy responses of the tower with SMA dampers. (a) Taft; (b) El Centro; (c) Kobe; (d) Hachinohe.
Materials 14 06987 g019aMaterials 14 06987 g019b
Figure 20. Variation of energy responses with damper stiffness. (a) Tower energy; (b) Energy by structural damping; (c) Energy by damper.
Figure 20. Variation of energy responses with damper stiffness. (a) Tower energy; (b) Energy by structural damping; (c) Energy by damper.
Materials 14 06987 g020
Table 1. Comparison of dynamic properties of the 3D and 2D models.
Table 1. Comparison of dynamic properties of the 3D and 2D models.
Freq.no3D Finite Element Model2D Dynamic ModelError
10.1290.1280.621%
20.4160.3867.211%
30.7290.781−7.133%
41.1711.1630.683%
51.5141.4563.831%
62.2212.1423.556%
72.5782.5690.349%
82.6322.783−5.737%
Table 2. Physical parameters of SMA dampers.
Table 2. Physical parameters of SMA dampers.
ParameterValueParameterValueParameterValue
Mf−46 °CCM10 MPa/°C σ s c r 120 MPa
As−18.5 °CCA15.8 MPa/°C σ f c r 190 MPa
Ms−37.4 °CDA75,000 MPaεL0.079
Af−6 °CDM29,300 MPa
Table 3. Damper installation schemes.
Table 3. Damper installation schemes.
Scheme 1Scheme 2Scheme 3
30 dampers around the low turret (Mass No. 4)
30 dampers around the upper turret (Mass No. 24)
60 dampers around the low turret (Mass No. 4)60 dampers around the upper turret (Mass No. 24)
Table 4. VRFs of the TV transmission tower under different earthquakes (%).
Table 4. VRFs of the TV transmission tower under different earthquakes (%).
Location TaftEl CentroKobeHachinoheNorthridge
Mass No. 30
(Tower top)
Displacement 39.6431.9522.3823.2429.03
Velocity25.5713.8718.7033.895.59
Acceleration24.756.4321.0913.8911.36
Mass No. 28
(Mast)
Displacement 30.7019.4623.6422.5610.70
Velocity26.59−5.2120.6919.6011.42
Acceleration22.48−9.2817.5218.5115.04
Mass No. 22
(Upper turret)
Displacement 9.729.248.318.9110.77
Velocity−10.1911.8823.94−1.05−0.78
Acceleration6.958.8610.48−4.565.55
Mass No. 4
(Low turret)
Displacement 31.88−2.8135.342.694.04
Velocity32.14−10.0523.301.080.42
Acceleration25.30−5.5615.608.74−1.64
Table 5. Comparison of VRFs with different control approaches under different earthquake intensities.
Table 5. Comparison of VRFs with different control approaches under different earthquake intensities.
Location SMA Damper
PGA = 4.0
Friction Damper
PGA = 4.0
SMA Damper
PGA = 1.0
Friction Damper
PGA = 1.0
Mass No. 30
(Tower top)
Displacement39.6439.8753.9240.69
Velocity25.5738.8421.3733.72
Acceleration24.7552.6327.4937.30
Mass No. 28
(Mast)
Displacement 30.7047.0730.998.61
Velocity26.5949.5131.1231.33
Acceleration22.4844.7520.3323.28
Mass No. 22
(Upper turret)
Displacement 9.728.619.688.57
Velocity−10.19−3.93−13.67−5.77
Acceleration6.959.113.5299.08
Mass No. 4
(Low turret)
Displacement 31.8839.4335.47133.37
Velocity32.1430.8031.3736.74
Acceleration25.3029.6129.0223.16
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wu, J.; Chen, B.; Zhi, L.; Song, X. Seismic Response Mitigation of a Television Transmission Tower by Shape Memory Alloy Dampers. Materials 2021, 14, 6987. https://0-doi-org.brum.beds.ac.uk/10.3390/ma14226987

AMA Style

Wu J, Chen B, Zhi L, Song X. Seismic Response Mitigation of a Television Transmission Tower by Shape Memory Alloy Dampers. Materials. 2021; 14(22):6987. https://0-doi-org.brum.beds.ac.uk/10.3390/ma14226987

Chicago/Turabian Style

Wu, Jingbo, Bo Chen, Lunhai Zhi, and Xinxin Song. 2021. "Seismic Response Mitigation of a Television Transmission Tower by Shape Memory Alloy Dampers" Materials 14, no. 22: 6987. https://0-doi-org.brum.beds.ac.uk/10.3390/ma14226987

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

Article Metrics

Back to TopTop