| public GetDecelCoast ( double speed, double acc, double gradient ) : double | ||
| speed | double | |
| acc | double | |
| gradient | double | |
| return | double | |
public double GetDecelCoast(double speed, double acc, double gradient)
{
//Declaration
int upperIndex;
int lowerIndex;
if (speed < Constants.SPEED_DCEL_MIN)
{
return speed / Constants.SPEED_DCEL_MIN * GetDecelCoast(Constants.SPEED_DCEL_MIN, acc, gradient);
}
double rotCoeff = GetRotationalCoeffecient(speed);
FindLowerUpperInPattern(out lowerIndex, out upperIndex, _speedPatternRotational, speed);
double iGear = Interpolate(speed,
_speedPatternRotational[lowerIndex],
_speedPatternRotational[upperIndex],
_gearTransmissionCurve[lowerIndex],
_gearTransmissionCurve[upperIndex]);
double iTot = iGear * _axleRatio;
double n = (30 * speed * iTot) / ((_effectiveWheelDiameter / 2) * Math.PI);
double nNorm = (n - _engineIdlingSpeed) / (_engineRatedSpeed - _engineIdlingSpeed);
FindLowerUpperInPattern(out lowerIndex, out upperIndex, _nNormTable, nNorm);
double fMot = 0;
if (speed >= 10e-2)
{
fMot = (-Interpolate(nNorm,
_nNormTable[lowerIndex],
_nNormTable[upperIndex],
_dragNormTable[lowerIndex],
_dragNormTable[upperIndex]) * _ratedPower * 1000 / speed) / 0.9;
}
double fRoll = (_resistanceF0
+ _resistanceF1 * speed
+ Math.Pow(_resistanceF2 * speed, 2)
+ Math.Pow(_resistanceF3 * speed, 3)
+ Math.Pow(_resistanceF4 * speed, 4)) * (_massVehicle + _vehicleLoading) * Constants.GRAVITY_CONST;
double fAir = _cWValue * _crossSectionalArea * 1.2 * 0.5 * Math.Pow(speed, 2);
double fGrad = (_massVehicle + _vehicleLoading) * Constants.GRAVITY_CONST * gradient / 100;
return -(fMot + fRoll + fAir + fGrad) / ((_massVehicle + _vehicleLoading) * rotCoeff);
}
#endregion static public VehicleResult CreateVehicleStateData(Helpers Helper,
CEP currCep,
double time,
double inputSpeed,
double inputAcc,
double Gradient = 0,
Correction DataCor = null)
{
//Declaration
double speed = Math.Max(inputSpeed, 0);
double acc;
double P_pos;
//Speed/Acceleration limitation
if (speed == 0)
{
acc = 0;
}
else
{
acc = Math.Min(inputAcc, currCep.GetMaxAccel(speed, Gradient, (Helper.pClass == Constants.strBEV | Helper.uClass == Constants.strHybrid)));
}
//Calculate the power
double power = currCep.CalcPower(speed, acc, Gradient, (Helper.pClass == Constants.strBEV | Helper.uClass == Constants.strHybrid));
double P_eng = currCep.CalcEngPower(power);
double Pwheel = 0;
if (Helper.uClass == Constants.strHybrid)
{
Pwheel = currCep.CalcWheelPower(speed, acc, Gradient);
}
//Power limitation
if (P_eng >= 0)
{
P_pos = power;
}
else
{
P_pos = 0;
}
//Calculate the result values (BEV)
if (Helper.pClass == Constants.strBEV)
{
return(new VehicleResult(Helper.gClass,
"",
time,
speed,
Gradient,
P_eng,
P_pos,
P_eng / currCep.RatedPower,
P_eng / currCep.DrivingPower,
acc,
currCep.GetAllEmission(P_eng, speed, Helper)));
}
//Calculate the decel costing
double decelCoast = currCep.GetDecelCoast(speed, acc, Gradient);
//Calculate the result values (Zero emissions by costing, Idling emissions by v <= 0.5m / s²)
if (acc >= decelCoast || speed <= Constants.ZERO_SPEED_ACCURACY)
{
if (Helper.uClass == Constants.strHybrid)
{
return(new VehicleResult(Helper.gClass,
"",
time,
speed,
Gradient,
P_eng,
P_pos,
P_eng / currCep.RatedPower,
P_eng / currCep.DrivingPower,
acc,
currCep.GetAllEmission(Pwheel, speed, Helper)));
}
else
{
return(new VehicleResult(Helper.gClass,
"",
time,
speed,
Gradient,
P_eng,
P_pos,
P_eng / currCep.RatedPower,
P_eng / currCep.DrivingPower,
acc,
currCep.GetAllEmission(P_eng, speed, Helper)));
}
}
else
{
if (Helper.uClass == Constants.strHybrid)
{
return(new VehicleResult(Helper.gClass,
"",
time,
speed,
Gradient,
P_eng,
P_pos,
P_eng / currCep.RatedPower,
P_eng / currCep.DrivingPower,
acc,
currCep.GetAllEmission(Pwheel, speed, Helper, true)));
}
else
{
return(new VehicleResult(Helper.gClass,
"",
time,
speed,
Gradient,
P_eng,
P_pos,
P_eng / currCep.RatedPower,
P_eng / currCep.DrivingPower,
acc,
currCep.GetAllEmission(P_eng, speed, Helper, true)));
}
}
}
