Control method for hybrid electromagnetic suspension

Abstract

A control method for hybrid electromagnetic suspension. The method provides four modes for hybrid electromagnetic suspension: a comfort mode, a sport mode, a combined mode, and an energy feedback mode. A driver can switch between the four modes as desired. For the comfort, sport, and combined modes, hybrid control is adopted, and two sub-modes are provided: an active control mode and a semi-active control mode. A switching condition between the two sub-modes is determined by using a novel parameter C.sub.act and comparing the same against a maximum equivalent electromagnetic damping coefficient C.sub.eqmax of a linear motor. The present invention solves the problem of achieving a balance between suspension comfort and tire traction, and meets the demands of different operating conditions and users by enabling manual mode switching. In addition, the hybrid control is employed to solve the problems of high energy consumption of active suspension and limited control performance of semi-active suspension, thus ensuring good kinematic performance of automobile suspension while reducing energy consumption. Furthermore, the energy feedback mode is designed to enable the suspension to perform energy recovery, meeting demands of energy conservation and emission reduction.

Claims

1. A control method of an electromagnetic hybrid suspension, comprising the following steps: Step A: providing a quarter electromagnetic hybrid suspension system model, established and expressed by a dynamic differential equation: $\begin{array}{cc}\{\begin{array}{c}{\stackrel{.Math.}{Z}}_b{m}_b=-k_s\left(Z_b-Z_w\right)-C_s\left({\stackrel{.}{Z}}_b-{\stackrel{.}{Z}}_w\right)+F_{\mathrm{act}}\\ {\stackrel{.Math.}{Z}}_w{m}_w=k_s\left(Z_b-Z_w\right)+C_s\left(Z_b-Z_w\right)-k_t\left(Z_w-Z_0\right)-F_{\mathrm{act}}\end{array}& \left(1\right)\end{array}$ in the formula, m.sub.b is spring-loaded mass, m.sub.w is non-spring-loaded mass, k.sub.s is spring stiffness, C.sub.s is passive adjustable damper damping coefficient, k.sub.t is tire equivalent stiffness, Z.sub.0 is road input displacement, and the mathematical model can be expressed as:
Ż.sub.0(t)=−2pf.sub.0Z.sub.0(t)+2p√{square root over (G.sub.0u)}w(t)  (2) in the formula, G.sub.0 is the road roughness, P is π, u is the vehicle speed, f.sub.0 is the lower cut-off frequency, w(t) is the Gauss white noise and the mean value is 0, which is used to generate random road surface; Step B: dividing the electromagnetic hybrid suspension into four working modes based on the quarter electromagnetic hybrid suspension model: a comfort mode, a sport mode, a comprehensive mode and an energy-feeding mode, each with a main improvement target; while the electromagnetic hybrid suspension is: in comfort mode, a Root Mean Square (“RMS”) value of vehicle body acceleration is the main improvement target; in sport mode, a RMS value of tire dynamic load is taken as the main improvement target; in the comprehensive mode, combined RMS values of the vehicle body acceleration and tire dynamic load are the main improvement target; in the energy-feeding mode, energy-feeding is maximized while guaranteeing the basic dynamic performance; Step C: while the electromagnetic hybrid suspension is in comfort mode, sport mode and comprehensive mode, the electromagnetic hybrid suspension comprises two sub-modes: active control mode and semi-active control mode: in the active control mode, calculating a required control force F.sub.des with an Electric Control Unit (“ECU”) through a control strategy, then dividing the required control force F.sub.des by a thrust coefficient K.sub.f of a linear motor to get a current I.sub.des needed by the linear motor, thereafter, the current I.sub.des is input to the ECU; wherein an on-board power supply is connected with the ECU, and the ECU is connected with the linear motor; while the current I.sub.des is input to the linear motor, the linear motor outputs an active control force F.sub.act to the electromagnetic hybrid suspension, at this time, the electric hybrid suspension consumes energy; in the semi-active control mode, calculating the required control force F.sub.des with the ECU through a control strategy, then dividing the required control force F.sub.des by a thrust coefficient K.sub.f of the linear motor to get the current I.sub.des needed by the linear motor, thereafter, the current I.sub.des is input to the ECU; the ECU is connected with the semi-active control loop, the duty cycle is input to the ECU to adjust the current I.sub.des in the ECU, thereby adjusting the equivalent damping coefficient C.sub.eq of the linear motor; a semi-active control loop is connected with the linear motor; the linear motor outputs an electromagnetic damping force to the electromagnetic hybrid suspension system; the linear motor is used as an electromagnetic damper generator, and electric energy recovered by the linear motor acting as a generator is stored in a plurality of super capacitors; introducing parameters: $\begin{array}{cc}{C}_{\mathrm{act}}=-\frac{F_{\mathrm{des}}}{v_{\mathrm{rel}}}& \left(3\right)\end{array}$ in the formula, F.sub.des represents the required control force to be provided by the linear motor and calculated by the control strategy; comparing C.sub.act with a maximum equivalent damping coefficient C.sub.eqmax of the linear motor, if 0<Cact<Ce.sub.eqmax, the required control force F.sub.des is in semi-active control mode, and if C.sub.act>C.sub.eqmax, the required control force F.sub.des is other cases-are in active control mode; Step D: causing the linear motor to act as a generator while the electromagnetic hybrid suspension is in the energy-feeding mode, which is equivalent to a nonadjustable electromagnetic damper acting in the electromagnetic hybrid suspension, and the electromagnetic damping coefficient is a maximum equivalent electromagnetic damping coefficient C.sub.eqmax to ensure that as much electric energy as possible can be recovered, the electric energy is stored in the plurality of super capacitors.

2. The control method of electromagnetic hybrid suspension according to claim 1, wherein the control force F.sub.des in step C is calculated from a Linear Quadratic Gaussian (“LQG”) control strategy; changing the selection of LQG weighting coefficients corresponds to the comfort mode, sport mode and comprehensive mode respectively; a LQG performance index function is designed as follows: $\begin{array}{cc}J=\underset{T.fwdarw.\infty }{\lim }\frac{1}{T}{\int }_0^T\left[{q_1\left(Z_w-Z_0\right)}^2+{q_2\left(Z_b-Z_w\right)}^2+q_3{\stackrel{.Math.}{Z}}_b^2\right]\mathrm{dt}& \left(4\right)\end{array}$ in the formula, Z.sub.w-Z.sub.0 is the dynamic displacement of the tire, wherein multiplying Z.sub.w-Z.sub.0 by an equivalent stiffness k.sub.t of the tire is the dynamic load of the tire; .sub.b-Z.sub.w is the suspension dynamic travel; Z.sub.b is the acceleration of the vehicle body; q.sub.1, q.sub.2 and q.sub.3 represent the weighting coefficients of tire dynamic load, suspension dynamic travel and vehicle body acceleration, respectively.

3. The control method of electromagnetic hybrid suspension according to claim 2, wherein, a specific process of F.sub.des calculation by the LQG control strategy is as follows: a state space equation of electromagnetic hybrid suspension is as follows:
{dot over (X)}=AX+BU  (5) in the formula, selecting X=[Ż.sub.b Ż.sub.w Z.sub.b Z.sub.w Z.sub.0].sup.T as the state variable, U=[F.sub.act] as the input control variable, the following can be obtained: $A=\left[\begin{array}{ccccc}-\frac{C_s}{m_b}& \frac{C_s}{m_b}& -\frac{k_s}{m_b}& \frac{k_s}{m_b}& 0\\ \frac{C_s}{m_w}& -\frac{C_s}{m_w}& \frac{k_s}{m_w}& -\frac{k_s+k_t}{m_w}& \frac{k_t}{m_w}\\ 1& 0& 0& 0& 0\\ 0& 1& 0& 0& 0\end{array}\right]B={\left[\begin{array}{cccc}\frac{1}{m_b}& -\frac{1}{m_w}& 0& 0\end{array}\right]}^T$ selecting Y=[Z.sub.w-Z.sub.0 Z.sub.b-Z.sub.w {umlaut over (Z)}.sub.b].sup.T as the output control variable, the output equation is as follows: $\begin{array}{cc}Y=\mathrm{CX}+\mathrm{DU}C=\left[\begin{array}{ccccc}0& 0& 0& 1& -1\\ 0& 0& 1& -1& 0\\ -\frac{C_s}{m_b}& \frac{C_s}{m_b}& -\frac{k_s}{m_b}& \frac{k_s}{m_b}& 0\end{array}\right]D={\left[\begin{array}{ccc}0& 0& \frac{1}{m_b}\end{array}\right]}^T& \left(6\right)\end{array}$ for the LQG control strategy, the objective is to find the input control variable U=[F.sub.act] so that a performance index function takes the minimum value; the input control variable U belongs to one of the quadratic optimal control problems, so the performance index function is written in the standard quadratic form as follows: $\begin{array}{cc}J=\underset{T.fwdarw.\infty }{\lim }\frac{1}{T}{\int }_0^T\left(X^T\mathrm{QX}+U^T\mathrm{RU}+2X^T\mathrm{NU}\right)\mathrm{dt}& \left(7\right)\end{array}$ the performance index can be obtained from formula (4) and formula (6): $\begin{array}{cc}J=\underset{T.fwdarw.\infty }{\lim }\frac{1}{T}{\int }_0^T{Y}^T{Q}_{0}Y\mathrm{dt}& \left(8\right)\end{array}$ in the formula, Q.sub.0 is a LQG weighting coefficient matrix, which is expressed as follows: $Q_0=\left[\begin{array}{ccc}{q}_1& 0& 0\\ 0& q_2& 0\\ 0& 0& q_3\end{array}\right]$ therefore, the weighting matrix Q of the state variable, the weighting matrix R of the control variable and the weighting matrix N of the crossover term can be obtained: $\begin{array}{c}Q=C^T{Q}_{0}C\\ =\left[\begin{array}{ccccc}\frac{C_s^2}{m_b^2}{q}_3& -\frac{C_s^2}{m_b^2}{q}_3& \frac{C_s{k}_s}{m_b^2}{q}_3& -\frac{C_s{k}_s}{m_b^2}{q}_3& 0\\ -\frac{C_s^2}{m_b^2}{q}_3& \frac{C_s^2}{m_b^2}{q}_3& -\frac{C_s{k}_s}{m_b^2}{q}_3& \frac{C_s{k}_s}{m_b^2}{q}_3& 0\\ \frac{C_s{k}_s}{m_b^2}{q}_3& -\frac{C_s{k}_s}{m_b^2}{q}_3& q_2+\frac{k_s^2}{m_b^2}{q}_3& -q_2-\frac{k_s^2}{m_b^2}{q}_3& 0\\ -\frac{C_s{k}_s}{m_b^2}{q}_3& \frac{C_s{k}_s}{m_b^2}{q}_3& -q_2-\frac{k_s^2}{m_b^2}{q}_3& q_1+q_2+\frac{k_s^2}{m_b^2}{q}_3& -q_1\\ 0& 0& 0& -q_1& q_1\end{array}\right]\end{array}R=D^T{Q}_{0}D=\frac{q_3}{m_b^2}\begin{array}{c}N=C^T{Q}_{0}D\\ D^T{Q}_{0}C\\ ={\left[\begin{array}{ccccc}-\frac{C_s}{m_b^2}{q}_3& \frac{C_s}{m_b^2}{q}_3& -\frac{k_s}{m_b^2}{q}_3& \frac{k_s}{m_b^2}{q}_3& 0\end{array}\right]}^T\end{array}$ the size of the state variable matrix Q is related to the values of q.sub.1, q.sub.2 and q.sub.3, so the weighting coefficient plays a decisive role in the control effect of the LQG control strategy; to obtain an optimal control force F.sub.act of the LQG control strategy, the feedback state gain matrix K is first obtained, namely:
F.sub.act=−KY  (9) the optimal control force F.sub.act can also be expressed as:
F.sub.act=−(k.sub.jŻ.sub.b+k.sub.2Ż.sub.w+k.sub.3Z.sub.b+k.sub.4Z.sub.w+k.sub.5Z.sub.0)  (10) according to LQG control theory, gain matrix K can be expressed as:
K=R.sup.−1(SB+N).sup.T  (11) matrix S is the solution of Ricatti equation, the form of Ricatti equation is as follows:
(SA).sup.T+SA−(SB+N)R.sup.−1(SB+N).sup.T+Q=0  (12) from matrix K, S can be obtained by using LQR function:
[K S E]=lqr(A,B,Q,R,N)  (13).

4. The control method of electromagnetic hybrid suspension according to claim 3, wherein, three state variables of vehicle body acceleration, tire dynamic load and suspension dynamic travel are collected and input into a LQG controller; wherein different mixing modes have different combinations of the LQG weighting coefficients, the LQG controller calculates F.sub.des based on the three state variables and the LQG weighting coefficients.

5. The control method of electromagnetic hybrid suspension according to claim 2, wherein, the weighting coefficients of the LQG control strategy are obtained by a genetic algorithm optimization, an optimization objective of the genetic algorithm optimization are a vehicle body acceleration and a tire dynamic load, which represent a vibration isolation of vehicle and a tire road holding, respectively; a performance parameter of suspension dynamic travel is used as a constraint of the optimization algorithm; the optimized design variable is a LQG weighting coefficient:
x=(q.sub.1 q.sub.2 q.sub.3).sup.T  (14) in the process of optimizing the LQG weighting coefficients by the genetic algorithm optimization, a fitness function needs to be calculated by the objective function; the electromagnetic hybrid suspension has three modes while using the LQG control strategy: comfort, sport and comprehensive, in which comfort mode requires a minimum vehicle body acceleration, sport mode requires a minimum tire dynamic load, and comprehensive mode requires both the minimum vehicle body acceleration and the minimum tire dynamic load; therefore, the calculation methods of the fitness function are different when the genetic algorithm optimization is carried out under three modes; in comfort mode, the objective function is the RMS value of vehicle body acceleration; since the root mean square of a vehicle body acceleration is required to take the minimum value, the fitness function should be proportional to a genetic probability, so the fitness function is required to take the maximum value, therefore, the fitness function is designed as follows: $\begin{array}{cc}F\left(x\right)=\{\begin{array}{cc}{A}_P-A_H,& \mathrm{if}{A}_H<A_P\\ 0,& \mathrm{if}{A}_H\ge A_P\end{array}& \left(15\right)\end{array}$ in the formula, A.sub.P is the RMS value of vehicle body acceleration of traditional passive suspension under the same simulation conditions, A.sub.H is the RMS value of vehicle body acceleration of the electromagnetic hybrid suspension under the LQG control strategy; in the sport mode, the fitness function is designed as follows: $\begin{array}{cc}F\left(x\right)=\{\begin{array}{cc}{T}_P-T_H,& \mathrm{if}{T}_H<T_P\\ 0,& \mathrm{if}{T}_H\ge T_P\end{array}& \left(16\right)\end{array}$ in the formula, T.sub.P is the RMS value of a tire dynamic load of a traditional passive suspension under the same simulation conditions, TH is the RMS value of the tire dynamic load of electromagnetic hybrid suspension under the LQG control strategy; in the comprehensive mode, an objective function is the RMS value of a vehicle body acceleration and a tire dynamic load, so the objective function belongs to a multi-objective genetic algorithm optimization, and the fitness function is designed as follows: $\begin{array}{cc}F\left(x\right)=\frac{A_P}{A_H}+\frac{T_P}{T_H}.& \left(17\right)\end{array}$

6. The control method of electromagnetic hybrid suspension according to claim 5, wherein, the specific steps of genetic algorithm to optimize the weighting coefficients of the LQG control strategy are as follows: (1) parameter initialization: q.sub.1, q.sub.2 and q.sub.3 are encoded by a binary system respectively, and then they are connected in series to form a chromosome string to form an individual; while decoding, the binary code is converted to a decimal value, then the electromagnetic hybrid suspension system is controlled; (2) generating population: an initial population with size of M=80 is generated by a random method, while a new population is generated by a genetic operator operation on the previous generation population; (3) fitness function evaluation: objective functions A.sub.H and T.sub.H of the electromagnetic suspension system controlled by a decoded weighting coefficients are extracted and then transformed into a fitness function F(x); (4) selection operation: using proportional selection method, assuming that the fitness of the i.sup.th individual is F(i), then a genetic probability p.sub.si is: $\begin{array}{cc}\begin{array}{cc}{P}_{si}=F\left(i\right)/\underset{i=1}{\overset{M}{.Math.}}F\left(i\right)& i=1,2,.Math.,M\end{array}& \left(18\right)\end{array}$ (5) cross operation: single point crossover is adopted; select a single crossover point for each pair of the paired individuals according to a crossover probability p.sub.c=0.6 and an exchange some of their chromosomes; (6) mutation operation: a basic bit mutation is used; according to the mutation probability p.sub.m=0.001, a part of a gene value in the coding string of an individual chromosome is inversely calculated; (7) termination condition judgment: a termination algebra T=300 is set, and an optimal individual in a 300.sup.th evolutionary algebraic population is output as an optimal solution; if an average fitness difference of 15 successive generations is less than 0.5%, the operation will be terminated ahead of time, and the optimal individual in a current population will be output as the optimal solution.

7. The control method of electromagnetic hybrid suspension according to claim 1, wherein, the electromagnetic hybrid suspension switching between the active control mode and the semi-active control mode in step C is realized by a first switching device comprising an insulated-gate bipolar transistor (“IGBT”): the on-board power supply is connected with an inverter, and a first end of the inverter is connected with a second switching device; a second end of the second switching device is connected with a third switching device and the linear motor; the third switching device is connected with a rectifier bridge and a DC/DC converter; the output of the DC/DC converter is connected with a positive end of a diode and a negative end of the diode is connected with the super capacitor.

8. The control method of electromagnetic hybrid suspension according to claim 1, wherein, the electromagnetic hybrid suspension used in the energy-feeding mode in step D is as follows: the linear motor is connected with a rectifier bridge, and a first end of the rectifier bridge is connected with a first end of the inductance; a second end of the inductance is connected with a first of a switch and a first end of the switch respectively; a second end of the switch is connected with a positive end of a diode, and a negative end of the diode is connected with a first end of the super capacitor; a second end of the super capacitor is connected with a second end of the switch and a second end of the rectifier bridge respectively.

9. The control method of electromagnetic hybrid suspension according to claim 1, wherein, the adjustable dampers in the comfort mode, sport mode, comprehensive mode and energy-feeding mode correspond to four different gears respectively.

Description

DESCRIPTION OF DRAWINGS

(1) FIG. 1 is an overall flowchart of the present invention;

(2) FIG. 2 is a flowchart of genetic algorithm optimization;

(3) FIG. 3 is a block diagram of the active control sub-mode system;

(4) FIG. 4 is a block diagram of the semi-active control sub-mode system;

(5) FIG. 5 is a topological graph of the semi-active control loop;

(6) FIG. 6 is a block diagram of the switching circuit;

(7) FIG. 7 shows the comparison of time domain simulation of instantaneous power of motor between general active LQG control and hybrid LQG control;

(8) FIG. 8 shows the energy consumption of active suspension and electromagnetic hybrid suspension within 10 seconds;

(9) FIG. 9 is a charging circuit diagram of energy feeding mode.

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS

(10) The present invention is further described in connection with the drawings and examples.

(11) An overall flowchart of a control method of the electromagnetic hybrid suspension of the present invention is shown in FIG. 1.

(12) First, a quarter electromagnetic hybrid suspension system model is established. The linear mathematical model can be expressed by dynamic differential equation:

(13) $\begin{array}{cc}\{\begin{array}{c}{\stackrel{.Math.}{Z}}_b{m}_b=-k_s\left(Z_b-Z_w\right)-C_s\left({\stackrel{.}{Z}}_b-{\stackrel{.}{Z}}_w\right)+F_{\mathrm{act}}\\ {\stackrel{.Math.}{Z}}_w{m}_w=k_s\left(Z_b-Z_w\right)+C_s\left({\stackrel{.}{Z}}_b-{\stackrel{.}{Z}}_w\right)-k_t\left(Z_w-Z_0\right)-F_{\mathrm{act}}\end{array}& \left(1\right)\end{array}$

(14) In the formula, m.sub.b is spring-loaded mass, m.sub.w is non-spring-loaded mass, k.sub.s is spring stiffness, C.sub.s is passive adjustable damper damping coefficient, k.sub.t is tire equivalent stiffness, Z.sub.0 is road input displacement, and the mathematical model can be expressed as:
Ż.sub.0(t)=−2pf.sub.0Z.sub.0(t)+2p√{square root over (G.sub.0u)}w(t)  (2)

(15) In the formula, G.sub.0 is the road roughness, u is the vehicle speed, f.sub.0 is the lower cut-off frequency, w(t) is the Gauss white noise and the mean value is 0, which is used to generate random road surface. In the present invention, the pavement grade is set to C, road roughness G.sub.0 is 256.Math.10.sup.−6 m.sup.3, and the vehicle speed is set to 20 m/s.

(16) In the present invention, electromagnetic hybrid suspension can be divided into four modes: comfort mode, sport mode, comprehensive mode and energy-feeding mode. The switch between the four modes is decided by the driver. In each mode, the damping coefficient of the adjustable damper is different. The adjustable damper is designed to be adjustable among four gears. When the driver switches modes, the controller controls the stepping motor to adjust the throttle size of the damper, thereby adjusting the damping coefficient C.sub.s of the damper.

(17) When the electromagnetic hybrid suspension is in comfort mode, sport mode and comprehensive mode, the suspension is controlled by LQG control strategy. The LQG performance index function is designed as follows:

(18) $\begin{array}{cc}J=\underset{T.fwdarw.\infty }{\lim }\frac{1}{T}{\int }_0^T\left[{q_1\left(Z_w-Z_0\right)}^2+{q_2\left(Z_b-Z_w\right)}^2+q_3{\stackrel{.Math.}{Z}}_b^2\right]\mathrm{dt}& \left(3\right)\end{array}$

(19) In the formula, q.sub.1, q.sub.2 and q.sub.3 represent the weighting coefficients of tire dynamic load, suspension dynamic travel and vehicle body acceleration, respectively. Different weighting coefficients have different effects on suspension dynamic performance. The requirements of the present invention are as follows: the improvement of vehicle body acceleration is emphasized in comfort mode, the improvement of tire dynamic load is emphasized in sport mode, the consideration of vehicle body acceleration and tire dynamic load is given in comprehensive mode, and the index of the suspension dynamic travel plays a limiting role. The genetic algorithm is used for optimization, and the optimization objectives are the vehicle body acceleration and tire dynamic load, and the suspension dynamic travel is taken as the optimum constraint condition. The optimization process is shown in FIG. 2. The weighting coefficients of the three modes after optimization are as follows:

(20) TABLE-US-00001 Working mode q.sub.1 q.sub.2 q.sub.3 Comfort mode 30146 4838 4.85 Sport mode 920358 65037 0.93 Comprehensive mode 41369 4603 1.08

(21) Because the active control of suspension has the disadvantage of high energy consumption and the semi-active control has the disadvantage of limited control effect, the present invention proposes a hybrid control method, which designs two sub-modes: active control mode and semi-active control mode when the electromagnetic hybrid suspension is in comfort mode, sport mode and comprehensive mode.

(22) In the active control mode, on-board power supply powers the linear motor and actively output the control force, which actually controls the input current I of the linear motor, the double-loop control system is designed in the present invention, as shown in FIG. 3. The state feedback variables measured in the suspension system are input to the external controller. The control force F.sub.des required to be provided by the motor is calculated by the controller, and the required current I.sub.des is input to the inner loop control. The inner loop control is the current tracking control, which enables the input current I of the linear motor to track the required current I.sub.des, thereby outputting the control force F.sub.act. External energy comes from direct power supply of on-board power supply.

(23) The control system in semi-active control mode is shown in FIG. 4. According to the feedback state variables in the suspension system, the required equivalent electromagnetic damping force F.sub.des is calculated, and the required current I.sub.des in the winding of linear motor at this time is calculated. The semi-active loop controller adjusts the duty cycle in the semi-active control loop so that the actual current in the winding of the linear motor I.sub.real tracks I.sub.des, so that the equivalent electromagnetic damping force F.sub.act provided by the linear motor tracks F.sub.des. FIG. 5 shows the topological structure of the semi-active control loop. The linear motor M is connected with the rectifier bridge, and the rectifier bridge is connected with the drain of the switch S1. The source of switch S1 is connected to the inductance L and the negative end of the diode D1 respectively. The other end of the inductance L is connected with the drain of the switch S2 and the positive end of the diode D2, and the negative end of the diode D2 is connected with the supercapacitor SC. The other end of the supercapacitor SC is connected with the positive end of the diode D1 and the source of the switch S2. When the terminal voltage of linear motor is less than that of supercapacitor, S1 switches on and S2 switches off, which makes the circuit generate current and avoids dead-zone phenomenon. In other cases, S1 switches off, S2 switches on, and electric energy is stored in the supercapacitor. The upward and downward movement of suspension will make the direction of current in the linear motor change constantly, so it is necessary to add the rectifier bridge in the loop to ensure the consistency of current direction. DC-DC converter plays the role of boosting and lowering voltage. After calculation, the semi-active loop controller outputs two pulse signals to S1 and S2 to adjust duty cycle. In semi-active control, the electric energy generated by linear motor can be stored in the supercapacitor through the semi-active control loop to recover energy.

(24) Introducing parameters:

(25) $\begin{array}{cc}{C}_{\mathrm{act}}=-\frac{F_{\mathrm{des}}}{v_{\mathrm{rel}}}& \left(4\right)\end{array}$

(26) In the formula, F.sub.des represents the control force required to be provided by the linear motor and calculated by the control strategy. Comparing C.sub.act with the maximum equivalent damping coefficient C.sub.eqmax of linear motor, if 0<C.sub.act<C.sub.eqmax, it is in semi-active control mode, and other cases are in active control mode. Switching between the two sub-modes is realized by switching device IGBT. The frequency of switching devices is much higher than that of suspension vibration, as shown in FIG. 6.

(27) The specific process of F.sub.des calculation by LQG control strategy is as follows:

(28) The state space equation of electromagnetic hybrid suspension is as follows:
{dot over (X)}=AX+BU  (5)

(29) Select X=[Ż.sub.b Ż.sub.w Z.sub.b Z.sub.w Z.sub.0].sup.T as the state variable, U=[F.sub.act] as the input control variables, the following can be obtained:

(30) $A=\left[\begin{array}{ccccc}-\frac{C_s}{m_b}& \frac{C_s}{m_b}& -\frac{k_s}{m_b}& \frac{k_s}{m_b}& 0\\ \frac{C_s}{m_b}& -\frac{C_s}{m_b}& \frac{k_s}{m_w}& -\frac{k_s+k_t}{m_w}& \frac{k_t}{m_w}\\ 1& 0& 0& 0& 0\\ 0& 1& 0& 0& 0\end{array}\right]B={\left[\begin{array}{cccc}\frac{1}{m_b}& -\frac{1}{m_w}& 0& 0\end{array}\right]}^T$

(31) Select Y=[Z.sub.w-Z.sub.0 Z.sub.b-Z.sub.w {umlaut over (Z)}.sub.b].sup.T as the output control variables, the output equation is as follows:

(32) $\begin{array}{cc}Y=\mathrm{CX}+\mathrm{DU}C=\left[\begin{array}{ccccc}0& 0& 0& 1& -1\\ 0& 0& 1& -1& 0\\ -\frac{C_s}{m_b}& \frac{C_s}{m_b}& -\frac{k_s}{m_b}& \frac{k_s}{m_b}& 0\end{array}\right]D={\left[\begin{array}{ccc}0& 0& \frac{1}{m_b}\end{array}\right]}^T& \left(6\right)\end{array}$

(33) For LQG control algorithm, the objective is to obtain the input control variables U=[F.sub.act] and make the performance index function take the minimum value. It belongs to one of the quadratic optimal control problems, so the performance index function is written in the standard quadratic form as follows:

(34) $\begin{array}{cc}J=\underset{T.fwdarw.\infty }{\lim }\frac{1}{T}{\int }_0^T\left(X^T\mathrm{QX}+U^T\mathrm{RU}+2X^T\mathrm{NU}\right)\mathrm{dt}& \left(7\right)\end{array}$

(35) It can be obtained from formula (4) and formula (6):

(36) 0$\begin{array}{cc}J=\underset{T.fwdarw.\infty }{\lim }\frac{1}{T}{\int }_0^T{Y}^T{Q}_{0}Y\mathrm{dt}& \left(8\right)\end{array}$

(37) In the formula, Q.sub.0 is LQG weighting coefficient matrix, which is expressed as follows:

(38) $Q_0=\left[\begin{array}{ccc}{q}_1& 0& 0\\ 0& q_2& 0\\ 0& 0& q_3\end{array}\right]$

(39) Therefore, the weighting matrix Q of the state variable, the weighting matrix R of the control variable and the weighting matrix N of the crossover term can be obtained:

(40) $\begin{array}{c}Q=C^T{Q}_{0}C\\ =\left[\begin{array}{ccccc}\frac{C_s^2}{m_b^2}{q}_3& -\frac{C_s^2}{m_b^2}{q}_3& \frac{C_s{k}_s}{m_b^2}{q}_3& -\frac{C_s{k}_s}{m_b^2}{q}_3& 0\\ -\frac{C_s^2}{m_b^2}{q}_3& \frac{C_s^2}{m_b^2}{q}_3& -\frac{C_s{k}_s}{m_b^2}{q}_3& \frac{C_s{k}_s}{m_b^2}{q}_3& 0\\ \frac{C_s{k}_s}{m_b^2}{q}_3& -\frac{C_s{k}_s}{m_b^2}{q}_3& q_2+\frac{k_s^2}{m_b^2}{q}_3& -q_2+\frac{k_s^2}{m_b^2}{q}_3& 0\\ -\frac{C_s{k}_s}{m_b^2}{q}_3& \frac{C_s{k}_s}{m_b^2}{q}_3& -q_2+\frac{k_s^2}{m_b^2}{q}_3& q_1+q_2+\frac{k_s^2}{m_b^2}{q}_3& -q_1\\ 0& 0& 0& -q_1& q_1\end{array}\right]\end{array}R=D^T{Q}_{0}D=\frac{q_3}{m_b^2}\begin{array}{c}N=C^T{Q}_{0}D\\ =D^T{Q}_{0}C\\ ={\left[\begin{array}{ccccc}-\frac{C_s}{m_b^2}{q}_3& \frac{C_s}{m_b^2}{q}_3& -\frac{k_s}{m_b^2}{q}_3& \frac{k_s}{m_b^2}{q}_3& 0\end{array}\right]}^T\end{array}$

(41) It can be concluded that the size of the state variable matrix Q is related to the values of q.sub.1, q.sub.2 and q.sub.3, so the weighting coefficient plays a decisive role in the control effect of the LQG control strategy.

(42) To obtain the optimal control force F.sub.act of LQG control strategy, it is required to obtain the feedback state gain matrix K, namely:
F.sub.act=−KY  (9)

(43) It can also be expressed as:
F.sub.act=−(k.sub.jŻ.sub.b+k.sub.2Ż.sub.w+k.sub.3Z.sub.b+k.sub.4Z.sub.w+k.sub.5Z.sub.0)  (10)

(44) According to LQG control theory, gain matrix K can be expressed as:
K=R.sup.−1(SB+N).sup.T  (11)

(45) Wherein, matrix S is the solution of Ricatti equation. The form of Ricatti equation is as follows:
(SA).sup.T+SA−(SB+N)R.sup.−1(SB+N).sup.T+Q=0  (12)

(46) In software Matlab, matrix K, S can be obtained by using LQR function.
[K S E]=lqr(A,B,Q,R,N)  (13)

(47) Weighting coefficients of LQG control strategy are obtained by genetic algorithm optimization. The optimization objectives of genetic algorithm are vehicle body acceleration and tire dynamic load, which represent vibration isolation of vehicles and tire road holding. The performance parameter of suspension dynamic travel is used as the constraints of the optimization algorithm. The optimized design variable is LQG weighting coefficient:
x=(q.sub.1 q.sub.2 q.sub.3).sup.T  (14)

(48) In the process of optimizing LQG weighting coefficients by genetic algorithm, the fitness function needs to be obtained by calculation through the objective function. The electromagnetic hybrid suspension has three modes when using LQG control strategy: comfort, sport and comprehensive, in which comfort mode requires minimum vehicle body acceleration, sport mode requires minimum tire dynamic load, and comprehensive mode requires both. Therefore, the calculation methods of fitness function are different when genetic algorithm optimization is carried out under three modes.

(49) In comfort mode, the objective function is the RMS value of vehicle body acceleration. Since the root mean square of vehicle body acceleration is required to take the minimum value, while the fitness function should be proportional to the genetic probability, so the fitness function is required to take the maximum value. Therefore, the fitness function is designed as follows:

(50) $\begin{array}{cc}F\left(x\right)=\{\begin{array}{cc}{A}_P-A_H,& \mathrm{if}{A}_H<A_P\\ 0,& \mathrm{if}{A}_H\ge A_P\end{array}& \left(15\right)\end{array}$

(51) In the formula, A.sub.P is the RMS value of vehicle body acceleration of traditional passive suspension under the same simulation conditions. A.sub.H is the RMS value of vehicle body acceleration of electromagnetic hybrid suspension under LQG control strategy. In the sport mode, the fitness function is designed as follows:

(52) $\begin{array}{cc}F\left(x\right)=\{\begin{array}{cc}{T}_P-T_H,& \mathrm{if}{T}_H<T_P\\ 0,& \mathrm{if}{T}_H\ge T_P\end{array}& \left(16\right)\end{array}$

(53) In the formula, T.sub.P is the RMS value of tire dynamic load of traditional passive suspension under the same simulation conditions. T.sub.H is the RMS value of tire dynamic load under LQG control strategy of electromagnetic hybrid suspension.

(54) In the comprehensive mode, the objective function is the RMS value of vehicle body acceleration and tire dynamic load, so it belongs to multi-objective genetic algorithm optimization, and the fitness function is designed as follows:

(55) $\begin{array}{cc}F\left(x\right)=\frac{A_P}{A_H}+\frac{T_P}{T_H}& \left(17\right)\end{array}$

(56) The specific steps of genetic algorithm to optimize the weighting coefficients of LQG control strategy are as follows:

(57) (1) Parameter initialization: q.sub.1, q.sub.2 and q.sub.3 are encoded by binary system respectively, and then they are connected in series to form a chromosome string to form an individual. When decoding, the binary code is converted to decimal value, then the suspension system is controlled;

(58) (2) Generating population: the initial population with size of M=80 is generated by random method, while the new population is generated by genetic operator operation on the previous generation population;

(59) (3) Fitness function evaluation: the objective functions A.sub.H and T.sub.H of the suspension system controlled by the decoded weighting coefficients are extracted and then transformed into fitness function F(x);

(60) (4) Selection operation: using proportional selection method, assuming that the fitness of the i.sup.th individual is F(i), then the genetic probability p.sub.si is:

(61) $\begin{array}{cc}\begin{array}{cc}{p}_{\mathrm{si}}=F\left(i\right)/\underset{i=1}{\overset{M}{.Math.}}F\left(i\right)& i=1,2,\end{array}.Math.,M& \left(18\right)\end{array}$

(62) (5) Cross operation: single point crossover is adopted. Select a single crossover point for each pair of the paired individuals according to the crossover probability p.sub.c=0.6 and exchange some of their chromosomes.

(63) (6) Mutation operation: basic bit mutation is used. According to the mutation probability p.sub.m=0.001, a part of the gene value in the coding string of an individual chromosome is inversely calculated.

(64) (7) Termination condition judgment: the termination algebra T=300 is set, and the optimal individual in the 300.sup.th evolutionary algebraic population is output as the optimal solution. If the average fitness difference of 15 successive generations is less than 0.5%, the operation will be terminated ahead of time, and the optimal individual in the current population will be output as the optimal solution.

(65) The advantage of hybrid control method is that the energy consumption is obviously reduced than active control suspension, and the control effect is similar to active suspension. The comparison of time domain simulation of instantaneous power of motor between general active LQG control and hybrid LQG control is shown in FIG. 7. The analysis shows that the instantaneous power of the two motors is the same for most of the time, and all the power in active control is the power consumed by the motor, while the part below line 0 in hybrid control is semi-active control, which is the power recovered by the motor. FIG. 8 shows the energy consumption of the motor in 10 seconds. It can be seen that the energy consumption of hybrid control is reduced by about 20% compared with that of active control.

(66) When the electromagnetic hybrid suspension is in the energy-feeding mode, the linear motor is required to recover as much vibration energy as possible. Therefore, the linear motor is used as an electromagnetic damper, and the equivalent damping coefficient is the maximum equivalent damping coefficient C.sub.eqmax of the linear motor. The charging circuit of energy recovery in energy-feeding mode is shown in FIG. 9. When the output voltage U.sub.m of linear motor is less than supercapacitor terminal voltage U.sub.c, there will be no induction current in the circuit, and there will be dead-zone phenomenon, which will affect the suspension performance. Therefore, at this time, control switch S1 to switch on and S2 to switch off, and the recovered vibration energy is stored in the inductor L. When the output voltage U.sub.m of linear motor is larger than the supercapacitor terminal voltage U.sub.c, control S1 to switch off and S2 to switch on, charge the supercapacitor at this time.

(67) It should be understood that the above examples are used only to illustrate the present invention and not to limit the scope of the present invention. After reading the present invention, the modifications of various equivalent forms of the invention by skilled persons in the art fall within the scope of the claims appended to the present application.