Sliding Mode Switch Control of Adjustable Hydro-Pneumatic Suspension based on Parallel Adaptive Clonal Selection Algorithm

: The hydro-pneumatic suspension, as a widely used suspension for heavy vehicles, has been taken seriously by researchers for a long time because it is crucial in terms of handling stability, riding comfort, and driving safety of these vehicles. Most previous studies only discussed the control of ride comfort or vehicle handling stability of the suspension system separately. This article proposes a dynamic switch control strategy which can switch between ride comfort and handling stability controllers under di ﬀ erent road surfaces and driving conditions. The load transfer ratio (LTR) is selected as the switch performance index, and it is calculated through a six-degrees-of-freedom (6-DOF) model. The ride comfort and handling stability controller of the hydro-pneumatic suspension are designed based on the sliding mode control theory. The objective functions of parameters optimization of the sliding mode controller (SMC) are obtained by means of analytic hierarchy process (AHP), and then the controller’s parameters are optimized by the parallel adaptive clonal selection algorithm (PACSA). The simulation results based on MATLAB / Simulink show that: (1) the PACSA performs better than a genetic algorithm in terms of the parameters optimization of the SMC; (2) the proposed switch control strategy can simultaneously improve the ride comfort and handling stability under several typical steering maneuvers and various road proﬁles compared with the conventional SMC-controlled suspension. RCSMC and HDSMC are obtained. Then the control of hydro-pneumatic suspension is simulated based on the two optimized controllers respectively, which reﬂects the optimizing ability for SMC of the PACSA and GA. The control results are obtained as follows. This paper presents a sliding mode switch control strategy for both the ride comfort and handling stability of the adjustable hydro-pneumatic suspension system of a rescue vehicle. The switch index of the Switched SMC is deﬁned according to the value of LTR. Three di ﬀ erent SMC controllers are designed respectively and tested through the simulations of the half-car model based on MATLAB / Simulink. AHP is applied to calculate the weight coe ﬃ cients of the objective functions, and the parameters of the controller are tuned by PACSA, the optimization ability of which proved better than that of GA. The performance of the switch strategy is proved through simulations with di ﬀ erent road input (random road and bump road), di ﬀ erent vehicle speeds, and di ﬀ erent steering maneuvers (J-turn steering and Fishhook steering).


Introduction
Rescue vehicles are important in the transportation of rescue goods as well as rescue workers. The hydro-pneumatic suspension is applied to rescue vehicles due to the great load carrying capacity and the complexity of rescue terrain. The hydro-pneumatic suspension system plays a decisive role in the ride comfort and handling stability of the vehicle, but the passive suspension cannot address various road profiles because of its fixed parameters. As a result, it is meaningful and feasible to study a control strategy for the adjustable hydro-pneumatic suspension system which can perform better in ride comfort and handling stability than the traditional passive suspension [1].
In previous studies, many scholars have studied the control strategy of the adjustable hydro-pneumatic suspension. Andres Riofrio [2] improved the vehicle rollover stability of active suspension through an LQR (linear quadratic regulator)-based controller with estimation of the road. Shuai Wang [3] discussed the active control of hydro-pneumatic suspension parameters of wheel loaders based on a fuzzy controller, which reduces the vertical vibration of the wheel loader. Sarel F [4] presented the possibility of using slow active suspension control to reduce the body roll and The rest of this paper is organized as follows. Section 2 establishes the road model and the 6-DOF half-car model of the suspension system. Section 3 focuses on the switch control strategy of the hydro-pneumatic suspension. Section 4 describes the parameter tuning methods of the SMC. Section 5 carries out the simulation of the designed controllers, and the simulation results are analyzed. Section 6 is the conclusions.

System Modelling
For the controller design and the validation of the switch control strategy, the models of the hydro-pneumatic suspensions system, nonlinear tire input, and road input are established in this section.

Modelling of the Hydro-Pneumatic Suspension System
The vehicle studied in this paper is a rescue vehicle with an adjustable hydro-pneumatic suspension system. The corresponding model has six degrees of freedom which represent the vertical, roll, and lateral dynamics of the vehicle. The rest of this paper is organized as follows. Section 2 establishes the road model and the 6-DOF half-car model of the suspension system. Section 3 focuses on the switch control strategy of the hydro-pneumatic suspension. Section 4 describes the parameter tuning methods of the SMC. Section 5 carries out the simulation of the designed controllers, and the simulation results are analyzed. Section 6 is the conclusions.

System Modelling
For the controller design and the validation of the switch control strategy, the models of the hydro-pneumatic suspensions system, nonlinear tire input, and road input are established in this section.

Modelling of the Hydro-Pneumatic Suspension System
The vehicle studied in this paper is a rescue vehicle with an adjustable hydro-pneumatic suspension system. The corresponding model has six degrees of freedom which represent the vertical, roll, and lateral dynamics of the vehicle. The 6-DOF model consists of two parts, one is a 4-DOF roll motion model (Figure 1a  For convenience, we give the linear approximation for the damping and stiffness of the hydropneumatic suspension by fitting, and the fitted stiffness and damping of the suspension are shown in Table 1. In addition, two tires are modeled by means of two linear springs. The parameters of the rescue vehicle's suspension system are displayed in Table 1.  For convenience, we give the linear approximation for the damping and stiffness of the hydro-pneumatic suspension by fitting, and the fitted stiffness and damping of the suspension are shown in Table 1. In addition, two tires are modeled by means of two linear springs. The parameters of the rescue vehicle's suspension system are displayed in Table 1.
According to Figure 2 and Newton's law, the dynamic motions of the rescue vehicle are modelled as follows.   γ = F f y l f _F ry l r where ϕ, γ, β represent the roll angle, yaw rate, and side-slip angle, respectively; v x is the speed of the vehicle; Vertical motion equation of the spring mass: where z s denotes vertical displacement of sprung mass; f ul and f ur indicate the adjustable control forces of the suspension.
Vertical motion equation of the unspring mass: where z sl = z s + Tϕz sr = z s − Tϕ, z ql , z qr indicate the vertical displacements of the wheels excited by the road surface, z ul , z ur is the vertical displacements of unspring mass. The subscript l and r of the symbols represent the left and right suspension.

Tire Model
The tire mechanics property plays a vital role in the vehicle's dynamics performance, handling stability and comfort. The vertical dynamics and lateral dynamics of the vehicle are coupled to each other by the tire force; hence, the accuracy of the tire model is essential. Tire models are usually divided into theoretical models and empirical or semi-empirical models according to their structural points. The empirical tire is generated based on a large number of experimental data fittings [18], and it is generally accurate, as for example, Pacejka's tire model. Fiala and UA tire belong to the theoretical tire model, which is inaccurate but physically meaningful [19]. WEI Yintao [20] proposed that the theoretical physical model is impossible to practically apply to automobile dynamics simulations because it lacks the accuracy of quantitative description. Yasheen Babulal [21] parameterized three existing tire models namely, Fiala, UA, and Pacejka tire models, and he drew a conclusion that the Pacejka tire model tire provided the closest fit for the side-force versus slip angle tests of all three tire models. The Pacejka tire model shows the relationship among the longitudinal force, the lateral force, the aligning torque and the side angle. It has a high fitting accuracy and can be used within a certain range except for extreme conditions. Therefore, the Pacejka tire model is selected in this paper to calculate the tire lateral force coupled with the vertical load. The relevant expression of the tire model can be found in Ref [19]. The Pacejka tire model is established in MATLAB/Simulink, and the correlative curves under the steering condition of Figure 2 are shown in Figure 3. It can be seen from the tire simulation results in Figure 3, the values and the trend of curves are consistent with the curves in Ref [22], which illustrates the correctness of the tire model established in this article. It can be seen from the tire simulation results in Figure 3, the values and the trend of curves are consistent with the curves in Ref [22], which illustrates the correctness of the tire model established in this article.

Road Input Modelling
All analyses and studies of the suspension system are implemented on the basis of the input of the stochastic road surface, because the excitation of road wave surface is the main vibration source of the vehicle vibration system [23]. Accurate road models are important to the research and optimization of the suspension system. Random road surface is generally selected as the most commonly used in simulation, and more results and significance can be obtained in Ref [24]. The expressions of z ql , z qr are shown in Equation (6).
where w(t) is white noise signal of the road, and f 0 = v x × n 00 .X 1 and X 2 are the state variables in the state equation. The random road model is established based on MATLAB/Simulink, and the simulated displacement curves of the left and right tires are shown in Figure 4,

Switch Control Strategy
Based on the conflicting relationship between them, it is difficult to optimize both aspects of ride comfort and handling stability of the rescue vehicle based on the control of the hydro-pneumatic suspension. However, when the vehicle is in different driving conditions or driving on different road surfaces, level of magnitude of the two aspects can be sorted. For instance, when the road surface is relatively flat and the driving state of the vehicle is unlikely to cause the vehicle to be dangerous, the ride comfort is the main goal pursued by the suspension system in this case. On the contrary, when the vehicle is driving on an uneven road or steering at a high speed, the possibility of rolling of the vehicle is relatively large at this time, and the handling stability of the vehicle is the major demand of the vehicle. As a result, it should be possible to analyze the weight of the vehicle's ride comfort and the handling stability under different conditions, so as to achieve the corresponding switch according to different requirements of working conditions. For the above-mentioned idea, the switch control strategy is presented in this paper, and the schematic diagram of the switch control strategy

Switch Control Strategy
Based on the conflicting relationship between them, it is difficult to optimize both aspects of ride comfort and handling stability of the rescue vehicle based on the control of the hydro-pneumatic suspension. However, when the vehicle is in different driving conditions or driving on different road surfaces, level of magnitude of the two aspects can be sorted. For instance, when the road surface is relatively flat and the driving state of the vehicle is unlikely to cause the vehicle to be dangerous, the ride comfort is the main goal pursued by the suspension system in this case. On the contrary, when the vehicle is driving on an uneven road or steering at a high speed, the possibility of rolling of the vehicle is relatively large at this time, and the handling stability of the vehicle is the major demand of the vehicle. As a result, it should be possible to analyze the weight of the vehicle's ride comfort and the handling stability under different conditions, so as to achieve the corresponding switch according to different requirements of working conditions. For the above-mentioned idea, the switch control strategy is presented in this paper, and the schematic diagram of the switch control strategy of the hydro-pneumatic suspension is shown in Figure 5.  Figure 5. Sliding mode switch control strategy of adjustable hydro-pneumatic suspension.
where f δ is the steering angle of the wheel, and y a is the lateral acceleration of the spring mass, For the switch control strategy, an indicator is needed to determine whether and when the switch occurs. In order to evaluate the possibility of the vehicle's danger, LTR is selected as the performance index to judge whether the switch action between the ride comfort and handling stability controllers is triggered, and LTR is defined as follows.
When the LTR is equal to 0, the vehicle has a stable roll dynamic, and the risk becomes higher with the increase of LTR [25]. The threshold of LTR is defined as Q, then the switch control strategy can be illustrated in Equation (8). In order to ensure sufficient time for the controller to react before the rollover danger occurs, Q is set to 0.8 in this paper. On this occasion, the controller has enough time to lower the LTR to avoid the rollover before obtaining high values of LTR. where δ f is the steering angle of the wheel, and a y is the lateral acceleration of the spring mass, For the switch control strategy, an indicator is needed to determine whether and when the switch occurs. In order to evaluate the possibility of the vehicle's danger, LTR is selected as the performance index to judge whether the switch action between the ride comfort and handling stability controllers is triggered, and LTR is defined as follows.
Appl. Sci. 2020, 10, 1852 When the LTR is equal to 0, the vehicle has a stable roll dynamic, and the risk becomes higher with the increase of LTR [25]. The threshold of LTR is defined as Q, then the switch control strategy can be illustrated in Equation (8). In order to ensure sufficient time for the controller to react before the rollover danger occurs, Q is set to 0.8 in this paper. On this occasion, the controller has enough time to lower the LTR to avoid the rollover before obtaining high values of LTR.
controller = ride com f ort controller, LTR ≤ Q handing stability controller, LTR > Q In order to easily understand the proposed control strategy, the control block diagram of the switch control strategy of the hydro-pneumatic suspension is shown in Figure 6. As shown in Figure 6, when the value of LTR is lower than Q, it indicates that the vehicle is not at risk of rollover at this time. As a result, the major task of the controller is to optimize the ride comfort of the suspension system, then the ride comfort SMC (RCSMC) is active. The handling stability SMC (HDSMC) is not activated until the LTR reaches the threshold Q.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 9 of 31 In order to easily understand the proposed control strategy, the control block diagram of the switch control strategy of the hydro-pneumatic suspension is shown in Figure 6. As shown in Figure 6, when the value of LTR is lower than Q, it indicates that the vehicle is not at risk of rollover at this time. As a result, the major task of the controller is to optimize the ride comfort of the suspension system, then the ride comfort SMC (RCSMC) is active. The handling stability SMC (HDSMC) is not activated until the LTR reaches the threshold Q.

Handling Stability Control
Once the gravity center of the heavy vehicle is laterally offset owing to severe maneuvers, a centrifugal moment ( c M ) is generated, and the vehicle is at risk of rolling over. The actuator of the adjustable hydro-pneumatic suspension can generate a torque ( t M ) against the centrifugal moment under the appropriate control of the controller, thereby preventing the vehicle from rolling over. Consequently, the following vehicle roll dynamic model is obtained.

( )
sin , where κ represents the allowable uncertainties, and For the purpose of effectively improving the asymptotic stability and decreasing chattering aroused from SMC's switch characteristic, the exponential reaching law is selected to converge towards 0 exponentially in this paper, which is rewritten in Equation (11).

Handling Stability Control
Once the gravity center of the heavy vehicle is laterally offset owing to severe maneuvers, a centrifugal moment (M c ) is generated, and the vehicle is at risk of rolling over. The actuator of the adjustable hydro-pneumatic suspension can generate a torque (M t ) against the centrifugal moment under the appropriate control of the controller, thereby preventing the vehicle from rolling over. Consequently, the following vehicle roll dynamic model is obtained.
ϕ . Both κ 1 κ 2 are positive constants. As we can see from Equation (9), the system input is M t , and the output is ϕ. From SMC theory, the sliding surface is defined as below.
where ρ(ρ > 0) indicates the performance of system converging to the equilibrium point. e is the difference between the ideal output and actual output, and the ideal output of the roll angle is zero, in other words, e = ϕ.
For the purpose of effectively improving the asymptotic stability and decreasing chattering aroused from SMC's switch characteristic, the exponential reaching law is selected to converge towards 0 exponentially in this paper, which is rewritten in Equation (11).
In order to satisfy the Hurwitz criteria [26], Equation (12) is selected as the Lyapunov candidate function, which is strictly negative. Equation (13) should be satisfied to ensure the robustness stability of the controller.
where µ and η are the coefficients that represent the speed reaching the sliding surface, µ > 0 η > 0. According to Equations (10)~(13) mentioned above, the control input can be rewritten as Equation (14).

Ride Comfort Control
Through many years of investigation and practice, the quarter-car suspension nonlinear model is suitable for the research on the ride comfort control [27]. In this paper, the ride comfort controller of the hydro-pneumatic suspension is built based on the quarter-car model, which makes it easy to compute. Through this quarter-car model, the controller output ( f ul ) of one side of the suspension can be calculated, and the other side controller output is expressed as f ur = ζ f ul , where ζ is the correlation coefficient that connects the control quantity of left and right sides' suspensions. The value of ζ can be obtained through the optimization algorithm which will be introduced in the Section 4. According to Newton's law, the state space representation can be written as follows. where Define the sliding surface as where X is the state vector, and the polynomial p 3 + n 3 p 2 + n 2 p + n 1 should satisfy the Hurwitz law based on the design criteria of the SMC. [28] On the basis of modern control theory, the variable structure systems' controllability should be constant for any nonsingular transformation. Because of the controllability of the matrices A and B, the nonsingular transformation is carried on X.
The values of the matrices in Equations (17) and (18) are shown as follows.
x 2 The sliding surface can be rewritten according to Equations (16) to (18). where , m u − m s n 3 0 For the purpose of satisfying the Hurwitz criteria to ensure the stability according to the SMC theory, Equation (21), the characteristic polynomial of Equation (20), should ensure the roots are negative. In other words, the value of p must be greater than zero. As a consequence, the bounds of n 1 , n 2 , and n 3 are confirmed, as shown in Equation (22).
Similarly, the exponential reaching law is also applied to the design of the ride comfort controller, and the controller output ( f ul ) can be obtained.

SMC Parameters Tuning Based on PACSA
During the process of designing the SMC, parameters of SMC should be chosen suitably through some optimization algorithms to obtain optimal performance. For the RCSMC, there are six parameters that can be tuned, i.e., n 1 , n 2 , n 3 , η, µ, and ζ. Likewise, five related parameters of the HDSMC need to be optimized, which are ρ, µ, η, κ 1 , and κ 2 . In this section, two objective functions of parameters optimization are designed through AHP, and then the PACSA is applied to optimize the above parameters of the two controllers based on the objective functions.

Objective Function
Appropriate performance evaluation indicators are the basis for judging the quality of control strategies. Ride comfort and handling stability are the main aspects of vehicle performance during the running process of vehicle. In order to achieve optimal SMC of vehicle ride comfort and handling stability, the objective function of the SMC' parameters optimization consists of vehicle vertical acceleration .. z s , dynamic deflection of left and right tire z ql − z ul and z qr − z ur , left and right side suspension dynamic deflection z ul − z sl and z ur − z sr , and roll angle acceleration of vehicle ..
ϕ. The first three items are the objective function factors often used in suspension ride control [29,30]. The roll angle acceleration is the main control target in the handling stability control of the suspension system [31,32]. In order to trade off among the various performance requirements for better ride comfort and handling stability, a quadratic objective function expression is taken as the performance index based on the quadratic optimal control theory, which is shown in Equation (24).
where τ 1 ∼ τ 6 are the weighting coefficients of the corresponding items. First of all, the six weighting indicators mentioned above should be quantified because of the big differences in magnitude between them. The root mean square (RMS) value of the above six indicators in the passive suspension model is calculated, and the quantized scale factor of the same scale can be obtained from Equation (25).
where γ 1 ∼ γ 6 are the quantized scale factors. The vertical displacement acceleration .. z s is used as a reference, and the value of γ 1 is set to 1. The calculation results are shown in Table 2.

Weighting Coefficient Optimization of the Objective Function based on AHP
The weighting coefficients have a great influence on the value of the objective function; therefore, it is essential to reasonably select the weighting coefficients of the objective function. In this study, the weighting coefficients are determined by AHP to avoid multifarious trial and error.
AHP is used for multi-criteria decision making in which the factors are arranged in a hierarchical structure [33]. The weight coefficients evaluation process of the objective function based on AHP is shown in Figure 7.
The weighting coefficients have a great influence on the value of the objective function; therefore, it is essential to reasonably select the weighting coefficients of the objective function. In this study, the weighting coefficients are determined by AHP to avoid multifarious trial and error.
AHP is used for multi-criteria decision making in which the factors are arranged in a hierarchical structure [33]. The weight coefficients evaluation process of the objective function based on AHP is shown in Figure 7.

Set up hierarchical analysis model
The hierarchical analysis model is displayed in Figure 8, where the first level evaluation feature is composed of ride comfort and handling stability during the vehicle driving. The second level evaluation feature includes vehicle vertical acceleration, dynamic deflection of left and right tire, leftand right-side suspension dynamic deflection, and roll angle acceleration of the vehicle. Therefore, the structure of the hierarchical analysis model is derived from the objective function to acquire the priority of the evaluation feature on each level.

Set up hierarchical analysis model
The hierarchical analysis model is displayed in Figure 8, where the first level evaluation feature is composed of ride comfort and handling stability during the vehicle driving. The second level evaluation feature includes vehicle vertical acceleration, dynamic deflection of left and right tire, leftand right-side suspension dynamic deflection, and roll angle acceleration of the vehicle. Therefore, the structure of the hierarchical analysis model is derived from the objective function to acquire the priority of the evaluation feature on each level.

Construct judgment matrix T
The judgment matrix T is applied to distinguish the relative emphasis of evaluation features per hierarchical layer. T is constructed by using a nominal scale from 0 to 9 with the value ij a , which is assigned to express the relative importance of evaluation element i a over j a . The matrix is established by the importance comparison between evaluation features.

Calculate weight vector corresponding to T
The weight vector can be written as

Construct judgment matrix T
The judgment matrix T is applied to distinguish the relative emphasis of evaluation features per hierarchical layer. T is constructed by using a nominal scale from 0 to 9 with the value a ij , which is assigned to express the relative importance of evaluation element a i over a j . The matrix T = a ij is established by the importance comparison between evaluation features.

Calculate weight vector corresponding to T
The weight vector can be written as Λ = (τ 1 , τ 2 , . . . , τ r ) T , where r denotes numbers of evaluation features per level. Λ is gained through the normalization process for judgment matrix T.

Consistency check of judgment matrix
Generally, CI (consistency index) is used to check the consistency deviation of the judgment matrix, whose expression is shown as follows.
where m represents the order of judgment matrix T. η max indicates the root of maximum eigenvector value. At the same time, the formula of η max is displayed as follows.
In addition, in order to check the consistency more effectively, CR (ratio of consistency check) is introduced as shown in the following equation.
where RI is the reference criterion, and Table 3 shows the value of RI. If the value of CR is less than 0.10, it means that the consistency check is qualified.

Calculate the weight of evaluation features per level
If the judgment matrix satisfies the consistency, then the weight of the evaluated feature on each level is the eigenvector obtained; otherwise, the eigenvector needs to be calculated. The weights of evaluation features on two levels can be calculated by means of the four steps mentioned above. If τ f irst = ρ i , τ sec ond = q j , then the weight of the lowest level feature is described as follows.
where τ ij represents the weight of jth feature at the second-level relative to the weight of ith feature at the first level. γ i is the quantization scale factor. Above all, the weight coefficients of the objective function of the SMC' parameters optimization in Equation (24) are shown in Table 4. It can also be seen from Table 4 that the vehicle vertical acceleration and the roll angle acceleration of vehicle are the two main factors which can reflect the ride comfort and handling stability of the vehicle, respectively. As a result, the two objective functions for the optimization of the RCSMC and HDSMC are shown in Equation (29).

Brief Introduction of PACSA
The natural biological immune system has greatly influenced the development of computational intelligence in recent years, and many scholars have proposed many intelligent computing technologies based on immune systems. The artificial immune algorithm is abstracted from the mechanism of the biological immune system. The biological principles like clone generation, proliferation, and maturation are mimicked and incorporated into an artificial immune-based algorithm termed as the clonal selection algorithm (CSA) [34]. Many concepts and operators in the algorithm have a corresponding relationship with the concepts and mechanisms in the immune system, which should be introduced. Antigen is a substance that can stimulate the body to produce specific antibodies, such as to viruses or bacteria. Antibodies can bind to an antigen that invades from the outside and then help destroy these antigens to keep the body healthy. However, the extent of binding of different antibodies to antigens is different, and the higher the matching extend, the stronger the ability of antibodies to destroy such antigens. Antibody-antigen affinity is used to describe the matching extent between the antibody and the antigen. In the immune algorithm, the corresponding meanings of the antigen, antibody, and antibody-antigen affinity correspond respectively to: (1) the objective function and constraints of the optimization problem; (2) the optimal solution of the optimization problem; (3) the matching between the solution of the optimization problem and the objective function.
PACSA is proposed on the basis of ACSA. 'P' stands for the parallel computing, and the first letter 'A' represents the tune of the parameters is adaptive, which improves the running speed and convergence of the CSA, respectively. Parallel computing can greatly reduce the running time of optimization operations and improve the performance of CSA through sacrificing hardware resources [16,17]. ACSA enables the cloning operator to adjust the population evolution strategy dynamically based on the evolution of the population.
PACSA refers to random mapping of antibodies induced by antibody-antigen affinity, where the meaning of cloning corresponds to the generation of new antibodies based on genetic operators. The main idea of the PACSA could be described as follows: only the antibodies recognized by the antigen are cloned during the process of optimization; in other words, the antibody is selected according to the antibody-antigen affinity in each generation for the cloning operation in each generation to generate a new population. The core operation of PACSA is the cloning operator, which ensures the algorithm converges quickly in the direction of the global optimal solution based on the cloning operation and selection mechanism. The clonal selection operation not only expands the spatial range for the algorithm to search for the optimal solution, but also increases the diversity of antibodies through the mutation operation, helps prevent premature convergence, avoids plunging into a local optimum, and accelerates the convergence of the group through selection operations. According to the principle of self-adjusting, the mutation rate of PACSA can be obtained from Equation (30) based on antibody-antigen affinity. It can be seen that the closer the number of iterations to the maximum set algebra, the smaller the mutation rate. As a consequence, the mutation rate is higher at the beginning of optimization, so that the antibody can jump out of the local extreme point and move to the global maximum point quickly. The mutation rate gradually decreases with the progress of iteration, making the antibody infinitely close to the global optimal solution, thereby enhancing the local search ability.
where A expresses the antibody-antigen affinity; A avg is average value of A; p r indicates the high frequency mutation rate of the antibody to be mutated; p r_max and p r_min are the maximum value and minimum value of the mutation rate, respectively; x n represents the current evolutional algebra; x max is the maximum evolutional algebra.
Based on antibody-antigen affinity, the clone scale of antibody is calculated as follows.
where A min is the minimum value of antibody-antigen affinity; n m represents the base set value of the clone scale; n indicates the number of antibodies that need to be cloned. The symbol 'int[]' is the rounded-down function.
As can be seen from Formula (31), the number of clones of an antibody is proportional to its own affinity value. In other words, if the antibody-antigen affinity is less than the average value (the antibody is poor), the antibody will be inhibited. On the contrary, if the antibody is superior, then the antibody will be propagated, thus realizing the survival of the fittest of antibody evolution.
According to the antibody-antigen affinity and evolutionary algebra to adjust the number of updates per generation, the PACSA's adaptive update number is calculated as follows.
where λ is number of updates per generation; λ max and λ min are the set maximum value and minimum value of the updates respectively; A max and A min are the optimal and worst values of antibody-antigen affinity in the antibody population. It can be seen from Equation (30) that the antibody-antigen affinity values have a great discrepancy in the initial stage of optimization, and the mean value is small that λ is larger and the probability of the antibody being updated is higher, which ensures the diversity of antibodies. With the evolution of antibodies, the optimal and worst antibody-antigen affinity are gradually approached, and the average value is gradually increased; hence, λ is smaller. As a result, the probability of the antibody being updated is lower in the last part of the iteration.

Optimization Process of PACSA
The procedure of PACSA for optimizing the SMC's parameters is shown in Figure 9, and the steps are described in detail as follows.
Step 1. Parameter setting. Set the parameters of the PACSA, including the clone scale n m , the population number of per generation initial antibody population n p , the number of updates per generation λ, and the selection number of original antibody per generation n s .
Step 2. Initialization. The initial value of the objective function (which is the initial antibody population in this article) Y consists of Z, a binary-encoded antibody generated at random.
Step 3. Calculation of antibody-antigen affinity. In this paper, the antibody-antigen affinity is calculated according to the controller's error integration criterion. The smaller the error integral value, the higher the antibody-antigen affinity. The error integral is calculated with the running of the simulation model. This process is the most computationally intensive part in the whole algorithm, and the error integral calculation for each simulation is independent from others. Therefore, the computing task is assigned to multiple threads of the computer for parallel computing.
Step 4. Clone. Cloning and proliferating the antibody in the antibody population Y to obtain an amplified antibody population R. Whether the antibody is cloned, as well as the number of clones, is determined by Equation (26).
Step 5. Somatic hyper-mutation. Somatic hyper-mutation of antibodies in the antibody population R to obtain R ⊗ . The rate of somatic hyper-mutation is defined by Equation (25).
Step 6. Selection. Select λ high-affinity antibodies from the antibody population R ⊗ to replace λ low-affinity antibodies in Y. The value of λ can be obtained from Equation (27).
Step 7. Judgment of the iteration termination condition. Whether the number of iterations meets the requirements, or whether the antibody-antigen affinity of the optimal antibody in the current population meets the requirements, and if it is satisfied, the procedure is terminated, and the result is the output. Otherwise, return to step (4).
an amplified antibody population R . Whether the antibody is cloned, as well as the number of clones, is determined by Equation (26).
Step 5. Somatic hyper-mutation. Somatic hyper-mutation of antibodies in the antibody population R to obtain R ⊗ . The rate of somatic hyper-mutation is defined by Equation (25).
Step 6. Selection. Select λ high-affinity antibodies from the antibody population R ⊗ to replace λ low-affinity antibodies in Y . The value of λ can be obtained from Equation (27).
Step 7. Judgment of the iteration termination condition. Whether the number of iterations meets the requirements, or whether the antibody-antigen affinity of the optimal antibody in the current population meets the requirements, and if it is satisfied, the procedure is terminated, and the result is the output. Otherwise, return to step (4).

Simulations Results and Analysis
The simulation model including the 6-DOF model, the ride comfort controller, the handling stability controller, the switch strategy, and the road input models are built in MATLAB/Simulink, as shown in Figure 10. In order to verify the quality of the proposed switch control strategy, the

Simulations Results and Analysis
The simulation model including the 6-DOF model, the ride comfort controller, the handling stability controller, the switch strategy, and the road input models are built in MATLAB/Simulink, as shown in Figure 10. In order to verify the quality of the proposed switch control strategy, the simulations are performed under the J-turn maneuver and fishhook maneuver to evaluate the ride comfort and handling stability of the vehicle.

SMC Controller Tuning Results of PACSA
The GA-based SMC controller is designed to contrast the performance of the proposed SMC controller optimized by PACSA. The relationships between the fitness functions value and the genetic generation of GA and PACSA are shown in Figure 11. Figure 10 displays the changes in the objective function of RCSMC and HDSMC during the PACSA tuning optimization. As can be seen from Figure 10, the changes in the objective function of RCSMC tuned by the PACSA are obvious in the first 40 generations and the corresponding convergence speed is fast. However, even until 100 generations, no convergence occurred during the GA tuning process. In addition, the whole tuning process of the GA lasted for 33 h, and the PACSA optimization procedure took only 6 h.
The optimized parameter results based on the GA and the PACSA are listed in the Table 5.
According to the optimized parameters in Table 5, the optimized RCSMC and HDSMC are obtained. Then the control of hydro-pneumatic suspension is simulated based on the two optimized controllers respectively, which reflects the optimizing ability for SMC of the PACSA and GA. The control results are obtained as follows.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 18 of 31 simulations are performed under the J-turn maneuver and fishhook maneuver to evaluate the ride comfort and handling stability of the vehicle.

SMC Controller Tuning Results of PACSA
The GA-based SMC controller is designed to contrast the performance of the proposed SMC controller optimized by PACSA. The relationships between the fitness functions value and the genetic generation of GA and PACSA are shown in Figure 11.

SMC Controller Tuning Results of PACSA
The GA-based SMC controller is designed to contrast the performance of the proposed SMC controller optimized by PACSA. The relationships between the fitness functions value and the genetic generation of GA and PACSA are shown in Figure 11.    Above all, it can be concluded that the convergence rate, running time, and optimized results of the PACSA approach are better than the GA approach, which indicates that PACSA is more suitable for tuning the parameters of the SMC controller.

Results of Switch Control Strategy
According to the optimization results mentioned in Table 5, the RCSMC and the HDSMC are obtained. A comparison between the proposed switch control SMC (Switched SMC) and two conventional SMCs (HDSMC and RCSMC) is discussed, which can better assert that the proposed control is superior. In order to demonstrate the effectiveness of the proposed switch control strategy, different steering maneuvers including J-turn and Fishhook are tested. Furthermore, to better illustrate the adaptability of control strategy, two different typical road conditions including random road and bump road are built. The related parameters of the vehicle driving state and the corresponding road input are listed in Table 6.  Above all, it can be concluded that the convergence rate, running time, and optimized results of the PACSA approach are better than the GA approach, which indicates that PACSA is more suitable for tuning the parameters of the SMC controller.

Results of Switch Control Strategy
According to the optimization results mentioned in Table 5, the RCSMC and the HDSMC are obtained. A comparison between the proposed switch control SMC (Switched SMC) and two conventional SMCs (HDSMC and RCSMC) is discussed, which can better assert that the proposed control is superior. In order to demonstrate the effectiveness of the proposed switch control strategy, different steering maneuvers including J-turn and Fishhook are tested. Furthermore, to better illustrate the adaptability of control strategy, two different typical road conditions including random road and bump road are built. The related parameters of the vehicle driving state and the corresponding road input are listed in Table 6.

J-Turn Maneuver
The random road input is shown in Figure 13a, and the vehicle steering angle of J-turn maneuver is displayed in Figure 13b. The dynamic responses of the HDSMC-controlled vehicle (blue lines), RCSMC-controlled vehicle (black lines), and the vehicle with the Switched SMC (red lines) are shown in Figure 14.

J-Turn Maneuver
The random road input is shown in Figure 13a, and the vehicle steering angle of J-turn maneuver is displayed in Figure 13b. The dynamic responses of the HDSMC-controlled vehicle (blue lines), RCSMC-controlled vehicle (black lines), and the vehicle with the Switched SMC (red lines) are shown in Figure 14. It can be seen from Figure 13b that the vehicle is going straight in the time range of [0, 5] s, and the LTR value is much smaller than the threshold. During this period, the main task for the Switched SMC is to pursue the ride comfort. Therefore, the vehicle's vertical acceleration curves in Figure 14a during this period by the RCSMC and Switched SMC are basically identical, and the ride comfort of HDSMC is noticeably worse than that of the other two controllers. Moreover, the roll angle of the vehicle is small with different controllers during the time of [0, 5] s. The J-turn begins at the 5th second; then the vehicle roll angle and the LTR increases rapidly. The sudden steering has an impact on the ride comfort of the vehicle (point I), which is can be seen in Figure 14a. As a result, the switch of Switched SMC is not active until the LTR reaches the threshold at the 6.6th s; then the main purpose of the Switched SMC changes from the ride comfort of the vehicle to the handling stability of the vehicle. As shown in Figure 14b-f, during [6.6, 10], the values of the roll angle are reduced significantly and the LTR is dragged towards to the set value with the control of the Switched SMC. However, the LTR of the vehicle with the RCSMC is increasing with the steering of the vehicle, which is likely to cause unsafe driving conditions that lead to accidents and delays. In addition, the switch action also can be confirmed by the switch value in Figure 14e. When the switch value changes from 0 to positive, then the switch occurs. In addition, the switch action of the Switched SMC causes to a certain extent an impact on the ride comfort of the vehicle, which is point II shown in Figure 14a. The frequency domain response of the body acceleration and roll angle are shown in Figure 14c,d, which provide further obvious vindication of the Switched SMC. On the one hand, the curve in Figure 14c indicates that the Switched SMC performs better than the HDSMC with respect to ride comfort. On the other hand, curves in Figure 14d validate the superior stability of the vehicle with Switched SMC than the RCSMC. In other words, Switched SMC can take care of both sides of ride comfort and handling stability rather than one aspect. It can be seen from Figure 13b that the vehicle is going straight in the time range of [0, 5] s, and the LTR value is much smaller than the threshold. During this period, the main task for the Switched SMC is to pursue the ride comfort. Therefore, the vehicle's vertical acceleration curves in Figure 14a during this period by the RCSMC and Switched SMC are basically identical, and the ride comfort of HDSMC is noticeably worse than that of the other two controllers. Moreover, the roll angle of the In order to validate the robustness of the suspension system, simulations considering the uncertainty of system parameters and external disturbances were performed. In this paper, road disturbance was used as the external uncertain disturbance factor. The spring mass, vehicle speed, and suspension damping were used as the internal uncertain disturbance of the system.
Taking the working conditions in Figure 13 as the initial conditions without disturbance, and the road disturbance signal is shown in Figure 15, which is imposed on the left tire. It can be seen from Figure 15, two road disturbances are imposed on the system respectively, and disturbance-2 is more serve than disturbance-1. On the one hand, the acceleration curves of the three conditions are roughly similar in the aspects of amplitude and trend in Figure 16a. On the other hand, the roll angle curves in Figure 16b show almost no difference compared with the initial system without disturbance. In other words, the ride comfort and the stability are stable, and the suspension system has excellent immunity to external disturbances. Taking into account the typical non-linear characteristics of the suspension system, the following three parameters of the suspension system, unspring mass, damping, and vehicle speed, are used as internal uncertainties to validate the system robustness. The corresponding simulation results are shown from Figures 17-19. It can be seen from Figure 15, two road disturbances are imposed on the system respectively, and disturbance-2 is more serve than disturbance-1. On the one hand, the acceleration curves of the three conditions are roughly similar in the aspects of amplitude and trend in Figure 16a. On the other hand, the roll angle curves in Figure 16b show almost no difference compared with the initial system without disturbance. In other words, the ride comfort and the stability are stable, and the suspension system has excellent immunity to external disturbances. It can be seen from Figure 15, two road disturbances are imposed on the system respectively, and disturbance-2 is more serve than disturbance-1. On the one hand, the acceleration curves of the three conditions are roughly similar in the aspects of amplitude and trend in Figure 16a. On the other hand, the roll angle curves in Figure 16b show almost no difference compared with the initial system without disturbance. In other words, the ride comfort and the stability are stable, and the suspension system has excellent immunity to external disturbances. Taking into account the typical non-linear characteristics of the suspension system, the following three parameters of the suspension system, unspring mass, damping, and vehicle speed, are used as internal uncertainties to validate the system robustness. The corresponding simulation results are shown from Figures 17-19. Taking into account the typical non-linear characteristics of the suspension system, the following three parameters of the suspension system, unspring mass, damping, and vehicle speed, are used as internal uncertainties to validate the system robustness. The corresponding simulation results are shown from Figures 17-19. Appl. Sci. 2020, 10, x FOR PEER REVIEW 24 of 31 The systems with three different unspring masses are compared in Figure 17. The black curve represents the initial condition with a sprung mass of 1525 kg, and the red and blue curves represent an increase and decrease of 20% of the sprung mass on the basis of the initial condition, respectively. As shown in Figure 16a, the differences between the three curves are small both in terms of trend and amplitude. The roll angle curves in Figure 17b show that the trends of the three conditions are basically the same, and the amplitude is slightly different. The systems with three different unspring masses are compared in Figure 17. The black curve represents the initial condition with a sprung mass of 1525 kg, and the red and blue curves represent an increase and decrease of 20% of the sprung mass on the basis of the initial condition, respectively. As shown in Figure 16a, the differences between the three curves are small both in terms of trend and amplitude. The roll angle curves in Figure 17b show that the trends of the three conditions are basically the same, and the amplitude is slightly different. The systems with three different unspring masses are compared in Figure 17. The black curve represents the initial condition with a sprung mass of 1525 kg, and the red and blue curves represent an increase and decrease of 20% of the sprung mass on the basis of the initial condition, respectively. As shown in Figure 16a, the differences between the three curves are small both in terms of trend and amplitude. The roll angle curves in Figure 17b show that the trends of the three conditions are basically the same, and the amplitude is slightly different. The systems with three different unspring masses are compared in Figure 17. The black curve represents the initial condition with a sprung mass of 1525 kg, and the red and blue curves represent an increase and decrease of 20% of the sprung mass on the basis of the initial condition, respectively. As shown in Figure 16a, the differences between the three curves are small both in terms of trend and Appl. Sci. 2020, 10, 1852 23 of 30 amplitude. The roll angle curves in Figure 17b show that the trends of the three conditions are basically the same, and the amplitude is slightly different.
As shown in Figure 18, the red and blue curves represent an increase and decrease of 20% of the damping comparing to the black curve of the initial condition, respectively. Obviously, changes in suspension damping in this range have almost no effect on the acceleration and roll angle of the vehicle system. As a result, the suspension system in this paper has a good anti-interference ability to the parameter uncertainty of suspension damping. Figure 19 shows the impact of vehicle speed changes on the suspension system. The vehicle speed under the initial condition is 70 km/h, which is the maximum speed of the rescue vehicle in this paper. Therefore, reducing vehicle speed by 20% and 30%, respectively, was used to verify the interference caused by the uncertainties of vehicle speed on the system. As can be seen from the two diagrams in Figure 19, the system has a strong anti-interference ability to vehicle speed uncertainty.
As a result, it can be seen from the above simulation results and analysis that the suspension system has strong ability to resist the disturbances, which were caused by both the external uncertainty and parametric uncertainty. In other words, the suspension system of the rescue vehicle built in this paper has good stability and robustness.
In addition, in order to better illustrate the adaptability of the proposed switch control strategy, a bumpy road was used as the road input, which is shown in Figure 20a, and the vehicle steering angle of the J-turn maneuver is displayed in Figure 20b. The dynamic responses of the HDSMC-controlled vehicle (blue lines), RCSMC-controlled vehicle (black lines), and the vehicle with the Switched SMC (red lines) are shown in Figure 21.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 25 of 31 As shown in Figure 18, the red and blue curves represent an increase and decrease of 20% of the damping comparing to the black curve of the initial condition, respectively. Obviously, changes in suspension damping in this range have almost no effect on the acceleration and roll angle of the vehicle system. As a result, the suspension system in this paper has a good anti-interference ability to the parameter uncertainty of suspension damping. Figure 19 shows the impact of vehicle speed changes on the suspension system. The vehicle speed under the initial condition is 70 km/h, which is the maximum speed of the rescue vehicle in this paper. Therefore, reducing vehicle speed by 20% and 30%, respectively, was used to verify the interference caused by the uncertainties of vehicle speed on the system. As can be seen from the two diagrams in Figure 19, the system has a strong anti-interference ability to vehicle speed uncertainty.
As a result, it can be seen from the above simulation results and analysis that the suspension system has strong ability to resist the disturbances, which were caused by both the external uncertainty and parametric uncertainty. In other words, the suspension system of the rescue vehicle built in this paper has good stability and robustness.
In addition, in order to better illustrate the adaptability of the proposed switch control strategy, a bumpy road was used as the road input, which is shown in Figure 20a, and the vehicle steering angle of the J-turn maneuver is displayed in Figure 20b. The dynamic responses of the HDSMCcontrolled vehicle (blue lines), RCSMC-controlled vehicle (black lines), and the vehicle with the Switched SMC (red lines) are shown in Figure 21. As shown in Figure 18, the red and blue curves represent an increase and decrease of 20% of the damping comparing to the black curve of the initial condition, respectively. Obviously, changes in suspension damping in this range have almost no effect on the acceleration and roll angle of the vehicle system. As a result, the suspension system in this paper has a good anti-interference ability to the parameter uncertainty of suspension damping. Figure 19 shows the impact of vehicle speed changes on the suspension system. The vehicle speed under the initial condition is 70 km/h, which is the maximum speed of the rescue vehicle in this paper. Therefore, reducing vehicle speed by 20% and 30%, respectively, was used to verify the interference caused by the uncertainties of vehicle speed on the system. As can be seen from the two diagrams in Figure 19, the system has a strong anti-interference ability to vehicle speed uncertainty.
As a result, it can be seen from the above simulation results and analysis that the suspension system has strong ability to resist the disturbances, which were caused by both the external uncertainty and parametric uncertainty. In other words, the suspension system of the rescue vehicle built in this paper has good stability and robustness.
In addition, in order to better illustrate the adaptability of the proposed switch control strategy, a bumpy road was used as the road input, which is shown in Figure 20a, and the vehicle steering angle of the J-turn maneuver is displayed in Figure 20b. The dynamic responses of the HDSMCcontrolled vehicle (blue lines), RCSMC-controlled vehicle (black lines), and the vehicle with the Switched SMC (red lines) are shown in Figure 21. The results for the bump road in Figure 21 are similar to those of Figure 14. In the time range of [0, 3] s, the vehicle vertical acceleration of the suspension system with RCSMC is significantly smaller than HDSMC and Switched SMC. The vehicle vertical acceleration and roll angle become bigger after 3 s because of the input of bump road; however, the vertical acceleration of the suspension system with RCSMC is still smaller than the vehicles with the other two controllers during the time of [3, 6.5] s. One can notice that at time 6.5 s, when the Switched SMC is switched, the performance of roll angle and LTR is better than the RCSMC-controlled vehicle. Consequently, the value of LTR is stabilized to its set value with the control of Switched SMC, whereas the value changes in an unsafe range with RCSMC. Similarly, there are also shocks that are induced by the switch action of the vehicle, which are not obvious as in Figure 14a because of the worse road condition. The frequency domain response of the body acceleration and roll angle are shown in Figure 21c,d, which provide further obvious vindication of the Switched SMC. It can also be seen that the Switched SMC can take care of both sides of ride comfort and handling stability rather than one aspect.

Fishhook Maneuver
The fishhook maneuver is applied to demonstrate the proposed switch control strategy of the hydro-pneumatic suspension, and the vehicle steering angle under the fishhook maneuver is shown in Figure 22b. The random road shown in Figure 21a and the bump road displayed in Figure 23a are used as the road input of simulation, respectively. The results for the bump road in Figure 21 are similar to those of Figure 14. In the time range of [0, 3] s, the vehicle vertical acceleration of the suspension system with RCSMC is significantly smaller than HDSMC and Switched SMC. The vehicle vertical acceleration and roll angle become bigger after 3 s because of the input of bump road; however, the vertical acceleration of the suspension system with RCSMC is still smaller than the vehicles with the other two controllers during the time of [3, 6.5] s. One can notice that at time 6.5 s, when the Switched SMC is switched, the performance of roll angle and LTR is better than the RCSMC-controlled vehicle. Consequently, the value of LTR is stabilized to its set value with the control of Switched SMC, whereas the value changes in an unsafe range with RCSMC. Similarly, there are also shocks that are induced by the switch action of the vehicle, which are not obvious as in Figure 14a because of the worse road condition. The frequency domain response of the body acceleration and roll angle are shown in Figure 21c,d, which provide further obvious vindication of the Switched SMC. It can also be seen that the Switched SMC can take care of both sides of ride comfort and handling stability rather than one aspect.

Fishhook Maneuver
The fishhook maneuver is applied to demonstrate the proposed switch control strategy of the hydro-pneumatic suspension, and the vehicle steering angle under the fishhook maneuver is shown in Figure 22b. The random road shown in Figure 21a  Once the fishhook steering maneuver is exerted at the 5th second in Figure 22b, the value of LTR goes up rapidly and reaches 0.8 at the 5.8th s, and then the controller is switched to drag the LTR to the set value. Until the 9th s, the value of LTR is reduced to less than the set value. Then the controller is switched to pursue the ride comfort of the vehicle. The switch actions are obvious, which can be observed and proved in Figure 23e. The steering angle changes rapidly in the opposite direction in the time rage of [9 9.4] s, but the value of LTR does not rise to the set value, and the control mode does not switch. It is easy to see that the roll angle and the LTR of the vehicle with the Switched SMC become bad and move toward the value of the RCSMC-controlled vehicle, because the Switched SMC has been pursuing the ride comfort when the vehicle is not in danger from the 9th s. The mutation points II and III in Figure 23a are caused by the switch action, and the mutation points I, IV, and V are induced by the steering actions. From the frequency domain response of curves shown in Figure  23c,d, the same conclusion can be reached that the Switched SMC can simultaneously improve the ride comfort and handling stability under this condition compared with the conventional SMC controlled suspension.
The rescue vehicle dynamic responses based on the switch strategy under the bump road and the fishhook maneuver ( Figure 24) are shown in Figure 25. Similarly, the analysis results are the same as those in Figure 23, and the switch actions of the Switched SMC come true according to LTR value.  Once the fishhook steering maneuver is exerted at the 5th second in Figure 22b, the value of LTR goes up rapidly and reaches 0.8 at the 5.8th s, and then the controller is switched to drag the LTR to the set value. Until the 9th s, the value of LTR is reduced to less than the set value. Then the controller is switched to pursue the ride comfort of the vehicle. The switch actions are obvious, which can be observed and proved in Figure 23e. The steering angle changes rapidly in the opposite direction in the time rage of [9 9.4] s, but the value of LTR does not rise to the set value, and the control mode does not switch. It is easy to see that the roll angle and the LTR of the vehicle with the Switched SMC become bad and move toward the value of the RCSMC-controlled vehicle, because the Switched SMC has been pursuing the ride comfort when the vehicle is not in danger from the 9th s. The mutation points II and III in Figure 23a are caused by the switch action, and the mutation points I, IV, and V are induced by the steering actions. From the frequency domain response of curves shown in Figure 23c,d, the same conclusion can be reached that the Switched SMC can simultaneously improve the ride comfort and handling stability under this condition compared with the conventional SMC controlled suspension.
The rescue vehicle dynamic responses based on the switch strategy under the bump road and the fishhook maneuver ( Figure 24) are shown in Figure 25. Similarly, the analysis results are the same as those in Figure 23, and the switch actions of the Switched SMC come true according to LTR value.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 28 of 31 Once the fishhook steering maneuver is exerted at the 5th second in Figure 22b, the value of LTR goes up rapidly and reaches 0.8 at the 5.8th s, and then the controller is switched to drag the LTR to the set value. Until the 9th s, the value of LTR is reduced to less than the set value. Then the controller is switched to pursue the ride comfort of the vehicle. The switch actions are obvious, which can be observed and proved in Figure 23e. The steering angle changes rapidly in the opposite direction in the time rage of [9 9.4] s, but the value of LTR does not rise to the set value, and the control mode does not switch. It is easy to see that the roll angle and the LTR of the vehicle with the Switched SMC become bad and move toward the value of the RCSMC-controlled vehicle, because the Switched SMC has been pursuing the ride comfort when the vehicle is not in danger from the 9th s. The mutation points II and III in Figure 23a are caused by the switch action, and the mutation points I, IV, and V are induced by the steering actions. From the frequency domain response of curves shown in Figure  23c,d, the same conclusion can be reached that the Switched SMC can simultaneously improve the ride comfort and handling stability under this condition compared with the conventional SMC controlled suspension.
The rescue vehicle dynamic responses based on the switch strategy under the bump road and the fishhook maneuver ( Figure 24) are shown in Figure 25. Similarly, the analysis results are the same as those in Figure 23, and the switch actions of the Switched SMC come true according to LTR value.  Based on the analysis of simulation results mentioned above, it is concluded that (1) the PACSA performs better than GA in terms of the parameters optimization of the SMC; (2) the switch strategy proposed in this paper is all effective under different vehicle speeds, road inputs, and the steering modes.

Conclusions
This paper presents a sliding mode switch control strategy for both the ride comfort and handling stability of the adjustable hydro-pneumatic suspension system of a rescue vehicle. The switch index of the Switched SMC is defined according to the value of LTR. Three different SMC Based on the analysis of simulation results mentioned above, it is concluded that (1) the PACSA performs better than GA in terms of the parameters optimization of the SMC; (2) the switch strategy proposed in this paper is all effective under different vehicle speeds, road inputs, and the steering modes.

Conclusions
This paper presents a sliding mode switch control strategy for both the ride comfort and handling stability of the adjustable hydro-pneumatic suspension system of a rescue vehicle. The switch index of the Switched SMC is defined according to the value of LTR. Three different SMC controllers are designed respectively and tested through the simulations of the half-car model based on MATLAB/Simulink. AHP is applied to calculate the weight coefficients of the objective functions, and the parameters of the controller are tuned by PACSA, the optimization ability of which proved better than that of GA. The performance of the switch strategy is proved through simulations with different road input (random road and bump road), different vehicle speeds, and different steering maneuvers (J-turn steering and Fishhook steering).
Compared with the conventional SMC controlled hydro-pneumatic suspension system, the simulation results show that the proposed switch control strategy has the following advantages: (1) The PACSA performs better than a genetic algorithm in terms of parameter optimization of the SMC. (2) The proposed switch control strategy can simultaneously address concerns on both ride comfort and handling stability under different road surfaces and driving conditions through switch actions.
In the future work, on the one hand, it is essential to reduce the shocks induced by the switch action, which will make the switch control strategy work better. On the other hand, we will apply the proposed control strategy based on related algorithms to a real fire rescue prototype to validate the proposed strategy when the prototype is completed in the future.