応力拡大係数の求め方・き裂モードと破壊じん性の基本

応力拡大係数の求め方・き裂モードと破壊じん性の基本

高張力鋼を使うほど、き裂に対して材料が敏感になり破壊リスクが上がります。

応力拡大係数の求め方における基本式とモードの種類

応力拡大係数（英語：stress intensity factor、略称：K）とは、き裂や欠陥が存在する金属材料において、き裂先端近傍に形成される応力場の強さを定量的に表すパラメータです。単位はSI単位系で「MPa·m^(1/2)」であり、現場の技術文書では「MPa√m」と表記されることが多い値です。

基本式はシンプルです。それが以下の形になります。

$$K = F \cdot \sigma \cdot \sqrt{\pi a}$$

ここで各記号の意味は次のとおりです。

• K：応力拡大係数（MPa√m）
• F：き裂形状や部材形状・負荷条件による無次元の補正係数
• σ：き裂がない場合を想定した公称応力（MPa）
• a：き裂の長さ（m）。中央き裂の場合は全長2aの半分を使用

最も基本となるのは、端のない無限平板（無限平面）に長さ2aの貫通き裂がある場合で、このときF＝1.0となり、式は $$K_I = \sigma\sqrt{\pi a}$$ と簡略化されます。現実の部材には「端」があり形状も複雑なため、この単純モデルをそのまま適用できるケースはほとんどありません。それが基本です。

き裂の変形様式（モード）は3種類に分類されます。これはき裂に対する力の作用方向によって区別されます。

• 🔺 モードⅠ（開口型）：き裂面に垂直な引張応力によってき裂が開口する様式。最も発生頻度が高く、破壊力学の計算で主に扱われるのがこのモードです。
• 🔺 モードⅡ（面内せん断型）：き裂面に平行な方向にすべり（面内せん断応力）が作用する様式。
• 🔺 モードⅢ（面外せん断型）：き裂面に垂直な方向にねじれるようなせん断（面外せん断）が作用する様式。

実務上の強度評価では、モードⅠの応力拡大係数「K_I」を求める場面が圧倒的に多いです。つまりモードⅠが原則です。ただし、複合荷重が作用する溶接構造やシャフト類では、複数モードが重なるケースもあるため、設計段階での負荷条件の整理が不可欠です。

応力拡大係数の計算で単位の取り扱いには注意が必要です。き裂長さ「a」はメートルm単位、応力「σ」はパスカルPaではなくMPa単位を使う場合は計算経路を正しく設定する必要があります。長さをメートルm、応力をMPa単位で計算すると、結果のKの単位はMPa√mになります。単位は要確認です。

参考：線形破壊力学と応力拡大係数の基礎理論について体系的に解説されています。

CAEと強度計算の森 — 線形破壊力学と応力拡大係数

応力拡大係数の求め方で使う補正係数Fとハンドブックの活用法

現実の部材に応力拡大係数を適用するとき、最初にぶつかる壁が「補正係数F」です。Fは、き裂の幾何形状・部材の寸法・負荷方式・境界条件などが複雑に絡み合って決まる無次元数で、世界中の研究者が解析や実験を通じてデータを蓄積してきました。これが実用上の核心です。

補正係数Fのデータが集積された書籍として最も権威があるのが「応力拡大係数ハンドブック（Stress Intensity Factors Handbook）」です。ASTM（米国材料試験協会）が発刊したもので、多くの工学系大学図書館に所蔵されています。日本語では、日本材料学会が発行する「応力拡大係数ハンドブック Vol.1 & Vol.2」電子版が広く活用されています。

• 📚 ASTM版ハンドブック：世界標準として参照される英語資料。多くの試験片形状・き裂パターンの補正係数Fが掲載。
• 📚 日本材料学会版（JSMS）：試験片名・材質・負荷形態・き裂形態などで検索できる電子版あり。国内向けに使いやすい構成。

具体的な使い方の手順は次のようになります。

1. 評価したい部材のき裂形状・負荷条件を整理する
2. ハンドブックの中から近い形状のケースを探す
3. 補正係数Fの式またはグラフ値を読み取る
4. 公称応力σとき裂長さaをメートルmで代入し、Kを計算する

代表的なケースをいくつか示します。まず「半無限板の片側き裂（表面き裂）」の場合、補正係数はF＝1.12が基本値として使われます。「有限幅板の中央き裂」の場合は、き裂長さ(2a)と板幅(W)の比がFの値を左右します。き裂が板幅に対して小さい（a/W ≪ 1）うちはF≒1.0ですが、き裂が板幅の半分程度（a/W ≈ 0.4）に達すると、Fは1.2〜1.3程度まで増加します。これは使えそうです。

ハンドブックにズバリ一致する形状が見つからない場合も多いです。そのような場合の実務的な対処法として、現在では有限要素法（FEM）を使った数値計算が一般的になっています。FEMを使った応力拡大係数の求め方については後の節で詳しく取り上げます。

参考：補正係数と応力拡大係数ハンドブックの具体的な計算方法について詳しく解説されています。

イプロス — 破壊工学の基礎知識7「応力拡大係数の計算方法」

参考：日本材料学会の応力拡大係数ハンドブック電子版の情報はこちらです。

日本材料学会 — 応力拡大係数ハンドブック Vol.1&Vol.2 電子版

応力拡大係数の求め方：FEM（有限要素法）による直接変位法

ハンドブックに一致する形状が見つからない場合や、複雑な応力場を持つ実機部品を評価する場面では、有限要素法（FEM）による応力拡大係数の算出が有効です。FEMで応力拡大係数を求める手法には大きく「直接法」「エネルギ法」「重ね合わせ法」の3系統があります。

• ⚙️ 直接変位法（直接法）：き裂開口部の変位をFEM解析で求め、理論式に代入してKを算出する方法。汎用FEMソフトでも実施可能。
• ⚙️ J積分法（エネルギ法）：き裂先端まわりの閉じた積分経路でエネルギを評価してKを求める方法。高精度かつ安定していて、多くのハイエンドFEMソフトに実装されています。
• ⚙️ VCE法（仮想き裂閉口法）：き裂を微小量だけ閉じる仮想操作でエネルギ解放率Gを求め、KとGの関係式でKに換算する方法。

ここでは現場でも実践しやすい「直接変位法」の手順を具体的に説明します。モードⅠの応力拡大係数K_Iを対象とします。

き裂近傍のy方向変位v(r)は理論的に次の式で表されます。

$$v(r) = \frac{K_I}{E} \cdot \frac{4(1-\nu^2)}{\sqrt{2\pi}} \cdot \sqrt{r}$$

ここで、rはき裂先端からの距離、Eはヤング率、νはポアソン比です。

この式を変形すると。

$$K_I = \frac{E}{4(1-\nu^2)} \cdot \sqrt{\frac{2\pi}{r}} \cdot v(r)$$

という形になり、FEM解析で得られたき裂面上の節点の座標rと変位v(r)を代入すれば、各節点でのK_I推定値が得られます。

実際の作業では、き裂面上のn個の節点についてそれぞれK_I推定値を計算し、横軸をr、縦軸をK_I推定値としてExcelで散布図を作成します。その後、1次式で近似したカーブフィットの「縦軸切片（r=0でのK_I値）」が最終的なK_Iとなります。FEMによる実測値と理論解の差が1%未満に収まるケースも多く、計算精度は十分実用的です。

FEM解析でき裂をモデリングする際の重要なポイントが2点あります。1点目は「き裂先端近傍の要素を細かく分割すること」で、理論的にはき裂先端付近のr≦a/10（き裂長さの1/10以内）の範囲を対象とすることが推奨されています。2点目は「き裂面を接触しないよう設定すること」です。き裂面を「接合」したままではき裂開口変位が正しく取得できません。要素分割が命です。

ハイエンドのFEMソフト（ANSYS、ABAQUS、NX Nastranなど）にはJ積分やVCE法によるK・G・J値を自動算出する機能が搭載されています。業務で使用しているFEMソフトがあれば、まずその機能を確認することをおすすめします。

参考：有限要素法を使った応力拡大係数の直接変位法と計算例が詳しく掲載されています。

CAEと強度計算の森 — 応力拡大係数の求め方（FEM直接変位法）

応力拡大係数と破壊じん性値・き裂進展判定の考え方

応力拡大係数Kを求めた後に、「では実際にそのき裂は危険なのか？」を判断するための指標が「破壊じん性値（Fracture Toughness）」です。破壊じん性値とは、材料がき裂の進展に対して持つ抵抗力を示す材料特性値で、記号はKc（または平面ひずみ条件下の厳密値としてKIc）で表されます。単位はKと同じくMPa√mです。

破壊の判定条件はシンプルです。

$$K \geq K_c \Rightarrow \text{き裂が進展して破壊}$$

応力拡大係数Kが材料のKcを超えた瞬間、き裂は急速に進展して一発破壊に至ります。KがKcを下回っている間は、静的には安全と判断できます。結論は「KとKcの大小比較」です。

KIcは試験によって求める材料定数です。代表的な試験規格として、ASTM E399（平面ひずみ破壊じん性試験）、JIS G 0564、ISO 12135などがあります。試験片に疲労予き裂を導入し、荷重を段階的に増加させながらき裂開口変位をクリップゲージで測定して算出します。

代表的な金属材料のKIc目安値（平面ひずみ条件）を参考として示します。

材料	引張強さの目安	KIcの目安
軟鋼（SS400）	約400 MPa	140〜200 MPa√m
高張力鋼（HT780）	約780 MPa	50〜100 MPa√m
アルミニウム合金（A7075）	約500 MPa	25〜30 MPa√m
チタン合金（Ti-6Al-4V）	約900 MPa	50〜80 MPa√m

ここに注目すべき事実があります。高張力鋼（HT780）は軟鋼（SS400）に比べて約2倍の引張強さを持ちますが、KIcは軟鋼の約1/3〜1/2程度にとどまります。つまり、高張力鋼は「強い」のに「き裂に対してはるかに脆い」という特性を持つわけです。高強度材料を採用したからといって、必ずしもき裂に対する安全性が向上するとは限りません。高強度材ほどKIcに注意が必要です。

この事実は現場において深刻な意味を持ちます。高張力鋼を使った軽量化設計では、同時に微細欠陥や加工傷に対するき裂進展しやすさも増大します。研磨傷・溶接止端部・孔端部などの加工プロセスに起因する微細欠陥の管理が、高張力鋼使用設計では特に重要になります。

参考：破壊靭性試験の試験規格ASTM・JISの概要と試験方法が説明されています。

日鉄テクノロジー — 靭性・破壊靭性試験の概要

参考：応力拡大係数と破壊じん性の関係がリバティ船の事故事例とともに解説されています。

イプロス — 破壊工学の基礎知識5「応力拡大係数と破壊じん性」

パリス則を使った応力拡大係数による疲労き裂進展寿命の推定

静的破壊の評価にK vs Kcの比較が有効であるのと同様に、繰り返し荷重下での疲労き裂の進展評価には「応力拡大係数範囲ΔK（デルタK）」を用います。ΔKは、変動荷重の1サイクルにおける応力拡大係数の最大値Kmax と最小値Kmin の差です。

$$\Delta K = K_{max} - K_{min} = F \cdot \Delta\sigma \cdot \sqrt{\pi a}$$

ΔKとき裂の1サイクルあたりの進展量（き裂進展速度）da/dn \μm/cycle の関係は、多くの金属材料で「パリス則（Paris' law）」と呼ばれるべき乗関係が成立することが知られています。

$$\frac{da}{dn} = C \cdot (\Delta K)^m$$

CとmはΔK–da/dn試験（疲労き裂進展試験）によって決定される材料定数で、鉄鋼材料ではm＝2〜4程度の値になることが多いです。これが基本です。

パリス則が成立するΔK中間領域（安定進展領域）では、き裂初期長さa1から最終き裂長さa2まで成長するのに必要な繰り返し数（疲労き裂進展寿命N）を積分によって推定できます。

$$N = \int_{a_1}^{a_2} \frac{da}{C \cdot (F\Delta\sigma\sqrt{\pi a})^m}$$

Fが一定とみなせる小さいき裂範囲では、この積分の解が得られます。m≠2のとき。

$$N = \frac{2}{(m-2) \cdot C \cdot (F\Delta\sigma\sqrt{\pi})^m} \left(a_1^{1-m/2} - a_2^{1-m/2}\right)$$

m=2のときは。

$$N = \frac{1}{C \cdot (F\Delta\sigma)^2 \cdot \pi} \cdot \ln\left(\frac{a_2}{a_1}\right)$$

これを活用することで、非破壊検査で発見された欠陥（例えば溶接部の0.5mmの表面き裂）が、使用条件下で何サイクル後に臨界き裂長さに達するかを定量的に予測できます。これは使えそうです。

パリス則を活用するために必要な入力データをまとめると次のようになります。

必要なデータ	取得方法
材料定数 C・m	ASTM E647準拠の疲労き裂進展試験、または文献データベース
応力拡大係数範囲ΔK	ハンドブック公式またはFEM解析で算出
初期き裂長さ a1	超音波・磁粉・浸透探傷などの非破壊検査で測定
最終き裂長さ a2（臨界値）	破壊靭性値KIcから K=KIc となるき裂長さを逆算

また、ΔKが非常に小さい領域では「下限界応力拡大係数範囲ΔKth」という重要な指標があります。ΔK ≦ ΔKthであれば疲労き裂はほぼ進展しないため、安全側に設計する際の判断基準になります。鋼材（オーステナイト系ステンレスを含む）の高応力比（R≧0.5）条件でのΔKthは約2.0 MPa√m、アルミニウム合金では約0.7 MPa√mが設計時の目安値として広く使われています。ΔKthが条件です。

非破壊検査の結果と応力拡大係数の計算を組み合わせることで、「今ある欠陥があとどれくらい使えるか」を数値で評価できます。これを「損傷許容設計（Damage Tolerant Design）」と呼びます。航空機・圧力容器・橋梁などでは義務化されており、金属加工の現場でも検査頻度や補修判断の根拠として活用が広がっています。

参考：パリス則・下限界ΔKthの詳細な解説と設計への応用法が掲載されています。

CAEと強度計算の森 — き裂進展速度・パリス則・下限界ΔKth

参考：疲労き裂進展速度と応力拡大係数範囲の関係についての技術資料です。

ねじ締結技術ナビ — 応力拡大係数と疲労寿命（パリス則）

応力拡大係数の求め方で見落とされがちな「残留応力」の影響

応力拡大係数の計算で、多くの現場技術者が見落としがちな要因が「残留応力」です。この視点は教科書的な解説にはほとんど登場しませんが、実際の金属加工品・溶接構造物の強度評価では無視できない重要因子です。意外ですね。

残留応力とは、外力が取り除かれた状態でも部材の内部に残存する応力のことです。溶接・プレス加工・旋削・研削・熱処理などほぼすべての金属加工プロセスで、多かれ少なかれ残留応力が発生します。残留応力は必ず発生します。

残留応力がき裂に与える影響は、応力拡大係数の「見かけの値」を変化させることで現れます。具体的には次のように考えます。

$$K_{total} = K_{applied} + K_{residual}$$

引張側の残留応力（＋）がき裂先端に存在すると、K_totalがK_appliedより大きくなり、き裂は外部荷重だけの評価より速く進展します。逆に圧縮側の残留応力（−）が存在すると、K_totalを低下させ、き裂進展を抑制する効果があります。ショットピーニング処理や表面圧延が疲労強度改善に有効なのはこの原理によるものです。

溶接構造物を例に挙げると、溶接止端部には高い引張残留応力（材料の降伏応力に近い値になる場合がある）が発生しているケースがほとんどです。この残留応力を無視してΔKを計算すると、き裂進展寿命を実態より大幅に楽観的に評価してしまう危険があります。実態より危険側かもしれません。

残留応力の影響を考慮した応力拡大係数の評価方法は次のとおりです。

• 🔧 残留応力測定：X線回折法・中性子回折法・磁気応力測定法など。X線法は表面の残留応力（深さ数μm〜数十μm）を高精度で測定可能。
• 🔧 FEM解析（加工・熱解析連成）：溶接や加工プロセスをFEMでシミュレーションし、残留応力分布を予測した上でき裂評価に組み込む方法。
• 🔧 FITNETやAPI 579などの評価基準の活用：残留応力を含む応力拡大係数の計算手順を規定した国際規格。原子力・石油化学プラントでは必須の評価文書です。

特に溶接後に応力除去焼なまし（SR処理）を実施することで、残留応力を大幅に低減できます。SR処理は「き裂進展の促進因子を除去する」という観点でも合理的な処置です。JIS Z 3700シリーズなどの溶接構造物の疲労設計基準では、残留応力状態の違いによるΔK補正が規定されています。残留応力まで見て初めて完全な評価です。

ここでのポイントを整理すると、「外部荷重から求めたΔKだけで疲労寿命を評価するのは、加工品・溶接品では不十分である」ということです。特に溶接止端部・孔加工端部・曲げ加工R部などの高残留応力が発生しやすい箇所では、残留応力を把握した上で応力拡大係数を評価することが、現場での破損防止に直結します。残留応力の確認が条件です。

参考：き裂の進展と応力拡大係数・疲労き裂進展に影響する因子についての詳細解説です。

エンジニア教育.com — き裂の進展と応力拡大係数について丁寧に解説