Quadrature Formulas Applicable to the Solution of Integral Equations on a Semi-Infinite Interval by the Mechanical Quadrature Method

Resumo Técnico: Fórmulas de Quadratura para Equações Integrais em um Intervalo Semi-Infinito

Enunciado do Problema
O artigo aborda o desafio de resolver numericamente equações integrais definidas em um intervalo semi-infinito $(0, \infty)$, especificamente aquelas que surgem na teoria matemática da elasticidade. Embora fórmulas de quadratura do tipo Gauss sejam bem estabelecidas para intervalos finitos, sua aplicação a domínios semi-infinitos é complicada pelo comportamento da solução no infinito. Métodos existentes baseados em polinômios de Laguerre impõem uma condição de decaimento exponencial ($x^\alpha e^{-x}$) sobre a solução, o que é frequentemente impreciso do ponto de vista físico para problemas de elasticidade, onde a solução tipicamente decai como uma função de potência. Os autores visam construir um método de quadratura mecânica capaz de lidar com integrais com funções de peso de lei de potência, incluindo integrais regulares, integrais singulares do tipo Cauchy e integrais com núcleos logarítmicos e de função degrau.

Metodologia
O núcleo da metodologia envolve a construção de uma nova família de funções racionais, $R^{(p,q)}_n(x)$, ortogonais no raio $(0, \infty)$ em relação à função de peso $(1+x)^{-p}x^{-q}$. Essas funções são derivadas de polinômios de Jacobi $P^{(p+q-2, -q)}_n$ por meio de uma transformação fracionária linear que mapeia o intervalo $(-1, 1)$ para $(0, \infty)$.

Os autores procedem através das seguintes etapas:

1. Interpolação: Uma função contínua e limitada $f(x)$ é interpolada usando as funções racionais $R^{(p,q)}_n(x)$.
2. Derivação da Quadratura:
   • Integrais Regulares: Ao substituir a fórmula de interpolação na integral e utilizar as condições de ortogonalidade, deriva-se uma fórmula de quadratura com nós $z_i$ (raízes do polinômio racional) e pesos específicos $w_i$. Esta fórmula é exata para polinômios em $\frac{x-1}{x+1}$ de grau até $2n-1$.
   • Integrais do Tipo Cauchy: Usando a identidade para frações parciais e as propriedades das funções racionais, uma fórmula de quadratura para a integral do valor principal de Cauchy é construída.
   • Operadores Singulares de Segunda Espécie: O método é adaptado para a parte característica de equações integrais singulares de segunda espécie. Os autores determinam os expoentes $p$ e $q$ com base na limitação da solução próximo à origem e na estrutura específica do núcleo, permitindo que o operador singular seja discretizado.
   • Singularidades Fracas: Fórmulas para integrais envolvendo núcleos logarítmicos ($\ln|x-y|$) e funções degrau ($\text{sign}(x-y)$) são derivadas usando a fórmula de Christoffel–Darboux e representações integrais específicas.

Aplicação e Validação
Para validar as fórmulas derivadas, os autores aplicam o método de quadratura mecânica a um problema de elasticidade bidimensional: a deformação de um semiplano elástico.

• Configuração Física: O semiplano ocupa $y \le 0$ e é rigidamente fixado ao longo da fronteira semi-infinita $x > 0$. Uma carga normal uniformemente distribuída atua sobre um segmento finito da fronteira livre.
• Redução Matemática: As condições de contorno reduzem o problema a uma equação integral singular de segunda espécie no intervalo semi-infinito. Busca-se as tensões de contato desconhecidas em uma forma consistente com as funções de peso derivadas.
• Discretização: A equação integral é convertida em um sistema de equações algébricas lineares usando os nós e pesos de quadratura propostos. O sistema é resolvido para os coeficientes de tensão desconhecidos.
• Pós-processamento: Uma vez determinada a distribuição de tensão, as fórmulas para integrais com núcleos logarítmicos e de degrau são usadas para calcular os componentes de deslocamento e reconstruir a forma deformada da fronteira.

Principais Resultados

• Convergência: A análise numérica demonstra uma excelente convergência do algoritmo. Para um coeficiente de Poisson $\nu = 0,27$ e um parâmetro de carga $a=1$, a solução estabiliza-se rapidamente, sendo que $n=6$ nós são frequentemente suficientes para alta precisão. Mesmo para valores pequenos de $a$ (onde a carga está próxima da região fixa), aumentar a ordem para $n=20$ produz resultados estáveis.
• Precisão: Os valores computados da função de tensão $\phi^*(x)$ em vários pontos mostram uma estabilidade notável conforme a ordem de aproximação $n$ aumenta, confirmando o alto grau algébrico de precisão das regras de quadratura propostas.
• Insight Físico: O método determina com sucesso as tensões de contato e a forma deformada da fronteira do semiplano elástico, ilustrando a eficácia da abordagem para problemas práticos de elasticidade.

Significância e Alegações
Os autores alegam que este trabalho apresenta, "aparentemente pela primeira vez", fórmulas de quadratura especificamente projetadas para avaliar integrais com funções de peso de lei de potência em um intervalo semi-infinito. Diferentemente de abordagens anteriores baseadas em polinômios de Laguerre, este método acomoda o decaimento de lei de potência típico dos problemas de elasticidade sem impor um decaimento exponencial irrealista. O artigo estabelece um framework unificado onde o mesmo conjunto de nós e pesos pode ser aplicado a integrais regulares, integrais singulares do tipo Cauchy e integrais com singularidades fracas (núcleos logarítmicos e de degrau). A aplicação bem-sucedida ao problema de contato de um semiplano elástico fixo serve como uma validação rigorosa das derivações teóricas.