Computation of stresses in jammed packings modeled with Tresca friction

Technische Samenvatting: Berekening van Spanningen in Verende Pakkingen Gemodelleerd met Tresca-Wrijving

Probleemstelling
Dit artikel behandelt de berekening van spanningsvelden in verende (jammed) pakkingen van rigide polygoonvormige cellen, zoals metselwerkstructuren of dichte granulaire materialen. In tegenstelling tot traditionele Discrete Element Methods (DEM), die vaak vertrouwen op elastische deformaties om penetratie te bestraffen, richt dit werk zich op systemen waarbij cellen volledig rigide zijn en niet elastisch vervormen. In dergelijke "verende" configuraties wordt het mechanisch evenwicht uitsluitend onderhouden door contactkrachten en wrijving, zonder significante elastische rek. De specifieke uitdaging ligt in het bepalen van de interne spanningsverdeling in deze statisch onbepaalde (hyperstatische) systemen, waarbij het aantal krachtvariabelen groter is dan het aantal evenwichtsbeperkingen, wat leidt tot een niet-unieke oplossingsruimte. De auteurs modelleren de contactfysica met een Tresca-wrijvingswet, die een vaste drempelwaarde voor tangentiële glijding postuleert, in tegenstelling tot de stressafhankelijke Coulomb-wrijving.

Methodologie
De voorgestelde methodologie volgt een driedelige computationele procedure:

1. Primal Energie-minimalisatie:
   De auteurs formuleren een beperkt minimalisatieprobleem waarbij de totale wrijvingsenergie wordt geminimaliseerd. De energiefunctie bestaat uit de som van de Tresca-wrijvingsarbeid over alle interne interfaces, gedefinieerd als $s_T |[u]_e \cdot t_e|$, waarbij $s_T$ de wrijvingscoëfficiënt is, $[u]_e$ de verplaatsingssprong en $t_e$ de tangentiële vector. Deze minimalisatie is onderworpen aan:

   • Niet-interpenetratie beperkingen: De normale component van de verplaatsingssprong moet niet-negatief zijn ($[u]_e \cdot n_e \geq 0$).
   • Opsluiting: Een Neumann-randvoorwaarde ($g$) wordt toegepast om te voorkomen dat cellen naar oneindig glijden, aangezien de energie geen elastische herstelkrachten bezit.
   • Translatie-invariantie: Een globale beperking wordt toegevoegd om het zwaartepunt vast te leggen.
     De cellen worden gemodelleerd als rigide lichamen met enkel translationele vrijheidsgraden; rotationele vrijheidsgraden worden expliciet genegeerd. Deze vereenvoudiging wordt gemotiveerd door de observatie dat de weerstand tegen rotatie (die dilatantie vereist) orde van grootte hoger is dan de weerstand tegen glijden in verende configuraties.

2. Duale Maximalisatie:
   Omdat het primaire probleem naar verplaatsingen streeft (die in een stabiele verende toestand effectief nul zijn), leiden de auteurs een duaal maximalisatieprobleem af. De duale variabelen komen overeen met de krachten die op de interne randen werken. Het oplossen van dit duale probleem levert de normale en tangentiële spanningen op de interfaces op. Het bestaan van minimizers voor de primaire energie en oplossingen voor het duale probleem wordt vastgesteld onder passende hypothesen met betrekking tot de opsluitingsterm. De duale formulering selecteert effectief een specifieke oplossing uit de hyperstatische manifold van mogelijke spanningsstaten.

3. Spanningsreconstructie:
   De output van het duale probleem levert discrete spanningswaarden op de randen. Om een continu spanningsveld $\sigma$ binnen elke cel te reconstrueren, gebruiken de auteurs lowest-order Raviart–Thomas ($RT_0$) eindige elementen. Deze elementen zijn $H(\text{div})$-conforme, wat betekent dat hun vrijheidsgraden overeenkomen met de normale spanningen op de barycentra van de randen. De reconstructie omvat het oplossen van een lokaal gemengd probleem in elke cel om een spannings-tensor $\sigma_c$ te vinden zodat:

   • $\text{div}(\sigma_c) = 0$ (evenwicht in de afwezigheid van volumelasten).
   • De normale spanning op de rand overeenkomt met de waarden verkregen uit de duale oplossing.
   • De symmetrie van de spannings-tensor ($\sigma = \sigma^T$) wordt zwak opgelegd via een Lagrange-multiplier. Dit is noodzakelijk omdat het onderliggende discrete model momentenbalansvergelijkingen verwaarloost, wat betekent dat de gereconstrueerde spanning een asymmetrisch component kan bezitten die residuele momenten vertegenwoordigt.

Belangrijkste Bijdragen

• Rigide Wrijvingsmodel: Het artikel introduceert een model voor verende pakkingen dat uitsluitend steunt op rigide wrijving (Tresca-criterium) zonder elastische parameters, wat het onderscheidt van standaard DEM-benaderingen.
• Afhandeling van Hyperstatische Onbepaaldheid: Het werk adresseert expliciet de statische onbepaaldheid van verende pakkingen. Het toont aan dat de optimalisatiebenadering een unieke oplossing selecteert uit de ruimte van zelf-spanningsstaten door de wrijvingsenergie-functionaal te minimaliseren.
• Consistente Spanningsreconstructie: Een nieuwe reconstructietechniek wordt gepresenteerd met behulp van $RT_0$-elementen om een consistent spanningsveld te genereren uit discrete interfacekrachten, waarbij lokaal krachtevenwicht binnen elke cel wordt gewaarborgd.
• Verwaarlozing van Rotatie: De auteurs stellen een vereenvoudigd model voor dat rotationele vrijheidsgraden negeert, en beargumenteren dat dit fysiek gerechtvaardigd is voor gecomprimeerde pakkingen waar translationeel glijden domineert, terwijl ze erkennen dat dit leidt tot niet-symmetrische spannings-tensoren.

Numerieke Resultaten
De methodologie werd getest op diverse 2D-polygoonmeshes (specifiek Voronoi-tessellaties) met de volgende bevindingen:

• Verificatie: In tests met uniforme randvoorwaarden kwam het gereconstrueerde spanningsveld overeen met de verwachte uniforme spanningsstaat, wat de consistentie van de methode bevestigt.
• Schuiftest: In een scenario met schuifbelasting vertoonde de gereconstrueerde spannings-tensor asymmetrie ($\sigma_{12} \neq \sigma_{21}$). De auteurs schrijven dit toe aan het gebrek aan afgedwongen momentenbalans in het discrete model, waarbij zij opmerken dat het asymmetrische deel de residuele momenten account.
• Homogene Compressie: Ondanks uniforme compressieve randvoorwaarden produceerde de solver een niet-homogeen spanningsveld. De auteurs verklaren dit als een gevolg van de hyperstatische onbepaaldheid van het systeem; de optimalisatie selecteert een specifieke energie-minimaliserende staat uit de nulruimte, wat van nature leidt tot de vorming van "force chains" (ijle netwerken van hoge spanning) en regio's met lage of nul spanning.
• Trekspanning in Compressie: De resultaten toonden cellen met trekspanning (positieve hoofdspanningen) zelfs onder globale compressie. De auteurs identificeren dit als een fysisch fenomeen dat voortkomt uit wrijvingsmatige inklemming en brugvorming (bridging), in plaats van een numeriek artefact.

Betekenis en Claims
Het artikel claimt dat de voorgestelde methodologie een robuuste en consistente manier biedt om spanningen in verende pakkingen te berekenen door gebruik te maken van de fysica van rigide wrijving en energie-minimalisatie. De auteurs beargumenteren dat hun benadering de onderliggende fysica van verende systemen beter vangt dan elastisch gebaseerde DEM's, met name in regimes waar elastische deformatie verwaarloosbaar is.

De betekenis van de resultaten ligt in de demonstratie dat het minimaliseren van wrijvingsenergie (een $L_1$-type objectief) van nature leidt tot ijle krachtnetwerken, wat overeenkomt met experimentele observaties van krachtketens in granulaire materialen. De auteurs erkennen bescheiden de beperkingen van hun werk: het negeren van rotationele evenwicht resulteert in niet-symmetrische spannings-tensoren, en het gebruik van Tresca-wrijving (vaste drempelwaarde) vereenvoudigt het probleem ten opzichte van de complexere, niet-convexe Quasi-Variational Inequalities die vereist zijn voor Coulomb-wrijving. De auteurs suggereren dat toekomstig werk het model kan uitbreiden door rotationele vrijheidsgraden en realistischere wrijvingswetten op te nemen om volledige mechanische consistentie te bereiken.