Code development platform for open source projects from the European Union institutions :large_blue_circle: EU Login authentication by SMS will be completely phased out by mid-2025. To see alternatives please check here

Skip to content
Snippets Groups Projects
Commit cc50fca8 authored by Markus Quaritsch's avatar Markus Quaritsch Committed by Raphael LUZ
Browse files

analytic ratedPower implementation

parent 7a542ac4
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
VectoCore/Models/SimulationComponent/Data/Engine/FullLoadCurve.cs
+ 304
− 264
View file @ cc50fca8
using System.Collections.Generic;
using System;
using System.Collections.Generic;
using System.Data;
using System.Diagnostics.Contracts;
using System.Linq;
using Common.Logging;
using Newtonsoft.Json;
using TUGraz.VectoCore.Exceptions;
using TUGraz.VectoCore.Models.SimulationComponent.Data;
using TUGraz.VectoCore.Utils;
namespace TUGraz.VectoCore.Models.SimulationComponent.Data.Engine
{
/// <summary>
/// Represents the Full load curve.
/// </summary>
public class FullLoadCurve : SimulationComponentData
{
[JsonProperty] private List<FullLoadCurveEntry> _entries;
public static FullLoadCurve ReadFromFile(string fileName)
{
var data = VectoCSVFile.Read(fileName);
//todo Contract.Requires<VectoException>(data.Columns.Count != 4, "FullLoadCurve Data File must consist of 4 columns.");
if (data.Columns.Count != 4) {
throw new VectoException("FullLoadCurve Data File must consist of 4 columns.");
}
//todo Contract.Requires<VectoException>(data.Rows.Count < 2, "FullLoadCurve must consist of at least two lines with numeric values (below file header)");
if (data.Rows.Count < 2) {
throw new VectoException(
"FullLoadCurve must consist of at least two lines with numeric values (below file header)");
}
List<FullLoadCurveEntry> entries;
if (HeaderIsValid(data.Columns)) {
entries = CreateFromColumnNames(data);
} else {
var log = LogManager.GetLogger<FullLoadCurve>();
log.WarnFormat(
"FullLoadCurve: Header Line is not valid. Expected: '{0}, {1}, {2}, {3}', Got: '{4}'. Falling back to column index.",
Fields.EngineSpeed, Fields.TorqueFullLoad, Fields.TorqueDrag, Fields.PT1,
string.Join(", ", data.Columns.Cast<DataColumn>().Select(c => c.ColumnName).Reverse()));
entries = CreateFromColumnIndizes(data);
}
return new FullLoadCurve { _entries = entries };
}
private static bool HeaderIsValid(DataColumnCollection columns)
{
Contract.Requires(columns != null);
return columns.Contains(Fields.EngineSpeed)
&& columns.Contains(Fields.TorqueDrag)
&& columns.Contains(Fields.TorqueFullLoad)
&& columns.Contains(Fields.PT1);
}
private static List<FullLoadCurveEntry> CreateFromColumnNames(DataTable data)
{
Contract.Requires(data != null);
return (from DataRow row in data.Rows
select new FullLoadCurveEntry {
EngineSpeed = row.ParseDouble(Fields.EngineSpeed).SI().Rounds.Per.Minute.To<PerSecond>(),
TorqueFullLoad = row.ParseDouble(Fields.TorqueFullLoad).SI<NewtonMeter>(),
TorqueDrag = row.ParseDouble(Fields.TorqueDrag).SI<NewtonMeter>(),
PT1 = row.ParseDouble(Fields.PT1).SI<Second>()
}).ToList();
}
private static List<FullLoadCurveEntry> CreateFromColumnIndizes(DataTable data)
{
Contract.Requires(data != null);
return (from DataRow row in data.Rows
select new FullLoadCurveEntry {
EngineSpeed = row.ParseDouble(0).SI().Rounds.Per.Minute.To<PerSecond>(),
TorqueFullLoad = row.ParseDouble(1).SI<NewtonMeter>(),
TorqueDrag = row.ParseDouble(2).SI<NewtonMeter>(),
PT1 = row.ParseDouble(3).SI<Second>()
}).ToList();
}
/// <summary>
/// [rad/s] => [Nm]
/// </summary>
/// <param name="angularFrequency">[rad/s]</param>
/// <returns>[Nm]</returns>
public NewtonMeter FullLoadStationaryTorque(PerSecond angularFrequency)
{
var idx = FindIndex(angularFrequency);
return VectoMath.Interpolate((double) _entries[idx - 1].EngineSpeed, (double) _entries[idx].EngineSpeed,
(double) _entries[idx - 1].TorqueFullLoad, (double) _entries[idx].TorqueFullLoad,
(double) angularFrequency).SI<NewtonMeter>();
}
/// <summary>
/// [rad/s] => [W]
/// </summary>
/// <param name="angularFrequency">[rad/s]</param>
/// <returns>[W]</returns>
public Watt FullLoadStationaryPower(PerSecond angularFrequency)
{
return Formulas.TorqueToPower(FullLoadStationaryTorque(angularFrequency), angularFrequency);
}
/// <summary>
/// [rad/s] => [Nm]
/// </summary>
/// <param name="angularFrequency">[rad/s]</param>
/// <returns>[Nm]</returns>
public NewtonMeter DragLoadStationaryTorque(PerSecond angularFrequency)
{
var idx = FindIndex(angularFrequency);
return VectoMath.Interpolate((double) _entries[idx - 1].EngineSpeed, (double) _entries[idx].EngineSpeed,
(double) _entries[idx - 1].TorqueDrag, (double) _entries[idx].TorqueDrag,
(double) angularFrequency).SI<NewtonMeter>();
}
/// <summary>
/// [rad/s] => [W].
/// </summary>
/// <param name="angularFrequency">[rad/s]</param>
/// <returns>[W]</returns>
public Watt DragLoadStationaryPower(PerSecond angularFrequency)
{
Contract.Requires(angularFrequency.HasEqualUnit(new SI().Radian.Per.Second));
Contract.Ensures(Contract.Result<SI>().HasEqualUnit(new SI().Watt));
return Formulas.TorqueToPower(DragLoadStationaryTorque(angularFrequency), angularFrequency);
}
/// <summary>
/// [rad/s] => [-]
/// </summary>
/// <param name="angularFrequency">[rad/s]</param>
/// <returns>[-]</returns>
public SI PT1(SI angularFrequency)
{
Contract.Requires(angularFrequency.HasEqualUnit(new SI().Radian.Per.Second));
Contract.Ensures(Contract.Result<SI>().HasEqualUnit(new SI()));
var idx = FindIndex(angularFrequency);
return VectoMath.Interpolate((double) _entries[idx - 1].EngineSpeed, (double) _entries[idx].EngineSpeed,
(double) _entries[idx - 1].PT1, (double) _entries[idx].PT1,
(double) angularFrequency).SI();
}
/// <summary>
/// [rad/s] => index. Get item index for engineSpeed.
/// </summary>
/// <param name="engineSpeed">[rad/s]</param>
/// <returns>index</returns>
protected int FindIndex(SI engineSpeed)
{
Contract.Requires(engineSpeed.HasEqualUnit(new SI().Radian.Per.Second));
int idx;
if (engineSpeed < _entries[0].EngineSpeed) {
Log.ErrorFormat("requested rpm below minimum rpm in FLD curve - extrapolating. n: {0}, rpm_min: {1}",
engineSpeed.To().Rounds.Per.Minute, _entries[0].EngineSpeed.To().Rounds.Per.Minute);
idx = 1;
} else {
idx = _entries.FindIndex(x => x.EngineSpeed > engineSpeed);
}
if (idx <= 0) {
idx = engineSpeed > _entries[0].EngineSpeed ? _entries.Count - 1 : 1;
}
return idx;
}
private static class Fields
{
/// <summary>
/// [rpm] engine speed
/// </summary>
public const string EngineSpeed = "engine speed";
/// <summary>
/// [Nm] full load torque
/// </summary>
public const string TorqueFullLoad = "full load torque";
/// <summary>
/// [Nm] motoring torque
/// </summary>
public const string TorqueDrag = "motoring torque";
/// <summary>
/// [s] time constant
/// </summary>
public const string PT1 = "PT1";
}
private class FullLoadCurveEntry
{
/// <summary>
/// [rad/s] engine speed
/// </summary>
public PerSecond EngineSpeed { get; set; }
/// <summary>
/// [Nm] full load torque
/// </summary>
public NewtonMeter TorqueFullLoad { get; set; }
/// <summary>
/// [Nm] motoring torque
/// </summary>
public NewtonMeter TorqueDrag { get; set; }
/// <summary>
/// [s] PT1 time constant
/// </summary>
public Second PT1 { get; set; }
#region Equality members
protected bool Equals(FullLoadCurveEntry other)
{
Contract.Requires(other != null);
return EngineSpeed.Equals(other.EngineSpeed)
&& TorqueFullLoad.Equals(other.TorqueFullLoad)
&& TorqueDrag.Equals(other.TorqueDrag)
&& PT1.Equals(other.PT1);
}
public override bool Equals(object obj)
{
if (ReferenceEquals(null, obj)) {
return false;
}
if (ReferenceEquals(this, obj)) {
return true;
}
return obj.GetType() == GetType() && Equals((FullLoadCurveEntry) obj);
}
public override int GetHashCode()
{
var hashCode = EngineSpeed.GetHashCode();
hashCode = (hashCode * 397) ^ TorqueFullLoad.GetHashCode();
hashCode = (hashCode * 397) ^ TorqueDrag.GetHashCode();
hashCode = (hashCode * 397) ^ PT1.GetHashCode();
return hashCode;
}
#endregion
}
#region Equality members
protected bool Equals(FullLoadCurve other)
{
return _entries.SequenceEqual(other._entries);
}
public override bool Equals(object obj)
{
if (ReferenceEquals(null, obj)) {
return false;
}
if (ReferenceEquals(this, obj)) {
return true;
}
return obj.GetType() == GetType() && Equals((FullLoadCurve) obj);
}
public override int GetHashCode()
{
return (_entries != null ? _entries.GetHashCode() : 0);
}
#endregion
}
public class FullLoadCurve : SimulationComponentData
{
[JsonProperty] private List<FullLoadCurveEntry> _entries;
public static FullLoadCurve ReadFromFile(string fileName)
{
var data = VectoCSVFile.Read(fileName);
//todo Contract.Requires<VectoException>(data.Columns.Count != 4, "FullLoadCurve Data File must consist of 4 columns.");
if (data.Columns.Count != 4) {
throw new VectoException("FullLoadCurve Data File must consist of 4 columns.");
}
//todo Contract.Requires<VectoException>(data.Rows.Count < 2, "FullLoadCurve must consist of at least two lines with numeric values (below file header)");
if (data.Rows.Count < 2) {
throw new VectoException("FullLoadCurve must consist of at least two lines with numeric values (below file header)");
}
List<FullLoadCurveEntry> entries;
if (HeaderIsValid(data.Columns)) {
entries = CreateFromColumnNames(data);
} else {
var log = LogManager.GetLogger<FullLoadCurve>();
log.WarnFormat(
"FullLoadCurve: Header Line is not valid. Expected: '{0}, {1}, {2}, {3}', Got: '{4}'. Falling back to column index.",
Fields.EngineSpeed, Fields.TorqueFullLoad, Fields.TorqueDrag, Fields.PT1,
string.Join(", ", data.Columns.Cast<DataColumn>().Select(c => c.ColumnName).Reverse()));
entries = CreateFromColumnIndizes(data);
}
return new FullLoadCurve { _entries = entries };
}
private static bool HeaderIsValid(DataColumnCollection columns)
{
Contract.Requires(columns != null);
return columns.Contains(Fields.EngineSpeed)
&& columns.Contains(Fields.TorqueDrag)
&& columns.Contains(Fields.TorqueFullLoad)
&& columns.Contains(Fields.PT1);
}
private static List<FullLoadCurveEntry> CreateFromColumnNames(DataTable data)
{
Contract.Requires(data != null);
return (from DataRow row in data.Rows
select new FullLoadCurveEntry {
EngineSpeed = row.ParseDouble(Fields.EngineSpeed).SI().Rounds.Per.Minute.To<PerSecond>(),
TorqueFullLoad = row.ParseDouble(Fields.TorqueFullLoad).SI<NewtonMeter>(),
TorqueDrag = row.ParseDouble(Fields.TorqueDrag).SI<NewtonMeter>(),
PT1 = row.ParseDouble(Fields.PT1).SI<Second>()
}).ToList();
}
private static List<FullLoadCurveEntry> CreateFromColumnIndizes(DataTable data)
{
Contract.Requires(data != null);
return (from DataRow row in data.Rows
select new FullLoadCurveEntry {
EngineSpeed = row.ParseDouble(0).SI().Rounds.Per.Minute.To<PerSecond>(),
TorqueFullLoad = row.ParseDouble(1).SI<NewtonMeter>(),
TorqueDrag = row.ParseDouble(2).SI<NewtonMeter>(),
PT1 = row.ParseDouble(3).SI<Second>()
}).ToList();
}
/// <summary>
/// [rad/s] => [Nm]
/// </summary>
/// <param name="angularFrequency">[rad/s]</param>
/// <returns>[Nm]</returns>
public NewtonMeter FullLoadStationaryTorque(PerSecond angularFrequency)
{
var idx = FindIndex(angularFrequency);
return VectoMath.Interpolate((double) _entries[idx - 1].EngineSpeed, (double) _entries[idx].EngineSpeed,
(double) _entries[idx - 1].TorqueFullLoad, (double) _entries[idx].TorqueFullLoad,
(double) angularFrequency).SI<NewtonMeter>();
}
/// <summary>
/// [rad/s] => [W]
/// </summary>
/// <param name="angularFrequency">[rad/s]</param>
/// <returns>[W]</returns>
public Watt FullLoadStationaryPower(PerSecond angularFrequency)
{
return Formulas.TorqueToPower(FullLoadStationaryTorque(angularFrequency), angularFrequency);
}
/// <summary>
/// [rad/s] => [Nm]
/// </summary>
/// <param name="angularFrequency">[rad/s]</param>
/// <returns>[Nm]</returns>
public NewtonMeter DragLoadStationaryTorque(PerSecond angularFrequency)
{
var idx = FindIndex(angularFrequency);
return VectoMath.Interpolate((double) _entries[idx - 1].EngineSpeed, (double) _entries[idx].EngineSpeed,
(double) _entries[idx - 1].TorqueDrag, (double) _entries[idx].TorqueDrag,
(double) angularFrequency).SI<NewtonMeter>();
}
/// <summary>
/// [rad/s] => [W].
/// </summary>
/// <param name="angularFrequency">[rad/s]</param>
/// <returns>[W]</returns>
public Watt DragLoadStationaryPower(PerSecond angularFrequency)
{
Contract.Requires(angularFrequency.HasEqualUnit(new SI().Radian.Per.Second));
Contract.Ensures(Contract.Result<SI>().HasEqualUnit(new SI().Watt));
return Formulas.TorqueToPower(DragLoadStationaryTorque(angularFrequency), angularFrequency);
}
/// <summary>
/// [rad/s] => [-]
/// </summary>
/// <param name="angularFrequency">[rad/s]</param>
/// <returns>[-]</returns>
public SI PT1(SI angularFrequency)
{
Contract.Requires(angularFrequency.HasEqualUnit(new SI().Radian.Per.Second));
Contract.Ensures(Contract.Result<SI>().HasEqualUnit(new SI()));
var idx = FindIndex(angularFrequency);
return VectoMath.Interpolate((double) _entries[idx - 1].EngineSpeed, (double) _entries[idx].EngineSpeed,
(double) _entries[idx - 1].PT1, (double) _entries[idx].PT1,
(double) angularFrequency).SI();
}
/// <summary>
/// [rad/s] => index. Get item index for engineSpeed.
/// </summary>
/// <param name="engineSpeed">[rad/s]</param>
/// <returns>index</returns>
protected int FindIndex(SI engineSpeed)
{
Contract.Requires(engineSpeed.HasEqualUnit(new SI().Radian.Per.Second));
int idx;
if (engineSpeed < _entries[0].EngineSpeed) {
Log.ErrorFormat("requested rpm below minimum rpm in FLD curve - extrapolating. n: {0}, rpm_min: {1}",
engineSpeed.To().Rounds.Per.Minute, _entries[0].EngineSpeed.To().Rounds.Per.Minute);
idx = 1;
} else {
idx = _entries.FindIndex(x => x.EngineSpeed > engineSpeed);
}
if (idx <= 0) {
idx = engineSpeed > _entries[0].EngineSpeed ? _entries.Count - 1 : 1;
}
return idx;
}
private static class Fields
{
/// <summary>
/// [rpm] engine speed
/// </summary>
public const string EngineSpeed = "engine speed";
/// <summary>
/// [Nm] full load torque
/// </summary>
public const string TorqueFullLoad = "full load torque";
/// <summary>
/// [Nm] motoring torque
/// </summary>
public const string TorqueDrag = "motoring torque";
/// <summary>
/// [s] time constant
/// </summary>
public const string PT1 = "PT1";
}
private class FullLoadCurveEntry
{
/// <summary>
/// [rad/s] engine speed
/// </summary>
public PerSecond EngineSpeed { get; set; }
/// <summary>
/// [Nm] full load torque
/// </summary>
public NewtonMeter TorqueFullLoad { get; set; }
/// <summary>
/// [Nm] motoring torque
/// </summary>
public NewtonMeter TorqueDrag { get; set; }
/// <summary>
/// [s] PT1 time constant
/// </summary>
public Second PT1 { get; set; }
#region Equality members
protected bool Equals(FullLoadCurveEntry other)
{
Contract.Requires(other != null);
return EngineSpeed.Equals(other.EngineSpeed)
&& TorqueFullLoad.Equals(other.TorqueFullLoad)
&& TorqueDrag.Equals(other.TorqueDrag)
&& PT1.Equals(other.PT1);
}
public override bool Equals(object obj)
{
if (ReferenceEquals(null, obj)) {
return false;
}
if (ReferenceEquals(this, obj)) {
return true;
}
return obj.GetType() == GetType() && Equals((FullLoadCurveEntry) obj);
}
public override int GetHashCode()
{
var hashCode = EngineSpeed.GetHashCode();
hashCode = (hashCode * 397) ^ TorqueFullLoad.GetHashCode();
hashCode = (hashCode * 397) ^ TorqueDrag.GetHashCode();
hashCode = (hashCode * 397) ^ PT1.GetHashCode();
return hashCode;
}
#endregion
}
#region Equality members
protected bool Equals(FullLoadCurve other)
{
return _entries.SequenceEqual(other._entries);
}
public override bool Equals(object obj)
{
if (ReferenceEquals(null, obj)) {
return false;
}
if (ReferenceEquals(this, obj)) {
return true;
}
return obj.GetType() == GetType() && Equals((FullLoadCurve) obj);
}
public override int GetHashCode()
{
return (_entries != null ? _entries.GetHashCode() : 0);
}
#endregion
private Tuple<PerSecond, Watt> FindMaxPower(FullLoadCurveEntry p1, FullLoadCurveEntry p2)
{
if (p1.EngineSpeed == p2.EngineSpeed) {
return new Tuple<PerSecond, Watt>(p1.EngineSpeed, Formulas.TorqueToPower(p1.TorqueFullLoad, p1.EngineSpeed));
}
if (p2.EngineSpeed < p1.EngineSpeed) {
var tmp = p1;
p1 = p2;
p2 = tmp;
}
// y = kx + d
var k = (p2.TorqueFullLoad - p1.TorqueFullLoad) / (p2.EngineSpeed - p1.EngineSpeed);
var d = p2.TorqueFullLoad - k * p2.EngineSpeed;
if (k == 0.0.SI()) {
return new Tuple<PerSecond, Watt>(p2.EngineSpeed, Formulas.TorqueToPower(p2.TorqueFullLoad, p2.EngineSpeed));
}
var engineSpeedMaxPower = (-1 * d / (2 * k)).To<PerSecond>();
if (engineSpeedMaxPower < p1.EngineSpeed || engineSpeedMaxPower > p2.EngineSpeed) {
if (k > 0) {
return new Tuple<PerSecond, Watt>(p2.EngineSpeed, Formulas.TorqueToPower(p2.TorqueFullLoad, p2.EngineSpeed));
}
return new Tuple<PerSecond, Watt>(p1.EngineSpeed, Formulas.TorqueToPower(p1.TorqueFullLoad, p1.EngineSpeed));
}
//return null;
var engineTorqueMaxPower = FullLoadStationaryTorque(engineSpeedMaxPower);
return new Tuple<PerSecond, Watt>(engineSpeedMaxPower, Formulas.TorqueToPower(engineTorqueMaxPower, engineSpeedMaxPower));
}
public PerSecond RatedSpeed()
{
var max = new Tuple<PerSecond, Watt>(new PerSecond(), new Watt());
for (var idx = 1; idx < _entries.Count; idx++) {
var currentMax = FindMaxPower(_entries[idx - 1], _entries[idx]);
if (currentMax.Item2 > max.Item2) {
max = currentMax;
}
}
return max.Item1;
}
}
}
\ No newline at end of file
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment

For further information about the platform and support, please see https://code.europa.eu/info/about