Fractional-Derivative Maxwell Kelvin Model for “5+4” Viscoelastic Damping Wall Subjected to Large Deformation
Considering the larger vibration amplitude and several viscoelastic material layers, a fractional-derivative Maxwell Kelvin (FDMK) viscoelastic mechanical model is proposed for “5+4” viscoelastic damping wall, which is used for vibration control of building structures. The development of the model is based on in-parallel combination of fractional Maxwell model and fractional Kelvin model. The proposed model is experimentally validated and very good agreement between predicted and experimental results was obtained. The results confirm that the FDMK model is accurate in simulating the hysteresis properties of the “5+4” viscoelastic damping wall under large deformation. From the areas of the experimental and theoretical hysteresis loops, under 300% strain, the predicted result is the most accurate in prediction of the energy dissipation and the second is the prediction under 450% strain. Moreover, from the comparisons of dynamic properties (storage modulus, loss modulus, etc.), the FDMK model works satisfactorily. The FDMK model is more sensitive in energy dissipation than in energy storage.
Protection of the mechanical response of the building structures subjected to earthquake has become an increasingly critical issue. Viscoelastic dampers have long been used in controlling the vibration in civil engineering . For example, 260 viscoelastic dampers were installed in the Columbia Center building in Seattle in order to decrease wind-induced vibration. A pioneer application can be traced back to New York in 1969 where around 10 000 viscoelastic dampers were installed in the towers of the World Trade Center to resist wind loads . Similar implementation of viscoelastic dampers for seismic mitigation can be also found in . However, general application of the structures with added viscoelastic dampers cannot be significantly developed without comprehensive understanding of the constitutive behavior of viscoelastic dampers. Current researches on constitutive modeling mainly focus on the classical and fractional-derivative models which can be found in [4–11].
Classical constitutive modeling is the diverse combination of springs and dashpots [12–16], for example, the Maxwell (M) and Kelvin (K) models, Burger model, generalized Maxwell model (GMM), and other models [17, 18]. The literature  studied shear-type buildings with Maxwell model-based brace-damper systems. The generalized Maxwell model (GMM) was experimentally verified by simulating the behavior of a fluid damper . Although these models are simple and can be easily applied in practical engineering, studies [19–21] have shown that the parameters obtained under low-frequency vibration could be inaccurate when vibration frequency increased to a higher extent. A two-parameter cooling law was proposed to offer a dependable estimate of the internal temperature of fluid dampers . Under this kind of situation, more complex model parameters should be included or otherwise the original model would lose its accuracy [6, 22]. However, this can make the classical model too cumbersome for general application. To solve the problem, fractional models, such as the fractional Kelvin model and fractional Maxwell model, were proposed and received increasingly considerable attention [23–25]. The literature  presented an approximate analytical solution to calculate the force of viscous dampers subjected to low-frequency motions. A family of parameters identification methods were proposed for the Maxwell and Kelvin fractional models . The fractional-derivative models can accurately capture the viscoelastic damping properties of the dampers by using even fewer model parameters.
The dynamic characteristics of a viscoelastic damper largely depend on the geometry of the device once the material is determined . The classical and fractional-derivative models [6, 27, 28] mentioned above were mainly verified based on the “3+2” (composed of two viscoelastic layers bonded with three parallel steel plates) viscoelastic dampers. The literature  discussed two specific examples by constitutive modeling of viscoelastic solid damper and viscoelastic liquids damper. Considering the larger vibration amplitude, the literature  proposed a “5+4” (composed of four viscoelastic layers bonded with five parallel steel plates) super large viscoelastic damping wall (VDW) under large deformation (at 225%, 300%, and 450% strain). Compared with the traditional viscoelastic damping device, “5+4” VDW has a substantial number of advantages. For example, “5+4” VDW can supply large space for doors or windows in a building. The placement of “5+4” VDW can be also equal to the interstory displacement so as to improve the efficiency of energy dissipation with a great deal of application potential.
However, the constitutive model of “5+4” VDW is not yet sufficiently validated by test results. Due to the large deformation and four viscoelastic material layers, highly nonlinear viscoelastic behavior may exist in such damper. The objective of this paper is to propose a reasonable mathematical model for “5+4” VDW. As will be shown, the proposed model, called FDMK, will be derived by parallel combination of the Maxwell and Kelvin components. It will be validated by characterizing various damping properties (storage modulus, loss modulus, storage stiffness, loss stiffness, equivalent damping coefficient, and loss factor) of the “5+4” VDW.
The paper is organized as follows. In Section 2, constitutive model of “5+4” VDW is derived. In Section 3, parameters of the derived model are identified with nonlinear least square method by using MATLAB (version 7.0). Then, the model behavior is investigated in the frequency domain and verified by fitting the experimental data.
2. Dynamic Characteristics and Modeling of “5+4” VDW
Experimental studies have demonstrated that the damping properties of viscoelastic dampers are related not only to the viscoelastic material, but also to the damper deformation. The structure diagram of the “5+4” VDW can be found in Figure 1(a) and the three-dimensional view of the corresponding loading device can be found in Figure 1(b). The upper steel plate is connected to the flange structures (green part in Figure 1(b)). The details are shown in Figure 1(c). The bottom steel plate is fixed at the test bed (purple part in Figure 1(b)). The details are shown in Figure 1(d). Experimental details, such as the test material, test setup, and test results, are available in . Figure 3 shows the force-displacement hysteresis loops measured in a sinusoidal test by PA30-240 (Tobul Accumulator Inc., as shown in Figure 2). It is quite evident that increasing the deformation generally results in larger cycles, that is, higher temporal shifting between stress and strain. It can be also observed that the force-deformation loops exhibit high nonlinearity.
(a) Geometry  (unit: mm)
(b) Three-dimensional view of loading device
(c) The connection between the upper steel plate and the flange
(d) The connection between the bottom steel plate and the test bed
(a) High pressure nitrogen gas bottle
(b) Piston energy storage system
(a) 0.5 Hz
(b) 1.0 Hz
(c) 1.5 Hz
(d) 2.0 Hz
The FDMK model proposed by this paper is considered as in-parallel combination of the fractional Maxwell model and the fractional Kelvin model, as schematically shown in Figure 4.
The Riemann-Liouville fractional derivative is defined as where is the fractional order of derivative with respect to time and . is the gamma function. And , .
The constitutive equation of the Abel dashpot is defined as
So, the fractional Kelvin model in Figure 4 can be written aswhere is the viscosity coefficient of the Abel dashpots and is the elastic coefficient in the fractional Kelvin model.
The fractional Maxwell model in Figure 4 can be written aswhere is the viscosity coefficient of the Abel dashpots and is the elastic coefficient in the fractional Maxwell model.
Due to parallel combination, (3) and (4) can be treated as a set of equations which can be written in the form of By rearranging (5), the relationship between stress and strain can be explicitly formulated aswhere
The Fourier transform of (6a) and (6b) is written as The symbol denotes the imaginary unit; denotes the circular frequency, where , and is the frequency. and are the Fourier transforms of the stress and strain histories, respectively.
By performing inverse Fourier transform of (7), the following relationships (, the storage modulus of FDMK model, and , the loss modulus of FDMK model) can be obtained:
So, the loss factor can be expressed as
The viscoelastic damper can be characterized by the storage stiffness and loss stiffness which are related to and aswhere is the total shear area and is the thickness of the viscoelastic damping material slab. The values of and are 1000000 mm2 and 10 mm, respectively. The equipment damping coefficient is defined by
3. Model Validation
In this study, the FDMK model characterized by (8), (9), (10), and (11) will be employed to simulate the viscoelastic damping behavior of “5+4” VDW under large deformation. As shown above, there are only five characteristic parameters (, , , , and ) to be determined. The test results of “5+4” VDW by Xu and Li  as shown in Tables 1, 2, and 3 are used to validate the proposed model.
An important problem of the FDMK model is the identification of the model parameters from experimental data. Pritz  discussed the method of parameters identification of five-parameter fractional-derivative model, which utilized the asymptotic properties of the storage and loss modulus functions. Makris  mentioned one more detailed description of the least square method for identification procedure. Park  determined the parameters of two mechanical models of viscoelastic dampers. In this paper, the parameters are obtained by using nonlinear least square method within MATLAB to analyze a set of storage modulus and loss modulus data (defined as the square root of the sum of squares of storage modulus and loss modulus data). In the least square method, the error between the model and the experimental complex modulus () is minimized in order to search the model parameters. The optimal parametric values of the FDMK model are obtained by the minimum of the following objective function (see (12)) which is defined by the error functional minimum between the theoretical and experimental data:
Note that , since there are a total of 6 excitation frequencies in the test, that is, 0.5, 1.0, 1.5, 2.0, 5.0, and 10.0 Hz. Symbols and denote the data point obtained from the experiment and the FDMK model for the th excitation frequency, respectively. The results of the parameters identification are shown in Table 4, together with their normalized error.
3.1. Comparisons of Hysteresis Loops between Test and Simulation
Figure 5 compares the experimental hysteresis loops with the simulation by the FDMK model for the excitation frequencies of 0.5 and 1.5 Hz under 225% strain, 300% strain, and 450% strain. The red solid lines in Figure 5 represent the experimental loops, while the dotted lines depict the theoretical loops simulated by using the parameters listed in Table 4. Notably, the experimental and theoretical loops have almost the same relative displacement and force. Figure 5 also shows that the hysteresis loops are Z-shaped, whereas both loops do not show an inclined elliptical shape and behaved in a nonlinear viscoelastic manner because the viscoelastic materials belonged to Si content . Figure 5(a) demonstrates that the simulated results obtained by the FDMK model underestimate the damper force in the third quadrant at 225% strain and frequency of 0.5 Hz but overestimate the damper force in the fourth quadrant at 225% strain and frequency of 1.5 Hz. Figures 5(c), 5(d), 5(e), and 5(f) illustrate that the theoretical results are in good agreement with the experimental loops under large deformation.
(a) At 225% strain, frequency 0.5 Hz
(b) At 225% strain, frequency 1.5 Hz
(c) At 300% strain, frequency 0.5 Hz
(d) At 300% strain, frequency 1.5 Hz
(e) At 450% strain, frequency 0.5 Hz
(f) At 450% strain, frequency 1.5 Hz
Moreover, the area of a hysteresis loop is equivalent to the ability of the energy dissipated by the damper. Table 5 compares the area of the theoretical loops to the experimental loops under frequencies 0.5 and 1.5 Hz at three strains, 225%, 300%, and 450%, respectively. Note that Δ in Table 5 denotes the error percentage of the dissipated energy simulated by the FDMK model. The error percentage is defined as where and denote the areas of the experimental and theoretical hysteresis loops, respectively. The areas of the experimental and theoretical hysteresis loops were integrated by MATLAB programs and the values are shown in Table 5. Based on the error percentage calculated by (13), under 300% strain, the predicted result is the most accurate in prediction of the energy dissipation and the second is the prediction under 450% strain. From the area values, at 225% strain and the frequency of 0.5 Hz, the model underestimates the experimental results, while, at 225% strain and the frequency of 1.5 Hz, the model overestimates the experimental results. This is consistent with the results presented in Figure 5.
3.2. Comparisons of Mechanical Properties
The mechanical properties of viscoelastic damping wall are important technical parameters used for the designing of vibration reduction in engineering structures. So, this section will demonstrate the model prediction of the mechanical properties, typically the storage modulus , the loss modulus , the storage stiffness , the loss stiffness , the equivalent damping coefficient , and the loss factor . These values from experiment are presented in Tables 1–3. The FDMK viscoelastic mechanical model representations for and are given in the form of (8). Additionally, and are related to and as (10).
The results obtained at 225% and 450% strains will be briefly described. A comparison of the storage modulus resulting from the FDMK model with the corresponding experimental results is presented in Figure 6. A similar comparison for the loss modulus is shown in Figure 7. Note that, in view of (10), the storage stiffness and the loss stiffness are linear functions of the storage modulus and the loss modulus , respectively. The curve shapes for and are similar to each other; likewise, and share the same trends. The presented results demonstrate that the FDMK model could reasonably describe the modulus of the damping wall under large deformation. The presented model works satisfactorily. Moreover, the values of and calculated by FDMK model are closer to the experimental test than the values of and . This indicates that the FDMK model is more sensitive in energy dissipation than in energy storage. As shown in Figure 7(a), the value under exhibits mutation. This may be explained as follows: the mutation possibility of the energy dissipation of the viscoelastic material used for VDW has occurred under . These general trends can provide a beneficial guideline for engineers to predict the behavior of the damping wall subjected to different excitation.
(a) At 225% strain
(b) At 450% strain
(a) At 225% strain
(b) At 450% strain
The plot of the equivalent damping coefficient versus frequency is shown in Figure 8. It can be observed that the equivalent damping coefficient of the viscoelastic damping wall can be well characterized by the FDMK model. The results of (Figure 9) show some departures from the FDMK representation, but the overall trend of the curve coincides with the experimental results. Therefore, it is suggested that the design values could be increased in actual design of “5+4” viscoelastic damping wall under large deformation.
(a) At 225% strain
(b) At 450% strain
(a) At 225% strain
(b) At 450% strain
This paper proposed a FDMK model for “5+4” viscoelastic damping wall subjected to large deformation. The FDMK model was composed of a fractional Kelvin model and a fractional Maxwell model, which were connected in a parallel way. Comparison between the experimental and predicted results showed that the theoretical results matched very well the experimental hysteresis loops under the large deformation. From the areas of the experimental and theoretical hysteresis loops, under 300% strain, the predicted result was the most accurate in prediction of the energy dissipation and the second was the prediction under 450% strain.
Moreover, the FDMK model was more sensitive in energy dissipation than in energy storage. The design values could be increased in actual design of “5+4” viscoelastic damping wall under large deformation.
The authors declare that they have no competing interests.
The authors wish to acknowledge the National Science Foundation of China (51278104) and Key Projects of Application Development Plan in Chongqing (cstc2014yykfB30003) in connection with this paper.
C. Christopoulos and A. Filiatrault, Principles of Passive Supplemental Damping and Seismic Isolation, IUSS Press, Pavia, Italy, 2006.
G. Zhao, P. Pan, J. Qian, and J. Lin, “Experimental study of viscoelastic dampers subjected to large deformation,” Journal of Building Structures, vol. 33, no. 10, pp. 126–133, 2012.View at: Google Scholar
R. M. Christensen, Theory of Viscoelasticity: An Introduction, Academic Press, New York, NY, USA, 1971.
J. D. Ferry, Viscoelastic Properties of Polymers, John Wiley & Sons, New York, NY, USA, 1980.
N. W. Tschoegl, The Phenomenological Theory of Linear Viscoelastic Behavior, Springer, Berlin, Germany, 1989.View at: Publisher Site
M. Lai, P. Lu, and D. Lunsford, “Viscoelastic damper: a damper with linear or nonlinear material,” in Proceedings of the 11th World Conference on Earthquake Engineering, Elsevier Science, Mexico City, Mexico, 1996.View at: Google Scholar
Z. D. Xu, Experiment and Study about the (Lead) Viscoelastic Structure, Xi'an University of Architecture and Technology, Xi'an, China, 2001.
I. Podlubny, Fractional Differential Equations, vol. 198 of Mathematics in Science and Engineering, Academic Press, 1999.View at: MathSciNet
The Japan Society of Seismic Isolation, Passhibu Seishin Kouzou Sekkei. Sekou Manyuaru Sono, The Japan Society of Seismic Isolation, Tokyo, Japan, 2008.