J積分（弾塑性破壊力学）
J積分（弾塑性破壊力学）の理論基礎
J積分とは
先生、J積分って何ですか？
$\Gamma$ は亀裂先端を囲む任意の経路。$W$ はひずみエネルギー密度、$\mathbf{T}$ は牽引力。
経路に依存しない！ だから「積分経路をどこに取っても同じ値」。
弾性体では経路独立性が厳密に成り立つ。弾塑性でも除荷がなければ（単調載荷）ほぼ経路独立。FEMでは亀裂先端を囲む複数のcontour（輪郭）でJを計算し、収束を確認する。
$J$ と $K$ の関係
線形弾性では：
$E' = E$（平面応力）、$E' = E/(1-\nu^2)$（平面ひずみ）。$J$ は $K$ の2乗に比例。
破壊条件
$J_{Ic}$ は臨界J積分値（材料特性）。ASTM E1820で試験方法が規定。
まとめ
要点：
- $J$ = 亀裂先端のエネルギー解放率 — 弾塑性に対応
- 経路独立 — 亀裂先端を囲む任意の経路で同じ値
- $J = K^2/E'$ — 線形弾性との関係
- $J \geq J_{Ic}$ で破壊 — ASTM E1820で$J_{Ic}$を測定
- FEMで*CONTOUR INTEGRAL — 亀裂先端のJを自動計算
Riceが9ページの論文で世界を変えた
J積分はJames Rice（ハーバード大）が1968年のJAppl Mech誌に発表した9ページの論文で提案した。き裂先端を囲む任意の積分経路で同じ値になる「経路独立性」が最大の特徴で、弾塑性条件でもエネルギー解放率を定義できる。この発見は破壊力学の弾塑性への拡張を可能にし、Riceは1983年のTimoshenko賞を受賞した。
J積分（弾塑性破壊力学）の数値計算手法
J積分のFEM
```
*CONTOUR INTEGRAL, CONTOURS=5, TYPE=J
crack_tip_node, direction_vector
```
5つの輪郭（contour）でJを計算。輪郭が亀裂先端から遠ざかるにつれて値が収束するはず。
5つの輪郭で値が収束しなかったら？
メッシュが粗い or 塑性域が大きい。メッシュを細分化するか、輪郭数を増やす。外側3〜4個の輪郭の値がほぼ一致すれば収束OK。
亀裂先端のメッシュ
亀裂先端には集中メッシュ（Spider web mesh）を配置。中心の1点（亀裂先端）から放射状に要素を配置。
- 2次要素（C3D20R）推奨 — 亀裂先端の特異性を正確に捕捉
- Quarter-Point要素 — 亀裂先端で中間節点を1/4に移動。$1/\sqrt{r}$特異場
まとめ
FEMでのJ積分計算：仮想き裂進展法
FEMでのJ積分計算には仮想き裂進展法（Domainアーカイブ法）が精度・効率ともに優れている。き裂先端から3〜5要素離れた領域でのみ積分し、FEMの特異要素（collapsed quarter-point element）を使う必要がない。ANSYSのFRACTURE TBコマンドは内部でこの方法を使い、経路1〜10の平均値を自動出力するため収束確認が容易だ。
J積分（弾塑性破壊力学）の実務適用
J積分の実務
圧力容器の亀裂評価（API 579 FFS-1）、パイプラインの欠陥評価、原子力の破壊力学評価（R6法）で使用。
ASTM E1820試験
J-R曲線（J vs. 亀裂進展量$\Delta a$）の試験。CT（Compact Tension）試験片で実施。$J_{Ic}$（亀裂開始の臨界値）と$J-R$曲線を取得。
実務チェックリスト
圧力容器ノズル部の弾塑性破壊評価
ASME Sec.XI Code では圧力容器ノズル隅肉部の亀裂評価にJ積分を使う。原子力級配管のSA-508 Cl.3鋼のJ-R曲線（J vs Δa）を入力に、初期き裂が不安定成長に転じる条件（Ji=安定成長開始点）を計算する。60年延長運転を検討する際に、照射脆化後のJIcが保守的に50%低下しても安全余裕があることを証明する解析が必要だ。
J積分（弾塑性破壊力学）のソフトウェア比較
J積分のツール
選定ガイド
Abaqus J積分計算の精度検証
AbaqusのJ積分計算精度はNIST標準問題での検証が公開されており、単辺ノッチ引張（SENT）試験体での誤差は0.5%以内とされる。IAEA原子力安全規制では各国の規制当局がAbaqus計算結果を受け入れる際にNISTベンチマーク結果の提出を要求する場合がある。最外積分経路（5〜10要素）の使用と粗いメッシュ回避が精度確保のポイントだ。
J積分（弾塑性破壊力学）の先端研究
J積分の先端
3次元J積分と面外変位の補正
工業構造物では板厚方向にJ積分値が変化し、中心部（平面ひずみ）で最大、表面（平面応力）で最小になる。3DFEMによる「修正J積分」では面外変位の貢献項を追加した「Js値」を計算し、薄板（t<4mm）への適用精度を向上させる。空力荷重を受ける航空機外板の亀裂評価では2DとJsの差が20〜40%に達することがある。
J積分（弾塑性破壊力学）のトラブル対応
J積分のトラブル
積分経路依存が生じる原因
J積分の経路独立性が成り立たない場合、主な原因は体積力の無視と弾性除荷の存在だ。熱応力がある場合は体積力項を追加した「Jth」を計算する必要があり、FABquis VFLUX法などで対応できる。また繰り返し塑性変形が生じるサイクル荷重下では単調負荷の前提が破れるため、ΔJやCyclic J積分概念を使う必要がある。
関連トピック
なった
詳しく
報告