破壊力学 (Fracture Mechanics) — CAE用語解説
破壊力学とは
破壊力学って何ですか？ 普通の強度計算——たとえば「応力が降伏応力以下ならOK」みたいな評価とは何が違うんですか？
いい質問だね。普通の強度計算は、材料が健全——つまりき裂がない前提で、応力が許容値以下かどうかを見る。でも現実の構造物には、溶接の欠陥、疲労き裂、腐食ピットなど、必ず何らかの不連続が存在するんだ。
つまり「き裂がある前提」で評価する学問ということですか？
そのとおり。破壊力学（Fracture Mechanics）は「き裂が既に存在する」ことを出発点にして、そのき裂が安定しているのか、いつ不安定破壊（急速破壊）に至るのかを定量的に評価する学問だ。例えば航空機の胴体に微小き裂が見つかったとき、「次の点検までに危険なサイズまで成長するか？」を判定できる。これが損傷許容設計の基盤だよ。
応力拡大係数 K（SIF: Stress Intensity Factor）
まったく別物だよ。応力集中係数 $K_t$ は無次元の倍率で、ノッチや孔の近くで公称応力の何倍になるかを表す。一方、応力拡大係数 $K$ は単位が $\mathrm{MPa}\sqrt{\mathrm{m}}$ で、き裂先端の応力場の「強さ」そのものを表すパラメータだ。定義式はこうなる：
ここで $\sigma$ は遠方応力、$a$ はき裂長さ、$Y(a/W)$ は形状補正係数だ。例えば無限板の中央き裂なら $Y=1$ になる。
$K_I$ の添字の "I" って何ですか？
それはき裂の変形モードを表している。ちょうどいいタイミングだから、き裂モードの話をしよう。
き裂モード I / II / III
き裂の開き方に種類があるんですか？
うん、3つある。実務でも非常に重要だから覚えておいて。
- モード I（開口モード）：き裂面に垂直な引張でき裂が「パカッ」と開く。実務で圧倒的に多いのがこれ。
- モード II（面内せん断モード）：き裂面に沿って上下の面がスライドする。地震の断層運動をイメージするとわかりやすい。
- モード III（面外せん断モード）：き裂面に沿って前後にずれる、はさみで紙を切るような動き。ねじり荷重で発生する。
それぞれに対応する応力拡大係数 $K_I$, $K_{II}$, $K_{III}$ がある。実際の構造物では複数のモードが混在する「混合モード（Mixed-mode）」が多いんだけど、設計で最も重要なのはモード I だ。破壊靱性値 $K_{Ic}$ もモード I で定義されている。
なるほど。じゃあ $K_I < K_{Ic}$ なら安全、$K_I \geq K_{Ic}$ になったら破壊が起きるということですか？
基本的にはそう。これがLEFM（線形弾性破壊力学）の破壊判定基準だ。ただしこれには大事な前提条件がある。それを理解するために、LEFMとEPFMの違いを話そう。
LEFM vs EPFM
LEFMとEPFMって何の略ですか？ どう使い分けるんですか？
LEFMは Linear Elastic Fracture Mechanics（線形弾性破壊力学）、EPFMは Elastic-Plastic Fracture Mechanics（弾塑性破壊力学） だ。ざっくり言うと、き裂先端の塑性域が小さいか大きいかで使い分ける。
塑性域が小さいっていうのは、具体的にどのくらいですか？
LEFMが有効なのは「小規模降伏（Small-Scale Yielding: SSY）」条件、つまりき裂先端の塑性域の大きさ $r_p$ がき裂長さ $a$ や試験片寸法に比べて十分小さいとき。目安としてはアーウィンの塑性域サイズ推定式がある：
$r_p$ がき裂長さの1/50以下くらいならLEFMで十分。でも延性材料——たとえばアルミ合金やオーステナイト系ステンレス鋼——だと塑性域が大きくなりやすい。そうなると $K$ だけでは評価できなくて、EPFMのパラメータであるJ積分やCTOD（き裂先端開口変位）が必要になるんだ。
じゃあ高強度鋼のような脆い材料はLEFM、柔らかい材料はEPFMということですか？
J積分
J積分って具体的にはどういうものですか？ 数式がとっつきにくくて……
直感的に言うと、き裂が単位面積だけ進展するときに解放されるエネルギーだ。Rice（1968）が定義した経路積分で、き裂先端を囲む任意の経路 $\Gamma$ に沿って計算する：
ここで $W$ はひずみエネルギー密度、$\mathbf{T}$ は経路上の表面力ベクトル、$\mathbf{u}$ は変位ベクトルだ。
経路をどこに取っても同じ値になるんですか？
そこがJ積分のすごいところで、経路独立（path-independent）なんだ。き裂先端のすぐ近くで計算しても、遠くの経路で計算しても同じ値になる。だからFEMでは、き裂先端の特異場の精度が多少悪くても、離れた領域で精度よくJ値を計算できる。AbaqusやANSYSでは「contour integral」として複数の経路でJを自動計算して、収束を確認できるようになっているよ。
J積分とKの関係ってあるんですか？
小規模降伏条件では、こういう関係がある：
つまりLEFMが適用できる範囲では $J$ と $K_I$ は等価だ。でも大規模降伏になると $K$ の概念は破綻して、$J$ だけが物理的に意味を持つ。だからEPFMでは $J_{Ic}$（臨界J値）が破壊靱性の指標になるんだ。
XFEM（拡張有限要素法）
従来のFEMでき裂を扱う場合、き裂面にメッシュの要素境界を一致させる必要がある。き裂先端には $1/\sqrt{r}$ の応力特異性があるから、特異要素（中間節点を1/4点に移動させたもの）を配置するのが定番だ。でも、き裂が進展するたびにリメッシュしなきゃいけない。3Dだとこれが本当に大変でね。
XFEMだとリメッシュが不要になるんですか？
そう、それが最大のメリット。XFEMでは通常のFEMの形状関数に、濃縮関数（enrichment function）を追加する。具体的にはこういう形だ：
第1項は通常のFEM、第2項はHeaviside関数 $H$ でき裂面の不連続を表現、第3項はき裂先端の特異性を表す枝関数 $F_\alpha$ だ。メッシュをき裂に合わせる必要がないから、き裂の位置や形状はレベルセット関数で自由に定義できる。
実務ではどのソルバーで使えるんですか？
Abaqus 6.9以降でXFEMが標準サポートされていて、き裂の挿入・進展が比較的簡単に設定できる。ANSYS MechanicalにもSmart Crack Growthという機能がある。オープンソースだとCode_AsterがX-FEMをサポートしている。ただしXFEMは万能じゃなくて、3Dの複雑なき裂分岐や接触を伴うケースではまだ課題が多い。そういう場合はCohesive Zone Model（CZM）やPhase-field法も検討するといい。
CAE実務での破壊力学
CAEで破壊力学の解析をするとき、一番注意すべきポイントは何ですか？
いくつかあるけど、最も重要なのはき裂先端のメッシュだね。従来FEMなら、き裂先端に同心円状の「スパイダーウェブメッシュ」を作って、先端節点を1点に集約した特異要素を使う。最内周の要素サイズはき裂長さの 1/20 以下が目安だ。J積分のコンター数は5〜10本取って、外側のコンターで値が収束していることを確認する。
疲労き裂進展の解析はどうやるんですか？
まず各荷重サイクルでの応力拡大係数範囲 $\Delta K = K_{\max} - K_{\min}$ を計算して、Paris則でき裂進展速度を求める：
$C$ と $m$ は材料定数で、鋼だとだいたい $m \approx 3$、アルミで $m \approx 3\text{--}4$ が典型値。これを数値積分すれば、初期き裂長さ $a_0$ から臨界長さ $a_c$（$K_{I\max} = K_{Ic}$ となる長さ）まで何サイクルかかるか——つまり残余寿命が計算できる。Abaqusの *CONTOUR INTEGRAL でKを求めて、外部スクリプトやfe-safeのようなツールで進展解析するのが実務の流れだ。
LEFM、EPFM、XFEM、CZM……選択肢が多くて迷いそうです。どう選べばいいんですか？
まとめるとこんな感じだ：
- 高強度鋼・脆性材料、き裂長さ既知 → LEFM（$K$ 評価）が最もシンプルで信頼性も高い
- 延性材料・大規模降伏 → EPFM（$J$ 積分 or CTOD）
- き裂の進展経路が未知・3D → XFEM または Phase-field
- 界面はく離・接着・複合材 → CZM + VCCT
最初の判断基準は「塑性域の大きさ」と「き裂経路が既知かどうか」だ。迷ったらまずLEFMで評価して、$r_p$ が大きすぎたらJ積分に移行する、という段階的なアプローチがいいよ。
CAE用語の正確な理解は、チーム内のコミュニケーションの基盤です。 — Project NovaSolverは実務者の学習支援も視野に入れています。
破壊力学の実務で感じる課題を教えてください
Project NovaSolverは、CAEエンジニアが日々直面する課題——セットアップの煩雑さ、計算コスト、結果の解釈——の解決を目指しています。あなたの実務経験が、より良いツール開発の原動力になります。
お問い合わせ（準備中）関連トピック
なった
詳しく
報告