En-Route Battery Management and a Mixed Network Equilibrium Problem Based on Electric Vehicle Drivers’ En-Route Recharging Behaviors

With the rapidly increasing number of electric vehicle users, in many urbans transport networks, there are mixed traffic flows (i.e., electric vehicles and gasoline vehicles). However, limited by driving ranges and long battery recharging, the battery electric vehicle (BEV) drivers’ route choice behaviors are inevitably affected. This paper assumes that in a transportation network, when BEV drivers are traveling between their original location and destinations, they tend to select the path with the minimal driving times and recharging time, and ensure that the remaining charge is not less than their battery safety margin. In contrast, gasoline vehicle drivers tend to select the path with the minimal driving time. Thus, by considering BEV drivers’ battery management strategies, e.g., battery safety margins and en-route recharging behaviors, this paper developed a mixed user equilibrium model to describe the resulting network equilibrium flow distributions. Finally, a numerical example is presented to demonstrate the mixed user equilibrium model. The results show that BEV drivers’ en-route recharging choice behaviors are significantly influenced by their battery safety margins, and under the equilibrium, the travel routes selected by some BEV drivers may not be optimal, but the total travel time may be more optimal.


Introduction
Thanks to the rapid development of battery technologies and increasing price of crude oil [1][2][3][4], in recently years, the number of battery electric vehicles (BEVs) is increasing rapidly in many cities in the world. However, for many BEV users, although the number of charging stations is increasing in the urban area, they still fear that their vehicle may run out of power when driving between their original location and destination. This fear, known as range anxiety [5][6][7], is not only experienced by BEV drivers, but also by alternative-fuel vehicles drivers [8,9]. Nevertheless, for alternative-fuel vehicle drivers, range anxiety can be easily eliminated by taking onboard spare fuel cans. However, for BEV drivers, it seems unrealistic to take onboard spare batteries for replacing. Furthermore, to ease range anxiety, many BEV drivers are likely to set a battery safety margin to ensure the remaining battery can finish the trip [10][11][12]. Moreover, range anxiety also makes some BEV drivers prefer the path on which they can easily access charging stations. In a word, range anxiety will inevitably affect BEV drivers' travel choice behaviors. Thus, in the future, with BEVs widely participating in urban transportation, the equilibrium distribution of the urban traffic flow will be influenced.
However, in the research of traffic assignment, the assumptions of travelers' choice behaviors play an important role in solving network equilibrium problems [13][14][15][16][17]. Thus, in order to understand how 2 of 14 BEVs influence the equilibrium distribution of the urban traffic flow, the assumptions of BEV drivers' travel choice behaviors and the assumptions of BEVs' structure and characteristics will be important. For example, by considering electric vehicles' limited driving range and drivers' recharging behaviors, He et al. [18] proposed a network equilibrium model. In their research, they assumed that there are only BEVs in the network and these BEVs are homogeneous with the same battery size and same initial state of battery, and that BEV drivers are rational. These assumptions are restrictive and may lead a discrepancy between the research results and reality. Then, by considering stochastic range anxiety, Xie et al. [19] developed a network equilibrium with electric vehicles. In their research, they assumed that the variation of driving distance limits BEV drivers' perception and causes battery depletion of BEV drivers with heterogeneous perception errors and risk-taking behaviors. In addition, Jiang and Xie [20] and Xu et al. [21] considered a mixed network equilibrium problem with two vehicle types, i.e., gasoline vehicles and electric vehicles. In the former, a path-constrained traffic assignment is formulated; in the latter, a nonlinear minimization model in terms of path flows is formulated. The above studies lack of consideration on BEV drivers' recharging psychological barriers.
Becker [22] suggested that the utility of time spent at a meaningful activity is larger than the utility of time spent at a non-meaningful activity. Since recharging activity is an additional activity and makes BEV drivers' total travel time become longer, for BEV drivers who need to rechange their BEVs at en-route charging stations, their perception of disutility of unit driving time may be different from their perception of disutility of unit recharging time. Thus, to be more in line with reality, in this paper, we assume that there are two kinds of vehicles in the network, i.e., gasoline vehicles and electric vehicles; moreover, BEV energy consumption is not influenced by traffic congestion and is independent of traffic flows [5,23,24]. Then, for battery electric vehicles, assuming that these BEVs have the same battery size but different initial state of battery, and for BEV drivers who need recharging en-route, it is assumed that they are heterogenous with different risk attitudes (this can be reflected by different battery safety margins) and that their perceptions of disutility of unit recharging time is smaller than the perception of disutility of unit travel time. That is, if the value of the travel time coefficient is equal to the value of the recharging time coefficient, there must be an error in some BEV drivers' minimal recharging time perception that makes their actual minimal recharging time longer. Based on the above assumptions and considering BEV drivers' range anxiety and recharging behavior, we formulate a mixed network equilibrium model to deal with the traffic assignment of gasoline vehicles and electric vehicles.
In summary, the contribution of this paper is that by considering different types of drivers with different risk attitudes, different perception errors, and different initial states of battery, we developed a mixed user equilibrium model to describe gasoline vehicle drivers' route choice behavior, BEV drivers' en-route recharging behaviors, the utilization of charge stations, and the distribution of traffic flow in the network under a long-term equilibrium. The results of this paper can help to understand how BEV drivers' battery management strategies (e.g., battery safety margin setting and en-route recharging behavior) influence urban transport systems.
The remainder of this paper is organized as follows: Section 2 formulates a base model that describes a mixed network equilibrium, and then a solution procedure is proposed. In Section 3, a numerical example is presented, and the numerical example results and some conclusions are also presented. Finally, Section 4 concludes the paper.

The Characterization of a Traffic Network
Let G = (N, A) denote a traffic network where N is set of nodes and A is set of links. Let W be the set of origin-destination pairs (hereinafter referred to as "OD pairs") and w ∈ W be the OD pair index. Let R w and D w be the set of paths and the travel demand between OD pair w, and r ∈ R w the path index. We denote a link as a ∈ A, and f a is the corresponding link flow. Let δ w a,r be the path-link incidence. If path r traverses link a, then δ w a,r = 1, otherwise δ w a,r = 0. Let f a be the traffic flow on link a and z a be the travel distance of link a. The travel time of link a is a strictly increasing function of the flow in the link, for example, in this study, we used the following Bureau of Public Roads (BPR) function: where t 0 a is link a's free-flow travel time, and c a is link a's capacity.

Definition and Formulation of a Mixed User Equilibrium
In this paper, it is assumed that there are two kinds of vehicle drivers in the network, i.e., BEV drivers and gasoline vehicle drivers. For BEV drivers, they many need rechange their BEVs at one en-route charging station to avoid running out of charge before reaching their destinations. It is worth mentioning that this paper assumes that in order to reach the destination as quickly as possible, the BEVs will not be fully charged at the en-route charging station; BEV drivers will recharge a small amount of electricity (this also suggests that the values of travel time are higher for these BEV drivers). In addition, due to the uncertainty of fuel consumption, BEV drivers may avoid using all of the power supplied by their batteries when they arrive at the destinations. Instead, they prefer to set a battery safety margin that they keep the remaining battery higher than [10]. Thus, we assumed that when BEV drivers arrive their destinations, the remaining battery is no less than the battery safety margin. This battery safety margin will ensure that they can find the next charging station to keep the BEV going. Note that in this paper, the minimal recharging amount of electricity depends on the remaining distance traveled and BEV drivers' battery safety margins. Figure 1 is an example used to illustrate the en-route recharging behavior. As shown in Figure 1, the BEV's initial state of battery is 6 kWh, the energy consumption rate is 0.3 kWh/km [18], and the BEV driver's battery safety margin is 3 kWh. Then, under this setting, this BEV driver will choose path 2 and the amount of recharging electricity will be 6.9 kWh. On the other hand, for gasoline vehicle drivers, it is assumed that they have enough fuel. increasing function of the flow in the link, for example, in this study, we used the following Bureau of Public Roads (BPR) function: where is link 's free-flow travel time, and is link 's capacity.

Definition and Formulation of a Mixed User Equilibrium
In this paper, it is assumed that there are two kinds of vehicle drivers in the network, i.e., BEV drivers and gasoline vehicle drivers. For BEV drivers, they many need rechange their BEVs at one en-route charging station to avoid running out of charge before reaching their destinations. It is worth mentioning that this paper assumes that in order to reach the destination as quickly as possible, the BEVs will not be fully charged at the en-route charging station; BEV drivers will recharge a small amount of electricity (this also suggests that the values of travel time are higher for these BEV drivers). In addition, due to the uncertainty of fuel consumption, BEV drivers may avoid using all of the power supplied by their batteries when they arrive at the destinations. Instead, they prefer to set a battery safety margin that they keep the remaining battery higher than [10]. Thus, we assumed that when BEV drivers arrive their destinations, the remaining battery is no less than the battery safety margin. This battery safety margin will ensure that they can find the next charging station to keep the BEV going. Note that in this paper, the minimal recharging amount of electricity depends on the remaining distance traveled and BEV drivers' battery safety margins. Figure 1 is an example used to illustrate the en-route recharging behavior. As shown in Figure 1, the BEV's initial state of battery is 6 kWh, the energy consumption rate is 0.3 kWh/km [18], and the BEV driver's battery safety margin is 3 kWh. Then, under this setting, this BEV driver will choose path 2 and the amount of recharging electricity will be 6.9 kWh. On the other hand, for gasoline vehicle drivers, it is assumed that they have enough fuel. When traveling between their original location and destination, for each driver, it is assumed that they select the route with minimal travel costs, which includes the cost of travel time and electricity. It is worth mentioning that the cost of travel time is much higher than the cost of electricity [25]. Thus, we only consider travel time cost when drivers select paths. Notably, for BEV drivers who need recharging en-route, their total travel time includes driving time and recharging time. Moreover, different from others' research, in our based model, we consider heterogenous BEV drivers with different risk attitudes and different levels of perception errors regarding the values of recharging time, and the BEVs have the same battery size but a different initial state of battery. Moreover, we further assumed that in the transport network, there is only a finite number of charging stations, and these charging stations are located at certain nodes of the network. Note that a vehicle traveling along When traveling between their original location and destination, for each driver, it is assumed that they select the route with minimal travel costs, which includes the cost of travel time and electricity. It is worth mentioning that the cost of travel time is much higher than the cost of electricity [25]. Thus, we only consider travel time cost when drivers select paths. Notably, for BEV drivers who need recharging en-route, their total travel time includes driving time and recharging time. Moreover, different from others' research, in our based model, we consider heterogenous BEV drivers with different risk attitudes and different levels of perception errors regarding the values of recharging time, and the BEVs have the same battery size but a different initial state of battery. Moreover, we further assumed that in the transport network, there is only a finite number of charging stations, and these charging stations are located at certain nodes of the network. Note that a vehicle traveling along a path Energies 2020, 13, 4061 4 of 14 may not pass by a charging station. Thus, based on the above assumptions and considerations, we can obtain the following network equilibrium: Definition 1. In a long-term equilibrium, for the same type of drivers, all the utilized paths are usable and the total travel time costs of all the utilized paths of one OD pair are the same, that is, they are less than or equal to those of any unutilized usable paths of the same OD pair.
Here, a mathematic model will be proposed to describe the above long-term equilibrium. For convenience, it is assumed that there are K types of vehicle drivers in the network where K − 1 types belong to BEV drivers and they have different risk attitudes and different levels of perception of rechanging time. The remaining type belongs to gasoline vehicle drivers. Moreover, because BEV energy consumption is independent of traffic flow, letR w denote the set of all usable paths between OD pair w, thus, we can obtain following network equilibrium with heterogenous users (NE-HU): where, s w r,k is the minimal actual time (minutes) that drivers need to spend on recharging activity. Moreover, for gasoline vehicle drivers, there is s w r,k = 0. β k (≥ 1) is a coefficient and relates to the drivers' risk attitudes. α k (≥ 1) is a coefficient and relates to BEV drivers' perception errors of the values of recharging time. Note that limited by cognitive ability, drivers may not be able to accurately estimate their recharging times, thus, the "perception error" in this study is defined as the difference between a BEV driver's perceived value of recharging time and the actual value of recharging time. Moreover, α k = 1 indicates that a driver's perception of disutility of unit travel time is equal to their perception of disutility of unit recharging time, and the greater the value of α k , the smaller the perception of disutility of unit recharging time. g w k is travel demand (i.e., total vehicles) of type k between OD pair w, and f w r,k (veh/h) is the traffic flow of type k ∈ K on path r between OD pair w. Then, the above NE-HU model's Karush-Kuhn-Tucker (KKT) conditions can be written as: where ∂ f a ∂ f w a,k = δ w a,r and a t a ( f a )δ w a,r = t w r . Moreover, t w r is the travel time (minutes) of path r. Thus, Equation (7) can be further written as: Finally, based on the Equations (10)-(12), we can obtain the following conditions: Here, Equation (14) is a complementary condition.
This complementary condition guarantees that the traffic flow is only distributed on the usable paths.

Solution Procedure
The above NE-HU model is a convex program with linear constraints. We can solve it by the enumeration method. However, the work of enumerating all the usable paths will be tremendous. Thus, in this paper, we adopt an iterative solution procedure which was proposed by He et al. [18]. However, different from He's work, in this paper, we consider the perception recharging time instead of actual recharging time and heterogenous BEV drivers with the different risk attitudes and different levels of perception of recharging time. Thus, the sub-problem, finding the shortest usable path by considering BEV drivers' perception of recharging time (SP-PT), can be formulated as follows: where, l w i,k (kWh) is the recharging amount of electricity in theory at charging stations i for type k drivers, and β k l w i,k (kWh) is the actual recharging amount of electricity for type k drivers when considering their different risk attitudes. Let c i β k l w i,k represent the recharging time (minutes) that at node i ∈ N, type k drivers spend recharging β k l w i,k amount of electricity. Then, α k c i β k l w i,k is the perception of time for type k drivers to recharge β k l w i,k amount of electricity. In this paper, the following perception of charging time function is used, i.e., α k c i β k l w i,k = α k c 1 i + α k c 2 i β k l w i,k , and here c 1 i (minutes) represents the fixed recharging time; c 2 i (min/kWh) is a variable of recharging time and relies on the class of chargers, for example, direct current charging (fast charging) and alternating current charging (slow charging). However, in this study, it is assumed that the en-route charging stations only supply direct current charging facilities and the charging facilities are adequate. Note that if there is no charging station on node i, the value of α k c i β k l w i,k will be zero, and for gasoline vehicle drivers, the value of α k c i β k l w i,k will also be zero.
Moreover, L w j,k (kWh) is the battery charge after recharging, L max (kWh) is the battery size, and L 0,k (kWh) is the initial state of battery charge. ∆ is the node-link incidence matrix, E w k is a vector with Energies 2020, 13, 4061 6 of 14 a length of |N| and consists of two nonzero components: one has a value of 1, which corresponds to the original location; the other has a value of −1 which corresponds to the destination. x w a,k is a binary variable, and if link a is used, then x w a,k = 1, otherwise x w a,k = 0; y w i,k is also a binary variable, and if charging station i is selected, then y w i,k = 1, otherwise, y w i,k = 0. ρ w a,k = 0 if link a is used by driver k and unrestriced otherwise. is the BEV's energy consumption, m is the minimum comfortable range, and H and M are sufficiently large constants.
In the above SP-PT model, the objective function is to minimize the total trip time, including driving time (i.e., k∈K a∈A t a f a x w a,k ) and BEV drivers' perceived recharging time (i.e., k∈K i∈N α k c 1 i y w i,k + c 2 i β k l w i,k . Constraint (16) ensures flow balance. Constraints (17) and (19) specify the relationship between the states of charge of BEV batteries at the starting and ending nodes of any utilized link. Constraint (18) ensures that the BEV driver will not fully deplete their battery on any utilized link, and the BEV driver prefers to keep the remaining battery range no less than a comfortable range.
Let . . . , x w a,k , . . . , y w i,k , . . . , l w i,k , . . . be the optimal solution to SP-PT for OD pair w ∈ W for each type of k ∈ K. By solving x w a,k , we can obtain the shortest usable path, i.e., r w k . Moreover, through solving y w i,k and l w i,k , we can obtain the shortest perception recharging time, i.e., k∈K i∈N α k c 1 i y w i,k + c 2 i β k l w i,k . Then, to solve NE-HU, we can use the following iterative procedures: Step 1: For each type of k ∈ K between OD pair w ∈ W solve SP-PT, i.e., . . . , f a , . . . = (. . . , 0, . . .).
Obtain the initial usable path set R w k = r w k and calculate the minimal actual recharging time s w r,k = c 1 i y w i,k + c 2 i β k l w i,k .
Step 2: Based on the initial usable path set R w k = r w k , solve the NE-HU model and obtain the optimal traffic flow distribution, i.e., . . . , f a,k , . . . and the Lagrange multiplier, i.e., . . . , λ w k , . . . . Step 3: Solve the SP-PT model. For each type of k ∈ K and each OD pair w ∈ W, if λ w k ≤ k∈K a∈A t a f a x w a,k + k∈K i∈N α k c 1 i y w i,k + c 2 i β k l w i,k , then stop and . . . , f a,k , . . . is the optimal equilibrium link flow distribution; otherwise, go to Step 1.

Numerical Example
To demonstrate the above mixed user equilibrium model, a numerical example will be shown in this section. We considered the following Nguyen-Dupius network [26] (see Figure 2). There are 13 nodes, 19 links, and four OD pairs in the Nguyen-Dupius network. We assume that the link distances are two times the link's free-flow travel times, and the link characteristics are shown in Table 1. According to the Nissan Leaf 2020 [27], in the numerical example, we set the battery capacity as L max = 40 kWh and the energy consumption rate as = 0.167 kWh/km. Since in this paper we assume the en-route charging stations only supply direct current charging facilities (fast charging facilities), then the fixed time c 1 i = 5 min the variable time c 2 i = 0.6 min/kWh if the current battery level is less than 70%, and c 2 i = 10 min/kWh if the current battery level is higher than 70%. Note that in this paper, it is assumed that in order to reach the destination as quickly as possible, BEV drivers do not fully charge their electric car at the en-route charging station.
As shown in Figure 2, in the network, the charge stations are located in the nodes 5 and 10. It is assumed that in the network, 30% of vehicle drivers are BEV drivers, and the rest are gasoline vehicle drivers. Because of the different risk attitudes towards en-route recharging, BEV drivers are further divided into three different types and the proportion of each type is 10%. The OD demands, the number of each type of driver, and the parameter setting of risk attitude β k and perception error α k are shown in Table 2. Moreover, the setting of the battery safety margin and the initial state of the battery are shown in Table 3. It is worth mentioning that when the battery safety margin is set to 0.10 L max , then according to the setting in this paper, if the initial state of the battery is equal to the battery safety margin, i.e., 0.10 L max , then the battery electric vehicle can have a range of about 24km in theory. In general, there are 27 classes of BEV drivers (3 × 3 × 3 = 27) with different risk attitudes, Energies 2020, 13, 4061 7 of 14 different perception errors, and different initial states of the battery used in the numerical example. The above settings ensure that all BEV drivers can reach their destinations.
13 nodes, 19 links, and four OD pairs in the Nguyen-Dupius network. We assume that the link distances are two times the link's free-flow travel times, and the link characteristics are shown in Table 1. According to the Nissan Leaf 2020 [27], in the numerical example, we set the battery capacity as = 40 kWh and the energy consumption rate as = 0.167 kWh/km. Since in this paper we assume the en-route charging stations only supply direct current charging facilities (fast charging facilities), then the fixed time = 5 min the variable time = 0.6 min/kWh if the current battery level is less than 70%, and = 10 min/kWh if the current battery level is higher than 70%. Note that in this paper, it is assumed that in order to reach the destination as quickly as possible, BEV drivers do not fully charge their electric car at the en-route charging station.    1  14  7  3000  11  18  9  5000  2  18  9  2000  12  20  10  5500  3  18  9  2000  13  18  9  2000  4  24  12  2000  14  12  6  4000  5  6  3  3500  15  18  9  3000  6  18  9  4000  16  16  8  3000  7  10  7  5000  17  14  7  2000  8  26  13  2500  18  20  10  5000  9  10  5  2500  19  18  9  2000  10  18  9 2000 Table 2. The setting of parameters β k and α k .   Table 4 lists the detailed path information under the long-term equilibrium, including path flows and actual path travel times. It is worth noting that according to the above setting, all BEV drivers need to rechange their BEVs at one en-route charging station to keep the remaining battery no less than their battery safety margin. From Table 4, we can see that under the equilibrium, there are 13 available paths in total. However, due to the fact that path 1, path 6, and path 13 do not have en-route charging stations, these paths are only available for gasoline vehicles. Moreover, limited by travel distance and the location of charging stations, path 3, path 5, path 9, path 10, and path 12 are not available for battery electric vehicles whose initial state of batteries are low. Therefore, under the equilibrium, for OD pair (1,11), path 2's travel time is higher than the other two paths' travel times, and for OD pair (1,13), path 4's travel time is higher than the other two paths' travel times.  Figure 3 shows the average utilization of charging station 5 and charging station 10 under the long-term equilibrium. It can be observed from Figure 3 that with increased risk aversion and the initial state of battery, charging station 10 is chosen by more BEV drivers. This is because as the level of risk aversion increases, so do the battery safety margins and initial states of batteries of such BEV drivers, as shown in Table 3. Figure 4 shows the recharging time information of charging station 5 and charging station 10. It can be observed from Figure 4 that the average recharging time of charging station 5 is longer than the average recharging time of charging station 10. This is because compared with charging station 5, charging station 10 is closer to the destination, and to ensure the remaining battery is higher than the battery safety margin, BEV drivers need to take more time to recharge their electric vehicles if they choose charging station 5. In addition, as the level of risk aversion increases, so does the recharging time of such BEV drivers, and this will also add to the feeling of more time spent recharging.

Origin-Destination (OD) Pair
Energies 2020, 13, 4061 9 of 14 station 5 is longer than the average recharging time of charging station 10. This is because compared with charging station 5, charging station 10 is closer to the destination, and to ensure the remaining battery is higher than the battery safety margin, BEV drivers need to take more time to recharge their electric vehicles if they choose charging station 5. In addition, as the level of risk aversion increases, so does the recharging time of such BEV drivers, and this will also add to the feeling of more time spent recharging.   To further characterize the impact of the battery safety margin (i.e., a BEV driver's risk attitude) on the mixed user equilibrium, we gradually reduce the values of BEV drivers' battery safety margins from the initial values of the state to zero. The corresponding results of the average utilization of charging station 5 and charging station 10 are shown in Figure 5.  To further characterize the impact of the battery safety margin (i.e., a BEV driver's risk attitude) on the mixed user equilibrium, we gradually reduce the values of BEV drivers' battery safety margins from the initial values of the state to zero. The corresponding results of the average utilization of charging station 5 and charging station 10 are shown in Figure 5.
It can be observed from Figure 5 that when β = 1.1 (i.e., when there is a lower risk aversion), as the battery safety margins decrease, there are only a few BEV drivers switching from charging station 5 to charging station 10, and when the value of the battery safety margin equals zero, there are very few BEV drivers choosing not to recharge their electric vehicles. Notably, these BEV drivers change their en-route recharging behaviors because their initial state of batteries are higher than others. When β = 1.3 (i.e., when there is a medium risk aversion), as the battery safety margins decrease, a growing number of BEV drivers switch to recharging station 10, and there are close to 350 BEV drivers who choose not to recharge their electric vehicles. When β = 1.5 (i.e., when there is a higher risk aversion), as the battery safety margins decrease, a lot of BEV drivers opt for not recharging their electric vehicles. Moreover, according to the Table 4, we can find that when some BEV drivers choose not to recharge their electric vehicles en-route, path 2 and path 4 are not the best choice. Thus, Table 5 lists the detailed path information under the equilibrium without the battery safety margin. Compared with Table 4, it can be observed that in Table 5, some paths' travel times are increased and some paths' travel times are decreased. However, the total travel time is 851,318 min in Table 4, and 907,994 min in Table 5. Figure 6 shows the recharging time information of charging station 5 and charging station 10 without the battery safety margin. Compared with Figure 4, it can be observed that the recharging times become smaller for charging station 5 and charging station 10 when β = 1.1, but the recharging times become larger for charging station 10 when β = 1.3 and β = 1.5. This is because more and more BEV drivers choose to recharge their electric vehicles at charging station 10 when they do not consider their battery safety margin. To further characterize the impact of the battery safety margin (i.e., a BEV driver's risk attitude) on the mixed user equilibrium, we gradually reduce the values of BEV drivers' battery safety margins from the initial values of the state to zero. The corresponding results of the average utilization of charging station 5 and charging station 10 are shown in Figure 5. It can be observed from Figure 5 that when = 1.1 (i.e., when there is a lower risk aversion), as the battery safety margins decrease, there are only a few BEV drivers switching from charging station 5 to charging station 10, and when the value of the battery safety margin equals zero, there are very few BEV drivers choosing not to recharge their electric vehicles. Notably, these BEV drivers change their en-route recharging behaviors because their initial state of batteries are higher than others. When = 1.3 (i.e., when there is a medium risk aversion), as the battery safety margins

Conclusions and Future Work
By considering electric vehicle drivers' battery management strategies, this paper investigated mixed user equilibrium problems with battery electric vehicle drivers' en-route recharging behaviors. Considering two kinds of vehicle drivers (i.e., BEV drivers and gasoline vehicle drivers) in the network, we assume that BEV drivers select paths with the shortest driving time and shortest perceived recharging time, and gasoline vehicle drivers select paths with the shortest driving travel time. Moreover, we also assume that the BEVs will not be fully charged at the en-route charging station, and in order to reach the destination as quickly as possible, BEV drivers will recharge a small amount of electricity and keep the remaining battery no less than the battery safety margin when they reach their destinations. We proposed a mixed user equilibrium model to describe the resulting network equilibrium flow distributions. A numerical example was used to demonstrate the proposed mixed user equilibrium model. It was observed that affected by the battery safety margin, under the equilibrium, for some OD pairs, some paths' travel times were higher than other available paths. Additionally, with increased risk aversion, charging station 10 was chosen by more BEV drivers, which may add to the feeling of more time spent recharging. Furthermore, in the numerical example, we further characterized the impact of the battery safety margin on the mixed user equilibrium, and provided the mixed user equilibrium information without considering the battery safety margin. It was observed that BEV drivers' en-route recharging choice behaviors are influenced by their battery safety margins, and the total travel time was 907,994 min when not considering BEV drivers' battery safety margins, and 851,318 min when considering BEV drivers' battery safety margins at the outset.
Note that the scale of BEV drivers and BEV drivers' risk attitudes, battery safety margins, perception errors, and the initial states of batteries are imaginary. Thus, our future study will attempt to investigate BEV drivers' risk attitudes, battery safety margins, and en-route recharging choice

Conclusions and Future Work
By considering electric vehicle drivers' battery management strategies, this paper investigated mixed user equilibrium problems with battery electric vehicle drivers' en-route recharging behaviors. Considering two kinds of vehicle drivers (i.e., BEV drivers and gasoline vehicle drivers) in the network, we assume that BEV drivers select paths with the shortest driving time and shortest perceived recharging time, and gasoline vehicle drivers select paths with the shortest driving travel time. Moreover, we also assume that the BEVs will not be fully charged at the en-route charging station, and in order to reach the destination as quickly as possible, BEV drivers will recharge a small amount of electricity and keep the remaining battery no less than the battery safety margin when they reach their destinations. We proposed a mixed user equilibrium model to describe the resulting network equilibrium flow distributions. A numerical example was used to demonstrate the proposed mixed user equilibrium model. It was observed that affected by the battery safety margin, under the equilibrium, for some OD pairs, some paths' travel times were higher than other available paths. Additionally, with increased risk aversion, charging station 10 was chosen by more BEV drivers, which may add to the feeling of more time spent recharging. Furthermore, in the numerical example, we further characterized the impact of the battery safety margin on the mixed user equilibrium, and provided the mixed user equilibrium information without considering the battery safety margin. It was observed that BEV drivers' en-route recharging choice behaviors are influenced by their battery safety margins, and the total travel time