# Modeling, Testing, and Characteristic Analysis of a Planetary Flywheel Inerter

Authors: Zheng Ge; Weirui Wang
Journal: Shock and Vibration
 (2018)

---

## Abstract

We propose the planetary flywheel inerter, which is a new type of ball screw inerter. A planetary flywheel consists of several planetary gears mounted on a flywheel bracket. When the flywheel bracket is driven by a screw and rotating, each planetary gear meshing with an outer ring gear generates a compound motion composed of revolution and rotation. Theoretical analysis shows that the output force of the planetary flywheel inerter is proportional to the relative acceleration of one terminal of the inerter to the other. Optimizing the gear ratio of the planetary gears to the ring gear allows the planetary flywheel to be lighter than its traditional counterpart, without any loss on the inertance. According to the structure of the planetary flywheel inerter, nonlinear factors of the inerter are analyzed, and a nonlinear dynamical model of the inerter is established. Then the parameters in the model are identified and the accuracy of the model is validated by experiment. Theoretical analysis and experimental data show that the dynamical characteristics of a planetary flywheel inerter and those of a traditional flywheel inerter are basically the same. It is concluded that a planetary flywheel can completely replace a traditional flywheel, making the inerter lighter.


---

## Body

## 1. Introduction

Mechanical vibration is widely found in nature as well as construction, vehicles, manufacture process, and other engineering fields. For a long time, researchers have deeply carried out studies on how to control the mechanical vibration. The traditional vibration control system is composed of two basic vibration control elements, spring and damper [1].In 2002, Smith proposed a new type of two-terminal mechanical device called inerter, whose output force is proportional to the relative acceleration of one of the inerter’s terminals to the other, as sketched in Figure1. x1 and x2 are the two terminals’ displacement of the inerter. The proportional coefficient b in the figure is called theinertance and has units of kilograms [2].(1)F=bx¨1-x¨2.Figure 1
Schematic diagram of an inerter.On the one hand, an inerter has two ends, which cannot be implemented by a mass block. On the other hand, the device can convert a linear motion into the rotational motion of a flywheel, which generates much more equivalent inertial mass than the gravitational mass of the flywheel [3]. For example, if the mass block in a dynamic vibration absorber is replaced by an inerter, the overall weight can be reduced by more than 90% [4]. The two characters of inerters have considerably enriched the design theory of mechanical vibration networks. As the effect of the inerter in the mechanical network is equivalent to that of the capacitor in the electrical network, the theoretical results of the electrical network synthesis can be applied to the mechanical network design. The realization problem of mechanical impedance with inerter was therefore studied in depth [5]. A number of special classes of admittances realization techniques have been applied to vehicle suspension designs [6, 7]. Various types of new passive suspension were proposed [8, 9]. Nowadays, the mechanical network design technology with inerter has also been applied to many fields such as construction [10] and robot [11].Currently, there are multiple structural variants of inerters, ball screw type [12], rack and pinion type [13], hydraulic motor type [14], fluid type [15, 16], and so on. The ball screw inerter has less friction. And the nut preload of the ball screw inerter helps to reduce the impact of the backlash, gaining a better comprehensive performance [17].Due to the unique advantages of ball screw inerters, the current theoretical research on them is relatively deep. The ball screw has become the structural basis of many new types of inerters. Papageorgiou and Smith established a dynamical model of a ball screw inerter with backslash and nonlinear elastic force and analyzed the influence of nonlinear factors on the dynamical characteristics of the inerter [18]. Hu et al. built a semiactive inerter by adding an electrically adjusted centrifugal block to a ball screw inerter, which provided a method of changing the inertance dynamically [19]. Pires et al. connected a motor to the end of the screw in an inerter and set up the electrical elements at the motor terminals, forming a flexible configuration of electromechanical network [20]. Sun et al. studied the suspension characteristics based on a nonlinear ball screw inerter model, and the parameters of the model were obtained by step-by-step identification of friction parameters and linear parameters [21]. Gonzalez-Buelga et al. researched the influence of nonlinear factors on the characteristics of tuned inerter damper (TID) and verified the theory by experiments. [22].The flywheel is the main inertial component of an inerter. An inerter can accomplish dynamical characteristics as a heavy object using a flywheel of small gravitational mass. The greater the ratio of the rotational inertia to the flywheel gravitational mass, the more evident the lightweight effect of the inerter benefits. A traditional flywheel is usually designed as a disk structure. Since the flywheel diameter is normally restricted by the surroundings, the traditional flywheel inertia is largely controlled by and is proportional to the gravitational mass of the flywheel. At present, the research on the flywheel for an inerter is primarily focused on how to dynamically change the inertance by diversification on flywheel structures, which is helpful to control the vibration by active or semiactive method through inerter [19, 23–25]. However, for passive inerter with fixed inertance, there are few designs for lightening the flywheel. In order to further improve the ratio of the rotational inertia to the flywheel gravitational mass (the inertia-to-mass ratio) of the inerter, we propose a planetary flywheel, which utilizes the revolution and rotation of the planetary gears to increase the inertia-to-mass ratio.This paper is organized as follows: Section2 illustrates the design and a prototype of the planetary flywheel inerter. An ideal model of the new type inerter is established, and the advantages of the planetary flywheel were revealed by comparing it with a traditional flywheel in Section 3. Section 4 analyzes the nonlinear factors of a planetary flywheel and a ball screw inerter, respectively. And then a dynamical model of the planetary flywheel inerter is established. Section 5 focuses on the experimental setup, the model parameter identification, and the model validation. The conclusions are drawn in Section 6.


## 2. Design and Prototype Production of a Planetary Flywheel Inerter

### 2.1. Introduction of the Design and Principle of the Inerter

As exhibited in Figure2, the piston rod and the housing are the two moving terminals of the planetary flywheel inerter. The left end of the piston rod is threaded, and the housing has a mounting flange on it for connection. When the relative movement occurs between the two terminals of the inerter, the piston rod moves the nut along the linear bearing and the linear guide rail. The nut then converts the linear motion into the rotary motion of the screw. The flywheel consists of a flywheel bracket and a plurality of planetary gears. The flywheel bracket is fixed at one end of the screw. And the planetary gears are meshed with the ring gear. During the movement, the screw and the flywheel bracket are rotating synchronously. The planetary gears, driven by the flywheel bracket and the ring gear, are governed by a compound motion composed of revolution and rotation.Figure 2
Structure of a planetary flywheel inerter:(1) piston rod, (2) linear bearing, (3) mounting flange, (4) housing, (5) linear guide rail, (6) nut, (7) ball screw, (8) screw bearings, (9) flywheel bracket, (10) pin, (11) ring gear, (12) planetary gears, and (13) planetary gear bearings.


### 2.2. Prototype Production

The prototypes of flywheels and an inerter are shown in Figures3 and 4. The main structural parameters of the prototype are listed in Table 1.Table 1
Structural parameters of the prototype.
Name
Notation
Specification
Modulus of gear
M
m
o
d
2
Pressure angle of gear
α
20°
Teeth number of ring gear
Z
2
50
Reference radius of ring gear
R
0.050 m
Teeth number of planetary gear
Z
1
17
Reference radius of planetary gear
r
0.017 m
Thickness of planetary gear
l
0.015 m
Mass of planetary gear
m
p
0.094 kg
Rotational inertia of rotation of planetary gear
J
p
1.66 × 10−5 kg·m2
Rotational inertia of revolution of planetary gear
J
r
1.02 × 10−4 kg·m2
Number of planetary gears
N
4
Radius of pin
s
0.0095 m
Thickness of flywheel bracket
L
0.0075 m
Mass of flywheel bracket (including pins)
m
s
0.426 kg
Rotational inertia of flywheel bracket (including pins)
J
s
1
2.89 × 10−4 kg·m2
Major diameter of screw
D
0.020 m
Lead of screw
P
0.010 m
Rotational inertia of screw
J
s
2
2.80 × 10−5 kg·m2Figure 3
Prototypes of flywheel: (a) planetary flywheel and (b) traditional flywheel.
(a)
(b)Figure 4
Prototype of a planetary flywheel inerter.


## 2.1. Introduction of the Design and Principle of the Inerter

As exhibited in Figure2, the piston rod and the housing are the two moving terminals of the planetary flywheel inerter. The left end of the piston rod is threaded, and the housing has a mounting flange on it for connection. When the relative movement occurs between the two terminals of the inerter, the piston rod moves the nut along the linear bearing and the linear guide rail. The nut then converts the linear motion into the rotary motion of the screw. The flywheel consists of a flywheel bracket and a plurality of planetary gears. The flywheel bracket is fixed at one end of the screw. And the planetary gears are meshed with the ring gear. During the movement, the screw and the flywheel bracket are rotating synchronously. The planetary gears, driven by the flywheel bracket and the ring gear, are governed by a compound motion composed of revolution and rotation.Figure 2
Structure of a planetary flywheel inerter:(1) piston rod, (2) linear bearing, (3) mounting flange, (4) housing, (5) linear guide rail, (6) nut, (7) ball screw, (8) screw bearings, (9) flywheel bracket, (10) pin, (11) ring gear, (12) planetary gears, and (13) planetary gear bearings.


## 2.2. Prototype Production

The prototypes of flywheels and an inerter are shown in Figures3 and 4. The main structural parameters of the prototype are listed in Table 1.Table 1
Structural parameters of the prototype.
Name
Notation
Specification
Modulus of gear
M
m
o
d
2
Pressure angle of gear
α
20°
Teeth number of ring gear
Z
2
50
Reference radius of ring gear
R
0.050 m
Teeth number of planetary gear
Z
1
17
Reference radius of planetary gear
r
0.017 m
Thickness of planetary gear
l
0.015 m
Mass of planetary gear
m
p
0.094 kg
Rotational inertia of rotation of planetary gear
J
p
1.66 × 10−5 kg·m2
Rotational inertia of revolution of planetary gear
J
r
1.02 × 10−4 kg·m2
Number of planetary gears
N
4
Radius of pin
s
0.0095 m
Thickness of flywheel bracket
L
0.0075 m
Mass of flywheel bracket (including pins)
m
s
0.426 kg
Rotational inertia of flywheel bracket (including pins)
J
s
1
2.89 × 10−4 kg·m2
Major diameter of screw
D
0.020 m
Lead of screw
P
0.010 m
Rotational inertia of screw
J
s
2
2.80 × 10−5 kg·m2Figure 3
Prototypes of flywheel: (a) planetary flywheel and (b) traditional flywheel.
(a)
(b)Figure 4
Prototype of a planetary flywheel inerter.


## 3. Ideal Model and Advantages of a Planetary Flywheel Inerter

### 3.1. Ideal Model of a Planetary Flywheel Inerter

The inertanceb of an ball screw inerter could be calculated using the equation given in [26] as(2)b=2πP2Je,where Je is the flywheel’s rotational inertia.According to the law of conservation of energy, the torqueT of a planetary flywheel satisfies(3)Tθs=12Jsθ˙s2+N12Jpθ˙p2+12Jrθ˙s2,where Js=Js1+Js2 is the total rotational inertia of the flywheel bracket and the screw, θs is the angular displacement of the flywheel about the axis of the screw, and θp is the angular displacement of the planetary gear rotation.Differentiating both sides of (3) yields(4)Tθ˙s=Jsθ˙sθ¨s+NJpθ˙pθ¨p+Jrθ˙sθ¨.In addition, the planetary gears are rotating in accordance with the following relationship:(5)θ˙p=Z2-Z1Z1θ˙s,θ¨p=Z2-Z1Z1θ¨s.Substituting (5) into (4), we get(6)T=Js+NJpZ2-Z1Z12+Jrθ¨s.Then the planetary flywheel’s rotational inertia could be calculated by(7)Je=Tθ¨s=Js+NJpZ2-Z1Z12+Jr.Substituting (7) into (2) gives(8)b=2πP2Je=2πP2Js+NJpZ2-Z1Z12+Jr.Equation (8) reveals that the inertance b of a planetary flywheel inerter is directly proportional to the planetary flywheel’s rotational inertia Je. With the compound motion of revolution and rotation, the planetary gear’s rotational inertia of rotation Jp is added to the inerter’s rotational inertia with a gain factor Z2-Z1/Z12. The inerter’s rotational inertia is thus conspicuously increased from the traditional case where the mass is the same but is revolving about the axis of the screw only.By substituting the parameters in Table1 into (7) and (8), we obtained the rotational inertia to be 9.752 × 10−4 kg·m2 for the planetary flywheel prototype shown in Figure 3(a) and the inertance to be 385 kg for the inerter prototype depicted in Figure 4, where the mass of the planetary flywheel was mpf=ms+4mp=0.802 kg.


### 3.2. Advantages of a Planetary Flywheel

As the modulus of the gear is relatively small in terms of the radial scale of the flywheel group, the planetary gears, the flywheel bracket, and the pins are considered homogeneous cylinders. The radius of the planetary flywheel cylinder is approximately equal to the reference circle’s radius of the ring gear. The radius of a planetary gear cylinder is approximately equal to the radius of its reference circle. Given that the flywheel bracket and planetary gears are made of identical material with densityρ, the mass and the rotational inertia of them are(9)ms=πρLR2+πNρls2,mp=πρlr2-s2,(10)Js=πρ12LR4+Nls2R-r2,Jp=12πρlr2-s2r2,Jr=πρlr2-s2R-r2.Substituting (10) into (7) and normalizing the result, the rotational inertia of a planetary flywheel could be rewritten as(11)Je=12πρR4L1+Nτλ23-β21-λ2,where(12)β=sr,0<β<1,λ=Z1Z2=rR,0<λ<sin⁡π/N1+sin⁡π/NAdjacent  condition,τ=lL,τ>0.Then the planetary flywheel’s inertia-to-mass ratioA could be determined by(13)A=Jems+Nmp=R221+Nτλ23-β21-λ2-11+Nτλ2.For the traditional flywheel inerter sketched in Figure5, its inertia-to-mass ratio A′ could be calculated by(14)A′=Je′m′=1/2πρH′R′4πρH′R′2=R′22.Figure 5
Traditional flywheel inerter.Equation (14) indicates that the inertia-to-mass ratio of a traditional inerter is determined by the radial size of the flywheel R′ and is independent of the flywheel thickness. A larger inertia-to-mass ratio can be accomplished only by a larger flywheel radial dimension.According to (13) and (14), given that the two types of flywheels share the identical radial size, the ratio between the two flywheels’ inertia-to-mass ratios is(15)AA′=1+Nτλ23-β21-λ2-11+Nτλ2.It can be deduced from (15) that when λ<1-1/3-β2, there are (3-β2)(1-λ)2-1>0 and A>A′. That is, given the gear ratio λ of the planetary gear to the ring gear being less than a certain value, the planetary flywheel’s inertia amplification capacity is greater than that of the traditional flywheel.As a reference, a prototype of the traditional flywheel was produced as shown in the Figure3(b), which has a diameter of 100 mm, same as the diameter of the planetary flywheel prototype. The mass of the traditional flywheel was 0.974 kg, and the rotational inertia was 9.752 × 10−4 kg·m2. Given that the traditional flywheel is installed to the end of the screw, the inertance of the inerter was 385 kg, equal to that of the planetary flywheel inerter.Therefore, the mass of the planetary flywheel was 17.7% less than that of the traditional flywheel while keeping the flywheel diameters as well as the inerter inertances the same. The weight of the inerter was reduced significantly.


## 3.1. Ideal Model of a Planetary Flywheel Inerter

The inertanceb of an ball screw inerter could be calculated using the equation given in [26] as(2)b=2πP2Je,where Je is the flywheel’s rotational inertia.According to the law of conservation of energy, the torqueT of a planetary flywheel satisfies(3)Tθs=12Jsθ˙s2+N12Jpθ˙p2+12Jrθ˙s2,where Js=Js1+Js2 is the total rotational inertia of the flywheel bracket and the screw, θs is the angular displacement of the flywheel about the axis of the screw, and θp is the angular displacement of the planetary gear rotation.Differentiating both sides of (3) yields(4)Tθ˙s=Jsθ˙sθ¨s+NJpθ˙pθ¨p+Jrθ˙sθ¨.In addition, the planetary gears are rotating in accordance with the following relationship:(5)θ˙p=Z2-Z1Z1θ˙s,θ¨p=Z2-Z1Z1θ¨s.Substituting (5) into (4), we get(6)T=Js+NJpZ2-Z1Z12+Jrθ¨s.Then the planetary flywheel’s rotational inertia could be calculated by(7)Je=Tθ¨s=Js+NJpZ2-Z1Z12+Jr.Substituting (7) into (2) gives(8)b=2πP2Je=2πP2Js+NJpZ2-Z1Z12+Jr.Equation (8) reveals that the inertance b of a planetary flywheel inerter is directly proportional to the planetary flywheel’s rotational inertia Je. With the compound motion of revolution and rotation, the planetary gear’s rotational inertia of rotation Jp is added to the inerter’s rotational inertia with a gain factor Z2-Z1/Z12. The inerter’s rotational inertia is thus conspicuously increased from the traditional case where the mass is the same but is revolving about the axis of the screw only.By substituting the parameters in Table1 into (7) and (8), we obtained the rotational inertia to be 9.752 × 10−4 kg·m2 for the planetary flywheel prototype shown in Figure 3(a) and the inertance to be 385 kg for the inerter prototype depicted in Figure 4, where the mass of the planetary flywheel was mpf=ms+4mp=0.802 kg.


## 3.2. Advantages of a Planetary Flywheel

As the modulus of the gear is relatively small in terms of the radial scale of the flywheel group, the planetary gears, the flywheel bracket, and the pins are considered homogeneous cylinders. The radius of the planetary flywheel cylinder is approximately equal to the reference circle’s radius of the ring gear. The radius of a planetary gear cylinder is approximately equal to the radius of its reference circle. Given that the flywheel bracket and planetary gears are made of identical material with densityρ, the mass and the rotational inertia of them are(9)ms=πρLR2+πNρls2,mp=πρlr2-s2,(10)Js=πρ12LR4+Nls2R-r2,Jp=12πρlr2-s2r2,Jr=πρlr2-s2R-r2.Substituting (10) into (7) and normalizing the result, the rotational inertia of a planetary flywheel could be rewritten as(11)Je=12πρR4L1+Nτλ23-β21-λ2,where(12)β=sr,0<β<1,λ=Z1Z2=rR,0<λ<sin⁡π/N1+sin⁡π/NAdjacent  condition,τ=lL,τ>0.Then the planetary flywheel’s inertia-to-mass ratioA could be determined by(13)A=Jems+Nmp=R221+Nτλ23-β21-λ2-11+Nτλ2.For the traditional flywheel inerter sketched in Figure5, its inertia-to-mass ratio A′ could be calculated by(14)A′=Je′m′=1/2πρH′R′4πρH′R′2=R′22.Figure 5
Traditional flywheel inerter.Equation (14) indicates that the inertia-to-mass ratio of a traditional inerter is determined by the radial size of the flywheel R′ and is independent of the flywheel thickness. A larger inertia-to-mass ratio can be accomplished only by a larger flywheel radial dimension.According to (13) and (14), given that the two types of flywheels share the identical radial size, the ratio between the two flywheels’ inertia-to-mass ratios is(15)AA′=1+Nτλ23-β21-λ2-11+Nτλ2.It can be deduced from (15) that when λ<1-1/3-β2, there are (3-β2)(1-λ)2-1>0 and A>A′. That is, given the gear ratio λ of the planetary gear to the ring gear being less than a certain value, the planetary flywheel’s inertia amplification capacity is greater than that of the traditional flywheel.As a reference, a prototype of the traditional flywheel was produced as shown in the Figure3(b), which has a diameter of 100 mm, same as the diameter of the planetary flywheel prototype. The mass of the traditional flywheel was 0.974 kg, and the rotational inertia was 9.752 × 10−4 kg·m2. Given that the traditional flywheel is installed to the end of the screw, the inertance of the inerter was 385 kg, equal to that of the planetary flywheel inerter.Therefore, the mass of the planetary flywheel was 17.7% less than that of the traditional flywheel while keeping the flywheel diameters as well as the inerter inertances the same. The weight of the inerter was reduced significantly.


## 4. Dynamical Model of a Planetary Flywheel Inerter

### 4.1. Dynamical Model of a Planetary Flywheel

Given that all the gears are standard cylindrical spur gears and the planetary gears have the same physical and geometric parameters, the concentrated mass method given in [27] was used to establish a planetary flywheel dynamical model, which is shown in Figure 6.Figure 6
Schematic dynamical model of planetary flywheel.Whererc=R-r is the center distance between the planetary gear and the ring gear, i=1,2,…,N, kri is the meshing stiffness of the meshing gear pair along the line of action, cri is the meshing damping coefficient of the gear pair, gri indicates the gear backslash, Fpi is the meshing force between the planetary gear and the ring gear, and Tfpi is the equivalent friction torque of the planetary gear about its rotation axis.Lagrange’s equation was used to derive the differential equation governing the dynamical planetary flywheel system as follows:(16)Jpθ¨pi=Fpirb-Tfpi,Js+NJrθ¨s=Ta-∑i=1NFpirccos⁡α,where rb is the radius of the planetary gear’s base circle and Ta is the total torque of the planetary flywheel.The generalized coordinates and the equivalent mass of gears could be defined as(17)xpi=θpirb,xs=θsrccos⁡α,mpe=Jprb2,mre=Jrrc2cos2α,mse=Jsrc2cos2α.Substituting (17) into (16), we get(18)mpex¨pi=Fpi-Tfpirb,mse+Nmrex¨s=Tarccos⁡α-∑i=1NFpi.The elastic meshing forcePi and the viscous meshing force Di of the gear pair could be expressed as(19)Pi=krixs-xpi,Di=crix˙s-x˙pi.The meshing forceFpi could be written using the mechanical clearance model given in [28] as(20)Fpi=krixs-xpi-gri+DiPi+Di≥krigri0Pi+Di<krigrikrixs-xpi+gri+DiPi+Di≤-krigri.Using the Coulomb viscous friction model, the equivalent friction torque of the planetary wheel about the rotation axisTfpi is(21)Tfpi=Tfpmsgnθ˙pi+Tfpvθ˙pi,where Tfpm is the Coulomb friction torque of the planetary wheel about the rotation axis; Tfpv is the viscous friction torque coefficient of the planetary wheel about the rotation axis.According to (18)–(21), a dynamical model of the planetary flywheel is shown in Figure 7, where Ffpi=Tfpi/rb. The torque of the planetary flywheel is generated by the revolution and rotation of the gears.Figure 7
Dynamical model of a planetary flywheel.


### 4.2. Dynamical Model of a Ball Screw Type Inerter

A traditional flywheel ball screw inerter model with nonlinear factors such as clearance and friction is presented in [3]. Given the traditional flywheel being replaced by a planetary flywheel, the mechanical principle of the planetary flywheel inerter is shown in Figure 8.Figure 8
Mechanical scheme of a ball screw type inerter with planetary flywheel.Whereks is the axial stiffness of the nut-screw pair, cs is the axial damping coefficient of the nut-screw pair, ε1 is the axial backslash between the nut and the screw, ε2 is the axial internal clearance of the bearing, mr is the mass of the nut and the piston rod, x3 is the axial displacement of the screw, Fd is the axial driving force applied to the screw by the nut, Ta is the torque transmitted to the flywheel by the screw, Tf1 is the friction torque between the nut and the screw, and Tf2 is the bearing friction torque.The output forceF of the inerter and the forces applied on the piston rod and the nut satisfy(22)F=F1-F2=Fd+mrx¨1.The axial elastic forcePi and the axial viscous force Di applied upon the nut-screw system could be given as(23)Ps=ksx1-x3,Ds=csx˙1-x˙3.Letε=ε1+ε2; Fd could be written as(24)Fd=ksx1-x3-ε+DsPs+Ds≥ksε0Ps+Ds<ksεksx1-x3+ε+DsPs+Ds≥-ksε.The torque applied to the flywheel by the screw is(25)Ta=Td-Tf,where Td is the driving torque generated by the nut-screw pair and Tf=Tf1+Tf2 is the total friction torque of the nut-screw pair and the bearing. Using the Stribeck friction model given in [29], the friction torque Tf can be expressed as(26)Tf=Tc+Ts-Tce-θ˙s/θ˙stribecksgnθ˙s+Tvθ˙s,where Tc is the total Coulomb friction torque of the nut-screw pair and the bearing; Ts is the total static friction torque; θ˙stribeck is the Stribeck velocity; Tv is the viscous friction torque coefficient.By the principle of nut-screw pair, the relationship between the axial forces and their corresponding torques is(27)Fd=2πPTd,Ff=2πPTf,Fa=2πPTa.Based on the dynamical model of a planetary flywheel established in Section4.1 and (22)~(27), a comprehensive dynamical model of a planetary flywheel inerter is presented in Figure 9, where θs=2π/P(x3-x2).Figure 9
Dynamical model of a planetary flywheel inerter.


## 4.1. Dynamical Model of a Planetary Flywheel

Given that all the gears are standard cylindrical spur gears and the planetary gears have the same physical and geometric parameters, the concentrated mass method given in [27] was used to establish a planetary flywheel dynamical model, which is shown in Figure 6.Figure 6
Schematic dynamical model of planetary flywheel.Whererc=R-r is the center distance between the planetary gear and the ring gear, i=1,2,…,N, kri is the meshing stiffness of the meshing gear pair along the line of action, cri is the meshing damping coefficient of the gear pair, gri indicates the gear backslash, Fpi is the meshing force between the planetary gear and the ring gear, and Tfpi is the equivalent friction torque of the planetary gear about its rotation axis.Lagrange’s equation was used to derive the differential equation governing the dynamical planetary flywheel system as follows:(16)Jpθ¨pi=Fpirb-Tfpi,Js+NJrθ¨s=Ta-∑i=1NFpirccos⁡α,where rb is the radius of the planetary gear’s base circle and Ta is the total torque of the planetary flywheel.The generalized coordinates and the equivalent mass of gears could be defined as(17)xpi=θpirb,xs=θsrccos⁡α,mpe=Jprb2,mre=Jrrc2cos2α,mse=Jsrc2cos2α.Substituting (17) into (16), we get(18)mpex¨pi=Fpi-Tfpirb,mse+Nmrex¨s=Tarccos⁡α-∑i=1NFpi.The elastic meshing forcePi and the viscous meshing force Di of the gear pair could be expressed as(19)Pi=krixs-xpi,Di=crix˙s-x˙pi.The meshing forceFpi could be written using the mechanical clearance model given in [28] as(20)Fpi=krixs-xpi-gri+DiPi+Di≥krigri0Pi+Di<krigrikrixs-xpi+gri+DiPi+Di≤-krigri.Using the Coulomb viscous friction model, the equivalent friction torque of the planetary wheel about the rotation axisTfpi is(21)Tfpi=Tfpmsgnθ˙pi+Tfpvθ˙pi,where Tfpm is the Coulomb friction torque of the planetary wheel about the rotation axis; Tfpv is the viscous friction torque coefficient of the planetary wheel about the rotation axis.According to (18)–(21), a dynamical model of the planetary flywheel is shown in Figure 7, where Ffpi=Tfpi/rb. The torque of the planetary flywheel is generated by the revolution and rotation of the gears.Figure 7
Dynamical model of a planetary flywheel.


## 4.2. Dynamical Model of a Ball Screw Type Inerter

A traditional flywheel ball screw inerter model with nonlinear factors such as clearance and friction is presented in [3]. Given the traditional flywheel being replaced by a planetary flywheel, the mechanical principle of the planetary flywheel inerter is shown in Figure 8.Figure 8
Mechanical scheme of a ball screw type inerter with planetary flywheel.Whereks is the axial stiffness of the nut-screw pair, cs is the axial damping coefficient of the nut-screw pair, ε1 is the axial backslash between the nut and the screw, ε2 is the axial internal clearance of the bearing, mr is the mass of the nut and the piston rod, x3 is the axial displacement of the screw, Fd is the axial driving force applied to the screw by the nut, Ta is the torque transmitted to the flywheel by the screw, Tf1 is the friction torque between the nut and the screw, and Tf2 is the bearing friction torque.The output forceF of the inerter and the forces applied on the piston rod and the nut satisfy(22)F=F1-F2=Fd+mrx¨1.The axial elastic forcePi and the axial viscous force Di applied upon the nut-screw system could be given as(23)Ps=ksx1-x3,Ds=csx˙1-x˙3.Letε=ε1+ε2; Fd could be written as(24)Fd=ksx1-x3-ε+DsPs+Ds≥ksε0Ps+Ds<ksεksx1-x3+ε+DsPs+Ds≥-ksε.The torque applied to the flywheel by the screw is(25)Ta=Td-Tf,where Td is the driving torque generated by the nut-screw pair and Tf=Tf1+Tf2 is the total friction torque of the nut-screw pair and the bearing. Using the Stribeck friction model given in [29], the friction torque Tf can be expressed as(26)Tf=Tc+Ts-Tce-θ˙s/θ˙stribecksgnθ˙s+Tvθ˙s,where Tc is the total Coulomb friction torque of the nut-screw pair and the bearing; Ts is the total static friction torque; θ˙stribeck is the Stribeck velocity; Tv is the viscous friction torque coefficient.By the principle of nut-screw pair, the relationship between the axial forces and their corresponding torques is(27)Fd=2πPTd,Ff=2πPTf,Fa=2πPTa.Based on the dynamical model of a planetary flywheel established in Section4.1 and (22)~(27), a comprehensive dynamical model of a planetary flywheel inerter is presented in Figure 9, where θs=2π/P(x3-x2).Figure 9
Dynamical model of a planetary flywheel inerter.


## 5. Experiment

### 5.1. Test-Rig Setup

As shown in Figures10 and 11, the inerter was mounted on the test-rig frame through the mounting flange. The hydraulic actuator’s output rod was fixed to the piston rod through the connector. The hydraulic actuator was controlled by a Moog™ portable controller to create a relative movement between the flange end and the piston rod end of the inerter. The NI Compact-RIO™ data acquisition system was deployed to collect the displacement and force data of the inerter with sampling frequency of 1000 Hz and send it to the PC system for subsequent data analysis.Figure 10
Scheme of the test-rig.Figure 11
Prototype of the inerter and the test-rig:(1) hydraulic actuator, (2) force transducer, (3) connector, (4) test-rig frame, (5) inerter, and (6) flywheel chamber.


### 5.2. Parameter Identification and Model Validation

In order to reduce the computational complexity of the parameter identification process, the parameter identification was carried out stepwise. Firstly, no flywheel was attached to the inerter. The frictional parameters of the inerter were identified. Next, a traditional flywheel was mounted on the inerter and the dynamical parameters of the ball screw were identified. Then, the traditional flywheel was replaced with the planetary flywheel. And the dynamical parameters of the planetary flywheel were identified. Finally, the determined parameters were substituted into the model, and the accuracy of the model was validated by various experiments. The model was solved using the MATLAB™ software with the ode45 solver for both the parameter identification and the model validation.


#### 5.2.1. Frictional Parameters of the Ball Screw Type Inerter

During the experiment phase for the frictional parameter identification, no flywheel was fixed to the inerter. The hydraulic actuator drove the piston rod in a sinusoidal style with amplitude of 20 mm and a frequency of 0.1 Hz. At this time, the frequency of motion was very low, the relative velocity and acceleration of one terminal to the other of the inerter were small, and the damping and inertia effect could be ignored. The output force of the inerter was basically due to friction.In the experiment, the output force data of the inerter prototypeFexp was collected. A population of candidate frictional parameter solutions was given to the Stribeck friction model described by (26) and (27). The resulting model output forces Fmod=Ff were then compared with the target Fexp using the following loss function:(28)min⁡Fexpt-FmodTc,Ts,θ˙s,Tv22.F
e
x
p and Fmod were both filtered by a 6-order low-pass Butterworth filter with a cutoff frequency of 0.5 Hz. As a result, the Stribeck friction model parameters of the overall ball screw type inerter were estimated and listed in Table 2.Table 2
Frictional parameters of the ball screw type inerter.
Parameter
Specification
T
c
1.975 × 10−1 N·m
T
s
2.231 × 10−1 N·m
θ
˙
s
t
r
i
b
e
c
k
5.437 × 108 rad/s
T
v
2.703 × 10−3 N·m·s/radAs shown in Figure12, given the identified inerter’s frictional parameters listed in Table 2, the model could basically describe the friction characteristics of the ball screw type inerter.Figure 12
Low-frequency response of the ball screw type inerter: displacement input amplitude = 20 mm and frequency = 0.1 Hz.


#### 5.2.2. Parameters of the Ball Screw

The traditional flywheel shown in Figure3(b) was installed at the end of the screw. The purpose of the flywheel was to amplify the magnitude of output for improving the measurement accuracy. The hydraulic actuator drove the piston rod in a sinusoidal style with amplitude of 10 mm and a frequency of 3 Hz.With the inerter’s frictional parameters given in Table2, a population of candidate dynamical parameter solutions of the ball screw were set to the model of the traditional flywheel inerter described in the literature [3]. The resulting model output forces Fmod were then compared with the target Fexp using the following loss function:(29)min⁡Fexpt-Fmodmr,b,ks,cs,ε22.As a result, the parameters of the ball screw were estimated and listed in Table 3.Table 3
Parameters of the ball screw.
Parameter
Specification
m
r
1.701 kg
k
s
1.168 × 107 N/m
c
s
1.223 × 103 N·s/m
ε
1.006 × 10−7 mAs displayed in Figure13, given that the identified inerter’s frictional parameters and the ball screw’s dynamical ones are shown in Table 3, the model output force agreed well with the experiment data.Figure 13
Time response of the traditional flywheel inerter: displacement input amplitude = 10 mm and frequency = 3 Hz.


#### 5.2.3. Parameters of the Planetary Flywheel

The traditional flywheel was then replaced with the planetary flywheel shown in Figure3(a). The hydraulic actuator drove the piston rod in a sinusoidal style with amplitude of 10 mm and a frequency of 3 Hz.Since the four planetary gears were manufactured together, their mutual difference was considered ignorable. It was assumed that their meshing stiffness, meshing damping coefficient, backslash, and friction torque were the same, respectively.With the inerter’s frictional parameters given in Table2 and the parameters of the ball screw given in Table 3, a population of candidate planetary flywheel parameter solutions were assigned to the model of planetary flywheel inerter shown in Figure 9. The resulting model output forces Fmod=F1-F2 were then compared with the target Fexp using the following loss function:(30)min⁡Fexpt-Fmodkr,cr,gr,Tpm,Tpv22.As a result, the parameters of the planetary flywheel were estimated and listed in Table 4.Table 4
Parameters of the planetary flywheel.
Parameter
Specification
k
r
1.353 × 109 N/m
c
r
4.671 × 104 N·s/m
g
r
1.838 × 10−6 mm
T
p
m
8.067 × 10−4 N·m
T
p
v
1.366 × 10−4 N·m·s/radAs illustrated in Figure14(b), the model output force agreed well with the experiment data as all the identified parameters listed in Table 4 are given.Figure 14
Time response of the inerter with planetary flywheel: (a) displacement input amplitude = 20 mm and frequency = 0.1 Hz; (b) amplitude = 10 mm and frequency = 3 Hz; (c) amplitude = 7 mm and frequency = 5 Hz; (d) amplitude = 1.5 mm and frequency = 10 Hz.
(a)
(b)
(c)
(d)


#### 5.2.4. Validation of the Dynamical Model

With all the key parameters of the planetary flywheel inerter determined, the frictional parameters, the parameters of the ball screw, and the planetary flywheel parameters were substituted into the inerter model shown in Figure9. Then, by inputting the sinusoidal displacement excitation of different frequencies to the inerter, the feasibility of the model was tested by comparing the experiment data with the model output.Figures14(a), 14(c), and 14(d) demonstrate that the output of the model agreed well with the output of the experimental inerter under the excitation of three frequencies (low frequency: frequency = 0.1 Hz and amplitude = 20 mm; Fexp and Fmod were both filtered by a 6-order low-pass Butterworth filter with a cutoff frequency of 0.5 Hz; intermediate frequency: frequency = 5 Hz and amplitude = 7 mm; high frequency: frequency = 10 Hz and amplitude = 1.5 mm). Experiments showed an outstanding accuracy of the theoretical dynamical model.


### 5.3. Random Input Test and Transfer Characteristics Analysis

To reveal the transfer characteristics of the inerter in the frequency domain, a white noise speed excitation of a range of 0.5~15 Hz as shown in Figure15 was applied to the planetary flywheel inerter and its traditional flywheel counterpart by the hydraulic actuator. The power spectral density of input was(31)Ginputf=6.4×10-6·f-2m2/Hz.Figure 15
Random input: (a) input signal and (b) power spectrum density of input signal.
(a)
(b)The transfer characteristics of an inerter is defined as(32)Hf=Gin-outfGinputf,where Gin-outf is the cross power spectral density of the input displacement and the output force.From Figure16, the transfer characteristics of the planetary flywheel inerter are similar to that of the traditional flywheel inerter with the same inertance. In the lower frequency band, the amplitude gain of the planetary flywheel inerter is slightly larger than that of the traditional flywheel counterpart, as a result of the introduction of the planetary gears’ friction to the former one. In the higher frequency band above 2 Hz, the transfer characteristics of the two types of inerter are basically the same, and both increase with a higher frequency, in line with the characteristics of an ideal inerter.Figure 16
Frequency response of the inerters.Moreover, the transfer characteristics of the dynamical model were in good agreement with the experimental results at a relatively high frequency, which further validates the accuracy of the model. When the frequency is low, the amplitude of the model is smaller than that of the experiment. This is mainly because the impact of friction on the output force is relatively large at low frequency, but the friction model is only a theoretical approximation to the friction phenomenon. In general, the change of model output with frequency was consistent with that of the experiment data, which proves the model’s validity of being used in the investigation and design work on mechanical network.Literature [18] gives the response of the ball screw inerter and the rack-pinion inerter with traditional flywheel. When excited at a frequency of 0.1 Hz, the output force of the rack-pinion inerter is a triangular wave with large fluctuating, and the output force of the ball screw inerter is similar to a square wave. Although the ball screw inerter with planetary flywheel and the rack-pinion inerter both use gear mechanisms, their dynamic characteristics are quite different. This is because, for the rack-pinion inerter, the conversion of the linear motion to the rotational motion is accomplished by the rack and pinion. The gear backlash will cause the output force of the inerter to have a greater fluctuation. For the planetary flywheel inerter, the conversion of the linear motion to the rotational motion is accomplished by the ball screw which can be preloaded and the simultaneous engagement of multiple planetary gears produces an average effect on the fluctuation of the output force. The gear backlash has little effect on the output force of the planetary flywheel inerter. The dynamic performance of the ball screw inerter with planetary flywheel is similar to that of the ball screw inerter with traditional flywheel. In addition, the ball screw inerter has a greater inertia-to-mass ratio than the rack-pinion inerter. The planetary flywheel inerter can reduce the flywheel weight more effectively. The rack-pinion inerter is more suitable for the applications where the flywheel cannot be mounted on the inerter axis due to space constraints.As presented in Figures13–16, the inerter’s dynamic characteristics of the planetary flywheel type and the traditional flywheel type are basically the same, indicating that the planetary flywheel can be the core inertial component of an inerter and be a replacement of the traditional flywheel.


## 5.1. Test-Rig Setup

As shown in Figures10 and 11, the inerter was mounted on the test-rig frame through the mounting flange. The hydraulic actuator’s output rod was fixed to the piston rod through the connector. The hydraulic actuator was controlled by a Moog™ portable controller to create a relative movement between the flange end and the piston rod end of the inerter. The NI Compact-RIO™ data acquisition system was deployed to collect the displacement and force data of the inerter with sampling frequency of 1000 Hz and send it to the PC system for subsequent data analysis.Figure 10
Scheme of the test-rig.Figure 11
Prototype of the inerter and the test-rig:(1) hydraulic actuator, (2) force transducer, (3) connector, (4) test-rig frame, (5) inerter, and (6) flywheel chamber.


## 5.2. Parameter Identification and Model Validation

In order to reduce the computational complexity of the parameter identification process, the parameter identification was carried out stepwise. Firstly, no flywheel was attached to the inerter. The frictional parameters of the inerter were identified. Next, a traditional flywheel was mounted on the inerter and the dynamical parameters of the ball screw were identified. Then, the traditional flywheel was replaced with the planetary flywheel. And the dynamical parameters of the planetary flywheel were identified. Finally, the determined parameters were substituted into the model, and the accuracy of the model was validated by various experiments. The model was solved using the MATLAB™ software with the ode45 solver for both the parameter identification and the model validation.


### 5.2.1. Frictional Parameters of the Ball Screw Type Inerter

During the experiment phase for the frictional parameter identification, no flywheel was fixed to the inerter. The hydraulic actuator drove the piston rod in a sinusoidal style with amplitude of 20 mm and a frequency of 0.1 Hz. At this time, the frequency of motion was very low, the relative velocity and acceleration of one terminal to the other of the inerter were small, and the damping and inertia effect could be ignored. The output force of the inerter was basically due to friction.In the experiment, the output force data of the inerter prototypeFexp was collected. A population of candidate frictional parameter solutions was given to the Stribeck friction model described by (26) and (27). The resulting model output forces Fmod=Ff were then compared with the target Fexp using the following loss function:(28)min⁡Fexpt-FmodTc,Ts,θ˙s,Tv22.F
e
x
p and Fmod were both filtered by a 6-order low-pass Butterworth filter with a cutoff frequency of 0.5 Hz. As a result, the Stribeck friction model parameters of the overall ball screw type inerter were estimated and listed in Table 2.Table 2
Frictional parameters of the ball screw type inerter.
Parameter
Specification
T
c
1.975 × 10−1 N·m
T
s
2.231 × 10−1 N·m
θ
˙
s
t
r
i
b
e
c
k
5.437 × 108 rad/s
T
v
2.703 × 10−3 N·m·s/radAs shown in Figure12, given the identified inerter’s frictional parameters listed in Table 2, the model could basically describe the friction characteristics of the ball screw type inerter.Figure 12
Low-frequency response of the ball screw type inerter: displacement input amplitude = 20 mm and frequency = 0.1 Hz.


### 5.2.2. Parameters of the Ball Screw

The traditional flywheel shown in Figure3(b) was installed at the end of the screw. The purpose of the flywheel was to amplify the magnitude of output for improving the measurement accuracy. The hydraulic actuator drove the piston rod in a sinusoidal style with amplitude of 10 mm and a frequency of 3 Hz.With the inerter’s frictional parameters given in Table2, a population of candidate dynamical parameter solutions of the ball screw were set to the model of the traditional flywheel inerter described in the literature [3]. The resulting model output forces Fmod were then compared with the target Fexp using the following loss function:(29)min⁡Fexpt-Fmodmr,b,ks,cs,ε22.As a result, the parameters of the ball screw were estimated and listed in Table 3.Table 3
Parameters of the ball screw.
Parameter
Specification
m
r
1.701 kg
k
s
1.168 × 107 N/m
c
s
1.223 × 103 N·s/m
ε
1.006 × 10−7 mAs displayed in Figure13, given that the identified inerter’s frictional parameters and the ball screw’s dynamical ones are shown in Table 3, the model output force agreed well with the experiment data.Figure 13
Time response of the traditional flywheel inerter: displacement input amplitude = 10 mm and frequency = 3 Hz.


### 5.2.3. Parameters of the Planetary Flywheel

The traditional flywheel was then replaced with the planetary flywheel shown in Figure3(a). The hydraulic actuator drove the piston rod in a sinusoidal style with amplitude of 10 mm and a frequency of 3 Hz.Since the four planetary gears were manufactured together, their mutual difference was considered ignorable. It was assumed that their meshing stiffness, meshing damping coefficient, backslash, and friction torque were the same, respectively.With the inerter’s frictional parameters given in Table2 and the parameters of the ball screw given in Table 3, a population of candidate planetary flywheel parameter solutions were assigned to the model of planetary flywheel inerter shown in Figure 9. The resulting model output forces Fmod=F1-F2 were then compared with the target Fexp using the following loss function:(30)min⁡Fexpt-Fmodkr,cr,gr,Tpm,Tpv22.As a result, the parameters of the planetary flywheel were estimated and listed in Table 4.Table 4
Parameters of the planetary flywheel.
Parameter
Specification
k
r
1.353 × 109 N/m
c
r
4.671 × 104 N·s/m
g
r
1.838 × 10−6 mm
T
p
m
8.067 × 10−4 N·m
T
p
v
1.366 × 10−4 N·m·s/radAs illustrated in Figure14(b), the model output force agreed well with the experiment data as all the identified parameters listed in Table 4 are given.Figure 14
Time response of the inerter with planetary flywheel: (a) displacement input amplitude = 20 mm and frequency = 0.1 Hz; (b) amplitude = 10 mm and frequency = 3 Hz; (c) amplitude = 7 mm and frequency = 5 Hz; (d) amplitude = 1.5 mm and frequency = 10 Hz.
(a)
(b)
(c)
(d)


### 5.2.4. Validation of the Dynamical Model

With all the key parameters of the planetary flywheel inerter determined, the frictional parameters, the parameters of the ball screw, and the planetary flywheel parameters were substituted into the inerter model shown in Figure9. Then, by inputting the sinusoidal displacement excitation of different frequencies to the inerter, the feasibility of the model was tested by comparing the experiment data with the model output.Figures14(a), 14(c), and 14(d) demonstrate that the output of the model agreed well with the output of the experimental inerter under the excitation of three frequencies (low frequency: frequency = 0.1 Hz and amplitude = 20 mm; Fexp and Fmod were both filtered by a 6-order low-pass Butterworth filter with a cutoff frequency of 0.5 Hz; intermediate frequency: frequency = 5 Hz and amplitude = 7 mm; high frequency: frequency = 10 Hz and amplitude = 1.5 mm). Experiments showed an outstanding accuracy of the theoretical dynamical model.


## 5.2.1. Frictional Parameters of the Ball Screw Type Inerter

During the experiment phase for the frictional parameter identification, no flywheel was fixed to the inerter. The hydraulic actuator drove the piston rod in a sinusoidal style with amplitude of 20 mm and a frequency of 0.1 Hz. At this time, the frequency of motion was very low, the relative velocity and acceleration of one terminal to the other of the inerter were small, and the damping and inertia effect could be ignored. The output force of the inerter was basically due to friction.In the experiment, the output force data of the inerter prototypeFexp was collected. A population of candidate frictional parameter solutions was given to the Stribeck friction model described by (26) and (27). The resulting model output forces Fmod=Ff were then compared with the target Fexp using the following loss function:(28)min⁡Fexpt-FmodTc,Ts,θ˙s,Tv22.F
e
x
p and Fmod were both filtered by a 6-order low-pass Butterworth filter with a cutoff frequency of 0.5 Hz. As a result, the Stribeck friction model parameters of the overall ball screw type inerter were estimated and listed in Table 2.Table 2
Frictional parameters of the ball screw type inerter.
Parameter
Specification
T
c
1.975 × 10−1 N·m
T
s
2.231 × 10−1 N·m
θ
˙
s
t
r
i
b
e
c
k
5.437 × 108 rad/s
T
v
2.703 × 10−3 N·m·s/radAs shown in Figure12, given the identified inerter’s frictional parameters listed in Table 2, the model could basically describe the friction characteristics of the ball screw type inerter.Figure 12
Low-frequency response of the ball screw type inerter: displacement input amplitude = 20 mm and frequency = 0.1 Hz.


## 5.2.2. Parameters of the Ball Screw

The traditional flywheel shown in Figure3(b) was installed at the end of the screw. The purpose of the flywheel was to amplify the magnitude of output for improving the measurement accuracy. The hydraulic actuator drove the piston rod in a sinusoidal style with amplitude of 10 mm and a frequency of 3 Hz.With the inerter’s frictional parameters given in Table2, a population of candidate dynamical parameter solutions of the ball screw were set to the model of the traditional flywheel inerter described in the literature [3]. The resulting model output forces Fmod were then compared with the target Fexp using the following loss function:(29)min⁡Fexpt-Fmodmr,b,ks,cs,ε22.As a result, the parameters of the ball screw were estimated and listed in Table 3.Table 3
Parameters of the ball screw.
Parameter
Specification
m
r
1.701 kg
k
s
1.168 × 107 N/m
c
s
1.223 × 103 N·s/m
ε
1.006 × 10−7 mAs displayed in Figure13, given that the identified inerter’s frictional parameters and the ball screw’s dynamical ones are shown in Table 3, the model output force agreed well with the experiment data.Figure 13
Time response of the traditional flywheel inerter: displacement input amplitude = 10 mm and frequency = 3 Hz.


## 5.2.3. Parameters of the Planetary Flywheel

The traditional flywheel was then replaced with the planetary flywheel shown in Figure3(a). The hydraulic actuator drove the piston rod in a sinusoidal style with amplitude of 10 mm and a frequency of 3 Hz.Since the four planetary gears were manufactured together, their mutual difference was considered ignorable. It was assumed that their meshing stiffness, meshing damping coefficient, backslash, and friction torque were the same, respectively.With the inerter’s frictional parameters given in Table2 and the parameters of the ball screw given in Table 3, a population of candidate planetary flywheel parameter solutions were assigned to the model of planetary flywheel inerter shown in Figure 9. The resulting model output forces Fmod=F1-F2 were then compared with the target Fexp using the following loss function:(30)min⁡Fexpt-Fmodkr,cr,gr,Tpm,Tpv22.As a result, the parameters of the planetary flywheel were estimated and listed in Table 4.Table 4
Parameters of the planetary flywheel.
Parameter
Specification
k
r
1.353 × 109 N/m
c
r
4.671 × 104 N·s/m
g
r
1.838 × 10−6 mm
T
p
m
8.067 × 10−4 N·m
T
p
v
1.366 × 10−4 N·m·s/radAs illustrated in Figure14(b), the model output force agreed well with the experiment data as all the identified parameters listed in Table 4 are given.Figure 14
Time response of the inerter with planetary flywheel: (a) displacement input amplitude = 20 mm and frequency = 0.1 Hz; (b) amplitude = 10 mm and frequency = 3 Hz; (c) amplitude = 7 mm and frequency = 5 Hz; (d) amplitude = 1.5 mm and frequency = 10 Hz.
(a)
(b)
(c)
(d)


## 5.2.4. Validation of the Dynamical Model

With all the key parameters of the planetary flywheel inerter determined, the frictional parameters, the parameters of the ball screw, and the planetary flywheel parameters were substituted into the inerter model shown in Figure9. Then, by inputting the sinusoidal displacement excitation of different frequencies to the inerter, the feasibility of the model was tested by comparing the experiment data with the model output.Figures14(a), 14(c), and 14(d) demonstrate that the output of the model agreed well with the output of the experimental inerter under the excitation of three frequencies (low frequency: frequency = 0.1 Hz and amplitude = 20 mm; Fexp and Fmod were both filtered by a 6-order low-pass Butterworth filter with a cutoff frequency of 0.5 Hz; intermediate frequency: frequency = 5 Hz and amplitude = 7 mm; high frequency: frequency = 10 Hz and amplitude = 1.5 mm). Experiments showed an outstanding accuracy of the theoretical dynamical model.


## 5.3. Random Input Test and Transfer Characteristics Analysis

To reveal the transfer characteristics of the inerter in the frequency domain, a white noise speed excitation of a range of 0.5~15 Hz as shown in Figure15 was applied to the planetary flywheel inerter and its traditional flywheel counterpart by the hydraulic actuator. The power spectral density of input was(31)Ginputf=6.4×10-6·f-2m2/Hz.Figure 15
Random input: (a) input signal and (b) power spectrum density of input signal.
(a)
(b)The transfer characteristics of an inerter is defined as(32)Hf=Gin-outfGinputf,where Gin-outf is the cross power spectral density of the input displacement and the output force.From Figure16, the transfer characteristics of the planetary flywheel inerter are similar to that of the traditional flywheel inerter with the same inertance. In the lower frequency band, the amplitude gain of the planetary flywheel inerter is slightly larger than that of the traditional flywheel counterpart, as a result of the introduction of the planetary gears’ friction to the former one. In the higher frequency band above 2 Hz, the transfer characteristics of the two types of inerter are basically the same, and both increase with a higher frequency, in line with the characteristics of an ideal inerter.Figure 16
Frequency response of the inerters.Moreover, the transfer characteristics of the dynamical model were in good agreement with the experimental results at a relatively high frequency, which further validates the accuracy of the model. When the frequency is low, the amplitude of the model is smaller than that of the experiment. This is mainly because the impact of friction on the output force is relatively large at low frequency, but the friction model is only a theoretical approximation to the friction phenomenon. In general, the change of model output with frequency was consistent with that of the experiment data, which proves the model’s validity of being used in the investigation and design work on mechanical network.Literature [18] gives the response of the ball screw inerter and the rack-pinion inerter with traditional flywheel. When excited at a frequency of 0.1 Hz, the output force of the rack-pinion inerter is a triangular wave with large fluctuating, and the output force of the ball screw inerter is similar to a square wave. Although the ball screw inerter with planetary flywheel and the rack-pinion inerter both use gear mechanisms, their dynamic characteristics are quite different. This is because, for the rack-pinion inerter, the conversion of the linear motion to the rotational motion is accomplished by the rack and pinion. The gear backlash will cause the output force of the inerter to have a greater fluctuation. For the planetary flywheel inerter, the conversion of the linear motion to the rotational motion is accomplished by the ball screw which can be preloaded and the simultaneous engagement of multiple planetary gears produces an average effect on the fluctuation of the output force. The gear backlash has little effect on the output force of the planetary flywheel inerter. The dynamic performance of the ball screw inerter with planetary flywheel is similar to that of the ball screw inerter with traditional flywheel. In addition, the ball screw inerter has a greater inertia-to-mass ratio than the rack-pinion inerter. The planetary flywheel inerter can reduce the flywheel weight more effectively. The rack-pinion inerter is more suitable for the applications where the flywheel cannot be mounted on the inerter axis due to space constraints.As presented in Figures13–16, the inerter’s dynamic characteristics of the planetary flywheel type and the traditional flywheel type are basically the same, indicating that the planetary flywheel can be the core inertial component of an inerter and be a replacement of the traditional flywheel.


## 6. Conclusion

(
1
) In this paper, we propose the planetary flywheel inerter as a new structural type of inerters. When a relative movement occurs between the two terminals of the inerter, the planetary flywheel produces a compound screw-driven motion composed of revolving and rotating. The torque of the flywheel is proportional to the rotational acceleration of the screw. Therefore, the output force is proportional to the relative linear acceleration between the inerter’s two terminals. In other words, the planetary flywheel inerter is compatible with the mechanical output requirements of an ideal inerter.(
2
) A nonlinear dynamical model of the planetary flywheel inerter is established by theoretical analysis. Based on the model parameter identification, the response characteristics of the model shows remarkable agreement with the experimental data, which proves the model’s validity of being used in the investigation and design work on mechanical network.(
3
) The dynamical transfer characteristics of the planetary flywheel inerter and that of the traditional flywheel inerter are basically the same. However, the former’s transmissibility is slightly larger than that of the latter in low-frequency range due to the friction torque of the planetary gears.(
4
) The planetary flywheel can be a good substitute for the traditional flywheel. While the same inertance is accomplished, the mass of a planetary flywheel is notably less than a traditional flywheel with the same outer diameter, which helps to make the inerter lightweight.



---

*Source: 2631539-2018-01-28.xml*