A robust mixed finite element formulation for third medium contact

Technische Samenvatting: Een Robuuste Mixed Finite Element Formulering voor Derde Medium Contact

Probleemstelling
Klassieke eindige elementen contactformuleringen vertrouwen op expliciete discretisatie van contactinterfaces, waarbij beperkingen worden afgedwongen via penalty-methoden, Lagrange-multiplicatoren of barrière-methoden. Hoewel deze volwassen zijn, kampen deze benaderingen met aanzienlijke algoritmische lasten in scenario's met grote glijbewegingen, zelfcontact, ernstige mesh-distorsie of topologieoptimalisatie waarbij de contactgrenzen vooraf niet bekend zijn. De "third medium contact" (TMC) methode adresseert dit door expliciete beperkingen te vervangen door een zeer compliant fictief medium dat de ruimte tussen potentieel contactmakende lichamen vult. Dit medium blijft zacht vóór contact, maar verstijft drastisch onder compressie, waardoor contactkrachten via het continuüm worden overgebracht.

Een centrale uitdaging in TMC is echter stabilisatie. Het medium moet zacht genoeg zijn om spurieuze krachten vóór contact te vermijden, maar stabiel genoeg om extreme compressie en afschuiving zonder mesh-distorsie te weerstaan. Bestaande stabilisatiestrategieën vereisen vaak hogere-orde eindige elementen (om tweede afgeleiden van verplaatsing te berekenen), introduceren globale hulpvelden die de systeemomvang vergroten, of vertrouwen op gespecialiseerde discretisaties. Er is behoefte aan een robuust stabilisatiemechanisme dat compatibel is met eenvoudige, computationeel efficiënte eerste-orde eindige elementen en dat de directe evaluatie van tweede verplaatsingsafgeleiden vermijdt.

Methodologie
De auteurs stellen een mixed eindige elementen formulering voor die een auxiliary-field stabilisatie introduceert voor het derde medium. De kern van de methode omvat:

1. Introductie van een Hulpveld: Een deformatiegradiënt-achtig veld, aangeduid als $\Theta$, wordt geïntroduceerd in het derde medium en behandeld als een onafhankelijke onbekende, ontkoppeld van het fysieke verplaatsingsveld $u$.
2. Zwakke Koppeling: Het hulpveld $\Theta$ wordt zwak gekoppeld aan de fysieke deformatiegradiënt $F$ via een penalty-term ($W^p_{tm} \propto |\Theta - F|^2$). Dit zorgt ervoor dat $\Theta$ de fysieke deformatie volgt.
3. Gradiënt Regularisatie: In plaats van de gradiënt van de fysieke deformatiegradiënt ($\nabla F$) te bestraffen, wat tweede afgeleiden van verplaatsing zou vereisen, bestraft de formulering de gradiënt van het hulpveld ($\nabla \Theta$). De regularisatie-energie wordt gedefefinieerd als $W^r_{tm} \propto |\nabla \Theta|^2$.
4. Mixed Discretisatie: De formulering staat onafhankelijke interpolatie van $u$ en $\Theta$ toe. De auteurs onderzoeken zowel continue als element-gewijs discontinue benaderingen voor $\Theta$.
   • Continue $\Theta$: Knopenwaarden worden gedeeld tussen elementen, wat zorgt voor inter-element koppeling.
   • Discontinue $\Theta$: Vrijheidsgraden zijn lokaal aan elk element toegewezen. Deze kunnen worden geëlimineerd via statische condensatie, waardoor het globale systeem alleen de verplaatsingsvrijheidsgraden behoudt.

Belangrijkste Bijdragen

• Compatibiliteit met Lage Orde: De formulering bereikt gradiënt-type regularisatie zonder tweede afgeleiden van het verplaatsingsveld te berekenen, wat het compatibel maakt met standaard eerste-orde (lineaire) en tweede-orde eindige elementen.
• Mechanisme van Stabilisatie: Het artikel verduidelijkt dat voor lage-orde verplaatsingselementen (waarbij $F$ element-gewijs constant is), een continue interpolatie van $\Theta$ essentieel is. Het draagt de element-gewijze sprongen in $F$ over naar niet-verdwijnende gradiënten van $\Theta$, waardoor de regularisatie wordt geactiveerd. In tegenstelling hiermee kunnen discontinue $\Theta$-velden geen inter-element gradiënten genereren, wat het stabilisatie-effect voor lage-orde elementen verzwakt.
• Statische Condensatie: Voor discontinue hulpvelden demonstreert de methode dat de extra onbekenden lokaal gecondenseerd kunnen worden, waardoor de globale systeemomvang niet toeneemt terwijl de lokale stabilisatie behouden blijft.

Numerieke Resultaten
De formulering werd beoordeeld via verschillende benchmark-problemen:

• Twee-blok Compressie: Toonde aan dat de schalingparameter $\gamma$ van het derde medium de afweging tussen de residuele kloof (contactnauwkeurigheid) en de nietlineaire convergentierobustheid regelt. Een kleinere $\gamma$ levert nauwer contact op, maar vereist meer Newton-iteraties. De methode bleef stabiel onder gecombineerde compressie en afschuiving.
• Zelfcontact binnen een Box:
  • Interpolatie Vergelijking: Continue hulpvelden ($T_u^1 T_\Theta^1$, $T_u^2 T_\Theta^2$) boden robuuste contactsluiting. Het mixed paar $T_u^2 T_\Theta^1$ (kwadratische verplaatsing, lineair hulpveld) bood resultaten die vergelijkbaar waren met $T_u^2 T_\Theta^2$, maar met minder hulpveld-vrijheidsgraden.
  • Discontinu versus Continu: Discontinue hulpvelden vertoonden een verminderde effectiviteit van de regularisatie, met name voor lineaire elementen waar de gradiëntterm verdween of de elementen niet kon koppelen.
  • Vergelijking met Deformatiegradiënt Middeling: Vergeleken met de methode van Faltus et al. [18], vereiste de voorgestelde mixed formulering aanzienlijk minder belastingsstappen en Newton-iteraties, wat een superieure nietlineaire robuustheid aantoont. Hoewel de mixed formulering een groter algebraïsch systeem heeft, compenseerde de verbeterde convergentie voor de kosten bij gemiddelde mesh-groottes, hoewel de kosten van de hulpveld-vrijheidsgraden dominant werden bij extreem fijne meshes ($>10^6$ onbekenden).
• C-vormige Box: De formulering legde grote rotaties, progressief zelfcontact en ernstige compressie van het derde medium succesvol vast. De kracht-verplaatsingsrespons kwam nauw overeen met klassieke node-to-segment contactformuleringen, met vloeiende overgangen naar contact.
• 3D Uitbreiding: Een driedimensionaal zelfcontactvoorbeeld bevestigde de toepasbaarheid van de methode buiten de plane strain-instellingen.

Betekenis en Claims
Het artikel claimt dat de voorgestelde hulpveld-stabilisatie een robuust, laag-orde-compatibel alternatief biedt voor klassieke contactalgoritmen en bestaande TMC-regularisaties. Door de directe evaluatie van tweede verplaatsingsafgeleiden te vermijden, maakt het het gebruik van eenvoudige, computationeel goedkope eindige elementen (lineair en kwadratisch) mogelijk voor complexe contactproblemen met grote vervormingen en zelfcontact.

De auteurs benadrukken dat de effectiviteit van de methode intrinsiek verbonden is met de continuïteit van het hulpveld; continue interpolaties zijn noodzakelijk om het gradiënt-gebaseerde stabilisatiemechanisme voor lage-orde verplaatsingselementen te genereren. De formulering is bijzonder geschikt voor toepassingen waar contactgrenzen onbekend zijn of evolueren, zoals in contact-ondersteunde mechanismen, soft robotics en topologieoptimalisatie, waarbij de continuüm-gebaseerde aanpak de noodzaak voor expliciete contactzoekalgoritmen en active-set strategieën elimineert. Het werk concludeert dat hoewel de methode een groter lokaal algebraïsch systeem ruilt voor verbeterde globale robuustheid, het een levensvatbaar pad biedt voor het simuleren van ernstige contactscenario's zonder de complexiteit van hogere-orde elementen of gespecialiseerde mesh-technologieën.