imported>Leyo: redundant zu DOI

2026-04-08T19:10:39Z

redundant zu DOI

Neue Seite

Der Begriff '''Diskrete-Elemente-Methode''' (engl. {{lang|en|discrete element method}}, '''DEM''') wird heutzutage für zwei numerische Berechnungsverfahren verwendet.

Die häufigste Verwendung findet die von Cundall<ref>P. A. Cundall: ''A computer model for simulating progressive large scale movements in blocky rock systems.'' In: ''Proceedings Symposium Int. Soc. Rock Mech.'' Nancy Metz, vol. 1, 1971, Paper II–8.</ref> im Jahre 1971 entwickelte [[Numerik|numerische]] Berechnungsmethode, mit der die Bewegung einer großen Zahl von Teilchen berechnet werden kann. Die Methode wird manchmal auch als {{lang|en|''Distinct Element Method''}} bezeichnet. Typischerweise erweitert die Diskrete-Elemente-Methode ''[[Molekulardynamik]]''-Simulationen um [[Rotationsfreiheitsgrad]]e. Seit ihrer Einführung hat sich ihr Einsatzgebiet ausgedehnt, wie z. B. auf die Simulation aus der Partikelverfahrenstechnik, der Geotechnik und des Maschinenbaus. Eine Weiterentwicklung des Verfahrens ist die [[erweiterte Diskrete-Elemente-Methode]].

Auf der anderen Seite wird der Begriff DEM auch für ein Stabgittermodell verwendet. Diese Betrachtungsweise – Abbildung eines Körpers durch Stäbe – geht auf Arbeiten von E. G. Kirsch<ref>E. G. Kirsch: ''Die Fundamentalgleichungen der Theorie der Elastizität fester Körper, hergeleitet aus der Betrachtung eines Systems von Punkten, welche durch elastische Streben verbunden sind.'' In: ''Zeitschrift des Vereines Deutscher Ingenieure. ''Band 7, Heft 8, 1868, S. 481–487, 553–570, 631–638.</ref> aus dem Jahr 1868 zurück und wurde zum Beispiel von [[Felix Klein (Mathematiker)|Felix Klein]] und [[Karl Wieghardt (Physiker)|Karl Wieghardt]] Anfang des 20. Jahrhunderts weiterentwickelt. Heute wird die Stabgittermethode unter anderem zur Simulation des Materialverhaltens von Verbundwerkstoffen, insbesondere Gewebestrukturen, eingesetzt.<ref>F. K. Wittel: ''Diskrete Elemente – Modelle zur Bestimmung der Festigkeitsevolution in Verbundwerkstoffen.'' Universität Stuttgart, Fakultät für Luft- und Raumfahrttechnik und Geodäsie, Institut für Statik und Dynamik der Luft- und Raumfahrtkonstruktionen, 2006.</ref><ref>D. Ballhause: ''Diskrete Modellierung des Verformungs- und Versagensverhaltens von Gewebemembranen.'' Universität Stuttgart, Fakultät für Luft- und Raumfahrttechnik und Geodäsie, Institut für Statik und Dynamik der Luft- und Raumfahrtkonstruktionen, 2006.</ref> Darüber hinaus zeigte sich, dass die DEM sich zur Lebensdauerabschätzung von metallischen und keramischen Werkstoffen eignet.<ref>M. Hahn, M. Bouriga, B.-H. Kröplin, T. Wallmersperger: ''Life time prediction of metallic materials with the Discrete-Element-Method.'' In: ''Computational Materials Science.'' Band 71, 2013, S. 146–156.</ref><ref>M. Hahn: ''Lebensdauerabschätzung von metallischen Strukturen mittels der Diskrete-Elemente-Methode im gekoppelten thermo-mechanischen Feld.'' Universität Stuttgart, Fakultät für Luft- und Raumfahrttechnik und Geodäsie, Institut für Statik und Dynamik der Luft- und Raumfahrtkonstruktionen, 2012.</ref><ref>{{Webarchiv |url=http://www.isd.uni-stuttgart.de/forschung/fsp/sum/diskrete_elemente_sum/index.html |wayback=20140917020939 |text=Diskrete-Elemente-Methode, Universität Stuttgart}}.</ref> DEM wird auch zur Beschreibung des Materialverhaltens in der Geomechanik angewandt.

Von der DEM zu unterscheiden ist die [[Finite-Elemente-Methode]] (FEM), ein numerisches Verfahren, das für physikalische Aufgabenstellungen, zum Beispiel bei Festkörpern, angewendet wird.

Der folgende Text beschränkt sich auf die DEM nach Cundall, da sie derzeit eine höhere Relevanz in der Forschung erfährt.

== Verfahren ==
=== Anwendungsbereiche ===
Die Grundannahme des Verfahrens beruht darauf, dass die zu berechnende [[Materie (Physik)]] sich aus einzelnen, abgeschlossenen Elementen zusammensetzt. Diese Elemente können unterschiedliche Formen und Eigenschaften haben.

Anwendung findet die Methode in folgenden Bereichen: [[Numerische Simulation|Simulation]] des Verhaltens von [[Atom]]en und [[Molekül]]en (siehe auch chemische [[Moleküldynamik]]), Entwicklung von Verarbeitungsmethoden für [[Schüttgut|Schüttgüter]] in [[Silo]]s, z. B. [[Getreide]], Entwicklung von Verarbeitungsmethoden für [[Massenrohstoff]]e, z. B. [[Sand]], Simulation von [[Staub]]entwicklung, z. B. [[Toner]] sowie Simulation von Prozessen der [[Geodynamik]], z. B. Dynamik von [[Akkretionskeil]]en.
=== Ablauf ===
Bei einer DEM-[[Simulation]] werden alle Teilchen in einer bestimmten Startgeometrie positioniert und mit einer Anfangs[[geschwindigkeit]] versehen. Aus diesen Anfangsdaten und den physikalischen Gesetzen, die für die Teilchen relevant sind, werden die Kräfte ausgerechnet, die auf jedes Teilchen wirken.

Kräfte, die hier in Frage kommen, sind zum Beispiel im [[makroskopisch]]en Fall: [[Reibung]]skräfte, wenn zwei Teilchen einander streifen, Rückstoßende Kräfte, wenn zwei Teilchen aufeinander treffen und dabei leicht reversibel deformiert werden, [[Gravitation]]skräfte, also die Anziehung der Teilchen aufgrund ihrer Massen (nur relevant bei astronomischen Simulationen).

Oder auf [[Molekül|molekularer]] Ebene: [[Coulombkraft|Coulomb-Kräfte]], also die [[Elektrostatik|elektrostatische]] Anziehung oder Abstoßung der Teilchen, falls diese eine [[elektrische Ladung]] tragen, [[Austauschwechselwirkung|Pauli Repulsion]], wenn zwei Atome nahe aneinandergeraten sowie [[Van-der-Waals-Kräfte]].

Alle diese Kräfte werden aufsummiert und danach mit Hilfe eines numerischen Integrationsverfahren aus der [[Isaac Newton|Newtonschen]] Bewegungsgleichung die Veränderung der Teilchengeschwindigkeit und -position berechnet, die sich in einem gewissen Zeitschritt ergibt. Danach werden mit den veränderten Positionen und Geschwindigkeiten erneut die Kräfte berechnet und diese [[Schleife (Programmierung)|Schleife]] so lange wiederholt, bis der Simulationszeitraum beendet ist.

=== Langreichweitige Kräfte ===
Wenn langreichweitige Kräfte (typischerweise Gravitationskräfte oder elektrostatische Kräfte) berücksichtigt werden, so muss grundsätzlich die Wechselwirkung von jedem Teilchen mit allen anderen Teilchen berechnet werden. Die Zahl der Interaktionen und damit auch der Rechenaufwand steigt dann quadratisch mit der Zahl der Teilchen. Bei hohen Teilchenzahlen steigt damit die Rechenzeit inakzeptabel an. Eine Möglichkeit, dies zu vermeiden, besteht darin, mehrere Teilchen, die weit entfernt vom aktuellen Teilchen liegen, zu einem Pseudoteilchen zusammenzufassen und nur eine Interaktion zwischen dem aktuellen Teilchen und dem Pseudoteilchen zu berechnen. Als Beispiel kann die Interaktion zwischen einem Stern und einer weit entfernten Galaxie dienen: Der Fehler, der entsteht, wenn alle Sterne der entfernten [[Galaxie]] zu einem einzigen Massepunkt zusammengefasst werden, ist bei normalen Anforderungen vernachlässigbar. Um zu entscheiden, welche Teilchen zu Pseudoteilchen zusammengefasst werden können, werden sogenannte Baumverfahren angewendet. Dabei werden die Teilchen in einem hierarchischen Baum, im zweidimensionalen Fall einem Quadtree, im dreidimensionalen Fall einem [[Octree]] angeordnet.
Bei Molekulardynamik-Simulationen wird dagegen der Raum, in dem die Simulation stattfinden soll, in Simulationszellen eingeteilt. Sowohl die Kräfte als auch die Teilchen werden, wenn sie über den Rand der Zelle hinausgehen, einfach auf der anderen Seite der Zelle wieder eingefügt ([[periodische Randbedingung]]). Um zu verhindern, dass ein Teilchen nun sowohl von der eigentlichen Kraft als auch von deren Spiegelbild auf der anderen Seite erfasst wird, wird diese Kraft ab der sogenannten Cutoff-Distanz (normalerweise die halbe Länge der Zelle) nicht mehr berücksichtigt. Um nun die Anzahl der beteiligten Teilchen zu erhöhen, wird einfach die Simulationszelle beliebig vervielfacht.

== Algorithmen ==
=== Integrationsalgorithmen ===
Integrationsalgorithmen sind der [[Verlet-Algorithmus]], der [[Velocity-Verlet-Algorithmus]], das [[Leapfrog-Verfahren]] und die [[Mehrschrittverfahren#Prädiktor-Korrektor-Methode|Prädiktor-Korrektor-Methode]].

=== Langreichweitige Kräfte ===
Hierzu zählen der [[Barnes-Hut-Algorithmus]], die [[Fast-Multipole-Methode]], die [[Ewaldsumme]] und die [[Particle-Mesh-Ewald-Methode]].

== Software ==

* Aspherix (DCS Computing, Weiterentwicklung von LIGGGHTS)
* Chute Maven (Hustrulid Technologies Inc.)
* Altair EDEM<ref>{{Internetquelle |url=https://www.altair.de/edem/ |titel=Diskrete-Elemente-Modellierung - DEM Software |hrsg=Altair Engineering Inc. |abruf=2022-11-02 |sprache=de}}</ref> ([[Altair Engineering|Altair Engineering Inc.]])
* ELFEN
* ESyS-Particle (Open-Source)
* eXtended Particle System (XPS)
* GROMACS
* GROMOS 96
* LAMMPS (Open-Source)
* LIGGGHTS<ref>{{cite journal |author=Christoph Kloss, Christoph Goniva, Alice Hager, Stefan Amberger, Stefan Pirker |title=Models, algorithms and validation for opensource DEM and CFD-DEM |journal=Progress in Computational Fluid Dynamics |volume=12 |issue=2–3 |pages=140–152 |date=2012-06 |url=https://www.inderscienceonline.com/doi/abs/10.1504/PCFD.2012.047457 |access-date=2021-12-15 |language=en}}</ref> (Open-Source)
* MUSEN<ref>{{Literatur |Titel=MUSEN: An open-source framework for GPU-accelerated DEM simulations |Sammelwerk=SoftwareX |Band=12 |Datum=2020-07-01 |Seiten=100618 |ISSN=2352-7110 |DOI=10.1016/j.softx.2020.100618 |Sprache=en}}</ref> (Open-Source)
* PASIMODO Partikelsimulationspaket
* PFC2D und 3D von ITASCA
* ROCKY DEM
* Simcenter STAR-CCM+
* SimPARTIX vom Fraunhofer IWM
* ThreeParticle (BECKER 3D GmbH)
* TinyDEM<ref>{{cite journal |first1=Roman |last1=Vetter |title=TinyDEM: Minimal open granular DEM code with sliding, rolling and twisting friction |journal=Computer Physics Communications |volume=320 |pages=109942 |doi=10.1016/j.cpc.2025.109942 |date=2025-11-19 |language=en}}</ref> (Open-Source)
* YADE (Open-Source)

== Literatur ==
* M. P. Allen, D. J. Tildesly: ''Computer Simulation of Liquids''. Oxford University Press, 1989, ISBN 0-19-855645-4.
* {{Literatur |Autor=Michael Griebel, Stephan Knapek, Gerhard Zumbusch, Attila Caglar |Titel=Numerische Simulation in der Moleküldynamik: Numerik, Algorithmen, Parallelisierung, Anwendungen |Verlag=Springer-Verlag |Ort=Berlin/Heidelberg |Datum=2004 |Reihe=Springer-Lehrbuch |ISBN=3-540-41856-3 |Sprache=de}}
* [[Nenad Bićanić]]: ''Discrete Element Methods.'' In: Erwin Stein, René de Borst, J. R. Hughes: ''Encyclopedia of Computational Mechanics, Vol. 1''. Wiley, Chichester/West Sussex 2004. ISBN 0-470-84699-2.

== Weblinks ==
* [http://www.ged.rwth-aachen.de/index.php?cat=Tools_n_Methods&subcat=Geomechanical_modeling&page=Discrete_Element_Modeling RWTH Aachen]
* [http://www.isd.uni-stuttgart.de/forschung/fsp/sum/diskrete_elemente_sum/index.html ISD Universität Stuttgart]
* [http://www.rcpe.at/en/xps-at/ Research Center for Pharmaceutical Engineering]

== Einzelnachweise ==
<references />

{{Normdaten|TYP=s|GND=4577038-4|LCCN=sh2013003076}}

[[Kategorie:Computerchemie]]
[[Kategorie:Computerphysik]]

Diskrete-Elemente-Methode - Versionsgeschichte

imported>Leyo: redundant zu DOI