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

技術要約：半無限区間における積分方程式のための求積公式

問題の定義
本論文は、弾性理論の数学的理論において特によく見られる、半無限区間 $(0, \infty)$ で定義される積分方程式を数値的に解く際の課題に対処している。ガウス型求積公式は有限区間においては確立されているが、半無限領域への適用は、解の無限遠における振る舞いによって困難が生じる。ラグーエル多項式に基づく既存の手法は、解に対して指数関数的な減衰条件（$x^\alpha e^{-x}$）を課すが、これは弾性学の問題において解が通常べき関数として減衰するという物理的な実態と一致しないことが多い。著者らは、べき乗型の重み関数を持つ積分（通常の積分、コーシー型の特異積分、および対数やステップ関数のカーネルを持つ積分を含む）を扱うことが可能な、新しい機械的な求積法を構築することを目的としている。

手法
手法の核心は、重み関数 $(1+x)^{-p}x^{-q}$ に関して $(0, \infty)$ 上で直交する、新しい有理関数族 $R^{(p,q)}_n(x)$ の構築にある。これらの関数は、区間 $(-1, 1)$ を $(0, \infty)$ へ写像する線形分数変換を通じて、ヤコビ多項式 $P^{(p+q-2, -q)}_n$ から導出される。

著者らは以下の手順を進める：

1. 補間： 連続かつ有界な関数 $f(x)$ を、有理関数 $R^{(p,q)}_n(x)$ を用いて補間する。
2. 求積法の導出：
   • 通常の積分： 補間公式を積分に代入し、直交条件を利用することで、ノード $z_i$（有理多項式の根）と特定の重み $w_i$ を持つ求積公式を導出する。この公式は、$\frac{x-1}{x+1}$ の次数が $2n-1$ までの多項式に対して正確である。
   • コーシー型積分： 部分分数分解の恒等式と有理関数の性質を用いて、コーシー主値積分に対する求積公式を構築する。
   • 第二種特異演算子： 本手法を第二種特異積分方程式の特性部分に適応させる。著者らは、原点付近での解の有界性と特定のカーネル構造に基づいて指数 $p$ および $q$ を決定し、特異演算子の離散化を可能にする。
   • 弱特異性： クリストッフェル・ダルブー公式および特定の積分表示を用いて、対数カーネル（$\ln|x-y|$）およびステップ関数（$\text{sign}(x-y)$）を含む積分の公式を導出する。

適用と検証
導出された公式を検証するために、著者らは二次元弾性学の問題である「弾性半平面の変形」にこの機械的求積法を適用する。

• 物理的設定： 半平面は $y \le 0$ に位置し、半無限境界 $x > 0$ に沿って剛体に固定されている。自由境界の有限セグメントに対して、一様に分布した垂直荷重が作用する。
• 数学的還元： 境界条件により、問題は半無限区間上の第二種特異積分方程式へと還元される。未知の接触応力は、導出された重み関数と整合する形式で求められる。
• 離散化： 提案された求積ノードと重みを用いて、積分方程式を線形代数方程式の系へと変換する。未知の応力係数を求めるために、この系を解く。
• 後処理： 応力分布が決定された後、対数およびステップカーネルを用いた積分公式を使用して、変位成分を計算し、境界の変形形状を再構成する。

主な結果

• 収束性： 数値解析により、アルゴリズムの優れた収束性が示された。ポアソン比 $\nu = 0.27$、荷重パラメータ $a=1$ の場合、解は急速に安定し、$n=6$ 個のノードで十分な精度が得られることが多い。$a$ が小さい値（荷重が固定領域に近い場合）であっても、次数を $n=20$ まで上げれば安定した結果が得られる。
• 精度： 様々な点における応力関数 $\phi^*(x)$ の計算値は、近似次数 $n$ の増加に伴い極めて高い安定性を示しており、提案された求積則の高い代数的精度を裏付けている。
• 物理的洞察： 本手法は、弾性半平面の接触応力と変形した境界形状を正常に決定しており、実用的な弾性学の問題に対する本アプローチの有効性を実証している。

意義と主張
著者らは、本研究が「おそらく初めて」となる、べき乗型の重み関数を持つ積分を評価するために特別に設計された求積公式を提示していると主張している。ラグーエル多項式に依存する従来のアプローチとは異なり、本手法は、弾性学に典型的なべき乗型の減衰を、非現実的な指数関数的減衰を課すことなく許容する。本論文は、同じノードと重みのセットを、通常の積分、コーシー型特異積分、および弱特異性（対数およびステップカーネル）を持つ積分に適用できる統一的な枠組みを確立している。固定された弾性半平面の接触問題への適用成功は、理論的導出の厳密な検証となっている。