Análisis de Tensiones Residuales Elasto-Plásticas en Materiales Porosos SLS mediante Simulaciones de Campo de Fase 3D Multicapa

2.1. Modelo 3D Multicapa de Campo de Fase

Un modelo de campo de fase simula la consolidación del polvo y la evolución de la microestructura. El parámetro de orden $\phi$ distingue entre las fases sólida ($\phi=1$) y polvo/vacío ($\phi=0$). La evolución se rige por la ecuación de Allen-Cahn, acoplada con la transferencia de calor:

$\frac{\partial \phi}{\partial t} = -L \frac{\delta F}{\delta \phi}$

donde $L$ es el coeficiente cinético y $F$ es el funcional de energía libre total, que incorpora energía de gradiente y un potencial de doble pozo. El modelo tiene en cuenta la fusión parcial y la formación de cuellos entre partículas.

2.2. Marco de Elementos Finitos Termo-Elasto-Plástico

El historial térmico de la simulación de campo de fase impulsa un análisis mecánico. El modelo constitutivo considera propiedades del material dependientes de la temperatura y la fase. La deformación total $\boldsymbol{\varepsilon}$ se descompone en componentes elástica $\boldsymbol{\varepsilon}^e$, plástica $\boldsymbol{\varepsilon}^p$ y térmica $\boldsymbol{\varepsilon}^{th}$:

$\boldsymbol{\varepsilon} = \boldsymbol{\varepsilon}^e + \boldsymbol{\varepsilon}^p + \boldsymbol{\varepsilon}^{th}$

La plasticidad se modela utilizando una teoría de flujo $J_2$ con endurecimiento isotrópico. La tensión se calcula mediante la ley de Hooke: $\boldsymbol{\sigma} = \mathbb{C}(T, \phi) : \boldsymbol{\varepsilon}^e$, donde $\mathbb{C}$ es el tensor de rigidez dependiente de la temperatura $T$ y la fase $\phi$.

2.3. Propiedades del Material y Condiciones de Contorno

Las propiedades del material (módulo de Young, límite elástico, coeficiente de expansión térmica) se definen como funciones de la temperatura y la fase. El dominio de simulación representa una sección del lecho de polvo con condiciones de contorno térmicas (convección, radiación) y mecánicas (inicialmente libres de tensión) apropiadas. El haz láser se modela como una fuente de calor móvil con una distribución de potencia gaussiana.

3.1. Evolución Microestructural y Concentración de Tensiones

Las simulaciones revelan que la concentración de tensiones no es uniforme. Se localiza principalmente en dos regiones críticas: (1) las zonas de cuello entre partículas parcialmente fundidas, y (2) las unión entre capas depositadas sucesivamente. Estas regiones actúan como puntos calientes para el inicio de la plasticidad.

3.2. Distribución de Tensiones Residuales y Deformación Plástica

La deformación plástica acumulada y la tensión residual resultante siguen el gradiente térmico. Se encuentran altas tensiones residuales de tracción en el núcleo de las regiones sinterizadas durante el enfriamiento, mientras que pueden desarrollarse tensiones de compresión cerca de la superficie o las interfaces. Este desajuste es un impulsor principal de la distorsión y la posible delaminación.

3.3. Efecto de los Parámetros del Proceso

El estudio evalúa el impacto de la potencia del haz y la velocidad de escaneo (a menudo combinadas como densidad de energía lineal). Una mayor entrada de energía generalmente aumenta el tamaño del baño de fusión, reduce la porosidad, pero puede exacerbar los gradientes térmicos, lo que lleva a una mayor tensión residual si no se optimiza. Se identifica una ventana para lograr una densidad aceptable con una tensión minimizada.

4.1. Idea Central y Flujo Lógico

La idea central del artículo es que las simulaciones mesoscópicas, resueltas a nivel de polvo, son no negociables para predecir los resultados del SLS. El flujo lógico es robusto: comenzar con un modelo de campo de fase basado en la física para obtener la "forma" del material, luego usar esa geometría y el historial térmico como entrada directa para un análisis mecánico. Este es un paso significativo frente a los modelos homogeneizados que difuminan los fenómenos críticos a escala de polvo. Refleja la filosofía en la ciencia computacional de materiales donde capturar la microestructura correcta es primordial, similar al trabajo sobre fractura de campo de fase de Miehe et al. (2010).

4.2. Fortalezas y Debilidades

Fortalezas: El acoplamiento es elegante y físicamente sólido. Identificar las regiones de cuello e intercapas como concentradores de tensiones es un hallazgo valioso y accionable. El intento de construir modelos de regresión a partir de datos de simulación es pragmático para la adopción industrial.

Debilidades: El costo computacional es sin duda masivo, limitando el tamaño del dominio simulado. Es probable que el modelo simplifique la distribución del tamaño del polvo y la aleatoriedad del empaquetamiento. La validación, aunque se compara con experimentos, necesita mediciones de tensión cuantitativas y punto por punto más detalladas (por ejemplo, usando difracción de rayos X de sincrotrón) para establecer plenamente su credibilidad.

4.3. Conclusiones Accionables

Para ingenieros de proceso: Enfocar la planificación de la trayectoria del láser y el desarrollo de parámetros no solo en la densidad, sino en minimizar la tensión en las regiones críticas identificadas (cuellos, uniones de capa). Usar los modelos de regresión como punto de partida para definir una ventana de proceso. Para investigadores: Este marco es una plantilla. El siguiente paso es incorporar más física (dinámica de esparcido de polvo, vaporización y formación explícita de grietas) y aprovechar la reducción de orden del modelo o sustitutos de aprendizaje automático para hacerlo más rápido.