An implicit octree-based adaptive Material Point Method

Resumen Técnico: Un Método de Punto Material Adaptativo basado en Octree Implícito

Planteamiento del Problema
El Método de Punto Material (MPM, por sus siglas en inglés) es un marco robusto para simular grandes deformaciones, particularmente en interacciones entre estructuras continuas y rígidas, como los problemas de interacción suelo-estructura. Sin embargo, las implementaciones estándar de MPM dependen de mallas de fondo estructuradas, lo que requiere una alta resolución en todo el dominio para capturar con precisión las interfaces evolutivas entre cuerpos rígidos y continuos deformables. Este requisito conduce a sistemas numéricos de gran tamaño y tiempos de ejecución computacional excesivos, especialmente en tres dimensiones. Además, los esquemas de integración temporal implícitos, que son ventajosos para problemas de ingeniería de larga duración (por ejemplo, la interacción suelo-estructura), sufren de una mala escalabilidad de los resolvedores lineales dispersos a medida que aumenta el número de grados de libertad. El método también enfrenta inestabilidades numéricas bien conocidas: la inestabilidad de cruce de celdas (cell-crossing instability) y el problema del "corte pequeño" (small-cut problem), donde los elementos son parcialmente cortados por las fronteras, provocando sistemas mal condicionados.

Metodología
Este artículo presenta el primer GIMPM (Método de Punto Material de Interpolación Generalizada) basado en un octree implícito. El núcleo del enfoque es una malla de fondo adaptativa que refina y coarce automáticamente basándose en la posición del cuerpo rígido que interactúa con el continuo.

• Discretización por Octree: La malla de fondo utiliza una estructura de octree para concentrar el esfuerzo computacional espacialmente. Los elementos se refinan solo en las regiones que requieren alta fidelidad (cerca de la interfaz), mientras que se mantienen elementos gruesos en otras zonas. Esto introduce nodos colgantes (hanging nodes) (nodos en la cara o el borde de un vecino más grueso), lo que rompe la continuidad $C^0$ estándar.
• Formulación de la Base GIMP Discreta: Para manejar la malla de octree no conforme, los autores formulan una base GIMP discreta. En lugar de construir funciones complejas por trozos para todas las posibles variaciones de intersección, las funciones de base se calculan mediante una única regla de cuadratura de puntos de Gauss sobre la intersección del dominio del punto material y los elementos de fondo. Este enfoque incorpora naturalmente las ecuaciones de restricción necesarias para vincular los nodos colgantes con sus nodos de vértice de soporte, manteniendo la continuidad $C^1$ de las funciones de base.
• Estabilización por Fantasma (Ghost Stabilisation): El problema del corte pequeño, que causa un mal condicionamiento en dominios no ajustados (unfitted domains), se mitiga extendiendo las técnicas de estabilización de fantasma a las mallas de octree. Esto implica añadir una forma bilineal de fantasma a la matriz de rigidez a través de "caras fantasma" (interfaces entre elementos activos e inactivos o de frontera), asegurando la estabilidad sin necesidad de técnicas de agregación que podrían complicar el marco de trabajo.
• Estrategia Adaptativa: La adaptividad es impulsada por la posición del cuerpo rígido. La estrategia consiste en tres pasos: identificar los elementos que intersectan el cuerpo rígido, refinar los elementos dentro de una distancia especificada ($dx_{region}$) y refinar los puntos materiales (GIMPs) mismos basándose en el tamaño local del elemento. Los datos del material se transfieren de los GIMPs gruesos a los refinados, aunque los autores señalan que, para propiedades de material dependientes de la posición, esta transferencia depende de coordenadas de referencia en lugar de un marco energético.

Contribuciones Clave

1. Primer GIMPM de Octree Implícito: El artículo introduce el primer GIMPM que utiliza una malla de fondo de tipo octree, permitiendo la solución eficiente de problemas de interacción continuo-estructura con grandes deformaciones.
2. Formulación de Base Discreta: Un nuevo método de integración discreta para las funciones de base GIMP en mallas no conformes, evitando la necesidad de acoplamientos de malla jerárquicos o definiciones complejas por trozos.
3. Estabilización por Fantasma Extendida: La adaptación de la estabilización por fantasma a mallas de octree con nodos colgantes para resolver la inestabilidad del corte pequeño.
4. Adaptividad Automática: Una estrategia sistemática que adapta tanto la malla computacional como la discretización del punto material basándose en la interfaz de contacto, eliminando la necesidad de algoritmos separados de desrefinamiento.

Resultados
El método fue validado mediante cinco ejemplos numéricos:

• Estudios de Convergencia: Una columna bajo su propio peso y una columna bajo esfuerzo de contacto constante demostraron la convergencia de las soluciones de esfuerzo con el refinamiento de la malla, validando la formulación de nodos colgantes y el enfoque de penalización de contacto.
• Contacto Dinámico: Una esfera rodando por una pendiente validó la implementación de fricción y la formulación de contacto en presencia de nodos colgantes, mostrando una excelente concordancia con las soluciones analíticas para diversos coeficientes de fricción.
• Aplicaciones Geotécnicas:
  • Ensayo de Penetración de Cono (CPT): La malla de octree adaptativa reprodujo con éxito los resultados experimentales del CPT. El estudio destacó que, si bien el refinamiento aumenta la rigidez en esta configuración específica (deb debido a la definición del módulo de Young en los centros de los GIMPs), los resultados convergieron con el refinamiento.
  • Penetración de Anclaje: Se simuló un anclaje arrastrado a través de un lecho de arena y se comparó con datos experimentales y una simulación previa de malla estructurada.

Desempeño e Impacto Ambiental
En el caso de la penetración del anclaje, el enfoque de octree adaptativo logró resultados comparables o mejores que una malla conforme de alta resolución.

• Velocidad: El método adaptativo logró una aceleración de hasta 5.5 veces en comparación con la simulación de malla estructurada.
• Eficiencia: La reducción de los grados de libertad y de los puntos materiales condujo a una disminución significativa de los recursos computacionales.
• Huella de Carbono: El artículo reporta una reducción de 21 veces en las emisiones de CO2 para la simulación adaptativa en comparación con su equivalente de malla estructurada, calculada en función del uso de recursos de Computación de Alto Rendimiento (HPC).

Significancia
El artículo afirma que este trabajo representa un paso significativo hacia la eficiencia del MPM implícito para problemas de grandes deformaciones. Al integrar la adaptatividad de octree con el GIMPM, el método concentra la resolución computacional estrictamente en la interfaz evolutiva, reduciendo el número de grados de libertad y de puntos materiales. Esto no solo acelera las simulaciones, sino que también ofrece beneficios ambientales sustanciales al reducir el consumo energético requerido para el modelado geotécnico de alta fidelidad. Los autores enfatizan que el método requiere una intervención mínima del usuario (definiendo el tamaño mínimo del elemento y la región de refinamiento) y proporciona un marco robusto para complejas interacciones suelo-estructura rígida.