Analytically optimal parameters of fractional-order dynamic vibration absorber

Shen, Yongjun; Peng, Haibo; Wen, Shaofang; Yang, Shaopu; Xing, Haijun

doi:10.21595/jve.2016.16617

Journal of Vibroengineering

Browse Journal

Submit article

Published: 15 August 2016

Check for updates

Analytically optimal parameters of fractional-order dynamic vibration absorber

Yongjun Shen¹

Haibo Peng²

Shaofang Wen³

Shaopu Yang⁴

Haijun Xing⁵

^{1, 4, 5}Department of Mechanical Engineering, Shijiazhuang Tiedao University, Shijiazhuang, 050043, China

²Department of Engineering Mechanics, Shijiazhuang Tiedao University, Shijiazhuang, 050043, China

³Transportation Institute, Shijiazhuang Tiedao University, Shijiazhuang, 050043, China

Corresponding Author:

Yongjun Shen

Cite the article Download PDF

Downloads 1656

WoS Core Citations 3

CrossRef Citations 6

Abstract

In this paper the optimal parameters of the fractional-order Voigt type dynamic vibration absorber (DVA) are analytically studied for two cases, named as $H_{\infty}$ and $H_{2}$ optimization criteria. At first the approximately analytical solution is obtained by the averaging method when the primary system is subjected to harmonic excitation. Then the optimal fractional coefficient and order are obtained based on $H_{\infty}$ optimization criterion, which is designed to minimize the maximum amplitude magniﬁcation factor of the primary system. Based on H₂ optimization criterion, the optimal fractional parameters are obtained to reduce the total vibration energy of the primary system over the whole-frequency range. The comparisons of the approximate solutions with the numerical ones in the two cases are fulfilled, and the results verify that the approximately analytical solutions are correct and satisfactorily precise. At last the control performance of the fractional-order Voigt type DVA is compared with the classical integer-order counterpart, and it could be concluded that the fractional-order DVA has superiority in vibration engineering, and fractional-order element could replace the traditional damper and spring simultaneously in some cases.

1. Introduction

The dynamic vibration absorber (DVA), also named as tuned mass damper (TMD), is a vibration control device which is attached to a vibrating primary system in order to reduce its response by appropriately designing the parameters of DVA. In 1909 Frahm [1] invented the first DVA without damping element, and it could only work in a narrow frequency range close to the natural frequency of the primary system. In 1928, Den Hartog and Ormondroyd [2] presented a DVA with damping element, and found that it could suppress the amplitude of the primary system in a broader frequency range. The DVA by Den Hartog and Ormondroyd was very typical and called as the Voigt type DVA. Because it is efficient, reliable and low-cost, the Voigt type DVA has been used in many fields of engineering practice. In order to improve the control performance, it is necessary to study the optimal parameters of the Voigt type DVA.

At present the study on Voigt type DVA can be mainly divided into the three kinds of optimization criteria, i.e. $H_{\infty}$ , $H_{2}$ , and stability maximization. The purpose of $H_{\infty}$ optimization criterion is to minimize the maximum amplitude magnification factor. This optimization criterion started from the research by Ormondroyd and Den Hartog [2], when they found that the amplitude of the primary system in the damped DVA would pass through two fixed points independent of the absorber damping. Hahnkamm [3] obtained the optimum tuning ratio in 1932, and later the optimum damping ratio was proposed by Brock in 1946 [4]. In 2002, Nishihara and Asami [5] presented the exact series solution for $H_{\infty}$ optimization criterion and compared it with the result by Den Hartog, and found that they are very close in control performance. The $H_{2}$ optimization criterion was first proposed by Crandall and Mark [6] in 1963 when they considered the primary system was subjected to random excitation. The object of this criterion is to minimize the total energy of the primary system over the whole-frequency range. In other words, this object is to minimize the area under the amplitude-frequency curve of the primary system. In the late 2000s Iwata [7] and Asami [8] also studied the optimum parameters respectively based on this optimization criterion. The stability maximization criterion was presented by Yamaguchi [9] in 1988 to stabilize the transient vibration of the primary system as quick as possible. All the three optimization criteria have been analytically solved for the undamped primary system. However, for the damped primary system it is difficult to obtain the analytical solution. Accordingly, numerical method was usually dopted to solve this problem. For example, Ioi and Ikeda [10], Randll [11], and Thompson [12] studied the optimum parameters and presented some important results by numerical methods. A detailed numerical study was conducted by Warburton [13] for the primary system with light damping, and he designed a table about the optimum parameters. In 1997 Nishihara and Matsuhisa [14] derived the optimum parameters for the primary system with light damping according to the stability maximization criterion. After that Asami and Nishihara [15, 16] presented the series solutions for the $H_{\infty}$ and $H_{2}$ optimization criteria.

Fractional-order calculus was first proposed by Hospital and Leibniz in the late 1700s. In the early development, the research was focused on the definition, properties and computation methods of the fractional-order calculus. Due to the lack of physical and mechanical meaning, it was slowly developed only as a mathematical branch [17-19]. Until Mandelbrot [20] proposed the fractal theory, the study on fractional-order calculus, especially the study on fractional-order differential equation started to develop rapidly. The research on fractional-order dynamical systems could be mainly divided into three groups, i.e. the qualitatively analysis, numerical study, and analytical research on the fractional-order system. The qualitatively analysis is focused on the number and stability of periodic solutions. Machado and Galhano [21], Li et al. [22], Wang and Hu [23], and Wang and Du [24], Rossikhin and Shitikova [25] had studied the composition of the solution for some fractional-order differential equations, and obtained some important results on its stability and properties. The numerical study is focused on the numerically investigation about the complicated nonlinear phenomenon in fractional-order dynamical system, such as bifurcation and chaos. Cao et al. [26], Sheu et al. [27] investigated the effect of the fractional-order parameters on the different nonlinear systems. The analytical research is focused on the approximate solution and the quantitative analysis of fractional-order differential equation. Wahi and Chatterjee [28], Shen and Yang [29-31], Chen and Zhu [32-35] studied the effects of fractional-order parameters on the dynamical response, and presented some important results to improve the control performance by approximately selecting fractional-order parameters.

Although the fractional-order has many advantages in engineering practice, little study on fractional-order DVA has been fulfilled. In engineering, the viscoelastic material could be modelled by fractional-order derivative, such as the air spring, metal rubber and magnetorheological elastomer. In this paper, fractional-order Voigt type DVA is introduced, and the optimum parameters of the presented DVA is studied in detail. In Section 2 the optimum fractional-order parameters, i.e. the fractional order and coefficient are obtained based on the $H_{\infty}$ and $H_{2}$ optimization criteria. The research shows that the linear damping and stiffness in the traditional Voigt type DVA can be replaced by a single fractional-order element completely in the two optimization cases. The comparisons of the approximately analytical solutions with the numerical results are fulfilled in Section 3, and the comparisons between the fractional-order and the traditional integer-order Voigt type DVA are also given in this section. Due to the advantages of fractional-order derivative, this research provides a theoretical basis to use only one fractional-order element to replace the linear spring and damper simultaneously.

2. Analytically investigation on fractional-order DVA

The model of fractional-order Voigt type DVA is shown in Fig. 1, which consists of the primary system $m_{1}$ and a damped DVA $m_{2}$ . The fractional-order element is put between the primary system and the DVA, which is in parallel with the spring and damper of the DVA. In vibration engineering, the fractional-order element may be arbitrary viscoelastic material and its force could be modelled as:

1a

F_{1} = K_{1} D^{P} (x_{1} - x_{2}),

where $x_{1}$ and $x_{2}$ are the displacement of the primary system and the DVA, $D^{P} (x_{1} - x_{2})$ is the $p$ -order derivative of $x_{1} - x_{2}$ to $t$ with the fractional coefficient $K_{1}$ ( $K_{1} \geq 0$ ) and the fractional order $p$ ( $0 \leq p \leq 1$ ). There are several definitions for fractional-order derivative, and under some conditions they are equivalent. Here the Caputo’s definition is adopted for simplicity:

1b

D^{p} [x (t)] = \frac{1}{Γ (1 - p)} \int_{0}^{t} \frac{x^{'} (u)}{(t - u)^{p}} d u,

where $Γ (z)$ is the Gamma function satisfying $Γ (z + 1) = z Γ (z)$ . In the fractional-order derivative, the initial condition is very important [36-38]. However, the steady-state responses are more meaningful in vibration engineering and both types of the derivatives will lead to the same steady-state solutions. Accordingly, the objective will be focused on the steady-state solution in the rest parts of this paper.

According to Newtonian second law, the motion equation could be established as:

2

\{\begin{array}{l} m_{1} {\ddot{x}}_{1} + c_{1} {\dot{x}}_{1} + k_{1} x_{1} + c_{2} ({\dot{x}}_{1} - {\dot{x}}_{2}) + k_{2} (x_{1} - x_{2}) + K_{1} D^{P} (x_{1} - x_{2}) = F_{0} c o s (ω t), \\ m_{2} {\ddot{x}}_{2} - c_{2} ({\dot{x}}_{1} - {\dot{x}}_{2}) - k_{2} (x_{1} - x_{2}) - K_{1} D^{P} (x_{1} - x_{2}) = 0, \end{array}

where $m_{1}$ , $m_{2}$ , $k_{1}$ , $k_{2}$ , $c_{1}$ , $c_{2}$ are the masses, linear stiffness coefficients, and linear viscous damping coefficients of the primary system and DVA respectively. $F_{0}$ and $ω$ are the amplitude and frequency of the external excitation.

Using the following parametric transformation:

y_{1} = x_{1}, y_{2} = x_{1} - x_{2}, c_{1} = 2 ζ_{1} ω_{1} m_{1}, c_{2} = 2 ζ_{2} ω_{2} m_{2}

ω_{1}^{2} = \frac{k_{1}}{m_{1}}, ω_{2}^{2} = \frac{k_{2}}{m_{2}}, μ = \frac{m_{2}}{m_{1}}, K = \frac{K_{1}}{m_{2}}, F = \frac{F_{0}}{m_{1}} .

Eq. (2) becomes:

3

\{\begin{array}{l} {\ddot{y}}_{1} + 2 ζ_{1} ω_{1} {\dot{y}}_{1} + ω_{1}^{2} y_{1} + 2 μ ζ_{2} ω_{2} {\dot{y}}_{2} + μ ω_{2}^{2} y_{2} + μ K D^{p} (y_{2}) = F c o s (ω t), \\ {\ddot{y}}_{1} - {\ddot{y}}_{2} - 2 ζ_{2} ω_{2} {\dot{y}}_{2} - ω_{2}^{2} y_{2} - K D^{p} (y_{2}) = 0 . \end{array}

Fig. 1The model of fractional-order DVA

$The model of fractional-order DVA$

2.1. The analytical solution based on the averaging method

According to the averaging method, one can suppose Eq. (3) has the solution as:

4

\{\begin{array}{l} y_{1} = a_{1} c o s φ_{1}, \\ {\dot{y}}_{1} = - ω a_{1} s i n φ_{1}, \end{array} \{\begin{array}{l} y_{2} = a_{2} c o s φ_{2}, \\ {\dot{y}}_{2} = - ω a_{2} s i n φ_{2}, \end{array}

where $φ_{1} = ω t + θ_{1}$ , $φ_{2} = ω t + θ_{2}$ are the generalized phases. By differentiating Eq. (4) to $t$ one can obtain:

5

\{\begin{array}{l} {\dot{y}}_{1} = {\dot{a}}_{1} c o s φ_{1} - a_{1} (ω + {\dot{θ}}_{1}) s i n φ_{1}, \\ {\ddot{y}}_{1} = - ω {\dot{a}}_{1} s i n φ_{1} - ω a_{1} (ω + {\dot{θ}}_{1}) c o s φ_{1}, \end{array}

\{\begin{array}{l} {\dot{y}}_{2} = {\dot{a}}_{2} c o s φ_{2} - a_{2} (ω + {\dot{θ}}_{2}) s i n φ_{2}, \\ {\ddot{y}}_{2} = - ω {\dot{a}}_{2} s i n φ_{2} - ω a_{2} (ω + {\dot{θ}}_{2}) c o s φ_{2} . \end{array}

Combining Eq. (3) and Eq. (4) with Eq. (5), one can get:

6

\{\begin{array}{l} {\dot{a}}_{1} c o s φ_{1} - a_{1} {\dot{θ}}_{1} s i n φ_{1} = 0, \\ {\dot{a}}_{1} s i n φ_{1} + a_{1} {\dot{θ}}_{1} c o s φ_{1} = - \frac{Δ_{1} - μ K D^{p} (a_{2} c o s φ_{2})}{ω}, \end{array}

\{\begin{array}{l} {\dot{a}}_{2} c o s φ_{2} - a_{2} {\dot{θ}}_{2} s i n φ_{2} = 0, \\ {\dot{a}}_{2} s i n φ_{2} + a_{2} {\dot{θ}}_{2} c o s φ_{2} = - \frac{Δ_{2} - (1 + μ) K D^{p} (a_{2} c o s φ_{2})}{ω}, \end{array}

where:

Δ_{1} = F c o s (ω t) + 2 ζ_{1} ω_{1} ω a_{1} s i n φ_{1} - ω_{1}^{2} a_{1} c o s φ_{1} + 2 μ ζ_{2} ω_{2} ω a_{2} s i n φ_{2}

- μ ω_{2}^{2} a_{2} c o s φ_{2} + a_{1} ω^{2} c o s φ_{1},

Δ_{2} = F c o s (ω t) + 2 ζ_{1} ω_{1} ω a_{1} s i n φ_{1} - ω_{1}^{2} a_{1} c o s φ_{1} + 2 (1 + μ) ζ_{2} ω_{2} ω a_{2} s i n φ_{2}

- (1 + μ) ω_{2}^{2} a_{2} c o s φ_{2} + a_{2} ω^{2} c o s φ_{2} .

Solving Eq. (6) with ${\dot{a}}_{1}$ , ${\dot{θ}}_{1}$ , ${\dot{a}}_{2}$ , ${\dot{θ}}_{2}$ as unknowns, one obtains:

7

\{\begin{array}{l} {\dot{a}}_{1} = - \frac{Δ_{1} - μ K D^{p} (a_{2} c o s φ_{2})}{ω} s i n φ_{1}, \\ a_{1} {\dot{θ}}_{1} = - \frac{Δ_{1} - μ K D^{p} (a_{2} c o s φ_{2})}{ω} c o s φ_{1}, \end{array}

\{\begin{array}{l} {\dot{a}}_{2} = - \frac{Δ_{2} - (1 + μ) K D^{p} (a_{2} c o s φ_{2})}{ω} s i n φ_{2}, \\ a_{2} {\dot{θ}}_{2} = - \frac{Δ_{2} - (1 + μ) K D^{p} (a_{2} c o s φ_{2})}{ω} c o s φ_{2} . \end{array}

Moreover, one could apply the standard averaging method to Eq. (7) in time interval $T$ [39-40], that means:

8

\{\begin{array}{l} {\dot{a}}_{1} = \frac{1}{T} \int_{0}^{T} [- \frac{Δ_{1} - μ K D^{p} (a_{2} c o s φ_{2})}{ω} s i n φ_{1}] d φ_{1}, \\ a_{1} {\dot{θ}}_{1} = \frac{1}{T} \int_{0}^{T} [- \frac{Δ_{1} - μ K D^{p} (a_{2} c o s φ_{2})}{ω} c o s φ_{1}] d φ_{1}, \end{array}

\{\begin{array}{l} {\dot{a}}_{2} = \frac{1}{T} \int_{0}^{T} [- \frac{Δ_{2} - (1 + μ) K D^{p} (a_{2} c o s φ_{2})}{ω} s i n φ_{2}] d φ_{2}, \\ a_{2} {\dot{θ}}_{2} = \frac{1}{T} \int_{0}^{T} [- \frac{Δ_{2} - (1 + μ) K D^{p} (a_{2} c o s φ_{2})}{ω} c o s φ_{2}] d φ_{2} . \end{array}

One could select the time terminal $T$ as $T = 2 π$ if the integrands are periodic or $T = \infty$ if the integrands are aperiodic. Then Eq. (8) can be divided into the following types:

9

\{\begin{array}{l} {\dot{a}}_{1} = {\dot{a}}_{11} + {\dot{a}}_{12}, \\ a_{1} {\dot{θ}}_{1} = a_{1} {\dot{θ}}_{11} + a_{1} {\dot{θ}}_{12}, \end{array} \{\begin{array}{l} {\dot{a}}_{2} = {\dot{a}}_{21} + {\dot{a}}_{22}, \\ a_{2} {\dot{θ}}_{2} = a_{2} {\dot{θ}}_{21} + a_{2} {\dot{θ}}_{22}, \end{array}

where:

10a

\{\begin{array}{l} {\dot{a}}_{11} = - \frac{1}{2 π} \int_{0}^{2 π} \frac{Δ_{1} s i n φ_{1}}{ω} d φ_{1}, \\ a_{1} {\dot{θ}}_{11} = - \frac{1}{2 π} \int_{0}^{2 π} \frac{Δ_{1} c o s φ_{1}}{ω} d φ_{1}, \end{array}

10b

\{\begin{array}{l} {\dot{a}}_{12} = \underset{T \to \infty}{l i m} \frac{1}{T} \int_{0}^{T} [\frac{μ K D^{p} (a_{2} c o s φ_{2})}{ω} s i n φ_{1}] d φ_{1}, \\ a_{1} {\dot{θ}}_{12} = \underset{T \to \infty}{l i m} \frac{1}{T} \int_{0}^{T} [\frac{μ K D^{p} (a_{2} c o s φ_{2})}{ω} c o s φ_{1}] d φ_{1}, \end{array}

10c

\{\begin{array}{l} {\dot{a}}_{21} = - \frac{1}{2 π} \int_{0}^{2 π} \frac{Δ_{2} s i n φ_{2}}{ω} d φ_{2}, \\ a_{2} {\dot{θ}}_{21} = - \frac{1}{2 π} \int_{0}^{T} \frac{Δ_{2} c o s φ_{2}}{ω} d φ_{2}, \end{array}

10d

\{\begin{array}{l} {\dot{a}}_{22} = \underset{T \to \infty}{l i m} \frac{1}{T} \int_{0}^{2 π} [\frac{(1 + μ) K D^{p} (a_{2} c o s φ_{2})}{ω} s i n φ_{2}] d φ_{2}, \\ a_{2} {\dot{θ}}_{22} = \underset{T \to \infty}{l i m} \frac{1}{T} \int_{0}^{T} [\frac{(1 + μ) K D^{p} (a_{2} c o s φ_{2})}{ω} c o s φ_{2}] d φ_{2} . \end{array}

Then the first part of Eq. (9) will become:

11

\{\begin{array}{l} {\dot{a}}_{11} = - \frac{1}{2 ω} \{F s i n θ_{1} + 2 ω a_{1} ζ_{1} ω_{1} + μ a_{2} ω_{2} [2 ω ζ_{2} \cos (θ_{1} - θ_{2}) - ω_{2} \sin (θ_{1} - θ_{2})]\}, \\ a_{1} {\dot{θ}}_{11} = - \frac{1}{2 ω} \{F c o s θ_{1} + a_{1} (ω^{2} - ω_{1}^{2}) - μ a_{2} ω_{2} [2 ω ζ_{2} s i n (θ_{1} - θ_{2}) + ω_{2} c o s (θ_{1} - θ_{2})]\}, \end{array}

\{\begin{array}{l} {\dot{a}}_{21} = - \frac{1}{2 ω} \{\begin{matrix} F s i n θ_{2} + a_{1} ω_{1} [2 ω ζ_{1} \cos (θ_{1} - θ_{2}) + ω_{1} \sin (θ_{1} - θ_{2})] \\ + 2 (1 + μ) ω a_{2} ζ_{2} ω_{2} \end{matrix}\}, \\ a_{2} {\dot{θ}}_{21} = - \frac{1}{2 ω} \{\begin{matrix} F c o s θ_{2} + a_{1} ω_{1} [2 ω ζ_{1} \sin (θ_{1} - θ_{2}) - ω_{1} \cos (θ_{1} - θ_{2})] \\ + a_{2} [ω^{2} - (1 + μ) ω_{2}^{2}] \end{matrix}\} . \end{array}

For the second part of Eq. (9), the integration region should be $[0, \infty]$ . Here $a_{12}$ is taken as an example, and the others could be solved similarly. In order to calculate this integration, two important formulae are introduced [29-31]:

12

\underset{T \to \infty}{l i m} \int_{0}^{T} \frac{s i n (ω t)}{t^{p}} d t = ω^{p - 1} Γ (1 - p) c o s (\frac{p π}{2}),

\underset{T \to \infty}{l i m} \int_{0}^{T} \frac{c o s (ω t)}{t^{p}} d t = ω^{p - 1} Γ (1 - p) s i n (\frac{p π}{2}) .

According to Eq. (12), one can get the integral in Eq. (10b):

13a

{\dot{a}}_{12} = \frac{μ K a_{2} ω^{p - 1}}{2} [s i n (θ_{1} - θ_{2}) c o s (\frac{p π}{2}) - c o s (θ_{1} - θ_{2}) s i n (\frac{p π}{2})] .

The detail derivation procedures are presented in Appendix A. Similarly, one can get:

13b

a_{1} {\dot{θ}}_{12} = \frac{μ K a_{2} ω^{p - 1}}{2} [s i n (θ_{1} - θ_{2}) s i n (\frac{p π}{2}) + c o s (θ_{1} - θ_{2}) c o s (\frac{p π}{2})],

13c

{\dot{a}}_{22} = - \frac{(1 + μ) K a_{2} ω^{p - 1}}{2} s i n (\frac{p π}{2}),

13d

a_{2} {\dot{θ}}_{22} = \frac{(1 + μ) K a_{2} ω^{p - 1}}{2} c o s (\frac{p π}{2}) .

Combining Eq. (13) with Eq. (11), one could obtain:

14

\{\begin{array}{l} {\dot{a}}_{1} = - \frac{1}{2 ω} \{F s i n θ_{1} + 2 ω a_{1} ζ_{1} ω_{1} + μ a_{2} ω_{2} [2 ω ζ_{2} c o s (θ_{1} - θ_{2}) - ω_{2} s i n (θ_{1} - θ_{2})]\} \\ + \frac{μ K a_{2} ω^{p - 1}}{2} [s i n (θ_{1} - θ_{2}) c o s (\frac{p π}{2}) - c o s (θ_{1} - θ_{2}) s i n (\frac{p π}{2})], \\ a_{1} {\dot{θ}}_{1} = - \frac{1}{2 ω} \{F c o s θ_{1} + a_{1} (ω^{2} - ω_{1}^{2}) - μ a_{2} ω_{2} [2 ω ζ_{2} s i n (θ_{1} - θ_{2}) + ω_{2} c o s (θ_{1} - θ_{2})]\} \\ + \frac{μ K a_{2} ω^{p - 1}}{2} [s i n (θ_{1} - θ_{2}) s i n (\frac{p π}{2}) + c o s (θ_{1} - θ_{2}) c o s (\frac{p π}{2})], \end{array}

\{\begin{array}{l} {\dot{a}}_{2} = - \frac{1}{2 ω} \{F s i n θ_{2} + a_{1} ω_{1} [2 ω ζ_{1} c o s (θ_{1} - θ_{2}) + ω_{1} s i n (θ_{1} - θ_{2})] + 2 (1 + μ) ω a_{2} ζ_{2} ω_{2}\} \\ - \frac{(1 + μ) K a_{2} ω^{p - 1}}{2} s i n (\frac{p π}{2}), \\ a_{2} {\dot{θ}}_{21} = - \frac{1}{2 ω} \{\begin{matrix} F c o s θ_{2} - a_{1} ω_{1} [- 2 ω ζ_{1} \sin (θ_{1} - θ_{2}) + ω_{1} \cos (θ_{1} - θ_{2})] \\ + a_{2} [ω^{2} - (1 + μ) ω_{2}^{2}] \end{matrix}\} \\ + \frac{(1 + μ) K a_{2} ω^{p - 1}}{2} c o s (\frac{p π}{2}) . \end{array}

The steady-state system response is more meaningful. Letting the right part of Eq. (14) be zero and simplifying the equation, one can obtain:

15

\{\begin{array}{l} 0 = - \frac{1}{2 ω} \{F s i n {\bar{θ}}_{1} + 2 ω {\bar{a}}_{1} ζ_{1} ω_{1} + μ {\bar{a}}_{2} [K ω^{p} s i n (\frac{p π}{2} - {\bar{θ}}_{1} + {\bar{θ}}_{2}) \\ + 2 ω ζ_{2} ω_{2} c o s ({\bar{θ}}_{1} - {\bar{θ}}_{2}) - ω_{2}^{2} s i n ({\bar{θ}}_{1} - {\bar{θ}}_{2})]}, \\ 0 = \frac{1}{2 ω} \{- F c o s {\bar{θ}}_{1} + {\bar{a}}_{1} (- ω^{2} + ω_{1}^{2}) + μ {\bar{a}}_{2} [K ω^{p} c o s (\frac{p π}{2} - {\bar{θ}}_{1} + {\bar{θ}}_{2}) \\ + 2 ω ζ_{2} ω_{2} s i n ({\bar{θ}}_{1} - {\bar{θ}}_{2}) + ω_{2}^{2} c o s ({\bar{θ}}_{1} - {\bar{θ}}_{2})]}, \\ 0 = - \frac{1}{2 ω} {F s i n {\bar{θ}}_{2} + {\bar{a}}_{1} ω_{1} [2 ω ζ_{1} c o s ({\bar{θ}}_{1} - {\bar{θ}}_{2}) + ω_{1} s i n ({\bar{θ}}_{1} - {\bar{θ}}_{2})] \\ + (1 + μ) {\bar{a}}_{2} [K ω^{p} s i n (\frac{p π}{2}) + 2 ω ζ_{2} ω_{2}]\}, \\ 0 = \frac{1}{2 ω} {- F c o s {\bar{θ}}_{2} + {\bar{a}}_{1} ω_{1} [- 2 ω ζ_{1} s i n ({\bar{θ}}_{1} - {\bar{θ}}_{2}) + ω_{1} c o s ({\bar{θ}}_{1} - {\bar{θ}}_{2})] \\ + {\bar{a}}_{2} [- ω^{2} + K (1 + μ) ω^{p} c o s (\frac{p π}{2}) + (1 + μ) ω_{2}^{2}]\} . \end{array}

The steady-state amplitudes ${\bar{a}}_{1}$ , ${\bar{a}}_{2}$ and phases ${\bar{θ}}_{1}$ , ${\bar{θ}}_{2}$ can be obtained by solving Eq. (15). Accordingly, one can get:

16

{\bar{a}}_{1} = \frac{F \sqrt{{[K ω^{p} s i n (\frac{p π}{2}) + 2 ω ζ_{2} ω_{2}]}^{2} + {[K ω^{p} c o s (\frac{p π}{2}) + ω_{2}^{2} - ω^{2}]}^{2}}}{\sqrt{Δ_{3}}}, {\bar{a}}_{2} = \frac{F ω^{2}}{\sqrt{Δ_{3}}},

where:

Δ_{3} = {\{[(1 + μ) ω^{2} - ω_{1}^{2}] [K ω^{p} s i n (\frac{p π}{2}) + 2 ω ζ_{2} ω_{2}] + 2 ω ζ_{1} ω_{1} [ω^{2} - K ω^{p} c o s (\frac{p π}{2}) - ω_{2}^{2}]\}}^{2}

+ \{- 2 ω ζ_{1} ω_{1} [K ω^{p} s i n (\frac{p π}{2}) + 2 ω ζ_{2} ω_{2}] + ω_{1}^{2} [- ω^{2} + K ω^{p} c o s (\frac{p π}{2}) + ω_{2}^{2}]

+ {ω^{2} [ω^{2} - K (1 + μ) ω^{p} c o s (\frac{p π}{2}) - (1 + μ) ω_{2}^{2}]\}}^{2} .

Substituting the parameters with the original ones, Eq. (16) could be transformed into:

17

{\bar{a}}_{1} = \frac{F \sqrt{{[K_{1} ω^{p} s i n (\frac{p π}{2}) + ω c_{2}]}^{2} + {[K_{1} ω^{P} c o s (\frac{p π}{2}) + k_{2} - ω^{2} m_{2}]}^{2}}}{\sqrt{Δ_{4}}}, {\bar{a}}_{2} = \frac{F ω^{2} m_{2}}{\sqrt{Δ_{4}}},

where:

Δ_{4} = {\{ω c_{1} [K_{1} ω^{p} c o s (\frac{p π}{2}) + k_{2} - ω^{2} m_{2}] + [K_{1} ω^{p} s i n (\frac{p π}{2}) + ω c_{2}] [k_{1} - ω^{2} (m_{1} + m_{2}]\}}^{2}

+ \{ω c_{1} [K_{1} ω^{p} s i n (\frac{p π}{2}) + ω c_{2}] - k_{1} [K_{1} ω^{p} c o s (\frac{p π}{2}) + k_{2} - ω^{2} m_{2}]

+ {ω^{2} m_{2} [K_{1} ω^{p} c o s (\frac{p π}{2}) + k_{2}] + ω^{2} m_{1} [K_{1} ω^{p} c o s (\frac{p π}{2}) + k_{2} - ω^{2} m_{2}]\}}^{2} .

Two new parameters $C_{2} (p)$ and $K_{2} (p)$ , defined as equivalent damping coefﬁcient and the equivalent stiffness coefﬁcient, can be introduced. Then Eq. (17) can be rewritten as:

18

{\bar{a}}_{1} = \frac{F_{0} \sqrt{ω^{2} C_{2} (p)^{2} + {[K_{2} (p) - ω^{2} m_{2}]}^{2}}}{\sqrt{Δ_{5}}}, {\bar{a}}_{2} = \frac{F_{0} ω^{2} m_{2}}{\sqrt{Δ_{5}}},

where:

C_{2} (p) = c_{2} + K_{1} ω^{p - 1} s i n (\frac{p π}{2}), K_{2} (p) = k_{2} + K_{1} ω^{p} c o s (\frac{p π}{2}),

Δ_{5} = {\{ω c_{1} [K_{2} (p) - ω^{2} m_{2}] + ω C_{2} (p) [k_{1} - ω^{2} (m_{1} + m_{2})]\}}^{2}

{+ \{ω^{2} c_{1} C_{2} (p) - k_{1} [K_{2} (p) - ω^{2} m_{2}] + ω^{2} m_{2} K_{2} (p) + ω^{2} m_{1} [K_{2} (p) - ω^{2} m_{2}]\}}^{2} .

2.2. $H_{\infty}$ optimization of the DVA

The calculation will be very complicated if the primary system contains viscous damping, so that the primary system without damping is considered. One could suppose $ω \approx ω_{1}$ in the equivalent damping coefﬁcient and the equivalent stiffness coefﬁcient because the main object of $H_{\infty}$ optimization criterion is to reduce the maximum resonance amplitude of the primary system. Here the amplitude amplification factor $A$ is defined based on Eq. (18) and the parameters in Eq. (3):

19

A = \frac{{\bar{a}}_{1}}{(F_{0} / k_{1})}

= \sqrt{\frac{(2 \bar{ν} λ {\bar{ζ}}_{2})^{2} + ({\bar{ν}}^{2} - λ^{2})^{2}}{4 λ^{2} [ζ_{1} ({\bar{ν}}^{2} - λ^{2}) + \bar{ν} {\bar{ζ}}_{2} (1 - λ^{2} - μ λ^{2})]^{2} + [{\bar{ν}}^{2} - (1 + 4 \bar{ν} ζ_{1} {\bar{ζ}}_{2} + {\bar{ν}}^{2} + μ {\bar{ν}}^{2}) λ^{2} + λ^{4}]^{2}}},

where:

{\bar{ζ}}_{2} = \frac{C_{2} (p)}{2 m_{2} Ω_{2}}, \bar{ν} = \frac{Ω_{2}}{ω_{1}}, λ = \frac{ω}{ω_{1}}, Ω_{2} = \sqrt{\frac{K_{2} (p)}{m_{2}}} .

Considering the undamped primary system, namely $c_{1} = 0$ , Eq. (19) can be:

20

\bar{A} = \frac{{\bar{a}}_{1}}{(F_{0} / k_{1})} = \sqrt{\frac{(2 \bar{ν} λ {\bar{ζ}}_{2})^{2} + ({\bar{ν}}^{2} - λ^{2})^{2}}{(2 λ \bar{ν} {\bar{ζ}}_{2})^{2} (1 - λ^{2} - μ λ^{2})^{2} + [{\bar{ν}}^{2} - (1 + {\bar{ν}}^{2} + μ {\bar{ν}}^{2}) λ^{2} + λ^{4}]^{2}}} .

According to the fixed-point theory, the amplitude of the primary system will pass through two fixed points independent of the damping ratio ${\bar{ζ}}_{2}$ . Assuming the two $\bar{A}$ to be equal when ${\bar{ζ}}_{2} = 0$ and ${\bar{ζ}}_{2} \to \infty$ , one can get:

21

\bar{A} = \frac{{\bar{ν}}^{2} - λ^{2}}{{\bar{ν}}^{2} - (1 + {\bar{ν}}^{2} + μ {\bar{ν}}^{2}) λ^{2} + λ^{4}} = \frac{\pm 1}{1 - λ^{2} - μ λ^{2}} .

One can find that there is no meaning if the right part of the equal sign is positive. Taking the negative one and simplifying the equation, one could obtain:

22

λ^{4} - \frac{2 λ^{2} (1 + {\bar{ν}}^{2} + μ {\bar{ν}}^{2})}{2 + μ} + \frac{2 {\bar{ν}}^{2}}{2 + μ} = 0 .

Solving Eq. (22), one can get:

23

λ_{1,2}^{2} = \frac{1 + (1 + μ) {\bar{ν}}^{2}}{2 + μ} \pm \frac{\sqrt{1 - 2 {\bar{ν}}^{2} + (1 + μ)^{2} {\bar{ν}}^{4}}}{2 + μ} .

2.2.1. The optimal equivalent frequency ratio

Substituting $λ_{1,2}^{2}$ into $\bar{A}$ respectively when ${\bar{ζ}}_{2} \to \infty$ , one can obtain:

24

\frac{1}{1 - {λ_{1}}^{2} - μ {λ_{1}}^{2}} = - \frac{1}{1 - {λ_{2}}^{2} - μ {λ_{2}}^{2}} .

Simplifying Eq. (24), one will get:

25

λ_{1}^{2} + λ_{2}^{2} = \frac{2}{1 + μ} .

According to Eq. (23), one can obtain:

26

λ_{1}^{2} + λ_{2}^{2} = \frac{2 (1 + {\bar{ν}}^{2} + μ {\bar{ν}}^{2})}{2 + μ} .

Comparing Eq. (25) with Eq. (26), one can get:

27

\frac{2 (1 + {\bar{ν}}^{2} + μ {\bar{ν}}^{2})}{2 + μ} = \frac{2}{1 + μ} .

The optimal equivalent frequency ratio can be obtained by solving Eq. (27):

28

\bar{ν} = \frac{1}{1 + μ} .

Substituting $\bar{ν}$ into $λ_{1,2}^{2}$ , then the abscissas of the two fixed points can be obtained:

29

λ_{1,2}^{2} = \frac{1}{1 + μ} (1 \pm \sqrt{\frac{μ}{2 + μ}}) .

2.2.2. The optimal equivalent damping ratio

Substituting the optimal equivalent frequency ratio $\bar{ν}$ into $\bar{A}$ , differentiating $\bar{A}$ to $t$ , and equating the obtained slopes to zero at the two optimal equivalent frequency ratio, one could get the best equivalent damping ratio:

30

({\bar{ζ}}_{2})_{1,2} = \sqrt{\frac{μ}{8 (1 + μ)} (3 \pm \sqrt{\frac{μ}{2 + μ}})} .

Based on Eq. (30), one could obtain an average value between the two given optimal equivalent damping ratios:

31

{\bar{ζ}}_{2} = \sqrt{\frac{({\bar{ζ}}_{2})_{1}^{2} + ({\bar{ζ}}_{2})_{2}^{2}}{2}} = \sqrt{\frac{3 μ}{8 (1 + μ)}} .

In order to determine the optimal equivalent stiffness and equivalent damping, the optimal frequency ratio and the optimal damping ratio are replaced by the equivalent stiffness and the equivalent damping. After simplification one could get:

32

\{\begin{array}{l} C_{2} (p) = \frac{2 m_{2} ω_{1}}{1 + μ} \sqrt{\frac{3 μ}{8 (1 + μ)}} = c_{2} + K_{1} ω_{1}^{p - 1} s i n (\frac{p π}{2}), \\ K_{2} (p) = \frac{m_{2} ω_{1}^{2}}{(1 + μ)^{2}} = k_{2} + K_{1} ω_{1}^{p} c o s (\frac{p π}{2}) . \end{array}

The optimal fractional order $p$ and the fractional coefficient $K_{1}$ can be obtained by solving Eq. (32):

33

\{\begin{array}{l} p = \frac{2}{π} a r c t a n [\frac{m_{2} ω_{_{1}}^{2} \sqrt{6 μ (1 + μ)} - c_{2} ω_{1} (1 + μ)^{2}}{2 m_{2} ω_{1}^{2} - 2 k_{2} (1 + μ)^{2}}], \\ K_{1} = \frac{\sqrt{{[k_{2} - \frac{m_{2} ω_{1}^{2}}{(1 + μ)^{2}}]}^{2} + {[c_{2} ω_{1} - \frac{m_{2} ω_{1}^{2} \sqrt{6 μ}}{2 (1 + μ)^{3 / 2}}]}^{2}}}{ω_{1}^{p}} . \end{array}

That means, the fractional-order element could replace the linear spring and damper simultaneously in $H_{\infty}$ optimization criterion, and reduce the maximum resonance amplitude of the primary system the same as traditional Voigt type DVA.

2.3. $H_{2}$ optimization of the DVA

$H_{2}$ optimization will be more desirable if the system is subjected to random excitation instead of harmonic excitation, because the object of this optimization criterion is to reduce the total vibration energy of the primary system over the whole-frequency range. It means the area under the amplitude-frequency curve of the primary system should be minimized.

The performance index is defined as follow for the $H_{2}$ optimization criterion:

34

I_{1} = \frac{E [x_{1}^{2}]}{2 π S_{f} ω_{1} / k_{1}^{2}} = \frac{⟨x_{_{1}}^{2}⟩}{2 π S_{f} ω_{1} / k_{1}^{2}},

where the symbol $E [\cdot]$ is statistical average, the symbol $⟨\cdot⟩$ is the temporal averages, and $S_{f}$ is the uniform power spectrum density of the excitation force respectively. The mean square value of the displacement of the primary system is calculated by the following equation:

35

⟨x_{1}^{2}⟩ = \frac{S_{f}}{k_{1}^{2}} \int_{- \infty}^{+ \infty} A^{2} d ω = \frac{S_{f} ω_{1}}{k_{1}^{2}} \int_{- \infty}^{+ \infty} A^{2} d λ,

where $A$ is the amplitude amplification factor. Thus the performance index can be simplified as:

36

I_{1} = \frac{1}{2 π} \int_{- \infty}^{+ \infty} A^{2} d λ .

Based on the residue theorem, Eq. (36) will be:

37

I_{1} = \frac{p_{1} {\bar{ζ}}_{2}^{3} + p_{2} {\bar{ζ}}_{2}^{2} + p_{3} {\bar{ζ}}_{2}^{} + p_{4}}{4 (q_{1} {\bar{ζ}}_{2}^{3} + q_{2} {\bar{ζ}}_{2}^{2} + q_{3} {\bar{ζ}}_{2}^{} + q_{4})},

where:

p_{1} = 4 (1 + μ) {\bar{ν}}^{2}, p_{2} = 4 \bar{ν} ζ_{1} [1 + (1 + μ) {\bar{ν}}^{2}],

p_{3} = 1 + (1 + μ)^{2} {\bar{ν}}^{4} - {\bar{ν}}^{2} (2 + μ - 4 ζ_{1}^{2}), p_{4} = μ {\bar{ν}}^{3} ζ_{1},

q_{1} = 4 {\bar{ν}}^{2} ζ_{1} (1 + μ), q_{2} = \bar{ν} \{μ + 4 ζ_{1}^{2} [1 + {\bar{ν}}^{2} (1 + μ)]\},

q_{3} = ζ_{1} [1 + {\bar{ν}}^{4} (1 + μ)^{2} - 2 {\bar{ν}}^{2} (1 - 2 ζ_{1}^{2})], q_{4} = μ {\bar{ν}}^{3} ζ_{1}^{2} .

The performance index has a minimum value at a certain combination of ${\bar{ζ}}_{2}$ and $\bar{ν}$ when the value of $μ$ and $ζ_{1}$ is given. This condition is achieved when both the partial derivatives of $I_{1}$ with respect to ${\bar{ζ}}_{2}$ and $\bar{ν}$ are zero. Considering the undamped primary system, one can get a pair of simultaneous equations after differentiation:

38

\{\begin{array}{l} {\bar{ζ}}_{2}^{4} (4 μ {\bar{ν}}^{2} + 4 {\bar{ν}}^{2}) + {\bar{ζ}}_{2}^{2} (3 μ^{2} {\bar{ν}}^{4} + 6 μ {\bar{ν}}^{4} - μ {\bar{ν}}^{2} + 3 {\bar{ν}}^{4} - 2 {\bar{ν}}^{2} - 1) = 0, \\ {\bar{ζ}}_{2}^{4} (4 μ {\bar{ν}}^{2} + 4 {\bar{ν}}^{2}) + {\bar{ζ}}_{2}^{2} (- μ^{2} {\bar{ν}}^{4} - 2 μ {\bar{ν}}^{4} + μ {\bar{ν}}^{2} - {\bar{ν}}^{4} + 2 {\bar{ν}}^{2} - 1) = 0 . \end{array}

Eliminating ${\bar{ζ}}_{2}$ for solving the functions, the equation about $\bar{ν}$ can be obtained:

39

2 {(1 + μ)}^{2} {\bar{ν}}^{4} - (2 + μ) {\bar{ν}}^{2} = 0 .

The optimal equivalent tuning ratio can be obtained by solving Eq. (39):

40

{\bar{ν}}_{o p t} = \frac{1}{1 + μ} \sqrt{\frac{2 + μ}{2}} .

Substituting ${\bar{ν}}_{o p t}$ into Eq. (38), the optimal equivalent damping ratio ${\bar{ζ}}_{2}$ can be obtained as:

41

{\bar{ζ}}_{2}_{o p t} = \sqrt{\frac{μ (4 + 3 μ)}{8 (1 + μ) (2 + μ)}} .

Replacing the parameters with the equivalent stiffness and the equivalent damping, one could get a pair of equations about $p$ and $K_{1}$ :

42

\{\begin{array}{l} C_{2} (p) = \frac{1}{2} m_{2} ω_{1} \sqrt{\frac{μ (4 + 3 μ)}{(1 + μ)^{3}}} = c_{2} + K_{1} ω_{1}^{p - 1} s i n (\frac{p π}{2}), \\ K_{2} (p) = \frac{(2 + μ) m_{2} ω_{1}^{2}}{2 {(1 + μ)}^{2}} = k_{2} + K_{1} ω_{^{1}}^{p} c o s (\frac{p π}{2}) . \end{array}

Solving the equations one can get the optimal fractional order and the fractional coefficient for $H_{2}$ optimization criterion:

43

\{\begin{array}{l} p = \frac{2}{π} a r c t a n [\frac{2 ω_{1} c_{2} (1 + μ)^{2} - m_{2} ω_{1} (1 + μ) \sqrt{μ (4 + 3 μ) / (1 + μ)}}{2 k_{2} (1 + μ)^{2} - m_{2} ω_{1}^{2} (2 + μ)}], \\ K_{1} = \frac{\sqrt{{[\frac{1}{2} m_{2} ω_{1}^{2} \sqrt{\frac{μ (4 + 3 μ)}{(1 + μ)^{3}}} - c_{2} ω_{1}]}^{2} + {[\frac{(2 + μ) m_{2} ω_{1}^{2}}{2 {(1 + μ)}^{2}} - k_{2}]}^{2}}}{ω_{1}^{p}} . \end{array}

That means, the fractional-order element could replace the linear spring and damper simultaneously in $H_{2}$ optimization criterion, and reduce the total energy of the primary system as traditional Voigt type DVA.

3. Numerical simulation and comparison

3.1. The comparison between the analytical and numerical solution

The numerical results are also presented to verify the precision of the analytical solution. Here the numerical formula presented in [17, 18] are adopted:

44

D^{p} [x (t_{l})] \approx h^{- p} \sum_{j = 0}^{l} C_{j}^{p} x (t_{l - j}),

where $t_{l} = l h$ is the sample points, $h$ is the sample step, $C_{j}^{p}$ is the binomial coefficient with the iterative relationship as:

45

C_{0}^{p} = 1, C_{j}^{p} = (1 - \frac{1 + p}{j}) C_{j - 1}^{p} .

The total computation time is selected as 400 s and the sample step is 0.001 s. The peak value of the later 50 s response is taken as the steady-state amplitude of the numerical results and the temporary response in frontal 350 s is omitted.

Based on the above optimized results, one could investigate the effect of the fractional-order parameters through the following four groups of parameters:

a) $k_{2} =$ 0, $c_{2} = 0$ .

b) $k_{2} =$ 0, $c_{2} = 100$ ,

c) $k_{2} =$ 6000, $c_{2} = 0$ ,

d) $k_{2} =$ 6000, $c_{2} = 100$ .

The other system parameters are selected as $m_{1} =$ 1000, $k_{1} =$ 100000, $c_{1} =$ 0, $m_{2} =$ 100, $F_{0} =$ 1000000.

3.1.1. $H_{\infty}$ optimization criterion

According to Eq. (19) and Eq. (33), the analytical normalized amplitude-frequency curves of the primary system are plotted in Fig. 2 and denoted by the solid lines. The corresponding numerical results are also shown in Fig. 2, and denoted by the circles. From the figure it could be found that the two kinds of curves agree well, that means the approximately analytical solution is satisfactorily precise. Moreover, the similarity of those curves shows that all the analytical amplitudes remain almost the same no matter how the stiffness coefficient $k_{2}$ and damping coefficient $c_{2}$ changes. That means the fractional order and fractional coefficient can affect the dynamic system characteristics by affecting the equivalent stiffness and the equivalent damping coefficient. Even it could replace the linear spring and damper of DVA completely. This result is different from the traditional results where the fractional-order derivative was generally regarded as damping device only.

3.2. $H_{2}$ optimization criterion

Based on the aforementioned system parameters, the normalized amplitude-frequency curves are plotted in Fig. 3 according to the optimal $p$ and $K_{1}$ in Eq. (43). It could also be found that the analytical solutions agree very well with the numerical results in $H_{2}$ optimization case.

Fig. 2The normalized frequency response curve based on H∞ optimization criterion

a) $k_{2} =$ 0, $c_{2} =$ 0

b) $k_{2} =$ 0, $c_{2} =$ 100

c) $k_{2} =$ 6000, $c_{2} =$ 0

d) $k_{2} =$ 6000, $c_{2} =$ 100

Fig. 3The normalized frequency response curve based on H2 optimization criterion

a) $k_{2} =$ 0, $c_{2} =$ 0

b) $k_{2} =$ 0, $c_{2} =$ 100

c) $k_{2} =$ 6000, $c_{2} =$ 0

d) $k_{2} =$ 6000, $c_{2} =$ 100

To study more realistic situation, 50 s random excitation are constructed, which is composed of 5000 random numbers with normal distribution. The excitation is normalized so that it is with mean value 0 and variance 1. The time history of the random excitation is shown in Fig. 4. When there is no fractional-order derivative in the Voigt type DVA, the time history of the primary system is given in Fig. 5. Here the optimum parameters are $k_{2} =$ 8686 and $c_{2} =$ 285, which are obtained according to reference [15]. For comparison the time histories of the primary system with the fractional-order DVA for the above different system parameters are presented in Fig. 6 to Fig. 9. The variances of the primary system in the above five systems are summarized in Table 1.

Table 1The variances of the displacement of the primary system

	The integer-order DVA	The fractional-order DVA
	The integer-order DVA	$k_{2} =$ 0, $c_{2} =$ 0	$k_{2} =$ 0, $c_{2} =$ 100	$k_{2} =$ 6000, $c_{2} =$ 0	$k_{2} =$ 6000, $c_{2} =$ 100
Variances	3.929e-11	2.8434e-11	2.8632e-11	2.8348e-11	2.8577e-11

Fig. 4The time history of random excitation

From the Figs. 4-9 and Table 1 it could be concluded that the fractional-order DVA can not only reduce the peak displacement of the primary system ( $H_{\infty}$ optimization criterion), but also reduce the total energy of the primary system in the whole-frequency range ( $H_{2}$ optimization criterion).

Fig. 5The time history of the primary system without fractional-order derivative (k2= 8686, c2= 285)

$The time history of the primary system without fractional-order derivative (k2= 8686, c2= 285)$

Fig. 6The time history of the primary system with fractional-order derivative (k2= 0, c2= 0)

$The time history of the primary system with fractional-order derivative (k2= 0, c2= 0)$

Fig. 7The time history of the primary system with fractional-order derivative (k2= 0, c2= 100)

$The time history of the primary system with fractional-order derivative (k2= 0, c2= 100)$

Fig. 8The time history of the primary system with fractional-order derivative (k2= 6000, c2= 0)

$The time history of the primary system with fractional-order derivative (k2= 6000, c2= 0)$

Fig. 9The time history of the primary system with fractional-order derivative (k2= 6000, c2= 100)

$The time history of the primary system with fractional-order derivative (k2= 6000, c2= 100)$

3.3. The comparison with the Voigt type DVA

In order to verify the performance of the fractional-order DVA, the comparison between the fractional-order and the traditional integer-order Voigt type DVA is presented in this subsection.

The amplification factor of the integer-order Voigt type DVA is:

46

A = \sqrt{\frac{(2 ν λ ζ_{2})^{2} + (ν^{2} - λ^{2})^{2}}{(2 λ ν ζ_{2})^{2} (1 - λ^{2} - μ λ^{2})^{2} + [ν^{2} - (1 + ν^{2} + μ ν^{2}) λ^{2} + λ^{4}]^{2}}} .

According to the reference [4], the optimal parameters of Voigt DVA for the $H_{\infty}$ optimization criterion is $μ = 0.1$ , $ν_{o p t} = 0.909$ , $ζ_{2 o p t} = 0.185$ .

The comparison of the fractional-order and the traditional integer-order Voigt type DVA is shown in Fig. 10. The parameters $k_{2}$ and $c_{2}$ are all selected as 0 in the fractional-order DVA so as to find out the effect of fractional-order derivative on the system clearly.

Similarly, one can choose the following parameters according to the existing literature [15] for the $H_{2}$ optimization criterion $μ = 0.1$ , $ν_{o p t} = 0.932$ , $ζ_{2 o p t} = 0.153$ .

Fig. 10The comparison of the fractional-order and integer-order Voigt type DVA for H∞ optimization criterion

$The comparison of the fractional-order and integer-order Voigt type DVA for H∞ optimization criterion$

Fig. 11The comparison of the fractional-order and integer-order Voigt type DVA for H2 optimization criterion

$The comparison of the fractional-order and integer-order Voigt type DVA for H2 optimization criterion$

To observe the effect of the fractional-order derivative, one can also select $k_{2}$ and $c_{2}$ to be zero. The comparison results are shown in Fig. 11. From Fig. 10 and Fig. 11, it could be found that in the two cases the amplitude-frequency curves of the fractional-order DVA agree well with those of the integer-order DVA, even when the linear stiffness and damping coefficients of the fractional-order DVA equal zero. That means the fractional-order derivative has a very good substitution for the stiffness and damping elements of DVA. The results also prove the correctness of the concept of equivalent stiffness and damping coefficients presented in references [29-31].

4. Conclusion

The optimization of the fractional-order DVA is studied in this paper, and the optimal fractional order and fractional coefficient are obtained for the $H_{\infty}$ and $H_{2}$ optimization criteria. The comparisons of the approximately analytical solutions with the numerical solutions certify the satisfactory precision of the approximately analytical solutions. The comparisons with the integer-order DVA show that the fractional-order derivative can almost replace the linear stiffness and linear damping element, which may make the dynamical system much simpler. As the fractional-order derivative can accurately model many viscoelastic materials for vibration control engineering, the linear spring and damping elements may be replaced by one fractional-order element. These results may be very useful to design or revise vibration control device in engineering. Of course, it is important to study the replacement in engineering practice, which may be the next study object.

References

Frahm H. Device for Damping Vibrations of Bodies. U.S. Patent, No. 989958, 30 October 1909.

Search CrossRef
Ormondroyd J., Den Hartog J. P. The theory of the dynamic vibration absorber. ASME Journal of Applied Mechanics, Vol. 1928, Issue 50, 7, p. 9-22.

Search CrossRef
Hahnkamm E. The damping of the foundation vibrations at varying excitation frequency. Ingenieur Archiv, Vol. 1932, 4, p. 192-201.

Publisher
Brock J. E. A note on the damped vibration absorber. ASME Journal of Applied Mechanics, Vol. 13, Issue 4, 1946, p. 284.

Publisher
Nishihara O., Asami T. Close-form solutions to the exact optimizations of dynamic vibration absorber (minimizations of the maximum amplitude magnification factors). ASME Journal of Vibration and Acoustics, Vol. 124, 2002, p. 576-582.

Publisher
Crandall S. H., Mark W. D. Random Vibration in Mechanical Systems. Academic Press, New York, 1963.

Search CrossRef
Iwata Y. On the construction of the dynamic vibration absorbers. Japan Society of Mechanical Engineers, Vol. 820, Issue 8, 1982, p. 150-152.

Search CrossRef
Asami T. Optimum design of dynamic absorbers for a system subjected to random excitation. JSME International Journal, Vol. 34, Issue 2, 1991, p. 218-226.

Publisher
Yamaguchi H. Damping of transient vibration by a dynamic absorber. Transactions of the Japan Society of Mechanical Engineers Series C, Vol. 54, Issue 499, 1988, p. 561-568.

Publisher
Ioi T., Ikeda K. On the dynamic vibration damped absorber of the vibration system. Bulletin of The Japan Society of Mechanical Engineers, Vol. 21, Issue 151, 1978, p. 64-71.

Publisher
Randall S. E., Halsted D. M., Taylor D. L. Optimum vibration absorbers for linear damped systems. ASME Journal of Mechanical Design, Vol. 103, 1981, p. 908-913.

Publisher
Thompson A. G. Optimum tuning and damping of a dynamic vibration absorber applied to a force excited and damped primary system. Journal of Sound and Vibration, Vol. 77, Issue 3, 1981, p. 403-415.

Publisher
Warburton G. B. Optimal absorber parameters for various combinations of response and excitation parameters. Earthquake Engineering and Structural Dynamics, Vol. 10, 1982, p. 381-401.

Publisher
Nishihara O., Matsuhisa H. Design and tuning of vibration control devices via stability criterion. Japan Society Mechanical Engineering, Vol. 97, Issue 10, 1997, p. 165-168.

Search CrossRef
Asami T., Nishihara O., Baz A. M. Analytical solutions to H_∞ and H₂ optimization of dynamic vibration absorbers attached to damped linear systems. ASME Journal of Vibration and Acoustics, Vol. 124, 2002, p. 284-295.

Publisher
Asami T., Nishihara O. Closed-form exact solution to H_∞ optimization of dynamic vibration absorbers application to different transfer functions and damping systems. ASME Journal of Vibration and Acoustics, Vol. 125, Issue 3, 2003, p. 398-405.

Publisher
Oldham K. B., Spanier J. The Fractional Calculus Theory and Applications of Differentiation and Integration to Arbitrary Order. Academic Press, New York, 1974.

Search CrossRef
Podlubny I. Fractional Differential Equations. Academic Press, London, 1999.

Search CrossRef
Petras I. Fractional-Order Nonlinear System. Higher Education Press, Beijing, 2011.

Publisher
Mandelbrot B. B. Les Objets Fractals: Forme, Hasard et Dimension. Flammarion, Paris, 1975.

Search CrossRef
Machado J. A. T., Galhano A. M. S. Fractional dynamics: a statistical perspective. ASME Journal of Computational and Nonlinear Dynamics, Vol. 3, Issue 2, 2008, p. 021201

Publisher
Li G., Zhu Z., Cheng C. Dynamical stability of viscoelastic column with fractional derivative constitutive relation. Applied Mathematics and Mechanics, Vol. 22, Issue 3, 2001, p. 294-303.

Publisher
Wang Z., Hu H. Stability of a linear oscillator with damping force of fractional order derivative. Science China Physics, Mechanics and Astronomy, Vol. 53, Issue 2, 2010, p. 345-352.

Publisher
Wang Z., Du M. Asymptotical behavior of the solution of a SDOF linear fractionally damped vibration system. Shock and Vibration, Vol. 18, 2011, p. 257-268.

Publisher
Rossikhin Y. A., Shitikova M. V. Application of fractional derivatives to the analysis of damped vibrations of viscoelastic single mass systems. Acta Mechanica, Vol. 120, 1997, p. 109-125.

Publisher
Cao J., Ma C., Xie H., Jiang Z. Nonlinear dynamics of Dufﬁng system with fractional order damping. ASME Journal of Computational and Nonlinear Dynamics, Vol. 5, 2010, p. 041012

Publisher
Sheu L. J., Chen H. K., Chen J. H., Tam L. M. Chaotic dynamics of the fractionally damped Dufﬁng equation. Chaos, Solitons and Fractals Vol. 32, 2007, p. 1459-1468.

Publisher
Wahi P., Chatterjee A. Averaging oscillations with small fractional damping and delayed terms. Nonlinear Dynamics, Vol. 38, 2004, p. 3-22.

Publisher
Shen Y. J., Yang S. P., Xing H. J. Dynamical analysis of linear single degree-of-freedom oscillator with fractional-order derivative. Acta Physica Sinica, Vol. 61, Issue 11, 2012, p. 110505.

Publisher
Shen Y. J., Yang S. P., Xing H. J., Ma H. X. Primary resonance of Duffing oscillator with two kinds of fractional-order derivatives. International Journal of Non-Linear Mechanics, Vol. 47, 2012, p. 975-983.

Publisher
Shen Y. J., Yang S. P., Xing H. J., Gao G. S. Primary resonance of Dufﬁng oscillator with fractional-order derivative. Communications in Nonlinear Science and Numerical Simulation, Vol. 17, 2012, p. 3092-3100.

Publisher
Chen L., Hu F., Zhu W. Stochastic dynamics and fractional optimal control of quasi integrable Hamiltonian systems with fractional derivative damping. Fractional Calculus and Applied Analysis, Vol. 16, Issue 1, 2013, p. 189-225.

Publisher
Chen L. C., Zhu W. Q. Stochastic jump and bifurcation of Duffing oscillator with fractional derivative damping under combined harmonic and white noise excitations. International Journal of Non-linear Mechanics, Vol. 46, Issue 12, 2011, p. 1324-1329.

Publisher
Chen L. C., Li H. F., Li Z. S., Zhu W. Q. First passage failure of single-degree-of-freedom nonlinear oscillators with fractional derivative. Journal of Vibration and Control, Vol. 19, Issue 14, 2013, p. 2154-2163.

Publisher
Chen L. C., Wang W. H., Li Z. S., Zhu W. Q. Stationary response of Duffing oscillator with hardening stiffness and fractional derivative. International Journal of Non-Linear Mechanics, Vol. 48, Issue 1, 2013, p. 44-50.

Publisher
Achar N., Lorenzo C. F., Hartley T. T. Initialization and the Caputo Fractional Derivative. NASA John H. Glenn Research Center at Lewis Field Report, 2003.

Search CrossRef
Ortigueira M. D., Coito F. J. Initial conditions: what are we talking about? 3rd IFAC Workshop on Fractional Differentiation and its Applications, Ankara, Turkey, 2008.

Search CrossRef
Sabatier J., Merveillaut M., Malti R., Oustaloup A. On a representation of fractional order systems: interests for the initial condition problem. 3rd IFAC Workshop on Fractional Differentiation and its Applications, Ankara, Turkey, 2008.

Search CrossRef
Burd V. Method of Averaging for Differential Equations on an Infinite Interval: Theory and Applications. Taylor and Francis Group, 2007.

Publisher
Shen Y. J., Zhao Y. X., Yang S. P., Xing H. J. A revised averaging method and general forms of approximate solution for nonlinear oscillator with only polynomial-type displacement nonlinearity. Journal of Vibroengineering, Vol. 16, Issue 4, 2014, p. 1864-1876.

Search CrossRef

Cited by

Torsional vibration suppression for a high-damped rotor using dynamic vibration absorber based on weighted dual estimation of equivalent linearization techniques

Vu Duc Phuc | Phan Ngoc Tuan | Van – The Tran

(2024)

Optimal design and dynamic performance analysis of a fractional-order electrical network-based vehicle mechatronic ISD suspension

(2023)

Recent Advances in Manufacturing, Automation, Design and Energy Technologies

Anant J. Sheth | Utkarsh A. Patel

(2022)

The effect of fractional viscoelastic supports on the response of a flexible rotor based on H∞ and H2 optimization methods

Hamidreza Heidari | Habib Ramezannejad Azarboni | Pedram Safarpour

(2021)

The missing mem-inerter and extended mem-dashpot found

Xiao-Liang Zhang | Can Geng | Jia-Mei Nie | Qiao Gao

(2020)

Vibration performance improvement of D/G-set employing inerter-rubber vibration isolator

Huabing Wen | Yang Li | Kun Zhang | Yue Liu | Chengwei Chang

(2019)

About this article

Received

09 November 2015

Accepted

17 February 2016

Published

15 August 2016

SUBJECTS

Mechanical vibrations and applications

DOI

https://doi.org/10.21595/jve.2016.16617

Keywords

dynamic vibration absorber

fractional derivative

parameters optimization

averaging method

Acknowledgements

The authors are grateful to the support by National Natural Science Foundation of China (No. 11372198), the Cultivation Plan for Innovation Team and Leading Talent in Colleges and Universities of Hebei Province (LJRC018), the Program for Advanced Talent in the Universities of Hebei Province (GCC2014053), and the Program for Advanced Talent in Hebei Province (A201401001).

Author Contributions

Yongjun Shen presented the whole structure of the article. Yongjun Shen and Haibo Peng wrote the draft manuscript. Shaofang Wen, Shaopu Yang, and Haijun Xing helped to improve the quality of the article in the preparation procedure.

This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Next article in issue Next

Research article

2022 12 12

A position optimization strategy of the dynamic vibration absorber for the noise control of the power equipment

Mingxin Geng, Jiangang Ma, Yan Jing, Yalin Zhao, Lv Wang

Research article

2022 04 21

Optimal analysis on low frequency line spectrum sound control of power equipment using dynamic vibration absorber

Jiangang Ma, Jian Wu, Mingxin Geng, Bo Ma, Chen Shen

Research article

2019 05 15

Stochastic P-bifurcation in a tri-stable Van der Pol system with fractional derivative under Gaussian white noise

Yajie Li, Zhiqiang Wu

Research article

2019 03 31

Parameter optimization design of rotor dynamic vibration absorber

Zi Liang Liu, Qin Zhang, Yan Bo Cao, Qi Xu, Hong Liang Yao, Bang Chun Wen

Y. Shen, H. Peng, S. Wen, S. Yang, and H. Xing, “Analytically optimal parameters of fractional-order dynamic vibration absorber,” Journal of Vibroengineering, Vol. 18, No. 5, pp. 2714–2734, Aug. 2016, https://doi.org/10.21595/jve.2016.16617

Copy Extrica

Copied to clipboard!

TY  - JOUR
DO  - 10.21595/jve.2016.16617
UR  - https://doi.org/10.21595/jve.2016.16617
TI  - Analytically optimal parameters of fractional-order dynamic vibration absorber
T2  - Journal of Vibroengineering
AU  - Shen, Yongjun
AU  - Peng, Haibo
AU  - Wen, Shaofang
AU  - Yang, Shaopu
AU  - Xing, Haijun
PY  - 2016
DA  - 2016/08/15
PB  - JVE International Ltd.
SP  - 2714-2734
IS  - 5
VL  - 18
SN  - 1392-8716
ER  - 

Copy Ris

Copied to clipboard!

@article{Shen_2016,
	doi = {10.21595/jve.2016.16617},
	url = {https://doi.org/10.21595/jve.2016.16617},
	year = 2016,
	month = {aug},
	publisher = {{JVE} International Ltd.},
	volume = {18},
	number = {5},
	pages = {2714--2734},
	author = {Yongjun Shen and Haibo Peng and Shaofang Wen and Shaopu Yang and Haijun Xing},
	title = {Analytically optimal parameters of fractional-order dynamic vibration absorber},
	journal = {Journal of Vibroengineering}
}

Copy Bibtex

Copied to clipboard!

[1]Y. Shen, H. Peng, S. Wen, S. Yang, and H. Xing, “Analytically optimal parameters of fractional-order dynamic vibration absorber,” Journal of Vibroengineering, vol. 18, no. 5, pp. 2714–2734, Aug. 2016, doi: 10.21595/jve.2016.16617.

Copy IEEE

Copied to clipboard!

Shen, Yongjun, Haibo Peng, Shaofang Wen, Shaopu Yang, and Haijun Xing. “Analytically Optimal Parameters of Fractional-Order Dynamic Vibration Absorber.” Journal of Vibroengineering 18, no. 5 (August 15, 2016): 2714–34. https://doi.org/10.21595/jve.2016.16617.

Copy Chicago

Copied to clipboard!

Analytically optimal parameters of fractional-order dynamic vibration absorber

Abstract

1. Introduction

2. Analytically investigation on fractional-order DVA

2.1. The analytical solution based on the averaging method

2.2. H∞ optimization of the DVA

2.2.1. The optimal equivalent frequency ratio

2.2.2. The optimal equivalent damping ratio

2.3. H2 optimization of the DVA

3. Numerical simulation and comparison

3.1. The comparison between the analytical and numerical solution

3.1.1. H∞ optimization criterion

3.2. H2 optimization criterion

3.3. The comparison with the Voigt type DVA

4. Conclusion

References

Cited by

About this article

Related Articles

2.2. $H_{\infty}$ optimization of the DVA

2.3. $H_{2}$ optimization of the DVA

3.1.1. $H_{\infty}$ optimization criterion

3.2. $H_{2}$ optimization criterion