# Dynamic Vibration Absorber with Negative Stiffness for Rotor System

Authors: Hongliang Yao; Zidong Chen; Bangchun Wen
Journal: Shock and Vibration
 (2016)

---

## Abstract

To suppress the vibration of a rotor system, a vibration absorber combining negative stiffness with positive stiffness together is proposed in this paper. Firstly, the negative stiffness producing mechanism using ring type permanent magnets is presented and the characteristics of the negative stiffness are analyzed. Then, the structure of the absorber is proposed; the principles and nonlinear dynamic characteristics of the absorber-rotor system are studied numerically. Finally, experiments are carried out to verify the numerical conclusions. The results show that the proposed vibration absorber is effective to suppress the vibration of the rotor system, the nonlinearity of the negatives stiffness affects the vibration suppression effect little, and the negative stiffness can broaden the effective vibration control frequency range of the absorber.


---

## Body

## 1. Introduction

Vibrations of rotor system are the source of many faults in rotating machinery such as rotor-to-stator rub and shaft crack growth. To maintain the safety operation of the rotating machinery, the excessive vibration of rotor system must be controlled.Many studies have been carried out about vibration suppression methods of rotor systems, which can be largely divided into three kinds: vibration control by applying external forces, vibration reduction by adjusting parameters, and vibration absorbing by using dynamic vibration absorbers.The first kind of method is by applying external forces, which are usually produced by active magnetic bearings or active magnetic exciters [1–4]. The advantage of magnetic forces is that they can be acted on the rotating parts contactless, thus changing the structure of the rotor little.The second kind of method is by adjusting parameters of the rotor system such as stiffness, damping, and eccentricity, which can be carried out by using shape memory alloy metals [5], discontinuous springs [6], and so forth. In all these methods, the damping adjusting methods are more popular and many approaches are used, such as squeeze film dampers [7], eddy current damping [8, 9], piezoceramic induced damping [10], magnetorheological elastomer damping [11–13], and friction controlled damping [14]. Eccentricity adjusting is another kind of popular method and is mainly realized by using an automatic ball balancer [15, 16], which can be very efficient when the rotating speed is above the critical speed.The third kind of method is using a vibration absorber or a nonlinear energy sink, by which the vibration of the rotor system can be “absorbed.” For example, the centrifugal pendulum vibration absorbers have long been widely used in torsional vibration suppression [17]. In lateral vibration attenuation, many studies are also carried out; the vibration absorbers are placed on the bearing or foundation [18, 19] or inside the rotating parts [20–22]. As the lateral vibrations of rotor system are in both the horizontal and vertical directions, ring type absorbers are widely used for vibration suppression of the rotor system. For example, in [23], a new kind of viscoelastic vibration absorber was designed for the rotating system and the experiment results verified the vibration suppression effect; in [24], a tunable vibration absorber is designed to suppress regenerative chatter in milling process; in [25], a vibration absorber for a washing machine is proposed, which is mounted on the rotating basket of the machine; in [26], a circular shape vibration absorber to reduce the vibration caused by the imbalance of an optical disk drive is presented; the vibration absorber is speed-dependent and can keep resonance in a specific frequency range.The vibration absorbers often become inefficient in the low frequency range, because it is difficult to reduce the stiffness of absorber to very low values, as the low stiffness will cause large static deflection. Negative stiffness is a good solution for this problem and has been widely used in vibration isolator design. Permanent magnets are a good way to produce negative stiffness and have been widely used in single-direction vibration suppression, such as vibration isolators [27] or dynamic vibration absorbers [28], but, to the authors’ knowledge, no permanent magnet has been used in dynamic vibration absorber design for rotor system.So, in this paper, a novel vibration absorber with negative stiffness produced by permanent magnet is proposed and its fundamental principles and characteristics are studied, and then experiments are carried out to verify the analytical results.


## 2. Negative Stiffness Spring Made Up of Permanent Magnets

### 2.1. Structure of the Negative Stiffness Magnetic Spring

The structure is made up of ring type magnets and is shown as Figure1. The outer magnets and central magnets are mounted in repulsive interaction, and the two central magnets are connected by a connecting part which contains a rolling bearing.Figure 1
Structure of the magnetic spring with negative stiffness.As the outer magnet and the central magnet are in repulsive interaction, a repulsive force in radial direction will occur; the central magnet will leave the equilibrium position and will not return without external force. So the central magnet acts as a spring with negative stiffness.


### 2.2. Magnet Force and Stiffness Calculation Using Equivalent Magnetic Charge Method

The magnet force generated between the two parallel ring type magnets can be calculated by using equivalent magnetic charge method [29, 30].The outer diameter, inner diameter, and thickness of the outer magnets areRD1, Rd1, and b1, and those of the inner magnets are RD2, Rd2, and b2, respectively. The distance between the outer magnet and the central magnet is h. The width of the rolling bearing is B. According to the equivalent magnetic charge method, the point charge of a point r2,α in plane 2 (as shown in Figure 1) is(1)q2=Brr2dr2dα,where Br is the residual flux density.The point charge of a pointr3,β in plane 3 is(2)q3=Brr3dr3dβ.The interaction force between the two points is(3)dF→23=Br24πμ0r2r3dr2dαdr3dβr232r→23,where μ0 is the relative permeability of the permanent magnets.When the central magnets leave the equilibrium point for distancee, the interaction force between planes 2 and 3 is(4)F23=Br24πμ0∫02π∫02π∫Rd1RD1∫Rd2RD2r2r3r3cos⁡β-r1cos⁡α-er3cos⁡β-r1cos⁡α-e2+r3sin⁡β-r1sin⁡α2+h23/2dr2dr3dαdβ.The interaction force between plane 2 and planei can be written as(5)F2i=-1i+1Br24πμ0∫02π∫02π∫Rd1RD1∫Rd2RD2r2riricos⁡β-r2cos⁡α-ericos⁡β-r2cos⁡α-e2+risin⁡β-r2sin⁡α2+h2i23/2dr2dridαdβi=3,…,6,where h23=h, h24=h+b2, h25=h+b2+B, and h26=2h+b2+B.Also, the interaction force between plane 1 and planei can be written as(6)F1i=-1iBr24πμ0∫02π∫02π∫Rd1RD1∫Rd2RD2r1riricos⁡β-r1cos⁡α-ericos⁡β-r1cos⁡α-e2+risin⁡β-r1sin⁡α2+b1+h2i23/2dr1dridαdβi=3,…,6.As the structure is symmetrical, the interaction forces between the right outer magnet and the central magnets are equal to those of the left ones. So the total force between the central magnets and the outer magnets is(7)Fy=2∑i=36F1i+F2i.Equations (4)–(7) can be solved numerically. When the different forces Fy with different e are obtained, the radial stiffness of the negative stiffness magnetic spring can be obtained by (8)km=dFyde.


### 2.3. Numerical Simulation and Stiffness Analysis

The parameters used for numerical simulation are shown in Table1.Table 1
Parameters for numerical simulation of magnetic force and stiffness.
Parameter Value R D 1 50 mm R d 1 25 mm b 1 6 mm R D 2 50 mm R d 2 25 mm b 2 6 mm B 15 mm μ 0 4 π × 1 0 - 7 N/A2 B r 1.34 TWhen the distance between the outer magnet and the central magnet changes from 1 to 5 mm and the distance from the center of central magnet to the equilibrium point changes from 0 to 1 mm, the magnetic force in the radial direction is obtained and is shown in Figure2. It can be seen from Figure 2 that the radial force changes violently when the magnet distance h is small.Figure 2
Magnetic force of the negative stiffness structure.The negative stiffness can be calculated by (8), and the stiffness of different h is shown in Figure 3. It can be seen from Figure 3 that the nonlinearity is weak when h is large, and the nonlinearity increases with the decrease of h and becomes suddenly evident when h is 2.5 mm.Figure 3
Stiffness of the magnetic spring when the distance is from 2.5 to 20 mm.When polynomial fitting is applied, the relationship between the total stiffness and the radial displacement can be written as(9)km=km1+km2e2+km3e4.The parameters of different magnet distances are shown in Table2 and the fitting effects are shown in Figure 4. It can be seen from Table 2 and Figure 4 that the polynomial fitting effect is rather good, and km2 and km3 are rather small when h is larger.Table 2
Parameters for approximate stiffness.
h (mm) k m 1 (N/mm) k m 2 (N/mm) k m 3 (N/mm) 20 −1.212 0.012 −0.001 15 −2.48 0.033 −0.002 10 −5.898 0.112 −0.008 5 −17.28 0.637 −0.054 2.5 −35.88 3.238 −0.434Figure 4
Approximate stiffness and exact stiffness comparison whenh is 15 mm, 10 mm, and 5 mm, respectively.
(a)
h = 15mm (b)
h = 10mm (c)
h = 5mm


## 2.1. Structure of the Negative Stiffness Magnetic Spring

The structure is made up of ring type magnets and is shown as Figure1. The outer magnets and central magnets are mounted in repulsive interaction, and the two central magnets are connected by a connecting part which contains a rolling bearing.Figure 1
Structure of the magnetic spring with negative stiffness.As the outer magnet and the central magnet are in repulsive interaction, a repulsive force in radial direction will occur; the central magnet will leave the equilibrium position and will not return without external force. So the central magnet acts as a spring with negative stiffness.


## 2.2. Magnet Force and Stiffness Calculation Using Equivalent Magnetic Charge Method

The magnet force generated between the two parallel ring type magnets can be calculated by using equivalent magnetic charge method [29, 30].The outer diameter, inner diameter, and thickness of the outer magnets areRD1, Rd1, and b1, and those of the inner magnets are RD2, Rd2, and b2, respectively. The distance between the outer magnet and the central magnet is h. The width of the rolling bearing is B. According to the equivalent magnetic charge method, the point charge of a point r2,α in plane 2 (as shown in Figure 1) is(1)q2=Brr2dr2dα,where Br is the residual flux density.The point charge of a pointr3,β in plane 3 is(2)q3=Brr3dr3dβ.The interaction force between the two points is(3)dF→23=Br24πμ0r2r3dr2dαdr3dβr232r→23,where μ0 is the relative permeability of the permanent magnets.When the central magnets leave the equilibrium point for distancee, the interaction force between planes 2 and 3 is(4)F23=Br24πμ0∫02π∫02π∫Rd1RD1∫Rd2RD2r2r3r3cos⁡β-r1cos⁡α-er3cos⁡β-r1cos⁡α-e2+r3sin⁡β-r1sin⁡α2+h23/2dr2dr3dαdβ.The interaction force between plane 2 and planei can be written as(5)F2i=-1i+1Br24πμ0∫02π∫02π∫Rd1RD1∫Rd2RD2r2riricos⁡β-r2cos⁡α-ericos⁡β-r2cos⁡α-e2+risin⁡β-r2sin⁡α2+h2i23/2dr2dridαdβi=3,…,6,where h23=h, h24=h+b2, h25=h+b2+B, and h26=2h+b2+B.Also, the interaction force between plane 1 and planei can be written as(6)F1i=-1iBr24πμ0∫02π∫02π∫Rd1RD1∫Rd2RD2r1riricos⁡β-r1cos⁡α-ericos⁡β-r1cos⁡α-e2+risin⁡β-r1sin⁡α2+b1+h2i23/2dr1dridαdβi=3,…,6.As the structure is symmetrical, the interaction forces between the right outer magnet and the central magnets are equal to those of the left ones. So the total force between the central magnets and the outer magnets is(7)Fy=2∑i=36F1i+F2i.Equations (4)–(7) can be solved numerically. When the different forces Fy with different e are obtained, the radial stiffness of the negative stiffness magnetic spring can be obtained by (8)km=dFyde.


## 2.3. Numerical Simulation and Stiffness Analysis

The parameters used for numerical simulation are shown in Table1.Table 1
Parameters for numerical simulation of magnetic force and stiffness.
Parameter Value R D 1 50 mm R d 1 25 mm b 1 6 mm R D 2 50 mm R d 2 25 mm b 2 6 mm B 15 mm μ 0 4 π × 1 0 - 7 N/A2 B r 1.34 TWhen the distance between the outer magnet and the central magnet changes from 1 to 5 mm and the distance from the center of central magnet to the equilibrium point changes from 0 to 1 mm, the magnetic force in the radial direction is obtained and is shown in Figure2. It can be seen from Figure 2 that the radial force changes violently when the magnet distance h is small.Figure 2
Magnetic force of the negative stiffness structure.The negative stiffness can be calculated by (8), and the stiffness of different h is shown in Figure 3. It can be seen from Figure 3 that the nonlinearity is weak when h is large, and the nonlinearity increases with the decrease of h and becomes suddenly evident when h is 2.5 mm.Figure 3
Stiffness of the magnetic spring when the distance is from 2.5 to 20 mm.When polynomial fitting is applied, the relationship between the total stiffness and the radial displacement can be written as(9)km=km1+km2e2+km3e4.The parameters of different magnet distances are shown in Table2 and the fitting effects are shown in Figure 4. It can be seen from Table 2 and Figure 4 that the polynomial fitting effect is rather good, and km2 and km3 are rather small when h is larger.Table 2
Parameters for approximate stiffness.
h (mm) k m 1 (N/mm) k m 2 (N/mm) k m 3 (N/mm) 20 −1.212 0.012 −0.001 15 −2.48 0.033 −0.002 10 −5.898 0.112 −0.008 5 −17.28 0.637 −0.054 2.5 −35.88 3.238 −0.434Figure 4
Approximate stiffness and exact stiffness comparison whenh is 15 mm, 10 mm, and 5 mm, respectively.
(a)
h = 15mm (b)
h = 10mm (c)
h = 5mm


## 3. Dynamic Vibration Absorber with Magnetic Spring

### 3.1. Structure of the Dynamic Vibration Absorber

The structure of the absorber is shown in Figure5(a) and that of the whole test rig is shown in Figure 5(b). The central magnets are fixed on the connecting part. The connecting part is mounted on the rolling bearing. With the rolling bearing, the absorber will not rotate with the shaft together, so the electromagnetic damping between the central magnet and the frame of the absorber can be avoided. The outer magnets are fixed on the aluminum plates, which are connected by beams with circular section. The circular section beams are connected to the supporting frame by linear bearings, which ensure that the beams have single-direction stiffness in the lateral vibration plane: the vertical beam has stiffness in horizontal direction and the horizontal beam has stiffness in vertical direction, as shown in Figure 5(c). The beams act as positive stiffness springs and the magnet structure acts as negative stiffness spring.Figure 5
Structure of the dynamic vibration absorber and the rotor-absorber system.
(a)
Structure of the absorber (b)
Structure of the whole test rig (c)
Schematic diagram of the beam installmentThe dynamic model of the unbalanced rotor-absorber system is shown in Figure6, and the dynamic equations of the system are (10)m100m2r¨1r¨2+c1+cm-cm-cmcm+c3r˙1r˙2+k1+km1-km1-km1km1+k3r1r2+km2r1-r23+km3r1-r25-km2r1-r23-km3r1-r25=meδω2eiωt0,where m1, k1, and c1 are the mass, stiffness, and damping of the rotor, respectively. km and cm are the stiffness and damping of the negative stiffness structure. m2 is the mass of the part of the absorber attached to the rotor and k3 is the stiffness of the circular beam. r1 and r2 are the displacement vectors of the rotor and the absorber in the radial direction, respectively. me and δ are the eccentric mass and eccentric distance of the disc, respectively. ω is the rotating speed of the rotor.Figure 6
Dynamic model of the rotor-absorber system.


### 3.2. Numerical Simulation of the Dynamic Responses of the Rotor-Absorber System

When the nonlinearity of the negative stiffness is considered, numerical simulation must be applied to study the responses of the rotor-absorber system. The Incremental Harmonic Balance (IHB) method [31] is applied in the following analysis. To trace the unstable periodic responses, the arc-length method [32] is combined with the IHB method.The parameters of the magnets are assumed the same as those in Table1, and the other parameters of the rotor-absorber system are shown in Table 3.Table 3
Parameters of the rotor-absorber system.
Parameter Value m 1 1 kg m 2 0.7 kg k 1 9.5 × 1 0 4 N/m c m 0 Nm/s c 1 10 Nm/s c 3 6 Nm/sThe stiffnesskm is adjusted by the distance h between the outer magnet and the central magnet, and k3 is adjusted by the lengths of the circular section beams. For the cantilever beam with length l and diameter d, the bending stiffness can be obtained by (11)k3=3πEd464l3,where E is the elastic modulus of material. When the diameter of the beam is 5 mm and the length is 67 mm, 64 mm, and 61 mm, the stiffness k3 is about 6×104 N/m, 7×104 N/m, and 8×104 N/m, respectively.Firstly, the amplitude-frequency response curve of the rotor system without absorber whenδ is 0.015 mm is shown in Figure 7. It can be seen from Figure 7 that the first-order critical speed is about 50 Hz and the largest amplitude is about 0.45 mm.Figure 7
Responses of the rotor system.Then the absorber is mounted on the rotor system, and the frequency-amplitude curves of the rotor system and the absorber whenδ is 0.015 mm are studied. As the nonlinearity of the negative stiffness differs with the distance h, here the cases with different distance are studied. Firstly, the responses of the system when h is 10 mm are calculated. Figures 8(a), 8(b), and 8(c) are the responses of the rotor and absorber when k3 is 6×104 N/mm, 7×104 N/mm, and 8×104 N/mm, respectively. Then, the responses of the system when h is 5 mm are calculated and shown in Figure 9.Figure 8
Responses of the rotor-absorber system whenh is 10 mm.
(a)
k 3 = 6 × 1 0 4 N/mm (b)
k 3 = 7 × 1 0 4 N/mm (c)
k 3 = 8 × 1 0 4 N/mmFigure 9
Responses of the rotor-absorber system whenh is 5 mm.
(a)
k 3 = 6 × 1 0 4 N/mm (b)
k 3 = 7 × 1 0 4 N/mm (c)
k 3 = 8 × 1 0 4 N/mmStudying Figures8 and 9, some conclusions can be obtained:(1)
The vibrations of the rotor can be suppressed when the absorber is added. The maximum amplitude of the rotor is much smaller than that of the rotor system without absorber. The antiresonant frequency changes with the change ofk3, so when the rotating speed of the rotor is determined, we can change k3 to keep the small vibration of the rotor.(2)
The vibration suppression effect is better when the distanceh is smaller. The residual vibrations in Figure 9 are much smaller than those in Figure 8 and the low vibration regions are much larger.(3)
The nonlinearity affects the suppression effect little, especially near the antiresonant frequency range. Although there exists nonlinearity in the magnetic stiffness, the response curves show little “hard” or “soft” characteristic, because the vibration of the rotor system is small.When the vibration is larger, strong nonlinearity phenomenon may appear in the resonance region, as shown in Figure10, in which k3 is 7×104 N/mm and the eccentricity of the rotor is δ=0.1mm. But this case is very rare, as the vibration amplitude cannot be so great in real rotor systems in normal condition. Also, although the nonlinearity becomes stronger, the vibration suppression effects near antiresonant frequency remain the same.Figure 10
Responses of the rotor-absorber system whenh is 5 mm and δ is 0.1 mm.


## 3.1. Structure of the Dynamic Vibration Absorber

The structure of the absorber is shown in Figure5(a) and that of the whole test rig is shown in Figure 5(b). The central magnets are fixed on the connecting part. The connecting part is mounted on the rolling bearing. With the rolling bearing, the absorber will not rotate with the shaft together, so the electromagnetic damping between the central magnet and the frame of the absorber can be avoided. The outer magnets are fixed on the aluminum plates, which are connected by beams with circular section. The circular section beams are connected to the supporting frame by linear bearings, which ensure that the beams have single-direction stiffness in the lateral vibration plane: the vertical beam has stiffness in horizontal direction and the horizontal beam has stiffness in vertical direction, as shown in Figure 5(c). The beams act as positive stiffness springs and the magnet structure acts as negative stiffness spring.Figure 5
Structure of the dynamic vibration absorber and the rotor-absorber system.
(a)
Structure of the absorber (b)
Structure of the whole test rig (c)
Schematic diagram of the beam installmentThe dynamic model of the unbalanced rotor-absorber system is shown in Figure6, and the dynamic equations of the system are (10)m100m2r¨1r¨2+c1+cm-cm-cmcm+c3r˙1r˙2+k1+km1-km1-km1km1+k3r1r2+km2r1-r23+km3r1-r25-km2r1-r23-km3r1-r25=meδω2eiωt0,where m1, k1, and c1 are the mass, stiffness, and damping of the rotor, respectively. km and cm are the stiffness and damping of the negative stiffness structure. m2 is the mass of the part of the absorber attached to the rotor and k3 is the stiffness of the circular beam. r1 and r2 are the displacement vectors of the rotor and the absorber in the radial direction, respectively. me and δ are the eccentric mass and eccentric distance of the disc, respectively. ω is the rotating speed of the rotor.Figure 6
Dynamic model of the rotor-absorber system.


## 3.2. Numerical Simulation of the Dynamic Responses of the Rotor-Absorber System

When the nonlinearity of the negative stiffness is considered, numerical simulation must be applied to study the responses of the rotor-absorber system. The Incremental Harmonic Balance (IHB) method [31] is applied in the following analysis. To trace the unstable periodic responses, the arc-length method [32] is combined with the IHB method.The parameters of the magnets are assumed the same as those in Table1, and the other parameters of the rotor-absorber system are shown in Table 3.Table 3
Parameters of the rotor-absorber system.
Parameter Value m 1 1 kg m 2 0.7 kg k 1 9.5 × 1 0 4 N/m c m 0 Nm/s c 1 10 Nm/s c 3 6 Nm/sThe stiffnesskm is adjusted by the distance h between the outer magnet and the central magnet, and k3 is adjusted by the lengths of the circular section beams. For the cantilever beam with length l and diameter d, the bending stiffness can be obtained by (11)k3=3πEd464l3,where E is the elastic modulus of material. When the diameter of the beam is 5 mm and the length is 67 mm, 64 mm, and 61 mm, the stiffness k3 is about 6×104 N/m, 7×104 N/m, and 8×104 N/m, respectively.Firstly, the amplitude-frequency response curve of the rotor system without absorber whenδ is 0.015 mm is shown in Figure 7. It can be seen from Figure 7 that the first-order critical speed is about 50 Hz and the largest amplitude is about 0.45 mm.Figure 7
Responses of the rotor system.Then the absorber is mounted on the rotor system, and the frequency-amplitude curves of the rotor system and the absorber whenδ is 0.015 mm are studied. As the nonlinearity of the negative stiffness differs with the distance h, here the cases with different distance are studied. Firstly, the responses of the system when h is 10 mm are calculated. Figures 8(a), 8(b), and 8(c) are the responses of the rotor and absorber when k3 is 6×104 N/mm, 7×104 N/mm, and 8×104 N/mm, respectively. Then, the responses of the system when h is 5 mm are calculated and shown in Figure 9.Figure 8
Responses of the rotor-absorber system whenh is 10 mm.
(a)
k 3 = 6 × 1 0 4 N/mm (b)
k 3 = 7 × 1 0 4 N/mm (c)
k 3 = 8 × 1 0 4 N/mmFigure 9
Responses of the rotor-absorber system whenh is 5 mm.
(a)
k 3 = 6 × 1 0 4 N/mm (b)
k 3 = 7 × 1 0 4 N/mm (c)
k 3 = 8 × 1 0 4 N/mmStudying Figures8 and 9, some conclusions can be obtained:(1)
The vibrations of the rotor can be suppressed when the absorber is added. The maximum amplitude of the rotor is much smaller than that of the rotor system without absorber. The antiresonant frequency changes with the change ofk3, so when the rotating speed of the rotor is determined, we can change k3 to keep the small vibration of the rotor.(2)
The vibration suppression effect is better when the distanceh is smaller. The residual vibrations in Figure 9 are much smaller than those in Figure 8 and the low vibration regions are much larger.(3)
The nonlinearity affects the suppression effect little, especially near the antiresonant frequency range. Although there exists nonlinearity in the magnetic stiffness, the response curves show little “hard” or “soft” characteristic, because the vibration of the rotor system is small.When the vibration is larger, strong nonlinearity phenomenon may appear in the resonance region, as shown in Figure10, in which k3 is 7×104 N/mm and the eccentricity of the rotor is δ=0.1mm. But this case is very rare, as the vibration amplitude cannot be so great in real rotor systems in normal condition. Also, although the nonlinearity becomes stronger, the vibration suppression effects near antiresonant frequency remain the same.Figure 10
Responses of the rotor-absorber system whenh is 5 mm and δ is 0.1 mm.


## 4. Experiments

### 4.1. Experiment Setup

The experiments are carried out on the Bently Nevada test rig; the absorber is mounted on the single-disc rotor system, as shown in Figure11. The lengths of the vertical and horizontal beam can be changed by changing the positions and heights of the manual tuning platforms. The parameters of the rotor system and the absorber are almost the same as those in Section 3. The vibration responses of the rotor are collected by eddy current transducer, and those of the absorber are collected by acceleration transducers. To prevent the rotation of the central magnet which is mounted on rolling bearings, a flexible copper wire is used, which holds the central magnet on one end and is fixed to the ground on the other end.Figure 11
Experimental setup: the photo.
(a)
Structure of the absorber (b)
Whole test rig


### 4.2. Experiment Results

To test the vibration suppression effect of the absorber, the vibration signals of the run-up course are recorded and then the amplitude-frequency response curves of the rotor and the absorber are obtained. As the rotor system is symmetry in horizontal and vertical direction, only the responses of the vertical direction are studied.Firstly, the amplitude-frequency response curve of the rotor system without absorber is obtained and shown in Figure12. It can be seen that the natural frequency of the primary rotor system is about 50 Hz.Figure 12
Response of the rotor system without the absorber.Then, the absorber is mounted on the rotor system and the amplitude-frequency response curves of rotor system and the absorber are obtained. The lengthl of the absorber beams can be adjusted. Figures 13(a)–13(c) present the experiment results of the rotor and absorber responses when h is 10 mm and k3 is about 6×104 N/mm, 7×104 N/mm, and 8×104 N/mm, respectively. Figures 14(a)–14(c) are corresponding results when h is 5 mm.Figure 13
Responses of the rotor-absorber system with differentk3 when h is 10 mm.
(a)
k 3 = 6 × 1 0 4 N/mm (b)
k 3 = 7 × 1 0 4 N/mm (c)
k 3 = 8 × 1 0 4 N/mmFigure 14
Responses of the rotor-absorber system with differentk3 when h is 5 mm.
(a)
k 3 = 6 × 1 0 4 N/mm (b)
k 3 = 7 × 1 0 4 N/mm (c)
k 3 = 8 × 1 0 4 N/mmComparing the experimental results with the numerical simulation results, it can be seen that there are still some differences between the experimental and the simulation results. There are additional resonance peaks in experimental results that are not present in the simulation results. The possible reason is that the additional resonance peaks are caused by rotational degree of freedom of the absorber. In the numerical simulations, the absorber is only modeled to move inx and y directions constrained by the beams, but the possibility of angular oscillations is not included. So the affections to the movement by the rotational degree of freedom are neglected.But the main principles of the absorber can be verified after comparing the experimental results with the numerical simulation results:(1)
The vibration suppression effect can be verified by the experimental results. In all the cases, the vibration amplitudes near the antiresonant frequency are very low. Also, the antiresonant frequency of the rotor system moves to low frequency when the length of the beam increases, so the primary system can remain in low vibration state so far as the beam length is adjustable.(2)
The vibration suppression effect improves when the distance between the central magnet and the outer magnet reduced, which can be verified by comparing Figure13 with Figure 14.(3)
The nonlinearity affects the suppression effect little, as the amplitude-frequency response curves show little “hard” or “soft” effect.As the negative stiffness is adopted, the natural frequency of the absorber can be reduced, as it is determined by the combination of the negative stiffness and positive stiffness. This means that the negative stiffness can broaden the effective frequency range of the absorber. Also, as the adjustable stiffness connects to the ground, it will reduce the mass of the primary system, and it will be easy to add motor driven equipment and make it an automatic adjustable absorber.


## 4.1. Experiment Setup

The experiments are carried out on the Bently Nevada test rig; the absorber is mounted on the single-disc rotor system, as shown in Figure11. The lengths of the vertical and horizontal beam can be changed by changing the positions and heights of the manual tuning platforms. The parameters of the rotor system and the absorber are almost the same as those in Section 3. The vibration responses of the rotor are collected by eddy current transducer, and those of the absorber are collected by acceleration transducers. To prevent the rotation of the central magnet which is mounted on rolling bearings, a flexible copper wire is used, which holds the central magnet on one end and is fixed to the ground on the other end.Figure 11
Experimental setup: the photo.
(a)
Structure of the absorber (b)
Whole test rig


## 4.2. Experiment Results

To test the vibration suppression effect of the absorber, the vibration signals of the run-up course are recorded and then the amplitude-frequency response curves of the rotor and the absorber are obtained. As the rotor system is symmetry in horizontal and vertical direction, only the responses of the vertical direction are studied.Firstly, the amplitude-frequency response curve of the rotor system without absorber is obtained and shown in Figure12. It can be seen that the natural frequency of the primary rotor system is about 50 Hz.Figure 12
Response of the rotor system without the absorber.Then, the absorber is mounted on the rotor system and the amplitude-frequency response curves of rotor system and the absorber are obtained. The lengthl of the absorber beams can be adjusted. Figures 13(a)–13(c) present the experiment results of the rotor and absorber responses when h is 10 mm and k3 is about 6×104 N/mm, 7×104 N/mm, and 8×104 N/mm, respectively. Figures 14(a)–14(c) are corresponding results when h is 5 mm.Figure 13
Responses of the rotor-absorber system with differentk3 when h is 10 mm.
(a)
k 3 = 6 × 1 0 4 N/mm (b)
k 3 = 7 × 1 0 4 N/mm (c)
k 3 = 8 × 1 0 4 N/mmFigure 14
Responses of the rotor-absorber system with differentk3 when h is 5 mm.
(a)
k 3 = 6 × 1 0 4 N/mm (b)
k 3 = 7 × 1 0 4 N/mm (c)
k 3 = 8 × 1 0 4 N/mmComparing the experimental results with the numerical simulation results, it can be seen that there are still some differences between the experimental and the simulation results. There are additional resonance peaks in experimental results that are not present in the simulation results. The possible reason is that the additional resonance peaks are caused by rotational degree of freedom of the absorber. In the numerical simulations, the absorber is only modeled to move inx and y directions constrained by the beams, but the possibility of angular oscillations is not included. So the affections to the movement by the rotational degree of freedom are neglected.But the main principles of the absorber can be verified after comparing the experimental results with the numerical simulation results:(1)
The vibration suppression effect can be verified by the experimental results. In all the cases, the vibration amplitudes near the antiresonant frequency are very low. Also, the antiresonant frequency of the rotor system moves to low frequency when the length of the beam increases, so the primary system can remain in low vibration state so far as the beam length is adjustable.(2)
The vibration suppression effect improves when the distance between the central magnet and the outer magnet reduced, which can be verified by comparing Figure13 with Figure 14.(3)
The nonlinearity affects the suppression effect little, as the amplitude-frequency response curves show little “hard” or “soft” effect.As the negative stiffness is adopted, the natural frequency of the absorber can be reduced, as it is determined by the combination of the negative stiffness and positive stiffness. This means that the negative stiffness can broaden the effective frequency range of the absorber. Also, as the adjustable stiffness connects to the ground, it will reduce the mass of the primary system, and it will be easy to add motor driven equipment and make it an automatic adjustable absorber.


## 5. Conclusions

A dynamic absorber combining negative stiffness with conventional positive stiffness for rotor system vibration suppression is proposed in this paper. The principles and fundamental characteristics of the absorber are studied analytically and experimentally and some conclusions are drawn as follows:(1)
The proposed dynamic absorber is effective, and the vibration suppression effect can be adjusted by changing the magnitude of negative stiffness and stiffnessk3. The antiresonant frequency increases with the increase of stiffness k3; the vibration suppression effect improves when the magnitude of the negative stiffness is greater.(2)
The nonlinearity of the negative stiffness is weak and it affects the vibration suppression effect little in normal condition.(3)
The absorber can keep efficiency at low frequency as the negative stiffness lowers the natural frequency of the absorber.



---

*Source: 5231704-2016-12-07.xml*