From 078719c134f04cd06a93ea954bad574d4e47cb09 Mon Sep 17 00:00:00 2001
From: Raphael Luz <luz@ivt.tugraz.at>
Date: Thu, 29 Jan 2015 12:07:40 +0100
Subject: [PATCH] Cleaner engine / clutch power calculation
---
VECTO/MODcalc/cPower.vb | 64 ++++++++++++++++++++---------------------
1 file changed, 31 insertions(+), 33 deletions(-)
diff --git a/VECTO/MODcalc/cPower.vb b/VECTO/MODcalc/cPower.vb
index 4b8497036e..76a5cb3d32 100644
--- a/VECTO/MODcalc/cPower.vb
+++ b/VECTO/MODcalc/cPower.vb
@@ -219,27 +219,27 @@ Public Class cPower
PaMot = fPaMot(nU, LastnU)
End If
- 'Aux Demand
- Paux = fPaux(i, nU)
- 'Engine Power (at Clutch)
- If Pplus Or Pminus Then
+ 'Calculate powertrain losses => power at clutch
+ If Pplus Or Pminus Then
PlossGB = fPlossGB(Pwheel, Vact, Gear, True)
PlossDiff = fPlossDiff(Pwheel, Vact, True)
PlossRt = fPlossRt(Vact, Gear)
PaGetr = fPaG(Vact, aact)
-
Pkup = Pwheel + PlossGB + PlossDiff + PaGetr + PlossRt
- P = Pkup + Paux + PaMot
-
Else
-
Pkup = 0
- P = Paux + PaMot
-
End If
+ 'Total aux power
+ '[kW]
+ Paux = fPaux(i, nU)
+
+ 'Internal Engine Power (Pclutch plus Aux plus Inertia)
+ P = Pkup + Paux + PaMot
+
+
'Full load / motoring
Pmin = FLD(Gear).Pdrag(nU)
@@ -933,26 +933,9 @@ lb_nOK:
'************************************ Determine Engine-state ************************************
' nU is final here!
- 'Determine Aux power demand (from VEH and DRI)
- Paux = fPaux(jz, Math.Max(nU, ENG.Nidle))
-
- 'ICE-inertia
- If jz = 0 Then
- PaMot = 0
- Else
- 'Not optimal since jz-1 to jz not the right interval
- PaMot = fPaMot(nU, MODdata.nU(jz - 1))
- End If
-
-
-
- 'Total Engine-power
- ' => Pantr
- ' => P
- ' => Pkup
+ 'Power at clutch
Select Case Clutch
Case tEngClutch.Opened
- P = Paux + PaMot
Pclutch = 0
PlossGB = 0
PlossDiff = 0
@@ -963,7 +946,7 @@ lb_nOK:
If GBX.TCon And GBX.IsTCgear(Gear) Then
- P = nMtoPe(nU, GBX.TCMin) + Paux + PaMot
+ Pclutch = nMtoPe(nU, GBX.TCMin)
If P >= 0 Then
PlossTC = Math.Abs(nMtoPe(GBX.TCnUin, GBX.TCMin) * (1 - GBX.TC_mu * GBX.TC_nu))
@@ -971,7 +954,6 @@ lb_nOK:
PlossTC = Math.Abs(nMtoPe(GBX.TCnUout, GBX.TCMout) * (1 - GBX.TC_mu * GBX.TC_nu))
End If
-
Else
PlossGB = fPlossGB(Pwheel, Vact, Gear, False)
@@ -980,7 +962,6 @@ lb_nOK:
PlossTC = 0
PaGbx = fPaG(Vact, aact)
Pclutch = Pwheel + PlossGB + PlossDiff + PaGbx + PlossRt
- P = Pclutch + Paux + PaMot
End If
Case Else 'tEngClutch.Slipping: never in AT mode!
@@ -990,9 +971,26 @@ lb_nOK:
PlossTC = 0
PaGbx = fPaG(Vact, aact)
Pclutch = (Pwheel + PlossGB + PlossDiff + PaGbx + PlossRt) / ClutchEta
- P = Pclutch + Paux + PaMot
End Select
+
+ 'Total aux power
+ '[kW]
+ Paux = fPaux(jz, Math.Max(nU, ENG.Nidle))
+
+
+ 'ICE-inertia
+ If jz = 0 Then
+ PaMot = 0
+ Else
+ 'Not optimal since jz-1 to jz not the right interval
+ PaMot = fPaMot(nU, MODdata.nU(jz - 1))
+ End If
+
+ 'Internal Engine Power (Pclutch plus Aux plus Inertia)
+ P = Pclutch + Paux + PaMot
+
+
'EngState
If Clutch = tEngClutch.Opened Then
@@ -2295,7 +2293,7 @@ lb10:
Return CSng(((VEH.Loading + VEH.Mass + VEH.MassExtra) * 9.81 * Math.Sin(Math.Atan(Grad * 0.01)) * v) * 0.001)
End Function
- '----------------Ancillaries(Nebenaggregate) ----------------
+ '----------------Auxillaries(Nebenaggregate) ----------------
Public Function fPaux(ByVal t As Integer, ByVal nU As Single) As Single
Return CSng(MODdata.Vh.Padd(t) + VEC.PauxSum(t, nU))
End Function
--
GitLab