圧壊フォーム材モデル
圧壊フォーム材の理論基礎
フォーム材の力学
先生、フォーム材（発泡体）の力学は金属とどう違いますか？
フォーム材（EPS、PUフォーム、金属フォーム）は圧縮で大きく体積が減少する。金属は非圧縮塑性（$\Delta V = 0$）だが、フォームは圧縮性塑性（$\Delta V \neq 0$）。セル壁の座屈と圧壊でエネルギーを吸収。
圧縮応力-ひずみ曲線
フォームの典型的な圧縮曲線：
1. 弾性域 — セル壁の弾性変形
2. プラトー域 — セル壁の座屈/圧壊。応力がほぼ一定で大変形
3. 緻密化域 — セル壁が密着。急激に応力上昇
FEMでのモデル化
まとめ
フォームモデル誕生の動機
Deshpande & Fleck（ケンブリッジ大）が2000年に発表した「Isotropic constitutive models for metallic foams」は、アルミニウム発泡金属（Alporas、Cymat）の衝撃吸収を定量化するために開発された。それ以前のvon Misesモデルではフォームの等方圧縮による大変形を全く表現できなかったため、体積塑性ひずみを許容する新たな降伏面が必要だった。この論文はAbaqus「Crushable Foam」モデルの理論的基盤となっている。
圧壊フォーム材の数値計算手法
フォームのFEM設定
```
*MAT_LOW_DENSITY_FOAM
$ 圧縮応力-ひずみテーブル
*DEFINE_CURVE
0., 0.
0.1, 0.5
0.5, 0.6
0.8, 5.0
```
Abaqus: *CRUSHABLE FOAM + テーブル入力。
まとめ
応力比パラメータk₀の実測法
Crushable Foamの降伏面形状を決めるパラメータk₀（初期静水圧降伏応力比）は、一軸圧縮試験と静水圧圧縮試験の2種類の実験から決定する。Alporas（住友電工製アルミフォーム、密度0.25 g/cm³）の場合、一軸圧縮降伏応力≈1.6 MPa、静水圧降伏≈1.9 MPaより k₀≈1.19が得られることがDeshpane自身の実験で報告されている。k₀が不明な場合は最初にk₀=1.1〜1.3を試すのが実務上の慣例である。
圧壊フォーム材の実務適用
フォームの実務
パッケージング（電子機器の緩衝材）、自動車のバンパーフォーム、ヘルメットのEPSライナー。
実務チェックリスト
自動車シートクッション衝撃吸収
ポリプロピレン（PP）発泡フォーム（密度30〜60 kg/m³）はヘッドレスト・膝パッドなどの乗員保護部品に使われる。Toyota・Honda・Volkswagenは2010年代からFMVSS201U（乗用車内装衝撃規格）準拠のためにAbaqusのCrushable Foamモデルによる仮想試験を活用しており、樹脂発泡材の圧縮-密実化曲線（応力〜100%ひずみ）を入力するだけで衝撃加速度-時刻歴を±10%以内に再現できるとメーカー各社が報告している。
圧壊フォーム材のソフトウェア比較
フォームのツール
LS-DYNAとAbaqusの発泡材実装の違い
LS-DYNAの代表的な発泡材モデルはMAT_LOW_DENSITY_FOAM（#57）で、ひずみ速度依存の積載-除荷曲線をそのまま入力できる準物理モデルである。一方AbaqusのCrushable Foamは降伏面の数学的定式化に基づくため、実験曲線からk₀やα等のパラメータを逆解析で同定する必要がある。落下衝撃や爆発など高速現象ではLS-DYNA MAT57、クリープ変形や長期疲労ではAbaqus CDP系が向いているとの実務的な使い分けが業界でコンセンサスを得ている。
圧壊フォーム材の先端研究
フォームの先端
異方性フォームへの拡張
Deshpande-Fleckモデルは等方性を前提とするが、押し出し成形フォームや繊維複合フォームは厚み方向と面方向で降伏応力が2〜5倍異なる場合がある。2006年にChenとLuは「Orthotropic Crushable Foam」を提案し、3方向独立の降伏応力を定義する形式を導入した。この定式化はLS-DYNA MAT_ORTHOTROPIC_CRUSHABLE_FOAM（Material #142）に実装されており、ハニカムパネルのOOB（面外圧縮）解析で用いられている。
圧壊フォーム材のトラブル対応
フォームのトラブル
密実化後の急峻硬化と収束問題
フォームは圧縮ひずみが70〜80%を超えると密実化（densification）が始まり、応力が指数的に増加する。この急峻な硬化は増分解析の収束を妨げることが多い。対策として、密実化領域をテーブル入力の最大ひずみ点でカットオフし、それ以上を線形外挿に切り替える方法が使われる。AbaqusのCrushable Foam入力では応力-ひずみテーブルの最終勾配が自動で外挿されるため、テーブルを密実化開始点で適切に打ち切ることが解析安定の要となる。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——圧壊フォーム材モデルの問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告