Gear Ratio Optimization along with a Novel Gearshift Scheduling Strategy for a Two-Speed Transmission System in Electric Vehicle

: A novel gearshift scheduling strategy has been framed for a two-speed transmission system in electric vehicles that can save energy during hilly driving and frequently changing driving conditions through e ﬃ cient electric motor operation. Unlike the traditional approach, the proposed gearshift strategy is based on the preferred vehicle speed range, vehicle acceleration, and road grade to ensure desired vehicle performances with minimum energy consumption. Meanwhile, the vehicle speed range is chosen around the electric motor rated speed, and two gearshift schedules in relation to vehicle acceleration and road grade are developed based on the motor torque generating capacity and e ﬃ ciency. Appropriate gear is selected through a combined assessment of the required vehicle speed, acceleration, and road grade information. A guideline is developed and explained for the primary gearshift schedule. Next, the gear ratios and gearshift schedules are optimized combinedly in a Simulink environment using the gradient descent method and pattern search method on three driving cycles separately. Depending on the driving scenarios, around 4% to 7.5% energy saving has been experienced through optimization, while the gear ratios and gearshift schedule in relation to the road grade are found to be major contributors to the vehicle economic driving compared to that with the gearshift schedule for vehicle acceleration.


Introduction
Incorporating the multi-speed transmission system in electric vehicles (EVs) could be a solution to achieve better drivetrain efficiency. In the review paper by Ahssan et al. [1], it has been shown how the multi-speed transmission system outperforms the single-speed transmission system in EVs. However, the multi-gear system introduces additional complexities such as the selection of optimal gear ratios, appropriate gearshift schedule and shifting control, torque interruption during gear shifting, increased mass, and transmission loss. In this study, an appropriate gearshift schedule strategy, and then, optimization of both gear ratios and the gearshift schedule will be investigated.
A higher value of gear ratio contributes to achieving higher gradeability and faster acceleration, while a lower value of gear ratio determines the maximum vehicle speed limit and reduces energy consumption [2,3]. The gearshift schedule dictates the initiation of gear change i.e., either upshift or downshift depending on the demand of vehicle performances. Vehicle driving economy or dynamic  [7] 5.9/3.11 Throttle Position Motor speed 2-Speed [4] 10.63/5.12 Throttle opening Vehicle speed 2-Speed [9] 9.39/4.83 Throttle angle Vehicle speed 2-Speed [13] 11.47/4.64 Throttle angle Vehicle speed 2-Speed [14] 9.8/3.56 Accelerator opening Vehicle speed 2-Speed [2] -Driver demand Vehicle speed 2-Speed [15] -Driver demand Vehicle speed In the above-mentioned research papers, maximum gradeability is set as one of the dynamic performance targets, but energy consumption on a route with frequently varied altitude has not been studied. Moreover, standard drive cycles i.e., New European Drive Cycle (NEDC), Urban Dynamometer Driving Schedule (UDDS), etc. are based on the flat road [16]. Recently, Han et al. [17] presented a systematic framework to optimize the gear numbers, gear ratios, and gearshift schedule in the transmission system for EVs considering several driving conditions. In another study, Li et al. [18] conducted parameter optimization of both the electric motor and a two-speed transmission system in an EV. Both research groups used the road inclination to set the dynamic performance target and the upper limit of the gear ratios during the optimization process. However, road inclination is ignored in the driving cycles. A single gearshift schedule may not be equally suitable or efficient for all drive cycles without compromising the vehicle performances. It is also difficult to come up with a single optimized gearshift schedule that would ensure efficient motor operation in all possible driving scenarios. The suggested gearshift schedule strategy will assess motor speed, vehicle acceleration and road grade to initiate the gearshift and will address the mentioned shortcomings. In the next section, typical models of major powertrain components in EV are briefly described along with a detailed explanation of shift schedule strategy for a two-speed transmission unit that can be further extended to a multi-speed transmission system in electric vehicle. In Section 3, a combined gear ratio and shift schedule optimization process is explained. The results are analyzed in the fourth section followed by conclusion in Section 5.

EV Powertrain
Major components of EV powertrain include a battery, an electric motor, and a transmission system. All the components are modeled within MATLAB/Simulink environment. In the following sub-sections, the functionality of the battery and electric motor blocks is briefly explained along with their specifications. The transmission unit is described in detail as the focus of this study. To evaluate the performances of the proposed gearshift schedule strategy for a two-speed transmission system in electric vehicles, the parameters of a large-size vehicle shown in Table 2 have been collected from a published research paper [6]. The inertial mass of the vehicle has been estimated as m inertia = 336.45 kg based on the information of inertia of parts and associated formula provided in research papers by Lucente et al. [19] and Saini et al. [12], respectively. The vehicle model receives the information of the input drive cycle speed, road grade, and wind speed as inserted in the Simulink constant block and generates an output vehicle speed based on the following equation of motion i.e., Equation (1), where η PT is the powertrain efficiency and GR i is the gear ratios [13]. Here, GR i is representing the product of actual gear ratios and the final drive ratio, i f . Meanwhile, wind speed is assumed to be insignificant or zero; drive cycles and road grade information will be further explained later in this study.

Battery
The battery unit stores the electric energy as presented in Table 3 and supplies it to the traction motor. However, the battery block is connected to the electric motor block as a backward-facing module in the EV model for simulation. This unit receives power demand P M from the electric motor and adjusts it with its energy capacity over time. When the battery state-of-charge, SOC, falls below a certain value i.e., 5%, it stops the simulation as an indication of a critical level of SOC. A single battery cell is modeled based on the lookup table data on cell voltage versus battery state-of-charge, SOC. Next, the battery output voltage is estimated according to the total number of battery cells considering both parallel and series arrangement.
Equations (2)-(5) explain how the battery state-of-charge, SOC, and battery voltage, V bat , are estimated. Here, I bat is the battery current, while N cell_s and N cell_p are the number of battery cells connected in series and parallel. The battery open-circuit voltage, V OC , and battery internal resistance, R int , are the function of battery state-of-charge, SOC, and the surrounding temperature, temp, as shown in Equations (6) and (7), and these are recalled from the look-up table data retained in the battery block in Simulink.

Electric Motor
A map-based motor model from the Simulink Powertrain Blockset toolbox is used for simulation. The motor model is integrated with two lookup tables in Simulink. One lookup table provides maximum torque vs. speed data while the other holds motor efficiency data at numerous torque-speed combinations as referred in Equations (8) and (9). Different types of electric motors, i.e., a permanent magnet motor, switched-reluctance motor, induction motor, etc. can be chosen for EVs depending on their robustness and availability. Characteristics of these types of electric motors are discussed in detail in the literature [20][21][22][23]. Interested readers are referred to [24][25][26] for in-depth knowledge on the control and modeling of electric motor for EVs. While some research can be found with more than one electric motor as the propulsion device in EVs [27,28], a single electric motor is considered in this study as the propulsion device. Table 4 shows the motor capacity that is collected from a published article by Zhu et al. [8]. The motor power, P M , comprises the required motor output power at the input side of the transmission system, i.e., mechanical power, P Mech , and loss of power, P Loss , to generate that required mechanical power at any combination of the required electric motor output torque, T M , and speed, n M , as shown in Equation (10). Being a secondary module in the EV model, the battery block shows the drop in the battery state-of-charge, SOC, over time based on the motor power, P M .
Efficiency data are extracted from the available motor efficiency map shown in Figure 1. During the simulation, the motor efficiency value at any torque-speed combination is determined through linear interpolation/extrapolation. The electric motor unit receives an acceleration command from the driver module, generates the required torque based on lookup table data, and transfers it to the transmission unit. A fully pressed acceleration pedal is linked to the maximum motor output torque generation at Energies 2020, 13, 5073 5 of 24 any speed, while any fraction of the acceleration pedal is decoded as how much motor output torque is needed to be scaled down from its maximum limit.  Table 4. Traction motor parameter [8]. Rated Power  40  kW  Rated Torque  127  Nm  Base Speed  3000  rpm  Maximum Power  80  kW  Maximum Torque  255  Nm  Maximum Speed  9000  rpm Efficiency data are extracted from the available motor efficiency map shown in Figure 1. During the simulation, the motor efficiency value at any torque-speed combination is determined through linear interpolation/extrapolation. The electric motor unit receives an acceleration command from the driver module, generates the required torque based on lookup table data, and transfers it to the transmission unit. A fully pressed acceleration pedal is linked to the maximum motor output torque generation at any speed, while any fraction of the acceleration pedal is decoded as how much motor output torque is needed to be scaled down from its maximum limit.

Transmission
In conventional internal combustion engine (ICE) vehicles, higher number of gears are used in the transmission system to reduce fuel consumption [29]. However, that is not the case with EV because of the torque-speed characteristics of the electric motor i.e., high torque at the lower speed region and a constant power region throughout the extended speed range [21], and fewer gears could serve the purpose of achieving the required vehicle performances [13]. A two-speed automatic manual/mechanical transmission (AMT) system is considered for the EV model in this study. For an EV powertrain, several types of multi-speed transmission system are available to choose such as automatic transmission (AT) [30,31], continuous variable transmission (CVT), dual-clutch transmission (DCT) [32,33], inverse automatic mechanical transmission (IAMT) [4,34], uninterrupted mechanical

Transmission
In conventional internal combustion engine (ICE) vehicles, higher number of gears are used in the transmission system to reduce fuel consumption [29]. However, that is not the case with EV because of the torque-speed characteristics of the electric motor i.e., high torque at the lower speed region and a constant power region throughout the extended speed range [21], and fewer gears could serve the purpose of achieving the required vehicle performances [13]. A two-speed automatic manual/mechanical transmission (AMT) system is considered for the EV model in this study. For an EV powertrain, several types of multi-speed transmission system are available to choose such as automatic transmission (AT) [30,31], continuous variable transmission (CVT), dual-clutch transmission (DCT) [32,33], inverse automatic mechanical transmission (IAMT) [4,34], uninterrupted mechanical transmission (UMT) [35], automatic mechanical transmission (AMT), etc. Each type of transmission system has its own benefits or limitations. Chai et al. [36] and Kim et al. [37] made enough reasonings in favor of using AMT in EVs. Some of those are features of ATs with higher-energy economy, higher controllability, converting existing manual transmission (MT) into AMT causing lower development cost, etc.
As shown in Figure 2, an AMT can generally be formed by adding an actuator (i.e., hydraulic actuator or small electric motor) and sensors to the traditional MT and by removing the existing clutch arrangement. The gearshift is conducted through an actuator according to the command from the transmission control unit (TCU). As a result of the high torque-speed controllability of the traction motor, it is possible to adjust the torque and speed between the rotor and input shaft of transmission during gear disengagement and engagement in the gearshift process to minimize the torque interruption at the vehicle wheels. Based on these criteria, AMT is considered as a reasonable choice for this study. actuator or small electric motor) and sensors to the traditional MT and by removing the existing clutch arrangement. The gearshift is conducted through an actuator according to the command from the transmission control unit (TCU). As a result of the high torque-speed controllability of the traction motor, it is possible to adjust the torque and speed between the rotor and input shaft of transmission during gear disengagement and engagement in the gearshift process to minimize the torque interruption at the vehicle wheels. Based on these criteria, AMT is considered as a reasonable choice for this study. The role of the transmission system is to ensure enough wheel torque to run the vehicle at maximum road grade and at desired maximum speed and to achieve faster acceleration with minimum energy consumption. Based on recent literature on the multi-gear transmission system in an EV, Table 5 demonstrates the targets of the vehicle dynamic performances that have been set to evaluate the performances of the proposed gearshift scheduling strategy for two-speed AMT in an EV. A comprehensive discussion can be found in [10,38], where the efficiency of a two-speed DCT may vary from 86% to 95% depending on the gear number and other factors i.e., viscous loss, wet clutch loss, synchronizer mechanism, differential, etc. However, in this study, the efficiency of the two-speed AMT system is considered as 90% on gear 1 and 95% on gear 2. In the following sub-sections, the novel gearshift strategy is explained about how this would meet up the target vehicle performances as well as save energy consumption at different driving conditions.  The role of the transmission system is to ensure enough wheel torque to run the vehicle at maximum road grade and at desired maximum speed and to achieve faster acceleration with minimum energy consumption. Based on recent literature on the multi-gear transmission system in an EV, Table 5 demonstrates the targets of the vehicle dynamic performances that have been set to evaluate the performances of the proposed gearshift scheduling strategy for two-speed AMT in an EV. A comprehensive discussion can be found in [10,38], where the efficiency of a two-speed DCT may vary from 86% to 95% depending on the gear number and other factors i.e., viscous loss, wet clutch loss, synchronizer mechanism, differential, etc. However, in this study, the efficiency of the two-speed AMT system is considered as 90% on gear 1 and 95% on gear 2. In the following sub-sections, the novel gearshift strategy is explained about how this would meet up the target vehicle performances as well as save energy consumption at different driving conditions.

Proposed Shift Scheduling Strategy
A three-parameter-based gearshift strategy has been formulated to choose appropriate gear at various driving conditions. These parameters are vehicle speed, road grade, and vehicle acceleration, which are estimated in real time. Vehicle speed information can be directly available through the speedometer, while acceleration can be calculated by differentiating the current velocity. Although real-time road grade information is not common in most of the present production vehicle control systems [39], road grade estimation is becoming significant for advanced driver assistance system (ADAS), fuel economy, intelligent energy management strategy (EMS), and autonomous cars [40,41]. In the literature, several approaches for real-time road grade estimation can be found such as recursive least squares (RLS) using Energies 2020, 13, 5073 7 of 24 vehicle longitudinal dynamic models [42,43], adaptive extended Kalman filter (AEKF) [44], parallel mass and grade (PNG) estimation [45], grade estimation through orientation filter from inertial measurement unit (IMU) data [40], etc. Within the scope of this study, it is assumed that any approach of road grade estimation would meet up the requirement for this study, and road grade information would be available as a feature of ADAS. The upper and lower bound of each parameter are set for upshift and downshift, respectively. Then, a set of shift logic initiates the change of gear considering the combined assessment of all three parameters.

Gearshift Schedule Parameter-Vehicle Speed
The vehicle speed range is set based on the rated motor speed, n b i.e., 3000 rpm, as shown in Table 4, because it is known that the motor runs most efficiently at its rated speed. For the preferred motor speed range, there will be two sets of vehicle speed range due to considering two gears. Motor speed is not considered directly in decision making because of the sudden fluctuation of motor speed during gearshift to avoid the gear hunting problem. When the vehicle speed goes above the upshift speed of gear 1, then gear 2 will be allowed to select. Next, when the vehicle speed goes below the downshift speed of gear 2, then gear 1 will be the preferred choice.

Gearshift Schedule Parameter-Road Grade and Vehicle Acceleration
Initially, two separate shifting maps are developed to meet up the torque demands at various road grades and accelerations. To generate the gearshift maps, the maximum and rated capacity of the electric motor i.e., torque, speed, power, and motor efficiency at each gear, as shown in Figures 1 and 3, as well as common driving scenarios, industrial practice, and knowledge are the considerations. through the speedometer, while acceleration can be calculated by differentiating the current velocity. Although real-time road grade information is not common in most of the present production vehicle control systems [39], road grade estimation is becoming significant for advanced driver assistance system (ADAS), fuel economy, intelligent energy management strategy (EMS), and autonomous cars [40,41]. In the literature, several approaches for real-time road grade estimation can be found such as recursive least squares (RLS) using vehicle longitudinal dynamic models [42,43], adaptive extended Kalman filter (AEKF) [44], parallel mass and grade (PNG) estimation [45], grade estimation through orientation filter from inertial measurement unit (IMU) data [40], etc. Within the scope of this study, it is assumed that any approach of road grade estimation would meet up the requirement for this study, and road grade information would be available as a feature of ADAS. The upper and lower bound of each parameter are set for upshift and downshift, respectively. Then, a set of shift logic initiates the change of gear considering the combined assessment of all three parameters.

Gearshift Schedule Parameter-Vehicle Speed
The vehicle speed range is set based on the rated motor speed, i.e., 3000 rpm, as shown in Table 4, because it is known that the motor runs most efficiently at its rated speed. For the preferred motor speed range, there will be two sets of vehicle speed range due to considering two gears. Motor speed is not considered directly in decision making because of the sudden fluctuation of motor speed during gearshift to avoid the gear hunting problem. When the vehicle speed goes above the upshift speed of gear 1, then gear 2 will be allowed to select. Next, when the vehicle speed goes below the downshift speed of gear 2, then gear 1 will be the preferred choice.

Gearshift Schedule Parameter-Road Grade and Vehicle Acceleration
Initially, two separate shifting maps are developed to meet up the torque demands at various road grades and accelerations. To generate the gearshift maps, the maximum and rated capacity of the electric motor i.e., torque, speed, power, and motor efficiency at each gear, as shown in Figures 1  and 3, as well as common driving scenarios, industrial practice, and knowledge are the considerations.  An imaginary schedule line or a separation line as shown in Figure 4 is developed to identify the efficient operation area of the electric motor at each gear among the two consecutives. Several representative points are identified on both the road grade vs. vehicle speed and acceleration vs. vehicle speed graph. These points can be named as transition points between two consecutive gears. The gearshift schedule line is achieved by connecting those points in relation to road grade and vehicle acceleration. Then, to avoid the gear hunting problem, this schedule line will be split up into two lines by imposing a weighting factor. The line on the left side represents the downshift line, while the line on the right side stands for the upshift. representative points are identified on both the road grade vs. vehicle speed and acceleration vs. vehicle speed graph. These points can be named as transition points between two consecutive gears. The gearshift schedule line is achieved by connecting those points in relation to road grade and vehicle acceleration. Then, to avoid the gear hunting problem, this schedule line will be split up into two lines by imposing a weighting factor. The line on the left side represents the downshift line, while the line on the right side stands for the upshift. The gearshift schedule achieved at this stage can be considered as the primary gearshift schedule to run the simulation of the EV model. Next, the primary shift schedule will be further optimized through an appropriate optimization process. A guideline is articulated here about how the primary gearshift schedule can be generated with less effort.

•
Considering vehicle acceleration at different vehicle speed ranges and various road grades at different constant vehicle speeds, numerical analysis on the motor torque requirements has been conducted based on Equation (1). For any condition, there are three possible choices for a two-speed transmission system; i.e., the first is that only gear 1 can provide the required wheel torque at this vehicle speed, the second choice is only gear 2 can be the option, and the third is any gear can be chosen to supply the necessary torque and speed at the vehicle wheel. For the third scenario, the gear with the higher motor efficiency can be chosen to save energy consumption. Therefore, it is convenient to identify the three-speed regions such that each speed region can be associated with any of the above three choices as stated. At the low-speed region i.e., zero to approximately 20-25 km/h, gear 1 is an obvious option because of the high torque demand from the motor to meet up with the requirements at the initial vehicle acceleration or uphill driving scenario. Practically, gear 2 is not an option in this speed range, because the required motor output torque exceeds the maximum motor output torque limit in many cases or the motor operates at a very low efficiency region.

•
Similarly, at a high-speed region i.e., 65-70 km/h or above, gear 1 cannot be selected, because the required motor speed goes over the maximum speed limit of the motor. Therefore, gear 2 can be set as the default selection in this speed range. • Now, the gear selection problem can be narrowed down to the vehicle speed range between 20-25 km/h and 65-70 km/h where any gear could be a choice to meet up the requirements to run the vehicle. Here, the motor efficiency and rated motor capacity play the decision-making role for gear selection. Transition points for gear change are identified where the required motor output torque and speed at both gears are close to its rated capacity. Gear one is preferred at the vehicle speed below the speed corresponding to the transition points, and gear two is selected The gearshift schedule achieved at this stage can be considered as the primary gearshift schedule to run the simulation of the EV model. Next, the primary shift schedule will be further optimized through an appropriate optimization process. A guideline is articulated here about how the primary gearshift schedule can be generated with less effort.

•
Considering vehicle acceleration at different vehicle speed ranges and various road grades at different constant vehicle speeds, numerical analysis on the motor torque requirements has been conducted based on Equation (1). For any condition, there are three possible choices for a two-speed transmission system; i.e., the first is that only gear 1 can provide the required wheel torque at this vehicle speed, the second choice is only gear 2 can be the option, and the third is any gear can be chosen to supply the necessary torque and speed at the vehicle wheel. For the third scenario, the gear with the higher motor efficiency can be chosen to save energy consumption. Therefore, it is convenient to identify the three-speed regions such that each speed region can be associated with any of the above three choices as stated. At the low-speed region i.e., zero to approximately 20-25 km/h, gear 1 is an obvious option because of the high torque demand from the motor to meet up with the requirements at the initial vehicle acceleration or uphill driving scenario. Practically, gear 2 is not an option in this speed range, because the required motor output torque exceeds the maximum motor output torque limit in many cases or the motor operates at a very low efficiency region.

•
Similarly, at a high-speed region i.e., 65-70 km/h or above, gear 1 cannot be selected, because the required motor speed goes over the maximum speed limit of the motor. Therefore, gear 2 can be set as the default selection in this speed range. • Now, the gear selection problem can be narrowed down to the vehicle speed range between 20-25 km/h and 65-70 km/h where any gear could be a choice to meet up the requirements to run the vehicle. Here, the motor efficiency and rated motor capacity play the decision-making role for gear selection. Transition points for gear change are identified where the required motor output torque and speed at both gears are close to its rated capacity. Gear one is preferred at the vehicle speed below the speed corresponding to the transition points, and gear two is selected when the vehicle speed goes over the speed at the transition points. Next, the imaginary separation line can be achieved by connecting the transition points.

•
An imaginary separation line needs to be split up into two lines i.e., upshift and downshift lines, to avoid frequent gearshift through applying an appropriate weighting factor. Zhang et al. [46] used a hysteresis strategy where the gearshift line is increased by 20% to set the upshift line and is decreased by 20% for the downshift line. Another common technique of implementing a weighting factor is to apply a penalty factor within the optimization algorithm to prevent the repeated gearshift [47,48]. Although these techniques are applied on throttle demand vs. vehicle speed-based gearshift schedule and cannot be directly linked to the proposed schedule in relation to road grade and vehicle acceleration, these could be considered as a general guideline. In this study, the weighting factor is set manually considering mainly two factors [49], i.e., avoiding frequent gearshift to ensure rider comfort and kick-downshifting during high acceleration demand. Therefore, the weighting factor is not necessarily the same at all vehicle speeds. However, a buffer zone of 10-40% has been maintained between the downshift line and upshift line that will be further optimized through applying an appropriate optimization method in the later part of this paper. • Following the above steps, primary gearshift schedules in relation to the vehicle acceleration and road grade have been shown in Figures 5 and 6, respectively.
separation line can be achieved by connecting the transition points. • An imaginary separation line needs to be split up into two lines i.e., upshift and downshift lines, to avoid frequent gearshift through applying an appropriate weighting factor. Zhang et al. [46] used a hysteresis strategy where the gearshift line is increased by 20% to set the upshift line and is decreased by 20% for the downshift line. Another common technique of implementing a weighting factor is to apply a penalty factor within the optimization algorithm to prevent the repeated gearshift [47,48]. Although these techniques are applied on throttle demand vs. vehicle speed-based gearshift schedule and cannot be directly linked to the proposed schedule in relation to road grade and vehicle acceleration, these could be considered as a general guideline. In this study, the weighting factor is set manually considering mainly two factors [49], i.e., avoiding frequent gearshift to ensure rider comfort and kick-downshifting during high acceleration demand. Therefore, the weighting factor is not necessarily the same at all vehicle speeds. However, a buffer zone of 10-40% has been maintained between the downshift line and upshift line that will be further optimized through applying an appropriate optimization method in the later part of this paper. • Following the above steps, primary gearshift schedules in relation to the vehicle acceleration and road grade have been shown in Figures 5 and 6, respectively.   lines, to avoid frequent gearshift through applying an appropriate weighting factor. Zhang et al. [46] used a hysteresis strategy where the gearshift line is increased by 20% to set the upshift line and is decreased by 20% for the downshift line. Another common technique of implementing a weighting factor is to apply a penalty factor within the optimization algorithm to prevent the repeated gearshift [47,48]. Although these techniques are applied on throttle demand vs. vehicle speed-based gearshift schedule and cannot be directly linked to the proposed schedule in relation to road grade and vehicle acceleration, these could be considered as a general guideline. In this study, the weighting factor is set manually considering mainly two factors [49], i.e., avoiding frequent gearshift to ensure rider comfort and kick-downshifting during high acceleration demand. Therefore, the weighting factor is not necessarily the same at all vehicle speeds. However, a buffer zone of 10-40% has been maintained between the downshift line and upshift line that will be further optimized through applying an appropriate optimization method in the later part of this paper.

•
Following the above steps, primary gearshift schedules in relation to the vehicle acceleration and road grade have been shown in Figures 5 and 6, respectively.

Gearshift Control
A flowchart is shown in Figure 7 that demonstrates the gear selection process. A set of control logics receives information related to the vehicle speed, road grade, and vehicle acceleration, and it evaluates these values against pre-defined upper and lower bounds of these parameters respectively. Next, it is checked whether the current values of the three parameters satisfy the upshift or downshift condition or hold the current gear.
It is assumed that the vehicle starts on gear 1. At every instance, the upper and lower range of the vehicle speed (i.e., up_speed and down_speed), road grade (i.e., up_grd and down_grd), and acceleration (i.e., up_acc and down_acc) are determined according to the current vehicle speed and gear number. Then, the current values of these parameters will be compared with either the upper range of all parameters respectively for upshift if the current gear is on gear 1 or the lower range of all parameters respectively for downshift if the current gear is on gear 2. The controller's decision could be either to change gear (upshift/downshift) or to stay on the current gear. However, in case of more than two gears in the transmission system, current parameter values are compared with both the upper and lower bounds for all middle gears i.e., the 2nd gear in a three-speed transmission system or the 2nd gear and 3rd gear in a four-speed transmission system, and so on. A minimum waiting period or shift interval is maintained between every consecutive gearshift to avoid frequent gear shifting as well as to ensure riding comfort. The shift interval must be greater than the maximum gearshift time, because the gearshift time is not necessarily the same for every gear change in a transmission system with more than two gears for both upshift and downshift [36,49]. Based on the study conducted by Zhang et al. [46], a minimum shift interval of 4 s is assumed to be sufficient for an acceptable level of riding comfort. It has been found through simulation that higher shift interval time leads to energy consumption. For example, increasing shift interval time from 1 to 4 s costs 1-1.5% additional energy consumption depending on different drive cycles.

Gearshift Control
A flowchart is shown in Figure 7 that demonstrates the gear selection process. A set of control logics receives information related to the vehicle speed, road grade, and vehicle acceleration, and it evaluates these values against pre-defined upper and lower bounds of these parameters respectively. Next, it is checked whether the current values of the three parameters satisfy the upshift or downshift condition or hold the current gear. It is assumed that the vehicle starts on gear 1. At every instance, the upper and lower range of the vehicle speed (i.e., up_speed and down_speed), road grade (i.e., up_grd and down_grd), and acceleration (i.e., up_acc and down_acc) are determined according to the current vehicle speed and gear number. Then, the current values of these parameters will be compared with either the upper range of all parameters respectively for upshift if the current gear is on gear 1 or the lower range of all parameters respectively for downshift if the current gear is on gear 2. The controller's decision could be either to change gear (upshift/downshift) or to stay on the current gear. However, in case of more than two gears in the transmission system, current parameter values are compared with both the upper and lower bounds for all middle gears i.e., the 2nd gear in a three-speed transmission system or the 2nd gear and 3rd gear in a four-speed transmission system, and so on. A minimum waiting period or shift interval is maintained between every consecutive gearshift to avoid frequent gear shifting as well as to ensure riding comfort. The shift interval must be greater than the maximum gearshift time, because the gearshift time is not necessarily the same for every gear change in a transmission system with more than two gears for both upshift and downshift [36,49]. Based on the study conducted by Zhang et al. [46], a minimum shift interval of 4 s is assumed to be sufficient for an acceptable level of riding comfort. It has been found through simulation that higher shift interval time leads to energy consumption. For example, increasing shift interval time from 1 to

Optimization Method
It is evident in the literature that each set of gear ratios is associated with a certain gearshift schedule that can ensure minimum energy consumption. In other words, it can be said that the gear ratios and gearshift strategy are mutually associated to achieve the desired vehicle dynamic and economic performances. In this study, both the gear ratios and gearshift schedule are optimized combinedly in the optimization process. Table 6 shows some approaches or methods of optimization of the multispeed transmission system in EVs in recent studies. All these methods are generally known as global optimization methods that are computationally intensive [13]. The global optimization method can generate more accurate results than the results achieved from the optimization method suitable for the local minimum value of the objective function. However, it is not shown how much more accurate results or benefits can be made through using global methods compared to those by using methods for local minimum values in this problem. In this study, both global and local optimization methods will be implemented to search for optimal gear ratios of the two-speed transmission system in EV, and then, the results will be compared. Table 6. Approaches for optimization of gear ratios in the multi-speed transmission system in EVs.

Optimization Method Cost Function Decision Variables Constraints
Multi-Objective Genetic Algorithm [7] Overall energy consumption Acceleration time, t 0-100 sec Gear ratios Step ratio: 1.9 Dynamic Programming [4] Energy consumption in terms of ∆SOC Gear ratios Max. speed, Grade ability, Step ratio: 1.7-1.8 Auto Search Method [9] Driving range Gear ratios and shift schedule Step ratio: 3.4 Genetic Algorithm [13] Mean motor efficiency Driving range Gear ratios and shift schedule Max. speed, Grade ability, Acceleration time Genetic Algorithm [14] Energy consumption in terms of driving range Gear ratios Max. speed, Grade ability, Acceleration time

Brute Force Iterative
Algorithm [2] Acceleration time, Grade ability, Top speed, Energy consumption Gear ratios Step ratio, Max. speed, Grade ability, Acceleration time

Brute Force Iterative
Algorithm [15] Acceleration time, Grade ability, Top speed, Energy consumption Gear ratios Step ratio, Max. speed, Grade ability, Acceleration time Gradient descent (GD) and pattern search (PS) methods are implemented to optimize the gear ratios and shift schedule of the two-speed transmission system of the EV model, and next, the results obtained from these methods are compared. Both methods can handle multi-variables and multiple objectives or cost functions. However, the GD method using the sequential quadratic programming algorithm is appropriate for optimization problems with the local minimum value of objective functions [50]. This method has the flexibility to consider any values of control variables within the bounds during each iteration, even the values on the boundary of variables. To the other end, the PS method is a derivative-free method to find the optimum solution and can handle discontinuous functions [51]. The PS method can be implemented to both global and local search. With a good initial start point, PS can provide a better global solution and may require a reduced computational time compared to the other traditional global optimization method i.e., evolutionary programming and genetic algorithm. To conduct the optimization, two objectives have been considered; one is the energy consumption in a complete drive cycle, and the other is the tracking error between the drive cycle input speed and vehicle output speed. The optimization problem is to find the minimum values of the objective functions for the desired set of gear ratios and optimized gearshift schedules in relation to vehicle acceleration and road grade.

Problem Formulation
In the optimization problem, state variables are represented by the objective functions i.e., energy consumption over the entire drive cycle and tracking error.

of 24
In Equations (11) and (12), f 1 or E dc stands for the total energy consumption in a complete drive cycle, where P M is the required motor output power at the time step, dt, and t is the end time of the drive cycle. The motor output power, P M , is the summation of the required mechanical power, P Mech , at the input side of the transmission system to run the vehicle and motor loss, P Loss , corresponding to the motor torque, T M , and speed, n M , as shown in Equation (10). The motor efficiency, η M , is recalled from the lookup table during the simulation. In Equations (11) and (13), f 2 represents the tracking error between the input drive cycle speed, V in , and output vehicle speed, V out . NEDC and UDDS drive cycles are used as the input to the EV model. The aim of the optimization is to minimize the values of these two objectives and to find out the corresponding values of gear ratios along with the optimal gearshift schedule.
The control variable u in Equation (14) includes gear ratios i.e., GR i , and shift factors i.e., k i , for both shifting maps in relation to vehicle acceleration and road grade for the upshift and downshift lines, respectively. Shift factors are introduced in such a way that these would optimize the gearshift lines within the middle region of both gearshift schedule maps i.e., around the 30 km/h to 70 km/h speed range of Figure 5 and around the 20 km/h to 80 km/h speed range of Figure 6. Since gear selection is obvious (i.e., either gear 1 or gear 2) on either side of the middle region of both schedules, these areas have been excluded in the optimization process. It is expected that this measure would eliminate some unnecessary computation during the optimization.

Boundary Conditions and Constraints
The upper and lower bound of gear ratios are set based on the target performances as presented in Table 5. The lower bound of the 1st gear ratio is determined by the maximum gradeability requirement and can be found using Equation (15).
The upper bound of the 2nd gear ratio is set based on the maximum vehicle speed requirement and is calculated according to Equation (16).

Maximum Motor Speed * Wheel Radius
Maximum Vehicle Speed GR 2max = (2 * π * n max /60) * r w V max /3.6 However, it is necessary to check whether the motor can generate enough torque required to run the vehicle at maximum speed. Equation (17) shows the torque ratio between the vehicle resistance torque at the desired top speed and the maximum motor torque at the maximum motor speed i.e., γ = Vehicle Resistance Torque at Maximum Speed Maximum Motor Torque at Maximum Motor Speed γ = (m * g * (µ * cos θ+sin θ)+0.5 * C d * ρ * A * (V max /3.6) 2 ) * r w T Mmax at n max GR 2max will be acceptable under the condition of γ ≤ GR 2max i.e., enough motor torque will be available at the vehicle's top speed. To set the upper bound of the 1st gear ratio, GR 1max , and the lower bound of the 2nd gear ratio, GR 2min , the idea of step ratio will be applied.
Step ratio (str), is the ratio between two consecutive gear ratios. The maximum step ratio is set to 3.4 based on a study conducted by Zhou et al. [9]. Now, GR 1max and GR 2min can be estimated by using Equation (18) i.e., Energies 2020, 13, 5073 13 of 24 Considering the above-mentioned assumptions and the vehicle parameters presented in Table 2, the range of both gear ratios, as shown in Equation (19), can be found with imposing the requirement of 8.23 ≤ GR 1 ≤ 20.51 2.42 ≤ GR 2 ≤ 6.03 (19) It has been mentioned in the earlier section that the traditional approach of imposing the weighting factor to create a buffer zone between the upshift and downshift lines has been considered as a primary guideline. In line with that, a 20% limit as shown in Equation (20) is set for each shifting line to move or shift either downward or upward from its primary position.
Apart from the boundary conditions, four additional constraints are set out in Equations (21)-(24) to prevent the generation of an unacceptable combination of control variables in the optimization, even if the values are within the defined boundary of each variable separately. As in Equation (21), a maximum limit of the step ratio, str, needs to be ensured to avoid an undesirable shift jerk or shock load on the gears during gearshift as well as riding discomfort. The maximum step ratio, str, is considered 3.4 as mentioned earlier in this section. The second constraint in Equation (22) ensures the difference between the vehicle upshift speed @GR 1 in km/h and downshift speed @GR 2 in km/h above a certain value, c, to prevent the gear hunting problem while c can be any small positive number i.e., 5 or less, etc. In other words, a higher value of c may have a negative impact on energy consumption. The third and fourth constraints in Equations (23) and (24) are about providing the flexibility of moving both upshift and downshift lines up to 20% downward or upward direction from its original position, while the buffer zone should not be less than that of the primary gearshift maps associated with the vehicle acceleration and road grade, respectively.

Simulation and Result Analysis
Simulation has been conducted on the EV model with a two-speed transmission system in the MATLAB/Simulink environment. All the powertrain components are connected within a forward-facing model as shown in Figure 8, while only the battery is connected to the electric motor as a backward-facing module to show the drop in battery state-of-charge. The EV model is developed based on the Powertrain Blockset library available in the Simulink library browser. The Stateflow library is employed to construct the gearshift control logic as outlined in Figure 7.
Optimization has been conducted in two phases. The gear ratios and the gearshift map in relation to vehicle acceleration are optimized in the first phase. Next, the gearshift schedule for the road grade is optimized separately because the available drive cycles are based on a flat road. To evaluate the vehicle performances, a standard driving scenario on a frequently changing road grade is not readily available. Tan et al. [16] presented a driving cycle with associated road grade information in their study. However, that driving scenario is applicable to dump trucks that operated on a very specialized structured route in the mining zone. In this study, customized road grade information is considered on a standard drive cycle to evaluate the primary gearshift schedule in relation to the road grade for the two-speed transmission system in an EV model.  Optimization has been conducted in two phases. The gear ratios and the gearshift map in relation to vehicle acceleration are optimized in the first phase. Next, the gearshift schedule for the road grade is optimized separately because the available drive cycles are based on a flat road. To evaluate the vehicle performances, a standard driving scenario on a frequently changing road grade is not readily available. Tan et al. [16] presented a driving cycle with associated road grade information in their study. However, that driving scenario is applicable to dump trucks that operated on a very specialized structured route in the mining zone. In this study, customized road grade information is considered on a standard drive cycle to evaluate the primary gearshift schedule in relation to the road grade for the two-speed transmission system in an EV model. Figure 9 displays NEDC and UDDS drive cycles that are used as the input to the EV model. In this phase, optimization is conducted considering two groups of control variables where only gear ratios are considered as the first group of control variables and both the gear ratios and shift factors of the shift map for acceleration are considered combinedly as the second group of control variables for each drive cycle. Being a global optimization method, for the second group of control variables, the PS method takes an excessively longer duration to conduct the optimization. To accelerate the optimization process, an exception is made to avoid the inconvenience of uncertain waiting time, and each set of control variables are optimized separately through the PS method A set of generic gear ratios i.e., 10.00 for and 5.20 for , have been considered to start the simulation. The initial value of each shift factor of the gearshift map for vehicle acceleration is set to 1.0 to ensure that the primary gearshift maps are considered at the beginning of the optimization process.  Figure 9 displays NEDC and UDDS drive cycles that are used as the input to the EV model. In this phase, optimization is conducted considering two groups of control variables where only gear ratios are considered as the first group of control variables and both the gear ratios and shift factors of the shift map for acceleration are considered combinedly as the second group of control variables for each drive cycle. Being a global optimization method, for the second group of control variables, the PS method takes an excessively longer duration to conduct the optimization. To accelerate the optimization process, an exception is made to avoid the inconvenience of uncertain waiting time, and each set of control variables are optimized separately through the PS method A set of generic gear ratios i.e., 10.00 for GR 1 and 5.20 for GR 2 , have been considered to start the simulation. The initial value of each shift factor of the gearshift map for vehicle acceleration is set to 1.0 to ensure that the primary gearshift maps are considered at the beginning of the optimization process. It has been noticed that the maximum tracking error is around 4 km/h for NEDC and 6.4 km/h for the UDDS drive cycle with optimized gear ratios as presented in Figure 10a,b respectively. Table  7 shows the energy consumption of the two-speed EV model before and after optimization on complete NEDC and UDDS drive cycles. It has been noticed that the maximum tracking error is around 4 km/h for NEDC and 6.4 km/h for the UDDS drive cycle with optimized gear ratios as presented in Figure 10a,b respectively. Table 7 shows the energy consumption of the two-speed EV model before and after optimization on complete NEDC and UDDS drive cycles. It has been noticed that the maximum tracking error is around 4 km/h for NEDC and 6.4 km/h for the UDDS drive cycle with optimized gear ratios as presented in Figure 10a,b respectively. Table  7 shows the energy consumption of the two-speed EV model before and after optimization on complete NEDC and UDDS drive cycles.   Several more observations can be made based on the results achieved through the optimization process.

Optimization of Gear Ratios and Gearshift Schedule for Vehicle Acceleration
• After optimization, both GD and PS methods offer almost similar improvements in terms of energy consumption i.e., 4-5% less energy consumption depending on the driving scenarios i.e., NEDC or UDDS, as plotted in Figure 11. Another reflection from the results in Table 7 is that compared to the impact on the energy consumption by the shift factors of the shifting map for acceleration, gear ratios play more dominating roles on vehicle economic performance.
• For both driving scenarios i.e., NEDC and UDDS, optimized gear ratios achieved through both GD and PS methods are practically the same with little difference (approximately 0.02%) in the second gear ratio. To the other end, the shift factors of the gearshift map for vehicle acceleration are same for the GD methods on either driving conditions, while some difference can be noticed in the results through the PS method. • A significant observation is that both optimization methods tried to move the gear ratios toward the lower bound of each gear ratio. This indicates that the lower value of gear ratios will contribute to economic driving. However, further lowering the value of gear ratios has been restricted by dynamic performance targets. Secondly, an upward move (20% for the GD method and nearly 15% for the PS method) based on the results in Table 7 of both upshift and downshift lines of the gearshift schedule for acceleration is experienced after optimization, while the buffer zone of primary gearshift schedule is maintained as illustrated in Figure 12. An exception with the PS method i.e., close to the primary gearshift schedule with a slightly wider buffer zone, is observed after optimization with the NEDC driving cycle. Although shift factors for vehicle acceleration have relatively less impact on energy consumption, the results could be viewed as a guideline to generate the primary gearshift schedule with more confidence.
Energies 2020, 13, x FOR PEER REVIEW 17 of 26 Figure 11. Energy-saving through optimization of gear ratios and shift factors of a gearshift schedule for acceleration combinedly for a two-speed transmission system in an EV model on NEDC and UDDS drive cycles.

Optimization of Gearshift Schedule for Road Grade
It is mentioned earlier that a standard driving scenario with a frequently varied road grade has not been found in the literature. Therefore, to evaluate the primary gearshift map for a road grade, the Economic Commission for Europe (ECE, Geneva, Switzerland) Extra-Urban driving cycle is chosen because this drive cycle covers a wider speed range from low to relatively high with reasonable speed variations, as shown in Figure 13a. The low-speed region of this drive cycle is

Optimization of Gearshift Schedule for Road Grade
It is mentioned earlier that a standard driving scenario with a frequently varied road grade has not been found in the literature. Therefore, to evaluate the primary gearshift map for a road grade, the Economic Commission for Europe (ECE, Geneva, Switzerland) Extra-Urban driving cycle is chosen because this drive cycle covers a wider speed range from low to relatively high with reasonable speed variations, as shown in Figure 13a. The low-speed region of this drive cycle is

Optimization of Gearshift Schedule for Road Grade
It is mentioned earlier that a standard driving scenario with a frequently varied road grade has not been found in the literature. Therefore, to evaluate the primary gearshift map for a road grade, the Economic Commission for Europe (ECE, Geneva, Switzerland) Extra-Urban driving cycle is chosen because this drive cycle covers a wider speed range from low to relatively high with reasonable speed variations, as shown in Figure 13a. The low-speed region of this drive cycle is associated with a higher road grade, and the high-speed region is associated with a lower road grade, as shown in Figure 13. In this phase of optimization, the first group of control variables comprises both generic gear ratios as well as the shift factors of both primary gearshift maps for vehicle acceleration and road grade at the beginning of optimization. The second group of control variables contains only the shift factors of a shifting map for road grade considering the fixed and previously optimized values of gear ratios and shift factors of a shifting map for vehicle acceleration. Figure 14a shows the deviation in vehicle speed tracking (i.e., around 5.5 km/h), which is similar or close to that deviation on flat road scenarios, as shown in Figure 10. Another observation in Figure 14b,c is that the first gear selection in the mid-region (i.e., between 120 s to 180 s) of the drive cycle is altered after optimization, and the second gear is chosen throughout most of the driving time. The optimization results are demonstrated in Table 8. In this phase of optimization, the first group of control variables comprises both generic gear ratios as well as the shift factors of both primary gearshift maps for vehicle acceleration and road grade at the beginning of optimization. The second group of control variables contains only the shift factors of a shifting map for road grade considering the fixed and previously optimized values of gear ratios and shift factors of a shifting map for vehicle acceleration. Figure 14a shows the deviation in vehicle speed tracking (i.e., around 5.5 km/h), which is similar or close to that deviation on flat road scenarios, as shown in Figure 10. Another observation in Figure 14b,c is that the first gear selection in the mid-region (i.e., between 120 s to 180 s) of the drive cycle is altered after optimization, and the second gear is chosen throughout most of the driving time. The optimization results are demonstrated in Table 8.
Energy saving through optimization is evident in Figure 15. The optimized gearshift schedule in relation to road grade shows a significant contribution to vehicle driving economy compared to the combined dominance through optimized gear ratios and gearshift schedule for vehicle acceleration.
An overall assessment of all the results of optimization in this study reflects that the gear ratios and gearshift schedule for the uphill driving scenario could be considered as the key parameters that have a major influence on the vehicle economic performance.    Energy saving through optimization is evident in Figure 15. The optimized gearshift schedule in relation to road grade shows a significant contribution to vehicle driving economy compared to the combined dominance through optimized gear ratios and gearshift schedule for vehicle acceleration. An overall assessment of all the results of optimization in this study reflects that the gear ratios and gearshift schedule for the uphill driving scenario could be considered as the key parameters that have a major influence on the vehicle economic performance. It is apparent from the results in Tables 7 and 8 that both the GD and PS methods generate almost similar optimized gear ratios and shift factors of gearshift schedule in relation to vehicle acceleration with some exceptions. In Table 8, optimized gear ratios through the PS method can be ignored because for all other cases, similar gear ratios as well as energy savings are experienced after optimization. For uphill driving conditions, some difference is found on the optimized gearshift schedule, although the equivalent amount of energy-saving is observed with both optimization methods. An optimized gearshift schedule through the GD method is shifted 20% in the upward direction from the primary position, while the optimized downshift and upshift lines through the PS method are shifted around 12.73% upward and 3.42% downward, respectively, from the primary position. In other words, the buffer zone is further widened through the PS optimization method, as shown in Figure 16, while having similar energy-saving achieved through the GD method.  It is apparent from the results in Tables 7 and 8 that both the GD and PS methods generate almost similar optimized gear ratios and shift factors of gearshift schedule in relation to vehicle acceleration with some exceptions. In Table 8, optimized gear ratios through the PS method can be ignored because for all other cases, similar gear ratios as well as energy savings are experienced after optimization. For uphill driving conditions, some difference is found on the optimized gearshift schedule, although the equivalent amount of energy-saving is observed with both optimization methods. An optimized gearshift schedule through the GD method is shifted 20% in the upward direction from the primary position, while the optimized downshift and upshift lines through the PS method are shifted around 12.73% upward and 3.42% downward, respectively, from the primary position. In other words, the buffer zone is further widened through the PS optimization method, as shown in Figure 16, while having similar energy-saving achieved through the GD method.
Energy saving through optimization is evident in Figure 15. The optimized gearshift schedule in relation to road grade shows a significant contribution to vehicle driving economy compared to the combined dominance through optimized gear ratios and gearshift schedule for vehicle acceleration. An overall assessment of all the results of optimization in this study reflects that the gear ratios and gearshift schedule for the uphill driving scenario could be considered as the key parameters that have a major influence on the vehicle economic performance. It is apparent from the results in Tables 7 and 8 that both the GD and PS methods generate almost similar optimized gear ratios and shift factors of gearshift schedule in relation to vehicle acceleration with some exceptions. In Table 8, optimized gear ratios through the PS method can be ignored because for all other cases, similar gear ratios as well as energy savings are experienced after optimization. For uphill driving conditions, some difference is found on the optimized gearshift schedule, although the equivalent amount of energy-saving is observed with both optimization methods. An optimized gearshift schedule through the GD method is shifted 20% in the upward direction from the primary position, while the optimized downshift and upshift lines through the PS method are shifted around 12.73% upward and 3.42% downward, respectively, from the primary position. In other words, the buffer zone is further widened through the PS optimization method, as shown in Figure 16, while having similar energy-saving achieved through the GD method.  It should be remembered that the road grade information associated with the ECE Extra-Urban driving cycle is not from any recognized source; rather, this driving scenario is implemented in this study to see how vehicle driving economy is impacted by the frequently changing road grade. Therefore, without focusing on the specific numerical results, it would be wise to emphasize what message these results are producing. On that note, an upward shift of the gearshift schedule for both the case with vehicle acceleration and road grade could have a positive impact on energy saving with a multi-speed transmission system in EV that would be further validated through the experiments to draw a more specific conclusion.

Comparison of Results with Conventional Approach
To have a comprehensive understanding, it is necessary to evaluate the proposed gearshift strategy with the conventional rule-based gearshift strategy. On that note, two gearshift schedules as shown in Figure 17a,b are chosen from two studies published by Ruan et al. [11] and Zhu et al. [8], respectively. It is understandable that the assumptions are not necessarily the same in every study. However, a similar vehicle size is assumed to be worth consideration as a basis to compare the performances of these gearshift strategies.
driving cycle is not from any recognized source; rather, this driving scenario is implemented in this study to see how vehicle driving economy is impacted by the frequently changing road grade. Therefore, without focusing on the specific numerical results, it would be wise to emphasize what message these results are producing. On that note, an upward shift of the gearshift schedule for both the case with vehicle acceleration and road grade could have a positive impact on energy saving with a multi-speed transmission system in EV that would be further validated through the experiments to draw a more specific conclusion.

Comparison of Results with Conventional Approach
To have a comprehensive understanding, it is necessary to evaluate the proposed gearshift strategy with the conventional rule-based gearshift strategy. On that note, two gearshift schedules as shown in Figure 17a,b are chosen from two studies published by Ruan et al. [11] and Zhu et al. [8], respectively. It is understandable that the assumptions are not necessarily the same in every study. However, a similar vehicle size is assumed to be worth consideration as a basis to compare the performances of these gearshift strategies. Significant findings can be drawn from Figure 18. Around 3.0-4.2% and 3.7-3.9% less energy is consumed on NEDC and UDDS drive cycles respectively with a conventional gearshift schedule compared to that with the proposed gearshift schedule. However, the proposed gearshift strategy on uphill driving conditions can save energy up to 17.7% compared to that with the conventional approach. Although these results may vary slightly based on other considerations (i.e., EV powertrain components specifications, efficiency etc.), it is revealed that the proposed gearshift strategy could respond effectively to changing driving conditions, especially in hilly driving conditions, and it can also offer better economic performance. Significant findings can be drawn from Figure 18. Around 3.0-4.2% and 3.7-3.9% less energy is consumed on NEDC and UDDS drive cycles respectively with a conventional gearshift schedule compared to that with the proposed gearshift schedule. However, the proposed gearshift strategy on uphill driving conditions can save energy up to 17.7% compared to that with the conventional approach. Although these results may vary slightly based on other considerations (i.e., EV powertrain components specifications, efficiency etc.), it is revealed that the proposed gearshift strategy could respond effectively to changing driving conditions, especially in hilly driving conditions, and it can also offer better economic performance.

Conclusions
A novel gearshift schedule has been presented and explained in the first part of this study. While meeting up the dynamic performance requirements, this gearshift schedule technique is

Conclusions
A novel gearshift schedule has been presented and explained in the first part of this study. While meeting up the dynamic performance requirements, this gearshift schedule technique is expected to select the appropriate gear in relation to acceleration demand and change in road grade. Next, the gear ratios of a two-speed transmission system in EV and a gearshift schedule for both vehicle acceleration and road grade are optimized combinedly through two optimization methods i.e., GD and PS methods. After optimization of the two-speed transmission system, around 4% to 7.5% energy saving has been found, depending on the different driving scenarios considering both flat road and uphill road conditions, according to the results in Tables 7 and 8. Contributions on energy economy through the gear ratios and gearshift schedule for both vehicle acceleration and road grade have been assessed separately. It is revealed that compared to the impact of the gearshift schedule for acceleration, gear ratios play more dominating roles on vehicle energy consumption on flat road conditions. In uphill driving scenarios, the impact of gearshift schedule in relation to road grade on vehicle economic performance can even be more than the combined effect of gear ratios and gearshift map for acceleration. Apart from this, a comparison is conducted on the performance between the proposed gearshift strategy and the conventional approach. It is found that the proposed gearshift strategy shows a better response on hilly road compared to that with a conventional rule-based gearshift strategy.
The next significant finding of this study is about how the gear ratios and gearshift schedules are tuned after the implementation of optimization methods. It is observed that the values of both gear ratios are moving toward their lower boundaries based on the results in Tables 7 and 8 after optimization, while the primary gearshift schedule for vehicle acceleration is shifted in the upward direction, extending the working area of the second gear, as revealed in Figure 12. For the uphill driving conditions, the primary gearshift schedule in relation to the road grade has been shifted in the upward direction after optimization through the GD method, similar to that observed with the shifting map for acceleration. However, the gearshift schedule with a wider buffer zone is observed after implementing the PS method with equivalent energy consumption as found through implementing the GD method. Through this analysis, it has been explained how the shifting of gearshift schedule lines has an impact on the energy consumption of EVs with a multi-speed transmission system.

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