構造信頼性解析 戻る
構造解析

構造信頼性解析・破損確率計算

FORM(一次信頼性法)とモンテカルロ法で信頼性指標β・破損確率Pf を計算。抵抗RとS荷重の分布干渉を可視化。

パラメータ設定
抵抗 R(強度)
分布型
平均 μ_R
変動係数 CoV_R
荷重 S (Load)
分布型
平均 μ_S
変動係数 CoV_S
目標信頼性指標 β_T
計算結果
信頼性指標 β
破損確率 P_f
安全率 γ = μ_R/μ_S
必要 μ_R (目標β)
分布
理論・主要公式

正規分布(FORM):

$$\beta = \frac{\mu_R - \mu_S}{\sqrt{\sigma_R^2 + \sigma_S^2}}, \quad P_f = \Phi(-\beta)$$

対数正規分布: $\zeta = \sqrt{\ln(1+\text{CoV}^2)}$, $\lambda = \ln\mu - \zeta^2/2$

$$\beta_{LN}= \frac{\lambda_R - \lambda_S}{\sqrt{\zeta_R^2 + \zeta_S^2}}$$

目標β_T に対する必要平均強度:

$$\mu_R^* = \mu_S + \beta_T\sqrt{\sigma_R^2 + \sigma_S^2}$$

構造信頼性解析・破損確率計算とは

🙋
「構造信頼性解析」って何ですか? 普通の安全率計算と何が違うんですか?
🎓
大まかに言うと、「ばらつき」を考慮した確率論的な安全評価だね。普通の安全率は「平均強度÷平均荷重」だけど、現実には材料強度も荷重もバラつくよね。例えば、同じ工場で作った鋼材でも強度は少しずつ違うし、建物にかかる積雪荷重も年によって違う。このツールでは、強度と荷重の分布がどれだけ重なるか(干渉するか)を可視化して、破損する確率そのものを計算するんだ。上のスライダーで「平均強度μ_R」と「平均荷重μ_S」を動かしてみると、二つの分布が離れたり近づいたりするのがわかるよ。
🙋
え、そうなんですか! で、出てくる「信頼性指標β」って、数字が大きいほど安全ということですか? 3.0とか4.0って、どんな意味があるんですか?
🎓
その通り! βは安全の度合いを表す無次元の指標で、大きいほど破損確率は低くなる。実務では、β=3.8なら破損確率Pfは約0.007%、β=4.2なら約0.0013%といった感じだ。建築基準では目標をβ=3.5以上とすることが多いし、航空機のエンジン部品などはもっと厳しいβ=6以上を要求されることもあるよ。シミュレーター右側の「目標信頼性指標β_T」を4.0に設定して「計算」ボタンを押すと、その目標を満たすために必要な平均強度が逆算される機能もあるんだ。
🙋
なるほど! で、「FORM」と「モンテカルロ法」ってどっちを使えばいいんですか? パラメータの「分布型」を「対数正規分布」に変えると結果が結構変わりますね。
🎓
良いところに気づいたね! FORM(一次信頼性法)は計算が速く、設計の初期段階でパラメータをいじりながら感度分析するのに向いてる。ただ、強度や荷重が対数正規分布のように正規分布から外れると、FORMだけだと誤差が出る場合がある。そこで、分布型を「対数正規」に変えて、「サンプル数N」を増やして「モンテカルロ法」で計算してみて。乱数シミュレーションで真の破損確率に近づく様子が体感できるよ。現場では、FORMでおおよその目安をつけて、最終確認でモンテカルロ法を使うことが多いんだ。

よくある質問

Rは材料試験データや設計基準値から平均と標準偏差を設定します。Sは想定する荷重条件(静的・動的)や過去の計測データに基づき設定してください。正規分布を仮定する場合、負の値が出ないよう注意が必要です。
FORMは分布が正規分布で限界状態関数が線形の場合に正確ですが、非線形性や非正規分布では誤差が生じます。モンテカルロ法はサンプル数が少ないとばらつきが大きく、収束には数万回以上の試行が必要です。両者の差異は分布形状や非線形性の影響を反映しています。
はい、あります。βが負の値は、荷重Sの平均が抵抗Rの平均より大きい場合に生じ、破損確率Pfが0.5を超えることを意味します。この場合、設計の見直し(材料変更や断面積増加など)が必要です。
材料強度や荷重は正の値のみをとり、分布が右に歪む特性があります。対数正規分布を用いると負の値を排除でき、実現象に近いモデル化が可能です。特に疲労強度や風荷重など、変動係数が大きいデータに適しています。

実世界での応用

建築・土木構造物のコード適合性確認: 日本の建築基準法の荷重指針や国際規格ISO 2394では、目標信頼性指標が定められています。FEM解析で求めた応力分布の統計量と材料強度の統計データを本ツールに入力し、コードが要求するβ(例:通常時β≥3.5、地震時β≥2.9など)を満たしているか確認する設計フローに活用されます。

航空宇宙・自動車部品の確率設計: 軽量化と高信頼性が求められるエンジン部品やシャシーでは、従来の安全率だけでは不十分です。疲労強度や破壊靭性のばらつき、運用荷重の変動を考慮し、部品ごとに許容破損確率(例:10⁻⁷)を設定。モンテカルロ法を用いてその目標を満たす材料仕様や寸法を決定します。

CAE/FEM解析との連携(確率的FEM): 有限要素法の結果(応力、ひずみ、変形)は入力条件(材料定数、荷重、境界条件)のばらつきに依存します。これらの入力パラメータに確率分布を設定し、多数のFEM実行結果から応力の分布を求め、本ツールで強度分布と照合することで、構造全体のシステム信頼性を評価できます。

既存構造物の残存寿命評価と保全計画: 橋梁やプラントなど老朽化が進むインフラでは、材料の劣化(CoVの増大)や荷重条件の変化(μ_Sの増加)をモデルに反映させ、現在の信頼性指標βを時系列で追跡します。これに基づいて、補修や交換の優先順位付け、経済的な保全計画を立案します。

よくある誤解と注意点

このツールを使い始めるとき、いくつか陥りがちな落とし穴があるから気をつけてね。まず「平均値と変動係数(CoV)の設定根拠」。例えば「鋼材の降伏強度」を入力するとき、カタログ値の「最小保証値」をそのまま平均値μ_Rに入れていませんか? それは誤りだ。平均値は、実際の製造バラツキを反映したデータから推定する必要がある。カタログ値が490MPaなら、実際の平均はそれより高い520MPaくらいかもしれない。この設定を間違えると、破損確率が実際より大幅に高く(または低く)出てしまう。実務では材料メーカーにバラツキデータを問い合わせるのが第一歩だ。

次に「分布型の選択の影響を過小評価」。とりあえず正規分布で計算して「β=4.0だから大丈夫」と安心するのは危険。特に強度は正の値しか取らず、低い値側に裾を引く(左に歪んだ)分布になりがち。例えばコンクリートの圧縮強度は対数正規分布に近い。この場合、同じ平均・標準偏差でも、正規分布仮定より破損確率Pfが数倍高くなることもある。重要な検討では、必ず背景データのヒストグラムを確認し、適切な分布型を選ぼう。

最後は「モンテカルロ法のサンプル数Nの盲信」。N=10万で計算すればいつも正確、と思いがちだが、破損確率が極めて小さい(例えば10^-6以下)場合、単純なモンテカルロ法では数億回のサンプリングでも結果が安定しない。こんな時は「重点サンプリング法」などの高度な手法が必要になる。このツールはあくまで概念理解と初期検討用。確率が極端に小さい領域の精密計算は、専門のソフトウェアに任せるのが現実的だ。