Revisiting Stress Analysis in a Three-Dimensional Elastic Hollow Sphere under Uniaxial Compression via the Inverse Laplace Transform Expressions within an Elastodynamic Framework

Resumen Técnico: Revisitando el Análisis de Esfuerzos en una Esfera Hueca Elástica Tridimensional bajo Compresión Uniaxial

Planteamiento del Problema
Este estudio aborda el análisis de esfuerzos de una esfera hueca elástica tridimensional sometida a compresión uniaxial. Si bien el comportamiento de las esferas sólidas bajo cargas similares está bien documentado en la teoría de la elasticidad clásica, y existen análisis dinámicos para cilindros y esferas, persiste una laguna en la literatura respecto a soluciones analíticas explícitas y cerradas para esferas huecas. Trabajos previos, como el de Kasano et al., abordaron esferas huecas pero determinaron los coeficientes de desplazamiento y esfuerzo de forma incremental, careciendo de la claridad y utilidad de las expresiones explícitas. Este artículo tiene como objetivo llenar ese vacío mediante la derivación de expresiones analíticas explícitas para los coeficientes de desplazamiento y esfuerzo, mejorando así la comprensión teórica del comportamiento de esfuerzos en estructuras esféricas porosas o huecas.

Metodología
Los autores emplean un marco elastodinámico para modelar el sistema. El problema se define por una esfera hueca elástica con un radio exterior $R_o$ y un radio interior $R_i$, caracterizada por una densidad de masa $\varrho_0$, un módulo de corte $G$ y una relación de Poisson $\nu$. La esfera está sometida a una compresión uniaxial definida por una presión $\sigma_0$ que actúa diametralmente.

El proceso de solución involucra los siguientes pasos:

1. Ecuaciones Gobernantes: La ecuación de Navier-Cauchy adimensional se resuelve utilizando el teorema de Helmholtz, el cual descompone el vector de desplazamiento en potenciales escalar ($\varphi$) y vectorial ($\chi$).
2. Transformada de Laplace: Los potenciales se transforman al dominio de Laplace. Las soluciones se expanden utilizando funciones de Bessel modificadas de primera y segunda clase ($i_n$ y $k_n$) y polinomios de Legendre ($P_n$).
3. Condiciones de Contorno: Los coeficientes de la expansión se determinan satisfaciendo las condiciones de contorno en las superficies interna y externa. Esto genera un sistema de ecuaciones simultáneas, cuya solución proporciona expresiones explícitas para los coeficientes $a_n, b_n, c_n$ y $d_n$ en términos de determinantes que involucran las funciones de Bessel y los parámetros del material.
4. Límite Estático: Para obtener soluciones estáticas, los autores aplican el teorema del valor final a las expresiones de la transformada inversa de Laplace. Esto permite la derivación de series infinitas explícitas para los componentes de desplazamiento y esfuerzo estáticos a medida que el parámetro de tiempo tiende a infinito.
5. Implementación Numérica: Se analiza la convergencia de la serie infinita, estableciendo que un límite de truncamiento de $N_{max} = 1000$ es suficiente para un cálculo preciso, incluso cerca de la superficie interna.

Resultos Clave
El estudio deriva expresiones de series explícitas para los desplazamientos radial y tangencial adimensionales ($\tilde{u}_\rho, \tilde{u}_\theta$) y los componentes de esfuerzo ($\tilde{\sigma}_{\rho\rho}, \tilde{\sigma}_{\rho\theta}, \tilde{\sigma}_{\theta\theta}, \tilde{\sigma}_{\phi\phi}$). Los hallazgos clave incluyen:

• Concentración de Esfuerzos: La magnitud de la deformación y la diferencia del segundo esfuerzo principal se ven incrementadas en las proximidades de los puntos de carga (polos).
• Comparación con Esferas Sólidas: En el eje de carga ($\theta=0$), el esfuerzo de compresión en una esfera hueca se comporta de manera similar a una esfera sólida en la región del volumen (bulk), pero disminuye drásticamente a medida que se aproxima a la superficie interna. Por el contrario, el esfuerzo de extensión es menor que el de una esfera sólida en la región intermedia, pero aumenta significativamente cerca del diámetro interno.
• Comportamiento de la Superficie Interna: En la superficie interna ($\rho = \rho_i$), el esfuerzo circunferencial $\tilde{\sigma}_{\theta\theta}$ exhibe un pico distintivo en un ángulo $\theta$ entre 0 y $\pi/2$ (específicamente donde el ángulo entre la fuerza de compresión y el punto de la superficie es perpendicular). Este pico se vuelve más pronunciado a medida que el radio interno aumenta.
• Límite de Capa Delgada: En el límite donde el ancho de la capa es pequeño ($\delta \ll 1$), los componentes de desplazamiento radial y de esfuerzo escalan aproximadamente como $(1 - \rho_i)^{-1}$, lo cual es consistente con las expectativas de la teoría de capas delgadas.

Significancia y Reivindicaciones
Los autores afirman que este trabajo proporciona un enfoque analítico robusto para comprender y predecir el comportamiento mecánico de esferas elásticas huecas bajo carga uniaxial. Al pasar de expresiones incrementales a expresiones cerradas explícitas, el estudio ofrece herramientas prácticas para:

• Simplificar el modelado computacional de estructuras porosas.
• Ayudar en el diseño de materiales con propiedades mecánicas a medida.
• Mejorar la precisión de las predicciones de falla en aplicaciones de ingeniería que involucran geometrías esféricas huecas o porosas.

El artículo concluye que estos hallazgos contribuyen a una comprensión más profunda de la distribución de esfuerzos en geometrías internas complejas. Los autores señalan que, si bien el enfoque actual se centra en soluciones estáticas, las expresiones de la transformada inversa de Laplace derivadas pueden extenderse a problemas dinámicos para analizar fenómenos de propagación de ondas (ondas longitudinales, transversales y de Rayleigh), aunque dicho análisis dinámico requeriría una evaluación compleja de residuos y una categorización de casos basada en el tiempo y la ubicación.