Nonlinear Incompressible Shear Wave Models in Hyperelasticity and Viscoelasticity Frameworks, with Applications to Love Waves

Este artigo apresenta equações gerais para ondas de cisalhamento não lineares em materiais hiperelásticos e viscoelásticos incompressíveis, aplicando-as à análise de ondas do tipo Love em interfaces e demonstrando, através de simulações numéricas, como a velocidade variável dessas ondas evolui ao longo do tempo em relação às propriedades mecânicas das camadas.

O essencial

Imagine que você está tentando entender como as ondas de um terremoto se movem através da Terra, ou como um gelatina treme quando você a cutuca. A maioria das pessoas e cientistas, historicamente, olhava para isso como se o mundo fosse feito de "elásticos perfeitos" ou "molas rígidas". Se você esticasse um pouco, eles voltavam ao normal de forma previsível e linear. É como puxar um elástico de escritório: quanto mais puxa, mais ele resiste, mas sempre de forma reta e simples.

Mas a realidade é mais complexa. Materiais como borracha, tecidos biológicos, rochas profundas e até o nosso próprio corpo não são elásticos perfeitos. Eles têm "personalidade". Quando você os estica muito, eles mudam de comportamento, ficam mais duros ou mais moles de formas inesperadas. Além disso, eles não são perfeitos; eles perdem um pouco de energia na forma de calor (como quando você esfrega as mãos e elas esquentam). Isso é chamado de viscoelasticidade.

Este artigo é como um manual de instruções avançado para entender como as ondas de cisalhamento (ondas que fazem o material "deslizar" lateralmente, como quando você sacode um tapete) se comportam nesses materiais "reais" e não perfeitos.

Aqui está a explicação passo a passo, usando analogias do dia a dia:

1. O Cenário: Duas Camadas de "Massa"

Imagine dois tipos de massa de bolo empilhados.

• A camada de cima: Uma massa mais leve e macia (como a crosta terrestre).
• A camada de baixo: Uma massa mais densa e dura (como o manto da Terra).
• A Interface: A linha onde elas se tocam.

Quando você dá um "soco" (uma explosão ou terremoto) em uma dessas camadas, ondas viajam. Se a camada de baixo for muito mais rígida que a de cima, as ondas podem ficar "presas" na interface, viajando como um trem em trilhos. Isso é o que chamamos de Ondas de Love.

2. O Problema: A Física "Linear" vs. A Realidade "Não Linear"

A física tradicional (linear) diz: "Se você dobrar a força, a onda dobra de tamanho e viaja na mesma velocidade". É como dirigir um carro em uma estrada vazia: você pisa no acelerador, ele vai mais rápido, mas a estrada não muda.

Os autores deste artigo dizem: "Espera aí! Materiais reais mudam a estrada enquanto você dirige".

• Não Linearidade: Se você esticar muito o material, ele pode ficar mais rígido ou mais macio. Isso faz com que a velocidade da onda mude dependendo de quão forte ela é. Uma onda gigante viaja em velocidade diferente de uma onda pequena.
• Viscoelasticidade: É como se o material fosse um melado. Ele não apenas estica, ele "gruda" um pouco e perde energia. Isso faz a onda diminuir de tamanho (amortecer) enquanto viaja.

3. A Grande Descoberta: Novas Equações para Ondas Reais

Os autores criaram novas equações matemáticas (como receitas de bolo complexas) para prever exatamente como essas ondas se comportam quando:

1. O material é "hiperelástico" (pode se deformar muito sem quebrar, como borracha).
2. O material tem "viscosidade" (perde energia, como melado).
3. A onda é forte o suficiente para mudar o comportamento do material (não linear).

Eles usaram um modelo específico chamado Yeoh (que é como uma receita matemática para descrever borrachas e tecidos) e adicionaram ingredientes de "atrito" (viscosidade).

4. O Que Eles Descobriram (A Simulação)

Como é difícil fazer esses experimentos na vida real (quem quer explodir camadas da Terra?), eles usaram computadores para simular o cenário.

• O Experimento: Eles simularam uma "explosão" (um ponto de origem) em diferentes lugares (perto da interface ou longe dela) e viram como a onda se espalhou.
• A Regra de Ouro: Eles descobriram que, mesmo com toda essa complexidade (ondas gigantes, materiais estranhos, perda de energia), as ondas de Love ainda seguem uma regra básica: elas só existem se a velocidade da onda estiver entre a velocidade da camada de cima e a da camada de baixo. É como se a onda precisasse de um "corredor" específico para viajar.
• O Comportamento a Longo Prazo: No início, a onda é bagunçada e rápida. Mas, com o tempo, ela se estabiliza e viaja na velocidade máxima permitida pelo material mais rápido das duas camadas. É como se a onda, cansada de tentar ir rápido demais, decidisse seguir o fluxo mais rápido disponível.

5. A Analogia do Trânsito

Pense nas ondas de Love como carros em uma rodovia de duas pistas:

• Pista 1 (Superfície): Carros lentos.
• Pista 2 (Profundidade): Carros rápidos.
• A Interface: A linha divisória.

Se você tentar entrar na pista rápida vindo da lenta, você só consegue se sua velocidade for compatível. Se você for muito rápido, você "reflete" de volta. Se for muito devagar, você fica preso na pista lenta.
Os autores mostraram que, mesmo que os carros sejam "elásticos" (podem esticar e mudar de forma) e tenham "freios" (viscosidade), eles ainda precisam respeitar essa regra de velocidade para conseguir viajar juntos na interface sem se desintegrar.

6. Por Que Isso Importa?

• Terremotos: Ajuda a entender melhor como a energia de um terremoto se propaga, o que é crucial para prever danos em edifícios.
• Medicina: Pode ajudar a entender como ondas sonoras se movem através de tecidos biológicos (como tendões ou músculos) para diagnósticos mais precisos.
• Exploração: Ajuda a encontrar petróleo ou minerais entendendo como as ondas se comportam em diferentes camadas de rocha.

Resumo Final

Este artigo é como atualizar o mapa de um GPS. O GPS antigo (física linear) dizia: "Vá em linha reta, a velocidade é constante". O novo GPS (física não linear e viscoelástica) diz: "Cuidado! A estrada está molhada, o asfalto está elástico e, se você acelerar demais, a estrada muda de forma. Mas, se você seguir a regra de velocidade entre as duas pistas, você chegará ao destino".

Os autores provaram matematicamente e mostraram no computador que, mesmo com toda essa complexidade, a natureza ainda segue padrões previsíveis, o que é uma ótima notícia para quem tenta prever desastres naturais ou entender o corpo humano.

Resumo técnico

1. Problema e Contexto

O artigo aborda a modelagem matemática de ondas de cisalhamento (ondas SH) não lineares que se propagam através de materiais hiperelásticos e hiper-viscoelásticos incompressíveis. O foco principal é a descrição de ondas do tipo Love, que ocorrem na interface entre dois meios com propriedades mecânicas distintas (por exemplo, a crosta terrestre e o manto superior).

Enquanto a teoria linear de ondas Love é bem estabelecida (baseada na Lei de Hooke e na equação de onda linear), ela falha em capturar efeitos importantes em grandes deformações, como os observados em sismos próximos ao epicentro, em tecidos biológicos ou em materiais poliméricos. O trabalho busca preencher essa lacuna derivando equações governantes para deslocamentos finitos, incorporando não linearidades geométricas e constitutivas, bem como efeitos de dissipação de energia (viscosidade).

2. Metodologia

A pesquisa segue uma abordagem rigorosa de mecânica do contínuo, dividida em etapas teóricas e numéricas:

• Formulação Teórica (Hiperelasticidade):

  • Os autores derivam as equações de movimento no referencial de Lagrange para materiais incompressíveis.
  • Utilizam uma função de densidade de energia de deformação ($W^H$) arbitrária, dependendo dos invariantes principais do tensor de deformação de Cauchy-Green.
  • Focam em materiais do tipo "Neo-Hookeano Generalizado" (onde a energia depende apenas do primeiro invariante $I_1$) para simplificar o sistema e garantir a existência de soluções de pressão hidrostática bem definidas.
  • Especificam o modelo para o Modelo Cúbico de Yeoh, que introduz termos polinomiais cúbicos e quárticos/quinticos na relação tensão-deformação.

• Extensão para Viscoelasticidade (Hiper-Viscoelasticidade):

  • Incorporam efeitos dissipativos através de um potencial viscoso ($W^V$) dependente da taxa de variação do tensor de deformação.
  • Derivam um termo de viscosidade não linear que inclui derivadas mistas de deslocamento (ex: $u_{xxt}$), representando a dispersão e amortecimento da onda.

• Análise Numérica:

  • Utilizam o Método das Retas (Method of Lines) para discretizar as equações diferenciais parciais (EDPs) não lineares em um domínio 2D (plano $X-Z$).
  • As derivadas espaciais são aproximadas por diferenças finitas de segunda ordem em uma malha não uniforme (com maior densidade de pontos perto da fonte da explosão e da interface).
  • O sistema de EDOs resultante é integrado no tempo usando o solucionador ode23 do MATLAB.
  • Simulações são realizadas para diferentes cenários: explosões gaussianas na camada superior, na camada inferior e em diferentes relações de velocidades de onda ($c_1$ vs $c_2$).

• Análise de Simetria:

  • Para o caso unidimensional, aplicam a análise de simetria de Lie para encontrar soluções exatas invariantes, comparando-as com os resultados numéricos.

3. Principais Contribuições

1. Derivação de Novas EDPs: Apresentação de equações de onda de cisalhamento não lineares gerais para materiais hiperelásticos incompressíveis com uma função de energia arbitrária, e sua especialização para o modelo cúbico de Yeoh.
2. Modelo Hiper-Viscoelástico: Desenvolvimento de um modelo que acopla a não linearidade hiperelástica com termos de viscosidade de ordem superior, resultando em equações que contêm termos polinomiais cúbicos e quárticos/quinticos, além de derivadas mistas de dispersão.
3. Validação Numérica em 2D: Realização de simulações completas (2+1 dimensões) de ondas Love não lineares, rastreando a evolução de frentes de onda na interface e na superfície livre.
4. Correção de Erros e Esclarecimento Teórico: O artigo inclui um adendo corrigindo uma suposição anterior sobre a dependência do segundo invariante ($I_2$), restringindo o modelo principal a materiais onde $W^H$ não depende de $I_2$ (Neo-Hookeano Generalizado) para garantir a consistência da pressão hidrostática.

4. Resultados Chave

• Condição de Existência Não Linear: As simulações demonstram que, mesmo no regime totalmente não linear, as ondas Love na interface e na superfície satisfazem a condição de existência linear clássica: $c_1 < |v| < c_2$, onde $c_1$ e $c_2$ são as velocidades de onda dos dois meios.
• Comportamento Assintótico: Observou-se que, para tempos longos, a velocidade de fase da onda não linear tende a convergir para a maior das velocidades de onda dos materiais ($\max(c_1, c_2)$), independentemente da posição inicial da fonte.
• Efeitos de Viscosidade: A introdução de termos viscosos no modelo hiper-viscoelástico resulta em uma dissipação rápida de energia e dispersão mais acentuada das ondas. Isso previne a formação de "quebras" de onda (wave breaking) que poderiam ocorrer em modelos puramente hiperelásticos sob certas condições, tornando o comportamento da onda mais previsível e estável.
• Soluções Exatas vs. Numéricas: A análise de simetria de Lie no caso 1D produziu soluções exatas que, embora ilimitadas, apresentam um perfil qualitativo muito próximo ao observado nas simulações numéricas com condições iniciais gaussianas, validando a precisão do método numérico.
• Ondas Quebrantes Numéricas: O apêndice do artigo relata a descoberta de soluções que "quebram" numericamente (formação de oscilações espúrias e descontinuidades na derivada espacial) em configurações radialmente simétricas, um fenômeno que depende da resolução da malha e da não linearidade.

5. Significado e Impacto

Este trabalho é significativo para várias áreas:

• Sismologia: Oferece um modelo mais realista para ondas Love geradas por terremotos, especialmente perto do epicentro onde as deformações são grandes e os efeitos não lineares não podem ser ignorados. A capacidade de modelar a interação entre camadas geológicas com propriedades distintas é crucial para a interpretação de sismogramas.
• Ciência dos Materiais e Biomecânica: O modelo é aplicável a materiais complexos como tecidos biológicos, polímeros e espumas, que exibem comportamento hiperelástico e viscoelástico sob grandes deformações.
• Fundamentação Teórica: A derivação rigorosa das equações de movimento para o modelo de Yeoh com viscosidade fornece uma base matemática sólida para futuras pesquisas em dinâmica de ondas não lineares em meios complexos.
• Métodos Numéricos: A aplicação bem-sucedida do Método das Retas com malhas não uniformes para este tipo de EDP não linear demonstra a viabilidade de simular fenômenos de ondas complexas em domínios estratificados.

Em resumo, o artigo avança a compreensão das ondas de cisalhamento em meios reais, transicionando da teoria linear simplificada para modelos não lineares robustos que capturam a física essencial da dissipação de energia e da dependência da amplitude na velocidade de propagação.