Informatik-Fachberichte 230 Herausgeber: W. Brauer im Auftrag der Gesellschaft fOr Informatik (GI)
Informatik-Fachberichte 230
Herausgeber: W. Brauer im Auftrag der Gesellschaft fOr Informatik (GI)
A. Bode R. Dierstein M. Gobel A. Jaeschke (Hrsg.)
Visualisierung von Umweltdaten in Supercomputersystemen 1. Fachtagung Karlsruhe, 8. November 1989
Proceedings
Springer-Verlag Berlin Heidelberg New York London Paris Tokyo Hong Kong Barcelona
Herausgeber
A. Bode Institut fur Informatik der Technischen Universitat Munchen Barerstr. 23, 0-8000 Munchen 2
R. Oierstein Deutsche Forschungsanstalt fur Luft- und Raumfahrt e. V. Zentrale Oatenverarbeitu ng 0-8031 Oberpfaffenhofen
M. Gobel Fraunhofer-Arbeitsgruppe Graphische Oatenverarbeitung Wilhelminenstr. 7, 0-6100 Darmstadt
A.Jaeschke Kernforschungszentrum Karlsruhe GmbH Institut fur Oatenverarbeitung in der Technik Postfach 3640, 0-7500 Karlsruhe 1
CR Subject Classifications (1987): B.2.1, B.4.4, C.1.1-2, C.5.1, 1.3, J.2, J.6., J.m
ISBN-13: 978-3-540-52746-6 e-ISBN-13: 978-3-642-75805-8 DOl: 10.1007/978-3-642-75805-8
CIP-Titelaufnahme der Deutschen Bibliothek. Visualisierung von Umweltdaten in Supercomputersystemen: ... Fachtagung ... ; proceedings. - Berlin; Heidelberg; New York; London; Paris; Tokyo: Springer, 1. Karlsruhe, 8. November 1989. - 1990
(lnformatik-Fachberichte: 230)
NE:GT
Dieses Werk ist urheberrechtlich geschOtzt. Die dadurch begrOndeten Rechte, insbesondere die der Obersetzung, des Nachdrucks, des Vortrags, der Entnahme von Abbildungen und Ta-bell en, der Funksendung, der Mikrovertilmung oder der Vervielfaltigung auf anderen Wegen und der Speicherung in Datenverarbeitungsanlagen, bleiben, bei auch nur auszugsweiser Verwer-tung, vorbehalten. Eine Vervieifaltigung dieses Werkes oder von Teilen dieses Werkes ist auch im Einzelfall nur in den Grenzen der gesetzlichen Bestimmungen des Urheberrechtsgesetzes der Bundesrepublik Deutschland vom 9. September 1965 in der Fassung vom 24.Juni 1985 zulassig. Sie ist grundsatzlich vergOtungspflichtig. Zuwiderhandlungen unterliegen den Straf-bestimmungen des Urheberrechtsgesetzes.
© Springer-Verlag Berlin Heidelberg 1990
2145/3140-543210 - Gedruckt auf saurefreiem Papier
Vorwort
Die Datenmengen, die bei der Modellierung, der Simulation und der Oberwaehung von Umweltsituationen an fallen, sind heute schon so groB, daB sic nur noch mit Supercomputersystemen in verantwortbaren Zeiten sinnvoll bearbeitet werden konnen. Was fur die Eingabedaten gilt, wiederholt sich bei der Ausgabe: Die Ergebnisse der Berechnungen oder Auswertungen sind in den meisten Fallen schier endlos scheinende Zahlenkolonnen, die, so hergenommen, wie sie entstanden sind, keine brauchbaren Aussagen gestatten. Die Ergehnisdaten milssen deshalb weiter analysiert und interpretiert werden, sei es durch erfahrenes und DV-kundiges Fachpersonal, sei es durch weitergehende Transformationen und geeignete Formen der Darstellung.
Die Visualisierung von Umweltdaten setzt in dies em Punkt ein. Sie versucht, das umfangreiche Datenmaterial mit Werkzeugen der Graphischen Datenverarbeitung in Form von Bildern darzustellen und dabei die Modellvorstellungen derjenigen Wissenschaftler und Ingenieure wiederzugeben, 'die in erster Linie Fachleute ihrer wissenschafUichen und technischen Disziplinen sind und erst danach - wenn ilberhaupt - DV-Spezialisten.
Darilber hinaus ist es Aufgabe und Ziel der Visualisicrung umfangreicher und komplcxcr Mcngen wisscnschaltlichcr Daten, die Modcllicrung und Stcucrung von UrnwcItprozcsscll nicht mehr nur rein numerisch anzugehen, sich also ganz auf die Datenstrukturen Ztl heziehen, sondern sie graphiseh-interaktiv auf der Grundlage bildlicher Darstellungen aufzubauen und durchzufuhren.
Del' Fachbereich 4 Irljormationstec/znik und Technische Nu/zul1g der In!orl1lafik der GeseIlschaft fur Informatik hat deshalb am 8. November 1989 im Kernforschungszentrum Karlsruhe (KfK) die 1. Fachtagung zum Thema
Visualisierung von Umweltdaten in Supercomputersystemell: Datenstrukturell - Modelle - Graphische Darstellung
veranstaltet, und zwar im Rahmen des 4. Symposiums bljormatik im Umweltschutz der Gesellschaft fur Informatik (GI) und der 10. Jahrestagung del' Gesellschqftfilr bljormatik ill dcr Landwirtschaft (GIL).
Schwerpunkte der Tagung waren die
• Visualisierung raumlicher Daten, • Anwendung von I Iochlcistungssystemen, • Visualisicrung und Interaktion in Umweltinformationssystemcn.
Viele haben dazu beigetragen, daB dicse l'achtagung durchgefuhrt werden konnte. Dank gilt den I nitiatoren und Organisatoren dieser Tagung, den IIerren Prof. Dr. Arndt Bode, Ttl Milnchen, Prof. Dr. Jose Encarnac;:ao, Til Darmstadt, Dr. Andreas Jaeschke, KfK Karlsruhe und Dipl.-Inform. Martin Gobel, FhG-AGD, Darmstadt, Dank auch allen Mitgliedern des Organisations- und des Programmkomitees - Dank vor allem aber den Vortragenden, die durch ihre Beitrage diese Tagung und diesen Tagungsband moglich machten.
OberpfafTenhofen, im Januar 1990 Rildigcr-Dicrstein
Inhaltsverzeichnis
Farbtafeln
Visualisierung diumlicher Daten
M. Fruhauf, K. Karlsson. FhG-AGD. Darmstadt Visualisierung von Volumendaten in verteilten Systemen
W. Kruger. ART & COM. Berlin Visualisierung 3-dimensionaler skalarer Datenfelder - Transporttheorie-M odell ........ 11
H ochleistungssysteme
H. Burfeindt. VAG Wolfsburg; H. Zimmer. SFE Berlin Schalldruckberechnung im PKW-Innenraum .............................................................. 23
W. Giloi. GMD-FIRST. Berlin Der zukunftige Parallelrechner - von der Superworkstation zum Supercomputer ...... 51
J. Hillmann. VAG W o(f~burg Einsatz von Supercomputern flir Crashberechnung ................. .................... ......... ...... 63
Visualisierung und Interaktion in Umwelt-Informationssystemen
R. Muller-Meernach. F.J. Prester. S.E.P.P. GmbH. Rotlenbach P. Weigert. Bundesgesundheitsamt ZEBS. Berlin
Kartographisch orientierte Darstellung von Mef3werten ............................................ 81
R. Riegger. EDV Studio Ploenzke. Milnchen Die verschiedenen Ebenen bei der Bearbeitung von Umweltdaten mit Geo-informationssystemen ................................................................................................... 91
F. Schmidt (IKE). R. LWkemeyer (BMU). J. Lieser (BGA-ZdB). A. Bohm. G.llehn. Th. Milller. W. Scheuermann. V. Sundararaman, IKE Universitiit Stuttgart
Ver<;trbeitung von Umweltdaten unter Real-Time-Bedingungen - Konzept und prototypische Realisierung ............................................................................... ........... 103
Farbtafeln
Farbige Abbildungen zu den Beitdigen
H. Burfeindt Schalldruckberechnung im PKW-Innenraum
Bild 5 und 10
J. Hillmann Einsatz von Supercomputern fur die Crashberechnung
Bild 11
P. Riegger Die verschiedenen Ebenen bei der Bearbeitung '"
Bild 8 - Bild 11
IX
1.e.C£-0S SYMM. MODE DES GOLF-INNENRAUMES
l. 3SE-0S
I. 01£-05
6.7e.E- 0e.
'I.
FREQUENZ 71 HZ
BILD 5
? --------------------------------------------~
5 . 00E-0l
0.
-5.
- 1.~+10
DISlU~rNT FI!1NG€S i---
BILD 10
'--- BETELIGUNG DER KAR088ERE-FLAECHEN AM 8CHALLDRUCK
-
1-O
YNA
301
t r --
-PA
MC
RA
SH I
. II
I
• ~
t{j \
tI
, , ~
~/,
~ ~
-.. ---
I--
r,t'''
~~ ~
--'
~
&
tI I
I
I I-
-
,VI
I '-
... ~~
~ L.
. J
.. , 1"'';'
I-
--.
--1
----
.. -f-
. ~
~
Zei1
Bild
11
: V
erg
leic
h d
er
Au
fpra
llkra
ft-Z
eit
Ku
rve
n
fUr
DY
NA
3D u
nd
PA
MC
RA
SH
> C -en n ~
x ji"
CD
'" ., ED - -
Bild 8.
Bild 9.
Altlastenkataster der Stadt Wuppertal
UMW(l1 A Tl.AS S(RLIN
XI
Bewertung der Belastungspotentiale von Altlasten
OIo,-OIo;ISC,," PJo""'"'9~s\rufM''''1 io1l.'UIClllIftOlt I Um.~lt
A.:tJo:lten\lfldoc:hlsrlOCMn
'*---'" .... .. . ................ -.
_ ......... ...... ~~
_ ..... _w ===.:.-~-
t'lt-f~r .. "'9 1f1"'I~-=<,* "''''t..~1 4I1i01.t~ot»(fI\J1«
XII
Bild 10. Landnutzungsklassifizierung im Raum Rosenheim mit LANDSAT-TM-Daten
Bild 11. Vegetationsverteilung der Stadt Munchen - Klassifizierung mit LANDSAT-TM- und SOJUSDaten
Visualisierung von Volumendaten in verteilten Systemen
Zusammenfassung
Martin FrUhauf, Kennet Karlsson
Fraunhofer-Arbeitsgruppe Graphische Datenverarbeitung
Wilhelminenstr. 7, D-6100 Darmstadt
Durch den Einsatz leistungsfahiger Rechner wird es in zunehmendem Maile maglich, auch die in
nere Struktur komplexer, heterogener Objekte zu visualisieren. Das erfordert die Verarbeitung von
umfangreichen Volumenmodellen. Es werden an Beispielen Maglichkeiten zur Generierung solcher
Volumenmodelle gezeigt. Danach werden die Verarbeitungsschritte zur Generierung von 3D-Dar
stellungen der Volumenmodelle und Maglichkeiten zu deren Beschleunigung durch schnelle Algo
rithmen und Datenreduktion vorgestellt. Zum SchluB wird auf die Moglichkeiten der graphisch-in
teraktiven Visualisierung solcher Daten in verteilten Systemen bestehend aus Supercomputern und
leistungsfahigen Graphik-Workstations eingegangen.
1. Einlei tung
Leistungsfahige Supercomputer ermaglichen die Lasung von immer komplexeren Problemen. Dabei
werden groBe Datenmengen verarbeitet oder als Berechnungsergebnisse produziert. Diese Daten
mengen werden in Zukunft noch weiter zunehmen. Man kann sagen, daB Forschungsergebnisse da
von abhangen, wieviele Daten und Simulationsergebnisse der Forscher begreifen und verarbeiten
kann [1]. Daher muB man geeignete Wege suchen, diese Informationen zu erfassen, zu analysieren
und sie in einer zweckmaBigen Form an Wissenschaftler und Ingenieure weiterzugeben. Diese
Technik wird Visualisierung genannt [2]. Die Computergraphik spielt hierbei eine zentrale Rolle. In
zunehmenden MaBe wird ein Bedarf an graphisch-interaktiver Visualisierung bestehen, d. h. der
Forscher muB die Maglichkeit haben, die Darstellungsweise der Berechnugsergebnisse interaktiv zu
steuern, zu manipulieren und seinen speziellen Bedtlrfnissen anzupassen. Interaktive Anwendungen
stellen aber besonders hohe Anforderungen an die Zeit, die die Generierung von Darstellungen
benotigen darf. Diese Zeitanforderungen einzuhalten, erfordert den Einsatz moglichst leisi-ungsfahi
ger Rechner auch fUr die graphischen Verfahren zur Visualisierung von MeBwerten oder Simula
tionsergebnissen. Andererseits bieten Graphik-Workstations exzellente Moglichkeiten zur qualitativ
2
hochwertigen Darstellung, so dafi der Einsatz emes verteilten Systems bestehend aus emem
Supercomputer und einer oder mehrerer Graphik-Workstations vorteilhaft erscheint.
Viele Ansatze zur Visualisierung beschranken sich auf die Verarbeitung oberflachenorientierter Da
tenstrukturen. Einhiillende Oberflachen von Mefipunkten im dreidimensionalen Raum werden da
bei durch Triangulierung erzeugt, d. h. die darzustellende Oberflache besteht aus einer - meistens
relativ geringen - Anzahl von dreieckigen Polygonen. In zunehmendem Mafie wird aber auch die
Visualisierung komplexer Volumen moglich und notwendig. Dabei interessiert den Wissenschaftler
nicht nur die Oberflache eines Objektes, Korpers oder allgemein eines Volumens, sondern auch des
sen inn ere Struktur. Es ist daher notwendig stets das gesamte Volumenmodell verfiigbar zu haben.
Zur Visualisierung werden dann Schnitte und Transparenzen im Volumen vom Benutzer definiert
und die entsprechenden Darstellungen erzeugt. Ein Volumenmodell besteht aus einer Anzahl von
einzelnen, regelmafiig angeordneten Volumenelementen. Diese Volumenelemente sind mit be
stimmten Wert en belegt. Ein Volumenelement in diesem Modell wird Voxel genannt. Beispiels
weise besteht ein Volumen aus 2563 Voxeln.
Die im folgenden vorgestellten Verfahren und 'Oberlegungen zur Visualisierung von Volumendaten
sind in den verschiedensten Anwendungsgebieten der Visualisierung technischer und naturwissen
schaftlicher Daten und Mefiwerte anwendbar. Insbesondere auch auf Umweltdaten aus der Mete
orologie oder Geophysik sowie bei der Visualisierung der Ergebnisse von Stromungsberechnungen.
2. Volumenmodelle
2.1 Generierung von Volumenmodellen aus Mefiwerten
Oft sind Mefipunkte in emem regelmafiigen Gitter im dreidimensionalen Raum angeordnet, Sle
konnen aber auch unregelmafiig im Raum verteilt sein. Die Daten an den Mefipunkten konnen
Skalare, Vektoren oder Funktionen mit ihren Parametern sein.
Bei der Visualisierung der Ergebnisse von Stromungsberechnungen besteht das Volumenmodell aus
dem berechneten Druck sowie dem Geschwindigkeitsvektor an den Knotenpunkten des Finite-Ele
ment-Gitters.
Bei medizinischen Anwendungen wird ein Volumenmodell aus einer Sene von Schichtbildern durch
Interpolation zwischen den Schichten erstellt. Diese Schichtbilder konnen z. B. Computer- oder
Kernspin-Tomogramme sein. Die Daten an den Mefipunkten sind Dichtewerte der untersuchten
Gewebe und sind skalar. In einem Volumenmodell ist stets die gesamte Information zur darzustel
lend en Szene enthalten. Dadurch lassen sich graphisch-interaktive Anderungen des Modells, wie sie
bei der Operationsplanung benotigt werden, ohne eine aufwendige Neuberechnung der Ob
jektoberflachen durchfiihren.
Auch die Visualisierung von Materialeigenschaften erfordert die Verarbeitung von Volumen
modellen. Dabei kann das Volumenmodell aus den geometrischen Abmessungen des Werkstiicks
3
generiert werden. Verschiedenen Materialeigenschaften oder -zustanden - z. B. Temperatur, Span
nung usw. - konnen dann durch verschiedene Werte der Voxel reprasentiert werden.
Ein Problem bei der Verwendung von Volumenmode11en stellt die grofle Anzahl der in einer Szene
enthaltenen Voxel dar. Bei der a11gemein iiblichen Auflosung von Computer-Tomogrammen in der
Medizin von 256 x 256 Pixeln ergibt sich beispielsweise ein V olumen von 256 x 256 x 256 ~ 16
Millionen Voxeln. Hier wird deutlich das solche Datenmengen nur von Rechnern mit entsprechen
der Hauptspeicherkapazitat effizient verarbeitet werden konnen. Die Datenmenge laflt sich zwar
durch die Erste11ung eines Octrees [3] komprimieren, ailerdings kann die Rotation des im Octtree
gespeicherten Volumen nicht mehr effizient ausgefiihrt werden [4].
2.2 Generierung von 3D-Darste11ungen aus Volumenmode11en
Zur Generierung verschiedener Ansichten der als Volumenmode11 vorliegenden Szene sind folgende
Operationen durchzufiihren:
1. A uswahl der darzuste11enden V oxel
2. Schneiden des Volumens
3. Rotation des Volumens
4. Projektion des Volumens
5. Berechnung der Oberflachennormale
6. Schattieren der Darste11ung
7. Skalierung der Darste11ung
Die Operationen 1 - 4 miissen immer im Objektraum durchgefiihrt werden, d. h. sie erfordern die
Verarbeitung ailer Volumenelemente. Einfache Verfahren berechnen die Oberflachennormale im
Bildraum mit Hilfe eines Z-Buffers. Anspruchsvo11ere Darste11ungen erfordern die Berechnung der
Oberflachennormalen im Objektraum. Die Oberflachennormale wird dabei aus den vorliegenden
Meflpunkten hergeleitet. Da die Oberflachennormale zusammen mit dem Farbwert an einem Mefl
punkt die Darste11ung des Meflwertes bestimmt, sind je nach der Art der vorliegenden Daten und
der gewiinschten Darste11ung unterschiedliche Verfahren zur Berechnung notwendig, auf die hier
nicht im einzelnen eingegangen werden so11.
Wahrend die Operationen 1 und 2 durch einfache Vergleiche ausgefiihrt werden konnen, erfordert
die Rotation der Szene um beliebige Winkel die Anwendung einer 3 x 3 Rotationsmatrix auf jedes
einzelne Volumenelement. Diese Operation ben5tigt einen groflen Anteil der zur Darste11ungsbe
rechnung notwendigen Rechenzeit. Daher so11 hier ein Verfahren vorgeste11t werden, mit dem der
Rechenaufwand fiir die Anwendung der Rotationsmatrix erheblich verringern laflt, in dem man sich
die regelmaflige Anordnung der Voxel im Volumen zu Nutze macht. Dennoch ist es sinnvo11, so
viele Schritte zur Visualisierung von Volumenmode11en wie moglich auf Supercomputern durch
zufiihren, um die Rechenzeit zu verkiirzen.
4
2.3 Rotation und Projektion von Volumenmodellen
Die Koordinaten der Voxel werden als ganzzahlig angenommen. Der Bildraum wird back-to-front
(BTF) traversiert und mit Hilfe der Inversen der Rotationsmatrix werden die Positionen der im
Bildraum einzutragenden Voxel im Objektraum bestimmt. Der Bildraum wird durch ein zwei
dimension ales Bildraster und einen z-Buffer, in dem die z-Koordinaten der Voxel im Bildraster ab
gelegt sind, reprasentiert. Die Projektion der Szene geschieht durch die Back-to-Front-Traversie
rung des Volumens [5) und Eintragen des rotierten Voxels in das Bildraster bei gleichzeitigem Ein
trag der z-Koordinate in den z-Buffer.
z
Abb. 1: Ergebnis der BTF-Traversierung
Soil die Projektion des V olumens nicht gleichzeitig mit der Rotation erfolgen, da noch weitere Be
rechnungsschritte im Objektraum durchgefiihrt werden miissen, kann mit dies em inkrementellen
Rotationsverfahren natiirlich auch ein rotiertes Volumenmodell erzeugt werden.
Durch die Anwendung des BTF-Verfahrens wird eine implizite hidden-surface Eliminierung durch
gefiihrt. Falls der Objektraum traversiert wird und die Voxel in den Bildraum rotiert und proje
ziert werden, k6nnen durch Rundung auf ganzzahlige Positionen Locher im Bildraum entstehen.
Bei der Traversierung des Bildraumes wird das vermieden, da fiir jede Position im Bildraum genau
ein Voxel aus dem Objektraum projeziert wird.
Sei k die Position eines Voxels im Objektraum, k' die Position des Voxels im Bildraum ]lnd R die
Rotationsmatrix.
Es gilt:
k' = k * R
mit T = R-1 folgt
k = k' * T
oder
Sei
[kx', ky', kz'J = [(k-1)x'+ I, (k-1)y', (k-1)z'J
dann gilt
daraus folgt
5
(2.1)
(2.2)
(2.3)
(2.4)
(2.5)
(2.6)
Dies gilt analog auch fUr die Inkrementierung der y- oder z-Koordinate SOWle fUr die Dekre
mentierung.
In diesem Algorithmus enthalten nur die Koeffizienten der Rotationsmatrix gebrochene Zahlen.
Wenn fUr diese Koeffizienten eine Festpunktdarstellung gewahlt wird, indem man sie mit 2n multi
pliziert und danach in eine ganze Zahl konvertiert, ben5tigt man zur Rotation eines Voxels drei
Ganzzahladditionen und drei Shift-Operationen urn n Stellen. Mit den Shift-Operationen wird die
Multiplikation der Matrixkoeffizienten egalisiert.
Ein anderes Rotationsverfahren benutzt Tabellen zur Rotation der Voxel [5J. In dies en Tabellen
sind Eintrage tij * k mit k = 0, I, 2, ... , 255 und i, j = I, 2, 3 abgelegt. Die Position eines Voxels
im Bildraum kann dann durch Addition der entsprechenden Tabelleneintrage berechnet werden.
Dieses Verfahren hat Vorteile, wenn eine gro£e Zahl der Voxel im Objektraum fUr die Darstellung
nicht relevant sind, so dan sie nicht in den Bildraum transformiert werden mUssen. Bei dem oben
vorgestellten inkrementellen Verfahren bringt das keine Vorteile, da jede Transformation von der
Transformation des Vorgangers abhangt.
Mit Hilfe des z-Buffers kann zu jedem Voxel im Bildraster der Gradient berechnet werden [6J. Der
Gradient ist ein MaE fUr die Neigung der durch die Voxel approximierten Objektoberflache. Er ist
die Grundlage fUr die Generierung schattierter Darstellungen mit Beleuchtungsmodellen [2L :Qi~ses Schattieren im Bildraum hat den Vorteil, dan das Berechnen der Neigungen nur fUr die Voxel
durchgefuhrt werden mu£, die auch tatsachlich dargestellt werden. Dies ist sinnvoll, denn die
6
einmalige Berechnung der Neigungen im Objektraum ist nicht moglich, wenn interaktiv Anderun
gen des Volumenmodells vorgenommen werden. Allerdings hangt die Qualitat der Darstellungen
entscheidend von der Prazision der Neigungsberechnung abo Die Berechnung der Neigung im Bild
raum ist daher nicht immer ausreichend. Auf die Verfahren zum Rendering von Voxelmengen soIl
hier nicht eingegangen werden. Es sei verwiesen auf [6 8 9 10]. Auch Teilberechnungen dieser Ver
fahren konnen auf einem Supercomputer durchgefilhrt werden.
2.4 Generierung eines oberflachenorientierten Datenmodells aus dem Volumenmodell
Ein Verfahren zur Reduktion der Datenmenge und der damit verbundenen Bearbeitungszeit, das
sogenannte Cuberille-Verfahren, ist von G. T. Herman [11] entwickelt worden. Ausgehend von ei
nem Oberflachenelement eines Objekts wird die ganze Oberflache des Objekts mit Hilfe eines Gra
phensuchverfahrens im segmentierten Volumen verfolgt. Ausgabe des Algorithmus ist eine Menge
von Oberflachenelementen, die die Oberflache des Objekts bilden. Ein Oberflachenelement ist in
diesem Modell eine Voxelseite, d.h. eine der sechs Seiten eines wilrfelformigen Volumenelements.
Das Datenmodell baut auf sauberen mathematischen Definitionen auf und hat u.a. den Vorteil, daJl
mehrere Objekte in einem Volumen unabhangig voneinander behandelt werden konnen.
Filr die Darstellung der Cuberille-Oberflache konnen dieselben Verfahren verwendet werden wie fiir
die direkte Darstellung des segmentierten Volumens, z.B. das Back-to-Front Verfahren, wobei die
relativ geringe Anzahl von zu zeichnenden Voxeln die Darstellung erheblich schneller macht.
Ein groller Nachteil dieses Modells und oberflachenorientierter Datenmodelle im allgemeinen, ist,
dall interaktive Manipulationen aufwendig werden, da die Information ilber die innere Struktur der
Objekte verlorengegangen ist. Jede Manipulation mull auf dem segmentierten Volumen ausgefiihrt
werden und die Oberflache mull dann neu generiert werden. Oberflachenorientierte Modelle werden
deshalb verwendet, urn schnelle Ansichten aus verschiedenen Richtungen zu erzeugen, werden aber
fiir interaktive Manipulationen kaum eingesetzt. Nach un serer Auffassung ist es aber durchaus
moglich, ein lokales "Stopfen" eines Teiles der Oberflache nach einer lokalen Manipulation in dem
Cuberille-Modell schnell durchzufiihren.
Ein Objekt 0 in dem segmentierten Volumen besteht aus einer Menge von Voxeln. Die Oberflache
des Objekts 0 besteht aus einer Menge F von Oberflachenelementen (Voxelseiten). Betrachten wir
ein Voxel V in dem Objekt O. Wenn das Voxel V geloscht wird (in dem segmentierten Volumen)
miissen die Seiten von V, die an der Oberflache lagen, aus der Menge F geloscht werden und neue
Voxelseiten, die jetzt neu zur Oberflache dazugekommen sind, miissen zu F hinzugefiigt werden.
Diese neuen Oberflachenelemente sind identisch mit den Seiten von V, die nicht an der Oberflache
lagen. Bei dem Loschen eines einzelnen Voxels sind damit die notwendigen Anderungen in der
Menge der Oberflachenelemente unmittelbar bestimmt.
Ein allgemeiner Ansatz fur die lokale Neugenerierung der Oberflache ist der lokale Einsatz des
Graphensuchalgorithmus. Betrachtet man die Oberflachenelemente als Knoten und die Rander der
7
Oberflachenelemente als Kanten, kann die Oberflache in emem gerichteten Graphen dargestellt
werden. In dem Graphen hat jeder Knoten genau zwei Vorganger und genau zwei Nachfolger. Be
trachten wir wieder das Voxel V im Objekt O. Wenn V im segmentierten Volumen geloscht wird,
mussen die Seiten von V, die an der Oberflache lagen, aus der Menge F geloscht werden. Der AI
gorithmus wird dann am Vorganger der geloschten Oberflachenelemente fortgesetzt. Wenn der AI
gorithmus terminiert, ist die Oberflache wieder vollstandig (Beweis dazu, siehe [11]). Mit diesem
Verfahren kann die Oberflache auch bei groBeren Manipulationen lokal neugeneriert werden. Der
gerichtete Graph muB dabei nicht gespeichert werden, da die Vorganger eines Oberflachenelements
sich sehr schnell berechnen lassen.
Das Hauptproblem bei der Implementierung des Graphensuchalgorithmus ist die Handhabung der
groBen Datenmengen. Das Volumen besteht beispielsweise aus 256 x 256 x 256 Voxeln. Wird jedes
Voxel mit vier Bit reprasentiert belegt das Volumen 8 MByte. Die Ausgabe des Algorithmus kann
in einem vier-dimensionalen Array mit 256 x 256 x 256 x 6 Voxelseiten gehalten werden, wobei
jede Voxelseite mit einem Bit reprasentiert wird. Dieses Bit gibt die Zugehorigkeit zur Oberflache
des Objekts an. Mit 12 MByte ist dieses Array viel grBBer als die in [11] beschriebene Datenstruk
tur, die fur eine Implementierung auf einem Minicomputer entwickelt wurde. Mit dem Array wird
aber der Algorithmus einfacher und schneller. Die Rechenzeit verhalt sich linear zu der Anzahl von
Oberflachenelementen. Die groBen Datenmengen set zen allerdings einen groBen Hauptspeicher vor
aus und das Verfahren eignet sich deshalb fur die Verarbeitung auf einem Supercomputer im
Hintergrund der Visualisierung.
3. Verteilte Systeme
Die oben erwahnten Schritte zur Bildberechnung 1 - 5 sind unabhangig von dem zur Darstellung
verwendeten Graphikpaket und der verwendeten Graphik-Workstation. Sie dienen zur Erzeugung
eines Bildraster sowie des Gradienten zu jedem Bildpunkt als MaB fur die Neigung der Objekt
oberflache. Diese Daten dienen bei der Berechnung auf einem Supercomputer als Eingabe zur Bild
generierung auf einer Graphik-Workstation.
Da die einzelnen Voxel voneinander unabhangig sind und unabhangig voneinander dargestellt wer
den, laBt sich die Bildberechnung aus Volumenmodellen auch leicht auf Parallelrechnern ausfuhren,
ohne daB dazu spezielle Algorithmen entwickelt werden mussen. Beim oben vorgestellten Algorith
mus zur Rotation der Voxelraume kann man zunachst die Rotation fur das jeweils erste Voxel ei
ner Schicht berechnen, und kann dann die Rotation der Schichten parallel ausfuhren. Wenn die
Projektion parallel durchgefuhrt wird, ist zu beachten, daB das Eliminieren verdeckter Voxel nicht
mehr durch einen einfachen Back-to-Front-Algorithmus durchgefuhrt werden kann. Allerdings
kann der BTF-Algorithmus in zwei Schritten durchgefuhrt werden. 1m ersten Schritt wird der
BTF-Algorithmus parallel auf bestimmte disjunkte Intervalle der Tiefe des Voxelraumes angewen
det. Auf die Ergebnisse dieser Berechnungen wird dann erneut ein BTF-Algorithmus angewendet
und das endgultige Bildraster erzeugt.
8
3.1 Kommunikation in verteilten Systemen
UNIX-Systeme bieten eine komfortable M<lglichkeit fiir die Kommunikation zwischen Prozessen.
Zwei Prozesse sind iiber "sockets" gekoppelt und kommunizieren iiber diese miteinander. Diese
Prozesse k<lnnen auch auf verschiedenen Rechnern ablaufen. Wenn die Verbindung zwischen zwei
Prozessen einmal eingerichtet wurde, k<lnnen beliebig grofie Datenpakete von einem Prozefi zu ei
nem anderen transferiert werden. Damit k<lnnen aufwendige Berechnungen zur Bilderzeugung auf
einem schnellen Hintergrundrechner durchgefiihrt werden und die Graphik-Workstation zur glei
chen Zeit ausschliefilich zur Generieqmg von Bildern verwendet werden. Voraussetzung dafiir ist,
dafi Formate zum Austausch der Bilddaten festgelegt werden. Solche Formate k<lnnen Bildraster
mit Z-Buffer oder den Gradient en zu jedem Bildpunkt wie oben beschrieben oder bereits vollstan
dig berechnete Rasterbilder sein [13].
Bei der Verteilll.ng der einzelnen Berechnungsoperationen ist der Zeitgewinn durch die Berechnung
auf einem schnellen Hintergrundrechner gegen den Zeitaufwand fiir den Datentransfer abzuwagen.
Ein Bildraster belegt z. B. 256 x 256 Byte (bei einem Byte pro Bildpunkt). Die gleiche Daten
menge wird noch einmal fiir den z-Buffer ben<ltigt. Diese Datenmengen k<lnnen nur auf einem LAN
effizient iibertragen werden, d. h. eine verteilte Berechnung ist nur sinnvoll, wenn z. B. UNIX-Sys
teme iiber Ethernet gekoppelt sind.
Es ist natiirlich auch m<lglich aile Verarbeitungsschritte zur 3D-Darstellung von Volumenmodellen
auf einem Supercomputer durchzufiihren und die dort berechneten Rasterbilder auf eine Graphik
Workstation zu iibertragen. Allerdings fiihren leistungsfahige Graphik-Workstations Schattierungen
heute bereits mit spezieller Hardware durch, so dafi die dort gebotene Leistung durch Software auf
Supercomputern i. a. nicht erreicht wird. Weiterhin werden in Zukunft Schattierungsfunktionen in
Graphikstandards [14] auf Workstations verfiigbar sein, die dann dort genutzt werden k<lnnen.
3.2 Interaktion in verteilten Systemen
Fails eine Graphik-Workstation und ein schneller Hintergrundrechner iiber ein schnelles Netz ge
koppelt sind, lassen sich in diesem System auch graphisch-interaktive Anwendungen realisieren.
Mit einer graphischen Benutzungsoberflache werden yom Benutzer Parameter zur Darstellung der
Szene spezifiziert. Parameter zur Darstellung sind z. B. Winkel zur Rotation der Szene oder
Schnittebenen. Sie parametrisieren die Rechenschritte zur Darstellung, die auf dem Hintergrund
rechner ablaufen. Die berechnete Darstellung wird dann wie oben beschrieben auf die Graphik
Workstation iibertragen. Nach unseren Erfahrungen mit Systemen, in denen Benutzungsoberflache
und Anwendungsprogramm auf unterschiedlichen Rechnern installiert sind und ablaufen, stellt die
Ubertragung der spezifizierten Parameter iiber ein Ethernet keine yom Benutzer wahrnehmbare
Verz<lgerung bei der Ausfiihrung des Anwendungsprogramms dar.
9
4. Ausblick
Es wurde ein Uberblick fiber die bei der interaktiven Visualisierung von komplexen, heterogenen
Volumenmodellen auftretenden Probleme gegeben, und es wurden Moglichkeiten der Be
schleunigung der Visualisierung vorgestellt. Die damit erreichten Antwortzeiten konnen interaktive
Anwendungen aber noch keineswegs befriedigen, sind aber in Anbetracht der Komplexitat des
Visualisierungsproblems beachtlich. Wie in allen Anwendungsbereichen der Computergraphik wird
sich in Zukunft ein weiterer Zeitgewinn durch schnellere Hardware erreichen lassen. Ferner ist es
moglich, Darstellungen minderer Qualitat - etwa als Preview - in kurzeren Zeiten zu erzeugen.
Literatur
[1]
[2]
[3]
[4]
[5]
[6]
[7]
[8]
[9]
[10]
[11]
[12]
Fangmeier, S. M.: "The Scientific Visualization Process" in: Proc. ACM Siggraph '88, Vol. 20, pp. 26, Atlanta, USA, 1988.
Encarnacao, J. L.; Schonhut, J.: "High Performance, Visualisation and Integration: The Computer Graphics Headlines for the 90's" Proc. IFIP TC 5 Conference on CAD/CAM Technology Transfer, Mexico City, Mexico, North-Holland, 1988.
Meagher, D.: "Geometric Modeling Using Octree Encoding" Computer Graphics and Image Processing 19, pp. 129, 1982.
Samet, H.; Webber, R. E.: "Hierachical Data Structures and Algorithms for Computer Graphics" IEEE GC&A, pp. 48, May 1988.
Frieder, G.; Gordon, D.; Reynolds, R. A.: "Back-to-Front Display of Voxel-Based Objects" IEEE CG&A, pp. 52, Jan 1985.
L. Chen, G. T. Herman, R. A. Reynolds, J. K. Udupa: "Surface Shading in the Cuberille Environment" IEEE CG&A, Dec 1985.
Foley, J. D.; van Dam, A.: "Fundamentals of Interactive Computer Graphics" Addison Wesley, 1984.
Kajiya, J. T.; "Von Herzen, B. P.: Ray Tracing Volume Densities" Computer Graphics, pp.165-174, VoLl8, No.3 (1984)
Sabella, P.: "Rendering Algorithm for Visualizing 3D Scalar Fields" Computer Graphics, pp.51-58, Vol.22, No.4 (1988)
Upson, C.; Keeler, M.: "V-BUFFER: Visible Volume Rendering" Computer Graphics, pp.59-64, Vol.22, No.4 (1988)
Artzy, E.; Frieder, G.; Herman, G. T.: "The Theory, Design, Implementation and Evaluation of a Three-Dimensional Surface Detection Algorithm" Computer Graphics and Image Processing 15, pp.1-24, (1981)
Boecker, F. R. P.; Hoehne, K.-H.; Tiede, U.; Riemer, M.: "Algorithmen und Datenstrukturen fur die Interaktion in dreidimensionalen medizinischen Szenen" in: Hommel, G.; Schindler, S. (Hrsg.): Proc. GI - 16. Jahrestagung, S. 659-675, SpringerVerlag, Berlin (1986)
10
[13] Hofmann, G. R.: "FTCRP - Ubertragungsformat fiir den farbtreuen, gerateunabhangigen Transfer von Rasterbildern" Fb,G-AGD, (1988)
[14] ISO: "Programmer's Hierachical Interactive Graphics System - Basic Rendering (PHIGS BR)" ISO/WD.
VISUALISIERUNG 3-DIMENSIONALER SKALARER DATEN FELDER:
TRANSPORTTHEORIE-MODELL
Zusammenfassung
Wolfgang Kriiger ART +COM e. V.
Hardenbergplatz 2 D-1000 Berlin 12
In dieser Arbeit wird eine allgemene Methode fur die Visualisierung von Volumendaten aus Wissenschaft und Technik vorgestellt. Es werden insbesondere die Modellgrundlagen, die moglichen Anwendungen und die Relationen zu vorhandenen Modellen diskutiert. Das verwendete Modell, die lineare Transporttheorie zur Beschreibung der Ausbreitung von virtuellen Lichtstrahlen, enthdlt unterschiedlichste Moglichkeiten, 3-dimensionale Datenfelder mit inneren Strukturen (z. B. punktformigen Details oder starken Inhomogenitiiten auf Fliichen) zu visualisieren.
1. Einfuhrung
Zu den klassischen Bereichen der Wissenschaft, Theorie und Experiment, kam mit der Entwicklung leistungsfiihiger Computer ein dritter hinzu, der der rechnergestutzten Simulation. Die Fiihigkeit der Supercomputer bzw. moderner Experimentaleinrichtungen, gewaltige Datenmengen in kurzer Zeit zu berechnen bzw. durch Messungen zu erzeugen, fuhrt zu einer nicht mehr auswertbaren Flut von "zweidimensionalen" Daten auf Papier und Biindern. Unter dem Stichwort "Visualization in Scientific Computing" [MDB87] hat eine Kommission der amerikanischen Supercomputerzentren ein Programm fur einen Ausweg aus diesem Dilemma entwickelt. Die Entwicklung von allgemein benutzbaren, interaktiven Visualisierungsprogrammen zur Darstellung groBer 3-dimensionaler Datenfelder und ihrer zeitlichen Veriinderung ist eine der vorgeschlagenen Methoden, da schon heute leistungsfiihige Graphik-Workstationen fur den allgemeinen Gebrauch in Forschungs- und Entwicklungseinrichtungen bereitstehen. Die bei interaktiver Nutzung der Visualisierungsmoglichkeiten auftretenden Datenflusse konnen mit Hilfe von Breitbandkommunikationsnetzen auf Glasfaserbasis bewrutigt werden.
Die, 3-dimensionale graphische Visualisierung in Wissenschaft und Technik- S-611-der
12
Information, dem besseren Verstandnis und der Erkenntnissteigerung dienen. Der interaktive Gebrauch der Visualisierungstechniken wahrend des Simulations- bzw. Me£prozesses wird gegenuber der blo£en Dokumentation immer wichtiger werden, d. h. sie mussen vor allem fur die Anwender selbst leicht handhabbar sein. Fur ein Softwarepaket zur Visualisierung 3-dimensionaler Datenfelder ergeben sich deshalb als wesentliche Anforderungen:
• Flexibilitat des Visualisierungsprogramms gegenuber unterschiedlichen 3-dimensionalen skalaren Datenarten und -groBen;
• leichte Handhabbarkeit auch fur Computergraphik-Nichtspezialisten durch Benutzung nur weniger, leicht interpretierbarer Steuerparameter;
• Ausnutzung der auf Graphik-Workstationen vorhandenen Routinen zur Datendarstellung, -transformation und -nachbehandlung;
• Moglichkeit der Berechnung von Erscheinungsbildern der Volumeneigenschaften auf beliebige Schnittebenen durch das Volumen;
• Moglichkeiten der Uberlagerung bzw. Differenzbildung von Datenfeldern durch Entkopplung der Berechnung von Lichtintensitat und -farbej
• Moglichkeit der Ankopplung von Methoden zur Kontrastverstarkung bei Datenfeldern mit schwachen lokalen Inhomogenitaten;
• effizienter, moglichst rekursiver Algorithmus zur Abbildung der Volumeneigenschaften auf eine 2-dimensionale Lichtintensitatsverteilung auf dem Bildschirm.
Beispiele fur 3-dimensionale skalare Felder, die in Wissenschaft und Technik untersucht werden, sind:
• Gravitationsfelder; • elektrische und andere Potentialfelder; • Dichteverteilungen von Flussigkeiten, Gasen, Plasmas; • Temperaturverteilungen; • Verteilungen von Verunreinigungen in der Atmosphare, in Fllissigkeiten, in fest en
Korpernj • einzelne Komponenten eines Vektorfeldes, z. B. die Vertikalgeschwindigkeit emer
Stromungj • Dichteverteilungen von Gewebe aus radiologischen Messungen.
Beispiele fur Visualiserungen von Volumeneigenschaften kommen aus so unterschiedlichen Bereichen wie der Astrophysik und Kosmologie [Ups86], der Meteorologie [Hus83), Geophysik [Sab88), Stromungslehre [UK88) und der' Humanmedizin [RM87, DCH88, Lev88). Die Begrenztheit von 2-dimensionalen, z. B. projektiven, Visualisierungsmodellen fur die Darstellungen von Messungen in der Meteorologie und Geophysik wird in [Hus83) diskutiert.
Das in dieser Arbeit vorgestellten Visualiserungsprogramm fur 3-dimensionale Datenfelder enthalt im wesentlichen drei Teile:
• Aufbereitung des vorgegebenen Datenfeldes (Definition und Normierung des 3-dimensionalen Gittersj Auswahl bzw. Erganzung, Normierung und Skalierung des gegebenen Volumendatenfeldesj optimale Speicherung);
13
• Abbildung der Volumeneigenschaften auf eine 2-dimensionale Lichtintensitatsverteilung auf den Rasterbildschirm (Transporttheorie-Modell fur die simulierten Lichtstrahlen; Abbildung der Volumeneigenschaften auf die Absorptions- und Streukoeffizienten und die verschiedenen Quellterme);
• N achbearbeitung und Darstellung der Bilddaten (Kontrastverstarkung durch Filterprozesse oder Falschfarben-Methoden; Glattung von Rauscheffekten hervorgerufen durch Datenfeldungenauigkeiten oder Diskretisierungsfehler; Darstellung des Erscheinungsbildes auf verschiedenen Schnittflachen).
Diese Arbeit konzentrierl sich auf den zweiten Punkt, die Erlauterung des Transporttheorie-Modells zur Visualisierung 3-dimensionaler skalarer Datenfelder. Es laBt sich zeigen, daB einige der schon vorhandenen Visualisierungsmodelle [RM87, Sab88, UK88, Lev88] als Spezialfalle enthalten sind. 1m folgenden Abschnitt wird die Methode des Volumen-Ray Tracing auf Grundlage der linearen Transporttheorie beschrieben. Die besonderen Vorteile dieses sehr allgemeinen Modells liegen dabei in den vielfaltigen Abbildungsmoglichkeiten der vorgegebenen Volumeneigenschaften auf die strukturell unterschiedlich wirkenden Parameterfelder Absorptions-, Streu- und Quellterm. Die Wirkungsweise dieser Terme innerhalb des Visualisierungsprozesses wird durch Testbilder dokumentiert. In einem Anhang werden kurz die Grundlagen der linearen Transporttheorie erlautert und Algorithmen fur die numerische Berechnung angegeben.
2. Visualisierung via Transportiheorie-Modell
Da der Visualisierungsalgorithmus moglichst universell einsetzbar sein solI, ist eine umfangreiche Vorbearbeitung des vorgegebenen Datenfeldes erforderlich. Einzelne Schritte sind:
• Wahrend des Einlesens werden die Daten entsprechend der Gitterbeschreibung sortiert und die Extremalwerte bestimmt;
• "zufallig" angeordnete Datenwerte werden durch Interpolationsalgorithmen auf ein angepaBtes Gitter abgebildet, urn "explite" Speicherung der Gitterkoordinaten zu vermeiden;
• die eingelesenen diskreten 3-dimensionalen Felddaten werden normiert, z. B. auf den Bereich [-1,1], d. h. das Datenfeld wird auf eine diskrete, normierle Volumendichte abgebildet.
Das Modell fur die Darstellung dieser 3-dimensionalen Volumendichte als 2-dimensionale Lichtintensitatsverteilung auf dem Bildschirm benutzt die Formalismen der linearen Transporttheorie fur die Lichtausbreitung in inhomogenen amorphen Korpern [Cha60, CZ67, DM79]. Auf dies em Modell beruhen auch moderne Algorithmen fur die Bildsynthese von Wolken [KV84] oder fur die Synthese komplexer "naturlicher" Szenen [Kaj86].
Fur die Anwendung dieses Modells auf das Problem der Visualisierung allgemeiner 3-dimensionaler Datenfelder sind drei Hauptbestandteile zu erlautern, das abstrakte Transporttheorie-Modell, die zugehorigen Auswertungsalgorithmen und die Moglich-
14
keiten der Abbildung des vorgegebenen Datenfeldes auf die Modellparametet. Die mathematisch-physikalischen Grundlagen des linearen Transporttheorie-Modells und der angepaBten Auswertungsalgorithmen werden im Anhang kurz erlautert.
Die lineare Transporttheorie beschreibt die Ausbreitung von Teilchen (Photonen, Elektronen, Neutronen, Ionen, usw.) in inhomogenen amorphen Korpern. Die Anzahl und Bewegungsrichtung der Teilchen in jedem Volumenelement geniigt dabei einer Integro-Differentialgleichung (Boltzmann Gleichung) (s. Formel Al im Anhang). Die Art der raumlichen Ausbreitung der Teilchen hangt dabei von den Parametern Absorption, Quellstarke und Streuung ab, die durch Eigenschaften der Volumendichte modifiziert werden konnen. Die Auswertung dieser Gleichung im Rechner erfolgt durch Diskretisierung der formalen Losung (A4) (s. Anhang).
Das hier vorgestellte Visualisierungsmodell simuliert die Ausbreitung "virtueller" Lichtstrahlen durch das 3-dimensionale Datenfeld, das auf eine Volumendichte abgebildet wurde. Dabei wird die Lichtintensitat und die Ausbreitungsrichtung an den Endpunkten eines jeden Wegelements gemaB Formel (A4) unter Benutzung von (AI3) und (AI4) in Abhangigkeit von der lokalen Absorption, Quellstarke und Streustarke berechnet (s. Abbildung 1). Ohne Beriicksichtigung von Streuungen ergibt sich der gerade Weg 1, wahrend bei Einschalten von Streuprozessen die Bahn 2 berechnet werden muB. Die Gesamtlichtintensitat auf dem Bildschirm ergibt sich dann aus der normierten Summe der Intensitat 18(1) und allen Wert en der Intensitat 18(2), die durch die Monte-Carlo Simulation von zusatzlichen Lichtstrahlen erhalten werden kann.
Die Anwendung des linearen Transporttheorie-Modells auf das Problem der Visualisierung "natiirlicher" oder abstrakter skalarer Datenfelder erfordert eine angepaBte Abbildung der vorgegebenen Information auf die charakteristischen Parameterfelder (Absorptionterm L;t, Streuterm L;8 und Quellterme q8 bzw. qv). Das vorgeschlagene Transporttheorie-Modell besitzt eine groBe Anzahl von Moglichkeiten, urn die Eigenschaften des normierten und skalierten Datenfeldes F(x) mit x = (x, y, z) auf diese drei Parameterfelder der Transporttheorie abzubilden.
1) Fiir nichtkomplexe, relative homogene Volumendichten bietet sich die einfache Abbildung auf den Absorptionsterm und einen Volumenquellterm mit Emissionsrichtung s
L;t(x; k) = cl(k)· F(x) qv(x, s; k) = c2(k) . F(x)
(1)
an [Sab88], wobei Cl 2:: 1 und C2 2:: 1 freie Parameter sind, die farbabhangig sein konnen. Diese beiden Terme konnen natiirlich auch proportional zum Betrag des Feldgradienten gewiihlt werden. Die Berechnung der Lichtintensitat auf dem Bildschirm erfolgt nach einem einfachen Rekursionsalgorithmus (s. Gleichung A13 im Anhang), wobei die Parameterwerte an den Wegelementenden durch entfernungsabhangige Interpolation iiber die benachbarten acht Dichtewerte ermittelt werden. Der Vorteil dieser Abbildung liegt in der moglichen einfachen Rechnerauswertung, nur die geraden Strahlen (s. Abbildung 1) werden verfolgt. Der Nachteil dieser Wahl besteht im Verlust an Tiefeninformation durch das Fehlen von Schattierungseffekten (s. auch Abbildung 2). Die Zuordnung (1) entspricht der Simulation einer Rontgen-Durchleuchtung fiir C2 = O.
15
I \
2 'I'-- -..... ......
----- r'- -
Bildschirm au Bere Quelle
Abb. 1: Beispiel fur mogliche Lichtwege innerhalb des Volumen-Ray Tracing
2) Fur das Hervorheben von starken Feldinhomogenitaten auf inneren Flachenelementen oder von z. B. Aquipotentialflachen bietet sich das Modellieren von Flachenladungen an. Sie sind definierl durch
(2)
wobei Xs die Flache beschreibt, es die lokale Oberflachennormale ist und die Intensitatsdifferenz zwischen den beiden Oberflachenseiten proportional zum VolumendichteGradienten senkrecht zur Oberflache gesetzt wird [Sab88]. Der Parameter C3 > 0 kann farbabhangig gewahlt werden.
Diese Abbildung eignet sich besonders fur Anwendungen in der Medizin, wo die Grenzen zwischen unterschiedlichen Gewebearten besonders hervorgehoben werden sollen [DCH88, Lev88]. Die Anwendung von Zuordnung (1) und (2) of CT-Daten in einem vereinfachten phanomenologischen Modell wird in [RM87] gezeigt. Auch interessante Aquipotentialflachen von komplexen Makromolekulen konnten mit dieser Methode besonders hervorgehoben werden, wenn man der entsprechnenden Potentialflache eine Flachenladung zuordnet (s. auch Abbildung 3).
3) Vorgegebene extreme Volumendichten oder Dichteschwankungen konnen auch durch Abbildung auf dem Streuterm in Gleichung (A3) besonders hervorgehoben werden.
16
Der Streuterm konnte in zwei wesentlichen Fiillen Bedeutung haben:
• Selektives Hervorheben von lokalen Volumendichte-Fluktuationen kann mit einem Volumenstreukoeffizienten 2;s erreicht werden durch die Identifizierung
2;s{x; k) = c4(k)· F(x), (3)
wobei C4 > 0 ist. Ein zusiitlicher Freiheitsgrad besteht in der Moglichkeit, den Winkelanteil des differentiellen Streuquerschnitts z. B. mit verstiirkter Vorwiirts- oder Riickstreuung zu wiihlen (s. Anhang). Dieser Fall ist besonders geeignet fiir die Visualisierung von atmosphiirischen Effekten (Wolken, Dunst, Verunreinigungen) [KVS4]. Aerosole (\Volken, Dunst) oder Regengebiete streuen und/oder brechen das durchfallen de Licht, so daf3 Punktlichtquellen diffus erscheinen oder Regenbogeneffekte auftreten. Auch Datenfelder mit Substrukturen in symmetrischen Anordnungen konnen durch ein Einbeziehen von Streuungen gemiif3 Formel (A6) und (A7) mit kleinem mittleren Streuwinkel deutlicher dargestellt werden ("channeling" -Effekt) [SUS9] .
• Eine Hervorhebung von Grenzfliichen zwischen Volumenregionen mit stark abweichender Dichte (z. B. bei CT-Anwendungen) kann durch Beriicksichtigung eines Reflexions- und/oder Transmissionskerns (s. Gleichungen AS-ll im Anhang) erreicht werden. Wird eine spiegelnde Reflexion auf einer interessanten Fliiche mit einer zusiitzliche Punktlichtquelle simuliert, ergeben sich yom Einfalls- und Ausfallswinkel abhiingige Schattierungseffekte, die den 3-dimensionalen Eindruck verstiirken (s. Abbildung 4). Die Definition einer Semitransparenz der Fliichen erlaubt die Abbildung von auch verdeckten Strukturen, der Tiefeneindruck wird ebenfalls verstiirkt.
Da die Gleichungen der Transporttheorie linear sind, konnen mehrere korrelierte Datenfelder (z. B. Volumendichte und -temperatur) gleichzeitig visualisiert werden, indem man die Modulation der Lichtintensitiit und die der Farbgebung entkoppelt. Elektrische Potentialfelder von verschiedenen Atomen in Molekiilen oder Datenfelder mit unterschiedlichen Vorzeichen konnen z. B. auf getrennte Frequenzbereiche bzw. RGBBereiche der Lichtintensitiit abgebildet werden. Ebenso liif3t sich auch die Differenz zweier Datenfelder (z. B. Ergebnisse aus Messungen und Simulationen) darstellen.
Als Beispiel fiir ein Datenfeld wurde die quantenmechanische Wahrscheinlichkeitsdichte fur den Aufenthaltsort des Elektrons in einem H-Atom (3d- und 4d-Konfiguration) gewiihlt. Abbildung 2 zeigt die Visualisierung mit Hilfe des Algorithmus (A13) und der Zuordnung (1). Abbildung 3 beriicksichtigt zusiitzlich die Wirkung einer Fliichenquelle gemiif3 Zuordnung (2). Die Rolle der spiegelnden Reflexion zur Hervorhebung des Tiefeneindrucks wird in Abbildung 4 dokumentiert, wobei in jedem Oberfliichenpunkt die Reflexion einer zusiitzlichen punktformigen Lichtquelle an einer "rauhen" Oberfliiche gemiif3 Formel (AlO) und den Algorithmen in [CTSl] simuliert wird.
Die Vorteile des vorgestellten Visualiserungsmodells fiir 3-dimensionale skalare Datenfelder liegen in seinem exakten mathematisch-physikalischen Fundament, dem schon vorhandenen Fundus an Auswertungsalgorithmen aus uber 40 Jahren Anwendung in unterschiedlichsten Bereichen und den vielfiiltigen Moglichkeiten, die vorgegebenen Daten
17
Abb. 2: Darstellung mit Zuordnung auf Absorptions und Volumenquellterm
Abb. 3: Darstellung mit Zuordnung auf Volumenquellterm und einen Flachenquellterm ("Visualisierung von Aquipotentialfiachen")
18
Abb. 4: Darstellung mit Zuordnung auf Fliichenquellterm und spiegelnden Oberflachenreflexionsterm
(Feldstiirke, Gradient, Kriimmung, usw.) aufunterschiedlich wirkende Modellparameter abzubilden. Viele der schon verwendeten "phiinomenologischen" Modelle sind in diesem Modell enthalten. Der groBte N achteil dieses aufwendigen Modells ist noch die erforderliche Rechenzeit, insbesondere bei Vorhandensein von feinkornigen Inhomogenitaten oder bei Benutzung von Streuungssimulationen.
A nhang: Lineare Transporttheorie fur die Lichtausbreitung
Mit Hilfe der linearen Transporttheorie kann die Ausbreitung von "Lichtstrahlen" in inhomogenen amorphen Korpern berechnet werden. 1m folgenden werden die wesentlichen Formeln und Algorithmen der linearen Transporttheorie erliiutert. Weitergehende Informationen sind in Textbiichern zu finden [Cha60, CZ67, DM79, Mar80].
1) Grundgleichungen der linearen Transporttheorie
Ein Lichtstrahl wird durch die lokale Lichtintensitat I(x,s;k) im Punkt x = (x,y,z)
19
mit der Ausbreitungsrichtung s und der Frequenz k (z. B. im einfachsten diskreten Fall k = R, G oder B) beschrieben.
Die riiumliche Ausbreitung der lokalen Intensitiit I folgt einer Bilanzgleichung in der Form einer linearen Integro-Differentialgleichung
(s· D) . I(x, s; k) + L;t(x; k) . I(x, s,; k) = Q(x, S; k), (AI)
wobei D der Nabla-Operator ist. Der Diimpfungs koeffizient L;t ist gegeben durch
(A2)
wobei L;a den lokalen Absorptionskoeffizienten pro Wegliingeneinheit und L;s den totalen Streuquerschnitt pro Volumeneinheit bezeichnet.
Der verallgemeinerte Quellterm Q hat die Form
Q(x,s;k) = q(x,s;k)
+ L;s(x; k) J dk' J ds' K(x, s' --+ S; k' --+ k)· I(x, s'; k'), (A3)
wobei q eine punkt-, fiiichen- und/oder volumenfarmige Lichtquelle beschreibt. Kist der differentielle Streuquerschnitt, der die Winkelstreuung und die Frequenzverschiebung des Lichts generiert.
Die Lasung der Gleichung (AI) kann in Form einer Integralgleichung
I(x, S; k) = Is(x - Rs,s; k)· exp [-peR)]
+ J dR' Q(x - R's, S; k) . exp [-peR')] (A4)
gegeben werden. Hierbei bezeichneit R die Wegliinge und Is ist die durch die Volumenoberfiiichen im Punkt Xs = x - Rs einfallende Lichtintensitiit. Die effektive optische Wegliinge wird definiert durch
R
peR) = J dR' . L;t(x - R's; k). o
(A5)
Der zweite Anteil des Quellterms in (A3) beschreibt die Streuung des einfallenden Lichts pro Wegliinge. Grundsiitzlich kann man zwei Arten von Streuung unterscheiden:
a) Volumenstreuung, wobei L;s(x; k) z. B. proportional zur Volumendichte gewiihlt werden kann. Der Kern K beschreibt die Winkelverteilung des gestreuten Lichts z. B. in den Formen
{
417r fur isotrope Streuung, K(x,s' --+ s) = 8(s + s') fur Ruckstreuung,
K HO ( cos 8) fur mittlere Streuwinkel. (A6)
20
Hierbei ist J(HG die Henyey-Greenstein Approximation [HG41] fiir die allgemeine Streufunktion J(
r 1- g2 li. HG ( cos 0) = ( 0 2)3/2 1 - 2g . cos ~ + g
(A7)
mit dem Streuwinkel cos 0 = (s . s') und g mit jgj ::; 1 entspricht dem Mittelwert von cos 0.
b) Oberfliichenstreuung an einem Fliichenpunkt Xs kann beschrieben werden durch die Relation
-(es · Sout)· jOut(xs,Sout; k) = L:s(xs; k)· J dk' J dsin ·
(e s . sin)I<AXs, Sin -+ Sout; k' -+ k)· in(xs, Sin; k'), (AS)
wobei e s die Fliichennormale ist. Die Fliichenstreufunktion Ie hiingt nur von den angenommenen lokalen Eigenschaften der betrachteten Fliiche abo J(s kann in der Form
J(s = J(spec + J(back + J(diff
modelliert werden. Hierbei bezeichnet J(spec den spiegelnden Reflexionsanteil
J(spec = J(s(x s; k' -+ k)· 8(Sin - ss)
(A9)
(A10)
mit der Reflexionsrichtung Ss= Sin - es · (Sin· e s). Urn ausgedehnte Glanzlichter auf den Oberfliichen zu erhalten wird die Delta-Funktion in (A10) z. B. durch eine gau:f3formige Verteilung geniihert, urn die Reflexion an rauhen Oberfliichen zu simulieren [eTS1]. Der Riickstreuanteil J(back kann durch
(All)
beschrieben werden. Der Anteil J(diff kann iiber aIle Raumwinkel gleichmii:f3ig verteilte Riickstreuung beschreiben.
Grundsiitzlich konnen immer auch farbverschiebende (inelastische) Streukeme eingefiihrt werden.
2) A 19orithmen zur rechnergestutzten A uswertung
Die Gleichung (A4) bildet die Grundlage fiir den Volumen-Ray Tracing Algorithmus, der das 3-dimensionale Datenfeld (Volumendichte) auf den Bildschirm abbilden solI. Eine einfache Ray-Tracing-Konfiguration besteht z. B. aus einer ebenen Lichtquelle im Unendlichen, von der parallele Lichtstrahlen ausgehen, die durch das betrachtete Volumen moduliert und auf dem Bildschirm als Intensitiitsverteilung visualisiert werden (s. Abb. 1).
Die lokalen Werte der Transporttheorie-Parameter (L:t, L:s und q) an den Endpunkten der Wegelemente konnen durch Interpolation der Werte auf den benachbarten acht Gitterpunkten mit entsprechender Abstandsgewichtung erhalten werden.
21
Als Lasungsansatz fur die Gleichung (A4) nimmt man die Entwicklung (Neumann'sche Reihe) [Mar80, Kaj86]
n
I = Io + I:: Ii, (A.12) i=l
wobei Io die Lasung von Gleichung (A4) ohne den Streuterm in (A3) ist. Die Lasung Io entspricht einem geraden Lichtweg durch das Volumengitter (s. Weg 1 in Abb. 1). Die Entwicklungsterme Ii entsprechen dann Lichtwegen mit i-facher Streuung (s. Weg 2 in Abb.1).
Die Berechnung von I o , das Verfolgen von geraden Lichtstrahlen durch die Volumendichte, kann bei Benutzung der Trapezregel durch die einfache Rekursionsregel
Io(s + ds) = [Io(s) +.5· q(s)· ds]
. exp[-.5(~t(s) + ~t(s + ds))· ds] +.5· q(s + ds)· ds (A13)
erfolgen, vrobei Io(O) gleich dem Anfangswert Is gesetzt wird. Die Lange des Wegelements ds mU£ dabei der GraBenordnung der Volumendichte-Inhomogenitaten angepaBt werden.
Die multiplen Streuterme Ii folgen gemaB (A3) und (A4) aus der Rekursionsformel
Ii(x, S; k) = J dR' . exp[-P(R')]· Qi-1(R'; k),
wobei Qi-1 gegeben ist durch
Qi-1 = ~s(x - R's; k)· J dk' J ds'· K(x - R's,s' --t S; k' --t k)
. Ii- 1(X - R's,s';k').
(A14)
(A15)
Die Auswertung dieser Rekursionsformel kann bei nichttrivialem Streuquerschnit K mit Hilfe der Monte-Carlo-Methode erfolgen [Mar80, Kaj86], wobei die Weglange zwischen zwei Streuungen z. B. durch die mittlere Volumendichte bestimmt wird.
22
Literaturverzeichnis
[CZ67]
[Cha60] [CT81]
[DCH88]
[DM79] [HG41] [Hus83]
[Kaj86]
[KV84]
[Lev88]
[Mar80]
[MDB87]
[RM87]
[Sab88]
[SU89]
[UK88]
[Ups86]
K.M. Case, P.F. Zweifel, Linear Transport Theory, Addison-Wesley, Reading, 1967 S. Chandrasekhar, Radiative Transfer, Dover, New York, 1960 R.L. Cook, K.E. Torrance, A reflectance model for computer graphics, Computer Graphics 13, No.3 (1981), 307-316 RA. Drebin, L. Carpenter, P. Hanrahan, Volume rendering, Computer Graphics 22, No.4 (1988), 65-74 J.J. Duderstadt, W.R Martin, Transport Theory, Wiley, New York, 1979 L.G. Henyey, J.L. Greenstein, in Astrophys. J. 93 (1941), 70 K.J. Hussey et al., Environmental data display, Environ. Sci. Technol. 17, No.2 (1983), 163-170 J.T. Kajiya, The rendering equation, Computer Graphics 20, No.4 (1986), 143-150 J.T. Kajiya, B.P Van Herzen, Ray tracing volume densities, Computer Graphics 18, No.3 (1984),165-174 M. Levoy, Display of surfaces from volume data, IEEE CG & A, May 1988, 29-37 G.l. Marchuk et al., The Monte Carlo Method in Atmospheric Optics, Springer-Verlag, Berlin 1980 B.H. McCormick, Th.A. DeFanti, M.D. Brown, eds., Visualization in scientific computing, Computer Graphics 21, No.6 (Nov. 1987) G. Russel, RB. Miles, Display and perception of 3-D space-filling data, Appl. Optics 26, No.6 (1987), 973-982 O. Sabella, A rendering algorithm for visualizing 3D scalar fields, Computer Graphics 22, No.4 (1988), 51-58 A.H. S¢rensen, E. Uggerh¢j, The channeling of electrons and positrons, Scientific American, (June 1989) C. Upson, M. Keeler, V-buffer: Visible volume rendering, Computer Gra-phics 22, No.4 (1988), 65-74 C. Upson, The visual simulation of amorphous phenomena, The Visual Computer 2, (1986), 321-326
Schalldruckberechnung im PKW-Innenraum
H. Burfeindt, Volkswagen AG Wolfsburg H. Zimmer, SFE GmbH Berlin
(Gesellschaft fur Strukturanalyse in Forschung und Entwicklung)
Es werden die Ergebnisse eines auf dem Finite Elemente Programm NASTRAN basierenden Berechnungsverfahrens zur PKW-Innenraumakustik im Vergleich zu Messungen dargestellt.
Das Verfahren berucksichtigt die Interaktion zwischen Blechstruktur und eingeschlossenem Luftvolumen des Innenraumes, die D5mpfung der Luftschallwellen Uber die meBbaren Impedanzen von absorbierender Innenauskleidung und die Wirkung von Entdrehnungsmaterialien.
1. Einleitung
Akustische Eigenschaften von Fahrzeugen sind zur zeit erst in einer sp5ten Prototypenphase experimentell bestimmbar.
Es ist wUnschenswert, durch geeignete Berechnungsverfahren in einer frUhen Entwicklungsphase des Fahrzeugs Aussagen und prognosen zur Innenraumakustik machen zu kennen. Ein meglicher Ansatz hierzu ist die Methode der Finiten Elemente, deren Anwendung auf komplexe Bauteile zur Ermittlung des strukturmechanischen Verhaltens heute Stand der Technik ist.
Ein solches Verfahren lieBe sich da£~ schon sinnvoll einsetzen, wenn berechnete Schallpegel die gleichen Charakteristiken wie gemessene zeigen und tendenzm5Big hinsichtlich Strukturver5nderungen qualitativ "richtige" Schalldruckpegel berechnet werden kennen.
Damit g5be es die Meglichkeit, TrendausSagen hinsichtlich des zu erwartenden Schalldrucks bei Parameter5nderungen an einzelnen Bauteilen der Karosserie rechnerisch zu prognostizieren.
2. Berechnungsverfahren
Die Berechnung des Schalldrucks im Innern eines PKW's simuliert zwei schwingungsf5hige Systeme, die miteinander gekoppelt sind (1, 2):
- die Fahrzeugstruktur, an der die Erregerkrafte angreifen und deren Sterung als Kerperschall weitergeleitet wird.
- den von der Fahrzeugstruktur eingeschlossenen Luftraum, an dessen Berandung der Kerperschall Luftschallwellen abstrahlt, die einen Schalldruckpegel erzeugen.
Bei Vorhandensein von Innenausstattung ist der Luftraum zusatzlich teilweise von schallabsorbierendem Material berandet.
I)Berechneter bzw. gemessener Schallpegel bedeuten in diesem Bericht die bezogene Ubertragungsfunktion vom Karosserie-Anregungspunkt zum Aufpunkt im Luftinnenraum.
24
Die L6sung des gekoppelten Gleichungssystems in den physikalischen Koordinaten
{pI - Druckvektor der Luftknoten
{wI - Verschiebungsvektor der Strukturknoten
ware wegen der groBen Zahl von Freiheitsgraden (insbesondere des Strukturmodellsl zu aufwendig.
Daher wird eine Transformation in den Modalraum des Struktur- bzw. Luftmodells durchgeftihrt tiber die Verkntipfungen:
{pI - [+L] {ll}
{wI - [+s] Ill}
Dabei bedeuten:
I llL}
Ills}
[+L]
[+LA]
- modaler Druckvektor
- modaler Verschiebungsvektor
L L - [VI' V2 ' ... ] - Modalmatrix der Luft
[ "'ILA, ..., LA ] I 12' .. Modalmatrix der Luft im
Bereich Luft/Absorber
s s - [VI' Y2 ' ... ] - Modalmatrix der Struktur
Nach Ausftihrung der Transformation erhalt man als gekoppeltes Gleichungssystem in den modalen Freiheitsgraden der Luft {llL} und der Struktur Ills}:
1 [S] - - [<t>Ll IS] [4lJ
f!
[Kl - [<t>~l IK] ItA;]
1 Iii - -,- !<t>[J[PJ 14lJ
C-f!
1M] - I<t>~] 1M] [tA;]
- Modale Sleifiglell der Luft
Modale Sleifiglell der Slrullur
- Modale Mass< der Lufl
- Modale Mass<: der Slrullur
B(w) [SA] - --!<t>L.J IS,J IlPul
Z, (w)
- Modale Sleifigkeil des Absorbers
til - I<t>~l 1f1
[e] - [<t>[lI8l1tA;]
- Modale Masse des AI>sorbers
- Modaler Laslveklor
- Modale Kopplungsma-tru.
25
Die GraBen Z (w) und B (w) beschreiben die.Absorptionswirkung schaildampfender Materialien der Innenausstattung.
ZA(W) ist die lokal wirkende, komplexe akustische Impedanz bei senkrechtem und B (w) die flachenhaft wirkende, komplexe akustische Impedanz bei schragem Einfall der Schallwellen (3). Beide GraBen sind meBtechnisch bestimmbar.
Zusatzliche Dampfungswirkung der Struktur kann als sog. strukturdampfung und modale Dampfung berucksichtigt werden.
Lokal dampfende Bereiche der Struktur, wie etwa die Gummi-Einfassungen der Scheiben, werden durch diskrete viskose Dampfer beschrieben.
Die karperschalldampfende wirkung von auf Blechflachen aufgeklebten Schwermatten kann berucksichtigt werden, indem der betreffende Bereich als Sandwich-Platte mit Strukturdampfung und gekoppelter Biege-Membranspannung aufgefaBt wird.
Mit dem in Gl. (1) dargestellten Gleichungssystem wird die Erregungsrechnung zur Ermittlung des Schalldrucks durchgefuhrt.
3. Finite Element Modellierung
Urn die anfallende Datenmenge zu begrenzen, werden die Schalldruckberechnungen an einem halben Karosseriemodell durchgefuhrt. Symmetrische und antimetrische tasungen werden zum SchluB uberlagert zum Gesamtergebnis.
Dabei wird eine symmetrische Fahrzeugstruktur vorausgesetzt, was nur naherungsweise erfullt ist.
Das Finite-Elemente-Modell der Rohkarosserie zeigt Bild 1. Es besteht aus ca. 11 000 Schalenelementen und 10 400 Knoten und hat ca. 62 000 Freiheitsgrade.
Turen und Heckklappe sind nicht als Finite-Elemente-Modell nachgebildet, sondern als modale Modelle aus experimenteller Modalanalyse in die Rechnung mit einbezogen.
Daraus ergab sich die Forderung nach Aufteilung der Fahrzeugstruktur in Superelemente. Deren Residua1struktur besteht aus ca. 3 900 Freiheitsgraden.
Das in Bild 2 dargestellte Finite-Elemente-Modell der tuft besteht aus ca. 1 280 Volumelementen mit ca. 1 500 Freiheitsgraden.
Bild 3 zeigt das in den Rechnungen verwendete Finite-Elemente-Modell der absorbierenden Flachen, bestehend aus ca. 900 Schalenelementen entsprechend ca. 1 000 Freiheitsgraden.
26
4. Ablauf einer.Schalldruckberechnung
Die auf NASTRAN basierende Berechnungskette lauft in drei Schritten ab:
- Erzeugung des Lastvektors {F} und des Koppellastvektors zur Bestimmung der Flachenmatrix [A] zwischen Luft/Struktur und Luft/Absorber in SOL 61.
- Reduktion der Superelemente von Struktur und Luft und Berechnung der Eigenwerte in SOL 63.
Bestimmung der Modalmatrizen [~L] und [~S] der Luft bzw. der Struktur und der Kopplungsmatrix [8].
Bei der Kopplung Struktur/Luft brauchen die Strukturknoten nicht auf korrespondierende Luftknoten zu fallen. Ihre wirkung auf umgebende Luftknoten wird durch EinfluBfaktoren entsprechender Gewichtung beschrieben.
- Aufstellen der modalen Steifigkeits- und Massenmatrizen [SA] bzw. [P ] des Absorbermaterials und Berechnung der erzwungenen ~chwingung des gekoppelten System Struktur/Luft (Gl. 1) einschlieBlich Dampfung durch Absorbermaterial und Struktur durch Erregung mittels der modalen Last.
Das Ergebnis ist der Modalvektor {nt } bzw. nach erfolgter Rticktransformation der Druck {pj im Hohlraum. All dies geschieht in SOL 71.
Die aus der symmetrischen bzw. antimetrischen Rechnung erzielten Ergebnisse werden anschlieBend tiberlagert, und man erhalt den Schalldruck im Innern des kompletten Fahrzeugmodells.
5. Ergebnisse der Rechnungen und Vergleich mit Messungen
Ftir den Vergleich Rechnung/Messung wurden eine Rohkarosserie des Golf schrittweise bis zu vollstandigen Innenausstattung aufgertistet und die Rechnungen an einem dem jeweiligen Ausbauzustand der Rohkarosserie entsprechenden Berechnungsmodell durchgeftihrt. Der Vergleich erfolgte mit Messungen an der Rohkarosserie im dazugehorigen Ausbauzustand (4).
Folgende Ausbaustufen wurden gewahlt:
(1) offene Rohkarosserie ohne Scheiben, Ttir und Heckklappe
(2) geschlossene Rohkarosserie mit Front-, Seiten- und Heckscheibe sowie Tliren und Heckklappe ohne sonstige Innenausstattung
(3) wie Stufe (2) und zusatzlich Vorder- und Hintersitze
(4) wie Stufe (3) und zusatzlich mit Schwermatten-Belegung im Querwand- und vorderen Bodenbereich
(5) geschlossene Rohkarosserie mit vollstandiger Innenausstattung
27
Rechnungen und Messungen erfolgten im Frequenzbereich 20 -200 Hz. Einleitungsort fUr die Krafterregung (Sinus-Kraft von 1 N) war vorn im linken Langstrager. Im folgenden werden die Ergebnisse der jeweiligen Ausbaustufe diskutiert.
5.1 Offene Rohkarosserie in Ausbaustufe (1)
Berechnete Eigenwerte, Schwingungsmoden sowie tibertragungsfunktionen wurden mit den Ergebnissen experimenteller Modalanalyse verglichen. Eine teilweise eindeutige Zuordnung der berechneten und gemessenen Schwingungsformen konnte bis ca. 60 Hz erzielt werden. FUr h6here Frequenzen wird die Zuordnung wegen der vermehrt auftretenden lokalen Schwingungsmoden und h6heren Modendichte immer schwieriger.
Tab. 1 zeigt eine GegenUberstellung berechneter und gemessener Eigenfrequenzen und Moden.
Ein Beispiel einer berechneten tibertragungsfunktion im Vergleich zur Messung ist in eild 4 dargestellt.
Die Rohkarosserie ist in der Rechnung modal gedampft, wobei die Dampfungswerte aus dem Experiment stammen.
5.2 Geschlossene Rohkarosserie in Ausbaustufe (2)
Die Rohkarosserie enthalt in diesem zustand nur Scheiben sowie TUren und Heckklappen, umhUllt damit einen geschlossenen Luftraum. Die elastische Lagerung der Scheiben ist durch Federelemente simuliert, Dampfung der Gummieinfassung wird durch viskose Dampfer modelliert. TUren und Heckklappe sind als modales Modell dargestellt. tiber SOL 43 (Component mode synthesis) werden die aus experimenteller Modalanalyse erhaltenen modalen Gr6Ben wie Masse, Steifigkeit, Dampfung unter Benutzung der Eigenformen, Gesamtmasse und Tragheitsmomente in Steifigkeits-, Massen- und Dampfungsmatrizen umgewandelt und in NASTRAN eingelesen. Die Modalmessungeh an TUr bzw. Heckklappe erfolgten am ungefesselten System.
Die Berechnung der Eigenfrequenzen und Moden des Luftraumes ist unproblematisch und bedarf keines besonderen Aufwandes.
Zu vergleichszwecken wurden diese auch experimentell bestimmt, die Ergebnisse sind in Tab. 2 den Resultaten der Rechnung gegentibergestellt. Bild 5 zeigt als Beispiel die niedrigste im Luftraum auftretende Schwingungsmode.
Zur UberprUfung des Rechenmodells wurden wiederum Eigenfrequenzen, Moden und tibertragungsfunktionen der geschlossenen Rohkarosserie ermittelt.
Die Messungen zeigen, daB durch die Ausstattung der Rohkarosserie mit TUren, Heckklappe und Scheiben bei tiefen Frequenzen hohe Verlustfaktoren durch System~ dampfung infolge von Gummidichtungen, Scharnieren etc. auftreten. Die Resonanzfrequenzen verlagern sich nur
28
geringfugig zu hoheren Frequenzen, jedoch kann sich die geanderte ortliche Steife- und Massenverteilung lokal durchaus in einer Veranderung der Modenform auswirken.
Tab. 3 enthalt einige berechnete Moden im Vergleich zur Messung, Ubertragungsfunktionen sind im Bild 6 und 7 dargestellt.
Bei der Berechnung des Schalldrucks im Fahrzeuginnenraum wird die vollstandige Kopplung zwischen Luftkorper und umgebender Blechstruktur berlicksichtigt, Ruckwirkungen der Luft auf die Struktur werden also mit erfaBt.
Tliren und Heckklappe existieren als modales Modell und sind uber Scharnier und SchloB mit der Struktur verbunden. Eine Ankopplung an den Luftraum besteht nicht.
Bild 8 und 9 zeigen den berechneten Schalldruck im vergleich zur Messung bei Anregung durch eine Sinuskraft am vorderen Langstrager.
Die Rechnung zeigt im oberen Frequenzbereich bei ca. 130, 140 und 170 Hz Uberhohungen gegenliber der Messung. Die Abweichung bei 170 Hz ist zum Teil in der zu weichen Anbindung des Daches an die Quertrager, modelliert durch elastische Federn zur Nachbildung einer Klebung, begrlindet.
Anhand von Polardiagrammen kann die Beteiligung der Karosseriemoden am Schalldruck nach Betrag und Phase graphisch sichtbar gemacht werden. Karosseriemoden, die den Schalldruck verstarken, sol len gezielt durch Strukturveranderungen an der Karosserie beeinfluBt werden, urn so den Schalldruck zu reduzieren.
AuBerdem ist es moglich, den EinfluB der Karosseriezonen am erzeugten Schalldruck graphisch zu verdeutlichen. Eine solche Darstellung liefert Hinweise uber Strukturbereiche, die als Schallquelle bzw. -senke wirken, wie in Bild 10 gezeigt.
5.3 Geschlossene Rohkarosserie in Ausbaustufe (3)
Diese Ausbaustufe entspricht der geschlossenen Rohkarosserie mit eingebautem vorder- und Hintersitz.
Zur Beschreibung der Absorptionswirkung des Sitzbelages ist nach Gl. (1) die Kenntnis der Impedanzwerte ZA(oo) bzw. B(oo) erforderlich. Diese wurden experimentell ermittelt.
Im Schwerpunkt der Sitze und Lehnen angebrachte konzentrierte Massen beschreiben im Rechenmodell die dynamische Massenwirkung von Vorder- und Hintersitz.
Wie aus Bild 11, 12 zu ersehen, ist der EinfluB der Sitze auf den Schallpegel gering.
Den berechneten Schallpegel im vergleich zur Messung zeigt Bild 13, 14. Die Zuverlassigkeit der gemessenen
29
Impedanzwerte, insbesondere fur die Bulkimpedanz B(oo), ist nicht bekannt.
Die bestehenden Diskrepanzen zwischen Rechnung und Messung sind zum Teil auf diese Unsicherheit zuruckzufuhren.
5.4 Geschlossene Rohkarosserie in Ausbaustufe (4)
Zusatzlich zu den Sitzen enthalt die Rohkarosserie k6rperschalldampfende Schwermatten im Querwand- und vorderen Bodenbereich.
Die Elemente der Strukturflachen, die mit diesem Material bedeckt sind, werden als "Sandwich" behandelt und uber NASTRAN-PCOMP-Karten beschrieben. NASTRAN erzeugt hieraus entsprechende PSHELL und MAT2 "balkdata" Karten fur die jeweiligen Elemente.
E-Modui und Verlustfaktor fur das Schwermatten-Material wurden experimentell ermittelt.
Der Dampfungseffekt ist auch hier recht gering, eine deutliche Reduktion des Schallpegels ist nur im oberen Frequenzbereich etwa bei 130 bis 140 und 170 Hz erkennbar.
Bild 17, 18 zeigen den berechneten Scha1lpegel im Vergleich zur Messung.
5.5 Geschlossene Rohkarosserie in Ausbaustufe (5)
Jetzt enthalt die geschlossene Rohkarosserie die gesamte innere Ausstattung einschlieBlich der Armaturentafel. Die Massenbeitrage des jeweiligen Ausstattungsmaterials werden durch Angabe der "nonstructural mass" auf den PSHELL "bulk data" Karten der entsprechenden Elemente berucksichtigt.
Die Impedanzen ZA (00) und B (00) der schallabsorbierenden Materialien wie Teppich und Poisterung wurden experimentell ermittelt.
Ihre dampfende Wirkung ist wesentlich starker als die des Materials zur Dach- und Seitenverkleidung, das eine glatte, harte Oberflache besitzt.
In der gegenwartigen Rechnung beeinfluBt der Teppich aufgrund seiner Masse und absorbierenden Oberflache den Schalldruck.
Einige Messungen scheinen darauf hinzudeuten, daB der Teppich auf einer dunnen Luftschicht liegend gedacht werden kann und so ein schwingungsfahiges System bildet, besonders im h6herfrequenten Bereich. Dadurch wird das akustische Verhalten zusatzlich beeinfluBt.
Die Dampfung des schallpegels bei Vorhandensein der gesamten Innenausstattung ist signifikant besonders im h6heren Frequenzbereich von 130 bis 200 HZ, wie aus Bild 19 und 20 ersichtlich wird. Hier ist der Schalldruck fur die FaIle mit/ohne vollstandige Innenaus-
30
stattung dargestellt. Sowohl Rechnung als auch Messung zeigen deutliche Pegelabsenkungen besonders bei 130, 170 und 190 Hz. TendenzmaEig ergibt sich aus der Rechnung also ahnliches (wenn auch starkeres) Dampfungsverhalten wie bei der Messung.
Bild 21, 22 zeigen den berechneten Schallpegel im FaIle vollstandiger Innenausstattung im Vergleich zur Messung an den positionen des Fahrer-/Beifahrerohres.
Auffallend ist die Diskrepanz zwischen Rechnung und Messung im Frequenzbereich 140 - 180 Hz in Bild 22 (Position Beifahrerohr).
Moglicherweise spielt hier die Tatsache eine Rolle, daE die Berechnung an einem Rechenmodell durchgeflihrt wurde, das Symmetrie voraussetzt, die Karosserie mit der Innenausstattung sich jedoch unsymmetrisch verhalt.
Die Messungen selbst zeigen in dem o. g. Frequenzbereich deutliche unterschiede an den Positionen Fahrer-/Beifahrerohr, wie aus Bild 23 zu entnehmen ist.
6. Zusammenfassung und Bewertung
Das beschriebene Berechnungsverfahren berlicksichtigt die wechselseitige Interaktion zwischen Karosserie und Luftinnenraum und die Dampfung des Luftschalles durch porose Absorberstoffe. Als Dampfungsmechanismen in der Struktur werden die modale Dampfung und lokal wirkende Dampfer (z. B. zur Beschreibung dampfender Gummieinfassungen) berlicksichtigt.
Auf dem Karosserieblech zur Korperschalldampfung aufgeklebte Schwermatten werden als Sandwichplatte mit Struktur-Dampfung behandelt.
Aufgrund bisheriger Erfahrungen hat sich gezeigt, daE Finite Element Berechnungen flir die Prognose akustischer Trendaussagen bei Parameterstudien an einzelnen Bauteilen heute schon unterstlitzend zum Experiment eingesetzt werden konnen.
Der sich jetzt bietenden Moglichkeit, rechnerisch die Beteiligung der Karosserieflachen und -moden am Schalldruck bestimmen zu konnen, urn daraus Rlickschllisse auf MaEnahmen zur Verminderung des Schalldruckpegels abzuleiten, wird dabei groEe Bedeutung beigemessen.
Dagegen ist der absolute Schalldruck im Fahrzeug zur zeit noch nicht zuverlassig im Sinne einer Prognose berechenbar.
Dies liegt unter anderem an den Grenzen der mathematischen Modellierbarkeit von zufallig variierenden Fertigungsdetails der Karosserie, die dem Rechenmodell nicht bekannt gemacht werden konnen, das Ubertragungsverhalten jedoch beeinflussen.
31
Hierzu gehoren auch beispielsweise Undichtigkeiten der Tlirgummileisten, die einen Druckausgleich nach au~en ermoglichen, der im Rechenmodell aber nicht berlicksichtigt werden kann.
Hinzu kommen Vereinfachungen bei der Erstellung des Finite Elemente Modells, wie zum Beispiel die Annahme kontinuierlicher Schwei~verbindungen an stelle von Schwei~punkten am realen Modell oder zu grober Vernetzung einzelner Bereiche und Vernachlassigung von Sicken oder Lochern.
Schlie~lich ist die Frage zu stellen, ob au~er den oben angegebenen Dampfungsmechanismen noch weitere mit in Betracht gezogen werden mlissen, wie etwa die vernachlassigte Schallabstrahlung nach au~en oder Dampfung durch doppelwandige Strukturen wie Dachhimmel oder Tliren.
Literatur:
(1) Burfeindt, Sarfeld, Zimmer:
(2) Burfeindt, Zimmer:
(3) Bliss:
Finite Elemente Berechnungsverfahren zur PKW-Innenraumakustik unter Berlicksichtigung von Dampfungsmechanismen durch Absorber, VDI-Bericht 537
Calculating sound pressure in car interiors, Proceedings of the MSC/NASTRAN European User's Conference, London 1989
A Study of Bulk Reacting Parous Sound Absorbers and a New Boundary Condition for Thin Porous Layers. 1981
(4) Die Messungen wurden an der TU Berlin (Institut flir Technische Akustik, Prof. Heckl) durchgeflihrt und sind in verschiedenen Berichten dokumentiert.
Die Berechnungen wurden mit Unterstlitzung der Firma SFE, Berlin, durchgeflihrt.
32
Eigenwert Frequenz Nr. Form Rechnung Messung
1 ANTI Dachscheren + Rotation 29.1 26.3
2 SYM Reserveradmu1de (loka1) 31.5
3 ANTI Querschwingung + TORSION 32.1 30.6
4 ANTI 1. Globale Torsion 38.5 35.8
5 SYM Globale Biegung, starkes 39.3 35.9 Abknicken des Vorderwagens
6 SYM 1. Globale Biegung 40.7
Tabelle 1
Berechnete und gemessene Eigenfrequenzen/Moden der offenen Rohkarosserie
Mode Experiment Rechnung Nr. Eigenfrequenz Hz Eigenfrequenz Hz
* 1 18.05 0
2 74.83 70.98
3 133.75 131.54 und 127.39
4 150.74 152.16
5 171.16 166.85
6 179.60 177 • 69
7 181. 78
Tabelle 2
Berechnete und gemessene Eigenfrequenz der Luftraummoden
* Die erste Mode stellt eine tiber dem ganzen Volumen konstante Druckverteilung dar.
Die Messung ergibt nicht Null, da in der Realitat die Berandung des Luftvolumens flexibel und nicht starr ist, wie in der Rechnung vorausgesetzt.
Eigenwert Nr.
1 ANTI
2 SYM
3 ANTI
4 SYM
Tabelle 3
33
Form
Torsion des Vorderwagens gegen Kabine
Reserveradmulde (lokal)
Globale Torsion des Fahr-zeugs
Globale Biegung und Dach-schwingung
Frequenz Rechnung Messung
30.1 29.6
31. 7
34.2 31. 6
41. 93 37.0
Berechnete und gemessene Eigenfrequenzen/Moden der geschlossenen Rohkarosserie ohne Ausstattung
34
A n han 9
- Rechnerkonfiguration und Rechenzeitbedarf
Die Berechnungen sind auf einer CRAY X-MP 14 mit VAX 8800 und CDC CYBER 990 als Frontendrechner durchgefuhrt worden. Diese sind mit der CRAY uber Hyperchannel verbunden (vgl. Bild 23). Die Rechenzeit fur eine vollstandige Analyse per NASTRAN V65, d. h. Berechnung des symmetrischen bzw. antimetrischen Lastfalles und uberlagerung betragt ca. elf bis zwolf cpu/h.
Die Verweilzeit fur eine Analyse belauft sich auf durchschnittlich 32 h.
Bei der Berechnung fallen relativ groBe Datenmengen an (Restartinformation) in der GroBenordnung von ca. 1,2 GB, die auf den CRAY-Platten uber langere Zeit nicht abgelegt werden konnen.
Diese werden stattdessen auf den platten der IBM 3090 gespeichert. Bild 24 zeigt benutzte Software und die Rechnerumwelt.
- Anforderungen aus Sicht des Anwenders
Die Berechnung lauft uber eine Anzahl verketteter Jobs ab, wobei der jeweils fertig werdende Job den darauffolgenden startet.
In Zukunft ware anzustreben, die "Turnaroundzeit" fur eine Analyse bei ModellgroBen wie sie momentan vorliegen (- 60 000 Freiheitsgrade) dahingehend zu reduzieren, daB das Ergebnis einer Berechnung, die heute gestartet wurde, morgen vorlage.
Momentan stehen maximal 65 500 Sektoren als SSD einem Job zur Verfugung. Jedoch wurde ein Job mit einer solchen SSDAnforderung lange in der Input-Queue stehen mussen. Daher wird der SSD-Bedarf fur einen Job der Berechnungskette auf max. 30 000 Sektoren begrenzt. Als Nachteil handelt man sich dann relativ hohe I/O-Zeiten, die in die GroBenordnung des CPU-Bedarfs kommen, ein. Daher ware eine Erweiterung der SSD-Kapazitat anzustreben. Eine Erhohung des Kernspeichers wurde ebenfalls den I/O-Bedarf reduzieren.
Weiterhin fallen, wie oben erwahnt, sehr groBe Datenmengen an, die als Restartdaten auf den IBM-Platten abgelegt werden. Bei Ausfuhrung von Restart jobs zur Datenauswertung ist i. a. sehr wenig CPU-Zeit notwendig, und die meiste Zeit wird zum Datentransfer IBM <--) CRAY benotigt.
35
Angestrebt wird daher eine Datenha1tung auf den CRAY-P1atten, wobei gleichzeitg versucht werden soll, durch spezie1-1es DMAP in NASTRAN den Antei1 der ausgeschriebenen Restartdaten weiter auf ein notwendiges Minimum zu reduziereno
Abbildungen
Abb. 1:
Abb. 2:
Abb. 3:
Abb. 4:
Abb. 5:
Abb. 6, 7:
Abb. 8, 9:
Abb. 10:
36
Finite E1emente Modell der Rohkarosserie
Finite Elemente Modell der Luft
Finite Elemente Modell der absorbierenden Flachen
Berechnete Karosserie-Ubertragungsfunktion bei offener Karosserie im Vergleich zur Messung
Schwingungsmode im Luftraum
Berechnete Ubertragungsfunktionen der geschlossenen Karosserie im Vergleich zur Messung
Berechneter schalldruck im Fahrzeuginnenraum im Vergleich zur Messung
Beteiligung der Karosseriezonen am Schalldruck
Abb. 11, 12: Einf1uE der Sitze auf den Schallpegel, berechnet und gemessen
Abb. 13, 14: Berechneter Schalldruck im Fahrzeuginnenraum mit Sitzen im vergleich zur Messung
Abb. 15, 16: Einf1uE der Schwermatten auf den Schal1pegel, berechnet und gemessen
Abb 17, 18: Berechneter Schal1druck im Fahrzeuginnenraum mit Schwermatten im Verg1eich zur Messung
Abb. 19, 20: EinfluE der gesamten inneren Ausstattung auf den Schallpegel, berechnet und gemessen
Abb. 21, 22: Berechneter Schalldruck im Fahrzeuginnenraum mit vollstandiger Innenausstattung im vergleich zur Messung
Abb. 23:
Abb. 24:
Gemessener Schalldruck im Fahrzeuginnenraum mit vollstandiger Innenausstattung an der Position Fahrerohr/Beifahrerohr
Benutzte Software und Rechnerumgebung
Generell bezeichnet in der Legende "NASTRAN" die berechnete Kurve und "TU/SFE" oder "vw" die entsprechende Messung.
FEM
-M
OD
ELL
DE
R G
OLF
RO
HK
AR
OS
SE
RIE
FUE
R A
KU
ST
IK R
EC
HN
UN
GE
N
KO
PP
LUN
G S
TR
UK
TU
R/L
UF
T
•
DA
EMPF
UN
G D
UR
CH
AB
SOR
BE
R
AN
RE
GU
NG
AM
VO
RD
ER
EN
LA
EN
GS
TR
BIL
D 1
W
-..J
38
tn W :E ::::» c( z a:: w z Cl w Z z ::::» z z - :c I LL. 0
N ...I W 0 a:: Cl
e ..J
~ -III tn -W l-e tn
::::» ...I ~ ...I < w Q a:: 0 w :E ::::»
LL. I
:E w LL.
39
II. ..I 0
" Z 0
W
W " -0 Z :l
..I Z
..I :c W 0 0 W 0 a: :e ~ a: -w I-
ew')
m 0
0
a: :l 0 ~
..I -0 c(
m
m a: c( w I :l
II. :E W I&.
40
110 '" UEBERTRAfiUNfiSnlliKTIOI'! OCR orFEN[N ROHI(AR05S (A2) GR WNIl. 45074 M[SSf'UNK TNR. 17 KO I1l"'OI1[NT[ Y
:r: u ,fJ llJOiIII -~ :c
:r: ~OOII ~ ~ --- ~
~U a< cr: ~ = ...J
, ~ Vl
, f\ '\ Ii A I. ,
ru ~ - ~OO VI .... 1::
I HI. if' lor--r\ ,,-" l/ \: V
I , 11.' \ot". - ~ -, , \ , := .... lOll
~ 58
J , " I
f I
=, '-" -, W ...J :r: 1a u lfJ
r- .... :1 , ,
-I
W
"" ...J 9.
I
i I ,
I--~ I,' :x
ill II <'4 (6 <8 36 3l 3' 36 lB 4H « •• <6 48 ~~
FRE QUENZ (HZ) LOGRR]THM]SCH
' IB ~ UEBERTRAGUNGSFUNKT[ON OER OFFENE ROHKAROSSE (A2) GR[ONR. 45074 MES5PUNKTNR, 17 KOMPONEN rE Y
:r: <-J U"l ~S~ ~--~~~r---r--+--+-----------------~ :r: f-
a< a: <.!> (:)
...J I~~~ __ ~--+H~+--;~~Il---------------~
z: .... 1118
<.!> z: 58 ~
::2 z: ~ W ...J :r: U U"l 18 w CO
...J
~ a:
58 68 78 88 'leI
FREQUENZ (HZ) LOGRR[THMISCH
BILD 4
VW/SPB
VW/SPB
41
J I UEBERTRAGUNGSrUNK nml DER GESCHLOSSErIEN ROHKAROSSE (A2) GR [ONR. '15fl71 SSPUN TNR. 17 OtflONENTE y
~ 1~;----7.~~~~-'~-r~~~~~~~~~~ ~
~ ~~ +-~~-r--~~~--+--r~~-r~R-. ::c
z:: ..... I~
<.::J z:: ~a ::::J <.::J
z:: ::::J L...J -' 18 ::t: lJ V'> L...J CD
-' 'f: a:
C\l lI! ~a S8 6e 1'6 88 ge 100
FREQUENZ (HZ) LO GAR [THM[SCH
BILO 6 JI0ut BERTRAGUNGSrUNKT [ON DER GESCHLOSSENEN ROHKAROSSE (A2)
::t: lJ V'>
~1et18 ..... 8IlII
~ 600 <.::J ~ 100
:CIla V'> ..... ::c ~
100 z:: 88 ..... 61!
<.::J z:: 18 ::::J 8 z:: 3 C\l ....J ::t: lJ V'> ....... 18 CD
GR [ONR . 45375 I"'(SSPUNK TNR. 23 KOtflONENTE Z I'
~ ~
I :~ I . , I n ,I I ~
.I ~~ ~ ,I
1\ i~ Ii II
/ I; ~I
~ '" , ; ... I ~ j
II III II , !: I
-v-J III .. pi
lI! 18 5e 61! 711 88 ge 100
FREQUE Z (HZ) LOGAR[THM[SCH
BILO 7
VW/SPB
VW/SPE
42
SCHALLPEGEL ESCHLOSSENE ROHKAROSSE (A2) AM f AHR[ROHR
I I I I .~
, ~ ::=~
, I,
~ , I • , I 11 II " ~ I~ i I , i ~ . \~ ~ ~ III' II
, ~ )I rt I. ~~ , " I '.
il ~ ~, I! • I~ " r.: , I
" " " ~~ i" , \I I I
I
,I I
" ,
L' , , , YW/SFE
~8 lB <6 58 66 18 811 'Ie 1811 116 I~ 136 146 156 16\1 118 1811 I!/II l8II
FREaU ENZ (HZ) LINERR
BILD 8 SCHALLPEGEL GESCHLOSSE N[ ROHKAROSSE (A2 ) AM B([fAHR[ROHR
65
~ G6 a> <:> ~ 5S -' <.oJ
"" 56 ..., n.. -' -' 4S a: :::t: L.> <.n 46
J8
,Iul :~i I I
1--- --- """'"' ,\ ~ .
IIY<Jl
I ~ ~ 1
~ ~ .. I ,: ' . ~ I~~ ~I Ia. II I I • I.e.
; Jl ,~, I I ;.
1;V ~ iAl , I II
I , ~, , 't{ • , , , " I
, " • . , • r~ ~' I I , ,
, I I: U , I 1,' , I' J '~ , II
- " -
,I
~ :: • f • I , YW/SFE
1 i 16 16 ~e 56 66 18 8\1 'Ie 1811 116 I~ 136 146 156 168 118 1811 1'Ie l8II
FREaUENZ (HZ) LINERR
BILD 9
43
SCW1LU'L l L C;(S(H OSS[N fWIII«WOSS[ (A2) AM BELrAHR[ROHR
: 6
i ..;
&3
&6
&l
~tK ~(.8
~ ~8 ..:-....... ~6 (L
....J S~ -' ~ ~? ;,: Sit
' 6
)8
_ Jltlll[ '..!:1lL S 11 Z( N
J I 1--r~<vo~
-_. nJ'''~H 1
1-~ nJt91.1
I -I I I I
- - I II I it t-I- I l~
'. ~ ~ n II: • ! ~
~ • - f--I I .
f : )I
til I
<D 40 ~o 66 7ll 011 9Il tOO ltO ti'll 136 148 156 tOO 1711 tEll 19I1 i'8H
FREQU[NZ (HZ) LINEAR
BILD 11 MESSUNG S(I-lALLf'[G( L C;[ SCHLOSS(N[ ROHKAROSS[ (A2) AM un rAHR[ROHR
.Qf.!!'!.f Pll ~; !1 zeN
-' ~JS ~~+-+P+-+-~~~~~~-4-4-4-+-+-+~ a:
VW/SPB
VW/SPB
<n 3iI 4~ SO 66 1D Iro 911 tllD 11 0 Ii'll t36 ''0 156 tOO 1711 Ilro 196
FREQU(NZ ( HZ) LIN[RR
BILD 12 RECHNUNG
<'S
44
SCHtlLLf'EGEl GESCHLOSS[r![ ROHI(AROSS( (fl2) flM rflHREROHR
>llT SITZ[N
! ~ :~ •
I ~ld !'O
~. I • II ••• I'IlS1liJrI
Ii:.::~ ,
", I
1 " ,I 'I ~ n~
III 1-lr ~~ I ~~ I ~ I
, ~ !
ltIt' ~ iT: t· • f r-~ III
~ " , JI :. I I I
" II
~ , I , I II
[ I, ,
, I " i ~ 'I • i [ ~
" • , II
~ , • I ,
(Il )d ~8 so 66 7& SO 98 106 llU lZ11 130 148 I~O 166 170 ISO 198 ?OO
FREOU[NZ (HZ) LIN[RR
BILD 13
VW/SFB
SCHnlU'(GEL {;(SClIl USS(tI[ RQlII(tlR05S( (fl<') AM llElrfl~IR(R(lHR
MIT S!TZ[N
;'8
I I J fr-" , It/£/U , ... -- I\W" II II ,I I
j' lJ ' I , I
I ~ I ,
: 'A ':' :t~ r:, 1 I' ~' 1 -.1'
~ I~ , ': , , , ~ I'" t" ,I 1 '" , I ~ ~n~ ~' , : ,I"
" I .. I , ~I I I,
, " 1 1
: ~ , • ~ I I • I , I
i III , I I Jr- ,
. ~ t
VW/SFB 3S
i'll 11l ~o ~o 6!l 70 au 98 1011 110 1(1l no 140 ISO 16\1 170 ISO 198
FREOUENZ (HZ) LIN[RR
BILD 14
45
SCllfllLf'(G L G[SCHLOSS[NL fWHI<l-HWSS[ (A2) AM rAIIR[ROW~
OHrIL~1T VW <;CHWLRMflTl[N
i'J -
GS -r r~ _I ~
'" JlI ~
-, 66 ri-c:<> ~, [' .::> , -
N" -J U
lj .: w I -=> II 2t' sa -.J I -.J
1 cr 3 04~ ,n I I I
-.J '8 ~ a:
:;!,
1 c8 16 'H ~6 6D 7lI ~u 96 1(16 1 m '<0 I)D '48 1~6 166 'Ill 'b'll '96 ,1m
FREQU NZ (HZ ) LI N[RR
BILD 15 MESSUNG SCI1ALLf'[G[L (;[SCll105S[N[ ROHI(IIR055 (A2) W1 rAHR[RQIIR
_ j)tJ11I 'l1l.l. VW SCHW[RI1AH[N
1 1~ • III Wil'f.l
I I
I ~~ ---It>IN.1
•• -=-~ • , AI I ..
~ ~r\
" J
, I i I ll. I~ ~ ll~ -1l
11111 ~II ~~ ,I I II I
I " ~ " , .' " II ,
II
~ is ~ N • ~ ~ III ro 9t! 'OOI~I~IW~I~ 1~lrol~ Iweoo
FREQU[NZ (HZ) LI NERR
BILD 16 RECHNUNG
VW/SFB
VW/SFB
-' L.J '-=' L.J 59 a... -' -' ~ 'S '-' V1
...,; "0 9.. a:
Y"
46
S(HnlLP[GEL G SCIILUSS(N R(lHI(flROSS[ (A2) AM rAHR(ROHR
MIT UW-SCHW[RMATTEN I I ~ _Ll " l{~
" I ~ ! ••• rc:6JtI;ft , I -f- nwr , ,
~-I
'r~' ~ , \ , I~ I
II" ~ I: I , , ) ~ " , , , . '.' : ~ " • • I • : 1 I: " ", II ~ ,
• '.1 , "
, " , , It . , . , ,r 11 " , , , , AI
~ ;jl , ~ ,
~ ~ ffl W ~ n ~ ~ l00I~I~I~lffllWl~lnl~l~lOO
FREGUENZ (HZ) LINEAR
BILD 17
VW/SFB
SCHA'LF'[ L G(SOI OSS(N( ROIII(mmSS[ (A2) AM 1l[lrAHR[ROIiR
-' ~ SO L.J a... -' ;:t . ~ :z: J; -'
fi 35
111 T UW-SCHW[RMATT[N
III J J ~ , LLIll«
" f- ----- II>Y[
~ " ~", jl , At: I [I It
IJrri ,,: , n Vi' \ I ~' , ,
if, I~: IH ~ . ~.I I:~ I " , I! , • lit I ,I ~! ir I
I I" " • , ' I ~ • , , .,
: : I I !~ I, , , I.
,: ' , I , I
"' • l1 I I I VW/SFB
:a Jll .0 S6 68 70 8U ~ 100 110 ,,0 IlO 148 lS6 I ~ In I~ IW cOO
FREGUENZ (HZ) LI NEAR
BILD 18
47
~,CHALU'[G[L RUSCi KL[ ID ,/(;ESC'HL. KAROSSE Ul2J AM rAHR[ROHR
rf~' nl TY'-~12 A J In
~ J!ll
'1 ! \ it ~ "
J i 1\ r1' ~ JII I I!
~ , Ji,l ~
~ ii ~ Ii 'N . z I
Ii r I~ rl ~ , V ., I I iii •
• ~ Ii : ; ~ ..
I~ ~I M I I
I
)5 VW/SFB
,8 Je ~6 ~6 6il 70 W 90 101'1 110 1<'0 l JO 1 4~ l~a 1611 17a 188 190 l Oll
FREQUENZ (HZ) LINERR
BILD 19 MESSUNG SCHRLU'[GEL AUS(;[KL ID./(;[SCHL . KAROS5[ (At') AM rAHR[ROHR
VW/SFB
<'0 Je ~6 S6 66 70 86 \16 101'1 110 1<'0 lJll 140 \50 1611 17a 188 \\16 ~OI'I
FREQUENZ (HZ) LINERR
BILD 20 RECHNUNG
OItiEJllWT ~~
AU88TATnJNG
OtflE/WT ----- -~tEN
AU88TATTUNG
,....,. 'l(' .:0 O~
--' ~? '-'-' so ....., c.. ~8 --' ~ <6
3 H ..n .~
~ "8
~ lS
Jj; ,. "J2
48
'.CIIHLLf'[ [l nUSG( KL[ IOl H Krll~(lSS( (flt') Al1 rflHU(ROHU
rv .J J 1\ - +~..:." . ~ \I" ~l ~nwr 'I , -j
- ;If l- I, - I r J I' II ~ \ ';1 ~ tr; ; ~ I I I" I I n j~ I , ] II I'l
.: ti: ~ , ~~ I I ~1V \" , r ,I L It ' :1~ 1r ~ T' ". I ': "
, " I
J I _i I ~eJ , I;, I ;
~! .~r r .. II! I :L ,+ "1 .. , ' 'I I I ~ ,
I I ~ I, 'I~ I Ii 'I I~ I , f
I 1\1
~ q ~ w ro w ~ loolm .~ I W lq.~~.m'~lw ~oo
FREQUENZ (HZ) LIN AR
BILD 21 ~,CHALLf'( (L AUSG(KL 10 T( KHROSS( (flo?) AM BE rrAHR ROfIR
65
--' .8 UJ <.:> ~3S --' --' a; ::I: ~ ,n
--' i'S
!j!. l O cr::
.0
~-
" ~~ u- ~ :~ I ~\ I
~ .. .. _- ~ .~ I '" ~ ~r I I - nwr ca. I I I , Ii! lL: nt ~ N , ,
I !JI ~ : ,.~ I\J , ~ I I 'J ,. r:f I' ~ I
I r ~ ~ , I • I ~ I I 1
I J . ~ I ~I : I: I 1 f
\I I
" I I"
II
: • I
: I
,9 )9 '0 50 611 m ~ W '00 110 ,;:0 130 ' .0 '56 1611 1m 186 '911 ,00 FREQUENZ (HZ) LINEAR
BILD 22
VW/SFE
VW/SFB
49
" ~ t-t-ee III
t-c:. It) ......
UJ ~ UJ >-:;:) a: ~
ee 5 g;:
z ~ w a: ~ " z W ~ ~ z I - "- ~
• 5 "':::.::: UJ .... < '-
! c U)
a. u. :.;:. .. ~ W '""";. ~= C
:Ii - "" ~: :Ii W lii 0 III .... =~ ~ ...... ~ ..... ~ 0. II 0 z:
t-...... 1i:= C") --' 0: '" Sl ~ - ..., N
:& ~
~ "" a: 1
;:: ..... C
W Wi ....J ~
a: 1 '--' .... - lj a:: :< -:z:' 5: III
W <I -en -11.: v en '-"
0 IE :;; <8 -, ~ ;: !;' :"l ~ -,
a:: (80 l 13~3dll~H)S · 1~~~
ee ~ = 0 a::
BE
NU
TZTE
SO
FT
WA
RE
UN
O R
EC
HN
ER
DA
TA
PU
NC
H
BA
SE
PL
OT
1.2
GB
M
OD
ELL
TR
AP
LO
IBM
I
I CR
AY
H
CD
C
I I
~; B
ETP
LO
(]1
30
90
X
-MP
14
CY
BE
R
0
MO
VIE
N
A.S
TRA
N
99
0
DA
TEN
I
I (A
NA
LYS
E)
PO
ST
L
••• .1
DR
UP
LO
VA
X 8
82
0 H
CYBER
D
RU
BE
T
PA
TR
AN
91
0 (P
RE
) -
"lO
AN
(A
NIM
AT
)
BIL
D 2
4
Der zukiinftige Parallelrechner von der Superworkstation zum Supercomputer
Wolfgang K. Giloi GMD-TUB
Forschungszentrum fUr innovative Rechnersysteme und -technologie Berlin
Zusammenfassung
Hochparallele MIMD-Architekturen bestehen aus einer groBen Zahl von kooperierenden Knotenrechnern. Durch Verwendung einer neuen Generation von Mikroprozessoren wird der einzelne Knotenrechner eine Leistung von einigen hundert MFLOPS haben. Das hierbei noch zu losende Problem ist das der Entwicklung von hochoptimierenden Compilern, die es ermoglichen, das enorme Leistungspotential dieser skalaren Pipeline-Prozessoren auch wirklich zunutzen. Werden so1che Architekturen als Systerne mit verteiltem Speicher ausgefUhrt, so werden sie in weiten Grenzen skalierbar, d.h. es lassen sich Parallelrechner mit einheitlicher Hardware und Systemsoftware im Leistungsbereich von einigen hundert MFLOPS (Superworkstation) bis zu einigen hundert GFLOPS (Supercomputer) realisieren. Fiir spezielle, komplexe Operationen der dreidimensionalen Graphik und der Bildanalyse wird es sich weiterhin lohnen, die universellen Knotenrechner des Systems durch·Speziaiknoten fUr diese Operationen (dies sind spezielle Vektorprozessoren oder systolische Arrays) zu ergiinzen. Die photorealistische Erzeugung dreidimensionaler Bilder und die dreidimensionale Bildanalyse konnen auf so1chen Superworkstations in Echtzeit durchgeftihrt werden.
1. Einleitung
In den letzten Jahren hat sich ein inflationarer Gebrauch des Worts "Supercomputer" breit gemacht. Jeder Rechner, der die Leistung erreicht, die die 1976 auf den Markt gebrachte CRAY-1 hatte, wird gerne als Supercomputer apostrophiert. Die kommende Generation von auf einem Chip untergebrachten 64-bit-Mikroprozessoren wie z.B. der Intel NIl wird aber auch diese Leistungsklasse erreichen, so daB nach dieser Definition in absehbarer Zukunft auch der einfache Einprozessor-Arbeitsplatzrechner ein Supercomputer ware.
Wir wollen daher den Gebrauch der Bezeichnung "Supercomputer" denjenigen Rechnersystemen vorbehalten, die das nach dem jeweiligen Stand der Technik mogliche H5chstrnaB an Leistung auf dem Anwendungsgebiet, fUr das sie konzipiert sind, erbringen. 1m Sinne dieser Definition kann man von "Superworkstations" fUr CAD oder CAE sprechen, die im Leistungsbereich von einigen 100 MFLOPS bis zu einigen GFLOPS liegen werden, dabei aber noch unter oder sogar auf dem Schreibtisch Platz finden, bzw. von "Supercomputern" fiir nichtnumerische oder numerische Anwendungen, deren Leistung gegenwartig im Bereich von einigen GFLOPS, in etwa 2-3 Jahren bei 10-100 GFLOPS, und zu Beginn des dritten Jahrtausends bei etwa einem TFLOPS liegen wird.
Bis vor wenigen Jahren war die Supercomputer-Architektur ausschlieBlich die der Vektormaschine und die Supercomputer-Technologie die der Hochstgeschwindigkeits-Schaltkreise (ECL). Nur wenige Rechnerarchitekten (wozu sich der Verfasser rechnen darf) hatten bisher die Weitsicht, zu erkennen, daB man mit hochparallelen MIMD-Architekturen nicht nur eine ebenso hohe oder auf die Dauer sogar hohere Leistung als mit den Vektormaschinen wiirde erzielen konnen, sondern durch Verwendung der hochintegrierten, wenig Leistung verbrauchenden und damit sehr viel billigeren Schaltkreis- und Aufbautechnik dieses auch zu einem wesentlich giinstigeren Preis. Die Masse der "Ungliiubigen" konnte den Preisvorteil zwar nicht bestreiten, war aber bisher der Meinung, daB die hochparallelen Systeme wesentlich schwieriger zu programmieren und letztlich aufgrund des Kommunikationsaufwands vielleicht auch weniger effizient als die einfacheren Pipeline-Vektorprozessoren waren.
52
Erst in allerjlingster Zeit beginnt sich der Durchbruch der hochparallelen MIMD-Architekturen abzuzeichnen. Ein Indiz hierflir ist die ACM Supercomputing Conference yom November 1989, wo namhaf,teste Fachleute (Rechnerarchitekten, Compiler- und Anwendungsexperten) der USA in einer dreistlindigen Panel Discussion das Thema "Teraflops Machines - Myth or Reality" diskutierten und dabei ohne Ausnahme die Sicht auGerten, daG in den nachsten Jahren die "killer micros" (wie es ein Panelist ausdrtickte) die Vektormaschinen zu verdrangen beginnen werden. Als Beispiele flir besonders richtungsweisende Entwicklungen in dieser Richtung wurden die Connection Machine [1], die iPSCRechner [1] und der SUPRENUM-Rechner [1], [21 genannt. Besonders illustriert wurde diese Entwicklung aber durch die Andeutungen des Vorstandsvorsitzenden der Firma Cray Research, ebenfalls in dieser Richtung tatig zu werden [3].
Dieser Beitrag versucht aufzuzeigen, we1che Entwicklungen auf diesem Gebiet bis zur Mitte der neunziger Jahre zu erwarten sind und we1che Problemkomplexe erst noch befriedigend ge16st werden mlissen. Dabei wird nicht nur auf Allzweck-Supercomputer flir numerische Anwendungen eingegangen (auf hochparallele nicht-numerische Rechner einzugehen zu wollen ware noch zu frlih), sondern auch auf Spezialrechner, z.B. flir photorealistische dreidimensionale Bilddarstellungen oder Echtzeit-Bilderkennung.
2. Technologie-Entwicklungsprognosen
In jlingster Zeit ist eine revolutionare Entwicklung auf dem Gebiet der Mikroprozessoren eingetreten, die sich insbesondere in einer exponentiellen Zunahme derr MIPS-Leistung auGert. Nach Joy gilt flir die Leistungszunahme seit 1984 die Formel
MIPS = 2YEIIR 1984
Dieses Joysche Gesetz wird durch Bild 1 veranschaulicht.
P/MIPS
128 Y-1984
P =2
64
32
16
8
4
2
1 Y '84 '85 '86 '87 '88 '89 '90 '91
Bild 1. Joysches Gesetz
Das Joysche Gesetz laGt den Eindruck entstehen, als handle es sich bei dieser Leistungssteigerung urn einen stetigen ProzeG. Dem ist jedoch nicht so. Das Joysche Gesetz - beliebig extrapoliert - erweckt die Vorstellung, daG 1994 die Leistung des einzelnen Mikroprozessors 1 GIPS erreicht. Wenn diese
53
exponentielle Leistungszunahme iiber die Jahre gel ten wiirde, dann brauchte man sich schlieBlich nicht mehr mit Parallelarchitekturen zu beschiiftigen. Spezielle Hardware fUr Funktionen der photorealistischen Bilderzeugung oder der Bildanalyse wiirde sich ebenfalls nicht mehr lohnen. In Wirklichkeit wird die in Bild 1 gezeigte Kurve aber irgendwann abknicken. Wir wollen hier untersuchen, wann dieser Knickpunkt zu erwarten ist. Dazu miissen wir betrachten, wo die starke Leistungssteigerung der letzten Jahre herriihrt. Wir werden weiterhin untersuchen, wie iiber die Leistung des einzelnen Prozessors hinaus weitere Leistungssteigerung erzielt werden kann.
Bild 2 zeigt die Zunahme der Taktfrequenz von Mikroprozessoren iiber die Jahre 1987-1993. Die Kurve zeigt eine Geschwindigkeitsverdopplung etwa aile drei Jahre an, die durch eine Halbierung der Kanalliingen erreicht wird. 1) Die treibende Kraft ist hierbei die Entwicklung bei den dynamischen Speicherbausteinen (DRAM). Dabei ist bemerkenswert, daB die DRAM-Bausteine selbst aufgrund gegenliiufiger Effekte nicht im gleichen MaBe schneller werden.
MHz=MIPS MFLOPS (2 operations chained)
100 200
80
60
1 00
40
20
10 ' 87 ' 88 ' 89 ' 90 '91 '92 ' 93
Bild 2. Zunahme der Taktfrequenz von Prozessoren 1987 - 1993
Tabelle 1 gibt die vorhergesagte Weiterentwicklung der DRAM-Bausteine an. Die Zugriffszeit der DRAM wird kaum mehr wesentlich unter den jetzt erreichten Bestwert von etwa 25 ns fUr den Zugriff im static column mode bzw. 60 ns fUr den wahlfreien Zugriff (die Zykluszeit betriigt mehr als das doppelte) gesenkt werden k6nnen. Man erkennt, daB die Diskrepanz zwischen Prozessor-Geschwindigkeit und DRAM-Geschwindigkeit immer gr6Ber wird, was Riickwirkungen auf die Prozessorarchitektur hat (Notwendigkeit von Cache-Speichem). Giinstiger sind die Verhiiltnisse bei statischen Speicherbausteinen (SRAM), die sich grundsiitzlich jeweils mit einem Viertel der Kapazitiit des entsprechenden SRAM-Bausteins realisieren lassen2), wobei die Zugriffszeit im Maximalfalle bei CMOS-SRAM inzwischen nahe bei 10 und bei BiCMOS-SRAM deutlich unter 10 Nanosekunden liegt.
Eine Geschwindigkeitsverdopplung bei den Prozessoren aile 3 Jahre entspricht einer Steigerung von 26 % pro Jahr. Der Unterschied zum Faktor 2, wie ihn das Joysche Gesetz angibt, kann somit nur durch architektonische MaBnahmen erzieit worden sein. Wie sie aile wissen, ist es die RISC-Bewegung, die diesen ProzeB in Gang gebracht hat.
1) Die Halbierung der Kanallange halbiert den Strom, reduziert aber die der Flache proportionale Kapazitlit urn den Faktor 4. Dadurch werden Kapazitliten doppelt so schnell umgeladen.
2) Fiir eine DRAM-Zelle werden 2 Transistoren benbtigt (einer als Ladungsspeicher, der andere zur AnkQPplung). SRAM-Zellen sind Flipflops, fiir die typischerweise 6 Transistoren aufgewendet werden.
54
Jahr Dimension (11) Kapazitiit (Mbit) Zahl der Transistoren
1983 2.5 ..... 2.0 0.25 0.5 Millionen 1986 1.25 ... 1.0 1 2 Millionen 1989 0.8 ..... 0.7 4 8 Millionen 1992 0.5 16 32 Millionen 1995 0.35 64 128 Millionen
Tabelle 1. Entwicklung der DRAM-Technologie 1983-1995
3. Der Pipeline-Prozessor
Das Ziel der RISC-Prozessor-Entwicklung war es, durch extreme Einfachheit der Befehle, Uberlappung der Transporte von Befehlen und Daten und der Befehlsausflihrung und fest verdrahteter Steuerung anstelle der Mikroprogrammierung zu Prozessoren zu gelangen, die mit jedem Takt einen Befehl ausflihren konnen. Dies mag flir Arbeitsplatzrechner, die nur der Programmentwicklung dienen, akzeptabel sein, flir Arbeitsplatzrechner im CAD- und CAE-Bereich benotigt man jedoch nach wie vor komplexere Funktionen wie zum Beispiel die Gleitpunktoperationen, die nicht innerhalb eines Taktzyklus von 20 ... 25 Nanosekunden ausgeflihrt werden konnen, zumindest nicht in CMOS-Technologie.
Da aber bereits das Holen und Interpretieren der Befehle, das Holen der Daten und die Ausflihrung der Befehle im Pipelineverfahren iiberlappt ausgeflihrt werden, hindert einen niemand daran, die Ausflihrung komplexerer Befehle ebenfaJls in das pipelining einzubeziehen. Dies hat zur Entwicklung der Pipeline-Skalarprozessoren geflihrt.
Das Pipeline-Verfahren ist nun nicht gerade neu. Es ist diese Methode, die den Vektorrechnem seit vielen Jahren ihren starken Leistungsvorsprung gegentiber konventionellen "Super Mainframes", die von der Hardware her etwa gleich aufwendig sind, verleiht. Neu ist die Anwendung dieses Prinzips auch im skalaren Prozessor, und dies wurde erst durch die enormen Fortschritte moglich, die man in den letzten Jahren beztiglich der Entwicklung hochoptimierender Compiler gemacht hat.
Was der Pipeline-SkaJarprozessor mit dem Vektorprozessor gemein hat, ist die Tatsache, daB beide Prozessoren im gtinstigsten Faile mit jedem Taktzyklus eine vollstandige Operation ausflihren kann. Beim Vorhandensein entsprechender Hardware schlieBt dies Gleitpunktoperationen ein. Hierdurch ist mehr als durch jede andere MaBnahme die plotzliche enonne Leistungssteigerung erzielt worden, wie sie durch die Aufstellung in Tabelle 2 zum Ausdruck gebracht wird3) [4]. Wie Tabelle 2 zeigt, steht zu erwarten, daB der skalare Pipeline-Prozessor in wenigen Jahren den in gleicher Technologie ausgeflihrten Vektorprozessor leistungsmaBig eingeholt haben wird - vorausgesetzt, das Compilerproblem kann gelost werden.
1m Gegensatz zum Vektorprozessor, bei der der ununterbrochene Strom von ausflihrbaren Operationen, der zur Ftillung der Pipeline benotigt wird, durch die natiirliche Ordnung der Vektorelemente und die a priori gegebene Datenunabhangigkeit der elementenweisen Ausflihrung der Operationen der Vektor-Datentypen gegeben ist, muB beim skalaren Pipeline-Prozessor ein so1cher Operationsstrom erst vom Compiler erzeugt werden, und zwar auf der Grundlage einer aufwendigen Datenabhangigkeitsanalyse, sowie komplizierter Optimierungsalgorithmen. Urn es vorwegzunehmen: Solche Compiler gibt es im Augenblick noch nicht, aber es wird sie in absehbarer Zeit geben.
Es gibt noch einen anderen Unterschied zwischen dem Pipeline-Vektorprozessor und dem Pipeline-SkaJarprozessor. Der Vektorprozessor erhalt seine Datenstrome aus dem Hauptspeicher oder einem eigenen, hinreichend groBen Vektorspeicher. Seine Leistung ist damit letztlich durch die Speicherbandbreite begrenzt, ganz gleich, wie leistungsfahig der Pipelineprozessor ist. Beispiel: der Intel-Prozessor i860, der im gtinstigsten Falle flir einen Speicherzugriff 2 Takte benOtigt, vorausgesetzt der Speicher ist ein im static column mode betriebener DRAM oder ein schneller SRAM, so daB keine
3) Die in Tarelle 2 angegebenen Werte entstammen zum Teil entsprechenden Firmen-Druckschriften, zum Teil sind es Werte, die von Herstellern in Gesprachen mit uns genannt wurden.
55
Prozessor-Typ 1990 1991 Zahl der Chips! (40 MHz) (50 MHz) Technologie
1. "Yon Neumann Type"
TRANSPUTER T840 4 5 I PGA/CMOS
2. Skalarer Pipeline-Prozessor
MOTOROLA MC68040 6.4* 8* 1 PGA/CMOS MOTOROLA MC88000 10 16 3 PGA/CMOS INTEL i860 40/20** 50/25** 1 PGA/CMOS INTEL NIl 100/100 1 PGA/CMOS
3. Yektor-Prozessor
ADSP711On120 80* ? 3 PGAJECL WEITEK WTL3364 80* 100** 1 PGA/CMOS GENESIS (Entwurf) 120** 2 PGA/BIGMOS
Erliuterungen Spitzenleistungen in MFLOPS, IEEE double precision
PGA pin grid array * geschiitzt (Motorola gibt keine Spitzenleistung an) ** "chaining" von 2 Operationen; der zweite Wert gibt die bei der Vektor-
verarbeitung erzielbare Spitzenleistung an
Tabelle 2. Spitzenleistung der neuesten Generation von Mikroprozessoren
zusatzlichen Wartezustande entstehen, wird bei der Verarbeitung langerer Vektoren niemals iiber 20 MFLOPS hinauskommen. 1m Gegensatz dazu fiihrt der Skalarprozessor seine Operationen so weit wie moglich auf Daten aus, die sich in seinen internen Registern oder dem on-chip cache befinden, und es ist Sache des Compilers, die im Gegensatz zur Vektorverarbeitung nun gegebenen Datenabhangigkeiten zwischen den Operationen dazu zu nutzen, moglichst viele Zwischenergebnisse intern zu speichern und damit die Anzahl der externen Speicherzugriffe zu minimieren. Das bedeutet aber auch, daB die angegebene Spitzenleistung nur im Burst-Modus allenfalls erreicht wird (Merke: Die Spitzenleistung ist die Leistung, von der der Hersteller garantiert, daB sie niemals iiberschritten wird).
4. Das Compiler-Problem
Die enormen Spitzenleistungs-Angaben - in der Regel so viele MIPS oder MFLOPS pro Funktionseinheit des Prozessors wie der Takt in MHz betragt - der neuesten Generation von Mikroprozessoren wie den Intel-Prozessoren i486 oder i860 haben zunachst eine Euphorie hervorgerufen, die sehr bald von einer ebenso wenig gerecht fertigten negativen Beurteilung abgelOst wurde, sobald echte Leistungsdaten wie z.B. LINPACK-Benchmarkwerte bekannt wurden. So hat der i860, der von Intel als 80 MFLOPSProzessor, bzw. (unter Zusarnrnenzahlung aller "Apfel und Birnen") als 1 20-MIPS-Prozessor eingefiihrt wurde, eine LINPACK-Leistung von 7 MFLOPS mit dem Green Hills Fortran Compiler (der C Compiler liefert etwa den gleichen Wert) oder 13 MFLOPS handcodiert [5]. Der Leistungsunterschied urn fast den Faktor 2 zeigt, wie sehr es noch beim Compiler hapert. In der Tat fiihrt der Green Hills Compiler alle heute iiblichen Optirnierungen wie colored register allocation und loop unrolling and com-paction' aus, kann jedoch weder den Pipelinebetrieb des Prozessors noch seinen dual instruction -,tlOde
56
noch seinen dual operation mode ausnutzen. Einen Compiler, der das konnen wird, wird es erst in etwa 2 bis 3 Jahren geben.
Was der i860 in bestimmten Fallen wirklich leisten kann, zeigt zum Beispiel Tabelle 3. Wir sehen hier ein (handcodiertes) Programm fUr die innerste Schleife eines Ray-Tracing-Algorithmus4). Durch Ausnutzung des dual instruction modes sowie dem vierfachen loop unrolling kann der gesamte AIgorithmus mit nur 8 Befehlen - d.h. in nur 200 Nanosekunden - ausgefUhrt werden. Man sieht, daB dabei die Pipeline sehr gut gefUllt werden kann. Bis es einen Compiler gibt, der einen solch optimalen Code zu erzeugen in der Lage ist, werden, wie gesagt, aber noch mindestens zwei Jahre vergehen.
Nun kann man, wie dies im obigen Beispiel geschehen ist, fUr Spezialanwendungen wie 3D-Graphik oder Bildanalyse die zeitkritischsten Algorithmen durchaus von Hand codieren, da diese in der Regel aus nur kurzen und damit Uberschaubaren Prozeduren bestehen. FUr komplexe Anwendungsprogramme gilt dies narurlich nicht, hier wird man auf die noch ausstehenden leistungsfahigen Compiler warten mUssen.
5. Der Multi-PE-Prozessor
Was Hillt sich in Zukunft tun, urn das Joysche Gesetz noch weiter erfUllen zu konnen? Mit dem Takt ist man bei CMOS mit 50 MHz am Ende angelangt. BiCMOS wird es erlauben, noch einmal. einen Faktor 4 an Geschwindigkeit zu gewinnen (Uber einen Zeitraum von 6 Jahren, d.h. bis etwa 1995). Damit sind etwa urn 1995 Mikroprozessoren mit einer Spitzenleistung von 200 MIPS bzw. 200 MFLOPS fUr einzelne Operationen und 400 MFLOPS fUr zwei gleichzeitig ausgefUhrte GleitpunktOperationen ("chaining") zu erwarten.
Weitere Leistungssteigerungen konnen daruber hinaus nur noch durch weitere Vervielfachung der Funktionseinheiten auf dem Chip erzielt werden. Intel hat mit dem i960 bereits einen Prozessor fUr "eingebettete Systeme" auf dem Markt, der zwei lnteger-Einheiten besitzt und daher als 66-MIPS-Prozessor (bei 33 MHz Taktfrequenz) gehandeit wird. Motorola wird mit dem 88110 das gleiche tun.
Die Vervielfachung von Funktionseinheiten verscharft natUrlich die Compilerproblematik noch weiter stark. Es ist ja bereits jetzt schon gerade der Dual Instruction Mode (AusfUhrung einer Integerund einer Gleitpunkt-Operation gleichzeitig) bzw. der Dual Operation Mode (AusfUhrung von 2 Gleitpunkt-Operationen gleichzeitig), die beim i860 das groEte Problem fUr den Compiler darstellt. Dies ist plausibel, wenn man bedenkt, daB jetzt Compiler jetzt nicht nur einen optimierten Operationsstrom erzeugen muB, sondern mehrere, aufeinander synchronisierte Strome.
Das gleiche Problem stellt sich auch bisher schon bei der Very Long Instruction Word (VUW) Machine [6], die ebenfalls durch das Vorhandensein mehrerer Funktionseinheiten gekennzeichnet ist. Der Compiler fUr eine VLIW -Maschine "kompaktisiert" den ursprunglichen, sequentiellen Code durch Ausnutzung der Parallelarbeitsmoglichkeiten mehrerer Funktionseinheiten und beschleunigt dadurch die ProgrammausfUhrung. Die Code-Kompaktisierung ist an sich ein Verfahren, welches man in der Mikrocode-Optimierung seit vielen Jahren anwendet. Entscheidend ist aber jetzt, daB man nicht wie beim Mikrocode nur basic blocks (das sind Code-Teile, die auBer am Eingang oder Ausgang keine Sprungbefehle aufweisen) betrachtet, sondern lange, mit Sprungbefehlen durchsetzte Code-Sequenzen. Dadurch ergibt sich ein hoherer Grad an nutzbarer Parallelitat. Ein Verfahren, das dies ermoglicht, ist unter dem Namen Trace Scheduling entwickelt worden [6]. Die Anzahl der Funktionseinheiten in VLIW-Maschinen muE dennoch aber klein bleiben (typisch 4, hochstens 8), urn noch sich lohnende Ausnutzungsgrade zu erhalten.
Auf diese Weise die von Joy fUr 1994 vorhergesagten Leistung von 1 GIPS wirklich zu erreichen (z.E. durch 8 parallel arbeitende Funktionseinheiten bei 120 MHz Takt) erscheint angesichts dieser Betrachtungen unrealistisch. Aber auch wenn man realistischerweise annimmt, daB die Joysche Kurve etwa bei einem halben GFLOPS abknicken wird, wird man immer noch einen Prozessor mit sehr hoher Leistung als Baustein fUr Multiprozessor-Architekturen zur VerfUgung haben, durch die man, wie es
4) Das Programm wurde von Michael Riisse geschrieben.
57
FB UNIT CORE UNIT ADDER PIPELINE 1 2 3 R
LOOP - fix(n) fix(e) -d.pfadd e de fe ; be DONE e+de - fix(n) fixe(e) d.pfadd n dn fn ; shl #2 in in n+dn e+de - fix(n) d.pfsub zr dz -- ; fld terrain(in) zl zr-dz n+dn e+de -d.fxfr fe -- ii ; nop zr-dz n+dn c+de d.pfix e -- e ; nop fix(c) zr-dz n+dn c+de d.fxfr fn -- in ; sht #10 ii ii fix(e) zr-dz n+dn d.pfix n -- n ; nop fix(n) fix(c) zr-dz n+dn d.pfgt zr zt zr ; add - fix(n) fix(e) zr-dz
d.pfadd e de fe ; be DONE e+dc - fix(n) fix(e) d.pfadd n dn fn ; shl #2 in in n+dn e+de - fix(n) d.pfsub zr dz -- ; fld telTain(in) zi zr-dz n+dn e+dc -
d.fxfr fe -- ii ; nop zr-dz n+dn e+de d.pfix c -- e ; nop fix(e) zr-dz n+dn c+dc d.fxfr fn -- in ; shl #10 ii ii fix(c) zr-dz n+dn d.pfix n -- n ; nop fix(n) fix(e) zr-dz n+dn d.pfgt zr zt zr ; add - fix(n) fix(c) zr-dz
d.pfadd e de fe ; be DONE e+de - fix(n) fix(e) d.pfadd n dn fn ; shl #2 in in n+dn c+dc - fix(n) d.pfsub zr dz -- ; fld tcrrain(in) zl zr-dz n+dn e+de -
d.fxfr fe -- ii ; nop zr-dz n+dn e+dc d.pfix c -- e ; nop fix(c) zr-dz n+dn e+de d.fxfr fn -- in ; shl #10 ii ii fix(c) zr-dz n+dn d.pfix n -- n ; nop fix(n) fix(e) zr-dz n+dn d.pfgt zr zt zr ; add - fix(n) fix(e) zr-dz
d.pfadd e de fe ; be DONE e+de - fix(n) fix(e) d.pfadd n dn fn ; shl #2 in in n+dn e+de - fix(n) d.pfsub zr dz -- ; fld terrain(in) zl zr-dz n+dn e+de -
d.fxfr fe -- ii ; nop zr-dz n+dn e+de d.pfix e -- c ; nop fix(e) zr-dz n+dn e+de d.fxfr fn -- in ; shl #10 ii ii fix(e) zr-dz n+dn d.pfix n -- n ; nop fix(n) fix(e) zr-dz n+dn d.pfgt zr zt zr ; add - fix(n) fix(e) zr-dz
DONE
Tabelle 3. i860-Programm fiir die innerste Schleife eines Ray Tracing-Algorithmus
inzwischen zum Stand der Technik geworden ist, die Leistung des einzelnen Prozessors multiplizieren kann.
6. Multiprozessor-Chips
1m Gegensatz zu den hochregularen Speicherstrukturen sind Prozessor-Layouts stark ilTegular. Vor aI-lem wird beim Prozessor ein erheblicher Teil der Chipflache fUr die internen Verbindungsstrukturen (Busse) benotigt. Dadurch laBt sich bei gleicher Technologie hochstens die Halfte der Transistoren des entsprechenden DRAM-Bausteins realisieren. Urn hier ein Musterbeispiel zunennen: den 2 Millionen Transistoren des I-Mbit-DRAM stehen beim i860 knapp 1 Million Transistoren gegen liber. Der i860 zeigt, daB man damit noch nicht einen voll befriedigenden Prozessor realisieren kann, sondern Kompromisse eingehen muB, die die Leistung beeintrachtigen (keine snoop logic zur Gewahrleisnmg .der cache coherence, keine Prozesskennung im Cache, die Caches sind zu klein, die doppeltcgenaue Mul-
58
tiplikation benotigt 2 Takte). AuBerdem scheinen auBer den Erfindern des Transputers die ProzessorArchitekten noch nicht erkannt zu haben, daB man bei einem Prozessor, der als Baustein flir Multiprozessor-Architekturen gedacht ist, spezieUe Hardwareeinrichtungen flir die Kommunikation mit anderen Prozessoren vorsehen muB, wenn man eine effiziente Kommunikation erzielen will [4], [7].
Mit der inzwischen existierenden Technologie des 4-Mbit-DRAMs lassen sich aber aUe diese Nachteile beheben. Solche Prozessoren werden nicht zuletzt mit Registerfiles, Instruktions- und Datencaches von angemessener Kapazitiit, sowie speziellen Funktionseinheiten z.B. flir Graphik oder Signalverarbeitung ausgestattet sein. Mit der ab 1992 zu erwanenden Technologie des 16-Mbit-DRAM wird man dann 4 solcher Prozessoren auf einem Chip unterbringen konnen. Da dies im wesentlichen nur die Replikation des Layouts des einzelnen Prozessors bedeutet, fiiUt der Entwurfsaufwand hierflir nicht sehr ins Gewicht. Entsprechende Gehiiuseformen mit einer hinreichend groBen Zahl von Anschliissen wird es in der Form der TAB-Bausteine (tape bounded devices) geben.
Wir sehen hierin einen Weg, den Fortschritt der GroBtintegration zu nutzen, ohne wie beim Multi-PE-Prozessor das Compilerproblem immer weiter zu komplizieren. Tatsiichlich verlagert sich das Software-Problem, und zwar vom Compiler zum Knoten-Betriebssystem, welches zwangsHiufig nun auch ein Multitasking-Betriebssystem sein muB, urn die Prozessoren geeignet einsetzen zu konnen. Entwicklungen wie MACH [8] oder PEACE [9] sind aber der Beweis daflir, daB sich solche Multitasking-Betriebssysteme sehr effizient implementieren lassen.
Mit solchen Multiprozessor-Bausteinen lassen sich hochparallele Systeme aufbauen, bei denen mehrere Knotenrechner mit ihren privaten, viele MByte groBen Speichern auf einer Platine untergebracht sind. Voraussetzung hierflir ist, daB diese Bausteine wirklich konsequent als Komponenten flir Multiprozessor-Architekturen konzipiert werden, wozu noch erhebliche Forschungsarbeit mit dem Ziel der Entwicklung hocheffizienter Kommunikationsmechanismen zu leisten ist.
Den kompletten Rechenknoten einschlieBlich seines Speichers auf einem Chip wird es hingegen nicht geben konnen, da man heute Speicherkapazitiiten forden, die ein Vielfaches der Kapazitiit eines einzelnen DRAM-Bausteins betragen.
7. Parallelrechner-Architekturen
Die Leistung im GFLOPS-Bereich, die zum Beispiel flir Echtzeit-Anwendungen auf den Gebieten der Simulation, der photorealistischen dreidimensionalen GraphikdarsteUungen oder der Bildanalyse benotigt werden, wird man auch in absehbarer Zeit nur durch die ParaUelarbeit einer entsprechenden Anzahl von Prozessoren erhalten, d.h. durch MIMD-Parallelrechnerarchitekturen.
Es gibt zwei Kategorien von MIMD-ParaUelrechnerarchitekturen, die speichergekoppelten Multi-prozessorsysteme (shared memory systems) und die aus einer Vielzahl von Knotenrechnern bestehenden Systeme mit verteiltem Speicher (distributed memory systems). Da beim System mit verteiltem Speicher die Kommunikation nicht tiber gemeinsam benutzte, allen Knoten zugiingliche Speicherbereiche erfolgen kann (da es nur "private" Speicher gibt), sondem nur durch den Austausch von Botschaften, wird dieses auch als botschaften-orientiertes System (message-passing system) bezeichnet.
Beide Kategorien unterscheiden sich wesentlich in der Granularitiit der ParaUelarbeit und damit auch in dem Programmiermodell flir den Benutzer.
In den speichergekoppelten Multiprozessorsystemen findet ParaUelarbeit auf der Ebene der einzelnen Anweisungen einer Programmiersprache mit skalaren Datentypen (in der Regel ist dies eine "von Neumann"-Sprache wie Fortran, C, usw.) statt (fine grain parallelism). Dies geschieht vorwiegend dadurch, daB die Iterationen einer Schleife gleichzeitig ausgeflihrt werden (DO ALL-Konstrukt). Berechnung wird im Sinne des Von Neumann-Modells als schrittweise Veriinderung des Speicherzustands der Maschine verstanden, wobei hier im Gegensatz zurn sequentieUen Rechner der Zustand einer Anzahl von Speicherobjekten gleichzeitig veriindert werden kann.
In den botschaften-orientierten Parallelrechnem findet die ParaUelarbeit hingegen in der Form konkurrenter, kommunizierender Prozesse statt (coarse grain parallelism). Dies entspricht dem ProgrammiermodeU des objektorientierten Programmierens.
59
Bisher sind die speichergekoppelten Multiprozessorsysteme dominierend. HierfUr gibt es zwei Griinde:
1. Ihr Programmiermodell entfemt sich weniger yom Gewohnten. 2. Die Kommunikation iiber einen gemeinsam zugiinglichen Speicher steHt die einfachste und effi-
zienteste Uisung dar, solange die Zahl der Prozessoren klein ist (4, 8 oder h5chstens 16).
Bei hochparallelen Systemen mit einer groBen Zahl von Knotenrechnem liiBt sich die Speicherkopplung jedoch nicht verniinftig anwenden, und zwar wegen der iiberhand nehmenden Zugriffskonflikte, die sehr schnell einen untragbaren "System-Flaschenhals" erzeugen konnen. Ein weiterer Vorteil der verteilten Systeme ist, daB sie im Gegensatz zu den speichergekoppelten Systemen in wei ten Grenzen skalierbar sind (s. unten). SchlieBlich ist der coarse grain parallelism prinzipiell effizienter als der fine grain parallelism, da zwischen den einzelnen Kommunikationsvorgiingen ein wesentlich hoheres MaB an Rechenarbeit geleistet wird.
Bei den botschaften-orientierten Systemen ist man dafUr mit einem potentiell starken Flaschenhals anderer An konfrontiert, der durch die Verlustzeiten gegeben ist, wiihrend der die Knotenrechner auf eine Botschaft warten miissen und dabei nichts anderes tun konnen. Diese Latenzzeit fUr einen Botschaften-Austausch wird durch die Summe zweier Werte gegeben, die startup time und die transmission time.
Die Transmission Time ist die fiir die Ubertragung der Botschaft selbst benotigte Zeit. Sie kann durch ein sehr breitbandiges Verbindungsnetzwerk und die Verwendung moglichst effizienter Verbindungsprotokolle beim heutigen Stand der Technik sehr gering gemacht werden; wie gering, ist letztlich eine Frage des Aufwands. Es liiBt sich zeigen, daB ein gleichzeitiges Optimum von Leistung und Kosteneffektivitat durch eine zweistufige Verbindungshierarchie erreicht wird, d.h. durch eine Struktur, in der auf der unteren Ebene eine gewisse Zahl von Knotenrechnern jeweils zu einem Cluster verbunden werden und auf der oberen Ebene dann die so geformten Cluster untereinander [4], [10].
'pie Startup Time ist die Zeit, die fiir die vorbereitenden Schritte benotigt wird, bis eine Botschaften-Ubertragung iiberhaupt erst begonnen werden kann. Sie hiingt sehr stark von der Art des fiir die Inter-ProzeB-Communication (!PC) verwendet en Protokolls ab (und damit letztlich davon, welche Sicherheitsanforderungen fiir den Botschaften-Austausch gestellt werden). Fiir den Fall einer yom Betriebssystem durchgefiihrten Kommunikation auf der Grundlage eines gesicherten Rendezvous-Protokolls liegt die Startup Time typischerweise bei etwa 0.5 - 2 ms. Durch eine extrem effiziente Implementierung, wie sie im PEACE-Betriebssystem vorgenommen wurde [9] (Teams von leichtgewichtigen Prozessen), liiBt sich dieser Wert auf 100 - 300 Ils reduzieren. Durch Aufwand eines speziellen Kommunikationsprozessors kann man die Startup-Zeit urn eine weitere Zehnerpotenz verringem [7]. Damit steht dem Aufbau sehr effizienter, hochparalleler MIMD-Architekturen mit verteiltem Speicher nichts mehr im Wege.
Will man auch beim System mit verteiltem Speicher am gewohnteren Programmierstil des gemeinsam benutzten Speichers (s. oben) festhalten, so erfordert dies, daB man auf das botschaftenorientiene System ein virtual shared memory aufsetzt. Wie man dies moglichst effizient und gleichzeitig kosteneffektiv durchfiihrt, ist ebenfalls noch ein Forschungsthema. Man kann aber vorweg sagen, daB auch hier nur eine geeignete Hardware-Unterstiitzung (z.B. wieder in der Form des dediziert en Kommunikationsprozessors) zu einer tragbaren Effizienz fiihren kann.
8. LohDt sich Doch die VektorverarbeituDg?
Bis vor etwa einem Jahr war die Antwort auf diese Frage ein eindeutiges Ja. Bevor 1988 die ersten RISC-Prozessoren mit Gleitpunkt-Hardware wie z.B. der Motorola MC88000 verfiigbar wurden, gab es als billige Uisung nur die G1eitpunkt-Coprozessoren von der Art der Motorola MC88881/88882 oder der Intel i80287/80387) mit einer Leistung in der GroBenordnung von 100 KFLOPS. Durch Hinzufiigen eines Vektorprozessors (z.B. mit den Bausteinen WTL2264/2265 oder WTL3364 von Weitek oder den Bausteinen ADSP7110n 120 von Analog Devices) konnte man die skalare Gleitpunkt -Leistung urn etwa
60
eine GroBenordnung und die Pipeline-Vektorleistung urn iiber 2 GroBenordnungen steigern. Gegenwiirtig ist bei einem Knotenrechner mit dem Intel i860 als CPU durch Hinzufligen eines Vektorprozessors (z.? WTL 3364) immer noch eine Steigerung der Vektorverarbeitungsleistung etwa urn den Faktor 3 zu erzielen. Da die dadurch verursachte Kostensteigerung des Knotenrechners nur etwa 40-50 % des Preises eines skalaren Knotens betragt, ist die "Vektorisierung" bei Rechnern flir numerische Anwendungen immer noch ein gutes Geschaft.
Mit dem etwa in einem Jahr verfiigbar gewordenen i860-Nachfolger (NIl) wird aber zumindest nominell die gleiche Spitzenleistung erreicht werden, die die dann verftigbaren Vektorprozessor-Bausteine haben werden (s. Tabelle 2). Wie der skalare Pipeline-Prozessor im Faile der Vektorverarbeitung gegeniiber einem speziellen Vektorprozessor beziiglich der echt erzielbaren Leistung abschneiden wird, wirft letztlich wieder die bereits diskutierte Compilerfrage auf. In den nachsten 2-3 Jahren zumindest ist immer noch ein Leistungsunterschied von immer noch einem Faktor 2 zu erwarten. Dieser Leistungsunterschied kann natiirlich aber auch durch Verdoppelung der Knotenzahl ausgeglichen werden (s. den Abschnitt "Multiprozessor-Chips"), und es wird noch sorgfaltiger Untersuchungen bediirfen, hier die optimale LOsung zu finden.
Eindeutiger zu beantworten ist die obige Fragestellung, wenn es urn spezielle Operationen der Graphik oder der Signal- und Bildverarbeitung geht. Tabelle 4 liefert einen Vergleich der Arbeitsgeschwindigkeit des Intel i860 mit der Arbeitsgeschwindigkeit spezieller Pipeline-Prozessoren flir reprasentative Operationen der Graphik und Bildverarbeitung.
Algorithm INTEL i860 Pipeline Processor (40 MHz) (30 MHz)
Double precision vector operations 20 MFLOPS 60 MFLOPS Ray tracing 5 million incremental steps ---Phong shading 0.2 million pixels 20 million pixels 4x4 matrix multiplication I million transfonnations 5 million transfonnations 3x3 convolution over small win- 20 million operations 60 million operations dow 3x3 convulution of unwindowed 10 million operations 340 million operations image (systolic array chip) Fast fourier transfonn 0.7 million complex values 5 million complex values
(TRW chip set)
Tabelle 4. Leistung des i860 im Vergleich zu dedizierten Prozessoren
Dieser Leistungsvergleich zeigt, daB es sich in jedem Faile lohnt, in einer hochparallelen MIMDArchitektur spezielle Knoten fiir spezielle Aufgaben vorzusehen. Wendet man das oben beschriebene Prinzip der Clusterung durch Aufbau einer mehrstufig-hierarchischen Verbindungsstruktur an, so wird man solche speziellen Knoten in die Cluster integrieren. Bild 3 zeigt dies am Beispiel eines speziellen Knotens fiir die Erzeugung dreidimensionaler Bilder mit schattierten Oberflachen in Echtzeit, wie er in die Cluster des in der Entwicklung befindlichen GENESIS-Supercomputers [4], integriert werden solI.
9. Skalierbarkeit der Architektur
Es wurde bereits erwahnt, daB botschaften-orientierte Systeme in wei ten Grenzen skalierbar sind. Das heiBt, man kann die Zahl der Knotenrechner in weiten Grenzen wahlen, ohne daB dies Riickwirkungen auf die Systemsoftware hat. Damit sind bei standardisierter Hardware und Systemsoftware zum Beispiel die in Tabelle 5 angeflihrten Konfigurationen moglich.
Jedes dieser Systeme kann durch Spezialknoten flir Signal- und Bildverarbeitung, Graphik und Netzwerkanschliisse (Ethernet, FDDI, ... ) erganzt werden. Platten speicher lassen sich sowohl auf der Ebene des einzelnen Knotenrechners (billige Festplatten (Winchester) mit SCSI-Interface) wieaticn auf
SEGMENT MEMORY
SEGMENT l ~vnrf'''E:I VERTEX GROUP COORD.
FilE FilE FilE
61
32 CROSSBAR SWITCH
CLUSTER OF NODES
20 SUB-SYSTEM
RENDERING PROCESSOR
RENDERING SUBSYSTEM
FRAME BUFFER SYSTEM
Z·BUFFER MEMORY
Bild 3. Integration eines Graphikknotens in ein Cluster [4]
System-Typ Leistungsbereich Aufbau
Super-Arbeitsplatzrechner 0.1 ..• 1 GFLOPS 1 •.. 4 Knotenrechner 1 Platine, Tischgehause
Super-Workstation 1 ... 10 GFLOPS 1 Cluster von (bis zu) 32 Knoten Supercomputer hochster Lei- 100 .•• 1000 GFLOPS bis zu 32 Cluster, entsprechend stung viele Schriinkc
Tabelle S. Konfiguration verschiedener Systeme aus gleicher Hardware und Software
der Clusterebene (Hochgeschwindigkeitsplatten mit eigenem Controller) in das System integrieren [4J. Als vorgeschalteter Rechner, der die Benutzeroberflache bereitstellt, dient ein Standard-Arbeitsplatz-
62
rechner mit seiner universellen Software-Ausstattung und Programmentwicklungs-Umgebung (UNIX mit X-Windows und entsprechenden Crosscompilern). (1989), ACM order no. 4 15891,313-321
10. Literaturverzeichnis
[1] Gonauser M., Mrva M.(ed.): Multiprozessor-Systeme, Springer-Verlag 1989
[2] Giloi W.K.: The SUPRENUM Architecture, Proc. CONPAR 88, Part A, British Computer Society 1988, 1-8
[3] Rollwagen J.: Keynote Speech at Supercomputing '89, ACM Internal. Conf., Reno, Nevade, Nov. 1989
[4] Gi10i W.K.(ed.): GENESIS - The Architecture and Its Rationale, Intem. Report, ESPRIT-II Project P2702, June 1989
[5] Anonymus: i860 Processor Performance, Release 1.0 (March 1989). Intel Corporation
[6] Colwell R.P. et al.: A VLfW Architecture for a Trace Scheduling Compiler, Proc. 2nd IntemaL Conf. on Architectural Support for Programming Languages and Operating Systems, IEEE Computer Society Press, Oct. 1987, 180-192
[7] Giloi W.K., Schroeder W.: Very High-Speed Communication in Large MIMD Super computers, Proc. 3rd. Internal. Conf. on Supercomputing
[8] Young M. et al.: The Duality of Memory and Communication in the Implementation Multiprocessor Ope-rating System, of a ACM Operating System Review 21.5, (Proe. of the Eleventh ACM Symposium on Operating Systems Principles, Austin, Texas 1987)
[9] Schroeder W.: The PEACE Operating System and Its Suitability for MIMD Message Passing Systems, Proe. CONPAR 88, Part A, British Computer Society 1988, 17-24
[l0] Giloi W.K., Montenegro S.: Super Interconnection Structures for Super Computers, Proe. ICS 89, Internal. Supercomputing Institute 1989,310-316
Einsatz von Supercomputern fUr die Crashberechnung
J. Hillmann
PKW-Berechnung
Volkswagen AG Wolfsburg
Wurden Crashsimulationen von Gesamtfahrzeugen noch vor wenigen Jahren als
'numerisches Abenteuer ' bewertet (1), so hat sich dieser Zweig der Berechnung
inzwischen so weit entwickelt, daB er heute als zielfindend eingestuft wird (2).
In diesem Beitrag wird der erreichte Stand der Crashberechnungen bei Volkswa-
gen an einigen Beispielen aufgezeigt. Neben der sUindigen Weiterentwicklung der
Finite-Element Verfahren kommt der grafischen Aufbereitung der Berechnungser-
gebnisse (Postprozessing) eine wachsende Bedeutung zu (3), denn sie fUhrt zu
einer verbesserten Dialogfahigkeit zwischen den Berechnungsingenieuren und den
Mitarbeitern in den Konstruktions- und Versuchsabteilungen. AbschlieBend werden
aus der Sicht des Berechnungsingenieurs zukunftige Anforderungen an die Hard-
und Software fur die Crashsimulation formuliert.
EINLEITUNG
Yom Beginn der Entwicklungsgeschichte eines neuen Fahrzeuges bis zum Serien-
anlauf mussen immer wieder Zielkonflikte gel6st werden. Um die optimale 'Crash '-
L6sung zu finden, wurden bisher fast ausschlieBlich Crashversuche an Komponen-
ten, Aggregatetragern und spater an Prototypen durchgefUhrt. Mit der VerfUgbar-
keit von Supercomputern und numerisch stabilen expliziten dynamischen Finite-Ele-
ment Verfahren sind Berechnungsingenieure seit einiger Zeit in der Lage, prog-
nosefahige Aussagen zum Crashverhalten komplexer Strukturen zu liefern, bevor
diese gebaut werden sind (4).
64
Der Wert der Berechnungsunterstutzung wird wesentlich durch den zeitkritischen
Ablauf der Entwicklung bestimmt. In der Vorentwicklungsphase werden ausgehend
yom Finite-Element-Modell (FE-Modell) des Vorgiingerfahrzeuges umfangreiche Grund-
satzuntersuchungen in enger Zusammenarbeit mit den Konstrukteuren durchge-
fUhrt. Die gesammelten Erkenntnisse fliel3en in das neue Fahrzeug ein. Zum Ende
der Konstruktionsphase mul3 in maglichst kurzer Zeit aus den CAD-Daten und
evtl. weiteren Zeichnungen das FE-Modell der I. und spiiter II. Baustufe gene-
riert werden. Bei Volkswagen ist fUr diesen Prozel3 ein in dem CAD-System
VWSURF integrierter Netzgenerator VWMESH entstanden. Ober das Konzept und
den erreichten Entwicklungsstand berichten U. Sorgatz und H. Deuter in (5). Je
schneller das Modell erstellt werden kann, urn so fruher kannen Fehlstellen auf-
gespurt und korrigiert werden.
2 CRASHBERECHNUNGEN
I m Jahre 1985 wurden die ersten Crashberechnungen an Fahrzeugkarosserien durch-
gefUhrt (vgl. (6)). Seit 1987 werden bei Volkswagen die Programme DYNA3D (8)
(Prof. Hallquist/Prof. Schweizerhof) und PAMCRASH (9) (Haug u.a., Fa. ESI)
aufeiner CRAY X-MP/14 verwendet (vgl. (10)).
Bis heute konnte durch die enge Zusammenarbeit mit den Softwareerstellern die
GENAUIGKEIT und die
ZUVERLASSIGKEIT gesteigert, das
ANWENDUNGSSPEKTRUM erweitert und die
RECHENZEITEN erheblich verkurzt werden.
Die gleichzeitige Verfugbarkeit beider Programme erlaubt eine absolute Beurteilung
des Leistungsvermagens bezuglich dieser Kriterien.
2.1 BERECHNUNGSVERFAHREN
Die Gleichgewichtsbedingung fUr den linearen dynamischen Fall lautet:
mit den bekannten GraBen und Pt als Lastvektor.
65
1m linearen Fall ohne Berucksichtigung von Diimpfungsanteilen k6nnen die inneren
Kriifte fi als Produkt aus Steifigkeitsmatrix und Verschiebungsvektor berechnet
werden:
Bei Berucksichtigung von Nichtlinearitiiten werden die inneren Krafte berechnet
aus:
fi S B TS dV. V
Da die Massenmatrix M nur auf der Hauptdiagonalen besetzt ist, sind aile Glei-
chungen entkoppelt:
Mit dem zentralen Differenzenverfahren werden die Geschwindigkeiten und Ver-
schiebungen zum Zeitpunkt t + At ermittel t (expl izites Verfahren).
Die Einfachheit des Verfahrens wird durch sehr kleine Zeitschritte 'erkauft'. Die
numerische Stabilitiit des expliziten Verfahrens ist nur dann hinreichend gesi-
chert, wenn fUr den kleinsten Zeitschritt min dt die Bedingung eingehalten wird:
mit
min A t ~
min A I
E
e.
min .0.1
{Eie
kleinste Elementkantenliinge im gesamten FE-Netz;
Elastizitiitsmodul;
Dichte des verwendeten Materials.
Fur Nutzanwendungen im Automobilbau liegt der kleinste Zeitschritt zwischen 2 I-ls
und 0,5 I-ls.
66
In DYNA3D und PAMCRASH stehen dem Anwender Balken-, Schalen-, SchweiB-
punkt- bzw. Nietelemente, unterschiedliche Werkstoffmodelle, verschiedene Kon-
taktdefinitionen usw. zur VerfUgung.
Kontaktbereiche konnen als Master-slave-contact oder Single-surface-contact de-finiert werden. Mit steigender Zahl von Kontaktdefinitionen steigt auch der Be-rechnungsaufwand deutlich an. Bis zu 75 % der CPU-Zeit kann von den Kontaktalgo-
rithmen benotigt werden.
2.2 BERECHNUNGSAUFWAND
Der Berechl;lungsaufwand steigt linear mit der Anzahl der verwendeten Elemente an. Er kann wie folgt abgeschatzt werden:
CPU-Sek. Ncycle x Nelem x Tcycle
mit
Ncycle - Anzahl der Zeitschritte fUr die Gesamtberechnungszeit T: T : Ncycle T/(1 ,5*4t);
Nelem - Anzahl der Elemente und Tcycle - Time per zone cycle
Der 'Time per zone cycle', d. h. die Berechnungszeit fUr ein Element ist stark abhangig vom verwendeten Rechner und von Art und Umfang der Kontaktberei-
che.
Fur die CRA Y X-MP 114 gilt:
) 20 \..IS - ohne Kontakt und elastisches Materialverhalten;
" 90 \..IS - 50 % Kontakt und plastisches Materialverhalten.
Fur die Cyber 990 VF liegen die Zeiten etwa doppelt, fur eine VAX 11-780 etwa 500 - 600 mal so hoch.
67
Zur Zeit wird bei VW mit speziellen Fahrzeugmodellen fUr den Front-, Seiten- und
Heckcrash gearbeitet. Jedes FE-Modell umfaBt etwa 15 000 bis 25 000 Elemente.
Fur etwa 40 - 50 % der Elemente sind Kontaktbedingungen definiert. Der mittlere
Zeitschritt liegt in der GroBenordnung At = 1 ~s. Wird die Crashsimulation mit
T = 100 ms festgelegt, so ergibt sich fUr die CRAY X-MP/14 folgende CPU-Zeit:
CPU-Sek. = 100 000 cycle * 20 000 Elemente * 90 ~s
180 000 Sekunden = 50 Stunden.
Eine VAX 11-780 ware also mit der Losung dieses Problems ca. 3 - 3 1/2 Jahre
beschilftigt. Doch 50 Stunden Rechenzeit fUr einen solchen Job auf der CRAY
konnen - auch an einem Wochenende - kaum zur Verfugung gestellt werden. Hier
muB also die Gesamtberechnungszeit gesenkt und es muB mit Hestarts gearbeitet
werden.
2.3 PRE- UND POSTPROZESSING
Fur die Geometrie- und Netzerstellung wird neben PAT RAN immer starker das
CAD-System VWSURF mit VWMESH eingesetzt. Die Einzelbauteile werden schlieBlich
mit PAT RAN zum FE-Gesamtmodell zusammengefUgt. Hier werden auch Materialeigen-
schaften, Kontaktbereiche usw. definiert. Mit Hilfe des VW eigenen Programms
TRECAN werden dann die Inputdecks fUr DYNA3D und PAMCRASH erstellt.
Fur die Auswertung der meist sehr umfangreichen Ergebnisdaten (Postprozessing-
Phase) stehen eine ganze Reihe weiterer Programmsysteme zur VerfUgung. Mit
TAURUS, MOVIE und PAT RAN konnen Verformungsbilder als Drahtmodell oder
'fringe plots' - auch mit verdeckten Linien - dargestellt werden. Einzelne Kur-
venverlaufe konnen mit den UN I RAS-Routinen prasentationsfahig aufbereitet wer-
den. Das Programm PASIR wurde von VW speziell fur die Animation von FE-Struk-
turen unter Ausnutzung der lokalen Fahigkeiten des PS390 von Evans and Suther-
land entwickelt. Die Verformungen fUr einzelne Zeitpunkte werden als Drahtmo-
delle im richtigen zeitlichen Ablauf nacheinander dargestellt. Schon mittels we-
niger Bilder kann der Betrachter den Ablauf der wichtigsten Ereignisse verfolgen
und bewerten. Falsche Randbedingungen, Durchdringungen, Modellierungsfehler
usw. werden erkannt.
68
Ober einen Converter kann ein Videorecorder direkt angeschlossen werden, so
daB visualisierte Crashberechnungsergebnisse auf diese Weise sehr kostengunstig
gespeichert werden k6nnen. Fur die I nsassensimulation stehen weitere Programme
zur Verfugung (vgl. 11).
Bild 1 gibt einen Oberblick uber die bei VW verwendete Hard- und Software fUr
Crashberechnungen.
3 CRASHBERECHNUNGSERGEBNISSE
3.1 LANGSTRAGERBERECHNUNG
Crashanalysen sind zu einem wichtigen Werkzeug in der Entwicklungsgeschichte
eines neuen Fahrzeuges geworden. Die Einbindung in die Sicherheitsentwicklung
ist in (4) ausfUhrlich dargestellt. Die dort angegebenen Beispiele werden mit den
aktuellsten Programmversionen fur DYNA3D und PAMCRASH (hier: Stand Februar
1989) erneut berechnet, dargestellt und verglichen, dabei greifen wir fUr die
Beurteilung von Verbesserungen in den Programmpaketen im Sinne eines Benchmark-
tests immer wieder auf die 'alten' Beispiele zuruck.
Bild 2 zeigt das FE-Modell eines Uingstdigerabschnitts, der im vorderen Bereich
durch den Kuhlertrager in Querrichtung abgestutzt ist. Die Anbindungspunkte am
Langstrager sind zusatzl ich verstarkt.
Die Abschatzung nach Beermann (12) fUr eine statische Belastung ergibt fur den
vorderen Bereich ei ne mittlere Faltenbeulkraft von 62,0 kN. Berucksichtigt man
den EinfluB der Geschwindigkeit mit dem Faktor Ct. = 1,2, so ergibt sich eine
mittlere Faltenbeulkraft von 74,4 kN. Um das rechnerisch ermittelte Energiepo-
tential des Tragers voll auszusch6pfen, wird fUr den Versuch mit einer 1100 kg
Barriere die Aufprallgeschwindigkeit mit Vo = 24,6 km/h festgelegt.
Das FE-Modell besteht aus 2928 Schalenelementen. Die Kontaktbereiche sind in
sechs Zonen mit insgesamt 1520 Elementen (Single-surface-contact) aufgeteilt.
69
Wie in der Realitat werden auch im FE-Modell die Flansche des Langstragers und
des Deckbleches iirtlich durch S<;:hweiBpunktelemente miteinander verknupft. Diese
Bereiche werden in die Kontaktdefinitionen mit aufgenommen. Aus Symmetriegrun-
den wird nur eine Fahrzeugseite berechnet.
Es werden die Ergebnisse der Programmversionen von August 1988 (V0888) und
Februar 1989 (V0289) verglichen. Die Berechnungsergebnisse bis 120 ms sind fUr
PAMCRASH in Tafel 1 und fUr DYNA3D in Tafel 2 zusammengestellt.
Dabei bedeutet:
volle Ite.
drei Ite,
ohne Ite.
Iterative ErfUllung der FlieBbedingung in jedem Zeit-
schritt;
- Annaherung an die FlieBbedingung mit drei Iterations-
schritten;
Keine Iteration im Werkstoffgesetz.
Mit max. U x wird die griiBte Verschiebung in x-Richtung bezeichnet.
In PAMCRASH wird das isoparametrische Schalenelement nach Belytschko/Linl
Tsay (13) verwendet. In DYNA3D ist auch ein Schalenelement nach Hughesl Liu
(14) implementiert. Damit werden groBe Dehnungen und groBe Rotationen korrekt
erfaBt.
Version Plastizitat Ncycle CPU-Sek. I max u
V0888 ohne Ite. 137 537 40 500 271 m; ---
--~~~::--------1--~:~:-:~:::---1-~~~--~~-~-----1-~~--:~~------1--~-~-::------ I Tafel 1: Zusammenstellung der Berechnungsergebnisse fUr PAMCRASH
Vergleicht man die max. Verschiebungen fUr PAMCRASH, so ergeben sich fUr
die neue Version griiBere Verschiebungen als bisher (s. Tafel 1). Hier wur-
den die Default-Werte der Hourglass-Koeffizienten geandert und es wurden
einige Anderungen in den Kontaktalgorithmen vorgenommen. Die Rechenge-
schwindigkeit konnte um ca. 13 % gesteigert werden.
70
~:::~on :,:m::: :~',~t;,::~t I ~::::4 I C:'";::k. :;::: __________________ 9 _________________________________________________________
V0289 Belytschko ohne Ite. 124 928 24 276 294 mm V0289 Belytschko volle Ite. 138 206 27 360 280 mm
V0289 Hughes volle Ite. 131 045 35 550 285 mm
Tafel 2: Zusammenstellung der Berechnungsergebnisse fUr DYNA3D
Die Berechnungsergebnisse fUr DYNA3D zeigen, daB die max. Verschiebungen
in der neuen Version kleiner sind als in der alten Version (5. Tafel 2). Hier
wurden die Kontaktalgorithmen weiter verbessert. Die Unterschiede in den
Verschiebungen fUr die Schalenelemente nach Hughes und Belytschko sind
recht gering, es ergeben sich aber doch deutliche Unterschiede zwischen den
verschiedenen Moglichkeiten zur ErfUliung der FlieBbedingung.
Auch in DYNA3D konnte eine deutliche Steigerung der Rechengeschwindigkeit
erreicht werden. Die Werte in Tafel 2 belegen, daB das Schalenelement nach
Hughes erheblich rechenaufwendiger ist als das Schalenelement nach Belytsch-
ko. Dabei stellt das Element nach Hughes mit voller Iteration im Werkstoffge-
setz theoretisch die bestmogliche Naherung dar.
Die berechneten Verformungen fur DYNA3D und PAMCRASH sind in den Bil-
dern 3 und 4 in sechs Phasen bis t = 120 ms nach dem Aufprall dargestellt.
Fur den Versuch wurde ein Uingstragerpaar mit dem Kuhlertrager als vordere
Querverbindung auf einer Halteplatte befestigt (5. Bild 5). Bild 6 zeigt einen
Langstragerabschnitt nach dem Versuch. Es wurde eine max. dynamische Verschie-
bung von 306 mm erreicht. Wird das kontaktgebende Element in die Auswer-
tung mit einbezogen, so ergibt sich ein etwas kleinerer Verformungsweg. Der
entsprechende Zeitversatz ist dann ebenfalls zu berucksichtigen.
In Bild 7 sind die von DYNA3D und PAMCRASH berechneten und im Versuch
gemessenen AufprallkrMte als Funktion der Zeit aufgetragen.
71
3.2 GESAMTFAHRZEUGBERECHNUNG
Fur die Erstellung des Gesamtfahrzeug-Modells wurden weitgehend Teile aus
bereits bestehenden FE-Netzen verwendet. Der Motor- und der Getriebeblock
sind uber spezielle Lager mit den tragenden Teilen verbunden. Anbauteile
wurden als Punktmassen berucksichtigt. 1m Bereich des Vorderwagens sind
insgesamt 16 Kontaktzonen (Single-surface-contact) mit 8153 Kontaktelementen
definiert.
Bi Id 8 zeigt das Gesamtfahrzeug-Modell. Es besteht aus 15267 Schalenelemen-
ten und mehr als 100 Balkenelementen. Fast 50 % der Elemente wurden fUr die
Modellierung der Uingstriiger verwendet. So kann das Faltenbeulphiinomen der
Liingstriiger auch im Gesamtfahrzeugmodell erfaGt werden.
Das Crashverhalten des Fahrzeugmodells wurde mit PAMCRASH und DYNA3D
berechnet. Die Daten sind in Tafel 3 zusammengestellt.
Version Element Plastizitiit Ncycle CPU-Sek.
DYNA3D V0289 Belytschko volle Ite. 123 302 128 840
PAMCRASH V0289 Belytschko ohne Ite. 95 258 128 100
Tafel 3: Zusammenstellung der Daten
Die verformte Struktur 90 ms nach dem Aufprall ist in Bild 9 bzw. Bild 10
dargestellt.
Vergleicht man beide Bilder, so sind nur geringe Unterschiede feststellbar.
Die mit PAMCRASH berechnete max. Verschiebung in x-Richtung weicht nur
um 12,6 mm yom DYNA3D Ergebnis abo
72
Fur eine Bewertung der Ergebnisse sollte aber nicht nur das globale Struk-
turverhalten mal3gebend sein. Es muf3 auch das lokale Verhalten einzelner
Strukturelemente und die Wechselwirkung miteinander betrachtet werden. Da-
fUr ist die Animation der Berechnungsergebnisse auf dem PS 300 von Evans
uns Sutherland mittels PAS I Rein wichtiges Werkzeug.
In Bild 11 sind die Aufprallkraft-Zeit Kurven dargestellt. Der Zeitpunkt des
Wandaufprall des Motors ist identisch. Die Kraftspitze ist in der Berechnung
jedoch immer grof3er als im Versuch, solange 'weiche' Anbauteile wie der Kuh-
ler nicht berucksichtigt werden.
Nach 30 ms sind bereits 60 % der kinetischen Energie in plastische Arbeit
umgewandelt. Den uberwiegenden Teil haben die Uingstdiger ubernommen.
Mit geringem Aufwand konnen auch die Beschleunigungs-, Geschwindigkeits-
und Verschiebungs-Zeit Kurven fUr verschiedene Knotenpunkte ermittelt wer-
den.
Die speziell gefilterten Beschleunigungs-Zeit Kurven werden dann als Einga-
bedaten fur das I nsassensimulationsprogramm MADYMO verwendet. Mit diesem
Werkzeug kann das Ruckhaltesystem optimal an das Verhalten der Fahrzeug-
struktur angepal3t werden, um so ein hohes MaB an Sicherheit fUr die Insas-
sen zu erreichen.
4. AUSBLICK
Supercomputer der nachsten Generation mit mehreren Prozessoren werden die
vielfache Leistungsfahigkeit heutiger Supercomputer haben. Diese Tatsache
wird nicht nur Auswirkungen auf die Software, sondern auch auf die GroBe
der FE-Modelle haben. Werden heute noch Modelle fUr jede spezielle Crashart
aufgebaut, um die Berechnungszeit in einer sinnvollen GroBenordnung zu hal-
ten, so wird es dann sicher nur noch ein Fahrzeugmodell geben, das fUr sehr
viele Berechnungslastfalle geeignet ist. Ein solches Gesamtfahrzeugmodell wird
dann aus mehr als 50 000 Elementen bestehen. Bis dahin mussen jedoch noch
eine ganze Reihe von 'Werkzeugen' weiterentwickelt werden, damit Modelle
solcher Gr6Benordnung sinnvoll aufgebaut, verwaltet und berechnet werden
k6nnen.
73
Die heute eingesetzte Software fUr das Postprozessing muB weiter verbessert
und erhebl ich komfortabler gestaltet werden. Besonders die Crashsimulation
kann nicht auf die Animation verzichten. Es muB moglich werden, flachenhaft
dargestellte Modelle (hidden line) in Echtzeit zu drehen. Ein solches Werkzeug
kann ganz erheblich zum Verstandnis der komplexen physikalischen Vorgange
beitragen. Dieses wird dann eine hervorragende Basis fUr die Diskussion mit
den Konstukteuren und Versuchsingenieuren sein.
Mit den heute zur VerfUgung stehenden Mitteln k6nnen jedoch bereits progno-
sefahige Aussagen zum Verhalten der Karosserie unter Crashbelastung ge-
macht werden. Die Gute der Ergebnisse wird mit feiner werdenden Netzen,
weiter verbesserten Berechnungsalgorithmen und noch leistungsfahigeren Su-
percomputern auch in Zukunft weiter steigen. Durch eine rechnerische Be-
gleitung von der Vorentwicklungsphase bis zum Serienanlauf wird es gelingen,
Fehlentwicklungen zu vermeiden und damit einen Beitrag zur Verkurzung der
Entwicklungszeiten und Senkung der Entwicklungskosten zu leisten.
74
LITERATURLISTE
(1) Fiala, E.; Sorgatz, U.: Moderne Methoden fUr die Karosserie-Ent-wicklung. Stahl und Eisen 106 (1986), Nr. 12, S. 29-35.
(2) Seiffert, U.; Scharnhorst, T.: Die Bedeutung von Berechnungen und Simulationen fur den Automobilbau. VDI-Bericht Nr. 699 (1988), S. 1 - 36.
(3) Spreng, H.P.; Pries, H.: Einsatz numerischer Verfahren bei Ent-wicklung und Optimierung des neuen VW-PASSAT und ihre Auswir-kungen. VDI-Bericht Nr. 699 (1988), S. 655-676.
(4) Hillmann, J.; Rabethge, W.: Rechnergestutzte Entwicklung einer Vorderwagenstruktur zur Verbesserung der Sicherheit von PKW-In-sassen beim Frontaufprall. Automobiltechnische Zeitschrift (1988), N r. 11, S. 641-646.
(5) Sorgatz, U.; Deuter, H.: Das System VWMESH zur Idealisierung von Tragstrukturen im CAE-Konzept. VDI-Z 131 (1989) Nr. 3, S. 26-32.
(6) Haug, E., Scharnhorst, T. and Dubois, P.: FEM-Crash, Berech-nung eines Fahrzeugfrontalaufpralls. VD I-Bericht N r. 613 (1986), S. 479-505.
(7) DYNA3D USER'S MANUAL (Version 13). Hallquist, J .0.; Benson, D.J.: (Nonlinear Dynamic Analysis of Structures in three Dimensions) Rev. Jan. 1989.
(8) Hallquist, J .0.; Schweizerhof, K.: CRASH/IMPACT-Berechnungen mit DYNA3D - Neue Entwicklungen, User-Seminar Munchen (Septem-ber 1987).
(9) PAMCRASH MANUAL (Version 10.2), Rev. 1-89/JCB Engineering System International. (1989).
(10) Hillmann, J., Konig, C., Schettler-Kohler, R. und Schneider, J.: Crashsimulation. I nterner Bericht VW-Wolfsburg (1988).
(11) Meyer-Prussner, R.: PAPAS-Animation of Passenger Simulation Results. 1. Madymo User's Meeting. TNO Delft (Sept. 1988).
(12) Beermann, H.-J.; Staisch, A.: Aufpralluntersuchungen mit verein-fachten Strukturmodellen. IfF-Tagung Braunschweig (1982), S. 205-237.
(13) Belytschko, T.; Lin, J.I.; Tsay, C.S.: Explicit Algorithms for the Nonlinear Dynamic of Shells. Compo Methods in Applied Mechanics and Engineering 42 (1984), pp. 225-251.
(14) Hughes, T.J.R.; Liu, W.K.: Nonlinear Finite Element Analysis of Shells. Three-dimensional Shells. Computational Methods in Applied Mechanics 27 (1981), pp. 331-362.
Geometrie-~ Erstellung ~'----""
"C
~ Tektronix Workstation PS 390
Preprozessing FE-Modell
Tektronix PS 390 Workstation
2 VWMESH '2 PATRAN .~ c .c. u en 4: o >
75
Crash-berechnungen
Postprozessing FE -Ergebnisse
CRAY X-MP/14 Tektronix PS 390 Workstation Video
z DYNA3D « ~ PAMCRASH 0:: I--
1---9 PATRAN a.. PASIR a: TAURUS t-=" • MOVIE ~:!4 UNIRAS z-oG ua..
MKS-Insassen-simulation
Tektronix Video
MADYMO • PASIR
PAPAS
Bild 1: Einsatz von Hard- und Software fur die Crashberechnung bei VW-Wolfsburg
Bild 2: Finite-Element-Modell des Uingstragerabschnitts
o ms
Bild 3:
76
....."
~~ ~~
~ -a;:.,
Sequenz der verformten Langstragerstruktur bis 120 ms DYNA3D - V0289, Hughes, volle Ite.
•
o ms
~
20 ms
~
40 ms
~
Bild 4:
77
~ ~ ~
~
"
~ ~ ~,GJ
.£ ;:.l.J-""
""
Sequenz der verformten Uingstragerstruktur bis 120 ms PAMCRASH - V0289, volle Ite.
78
Bild 5: Langstragerpaar vor dem Versuch
Bild 6: Langstrager nach dem Aufprall
200
110
..,
4D
2t)
0
79
.. ·· .. · VERSUCH :" ., --_. PAMCRASH
-DYNA3D
.. : : :',: t:'" : : ! 1.1 1;-
; , : ~K'" ,: \J ....... l: I r i
:W \l If' '\~ ~.
.... r- ' . ':.. .,
',\: ~' ~-'.
\ "\.
'" ~ .. I.:.. •• , o m W 00 ~ w ro ro 00 W m ~ ~
Zeit [ms]
Bild 7: Vergleich der Kraft-Zeit Kurven aus Berechnung und Versuch
B i Id 8: Finite-Element Modell des Gesamtfahrzeugs
80
Bild 9: PAMCRASH Verformte Fahrzeugstruktur 90 ms nach dem Aufprall
Bild 10: DYNA3D Verformte Fahrzeugstruktur 90 ms nach dem Aufprall
KARTOGRAPHISCH ORIENTIERTE DARSTEILUNG VON MESSWERTEN
Reinhold Millier-Meernach, Franz Josef Prester S.E.P.P. GmbH, Lohmiihlweg 4,8551 Rottenbach
Dr. Peter Weigert Bundesgesundheitsamt, Zentrale Erfassungs- und Bewertungstelle fur Umweltchemikalien (ZEBS), Werner-VoB-Damm 62, 1000 Berlin 42
1. Zusammenfassung
In einem bundesweiten Projekt werden Lebensmittelproben genommen und auf Schadstoffe untersucht. Die gemessenen Daten werden auf einem Mainframe gesammelt und verwaltet. Es wird ein Arbeitsplatz beschrieben, welcher diese Daten yom Host-Rechner holt und graphisch visualisiert. Hierbei steht die kartographische Darstellung im Vordergnmd. Die Aufgabenschwerpunkte des Systems sind erstens das interaktive Arbeiten am Bildschirm urn die erfaBten Mefiwerte schnell zu sichten und zu bewerten, zweitens die graphische Verkniipfung von MeBreihen urn auffallige Korrelationen zwischen den MeBreihen festzustellen, sowie Gebiete mit hoher oder niedriger Schadstoffkonzentration zu lokalisieren und drittens die Erstellung von veriifi"entlichungsreifem Kartenmaterial. Die Realisierung erfolgte auf einem modernen Arbeitsplatzrechner mit aktueller Graphikhardware, deren integrierter Graphikprozessor schnelle Bildaufbauzeiten ermoglicht. Standardisierte Graphiksoftware GKS und das mehrtaskfahige Betriebssystem OS/2 halten die Softwarekosten niedrig und ermoglichen den Einsatz des Arbeitsplatzes auch flir andere Aufgabenstellungen.
2. Daten
In etwa 70 Regionen der Bundesrepublik Deutschland werden zusammen etwa 1500 Lebensmittel auf 20-40 Stoffe untersucht. Dies fiihrt zu etwa 30.000 bis 40.000 MeBwerten pro Quartal, ab 1991 verdoppelt sich diese Zahl (11/, /2/). Die MeBwerte werden auf einem BS2000 Rechner tiber das Datenbanksystem ADABAS verwaltet. Uber eine Standard-PC-Kopplung werden ausgewahlte MeBwerte und Aggregationen aus ihnen auf den graphischen Arbeitsplatz gekoppelt und dort mit Hilfe der Datenbank ORACLE organisiert. Hier werden die darzustellenden MeBreihen zusammengestellt.
82
3. Regionen
Bezugsgebiet ist die Bundesrepublik Deutschland einschlieBlich West-Berlin. Es stehen aus Grenzabschnittsdateien generierte Grenzpolynomdateien fur BundesHinder, Regierungsbezirke, Landkreise und Gemeinden zur Verfiigung. Die Identifizierung erfolgt tiber den Gemeindekennziffernschltissel. Die Regionengrenzdaten und -polygone wurden von /3/ erfaBt. Darzustellende Regionen konnen aus den MeBwertdateien generiert werden, interaktiv bestimmt werden oder aus vorhergehenden Sitzungen tibernommen werden. Die Regionen konnen sowohl zusammenhangend als auch explodiert dargestellt werden, auch konnen interessante Regionen mit einer Urnrandung in einstellbarer Strichstarke versehen werden.
4. Datendarstellung
Die Darstellung basiert in erster linie auf der Relation 'Gemeindekennziffer-MeBwert'. Zur kartographischen Einfarbung ist eine Klassenbildung erforderlich. Diese kann automatisch, halbautomatisch oder manuell erfolgen. Die Klasseneinteilung wird als Vorschlag fur eine Legendengenerierung verwendet. Den Klassen wird ein Farbverlauf zugeordnet. Neben der Erstellung von veroffentlichungsreifem Kartenmaterial ist die gleichzeitige Darstellung von mehreren MeBreihen fur gleiche Regionen Aufgabenschwerpunkt des Systems. Flir diese Verknlipfung von MeBwerten stehen zwei Verfahren zur Verfiigung:
1. Normierte Addition der Klassenran~ Die addierten Klassenrange der MeBwerte fur eine Region werden beziiglich der Anzahl der MeBreihen normiert. Der Farbverlauf tiber die Klassen bleibt erhalten. Diese Darstellung ermoglicht das Erkennen von 'hoch-hoch' und 'niedrigniedrig' Korrelationen.
2. Farborientierte Darstellung;. Diese Darstellung wirkt wie das Ubereinanderlegen von Folien auf dem Overheadprojektor. Es entstehen neue Farbverlaufe. Die Anzahl der resultierenden Klassen ist das Produkt der Ausgangsklassen. Hier schranken die Hardwaremoglichkeiten die Anzahl der Ausgangsklassen ein. Diese Darstellung ermoglicht dem Anwender, die resultierende Karte zu be-urteilen und laBt Rtickschllisse auf die konstituierenden MeBwerte zu.
Durch eine Mehrebenentechnik konnen Fllisse und die Hauptverkehrswege (Schiene und StraBe) ein- und ausgeblendet werden. Dies gilt auch fur die Namen der Regionen. FUr ein sinnvolles Arbeiten in diesem Sinne sind sehr schnelle Bildaufbauzeiten erforderlich. Dies betrifft sowohl den Gesamtbildaufbau, als auch insbesondere die Darstellung der Anderung der MeBwerte beim Visualisieren von neuen oder weiteren MeJ3-reihen. Die verwendetete Hardware wird durch geeignete Softwareuntersrutzung diesen Anforderungen gerecht.
83
3. Farben
Die Farbzuordnung erfolgt tiber sogenannte Paletten. Dies sind Farbverlaufe, welche yom Anwender interaktiv bestimmt werden und unter einem Namen abgespeichert werden konnen. So kann zu einen ein einheitliches farbliches Erscheinungsbild des Kartenmaterials gesichert werden, zum anderen kann das System an unterschiedliche Farbfahigkeiten von z.B. Bildschirm und Plotter angepaf3t werden.
5. Layout
Zum Gestalten der Karten stehen clem Anwender eine Reihe von Moglichkeiten zur Verfugung: Diagramme in Balken- und Kuchenform, Notchbox- und Box-Whiskerdiagramme sowie entsprechende Legenden hierzu erweitern das Darstellungspektrum. AuBerdem konnen interaktiv freie Texte mit verschiedenen Textparametern plaziert werden. Eine Vektorgraphik ermoglicht es, im Dialog mit Maus oder Tablett Bezugspfeile zu setzen und Objekte zu umranden oder zu unterstreichen. Ein Logo mit variablem Text ermoglicht eine leichte Identifizierung der Karten. Verschiedene Plottformate sind einstellbar (GroBe und Orientierung).
6. Hardware
Das System wurde auf einen Arbeitsplatzrechner mit 80386 Prozessor und mathemathischem Co-Prozessor realisiert. Das System ist mit einem Hochleistungsgraphikboard ausgestattet, welches auf einem der modernsten integrierten Graphikkontroller, dem QPDM von AMD, basiert. Die AuflOsung betragt 1024 x 1268 Bildpunkte bei 256 Farben gleichzeitig auf einem 19" Monitor.
7. Software
Die graphische Implementierung basiert auf dem Standard system GKS. So konnten die Softwarekosten gering gehalten werden und das System ist offen fur die Zukunft. Dies betrifft sowohl die allgemeine Portierbarkeit als auch den technischen Fortschritt der Hardware. Ebenso konnen problemlos verschiedene Plotter- oder Druckermodelle eingesetzt werden. Das Betriebssystem OS/2 wurde aus mehreren Grunden eingesetzt: Es fehlen einige liistige MS-DOS Beschriinkungen wie z.E. die 640 kB Speicherplatzgrenze, die Softwareentwicklung ist komfortabel. Datenbank, Hostkopplung und Visualisierung sind durch die Mehrtaskfahigkeit parallel betreibbar. Weiterhin sind neben dem beschriebenen System weitere preiswerte Standard-Softwareprodukte einsetzbar. OS/2 erlaubt gegentiber anderen Mehrtaskbetriebssystemen eine effiziente und schnelle Ansteuerung der graphischen Hardware. In /4/ befindet sich eine ausfuhrliche Darstellung tiber die technischen Datails der Implementierung.
84
8. AIternativen
Neben dem Einsatz eines modernen PCs unter MS-OS/2 wurde auch die Verwendung von sogenannten Grafikworkstations unter UNIX (Warenz. v. AT&T) erwogen. Jedoch lagen bier die durchschnittlichen, vergleichbaren Hardwarekosten ungefahr doppelt so hoch, wobei oft im Grafikbereich nicht die gleiche Leistung (Vektoren/sec, Auf16sung) erreicht werden konnnte. AuBerdem ware ein Unix-Rechner am Einsatzort in einer MS-DOS-Umgebung ein Fremdk6rper gewesen.
9. Erfahrungen
Gegeniiber anderer Software, welche den Entwurf von Karten in den Vordergrund stellt, bietet das beschriebene System die M6glichkeit zum interaktiven Arbeiten, d.h. zum schnell en Sichten und Arbeiten mit MeBwerten und Regionen. Erm6glicht wurden die sehr kurzen Bildaufbauzeiten durch Mehrebenentechnik sowie die intelligente Ausnutzung der Fabigkeiten moderner Graphikhardware.
Der Einsatz von Standardkomponenten halt die Beschaffungskosten niedrig, und macht den Arbeitsplatz auch iiber den vorgesehenen Zweck binaus nutzbar.
Die Erfahrungen mit Softwareentwicklung unter MS-OS/2 k6nnen als durchweg positiv bezeichnet werden.
10. Literatur/lnfonnation
Ein Oberblick iiber das System befindet sich in /5/, /6/ ist eine umfangreiche Produktbeschreibung mit vielen Einsatzbeispielen.
/1/ BGA Tatigkeitsbericht 1988, Berlin 1988.
/2/ Weigert, P., Miiller, J.: "Analytical Quality Assurance in the German Food Contamination Monitoring Programme", Fresenius Zeitschrift fur analytische Chemie, 6/1988, pp.736-737.
/3/ GH Kassel, FB Wirtschaftwissenschaften, Dr. Hans-Friedrich Eckey, Nora-Platiel-Str. 4, 35 Kassel.
/4/ Woznik, W: "Kundenspezifische Entwicklung unter OS/2", CAD-CAM Report, Juni 1989, pp.70-73.
/5/ Miiller-Meernach, R.: "Interaktives MeBwertdarstellungssystem", Infografik, Feb. 1989, pp.26-28.
/6/ KoDaM Produktbeschreibung, 22 Seiten, S.E.P.P. GmbH, Lohmiihlweg 4, 8551 R6ttenbach.
85
Bundesland Bayern
Urtterfronken
Rottenboch
Landk ...... 1n Boy.m Re<jie<unqstJ..zi<1<e ""p!odiert dorgeste'lI.
Schwo~
Oberfmnken
Londtogswohlen 1982
BAYERN
• CSU .SPO • FOP
Sonstige
The mat i s c h e Kart e n mit KoDaM von S.E.P.P.
KoDoM von S.E.P.P.
86
Th em ati sc h e Ka rte n KoDaM von S.E.P.P.
~ )c~le~wi~-~ol~tein [J ~om~~I~ ~ ~ie~elsoc~sen
, ~Iemen
Anfoll und 8eseitigung von Klarschlomm in 6ffentlichen Klaronlogen 1979/1983 Quell • • Ooten LU,Umw.11 1986/87
KoDoM von S.E.P.P.
87
Th em oti sc h e Ka rte n KoDaM
280.0
• A~ID~m • ~OOI~OSt ieren III Yer~rmen ,I lDn~,i rtsc~Dltl. Yerrertun~ o )onsti~es
von
Anfall und Beseitigung
S.E.P.P.
von Klarsch lamm in ()ffentlichen K lora n I age n 1 979 Quell • • Dolen '". Umwell 1986/87
KoDoM von S.E.P.P
88
The mat i s c h e Kart e n
KoDaM von S.E.P.P.
lx.P (c) Copyright. 1988: S.E.P.P. GmbH 8551 ROttenboch. LohmOhlweg 4 Tel. : 09195/3025 Fox.: 09195/1739 -
• - 1500.950 • 1500 .950 - 2000 .900 • 2000.900 - 2500.m • mO.8SO - lOOO.800 • 5000 .800 - mO.150
5500.150 - 4000 .100 4000.100 - 4500 .650
D 4500.650 - 5000 .600 D 500UOO - 5500 .550
~ 0 5500.m - 6000.500 o 6000 .500 - 6500. 451 • 6500.m - 1000 .401 • 1000AOI - 1500.351 • 1500.351 - 8000.301 • 8000 .301 - 8500 .251 • 8500 .251 - 9000.201
9000.201 - 9500.151 9500.151 - 10000 .1 02
o 10000.102 - 10500.052 D 10500.052 -
KoDoM von S.E.P.P.
IKoD
aM:
them
ati
sch
e K
art
en
I O
logr
omm
1
f! .6
.. • 2
-.2 -.'
-.8
Ja
n Fe
b M
it Ol
Hri«:
m\;c
"
Olo
grom
m 2
:rIJ
ij O
J <0
-.2
-.f
II ~~
erf ra
n~en
-.. D
~~rto
ll~l~
Jc.1'I
f"e
b va
, lk
It"r
lr.all
,lcu"
O
logr
omm
J
II ~i
ltell
ran~e
n D
~~lr~
i~t
I .. . f
II ~n
terlra
n~en
.2
t==]
~~i~
I -.
2
-.f
-.. Jo
n f~
6 "i
r Itl
Ueltr
onke
n
Ob
er-
, U
nte
r-un
d M
ittel
fran
ken
KoO
ot.l
von
S.E.
P.P
.
Bu
nd
esla
nd
N o
rdrh
ein
-W
estf
alen
•
-51
5435
2.00
0 .
5154
352.
000
-51
9770
4.00
0 .5
19
77
04.0
00
-52
4105
6.00
0 .5
24
1056
.000
-
5284
408.
000
.52
84
408
.000
-
5327
760.
000
5327
760.
000
-53
7111
2.00
0 ~ 5
3711
12.0
00
-54
1446
4.00
0 05
4144
64.0
00
-54
5781
6.00
0 05
4578
16.0
00
-55
0116
8.00
0 05
5011
68.0
00
-55
4452
0.00
0 05
5445
20.0
00
-55
8787
2.00
0 .
5587
872.
000
-56
3122
4.00
0 .
5631
224.
000
-56
7457
6.00
0 .5
6745
76.0
00 -
5717
928.
000
.571
7928
.000
-57
6128
0.00
0 .5
76
1280
.000
-
5804
632.
000
5804
632.
000
-58
4798
4.0
00
~ 5
8479
84.0
00
-58
9133
6.00
0 0
58
913
36.0
00
-59
3468
8.00
0 05
9346
88.0
00
-
Oo,
.slell
unq
bi's
ouf
Cem
einC
leeb
enc
The
ma
tisc
he K
art
en m
it K
oDaM
vo
n S
.E.P
.P.
KoDo
M v
on S
.E;P
.P.
<0 a
Die verschiedenen Ebenen bei der Bearbeitung von Umweltdaten mit Geoinformationssystemen
Zusammenfassung
Peter Riegger EDV Studio Ploenzke
MUnchen
Projekterfahrungen bei UmweltbehOrden haben eine in den Anforderungen ganzlich unterschiedliche Behandlung von Umweltdaten in Geoinformationssystemen aufgezeigt. Diese Erkenntnis hat zu einer strategischen Abgrenzung geftihrt, die sich auf drei Ebenen beschreiben laBt:
1. UnterstUtzung von Umweltkatastem auf der Verwaltungsebene, 2. der Umweltarbeitsplatz zur analytischen Bearbeitung von Umweltdaten; 3. die hybride Grafik flir die Integration von Raster- und Vektordaten.
Vom Umweltkataster wird eine exakte Nachbildung von Grenzverlaufen erwartet, Systemvoraussetzung ist ein Vektorgrafiksystem mit Sachdatenhaltung. Der Umweltarbeitsplatz stellt geringere Anforderungen an die Kartographie, verlangt dagegen viele und sehr leistungsfahige methodische Komponenten. FUr die Sachdatenhaltung sollte eine Datenbank-Schnittstelle existieren. Bei den hybriden Grafiksystemen wird die Leistungsfahigkeit von vektorbasierten und rasterbasierten Systemen ausgenutzt. Ihr wichtigstes Einsatzgebiet ist derzeit die Integration von Femerkundungsdaten.
1. Einleitung
Die Forderungen nach wirksamen Instrumentarien zur Bewaltigung der Umweltkrise hat sich mittlerweile auf allen Verwaltungsebenen ausgebreitet. Eine konkrete Vorstellung dessen, was wirklich gewUnscht wird, ist allerdings in den wenigsten Fallen gegeben. Am liebsten ware dem Kunden ein Instrumentarium, das all seine Probleme bewaltigen konnte - jedoch nach seinen Vorstellungen. Diese konnen - soweit keine Rechtsnormen und Verwaltungsvorschriften den Rahmen bestimmen - recht subjektiv sein und werden tei!weise durch regionale Besonderheiten und politische Zielsetzungen gepragt.
Dies trifft besonders auf die Integration von Geo-Informationssystemen (GIS) in die Umweltarbeit - vorzUglich in Umweltinfommtionssysteme (UIS) - zu. Hier konnen sogar schon Unterschiede in der Verwaltungsstruktur zu verschiedenen Losungswegen fUhren. So ist es bei Kommunen, die Uber eine l?~gene Vermessungsverwaltung verfligen, selbstverstandlich, die dort aufgebauten Systeme mit in die Uberlegung einzubeziehen. Dies geht teilweise soweit, daB die Vermessungsverwaltung - wei! in ihr die Know-how-Trager sitzen - mafigeblich fUr die Integration der Umweltdaten in ein GIS verantwortlich ist.
Vielerorts hat diese Vorgehensweise zu Spannungen gefUhn, denn die Umweltamter sind oft weniger an groBen Verwaltungssystemen interessiert, sondem wUnschen sich vielmehr fUr sie leichter verstandliche - und das heiBt ihren Problemen (Analyse der geographischen Information) angepaBte -GIS. Dieser Ansatz verlangt vor allem methodische Bausteine flir die Auswertung der alphanumerischen (Statistik) und graphischen (Verschneidung) Daten. Der Anspruch auf die geometrische Genauigkeit der I?aten ist weniger hoch, denn auch in der Natur gibt es kaum exakte Abgrenzungen, sondem meist nur Ubergange von einen Zustand in den anderen.
Auch fUr diese Anforderungen bietet der Markt Produkte an. Oft passen sie sich dem allgemeinen Zustand der Umweltamter an: Auf dem eventuell sogar von zu Hause mitgebrachten PC werden die ersten Programme geschrieben, weitere Ausbaustufen sollen nun auf dem angeeigneten Computerver-
92
standnis aufbauen. Kommen die Umweltexperten gerade von der Universitat, dann werden vielleicht noch die dort eingesetzten Systeme adaptiert - wobei man schon von Gliick sprechen kann, wenn jemand in der landschaftsplanerischen, geographischen oder okologischen Ausbildung iiberhaupt an einer gescheiten DV-Ausbildung teilhaben konnte.
Festzuhalten ist die Tatsache, daB diese kleinen Systeme bei den groBen Datenmengen und komplexen Aufgabenstellungen, die im Umweltschutz anzutreffen sind, schnell an die Kapazitatsgrenze stoBen. Selbst wenn die Entwicklung von PC's in dem Bereich der zu verarbeitenden Datenmenge rasch voranschreitet, so kann sie dennoch nicht mit den jetzt schon sehr dringend anstehenden Notwendigkeiten an LOsungsansatzen im Umweltschutz Schritt halten. GIS miissen integraler Bestandteil eines vernetzten Systems sein, und dementsprechend sind auch die Standards zu definieren.
Es sollte als eine selbstverstandliche Tatsache angesehen werden, daB die Belange der Vermessung nicht unbedingt denen des Umweltschufzes entsprechen. Dies muB auch in der Systembetrachtung beriicksichtigt werden. Das heiBt nicht, daB ein GIS, das fUr die Vermes sung gut ist, deshalb fUr den Umweltschutz nicht taugt. Bei ein und demselben System konnen ja durch das modulare Angebot die verschiedenen Aufgabenbereiche beriicksichtigt werden. Kernaussage solI sein, daB die Systeme nicht unbedingt identisch sein miissen. Wichtig ist vielmehr, daB der Transfer von Daten gewahrleistet ist, urn den Aufbau von redundanten Datensatzen zu vermeiden.
Damit waren zwei Systemansatze beschrieben, wobei der erste, der sich mit der Datenverwaltung beschaftigt, als Umweltkataster und der zweite, der sich auf die geographische Analyse bezieht, als Umweltarbeitsplatz bezeichnet werden solI. Bliebe noch das hybride graphische Informationssystem, das als reale Zukunftsperspektive mit in die Betrachtung einbezogen werden solI. Bisher widerstrebt die preisliche Kalkulation, als auch die Frage nach der Niitzlichkeit von Bildverarbeitungsfunktionen, dem Einsatz solcher Systeme. In einem integrierten Konzept werden sie aber dennoch bald Beriicksichtigung finden - wenn auch vorerst als Einzelerscheinung fUr besondere Aufgabenstellung (Bild 1).
Umweltkataster (MERKIS)
• Raumbezugsbasis • Fachdatenbasis • Geographische
Datenbasis
\
Umweltarbeitsplatz
• Datenbankkonzept • Methodeneinbindung • kartographische
.. Grund- und Dar-steliungsfunktion
/ hybride Graphik
• Einbeziehung von Fernerkundungsdaten
• Konvertierung von gescannter Information
Bild 1. Das Drei-Ebenen-Konzept fiir den Einsatz von GIS im Umweltschutz
2. Datenstruktur
Stiefel /STIE87/ unterscheidet GIS nach den zwei verschiedenen Ansatzen in vektorbasiene und rasterbasierte Systeme. Damit waren durchaus die beiden wichtigen Graphikkomponenten angesprochen, jedoch sind GIS noch durch spezifische Methodenansatze sowie die Sachdatenhaltung charakterisiert. Geospezifische Methodenansatze sind z.B. das digitale Gelandemodell (DGM) und die Rasterzellenverarbeitung (Griding). Nicht alle verwendete Methoden mlissen integraler Bestandteil eines GIS sein (Statistikverfahren, Zeitreihenanalyse u.a.). Ebenso ist die Sachdatenhaltung nicht auf das GIS beschrankt. Deshalb mlissen von einem GIS Schnittstellen zu extemen Methoden- und Datenbanken bestehen (Bild 2).
2.1 Vektorgraphik
Der GroBteil der heutigen GIS hat sich aus den vektorbasierten Systemen entwickelt. Vektordaten werden aus vermessungstechnischen Berechnungsprogrammen, der photogrammetrischen Auswertung von Luftbildem oder durch Digitalisieren analoger Vorlagen gewonnen. Die einzelnen Objekte (Punkt, Linie, Flache) werden dabei durch die Speicherung von Deskriptoren (geometrische, qualitative, quantitative und evt. Namensinformationen) eindeutig beschrieben.
Flir die Speicherungskonzepte nimmt der Punkt eine Sonderstellung ein, da seine geometrische Information schon durch ein einziges Koordinatenpaar eindeutig festgelegt ist. Linien entstehen durch die Verkettung von Punkten mitte1s Zeichenvorschrift. Flachen konnen durch Verweis von Flachen auf Linien oder relational durch den Verweis von Linien auf Flachen abge1egt werden. Modeme kartographische Datenhaltungskonzepte zeichnen sich durch eine horizon tale (blattschnittlose) und venikale (Zeichenebenen) Speicherung der geometrischen Information aus /SCHI85/.
2.2 Rastergrafik
Unter dem Begriff Rastergraphik wird die Ubemahme von Bildem in digitaler Form verstanden. Rasterdaten konnen zum einen durch Femerkundungssensoren, zum anderen durch Scannen analoger Vorlagen erfaBt werden. Die Rasterbilder setzen sich aus Pixeln zusammen, deren qualitative Information sich aus den Farbwerten ergibt. Flir die Manipulation der Rasterdatem stehen die Methoden der Bildverarbeitung (Grauwertoperationen, Filterverfahren u.a.) zur Verfligung /GOEP87/.
Je nach Informationstiefe (Frequenz) und Bildintensitat (Kontrast) benotigen Rasterdaten sehr viel Speicherplatz. Kompressionsverfahren, wie runlength-Codierung oder Quad-Tree-Strukturen, konnen bei binaren und homogenen Bildem zu einer erheblichen Einsparung von Speicherplatz flihren - flir Satellitenbilder mit einer 8 Bit Tiefe benotigen sie unter Umstanden mehr Speicherplatz. Deshalb wird die Entwicklung von groBen lokalen Speicherkapazitaten flir die Integration von Rasterbildem in komplexe GIS von immenser Wichtigkeit sein.
2.3 DigitaJes GeJiindemodel\
Mark /MARK79/ zeigt die Vielfalt der Ansatze zur digitalen Abbildung der Gelandeoberflache. Die beiden wichtigsten Typen sind die Raster-Gelandemodelle und die triangulienen irregularen Netzwerke (TIN). Raster-GeIandemodell werden aus beliebig im Raum verteilten Hohenangaben interpoliert und ergeben ein dichtes liber das Relief gelegtes Gitter mit den Hohen an den Schnittpunkten. TIN's werden aus den meist eigens aufgenommenen Gelandepunkten mit dem hochsten Informationsgehalt durch Dreiecksvermaschung gebildet.
Die Beziehung von TIN's zur Vektorgraphik ist zwar naherliegend und damit ware die Ubemahme ,:,.on Flacheninformationen 1eichter moglich, andererseits beglinstigt die Analogie der Datenstruktur die Ubemahme der Rastermodelle in Bildverarbeitungssysteme. Ausserdem kann bei Rastermodellen die
94
Geogr aphisches Informationssys tern Dig itales Gelandemode l/
Sachdaten
Tools Prozessebene DB/DC
Bild 2. Datenstruktur in einem Geoinformationssystem
Interpolation von verstreut vorliegenden Hohenwerten auch auf andere Z-Werte bezogen werden. So konnen z.B. die ImmissionsmeBwerte verschiedener Stationen flir eine Isoliniendarstellung ausgewertet werden. Allerdings zeigt die Praxis, daB hierbei noch spezifische Parameter zu beriicksichtigen sind (im Beispiel die Windrichtung und -geschwindigkeit), die in einer vorgelagerten zusatzlichen Modellkomponente zu berechnen sind.
2.4 Rasterzellenverarbeitung
Ahnlich wie bei der Rasterverarbeitung wird auch bei der Rasterzellenverarbeitung eine Matrix angeiegt, in der Sachinfortnationen verschllisselt umgesetzt werden. Hierbei handelt es sich urn ein alphanumerisches System, mit dem unter Verzicht auf graphische Komponenten verschiedene Rauminfortnationen
95
miteinander aggregiert werden k6nnen. Der wesentliche Tell eines Gridpaketes besteht aus Verfahrenskomponenten zur Verkniipfung der verschiedenen Matrizen.
Gridingmethoden sind eigentlich die ersten Ansatze zur DV-unterstiitzten Verarbeitung raumbezogener Daten und vor allem in den geowissenschaftlichen Berekhen angesiedelt. Die Verfahren stammen meist aus Amerika, haben sich aber in Deutschland bisher kaum durchgesetzt /siehe u.a. BF AL82/. Eine Folge mogen die schlecht lesbaren Druckausgaben gewesen sein, die oft nur mit dem alphanumerischen Zeichensatz erstellt wurden. Mittlerweile werden Datenaufnahme und -ausgabe durch Schnittstetlen zu Vektorsystemen untersttitzt. Da auch Grids eine Analogie zu den anderen Rasterstrukturen aufweisen, konnten sich die Verfahren in der naheren Zukunft vielleicht besser durchsetzen.
2.5 Sachdaten
Graphischen Objekten konnen Sachdaten zugeordnet werden, die qualitative und quantitative Merkmale beinhalten. Die Verbindung wird entweder in Form von Pointern oder tiber Objektnamen herge- steHt. Diese Zuordnung wurde bisher nur flir die Vektorgraphik im kommerziellen Rahmen realisiert. Die Moglichkeiten der Objektbildung in der Rastergraphik werden derzeit noch untersucht bzw. sollen in hybriden graphischen Systemen umgesetzt werden /FRIT88/.
Die Sachdatenhaltung wird in den meisten Fallen auf einer relationalen Datenbank verwirklicht. Diese kann integraler Bestandteil des GIS sein, sie kann aber ebenso extern tiber Schnittstellen angeschlossen werden. Der Vorteil der internen Datenbank liegt in den schnelleren Zugriffen in beiden Richtungen, wohingegen eine externe Datenbank bessere AnschluBmoglichkeiten fUr andere Anwendungen bietet. Nouhuys /NOUH86/ empfiehlt aus Grunden der Praktikabilitiit eine Trennung der Speicherarten von geometrischen Sachdaten (gemeinsame Verwaltung mit Graphikdaten) und disaggregier-ten Daten (extern in Dateien oder Datenbanken). .
3. Umweltkataster
Verwaltungsbereiche, die sich in Umweltplanung und Umweltschutz auf genaue Abgrenzungen sttitzen, sind auf einen exakten Raumbezug angewiesen. Dieser Raumbezug kann nur durch ein vermessungstechnisch ausgereiftes GIS realisiert werden. Basis flir ein Umweltkataster muB deshalb ein Vektorsystem sein. Die Struktur fi.ir solch ein System wird durch die MaBstaborientierte Einheitliche Raumbezugsbasis fi.ir Kommunale Informationssysteme (MERKIS) beschrieben IDST88/.
Der einheitliche Raumbezug soIl schrittweise in verschiedenen MaBstabsebenen, die auf dem GauB-Kriiger-Koordinatensystem basieren, realisiert werden. Die MaBstabsebenen entsprechen den bisherigen analogen Kartenwerken. Auch die Uberlagerung mit den fachbezogenen Geometrien folgt den bisherigen Planungs- und Verwaltungsstandards. So setzen Leitungskataster und verbindliche Bauleitplanung auf die Stadtgrundkarte (l :500), das Altlastenkataster und die Flachennutzungsplanung auf die Deutsche Grundkarte (1 :5000) auf.
In der Fachdatenbasis werden die Sachdaten gehalten. Wahrend die Raumbezugsbasis zentral in einem Fachbereich vorliegt, werden die Fachdaten in den verantwortlichen Fachbereichen fortgeftihrt. Die Abfrage der Daten kann durch geometrische oder alphanumerische Sachfragen erfolgen (Bild 3). Obwohl Geometrie- und Sachdaten jeweils in eigenen Datenbanken verwaltet werden, sollte deshalb aus der Sicht des Anwenders nur ein logischer Zugriff auf die verschiedenen Datenbestande notwendig sein. Eine Umsetzung der in MERKIS formulierten Zielsetzung zeigt der Umweltschutzbericht der Stadt Wuppertal IWUPP87/. Digitale fachbezogene Geometriedaten werden flir die Umweltbereiche Altablagerung, Bodenproben, Landschafts-, Natur- und Wasserschutz erhoben. Andere Bereiche, wie Luftbelastung und Gewassergtite werden tiber thematische Bausteine ausgewertet. Durch die Novellierung des Bundesbaugesetzes hat die Aufnahme der Altlasten besonders an Bedeutung gewonnen (Bild 8, Farbtafeln). Die genaue Abgrenzung der Altlastengebiete ist nun flir die Aufstellung des BebauungspJa,ns unerliiBlich. ... -
96
Raumbezogenes Informal ionssyslem
FachdatenbaSls
Bild 3. Zugriffsmoglichkeiten in MERKIS
Einen Schritt weiter geht der okologische Planungsansatz von Berlin /SENA87/. Uber Auswerteverfahren wird hier eine Bewertung der Belastungspotentiale von Altlasten erstellt (Bild 9, Farbtafeln). Dieser Ansatz sieht den Weg zu einem raumlichen Umweltinformationssystem erst durch die Datenanalyse, -bewertung und Wertaggregation gegeben. Die Verfahren der Bewertung und Wertesynthese sind Bestandteil einer Methodenbank, die von der Graphik aus interaktiv gesteuert werden kann. Dies erzeugt Transparenz, Nachvollziehbarkeit und Flexibilitat.
Die yom Anwender erzeugten Nutzenfunktionen und Synthesealgorithmen werden in Modulbibliotheken verwaltet und konnen flir weitere Anwendung unter veranderten Bedingungen als Simulationsmodelle eingesetzt werden. Gesteuert werden sollen die Verfahren durch eine sich standig erweitemde Benutzeroberflache, die in Form von Masken, Dialogprozeduren und Mentis aufgebaut wird. Ziel dieses relativ aufwendigen Ansatzes ist die moglichst einfache Zuganglichkeit und Anwendbarkeit des Instrumentariums /BOCK89!.
4. Umweltarbeitsplatz
Durch die Dringlichkeit von MaBnahmen im Umweltschutz, einem wachsenden UmweltbewuBtsein sowie der raschen Entwicklung neuer Umwelttechnologien fallen immer groBere Mengen von umweltrelevanten Daten an, die vorerst volIig unstrukturiert sind. Diese Daten bilden Planungs- und Arbeitsgrundlage flir eine Vielzahl von Bereichen (Politik, Verwaltung, Wissenschaft, OffentIichkeitsarbeit
97
usw.). Um die Daten auf versHindliche Weise zuganglich zu machen, sind Hilfsmittel, wie graphische Datenverarbeitung, Textverarbeitung oder Methodenbausteine unerHiBlich.
Die'se Anforderungen verlangen ein spezifisches Instrumentarium, das hier durch das Konzept "Umweltarbeitsplatz" beschrieben wird (Bild 4). Der Umweltarbeitsplatz solI folgende Aufgaben iibernehmen:
>- externe Kommunikation >- Sachdatenverwaltung >- Geometriedatenverwaltung und -darstellung >- Datenanalyse >- graphische Darstellung von Auswertungen >- Dokumentenerstellung und -verwaltung
DATEN -
BANK
Bild 4. Konzept eines Umweltarbeitsplatzes
BORD -
ANWEND.
Damit ist die Funktionalitat eines pes schnell liberlastet. Zunehmend zie1en Ausschreibungen auf UNIX-LOsungen abo UNIX-Systeme bieten mitlerweile eine umfangreiche Software flir Blirofunktionen, Datenverwaltung, Kommunikationsnetze und wissenschaftliche Analyse. Damit der Anwender sich nicht intensiv einarbeiten muG, ist eine leicht zugangliche Benutzeroberflache (Fenstertechnik) mit Moglichkeiten zu Interprozesskommunikation und Multitasking gefordert.
Die kartographischen Anforderungen werden auf ein weniger komplexes Vektorsystem beschrankt. Flir die projektbezogene Analyse sowie Ausgabe von Ubersichtskarten reicht die blattschnittbezogene Speicherung der Geometriedaten. Die Sachdaten werden in einer relationa1en Datenbank gehalten~ Die
98
Verkniipfung der Geometriedaten mit den Sachsatzen erfolgt liber Graphik- und Sachsatzschllissel. AuBerdem ist eine methodenbezogene Datenbankkopplung vorgesehen.
Zlir Diskussion steht die Anbmdung von Rastergraphik und Griding. Durch die Rastergraphik soli zum einen die Unterlegung mit gescannten Kartenvorlagen realisiert, zum anderen die Ubernahme von interpretierten Fernerkundungsdaten ermoglicht werden. Eine Rasterdatenbasis kann auch Grundlage flir den Einsatz von Gridverfahren sein. Ein sinnvoller Ansatz ware, die Verfahren auf einer gemeinsamen Rasterdatenbasis aufzusetzen. 1m einzelnen ist diese Entscheidung jedoch von dem eingesetzten Gridpaket abhangig.
is 18 17 16 15 14 13 12 11 10 ~ 8 7 6 S 4 3 2 1 o
Lu.fte.per-otur Ucsse,.tc.pcro'tur
VERDEN VERDEN
1~87
1S87
_1 ~~dr--~ __ ~ __ ~ __ L-__ L-__ ~~ __ ~ __ ~ __ ~ __ ~ __ J
-2 -3 -4 -5 - 6
Bild 5. Jahresgang von Luft- und Wassertemperatur an der Me6station Verden
In einer ersten Erprobung des Umweltarbeitsplatzes wurden in die Datenbank flir 1987 die 1/2 Stunden Mittelwerte von Oz-Gehalt, ph-Wert, Temperatur Luft/Wasser und elektrische Leitfahigkeit aus 34 MeBstationen des automatisierten Gewasserliberwachungsnetz Niedersachsen (GUN) einge1esen. Damit ist ein Datensatz von 40 MB zu verwalten. Die ausgewerteten Daten konnen wahl weise in einer Ubersichtskarte oder einer Geschaftsgraphik (Bild 5) betrachtet werden.
5. Hybride grafische Informationssysteme
Kammerer u.a. /KAMM881 beschreiben die Entwicklung von hybriden graphischen Informationssystemen in drei Stufen:
99
Erste Stufe: Auskunftssystem
Uberlagerung von Vektor- und Rasterdaten, wobei der Benutzer dafUr zu sorgen hat, daB beide Datentypen in demselben Koordinatensystem vorliegen und der identische Ausschnitt zur Darstellung kommt.
Zweite Stufe: Integrierte Datenverarbeitung
Gemeinsame Verarbeitung auf der Basis von Schnittstellen sowie Konvertierungsm6glichkeiten. Koordinaten k6nnen fUr den gegenseitigen Zugriff ausgetauscht werden, Verarbeitungsfunktionen werden unter einer einheitlichen Benutzeroberflache zusammengefaBt.
Dritte Stufe: Integrierte Datenverwaltung
Gemeinsame Verwaltung von Vektor-, Sach- und Rasterdaten. Unterschiedliche Datenhaltung, aber gemeinsame Verwaltungsoberflache, so daB sich fUr den Anwender die drei Welten als logische Informationsebene eines Beschreibungsmodells darstellen.
Dies fUhrt zu einem Instrumentarium, mit dem erstmals aile Vorteile der verschiedenen Ansatze der Datenaufnahme, -haJtung, -verarbeitung und -ausgabe in einem System vereint werden k6nnen (Bild 6).
[IJ~~ Ql 19J 0~§ Gel.,jndeo· lullbllCl · Olgll,ll ' Alph;)num Vo,h .. nd~ne ceo· 'ilcannt'f \alell.l ,)ufn~"'me me,\u('Il~ \'t" un9 TlI:fmln.ll ~dcnd,uen Kamer"
l- t I-Ve'l.CIfd~IM I Sathd~ICiI I I R"l'tf~LtII
I I I ~ l- I-
Vt 101°. Sach- und R~'ltr~ttfWttw"~
HrrYd.Vetll~
K=) o Ubo''''J'''"'I v·R W Vcll.Cf ... e-riol~1VII9 ov R RI!o[fI"Wtf"'Dl:l~
O.-v DAII\'Wtf"""9
Hybrid •• G,.fik.y.tem
Sllll. I M<10- I P.lllSolOMo I I Il,ulitf " I I m·o- I oIotlet , .... lc<Attlllsdt """I<' IIttcOl(ft co;.,
Bild 6. Aufbau eines hybriden graphischen Informationssystems
Hybride Grafik bietet sich im Umweltbereich vor aHem bei der Integration von Fernerkundungsdaten zur Umweltiiberwachung und zum Scannen von topographischen Karten als Hintergrundinformation an. Auf3erdem ermoglicht die Konvertierung von Vektor- in Rasterdaten vereinfachte und Schnellere Ver-
100
schneidungsroutinen, die durch entsprechende Datenkomprimierung (Quad-Trees) noch komfortabler gestalltet werden konnen. Dber die Organisation von Vektor- und Rasterdaten auf Objektebene wird vielleicht schon in naher Zukunft der gemeinsame Zugriff auf Attributdaten verwirklicht. Damit erlauben hybride grafische Systeme den umfangreichsten Losungsansatz zur Bearbeitung und Visualisierung umweltrelevanter Daten.
Dber den Einsatz von Femerkundungsdaten zur Umwelttiberwachung gibt es mittlerweile eine Menge an Literatur /HA YD87 /. Neue Arbeiten beschaftigen sich mit dem Einsatz von Femerkundungsdaten in hybriden graphischen Informationssystemen. Gegg /GEGG89/ untersucht die Moglichkeit zum Einsatz von Klassifizierungsverfahren in einem hybriden graphischen Informationssystem. Das Beispiel bezieht sich auf versiegelte Flachen (Siedlungsflachen und Verkehrsnetz). Dabei konnte die automatische Klassifizierung von "Siedlungen" und "Autobahnen" erfolgreich durchgeflihrt werden, wohingegen sich die Erfassung von StraBen wegen der geringen Breite als problematisch erwies. Optimiert wurde die visuelle Erkennbarkeit von Siedlungen und Verkehrswegen durch die Verwendung eines Winterszenarios (Bild 10, Farbtafeln), das sich allerdings flir die Klassifizierung nicht einsetzen lieS. Selbst wenn das Ziel dieser Arbeit die automatische Klassifizierung ist, so kommt die Informationsgewinnung bisher nicht ohne das Expertenwissen des Bearbeiters aus.
Den Schritt zur Integration von Satellitenbildauswertungen in ein GIS geht Kerl /KERL88/. Dber die Kombination von spektral hochauflosenden Landsat -TM -Daten mit den geometrisch hochauflosenden Sojuzbildetn erstellt sie eine Karte der Vegetationsverteilung in der Stadt Milnchen. Aus einer geographischen Datenbasis ilbemimmt sie die als Vektordaten hinterlegten Stadt- und Bezirksgrenzen. Aus den geschlossenen Polygonzilgen konnen Masken erstellt und mit der Vegetationsverteilung verschnitten werden. So lassen sich zu den einzelnen Bezirken Flachenbilanzen ausgeben (Bild 11, Farbtafeln). Fi.ir Breadt /BRAE89/ ist die Einbeziehung von Femerkundungsdaten in ein GIS nicht nur aufgrund der wachsenden Informationsgrundlage notwendig. Auch werden zur Interpretation von Satellitenbildem in digitaler Form vorliegende Bodenkarten benotigt. Nach seiner Ansicht mils sen hybride grafische Informationssysteme integraler Bestandteil in ein Umweltinformationssystem sein.
6. Das GIS als Bestandteil eines Umweltinformationssystems
Dberlegungen zum Aufbau von UIS werden zur Zeit in verschiedenen Kommunen und Landesbehorden angestrengt. Der konkreteste Ansatz wird derzeit in Baden-Wilrttemberg umgesetzt. Hier solI ein, in seiner Komplexitat einmaliges System realisiert werden, an dem 8 Ressorts der Landesverwaltung beteiligt sind, das aile Umweltbereiche integriert, Aufgaben von der Frilherkennung bis hin zur Fachplanung bewaItigt, samtliche Verwaltungsebenen unterstiitzt und dabei noch eine flachendeckende Aufnahme gewahrleistet. Der Umfang eines solchen Systems ist schwer zu modellieren. Bild 7 soll zumindest einen Eindruck von den wesentlichen Komponenten geben. Dabei muS sich das UIS mit einer groBen Menge heterogener Daten auseinander setzen. Die Daten zu den Umweltmedien werden in Listen, Katastem und Karten geflihrt, aus Luft- und Satellitenbildem gewonnen oder ilber manuelle bzw. automatisierte MeBnetze erfaBt. Ob die Daten zentral oder dezentral gehalten werden - die Organisation wird immer von einer Zentrale aus geregelt. Dort werden die Kataloge zu den Umweltdaten, den Methoden (Tools) und den Wissensbereichen (Pools) geflihrt sowie die Kommunikation im Netz und an den Bildschirmen bestimmt.
Die Tools lei ten sich aus den zu bewaItigenden Aufgaben abo So sind neben den normalen Schreibtischtatigkeiten und der Dokumentation besonders Prognosemodelle und Planungsinstrumentarien gefragt. Optimale Unterstiltzung im Verwaltungsablauf konnen Expertensysteme beisteuem. Gerade bei den Expertensystemen ist der Zugriff auf Pools von Bedeutung. Damit werden zusatzliche interne und exteme Datenhaltungen bezeichnet, in denen fachiibergreifende und fachspezifische Informationen abgelegt sind.
Die kartographischen Anforderungen werden besonders aus der Dateneingabe und -ausgabe ersichtlich. So lassen sich zu der Beschreibung der Umweltmedien folgende Beispiele herausnehmen:
101
Umweltmedien Natur und I W Id I Landschaftj a Luft I L!lrm I Wasser 1 Abfall I Boden I Nahrung 1 =~~~?t~t Tools Textverarbeitung
Kommunikationssystem D
S a t
y e
Pools Dokumentendatenbank
Desktop Publishing s n t b Statistiken
Geschattsgrafik
Vektorgrafik
e a m Kern n v k e System r
Gefahrstoffdatei
Geoinformationssystem
Griding w Bibliographie a
Bildverarbeitung I Umweltforschungskatalog t
Methodenbank u Umweltvorschriften n
Exportsystem 9 Benutzeroberfl!lche Wissensbasis
Verwaltung I Politik I Alarmsystem I Prognosen lDokumentati°1 UVP 1 Forschung
Aufgaben
Bild 7. Komponenten eines Umweltinformationssystems
Natur und Land-schaft Wald Luft Larm Wasser Abfall Boden Nahrung Radioaktivitiit
Biotopkartierung, Artenschutzkartierung, NSG, LSG, Naturraumliche Gliederung Infrarotbefliegung, forstliche Standortkartierung Emissions- und Immissionskataster, FJechtenkartierung Larrnschutzzonen, Schallimmissionsplan Gewassergiitekarte, GrundwasserabfluBmodelle, WSG Altlastenkataster, Deponienverzeichnis Bodenversiegelung, Reichsbodenschatzung, Bodenerosion Viehbesatz, Diingemitteleinsatz, landwirtschaftliche Nutzflache Trajektoriendarstellung
Die Aufgabenbewaltigung verlangt besondere DV-Ausstattungen, die sich flir die GIS mit einigen Schlagworten anreiBen lassen:
Verwaltung PoUtik Alarmsystem Prognosen Dokumentation
UVP
Expertensysteme zur Erleichterung von Suchanfragen Auskunftsbildschinne mit aktueller Lageiibersicht Ausbreitungsberechnungen und Katastrophenplan Zeitkarten zum Test von Hypothesen und GegenmaBnahmen Schnittstellenforrnat zur Einbindung in die Umweltberichterstattung (DTP) Aggregation raumlicher Informationen zur Ermittlung von Ressourcenund Konfliktpotential
102
Forschung Weiterentwicklung der Klassifizierungsansatze in der Satellitenbildauswertung
Dieses Aufgabenspektrum kann vorerst nur als Zielvorstellung angesehen werden. Dennoch soli ten diese Uberlegung bei der Planung eines GIS in einem UIS mit einflieBen. In einer Umweltverwaltung mit umfassendem Aufgabenspektrum konnten diese Anforderungen den Ausbau aller drei Ebenen bedeuten. Stark verwaltungsorientierte Behorden sollten sich auf das "Umweltkataster" konzentrieren, fUr die mehr analytisch arbeitenden lnstitutionen geht der Weg zum "Umweltarbeitsplatz". Ais iibergeordneter Ausbau empfiehlt sich fUr beide Bereiche der schrittweise Ubergang zu einem hybriden grafischen Informationssystems.
7. Literatur
/BFAL82/ Bundesforschungsanstalt fur Naturschutz und Landschaftspflege (Hrsg.) Landschafts-Infonnationssysteme. NuL 57 (1982) H.12.
/BOCK89/ Bock, M. Umweltatlas Berlin - Autbau eines okologischcn Planungsinstruments. In: Schilcher, M. und Fritsch, D. (Hrsg.) Geo-Infonnationssysteme. Karlsruhe 1989, S.191-208.
/BRAE89/ Braedt,l. Satellitenbilder als Baustein eines Infonnationssystcms fur Landesplanung und Umweltschutz. In: Schilcher, M. und Fritsch, D.(Hrsg.) Geo-Infonnationssystemc. Karlsruhe 1989, S.247-259.
/DST88/ Deutscher Stiidtetag (Hrsg.) MaBstabsorientierte Einheitlichc Raumbezugsbasis fur Kommunale Infonnationssysteme (MERKIS). KOln 1988.
/FRIT88/ Fritsch, D., Schilcher, M. u. Yang,H. Object Oriented Management of Raster Data in Geographic Infonnation Systems. Int. Arch. Phot. Rem. Sens., 27 (1988), B4, S.538-546.
/GEGG89/ Gegg, G. Untersuchungen zum Einsatz von Klassifizierungsverfahren in einem hybriden grafischen Infonnationssystem. Diplomarbeit, Miinchen 1989.
/GOEP87/ Goepfert, W. Raumbezogene Infonnationssysteme. Karlsruhe 1987.
/HA YD87/ Haydn, R. u. Volk, P. Erkennen von Umweltproblemen im Lufl- und Satellitenbild. GR 39 (1987), H.3, S.316-323.
/KAMM88/ Kammerer, 1. u.a. Hybride Grafik in Geoinfonnationssystemen. Beitrage zum X. Intemationalen Kurs fur Ingenieurvennessung, Sonderdruck, Miinchen 1988.
/KERL88/ Kerl, l. Die Integration von Femerkundungsdaten zum Autbau von Geo-Infonnationssystemen fur stiidtische Gebiete. Diplomarbeit, Trier 1988.
/MARK79/ Mark, D.M. Phenomenon-based data-structuring and digital terrain modelling. Geo-Processing 1 (1979), S.27-36.
/NOUH86/ Nouhuys, J.v. Okologischer Planungsansatz Berlin - DV-Konzept. Forschungsbericht 101 04 048 (1986) UBA-FB.
/SCHI85/ Schilcher, M. (Hrsg.) CAD-Kartographie. Karlsruhe 1985.
/SENA87/ Senator fijr Stadtentwicklung und Umweltschutz (Hrsg.) Umweltatlas - Okologisches Planungsinstrument. Berlin 1987.
/STIE87/ Stiefel, M. Mapping Out the Differences Among Geographic Infonnation Systems. The S. Klein Computer Graphics Review, Fall 1987, S.73-87.
/WUPP87/ Stadt Wuppertal (Hrsg.) Umweltschutzbericht 1987.
Verarbeitung von Umweltdaten unter Real-Time-Bedingungen -
Konzept und prototypische Realisierung
F. Schmidt (IKE), R. Lutkemeyer (BMU), J. Lieser (BGA-ZdB),
A. Bohm, G. Hehn, Th. Muller, W. Scheuermann, V. Sundararaman (aIle IKE)
IKE Universitiit Stuttgart, Pfaffenwaldring 31 7000 Stuttgart-80
Einleitung Am 31. Dezember 1986 trat das 2. Strahlenvorsorgeschutzgesetz in Kraft. Zweck dieses Gesetzes ist es, zum Schutz der Bevolkerung die Radioaktivitat in der Umwelt zu uberwachen und die Strahlenexposition der Menschen im Falle von Ereignissen durch angemessene Mafinahmen so gering wie moglich zu halten. Zur Erreichung dieser Ziele wird ein Integriertes Mefi- und Informations-System (IMIS) [1] eingerichtet und im Auf trag des BMU durch die Zentralestelle des Bundes (ZdB) betrieben. Wir haben ein Konzept fUr solch ein System entwickelt und eine prototypische Realisierung wichtiger Komponenten vorgenommen. Konzept und Prototyp werden anhand der Komponenten zur Behandlung von Daten, die in sogenannten Bundesmefinetzen erhoben werden, vorgestellt. Insbesondere wird auf Datenmodelle und Funktionalitaten eingegangen. Die Struktur von IMIS ist in Abb. 1 zu sehen.
Die Entscheidungszentrale ist beim Bundesminister fur Umwelt (BMU) angesiedelt. Dort wird uber mogliche Mafinahmen und Empfehlungen entschieden. Dies geschieht auf der Grundlage von Vorarbeiten durch die Zentralstelle des Bundes (ZdB), deren Hauptaufgaben das Sammeln, Zusammenfassen, Aufbereiten und Dokumentieren der Daten sind. In diesem Beitrag konzentrieren wir uns auf den Teilaspekt Aufbereiten von Daten und dort wieder auf die Daten, die aus dem Bundesmefinetz kommen. Sie spielen insofern eine zentrale Rolle, da aus ihnen zum einen erste Hinweise auf Storfiille (Ereignisfall) erwartet werden und sie zum anderen zur Validierung der restlichen Daten herangezogen werden [2].
Das Informationssystem unterscheidet im wesentlichell zwei Betriebsarten:
1. Normalbetrieb,
2. Ereignisfall (Intensivbetrieb).
1m Normalbetrieb sind keine Empfehlungen notig. Die anfallenden Daten 1l1iissell dann lediglich bewertet, archiviert und in einem Lagebericht zusa1l11l1engefafit werden.
104
1m Ereignisfall steigt die Menge der MeBdaten gegenuber dem Normalbetrieb stark an. Inhalt und Struktur der zu ubertragenden Daten andern sich aber nicht. Zusatzlich werden zur Bewertung der Daten weitere Informationen notig. 1m Ereignisfall mussen daher eine Vielzahl von Daten zwischen dem BMU und der ZdB ausgetauscht werden. Sie umfassen MeBdaten, prognostizierte Daten und abgeleitete (errechnete) Daten. Diese Daten bilden die Grundlage des Bewertungsprozesses und der daraus abgeleiteten Empfehlungen.
Die Bewertung erfolgt unter zeitlichem Druck und muB so schnell erfolgen, daB die empfohlenen MaBnahmen auch rechtzeitig umgesetzt werden konnen. Dies erfordert die gleichzeitige VerfUgbarkeit einer Vielzahl von Daten, Programmen und Erfahrungen. Sie mussen fUr den Ereignisfall in einer konsistenten Weise integriert sein. U nseren Vorschlag fur solch ein integriertes System zur Analyse von Umweltdaten stellen wir im nachsten Abschnitt zur Diskussion. Wir detaillieren ihn durch Aussagen uber mogliche Datenmodelle und Funktionalitaten und durch die Kurzbeschreibung unserer prototypischen U mweltanalysestation.
Grundlagen des DV -technischen Designs eines integrierten Sy-stems zur Verarbeitung von Umweltdaten
Ziel der Entwicklung integrierter Systeme ist es, konsistente Analysen komplexer Situationen unter ungewohnlichen Randbedingungen durchfUhren zu konnen. Voraussetzung dazu sind die Verfugbarkeit
• konsistenter Daten in wohliiberlegten Strukturierungen und mit problemange-paBten Datenmodellen;
• konsistenter Programme zur Verarbeitung, Diagnose und Prognose der Daten;
• konsistentes Wissen zur DurchfUhrung der Analysen;
• einer einheitlichen Benutzeroberflache zur automatischen Erledigung von Routineaufgaben und zur flexiblen Formulierung von Losungsstrategien von Spezialproblemen;
• flexibler Methoden zur konsistenten Erneuerung von Daten und Wissen und der sie verarbeitenden Methoden.
Eine mogliche Sicht solch eines Systems zeigt Abb. 2. Drei Hauptkomponenten sind erkennbar.
• Daten und Wissensbanken
Sie enthalten zunachst die gemessenen Daten, die sie beschreibenden Informationen und die zugehorigen Fehlerangaben.
Die zweite wichtige Komponente sind allgel11ein vorzuhaltende Infol'l11ationen wie Basiskarten, Stal11l11daten oder Grundinformation zur Interpretation einzelner MeBwerte.
105
Schliei31ich ist eine Bibliothek von Regeln zur Veral'beitung und Verteilung del' Daten aufzubauen. Sie enthalt nicht nur die zu beachtenden Vorschriften, sondern mittelfristig auch Wissen iiber die Auswahl und entscheidungsunterstiitzende Darstellung aktueller Lagekarten .
• Methoden zum Erwerb und Verarbeitung von Wissen
Aufgabe diesel' Systemkomponente sind Pflege und Verarbeitung des in del' ersten Komponente gespeicherten Wissens. Da sich das Wissen aus Daten, Regeln und Beschreibungen zusammensetzt, sind auch in diesel' Komponente entsprechend komplexe Methoden anzustreben. Die Wissensaquisition geschieht zunachst iiber Abfrage, Verteilung und Archivierung del' Daten aus den verschiedenen Netzen. Dies sollte weitmoglichst nach einheitlichen Methoden geschehen.
Die Verarbeitung des Wissens geschieht mit Modellen, wie sie etwa durch die Computerprogramme PARK [3] oder SPEED! [4] bereitgestellt werden. Es ist wichtig, solche Modelle nicht nur an das Gesamtsystem anzuschlieBen, sondern sie in das System zu integrieren, urn dadurch MeB- und Rechendaten in konsistenter Form weiterverarbeiten zu konnen.
Ein zweites Modell, das schon in del' Anfangsphase bereitgestellt werden muE, ist das Modell zur Darstellung von Lageinformationen. Typischerweise erfordert es folgende Schritte:
o Auswahl del' darzustellenden Punkte,
o Interpolation del' MeB- oder Rechenwerte,
o Umsetzung in aussagekraftige Farbskalen,
o Darstellung mit Standard-Graphik-Paket .
• BenutzeroberfHiche
Del' Benutzeroberflache kommt besondere Bedeutung zu. Sie muB die Basisfunktionalitaten so unterstiitzen, daB sie auch von nichtwissenschaftlichem Personal routinemiiBig durchgefiihrt werden konnen. Gleichzeitig muB sie es dem Wissenschaftler erlauben, im Ereignisfall flexible und situationsangepaBte Fragestellungen zu losen. Dazu sind folgende Komponenten notig:
o Ein Optionen-Generator, der zeigt, welche Aktionen im gegenwartigen SystemStatus moglich sind, bei Bedarf erklart, welche Annahmen dazu gemacht wurden und welche Datenanforderungen zu befriedigen sind. Del' OptionenGenerator soUte durch Menii- und Help-Systeme realisiert werden. Das Help-System sollte auch Unterstiitzung bei der Formulierung eigener Problemlosungen bieten. Besondere Aufmerksamkeit verdient dabei die Aktion "GegenmaBnahmen". Hier muB das System in del' Lage sein, die Voraussetzungen zu analysieren, danach mogliche MaBnahmen auszuwahlen und schlieBlich die Abschiitzung del' Folgen zu unterstiitzen. Eine Empfehlung durch das System ist mindestens im gegenwartigen Zeitpunkt nicht vert retbar. tiber erste Ansatze zu solch einer U nterstiitzung haben wir in einem extra Paper [5] berichtet.
106
o Ein Report-Writer, der insbesondere in der Lage ist, statistische Auswertungen fiir Routineberichte in Veri:iffentlichungen und Mitteilungen zu integrieren. Die DarsteUung der statistischen Daten soUte tabeUarisch und iiber Business-Graphik mi:iglich sein.
o Ein Konsistenz-Priifer, dessen erste Aufgabe es ist, ankommende Daten auf VollsHindigkeit zu iiberpriifen, bei unvoUstandigen Daten entsprechende Anforderungen zu erstellen und Statistiken zu fiihren. AuBerdem soll er in der Lage sein, Modifikationen von Daten zu iiberwachen und Vergleiche von Daten zu unterstiitzen .
o Ein Karten-Generator zur Unterstiitzung von Lagedarstellungen.
o Ein Daten-Verdichter, des sen Aufgabe es ist, Daten fiir die einzelnen Aufgaben zu selektieren und in der aufgabenspezifischen Art zusammenzustellen.
o Eine Erklarungskomponente, die in der erst en Ausbaustufe ein von jedem Schritt des Arbeitsprozesses aufrufbares Help-System enthalt, das sowohl aufgetretene Fehler als auch weitere Schritte erHiutert.
Diese Aufgaben sind durch eine Reihe von Basispaketen zu unterstiitzen. Solche Basispakete sind
1. Datenverwaltungssystem,
2. Modellierungssysteme zur Erstellung von Prognosen,
3. Auswertesystem zur statistischen Analyse,
4. Graphikpaket zur Darstellung der Ergebnisse,
5. Report-System zur Dokumentation,
6. Kommunikationssystem zum Austausch von Informationen,
7. Bewertungssystem zum Vergleich von Alteruativen und zur Erstellung von Empfehl ungen,
8. U nterstiitzung bei der Alarmausli:isung.
Die Handhabung dieser Systeme erfordert
1. die Integration der Systeme in eill Gesamtsystem;
2. eille Bellutzeroberfiache, die Illteraktionell wahlweise iiber Meniis, Formulare und Eingabe erlaubt und vom Benutzer selbst bei Bedarf modifiziert werden kann;
3. Kommunikation iiber Maus oder Tastastur;
4. ein Fenstersystem, das es erlaubt, gleichzeitig Ergebnisse darzustellen und mit dem System im Dialog zu arbeiten;
5. efn intelligentes Help-System, das zum einen Grundfehler verhindert und zum anderen Erklarungen zur Benutzung des Systems anbietet.
107
Aus dem Vergleich del' Abb. 2 und del' Liste del' Basissysteme geht hervor, daB ein solches System zur Zeit nicht als geschlossenes System gekauft werden kann. Daher ist es erforderlich, daB bei sol cheri Systemen Schnittstellen geeignet definiert und offengelegt werden, Anderungsmoglichkeiten bestehen und neue Funktionalitaten zugefiigt werden konnen. Grundprinzipien del' modernen Informatik, wie die del' Strukturierung, Modularisierung, Abstraktion und Lokalitat sind zwingend zu beachten. Nur mit ihrer Hilfe ist es moglich, die Komplexitat des Systems zu bewaJtigen.
1m folgenden werden Anforderungen an solch ein System beschrieben, wie sie durch die Verarbeitung von Daten aus dem BundesmeBnetz vorgegeben werden. Dies geschieht in zwei Schritten. Zunachst werden die im Rahmen dieses MeBnetzes anfallenden Daten beschrieben. Die Funktionalitaten sind dem Benutzer in Form von Meniis anzubieten, die mit Standard-Menii-Generatoren erstellt werden konnen. Darauf wird nicht weiter eingegangen.
Daten und Datenmodelle
Daten aus dem BundesmeBnetz Durch die am BundesmeBnetz beteiligten Institutionen werden drei Klassen von Daten an die ZdB iibertragen:
1. MeBdaten bzw. gepriifte MeBwerte entsprechend §2 StrVG (alle Daten sind plausibel ),
2. prognostizierte Daten,
3. Friihwarndaten (z.B. yom Deutschen Wetterdienst - DWD) von Ereignissen auBerhalb del' Grenzen von IMIS.
Die Daten miissen wedel' iiberpriift noch mit Prognosen verglichen werden. Die Ubertragung del' Daten an die ZdB solI grundsatzlich auf Veranlassung des Erzeugers und unmittelbar nach Fertigstellung geschehen.
MeBdaten Alle MeBdaten des BundesmeBnetzes werden in einem einheitlichen Format iibertragen. Charakteristisch fUr die Daten ist:
1. MeBort, gemessene GroBe und MeBverfahren sind je eindeutig vorgegeben.
2. Alle MeBdaten sind mindestens gepriifte Me£werte und miissen keinen weiteren Plausibilitatskontrollen unterzogen werden.
3. MeBdaten aus einem Teilnetz sind fiir weitere Bearbeitung lediglich in plausible und nichtplausible Daten zu unterteilen.
Die Mefldaten des BUlldesmeflnetzes konnell bei del' ZdB dahel' unmittelbal' nach ihl'em Eintrefl'en archiviert werden.
Prognostizierte Daten An Prognosedaten fallen an
108
1. Prognosen etwa der Bundesanstalt fiir Gewasserkunde,
2. Trajektoriendaten des Deutschen Wetterdienstes,
3. Prognose- und Diagnosedaten des Deutschen Wetterdienstes (DWD).
Prognose- und Diagnose-Daten des DWD werden im selben Format iibertragen, da sie aus analogen Modellen abgeleitet werden. Sie unterscheiden sich wie folgt:
A usbreitungsprognosen sind Berechnungen der Ausbreitung anhand prognostizierter meterologischer Parameter.
Ausbreitungsdiagnosen sind Berechnungen der Ausbreitung anhand bekannter metereologischer Parameter.
Trajektorien und Diagnosen konnen auch auf Anfrage durch die ZdB erstellt und iibertragen werden. Ausbreitungsprognosen werdn nur routinemaBig erstellt und nach ihrer Erstellung durch den DWD direkt an die ZdB iibertragen.
Friihwarndaten von Ereignissen auBerhalb des Geltungsbereichs von IMIS Der DWD ist an das GTS (Global Telecommunication System) der WMO (World Meterological Organisation) angeschlossen und erhalt von dort
1. Friihwarnungen iiber Ereignisse mit radiologischer Folge,
2. Ausbreitungsprognosen internationaler Institute.
Bei Eintreffen solcher Warnungen werden diese an die ZdB weitergeleitet. In diesem Fall ist der Intensivbetrieb von IMIS vorzubereiten (Vorwarnstufe).
Datenmodelle Das Datenmodell ist die funktionale Schnittstelle, die das Datenbanksystem dem Benutzer zur Verfiigung stellt und definiert somit alle Moglichkeiten der Strukturierung der Daten und des Zugriffs auf die Daten (Speichern, Andern und Wiederauffinden).
Das Datenmodell soIl es ermoglichen, den betrachteten Ausschnitt der Umwelt moglichst einfach und naturgetreu zu modellieren. Die im letzten Abschnitt genallllten Dateninhalte (MeBdaten, Rechendaten, Progllosedatell) sind in ihrer Struktur breit gefiichert, so daB kein einheitliches Datenmodell alle Allwendungen gleich gut ullterst.iit.zen kallll. Insbesondere lassen sich die Rechenergebnisse (Matrizen und Vektoren), die Beschreibungen (Texte) sowie die geographischen Dat.en nur widerstrebend und in sehr unnatiirlicher Weise auf die Attribute und Tabellen eines relationen Dat.enballksystems abbilden. Erweiterungen am relationalen Datenmodell durch Int.egration der Datenstrukturell Text und Langes-Feld erscheinen als Notli:ist1ng, die das Modell an Anwendungsbereiche anpassen soIl, fUr die es zunachst nicht geeignet ist.
109
Objektorientierte Datenbanksysteme konnen dagegen durch ihre abstrakten Datentypen stark divergierende AnwendungsfaIle wesentlich besser integrieren: Fur jeden der genannten Dateninhalte wird ein eigener abstrakter Datentyp definiert, dessen Operatoren auf einer jeweils fUr sie zugeschnittenen und nur ihnen bekannten Datenstruktur arbeiten. Beispiele fUr abstrakte Datentypen sind u.a. Wert, Text, Vektor, Matrix, Polygon, Klasse. (Der Datentyp Vektor eignet sich zur Beschreibung zeitlicher Verlaufe. Der Datentyp Matrix ist Grundlage der DarsteIlung von Lagekarten. Der Datentyp Klasse erlaubt es, MeBdaten zu erfassen und Selektionen vorzunehmen.) Die Operatoren der ADT-Module bieten dem Anwendungsprogrammierer die dem Problem jeweils am besten angepaBte Sicht auf seine Daten, bzw. Operationen zur Manipulation seiner Daten. Die innere Struktur der Datenobjekte sollte so frei definierbar sein, daB auch existierende Software-Werkzeuge mit geringem Aufwand als Operatoren eines neuen ADT -Moduls integriert werden konnen.
Fur den Prototyp stand mit RSYST [6] eine Datenbasis zur VerfUgung, die eine hierarchische Strukturierung von Datenobjekten beliebiger interner Struktur erlaubt.
Strukturierung der Daten AIle gemessenen und gerechneten Daten konnen durch eine Reihe von Parametern unterschieden werden. Darunter fallen
Quelle der Daten (Bund, Lander, Rechner);
Zeitpunkt der Erstellung (Jahr, Monat, Tag);
Art der Daten (MeBnetz, Mefistelle, Mefimethode, Art der Messung).
Diese Parameter konnen zur hierarchischen Strukturierung der Daten verwendet werden. Abb. 3 zeigt einen Vorschlag fUr solch eine Strukturierung der Basisdaten.
Funktionalitaten
Beim Arbeiten mit IMIS sollte zwischen zwei Stati unterschieden werden:
Status 1 Systemmanager oder DB-Administrator,
Status 2 Systemanwender.
Die Hauptaufgaben des Systemmangers sind
• Verwaltung der Strukturen der DB (Archivierung von Daten, Strukturierung von Daten, Festlegung von Datenobjektbeschreibungen);
• Aktualisierung der Daten (Verwaltung der Stammclaten, Verwaltung cler Verarbei tungsregeln, Verwalt ung der Anwencler berechtigungen).
110
Del' Systemmanager hat alle Privilegien des Systems. Er kann Teilaufgaben und die zugeharigen Privilegien an Teilsystemmanager delegieren. Die Aufgabenverteilung soU anharid einer besonders zu schiitzenden ID erkannt werden.
Systemanwender sind alle, die mit del' DarsteUung, Auswertung und Beurteilung von Daten zu tun haben. Entsprechend ihres Aufgabenbereichs sollen sie nul' zu Teilen del' DB Zugriff haben.
Folgende Funktionalitaten sind auf die aus dem BundesmeBnetz an die ZdB iibertragenen Daten durchzufiihren: Archivierung, Selektion, Statistische Auswertung, Darstellung, Transformationen, Modifikationen, Vergleiche.
Alle Funktionalitaten werden routinemaf3ig nur auf plausible Daten durchgefUhrt. Sollen auch als nichtplausibel eingestufte Daten behandelt werden, so ist dies yom Benutzer anzugeben. Einmal getroffene Wahlen sind solange giiltig, bis der Benutzer sie explizit andert.
Archivierung
Die im Rahmen des BundesmeBnetzes zur ZdB iibertragenen Daten kannen direkt archiviert werden. Die Archivierung sonte strukturiert erfolgen. Ein Beispiel fiir eine magliche Struktur zeigt Abb. 2. MeBsatz meint dann die iibertragenen Daten. Konfliktieren sie mit del' Struktur, miissen sie aufgespalten werden. Es ist notwendig, zu priifen, ob alle erwarteten Daten angekommen sind. Nicht eingetroffene Datensatze soUten iiber elektronic mail angefordert werden.
Das DB-System zur Archivierung der Daten soUte mit dem zur Verarbeitung der Daten identisch sein, so daB Daten aus dem Archivierungssystem leicht iibertragen werden kannen (siehe auch Transformation). Bei Verwendung optischer Disks entfallen aufwendige Daten-Refresh-Mechanismen.
Selektionen
Selektionen sollen es erlauben, Daten aus raumlichen und zeitlichen Fenstern auszuwahlen, weiterzuverarbeiten oder an andere Partner im IMIS weiterzureichen.
Magliche Selektionen im Rahmen del' hier behandelten Daten sind Daten fUr Lander, Zeitfenster fUr Auswertung von Verlaufen, Auswahl nach vorgegebenen Kriterien (SQL) oder Zusammenfassung von Daten aus verschiedenen MeBnetzen.
Durch die Funktionalitat Selektion solI eine neue Teilstruktur aufgebaut werden, die als Blatter gerade die selektierten Datensatze ellthalt. Fiir den Fall del' Auswahl raumlicher Teilmengen aus del' Arbeitsdatenbank sollen zusatzlich auch reduzierte Datenobjekte generiert werden.
111
Statistische Auswertung
Die Aufgaben des Statistikpaketes sind zeitliehe und raumliehe Aggregierung von Daten, Angabe von Streuungen, Bereehnung von Korrelationen raumlieh und zeitlieh verteilter GraBen (eventuell unter Beriieksiehtigung einer zeitliehen Versehiebung), Bereehnung von Regressionen, Durehfiihrung einfaeher Trelldanalysen, Verteilungsuntersuehungen.
Die Funktionen des Statistikpaketes sollen auf aIle Datensatze del' ausgewahlten Datenbank angewandt werden, solange sie von derselben Art sind. Es ist deswegen natig, ein verfiigbares Statistikpaket dureh entspreehende Sehnittstellen zu DB und Benutzeroberftaehe in das Gesamtsystem zu integrieren.
Darstellung
Die Anforderungen an das Graphiksystem sind:
1. Das Graphiksystem solI hardware-unabhangig sein, d.h. auf einem Standard wie etwa GKS aufsetzen kannen.
2. Das Graphiksystem muB Zugang zur Datenbank haben, urn zum einen die dort abgelegten Strukturen veransehaulichen zu kannen und zum anderen die Generierung neuer Strukturen zu unterstiitzen.
3. Das Graphiksystem sollte objektorientiert arbeiten. Ein graphisehes Objekt sollte dureh seinen N amen gekennzeiehnet werden. Del' Name sollte indizierbar sein.
4. Grundfunktionen sollten benutzerfreundlieh angeboten werden.
5. Die Vorauswahl del' Grundfunktionen sollte in Abhangigkeit del' darzustellenden Werte gesehehen.
6. Folgende Grunddarstellungen sollten verfiigbar sein:
• Grenzen (etwa Bundesgebiet, Lander, Regierungsbezirke, Kreise);
• MeBstellen (etwa Orte, Werte numeriseh, Werte graphiseh);
• Interpolationen (etwa N etzaufbau anhand aktueller MeBdaten, Netzinterpolation, Kreisinterpolation);
• erreehnete und prognostizierte Daten;
• Verlaufe (etwa zum Vergleieh gleieher odeI' versehiedener GraBen).
7. Die Software sollte auBerdem folgende Operationen erlauben:
• Kombinationen graphiseher Objekte und daraus erzeugter Darstellungen;
• Vergleiche etwa von errechneten und gemessenen Daten;
• Standardprozeduren zur Generierung von Standardbildern.
112
Transformationen
Das Arbeiten mit den Daten sollte in der Regel nicht liber die Archivierungsdatenbank, sondern im Rahmen temporarer Arbeitsdatenbanken geschehen. Dazu ist es notwendig, die Daten zu transformieren. Bei der Transformation soIl die Grundstruktur erhalten bleiben. Die MeBsatze sollen aber den Arbeitsanforderungen angepaBt werden. Es ist sinnvoll, in einem Blatt Daten aller MeBstellen zusammenzufassen (Lagedarstellungen).
Diese Daten sollen auch als Grundlage der Ubermittlung von Daten anderer Partner dienen. Das Selektionsmenli soll deswegen auch auf diese DB anwendbar sein.
Modifikation
Die Aufgaben der Funktionalitat Modifikation sind
Korrigieren einzelner Werte in Datenobjekten und
Manipulation von Netzen, Feldern und Vektoren zwecks einheitlicher Darstellung.
Zur ErfUllung der zweiten Aufgabe ist ein Paket mathematischer Grundoperationen auf diese Datenobjekte bereitzustellen.
Die Anderung von Werten in einzelnen Datenobjekten muB in der DB reproduzierbar dokumentiert werden.
Vergleiche
Vergleiche zwischen verschiedenen Arten sollen moglich sein
• numerisch - etwa durch Angabe von Tabellen, Difl'erenzen oder Verhaltnissen;
• statistisch - etwa durch Angabe von Korrelationen, Regressionen und Verteilungen;
• graphisch - etwa durch Angabe von Lagebildern, Bildern mit Verteilungen an MeBstellen, Histogrammen, Verlaufe von verschiedenen GraBen.
Die Aufgabe der Funktionalitat Vergleiche ist es, die entsprechenden Daten bereitzustellen, eventuelle Anpassungen wie zeitliche oder raumliche Normierungen durchzufUhren und die gewahlte Darstellung zu initialisieren.
Der IKE-Prototyp zur Real-Time-Analyse von Umweltdaten
Zum Nachweis der Realisierbarkeit des beschriebenen Konzeptes haben wir einen Prototyp fUr eine Umwelt-Workstation entwickelt, mit dessen Hilfe wir ForderUllgen an Software und Hardware von Umwelt-Informationssystemen untersuchen uncI Teilprodukte realer Systeme auf ihre Funktionalitat iiberprlifen konnen. N aturgemaB fehlen
113
dem Prototyp wichtige Eigenschaften eines endgiiltigen Systems. Insbesondere erfiillen seine Benutzeroberfliiche und damit die Bedienbarkeit und die Zahl der angeschlossenen Simulationsprogramme nicht alle Forderungen. Trotzdem gibt er eine gute Basis zur Veranschaulichung unserer Vorstellungen und zur Diskussion alternativer Losungen.
Die Konfigurierung dieses Prototyps ist in Abb. 4 dargestellt.
Begleitend zum Vortrag sollen auf diesem Prototyp folgende Schritte gezeigt werden:
1. Ubermittlung von Rohdaten aus MeBnetz nach Eliminierung von AusreiJlern
2. Lagedarstellung mittels der Rohdaten
3. Darstellung der aktuellen Lage
4. Ausblenden eines Datensatzes fiir ein Bundesland
5. Vergleich der aktuellen Lage mit Daten eines Ereignisfalles
6. Simulation eines Ereignisfalles
7. Darstellung des Quell-Terms
8. Darstellung des Windfeldes fiir Ausbreitungsrechnungen
9. Durchfiihrung von Ausbreitungsrechnungen auf der CRAY-2 in Stuttgart und Erstellen einer Prognose
10. Darstellung des zeitlichen Verlaufs der Prognosedaten
Literatur
[1] Biihling, A.j et al.: Integrated Measurement and Information System for Surveillance of Environmental Radioactivity in F.R.G. Nuclear Europe 3-4 (1988), p. 22
(2) Schmidt, F.; et al.: Konzepte zur DV-gestiitzten Uberwachung der Umweltradioaktivitat im Rahmen des Integrierten MeJl- und Informationssystems (IMIS) des BMU. Stuttgart: IKE, Oktober 1988 (IKE 4-128)
[3] Miiller, H.; et al.: Entwicklung und Einsatz verbesserter, zeitabh8.ngiger Modelle zur Berechnung der potentiellen Strahlenexposition nach Storf8.llen. AbschluBbericht zu Forschungsvorhaben St.Sch. 798 des BMI (1986)
[4] Masamichi Chino; Hirohiko Ishikawa; Michiaki Kai: SPEEDI, Manual of a Suite of Computer Codes, Code System for Real-time Prediction of Radiation Dose to the Public due to an Accidental Release from a Nuclear Power Plant. RSIC Computer Code Collection, CCC-507 (1987)
(5) Schmidt, F.; Tischendorf, M.: Moglichkeiten der Kontrolle und Analyse von Umweltdaten durch Kopplung von Datenbank- und Expertensystemtechniken. 4. Symposium 'Informatik im Umweltschutz', Karlsruhe, November 1989
[6] Riihle, R.: RSYST als Software-Basis eines einheitlichen Systems {iir 'ComputerAided Aerospace Analysis and Design'. Stuttgart. RUS, 1989, Studie.
N .... Q)
C en en Q)
E en Q)
"0 C :J
CD
Luft
Orts-dosis-leistung
Flusse
Nord/ Ostsee
Modell-annahmen
114
BMU
- ZdB -'---
Leitstellen
Abb. 1: Struktur von IMIS
Wissensbasen
Messwerte
/ ~
Regel n
r o ::J C. Cb en C. o "" Cb ::J N Cb :J
"" ., o Cb :J
und Beziehungen
Wissensbeschaffung Wissensverarbeitung w5hrend durch
Planung und Betrieb Simulation und Inferenz
Ben utzero berflache
Optionen Berichts Konsitenz Karten Daten Erklorungs Generierer Editor PrO fer Generierer Verdichter Komponente
Abb. 2: IMIS als integriertes Syst,em
115
Abb. 3: Basisstruktur fiir Daten in del' IMIS-Datenbank
Men·
netze
z.B.
IMIS
elrtended LAN, z.B. BelWii
116
OEC: Simulation·Output
S02·Emiss lonen der Haushalte und Klelnverbraucher
24 .3.1 2 .00-3.00 Uhr
Workstat ion
Abb.4:
CRAY·2: Sim.
Simulalion Deposition NeU,. SO ,,50 Hbhe,. 20 Ze ll ,. 24
Breitband Netz, z.B. Ethernet Vektor.
Rechner
z.B .
CRAV·2
IKE-Prototyp eines Arbeitsplatzes zur Analyse von Umweltdaten unter Real-Time-Bedingungen
Informatik - Fachberichte
Band 150: J. Halin (Hrsg.), Simulationstechnik. 4. Symposium, ZOrich, September 1987. Proceedings. XIV, 690 Seiten. 1987.
Band 151: E. Buchberger, J. Retti (Hrsg.), 3. Osterreichische Arti-ficial-Intelligence-Tagung. Wien, September 1987. Proceedings. VIII, 181 Seiten. 1987.
Band 152: K. Marik (Ed.), GWAI-87. 11th German Workshop on Artificial Intelligence. Geseke, SeptlOkl. 1987. Proceedings. XI, 405 Seiten. 1987.
Band 153: D. Meyer-Ebrecht (Hrsg.), ASST'87. 6. Aachener Sym-posium far Signaltheorie. Aachen, September.1987. Proceedings. XII, 390 Seiten. 1987. Band 154: U. Herzog, M. Paterok (Hrsg.), Messung, Modellierung und Bewertung von Rechensystemen. 4. GI/ITG-Fachtagung, Er-langen, SeptlOkl.1987. Proceedings. XI, 388 Seiten. 1987.
Band 155: W. Brauer, W. Wahlster (Hrsg.), Wissensbasierte Sy-steme. 2. Internationaler GI-KongreB, MOnchen, Oklober 1987. XIV, 432 Seiten. 1987.
Band 156: M. Paul (Hrsg.), GI-17. Jahrestagung. Computerinte-grierter Arbeitsplatz im BOro. MOnchen, Oktober 1987. Proceed-ings. XIII, 934 Seiten. 1987.
Band 157: U. Mahn, Attributierte Grammatiken und Attributierungs-algorithm en. IX, 272 Seiten. 1988.
Band 158: G. Cyranek, A. Kachru, H. Kaiser (Hrsg.), Informatik und "Dritte Welt". X, 302 Seiten. 1988. Band 159: Th. Christaller, H.-W. Hein, M. M. Richter (Hrsg.), KOnst-Iiche Intelligenz. FrOhjahrsschulen, Dassel, 1985 und 1986. VII, 342 Seiten. 1988.
Band 160: H. Miincher, Fehlertolerante dezentrale ProzeBauto-matisierung. XVI, 243 Seiten. 1987.
Band 161: P. Peinl, Synchronisation in zentralisierten Datenbank-systemen. XII, 227 Seiten. 1987.
Band 162: H. Stoyan (Hrsg.), BegrOndungsverwaltung. Proceed-ings, 1986. VII, 153 Seiten. 1988.
Band 163: H. Maller, Realistische Computergraphik. VII, 146 Seiten. 1988.
Band 164: M. Eulenstein, Generierung portabler Compiler. X, 235 Seiten. 1988.
Band 165: H.-U. Heia, Oberlast in Rechensystemen. IX, 176 Seiten. 1988.
Band 166: K. HOrmann, Kollisionsfreie Bahnen far Industrieroboter. XII, 157 Seiten. 1988.
Band 167: R. Lauber (Hrsg.), ProzeBrechensysteme '88. Stuttgart, Miirz 1988. Proceedings. XIV, 799 Seiten. 1988.
Band 168: U. Kastens, F. J. Rammig (Hrsg.), Architeklur und Betrieb von Rechensystemen. 10. GI/ITG-Fachtagung, Paderborn, Miirz 1988. Proceedings. IX, 405 Seiten. 1988.
Band 169: G. Heyer, J. Krems, G. GOrz (Hrsg.), Wissensarten und ihre Darstellung. VIII, 292 Selten. 1988.
Band 170: A. Jaeschke, B. Page (Hrsg.), Informatikanwendungen im Umweltbereich. 2. Symposium, Karlsruhe, 1987. Proceedings. X, 201 Seiten. 1988.
Band 171: H. Lutterbach (Hrsg.), Non-Standard Datenbanken far Anwendungen der Graphischen Datenverarbeitung. GI-Fachge-spriich, Dortmund, Milrz 1988, Proceedings. VII, 183 Seiten. 1988.
Band 172: G. Rahmstorf (Hrsg.), Wissensreprllsentation in Ex-pertensystemen. Workshop, Herrenberg, Milrz 1987. Proceedings. VII, 189 Seiten. 1988.
Band 173: M. H. Schulz, Testmustergenerierung und Fehlersimu-lation in digitalen Schaltungen mit hoher Komplexitat. IX, 165 Sei-ten. 1988.
Band 174: A. EndrOs, Rechtsprechung und Computer in den neun-ziger Jahren. XIX, 129 Seiten. 1988.
Band 175: J. HOlsemann, Funklioneller Test der·AuflOsung von Zu-griffskonfliklen in Mehrrechnersystemen. X, 179 Seiten. 1988.
Band 176: H. Trost (Hrsg.), 4.0sterreichische Artificial-Intelli-gence-Tagung. Wien, August 1988. Proceedings. VIII, 207 Seiten. 1988.
Band 177: L. Voelkel, J. Pliquett, Signaturanalyse. 223 Seiten.1989.
Band 178: H. GOttler, Graphgrammatiken in der Softwaretechnik. VIII, 244 Seiten. 1988.
Band 179: W.Ameling (Hrsg.), Simulationstechnik. 5. Symposium. Aachen, September 1988. Proceedings. XIV, 538 Seiten. 1988.
Band 180: H. Bunke, O. KObler, P. Stucki (Hrsg.), Mustererkennung 1988. 10. DAGM-Symposium, ZOrich, September 1988. Proceed-ings. XV, 361 Seiten. 1988.
Band 181: W. Hoeppner (Hrsg.), KOnstliche Intelligenz. GWAI-88, 12. Jahrestagung. Eringerfeld, September 1988. Proceedings. XII, 333 Seiten. 1988.
Band 182: W. Barth (Hrsg.), Visualisierungstechniken und Algorith-men. Fachgesprilch, Wien, September 1988. Proceedings. VIII, 247 Seiten. 1988.
Band 183: A. Clauer, W. Purgathofer (Hrsg.), AUSTROGRAPHICS '88. Fachtagung, Wien, September 1988. Proceedings. VIII, 267 Seiten. 1988.
Band 184: B. Gollan, W. Paul, A. Schmitt (Hrsg.), Innovative Infor-mations-Infrastrukluren. 1.1.1. - Forum, SaarbrOcken, Oklober 1988. Proceedings. VIII, 291 Seiten. 1988.
Band 185: B. Mitschang, Ein MolekOI-Atom-Datenmodell fOr Non-Standard-Anwendungen. XI, 230 Seiten. 1988.
Band 186: E. Rahm, Synchronisation in Mehrrechner-Datenbank-systemen. IX, 272 Seiten. 1988.
Band 187: R. Valk (Hrsg.), GI-18. Jahrestagung I. Vernetzte und komplexe Informatik-Systeme. Hamburg, Oklober 1988. Proceed-ings. XVI, 776 Seiten.
Band 188: R. Valk (Hrsg.), GI-18. Jahrestagung II. Vernetzte und komplexe Informatik-Systeme. Hamburg, Oklober 1988. Proceed-ings. XVI, 704 Seiten.
Band 189: B. Wolfinger (Hrsg.), Vernetzte und komplexe Informatik-Systeme. Industrieprogramm zur 18. Jahrestagung der GI, Ham-burg, Oklober 1988. Proceedings. X, 229 Seiten. 1988.
Band 190: D. Maurer, Relevanzanalyse. VIII, 239 Seiten. 1988.
Band 191: P. Levi, Planen far autonome Montageroboter. XIII, 259 Seiten. 1988.
Band 192: K. Kansy, P.WiBkirchen (Hrsg.), Graphik im BOrobe-reich. Proceedings, 1988. VIII, 187 Seiten. 1988.
Band 193: W. Gotthard, Datenbanksysteme fOr Software-Prod uk-tionsumgebungen. X, 193 Seiten. 1988.
Band 194: C. Lewerentz, Interaklives Entwerfen groBer Programm-systeme. VII, 179 Seiten. 1988.
Band 195: I. S. Batori, U. Hahn, M. Pinkal, W. Wahlster (Hrsg.), Com-puterlinguistik und ihre theoretischen Grundlagen. Proceedings. IX, 218 Seiten. 1988.
Band 197: M. Leszak, H. Eggert, Petri-Netz-Methoden und -Werk-zeuge. XII, 254 Seiten. 1989.
Band 198: U. Reimer, FRM: Ein Frame-Reprilsentationsmodell und seine formale Semantik. VIII, 161 Seiten. 1988.
Band 199: C. Beckstein, Zur Logik der Logik-Programmierung. IX, 246 Seiten. 1988.
Informatik - Fachberichte
Band 200: A. Reinefeld, Spielbaum-Suchverfahren. IX, 191 Seiten. 1989. Band 201: A. M. Kotz, Triggermechanismen in Datenbanksystemen. VIII, 187 Seiten. 1989.
Band 202: Th. Christaller (Hrsg.), Kanstliche Intelligenz. 5. Frah-jahrsschule, KIFS-87, GOnne, MarzlApril1987. Proceedings. VII, 403 Seiten. 1989.
Band 203: K. v. Luck (Hrsg.), KOnstliche Intelligenz. 7. FrOhjahrs-schule, KIFS-89, GOnne, Marz 1989. Proceedings. VII, 302 Seiten. 1989. Band 204: T. Harder (Hrsg.), Datenbanksysteme in BOro, Technik und Wissenschaft. GI/SI-Fachtagung, ZOrich, Marz 1989. Pro-ceedings. XII, 427 Seiten. 1989. Band 205: P. J. KOhn (Hrsg.), Kommunikation in verteilten Sy-stemen. ITG/GI-Fachlagung, Stuttgart, Februar 1989. Proceed-ings. XII, 907 Seiten. 1989. Band 206: P. Horster, H.lsselhorst, Approximative Public-Key-Kryptosysteme. VII, 174 Seiten. 1989. Band 207: J. Knop (Hrsg.), Organisation der Datenverarbeitung an der Schwelle der 90er Jahre. 8. GI-Fachgesprach, DOsseldorf, Marz 1989. Proceedings. IX, 276 Seiten. 1989. Band 208: J. Retti, K. Leidlmair (Hrsg.), 5. Osterreichische Artificial-Intelligence-Tagung, IglslTirol, Marz 1989. Proceedings. XI, 452 Seiten. 1989.
Band 209: U. W. Lipeck, Dynamische Integritat von Datenbanken. VIII, 140 Seiten. 1989.
Band 210: K. Drosten, Termersetzungssysteme. IX, 152 Seiten. 1989. Band 211: H. W. Meuer (Hrsg.), SUPERCOMPUTER '89. Mannheim, Juni 1989. Proceedings, 1989. VIII, 171 Seiten. 1989. Band 212:W.-M. Lippe (Hrsg.), Software-Entwicklung. Fachtagung, Marburg, Juni 1989. Proceedings. IX, 290 Seiten. 1989. Band 213: I. Walter, DatenbankgestOtzte Reprasentation und Ex-traktion von Episodenbeschreibungen aus Bildfolgen. VIII, 243 Seiten. 1989. Band 214: W. GOrke, H. SOrensen (Hrsg.), Fehlertolerierende Rechensysteme I Fault-Tolerant Computing Systems. 4.lnternatio-nale GIlITG/GMA-Fachtagung, Baden-Baden, September 1989. Proceedings. XI, 390 Seiten. 1989. Band 215: M. Bidjan-Irani, Qualitat und Testbarkeit hochinte-grierter Schaltungen. IX, 169 Seiten. 1989. Band 216: D. Metzing (Hrsg.), GWAI-89. 13th German Workshop on Artificial Intelligence. Eringerfeld, September 1989. Proceed-ings. XII, 485 Seiten. 1989. Band 217: M. Zieher, Kopplung von Rechnenoetzen. XII, 218 Seiten. 1989. Band 218: G. Stiege, J. S. Lie (Hrsg.), Messung, Modellierung und Bewertung von Rechensystemen und Netzen. 5. GIIITG-Fachta-gung, Braunschweig, September 1989. Proceedings. IX, 342 Seiten. 1989. Band 219: H. Burkhardt, K. H. HOhne, B. Neumann (Hrsg.), Muster-erkennung 1989. 11. DAGM-Symposium, Hamburg, Oktober 1989. Proceedings. XIX, 575 Seiten. 1989 Band 220: F. Stetter, W. Brauer (Hrsg.), Informatik und Schule 1989: Zukunftsperspektiven der Informatik far Schule und Ausbildung. GI-Fachtagung, Manchen, November 1989. Proceedings. XI, 359 Seiten. 1989. Band 221: H. Schelhowe (Hrsg.), Frauenwelt - Computerraume. GI-Fachtagung, Bremen, September 1989. Proceedings. XV, 284 Seiten. 1989. Band 222: M. Paul (Hrsg.), GI-19. Jahrestagung I. MOnchen, Oktober 1989. Proceedings. XVI, 717 Seiten. 1989.
Band 223: M. Paul (Hrsg.), GI-19. Jahrestagung II. Manchen, Oktober 1989. Proceedings. XVI, 719 Seiten. 1989.
Band 224: U. Voges, Software-Diversitat und ihre Modellierung. VIII, 211 Seiten.1989
Band 225: W. Stoll, Test von OSI-Protokollen. IX, 205 Seiten. 1989.
Band 226: F. Mattern, Verteilte Basisalgorithmen. IX, 285 Seiten. 1989.
Band 227: W. Brauer, C. Freksa (Hrsg.), Wissensbasierte Systeme. 3. Internationaler GI-KongreB, MOnchen, Oktober 1989. Proceed-ings. X, 544 Seiten. 1989.
Band 228: A. Jaeschke, W. Geiger, B. Page (Hrsg.), Informatik im Umweltschutz. 4. Symposium, Karlsruhe, November 1989. Proceedings. XII, 452 Seiten. 1989.
Band 229: W. Coy, L. Bonsiepen, Erfahrung und Berechnung. Kritik der Expertensystemtechnik. VII, 209 Seiten. 1989.
Band 230: A. Bode, R. Dierstein, M. GObel, A. Jaeschke (Hrsg.), Visualisierung von Umweltdaten in Supercomputersystemen. Karlsruhe, November 1989. Proceedings, 1989. XII, 116 Seiten. 1990.
Band 231: R. Henn, K. Stieger (Hrsg.), PEARL 89 - Workshop aber Realzeitsysteme. 10. Fachtagung, Boppard, Dezember 1989. Pro-ceedings. X, 243 Seiten. 1989.
Band 232: R. Loogen, Parallele Implementierung funktionaler Pro-grammiersprachen. IX, 385 Seiten. 1990.
Band 233: S. Jablonski, Datenverwaltung in verteilten Systemen. XIII, 336 Seiten. 1990.
Band 234: A. Pfitzmann, Diensteintegrierende Kommunikations-netze mit teilnehmerOberprOfbarem Datenschutz. XII, 343 Seiten. 1990.
Band 235: C. Feder, Ausnahmebehandlung in objektorientierten Programmiersprachen. IX, 250 Seiten. 1990.
Band 236: J. Stoll, Fehlertoleranz in verteilten Realzeitsystemen. IX, 200 Seiten. 1990.
Band 237: R. Grebe (Hrsg.), Parallele Datenverarbeitung mit dem Transputer. Aachen, September 1989. Proceedings, 1989. VIII, 241 Seiten. 1990.
Band 238: B. Endres-Niggemeyer, T. Hermann, A. Kobsa, D. ROs-ner (Hrsg.), Interaktion und Kommunikation mit dem Computer. Ulm, Miirz 1989. Proceedings, 1989. VIII, 175 Seiten. 1990.
Band 239: K. Kansy, P. WiBkirchen (Hrsg.), Graphik und KI. KOnigs-winter, April 1990. Proceedings, 1990. VII, 125 Seiten. 1990.
Band 240: D. Tavangarian, Flagorientierte Assoziativspeicher und -prozessoren. XII. 193 Seiten. 1990.
Band 241: A. Schill, Migrationssteuerung und Konfigurationsver-waltung fOr verteilte objektorientierte Anwendungen.IX, 174 Seiten. 1990.
Band 242: D. Wybranietz, Multicast-Kommunikation in verteilten Systemen. VIII, 191 Seiten. 1990.
Band 244: B. R. Kammerer, Sprecherunabhangigkeit und Spre-cheradaption. VIII, 110 Seiten. 1990.
Band 246: Th. Braunl, Massiv parallele Programmierung mit dem Parallaxis-Modell. XII,.168 Seiten. 1990
Band 247: H. Krumm, Funktionelle Analyse von Kommunikations-protokollen. IX, 122 Seiten. 1990.
Band 250: H. W. Meuer (Hrsg.), SUPERCOMPUTER'90. Mannheim, Juni 1990. Proceedings, 1990. VIII, 209 Seiten. 1990.