Research Article  Open Access
FractionalDerivative Maxwell Kelvin Model for “5+4” Viscoelastic Damping Wall Subjected to Large Deformation
Abstract
Considering the larger vibration amplitude and several viscoelastic material layers, a fractionalderivative 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 inparallel 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.
1. Introduction
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 [1]. For example, 260 viscoelastic dampers were installed in the Columbia Center building in Seattle in order to decrease windinduced 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 [2]. Similar implementation of viscoelastic dampers for seismic mitigation can be also found in [3]. 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 fractionalderivative 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 [17] studied sheartype buildings with Maxwell modelbased bracedamper systems. The generalized Maxwell model (GMM) was experimentally verified by simulating the behavior of a fluid damper [18]. Although these models are simple and can be easily applied in practical engineering, studies [19–21] have shown that the parameters obtained under lowfrequency vibration could be inaccurate when vibration frequency increased to a higher extent. A twoparameter cooling law was proposed to offer a dependable estimate of the internal temperature of fluid dampers [20]. 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 [23] presented an approximate analytical solution to calculate the force of viscous dampers subjected to lowfrequency motions. A family of parameters identification methods were proposed for the Maxwell and Kelvin fractional models [24]. The fractionalderivative 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 [26]. The classical and fractionalderivative 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 [6] discussed two specific examples by constitutive modeling of viscoelastic solid damper and viscoelastic liquids damper. Considering the larger vibration amplitude, the literature [29] 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 threedimensional 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 [29]. Figure 3 shows the forcedisplacement hysteresis loops measured in a sinusoidal test by PA30240 (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 forcedeformation loops exhibit high nonlinearity.
(a) Geometry [29] (unit: mm)
(b) Threedimensional 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 inparallel combination of the fractional Maxwell model and the fractional Kelvin model, as schematically shown in Figure 4.
The RiemannLiouville fractional derivative is defined as [30]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 mm^{2} 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 [29] 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 [4] discussed the method of parameters identification of fiveparameter fractionalderivative model, which utilized the asymptotic properties of the storage and loss modulus functions. Makris [31] mentioned one more detailed description of the least square method for identification procedure. Park [6] 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 Zshaped, 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 [32]. 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
4. Conclusions
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.
Competing Interests
The authors declare that they have no competing interests.
Acknowledgments
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.
References
 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
 S. S. Tezcan and O. Uluca, “Reduction of earthquake response of plane frame buildings by viscoelastic dampers,” Engineering Structures, vol. 25, no. 14, pp. 1755–1761, 2003. View at: Publisher Site  Google Scholar
 T. Pritz, “Fiveparameter fractional derivative model for polymeric damping materials,” Journal of Sound and Vibration, vol. 265, no. 5, pp. 935–952, 2003. View at: Publisher Site  Google Scholar
 R. M. Christensen, Theory of Viscoelasticity: An Introduction, Academic Press, New York, NY, USA, 1971.
 S. W. Park, “Analytical modeling of viscoelastic dampers for structural and vibration control,” International Journal of Solids and Structures, vol. 38, no. 4445, pp. 8065–8092, 2001. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 A. Aprile, J. A. Inaudi, and J. M. Kelly, “Evolutionary model of viscoelastic dampers for structural applications,” Journal of Engineering Mechanics, vol. 123, no. 6, pp. 551–560, 1997. View at: Publisher Site  Google Scholar
 A. Palmeri, F. Ricciardelli, A. De Luca, and G. Muscolino, “State space formulation for linear viscoelastic dynamic systems with memory,” Journal of Engineering Mechanics, vol. 129, no. 7, pp. 715–724, 2003. View at: Publisher Site  Google Scholar
 T. Hatada, T. Kobori, M. Ishida, and N. Niwa, “Dynamic analysis of structures with Maxwell model,” Earthquake Engineering & Structural Dynamics, vol. 29, no. 2, pp. 159–176, 2000. View at: Publisher Site  Google Scholar
 K.W. Min, J. Kim, and S.H. Lee, “Vibration tests of 5storey steel frame with viscoelastic dampers,” Engineering Structures, vol. 26, no. 6, pp. 831–839, 2004. View at: Publisher Site  Google Scholar
 K. C. Chang, T. T. Soong, S.T. Oh, and M. L. Lai, “Seismic behavior of steel frame with added viscoelastic dampers,” Journal of Structural Engineering, vol. 121, no. 10, pp. 1418–1426, 1995. View at: Publisher Site  Google Scholar
 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
 V. A. Matsagar and R. S. Jangid, “Viscoelastic damper connected to adjacent structures involving seismic isolation,” Journal of Civil Engineering and Management, vol. 11, no. 4, pp. 309–322, 2005. View at: Publisher Site  Google Scholar
 M. P. Singh and T.S. Chang, “Seismic analysis of structures with viscoelastic dampers,” Journal of Engineering Mechanics, vol. 135, no. 6, pp. 571–580, 2009. View at: Publisher Site  Google Scholar
 T.S. Chang and M. P. Singh, “Mechanical model parameters for viscoelastic dampers,” Journal of Engineering Mechanics, vol. 135, no. 6, pp. 581–584, 2009. View at: Publisher Site  Google Scholar
 Y.T. Chen and Y. H. Chai, “Effects of brace stiffness on performance of structures with supplemental Maxwell modelbased bracedamper systems,” Earthquake Engineering & Structural Dynamics, vol. 40, no. 1, pp. 75–92, 2011. View at: Publisher Site  Google Scholar
 L.Y. Lu, G.L. Lin, and M.H. Shih, “An experimental study on a generalized Maxwell model for nonlinear viscoelastic dampers used in seismic isolation,” Engineering Structures, vol. 34, pp. 111–123, 2012. View at: Publisher Site  Google Scholar
 M. Sjöberg and L. Kari, “Nonlinear behavior of a rubber isolator system using fractional derivatives,” Vehicle System Dynamics, vol. 37, no. 3, pp. 217–236, 2002. View at: Publisher Site  Google Scholar
 C. J. Black and N. Makris, “Viscous heating of fluid dampers under small and large amplitude motions: experimental studies and parametric modeling,” Journal of Engineering Mechanics, vol. 133, no. 5, pp. 566–577, 2007. View at: Publisher Site  Google Scholar
 R. W. Wolfe, H.B. Yun, S. Masri, F. Tasbihgoo, and G. Benzoni, “Fidelity of reducedorder models for largescale nonlinear orifice viscous dampers,” Structural Control and Health Monitoring, vol. 15, no. 8, pp. 1143–1163, 2008. View at: Publisher Site  Google Scholar
 S. Dzierzek, “Experimentbased modeling of cylindrical rubber bushings for the simulation of wheel suspension dynamic behavior,” SAE Technical Paper 2000010095, 2000. View at: Publisher Site  Google Scholar
 N. Makris, M. C. Constantinou, and G. F. Dargush, “Analytical model of viscoelastic fluid dampers,” Journal of Structural Engineering, vol. 119, no. 11, pp. 3310–3325, 1993. View at: Publisher Site  Google Scholar
 R. Lewandowski and B. Chorazyczewski, “Identification of the parameters of the KelvinVoigt and the Maxwell fractional models, used to modeling of viscoelastic dampers,” Computers and Structures, vol. 88, no. 12, pp. 1–17, 2010. View at: Publisher Site  Google Scholar
 R. Lewandowski and Z. Pawlak, “Dynamic analysis of frames with viscoelastic dampers modelled by rheological models with fractional derivatives,” Journal of Sound and Vibration, vol. 330, no. 5, pp. 923–936, 2011. View at: Publisher Site  Google Scholar
 N. Makris and M. C. Constantinou, “Fractionalderivative maxwell model for viscous dampers,” Journal of Structural Engineering, vol. 117, no. 9, pp. 2708–2724, 1991. View at: Publisher Site  Google Scholar
 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.
 J. H. Xu and A. Q. Li, “Experiments of mechanical properties of super large viscoelastic damping wall under shear loading along transverse direction,” Journal of Southeast University (Natural Science Edition), vol. 45, no. 1, pp. 133–138, 2015. View at: Publisher Site  Google Scholar
 I. Podlubny, Fractional Differential Equations, vol. 198 of Mathematics in Science and Engineering, Academic Press, 1999. View at: MathSciNet
 N. Makris, “Complexparameter Kelvin model for elastic foundations,” Earthquake Engineering & Structural Dynamics, vol. 23, no. 3, pp. 251–264, 1994. View at: Publisher Site  Google Scholar
 The Japan Society of Seismic Isolation, Passhibu Seishin Kouzou Sekkei. Sekou Manyuaru Sono, The Japan Society of Seismic Isolation, Tokyo, Japan, 2008.
Copyright
Copyright © 2016 Junhong Xu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.