An essay on $d_0$ in neutron-based strain measurement techniques: equilibrium-based methods for non-destructive $d_0$ estimation

Resumen Técnico: Métodos Basados en el Equilibrio para la Estimación No Destructiva de $d_0$ en la Medición de Deformación por Neutrones

El Problema: La Ambigüedad de $d_0$
En la medición de deformación mediante difracción de neutrones, el cálculo de la deformación elástica ($\epsilon$) se basa en la ecuación $\epsilon = (d - d_0)/d_0$, donde $d$ es el espaciado de red medido y $d_0$ es el espaciado de red de referencia libre de tensiones. Aunque a menudo se trata como una constante del material, $d_0$ frecuentemente no es constante ni intrínseco. Varía espacial y direccionalmente debido a gradientes de composición, transformaciones de fase, textura y tensiones residuales a escalas micro e intergranulares (tensiones de Tipo II y III).

La dificultad fundamental, denominada "el problema de $d_0$", es que las mediciones de difracción por sí solas no pueden distinguir entre un cambio en el espaciado de red de referencia ($d_0$) y una deformación elástica hidrostática uniforme. Matemáticamente, un conjunto de mediciones de difracción es insuficiente para resolver de forma única tanto el tensor de deformación como $d_0$ simultáneamente porque el sistema es deficiente en rango; una deformación hidrostática uniforme es observacionalmente indistinguible de un desplazamiento en $d_0$. En consecuencia, determinar un $d_0$ apropiado es una tarea crítica y no trivial. La solución convencional implica pruebas destructivas —cortar la muestra en pequeños cupones con alivio de tensiones para medir $d_0$ directamente—. Sin embargo, esto suele ser impráctico o imposible para artefactos valiosos, grandes o únicos donde el muestreo destructivo está prohibido.

Metodología: Una Jerarquía de Restricciones de Equilibrio
Este artículo propone una jerarquía de alternativas no destructivas basadas en el principio de equilibrio mecánico para estimar $d_0$ sin muestras de referencia destructivas. El enfoque aprovecha el hecho de que los campos de tensión residual dentro de una muestra deben satisfacer las ecuaciones de equilibrio ($\nabla \cdot \sigma = 0$) y las condiciones de contorno (por ejemplo, superficies libres de tracción).

La metodología progresa a través de tres niveles de restricción:

1. Balance de Fuerzas (Equilibrio Global):
   Para una sección transversal planar dentro de una muestra, la fuerza neta integrada sobre el área debe ser cero (en ausencia de cargas externas). Al integrar las deformaciones normales medidas a través de una sección transversal, se puede derivar una ecuación integral (Ecuación 15) para resolver un único valor de $d_0$ constante que satisfaga el balance de fuerzas global. Este método asume un $d_0$ isotrópico pero utiliza el conjunto completo de datos de medición de deformación para resolver la ambigüedad hidrostática.

2. Condiciones de Contorno (Equilibrio Local):
   En superficies libres donde el vector de tracción es cero ($\sigma \cdot n = 0$), los componentes de tensión deben anularse. Al extrapolar los datos de deformación medidos hacia el límite, se pueden aplicar estas condiciones de tensión cero para resolver valores de $d_0$ locales. Esto añade efectivamente ecuaciones adicionales al sistema de incógnitas, permitiendo la determinación única de $d_0$ en el contorno, que luego puede usarse para restringir el interior.

3. Equilibrio Punto a Punto mediante el Análisis de Eigenstrain:
   Para escenarios más complejos, el artículo utiliza el marco de trabajo de eigenstrain (deformación propia). La deformación total se descompone en deformación elástica ($\epsilon$) y deformación propia inelástica ($\epsilon^*$). Los autores formulan un problema inverso donde se reconstruyen simultáneamente una distribución de $d_0$ espacialmente variable y un campo de eigenstrain para ajustar los espaciados-$d$ observados.

   • Sistemas Axisimétricos: Para geometrías cilíndricas, el equilibrio se reduce a una ecuación diferencial ordinaria de primer orden (Ecuación 18). Esto permite la reconstrucción explícita de una distribución radial de $d_0(r)$ a partir de mediciones de espaciado de red radial, de aro y axial, siempre que se satisfaga la condición de contorno en la superficie.
   • Planos de Simetría y Sistemas 2D: Para mediciones planas (por ejemplo, planos de simetría en cubos o secciones transversales de palas), los autores construyen una base de campos de eigenstrain (usando funciones polinómicas). Resuelven un problema de optimización de mínimos cuadrados regularizado para encontrar los coeficientes de estas funciones base y los parámetros de una distribución de $d_0(x,y)$ que mejor se ajusten a los datos mientras satisfacen estrictamente el equilibrio y las condiciones de contorno. Este enfoque utiliza el modelo de eigenstrain como un regularizador físico para manejar la naturaleza mal planteada del problema inverso.

Resultados Clave y Casos de Estudio
El artículo valida estos métodos utilizando tres ejemplos distintos:

• Cubo de Inconel Fabricado Aditivamente:
  Se midió un cubo de Inconel fabricado mediante fusión selectiva por láser (SLM) mediante difracción de neutrones. Un cubo de muestra tipo "peine" con alivio de tensiones, cortado de un cubo idéntico, proporcionó mediciones directas de $d_0$ para comparación.

  • Resultados: Los valores de $d_0$ derivados del balance de fuerzas (1.085808 Å) y la extrapolación de la condición de contorno (1.085810 Å) fueron consistentes con las mediciones directas (promedio 1.085869 Å) y entre sí. La reconstrucción de equilibrio punto a punto arrojó un promedio de 1.085792 Å con una pequeña variación espacial. Esto confirmó que los métodos basados en el equilibrio pueden recuperar con precisión $d_0$ incluso cuando la medición directa es posible, sirviendo como una validación robusta.

• Sondas Médicas de Bronce de la Antigua Roma:
  Se analizaron cuatro sondas de bronce antiguas donde el muestreo destructivo era imposible.

  • Resultados: La suposición de un $d_0$ constante falló al satisfacer la condición de contorno superficial (prediciendo ~100 MPa de tensión radial en la superficie). La aplicación directa de la ecuación diferencial axisimétrica produjo una distribución de $d_0$ que satisfacía el equilibrio y las condiciones de contorno, pero resultó en una fuerza neta desequilibrada. Finalmente, la reconstrucción basada en eigenstrain satisfizo tanto la condición de contorno como el balance de fuerzas global simultáneamente. El $d_0$ resultante mostró una pequeña pero medible variación (rango ~3.9 × 10⁻⁴ Å), lo que sugiere que imponer las restricciones de equilibrio puede revelar variaciones espaciales en $d_0$ que las suposiciones constantes pasan por alto.

• Daga de Bronce Persa Antigua:
  Se escaneó una daga de bronce (c. 1000 a.C.) utilizando una combinación de imágenes y difracción de neutrones.

  • Resultados: Similar a las sondas, un $d_0$ constante fue insuficiente. La reconstrucción de eigenstrain de un $d_0$ espacialmente variable (rango ~±1.9 × 10⁻⁴ Å) produjo un campo de tensiones consistente con los procesos de fundición/temple y posibles trabajos en frío. El $d_0$ reconstruido mostró valores ligeramente mayores cerca de las superficies en comparación con el núcleo. Los residuos entre el modelo y los datos estuvieron dentro de la incertidancia típica del ajuste de picos, demostrando la viabilidad del enfoque para el análisis no destructivo de artefactos únicos.

Significado y Reivindicaciones
El artículo no pretende demostrar que $d_0$ varía definitivamente en estas muestras específicas, reconociendo que las variaciones inferidas son una de las muchas soluciones posibles consistentes con los datos. Sin embargo, los autores sostienen que su trabajo demuestra un cambio crítico de metodología:

1. Viabilidad No Destructiva: Las suposiciones de equilibrio proporcionan una restricción física práctica para estimar $d_0$ cuando las mediciones de referencia destructivas son totalmente inapropiadas.
2. Consistencia Física: Al tratar a $d_0$ no como una constante de calibración independiente, sino como una variable restringida por la mecánica de la muestra, los campos de deformación y tensión resultantes están garantizados para ser físicamente admisibles (satisfaciendo el equilibrio y las condiciones de contorno).
3. Jerarquía de Restricciones: El artículo establece un marco escalable que va desde balances de fuerzas simples hasta reconstrucciones complejas de eigenstrain punto a punto. Esta jerarquía permite a los investigadores suplementar, probar o reemplazar las mediciones directas dependiendo de la geometría de la muestra y la disponibilidad de simetrías o condiciones de contorno.

En conclusión, los autores argumentan que $d_0$ debe ser reconstruido utilizando la mecánica de la muestra siempre que sea posible, yendo más allá de la visión simplista de considerarlo una constante de material fija. Este enfoque transforma el problema de $d_0$ de un obstáculo de calibración en un problema inverso resoluble basado en las leyes fundamentales de la elasticidad.