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

Technische Samenvatting: Kwadratuurformules voor Integrale Vergelijkingen op een Semi-Oneindig Interval

Probleemstelling
Het artikel behandelt de uitdaging van het numeriek oplossen van integrale vergelijkingen gedefinieerd op een semi-oneindig interval $(0, \infty)$, specifiek die voortvloeien uit de wiskundige elasticiteitstheorie. Hoewel Gauss-type kwadratuurformules goed gevestigd zijn voor eindige intervallen, wordt hun toepassing op semi-oneindige domeinen bemoeilijkt door het gedrag van de oplossing bij oneindig. Bestaande methoden gebaseerd op Laguerre-polynomen leggen een voorwaarde van exponentiële afname ($x^\alpha e^{-x}$) op aan de oplossing, wat vaak fysiek onjuist is voor elasticiteitsproblemen waarbij de oplossing doorgaans een machtsfunctie-afname vertoont. De auteurs beogen een mechanische kwadratuurmethode te construeren die in staat is integralen met machtsfunctie-gewichtsfuncties te verwerken, inclusclusief reguliere integralen, Cauchy-type singuliere integralen en integralen met logaritmische en stapfunctie-kernels.

Methodologie
De kern van de methodologie betreft de constructie van een nieuwe familie van rationale functies, $R^{(p,q)}_n(x)$, die orthogonaal zijn op de straal $(0, \infty)$ met betrekking tot de gewichtsfunctie $(1+x)^{-p}x^{-q}$. Deze functies worden afgeleid van Jacobi-polynomen $P^{(p+q-2, -q)}_n$ via een lineaire fractionele transformatie die het interval $(-1, 1)$ naar $(0, \infty)$ brengt.

De auteurs volgen hierbij de volgende stappen:

1. Interpolatie: Een continue, begrensde functie $f(x)$ wordt geïnterpoleerd met behulp van de rationale functies $R^{(p,q)}_n(x)$.
2. Afleiding van Kwadratuur:
   • Reguliere Integralen: Door de interpolatieformule in de integraal te substitueren en gebruik te maken van de orthogonaliteitsvoorwaarden, wordt een kwadratuurformule afgeleid met knopen $z_i$ (wortels van het rationale polynoom) en specifieke gewichten $w_i$. Deze formule is exact voor polynomen in $\frac{x-1}{x+1}$ van graad tot en met $2n-1$.
   • Cauchy-type Integralen: Met behulp van de identiteit voor partiële breuken en de eigenschappen van de rationale functies wordt een kwadratuurformule voor de Cauchy-hoofdwaarde integraal geconstrueerd.
   • Singuliere Operatoren van de Tweede Soort: De methode wordt aangepast aan het karakteristieke deel van singuliere integrale vergelijkingen van de tweede soort. De auteurs bepalen de exponenten $p$ en $q$ op basis van de begrensdheid van de oplossing nabij de oorsprong en de specifieke kernelstructuur, waardoor de singuliere operator gediscretiseerd kan worden.
   • Zwakke Singulariteiten: Formules voor integralen met logaritmische kernels ($\ln|x-y|$) en stapfuncties ($\text{sign}(x-y)$) worden afgeleid met behulp van de Christoffel–Darboux-formule en specifieke integraalrepresentaties.

Toepassing en Validatie
Om de afgeleide formules te valideren, passen de auteurs de mechanische kwadratuurmethode toe op een tweedimensionaal elasticiteitsprobleem: de vervorming van een elastisch halfvlak.

• Fysieke Opstelling: Het halfvlak beslaat $y \le 0$ en is star bevestigd langs de semi-oneindige grens $x > 0$. Een uniform verdeelde normale last werkt op een eindig segment van de vrije grens.
• Wiskundige Reductie: De randvoorwaarden reduceren het probleem tot een singuliere integrale vergelijking van de tweede soort op het semi-oneindige interval. De onbekende contactspanningen worden gezocht in een vorm die consistent is met de afgeleide gewichtsfuncties.
• Discretisatie: De integrale vergelijking wordt omgezet in een stelsel lineaire algebraïsche vergelijkingen met gebruik van de voorgestelde kwadratuurknopen en gewichten. Het stelsel wordt opgelost voor de onbekende spanningscoëfficiënten.
• Nabewerking: Zodra de spanningsverdeling is bepaald, worden de formules voor integralen met logaritmische en stapkernels gebruikt om de verplaatsingscomponenten te berekenen en de vervormde vorm van de grens te reconstrueren.

Belangrijkste Resultaten

• Convergentie: Numerieke analyse toont een uitstekende convergentie van het algoritme aan. Voor een Poisson-ratio $\nu = 0,27$ en een lastparameter $a=1$ stabiliseert de oplossing snel, waarbij $n=6$ knopen vaak voldoende zijn voor een hoge nauwkeurigheid. Zelfs voor kleine waarden van $a$ (waarbij de last dicht bij de vaste regio ligt), levert het verhogen van de orde naar $n=20$ stabiele resultaten op.
• Nauwkeurigheid: De berekende waarden van de spanningsfunctie $\phi^*(x)$ op diverse punten vertonen een opmerkelijke stabiliteit naarmate de benaderingsorde $n$ toeneemt, wat de hoge algebraïsche precisiegraad van de voorgestelde kwadratuurregels bevestigt.
• Fysisch Inzicht: De methode bepaalt succesvol de contactspanningen en de vervormde grens vorm van het elastische halfvlak, wat de effectiviteit van de aanpak voor praktische elasticiteitsproblemen illustreert.

Betekenis en Claims
De auteurs stellen dat dit werk, "blijkbaar voor het eerst", kwadratuurformules presenteert die specifiek zijn ontworpen voor het evalueren van integralen met machtsfunctie-gewichtsfuncties op een semi-oneindig interval. In tegenstelling tot eerdere benaderingen die steunen op Laguerre-polynomen, faciliteert deze methode de machtsfunctie-afname die typerend is voor elasticiteitsproblemen zonder een onrealistische exponentiële afname op te leggen. Het artikel vestigt een verenigd kader waarbij dezelfde set knopen en gewichten kan worden toegepast op reguliere integralen, Cauchy-type singuliere integralen en integralen met zwakke singulariteiten (logaritmische en stapkernels). De succesvolle toepassing op het contactprobleem van een vast elastisch halfvlak dient als een rigoureuze validatie van de theoretische afleidingen.