#### Abstract

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 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.

#### 1. 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–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.

#### 2. 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 bearing, is the mass center of the flywheel, is the distance between 1 and 2, is the distance from the mass center to 1, is the distance from the mass center to 2, and are, respectively, the angular displacements around the -axis and -axis, , , , , 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]where , , , , and are the displacements of the flywheel in five freedom degrees, , , , , and are, respectively, the control current corresponding to each freedom degree, , , , , and are, respectively, the displacement stiffness of each freedom degree, and , , , , and are, respectively, the current stiffness of each freedom degree.

For the convenience of analysis, type (2) can be expressed by type (3):where is the mass matrix, is the damping coefficient matrix, is the displacement stiffness matrix, is the current stiffness matrix, is the unit matrix, is the displacement vector, is the control current vector, and is the unbalanced inertial force vector. Type (3) can be rewritten as

The state space of the flywheel system transfer function matrix iswhere is the state vector, is the vector control, is the measured quantity, namely, the sensor output signal vector of the five freedom degrees. Moreover, type (5) can also be expressed aswhere 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.

#### 3. 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 dual-mode 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.

##### 3.1. 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, is the transfer function matrix of the controller, is the input signal matrix, , , and are the weighted functions, is the state space equation of the maglev flywheel (it contains the power amplifier transfer function , the sensor transfer function , and the control current , and so forth), is the sensor output signal vector, is the system reference signal vector, and is the error signal matrix.

is the closed transfer function in Figure 3 can be written as [11]where type (8) is the linear fractional transformation of the controller . The standard of control problem is to find a real rational , to make the controlled object stable work, and to make the minimal norm of in the whole frequency range [12], should be satisfied bywhere is the biggest singular value of , “sup” is the supremum of in the whole frequency range, and is a given positive number.

The corresponding sensitivity and complementary sensitivity matrix functions and combined with Figure 2 can be deduced: is the weighted function of ; its main purpose is to limit the amplitude of in a specified frequency range. is the weighted function of , whose purpose is to limit the amplitude of . In addition, the transfer function matrix of the controller output is also has a weighted function ; its main purpose is to limit the controller output. 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 kg, N/A, N/A, N/m, and N/m. Based on controller problem solving limit and MATLAB robust control instructions, the solved parameters are involved in the solution procedure, and the weighting functions and are determined, respectively, by

In addition, is selected as through the simulation analysis. Then, the robust performance index (0.4232) is calculated based on the controlled object transfer function matrix 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.

##### 3.2. 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 . 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, is the learning gain of AILC, and is the impacting factor of . 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]

The discrete transfer function of the learning law of (14) can be calculated bywhere 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 bywhere 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 when infinite iteration is operated, and the error signal will become zero. According to (22) and (26), the error signal can be shown as

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 aswhere 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.

#### 4. Simulation and Experimental

##### 4.1. Simulation Analysis

First, the stability of controller is attested through analyzing the singular values relationship between and , as well as between and . 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 , , and , and so on.

The choice principle of is guaranteeing the anti-interference 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 ; if it does not conform to the demand, the weighted function must be chosen again. The choice principle of 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 ; if it does not conform to the demand, the weighted function must be chosen again. Figure 5 shows that the curve of is under the curve of , and the curve of is under the curve of , which demonstrates that the selection of weighted function in types (14) and (15) can meet the design requirements, and the solved controller is appropriate.

Second, the simulation parameters of AILC are selected as Hz, , , , , and the definition of “” indicates that the AILC algorithm has 100 memory points in one control period. Figure 6 shows the rotor simulation orbits and radial run-out displacements with AILC compensation.

**(a) Rotor position orbit**

**(b)**Rotor run-out in axisThe flywheel has a regular circular trajectory in 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 and ensuring the flywheel to rotate around its collection center.

##### 4.2. 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 dual-model 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.

#### 5. 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 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, help to prolong its service life, and accelerate the development scale of electric vehicles.

#### Conflict of Interests

The authors declare that they have no conflict of interests regarding this paper.

#### Acknowledgments

This work is supported by the National Natural Science Foundation of China (51405244), the China Postdoctoral Science Foundation (2014M551634), and the Jiangsu Province Natural Science Foundation of China (BK20140880). The authors would like to thank the editor and the reviewers for their constructive comments and suggestions to improve the quality of the paper.