Next Article in Journal
A Review on Rail Defect Detection Systems Based on Wireless Sensors
Previous Article in Journal
Mapping Local Climate Zones in the Urban Environment: The Optimal Combination of Data Source and Classifier
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Vibration Analysis of a 1-DOF System Coupled with a Nonlinear Energy Sink with a Fractional Order Inerter

1
College of Mechanical and Electronic Engineering, Nanjing Forestry University, Nanjing 210037, China
2
College of Intelligent Equipment Engineering, Wuxi Taihu University, Wuxi 214151, China
3
College of Automobile and Traffic Engineering, Nanjing Forestry University, Nanjing 210037, China
4
Anhui Vocational and Technical College, College of Mechanical Engineering, Hefei 230009, China
*
Author to whom correspondence should be addressed.
Sensors 2022, 22(17), 6408; https://doi.org/10.3390/s22176408
Submission received: 9 July 2022 / Revised: 22 August 2022 / Accepted: 23 August 2022 / Published: 25 August 2022
(This article belongs to the Section Physical Sensors)

Abstract

:
The fluid inerter described by the fractional derivative model is introduced into the traditional nonlinear energy sink (NES), which is called fractional-order NES in this paper. The slowly varying dynamic equation (SVDE) of the system coupled with fractional-order NES is obtained by the complex averaging method, in which the fractional derivative term is treated using the fractional Leibniz theorem. Then, the discriminants (Δ, Δ1, and Δ2) of the number of equilibrium points are derived. By using the variable substitution method, the characteristic equation for judging the stability is established. The results show: (1) the approximate SVDE is sufficient to reflect the slowly varying characteristics of the system, which shows that the mathematical treatment of the fractional derivative term is reliable; (2) the discriminant conditions (Δ1, Δ2) can accurately reflect the number of equilibrium points, and the corresponding range of nonlinear parameter κ can be calculated when the system has three equilibrium points. The expressions of Δ1, Δ2 are simpler than Δ, which is suitable for analysis and design parameters; (3) the stability discrimination methods of schemes 1 and 2 are accurate. Compared with scheme 2, scheme 1 is more prone to various responses, especially various strongly and weakly modulated responses. In scheme 2, the inertia effect of mass can be completely replaced by integer order inerter. Compared with integer order inerter, the introduction of fractional order inerter, whether in series or in parallel, means that the amplitude of the equilibrium point on the NES vibrator is smaller, but it is also for this reason that it is not easy to produce a modulated response with scheme 2, and the vibration suppression effect of the main structure is not good.

1. Introduction

When a fluid inerter works [1,2,3,4,5], it produces inertia force, usually accompanied by large parasitic damping force, which cannot be ignored. Therefore, the traditional mathematical model of an ideal inerter, such as a mechanical inerter [6,7,8], cannot describe the characteristics of a fluid inerter. The approximate model of its output force is mainly obtained by referring to the empirical formula. For example, the output force model of the fluid inerter proposed in [9,10,11] is obtained using the laminar flow empirical formula, and the output force model is established using the turbulence empirical formula in [4,5]. Some of these models consider laminar flow, and others consider the form of turbulence, but the actual situation may be both. There is still room for discussion on the model of fluid inerters. In reference [12], the author noted that describing the multiphase mechanical properties by simply adding the integer derivatives of different derivative orders is not satisfactory. From the perspective of a modeling fractional differential equation, Westerlund [13] proposed the use of a unified fractional model containing multiphase mechanical properties. In reference [14], a fractional derivative model is proposed to describe the output force model of a fluid inerter. It is verified both mathematically and physically that the model can reflect the multiphase mechanical characteristics of the inerter. Compared with the integer order, an inerter, although it brings difficulties in mathematical processing [15], can have more parameter combinations, to ensure that the system achieves the required performance.
In recent years, NES [16,17] became a research hotspot of nonlinear dynamic vibration absorbers. The typical NES is composed of additional mass, a nonlinear spring, and linear damping, which is far superior to the traditional linear vibration absorber. The effective vibration suppression of NES is mainly through the TET [18,19,20], and its internal mechanism is resonance capture. Only when the external excitation intensity reaches the trigger threshold does the TET become excited. Under harmonic excitation, the study of NES systems is accompanied by the discussion of the mechanism of the strongly modulated response (SMR) generation [16,21,22,23,24]. These research results show that the vibration system has better energy transfer efficiency and vibration suppression effect when SMR appears. The SMR is due to the coupling system slowly varying power flow caused by the saddle node bifurcation of limit cycles, and is a real phenomenon in engineering. The saddle node bifurcation is directly related to the stability of the equilibrium point of the slowly varying power flow. Scholars [24,25,26] used a complex average method, multi-scale method, harmonic balance method, and so on, to study the stability of the equilibrium points. References [27,28,29,30] introduce a new idea that parallel or series inertial ideal integer order inerter based on traditional NES can reduce the added mass and improve the vibration suppression performance. However, the fractional order inerter has not been applied to nonlinear vibration absorbers yet. In this paper, we will introduce the fractional order inerter to NES structure and discuss the influence of nonlinear parameters and inerter parameters on the vibration system.
The remainder of this paper is organized as follows. In Section 2, the dynamic models of the two schemes are established. In Section 3, the SVDE, the discriminant of the number of equilibrium points, and the stability discrimination of the equilibrium points are derived for the two schemes. Section 4 is the numerical verification of the theoretical analysis in Section 3. Section 5 contains the conclusions.

2. Dynamic Model of a Linear SDOF System Coupled with Fractional-Order NES

The two models of the SDOF system coupled with fractional-order NES are shown in Figure 1, where scheme 1 (S.1) is the NES with an inerter in parallel (a); and scheme 2 (S.2) is the NES with one inerter in series (b). In Figure 1, m1, ks, and cs are the main structural mass, stiffness, and damping coefficient, respectively; m2, kn, cn, and bn are the mass, ideal nonlinear cubic stiffness, damping coefficient of the NES, and the inertance of the inerter, respectively. Since this paper discusses the influence of the fractional derivative of the fractional order inerter on the research, the influence of nonlinear factors and Coulomb friction is ignored [10,11,14]. The fractional order inerter model is
F b = b n D μ z  
where μ between 1 and 2, especially when μ = 2, is an ideal integer order inerter. In recent decades, the most used definitions of fractional calculus are the Riemann–Liouville, Caputo, and Grunwald–Letnikov derivatives [31]. Among them, the Caputo definition is suitable for the application of various engineering scenarios. Therefore, Caputo’s definition is adopted for the fractional derivatives in this paper.
Therefore, the dynamic models of the two schemes can be obtained from Newton’s second law, as shown in Equations (2) and (3).
S.1:
{ m s z ¨ 1 + m 2 z ¨ 2 + c s z ˙ 1 + k s z 1 = f cos ω t m 2 z ¨ 2 [ b n D μ ( z 1 z 2 ) + c n ( z ˙ 1 z ˙ 2 ) + k n ( z 1 z 2 ) 3 ] = 0 ,
and
S.2:
{ m s z ¨ 1 + m 2 z ¨ 2 + b n D μ z 2 + c s z ˙ 1 + k s z 1 = f cos ω t m 2 z ¨ 2 + b n D μ z 2 c n ( z ˙ 1 z ˙ 2 ) k n ( z 1 z 2 ) 3 = 0
Obviously, when bn = 0, it is the vibration dynamic model coupled with the traditional NES.
According to reference [32], when the generalized harmonic analysis is performed, the fractional derivative terms in vibration Equations (2) and (3) can be approximately equivalent in the following integer order as follows:
{ D μ z 2 = ω μ 1 cos ( ( μ 1 ) π / 2 ) z ˙ 2 + ω μ 2 sin ( ( μ 1 ) π / 2 ) z ¨ 2 D μ ( z 1 z 2 ) = ω μ 1 cos ( ( μ 1 ) π / 2 ) ( z ˙ 1 z ˙ 2 ) + ω μ 2 sin ( ( μ 1 ) π / 2 ) ( z ¨ 1 z ¨ 2 )
Here, the first term on the right of the above formula plays a damping role, and the second term plays an inertia role. After substituting (4) into Equations (2) and (3), m2 cannot be 0 for S.1, but m2 can be equal to 0 for S.2, and its inertia is completely provided by the inerter.
By introducing a new time scale τ = ω 0 t , ω 0 = k s / m s , Ω = ω / ω 0 , the above equation is dimensionless, and the following equation is obtained:
S.1:
{ z 1 + ε z 2 + ξ s z 1 + z 1 = f e cos Ω τ ε z 2 [ ε 1 ω 0 μ 2 D μ ( z 1 z 2 ) + ξ n ( z 1 z 2 ) + κ ( z 1 z 2 ) 3 ] = 0 ,
and
S.2:
{ z 1 + ε z 2 + ε 1 ω 0 μ 2 D μ z 2 + ξ s z 1 + z 1 = f e cos Ω τ ε z 2 + ε 1 ω 0 μ 2 D μ z 2 ξ n ( z 1 z 2 ) κ ( z 1 z 2 ) 3 = 0
where Ω = ω ω 0 ,   ε = m 2 m s ,   ε 1 = b n m s ,   ξ s = c s m s k s ,   ξ n = c n m s k s ,   κ = k n k s ,   f e = f k s .

3. Equilibrium Point Analysis

3.1. The SVDE of the System Coupled with Fractional-Order NES (S.1)

In this paper, the nonlinear term is not limited to a small quantity, so system (5) can be a strongly nonlinear system. The complex averaging method (CA-X) [25,26] is used to derive the SVDE of system (5). Internal resonance of the system is one of the preconditions to ensure TET. Therefore, this paper studies the characteristics of the equilibrium point at 1:1 internal resonance i.e., Ω =1, and introduces complex variables as follows:
φ j e i τ = z j + i z j , φ j * e i τ = z j i z j ( j = 1 , 2 ) .
The following relations can be obtained by further derivation of the (7):
z j = φ j e i τ φ j * e i τ 2 i , z j = φ j e i τ + φ j * e i τ 2 , z j = φ ˙ 1 e i τ + i φ j e i τ φ j * e i τ 2 .
Equation (8) can be further transformed to obtain the following relationship:
φ ˙ j e i τ φ ˙ j * e i τ = 0 , φ ˙ j e i τ + φ ˙ j * e i τ = 2 φ ˙ 1 e i τ
According to the properties of fractional calculus [31], D μ z j = D μ 1 [ d z j d t ] , that is, D μ z j = D α [ z ˙ j ] , α = μ − 1, α ∈ (0, 1). The fractional Leibniz theorem [31] is introduced as follows:
D t α t 0 [ f ( t ) g ( t ) ] = k = 0 n ( α k ) f ( k ) ( t ) g ( α k ) ( t ) R n α ( t ) ,  
where ( α k ) = Γ ( α + 1 ) k ! Γ ( α k + 1 ) , Γ ( ) is the gamma function, and the functions f(t) and g(t) and their derivatives are continuous on the interval [t0, t]. Since 0 < α < 1, it can be seen from [31] that when n = 2 and k < 2, the remainder R n α ( t ) can be ignored.
Therefore, according to Equations (8)–(10), D μ z j can be reduced to:
D μ z j = D α ( φ j e i τ + φ j * e i τ 2 ) = ( φ j e i τ + φ j * e i τ ) cos α π 2 + ( φ j e i τ φ j * e i τ ) sin α π 2 i 2 + α φ ˙ j e i τ sin α π 2 .
Then, the Equations (7), (8) and (11) are substituted into Equation (5), and averaged to obtain the SVDE of S.1:
{ φ 1 + ε φ 2 + ξ s φ 1 2 + i ε φ 2 2 = f e 2 ( ε + α b ) φ 2 α b φ 1 + i ε φ 2 2 ( c + b i + ξ n 3 κ 4 i | φ 1 φ 2 | 2 ) φ 1 φ 2 2 = 0 .
where b = ε 1 ω 0 μ 2 sin α π 2 , c = ε 1 ω 0 μ 2 cos α π 2 , and the following b and c express the same meaning.
To compare the slow variable model (12) with the zj obtained from the original system (5), the slow variable obtained from (12) is φj, to be converted into a response in system (5). According to (7), zj can be determined by φj, the formula is zj = imag(φje), zj = real(φje), j = 1,2. The fractional derivative term in the original system (5) is solved by using the improved Oustaloup’s method [33].
Figure 2 (ε1 = 0) and Figure 3 (ε1 = 0.05) are displacement responses: (a) main structure (z1), (b) NES (z2), respectively (other parameters: ε = 0.05, μ = 1.9, ω0 = 2.5, ξs = 0.01, ξn = 0.005, fe = 0.01, κ = 1). It can be seen from Figure 2 and Figure 3 that the responses of the slow varying model are in good agreement with the numerical solution in an abbreviated period at the beginning, and has a small deviation from the numerical solution after a long period of time. In general, the SVDE (12) can reflect the change trend of the original system response, and can be used for the comparison of the later stable responses.

3.2. The Discriminant Derivation of the Number of Equilibrium Points of S.1

To discuss the steady-state characteristics of the SVDE (S.1), it is difficult to calculate the equilibrium point of Equation (12) by directly making the derivative term 0. Therefore, the following variables are substituted for Equation (12):
Ψ 1 = φ 1 , Ψ 2 = φ 1 φ 2
Substituting Equation (13) into Equation (12), the following relation can be obtained:
{ Ψ ˙ 1 + ε ( Ψ ˙ 1 Ψ ˙ 2 ) + ξ s Ψ 1 2 + i ε Ψ 1 Ψ 2 2 = f e 2 ( ε + α b ) ( Ψ ˙ 1 Ψ ˙ 2 ) α b Ψ ˙ 1 + i ε Ψ 1 Ψ 2 2 ( c + b i + ξ n 3 κ 4 i | Ψ 2 | 2 ) Ψ 2 2 = 0 .
If the derivative term in the above equation is zero, we obtain the equilibrium equation of the system as follows:
9 16 κ 2 | Ψ 2 | 6 γ 3 3 κ 2 | Ψ 2 | 4 + ( γ 2 2 + γ 3 2 ) | Ψ 2 | 2 = γ 1 2 ,
Ψ 2 = γ 1 ( γ 2 + γ 3 i 3 κ 4 i | Ψ 2 | 2 ) ,
and
Ψ 1 = f e ξ s + i ε + i ε ξ s + i ε Ψ 2 ,
where
γ 1 = ε ξ s 2 + ε 2 f e , γ 2 = ( ξ s ε 2 ξ s 2 + ε 2 + c + ξ n ) , γ 3 = ( ξ s 2 ε ξ s 2 + ε 2 + b ) .
Equations (15a)–(15c) are a set of uncoupled nonlinear equations, and the equilibrium points of the system can be calculated from top to bottom. It can be seen from Equations (15a)–(15c) that (15b) and (15c) have only a single real root, so the number of equilibrium points of the system is determined by the number of real roots of (15a). The number of real roots of the polynomial can be determined by the Cardano discriminant. For equation a 2 x 3 + a 1 x + a 0 = 0 , the Cardano discriminant is:
Δ = ( a 0 2 a 2 ) 2 + ( a 1 3 a 2 ) 3 .
When Δ > 0, (15a) has a real root. When Δ = 0, (15a) has three real roots, of which at least two are equal. When Δ < 0, (15a) has three unequal real roots.
To determine the number of real roots of (15a), it is necessary to replace them with variables, as shown below:
| Ψ 2 | 2 = β v + γ 3 3 κ 2 3 9 16 κ 2 = β v + 8 γ 3 9 κ .
Then, (15a) can be rewritten as:
9 16 κ 2 β v 3 + ( γ 2 2 γ 3 2 3 ) β v + 8 9 κ γ 3 ( γ 2 2 + 1 9 γ 3 2 ) γ 1 2 = 0 .
Therefore:
Δ = ( 64 γ 3 ( γ 2 2 + 1 9 γ 3 2 ) 72 γ 1 2 κ 81 κ 3 ) 2 + ( 16 ( γ 2 2 γ 3 2 3 ) 27 κ 2 ) 3 .
Particularly, when Δ = 0, the following relation can be obtained:
κ 1 , 2 = 8 ( ( 9 γ 2 2 + γ 3 2 ) γ 3 ± ( γ 3 2 3 γ 2 2 ) 3 ) 81 γ 1 2 .
Obviously, the number of κ satisfying Δ = 0 is determined by the positive and negative of (21) as follows:
Δ 1 = γ 3 2 3 γ 2 2
Δ1 is the judgment condition for the quantity whose Δ is equal to 0. When Δ1 = 0, Δ = 0, the system has three equilibrium points, and in other positions when Δ > 0, the system has only one equilibrium point. When Δ1 > 0, Δ has two zero points, and there are three equilibrium points in the system between the two zero points, and there is only one equilibrium point at other positions. When Δ1 < 0, Δ has no zero points, and there is only one equilibrium point at all positions.

3.3. The Discriminant Derivation of the Number of Equilibrium Points of S.2

The derivation method of S.2 is the same as S.1, so the derivation process is simplified, and the complex variables are substituted into (6) to obtain the SVDE of S.2 at 1:1 internal resonance, as shown below:
{ φ 1 + ( ε + α b ) φ 2 + ξ s φ 1 2 + ( c + ( ε + b ) i ) φ 2 2 = f e 2 ( ε + α b ) φ 2 + ( c + ( ε + b ) i ) φ 2 2 ( ξ n 3 κ 4 i | φ 1 φ 2 | 2 ) φ 1 φ 2 2 = 0
Selecting the same variable as in the previous section to replace (Ψ1 = φ1, Ψ2 = φ1φ2), and the SVDE of S.2 is rewritten as:
{ Ψ ˙ 1 + ( ε + α b ) ( Ψ ˙ 1 Ψ ˙ 2 ) + ξ s Ψ 1 2 + ( c + ( ε + b ) i ) ( Ψ 1 Ψ 2 ) 2 = f e 2 ( ε + α b ) ( Ψ ˙ 1 Ψ ˙ 2 ) + ( c + ( ε + b ) i ) ( Ψ 1 Ψ 2 ) 2 ( ξ n 3 κ 4 i | Ψ 2 | 2 ) Ψ 2 2 = 0
If the derivative term in (23) is 0, the equilibrium point equations of S.2 can be obtained as follows:
9 16 κ 2 | Ψ 2 | 6 3 2 γ ¯ 3 κ | Ψ 2 | 4 + γ ¯ 2 | Ψ 2 | 2 = γ ¯ 1
Ψ 2 = K 1 f e ξ s ( ξ n 3 κ 4 | Ψ 2 | 2 i ) ( c B 1 i ) + ( C 1 3 κ 4 | Ψ 2 | 2 i ) K 1
Ψ 1 = ξ n + c + ( B 1 3 κ 4 | Ψ 2 | 2 ) i c + B 1 i Ψ 2
where B 1 = ( ε + b ) , C 1 = ξ s + ξ n , K 1 = c 2 + B 1 2 , γ ¯ 1 = K 1 2 ξ s 2 B 1 2 + ( K 1 + c ξ s ) 2 f e 2 , γ ¯ 2 = ξ s 2 ξ n 2 B 1 2 + ( ξ s ξ n c + K 1 C 1 ) 2 ξ s 2 B 1 2 + ( K 1 + c ξ s ) 2 , γ ¯ 3 = ξ s 2 B 1 K 1 ξ s 2 B 1 2 + ( K 1 + c ξ s ) 2 .
To judge the number of real roots of (24a), comparing (24a) with (15a) to find that the same variable can be used to replace. Substituting (17) into (24a), (24a) can be rewritten as follows:
9 16 κ 2 β v 3 + ( γ ¯ 2 4 γ ¯ 3 2 3 ) β v + 8 γ ¯ 2 γ ¯ 3 9 κ 64 γ ¯ 3 3 81 κ γ ¯ 1 = 0
The discriminant is obtained by substituting the coefficient of (25) into (16):
Δ = ( 8 ( 72 γ ¯ 2 γ ¯ 3 81 κ γ ¯ 1 64 γ ¯ 3 3 ) 729 κ 3 ) 2 + ( ( 48 γ ¯ 2 64 γ ¯ 3 2 ) 81 κ 2 ) 3
Similarly, assuming Δ = 0, the expression of κ can be simplified as follows:
κ 1 , 2 = 8 ( 9 γ 2 8 γ 3 2 ) γ 3 ± 8 ( 4 γ 3 2 3 γ 2 ) 3 81 γ 1
Therefore, the discriminant condition for the number of κ is:
Δ 2 = 4 γ ¯ 3 2 3 γ ¯ 2
To understand the influence of various parameters on the number of equilibrium points when there is no inerter (i.e., bn = 0), according to Equations (2) and (3), the dynamic equations of the two schemes are the same. Therefore, taking S.2 as an example, this paper gives the discriminant Δ2 when bn = 0 by replacing parameters, as shown below:
Δ 2 = 4 ξ s 4 ε 2 3 [ ( ξ s ξ n ) 2 + ( ξ n + ξ s ) 2 ε 2 ] ( ξ s 2 + ε 2 ) ε 2 ( ξ s 2 + ε 2 ) 2
Observing the above formula, the value of Δ2 is determined by ξnξs, and ε, and the clear relationship is that the smaller ξn is, the larger Δ2 is.

3.4. Stability Analysis Based on Lyapunov Theory

The equilibrium obtained from the equations derived in the previous section may be unstable. If the equilibrium point is unstable, the system may have a weakly modulated response (WMR) or SMR. As the response type of the system has a very important impact on the vibration suppression effect, it is necessary to analyze the stability of the equilibrium point of the system.
For S.1, the Equation (14) is expanded near the equilibrium point, i.e.,
Ψ 1 = Ψ 10 + δ 1 , Ψ 2 = Ψ 20 + δ 2
where Ψ10 and Ψ20 are the equilibrium points of the system, and δ1 and δ2 are the small increments near the corresponding equilibrium points. After algebraic operation, the linear differential equations about the increment can be obtained as follows:
{ δ ˙ 1 = ξ s ( ε + α b ) ϒ δ 1 2 α b ϒ i ε δ 1 δ 2 2 ε ϒ ( c + b i + ξ n ) δ 2 2 + ε ϒ 3 κ 8 i ( 2 | Ψ 20 | 2 δ 2 + Ψ 20 2 δ 2 * ) δ ˙ 2 = ξ s ε ϒ δ 1 2 + i ε ϒ δ 1 δ 2 2 1 + ε ϒ ( c + b i + ξ n ) δ 2 2 + 1 + ε ϒ 3 κ 8 i ( 2 | Ψ 20 | 2 δ 2 + Ψ 20 2 δ 2 * )
By substituting the following formula into (31):
Ψ 20 = N 20 e i θ 20 , δ 1 = p 1 + q 1 i , δ 2 = p 2 + q 2 i , δ 2 * = p 2 q 2 i
and separating the real part from the imaginary part, the following relationship can be obtained:
{ p ˙ 1 = ε 2 ϒ ( ξ s ( 1 + b α / ε ) p 1 + b α q 1 ϒ 2 p 2 + ( ϒ 3 b α ) q 2 ) q ˙ 1 = ε 2 ϒ ( b α p 1 ξ s ( 1 + b α / ε ) q 1 + ( b α ϒ 4 ) p 2 ϒ 1 q 2 ) p ˙ 2 = ε 2 ϒ ( ξ s p 1 q 1 ( 1 + 1 / ε ) ϒ 2 p 2 + ( 1 + ( 1 + 1 / ε ) ϒ 3 ) q 2 ) q ˙ 2 = ε 2 ϒ ( p 1 ξ s q 1 ( 1 + ( 1 + 1 / ε ) ϒ 4 ) p 2 ( 1 + 1 / ε ) ϒ 1 q 2 )
where ϒ 1 = c + ξ n 3 N 2 κ sin 2 θ 20 4 , ϒ 2 = c + ξ n + 3 N 2 κ sin 2 θ 20 4 ,   ϒ 3 = b 3 4 ( 2 cos 2 θ 20 ) κ N 20 2 , ϒ 4 = b 3 4 ( 2 + cos 2 θ 20 ) N 20 2 κ .
θ20 in the above equation can be obtained from the following equation:
θ 20 = i ln Ψ 2 | Ψ 2 |
Obviously, Equation (33) is a linear constant system, the stability can be directly determined by Lyapunov theory, and the characteristic factor is defined:
Ξ = max ( Re ( λ i ) )
λi are the eigenvalues of the coefficient matrix on the right side of (33). According to Lyapunov theory, if Ξ < 0, the equilibrium point of the system is asymptotically stable; if Ξ = 0, the equilibrium point is critically stable; if Ξ > 0, the equilibrium point is unstable. If a pair of eigenvalues cross the imaginary axis when they reach a certain point with a certain parameter, then this point is also called the Hopf bifurcation point of the system.
The characteristic matrix of the system is as follows:
A 1 = ε 2 ϒ ( ξ s ( 1 + b 1 / ε ) b 1 ϒ 2 ϒ 3 b 1 b 1 ξ s ( 1 + b 1 / ε ) b 1 ϒ 4 ϒ 1 ξ s 1 ( 1 + 1 / ε ) ϒ 1 1 + ( 1 + 1 / ε ) ϒ 3 1 ξ s 1 ( 1 + 1 / ε ) ϒ 4 ( 1 + 1 / ε ) ϒ 2 )
Therefore, according to the following formula:
| λ I 4 × 4 A 1 | = 0
its characteristic polynomial can be obtained as follows:
a 1 λ 4 + a 2 λ 3 + a 3 λ 2 + a 4 λ 1 + a 5 = 0
Since the forms of ai (i = 1, 2, …5) are too complex, and the eigenvalues in (38) can be easily calculated in MATLAB, the specific forms are not given here. For S.2, using the same method above, the system matrix A2 of its incremental state equation can be derived from (23) as follows:
A 2 = ( 1 2 ξ s 0 2 3 0 ξ s 2 4 1 c ξ s 2 b + ε 2 c 2 1 + 2 b + ε 2 + 1 + 3 b + ε 2 c ξ s 2 b + ε 2 + 1 + 4 1 2 + 1 + 1 )
where = ( b α + ε ) , 1 = 3 N 2 κ sin 2 θ 4 ξ n 8 , 2 = 3 N 2 κ sin 2 θ + 4 ξ n 8 , 3 = 3 N 2 κ ( cos 2 θ 2 ) 8 , 4 = 3 κ N 2 ( cos 2 θ + 2 ) 8 .
The stability of the equilibrium point of S.2 is determined by the eigenvalue of the characteristic equation corresponding to the system matrix A2.

4. Numerical Simulation and Discussion

Since the cubic stiffness of the NES is not much different from the linear stiffness of the actual structure, this section only studies the characteristics of the system equilibrium point when κ changes within a certain range, but the method used is not limited by parameter changes. The system parameters discussed in this section are normalized parameters, in which the constant parameters ms = 10, ks = 0.6, and fe = 0.01.

4.1. Influence of Parameters on the Number of Equilibrium Point

To find the parameter combination of three equilibrium points in the system without inerter (bn = 0), Figure 4a–c are drawn according to Δ2 (28), where Figure 4a shows that the larger cn, the smaller Δ2, which is consistent with the previous prediction. Figure 4b shows that when m2 increases, Δ2 first increases and then decreases, while the Figure 4c shows that the influence of cs is just opposite to m2. The greater the cs, Δ2 first decreases and then increases. To make Δ2 > 0, m2 should choose a small value, and cs should choose a relatively large value.
Then, selecting the appropriate parameters, and considering the NES of parallel and series inerter, the effects of parameters μ and bn on Δ1 and Δ2, respectively, are discussed. Considering the influence of inerter (bn ≠ 0), for S.1, Figure 5a–c are drawn according to Δ1 (21), as shown below.
Figure 5a shows that the two surfaces almost coincide, and m2 has an insignificant effect on Δ1, which is in line with the prediction of the curve when cs = 0.003 in Figure 4b. In Figure 5b, the change of μ directly changes the slope of Δ1 with bn, especially when μ = 5/3, the slope is 0. As shown in Figure 5c, bn has no effect on the μ range of Δ1 > 0, which is always between μ > 5/3 and 2. However, the larger bn is, the larger Δ1 is when μ > 5/3. In addition, when there is no inerter (bn = 0), the slope is 0, and Δ1 remains unchanged and less than 0, the system does not have multiple equilibrium points. Therefore, in the following discussion, μ should be selected between (5/3,2).
When S.2 adopt the same parameters as S.1, Δ2 is always negative, indicating that it is not easy to make Δ2 > 0 in S.2. To highlight the role of the inerter, for S.2, m2 = 0, and the influence of the inerter parameters on the value of Δ2 is analyzed. Under the premise of m2 = 0, the surfaces of Δ2 versus μ and bn is drawn for different cs (cs = 0.1, cs = 0.3) in Figure 6a; Figure 6b shows the curve of Δ2 versus bn for different μ (cs = 0.3); Figure 6c shows the curve of Δ2 versus μ for different bn (m2 = 0, cn = 0.002, bn = 0.4), and other parameters are the same as those in S.1.
As can be seen from Figure 6a,c, in the case of other parameter phase diagrams, the larger the cs, the larger the range of Δ2 greater than zero. It can be seen from Figure 6b that the change law of the curve is consistent with that of Figure 4b, which shows that the inertance (bn) of the series inerter can nominally replace m2. However, for S.1 the change law is different from that of the curve in Figure 4b, which shows that bn cannot simply replace m2 when the inerter is connected in parallel.

4.2. Verification of Equilibrium Point Discriminant

To verify the positive and negative correspondence between the number of equilibrium points and Δ, and the prediction results of Δ1 and Δ2 in the previous section, the curves of Δ and |Ψ2| with κ are drawn.
For S.1, the values of m2 = 0.3 and m2 = 4 are selected, and the other parameters are consistent with those in Figure 5b (cs = 0.003, cn = 0.002, bn = 0.4). The curves of Δ (Equation (19)) and |Ψ2|of S.1 versus μ under different κ are drawn, as shown in Figure 7 and Figure 8.
In S.2, the parameters of Figure 6b (m2 = 0, cs = 0.3, bn = 0.4, cn = 0.002) and Figure 6c (m2 = 0, cs = 0.1, bn = 0.4, cn = 0.002) in the previous section are selected, and the curves of (a) Δ (Equation (26)) and (b) |Ψ2| of S.2 versus κ for different μ are drawn, as shown in Figure 9 and Figure 10.
Whether in S.1 or S.2, the positive and negative of Δ correspond to the number of equilibrium points of the system one by one, and can also correspond to the positive and negative prediction of Δ1 or Δ2 in the previous section. Under the same parameters, the smaller μ, the smaller the amplitude of the equilibrium point, and the less likely it is to switch the number of equilibrium points, indicating that the inertia and damping effects of the inerter exist simultaneously when μ is between 1 and 2. Here, κ1,2 = (0.22, 0.34) is calculated by substituting the parameters into (20) when μ = 1.8 in Figure 7b, and κ1,2 = (0.77, 1.26) is calculated by substituting the parameters into (27) when μ = 2 in Figure 9b. Obviously, the calculated results of the two groups of parameters are consistent with the positions of the actual simulation curves. In addition, as predicted in the previous section, it is not easy to produce multiple equations with S.2. Particularly, when Δ1 = 0 or Δ2 = 0, the system is always only one equilibrium point. It can be considered that the system is just in a critical state.

4.3. Stability Discrimination of Equilibrium Points

In this section, the Lyapunov method in Section 3.3 is used to judge the stability of the equilibrium points of S.1 and S.2, and to verify whether the calculation results of the equilibrium point are consistent with those calculated by the numerical algorithm. Of course, this is meaningful only when the equilibrium point is asymptotically stable.
Taking S.1 as an example, the verification method for slowly varying parameters is established. According to (7), after the equilibrium point of system (14) is calculated by (15), the state of the original system (5) corresponding to the equilibrium point when t = 0 can be expressed as:
[ z 1 , z 1 , z 2 , z 2 ] = [ Im ( φ 1 ) , Re ( φ 1 ) , Im ( φ 2 ) , Re ( φ 2 ) ]
A small disturbance (δ1, δ2, δ3, δ4) is added to (40) and brought into the original system (5) as an initial condition. If the equilibrium point of the system is stable, the time response is stable after a period.
If Ψ2 is expanded in the complex plane with N2e, then |Ψ2| = N2. The steady-state amplitude |Ψ2| is obtained from (15a). Since Ψ2 = φ1φ2, the response of the original system can be converted to Ψ2, that is, Ψ2 = −z′1z′2 + (z1z2)i. The original system (5) is calculated by the fourth order Runge–Kutta (R–K) method, and its steady-state value N2 can be obtained through conversion, as shown below:
N 2 = ( z 1 z 2 ) 2 + ( z ˙ 1 z ˙ 2 ) 2 ( τ 1 , τ 2 )
where < >(τ1, τ2) represents the weighted average in the (τ1, τ2) period. Obviously, S.2 can also be treated in the same way.
Adding a disturbance δ1 = 0.01 to z1(0), the stability of the equilibrium point of the original system is analyzed. For S.1, the different stabilities of the equilibrium points are drawn, as shown in Figure 11, and the parameters of Figure 7b and Figure 8b are used in Figure 11a,b, respectively. For S.2, the different stabilities of the equilibrium points are drawn, as shown in Figure 12, and the parameters of Figure 9b and Figure 10b are used in Figure 12a,b, respectively. In Figure 11 and Figure 12, the red squares represents the stable equilibrium points obtained by the R–K method and the black dots represents the stable equilibrium points obtained by (15a). The other represent unstable equilibrium points are also represented.
Firstly, from the above four pictures, the stable equilibrium points obtained by the R–K method are the same as that obtained by the calculation in (15a) or (24a), which verifies the accuracy of the analytical method. No matter which scheme, when the system has three equilibrium points, the switching of the stability always occurs near the turning point of the equilibrium point curve, and Δ1 and Δ2 proposed in this paper can only predict the turning point position, and the parameters are few and easy to analyze.
Figure 11a,b show that both schemes have three equilibrium points, but Figure 11a shows two unstable equilibrium points with large amplitude and one stable equilibrium point with small amplitude, and Figure 11b shows an unstable equilibrium point sandwiched between the two stable equilibrium points, indicating that there are two different bifurcation phenomena in the system. When bn is calculated according to Δ1 = 0, or Δ2 = 0 is taken as the inertance of the system, the equilibrium point of the system is always stable with the change of κ. It is more difficult to produce multiple equilibrium points in S.2 (Figure 12) than S.1 (Figure 11), and the range of unstable parameters is relatively small.

4.4. The Time-Domain Response

Next, the time-domain responses are drawn to verify the accuracy of the above equilibrium stability analysis results, and the influence of equilibrium stability on the response type when the slow varying system has multiple solutions is explored. Particularly, to maintain the consistency between the steady-state amplitude N2 of the system and |Ψ2| in the previous section, “z1z2” is drawn instead of “z2”.
For S.1, we selected the parameters when k = 0.11 in Figure 11a (μ = 2 and μ = 1.9), and plotted the time-domain response of the original system (Figure 13 (μ = 2) and Figure 14 (μ = 1.9)) under different initial interference (δ1).
In Figure 13a–d, the corresponding initial amplitude|Ψ2|(0) is (0.11, 0.41, 0.56, 0.91), just between the three equilibrium points (|Ψ2|1,2,3 = 0.3109, 0.49, 0.79). The results show that the response can be stabilized at the small stable equilibrium point only when it is between the small equilibrium point and the intermediate equilibrium point, while the system response is in the SMR when other disturbances occur, which is consistent with the conclusion when μ = 2, κ = 0.11 in Figure 11a. The four δ1 are shown in Figure 14a–d, and the corresponding initial amplitudes (|Ψ2|(0)) are (0.07, 0.43, 0.73, 0.93) are just between the three equilibrium points (|Ψ2|1,2,3 = 0.24, 0.63, 0.81). If the initial amplitude is less than the second largest equilibrium point, the system response asymptotically tends to the small stable equilibrium point after a long time of WMR (as shown in Figure 14a,b), while when it is greater than the second largest equilibrium point, the system response is in the SMR (as shown in Figure 14c,d).
Then, the parameters (m2 = 4, κ = 0.1, μ = 2) and (m2 = 4, κ = 0.1, μ = 1.9) in Figure 11b are selected to draw the time-domain response of the original system (5) (as shown in Figure 15 and Figure 16). Figure 15a–d show that no matter what kind of initial interference, the system response tends to be stable after a long time of WMR. When the initial value is between two stable equilibrium points, it is always stable at the small equilibrium point, otherwise, it is stable at the large stable equilibrium point. The four δ1 are shown in Figure 16a–d, and the corresponding initial amplitudes (|Ψ2|(0)) are (0.06, 0.33, 0.73, 0.93) are just between the three equilibrium points (|Ψ2|1,2,3 = 0.24, 0.67, 0.84). Basically consistent with the situation of μ = 2 (Figure 15), the system response can be stabilized at two stable equilibrium points, but the time to reach the stable equilibrium point is much shorter, and there is no obvious modulated process.
For S.2: from Figure 10a,b in Section 4.2 and Figure 12a,b in Section 4.3, we can see that it is not easy to maintain multiple solutions, and the stability of the three equilibrium points is the same as that of S.1 (Figure 11b). To verify the special case of S.2, that is, whether the inerter can replace the mass function of traditional NES, we selected three groups of parameter schemes for comparison: T1 (μ = 2, m2 = 0, bn = 0.4), T2 (bn = 0, m2 = 0.4), T3 (μ = 1.9, m2 = 0, bn = 0.4) and the other parameters are selected when κ = 1.1 in Figure 12a.
Figure 17a,b, and c show the phase diagrams of the three parameter schemes of the original system (6). The initial interference selection method is the same as that of S.1, and the four initial interferences (δ1) are (−0.01, 0.0, 8, 0.14, 0.3). Figure 18a,b show the displacement responses (z1, z1z2) of T1 and T2 when the interference is −0.01 and 0.14, respectively. Figure 18c,d show the displacement (z1, z1z2) response of T3 and T2 when the interference is −0.01 and 0.14, respectively.
It can be seen from Figure 18c,d that the time-domain responses of T3 and T2 are basically coincidental, which shows that the bn of S.2 can replace the role of m2 equivalently when the integer order inerter is connected in series. From Figure 17a and Figure 18a,b, T1 and T2 can be stabilized at small stable equilibrium points no matter how large the interference (δ1) is. As can be seen from Figure 18b, when the initial amplitude is above the small equilibrium point, the amplitude of NES jumps from the large equilibrium point to the small equilibrium point due to the fractional order inerter in series in S.2, resulting in the inability to effectively suppress the displacement of the main structure. At this time, the effect of vibration suppression is not as good as that of series integer order inerter.
The analysis of different situations of the above two schemes shows that the stability discrimination in the previous section is accurate. No matter the scheme, μ has a great impact on response, it is sensitive to the initial value when only the inertia effect (i.e., μ = 2) exists, and it is easy for various modulated responses to appear in the time-domain response. Under the same conditions, S.1 is more sensitive to the nonlinear parameter κ than S.2, and it is easy to produce different system responses.

5. Conclusions

In this paper, the inerter described by fractional derivative is used in traditional NES for the first time. Combined with the CA-X method and fractional properties, especially fractional Leibniz theorem, the difficulty of deriving SVDE from fractional nonlinear dynamic equations is solved. Through the verifications of analytical and numerical methods, the accuracy of SVDE ensures that the analysis is reliable.
The discriminants Δ1 (S.1) and Δ2 (S.2) for the number of equilibrium points are proposed. Compared with Δ, their expressions are simpler, which can directly and accurately predict the number of equilibrium points of the equilibrium point equation, and the range of κ corresponding to the system with three equilibrium points can be calculated by (20) and (27).
The verifications in Figure 11, Figure 12, Figure 13, Figure 14, Figure 15, Figure 16, Figure 17 and Figure 18 show that the Lyapunov theory given in this paper is accurate in judging the stability of the equilibrium point. For S.2, when m2 = 0, the introduction of the inerter does not increase the complexity of the structure, and the inerter can not only provide inertia, but also provide damping. S.1 has two unstable equilibrium points or two stable equilibrium points at the same time, so S.1 is more prone to various responses, especially various WMR and SMR. Compared with integer order inerter, the introduction of fractional order inerter, whether in series or in parallel, means that the amplitude of the equilibrium point on the NES vibrator is smaller, but it is also for this reason that it is not easy to produced modulated responses in S.2, and the vibration suppression effect of the main structure is not good.

Author Contributions

Writing—original draft preparation, Y.C.; writing—review and editing, Y.C., Y.T. and J.X.; supervision, X.X.; supervision, N.C. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China (Grant no. 11272159) and (Grant no. 51605228).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare that there are no conflict of interest regarding the publication of this paper.

References

  1. Swift, S.J.; Smith, M.C.; Glover, A.R.; Papageorgiou, C.; Gartner, B.; Houghton, N.E. Design and modelling of a fluid inerter. Int. J. Control 2013, 86, 2035–2051. [Google Scholar] [CrossRef]
  2. Pei, D.J.W.J.-S. Modeling a helical fluid inerter system with time-invariantmem-models. Struct. Control Health Monit. 2020, 27, e2579. [Google Scholar] [CrossRef]
  3. Shen, Y.J.; Shi, D.H.; Chen, L.; Liu, Y.L.; Yang, X.F. Modeling and experimental tests of hydraulic electric inerter. Sci. China Technol. Sci. 2019, 62, 2161–2169. [Google Scholar] [CrossRef]
  4. De Domenico, D.; Ricciardi, G.; Zhang, R.F. Optimal design and seismic performance of tuned fluid inerter applied to structures with friction pendulum isolators. Soil Dyn. Earthq. Eng. 2020, 132, 106099. [Google Scholar] [CrossRef]
  5. De Domenico, D.; Deastra, P.; Ricciardi, G.; Sims, N.D.; Wagg, D.J. Novel fluid inerter based tuned mass dampers for optimised structural control of base-isolated buildings. J. Franklin I. 2019, 356, 7626–7649. [Google Scholar] [CrossRef]
  6. Zhao, Z.P.; Chen, Q.J.; Zhang, R.F.; Pan, C.; Jiang, Y.Y. Energy dissipation mechanism of inerter systems. Int. J. Mech. Sci. 2020, 184, 105845. [Google Scholar] [CrossRef]
  7. Xu, K.; Bi, K.M.; Ge, Y.J.; Zhao, L.; Han, Q.; Du, X.L. Performance evaluation of inerter-based dampers for vortex-induced vibration control of long-span bridges: A comparative study. Struct. Control Health 2020, 27, e2529. [Google Scholar] [CrossRef]
  8. Shen, Y.J.; Jiang, J.Z.; Neild, S.A.; Chen, L. Vehicle vibration suppression using an inerter-based mechatronic device. Proc. Inst. Mech. Eng. Part D J. Automob. Eng. 2020, 234, 095440702090924. [Google Scholar] [CrossRef]
  9. Deastra, P.; Wagg, D.J.; Sims, N.D. The Realisation of an Inerter-Based System Using Fluid Inerter. Dyn. Civ. Struct. 2019, 2, 127–134. [Google Scholar] [CrossRef]
  10. Shen, Y.J.; Chen, L.; Liu, Y.L.; Zhang, X.L.; Yang, X.F. Optimized modeling and experiment test of a fluid inerter. J. Vibroeng. 2016, 18, 2789–2800. [Google Scholar] [CrossRef] [Green Version]
  11. Shen, Y.J.; Chen, L.; Liu, Y.L.; Zhang, X.L. Influence of fluid inerter nonlinearities on vehicle suspension performance. Adv. Mech. Eng. 2017, 9, 168781401773725. [Google Scholar] [CrossRef]
  12. Chen, W. Fractional Derivative Modeling of Mechanical and Engineering Problems; Science Press: Beijing, China, 2010. [Google Scholar]
  13. Westerlund, S. Report No. 940426; University of Kalmar: Kalmar, Sweden, 1994. [Google Scholar]
  14. Chen, Y.; Xu, J.; Tai, Y.; Xu, X.; Chen, N. Critical damping design method of vibration isolation system with both fractional-order inerter and damper. Mech. Adv. Mater. Struct. 2022, 29, 1348–1359. [Google Scholar] [CrossRef]
  15. Mu’lla, M.A.M. Fractional Calculus, Fractional Differential Equations and Applications. OALib 2020, 7, 1–9. [Google Scholar] [CrossRef]
  16. AL-Shudeifat, M.A.; Saeed, A.S. Comparison of a modified vibro-impact nonlinear energy sink with other kinds of NESs. Meccanica 2021, 56, 735–752. [Google Scholar] [CrossRef]
  17. Ding, H.; Chen, L.Q. Designs, analysis, and applications of nonlinear energy sinks. Nonlinear Dynam. 2020, 100, 3061–3107. [Google Scholar] [CrossRef]
  18. Li, X.; Zhang, Y.W.; Ding, H.; Chen, L.Q. Dynamics and evaluation of a nonlinear energy sink integrated by a piezoelectric energy harvester under a harmonic excitation. J. Vib. Control 2019, 25, 851–867. [Google Scholar] [CrossRef]
  19. Chen, J.E.; Zhang, W.; Yao, M.H.; Liu, J. Vibration suppression for truss core sandwich beam based on principle of nonlinear targeted energy transfer. Compos. Struct. 2017, 171, 419–428. [Google Scholar] [CrossRef]
  20. Phan, T.N.; Bader, S.; Oelmann, B. Performance of An Electromagnetic Energy Harvester with Linear and Nonlinear Springs under Real Vibrations. Sensors 2020, 20, 5456. [Google Scholar] [CrossRef]
  21. Lamarque, C.H.; Thouverez, F.; Rozier, B.; Dimitrijevic, Z. Targeted energy transfer in a 2-DOF mechanical system coupled to a non-linear energy sink with varying stiffness. J. Vib. Control 2017, 23, 2567–2577. [Google Scholar] [CrossRef]
  22. Li, T.; Seguy, S.; Berlioz, A. Optimization mechanism of targeted energy transfer with vibro-impact energy sink under periodic and transient excitation. Nonlinear Dynam. 2017, 87, 2415–2433. [Google Scholar] [CrossRef] [Green Version]
  23. Li, H.Q.; Li, A.; Zhang, Y.F. Importance of gravity and friction on the targeted energy transfer of vibro-impact nonlinear energy sink. Int. J. Impact Eng. 2021, 157, 104001. [Google Scholar] [CrossRef]
  24. Shao, J.W.; Luo, Q.M.; Deng, G.M.; Zeng, T.; Yang, J.; Wu, X.; Jin, C. Experimental study on influence of wall acoustic materials of 3D cavity for targeted energy transfer of a nonlinear membrane absorber. Appl. Acoust. 2021, 184, 108342. [Google Scholar] [CrossRef]
  25. Manevitch, L.I.; Musienko, A.I.; Lamarque, C.H. New analytical approach to energy pumping problem in strongly nonhomogeneous 2dof systems. Meccanica 2007, 42, 77–83. [Google Scholar] [CrossRef]
  26. Gendelman, O.V.; Starosvetsky, Y.; Feldman, M. Attractors of harmonically forced linear oscillator with attached nonlinear energy sink I: Description of response regimes. Nonlinear Dynam. 2008, 51, 31–46. [Google Scholar] [CrossRef]
  27. Zhang, Z.; Lu, Z.Q.; Ding, H.; Chen, L.Q. An inertial nonlinear energy sink. J. Sound Vib. 2019, 450, 199–213. [Google Scholar] [CrossRef]
  28. Zhang, Y.W.; Lu, Y.N.; Zhang, W.; Teng, Y.Y.; Yang, H.X.; Yang, T.Z.; Chen, L.Q. Nonlinear energy sink with inerter. Mech. Syst. Signal Pr. 2019, 125, 52–64. [Google Scholar] [CrossRef]
  29. Wang, J.J.; Wang, B.; Zhang, C.; Liu, Z.B. Effectiveness and robustness of an asymmetric nonlinear energy sink-inerter for dynamic response mitigation. Earthq. Eng. Struct. D 2021, 50, 1628–1650. [Google Scholar] [CrossRef]
  30. Javidialesaadi, A.; Wierschem, N.E. An inerter-enhanced nonlinear energy sink. Mech. Syst. Signal Pr. 2019, 129, 449–454. [Google Scholar] [CrossRef]
  31. Podlubný, I. Fractional Differential Equations; Academic Press: New York, NY, USA, 1999. [Google Scholar]
  32. Shen, Y.J.; Yang, S.P.; Xing, H.J.; Ma, H.X. Primary resonance of Duffing oscillator with two kinds of fractional-order derivatives. Int. J. Nonlin Mech. 2012, 47, 975–983. [Google Scholar] [CrossRef]
  33. Xue, D.Y.; Zhao, C.N.; Chen, Y.Q. A modified approximation method of fractional order system. In Proceedings of the 2006 IEEE International Conference on Mechatronics and Automation, Luoyang, China, 25–28 June 2006; Volume 1–3, pp. 1043–1048. [Google Scholar]
Figure 1. Two vibration models of the SDOF system coupled with NES: (a) (S.1); (b) (S.2).
Figure 1. Two vibration models of the SDOF system coupled with NES: (a) (S.1); (b) (S.2).
Sensors 22 06408 g001
Figure 2. Displacement responses: (a) main structure, (b) NES (ε1 = 0).
Figure 2. Displacement responses: (a) main structure, (b) NES (ε1 = 0).
Sensors 22 06408 g002
Figure 3. Displacement responses: (a) main structure, (b) NES (ε1 = 0.05).
Figure 3. Displacement responses: (a) main structure, (b) NES (ε1 = 0.05).
Sensors 22 06408 g003
Figure 4. (a) Δ2 vs. cs and m2 at different cn; (b) Δ2 vs. m2 at different cs (bn = 0, cn = 0.002); (c) Δ2 vs. cs at different m2 (bn = 0, cn = 0.002).
Figure 4. (a) Δ2 vs. cs and m2 at different cn; (b) Δ2 vs. m2 at different cs (bn = 0, cn = 0.002); (c) Δ2 vs. cs at different m2 (bn = 0, cn = 0.002).
Sensors 22 06408 g004
Figure 5. (a) Δ1 vs. μ and bn for different m2; (b) Δ1 vs. bn for different μ (m2 = 0.3, cs = 0.003, cn = 0.002); (c) Δ1 vs. μ for different bn (m2 = 0.3, cs = 0.003, cn = 0.002).
Figure 5. (a) Δ1 vs. μ and bn for different m2; (b) Δ1 vs. bn for different μ (m2 = 0.3, cs = 0.003, cn = 0.002); (c) Δ1 vs. μ for different bn (m2 = 0.3, cs = 0.003, cn = 0.002).
Sensors 22 06408 g005
Figure 6. (a) Δ2 vs. μ and bn for different cs; (b) Δ2 vs. bn for different μ (cs = 0.3); (c) Δ2 vs. bn for different cs (m2 = 0, cn = 0.002, bn = 0.4).
Figure 6. (a) Δ2 vs. μ and bn for different cs; (b) Δ2 vs. bn for different μ (cs = 0.3); (c) Δ2 vs. bn for different cs (m2 = 0, cn = 0.002, bn = 0.4).
Sensors 22 06408 g006
Figure 7. (a) Δ and (b)|Ψ2| of S.1 vs. κ for different μ (m2 = 0.3).
Figure 7. (a) Δ and (b)|Ψ2| of S.1 vs. κ for different μ (m2 = 0.3).
Sensors 22 06408 g007
Figure 8. (a) Δ and (b)|Ψ2| of S.1 vs. κ for different μ (m2 = 4).
Figure 8. (a) Δ and (b)|Ψ2| of S.1 vs. κ for different μ (m2 = 4).
Sensors 22 06408 g008
Figure 9. (a) Δ and (b) |Ψ2| of S.2 vs. κ for different μ (cs = 0.3).
Figure 9. (a) Δ and (b) |Ψ2| of S.2 vs. κ for different μ (cs = 0.3).
Sensors 22 06408 g009
Figure 10. (a) Δ and (b) |Ψ2| of S.2 vs. κ for different μ (cs = 0.1).
Figure 10. (a) Δ and (b) |Ψ2| of S.2 vs. κ for different μ (cs = 0.1).
Sensors 22 06408 g010
Figure 11. The different stabilities of the equilibrium points (S.1): (a) m2 = 0.3; (b) m2 = 4.
Figure 11. The different stabilities of the equilibrium points (S.1): (a) m2 = 0.3; (b) m2 = 4.
Sensors 22 06408 g011
Figure 12. The different stabilities of the equilibrium points (S.2) (a) cs = 0.3; (b) cs = 0.1.
Figure 12. The different stabilities of the equilibrium points (S.2) (a) cs = 0.3; (b) cs = 0.1.
Sensors 22 06408 g012
Figure 13. Time-domain response of (5) (parameters in Figure 11a m2 = 0.3, μ = 2, κ = 0.11): (a) δ1 = 0.2, (b) δ1 = −0.1, (c) δ1 = 0.25, (d) δ1 = −0.6.
Figure 13. Time-domain response of (5) (parameters in Figure 11a m2 = 0.3, μ = 2, κ = 0.11): (a) δ1 = 0.2, (b) δ1 = −0.1, (c) δ1 = 0.25, (d) δ1 = −0.6.
Sensors 22 06408 g013
Figure 14. Time-domain response of (5) (parameters in Figure 11a m2 = 0.3, μ = 1.9, κ = 0.11): (a) δ1 = 0.2, (b) δ1 = −0.2, (c) δ1 = −0.5, (d) δ1 = −0.7.
Figure 14. Time-domain response of (5) (parameters in Figure 11a m2 = 0.3, μ = 1.9, κ = 0.11): (a) δ1 = 0.2, (b) δ1 = −0.2, (c) δ1 = −0.5, (d) δ1 = −0.7.
Sensors 22 06408 g014
Figure 15. Time-domain response of (5) (parameters in Figure 11b m2 = 4, μ = 2, κ = 0.1): (a) δ1 = 0.2, (b) δ1 = −0.1, (c) δ1 = 0.25, (d) δ1 = −0.6.
Figure 15. Time-domain response of (5) (parameters in Figure 11b m2 = 4, μ = 2, κ = 0.1): (a) δ1 = 0.2, (b) δ1 = −0.1, (c) δ1 = 0.25, (d) δ1 = −0.6.
Sensors 22 06408 g015
Figure 16. Time-domain response of (5) (parameters in Figure 11b m2 = 4, μ = 1.9, κ = 0.1): (a) δ1 = 0.2, (b) δ1 = −0.1, (c) δ1 = 0.5, (d) δ1 = −0.7.
Figure 16. Time-domain response of (5) (parameters in Figure 11b m2 = 4, μ = 1.9, κ = 0.1): (a) δ1 = 0.2, (b) δ1 = −0.1, (c) δ1 = 0.5, (d) δ1 = −0.7.
Sensors 22 06408 g016
Figure 17. Phase diagrams (κ = 1.1) of (6) at different δ1: (a) T1, (b) T2, (c) T3.
Figure 17. Phase diagrams (κ = 1.1) of (6) at different δ1: (a) T1, (b) T2, (c) T3.
Sensors 22 06408 g017
Figure 18. Time-domain responses (κ = 1.1) of (6) with different parameter schemes (a) T1 and T2 (δ1 = −0.01), (b) T1 and T2 (δ1 = −0.14), (c) T3 and T2 (δ1 = −0.01), (d) T3 and T2 (δ1 = −0.14).
Figure 18. Time-domain responses (κ = 1.1) of (6) with different parameter schemes (a) T1 and T2 (δ1 = −0.01), (b) T1 and T2 (δ1 = −0.14), (c) T3 and T2 (δ1 = −0.01), (d) T3 and T2 (δ1 = −0.14).
Sensors 22 06408 g018
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Chen, Y.; Tai, Y.; Xu, J.; Xu, X.; Chen, N. Vibration Analysis of a 1-DOF System Coupled with a Nonlinear Energy Sink with a Fractional Order Inerter. Sensors 2022, 22, 6408. https://doi.org/10.3390/s22176408

AMA Style

Chen Y, Tai Y, Xu J, Xu X, Chen N. Vibration Analysis of a 1-DOF System Coupled with a Nonlinear Energy Sink with a Fractional Order Inerter. Sensors. 2022; 22(17):6408. https://doi.org/10.3390/s22176408

Chicago/Turabian Style

Chen, Yandong, Yongpeng Tai, Jun Xu, Xiaomei Xu, and Ning Chen. 2022. "Vibration Analysis of a 1-DOF System Coupled with a Nonlinear Energy Sink with a Fractional Order Inerter" Sensors 22, no. 17: 6408. https://doi.org/10.3390/s22176408

APA Style

Chen, Y., Tai, Y., Xu, J., Xu, X., & Chen, N. (2022). Vibration Analysis of a 1-DOF System Coupled with a Nonlinear Energy Sink with a Fractional Order Inerter. Sensors, 22(17), 6408. https://doi.org/10.3390/s22176408

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop