Homogenization of rod-like metamaterials as a special Cosserat rod

Este artículo presenta un esquema de homogeneización variacional basado en la teoría de varillas de Cosserat especial geométricamente exacta para derivar la respuesta constitutiva no lineal y la rigidez de metamateriales tipo varilla ensamblados periódicamente, validado mediante ejemplos numéricos que van desde retículas simples hasta estructuras complejas auxéticas y de músculo artificial.

Lo esencial

Imagina que tienes un tubo largo y flexible hecho no de un material sólido único, sino de un patrón complejo y repetitivo de pequeños palos conectados entre sí, como una escalera microscópica o una cerca de malla metálica enrollada en forma de cilindro. Esto es lo que los autores denominan un "metamaterial de tipo varilla".

El problema que abordan es el siguiente: si quieres saber cómo se dobla, estira o tuerce todo este tubo largo, no puedes simplemente observar un solo palito. Tienes que examinar cómo interactúa toda la red de miles de palos. Simular cada palito individual para un tubo largo es como intentar contar cada grano de arena en una playa para entender cómo se mueve la playa con el viento: requiere demasiada potencia de computación.

Los autores proponen un atajo ingenioso, una "receta" para predecir cómo se comporta todo el tubo estudiando solo una pieza diminuta y representativa de él. Así es como lo hacen, explicado con analogías sencillas:

1. El "Zoom Mágico" (Homogeneización)

Piensa en el metamaterial como un patrón de papel tapiz gigante y repetitivo. En lugar de analizar toda la pared, solo miras un solo cuadrado del papel tapiz (llamado RVE, o Volumen Elemento Representativo).

El truco de los autores consiste en asumir que si estiras o tuerces todo el tubo largo, ese pequeño trozo cuadrado también se estira o tuerce, pero de una manera muy específica y en espiral. Lo llaman una deformación "helical". Imagina tomar un resorte y tirarlo; las espiras no solo se alargan, sino que también giran ligeramente. Los autores se dieron cuenta de que, al obligar a ese pequeño trozo cuadrado a imitar exactamente ese movimiento en espiral, pueden determinar cómo reaccionaría todo el tubo largo sin simularlo por completo.

2. Las Varillas "Perfectamente Flexibles"

La mayoría de los modelos informáticos tratan las varillas como rígidas e inmutables, como una regla de acero. Pero en la vida real, especialmente con estos metamateriales diminutos, las varillas pueden doblarse, estirarse y cortarse (deslizarse lateralmente) todo a la vez, incluso cuando la deformación es enorme.

Los autores utilizan un modelo matemático especial llamado "Varilla Cosserat Especial".

• Analogía: Imagina un trozo de espagueti cocido. Puede doblarse, puede estirarse un poco y puede torcerse. Ahora imagina que ese espagueti está hecho de un material que puede hacer todas estas cosas de manera perfecta y precisa, incluso si lo doblas en un círculo o lo estiras hasta el doble de su longitud. Eso es lo que hace su modelo. No solo aproxima; captura la geometría exacta de la flexión y la torsión.

3. Las Reglas de la "Pista de Baile" (Condiciones de Contorno)

Para hacer que el pequeño trozo cuadrado se comporte como si fuera parte de un tubo gigante y repetitivo, los autores tuvieron que inventar un conjunto de reglas sobre cómo los bordes de ese cuadrado se comunican entre sí.

• El Problema: Si cortas un trozo de una escalera de caracol, el borde superior no coincide perfectamente con el borde inferior.
• La Solución: Crearon una "condición de contorno helical". Imagina que el lado izquierdo de tu pequeño trozo cuadrado se da la mano con el lado derecho, pero el lado derecho está ligeramente girado y desplazado, al igual que los peldaños de una escalera de caracol.
• La Innovación: Los métodos anteriores solo podían manejar movimientos pequeños y suaves. La nueva regla de los autores funciona incluso si el tubo se tuerce en forma de pretzel o se estira hasta quedar tan delgado como un hilo. Es "geométricamente exacta", lo que significa que nunca pierde precisión, sin importar cuán salvaje se vuelva la forma.

4. Las "Articulaciones" y el "Pegamento"

Dentro de ese pequeño trozo cuadrado, las varillas están conectadas en articulaciones.

• Articulaciones Rígidas: Algunas articulaciones son como un pegamento superfuerte; las varillas no pueden moverse una respecto a la otra en el punto de conexión.
• Las Matemáticas: Los autores establecieron un sistema donde la computadora calcula la mejor posición de cada varilla en ese pequeño cuadrado, asegurando que las articulaciones permanezcan conectadas y que se sigan las reglas de la "escalera de caracol", mientras se utiliza la mínima cantidad de energía posible.

5. Lo Que Encontraron (Los Resultados)

Una vez que resolvieron las matemáticas para la pieza diminuta, pudieron predecir cómo actuaría todo el tubo. Lo probaron con diferentes formas:

• Las Formas de Cruz y Cuadrado: Comenzaron con formas simples (como una cruz o un cuadrado hecho de varillas) para demostrar que sus matemáticas funcionaban. Descubrieron que si las varillas diminutas son gruesas y cortas, importa mucho si pueden estirarse o cortarse. Si son muy delgadas y largas, las matemáticas antiguas y más simples funcionan bien.
• Las Varillas Helicales (Resortes): Observaron un cuadrado hecho de varillas que ya están curvadas como resortes (hélices).
  • El Estiramiento en "J": Cuando estiraron este material, al principio era blando (como un resorte desenrollándose) pero se volvió muy rígido a medida que se enderezaba. Esto crea una curva en forma de "J". Esto es exactamente cómo se comportan los tejidos biológicos (como los músculos), razón por la cual los autores mencionan que podrían usarse para músculos artificiales.
  • El Ablandamiento al Doblar: Cuando lo doblaron, el material se volvió más blando cuanto más lo doblaban. Esto sucedió porque la varilla-resorte conectada comenzó a torcerse fuera del plano, actuando como una bisagra.
• El Tubo Auxético: Modelaron un tubo hueco que se ensancha cuando lo tiras (como un panal).
  • Mostraron que al cambiar el ángulo de las varillas, puedes ajustar el tubo para que sea muy flexible de lado a lado (bueno para doblarse) pero muy rígido contra ser aplastado (bueno para mantener abiertos los vasos sanguíneos).
  • Señalaron que estas estructuras pueden ajustarse para evitar el "acortamiento" (hacerse más corto al expandirse), que es un problema común en stents cardiovasculares (tubos de malla utilizados para mantener abiertas las arterias).

Resumen

Los autores construyeron un "traductor universal" para metamateriales. Crearon un método que toma una red compleja y tridimensional de varillas diminutas y la traduce en una descripción matemática simple y suave de una sola varilla. Esto permite a los ingenieros diseñar materiales complejos y flexibles para cosas como brazos robóticos, músculos artificiales y stents médicos ajustando los patrones internos diminutos, sabiendo exactamente cómo se doblará y estirará el producto final, sin necesidad de ejecutar una simulación de supercomputadora para cada cambio de diseño.

Resumen técnico

Resumen Técnico: Homogeneización de Metamateriales Tipo Varilla como una Varilla Cosserat Especial

Enunciado del Problema
Los metamateriales tipo varilla son estructuras formadas por el ensamblaje periódico de unidades microestructurales (típicamente redes de varillas) en una sola dirección, lo que resulta en una geometría macroscópica significativamente más larga en una dimensión que en las otras dos. Estas estructuras se utilizan en aplicaciones que van desde stents cardiovasculares y músculos artificiales hasta actuadores robóticos. Si bien la simulación numérica directa a escala completa de tales estructuras heterogéneas es computacionalmente prohibitiva, las técnicas de homogeneización estándar a menudo enfrentan limitaciones. Específicamente, los esquemas de homogeneización de primer orden son insuficientes para capturar grandes deformaciones que involucran gradientes de deformación significativos (como flexión y torsión) típicos de las estructuras tipo varilla. Además, los enfoques estándar de homogeneización tridimensional no logran tener en cuenta los efectos de tamaño cuando la separación de escalas de longitud no se cumple en las direcciones transversales. Existe la necesidad de un marco que pueda modelar con precisión la respuesta constitutiva no lineal de estos metamateriales bajo deformaciones grandes arbitrarias, teniendo en cuenta el cizallamiento, estiramiento, flexión y torsión finitos tanto a escala micro como macro.

Metodología
Los autores proponen un esquema de homogeneización computacional que modela los metamateriales tipo varilla como una "varilla Cosserat especial" a escala macroscópica, mientras que la unidad microestructural (Volumen Elemento Representativo o RVE) se modela como una red de varillas Cosserat geométricamente exactas.

1. Cinemática Geométricamente Exacta: Tanto la varilla efectiva macroscópica como las varillas constituyentes microscópicas se modelan utilizando la teoría de la varilla Cosserat especial geométricamente exacta. Esta formulación tiene en cuenta la flexión, torsión, cizallamiento y estiramiento axial finitos, evitando las suposiciones de pequeñas deformaciones de las teorías clásicas de vigas.
2. Extensión de la Regla de Cauchy-Born Helical (HCB): El método extiende la regla de Cauchy-Born helical. Asume que el metamaterial está sometido a un campo de deformación uniforme a lo largo de su longitud de arco a escala macroscópica. Esta suposición permite reducir el problema estructural completo al de un solo RVE.
3. Condiciones de Contorno Periódicas Helicoidales: Una innovación clave es la derivación de un mapeo no lineal geométricamente exacto para aplicar condiciones de contorno periódicas helicoidales al RVE. Este mapeo relaciona los grados de libertad traslacionales y rotacionales del límite derecho del RVE con su límite izquierdo basándose en los parámetros de deformación macroscópica ($v_0, k_0$). A diferencia de los mapas de transición lineales anteriores, esta formulación permanece válida para deformaciones macroscópicas arbitrariamente grandes.
4. Minimización de Energía Restringida: El problema del RVE se formula como una minimización restringida de la energía elástica. Las restricciones incluyen:
   • Restricciones de Articulación Interna: Imponer una conexión rígida entre las varillas dentro del RVE utilizando un enfoque maestro-esclavo.
   • Restricciones de Contorno Helicoidales: Imponer la periodicidad derivada de la regla HCB.
   • Restricciones Globales: Restringir la traslación y rotación de cuerpo rígido del RVE para asegurar una solución única y mantener fijos los parámetros de deformación macroscópica durante la minimización.
5. Discretización y Solución: La energía del RVE se discretiza utilizando "elementos helicoidales", donde se asume que cada segmento experimenta un campo de deformación uniforme, haciendo que el segmento en sí sea helicoidal. La forma débil no lineal resultante se resuelve utilizando un código interno en la plataforma deal.II.
6. Derivación de Propiedades Macroscópicas: Se derivan expresiones en forma cerrada para la fuerza de contacto interna, el momento y el tensor de rigidez de la varilla homogeneizada mediante la diferenciación de la densidad de energía del RVE minimizada con respecto a los parámetros de deformación macroscópica.

Contribuciones Clave
El artículo identifica tres contribuciones novedosas principales:

1. Exactitud Geométrica Dual-Escala: Tanto las varillas a escala micro como macro se modelan como varillas Cosserat especiales geométricamente exactas, capturando el cizallamiento y estiramiento finitos además de la flexión y torsión en ambas escalas.
2. Mapa Macro-a-Micro Geométricamente Exacto: Se propone un mapeo completamente no lineal para aplicar condiciones de contorno periódicas helicoidales a los grados de libertad traslacionales y rotacionales. Este mapa es válido incluso a grandes deformaciones macroscópicas, abordando una limitación de los mapas de transición lineales anteriores.
3. Comportamiento Constitutivo Microscópico Arbitrario: El marco permite un comportamiento constitutivo arbitrario de las varillas constituyentes dentro del RVE, yendo más allá de las suposiciones de elasticidad lineal.

Resultados y Ejemplos Numéricos
Los autores validan el esquema utilizando varios ejemplos numéricos con diferentes complejidades de RVE:

• RVEs Cruz y Cuadrado: Se utilizaron RVEs planares simples (formas de cruz y cuadrado) para validar los resultados contra fórmulas analíticas existentes y literatura. El estudio demostró que incluir cizallamiento y extensión en las varillas a escala micro afecta significativamente la rigidez macroscópica, particularmente para varillas con mayores relaciones de esbeltez. Los resultados coincidieron con los límites analíticos para varillas Kirchhoff (inextensibles/indeformables).
• RVEs de Varillas Helicoidales: Se analizó un RVE cuadrado más complejo compuesto por varillas constituyentes helicoidales. Los resultados mostraron una respuesta de esfuerzo-estiramiento macroscópica en forma de J, transitando de un comportamiento dominado por flexión a uno dominado por estiramiento. El estudio destacó cómo parámetros como el número de vueltas helicoidales y el paso podían ajustar esta respuesta. El análisis de flexión reveló un comportamiento de ablandamiento desencadenado por la flexión fuera del plano de las varillas helicoidales de entrecruzamiento.
• Estructuras Tubulares Auxéticas: El esquema se aplicó a metamateriales tubulares auxéticos (relevantes para stents). El modelo capturó con éxito la dependencia de la rigidez de flexión macroscópica, el coeficiente de Poisson (comportamiento de acortamiento) y la rigidez radial con respecto al ángulo de diseño $\theta_a$. También demostró la captura del efecto Bazier (ovalización de la sección transversal) bajo grandes flexiones.
• Validación FE2: Una comparación entre el esquema de homogeneización propuesto (enfoque FE2) y simulaciones a escala completa de un metamaterial de RVE cruzado mostró un buen acuerdo, validando la precisión del método incluso bajo grandes momentos de torsión donde los modelos constitutivos lineales fallaron.

Significado y Afirmaciones
El artículo afirma proporcionar un marco robusto para la homogeneización de metamateriales tipo varilla que supera las limitaciones de los esquemas de homogeneización de primer orden y lineales. Al utilizar una teoría de varilla Cosserat especial geométricamente exacta en ambas escalas y derivar un mapa de transición macro-a-micro no lineal, el método predice con precisión la respuesta constitutiva no lineal de estas estructuras bajo grandes deformaciones.

Los autores enfatizan que su enfoque es general y puede extenderse a otros tipos de RVE (por ejemplo, sólidos porosos, estructuras de origami) y puede incorporar fenómenos microscópicos complejos como elastoplasticidad, viscoelasticidad y acoplamiento multifísico (termo-, quimio- o electroelasticidad). El trabajo sirve como base para diseñar metamateriales a medida con propiedades macroscópicas específicas, como músculos artificiales o stents desplegables, mediante el ajuste de parámetros microestructurales. Los autores notan modestamente que, aunque el trabajo actual asume que el RVE es una red de varillas, las condiciones de contorno derivadas podrían adaptarse para otras representaciones microestructurales.