Hashin破壊基準
Hashin破壊基準の理論基礎
Hashin基準とは
先生、Hashin基準はTsai-Wuとどう違うんですか？
最大の違いは破壊モードを区別することだ。Tsai-Wuは「破壊した/しない」の二値判定だが、Hashin基準（1980）は4つの独立した破壊モードを識別する。
4つの破壊モード
Hashin基準の4モード：
1. 繊維引張破壊（Fiber Tension）
2. 繊維圧縮破壊（Fiber Compression）
3. マトリクス引張破壊（Matrix Tension）
4. マトリクス圧縮破壊（Matrix Compression）
繊維とマトリクスの破壊を区別し、さらに引張と圧縮でも分ける。4つの式で4モードを判定するんですね。
そう。これがHashin基準の最大の利点だ。破壊モードがわかれば、どの強度を上げればいいかが明確になる。繊維破断なら繊維量を増やす、マトリクスクラックならマトリクスの樹脂を変える…と対策が具体的になる。
プログレッシブ損傷への対応
Hashin基準はプログレッシブ損傷解析に使えるんですか？
使える。各モードの破壊指標が1を超えたとき、そのモードに対応する剛性を低減する：
- 繊維破壊 → $E_1$ を低減
- マトリクスクラック → $E_2, G_{12}$ を低減
- 両方の組み合わせ → 全剛性を低減
AbaqusのBuilt-in Hashin Damageは、破壊判定＋剛性低減＋エネルギー散逸を一体で実装している。複合材のプログレッシブ損傷解析の事実上の標準だ。
Tsai-Wuにはこの機能がないんですね。
Tsai-Wuは1つの式で全モードを混ぜているため、「どの剛性を低減すべきか」が決まらない。モード分離ができるHashin基準がプログレッシブ損傷の基盤になっている。
Hashin基準の限界
Hashin基準にも限界がある：
| 限界 | 説明 |
|---|---|
| 繊維圧縮モードが単純すぎる | 実際はキンクバンド（繊維の局所座屈）が支配 |
| マトリクス圧縮で破壊面角度を考慮しない | Puck基準やLaRC基準のほうが正確 |
| 層間剥離を扱えない | 別途CZM（Cohesive Zone Model）が必要 |
Hashinでもまだ改善の余地があるんですね。
Hashin基準は「モード分離の先駆者」だが、各モードの物理描写はまだ粗い。LaRC05基準（NASAランリー研究所、2005年）はHashinの限界を改善した最新基準だ。
まとめ
Hashin基準を整理します。
要点：
- 4つの破壊モードを区別 — 繊維引張/圧縮、マトリクス引張/圧縮
- プログレッシブ損傷に対応 — モード別の剛性低減が可能
- Abaqusに標準実装 — Built-in Hashin Damage
- Tsai-Wuより物理的に正確 — ただし計算が複雑
- 限界もある — 繊維圧縮のキンクバンド、マトリクス圧縮の破壊面角度
Tsai-Wuがスクリーニング、Hashinが詳細評価、という使い分けが明確になりました。
その通り。設計の初期段階でTsai-Wuで概算し、詳細設計でHashin（またはPuck/LaRC）で精密評価する。段階的なアプローチが実務的だ。
HashinがCFRP破損の2軸相互作用を発見
Hashin破損基準は1973〜1980年にZvi Hashin（テクニオン大学）が開発した。従来のMaximum Stress基準では多軸応力場での複合材破損が不正確だったため、繊維破断・マトリクス破断・繊維圧縮・マトリクス圧縮の4モードをそれぞれ独立した相互作用式で定義した。その後WWFE（World Wide Failure Exercise、2002〜2004年）でCFRP設計で最も精度の高い基準として評価された。
Hashin破壊基準の数値計算手法
AbaqusのHashin損傷モデル
AbaqusのHashin損傷モデルの設定を教えてください。
AbaqusのBuilt-in Hashin Damageは2段階で設定する。
Step 1: 損傷開始基準（Damage Initiation）
```
*DAMAGE INITIATION, CRITERION=HASHIN
X_t, X_c, Y_t, Y_c, S_L, S_T, alpha
```
7つのパラメータ: 繊維引張/圧縮強度、マトリクス引張/圧縮強度、縦/横せん断強度、$\alpha$（せん断応力の寄与係数）。
Step 2: 損傷進展則（Damage Evolution）
```
*DAMAGE EVOLUTION, TYPE=ENERGY
G_ft, G_fc, G_mt, G_mc
```
4つの破壊エネルギー: 繊維引張/圧縮、マトリクス引張/圧縮の破壊エネルギー（$G_c$）。
破壊エネルギーが必要なんですね。強度値だけでは足りない？
ソルバー別のHashin実装
| 機能 | Abaqus | Nastran | Ansys |
|---|---|---|---|
| Hashin破壊判定 | ○（標準） | △（USDFLD） | ○（ACP Post） |
| プログレッシブ損傷 | ○（DAMAGE EVOLUTION） | △（SOL 400 + ユーザー） | △（ACP + APDL） |
| 破壊エネルギー法 | ○ | — | 限定的 |
| 要素削除 | ○（STATUS） | ○（PARAM,ERODEL） | ○（EKILL） |
Abaqusが圧倒的に充実していますね。
Hashinベースのプログレッシブ損傷解析はAbaqusが事実上の標準。論文で「Hashin damage」と言えばほぼAbaqusの実装を指す。NastranやAnsysでは同等のことがユーザーサブルーチンやスクリプトで可能だが、手間がかかる。
メッシュ依存性と正則化
破壊エネルギーによる正則化について教えてください。
損傷が局所化（1要素に集中）すると、エネルギー散逸が要素サイズに依存する。特性長さ（characteristic length）を用いて破壊エネルギーを要素サイズで正規化する：
$L_c$ は要素の特性長さ（概ね要素サイズ）。これにより、メッシュを変えても散逸エネルギーが一定になる。
要素サイズが小さいと破壊ひずみが大きくなり、大きいと小さくなる。エネルギーが保存される。
完璧な理解だ。ただし要素が極端に大きい（特性長さ > $2G_c/\sigma_c^2 \cdot E$）と、スナップバックが起きて数値的に不安定になる。メッシュサイズの上限に注意。
まとめ
Hashin基準の数値手法、整理します。
要点：
- AbaqusのDAMAGE INITIATION + DAMAGE EVOLUTIONが標準 — Hashin損傷の事実上の業界標準
- 7つの材料強度＋4つの破壊エネルギー — 合計11パラメータ
- 破壊エネルギーによる正則化 — メッシュ依存性を排除
- 要素サイズの上限に注意 — 大きすぎるとスナップバック
- NastranやAnsysではユーザーサブルーチンが必要 — Abaqusが最も手軽
Hashin基準の4モード計算手順
Hashin基準では①繊維引張(σ11/XT)²+(τ12/Xt)²≥1で繊維破断②繊維圧縮(σ11/XC)²≥1でkinking③マトリクス引張(σ22/YT)²+(τ12/S12)²≥1で横クラック④マトリクス圧縮計算の4モードを全て評価し、最初に基準値1を超えたモードが初期破損モードとなる。CFRPのXT=2000MPa・XC=1500MPa・YT=50MPa・S12=80MPaが代表的な入力値だ。
Hashin破壊基準の実務適用
Hashin基準の実務適用
Hashin基準は実務でどう使われていますか？
主にプログレッシブ損傷解析（PDA: Progressive Damage Analysis）の基盤として使われる。
PDAの適用場面
| 場面 | 目的 |
|---|---|
| 航空機のCAI（Compression After Impact） | 衝撃後の圧縮強度を予測 |
| OHC/OHT（Open Hole Compression/Tension） | 穴あき板の引張/圧縮強度 |
| ノッチ効果 | 切り欠きでの破壊進展 |
| 複合材ボルト接合 | ベアリング破壊のシミュレーション |
CAI（衝撃後圧縮）は航空機認証で最も重要な試験の一つですよね。
その通り。複合材構造の設計許容値はCAI強度で決まることが多い。Hashin＋CZM（層間剥離）のPDAでCAI強度を予測し、試験数を削減する「ビルディングブロックアプローチ」が航空宇宙の標準ワークフローだ。
解析のワークフロー
2. 損傷状態の転写 — 衝撃後の損傷変数をAbaqus/Standardに転写
3. 圧縮解析（陰解法） — 損傷を含んだ状態で圧縮荷重を与え、残留圧縮強度を評価
陽解法と陰解法を切り替えるんですね。
衝撃は短時間の動的現象だから陽解法、圧縮は準静的だから陰解法。Abaqusの*IMPORTで陽解法→陰解法の状態転写が可能。
材料パラメータの取得
破壊エネルギー $G_c$ はどうやって測定しますか？
| パラメータ | 試験方法 | 規格 |
|---|---|---|
| $G_{Ic}$（モードI層間） | DCB（Double Cantilever Beam） | ASTM D5528 |
| $G_{IIc}$（モードII層間） | ENF（End-Notched Flexure） | ASTM D7905 |
| $G_{ft}$（繊維引張） | CT試験や推定 | 直接測定困難 |
| $G_{mc}$（マトリクス圧縮） | Compact Compression | 研究段階 |
繊維の破壊エネルギーは直接測定が困難？
繊維破壊のエネルギーは非常に大きく（$G_{ft} \approx 50 \sim 100$ kJ/m²）、標準試験法がない。文献値や逆解析（FEMの結果と試験結果を合わせる）で推定することが多い。
実務チェックリスト
Hashin基準のチェックリストをお願いします。
破壊エネルギーの設定が最も重要で、最も難しいパラメータですね。
強度値は標準試験で比較的容易に取得できるが、破壊エネルギーは試験が困難または高コスト。パラメータの不確かさに対する感度分析が不可欠だ。
CFRP翼パネルの破損包絡線設計
航空機CFRP翼パネルの多軸荷重下での安全マージンを「破損包絡線（Failure Envelope）」で管理する。FEM + Hashin基準でσ11-σ22の荷重空間上に安全域を描き、飛行荷重の全ケースが包絡線内に収まることを確認する。Bombardier C-Series（現A220）のCFRP翼スキンはこの設計で軽量化達成し、ALCF（許容設計荷重比率）0.85以下を全荷重ケースで保証している。
Hashin破壊基準のソフトウェア比較
Hashin基準のツール
Hashin基準を使えるツールを比較してください。
AbaqusとLS-DYNAが強いですね。
LS-DYNAのMAT58は Chang-Chang 基準（Hashinの変形版）をベースにしたプログレッシブ損傷モデル。衝突時の複合材破壊に広く使われている。
選定ガイド
Hashin基準はAbaqusで最も完成度が高いということですね。
AbaqusのHashin Damage実装は、論文で最も引用される複合材損傷モデルだ。同じ解析を他のソルバーで行うにはかなりの開発努力が必要。
Abaqus HashinDamage材料モデル
Abaqusの*DAMAGE INITIATION, CRITERION=HASHINはHashin基準を自動的に4モードで評価し、初期破損後に材料劣化（剛性低減）を進行させるプログレッシブ解析を一括実行できる。DAMAGE EVOLUTIONカードでエネルギーベースの剛性低減率を設定する。GKN Aerospace社はCFRP機体部品の認証解析にAbaqus + Hashinを使用し、Boeing 787の翼形成型パネルの静的破壊試験を±8%精度で予測した。
Hashin破壊基準の先端研究
Hashin基準を超える破壊理論
Hashin基準の限界を改善する研究はありますか？
Hashin基準の各モードをより物理的に正確にした次世代基準が開発されている。
LaRC05基準（NASAランリー）
NASAランリー研究所が開発したLaRC05（Langley Research Center Criteria, 2005）は、Hashinの4モードを物理ベースで改善：
- 繊維圧縮 — キンクバンド理論（繊維の微小座屈）。Hashinの単純な圧縮基準より正確
- マトリクス圧縮 — Mohr-Coulombベースの破壊面理論。Puck基準と同等の精度
- 繊維引張 — Hashinと同様
- マトリクス引張 — Hashinと同様だが横せん断の扱いを改善
キンクバンドって何ですか？
繊維が圧縮を受けると、繊維の初期不整（微小なミスアライメント）が増幅されて局所的な座屈が起きる。これがキンクバンドだ。圧縮強度は繊維自体の強度ではなく、繊維の安定性（初期不整と樹脂のせん断剛性）で決まる。
Hashinの $\sigma_1 / X_c$ では表現できない物理なんですね。
キンクバンドの圧縮強度は：$X_c \approx G_{12} / (1 + \phi_0 / \gamma_c)$。樹脂のせん断剛性 $G_{12}$ と初期繊維ミスアライメント $\phi_0$ で決まる。Hashin基準では見えない物理だ。
Puck基準
Puck基準（ドイツ、2002年）はマトリクス破壊の破壊面角度を考慮する。マトリクス圧縮では破壊が板厚方向の斜め面に沿って起きるため、破壊面の角度を変数として最も危険な面を探索する。
WWFE（World-Wide Failure Exercise）での比較では、Puck基準が最も実験結果に近い予測を示した。ただし計算コストが高い（破壊面角度の最適化が必要）。
Phase-Field法による複合材破壊
Phase-Field法を複合材に適用し、繊維破断とマトリクスクラックの進展を連続場で追跡する研究が急増している。Hashinのようなモード分離を事前に仮定せず、破壊経路が自然に決定される。
Phase-Fieldなら破壊モードの仮定が不要？
理想的にはそう。ただし複合材の異方性をPhase-Field法に組み込む定式化がまだ研究段階だ。
まとめ
Hashin基準を超える破壊理論、まとめます。
Hashin基準は複合材破壊の「出発点」であり、より精密な予測にはLaRC05/Puck/Phase-Fieldに進む。
Hashin基準の3D拡張：三次元複合材
Hashin基準の2D版は面外応力成分σ33・τ23・τ13を扱えないため、厚肉積層板や3D繊維強化材には不足だ。Puck（2001年）やCuntze（2004年）が三次元破損基準を提案しており、特に斜め衝突荷重を受ける3次元CFRP（炭素繊維織物・3D編物）の破損予測精度を改善した。宇宙機の分厚いCFRP結合部設計でこれらの3D拡張基準が2010年代から採用されている。
Hashin破壊基準のトラブル対応
Hashin損傷解析のトラブル
Hashin損傷解析でよくあるトラブルを教えてください。
プログレッシブ損傷解析は非線形解析の中でも最も収束が難しい分野だ。
収束困難
損傷が始まると解析が収束しなくなります。
損傷による急激な剛性低下で平衡経路が不安定になる。対策：
1. 粘性正則化 — *DAMAGE STABILIZATION で微小な粘性を追加。Abaqusの推奨値は $10^{-5}$ 程度
2. 増分を小さくする — 最小増分を $10^{-10}$ 程度に
3. 安定化法 — *STATIC, STABILIZE
4. 陽解法に切り替え — 収束問題を回避
粘性正則化の値はどう決めますか？
小さいほど物理的に正確だが収束しにくい。大きいと収束するが結果が不正確。$10^{-5}$ から始めて、結果が粘性値に依存しないことを確認する（$10^{-4}$ と $10^{-6}$ で同じ結果か）。
メッシュ依存性
メッシュを変えると損傷パターンが変わります。
破壊エネルギーの正則化が正しく機能しているか確認。特に：
- 特性長さが正しいか — Abaqusは自動計算するが、歪んだ要素では不正確
- 要素サイズが破壊プロセスゾーンに対して適切か — マトリクスクラックのプロセスゾーンは0.5〜2 mm程度。要素サイズはこのオーダーにすべき
0.5〜2 mmの要素って、かなり細かいメッシュですね。
プログレッシブ損傷解析はメッシュ要件が厳しい。OHC試験のシミュレーションで穴の周囲に0.5 mmの要素が必要。計算コストが大きくなるが、精度のために避けられない。
破壊エネルギーの不足
破壊エネルギーが正しいかどうかどう確認しますか？
まとめ
Hashin損傷のトラブル対処、整理します。
プログレッシブ損傷解析は「FEM上級者の領域」ですね。
その通り。材料特性の理解、非線形解析のスキル、実験データとの整合…全てが揃って初めてプログレッシブ損傷解析が意味を持つ。
Hashin基準が実験破損より早期に予測する問題
Hashin基準が実験のFirst Ply Failure（FPF）より20〜30%低い荷重で破損を予測する場合、マトリクス引張強度YTまたは面内せん断強度S12の値が低く設定されていることが多い。特に湿熱吸収後の強度低下（YTで最大30%、S12で20%程度）を設計値に過剰に適用していると保守的すぎる予測になる。ドライ条件の実験値を基準に、使用環境に応じた低減係数を別途適用するのが正しいアプローチだ。
関連トピック
なった
詳しく
報告