Atom-dynamische Betrachtung zum funktionellen und strukturellen Verhalten von Formgedächtnislegierungen
Zusammenfassung der Projektergebnisse
Metall-Legierungen mit der Fähigkeit, eine mechanisch aufgeprägte Verformung unter Temperaturerhöhung in den Ursprungszustand zurückzustellen, nennt man Formgedächtnislegierungen (FGL). Der Formgedächtnis-Effekt wird durch fest-fest Phasenübergänge der metallischen Kristallgitter verursacht. Diese Phasenübergänge treten zwischen den zwei generischen Gitterstrukturen Austenit und Martensit auf, welche unterscheidbare Feststoff-Phasen darstellen, deren Stabilität sowohl von Belastungszuständen, als auch von der Temperatur abhängt. Der Übergang zwischen diesen Kristallstrukturen werden allgemein martensitische Transformationen genannt. Es handelt sich um diffusionslose Umstrukturierungen, die die atomare Nahordnung erhalten. Austenit ist bei hoher Temperatur und kleiner Last stabil, während Martensit bei niedriger Temperatur und/oder bei großer Last stabil ist. Martensit liegt in verschiedenen Varianten vor, welche geometrisch kongruent und energetisch äquivalent sind. Der Austenit-Martensit Phasenübergang wird durch Temperaturwechsel (Temperatur-induzierter Übergang) oder durch mechanische Belastung (last-induzierter Übergang) ausgelöst. Im Laborexperiment zeigen die Prozessdiagramme charakteristische Hysteresen, die auf Prozess-Irreversibilitäten hinweisen. Mikroskopische Befunde zeigen, dass unter martensitischen Transformationen Makrostrukturen erzeugt werden. Charakteristisch ist die Ausbildung von Domänen, welche aus regelmäßig gefalteten („verzwillingten") Martensitplatten bestehen. Das Ziel des Projekts bestand in einer Modellstudie dieser Transformationen auf der atomistischen Längen- und Zeitskala. Die verwendete Methode der atom- oder molekulardynamischen Simulation ist dabei in der Lage, den Mechanismus der Phasentransformationen zwischen Austenit und Martensit dynamisch nachzustellen und dabei die treibenden physikalischen Kräfte zu analysieren. Es handelt sich dabei um ein numerisches Verfahren zur Berechnung der klassischen Bewegungsgleichungen der Atome in endlichen Kristallen. Dabei kommt der Modellierung der atomaren Wechselwirkungskräfte die wichtigste Bedeutung zu. Trotz wesentlicher Fortschritte auf dem Gebiet der Potentialmodellierung in den vergangenen Jahrzehnten gibt es bis heute keine Modellpotentiale, welche alle Eigenschaften von FGL qualitativ und quantitativ zufriedenstellend abdecken. Stets müssen Kompromisse bei der Anpassen der Potentiale eingegangenen werden. Ein weiteres Problem realistische, drei-dimensionale Modelle besteht in ihrer atomaren Größe. Diese Arbeit vermeidet Probleme dieser Art, in dem sie a priori ein künstliches Modellsystem betrachtet. Es hat sich herausgestellt, dass das verwendete Modellsystem trotz seiner konzeptionellen Vereinfachungen die wesentlichen Merkmale der martensitischen Transformationen qualitativ abbilden karm. Das Projekt verfolgte damit das Ziel einer Prinzipstudie, in der die den martensitischen Transformationen auf der atomaren Skala zugrundeliegenden physikalische Mechanismen untersucht wurden. Das verwendete Modell basiert auf Lennard-Jones (L-J) Wechselwirkungen in zwei Raumdiraensionen. Dazu werden die Paar-Potentiale in binären Legierungen so an gepasst, dass Transformationen in Simulations-Experimenten beobachtet werden können. Unter Berücksichtigung von zwei Atomsorten A und B können stabile Kristallgitter konstruiert werden. Die Wechselwirkungsparameter können so identifiziert werden, dass temperaturabhängige Gittertransformationen in MD Simulationen abgebildet werden können. Austenit, im 2D Modell durch quadratische, ineinander gestellte Komponentengitter repräsentiert, ist stabil bei hoher Temperatur. Durch eine shear/shuffle Transformation kann der Modell-Austenit in Varianten von Modell-Martensit transformieren, welche bei kleiner Temperatur stabil sind. Der Modell-Martensit hat eine hexagonale Struktur. Ähnlich dem 3D Fall kommt es auch in 2D zur Ausprägung von Martensitvarianten, die durch die Scherrichtung und die Verschiebungsrichtung der Komponentengitter relativ zueinander definiert sind. Im 2D Modell können vier prinzipielle Varianten von Martensit unterschieden werden. Die Gitterstabilität hängt von der Temperatur und von Lastzuständen ab. Basierend auf diesem Modell wurden Simulationsstudien in Testkörpern von einigen hunderttausend Atomen unter den folgenden Fragestellungen untersucht: 1. Untersuchung des Modellmaterials in isothermen Zugversuche (last- und deformationskontrolliert); 2. Untersuchung der martensitischen Umwandlung in den Modellkörpern im Hinblick auf die Entstehung martensitische Mikrostruktur; 3. Einstellung des thermodynamische Gleichgewichts; 4. Einfluss der Mikrostruktur auf das Transformationsverhalten hinsichtlich der Hysteresenbildung. Die Simulationsexperimente zeigen die dynamischen Abläufe der Ausbildung einer charakteristischen Mikrostruktur detaiUiert. Wir beobachten die Nukleation und das Wachs tum keilförmiger Martensit domänen, welche eine charakteristische Zwillingsstruktur aufweisen. Diese Wachstumsprozesse führen auf der mesoskopischen Längenskala zur Ausbildung ausgeprägter Domänstrukturen die durch signifikante Domängrenzen separiert werden. Hier werden Gitterdefekte generiert, welche teilweise durch Versetzungsbewegungen mobilisiert werden können oder immobil im Gitter verbleiben. Solche Defekte können selbst unter temperatur- oder lastinduzierter Rücktransformation erhalten bleiben. Es zeigt sich, dass solche Defektresiduen als Nukleationszentren in nachfolgenden Transformationen dienen und akkumulieren können. In zyklischen Transformationsprozessen wird auf diese Weise eine permanente Defekt-Mikrostruktur erzeugt, welche ein „Morphologie- Gedächtnis"auf der mikrostrukturellen Längenskala repräsentiert. Da Gitterdefekte stets eine potentielle Exzessenergie tragen, bewirkt eine solche Defekt akkumulation andererseits auch eine irreversible Veränderung der Energie eines gegebenen Kristalls. Auf diese Weise wird auch das thermodynamische Phasengleichgewicht verändert und mit der Transformationsgeschichte eines gegebenen Testkörpers in Beziehung gesetzt. Die Methode der molekulardynamischen Simulation erfordert beträchtliche Rechnerkapazitäten. Die hier dargestellten Simulationsexperimente konnten nur auf spezialisierten Parallelrechnern abgebildet werden. Dank der konzeptionellen Modellannahmen konnten die Rechenzeitressourcen besonders effektiv eingesetzt werden. Hier Spielte die Beschränkung auf zwei Raumdimensionen die wichtigste Rolle. Sie ermöglichte eine Betrachtung sowohl zugunsten der simulierten Zeitskalen, deren Behandlung für dreidimensionale Modelle außerhalb der verfügbaren Mittel läge, als auch zugunsten der Bandbreite der simulierten thermo-mechanischen Prozesse.
Projektbezogene Publikationen (Auswahl)
- 2005 Workshop on Macroscopic Limits and Modulation Theory for PDEs and discrete lattice models, Weierstraß Institut für angewandte Mathematik und Stochastik Berlin, 3. 11. 2005
- 2006 Third International Conference Multiscale Materials Modeling (MMM2006), Freiburg, 18. - 22. 9. 2006
- Fachausschusses Computersimulation der DGM, Max-Planck-Istitut für Eisenforschung, 15. 11. 2006
- 6. Materialwissenschaftlicher Tag and der Ruhr-Universität Bochum, 18. 10. 2007
- 6th International Congress on Industrial and Applied Mathematics (ICIAM 07), Zürich/Schweiz, 16. - 20. Juli 2007. (Eingeladener Vortrag)
- Atomistic Simulations of martensitic transformations in Lennrad-Jones solids. Part I. In Science and Supercomputing in Europe, pages 950–959. HPC-Europa, 2007
O. Kastner and G.J. Ackland
- Atomistic Simulations of martensitic transformations in Lennrad-Jones solids. Part II. In Science and Supercomputing in Europe, pages 960–962. HPC-Europa, 2007
O. Kastner and G.J. Ackland
- Gastvortrag an der School of Physics, Universität Edinburgh, Januar 2007
- 2008 International Conference on Martensitic Transformations (ICOMAT), Santa Fe, NM, USA, 29.6. - 5.7. 2008
- Conference on Smart Materials, Adaptive Structures and Intelligent Systems (SMASIS 2008), Ellicott City, MD, USA, 28.-30. 10. 2008
- Gastvortrag am Department of Mechanical and Aerospace Engineering, North Carolina State University, USA,November 2008
- Gastvortrag im Mathematischen Institut der Universität Oxford, 3. 10. 2008
- Load-induced martensitic transformations in pseudo-elastic Lennrad-Jones cyrstals. In Proceedings of the SMASIS 2008, Ellicott City, MD, USA, October 28-30 2008. ASME
O. Kastner and G.J. Ackland
- Martensitic transformations in 2D Lennard-Jones crystals. In Proceedings of ICOMAT 2008, International conference on martensitic transformations, Santa Fe, NM, USA, June 29 – July 5 2008
O. Kastner and G.J. Ackland
- 2. Forum Metallplastizität — Thermo-mechanisch gekoppelte Phänomene, Universität Karlsruhe, 26. 5. 2009
- Conference on Smart Materials, Adaptive Structures and Intelligent Systems (SMASIS 2009), Oxnard, CA, USA, 21.-23. 9. 2009
- European Symposium of Martesitic Transformations (ESOMAT 2009), Prag/Tschechien, 7.-11. September 2009
- Gastvortrag an der Ben-Gurion Universität, Beer-Sheva/Israel, 8. 3. 2009
- Gastvortrag an der Fudan-Universität, Shanghai/China, 3. 7. 2009
- Kolloquium im Institut für Nanotechnologie des Technion, Haifa/Israel, 11. 3. 2009
- OxMOS workshop on atomistic models of solids, Universität Oxford, 7.-8. Dezember 2009
G. Ackland
- Second International Conference on Smart Materials and Nanotechnology in Engineering (SMN2009), Weihai/China, 8.-11. Juli 2009
- Seminarvortrag am Weierstraß Institut für angewandte Mathematik und Stochastik, Berlin, 24. 11. 2009
- Seminarvortrag im Interdisciplinary Centre for Advanced Materials Simulation (ICAMS), Ruhr-Universität Bochum, April 2009
- Institutsvortrag am Weierstraß Institut für angewandte Mathematik und Stochastik, Berlin, 26. 4. 2010. (Eingeladener Vortrag)
- Spring Meeting of the European Materials Research Society (E-MRS 2010), Strasburg/Frankreich, 7.-11. Juni 2010. (Eingeladener Vortrag)
- VDI Arbeitskreis Werkstofftechnik, Fachhochschule Würzburg-Schweinfurt, Fakultät Maschinenbau, 30. 3. 2010
- Vortrag auf der Fachtagung “phase transformations”, Mathematisches Forschungszentrum Oberwolfach, 30. 5. - 4. 6. 2010
