ブシネスク問題（半無限弾性体の点荷重） — 実践ガイド

以下の手順で進める。

1. 参照解プログラムの作成: PythonでBoussinesq解を計算する関数を用意。10行程度で書ける。この関数自体をNumPyの数値微分で検算しておく

2. 標準メッシュテンプレートの作成: 3水準のバイアスメッシュをGmsh等で自動生成するスクリプトを用意

3. 自動比較スクリプト: ソルバー実行→結果抽出→理論解比較→GCI算出→合否判定を一気通貫で行うPythonスクリプト

4. pytestへの統合: pytest test_boussinesq.py で全チェックが走るようにする

評価点($z = 0.1$ m)での$\sigma_{zz}$と$u_z$が理論値の1%以内。メッシュ「細」での値を基準にGCI < 5%。この基準は片持ち梁ベンチマークと同じレベルだから、一貫性がある。

Hertz接触理論は、Boussinesq解の面圧分布に対する畳み込み積分として導出される。球と平面の接触で、接触面の半楕円面圧分布

を積分すると接触半径$a$と最大面圧$p_0$が得られる。Boussinesq検証→Hertz接触検証→一般接触問題という段階的Validationがベストプラクティスだ。

その通り。V&Vの基本原則は「単純から複雑へ」だ。Boussinesq（線形弾性、軸対称）→Hertz（接触だが解析解あり）→一般接触（非線形、摩擦あり）と進む。各段階で参照解との一致を確認し、次の段階の信頼性基盤とする。

基礎設計の第一近似として使われる。地中の増加応力を計算し、圧密沈下量を推定する。具体的にはNewmarkチャートやFadum積分表を使って矩形分布荷重下の応力増分を求める手法が広く普及している。

FEAとの比較検証にも有用で、特にFEAモデルの境界条件（底面固定深度や側面位置）の妥当性を確認するのにBoussinesq解が参照になる。FEA結果がBoussinesq解から大きく外れている場合、モデル境界が近すぎるか境界条件の設定に問題がある。

作用応力が地盤の降伏応力に比べて十分小さい領域（常時荷重下の地盤）では弾性近似で実用的な精度が得られる。大変形や液状化を扱う場合は弾塑性モデルや有効応力解析に進む必要がある。その場合でもBoussinesq解は初期応力状態の推定やオーダー推定に使える。