📑 Inhaltsverzeichnis
Diese Arbeit stellt eine Finite-Elemente-Formulierung für gekoppelten chemo-mechano-thermodynamischen Kontakt bei großen Deformationen vor. Basierend auf der Kontakttheorie von Sauer et al. (2022) werden sechs Felder betrachtet: Deformation und Temperatur $T$ der beiden Kontaktkörper sowie ein Grenzflächenbindungsfeld $\phi$ und eine Grenzflächentemperatur. Letztere wird durch chemische und mechanische Energiedissipation an der Grenzfläche bestimmt. Der Fokus liegt auf der Entwicklung von Bindung und Entbindung (Bonding/Debonding) und deren Kopplung an den mechanischen und thermischen Kontaktzustand.
Warum sechs Felder? Kontakt ist ein Mehrfeldproblem
Wenn zwei Körper einander berühren, treffen dort drei Physiken aufeinander: Die Mechanik entscheidet über Verformung, Kontaktdruck $p$ und Spaltweite $g$, die Thermodynamik über Wärmeleitung und Temperatur, und die Chemie darüber, ob an der Grenzfläche neue Bindungen entstehen oder alte brechen. Jeder der beiden Kontaktkörper trägt zwei Felder – seine Deformation und seine Temperatur $T$ –, das ergibt vier. Die Grenzfläche selbst fügt zwei weitere hinzu: das Grenzflächenbindungsfeld $\phi$ und die Grenzflächentemperatur. Erst zusammen ergeben diese sechs Felder ein geschlossenes Bild davon, warum eine Verbindung haftet, sich erwärmt oder sich wieder löst.
Das Bindungsfeld $\phi$ hat dabei eine vertraute chemische Gestalt: Es ist ein Reaktionsfortschritt an der Oberfläche – der Anteil der Grenzfläche, der bereits gebunden ist, mit Werten zwischen $\phi = 0$ (völlig ungebunden) und $\phi = 1$ (vollständig gebunden). Wie ein Umsatzgrad wächst oder schwindet $\phi$ im Zeitverlauf und wird nicht als Konstante, sondern über eine Evolutionsgleichung beschrieben.
Das quadratische Kontaktpotential als Baukasten
Ausgangspunkt der Formulierung ist ein quadratisches Kontaktpotential, aus dem sich die Kontakttraktionen als Funktionen von Druck, Spalt und Bindungszustand ergeben. Darauf aufbauend schlagen die Autoren mehrere elementare Modelle vor, etwa für druck- und spaltabhängige Bindung: Die Rate, mit der neue Bindungen entstehen, hängt davon ab, wie fest die Flächen aufeinandergepresst werden ($p$) und wie nah sie einander sind ($g$). In elementarer, ratenbasierter Form lässt sich das wie eine kinetische Gleichung lesen,
$$\dot{\phi}_{\text{Bond}} = k(T)\,f(p,g)\,(1-\phi), \qquad k(T) = k_0\,\exp\!\left(-\frac{E_a}{R\,T}\right),$$mit einer Arrhenius-artigen Temperaturabhängigkeit, wie sie aus der Reaktionskinetik bekannt ist: Aktivierungsenergie $E_a$, Gaskonstante $R$, Vorfaktor $k_0$. Der Faktor $(1-\phi)$ sorgt dafür, dass die Bindungsbildung zum Sättigungszustand hin ausläuft – das freie (noch ungebundene) Grenzflächenareal wird knapp, genau wie ein Edukt in einer Reaktion erster Ordnung. Gleichzeitig kann Entbindung als gegenläufiger Prozess modelliert werden, sodass Binden und Lösen auch gleichzeitig auftreten können.
Chemie an der Grenzfläche: Wärme aus Reibung und Reaktion
Besonders lehrreich ist die Behandlung der Grenzflächentemperatur. Sie wird durch zwei Dissipationsquellen gespeist: mechanische Energie (etwa aus Reibungs- und Separationsprozessen beim Kontakt) und chemische Energie aus exothermen Bindungsreaktionen. In einer Energiebilanz an der Grenzfläche $\Gamma$ steht das nebeneinander,
$$\rho\,c\,\dot{T}_\Gamma \;=\; \underbrace{q_{\text{chem}}(\dot{\phi})}_{\text{Reaktionswärme}} \;+\; \underbrace{\mathcal{D}_{\text{mech}}}_{\text{mechanische Dissipation}} \;-\; \underbrace{q_{\text{cond}}}_{\text{Wärmeleitung in die Körper}},$$und die Wärmeleitung $q_{\text{cond}}$ koppelt die Grenzflächentemperatur zurück an die Volumenkörper. Hier entsteht eine echte Rückkopplungsschleife, die Chemikern vertraut ist: Eine exotherme Bindungsreaktion ($\Delta H < 0$) erwärmt die Grenzfläche, die höhere Temperatur erhöht über $k(T)$ die Reaktionsrate, diese setzt erneut Wärme frei – ein selbstverstärkender Prozess, wie man ihn vom Aushärten zweikomponentiger Klebstoffe kennt. Umgekehrt zeigt das Beispiel der thermischen Härtung, dass Wärmezufuhr die Bindung erst ermöglicht, während thermische Ausdehnung der Körper den Spalt $g$ verändert und damit die Mechanik zurückspielt.
Die Thermodynamik liefert dabei die Ordnungsmacht: Die Modelle werden so formuliert, dass die gesamte Energiedissipation – chemisch wie mechanisch – der Clausius-Duhem-Bedingung genügt, also nie negativ wird. Ein Bonding-Modell, das Energie aus dem Nichts erzeugen würde, ist damit ausgeschlossen; Bindungsenergie, die bei der Entbindung „freigesetzt" wird, muss messbar dissipieren und erwärmen.
Numerik: monolithisch, linearisiert, kondensiert
Die Implementierung ist monolithisch: Alle sechs Felder werden simultan in einem gekoppelten nichtlinearen Gleichungssystem gelöst, statt die Physiken nacheinander („gestaffelt") zu behandeln. Für das Newton-Raphson-Verfahren wird dafür die vollständige konsistente Linearisierung hergeleitet – die Ableitungen der Gleichungen nach allen Feldern, einschließlich der Kopplungsterme. Das ist aufwendig, zahlt sich aber in robuster und schneller Konvergenz aus. Die Zeitintegration erfolgt implizit, was für steife, ratenbasierte Evolutionsgleichungen wie $\dot{\phi}$ essenziell ist.
Als Ansatzräume dienen klassische Formfunktionen ebenso wie isogeometrische (NURBS-basierte). Letztere beschreiben die Geometrie der Körper direkt mit glatten Funktionen; die C¹-stetigen Oberflächen machen die Kontaktbeschreibung besonders maßhaltig. Ein weiterer Kniff: Sind die Bindungsstellen Materialpunkte, lässt sich die Bindungsvariable lokal kondensieren – sie wird auf der Ebene einzelner Integrationspunkte aus dem globalen System eliminiert (mathematisch ein Schur-Komplement-artiger Schritt), ohne die Lösung zu verändern. Das globale System schrumpft, die Simulation wird schneller.
Die Formulierung ist vielseitig einsetzbar und wird durch Beispiele illustriert: druck- und spaltabhängige Bindung, exotherme Bindungsreaktionen, thermische Härtung und Ausdehnung sowie gleichzeitiges Binden und Lösen. Hintergrund: Die Simulation von Haftungsprozessen ist entscheidend für die Entwicklung verbundener Materialien in der Werkstofftechnik – von Klebverbindungen über Verbundwerkstoffe bis hin zu Fügeprozessen, bei denen Chemie, Druck und Wärme untrennbar zusammenwirken.
Lernimpuls
Kernidee: An einer Kontaktfläche bilden Chemie, Mechanik und Wärmeleitung ein einziges gekoppeltes System. Die Bindungsvariable $\phi$ spielt dabei die Doppelrolle eines Reaktionsfortschritts und einer mechanischen Zustandsgröße: Sie wächst nach ratenbasierten Evolutionsgleichungen (mit Arrhenius-artiger Temperaturabhängigkeit), sie erzeugt Zugkräfte, die den mechanischen Kontakt verändern, und sie setzt bei exothermen Bindungsreaktionen Wärme frei, die über die Grenzflächentemperatur auf die eigene Reaktionsrate zurückwirkt. Wer diese Rückkopplungen versteht, kann Haften und Lösen nicht nur beobachten, sondern berechnen – und gezielt Werkstoffe und Prozesse so gestalten, dass Verbindungen halten, wo sie sollen, und sich lösen, wo sie sollen.
Häufige Fehlvorstellung: „Haftung ist eine Stoffkonstante – ein Material haftet oder haftet nicht." Richtig ist: Haftung ist ein Zustand in Entwicklung. $\phi$ wächst oder schwindet je nach Druck, Spaltweite und Temperatur; dieselbe Fläche kann gleichzeitig binden und sich lösen. Die zweite Fehlvorstellung: „An einer Kontaktfläche entsteht Wärme höchstens durch Reibung." Tatsächlich heizt auch die Chemie: Exotherme Bindungsreaktionen tragen zur Grenzflächentemperatur bei – und beschleunigen damit ihre eigene Kinetik. Wer nur eine der drei Physiken betrachtet, versteht Haftprozesse nur zur Hälfte.
Prüfe dich:
- Erkläre, warum das Grenzflächenbindungsfeld $\phi$ über eine Evolutionsgleichung $\dot{\phi}$ modelliert wird – welche Analogie besteht zum Umsatzgrad einer chemischen Reaktion, und warum läuft die Bindungsbildung beim Sättigungswert $\phi = 1$ aus?
- Eine exotherme Bindungsreaktion setzt Wärme frei ($\Delta H < 0$). Skizziere die Rückkopplungsschleife zwischen Reaktionswärme, Grenzflächentemperatur $T_\Gamma$ und Arrhenius-Rate $k(T) = k_0\,e^{-E_a/(RT)}$ – und erkläre, warum dieser Prozess sich selbst verstärken kann.
- Begründe mit der Clausius-Duhem-Bedingung, warum ein thermodynamisch konsistentes Bonding-Modell dissipativ sein muss. Wohin fließt die Bindungsenergie, wenn eine Grenzfläche entbunden wird?
- Vergleiche klassische Lagrange-Formfunktionen mit isogeometrischen NURBS-Ansätzen: Warum sind glatte, C¹-stetige Oberflächen für die Kontaktmechanik vorteilhaft?
- Was bedeutet es, dass die Bindungsvariable „lokal kondensiert" werden kann, wenn Bindungsstellen Materialpunkte sind? Erkläre, warum das das globale Gleichungssystem verkleinert, ohne die Lösung zu verändern.
Bei Schwierigkeiten: Stell dir Klettverschluss vor: $\phi$ ist der Anteil der geschlossenen Häkchen. Drücken (Druck $p$) und enger Kontakt (kleiner Spalt $g$) schließen mehr Häkchen – aber nur, solange noch freie Häkchen da sind ($1-\phi$). Für die Rückkopplung hilft das Alltagsbild vom zweikomponentigen Kleber: Beim Aushärten wird der Klebstoff warm, die Wärme beschleunigt das Aushärten – Chemie heizt sich selbst. Und für die Thermodynamik: Energie verschwindet nie – Bindungsenergie, die beim Lösen „frei" wird, muss als Wärme wieder auftauchen und die Grenzfläche erwärmen. Genau das zwingt die Clausius-Duhem-Bedingung dem Modell auf.
Entitäten: finite-elemente-methode | methode chemo-mechano-thermodynamik | konzept kontaktmechanik | konzept grenzflächenbindungsfeld | konzept bonding-debonding | konzept arrhenius-gleichung | konzept exotherme reaktion | reaktion thermische härtung | konzept newton-raphson-verfahren | methode isogeometrische analyse | methode implizite zeitintegration | methode lokale kondensation | methode
📄 Quelle
🧪 Verwandte Rechner
Mit diesen interaktiven Werkzeugen können Sie das Thema vertiefen:
🔬 bindungspotential → 🔬 druck-flaechen-rechner → 🔬 waermeleitung →