# Fractional-Order Modeling and Analysis of a Variable Structure Hybrid Energy Storage System for EVs

Authors: Jianlin Wang; Dan Xu; Jiahui Zhou; Jinlu Mao
Journal: Complexity
 (2020)

---

## Abstract

Hybrid energy storage system has been widely studied as an important technology for electric vehicles. Since the hybrid energy storage system is a nonlinear and complex system, the modeling of the system and the high-precision nonlinear control strategy are technical difficulties for research. The establishment of a high-precision mathematical model of the hybrid energy storage system is the basis for the study of high-quality nonlinear control algorithms. Fortunately, the theory of fractional calculus can help build accurate mathematical models of hybrid energy storage systems. In order to obtain the high-quality nonlinear control strategy of this complex system, this paper, respectively, carried out fractional-order modeling and analysis on the three basic equivalent working states of the hybrid energy storage system of electric vehicles. Among them, the fractional-order average state space model is carried out for the equivalent Buck and Boost mode. Also, the steady-state analysis of the equivalent Dual-Boost mode is carried out by combining the fractional-order calculus theory with the equivalent small parameter variable method. Finally, the effectiveness and precision of the fractional-order model are proved by simulation and experiment.


---

## Body

## 1. Introduction

The hybrid energy storage system (HESS) of electric vehicles, that is, the composite power supply technology, is a combination of two or more energy sources to form a power system for electric vehicles [1, 2] . Generally, the overall performance of electric vehicles is improved by complementing the advantages of various energy sources, which is a relatively advanced technology in the field of electric vehicles. As we all know, the hybrid energy storage system contains two or more power supplies, and there are also multiple DC-DC converters with strong nonlinear characteristics, which make the design of the control strategy for the hybrid energy storage system difficult. In order to be able to design high-performance control strategies, it is necessary to carry out topology design and high-precision system modeling for hybrid energy storage systems.The topology of a HESS is generally passive, semiactive, active, or cascaded. It is all well known, and each topology of HESS has obvious advantages and disadvantages [3]. In order to integrate the advantages of various topologies, a new variable structure HESS has been proposed, which can realize the function of battery/supercapacitor (SC) semiactive structure, SC/battery semiactive structure, and fully active structure. The topology of the new variable structure HESS is shown in Figure 1. It mainly has three working states; that is, it works in the Buck mode when the SC absorbs the braking energy; it works in the Boost mode when the SC or the battery is outputting separately; and it works in the Dual-Boost mode when the SC and the battery are outputting together.[4, 5].Figure 1
Topology of the variable structure HESS.The control problem of this nonlinear complex system is a difficult point, but there are some research results that can be used. For example, the sliding mode control method is used to control the DC-DC converter to implement power distribution, and the convergence time of the positive Markov neural network is optimized [6–8], and the fractional calculus theory is used to perform fractional control to obtain high performance [9]. However, the accuracy of these control strategies needs to be based on a high-precision system model, and the design of HESS for electric vehicles is diverse so that there are relatively few studies on the modeling of variable structure HESS for electric vehicles.An accurate mathematical model for the equivalent DC-DC converter system corresponding to each working state is the basis for applying the variable structure HESS use in electric vehicles. Considering the nondualities of commercial inductors, it is difficult to describe such circuit components with long-memory characteristics using conventional integer-order calculus. Fortunately, the fractional-order theory is gradually maturing, and in the past decade, fractional-order calculus has been applied in almost all fields of science, engineering, and mathematics. Concepts from fractional-order circuits and systems have attracted much attention from the electrical engineering community. Some research results prove that the fractional-order model has higher precision than the integer-order model in mathematical modeling of nonlinear circuit systems, and the integer-order model is a special case of the fractional-order model [10]. Therefore, this article focuses on the fractional-order modeling research on the three working states of the variable structure HESS.The volt-ampere characteristics of fractional inductor and capacitor shown in Figure1 can be described as(1)vL1=LdαiL1dtα,vL2=LdβiL2dtβ,iC=CdγvCdtγ,where α is the order of the fractional inductor L1, β is the order of the fractional inductor L2, and γ is the order of the fractional capacitor C, respectively [11].In this article, the fractional-order modeling of equivalent DC-DC converters corresponding to the three working states of the variable structure HESS is mainly based on the fractional-order characteristics of inductors and capacitors. The rest of this manuscript is organized as follows: Section2 introduces the fractional-order modeling of the equivalent Buck converter for braking energy recovery. Section 3 introduces the fractional-order modeling of the equivalent Boost converter for power output. Section 4 introduces the steady-state analysis of the dual-Boost output mode based on the fractional-order equivalent small parameter variable method. In Section 5, simulations and experiments are performed to verify the fractional-order mode. Finally, conclusions are given in Section 6.


## 2. Fractional-Order Modeling of Equivalent Buck Converter for Braking Energy Recovery

From the working principle of the variable structure HESS, it is known that when the electric vehicle brakes to produce energy, it is absorbed by the SC alone. In this mode of operation, the equivalent Buck converter of the variable structure HESS is shown in Figure2. The input voltage is equivalent to the voltage of the motor inverter, the load is the internal resistance of the SC, and the capacitance is the equivalent capacitance of the SC.Figure 2
Equivalent Buck converter of the variable structure HESS.The equivalent Buck converter has two main operating states. The operating state 1 is that the switch SW is closed; that is, the capacitor is charged; the operating state 2 is that the switch SW is open; that is, the capacitor is discharged. Throughout one duty cycle period, the inductor voltage can be defined as(2)vL=vin−vo,0<t<dT,−vo,dT<t<T.The average current of the inductor current over the entire period can be described as(3)iL1t,αT=1T∫tt+TiL1τ,αdτ.The average inductor current also can be converted into a form by fractional-order calculation as follows:(4)dαiL1t,αTdtα=1T∫tt+TdαiL1t,αdταdτ=dαiL1t,αdtα,where iL1t,α is the superposition of the AC component under the influence of the DC component and the high-frequency switching harmonics. Similarly, the input voltage, output voltage, inductor current, and switching duty cycle in an equivalent Buck converter can be considered as a superposition of the DC component and the interference component [12, 13].Assume that the equivalent capacitance of the SC also has a fractional-order characteristic, and the order is defined asρ. A fractional-order average state space model can be established for equivalent Buck converter.Whent∈0,dT, the equation of state for the equivalent Buck converter is(5)dαiL1dtα=1L1vin−vo,dρvodtρ=1CsiL1−voRs.Whent∈dT,T, the equation of state for the equivalent Buck converter is(6)dαiL1dtα=−1L1vo,dρvodtρ=1CsiL1−voRs.Introducing the duty cycle function in equations (5) and (6), the fractional-order average state model of the equivalent Buck circuit is obtained, as(7)dαiL1dtαdρvodtρ=0−1L11Cs−1RsCsiL1vo+dL10vin.Based on the fractional-order Laplace transform method [6], the Laplace transform of equation (7) is obtained:(8)sαiL1ssρvos=0−1L11Cs−1RsCsiL1svos+dsL10vins.The transfer function of the output voltage to the input voltage is obtained by equation (8) as(9)GVs=vosvins|d^s=0=D/L1CSsαρ+1/RSCSsα+1/L1CS.Assuming that the input and output voltages are constant during the cycle, we can obtain the input and output energies as(10)Ein=vin∫0dTiL1τ,αdτ,=vindTiR+vin−voL1Γα+1Tα+1ddα−d,Eo=vo∫0dTiL1τ,αdτ+∫dTTiL2τ,αdτ,=voTiR+dα+1−dvin+1−dαL1Γα+2Tα+1.At steady state, the average energy efficiency is(11)η=EoEm=d1−1−dαiRiR−vinTα1−dαdα−dL1Γα+2.


## 3. Fractional-Order Modeling of Equivalent Boost Converter for Power Output

When the power demand of the electric vehicle is positive, the variable structure HESS works in the Boost mode and boosts the voltage of the battery or the SC for output. In this working state, the equivalent Boost converter of the variable structure HESS is shown in Figure3. The input voltage vin is equivalent to the voltage of the battery or SC, and the load is the motor inverter. In this section, the battery boost output is chosen as an example to illustrate the fractional-order modeling and analysis of the Boost mode.Figure 3
Equivalent Boost converter of the variable structure HESS.Similar to the equivalent Buck converter, the equivalent Boost converter also contains two operating states. In one duty cycle, the inductor voltage can be described as(12)vL=vin,0<t<dT,vin−vo,dT<t<T.Whent∈0,dT, the equation of state for the equivalent Boost converter is(13)dβiL2dtβ=1L2vin,dγvodtγ=1RCvo.Whent∈dT,T, the equation of state for the equivalent Boost converter is(14)dβiL2dtβ=1L2vin−vo,dγvodtγ=1CiL2−voR.Introducing the duty cycle function in equations (13) and (14), the fractional-order average state model of the equivalent Buck circuit is obtained, as(15)dβiLdtβdγvodtγ=0−1−dL21−dC−1RCiL2vo+1L20vin.Based on the fractional-order Laplace transform method, the Laplace transform of equation (15) is obtained:(16)sβiL2ssγvos=0−1−dsL21−dsC−1RCiL2svos+1L20vins.The transfer function of the output voltage to the input voltage is obtained by equation (16) as(17)GVs=vosvins|d^s=0=1−D/L2Csβ+γ+1/RCsβ+1−Dβ/LC.Assuming that the input and output voltages are constant during the cycle, we can obtain the input and output energies as(18)Ein=vin∫0dTiL2′τ,βdτ+∫aTTiL2″τ,βdτ,=vinT1−diR−vinTβ1−d+d2−dβ+1−voTβ1−dβ+1L2Γβ+2,Eo=vo∫aTTiL2″τ,βdτ=voiRT=viniRT1−dβ.At steady state, the average energy efficiency is(19)η=EoEin=1−dβ1−d1−vinTβ1−d+d2−dβ+1−vo1−dβ+1L2Γβ+2iR.


## 4. Fractional-Order Modeling and Steady-State Analysis for Dual-Boost Output Mode

When the state of the battery and the SC are not in an extreme state, the variable structure HESS is often operating in a common output mode. In this mode, the switch SW3 is turned off, and the variable structure HESS can be equivalent to the dual-Boost parallel output. The equivalent DC-DC converter is shown in Figure4. Such a parallel output circuit system has high requirements for the design of the control algorithm and the implementation of the power allocation strategy.Figure 4
Equivalent dual-Boost converter of variable structure HESS.The state space averaging has many advantages in modeling, but not enough precision. It cannot perform steady-state analysis of this type of circuitry, especially ripple analysis. The equivalent small parameter variable method is a high-precision symbol analysis method combining the disturbance analysis method and the harmonic balance method [4]. This method does not need to introduce small parameter quantities, which has certain advantages for the ripple analysis of nonlinear high-order systems [14, 15].Considering the fractional-order characteristics of the capacitor and inductor in the circuit system shown in Figure4, the results of the circuit analysis are(20)dαiL1dtα=VSCL1−vCL1,dβiL2dtβ=VBatL2−vCL2,dγvCdtγ=iL1+iL2C−vCR.The state equation of the dual-PWM switch Boost equivalent circuit system is described as(21)G0pαpβpγx+G1f1x+G2f2x=U,where x=iL1iL2vCT, G0=pα01/L10pβ1/L2−1/C−1/Cpγ+1/RC,G1=00−1/L10001/C00, G2=00000−1/L201/C0,U=VSC/L1VBat/L20,pα=Dα is fractional calculus operator, f1x=δ1tx, f2x=δ2tx, and δt is the pulse function.State variablex and δt pulse function are, respectively, expanded into a series form of a sum of a main component and a small component, that is,(22)x=x0+∑i=1∞εixi,(23)δnt=δ0n+∑i=1∞εiδin.Substituting equations (22) and (23) into equation (21), we get(24)fnδn,x=f0n+∑i=1∞εifin,where i is the order of the harmonics and xi, δi, and fi are the i-th order components of x, δx, and fx. ε is a small component symbol, which only means that xi is much smaller than x.The Fourier series ofx0 and xii≥1 are expanded to(25)x0=∑k∈E0xk0=a00,xi=∑k∈Eirxkj=a0i+∑k∈Eirakiejkτ+a¯kie−jkτ,where τ=ωt=2πt/T, a0t is the DC component of xi, and aki is the k-th order amplitude of the i-th harmonic. E0 is the spectrum set of the state vector x0, and Eir is the spectrum set of the correction xi during the iterative calculation.Furthermore, the Fourier series expansion of theδ0n and δii≥1n is(26)δ0n=b0n+b1nejr+b¯1ne−jt,δin=b2inej2iτ+b2i+1nej2i+1τ+b¯2ine−j2iτ+b¯2i+1ne−j2i+1τ,where b0n=1/T∫0Tδntdt=dnandbin=1/2ain−jβin.Furthermore, equation (20) can be rewritten as(27)f0n=δ0nx0,f1n=δ0nx1+δ1nx0,f2n=δ0nx2+δ1nx1+δ2nx0,⋯.Expand the functionfinx into the series form of the sum of the main and small components, which can be expressed as(28)fin=fi=1n+εRi+1,i=0,1,2,…,.Wheni=0, f0mn is the main component of f0n, which is the same as the spectrum of x0, and R1 is the remaining high-order small components. According to the spectrum of R1, the frequency set of can be determined. Similarly, the frequency set of f0mn can be determined according to the spectrum of fimn. f0mn, fimn, and Ri in equation (28) are(29)f0mn=h00+∑k∈E0hΣ0ejkτ+h¯k0e−jkz,fimn=h0i+∑k∈Eirhkiejkr+h¯kie−jkτ,Ri=g0i+∑k∈Eirgkiejkτ+g¯kie−jkτ,where h00,h0i,andg0i are all DC components. Substitute the above equation into equation (28) to get the following equation:(30)fn=f0mn+εf1mn+ε2f2mn+…+εR1+ε2R2+….According to the above analysis, the principal component of the steady-state solution is DC, and the principal component of the state variable is defined as(31)x0=a00=IL100IL200VC00T.Considering that there is no derivative change in the DC component of the system, letG00=G00α0β0γ, and we can get(32)001−d1L1001−d2L2−1−d1C−1−d2C1RC⋅IL1∞IL200VC00=VSCL1VBatL20.Solving equation (32) can obtain the steady-state solution of the main oscillation component, and we can get(33)VC00=VSC1−d1=VBat1−d2,IL,00=VSC2R1−d12,IL200=VBat2R1−d22.The first-order correction contains only the fundamental component, so suppose(34)x1=a11ejτ+a¯11e−jτ,where a11=IL111IL211VC11T and IL111, IL211, and VC11 are the first-order amplitude small components of the first harmonics of iL111, iL211, and vC11, respectively. Substituting equation (34) into equation (30), we get(35)f1m1=b01a11ejτ+a¯11e−jτ,f1m2=b02a11ejτ+a¯11e−jτ,R11=a11b11ejτ+b¯11e−jτ,R12=a11b12ejτ+b¯12e−jτ.According to the calculation principle of fractional calculus,pαejωt=jωαejωt, pβejωt=jωβejωt, and pγejωt=jωγejωt. [16] The first-order correction amount that can be calculated is(36)VC11≈b1jωαL1IL,00+jωβL2IL200−1−d′VC00Δjω,IL11=b1VC00−1−d1VC11jωαL1,IL211=b1VC00−1−d2VC11jωβL2.According to the above calculation principle,Eir=0,2,3 is known; then, we suppose(37)x2=a02+a22ej2τ+a¯22e−j2τ+a32ej3τ+a¯32e−j3τ,where a02=IL102IL202VC02T is the second-order correction of the DC component and a22=IL122IL222VC22T and a32=IL132IL232VC32T are the higher-order components of the second harmonic, respectively.Substituting equation (37) into equation (30), after ignoring the high-order small component of the second harmonic, we can get(38)f2m1≈b01a02+b01a22+b31a¯11ej2τ+b01a32+b11a22+b21a11ej3t+c.c,f2m2≈b02a02+b02a22+b32a¯11ej2t+b02a32+b12a22+b22a11ej3t+c.c.Furthermore, the second-order correction of the DC component and the higher-order component of the second-order harmonic can be obtained:(39)VC02=b¯1VC11+b1V¯C111−d′,IL102=VC022R1−d1,IL202=2VC022R1−d2,VC22=1−db1VC11+b2VC00+b3V¯C11Δj2ω−j2ωαL1b1IL111+b2IL100+b3I¯L111+j2ωβL2b1IL211+b2IL200+b3I¯L211Δj2ω,IL122=b1VC11+b2VC00+b3V¯C11−1−d1VC22j2ωαL1,IL222=b1VC11+b2VC00+b3V¯C11−1−d2VC22j2ωβL2,VC32=1−db1VC22+b2VC11+b3V¯C00Δj3ω−j3ωαL1b1IL122+b2IL111+b3I¯L100+j3ωβL2b1IL222+b2IL211+b3I¯L200Δj3ω,IL132=b1VC22+b2VC11+b3V¯C00−1−d1VC23j3ωαL1,IL232=b1VC22+b2VC11+b3V¯C00−1−d2VC23j3ωβL2.Using the principle of harmonic balance, it is possible to approximatea00,a11,a02,a22,a32. Thereby, a steady-state periodic solution of the system of equation (25) is obtained, which includes a DC solution and a ripple solution, as(40)iL1≈IL100+IL102+2ReIL111cosτ−ImIL111sinτ,+ReIL122cosτ−ImIL122sinτ,+ReIL132cosτ−ImIL132sinτ,iL2≈IL200+IL202+2ReIL211cosτ−ImIL211sinτ,+ReIL222cosτ−ImIL222sinτ,+ReIL232cosτ−ImIL232sinτ,vC≈VC00+VC02+2ReVC11cosτ−ImVC11sinτ,+ReVC22cosτ−ImVC22sinτ,+ReVC32cosτ−ImVC32sinτ.


## 5. Simulation

In order to verify the validity and accuracy of the fractional-order average state space model and the fractional-order equivalent small parameter variable method, simulation studies are implemented in MATLAB2018b [17]. The values of the main parameters of the variable structure HESS L1, L2, C, and R are set to 100μH, 100μH, 470μF, and 10Ω, respectively.According to the fractional average state space model of the equivalent Buck converter and the equivalent Boost converter, the corresponding mathematical simulation models were established in MATLAB2018b/Simulink, as shown in Figure5. In the open-loop simulation of the equivalent buck converter fractional-order model, the input voltage is set to 24V, the duty cycle is set to 0.7, and the order is chosen to be 0.8, 0.9, and 1, respectively. The simulation results are shown in Figure 6(a) and Figure 6(b). In the open-loop simulation of the equivalent Boost converter fractional-order model, the input voltage is set to 24V, the duty cycle is set to 0.5, and the order is chosen to be 0.8, 0.9, and 1, respectively. The simulation results are shown in Figure 6(c) and Figure 6(d). In addition, experiments were carried out on the Buck mode and Boost mode, respectively. It was found that the experiment results of output voltage are close to the simulation results with the order of 0.95. This shows that the inductor and capacitor used in the circuit have fractional-order characteristics, and the order is close to 0.95.Figure 5
Fractional-order simulation of equivalent DC-DC converter: (a) Buck converter; (b) Boost converter.
(a)(b)Figure 6
Simulation results of fractional-order model of equivalent DC-DC converter: (a) output voltage of fractional-order model of equivalent Buck converter; (b) inductor current of fractional-order model of equivalent Buck converter; (c) output voltage of fractional-order model of equivalent Boost converter; (d) inductor current of fractional-order model of equivalent Boost converter.
(a)(b)(c)(d)In order to verify the effectiveness of the fractional-order small parameter variable method for the steady-state analysis of the equivalent dual-switch converter, a simulation program of the fractional-order small parameter variable method was written in MATLAB2018b. The calculation process of programming can be divided into the following 5 parts:①
Establish the fractional-order state equation based on the equivalent circuit and determine the sign row vector.②
Calculate the Fourier series of the single pulse function.③
Calculate the DC componenta00.④
Calculate the ripple componentsa11,a02,a22,a32.⑤
Get the steady-state periodic solution of the equivalent circuit.Write②-⑤ as a subprogram, the main program entry is the circuit parameter matrix and duty cycle row vector, and the exit is the steady-state periodic solution of the circuit. In the simulation process, the SC voltage is set to 38 V, the lithium battery voltage is 24 V, and the duty cycle of SW1 and SW2 is set to 0.4 and 0.5 (synchronously boosted to 48 V output). The calculation results are obtained in the time domain, and the steady-state waveforms that generate the inductor current and output voltage are shown in Figure 7.Figure 7
Simulation result of fractional-order model of equivalent dual-switch Boost converter: (a) inductor L1 current; (b) inductor L2 current; (c) output voltage.
(a)(b)(c)From the simulation results, it can be seen that the fractional-order model of the equivalent DC-DC converter has more freedom to choose when describing the DC-DC operating characteristics. As the fractional-order decreases, the overshoot of the dynamic characteristics decreases. The amount of ripple in the steady-state characteristic increases. In practical applications, the voltage response of the equivalent Boost converter of the variable structure HESS has a certain degree of ripple component. However, the ripple of the steady-state characteristics described by the integer-order model is very small, which is obviously in error with the actual situation, and it also shows that the accuracy of integer-order modeling is insufficient. That is to say, when the order of the fractional model is selected to an appropriate value, the dynamic and steady-state characteristics of the DC-DC converter can be accurately described.


## 6. Conclusions

As an important technology for electric vehicles, hybrid energy storage systems have been extensively studied. However, there are still many problems, and this technology cannot be made mature and widely used. Especially the high-precision modeling and control strategy for this complex system is the bottleneck of this technology. In this paper, a fractional-order modeling of a variable structure HESS is carried out. For the three working modes of variable structure HESS, different fractional modeling methods are used to perform fractional modeling and analysis. The effectiveness and precision of the fractional-order model are proved by simulation and experiment.



---

*Source: 7643812-2020-12-09.xml*