Computation of stresses in jammed packings modeled with Tresca friction

Resumo Técnico: Computação de Tensões em Empacotamentos Bloqueados Modelados com Atrito de Tresca

Enunciado do Problema
Este artigo aborda a computação de campos de tensão em empacotamentos bloqueados (jammed packings) de células poligonais rígidas, como estruturas de alvenaria ou materiais granulares densos. Diferente dos métodos tradicionais de Elementos Discretos (DEM), que frequentemente dependem de deformações elásticas para penalizar a interpenetração, este trabalho foca em sistemas onde as células são totalmente rígidas e não se deformam elasticamente. Em tais configurações "bloqueadas", o equilíbrio mecânico é mantido exclusivamente através de forças de contato e atrito, sem deformação elástica significativa. O desafio específico reside em determinar a distribuição de tensão interna nesses sistemas estaticamente indeterminados (hiperestáticos), onde o número de variáveis de força excede o número de restrições de equilíbrio, levando a um espaço de solução não único. Os autores modelam a física de contato usando uma lei de atrito de Tresca, que postula um limiar fixo para o deslizamento tangencial, distinto do atrito de Coulomb dependente da tensão.

Metodologia
A metodologia proposta segue um arcabouço computacional de três estágios:

1. Minimização de Energia Primal:
   Os autores formulam um problema de minimização restrita onde a energia de atrito total é minimizada. O funcional de energia consiste na soma do trabalho de atrito de Tresca sobre todas as interfaces internas, definido como $s_T |[u]_e \cdot t_e|$, onde $s_T$ é o coeficiente de atrito, $[u]_e$ é o salto de deslocamento e $t_e$ é o vetor tangente. Esta minimização é sujeita a:

   • Restrições de não-interpenetração: A componente normal do salto de deslocamento deve ser não-negativa ($[u]_e \cdot n_e \geq 0$).
   • Confinamento: Uma condição de contorno de Neumann ($g$) é aplicada para evitar que as células deslizem ao infinito, uma vez que a energia carece de forças restauradoras elásticas.
   • Invariância de translação: Uma restrição global é adicionada para fixar o centro de massa.
     As células são modeladas como corpos rígidos com apenas graus de liberdade de translação; os graus de liberdade rotacionais são explicitamente negligenciados. Esta simplificação é motivada pela observação de que a resistência à rotação (que requer dilatância) é ordens de magnitude superior à resistência ao deslizamento.

2. Maximização Dual:
   Uma vez que o problema primal busca deslocamentos (que são efetivamente zero em um estado bloqueado estável), os autores deduzem um problema de maximização dual. As variáveis duais correspondem às forças que atuam nas arestas internas. A resolução deste problema dual fornece as tensões normal e tangencial nas interfaces. A existência de minimizadores para a energia primal e soluções para o problema dual é estabelecida sob hipóteses apropriadas em relação ao termo de confinamento. A formulação dual seleciona efetivamente uma solução específica do manufoldo de tensões possíveis hiperestáticas.

3. Reconstrução de Tensão:
   O resultado do problema dual fornece valores de tensão discretos nas arestas. Para reconstruir um campo de tensão contínuo $\sigma$ dentro de cada célula, os autores empregam elementos finitos de ordem mais baixa Raviart–Thomas ($RT_0$). Estes elementos são conformes a $H(\text{div})$, o que significa que seus graus de liberdade correspondem às tensões normais nos baricentros das arestas. A reconstrução envolve a resolução de um problema misto local em cada célula para encontrar um tensor de tensão $\sigma_c$ tal que:

   • $\text{div}(\sigma_c) = 0$ (equilíbrio na ausência de cargas volumétricas).
   • A tensão normal na fronteira coincide com os valores obtidos da solução dual.
   • A simetria do tensor de tensão ($\sigma = \sigma^T$) é imposta fracamente via multiplicador de Lagrange. Isto é necessário porque o modelo discreto subjacente negligencia equações de balanço de momento, o que significa que a tensão reconstruída pode possuir uma componente antissimétrica representando torques residuais.

Principais Contribuições

• Modelo de Atrito Rígido: O artigo introduz um modelo para empacotamentos bloqueados que depende exclusivamente de atrito rígido (critério de Tresca) sem parâmetros elásticos, distinguindo-se das abordagens DEM padrão.
• Tratamento da Hiperestaticidade: O trabalho aborda explicitamente a indeterminação estática dos empacotamentos bloqueados. Demonstra que a abordagem de otimização seleciona uma solução única do espaço de estados de auto-tensão ao minimizar o funcional de energia de atrito.
• Reconstrução de Tensão Consistente: Uma técnica de reconstrução inovadora é apresentada usando elementos $RT_0$ para gerar um campo de tensão consistente a partir de forças de interface discretas, garantindo o equilíbrio de forças local dentro de cada célula.
• Negligência de Rotação: Os autores propõem um modelo simplificado que ignora os graus de liberdade rotacionais, argumentando que isso é fisicamente justificado para empacotamentos comprimidos onde o deslizamento translacional domina, embora reconheçam que isso leva a tensores de tensão não simétricos.

Resultados Numéricos
A metodologia foi testada em várias malhas poligonais 2D (especificamente tesselações de Voronoi) com as seguintes conclusões:

• Verificação: Em testes com condições de contorno uniformes, o campo de tensão reconstruído coincidiu com o estado de tensão uniforme esperado, confirmando a consistência do método.
• Teste de Cisalhamento: Em um cenário de carregamento por cisalhamento, o tensor de tensão reconstruído exibiu assimetria ($\sigma_{12} \neq \sigma_{21}$). Os autores atribuem isso à falta de imposição do balanço de momentos no modelo discreto, observando que a parte antissimétrica contabiliza os torques residuais.
• Compressão Homogênea: Apesar das condições de contorno compressivas uniformes, o solver produziu um campo de tensão não homogêneo. Os autores explicam isto como uma consequência da hiperestaticidade do sistema; a otimização seleciona um estado específico de minimização de energia do espaço nulo, o que naturalmente leva à formação de "cadeias de força" (redes esparsas de alta tensão) e regiões de baixa ou zero tensão.
• Tração em Compressão: Os resultados mostraram células em tração (tensões principais positivas) mesmo sob compressão global. Os autores identificam isto como um fenômeno físico decorrente do travamento por atrito (frictional interlocking) e efeitos de ponte (bridging), e não como um artefato numérico.

Significância e Alegações
O artigo afirma que a metodologia proposta oferece uma maneira robusta e consistente de computar tensões em empacotamentos bloqueados, aproveitando a física do atrito rígido e a minimização de energia. Os autores argumentam que sua abordagem captura melhor a física subjacente dos sistemas bloqueados em comparação com os DEMs baseados em elasticidade, particularmente em regimes onde a deformação elástica é desprezível.

A significância dos resultados reside na demonstração de que a minimização da energia de atrito (um objetivo do tipo $L_1$) leva naturalmente a redes de força esparsas, alinhando-se com as observações experimentais de cadeias de força em materiais granulares. O artigo reconhece modestamente suas limitações: a negligência do equilíbrio rotacional resulta em tensores de tensão não simétricos, e o uso do atrito de Tresca (limiar fixo) simplifica o problema em comparação com as Desigualdades Variacionais Quasi-Não Convexas mais complexas exigidas pelo atrito de Coulomb. Os autores sugerem que trabalhos futuros podem estender o modelo para incluir graus de liberdade rotacionais e leis de atrito mais realistas para alcançar plena consistência mecânica.