3D Finite Element-Based Multiphysics Simulation of a Shape Memory Alloy Hybrid Composite Module

Título: Simulação Multiphysics 3D Baseada em Elementos Finitos de um Módulo Híbrido de Compósito com Liga com Memória de Forma (SMAHC)

Autores: L. Handl, M. Kaiser, M. Duhovic e M. Gurka (Leibniz-Institut für Verbundwerkstoffe GmbH, Alemanha)

---

1. Problema e Contexto

Os atuadores baseados em Ligas com Memória de Forma (SMAs) integrados em compósitos híbridos (SMAHCs) são promissores para aplicações de adaptação de forma, como bordas de fuga de aeronaves, geradores de vórtices e robótica. Eles funcionam através de uma transformação de fase termoelástica (martensita para austenita) induzida por calor (efeito Joule), gerando deformação mecânica.

Apesar da existência de diversos modelos analíticos e de elementos finitos (FE) para prever a resposta desses materiais, a literatura atual apresenta limitações significativas:

• Muitos modelos não consideram o acoplamento eletro-termo-mecânico completo.
• Falta validação experimental abrangente para simulações 3D detalhadas.
• Modelos existentes muitas vezes ignoram fatores críticos como condutividade térmica/elétrica, temperatura ambiente e cargas mecânicas aplicadas.
• Não há um método validado que simule a geometria 3D completa de um SMAHC, incluindo a microestrutura e as condições de contorno reais.

2. Metodologia

Os autores desenvolveram uma abordagem de simulação multiphysics tridimensional (3D) utilizando o código de elementos finitos ANSYS LS-DYNA.

• Modelo de Material: Foi implementado um modelo constitutivo micromecânico (baseado em Kelly et al. [25]) que descreve a cinética de crescimento e saturação das fases de Martensita e Austenita. O modelo minimiza a energia livre de Helmholtz, considerando variáveis de estado internas como a fração volumétrica de martensita ($\lambda$) e a tensão de martensita.
• Acoplamento Multiphysics: O modelo integra três solucionadores:
  1. Mecânico: Para deformação e tensão.
  2. Térmico: Para transferência de calor e convecção.
  3. Eletromagnético (EM): Um solver simplificado que calcula o aquecimento Joule (resistência ôhmica) sem efeitos de indução, acoplado ao solver térmico como fonte de calor.
• Procedimento de Inicialização Física: Um desafio crítico abordado foi a definição do estado inicial da martensita. Como o software inicializa o material puramente em estado austenítico, os autores introduziram uma etapa de pré-estiramento na simulação. Neste passo, o fio SMA é esticado mecanicamente a 23°C para induzir o "desembaralhamento" (detwinning) parcial da microestrutura de martensita, replicando o estado físico real do atuador antes da ativação térmica.
• Geometria e Malha: O modelo representa um atuador comercial (CompActive GmbH, Tipo A3950) com diferentes espessuras de substrato e intercamadas. Foram utilizados elementos hexaédricos sólidos (32.360 elementos) para garantir precisão numérica e estabilidade em grandes deformações.
• Validação: Os resultados foram comparados com:
  1. Dados experimentais reais (deflexão, temperatura, histerese).
  2. Um modelo de esquema escalonado (SSM) totalmente acoplado e validado anteriormente.

3. Contribuições Chave

• Implementação 3D Completa: Primeira simulação 3D detalhada de um SMAHC que integra acoplamento eletro-termo-mecânico completo em um ambiente de engenharia assistida por computador (CAE) padrão (LS-DYNA).
• Protocolo de Inicialização: Desenvolvimento de um método de pré-processamento para estabelecer um estado de martensita desembaralhada fisicamente motivado, superando a limitação de inicialização padrão do software.
• Validação Multiescala: Comparação rigorosa entre simulação 3D, modelo 1D (SSM) e dados experimentais, cobrindo diferentes correntes elétricas, temperaturas ambiente e cargas mecânicas.
• Visualização de Campo: Capacidade de visualizar a distribuição de temperatura e deformação da malha em tempo real, permitindo a análise de gradientes térmicos ao longo do fio.

4. Resultados

• Aquecimento e Temperatura: A simulação reproduziu qualitativamente bem a taxa de aquecimento e a temperatura máxima (403 K). Observou-se uma distribuição de temperatura inhomogênea (o centro aquece mais rápido que as extremidades), consistente com os experimentos.
• Resfriamento: Houve uma discrepância significativa na fase de resfriamento. A simulação resfria mais rápido do que o experimento real, especialmente em temperaturas ambiente mais altas. O modelo não capturou totalmente a inércia térmica observada experimentalmente durante o resfriamento.
• Deflexão e Histerese:
  • As curvas de deflexão mostraram boa concordância qualitativa, reproduzindo a histerese característica dos SMAs.
  • Os valores de deflexão máxima estão na ordem de grandeza correta, embora a simulação tenda a subestimar ligeiramente a deflexão máxima em comparação com os dados experimentais.
  • Em cargas mais altas, o atuador simulado mostrou-se mais rígido que o físico, resultando em desvios maiores.
• Comparação com SSM: Em geral, o modelo de esquema escalonado (SSM) apresentou melhor concordância quantitativa com os dados experimentais do que a simulação FE 3D, especialmente nas fases de resfriamento e em cargas elevadas. No entanto, a abordagem 3D capturou corretamente as tendências gerais e os fenômenos físicos.

5. Significância e Conclusão

O estudo demonstra que é possível simular com sucesso a resposta dinâmica de atuadores SMAHC complexos em 3D, integrando efeitos elétricos, térmicos e mecânicos. Embora existam desvios quantitativos (principalmente relacionados à modelagem da rigidez das camadas intermédias e à taxa de resfriamento), a abordagem oferece uma ferramenta de engenharia robusta para prever o comportamento de sistemas mais complexos.

Conclusões Principais:

1. O aquecimento Joule e a temperatura máxima são modelados com boa concordância qualitativa.
2. O resfriamento convectivo apresenta discrepâncias significativas, indicando a necessidade de refinamento nos modelos de transferência de calor ou propriedades do material.
3. A deflexão é modelada de forma razoável, com tendências consistentes entre simulação e experimento.
4. A abordagem permite a exploração futura de sistemas 3D mais complexos, com planos para melhorar o modelo da camada intermédia de elastômero e realizar análises de sensibilidade.

Este trabalho representa um passo importante na transição de modelos simplificados para simulações multiphysics 3D completas para o projeto de atuadores inteligentes integrados.