Published: 30 June 2014

Design and analysis of a quasi-zero stiffness isolator using a slotted conical disk spring as negative stiffness structure

Fu Niu1
Lingshuai Meng2
Wenjuan Wu3
Jinggong Sun4
Wenming Zhang5
Guang Meng6
Zhushi Rao7
1, 5, 6, 7State Key Laboratory of Mechanical System and Vibration, School of Mechanical Engineering Shanghai Jiao Tong University, 800 Dongchuan Road, Shanghai, 200240, China
2, 3, 4, 1Institute of Medical Equipment, Academy of Military Medical Sciences, Tianjin, 300161, China
Corresponding Author:
Guang Meng
Views 445
Reads 239
Downloads 1958

Abstract

This paper concerns the characteristics of a novel quasi-zero stiffness (QZS) isolator developed by parallelly combining a slotted conical disk spring with a vertical linear spring. The static characteristics of the slotted conical disk spring as well as the QZS isolator are presented. The configurative parameters are optimized to achieve a wide displacement range around the equilibrium position for which the stiffness has a low value and changes slightly. The overload and underload conditions are taken into account, resulting in a Helmoholtz-Duffing equation. The primary resonance response of the nonlinear system composed by a loaded mass and the QZS isolator are determined by employing the Harmonic Balance Method (HBM) and confirmed with the results of numerical simulation. The frequency response curves (FRCs) are obtained for both force and displacement excitations. The force transmissibility, the absolute displacement and acceleration transmissibility are defined and investigated. The study shows that the overloaded or underloaded system can exhibit linear stiffness, softening stiffness, softening-hardening stiffness and hardening stiffness with the increasing excitation amplitude. The response and the resonance frequency of the system are affected by the excitation amplitude and the offset displacement to the position at which the dynamic stiffness is zero. To enlarge the isolation frequency range and improve the isolation performance, the loaded mass and the excitation amplitude should be suitably controlled.

1. Introduction

Low-frequency vibrations are more than often thought to induce harmful effects affecting the health of drivers and reducing the accuracy of high-precision machinery. In the ideal case of a mass supported by a linear stiffness on a rigid foundation, the vibration attenuation is only obtained for the excitation frequency is greater than 2 times the natural frequency of the passive linear isolator. It is evident that a smaller stiffness results in a border frequency band of isolation but leads to a larger static displacement of the supported mass as well [1]. In recent years, passive nonlinear vibration isolators, owning a high static stiffness resulting in a small static displacement and a low dynamic stiffness resulting in a wide frequency range of isolation, have drawn increasing attention since it can overcome the disadvantage. By choosing the appropriate system parameters of nonlinear isolators, a quasi-zero stiffness (QZS) isolator possessing zero dynamic stiffness at equilibrium position can be realized [2].

QZS isolators are mainly achieved by parallelly combining a negative stiffness structure with a positive stiffness structure. Ibrahim [3] reviewed the recent advances in passive nonlinear isolators with good ultra-low frequency isolation performance in detail. The comprehensive book by Alabuzev [4] et al. covered the fundamental theory and many prototypes of vibration protecting systems characterized with QZS. Platus [5] combined positive and negative stiffness with two compressed bars hinged at the center. Zhang et al. [6] added a beam under axial force to a positive stiffness spring. Carrella et al. [2] proposed a High-Static-Low-Dynamic stiffness (HSLDS) isolator with a vertical linear spring in parallel with two oblique linear springs. Le and Ahn [7] studied a vibration isolator for vehicle seat exposed to low frequency excitation theoretically and experimentally, in which the negative stiffness structure is configured by a horizontal spring in series with a bar. Liu et al. [8] investigated the characteristics of a QZS isolator using Euler buckled beam as negative stiffness corrector.

In this paper, a new QZS isolator including a vertical linear spring and a slotted conical disk spring is presented as shown in Fig. 1. The reason for taking a slotted conical disk spring as the negative structure is that it can supply a certain restoring force at the flatten state and produce axial nonlinear restoring force, which enables the isolator to support greater load and achieve the QZS property at the equilibrium position.

Most studies of QZS isolators assume that the system keep balance at the equilibrium where the dynamic stiffness equals zero. However, the overload or underload condition is more common in practical engineering, which motivates this paper to some extent. The aim of this paper is to investigate the characteristics of the QZS isolator and the influence of overload or underload on the isolation performance of the QZS isolator.

Fig. 1Prototype model of the proposed QZS isolator

Prototype model of the proposed QZS isolator

2. Static characteristics of the QZS isolator

2.1. The slotted conical disk spring

Consider a slotted conical disk spring loaded axially as shown in Fig. 2. It can be divided into a coned disk segment and a number of lever arm segments. The solution for the straight slotted disk spring is fully stated by Schremmer based on the theory of Almen and Laszlo [9]. Note that the total displacement δ includes the rigid displacement δ1 and bending displacement δ2. The formulations between the axial applied force Fn and the displacement are given by:

1a
Fn=E1-μ2t3D2K1δ11+ht-δ1tht-δ12t1-DmD/1-dD,
1b
δ2=C(Dm-d)3(1-μ2)Fn2Et3b2Z,
1c
δ=1-dD/1-DmDδ1+δ2,

where E is the Young’s modulus of the slotted conical disk spring, μ is the Poisson’s ratio and the other parameters are presented in Fig. 2. While the constants K1 and C are defined as:

2a
K1=23πDDm2lnDDmDDm-12,
2b
C=3(1-b1/b2)312-2b1b2+b1b2232-lnb1b2.

The width ratio of slot tip b1 to slot datum b2 which is used to express the slot profile of a slotted conical disk spring does not significantly affect the nonlinear force-displacement [10]. Thus, simplified formula defining the slot profile with respect to the constant slot number Z, is introduced as follows:

3
b1=πd2Z, b2=πDm2Z.

Based on Eqs. (1) and (3), one can get the relationship between the total displacement δ and the rigid displacement δ1:

4
δ=1-dD/1-DmDδ1+CDm-d3πDmK1D21-DmD/1-dDδ11+ht-δ1tht-δ12t.

When the coned disk segment of the slotted conical disk spring is in a horizontal line which means δ1=h, the total displacement and the restoring force is (δe, Fne). It is worthy of note that the force-displacement curve is symmetric about the point as δ1 changes in the range [0, 2h], i.e. δ changes in the range [0, 2δe]. It is convenient to define the following non-dimensional parameters:

5a
δ^=δD,
5b
δ^1=δ1D,
5c
D^m=DmD,
5d
d^=dD,
5e
h^=hD,
5f
t^=tD,
5g
F^n=FnED2/(1-μ2),
5h
Γ=K1t^1-D^m1-d^,
5i
Μ=1-d^1-D^m,
5j
Ν=K1Cπ(D^m-d^)3D^mt^21-D^m1-d^.

Let x^=δ^1-h^, i.e. y^=δ^-δ^e. By substituting the two equations and Eq. (5) given above into Eqs. (1a) and (4), the non-dimensional restoring force and the non-dimensional total displacement can be derived:

6a
F^n=Γ12x^3+t^2-h^22x^+t^2h^,
6b
y^=Μx^+Ν12x^3+t^2-h^22x^,

where F^n is the non-dimensional force, x^ is the non-dimensional displacement of the coned disk segment, and y^ is the non-dimensional displacement of the slotted conical disk spring. In addition, δ^e can be derived as δ^e=(Μ+Νt^2)h^.

By differentiating Eq. (6a, b) with respect to the non-dimensional displacement x^ separately, the non-dimensional stiffness of the slotted conical disk spring is obtained as:

7
k^n=Γ32x^2+t^2-12h^2Μ+Ν32x^2+t^2-12h^2.

Fig. 2The slotted conical disk spring under axial force

The slotted conical disk spring under axial force

Considering the parameters Γ>0, Μ>0 and Ν>0, it can be found obviously that the stiffness is symmetric about x^=0 and gets the minimum value at this position as x^ changes in the range [-h^,h^], i.e. y^ changes in the range [-δ^e,δ^e]. Meanwhile, the slotted conical disk spring possesses continuous negative stiffness region when the parameters meet the condition h^(2t^,2(t^2+Μ/Ν)).

2.2. Design and optimization of the QZS isolator

The schematic of the QZS isolator with a vertical linear spring and a viscous damper in parallel with a slotted conical disk spring acting as the negative stiffness structure is shown in Fig. 3. Here, the weight of the isolated mass is ignored and the stiffness of the linear spring is kv. As displayed in Fig. 3(a), the mass moves a certain total displacement δ from the initial position by the force F. The restoring force can be expressed by:

8
F=kvδ+Fn=kv1-dD/1-DmDδ1+CDm-d3πDmK1D21-DmD/1-dDδ11+ht-δ1tht-δ12t+E1-μ2t3D2K1δ11+ht-δ1tht-δ12t1-DmD/1-dD.

By introducing the non-dimensional restoring force F^=F/kvD, and transforming the initial point to the symmetric point of the isolator following the procedure in Section 2.1. Then the non-dimensional restoring force can be obtained as:

9
F^=Μx^+h^+Ν+λΓ12x^3+t^2-h^22x^+t^2h^,

where λ=ED/kv(1-μ2) is defined as the stiffness ratio between the slotted conical disk spring and the linear spring, and the other parameters have the same meaning with that in Eq. (3). Differentiating Eqs. (9) and (6b) with respect to the non-dimensional displacement x^ separately, one can get the non-dimensional stiffness of the isolator:

10
k^=1+λΓ32x^2+t^2-h^22Μ+Ν32x^2+t^2-h^22,

where k^=k/kv is the non-dimensional stiffness.

Fig. 3Schematic representation of the QZS isolator. a) System balance at the position at which the dynamic stiffness is zero; b) overloaded system balance at a lower position; c) underloaded system balance at a higher position

Schematic representation of the QZS isolator. a) System balance at the position at which the dynamic stiffness is zero; b) overloaded system balance at a lower position;  c) underloaded system balance at a higher position

a)

Schematic representation of the QZS isolator. a) System balance at the position at which the dynamic stiffness is zero; b) overloaded system balance at a lower position;  c) underloaded system balance at a higher position

b)

Schematic representation of the QZS isolator. a) System balance at the position at which the dynamic stiffness is zero; b) overloaded system balance at a lower position;  c) underloaded system balance at a higher position

c)

In operation, it is prospected that the isolator can reach the static equilibrium at the state the coned disk segment of the slotted conical disk spring is horizontal after loaded with a mass, in which x^=x^e=0, i.e. y^=y^e=0, and the stiffness of the isolator has a minimum value. Then setting the stiffness of the isolator to be zero at the equilibrium position provides that:

11
λQZS=ΜΓ(h^2/2-t^2)-ΝΓ.

In addition to the isolator possessing the QZS property, it is desirable for it to own a wide range of non-dimensional total displacement R^ from the equilibrium position for which the non-dimensional stiffness has a low value as well as changes slightly. By substituting y^=y^e±R^=±R^ and x^=x^e±r^=±r^ into Eqs. (6b) and (10) separately, these displacements are found to satisfy:

12a
R^=Μr^+Ν12r^3+t^2-h^22r^,
12b
r^=23(k^-1)Μ+1-k^Ν+λΓh^22-t^2(1-k^)Ν+λΓ.

In this paper, the 50CrVA is used to be the material of the slotted conical disk spring, which has the following material properties: Young’s modulus E=206 GPa, Poisson’s ratio μ=0.3 and the ultimate stress S=1275 MPa. The configurative parameters D^m, d^, t^ and h^ are chosen from the set D^m[0.6, 0.9], d^[0.1, 0.5], t^[0.01, 0.1] and h^(2t^,2(t^2+Μ/Ν)) according to the practical engineering conditions [11]. Note that the stiffness ratio λ calculated using Eq. (11) should be positive and the non-dimensional stiffness k^ should not be negative. The optimization criteria includes the condition that the largest displacement is achieved at which the stiffness of the isolator is equal to that of the linear spring, i.e. k^=1, and the requirement that the non-dimensional stiffness changes slightly with the tolerance of Δk^=0.0025 for Δy^=0.01 in the region around the equilibrium position. In addition, the maximum stress of the slotted conical disk spring with the configurative parameters given above occurs at the lower outer edge and should not be larger than the ultimate stress [9].

The results of the optimization are D^m=0.83, d^=0.1, t^=0.01and h^=0.022. The relationship between the non-dimensional stiffness k^ and the non-dimensional displacement y^ for the optimal parameters (Case 1) and other parameters satisfying the optimization criteria listed in Table 1 can be plotted in Fig. 4. The circles shown in Fig. 4 denote the maximum range R^ given by Eq. (12) calculated at k^=1. It can be seen obviously that the optimized isolator possesses a very small stiffness in the region around the equilibrium position and a smaller stiffness for a larger displacement from the static equilibrium position.

Fig. 4Non-dimensional stiffness characteristics of the isolator for different parameters: ‘red line’ Case 1; ‘blue line’ Case2; ‘brown line’ Case 3; ‘green line’ Case 4

Non-dimensional stiffness characteristics of the isolator for different parameters:  ‘red line’ Case 1; ‘blue line’ Case2; ‘brown line’ Case 3; ‘green line’ Case 4

Table 1The configurative parameters of the QZS isolator

Case
D^m
d^
t^
h^
1
0.83
0.1
0.01
0.022
2
0.83
0.1
0.02
0.031
3
0.7
0.18
0.01
0.031
4
0.61
0.25
0.01
0.036

3. Dynamic characteristics of the QZS isolator

3.1. Approximation of the restoring force

The restoring force of the isolator F^ with respect to the total displacement y^ could be derived by Eqs. (6) and (9), which is very complicated and can be expanded to the third order Taylor series at the zero stiffness position y^=y^e=0. The approximate restoring force is given by:

13
F^=η+χy^3,

where η=[Μ+(Ν+λΓ)t^2]h^ and χ=(λΓ+Ν)4/2(λΓΜ)3. The exact restoring force expressed by Eqs. (6) and (9) and the approximate restoring force expressed by Eq. (13) are plotted in Fig. 5 both for the optimal parameters obtained in Section 2.2. Note that the error of the approximation depends on the non-dimensional displacement and less than 5 percent within the region |y^|0.06.

Fig. 5Non-dimensional force-displacement of the QZS isolator: ‘red line’ exact expression; ‘blue line’ approximate expression

Non-dimensional force-displacement of the QZS isolator:  ‘red line’ exact expression; ‘blue line’ approximate expression

3.2. Dynamic modeling and solution

As shown in Fig. 3(a), the QZS isolator will keep balance at the equilibrium position y^=y^e=0 when loaded with an appropriate mass, at which the dynamic stiffness is zero. Considering the practical applications, it is more possible that the isolator will balance at y^=+y^0 for overload or y^=-y^0 for underload as shown in Fig. 3(b) and 3(c) separately, owning an offset displacement y^0. The static equilibrium equation of the QZS isolator can be derived as:

14
kvDη±χy^03=mg.

The non-dimensional force-displacement and stiffness characteristics of the isolator for overload and underload can be plotted in Fig. 6. It is evident that the equilibrium position for overload and underload will moves from the original position denoted by ‘o’ to a new position denoted by ‘*’ and ‘+’ separately, and the non-dimensional stiffness will not be zero at the new equilibrium position. The influence of overload and underload on the isolation performance of the isolator cannot be ignored. Go back to Fig. (3), the mass m can represent a machine excited by the harmonic force F=fecos(ωt), or an equipment isolated from the base excitation Z=zecos(ωt). By applying the Newton’s second law of motion, one can get the dynamic equations of the system separately for the force and displacement excitations given above:

15a
my¨+cy˙+kvDη+χy^±y^03-mg=fecosωt,
15b
mu¨+cu˙+kvDη+χu^±y^03-mg=mω2zecosωt,

where u=y-Z is the relative displacement between the mass and the base, u^=u/D is the non-dimensional relative displacement. Introducing the non-dimensional parameters as follows:

ωn=kvm, τ=ωnt, Ω=ωωn, ξ=c2mωn, f^e=femωn2D, z^e=zeD,

and combing Eq. (14), Eqs. (15a) and (15b) can be rewritten separately as below:

16a
y^''+2ξy^'+3χy^02y^±3χy^0y^2+χy^3=f^ecosΩτ,
16b
u^''+2ξu^'+3χy^02u^±3χy^0u^2+χu^3=Ω2z^ecosΩτ.

Eq. (16) can be expressed by a uniform equation for simplicity:

17
v^''+2ξv^'+α1v^±α2v^2+χv^3=γβcosΩτ,

where α1=3χy^02, α2=3χy^0, β is the amplitude of the harmonic excitation, and γ=1 for the force excitation; γ=Ω2 for the displacement excitation. By applying the transformation e^=v^±α2/(3χ)=v^±y^0 into Eq. (17) [12], one can get a Duffing equation under asymmetric excitation:

18
e^''+2ξe^'+χe^3=±α+γβcosΩτ,

where the constant term α=χy^03. Using the Harmonic Balance Method (HBM), the approximate solution corresponding to the steady-state response in the region of the primary resonance is assumed to be:

19
e^=A0+A1cosΩτ+φ.

Fig. 6Non-dimensional force-displacement and stiffness characteristics of the overloaded or underloaded system. ‘*’ Equilibrium position of the overloaded system shown in Fig. 3(b); ‘+’ equilibrium position of the underloaded system shown in Fig. 3(c)

Non-dimensional force-displacement and stiffness characteristics of the overloaded or  underloaded system. ‘*’ Equilibrium position of the overloaded system shown in Fig. 3(b);  ‘+’ equilibrium position of the underloaded system shown in Fig. 3(c)

Substituting Eq. (19) into Eq. (18), equating constant terms, and the coefficients of the terms containing sin(Ωτ) and cos(Ωτ) separately to zero, the steady-state response can be expressed by the constant term A0, the amplitude of the harmonic term A1 and the phase φ as:

20a
χA03+32χA0A12=±α,
20b
-Ω2A1+3χA02A1+34χA13=γβcosφ,
20c
-2ξΩA1=γβsinφ.

Combing Eqs. (20a), (20b) and (20c) gives the implicit equation for the amplitude of the constant term A0:

21
25χ3A09-20χ2Ω2A07-15χ2αA06+4χΩ2(Ω2+4ξ2)A05+16χαΩ2A04+3χ(2γ2β2-3α2)A03-4αΩ2(Ω2+4ξ2)A02+4α2Ω2A0-α3=0.

The implicit equations for the peak amplitude of the constant term for the force and displacement excitation, i.e. Aopf and A0pd, can be derived separately from Eq. (21):

22a
-20χ3ξ2A0pf9+8χ2ξ4A0pf7+36χ2αξ2A0pf6-3χ2β2A0pf5-16χαξ4A0pf4-12χα2ξ2A0pf3+3χαβ2A0pf2+8α2ξ4A0pf-4α3ξ2=0,
22b
(-75χ4β2-80χ3ξ2)A0pd9+32χ2ξ4A0pd7+(45χ3αβ2+144χ2αξ2)A0pd6-64χαξ4A0pd4+(27χ2α2β2-48χα2ξ2)A0pd3+32α2ξ4A0pd+(3χα3β2-16α3ξ2)=0.

Then the excitation frequency ratios that corresponding to the maximum responses can be expressed by:

23a
Ω0pf2=5χ2A0pf6-4χξ2A0pf4-4χαA0pf3+4αξ2A0pf-α22χA0pf4-2αA0pf,
23b
Ω0pd2=5χ2A0pd6-4χξ2A0pd4-4χαA0pd3+4αξ2A0pd-α22χA0pd4+3χβ2A0pd2-2αA0pd.

The amplitude A1, the peak amplitude A1pf and A1pd of the harmonic term for the force and displacement excitation can be calculated separately from Eq. (20a). Eqs. (21), (22) and (23) are valid for the overloaded system balance at the position y^=+y^0. The solution for the response of the underloaded system balance at another position y^=-y^0 can be obtained by transforming y^0 to -y^0 in Eq. (21). According to the investigation of Kovacic et al. [13], it was found that the system excited by the asymmetric force may have a maximum number of one, three of five steady-state responses and multiple jumps for different combinations of excitation amplitudes.

For the system balance at the equilibrium position y^=y^e=0 after loaded with an appropriate mass, the steady-state solutions of the system for both the force and displacement excitation can be derived by following the procedure above and setting y^0=0. The dynamic equation is:

24
v^''+2ξv^'+χv^3=γβcosΩτ.

The implicit amplitude frequency equation is:

25
916χ2A16-32χΩ2A14+Ω2(Ω2+4ξ2)A12-γ2β2=0.

For the two types of excitation, the peak amplitude of the response A1pf and A1pd, and the excitation frequency ratios corresponding to the maximum responses can be obtained as follows:

26a
A1pf=2ξ3+4ξ6+3χβ23χξ,
26b
A1pd=8ξ248χξ2-9χ2β2,
27a
Ω1pf=4ξ6+3χβ2-6ξ34ξ,
27b
Ω1pd=92χβ2ξ6-12ξ812ξ6-454χβ2ξ4+278χ2β4ξ2-81256χ3β6.

The single degree of freedom (SDOF) system without the slotted conical disk spring can be regarded as the equivalent linear system (ELS). The dynamic equation for the ELS is:

28
v^''+2ξv^'+v^=γβcosΩτ.

The amplitude of the steady-state response for Eq. (28) can be derived as:

29
A1=γβ(1-Ω2)2+(2ξΩ)2 .

The optimal parameters D^m=0.83, d^=0.1, t^=0.01, h^=0.022 and the damping ratio ξ=0.03 are chosen to conduct the following investigation. Based on Descartes’ Rule of Signs [14] and the analysis above, the ways in which the maximum number of the steady-state values of A0 and A1 depend on the offset displacement y^0 and the excitation amplitude β are shown in Fig. 7.

Fig. 7The maximum number of the steady-state amplitudes: one, three or five, as a function of the non-dimensional offset displacement y^0 and non-dimensional excitation amplitude β

The maximum number of the steady-state amplitudes: one, three or five, as a function of the non-dimensional offset displacement y^0 and non-dimensional excitation amplitude β

As shown in Fig. 7, there are a maximum number of one, three and five steady-state values for both force excitation and displacement excitation. For the force excitation, the FRCs with a maximum number of one, three and five steady-state values corresponding to y^0=1.7×10-2, β=10-4; y^0=1.7×10-2, β=10-4; y^0=2.7×10-2, β=1.2×10-3 respectively, and for the displacement excitation, the FRCs with a maximum number of one, three and five steady-state values corresponding to y^0=1.7×10-2,β=10-3;y^0=1.7×10-2,β=4.8×10-3; y^0=2.7×10-2, respectively are shown in Fig. 8 to illustrate these cases. In Fig. 8, the approximate solutions obtained by HBM are compared with the exact solutions obtained by numerical simulation. By using the MATLAB ode45 function and calculating the first harmonic from the Fourier series coefficients of the steady-state response, the exact solutions are obtained as shown in Fig. 8 with the symbol ‘*’. It can be obviously seen that both the constant term and the first harmonic term of the response are predicted reasonably well in the frequency range using the HBM.

Fig. 8FRCs of the constant term A0 and harmonic term A1 for the optimal parameters and ξ=0.03. a) For the force excitation: ‘black line’ y^0=1.7×10-2, β=1×10-4; ‘red line’ y^0=1.7×10-2, β=8×10-4; ‘blue line’ y^0=2.7×10-2, β=1.2×10-3; b) For the displacement excitation: ‘black line’ y^0=1.7×10-2, β=1×10-3; ‘red line’ y^0=1.7×10-2, β=4.8×10-3; ‘blue line’ y^0=2.7×10-2, β=4.8×10-3; ‘green dashed line’ unstable solution, ‘o’ peak amplitude and ‘*’ numerical solution

FRCs of the constant term A0 and harmonic term A1 for the optimal parameters and ξ=0.03.  a) For the force excitation: ‘black line’ y^0=1.7×10-2, β=1×10-4; ‘red line’ y^0=1.7×10-2, β=8×10-4; ‘blue line’ y^0=2.7×10-2, β=1.2×10-3; b) For the displacement excitation: ‘black line’ y^0=1.7×10-2, β=1×10-3; ‘red line’ y^0=1.7×10-2, β=4.8×10-3; ‘blue line’ y^0=2.7×10-2, β=4.8×10-3; ‘green dashed line’ unstable solution, ‘o’ peak amplitude and ‘*’ numerical solution

4. Effects of the offset displacement and excitation amplitude on the QZS isolator

To investigate the influence of the offset displacement y^0 and excitation amplitude β on the FRCs and the transmissibility, the increasing values of the excitation amplitude β are chosen to be 10-4, 10-4, 10-4, 1.2×10-3 and 10-3 when the value of the offset displacement y^0 is fixed to be 1.7×10-2 or 2.7×10-2 separately for the system excited by the harmonic force. And then the increasing values of the excitation amplitude β are chosen to be 10-3, 10-3, 4.8×10-3, 10-3 and 10-2 when the value of the offset displacement y^0 is also fixed to be 1.7×10-2 or 2.7×10-2 separately for the system excited by the harmonic displacement. For the illumination convenience in the following sections, the system not being or being balance at the position at which the dynamic stiffness of the system is zero, i. e. the offset displacement y^0>0 or y^0=0, are called system I and system II separately.

4.1. Effects on the FRCs for the force excitation and displacement excitation

As shown in Fig. 9 and Fig. 10, the FRCs of system I and II are plotted separately for the both types of excitation, and that of the ELS is also plotted on the same figure for comparison. Note that the same values of the excitation amplitude β and the damping ratio ξ are chosen in the ELS.

Fig. 9FRCs of system I, system II and their ELS for force excitations with different offset displacements and excitation amplitudes. ‘red line’system I with y^0=1.7×10-2,‘blue line’ system I with y^0=2.7×10-2, ‘black line’ system II, ‘magenta line’ ELS, ‘green dotted line’ unstable solutions, ‘o’ peak amplitude of response

FRCs of system I, system II and their ELS for force excitations with different offset  displacements and excitation amplitudes. ‘red line’system I with y^0=1.7×10-2,‘blue line’  system I with y^0=2.7×10-2, ‘black line’ system II, ‘magenta line’ ELS, ‘green dotted  line’ unstable solutions, ‘o’ peak amplitude of response

a) β=10-4

FRCs of system I, system II and their ELS for force excitations with different offset  displacements and excitation amplitudes. ‘red line’system I with y^0=1.7×10-2,‘blue line’  system I with y^0=2.7×10-2, ‘black line’ system II, ‘magenta line’ ELS, ‘green dotted  line’ unstable solutions, ‘o’ peak amplitude of response

b) β=10-4

FRCs of system I, system II and their ELS for force excitations with different offset  displacements and excitation amplitudes. ‘red line’system I with y^0=1.7×10-2,‘blue line’  system I with y^0=2.7×10-2, ‘black line’ system II, ‘magenta line’ ELS, ‘green dotted  line’ unstable solutions, ‘o’ peak amplitude of response

c) β=10-4

FRCs of system I, system II and their ELS for force excitations with different offset  displacements and excitation amplitudes. ‘red line’system I with y^0=1.7×10-2,‘blue line’  system I with y^0=2.7×10-2, ‘black line’ system II, ‘magenta line’ ELS, ‘green dotted  line’ unstable solutions, ‘o’ peak amplitude of response

d) β=1.2×10-3

FRCs of system I, system II and their ELS for force excitations with different offset  displacements and excitation amplitudes. ‘red line’system I with y^0=1.7×10-2,‘blue line’  system I with y^0=2.7×10-2, ‘black line’ system II, ‘magenta line’ ELS, ‘green dotted  line’ unstable solutions, ‘o’ peak amplitude of response

e) β=10-3

By observing Fig. 9, for system I, a decrease in the offset displacement y^0 yields a decrease in both the amplitude of the constant term A0 and the resonance frequency, and an increase in the peak amplitude of the harmonic term A1 when the excitation amplitude is fixed. Note that the constant term A0 disappears in system II. With the increase of the excitation amplitude, the FRCs of the ELS is shifted upwards in the whole frequency region and the peak amplitudes always occur at Ω=1. For system I, the amplitude of the harmonic term A1 is also increases. But the effect of the excitation amplitude on the constant term A0 is only obvious around the resonance frequency, and larger excitation amplitude results in smaller peak amplitude of the constant term A0. For the frequency far away from the resonance frequency, the amplitude of the constant term A0 does not change much and will approach to the value of the offset displacement y^0. It is also worthy of note that the resonance frequency of system I decreases at first, increases later, and finally becomes larger than that of the ELS with the increasing excitation amplitude. For system II, larger excitation amplitude yields both larger amplitude of the response and the resonance frequency.

The excitation amplitude also effects the stiffness characteristics of system I and II obviously. Both system I and II approach to be a linear system if the excitation amplitude is small. When the excitation amplitude increases, system II always exhibits the hardening stiffness. However, with the increasing excitation amplitude, system I exhibits the softening stiffness, and then enters a region with the softening stiffness at first and the hardening stiffness later on. If the excitation amplitude keeps increasing, system I will only exhibit the hardening stiffness.

Fig. 10FRCs of system I, system II and their ELS for displacement excitations with different offset displacements and excitation amplitudes. ‘red line’ system I with y^0=1.7×10-2, ‘blue line’ system I with y^0=2.7×10-2, ‘black line’ system II, ‘magenta line’ ELS, ‘green dotted line’ unstable solutions, ‘o’ peak amplitude of response

FRCs of system I, system II and their ELS for displacement excitations with different  offset displacements and excitation amplitudes. ‘red line’ system I with y^0=1.7×10-2, ‘blue line’  system I with y^0=2.7×10-2, ‘black line’ system II, ‘magenta line’ ELS, ‘green dotted line’  unstable solutions, ‘o’ peak amplitude of response

a) β=10-3

FRCs of system I, system II and their ELS for displacement excitations with different  offset displacements and excitation amplitudes. ‘red line’ system I with y^0=1.7×10-2, ‘blue line’  system I with y^0=2.7×10-2, ‘black line’ system II, ‘magenta line’ ELS, ‘green dotted line’  unstable solutions, ‘o’ peak amplitude of response

b) β=10-3

FRCs of system I, system II and their ELS for displacement excitations with different  offset displacements and excitation amplitudes. ‘red line’ system I with y^0=1.7×10-2, ‘blue line’  system I with y^0=2.7×10-2, ‘black line’ system II, ‘magenta line’ ELS, ‘green dotted line’  unstable solutions, ‘o’ peak amplitude of response

c) β=4.8×10-3

FRCs of system I, system II and their ELS for displacement excitations with different  offset displacements and excitation amplitudes. ‘red line’ system I with y^0=1.7×10-2, ‘blue line’  system I with y^0=2.7×10-2, ‘black line’ system II, ‘magenta line’ ELS, ‘green dotted line’  unstable solutions, ‘o’ peak amplitude of response

d) β=10-3

FRCs of system I, system II and their ELS for displacement excitations with different  offset displacements and excitation amplitudes. ‘red line’ system I with y^0=1.7×10-2, ‘blue line’  system I with y^0=2.7×10-2, ‘black line’ system II, ‘magenta line’ ELS, ‘green dotted line’  unstable solutions, ‘o’ peak amplitude of response

e) β=10-2

In addition, another conclusion about the three systems can be drawn. For the same excitation amplitude, the amplitude of the harmonic term A1 in system II is larger than that in system I and the amplitude of the harmonic term A1 in system I is larger than that in ELS at lower frequencies; but when the excitation frequency is around the resonance frequency of the ELS, the amplitude of the harmonic term A1 in system II becomes smaller than that in system I and the amplitude of the harmonic term A1 in system I becomes smaller than that in ELS; and at higher frequencies the amplitudes of the harmonic term A1 in the three systems are almost the same.

Fig. 10 is for the displacement excitation, of which the vertical coordinate represents the relative displacement. Unlike the FRCs for the force excitation, the amplitudes of the harmonic term A1 in system I, II and the ELS will approach to the excitation amplitude with the increasing excitation frequency. The amplitude of the harmonic term A1 in system II reaches the excitation amplitude at the lower frequency than that in system I, and the amplitude of the harmonic term A1 in system I reaches the excitation amplitude at the lower frequency than that in ELS. When the excitation amplitude is relatively large, unbounded responses will occur for both system I and II which is not observed in Fig. 9.

4.2. Definition of the transmissibilities

Force transmissibility for the force excitation and absolute displacement transmissibility for the displacement excitation are the key indexes to evaluate the vibration isolation performance of a vibration isolator [15]. According to the study of Ravindra and Mallik [12], it is concluded that the nonlinear system with asymmetric restoring force may not perform satisfactorily for the displacement excitation by using the absolute displacement transmissibility to evaluate the isolation performance. Therefore, the absolute acceleration transmissibility for the displacement excitation is also introduced to evaluate the isolation performance for the displacement excitation.

4.2.1. Force transmissibility

The force transmissibility is defined as the ratio of the amplitude of the non-dimensional dynamic force transmitted to the base, to the amplitude of the non-dimensional excitation force. It is given by:

30
Tf=F^tβ,

where F^t=F^te2+F^td2, F^te is the non-dimensional elastic force and F^td is the non-dimensional damping force.

For system I, the solution of Eq. (17) can be expressed by [16]:

31
v^=A0'+A1cosΩτ+φ,

where the constant term A0'=A0y^0.

The non-dimensional elastic force is:

32
F^te=α1v^±α2v^2+χv^3.

By substituting Eq. (31) into Eq. (32), one can get the implicit expression of the non-dimensional elastic force:

33
F^te=F^t0+F^t1cosΩτ+φ,

where:

F^t0=α1A0'+α2A0'2+α2A122+χA0'3+32χA0'A12,
F^t1=α1A1+2α2A0'A1+34χA13+3χA0'2A1.

Only considering the dynamic force here, then the force transmissibility of system I can be expressed by:

34
TfI=F^t12+(2ξΩA1)2β.

For system II, the force transmissibility is determined by [2]:

35
TfII=(χA13)2+(2ξΩA1)2β.

Using Eq. (29), the force transmissibility of ELS can be derived as:

36
Tfl=1+(2ξΩ)2(1-Ω2)2+(2ξΩ)2.

4.2.2. Absolute displacement and acceleration transmissibility

The absolute displacement transmissibility is defined as the ratio between the amplitude of the non-dimensional absolute displacement and amplitude of the non-dimensional excitation displacement. It is given by:

37
Td=v0β.

For system I, one can get the expression of the absolute displacement of the mass by employing Eq. (18):

38
v^=e^y^0+z^=A0y^0+A1cosΩτ+φ+βcosΩτ.

The absolute displacement transmissibility can be defined as:

39
TdI=|A0y^0|+A12+β2+2A1βcos(φ)β.

Then the absolute acceleration transmissibility of system II can be derived as:

40
TacI=|v¨||z¨|=|v^''||z^''|=|u^''+z^''||z^''|=|e^''+z^''||z^''|=A12+β2+2A1βcos(φ)β,

where cos(φ) in Eqs. (39) and (40) can be obtained by Eq. (20b).

For system II, the absolute displacement and acceleration transmissibility have the same expression as:

41
TdII=TacII=A12+β2+2A1βcos(φ)β,

noting that the cos(φ) in Eq. (41) should be obtained by setting A0=0 in Eq. (20b).

For the ELS, the absolute displacement and acceleration transmissibility are also the same with the force transmissibility and given by:

42
Tdl=Tacl=1+(2ξΩ)2(1-Ω2)2+(2ξΩ)2.

4.3. Effects on the tranmissibilities

The force transmissibility is plotted in Fig. 11. The absolute displacement and acceleration transmissibility are plotted in Fig. 12. All the transmissibility results are plotted in dB, i.e. as 20log10Tf, 20log10Td and 20log10Tac. As shown in Fig. 11 and Fig. 12, both system I and II will exhibit superior or inferior isolation performance to the ELS depending on the excitation frequency range and the excitation amplitude.

Fig. 11Force transmissibility of system I, system II and their ELS for displacement excitations with different offset displacements and excitation amplitudes. ‘red line’ system I with y^0=1.7×10-2, ‘blue line’ system I with y^0=2.7×10-2, ‘black line’ system II, ‘magenta line’ ELS, ‘green dotted line’ unstable solutions, ‘o’ peak amplitude of transmissibility

Force transmissibility of system I, system II and their ELS for displacement excitations with different offset displacements and excitation amplitudes. ‘red line’ system I with y^0=1.7×10-2,  ‘blue line’ system I with y^0=2.7×10-2, ‘black line’ system II, ‘magenta line’ ELS,  ‘green dotted line’ unstable solutions, ‘o’ peak amplitude of transmissibility

a) β=10-4

Force transmissibility of system I, system II and their ELS for displacement excitations with different offset displacements and excitation amplitudes. ‘red line’ system I with y^0=1.7×10-2,  ‘blue line’ system I with y^0=2.7×10-2, ‘black line’ system II, ‘magenta line’ ELS,  ‘green dotted line’ unstable solutions, ‘o’ peak amplitude of transmissibility

b) β=10-4

Force transmissibility of system I, system II and their ELS for displacement excitations with different offset displacements and excitation amplitudes. ‘red line’ system I with y^0=1.7×10-2,  ‘blue line’ system I with y^0=2.7×10-2, ‘black line’ system II, ‘magenta line’ ELS,  ‘green dotted line’ unstable solutions, ‘o’ peak amplitude of transmissibility

c) β=10-4

Force transmissibility of system I, system II and their ELS for displacement excitations with different offset displacements and excitation amplitudes. ‘red line’ system I with y^0=1.7×10-2,  ‘blue line’ system I with y^0=2.7×10-2, ‘black line’ system II, ‘magenta line’ ELS,  ‘green dotted line’ unstable solutions, ‘o’ peak amplitude of transmissibility

d) β=1.2×10-3

Force transmissibility of system I, system II and their ELS for displacement excitations with different offset displacements and excitation amplitudes. ‘red line’ system I with y^0=1.7×10-2,  ‘blue line’ system I with y^0=2.7×10-2, ‘black line’ system II, ‘magenta line’ ELS,  ‘green dotted line’ unstable solutions, ‘o’ peak amplitude of transmissibility

e) β=10-3

By observing Fig. 11, for system I, a decrease in the offset displacement y^0 yields a decrease in both the peak amplitude of the force transmissibility and the resonance frequency, and the transmissibility curve approaches to the transmissibility curve of system II. If the force excitation amplitude is relatively large, the peak amplitude of the force transmissibility and the resonance frequency in system II are larger than that in system I. However, system II will present better isolation performance compared with system I when the excitation frequency exceeds the resonance frequency. For system II, larger excitation amplitude results in both larger peak amplitude of the force transmissibility and the resonance frequency. By combing the characteristics of the stiffness and resonance frequency for different excitation amplitudes in Section 4.2, it is found that the peak amplitude of the force transmissibility of system I decreases at first, increases later, and finally becomes larger than that of the ELS with the increasing excitation amplitude. For the same excitation amplitude, the force transmissibility of system II is larger than that of system I and the force transmissibility of system I is larger than that of the ELS at lower frequencies; but when the excitation frequency is around the resonance frequency of the ELS, the force transmissibility of system II becomes smaller than that of system I and the force transmissibility of system I becomes smaller than that of ELS; and at higher frequencies the force transmissibility of the three systems are almost the same.

Different from the force transmissibility for the force excitation, both system I and II have unbounded absolute displacement and acceleration transmissibility as shown in Fig. 12. For system I, the absolute acceleration transmissibility is always better than the absolute displacement transmissibility. For the same excitation amplitude, the interactive characteristics of both the absolute displacement and acceleration transmissibility among system I, II and their ELS are same with that of the force transmissibility when the excitation frequency is lower or around the resonance frequency of the ELS. However, at higher frequencies the absolute acceleration transmissibility of the three systems are almost the same while the absolute displacement transmissibility of system I is larger than that of system II and the ELS. And the absolute displacement transmissibility of system I increases with the increasing excitation amplitude at higher excitation frequencies.

Fig. 12Absolute displacement and acceleration transmissibility of system I, system II and their ELS for offset displacements with different equilibrium positions and excitation amplitudes. ‘red line’ absolute displacement transmissibility of system I with y^0=1.7×10-2, ‘cyan line’ absolute acceleration transmissibility of system I with y^0=1.7×10-2, ‘blue line’ absolute displacement transmissibility of system II with y^0=2.7×10-2, ‘brown line’ absolute acceleration transmissibility of system II with y^0=2.7×10-2, ‘black line’ system II, ‘magenta line’ ELS, ‘green dotted line’ unstable solutions, ‘o’ peak amplitude of transmissibility

Absolute displacement and acceleration transmissibility of system I, system II and their ELS for offset displacements with different equilibrium positions and excitation amplitudes. ‘red line’ absolute displacement transmissibility of system I with y^0=1.7×10-2, ‘cyan line’ absolute acceleration transmissibility of system I with y^0=1.7×10-2, ‘blue line’ absolute displacement transmissibility  of system II with y^0=2.7×10-2, ‘brown line’ absolute acceleration transmissibility of system II  with y^0=2.7×10-2, ‘black line’ system II, ‘magenta line’ ELS, ‘green dotted line’ unstable  solutions, ‘o’ peak amplitude of transmissibility

a) β=10-3

Absolute displacement and acceleration transmissibility of system I, system II and their ELS for offset displacements with different equilibrium positions and excitation amplitudes. ‘red line’ absolute displacement transmissibility of system I with y^0=1.7×10-2, ‘cyan line’ absolute acceleration transmissibility of system I with y^0=1.7×10-2, ‘blue line’ absolute displacement transmissibility  of system II with y^0=2.7×10-2, ‘brown line’ absolute acceleration transmissibility of system II  with y^0=2.7×10-2, ‘black line’ system II, ‘magenta line’ ELS, ‘green dotted line’ unstable  solutions, ‘o’ peak amplitude of transmissibility

b) β=10-3

Absolute displacement and acceleration transmissibility of system I, system II and their ELS for offset displacements with different equilibrium positions and excitation amplitudes. ‘red line’ absolute displacement transmissibility of system I with y^0=1.7×10-2, ‘cyan line’ absolute acceleration transmissibility of system I with y^0=1.7×10-2, ‘blue line’ absolute displacement transmissibility  of system II with y^0=2.7×10-2, ‘brown line’ absolute acceleration transmissibility of system II  with y^0=2.7×10-2, ‘black line’ system II, ‘magenta line’ ELS, ‘green dotted line’ unstable  solutions, ‘o’ peak amplitude of transmissibility

c) β=4.8×10-3

Absolute displacement and acceleration transmissibility of system I, system II and their ELS for offset displacements with different equilibrium positions and excitation amplitudes. ‘red line’ absolute displacement transmissibility of system I with y^0=1.7×10-2, ‘cyan line’ absolute acceleration transmissibility of system I with y^0=1.7×10-2, ‘blue line’ absolute displacement transmissibility  of system II with y^0=2.7×10-2, ‘brown line’ absolute acceleration transmissibility of system II  with y^0=2.7×10-2, ‘black line’ system II, ‘magenta line’ ELS, ‘green dotted line’ unstable  solutions, ‘o’ peak amplitude of transmissibility

d) β=10-3

Absolute displacement and acceleration transmissibility of system I, system II and their ELS for offset displacements with different equilibrium positions and excitation amplitudes. ‘red line’ absolute displacement transmissibility of system I with y^0=1.7×10-2, ‘cyan line’ absolute acceleration transmissibility of system I with y^0=1.7×10-2, ‘blue line’ absolute displacement transmissibility  of system II with y^0=2.7×10-2, ‘brown line’ absolute acceleration transmissibility of system II  with y^0=2.7×10-2, ‘black line’ system II, ‘magenta line’ ELS, ‘green dotted line’ unstable  solutions, ‘o’ peak amplitude of transmissibility

e) β=10-2

5. Conclusions

In this paper, a QZS isolator fabricated by adding a slotted conical disk spring in parallel with a vertical linear spring is presented. Firstly, the geometrical configuration for designing the unique feature of quasi-zero stiffness is described and the corresponding mathematical modeling is formulated. The configurative parameters are optimized to achieve a wide range of the displacement from the equilibrium position for which the stiffness has a low value and changes slightly.

Secondly, the system of an offset equilibrium position for overloaded or underloaded, and the system of an equilibrium position in which the dynamic stiffness is zero for loaded with an appropriate mass are studied. The primary resonance response for the force excitation and displacement excitation is separately derived by employing the HBM and confirmed by the results of numerical simulation. The ways in which the maximum number of the steady-state values depends on the offset displacement and the amplitude of the two types of excitation are also separately reported.

Thirdly, the frequency response curves of the two kinds of system for the two types of excitation have been plotted. The frequency response curves of their equivalent linear system are plotted in the same figure for comparison. The overloaded system can exhibits purely softening, mixed softening-hardening and purely hardening with the increasing excitation amplitude. It can be summarized that both larger offset displacement and excitation amplitude result in a higher resonance frequency.

Finally, the force transmissibility, the absolute displacement and acceleration transmissibility are chosen to evaluate the isolation for the two kinds of systems and compared with their equivalent linear system. It can be concluded that adding the slotted conical disk spring acting as the negative structure to a linear spring is a feasible way to achieve the wider isolation frequency region and smaller transmissibility. Additionally, an effective way to achieve better isolation performance is applying the QZS isolator in the conditions that the mass differs less with its load capability and the excitation amplitude is not too large.

References

  • Rivin E. I. Passive Vibration Isolation. ASME Press, New York, 2001.
  • Carrella A. Passive vibration isolators with high-static-low-dynamic-stiffness. Ph. D. Thesis, ISVR, University of Southampton, 2008.
  • Ibrahim R. A. Recent advances in nonlinear passive vibration isolators. Journal of Sound and Vibration, Vol. 314, 2008, p. 371-452.
  • Alabuzhev P., Gritchin A., Kim L., Migirenko G., Chon V., Stepanov P. Vibration Protecting and Measuring Systems with Quasi-Zero Stiffness. Hemisphere, New York, 1989.
  • Platus D. L. Negative-stiffness-mechanism vibration isolation systems. Proceedings of SPIE, Vol. 3786, 1999, p. 98-105.
  • Zhang J., Li D., Dong S. An ultra-low frequency parallel connection nonlinear isolator for precision instruments. Key Engineering Materials, Vol. 257-258, 2004, p. 231-236, (in Chinese).
  • Le T. D., Ahn K. K. A vibration isolation system in low frequency excitation region using negative stiffness structure for vehicle seat. Journal of Sound and Vibration, Vol. 330, 2011, p. 631-6335.
  • Liu X., Huang X., Hua H. On the characteristics of a quasi-zero stiffness isolator using Euler buckled beam as negative stiffness corrector. Journal of Sound andVibration, Vol. 332, 2013, p. 3359-3376, (in Chinese).
  • Schremmer G. The slotted coned disk spring. Journal of Engineering for Industry, Trans. of the ASME, 1973, p. 765-770.
  • Frawazi N., Lee J. Y., Oh J. E. A load-displacement prediction for a bended slotted disc using the energy method. Proceedings of the Institution of Mechanical Engineers, Part C: Journal of Mechanical Engineering Science, Vol. 2126, 2012, p. 2126-2137.
  • Zhang Y., Liu H., Wang D. Spring Manual. China Machine Press, Beijing, 2008, (in Chinese).
  • Ravindra B., Mallik A. K. Performance of non-linear vibration isolators under harmonic excitation. Journal of Sound and Vibration, Vol. 170, Issue 3, 1994, p. 325-337.
  • Kovacic I., Brennan M. J., Lineton B. On the resonance response of an asymmetric Duffing oscillator. International Journal of Non-Linear Mechanics, Vol. 43, 2008, p. 858-867.
  • Korn G. A., Korn T. M. Mathematical Handbook for Scientists and Engineers. McGraw-Hill, New York, 1961.
  • Guo P., Lang Z., Peng Z. Analysis and design of the force and displacement transmissibility of nonlinear viscous damper based vibration isolation systems. Nonlinear Dynamics, Vol. 67, 2012, p. 2671-2687.
  • Liu X., Huang X., Zhang Z., Hua H. Influence of excitation amplitude and load on the characteristics of quasi-zero stiffness isolator. Journal of Mechanical Engineering, Vol. 49, Issue 6, 2013, p. 89-94, (in Chinese).

About this article

Received
16 January 2014
Accepted
25 April 2014
Published
30 June 2014
Keywords
quasi-zero stiffness (QZS)
isolator
optimization
overload
frequency response curve (FRC)
transmissibility
Acknowledgements

This research was supported in part by the National Science and Technology Major Project of China (Grant No. 2012ZX10004801-002-013).