Solving the Dissipation Inequality not as a constitutive restriction

Technische Samenvatting: Het oplossen van de dissipatie-ongelijkheid als een randvoorwaardelijke vergelijking

Probleemstelling
Het artikel behandelt een fundamentele uitdaging in de continuümmechanica: de moeilijkheid om de Tweede Wet van de Thermodynamica (de dissipatie-ongelijkheid) te accommoderen als een strikte beperking wanneer het constitutieve gedrag van een materiaal niet volledig bekend is of op een wijze is gespecificeerd die onbedoeld de thermodynamische principes kan schenden. Specifiek onderzoeken de auteurs een scenario waarin een voorgeschreven constitutieve vergelijking voor een snelheid-afhankelijke, elastisch-plastische staaf opzettelijk wordt gekozen om de Dissipatie-ongelijkheid over een specifiek tijdsinterval te schenden. Het doel is om een zelfconsistent procedure te ontwikkelen die een dergelijke "gebrekkige" specificatie als input accepteert, maar een oplossing produceert die strikt voldoet aan de Tweede Wet met minimale afwijking van het voorgeschreven materiaalgedrag.

Methodologie
De auteurs formuleren het probleem als een niet-lineaire geoptimaliseerd constrained optimalisatieprobleem. In plaats van de Dissipatie-ongelijkheid te behandelen als een restrictie op de keuze van constitutieve functies (zoals traditioneel het geval is), behandelen zij het als een randvoorwaardelijke vergelijking binnen het oplossingsproces.

1. Wiskundige Formulering:

   • De evolutie van de plastische rek-snelheid wordt ontleed in een voorgeschreven deel, $f_c(\sigma, t)$, en een correctief veld, $a(x, t)$.
   • Het systeem zoekt naar de velden voor verplaatsing ($u$), plastische rek ($p$), het correctieve veld ($a$), en een dissipatiedichtheidsveld ($s$).
   • Het doel is om de integraal van een functionaal $H(a, s) = \frac{1}{2}c_a a^2 + \frac{1}{2}c_s s^2$ over de tijd te minimaliseren. Deze minimalisatie zorgt ervoor dat het correctieve veld $a$ en de dissipatie $s^2$ zo klein mogelijk zijn, wat een "vertrouwen" in het voorgeschreven constitutieve model weerspiegelt, terwijl de niet-negatieve dissipatie strikt wordt afgedwongen.
   • De randvoorwaarden omvatten krachtevenwicht, de gewijzigde constitutieve relatie en de Dissipatie-ongelijkheid.

2. Analytische Oplossing:

   • Een gesloten vorm oplossing wordt afgeleid voor het 1D-staafprobleem. De optimalisatie levert een specifieke selectie van het correctieve veld $a(t)$ en de dissipatie $s(t)$ die de randvoorwaarden bevredigt.
   • In het limiet waar de ratio van de materiaalconstanten naar nul nadert, selecteert de oplossing natuurlijk een toestand waarin de plastische rek-snelheid nul wordt telkens wanneer de voorgeschreven evolutie de Tweede Wet zou schenden, wat effectief een "elastische kloof" creëert.

3. Computationeel Schema:

   • De auteurs maken gebruik van een Dual-to-Primal (DtP) variationele formulering gebaseerd op een convex duaal variatieprincipe.
   • De methode maakt gebruik van een sequentie van convexe optimalisatieproblemen om de sturende vergelijkingen op te lossen.
   • Het numerieke algoritme verloopt in twee fasen:
     • Gradient Flow: Een gradiëntafdaling-schema in een fictieve tijdsvariabele stuurt de initiële benadering naar de oplossing.
     • Newton-Raphson: Een Newton-Raphson-schema met stapgroottecontrole verfijnt de oplossing.
   • Cruciaal is dat het algoritme het oplossingsdomein (de "DtP-zone") beperkt om de convexiteit van de duale functionaal te waarborgen, waardoor wordt voorkomen dat de oplossing de Tweede Wet schendt.
   • De methode maakt gebruik van eindige-elementen-discretisatie in de tijd voor de duale velden.

Belangrijkste Resultaten
Het schema werd getest op een 1D-staaf onder quasi-statische rek met een voorgeschreven plastische evolutiefunctie die negatief wordt (waardoor de Tweede Wet wordt geschonden) tijdens een specifiek tijdsinterval.

• Correctie van Constitutieve Schending: Het computationele schema slaagde erin de opzettelijk gebrekkige constitutieve specificatie te "corrigeren". Wanneer de voorgeschreven evolutie $f_c$ dreigde een negatieve dissipatie te produceren, werd het correctieve veld $a(t)$ geactiveerd om te garanderen dat de totale dissipatie niet-negatief bleef.
• Elastische Kloven: De activatie van het correctieve veld resulteerde in een plastische rek-snelheid van nul tijdens het interval van de schending. Bij consequentie vertoonde de spanning-rek respons een "elastische kloof" waarbij het materiaal puur elastisch gedrag vertoonde ($d\sigma/d\epsilon = E$), ondanks de toenemende externe belasting. Dit gedrag kwam natuurlijk voort uit de optimalisatie in plaats van te worden voorgeschreven.
• Nauwkeurigheid: De numerieke oplossingen (Eindige Elementen Methode) vertoonden een uitstekende overeenstemming met de analytische gesloten vorm oplossingen.
  • Voor het snelheid-gevoelige geval ($m=1$), was de fout in de spanning-rek relatie lager dan 0,05%.
  • Voor het snelheid-ongevoelige geval ($m=0,1$), was de fout lager dan 0,8%.
  • Fouten waren gelokaliseerd in regio's waar de modulerende functie $g(t)$ snelle transities onderging, wat suggereert dat de fouten voortkomen uit de mesh-resolutie in plaats van het algoritmische schema zelf.
• Rol van Initiële Condities: De studie benadrukte dat de initiële basis-toestand in het numerieke algoritme fungeert als een selectieparameter. Een kleine initiële waarde voor het dissipatieveld leidde de numerieke oplossing naar de specifieke gesloten vorm oplossing die in het analytische gedeelte is afgeleid. Grotere initiële waarden convergeerden naar andere, doch fysiek geldige oplossingen binnen de oneindige familie van oplossingen die door de randvoorwaarden worden toegestaan.

Betekenis en Claims
Het artikel beweert een levensvatbaar schema te hebben aangetoond voor het behandelen van de Dissipatie-ongelijkheid niet als een voorafgaande voorwaarde voor constitutieve modellering, maar als een randvoorwaarde die moet worden voldaan tijdens het oplossen van randwaardeproblemen.

• Robuustheid tegen Onzekerheid: De aanpak staat het gebruik toe van constitutieve modellen die mogelijk incompleet of imperfect zijn, door ze automatisch te corrigeren naar de fundamentele thermodynamische postulaten met een "minimale" afwijking.
• Computationele Haalbaarheid: De auteurs bieden de eerste numerieke bevestiging dat het duale variatieprincipe en het bijbehorende DtP-schema de Tweede Wet kunnen afdwingen in scenario's waar de wet anders mogelijk door het gespecificeerde proces zou worden geschonden.
• Toepassingsbereik: De auteurs suggereren dat de techniek bijzonder relevant is voor complexe mechanische modellen waarbij de Tweede Wet-restricties op constitutieve functies ingewikkeld zijn, of voor het koppelen van uiteenlopende fenomenen waarbij standaard constitutieve restricties moeilijk te formuleren zijn. Zij merken ook op dat de opkomst van "elastische kloven" in de spanning-rek respons overeenkomt met gedragingen waargenomen in strain-gradient plasticiteitstheorieën, hoewel dit wordt gepresenteerd als een emergent kenmerk van het specifieke voorbeeld in plaats van een nieuwe fysieke theorie.

Het werk stelt geen nieuwe experimentele opstellingen voor of claimt alle constitutieve modelleringsproblemen op te lossen, maar biedt eerder een specifiek wiskundig kader om de randvoorwaarde van thermodynamische toelaatbaarheid te beheren in de aanwezigheid van onzekere of potentieel gebrekkige materiaalspecificaties.