Resumo Técnico: SoliDualSPHysics

Definição do Problema

A Hidrodinâmica de Partículas Suavizadas (SPH) estabeleceu-se como um poderoso método livre de malha para a dinâmica de fluidos e mecânica de sólidos com grandes deformações devido à sua natureza lagrangiana. No entanto, o cenário atual de frameworks SPH de alto desempenho e código aberto está fortemente inclinado para a dinâmica de fluidos computacional (CFD) e interação fluido-estrutura (FSI) envolvendo corpos rígidos ou elasticidade de pequenas deformações. Existe uma lacuna significativa em frameworks unificados capazes de lidar com mecânica de sólidos deformáveis autônomos com comportamentos constitutivos avançados, especificamente plasticidade de deformação finita e fratura frágil fisicamente consistente.

As formulações de fratura SPH existentes frequentemente dependem de modelos de dano empíricos ad hoc, abordagens de partículas de trinca ou modelos de zona coesiva, que sofrem com dependência de malha, resultados não convergentes ou a necessidade de rastreamento explícito de trincas e enriquecimento local. Além disso, integrar essas capacidades avançadas de mecânica dos sólidos em arquiteturas escaláveis CPU/GPU adequadas para simulações de grande escala continua sendo um desafio. A necessidade específica identificada é um formulário SPH unificado que suporte hiperelasticidade, plasticidade de deformação finita e fratura frágil sem rastreamento explícito de trincas, mantendo a eficiência computacional através da paralelização.

Metodologia

Os autores apresentam o SoliDualSPHysics, uma extensão de código aberto e acelerada por GPU do framework DualSPHysics. A metodologia é construída sobre uma formulação Total Lagrangian SPH (TLSPH), que permite a aplicação direta de cargas externas e condições de contorno na configuração não deformada.

Estrutura Teórica

1. Balanço de Momento: O solver emprega uma discretização TLSPH de forma forte para a equação de balanço de momento. Utiliza uma correção de gradiente de kernel de primeira ordem para restaurar a consistência e inclui um termo de viscosidade artificial para mitigar instabilidades numéricas e modos de energia zero espúrios.
2. Modelagem Constitutiva:
   • Hiperelasticidade: O framework suporta os modelos de St. Venant–Kirchhoff e neo-Hookean compressível. Estes são acoplados a uma formulação de campo de fase para simular a fratura frágil.
   • Plasticidade de Deformação Finita: Uma extensão para plasticidade $J_2$ independente de taxa com endurecimento isotrópico é implementada. Esta utiliza uma decomposição multiplicativa do gradiente de deformação ($F = F^e F^p$) e um algoritmo de retorno de mapa exponencial para garantir a objetividade de referencial e consistência termodinâmica sob grandes rotações.
   • Nota de Acoplamento: Na implementação atual, a plasticidade e a fratura por campo de fase não são acopladas; o modelo de plasticidade descreve a resposta elastoplástica sem fratura, enquanto o modelo de campo de fase lida com a fratura frágil em materiais hiperelásticos.
3. Fratura Frágil (Campo de Fase): A fratura é modelada usando uma formulação de campo de fase hiperbólica. Diferente das formulações elípticas que exigem a resolução de sistemas lineares globais (incompatíveis com integração de tempo explícita), esta abordagem utiliza uma equação de evolução hiperbólica para a variável de dano $s$. Isso permite que a iniciação, propagação, ramificação e coalescência de trincas surjam naturalmente sem rastreamento ou enriquecimento explícito. O modelo impõe irreversibilidade via um funcional de histórico e restringe a propagação de trincas a regiões de tração através da decomposição da energia de deformação elástica.
4. Integração Temporal: O sistema é avançado usando esquemas de integração de tempo explícitos (Verlet ou Simplético/Leapfrog), com passos de tempo restritos por condições do tipo CFL baseadas em propagação de ondas e aceleração.

Arquitetura de Software e Implementação

• Plataforma: O código é escrito em C++ e CUDA, aproveitando a arquitetura paralela nativa de CPU (OpenMP) e GPU (CUDA) do DualSPHysics.
• Entrada/Saída: A configuração é gerenciada via arquivos de entrada XML. A extensão introduz novos blocos (<deformstrucs>, <mathexpressions>) para definir corpos deformáveis, propriedades de materiais, condições de contorno e expressões definidas pelo usuário.
• Expressões Definidas pelo Usuário: Uma característica fundamental é a inclusão de um analisador de expressões matemáticas que permite aos usuários prescrever quantidades dependentes de tempo e espaço (ex: campos de velocidade, cargas de tração) usando variáveis integradas (posição, tempo, deslocamento) e operadores matemáticos/lógicos padrão.
• Condições de Contorno: Condições de contorno flexíveis de velocidade e tração podem ser aplicadas globalmente ou a grupos de superfície específicos definidos por geometrias auxiliares.

Principais Contribuições

1. Framework Unificado: O SoliDualSPHysics fornece um ambiente único e de código aberto para simular hiperelasticidade, plasticidade de deformação finita e fratura frágil dentro do paradigma SPH.
2. Fratura por Campo de Fase em SPH: Implementa uma formulação de campo de fase hiperbólica para fratura frágil que evita o custo computacional de solvers elípticos, permitindo a simulação eficiente de ramificação e coalescência de trincas sem critérios ad hoc.
3. Plasticidade de Deformação Finita: A integração de um modelo de plasticidade $J_2$ multiplicativa com mapeamento de retorno de exponencial estende a aplicabilidade do código para problemas plásticos de grande deformação.
4. Escalabilidade: Ao utilizar a infraestrutura CPU/GPU do DualSPHysics, o software alcança aceleração substancial para simulações de grande escala envolvendo milhões de partículas.
5. Flexibilidade: A inclusão de um analisador de expressões definido pelo usuário e definições flexíveis de condições de contorno aumenta a adaptabilidade do código para cenários de carregamento complexos e aplicações de pesquisa.
6. Disponibilidade de Código Aberto: Todo o código fonte, documentação e casos de teste foram disponibilizados gratuitamente via um repositório público no GitHub para garantir a reprodutibilidade.

Resultados e Validação

O software foi verificado e validado contra uma série de problemas numéricos de referência e dados experimentais:

• Elasto-Dinâmica: Oscilações livres de vigas e placas cantilever (2D e 3D) demonstraram excelente concordância com soluções analíticas e resultados de referência do Método de Elementos Finitos (FEM). As simulações mostraram troca de energia quase conservativa e dissipação numérica mínima.
• Grande Deformação: Uma viga cantilever 3D sujeita a uma grande tração superficial distribuída confirmou a capacidade de lidar com não linearidade geométrica e condições de contorno complexas sem oscilações espúrias.
• Dinâmica de Torção: Um caso de coluna de torção 3D validou a capacidade do framework de resolver a propagação de ondas torcionais e recuperação elástica, com resultados correspondentes a estudos de referência.
• Fratura: Um teste de referência de ramificação de trinca dinâmica (placa pré-entalhada sob tensão) demonstrou a capacidade do modelo de capturar padrões de fratura complexos, incluindo ramificação, consistente com a literatura estabelecida de campo de fase.
• Desempenho: A implementação em GPU mostrou entregar aceleração significativa em comparação com execuções apenas em CPU, suportando simulações escaláveis.

Significância e Escopo

O artigo posiciona o SoliDualSPHysics como um passo significativo para preencher a lacuna entre solvers de fluidos SPH de alto desempenho e capacidades avançadas de mecânica dos sólidos. Sua principal significância reside em fornecer uma ferramenta unificada, de código aberto e acelerada por GPU que elimina a necessidade de rastreamento ad hoc de trincas em simulações de fratura frágil e suporta plasticidade de deformação finita dentro de um framework livre de malha.

Os autores explicitamente observam as limitações da versão atual:

• É restrita a resposta de material isotrópica (hiperelasticidade isotrópica, plasticidade $J_2$ isotrópica).
• Plasticidade e fratura não são acopladas; mecanismos de fratura dúctil não estão incluídos.
• Modelos constitutivos anisotrópicos, materiais compostos e operadores dedicados de controle de hourglass ainda não foram implementados.

Apesar dessas limitações, este trabalho estabelece uma base robusta para desenvolvimentos futuros em mecânica dos sólidos e modelagem de fratura dentro da comunidade SPH, oferecendo uma plataforma reprodutível e extensível para pesquisadores.