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.
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 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 [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 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 [17] 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 [18]. 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 [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 low-frequency motions. A family of parameters identification methods were proposed for the Maxwell and Kelvin fractional models [24]. 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 [26]. 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 [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 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 [29]. 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.
“5+4” super large viscoelastic damping wall.
Geometry [29] (unit: mm)
Three-dimensional view of loading device
The connection between the upper steel plate and the flange
The connection between the bottom steel plate and the test bed
Test power source.
High pressure nitrogen gas bottle
Piston energy storage system
F-u hysteresis curves of VDW at different frequencies and displacement values.
0.5 Hz
1.0 Hz
1.5 Hz
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.
A schematic show of FDMK viscoelastic mechanical model.
The Riemann-Liouville fractional derivative is defined as [30](1)Drft=1Γ1-rddt∫0tft-ττrdτ=ddt∫0tIrτft-τdτ=Irtf0+∫0tIrτddtft-τdτ,where r is the fractional order of derivative with respect to time t and 0<r<1. Γ is the gamma function. And Γ(z)=∫0∞tz-1e-tdt, z>0.
The constitutive equation of the Abel dashpot is defined as(2)σt=ηDrεt=ηIrt∗dεt.
So, the fractional Kelvin model in Figure 4 can be written as(3)σ1t=η1Drε1t+G1ε1t,where η1 is the viscosity coefficient of the Abel dashpots and G1 is the elastic coefficient in the fractional Kelvin model.
The fractional Maxwell model in Figure 4 can be written as(4)σ2t+η2G2Drσ2t=η2Drε2t,where η2 is the viscosity coefficient of the Abel dashpots and G2 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 (5)σ1t=η1Drε1t+G1ε1t,σ2t+η2G2Drσ2t=η2Drε2t,εt=ε1t=ε2t,σt=σ1t+σ2t.By rearranging (5), the relationship between stress and strain can be explicitly formulated as(6a)σt=Θεt,where(6b)Θ=G2η1η2DrG1+η1DrG1η2+η1η2DrG2η2+η22Dr+G2η22Dr.
The Fourier transform of (6a) and (6b) is written as (7)σ~ω=G2η1η2iωrG1+η1iωrε~ωG1η2+η1η2iωrG2η2+η22iωr+G2η22iωr.The symbol i denotes the imaginary unit; ω denotes the circular frequency, where ω=2πf, and f 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 (G′, the storage modulus of FDMK model, and G′′, the loss modulus of FDMK model) can be obtained:(8)G′=G1G2η22ωrsinrπ/2+G2η2η12ωrsinrπcosrπ,G′′=G1G2η2η1ωrcosrπ/2+G2η1η22ωrcosrπ/2+ω2rcosrπ/2sinrπ.
So, the loss factor η can be expressed as(9)η=G′′G′=G1G2η2η1ωrcosrπ/2+G2η1η22ωrcosrπ/2+ω2rcosrπ/2cosrπG1G2η22ωrsinrπ/2+G2η2η12ωrsinrπsinrπ.
The viscoelastic damper can be characterized by the storage stiffness Kd′ and loss stiffness Kd′′ which are related to G′ and G′′ as(10)Kd′=G′Ah=G1G2η22ωrsinrπ/2+G2η2η12ωrsinrπAhcosrπ,Kd′′=G′′Ah=G1G2η2η1ωrcosrπ/2+G2η1η22ωrcosrπ/2+ω2rcosrπ/2Ahsinrπ,where A is the total shear area and h is the thickness of the viscoelastic damping material slab. The values of A and h are 1000000 mm2 and 10 mm, respectively. The equipment damping coefficient Cd is defined by(11)Cd=Kd′′ω=G1G2η2η1ωrcosrπ/2+G2η1η22ωrcosrπ/2+ω2rcosrπ/2Ahωsinrπ.
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 (η1, G1, η2, G2, and r) 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.
Properties of “5+4” viscoelastic damping wall (at 225% strain and 21.2°C).
Frequency (Hz)
G′ (MPa)
G′′ (MPa)
η (MPa·s)
Kd′ (kN/mm)
Kd′′ (kN/mm)
Cd (kN·s/mm)
f=0.5
2.65
0.76
0.29
2.65
0.76
1.52
f=1.0
1.79
1.05
0.59
1.79
1.05
0.70
f=1.5
1.34
0.85
0.63
1.34
0.85
0.57
f=2.0
1.16
0.89
0.77
1.16
0.89
0.45
f=5.0
1.00
0.91
0.91
1.00
0.91
0.40
f=10.0
0.96
0.94
0.98
0.96
0.94
0.38
Properties of “5+4” viscoelastic damping wall (at 300% strain and 21.2°C).
Frequency (Hz)
G′ (MPa)
G′′ (MPa)
η (MPa·s)
Kd′ (kN/mm)
Kd′′ (kN/mm)
Cd (kN·s/mm)
f=0.5
2.36
0.67
0.29
2.36
0.67
1.34
f=1.0
2.69
0.63
0.23
2.69
0.63
0.42
f=1.5
1.89
0.63
0.33
1.89
0.63
0.42
f=2.0
2.43
0.62
0.26
2.43
0.62
0.31
f=5.0
2.20
0.61
0.28
2.20
0.61
0.28
f=10.0
2.10
0.50
0.24
2.10
0.50
0.25
Properties of “5+4” viscoelastic damping wall (at 450% strain and 21.2°C).
Frequency (Hz)
G′ (MPa)
G′′ (MPa)
η (MPa·s)
Kd′ (kN/mm)
Kd′′ (kN/mm)
Cd (kN·s/mm)
f=0.5
3.54
0.29
0.08
3.54
0.29
0.58
f=1.0
3.20
0.29
0.09
3.20
0.29
0.19
f=1.5
2.89
0.32
0.11
2.89
0.32
0.21
f=2.0
2.72
0.33
0.12
2.72
0.33
0.17
f=5.0
2.71
0.34
0.13
2.71
0.34
0.16
f=10.0
2.65
0.36
0.14
2.65
0.36
0.11
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 five-parameter fractional-derivative 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 G∗ (G∗=G′+iG′′) 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:(12)minI~G1,G2,η1,η2,r=∑i=16ri-si2.
Note that i=1–6, 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 ri and si denote the data point obtained from the experiment and the FDMK model for the ith excitation frequency, respectively. The results of the parameters identification are shown in Table 4, together with their normalized error.
Best-fitting parameters of the FDMK model.
G1 (MPa)
G2 (MPa)
η1 (MPa·s)
η2 (MPa·s)
r
Value
8.69
3.58
0.75
0.98
0.60
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 [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.
Comparison between the experimental and numerically reconstructed hysteresis loops.
At 225% strain, frequency 0.5 Hz
At 225% strain, frequency 1.5 Hz
At 300% strain, frequency 0.5 Hz
At 300% strain, frequency 1.5 Hz
At 450% strain, frequency 0.5 Hz
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 (13)Δ=AExp-ATheAExp×100%,where AExp and AThe 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.
Comparison of the areas of the experimental and theoretical hysteresis loops.
Frequency (Hz)
Theoretical value (kN·mm)
Experimental value (kN·mm)
Δ
0.5 (225%)
2279
2413
5.6%
1.5 (225%)
2104
2249
6.9%
0.5 (300%)
2849
2826
0.8%
1.5 (300%)
2445
2540
3.7%
0.5 (450%)
5240
5394
2.9%
1.5 (450%)
4427
4530
2.3%
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 G′, the loss modulus G′′, the storage stiffness Kd′, the loss stiffness Kd′′, the equivalent damping coefficient Cd, and the loss factor η. These values from experiment are presented in Tables 1–3. The FDMK viscoelastic mechanical model representations for G′ and G′′ are given in the form of (8). Additionally, Kd′ and Kd′′ are related to G′ and G′′ 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 Kd′ and the loss stiffness Kd′′ are linear functions of the storage modulus G′ and the loss modulus G′′, respectively. The curve shapes for G′ and Kd′ are similar to each other; likewise, G′′ and Kd′′ 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 G′′ and Kd′′ calculated by FDMK model are closer to the experimental test than the values of G′ and Kd′. This indicates that the FDMK model is more sensitive in energy dissipation than in energy storage. As shown in Figure 7(a), the G′′ value under ω=1.0 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 ω=1.0. These general trends can provide a beneficial guideline for engineers to predict the behavior of the damping wall subjected to different excitation.
Model simulation of the storage modulus.
At 225% strain
At 450% strain
Model simulation of the loss modulus.
At 225% strain
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 Cd 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.
Model simulation of the equivalent damping coefficient.
At 225% strain
At 450% strain
Model simulation of the loss factor.
At 225% strain
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.
ChristopoulosC.FiliatraultA.2006Pavia, ItalyIUSS PressZhaoG.PanP.QianJ.LinJ.Experimental study of viscoelastic dampers subjected to large deformation201233101261332-s2.0-84869130544TezcanS. S.UlucaO.Reduction of earthquake response of plane frame buildings by viscoelastic dampers200325141755176110.1016/j.engstruct.2003.07.0012-s2.0-0242406957PritzT.Five-parameter fractional derivative model for polymeric damping materials2003265593595210.1016/S0022-460X(02)01530-42-s2.0-0042884059ChristensenR. M.1971New York, NY, USAAcademic PressParkS. W.Analytical modeling of viscoelastic dampers for structural and vibration control20013844-458065809210.1016/S0020-7683(01)00026-9ZBL1090.746362-s2.0-0035850692AprileA.InaudiJ. A.KellyJ. M.Evolutionary model of viscoelastic dampers for structural applications1997123655156010.1061/(asce)0733-9399(1997)123:6(551)2-s2.0-0000067881PalmeriA.RicciardelliF.De LucaA.MuscolinoG.State space formulation for linear viscoelastic dynamic systems with memory2003129771572410.1061/(ASCE)0733-9399(2003)129:7(715)2-s2.0-0042313959HatadaT.KoboriT.IshidaM.NiwaN.Dynamic analysis of structures with Maxwell model200029215917610.1002/(sici)1096-9845(200002)29:2<159::aid-eqe895>3.0.co;2-12-s2.0-0034143698MinK.-W.KimJ.LeeS.-H.Vibration tests of 5-storey steel frame with viscoelastic dampers200426683183910.1016/j.engstruct.2004.02.0042-s2.0-16544395354ChangK. C.SoongT. T.OhS.-T.LaiM. L.Seismic behavior of steel frame with added viscoelastic dampers1995121101418142610.1061/(asce)0733-9445(1995)121:10(1418)2-s2.0-0029391851FerryJ. D.1980New York, NY, USAJohn Wiley & SonsTschoeglN. W.1989Berlin, GermanySpringer10.1007/978-3-642-73602-5MatsagarV. A.JangidR. S.Viscoelastic damper connected to adjacent structures involving seismic isolation200511430932210.1080/13923730.2005.96363622-s2.0-28844450374SinghM. P.ChangT.-S.Seismic analysis of structures with viscoelastic dampers2009135657158010.1061/(asce)0733-9399(2009)135:6(571)2-s2.0-66349095669ChangT.-S.SinghM. P.Mechanical model parameters for viscoelastic dampers2009135658158410.1061/(asce)0733-9399(2009)135:6(581)2-s2.0-66349089024ChenY.-T.ChaiY. H.Effects of brace stiffness on performance of structures with supplemental Maxwell model-based brace-damper systems2011401759210.1002/eqe.10232-s2.0-78650562064LuL.-Y.LinG.-L.ShihM.-H.An experimental study on a generalized Maxwell model for nonlinear viscoelastic dampers used in seismic isolation20123411112310.1016/j.engstruct.2011.09.0122-s2.0-80155154094SjöbergM.KariL.Non-linear behavior of a rubber isolator system using fractional derivatives200237321723610.1076/vesd.37.3.217.35322-s2.0-0036494719BlackC. J.MakrisN.Viscous heating of fluid dampers under small and large amplitude motions: experimental studies and parametric modeling2007133556657710.1061/(asce)0733-9399(2007)133:5(566)2-s2.0-34247232653WolfeR. W.YunH.-B.MasriS.TasbihgooF.BenzoniG.Fidelity of reduced-order models for large-scale nonlinear orifice viscous dampers20081581143116310.1002/stc.2562-s2.0-56849122160DzierzekS.Experiment-based modeling of cylindrical rubber bushings for the simulation of wheel suspension dynamic behavior20002000-01-009510.4271/2000-01-0095MakrisN.ConstantinouM. C.DargushG. F.Analytical model of viscoelastic fluid dampers1993119113310332510.1061/(asce)0733-9445(1993)119:11(3310)2-s2.0-0027789968LewandowskiR.ChorazyczewskiB.Identification of the parameters of the Kelvin-Voigt and the Maxwell fractional models, used to modeling of viscoelastic dampers2010881-211710.1016/j.compstruc.2009.09.0012-s2.0-70350770826LewandowskiR.PawlakZ.Dynamic analysis of frames with viscoelastic dampers modelled by rheological models with fractional derivatives2011330592393610.1016/j.jsv.2010.09.0172-s2.0-78649448161MakrisN.ConstantinouM. C.Fractional-derivative maxwell model for viscous dampers199111792708272410.1061/(asce)0733-9445(1991)117:9(2708)2-s2.0-0026219915LaiM.LuP.LunsfordD.Viscoelastic damper: a damper with linear or nonlinear materialProceedings of the 11th World Conference on Earthquake Engineering1996Mexico City, MexicoElsevier ScienceXuZ. D.2001Xi'an, ChinaXi'an University of Architecture and TechnologyXuJ. H.LiA. Q.Experiments of mechanical properties of super large viscoelastic damping wall under shear loading along transverse direction201545113313810.3969/j.issn.1001-0505.2015.01.0242-s2.0-84923763515PodlubnyI.1999198Academic PressMathematics in Science and EngineeringMR1658022MakrisN.Complex-parameter Kelvin model for elastic foundations199423325126410.1002/eqe.42902303032-s2.0-0028184349The Japan Society of Seismic Isolation2008Tokyo, JapanThe Japan Society of Seismic Isolation