CTOD（Crack Tip Opening Displacement）破壊力学
理論と物理
概要
先生、破壊力学でKIcとかKIIIcって応力拡大係数は学んだんですが、CTODって何が違うんですか？
線形弾性破壊力学（LEFM）の $K_{Ic}$ は亀裂先端近傍が「完全弾性体」という仮定に基づいている。でも鋼材のように靭性が高い材料では亀裂先端で大きな塑性変形が起きる——この「塑性域」が小さければLEFMが使えるが、大きくなると $K_{Ic}$ だけでは記述できない。CTODは亀裂先端がどれだけ「開いた（変位した）」かを直接計測する弾塑性破壊力学（EPFM）のパラメータで、船舶・橋梁・海洋構造物など靭性重視の設計で使われる指標だ。1963年にWells が提案した。
「亀裂先端が開く」というのをFEMで計算するんですか？ どうやって「亀裂先端」を特定するんでしょう？
FEMでは亀裂は「亀裂面上の2つの重なった節点（crack tip nodes）」として表現する。亀裂先端節点の対称面上での変位差が $\delta$（CTOD）になる。亀裂先端要素には特殊な四分の一点要素（Quarter-Point Element, QPE）を使って応力場の $r^{-1/2}$ 特異性を正確に再現するのが標準的だよ。
どんなケースでCTODを使い、どんなときにKIcを使うという判断基準はありますか？
塑性域の大きさで判断するのが原則だ。塑性域半径の推定式：
平面応力（薄板）では $r_p$ が試験片や亀裂長さに比べて無視できない大きさになる。目安として「亀裂長 $a$ に対して $r_p / a > 0.02$」になったらLEFMの精度が疑われ、CTODまたはJ積分ベースのEPFMを使う。実務では溶接構造の熱影響部（脆性転移温度付近）、極低温のLNGタンク、高エネルギー管路などがCTOD評価の典型対象だ。
CTODの定義と物理的意味
CTODの値が大きいほど「靭性が高い」ということですか？
そうだよ。CTODが大きいほど亀裂先端がより多く開いてから破壊が始まる——つまり吸収エネルギーが大きく靭性が高い。設計では材料試験から得られる臨界CTOD（$\delta_c$）が、FEM解析で得られる実際の構造物CTOD（$\delta_{appl}$）より大きければ「安全」と判定する：
ABS（米国規格協会）や英国規格 BS 7448 が主要な試験・評価規格だ。特に低温環境（極低温の LNG タンク、極地海洋構造）ではシャルピー試験だけでなくCTOD試験が必須になる。
$\delta_c$ は実験でどう測るんですか？ CTOD試験の試験片ってどんな形状ですか？
CTOD試験（BS 7448 / ASTM E1820）では三点曲げ試験片（SENB: Single Edge Notch Bend）またはコンパクト試験片（CT）を使う。試験片には疲労き裂を導入し、CODゲージ（クリップゲージ）で亀裂口の開口変位を計測する。き裂先端の $\delta_c$ はき裂口変位から幾何学的に換算する：
$V_p$ は塑性成分のき裂口変位、$r_p$ は塑性回転係数（通常0.4〜0.45）、$W$ は試験片幅、$a$ は亀裂長さ。試験温度を変えて $\delta_c$ の温度依存性（ダクタイル-脆性遷移）を評価するのが実務の標準的なアプローチだ。
J積分との関係
J積分とCTODって関係があるんですか？
弾性域では $K_I, J, \delta$（CTOD）は互いに換算できる：
$\sigma_{ys}$ は降伏応力、$m$ は拘束係数（平面応力 $m \approx 1$、平面ひずみ $m \approx 2$）、$E' = E$（平面応力）、$E' = E/(1-\nu^2)$（平面ひずみ）。塑性変形が大きい場合はJ積分の方が一般的だが、CTODは実験的に直接測定しやすい（亀裂開口変位計 COD Gauge を使用）という利点がある。
数値解法と実装
FEMによるCTOD計算
AbaqusやAnsysでCTODを計算するには何が必要ですか？
主要なポイントを挙げよう。
- 亀裂モデリング：亀裂面を seam（縫い目）として定義し、亀裂先端に四分の一点要素（QPE）を配置。Abaqusなら contour integral オプション、Ansysなら CINT コマンド。
- 亀裂先端メッシュ：亀裂先端から 0.01〜0.1mm の細かい扇状メッシュ（rosette mesh）。塑性域が大きい場合はそれ以上細かくする。
- 材料モデル：弾塑性解析が必要。等方硬化（isotropic）または移動硬化（kinematic）の降伏モデル。
- CTOD値の抽出：亀裂先端後方の特定距離（通常 $r = CTOD/2$ 相当位置）での亀裂面変位差を計測。または $J$ 値からCTODに換算。
Abaqusの contour integral って何を出力するんですか？ 複数の積分経路での $J$ 値が出るとのことですが。
Abaqusでは亀裂先端から外側へ向かって同心円状に複数の積分経路（Contour 1〜N）を自動設定し、各経路での $J$ 値・応力拡大係数 $K_I, K_{II}, K_{III}$ および $T$-応力を同時出力する。重要なのは「経路に依存しない（Path-independent）」ことの検証——弾性域では全経路で $J$ 値が一致するはずだが、塑性が大きいと経路依存が出る。この場合は複数経路の収束値または最外経路の値を採用するのが実務慣行だ。
四分の一点要素（QPE）と亀裂先端メッシュ
四分の一点要素（QPE）って何が特別なんですか？ 通常の要素と何が違うんですか？
通常の8節点四辺形要素や20節点六面体要素では節点が辺の中点に配置されるが、QPEでは亀裂先端に隣接する辺の中点節点を先端から1/4の位置（四分の一点）にずらす。これにより要素内の変位場が $r^{1/2}$（$r$ は亀裂先端からの距離）に比例するようになり、$\varepsilon \propto r^{-1/2}$、$\sigma \propto r^{-1/2}$ の特異性を正確に再現できる。通常要素では亀裂先端の応力が有限値に収束してしまうため、メッシュをいくら細かくしてもK値の精度が出ない。
QPEの代わりにXFEM（拡張有限要素法）を使う方法もありますか？
XFEMはメッシュに依存せず亀裂を表現できる方法で、亀裂伝播解析（伝播方向・速度の予測）に非常に有用だ。ただしCTODのような「亀裂先端の精密な変位計算」には、現状ではQPEの方が精度が高い傾向がある。XFEMの亀裂先端エンリッチメント関数が $r^{1/2}$ 特異性を近似的に表現するのに対し、QPEは要素形状で厳密に実現するためだ。実務では定常亀裂のCTOD評価にはQPE＋ rosette mesh、亀裂伝播の経路予測にはXFEM、と使い分けるのが現在の標準だ。
実践ガイド
FFS適合性評価フロー
海洋構造物の安全評価でCTODを使う場合のフローを教えてください。
適合性評価（Fitness for Service, FFS）のフローを示すよ。ECA（Engineering Critical Assessment）とも呼ばれる。
- 欠陥調査：超音波探傷（UT）や磁粉探傷（MT）で亀裂サイズ・形状を特定。
- 材料特性取得：使用温度でのCTOD試験（BS 7448またはASTM E1820）を実施。$\delta_c$ を取得。
- 応力解析：FEMで一次（荷重）+ 二次（熱・残留）応力を計算。
- 亀裂先端CTOD計算：FEMまたはBS 7910の簡易式でCTOD適用値 $\delta_{appl}$ を計算。
- 判定：$\delta_{appl} / \delta_c \le 1$（安全係数考慮後）ならば合格。
BS 7910の「簡易式」って具体的にどんな式ですか？ FEM解析なしで評価できるんですか？
BS 7910ではFAD（Failure Assessment Diagram）を用いた簡易評価が提供されている。FADは縦軸に $K_r = K_I / K_{mat}$（応力拡大係数比）、横軸に $L_r = \sigma / \sigma_Y$（荷重比）をとり、評価点がFAD曲線の内側に入れば安全と判定する。CTODベースなら $\delta_r = \delta_{appl} / \delta_c$ を横軸に置き換えた形式も使える。Level 1（保守的な簡易評価）からLevel 3（詳細FEM解析）まで段階的な評価レベルが定義されており、高リスク部位にはLevel 3のFEM解析を推奨する。
実践チェックリスト
- 材料特性：使用温度でのCTOD試験データ（δc）を取得、温度依存性を確認
- 亀裂モデリング：seam定義、QPE（四分の一点要素）配置を確認
- 亀裂先端メッシュ：rosette mesh、先端要素サイズ ≤0.1mm 以下
- 材料モデル：弾塑性（降伏曲線データ）を設定
- J積分の経路独立性確認：複数経路で値がほぼ一致（差 ≤5%）することを検証
- 残留応力の考慮：溶接部は残留応力分布をFEMまたは試験値から取得
- 安全係数の適用：規格（BS 7910 / API 579）指定の安全係数を乗じて判定
- 感度解析：亀裂サイズ・荷重・材料特性のバリエーションで評価の頑健性を確認
ソフトウェア比較
破壊力学解析に使える主要ツールを教えてください。
代表的なツールの比較を見てみよう。
| ツール | CTOD計算 | 亀裂モデリング | 特徴 |
|---|---|---|---|
| Abaqus/Standard | Contour Integral（J・SIF・CTOD） | Seam crack + QPE | J積分・SIF・CTODを多経路で自動計算、業界標準 |
| Ansys Mechanical | Fracture Mechanics（CINT） | SMART法（自動亀裂伝播） | 亀裂伝播（VCCT/CZM）も可能 |
| MSC Nastran（SOL 401） | CTOD（*CTOD） | CRACK要素 | 非線形破壊力学 |
| FRANC3D | SIF + J積分 | 自動亀裂挿入・伝播 | 亀裂伝播専用、疲労き裂成長予測に特化 |
| Code_Aster | K-J-CTOD | XFEM対応 | オープンソース、EDF開発、XFEM亀裂伝播対応 |
先端技術
CTODや弾塑性破壊力学の最新研究トレンドを教えてください。
いくつかの注目テーマがある。
- XFEM（拡張有限要素法）：メッシュに依存せずに亀裂を表現できるXFEMとCTODの組み合わせが普及。亀裂伝播経路をメッシュ再生成なしに追跡。
- 位相場破壊モデル（Phase-Field）：連続体として亀裂を表現。CTODに相当する情報を陰的に含む。複雑な亀裂分岐・合体に対応。
- 機械学習ベースの破壊靭性予測：材料組成と加工履歴からCTOD値をNNで予測する研究が急増。
- 水素脆化とCTOD：水素を含む環境（燃料電池スタック、水素タンク）での $\delta_c$ 低下予測にCTOD解析が活用。
水素脆化でCTODが低下するメカニズムはどういうものですか？
水素は亀裂先端の塑性域に拡散・濃縮し、転位の移動を妨げて（HELP機構：Hydrogen Enhanced Localized Plasticity）局所的な塑性変形を抑制する。その結果、CTODの値（すなわち亀裂先端が開くのに必要な変位量）が低下する——つまりより少ないエネルギーで亀裂が進展する「脆化」が起きる。高強度鋼（YS > 700MPa）ほど影響が大きく、水素タンク・パイプラインでの設計では水素環境での $\delta_c^H$（水素環境CTOD）を使った評価が必要になる。
リバティ船の脆性破壊とCTOD研究の始まり
1940年のタコマナローズ橋崩壊は空力振動の有名な事例だが、船舶や橋梁の突然の脆性破壊（特に低温下）も当時の工学者を悩ませていた。1943〜1946年に多数のリバティ船（溶接鋼製輸送船）が突然破断した事故の調査がCTOD研究の出発点のひとつ。「溶接熱影響部の脆化」という概念が確立し、CTODを用いた溶接継手の品質評価が現代の造船・橋梁規格の基礎になった。
トラブルシューティング
FEMでCTODを計算したら値がメッシュサイズに依存してしまいます。どうすればいいですか？
それはよくあるトラブルだよ。亀裂先端要素サイズの収束確認、QPE配置の確認、J積分の経路依存性確認など対策がある。
| 症状 | 主な原因 | 対策 |
|---|---|---|
| J値がメッシュサイズに依存 | QPEなし、または通常要素のサイズが粗すぎ | 亀裂先端にQPE配置、rosette meshサイズを1/2〜1/4に縮小してメッシュ収束確認 |
| J値が経路依存（弾性域でも） | 塑性変形の大きさ、または解析誤差 | 弾性域で経路独立性を確認してから塑性解析へ、最外経路の値を採用 |
| CTODが過大評価される | 亀裂先端変位の計測位置が不適切 | J値からCTODへの換算を使う（δ = J/(mσ_ys)） |
| 非線形解析が収束しない | 亀裂先端の塑性域で変位ジャンプ | 亀裂先端メッシュをさらに細化、荷重増分を小さくする |
| 3D解析でのK3（モードIII）の影響 | 実際の溶接部など混合モード亀裂 | T-応力・混合モードSIFを出力してMTS（最大接線応力）基準で評価方向を決定 |
メッシュ依存性、亀裂先端特異性の再現失敗、J積分の経路依存性など詳細解説