Abstract
Earthquake excitation is a kind of complicated multi-dimensional ground motion. Therefore, anti-earthquake measuresfor complex and important engineering structures, such as long-span structures and high rise buildings under multi-dimensional earthquake excitations, have become an extremely important research topic. In this paper, a new multi-dimensional earthquake isolation and mitigation device (MEIMD) is proposed for protecting long-span reticulated structures due to strong earthquake motions. It consists of a viscoelastic bearing and several viscoelastic dampers. Property tests on the MEIMD under different excitation frequencies, excitation amplitudes and vertical pressure forces are carried out. Then, a mathematical simulation is carried out for the equivalent stiffness and the equivalent damping of the device, and the numerical results agree reasonably well with the experimental data. Experimental and numerical results show that the MEIMD has a fine energy dissipation capacity, and can be used for anti-earthquake purposes for long-span reticulated structures.
1. Introduction
Conventional seismic design cannot usually ensure that some important structures are fully functional to avoid great economic losses and seismic casualties during a major earthquake (e.g., tall buildings, and long-span structures). Since 1972, when the vibration control theory was first proposed (Yao, 1972), the idea that reducing the seismic responses of structures by using vibration control devices has been widely accepted by structural engineers. The seismic isolation system is considered as the most effective method of decoupling seismic responses of structures. For example, at present, laminated rubber bearing has been extensively used in seismic isolation structures and it performs well (Ibrahim, 2008). However, conventional laminated rubber bearing can dissipate little earthquake energy during an earthquake (Narasimhan et al., 2006). Therefore new types of isolation systems with energy dissipation capacity were studied, such as laminated rubber bearing with a lead-core, frictional type sliding isolators and hybrid seismic isolation systems (Kasai et al., 1998; Roussis and Constantinou, 2006; Shrimali and Jangid, 2003). These isolation and mitigation devices are mostly designed to reduce only horizontal earthquake responses of structures. In fact, the earthquake excitations are complicated multi-dimensional ground motions, which include both horizontal and vertical components. Vertical earthquake-induced ground motion must be also be considered for important or complex structures, such as long-span reticulated structures and high-rise buildings. However, nowadays few devices have been studied and designed for reducing both horizontal and vertical earthquake responses of structures simultaneously. Laminated rubber bearings and viscous dampers were used to isolate the horizontal earthquake for a nuclear power station, while helical spring and viscous dampers were designed to control the vertical earthquake component (Huang et al., 2007). A new three-dimension isolation device composed of a rubber isolation device and sealed aircushion was proposed, and tests indicated that the device performed well in both horizontal and vertical directions (Suhara et al., 2002). On the basis of Suhara’s research results, a hydraulic cylinder and sealed aircushion were used to isolate a vertical earthquake while a rubber isolation device was used to reduce the horizontal response (Shimada et al., 2005). A device using a butterfly spring, rubber shock isolation blankets and MR damper was proposed to reduce the multi-dimensional earthquake responses (Qu, 2003).
Here, a new multi-dimensional earthquake isolation and mitigation device (MEIMD) which has not only earthquake isolation ability but also earthquake mitigation ability is proposed to reduce both horizontal and vertical earthquake responses of long-span reticulated structures simultaneously (Approved National Patent in China with granted number: ZL2006 1 0097219.3). The device makes full use of high energy absorption capacity of viscoelastic material. The excellent performance of passive control using a viscoelastic damper has been proved by both shaking table tests (Chang et al., 1995; Hayes et al., 1999; Xu, 2007) and real applications under strong earthquakes (Aiken et al., 1990; Crosby et al., 1994; Soong and Shen, 1995). The device consists of a viscoelastic bearing and several viscoelastic dampers, which can isolate and dissipate both horizontal and vertical vibration energy of structures. Property tests on the MEIMD under different excitation frequencies, excitation amplitudes and vertical pressure forces are carried out. Then, a mathematical simulation is carried out to simulate the equivalent stiffness and the equivalent damping of the device. The numerical results agree reasonably well with the experimental data. Experimental and numerical results show that the MEIMD has a fine energy dissipation capacity, and can be used for anti-earthquake purposes for long-span reticulated structures.
2. Test setup
2.1. Description of the multi-dimensional earthquake isolation and mitigationdevice (MEIMD)
The newly proposed MEIMD for reducing the dynamic responses of long-span reticulated structures is shown in Figure 1. The device consists of a viscoelastic core bearing and two viscoelastic dampers. For the horizontal earthquake excitation, the viscoelastic core bearing can isolate the transfer of vibration energy. At the same time, the viscoelastic core bearing, under vertical pressure due to the weight of the upper structure and two viscoelastic dampers, will dissipate vibration energy because of the high energy dissipation capacity of the viscoelastic material. For the vertical earthquake excitation, viscoelastic bearing and dampers also form a stiffness softening layer, i.e. isolation layer, which can isolate the transfer of vibration energy in a vertical direction.
Diagram and photograph of the multi-dimensional earthquake isolation and mitigation device.
Simultaneously, the viscoelastic dampers and bearing can dissipate vibration energy when the device produces a relative displacement. The energy dissipation is small under a minor earthquake due to a small relative vertical displacement of the device, while it is large under a strong earthquake which will lead to obvious vertical tension displacement of the device. If the MEIMD is tensioned, the upper steel plate will detach from the core bearing, therefore the viscoelastic dampers will dissipate most of the vibration energy, and the core bearing will be prevented from failure caused by tension. The tension failure is the important shortcoming of the traditional laminated rubber bearing. If the MEIMD is pressed, the pressure will be undertaken by the core bearing and the viscoelastic dampers. This is sufficient to carry the load of a superstructure just like the traditional laminated rubber bearing. It must be noted that the device cannot only permit a large deformation in a horizontal direction, but should also permit large tension deformation in a vertical direction. Above all, the viscoelastic damper and viscoelastic core bearing can work together to reduce the seismic responses effectively under both horizontal and vertical directions.
Detail design parameters of the MEIMD
MEIMD: multi-dimensional earthquake isolation and mitigation device.
2.2. Test procedure
The most important mechanical characteristics of the MEIMD are force-displacement energy dissipation hysteresis loop, the equivalent stiffness, and the equivalent damping. In this study, the effects of excitation frequency, excitation amplitude and vertical pressure force on these characteristics of the MEIMD under horizontal direction are considered. The horizontal performance tests have been done in Key Laboratory of C&PC Structures of the Ministry of Education, People’s Republic of China. As shown in Figure 2, the MEIMD is mounted to a 30 ton MTS closed-loop feedback hydraulic testing machine which works through a displacement control mode. The test system can immediately return the actual output force, displacement, velocity and acceleration of the actuator and can show the force-displacement behavior simultaneously. The MEIMD is subjected to ten cycles of sinusoidal excitation with constant displacement amplitude, excitation frequency and vertical pressure load. The vertical load on the MEIMD is provided by a 50 ton oil jack. Outputs from the actuator load cell and displacement transducer are recorded by a high speed digital data acquisition system. The corresponding high speed digital data acquisition system is provided by MTS Systems Corporation in the USA, which includes a high speed computer server and some corresponding software. The system can immediately return actual output force, displacement, velocity and acceleration data, and can plot the force-displacement hysteresis curves simultaneously. Tests are conducted at displacements and frequencies selected to be appropriate for seismic applications. Tests are carried out under the ambient temperature 22°C.
Loading device of the multi-dimensional earthquake isolation and mitigation device in horizontal performance test.
Horizontal property test conditions of the MEIMD
MEIMD: multi-dimensional earthquake isolation and mitigation device.
3. Tests results and analysis
During the entire process of the tests, the MEIMD isnot damaged. When the vertical load is 0 ton, theviscoelastic core bearing does not deform at all, andonly the dampers cause shear deformation and dissipate energy. While the vertical load is applied onthe MEIMD, the shear deformation of the core bearing can be found because the horizontal shear force can be transmitted to the bearing by the friction force between the upper steel plate and the core bearing.
The force-displacement hysteresis curves are plotted to understand thoroughly the energy dissipation capacity and the horizontal seismic performance of the proposed MEIMD under different excitation frequencies and excitation amplitudes without the vertical pressure force, as shown in Figures 3 and 4. Figure 3 shows force-displacement hysteresis curves under different displacement amplitudes at the excitation frequencies of 0.1 Hz, 0.2 Hz, 0.5 Hz and 1.0 Hz. The hysteresis loops shown in Figure 3 are smooth and plump ellipses, and this means that the MEIMD has fine energy dissipation ability. With the increase of the excitation amplitude, the damping force amplitude and the area of the hysteresis curve both increase, while the obliquity of the center inclined line decreases, and this decrease of the obliquity means a decrease of the equivalent stiffness of the MEIMD. These characteristic of hysteresis curves are the same under the frequencies of 0.1 Hz, 0.2 Hz, 0.5 Hz and 1 Hz. It is also found that the real displacement amplitude of the MEIMD cannot reach the required values under the frequency of 1 Hz, which is due to the limitation of the output power of servo-hydraulic test machine (MTS) under high frequency. Figure 4 shows force-displacement hysteresis curves under different excitation frequencies at the excitation amplitudes of 2 mm, 4 mm, 6 mm and 8 mm. It can be seen that the obliquity of the center inclined line increases with the increase of excitation frequency, which means that there is an increase in equivalent stiffness of the MEIMD when the excitation frequency increases. In order to acquire the effect of the vertical pressure on properties of the MEIMD, the force-displacement hysteresis curves under different pressure loads at the excitation amplitudes of 2 mm, 4 mm, 6 mm and 8 mm and the excitation frequency of 0.1 Hz are plotted as Figure 5. As seen, vertical pressure loads have no obvious effect on properties of theMEIMD when the excitation amplitudes are 2 mm, 4 mm, 6 mm, and the inclination of the hysteresis curves increases slightly with the increasing vertical pressure loads at the excitation amplitude of 8 mm.
Force-displacement hysteretic loops under different displacement. Force-displacement hysteretic loops under different frequency. Force-displacement hysteretic loops under different vertical load.


As mentioned above, the equivalent stiffness Ke and the equivalent damping Ce are the most important factors that reflect the horizontal characteristics of the MEIMD. The two factors can be obtained from the tested hysteresis curves. Under the sinusoidal excitation of Force-displacement curve of a damper. The equivalent stiffness and damping of the MEIMD without vertical loads MEIMD: multi-dimensional earthquake isolation and mitigation device.

3.1. Excitation frequency effect
In order to reveal the effect of the excitation frequency on the equivalent stiffness and the equivalent damping of the device, the calculated data listed in Table 3 are plotted in Figure 7. It can be seen from Figure 7(a) that the equivalent stiffness Ke increases slowly with increasing frequency. The equivalent stiffness Ke mainly comes from the stiffness of two viscoelastic plate dampers when the pressure load is 0 ton. Take a 4 mm excitation displacement amplitude condition for example, when the frequency increases from 0.1 Hz to 0.2 Hz, 0.2 Hz to 0.5 Hz, 0.5 Hz to1 Hz, and by every 0.1 Hz increasing frequency, the corresponding equivalent stiffness Keis increased by 7.3%, 4.1%, 3%, respectively, and the detailed data can be seen in Table 3. Adopting thesame data disposal measures, under the excitation displacement amplitude of 4 mm and the vertical pressure load of 30 ton, and by every 0.1 Hz of increasing frequency, the corresponding increase percentages of the equivalent stiffness Ke are 6.4%, 5.8% and 3.2%. The equivalent stiffness is increased with more frequency which is due to the storage modulus of viscoelastic materials increasing with more frequency. The storage modulus and loss factor are the main indices of properties for viscoelastic materials, and they mean the stored energy representing the elastic portion and the energy dissipated as heat representing the viscous portion, respectively.
Influence of excitation frequency on properties of the multi-dimensional earthquake isolation and mitigation device.
Variation of the equivalent damping Cewith excitation frequency can be seen in Figure 7(b). It can be found that the equivalent damping Ce decreases sharply when the frequency changes from 0.1 Hz to 0.2 Hz, and it decreases slowly when the frequency changes from 0.2 Hz to 1 Hz. Take 4 mm excitation displacement amplitude condition for example, when the frequency is increased from 0.1 Hz to 0.2 Hz, 0.2 Hz to 0.5 Hz, 0.5 Hz to 1 Hz, and by every 0.1 Hz increasing frequency, the corresponding equivalent damping Ce is decreased by 44.9%, 16.2%, 12.7%, respectively. In the same way, the equivalent damping under the excitation displacement amplitude of 4 mm and the vertical pressure load of 30 ton, and by every 0.1 Hz increasing frequency, the corresponding decrease percentages of the equivalent damping Ce are 42.5%, 16.7%, 8.6%. Decrease of the damping coefficient Ce with increasing frequency is primarily caused for two reasons, one is the energy dissipation per cycle will decrease slightly in higher frequency because of the power limitation of MTS, and the other is that the loss factor of viscoelasticmaterials decreases with increasing frequency. Inaddition, as shown in Figure 7(b), the excitation frequency has a more significant effect on the equivalent damping than high frequency, because the equivalent damping Ce and the frequency ω are inversely proportional according to equation (3), so the equivalent damping is more sensitive at low frequency.
3.2. Excitation displacement amplitude effect
For the purpose of showing the excitation displacement amplitude effect on the equivalent stiffness and damping of the device, the data listed in Table 3, obtained from the tests, are plotted in Figure 8. It can be seen in Figure 8(a) that the equivalent stiffness Ke decreases slowly with increasing displacement amplitude. Take 1.0 Hz excitation frequency as an example, when the displacement amplitude increases from 2 mm to 4 mm, 4 mm to 6 mm, 6 mm to 8 mm, the corresponding equivalent stiffness Ke decreases by 10.6%, 13.2%, and 28.2%, respectively. When the excitation frequency is 1.0 Hz and the vertical load is 30 ton, the corresponding decreasing percentages of the equivalent stiffness Ke are 10.3%, 11.7% and 9.4%. This phenomenon of the decreasing equivalent stiffness is due to the fact that the energy dissipation per cycle increases with increasing displacement amplitude. This dissipation energy transfers into the heat energy of the viscoelastic material, which leads to a temperature rise. This rise of temperature will in turn cause a decrease in the storage modulus of the viscoelastic material, and the decrease in storage leads to the decrease in equivalent stiffness Ke. Variation of the equivalent damping Ce with excitation amplitude can be seen in Figure 8(b). It can be found that the equivalent damping coefficient Ce decreases slowly with an increasing displacement amplitude.
Influence of displacement amplitude on equivalent stiffness and equivalent damping of the multi-dimensional earthquake isolation and mitigation device.
Take 0.2 Hz excitation frequency condition as an example, when the displacement amplitude is increased from 2 mm to 4 mm, 4 mm to 6 mm, 6 mm to 8 mm, the corresponding decrease percentages of the equivalent damping Ceare 10.5%, 14.2% and 26.3%, respectively. In the same way, considering the equivalent damping under the excitation frequency of 0.2 Hz and the vertical load of 30 ton, the corresponding decrease percentages of the equivalent damping Ce are 11.4%, 9.2% and 7.4%. Decrease of the damping with increasing amplitude is also caused by the temperature rise of viscoelastic material under large amplitudes. The transferred heat during vibration under large amplitudes in viscoelastic material is difficult to be radiated into the air space very quickly, which will lead to a rise in the temperature of viscoelastic material and cause the loss factor of viscoelastic materials decreases with increasing frequency.
3.3. Vertical pressure load effect
In order to reveal the vertical pressure load effect on the equivalent stiffness and damping of the device, the data obtained from tests are plotted in Figure 9. It can be seen from Figure 9(a) that the influence of vertical load on Ke is not as obvious as the excitation frequency anddisplacement amplitude. The equivalent stiffnessKe increases slowly, and the increasing percentage decreases with an increasing vertical load. Take 0.2 Hz excitation frequency and 6 mm displacement amplitude for example, when the vertical load increases from 0 ton to 10 ton, 10 ton to 20 ton and 20 ton to 30 ton, the corresponding equivalent stiffness is increased by 2.26% and 0.08% and 1.8%. In the same way, when the equivalent stiffness under the excitation frequency of 1 Hz and the excitation amplitude of 6 mm, the corresponding increase percentages of the equivalent stiffness are 5.4%, 1.1% and -0.4%, respectively. The equivalent stiffness under the excitation frequency of 0.5 Hz and the excitation amplitude of 8 mm are 7027, 8664, 8975, 9225 kN/m. When the vertical load increases from 0 ton to 10 ton, 10 ton to 20 ton and 20 ton to 30 ton, the corresponding equivalent stiffness is increased by 23.3%, 3.6% and 2.8% respectively. The law of variation of the equivalent damping Cewith the vertical load is plotted in Figure 9(b). It can be found that the effect of vertical load on Ce is also not as obvious as the effects of the excitation frequency and displacement amplitude. The equivalent damping coefficient Ce increases slowly and the increasing percentage decreases with an increasing vertical load. Take 0.2 Hz excitation frequency and 6 mm displacement amplitude condition for example, when the vertical load increases from 0 t to 10 t, 10 t to 20 t and 20 t to 30 t, the corresponding equivalent damping coefficient is increased by 20.1%, 3.6%, 2.98%. The equivalent damping under the excitation frequency of 0.5 Hz and the excitation amplitude of 8 mm are 625, 869, 916, 916 kN…s/m. When the vertical load increases from 0 ton to 10 ton, 10 ton to 20 ton and 20 ton to 30 ton, the corresponding equivalent stiffness is increased by 39.0%, 5.4% and 0%, respectively.
Influence of vertical load on equivalent stiffness and equivalent damping of the multi-dimensional earthquake isolation and mitigation device.
4. Numerical analysis
4.1. Mathematical simulation
As mentioned before, the MEIMD consists of a viscoelastic core bearing and several viscoelastic dampers, which makes the device a newly complex and nonlinear system.
It is necessary to develop mathematical equations to simulate the characteristics of the MEIMD. The horizontal characteristic parameters of the MEIMD include the equivalent stiffness Ke and equivalent damping Ce. As known from the above tests, the excitation frequency and displacement amplitude have influence on Ke and Ce, while the vertical load has a slight influence. Therefore, when developing mathematical equations only the excitation frequency and excitation displacement amplitude are considered. It can be seen from experimental data that the equivalent stiffness increases with increasing excitation frequency or decreasing displacement amplitude, and the equivalent damping decreases with increasing excitation frequency or displacement amplitude. The equivalent stiffness Ke andthe equivalent damping Ce can be expressed as follows
4.2. Comparison analysis
Comparison of experimental and numerical values
MEIMD: multi-dimensional earthquake isolation and mitigation device.
Figure 10 shows the experimental and numerical results comparison of the MEIMD under different frequencies with the fixed displacement amplitude of6 mm.
Influence of excitation frequency on equivalent stiffness and equivalent damping of the multi-dimensional earthquake isolation and mitigation device.
This group of experimental data is not included in those fitting data for obtaining equations (6) and (7). The experimental data for the equivalent stiffness Ke are 8352, 8504, 9604, 11092 kN/m under frequencies of 0.1, 0.2, 0.5, 1.0 Hz excitations, and the corresponding numerical data are 8259, 8593, 9555, 11174 kN/m, respectively. The experimental equivalent damping Ce are 3206, 1672, 865, 506 kN·s/m and the corresponding numerical results are 2010, 1623, 855, 294. As shown in Figure 10(a), there is little difference between the numerical and experimental results, and the maximum error is 1.1%. As shown in Figure 10(b), when loading frequency is between 0.2 Hz and 0.5 Hz, the numerical results agreewell with experimental values, and the maximum error is 2.9%; however, when the loading frequency is at 0.1 Hz or 1.0 Hz, there are some differences between experimental and numerical results.
Figure 11 shows the comparison results between experimental and numerical results of the MEIMD under different displacement amplitudes with a fixed frequency of 0.5 Hz. The experimental data for the equivalent stiffness are 12795, 11116, 9604 and 7027 kN/m under 2 mm, 4 mm, 6 mm and 8 mm displacement amplitudes, and the corresponding numerical data are 12632, 10986, 9555 and 8309 kN/m, respectively.
Comparison of equivalent stiffness and equivalent damping in calculation results and test results under 0.5 Hz frequency of the multi-dimensional earthquake isolation and mitigation device.
The experimental equivalent dampings are 1201, 1035, 865 and 625 kN…s/m, and the corresponding numerical data are 1215, 1035, 855 and 674 kN…s/m. It is shown from Figure 11(a) that the numerical results agree well with experiment results for the equivalent stiffness when the excitation displacement amplitude is between 2 mm and 6 mm, and the maximum error is 1.2%, while the maximum error is 18.26% under 8 mm displacement amplitude working condition. As shown in Figure 11(b), when the displacement amplitude is between 2 mm and 6 mm, the numerical results agree well with the experimental results, and the maximum error of the equivalent damping is 1.2%, and the maximum error is 7.8% under 8 mm working condition. It can be found from the analysis that the mathematical equations can describe the equivalent stiffness and damping of the MEIMD under different excitation frequency and displacement amplitudes.
5. Conclusions
In this paper, a new MEIMD has been developed and horizontal property tests are carried out. In accordance with the experimental data, the numerical simulation on the equivalent stiffness and the equivalent damping of the device is carried out. The following conclusions can be obtained from the experimental and numerical analysis.
The force-displacement hysteretic loops of the MEIMD are full smooth ellipse, which means that the device has a good performance of energy dissipation. The properties of the MEIMD are affected by the excitation frequency, the excitation amplitude, and the vertical pressure load. Excitation frequency has obvious influence on the equivalent stiffness and damping, excitation amplitude has a smaller influence, and vertical pressure load has a slight influence. The mathematical equations can describe the equivalent stiffness and damping of the MEIMD under different excitation frequency and displacement amplitudes.
Footnotes
Acknowledgements
The authors are grateful for the support of their funding bodies.
Funding
This work was supported by the National Natural Science Foundation of China (grant number 90915004), Jiangsu Province Natural Science Foundation (grant number BK2011282), Doctoral Fund of Ministry of Education of China (grant number 20090092110012) and the Priority Academic Program Development of Jiangsu Higher Education Institutions.
