An Optimal Slip Ratio-Based Revised Regenerative Braking Control Strategy of Range-Extended Electric Vehicle

The energy recovered with regenerative braking system can greatly improve energy efficiency of range-extended electric vehicle (R-EEV). Nevertheless, maximizing braking energy recovery while maintaining braking performance remains a challenging issue, and it is also difficult to reduce the adverse effects of regenerative current on battery capacity loss rate (Qloss,%) to extend its service life. To solve this problem, a revised regenerative braking control strategy (RRBCS) with the rate and shape of regenerative braking current considerations is proposed. Firstly, the initial regenerative braking control strategy (IRBCS) is researched in this paper. Then, the battery capacity loss model is established by using battery capacity test results. Eventually, RRBCS is obtained based on IRBCS to optimize and modify the allocation logic of braking work-point. The simulation results show that compared with IRBCS, the regenerative braking energy is slightly reduced by 16.6% and Qloss,% is reduced by 79.2%. It means that the RRBCS can reduce Qloss,% at the expense of small braking energy recovery loss. As expected, RRBCS has a positive effect on prolonging the battery service life while ensuring braking safety while maximizing recovery energy. This result can be used to develop regenerative braking control system to improve comprehensive performance levels.


Introduction
With increasingly prominent energy and environmental issues, electric vehicles have developed rapidly in recent years due to their advantages of low pollution-free emission. However, due to the low energy density of battery, short driving range and long charging time, the marketing of battery electric vehicles (BEV) is greatly limited [1,2]. Range-extended electric vehicles (R-EEVs) add auxiliary power units, it can solve the range anxiety from consumers by increasing driving range. In addition, a R-EEV has a certain commercial competitiveness in the future market [3]. The regenerative braking system (RBS) provides an effective way to greatly improve overall electrical performance of R-EEVs [4,5]. Thereby, the regenerative braking control strategy (RBCS) is worthwhile research, and is receiving a good deal of attention.
The existing research has focused on RBCS in order to achieve better capacity of the regenerative braking energy. Literature [6] has analyzed the braking force distribution strategy and studied the energy saving potential. Different control strategies are studied to achieve the goals of regeneration efficiency and braking safety in [7][8][9]. Similar researches also include Literature [10][11][12][13], which develop braking torque allocating method to achieve improved braking performance and energy regeneration. In terms of the currently available approaches, it is mainly focused on how to coordinate recovery while ensuring the braking efficiency. It should be noted that the impact of regenerative braking current (RCC) on battery health should also be given sufficient attention and research.
The other focus is optimization and the analysis of braking torque distributions strategies under extreme adhesion conditions or an emergency braking process. A revised control strategy is presented with the front wheel slip ratio consideration in [14] to prevent a front wheel lock-up and maximize the regenerative braking efficiency under low tire-road friction conditions. Literature [15] proposed an efficient RBCS based on the modified nonlinear model predictive control method to ensure braking safety. In addition, a study team [16] built a high-precision predictive model based on the off-line optimization data of the combined model to solve the poor real-time problem of the optimization. The artificial neural network-based control mechanism was utilized to optimize the switching scheme of the vehicular braking force distribution was proposed [17].
The previous studies have focused more on the braking efficiency and maximization of regenerative energy, but the impact of RCC on battery service life is not considered sufficiently. It means there is still much room for optimization of RBCS. Therefore, a revised regenerative braking control strategy (RRBCS) is proposed in this paper. The strategy based on optimal slip ratio (OSR) to ensure braking performance while maximizing braking energy recovery and reducing battery capacity loss rate (Qloss,%). The rest of this paper is organized as follows: In Section 2, the initial regenerative braking control strategy (IRBCS) is given, and the distribution of the front and rear slip ratios and allocation logic of the braking work-point are introduced. A simulation model is built in Section 3, and the tire-road adhesion coefficient recognition module (TACRM) and double fuzzy logic controller (DFLC) are verified by simulation test. In section 4, battery capacity loss module (BCLM) is established. The BCLM uses test results of battery life decay test. Eventually, RRCBS is proposed. In section 5, the effectiveness and feasibility of the proposed RRBCS are validated by test under the world light vehicle test procedure (WLTP). This is followed by the conclusion in the final section.

Control Principle Overview of IRBCS
The research object is a medium-sized R-EEV. As shown in Figure 1, the R-EEV is driven by the front axle. The regenerative braking system consists of two systems: hydraulic braking system and an electric braking system. Four wheels are equipped with regulating valves, so each wheel can be controlled independently by braking controller. The optimal control of braking efficiency and energy recovery cannot be guaranteed with fixed slip ratio on different tire-road adhesion conditions. Through monitoring and identification of tireroad adhesion conditions in real time, the OSR is output as the control target value to ensure braking performance [18,19]. Process of IRBCS based on OSR is shown in Figure 2.  The optimal control of braking efficiency and energy recovery cannot be guaranteed with fixed slip ratio on different tire-road adhesion conditions. Through monitoring and identification of tire-road adhesion conditions in real time, the OSR is output as the control target value to ensure braking performance [18,19]. Process of IRBCS based on OSR is shown in Figure 2. When braking starts, hydraulic brake control unit judges the battery state of charge (SoC). If SoC ≥0.8, it means battery energy is high, the RBS is not carried out and the wheels are braked by hydraulic braking system. If SoC<0.8 and Vvehic>0, the RBS start to work and braking torque distributed between front wheels with regenerative braking torque and rear wheels with hydraulic braking torque provided by the motor. If the regenerative braking torque can provide sufficient required braking torque, the front wheel braking torque is all provided by the regenerative braking torque. If the regenerative braking torque cannot provide sufficient required braking torque, the front wheel braking torque is provided by the regenerative braking torque and the hydraulic braking torque.

Tire-Road Adhesion Coefficient Recognition Module (TACRM)
Through TACRM, the tire-road adhesion coefficient (TAC) is estimated in real time and the realtime optimal slip ratio (sosr (t)) and real-time tire-road adhesion coefficient (μosr (t)) are outputted. Several mature tire-road models based on empirical and analytical can be selected [13]. A braking force distribution strategy is proposed in [19] based on the estimation of the tire-road friction coefficient using a fuzzy logic estimation approach. In this paper, the empirical model 'Burckhardt tire model' is selected because of its suitability for tire behavior simulation, as shown in Figure 3(a) and Table 1. The equation of the 'Burckhardt tire model' is expressed as follows:    When braking starts, hydraulic brake control unit judges the battery state of charge (SoC). If SoC ≥0.8, it means battery energy is high, the RBS is not carried out and the wheels are braked by hydraulic braking system. If SoC<0.8 and V vehic >0, the RBS start to work and braking torque distributed between front wheels with regenerative braking torque and rear wheels with hydraulic braking torque provided by the motor. If the regenerative braking torque can provide sufficient required braking torque, the front wheel braking torque is all provided by the regenerative braking torque. If the regenerative braking torque cannot provide sufficient required braking torque, the front wheel braking torque is provided by the regenerative braking torque and the hydraulic braking torque.

Tire-Road Adhesion Coefficient Recognition Module (TACRM)
Through TACRM, the tire-road adhesion coefficient (TAC) is estimated in real time and the real-time optimal slip ratio (s osr (t)) and real-time tire-road adhesion coefficient (µ osr (t)) are outputted. Several mature tire-road models based on empirical and analytical can be selected [13]. A braking force distribution strategy is proposed in [19] based on the estimation of the tire-road friction coefficient using a fuzzy logic estimation approach. In this paper, the empirical model 'Burckhardt tire model' is selected because of its suitability for tire behavior simulation, as shown in Figure 3a and Table 1. The equation of the 'Burckhardt tire model' is expressed as follows: Energies 2019, 12, x FOR PEER REVIEW 3 of 21 Figure 2. Process of regenerative braking control strategy based on optimal slip ratio.
When braking starts, hydraulic brake control unit judges the battery state of charge (SoC). If SoC ≥0.8, it means battery energy is high, the RBS is not carried out and the wheels are braked by hydraulic braking system. If SoC<0.8 and Vvehic>0, the RBS start to work and braking torque distributed between front wheels with regenerative braking torque and rear wheels with hydraulic braking torque provided by the motor. If the regenerative braking torque can provide sufficient required braking torque, the front wheel braking torque is all provided by the regenerative braking torque. If the regenerative braking torque cannot provide sufficient required braking torque, the front wheel braking torque is provided by the regenerative braking torque and the hydraulic braking torque.

Tire-Road Adhesion Coefficient Recognition Module (TACRM)
Through TACRM, the tire-road adhesion coefficient (TAC) is estimated in real time and the realtime optimal slip ratio (sosr (t)) and real-time tire-road adhesion coefficient (μosr (t)) are outputted. Several mature tire-road models based on empirical and analytical can be selected [13]. A braking force distribution strategy is proposed in [19] based on the estimation of the tire-road friction coefficient using a fuzzy logic estimation approach. In this paper, the empirical model 'Burckhardt tire model' is selected because of its suitability for tire behavior simulation, as shown in Figure 3(a) and Table 1. The equation of the 'Burckhardt tire model' is expressed as follows:     As shown in Figure 3b, the peak tire-road adhesion coefficient on j-th road (µ j (s(t))) can be calculated by substituting the real-time slip ratio (s(t)). The utilization adhesion coefficient (f(t)) is calculated as follows:

C1 C2 C3
where F x_adhesion is the wheel longitudinal adhesion force; F z_vertical is the wheel vertical force. When µ j (s(t))>f(t), the upper limit value of TAC (µ J_ul (t)) is calculated as follows: When µ j (s(t))≤f(t), the lower limit value of TAC (µ J_ll (t)) is calculated as follows: The upper and lower limit values of OSR (s J_osr_ul (t), s J_osr_ll (t)), upper and lower limit value of real-time tire-road adhesion coefficient (µ J_osr_ul (t), µ J_osr_ll (t)) on the j-th road are obtained using a look-up table. Then, the distribution proportion coefficient (k r ) is defined as follows: Real-time optimal slip ratio (s osr (t)) and corresponded adhesion coefficient (µ osr (t)) can be obtained as follows: s osr (t) = s J_osr_ul (t) + k r s J_osr_ul (t) − k r s J_osr_ll (t) (6) µ osr (t) = µ J_osr_ul (t) − k r µ J_osr_ul (t) + k r µ J_osr_ll (t)

Distribution of Front and Rear Slip Ratio
The braking working point is located on the ideal braking force distribution curve (I curve). Making the best use of the tire-road adhesion conditions can maximize the braking safety [20]. It is certain that the Economic Commission of Europe Regulation 13 needs to be considered for the braking force allocation to ensure braking safety in design process of RBCS. As shown in Figure 4a, the left intersection of the ECE regulation line (M curve) and the front wheel braking force axis (x-axis) is the maximum braking force value of working point A, which is also limited by the motor maximum torque (T mot_max ). It can ensure that all the braking working points obtained in RBCS are above the M curve to ensure braking safety.
When the braking working point is at point A, the braking strength is defined as the threshold of the low braking intensity (Z a_1 ), and it is calculated as follows: where m is the total mass of vehicle; R is the wheelbase; i is the transmission ratio. W is the correction factor with the tire-road adhesion condition consideration, and it is calculated as follows: where µ w_h is the threshold value of the high TAC; µ w_l is the threshold value of the low TAC; µ max is the peak tire-road adhesion coefficient (PTAC).

Distribution of Front and Rear Slip Ratio
The braking working point is located on the ideal braking force distribution curve (I curve). Making the best use of the tire-road adhesion conditions can maximize the braking safety [20]. It is certain that the Economic Commission of Europe Regulation 13 needs to be considered for the braking force allocation to ensure braking safety in design process of RBCS. As shown in Figure 4(a), the left intersection of the ECE regulation line (M curve) and the front wheel braking force axis (xaxis) is the maximum braking force value of working point A，which is also limited by the motor maximum torque (Tmot_max). It can ensure that all the braking working points obtained in RBCS are above the M curve to ensure braking safety.  When the required braking strength (Z req ) is less than Z a_l , the front wheel braking torque is provided by the regenerative braking torque totally. If the regenerative braking torque provides insufficient required braking torque, the motor remains at the maximum torque, and the insufficient braking force would be compensated by the hydraulic braking system until reaching point B on I curve. As shown in Figure 4a, When the braking working point is at point B, the braking strength is defined as the threshold of the high braking intensity (Z b_h ), which is calculated as follows: where m f is the front axle mass distribution coefficient; hg is the centroid height. If Z req > Z b_h , the front and rear wheel braking force distribution is performed according to the I curve, that is, the BC segment coincides with the I curve.
The ideal braking force distribution control is realized by controlling the front and rear slip ratio under braking process [6]. As shown in Figure 4b, the relationship between braking force and the slip ratio is derived according to the braking dynamics theory. The distribution values of front wheel slip rate slip ratio (s f ) and rear wheel slip rate slip ratio (s r ) in three braking stages are described as follows: (1) OA stage: full regenerative braking (2) AB stage: hydraulic braking compensation (3) BC stage: coordination braking on the I curve Energies 2020, 13, 1526 6 of 21

Allocation of Braking Work-Point
The allocation logic of braking work-point is a major part of IRBCS, and it determines the efficiency of IRBCS. The braking work-point in IRBCS is divided into five types in this paper, as shown in Figure 5. The allocation logic of braking work-point is a major part of IRBCS, and it determines the efficiency of IRBCS. The braking work-point in IRBCS is divided into five types in this paper, as shown in Figure 5. As shown in Figure 5 (a), μμ_h and μμ_l are the high tire-road adhesion threshold value and the low tire-road adhesion threshold value, respectively. The allocation logic of braking work-point in IRBCS is as follows: (1) Mild braking (Zreq_l< Za_l): when the PTAC is small(μmax<μμ_l), such as ice road, the wheel has the risk of locking and slipping in braking process. Therefore, the braking work-point should follow the I curve on work-point 1 (WC1) to ensure the braking efficiency. If the PTAC is better (μmax>μμ_l), the braking work-point should follow the "OA" curve on work-point 2 (WC2). The front wheel braking torque is provided by the motor and the rear wheel hydraulic force is 0. (2) Moderate braking (Zb_h>Zreq_m>Za_l), when the PTAC is not high enough (μmax<μμ_h), there is also a risk of locking slip due to the significant braking strength. Therefore, the braking work-point should follow the I curve on work-point 3 (WC3) to ensure the braking efficiency. If the PTAC is high (μmax>μμ_h), the braking work-point should follow the "AB" curve on work-point 4 (WC4) to achieve good braking efficiency while maximizing braking energy recovery. (3) In the emergency braking condition (Zreq_h>Zb_h), there is a risk of locking slip due to the significant braking strength. Therefore, the braking work-point should follow the I curve on work-point 5 (WC5) to ensure the braking efficiency. The rear wheel braking torque is provided by hydraulic braking system. The motor provides its maximum regenerative torque applied at the front wheels.

Simulation Model of R-EEV
In this section, simulation experiments are performed in MATLAB/Simulink and AVL/Cruise software (version 2014, AVL List GmbH, Austria) to evaluate the control effect of IRBCS. As shown in Figure 6, the vehicle dynamics model of R-EEV, including regenerative braking system is established in Cruise software. The vehicle performance simulation calculations are performed in Cruise software. The IRBCS control strategy is performed on the Simulink software through the API function software provided by Cruise. Signals such as traction and braking, battery SoC and Vvehic are transmitted between Cruise and Simulink. The main parameters used in the simulation are shown in Table 2.  As shown in Figure 5a, µ µ_h and µ µ_l are the high tire-road adhesion threshold value and the low tire-road adhesion threshold value, respectively. The allocation logic of braking work-point in IRBCS is as follows: (1) Mild braking (Z req_l < Z a_l ): when the PTAC is small(µ max <µ µ_l ), such as ice road, the wheel has the risk of locking and slipping in braking process. Therefore, the braking work-point should follow the I curve on work-point 1 (W C1 ) to ensure the braking efficiency. If the PTAC is better (µ max >µ µ_l ), the braking work-point should follow the "OA" curve on work-point 2 (W C2 ). The front wheel braking torque is provided by the motor and the rear wheel hydraulic force is 0. (2) Moderate braking (Z b_h >Z req_m >Z a_l ), when the PTAC is not high enough (µ max <µ µ_h ), there is also a risk of locking slip due to the significant braking strength. Therefore, the braking work-point should follow the I curve on work-point 3 (W C3 ) to ensure the braking efficiency. If the PTAC is high (µ max >µ µ_h ), the braking work-point should follow the "AB" curve on work-point 4 (W C4 ) to achieve good braking efficiency while maximizing braking energy recovery. (3) In the emergency braking condition (Z req_h >Z b_h ), there is a risk of locking slip due to the significant braking strength. Therefore, the braking work-point should follow the I curve on work-point 5 (W C5 ) to ensure the braking efficiency. The rear wheel braking torque is provided by hydraulic braking system. The motor provides its maximum regenerative torque applied at the front wheels.

Simulation Model of R-EEV
In this section, simulation experiments are performed in MATLAB/Simulink and AVL/Cruise software (version 2014, AVL List GmbH, Austria) to evaluate the control effect of IRBCS. As shown in Figure 6, the vehicle dynamics model of R-EEV, including regenerative braking system is established in Cruise software. The vehicle performance simulation calculations are performed in Cruise software. The IRBCS control strategy is performed on the Simulink software through the API function software provided by Cruise. Signals such as traction and braking, battery SoC and V vehic are transmitted between Cruise and Simulink. The main parameters used in the simulation are shown in Table 2

Verification of TACRM
The TACRM proposed is verified by simulation test in this section. A co-simulation test is carried out based on the braking system model of R-EEV in simulation environment. In simulations, the initial braking speed is set at 80 km/h, and four tire-road conditions were set with PTAC of 0.6, 0.45, 0.2, and 0.1, respectively, as shown in Figure 7. The braking strength increases from 0 to 0.7 quickly to simulate emergency braking conditions, and it stabilizes at 0.5 until the vehicle speed is 0. The estimated value of the PTAC will be output in real time through the TACRM. The effectiveness and accuracy of the TACRM can be seen by comparing with the preset value of PTAC. The interpolation calculation method is used to estimate the tire-road adhesion condition in the four types of tire-roads, and it can effectively and accurately estimate the PTAC in real time.

OSR Control Method By DFLC
The braking force distribution of the front and rear wheels is different. The membership function of the fuzzy controller and the formulation of the fuzzy rules are also different, which will affect the target slip ratio control result. Considering the difference of braking torque between front and rear wheels, the DFLC for front and rear wheels is designed to ensure OSR tracking control target values effectively, which is different from the fuzzy controllers in [9,13,[21][22][23][24][25]. The requirement slip ratio

Verification of TACRM
The TACRM proposed is verified by simulation test in this section. A co-simulation test is carried out based on the braking system model of R-EEV in simulation environment. In simulations, the initial braking speed is set at 80 km/h, and four tire-road conditions were set with PTAC of 0.6, 0.45, 0.2, and 0.1, respectively, as shown in Figure 7. The braking strength increases from 0 to 0.7 quickly to simulate emergency braking conditions, and it stabilizes at 0.5 until the vehicle speed is 0. The estimated value of the PTAC will be output in real time through the TACRM. The effectiveness and accuracy of the TACRM can be seen by comparing with the preset value of PTAC.

Verification of TACRM
The TACRM proposed is verified by simulation test in this section. A co-simulation test is carried out based on the braking system model of R-EEV in simulation environment. In simulations, the initial braking speed is set at 80 km/h, and four tire-road conditions were set with PTAC of 0.6, 0.45, 0.2, and 0.1, respectively, as shown in Figure 7. The braking strength increases from 0 to 0.7 quickly to simulate emergency braking conditions, and it stabilizes at 0.5 until the vehicle speed is 0. The estimated value of the PTAC will be output in real time through the TACRM. The effectiveness and accuracy of the TACRM can be seen by comparing with the preset value of PTAC. The interpolation calculation method is used to estimate the tire-road adhesion condition in the four types of tire-roads, and it can effectively and accurately estimate the PTAC in real time.

OSR Control Method By DFLC
The braking force distribution of the front and rear wheels is different. The membership function of the fuzzy controller and the formulation of the fuzzy rules are also different, which will affect the target slip ratio control result. Considering the difference of braking torque between front and rear wheels, the DFLC for front and rear wheels is designed to ensure OSR tracking control target values effectively, which is different from the fuzzy controllers in [9,13,[21][22][23][24][25]. The requirement slip ratio The interpolation calculation method is used to estimate the tire-road adhesion condition in the four types of tire-roads, and it can effectively and accurately estimate the PTAC in real time.

OSR Control Method By DFLC
The braking force distribution of the front and rear wheels is different. The membership function of the fuzzy controller and the formulation of the fuzzy rules are also different, which will affect the target slip ratio control result. Considering the difference of braking torque between front and rear Energies 2020, 13, 1526 8 of 21 wheels, the DFLC for front and rear wheels is designed to ensure OSR tracking control target values effectively, which is different from the fuzzy controllers in [9,13,[21][22][23][24][25]. The requirement slip ratio (s req (t)) can be obtained from the requirement braking strength calculation and real-time OSR (s osr (t)) can be outputted. The determination process of the real-time OSR control target (s csr (t)) of the control system is as follows: If s req (t) <s osr (t), s csr (t) = s req (t). Otherwise, s csr (t) = s osr (t). As shown in Figure 8, the double fuzzy logic controller is designed in this section. (sreq(t)) can be obtained from the requirement braking strength calculation and real-time OSR (sosr(t)) can be outputted. The determination process of the real-time OSR control target (scsr(t)) of the control system is as follows: If sreq(t) <sosr(t), scsr(t) = sreq(t). Otherwise, scsr(t) = sosr(t). As shown in Figure 8, the double fuzzy logic controller is designed in this section. The Mamdani method is used to perform the fuzzy logic calculation in this paper. Triangular shapes are selected for the membership functions of the inputs and outputs and five fuzzy sets are selected. As shown in Figure 9, in the fuzzy logic controller of front wheels, two inputs are the difference between the real slip ratio and the target slip ratio (ef(t)) and slip ratio change rate (s f(t)). ef(t) ∈ (−0.05, 0.05), s f(t) ∈ (−0.2, 0.2). The input universe is divided into 5 membership functions, which are negative large (NB), negative small (NS), zero (ZO), positive small (PS), positive large (PB). The outputs are the regenerative braking force adjustment value (ΔTf _e(t)) and the hydraulic braking force adjustment value (ΔTf_m (t)). ΔTf_m (t) ∈ (−500, 500), ΔTf _e(t) ∈ (−500, 500). The output universe is divided into 7 membership functions, which are negative large (NB), negative medium (NM), negative small (NS), zero (ZO), positive small (PS), positive medium (PM), and positive large (PB).  The Mamdani method is used to perform the fuzzy logic calculation in this paper. Triangular shapes are selected for the membership functions of the inputs and outputs and five fuzzy sets are selected. As shown in Figure 9, in the fuzzy logic controller of front wheels, two inputs are the difference between the real slip ratio and the target slip ratio (e f (t)) and slip ratio change rate (s f (t)). e f (t) ∈ (−0.05, 0.05), s f (t) ∈ (−0.2, 0.2). The input universe is divided into 5 membership functions, which are negative large (NB), negative small (NS), zero (ZO), positive small (PS), positive large (PB). The outputs are the regenerative braking force adjustment value (∆T f _e (t)) and the hydraulic braking force adjustment value (∆T f_m (t)). ∆T f_m (t) ∈ (−500, 500), ∆T f _e (t) ∈ (−500, 500). The output universe is divided into 7 membership functions, which are negative large (NB), negative medium (NM), negative small (NS), zero (ZO), positive small (PS), positive medium (PM), and positive large (PB).
The fuzzy rules are the key part of the fuzzy controller. Taking the front wheel fuzzy controller as an example, its design principles are as follows: (1) When e f (t)<0: If s f (t)>0, it means that s act (t) has approached s csr (t), and ∆T f should be PS or ZO. If s f (t) = 0, it means that s act (t) is relatively small relative to s csr (t). There is no risk of slipping during braking due to the slip ratio is small, ∆T f should be PS or PM. if s f (t)>0, it means that s act (t) has deviated greatly from s csr (t). In order to eliminate the target deviation as soon as possible, ∆T f is PB. (2) When e f (t) = 0: if s f (t) ≤0, it means that s act (t) is closer to s csr (t), ∆T f should be PS or ZO. If s f (t)>0, it means that s act (t) becomes larger, and the front wheel has the risk of slipping. ∆T f should be NS or NM. (3) When e f (t)>0: If s f (t)>0, it means that s act (t) has deviated seriously from s csr (t), and the wheel has the risk of slipping. ∆T f should be NB or NM to reduce the slip ratio as soon as possible. If s f (t) = 0, it means that s act (t) s relatively large relative to s csr (t). ∆T f should be NM or NS. If s f (t)>0, it means that s act (t) is closer to s csr (t). ∆T f should be NS or ZO. The fuzzy control rules base of front wheel fuzzy controller is shown in Table 3. ef(t) ∈ (−0.05, 0.05), s f(t) ∈ (−0.2, 0.2). The input universe is divided into 5 membership functions, which are negative large (NB), negative small (NS), zero (ZO), positive small (PS), positive large (PB). The outputs are the regenerative braking force adjustment value (ΔTf _e(t)) and the hydraulic braking force adjustment value (ΔTf_m (t)). ΔTf_m (t) ∈ (−500, 500), ΔTf _e(t) ∈ (−500, 500). The output universe is divided into 7 membership functions, which are negative large (NB), negative medium (NM), negative small (NS), zero (ZO), positive small (PS), positive medium (PM), and positive large (PB).   As shown in Figure 10, in the fuzzy controller of rear wheels, two inputs are the difference between the real slip ratio and the target slip ratio (e r (t)) and slip ratio change rate (s r (t)). e r (t) ∈ (−0.05, 0.05), s r (t) ∈ (−0.2, 0.2). The output is the hydraulic braking force adjustment value of rear wheel ∆T r_m (t). The fuzzy rules are the key part of the fuzzy controller. Taking the front wheel fuzzy controller as an example, its design principles are as follows: (1) When ef (t)<0: If sf (t)>0, it means that sact(t) has approached scsr(t), and ∆Tf should be PS or ZO. If sf (t) = 0, it means that sact(t) is relatively small relative to scsr(t). There is no risk of slipping during braking due to the slip ratio is small, ∆Tf should be PS or PM. if sf (t)>0, it means that sact(t) has deviated greatly from scsr(t). In order to eliminate the target deviation as soon as possible, ∆Tf is PB. (2) When ef (t) = 0: if sf (t) ≤0, it means that sact(t) is closer to scsr(t), ∆Tf should be PS or ZO. If sf (t)>0, it means that sact(t) becomes larger, and the front wheel has the risk of slipping. ∆Tf should be NS or NM. (3) When ef (t)>0: If sf (t)>0, it means that sact(t) has deviated seriously from scsr(t), and the wheel has the risk of slipping. ∆Tf should be NB or NM to reduce the slip ratio as soon as possible. If sf (t) = 0, it means that sact(t) s relatively large relative to scsr(t). ∆Tf should be NM or NS. If sf (t)>0, it means that sact(t) is closer to scsr(t). ∆Tf should be NS or ZO. The fuzzy control rules base of front wheel fuzzy controller is shown in Table 3.  Figure 10, in the fuzzy controller of rear wheels, two inputs are the difference between the real slip ratio and the target slip ratio (er(t)) and slip ratio change rate (s ŕ(t)). er(t) ∈ (−0.05, 0.05), s ŕ(t) ∈ (−0.2, 0.2). The output is the hydraulic braking force adjustment value of rear wheel ΔTr_m(t). ΔTr_m (t) ∈ (−500, 500). The selection and definition of the membership function are different with the front wheel fuzzy controller obviously. The possibility of sudden changes in rear wheel braking torque and the risk of overshoot fluctuations should be considered. The principles are roughly the same in the formulation of fuzzy rules. Compared with the front wheel controller, the selection of fuzzy rules should be softer. The fuzzy control rules base of rear wheel fuzzy controller after repeated design adjustments is shown in Table 4.  The fuzzy control rules base of rear wheel fuzzy controller after repeated design adjustments is shown in Table 4. Defuzzification is an important step in fuzzy inference system. There are many methods of defuzzification. The most common used methods are the maximum membership method, the center of gravity method, and the weighted average method. The weighted average method is selected in fuzzy controller for defuzzification, which is widely used in industrial control. The off-line calculation method is adopted to ensure the control accuracy and calculation speed, and the corresponding relationship between the observed value and the actual control value can be calculated using following equation: where X is the final value; n is the number of elements, ω i is the membership.

Simulation Analysis of IRBCS
Based on the simulation model of the regenerative braking system, the braking performance was tested to verify the effectiveness of IRBCS and detect the control effect of the DFLC on the target slip ratio. In simulations, three different roads are set up, namely dry asphalt surface (µ max _ h = 0.8), dry cement surface (µ max _ m = 0.6) and wet cobblestone surface (µ max _ l = 0.35). The initial braking speed is set at 80 km/h. The braking strength increases from 0 to 0.2 and then increases to 0.7 to simulate emergency braking conditions, and it stabilizes at 0.5 until the V vehic = 0. Two evaluation indexes are used to evaluate the energy recovery performance of IRBCS. Evaluation parameter of regenerative braking efficiency (η b ) and energy recovery ratio (C b_cyc ) can be calculated using following equation: where F f_e is the regenerative braking force; F total is the total demand force; E regen_on is the energy consumed with RBS under the single braking test; E regen_off is the energy consumed without RBS under the single braking test. C b_cyc can be expressed as the contribution rate of regenerative braking to the reduction of energy consumption. In order to make the simulation more realistic, the research uses the use of real-like sensor signals affected by errors. During the simulation process, the effects of errors are simulated by adding noise signals.

High Tire-Road Adhesion Condition
As shown in the Figure 11a, IRBCS tries to stop the car and the actual braking strength (Z act ) can reach the required braking strength (Z req ). The battery SoC increases from 0.4 to 0.4053, and C b-cyc is about 61.325%. The speed and slip ratio of front and rear wheel with IRBCS are shown in Figure 11b,c, respectively. We can see that the IRBCS keeps the slip ratio at the optimal value. The actual slip ratio sact(t) responses in time with no overshooting, and the fluctuation is small. It means that the DFLC designed in IRBCS can achieve good follow-up of the OSR. The actual vehicle speed overlaps with the desired vehicle speed. The difference between the vehicle speed and the wheel speed is very small. Figure  11(d) shows that the regenerative braking efficiency and the variations of the braking torque include hydraulic braking torque and regenerative braking torque, the hydraulic braking torque and regenerative braking torque can cooperate well under IRBCS. Initially, the braking torque is all provided by the regenerative braking system and ηb is 1 from 0s to 2.4s to increase regenerative braking energy. When the braking demand is increased with braking torque demand value exceeding the maximum motor torque Tmot_max, the hydraulic braking torque participates in the braking process and ηb∈[0.65, 0.81]; it shows that a large regenerative braking efficiency can still be achieved.

Medium Tire-Road Adhesion Condition
Similar analysis with the previous Section 3.4.1 is as follows: Zact can reach the required braking strength well on medium PTAC, the SoC is increased from 0.4 to 0.4049, and Cb-cyc is about 56.677%. IRBCS keeps the slip ratio at the optimal value, as shown in Figure 12(b) and 12(c). Figure 12(d) shows that variations of the braking torque limiting the follow-up of the target braking intensity when Zreq>0.6 after 3s-3.8s with μmax_h reaching the upper limitation. The target slip ratio is equal to the ground best and provides the maximum braking strength. During the braking process, ηb is in the range of 0.58 to 1; it shows that a large regenerative braking efficiency can still be achieved. We can see that the IRBCS keeps the slip ratio at the optimal value. The actual slip ratio s act (t) responses in time with no overshooting, and the fluctuation is small. It means that the DFLC designed in IRBCS can achieve good follow-up of the OSR. The actual vehicle speed overlaps with the desired vehicle speed. The difference between the vehicle speed and the wheel speed is very small. Figure 11d shows that the regenerative braking efficiency and the variations of the braking torque include hydraulic braking torque and regenerative braking torque, the hydraulic braking torque and regenerative braking torque can cooperate well under IRBCS. Initially, the braking torque is all provided by the regenerative braking system and η b is 1 from 0s to 2.4s to increase regenerative braking energy. When the braking demand is increased with braking torque demand value exceeding the maximum motor torque T mot_max , the hydraulic braking torque participates in the braking process and η b ∈[0.65, 0.81]; it shows that a large regenerative braking efficiency can still be achieved.

Medium Tire-Road Adhesion Condition
Similar analysis with the previous Section 3.4.1 is as follows: Z act can reach the required braking strength well on medium PTAC, the SoC is increased from 0.4 to 0.4049, and C b-cyc is about 56.677%. IRBCS keeps the slip ratio at the optimal value, as shown in Figure 12b,c. Figure 12d shows that variations of the braking torque limiting the follow-up of the target braking intensity when Z req >0.6 after 3-3.8 s with µ max_h reaching the upper limitation. The target slip ratio is equal to the ground best and provides the maximum braking strength. During the braking process, η b is in the range of 0.58 to 1; it shows that a large regenerative braking efficiency can still be achieved.

Low Tire-road Adhesion Condition
As shown in the Figure 13a, Z act can reach the required braking strength well on low adhesion coefficient condition, the SoC is increased from 0.4 to 0.4031, and C b-cyc is about 32.841%. IRBCS keeps the slip ratio at the optimal value, as shown in Figure 13b,c. Figure 13d shows that variations of the braking torque limiting the follow-up of the target braking intensity when Z req >0.35 after 2.2 s with µ max_l reaching the upper limitation. The control target slip ratio is equal to the ground best and provides the maximum braking strength. During the braking process, η b is always 1. Under the Energies 2020, 13, 1526 12 of 21 premise of ensuring the braking efficiency, the braking energy is recovered to the greatest extent, and the control system has excellent performance.

Low Tire-road Adhesion Condition
As shown in the Figure 13(a), Zact can reach the required braking strength well on low adhesion coefficient condition, the SoC is increased from 0.4 to 0.4031, and Cb-cyc is about 32.841%. IRBCS keeps the slip ratio at the optimal value, as shown in Figure 13(b) and 13(c). Figure 13(d) shows that variations of the braking torque limiting the follow-up of the target braking intensity when Zreq>0.35 after 2.2s with μmax_l reaching the upper limitation. The control target slip ratio is equal to the ground best and provides the maximum braking strength. During the braking process, ηb is always 1. Under the premise of ensuring the braking efficiency, the braking energy is recovered to the greatest extent, and the control system has excellent performance. In summary: as the braking control intervenes, the wheel speed decreases smoothly, and the slip ratio maintains the ideal value while in deceleration process. As shown in Figure 11- Figure 13, the wheels actual slip ratio can achieve a good tracking of the target slip ratio on various TAC conditions. It means that the IRBCS can prevent the wheel from locking to ensure braking safety. By using the proposed IRBCS, the accumulated Cb_cyc on the three types of roads are 61.325%, 56.677%, 32.841%, respectively. Thereby, tire-road adhesion is fully utilized and the regenerative braking performance of R-EEV is assured, testifying the effectiveness and adaption of the developed control strategy. In summary: as the braking control intervenes, the wheel speed decreases smoothly, and the slip ratio maintains the ideal value while in deceleration process. As shown in Figures 11-13, the wheels actual slip ratio can achieve a good tracking of the target slip ratio on various TAC conditions. It means Energies 2020, 13, 1526 13 of 21 that the IRBCS can prevent the wheel from locking to ensure braking safety. By using the proposed IRBCS, the accumulated C b_cyc on the three types of roads are 61.325%, 56.677%, 32.841%, respectively. Thereby, tire-road adhesion is fully utilized and the regenerative braking performance of R-EEV is assured, testifying the effectiveness and adaption of the developed control strategy.

Battery Capacity Loss Model (BCLM)
Relevant research shows that the better charging capacity and the better health state of the battery can be obtained by controlling the charging current profiles [26]. In this paper, the Li-ion cylindrical batteries with model M18650 have been selected to study the impact of charging profiles on battery performance. The battery storage systems were subjected to repeat operations of charging and discharging profiles, and the rate and shape of RCC indubitably affected the charging time and the aging rate of a battery. In this section, we reveal the relationship and establish the coupling relationship model between the battery current (I(t)) and the battery capacity loss rate (Q loss,% ). The battery model and the measurements of battery capacity fading are proposed in [26][27][28]. Under a comprehensive driving cycle, the BCLM can be established as follows: where α is the capacity attenuation coefficient; DOD is the depth of discharge; Ah cell is the cumulative capacity of the battery; B 1 and B 2 are fitting parameters, which can be calculated through the least square method; C n is the charging rate; I(t) is the battery current; T cyc is the total cycle time.
It is difficult to evaluate the impact of discharge profiles on battery performance, because most discharge conditions are multifarious based on the load itself [26], and it is necessary to correct the original equation of BCLM with rate and current shape characteristic index. In order to evaluate the capacity fading of the cells during their service life, multi-stage constant current-constant voltage with a negative pulse (MCC-CVNP) is selected to simulate the battery current shape during operation, and the reference performance test was carried out regularly (every 1700 cycles). As shown in Figure 14, six cases have been selected in this paper to investigate the impact of the battery current curve on the lithium-ion battery: Where α is the capacity attenuation coefficient; DOD is the depth of discharge; Ahcell is the cumulative capacity of the battery; B1 and B2 are fitting parameters, which can be calculated through the least square method; Cn is the charging rate; I(t) is the battery current; Tcyc is the total cycle time.
It is difficult to evaluate the impact of discharge profiles on battery performance, because most discharge conditions are multifarious based on the load itself [26], and it is necessary to correct the original equation of BCLM with rate and current shape characteristic index. In order to evaluate the capacity fading of the cells during their service life, multi-stage constant current-constant voltage with a negative pulse (MCC-CVNP) is selected to simulate the battery current shape during operation, and the reference performance test was carried out regularly (every 1700 cycles). As shown in Figure 14, six cases have been selected in this paper to investigate the impact of the battery current curve on the lithium-ion battery: Case 1 (CC-CVNP): constant current discharging at current rate Id for Td and constant current charging at current rate Ich for Tch, followed by rest time Ts were applied.
Cases 2-6 (CC-CVNP): initial constant current discharging at current rate I ⃰ d for Td and initial constant current charging at current rate I ⃰ ch for Tch, followed by rest time Ts were applied. Each pulse decreases in sequence of ∆Id. Under the same discharge rate Cn, the parameters of six cases are designed as follows: Where ∆Id is the decrement of discharging current; ∆Ich is the decrement of charging current. As shown in Figure 14 (b), each pulse decreases in sequence of ΔId. The reduced current value Case 1 (CC-CVNP): constant current discharging at current rate I d for T d and constant current charging at current rate I ch for T ch , followed by rest time T s were applied. Cases 2-6 (CC-CVNP): initial constant current discharging at current rate I * d for T d and initial constant current charging at current rate I * ch for T ch , followed by rest time T s were applied. Each pulse decreases in sequence of ∆I d .
Under the same discharge rate C n , the parameters of six cases are designed as follows: where ∆I d is the decrement of discharging current; ∆I ch is the decrement of charging current. As shown in Figure 14b, each pulse decreases in sequence of ∆I d . The reduced current value becomes equal difference series size relationship, which increases in sequence (∆I d , 2∆I d , 3∆I d . . . ). The relationship between ∆I d , ∆I ch , I * d , and I * ch needs to satisfy the Equation (20) to ensure that maintain the same charging rate under the shape of pulse current.
S c 2 is defined to characterize the fluctuation of current shape under cycle condition and is calculated as follows: The related research results reveal that the dynamic discharging-charging profile has a significant impact on reducing the capacity fade of the lithium-ion battery cells compared with the static profile [26]. Therefore, case 1 is selected as reference scheme because of their minimum S c 2 in this paper, and k s2 is defined as follows: To characterize the fluctuation of current shape under cycle condition, S wc 2 (t) can be calculated as follows: Commercial lithium-ion M18650 cells has been used to study the battery capacity loss rate Q loss,% . Table 5 shows the electrical parameters of battery cells. Six cases of the applied charging process are shown in Table 6. Table 6. Six cases of the applied discharging process of C n = 0.5C.  The same test procedure can be performed in the case of C n = 1.0C, C n = 2.0C and C n = 3.0C. The relationship of correction factor k s2 (t) and S wc 2 (t) can be computed and concluded as shown in Figure 15. 1.0 0.5 0.5 21 The same test procedure can be performed in the case of Cn = 1.0C, Cn = 2.0C and Cn = 3.0C. The relationship of correction factor ks2(t) and Swc 2 (t) can be computed and concluded as shown in Figure  15. According to equation (21) and equation (22), Correction factor ks2 can be calculated as follows: According to Equation (21) and Equation (22), Correction factor k s2 can be calculated as follows: By fitting the experimental data, the following can be acquired: C 1 = 1.0419 and C 2 = 0.0363. Equation (23) can be calculated as follows: Finally, according to the above analysis of the relationship between k s2 and S wc 2 , the BCLM with current curve shape considerations under the cycle condition is obtained as follows:

Regenerative Braking Work-Point Switching
The goal of R-EEV's regenerative braking control strategy is to ensure comprehensive regenerative braking performance and the OSR-based RRCCS is proposed. The shape and rate of the RCC characterize the details of the battery being charged and discharged during vehicle operation, which can be seen in the battery SoC. In order to control the shape and rate of the battery RCC and reduce the adverse effect of RCC on battery life. The operating point switching mechanism in RRBCS is to conservatively adjust the operating point by referring to the battery SoC characteristics. Through the "conservative braking mode", the actual SoC will be close to the pre-set value to make full use of the power of the battery pack, while protecting the battery health. As shown in Figure 16, the switching of the braking working point is adjusted in real time according to the vehicle state. W C1 is switched to W C2 , and W C4 is switched to W C3 .
The reference index for braking work-point switching logic is calculated as follows: where k SoC is the SoC change rate; k cr is the change rate discrimination value; SoC (k) and SoC 0 (k) are the battery real time SoC and the SoC controls target value at the k instant; SoC (k + 1)) and SoC 0 (k+1) are the battery real time SoC and the SoC controls target value at the k+1 instant; k 0 and k 0_CD are the SoC control target change rate and preset value in charge depletion, and the k 0_CS is the SoC preset value in charge sustaining with the value assigned is 0. ζ 0 is a correction coefficient and the value assigned to the study is 0.5.
regenerative braking performance and the OSR-based RRCCS is proposed. The shape and rate of the RCC characterize the details of the battery being charged and discharged during vehicle operation, which can be seen in the battery SoC. In order to control the shape and rate of the battery RCC and reduce the adverse effect of RCC on battery life. The operating point switching mechanism in RRBCS is to conservatively adjust the operating point by referring to the battery SoC characteristics. Through the "conservative braking mode", the actual SoC will be close to the pre-set value to make full use of the power of the battery pack, while protecting the battery health. As shown in Figure 16, the switching of the braking working point is adjusted in real time according to the vehicle state. WC1 is switched to WC2, and WC4 is switched to WC3. The reference index for braking work-point switching logic is calculated as follows: Where kSoC is the SoC change rate; kcr is the change rate discrimination value; SoC (k) and SoC0 (k) are the battery real time SoC and the SoC controls target value at the k instant; SoC (k + 1)) and SoC0 (k+1) are the battery real time SoC and the SoC controls target value at the k+1 instant; k0 and k0_CD are the SoC control target change rate and preset value in charge depletion, and the k0_CS is the SoC preset value in charge sustaining with the value assigned is 0. ζ0 is a correction coefficient and the value assigned to the study is 0.5.
The details are explained as follows: (1) In charge depletion (CD), when the SoC＞0, if the ksoc is positive and greater than the discrimination value kcr, it means that the SoC has a greater increasing trend and the possibility of deviation from the SoC0, "conservative braking mode" should be adopted.
(2) If kSoC is negative and |kSoC| <|kcr|, it means that the SoC has a trend of approaching the SoC0, but the "fallback" is slow. At this time, a "conservative braking mode" should be adopted. (3) In charge sustaining (CS), when SoC> 0, it means that the difference between SoC and SoC0 is too large. There The details are explained as follows: (1) In charge depletion (CD), when the SoC>0, if the k soc is positive and greater than the discrimination value k cr , it means that the SoC has a greater increasing trend and the possibility of deviation from the SoC 0 , "conservative braking mode" should be adopted.
(2) If k SoC is negative and |k SoC | <|k cr |, it means that the SoC has a trend of approaching the SoC 0 , but the "fallback" is slow. At this time, a "conservative braking mode" should be adopted. (3) In charge sustaining (CS), when SoC> 0, it means that the difference between SoC and SoC 0 is too large. There is a risk of increasing the deviation of the recovered braking energy at this time; thus, the "conservative braking method" should be adopted.

Comparison and Analysis of Three Control Strategies
The world light vehicle test procedure (WLTP) is close to actual driving behavior. In this paper, the WLTP driving cycle is adopted to carry out the tests for studying the comprehensive braking performance of R-EEV. Under the WLTP driving cycle, the vehicle braking intensity Z req ∈ [0, 0.15], it means that the regenerative braking work-point is W C1 or W C2 . In addition, no braking work condition is defined as the work-point 0 (W C0 ).
Control strategy 1 (the best-braking-performance strategy) is adopted to guarantee the best braking performance. The regenerative braking work-point is always W C1 in Control strategy 1. Control strategy 2 (the maximum-regeneration-efficiency strategy) is the IRBCS proposed mentioned before. The regenerative braking work-point is W C2 in Control strategy 2 to maximize the use of the regenerative braking torque and achieve the maximum regeneration efficiency theoretically. Control strategy 3 (the switchable-work-point strategy) is RRBCS through adding "conservative condition" based on IRBCS. The regenerative braking work-point can be switched from W C2 to W C1 in Control strategy 3 to balance the recovery energy and battery capacity loss rate. Comparative simulations are carried out among three control strategies, as shown in Figure 17. It should be noted that the braking energy recovery will be limited by the motor maximum torque. The motor torque is assumed to provide sufficient required braking torque, and the braking energy is recovered at all braking strengths in theory. This braking recovery is defined as the theoretical maximum regenerative braking to illustrate the contribution of the developed strategy to energy recovery.
As depicted in Figure 17a, the vehicle speed can be in good agreement with the target speed under the control strategy 1 and the control strategy 2. It indicated that the IRBCS proposed in this paper can effectively ensure braking efficiency and braking safety. Figure 17b shows that the regenerative braking efficiency of strategy 1 has decreased from 100% to about 56.7%, and the regenerative braking efficiency of strategy 2 has always been 100%. As shown in Figure 17c,d, after a WLTP driving cycle, battery SoC drops from the initial value of 0.8 to 0.671 and 0.683, respectively. The energy recovery ratio C b-cyc of strategy 1 is above 13.6%. In particular, the energy recovery ratio C b-cyc even reaches 24.8% for strategy 2, which is a relatively high level. In addition, there are about 37% of the regenerative braking work-points switched from W c2 to W c1 in strategy 3 during the driving cycle; it means that some of the braking recovery energy will not be fully utilized, as shown in Figure 18.
before. The regenerative braking work-point is WC2 in Control strategy 2 to maximize the use of the regenerative braking torque and achieve the maximum regeneration efficiency theoretically. Control strategy 3 (the switchable-work-point strategy) is RRBCS through adding "conservative condition" based on IRBCS. The regenerative braking work-point can be switched from WC2 to WC1 in Control strategy 3 to balance the recovery energy and battery capacity loss rate. Comparative simulations are carried out among three control strategies, as shown in Figure 17. It should be noted that the braking energy recovery will be limited by the motor maximum torque. The motor torque is assumed to provide sufficient required braking torque, and the braking energy is recovered at all braking strengths in theory. This braking recovery is defined as the theoretical maximum regenerative braking to illustrate the contribution of the developed strategy to energy recovery. As depicted in Figure 17 (a), the vehicle speed can be in good agreement with the target speed under the control strategy 1 and the control strategy 2. It indicated that the IRBCS proposed in this paper can effectively ensure braking efficiency and braking safety. Figure 17(b) shows that the regenerative braking efficiency of strategy 1 has decreased from 100% to about 56.7%, and the regenerative braking efficiency of strategy 2 has always been 100%. As shown in Figure 17(c), (d), after a WLTP driving cycle, battery SoC drops from the initial value of 0.8 to 0.671 and 0.683, respectively. The energy recovery ratio Cb-cyc of strategy 1 is above 13.6%. In particular, the energy recovery ratio Cb-cyc even reaches 24.8% for strategy 2, which is a relatively high level. In addition, there are about 37% of the regenerative braking work-points switched from Wc2 to Wc1 in strategy 3 during the driving cycle; it means that some of the braking recovery energy will not be fully utilized, as shown in Figure 18. The battery SoC and current curve in the three braking modes based on the WLTP cycle is shown in Figure 19 and Figure 20. Through the SoC comparison of the three strategies, it is not hard to discover that the IRBCS (strategy 2) and RRBCS (strategy 3) possess a more superior economy and reduces electricity consumption as much as possible. The battery SoC and current curve in the three braking modes based on the WLTP cycle is shown in Figures 19 and 20. recovery ratio Cb-cyc even reaches 24.8% for strategy 2, which is a relatively high level. In addition, there are about 37% of the regenerative braking work-points switched from Wc2 to Wc1 in strategy 3 during the driving cycle; it means that some of the braking recovery energy will not be fully utilized, as shown in Figure 18. The battery SoC and current curve in the three braking modes based on the WLTP cycle is shown in Figure 19 and Figure 20. Through the SoC comparison of the three strategies, it is not hard to discover that the IRBCS (strategy 2) and RRBCS (strategy 3) possess a more superior economy and reduces electricity consumption as much as possible.  Through the SoC comparison of the three strategies, it is not hard to discover that the IRBCS (strategy 2) and RRBCS (strategy 3) possess a more superior economy and reduces electricity consumption as much as possible. Compared with Control strategy 2, the change of energy recovery ratio ∆Cb_cyc and the change of battery capacity loss rate ∆Qloss,% are calculated as follows: Through the SoC comparison of the three strategies, it is not hard to discover that the IRBCS (strategy 2) and RRBCS (strategy 3) possess a more superior economy and reduces electricity consumption as much as possible.
Compared with Control strategy 2, the change of energy recovery ratio ∆C b_cyc and the change of battery capacity loss rate ∆Q loss,% are calculated as follows: where C i b_cyc is the energy recovery ratio C b_cyc of control strategy i (I = 1, 2, 3); Q i loss,% is the battery capacity loss rate Q loss,% of control strategy i (i = 1, 2, 3).
The comparison results of regenerative braking performance among the three strategies are listed in Table 7. As depicted in Figure 21, regenerative braking efficiency η b under three control strategies are 49.59%, 84.43% and 76.51%, respectively. The energy recovery ratio C b_cyc under three control strategies are 8.5%, 24.8% and 22.1%, respectively.
The C b_cyc slightly reduced by 16.6% and the Q loss,% reduced by 79.2%. It infers that the strategy 3 can improve battery life to a certain extent with less loss of braking energy recovery. This means that strategy 3 balances the regenerative braking energy and the battery capacity attenuation as much as possible. The mitigation effect of the strategy 3 on the battery capacity loss should be more clearly presented, especially during the entire battery life. Therefore, the results of Q loss,% after 1000-10,000 times under WLTP driving cycles are shown in Table 8.  Figure 21, regenerative braking efficiency ηb under three control strategies are 49.59%, 84.43% and 76.51%, respectively. The energy recovery ratio Cb_cyc under three control strategies are 8.5%, 24.8% and 22.1%, respectively. The Cb_cyc slightly reduced by 16.6% and the Qloss,% reduced by 79.2%. It infers that the strategy 3 can improve battery life to a certain extent with less loss of braking energy recovery. This means that strategy 3 balances the regenerative braking energy and the battery capacity attenuation as much as possible. The mitigation effect of the strategy 3 on the battery capacity loss should be more clearly presented, especially during the entire battery life. Therefore, the results of Qloss,% after 1000-10000 times under WLTP driving cycles are shown in Table 8.  Table 8, the battery Qloss,% gradually increases with the number of WLTP driving cycles. After 10,000 WLTP driving cycles, Qloss,% under control strategy 1 is 19.64%. By comparison, Qloss,% under control strategy 2 and control strategy 3 are 36.52% and 24.97%, respectively. It can indicate that the proposed control strategy 3 avoids the negative impact of the large braking feedback of strategy 2 on battery life and inherits the battery friendliness of control strategy 1. Therefore, the  As shown in Table 8, the battery Q loss,% gradually increases with the number of WLTP driving cycles. After 10,000 WLTP driving cycles, Q loss,% under control strategy 1 is 19.64%. By comparison, Q loss,% under control strategy 2 and control strategy 3 are 36.52% and 24.97%, respectively. It can indicate that the proposed control strategy 3 avoids the negative impact of the large braking feedback of strategy 2 on battery life and inherits the battery friendliness of control strategy 1. Therefore, the proposed control strategy 3, namely the switchable-work-point strategy, has better comprehensive braking performance, which can provide a feasible reference for regenerative braking control strategy development.

Conclusions
In this paper, a RRBCS with the rate and shape of battery current consideration is proposed to balance the comprehensive regenerative braking performance and battery service life. The test results validate the effectiveness and feasibility of the proposed control strategy. Firstly, The IRBCS based on OSR is introduced, and the aim anti-lock control braking torque based on OSR is provided to improve braking performance while maximizing braking energy recovery. It has been shown that, as expected, IRBCS can achieve good follow-up control of the target slip ratio by DFLC. During the control process, the actual slip ratio responds quickly without overshooting. At the same time, under the three types of tire-road adhesion conditions, IRBCS can maximize the braking energy regenerative efficiency while guaranteeing good braking efficiency with the energy regenerative ratio is 82.46%, 86.43% and 100%, respectively. Moreover, three different braking energy recovery strategies, namely the best-braking-performance strategy (Control strategy 1), the maximum-regeneration-efficiency strategy (Control strategy 2, IRBCS), and the switchable-work-point strategy (Control strategy 3, RRCBS) are implemented and analyzed during WLTP driving cycle. The experimental results indicate the effectiveness and adaptability of the RRBCS. About 37% of work-points are switched from W c2 to W C1 under strategy 3, with regenerative brake efficiency slightly reduced by 16.6% and the decline in the battery's capacity reduced by 79.2%, it infers that the RRBCS can reduce Q loss,% to extend its service life significantly with less loss of braking energy recovery. The research methods and conclusions proposed in this paper can provide a powerful reference for regenerative braking control strategy design of R-EEVs.
Author Contributions: Y.L. managed the project; Y.F. conceptualized the scheme; H.L. conceived the control method, completed the modeling, performed the simulation experiments, and finished the manuscript; X.L. collected the data and reviewed the paper. All authors have read and agreed to the published version of the manuscript.