確率論的破壊力学
確率論的破壊力学の理論基礎
確率論的破壊力学とは
先生、確率論的破壊力学って何ですか？
決定論的な破壊力学は「$K \geq K_{IC}$ で破壊」という二値判定。確率論的破壊力学は亀裂サイズ、材料特性、荷重のばらつきを確率変数として扱い、破壊確率を計算する。
「破壊するかどうか」ではなく「破壊する確率は何%か」を評価。
原子力の確率論的リスク評価（PRA）、航空機の損傷許容、パイプラインの信頼性設計で使用。
確率変数
| パラメータ | ばらつきの源泉 |
|---|---|
| 亀裂サイズ $a$ | 検査の不確かさ、初期欠陥分布 |
| $K_{IC}$ | 材料ロット間のばらつき |
| 荷重 $\sigma$ | 運転条件のばらつき |
| Paris定数 $C, m$ | 材料試験のばらつき |
計算手法
まとめ
破壊確率1/1000の意味
確率論的破壊力学では許容破壊確率Pf=10⁻⁶〜10⁻⁴を設定し、欠陥サイズや材料靭性のばらつきを考慮して安全裕度を評価する。原子力圧力容器のIAEA規格ではPf<10⁻⁶/年以下を要求するが、これは「100万容器×1年運転でも1件の破壊が起きない」という厳しい基準だ。Monte Carlo法で10⁷サンプルを使って評価するのが標準的な手順だ。
確率論的破壊力学の数値計算手法
確率論的破壊のFEM
1. FEMで$K$ or $J$を亀裂サイズの関数として計算 — パラメトリックに
2. モンテカルロシミュレーション — 亀裂サイズ、荷重、$K_{IC}$をランダムサンプリング
3. 各サンプルで破壊条件を判定 — $K \geq K_{IC}$?
4. 破壊確率を計算 — 破壊したサンプル数 / 全サンプル数
ツール
まとめ
モンテカルロ法とラテンハイパーキューブ
確率論的破壊解析の数値手法として、ランダムサンプリング（モンテカルロ）と分散削減技術（ラテンハイパーキューブ）がある。モンテカルロは10⁴〜10⁶回の試行が必要なのに対し、ラテンハイパーキューブは同じ精度を10²〜10³回で達成できる。重要サンプリング法と組み合わせれば、低確率破壊（Pf<10⁻⁶）の評価も計算効率よく行える。
確率論的破壊力学の実務適用
確率論的破壊の実務
実務チェックリスト
原子炉圧力容器の確率論的健全性評価
米国NRCはFAVOR（Fracture Analysis of Vessels Oak Ridge）コードを使い原子力圧力容器の確率論的破壊評価を実施している。圧力容器の溶接部に潜在する欠陥サイズをワイブル分布でモデル化し、緊急炉心冷却（ECCS）時の熱衝撃に対するPfを計算する。照射脆化後の10万時間運転での設計寿命評価がこれにより標準化されている。
確率論的破壊力学のソフトウェア比較
確率論的破壊のツール
DARWIN確率論的破壊評価ソフト
SwRI（サウスウエスト研究所）のDARWINは航空機エンジンタービンディスクの確率論的破壊評価専用ソフトだ。FAA/EPRI認定のモンテカルロエンジンを持ち、1ディスクあたり10⁷サンプルの計算を数時間で処理する。GE・P&W・RRの全主要エンジンメーカーがFAA認証プロセスに使用しており、DARWINの計算結果がそのままFAA提出書類の根拠になる。
確率論的破壊力学の先端研究
確率論的破壊の先端
ベイズ更新による破壊確率の修正
超音波探傷試験（UT）の結果を使ってベイズ更新により破壊確率を修正する手法が2010年代に整備された。検査での「検出されなかった」という情報も使え、欠陥サイズ分布の上限を絞り込める。SNC-Lavalinは加圧水炉のPRPD計算にベイズ更新を組み込み、検査間隔を15年から20年に延長する規制申請の根拠データを得た。
確率論的破壊力学のトラブル対応
確率論的破壊のトラブル
確率分布の尾部推定が難しい理由
確率論的破壊評価で最も難しいのは低確率領域（分布の尾部）の推定だ。材料靭性がワイブル分布に従うとしても、尾部のパラメータ推定には大量の実験データ（最低50〜100点）が必要だ。データが少ない段階で極値分布を使った外挿を行うと、設計破壊確率が10倍以上変わることがある。データ数と確率推定の不確実性の関係を常に意識すること。
関連トピック
なった
詳しく
報告