全体マップ

本サイトは3部構成です。基礎編で一般の構造解析と共通する土台を固め、衝突コア編（本サイトの主役）で衝突解析ならではの理論に踏み込み、実務・検証編で「正しく動かし、正しく疑う」方法を学びます。

読み進めルート

このサイトの読み方

• KEY EQUATION ボックス — その章で最も重要な式。式の「読み方」を言葉で添えています。
• インタラクティブ図解 — スライダーやボタンで自分の手を動かして数式の意味を確かめられます。
• ミニ演習 — 各章末の理解確認。クイズは選択肢を押すと解説が開き、計算問題は「解答を見る」で答え合わせできます。
• 用語に迷ったら付録の用語集・記号一覧へ。

学習目標

衝突解析が「何を予測する技術」かを知り、実車試験との関係を理解する。

前提知識

特になし。高校物理の力学のイメージがあれば十分。

読了目安

約15分

1.1 衝突現象の特徴 —— 短時間・大変形・非線形・接触

自動車の前面衝突は、車体が壁に触れてから乗員の運動が落ち着くまで、わずか100〜150ミリ秒で終わります。まばたき（約300ms）よりも短いこの時間の中で、次のことが同時に起こります。

• 短時間 —— 現象全体がミリ秒（1/1000秒）のスケール。荷重は静かに釣り合うのではなく、衝撃として構造を駆け抜けます。
• 大変形 —— フロント部材は元の長さの半分以下まで潰れます。「変形は小さい」という材料力学の仮定は最初の1msで破綻します。
• 非線形 —— 材料は降伏して塑性変形し（材料非線形）、形状が大きく変わって剛性そのものが変化し（幾何非線形）、部品同士が触れたり離れたりします（境界非線形）。3種類の非線形性が全部同時に現れます。
• 接触 —— バンパーと壁、折り重なるパネル同士、乗員とエアバッグ。「どこが何に触れるか」が事前に分からないのが衝突の本質です。

この4つの特徴こそが、後の章で学ぶ陽解法（第3・4章）・接触アルゴリズム（第5章）・ひずみ速度依存の材料モデル（第6章）という「衝突解析ならではの技術」を要求します。本章ではまず現象の時間感覚を掴みましょう。

1.2 応用分野 —— 衝突解析はどこで使われているか

共通するのは「短時間に大きく変形し、接触を伴う」こと。つまり本サイトで学ぶ理論は、自動車に限らず「モノがぶつかる現象」一般に使えます。

1.3 実車試験 vs 仮想試験

「実物をぶつけて測ればいいのでは？」——もちろん最終確認は実車試験で行います。しかし開発の途中では、シミュレーション（仮想試験）に大きな利点があります。

1.4 主要ソフトの位置づけ

衝突解析の商用ソルバーとして広く使われているのは次の4つです。名前だけ覚えておけば、文献や求人票が読めるようになります。

• LS-DYNA —— 衝突解析のデファクトスタンダード。自動車業界で最も広く使用。
• PAM-CRASH —— 欧州系メーカーで実績の長い衝突専用ソルバー。
• Altair Radioss —— 同じく衝突・衝撃に強い陽解法ソルバー。ソースコードが OpenRadioss としてオープンソース公開されており、無料で本物の衝突ソルバーを動かせます。本サイトでは各章に「OpenRadiossで確かめる」ボックスを置き、第11章で実行から結果分析まで実践します。
• Abaqus/Explicit —— 汎用FEMソフトAbaqusの陽解法モジュール。

いずれも中身は本サイトで学ぶ陽解法・接触・材料モデルの組み合わせです。理論を一度理解すれば、どのソフトに触れても「設定項目が何を意味するか」が読み解けるようになります。

学習目標

衝突で効いてくる「大変形」と「保存則」を再確認し、運動方程式の各項の意味を言えるようになる。

前提知識

材料力学の初歩（応力・ひずみ・ヤング率）、微分の基本。

読了目安

約25分

2.1 応力とひずみ —— 公称値と真値の違い

材料力学で習う応力・ひずみは、多くの場合公称応力・公称ひずみ（工学応力・工学ひずみ）です。長さ L₀、断面積 A₀ の棒を力 F で引っ張って長さ L になったとき：

これは「変形前の寸法で割った」量です。変形が数％以下なら問題ありませんが、衝突では部材が50％も潰れ、断面積も大きく変わります。そこで「いまこの瞬間の寸法」を基準にした真応力・真ひずみを使います。

数値を入れると差が実感できます。公称ひずみ 5% なら真ひずみは ln(1.05) ≈ 4.9% でほぼ同じ。しかし公称ひずみ 100%（長さ2倍）では真ひずみは ln 2 ≈ 69% と大きくずれます。衝突解析の材料データ（第6章）は基本的に真応力−真ひずみで入力するので、この換算は実務でも毎回登場します。

2.2 大変形・大回転 —— 微小変形の仮定はどこで破れるか

微小変形理論は、次の2つを暗黙に仮定しています。

1. 釣り合い式を変形前の形状で立ててよい（形が変わっても力の向きや腕の長さは変わらない）。
2. ひずみは変位の1次の項だけで表せる（回転しても伸びていなければひずみゼロ、が近似的に成り立つ）。

衝突ではどちらも成り立ちません。潰れたフロントサイドメンバーは折れ曲がり（大回転）、力の伝達経路そのものが刻々と変わります。剛体回転しただけの要素に「見かけのひずみ」が生じないように、変形を追いかけながら座標系を更新する定式化（有限変形理論）が必要になります。

2.3 運動方程式 —— すべての章がこの式に帰着する

FEMで離散化した構造の運動は、次の式で表されます。本サイトで最も重要な式です。

この式の見方で、解析の種類が整理できます。

• 静解析：時間変化がないので ü = u̇ = 0。残るのは K u = F、つまり連立一次方程式。
• 衝突解析：慣性項 M ü が主役。壁に当たった瞬間、構造を潰す力の正体は「減速する質量の慣性力」です。なお非線形問題では K u はより正確には内力ベクトル F^int(u) と書きます（剛性が変形とともに変わるため）。この書き換えが第4章で効いてきます。

2.4 運動量保存とエネルギー保存 —— 検証の伏線

衝突の前後で、外力（タイヤと路面の摩擦など）が無視できるなら運動量は必ず保存します。一方、運動エネルギーは保存するとは限りません。塑性変形や摩擦で熱に変わるからです。

実車の衝突は「ほぼ完全非弾性」に近く、失われた運動エネルギーこそが車体の変形仕事＝乗員を守るために吸収すべきエネルギーです。そして第9章では、シミュレーション結果の中でこのエネルギー収支が数値的に保存されているかを検証します。ここが崩れている解析は信用できません——保存則は「後で結果を疑うためのものさし」でもあるのです。

学習目標

FEMの離散化のイメージを掴み、なぜ衝突は陽解法で解くのかを自分の言葉で説明できるようになる（本サイト最大のハブ）。

前提知識

第2章（運動方程式 M ü + C u̇ + K u = F）。

読了目安

約30分

3.1 離散化とは —— 連続体を節点と要素に刻む

実物の車体は無限の点からなる連続体で、その変形は偏微分方程式で記述されます。これをそのまま解くことはできないので、FEMでは構造を要素（element）という小片に分割し、節点（node）の変位だけを未知数にします。要素内部の変位は、節点の値から補間して表す——これが離散化です。

刻んだ結果、未知数は「全節点の変位ベクトル u」という有限個の数字になり、偏微分方程式は第2章の運動方程式（連立の常微分方程式）に姿を変えます。フルカーモデルなら未知数は数百万〜数千万自由度です。

3.2 形状関数・質量行列・剛性行列のイメージ

要素内の変位を節点値からどう補間するかを決めるのが形状関数 N です。例えば2節点の棒要素なら、要素内の変位は両端の値を直線でつないだもの：

形状関数が決まると、エネルギーの釣り合い（仮想仕事の原理）から各行列が機械的に導かれます。イメージで押さえましょう。

• 剛性行列 K —— 「節点をこれだけ動かすと、どの節点にどれだけの内力が生じるか」。ばね定数の巨大な表。
• 質量行列 M —— 「各節点がどれだけの質量を背負っているか」。衝突解析では要素の質量を節点に山分けした集中質量行列（lumped mass）を使います。このとき M は対角行列になる——この一見地味な事実が、次章で陽解法の速さの源になります。
• 減衰行列 C —— 振動を減衰させる項。衝突ではエネルギー吸収の大半を塑性変形が担うため、脇役です。

3.3 陰解法（Implicit）—— 釣り合いを毎ステップ「解く」

運動方程式を時間方向に進める方法は大きく2系統あります。まず陰解法。時刻 t_n+1 の状態を求めるのに、時刻 t_n+1 自身の釣り合い式を使います：

未知数 u_n+1 が式の中に入っているので、毎ステップ大規模な連立方程式を解く必要があります。さらに非線形問題では剛性が変位に依存するため、Newton-Raphson法で「解いては修正」を収束するまで繰り返します。

• 長所：無条件安定。時間刻み Δt を大きく取れる（静解析なら1ステップでもよい）。
• 短所：1ステップが重い（行列の組み立てと分解）。接触や座屈で剛性が激変すると収束しなくなる——計算が「解けずに止まる」。

3.4 陽解法（Explicit）—— 加速度から直接「前進」する

陽解法は発想を変えます。いま分かっている時刻 t_n の情報だけで加速度を求め、次の時刻へ進みます：

• 長所：1ステップが圧倒的に軽い。収束計算が存在しないので、接触・座屈・破壊が暴れても計算は必ず前に進む。
• 短所：条件付き安定。時間刻み Δt をある臨界値（第4章のCFL条件）より小さくしないと解が発散する。衝突ではΔt ≈ 1マイクロ秒程度、100msの現象に約10万ステップが必要。

3.5 使い分けの原則

2つの方式の得意分野は歩幅の性質から自然に決まります。

判断の目安は「現象の時間」と「非線形の強さ」です。ミリ秒の現象なら、小さなΔtでも総ステップ数は現実的で、陽解法の頑健さが最大限生きます。逆に数秒〜数分の現象を陽解法で解くとステップ数が天文学的になります（準静的問題への流用テクニックは第4章・第8章で触れます）。

学習目標

中心差分法の式を追い、臨界時間刻み（CFL条件）の物理的意味を説明でき、質量スケーリングの効果と副作用を述べられる。

前提知識

第3章（陽解法と陰解法の違い）、テイラー展開の初歩。

読了目安

約30分

4.1 中心差分法 —— 差分近似で次の一歩を求める

陽解法ソルバーの標準的な時間積分が中心差分法です。速度を「半ステップずらした時刻」で定義するのがミソです。時刻 t_n で加速度を求めたら：

新しい変位 u_n+1 が決まれば、各要素のひずみ → 応力（材料モデル：第6章）→ 内力 F^int が計算でき、また①に戻ります。このループを何十万回も回すのが陽解法ソルバーの一生です。

4.2 条件付き安定性 —— なぜ歩幅に上限があるのか

中心差分法には決定的な弱点があります。Δt を大きくしすぎると、解が物理と無関係に発散するのです。直感的にはこう理解できます：構造の中を情報（応力波）は有限の速さで伝わります。1ステップの間に波が1要素分より遠くへ進んでしまうと、差分式は波を「追い越されて」しまい、つじつまの合わない答えを作り出します。数値のつじつまが破れると、誤差は1ステップごとに増幅され、数十ステップで変位が爆発します。

4.3 臨界時間刻み —— CFL（Courant）条件

重要なのは、Δt をユーザーが決めるのではなく、メッシュと材料が決めることです。ソルバーは毎ステップ全要素の Δt_cr を計算し、その最小値で全体を進めます。

4.4 「最小要素が全体を支配する」問題

CFL条件はモデル全体で最も厳しい（＝最小の）要素で決まります。100万要素のうち999,999個が5mmでも、たった1個の0.5mm要素が混じれば、全体のΔtは1/10になり、計算時間は10倍になります。

4.5 質量スケーリング —— 密度を上げてΔtを稼ぐ

微小要素を作り直せないとき、実務では質量スケーリングを使います。狙いはこうです：

密度を4倍にすれば音速が半分になり、Δt は2倍——計算時間は半分です。ソルバーには「目標Δtを下回る要素にだけ自動で質量を足す」機能があり、これが広く使われます。

ただしタダ飯ではありません。質量を足すということは慣性力 M ü を偽ることです。速く動く部位に質量を足せば、運動エネルギーも接触力も変わってしまいます。目安として：

• 追加質量は部品質量の数%以内に抑える（全体で監視）。
• 激しく運動・変形する部位（クラッシュボックスなど）への追加は特に危険。
• 準静的解析（ゆっくりした現象）では慣性の影響がもともと小さいため、比較的大胆に使える。

「どれだけ質量が足されたか」はソルバーが出力します。added mass を確認しない質量スケーリングは、結果を偽造しているのと同じ——第9章の検証項目に直結します。

学習目標

接触がペナルティ法でどう数値化されるかを理解し、接触剛性の大小が結果と時間刻みに与える影響を説明できる。

前提知識

第4章（時間刻みの制約）。ばねの力 F = k δ のイメージ。

読了目安

約25分

5.1 接触の物理 —— 数値モデルに求められる3つの仕事

FEMの要素と要素は、放っておくとすり抜けます。メッシュは互いの存在を知らないからです。「ぶつかる」という当たり前の物理を実現するために、接触アルゴリズムは毎ステップ次の仕事をします。

1. 貫入させない —— 面と面が重なったら押し戻す（法線方向の反力）。
2. 反発・分離を許す —— 押し付けている間だけ力を伝え、離れたらゼロに戻す（引っ張りは伝えない）。
3. 摩擦を伝える —— 接触面に沿った滑りに抵抗する（接線方向の力）。

第1章で見たとおり「どこがどこに触れるか」は事前に分かりません。だから接触は探索（どこが当たったか）と力の計算（どれだけ押し戻すか）の2段構えになります。まず力の計算から見ましょう。

5.2 ペナルティ法 —— 貫入をばねで罰する

ポイントは、ペナルティ法では微小な貫入を意図的に許していることです。貫入ゼロを厳密に強制する方法（ラグランジュ乗数法）もありますが、連立方程式が必要になり陽解法の軽さを殺すため、衝突解析ではばね1本で済むペナルティ法が標準になりました。「貫入は誤差だが、小さければ実用上問題ない」という割り切りです。

5.3 接触剛性の大小 —— 第4章との危険な連結

接触剛性 k はどう選ぶべきでしょうか。両極端を考えると性質が見えます。

• k が小さすぎる：貫入が深くなり、部材同士がめり込んだ嘘の変形モードになる。エネルギーの受け渡しも不正確に。
• k が大きすぎる：貫入は小さくなるが、硬いばねは高い振動数を持ちます。硬いばねの振動周期を解像するにはΔtを小さくしなければならず、第4章の安定条件が接触ばねで律速されてしまいます。Δtを下げなければ接触力が発振し、ガタガタしたノイズだらけの結果になります。

実務のソルバーは、接触する要素の剛性・寸法から「相手の要素と同程度の硬さ」になるよう k を自動計算し、ユーザーはスケール係数（既定値1.0前後）で微調整します。既定値から大きく外す前に、なぜ貫入が起きているか（メッシュ・板厚・速度）を疑うのが正しい順序です。

5.4 接触探索 —— 毎ステップ「どこが当たるか」を探す

力の計算より計算機に厳しいのが探索です。100万節点のモデルで「どの節点がどの面に触れそうか」を毎ステップ総当たりで調べると、組み合わせは10¹²のオーダー——不可能です。そこで2段階に分けます。

1. 大域探索：空間を格子（バケット）に区切り、同じ・隣の格子にいる相手だけを候補にする（バケットソート）。候補は数個〜数十個に絞られます。
2. 局所探索：候補の面に対して、節点がどの位置に投影されるか・貫入しているかを幾何計算で厳密に判定する。

大域探索は毎ステップやる必要はなく、数ステップに1回で済ませて高速化します（その間に移動しうる距離を見込んでマージンを取る）。フルカーモデルでは接触処理が全計算時間の3〜5割を占めることも珍しくなく、探索の効率はソルバーの性能そのものです。

5.5 自己接触 —— 自分が自分にぶつかる

衝突特有なのが自己接触です。クラッシュボックスが蛇腹状に座屈すると、同じ部品の面同士が折り重なって接触します。どの面とどの面が触れるか本当に予測不能なので、「モデル中の（指定した範囲の）すべての面を、すべての面に対する接触候補にする」単一面接触（single surface contact）が使われます。現代の衝突解析では、車体全体を1つの単一面接触で包むのが標準的な作法です。

5.6 摩擦 —— クーロン摩擦モデル

接線力 F_t は摩擦係数 μ × 法線力（接触反力）を超えられず、超えようとすると滑りが始まる——高校物理のクーロン摩擦がそのまま使われます（静摩擦と動摩擦を分け、滑り速度で遷移させる拡張が一般的）。摩擦で散逸したエネルギーはすべりエネルギーとして集計され、第9章のエネルギー収支の1項目になります。

学習目標

衝突特有の材料挙動（ひずみ速度依存・破壊）を理解し、Johnson-CookとCowper-Symondsの式の各項が何を表すか説明できる。

前提知識

第2章（真応力・真ひずみ）、応力−ひずみ線図の読み方。

読了目安

約30分

6.1 弾塑性の基礎 —— 降伏と硬化

金属は降伏応力 σ_y までは弾性（除荷すれば元に戻る）、それを超えると塑性変形が始まり、変形が永久に残ります。衝突でエネルギーを吸収するのは、まさにこの塑性変形です。降伏後も応力が上がり続ける現象を加工硬化と呼び、材料データでは「真応力−塑性ひずみ」の硬化曲線として与えます。よく使われる近似が指数硬化則です：

多次元の応力状態では、相当応力（von Mises応力）が降伏応力に達したかどうかで塑性化を判定します。ここまでは静的な弾塑性解析と共通です。衝突ならではの話はここからです。

6.2 ひずみ速度依存性 —— 速く変形するほど強くなる

金属をゆっくり引張ると 200 MPa で降伏する材料が、衝突のような高速変形では 300 MPa 出る——多くの金属（特に軟鋼）は変形速度が上がると降伏応力が上がる性質を持ちます。変形の速さはひずみ速度 ε̇ = dε/dt（単位 1/s = s⁻¹）で測ります。目安：

つまり衝突時の材料は、カタログの静的試験値より数十%も強いことがあります。これを無視すると、車体は実際より弱く潰れすぎる解析結果になります。逆に言えば、静的データしか入れていない衝突解析は系統的に間違う——ひずみ速度依存性は衝突解析の材料モデルに必須の要素です。

6.3 構成則の例 —— Johnson-Cook と Cowper-Symonds

ひずみ速度（と温度）の効果を取り込んだ代表的な構成則を2つ紹介します。どちらも「静的な硬化曲線に、速度に応じた倍率を掛ける」という構造です。

6.4 破壊・損傷と要素削除

変形が限界を超えると材料は破断します。FEMのメッシュは「裂ける」ことができないので、衝突解析では破壊条件を満たした要素をモデルから削除する（要素削除）のが最も一般的な扱いです。最も簡単な条件は「相当塑性ひずみが限界値 ε_f に達したら削除」。より精密には、応力状態（三軸度）によって破断ひずみが変わることを考慮した損傷モデル（Johnson-Cook破壊則、GISSMOなど）が使われます。

6.5 エネルギー吸収材 —— フォームとハニカム

金属以外にも、衝突解析には「潰れることが仕事」の材料が登場します。バンパー内のエネルギー吸収フォーム、アルミハニカム（衝突試験のバリアそのもの）、シートのウレタンなどです。これらは長い応力プラトー（ほぼ一定の応力で大きく圧縮変形する領域）を持ち、少ない反力で大きなエネルギー E = ∫σ dε を吸収します。圧縮試験の応力−ひずみ曲線をそのまま入力する専用材料モデルが各ソルバーに用意されています。密度が低く音速も遅いため、意外にも時間刻みを律速しにくいという実務上の性質も覚えておくと便利です。

学習目標

クラッシュ解析で多用される要素の特徴と、低減積分が生むアワーグラスモードの正体・対策を理解する。

前提知識

第3章（形状関数・剛性行列）、第4章（計算コストの感覚）。

読了目安

約20分

7.1 シェル要素 —— 板金モデルの主役

車体はほとんどが薄板（板厚0.6〜2mm程度）でできています。薄板を厚み方向にもソリッド要素で刻むと要素数が爆発し、しかも極薄の要素がΔtを殺します（第4章）。そこで「面としてモデル化し、板厚は要素の属性として持つ」シェル要素が主役になります。面内の伸縮＋曲げ＋横せん断を扱い、板厚方向には積分点を数点（多くは5点）置いて、曲げによる表裏の応力差を表現します。

衝突解析の標準は Belytschko-Tsay（BT）シェル。面内1点積分（後述の低減積分）で計算が非常に軽く、大変形に強い定式化です。「速いが、アワーグラス制御が前提」という性格を持ちます。精度が欲しい部位には完全積分に近いタイプ（Bathe-Dvorkin系など）を使い分けます。

7.2 低減積分とアワーグラスモード —— 偽の変形

要素の剛性・内力は要素内部の数点（積分点）で応力を評価して作ります。積分点を減らす（低減積分、4節点要素なら中心の1点だけ）と計算は数倍速くなり、しかも過剛性（せん断ロッキング）も緩和される——良いことずくめに見えます。しかし代償があります。

中心1点だけで見ていると、「積分点ではひずみがゼロなのに、要素は変形している」モードが存在してしまいます。4節点が交互に互い違いに動く、砂時計（hourglass）のような変形です。このモードは内力を生まない＝剛性ゼロなので、ひとたび励起されると抵抗なく成長し、メッシュがジグザグに崩れていきます。これがアワーグラスモード（ゼロエネルギーモード）です。

対策がアワーグラス制御です。アワーグラス型の変形パターンだけに反応する人工的な剛性または粘性の力を加え、モードの成長を抑えます。ただしこの力は物理には存在しない「作り物」なので、それがした仕事＝アワーグラスエネルギーが計算されます。制御を強くしすぎると要素が不自然に硬くなるため、「モードは抑えるが、エネルギーは小さく」が理想です。

7.3 ソリッド要素と要素品質

エンジンブロックのような塊、フォームやゴムのような厚みのある材料にはソリッド要素（8節点六面体が基本）を使います。ここでも低減積分＋アワーグラス制御の組み合わせが標準です。四面体は自動メッシュしやすい反面、1次要素は硬すぎる（ロッキング）ため、クラッシュでは六面体が好まれます。

要素の「形の良さ」も結果と Δt に直結します。チェックすべき代表指標：

• アスペクト比 —— 細長い要素は特定方向の精度が落ち、最短辺がΔtを律速する。
• ゆがみ（skew / warpage） —— 平行四辺形・ねじれ形状は積分精度を落とす。
• ヤコビアン —— 要素座標変換の健全性。負になると要素が「裏返って」おり計算不能。

静解析なら多少の悪形状は局所誤差で済みますが、陽解法では悪い要素1個がΔtと安定性を通じて全体を巻き込む——これが衝突解析でメッシュ品質にうるさい理由です。

学習目標

衝突解析モデルを組み立てる一連の流れ（メッシュ→境界条件→ソルバー設定→コスト見積り）を掴む。

前提知識

第4章（Δtとコスト）、第5章（接触）、第7章（要素品質）。

読了目安

約20分

8.0 全体の流れ

実務のモデル構築は、おおよそ次の一本道です。各工程がこれまでの章のどこに対応するかを添えます。

8.1 メッシュ作成 —— サイズは精度とコストの通貨

要素サイズの決定は衝突解析で最も重い意思決定です。理由は第4章のとおり、サイズを半分にすると要素数が（シェルで）4倍、Δtが半分、計算時間は約8倍になるからです。実務の相場観：

• フルカーの衝突モデル：基本サイズ4〜8mm程度（世代とともに細かくなってきた）。
• 変形が集中する部位（クラッシュボックス、サイドメンバー先端）：座屈の波長を拾えるよう局所的に細かく。目安として蛇腹のひと山を5要素以上で解像したい。
• 変形しない遠方（リアまわり等、前突の場合）：粗くてよい。質量と剛性が合っていれば十分。

「知りたい変形の波長より十分細かく、それ以外は思い切って粗く」——メッシュは均一である必要はなく、解きたい物理に予算を集中させるのが上手なメッシュです。

8.2 境界条件 —— 拘束・初速度の与え方

衝突解析の境界条件は静解析と発想が異なります。構造はどこにも固定されていない（車は宙に浮いた自由体）のが基本で、代わりに初期条件が主役になります。

• 初速度 —— 車両側の全節点に衝突速度（例：56 km/h = 15.6 m/s）を与える。回転を伴う場合は角速度も。
• 剛体壁・バリア —— 相手が変形しない場合は剛体壁（rigid wall）として定義。変形するバリア（アルミハニカム）は第6章の材料でモデル化。
• 重力 —— 100msの現象では落下距離は数cm以下でしばしば省略されるが、ロールオーバー等の長時間現象では必須。
• 拘束の濫用に注意 —— 実物で固定されていない場所を「計算の都合」で固定すると、そこが偽の荷重経路になります。部分モデル（コンポーネント解析）の切断面の扱いは慎重に。

8.3 ソルバー設定 —— 終了時間と出力の設計

• 終了時間 —— 前突なら120〜150ms程度。「変形が終わり、評価したい量（加速度・侵入量）が出そろうまで」。長いほど計算時間は線形に増えるだけなので、初回はやや長めが安全。
• 出力頻度 —— アニメーション用の全変形状態は1ms間隔程度（120コマ）、時刻歴（加速度・エネルギー等）は0.01〜0.1ms間隔など細かく。全状態を細かく吐くとファイルが容易に数十GBになるため、「絵は粗く、グラフは細かく」が定石。
• 質量スケーリング許容値 —— 目標Δtと追加質量の上限を設定（第4章）。
• 単位系 —— ソルバーは単位を知りません。mm-ms-kg-kN系など、一貫した単位系を最初に決め、全データで貫きます。単位系ミスは衝突解析の古典的事故です。

8.4 計算コストの見積もり

計算時間はおおよそ次の積に比例します：

（N_要素：要素数、T_終了/Δt：総ステップ数、N_CPU：並列コア数。並列化効率は100%ではないので後半は頭打ちになります。）

例えば200万要素・Δt = 0.8µs・終了時間120ms・128コア並列なら、総ステップ数は15万。1ステップあたりの処理が要素あたり〜µsオーダーとして、数時間〜十数時間のレンジに収まる——という概算が立ちます。ここでも第4章の結論が効いています：コストを決めるのは最小要素。見積もりが合わないときは、まずΔtを支配している要素を探します。

学習目標

「計算が回った」と「結果が正しい」の間にある検証プロセス、特にエネルギーバランスの読み方を身につける。

前提知識

第2章（保存則）、第4章（質量スケーリング）、第5章（接触）、第7章（アワーグラス）。

読了目安

約25分

9.1 エネルギーバランス —— 最強の健全性チェック

陽解法の結果を受け取ったら、変形アニメーションを眺める前にエネルギーの時刻歴を見ます。系の全エネルギーは、種類を変えながらも総和が保存するはずです：

健全な前面衝突解析なら、グラフはこう読めます：運動エネルギーが減った分だけ内部エネルギーが立ち上がり（エネルギーの引っ越し）、アワーグラス・接触エネルギーは床を這う程度、総和は水平な直線。

チェックの定量目安（いずれも経験則としてよく使われる値）：

• 総エネルギーのドリフト：初期値から±5%以内
• アワーグラスエネルギー：内部エネルギーの10%未満（部品別にも確認：第7章）
• 接触エネルギー：小さく、かつ大きな負値にならないこと（負に沈むのは貫入の異常＝接触の抜けのサイン）
• 質量スケーリングによる追加質量：数%以内、かつ運動部位に集中していないこと（第4章）

9.2 主要な出力量 —— 何を評価するのか

• 変形モード —— 意図した順序で潰れたか。クラッシュボックスが蛇腹に潰れる前にサイドメンバーが折れていないか。設計意図（荷重経路）との照合が第一。
• 加速度（クラッシュパルス） —— 車体（Bピラー下端など）の減速度時刻歴。乗員が受ける衝撃の入力そのもの。生波形はノイズを含むため、規定のフィルタ（SAE J211のCFC 60/180等）をかけて評価します。フィルタ条件を揃えずに波形を比較してはいけません。
• 侵入量（イントルージョン） —— 乗員室への構造の侵入。ステアリング後退量、ペダル後退量、フットウェル変形など。生存空間の確保を測る量です。
• 傷害値 —— ダミー（第10章）の頭部合成加速度から計算するHIC、胸部変位、大腿部荷重など。「構造の応答」を「人体へのリスク」に翻訳する指標です。

9.3 実験との相関 —— 何をもって「合った」と言うか

モデルの信頼性は実験との突き合わせで決まります。相関にはレベルがあります：

1. 定性的相関 —— 変形モードが同じか。潰れの順序・折れの位置を高速度カメラ映像・事後の車体と見比べる。まずここから。
2. 定量的相関 —— 加速度波形のピーク値・タイミング、侵入量、エネルギー吸収量の数値比較。波形類似度の指標（相関係数やISO/TS 18571のようなスコア）を使うこともあります。

相関が取れたモデルは、法規試験（各国の衝突安全基準）やアセスメント（NCAPなど市販車の安全性能評価プログラム）に向けた設計検討の土台になります。相関の取れていないモデルでの絶対値議論は禁物——ただし設計案AとBの相対比較なら、同じ癖を共有するぶん意味を持つことが多い、という使い分けも実務の知恵です。

9.4 よくある失敗カタログ

学習目標

これまでの各章の要素技術が、実際の開発現場でどう組み合わさるかを俯瞰し、次の学びへの道標を得る。

前提知識

第1〜9章（特に4・5・6・9章）。

読了目安

約20分

10.1 自動車の衝突モード —— 前面・側面・ロールオーバー

• 前面衝突 —— エネルギー吸収の主戦場。クラッシュボックス→フロントサイドメンバーの順に蛇腹座屈させて運動エネルギーを吸収し（第2・6章）、乗員室は変形させない、という「潰す所と守る所」の設計。フルラップ（壁全面）とオフセット（変形バリアに部分的に当てる）で問われる性能が異なります。
• 側面衝突 —— 潰し代がドア1枚分しかない最も苦しいモード。Bピラー・ドアビームの断面設計と超高張力鋼の使いこなし（第6章）、そして侵入量・侵入速度（第9章）が主戦場になります。
• ロールオーバー・ルーフ強度 —— 横転時の生存空間確保。現象時間が数秒と長く、陽解法では計算コストが厳しいため、準静的な押し込み試験（ルーフクラッシュ）で代替評価されることも多い領域です（第3章の使い分けの好例）。

どのモードでも仕事の骨格は同じです：荷重経路を設計し（第8章）、部材を意図どおり潰し（第6・7章）、エネルギー収支と実験相関で検証する（第9章）。

10.2 乗員保護 —— ダミーとエアバッグ

車体だけ守っても意味がありません。最終的な評価対象は乗員です。

• ダミー（ATD） —— 実験用の人体模型（Hybrid III、THORなど）を丸ごとFEMモデル化したものが市販されており、車体モデルに「乗せて」計算します。ダミー自体が数十万要素の精巧な接触・材料モデルの塊で、本サイトで学んだ技術の集大成です。近年は実人体の骨格・臓器まで模した人体FEモデル（THUMSなど）も使われます。
• エアバッグ —— 布（ファブリック膜要素）の袋にガスを噴き込む折り畳み展開の解析。ガスの扱いは、圧力を一様とみなす簡便法から、ガス流れを粒子や流体で解く方法（コーポレートパスクル法・ALEなど）まで段階があります。展開途中の乗員との接触（OOP：正規着座位置外）が難所。
• シートベルト —— 1次元要素とスリップリング・プリテンショナーの専用機能でモデル化。ダミー・エアバッグ・ベルトの3点セットで「拘束系（restraint system）」を構成し、車体のクラッシュパルス（第9章）を入力に傷害値を最適化します。

10.3 他分野への展開 —— 落下試験・鳥衝突

• 電子機器の落下試験 —— スマートフォンの角落ち・面落ち。ガラス・はんだ接合部の破壊（第6章）が焦点で、落下姿勢のばらつきをパラメトリックに振るのが定石。現象は数msと短く、陽解法の得意領域です。
• 鳥衝突（バードストライク） —— 鳥はほぼ水の塊として振る舞うため、通常のラグランジュ要素では変形に追従できず、SPH（粒子法）やALEで衝突体をモデル化します。「メッシュが極限まで歪む問題への拡張」という視点で眺めると、本サイトの知識の自然な延長線上にあります。
• バッテリーの圧壊・釘刺し —— 電動化時代の新しい衝突安全。セルの内部短絡をどう構造解析と連成させるかが現在進行形の研究テーマです。

10.4 これからの学びの道標

1. 手を動かす —— オープンソースのOpenRadiossなら今日から無料で始められます（→第11章で実行から結果分析まで実践します）。商用ではLS-DYNAやAltair Radiossの学生ライセンス・無償枠も。まず「箱を壁にぶつける」「筒を軸圧壊させる」だけのモデルで、Δt・接触力・エネルギー収支を自分の目で確認するのが最短の学びです。
2. 理論を深める —— 非線形有限要素法の教科書（Belytschko らの『Nonlinear Finite Elements for Continua and Structures』が定番）で、本サイトが省いた定式化（更新ラグランジュ、応力の客観性）を補う。
3. 規格・評価を知る —— 各国の衝突安全法規、NCAP系アセスメントのプロトコルを読むと、「解析で何を出力すべきか」の解像度が上がる。
4. 検証文化を身につける —— V&V（Verification & Validation）の考え方。第9章の内容を体系化した分野です。

総合ミニ演習 —— 全章横断クイズ

学習目標

オープンソース陽解法ソルバー OpenRadioss で、入力デッキの構造を理解し、解析を実行し、ログ・時刻歴・アニメーションから結果の妥当性を自分で判定できるようになる。

前提知識

第4章（Δt・質量スケーリング）、第5章（接触）、第6章（材料）、第9章（エネルギー検証）。コマンドライン操作の初歩。

読了目安

約30分（＋実習は半日）

11.1 OpenRadiossとは

OpenRadioss は、商用衝突ソルバー Altair Radioss のソースコードが2022年にオープンソース化されたものです。つまり第1章で名前を挙げた「本物の衝突解析ソルバー」を、無料で・ライセンスの制約なく・今日から動かせます。本サイトで学んだ理論——中心差分法、CFL条件、ペナルティ接触、Johnson-Cook材料、アワーグラス制御——がすべて実装されており、理論を実地で確かめる教材として最適です。

• 入手先：GitHub の OpenRadioss/OpenRadioss リポジトリ。ビルド済みバイナリがリリースページで配布されています（Linux / Windows。WindowsではWSL利用も便利）。
• 構成：Starter（モデルを読み込みチェックする前処理実行体）とEngine（時間積分を回す本体）の2段構成。後処理用の変換ツール（anim_to_vtk、th_to_csv）も同梱されています。
• サンプル：公式のモデル集（チュートリアル・検証モデル）が公開されており、最初の1本は自作せずサンプルの実行から始めるのが定石です。

11.2 入力デッキの構造 —— 理論と設定項目がつながる瞬間

OpenRadioss の入力は、/キーワード で始まるブロックを並べたテキストファイル（拡張子 .rad）です。Starterファイル（MODEL_0000.rad：モデル定義）とEngineファイル（MODEL_0001.rad：実行制御）に分かれます。構造を眺めてみましょう（書式は簡略化しています。正確な桁位置・パラメータはマニュアルを参照）。

#---- MODEL_0000.rad ─ Starter（モデル定義）----------------
/BEGIN                 ← タイトル・単位系（kg, m, s など）の宣言
/NODE                  ← 節点座標（第3章：離散化の実体）
/SHELL/1               ← シェル要素の結合情報（第7章）
/PART/1                ← 部品＝材料＋プロパティの組
/MAT/LAW2/1            ← 弾塑性材料。Johnson-Cook型 A, B, n, C（第6章）
/PROP/SHELL/1          ← 板厚・積分点数・Ishell（要素定式化、第7章）
/INTER/TYPE7/1         ← ペナルティ接触。自己接触対応（第5章）
/RWALL/1               ← 剛体壁（第8章）
/INIVEL/TRA/1          ← 初速度（第8章）
/BCS/1                 ← 拘束条件（第8章）
/END

#---- MODEL_0001.rad ─ Engine（実行制御）-------------------
/RUN/MODEL/1
  0.150                ← 終了時間 150 ms（第8章）
/DT/NODA/CST
  0.9  1.0E-6          ← 安全係数0.9・目標Δt。質量スケーリング（第4章）
/TFILE
  1.0E-5               ← 時刻歴（Tファイル）の出力間隔：細かく
/ANIM/DT
  0.0  1.0E-3          ← アニメーションの出力間隔：粗く（第8章）

本サイトの各章と設定項目の対応を一覧にします。この表が「理論→実務」の翻訳表です。

11.3 主要キーワードの設定項目リファレンス

サンプルデッキを開いたとき「この数字は何を意味するのか」が読めるように、結果を左右する主要パラメータを理論との対応付きで押さえます（全項目の正確な定義は OpenRadioss Reference Manual を参照。カード内の桁位置・省略時の既定値はバージョンで異なることがあります）。

(1) 材料：/MAT/LAW2（Johnson-Cook型 弾塑性）

実測の硬化曲線がある場合は /MAT/LAW36（PLAS_TAB）が素直です：真応力−塑性ひずみ曲線を関数（/FUNCT）として渡し、ひずみ速度ごとに複数の曲線を与えれば表引きで速度依存を表現できます。破壊をより精密に扱うときは /FAIL/JOHNSON 等の破壊カードを材料に重ねます。

(2) シェル要素：/PROP/SHELL

(3) 接触：/INTER/TYPE7

板のフチ同士が当たる問題（エッジ接触）では /INTER/TYPE11 を併用します。TYPE7は節点対面の接触なので、エッジ対エッジのすり抜けは拾えません。

(4) 実行制御と出力（Engineファイル）

11.4 実行とログの読み方 —— 理論が数字になって流れる

実行は2段階です：

# ① Starter：モデル読み込み・整合性チェック・初期貫入チェック
starter_linux64_gf -i MODEL_0000.rad -np 1
# ② Engine：時間積分の実行（-np で並列数を指定）
engine_linux64_gf -i MODEL_0001.rad -np 4

Starterの出力（MODEL_0000.out）では、部品ごとの質量、警告（初期貫入！）、エラーを確認します。WARNINGを読み飛ばさないのが第8章で学んだ試走の作法です。

Engineの出力（MODEL_0001.out）には、計算の「心電図」ともいえるサイクル表が流れます：

 CYCLE    TIME      TIME-STEP  ELEMENT        ERROR  I.ENERGY   K.ENERGY   ...  MAS.ERR
 40000  4.002E-02  8.91E-07  SHELL   3021    -1.2%  1.85E+04   6.11E+03        2.1E-03

• TIME-STEP と ELEMENT —— 現在のΔtと、それを支配している要素の番号。毎サイクル同じ要素が出続けるなら、その要素がモデル全体の計算時間を人質に取っています（第4章4.4そのもの）。
• ERROR —— エネルギー誤差（収支のずれ）。0%付近から小さな負側へゆっくり推移するのが健全（数値散逸）。正側へ増え続けるのは「エネルギーが湧いている」状態＝第9章の要調査サインです。目安±5%以内。
• I.ENERGY / K.ENERGY —— 内部エネルギーと運動エネルギー。衝突が進むと運動E→内部Eへ引っ越していく様子（第9章 図9-1）が数字で見えます。
• MAS.ERR —— 質量スケーリングによる追加質量の割合（ΔM/M）。数%以内に収まっているか毎回確認（第4章4.5）。

11.5 結果ファイルと後処理 —— グラフと動画にする

Engineは主に2種類の結果を吐き出します。

Pythonでの時刻歴プロットは、たとえば uv を使えば環境構築なしで1コマンドです：

th_to_csv MODELT01                       # T01 → CSV
uv run --with pandas --with matplotlib plot_energy.py
# plot_energy.py：CSVを読み、運動E・内部E・アワーグラスE・総和を1枚に描く
# （第9章 図9-1を、自分の解析結果で再現する）

最初に描くグラフはエネルギーバランス——第9章の結論を、そのまま自分の解析で実行します。総和が水平か、アワーグラスEが内部Eの10%未満か、接触Eが負に沈んでいないか。この3チェックが習慣になれば、あなたの後処理はもう実務レベルです。

11.6 演習：箱型ビームの軸圧壊 —— 最初の1本

公式サンプルにもある定番題材「薄肉箱型ビームの軸圧壊（クラッシュボックスの単品版）」を最初の1本に推奨します。進め方：

1. サンプルをそのまま実行し、正常終了・所要時間・出力ファイル一式を確認する。
2. ログで理論を確認：Δt支配要素はどこか（第4章）／MAS.ERRはいくらか／ERRORの推移は健全か（第9章）。
3. エネルギーバランスを描く：運動E→内部Eの引っ越し、アワーグラスEの割合（第7章の10%基準で部品別にも）。
4. ParaViewで変形モードを見る：蛇腹座屈がきれいに出ているか。折り重なる面がすり抜けていないか（第5章：自己接触）。
5. 1つだけ変えて再実行：板厚±20%、初速度×2、Ishellの変更、ひずみ速度項のON/OFF（第6章）——結果の変化を予想してから流し、予想と比べる。この「予想→実行→答え合わせ」のループが、教科書を実力に変えます。

A. 用語集

陽解法（Explicit method）

現在時刻の情報だけで加速度を求め、連立方程式を解かずに時間を前進させる解法。1ステップが軽く強非線形に頑健だが、時間刻みに上限（CFL条件）がある。→第3・4章

陰解法（Implicit method）

次の時刻の釣り合い式を連立方程式として解く解法。時間刻みを大きく取れるが、1ステップが重く、強非線形で収束困難に陥りやすい。→第3章

CFL条件（Courant条件）

陽解法の安定条件。時間刻みは「応力波が最小要素を横切る時間」以下でなければならない：Δt ≤ L_e/c。→第4章

中心差分法

速度を半ステップずらした時刻で定義する2次精度の時間積分法。陽解法ソルバーの標準。→第4章

質量スケーリング

時間刻みを支配する要素の密度を人為的に増やしてΔtを稼ぐ手法。慣性力を偽る副作用があり、追加質量の監視が必須。→第4章

ペナルティ法

接触面への貫入量δに比例した反力 F = kδ で押し戻す接触処理。微小貫入を許す代わりに連立方程式が不要。→第5章

単一面接触（自己接触）

指定範囲のすべての面同士を接触候補とする定義。座屈で自分自身に折り重なる衝突解析では事実上必須。→第5章

ひずみ速度依存性

変形速度が速いほど降伏応力が上がる材料の性質。Johnson-CookやCowper-Symondsの構成則で表現する。→第6章

要素削除（element erosion）

破壊条件を満たした要素をモデルから取り除く破断の表現法。質量・エネルギーが消えるため使用量の監視が必要。→第6章

低減積分

要素の積分点を減らして計算を高速化する手法。副作用としてアワーグラスモードを許す。→第7章

アワーグラスモード

低減積分要素で生じる、積分点にひずみを生まない砂時計状の偽変形（ゼロエネルギーモード）。人工的な制御力で抑える。→第7章

エネルギーバランス

運動・内部・アワーグラス・接触・すべり等のエネルギーの総和が保存しているかの検証。陽解法結果の健全性チェックの第一歩。→第9章

クラッシュパルス

衝突中の車体減速度の時刻歴波形。乗員保護設計への入力となる。→第9章

HIC（頭部傷害基準値）

頭部合成加速度の時刻歴から計算される傷害指標。傷害値の代表例。→第9・10章

B. 記号一覧

C. よくある質問（FAQ）

D. 参考文献・次の一歩

• 非線形FEMの理論 —— T. Belytschko, W. K. Liu, B. Moran, K. Elkhodary, Nonlinear Finite Elements for Continua and Structures（Wiley）。陽解法・大変形定式化の定番教科書。
• 陽解法・衝撃の古典 —— N. Jones, Structural Impact（Cambridge University Press）。構造衝撃力学の体系書。Cowper-Symonds則の背景もここに。
• ソルバー理論マニュアル —— LS-DYNA Theory Manual（無償公開）、OpenRadioss Theory Manual。本サイトの第4〜7章の内容が、実装レベルの詳細とともに網羅されています。理論と実務の橋渡しに最適。
• OpenRadioss —— github.com/OpenRadioss/OpenRadioss（ソルバー本体・後処理ツール・サンプルモデル）、openradioss.org（ドキュメント・コミュニティ）。第11章の実践の出発点。
• 和書 —— 非線形有限要素法・衝突安全の入門書、自動車技術会の衝突安全関連のテキスト・論文集。
• 規格・評価プロトコル —— 各国の衝突安全法規、NCAP系アセスメントの公開プロトコル、SAE J211（計測データのフィルタ処理）。
• 手を動かす環境 —— LS-DYNA、Altair Radioss などの学生ライセンス・無償枠。まず数百要素の「筒の軸圧壊」から始めるのがおすすめです（第10章）。