This paper establishes the simulation model of fuel cell hybrid golf cart (FCHGC), which applies the non-GUI mode of the Advanced Vehicle Simulator (ADVISOR) and the genetic algorithm (GA) to optimize it. Simulation of the objective function is composed of fuel consumption and vehicle dynamic performance; the variables are the fuel cell stack power sizes and the battery numbers. By means of simulation, the optimal parameters of vehicle power unit, fuel cell stack, and battery pack are worked out. On this basis, GUI mode of ADVISOR is used to select the rated power of vehicle motor. In line with simulation parameters, an electrical golf cart is refitted by adding a 2 kW hydrogen air proton exchange membrane fuel cell (PEMFC) stack system and test the FCHGC. The result shows that the simulation data is effective but it needs improving compared with that of the real cart test.
1. Introduction
The problem with fossil fuel combustion processes of atmospheric pollution and greenhouse effect has plagued the earth. In order to protect the environment and avoid energy crises, fuel cells are employed as power to drive vehicles. PEMFC due to the high power density, low temperature operating and quick start up, is very suitable for the vehicle power source, whose raw material hydrogen is ubiquitous and its reactant water is environmentally friendly; for this reason, it is an optimal resource [1–3]. However, the production and using costs of fuel cells are extremely expensive, though large fuel cell stacks are used in some buses at home and abroad, most of which are used in testing [4, 5]. It is of high requirements in both mass and heat transfer of a large stack, resulting in the short using lifespan of the stack. Mini stacks are selected instead of large ones to test the FCHGC so as to reduce costs. They can also be applied in some mini electrical vehicles to reduce the pollution caused by lead-acid battery.
There are some control strategies in fuel cell hybrid vehicle (FCHV), such as fuzzy control and state of charge (SOC) control [6–8], which refer to energy management and SOC balance. This paper employs control strategy based on power following, applying vehicle simulation software ADVISOR and GA to conduct vehicle simulation, comparing its results with the real vehicle testing to analyze the mini FCHV; it aims to find out the appropriate optimization parameters.
The layout of this paper is as follows. Section 2 demonstrates the use of ADVISOR and GA, as well as the vehicle simulation mode establishment. Section 3 analyses the vehicle structure mode, mainly elaborating the stack size, the batteries mode, and the vehicle out-looking structure that the vehicle matches. Section 4 conducts road test of the vehicle. Section 5 analyses the result of the simulation and the actual testing. Section 6 draws a conclusion of the whole paper.
2. The Establishment of the Simulation Mode2.1. ADVISOR Simulation Model
Vehicle simulation tools are mainly commercial software, such as ADVISOR and CRUISE, as well as simulation tools developed by researchers [9–12]. ADVISOR can match components mode selection of the vehicles beforehand; it can also estimate the energy control strategy by running the vehicles, optimizing all sorts of using parameters [13, 14]; meanwhile, this software is applied by vehicle manufacture enterprises and vehicle research institutions. ADVISOR employs both the forward and the backward hybrid simulation method, mainly the backward one and the forward simulation as its auxiliaries, in accordance with the energy demands to assign vehicle power energy.
The internal module of ADVISOR involving parameters is primarily from large business vehicles, while the mini ones need to be developed again, changing some parameters and structure so as to launch the FCHGC simulation. M-file of the main components including FUEL_CELL_defaults_in.m, block_diagram_name.m, VEH_SMART.m, pb.m in the toolkits of ADVISOR, and the driving cycle of the road testing are altered; the block diagram is demonstrated in Figure 1 [15].
ADVISOR block diagram.
2.2. The Objective Functions and Constraint
One of the vehicle objectives is to reduce the fuel consumption; namely, it is to reduce hydrogen consuming as far as the FCHV is concerned. FCHGC energy consumption is to provide power for electric motor, whose power is mainly composed of DC/DC converter (fuel cell stack) and batteries. In order to allocate power of batteries and DC/DC converter rationally, the control strategy necessitates the fuel cell stack work in the high efficient area and maintains the battery’s SOC in a reasonable range so as to prevent the batteries from discharging too much power or charging too little and to protect the stack and the batteries, which are introduced in [16].
The calculating function of the fuel cell consuming
(1)minQH2=∫0tf(PDC/DC,Pbat,Pacc)dt,
where QH2 is the consuming of hydrogen in the vehicle operating process, PDC/DC is the output power of the converters, Pbat is the output power of battery, and Pacc is accessory power.
Equation (2) is the power calculation of the batteries and the converter [6, 17]:
(2)f(PDC/DC,Pbat)=f(PDC/DC)+f(Pbat)+f(Pacc),f(PDC/DC)=Pfc*ηDC/DC,f(Pbat)={Pbat.disηdisηchg¯,Pbat≥0,ηdis¯ηchgPbat.chg,Pbat<0,f(Pacc)=K,ηdis=12+121-4RdisPbat.disUoc,ηchg=(12+121-4RchgPbat.chgUoc)-1,
where Pfc represents stack power; Pbat.dis and Pbat.chg are transient discharging and charging power of batteries; ηchg, ηchg¯ are the charging efficiency and average charging efficiency of batteries; ηdis, ηdis¯ are the discharging efficiency and average discharging efficiency of batteries; ηDC/DC is converter efficiency; K is set to be a constant accessories consuming power; Uoc is battery open-circuit voltage; and Rchg, Rdis are batteries charging internal resistance and discharging internal resistance.
Equations (3) are constraining ones:
(3)PDC/DCmin≤PDC/DC≤PDC/DCmaxPbatmin≤Pbat≤PbatmaxSOC_H≤SOC≤SOC_LPreq=PDC/DC+Pbat+Pacc0≤Vmax≤32km/h,
where PDC/DCmin and PDC/DCmax are the minimum and maximum DC/DC power output; Pbatmin and Pbatmax are the minimum and maximum batteries power output; SOC_H and SOC_L are the maximum and minimum SOC; Preq represents the power requirements of the motor; and Vmax is the maximum speed of vehicle.
While calculating simulation, DC/DC converter efficiency needs to be added to the ADVISOR mode which has no DC/DC module.
In order to better reflect the features of vehicle dynamic and fuel cell economy, the objective function is revised to be in line with the weighting of the vehicle maximum speed and fuel consuming. The maximum vehicle speed is in accordance with the vehicle designing demands; the estimated maximum vehicle speed is 32 km per hour.
In accordance with GA to seek the minimum, the maximum vehicle speed is revised to be
(4)Vmodify=32-accel.max_speed,
where accel.max_speed is the maximum speed of vehicle accelerating test.
The objective function of the joint simulation is
(5)Objf=λ1·QH2+λ2·Vmodify,
where λ1, λ2 is weighting, and λ1>0, λ2>0, λ1+λ2=1.
2.3. The Combination of GA and ADVISOR
Based on the concept of chromosome in the biology and genes, GA simulates the mechanism of genes and evolution of the natural biology to optimize, it is especially appropriate for the multiobjective, constrained nonlinear and multimode optimization simulation. Currently, GA is largely applied in the vehicle engineering to conduct optimization simulation, especially the vehicle calculation with features of multiobjective, constrained nonlinear, which are introduced in references [18, 19]. The basic steps in GA are mentioned in references [19].
The GUI mode of ADVISOR [20], whose mode is visual, easy to operate, and convenient to interact, is usually used, whereas, it is difficult for the GUI mode to carry through combination simulation with other algorithms, whose data cannot be exchanged; in this case, non GUI mode of ADVISOR is employed. A command function [error_code, resp] = adv_no_gui (action, input) is applied in using non-GUI mode.
In accordance with the calculating characteristics of the vehicle model selection and evaluation, Pfc_max and numbat are selected as variables, by revising the commands of non-GUI to simulate, mainly revising the command line [15]:
The essential steps of combination simulation are as follows.
Initialize ADVISOR and GA workspace;
Calculate initial population and lead in ADVISOR to work out objective function value. To judge whether terminate the condition or not, quit if feasible, or enter the next step;
The procedure of GA: the objective function value is ranking, operators (select, cross connection or copy, mutation) are selected, obtaining the new population after calculating;
ADVISOR is used to calculate objective function. To judge whether terminate the condition or not, quit if feasible, or enter the next step. The overall procedure is shown in Figure 2.
Flowchart for combination simulation.
The GUI functions of ADVISOR are fully used besides the joint simulation. It mainly utilizes the GUI functions to select the motor rated power after making sure the maximum power of the optimal stack and the battery numbers. The ultimate purpose is to make the stack, the batteries, and the motor work coordinately in order to guarantee the dynamic performance and the economy of the vehicle.
3. The Structural Models of an FCHGC
The FCHGC is refitted from the golf cart electric vehicle, retaining their fundamental structure, adapting the dynamic components, and adding 2 kW fuel cell stack and air supplying system; the passengers are cut down from eight to five, the vehicle parameters and system configuration; see Table 1 and Figure 3.
Parameters of the vehicle structure and power [21].
Parameters
Data
Vehicle structure
Size (m)
3.75 (L) × 1.35 (W) × 1.88 (H)
Wheelbase (m)
2.58
Minimum turning radius (m)
4.5
Maximum speed (km/h)
30
Motor voltage (V)
48
Stack parameters
Number of cells
110
Rated voltage (V)
66
Rated power (kW)
2 (0.6 V/cell)
Rated current (A)
30
Lead-acid battery pack
Number of batteries
8
Nominal voltage of a single battery (V)
6
Electric capacity of a single battery (A h)
210 (5 hr)
Maximum charging voltage (V)
53
Electrical schematic of the vehicle dynamic control.
The FCHGC uses 2 kW half integration PEMFC stack as main energy provider, mainly providing operating power to the vehicle and charging the batteries once running in the low load road according to demands. For the detailed parameters see Table 1.
The 3-D-210 specialized in electrical vehicle is applied as the power batteries, whose series are composed of eight batteries as auxiliary power system, supplying auxiliary power as soon as the vehicle starts, the stack activated, or when it has insufficient power.
4. FCHGC Road Test
Considering the maximum speed of the vehicle does not exceed 32 km per hour, the maximum climbing gradient does not exceed 11% and the vehicle cannot do high strength test, such as CYC_UDDS, CYC_ECE, and CYC_HWFET, also taking into account the FCHGC runs on a flat road. In this case, the vehicle testing group selected a flat road to accelerate and decelerate the cycle test and the relevant climbing test. The test was carried out in the South China University of Technology, North Campus, where the road is concrete pavement and flat. Furthermore, before the real road test, the FCHGC is carried out on the test bench so as to subtract test costs. Some small items are divided during testing, say, discontinuously accelerating, decelerating, and climbing as well as fixing load [22]. The testing is shown in Figure 4.
Vehicle test.
The test purposes are to obtain information of the stack dynamic performances, such as the maximum power, the corresponding time, and the durability of the stack, as well as the matching condition of the batteries and the stack. It reflects the dynamic performance, the speeding performance, and the maximum speed.
5. Analysis of Results
The FCHGC used ADVISOR for simulation, which has no proper driving cycle; therefore, the corresponding revision is made based on CYC_UDDS to form a new driving cycle diagram CYC_UDDS_golf, as shown in Figure 5. The subsequent simulation is carried out based on the new driving cycle.
CYC_UDDS_golf drive cycle.
The weighting λ1=0.7 and λ2=0.3 are selected in the simulation to iterate; the result is shown in Figure 6. After 15 generations of iteration, the optimal target and the average target vary to achieve stability, whose fuel consumption is 76.7 (L/100 km). By this time, variable value Pfc_max is 2 kW and numbat is 4.7.
Convergence of combination simulation.
Table 2 is the simulation results; line one demonstrates different parameter symbols, line two to six indicate the simulation results of different weighting, and the last line represents simulation results that simulate through GUI mode after appropriate dynamic component is selected. Table 2 shows the impact on the simulation results caused by different weighting of the objective function. With the increase of fuel economy weight λ1, the number of batteries increases, whilst the stack power remains 2 kW; meanwhile, fuel consumption is reduced. This testing forms better economy.
The data of different weighting after the combination simulation.
Weighting
numbat
Fuel Consumption (L/100 km)*
Distance (km)
λ1=0.8, λ2=0.2
5.2696
55.4
2.1
λ1=0.7, λ2=0.3
4.7304
76.7
1.6
λ1=0.6, λ2=0.4
4.7274
76.5
1.6
λ1=0.5, λ2=0.5
4.7304
76.7
1.6
λ1=0.4, λ2=0.6
4.1657
88.7
1.4
λ1=0.3, λ2=0.7
4.1667
88.7
1.4
After combination simulation
8
46.4
2.3
*ρH2=18 (g/L), density of the fuel; lhv = 120000 (J/g), lower heating value of the fuel.
λ1 represents the weighting of fuel consumption, λ2 represents the weighting of dynamic performance.
Figures 7 to 10 are the comparison between the speed of driving cycle and that of simulation; cyc_mph_r is the speed curve of driving cycle; mpha is the speed curve of simulation. Figures 7 and 8 demonstrate the speed curve of simulation, which cannot reach the speed demands of the driving cycle; the two curves are obviously different, indicating the insufficiency of vehicle power. The two curves in Figure 9 are slightly different, while the two curves in Figure 10 are completely coincident.
λ1=0.4, λ2=0.6, the speed of simulation & drive cycle.
λ1=0.6, λ2=0.4, the speed of simulation & drive cycle.
λ1=0.8, λ2=0.2, the speed of simulation & drive cycle.
The speed of simulation and drive cycle.
In accordance with the simulation, a stack of 2 kW is selected, while the battery numbers are required to be rounded off, six is selected as numbat. But according to Figure 9, there is a little difference between the two curves; it shows the insufficiency of vehicle power, battery numbers need increasing, for this reason, the battery numbers of the real vehicle are eight. By selecting 2 kW as Pfc_max, eight as numbat, simulating GUI mode of ADVISOR, whose results are shown in Figure 10 and Table 2, the data indicates the minimum fuel consuming of the power matching, satisfying the maximum dynamic demands, achieving the prospective effect.
The above joint simulation estimates the motor rated power to be 3 kW, whether it meets the requirements or not, the following simulation is simulated by means of GUI mode of ADVISOR, whose results are shown from Figure 11 to Figure 13. Figure 11 indicates that in the overall driving cycle, the vehicle driving distance is less than 2.3 km with a motor under a 2 kW rated power; moreover, the minimum fuel consumption of the motor power is 34.6 L/100 km with a motor power of 1.5 kW, and then it increases gradually. Figures 12 and 13 demonstrate that there is a difference between the simulation speed curve and the input driving cycle to the motor rated power less than 2 kW. It indicates that the selection of the undersize of the motor rated power consumes little, but the dynamic performance of vehicle is bad; however, selection with oversize increases the fuel consumption, and in this case, a 3 kW is optimal to be selected as the rate power of the motor.
Fuel consumption and vehicle versus motor power.
Pmotor=1kW, the speed of simulation and drive cycle.
Pmotor=2kW, the speed of simulation and drive cycle.
Road tests of climbing and flat ground are to verify the using status of the stack power and the batteries, comparing with the simulation results. Figure 11 indicates the hybrid test of the FCHGC, which discontinuously accelerates and decelerates on the road. When motor needs large power, the stack power of real testing is working fluctuation in 2 kW attachments, however, the maximum stack power of the simulation does not exceed 2 kW. The simulation software also makes the stack supply the accessories with 300 W power, while the real vehicle is supplied with batteries.
With a gradient of 11%, the vehicle climbing status is described in Figure 12, which shows the stack power ascends and fluctuates in the road test, not reaching the maximum power 2 kW, while the stack power of the simulation immediately reaches 2 kW and remains in this stage until the climbing ends. The cause is the vehicle motor power demand increases; the stack is required to provide more power, but it cannot respond immediately, because its reactant reaction and by-products and heat transfer need time; these factors cause the stack power oscillation. Simultaneously, in order to protect the stack, DC/DC do not allow too much change of the stack power; therefore, the batteries need to provide more power [23–25]. However, ADVISOR involving mass and heat transfer is static data; therefore, the stack provides the maximum power.
It reflects the discrepancy between the simulation and the real test; simulation accuracy needs improving in practice so as to better serve the vehicle development.
6. Conclusion
By analyzing the structure features and the using condition of the FCHGC and the problems emerged in the process of ADVISOR simulation, some viewpoints are concluded as follows.
The joint simulation of GA and ADVISOR can fast and accurately solve the problems of the vehicle power selection, mitigating the product research and development, saving the development cost.
By analyzing the optimization of the FCHGC, the appropriate fuel consumption and dynamic performance of mini PEMFC are found out, which can be widely applied and be used as a research reference for fuel cell vehicles with high-efficiency.
Simulation is one of the necessary work of the vehicle development, while simulation cannot replace real testing because the simulation software has some static data that influences the simulation result; in consequence, there still exists a large improvement.
NomenclatureP:
Power
R:
Resistance
U:
Voltage
K:
Constant accessories consuming
SOC_H:
Maximum SOC
SOC_L:
Minimum SOC
V:
Speed
num:
Battery number
GA:
Genetic algorithm
FCHGC:
Fuel cell hybrid golf cart
ADVISOR:
Advanced Vehicle Simulator
FCHV:
Fuel cell hybrid vehicle
SOC:
State of charge
PEMFC:
Proton exchange membrane fuel cell
Q:
The consuming of hydrogen.
Subscriptsλ:
Weighting
η:
Efficiency
bat:
Battery
acc:
Accessory
DC/DC:
Direct current converter
fc:
Fuel cell
dis:
Discharging
chg:
Charging
oc:
Open-circuit
min:
Minimum
max:
Maximum
req:
Requirement.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgment
The research was supported by Guangdong Science and Technology Department, China (Funding: no. 2011A090200075, no. 2012B091100330).
TirnovanR.GiurgeaS.Efficiency improvement of a PEMFC power source by optimization of the air management2012379774577562-s2.0-8486029297710.1016/j.ijhydene.2012.02.029CorboP.MigliardiniF.VeneriO.Performance investigation of 2.4 kW PEM fuel cell stack in vehicles20073217434043492-s2.0-3654907551910.1016/j.ijhydene.2007.05.043AkyalçinL.KaytakoǧluS.Optimization of structural combinations on the performance of a PEMFC's MEA200818027677722-s2.0-4274909537510.1016/j.jpowsour.2008.02.036AhluwaliaR. K.WangX.KwonJ.RousseauA.KalinoskiJ.JamesB.MarcinkoskiJ.Performance and cost of automotive fuel cell systems with ultra-low platinum loadings201119610461946302-s2.0-7995228205510.1016/j.jpowsour.2011.01.059BaptistaP.TomásM.SilvaC.Plug-in hybrid fuel cell vehicles market penetration scenarios2010351810024100302-s2.0-7795617660110.1016/j.ijhydene.2010.01.086BarelliL.BidiniG.OttavianoA.Optimization of a PEMFC/battery pack power system for a bus application2012977777842-s2.0-8265517198410.1016/j.apenergy.2011.11.043KisacikogluM. C.UzunogluM.AlamM. S.Load sharing using fuzzy logic control in a fuel cell/ultracapacitor hybrid vehicle2009343149715072-s2.0-5864911237710.1016/j.ijhydene.2008.11.035CauxS.HankacheW.FadelM.HisselD.On-line fuzzy energy management for hybrid fuel cell systems2010355213421432-s2.0-7554908903110.1016/j.ijhydene.2009.11.108BrownD.AlexanderM.BrunnerD.AdvaniS. G.PrasadA. K.Drive-train simulator for a fuel cell hybrid vehicle200818312752812-s2.0-4674914915710.1016/j.jpowsour.2008.04.089MusioF.TacchiF.OmatiL.Gallo StampinoP.DotelliG.LimontaS.BrivioD.GrassiniP.PEMFC system simulation in MATLAB-Simulink environment20113613804580522-s2.0-7995761710510.1016/j.ijhydene.2011.01.093HosseinzadehE.RokniM.AdvaniS. G.PrasadA. K.Performance simulation and analysis of a fuel cell/battery hybrid forklift truck2013381142414249/10.1016/j.ijhydene.2013.01.168WuZ.ZhangJ.JiangL.WuH.YinC.The energy efficiency evaluation of hybrid energy storage system based on ultra-capacitor and LiFePO4 battery201211395105XuJ.-W.ZhengL.Simulation and analysis of series hybrid electric vehicle (SHEV) based on ADVISORProceedings of the International Conference on Measuring Technology and Mechatronics Automation (ICMTMA '10)March 20103543572-s2.0-7795303921910.1109/ICMTMA.2010.678SameA.StipeA.GrossmanD.ParkJ. W.A study on optimization of hybrid drive train using advanced vehicle simulator (ADVISOR)201019519695469632-s2.0-7795312771910.1016/j.jpowsour.2010.03.057BrookerA.HaraldssonK.HendricksT.JohnsonV.KellyK.KramerB.ADVISOR Documentation (Version 2002)National Renewable Energy Lab., April 2002RaveyA.BlunierB.MiraouiA.Control strategies for fuel-cell-based hybrid electric vehicles: from offline to online and experimental results20126162452245710.1109/TVT.2012.2198680XuL.LiX.HuaJ.LiJ.OuyangM.Parameter identification and control strategy optimization of hybrid fuel cell powertrain200945256612-s2.0-6164910289110.3901/JME.2009.02.056FangL.QinS.XuG.LiT.ZhuK.Simultaneous optimization for hybrid electric vehicle parameters based on multi-objective genetic algorithms2011435325442-s2.0-7995321983110.3390/en4030532ArsalisA.NielsenM. P.KærS. K.Modeling and optimization of a 1 kWe HT-PEMFC-based micro-CHP residential system2012373247024812-s2.0-8485584507710.1016/j.ijhydene.2011.10.081LiQ.ChenW.LiY.LiuS.HuangJ.Energy management strategy for fuel cell/battery/ultracapacitor hybrid vehicle based on fuzzy logic201243514525TangY.YuanW.PanM.WanZ.Experimental investigation on the dynamic performance of a hybrid PEM fuel cell/battery system for lightweight electric vehicle application201188168762-s2.0-7795730151110.1016/j.apenergy.2010.07.033LiangJ.JianQ.Decay analysis of a 5kW PEMFC stack in sightseeing busProceedings of the 3rd International Conference on Frontiers of Manufacturing Science and Measuring Technology201318261829MatsuuraT.SiegelJ. B.ChenJ.StefanopoulouA. G.Multiple degradation phenomena in polymer electrolyte fuel cell operation with dead-ended anodeProceedings of the 9th Fuel Cell Science, Engineering and Technology ConferenceAugust 2013Washington, DC, USA19HwangJ. J.ChangW. R.WengF. B.SuA.ChenC. K.Development of a small vehicular PEM fuel cell system20083314380138072-s2.0-4894911560510.1016/j.ijhydene.2008.04.043CorboP.MigliardiniF.VeneriO.Experimental analysis of a 20 kWe PEM fuel cell system in dynamic conditions representative of automotive applications20084910268826972-s2.0-4884908651410.1016/j.enconman.2008.04.001