Deformation gradient averaging regularization for third medium contact

Titel: Regularisatie door gemiddelde vervormingsgradiënt voor contact van het derde medium

1. Probleemstelling

In de moderne materiaalwetenschap en structuuroptimalisatie (zoals topologie-optimalisatie) is het simuleren van contactmechanica onder grote vervormingen essentieel. Traditionele contactmethoden (zoals straf- of Lagrange-multiplicator-methoden) vereisen complexe zoekalgoritmen om contactoppervlakken te detecteren, wat problematisch is bij veranderende topologieën.

De derde-methode (third medium method) is een alternatief waarbij de holle ruimte tussen lichamen wordt gevuld met een fictief, zeer flexibel materiaal. Dit materiaal is zacht voordat contact optreedt, maar wordt stijf bij extreme compressie (bijna nul volume), waardoor contact wordt afgedwongen zonder een zoekalgoritme.

Echter, bestaande modellen voor het derde medium lijden onder instabiliteit in de pre-contactfase (overmatige vervorming). Om dit op te lossen, gebruiken eerdere modellen regularisatie via de straling van tweede afgeleiden van de verplaatsing (second gradients). Dit heeft twee grote nadelen:

1. Het vereist kwadratische elementen (tweede orde) of de toevoeging van extra vrijheidsgraden aan lineaire elementen, wat de implementatie complex maakt.
2. Het straalt vaak meer vervormingsmodes uit dan strikt noodzakelijk is.

2. Methodologie

De auteurs stellen een nieuwe regularisatiemethode voor die volledig afziet van het berekenen van tweede afgeleiden van de verplaatsing. In plaats daarvan gebruiken ze een elementgewijze gemiddelde vervormingsgradiënt.

De totale energie van het derde medium bestaat uit drie componenten:

1. Contactterm: Gebaseerd op het volumetrische deel van de neo-Hookean wet. Deze term zorgt ervoor dat het materiaal oneindig stijf wordt bij extreme compressie ($J \to 0$), waardoor penetratie wordt voorkomen. Het isochorische deel wordt weggelaten om stabiliteit te vergroten.
2. Lineair elastische term (constante stijfheid): Een term gebaseerd op de kleine-rek theorie ($W_{ss} = \frac{1}{2}(F-I):D_{3M}:(F-I)$). Deze zorgt voor een constante stijfheid in de pre-contactfase en voorkomt overmatige vervorming van de elementen voordat contact optreedt.
3. Regularisatieterm (Deformation Averaging): Dit is de kerninnovatie. In plaats van de gradiënt van de vervormingsgradiënt ($\nabla F$) te straffen, wordt de afwijking van de vervormingsgradiënt ($F$) op de integratiepunten ten opzichte van de gemiddelde vervormingsgradiënt over het element (berekend in het elementcentrum, $\bar{F}$) gestraft.
   • Formule: $W_{\bar{F}} = \frac{1}{2}\kappa_{\bar{F}} \| F - \bar{F} \|^2$.
   • Hierbij is $\kappa_{\bar{F}}$ een straffingsparameter.

Voordeel: Deze aanpak maakt het mogelijk om standaard lineaire (eerste-orde) elementen te gebruiken zonder extra vrijheidsgraden. De regularisatie werkt door de uniformiteit van de vervormingsgradiënt over de Gauss-integratiepunten van een element te controleren.

3. Belangrijkste Bijdragen

• Implementatievriendelijkheid: De methode vereist geen kwadratische elementen en geen extra vrijheidsgraden, wat de implementatie in bestaande FEM-codes (zoals OOFEM) aanzienlijk vereenvoudigt.
• Robuustheid: Het combineert een constante stijfheid in de pre-contactfase met een asymptotische stijfheid bij contact, wat de numerieke stabiliteit verbetert.
• Universeel toepasbaar: De methode is geldig voor zowel 2D als 3D problemen en werkt goed in combinatie met topologie-optimalisatie.
• Energiepotentiaal: Alle termen in het model hebben een energiepotentiaal, wat de symmetrie van de globale stijfheidsmatrix behoudt en belangrijk is voor stabiliteitsanalyses.

4. Resultaten en Validatie

De methode is getest op vier bekende benchmarks:

1. C-vorm benchmark:

   • Vergelijking met traditionele node-to-segment contactmethoden.
   • De resultaten tonen aan dat de nieuwe methode robuust is voor verschillende mesh-dichtheden. De contacten worden correct afgedwongen, hoewel de scherpte van het contactpunt iets minder is dan bij traditionele methoden (een kenmerk van de derde-methode).
   • De parameter $\kappa_{\bar{F}}$ moet schalen met de meshgrootte voor optimale resultaten.

2. Gesloten doos benchmark (Zelfcontact):

   • Simulatie van een doos die onder druk buigt en zelfcontact maakt (binnen- en buitenkant).
   • De kracht-verplaatsingsdiagrammen komen zeer goed overeen met traditionele contactmethoden.
   • De methode slaagt erin complexe zelfcontactsituaties te simuleren zonder dat het contact zoekalgoritme faalt.

3. Beperkte knik en zelfcontact:

   • Een instabiliteitsprobleem waarbij een balk wordt gecomprimeerd, knikt, en vervolgens contact maakt met wanden en zichzelf.
   • De methode kan de verschillende fasen van het knikgedrag en de overgang naar contact correct volgen, inclusief de overdracht van krachten via het derde medium.

4. Configuratieve krachten:

   • Validatie tegen een analytische oplossing voor de configuratieve kracht (J-integraal) bij een elastisch blok in contact met een ronde indenter.
   • De resultaten tonen een uitstekende overeenkomst met de analytische oplossing, zelfs bij grote vervormingen en schuifcontact, wat de nauwkeurigheid van de methode bevestigt.

5. Betekenis en Conclusie

Dit artikel introduceert een elegante en efficiënte oplossing voor een langdurig probleem in de contactmechanica: het stabiliseren van het derde-middel contact zonder de complexiteit van hogere-orde elementen of extra vrijheidsgraden.

De deformation gradient averaging regularisatie maakt het mogelijk om simpele, lineaire elementen te gebruiken voor complexe contactproblemen met grote vervormingen. Dit opent de deur voor bredere toepassing in:

• Topologie-optimalisatie: Waar contact vaak voorkomt en zoekalgoritmen problematisch zijn.
• Multiphysica simulaties: Waar lucht-achtige ruimtes moeten worden gemeshed.
• Stabiliteitsanalyses: Vanwege het behoud van de symmetrische stijfheidsmatrix.

De auteurs concluderen dat hun model robuust, universeel en eenvoudig te implementeren is, waardoor het een waardevolle toevoeging vormt aan de toolbox van computergestuurde mechanica.