CH686766A5

CH686766A5 - Temperaturregelung einer Strangpresse.

Info

Publication number: CH686766A5
Application number: CH00561/93A
Authority: CH
Inventors: Madhukar Pandit; Karlheinz Buchheit
Original assignee: Alusuisse Lonza Services Ag; Mannesmann Ag
Priority date: 1993-02-24
Filing date: 1993-02-24
Publication date: 1996-06-28
Also published as: EP0615795B1; DE59402183D1; JPH06277750A; CZ284315B6; CZ41294A3; US5614137A; EP0615795A1; CA2119941A1

Description

5

10

15

20

25

30

35

40

45

50

55

60

CH 686 766 A5

Beschreibung

Vorliegende Erfindung betrifft ein Verfahren für die Regelung von Strangpressen und die Verwendung des Verfahrens zur Herstellung von Strangpressprofilen.

Das Strangpressen ist ein bekanntes und vielseitig verwendbares Verfahren zum Herstellen von Profilen mittels Extrudieren von Werkstoffen, wie beispielsweise Metall, Glas oder Kunststoff durch eine Matrize, wobei die Matrizenöffnung einen beinahe beliebigen Querschnitt von kreisrund bis hin zu verwinkelten Mustern mit einem oder mehreren Hohlräumen aufweisen kann.

Eine Strangpressvorrichtung enthält im wesentlichen einen Aufnehmer mit einer zylinderförmigen Ausnehmung beliebigen Querschnittes, der das Pressgut, üblicherweise in Form eines zylinderförmigen Blockes, sowie einen mit einer Pressscheibe versehenen Stempel aufnehmen kann, wobei an der einen Öffnung der zylinderförmigen Ausnehmung des Aufnehmers eine Matrize angebracht werden kann.

Bei der Herstellung von Strangpressprofilen wird das Pressgut in die zylinderförmige Ausnehmung des Aufnehmers geführt und unter Anwendung eines über die Pressscheibe vermittelten hohen axialen Druckes an die Matrize gedrückt, so dass sich das Pressgut unter den gegebenen Temperaturen plastisch verformen kann und somit durch die Matrizenöffnung extrudiert wird.

Beim Strangpressen von kristallinen oder glasartigen Werkstoffen entspricht der Querschnitt der entstehenden Strangpressprofile im wesentlichen dem Querschnitt der Matrizenöffnung. Dies gilt jedoch nicht für das Strangpressen von Polymeren mit strukturviskosen (Abnahme der Viskosität mit zunehmender mechanischer Beanspruchung), entropie-elastischen (Strangaufweitung) und viskoelastischen (zeitabhängige Kopplung der Viskosität und Elastizität) Eigenschaften.

Die plastische Verformbarkeit des Pressgutes und damit die durch die Matrize extrudierte Pressgutmenge pro Zeiteinheit hängt neben der Materialzusammensetzung des Pressgutes und des angewendeten Druckes vorwiegend von der Prozesstemperatur ab. Um bei diesem Warmumformprozess eine möglichst hohe Ausstossgeschwindigkeit zu erreichen, wird die Strangaustrittstemperatur möglichst hoch gehalten. Die maximal zulässige Strangaustrittstemperatur liegt einerseits unterhalb der Schmelztemperatur des Pressgutes und ist andererseits durch die Bedingung gegeben, dass sich der aus der Matrizenöffnung austretende Strang im heissen Zustand unter seinem Eigengewicht nicht verformen darf. Des weiteren hat die Strangaustrittstemperatur einen erheblichen Einfluss auf die stofflichen Eigenschaften der Strangpressprofile und damit auf die Produktequalität (Homogenität, mechanische Spannungen usw.). Deshalb besteht auch aus Gründen der Qualitätssicherung ein erhebliches Interesse, die Strangaustrittstemperatur definiert vorzugeben und während dem Prozess konstant zu halten. Ein solches Verfahren mit einer vordefinierten und konstant gehaltenen Strangaustrittstemperatur wird als isothermes Strangpressen bezeichnet.

Die Bilanz des Energie-Haushaltes ergibt sich aus der Differenz aller zugeführten Energie (mechanische Arbeit und Wärme) und der abgeführten Energie (plastische Verformung, Wärmeleitung). Dabei bezieht sich der für den Warmumformprozess wesentliche Energie-Haushalt auf den Teil des Pressgutblockes, der plastisch verformt wird. Die resultierende Temperatur der Profile beim Austritt aus der Matrize kann beispielsweise durch die Vorheiztemperatur der Pressgutblöcke oder Barren und die Prozessgeschwindigkeit gezielt beeinflusst werden.

Die praktische Realisierung des isothermen Strangpressens bedingt jedoch die vollständige Kenntnis und Beherrschung aller Prozessparameter und insbesondere sämtlicher thermischer Prozessgrössen, weswegen dieses Verfahren mit erheblichen, technologisch unzureichend gelösten Problemen verbunden ist. Solche Probleme können durch Anwendung bekannter Methoden der Regelungstechnik wie beispielsweise das simulierte bzw. geregelte isotherme Strangpressen verhindert werden.

Beim simulierten Strangpressen wird die Strangaustrittstemperatur mittels einem Simulationsmodell vorausberechnet, wobei die Pressgeschwindigkeit der regelungstechnisch relevante Prozessparameter darstellt. Der Strangpressprozess stellt jedoch ein kompliziertes thermomechanisches System mit vielen nicht leicht kontrollierbaren Parametern dar, so dass der Gesamtprozess des Strangpressens analytisch nicht vollständig und mit numerischen Methoden nur ungenau beschrieben werden kann. Deshalb eignet sich dieses Verfahren für die Regelung einer Strangpresse nicht.

Beim geregelten Strangpressen wird das Herstellen und Aufrechterhalten der als Regelgrösse bezeichneten Strangaustrittstemperatur mittels einem geschlossenen Regelkreis erreicht, der durch permanenten Vergleich von Soll- und Ist-Wert der Regelgrösse, die zur Korrektur notwendige Pressgeschwindigkeit als Stellgrösse berechnet. Zur Messung der Strangaustrittstemperatur werden dabei üblicherweise Strahlungspyrometer eingesetzt.

Die pyrometrische Temperaturmessung geschieht dabei unter Ausnutzung der Planck'schen Strahlungsgesetze, die jedoch nur für ideale schwarze Körper gelten. Kennt man die Gesamtenergie der ausgesandten Strahlung, so kann aus der Messung der Energie eines gewissen Spektralbereichs unter Zuhilfenahme der Planck'schen Strahlungsgesetze die Temperatur berechnet werden, die der Körper hätte, wenn er ein schwarzer Körper wäre. Da die meisten Körper nicht ideal schwarz sind, ist die wahre Temperatur höher als die auf diese Weise errechnete. Um die Temperatur eines realen Körpers zu bestimmen, muss die Emissivität, d.h. das Strahlungsvermögen des betrachteten Körpers, bekannt sein. Die Emissivität eines opakten Körpers wird durch den Quotienten der emittierten Strahlung des Körpers und der Strahlung eines schwarzen Körpers mit derselben Temperatur definiert. Die Emissivität kann

2

5

10

15

20

25

30

35

40

45

50

55

60

65

CH 686 766 A5

physikalisch durch einen auf die Planck'schen Strahlungsgesetze multplikativ wirkenden Emissionsgrad e beschrieben werden. Ein idealer schwarzer Körper hat den Emissionsgrad e gleich eins.

Die kontaktlose, pyrometrische Temperaturmessung führt jedoch an Materialien mit kleinem und/oder wellenlängenabhängigem Emissionsgrad (e < 0,1) und/oder veränderlichen Oberflächencharakteristiken, wie beispielsweise an Materialien aus Aluminium oder Aluminiumlegierungen oft zu einer ungenauen Temperaturbestimmung. Deshalb eignet sich das geregelte Strangpressen für solche Materialien nicht.

Aufgabe vorliegender Erfindung ist es, ein Verfahren zu schaffen, das obige Nachteile überwindet und die präzise Steuerung einer Strangpresse zur Erreichung eines maximalen Durchsatzes bei gleichzeitig optimaler Qualität der Strangpressprofile erlaubt.

Erfindungsgemäss wird dies dadurch erreicht, dass die Pressgeschwindigkeit v(t) von Strangpressen derart geregelt wird, dass die Strangaustrittstemperatur ea(t) möglichst konstant und gleich einem vorgegebenen Sollwertverlauf eaw(t) ist und a) die Temperaturregelung zyklisch erfolgt;

b) der zeitliche Verlauf der Pressgeschwindigkeit v«(t) und der Strangaustrittstemperatur ea«(t) während jedem Zyklus k erfasst wird;

c) die Abhängigkeit der Strangaustrittstemperatur eaK(t) von der Pressgeschwindigkeit v«(t) während des ganzen Zyklus k bestimmt wird;

d) mit Hilfe dieser Abhängigkeit und dem zeitlichen Verlauf von v«(t) der Verlauf der Pressgeschwindigkeit VK+i(t) für den nächsten Zyklus k+1 derart bestimmt wird, dass der Regelfehler

(1) ek+l(t) = 03w(t) — 03k+l(t)

und der Stellaufwand

(2) dvk+i(t) = Vk+1 (t) - vk(t)

möglichst gering werden, wobei der Sollwertverlauf 9aw(t) für jeden Zyklus k definiert werden kann;

e) Stellgrössenbeschränkungen vmin, k < Vk(t) < vmax, k berücksichtigt werden können;

f) die Pressgeschwindigkeit vk+i(t) vor Beginn des Presszyklus k+1 berechnet wird;

g) Vk+i(t) während dem Presszyklus k+1 nicht verändert wird;

h) nach Beendigung des Presszyklus k+1 die Verfahrensschritte b) - g) in rekursiver Weise für jeden weiteren Presszyklus wiederholt werden bis der Pressvorgang beendet ist.

Dieses erfindungsgemässe Verfahren unterscheidet sich von den bekannten Festwertregelungen dadurch, dass nicht wie bei einem geschlossenen Regelkreis nur der lokale Arbeitspunkt, sondern stets der gesamte Zyklus optimiert wird. Wegen der Zyklizität, d.h. der repetitiven Natur des Regelprozesses, wird die Erfahrung aus einem Zyklus k für die Generierung der Pressgeschwindigkeitskurve k+1 verwendet, womit eine Rückkopplung von einem auf den anderen Zyklus vorhanden ist. Damit ist dieses Regelverfahren weniger störanfällig auf die pyrometrische Messung der Strangaustrittstemperatur und eignet sich bevorzugt für die Temperaturregelung von Strangpressen zur Herstellung von Strangpressprofilen mit kleinem und/oder wellenlängenabhängigem Emissionsgrad (e < 0,1) und/oder veränderlichen Oberflächencharakteristiken und somit insbesondere für die Herstellung von Strangpressprofilen aus Aluminium oder Aluminiumlegierungen.

Beim Strangpressen von Aluminium oder dessen Legierungen wird das Pressgut in einem Ofen auf eine Temperatur von 400 bis 500°C aufgeheizt und anschliessend in einen Aufnehmer oder Rezipienten geladen. Dieser ist einseitig von einer Matrize verschlossen, deren Öffnung oder Durchbruch dem Querschnitt des entstehenden Profilstranges entspricht. Von der der Matrize gegenüberliegenden Seite des Pressgutes wird mit einem Stempel unter Einwirkung eines hohen Druckes von mehr als 10 MN (Mega Newton) der Pressgutblock bis auf einen kleinen Rest durch die Matrize gedrückt. Nach Beendigung eines Zyklus wird ein neuer Block in den Rezipienten geladen und der Prozessvorgang kann wiederholt werden.

Zur Verdeutlichung des Strangpressprozesses zeigt Fig. 1 die für den Pressbetrieb wesentlichen Komponenten einer Strangpresse sowie die während dem Prozess auftretenden Wärmeflüsse.

Unter regelungstechnischen Gesichtspunkten gilt für eine Strangpresse mit einem Strahlungspyrometer als Messgerät für die Regelgrösse folgendes Systemverhalten:

- Die Sollwertkurve der Strangaustrittstemperatur eaw(t) ist vor Beginn eines Zyklus bekannt.

- Die Dauer eines Zyklus Tzyk hat stets die gleiche Grössenordnung, wobei die Zyklusdauer Tzyk je nach Maschinentyp, Werkzeug und Legierung zwischen 60 bis 1000 s variieren kann. Bei der Verwendung derselben Maschine, desselben Werkzeugs und derselben Legierung ist die Systemänderung während eines Zyklus jedoch im wesentlichen auf +/- 20% beschränkt.

- Das thermische Systemverhalten ist nur langsam zeitveränderlich und wird im wesentlichen durch den Rezipienten bestimmt, dessen thermische Zeitkonstante typischerweise zwischen 3 und 5 Stunden beträgt.

- Der Prozess ist nicht linear und mit analytischen Methoden kaum beschreibbar.

- Das Prozessverhalten ist determiniert, d.h. relevante Prozessparameter wie beispielsweise die Rezipienten-, Werkzeug- oder Blocktemperatur oder die geometrischen Abmessungen des Rezipienten

3

5

10

15

20

25

30

35

40

45

50

55

CH 686 766 A5

bzw. Werkzeuges ändern sich nicht willkürlich; somit unterliegt der Prozess keinen stochastischen Parameterschwankungen und ist stets reproduzierbar.

- Jeder Zyklus hat denselben Anfangszustand.

- Die hier betrachtete Stellgrösse des Prozesses (Pressgeschwindigkeit) ist betragsmässig begrenzt und nicht beliebig schnell veränderbar.

- Die Erfassung der Regelgrösse (Strangaustrittstemperatur ea) ist mit erheblichen Fehlern, Messstörungen und mit einer Totzeit (verzögerte Reaktion) behaftet, weswegen eine offline-Verarbeitung der Messsignale vorteilhaft ist. Im Gegensatz zu einer online-Verarbeitung der Messwerte, bei der die Auswertung der gemessenen Grössen während dem Pressvorgang vorgenommen wird, geschieht die Verarbeitung der Messgrössen beim offline-Verfahren in der Nebenzeit zwischen zwei Presszyklen.

Die Struktur des erfindungsgemässen Verfahrens wird durch das in Fig. 2 dargestellte Wirkungsschaltbild verdeutlicht, welches eine zyklische Regeleinrichtung darstellt, die das Herstellen und Aufrechterhalten einer möglichst konstanten und einem vorgegebenen Sollwertverlauf 0aw(t) entsprechenden Strangaustrittstemperatur ea(t) ermöglicht. Die Regeleinrichtung ist dabei der beeinflussende Teil und die Regelstrecke der beeinflusste Teil des Regelkreises. Nach Ablauf eines Presszyklus wird aus dem Verlauf der Pressgeschwindigkeit Vk(t) und der Austrittstemperatur eak(t) die Stellgrösse berechnet. Dies geschieht durch Identifikation oder Berechnung einer Sprungantwort hk(t) der Regelstrecke, wobei 0 < t < Tzyk-

Unter Identifikation wird allgemein die Berechnung bzw. Abschätzung von Parametern eines vorgegebenen Gleichungssystems verstanden wie dies beispielsweise für die Bestimmung der Koeffizienten von Differentialgleichungen oder zur Berechnung von Stützstellen der nachfolgend vorgeschlagenen Sprungantwort benützt wird. In einem Optimierungsverfahren wird dann aus der Sprungantwort hk(t) und dem Regelfehler ek(t) eine Korrekturkurve oder Korrekturtrajektorie dvk+i berechnet und zur Bahnkurve vk(t) hinzuaddiert. Die so entstehende Kurve Vk+i(t) wird in einem Speicher abgelegt, aus dem sie während dem nächsten Zyklus abgerufen werden kann.

Das erfindungsgemässe Verfahren erleichtert die Unterdrückung von Messstörungen, da im Gegensatz zu den bekannten Regelungen leistungsfähige, nichtkausale Filter eingesetzt werden können. Dabei ist der Ausgangswert y(to) eines nichtkausalen Filters zu einem Zeitpunkt to nicht nur - wie bei einem kausalen Filter - von Eingangswerten x(to-At) mit At > 0, sondern auch von Eingangswerten x(to+At) abhängig. Damit führt das erfindungsgemässe Verfahren trotz schwierigen Randbedingungen zu einem sicheren und bezüglich der Messwerte robusten Regelsystem.

Aufgrund der thermischen Trägheit der Strangpresse sind die Systemveränderungen, wie beispielsweise die Werkzeug-, Rezipienten-, Block- oder Stempeltemperatur aufeinanderfolgender Zyklen vernachlässigbar klein, wodurch die zyklische Regelung solchen Veränderungen schnell genug folgen kann, so dass ein optimaler Prozessverlauf gewährleistet ist. Zudem trägt die Identifikation der Regelstrecke zur Erhöhung der Konvergenzgeschwindigkeit bei, so dass bereits nach wenigen Zyklen der stationäre Presszustand erreicht wird.

Die Erfassung und Verarbeitung der Messgrössen geschieht üblicherweise mittels Datenerfassungsgeräten mit beschränkter Rechenkapazität, wie beispielsweise Mikrorechnern. Um den Rechenaufwand für die zyklische Regelung zu beschränken, werden deshalb die Zeitfunktionen (Strangaustrittstemperatur, Pressgeschwindigkeit) diskret abgetastet.

Eine Möglichkeit, das erfindungsgemässe Verfahren zweckmässig durchzuführen, besteht darin, dass a) das kontinuierliche Zeitverhalten in diskrete Zeitintervalle Ta

(3) t = i*TA, i = 0, 1, 2 ...

unterteilt wird;

b) finite Zustandsänderungen für die Pressgeschwindigkeit und die Strangaustrittstemperatur verwendet werden;

c) zur Begrenzung des Rechenaufwandes und zur Dämpfung des Regelsystems der Verlauf der Pressgeschwindigkeit nicht zu jedem beliebigen Zeitpunkt veränderbar ist, sondern abschnittweise während einem Zeitintervall j von der Dauer itTTa konstant ist, wobei j = 0, 1, 2, ... n-1, n und m eine natürliche Zahl ist, so dass für jeden Zyklus gilt: i = 0, 1,2, .., n*m-1;

d) der Pressgeschwindigkeitsverlauf nach Gleichung (4) durch Elementarfunktionen dargestellt wird ri - l

(4) vk(i*TA) =2 Avk . * a( ( i-j*m)*TÄ )

: _ n * ' J

und

4

5

10

15

20

25

30

35

40

45

50

55

60

65

CH 686 766 A5

(6) AVk, i = Vk(j*m*TA) - Vk((i*m-1)*TA)) mit j = 0, 1, 2, n-1

die Sprunghöhen des Pressgeschwindigkeitsverlaufes für den Zeitpunkt j*m*TA bezeichnen;

e) unter der Voraussetzung von Linearität und Zeitinvarianz, Annahmen, die in der Umgebung der Ar-beitstrajektorien zulässig sind, die Strangaustrittstemperatur nach Gleichung (7) dargestellt wird:

n-1

(7)0ak(i*TA) = 2 Avk>j * M(i-j*m)*TA),

j-0

wobei h(i*TA) die Reaktion der Strangpresse auf einen Einheitssprung o(i*TA) ist;

f) durch Invertierung von Gleichung (7) und nach Messung der Verläufe eak(i*ÏA) und Vk(i*TA) die Sprungantwort h(i'TA) identifiziert wird

0ak(i*TA) - \ Avk)j*hk((i-j*m)*TAj

1

(8) hk(i*TA) = * (

j-0

für l*m < i < (1+1 )*m, wobei wegen der Kausalität des Systems

(9) hk(i Ta) = 0 für i < 0;

g) der Verlauf der Pressgeschwindigkeitskurve Vk+i(i*TA) aus dem rekursiven Regelgesetz (10) gebildet wird:

(10) Vk+i(i*TA) = Vk(i*TA) + dvk+i(i*TA)

und

(11) eak+i(i*TA) = eak(i*TA) + deak+i(i*TA)

gilt;

h) durch Minimierung der Gütefunktion Q

(12)

Min f n " 1 1 n - 1 ^

ûdvk+ltJ Q= - * 2 (Advktl>j) + — * 2 [ek(i*TA)-d0akt;(i*TA)]2

n*nji*0 J

bezüglich der Stellgrössensprünge Advk+i, j der optimale Pressgeschwindigkeitsverlauf gefunden wird, wobei X eine frei wählbare Gewichtung beschreibt und

(13) ek(i*TA) = eaw(i*TA) - eak(i*TA)

der gemessene Regelfehler aus dem gerade vergangenen Zyklus k und n

(14) d0ak + 1(i*TA) = 2 ûdvk + 1 j*hk((i-j*m)*TA)

j -°

die durch Einwirkung von Advk+i, j vorausberechnete Änderung des Temperaturverlaufes deak+i(i*TA) ist;

und i) Steligrössenbeschränkungen von der Form j

(15) Ï 4vk,r < j=0, i, 2,.. ,n-l

5

10

15

20

25

30

35

40

45

50

55

60

65

CH 686 766 A5

j >

<16> 5 ûvk r 4 vBin. 3 >0,1,2,.. ,n-l r-0

berücksichtigt werden können.

Eine schematische Darstellung des Pressgeschwindigkeitsverlaufes eines Presszyklus k findet sich in Fig. 3. Der Zähler i stellt dabei den Index der diskreten Zeitintervalle Ta dar und j den Zählindex für die Stellgrösse v(t), welche jeweils während mindestens der Dauer iti*Ta konstant ist, und die Änderung der Stellgrösse mit aV| bezeichnet wird.

Unter der Voraussetzung von Zeitinvarianz gilt für ein System, das auf einen beliebigen Eingangsverlauf einer Systemgrösse x(t) mit einem Ausgangsverlauf y(t) reagiert, der Verlauf der Ausgangsgrösse y*(t) = y(t+x) für x*(t) = x(t+t).

Die Zeitinvarianz des hier betrachteten Systems ergibt sich aus der Konstanz der Parameter. Somit kann unter der Voraussetzung von Linearität und Zeitinvarianz in der Umgebung der Arbeitstrajektorien oder Bahnkurve vk(t) und 9k(t), d.h. in der Umgebung

(17) Vk+i(t) = Vk(t) + dvk+1 (t)

(18) 0k+i (t) = 8k(t) + dek+i(t),

die Strangaustrittstemperatur durch Gleichung (7) dargestellt werden. Dies obwohl das Systemverhalten der Strangpresse im Prinzip nichtlinear ist; für kleine Änderungen der Stellgrösse Vk(t) jedoch näherungsweise als linear betrachtet werden kann, so dass der Fehler des linearisierten Modells vernachlässigbar gering ist. Das durch Gleichung (7) beschriebene Systemverhalten ergibt durch Invertierung dieser Gleichung, d.h. durch Auflösung von Gleichung (7) nach hk(i*TA), das lineare Gleichungssystem (8), womit nach Messung der Verläufe eak(i*TA) und vk(i*TA) die Sprungantwort hk(i*TA) identifiziert werden kann. Der Wert I in Gleichung (8) kann auch durch (n*m-1) ersetzt werden, da die Summanden für j > I identisch Null sind. Wegen der Kausalität des Systems, d.h. das System reagiert gemäss Gleichung (9) auf eine Änderung der Pressgeschwindigkeit erst nachdem diese Änderung stattgefunden hat, kann der Verlauf der Pressgeschwindigkeitskurve Vk+i(i*TA) bzw. der Strangaustrittstemperatur eak+i(i*TA) nach dem rekursiven Regelgesetz (10) bzw. (11) berechnet werden.

Die gesuchte Grösse ist somit der Stellwertverlauf vk+i(t) für den Presszyklus k+1, wobei der Verlauf Vk(t) aus dem vorherigen Zyklus bekannt ist, so dass dvk+i(t) mittels der Gleichungen (4) und (10) gemäss Gleichung (19) dargestellt werden kann.

n - 1

(19) dvktl(i*TA) = 2 ûdvk+l j * a((i-j*ra)*TA)

j-0

Die gemittelten Veränderungen der Stell- und Regeigrössen werden somit durch die Gütefunktion Q gemäss Gleichung (12) beschrieben, welche entsprechend dem erfindungsgemässen Verfahren minimiert werden soll.

Typische Werte der Parameter für das erfindungsgemässe Verfahren liegen für Tzyk zwischen 60 und 1000 s, für Ta zwischen 0,5 und 3 s, für m zwischen 10 und 20 und diejenigen für n zwischen 10 und 15. Der Wert für den Gewichtsfaktor X liegt typischerweise bei 0,05*m*h((n*m-1)*TA), wobei h((n*m-1)*Ta) der stationäre Endwert der Streckensprungantwort bedeutet.

Wenn keine Pressgeschwindigkeitsbeschränkungen gelten, kann die Minimierung der Gütefunktion Q in Gleichung (12) mittels dem Gradienten-, dem konjugierten Gradienten-, dem Quasi-Newton-, dem Newton-Raphson- oder dem Newton-Verfahren durchgeführt werden.

Wenn hingegen Stellgrössenbeschränkungen für die Pressgeschwindigkeit gelten, kann die Minimierung der Gütefunktion Q in Gleichung (12) bevorzugt durch Verwendung des Kuhn-Tucker-Verfahrens durchgeführt werden.

Die Gütefunktion in Gleichung (12) kann in erfindungsgemässem Verfahren auch durch die Betragsgütefunktion (20) ersetzt werden:

n-1 n * m - 1

(20) Q = x* 2 |ûdvk .| + 2 |ek ( i*TA ) - d0ak+t(i*TA)|

j-0 i - 0

oder durch eine der beiden folgenden gewichteten Gütefunktionen:

6

5

10

15

20

25

30

35

40

45

50

55

60

65

CH 686 766 A5

n - 1 n * m - 1

(21) Q = 2 xj*(Advk+ )2+ 2 Hi*|ek(i*TA) - d0ak+1(i*TA)|2

j-0 J 1-0

ci — i n * ta - 1

(22) Q = 2 Xj*|ûdvk+1 . 1+ 2 M-i*|ek(i*TA) - d0aktl(i*TA)|

j-0 i-0

wobei bzw. m Gewichtsfaktoren darstellen, die vorgängig für jedes Zeitintervall j definiert werden müssen.

Der Gewichtsfaktor X in Gleichung (20) weist typischerweise einen Wert in der Grössenordnung von X = 0,1*m*h((n*m-1)*TA) auf. Die Werte der Gewichtsfaktoren X\ und ^ in Gleichung (21) betragen typischerweise

1

Hi- —*i bzw. Xj = 0,05*|i,( j*m, *h((n*m-l)*TA )

n * m und diejenigen für Gleichung (22), beispielsweise

1

Hi« bzw. Xj =* 0,liJLt j.mJ *h( (n*m-l)*TA)

n * m

Die direkte Berechnung der Sprungantwort nach Gleichung (8), kann zur Erreichung besser dämpfender Eigenschaften bei starken Messstörungen, durch ein Fehlerminimierungsverfahren ersetzt werden, wobei a) eine durch Gleichung (23) definierte Streckenimpulsantwort gk(i*TA), d.h. die Reaktion der Strecke auf einen Impuls nach Gleichung (24), eingeführt wird:

i

(23) 0k ( i*Ta ) = 2 vk (r*TA )*gk ( ( i-r)*TA )

r-0

fl für i=0 (24) 8( i *TA ) = ^0 sonst b) bei der Identifikation zur Begrenzung der Freiheitsgrade nur die ersten N-Werte berücksichtigt werden, so dass für gk(i*TA) gilt:

f*0 für 0<i<N-lì (25) gK ( i *Ta ) = < 0 sonst >

c) die bezüglich der Impulsantwort gk(i*TA) zu minimierende Gütefunktion F somit folgende Gestalt annimmt:

7

5

10

15

20

25

30

35

40

45

50

55

60

CH 686 766 A5

n * m - 1 "

(26) F= 2

6ak(i*TA)-2 vk(r*TA)*gk((i-r)*TA)

r - 0

+ Y[gk(s*TA)-gk((s-l)*TA)]2

s - 0

und d) die Sprungantwort gemäss

(27) MI'TJ - ì gk(r*TJ

r -0

aus dem numerischen Integral der Impulsantwort hervorgeht.

Die Identifikation von Gleichung (23) ist in Gleichung (26) formuliert. Damit wird gk(i*TA) so bestimmt, dass der Modellfehler möglichst klein und die Kurve gk(i*TA) möglichst glatt wird. Die in Gleichung (26) formulierte Gütefunktion F dient somit zur Identifikation der Funktion gk(i'ÏA) und weist keinen Bezug zur Gütefunktion Q auf, wobei letztere unverändert erhalten bleibt.

Der Wert des Parameters N liegt typischerweise bei N = 100-150, wobei er jedoch maximal den Wert n*m-1 annehmen kann.

Die Bestimmung der Sprungantwort kann auch über ein Fehlerminimierungsverfahren im Frequenzbereich (z-Bereich) gefunden werden, wobei a) die Streckenübertragungsfunktion im z-Bereich gemäss Gleichung (28) gegeben ist

N

2 b*z<-r)

/ v 0(z)

(28) G ( Z ) = = ■ wobei

V(z) «

l+2as*z(*s}

s - 0

e(z) und V(z) die z-Transformierten der zeitdiskreten Funktionen 0(i'*Ta) bzw. v(i'*Ta) sind und die Koeffizienten der Übertragungsfunktion as und br durch Minimierung des quadratischen Modellfehlers bestimmt werden;

b) die Rücktransformierte der z-Übertragungsfunktion Gs(z) die Impulsantwort gemäss Gleichung (29) ergibt:

(29) gk(i*TA) = r-i [Gs, k(z)]

und c) die Sprungantwort hiermit nach Gleichung (27) bestimmt wird.

Dieses Verfahren dient zur Minimierung des Modellfehlers:

n * m - 1 2

(30) F = 2 [0ak(i*TA)-0mk(i*TA)]l i-0

wobei emk(i*ÏA) den Modellwert darstellt, so dass

N n

(31) 0mk(i*TA)+ 2 as 0^ ( ( i-s )*TA ) = 2 br*vk((i-r)*TA)

s - 1 r- 1

und N die Modellordnung bezeichnet, welche für das erfindungsgemässe Verfahren typische Werte von 1 < N < 5

annehmen kann. Die in Gleichung (28) mit as und br bezeichneten Grössen sind die Koeffizienten der

8

5

10

15

20

25

30

35

40

45

50

55

60

65

CH 686 766 A5

Differenzengleichung (31). Die z-transformierten Funktionen Gs(z), ©(z) und V(z) in Gleichung (28) sind durch die folgenden Gleichungen (32-34) definiert, wobei z die komplexe Frequenz darstellt.

n * m - 1

(32) Gs (z) = 2 gk( i*Tj*z-i i-0 n * ta - 1

(33) e(z) = 2 0( i*Ta )*z" i i -0 n * m - 1

(34) V(z) = 2 v(i*Ta)*z"i i-0

Die Rücktransformation bedeutet das Auffinden einer diskreten Funktion im Zeitbereich, die die entsprechende z-Funktion als z-Transformierte hat. Im erfindungsgemässen Verfahren bedeutet die Rücktransformation der z-Übertragungsfunktion Gs(z) somit die Berechnung der Impulsantwort gk(i*ÏA), die gemäss Gleichung (23) definiert ist.

Die Messung des Verlaufes der Strangaustrittstemperatur und der Pressgeschwindigkeit sowie deren Auswertung für jeden Presszyklus k und die darauffolgende vorgängige Berechnung des Pressgeschwindigkeitsverlaufes für den nächsten Zyklus k+1 führt zum erfindungsgemässen Verfahren, welches gegenüber den bekannten Regelverfahren wesentlich weniger störanfällig auf die pyrometrische Messung der Strangaustrittstemperatur ist.

Das vorliegende Verfahren erlaubt somit die Temperaturregelung von Strangpressen zur Herstellung von Strangpressprofilen mit kleinem und/oder wellenlängenabhängigem Emissionsgrad (e < 0,1) und/ oder veränderlichen Oberflächencharakteristiken. Bevorzugt wird das Verfahren zur Temperaturregelung von Strangpressen für die Herstellung von hochreflektierenden, metallischen Strangpressprofilen verwendet. Das Verfahren eignet sich damit insbesondere für die Herstellung von Strangpressprofilen aus Aluminium oder Aluminiumlegierungen.

Das erfindungsgemässe Verfahren erlaubt die präzise Steuerung einer Strangpresse und ermöglicht damit die Erreichung eines maximalen Durchsatzes bei gleichzeitig optimaler Qualität der Strangpressprofile und findet weiter überall dort Verwendung, wo die Verarbeitungs- oder Betriebstemperaturen des Messgutes kritisch sind.

Claims

Patentansprüche

1. Verfahren für die Regelung von Strangpressen, dadurch gekennzeichnet, dass die Pressgeschwindigkeit v(t) derart geregelt wird, dass die Strangaustrittstemperatur ea(t) möglichst konstant und gleich einem vorgegebenen Sollwertverlauf eaw(t) ist, und a) die Temperaturregelung zyklisch erfolgt;

b) der zeitliche Verlauf der Pressgeschwindigkeit v«(t) und der Strangaustrittstemperatur eax(t) während jedem Zyklus k erfasst wird;

c) die Abhängigkeit der Strangaustrittstemperatur eax(t) von der Pressgeschwindigkeit vx(t) während des ganzen Zyklus k bestimmt wird;

d) mit Hilfe dieser Abhängigkeit und dem zeitlichen Verlauf von vx(t) der Verlauf der Pressgeschwindigkeit VK+1 (t) für den nächsten Zyklus k+1 derart bestimmt wird, dass der Regelfehler

(1) ek+i(t) = eaw(t) - eak+i(t)

und der Stellaufwand

(2) dvk+1 (t) = Vk+i(t) - Vk(t)

möglichst gering werden, wobei der Sollwertverlauf eaw(t) für jeden Zyklus k definierbar ist;

e) die Pressgeschwindigkeit Vk+i(t) vor Beginn des Presszyklus k+1 berechnet wird;

f) Vk+i(t) während dem Presszyklus k+1 nicht verändert wird;

g) nach Beendigung des Presszyklus k+1 die Verfahrensschritte b) -

h) in rekursiver Weise für jeden weiteren Presszyklus wiederholt werden bis der Pressvorgang beendet ist.

2. Verfahren gemäss Anspruch 1, dadurch gekennzeichnet, dass Stellgrössenbeschränkungen vmin, k < VK(t) < vmax, k berücksichtigt werden.

g

5

10

15

20

25

30

35

40

45

50

55

60

CH 686 766 A5

3. Verfahren nach Anspruch 1, dadurch gekennzeichnet, dass a) das kontinuierliche Zeitverhalten in diskrete Zeitintervalle Ta

(3) t = i*TA, i = 0, 1, 2...

unterteilt wird;

c) zur Begrenzung des Rechenaufwandes und zur Dämpfung des Regelsystems der Verlauf der Pressgeschwindigkeit nicht zu jedem beliebigen Zeitpunkt veränderbar ist, sondern abschnittweise während einem Zeitintervall j von der Dauer m*TA konstant ist, wobei j = 0, 1, 2, n-1, n und m eine natürliche Zahl ist, so dass für jeden Zyklus gilt: i = 0, 1, 2, n*m-1;

d) der Pressgeschwindigkeitsverlauf nach Gleichung (4) durch Elementarfunktionen dargestellt wird n-1

(4) vk(i*TA) =2 av * o((i-j*m)*TA)

j ■ 0

wobei mit o(ì*Ta) die Heaviside oder Einheitssprungfunktion fl für i>0 )

(5) a(i*TA) = jo sonst j-

und

(6) AVk, j = Vk(j*m*ta) - Vk((j*nr-1) Ta)) mit j = 0, 1, 2, n-1

e) unter der Voraussetzung von Linearität und Zeitinvarianz, Annahme, die in der Umgebung der Ar-beitstrajektorien zulässig sind, die Strangaustrittstemperatur nach Gleichung (7) dargestellt wird:

n-1

,7) ea„(in„) ■ 4Vj * hk ( ( i-j*m)*TA ),

wobei h(i*TA> die Reaktion der Strangpresse auf einen Einheitssprung a(i*TA) ist;

f) durch Invertierung von Gleichung (7) und nach Messung der Verläufe eak(i*TA) und vk(i*TA) die Sprungantwort h(i*ta> identifiziert wird

1

(8) = —

AVk , 0

1

6ak(i*Tj - 2

j-0

für l*m < i < (1+1 )*m, wobei wegen der Kausalität des Systems

(9) hk(i*TA) = 0 für k0)

(10) Vk+i (ì*Ta) = Vk(i*ta) + dvk+i(i*ta)

und

(11) ak+1 (ì*Ta) = eak(i*ta) + deak+i(i*ta)

gilt; und h) durch Minimierung der Gütefunktion Q

10

5

10

15

20

25

30

35

40

45

50

55

60

65

CH 686 766 A5

02) X-iJq =-*7(Advktlij)2+— n ï '[eji'Tj-dea^dnj]2}

[nj-0 n * n i - 0 J

(13) ek(i*TA) = eaw(i*TA) - 0ak(i*TA)

der gemessene Regelfehler aus dem gerade vergangenen Zyklus k und n

(14) d0ak+1(i*TA) = 2 Advk+1> j*hk((i-j*m)*TA)

j - 0

die durch Einwirkung von Advk+i, j vorausberechnete Änderung des Temperaturverlaufes deak+i(i*TA) ist.

4. Verfahren nach Anspruch 3, dadurch gekennzeichnet, dass Stellgrössenbeschränkungen von der Form

(15) £ A vk(r £ vmaxj j= 0,1,2 n-1

r»0

(16) ^ A Vj^.r > Vjjjjjjj j—0,1,2,...,n-1

r=. 0

berücksichtigt werden.

5. Verfahren gemäss Anspruch 3, wobei für die Pressgeschwindigkeit keine Stellgrössenbeschrän-kung gilt, dadurch gekennzeichnet, dass zur Minimierung der Gütefunktion Q das Gradienten-, das konjugierte Gradienten-, das Quasi-Newton-, das Newton-Raphson- oder das Newton-Verfahren verwendet wird.

6. Verfahren gemäss Anspruch 4, dadurch gekennzeichnet, dass zur Minimierung der Gütefunktion Q das Kuhn-Tucker-Verfahren verwendet wird.

7. Verfahren gemäss Anspruch 3, dadurch gekennzeichnet, dass die quadratische Gütefunktion gemäss Gleichung (12) durch die Betragsgütefunktion n-1 n * m - 1

(20) q = X* i » 2 - dealM(i*T4)|

j-0 J i-0

ersetzt wird.

8. Verfahren gemäss Anspruch 3, dadurch gekennzeichnet, dass die quadratische Gütefunktion gemäss Gleichung (12) durch eine der nachstehend gewichteten Gütefunktionen n-1 n * m - 1

(21) Q = Ï Xj*(4dvk. )'♦ 2 Hi*k(i*TA) - deak>1(i»TA)l!

j-0 J J i-0

n-1 n * in - 1

(22) Q = 2 x*|Advk J + 2 ^i*|ek(i*TÄ) - d0akM(i*TA)|

j-0 i-0

ersetzt wird, wobei bzw. m Gewichtsfaktoren darstellen, die vorgängig für jedes Zeitintervall j definiert werden müssen.

9. Verfahren gemäss Anspruch 3, dadurch gekennzeichnet, dass die direkte Berechnung der Sprungantwort nach Gleichung (8), zur Erreichung besser dämpfender Eigenschaften bei starken Messstörungen, durch ein Fehlerminimierungsverfahren ersetzt wird, wobei

11

5

10

15

20

25

30

35

40

45

50

55

CH 686 766 A5

a) eine durch Gleichung (23) definierte Streckenimpulsantwort gk(i*ÎA), d.h. die Reaktion der Strecke auf einen Impuls nach Gleichung (24), eingeführt wird:

1

(23) 9k(i*TA) = 2 vk(r*TA)*gk((i-r)*TA

r -0

(1 für i=0 1

(24) 5( i*TA ) = jo sonst ^

b) bei der Identifikation zur Begrenzung der Freiheitsgrade nur die ersten N-Werte berücksichtigt werden, so dass für gk(i*ÎA) gilt:

f/0 für 0<i<N-li

(25) 9k(i*TA) = \ 0 sonst j c) die bezüglich der Impulsantwort gk(i*ÎA) zu minimierende Gütefunktion F somit folgende Gestalt annimmt:

2

n * m - Ì

(26) F= 2

i - o

9ak(i*TA)-2 vk(r*TA)*gk((i-r)*TA)

o

♦ 'Y[gk(s*TA)-gk((s-l)*TA)]2

s - 0

und d) die Sprungantwort gemäss

(27) hk(i*TA) = Ì gk (r*TA )

r - 0

aus dem numerischen Integral der Impulsantwort hervorgeht.

10. Verfahren gemäss Anspruch 9, dadurch gekennzeichnet, dass die Bestimmung der Sprungantwort über ein Fehlerminimierungsverfahren im Frequenzbereich (z-Bereich) gefunden wird, wobei a) die Streckenübertragungsfunktion im z-Bereich gemäss Gleichung (28) gegeben ist

N

2 b_*z< -e(z) ...

(28) Gs(z) = = ,

V(z) N

l+2as *z(-3'

s - 0

e(z) und V(z) die z-Transformierten der zeitdiskreten Funktionen 8(ì*Ta) bzw. v(ì*Ta) sind und die Koeffizienten der Übertragungsfunktion as und br durch Minimierung des quadratischen Modellfehlers bestimmt werden;

(29) gk(j*TA) = r-i [Gs, k(z)]

und c) die Sprungantwort hiermit nach Gleichung (27) bestimmt wird.

11. Anwendung des Verfahrens nach einem der Ansprüche 1 bis 10 zur Herstellung von Strangpressprofilen mit kleinem und/oder wellenlängenabhängigem Emissionsgrad (e < 0,1) und/oder veränderlichen Oberflächencharakteristiken.

12. Anwendung des Verfahrens nach einem der Ansprüche 1 bis 10 zur Herstellung von Strangpressprofilen aus Aluminium oder Aluminiumlegierungen.

12