超弾性 — CAE用語解説

根本的な違いは、変形のメカニズムとその大きさです。金属などの線形弾性体は、結晶格子の歪みによる変形で、通常ひずみ1%未満が範囲です。一方、ゴムなどの超弾性体は、ポリマー鎖の絡み合いがほどけ、鎖自体が伸び縮みすることで大変形を生みます。ひずみ100%を超えることも珍しくありません。応力-ひずみ関係が非線形で、加えて「体積がほとんど変化しない（非圧縮性）」という強い制約があります。

その通りです。変形前後の体積比

実務で最もよく使われるのは「Mooney-Rivlinモデル」と「Ogdenモデル」です。Mooney-Rivlinの2項モデルでは、

良い着眼点です。厳密な非圧縮性は、体積ひずみが0という条件が圧力pに対する拘束方程式となり、標準的な変位法では「ロッキング」という数値的不安定を引き起こします。これを回避する主要な手法が2つ。「混合定式化」と「減積分」です。混合定式化では、変位と圧力を独立した自由度として扱い、例えばAbaqusの「C3D8H」要素（HはHybrid）のように、圧力自由度を追加して条件を弱形式で満たします。

体積ひずみ（＝変位の微分のトレース）を計算する積分点の数を、変形の計算に使う積分点より少なくすることです。例えば、Abaqusの「C3D8RH」要素は、変形の計算に8点（低減積分）、体積ひずみの計算に1点（中心）を使います。これにより、過剰な体積拘束が緩和され、ロッキングを防ぎます。ただし、今度は「アワーグラッシング」という別の不安定モードが発生する可能性があるので、安定化技術が必要です。

必須です。大きな変形と強い材料非線形性を扱うため、反復解法（Newton-Raphson法）が基本です。しかし、初期剛性が不適切だと発散します。そのため「ステップ制御」が重要。Ansysでは「Stabilization」機能、Abaqusでは「Automatic stabilization」を有効にし、不安定な接触や大変形の初期段階で仮想的な減衰力を加えて収束を助けます。また、荷重や変位を小さな増分（インクリメント）に分割し、1インクリメントあたりの最大反復回数を15〜20回に設定するのが一般的です。

載っていることは稀で、通常は実験データから「パラメータ同定（フィッティング）」を行います。最低限必要なのは「単軸引張試験」のデータです。より精度を求めるなら、「単軸圧縮」「平面引張（二軸引張り）」「せん断」のデータを追加します。実験データが異なる変形モードをカバーするほど、選択した材料モデルの予測精度が広い範囲で保証されます。例えば、シリコーンゴムのOリングを解析するなら、単軸引張と平面引張のデータを使い、Mooney-Rivlin 2項モデルでフィッティングするのが一つの定石です。

その場合は、Neo-Hookeanモデル（Mooney-Rivlinの特殊形でC01=0）か、Yeohモデルが現実的です。Yeohモデル

解析で想定される最大ひずみを少なくともカバーし、さらに余裕を持たせるべきです。例えば、シール部品で最大公称ひずみが50%と予想される場合、実験データは少なくともひずみ100%まで取得します。なぜなら、材料モデルはフィッティングに使ったデータ範囲内でのみ信頼性が高く、範囲外は外挿になり、予測が大きく外れるリスクがあるからです。また、ひずみ0%近傍の初期剛性を正確に捉えるため、ひずみ10%以下のデータ点も密に取ることが推奨されます。

基本的な材料モデル（Mooney-Rivlin, Ogden, Yeoh等）は全て備えていますが、実装の細部とユーザビリティに差があります。Abaqusは自動車タイヤ解析など、ゴム製品の産業応用に強く、超弾性要素（C3D8H, C3D10H等）や安定化アルゴリズムが非常に成熟しています。Ansys Mechanicalも同様の機能を持ちますが、材料データ入力インターフェースが直感的で、実験データを直接入力して複数のモデルで同時フィッティング比較できるツールが優れています。

COMSOL Multiphysicsは「超弾性」を単体の物理としてではなく、構造力学モジュールの一部、かつ他の物理場（熱、電気、流体）と容易に連成できる点が最大の強みです。例えば、発熱を伴うゴムの履歴現象（マリンの効果）や、電場による誘電エラストマーの大変形などを扱いやすい。材料モデルも豊富で、ユーザ定義のひずみエネルギー関数を直接入力できる柔軟性があります。ただし、純粋な大規模な構造接触問題のみの求解効率では、Abaqusに一歩譲る印象です。

可能です。例えば、Abaqusと互換性のある入力形式を扱うCalculiXは、超弾性材料モデル（HYPERELASTIC）と混合要素（C3D8H等）をサポートしています。Code_AsterもMooney-RivlinやOgdenモデルを実装。ただし、商用ソフトのような高度な自動収束制御や、グラフィカルな材料フィッティングツールはほぼありません。実験データから材料定数を求める作業は、Pythonスクリプトや別のツール（例えばMCalibrationのような専用ソフト）で事前に行う必要があります。学習コストは高いですが、原理を深く理解するには良い環境です。

最も多い原因は2つ。「不適切な初期条件」と「不安定な接触」です。超弾性体は初期状態でほとんど剛性がありません。部品同士が最初からめり込んでいたり（初期貫通）、大きな初期荷重がかかっていると、ソルバーが適切な変位増分を見つけられず発散します。まずは、初期接触状態を厳密に調整し、最初の1ステップはごく小さな変位（例えば0.1mm）または荷重を適用し、部品が安定して接触するようにします。

要素が極端に歪み、数学的に体積が負（つまり、要素が裏返しになる）と判定された状態です。大変形解析で特に低次要素（1次要素）で発生しやすい。直接的な原因は、①メッシュが粗すぎて変形を追従できない、②接触条件が急峻で要素が潰される、③材料モデルが不安定（例えばOgdenモデルのパラメータαが負で、あるひずみで応力が低下する）などです。対策は、メッシュを細かくする、要素次数を上げる（2次要素へ）、接触剛性を適切に設定する、材料モデルの安定性を確認する、です。

その通りです。これは「アワーグラッシングモード」の典型的な症状で、特に低減積分要素を使用した際に発生します。要素内部で変形モードがゼロエネルギー（ひずみエネルギーを生まない変形）を持ち、それが数値的に増幅される現象です。対策は、①完全積分要素（例えばAbaqusのC3D8HではなくC3D8）に切り替える、②「アワーグラッシング制御」を有効にする（AnsysやAbaqusのオプション）、③メッシュを細かくする（現象を軽減するが根絶は難しい）、です。完全積分要素は計算コストが上がりますが、結果の安定性は格段に向上します。

バグではなく、「材料不安定性」です。Ogdenモデルのパラメータ