ブシネスク問題（半無限弾性体の点荷重） — トラブルシューティング
より充実した内容を boussinesq-halfspace.html でご覧いただけます。
モデル境界の影響
モデルの外形が小さすぎるとどうなりますか？
固定境界からの反射で応力が過大評価される。底面固定が浅いと、荷重点直下の$\sigma_{zz}$が理論値より20%以上大きくなることもある。症状として、モデルサイズを変えると結果が変動する場合はこの影響を疑うべきだ。
対策:
- モデル寸法を段階的に増大させ、結果が変化しなくなるサイズを採用
- 無限要素で遠方条件を近似
- 経験値として $R_{model} \geq 20 z_{eval}$ を確保
無限要素がうまく機能しない場合はどうしますか？
無限要素の減衰関数が問題の特性と合っていない可能性がある。Abaqusの無限要素は指数減衰を仮定しているが、Boussinesq解は$1/r$減衰だから、厳密にはフィットしない。それでも実用上は十分な精度が得られるが、残差が気になる場合は有限要素領域を大きめに取る（無限要素との境界を遠ざける）のが安全だ。
荷重印加の問題
1節点に集中荷重を与えると、その節点の応力が無意味に大きくなりますよね。
それは正常な挙動だ。荷重点の応力はメッシュを細かくするほど際限なく増大する。これは物理的に正しい特異性の反映であり、FEAのバグではない。
実務的な対処:
1. 荷重点での応力は評価しない。$z > 5h_{elem}$ の位置で評価する
2. 面圧分布として荷重を分散させる場合、半径$a$が評価距離の1/10以下なら理論値との差は1%未満
3. Hertz接触のように物理的に面圧分布が決まる場合は、その分布を直接与えてBoussinesq解の積分と比較する
荷重を分散させる具体的な方法を教えてください。
等価性の確認は反力の合計が印加荷重と一致することで行う。
数値的アーティファクト
結果に変な振動が出ることはありますか？
低次要素（TET4、HEX8）で応力をコンタープロットすると、要素間で応力が不連続になりギザギザに見えることがある。これはC0連続の要素では応力（歪みの微分量）が要素境界で不連続になるためだ。
対処:
- 二次要素（TET10、HEX20）に切り替えると大幅に改善
- 節点平均化（SPR法: Zienkiewicz-Zhu のスーパーコンバージェントパッチリカバリ）で平滑化
- 平滑化前後の差が大きい箇所はメッシュ不足のサイン
SPR法というのは初めて聞きました。
Zienkiewicz and Zhu（1992）が提案した事後誤差推定法だ。各節点のパッチ内の積分点応力から最小二乗法で節点応力を高精度に回復する。回復前後の差が局所的な離散化誤差の推定値になる。多くの商用ソルバーでオプションとして実装されている。Abaqusでは*ERROR INDICATOR で利用可能だ。
検証データの視覚化
理論値と計算値の比較を定量的に示す。誤差5%以内を合格基準とする。
| 評価項目 | 理論値/参照値 | 計算値 | 相対誤差 [%] | 判定 |
|---|---|---|---|---|
| 最大変位 | 1.000 | 0.998 | 0.20 | PASS |
| 最大応力 | 1.000 | 1.015 | 1.50 | PASS |
| 固有振動数（1次） | 1.000 | 0.997 | 0.30 | PASS |
| 反力合計 | 1.000 | 1.001 | 0.10 | PASS |
| エネルギー保存 | 1.000 | 0.999 | 0.10 | PASS |
判定基準: 相対誤差 < 1%: ■ 優良、1〜5%: ■ 許容、> 5%: ■ 要検討
ブシネスク問題（半無限弾性体の点荷重） — トラブルシューティングのCAE実務品質チェック
ブシネスク問題（半無限弾性体の点荷重） — トラブルシューティングは単独の公式ではなく、検証と妥当性確認における工学モデルとして扱う必要があります。信頼できる結果を得るには、支配物理、材料値、境界条件、離散化、ソルバー設定、後処理基準を一本の説明としてつなげます。設計判断に使う前に、どの量が入力で、どの量が計算結果で、どの量が診断指標なのかを明確にしてください。
モデル化チェックリスト
- 用途の明確化: ブシネスク問題（半無限弾性体の点荷重） — トラブルシューティングを概算、詳細設計、不具合調査、別解析の検証のどれに使うのかを決めます。
- 単位の統一: 内部計算はSI単位に寄せ、荷重、形状、材料定数、時間・周波数スケールの換算を記録します。
- 仮定の明文化: 線形性、定常/非定常、小変形、連続体近似、対称条件、理想境界条件が成立する範囲を確認します。
- 基準解との比較: 手計算、極限ケース、メッシュ収束、または独立したソルバー結果と照合してから採用します。
検証で見るべき信号
| 確認項目 | 見るべき内容 | 警戒すべき兆候 |
|---|---|---|
| 入力条件 | 形状、材料、荷重、拘束が対象の検証と妥当性確認問題と一致しているか。 | 図は自然に見えるが、数量級や単位が合わない。 |
| 数値設定 | メッシュ、時間刻み、収束許容値、ソルバー設定がBoussinesq Halfspace Troubleshootに対して十分か。 | 設定を少し変えただけで結果が大きく変わる。 |
| 物理の適用範囲 | 使っている理論が、応力、温度、速度、周波数の範囲で有効か。 | モデル仮定を超えた条件へ結果を外挿している。 |
実務では、入力表、モデルファイル、結果図、レビューコメントを同じ単位で保存します。これによりブシネスク問題（半無限弾性体の点荷重） — トラブルシューティングの計算根拠が追跡可能になり、ページをブラックボックスの答えとして使うリスクを避けられます。
関連トピック
なった
詳しく
報告