Next Article in Journal
Speed-Sensorless Control of Induction Machines with LC Filter for Geothermal Electric Submersible Pumping Systems
Next Article in Special Issue
Investigation of Eddy Current Loss and Structure Design with Magnetic-Thermal Coupling for Toothless BLDC High-Speed PM Motor
Previous Article in Journal
A 3D-Printed Continuously Variable Transmission for an Electric Vehicle Prototype
Previous Article in Special Issue
Design and Analysis of a Novel Double-Stator Double-Rotor Motor Drive System for In-Wheel Direct Drive of Electric Vehicles
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Parameter Matching Methods for Li Battery–Supercapacitor Hybrid Energy Storage Systems in Electric Buses

1
College of Electrical and Electronic Engineering, Harbin University of Science and Technology, Harbin 150080, China
2
School of Electrical Engineering and Automation, Harbin Institute of Technology, Harbin 150001, China
3
Harbin Electrical Machinery Company Limited, Harbin 150040, China
*
Author to whom correspondence should be addressed.
Machines 2022, 10(2), 85; https://doi.org/10.3390/machines10020085
Submission received: 9 December 2021 / Revised: 12 January 2022 / Accepted: 18 January 2022 / Published: 25 January 2022
(This article belongs to the Special Issue Electric Vehicle Powertrains: Design, Development and Technology)

Abstract

:
The parameter matching of composite energy storage systems will affect the realization of control strategy. In this study, the effective energy and power utilizations of an energy storage source were defined. With the miniaturization of a composite energy storage system as the optimization goal, the linear programming simplex method was employed to obtain the optimized masses of Li batteries and supercapacitors under the constraints of maximum speed, climbing gradient, acceleration time and cost-effectiveness. As the module numbers shall be integers, the matching results were modified in combination with the graphical method. Owing to the influences of parameter matching schemes on the overall performance and battery life, the critical points of constraints were analyzed and the most appropriate matching numerical points of the composite energy storage system were identified. Simulation and experimental analysis were conducted under practical urban road conditions in China. The results show that the proposed matching method delivers accurate results. Compared with conventional electric buses, the mileage and overall performance of the prototype bus are improved.

1. Introduction

The parameter design of hybrid energy storage systems (HESS) includes power capacity and energy capacity. Currently, the formula method and accurate simulation method have been widely applied. According to the power performance indexes of vehicles and power and energy requirements under typical conditions, the energy and power that the energy storage source system is supposed to provide were separated based on the designed electric transmission system. Then, the specific scheme was designed according to specific energy, specific power, cost, temperature range, reliability and safety of current energy storage source and optimized with a parameter of HESS as the optimization objective. Since the 1980s, studies of composite energy storage sources have been proposed and applied in the energy systems of earth satellites. With the development of battery and supercapacitors, researchers have committed to the matching and control strategy of composite energy storage sources consisting of supercapacitor packs and theoretical and practical explorations of applications of supercapacitors in electric vehicles. Tehrani et al. published a series of articles on supercapacitors and composite energy storage sources and summarized the principle and performance of various supercapacitor energy storages. It was pointed out that supercapacitors are suitable for cases with large power demand but short time [1,2,3,4]. Li et al. investigated simulation and control of composite energy storage sources in Advanced Vehicle Simulator (ADVISOR) and proposed general methods for parameter matching of composite energy storage sources [5,6,7]. Giannoutsos et al. achieved power distribution of the battery and supercapacitor in a composite energy source using low-pass filtering [8,9]. Composite energy storage sources with supercapacitors have been investigated [10,11]. Cao et al. connected DC/DC with a supercapacitor and in parallel with a battery, followed by parameter matching using the formula method. The results showed that the composite energy storage system can effectively provide driving power for motors [12,13]. Shen et al. used the DIRECT (dividing rectangles) method to optimize the multiobjective function based on the total mass of the HESS and the battery capacity. The simulation results showed that the mass was reduced by 66 kg after optimization [14,15]. Hu et al. investigated the impacts of price fluctuation of a single energy storage component on the overall performance of vehicles using the convex optimization method, with cost as the optimization objective [16,17]. Araujo et al. combined parameter matching and energy management based on filtering and nonlinear optimization. The results showed that the total cost was reduced by 20% and the energy utilization was improved by 7.8%. Ravey et al. applied a statistical method to parameter matching of HESSs and achieved good results [18,19].
In China, the methods and principles of parameter matching of composite energy storage systems have been investigated [20,21]. Hong Kong university proposed optimization of fuel cost-effectiveness of electric vehicles by adjusting the ratio of the total mass of energy storage source to the vehicle mass, as well as the ratio of the mass of high specific power energy storage source to the total mass of composite energy storage source, without degradations in mileage and other performances. Harbin Institute of Technology proposed matching strategies based on optimization of the mass ratio and mixing ratio of electric vehicles. A model was established to quantitatively analyze the impact of the mass ratio and mixing ratio on the overall performance of vehicles [22]. Beijing Institute of Technology proposed that other parameters can be determined based on the time constant of the supercapacitor. The number of cells in parallel and maximum output power of the battery pack were optimized with the physical parameters of the composite energy storage source and the charge/discharge rate of the battery as the optimization objectives. The results showed that the acceleration time of the electric bus from 0 to 70 km/h was reduced by 9 s, while the fluctuation range of the output current of the battery decreased significantly [23]. Chongqing University proposed linear optimization of HESS parameters from the perspective of the overall performance of the vehicle. The results showed that compared with the conventional electric vehicles, the battery was protected and the power consumption per 100 km was reduced by 8.55% [24]. Jilin University proposed parameter optimization based on the measured data of vehicles, with charging/discharging capacity and service life of the HESS energy storage component as the optimization objectives [25,26].
In this study, parameter optimization was achieved using the linear programming method from the perspective of power and energy utilizations. The minimum mass of the composite energy storage system was regarded as the optimization objective, while the minimum mass of the energy storage system was regarded as the objective function. The integer solution of minimum mass was obtained in combination with the graphical method. From the perspective of requirements and evaluation indexes of power and cost-effectiveness of vehicles, the influences of characteristic factors of the evaluation indexes, mass ratio of the HESS and conditions on vehicle performances were analyzed.

2. Theories and Analysis of Parameter Matching

The basic parameters of a Cens bus (the target model) are shown in Table 1. Herein, Cens ferric phosphate Li battery (CBP2450) was employed as the Li battery module unit. Its nominal voltage, specific power, specific energy, nominal capacitance and mass are 25.6 V, 1000 W/kg, 215.49 W∙h/kg, 50 A∙h and 19 kg, respectively. BMOD0165 (Maxwell, USA) was employed as a supercapacitor module unit. Its rated voltage, peak voltage, nominal capacitance, peak current and nominal mass are 48 V, 51 V, 165 F, 1900 A and 13.5 kg, respectively. The rated voltage, rated power and rated revolution speed of the motor are 300 V, 120 kW and 1860 r/min, respectively.

3. Theories and Analysis of Parameter Matching

3.1. Composite Energy Storage Device

In composite energy storage devices, two or more energy storage devices are combined so that the advantages of each energy source can be exploited, while its drawbacks can be compensated for by other energy storage devices. The combination of chemical battery and supercapacitor can overcome the low specific power of the battery and low specific energy of the supercapacitor while achieving high specific energy and high specific power [15,16]. It is assumed that the composite energy storage source consists of n (typically 2 or 3) basic energy storage devices:
{ P load = η 1 × P 1 + η 2 × P 2 + + η n × P n = η 1 × ρ P 1 × m 1 + η 2 × ρ P 2 × m 2 + + η n × ρ P n × m n W load = η 1 × W 1 + η 2 × W 2 + + η n × W n = η 1 × ρ W 1 × m 1 + η 2 × ρ W 2 × m 2 + + η 3 × ρ W n × m n
where n = { 1 , 2 , , n } ; P load and W load represent the power and energy needs of the load, respectively; P n and W n represent the power and energy provided by nth device, respectively; η n is the efficiency of the corresponding power converter; and ρ P n and ρ W n are power and energy densities, respectively.
In this study, the combination of a Li battery and a supercapacitor was considered:
{ P load = η 1 × P batt + η 2 × P capa = η 1 × ρ P _ batt × m batt + η 2 × ρ P _ capa × m capa W load = η 1 × W batt + η 2 × W capa = η 1 × ρ W _ batt × m batt + η 2 × ρ W _ capa × m capa

3.1.1. Technical Indicators of Typical Battery

For electric vehicles, the battery is the main energy source and converts chemical energy into electric energy for driving. In terms of the development of battery technology, current batteries include Pb–acid batteries, Ni–Cd batteries, Ni–Fe batteries, Ni–Zn batteries, Ni–metal hydride batteries, Al–air batteries, Fe–air batteries, Zn–air batteries and Li-ion batteries. Table 2 summarizes the specific energy, specific power, cycle life and price of typical batteries. Among them, the Li battery has been widely applied in electric vehicles by virtue of its excellent performance and environmentally friendly nature. However, it is limited by some drawbacks, including high price and poor low-temperature performance.

3.1.2. Parameters of Typical Supercapacitors

According to the U.S. Department of Energy, specific energy and specific power of supercapacitors used in electric and hybrid electric vehicles shall be no less than 5 W·h/kg and 500 W/kg, respectively, and the improved values should exceed 15 W·h/kg and 1600 W/kg, respectively. To date, no applicable supercapacitors that can fully meet these requirements have been reported. Nevertheless, high storage and transient large storage characteristics of supercapacitors play a dominant role in their applications, and the cost has dropped to 2.85 USD/kJ or 0.01 USD/F. The comparison of characteristics of different supercapacitors is shown in Table 3.

3.2. HESS Constraints

3.2.1. Technical Indexes of Typical Batteries

In addition to providing power and energy for the entire system, the composite energy storage device also consumes power and energy. Herein, the proposed effective energy utilization and power utilization constraints were analyzed.
The power required by the power source to drive the motor can be obtained through the Equations (A1)–(A3) in Appendix A (1). In addition to providing power for the power-source-driven motor, batteries also consume energy during driving:
P batt = ρ P _ batt m batt
η P _ batt = P batt P loss _ batt P batt
where ρ W _ capa and ρ W _ batt are energy densities of supercapacitor and battery, respectively; ρ P _ capa and ρ P _ batt are power densities of supercapacitor and battery, respectively; and η P _ capa and η P _ batt are effective power utilizations of supercapacitor and battery, respectively.
η = { 1 1 3600 ρ P _ batt η T ( f g cos α + g sin α ) v ( Consider   the   battery ) 1 1 3600 ρ P _ capa η T ( f g cos α + g sin α ) v ( Consider   the   supercapacitor )
With constant slope angle (α), power utilization is related to speed change and power density, and it is a positive number between 0 and 1. Therefore, the constraints of effective power utilization can be obtained by:
( f g cos α + g sin α ) v 3600 η T m capa ρ P _ capa ( m capa + m batt ) m batt ρ P _ batt ( m capa + m batt ) 0
Likewise, the constraints of effective energy utilization can be obtained by:
( f g cos α + g sin α ) S 3600 η T - m capa ρ W _ capa ( m capa + m batt ) m batt ρ W _ batt ( m capa + m batt ) 0

3.2.2. Constraint of Mileage

The most intuitive indicator of cost-effectiveness is single charge mileage or power consumption under cyclic conditions. Herein, the single charge mileage was employed. In this study, the Li battery in the composite power supply structure is directly connected to the brushless DC motor via the controller, and the number of battery packs in series is related to the level of the motor controller. U motor and U batt are defined as the voltage level of the motor controller and the nominal voltage of the battery unit, respectively. Hence, the number of battery packs in series can be determined by:
n ser _ batt = int ( U motor U batt )
where int indicates rounding up.
Assuming that m batt is the mass of a single battery, m batt is the total mass of battery and p is the number of battery packs in parallel, which is determined according to the design,
m batt = m batt int ( U motor U batt ) p
Analysis of bidirectional DC/DC efficiency shows that the working efficiency is high when the rated voltages of the high-voltage terminal and the low-voltage terminal are consistent. As a supercapacitor provides a large transient current for initiation, its terminal voltage should be slightly higher than that of the battery terminal. Assuming that m capa is the mass of single supercapacitor, m capa is the total mass of supercapacitor and q is the number of battery packs in parallel,
m capa = m capa int ( U motor U capa ) q
The total mass of the HESS is:
m HESS = m batt int ( U motor U batt ) p + m capa int ( U motor U capa ) q + m DCDC
When the vehicle travels at a constant speed ( v ys ), the power demand is smooth and the energy of the driving motor is provided by the battery. Hence, the constraints of mileage can be obtained by:
( ρ W _ batt m batt ) ( S / v ys ) P req
The power demand for driving at a constant speed can be expressed as:
P req = v ys 3600 η T ( m g f cos α max + C d A v ys 2 21.15 )
The constraints of mileage can be expressed as:
( ρ W _ batt m batt ) ( S / v ) v ys 3600 η T ( m g f + C d A v ys 2 21.15 )

3.2.3. Constraint of Maximum Power of Motor

When the maximum power of the motor is known, from the perspective of energy supply by energy storage source, it can be obtained that:
P max ( m batt ρ P _ batt + m capa ρ P _ capa ) 0
Herein, power and energy densities of the supercapacitor are 14,080 W/kg and 3.91 W·h/kg, respectively; power and energy densities of the battery are 1000 W/kg and 215.49 W·h/kg, respectively.

3.2.4. Constraint of Supercapacitor Capacitance

The energy provided the supercapacitor module can be expressed as:
W UC = 1 2 C ( U max 2 U min 2 )
where U max is the maximum discharge voltage of the supercapacitor, U min is the discharge cut-off voltage and W UCmax is the maximum energy released by the supercapacitor in road cycles and is related to the corresponding conditions. The constraint equation is:
1 2 C N ( U max 2 U min 2 ) W UC max

4. Linear Programming Parameter Optimization Method

4.1. Establishment of Optimized Objective Function

In this study, the masses of the battery ( m batt ) and supercapacitor ( m capa ) were selected as the optimization variables. Related to energy loss and cost during operation, quality is a key factor at the early stage of design. The overall objective function was the output. The HESS device includes Li battery, supercapacitor and DC/DC. The mass of DC/DC with fixed power is ignored in the optimization process as it is constant. Hence, the optimized objective function is as follows:
F ( m 1 , m 2 ) = m batt + m capa

4.2. Parameter Optimization of Linear Programming

Determination of the optimal solution using linear programming is essentially based on the assumption that only one feasible solution is available. m batt _ main = 2581.98   kg and m capa _ min = 628.65   kg . The number of units under the constraints in this study was considered. The number of battery modules was calculated to be 136. Based on the rated voltage of the driving motor and the integer principle, it was determined that 12 battery modules shall be connected in series in one branch. As energy is provided by the supercapacitor during starting acceleration, from the perspective of cost saving, 3.84 battery modules are absent with 11 branches in parallel selected, and an additional energy storage source is required to provide 74 kW power and 4.99 kW∙h energy. It was calculated that m batt = 2508   kg .
The selected supercapacitor module had a rated voltage of 48 V, discharge cut-off voltage of 24 V and maximum released energy of 1.43 × 10 5 J (39.6 W∙h). Likewise, seven modules in series are required for one supercapacitor. The mass of one supercapacitor mass was 13.5 × 7 = 94.5 kg and the mass of two supercapacitors was 13.5 × 7 × 2 = 189 kg; the energy provided by one supercapacitor was 277.2 W∙h and the energy provided by seven supercapacitors was 1940.4 W∙h (1.94 kW∙h). When all constraints are satisfied, m capa = 661.5   kg and the mass of DC/DC was assumed to be 100 kg. In this case, the total mass of the HESS was m HESS = 3269.5   kg .
The requirements of energy and power by the vehicle for the composite energy storage device can be described by the energy/power ratio ( R W / P ). R W / P is also an important parameter of the composite energy storage device. A large R W indicates that the HESS is a high-energy component, while a small R W indicates that the HESS is a high-power component.
R W / P = 3.91 × 661.5 + 67.37 × 2508 14,080 × 661.5 + 1000 × 2508 = 1.45 %
The graphic analysis method was employed to determine the constraints of the HESS. k is defined as the mass ratio of the supercapacitor and Li battery; k 1 is the mass ratio of the Li battery and the residual mass of the vehicle. Additionally, k p and k w are defined as power and energy coefficients, respectively; η P _ HESS and ρ P _ HESS are power utilization and density of the energy storage system, respectively. Constraints are as follows:
{ k W k k 1 m qt ρ W _ capa ( k k 1 m qt + k 1 m qt ) k 1 m qt ρ W _ batt ( k k 1 m qt + k 1 m qt ) 0 k P k k 1 m qt ρ P _ capa ( k k 1 m qt + k 1 m qt ) k 1 m qt ρ P _ batt ( k k 1 m qt + k 1 m qt ) 0 P max ( k 1 m qt ρ P _ batt + k k 1 m qt ρ P _ capa ) 0 ( ρ W _ batt k 1 m qt ) ( S / v ) v 3600 η T ( m g f + C d A v 2 21.15 ) 1 2 C N ( U max 2 U min 2 ) W UCmax
The constraints of power and energy utilization have a dominant effect on the size and shape of the reasonably selected area. If ρ P _ batt < k P < ρ P _ capa and ρ P _ batt - k P k P - ρ P _ capa < ρ W _ batt - k W k P - ρ W _ capa , Point E (intersection of the effective power utilization constraint curve and the energy feedback constraint curve) is higher than Point F (intersection of the energy feedback constraint curve and the dynamic performance constraint curve) but lower than Point A (intersection of the mileage constraint curve and the energy feedback constraint curve). The constraints obtained are shown in Figure 1. Herein, the horizontal and vertical coordinates are k 1 and k , respectively. As observed, Points B and C exhibited optimized mileage. Located at the boundary of dynamic performance and energy feedback constraints, Point B exhibited limited dynamic performance and insignificant energy feedback effect; its role as auxiliary energy storage source was also hindered. Point C was at the intersection of the effective power utilization and dynamic performance constraints, resulting in low effective power utilization of the energy storage source; a surplus was observed and the level of dynamic performance was not high. Similarly, Points A and D are at the intersection of the energy feedback constraint and mileage constraint and the intersection of the mileage constraint and effective power utilization constraint, respectively. Point A exhibited the lowest energy feedback and shortest mileage; Point D exhibited limited mileage and the lowest effective power utilization. Point E exhibited the lowest effective power utilization and energy feedback efficiency. If the mass ratio of the energy storage system corresponding to effective energy utilization constraints is higher than Point G (the intersection of the dynamic performance constraint and mileage constraint) or the mass ratio of the energy storage system corresponding to effective power utilization constraints is lower than Point F, the constraint shape developed would change.
Based on the parameters of the model vehicle selected, k and k 1 should not be over-large with all constraints satisfied. In this study, Point O (0.251, 0.264) was identified, which is consistent with the results of previous linear programming. In order to analyze the correlation of added auxiliary energy storage component ( m capa ) with m HESS and the relative relationship between the composite energy storage system and the mass of other components, the mass ratio k zlb and the mixing ratio k hhb were introduced. In this study, k zlb = 20.23 % and k hhb = 27.44 % .
An energy storage source provides energy consumed by itself and power and energy required for the entire vehicle during operation:
P batt η P _ batt + P capa η P _ capa > P other
ρ P _ batt m batt η P _ batt + ρ P _ capa m batt η P _ capa > k p m other + C d A v 3 76,140 η T
The quality relationships between each other were introduced:
k 1 > k P + C d A v 3 76,140 η T m ch η P _ batt k P
As observed, the mass of the Li battery was reduced by 73.98 kg. For the selected Li cell module unit, the mass of the Li battery was reduced by 228 kg from the perspective of the equivalent total terminal voltage of a single series battery pack. The capacitance was reduced by 50 A·h if one battery branch was removed. Due to the addition of the supercapacitor, R W / P decreased from 21.55% to 1.45% and the proportion of HESS as a functional component increased. In subsequent simulations, it was demonstrated that the component parameters of the selected driving system and energy storage system can reflect the overall performance of the vehicle; the accuracy and validity of the matching results were verified. The mass ratio of HESS and energy storage source and the reasonable range of the mass ratio of entire vehicle and components other than the energy storage source that satisfy the constraints were determined using the proposed parameter matching method. Additionally, the utilization can be improved by reasonable selection of the energy storage source unit if the initial matching conditions are satisfied.

5. Modeling and Experimental Verification of Electric Buses

5.1. Simulations of Consumption under Cycle Conditions

With power demand as the input signal, a dynamic model of the vehicle was established, as shown in Figure 2. In CRUISE, the internal electrical and mechanical parts of the composite energy storage system were connected.
In the full-load climbing experiment, the climbing performance was added to the Project-Task Folder during simulations. The maximum first-gear climbing slopes of the proposed electric buses and conventional electric buses with battery as single energy storage source were 29.4% and 18%, respectively, indicating a significant improvement in the climbing performance. It was significantly higher than the preset requirement on climbing performance (see Figure 3).
The mileage was calculated. Specifically, a vehicle traveling at a constant speed of 40 km/h was added in CRUISE and the corresponding driving distance when the SOC of the Li battery dropped 90% to 20% was defined as the mileage. Compared with conventional electric buses, the mileage of the proposed electric bus increased from 265 to 373 km, as shown in Figure 4.
As shown in Figure 5, the acceleration times of 0~50 km/h of conventional electric buses and those with HESS were 23.8 and 15.88 s, respectively. Therefore, it can be concluded that the speed changing rate of conventional electric buses is small and cannot satisfy the large transient energy demand.
As shown in Figure 6, the target speed and the actual speed almost coincide under the standard single road condition and the deviation was far less than 5%, demonstrating that the driving capacity of the motor and the energy provided by the energy storage source meet the performance requirements of buses owing to the parameter matching of the energy storage system.

5.2. Experimental Analysis

A test setup was established according to the similarity principle to investigate the performance of a matched HESS. In the test of the driving part, the power demand of the vehicle during operation was simulated by pair trawling of the brushless DC motor and the electric dynamometer. A brushless DC motor with a rated power of 5 kW and rated voltage of 48 V was employed as the motor to be tested. A DC dynamometer with a rated power of 7.5 kW was employed as the loading motor. Communication was established between controllers of the two motors through USB-to-CAN. Meanwhile, a vehicle demand model was developed in LabVIEW according to the actual vehicle speed demand in typical conditions to calculate the actual power demand. Additionally, the power was distributed between the Li battery and supercapacitor in the composite energy storage system according to the control algorithms. The experimental setup is shown in Figure 7.

5.2.1. Constant-Speed Climbing Performance

With motor speed of 970 rpm, the corresponding vehicle speed was calculated to be 8.9 km/h. For the dynamometer, the load resistance at the loading point increased from 2 N·m by 1 N·m/s Newton resistance. Figure 8 shows the measured waveform. As observed, the current increased and fluctuated at each turning point of slope increase. The current increased from 7.6 to 21.4 A as the slope increased. The simulated slope angle corresponding to the maximum torque was about 40%, which satisfies the requirements of dynamic performance.

5.2.2. Accelerating Capability

As shown in Figure 9, the vehicle speed increased from 8.3 to 17.5 km/h and the motor speed increased from 907.62 to 1907.33 rpm in simulations according to the set point. During acceleration, the vehicle delivers a driving force greater than the resistance; torque, speed and output current of the driving motor exhibit upward trends.

5.2.3. Regenerative Braking Feedback of Electric Vehicles

In this test, a standard urban deceleration road in China was employed for simulations. The motor speed dropped from 2940.5 (vehicle speed = 27.0 km/h) to 272.3 rpm (vehicle speed of 2.5 km/h). Forty sampling points were arranged and the relevant parameter information of one point was recorded every 1 s. Table 4 shows the experimental data obtained. The braking feedback test simulates the deceleration and braking of vehicles. During the test, the driving motor is decelerating or braking. Since two-way loading is possible for the dynamometer, the dynamometer drives the motor to simulate the transient process of deceleration upon braking. Besides decreasing speed, current and power of the motor fed back to the battery also decrease accordingly.

5.2.4. Continuous Road Condition Simulation

At the beginning of the test, the driving motor was started in the positive direction and the dynamometer provided a positive loading resistance to the motor. At the next sampling time, the dynamometer applied a reverse resistance torque to the motor (so that the driving operated). When the motor was exposed to reverse driving, it automatically entered the braking state; the dynamometer drove the motor until the speed was stabilized at 10 m/s. In the braking feedback test, three consecutive sampling points were set for reverse loading, and driving points were inserted between every two braking points. Figure 10 shows the trend of current on the DC side in the test. As observed, the output waveform of battery current measured by the current probe was positive during driving and negative during braking feedback.

6. Conclusions

With the minimum mass of the energy storage system as the objective function, a method combining linear programming and graphical method for parameter matching of a composite energy storage system is proposed. From the perspective of evaluation indexes and based on the rounding principle of unit quantity, the proportion of HESS as a power component increases to meet the energy demand for frequent starting and stopping of buses and buses; this method provides a novel solution to parameter matching of composite energy storage systems. Therefore, the proposed method is of great significance both theoretically and practically.

Author Contributions

Y.Z.: Conceptualization, Supervision, Writing—Review and Editing, Project Administration. J.L.: Writing—original draft, Investigation, Methodology, Validation. S.C.: Visualization. M.Z.: Software. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the Fundamental Research Foundation for Universities of Heilongjiang Province (2018-KYYWF-1672).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Simulation System Parameters

(1)
When the vehicle travels at a speed of v on a road with a slope of α, the power required by this vehicle at this moment is P req , the power consumed by the energy storage system is P HESS and the power consumed by other parts is P other :
P req = 1 3600 η T ( f m g cos α + C d A v 2 21.15 + m g sin α ) v
P HESS = 1 3600 η T ( f g cos α + g sin α ) ( m capa + m batt + m DCDC ) v
P other = 1 3600 η T ( f m other g cos α + C d A v 2 21.15 + m other g sin α ) v
where P other is the power demand of the entire vehicle except the energy storage device, P HESS is the power demand of the composite energy storage device, m is the vehicle mass, m HESS is the mass of the composite energy storage device, m DC / DC is the mass of DC/DC, m batt is the mass of the Li battery, m capa is the mass of the supercapacitor and m q t is the mass of other components.
(2)
The matched HESS energy/power ratio is defined as:
R W / P = E req P req = m batt ρ W _ batt + m capa ρ W _ capa m batt ρ P _ batt + m capa ρ P _ capa
In the formula, E req and P req are the total energy and total power of the vehicle during the traveling process.
k = m capa m batt ,   k 1 = m batt m qt
(3)
That is, m capa = k × m batt and m batt = k 1 × m qt ; therefore, m capa = k × k 1 × m qt .
k P = ( f × g × cos α + g × sin α ) × v 3600 η T ;   k W = ( f × g × cos α + g × sin α ) × S 3600 η T .
(4)
From this, η P _ HESS = 1 k P ρ P _ HESS .
According to the formula, the power coefficient is determined by the rolling friction coefficient, the slope angle and the vehicle speed. In the case of a certain slope angle, the higher the speed, the greater the corresponding power coefficient; the longer the distance traveled, the higher the corresponding energy coefficient.
k zlb = m capa m capa + m batt + m DCDC ;   k hhb = m capa + m batt + m DCDC m other .
(5)
The mass ratio k zlb is the percentage of the total mass of the supercapacitor in the HESS, and the mixing ratio k hhb is the ratio of the total mass of the HESS to the remaining mass of the vehicle.

References

  1. Masih-Tehrani, M.; Ha’Iri-Yazdi, M.R.; Esfahanian, V.; Safaei, A. Optimum sizing and optimum energy management of a hybrid energy storage system for lithium battery life improvement. J. Power Sources 2013, 244, 2–10. [Google Scholar] [CrossRef]
  2. Masih-Tehrani, M.; Dahmardeh, M. A Novel Power Distribution System Employing State of Available Power Estimation for a Hybrid Energy Storage System. IEEE Trans. Ind. Electron. 2018, 65, 6676–6685. [Google Scholar] [CrossRef]
  3. Li, M.; Wang, L.; Wang, Y.; Wang, Y.J.; Chen, Z.H. Sizing Optimization and Energy Management Strategy for Hybrid Energy Storage System Using Multiobjective Optimization and Random Forests. IEEE Trans. Power Electron. 2021, 36, 11421–11430. [Google Scholar] [CrossRef]
  4. Tostado-Veliz, M.; Arevalo, P.; Jurado, F. An optimization framework for planning wayside and on-board hybrid storage systems for tramway applications. J. Energy Storage 2021, 43, 103207. [Google Scholar] [CrossRef]
  5. Li, Q.; Chen, W.R.; Li, Y.K.; Liu, S.K.; Huang, J. Energy management strategy for fuel cell/battery/ultracapacitor hybrid vehicle based on fuzzy logic. Int. J. Electr. Power Energy Syst. 2012, 43, 514–525. [Google Scholar] [CrossRef]
  6. Li, Q.; Wang, T.X.; Li, S.H.; Chen, W.R.; Liu, H.; Breaz, E.; Gao, F. Online Extremum Seeking-based Optimized Energy Management Strategy for Hybrid Electric Tram Considering Fuel Cell Degradation. Appl. Energy 2021, 285, 116505. [Google Scholar] [CrossRef]
  7. Yu, Y.; Li, Q.; Huang, W.Q.; Chen, W.E. Operation Optimization and Control Method Based on Optimal Energy and Hydrogen Consumption for the Fuel Cell/Supercapacitor Hybrid Tram. IEEE Trans. Ind. Electron. 2021, 68, 1342–1352. [Google Scholar]
  8. Giannoutsos, S.V.; Manias, S.N. A Data-Driven Process Controller for Energy-Efficient Variable-Speed Pump Operation in the Central Cooling Water System of Marine Vessels. IEEE Trans. Ind. Electron. 2015, 62, 587–598. [Google Scholar] [CrossRef]
  9. Song, F.Z.; Liu, Y.; Xu, J.X.; Yang, X.F.; Zhu, Q. Data-Driven iterative feedforward tuning for a wafer stage: A high-order approach based on instrumental variables. IEEE Trans. Ind. Electron. 2019, 66, 3106–3116. [Google Scholar] [CrossRef]
  10. Montazeri-Gh, M.; Mahmoodi-K, M. Optimized predictive energy management of plug-in hybrid electric vehicle based on traffic condition. J. Clean. Prod. 2016, 139, 935–948. [Google Scholar] [CrossRef]
  11. Guo, L.X.; Zhang, X.D.; Zou, Y.; Guo, N.Y.; Li, J.W.; Du, G.D. Cost-optimal energy management strategy for plug-in hybrid electric vehicles with variable horizon speed prediction and adaptive state-of-charge reference. Energy 2021, 232, 120993. [Google Scholar] [CrossRef]
  12. Cao, J.; Emadi, A. A New Battery/UltraCapacitor Hybrid Energy Storage System for Electric, Hybrid, and Plug-In Hybrid Electric vehicles. IEEE Trans. Power Electron. 2012, 27, 122–132. [Google Scholar]
  13. Bayati, M.; Abedi, M.; Farahmandrad, M.; Gharehpetian, G.B.; Tehrani, K. Important Technical Considerations in Design of Battery Chargers of Electric Vehicles. Energies 2021, 14, 5878. [Google Scholar] [CrossRef]
  14. Shen, J.Y.; Dusmez, S.; Khaligh, A. Optimization of Sizing and Battery Cycle Life in Battery/Ultracapacitor Hybrid Energy Storage Systems for Electric Vehicle Applications. IEEE Trans. Ind. Inform. 2014, 10, 2112–2121. [Google Scholar] [CrossRef]
  15. Prasanthi, A.; Shareef, H.; Asna, M.; Ibrahim, A.A.; Errouissi, R. Optimization of hybrid energy systems and adaptive energy management for hybrid electric vehicles. Energy Conv. Manag. 2021, 243, 114357. [Google Scholar] [CrossRef]
  16. Hu, X.S.; Murgovski, N.; Johannesson, J.M.; Egardt, B. Comparison of three electrochemical energy buffers applied to a hybrid bus powertrain with simultaneous optimal sizing and energy management. IEEE Trans. Intell. Transp. Syst. 2014, 15, 1193–1205. [Google Scholar] [CrossRef]
  17. Li, Y.P.; Tang, X.L.; Lin, X.K.; Grzesiak, L.; Hu, X.S. The role and application of convex modeling and optimization in electrified vehicles. Renew. Sustain. Energ. Rev. 2022, 153, 111796. [Google Scholar] [CrossRef]
  18. Ravey, A.; Watrin, N.; Blunier, B.; Bouquain, D.; Miraoui, A. Energy-Source-Sizing Methodology for Hybrid Fuel Cell Vehicles Based on Statistical Description of Driving Cycles. IEEE Trans. Veh. Technol. 2011, 60, 4164–4174. [Google Scholar] [CrossRef]
  19. Kadyk, T.; Schenkendorf, R.; Hawner, S.; Yildiz, B.; Romer, U. Design of Fuel Cell Systems for Aviation: Representative Mission Profiles and Sensitivity Analyses. Front. Energy Res. 2019, 7, 35. [Google Scholar] [CrossRef] [Green Version]
  20. Gao, F.Y.; Zhang, H.R.; Wang, W.X. Energy Saving Operation Optimization of Hybrid Energy Storage System for Hydrogen Fuel Cell Tram. Trans. China Electrotech. Soc. 2020, 1, 293–299. [Google Scholar]
  21. Du, X.W.; Shen, Y.X.; Li, J. Energy management strategy of DC microgrid hybrid energy storage based on model predictive control. Power Syst. Prot. Control. 2020, 48, 69–75. [Google Scholar]
  22. Zhou, M.L.; Feng, J.F.; Zhang, Y. Composite energy storage system and its energy control strategy for electric vehicles. Electr. Mach. Control. 2019, 23, 51–59. [Google Scholar]
  23. Hu, J.J.; Xiao, J.; Yan, J.J. Control strategy and parameter optimization of hybrid energy storage device for electric vehicles. J. Chongqing Univ. 2016, 39, 1–11. [Google Scholar]
  24. Zhou, M.L.; Feng, J.F.; Zhang, Y. Research on power allocation control strategy for compound electric energy storage system of pure electric bus. Trans. China Electrotech. Soc. 2019, 34, 5001–5013. [Google Scholar]
  25. Lu, X.Y.; Wang, H.Y. Optimal Sizing and Energy Management for Cost-Effective PEV Hybrid Energy Storage Systems. IEEE Trans. Ind. Inform. 2020, 16, 3407–3416. [Google Scholar] [CrossRef]
  26. Wang, Y.X.; Su, J.; Hu, J.G. Research on on-board charging system for electric vehicle. J. Electr. Technol. 2019, 7, 62–66. [Google Scholar]
Figure 1. Conditional constraints of parameter matching.
Figure 1. Conditional constraints of parameter matching.
Machines 10 00085 g001
Figure 2. Pure electric bus model.
Figure 2. Pure electric bus model.
Machines 10 00085 g002
Figure 3. Comparison of full-load climbing performance.
Figure 3. Comparison of full-load climbing performance.
Machines 10 00085 g003
Figure 4. Comparison of driving mileage performance.
Figure 4. Comparison of driving mileage performance.
Machines 10 00085 g004
Figure 5. Comparison of 0~50 km/h acceleration curve.
Figure 5. Comparison of 0~50 km/h acceleration curve.
Machines 10 00085 g005
Figure 6. The comparison between the desired speed and current speed.
Figure 6. The comparison between the desired speed and current speed.
Machines 10 00085 g006
Figure 7. The experimental platform construction.
Figure 7. The experimental platform construction.
Machines 10 00085 g007
Figure 8. DC bus current when climbing with uniform speed.
Figure 8. DC bus current when climbing with uniform speed.
Machines 10 00085 g008
Figure 9. Speed, torque and current curve in accelerating test.
Figure 9. Speed, torque and current curve in accelerating test.
Machines 10 00085 g009
Figure 10. DC side current curve of continuous road test.
Figure 10. DC side current curve of continuous road test.
Machines 10 00085 g010
Table 1. Vehicle performance requirements and parameters.
Table 1. Vehicle performance requirements and parameters.
System ParametersValue
Overall vehicle performance requirementsMaximum speed (km/h)≥80
Maximum climbing degree (%)≥16
Acceleration time (s)≤25 (0~50 km/h)
Driving range (km)≥250 km (40 km/h uniform speed)
Overall vehicle parametersSize (mm)12,000 × 2530 × 3080
Unladen mass (kg)12,000
Full load mass (kg)17,800
Air resistance coefficient Cd0.65
Rolling resistance coefficient ff0.0165
Rolling mass conversion factor1.08
Windward area A (m2)6.385
Wheel radius (m)0.469
Table 2. Comparison of characteristics of different batteries.
Table 2. Comparison of characteristics of different batteries.
Battery TypeSpecific Energy
(W·h/kg)
Specific Power (W/kg)Cycle Life
(Times)
Price
(USD/kW·h)
Pb–acid20~10050~400500~2000120~150
Ni–Cd40~6080~350600~3000250~350
Ni–Fe50~6080~1501500~2000200~400
Ni–Zn55~75150~300600~1200100~300
Ni–metal hydride70~95200~300750~1200200~350
Al–air200~300160----
Fe–air80~12090500+50
Zi–air100~22030~80600+90~120
Li-ion100~130150~2501000+About 200
Table 3. Comparison of characteristics of different supercapacitors.
Table 3. Comparison of characteristics of different supercapacitors.
Brand NameRated Voltage (V)Rated Capacitance (F)Series Equivalent Resistance (mΩ)Energy
Density (W·h/kg)
Power Density (W/kg)Weight (kg)
Maxwell2.730000.295.5254000.55
Nesscap2.750000.255.4450630.93
Panasonic2.525000.433.7010350.395
EPCOS2.734000.454.37600.60
ESMA1.780,0000.513.385832.4
Spscap2.795000.287.4050101.3
BMOD011542145102.22290016
BMOD01171443541.8219006.5
BCAP00102.526000.74.343000.525
Company A1.6140,0000.312.02055.0
Company B2.720499.8810000.116
Company C2.7100154.916770.021
Company D2.730000.43.8510200.41
Company E1.215002.43.338000.36
Table 4. Data of regenerative braking test.
Table 4. Data of regenerative braking test.
Output Voltage (V)Output Current (A)Input Torque (N·m)Input Power (W)Output Power (W)Efficiency (%)
53.39−26.994−6.6341860.56−1441.4277.47
53.62−26.856−6.6131866.72−1440.2377.15
53.77−26.235−6.5871857.65−1410.8575.95
53.87−26.172−6.5761838.85−1409.8876.67
53.93−25.784−6.5651823.40−1390.7176.27
54.99−25.744−6.5591811.09−1390.0876.75
54.03−25.496−6.5491800.30−1377.7776.53
54.08−24.907−6.5271746.41−1347.1777.14
54.12−24.893−6.5241748.61−1347.2277.05
54.12−24.404−6.5161736.84−1320.7676.04
54.14−24.380−6.5121720.43−1320.0476.73
54.16−24.187−6.5031709.27−1309.9776.64
54.18−24.194−6.5051704.14−1310.9476.93
54.15−23.273−6.4841654.57−1260.3676.17
54.17−23.284−6.4821644.83−1261.3776.69
54.10−21.677−6.4501541.27−1172.7376.09
54.11−21.682−6.4491534.07−1173.2176.48
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhang, Y.; Liu, J.; Cui, S.; Zhou, M. Parameter Matching Methods for Li Battery–Supercapacitor Hybrid Energy Storage Systems in Electric Buses. Machines 2022, 10, 85. https://doi.org/10.3390/machines10020085

AMA Style

Zhang Y, Liu J, Cui S, Zhou M. Parameter Matching Methods for Li Battery–Supercapacitor Hybrid Energy Storage Systems in Electric Buses. Machines. 2022; 10(2):85. https://doi.org/10.3390/machines10020085

Chicago/Turabian Style

Zhang, Yu, Jiahong Liu, Shumei Cui, and Meilan Zhou. 2022. "Parameter Matching Methods for Li Battery–Supercapacitor Hybrid Energy Storage Systems in Electric Buses" Machines 10, no. 2: 85. https://doi.org/10.3390/machines10020085

APA Style

Zhang, Y., Liu, J., Cui, S., & Zhou, M. (2022). Parameter Matching Methods for Li Battery–Supercapacitor Hybrid Energy Storage Systems in Electric Buses. Machines, 10(2), 85. https://doi.org/10.3390/machines10020085

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop