Mapping Fuzzy Energy Management Strategy for PEM Fuel Cell–Battery–Supercapacitor Hybrid Excavator

By replacing conventional supplies such as fossil fuels or internal combustion engines (ICEs), this paper presents a new configuration of hybrid power sources (HPS) based on the integration of a proton-exchange membrane fuel cell (PEMFC) with batteries (BATs) and supercapacitors (SCs) for hydraulic excavators (HEs). In contrast to conventional architectures, the PEMFC in this study functions as the main power supply, whereas the integrated BAT–SC is considered as an auxiliary buffer. Regarding shortcomings existing in the previous approaches, an innovative energy management strategy (EMS) was designed using a new mapping fuzzy logic control (MFLC) for appropriate power distribution. Comparisons between the proposed strategy with available approaches are conducted to satisfy several driving cycles with different load demands and verify the strategy’s effectiveness. Based on the simulation results, the efficiency of the PEMFC when using the MFLS algorithm increased up to 47% in comparison with the conventional proposed EMS and other approaches. With the proposed strategy, the HPS can be guaranteed to not only sufficiently support power to the system even when the endurance process or high peak power is required, but also extend the lifespan of the devices and achieves high efficiency.


Introduction
Hydraulic excavators (HEs) are essential to construction and agriculture [1][2][3][4][5][6][7][8][9][10][11], however, existing issues such as low energy efficiency and carbon emissions in HEs using an internal combustion engine (ICE) have urgently imposed alternative energy requirements with the aim of reducing environmental and economic problems. In the trend of exploring and developing alternative power sources to take the place of fossil fuels, the proton-exchange membrane fuel cell (PEMFC) has become a potential candidate due to its high performance in high energy conversion efficiency, low chemical pollution, quiet operation, low weight, and low volume, especially in producing zero-emissions. However, low power density and slow power response are drawbacks of the PEMFC. Moreover, transient change in the load demand due to complicated requirements or working terrains is difficult to adequately satisfy using the standalone PEMFC. Furthermore, when the HEs are in regeneration modes or when the PEMFC releases more power than the demand, the PEMFC itself is not able to store the excess energy. A power deficiency when increasing the operating load and the waste power in the case of regeneration mode require at least one auxiliary device to enhance and optimize performance. employed the FLC to manage power-sharing between three different sources and guarantee their performance. The simulations were implemented in four standard driving cycles and the results proved that the suggested algorithm could satisfy the power demanded in comparison with other methods. Using the advantages of FLC, Issam et al. applied this technique for the EMS used in stand-alone applications [27]. Another application was introduced by Fei et al. when a master-slave EMS-based fuzzy logic hysteresis state machine was exploited for power split management applied to a PEMFC-lithium battery-supercapacitor (PEMFC-LIB-SC) hybrid tramway powertrain system [28]. Saman et al. [29] utilized fuzzy rules combined with a genetic algorithm (GA) to improve the vehicle's dynamic performance of the hybrid FC-BAT-SC system and fuel economy. As the author claimed, the study contributed to related research from different aspects. Summarizing from previous approaches, Yakup et al. realized an analysis and comparison of several control methodologies for fuel-saving in hydrogen fuel cell vehicles integrated with BAT and SC [30]. In [31], Zhang et al. conducted an investigation on unplugged electric vehicles (UEVs) by employing an online fuzzy for the power management of the HPS composed of the PEMFC and BAT, and presented the procedure for power split management. The power required from the FC was produced after being processed in the FLC. The pulsed-power profile was used to verify its effectiveness and stability to the rapid change of the instantaneous high-power demand. With a similar rule-based approach used in previous works, Cong et al. established a new on-off power following a strategy-based fuzzy algorithm for FC extended-range vehicles in which the BAT was employed as the main source while the FC was considered as an extender for the system [32]. The simulation demonstrated that good results and improvement could be accomplished in both system dynamics and economy. Eckert et al. [33] utilized GA to optimize the lower and upper values of the membership function to maximize system performance. In [34], Ameur et al. proposed a master-slave fuzzy algorithm to design the EMS to effectively manage a hybrid power system and prolong all component lifespans.
Based on the literature, the FLC that mimics human knowledge [35][36][37] is an appropriate tool for establishing the EMS and applying it to the complicated hybrid system. Despite achieving good performance, few studies have reported on an application of these accomplishments for HEs. Furthermore, due to mimicking designer knowledge, it is challenging to achieve a good performance without setting an appropriate value and condition for the fuzzy system. Optimization methods have been successfully applied for the experimental system; nevertheless, most configurations have been conducted in two-device configurations such as FC-BAT, FC-SC, and BAT-SC. When configurations of three devices are considered, the power-sharing strategies are more complicated. Although the FLC, known to mimic designer knowledge, is difficult to design, it is still a suitable solution for a power-sharing design of a hybrid FC-BAT-SC. This technique can cover all scenarios of the working operation and mutual influences of charging-discharging mode in BAT-SC to maintain the high performance of all components. Considering the advantages of the integration of FC-BAT-SC and a lack of algorithms when designing the EMS for HEs, this study proposed a novel algorithm in which the FLC is exploited as the EMS to appropriately distribute power to each component. Building on our previous work in [38], we focused on how to operate all components in an optimal condition as the prior criterion to increase efficiency, thus prolonging lifespan and reducing costs. In order to handle these problems, an optimal mapping FLC (MFLC) was systematically designed based on a self-tuning methodology to generate the optimized parameters of all devices, thus adjusting power-sharing. The contributions of this research are as follows: (1) With our experience in fluid power and construction HEs, we proposed a novel configuration for the integration of FC-BAT-SC and attempted to apply it to HEs. In contrast to the conventional configuration in which the BATs are considered as the main supply, in this configuration, the FC functions as the main power source, and the BAT-SC is attached as supplements. (2) A novel EMS strategy was introduced in which the MFLC was designed to match the optimal condition during operation. While FLC was employed to distribute sufficient power to each component under different scenarios, the mapping condition was first introduced to calculate Energies 2020, 13, 3387 4 of 27 a suitable fuel cell power. This control scheme is the key point to addressing problems associated with HE power distribution, which are considered as constrained multi-objective problems. The effectiveness of the proposed algorithm was validated by the standard driving cycle in which all working operations of the HE were investigated. (3) The regenerative mode of the HEs is mentioned and the difficulty in designing power-saving transmission for regeneration is explained. (4) The dynamic model of the entire system comprising the HE and integrated power sources were derived in detail. This model was simulated in a co-simulation AMESim-MATLAB/Simulink environment. The HE model was simulated in the AMESim software, whereas the models of FC, BAT, and SC were derived and performed in a MATLAB/Simulink software. The goal of this study was how to establish a real-time EMS, achieve the demand of the powertrain, and stabilize the entire platform when highly-fluctuating power occurred. (5) Finally, comparisons between the proposed algorithm with other conventional approaches are discussed to verify the effectiveness of the new configuration compared to previous conventional approaches.
The rest of this paper is organized as follows: Section 2 expresses the modeling of all components. Based on the characteristics of each device, the new configuration and the proposed EMS with regeneration mode are introduced in Section 3. Section 4 describes the MFLC based on the proposed EMS to achieve high efficiency and performance. To verify its advantages, comparative simulations between the proposed and conventional EMSs are discussed in Section 5. Finally, Section 6 summarizes and presents potential future applications.

Hydraulic Excavator Configuration
The power demand for running the HE is simply based on the total power needed for moving each element as requirements, driving crawlers, and running the hydraulic circuit. In the entire hydraulic system as depicted in Figure 1, the total power is a product of the outlet pressure and the flow rate of the hydraulic pump [39]: where p, D, n, and η are the pressure (bar), displacement (l/rev), rotational speed (rev/min), and volumetric efficiency of the hydraulic pump, respectively. P M is the computed motor power to drive the hydraulic pump (kW); η M is an efficiency dropped when converting electric power to mechanical power. The hydraulic circuit of the HE is inherited from the AMESim library and described in Figure 1 [40].

Fuel Cell Modeling
The PEMFC model is referred to in [41][42][43][44][45]. The dynamics of a single fuel cell is calculated as follows: where Nernst E , a c t V , conc V , and o h m ic V are voltage losses of the thermodynamic potential, activation process, concentration, and ohmic voltage loss, respectively.
where Rd is the activation resistance and concentration resistance, and Vd is the voltage drop.

Fuel Cell Modeling
The PEMFC model is referred to in [41][42][43][44][45]. The dynamics of a single fuel cell is calculated as follows: where E Nernst , V act , V conc , and V ohmic are voltage losses of the thermodynamic potential, activation process, concentration, and ohmic voltage loss, respectively. p H 2 , p O 2 , T, R, and F represent the hydrogen and oxygen partial pressure, cell temperature, universal gas constant, and Faraday constant, respectively. R int is the internal resistance of the electrolyte membrane; c O2 is the oxygen concentration at the cathode/membrane interface; ξ 1 . . . ξ 4 are parametric coefficients; and i is the cell current [38].
Due to the double capacitor layer effect at the electrode-electrolyte interface, the voltage drop can be computed as [46] dV d dt where R d is the activation resistance and concentration resistance, and V d is the voltage drop. The voltage of a single cell and total voltage when combining the number of cells N can be computed as V stack = NV cell (11) Energies 2020, 13, 3387 6 of 27 On the anode-side, the reactant flow model is given by the following: . .
m H 2 ,out are the anode volume, hydrogen inlet, and outlet flow rates through the fuel cell (FC) stack, respectively. k a is a flow constant for the anode, and p tan k is the pressure of the hydrogen tank. V c , m O 2 ,out are the cathode volume, oxygen inlet, and oxygen outlet flow rate through the FC stack, respectively. k c and p BPR are the flow constant in the cathode and oxygen pressure at the outlet, respectively. P tot denotes a total power input of the system, which changes linearly with hydrogen consumed, and ∆H is the hydrogen enthalpy of combustion.
Finally, the electrical output power can be obtained as

Supercapacitor Modeling
The SC was selected as the first support unit to supply different power between the demand and power generated from the FC and BATs in the case of sudden change due to its advantages such as high power density, fast charge, and high release power. Without using a SC, the FC and BATs have to manage the entire required workload even when high peak power happens, thus degrading lifespan, or increasing the size and costs in a trade-off [47]. Moreover, as Phatiphat stated in [48], when the time constant is less than 0.1 s, energy cannot be obtained from the same sized BAT, but the SC can totally provide energy at a very high rate. This is the highlighted advantage of the SC. Many different models have been suggested using the RC circuit. The equivalent SC model can be referred to in [48][49][50]. Naturally, the SC consists of capacitors and resistors that represent the charging and discharging units, and an equivalent parallel resistor is the self-discharging loss [51]. The unit cell of the SC was constructed with two RC branches in a parallel manner, as presented in Figure 2 [52].
On the anode-side, the reactant flow model is given by the following:

Supercapacitor Modeling
The SC was selected as the first support unit to supply different power between the demand and power generated from the FC and BATs in the case of sudden change due to its advantages such as high power density, fast charge, and high release power. Without using a SC, the FC and BATs have to manage the entire required workload even when high peak power happens, thus degrading lifespan, or increasing the size and costs in a trade-off [47]. Moreover, as Phatiphat stated in [48], when the time constant is less than 0.1 s, energy cannot be obtained from the same sized BAT, but the SC can totally provide energy at a very high rate. This is the highlighted advantage of the SC. Many different models have been suggested using the RC circuit. The equivalent SC model can be referred to in [48][49][50]. Naturally, the SC consists of capacitors and resistors that represent the charging and discharging units, and an equivalent parallel resistor is the self-discharging loss [51]. The unit cell of the SC was constructed with two RC branches in a parallel manner, as presented in Figure 2 [52].

R1
R2  As depicted in Figure 2, the main cell, R 1 C 1 , specifies the immediate response during the charge or discharge process within a short time duration [53]. All charged energy is in the capacitor C 1 of the immediate branch at the end of the charging mode. Then, the charge re-splits itself to the second slow cell, R 2 C 2 . The resistor R f denotes a leakage behavior based on the fact that there always exists Energies 2020, 13, 3387 7 of 27 a leakage current flow of the SC during a self-discharge phenomenon [54]. This current leakage is approximately several milliamps in a large SC. The model of SC is derived as where U sc and I sc are the voltage and current of the pack SCs; and v sc and i sc are the voltage and current of an elementary SC. N s_sc and N p_sc denote the number of SCs in a serial connection and the number of parallel branches, respectively. i 1 and i 2 are the current through the main cell and the slow cell. Q 1 and C 1 are the instantaneous charge state of the main cell. v 1 , v 2 is the voltage of the first and secondary branches. Q 2 and C 2 are the instantaneous charge states of the slow cell. Finally, the SoC SC is an important parameter to evaluate the state of the SC bank for designing the EMS. The change rate of the SoC SC is proportional with the charging current i sc [55]: where Q SCmax denotes the maximum capacity of the SC.

Battery Modeling
The BAT is the second buffer supply for the system when the FC and/or the SC cannot maintain the workload demand during the endurance process. The BAT model is constructed with specific invoked parameters to evaluate its working status, thereby establishing criteria for effectively initializing the EMS. The BAT model used in this study follows the circuit proposed in [56].
where E is the controlled voltage source; V BAT and P BAT are the BAT voltage and output power, respectively; Q BAT and Q BATmax denote the instant and maximum BAT capacity, respectively; i BAT is the load current; and t is time parameterized. A BAT is an exponential zone amplitude (V).
The SoC BAT can be derived from the current charge and the maximum charge of it.
Energies 2020, 13, 3387 The SoC BAT is another important parameter, along with the SoC SC , that reveals the device status and is invoked to assess the available energy remaining. The direction of the current determines the charge or discharge mode of the BAT. Compared with the BAT model in [47,57,58], this model gives the shortened form and can be used in both charge and discharge mode without losing its characteristics.

DC/DC Converter Modeling
For the studied hybrid electric excavator system, the BAT and SC were interfaced with the DC bus by using a bidirectional DC/DC converter, which can allow energy flow in both directions. The PEMFC system was linked with the DC bus via a buck DC/DC converter, which can adapt the voltage level of the DC bus. These DC/DC converters are the crucial execution devices of the energy management process, which can regulate an output voltage and current using a local controller. The energy management process is at a higher level than a local controller of the DC/DC converters and it is assumed that the time constant of the inductors in the DC/DC converters are much greater than the switching period [59]. Therefore, once the inner-loop subsystem is well controlled, it can respond immediately to the reference. Thus, it is reasonable to reduce the fast dynamics of the DC/DC converter using the following equivalent static model [60,61].
where V I , V O are the DC/DC converter input and output voltage, respectively; R L is the resistor of the inductor and L is its inductance; κ DC is the ratio of converter output and input voltage; i L and i O are currents through the inductor and the output current of the converter, respectively; η DC is the converter efficient; ε = 1 for boost mode or the bidirectional mode with i O V O ≥ 0 and ε = -1 for bidirectional mode when i O V O < 0.

Hybrid Power Hydraulic Excavators Configuration
The hybrid power HE (HP-HE) was constructed as depicted in Figure 3a, in which the ICEs were alternatively replaced by the hybrid FC-BAT-SC power source and electric motor. The entire system can be considered as the HPS associated with the hydraulic actuators (HAs) through the DC bus. The HPS comprises the PEMFC functioning as the primary supply and the BAT-SC functioning as auxiliary buffers. The goal of the HPS is to supply power as the load power required and store energy in the case of regeneration. The power induced from the HPS provides power to the electric motor through DC/DC converters, and the electric motor drives the practical hydraulic system. The HAs consisted of three hydraulic cylinders for the driving boom, arm, and bucket, and one swing hydraulic motor for driving the HE body. The movement of the three cylinders was distinguished by using a main control valve block to control the flow rate to each actuator. Furthermore, the HAs were extended by the regeneration part to recover the energy-saving when the boom automatically moves down due to gravity, as shown in Figure 3b.
Energies 2020, 13, x FOR PEER REVIEW 9 of 27 extended by the regeneration part to recover the energy-saving when the boom automatically moves down due to gravity, as shown in Figure 3b.

Proposed Energy Management Strategy (EMS)
Based on the fact that HEs perform various functions such as digging, lifting soil, spinning, ground leveling, and so on, an EMS was proposed to productively achieve the requirements. Improving the system performance and considering regeneration mode for saving energy were taken into consideration. From the requested motions of the excavator, the power required (Pred) was determined as a reference for running the power from the HPS. The SoCBAT and the SoCSC were used to evaluate the state of the components. Furthermore, the minimum SoC signed by SoCBAT_min and SoCSC_min was used to monitor when the two devices needed to be charged. The flowchart of the proposed EMS is depicted in Figure 4.
As depicted in Figure 4, in the beginning, the power required was examined to determine if its value was negative. In this case, the HE is moving downhill, when the swing is braking after spinning the entire body, or when the boom is moving down, whereas other elements do not operate. When these events happen, the BAT-SC is not required to release power. Instead, depending on the level of the power-saving, the FC is switched OFF or runs in an optimal value in which the highest efficiency can be achieved to charge the BAT-SC.
Otherwise, when Preq is positive, the following circumstances should be considered. If the power demand is less than the nominal power of the FC (Pfc_n), the FC primarily supplies power to the system and charges the auxiliary supplements for later use. The charging is stopped if the SoCBAT/SoCSC reaches the maximum value (SoCBAT/SoCSC = SoCMax). If the load demand exceeds the FC nominal power, the combination between the FC and BAT or SC is considered, if their power can sufficiently satisfy the load demand. This scenario is highlighted in light blue and light green in Figure 4. The combination of FC-SC is considered in the case when the BAT cannot satisfy the sudden change in power required, as explained in the SC modeling part. Otherwise, the prior combination of FC-BAT is used. During this process, if the SoCBAT (or SoCSC) drops to the minimum level, the SC (or BAT) is alternatively used, and the other device is alarmed and switched to charge mode. If the Preq exceeds the combination of FC-BAT (or FC-SC), all devices enter together and provide power to the system. During these processes, the SoCBAT and SoCSC should be monitored and the SC is always charged first.

Proposed Energy Management Strategy (EMS)
Based on the fact that HEs perform various functions such as digging, lifting soil, spinning, ground leveling, and so on, an EMS was proposed to productively achieve the requirements. Improving the system performance and considering regeneration mode for saving energy were taken into consideration. From the requested motions of the excavator, the power required (P red ) was determined as a reference for running the power from the HPS. The SoC BAT and the SoC SC were used to evaluate the state of the components. Furthermore, the minimum SoC signed by SoC BAT_min and SoC SC_min was used to monitor when the two devices needed to be charged. The flowchart of the proposed EMS is depicted in Figure 4.
As depicted in Figure 4, in the beginning, the power required was examined to determine if its value was negative. In this case, the HE is moving downhill, when the swing is braking after spinning the entire body, or when the boom is moving down, whereas other elements do not operate. When these events happen, the BAT-SC is not required to release power. Instead, depending on the level of the power-saving, the FC is switched OFF or runs in an optimal value in which the highest efficiency can be achieved to charge the BAT-SC.
Otherwise, when P req is positive, the following circumstances should be considered. If the power demand is less than the nominal power of the FC (P fc_n ), the FC primarily supplies power to the system and charges the auxiliary supplements for later use. The charging is stopped if the SoC BAT /SoC SC reaches the maximum value (SoC BAT /SoC SC = SoC Max ). If the load demand exceeds the FC nominal power, the combination between the FC and BAT or SC is considered, if their power can sufficiently satisfy the load demand. This scenario is highlighted in light blue and light green in Figure 4. The combination of FC-SC is considered in the case when the BAT cannot satisfy the sudden change in power required, as explained in the SC modeling part. Otherwise, the prior combination of FC-BAT is used. During this process, if the SoC BAT (or SoC SC ) drops to the minimum level, the SC (or BAT) is alternatively used, and the other device is alarmed and switched to charge mode. If the P req exceeds the combination of FC-BAT (or FC-SC), all devices enter together and provide power to the system. During these processes, the SoC BAT and SoC SC should be monitored and the SC is always charged first.  Figure 4. The proposed strategy for energy management strategy (EMS) with the regeneration part.

Regeneration Mode
In electric vehicles (EVs), the charge mode occurs when the cars are moving downhill or braking, while the FC is still supplying power to the system. One important point to clarify is that the charge and discharge mode cannot simultaneously occur as analyses in conventional algorithms. Since the BAT-SC is connected with only one engine and the engine is directly coupled with two-back-wheels, by using bi-directional DC/DC converters in the hybrid source, the state of the devices can be easily switched to capture or release power, depending on the state of the system and can be used in the same transmission line.
In contrast to the EVs, the powertrain of the HEs consists of various links (boom, arm, bucket, swing, and crawlers). Therefore, designing the regeneration mode in the HEs is more complicated than that of the EVs. We need to consider a case when one link is in a power-saving procedure (for instance, the boom is moving down automatically due to gravity), while the others continuously operate as usual. This implies that the power sources supplying progress for the system and powersaving progress from one or more elements occur at the same time. Hence, we cannot apply the same design as that of the EVs. For general hydraulic circuits, one sub-system should be connected as the second line for saving energy. Consequently, in this design, the circuit diagram for hybrid sources is

Regeneration Mode
In electric vehicles (EVs), the charge mode occurs when the cars are moving downhill or braking, while the FC is still supplying power to the system. One important point to clarify is that the charge and discharge mode cannot simultaneously occur as analyses in conventional algorithms. Since the BAT-SC is connected with only one engine and the engine is directly coupled with two-back-wheels, by using bi-directional DC/DC converters in the hybrid source, the state of the devices can be easily switched to capture or release power, depending on the state of the system and can be used in the same transmission line.
In contrast to the EVs, the powertrain of the HEs consists of various links (boom, arm, bucket, swing, and crawlers). Therefore, designing the regeneration mode in the HEs is more complicated than that of the EVs. We need to consider a case when one link is in a power-saving procedure (for instance, the boom is moving down automatically due to gravity), while the others continuously operate as usual. This implies that the power sources supplying progress for the system and power-saving progress from one or more elements occur at the same time. Hence, we cannot apply the same design as that of the EVs. For general hydraulic circuits, one sub-system should be connected as the second line for saving energy. Consequently, in this design, the circuit diagram for hybrid sources is separated into two distinguished transmission lines: one for supply and one for regeneration. The power captured from the generation is transmitted through the second transmission line.

Fuzzy EMS for the Integrated System
According to the working principle presented by the above flow chart, the power obtained from the driving cycle can be separated into four cases: high power, medium power, low power, and regenerative mode. The conventional EMS can simply be considered as an on-off switch between devices. However, this method may generate a delay to the system, as explained in the Introduction. In practice, the change in the status of devices requires a certain time to achieve and adapt to the required tasks. This problem can be handled by employing a FLC, which is known to be a useful tool to split and distribute power. This method was previously suggested by [18]. In this approach, the prior SoC SC is retained so that the stability of the SoC SC is maintained in an acceptable condition. However, in practice, maintaining the SoC SC in an acceptable condition is not as important as that of BAT due to its specification. In large systems that always work under high power for an extended period, the BAT operation and its status should be regulated so that its SoC exhibits variation or is prevented from dropping down to a very low value as a solution to prolong life-time. The efficiency curve of the PEMFC is shown in Figure 5.
BAT due to its specification. In large systems that always work under high power for an extended period, the BAT operation and its status should be regulated so that its SoC exhibits variation or is prevented from dropping down to a very low value as a solution to prolong life-time. The efficiency curve of the PEMFC is shown in Figure 5.
As shown in Figure 5, efficiency does not increase as power increases. High power achievement requires more energy from the auxiliary devices (cooling or heating system), thus reducing the overall efficiency of the system. Therefore, designing an algorithm to maintain high efficiency when operating the FC is practically prioritized as the main criterion. Thereby, we can reduce hydrogen consumption and limit the on/off cycles of the system to extend lifespan. The FLS was set up with the following criteria to satisfy the above analysis: • Increase the FC efficiency and minimize hydrogen consumption.

•
The SoCBAT should be frequently maintained within the range of 0.5~0.9 as a solution to prolong lifespan. • The SoCSC should be maintained at a high level to boost the power in the case of an emergency.
Therefore, the fuzzy control rules can serve to map the input linguistic variables Preq, SoCBAT, and   As shown in Figure 5, efficiency does not increase as power increases. High power achievement requires more energy from the auxiliary devices (cooling or heating system), thus reducing the overall efficiency of the system. Therefore, designing an algorithm to maintain high efficiency when operating the FC is practically prioritized as the main criterion. Thereby, we can reduce hydrogen consumption and limit the on/off cycles of the system to extend lifespan. The FLS was set up with the following criteria to satisfy the above analysis: • Increase the FC efficiency and minimize hydrogen consumption.

•
The SoC BAT should be frequently maintained within the range of 0.5~0.9 as a solution to prolong lifespan.

•
The SoC SC should be maintained at a high level to boost the power in the case of an emergency.
Therefore, the fuzzy control rules can serve to map the input linguistic variables P req , SoC BAT , and SoC SC to the output linguistic variables P fc , P BAT , and P SC , respectively. The input linguistic  Tables 1-4. The membership function of the inputs and outputs are described in Figure 6.  For the regeneration mode, because the excavator does not require any power from the power supply, the power of the FC is adjusted to operate at maximum efficiency in the case of charging auxiliary devices when the SoC drops down to a low value and off when the SoC level is high. Furthermore, the power of BAT and SC power are assigned as outputs of the FLC. The fuzzy rule of the regenerative mode is presented in Table 4. For the regeneration mode, because the excavator does not require any power from the power supply, the power of the FC is adjusted to operate at maximum efficiency in the case of charging auxiliary devices when the SoC drops down to a low value and off when the SoC level is high. Furthermore, the power of BAT and SC power are assigned as outputs of the FLC. The fuzzy rule of the regenerative mode is presented in Table 4.  In this study, triangular membership functions are exploited as depicted in Figure 6. The left side illustrates the inputs of the FLS in which the power required P req , SoC BAT , and SoC SC are considered; the right side is the product including the FC power and gains of the BAT for charging or discharging mode. Additionally, to match the aforementioned criteria, the four following requirements should be satisfied:

•
In the case of low devices SoC (SoC BAT and SoC SC are low), the FC power can be set up to a high value for quick charging, even when the system is operating with medium or low power required. • For medium and low power required, if the SoC BAT is greater than medium level, the BAT charges the SC instead of using FC; therefore, the FC does not need to run at a high value, and the efficiency can be increased consequently.

•
In the case of charging, the SC is always charged so that a good condition of the SoC SC can be maintained for later use.

•
The final goal is to force the FC power to the highest efficiency point, as shown in Figure 5.
The completed control schematic diagram is depicted in Figure 7.
Since the power of the system in some cases can be referred to as the mixture between two cases, (for instance, if the required power is 85 kW, then this value can belong to both high power and medium power), then the scheduling gains is used to calculate the exact output power of the FLC.
Let us define the vector ( ) T K α β γ σ = as the gains of the high power, medium power, low power, and regenerative power, respectively. Then, the output power is calculated as: Consequently, the FC power used for the net is computed as

Parameters Setup for Simulation
In this section, the simulations and comparisons between the proposed control strategy with previous approaches are examined to evaluate the effectiveness of the proposed HPS under different Since the power of the system in some cases can be referred to as the mixture between two cases, (for instance, if the required power is 85 kW, then this value can belong to both high power and medium power), then the scheduling gains is used to calculate the exact output power of the FLC. Let us define the vector K = α β γ σ T as the gains of the high power, medium power, low power, and regenerative power, respectively. Then, the output power is calculated as: where P * FC = P * FC1 P * FC2 P * FC3 P *

FC4
T is an output of the four FLC, respectively, K = α β γ σ T is calculated as: Consequently, the FC power used for the net is computed as

Parameters Setup for Simulation
In this section, the simulations and comparisons between the proposed control strategy with previous approaches are examined to evaluate the effectiveness of the proposed HPS under different working conditions. The simulations were conducted in a co-simulation between LMS AMESim 15.2 software and MATLAB/Simulink 2019a as depicted in Figure 8. The LMS AMESim, known as specialized software for dynamic simulation with various applications, was employed to express the HE configuration. The parameters for simulating the HE are described in Table 5. The modeling of the HPS was deployed in MATLAB/Simulink, in which all mathematical equations were embedded, and the results were displayed with a sampling time of 10 ms. The component sizes of the HPS were reasonably chosen to satisfy the demand and verify the proposed methodology as listed in Tables 6-8. The testing trajectories of the boom, arm, and bucket of the excavator are described in Figure 9, and the required power from those motions is performed in Figure 10.  Table 7. Supercapacitor parameters [52]. The LMS AMESim, known as specialized software for dynamic simulation with various applications, was employed to express the HE configuration. The parameters for simulating the HE are described in Table 5. The modeling of the HPS was deployed in MATLAB/Simulink, in which all mathematical equations were embedded, and the results were displayed with a sampling time of 10 ms. The component sizes of the HPS were reasonably chosen to satisfy the demand and verify the proposed methodology as listed in Tables 6-8. The testing trajectories of the boom, arm, and bucket of the excavator are described in Figure 9, and the required power from those motions is performed in Figure 10.    Regarding the working cycles in Figure 9, the characteristic line sloped down when the element moved down (i.e., the cylinder retraces and increases when the element moves up, i.e., the cylinder is extending. With the driving cycles shown in Figure 9, changes in the characteristic lines result in changes in the power required, as depicted in Figure 10, for instance, referring to the driving cycle (a) at the time of the 5th to 10th second when all three elements are moving, the power required in Figure 10 increases. Through the power curve, the power needed to run the boom was the highest compared with that of the other two; in particular, the power curve increased significantly when the boom moved up. After that, at the time of the 25th second, all elements were kept in a fixed position and the power decreased. Therefore, the buffer supplements were charged at these times as in this scenario. Regarding the working cycles in Figure 9, the characteristic line sloped down when the element moved down (i.e., the cylinder retraces and increases when the element moves up, i.e., the cylinder is extending. With the driving cycles shown in Figure 9, changes in the characteristic lines result in changes in the power required, as depicted in Figure 10, for instance, referring to the driving cycle (a) at the time of the 5th to 10th second when all three elements are moving, the power required in Figure 10 increases. Through the power curve, the power needed to run the boom was the highest compared with that of the other two; in particular, the power curve increased significantly when the boom moved up. After that, at the time of the 25th second, all elements were kept in a fixed position and the power decreased. Therefore, the buffer supplements were charged at these times as in this scenario.

Simulation Results and Discussion
In this section, the two simulations comprising the proposed EMS in [38] and the proposed fuzzy EMS (F-EMS), as expressed in Section 3, were conducted and compared with the three other strategies: (1) conventional EMS (C-EMS) [62], (2) HFS-CS, and (3) LCS-CS in [16].
The responses of the total power supply, FC, BAT, and SC are described in Figures 11-14, respectively. In each figure, the responses when embedding the proposed EMS, proposed F-EMS, C-EMS, HFS-CS, and LCS-CS are depicted in sub-figures (a), (b), (c), (d), and (e), respectively. The total power supplied from the proposed algorithm is depicted in Figure 11a and proposed F-EMS is depicted in Figure 11b. The other results when embedding the C-EMS, HFS-CS, and LCS-CS into the system are depicted in Figure 11c,d,e, respectively. These results reveal that the proposed strategies can provide sufficient power to the system and maintain stable levels during an endurance process; consequently, good performance can be achieved.
In contrast, the other EMS strategies were not capable of maintaining sufficient energy. For this reason, the C-EMS did not consider the charging process for auxiliary devices when their SoCBAT/SoCSC dropped to a low level. When the SoCBAT/SoCSC reached the minimum level as shown in Figure 15, the supplements were out of power and could not support the system.

Simulation Results and Discussion
In this section, the two simulations comprising the proposed EMS in [38] and the proposed fuzzy EMS (F-EMS), as expressed in Section 3, were conducted and compared with the three other strategies: (1) conventional EMS (C-EMS) [62], (2) HFS-CS, and (3) LCS-CS in [16].
The responses of the total power supply, FC, BAT, and SC are described in Figures 11-14 Figure 11a and proposed F-EMS is depicted in Figure 11b. The other results when embedding the C-EMS, HFS-CS, and LCS-CS into the system are depicted in Figure 11c-e, respectively. These results reveal that the proposed strategies can provide sufficient power to the system and maintain stable levels during an endurance process; consequently, good performance can be achieved.
In contrast, the other EMS strategies were not capable of maintaining sufficient energy. For this reason, the C-EMS did not consider the charging process for auxiliary devices when their SoC BAT /SoC SC dropped to a low level. When the SoC BAT /SoC SC reached the minimum level as shown in Figure 15, the supplements were out of power and could not support the system.
As a result, the system can not complete the requirement due to the shortage of power, starting at the time of the 140th second. For the two HFS-CS and LCS-CS strategies, charging procedures were considered in which the positive power of the auxiliary devices indicates that they are in discharge mode, whereas the negative value indicates that they are in charge mode. However, instability occurs due to the inappropriate order of using devices. Starting at the time of the 210th second, when the SoC BAT /SoC SC dropped down to the minimum level, the chattering phenomenon took place because the devices immediately switched from releasing to charging mode when the SoC BAT /SoC SC hit the minimum level. Subsequently, when the SoC BAT /SoC SC is greater than the minimum level, the states of the devices instantly shift from charging to releasing the power, despite low remaining power. These processes repeatedly occur and generate high-frequency fluctuations. Therefore, these paradigms need to be noticeably improved based on the fact that either the system has to stop for charging like in conventional algorithms, or employ another strategy to maintain performance.     As a result, the system can not complete the requirement due to the shortage of power, starting at the time of the 140 th second. For the two HFS-CS and LCS-CS strategies, charging procedures were considered in which the positive power of the auxiliary devices indicates that they are in discharge mode, whereas the negative value indicates that they are in charge mode. However, instability occurs due to the inappropriate order of using devices. Starting at the time of the 210th second, when the SoCBAT/SoCSC dropped down to the minimum level, the chattering phenomenon took place because the devices immediately switched from releasing to charging mode when the SoCBAT/SoCSC hit the minimum level. Subsequently, when the SoCBAT/SoCSC is greater than the minimum level, the states of the devices instantly shift from charging to releasing the power, despite low remaining power. These processes repeatedly occur and generate high-frequency fluctuations. Therefore, these paradigms need to be noticeably improved based on the fact that either the system has to stop for charging like in conventional algorithms, or employ another strategy to maintain performance.
Compared with the three previous strategies, the proposed algorithm, which considers all circumstances for stable power-sharing, has greater endurance and the supplements are considerably maintained in good condition to extend the duration for the requirements. As expressed in Figure 15, the SoCBAT varied negligibly from 0.769 at the start to 0.7685 at the end of the process. The SoCSC had a periodic oscillation around the value of 0.6 with a high magnitude due to the high power required. Moreover, for the proposed EMS, the Pfc was maintained at constant power and was only switched off in the case that no power was required. Consequently, the BAT and SC are charged every time when high power is not required. However, this is not an optimal solution for the EMS since overcharge happened in the SC at the time of the 50th second, as shown in Figure 15. Moreover, the SoCBAT increased from the beginning to the end of the process because when the SC reached the maximum level, no additional power could be captured and the excess power from the FC was transferred to the BAT. Due to maintaining the constant level, the FC required more hydrogen consumption and the FC efficiency was low as a result. This issue was resolved by using the F-EMS employing the MFLC as explained in Section 4. Based on the comparison between the two proposed EMSs, the FC in the F-EMS consumes less fuel than that of the proposed EMS and is forced to the rated power to increase efficiency. The BAT takes the place of the FC in charging energy for the SC and maintaining it in good condition, instead of using the FC all the time.
Referring to the FC model in Section 2, the efficiency of the FC can be calculated as [63] max 0,1 elec aux tot elec The comparison of the efficiency for the two proposed algorithms is depicted in Figure 16. Compared with the three previous strategies, the proposed algorithm, which considers all circumstances for stable power-sharing, has greater endurance and the supplements are considerably maintained in good condition to extend the duration for the requirements. As expressed in Figure 15, the SoC BAT varied negligibly from 0.769 at the start to 0.7685 at the end of the process. The SoC SC had a periodic oscillation around the value of 0.6 with a high magnitude due to the high power required. Moreover, for the proposed EMS, the P fc was maintained at constant power and was only switched off in the case that no power was required. Consequently, the BAT and SC are charged every time when high power is not required. However, this is not an optimal solution for the EMS since overcharge happened in the SC at the time of the 50th second, as shown in Figure 15. Moreover, the SoC BAT increased from the beginning to the end of the process because when the SC reached the maximum level, no additional power could be captured and the excess power from the FC was transferred to the BAT. Due to maintaining the constant level, the FC required more hydrogen consumption and the FC efficiency was low as a result. This issue was resolved by using the F-EMS employing the MFLC as explained in Section 4. Based on the comparison between the two proposed EMSs, the FC in the F-EMS consumes less fuel than that of the proposed EMS and is forced to the rated power to increase efficiency. The BAT takes the place of the FC in charging energy for the SC and maintaining it in good condition, instead of using the FC all the time.
Referring to the FC model in Section 2, the efficiency of the FC can be calculated as [63] E f f = P elec P tot × max 0, 1 − P aux P elec The comparison of the efficiency for the two proposed algorithms is depicted in Figure 16. As revealed in Figure 16, the maximum efficiency that the FC can reach when applying the proposed EMS was approximately 37% due to maintaining the constant value of 65 kW. Conversely, the efficiency when embedding the F-EMS was higher than that of the proposed EMS due to reducing the hydrogen consumption. With the response of the FC shown in Figure 12a, the efficiency of the FC could reach up to 47%. The average efficiency of the FC increased from 32% when using the proposed EMS to 42% when applying the F-EMS under the same conditions. Furthermore, the number of As revealed in Figure 16, the maximum efficiency that the FC can reach when applying the proposed EMS was approximately 37% due to maintaining the constant value of 65 kW. Conversely, the efficiency when embedding the F-EMS was higher than that of the proposed EMS due to reducing the hydrogen consumption. With the response of the FC shown in Figure 12a, the efficiency of the FC could reach up to 47%. The average efficiency of the FC increased from 32% when using the proposed EMS to 42% when applying the F-EMS under the same conditions. Furthermore, the number of ON/OFF switches decreases as a solution for extending the FC lifespan.

Conclusions
The new configuration for the hydraulic excavators equipped with a hybrid power source comprising the fuel cell, batteries, and supercapacitors with the proposed EMS is the merit of this paper. Instead of using one auxiliary supplement until it runs out of power and using the other as a compensation supply later, the EMS was exploited to distribute the energy from the powertrain to each device so that good performance is achieved and the lifespan of all components can be extended. The simulation results displayed that the proposed EMS can satisfy the requirements in comparison with other previous approaches under different working conditions. Moreover, the generation procedure was considered as the time for healing buffer supplements. In a going effect, the use of the proposed F-EMS confirmed its effectiveness when reducing the fuel consumption in the FC while maintaining the SC in a good condition for use in the event of an emergency. The achievements in this study can not only guarantee power supply performance and prolong the lifespan of all devices, but also reduce the fuel consumption and increase the efficiency of the system. Therefore, this study is a premise for developing a hybrid power source for construction machinery in the future.

Conflicts of Interest:
The authors declare no conflict of interest.