Caveats on formulating finite elasto-plasticity in curvilinear coordinates

Este artigo apresenta uma metodologia prática e passo a passo para formular a elastoplasticidade de grandes deformações em coordenadas curvilíneas, esclarecendo o tratamento rigoroso de termos adicionais como o gradiente de deformação, o jacobiano e o deslocador em bases não-cartesianas para garantir a implementação robusta em elementos finitos de problemas axisimétricos.

O essencial

Imagine que você é um arquiteto projetando um arranha-céu. Se o prédio for um cubo perfeito, calcular como ele se deforma sob o vento é relativamente fácil: você usa réguas e ângulos retos (como num sistema de coordenadas cartesianas, o "x, y, z" clássico).

Mas e se o prédio for um silos de grãos, um tanque de água ou uma torre de resfriamento? Eles são redondos, cilíndricos. Se você tentar usar as mesmas réguas retas para medir um objeto redondo, as coisas ficam estranhas. É como tentar medir a circunferência de uma pizza usando apenas uma régua reta: você vai ter que fazer muitas contas extras e, se não tomar cuidado, vai errar o tamanho da fatia.

Este artigo científico é basicamente um manual de instruções para engenheiros que precisam fazer esses cálculos complexos em objetos redondos (simetria axial) que sofrem grandes deformações (como um pneu sendo espremido ou um solo sendo perfurado), mas que querem economizar tempo de computador.

Aqui está a explicação "traduzida" para o português, usando analogias do dia a dia:

1. O Problema: A Armadilha da "Régua Curva"

Os computadores de engenharia (como o Método dos Elementos Finitos) são ótimos em resolver problemas usando matrizes (tabelas de números). Normalmente, eles funcionam perfeitamente em sistemas "retos" (Cartesianos).

O problema surge quando tentamos aplicar essa lógica em sistemas "curvos" (como cilindros).

• A Analogia: Imagine que você está desenhando um mapa. Num mapa plano (Cartesiano), se você andar 1 km para o norte e 1 km para o leste, você sabe exatamente onde está. Mas num globo terrestre (Curvilíneo), se você andar 1 km para o norte e 1 km para o leste, a distância que você percorreu "para o leste" depende de quão perto você está do Polo Norte.
• O Erro Comum: Muitos programas de computador, ao tentarem simular objetos redondos, simplesmente copiam as fórmulas do mundo "reto" e as aplicam no mundo "curvo". Isso é como usar um mapa plano para navegar num globo: você acaba no lugar errado. O artigo mostra que, se você não corrigir essas fórmulas, o computador pode achar que o material está esticando ou comprimindo de formas que não existem na realidade.

2. As Ferramentas Secretas: O "Deslocador" e o "Jacobian"

Para consertar isso, os autores introduzem dois conceitos que são como "óculos corretivos" para o computador:

• O Deslocador (Shifter):
  • O que é: É uma ferramenta que ajusta a "régua" quando você muda de um ponto para outro num espaço curvo.
  • A Analogia: Imagine que você está andando em uma esteira rolante que gira. Se você segura um objeto e tenta movê-lo para frente, a esteira gira e muda a direção do objeto. O "Deslocador" é como um assistente que gira o objeto de volta para a direção correta antes de você medir. Sem ele, o computador acha que o objeto mudou de direção magicamente.
• O Jacobiano (O Medidor de Volume):
  • O que é: É um número que diz quanto o volume de um pedaço de material mudou.
  • A Analogia: Pense em uma massa de pão. Se você esticar o pão, ele fica mais fino, mas o volume total (a quantidade de massa) deve ser preservado (ou mudar de forma previsível). Em coordenadas curvas, o "cubo" de massa não é um cubo perfeito, é um pedaço de fatia de pizza. O Jacobiano é o fator que diz ao computador: "Ei, esse pedaço de pizza tem um tamanho diferente do cubo que você imaginou, ajuste o cálculo!".

3. O Desafio da "Plasticidade" (Quando o material não volta ao normal)

O artigo foca em materiais que sofrem deformação plástica.

• A Analogia: Imagine amassar uma lata de refrigerante. Se você soltar, ela não volta ao formato original (isso é plástico). Se você esticar um elástico, ele volta (isso é elástico).
• O Problema: Quando o material é amassado (plástico) e depois esticado (elástico), o computador precisa saber exatamente quanto do volume mudou por causa do amassado e quanto por causa do esticamento. Em coordenadas curvas, separar essas duas coisas é como tentar separar a manteiga derretida da água em um copo: é difícil e requer uma receita muito específica. O artigo fornece essa "receita" matemática para que o computador não se confunda.

4. A Solução: O "Corte" Inteligente (Simetria Axial)

A grande vantagem que o artigo propõe é usar a simetria axial.

• A Analogia: Se você tem um cilindro perfeito e o empurra de forma uniforme, ele continua sendo um cilindro. Você não precisa calcular o que acontece em todo o cilindro (360 graus). Você só precisa calcular uma "fatia" de 2D (como cortar uma fatia de bolo) e dizer ao computador: "Repita isso em volta do centro".
• O Ganho: Isso economiza um tempo computacional enorme. Em vez de calcular milhões de pontos em 3D, o computador calcula apenas milhares em 2D. Mas, para fazer isso funcionar sem erros, você precisa usar as "ferramentas secretas" (Deslocador e Jacobiano) explicadas acima.

5. O Resultado: O Teste do "Cilindro Espesso"

Os autores testaram sua teoria simulando um cilindro oco sendo expandido (como um balão sendo inflado dentro de um tubo).

• Eles compararam o resultado do "Corte Inteligente" (2D com as correções certas) com uma simulação completa e pesada em 3D.
• A Conclusão: Os resultados foram idênticos! O método "cortado" foi tão preciso quanto o método completo, mas muito mais rápido. Isso prova que, se você usar as fórmulas corretas para coordenadas curvas, pode simular grandes deformações em objetos redondos com muito menos esforço de computador.

Resumo Final

Este artigo é um aviso e um guia para engenheiros: "Não use fórmulas de mundo reto para objetos redondos sem ajustes!"

Ele ensina como adicionar os "óculos corretivos" matemáticos necessários para que os computadores entendam a geometria curva. Se você ignorar esses ajustes, seus cálculos de tensão, deformação e volume estarão errados, o que pode levar a falhas no projeto de pontes, túneis ou até implantes médicos. Com as correções certas, podemos simular o mundo real (que é cheio de curvas) com precisão e eficiência.

Resumo técnico

Resumo Técnico: Cuidados na Formulação de Elasto-Plasticidade de Grandes Deformações em Coordenadas Curvilíneas

1. Problema e Contexto

A mecânica dos meios contínuos é fundamentalmente baseada na análise tensorial, que é invariante ao referencial. No entanto, as implementações numéricas (como o Método dos Elementos Finitos - MEF) dependem de representações matriciais em bases escolhidas pelo utilizador.

• O Desafio: Quando se utilizam bases não-Cartesianas e não-ortonormais (como coordenadas curvilíneas, ex: cilíndricas), surgem termos adicionais que são normalmente ausentes ou implícitos nas formulações Cartesianas.
• A Lacuna: A extensão ingênua de formulações Cartesianas para sistemas curvilíneos, especialmente em problemas de grandes deformações (finite strain) e elasto-plasticidade, pode levar a erros fundamentais na definição de quantidades cinemáticas (como o gradiente de deformação e o Jacobiano) e, consequentemente, na tensão de Cauchy e na linearização consistente necessária para a convergência do MEF.
• Aplicações: Problemas com simetria axial são comuns em engenharia geotécnica (fundações, túneis, ensaios CPT), engenharia civil (torres de arrefecimento, silos) e metalurgia (estampagem). Explorar essa simetria reduz drasticamente o custo computacional (de 3D para 2D), mas exige uma formulação rigorosa.

2. Metodologia

O artigo não adota uma abordagem de variedades diferenciais (Riemannianas), optando por trabalhar dentro de representações Cartesianas padrão e utilizando mudanças de base explícitas para obter as formas curvilíneas. O estudo foca na implementação prática de códigos MEF.

• Problema Modelo: O trabalho utiliza o problema de expansão de cavidade (expansão radial de um cilindro oco) como caso de teste inicial para expor as definições e o papel de quantidades críticas.
• Quantidades Cinemáticas Chave:
  • Gradiente de Deformação ($F$): O artigo clarifica que, em coordenadas curvilíneas, os componentes de $F$ não são simplesmente derivadas parciais. É necessário distinguir entre o gradiente de deformação e o gradiente de um vetor.
  • Tensor Deslocador (Shifter - $S$): Introduzido para mapear vetores entre configurações de referência e atual quando as bases mudam. Em coordenadas Cartesianas, este tensor é trivial (identidade), mas em coordenadas cilíndricas, possui componentes não triviais que afetam a definição de deslocamentos e o gradiente de deformação.
  • Jacobiano ($J$): A relação entre volumes infinitesimais em coordenadas curvilíneas requer o uso dos coeficientes métricos (determinante da métrica). O artigo deriva a relação correta: $J = \det(F) \sqrt{\det(g)} / \sqrt{\det(G)}$.
• Elasto-Plasticidade de Grandes Deformações:
  • Utiliza-se a decomposição multiplicativa do gradiente de deformação ($F = F_e \cdot F_p$).
  • Discute-se a definição de Jacobianos elástico ($J_e$) e plástico ($J_p$). O artigo enfatiza que, para deformações plásticas isocóricas, a tensão de Cauchy depende da configuração intermediária livre de tensões, o que complica a definição de $J_p$ em coordenadas curvilíneas.
  • Propõe-se uma formulação para a linearização consistente do Jacobiano elástico baseada na tensão logarítmica elástica ($J_e = \exp(\text{tr}(\epsilon^e))$), minimizando componentes extras na matriz tangente.
• Formulação Fraca e Discretização:
  • Derivação das formas fortes e fracas do equilíbrio de momento linear em coordenadas curvilíneas gerais.
  • Particularização para o caso axissimétrico, reduzindo a integral de volume 3D para uma integral 2D no plano R-Z, com tratamento cuidadoso dos termos fora do plano (componentes $\Theta$).
  • Desenvolvimento de fluxos de trabalho algorítmicos para formulações Lagrangiana Atualizada (UL) e Lagrangiana Total (TL).

3. Contribuições Principais

1. Clarificação Cinemática: Demonstra que o componente "fora do plano" do gradiente de deformação em coordenadas cilíndricas para um problema de expansão de cavidade é unitário (e não zero ou outro valor), corrigindo mal-entendidos comuns na literatura.
2. Tratamento do Tensor Deslocador: Mostra explicitamente como o tensor deslocador ($S$) deve ser utilizado para definir corretamente os deslocamentos e o gradiente de deformação em bases curvilíneas, evitando erros de indexação.
3. Definição de Jacobianos: Fornece uma metodologia clara para separar e calcular os Jacobianos elástico e plástico em contextos de grandes deformações, essencial para materiais que não são perfeitamente incompressíveis.
4. Linearização Consistente: Apresenta uma derivação passo a passo para a linearização das equações de equilíbrio, crucial para a convergência quadrática do método de Newton-Raphson em implementações MEF.
5. Validação Numérica: Desenvolve e valida uma implementação MEF axissimétrica completa, comparando-a com uma análise 3D completa.

4. Resultados

O estudo foi validado através de um exemplo numérico de um cilindro de parede espessa submetido a deslocamentos radiais prescritos na superfície interna, utilizando um material hiperelástico (Neo-Hookean) com endurecimento plástico (modelo Cam-Clay Modificado).

• Comparação 3D vs. 2D Axissimétrico: Os resultados das quantidades invariantes (pressão de Cauchy, tensão desviadora, Jacobianos elástico e plástico) obtidos pela formulação axissimétrica proposta foram idênticos aos obtidos pela análise 3D completa, dentro da precisão numérica.
• Convergência: O esquema de Newton-Raphson demonstrou convergência robusta, com o número de iterações necessárias sendo muito semelhante entre as formulações 3D e 2D (máximo de 7 iterações, média de ~3.5).
• Eficiência: A formulação axissimétrica atingiu resultados de qualidade 3D com uma discretização 2D, confirmando a redução significativa de custo computacional sem perda de precisão.
• Detecção de Erros: O estudo mostrou que definições incorretas das quantidades cinemáticas levariam a inconsistências entre as formulações Lagrangiana Total e Atualizada, comprometendo a solução.

5. Significância

Este trabalho é fundamental para pesquisadores e engenheiros que desenvolvem códigos de elementos finitos para problemas de grandes deformações com simetria axial.

• Rigor Matemático: Elimina ambiguidades sobre como tratar gradientes, transpostos e Jacobianos em bases não-ortonormais, prevenindo erros sutis que podem invalidar simulações complexas.
• Aplicabilidade Prática: Oferece um guia passo a passo e fluxos de algoritmos (Algoritmos 1 e 2 no artigo) para a implementação direta em softwares de MEF.
• Confiança em Simulações: Garante que a exploração da simetria axial para economizar recursos computacionais não introduza erros físicos ou numéricos, especialmente em problemas críticos de engenharia geotécnica e metalurgia onde grandes deformações e plasticidade são predominantes.

Em suma, o artigo estabelece que, embora a mecânica dos meios contínuos seja invariante, a sua implementação numérica em coordenadas curvilíneas exige cuidados extras com a métrica e a mudança de base, e fornece as ferramentas necessárias para realizar essas implementações de forma correta e robusta.