Overview#
Since version 0.24.0 SUMO includes a model for electric vehicles. It was implemented by Tamás Kurczveil and Pablo Álvarez López from the TU-Braunschweig. The core of the model is implemented in the vehicle device device.battery. Additional features are a charging station (which can be placed on any lane in the network) and a new output option --battery-output <FILE>.
You can find a test case for these implementations at [1]
Defining Electric Vehicles#
Different aspects of electric vehicles are modeled separately. This page puts the focus on modeling the battery and how it is charged and discharged. The actual energy consumption values themselves are part of the emission modelling because SUMO can use different models for that.
To track the charging status of a vehicle, it must be equipped with a battery device. This may be done using the option --device.battery.explicit <vehID1,vehID2,...> or simply setting --device.battery.probability 1 to equip all vehicles. Alternatively, the device may be specified using Generic vehicle parameters.
Additional properties of the vehicle and its electrical components must then be defined via parameters of the vehicle or its type. Some property can only be defined for the vehicle type.
These values have the following meanings (the defaults are from the Kia below):
key | Value Type | Default | Description |
---|---|---|---|
device.battery.capacity | float | 35000 (Wh) | Maximum battery capacity Emax |
maximumPower | float | 150000 (W) | Maximum power which the vehicle can reach (unused) |
vehicleMass | float | 1830 (kg) | Vehicle mass mveh (deprecated) |
loading | float | 0 (kg) | Additional mass (to be defined in the vehicle type) |
frontSurfaceArea | float | 2.6 (m2) | Front surface area Aveh |
airDragCoefficient | float | 0.35 | Air drag coefficient cw |
rotatingMass | float | 40 (kg) | (Equivalent) mass of internal rotating elements |
radialDragCoefficient | float | 0.1 | Radial drag coefficient crad |
rollDragCoefficient | float | 0.01 | Rolling resistance coefficient croll |
constantPowerIntake | float | 100 (W) | Avg. (constant) power of consumers Pconst |
propulsionEfficiency | float | 0.98 | Drive efficiency ηprop |
recuperationEfficiency | float | 0.96 | Recuperation efficiency ηrecup |
stoppingThreshold | float | 0.1 (m/s) | Maximum velocity to start charging |
device.battery.maximumChargeRate | float | 150000 (W) | Maximum charging rate of the battery |
device.battery.chargeLevelTable | float list | Ordered list of state of charge values (from 0 to 1) for which maximum charge rates are defined in device.battery.chargeCurveTable |
|
device.battery.chargeCurveTable | float list | Corresponding maximum charge rates to each state of charge value in device.battery.chargeLevelTable |
Note
Before SUMO 1.20.0 the rotatingMass
was called internalMomentOfInertia
but it has been renamed to make clear
that it is a mass and not a moment of inertia. The old parameter is considered deprecated.
Also the vehicleMass
has been deprecated in favor of the new mass
attribute of vehicles / vehicle types.
An example of a vehicle with electric attribute (those are the values for a city bus from the original publication):
<routes>
<vType id="ElectricBus" accel="1.0" decel="1.0" length="12" maxSpeed="100.0" sigma="0.0" minGap="2.5" mass="10000" color="1,1,1">
<param key="has.battery.device" value="true"/>
<param key="device.battery.capacity" value="2000"/>
<param key="maximumPower" value="1000"/>
<param key="frontSurfaceArea" value="5"/>
<param key="airDragCoefficient" value="0.6"/>
<param key="rotatingMass" value="100"/>
<param key="radialDragCoefficient" value="0.5"/>
<param key="rollDragCoefficient" value="0.01"/>
<param key="constantPowerIntake" value="100"/>
<param key="propulsionEfficiency" value="0.9"/>
<param key="recuperationEfficiency" value="0.9"/>
<param key="stoppingThreshold" value="0.1"/>
<param key="device.battery.maximumChargeRate" value="150000"/>
</vType>
</routes>
If a vehicle has a battery device (and is not tracking fuel)
and no explicit emissionClass
is defined, it will be assigned the emission class Energy/unknown
.
It will not use the default emission class derived from the vehicle class then.
This is for backward compatibility and will issue a warning because in general it is preferable to set the emission class explicitly.
Most of the parameters above do actually apply to this emission class and not to the battery device itself.
The initial energy content of the battery (by default 0.5*device.battery.capacity
) can
be set in the vehicle definitions
<routes>
<vehicle id="0" type="type1" depart="0" color="1,0,0">
<param key="device.battery.chargeLevel" value="500"/>
</vehicle>
</routes>
The charging rate of the battery at a charging station is limited to model the effects of battery management controllers (e.g. charge a nearly full battery less than an nearly empty one).
There are two ways to define the maximum charge rate: For a constant rate set the attribute device.battery.maximumChargeRate
.
If instead a maximum charge rate depending on the state of charge is wanted, it can be defined through data points between which the rate will get interpolated.
The states of charge have to be given in device.battery.chargeLevelTable
and the corresponding charge rates in device.battery.chargeCurveTable
.
If defined, the maximum charge curve takes precedence over the constant maximum charge rate device.battery.maximumChargeRate
. An example definition where the charge rate decreases above
50% state of charge looks like the following:
<routes>
<vehicle id="0" type="type1" depart="0" color="1,0,0">
<param key="device.battery.chargeLevelTable" value="0 0.5 1"/>
<param key="device.battery.chargeCurveTable" value="45000 45000 20000"/>
</vehicle>
</routes>
Vehicle behavior#
By default, vehicle behavior will not be affected by battery level. Car will keep driving even when their battery capacity is at 0. To avoid this, either TraCI must be used to change speed or route based on the current battery level or the stationfinder device can be configured to monitor the battery capacity.
Charging Stations#
A charging station is a surface defined on a lane in which the vehicles equipped with a battery are charged. The basic structure and parameters of bus stops were used for the implementation of charging stations.
key | Value Type | Value range | Default | Description |
---|---|---|---|---|
id | string | id | Charging station ID (Must be unique) | |
name | string | simple String | Charging station name. This is only used for visualization purposes. | |
lane | string | valid lane id | Lane of the charging station location | |
startPos | float | lane.length < x < lane.length (negative values count backwards from the end of the lane) | 0 | Begin position in the specified lane |
endPos | float | lane.length < x < lane.length (negative values count backwards from the end of the lane) | End position in the specified lane | |
friendlyPos | bool | true or false | false | Whether invalid charging station positions should be corrected automatically |
power | float (W) or (mg/s) | power > 0 | 22000 | Charging power Pchrg (If the battery device being charged is configured to track fuel, charging power will be interpreted as mg/s) |
efficiency | float | 0 <= efficiency <= 1 | 0.95 | Charging efficiency ηchrg |
chargeInTransit | bool | true or false | false | Enable or disable charge in transit, i.e. vehicle is forced/not forced to stop for charging |
chargeDelay | float | chargeDelay > 0 | 0 | Time delay after the vehicles have reached / stopped on the charging station, before the energy transfer (charging) is starting |
chargeType | string | normal | Charging type (normal, electric, fuel) | |
parkingArea | string | valid parkingArea id | id of the parking the charging station should be positioned on (optional) - vehicles will only charge after reaching the parking |
Charging stations are defined in additional using the following format:
<additional>
<chargingStation chargeDelay="2" chargeInTransit="0" power="200000" efficiency="0.95" endPos="25" id="cS_2to19_0a" lane="2to19_0" startPos="10"/>
</additional>
And are represented in the simulation as shown:
Representation of chargingStation in GUI
Color of chargingStation during charge
Stopping at a Charging Station#
A stop at a charging station may occur due to traffic conditions, stopping at a defined location or stopping at an explicit chargingStation as defined below:
<routes>
<vehicle id="v0" route="route0" depart="0">
<stop chargingStation="myChargingStationID" until="50"/>
</vehicle>
</routes>
Charging Station output#
There are two variants of the charging station output. One lists values for all time steps whereas the other aggregates the data into charging events (thus related to the whole time a vehicle was connected to the charging station).
Full report#
Option --chargingstations-output chargingstations.xml generates a full report of energy charged by charging stations:
<output>
<battery-output value="battery.xml"/>
<chargingstations-output value="chargingstations.xml"/>
</output>
File chargingstations.xml has the following structure:
<chargingstations-export>
<chargingStation id="CS1" totalEnergyCharged="71.25" chargingSteps="27">
<vehicle id="veh0" type="ElectricVehicle1" totalEnergyChargedIntoVehicle="15.83" chargingBegin="12.00" chargingEnd="17.00">
<step time="12.00" chargingStatus="chargingStopped" energyCharged="2.64" partialCharge="2.64" power="10000.00" efficiency="0.95" actualBatteryCapacity="12962.97" maximumBatteryCapacity="35000.00"/>
<step time="13.00" chargingStatus="chargingStopped" energyCharged="2.64" partialCharge="5.28" power="10000.00" efficiency="0.95" actualBatteryCapacity="12965.59" maximumBatteryCapacity="35000.00"/>
<step time="14.00" chargingStatus="chargingStopped" energyCharged="2.64" partialCharge="7.92" power="10000.00" efficiency="0.95" actualBatteryCapacity="12968.22" maximumBatteryCapacity="35000.00"/>
</vehicle>
<vehicle id="veh1" type="ElectricVehicle2" totalEnergyChargedIntoVehicle="5.28" chargingBegin="17.00" chargingEnd="18.00">
<step time="17.00" chargingStatus="chargingStopped" energyCharged="2.64" partialCharge="18.47" power="10000.00" efficiency="0.95" actualBatteryCapacity="11967.35" maximumBatteryCapacity="35000.00"/>
<step time="18.00" chargingStatus="chargingStopped" energyCharged="2.64" partialCharge="21.11" power="10000.00" efficiency="0.95" actualBatteryCapacity="12978.72" maximumBatteryCapacity="35000.00"/>
</vehicle>
...
</chargingStation>
...
...
</chargingstations-export>
For the entire ChargingStation:
Name | Type | Description |
---|---|---|
id | string | ChargingStation ID |
totalEnergyCharged | float | Total energy charged in Wh during the entire simulation |
chargingSteps | int | Number of steps in which chargingStation was charging energy |
For the current charging vehicle
Name | Type | Description |
---|---|---|
id | string | ID of vehicle that is charging in these charging stop |
type | string | type of vehicle |
totalEnergyChargedIntoVehicle | double | Energy charged (in Wh) during these charging stop |
chargingBegin | float | TimeStep in which vehicle starts with the charge (in seconds) |
chargingEnd | float | TimeStep in which vehicle ends with the charge (in seconds) |
For every charging timeStep:
Name | Type | Description |
---|---|---|
time | float | Current timestep (s) |
chargingStatus | string | Current charging status (Charging, waiting to charge o not charging) |
energyCharged | float | Energy charged in current timeStep |
partialCharge | float | Energy charged by ChargingStation from begin of simulation to this timeStep |
power | float | Current power of ChargingStation |
efficiency | float | Current efficiency of ChargingStation |
actualBatteryCapacity | string | Current battery capacity of vehicle |
maximumBatteryCapacity | string | Current maximum battery capacity of vehicle |
Charging events#
Option --chargingstations-output chargingevents.xml --chargingstations-output.aggregated true generates a report of the charging events at any charging station in the network. The generated output file chargingevents.xml has the following structure:
<chargingstations-export>
<chargingEvent chargingStation="CS1" vehicle="veh0" type="ElectricVehicle1" totalEnergyChargedIntoVehicle="15.83" chargingBegin="12.00" chargingEnd="17.00" actualBatteryCapacity="12968.22" maximumBatteryCapacity="35000.00" minPower="10000.00" maxPower="10000.00" minCharge="2.64" maxCharge="2.64" minEfficiency="0.95" maxEfficiency="0.95" />
<chargingEvent chargingStation="CS1" vehicle="veh1" type="ElectricVehicle1" totalEnergyChargedIntoVehicle="5.28" chargingBegin="17.00" chargingEnd="18.00" actualBatteryCapacity="12978.72" maximumBatteryCapacity="35000.00" minPower="10000.00" maxPower="10000.00" minCharge="2.64" maxCharge="2.64" minEfficiency="0.95" maxEfficiency="0.95" />
...
</chargingstations-export>
If the option --chargingstations-output.aggregated.write-unfinished true is set as well, vehicles still charging at the end of the simulation are included as well (but do not contain the chargingEnd attribute). Attributes with the same name can be looked up in the table above. The remaining attributes are:
Name | Type | Description |
---|---|---|
minPower | float | Minimum charging power during the charging event |
maxPower | float | Maximum charging power during the charging event |
minCharge | float | Minimum charged energy in one time step of the charging event |
maxCharge | float | Maximum charged energy in one time step of the charging event |
minEfficiency | float | Minimum charging efficiency during the charging event |
maxEfficiency | float | Maximum charging efficiency during the charging event |
battery-output#
There are three output parameters to be set in the SUMO configuration to use the battery device:
<configuration>
<input>
<net-file value="netFile.xml"/>
<route-files value="routeFile.xml"/>
<additional-files value="additionalFile.xml"/>
</input>
<output>
<battery-output value="Battery.out.xml"/>
<battery-output.precision value="4"/>
<device.battery.probability value="1"/>
<summary-output value="summary_100.xml"/>
</output>
</configuration>
The battery-output generates a file with this structure:
<battery-export>
<timestep time="0.00">
<vehicle id="vehicle01" energyConsumed="0.00" totalEnergyConsumed="0.00" totalEnergyRegenerated="0.00"
actualBatteryCapacity="17500.00" maximumBatteryCapacity="35000.00"
chargingStationId="NULL" energyCharged="0.00" energyChargedInTransit="0.00" energyChargedStopped="0.00"
speed="12.92" acceleration="0.00" x="1428.27" y="25.57" lane="01to02_0"
posOnLane="0.00" timeStopped="0"/>
<vehicle id=..... */
</timestep>
<timestep time="1.00">
<vehicle id=.....
</timestep>
<timestep time=...
...
</timestep>
</battery-export>
Name | Type | Description |
---|---|---|
time | int | Current timestep |
id | string | id of vehicle |
energyConsumed | double | energy consumption in this timestep in Wh |
totalEnergyConsumed | double | cumulative sum of energy consumption up to this timestep in Wh |
totalEnergyRegenerated | double | cumulative sum of regenerated energy up to this timestep in Wh |
actualBatteryCapacity | double | energy content of the battery in this timestep |
maximumBatteryCapacity | double | Max energy capacity of the battery |
chargingStationId | string | If vehicle is exactly at a charging station, this value is the id of the charging station, in other case, is NULL |
energyCharged | double | Charge received in the current time step from a charging station (Only != 0 if vehicle is exactly at a charging station) |
energyChargedInTransit | double | Charge that a vehicle in transit received in the current time step from a charging station |
energyChargedStopped | double | Charge that a stopped vehicle received in the current time step from a charging station |
speed | double | Speed of vehicle in this timestep |
acceleration | double | Acceleration of vehicle in this timestep |
x | double | absolute position x of vehicle in the map |
y | double | absolute position y of vehicle in the map |
lane | string | id of the lane that the vehicle is currently on |
posOnLane | double | Position of vehicle on its current lane |
timeStopped | int | Counter with the number of timesteps that the vehicle has remained standing |
Emission Output#
The Emission model-outputs of
SUMO can be used together with the battery device when setting the
<vType>
-parameter emissionClass="Energy/unknown"
.
Tracking fuel consumption for non-electrical vehicles#
By setting option --device.battery.track-fuel, equipped vehicles with a conventional drive train (emissionClass other than Energy
) will monitor their fuel level based on the fuel consumption of their respective emission class. All capacity values are then interpreted as mg instead of Wh. Also, the chargingStation power is re-interpreted as mg/s when charging fuel.
TraCI#
The internal state of the battery device can be accessed directly using traci.vehicle.getParameter and traci.vehicle.setParameter. Charging stations can be inspected and updated using the respective getter and setter functions inside traci.chargingstation.
Furthermore, the function
traci.vehicle.getElectricityConsumption()
can be used to access the consumption of the vehicle if the emissionClass="Energy/unknown"
is
declared.
Query of the current state of charge#
The current state of charge of the battery can be computed as the quotient of the actual battery charge and the maximum battery capacity:
capacity = float(traci.vehicle.getParameter(vehID, "device.battery.capacity"))
currentCharge = float(traci.vehicle.getParameter(vehID, "device.battery.chargeLevel"))
stateOfCharge = currentCharge / capacity
Calculating the remaining Range#
After the vehicle has been driving for a while, the remaining range can be computed based on previous consumption and distance:
mWh = traci.vehicle.getDistance(vehID) / float(traci.vehicle.getParameter(vehID, "device.battery.totalEnergyConsumed"))
remainingRange = float(traci.vehicle.getParameter(vehID, "device.battery.chargeLevel")) * mWh
To compute the remaining range on departure, the value of mWh
(meters per Watt-hour) should be calibrated from a prior simulation.
Reducing the power of a charging station#
If too many consumers connect to the eletrical grid, it may not be able to supply the nominal power of the charging station. A temporary power drop of 20% can be modeled using the following sample code:
prevPower = traci.chargingstation.getPower(chargingStationID) # remember for restoring the full power later
traci.chargingstation.setPower(chargingStationID, prevPower * 0.8)
Model Details#
All information about the implemented device (including details on the vehicle energy consumption and charging model) can be found in the following publication.
Publications#
Example Configurations#
Kia Soul EV 2020#
The values are provided by courtesy of Jim Div based on his own calibration.
<!-- vehicle mass = 1682kg curb wt + average 2 passengers / bags -->
<vType id="soulEV65" minGap="2.50" maxSpeed="29.06" color="white" accel="1.0" decel="1.0" sigma="0.0" emissionClass="Energy/unknown" mass="1830">
<param key="has.battery.device" value="true"/>
<param key="airDragCoefficient" value="0.35"/> <!-- https://www.evspecifications.com/en/model/e94fa0 -->
<param key="constantPowerIntake" value="100"/> <!-- observed summer levels -->
<param key="frontSurfaceArea" value="2.6"/> <!-- computed (ht-clearance) * width -->
<param key="rotatingMass" value="40"/> <!-- guesstimate, inspired by PHEMlight5 PC_BEV -->
<param key="device.battery.capacity" value="64000"/>
<param key="maximumPower" value="150000"/> <!-- website as above -->
<param key="propulsionEfficiency" value=".98"/> <!-- guesstimate value providing closest match to observed -->
<param key="radialDragCoefficient" value="0.1"/> <!-- as above -->
<param key="recuperationEfficiency" value=".96"/> <!-- as above -->
<param key="rollDragCoefficient" value="0.01"/> <!-- as above -->
<param key="stoppingThreshold" value="0.1"/> <!-- as above -->
</vType>
Observations: - Simulation efficiencies of 6.3 - 6.7 km driven per kWh consumed agree with measured efficiencies of 6.4 - 6.8 (stddev 0.4 - 0.8) in short and medium range simulations with realistic traffic - abstract scenarios without junctions and other cars overestimate the efficiency by a large factor (~ twice as many km/kWh)