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

Resumo Técnico: Revisitando a Análise de Tensão em uma Esfera Oca Elástica Tridimensional sob Compressão Uniaxial

Definição do Problema
Este estudo aborda a análise de tensão de uma esfera oca elástica tridimensional sujeita à compressão uniaxial. Embora o comportamento de esferas sólidas sob carregamentos semelhantes seja bem documentado na teoria da elasticidade clássica, e existam análises dinâmicas para cilindros e esferas, permanece uma lacuna na literatura quanto a soluções analíticas explícitas e fechadas abrangentes para esferas ocas. Trabalhos anteriores, como o de Kasano et al., abordaram esferas ocas, mas determinaram os coeficientes de deslocamento e tensão de forma incremental, carecendo da clareza e utilidade de expressões explícitas. Este artigo visa preencher essa lacuna ao derivar expressões analíticas explícitas para os coeficientes de deslocamento e tensão, aumentando assim a compreensão teórica do comportamento de tensão em estruturas esféricas porosas ou ocas.

Metodologia
Os autores empregam uma estrutura elastodinâmica para modelar o sistema. O problema é definido por uma esfera oca elástica com raio externo $R_o$ e raio interno $R_i$, caracterizada por densidade de massa $\varrho_0$, módulo de cisalhamento $G$ e razão de Poisson $\nu$. A esfera é submetida a uma compressão uniaxial definida por uma pressão $\sigma_0$ atuando diametralmente.

O processo de solução envolve as seguintes etapas:

1. Equações Governantes: A equação de Navier-Cauchy adimensional é resolvida utilizando o teorema de Helmholtz, que decompõe o vetor deslocamento em potenciais escalar ($\varphi$) e vetorial ($\chi$).
2. Transformada de Laplace: Os potenciais são transformados para o domínio de Laplace. As soluções são expandidas usando funções de Bessel modificadas de primeira e segunda espécie ($i_n$ e $k_n$) e polinômios de Legendre ($P_n$).
3. Condições de Contorno: Os coeficientes da expansão são determinados satisfazendo as condições de contorno nas superfícies interna e externa. Isso resulta em um sistema de equações simultâneas, cuja solução fornece expressões explícitas para os coeficientes $a_n, b_n, c_n$ e $d_n$ em termos de determinantes envolvendo as funções de Bessel e parâmetros de material.
4. Limite Estático: Para obter soluções estáticas, os autores aplicam o teorema do valor final à expressão da transformada inversa de Laplace. Isso permite a derivação de séries infinitas explícitas para os componentes de deslocamento e tensão estáticos conforme o parâmetro de tempo se aproxima do infinito.
5. Implementação Numérica: A convergência das séries infinitas é analisada, estabelecendo que um limite de truncamento de $N_{max} = 1000$ é suficiente para um cálculo preciso, mesmo próximo à superfície interna.

Principais Resultados
O estudo deriva expressões de série explícitas para os deslocamentos radial e tangencial adimensionais ($\tilde{u}_\rho, \tilde{u}_\theta$) e componentes de tensão ($\tilde{\sigma}_{\rho\rho}, \tilde{\sigma}_{\rho\theta}, \tilde{\sigma}_{\theta\theta}, \tilde{\sigma}_{\phi\phi}$). As principais descobertas incluem:

• Concentração de Tensão: A magnitude da deformação e a diferença da segunda tensão principal são intensificadas nas proximidades dos pontos de carregamento (polos).
• Comparação com Esferas Sólidas: No eixo de carregamento ($\theta=0$), a tensão de compressão em uma esfera oca comporta-se de forma semelhante a uma esfera sólida na região do volume (bulk), mas diminui bruscamente à medida que a superfície interna é aproximada. Por outro lado, a tensão de extensão é inferior à de uma esfera sólida na região intermediária, mas aumenta significativamente próximo ao diâmetro interno.
• Comportamento da Superfície Interna: Na superfície interna ($\rho = \rho_i$), a tensão circunferencial $\tilde{\sigma}_{\theta\theta}$ exibe um pico distinto em um ângulo $\theta$ entre 0 e $\pi/2$ (especificamente onde o ângulo entre a força de compressão e o ponto da superfície torna-se perpendicular). Este pico torna-se mais pronunciado conforme o raio interno aumenta.
• Limite de Casca Fina: No limite em que a largura da casca é pequena ($\delta \ll 1$), os componentes de deslocamento radial e tensão escalam aproximadamente como $(1 - \rho_i)^{-1}$, consistente com as expectativas da teoria de cascas finas.

Significância e Alegações
Os autores afirmam que este trabalho fornece uma abordagem analítica robusta para compreender e prever o comportamento mecânico de esferas ocas elásticas sob carregamento uniaxial. Ao passar de formas incrementais para expressões explícitas de forma fechada, o estudo oferece ferramentas práticas para:

• Simplificar a modelagem computacional de estruturas porosas.
• Auxiliar no design de materiais com propriedades mecânicas customizadas.
• Melhorar a precisão das previsões de falha em aplicações de engenharia envolvendo geometrias esféricas ocas ou porosas.

O artigo conclui que esses achados contribuem para uma compreensão mais profunda da distribuição de tensão em geometrias internas complexas. Os autores observam que, embora o foco atual sejam as soluções estáticas, as expressões da transformada inversa de Laplace derivadas podem ser estendidas para problemas dinâmicos para analisar fenômenos de propagação de ondas (ondas longitudinais, transversais e de Rayleigh), embora tal análise dinâmica exija uma avaliação complexa de resíduos e categorização de casos baseada no tempo e na localização.