A Projection Method for an Elasto-plasticity Model with Linear Kinematic Hardening

1. Problema Investigado

O artigo aborda um sistema dinâmico de elastoplasticidade que incorpora dois mecanismos físicos fundamentais:

• Viscosidade de Kelvin-Voigt: Introduzida para regularizar o problema e garantir controle espacial da velocidade.
• Endurecimento Cinemático Linear (Tipo Melan-Prager): Um modelo que descreve o fenômeno de endurecimento do material, onde o limite de escoamento (superfície de escoamento) se traduz no espaço de tensões devido à deformação plástica acumulada.

Formulação Matemática:
O modelo é formulado em um quadro variacional onde o conjunto de restrições para a tensão $\sigma$ evolui no tempo e é traduzido por uma variável interna chamada tensão de retorno (backstress), denotada por $\alpha$.

• A relação constitutiva envolve a decomposição da deformação total $E(u)$ em partes elástica ($\varepsilon^e$) e plástica ($\varepsilon^p$).
• A regra de fluxo plástico é acoplada à equação de movimento através de uma estrutura quasi-variacional, pois o conjunto admissível $K(t; \alpha) = \tilde{K}(t) + \alpha(t)$ depende da variável desconhecida $\alpha$.
• O critério de escoamento é do tipo von Mises, definido por $|\tau^D| \leq g(t, \cdot)$, onde $g$ é um limite de escoamento que pode variar no tempo e no espaço.

O sistema de equações diferenciais parciais (EDPs) acopladas inclui:

1. Equação de movimento (conservação do momento) com termos viscosos e elásticos.
2. Inclusão variacional para a evolução da tensão (regra de fluxo).
3. Equação de evolução para a variável interna $\alpha$ (relacionada à taxa de deformação plástica).

2. Metodologia

Os autores empregam o Método de Rothe (discretização temporal primeiro) combinado com uma discretização baseada em projeção.

Esquema Numérico Proposto (P):
Para cada passo de tempo $n$, o problema é dividido em três etapas sequenciais:

1. Etapa de Tentativa (Trial Step): Resolve-se um subproblema linear visco-elástico para obter uma tensão de tentativa $\sigma^*_n$ e a velocidade $v_n$, ignorando temporariamente a restrição de plasticidade.
2. Atualização da Tensão de Retorno: A variável interna $\alpha_n$ é atualizada baseada na taxa de deformação plástica estimada.
3. Etapa de Projeção: A tensão de tentativa, deslocada pela tensão de retorno ($\sigma^*_n - \alpha_n$), é projetada no conjunto admissível transladado $\tilde{K}(t_n)$. Isso é feito utilizando um operador de projeção pontual $P_{g_n}$, que garante que a tensão final $\sigma_n$ satisfaça o critério de escoamento.

Análise de Solvabilidade:
O artigo demonstra que, apesar da não linearidade acoplada entre $\sigma_n$ e $\alpha_n$ na etapa de projeção, o sistema pode ser resolvido explicitamente ponto a ponto (quase em todo $\Omega$) devido à estrutura geométrica do problema no espaço das tensões desviadoras.

3. Principais Contribuições

• Abordagem de Projeção para Endurecimento Cinemático: Diferente de trabalhos anteriores que tratavam o conjunto de restrições como dado externamente, este trabalho lida com o caso onde o conjunto depende da solução (estrutura quasi-variacional) e propõe um esquema de projeção específico para lidar com a tradução do conjunto admissível.
• Relaxamento de Hipóteses sobre o Limite de Escoamento ($g$): O resultado de existência é estabelecido sob condições muito gerais para a função $g(t, x)$:
  • Não é necessária continuidade espacial.
  • Não é assumido um limite inferior estritamente positivo (ou seja, $g$ pode ser zero em algumas regiões).
• Fundação para Aproximação Numérica: A estrutura do esquema proposto é consistente com a discretização espacial. Os autores destacam que, ao substituir os espaços de funções por espaços de elementos finitos, o mesmo esquema fornece uma base direta para a construção de algoritmos numéricos implementáveis (FEM).

4. Resultados Principais

• Estabilidade do Esquema Discreto (Teorema 4.1):
  Os autores provam estimativas de estabilidade uniformes em relação ao passo de tempo $\Delta t$. As soluções discretas são limitadas em normas apropriadas (incluindo $L^\infty(0,T; H)$ para tensões e velocidades, e $L^2(0,T; V)$ para gradientes de velocidade). Isso garante que o esquema não exploda numericamente.
• Existência de Solução Fraca (Teorema 4.2):
  Utilizando as estimativas de estabilidade, compacidade (Teorema de Aubin-Lions e teorema de Ascoli-Arzelà) e passando ao limite quando $\Delta t \to 0$, prova-se a existência de uma solução fraca para o Problema 1.1 no sentido variacional.
• Unicidade da Solução:
  A unicidade é demonstrada através de um argumento de energia, comparando duas soluções possíveis e mostrando que a diferença entre elas deve ser nula, graças à monotonicidade do operador de subdiferencial e à estrutura do endurecimento linear.

5. Significado e Conclusão

Este trabalho é significativo por fornecer uma base teórica rigorosa para modelos de elastoplasticidade dinâmica com endurecimento cinemático, que são cruciais para simular o comportamento de metais sob carregamento cíclico (efeito Bauschinger).

• Rigor Matemático: A prova de existência e unicidade sob hipóteses fracas para o limite de escoamento preenche uma lacuna na literatura, permitindo a modelagem de cenários mais realistas onde o material pode não ter um limite de escoamento definido em certas regiões ou instantes.
• Ponte para a Computação: Ao utilizar o método de Rothe com projeção, o artigo não apenas prova a existência teórica, mas também oferece um algoritmo construtivo. Isso facilita a transição direta para simulações numéricas usando Elementos Finitos, algo que é frequentemente desafiador em problemas de plasticidade com restrições variáveis.
• Futuro: O artigo sugere que a próxima etapa natural é a análise de convergência de esquemas totalmente discretos (tempo e espaço) e a extensão para outros critérios de escoamento (como Tresca) ou para o caso não regularizado ($\eta = 0$), que apresenta desafios analíticos maiores devido à perda de controle espacial da velocidade.

Em resumo, o artigo estabelece um método robusto e matematicamente fundamentado para resolver sistemas de elastoplasticidade dinâmica complexos, combinando análise funcional avançada com uma estratégia numérica prática.