Abstract
To improve the energy recovery rate and increase the driving range of load-isolated electric buses, a composite brake control strategy based on fuzzy logic optimization is proposed. Considering the influence of the battery SOC value and the braking intensity on the braking ratio, a fuzzy controller is designed for small and medium braking strengths to optimize the control strategy. The vehicle simulation model of the load-isolated electric bus was built with the AVL CRUISE software platform. The AVL CRUISE-Simulink cosimulation was carried out under the original control strategy and the optimized control strategy. The simulation results show that, under the premise of ensuring the stability of braking, the driving range of the vehicle with the optimized control strategy is increased by 7.74% and the energy recovery rate is increased by 11.05%.
1. Introduction
Regenerative braking, also known as braking energy recovery, is a new technology that is applied to new energy vehicles [1, 2]. Through the regenerative braking technology, part of the kinetic energy of an electric vehicle can be converted into electric energy. The electric energy is stored in the battery to improve energy utilization and increase the driving range [3, 4]. The load-isolated electric bus is an electric bus that isolates the engine from the load, ensuring that the engine is operating in the optimum efficiency zone. The regenerative braking system and the original hydraulic brake system of load-isolated electric buses constitute a composite braking system. Therefore, it is of research significance to study the composite braking energy recovery of electric buses.
Based on the fuzzy control theory, Li et al. [5] designed a regenerative braking control strategy for the power system of a fuel cell hybrid vehicle. The simulation verified the strategy to recover the braking energy and maintain the battery SOC (state of charge). To solve the hysteresis problem of hydraulic brake systems [6], a brake feel consistency control strategy was designed. A series control strategy was proposed for regenerative braking and pneumatic braking systems, and a road test was conducted [7]. A brake system with only motor/generator braking was proposed, and the effectiveness of the braking system was verified by simulations and experiments [8]. Using the fuzzy control method, Amirabadi and Farhangi [9] designed a control strategy for one-way and bidirectional converters in fuel cell hybrid vehicles. Considering the acceleration, jerk, and inclination of the vehicle, a regenerative braking fuzzy logic model was proposed [10]. For the regenerative braking system of fuel cell hybrid vehicles, Li et al. [11] proposed a fuzzy logic controller to improve the fuel economy of the vehicle. Considering the effects of the battery SOC, braking strength, and motor speed, Xiao et al. [12] proposed a regenerative braking strategy based on fuzzy control logic. Yuan et al. [13] proposed a new type of regenerative electrohydraulic composite braking system and proved its practicability through simulations and experiments.
The economics of the load-isolated electric bus are very good. If the regenerative braking technology is applied and the composite braking method is combined with the original mechanical brake, the energy-saving potential of the electric bus will be largely realized [14, 15]. However, only a few studies have investigated the application of regenerative braking technology to load-isolated electric vehicles. This paper draws on the above research methods, combined with the actual situation of a load-isolated electric bus, and applies the regenerative braking method to the vehicles. This is a method improvement.
At the same time, the current research on control strategy optimization is mainly divided into two aspects: one is optimized by the Isight software and the other is optimized by the fuzzy toolbox of the MATLAB software. Using the MATLAB software optimization is to use the fuzzy toolbox and to set the fuzzy control rules for optimization. This paper uses MATLAB to optimize the control strategy and to propose a new fuzzy control rule and a fuzzy controller construction method.
Combined with the structure of the load-isolated electric bus, this paper proposes a series control strategy for the composite brake system. The fuzzy controller is designed for small braking and medium braking strengths to optimize the control strategy. Finally, the AVL CRUISE software and the Simulink software are used to carry out a cosimulation analysis on the proposed control strategy and optimization method to improve energy utilization and increase the driving range.
2. Load-Isolated Electric Bus Composite Brake Control Strategy
The experimental model in this paper is a load-isolated electric bus with rear-wheel drive. The load-isolated electric bus is a new type of electric car proposed by Qingdao University. It uses proprietary topology circuits and control strategies to protect the engine and energy storage batteries from changes in the vehicle load. At the same time, the engine is guaranteed to work in the best economic conditions, to improve the energy utilization rate of the vehicle and to increase the driving range. The load-isolated electric bus power source system consists of two roads: one is the energy storage battery-drive motor and the other is the engine-generator-power battery-drive battery. The energy storage battery is used as the main energy supply module to keep the bus run normally. When the energy storage battery SOC is reduced to the lowest power supply limit, the power battery system as an auxiliary energy module continues to provide energy to the bus. The design of the two batteries in the power system avoids the problem of overcharging and overdischarging of a single battery and ensures the running state of the electric bus.
A regenerative braking system is introduced in the brake system of the load-isolated electric bus, and a composite braking system is formed with the original mechanical friction braking system. The purpose of the regenerative braking system is to recover the braking energy as much as possible while improving the braking stability of the electric bus and to improve the energy utilization rate of the electric vehicle [16, 17]. The regenerative braking method only acts on the drive wheels, which in this paper, are the rear wheels of the electric bus.
2.1. Composite Brake System Structure
The brake system of the load-isolated electric bus includes a mechanical friction brake system and a motor regenerative brake system called the composite brake system. It is a brake energy recovery system. According to the coupling relationship between the feedback braking force and the friction braking force, the braking energy recovery system can be divided into a superimposed type (parallel type) and a coordinated type (series type).
The parallel brake energy recovery system superimposes the motor feedback braking force directly on the original friction braking force. It does not adjust the original friction braking force. The advantage is that the implementation is convenient, while the disadvantage is that the feedback efficiency is low and the braking feeling is poor. The principle of the series brake energy recovery system is to preferentially use the feedback braking force and then adjust the hydraulic braking force (or friction braking force) accordingly. It aligns the sum of the two braking forces with the total braking demand. The advantage is that the feedback efficiency is higher and the braking feeling is better. The disadvantage is that the traditional hydraulic brake system must be modified and the implementation is complicated. After comprehensive consideration, this paper uses a series braking energy recovery system, as shown in Figure 1, and carries out simulation modeling and analysis.

2.2. Braking Force Distribution Curve and Braking Force Distribution Method
In the composite brake system of the load-isolated electric bus, there are only a mechanical friction braking force on the front wheel and a mechanical friction braking force and a motor regenerative braking force on the rear wheel. The distribution relationship of each braking force is shown in Figure 2.

According to the ECE brake regulations proposed by the United Nations Economic Commission for Europe, this paper proposes a braking force distribution curve for the brake system of the load-isolated electric bus. As shown in Figure 3, the solid red line is the braking force distribution curve of the front and rear wheels of the load-isolated electric bus.

The dashed colored lines are r-line groups. The r-line group is a curve of the front and rear ground braking force when the front wheel is not locked and the rear wheel is locked.
The dashed green lines are the brake strength curves with different adhesion coefficients.
The purple line group is the f-line group. The f-line group is the relationship between the front and the rear ground braking force when the rear wheel is not locked and the front wheel is locked on various φ values.
The black line is the standard line of the ECE regulations.
The braking force distribution curve equation of the rear drive iswhere is the vehicle’s centroid height, G is the gravity of the vehicle, L is the wheelbase of the vehicle, Fxb1 is the ground braking force of the front wheel, Fxb2 is the ground braking force of the rear wheel, and b is the distance from the vehicle’s center of mass to the rear axle.
The braking force distribution method of the front and rear wheels during the braking process of the load-isolated electric bus is as follows:(1)When the braking strength of the electric bus z ≤ 0.1, the braking strength is small. The braking force distribution curve of the small braking strength is shown in section OA in Figure 3. At this braking strength, the braking code does not restrict the braking force distribution of the front and rear wheels. To maximize the recovery of braking energy, only regenerative braking is used. Therefore, only the motor braking force of the rear wheel works, and the mechanical friction braking force of the front wheel and the rear wheel does not work. The braking strength (z) is defined as follows, where is the deceleration of the vehicle and is gravitational acceleration: At small braking strengths, the braking force distribution of the front and rear wheels of the electric bus is as follows: The specific distribution of the braking force between the front and the rear wheels is as follows: In equations (3) and (4), Fbf is the braking force of the front wheels, Fbr is the braking force of rear wheels, m is the quality of the vehicle, is the acceleration of gravity, z is the braking strength, Fbf-fric is the mechanical friction braking force provided by the front wheels, Fbr-fric is the mechanical friction braking force provided by the rear wheels, and Freg is the regenerative braking force provided by the motor.(2)When the braking strength of the electric bus 0.1 < z ≤ 0.234, the braking strength is considered to be in the middle range. The braking force distribution curve of the middle braking strength is shown in section AB in Figure 3. In this braking strength, the braking force of the front and rear wheels of the electric bus is distributed according to the curve of r = 0.234. The r-line group is a curve of the front and rear ground braking force when the front wheel is not locked and the rear wheel is locked. At the middle braking strength, the braking force distribution of the front and rear wheels of the electric bus is as follows: The specific distribution of the braking force between the front and the rear wheels is as follows: where is the total demand braking force of the vehicle, φ is the adhesion coefficient, and a is the distance from the vehicle’s center of mass to the front axle.(3)When the braking strength of the electric bus 0.234 < z ≤ 0.61, the braking strength is large. The braking force distribution curve for the large braking strength is shown in section BC in Figure 3. At this braking strength, the braking force of the front and rear wheels of the electric bus is distributed according to an ideal braking force distribution curve (I curve). At large braking strengths, the braking force distribution of the front and rear wheels of the electric bus is as follows: The specific distribution of the braking force between the front and the rear wheels is as follows: In equation (7), Fz1 is the normal reaction force of the ground to the front wheel and Fz2 is the normal reaction force of the ground to the rear wheel.(4)When the braking strength of the electric bus z > 0.61, the braking strength is in the emergency range. To ensure the braking stability of the electric bus, regenerative braking is not used. The braking force requirement of the electric bus is satisfied by the mechanical friction braking force of the front and rear wheel brakes. At the emergency braking strength, the braking force distribution of the front and rear wheels of the electric bus is as follows: The specific distribution of the braking force between the front and the rear wheels is as follows:
2.3. Compound Braking Control Strategy Logic
At present, the widely used braking force distribution control strategies are of series and parallel types. The series control strategy can be divided into the energy maximization control strategy, the sensory optimal control strategy, and the fuzzy logic control strategy. The parallel control strategy is divided into a simple parallel control strategy and an asynchronous parallel control strategy [18–20]. The advantages of the parallel control strategy are its simple structure, low cost, and good braking stability. Therefore, it is convenient to promote it in a wide range. The disadvantage is that the brake energy recovery rate is relatively low [21–23]. The energy recovery rate can be improved by optimizing the control strategies.
The serial control strategy is adopted in this paper. According to the proposed braking force distribution curve and the braking force distribution mode under each braking strength, the control strategy logic diagram of the compound braking system is shown in Figure 4. In the figure, Y represents YES and N represents NO in the control strategy proposed in this paper, except for the case in which the braking strength z takes into account “equal.” The battery SOC signal (including the Csoc and the Dsoc), the vehicle speed signal V, and the maximum regenerative braking force signal Fm-act do not take into account “equal.”

During the braking of the load-isolated electric bus, the sensor detects the brake pedal opening signal. After receiving the signal, the vehicle controller VCU determines the braking intention of the electric bus driver and the braking strength of the electric bus braking system. As shown in Figure 4, the magnitude of the braking force is determined by the control strategy logic.
When the braking strength z ≤ 0.1, the braking strength is small. For the purpose of protecting the battery, the energy storage battery SOC (Csoc) and the power battery SOC (Dsoc) signals are comprehensively considered. Regenerative braking will not be used if both signals are greater than 80%. Mechanical brakes are used on the front and rear wheels. It is judged whether or not there is a signal of less than 80% among the two signals. If there is, then the velocity signal needs to be judged. It is not possible to use regenerative braking when the velocity signal V is lower than the minimum velocity Vmin. Vmin is the minimum velocity to perform regenerative braking. If at least one battery SOC set is less than 80% and the velocity signal V is also greater than the minimum velocity Vmin, the actual maximum regenerative braking force Fm-act of the motor is calculated by the control strategy. If Fm-act is greater than the rear wheel demand braking force Fbr, then the front wheel demand braking force is provided by the front wheel mechanical friction braking force, and the rear wheel demand braking force is provided by the motor regenerative braking force. The rear wheel mechanical friction braking force does not work. If Fm-act is smaller than the rear wheel demand braking force Fbr, then the front wheel demand braking force is provided by the front wheel mechanical friction braking force. The rear wheel demand braking force is provided by the rear wheel mechanical friction braking force and the motor regenerative braking force.
When the braking strength 0.1 < z ≤ 0.61, the braking strength belongs to the middle and the large range. When 0.1 < z ≤ 0.234, the braking strength belongs to the middle range. When 0.234 < z ≤ 0.61, the braking strength is large. After receiving the SOC signal detected by the battery SOC sensor, the control strategy makes a logical determination. If both Csoc and Dsoc are greater than 80%, all mechanical friction braking is used, and regenerative braking is not used. If at least one SOC signal is less than 80%, the control strategy further determines whether the vehicle speed signal V is less than Vmin. If V is less than Vmin, all mechanical friction braking is used.
If V is greater than Vmin, the control strategy calculates the total demand braking force Fz and Fm-act. The front wheel demand braking force is provided by the front wheel mechanical friction braking force. The rear wheel demand braking force is determined by the size of Fm-act and Fbr. If Fm-act > Fbr, the rear wheel demand braking force is provided by the motor regenerative braking force. If Fm-act < Fbr, the rear wheel demand braking force is provided by the rear wheel mechanical friction braking force and the motor regenerative braking force.
When the braking strength z > 0.61, the braking strength is in the emergency range. To ensure the braking stability of the electric bus, regenerative braking is not used. The braking force is entirely provided by the mechanical friction braking force.
This paper models the control strategy of the composite brake system based on the proposed control logic.
3. Fuzzy Controller Design
In classical set theory, the membership of an element to a set can only take values of 0 and 1 to represent “belonging” and “not belonging.” The fuzzy set theory shows that the membership degree can be any value in the interval [0, 1]. In the case of not accurately understanding the controlled object model, the fuzzy control theory optimizes the system through fuzzy rules. Because it does not need to know the exact mathematical model to get the output that meets a designer’s vision, it has been widely used in the field of system control [24–27]. The fuzzy controller is the key technology in fuzzy control theory. It includes membership functions, fuzzy control rules, and influence surfaces. Designers make fuzzy rules according to the system environment and blur the input signal. Then, it is optimized by fuzzy rules. Finally, the output signal is obtained by defuzzifying [28–30].
In this paper, two different fuzzy controllers are designed for small and middle braking strengths. The small braking strength controller is a fuzzy controller with dual inputs and a single output. The middle braking strength controller is a three-input and dual-output fuzzy controller.
3.1. Fuzzy Controller for Optimizing Small Braking Strength Control Strategy
3.1.1. Input Variable and Output Variable
In a fuzzy controller with a small braking strength, there are two inputs and one output. The two inputs are the energy storage battery SOC (Csoc) and the power battery SOC (Dsoc). The one output is the regenerative braking ratio of the motor (REG), as shown in Figure 5.

3.1.2. Fuzzifying
The fuzzy subset of the input signal energy storage battery SOC (Csoc) is {L, M, H}, and the fuzzy domain is [0, 1]. The fuzzy subset of the power battery SOC (Dsoc) is {L, M, H}, and the fuzzy domain is [0, 1]. According to the two input quantities and the one output quantity in the fuzzy controller, the corresponding membership functions are designed in the fuzzy toolbox in MATLAB. The input membership function is shown in Figure 6. Figures 6(a) and 6(b) show the membership functions of the input signals Csoc and Dsoc, respectively.

(a)

(b)
3.1.3. Fuzzy Control Rules
Fuzzy control rules are the core content of fuzzy controllers and are important factors affecting output signals. The fuzzy control rules established are shown in Table 1.
The fuzzy inference surface of the input and output is shown in Figure 7. The inference surface is observed by the surface observer. It visually reflects the three-dimensional relationship between input and output.

Figure 7 describes the effect of Csoc and Dsoc on the ratio of regenerative braking from the perspective of three-dimensional relationships.
3.1.4. Defuzzifying
The fuzzy subset of the output signal motor regenerative braking ratio REG is {L, M, H}, and the fuzzy domain is [0, 1]. Figure 8 shows a membership function of the output signal REG.

3.2. Fuzzy Controller for Optimizing Middle Braking Strength Control Strategy
3.2.1. Input Variable and Output Variable
In the middle braking strength fuzzy controller, there are three inputs and two outputs. The three inputs are the energy storage battery SOC (Csoc), the power battery SOC (Dsoc), and the braking intensity (z). The two outputs are the motor regenerative braking ratio (REG) and the mechanical friction braking ratio (FRIC), as shown in Figure 9.

3.2.2. Fuzzifying
The fuzzy subset of the input signal energy storage battery SOC (Csoc) is {L, M, H}, and the fuzzy domain is [0, 1]. The fuzzy subset of the power battery SOC (Dsoc) is {L, M, H}, and the fuzzy domain is [0, 1]. The braking intensity z fuzzy subset is {L, M, H}, and the fuzzy domain is [0, 1]. The input membership functions are shown in Figure 10. Figures 10(a)–10(c) show the membership functions of the input signals Csoc, Dsoc, and z, respectively.

(a)

(b)

(c)
3.2.3. Fuzzy Control Rules
The effects of Csoc, Dsoc, and z on the motor regenerative braking ratio REG are shown in Table 2. The effects of Csoc, Dsoc, and z on the mechanical friction braking ratio FRIC are shown in Table 3.
Figure 11(a) shows a fuzzy inference surface diagram of REG as a function of Csoc and Dsoc for medium braking intensity. Figure 11(b) shows a curved surface diagram showing the influence of Csoc and z on REG in the middle braking strength. The effects of Csoc and z on REG are described from the perspective of three-dimensional relationships. The larger the Csoc and z, the smaller the motor regenerative braking ratio REG. This is also in line with the development of control strategies. It can be seen from Figure 11(b) that as the braking strength increases, the regenerative braking ratio of the motor gradually decreases, which is in accordance with the set rules of the control strategy.

(a)

(b)
3.2.4. Defuzzifying
The fuzzy subset of the output signal motor regenerative braking ratio REG is {L, M, H}, and the fuzzy domain is [0, 1]. The fuzzy subset of the mechanical friction braking ratio FRIC is {L, M, H}, and the fuzzy domain is [0, 1]. The output membership functions shown in Figures 12(a) and 12(b) are the membership functions of the output signals REG and FRIC, respectively.

(a)

(b)
4. Control Strategy Fuzzy Logic Optimization Modeling
Based on the original control strategy, the fuzzy controller module in the MATLAB fuzzy toolbox is added. The fuzzy controller designed in Section 2 is added to the small and middle braking strength strategies to optimize the signal of the regenerative braking ratio and the friction braking ratio. Through the optimization of the fuzzy controller, the regenerative braking ratio and the energy utilization rate are increased [31, 32]. Figure 13 shows the model principle of the fuzzy controller. When the braking strength is small, the Csoc signal and the Dsoc signal are input into the small braking strength fuzzy controller. After the calculation of the fuzzy rules and the fuzzy reasoning, the REG signal is outputted. When the braking strength is in the middle range, the Csoc signal, Dsoc signal, and z signal are input into the middle braking strength fuzzy controller. Then, the output signals REG and FRIC are obtained by fuzzy reasoning.

The motor speed and torque characteristics are shown in Figure 14. When the battery is being charged, the relationship between the SOC and the voltage is shown in Figure 15.


5. Simulation Analysis of Optimal Control Strategy
5.1. Braking Stability Simulation Analysis
To verify the effectiveness of the optimal control strategy, the City_Cycle cycle condition is selected for simulation analysis. The City_Cycle cycle condition is the typical urban bus cycle-driving condition in China. Table 4 shows the main parameters of the City_Cycle cycle.
Figure 16 shows the vehicle speed curve of a load-isolated electric bus.

It can be seen that the optimal control strategy satisfies the requirements of the working conditions very well, which indicates that the strategy has good working condition following characteristics. This is in line with the brake stability requirements.
5.2. Simulation Comparison between Original Strategy and Optimization Strategy under Single Working Conditions
In the single-mode simulation comparison, the initial SOC of the energy storage battery and the power battery is set to 50%. The vehicle runs a City_Cycle cycle condition. After the vehicle operating conditions are over, the motor torque, power battery input energy, power battery SOC, and energy recovery rate of the original control strategy and the optimized control strategy are compared. The energy flow during operation is shown in Figure 17.

5.2.1. Torque Comparison of the Drive Motor
The torque of the drive motor is shown in Figure 18. It can be seen from the figure that there is a greater increase in the torque of the optimal control strategy during braking. This shows that the control strategy is optimized to distribute more motor torque during braking.

5.2.2. Input Energy Comparison of Power Battery
The control strategy stipulates that the power battery is given priority to recover braking energy, as shown in Figure 19. The vehicle power battery embedded with the original control strategy recovered 7808.25 kJ of energy. The vehicle power battery embedded with the optimized control strategy recovered 11444.64 kJ of energy. The results show that the optimal control strategy can recover 3636.39 kJ more under the same road conditions.

5.2.3. SOC Comparison of Power Battery
Figure 20 shows a comparison of the power battery SOC. It can be seen from the diagram that the SOC of the vehicle power battery embedded with the optimal control strategy is 52.12% and the SOC of the vehicle power battery embedded with the original control strategy is 52.03%. The energy storage battery is the main power source, and the power battery is the auxiliary power source. Because the power battery preferentially recovers braking energy, the power battery SOC will rise. The results show that the power battery SOC of the optimized control strategy is 0.09% more than that of the original control strategy. This indicates that more regenerative braking ratios are allocated after optimization.

5.2.4. Energy Recovery Ratio Comparison
The energy recovery rates of the original control strategy and the optimized control strategy are shown in Table 5. The original control strategy was 35.79%, and the optimal control strategy was 52.88%. In terms of energy recovery, the optimal control strategy increased by 17.09%.
5.3. Simulation Analysis of Multicycle Conditions
The initial SOC of the energy storage battery and the power battery is set to 100%. Then, the load-isolated electric bus continues to operate in the City_Cycle condition until the set battery discharges the minimum SOC limit (Csoc is 20% and Dsoc is 10%). After such testing, the energy recovery, economy, battery SOC, and charged energy of the control strategy can be detected.
5.3.1. Energy Recovery Rate Simulation Comparison
Under the original control strategy, the vehicle recovered 202121.47 kJ of energy, the energy consumption was 713712.24 kJ, and the energy recovery rate was 28.32%. Under the optimized control strategy, the vehicle recovered 141605.15 kJ of energy, the energy consumption was 359696.05 kJ, and the energy recovery rate was 39.37%, as shown in Table 6. The results show that the vehicle energy recovery rate of the embedded optimization control strategy increased by 11.05%.
5.3.2. Optimization Control Strategy Economic Simulation
The economic simulation results are shown in Table 7. The original control strategy has a driving range of 237.50 km and a power consumption of 61.48 kWh per 100 kilometers. At the same time, the driving range of the optimized control strategy is 255.89 km, and the power consumption per 100 kilometers is 51.01 kWh. Compared with the original strategy, the driving range of the embedded optimization control strategy increased by 44 km, which is an increase of 7.74%. Moreover, the power consumption per 100 kilometers was reduced by 10.47 kWh, which is a reduction of 17.03%
5.3.3. Battery SOC Change of Optimized Control Strategy
It can be seen from Figure 21 that when the energy storage battery is discharged until the SOC is reduced to 10%, the power is supplied by the power battery. Figure 20 shows a graph of the battery SOC variation for an optimized control strategy. When the power battery is discharged until the SOC is reduced to 20%, it is no longer discharged. This result is in line with the expected control objectives.

5.3.4. Charging Energy of Optimized Control Strategy’s Battery
The charging energy of the optimized control strategy’s battery is shown in Figure 22. It can be seen from Figure 22 that both the energy storage battery and the power battery SOC are greater than 80% and that no energy is recovered. At 5000 s, the energy storage battery SOC is lower than 80%, and energy recovery begins. At 35000 s, the power battery SOC is lower than 80%, and energy recovery begins.

6. Conclusions
In this paper, a regenerative braking system for the load-isolated electric bus was designed, and it provides a basis for research on the vehicle braking system. A series composite braking control strategy was proposed as the basis for further research on the composite braking of load-isolated electric buses.
This research designed two fuzzy controllers with small and medium braking strengths to optimize the series composite braking control strategy. The simulation analysis results show that the energy recovery rate of the optimized control strategy increased by 11.05% and the driving range increased by 7.74%. The effectiveness and rationality of the designed fuzzy controller are well verified. The results lay a foundation for the study of regenerative braking in the future.
Data Availability
The data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
This research work was supported by the National Key Research and Development Program of China (no. 2017YFB0102004), Shandong Province Major Innovation Engineering Project Subject (no. 2017CSGC0502), and Shandong Province Key Research and Development Program (no. 2017GGX50106).