# Fractional Critical Damping Theory and Its Application in Active Suspension Control.

1. Introduction

The vibrations of linear 1 DOF systems with ordinary damping can be classified as underdamped, critically damped, and overdamped according to the magnitude of the damping coefficient. Critical damping is defined as the threshold between overdamping and underdamping. In the case of critical damping, the oscillator returns to the equilibrium position as quickly as possible, without oscillating, and passes it once at most [1]. Considering the particularity of critical damping, it is frequently studied in other systems. The criterion for critical damping of viscously damped multi-degree-of-freedom systems is provided by Bulatovic [2]. The existence conditions for the critical damping in second-order pendulum-like systems are established by Li et al. [3]. A general method that determines the "critical damping surfaces" of a certain linear continuous dynamic system is proposed by Beskos and Boley [4]. However, so far, there are only a few researches on the critical damping in fractionally damped systems. In 1984, Torvik and Bagley [5] proposed a mechanical model with fractional derivatives in the study of the motion of a rigid plate immersed in a Newton fluid, and the study results in [6, 7] make the fractional calculus attractive for many engineers and technicians [8].

Vehicle suspension is an important component for improving the driving comfort and the handling performance [9], the research on its control strategy is a hot spot. In these control approaches, the skyhook control strategy proposed by Karnopp et al. [10, 11] is widely applied because of its simple algorithm and good control performance. The classical skyhook control principle is based on a SDOF vibration system, which is suitable for the vertical vibration control of two DOFs quarter-car models. In recent years, many scholars have studied the application of skyhook algorithm in full-car suspension model. The mainstream skyhook control strategies for full-car suspension systems are based on physical thinking; these strategies are the application extensions of the classical skyhook method that is widely used to control the quarter-vehicle suspension systems. The full-car suspension model is regarded as a simple combination of four 1/4 subsuspension models, and it is assumed that there is a "skyhook" connected with each 1/4 car body by four skyhook dampers to control the vibration of the car body, whereas the full-vehicle suspension has requirements of multiobjective suspension performances involving the vertical, pitch, and roll motions [12]. Therefore, the problem of how to coordinate the forces of the four independent controllers to keep a good body posture [13] needs to be solved and the typical solution is adding a decision-making system.

Although mainstream algorithms can achieve a good control performance, it is inconsistent with the original skyhook control principle. From the perspective of mathematical principles, the classic skyhook control principle is used to control a SDOF system with one skyhook damper. Whereas the vehicle suspension is a system with multi-DOFs (the existing models have seven or more DOFs), thus the same number of controllers is required. However, in reality, there are only four controllers. How to tackle this problem?

This work is divided into two parts. In the first part, the critical damping in fractional order system is studied. The existence conditions of the critical damping are given, and the relationship order is derived. Then the vibration attenuation characteristics of fractional critical damping systems with different order are discussed. In the second part, the fractional critical damping is applied to the control strategy of the vehicle suspension system. The method of modal decoupling is used to solve the problem that the number of required controllers is not consistent with that of the actual ones. In the modal space, the classical skyhook control strategy is used for depressing the decoupled single mode vibration. Here, the fractional critical damping coefficients are chosen as the skyhook damping coefficients. In this way, the number of designed controllers is consistent with that of DOFs of the system, then these modes are recoupled and the actual controllers are used to control the suspension. A four-wheel-correlated random road time domain model is built to test the effect of fractional derivative skyhook control strategy; a road bump is especially designed to demonstrate the advantages of the fractional derivative critical damping.

The organization of the paper is as follow. In Section 2, the conditions of fractional damped systems being in critical damping case are given first. Then the properties of the vibration with critical damping are studied. In Section 3, a new fractional skyhook control algorithm for full-car suspension systems is proposed. In Section 4, the simulation results are discussed. Conclusions are given in Section 5.

2. Critical Damping of the System with Fractional Derivative Damping

2.1. Formula Derivation. The free vibration differential equation of a SDOF system with fractional derivative damping has the form

m[??] (t) + [c.sub.0][D.sup.[alpha]]x (t) + kx (t) = 0, (1)

where x(t) is the displacement, [sub.0][D.sup.[alpha]]x(t) is the fractional time derivative of x(t), and m, c, and k are the mass, damping, and stiffness coefficient, respectively.

There are many definitions for fractional derivatives [14], among which Riemann-Liouville definition and Caputo definitions are most widely used [15]. The former is frequently used for problem description because of its demand moderately for the continuity of the function. The latter has the same

Laplace transform as the integer order one, so it is widely used in control theory. In this paper, the fractional derivative damping force is regarded as a control force to study the properties of free damped vibration of the system, so the Caputo definition is used here.

By the Laplace transform method, the characteristic equation of the system takes the form

[ms.sup.2] + c[s.sup.[alpha]] + k = 0, (2)

where s is the complex variable. By substituting its polar form s = r[[theta].sup.ie] into (2), we have

[mr.sup.2] [e.sup.i(2[theta])] + [cr.sup.[alpha]] [e.sup.i([alpha][theta])] + k = 0, (3)

Considering the Euler formula [e.sup.i[theta]] = cos [theta] + i sin [theta], (3) takes the form

[mr.sup.2] (cos 2[theta] + i sin 2[theta]) + [cr.sup.[alpha]] (cos [alpha][theta] + i sin [alpha][theta]) + k = 0, (4)

The establishment-condition of (4) is that both the real and imaginary parts are equal to zero, so we obtain

[mr.sup.2] cos 2[theta] + [cr.sup.[alpha]] cos [alpha][theta] + k = 0

[mr.sup.2] sin 2[theta] + [cr.sup.[alpha]] sin [alpha][theta] = 0. (5)

It is known that when the imaginary part of the roots of (2) is non-zero, the damped free motion of the system is always oscillating. To avoid the oscillation, the characteristic roots must lie in the negative real axis. Assume that [theta] = (2[k.sub.1] + 1)n where [k.sub.1] is an integer, so cos 2[theta] = 1 and sin 2[theta] = 0 are obtained, then (5) can be simplified as

[mr.sup.2] + [cr.sup.[alpha]] cos [alpha][theta] + k = 0, (6)

[cr.sup.[alpha]] sin [alpha][theta] = 0. (7)

The establishing condition for (7) is sin [alpha][theta] = 0, which means that cos [alpha][theta] = [+ or -] 1. Therefore, it can be obtained that

[alpha][theta] = [alpha] (2[k.sub.1] + l)[pi] = [k.sub.2][pi], (8)

where [k.sub.2] is an integer. As a result, we have

[alpha] = [k.sub.2]/2[k.sub.1] + 1. (9)

We find that the set of a is dense, but the probability density of any [alpha] locating in this domain is small, so the existence condition of critical damping is strict.

From (6), a negative damping coefficient c is obtained when cos ad = 1, which represents an energy input to the system. In this case, the system oscillation is strengthened, and there is no critical damping, while it is the opposite when cos [alpha][theta] = -1; that is, [k.sub.2] is odd, so substituting cos [alpha][theta] = -1 into (6) and then (10) is obtained. In summary, in (9) [k.sub.1] is an integer, [k.sub.2] is odd, and [alpha] [member of] (0,2). The existence conditions of critical damping in the vibration systems with fractional derivative damping and its calculation formula are presented.

For linear 1 DOF fractionally damped systems, only when (9) is satisfied by the order of fractional operator, there is a critical value of damping coefficients. To make the solutions of (1) be without oscillation, the relation between the damping coefficient and the order is

c = [mr.sup.2] + k/[r.sup.[alpha]] (10)

where r [member of] [R.sup.+]. In (10), when dc/dr = 0, i.e. r = [square root of k[alpha]/m(2 - [alpha])], we have the minimum value of damping coefficient c which represents the critical damping coefficient cc.

The curves that represent the relation between the variables in (10) are plotted in Figure 1. Take [alpha] = 3/5, for example, the lowest point of the curve represents the critical damping point and its corresponding damping coefficient is the critical value of damping coefficient. It is worth noting that many previous researches on 1 DOF fractionally damped systems focus on the solutions of the characteristic equations. From this perspective, we find when c < [c.sub.c], the characteristic equations only have complex or conjugate roots and they have negative real roots when c [??] [c.sub.c]. Therefore, when c > [c.sub.c], it represents the overdamping coefficient, and when c < [c.sub.c], it is the underdamping coefficient. In the case of critical damping, the characteristic equation has the root s = -r, which represents the convergence rate. When a increases from 0 to 2, the critical damping point is shifted to the lower right in the figure, which indicates that with larger a, it turns out a smaller cc and larger r; that is, with a smaller eigenvalue, the system is a faster convergent.

It should be noted also that Sakakibara [16] studied the properties of vibration with fractional derivative damping of order 1/2. By the analysis of the solutions of (1), it is concluded that there is no critical value of damping coefficient, which is not against the conclusions of this paper because [alpha] = 1/2 is not located in the set represented by (9). In fact, it is easy to understand that by reduction to absurdity, that is, when the roots s are negative real, they do not hold by substituting [alpha] = 1/2 into (2). This means that when [alpha] = 1/2, the eigenvalues cannot be negative real and always contain an imaginary part. Furthermore, we find that when [alpha] = 1, the critical damping coefficients [c.sub.c] = 2[square root of mk], r = [square root of k/m], and s = -[square root of k/m] = -[[omega].sub.n] are obtained, which are consistent with the critical damping in an integer order system. Because it is not our main objective to solve the equation and the critical damping coefficients can be obtained without analyzing the solutions, we will not return to these questions here and refer the interested reader to [17, 18]. As is shown in Figure 2, when [alpha] = 9/11, the critical damping coefficient is figured out according to the above analysis.

2.2. Properties of Vibration with Fractional Derivative Critical Damping. When [alpha] [member of] (0,1), the fractional damping plays not only the role of a conventional damping, but also the role of a supplementary spring [19]. If [alpha] [right arrow] 0 or [alpha] [right arrow] 2, the damping effect of the system will be weakened, and there is a typical behavior of the oscillation. Furthermore, the fractional order systems are easily affected by the initial state. Therefore, in practice, [alpha] should lie within the range of engineering interest.

Figure 3 shows the curves of decaying free motions of critical damping systems with different orders under the initial state [x.sub.0] = 0.1, [x.sub.0] = 0. It shows that in the case of the same other parameters, the systems with a large a return back to equilibrium position faster. When [alpha] [member of] (0,1), the systems are relatively slow as they go back to balance position and do not cross it. Otherwise when [alpha] [member of] (1,2), the systems are relatively fast and cross through the static equilibrium position once (overshoot occurs), which is different from the ordinary critical damping. Although the systems with a large a return back to the equilibrium position at a faster speed, it is easily to be aroused by external excitation such as step input; the response curves are shown in Figure 4.

It is expected that under the premise of nonoscillatory, the system is not easy to be aroused by external excitation and can return back to the equilibrium position as quickly as possible when there is no external force. A switch control law is designed to make the displacement as small as possible when the system is away from the equilibrium position and to limit the time it takes to reach the asymptotically stable position when there is no external force. The designed control law is

[mathematical expression not reproducible], (11)

where u is the control force, [[alpha].sub.1] and [[alpha].sub.2] are the orders of fractional derivative, and [c.sub.1] and [c.sub.2] are the corresponding fractional derivative critical damping coefficients, x is the displacement. The effectiveness of the proposed control strategy is tested by a pulse excitation. Figure 5 shows that, under impulse input, the switching control law makes the vibration performance of the fractional order system better than that of the integer order one.

3. Vehicle Skyhook Control Strategy

According to vehicle dynamics theory, the dynamic model of the vehicle with seven DOFs is established. The seven DOFs [Z.sub.b], [theta], [Z.sub.wA], [Z.sub.wB], [Z.sub.wC], and [Z.sub.wD] are the heave, pitch, roll displacement of the body, and the four wheels displacement, respectively. This model is similar to those used by [20, 21], here the matrix differential equation of the model can be described as

[mathematical expression not reproducible], (12)

where X [member of] [R.sup.(7x1)] is a vector consisting of [Z.sub.b], [theta], [Z.sub.wA], [Z.sub.wB], [Z.sub.wC], and [Z.sub.wD]. M [member of] [R.sup.(7x7)], C [member of] [R.sub.(7x7)] and K [member of] [R.sup.(7x7)] are the mass, damping, and stiffness matrix, respectively. [K.sub.t] [member of] [R.sup.(7x4)] is the input matrix and [Z.sub.g(4x1)] is a vector that stands for the four-wheel-related road excitation. R [member of] [R.sup.(7x4)] is the control vector and [U.sub.(4x1)] is the vector of active control force. Equation (12) represents a passive suspension when [U.sub.(4x1)] is a zero vector.

According to the linear vibration theory, the decoupled suspension system turns into isolated linear subsystems that can be controlled independently [22]. Therefore, a systematic modal decoupling method [23] is considered, with which the mass and stiffness matrix can be completely decoupled; however, the damping matrix cannot be completely decoupled generally. Here only the diagonal elements of the damping matrix are controlled to verify the effectiveness of the control strategy. The matrix differential equation of the fully decoupled system is considered

[mathematical expression not reproducible], (13)

where

[mathematical expression not reproducible], (14)

[eta] [member of] [R.sup.(7x1)] is the vector of principal coordinates, [eta] = [Q.sup.-1] X, Q is the feature matrix, and [[bar.C].sub.Di] is the diagonal matrix whose diagonal elements are equal to those in vector [bar.C]. In (13), [bar.M], [bar.K], and [[bar.C].sub.Di] are seven-order diagonal matrices and, assuming that [bar.U] is also a seven-order diagonal matrices, seven differential equations of independent scalar function [[eta].sub.i] are obtained; fractional skyhook control is used here to depress each independent modal vibration. Let [u.sub.isky] = -[c.sub.isky][[eta].sub.i.sup.([alpha])] (i = 1 ~ 7), so the seven independent differential equations have the form

[mathematical expression not reproducible], (15)

where [F.sub.i] is the external excitation for modal vibration systems, [c.sub.isky][[eta].sup.([alpha])] represents the fractional skyhook damping force.

The free vibration equations of the modal systems are considered, namely,

[mathematical expression not reproducible], (16)

where the control force [c.sub.isky][[eta].sup.([alpha])] used to keep the system in the case of critical damping. According to the method in Section 2, the relation between the damping coefficient and the order is obtained

[mathematical expression not reproducible]. (17)

When [r.sub.i] = ([c.sub.i](1 - [alpha]) + [square root of [([c.sub.i](1 - [alpha])).sup.2] + 4[alpha](2 - a)[m.sub.i][k.sub.i])]/ 2[m.sub.i]i (2 - [alpha]), the fractional derivative skyhook damping coefficient [c.sub.ishk] is equal to the fractional derivative critical damping coefficient. In the same way, it is hoped that with the fractional damping force, the modal system is not easily aroused by external force and returns to equilibrium position as fast as possible without oscillating when there is no force. Here, a switching control law is given as follows:

[mathematical expression not reproducible]. (18)

In practice, with a larger or smaller [alpha], these problems, such as the limitation of actuator force and the work efficiency of the actuator, arise. In order to achieve a relatively good control effect, only the limitation of actuator force is considered. Seven skyhook damping coefficients [c.sub.ishk] of the system are obtained. By coordinate reduction, the final control force vector is

U =[([Q.sup.T]R).sup.-] x [Q.sup.-1] x [C.sub.shk] x [X.sup.([alpha])], (19)

where [C.sub.shk] = diag([[c.sub.ishk]]) (i = 1 ~ 7). Equation (19) represents the force of integer order skyhook damping control strategy when a =1. The generalized inverse matrix of [Q.sup.T]R is used here because it is not a square matrix.

4. Simulation Results and Discussions

A four-wheels-correlated random road time domain model [24] is used here and the road profile is C grade. To verify the characteristics of fractional critical damping, a work condition is designed as follows: when the simulation goes to 5s, on the left side of the vehicle, the front and rear wheels have been raised successively by road bump shaped like a sine wave with a height of 0.1m. Vehicle suspension parameters are shown in Notations. For validating the superiority of the fractional derivative critical damping, meanwhile avoiding the following negative effects with a large or small [alpha], in the switching control law, the orders are chosen as [[alpha].sub.1] = 9/11 and [[alpha].sub.2] = 13/11.

Figures 6 and 7 show that the proposed vehicle skyhook control strategy can effectively suppress the vibration of the body; both vibration amplitude and acceleration are decreased significantly; the performance especially is good after crossing the road bump. Figure 6 shows that the vibration with fractional derivative critical damping has a better performance on amplitude responses than that with integer one. And Figure 7 shows that fractional order skyhook damping control strategy has no significant deterioration in acceleration response. But for a large or small a, the acceleration responses become worse than those in integer order control strategy, and that is why the order should locate within a reasonable domain in engineering application.

Compared with many other full-car suspension control strategies, there are two main advantages for the method in this paper. Firstly, the proposed method is much more simple than most of the control methods. For example, these methods presented in [25] are also tested by a road bump and can improve the vibration performance of the vehicle, but they are too complicated. Actually, the skyhook control strategy is one of several simple and practical methods which is widely applied. Among the full-car skyhook control algorithms, a skyhook-based asynchronous semiactive controller proposed by Zhang et al. can control each subsystem independently; the results show that the peak amplitudes of body accelerations increase more than those in the passive suspension when they are tested by a pulse excitation. Therefore, it is not easy to keep a good body posture particularly when a car crosses through a road bump. The existing solution is to introduce new controls such as modularize parallel fuzzy control in [26] and human-like intelligent control in [27] and this makes the strategies complex and difficult to apply. Secondly, there are many active suspension control strategies which are designed based on a more comprehensive usage of road preview information facilitated by utilising on-board cameras and global positioning systems [28]. For example, it is provided that the road preview is available in the control method in [29]. However, our control method does not need such facilities.

In a word, the proposed skyhook control has a simple algorithm and is consistent with the original skyhook damping scheme in principle. The strategy with integer order critical damping coefficients has a good effect, and the fractional one is seen as a supplement, which provides more parameter selection and has a better performance on amplitude responses.

5. Conclusions

(1) The free damped motion of SDOF systems with fractional derivative damping is firstly studied. Conditions of existing critical damping are given and the relation between the critical damping coefficient and the order fractional derivative is derived. It is also found that when the order increases from 0 to 2, the critical damping coefficient is getting small, but it is faster to return back to equilibrium position.

(2) Based on the mathematical thinking, a new fullcar skyhook damping control strategy is proposed, which is different from the logical thinking of most scholars. The mainstream algorithm can also achieve a good performance; here, it is not the purpose to deny its effectiveness, but to give a new perspective for scholars to re-examine the intrinsic mathematical logic of classic skyhook damping principle. The fractional order critical damping coefficient is selected as the skyhook damping coefficient to clarify the superiority of proposed fractional order critical damping in practical application.

(3) Simulation results show that compared with the passive suspension, the skyhook controlled active suspension has a better performance on vibration suppression. Furthermore, the fractional skyhook controlled suspension has better responses of the body vibrating, especially when the vehicle passes the road bump. The results not only confirm the superiority of fractional critical damping, but also validate the effectiveness of this control strategy.
```Abbreviations

Vehicle Parameters

[M.sub.b]:       Sprung mass, 810 kg
[I.sub.p]:       Inertia moment of vehicle pitch, 300 kg x [m.sup.2]
[I.sub.r]:       Inertia moment of vehicle roll, 1058 kg x [m.sup.2]
A:               Distance from axle to 1.14 m
B:               Center of gravity, 1.22 m
[K.sub.s/1/2]:   Front suspension stiffness, 20600 N/m
[K.sub.sr1/2]:   Rear suspension stiffness, 15200 N/m
[C.sub.f1/2]:    Front suspension damping, 1570 N/m
[C.sub.r1/2]:    Rear suspension damping, 1760 N/m
[K.sub.w]:       Tire stiffness, 138000 N/m
[M.sub.wf]:      Front tire mass, 26.5 kg
[M.sub.wr]:      Rear tire mass, 24.4 kg
[B.sub.1]:       Distance between two tires, 1.3 m
V:               Vehicle speed, 50 km/h.
```

https://doi.org/10.1155/2017/2738976

Conflicts of Interest

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

Acknowledgments

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

References

[1] D. E. Beskos, "Critical Damping Surfaces of Linear Dynamic Systems," Engineering Mechanics: ASCE, pp. 722-725, 1979.

[2] R. M. Bulatovic, "On the critical damping in multi-degree-of-freedom systems," Mechanics Research Communications, vol. 29, no. 5, pp. 315-319, 2002.

[3] X. B. Li, Y. N. Huang, Y. Yang, and L. Huang, "Critical damping of the second-order pendulum-like systems," Applied Mathematics and Mechanics. Yingyong Shuxue He Lixue, vol. 26, no. 1, pp. 7-15, 2005.

[4] D. E. Beskos and B. A. Boley, "Critical damping in certain linear continuous dynamic systems," International Journal of Solids and Structures, vol. 17, no. 6, pp. 575-588, 1981.

[5] P. J. Torvik and R. L. Bagley, "On the appearance of the fractional derivative in the behavior of real materials," Journal of Applied Mechanics, vol. 51, no. 2, pp. 294-298, 1984.

[6] R. L. Bagley and P. J. Torvik, "Fractional calculus in the transient analysis of viscoelastically damped structures," AIAA Journal, vol. 23, no. 6, pp. 918-925, 1985.

[7] R. L. Bagley, "Power law and fractional calculus model of viscoelasticity," AIAA Journal, vol. 27, no. 10, pp. 1412-1417, 1989.

[8] Z. H. Wang and H. Y. Hu, "Stability of a linear oscillator with damping force of the fractional-order derivative," Science China Physics, Mechanics & Astronomy, vol. 53, no. 2, pp. 345-352, 2010.

[9] P. W. Claar II and J. M. Vogel, "A review of active suspension control for on and off-highway vehicles," SAE Technical Papers, 1989.

[10] D. Karnopp, M. J. Crosby, and R. A. Harwood, "Vibration control using semi-active force generators," Journal of Manufacturing Science and Engineering, vol. 96, no. 2, pp. 619-626, 1974.

[11] D. Karnopp, "Design principles for vibration control systems using semi-active dampers," Journal of Dynamic Systems, Measurement, and Control, vol. 112, no. 3, pp. 448-455, 1990.

[12] S. Ikenaga, F. L. Lewis, J. Campos, and L. Davis, "Active suspension control of ground vehicle based on a full-vehicle model," in Proceedings of the American Control Conference, vol. 6, pp. 4019-4024, IEEE, 2000.

[13] S. B. A. Kashem, S. Roy, and R. Mukharjee, "A modified skyhook control system (SKDT) to improve suspension control strategy of vehicles," in Proceedings of the International Conference on Informatics, Electronics and Vision, (ICIEV '14), pp. 1-8, IEEE, May 2014.

[14] S. Das, "Functional fractional calculus for system identification and controls," Functional Fractional Calculus for System Identification and Controls, pp. 1-240, 2008.

[15] J. A. Machado, I. S. Jesus, R. Barbosa, M. Silva, and C. Reis, "Application of fractional calculus in engineering," in Dynamics, games and science I, vol. 1 of Springer Proc. Math., pp. 619-629, Springer, Heidelberg, 2011.

[16] S. Sakakibara, "Properties of Vibration with Fractional Derivative Damping of Order 1/2," JSME International Journal Series C Mechanical Systems, Machine Elements and Manufacturing, vol. 40, no. 3, pp. 393-399, 1997.

[17] Y. A. Rossikhin and M. V. Shitikova, "Application of fractional calculus for dynamic problems of solid mechanics: novel trends and recent results," Applied Mechanics Reviews, vol. 63, no. 1, Article ID 010801, 2010.

[18] R. Lewandowski and M. Baum, "Dynamic characteristics of multilayered beams with viscoelastic layers described by the fractional Zener model," Archive of Applied Mechanics, vol. 85, no. 12, pp. 1793-1814, 2015.

[19] P. Wahi and A. Chatterjee, "Averaging oscillations with small fractional damping and delayed terms," Nonlinear Dynamics, vol. 38, no. 1-4, pp. 3-22, 2004.

[20] S. Ikenaga, F. L. Lewis, J. Campos, and L. Davis, "Active suspension control of ground vehicle based on a full-vehicle model," in Proceedings of the American Control Conference, pp. 4019-4024, 2010.

[21] R. Tchamna, M. Lee, and I. Youn, "Attitude control of full vehicle using variable stiffness suspension control," Optimal Control Applications and Methods, vol. 36, no. 6, pp. 936-952, 2015.

[22] T. K. Caughey and M. E. O'Kelly, "Classical normal modes in damped linear dynamic systems," Journal of Applied Mechanics, vol. 32, p. 583, 1965.

[23] C. Acar and A. M. Shkel, "Nonresonant micromachined gyroscopes with structural mode-decoupling," IEEE Sensors Journal, vol. 3, no. 4, pp. 497-506, 2003.

[24] Y. S. Wang, C.-M. Lee, and L. J. Zhang, "Wavelet analysis of vehicle nonstationary vibration under correlated four-wheel random excitation," International Journal of Automotive Technology, vol. 5, no. 4, pp. 257-268, 2004.

[25] L. Khan, S. Qamar, and M. U. Khan, "Comparative analysis of adaptive neurofuzzy control techniques for full car active suspension system," Arabian Journal for Science and Engineering, vol. 39, no. 3, pp. 2045-2069, 2014.

[26] S. Zhao and Y. Li, "Full vehicle semi-active suspension modularize parallel fuzzy control," Computer Simulation, vol. 5, no. 036, 2013.

[27] X. Dong M, M. Yu, W. Chen M et al., "Research on application of human-like intelligent control used to vehicle magneto rheological semi-active suspension," China Academic Journal Electronic Publishing House, vol. 18, no. 7, pp. 764-769, 2007.

[28] H. E. Tseng and D. Hrovat, "State of the art survey: active and semi-active suspension control," Vehicle System Dynamics, vol. 53, no. 7, pp. 1034-1062, 2015.

[29] M. Coric, J. Deur, L. Xu, H. E. Tseng, and D. Hrovat, "Optimisation of active suspension control inputs for improved vehicle ride performance," Vehicle System Dynamics, vol. 54, no. 7, pp. 1004-1030, 2016.

Peng Wang, (1) Qingyun Wang, (2) Xiaomei Xu, (2) and Ning Chen (1)

(1) College of Mechanical and Electronic Engineering, Nanjing Forestry University, Nanjing 210037, China

(2) College of Automobile and Traffic Engineering, Nanjing Forestry University, Nanjing 210037, China

Correspondence should be addressed to Ning Chen; chenning@njfu.com.cn

Received 9 June 2017; Revised 8 October 2017; Accepted 23 October 2017; Published 19 November 2017

Caption: FIGURE 1: The relation between the three parameters in (12).

Caption: FIGURE 2: When [alpha] = 9/11, the three cases of free decay vibration.

Caption: FIGURE 3: The curves of free damped motion of critical damping system with different orders.

Caption: FIGURE 4: Step response curves of critical damping systems with different orders.

Caption: FIGURE 5: Impulse response curves of integer order and fractional orders switched systems.

Caption: FIGURE 6: Body motion amplitude response of three kinds of suspension.

Caption: FIGURE 7: Body motion acceleration response of three kinds of suspension.