Research on Adaptive Dual-Mode Switch Control Strategy for Vehicle Maglev Flywheel Battery

Because of the jamming signal is real-time changeable and control algorithm cannot timely tracking control flywheel rotor, this paper takes vehicle maglev flywheel battery as the research object. One kind of dual-model control strategy is developed based on the analysis of the vibration response impact of the flywheel battery control system. In view of the complex foundation vibration problems of electric vehicles, the nonlinear dynamic simulation model of vehicle maglev flywheel battery is solved. Through analyzing the nonlinear vibration response characteristics, one kind of dual-mode adaptive hybrid control strategy based on H∞ control and unbalance displacement feed-forward compensation control is presented and a real-time switch controller is designed. The reliable hybrid control is implemented, and the stability in the process of real-time switch is solved. The results of this project can provide important basic theory support for the research of vehicle maglev flywheel battery control system.


Introduction
As the future main traffic tools, electric vehicle (EV) is required in the performance of starting, acceleration, and climbing; however, this performance depends largely on the power battery performance [1].But, the faults of costly, short range and short service life become the bottleneck to restrict EV development scale.Specifically in the process of frequent start-stop or climbing, the chemical battery life is more shorten because of fast and deep discharging [2,3].How to eliminate these shortcomings becomes the key EV to be quickly developed.
Maglev flywheel can be applied in EV electric power system, aerospace, and other fields because of having high specific energy, high power, fast charge and discharge, long service life, no waste gas pollution, environment-friendly advantages, and so on [4][5][6].In the field of EV, the maglev flywheel either can be as an independent power driving EV [3], or can be used as auxiliary power assisting the motive power batteries work [7,8].However, the maglev flywheel control stability will be affected because of existing start-stop, acceleration and deceleration, steering, and road random vibration in the process of the EV driving and even cause instability.So, the efficiency of magnetic suspension flywheel must be reduced.
For the scientific research and practical application of maglev flywheel, a dual-mode adaptive hybrid control strategy is studied based on ∞ control and AILC algorithms, and the state space equation of maglev flywheel was analyzed.To improve the robust stability of flywheel control system and reduce the real-time interference, one ∞ controller based on the state space equation was solved.To reduce maglev flywheel radial run-out, one adaptive iterative learning control theory was deduced and unbalance displacement compensation was implemented.The experimental result shows the dual-model control strategy has better interference capability in maglev flywheel start-up process and has stronger active control ability relative to only PID control.The maglev flywheel based on the dual-model control can help the EV primary battery, improve its discharge characteristics, and help to prolong its service life.

Maglev Flywheel Nonlinear Dynamic Model
The below force and movement differential equation of the flywheel are discussed in order to solve the maglev flywheel nonlinear dynamic model.Figure 1 shows the below force of the flywheel.
In Figure 1, - is the space coordinate, 1 and 2 are, respectively, the left and right position of the radial magnetic Figure 1: Flywheel below force analysis.
bearing,   (  ,   ,   ) is the mass center of the flywheel,  is the distance between 1 and 2,  1 is the distance from the mass center to 1,  2 is the distance from the mass center to 2,  and  are, respectively, the angular displacements around the -axis and -axis, Δ 1 , Δ 2 , Δ 1 , Δ 2 , and Δ  are the electromagnetic force of the flywheel in the three axes, and   ,   , and   are the disturbing force and the unbalance force of the flywheel in the three axes.The flywheel movement differential equations are deduced according to the below force analysis in Figure 1 and the motion laws of particles [9]: where  is the flywheel quality,   and   are, respectively, the radial inertia moment and axial inertia moment, and  is the angular velocity of the flywheel.The motion differential equation of matrix form can be written by [10] ( 1 0 0 0 0 0 1 0 0 0 0 0 1 0 0 0 0 0 1 0 0 0 0 0 1 where  1 , where M is the mass matrix, C is the damping coefficient matrix, K is the displacement stiffness matrix, B is the current stiffness matrix, I is the unit matrix, x is the displacement vector, i c is the control current vector, and f is the unbalanced inertial force vector.Type (3) can be rewritten as  The state space of the flywheel system transfer function matrix () is where x = [x ẋ ] T is the state vector, u = i c is the vector control, y = x is the measured quantity, namely, the sensor output signal vector of the five freedom degrees.Moreover, type (5) can also be expressed as where  is the five-order zero matrix; it can be solved out through types ( 4) and ( 5): So far, the nonlinear dynamic model () of the maglev flywheel is solved, and the precise model is provided for the design of the following dual-mode switch controller.

Dual-Mode Switch Control Strategy
To improve the control stability and the energy storage density of vehicle maglev flywheel, a dual-mode switch control strategy is study based on ∞ control algorithm and adaptive iterative learning control algorithm.The control strategy is shown in Figure 2. The control strategy includes generalized control elephant and dual-mode switch controller.Among them, the generalized control elephant is composed of power amplifiers, active magnetic bearing (AMB), flywheel rotor, and displacement and the dualmode switch controller is composed of ∞ controller and unbalance compensation controller.
The basic response is a random signal when EV is driving on the bumpy road or start-stop, acceleration-deceleration, and steering.To reduce the random impact and improve the control robustness, a ∞ controller is solved.Besides, to limit the flywheel radial run-out of the maglev flywheel battery in the charging process and improve the energy storage density of the maglev flywheel, an unbalance displacement feed-forward compensation controller based on adaptive iterative learning control (AILC) algorithm is adopted [10].The switching work between ∞ control and AILC is implemented by the real-time switching controller by judging the status of the EV and the maglev flywheel.

𝐻∞ Controller Design.
A kind of ∞ control strategy with mixed sensitivity is designed according to the characteristic of the maglev flywheel state space equation; it is shown in Figure 3.
In Figure 3, K H is the transfer function matrix of the ∞ controller, u is the input signal matrix,  1 ,  2 , and  3 are the weighted functions, G is the state space equation of the maglev flywheel (it contains the power amplifier transfer function   , the sensor transfer function   , and the control current i c , and so forth), y is the sensor output signal vector, r is the system reference signal vector, and e is the error signal matrix.
T c is the closed transfer function in Figure 3 can be written as [11]  where type ( 8) is the linear fractional transformation of the ∞ controller K  .The standard of ∞ control problem is to find a real rational K  , to make the controlled object G stable work, and to make the minimal ∞ norm of T c in the whole frequency range [12], should be satisfied by where (T c ()) is the biggest singular value of T c , "sup" is the supremum of (T c ()) in the whole frequency range, and  is a given positive number.
The corresponding sensitivity and complementary sensitivity matrix functions S and T combined with Figure 2 can be deduced: 1 () is the weighted function of S; its main purpose is to limit the amplitude of S in a specified frequency range. 3 () is the weighted function of T, whose purpose is to limit the amplitude of T. In addition, the transfer function matrix of the ∞ controller output is R also has a weighted function  2 (); its main purpose is to limit the K  controller output. 2 () should be chosen as a relatively small value in order to reduce the compensators and shorten operation cycle.The structural parameters of the maglev flywheel are calculated and got as  = 10.8 kg,  1 =  2 =  1 =  2 = 954.6N/A,   = 1543.8N/A,  1 =  2 =  1 =  2 = 1.517 × 10 6 N/m, and   = 3.216 × 10 6 N/m.Based on ∞ controller problem solving limit and MATLAB robust control instructions, the solved parameters are involved in the G() solution procedure, and the weighting functions  1 () and  3 () are determined, respectively, by In addition,  2 () is selected as 2.0 × 10 7 through the simulation analysis.Then, the robust performance index (0.4232) is calculated based on the controlled object transfer function matrix G and the three weighted functions; if the index greater than 1, the three weighted functions need to be selected.Finally, the four radial and one axial discrete transfer functions of the ∞ controller are concluded: Because the maglev flywheel is axially symmetric distribution in four radial freedom degrees, the four radial discrete transfer functions are the same.In order to make the control cycle consistent with the actual flywheel control system, the simulation sampling frequency is selected as 20000 Hz.

Feed-Forward Compensation Control Analysis.
To reduce the radial run-out of the flywheel and solve the problem of variable speed, the AILC algorithm is adopted as the feed-forward controller to implement vibratory displacement compensation [10].In Figure 4, the control scheme consists of PID feedback control system, AILC feed-forward compensation controller, and generalized plant.The PID controller can steady the whole system and improve the anti-interference ability, and the action of AILC is to make the learning gain accurately track the expectant orbit.The generalized plant includes power amplifier, electromagnetic coils, and rotor system.
To improve control performance and enhance the convergence rate of the learning law, there are two modifications in AILC.The first one is enhancing the error information action of previous control period and the second one is proposing a novel impacting factor  as the coefficient of the learning gain V  . can reduce the effects of learning gain to the control system when rotor speed does not coincide with the learning cycle of AILC.AILC can implement vibratory displacement compensation without any information of the generalized plant, and it will not increase the interference of the feedback controller.Only the expectant signal   and the output of the sensor   are needed in AILC; here,   is 2.5 V witch is defined as the balance position of rotor during static suspension.The error signal between   and   is iteratively learned, then, the perfect and unknown control signal   is obtained as the input signal of the power amplifier.The functions of AILC can be introduced by the iterative formulas in discrete domain.The error formula is given.
The update learning law of AILC is summarized as To understand the action of AILC better, the control input signal of generalized plant is written as where   () is the controller input, () is PID controller output, V  () is the learning gain of AILC, and  is the impacting factor of V  ().The last objective is to obtain perfect controller signal   when having infinitely iterative learning and then make the error signal become zero.It should be satisfied as [13] lim The discrete transfer function of the learning law of ( 14) can be calculated by where  −1 represents the lag operator in time domain, and it could make the sampled signal lag one period.This is the reason why the former error information is added into the modified learning law of AILC.
The PID controller with incomplete differential part is given in time domain: The discrete function of ( 18) is deduced as The -transform of the control signal function of ( 15) is deduced by where defining The discrete function of error signal can be obtained as Putting ( 22) into (17), it has Equation (24) gives the transformation in limit calculation at the two sides of (23) If  and  all are 1, at the same time, the following inequality is satisfied by [14]: Equation ( 24) can be simplified by That is, the controller signal   () is replaced by V  () when infinite iteration is operated, and the error signal will become zero.According to ( 22) and ( 26), the error signal can be shown as lim The convergence of AILC has been demonstrated according to (26) and ( 27), and a perfect controller signal   has been obtained as power amplifier input.Therefore, the AILC algorithm as the feed-forward compensation controller can be adopted in the application of maglev flywheel unbalance vibratory compensation.
According to the start-up time and the variability of the rotor frequency from static suspension to one fixed speed, the equation of  is given as where  is the flywheel frequency,   is a given frequency, and the action of  is to reduce the value of  when  is far away from   (usually,  is greater than 2).In the start-up process, due to  ≪   , the value of  should be very small; it can weaken the influence of   on the generalized plant.However,  will be close to 1 if  ≈   ; it can enhance the effect of repetitive learning and can make the error signal to converge toward zero.

Simulation and Experimental
4.1.Simulation Analysis.First, the stability of ∞ controller is attested through analyzing the singular values relationship between () and  1 (), as well as between () and  3 ().
The sensitivity and complementary sensitivity function and the corresponding weighted function singular value relations are shown in Figure 5.The curves in Figure 5 are created in the "M" file of MATLAB through some command functions, the transfer functions of T c , K  , and G(), and so on.The choice principle of  1 () is guaranteeing the antiinterference and tracking ability of the flywheel system, and the smaller the singular value of () is, the better the system tracking ability is.The curve of () should be under the curve of the 1/ 1 (); if it does not conform to the demand, the weighted function  1 () must be chosen again.The choice principle of  3 () is guaranteeing the flywheel system output to recurrence the input, and the smaller the singular value of () is, the smaller the system impact by compound disturbance because of model uncertainty is.The curve of () should also be under the curve of the 1/ 3 (); if it does not conform to the demand, the weighted function  3 () must be chosen again.Figure 5 shows that the curve of () is under the curve of 1/ 1 (), and the curve of () is under the curve of 1/ 3 (), which demonstrates that the selection of weighted function in types ( 14) and (15) can meet the design requirements, and the solved ∞ controller is appropriate.
The flywheel has a regular circular trajectory in    and ensuring the flywheel to rotate around its collection center.

Experimental Results
. In Figure 7, the radial displacement voltage curve and the control current curve of the flywheel are shown in start-up process by the solved ∞ controller control; the displacement voltage quickly drops to the equilibrium position (2.5 V) from the sensor calibration position (4 V).It illustrates that the ∞ controller has good robust stability and can be used in the control system of vehicle maglev flywheel.Figure 8 includes the radial displacement curve, the control current curve, and the speed measuring pulse when flywheel normal is rotating by ∞ control and flywheel rotational frequency is 600 Hz.The displacement curve is a sine wave, and it can show the flywheel has mass unbalance.
To restrict flywheel radial run-out, the control current also is sine having basic consistent phase with the displacement voltage.The speed measuring pulse has a voltage range from 0 V to 3.3 V; the purpose is to guarantee the pulse be captured by DSP acquisition circuit.The compensated effect by AILC algorithm at 600 Hz is shown in Figure 9, including the radial displacement curve and the control current curve after compensation.Compared with the radial sine displacement in Figure 8, the radial displacement is balance in the position of 2.5 V, which indicates the flywheel is limited rotating around its geometric center.To limit the radial run-out and increase the active control effect, the control current amplitude is significantly larger than without compensation.Figure 10 gives the output power curves of EV primary battery in three different modes.If flywheel is not working, the output power is about 35 kW when EV starts up; if the flywheel battery is working and only PID controlling, the output power of the primary battery is about 25 kW in EV start-up process; if the dual-model switch control strategy is acting on the flywheel, the output power reduces to 23 kW.
The power change shows the flywheel has better auxiliary effect to EV primary battery when it is controlled by the dualmodel switch control strategy.Similarly, the output power of EV primary battery is the smallest in the acceleration process when the dual-model switch control strategy controls the flywheel.Moreover, because the flywheel rotate speed is higher when the dual-model switch control strategy is working, EV needs lower motive power when restarting.Therefore, when maglev flywheel participates in discharge and is controlled by dual-model switch control strategy, EV primary battery output power is obviously decreased in the whole driving process, the influence on chemical characteristics of primary battery is reduced because the instantaneous discharge depth is small, and the service life of EV primary battery can be improved.

Conclusion
Through analyzing the nonlinear dynamic characteristics and the application on EV, one kind of dual-model control strategy based on ∞ control and AILC algorithms has been studied in this paper, and the state space equation of maglev flywheel was analyzed.To improve the robust stability of flywheel control system and reduce the realtime interference, one ∞ controller based on the state space equation was solved.To reduce maglev flywheel radial run-out, one adaptive iterative learning control theory was deduced, and unbalance displacement compensation was implemented.The experimental result shows the dual-model control strategy has better interference capability in maglev flywheel start-up process and has stronger active control ability relative to only PID control.The maglev flywheel based on the dual-model control can help the EV primary battery, improve its discharge characteristics, help to prolong its service life, and accelerate the development scale of electric vehicles.

Figure 2 :
Figure 2: The dual-mode switch control strategy diagram.

Figure 3 :
Figure 3: The diagram of ∞ control strategy.

Figure 5 :
Figure 5: Sensitivity, complementary sensitivity, and the corresponding weighted function singular value relations.
Figure 6(a) and a sine radial run-out displacement in Figure 6(b) without adding AILC compensation.When the AILC algorithm starts to work at 0.05 s, the circular trajectory is gradual convergence to a point and the amplitude of the radial run-out displacement is obvious attenuation.The curves change simulation result in Figure 6 testifies the AILC algorithm having good displacement compensation effect Rotor run-out in axis  1

Figure 7 :
Figure 7: Flywheel displacement and control current in start-up process by ∞ control.

Figure 8 :
Figure 8: Flywheel displacement and control current in rotating process by ∞ control and the speed measuring pulse at 600 Hz.

Figure 9 :
Figure 9: Flywheel displacement and control current in AILC compensation process at 600 Hz.

Figure 10 :
Figure 10: EV primary battery output power in different modes.
2 ,  1 ,  2 , and  are the displacements of the flywheel in five freedom degrees,  1 ,  2 ,  1 ,  1 , and   are, respectively, the control current corresponding to each freedom degree,  1 ,  2 ,  1 ,  2 , and   are, respectively, the displacement stiffness of each freedom degree, and  1 ,  2 ,  1 ,  2 , and   are, respectively, the current stiffness of each freedom degree.