Modellbildung und Simulation eines Dieselmotors unter Dymola und ...

diese genutzt. Rein algebraische Gleichungen, die ggf. durch empirische Ansätze entstanden
oder ergänzt worden sind, werden hingegen in c-Code in s-functions gekapselt. in Bild
5-1 ist dies für P MIHP zu sehen.
Bild 5-1: P MIHP in Simulink
Andere Ansätze, wie z.B. alle Komponenten der Kühler, werden durchgängig mit Simulink-
Blöcken umgesetzt. Auch hier könnte man auf Gleichungsebene s-functions eingesetzen; im
vorliegenden Fall wurden diese zur Abgrenzung nur für Modelle mit datenbasiert optimierten
Parametern verwendet.
#define S_FUNCTION_LEVEL 2
#define S_FUNCTION_NAME P_MIHP_v14_dll
//[...] andere s -function Unter-Funktionen
static void mdlOutputs(SimStruct *S, int_T tid){
//Deklarationen etc.[...]
qstern=MF_FUEL*hu/(MF_AE*T_BE*cp);
eta=1-(1/(qstern*kappa))
*((pow((qstern-(1/(kappa*epsylon))
*(x1-(pow(epsylon,kappa)))
+x1/epsylon),kappa))
*(pow(x1,(1-kappa)))-1)
*(1+x2*MF_FUEL_stroke);
P_MIHP=(eta*MF_FUEL*hu)/(I*N_ENGINE*VH);
y[0]= P_MIHP;
}
Bild 5-2: Kühlermodul in Simulink
Wie zu sehen ist, wird die Drossel in zwei verschiedenen Ausführungen benötigt; für die signalorientierte
Darstellung in Simulink bedeutet dies automatisch zwei separate Umsetzungen.
Die Drossel mit dem Massenstrom als Ausgangsgröße kann exakt gemäß Formel (4.4)
umgesetzt werden. Diejenige mit dem Druck als Ausgangsgröße benötigt jedoch eine Umstellung
dieser Gleichung, die rein analytisch nicht möglich ist. Für den vorliegenden Fall
wurde dies relativ einfach gelöst, da ein konstantes Verhältnis der spezifischen Wärmekapazitäten
κ angenommen wurde, so dass der Verlauf von c p zu Beginn der Simulation für dieses
feste κ berechnet und zur Laufzeit tabellarisch invertiert wird. Dies führt zwar zu einer
praktikablen Umsetzung, entfernt sich aber aufgrund der angenommenen Voraussetzungen
von einer verallgemeinerten Einsetzbarkeit.