Hashin破壊基準 — トラブルシューティングガイド
より充実した内容を hashin-criterion.html でご覧いただけます。
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%程度）を設計値に過剰に適用していると保守的すぎる予測になる。ドライ条件の実験値を基準に、使用環境に応じた低減係数を別途適用するのが正しいアプローチだ。
トラブル解決の考え方
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——Hashin破壊基準の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告