Eclipse SUMO - Simulation of Urban MObility
Loading...
Searching...
No Matches
HelpersHBEFA.h
Go to the documentation of this file.
1/****************************************************************************/
2// Eclipse SUMO, Simulation of Urban MObility; see https://eclipse.dev/sumo
3// Copyright (C) 2001-2024 German Aerospace Center (DLR) and others.
4// This program and the accompanying materials are made available under the
5// terms of the Eclipse Public License 2.0 which is available at
6// https://www.eclipse.org/legal/epl-2.0/
7// This Source Code may also be made available under the following Secondary
8// Licenses when the conditions for such availability set forth in the Eclipse
9// Public License 2.0 are satisfied: GNU General Public License, version 2
10// or later which is available at
11// https://www.gnu.org/licenses/old-licenses/gpl-2.0-standalone.html
12// SPDX-License-Identifier: EPL-2.0 OR GPL-2.0-or-later
13/****************************************************************************/
19// Helper methods for HBEFA-based emission computation
20/****************************************************************************/
21#pragma once
22#include <config.h>
23
24#include <vector>
25#include <limits>
26#include <cmath>
30#include "EnergyParams.h"
31#include "PollutantsInterface.h"
33
34// ===========================================================================
35// class definitions
36// ===========================================================================
46private:
47 static const int HBEFA_BASE = 1 << 16;
48
49public:
53
54
68 inline double compute(const SUMOEmissionClass c, const PollutantsInterface::EmissionType e, const double v, const double a, const double slope, const EnergyParams* param) const {
69 if (e == PollutantsInterface::ELEC || (param != nullptr && param->isEngineOff())) {
70 return 0.;
71 }
72 if (v > ZERO_SPEED_ACCURACY && a < getCoastingDecel(c, v, a, slope, param)) {
73 return 0.;
74 }
75 const int index = (c & ~PollutantsInterface::HEAVY_BIT) - HBEFA_BASE;
76 const double kmh = v * 3.6;
77 const double scale = (e == PollutantsInterface::FUEL && myVolumetricFuel) ? 3.6 * 790. : 3.6;
78 if (index >= 42) {
79 const double* f = myFunctionParameter[index - 42] + 6 * e;
80 return MAX2((f[0] + f[3] * kmh + f[4] * kmh * kmh + f[5] * kmh * kmh * kmh) / scale, 0.);
81 }
82 const double* f = myFunctionParameter[index] + 6 * e;
83 const double alpha = RAD2DEG(asin(a / GRAVITY));
84 return MAX2((f[0] + f[1] * alpha * kmh + f[2] * alpha * alpha * kmh + f[3] * kmh + f[4] * kmh * kmh + f[5] * kmh * kmh * kmh) / scale, 0.);
85 }
86
87
88private:
90 static double myFunctionParameter[42][36];
91
92};
#define GRAVITY
Definition GeomHelper.h:37
#define RAD2DEG(x)
Definition GeomHelper.h:36
int SUMOEmissionClass
T MAX2(T a, T b)
Definition StdDefs.h:82
An upper class for objects with additional parameters.
bool isEngineOff() const
Returns the state of the engine when the vehicle is not moving.
Helper methods for HBEFA-based emission computation.
static const int HBEFA_BASE
double compute(const SUMOEmissionClass c, const PollutantsInterface::EmissionType e, const double v, const double a, const double slope, const EnergyParams *param) const
Computes the emitted pollutant amount using the given speed and acceleration.
static double myFunctionParameter[42][36]
The function parameter.
HelpersHBEFA()
Constructor (initializes myEmissionClassStrings)
zero emission model, used as superclass for the other model helpers
virtual double getCoastingDecel(const SUMOEmissionClass c, const double v, const double a, const double slope, const EnergyParams *param) const
Returns the maximum deceleration value (as a negative number), which can still be considered as non-b...
static const double ZERO_SPEED_ACCURACY
the lowest speed which allows reliable coasting calculations
bool myVolumetricFuel
return fuel consumption in l instead of mg
Helper methods for PHEMlight-based emission computation.
EmissionType
Enumerating all emission types, including fuel.