#!/usr/bin/env python
# -*- coding: UTF-8 -*-
#
# Copyright 2015-2022 European Commission (JRC);
# Licensed under the EUPL (the 'Licence');
# You may not use this work except in compliance with the Licence.
# You may obtain a copy of the Licence at: http://ec.europa.eu/idabc/eupl
"""
Constants for the CO2MPAS physical model.
"""
import schedula as sh
import co2mpas.utils as co2_utl
#: Container of node default values.
# noinspection PyMissingOrEmptyDocstring,PyPep8Naming
[docs]class Values(co2_utl.Constants):
#: Minimum distance between points for the elevation interpolation [m].
minimum_elevation_distance = 30
#: Is the vehicle plugin hybrid?
is_plugin = False
#: Is the vehicle serial hybrid?
is_serial = False
#: Apply RCB correction?
rcb_correction = True
#: Family correction factor for representative regional temperatures [-].
atct_family_correction_factor = 1.0
#: Drive battery technology type.
drive_battery_technology = 'unknown'
#: Belt efficiency [-].
belt_efficiency = 0.8
#: Starter efficiency [-].
starter_efficiency = 0.7
#: Tyre state (i.e., new or worm).
tyre_state = 'new'
#: Road state (i.e., dry, wet, rainfall, puddles, ice).
road_state = 'dry'
#: Number of engine cylinders [-].
engine_n_cylinders = 4
#: Default additive correction for vehicles with periodically regenerating
#: systems [CO2g/km].
ki_additive = 0
#: Number of passengers including driver [-].
n_passengers = 1
#: Average passenger mass [kg].
passenger_mass = 75
#: Cargo mass [kg].
cargo_mass = 0
#: Fuel mass [kg].
fuel_mass = 25
#: Does the engine have selective catalytic reduction technology?
has_selective_catalytic_reduction = False
#: Does the engine have lean burn technology?
has_lean_burn = False
#: Does the gear box have some additional technology to heat up faster?
has_gear_box_thermal_management = False
#: Does the vehicle has periodically regenerating systems? [-].
has_periodically_regenerating_systems = False
#: Possible percentages of active cylinders [-].
active_cylinder_ratios = (1.0,)
#: Does the engine feature variable valve actuation? [-].
engine_has_variable_valve_actuation = False
#: NEDC cycle time [s].
max_time_NEDC = 1180.0
#: WLTP cycle time [s].
max_time_WLTP = 1800.0
#: Maximum velocity to consider the vehicle stopped [km/h].
stop_velocity = 1.0 + 1.1920929e-07
#: Maximum acceleration to be at constant velocity [m/s2].
plateau_acceleration = 0.1 + 1.1920929e-07
#: Does the vehicle have start/stop system?
has_start_stop = True
#: Does the engine have cylinder deactivation technology?
engine_has_cylinder_deactivation = False
#: Minimum vehicle engine speed [RPM].
min_engine_on_speed = 10.0
#: Minimum time of engine on after a start [s].
min_time_engine_on_after_start = 4.0
#: Time window used to apply gear change filters [s].
change_gear_window_width = 4.0
#: Service battery start window width [s].
service_battery_start_window_width = 4.0
#: Threshold vehicle velocity for gear correction due to full load curve
#: [km/h].
max_velocity_full_load_correction = 100.0
#: Air temperature [°C].
air_temperature = 20
#: Atmospheric pressure [kPa].
atmospheric_pressure = 101.325
#: Has the vehicle a roof box? [-].
has_roof_box = False
#: Tyre class (C1, C2, and C3).
tyre_class = 'C1'
#: Angle slope [rad].
angle_slope = 0.0
#: A different preconditioning cycle was used for WLTP and NEDC?
correct_f0 = False
#: Final drive ratio [-].
final_drive_ratio = 1.0
#: Wheel drive (i.e., front, rear, front+rear).
wheel_drive = 'front'
#: Apply the eco-mode gear shifting?
fuel_saving_at_strategy = True
#: Cold and hot gear box reference temperatures [°C].
gear_box_temperature_references = (40.0, 80.0)
#: Constant torque loss factors due to engine auxiliaries [N/cc, N*m].
auxiliaries_torque_loss_factors = (0.175, 0.2021) # m, q
#: Constant power loss due to engine auxiliaries [kW].
auxiliaries_power_loss = 0.0213
#: If the engine is equipped with any kind of charging.
engine_is_turbo = True
#: Start-stop activation time threshold [s].
start_stop_activation_time = 30
#: Standard deviation of idle engine speed [RPM].
idle_engine_speed_std = 100.0
#: Is an hot cycle?
is_cycle_hot = False
#: CO2 emission model params.
co2_params = {}
#: Enable the calculation of Willans coefficients for all phases?
enable_phases_willans = False
#: Enable the calculation of Willans coefficients for the cycle?
enable_willans = False
#: Alternator efficiency [-].
alternator_efficiency = 0.67
#: Time elapsed to turn on the engine with electric starter [s].
delta_time_engine_starter = .5
#: If to use decision tree classifiers to predict gears.
use_dt_gear_shifting = False
#: Does the vehicle have energy recuperation features?
has_energy_recuperation = True
#: A/T Time at cold hot transition phase [s].
time_cold_hot_transition = 300.0
#: Time frequency [1/s].
time_sample_frequency = 1.0
#: Initial temperature of the test cell of NEDC [°C].
initial_temperature_NEDC = 25.0
#: Initial temperature of the test cell of WLTP [°C].
initial_temperature_WLTP = 23.0
#: K1 NEDC parameter (first or second gear) [-].
k1 = 1
#: K2 NEDC parameter (first or second gear) [-].
k2 = 2
#: K5 NEDC parameter (first or second gear) [-].
k5 = 2
#: WLTP base model params.
wltp_base_model = {}
#: Velocity downscale factor threshold [-].
downscale_factor_threshold = 0.01
#: Empirical value in case of CVT [-].
tyre_dynamic_rolling_coefficient = 3.05 / 3.14
#: Container of internal function parameters.
# noinspection PyMissingOrEmptyDocstring,PyPep8Naming
[docs]class Functions(co2_utl.Constants):
ENABLE_ALL_FUNCTIONS = False
[docs] class default_after_treatment_warm_up_duration(co2_utl.Constants):
#: After treatment warm up duration for conventional vehicles [s].
duration = 60
[docs] class default_after_treatment_cooling_duration(co2_utl.Constants):
#: After treatment cooling duration for conventional vehicles [s].
duration = float('inf')
# noinspection PyMissingOrEmptyDocstring,PyPep8Naming
[docs] class parse_solution(co2_utl.Constants):
#: Use all calibration outputs as relative prediction targets.
CALIBRATION_AS_TARGETS = False
# noinspection PyMissingOrEmptyDocstring,PyPep8Naming
[docs] class select_best_model(co2_utl.Constants):
#: Model selector mapping.
MAP = {None: {
'wltp_h': ('nedc_h', 'wltp_h'),
'wltp_l': ('nedc_l', 'wltp_l'),
'wltp_m': ('wltp_m',)
}}
# noinspection PyMissingOrEmptyDocstring,PyPep8Naming
[docs] class calculate_cylinder_deactivation_valid_phases(co2_utl.Constants):
#: Engine inertia [kW].
LIMIT = 0.0001
# noinspection PyMissingOrEmptyDocstring,PyPep8Naming
[docs] class calculate_drive_battery_n_parallel_cells_v1(co2_utl.Constants):
#: Voltage for calculating the number of parallel battery cells [V].
reference_volt = {
'unknown': 2.9,
'NiMH': 1.1,
'Li-NCA (Li-Ni-Co-Al)': 2.9,
'Li-NCM (Li-Ni-Mn-Co)': 2.9,
'Li-MO (Li-Mn)': 2.9,
'Li-FP (Li-Fe-P)': 2.4,
'Li-TO (Li-Ti)': 1.7
} # source: batteryuniversity.com
# noinspection PyMissingOrEmptyDocstring,PyPep8Naming
[docs] class default_planetary_mean_efficiency(co2_utl.Constants):
#: Planetary mean efficiency [-].
efficiency = .97
# noinspection PyMissingOrEmptyDocstring,PyPep8Naming
[docs] class default_planetary_ratio(co2_utl.Constants):
#: Fundamental planetary speed ratio [-].
ratio = 2.6
# noinspection PyMissingOrEmptyDocstring,PyPep8Naming
[docs] class identify_after_treatment_warm_up_phases(co2_utl.Constants):
#: After treatment cooling time [s].
cooling_time = 400
# noinspection PyMissingOrEmptyDocstring,PyPep8Naming
[docs] class define_service_battery_electric_powers_supply_threshold(
co2_utl.Constants):
#: Minimum SOC variation to define service battery charging status [%].
min_soc = 0.1
#: Maximum allowed negative current for the service battery being
#: considered not charging [A].
min_current = -1.0
# noinspection PyMissingOrEmptyDocstring,PyPep8Naming
[docs] class default_final_drive_efficiency(co2_utl.Constants):
#: Formula for the default final drive efficiency [function].
formula = '1 - (n_wheel_drive - 2) / 100'
# noinspection PyMissingOrEmptyDocstring,PyPep8Naming
[docs] class predict_clutch_tc_speeds_delta(co2_utl.Constants):
#: Enable prediction of clutch or torque converter speeds delta?
ENABLE = False
# noinspection PyMissingOrEmptyDocstring,PyPep8Naming
[docs] class default_tc_normalized_m1000_curve(co2_utl.Constants):
#: Normalized m1000 curve [-].
curve = dict(
x=[
0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.65, 0.7, 0.75, 0.8, 0.87,
0.9, 0.95, 1, 1.1, 1.222, 1.375, 1.571, 1.833, 2.2, 2.5, 3, 3.5,
4, 4.5, 5
],
y=[
1, 0.97, 0.93, 0.9, 0.87, 0.83, 0.8, 0.75, 0.7, 0.65, 0.6, 0.55,
0.5, 0.25, 0, -0.099, -0.198, -0.336, -0.535, -0.828, -1.306,
-1.781, -2.772, -4.071, -5.746, -7.861, -10.480
]
)
# noinspection PyMissingOrEmptyDocstring,PyPep8Naming
[docs] class default_n_wheel(co2_utl.Constants):
#: Total number of wheels [-].
n_wheel = 4
# noinspection PyMissingOrEmptyDocstring,PyPep8Naming
[docs] class default_static_friction(co2_utl.Constants):
#: Static friction coefficient [-].
coeff = dict(
new=dict(dry=.85, wet=.65, rainfall=.55, puddles=.5, ice=.1),
worm=dict(dry=1, wet=.5, rainfall=.4, puddles=.25, ice=.1)
)
# noinspection PyMissingOrEmptyDocstring,PyPep8Naming
[docs] class calculate_velocities(co2_utl.Constants):
#: Time window for the moving average of obd velocities [s].
dt_window = 5
# noinspection PyMissingOrEmptyDocstring,PyPep8Naming
[docs] class calculate_engine_temperature_derivatives(co2_utl.Constants):
#: Derivative spacing [s].
dx = 4
#: Degree of the smoothing spline [-].
order = 7
#: Time window for smoother [s].
tw = 20
# noinspection PyMissingOrEmptyDocstring,PyPep8Naming
[docs] class default_clutch_window(co2_utl.Constants):
#: Clutching time window [s].
clutch_window = (0, 0.95384615)
# noinspection PyMissingOrEmptyDocstring,PyPep8Naming
[docs] class define_fuel_type_and_is_hybrid(co2_utl.Constants):
#: Is the vehicle hybrid?
is_hybrid = {
# 0: False, # Not available
1: False, # Gasoline
2: False, # Methanol
3: False, # Ethanol
4: False, # Diesel
5: False, # LPG
6: False, # CNG
7: False, # Propane
8: True, # Electric
9: False, # Bifuel running Gasoline
10: False, # Bifuel running Methanol
11: False, # Bifuel running Ethanol
12: False, # Bifuel running LPG
13: False, # Bifuel running CNG
14: False, # Bifuel running Propane
15: True, # Bifuel running Electricity
16: True, # Bifuel running electric and combustion engine
17: True, # Hybrid gasoline
18: True, # Hybrid Ethanol
19: True, # Hybrid Diesel
20: True, # Hybrid Electric
21: True, # Hybrid running electric and combustion engine
22: True, # Hybrid Regenerative
23: False, # Bifuel running diesel
}
#: The vehicle fuel type.
fuel_type = {
# 0: 'gasoline', # Not available
1: 'gasoline', # Gasoline
2: 'methanol', # Methanol
3: 'ethanol', # Ethanol
4: 'diesel', # Diesel
5: 'LPG', # LPG
6: 'NG', # CNG
7: 'propane', # Propane
8: None, # Electric
9: 'gasoline', # Bifuel running Gasoline
10: 'methanol', # Bifuel running Methanol
11: 'ethanol', # Bifuel running Ethanol
12: 'LPG', # Bifuel running LPG
13: 'NG', # Bifuel running CNG
14: 'propane', # Bifuel running Propane
15: None, # Bifuel running Electricity
16: 'gasoline', # Bifuel running electric and combustion engine
17: 'gasoline', # Hybrid gasoline
18: 'ethanol', # Hybrid Ethanol
19: 'diesel', # Hybrid Diesel
20: None, # Hybrid Electric
21: 'gasoline', # Hybrid running electric and combustion engine
22: None, # Hybrid Regenerative
23: 'diesel', # Bifuel running diesel
}
# noinspection PyMissingOrEmptyDocstring,PyPep8Naming
[docs] class CMV(co2_utl.Constants):
#: Enable optimization loop?
ENABLE_OPT_LOOP = False
# noinspection PyMissingOrEmptyDocstring,PyPep8Naming
[docs] class default_start_stop_activation_time(co2_utl.Constants):
#: Enable function?
ENABLE = False
#: Start-stop activation time threshold [s].
threshold = 30
# noinspection PyMissingOrEmptyDocstring,PyPep8Naming
[docs] class calculate_last_gear_box_ratio_v1(co2_utl.Constants):
#: Maximum admissible ratio for the last gear [-].
MAX_RATIO = 2
#: Minimum admissible ratio for the last gear [-].
MIN_RATIO = 0.2
#: Delta ratio for calculating the last gear [-].
DELTA_RATIO = 0.1
# noinspection PyMissingOrEmptyDocstring,PyPep8Naming
[docs] class calculate_maximum_velocity(co2_utl.Constants):
#: Maximum admissible velocity for the vehicle maximum velocity [km/h].
MAX_VEL = 1000
#: Minimum admissible velocity for the vehicle maximum velocity [km/h].
MIN_VEL = 1
#: Delta ratio for calculating the vehicle maximum velocity [km/h].
DELTA_VEL = 1
#: Full load curve percentage fro calculating the available power [-].
PREC_FLC = 0.9
# noinspection PyMissingOrEmptyDocstring,PyPep8Naming
[docs] class MGS(co2_utl.Constants):
#: Maximum admissible velocity for the vehicle maximum velocity [km/h].
MAX_VEL = 1000
#: Minimum admissible velocity for the vehicle maximum velocity [km/h].
MIN_VEL = 1
#: Delta ratio for calculating the vehicle maximum velocity [km/h].
DELTA_VEL = 0.1
#: Full load curve percentage fro calculating the available power [-].
PREC_FLC = 0.9
# noinspection PyMissingOrEmptyDocstring,PyPep8Naming
[docs] class calculate_maximum_velocity_v2(co2_utl.Constants):
#: Maximum admissible velocity for the vehicle maximum velocity [km/h].
MAX_VEL = 1000
#: Minimum admissible velocity for the vehicle maximum velocity [km/h].
MIN_VEL = 1
#: Delta ratio for calculating the vehicle maximum velocity [km/h].
DELTA_VEL = 1
# noinspection PyMissingOrEmptyDocstring,PyPep8Naming
[docs] class calculate_first_gear_box_ratio(co2_utl.Constants):
#: Starting slope [-].
STARTING_SLOPE = 0.5 # --> 50%
#: Percentage of maximum engine torque to calculate the first gear [-].
MAX_TORQUE_PERCENTAGE = 0.95
# noinspection PyMissingOrEmptyDocstring,PyPep8Naming
[docs] class design_gear_box_ratios(co2_utl.Constants):
#: Two factor to design the gear box ratios [-].
f_two = [
1, 1.01, 1.02, 1.03, 1.04, 1.05, 1.06, 1.07, 1.08, 1.09, 1.1, 1.11,
1.12, 1.13, 1.14, 1.15, 1.16, 1.17, 1.18, 1.19, 1.2
]
#: Tuning factor to design the gear box ratios [-].
f_tuning = [
1, 1.01, 1.02, 1.03, 1.04, 1.05, 1.06, 1.07, 1.08, 1.09, 1.1
]
# noinspection PyMissingOrEmptyDocstring,PyPep8Naming
class _filter_temperature_samples(co2_utl.Constants):
#: Max abs val of temperature derivatives during the cold phase [°C/s].
max_abs_dt_cold = 0.7
# noinspection PyMissingOrEmptyDocstring,PyPep8Naming
class _rescaling_matrix(co2_utl.Constants):
#: Percentage width top base (i.e., short base) [-].
a = 0.9902
#: Percentage to define the bottom base (i.e., long base) from the phase
#: corner [-].
b = 0.1699
# noinspection PyMissingOrEmptyDocstring,PyPep8Naming
[docs] class calculate_aerodynamic_drag_coefficient_v1(co2_utl.Constants):
#: Aerodynamic drag coefficients function of vehicle body [-].
cw = {
'cabriolet': 0.28, 'sedan': 0.27, 'hatchback': 0.3,
'stationwagon': 0.28, 'suv/crossover': 0.35, 'mpv': 0.3,
'coupé': 0.27, 'bus': 0.35, 'bestelwagen': 0.35,
'pick-up': 0.4 # estimated.
}
# noinspection PyMissingOrEmptyDocstring,PyPep8Naming
[docs] class calculate_aerodynamic_drag_coefficient(co2_utl.Constants):
#: Aerodynamic drag coefficients function of vehicle category [-].
cw = {
'A': 0.34, 'B': 0.31, 'C': 0.29, 'D': 0.30, 'E': 0.30, 'F': 0.28,
'S': 0.29, 'M': 0.32, 'J': 0.38
}
# noinspection PyMissingOrEmptyDocstring,PyPep8Naming
[docs] class calculate_f2(co2_utl.Constants):
#: Deteriorating coefficient of the aerodynamic drag and frontal area
#: due to the roof box [-].
roof_box = 1.2
# noinspection PyMissingOrEmptyDocstring,PyPep8Naming
[docs] class calculate_rolling_resistance_coeff(co2_utl.Constants):
#: Rolling resistance coeff, function of tyre class and category [-].
#: http://eur-lex.europa.eu/legal-content/EN/TXT/PDF/?uri=CELEX:
#: 02009R1222-20120530&from=EN
#: (table A4/1 of eu legislation not world wide)
coeff = {
'C1': {
'A': 5.9, 'B': 7.1, 'C': 8.4, 'D': 9.8, 'E': 9.8, 'F': 11.3,
'G': 12.9
},
'C2': {
'A': 4.9, 'B': 6.1, 'C': 7.4, 'D': 8.6, 'E': 8.6, 'F': 9.9,
'G': 11.2
},
'C3': {
'A': 3.5, 'B': 4.5, 'C': 5.5, 'D': 6.5, 'E': 7.5, 'F': 8.5,
'G': 8.5
}
}
coeff = {k: {i: j / 1000.0 for i, j in v.items()}
for k, v in coeff.items()}
# noinspection PyMissingOrEmptyDocstring,PyPep8Naming
[docs] class calculate_f1(co2_utl.Constants):
#: Linear model coefficients.
qm = 2.7609 / 2, -71.735 / 2
# noinspection PyMissingOrEmptyDocstring,PyPep8Naming
[docs] class calculate_raw_frontal_area_v1(co2_utl.Constants):
#: Frontal area formulas function of vehicle_mass [function].
formulas = dict.fromkeys(
'ABCDEFS', '0.4041 * np.log(vehicle_mass) - 0.338'
)
formulas['J'] = '0.0007 * vehicle_mass + 1.8721'
formulas['M'] = formulas['J']
# noinspection PyMissingOrEmptyDocstring,PyPep8Naming
[docs] class calculate_frontal_area(co2_utl.Constants):
#: Projection factor from the row frontal area (h * w) [-].
projection_factor = 0.84
# noinspection PyMissingOrEmptyDocstring,PyPep8Naming
[docs] class select_prediction_data(co2_utl.Constants):
#: WLTP calibration data used to predict.
prediction_data = [
['calibration', [
'angle_slope', 'alternator_nominal_voltage',
'alternator_efficiency', 'battery_capacity', 'cycle_type',
'cycle_name', 'engine_capacity', 'engine_stroke',
'final_drive_efficiency', 'final_drive_ratios', 'frontal_area',
'final_drive_ratio', 'engine_thermostat_temperature',
'aerodynamic_drag_coefficient', 'fuel_type', 'ignition_type',
'gear_box_type', 'engine_max_power',
'engine_speed_at_max_power', 'rolling_resistance_coeff',
'time_cold_hot_transition', 'engine_idle_fuel_consumption',
'engine_type', 'has_start_stop', 'engine_is_turbo',
'engine_fuel_lower_heating_value', 'f0',
'has_energy_recuperation', 'fuel_carbon_content_percentage',
'f1', 'f2', 'full_load_speeds', 'plateau_acceleration',
'full_load_powers', 'fuel_saving_at_strategy',
'stand_still_torque_ratio', 'lockup_speed_ratio',
'change_gear_window_width', 'alternator_start_window_width',
'stop_velocity', 'min_time_engine_on_after_start',
'min_engine_on_speed', 'max_velocity_full_load_correction',
'is_hybrid', 'tyre_code', 'engine_has_cylinder_deactivation',
'active_cylinder_ratios', 'engine_has_variable_valve_actuation',
'has_torque_converter', 'has_gear_box_thermal_management',
'has_lean_burn', 'ki_multiplicative', 'n_wheel_drive',
'ki_additive', 'has_periodically_regenerating_systems',
'n_dyno_axes', 'has_selective_catalytic_reduction',
'start_stop_activation_time', 'has_exhausted_gas_recirculation',
'engine_n_cylinders', 'initial_drive_battery_state_of_charge',
'motor_p0_speed_ratio', 'motor_p1_speed_ratio',
'motor_p2_speed_ratio', 'motor_p2_planetary_speed_ratio',
'motor_p3_front_speed_ratio', 'motor_p3_rear_speed_ratio',
'motor_p4_front_speed_ratio', 'motor_p4_rear_speed_ratio',
'rcb_correction', 'vehicle_mass', 'speed_distance_correction',
'atct_family_correction_factor', 'is_plugin'
]],
['models', 'all'],
['user', 'all']
]
# noinspection PyMissingOrEmptyDocstring,PyPep8Naming
class _tech_mult_factors(co2_utl.Constants):
#: Multiplication factors of the engine parameters [-].
factors = {
'vva': {
# 0: {},
1: {'a': 0.98, 'l': 0.92},
},
'lb': {
# 0: {},
1: {'a': 1.1, 'b': 0.72, 'c': 0.76, 'a2': 1.25, 'l2': 2.85}
},
'egr': {
# 0: {},
1: {'a': 1.02, 'b': 1.1, 'c': 1.5, 'a2': 1.1}, # positive turbo
2: {'a': 1.02, 'b': 1.1, 'c': 1.5, 'a2': 1.1},
# positive natural aspiration
3: {'b': 1.08, 'c': 1.15, 'a2': 1.1}, # compression
4: {'b': 1.08, 'c': 1.15, 'a2': 1.1} # compression + scr
}
}
# noinspection PyMissingOrEmptyDocstring,PyPep8Naming
[docs] class FMEP_egr(co2_utl.Constants):
#: Exhausted gas recirculation multiplication factors ids [-].
egr_fact_map = {
('positive turbo', False): 1,
('positive natural aspiration', False): 2,
('compression', False): 3,
('compression', True): 4
}
# noinspection PyMissingOrEmptyDocstring,PyPep8Naming
[docs] class identify_co2_emissions(co2_utl.Constants):
#: Number of perturbations to identify the co2_emissions [-].
n_perturbations = 100
#: Enable first step in the co2_params calibration? [-]
enable_first_step = True
#: Enable second step in the co2_params calibration? [-]
enable_second_step = True
#: Enable third step co2_params calibration in perturbation loop? [-]
enable_third_step = False
#: Absolute error in k_refactor between iterations that is acceptable
#: for convergence in perturbation loop [-].
xatol = 1e-4
# noinspection PyMissingOrEmptyDocstring,PyPep8Naming
[docs] class calibrate_co2_params(co2_utl.Constants):
#: Enable first step in the co2_params calibration? [-]
enable_first_step = False
#: Enable second step in the co2_params calibration? [-]
enable_second_step = False
#: Enable third step in the co2_params calibration? [-]
enable_third_step = True
# noinspection PyMissingOrEmptyDocstring,PyPep8Naming
[docs] class BatteryStatusModel(co2_utl.Constants):
#: Minimum delta time to consider valid a charging state to fit charges
#: boundaries [s].
min_delta_time_boundaries = 5
#: Minimum acceptance percentile to fit the bers threshold [%].
min_percentile_bers = 90
#: Minimum delta soc to set the charging boundaries [%].
min_delta_soc = 8
# noinspection PyMissingOrEmptyDocstring,PyPep8Naming
[docs] class default_ki_multiplicative(co2_utl.Constants):
#: Multiplicative correction for vehicles with periodically regenerating
#: systems [-].
ki_multiplicative = {True: 1.05, False: 1.0}
# noinspection PyMissingOrEmptyDocstring,PyPep8Naming
[docs] class define_fmep_model(co2_utl.Constants):
#: Percentage of max full bmep curve used as limit in cylinder
#: deactivation strategy [-].
acr_full_bmep_curve_percentage = 0.72
#: Percentage of max mean piston speeds used as upper limit in cylinder
#: deactivation strategy [-].
acr_max_mean_piston_speeds_percentage = 0.45
#: Percentage of min mean piston speeds used as lower limit in cylinder
#: deactivation strategy [-].
acr_min_mean_piston_speeds_percentage = 1.5
#: Percentage of max full bmep curve used as limit in lean burn
#: strategy [-].
lb_full_bmep_curve_percentage = 0.4
#: Percentage of max mean piston speeds used as limit in lean burn
#: strategy [-].
lb_max_mean_piston_speeds_percentage = 0.6
#: Percentage of max full bmep curve used as limit in exhausted gas
#: recirculation strategy [-].
egr_full_bmep_curve_percentage = 0.5
#: Percentage of max mean piston speeds used as limit in exhausted gas
#: recirculation strategy [-].
egr_max_mean_piston_speeds_percentage = 0.5
# noinspection PyMissingOrEmptyDocstring,PyPep8Naming
[docs] class define_idle_model_detector(co2_utl.Constants):
#: eps parameter of DBSCAN [RPM].
EPS = 100.0
# noinspection PyMissingOrEmptyDocstring,PyPep8Naming
[docs] class identify_idle_engine_speed_std(co2_utl.Constants):
#: Min standard deviation value [RPM].
MIN_STD = 100.0
#: Max standard deviation percentage of median value [-].
MAX_STD_PERC = 0.3
# noinspection PyMissingOrEmptyDocstring,PyPep8Naming
[docs] class StartStopModel(co2_utl.Constants):
#: Maximum allowed velocity to stop the engine [km/h].
stop_velocity = 2.0
#: Minimum acceleration to switch on the engine [m/s2].
plateau_acceleration = 1 / 3.6
# noinspection PyMissingOrEmptyDocstring,PyPep8Naming
[docs] class correct_constant_velocity(co2_utl.Constants):
#: Constant velocities to correct the upper limits for NEDC [km/h].
CON_VEL_UP_SHIFT = (15.0, 32.0, 50.0, 70.0)
#: Window to identify if the shifting matrix has limits close to
# `CON_VEL_UP_SHIFT` [km/h].
VEL_UP_WINDOW = 3.5
#: Delta to add to the limit if this is close to `CON_VEL_UP_SHIFT`
# [km/h].
DV_UP_SHIFT = -0.5
#: Constant velocities to correct the bottom limits for NEDC[km/h].
CON_VEL_DN_SHIFT = (35.0, 50.0)
#: Window to identify if the shifting matrix has limits close to
# `CON_VEL_DN_SHIFT` [km/h].
VEL_DN_WINDOW = 3.5
#: Delta to add to the limit if this is close to `CON_VEL_DN_SHIFT`
# [km/h].
DV_DN_SHIFT = -1
# noinspection PyMissingOrEmptyDocstring,PyPep8Naming
[docs] class define_initial_co2_emission_model_params_guess(co2_utl.Constants):
#: Initial guess CO2 emission model params.
CO2_PARAMS = {
'positive turbo': [
('a', {'value': 0.468678, 'min': 0.0}),
# 'min': 0.398589, 'max': 0.538767},
('b', {'value': 0.011859}), # 'min': 0.006558, 'max': 0.01716},
('c', {'value': -0.00069}),
# 'min': -0.00099, 'max': -0.00038},
('a2', {'value': -0.00266, 'max': 0.0}),
# 'min': -0.00354, 'max': -0.00179},
('b2', {'value': 0, 'min': -1, 'max': 1, 'vary': False}),
('l', {'value': -2.14063, 'max': 0.0}),
# 'min': -3.17876, 'max': -1.1025}),
('l2', {'value': -0.0025, 'max': 0.0}),
# 'min': -0.00796, 'max': 0.0}),b
('t1', {'value': 3.5, 'min': 0.0, 'max': 8.0}),
('dt', {'value': 1.0, 'min': 0.0}),
('t0', {'expr': 't1 + dt', 'value': 4.5, 'max': 8.0}),
],
'positive natural aspiration': [
('a', {'value': 0.4851, 'min': 0.0}),
# 'min': 0.40065, 'max': 0.54315},
('b', {'value': 0.01193}), # 'min': -0.00247, 'max': 0.026333},
('c', {'value': -0.00065}),
# 'min': -0.00138, 'max': 0.0000888},
('a2', {'value': -0.00385, 'max': 0.0}),
# 'min': -0.00663, 'max': -0.00107},
('b2', {'value': 0, 'min': -1, 'max': 1, 'vary': False}),
('l', {'value': -2.39882, 'max': 0.0}),
# 'min': -3.27698, 'max': -1.72066},
('l2', {'value': -0.00286, 'max': 0.0}),
# 'min': -0.00577, 'max': 0.0},
('t1', {'value': 3.5, 'min': 0.0, 'max': 8.0}),
('dt', {'value': 1.0, 'min': 0.0}),
('t0', {'expr': 't1 + dt', 'value': 4.5, 'max': 8.0}),
],
'compression': [
('a', {'value': 0.391197, 'min': 0.0}),
# 'min': 0.346548, 'max': 0.435846},
('b', {'value': 0.028604}),
# 'min': 0.002519, 'max': 0.054688},
('c', {'value': -0.00196}),
# 'min': -0.00386, 'max': -0.000057},
('a2', {'value': -0.0012, 'max': 0.0}),
# 'min': -0.00233, 'max': -0.000064},
('b2', {'value': 0, 'min': -1, 'max': 1, 'vary': False}),
('l', {'value': -1.55291, 'max': 0.0}),
# 'min': -2.2856, 'max': -0.82022},
('l2', {'value': -0.0076, 'max': 0.0}),
# 'min': -0.01852, 'max': 0.0},
('t1', {'value': 3.5, 'min': 0.0, 'max': 8.0}),
('dt', {'value': 1.0, 'min': 0.0}),
('t0', {'expr': 't1 + dt', 'value': 4.5, 'max': 8.0}),
]
}
# noinspection PyMissingOrEmptyDocstring,PyPep8Naming
[docs] class default_specific_gear_shifting(co2_utl.Constants):
#: Specific gear shifting model.
SPECIFIC_GEAR_SHIFTING = 'ALL'
# noinspection PyMissingOrEmptyDocstring,PyPep8Naming
[docs] class default_clutch_k_factor_curve(co2_utl.Constants):
#: Torque ratio when speed ratio==0 for clutch model.
STAND_STILL_TORQUE_RATIO = 1.0
#: Minimum speed ratio where torque ratio==1 for clutch model.
LOCKUP_SPEED_RATIO = 0.0
# noinspection PyMissingOrEmptyDocstring,PyPep8Naming
[docs] class default_tc_k_factor_curve(co2_utl.Constants):
#: Torque ratio when speed ratio==0 for torque converter model.
STAND_STILL_TORQUE_RATIO = 1.9
#: Minimum speed ratio where torque ratio==1 for torque converter model.
LOCKUP_SPEED_RATIO = 0.87
# noinspection PyMissingOrEmptyDocstring,PyPep8Naming
[docs] class select_default_n_dyno_axes(co2_utl.Constants):
#: Number of dyno axes [-].
DYNO_AXES = {
'WLTP': {2: 1, 4: 2},
'NEDC': {2: 1, 4: 1}
}
# noinspection PyMissingOrEmptyDocstring,PyPep8Naming
[docs] class select_phases_integration_times(co2_utl.Constants):
#: Cycle phases integration times [s].
INTEGRATION_TIMES = {
'WLTP': (0.0, 590.0, 1023.0, 1478.0, 1800.0),
'NEDC': (0.0, 780.0, 1180.0)
}
# noinspection PyMissingOrEmptyDocstring,PyPep8Naming
[docs] class get_gear_box_efficiency_constants(co2_utl.Constants):
#: Vehicle gear box efficiency constants (gbp00, gbp10, and gbp01).
PARAMS = {
True: {
'gbp00': {'m': 0.0043233434399999994,
'q': {'hot': 0.29823614099999995,
'cold': 0.695884329}},
'gbp10': {'m': 2.4525999999999996e-06,
'q': {'hot': 0.0001547871,
'cold': 0.0005171569}},
'gbp01': {'q': {'hot': 0.9793688500000001,
'cold': 0.96921995}},
},
False: {
'gbp00': {'m': 0.0043233434399999994,
'q': {'hot': 0.29823614099999995,
'cold': 0.695884329}},
'gbp10': {'m': 2.4525999999999996e-06,
'q': {'hot': 5.15957e-05,
'cold': 0.00012958919999999998}},
'gbp01': {'q': {'hot': 0.9895177500000001,
'cold': 0.9793688500000001}},
}
}
# noinspection PyMissingOrEmptyDocstring,PyPep8Naming
[docs] class calculate_engine_mass(co2_utl.Constants):
#: Equivalent gear box heat capacity parameters.
PARAMS = {
'mass_coeff': {
'compression': 1.1,
'positive': 1.0
},
'mass_reg_coeff': (0.4208, 60)
}
# noinspection PyMissingOrEmptyDocstring,PyPep8Naming
[docs] class calculate_engine_heat_capacity(co2_utl.Constants):
#: Equivalent gear box heat capacity parameters.
PARAMS = {
'heated_mass_percentage': {
'coolant': 0.04, # coolant: 50%/50% (0.85*4.186)
'oil': 0.055, # oil: lubricant
'crankcase': 0.18, # crankcase: cast iron
'cyl_head': 0.09, # cyl_head: aluminium
'pistons': 0.025, # pistons: aluminium
'crankshaft': 0.08, # crankshaft: steel
'body': 0.1 # body: cast iron
},
# Cp in (J/kgK)
'heat_capacity': {
'coolant': 0.85 * 4186.0,
'oil': 2090.0,
'crankcase': 460.0,
'cyl_head': 910.0,
'pistons': 910.0,
'crankshaft': 490.0,
'body': 460.0
}
}
# noinspection PyMissingOrEmptyDocstring,PyPep8Naming
[docs] class calculate_equivalent_gear_box_heat_capacity(co2_utl.Constants):
#: Equivalent gear box heat capacity parameters.
PARAMS = {
'gear_box_mass_engine_ratio': 0.25,
# Cp in (J/kgK)
'heat_capacity': {
'oil': 2090.0,
},
'thermal_management_factor': 0.5
}
# noinspection PyMissingOrEmptyDocstring,PyPep8Naming
[docs] class default_full_load_speeds_and_powers(co2_utl.Constants):
#: Vehicle normalized full load curve.
FULL_LOAD = {
'positive': (
[0., 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1., 1.1, 1.2],
[0.1, 0.198238659, 0.30313392, 0.410104642, 0.516920841,
0.621300767, 0.723313491, 0.820780368, 0.901750158,
0.962968496, 0.995867804, 0.953356174, 0.85]
),
'compression': (
[0., 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1., 1.1, 1.2],
[0.1, 0.278071182, 0.427366185, 0.572340499, 0.683251935,
0.772776746, 0.846217049, 0.906754984, 0.94977083, 0.981937981,
1, 0.937598144, 0.85]
)
}
# noinspection PyMissingOrEmptyDocstring,PyPep8Naming
[docs] class calculate_engine_max_torque(co2_utl.Constants):
#: Engine nominal torque params.
PARAMS = {
'positive': 1.25,
'compression': 1.1
}
# noinspection PyMissingOrEmptyDocstring,PyPep8Naming
[docs] class calculate_engine_moment_inertia(co2_utl.Constants):
#: Engine moment of inertia params.
PARAMS = {
'positive': 1,
'compression': 2
}
# noinspection PyMissingOrEmptyDocstring,PyPep8Naming
[docs] class calculate_co2_emissions(co2_utl.Constants):
# idle ratio to define the fuel cutoff [-].
cutoff_idle_ratio = 1.1
# noinspection PyMissingOrEmptyDocstring,PyPep8Naming
[docs] class default_initial_drive_battery_state_of_charge(co2_utl.Constants):
# default initial state of charge of the drive battery [%].
initial_state_of_charge = {
'none': sh.NONE,
'mild': 50,
'full': 60,
'plugin': 70,
'electric': 80
}
# noinspection PyMissingOrEmptyDocstring,PyPep8Naming
[docs] class default_initial_service_battery_state_of_charge(co2_utl.Constants):
# default initial state of charge of the service battery [%].
initial_state_of_charge = {
'WLTP': 90,
'NEDC': 99
}
# noinspection PyMissingOrEmptyDocstring,PyPep8Naming
[docs] class default_fuel_density(co2_utl.Constants):
#: Fuel density [g/l].
FUEL_DENSITY = {
'gasoline': 745.0,
'diesel': 832.0,
'LPG': 43200.0 / 46000.0 * 745.0, # Gasoline equivalent.
'NG': 43200.0 / 45100.0 * 745.0, # Gasoline equivalent.
'ethanol': 794.0,
'methanol': 791.0,
'propane': 510.0,
'biodiesel': 890.0,
}
# noinspection PyMissingOrEmptyDocstring,PyPep8Naming
[docs] class default_fuel_lower_heating_value(co2_utl.Constants):
#: Fuel lower heating value [kJ/kg].
LHV = {
'gasoline': 43200.0,
'diesel': 43100.0,
'LPG': 46000.0,
'NG': 45100.0,
'ethanol': 26800.0,
'methanol': 19800.0,
'propane': 49680.0,
'biodiesel': 37900.0,
}
# noinspection PyMissingOrEmptyDocstring,PyPep8Naming
[docs] class default_fuel_carbon_content(co2_utl.Constants):
#: Fuel carbon content [CO2g/g].
CARBON_CONTENT = {
'gasoline': 3.17,
'diesel': 3.16,
'LPG': 1.35,
'NG': 3.21,
'ethanol': 1.91,
'methanol': 1.37,
'propane': 2.99,
'biodiesel': 2.81,
}
# noinspection PyMissingOrEmptyDocstring,PyPep8Naming
class _rcb_correction(co2_utl.Constants):
#: Willans factors [gCO2/MJ].
WILLANS = {
'gasoline': {
'positive turbo': 184, 'positive natural aspiration': 174
},
'diesel': {'compression': 161},
'LPG': {'positive turbo': 164, 'positive natural aspiration': 155},
'NG': {'positive turbo': 137, 'positive natural aspiration': 129},
'ethanol': {
'positive turbo': 179, 'positive natural aspiration': 169
},
'methanol': {
'positive turbo': 184 * 1.37 / 3.17 * 432.0 / 198.0,
'positive natural aspiration': 174 * 1.37 / 3.17 * 432.0 / 198.0
},
'propane': {
'positive turbo': 184 * 2.99 / 3.17 * 432.0 / 496.8,
'positive natural aspiration': 174 * 2.99 / 3.17 * 432.0 / 496.8
},
'biodiesel': {'compression': 161 * 2.81 / 3.16 * 431.0 / 379.0},
}
# noinspection PyPep8Naming,PyMissingOrEmptyDocstring
[docs]class Defaults(co2_utl.Constants):
values = Values()
functions = Functions()
#: Machine error.
EPS = 1.1920929e-07
#: Infinite value.
INF = 10000.0
dfl = Defaults()
