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設計で最も精度の高い基準として評価された。
各項の物理的意味
- 慣性項（質量項）：$\rho \ddot{u}$、つまり「質量×加速度」。急ブレーキで体が前に投げ出された経験はありませんか？ あの「持っていかれる感じ」がまさに慣性力です。重い物体ほど動き出しにくく、動き出したら止まりにくい。地震で建物が揺れるのも、地面が急に動いたのに建物の質量が「置いていかれる」から。静解析ではこの項をゼロにしますが、それは「ゆっくり力をかけるから加速度は無視できる」という仮定です。衝撃荷重や振動問題では絶対に省略できません。
- 剛性項（弾性復元力）：$Ku$ や $\nabla \cdot \sigma$。ばねを引っ張ると「戻ろうとする力」を感じますよね？ あれがフックの法則 $F=kx$ であり、剛性項の本質です。では質問——鉄の棒とゴム紐、同じ力で引っ張るとどちらが伸びるでしょうか？ 当然ゴムです。この「伸びにくさ」がヤング率 $E$ であり、剛性を決めます。よくある勘違い：「剛性が高い＝強い」ではありません。剛性は「変形しにくさ」、強度は「壊れにくさ」で、別の概念です。
- 外力項（荷重項）：体積力 $f_b$（重力など）と表面力 $f_s$（圧力、接触力など）。こう考えてみてください——橋の上のトラックの重さは「中身全体にかかる力」（体積力）、タイヤが路面を押す力は「表面だけにかかる力」（表面力）。風圧、水圧、ボルトの締付力…すべて外力です。ここでありがちな失敗：荷重の方向を間違える。「引張」のつもりが「圧縮」になっていた——笑い話に聞こえますが、3D空間で座標系が回転していると実際に起こります。
- 減衰項：レイリー減衰 $C\dot{u} = (\alpha M + \beta K)\dot{u}$。ギターの弦を弾いてみてください。音は鳴り続けますか？ いいえ、徐々に小さくなりますよね。振動エネルギーが空気抵抗や弦の内部摩擦で熱に変わるからです。車のショックアブソーバーも同じ原理——わざと振動エネルギーを吸収して乗り心地を良くしています。もし減衰がゼロだったら？ 建物は地震の後いつまでも揺れ続けることになります。実際にはそうならないので、適切な減衰の設定が重要です。
仮定条件と適用限界
次元解析と単位系
| 変数 | SI単位 | 注意点・換算メモ |
|---|---|---|
| 変位 $u$ | m（メートル） | mm入力時は荷重・弾性率もMPa/N系に統一すること |
| 応力 $\sigma$ | Pa（パスカル）= N/m² | MPa = 10⁶ Pa。降伏応力との比較時に単位系の不一致に注意 |
| 歪み $\varepsilon$ | 無次元（m/m） | 工学歪みと対数歪みの区別に注意（大変形時） |
| 弾性率 $E$ | Pa | 鋼: 約210 GPa、アルミ: 約70 GPa。温度依存性に注意 |
| 密度 $\rho$ | kg/m³ | mm系ではtonne/mm³（= 10⁻⁹ tonne/mm³ for 鋼） |
| 力 $F$ | N（ニュートン） | mm系ではN、m系ではNで統一 |
数値解法と実装
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が代表的な入力値だ。
線形要素（1次要素）
節点間を線形補間。計算コストは低いが、応力の精度が低い。せん断ロッキングに注意（低減積分やB-bar法で緩和）。
2次要素（中間節点付き）
曲線的な変形を表現可能。応力精度が大幅に向上するが、自由度は約2〜3倍に増加。推奨：応力評価が重要な場合。
完全積分 vs 低減積分
完全積分：過剰拘束（ロッキング）のリスク。低減積分：アワーグラスモード（零エネルギーモード）のリスク。適材適所で選択。
アダプティブメッシュ
誤差指標（ZZ推定量等）に基づく自動細分化。応力集中部の精度を効率的に向上。h法（要素分割）とp法（次数増加）がある。
ニュートン・ラフソン法
非線形解析の標準的手法。接線剛性マトリクスを毎反復更新。収束半径内で2次収束するが、計算コストが高い。
修正ニュートン・ラフソン法
接線剛性マトリクスを初期値または数反復毎に更新。各反復のコストは低いが、収束速度は線形的。
収束判定基準
力の残差ノルム: $||R|| / ||F_{ext}|| < \epsilon$（一般に $\epsilon = 10^{-3}$〜$10^{-6}$）。変位増分ノルム: $||\Delta u|| / ||u|| < \epsilon$。エネルギーノルム: $\Delta u \cdot R < \epsilon$
荷重増分法
全荷重を一度に負荷せず、小刻みに増加させる。弧長法（Riks法）は荷重-変位関係の極値点を越えて追跡可能。
直接法 vs 反復法のたとえ
直接法は「連立方程式を筆算で正確に解く」方法——確実だが大規模問題では時間がかかりすぎる。反復法は「当て推量を繰り返して正解に近づく」方法——最初は大雑把な答えだが、反復するたびに精度が上がる。辞書で言葉を探すとき、最初のページから順番に探す（直接法）より、見当をつけて開き、前後に調整する（反復法）方が効率的なのと同じ原理。
メッシュの次数と精度の関係
1次要素は「定規で曲線を近似する」——直線の折れ線で表現するため精度に限界がある。2次要素は「フレキシブルカーブ」——曲線的な変化を表現でき、同じメッシュ密度でも格段に精度が向上する。ただし、1要素あたりの計算コストは増えるため、トータルのコスト対効果で判断する。
実践ガイド
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以下を全荷重ケースで保証している。
解析フローのたとえ
解析の流れは、実は料理とそっくりです。まず材料を買い出し（CADモデルの準備）、下ごしらえをして（メッシュ生成）、火にかけて（ソルバー実行）、最後に盛り付ける（後処理で可視化）。ここで大事な問いかけ——料理で一番失敗しやすい工程はどこでしょう？ 実は「下ごしらえ」なんです。メッシュの品質が悪いと、どんなに優秀なソルバーを使っても結果はめちゃくちゃになります。
初心者が陥りやすい落とし穴
あなたはメッシュ収束性を確認していますか？ 「計算が回った＝結果が正しい」と思っていませんか？ これ、実はCAE初心者が最も陥りやすい罠です。ソルバーは与えられたメッシュで「それなりの答え」を必ず返します。でもメッシュが粗すぎれば、その答えは現実から大きくずれている。最低3段階のメッシュ密度で結果が安定することを確認する——これを怠ると「コンピュータが出した答えだから正しいはず」という危険な思い込みに陥ります。
境界条件の考え方
境界条件の設定は、試験の「問題文を書く」のと同じです。問題文が間違っていたら？ どんなに正確に計算しても答えは間違いますよね。「この面は本当に完全固定なのか」「この荷重は本当に一様分布なのか」——現実の拘束条件を正しくモデル化することが、実は解析全体で最も重要なステップだったりします。
ソフトウェア比較
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%精度で予測した。
選定で最も重要な3つの問い
- 「何を解くか」：Hashin破壊基準に必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
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損傷解析でよくあるトラブルを教えてください。
プログレッシブ損傷解析は非線形解析の中でも最も収束が難しい分野だ。
収束困難
損傷が始まると解析が収束しなくなります。
損傷による急激な剛性低下で平衡経路が不安定になる。対策：
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の要素が必要。計算コストが大きくなるが、精度のために避けられない。
破壊エネルギーの不足
破壊エネルギーが正しいかどうかどう確認しますか？
メッシュを変えると損傷パターンが変わります。
破壊エネルギーの正則化が正しく機能しているか確認。特に：
0.5〜2 mmの要素って、かなり細かいメッシュですね。
プログレッシブ損傷解析はメッシュ要件が厳しい。OHC試験のシミュレーションで穴の周囲に0.5 mmの要素が必要。計算コストが大きくなるが、精度のために避けられない。
破壊エネルギーが正しいかどうかどう確認しますか？
まとめ
Hashin損傷のトラブル対処、整理します。
プログレッシブ損傷解析は「FEM上級者の領域」ですね。
その通り。材料特性の理解、非線形解析のスキル、実験データとの整合…全てが揃って初めてプログレッシブ損傷解析が意味を持つ。
Hashin基準が実験破損より早期に予測する問題
Hashin基準が実験のFirst Ply Failure（FPF）より20〜30%低い荷重で破損を予測する場合、マトリクス引張強度YTまたは面内せん断強度S12の値が低く設定されていることが多い。特に湿熱吸収後の強度低下（YTで最大30%、S12で20%程度）を設計値に過剰に適用していると保守的すぎる予測になる。ドライ条件の実験値を基準に、使用環境に応じた低減係数を別途適用するのが正しいアプローチだ。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——Hashin破壊基準の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告