DE102009007808A1

DE102009007808A1 - Verfahren zum Betrieb einer Verbrennungskraftmaschine

Info

Publication number: DE102009007808A1
Application number: DE102009007808A
Authority: DE
Inventors: Robert Beckmann; Elmar Dr.Rer.Nat. Millich; Wolfgang Dr.-Ing. Drewelow
Original assignee: IAV GmbH Ingenieurgesellschaft Auto und Verkehr
Current assignee: Volkswagen AG
Priority date: 2009-02-04
Filing date: 2009-02-04
Publication date: 2010-08-05
Anticipated expiration: 2029-02-05
Also published as: DE102009007808B4

Abstract

Es ist Aufgabe der vorliegenden Erfindung, für den Betrieb beziehungsweise die Steuerung und Regelung einer Verbrennungskraftmaschine ein Verfahren zur stabilen Berechnung eines Modells einer Verbrennungskraftmaschine mittels eines numerischen Lösungsverfahrens mit einer festen Schrittweite bereitzustellen, so dass auch in Betriebspunkten der Verbrennungskraftmaschine, in denen zwei über eine Trennstelle, wie beispielsweise eine Drosselklappe, ein Turbinen- oder ein Verdichterrad miteinander verbundene Behälter wie ein einziger Behälter wirken, das jeweilige Modell sicher mittels eines gängigen Steuergerätes berechnet werden kann. Diese Aufgabe wird durch die vorliegende Erfindung dadurch gelöst, dass zur stabilen Berechnung eines zum Betrieb einer Verbrennungskraftmaschine dienenden Modells einer Verbrennungskraftmaschine mittels eines numerischen Lösungsverfahrens mit einer festen Schrittweite sämtliche Eigenwerte des linearisierten Modells in Arbeitspunkten der Verbrennungskraftmaschine, die hinsichtlich der Stabilität der Berechnung kritisch sind, in einen stabilen Bereich des numerischen Lösungsverfahrens verschoben werden.

Description

Die vorliegende Erfindung betrifft ein Verfahren zum Betrieb einer Verbrennungskraftmaschine mit den Merkmalen des Patentanspruches 1.
Es ist allgemein bekannt, der Steuerung und Regelung einer Verbrennungskraftmaschine physikalische Modellansätze zu Grunde zu legen, die mit numerischen Lösungsverfahren berechnet werden. Beispielsweise kann die so genannte Füll- und Entleermethode zur Beschreibung des Ansaug- und Abgassystems einer Verbrennungskraftmaschine herangezogen werden. Bei derartigen physikalischen Ansätzen handelt es sich allgemein um nichtlineare zeitkontinuierliche Modelle. Diese Modelle haben von Arbeitspunkt zu Arbeitspunkt eine verschiedene Dynamik und Bandbreite. Verfügt die zu Grunde liegende Verbrennungskraftmaschine beispielsweise über eine in der Ansauganlage angeordnete Drosselklappe zur Einstellung der den Zylindern zuzuführenden Frischluftmasse und ist die Drosselklappe nahezu geschlossen, ändert sich das Druckverhältnis über der Drosselklappe nur sehr langsam, da sich nur ein geringer Massenstrom über der Drosselklappe ausbildet. Ist dem hingegen die Drosselklappe weit geöffnet, kann sich ein großer Massenstrom über der Drosselklappe ausbilden, wodurch sich das Druckverhältnis sehr schnell ändert. Mit anderen Worten ist bei vollständig geöffneter Drosselklappe das Druckverhältnis über der Drosselklappe praktisch eins und die Behälter vor und hinter der Drosselklappe wirken wie ein einziger Behälter. Derartige Effekte können auch bei einer Verbrennungskraftmaschine auftreten, die über einen Abgasturbolader verfügt, wobei sich ebenfalls große Massenströme über dem Verdichter beziehungsweise der Turbine ausbilden und die Behälter vor und hinter dem Verdichter oder der Turbine wie ein einziger Behälter wirken. Zur Berechnung langsamer Ausgleichsvorgänge genügen große Integrationsschrittweiten, wohingegen schnelle Zustandsänderungen kleine Schrittweiten erfordern. Mit anderen Worten kommt es zwischen den zwei Behältern zu sehr schnellen Druckänderungen, für deren Berechnung eine bestimmte feste Schrittweite unter Umständen nicht genügt. Da die zwei genannten Drosselklappenszenarien beziehungsweise Massenstromänderungen über dem Verdichter und der Turbine im Betrieb der Verbrennungskraftmaschine permanent auftreten, kann zur Berechnung des jeweiligen Mo dells ein numerisches Lösungsverfahren mit variabler Schrittweite eingesetzt werden. Dieses ist in der Lage, die Integrationsschrittweite während der Berechnung an das zu berechnende Modell anzupassen. Als Alternative kann ein Verfahren mit einer festen Schrittweite eingesetzt werden, wobei diese so klein zu wählen ist, dass die schnellen Ausgleichsvorgänge stabil berechnet werden können. Das bedeutet jedoch, dass die Berechnung des Modells je nach Bandbreite des Modells vergleichsweise viel Rechenzeit konsumiert. Während die Wahl des Berechnungsverfahrens in einer Entwicklungsumgebung noch relativ frei ist, ist sie das bei einer Implementierung in dem Steuergerät einer Verbrennungskraftmaschine nicht. Hier kann nicht gewartet werden, bis das Verfahren eine Lösung in einem sehr steifen Dynamikbereich berechnet hat. Vielmehr muss zu definierten Zeitpunkten die Lösung vorliegen. Aus diesem Grund ist die Integrationsschrittweite im Steuergerät festgelegt. Für die Berechnung des Ansaugsystems ist beispielsweise eine Schrittweite von 10 ms vorgesehen. Unter den Verfahren mit einer festen Schrittweite gibt es zahlreiche Methoden, die sich hinsichtlich ihres Konvergenz- und Stabilitätsverhaltens zum Teil deutlich unterscheiden. Aufwendige Verfahren besitzen bei gleicher Schrittweite häufig ein besseres Stabilitätsverhalten als einfache Ansätze, benötigen dafür jedoch auch mehr Rechenschritte und damit mehr Rechenzeit. Da diese im Steuergerät einer Verbrennungskraftmaschine sehr begrenzt ist, ist es erforderlich, ein einfaches numerisches Verfahren zur näherungsweisen Lösung zu nutzen, beispielsweise das so genannte Euler-Verfahren zur numerischen Lösung von Differentialgleichungssystemen. Praktische Versuche haben jedoch gezeigt, dass, wenn bei einem in einer Entwicklungsumgebung implementierten Modell einer Verbrennungskraftmaschine das Lösungsverfahren von einem Verfahren mit variabler Schrittweite auf das Euler-Verfahren mit einer festen Integrationsschrittweite von 10 ms gewechselt wird, die Simulation in den meisten der Arbeitspunkte instabil ist.
Aufgabe
Es ist Aufgabe der vorliegenden Erfindung, für den Betrieb beziehungsweise die Steuerung und Regelung einer Verbrennungskraftmaschine ein Verfahren zur stabilen Berechnung eines linearisierten Modells einer Verbrennungskraftmaschine mittels eines numerischen Lösungsverfahrens mit einer festen Schrittweite bereitzustellen, so dass auch in Betriebspunkten der Verbrennungskraftmaschine, in denen zwei über eine Trennstelle wie beispielsweise eine Drosselklappe, ein Turbinen- oder ein Verdichterrad miteinander verbundene Behälter wie ein einziger Behälter wirken, das jeweilige Modell sicher mittels eines gängigen Steuergerätes berechnet werden kann.
Lösung
Diese Aufgabe wird durch die vorliegende Erfindung dadurch gelöst, dass zur stabilen Berechnung eines zum Betrieb einer Verbrennungskraftmaschine dienenden Modells einer Verbrennungskraftmaschine mittels eines numerischen Lösungsverfahrens mit einer festen Schrittweite sämtliche Eigenwerte des linearisierten Modells in Arbeitspunkten der Verbrennungskraftmaschine, die hinsichtlich der Stabilität der Berechnung kritisch sind, in einen stabilen Bereich des numerischen Lösungsverfahrens verschoben werden. Die Erfindung baut auf der Erkenntnis auf, dass die Eigenwerte eines linearisierten Modells beziehungsweise eines linearisierten Systems, bestehend aus zwei durch eine Trennstelle gekoppelten Behältern insbesondere durch den durch die Trennstelle erzeugten Eigenwertbeitrag bestimmt werden. Eine Trennstelle zwischen zwei Behältern ist dabei bevorzugt eine Drosselklappe oder ein Verdichter in der Ansaugleitung oder eine Turbine in der Abgasleitung einer Verbrennungskraftmaschine. Eine Entscheidung, ob das numerische Lösungsverfahren stabil ist, erfolgt bevorzugt anhand einer Untersuchung des Stabilitätskriteriums für zeitdiskrete Systeme in Zustandsdarstellung, wie im Ausführungsbeispiel näher beschrieben ist. Im Sinne der vorliegenden Erfindung entspricht der insbesondere von dem Massenstrom über eine Trennstelle zwischen zwei Behältern erzeugte Eigenwertbeitrag mindestens einem Eigenwertbeitrag mehrerer möglicher Eigenwertbeiträge der Systemmatrix eines linearen Differentialgleichungssystems, wobei das numerische Lösungsverfahren dann stabil arbeitet, wenn alle Eigenwerte der Systemmatrix eine von der Schrittweite des numerischen Lösungsverfahrens abhängige Bedingung erfüllen. Das Differentialgleichungssystem beschreibt dabei zumindest ein Teilmodell einer Verbrennungskraftmaschine und die Systemmatrix wird durch Linearisierung des zu Grunde liegenden Differentialgleichungssystems gebildet, wie im Ausführungsbeispiel näher beschrieben ist. Die Verschiebung von mindestens einem Eigenwert, der insbesondere durch den Eigenwertbeitrag des Massenstroms über eine Trennstelle zwischen zwei Behältern bestimmt ist, in einen stabilen Bereich erfolgt erfindungsgemäß mittels eines Skalierungsfaktors. Insbesondere beschreibt der Skalierungsfaktor einen Korrekturmassen strom. In einer Ausführung der vorliegenden Erfindung wird der Skalierungsfaktor als Funktion des Eigenwertbeitrages des Massenstromes über die jeweilige Trennstelle gewählt, wobei in einer besonders vorteilhaften Ausführung der vorliegenden Erfindung das Kriterium dabei ist, dass der Eigenwertbeitrag gleich einem Grenzwert ist, wobei nur dann in das zugrunde liegende Modell eingegriffen wird, wenn der Eigenwertbeitrag des Massenstroms über die jeweilige Trennstelle kleiner als dieser Grenzwert ist.
Durch die vorliegende Erfindung wird demnach der Vorteil erreicht, dass Modelle zur Beschreibung des Betriebsverhaltens einer Verbrennungskraftmaschine, die insbesondere zur Steuerung und Regelung einer Verbrennungskraftmaschine verwendet werden, mit numerischen Lösungsverfahren, wie dem Euler-Verfahren und einer festen Schrittweite von beispielsweise 10 ms berechnet und in dem Steuergerät einer Verbrennungskraftmaschine echtzeitfähig implementiert werden können. Da kritische Komponenten, wie die Drosselklappe oder der Verdichter in der Ansaugleitung beziehungsweise die Turbine in der Abgasleitung in einigen Arbeitspunkten zu große Eigenwertbeiträge in den umgebenden Zuständen liefern, ist nicht durchgängig eine echtzeitfähige Berechnung der dazugehörigen Modelle, beispielsweise mit dem Euler-Verfahren, möglich. Erfindungsgemäß vorteilhaft wird ein Ansatz zur Beeinflussung der Dynamik in den betreffenden Volumina und Arbeitspunkten angewendet. Dieser Ansatz erzeugt keinen stationären Fehler und insbesondere bei der Drosselklappe nur einen sehr kleinen dynamischen Fehler. Insgesamt wird mit diesen Maßnahmen ein echtzeitfähiges Modell erstellt, welches keinen zusätzlichen stationären Fehler erzeugt.
Ausführungsbeispiel
Weitere vorteilhafte Ausgestaltungen der vorliegenden Erfindung sind dem nachfolgenden Ausführungsbeispiel sowie den abhängigen Patentansprüchen zu entnehmen.
Hierbei zeigen:
1: Modell der Ansauganlage einer Verbrennungskraftmaschine,
2: Darstellung der Eigenwerte der Systemmatrix,
3: Darstellung einer Funktionsstruktur,
4: Darstellung eines Details der Funktionsstruktur.
Ein allgemein bekannter physikalischer Modellansatz zur Steuerung und Regelung einer Verbrennungskraftmaschine beschreibt die Ansauganlage. Wie in 1 dargestellt, beschreibt ein einfaches Modell der Ansauganlage einen ersten Behälter 1 und einen zweiten Behälter 2, die durch eine Drosselklappe 3 voneinander getrennt sind beziehungsweise durch das Schließen der Drosselklappe 3 voneinander getrennt werden können. Dem Behälter 2 fließt über einen Luftfilter 4 Umgebungsluft m ._LF zu, wie durch den Pfeil angedeutet. In dem Behälter 2 herrscht die Temperatur T_2i und der Druck p_2i und in dem Behälter 1 die Temperatur T_1i und der Druck p_1i. Vereinfachend wird bei diesem isothermen Modellansatz davon ausgegangen, dass die Temperaturen T_1i und T_2i in den Behältern 1 und 2 sowie die Umgebungstemperatur gleiche Werte haben und konstant sind. Außerdem wird angenommen, dass die Temperatur T_1i, T_2i in beiden Behältern 1 und 2 gleich der Temperatur der einströmenden Umgebungsluft ist und es zu keinem Wärmeaustausch mit der Umgebung kommt. Über die Drosselklappe 3 fließt dem Behälter 1 ein Massenstrom m ._DK zu. Aus dem Behälter 1 fließt weiterhin ein Massenstrom m ._EV zu den Brennräumen der Verbrennungskraftmaschine aus dem Behälter 1 ab. Aufbauend auf dem allgemeinen Gasgesetz kann unter den genannten Annahmen und definierten Größen die Zustandsdifferentialgleichung des Druckes p_1i im Behälter 1 als Funktion der Massenströme m ._DK und m ._EV gemäß Gleichung (1) dargestellt werden.
Der Massenstrom m ._DK über die Drosselklappe 3 kann wie allgemein bekannt gemäß Gleichung (2) dargestellt werden, wobei A der durchströmten Fläche als Funktion des Drosselklappenwinkels α_DK und Ψ der Durchflussfunktion entspricht, welche die Wirkung der Drosselklappe 3 in Abhängigkeit des Druckverhältnisses über der Drosselklappe 3 beschreibt.
Die Herleitung der Durchflussfunktion Ψ ist beispielsweise der Literaturstelle Merker, G. P.; Schwarz, C.; Stiesch, G.: Verbrennungsmotoren. Simulation der Verbrennung und Schadstoffbildung. Teuber, 2004 zu entnehmen.
Der Massenstrom m ._EV, der zu den Brennräumen aus dem Behälter 1 abfließt, kann weiterhin gemäß Gleichung (3) beschrieben werden, wobei K_Mot einer motorspezifischen Konstante mit der Einheit
und n_Mot der Drehzahl der Kurbelwelle der Verbrennungskraftmaschine entspricht. Die beiden Koeffizienten c_EV,1 und c_EV,2 beschreiben ferner den im Wesentlichen linearen funktionalen Zusammenhang zwischen dem Druck p_1i im Behälter 1 und einer so genannten relativen Füllung, wobei die relative Füllung als Quotient der aktuellen Luftfüllung zu einer Luftfüllung unter bestimmten Normbedingungen, multipliziert mit 100%, beschrieben wird. m .EV = KMot·nMot·(cEV,1·p1i + cEV,2) (3)
Zusätzlich wird der Behälter 2 gemäß Gleichung (4) betrachtet.
Der Behälter 2 wird über die Drosselklappe 3 entleert und über einen linearen Strömungswiderstand gegen den Druck der Umgebungsluft p_Umg abgegrenzt, wobei der Strömungswiderstand als linearisierte Darstellung des Luftfilters 4 gemäß Gleichung (5) interpretiert wird. m .LF = cLF(pUmg – p2i) (5)
Der Koeffizient c_LF beschreibt dabei den Strömungswiderstand, der durch den Luftfilter 4 bedingt ist.
Zum weiteren Verständnis wird weiterhin vorerst auf das Euler-Verfahren zur numerischen Lösung von gewöhnlichen Differentialgleichungssystemen erster Ordnung der Form gemäß Gleichung (5.1) eingegangen,
das mit der Anfangsbedingung gemäß (5.2)
beschrieben werden kann. Dafür werden (5.1) und (5.2) zunächst in die Integralform gemäß (5.3) überführt.
Die Lösung von (5.3) zum Zeitpunkt t₀ + Δt entspricht demnach wie (5.4) in dargestellt.
Da der Integrand in vielen Fällen nicht numerisch integrierbar ist, wird angenommen, dass dieser im Intervall [t₀, t₀ + Δt] für eine genügend kleine Integrationsschrittweite Δt gemäß (5.5) beziehungsweise (5.6) konstant ist. f(x(s)) = f(x0) (5.5)
Aus der Gleichung für den ersten Schritt (5.6) wird die Rekursionsgleichung (5.7) für die approximierte Lösung von (5.1) und (5.2) im gesamten Intervall [t₀, T] abgeleitet. x(ti + tΔ) = x(ti) + Δt·f(x(ti)) (5.7)
Für eine bessere Übersicht wird im folgenden die Schreibweise gemäß (5.8) verwendet. xi+1 = xi + Δt·f(xi) (5.8)
Die Konvergenz des Euler-Verfahrens hängt wesentlich von der gewählten Schrittweite Δt ab. Kommt es innerhalb eines Integrationsschrittes zu ausgeprägten Änderungen in den Zuständen x, ist das Euler-Verfahren nicht in der Lage diese abzubilden. Nicht nur die Konvergenz des Verfahrens, sondern auch die Stabilität hängt von der gewählten Schrittweite Δt ab. Ist die Schrittweite zur Berechnung einer bestimmten Eigenbewegung in (5.1), (5.2) zu groß, kommt es zu einem exponentiellen Wachstum des Verfahrensfehlers und das Euler-Verfahren wird instabil. Die Stabilität des Euler-Verfahrens in Verbindung mit dem zu Grunde liegenden Modell der Verbrennungskraftmaschine wird nur in ausgewählten Arbeitspunkten des Systems nachgewiesen. Ist das Euler-Verfahren in allen am Prüfstand gemessenen Arbeitspunkten stabil und sind diese eng genug gewählt, wird angenommen, dass auch alle Arbeitspunktübergänge stabil berechnet werden können und das Verfahren im interessierenden Arbeitsbereich insgesamt stabil arbeitet. Als Arbeitsbereich wird dabei der von den Arbeitspunkten abgedeckte Raum verstanden. Die Verfahrensstabilität selbst wird anhand des Stabilitätskriteriums für zeitdiskrete Systeme in Zustandsdarstellung untersucht.
Zur Beschreibung der Systemdynamik wird die Änderung der Systemzustände f(x(t)) in (5.1) im Arbeitspunkt x_R linearisiert. Die Bezeichnung Arbeitspunkt meint dabei, dass sich das System in diesem Punkt in einer Ruhelage befindet, die Zustände somit konstant sind. Die Reihenentwicklung von f(x) gemäß (5.9)
wird nach dem Linearglied abgebrochen, vergleiche (5.10).
Die Linearisierung von (5.1) soll die Systemdynamik nur in der Umgebung von x_R beschreiben, da x_R eine Ruhelage des Systems (5.1) ist, gilt f(x_R) = 0, wodurch sich (5.10) gemäß (5.11) vereinfacht.
Die linearisierte Darstellung von (5.1) im Arbeitspunkt x_R entspricht damit der Darstellung gemäß (5.12).
Üblicherweise wird die Verschiebung x(t) – x_R in der Schreibweise vernachlässigt, so dass die im Arbeitspunkt linearisierte Darstellung des Systems (5.12) mit ddt xR = 0 gemäß (5.13)
beziehungsweise gemäß (5.14)
lautet. Damit kann die Gleichung des Euler-Verfahrens (5.8) gemäß (5.15) beziehungsweise gemäß (5.16) und mit Φ = (I_n + Δt·A) gemäß (5.17) geschrieben werden.
xi+1 = (In + Δt·A)xi (5.16) xi+1 = Φxi (5.17)
Das System (5.17) ist stabil, wenn alle Eigenwerte λ_Φ,k der Matrix Φ die Bedingung gemäß (5.18) erfüllen, wie beispielsweise der Literaturstelle Lunze, D. I. J.: Regelungstechnik 2. Bd. 2. 2. Springer, 2002 zu entnehmen ist. |λΦ,k| ≤ 1 (k = 1, 2, ..., n) (5.18)
Daraus folgt, dass die Eigenwerte λ_A,k der Systemmatrix A für die Stabilität von (5.15) der Bedingung gemäß (5.19) genügen müssen. |λA,k·Δt + 1| ≤ 1 (k = 1, 2, ..., n) (5.19)
Mit diesen Ungleichungen (5.18, 5.19) wird der Zusammenhang zwischen der Systemdynamik des Modells und der für eine stabile Berechnung erforderlichen Abtastzeit formuliert. Um eine Aussage zu treffen, ob ein gegebenes System (5.1, 5.2) mit der Abtastzeit Δt im Ruhepunkt x_R beziehungsweise einer gewissen Umgebung um diesen Punkt stabil mit dem Euler-Verfahren berechnet werden kann, wird das System in diesem Arbeitspunkt linearisiert, die Systemmatrix A berechnet und anhand der Ei genwerte λ_A,k der Systemmatrix A entschieden, ob die Ungleichungen aus (5.19) erfüllt sind. Aus (5.19) folgt, dass die reellen Eigenwerte λ_A,k der Systemmatrix A bei einer Abtastzeit von 10 ms im Bereich gemäß (5.20) liegen müssen. –200 s–1 ≤ λA,k ≤ 0 s–1 (5.20)
Sind die Eigenwerte komplex, muss (5.19) komplexwertig ausgewertet werden.

Um den Einfluss des Massenstroms m ._DK über die Drosselklappe 3 auf die Dynamik des jeweiligen Modells beispielhaft zu beschreiben, werden alle Prozessgrößen auf Messungen gesetzt, die in einem Arbeitspunkt unter Einsatz eines Motorprüfstandes und geeigneter messtechnischer Mittel bestimmt wurden. Dabei wird bevorzugt ein Arbeitspunkt gewählt, bei dem das Druckverhältnis p_1i/p_2i über die Drosselklappe 3 nahe eins ist. Die Parameter des Modells und die Prozessgrößen in dem gemessenen Arbeitspunkt sind gemäß Tabelle (a) dargestellt. Tabelle (a)

p_1i = 1.2921·10⁵ Pa	T_1i = 300 K	V_1i = 2 dm³	n_Mot = 1500 1/min
p_2i = 1.3133·10⁵ Pa	T_2i = 300 K	V_2i = 1 dm³	α_DK = 34%
A(α_DK) = 3.5·10^–4 m²	R_Luft = 287 J/K·kg	c_LF = 1.5·10^–7 kg/s·Pa	K_Mot = 2.155·10^–7
c_EV,1 = 7.4051·10^–4	c_EV,2 = –4.8540	κ = 1.4

Zur weiteren Beschreibung des Einflusses des Massenstroms m ._DK über die Drosselklappe 3 auf die Dynamik des jeweiligen Modells wird im weiteren Verlauf zunächst die Dynamik der Differentialgleichung gemäß (1) untersucht. Dazu wird gemäß (6, 6.1, 6.2) eine Systemmatrix A von Gleichung (1) durch Linearisierung gebildet.
Der Eigenwertbeitrag λ_EV des Massenstroms m ._EV durch die Einlassventile kann daher gemäß (7, 7.1, 7.2) auf Grundlage von Gleichung (3) bestimmt werden.
Der Eigenwertbeitrag λ_EV ist betragsmäßig klein und hängt nicht vom Druck p_1i im Behälter 1 oder vom Druck p_2i im Behälter 2 ab.
Der Eigenwertbeitrag λ_DK des Massenstroms m ._DK über die Drosselklappe 3 hängt demgegenüber durch die Durchflusskennlinie
sehr stark vom Druckverhältnis über die Drosselklappe 3 ab, wie gemäß (8, 8.1, 8.2) dargestellt.
Ursache dafür ist, dass die Ableitung der Durchflusskennlinie
theoretisch für Druckverhältnisse nahe eins gegen minus unendlich geht. Das bedeutet, dass auch der Eigenwertbeitrag λ_DK des Massenstroms m ._DK über die Drosselklappe 3 gegen minus unendlich geht.
Wird das Modell gemäß den Gleichungen (1) und (4) betrachtet, so wird dessen Dynamik durch die Eigenwerte λ₁ und λ₂ der Systemmatrix A gemäß (9, 9.1) charakterisiert.
Es wird deutlich, dass das Modell bestehend aus den Behältern 1 und 2 und der Drosselklappe 3 mit dem Euler-Verfahren und einer Schrittweite von 10 ms nicht berechnet werden kann. Dieses Ergebnis kann auf ein übergeordnetes Gesamtmodell übertragen werden.
Um das Modell gemäß (1) und (4) mit einer Schrittweite von 10 ms berechnen zu können, wird erfindungsgemäß ein Skalierungsfaktor γ eingeführt, Gleichung (10).
Anschließend wird ein zusätzlicher Korrekturmassenstrom m ._DKK eingeführt. Weiterhin wird ein Zusammenhang zwischen dem Skalierungsfaktor γ und dem Korrekturmassenstrom m ._DKK gesucht, Gleichung (11).
Werden die Gleichungen (10) und (11) gleichgesetzt, kann m ._DKK durch den Skalierungsfaktor γ ausgedrückt werden, Gleichungen (12) und (13).
Der Skalierungsfaktor γ wird als Funktion des Eigenwertbeitrages λ_DK, der nach Gleichung (8) definiert sei, gewählt. Kriterium ist dabei, dass der Eigenwertbeitrag λ_DK gleich einem Grenzwert λ_G sein soll, Gleichung (14).
Gleichung (14) wird nach dem Skalierungsfaktor γ umgestellt, Gleichung (15).
Weiterhin wird erfindungsgemäß nur dann in die Dynamik eingegriffen, wenn der Eigenwertbeitrag λ_DK des Drosselklappenmassenstromes m ._DK kleiner als der Grenzwert λ_G ist. Die Bestimmungsgleichung für den Skalierungsfaktor γ lautet daher gemäß Gleichung (16).
Ist also der Eigenwertbeitrag λ_DK des Drosselklappenmassenstromes m ._DK größer als der Grenzwert λ_G, bleibt die Dynamik mit γ = 1 unverändert. Erst wenn der Eigenwertbeitrag λ_DK kleiner als der Grenzwert λ_G ist, erfolgt ein Eingriff in die Dynamik des Modells.
Um einen insgesamt geschlossenen Ansatz zu erhalten, wird der Korrekturmassenstrom m ._DKK auch in die Massenbilanz des Behälters 2 aufgenommen. Insgesamt besitzt das gekoppelte System dadurch die Form gemäß Gleichungen (17) bis (20).
Um zu untersuchen, wie der beschriebene Ansatz die Dynamik des Modells im geschlossenen Kreis beeinflusst, wird das System gemäß den Gleichungen (17) und (18) linearisiert und es werden die Eigenwerte λ₁ und λ₂ der Systemmatrix A als Funktion des Grenzwertes λ_G berechnet, Gleichung (21).
In 2 sind die beiden Eigenwerte λ₁ und λ₂ der Systemmatrix A über dem Grenzwert λ_G dargestellt. Es ist zu erkennen, dass ein Eigenwert ab λG = –10 1s nahe rungsweise konstant –10 1 / s ist. Der zweite Eigenwert wird hingegen ab λG = –10 1s annähernd linear mit λ_G skaliert. Durch die Kopplung der beiden Teilsysteme (Behälter 2 und 1 über die Drosselklappe 3) legt λ_G keinen Eigenwert direkt fest, sondern wirkt als Skalierungsfaktor auf den kritischen Eigenwert des Systems.
3 zeigt weiterhin eine Funktionsstruktur 5, die im Rahmen des Betriebs, also der Steuerung und Regelung einer Verbrennungskraftmaschine, den Massenstrom m ._DK' über die Drosselklappe 3 berechnet. Das Hochkomma ' symbolisiert dabei, dass es sich dabei um einen mittels des erfindungsgemäßen Verfahrens ermittelten Wert handelt, da bei dessen Berechnung mittels eines numerischen Lösungsverfahrens in kritischen Arbeitspunkten eine Verschiebung der Eigenwerte der Systemmatrix des linearisierten Systems gemäß Gleichung (1) und (4) in einen stabilen Bereich des numerischen Lösungsverfahrens erfolgt. In einem ersten Block 6 der Funktionsstruktur 5 erfolgt, wie allgemein bekannt, die Berechnung des Massenstroms m ._DK über die Drosselklappe 3 gemäß Gleichung (2). Dazu werden Block 6 die in Gleichung (2) einfließenden Größen zugeführt. In einem weiteren Block 7 erfolgt die erfindungsgemäße Berechnung des Korrekturmassenstroms m ._DKK. Mit dem Block 7 ist weiterhin ein Block 8 verbunden. In Block 8 erfolgt in Abhängigkeit des jeweils vorliegenden Druckverhältnisses p_1i/p_2i über die Drosselklappe 3 eine Entscheidung, ob eine Berechnung des Korrekturmassenstroms m ._DKK erfolgen soll oder nicht. Soll keine solche Berechnung erfolgen, wird der Korrekturmassenstrom m ._DKK auf null gesetzt. Erfolgt jedoch eine Berechnung des Korrekturmassenstroms m ._DKK, werden der gemäß Gleichung (2) berechnete Massenstrom m ._DK und der Korrekturmassenstrom m ._DKK in Block 9 addiert, so dass der Massenstrom m ._DK' über die Drosselklappe 3 einer weiteren Verarbeitung bereitsteht. Das Resultat einer ersten Berechnung des Massenstroms m ._DK über die Drosselklappe 3 gemäß Gleichung (2) wird wie in 3 gezeigt, dem Block 7 zugeführt, wie auch weitere Größen, die im folgenden näher beschrieben werden.
4 zeigt den Block 7 im Detail. Block 7 umfasst zunächst einen weiteren Block 10, in dem gemäß Gleichung (8.1) die Berechnung des Eigenwertbeitrages λ_DK erfolgt. Dazu werden Block 10 die in Gleichung (8.1) einfließenden Größen zugeführt. Im wei teren Verlauf wird ein Quotient aus dem in Block 10 berechneten Eigenwertbeitrag λ_DK und dem Grenzwert λ_G beziehungsweise der Skalierungsfaktor γ gebildet. Der Skalierungsfaktor γ wird einem Block 11 zugeführt. Dem Block 11 wird außerdem der Zahlenwert 1 zugeführt. Block 11 umfasst eine Funktion, die den Eingang weitergibt, der jeweils der kleinere ist. Ist demnach der Skalierungsfaktor γ kleiner 1, wird der Skalierungsfaktor γ weitergegeben. Ist der Skalierungsfaktor γ größer als 1, wird der Wert 1 weitergegeben. In einem noch weiteren Block 12 wird von dem jeweils von Block 11 weitergegebenen Wert der Wert 1 abgezogen beziehungsweise die Differenz (γ - 1) gebildet. In einem noch weiteren Block 13 werden die Differenz (γ – 1) und die Differenz (m ._DK – m ._EV) multipliziert, so dass Gleichung (19) zur Berechnung des Korrekturmassenstroms m ._DKK gebildet wird.

1: Behälter
2: Behälter
3: Drosselklappe
4: Luftfilter
5: Funktionsstruktur
6–13: Blöcke der Funktionsstruktur

ZITATE ENTHALTEN IN DER BESCHREIBUNG
Diese Liste der vom Anmelder aufgeführten Dokumente wurde automatisiert erzeugt und ist ausschließlich zur besseren Information des Lesers aufgenommen. Die Liste ist nicht Bestandteil der deutschen Patent- bzw. Gebrauchsmusteranmeldung. Das DPMA übernimmt keinerlei Haftung für etwaige Fehler oder Auslassungen.
Zitierte Nicht-Patentliteratur

- Merker, G. P.; Schwarz, C.; Stiesch, G.: Verbrennungsmotoren. Simulation der Verbrennung und Schadstoffbildung. Teuber, 2004 [0014]
- Lunze, D. I. J.: Regelungstechnik 2. Bd. 2. 2. Springer, 2002 [0029]

Claims

Verfahren zum Betrieb einer Verbrennungskraftmaschine, wobei ein zur Steuerung und Regelung der Verbrennungskraftmaschine dienendes Modell der Verbrennungskraftmaschine mittels eines numerischen Lösungsverfahrens mit einer festen Schrittweite berechnet wird, wobei sämtliche Eigenwerte (λ) des linearisierten Modells in Arbeitspunkten der Verbrennungskraftmaschine, die hinsichtlich der Stabilität der Berechnung kritisch sind, in einen stabilen Bereich des numerischen Lösungsverfahrens verschoben werden.
Verfahren nach Patentanspruch 1, wobei die Entscheidung, ob das numerische Lösungsverfahren stabil ist, anhand einer Untersuchung des Stabilitätskriteriums für zeitdiskrete Systeme in Zustandsdarstellung erfolgt.
Verfahren nach Patentanspruch 1 oder 2, wobei die Verschiebung der Eigenwerte (λ) in einen stabilen Bereich des numerischen Lösungsverfahrens mittels eines Skalierungsfaktors (γ) erfolgt.
Verfahren nach Patentanspruch 3, wobei der Skalierungsfaktor (γ) einen Korrekturmassenstrom (m ._DKK) beschreibt.
Verfahren nach Patentanspruch 3 oder 4, wobei der Skalierungsfaktor (γ) als Funktion des Eigenwertbeitrages (λ_DK) des Massenstromes (m ._DK) über die jeweilige Trennstelle (3) gewählt wird, wobei das Kriterium dabei ist, dass der Eigenwertbeitrag (λ_DK) gleich einem Grenzwert (λ_G) ist, wobei nur dann in die Dynamik eingegriffen wird, wenn der Eigenwertbeitrag (λ_DK) kleiner als der Grenzwert (λ_G) ist.