Model-Informed Joint Material-Structural Optimization of Hard-Magnetic Soft Materials

Resumo Técnico: Otimização Conjunta Modelo-Informada de Material-Estrutura de Materiais Magnéticos de Dureza Suave

Definição do Problema
Materiais magnéticos de dureza suave (hMSMs) são compósitos consistindo em partículas magnéticas de dureza elevada incorporadas em uma matriz elastomérica complacente. Eles permitem deformação livre de contato, impulsionada por campo, adequada para robótica suave e estruturas adaptativas. No entanto, o design preditivo de hMSMs enfrenta dois desafios interconectados:

1. Modelagem Constitutiva: A previsão precisa requer relações estrutura-propriedade eficazes que descrevam como a fração volumétrica de partículas magnéticas modifica os módulos elásticos efetivos. A literatura existente oferece vários modelos (inclusão rígida, autoconsistente, cinemática restrita) que podem prever valores de rigidez significativamente diferentes para a mesma concentração de partículas. Além disso, a escolha da função de densidade de energia de deformação pode influenciar a robustez numérica e a deformação prevista.
2. Otimização Conjunta: O design otimizado requer o controle simultâneo da densidade estrutural, da distribuição de partículas magnéticas e da direção da magnetização remanente. As abordagens de otimização topológica prévias para hMSMs frequentemente tratam a presença de material magnético como uma escolha discreta (magnetizado vs. não magnetizado) em vez de um campo de fração volumétrica contínuo que modula tanto a rigidez quanto a força de atuação. Adicionalmente, acoplar esses campos fisicamente é difícil porque as variáveis magnéticas só fazem sentido onde existe material estrutural.

Metodologia
Os autores desenvolvem um framework unificado que aborda tanto a seleção do modelo constitutivo quanto a otimização conjunta.

• Framework Constitutivo Unificado: Os autores reestruturam relações clássicas de inclusão rígida (Guth, Mooney, Kerner), a teoria de autoconsistência de Hill (especializada para inclusões rígidas e matrizes quase incompressíveis) e modelos de cinemática restrita (Lopez-Pamies, Alshammari) em uma forma comum de módulo de cisalhamento efetivo, $G(\phi)$, onde $\phi$ é a fração volumétrica de partículas. Isso resulta em sete relações de módulo de cisalhamento candidatas. Estas são combinadas com três funções de densidade de energia de deformação (duas formas neo-Hookeanas compressíveis, NH1 e NH2, e o modelo de St. Venant–Kirchhoff, SVK) para criar 21 modelos constitutivos candidatos.
• Análise de Sensibilidade do Modelo: Benchmarks numéricos (vigas em balanço, rodas, garras e vigas 3D) são usados para avaliar a sensibilidade da deformação prevista em relação à escolha da função de densidade de energia de deformação versus a escolha de $G(\phi)$.
• Seleção de Modelo Constitutivo: Dados experimentais de tensão-deformação uniaxial de Garcia-Gonzalez et al. (2021) são usados para selecionar a relação $G(\phi)$ mais representativa. Duas métricas são empregadas: (1) erro relativo contra valores de módulo de cisalhamento extraídos, e (2) erro de raiz quadrada média relativa (rRMSE) contra curvas completas de tensão-deformação reais.
• Formulação de Otimização Conjunta Material-Estrutura: Um framework de otimização topológica é desenvolvido onde três campos de design são parametrizados: densidade estrutural ($\rho$), fração volumétrica de partículas magnéticas ($\phi$) e direção de magnetização remanente ($\theta_{Br}$).
  • Mecanismo de Acoplamento: Para garantir a consistência física, a energia magnética e os módulos efetivos são interpolados de modo que as contribuições magnéticas desapareçam em regiões de vazio ($\rho \approx 0$). A fração volumétrica de partículas magnéticas $\phi$ atua como um campo contínuo que escala tanto a energia de interação magnética quanto o módulo de cisalhamento efetivo via a relação $G(\phi)$ selecionada.
  • Implementação: O framework utiliza a biblioteca FEniTop e FEniCSx para implementação de elementos finitos, empregando uma penalização do tipo SIMP para a densidade, filtragem de Helmholtz para todos os campos, e um cálculo de gradiente baseado em adjunção para o otimizador Method of Moving Asymptotes (MMA).

Principais Resultados

• Sensibilidade do Modelo Constitutivo: A função de densidade de energia de deformação (NH1, NH2, SVK) tem um efeito relativamente pequeno na deformação prevista para os problemas de atuação considerados. Em contraste, a escolha da relação de módulo de cisalhamento efetivo $G(\phi)$ afeta significativamente a deformação, particularmente quando o material magnético sobrepõe-se a regiões de alta deformação. Quando o material magnético é localizado longe das principais regiões de flexão (ex: vigas magnetizadas na extremidade), a dependência do modelo é negligenciável.
• Seleção de Modelo: Com base nos dados experimentais, a relação de Mooney para $G(\phi)$ demonstrou o melhor acordo geral tanto com os valores de módulo de cisalhamento extraídos quanto com as respostas completas de tensão-deformação. A função de densidade de energia de deformação NH2 foi selecionada por sua robustez numérica.
• Desempenho de Otimização: O framework de otimização conjunta produziu designs não intuitivos para três casos distintos:
  1. Atuador Rotacional: A otimização de $\phi$ e $\theta_{Br}$ em uma geometria de roda fixa aumentou o deslocamento tangencial em 59,7% comparado a um design uniforme intuitivo.
  2. Atuador Translacional: A otimização de $\phi$ e $\theta_{Br}$ em uma geometria tipo tesoura fixa alcançou o movimento horizontal alvo enquanto suprimiu o deslocamento parasitário vertical.
  3. Viga Restauradora: A otimização simultânea de $\rho$, $\phi$ e $\theta_{Br}$ criou uma viga capaz de restaurar sua forma sob cargas mecânicas opostas usando atuação magnética. O design confinou naturalmente o material magnético às regiões estruturalmente ocupadas.

Significância e Alegações
O artigo afirma fornecer uma rota sistemática e informada por modelos para o design de hMSMs. Suas principais contribuições são:

1. Um framework de módulo de cisalhamento efetivo unificado que permite a comparação direta de diversas relações estrutura-propriedade.
2. Evidência de que as relações de módulo de cisalhamento efetivo são críticas para a precisão preditiva quando os campos magnético e mecânico se sobrepõem, enquanto as formas de densidade de energia de deformação são menos sensíveis para os casos testados.
3. Uma formulação de otimização conjunta material-estrutura que trata a fração volumétrica de partículas magnéticas como uma variável de design contínua acoplada à densidade estrutural e à direção de magnetização, resolvendo o problema do acoplamento físico em regiões de baixa densidade através da interpolação de energia.
4. Demonstração de que a otimização simultânea da composição do material e da topologia estrutural produz um desempenho superior ao design magnético de geometria fixa.

Os autores observam que, embora o framework seja demonstrado em hMSMs, a estrutura computacional é aplicável a outros elastômeros responsivos a campos (ex: dielétricos, magnetorreológicos ou compósitos de memória de forma) substituindo os termos de energia relevantes. A implementação está disponível no repositório de código aberto CEADpx/top_optim para apoiar a reprodutibilidade e extensão.