A study of a nonlinear magnetic vibration isolator with quasi-zero-stiffness

A novel quasi-zero stiffness nonlinear magnetic isolator is designed by using magnetic springs in parallel with linear positive stiffness spring. Through the static analysis, the mathematical expressions of force-displacement-current and stiffness-displacement-current of the system are established, and the necessary conditions for the normal function and geometric parameters of the system are obtained. The nonlinear dynamic equations of the system under external excitation force are established. The amplitude and frequency characteristics of the system are deduced by harmonic balance method. The influence of system parameters and external excitation amplitude on the dynamic characteristics of the system is analyzed. Results showed that the vibration of the system can be controlled by controlling the electromagnet current or increasing the damping coefficient while the external excitation amplitude is kept in a certain of range. The researches provided a theoretical guidance for the design and application of the new quasi-zero stiffness nonlinear magnetic isolation system.


Introduction
At present, low frequency vibration isolation is a major research hotspot and challengeable problem in vibration engineering field.The vibration isolation modes include passive vibration and active vibration isolation.Due to the complex structure and high cost of active vibration isolation, the passive vibration isolation is a common method.According to the theory of vibration isolation, when the excitation frequency is greater than √2 times of the natural frequency, the system starts to isolate the harmful vibration.Therefore, the vibration isolation frequency band is wider when the natural frequency is lower.However, when the natural frequency is excessively lower, the greater static displacement will be produced.Then the contradiction of the lower natural frequency and the poorer stability resides in the system.In order to solve the problem, an isolator so-called quasi-zero stiffness (QZS) vibration isolator is proposed [1,2].The QZS isolator can achieve ultra-low stiffness, zero stiffness, or negative stiffness characteristics [3,4] by designing the appropriate structural parameters.The quasi-zero stiffness (QZS) vibration isolators are usually realized by connecting a positive-stiffness mechanism with a negative-stiffness mechanism.In previous researches [5,6], a typical configuration of the QZS isolator is shown in Fig. 1.The vertical spring acts as a positive-stiffness and the oblique springs act as a negativestiffness.When an appropriate mass is loaded on the isolator, the springs begin to compress until the oblique springs reach the equilibrium position.At this point, the dynamic stiffness is zero if the system parameters are appropriate.
Quasi-zero stiffness vibration isolator is widely studied by a large number of scholars.Carrella et al. [7] designed a vibration isolation system consisting of three springs, obtained the nonlinear force transfer characteristics of the vibration isolation system, and explained the concept of high static stiffness and low dynamic stiffness.Kovacic et al. [8][9][10] optimized the structural parameters of the three-spring quasi-zero stiffness isolator and studied its bifurcation characteristics.Zhou et al. [11,12] used the electromagnet, instead of the permanent magnet, to design a tunable magnetic isolator, which can adjust the system characteristic parameters by changing the electromagnet current.Thanh et al. [13] studied the horizontal spring by connecting the rod parallel to the vertical

Static analysis
The typical design of a quasi-zero stiffness mechanism generally consists of a vertical spring with positive stiffness and two oblique springs with negative stiffness [5][6][7].Considering a new model of the QZS isolator shown in Fig. 2(a), the quasi-zero stiffness of system is achieved by combining a vertical linear spring with two symmetrically inclined magnetic springs.The magnetic spring set consists of electromagnet and armature.The hollow ball is connected with the armature and the connecting rod can move freely in the guide.The two connecting rods and mechanical spring are hinged at the position .The vertical linear spring's stiffness is .The length of the connecting rod is .The compressed length of the vertical spring is ℎ when the system reaches the static equilibrium position.The coordinate defines the displacement from the point .
The magnetic force, is given by [15-17]: where 0 is air permeability, is coil number, is magnet pole area, 0 is air gap between armature and electromagnet, ∆ is armature displacement and is coil current.The relational expression between the applied force and the displacement is given by: Noting that ∆ = − (ℎ − ) − √ − ℎ , Eq. ( 2) can be written as: Eq. ( 3) becomes: Differentiating Eq. ( 4), the non-dimensional stiffness of the system can be obtained: If Eq. ( 5) is considered at the static equilibrium position = 1 and set to zero, the value of that ensures quasi-zero stiffness is: According to the Fig. 1 and Eq. ( 3), the geometric parameters of the system must be satisfied: Based on the analysis above, the parameters satisfying the Eq. ( 7) are listed in Table 1.The parameters designed in this paper are as follows: 0 = 4 ×10 -7 H/m, = 470, = 0.009 m 2 , = 30000 N/m, ℎ = 0.03 m, = 0.07 m, 0 = 0.008 m.The current range of the coil is: 0 ≤ ≤ 5 A, then the non-dimensional current range is: 0 ≤ ≤ 0.2.The relationship between the non-dimensional force, current and the non-dimensional displacement is plotted using Eq. ( 4), which is illustrated in Fig. 3.As shown in the Fig. 3, the nonlinearity of the system is getting stronger with the increase of the current.The relationship between the non-dimensional stiffness, current and the non-dimensional displacement of the system given by Eq. ( 5) is shown in Fig. 4, which indicates that the magnitude of the stiffness varies with the current.When the current is small, the system always exhibits positive stiffness.As the current increases to , the stiffness of the system at the static equilibrium position is zero, but in the other range, the stiffness is always positive.When the current is too large, the magnetic spring plays a leading role, and thus the system shows negative stiffness in some areas.
The non-dimensional stiffness is obtained by differentiating Eq. ( 8): where: The comparison between exact expressions given by Eq. ( 4) and Eq. ( 5) and approximate expressions given by Eq. ( 8) and Eq. ( 9) of non-dimensional force and non-dimensional stiffness is shown in Fig. 5.With the increase of the displacement, the error between the exact and the approximate expression increases.When the displacement of the system at the static equilibrium position is small, the error between the exact expression and the approximate one is small, and thus the Taylor series expansion can be applied to simulate the exact expression.

Dynamic modeling
The dynamical model of the system under external excitation is shown in Fig. 6, the nondimensional motion differentiate equation of the system is: where: The solution of the system is: For convenience, let = 1: Neglect the operation of the higher order term, the coefficients of sine and cosine in Eq. ( 12) should be zero: Introduce = cos( ), = sin( ), where is the response amplitude of the system.Eq. ( 13) and Eq. ( 14) become: sin( ) The amplitude-frequency relation expression can be obtained by combining Eq. ( 15) with Eq. (1š) and using sin ( ) + cos ( ) = 1: The transmitted force, shown in Fig. 6 is given by: Using the Harmonic theory, suppose the response of the system has the form ( ) = cos( + ) , substitute it to Eq. ( 18) and omit the higher harmonic terms.The magnitude of the transmitted force is obtained: Thus, the force transmissibility is given:

Effects of system parameters and excitation amplitude on dynamical characteristics
Based on the analysis above, the influences of system parameters and excitation amplitude on dynamic characteristics are shown in Fig. 7-9.Three cases are analyzed in detail as follows: (1) The damping ratio is varied.When the system parameters ̂= 0.04, = 0.1, the amplitude-frequency characteristic given by Eq. ( 19) and the force transmissibility given by Eq. ( 20) for several values of the damping ratio are plotted in Fig. 7.The nonlinear characteristic of the system is gradually weakened, as shown in Fig. 7(a).With the increase of damping coefficient, the resonance peak and the jumping frequency decreases, but the jump-up frequency is not substantially varied.As seen from the Fig. 7(b), the vibration isolation effect of the system is enhanced in the low frequency band, and the effect is weakened in the high frequency range.(2) The excitation amplitude parameter is varied.When the system parameters = 0.04, = 0.4, the amplitude-frequency characteristic given by Eq. ( 17) and the force transmissibility given by Eq. ( 22) for several values of the excitation amplitude parameter are plotted in Fig. 8.The system exhibits obvious nonlinear characteristics, as shown in Fig. 8(a).With the increase of excitation amplitude, the resonance peak, the jump-up and the jump-down frequency increase.As can be seen from the Fig. 8(b), the range of low frequency vibration isolation becomes narrower.When the external excitation frequency is higher than the jump-down frequency, the force transmissibility curve coincides.So, the vibration isolation effect is the same in this bandwidth.(3) The current is varied.When the system parameters = 0.1, = 0.4,When the system parameters = 0.04, = 0.1, the amplitude-frequency characteristic given by Eq. ( 17) and the force transmissibility given by Eq. ( 20) for several values of the current are plotted in Fig. 9.The nonlinear characteristic of the system is gradually enhanced, as shown in Fig. 9(a).With the increase of the current, the resonance peak decreases and the vibration isolation starting frequency decreases.As seen from the Fig. 9(b), the vibration isolation performance decreases gradually.According to the analysis above, based on the proper control of the external excitation amplitude, it is possible to obtain a good vibration isolation performance by controlling the current of the system or increasing the damping coefficient of the system.

Conclusions
In this paper, a new nonlinear electromagnetic vibration isolator with quasi-zero-stiffness is designed.Through the static analysis, the mathematical model of force-displacement-current and stiffness-displacement-current are derived.The conditions that the geometrical parameters of the system should be satisfied to ensure proper operation of the isolator have been analyzed.The geometric dimension of the isolator is designed under these constraints.It is found that changing the current can make the system at the equilibrium position to achieve quasi-zero stiffness characteristics.Through the dynamic analysis, the study shows that under the condition of certain damping coefficient and external excitation, the appropriate current control range can be selected to restrain the amplitude of vibration effectively.

2 .
a) Structural diagram of the new QZS isolator in the initial position b) Schematic of the QZS isolator when system activated c) The model of the proposed isolator Fig. Schematic of the QZS isolator Introducing the non-dimensional parameters:

Fig. 3 .Fig. 4 .
Fig. 3. Non-dimensional force of the system: a) the surface graph of the non-dimensional force-currentdisplacement, b) the non-dimensional force-displacement curves for different values of ̂

Fig. 5 .
Fig. 5. Comparison between exact and approximate expressions of non-dimensional force and non-dimensional stiffness ( = ): a) the non-dimensional force-displacement curves, b) the non-dimensional stiffness-displacement curves

7 .
force transmissibility curve Fig.The dynamical characteristic of the system when damping ratio is varied

Fig. 8 .
The dynamical characteristic of the system when parameter is varied

Fig. 9 .
The dynamical characteristic of the system when current is varied

Table 1 .
The parameters of the QZS isolator A STUDY OF A NONLINEAR MAGNETIC VIBRATION ISOLATOR WITH QUASI-ZERO-STIFFNESS. PAN SU, JIECHANG WU, SHUYONG LIU, YUAN FANG