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

Resumen Técnico: Optimización Conjunta de Material y Estructura Informada por Modelos de Materiales Blandos con Magnetismo Duro

Planteamiento del Problema
Los materiales blandos con magnetismo duro (hMSM, por sus siglas en inglés) son compuestos consistentes en partículas de magnetismo duro embebidas en una matriz elastomérica flexible. Estos permiten la deformación impulsada por campos sin contacto, adecuada para la robótica blanda y estructuras adaptativas. Sin embargo, el diseño predictivo de los hMSM enfrenta dos desafíos interconectados:

1. Modelado Constitutivo: La predicción precisa requiere relaciones efectivas estructura-propiedad que describan cómo la fracción de volumen de las partículas magnéticas modifica los módulos elásticos efectivos. La literatura existente ofrece diversos modelos (inclusión rígida, autoconsistente, de cinemática restringida) que pueden predecir valores de rigidez significativamente diferentes para la misma concentración de partículas. Además, la elección de la función de densidad de energía de deformación puede influir en la robustez numérica y en la deformación predicha.
2. Optimización Conjunta: El diseño óptimo requiere el control simultáneo de la densidad estructural, la distribución de las partículas magnéticas y la dirección de la magnetización remanente. Los enfoques previos de optimización topológica para hMSM suelen tratar la presencia de material magnético como una elección discreta (magnetizado vs. no magnetizado) en lugar de un campo de fracción de volumen continuo que modula tanto la rigidez como la fuerza de actuación. Además, acoplar estos campos físicamente es difícil porque las variables magnéticas solo tienen sentido donde existe el material estructural.

Metodología
Los autores desarrollan un marco unificado que aborda tanto la selección del modelo constitutivo como la optimización conjunta.

• Marco Constitutivo Unificado: Los autores reformulan las relaciones clásicas de inclusión rígida (Guth, Mooney, Kerner), la teoría de auticonsistencia de Hill (especializada para inclusiones rígidas y matrices casi incompresibles) y los modelos de cinemática restringida (Lopez-Pamies, Alshammari) en una forma común de módulo de corte efectivo, $G(\phi)$, donde $\phi$ es la fracción de volumen de partículas. Esto genera siete relaciones candidatas de módulo de corte. Estas se combinan con tres funciones de densidad de energía de deformación (dos formas neo-Hookeanas compresibles, NH1 y NH2, y el modelo de St. Venant–Kirchhoff, SVK) para crear 21 modelos constitutivos candidatos.
• Análisis de Sensibilidad del Modelo: Se utilizan pruebas numéricas de referencia (vigas en voladizo, ruedas, pinzas y vigas 3D) para evaluar la sensibilidad de la deformación predicha a la elección de la función de densidad de energía de deformación frente a la elección de $G(\phi)$.
• Selección del Modelo Constitutivo: Se utilizan datos experimentales de esfuerzo-deformación uniaxial de Garcia-Gonzalez et al. (2021) para seleccionar la relación $G(\phi)$ más representativa. Se emplean dos métricas: (1) el error relativo contra los valores de módulo de corte extraídos, y (2) el error cuadrático medio relativo (rRMSE) contra las curvas completas de esfuerzo-deformación reales.
• Formulación de la Optimización Conjunta de Material y Estructura: Se desarrolla un marco de optimización topológica donde tres campos de diseño están parametrizados: la densidad estructural ($\rho$), la fracción de volumen de partículas magnéticas ($\phi$) y la dirección de magnetización remanente ($\theta_{Br}$).
  • Mecanismo de Acoplamiento: Para asegurar la consistencia física, la energía magnética y los módulos efectivos se interpolan de tal manera que las contribuciones magnéticas desaparecen en las regiones de vacío ($\rho \approx 0$). La fracción de volumen de partículas magnéticas $\phi$ actúa como un campo continuo que escala tanto la energía de interacción magnética como el módulo de corte efectivo mediante la relación $G(\phi)$ seleccionada.
  • Implementación: El marco utiliza la biblioteca FEniTop y FEniCSx para la implementación de elementos finitos, empleando una penalización de tipo SIMP para la densidad, un filtrado de Helmholtz para todos los campos y un cálculo de gradiente basado en el adjunto para el optimizador Method of Moving Asymptotes (MMA).

Resultados Clave

• Sensibilidad Constitutiva: La función de densidad de energía de deformación (NH1, NH2, SVK) tiene un efecto relativamente pequeño en la deformación predicha para los problemas de actuación considerados. En contraste, la elección de la relación del módulo de corte efectivo $G(\phi)$ afecta significativamente la deformación, particularmente cuando el material magnético se superpone con regiones de alta deformación. Cuando el material magnético se localiza lejos de las regiones de flexión primaria (por ejemplo, vigas magnetizadas en el extremo), la dependencia del modelo es despreciable.
• Selección del Modelo: Basándose en los datos experimentales, la relación de Mooney para $G(\phi)$ demostró la mejor concordancia general tanto con los valores de módulo de corte extraídos como con las respuestas completas de esfuerzo-deformación. Se seleccionó la función de densidad de energía NH2 por su robustez numérica.
• Desempeño de la Optimización: El marco de optimización conjunta produjo con éxito diseños no intuitivos para tres casos distintos:
  1. Actuador Rotacional: La optimización de $\phi$ y $\theta_{Br}$ en una geometría de rueda fija aumentó el desplazamiento tangencial en un 59.7% en comparación con un diseño uniforme intuitivo.
  2. Actuador Traslacional: La optimización de $\phi$ y $\theta_{Br}$ en una geometría tipo tijera fija logró un movimiento horizontal objetivo mientras se suprimía el desplazamiento vertical parásito.
  3. Viga Restauradora: La optimización simultánea de $\rho$, $\phi$ y $\theta_{Br}$ creó una viga capaz de restaurar su forma bajo cargas mecánicas opuestas utilizando actuación magnética. El diseño confinó naturalmente el material magnético a las regiones estructuralmente ocupadas.

Significancia y Reivindicaciones
El artículo afirma proporcionar una ruta sistemática e informada por modelos para diseñar hMSM. Sus principales contribuciones son:

1. Un marco de módulo de corte efectivo unificado que permite la comparación directa de diversas relaciones estructura-propiedad.
2. Evidencia de que las relaciones del módulo de corte efectivo son críticas para la precisión predictiva cuando los campos magnéticos y mecánicos se superponen, mientras que las formas de la densidad de energía de deformación son menos sensibles para los casos probados.
3. Una formulación de optimización conjunta de material-estructura que trata la fracción de volumen de partículas magnéticas como una variable de diseño continua acoplada con la densidad estructural y la dirección de magnetización, resolviendo el problema del acoplamiento físico en las regiones de baja densidad mediante la interpolación de energía.
4. Demostración de que la optimización simultánea de la composición del material y la topología estructural produce un rendimiento superior al del diseño magnético de geometría fija.

Los autores señalan que, aunque el marco se demuestra en hMSM, la estructura computacional es aplicable a otros elastómeros sensibles a campos (por ejemplo, compuestos dieléctricos, magnetorreológicos o de memoria de forma) sustituyendo los términos de energía correspondientes. La implementación está disponible en el repositorio de código abierto CEADpx/top_optim para apoyar la reproducibilidad y la extensión.