On the passage from nonlinear to linearized viscoelastodynamics

1. Problema e Contexto

O artigo aborda a derivação rigorosa das equações da viscoelasticidade linearizada (na reologia de Kelvin-Voigt) a partir de um modelo não linear completo que satisfaz a indiferença de quadro dependente do tempo (time-dependent frame indifference), conforme definido por Antman.

O problema central é justificar matematicamente a transição de modelos de grandes deformações (não lineares) para modelos de pequenas deformações (linearizados) em regimes dinâmicos, considerando:

• Forças de inércia: Problemas hiperbólicos completos.
• Quase-estáticos: Onde a inércia é negligenciada.
• Estáticos: Onde o tempo não é uma variável.

O objetivo é demonstrar que a convergência das soluções do sistema não linear para o sistema linearizado ocorre consistentemente em todas as escalas de tempo (dinâmica, quase-estática e estática) e através de diferentes esquemas de aproximação (discretização temporal e limites de parâmetros).

2. Metodologia

Os autores utilizam uma abordagem multifacetada que combina análise variacional, teoria de $\Gamma$-convergência e estimativas de energia a priori.

A. Formulação do Modelo Não Linear

O modelo base considera materiais padrão generalizados (generalized standard materials) com:

• Energia Elástica ($E$): Depende do gradiente de deformação ($\nabla y$) e do gradiente de segunda ordem ($\nabla^2 y$) para garantir regularidade (materiais "não simples" de Toupin).
• Funcional de Dissipação ($R$): Depende da taxa de deformação e satisfaz a indiferença de quadro.
• Equação de Movimento:
  $$\varrho \partial_{tt} y - \text{div } T = \tilde{f}$$
  onde o tensor de tensão $T$ é a soma de partes elásticas e viscosas.

B. Escalonamento e Linearização

Introduz-se um parâmetro pequeno $\delta > 0$ para medir a magnitude das deformações. A deformação é reescalonada como $y(t, x) = x + \delta u(t, x)$, onde $u$ é o deslocamento.
Os funcionais de energia e dissipação são reescalonados ($E_\delta, R_\delta$) de forma que, no limite $\delta \to 0$, os termos de ordem superior (não lineares) desapareçam ou se tornem lineares, resultando nas equações clássicas da elasticidade linear com coeficientes elásticos ($C_W$) e viscosos ($C_D$).

C. Esquema de Aproximação Temporal (Duas Escalas)

Para conectar os regimes dinâmico, quase-estático e estático, os autores empregam um esquema de discretização temporal de duas escalas (inspirado em [BKS24] e De Giorgi):

1. Escala $h$ (Atraso de Tempo): A segunda derivada temporal (inércia) é aproximada por um quociente de diferenças atrasado: $\partial_{tt} y \approx \frac{\partial_t y(t) - \partial_t y(t-h)}{h}$. Isso transforma o problema hiperbólico em um problema de equilíbrio quase-estático com um termo de força dependente do passado.
2. Escala $\tau$ (Discretização Total): O problema contínuo é discretizado em passos de tempo $\tau$, levando a um problema variacional de minimização incremental (movimentos minimizantes).

D. Ferramentas Analíticas

• $\Gamma$-convergência: Usada para provar a convergência dos funcionais de energia e dissipação.
• Desigualdades de Energia: Estabelecimento de estimativas a priori uniformes para garantir a compacidade das soluções.
• Técnicas de Rigidez: Uso de estimativas de rigidez de Friesecke-James-Müller e desigualdades de Korn generalizadas para controlar a norma $W^{1,2}$ dos deslocamentos.
• Convergência EDP (Energy-Dissipation Principle): Uma abordagem variacional que utiliza a desigualdade de Fenchel-Young para caracterizar a solução como um fluxo de gradiente, permitindo provar a convergência sem depender exclusivamente da forma diferencial da equação.

3. Principais Contribuições e Resultados

1. Derivação Rigorosa da Linearização Dinâmica

O artigo prova que, sob condições de "dados bem preparados" (initial data e forças escalonados apropriadamente com $\delta$), as soluções do problema não linear dinâmico convergem para a solução única do problema linearizado dinâmico quando $\delta \to 0$.

2. Comutatividade dos Limites

Um dos resultados mais significativos é a demonstração de que a linearização comuta com a aproximação temporal. O diagrama de convergência apresentado no artigo mostra que:

• Pode-se primeiro passar para o limite $\delta \to 0$ (linearizar) e depois $\tau \to 0$ ou $h \to 0$.
• Ou pode-se primeiro discretizar ($\tau, h$) e depois linearizar ($\delta \to 0$).
• O resultado final é o mesmo: a solução do sistema linearizado correspondente.

3. Convergência em Múltiplos Regimes

O trabalho estabelece a convergência não apenas para o problema dinâmico completo, mas também para:

• O problema de atraso de tempo (time-delayed), que serve de ponte entre o dinâmico e o quase-estático.
• O problema quase-estático.
• O problema estático.
  Isso valida a consistência física e matemática de usar modelos linearizados em diferentes regimes de carregamento.

4. Convergência EDP com Forças

Os autores desenvolvem uma teoria de EDP-convergência com forças (EDP-convergence with forces). Eles mostram que a convergência dos funcionais de energia e dissipação (incluindo o trabalho das forças externas) é suficiente para garantir a convergência das soluções, mesmo na presença de forças dependentes do tempo. Isso oferece uma prova alternativa e robusta baseada puramente na estrutura variacional, sem necessidade de manipular diretamente as equações diferenciais parciais em todos os passos.

5. Tratamento de Materiais "Não Simples"

O modelo inclui termos de gradiente de segunda ordem ($\nabla^2 y$) para garantir a existência de soluções fracas no regime não linear. O artigo demonstra rigorosamente como esses termos de regularização desaparecem ou se tornam irrelevantes no limite linearizado, justificando sua omissão nos modelos clássicos de pequenas deformações.

4. Significado e Impacto

Este trabalho é fundamental para a mecânica dos sólidos e análise matemática por várias razões:

• Validação Teórica: Fornece a primeira prova rigorosa completa que conecta modelos não lineares hiperbólicos complexos com suas contrapartes linearizadas, cobrindo todo o espectro de escalas de tempo (dinâmica, quase-estática, estática).
• Justificativa de Modelos Reduzidos: Confirma matematicamente a prática comum de usar equações linearizadas para pequenas deformações, mesmo em regimes dinâmicos onde a inércia é significativa.
• Metodologia Variacional: A aplicação da $\Gamma$-convergência e da EDP-convergência a problemas hiperbólicos (dinâmicos) é uma contribuição metodológica importante, estendendo técnicas que eram predominantemente usadas em problemas de evolução rate-independent ou gradientes.
• Robustez Numérica: Ao provar a convergência em esquemas de discretização temporal (minimizing movements), o trabalho oferece uma base teórica sólida para métodos numéricos que simulam a transição de grandes para pequenas deformações.

Em resumo, o artigo estabelece que a estrutura variacional subjacente (a convergência dos funcionais de energia e dissipação) é o princípio fundamental que governa a transição de modelos não lineares para lineares em viscoelasticidade, independentemente da escala de tempo ou do método de aproximação utilizado.