Estimation of damping capacity of rubber vibration isolators under harmonic excitation

Svetlana Polukoshko1

1Ventspils University College, Engineering Research Institute “VSRC”, Ventspils, Latvia

1Corresponding author

Vibroengineering PROCEDIA, Vol. 8, 2016, p. 50-56.
Received 3 August 2016; accepted 6 August 2016; published 7 October 2016

Copyright © 2016 JVE International Ltd. 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.
Creative Commons License

The elastomeric compensation device mounted between the vibrating base and protecting object is the main element of any the passive vibration protection system. Rubber and rubber-like materials (elastomers) are the most suitable material for such devices due to their capability of absorbing input energy much better than other construction materials. In this paper rubber absorber of prismatic form with parallel flat ends under axial harmonic kinematic excitation is considered. For the estimation of damping properties of the rubber vibroabsorber two approaches are applied: using the mechanical models based on combination of elastic and viscous elements and using analytical representation of integral equations of creep and relaxation. In this article Maxwell and Burgers mechanical models and Rabotnov’s kernel for analytical representation are used. Damping properties are expressed by the ratio of amplitude of the driving vibration to the amplitude of the forced oscillations of object. Numerical calculations for each method are fulfilled and results are compared.

Keywords: elastomers, viscoelastic behaviour, relaxation kernel, dynamic modulus, absorption.

1. Introduction

Nowadays the questions of preventing the harmful effects of vibration and impact on machines and construction elements requires special attention. In most cases it is necessary to provide of shock and vibration protection at the design stage and best of all to use rubber-like materials for designing such devices [1-5].

The elastic properties of rubber are characterized by a large difference between bulk and shear moduli (K/G ratio reaches 500-5000), the high elasticity of rubbers allows them to withstand the large elongation without breaking; consequently, the rubber has a greater capacity to store energy per unit volume compared with other materials [3-7]. Along with the immediate elastic deformation rubber exhibit a retarded elastic deformation, viscous flow (creep) and relaxation, so its mechanical properties strongly depend on the time of application of external loads, that is, the mode of deformation [6-8]. If the duration of the force causing a deformation equal to the time of relaxation processes in the rubber, the elastic deformation is fully manifested, and this mode of deformation will correspond to well-defined modulus of rubber. When rate of deforming load application increases modulus of elasticity will be changing, and its magnitude will be larger than elastic modulus for the mode with full implementation of the elastic process, this fact must be taking into account when designing shock and vibration absorber. For the optimal design the rubber antivibration devices it is necessary to know mechanical characteristics of elastomers: bulk modulus of compression, dynamic and static shear modulus, damping properties, energy dissipation factor, relationship between stress and strain accounting viscoelastic properties.

In this paper the vibration absorber operating under uniaxial compression – tenion deformation is considered. In Fig. 1 the scheme of mounting and work of this device are shown: elastomeric absorber in the form of rectangular cylindrical or prism is placed between protected object and vibrating base. Rubber element 1 is rigidly attached to the undeformable plates 2. The aim of this work is to determine the dependence of damping properties of rubber absorber on the viscoelastic behaviour of elastomer, on the forced vibration frequency and natural frequency of the system for the model with one degree of freedom (ODF) under the action of harmonic excitation. Mass of absorber is not taken into account because of its smallness. The damping properties are determined by the ratio of the amplitude of the upper plate and mass forced oscillations to the amplitude of the forcing oscillation of base and lower plate. Such approach allows to appoint the dimentions of antivibration devices at design stage.

Fig. 1. Schemes of the installed absorbers: 1 – rubber element of height h, 2 – rigid plates, Δ – static deformation caused by Pz

Fig. 2. Models of viscoelastic medium: a) Maxwell’s model, b) Voigt model, c) Kelvin’s standard linear medium model, d) Burgers model

2. Simulation versions of viscoelastic properties of rubber

In the mechanical properties of rubber a set of features manifest itself: alongside with the instantaneous elastic deformation these materials exhibit retarding elastic deformation, viscous flow (creep) and relaxation. For clear presentation of the role of elastic and dissipative processes in the mechanical behaviour of materials the mechanical models are used which simplify the understanding of physical meaning of components of a differential equation, describing the relationship between stress and strain in time. It is assumed that for the most materials this equation is linear, and all features of mechanical behaviour may be described by system combining of two basic types of mechanical elements: elastic and viscous [8]. The elastic element is controlled by Hooke’s law σ=Eε, viscous element – by Newton’s law of flow σ=ηdε/dt. In the elastic element stresses changing is accompanied by an instantaneous directly proportional deformation changing. At constant stress the elastic element is in equilibrium state. The viscous element, in contrast, at constant stress is in a steady flow (at a constant strain rate). Two main elements may be connected together in series or in parallel joint forming two historically the first models of viscoelastic medium: Maxwell’s model (Fig. 2(a)) and Voigt model (Fig. 2(b)), which approximately describe the visco-elastic behaviour of the material. Stress – strain relationship for unilateral deformation are: for Voigt model: σ=Eε+ηε˙, for Maxwell’s model:

ε ˙ = σ ˙ E + σ η ,         σ ˙ = E ε ˙ - σ τ ,

where σ – normal stress, ε – relative deformation. E – Young’s modulus, ηviscosity factor, τ=η/Erelaxation time, which is necessary to establish a new equilibrium state of the material.

Both models have the certain disadvantages, using only them it is impossible to describe all aspects of the mechanical behavior of rubber-like materials. Thus, the Maxwell’s model allows instantaneous elastic deformation and stress relaxation, but does not simulate the equilibrium state of the material under load. Voigt’s model takes into account the retarded elastic deformation, but does not allow deformation jumps and stress relaxation. For adequate describing all features of the mechanical behavior of rubber-like materials it is necessary to complicate the model by increasing the number of elements used. There are four possible schemes of three-elements models and eight schemes of four-elements models. No one of these models is universal, since each of them has its limits of application. The example of three-elements Kelvin’s model, or standard linear medium model, is given in Fig. 2(c), example of four-elements Burgers model – in Fig. 2(d). The relationships σtεt for Kelvin’s and Burgers model under unilateral deformation:

σ ˙ + σ E 1 + E 2 η = ε E 1 E 2 η + ε ˙ E 2 ,           1 E 2 σ ¨ + 1 η 1 + 1 η 2 + E 1 E 2 η 1 σ ˙ + E 1 η 1 η 2 σ = ε ¨ + E 1 η 1 ε ˙ .

The blocks of elements may be, in turn, joined in series. Increasing number of elements (more than four) does not lead to a qualitative refinement time picture of stress – strain.

The most common relationship between stress and strain may be expressed by means of integral equations of relaxation and creep. In this case accounting of the time factor influence is reduced to solving systems of integro-differential equations. Methods of solving such systems are very diverse, the most widely used methods based on the use of integral transforms, especially Laplace and Fourier transforms, and also integral operators of Volterra’s type [5, 8]. In the case of axial deformation for weakly filled rubbers and small deformation (ε 20 %), when the linear dependence between stress and strain exists, the dependence is valid: σt=Etεt, where Et– time operator, not depended on σ and ε and reflects the viscoelastic behaviour of system.

This dependence may be written in the form of an integral equation of Volterra:

σ t = E 0 ε t - - t R t - τ ε τ d τ ,

where: E0 – instantaneous compressive modulus, which specifies the instantaneous stiffness in the classical sense; Rt-τ – relaxation kernel describing the viscoelastic properties of the rubber. Using Fourier – images transform of σω, εω and Rω Eq. (3) may be written as:

σ ω = E 0 1 - R ω ε ω = E ω ε ω ,           E ω = E ' ω + i E ' ' ω ,
E ' ω = E 0 R e 1 - R ω ,           E ' ' ω = - E 0 I m R ω ,
R e 1 - R ω = B 1 ω ,           - I m R ω = B 2 ω ,           E ω = E 0 B 1 ω + i E 0 B 2 ω ,

where Eω – dynamic compression modulus of elasticity expressed as complex modulis.

Proposed by A. Rzhanitsin, Y. Rabotnov and M. Koltunov kernels are in most common use. These kernels allow satisfactory describing mechanical properties of rubber in medium frequency range 10-3 <ω< 103 s-1. In this work Rabotnov’s kernel is considered, which represents fractional exponential function Rtτ=Эαβ,tτ. Rabotnov’s kernel and its Fourier image:

R t - τ = 1 - τ α n = 0 - β n t - τ n 1 - α Γ n + 1 1 + α ,           R ω = 1 i ω 1 + α + β ,

where 1 <α< 0, β> 0, α – parameter of fractionality, expressed in terms of the maximum value of the spectrum of relaxation times and modulus defect, β – rheological characteristics of rubber, Γ – gamma-function. Parameters α and β are in detail tabulated for many types rubbers.

Taking into account Eqs. (4) real and imaginary parts we obtain from Eq. (5):

B 1 ω = ω 1 + α c o s δ + β ω 2 1 + α + 2 ω 1 + α β c o s δ + β 2 ,           B 2 ω = ω 1 + α s i n δ ω 2 1 + α + 2 ω 1 + α β c o s δ + β 2 ,
δ = π 2 1 + α .

The coefficients B1 and B2 are Fourier cosine and sinus transform of the function Rtτ. The elastic modulus dependence on frequency should be considered as the conditional characteristics, as they are not applicable to any of a periodic process with a known frequency, but only to monochromatic described by harmonic functions sin(ωt), cos(ωt), exp(iωt).

3. Models of damping system with one degree of freedom

The dynamic properties of rubber absorber it is possible to estimate most clearly in the frame of the simplified one-degree-of-freedom model (Fig. 1), which allows to obtain explicit analytical dependences for damping characteristics. Protected object is represented as a point mass, vibration isolator is placed between protected object and vibrating base with kinematic excitation in accordance with time dependence of ζt. The time dependence of the protected object motion is zt. Since ζ and z are defined in the inertial system, the relative motion of the object and base will be equal to zt-ζt, and deformation of the elastomer layers thickness h is equal to:

ε t = z t - ζ t / h .

Let function ζt is periodic with amplitude ζ0:ζ=ζ0sin(ωt). Since the system is linear, for the harmonic input function all caused output processes will also be harmonic functions with the same frequency but with phase delay. Then it is natural to assume that function zt will be as follow: z=z0cos(ωt+φ), where z0 – amplitude of steady-state vibrations of mass, φ – angle of lag. The force imposed on the object with mass m from the elastic support is equal F=Aσt, where A – area of the upper face of the antivibration device.

The equation of motion of mass m in z – direction is:

m z ¨ t = - A σ t .

For Maxwell’s model the equation of mass m motion is received accounting Eqs. (1, 7, 8) and denoting natural frequency of mass vibration ω0=AE/mh:

m z ¨ t = - A ε ˙ η - σ ˙ η E = - A ε ˙ η + m η A E z t ,         z t + E η z ¨ t + ω 0 2 z ˙ t = ω 0 2 ζ ˙ t .

For harmonic deformation taking into account linear independency of the functions sin(ωt) and cos(ωt), the system of the equations is obtained from Eq. (9):

c o s φ - s i n φ s i n φ c o s φ 1 - b 2 E η b 2 ω z 0 = ζ 0 1 0 ,           b = ω ω 0 .

Received from system Eq. (10) forced and forcing vibration amplitudes ratio z0/ζ0 and lag angle φ:

z 0 ζ 0 ω = η ω 0 2 E 2 ω 2 + η 2 ω 2 - ω 0 2 2 ,         φ ω = a r c t a n E ω η ω 2 - ω   0 2 .

For Burgers model taking into account Eqs. (2, 7, 8) and denoting ω0=AE2/mh (since E2 is accounted initial modulus) equation of motion, amplitudes ratio z0/ζ0and lag angle φ are received:

z I V t + E 1 + E 2 η 1 + E 2 η 2 z t + E 1 E 2 η 1 η 2 + ω 0 2 z ¨ t + E 1 η 1 ω 0 2 z ˙ t = ω 0 2 ζ ¨ t + E 1 η 1 ω 0 2 ζ ˙ t ,
z 0 ζ 0 ω = ω 0 2 E 1 2 + η 1 2 ω 2 ω η 1 ω 2 - ω 0 2 - E 1 E 2 η 1 η 2 2 + E 1 ω 0 2 η 1 ω - ω E 2 η 2 + E 1 + E 2 η 1 2 ,
φ ω = a r c t a n E 1 ω 2 - ω 0 2 - E 1 E 2 η 1 η 2 + E 1 ω 0 2 η 1 ω - ω E 2 η 2 + E 1 + E 2 η 1 ω η 1 E 1 E 1 ω 0 2 η 1 ω - ω E 2 η 2 + E 1 + E 2 η 1 - ω 2 - ω 0 2 - E 1 E 2 η 1 η 2 ω η 1 .

Mass motion equation depending on the viscoelastic properties of the rubber is received taking into account Eqs. (3, 8):

m z ¨ t = - A E 0 ε t - - t R t - τ ε τ d τ .

Denoting the natural frequency of vibrating mass as ω0 and taking into account Eq. (7), the equation of object motion is received in the form:

z ¨ t + A E 0 m h z t - ζ t - - t R t - τ z τ - ζ τ d τ = 0 ,             ω 0 = A E 0 h m ,
z ¨ ( t ) + ω 0 2 z 0 s i n ω t + φ - ζ 0 s i n ω t
          - - t R t - τ z 0 sin ω τ + φ - ζ 0 sin ω τ d τ = 0 .

Using the property Rt-τ= 0 if t<τ, for harmonic deformation the system of the equations is obtained from Eq. (15):

c o s φ - s i n φ s i n φ c o s φ 1 - B 1 ω - b 2 B 2 ω z 0 = ζ 0 1 - B 1 ω B 2 ω ,         b = ω ω 0 .

From system Eq. (16) the amplitudes ratio or of the forced vibration amplitude of mass z0ω and the angle of lag φω may be found:

z 0 ζ 0 ω = ω 0 2 B 1 ω - 1 2 + B 2 2 ω B 2 2 ω ω 4 + B 1 ω - 1 2 ω 0 2 + B 1 ω - 1 ω 2 + B 2 2 ω ω 0 2 2 ,
φ ω = a r c t a n - B 2 ω ω 2 B 1 ω - 1 2 ω 0 2 + B 1 ω - 1 ω 2 + B 2 2 ω ω 0 2 .

By the same way parameters of shear vibration may be determined.

4. Numerical examples

The results of calculation of damping properties of rubber vibroisolator in accordance with above described algorithms and equations are given below. Numerical examples was solved with help of Mathcad-15 program. The curves of ratio z0/ζ0 dependence functions on ω, as well as φ on ω functions are constructed according to Eqs. (11-13, 17, 18); in all cases solution was fulfilled for five natural frequencies ω0= 1, 5, 10, 15 and 20 s-1.

Solution based on Rabotnov’s kernel is presented in Fig. 3, where plots of dependence the amplitudes ratio (a) and angle of lag (b) on excitation frequency are given for low-filled rubber 2959. Rheological parameters of this rubber α= –0.6 and factor β= 1.062 [6].

In Fig. 4(a) plots of dependence of the amplitudes ratio z0/ζ0 on excitation frequency ω calculated in accordance with Maxwell model, in Fig. 4(b) in accordance with Burgers model are presented.

In both cases calculations were performed for rubber 8157. Maxwell model parameters E= 2.55 MPa, η= 1.04 MNs/m2 and Burgers model parameters E1= 8.93 MPa, E2= 3.58 MPa and η1= 0.057 MNs/m2, η2= 1.04 MNs/m2 were taken from [5].

It is seen from the plots that simple Maxwell model describes the behaviour of rubber absorber inaccurately, but it may be used in a first approximation.

Fig. 3. Plots of dependence of a) amplitude ratio z0/ζ0 on excitation frequency ω and b) the lag angle φ on excitation frequency ω, calculated based on Rabotnov’s kernel for ω0= 1, 5, 10,15 and 20 s-1



Fig. 4. Plots of dependence of a) amplitude ratio z0/ζ0 on excitation frequency ω calculated with accordance to Maxwell mode, b) calculated in accordance to Burgers model for ω0= 1, 5, 10,15 and 20 s-1



5. Conclusions

This paper presents the procedure of estimation of damping properties of rubber absorber for ODF system under harmonic excitation. Definition of damping properties is fulfilled for system with kinematic excitation based on two approaches: using Maxwell and Burgers mechanical models combining of elastic and viscous elements and using Rabotnov’s kernel of relaxation for analytical representation of viscoelastic properties of rubber. Damping properties are expressed by the ratio of the amplitude of the driving lower base to the amplitude of the forced oscillations of the upper base and the object. Rubber absorbers of rectangular prismatic form, compressed in the normal to flat surface direction, is calculated as the examples.

Symbolic and numerical computations was fulfilled on Mathcad-15 program. Calculations confirm the well-known fact from the theory of vibration protection: if natural frequency of the object is three times lower than frequency of the forcing oscillations, the forced oscillation amplitude is attenuated approximately in ten times. Such vibration isolators show good results on above resonance frequencies. It may be concluded that proposed approach is suitable for dimensioning of rubber absorber on design stage.

Form factor accounting, i.e. shape influence on object natural frequency and experimental study of damping properties of rubber vibroisolator will be the next step of current research work.


  1. Kelly J. M., Konstantinidas D. A. Mechanics of Rubber Bearings for Seismic and Vibration Isolation. John Wiley and Sons, UK, 2011. [CrossRef]
  2. GentA. N. Engineering with Rubber: How to Design Rubber Components. Carl Hanser Verlag, Munich, 2011. [CrossRef]
  3. Bauman J. T. Fatigue, Stress and Strain of Rubber Components: Guide for Design Engineers. Carl Hanser Verlag, Munich, 2008. [CrossRef]
  4. LepetovV. A. Rubber Technical Products. Chemistry, Moscow, 1972, (in Russian). [CrossRef]
  5. Lyapunov V. T., Lavendel E. E., Shlyapochnikov S. A. Rubber Vibration Isolators. Handbook, Sudostrojenie, Leningrad, 1988, p. 216, (in Russian). [CrossRef]
  6. Ormonbekov T. O., Begaliev U. T. Use of Thin Rubber-Metal Elements in Earthquake Protection of Buildings, Engineering Constructions and Equipment. Ilim, Bishkek, Kyrgyzstan, 1996. [CrossRef]
  7. Poturaev V. N., Dyrda V. I. Rubber Machine Elements. Mashinostrojenie, Moscow, 1977. [CrossRef]
  8. Iljushin A. A., Pobedrja B. E. Fundamentals of Mathematical Theory Thermo-Viscoelasticity. Nauka, Moscow, 1970. [CrossRef]