AT522290B1 - Verfahren und Regelungseinheit zur Regelung eines nichtlinearen technischen Prozesses - Google Patents

Abstract

Für eine effiziente und mit wenig Rechenaufwand lösbare modellprädiktive Regelung ist vorgesehen, dass jedes als lineares Zustandsraummodell mit gleicher Systemordnung vorliegende Teilmodell (Mi) in die Regelungsnormalform transformiert wird, wobei jedes transformierte Teilmodell (Mi) dieselbe zumindest eine virtuelle Eingangsgröße v und denselben virtuellen Zustandsvektor z aufweist, sodass mit der modellprädiktiven Regelung (MPC) der zeitliche Verlauf der zumindest einen, für alle transformierten Teilmodelle (Mi) gültigen, virtuellen Eingangsgröße v über einen vorgegebenen Prädiktionshorizont NP zu berechnen, aus dem dann die zumindest eine Eingangsgröße (u) des technischen Prozesses (1) des aktuellen Zeitschrittes berechnet wird, die als Stellgröße für die Regelung des technischen Prozesses (1) verwendet wird.

Description

Beschreibung

VERFAHREN UND REGELUNGSEINHEIT ZUR REGELUNG EINES NICHTLINEAREN TECHNISCHEN PROZESSES

[0001] Die gegenständliche Erfindung betrifft ein Verfahren zur Regelung eines nichtlinearen technischen Prozesses mit zumindest einer Eingangsgröße und zumindest einer Ausgangsgröße, mit einer modellprädiktiven Regelung, wobei mit der modellprädiktiven Regelung in jedem Zeitschritt der Regelung die zumindest eine Eingangsgröße des nichtlinearen technischen Prozesses berechnet wird, die dem technischen Prozess im aktuellen Zeitschritt der Regelung zur Regelung aufgeschaltet wird, und der nichtlineare technische Prozess mit einem Modell modelliert wird, das die zumindest eine Eingangsgröße auf eine Schätzung der zumindest einen Ausgangsgröße abbildet, wobei das Modell aus einer Anzahl von Teilmodellen besteht und die Ausgangsgröße des Modells aus einer gegebenen Kombinationsvorschrift der Modellausgänge der Teilmodelle gebildet wird. Die Erfindung betrifft ebenso eine entsprechende Regelungseinheit zur Regelung eines nichtlinearen technischen Prozesses.

[0002] In technischen Maschinen und Anlagen gibt es üblicherweise eine Fülle von technischen Funktionen oder Prozessen, die für den Betrieb der Maschine oder Anlage geregelt werden müssen. Regeln bedeutet in diesem Zusammenhang in der Regel, dass zumindest eine Ausgangsgröße der Maschine oder Anlage vorgegebenen Sollwerten mit möglichst geringem Fehler und mit möglichst guter Dynamik folgen sollen. Der Fehler ist dabei die Abweichung zwischen den Sollwerten und den Istwerten der Ausgangsgröße. Die Dynamik der Regelung bestimmt im Wesentlichen wie gut die Ausgangsgröße bei raschen zeitlichen Anderungen der Sollwerte (z.B. bei einem Sprung oder einem steilen Anstieg) den Sollwerten folgen kann. Dabei werden von der Regelung gewisse Stellgrößen (Eingangsgrößen der Maschine oder Anlage) berechnet, die verwendet werden, um über Aktuatoren an der Maschine oder Anlage Regeleingriffe an der Maschine oder Anlage vorzunehmen. Die Regelung verarbeitet dabei in der Regel auch Istwerte der Ausgangsgröße. Beispielsweise wird bei der Regelung eine Abweichung zwischen den Sollwerten und den Istwerten ermittelt, um die Abweichung mit der Regelung auszuregeln.

[0003] Ein Beispiel für eine technische Maschine ist ein Kraftfahrzeug, in dem eine Vielzahl von geregelten technischen Funktionen realisiert sind, beispielsweise eine Motorsteuerung, eine Getriebesteuerung, die Steuerung eines Hybridantriebsstranges, ein Batteriemanagementsystem, die Regelung von Hilfsaggregaten des Fahrzeugs, eine Wärmerückgewinnung aus Abgasen (Waste Heat Recovery), die Regelung einer Abgasrückführung, die Regelung eines Turboladers, usw. Ein anderes Beispiel ist eine Brennstoffzelle mit einer Regelung der Konditionierung (z.B. Druck, Temperatur, Feuchtigkeit, Massenstrom) eines Prozessgases (z.B. Luft oder Sauerstoff). Daneben gibt es natürlich noch unzählige weitere Beispiele für eine technische Maschine oder Anlage mit zu regelnden technischen Funktionen / Prozessen.

[0004] Es ist natürlich auch eine Fülle von Regelungsmethoden bekannt, um eine technische Funktion der Maschine oder Anlage zu regeln. Für alle gängigen Regelungsmethoden ist die technische Funktion, die geregelt werden soll, zu modellieren, um anhand des Modells oder für das Modell dann einen geeigneten Regler für die technische Funktion zu entwerfen. Das Modell kann ein physikalisches Modell sein, beispielsweise in Form von Differentialgleichungen oder Differenzengleichungen, oder ein mathematisches, das anhand von Messdaten an der Maschine oder Anlage trainiert wird. Bekannte mathematische Modelle sind beispielsweise ein neuronales Netzwerk oder ein lokales Modellnetzwerk (LMN). Das Modell kann linear oder nichtlinear sein. Ein Modell soll die technische Funktion über den gewünschten Betriebsbereich der Maschine oder Anlage abbilden. Dazu stellt das Modell in Abhängigkeit von gewissen Betriebsparametern, die den Betriebsbereich aufspannen, einen Zusammenhang zwischen zumindest einer Eingangsgröße, beispielsweise einer Stellgröße eines Aktuators, und zumindest einer Ausgangsgröße dar. Ublicherweise werden mit dem Modell mehrere verschiedene Eingangsgrößen auf mehrere Ausgangsgrößen abgebildet (Mehrgrößensystem). Zur Regelung sind einfache Regler wie ein Pl oder PID Regler oder Zustandsraumregler bekannt, die oftmals anhand eines physikalischen Modells

entworfen werden, aber auch komplexere Regler wie ein Modellprädiktiver Regler (MPC, Modell Predictive Controller) oder evolutionäre Regelalgorithmen, die oftmals anhand eines trainierten mathematischen Modells entworfen werden. Insbesondere die Verwendung eines MPC ist für viele Anwendungen vorteilhaft, weil sich damit auf einfache Weise Stellgrößenbeschränkungen berücksichtigen lassen und sich gute Regelungsergebnisse erzielen lassen. Ein MPC eignet sich insbesondere dann, wenn die Sollwerte über die Zeit bekannt sind, weil ein MPC um einen gegebenen Prädiktionshorizont in die Zukunft rechnet, um den nächsten Stelleingriff für die Aktuatoren zu ermitteln.

[0005] Eine modellprädiktive Regelung basiert bekanntermaßen auf der Idee die zukünftigen Eingangsgrößen u zu ermitteln, die eine Abweichung zwischen bekannten zukünftigen Ausgangsgrößen yıer (als Sollwerte der Regelung) und den Ausgangsgrößen y, die sich als Reaktion auf die Eingangsgrößen u ergeben (was aus einem Modell des technischen Prozesses berechnet werden kann), minimiert. Dabei wird um einen vorgegebenen Prädiktionshorizont in die Zukunft gerechnet. Das wird üblicherweise als Optimierungsproblem formuliert, bei dem eine Kostenfunktion J durch Variation der Eingangsgrößen u minimiert wird. Oftmals wird die quadratische Abweichung zwischen yre: und y beurteilt, sodass sich die Kostenfunktion J = (yrer - y)“-Qy ergibt, mit einem Gewichtungsfaktor Q,. Ublicherweise wird auch gefordert, dass sich die Eingangsgröße u von einem Zeitschritt zum nächsten nicht zu stark ändert, was negative Auswirkungen auf die Regelung haben kann (z.B. Stabilität). Daher wird in der Kostenfunktion auch oftmals ein zweiter Term berücksichtigt der die Änderung der Eingangsgröße Au (im zeitdiskreten Fall) oder den Gradienten der Eingangsgröße u (im zeitkontinuierlichen Fall) bewertet, beispielsweise wieder als quadratische Anderung. Daneben können natürlich noch beliebige andere oder weitere Kostenterme in der Kostenfunktion berücksichtigt werden. Damit ergibt sich die Kostenfunktion beispielsweise zu J = (yıer - y)*-Qy + Au?-Qı, mit einem Gewichtungsfaktor Q,. Über die Gewichtungsfaktoren Q,, Qu hat man dabei zwei Einstellbarameter. Durch Minimierung dieser Kostenfunktion J erhält man eine optimale Folge von Eingangsgrößen u über den Prädiktionshorizont. Ublicherweise wird nur die nächste Eingangsgröße für die Regelung herangezogen und alle anderen weiter in der Zukunft liegenden Eingangsgrößen verworfen (was unter receding horizon bekannt ist). Das wird in jedem Zeitschritt der Regelung wiederholt.

[0006] Beispielsweise zeigt die WO 2014/131661 A1 einen nichtlinearen Regler für einen nichtlinearen Prozess. Dazu wird der Prozess mit einem LMN modelliert und anschließend für jedes lokale Modell des LMN mit einem multikriteriellen evolutionären Algorithmus ein lokaler PID Regler ermittelt. Der Regler ergibt sich dann aus einer Linearkombination der einzelnen lokalen Regler.

[0007] Die Dissertation Novak J., „Nonlinear System Identification and Control Using Local Model Networks“, Tomas Bata Universität, Feb. 2007 beschreibt den Entwurf von lokalen MPC Reglern für ein LMN. Dabei wird für jedes lokale Modell des LMN ein MPC Regler erstellt. Dieser Ansatz, für jedes lokale lineare Modell einen lokalen linearen MPC zu erstellen, ist wohl der intuitivste Ansatz. Allerdings bedingt dieser Ansatz einen sehr hohen Rechenaufwand, da in jedem Zeitschritt der Regelung für jedes lokale MPC eine Optimierung gelöst werden muss, was insbesondere bei erwünschten großen Prädiktionshorizonten der MPC zu sehr hohen Rechenaufwänden führen kann. Zusätzlich werden in einer MPC oftmals Zustandsbeobachter benötigt, um benötigte nicht gemessene oder direkt nicht messbare Zustände des technischen Systems oder der technischen Anlage aus verfügbaren Messkanälen zu schätzen, was wiederum eine aufwändige Parametrisierung des Zustandsbeobachters erforderlich macht. Damit eignet sich dieser Ansatz mit lokalen MPC nur bedingt für eine praktische Anwendung oder eine Echtzeitregelung mit sehr kurzen Zeitschritten der Regelung. Insbesondere auf Anlagen oder Maschinen mit eingeschränkten Rechnerkapazitäten lässt sich dieser Ansatz daher nicht einsetzen. Beispielsweise lässt sich ein solcher Regler mit den herkömmlichen Steuereinheiten eines Fahrzeuges und mit den Anforderungen an die Regelung in einem Fahrzeug nicht implementieren, zumindest nicht mit hinreichender Regelungsgüte (was einen hinreichend langen Prädiktionshorizont erforderlich macht).

[0008] Es ist daher eine Aufgabe der Erfindung eine modellprädiktive Regelung eines nichtlinearen technischen Prozesses anzugeben, wobei die Regelung effizient und mit wenig Rechenauf-

wand lösbar sein soll.

[0009] Diese Ausgabe wird durch die Merkmale der unabhängigen Ansprüche gelöst. Durch die Verwendung von Zustandsraummodellen mit gleicher Systemordnung für die Teilmodelle ist es möglich, die Zustandsraummodelle durch eine reguläre Zustandsraumtransformation in ein transformiertes Zustandsraummodell zu transformieren, für das dann die modellprädiktive Regelung erstellt werden kann. Der Vorteil ist der, dass dabei alle transformierten Teilmodelle dieselbe virtuelle Eingangsgröße und denselben neuen Zustandsvektor aufweisen und es damit ausreicht, nur eine einzige modellprädiktive Regelung in der virtuellen Eingangsgröße und dem neuen Zustandsvektor zu implementieren. Es ist damit nicht mehr erforderlich so wie im Stand der Technik, für jedes einzelne Teilmodell jeweils eine modellprädiktive Regelung zu implementieren und die Ergebnisse der einzelnen modellprädiktiven Regelungen zu kombinieren, was mit hohem Rechenaufwand verbunden ist. Durch die erfindungsgemäße Vorgehensweise kann der Rechenaufwand somit deutlich reduziert werden.

[0010] Als zusätzliches Ergebnis der erfindungsgemäßen modellprädiktiven Regelung kann aus dem neuen Zustandsvektor zusätzlich ein Schätzwert der zumindest einen Ausgangsgröße ermittelt werden. Dieser Schätzwert kann auch als Trajektorie der Ausgangsgröße in der Optimierung verwendet werden, womit sich eine Vorsteuerung realisieren lässt.

[0011] In einer besonders vorteilhaften Implementierung wird die Kostenfunktion in eine Funktion des zeitlichen Verlaufs der virtuellen Eingangsgröße, des vorgegebenen zeitlichen Verlauf des Sollwertes der Ausgangsgröße und von Matrizen, die sich aus der Zustandstransformation, den Systemmatrizen der Teilmodelle und der Kombinationsvorschrift ergeben, umgewandelt, wobei die Matrizen in jedem Zeitschritt der Regelung neu berechnet werden, und dass die Kostenfunktion zur Ermittlung des optimalen Verlaufs der zumindest einen virtuellen Eingangsgröße durch einen in der Regelungseinheit implementierten Lösungsalgorithmus durch Minimieren oder Maximieren optimiert wird. Damit kann das Modell des technischen Prozesses direkt bei der modellprädiktiven Optimierung berücksichtigt werden. Als Startwert für die Optimierung kann ein aktueller Messwert der Ausgangsgröße oder eine in einem vorhergehenden Zeitschritt prädizierte Ausgangsgröße des technischen Prozesses verwendet werden, womit entweder ein geschlossene Regelung oder eine Vorsteuerung erzielt werden kann. Besonders vorteilhaft ist dabei die Möglichkeit der Vorsteuerung unter gleichzeitiger Berücksichtigung von Randbedingungen der Eingangsgröße und/oder Ausgangsgröße.

[0012] Weitere Vorteile der erfindungsgemäßen nichtlinearen Regelung ergeben sich aus der nachfolgenden Beschreibung.

[0013] Die gegenständliche Erfindung wird nachfolgend unter Bezugnahme auf die Figuren 1 bis 4 näher erläutert, die beispielhaft, schematisch und nicht einschränkend vorteilhafte Ausgestaltungen der Erfindung zeigen. Dabei zeigt

[0014] Fig.1 die Regelung eines technischen Prozesses mit einer Regelungseinheit, [0015] Fig.2 eine Aufteilung eines betrachteten Betriebsbereichs in Teilbereiche, [0016] Fig.3 ein mögliches Modell des technischen Prozesses und

[0017] Fig.4 eine erfindungsgemäße modellprädiktive Regelung.

[0018] Die gegenständliche Erfindung geht allgemein von einem nichtlinearen technischen Prozess 1 einer technischen Maschine oder Anlage aus, der mit einem Regler 2 geregelt werden soll, wie in Fig.1 dargestellt. Der technische Prozess 1 weist üblicherweise eine Anzahl von Ausgangsgrößen y und eine Anzahl m von Eingangsgrößen u auf, sodass man im Allgemeinen von einem Mehrgrößensystem spricht. Die Eingangsgrößen u bewirken die Ausgangsgrößen y und beeinflussen auch eine Anzahl von Betriebsparametern S des technischen Prozesses 1, die sich im Betrieb einstellen. Die Eingangsgrößen u werden beispielsweise als Stellgrößen an Aktuatoren 3 des technischen Prozesses 1 eingestellt, sodass sich die Anzahl der Ausgangsgrößen y, und die Betriebsparameter s, ergeben. Natürlich sind die Eingangsgrößen u aus der physikalischen Implementierung des technischen Prozesses 1 oder der Aktuatoren beschränkt, sodass die Ein-

gangsgrößen u nur bestimmte Werte annehmen können. Damit ergibt sich auch ein möglicher Betriebsbereich des technischen Prozesses 1, der von den Betriebsparametern s aufgespannt wird. Bei einem Verbrennungsmotor sind typische Betriebsparameter s beispielsweise die Drehzahl und das Drehmoment des Verbrennungsmotors, aber auch Größen wie die Einspritzmenge, die Abgastemperatur oder ein Abgasmassenstrom. Typische zu regelnde Ausgangsgrößen sind beispielsweise emissionsrelevante Größen im Abgas oder im Ansaugtrakt (z.B. NOx, EGR-Rate, Verbrennungsluftverhältnis oder Ansaugluftmassenstrom) aber auch momentenrelevante Größen wie Ladedruck, Druck vor Turbolader oder gesamter Luftmassenstrom in den Zylinder. Diese Größen können beispielsweise mit einer Drosselklappe, einer Abgasrückführung (EGR), einem Turbolader mit variabler Geometrie (VTG), mit einem elektrischen Turbolader (Turbolader mit Elektromotor zum aktiven Antreiben des Turboladers, wobei der Elektromotor auch generatorisch zur Energierückgewinnung betrieben werden kann), usw. als Aktuatoren 3, beeinflusst werden. In verschiedenen technischen Prozessen 1 kommen natürlich verschiedene Aktuatoren 3 zum Einsatz. Die Regelung ermittelt damit die Stellgrößen (Eingangsgrößen) beispielsweise für die Abgasrückführung und den Turbolader, um die Ausgangsgröße (z.B. NOx und Luftmassenstrom) zu regeln. Natürlich gibt es, insbesondere in einem Fahrzeug, noch eine Fülle anderer technischer Prozesse, die geregelt werden, und die nicht alle einzeln aufgeführt werden können. Beispiele sind die Regelung der Verbrennung eines Verbrennungsmotors, eines Hybridantriebsstranges, einer Abgasnachbehandlung, eines Getriebes, usw. Je nach technischem Prozess 1 und je nachdem was geregelt werden soll, werden in der Regelung unterschiedliche Eingangsgrößen u, Ausgangsgrößen y und Betriebsparameter s zu berücksichtigen sein. Der technische Prozess 1 wird von einem Regler 2 geregelt, der in jedem Zeitschritt der Regelung anhand der implementierten Regelung die Eingangsgrößen u berechnet, die an den Aktuatoren 3 als Stellgrößen eingestellt werden müssen, damit die Ausgangsgröße y einem vorgegebenen Sollwert yıert, bzw. einen ZUkünftigen Verlauf des Sollwerts y‚e: über den Prädiktionshorizont bei einer modellprädiktiven Regelung, folgt.

[0019] Die Erfindung geht weiter davon aus, dass für den nichtlinearen technischen Prozess 1 ein Modell M vorhanden ist, das die zumindest eine Eingangsgröße u in Abhängigkeit von möglichen Betriebsparametern s auf die zumindest eine Ausgangsgröße y abbildet, also y = f(u, s). Das Modell M besteht aus einer Mehrzahl i, mit i=> 1, von Teilmodellen Mi (Fig.3). Der betrachtete Betriebsbereich 4 (was nicht dem gesamten möglichen Betriebsbereich 4 entsprechen muss) des technischen Prozesses 1 wird in | Teilbereiche TBi unterteilt, wie in Fig.2 dargestellt, und für jeden der i Teilbereiche gibt es ein Teilmodell Mi. Wieviele Teilmodelle Mi; es gibt und wie der betrachtete Betriebsbereich 4 in die | Teilbereiche aufgeteilt wird, ist für die Erfindung nicht entscheidend. Es erübrigt sich zu erwähnen, dass für verschiedene technische Prozesse 1 verschiedene Teilmodelle Mi zum Einsatz kommen und dass auch die Anzahl i der Teilmodelle Mi; unterschiedlich sein kann. Im Beispiel der Fig.2 wird der Betriebsbereich 4 von zwei Betriebsparametern s1, s2, beispielsweise die Drehzahl und das Drehmoment eines Verbrennungsmotors, aufgespannt. Selbstverständlich sind aber auch mehr oder weniger Betriebsparameter s möglich.

Der Ausgang des Modells M ergibt sich aus den Modellausgängen y; der Teilmodelle Mi nach einer vorgegebenen Kombinationsvorschrift, beispielsweise als gewichtete Summe der Modellausgänge yı der Teilmodelle M;. Dabei ist für jedes Teilmodell Mi eine Gewichtungsfunktion ©; gegeben, die den Einfluss des jeweiligen Teilmodells M; auf die zumindest eine Ausgangsgröße y des Modells M über den Betriebsbereich 4 beschreibt. Die Gewichtungsfunktion ©; ist üblicherweise eine Funktion der Betriebsparameter s, kann aber auch von den Eingangsgrößen u und/oder Ausgangsgrößen y beeinflusst werden. Beispielsweise kann ein Teilmodell Mi; nur in einem sehr engen Bereich um einen bestimmten Betriebspunkt (definiert durch eine Kombination der Betriebsparameter s) Einfluss auf eine Ausgangsgröße $ des Modells M haben, oder kann in einem größeren oder ganzen Bereich des betrachteten Betriebsbereichs Einfluss ausüben. Ein Beispiel einer Gewichtungsfunktion ©; ist eine Sigmoidfunktion. Allgemein ergibt sich die zumindest eine Ausgangsgröße $ des Modells M in der Form 5=X ii, wobei ein Teilmodell M; L

auch vergangene Eingangsgrößen und/oder Ausgangsgrößen verarbeiten kann, wie in Fig.3 angedeutet. Bei mehreren Ausgangsgrößen sind die Ausgangsgröße $ und der Modellausgang $ı

Vektoren mit den einzelnen Ausgangsgrößen als Vektorelemente. Ein typisches Beispiel eines solchen Modells M ist ein lokales Modellnetzwerk LMN. Es ist aber durchaus denkbar, dass die einzelnen Modellausgänge $%; der Teilmodelle Mi auch anders zur Ausgangsgröße # des Modells M kombiniert werden könnten. Die Ausgangsgröße 9 des Modells M ist damit eine Schätzung der tatsächlichen Ausgangsgröße y des technischen Prozesses 1.

[0020] Entscheidend für die Erfindung ist die Modellstruktur der Teilmodelle Mi. Jedes Teilmodell Mi ist als lineares Zustandsraummodell modelliert. Das Zustandsraummodell kann beispielsweise aus einer physikalischen Beschreibung des technischen Prozesses 1, beispielsweise mittels Differentialgleichungen oder Differenzengleichungen n-ter Ordnung, vorzugsweise maximal 2. oder 3.Ordnung, abgeleitet werden. In diesem Fall haben die Zustände x im Zustandsraummodell eine physikalische Entsprechung. Das Zustandsraummodell muss aber nicht zwingend eine physikalische Entsprechung haben, sondern kann auch als Black Box Modell implementiert sein. Jedes Teilmodell Mi ist als Zustandsraummodell in bekannter Zustandsraumdarstellung gegeben, also Fick) = Cixi(k)

. Üblicherweise sind die Matrizen

allgemein in zeitdiskreter Schreibweise mit Zeitschritt k in der Form

Xi = ArXi + Bu Yı = CiXi

Ai, Bi, Ci in den verschiedenen Darstellungsformen aber nicht gleich.

Darin sind A;, Bi, Ci die Systemmatrizen, die den technischen Prozess 1 beschreiben, x; die Zustände, %i die Ausgangsgröße(n) und u bezeichnet die Eingangsgröße(n). Je nach Anzahl der Eingangsgrößen u und Ausgangsgrößen $; kann das Teilmodell Mi als Eingrößensystem (SISO; single input single output), in diesem Fall wären B und C Vektoren, oder Mehrgrößensystem (MIMO); multiple input multiple output) vorliegen. Die Ordnung n des Zustandsraummodells ergibt sich aus der Größe der Systemmatrizen A;, Bi, Ci, beispielsweise ist die Systemmatrix A; einen x n Matrix und die Systemmatrix Bi; eine n x Anzahl der Eingangsgrößen Matrix und Ci; eine Anzahl der Ausgangsgrößen x n Matrix. Bei einem SISO System reduzieren sich die Matrizen Bi, Ci auf Vektoren. Zur Bestimmung des Modellausgangs $; des Modells M des technischen Prozesses 1 werden die Modellausgänge %; der Teilmodelle Mi mit der jeweiligen gegebenen Gewichtungsfunktion @i, des Teilmodells Mi in Zustandsraumdarstellung gewichtet und aufsummiert (oder anderweitig kombiniert), also beispielsweise 5(k)= X o;Cixi(k). Ublicherweise gilt dabei

L

[0021] Ein derartiges Modell M, oder die Teilmodelle Mi, kann bekannt oder vorgegeben sein oder kann trainiert werden. Zum Training kann der technische Prozess 1 beispielsweise auf einem Prüfstand (z.B. ein Motorprüfstand für einen Verbrennungsmotor) gemäß gewisser Vorgaben für Betriebsparameter s und Eingangsgrößen u betrieben werden und dabei die Ausgangsgrößen y gemessen oder aus anderen erfassten Messgrößen berechnet werden. Der betrachtete Betriebsbereich 4 wird dann in i Teilbereiche unterteilt und aus den diesen i Teilbereichen zugeordneten Eingangsgrößen, Ausgangsgrößen und Betriebsparametern s können mit hinlänglich bekannten mathematischen Trainingsmethoden die Teilmodelle Mi trainiert werden, also die Systemmatrizen A;, Bi, Ci der Teilmodelle Mi bestimmt werden. Beispielsweise wird in jedem Teilbereich TBi des Betriebsbereichs 4 zumindest ein Betriebspunkt BPi vermessen (Eingangsgröße, Ausgangsgröße, Betriebsparameter) und das Teilmodell Mi; dafür parametrisiert, sodass der Modellausgang %i des zugeordneten Teilmodells Mi den in diesem Betriebspunkt BPi gemessenen Ausgang bestmöglich approximiert. Solche Trainingsalgorithmen sind hinlänglich bekannt, weshalb nicht näher darauf eingegangen werden muss. Die Gewichtungsfunktionen @;, können vorgegeben werden oder vorgegeben sein, um die Gültigkeit oder den Einfluss eines Teilmodells M; in anderen Teilbereichen TBi festzulegen. Die Gewichtungsfunktionen @;, können aber genauso im Zuge des Trainings des Modells M mitbestimmt oder optimiert werden. Für die Erfindung kann davon ausgegangen werden, dass ein derartiges Modell M vorliegt. Das Training eines solchen Modells M ist nicht Teil der Erfindung, sondern bekannt.

oder äquivalent in zeitkontinuierlicher Form als

[0022] Für die Erfindung wird die Modellstruktur der Teilmodelle Mi; in der Zustandsraumdarstellung mit dem Zustandsvektor x zuerst mit einer regulären Zustandstransformation in ein System

mit gewählten neuen Zustandsgrößen in einem neuen Zustandsvektor z überführt. Die reguläre Zustandstransformation erfolgt allgemein mittels einer Transformationsmatrix T in der Form z = Tx. Das Ziel der Zustandstransformation ist bekanntermaßen das zugrundeliegende dynamische System in Form des Zustandsraummodells in äquivalenter, aber anderer mathematischer Weise darzustellen. Die Dynamik bleibt damit mit dieser Transformation aber erhalten. Es gibt verschiedene Wege, um ein Zustandsraummodell eines linearen SISO oder MIMO Systems in die Regelungsnormalform überzuführen. Diese Transformation in die Regelungsnormalform ist aber grundsätzlich hinlänglich bekannt und wird am Beispiel eines Mehrgrößensystems (MIMO) im zeitkontinuierlichen Fall kurz erläutert.

[0023] Die Transformation in die Regelungsnormalform erfolgt für ein lineares Zustandsraummodell vorzugsweise mit einer Transformationsmatrix T, die sich aus der Steuerbarkeitsmatrix Qs des linearen Zustandsraummodells ergibt. Die Transformation ist z.B. aus J. Lunze, Regelungstechnik 2, Springer Verlag, 9.Auflage, 2016 bekannt. Aus dem Zustandsraummodell

m

j=1 ergibt sich die reduzierte Steuerbarkeitsmatrix 0, bekanntermaßen zu Q; =

P=n [b,,..., AP171b.,-,bm-..., APm71 bb] mit den sogenannten Steuerbarkeitsindizes p;, mit 1 mit der Systemordnung n. Die Steuerbarkeitsindizes p; geben die Anzahl der linear unabhängigen Spalten in der Steuerbarkeitsmatrix Qs = [b4,..., bm, Ab, ..., Abm... , AP!b,--, AP 1b] des Zustandsraummodells an und ergeben sich aus der Modellierung des technischen Prozesses 1 durch das Zustandsraummodell. Die Transformation folgt aus dem Zustandsvektor x im Zustandsraummodell aus z=T-x mit dem neuen Zustandsvektor z in Regelungsnormalform.

[0024] Die Transformationsmatrix T ergibt sich zu T = [t,",...,t4TA?17* tm ",..., tm" APM71T. Der

pirte Vektor ti stellt die #1 Zeil der inversen reduzierten Steuerbarkeitsmatrix 05? dar.

[0025] Das wird am Beispiel eines MIMO Systems mit Systemordnung n=3 und mit zwei Eingangsgrößen und zwei Ausgangsgrößen und mit den Steuerbarkeitsindizes p1 = 1, p2 = 2 verdeutlicht. Die reduzierte Steuerbarkeitsmatrix ergibt sich dann zu Ö; = [b, bz Ab„]7, woraus die inverse reduzierte Steuerbarkeitsmatrix Ö7? folgt. Die Transformationsmatrix T ergibt sich dann zuT = [t” tt," t,TA]”, mitt," als erste Zeile t,” als dritten Zeile der inversen reduzierten Steuerbarkeitsmatrix Ö7*.

[0026] Die Transformation in die Regelungsnormalform ergibt sich dann aus dem Zustandsvektor X im Zustandsraummodell aus z=T-x bzw. x=T"z mit dem neuen Zustandsvektor z = [Zu1 Zip, Zma * Zmpm]' in Regelungsnormalform. Daraus lässt sich die folgende Transformationsvorschrift ableiten, mit ] = 1 ....m (Anzahl der Eingangsgrößen) und der virtuellen Eingangsgröße v;.

Zi1 = Zj,2 Zi,pj-1 7 Zip; Zip; = t} APIT71Z + t7APIT1Bu = vw;

[0027] In Matrixdarstellung folgt daraus z = Äz + Bv mit den virtuellen Eingangsgrößen

tTAPı tTAPı71 =vV= T-1z+ : Bu. tTAPm et APm-1 H

[0028] Die virtuellen Eingangsgrößen v ergeben sich also aus der Zustandstransformation in die Regelungsnormalform aus den Systemmatrizen A, B des Zustandsraumodells. Es ergeben sich

genauso viele virtuelle Eingangsgrößen v als es physikalische Eingangsgrößen u gibt. Weiters ist erkennbar, dass die Abbildung der virtuellen Eingangsgrößen v auf die neuen Zustandsgrößen z in Regelungsnormalform durch eine Integratorenkette gegeben ist (in zeitdiskreter Darstellung eine Kette von Verzögerungsgliedern). Außerdem kann die physikalische Eingangsgröße u aus der virtuellen Eingangsgröße v rekonstruiert werden in der Form u=H"v-H'Kz. Die Ausgangsgröße y ergibt sich aus der Transformation zu $ = CT! z

W ZB = tTAT 717 + tTBu = Zi [0029] Für das obige Beispiel ergäbe sich Z21 = Z2,2 bzw. in Matrixdarstel2,2 = AT 1z+t7ABu= v, lung 000 1 0 } 4 A tT Z1,1 0 0 1|z+10 0/vmitv = ı] 17+| : Buundz= [zu Die ZU0 0 0 Lo__ 11 lea] 2] LAT A Z2,2 aa

standsgrößen z in Regelungsnormalform sind “damit im zeitdiskreten Fall vergangene virtuelle Eingangsgrößen v und im zeitkontinuierlichen Fall zeitliche Integrale der virtuellen Eingangsgrößen v.

[0030] Für ein SISO System ergibt sich die Transformation in die Regelungsnormalform in analoge Weise.

[0031] Diese Transformation ist für jedes der i Teilmodelle M; durchzuführen, wobei sich für die Teilmodelle Mi aus der Transformation unterschiedliche Matrizen T;, Ki, Hi, Wi ergeben. Die Transformation kann natürlich in äquivalenter Weise auch in zeitdiskreter Darstellung durchgeführt werden.

[0032] Auf dasselbe Ergebnis kommt man, wenn man beispielsweise den Ausgang y im Zustandsraummodell jedes Teilmodells Mi solange nach der Zeit ableitet, bis in der Ableitung das erste Mal der Eingang u vorkommt. Auch das führt zur Regelungsnormalform für die Zustandsraumdarstellung und zu Matrizen H, K, W, die sich nun aber nicht direkt aus Matrixprodukten einer Transformationsmatrix T und den Systemmatrizen A;, Bi, Ci ergeben. Auch hierbei ergibt sich ein transformiertes Zustandsraummodell in der Form z = Az + Bv.

[0033] Diese Zustandstransformation in die Regelungsnormalform wird für jedes der i Teilmodelle Mi durchgeführt. Gemäß der Erfindung hat jedes Teilmodell Mi dieselbe Systemordnung n. Dieser Ansatz ist für ein Modell M eines technischen Prozesses 1, das aus Teilmodellen Mi zusammengesetzt ist, sinnvoll, weil die Teilmodelle Mi den technischen Prozess 1 nur in unterschiedlichen Teilbereichen TBi des Betriebsbereiches 4 modellieren. Damit sind die Strukturen der in die Regelungsnormalform transformierten Zustandsraummodelle der Teilmodelle Mi gleich und es treten lediglich unterschiedliche Matrizen Ti, Ki, Hi, Wi auf. Die virtuellen Eingangsgrößen v und die transformierten Zustände z gelten aber für alle Teilmodelle Mi, weil man bei der Transformation jedes Teilmodells Mi auf dieselbe Integratorenkette (Kette von Zeitverzögerungsgliedern) kommt. Dieser Umstand kann nun für eine erfindungsgemäße modellprädiktive Regelung des nichtlinearen technischen Prozesses 1, der wie beschrieben mit linearen Teilmodellen Mi modelliert wird, ausgenutzt werden.

[0034] Jedes Teilmodell Mi wird wie beschrieben in das transformierte Zustandsraummodell, beispielsweise in die Regelungsnormalform, transformiert, wobei die virtuellen Eingangsgrößen v und die transformierten Zustände z für jedes Teilmodell Mi gültig sind. Es kann somit nur eine einzige modellprädiktive Regelung entworfen werden, die die virtuelle Eingangsgrößen v regelt und das transformierte Zustandsraummodell z = Äz + Bv, also die Integratorenkette (Kette von Verzögerungsgliedern) bei einer Transformation in die Regelungsnormalform, als das der modellprädiktiven Regelung zugrunde liegende Modell verwendet, womit sich auch die transformierten Zustände z ergeben.

71717

[0035] Für eine vorteilhafte Implementierung der erfindungsgemäßen modellprädiktiven Regelung kann man sich den Umstand der Modellierung durch das Modell M des technischen Prozesses 1 mit Teilmodellen Mi in Form von Zustandsraummodellen zu Nutze machen, was nachfolgend in zeitdiskreter Formulierung erläutert wird. Bei dieser Vorgehensweise wird die Rücktransformation der virtuellen Eingangsgrößen in die physikalischen Eingangsgrößen und die Summierung über die Teilmodelle Mi; in der Optimierung mitberücksichtigt und somit in einem Rechenschritt berechnet.

[0036] Der nichtlineare technische Prozess 1 ist wie beschrieben mit einer Anzahl i Teilmodellen Mi in der Form von Zustandsraummodellen (nun in zeitdiskreter Darstellung)

Ik) = DIA) i modelliert. Aus der Zustandstransformation in die Regelungsnormalform der Teilmodelle Mi erge-

ben sich wie oben ausgeführt die Eingangsgröße u(k) und die Ausgangsgröße y(k) des technischen Prozesses 1 in der Form

u(k) = ) YHF* uk) DH, 2C0)

Nx Me

Ik) = YkIGT* 20)

NL Wr

[0037] Für die modellprädiktive Regelung werden die zukünftigen Ausgangsgrößen y(k+1), ..., y(k+Ng-1) und die zukünftigen Eingangsgrößen u(k+1), ..., u(k+N,-1) für den Prädiktionshorizont N; benötigt. Ferner wird wieder eine Kostenfunktion J für die Optimierung benötigt. Mit der Kostenfunktion J wird die Abweichung der Ausgangsgröße y zu gegebenen Sollwerten der Ausgangsgröße yıer ZUm jeweiligen Zeitpunkt bewertet. Das erfolgt beispielsweise wieder über die quadratische Abweichung und mit einem Einstellbarameter Q, der Regelung, wobei Qy für ein MIMO System eine Matrix mit mehreren Einstellbarametern ist. Das führt in der allgemeinen Formulierung (für SISO und MIMO) zu J = (Yıer- Y)' Qy(Yıer- Y). Vorzugsweise wird auch die Differenz zwischen zwei Stelleingriffen, also die Anderung der Eingangsgrößen u, bewertet, beispielsweise wieder als quadratische Abweichung. Das kann beispielsweise mit einem zusätzlichen Term (AU" Qu AU) in der Kostenfunktion J berücksichtigt werden. Q ist darin wieder ein Einstellbarameter der Regelung, wobei Qı für ein MIMO System eine Matrix mit mehreren Einstellbarametern ist. Es sei aber angemerkt, dass die Kostenfunktion J die Abweichung der Ausgangsgrößen y von Sollwerten y‚e: natürlich auch anders bewerten kann, was in gleicher Weise auch für die Änderung (oder den Gradienten) der Eingangsgrößen u zutrifft. Ebenso können in der Kostenfunktion J auch weitere Terme berücksichtigt werden. Beispielsweise könnte in einem weiteren Term die virtuelle Eingangsgröße V bewertet werden, z.B. wie die Eingangsgröße U. Es könnten ferner auch andere mit dem technischen Prozess 1 in Zusammenhang stehenden Terme in der Kostenfunktion J bewertet werden. Beispielsweise könnte für einen Verbrennungsmotor als technischen Prozess 1 auch der Kraftstoffverbrauch oder eine Emissionsgröße als weiterer Term in der Kostenfunktion J berücksichtigt werden, gegebenenfalls wieder mit einem Einstellbarameter. Damit könnten in der Kostenfunktion ein hoher Kraftstoffverbrauch oder zu hohe Emissionen bestraft werden.

[0038] Die vorteilhafte Ausgestaltung der Erfindung wird im Folgenden am Beispiel der Kostenfunktion J = (Yıer- Y)' Qy(Yıer - Y) + AU" Qu AU beschrieben.

[0039] Für die modellprädiktive Regelung werden die zukünftigen Ausgangsgrößen $ (k+1), ..., J (k+N;-1) über den Prädiktionshorizont N, benötigt. Diese können für den Prädiktionshorizont Np angeschrieben werden.

Pk) = WrZCK)

FI +1) = We+1Äz0k) + y(K) + Wer Bock) — WezCk)

Y(k +2) = Wy+42Ä?zCk) + y(k) + Wr42AÄBvCk) + Wr+2BoCk + 1) — WrzCk)

J(k +3) = Wr43Ä®zCk) + Wr+3Ä?BuCk) + Wr43ÄBvCk + 1) + Wr+3BuCk +2) — WrzCk)

[0040] In Matrixform ergibt sich das mit der Einheitsmatrix | und der oben definierten Matrix Wx ZU

Yk+1) Wer d-We I We+1B 0 0 v(K)

Fk +2) _ | Wer22? Wr I ZO] + | Wx+24B Wer2B 0 [| UCk + DI]

FIKA3] |[We43 A Wr 1|0OL WE, 3A2B WeyzAB WegzB |] 0 +2) Dy

[0041] In gleicher Weise kann man die Abweichungen Au(d(k) = u(k) - u(k-1) bis Au(k+Np-1) = u(k+Ng-1) - u(k+N,-2), usw. anschreiben und man erhält in Matrixform mit den oben definierten Matrizen Mx, N

Au(k) My 1! Au(k +1) | _ |[Myx41Ä— Mr, 0 [ z(k) Au(k + 2) My+2A? — My+2A 0 ulk 1) EL : : Nx 0 0 0717 v(k) My+1B — Nx Ne+1 0 0llv(k+1)

My424AB — My41B Mp4r2B — Neys Ne4y2 0||9Ck +2)

Dy V

Wy+1Ä-Wy 10 A2

[0042] Definiert man die Matrizen FR, = Wr+24 Wr I Olund Fı= Wy+34° -Wr I 0

My 0 I ; ©

My+1Ä — Mg 00

My+2Ä?* — My42AÄ 0 0

Y E.@)]| 709 Dy(@) erhält man zusammengefasst | AU = N (©) y@k) |+ lo Cl V. Die Matrizen Fy, Fu, ®y, O1 uNP)! uk — 1) u

sind von den Gewichtungsfunktionen ©; der Teilmodell M; abhängig. Die Matrizen 4,3 ergeben sich aus der obigen Beschreibung der Integratorketten oder Ketten von Verzögerungsgliedern als Zustandsraumsystem in Regelungsnormalform. y(k) bezeichnet die aktuell gemessene Ausgangsgröße und u(k-1) im vorherigen Zeitschritt ermittelte Eingangsgröße und stellen die Startwerte der Prädiktion in der Optimierung dar.

[0043] Setzt man das in die definierte Kostenfunktion J = (Yıer- Y)' Qy (Yrer - Y) + AU" Qu AU ein erhält man

J =vr05 o0[$ Ole 16 6_oultex)

205 01[% 011-8 eV

f

[0044] Weitere sich ergebende, von V unabhängige Terme können in der Kostenfunktion J, die nach V optimiert wird, unberücksichtigt bleiben, weil diese in der Kostenfunktion J nur konstante Terme wären. Die Kostenfunktion J für die Optimierung ergibt sich dann zu J = V'EV- 2V'f'. Diese Kostenfunktion J ist der Input in die Optimierung, insbesondere Input in einen Lösungsalgorithmus für das Optimierungsproblem

Vopt = min VTEV —-2VT FT.

[0045] Die Lösung dieses Optimierungsproblems ist ein Quadratisches Programm, das in der Literatur gut untersucht ist und für das es geeignete Lösungsalgorithmen gibt.

[0046] In der Optimierung wird die Kostenfunktion J üblicherweise durch Variation der virtuellen Eingangsgrößen v minimiert wird, wobei auch eine Maximierung denkbar ist. Das kann iterativ erfolgen, mit einem vorgegebenen Abbruchkriterium, beispielsweise eine Anzahl von Iterationsschritten oder ein gegebenes Abbruchkriterium. Es kann für bestimmte Kostenfunktionen J aber auch eine analytische Lösung der Optimierung gefunden werden.

[0047] Mit der derartigen modellprädiktiven Regelung wird der virtuelle Eingangsgrößenvektor v (bei SISO System nur eine Eingangsgröße) derart ermittelt, dass die gewichtete Summe der Modellausgänge %;, also der Ausgangsgrößenvektor $ (bei einem SISO System nur eine Ausgangsgröße) des Modells M, dem vorgegebenen Referenzwert yıe: gemäß dem Optimierungskriterium, das in der Kostenfunktion J festgelegt ist, folgt. Dabei kann bei entsprechender Wahl der Kostenfunktion J gleichzeitig auch die gewichtete Summe der Eingangsgrößen berücksichtigt werden. Da die Bestimmung der virtuellen Eingangsgrößen v auf der gewichteten Summe von %;,, und gegebenenfalls auch u;, basiert, wird immer dieselbe Kostenfunktion betrachtet, womit in der Optimierung trotz der Anzahl i von Teilmodellen Mi nur eine Kostenfunktion J optimiert werden muss. Das ermöglicht eine besonders effiziente Realisierung einer modellprädiktiven Regelung.

[0048] Der Vorteil dieses Ansatzes liegt auch darin, dass unmittelbar Randbedingungen und Nebenbedingungen auf die Eingangsgrößen u und die Ausgangsgrößen y bedarfsweise berücksichtigt werden können, beispielsweise Au" < Au <= Au"* oder y7" << y"*, wobei die mit min, max bezeichneten Grenzwerte in der Optimierung vorgegeben werden, beispielsweise durch physikalische Vorgaben des technischen Prozesses 1 oder der Aktuatoren 3.

[0049] Die Einstellbarameter Q,, Qu müssen in der Optimierung nicht konstant sein. Diese können sich auch von Zeitschritt zu Zeitschritt ändern, können aber auch abhängig sein vom jeweiligen Betriebspunkt BP und könnten sogar über den Prädiktionshorizont Nr veränderlich sein.

[0050] Diese Optimierung ist in jedem Zeitschritt k der Regelung zu lösen, wobei in jedem Zeitschritt auch die Matrizen E, f, bzw. die darin enthaltenen Matrizen, zu bestimmen sind, da diese von der Kombinationsvorschrift der Teilmodelle Mi, beispielsweise den Gewichtungsfunktionen @;, abhängig sind. Solche Matrizenoperationen lassen sich aber in jedem Zeitschritt k lösen.

[0051] Die Optimierung liefert somit den optimalen zeitlichen Verlauf (zeitdiskret oder zeitkontiniuierlich) der virtuellen Eingangsgrößen V über den Prädiktionshorizont Nr. Wie in der modellprädiktiven Regelung üblich werden vom ermittelten Verlauf zumindest die für den aktuellen Zeitschritt k gültigen Größen verwendet, und die anderen verworfen. Für ein SISO System enthält V nur den Verlauf einer virtuellen Eingangsgröße über den Prädiktionshorizont Np, für ein MIMO System den Verlauf einer virtuellen Eingangsgröße v über den Prädiktionshorizont Np für jede physikalische Eingangsgröße u. Damit kann die physikalische Eingangsgröße u(k) wie oben aus-

geführt rekonstruiert werden. Ebenso kann eine Schätzung der Ausgangsgröße $(k) ermittelt werden, wie oben ausgeführt. Die derart ermittelte Ausgangsgröße $(k) ist dabei eine Schätzung der sich ergebenden tatsächlichen Ausgangsgrößen y, wenn an das technische System 1 die Eingangsgrößen u angelegt werden, beispielsweise als Stellgrößen an den Aktuatoren 3.

[0052] Eine erfindungsgemäß ausgeführte Regelungseinheit 2 ist in Fig.4 dargestellt. Die Regelungseinheit 2 ist beispielsweise als Recheneinheit (Hardware) ausgeführt und die modellprädiktive Regelung MPC ist in Software auf der Recheneinheit implementiert. In der modellprädiktiven Regelung MPC wird in jedem Zeitschritt k der Regelung die oben beschriebene Optimierung der gewählten Kostenfunktion J gelöst. Daraus erhält man die im aktuellen Zeitschritt k gültigen virtuellen Eingangsgrößen v(k) und die virtuellen Zustandsgrößen z(k), aus denen die physikalischen Eingangsgrößen u(k) des technischen Prozesses 1 rekonstruiert werden, beispielsweise aus u(k) = N«v(k) + Mxz(k). Zusätzlich kann auch ein Schätzwert $(k) der Ausgangsgrößen y ermittelt werden, beispielsweise aus j(k) = Wız(k).

[0053] Zur Berechnung wird die aktuelle Ausgangsgröße y(k) am technischen Prozess 1 gemessen und der Regelungseinheit 2 rückgeführt, womit ein geschlossener Regekkreis realisiert werden kann.

[0054] Anstelle der Verwendung einer gemessenen Ausgangsgröße y(k) als Startwert könnte für die modellprädiktive Regelung MPC auch der ermittelte Schätzwert der Ausgangsgröße $(k) verwendet werden. Dazu kann der ermittelte Schätzwert der Ausgangsgröße $(k) der Regelungseinheit2 rückgeführt werden, wie in Fig.4 strichliert angedeutet. Dabei könnte man immer die im vorherigen Zeitschritt (k-1) ermittelte, geschätzte Ausgangsgröße $(k+1) als Ausgangsgröße $(k) im aktuellen Zeitschritt (k) verwenden, weil im Zeitschritt (k) die aktuelle Ausgangsgröße $(k) erst berechnet wird und somit nicht vorliegt. Damit erhält man eine Art Vorsteuerung des technischen Prozesses 1, die zusätzlich Begrenzungen der Eingangsgröße u (Stellgrößenbegrenzungen) einhält.

[0055] Die gemessene Ausgangsgröße y(k) könnte aber auch verwendet werden, um die Teilmodelle Mi zu verbessern, beispielsweise indem das Training der Teilmodelle Mi in gewissen Abständen mit aktuellen Werten der Ausgangsgröße y(k) und zugehörigen Eingangsgrößen u wiederholt wird, was auch online erfolgen könnte.

[0056] Nachfolgend werden als Beispiele von technischen Prozessen 1, die mit einer erfindungsgemäßen modellprädiktiven Regelung geregelt werden könnten, die Luftpfadregelung eines Verbrennungsmotors, die Kühlkreisregelung eines Personen- oder Lastkraftwagens, die Regelung einer Brennstoffzelle, die Regelung der Wärmerückgewinnung aus dem Abgas, die Regelung eines Hybridantriebstranges und eine Adaptive Fahrzeuggeschwindigkeitsregelung kurz erläutert.

[0057] Bei der Luftpfadregelung eines Verbrennungsmotors werden typischerweise die Ausgangsgrößen Abgasdruck oder Ladedruck und NOx geregelt. Als Stellgrößen der Regelung (Eingangsgrößen) kommen die EGR Position, die VGT Position, ein Sollwert eines elektrischen Turboladers, eine Abgas- oder Ansaugdrosselklappenstellung oder der Winkel Einspritzbeginn in Frage. Beschränkungen des technischen Prozesses 1 sind beispielsweise eine Turboladerdrahzahl, eine Abgastemperatur oder ein Verbrennungsluftverhältnis. Als zusätzlicher Kostenterm in der Kostenfunktion J kann auch der spezifische Kraftstoffverbrauch bewertet werden.

[0058] Bei einer Kühlkreisregelung eines Personen- oder Lastkraftwagens werden beispielsweise optimale Stellgrößen (Eingangsgrößen) für eine Kühlwasserpumpe, einen Lüfter, einer Durchflussmenge-Ventilöffnung oder eine elektrisch verstellbare Motorgrill-Jalousie-Position (Variable Grill Shutter Position) geregelt, sodass die Ausgangsgrößen wie Motortemperatur und Ladelufttemperatur geregelt werden. Dabei kann in der Kostenfunktion ebenfalls wieder ein zusätzlicher Kostenterm, der den Kraftstoffverbrauch bewertet, enthalten sein.

[0059] Ein anderes Beispiel eines technischen Prozesses 1 ist eine Brennstoffzelle. Dabei können die Ausgangsgrößen Luftmassenstrom der Kathodenseite und Kathodengegendruck mittels der Eingangsgrößen (Stellgrößen) Kompressorleistung oder Gegendruckventilstellung und Ven-

tilstellung der Wasserstoffzufuhr geregelt werden.

[0060] Eine andere Funktion eines Fahrzeugs, die als technischer Prozess 1 geregelt werden könnte ist die Wärmerückgewinnung aus dem Abgas (Waste Heat Recovery, WHR). Hierbei könnten als Ausgangsgrößen der Regelung Kühlwasser Temperaturanforderung, ein Wärmefluss in das WHR System oder eine Pumpenleistung im Wärmerückgewinnungskreis verwendet werden. Als Stellgrößen (Eingangsgrößen) kämen ein Massenfluss, der Druck oder die Temperatur (vor allem vor einem Verdampfer) eines Arbeitsmittels im Wärmerückgewinnungskreis in Frage. Als Beschränkungen könnten in der Optimierung Druck und Temperaturbeschränkungen des Arbeitsmittels, auch an verschiedenen Positionen, vor allem vor dem Verdampfer), ein maximal zulässige Wärmeeintrag in den Kühlkreis oder ein maximale zulässige Kühlwassertemperatur verwendet werden. Als zusätzliche Kostenterme in der Kostenfunktion könnte eine WHR System Effizienz oder eine WHR Turbinenausgangsleistung berücksichtigt werden.

[0061] Eine modellprädiktive Regelung ist im automotiven Umfeld besonders vorteilhaft, weil in modernen Fahrzeugen viele Funktionen anhand einer vorausschauenden Beurteilung von Fahrsituationen geregelt werden, insbesondere anhand einer vorausschauenden Beurteilung einer Fahrstrecke. Daraus können die zukünftigen zeitlichen Verläufe der Sollwerte Yref ermittelt werden, die in der modellprädiktiven Regelung benötigt werden. Das gilt insbesondere auch für die obigen Beispiele.

[0062] Beispielsweise kann das auch besonders vorteilhaft bei der Regelung eines Hybridantriebstranges eines Fahrzeuges verwendet werden, um die Momenten-Aufteilung zwischen Elektromotor und Verbrennungsmotor (Torque Split) zu optimieren. Die modellprädiktive Regelung basiert beispielsweise auf Information der davorliegenden Stecke, wie z.B. Steigung und Kurvenradius, und der zu erwartenden Fahrzeuggeschwindigkeit.

[0063] Auch die Prädiktive Adaptive Fahrzeuggeschwindigkeitsregelung (Predictive Adaptive Cruise Control) ist ein Beispiel für eine derartige Anwendung. Beispielsweise kann basierend auf einer Sollwertvorgabe des Fahrers (Cruise Control Set Speed) und einem erlaubtem Toleranzband (Stellgrößen-Verstellbereich) unter Berücksichtigung prädiktiver Information über die davorliegende Strecke eine optimale Geschwindigkeits-Trajektorie berechnet werden.

Claims (1)

1. Patentansprüche

   1. Verfahren zur Regelung eines nichtlinearen technischen Prozesses (1) mit zumindest einer Eingangsgröße (u) und zumindest einer Ausgangsgröße (y), mit einer modellprädiktiven Regelung (MPC), wobei mit der modellprädiktiven Regelung (MPC) in jedem Zeitschritt der Regelung die zumindest eine Eingangsgröße (u) des nichtlinearen technischen Prozesses (1) berechnet wird, die dem technischen Prozess (1) im aktuellen Zeitschritt der Regelung zur Regelung aufgeschaltet wird, und der nichtlineare technische Prozess (1) mit einem Modell (M) modelliert wird, das die zumindest eine Eingangsgröße (u) auf eine Schätzung der zumindest einen Ausgangsgröße (%) abbildet, wobei das Modell (M) aus einer Anzahl von Teilmodellen (Mi) besteht und die Ausgangsgröße ($) des Modells (M) aus einer gegebenen Kombinationsvorschrift der Modellausgänge ($;) der Teilmodelle (M;) gebildet wird, dadurch gekennzeichnet, dass für jedes Teilmodell (Mi) ein lineares Zustandsraummodell der Form xi(k + 1) = Ayxi(k) + Biu(k) oder Xi = AyXi + B; Yilk) = Cixilk) Yı = CiXi tor des Teilmodells (Mj), $i der Modellausgang des Teilmodells (M;), u die Eingangsgröße des technischen Prozesses und A,, Bi, Ci sich aus der Modellierung des technischen Prozesses (1) ergebende Systemmatrizen sind, wobei jedes Zustandsraummodell dieselbe Systemordnung n hat, dass die Zustandsraummodelle der Teilmodelle (Mi) jeweils mit einer regulären Zustandstransformation in ein transformiertes Zustandsraummodell in Regelungsnormalform der Form z = Az + Bv transformiert werden, mit Matrizen A, B, die sich aus der Transformationsvorschrift ergeben, mit gewählten neuen Zustandsgrößen in einem neuen Zustandsvektor z und mit zumindest einer sich aus der Zustandstransformation ergebenden virtuellen Eingangsgröße v, wobei jedes transformierte Teilmodell (Mi) dieselbe zumindest eine virtuelle Eingangsgröße v und denselben virtuellen Zustandsvektor z aufweist, dass eine auf einer Regelungseinheit (2) implementierte modellprädiktive Regelung (MPC) für die einzige virtuelle Eingangsgröße v verwendet wird, um den zeitlichen Verlauf der zumindest einen, für alle transformierten Teilmodelle (M;) gültigen, virtuellen Eingangsgröße v über einen vorgegebenen Prädiktionshorizont Np zu berechnen, dass aus dem zeitlichen Verlauf der derart berechneten zumindest einen virtuellen Eingangsgröße v die virtuelle Eingangsgröße v zumindest des aktuellen Zeitschrittes der Regelung verwendet wird, um daraus in der Regelungseinheit (2) die zumindest eine Eingangsgröße (u) des technischen Prozesses (1) des aktuellen Zeitschrittes zu berechnen und dass die berechnete Eingangsgröße (u) als Stellgröße für die Regelung des technischen Prozesses (1) verwendet wird.

   verwendet wird, wobei xi ein Zustandsvek-

   2, Verfahren nach Anspruch 1, dadurch gekennzeichnet, dass für die modellprädiktive Regelung (MPC) eine Kostenfunktion J verwendet wird, die in der Regelungseinheit (2) in jedem Zeitschritt der Regelung über den Prädiktionshorizont Ne durch Minimieren oder Maximieren optimiert wird.

   3. Verfahren nach Anspruch 2, dadurch gekennzeichnet, dass die die Abweichung eines prädizierten zeitlichen Verlaufs der Ausgangsgröße Y über den Prädiktionshorizont Np von einem vorgegebenen zeitlichen Verlauf eines Sollwertes der Ausgangsgröße Y,e: bewertet.

   4. Verfahren nach Anspruch 2, dadurch gekennzeichnet, dass in der Kostenfunktion J zusätzlich ein Kostenterm verwendet wird, der über den Prädiktionshorizont Ne die Anderung der zumindest einen Eingangsgröße u des technischen Prozesses (1) von einem Zeitschritt zum nächsten Zeitschritt bewertet.

   5. Verfahren nach einem der Ansprüche 2 bis 4, dadurch gekennzeichnet, dass die Kostenfunktion J in eine Funktion des zeitlichen Verlaufs V der virtuellen Eingangsgröße v, des vorgegebenen zeitlichen Verlauf Yıe: des Sollwertes der Ausgangsgröße y und von Matrizen, die sich aus der Zustandstransformation, den Systemmatrizen A;, Bi, Ci der Teilmodelle (M)) und der Kombinationsvorschrift ergeben, umgewandelt wird, wobei die Matrizen in jedem Zeitschritt der Regelung neu berechnet werden, und dass die Kostenfunktion J zur Ermittlung

   10.

   11.

   12.

   Österreichischer AT 522 290 B1 2021-01-15

   des optimalen Verlaufs der zumindest einen virtuellen Eingangsgröße v durch einen in der Regelungseinheit (2) implementierten Lösungsalgorithmus durch Minimieren oder Maximieren optimiert wird.

   Verfahren nach einem der Ansprüche 2 bis 5, dadurch gekennzeichnet, dass bei der Optimierung der Kostenfunktion J vorgegebene Randbedingungen für die zumindest eine Eingangsgröße u und/oder für die zumindest eine Ausgangsgröße y berücksichtigt werden.

   Verfahren nach einem der Ansprüche 2 bis 6, dadurch gekennzeichnet, dass als Startwert für die Optimierung ein aktueller Messwert y(k) der Ausgangsgröße y oder eine in einem vorhergehenden Zeitschritt prädizierte Ausgangsgröße ($(k)) des technischen Prozesses verwendet wird.

   Verfahren nach Anspruch 1, dadurch gekennzeichnet, dass die Ausgangsgröße (9) des Modells aus einer mit einer Gewichtungsfunktion @;, gewichteten Summe der Modellausgänge (%) der Teilmodelle (Mi) zusammengesetzt wird.

   Verfahren nach Anspruch 1, dadurch gekennzeichnet, dass ein betrachteter Betriebsbereich (4) des technischen Prozesses (1) in Teilbereiche (TBi) unterteilt wird und für jeden Teilbereich (TBi) ein Teilmodell (Mi) verwendet wird.

   Verfahren nach Anspruch 9, dadurch gekennzeichnet, dass in einem Teilbereich (TBi) zumindest eine Kombination der zumindest einen Ausgangsgröße y in Abhängigkeit von der zumindest einen Eingangsgröße u verwendet wird, um das für diesen Teilbereich (TBi) gültige Teilmodell (M;) zu trainieren.

   Verfahren nach Anspruch 1, dadurch gekennzeichnet, dass in der modellprädiktiven Regelung (MPC) aus dem neuen Zustandsvektor z zusätzlich ein Schätzwert (9) der zumindest einen Ausgangsgröße y ermittelt wird.

   Regelungseinheit zur Regelung eines nichtlinearen technischen Prozesses (1) mit zumindest einer Eingangsgröße (u) und zumindest einer Ausgangsgröße (y), mit einer modellprädiktiven Regelung (MPC), wobei die modellprädiktive Regelung (MPC) in der Regelungseinheit (2) implementiert ist und in jedem Zeitschritt der Regelung einen Wert der zumindest einen Eingangsgröße (u) des nichtlinearen technischen Prozesses (1) berechnet, und der nichtlineare technische Prozess (1) in der Regelungseinheit (2) mit einem Modell (M) modelliert ist, das die zumindest eine Eingangsgröße (u) auf eine Schätzung der zumindest einen Ausgangsgröße (y) abbildet, wobei das Modell (M) aus einer Anzahl von Teilmodellen (Mi) besteht und die Ausgangsgröße (%) des Modells (M) aus einer gegebenen Kombinationsvorschrift der Modellausgänge ($;) der Teilmodelle (Mi) gebildet wird, dadurch gekennzeichnet, dass jedes Teilmodell (Mi) als Ilineares Zustandsraummodell der Form x;(k +1) = Ayxi(k) + Bru(k) oder Xi = Aixi + Biu Jilk) = Cixi(k) Yı= CiXi

   vektor des Teilmodells (Mi), $i der Modellausgang des Teilmodells (Mi), u die Eingangsgröße des technischen Prozesses und A,, Bi, Ci sich aus der Modellierung des technischen Prozesses (1) ergebende Systemmatrizen sind, wobei jedes Zustandsraummodell dieselbe Systemordnung n hat, dass eine Zustandstransformation in Regelungsnormalform implementiert ist, die die Zustandsraummodelle der Teilmodelle (Mi) jeweils in ein transformiertes Zustandsraummodell in Regelungsnormalform der Form z = Az + Bv transformiert, mit Matrizen A,B, die sich aus der Transformationsvorschrift ergeben, mit gewählten neuen Zustandsgrößen in einem neuen Zustandsvektor z und mit zumindest einer sich aus der Zustandstransformation ergebenden virtuellen Eingangsgröße v, wobei jedes transformierte Teilmodell (Mi) dieselbe zumindest eine virtuelle Eingangsgröße v und denselben virtuellen Zustandsvektor z aufweist, dass die Regelungseinheit (2) mit der implementierten modellprädiktiven Regelung (MPC) einen zeitlichen Verlauf der zumindest einen, für alle transformierten Teilmodelle (M;) gültigen, virtuellen Eingangsgröße v über einen vorgegebenen Prädiktionshorizont Np berechnet, dass die Regelungseinheit (2) aus dem zeitlichen Verlauf der derart berechneten zumindest einen virtuellen Eingangsgröße v die virtuelle Eingangsgröße

   implementiert ist, wobei xi ein Zustands-

   v zumindest des aktuellen Zeitschrittes der Regelung verwendet, um daraus die zumindest eine Eingangsgröße (u) des technischen Prozesses (1) des aktuellen Zeitschrittes zu berechnen, und dass die Regelungseinheit (2) die berechnete Eingangsgröße (u) des technischen Prozesses (1) des aktuellen Zeitschrittes zur Regelung des technischen Prozesses (1) ausgibt.

   Hierzu 2 Blatt Zeichnungen