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

Authors: Hui Gao; Ying-Jun Wu; Jing-Jin Shen
Journal: Mathematical Problems in Engineering
 (2015)

---

## 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 onH∞ 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.


---

## Body

## 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 onH∞ 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 H∞ 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. Figure1 shows the below force of the flywheel.Figure 1
Flywheel below force analysis.In Figure1, O-XYZ is the space coordinate, 1 and 2 are, respectively, the left and right position of the radial magnetic bearing, Og(xg,yg,zg) is the mass center of the flywheel, l is the distance between 1 and 2, l1 is the distance from the mass center to 1, l2 is the distance from the mass center to 2, φ and θ are, respectively, the angular displacements around the x-axis and y-axis, ΔFx1, ΔFx2, ΔFy1, ΔFy2, and ΔFz are the electromagnetic force of the flywheel in the three axes, and fx, fy, and fz 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]:(1)mx¨g=ΔFx1+ΔFx2-fx,my¨g=ΔFy1+ΔFy2-fy,Irφ¨-ωIaθ˙=ΔFy1·l1-ΔFy2·l2,Irθ¨+ωIaφ˙=-ΔFx1·l1+ΔFx2·l2,mz¨g=ΔFz-fz,where m is the flywheel quality, Ir and Ia 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](2)
l
2
l
m
l
1
l
m
0
0
0
0
0
l
2
l
m
l
1
l
m
0
0
0
-
I
r
l
I
r
l
0
-
I
r
l
I
r
l
0
0
0
0
0
0
0
m
x
¨
1
x
¨
2
y
¨
1
y
¨
2
z
¨
+
0
0
0
0
0
0
0
0
0
0
ω
I
a
l
-
ω
I
a
l
0
0
0
0
0
-
ω
I
a
l
ω
I
a
l
0
0
0
0
0
0
x
˙
1
x
˙
2
y
˙
1
y
˙
2
z
˙
=
k
x
x
1
k
x
x
2
0
0
0
0
0
k
x
y
1
k
x
y
2
0
0
0
l
1
k
x
y
1
-
l
2
k
x
y
2
0
-
l
1
k
x
x
1
l
2
k
x
x
2
0
0
0
0
0
0
0
k
x
z
x
1
x
2
y
1
y
2
z
+
k
i
x
1
k
i
x
2
0
0
0
0
0
k
i
y
1
k
i
y
2
0
0
0
l
1
k
i
y
1
-
l
2
k
i
y
2
0
-
l
1
k
i
x
1
l
2
k
i
x
2
0
0
0
0
0
0
0
k
i
z
i
c
x
1
i
c
x
2
i
c
y
1
i
c
y
2
i
c
z
+
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
f
x
f
y
0
0
f
z
,where x1, x2, y1, y2, and z are the displacements of the flywheel in five freedom degrees, icx1, icx2, icy1, icy1, and icz are, respectively, the control current corresponding to each freedom degree, kxx1, kxx2, kxy1, kxy2, and kxz are, respectively, the displacement stiffness of each freedom degree, and kix1, kix2, kiy1, kiy2, and kiz are, respectively, the current stiffness of each freedom degree.For the convenience of analysis, type (2) can be expressed by type (3):(3)Mx¨+Cx˙+Kx=Bic+If,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, ic is the control current vector, and f is the unbalanced inertial force vector. Type (3) can be rewritten as(4)x¨=-CMx˙-KMx+BMic+IMf.The state space of the flywheel system transfer function matrixG(s) is(5)x¯=Ax¯+B1f+B2u,y=Cx¯,where x¯=[xx˙]T is the state vector, u=ic 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(6)
G
s
=
G
11
s
G
12
s
G
21
s
G
22
s
=
A
B
1
B
2
C
O
O
,where O is the five-order zero matrix; it can be solved out through types (4) and (5):(7)A=OI-KM-CM,B1=OIM,B2=O-IM,C=IO.So far, the nonlinear dynamic modelG(s) 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 onH∞ 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 H∞ controller and unbalance compensation controller.Figure 2
The dual-mode switch control strategy diagram.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, aH∞ 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 H∞ control and AILC is implemented by the real-time switching controller by judging the status of the EV and the maglev flywheel.


### 3.1.H∞ Controller Design

A kind ofH∞ control strategy with mixed sensitivity is designed according to the characteristic of the maglev flywheel state space equation; it is shown in Figure 3.Figure 3
The diagram ofH∞ control strategy.In Figure3, KH is the transfer function matrix of the H∞ controller, u is the input signal matrix, W1, W2, and W3 are the weighted functions, G is the state space equation of the maglev flywheel (it contains the power amplifier transfer function GA, the sensor transfer function GS, and the control current ic, 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](8)Tc=Fl(G,KH)=G11+G12KHI-G22KH-1G21,where type (8) is the linear fractional transformation of the H∞ controller KH. The standard of H∞ control problem is to find a real rational KH, to make the controlled object G stable work, and to make the minimal H∞ norm of Tc in the whole frequency range [12], should be satisfied by(9)Tc∞=supωσ-Tcjω<γ,where σ-(Tc(jω)) is the biggest singular value of Tc, “sup” is the supremum of σ-(Tc(jω)) in the whole frequency range, and γ is a given positive number.The corresponding sensitivity and complementary sensitivity matrix functionsS and T combined with Figure 2 can be deduced: (10)S=er=1I+KH(s)G(s),T=yr=KH(s)G(s)I+KH(s)G(s)=I-S.W1(s) is the weighted function of S; its main purpose is to limit the amplitude of S in a specified frequency range. W3(s) is the weighted function of T, whose purpose is to limit the amplitude of T. In addition, the transfer function matrix of the H∞ controller output is(11)R=ur=KH(s)I+KH(s)G(s).R also has a weighted function W2(s); its main purpose is to limit the KH controller output. W2(s) 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 asm=10.8 kg, kix1=kix2=kiy1=kiy2=954.6 N/A,  kiz=1543.8 N/A, kxx1=kxx2=kxy1=kxy2=1.517×106 N/m, and KXZ=3.216×106 N/m. Based on H∞ controller problem solving limit and MATLAB robust control instructions, the solved parameters are involved in the G(s) solution procedure, and the weighting functions W1(s) and W3(s) are determined, respectively, by (12)W1s=98s+2,W3s=1.5ss+3200.In addition,W2(s) is selected as 2.0×107 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 H∞ controller are concluded:(13)KHx1=KHx2=KHy1=KHy2=-21.8z3+15.3z2-25.6z+21.32z3-10.61z2+8.69z-1.73,KHz=-101.3z3+80.5z2+4.68z3-5.87z2+2.32z-0.21.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.Figure 4
AILC compensation principle of AMB 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 vk. β 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 yd and the output of the sensor yk are needed in AILC; here, yd is 2.5 V witch is defined as the balance position of rotor during static suspension. The error signal between yd and yk is iteratively learned, then, the perfect and unknown control signal uk 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(14)vk+1n=βPvkn+Qekn+Q1ekn-1.To understand the action of AILC better, the control input signal of generalized plant is written as(15)ukn=cn+βvkn,where uk(n) is the controller input, c(n) is PID controller output, vk(n) is the learning gain of AILC, and β is the impacting factor of vk(n). The last objective is to obtain perfect controller signal ud when having infinitely iterative learning and then make the error signal become zero. It should be satisfied as [13] (16)limk→∞uk(n)=ud,limk→∞ek(n)=0.The discrete transfer function of the learning law of (14) can be calculated by(17)Vk+1z=βPVkz+QEkz+Q1Ekzz-1=βPVkz+Q+Q1z-1Ekz,where z-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:(18)ct=GPt=Kpet+1Ti∫0tetdt+Td1+Tfdetdt.The discrete function of (18) is deduced as(19)Cz=Kp1+T0Ti11-z-1+Td1+TfT01-z-1Ez.Thez-transform of the control signal function of (15) is deduced by(20)Ukz=Cz+βVkz=kp1+T0Ti11-z-1+Td1+TfT01-z-1Ez+βVkz,where defining(21)Hz=1kp1+T0Ti11-z-1+Td1+TfT01-z-1-1.The discrete function of error signal can be obtained as(22)Ez=Ukz-βVkzHz.Putting (22) into (17), it has(23)Vk+1z=βPVkz+Q+Q1z-1Ukz-βVkzHz=βP-Q+Q1z-1HzVkz+Q+Q1z-1HzUkz.Equation (24) gives the transformation in limit calculation at the two sides of (23)(24)V∞z=limk→∞Vkz=limk→∞Q+Q1z-1HzUkz1-βP+Q+Q1z-1Hz.IfP and β all are 1, at the same time, the following inequality is satisfied by [14]:(25)βP-β[Q+Q1z-1]H(z)∞<1.Equation (24) can be simplified by(26)V∞z=limk→∞Ukz=Udz.That is, the controller signaluk(z) is replaced by βvk(z) when infinite iteration is operated, and the error signal will become zero. According to (22) and (26), the error signal can be shown as(27)limk→∞Ez=Hzlimk→∞Ukz-βlimk→∞Vkz=HzUdz-βV∞z=H(z)[Ud(z)-βUd(z)]=0.The convergence of AILC has been demonstrated according to (26) and (27), and a perfect controller signal ud 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(28)β=ffdn,where f is the flywheel frequency, fd is a given frequency, and the action of n is to reduce the value of β when f is far away from fd (usually, n is greater than 2). In the start-up process, due to f≪fd, the value of β should be very small; it can weaken the influence of Vk on the generalized plant. However, β will be close to 1 if f≈fd; it can enhance the effect of repetitive learning and can make the error signal to converge toward zero.


## 3.1.H∞ Controller Design

A kind ofH∞ control strategy with mixed sensitivity is designed according to the characteristic of the maglev flywheel state space equation; it is shown in Figure 3.Figure 3
The diagram ofH∞ control strategy.In Figure3, KH is the transfer function matrix of the H∞ controller, u is the input signal matrix, W1, W2, and W3 are the weighted functions, G is the state space equation of the maglev flywheel (it contains the power amplifier transfer function GA, the sensor transfer function GS, and the control current ic, 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](8)Tc=Fl(G,KH)=G11+G12KHI-G22KH-1G21,where type (8) is the linear fractional transformation of the H∞ controller KH. The standard of H∞ control problem is to find a real rational KH, to make the controlled object G stable work, and to make the minimal H∞ norm of Tc in the whole frequency range [12], should be satisfied by(9)Tc∞=supωσ-Tcjω<γ,where σ-(Tc(jω)) is the biggest singular value of Tc, “sup” is the supremum of σ-(Tc(jω)) in the whole frequency range, and γ is a given positive number.The corresponding sensitivity and complementary sensitivity matrix functionsS and T combined with Figure 2 can be deduced: (10)S=er=1I+KH(s)G(s),T=yr=KH(s)G(s)I+KH(s)G(s)=I-S.W1(s) is the weighted function of S; its main purpose is to limit the amplitude of S in a specified frequency range. W3(s) is the weighted function of T, whose purpose is to limit the amplitude of T. In addition, the transfer function matrix of the H∞ controller output is(11)R=ur=KH(s)I+KH(s)G(s).R also has a weighted function W2(s); its main purpose is to limit the KH controller output. W2(s) 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 asm=10.8 kg, kix1=kix2=kiy1=kiy2=954.6 N/A,  kiz=1543.8 N/A, kxx1=kxx2=kxy1=kxy2=1.517×106 N/m, and KXZ=3.216×106 N/m. Based on H∞ controller problem solving limit and MATLAB robust control instructions, the solved parameters are involved in the G(s) solution procedure, and the weighting functions W1(s) and W3(s) are determined, respectively, by (12)W1s=98s+2,W3s=1.5ss+3200.In addition,W2(s) is selected as 2.0×107 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 H∞ controller are concluded:(13)KHx1=KHx2=KHy1=KHy2=-21.8z3+15.3z2-25.6z+21.32z3-10.61z2+8.69z-1.73,KHz=-101.3z3+80.5z2+4.68z3-5.87z2+2.32z-0.21.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.Figure 4
AILC compensation principle of AMB 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 vk. β 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 yd and the output of the sensor yk are needed in AILC; here, yd is 2.5 V witch is defined as the balance position of rotor during static suspension. The error signal between yd and yk is iteratively learned, then, the perfect and unknown control signal uk 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(14)vk+1n=βPvkn+Qekn+Q1ekn-1.To understand the action of AILC better, the control input signal of generalized plant is written as(15)ukn=cn+βvkn,where uk(n) is the controller input, c(n) is PID controller output, vk(n) is the learning gain of AILC, and β is the impacting factor of vk(n). The last objective is to obtain perfect controller signal ud when having infinitely iterative learning and then make the error signal become zero. It should be satisfied as [13] (16)limk→∞uk(n)=ud,limk→∞ek(n)=0.The discrete transfer function of the learning law of (14) can be calculated by(17)Vk+1z=βPVkz+QEkz+Q1Ekzz-1=βPVkz+Q+Q1z-1Ekz,where z-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:(18)ct=GPt=Kpet+1Ti∫0tetdt+Td1+Tfdetdt.The discrete function of (18) is deduced as(19)Cz=Kp1+T0Ti11-z-1+Td1+TfT01-z-1Ez.Thez-transform of the control signal function of (15) is deduced by(20)Ukz=Cz+βVkz=kp1+T0Ti11-z-1+Td1+TfT01-z-1Ez+βVkz,where defining(21)Hz=1kp1+T0Ti11-z-1+Td1+TfT01-z-1-1.The discrete function of error signal can be obtained as(22)Ez=Ukz-βVkzHz.Putting (22) into (17), it has(23)Vk+1z=βPVkz+Q+Q1z-1Ukz-βVkzHz=βP-Q+Q1z-1HzVkz+Q+Q1z-1HzUkz.Equation (24) gives the transformation in limit calculation at the two sides of (23)(24)V∞z=limk→∞Vkz=limk→∞Q+Q1z-1HzUkz1-βP+Q+Q1z-1Hz.IfP and β all are 1, at the same time, the following inequality is satisfied by [14]:(25)βP-β[Q+Q1z-1]H(z)∞<1.Equation (24) can be simplified by(26)V∞z=limk→∞Ukz=Udz.That is, the controller signaluk(z) is replaced by βvk(z) when infinite iteration is operated, and the error signal will become zero. According to (22) and (26), the error signal can be shown as(27)limk→∞Ez=Hzlimk→∞Ukz-βlimk→∞Vkz=HzUdz-βV∞z=H(z)[Ud(z)-βUd(z)]=0.The convergence of AILC has been demonstrated according to (26) and (27), and a perfect controller signal ud 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(28)β=ffdn,where f is the flywheel frequency, fd is a given frequency, and the action of n is to reduce the value of β when f is far away from fd (usually, n is greater than 2). In the start-up process, due to f≪fd, the value of β should be very small; it can weaken the influence of Vk on the generalized plant. However, β will be close to 1 if f≈fd; 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 ofH∞ controller is attested through analyzing the singular values relationship between S(s) and W1(s), as well as between T(s) and W3(s). 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 Tc, KH, and G(s), and so on.Figure 5
Sensitivity, complementary sensitivity, and the corresponding weighted function singular value relations.The choice principle ofW1(s) is guaranteeing the anti-interference and tracking ability of the flywheel system, and the smaller the singular value of S(s) is, the better the system tracking ability is. The curve of S(s) should be under the curve of the 1/W1(s); if it does not conform to the demand, the weighted function W1(s) must be chosen again. The choice principle of W3(s) is guaranteeing the flywheel system output to recurrence the input, and the smaller the singular value of T(s) is, the smaller the system impact by compound disturbance because of model uncertainty is. The curve of T(s) should also be under the curve of the 1/W3(s); if it does not conform to the demand, the weighted function W3(s) must be chosen again. Figure 5 shows that the curve of S(s) is under the curve of 1/W1(s), and the curve of T(s) is under the curve of 1/W3(s), which demonstrates that the selection of weighted function in types (14) and (15) can meet the design requirements, and the solved H∞ controller is appropriate.Second, the simulation parameters of AILC are selected asfd=600 Hz, P=0.995, Q=0.75, Q1=0.25, n=0,1,…,99, and the definition of “n” 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.Figure 6
Flywheel trajectory with AILC compensation.
(a)
Rotor position orbit
(b)
Rotor run-out in axisx1The flywheel has a regular circular trajectory in Figure6(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 Figure7, the radial displacement voltage curve and the control current curve of the flywheel are shown in start-up process by the solved H∞ 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 H∞ controller has good robust stability and can be used in the control system of vehicle maglev flywheel.Figure 7
Flywheel displacement and control current in start-up process byH∞ control.Figure8 includes the radial displacement curve, the control current curve, and the speed measuring pulse when flywheel normal is rotating by H∞ 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.Figure 8
Flywheel displacement and control current in rotating process byH∞ control and the speed measuring pulse at 600 Hz.The compensated effect by AILC algorithm at 600 Hz is shown in Figure9, 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 9
Flywheel displacement and control current in AILC compensation process at 600 Hz.Figure10 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.Figure 10
EV primary battery output power in different modes.


## 4.1. Simulation Analysis

First, the stability ofH∞ controller is attested through analyzing the singular values relationship between S(s) and W1(s), as well as between T(s) and W3(s). 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 Tc, KH, and G(s), and so on.Figure 5
Sensitivity, complementary sensitivity, and the corresponding weighted function singular value relations.The choice principle ofW1(s) is guaranteeing the anti-interference and tracking ability of the flywheel system, and the smaller the singular value of S(s) is, the better the system tracking ability is. The curve of S(s) should be under the curve of the 1/W1(s); if it does not conform to the demand, the weighted function W1(s) must be chosen again. The choice principle of W3(s) is guaranteeing the flywheel system output to recurrence the input, and the smaller the singular value of T(s) is, the smaller the system impact by compound disturbance because of model uncertainty is. The curve of T(s) should also be under the curve of the 1/W3(s); if it does not conform to the demand, the weighted function W3(s) must be chosen again. Figure 5 shows that the curve of S(s) is under the curve of 1/W1(s), and the curve of T(s) is under the curve of 1/W3(s), which demonstrates that the selection of weighted function in types (14) and (15) can meet the design requirements, and the solved H∞ controller is appropriate.Second, the simulation parameters of AILC are selected asfd=600 Hz, P=0.995, Q=0.75, Q1=0.25, n=0,1,…,99, and the definition of “n” 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.Figure 6
Flywheel trajectory with AILC compensation.
(a)
Rotor position orbit
(b)
Rotor run-out in axisx1The flywheel has a regular circular trajectory in Figure6(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 Figure7, the radial displacement voltage curve and the control current curve of the flywheel are shown in start-up process by the solved H∞ 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 H∞ controller has good robust stability and can be used in the control system of vehicle maglev flywheel.Figure 7
Flywheel displacement and control current in start-up process byH∞ control.Figure8 includes the radial displacement curve, the control current curve, and the speed measuring pulse when flywheel normal is rotating by H∞ 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.Figure 8
Flywheel displacement and control current in rotating process byH∞ control and the speed measuring pulse at 600 Hz.The compensated effect by AILC algorithm at 600 Hz is shown in Figure9, 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 9
Flywheel displacement and control current in AILC compensation process at 600 Hz.Figure10 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.Figure 10
EV primary battery output power in different modes.


## 5. Conclusion

Through analyzing the nonlinear dynamic characteristics and the application on EV, one kind of dual-model control strategy based onH∞ 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 H∞ 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.



---

*Source: 327347-2015-09-16.xml*