加速試験とは｜アレニウス則で寿命を見積もる

この記事でわかること

アレニウス則による加速係数の計算手順（Excel対応）
目標寿命を確認するのに必要な試験時間の求め方
活性化エネルギーを2水準以上の温度から実測する方法
文献値を使うと寿命見積もりが何倍ぶれるか

📌 前提知識：信頼性工学の基礎｜故障率とバスタブ曲線を読んでいると理解しやすくなります

新しい樹脂部品に「使用温度40℃で10年保証」という仕様がつきました。設計審査で「根拠データを出してください」と言われましたが、10年待つわけにはいきません。上司からは「恒温槽で高温試験をやればいいだろう」と言われたものの、何℃で何時間やれば10年分になるのかがわかりません。

加速試験は、実使用より厳しい条件で故障を早め、短い試験時間から実使用での寿命を推定する手法です。温度が関わる劣化なら、アレニウス則という経験則を使って換算できます。

この記事では、125℃で1,000時間の試験が実使用の何年に相当するかを計算し、逆に10年を確認するには何時間必要かを求めます。あわせて、この手法でいちばん危ないところも扱います。

加速試験を使う場面と前提条件
アレニウス則と加速係数
必要な試験時間を求める
⭐ 活性化エネルギーを文献値で決めてはいけない
1. 2水準以上からEaを求める手順
ワイブル解析と組み合わせる
よくある質問
まとめ

加速試験を使う場面と前提条件

加速試験が成り立つのは、次の条件を満たすときです。

劣化の原因が温度で決まる化学反応であるとき。酸化、加水分解、拡散などが該当します
加速しても故障のしかたが変わらないとき。高温にしたせいで別の壊れ方をするなら、その結果は使えません
実使用条件が特定できているとき。使用温度が決まらないと換算先がありません

いちばん重要なのは2つ目です。NG例（△）として、樹脂の耐熱温度を超える温度で試験するのは避けてください。常温では起きない軟化や溶融が始まり、実使用とは別のメカニズムを測ってしまいます。この状態で得た加速係数を使っても、実使用の寿命は当たりません。

温度以外が効く場合は別の式を使います。電圧なら逆べき則、湿度が絡むならアレニウス則と湿度項を組み合わせたモデルが使われます。この記事では温度だけを扱います。

アレニウス則と加速係数

アレニウス則は、化学反応の速度が温度とともに指数関数的に上がることを表した式です。寿命は反応速度の逆数に比例すると考えると、2つの温度での寿命の比が加速係数です。

\[ AF = \exp\left[\frac{E_a}{k}\left(\frac{1}{T_{use}} – \frac{1}{T_{test}}\right)\right] \]

ExcelではEXP関数で書けます。B1に活性化エネルギー、B2に使用温度、B3に試験温度を摂氏で入れる形です。

=EXP((B1/0.00008617)*(1/(B2+273.15) - 1/(B3+273.15)))

\( E_a \) は活性化エネルギー（eV）、\( k \) はボルツマン定数（8.617×10⁻⁵ eV/K）、温度は絶対温度（K）です。摂氏のまま入れると答えが合わないので、273.15を足すのを忘れないでください。

使用温度40℃、試験温度125℃、活性化エネルギー0.7eVで計算します。

\[ \frac{E_a}{k} = \frac{0.7}{8.617 \times 10^{-5}} = 8123.5 \]

\[ \frac{1}{313.15} – \frac{1}{398.15} = 0.00319336 – 0.00251162 = 0.00068174 \]

\[ AF = \exp(8123.5 \times 0.00068174) = \exp(5.538) = 254.2 \]

Excelでは1行で書けます。

=EXP((0.7/0.00008617)*(1/(40+273.15) - 1/(125+273.15)))  → 254.2

125℃での1時間が、40℃での254.2時間に相当するという意味です。

必要な試験時間を求める

加速係数が出れば、試験時間と実使用時間は掛け算と割り算で行き来できます。

まず、125℃で1,000時間の試験を実施した場合。

\[ 1000 \times 254.2 = 254{,}200\ \text{時間} = 29.0\ \text{年} \]

=1000*254.2/(24*365)  → 29.0

1,000時間（約42日）の試験で、実使用29年分を見たことになる計算です。

逆に、目標の10年を確認するのに必要な試験時間を求めます。10年は 10×365×24 = 87,600時間です。

\[ \frac{87600}{254.2} = 344.6\ \text{時間} \approx 14.4\ \text{日} \]

=10*365*24/254.2  → 344.6

約15日の試験で10年保証の根拠が作れる計算です。設計審査に間に合う時間軸に落ちました。

ただしこれは「1個が344.6時間もった」という話ではありません。何個試験して何個壊れたかまで含めて評価する必要があります。必要な試験数の考え方はサンプルサイズの決め方が参考です。

⭐ 活性化エネルギーを文献値で決めてはいけない

ここまでの計算で、答えを決めているのは実質 \( E_a \) の一択です。そして \( E_a \) は材料や劣化モードによって変わります。同じ0.7eVを使い続けてよいのか、確かめてみます。

活性化エネルギー Ea	加速係数 AF	1,000時間試験の相当年数
0.5 eV	52.2	6.0年
0.7 eV	254.2	29.0年
0.9 eV	1,237	141.2年

Eaを0.5と見るか0.9と見るかで、加速係数は52倍から1,237倍まで、23.7倍の開きが出ます。1,000時間の試験結果が「6年分」にも「141年分」にもなってしまう。

文献から拾った値をそのまま使うと、この幅がまるごと見積もりの誤差として残ります。しかも危険なのは、Eaを大きめに見積もると寿命が長く出ることです。合格しやすい方向に誤るので、間違いに気づきにくくなります。

対策は1つで、温度を2水準以上に振って、Eaを自分で実測することです。

2水準以上からEaを求める手順

アレニウス則の式で両辺の対数を取ると、寿命の対数と温度の逆数が直線で結ばれます。

\[ \ln L = \frac{E_a}{k} \cdot \frac{1}{T} + \text{定数} \]

つまり横軸に \( 1/T \)、縦軸に \( \ln L \) を取って直線を引けば、傾きが \( E_a / k \) です。3水準で試験した結果が次のようになったとします。

試験温度	1/T（1/K）	平均寿命 L（時間）
85℃	0.00279213	9,764
105℃	0.00264445	2,942
125℃	0.00251162	1,000

ExcelのSLOPE関数で傾きを求め、ボルツマン定数を掛けます。

=SLOPE(LN(C2:C4), B2:B4)          → 8123.5
=SLOPE(LN(C2:C4), B2:B4)*0.00008617 → 0.700

傾き8,123.5にボルツマン定数を掛けて \( E_a = 0.700 \) eV が得られました。この値なら根拠があります。直線のあてはまりは決定係数R²で確認し、3点が直線に乗らない場合は温度によって劣化メカニズムが変わっている疑いがあります。

回帰の操作手順そのものはエクセル分析ツールを使った回帰分析で解説しています。

ワイブル解析と組み合わせる

加速試験のデータは、平均寿命だけでなくワイブル解析にかけると情報が増えます。

形状パラメータ \( m \) は故障のしかたを表します。初期故障なら1より小さく、偶発故障なら1前後、摩耗故障なら1より大きくなります。ここで見るべきなのは、加速条件と実使用に近い条件で m が揃っているかです。

85℃と125℃で m が大きく違っていれば、それは温度を上げたことで壊れ方が変わった証拠です。前提条件の2つ目が崩れているので、加速係数による換算は使えません。逆に m が揃っていれば、「同じメカニズムのまま速度だけが上がった」という前提を数字で裏づけられます。

設計段階で想定した故障モードと実際の故障モードが一致しているかは、FMEAの故障モード一覧と突き合わせて確認してください。加速試験で出た壊れ方がFMEAに載っていなければ、そのモード自体が見落としだったと判断できます。

よくある質問

Q. 温度は高いほど試験期間が短くて済みますか？

A. 計算上はそうですが、上げすぎると別の壊れ方が始まります。樹脂ならガラス転移温度や軟化点、はんだなら融点の手前が目安で、一般には材料の耐熱温度から余裕を取った温度を選びます。試験期間を1日短くするために前提を壊しては本末転倒です。迷ったら温度を1水準増やし、直線性で前提を確認するほうが確実です。