Design and Development of an Active Suspension System Using Pneumatic-Muscle Actuator and Intelligent Control

: A pneumatic muscle is a cheap, clean, and high-power active actuator. However, it is di ﬃ cult to control due to its inherent nonlinearity and time-varying characteristics. This paper presents a pneumatic muscle active suspension system (PM-ASS) for vehicles and uses an experimental study to analyze its stability and accuracy in terms of reducing vibration. In the PM-ASS, the pneumatic muscle actuator is designed in parallel with two MacPherson struts to provide a vertical force between the chassis and the wheel. This geometric arrangement allows the PM-ASS to produce the maximum force to counter road vibration and make the MacPherson struts generate signiﬁcant improvement. In terms of the controller design, this paper uses an adaptive Fourier neural network sliding-mode controller with H ∞ tracking performance for the PM-ASS, which confronts nonlinearities and time-varying characteristics. A state-predictor is used to predict the output error and to provide the predictions for the controller. Experiments with a rough concave-convex road and a two-bump excitation road use a quarter-car test rig to verify the practical feasibility of the PM-ASS, and the results show that the PM-ASS gives an improvement the ride comfort.


Introduction
Active suspension systems (ASS) [1][2][3][4] have been the subject of much development in the car industry over decades. To provide comfort ride and road holding for driving, ASS actively controls the vertical movement of the wheels when the road terrain changes. Therefore, ASS actively reduces road vibrations and gives a smoother ride experience for passengers than passive/conventional suspension systems because the design of a passive suspension system compromises between three conflicting criteria: road handling, road carrying, and passenger comfort [5,6].
To improve road handling and passenger comfort, ASS uses actuators to provide additional force between the sprung mass and the wheels to control the attitude of the vehicle. ASS also reduces vertical vibration in cameras that are mounted on a vehicle, thus, stable camera sequences are ensured for image processing [1,2]. Therefore, ASS is eminently suited to the needs of self-driving ground vehicles. The most-commonly used actuators for an ASS are pneumatic, linear electromagnetic, or hydraulic actuators. For a classical and conventional configuration, these actuators are commonly designed to be used in conjunction spring dampers to create an active suspension, which gives an additional vertical Figure 1 shows the design of the PM. It is made of natural rubber with fibers wrapped inside and (see Figure 1a) metal fittings attached at each end. The PM changes shape when the pressure of the air that flows into the interior of the rubber tube is changed. Injecting compressed air into the rubber tube expands the PM and produces a force that drags the sprung mass downward (see Figure 1b). Otherwise, the sprung mass returns to its original shape if the PM releases air. For simplicity, the reversible physical deformation during contraction and expansion produces a linear motion. According to the standard static physical model [22,23], (see Figure 2), the length and the diameter of the PM are formulated as: and: where L represents the actual length of PM, L 0 represents the original length of PM, l represents the length of the thread, n represents the number of turns in the thread, D represents the diameter, and θ represents the angle that the threads make with a horizontal axis.
Appl. Sci. 2019, 9, x FOR PEER REVIEW 3 of 22 desired H ∞ tracking performance for the overall system. It also significantly reduces the control chatter that is inherent in a conventional sliding mode control system. The main contributions of this paper are as follows: 1. This paper uses a pneumatic muscle actuator with two MacPherson struts to increase the isolation of vibration and show its feasibility on a quarter-car test rig. 2. A fully-functional test-rig for a quarter car is developed for examining the PM-ASS which has a road profile generator to provide real road conditions. 3. Using the orthogonal activation basis functions of Fourier series, Fourier-NN has a clear physical meaning and readily determined structure, which gives it an advantage over conventional NNs.
Because of the orthogonality of the basis functions, the Fourier-NN converges very fast, which means it is suited to real-time implementation. 4. A state-predictor M (1,1) G is used to predict the output error and to provide the predictions for the controller. 5. Four different control strategies are conducted to analyze the stability and accuracy in terms of reducing vibration, which are MacPherson Strut, FSASMC, FSB-ASMC+ H ∞ , and GP-FSB- Figure 1 shows the design of the PM. It is made of natural rubber with fibers wrapped inside and (see Figure 1a) metal fittings attached at each end. The PM changes shape when the pressure of the air that flows into the interior of the rubber tube is changed. Injecting compressed air into the rubber tube expands the PM and produces a force that drags the sprung mass downward (see Figure  1b). Otherwise, the sprung mass returns to its original shape if the PM releases air. For simplicity, the reversible physical deformation during contraction and expansion produces a linear motion. According to the standard static physical model [22,23], (see Figure 2), the length and the diameter of the PM are formulated as: (2)

Mathematical Model of the PM
where L represents the actual length of PM, 0 L represents the original length of PM, l represents the length of the thread, n represents the number of turns in the thread, D represents the diameter, and θ represents the angle that the threads make with a horizontal axis. Considering the PM as a cylindrical shape, its volume is defined as: Substituting Equations (1) and (2) into (3) yields: Appl. Sci. 2019, 9, x FOR PEER REVIEW 4 of 22 Considering the PM as a cylindrical shape, its volume is defined as: Substituting Equations (1) and (2) According to the Principle of the Conservation of Energy [24], the geometric force [22] is the product of the pressure and the change in the volume with respect to the length, which is: where P represents the pressure inside the PM and e P represents atmospheric pressure.
Substituting Equation (5) into (6) yields: A proportional pressure regulator (PPR) is used to regulate the pressure P inside the PM. This is modeled by a first order transfer function, which is: where v K and s T respectively represent the constant gain and the time constant for the PM, and u represents the input voltage. Substituting Equation (8) into (7) yields:  Because the active area of the cylinder changes over time, the volume of the PM changes with time. The active area is expressed as A = dV/dL, which is: According to the Principle of the Conservation of Energy [24], the geometric force [22] is the product of the pressure and the change in the volume with respect to the length, which is: where P represents the pressure inside the PM and P e represents atmospheric pressure. Substituting Equation (5) into (6) yields: A proportional pressure regulator (PPR) is used to regulate the pressure P inside the PM. This is modeled by a first order transfer function, which is: where K v and T s respectively represent the constant gain and the time constant for the PM, and u represents the input voltage. Substituting Equation (8) into (7) yields: Remark 1. In general, a PPR exhibits complex dynamic behavior and must be modeled by a high order nonlinear differential equation. Fortunately, for the PM-ASS, the bandwidth for a PPR is much wider than that for a PM; thus, the dynamics of a PPR can be modeled as a first order linear differential equation, as shown in Equation (8). Figure 3 shows an illustration of a quarter-car with a PM-ASS, in which the PM-ASS has two MacPherson struts and a PM, as shown in Figure 3a. For simplicity, the mathematical models of the two MacPherson struts are combined in one model, thus, they are expressed as one damper, one spring, and an internal frictional force. The dynamics of a tire are modeled as a spring if the tire fully contacts with the road surface. In this condition, the quarter-car is modeled as a system with 2-DOF (two-degrees-of-freedom). The frictional force of the piston in the shock absorber prevents the smooth travel of the MacPherson struts so the dynamic behavior of the PM-ASS must consider the frictional force. The dynamic equation for a quarter-car with a PM-ASS is:

Mathematical Model of a Quarter Car
and: where F µ is the frictional force, F a is the force that is exerted by the PM, Z r is the variation in the position of the road surface, Z s is the displacement of the sprung mass, Z u is the displacement of the un-sprung mass, M s is the weight of the rigid body (sprung mass), M u is the weight of the un-sprung mass, K t is a spring constant, and B s is a damping constant. If Z u , and x 5 = P, then manipulating Equations (8)-(11) with x 1 to x 5 yields a dynamic equation: Equation (12) can further be expressed as its matrix form, as given by: . where: Appl. Sci. 2019, 9, 4453 6 of 22 and: , the relative degree of the system is 3. Using input-output Feedback linearization, Equation (13) is transformed into a canonical form after simple manipulation, which is: where: and: The states are selected as ..
where the output z 1 is the displacement of the sprung mass, i.e., z 1 = x 1 = y, F(z) and G(z) cannot be assumed to be given continuous functions because they might contain unknown parameters, nonlinearities and uncertainties. It is noted that the functions F(z) and G(z) are unknown so it is impossible to determine a stable ideal control law for a quarter-car with a PM-ASS, as given by: Therefore, a Fourier neural network (FNN) is adapted to approximate the unknown functions F(z) and G(z); thus, a stable and robust control law can be derived for a quarter-car with a PM-ASS. Appl. Sci. 2019, 9, x FOR PEER REVIEW 7 of 22

Adaptive Fourier Neural Network Sliding-Mode Controller (FSB-ASMC) with H ∞ Tracking Performance
The PM-ASS features internal parametric uncertainties and external disturbances ( ) d t that arise from different unknown road terrains while driving; thus, the FSB-ASMC with H ∞ tracking performance was used for this study. While considering external disturbances ( ) d t , the mathematical model for the Quarter Car with PM-ASS (see Equation (19)) is changed to: is a vector of state s, which are assumed to be measurable; It is assumed that there exists a solution for (21) and the order of the nonlinear system (21) is known.

Sliding Mode Control
Sliding mode control can be used to reject uncertainties and external disturbance. Generally, it assumes that ( ) F z and ( ) G z are given. Here, we define a sliding mode control law for the quartercar with a PM-ASS while the system contains external disturbances. The output tracking error is written as: where ( ) d y t is a reference signal, which represents road conditions in this paper. A switch surface for an nth-order system is defined as: 1 2 s c e c e e = + +   (23) where i c are chosen such that Equation (23) shows that:

Adaptive Fourier Neural Network Sliding-Mode Controller (FSB-ASMC) with H ∞ Tracking Performance
The PM-ASS features internal parametric uncertainties and external disturbances d(t) that arise from different unknown road terrains while driving; thus, the FSB-ASMC with H ∞ tracking performance was used for this study. While considering external disturbances d(t), the mathematical model for the Quarter Car with PM-ASS (see Equation (19)) is changed to: where ..

y(t) ]
T ∈ R 3 is a vector of state s, which are assumed to be measurable; u(t) ∈ R and y(t) ∈ R are the control input and system output, respectively, d(t) is the bounded external disturbance, i.e., |d(t)| ≤ d u , F(z) and G(z) are smooth unknown functions with uncertain time-varying parameter. Without loss of generality, G(z) can be assumed to be strictly positive; i.e., G(z) ≥ G l > 0. It is assumed that there exists a solution for (21) and the order of the nonlinear system (21) is known.

Sliding Mode Control
Sliding mode control can be used to reject uncertainties and external disturbance. Generally, it assumes that F(z) and G(z) are given. Here, we define a sliding mode control law for the quarter-car with a PM-ASS while the system contains external disturbances. The output tracking error is written as: where y d (t) is a reference signal, which represents road conditions in this paper. A switch surface for an nth-order system is defined as: where c i are chosen such that 3 i=1 c i l i−1 is a Hurwitz polynomial in which l is a Laplace operator.
If e = [e, . e] T = [e 1 , e 2 ] T , the error dynamics are: . where Appl. Sci. 2019, 9, 4453 8 of 22 Lemma 1. For a nonlinear system (21) with the given nonlinear functions F(z) and G(z), if the control input is: and P 0, P ∈ R 2×2 satisfies the Lyapunov matrix equation: where s is the sliding surface defined in Equation (23), p 2i are elements of P, k > 0, and Q l 0 is given, then s → 0 and e → 0 as t → ∞ .
However, F(z) and G(z) for a quarter-car with a PM-ASS are generally uncertain rather than given, thus, there is no control input for Equation (26). The constant k in Equation (26) is designed to attenuate the uncertainties. However, a larger value for k produces increased chatter. To relax these two constraints, a Fourier Neural Network (Fourier-NN) is used to approximate the unknown functions F(z) and G(z) to eliminate model dependency. H ∞ tracking design is also used to compensate for the approximation error, which gives improved control performance and reduced chatter.

Description of Fourier Neural Networks
In this study, a Fourier-NN is used to approximate the unknown nonlinear functions with uncertainties, F(z) and G(z), in order to provide a stable control law for the PM-ASS. A Fourier-NN is eminently suited to the modeling and decomposition of a nonlinear function. It uses a Fourier series as the basic functions and allows much better convergence for the approximation. It is noted that the number of the hidden layers of the Fourier-NN depends on the bandwidth of the system. The weigh vector for a Fourier-NN is acquired used the Lyapunov stability theorem. The output function for a Fourier-NN is: where χ and p(χ), respectively, represent the input and the output for the Fourier-NN, q(x) ≡ cos(w 0 χ), cos(w 1 χ), sin(w 1 χ), . . . , cos(w M χ), sin(w M χ) T is a family of orthogonal Fourier T is a vector for network weights.

Remark 2.
The approximation error for a Fourier-NN is bounded and is expressed as: Clearly, the approximation error ε n (t) is eliminated as n → ∞ . This implies that f (χ) is approximated as long as n is sufficiently large.

Remark 3.
If a nonlinear function f (t) is a non-periodical function, which is written as: where F(w) = ∞ −∞ f (t)e −jwt dt, then by applying Shannon's theory in the frequency domain, Equation (30) is equivalent to the partition of a non-periodical function with an appropriate window in the time domain. This is written in discrete form as: where Ω = w n /n is the base frequency that is acquired by discretizing w. From an engineering point of view, the control objective is defined as a limited time interval of [0, T], thus, the base frequency is simply selected as ∆w = 2π/T.

Design of the FSB-ASMC with H ∞ Tracking Performance
Because the control voltage for the PM-ASS is specified and bounded, all of the states, which are real-world physical signals, are bounded. This implies that F(z(t)) and G(z(t) are bounded and meet the Dirichlet conditions. If time is the input for the FNN, i.e., χ = t in (28), and two Fourier-NNs are used to respectively approximate F(z(t)) and G(z(t) , Equation (21) is rewritten as: where q F (t) and q G (t) are the family of orthogonal Fourier activation function vector, W F and W G are the weight vectors, and q G produce an optimal mean-square approximation to the uncertain time-varying functions F(z(t)) and G(z(t)). (19) is assumed such that w ∈ L 2 [0, T], ∀T ∈ [0, ∞).

Assumption 2.
The internal dynamics of the system (19) that uses adaptive sliding-mode control law are stable.
Theorem 1. If a system can be expressed as a form of Equation (21) with unknown nonlinear time-varying functions F(z(t)) and G(z(t)) and has finite bandwidth. The sliding-mode control input is: whereŴ T F andŴ T G are the estimates of W T F and W T G , respectively. According to Lemma 1, P 0, P ∈ R 2×2 is selected to satisfy the Lyapunov matrix equation: where s is the sliding surface that is defined in Equation (23), p 2i are elements of P, σ is a designed constant, Q l 0 is given and the adaptive laws are: . and: .
where η 1 and η 2 (η 1 > 0 and η 2 > 0) are the adaptation gain matrices. Therefore, the H ∞ tracking performance for the overall system must satisfy: (see (37)) Remark 4. For a specified set of initial conditions, e(0) = 0, , and if the matrix Q l is set to an identity matrix, the control performance satisfies: T 0 e T (τ)e(τ)dτ. In other words, an arbitrary attenuation level is achieved if σ is properly chosen. Equation (38) also shows that e 2 ≤ σ w t 2 , thus, the bearing of the system (21) on the lumped uncertainty is stable. According to Assumption 2, the system (13) is stable. Figure 4 shows the overall scheme for a quarter-car with a PM-ASS that is controlled by the GP-FSB-ASMC+H ∞ , where the control law is shown in (33), the adaptive laws are Equations (34) and (35), and the sliding surface is Equation (23). The block "Road generator" simulates a reference road condition y m from a given terrain R r for the system. A Fuzzy logic controller (FLC) is used to track given terrain trajectories R r . To compensate for the delay in the control input, a "Grey Model First-Order One Variable" state-predictor GM(1, 1) predicts the future errorê for the FSB-ASMC+H ∞ , where GM(1, 1) is a time-series-forecasting model that predicts the possible error for the next sampling time. The steps for GM(1, 1) are:

Design Procedures
Step 1: Define the original data set: Step 2: Do Accumulating Generation Operator (AGO) to generate an accumulated data set: Step 3: Construct a differential equation for GM(1, 1): Step 4: Determine coefficients a and b using a Least Square Method (LSM): Step 5: Predict the next state as: whereê (1) (k + 1) is calculated using a whitening equation: Step 5: Predict the next state as: where (1) ( 1) e k + is calculated using a whitening equation:

Experiments
Four different control strategies were applied on the test-rig: MacPherson struts, FSASMC, FSB-ASMC+H ∞ and GP-FSB-ASMC+H ∞ . The ability of each strategy to reduce vibration was determined. When two Macpherson struts were used to suppress vibration, there was no additional force from the PM. FSB-ASMC means that the PM-ASS is operated on the FSB-ASMC; FSB-ASMC+H ∞ means that the PMASS is operated on the FSB-ASMC with a H ∞ performance guarantee and GP-FSB-ASMC+H ∞ means that the PMASS is operated on the FSB-ASMC with a H ∞ performance guarantee and GM(1, 1) error prediction. Two validations of the PM-ASS test-rig were performed: (1) It was shown that the techniques for the H ∞ performance guarantee and the state-predictor GM(1, 1) significantly improved the PM-ASS in terms of attenuating the displacement of the sprung mass and decreasing the acceleration of the sprung mass. (2) The performance of the Macpherson struts and the FSB-ASMC with/without the H ∞ performance guarantee and the state-predictor GM(1, 1) was compared. Figure 5 shows the test-rig, which has a PM and two MacPherson struts installed in parallel between the vehicle body and the support frame, a wheel firmly installed on one side of the support frame and a road profile generator that generates road profiles for various road conditions. The road profile generator has rollers, an induction motor and a pneumatic cylinder, and is in close contact with the wheel. Figure 6 illustrates the system blocks for a quarter-car with a PM-ASS. The components in Figure 6 are detailed as follows:

Test-rig for the PM-ASS
(1) The sprung mass is the portion of the vehicle's total mass that is supported by the suspension, i.e., a chassis.  Figure 5 shows the test-rig, which has a PM and two MacPherson struts installed in parallel between the vehicle body and the support frame, a wheel firmly installed on one side of the support frame and a road profile generator that generates road profiles for various road conditions. The road profile generator has rollers, an induction motor and a pneumatic cylinder, and is in close contact with the wheel. Figure 6 illustrates the system blocks for a quarter-car with a PM-ASS. The components in Figure 6 are detailed as follows:

Test-rig for the PM-ASS
(1) The sprung mass is the portion of the vehicle's total mass that is supported by the suspension, i.e., a chassis.      Figure 7 shows a control box, which has a NI cRIO-9074 embedded system, a NI-9263 D/A card, a NI-9215 A/D card and a NI-9411 encoder. The NI-9215 A/D card receives information (the green dashed line) about the displacement of the unsprung mass and the pneumatic cylinder and the acceleration of the sprung mass. The NI-9411 encoder receives a digital signal that shows the displacement of the sprung mass (the blue dashed line). The NI-9263AO sends the control signals (the red dashed line) to the PDVC, the invertor and the PPR to drive the PM, the pneumatic cylinder and the induction motor. All of the measured data was sent to a personal computer through an Ethernet interface and this is displayed on the screen. The system parameters can be modified in the PC and were sent to the NI cRIO-9074 embedded system. The Butter Worth filter was adopted to attenuate the sensor noise for the linear encoders, which is expressed as:    Figure 7 shows a control box, which has a NI cRIO-9074 embedded system, a NI-9263 D/A card, a NI-9215 A/D card and a NI-9411 encoder. The NI-9215 A/D card receives information (the green dashed line) about the displacement of the unsprung mass and the pneumatic cylinder and the acceleration of the sprung mass. The NI-9411 encoder receives a digital signal that shows the displacement of the sprung mass (the blue dashed line). The NI-9263AO sends the control signals (the red dashed line) to the PDVC, the invertor and the PPR to drive the PM, the pneumatic cylinder and the induction motor. All of the measured data was sent to a personal computer through an Ethernet interface and this is displayed on the screen. The system parameters can be modified in the PC and were sent to the NI cRIO-9074 embedded system. The Butter Worth filter was adopted to attenuate the sensor noise for the linear encoders, which is expressed as: where y out (k) is the output of the filter at the kth sampling and y in (k) is the displacement of the sprung mass measured by the linear encoder at the kth sampling.

Initial Setups for Experiments
It is critically important that the PM-ASS reacts quickly to changes in road conditions, thus, the PM was initially pre-set in the half-expanded state. The initial voltage of the PPR was 5V, in the range of [0V 10V] to maintain the PM in the half-expanded condition. Figure 8 shows three states for the PM: no-expansion (Figure 8a), the half-expanded state (Figure 8b) and the fully expanded state (Figure 8c). A rough concave-convex road profile and a two-bump excitation road profile were used to verify stability and robustness of the proposed PM-ASS. A fuzzy logic control (FLC) tracked the road profiles. Five membership functions were used for the error, the change in the error and the fuzzy output, as shown in Figure 9. The fuzzy rule refers to the rule template in Chapter of [Error! Reference source not found.], and the center of maximum defuzzification was used to calculate the fuzzy output. Table 1 shows the parameters for the FSB-ASMC, the FSB-ASMC+ H ∞ and the GP-FSB-

Initial Setups for Experiments
It is critically important that the PM-ASS reacts quickly to changes in road conditions, thus, the PM was initially pre-set in the half-expanded state. The initial voltage of the PPR was 5V, in the range of [0V 10V] to maintain the PM in the half-expanded condition. Figure 8 shows three states for the PM: no-expansion (Figure 8a), the half-expanded state (Figure 8b) and the fully expanded state (Figure 8c). A rough concave-convex road profile and a two-bump excitation road profile were used to verify stability and robustness of the proposed PM-ASS. A fuzzy logic control (FLC) tracked the road profiles. Five membership functions were used for the error, the change in the error and the fuzzy output, as shown in Figure 9. The fuzzy rule refers to the rule template in Chapter of [Error! Reference source not found.], and the center of maximum defuzzification was used to calculate the fuzzy output. Table 1 shows the parameters for the FSB-ASMC, the FSB-ASMC+ H ∞ and the GP-FSB-

Initial Setups for Experiments
It is critically important that the PM-ASS reacts quickly to changes in road conditions, thus, the PM was initially pre-set in the half-expanded state. The initial voltage of the PPR was 5V, in the range of [0V 10V] to maintain the PM in the half-expanded condition. Figure 8 shows three states for the PM: no-expansion (Figure 8a), the half-expanded state (Figure 8b) and the fully expanded state (Figure 8c). A rough concave-convex road profile and a two-bump excitation road profile were used to verify stability and robustness of the proposed PM-ASS. A fuzzy logic control (FLC) tracked the road profiles. Five membership functions were used for the error, the change in the error and the fuzzy output, as shown in Figure 9. The fuzzy rule refers to the rule template in Chapter of [25], and the center of maximum defuzzification was used to calculate the fuzzy output. Table 1 shows the parameters for the FSB-ASMC, the FSB-ASMC+H ∞ and the GP-FSB-ASMC+H ∞ . The rough concave-convex road profile is composed of a bump and a hollow with a sinusoidal disturbance, which is represented Equation 46.  1.554( 11.5) 3.5( 11.5) ( ) [10,11.5) ( )

Results and Disscussion of Experiments
where the sinusoidal disturbance is ( ) 0.175sin (2 ) d t πt = 0.07sin(7.5 ) πt + . Figure 10 shows the variation in the displacement of the sprung mass when the quarter car travels on a rough concaveconvex road with immediate 3 cm changes. In Figure 10, the blue line shows that the road profiles that are generated by the Road Generator have a 0.5 ± cm error in the rough concave-convex road profile. The pink dotted line shows the displacement of the sprung mass, which is controlled by the MacPherson struts. Because no additional force is supplied from PM in this case, the maximum displacement of the sprung mass is 2.4 cm, thus, only a 0.6 cm-high vibration is suppressed. Figure  10 shows that the FSB-ASMC, the FSB-ASMC+ H ∞ and the GP-FSB-ASMC+ H ∞ better suppress vibration than MacPherson struts under the same conditions. The black line in Figure 10 shows that FSB-ASMC suppresses a vibration in the sprung mass of less than 1.2 cm. The green dashed line shows that the FSB-ASMC+ H ∞ suppresses a maximum vibration of less than 0.65 cm because it  Table 1. Parameters for the proposed method for a rough concave-convex road profile.

Initial Weights
Control Parameters   The rough concave-convex road profile is composed of a bump and a hollow with a sinusoidal disturbance, which is represented Equation (46).

Results and Disscussion of Experiments
f or t ∈ [3.5, 5) 1.554(t − 6.5) 3 + 3.5(t − 6.5) 2 + d(t) f or t ∈ [5, 6.5) 1.554(t − 8.5) 3 − 3.5(t − 8.5) 2 + d(t) f or t ∈ [8.5, 10) −1.554(t − 11.5) 3 − 3.5(t − 11.5) 2 + d(t) f or t ∈ [10, 11.5) where the sinusoidal disturbance is d(t) = 0.175 sin(2πt) +0.07 sin(7.5πt). Figure 10 shows the variation in the displacement of the sprung mass when the quarter car travels on a rough concave-convex road with immediate 3 cm changes. In Figure 10, the blue line shows that the road profiles that are generated by the Road Generator have a ±0.5 cm error in the rough concave-convex road profile. The pink dotted line shows the displacement of the sprung mass, which is controlled by the MacPherson struts. Because no additional force is supplied from PM in this case, the maximum displacement of the sprung mass is 2.4 cm, thus, only a 0.6 cm-high vibration is suppressed. Figure 10 shows that the FSB-ASMC, the FSB-ASMC+H ∞ and the GP-FSB-ASMC+H ∞ better suppress vibration than MacPherson struts under the same conditions. The black line in Figure 10 shows that FSB-ASMC suppresses a vibration in the sprung mass of less than 1.2 cm. The green dashed line shows that the FSB-ASMC+H ∞ suppresses a maximum vibration of less than 0.65 cm because it attenuates lumped uncertainties. The GP-FSB-ASMC+H ∞ limits the maximum vibration to less than 0.4 cm, which is better than either the FSB-ASMC or the FSB-ASMC+H ∞ . Table 2 compares the maximum vibration and the root-mean-square error (RMSE) for the vibration, where the RMSE is defined as: where value vibration is the measured vibration, value zero_level is the zero displacement of the sprung mass and N is the total number of the measured vibration data points. Table 3 shows the weights after training for the GP-FSB-ASMC+H ∞ . In Figure 10, there is a slight oscillation on the sprung mass because of tire deformation, which is a small variation that depends on the road profile. Using the FLC, the variation in the output of the road generator is shown as a blue line in Figure 10, for which the tracking error is limited to ±0.5 cm. The acceleration and the control voltage and the sprung mass for all methods are respectively shown in Figures 11 and 12. It is noted that no additional force is applied to MacPherson struts, thus, Figure 12 shows no control voltage for them. is the zero displacement of the sprung mass and N is the total number of the measured vibration data points. Table 3 shows the weights after training for the GP-FSB-ASMC+ H ∞ . In Figure 10, there is a slight oscillation on the sprung mass because of tire deformation, which is a small variation that depends on the road profile. Using the FLC, the variation in the output of the road generator is shown as a blue line in Figure 10, for which the tracking error is limited to 0.5 ± cm. The acceleration and the control voltage and the sprung mass for all methods are respectively shown in Figures 11 and 12. It is noted that no additional force is applied to MacPherson struts, thus, Figure 12 shows no control voltage for them.   attenuates lumped uncertainties. The GP-FSB-ASMC+ H ∞ limits the maximum vibration to less than 0.4 cm, which is better than either the FSB-ASMC or the FSB-ASMC+ H ∞ . is the zero displacement of the sprung mass and N is the total number of the measured vibration data points. Table 3 shows the weights after training for the GP-FSB-ASMC+ H ∞ . In Figure 10, there is a slight oscillation on the sprung mass because of tire deformation, which is a small variation that depends on the road profile. Using the FLC, the variation in the output of the road generator is shown as a blue line in Figure 10, for which the tracking error is limited to 0.5 ± cm. The acceleration and the control voltage and the sprung mass for all methods are respectively shown in Figures 11 and 12. It is noted that no additional force is applied to MacPherson struts, thus, Figure 12 shows no control voltage for them. Figure 10. Displacement of the sprung mass for a rough concave-convex road profile. Figure 11. Acceleration of the sprung mass for the rough concave-convex road profile. Figure 11. Acceleration of the sprung mass for the rough concave-convex road profile.   (48) Figure 13 shows the variation in the displacement of the sprung mass when the quarter car moves along a rough sine wave excitation with a 1.5 cm amplitude and 2.5-sec wavelength. In Figure   13, the blue line shows the road profiles that are generated by the Road Generator, which gives a 0.5 ± cm error in the rough sine wave excitation road profile. The pink dotted line shows the displacement of the sprung mass, which is controlled by the MacPherson struts. There is almost no vibration-reduction. The black line in Figure 13 shows that FSB-ASMC suppresses the vibration of the sprung mass to less than 0.6cm. The green dashed line shows that the FSB-ASMC+ H ∞ suppresses the maximum vibration to less than 0.6 cm. The GP-FSB-ASMC+ H ∞ limits the maximum vibration to less than 0.25 cm, which is better than either the FSB-ASMC or the FSB-ASMC+ H ∞ . Table 4 compares the maximum vibration and the root-mean-square error (RMSE) of the vibration. Table 5 shows the weights after training for the GP-FSB-ASMC+ H ∞ . Using the FLC, the variation in the output of the  A sine wave excitation road profile is represented as: Figure 13 shows the variation in the displacement of the sprung mass when the quarter car moves along a rough sine wave excitation with a 1.5 cm amplitude and 2.5-sec wavelength. In Figure 13, the blue line shows the road profiles that are generated by the Road Generator, which gives a ±0.5 cm error in the rough sine wave excitation road profile. The pink dotted line shows the displacement of the sprung mass, which is controlled by the MacPherson struts. There is almost no vibration-reduction. The black line in Figure 13 shows that FSB-ASMC suppresses the vibration of the sprung mass to less than 0.6cm. The green dashed line shows that the FSB-ASMC+H ∞ suppresses the maximum vibration to less than 0.6 cm. The GP-FSB-ASMC+H ∞ limits the maximum vibration to less than 0.25 cm, which is better than either the FSB-ASMC or the FSB-ASMC+H ∞ . Table 4 compares the maximum vibration and the root-mean-square error (RMSE) of the vibration. Table 5 shows the weights after training for the GP-FSB-ASMC+H ∞ . Using the FLC, the variation in the output of the road generator is shown as the blue line in Figure 13 for which the tracking error is limited to ±0.1 cm. The acceleration and the control voltage and the sprung mass for all methods are, respectively, shown in Figures 14 and 15.
Appl. Sci. 2019, 9, x FOR PEER REVIEW 18 of 22 road generator is shown as the blue line in Figure 13 for which the tracking error is limited to 0.1 ± cm. The acceleration and the control voltage and the sprung mass for all methods are, respectively, shown in Figures 14 and 15.    road generator is shown as the blue line in Figure 13 for which the tracking error is limited to 0.1 ± cm. The acceleration and the control voltage and the sprung mass for all methods are, respectively, shown in Figures 14 and 15.    road generator is shown as the blue line in Figure 13 for which the tracking error is limited to 0.1 ± cm. The acceleration and the control voltage and the sprung mass for all methods are, respectively, shown in Figures 14 and 15.

Conclusions
This study developed an experimental ASS, which was termed the PM-ASS, on a quarter-car test rig. The design of the PM-ASS used a PM that is installed in parallel with two MacPherson struts and can provide a stable vertical force to reduce vibrations. In terms of the controller design, a GP-FSB-ASMC+H ∞ was used for the PM-ASS to address the nonlinearity and time-varying characteristics. Benefiting from the use of the orthogonal Fourier basis function, the FSB-ASMC+H ∞ has a clear physical meaning and readily determined structure for approximating the optimal control law to attenuate disturbances and uncertainties. In addition, the state-predictor GM(1, 1) predicted future errors to enhance the control performance for the FSB-ASMC+H ∞ . The experimental results demonstrated that the PM-ASS with the GP-FSB-ASMC+H ∞ significantly improved the Macpherson Strut suspension in terms of suppression of the displacement and acceleration in the vehicle body for improving the ride comfort.
To enhance the ability of vibration-reduction for the PM-ASS, we advise future researchers to use two air intakes on both side of the PM. A faster control response will then be produced for the PM-ASS. In terms of the GP-FSB-ASMC+H ∞ , a systemic method can be further developed for determining the initial weights and the controller parameters.

Patents
There are three Taiwan utility model patents resulting from the developed active suspension system, which are (1) a double-layer vibration isolation system using the pneumatic muscles with the patent number M437402; (2) a support system for the pneumatic muscles with the patent number M437403; (3) active vehicle suspension system using the pneumatic muscles with the patent number M438403.