Übertragungsfunktionen werden durch Linearisierung der transienten Energiebilanz des Rührkesselheizers um einen stationären Betriebspunkt abgeleitet. Für einen ideal durchmischten Tank mit konstanter Flüssigkeitsmasse, Einlasstemperatur (T_i), Dampftemperatur (T_{st}), Durchflussrate (F), Wärmeübertragungsfläche (A) und Wärmedurchgangskoeffizient (U) ist das resultierende Modell von erster Ordnung. Es erzeugt separate Übertragungsfunktionen von der Einlassflüssigkeitstemperatur und der Dampftemperatur zur Tanktemperatur.

Die Tanktemperatur reagiert auf beide Störungen durch dieselbe Wärmespeicher- und Wärmeabfuhrdynamik. Der Pfad der Einlasstemperatur wird durch den Flüssigkeitsaustausch bestimmt, während der Pfad der Dampftemperatur durch die Wärmeübertragungskapazität der Heizschlange (UA) gesteuert wird.

Beginnen Sie mit der instationären Energiebilanz

Definition der Tank- und Durchflussannahmen

Nehmen Sie an, dass der Tank ideal gerührt ist, sodass die Auslasstemperatur der Tanktemperatur (T) entspricht. Gehen Sie zudem von konstantem Flüssigkeitsvolumen, konstanter Dichte, Wärmekapazität, Durchflussrate und Wärmeübertragungseigenschaften aus.

Der Tank enthält eine Flüssigkeitsmasse (M), und das Prozessfluid hat eine Wärmekapazität (C_p). Wärme tritt durch die Dampfschlange gemäß folgender Gleichung ein:

[ Q = UA(T_{st}-T) ]

wobei (U) der Wärmedurchgangskoeffizient und (A) die effektive Wärmeübertragungsfläche ist.

Schreiben Sie den Akkumulationsterm

Die im Tank gespeicherte thermische Energie ändert sich mit der Rate

[ MC_p\frac{dT}{dt} ]

Dieser Term repräsentiert das dynamische Verhalten: Eine größere Flüssigkeitsmasse oder Wärmekapazität erfordert mehr Energie, um die Tanktemperatur zu ändern.

Schreiben Sie die Durchfluss- und Heizterme

Die Nettoenergie, die durch den Prozessstrom eintritt, beträgt

[ FC_p(T_i-T) ]

Der erste Teil, (FC_pT_i), stellt die Energie dar, die mit dem Zulauf eintritt. Der zweite Teil, (FC_pT), stellt die Energie dar, die mit dem ideal durchmischten Auslass austritt.

Die Addition der Dampfschlangenleistung ergibt die vollständige Bilanz:

[ MC_p\frac{dT}{dt}

FC_p(T_i-T)+UA(T_{st}-T) ]

Umstellen des dynamischen Modells

Erweitern der Wärmeübertragungsterme

Das Erweitern der rechten Seite ergibt

[ MC_p\frac{dT}{dt}

FC_pT_i-FC_pT+UAT_{st}-UAT ]

Das Zusammenfassen der Tanktemperaturterme führt zu

[ MC_p\frac{dT}{dt}+(FC_p+UA)T

FC_pT_i+UAT_{st} ]

Normalisieren der Gleichung

Die Division durch (MC_p) ergibt

[ \frac{dT}{dt} + \left( \frac{F}{M}+\frac{UA}{MC_p} \right)T

\frac{F}{M}T_i + \frac{UA}{MC_p}T_{st} ]

Definieren Sie

[ a=\frac{F}{M}+\frac{UA}{MC_p} ]

[ b_i=\frac{F}{M} ]

[ b_{st}=\frac{UA}{MC_p} ]

Das Modell wird zu

[ \frac{dT}{dt}+aT=b_iT_i+b_{st}T_{st} ]

Der Koeffizient (a) bestimmt die Geschwindigkeit, mit der sich der Tank einem neuen Gleichgewicht nähert. Sein Kehrwert,

[ \tau=\frac{1}{a}

\frac{MC_p}{FC_p+UA} ]

ist die effektive thermische Zeitkonstante.

Umwandlung in Abweichungsvariablen

Definition des stationären Betriebspunkts

Seien die anfänglichen stationären Werte (T_s), (T_{i,s}) und (T_{st,s}). Im stationären Zustand gilt

[ \frac{dT}{dt}=0 ]

sodass die stationäre Bilanz lautet

[ aT_s=b_iT_{i,s}+b_{st}T_{st,s} ]

Dieser Betriebspunkt ist die Referenz, die zur Untersuchung von Störungen verwendet wird.

Definition der Temperaturabweichungen

Definieren Sie Abweichungsvariablen als

[ T'=T-T_s ]

[ T_i'=T_i-T_{i,s} ]

[ T_{st}'=T_{st}-T_{st,s} ]

Da angenommen wird, dass die Systemparameter konstant sind, ergibt die Subtraktion der stationären Gleichung von der dynamischen Gleichung

[ \frac{dT'}{dt}+aT'

b_iT_i'+b_{st}T_{st}' ]

Dieser Schritt entfernt die konstanten Betriebspunktterme und lässt nur die Reaktion auf Änderungen der beiden Eingänge übrig.

Anwendung der Laplace-Transformation

Transformation der Differentialgleichung

Nehmen Sie eine anfängliche Abweichung von Null an, was bedeutet, dass der Tank bei seiner stationären Temperatur beginnt. Die Anwendung der Laplace-Transformation ergibt

[ sT'(s)+aT'(s)

b_iT_i'(s)+b_{st}T_{st}'(s) ]

Das Ausklammern der Tanktemperatur-Transformierten ergibt

[ (s+a)T'(s)

b_iT_i'(s)+b_{st}T_{st}'(s) ]

Daher gilt

[ T'(s)

\frac{b_i}{s+a}T_i'(s) + \frac{b_{st}}{s+a}T_{st}'(s) ]

Identifizierung der beiden Übertragungsfunktionen

Die Übertragungsfunktion von der Einlassflüssigkeitstemperatur zur Tanktemperatur lautet

[ G_i(s)

\frac{T'(s)}{T_i'(s)}

\frac{b_i}{s+a}

\frac{F/M}{s+F/M+UA/(MC_p)} ]

wobei die Abweichung der Dampftemperatur auf Null gehalten wird.

Die Übertragungsfunktion von der Dampftemperatur zur Tanktemperatur lautet

[ G_{st}(s)

\frac{T'(s)}{T_{st}'(s)}

\frac{b_{st}}{s+a}

\frac{UA/(MC_p)}{s+F/M+UA/(MC_p)} ]

wobei die Abweichung der Einlassflüssigkeitstemperatur auf Null gehalten wird.

Die vollständige Eingangs-Ausgangs-Beziehung ist daher

[ T'(s)=G_i(s)T_i'(s)+G_{st}(s)T_{st}'(s) ]

Physikalische Interpretation des Modells

Reaktion auf die Einlasstemperatur

Der Pfad der Einlasstemperatur hat den Zähler (F/M). Eine höhere Durchflussrate bewirkt, dass Einlassstörungen den Tank stärker beeinflussen und reduziert die Verweilzeit der Flüssigkeit.

Im stationären Zustand ist die Verstärkung der Einlasstemperatur

[ K_i=G_i(0)=\frac{b_i}{a} ]

oder

[ K_i

\frac{F/M}{F/M+UA/(MC_p)} ]

Reaktion auf die Dampftemperatur

Der Pfad der Dampftemperatur hat den Zähler (UA/(MC_p)). Eine Vergrößerung der Wärmeübertragungsfläche oder des Wärmedurchgangskoeffizienten erhöht den Einfluss von Änderungen der Dampftemperatur auf den Tank.

Die stationäre Verstärkung ist

[ K_{st}=G_{st}(0)=\frac{b_{st}}{a} ]

oder

[ K_{st}

\frac{UA/(MC_p)}{F/M+UA/(MC_p)} ]

Für dieses vereinfachte Modell gilt

[ K_i+K_{st}=1 ]

Dies spiegelt die Tatsache wider, dass die stationäre Tanktemperatur ein gewichtetes Gleichgewicht zwischen der Einlassflüssigkeitstemperatur und der Dampftemperatur darstellt.

Zeitkonstante und Polstelle

Beide Übertragungsfunktionen haben denselben Nenner und daher dieselbe Polstelle:

[ s=-a ]

Die gemeinsame Zeitkonstante ist

[ \tau=\frac{1}{a} ]

Unter Verwendung von (\tau) können die Übertragungsfunktionen auch geschrieben werden als

[ G_i(s)=\frac{K_i}{\tau s+1} ]

[ G_{st}(s)=\frac{K_{st}}{\tau s+1} ]

Die gemeinsame Polstelle bedeutet, dass beide Eingänge durch denselben thermischen Speicher des Tanks gefiltert werden. Ihre stationären Auswirkungen unterscheiden sich, aber ihre grundlegende Reaktionsgeschwindigkeit erster Ordnung ist unter den genannten Annahmen identisch.

Vorhersage von Sprungantworten

Sprungänderung der Einlasstemperatur

Für eine Sprungänderung der Größe (\Delta T_i) in der Einlasstemperatur, bei unveränderter Dampftemperatur,

[ T'(t)

K_i\Delta T_i \left(1-e^{-t/\tau}\right) ]

Die Tanktemperatur nähert sich exponentiell ihrem neuen stationären Wert.

Sprungänderung der Dampftemperatur

Für einen Dampftemperatursprung der Größe (\Delta T_{st}), bei unveränderter Einlasstemperatur,

[ T'(t)

K_{st}\Delta T_{st} \left(1-e^{-t/\tau}\right) ]

Die Form ist identisch, aber die Verstärkung wird durch die relative Stärke der Wärmeübertragung der Heizschlange bestimmt.

Gleichzeitige Änderungen

Wenn sich beide Temperaturen gleichzeitig ändern,

[ T'(t)

K_i\Delta T_i \left(1-e^{-t/\tau}\right) + K_{st}\Delta T_{st} \left(1-e^{-t/\tau}\right) ]

Für allgemeinere zeitabhängige Eingänge wird das Laplace-Modell direkt verwendet, oder die ursprüngliche Differentialgleichung kann numerisch gelöst und mit Sensordaten der Pilotanlage verglichen werden.

Verknüpfung der Parameter mit Messungen der Pilotanlage

Schätzung der Wärmeübertragungskapazität

Die Größe (UA) steuert, wie effektiv Änderungen der Dampftemperatur auf den Tank übertragen werden. Sie kann geschätzt werden aus

[ Q=UA\Delta T_m ]

wobei (\Delta T_m) die geeignete mittlere Temperaturdifferenz ist, oft die logarithmische mittlere Temperaturdifferenz für Berechnungen an Wärmetauschern.

Der gemessene Wert von (U) beinhaltet Konvektion, Wandleitung und Verschmutzungswiderstände.

Bewertung des Rühreffekts

Rühren verbessert die Durchmischung und kann den wärmeseitigen Wärmeübergangskoeffizienten erhöhen. Es unterstützt auch die Annahme der idealen Durchmischung, die im Modell verwendet wird.

Schlechte Durchmischung kann dazu führen, dass gemessene Reaktionen vom Ein-Temperatur-Modell abweichen, da verschiedene Bereiche des Tanks unterschiedlich schnell reagieren können.

Sorgfältiger Umgang mit dimensionslosen Parametern

Das Verhältnis

[ B=\frac{UA}{FC_p} ]

vergleicht die Wärmeübertragungskapazität der Schlange mit der Wärmekapazitäts-Durchflussrate des Prozessstroms.

Ein großes (B) zeigt an, dass die Dampfheizung im Verhältnis zu den Auswirkungen des Flüssigkeitsdurchsatzes stark ist. Ein kleines (B) zeigt an, dass der Einlassstrom einen größeren Einfluss auf die Auslasstemperatur hat.

Verständnis der Kompromisse

Ein Modell erster Ordnung ist eine Näherung

Das Modell nimmt konstante (U), (A), (F), (M) und (C_p) an. In einer realen Pilotanlage können diese Parameter mit Durchflussrate, Temperatur, Viskosität, Verschmutzung und Rührgeschwindigkeit variieren.

Die Übertragungsfunktionen sind daher in der Nähe des für die Linearisierung verwendeten stationären Betriebspunkts am zuverlässigsten.

Dampftemperatur entspricht möglicherweise nicht der Wärmeleistung

Das Modell behandelt die Dampftemperatur als Eingang und verwendet

[ Q=UA(T_{st}-T) ]

Dies ist angemessen, wenn der Dampf bei einer annähernd konstanten Temperatur kondensiert und die Dynamik der Heizschlange im Vergleich zur Tankdynamik schnell ist.

Wenn Dampfdruck, Ventilposition, Kondensatabfuhr oder die Temperatur der Schlange eine signifikante Dynamik aufweisen, repräsentiert die Dampftemperatur allein den Heizeingang möglicherweise nicht ausreichend.

Transport- und Sensordynamik können wichtig sein

Die ideale Bilanz hat keine explizite Verzögerung zwischen einer Änderung der Einlasstemperatur und ihrer Auswirkung auf den Tank. Echte Rohrleitungen, Probenahmestellen, Schutzrohre und Temperatursensoren können Transportverzögerungen und zusätzliche Zeitkonstanten einführen.

Diese Effekte sollten einbezogen werden, wenn angepasste Übertragungsfunktionen nicht mit experimentellen Sprungtestdaten übereinstimmen.

Linearisierung begrenzt den gültigen Bereich

Das Modell mit Abweichungsvariablen setzt Änderungen voraus, die klein genug sind, damit die Systemparameter annähernd konstant bleiben. Große Temperatur- oder Durchflussänderungen erfordern möglicherweise eine nichtlineare Energiebilanz oder eine numerische Integration.

Eine zuverlässige Validierung überprüft die Anfangsbedingung, den langfristigen stationären Zustand und ob die vorhergesagte transiente Kurve mit den gemessenen Daten übereinstimmt.

Die richtige Wahl für Ihr Ziel treffen

Verwenden Sie das aus der Bilanz abgeleitete Modell als Ausgangspunkt für Analyse und Regelungsentwurf.

Wenn Ihr Hauptaugenmerk auf Einlasstemperaturstörungen liegt: Verwenden Sie (G_i(s)) und betonen Sie das Verhältnis von Durchfluss zu Füllmenge (F/M) sowie die Verstärkung der Einlass-Sprungantwort (K_i).
Wenn Ihr Hauptaugenmerk auf der Dampftemperaturregelung liegt: Verwenden Sie (G_{st}(s)) und betonen Sie (UA/(MC_p)), Wärmeübertragungsfläche, Rühren und den Gesamtwärmedurchgangskoeffizienten (U).
Wenn Ihr Hauptaugenmerk auf der Validierung der Pilotanlage liegt: Führen Sie Sprungtests durch, schätzen Sie die Verstärkungen und die Zeitkonstante und vergleichen Sie die gemessene Reaktion mit den Vorhersagen erster Ordnung.
Wenn Ihr Hauptaugenmerk auf einer High-Fidelity-Simulation liegt: Behalten Sie die nichtlineare Energiebilanz bei und fügen Sie Transport-, Schlange-, Ventil- und Sensordynamiken hinzu, wo die Daten zeigen, dass das Modell erster Ordnung nicht ausreicht.

Eine Übertragungsfunktion ist die kompakte dynamische Konsequenz einer physikalischen Energiebilanz, die messbare Prozessänderungen mit der Reaktion der Tanktemperatur verknüpft.

Zusammenfassungstabelle:

Element	Ergebnis	Praktische Bedeutung
Dynamisches Modell	`MCp(dT/dt) = FCp(Ti - T) + UA(Tst - T)`	Beschreibt Tankwärmeakkumulation, Einlassstrom und Dampfheizung
Übertragungsfunktion Einlasstemperatur	`Gi(s) = (F/M) / [s + F/M + UA/(MCp)]`	Sagt die Tankreaktion auf Änderungen der Einlassflüssigkeitstemperatur voraus
Übertragungsfunktion Dampftemperatur	`Gst(s) = [UA/(MCp)] / [s + F/M + UA/(MCp)]`	Sagt die Tankreaktion auf Änderungen der Heizdampftemperatur voraus
Stationäre Verstärkung Einlass	`Ki = (F/M) / [F/M + UA/(MCp)]`	Misst den Einfluss der Einlasstemperatur auf die Tanktemperatur
Stationäre Verstärkung Dampf	`Kst = [UA/(MCp)] / [F/M + UA/(MCp)]`	Misst den Heizeinfluss der Dampftemperatur
Gemeinsame Zeitkonstante	`tau = MCp / (FCp + UA)`	Gibt an, wie schnell sich der Tank einer neuen Temperatur nähert
Wichtige Validierungsfaktoren	`U`, `A`, `F`, `M`, `Cp`, Mischung, Verzögerung und Sensordynamik	Bestimmt, wie genau das Modell mit den Daten der Pilotanlage übereinstimmt

Modellieren und validieren Sie Ihre Pilotanlage mit LABPARK

Benötigen Sie zuverlässige dynamische Daten für die chemische Verfahrenstechnik oder die Forschung zur Prozessregelung? Kontaktieren Sie LABPARK, um eine für Ihre Anwendung entwickelte Pilotanlage für Ausbildung oder Berufsbildung zu besprechen.

LABPARK unterstützt Universitäten, Forschungsinstitute und Unternehmen mit Pilotanlagen für chemische Verfahrenstechnik, Bioprozesse und Biotechnologie sowie Umwelt- und Wasseraufbereitung. Unser Rührkesselheizer und die zugehörige Ausrüstung für Verfahrenstechnik helfen Ihnen dabei:

Energiebilanzen, Übertragungsfunktionen, Prozessdynamik und Regelungsstrategien zu demonstrieren.
Wiederholbare Sprungtests für Einlass- und Dampftemperaturen durchzuführen.
Wärmedurchgangskoeffizienten, Verweilzeit, Rührwirkung und thermische Reaktion zu bewerten.
Praktische Daten für Lehre, Forschung, Prozessentwicklung und Bedienerschulung zu generieren.
Eine konfigurierbare Pilotplattform aufzubauen, die auf Ihre Labor-, Lehrplan- oder Produktionsanforderungen abgestimmt ist.

Teilen Sie uns Ihre erforderliche Kapazität, Betriebstemperatur, Instrumentierung sowie Lehr- oder Forschungsziele mit. LABPARK kann Ihnen helfen, die richtige Pilotanlagenkonfiguration zu identifizieren und praktische technische Unterstützung von der Spezifikation bis zur Anwendung zu bieten.