Dynamic Properties and Energy Dissipation Study of Sandwich Viscoelastic Damper Considering Temperature Inﬂuence

: Viscoelastic dampers are a kind of classical passive energy dissipation and vibration control devices which are widely utilized in engineering ﬁelds. The mechanical properties and energy dissipation capacity of the viscoelastic damper are signiﬁcantly affected by ambient temperature. In this work, dynamic properties tests of the sandwich type viscoelastic damper at different environmental temperatures are carried out. The equivalent fractional Kelvin model which can characterize the mechanical behavior of the viscoelastic damper with varying frequencies and temperatures is introduced to describe the dynamic properties and energy dissipation capability of the sandwich viscoelastic damper. The self-heating phenomenon of the sandwich viscoelastic damper is studied with a numerical simulation, and the dynamic properties and energy dissipation variation of the viscoelastic damper with self-heating processes are also analyzed. The results show that the dynamic properties of the viscoelastic damper are signiﬁcantly affected by temperature, excitation frequency and the internal self-generated heating.

In seismic structural design and retrofit engineering, viscoelastic dampers are often added to structures, working together to protect the building and reduce the damage caused by earthquakes. So far, a large number of investigations on dynamic properties of viscoelastic dampers and their applications in civil engineering have been conducted [7][8][9][10][11]. Asano et al. [12] studied the dynamic properties of the viscoelastic damper with different materials and investigated the influence of external loading circles, displacement amplitudes and excitation frequencies on the mechanical properties of the viscoelastic damper. Shen et al. [13] conducted a dynamic experimental investigation of the viscoelastic damper considering the frequency, strain amplitude and temperature effect and used the fractional derivative for modeling the dynamic behavior of the viscoelastic damper. Mazza et al. [14] used diagonal steel braces and viscoelastic dampers to retrofit a ten-story steel frame building with and without fire damages. A combination of two Maxwell models and one Kelvin model in parallel is employed to characterize the viscoelastic behavior of the damper. The dynamic response analysis of the structure shows that the viscoelastic dampers can effectively control the vibration of the fire-damaged steel structure after retrofitting. Christopoulos et al. [15] developed a viscoelastic coupling damper for application in vibration controls for tall buildings. It has been proved that the viscoelastic coupling damper can reliably reduce both dynamic responses caused by wind and earthquake for all loading amplitudes and frequencies. Xu et al. [16,17] utilized the viscoelastic damper to improve the seismic performance of the concrete frame structures. The design parameters and installation positions of the viscoelastic dampers are optimized, and shaking table tests of the structure models with and without dampers are carried out to validate the effectiveness of the dampers.
It can be seen that the mechanical properties and energy dissipation of the viscoelastic damper are greatly influenced by the excitation frequency, displacement amplitude and ambient temperature. It is important to select the appropriate mathematical model to describe the dynamic behavior of the viscoelastic damper in design and application procedures. Generally, the traditional viscoelastic mathematical models including the Kelvin model, Maxwell model, the generalized viscoelastic models, etc. are used to capture the properties of the viscoelastic materials [18][19][20]. Bagley et al. [21] and Meral et al. [22] employed the fractional calculus to model the constitutive relations of the viscoelastic materials. Gandhi et al. [23] used the standard linear solid model to simulate the nonlinear viscoelastic behaviors of the viscoelastic damper with dual-frequency excitations. Tsai et al. [7] proposed the finite element model which can consider the strain/displacement amplitude, frequency and ambient temperature influence on the mechanical properties of the viscoelastic dampers. Lewandowski et al. [24] utilized the generalized Kelvin and Maxwell model together with the W-L-F equation to characterize the properties of the viscoelastic damper considering the frequency and temperature influence. Xu et al. [25] theoretically and experimentally studied the temperature and displacement affection of viscoelastic dampers.
It also can be found that many scholars have paid attention to the influence of ambient temperature on the properties of viscoelastic dampers. However, the self-generated heating and temperature increment influence have received little attention [26][27][28][29][30]. Drozdov et al. [31] discussed the molecular chain fracture and deformation phenomenon of viscoelastic materials and studied the self-heating effects on dynamic properties of the viscoelastic material. Inaudi et al. [32] experimentally studied the mechanic properties change of the viscoelastic damper caused by heat generation during cyclic external loading, and the seismic responses of the structure with supplemental viscoelastic dampers considering self-generated heating are investigated by numerical simulation. Rodas et al. [33] formulated a thermo-visco-hyperelastic constitutive model to describe the self-heating process of viscoelastic materials with cyclic loadings. The finite-element-based simulations are conducted and compared with the experimental results.
In the present paper, dynamic properties tests of the sandwich viscoelastic damper at different environmental temperatures are carried out. The equivalent fractional Kelvin model is introduced to describe the dynamic properties and energy dissipation capability of the viscoelastic damper. The proposed model can characterize the mechanical behavior of the viscoelastic damper with varying frequencies and temperatures. The self-heating phenomenon of viscoelastic dampers is studied with numerical simulations. The dynamic properties and energy dissipation variations of the viscoelastic damper with different sizes and self-heating processes are analyzed. This demonstrates that the dynamic properties of the viscoelastic damper are significantly affected by temperature, excitation frequency and the internal self-generated heating.

Performance Tests
In order to better study the mechanical properties and energy dissipation of the sandwich viscoelastic damper with different temperatures, the dynamic properties tests of the sandwich viscoelastic damper at different temperatures (from 5 • C to 20 • C) are conducted. The results show that the viscoelastic damper has better performance of dynamic properties and energy dissipation capacities at low temperature. The dynamic parameters, including the storage modulus, loss factor, equivalent stiffness and equivalent damping, decrease with increasing ambient temperature.

Test Process
The sandwich viscoelastic damper consists of two viscoelastic material layers and three steel plates. The shear area of the viscoelastic damper is 60 mm × 50 mm, and the thickness of shear layer is 10 mm. The viscoelastic material is developed independently based on the nitrile rubber [11], and the steel type for the steel plates is Q235. The steel plates and viscoelastic layers of the damper are closely fixed with each other by chemical bonding in the vulcanization process and deform in different directions under cyclic external loading. The mechanical energy generated from external loading can be dissipated by the shear deformation of viscoelastic layers, and some parts of the mechanical energy are dissipated into the air in heat form. The schematic picture and photo of the viscoelastic damper are given in Figure 1a,b. properties and energy dissipation capacities at low temperature. The dynamic parameters, including the storage modulus, loss factor, equivalent stiffness and equivalent damping, decrease with increasing ambient temperature.

Test Process
The sandwich viscoelastic damper consists of two viscoelastic material layers and three steel plates. The shear area of the viscoelastic damper is 60 mm× , and the thickness of shear layer is 10 mm. The viscoelastic material is developed independently based on the nitrile rubber [11], and the steel type for the steel plates is Q235. The steel plates and viscoelastic layers of the damper are closely fixed with each other by chemical bonding in the vulcanization process and deform in different directions under cyclic external loading. The mechanical energy generated from external loading can be dissipated by the shear deformation of viscoelastic layers, and some parts of the mechanical energy are dissipated into the air in heat form. The schematic picture and photo of the viscoelastic damper are given in Figure 1a  To study the influence of ambient temperature on the dynamic properties and energy dissipation capability of the sandwich viscoelastic damper, the dynamic properties tests are conducted with a multipurpose universal servohydraulic LFV 100kN test series manufactured by w + b company, Switzerland, as seen in Figure 2a. Each test of the viscoelastic damper is carried out under 5 cycles of sinusoidal displacement loading = ( ), where represents the amplitude of the displacement loading, and denotes the loading frequency. The displacement and frequency in each test can be recorded directly in the control and data acquisition system with the DION 7 software. All the test conditions are given in Table 1. A temperature controlling device, as shown in Figure 2a  To study the influence of ambient temperature on the dynamic properties and energy dissipation capability of the sandwich viscoelastic damper, the dynamic properties tests are conducted with a multipurpose universal servohydraulic LFV 100kN test series manufactured by w + b company, Switzerland, as seen in Figure 2a. Each test of the viscoelastic damper is carried out under 5 cycles of sinusoidal displacement loading u d = u 0 sin(2π f t), where u 0 represents the amplitude of the displacement loading, and f denotes the loading frequency. The displacement and frequency in each test can be recorded directly in the control and data acquisition system with the DION 7 software. All the test conditions are given in Table 1. A temperature controlling device, as shown in Figure 2a,b, is used to adjust the environmental temperature of the damper.

Test Esults
In order to obtain the dynamic performance and energy dissipation indices of the sandwich viscoelastic damper, the fourth loop of the force-displacement curves at each test condition is selected and analyzed. The representative hysteresis curves of the sandwich viscoelastic damper are shown in Figure 3. Figure 3a,b presents the variations of the force-displacement hysteresis curves of the damper with changing ambient temperature (from −5 ℃ to 20 ℃) at the test conditions of displacement 1mm and frequency 0.2Hz and displacement 2mm and frequency 0.1 Hz, respectively. The sandwich viscoelastic damper has fuller hysteresis curves and larger elliptical slopes at low temperatures. As the temperature increases from −5 ℃ to 20 ℃, the fullness and slope of the hysteresis curves gradually decrease. Figure 3c,d gives the variations of the force-displacement curves of the damper with different loading frequencies (from 0.1Hz to 1.0Hz) with the test conditions of temperature −5 ℃ and displacement 1.0 mm, and temperature 15 ℃ and displacement 2.0 mm, respectively. The fullness and slope of the hysteresis curves gradually increase with increasing frequencies. As the fullness and slope of the hysteretic curves are greatly related to the stiffness and energy dissipation capacity of the viscoelastic damper, it can be concluded that the stiffness and energy dissipation capacity of the sandwich viscoelastic damper decrease with increasing temperature and increase with increasing frequency. The sandwich viscoelastic damper has better dynamic properties and energy dissipation ability at low temperature than that at high temperature.

Test Esults
In order to obtain the dynamic performance and energy dissipation indices of the sandwich viscoelastic damper, the fourth loop of the force-displacement curves at each test condition is selected and analyzed. The representative hysteresis curves of the sandwich viscoelastic damper are shown in Figure 3.  Under the sinusoidal displacement loading = 0 sin (2 ), the force-displacement hysteretic curve of the viscoelastic damper at each test condition can be treated as a full ellipse [8], as seen in Figure 4.  Figure 3a,b presents the variations of the force-displacement hysteresis curves of the damper with changing ambient temperature (from −5 • C to 20 • C) at the test conditions of displacement 1mm and frequency 0.2Hz and displacement 2mm and frequency 0.1 Hz, respectively. The sandwich viscoelastic damper has fuller hysteresis curves and larger elliptical slopes at low temperatures. As the temperature increases from −5 • C to 20 • C, the fullness and slope of the hysteresis curves gradually decrease. Figure 3c,d gives the variations of the force-displacement curves of the damper with different loading frequencies (from 0.1Hz to 1.0Hz) with the test conditions of temperature −5 • C and displacement 1.0 mm, and temperature 15 • C and displacement 2.0 mm, respectively. The fullness and slope of the hysteresis curves gradually increase with increasing frequencies. As the fullness and slope of the hysteretic curves are greatly related to the stiffness and energy dissipation capacity of the viscoelastic damper, it can be concluded that the stiffness and energy dissipation capacity of the sandwich viscoelastic damper decrease with increasing temperature and increase with increasing frequency. The sandwich viscoelastic damper has better dynamic properties and energy dissipation ability at low temperature than that at high temperature.
Under the sinusoidal displacement loading u d = u 0 sin(2π f t), the force-displacement hysteretic curve of the viscoelastic damper at each test condition can be treated as a full ellipse [8], as seen in Figure 4. Under the sinusoidal displacement loading = 0 sin (2  ), the force-displacement hysteretic curve of the viscoelastic damper at each test condition can be treated as a full ellipse [8], as seen in Figure 4.  The full ellipse can be expressed as In Equation (1), F d and u d denotes the damping force and corresponded displacement of the damper, respectively; K e = F 1 u o , which represents the equivalent stiffness, u o , and F 1 denotes the maximum displacement and the damping force at u o ; F 2 is the damping force of the damper when the displacement is zero, and F represents the biggest damping force of the viscoelastic damper on one hysteretic curve. Based on the hysteretic curve simplification of the viscoelastic damper, the dynamic properties and energy dissipation indices, such as the storage modulus G 1 , the loss factor η, the equivalent stiffness K e and equivalent damping C e of the viscoelastic damper can be expressed as [34] where n v denotes the number of viscoelastic layers, A v and h v represent the shear area and thickness of each viscoelastic layer, and ω = 2πf.
For the sandwich type viscoelastic damper utilized in this test, n v = 2, A v = 50 × 60 = 3000 mm 2 and h v =10 mm. Then, together with Equations (2)-(5), the storage modulus G 1 , the loss factor η, the equivalent stiffness K e and equivalent damping C e at each test condition can be calculated. To briefly show the dynamic properties and energy dissipation performance of the viscoelastic damper, only the dynamic parameters with displacement 1.0 mm are given and clearly pictured in Figure 5. The details of the dynamic parameters are shown in Table 2. This phenomenon indicates that the dynamic performance of the sandwich viscoelastic damper greatly decreases with increasing temperature, and the damper has better energy dissipation capacity at lower temperature. It can also be found that the storage modulus, loss factor, and equivalent stiffness increase when the excitation frequency increases, and equivalent damping decreases with the increasing frequency.   Figure 5 presents the characteristic parameters of the sandwich viscoelastic damper with changing temperature from −5 • C to 20 • C when the displacement amplitude is 1.0 mm. It is easy to see that the storage modulus, loss factor, equivalent stiffness and equivalent damping decrease significantly by increasing the temperature, especially at low temperatures (from −5 • C to 10 • C). Taking 1.0 Hz as an example, the storage modulus decreases by 15.91% in the range of −5 • C to 0 • C, decreases by 43.25% in the range of 0 • C to 10 • C, decreases by 12.14% in the range of 10 • C to 15 • C and decreases by 6.8% in the range of 15 • C to 20 • C. The loss factor decreases by 8.63% in the range of −5 • C to 0 • C, decreases by 34.86% in the range of 0 • C to 10 • C, decreases by 29.62% in the range of 10 • C to 15 • C, and decreases by 14.91% in the range of 15 • C to 20 • C. The equivalent stiffness decreases by 13.77% in the range of −5 • C to 0 • C, decreases by 44.66% in the range of Buildings 2021, 11, 470 7 of 16 0 • C to 10 • C, decreases by 12.14% in the range of 10 • C to 15 • C, and decreases by 6.8% in the range of 15 • C to 20 • C. The equivalent damping decreases by 23.17% in the range of −5 • C to 0 • C, decreases by 63.03% in the range of 0 • C to 10 • C, decreases by 38.17% in the range of 10 • C to 15 • C, and decreases by 20.7% in the range of 15 • C to 20 • C. This phenomenon indicates that the dynamic performance of the sandwich viscoelastic damper greatly decreases with increasing temperature, and the damper has better energy dissipation capacity at lower temperature. It can also be found that the storage modulus, loss factor, and equivalent stiffness increase when the excitation frequency increases, and equivalent damping decreases with the increasing frequency.

Mathematical Modeling
As seen in dynamic properties tests, the mechanical properties and energy dissipation capability of the sandwich viscoelastic damper are significantly affected by the ambient temperature and loading frequency. The equivalent fractional Kelvin model is employed in this section to describe the mechanical and damping properties variation of the sandwich viscoelastic damper with different ambient temperatures and loading frequencies. The fractional Kelvin model consists of one spring element and a fractional dashpot in parallel [25], as seen in Figure 6.

Mathematical Modeling
As seen in dynamic properties tests, the mechanical properties and energy dissipation capability of the sandwich viscoelastic damper are significantly affected by the ambient temperature and loading frequency. The equivalent fractional Kelvin model is employed in this section to describe the mechanical and damping properties variation of the sandwich viscoelastic damper with different ambient temperatures and loading frequencies. The fractional Kelvin model consists of one spring element and a fractional dashpot in parallel [25], as seen in Figure 6.
where 0 denotes the modulus of the spring element, and 1 denotes the damping of fractional dashpot; is the shear strain, and is the shear stress; is the fractional order derivative, and 0 < < 1.0. Applying the Fourier transform to Equation (6), we can  where q 0 denotes the modulus of the spring element, and q 1 denotes the damping of fractional dashpot; γ is the shear strain, and τ is the shear stress; α is the fractional order derivative, and 0 < α < 1.0. Applying the Fourier transform to Equation (6), we can obtain that where G * is the shear modulus of fractional Kelvin model in complex form, and j represents the unit complex number. By decomposing the complex modulus into two parts, the real part and the imaginary part, the storage modulus, shear modulus and loss factor of the viscoelastic material can be obtained as In the above equations, G 1 , G 2 and η represent the storage modulus, loss modulus and loss factor of viscoelastic materials, respectively.
When the ambient temperature is around T g to T g + 100 • C, the temperature and frequency influences on the dynamic properties of the viscoelastic damper have an equivalent relationship, and the dynamic properties of viscoelastic damper at high temperature are equivalent to those at low frequency, which can be described by the temperature-frequency equivalent theory [11,20,25], as shown in Equation (11) where T g is the glass transaction temperature, T 0 is the reference temperature, and α T can be written as From Equation (12), it can be seen that the temperature-frequency equivalent theory is an application of the W-L-F equation [24,35] in viscoelastic damping device studies. By considering the temperature-frequency equivalent theory together with Equations (8)-(10), the equivalent fractional kelvin model can be obtained as Equations (13) to (15) are expressions of the equivalent fractional Kelvin model. It can effectively characterize the affection of ambient temperatures and loading frequencies on dynamic properties and energy dissipation ability of viscoelastic materials and dampers.
To further reveal the temperature influence on dynamic properties and energy dissipation of the damper with equivalent fractional Kelvin model, the experimental results of the dynamic properties tests are used to determine the model parameters. As Equations (12)-(15) could not reflect the displacement amplitude influence and the fact that only the tests with displacement 1.0 mm and 2.0 mm are conducted, the test results with displacement 1.0mm and 2.0mm are treated separately. Only parts of the test data are selected for parameters determination. When the displacement is 1.0mm, the parameters of the equivalent fractional Kelvin model can be obtained as, q 0 = 6.69 × 10 5 , q 1 = 98.97, α = 0.4781, and T 0 = 490.52 K. When the displacement is 2.0 mm, q 0 = 6.2 × 10 5 , q 1 = 420.52, α = 0.4355, and T 0 = 479.02 K. The experimental and numerical results comparisons of the equivalent fractional Kelvin model of the sandwich viscoelastic damper with displacement 1.0 mm and 2.0 mm are plotted in Figure 7. The equivalent fractional Kelvin model can well characterize the mechanical properties and energy dissipation capacity of the sandwich viscoelastic damper with varying environmental temperatures. Taking a displacement of 2.0 mm as an example, the experimental and numerical data are given in Table 3. The numerical results agree well with the experimental results at different temperatures, and the errors do not exceed 21%. For storage modulus, the maximum error is 20.87% with an ambient temperature of −5 ℃ and a frequency 1.0 Hz; the experimental and numerical results are 2.45 MPa and 2.96 MPa, respectively. For loss factor, the maximum error is 19.75% with an ambient temperature of 20 ℃ and a frequency of 1.0 Hz; the experimental and numerical results are 0.314 and 0.376, respectively.  The equivalent fractional Kelvin model can well characterize the mechanical properties and energy dissipation capacity of the sandwich viscoelastic damper with varying environmental temperatures. Taking a displacement of 2.0 mm as an example, the experimental and numerical data are given in Table 3. The numerical results agree well with the experimental results at different temperatures, and the errors do not exceed 21%. For storage modulus, the maximum error is 20.87% with an ambient temperature of −5 • C and a frequency 1.0 Hz; the experimental and numerical results are 2.45 MPa and 2.96 MPa, respectively. For loss factor, the maximum error is 19.75% with an ambient temperature of 20 • C and a frequency of 1.0 Hz; the experimental and numerical results are 0.314 and 0.376, respectively.

Self-Generated Heating of the Sandwich Viscoelastic Damper
The dynamic properties and energy dissipation capability of viscoelastic materials are significantly affected by ambient temperature. In the design and seismic analysis of building structures installed with viscoelastic dampers, the influence of ambient temperature is always ignored. During the service life of the viscoelastic damper, the ambient temperature is not stable. At the same time, the external mechanical energy could be converted into heat through the shear deformation of the viscoelastic layers, which causes the internal temperature to rise and affects the dynamic performance of the damper. In this section, the self-heating phenomenon of the sandwich viscoelastic damper is investigated. The initial temperature of the whole damper and ambient temperature is taken as 10 • C.
Considering that the heat generation rate of the viscoelastic material is affected by ambient temperature and ignoring the thermoelastic coupling effect (thermal expansion), the heat transfer equation of the viscoelastic layer in the damper can be written as [31] q g + k∇ 2 T(t, x, y, z) = ρc p . T(t, x, y, z) (16) with the boundary conditions where, k, ρ, and c p represent the thermal conductivity, density and specific heat capacity of the materials, respectively. The value of k∇ 2 T(t, x, y, z) can be obtained by applying Fourier's law to a certain volume; ρc p . T(t, x, y, z) represents the heat restored inside the material. T 0 denotes the temperature loading on the boundary of the domain Ω D , q 0 denotes the heat flow on the boundary domain Ω N , h represents the heat convection coefficient between the boundary domain Ω C and the surroundings, T ∞ is the ambient temperature, q g is the heat generation rate of the viscoelastic material and has the form where β t is the thermal conversion ratio, which means the fraction of heat energy converted from the dissipated mechanical energy through shear deformation of the viscoelastic materials, . w m has the form [31] . w m (t, ω, T) = −ωG 2 (ω, T)Cε T 0 (t, ω, T)sin 2 (ωt + δ) (19) where C(t, ω, T) = G(t, ω, T)C, G(t, ω, T) = G 1 (t, ω, T) + G 2 (t, ω, T)j; ε 0 denotes the strain of the viscoelastic materials, and δ presents the phrase angle of the stress and strain. Based on the abovementioned heat generate and transfer theories, the finite-element simulation with a two-dimensional (2-D) model is conducted to study the self-heating phenomenon. The 2-D model of the viscoelastic damper is shown in Figure 8. Based on the abovementioned heat generate and transfer theories, the finite-element simulation with a two-dimensional (2-D) model is conducted to study the self-heating phenomenon. The 2-D model of the viscoelastic damper is shown in Figure 8. In the finite element simulation, the displacement amplitude of the viscoelastic damper is 1.0 mm, and the 'Coupled-Temperature-Displacement' time step is employed to simulate the thermo-viscoelastic behaviors of the material. The 'CPE4T'planar element type is used for meshing of viscoelastic damper. In order to avoid the element distortion with a large number of loading cycles, the element size of the viscoelastic layer is selected as 0.1 mm, while the element size of the steel plate is getting larger and larger when it is gradually getting away from the viscoelastic materials. The total number of finite elements of the 2-D model is 255428. In order to better understand the temperature rising at different positions of the viscoelastic layer, four reference points 'A', 'B', 'C' and 'D' are selected in the viscoelastic layer, as seen in Figure 8b. A local coordinate system is defined in the upper viscoelastic layer, and the coordinates of the reference points 'A', 'B', 'C' and 'D' are (30,7.5), (30,5.0), (30,2.5) and (55,5.0), respectively.
In this section, the temperature influence on the material properties should be considered to simulate the thermo-viscoelastic behavior of the viscoelastic damper. In ABAQUS, the Prony series by default defines the viscoelastic properties of the viscoelastic material. The W-L-F equation = − 1 ( − 0 )( 2 + − 0 ) [24,35] can be utilized to describe the temperature influence. 1 and 2 are material constants. With the Prony series and W-L-F equation, the viscoelastic properties of the viscoelastic material can be obtained and shown in Table 4 [35]. The Poisson's ratio of the viscoelastic material is 0.49, and the modulus and Poisson's ratio of the steel plate is 210GPa and 0.3. The thermal properties of the viscoelastic layer and steel plate are shown in Table 5 [36,37].  In the finite element simulation, the displacement amplitude of the viscoelastic damper is 1.0 mm, and the 'Coupled-Temperature-Displacement' time step is employed to simulate the thermo-viscoelastic behaviors of the material. The 'CPE4T'planar element type is used for meshing of viscoelastic damper. In order to avoid the element distortion with a large number of loading cycles, the element size of the viscoelastic layer is selected as 0.1 mm, while the element size of the steel plate is getting larger and larger when it is gradually getting away from the viscoelastic materials. The total number of finite elements of the 2-D model is 255428. In order to better understand the temperature rising at different positions of the viscoelastic layer, four reference points 'A', 'B', 'C' and 'D' are selected in the viscoelastic layer, as seen in Figure 8b. A local coordinate system is defined in the upper viscoelastic layer, and the coordinates of the reference points 'A', 'B', 'C' and 'D' are (30,7.5), (30,5.0), (30,2.5) and (55,5.0), respectively.
In this section, the temperature influence on the material properties should be considered to simulate the thermo-viscoelastic behavior of the viscoelastic damper. In ABAQUS, the Prony series by default defines the viscoelastic properties of the viscoelastic material. [24,35] can be utilized to describe the temperature influence. C 1 and C 2 are material constants. With the Prony series and W-L-F equation, the viscoelastic properties of the viscoelastic material can be obtained and shown in Table 4 [35]. The Poisson's ratio of the viscoelastic material is 0.49, and the modulus and Poisson's ratio of the steel plate is 210GPa and 0.3. The thermal properties of the viscoelastic layer and steel plate are shown in Table 5 [36,37]. Table 4. Viscoelastic parameters determined with Prony series and W-L-F equation in ABAQUS (n = 4). when d = 1.0 mm.
13.104 486.9721 0.3485 0.0357 0.0546 0.502 1.9055 × 10 −11 7.0947 × 10 −8 4.3061 × 10 −9 2.5573 × 10 −10 12 525 Table 5. Thermal properties parameters of the viscoelastic layers and steel plates. The heat convection coefficient h between the surfaces of the viscoelastic damper and surrounding environment is taken as 6 w·m 2 /K [38]. The model of the viscoelastic damper is subjected to 200 cycles of displacement loading u d = u 0 sin(2π f t). When the displacement amplitude is 1.0 mm and the loading frequency is 5.0 Hz, the temperature contours of the viscoelastic damper at different loading circles are shown in Figure 9. The heat convection coefficient ℎ between the surfaces of the viscoelastic damper and surrounding environment is taken as 6 ( • 2 ⁄ ) [38]. The model of the viscoelastic damper is subjected to 200 cycles of displacement loading = 0 (2 ). When the displacement amplitude is 1.0 mm and the loading frequency is 5.0 Hz, the temperature contours of the viscoelastic damper at different loading circles are shown in Figure 9.
The temperature of the viscoelastic materials increases gradually as the number of loading circles increases. The temperature of the upper and lower layers of the viscoelastic material is symmetrically distributed with respect to the intermediate steel plate. A higher temperature area is formed within a certain region inside the viscoelastic layers. After 200 cycles of sinusoidal displacement loading, the local temperature at the center of the viscoelastic layer reaches to 14.52 ℃, and the maximum increment is 4.52 ℃.   Figure 10 shows the temperature increment at the reference points with increasing loading circle numbers. The temperature rising processes at reference points 'A', 'B', 'C' and 'D' are different, which is caused by the difference in heat generation and transfer due to the change of position and boundary conditions. With the same loading circle, the temperature increment increases when the loading frequency goes up. The temperature with larger frequency is obviously higher than that with lower frequency. The reason is that the mechanical properties such as the modulus of the viscoelastic material increases with increasing frequency. The mechanical energy dissipation and heat generation rates in each single loading cycle of the viscoelastic material increase with increasing frequency. Take reference point 'B' as an example, at the 140th loading circle, the temperature is 11.06, 13.27, 14.88 and 16.08 ℃ with frequency 1.0, 5.0, 10.0 and 20.0 Hz. The temperature increment compared to initial temperature is 1.06, 3.27, 4.88 and 6.08 ℃, respectively. The temperature of the viscoelastic materials increases gradually as the number of loading circles increases. The temperature of the upper and lower layers of the viscoelastic material is symmetrically distributed with respect to the intermediate steel plate. A higher temperature area is formed within a certain region inside the viscoelastic layers. After 200 cycles of sinusoidal displacement loading, the local temperature at the center of the viscoelastic layer reaches to 14.52 • C, and the maximum increment is 4.52 • C. Figure 10 shows the temperature increment at the reference points with increasing loading circle numbers. The temperature rising processes at reference points 'A', 'B', 'C' and 'D' are different, which is caused by the difference in heat generation and transfer due to the change of position and boundary conditions. With the same loading circle, the temperature increment increases when the loading frequency goes up. The temperature with larger frequency is obviously higher than that with lower frequency. The reason is that the mechanical properties such as the modulus of the viscoelastic material increases with increasing frequency. The mechanical energy dissipation and heat generation rates in each single loading cycle of the viscoelastic material increase with increasing frequency. Take reference point 'B' as an example, at the 140th loading circle, the temperature is 11.06, 13.27, 14.88 and 16.08 • C with frequency 1.0, 5.0, 10.0 and 20.0 Hz. The temperature increment compared to initial temperature is 1.06, 3.27, 4.88 and 6.08 • C, respectively.  As the mechanical properties and energy dissipation capacity of the viscoelastic damper are greatly affected by ambient temperature, it is necessary to study the self-heating influence on the dynamic behaviors of the viscoelastic damper. Figure 11 gives the force-displacement hysteretic curve, damping force, storage modulus and loss factor variations of the sandwich viscoelastic damper at frequencies 1.0 , 5.0 , 10.0 and 20.0 Hz with loading circle numbers within 100. Figure 11a,b reveals that the fullness and slope of the force-displacement hysteretic curves change apparently when the loading circles increase. The damping force with circles from 95 to 100 is much less than that with circles from 1 to 5, which can be attributed to the modulus and stiffness decrease caused by the self-heating inside the viscoelastic layer. Taking the 5th and 95th loading circles at 5.0 Hz as an example, the peak force at unit width of the viscoelastic damper is 38.8438N and 34.7891N, respectively. The peak force at the 95th circle is 10.44% lower than that at the 5th circle. As the mechanical properties and energy dissipation capacity of the viscoelastic damper are greatly affected by ambient temperature, it is necessary to study the selfheating influence on the dynamic behaviors of the viscoelastic damper. Figure 11 gives the force-displacement hysteretic curve, damping force, storage modulus and loss factor variations of the sandwich viscoelastic damper at frequencies 1.0, 5.0, 10.0 and 20.0 Hz with loading circle numbers within 100. Figure 11a,b reveals that the fullness and slope of the force-displacement hysteretic curves change apparently when the loading circles increase. The damping force with circles from 95 to 100 is much less than that with circles from 1 to 5, which can be attributed to the modulus and stiffness decrease caused by the self-heating inside the viscoelastic layer. Taking the 5th and 95th loading circles at 5.0 Hz as an example, the peak force at unit width of the viscoelastic damper is 38.8438N and 34.7891N, respectively. The peak force at the 95th circle is 10.44% lower than that at the 5th circle. The FEA results of the storage modulus and loss factor of the viscoelastic damper are shown in Figure 11 c,d, and the storage modulus decreases when the loading cycles increase. The loss factor decreases as the number of loading circle increases at low frequencies (1.0 Hz and 5.0 Hz); and increases with the increasing loading cycles at high frequencies (10.0 Hz and 20.0 Hz). Taking the 10th, 30th, 50th and 80th circles with 10 Hz as examples, the FEA results of storage modulus are 3.6793, 3.4528, 3.2663 and 3.0417 MPa, respectively. It decreases by 6.16% in the range of the 10th to the 30th loading circles, decreases by 5.4% in the range of the 30th to the 50th loading circles and decreases by 6.88% in the range of the 50th to the 80th loading circles.

Conclusions
In the present work, dynamic properties tests of the sandwich viscoelastic damper at different ambient temperatures were conducted. The equivalent fractional Kelvin model was introduced to describe the thermo-viscoelastic behaviors of the viscoelastic damper. The self-heating phenomenon of the viscoelastic damper was investigated, and the finite element simulation was utilized to verify the self-heating influence on dynamic properties and energy dissipation capacity of the sandwich type viscoelastic damper.
Based on these studies some notable conclusions can be obtained as follows: ① The sandwich viscoelastic damper has perfect energy dissipation capacity, and the performance of the viscoelastic damper at low temperature is better than that at high temperature. Almost all the hysteretic curves of the sandwich viscoelastic damper are full ellipses, and the key energy dissipation parameter, loss factor η, reaches its maximum value near The FEA results of the storage modulus and loss factor of the viscoelastic damper are shown in Figure 11c,d, and the storage modulus decreases when the loading cycles increase. The loss factor decreases as the number of loading circle increases at low frequencies (1.0 Hz and 5.0 Hz); and increases with the increasing loading cycles at high frequencies (10.0 Hz and 20.0 Hz). Taking the 10th, 30th, 50th and 80th circles with 10 Hz as examples, the FEA results of storage modulus are 3.6793, 3.4528, 3.2663 and 3.0417 MPa, respectively. It decreases by 6.16% in the range of the 10th to the 30th loading circles, decreases by 5.4% in the range of the 30th to the 50th loading circles and decreases by 6.88% in the range of the 50th to the 80th loading circles.

Conclusions
In the present work, dynamic properties tests of the sandwich viscoelastic damper at different ambient temperatures were conducted. The equivalent fractional Kelvin model was introduced to describe the thermo-viscoelastic behaviors of the viscoelastic damper. The self-heating phenomenon of the viscoelastic damper was investigated, and the finite element simulation was utilized to verify the self-heating influence on dynamic properties and energy dissipation capacity of the sandwich type viscoelastic damper.
Based on these studies some notable conclusions can be obtained as follows: 1 The sandwich viscoelastic damper has perfect energy dissipation capacity, and the performance of the viscoelastic damper at low temperature is better than that at high temperature. Almost all the hysteretic curves of the sandwich viscoelastic damper are full ellipses, and the key energy dissipation parameter, loss factor η, reaches its maximum value near 1.0 with a temperature of −5 • C and a frequency of 1.0 Hz. 2 The storage modulus, loss factor, equivalent stiffness and equivalent damping decrease when the temperature increases, and the storage modulus, loss factor, and equivalent stiffness increase with an increase in frequency, while the equivalent damping decreases with an increase in frequency. 3 The temperature inside the viscoelastic layers could be obviously increased during the cyclic displacement loading due to the self-generated heating. The parameters decrease importantly because of the internal temperature increment. When the displacement amplitude is 1.0 mm and the loading frequency is 5.0 Hz, the maximum temperature increment is 4.52 • C after 200 cyclic loadings with initial temperature 10 • C. When displacement is 1mm, and frequency is 10 Hz, the storage modulus changes from 3.6793 to 3.0417 MPa at the 10th and the 80th loading circles, decreasing by 17.33%.
In this paper, dynamic properties tests of the sandwich viscoelastic damper with different temperatures were carried out. The equivalent fractional order kelvin model was formulated to capture the temperature influence on dynamic properties. The Finite Element Analysis method is used to investigate the self-heating phenomena. The results reveal that the dynamic properties of the damper are significantly affected by temperature, excitation frequency and the internal self-generated heating.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author. The data are not publicly available due to the confidentiality requirements of the lab.