Computation of stresses in jammed packings modeled with Tresca friction

Sintesi Tecnica: Computo degli Sforzi in Impaccamenti Intasati Modellati con Attrito di Tresca

Enunciato del Probleamento
Questo articolo affronta il calcolo dei campi di sforzo all'interno di impaccamenti intasati (jammed) di celle poligonali rigide, come strutture murarie o materiali granulari densi. A differenza dei tradizionali Metodi degli Elementi Discreti (DEM) che spesso si basano su deformazioni elastiche per penalizzare l'interpenetrazione, questo lavoro si concentra su sistemi in cui le celle sono completamente rigide e non si deformano elasticamente. In tali configurazioni "intasate", l'equilibrio meccanico è mantenuto esclusivamente attraverso forze di contatto e attrito, senza significative deformazioni elastiche. La sfida specifica risiede nel determinare la distribuzione interna degli sforzi in questi sistemi staticamente indeterminati (iperstatici), dove il numero di variabili di forza supera il numero di vincoli di equilibrio, portando a uno spazio di soluzioni non univoco. Gli autori modellano la fisica del contatto utilizzando una legge di attrito di Tresca, la quale postula una soglia fissa per lo scorrimento tangenziale, distinta dall'attrito di Coulomb dipendente dallo sforzo.

Metodologia
La metodologia proposta segue un framework computazionale in tre fasi:

1. Minimizzazione dell'Energia Primale:
   Gli autori formulano un problema di minimizzazione vincolato in cui viene minimizzata l'energia d'attrito totale. Il funzionale dell'energia consiste nella somma del lavoro d'attrito di Tresca su tutte le interfacce interne, definito come $s_T |[u]_e \cdot t_e|$, dove $s_T$ è il coefficiente d'attrito, $[u]_e$ è il salto di spostamento e $t_e$ è il vettore tangente. Questa minimizzazione è soggetta a:

   • Vincoli di non-interpenetrazione: La componente normale del salto di spostamento deve essere non negativa ($[u]_e \cdot n_e \geq 0$).
   • Confinamento: Viene applicata una condizione al contorno di Neumann ($g$) per evitare che le celle scivolino all'infinito, poiché l'energia manca di forze di ripristino elastico.
   • Invarianza per traslazione: Viene aggiunto un vincolo globale per fissare il centro di massa.
     Le celle sono modellate come corpi rigidi con soli gradi di libertà traslazionali; i gradi di libertà rotazionali sono esplicitamente trascurati. Questa semplificazione è motivata dall'osservazione che, nelle configurazioni intasate, la resistenza alla rotazione (che richiede dilatanza) è di ordini di grandezza superiore alla resistenza allo scorrimento.

2. Massimizzazione Duale:
   Poiché il problema primale cerca gli spostamenti (che sono effettivamente nulli in uno stato intasato stabile), gli autori deducono un problema di massimizzazione duale. Le variabili duali corrispondono alle forze che agiscono sulle interfacce interne. La risoluzione di questo problema duale fornisce gli sforzi normali e tangenziali alle interfacce. L'esistenza di minimizzatori per l'energia primale e di soluzioni per il problema duale è stabilita sotto ipotesi appropriate riguardanti il termine di confinamento. La formulazione duale seleziona efficacemente una soluzione specifica dal manifold iperstatico delle possibili configurazioni di sforzo.

3. Ricostruzione dello Sforzo:
   L'output del problema duale fornisce valori di sforzo discreti agli spigoli. Per ricostruire un campo di sforzo continuo $\sigma$ all'interno di ogni cella, gli autori utilizzano elementi finiti di ordine inferiore Raviart–Thomas ($RT_0$). Questi elementi sono conformi a $H(\text{div})$, il che significa che i loro gradi di libertà corrispondono agli sforzi normali ai baricentri degli spigoli. La ricostruzione prevede la risoluzione di un problema misto locale in ogni cella per trovare un tensore di sforzo $\sigma_c$ tale che:

   • $\text{div}(\sigma_c) = 0$ (equilibrio in assenza di carichi volumetrici).
   • Lo sforzo normale sul bordo corrisponde ai valori ottenuti dalla soluzione duale.
   • La simmetria del tensore di sforzo ($\sigma = \sigma^T$) è imposta debolmente tramite un moltiplicatore di Lagrange. Ciò è necessario perché il modello discreto sottostante trascura le equazioni di equilibrio dei momenti, il che significa che lo sforzo ricostruito può possedere una componente antisimmetrica che rappresenta i momenti residui.

Contributi Chiave

• Modello di Attrito Rigido: Il lavoro introduce un modello per impaccamenti intasati che si basa esclusivamente sull'attrito rigido (criterio di Tresca) senza parametri elastici, distinguendosi dagli approoli DEM standard.
• Gestione dell'Iperstaticità: Il lavoro affronta esplicitamente l'indeterminazione statica degli impaccamenti intasati. Dimostra che l'approccio di ottimizzazione seleziona una soluzione univoca dallo spazio degli stati di auto-sforzo minimizzando il funzionale dell'energia d'attrito.
• Ricostruzione Coerente dello Sforzo: Viene presentata una nuova tecnica di ricostruzione utilizzando elementi $RT_0$ per generare un campo di sforzo coerente a partire dalle forze di interfaccia discrete, garantendo l'equilibrio delle forze locali all'interno di ogni cella.
• Trascuratezza della Rotazione: Gli autori propongono un modello semplificato che ignora i gradi di libertà rotazionali, sostenendo che ciò è fisicamente giustificato per gli impaccamenti compressi dove lo scorrimento traslazionale domina, pur riconoscendo che ciò porta a tensori di sforzo non simmetrici.

Risultati Numerici
La metodologia è stata testata su varie mesh poligonali 2D (specificamente tassellazioni di Voronoi) con le seguenti evidenze:

• Verifica: In test con condizioni al contorno uniformi, il campo di sforzo ricostruito corrispondeva allo stato di sforzo uniforme atteso, confermando la coerenza del metodo.
• Test di Taglio: In uno scenario di carico di taglio, il tensore di sforzo ricostruito ha mostrato asimmetria ($\sigma_{12} \neq \sigma_{21}$). Gli autori attribuiscono ciò alla mancanza di un equilibrio dei momenti imposto nel modello discreto, notando che la parte antisimmetrica tiene conto dei momenti residui.
• Compressione Omogenea: Nonostante le condizioni di compressione uniforme, il solutore ha prodotto un campo di sforzo non omogeneo. Gli autori spiegano questo come una conseguenza dell'iperstaticità del sistema; l'ottimizzazione seleziona uno stato specifico che minimizza l'energia dal nucleo (null space), il che porta naturalmente alla formazione di "catene di forza" (reti sparse di alto sforzo) e regioni di basso o nullo sforzo.
• Tensione in Compressione: I risultati hanno mostrato celle in tensione (sforzi principali positivi) anche sotto compressione globale. Gli autori identificano questo come un fenomeno fisico derivante dall'incastro (interlocking) frizionale e dagli effetti di ponte (bridging), piuttosto che come un artefatto numerico.

Significatività e Rivendicazioni
Il lavoro sostiene che la metodologia proposta offre un modo robusto e coerente per calcolare gli sforzi negli impaccamenti intasati, sfruttando la fisica dell'attrito rigido e la minimizzazione dell'energia. Gli autori sostengono che il loro approccio cattura meglio la fisica sottostante dei sistemi intasati rispetto ai DEM basati sull'elasticità, particolarmente nei regimi in cui la deformazione elastica è trascurabile.

La significatività dei risultati risiede nella dimostrazione che la minimizzazione dell'energia d'attrito (un obiettivo di tipo $L_1$) porta naturalmente a reti di forza sparse, allineandosi alle osservazioni sperimentali delle catene di forza nei materiali granulari. Il lavoro riconosce modestamente i propri limiti: la trascuratezza dell'equilibrio rotazionale produce tensori di sforzo non simmetrici, e l'uso dell'attrito di Tresca (soglia fissa) semplifica il problema rispetto alle più complesse Disuguaglianze Quasi-Variazionali richieste dall'attrito di Coulomb. Gli autori suggeriscono che lavori futuri potrebbero estendere il modello includendo i gradi di libertà rotazionali e leggi d'attrito più realistiche per raggiungere una piena coerenza meccanica.