Ogden — CAE用語解説

カテゴリ: 用語集 | 2026-01-15

理論と物理

Ogdenの基本概念

🧑‍🎓

Ogdenモデルって、よく聞くけど、フックの法則やネオ・フックモデルと何が根本的に違うんですか?

🎓

根本的な違いは、ひずみエネルギーの表現形式です。フックの法則は微小変形の線形弾性、ネオ・フックは主伸びの2乗和で表される最も単純な超弾性モデルです。一方、Ogdenモデルは主伸び

$$ \lambda_i $$
のべき乗項の和でひずみエネルギー
$$ W $$
を直接定義します。式で書くと、
$$ W = \sum_{k=1}^{N} \frac{\mu_k}{\alpha_k} (\lambda_1^{\alpha_k} + \lambda_2^{\alpha_k} + \lambda_3^{\alpha_k} - 3) $$
となります。ここで
$$ \mu_k $$
$$ \alpha_k $$
は材料定数で、
$$ \alpha_k $$
は整数である必要がなく、例えばゴム材料では負の値も取り得ます。これが、広い範囲のひずみ(100%を超える伸びや圧縮)を高い精度でフィッティングできる理由です。

🧑‍🎓

「主伸びのべき乗」と言われてもイメージが湧きません。具体的に、

$$ \alpha_k $$
が例えば2.5とか-5.0とか、そんな中途半端な値になることがあるんですか?

🎓

はい、実際にあります。例えば、天然ゴムの単軸引張試験データをフィッティングすると、3項Ogdenモデル (

$$ N=3 $$
) で、
$$ \alpha_1 = 1.3, \alpha_2 = 4.0, \alpha_3 = -2.0 $$
といった値が得られることが報告されています。この柔軟性が強みで、
$$ \alpha_k=2 $$
の項だけにすればネオ・フックモデルと一致します。つまり、Ogdenモデルはネオ・フックを含むより一般的な形式と言えます。

🧑‍🎓

材料定数

$$ \mu_k $$
の単位は何ですか?また、初期せん断弾性係数
$$ \mu_0 $$
とはどう関係するんですか?

🎓

$$ \mu_k $$
の単位は応力(MPaやPa)です。重要なのは、これらの定数単体では物理的な弾性率として解釈できないことです。初期せん断弾性係数
$$ \mu_0 $$
は、ひずみがゼロの時の接線剛性で、Ogden定数から次のように計算されます:
$$ \mu_0 = \frac{1}{2} \sum_{k=1}^{N} \mu_k \alpha_k $$
。実務では、材料試験データからフィッティングして得られた
$$ \mu_k, \alpha_k $$
を使って、この式で初期剛性を確認します。例えば、
$$ \mu_0 $$
が1 MPa程度なら、比較的柔らかいゴムと言えます。

数値解法と実装

実装と安定性

🧑‍🎓

Ogdenモデルを有限要素法で実装する時、特に気をつけるべき数値的な問題はありますか?

🎓

主な課題は「体積ロック」と「非凸性」です。Ogdenモデルは非圧縮性を仮定することが多いため、混合形式(変位-圧力)の要素(例えばAbaqusのC3D8HやC3D10H)を使わないと、過剰剛性(体積ロック)が発生します。また、

$$ \alpha_k $$
が負の値を持つ項があると、主伸び
$$ \lambda_i $$
が1より小さくなる圧縮領域で、ひずみエネルギー関数が数学的に「凸」でなくなり、局所的な不安定性を引き起こすことがあります。これが収束不良の原因になります。

🧑‍🎓

「非凸性」による不安定性をソフトウェア側でどうやって検出・防止しているんですか?

🎓

主要なソルバーは、各積分点で接線剛性マトリックスの正定値性をチェックしています。例えば、Abaqus/Standardでは、材料ヤコビアン(接線剛性)が正定値でなくなると、「負の固有値」警告を出し、ステップを縮小したりします。また、実装上の工夫として、主伸び

$$ \lambda_i $$
に極小値(例えば0.01)を設けて、ゼロ除算や数値的爆発を防いでいます。ユーザー側でできる対策は、圧縮データを含まない試験データから定数を決定するか、
$$ \alpha_k $$
が負の項の影響を抑えるために、より多くの項(N=4以上)を使ってフィッティング精度を上げることです。

🧑‍🎓

ひずみエネルギー関数の第2変分、つまり材料ヤコビアン

$$ \mathbf{C} $$
を解析的に導出するのは大変じゃないですか?ソフトウェアはどうやって計算してるんですか?

🎓

その通りで、手計算は非常に煩雑です。実用的には、自動微分(Automatic Differentiation, AD)や数値微分を活用しています。例えば、COMSOL Multiphysicsの「非線形構造材料」モジュールでは、ユーザーが入力したひずみエネルギー式をシンボリックに処理し、必要な1次および2次偏微分を自動生成しています。また、Abaqusのような商用コードでは、Ogdenモデル用に最適化された閉じた形のヤコビアン式があらかじめハードコードされています。その式は、主伸びと主方向を使ったスペクトル分解に基づいていて、効率的に計算できるようになっています。

実践ガイド

材料定数の決定

🧑‍🎓

実際にゴム材料のOgden定数を決めたい時、どんな試験データが必要で、どうやってフィッティングするんですか?

🎓

最低限、単軸引張試験データは必須です。理想的には、単軸引張、等軸引張(または圧縮)、平面引張(純せん断)の3種類の試験データがあると、材料の挙動を多軸状態でよく規定できます。フィッティングは、Abaqus/CAEに付属する「材料キャリブレーション」機能や、MSC Marcの「Material Fitting」ツール、あるいは専用ソフトであるMCalibration(PolymerFEM社)を使うのが一般的です。これらのツールは、実験データと理論応力-伸び曲線の誤差(通常二乗誤差)を最小化するように、

$$ \mu_k, \alpha_k $$
を非線形最適化アルゴリズム(Levenberg-Marquardt法等)で自動計算します。

🧑‍🎓

項数

$$ N $$
は、どうやって決めればいいですか?多ければ多いほど精度が上がるんですか?

🎓

良い質問です。項数を増やしすぎると「過学習」の問題が起きます。実験データにはノイズが含まれるため、N=6も使うとノイズまで完璧にフィッティングし、見かけ上の精度は上がっても、実際の多軸状態での予測精度が落ち、数値的不安定性も増します。実務的な指針としては、まずN=2または3から始めます。ISO 17200などのゴム材料試験規格に基づくデータであれば、N=3で十分な場合がほとんどです。フィッティング後、未使用の変形モード(例えば平面引張データだけでフィッティングした場合の等軸圧縮挙動)をシミュレーションで予測し、妥当かどうかを確認する「検証」プロセスが不可欠です。

🧑‍🎓

フィッティング結果の「良さ」を判断する定量的な指標はありますか?

🎓

主に2つあります。1つは決定係数

$$ R^2 $$
で、1に近いほどフィッティング精度が高いことを示します。もう1つは、各データ点における応力の相対誤差の平均値や最大値です。実務では、
$$ R^2 > 0.99 $$
かつ、最大相対誤差が10%以内(厳しい要求なら5%以内)を目安にすることが多いです。ただし、低ひずみ域(20%以下)は剛性が設計上重要なので、この領域の誤差を特に重視します。フィッティングツールはこれらの指標を通常レポートしてくれます。

ソフトウェア比較

各ソフトウェアでの実装

🧑‍🎓

Ansys、Abaqus、COMSOLでOgdenモデルを定義する時、入力方法やオプションに違いはありますか?

🎓

大きな違いがあります。Abaqusでは、「Hyperelastic」材料モデルのサブオプションとして「Ogden」を選び、項数N(最大6)と各項の

$$ \mu_i, \alpha_i $$
を直接入力します。圧縮性を考慮する「Volumetric」部分も別途定義可能です。Ansys Mechanical APDLでは、「TB,HYPER」コマンドに「OGDEN」を指定し、同様に定数を入力します。一方、COMSOLは「超弾性材料」ノードで「Ogden」を選択しますが、その式が「次数」で表現され、N=1は1次Ogden(実質的にネオ・フック)、N=2は3次Ogdenモデル(
$$ N=3 $$
の特殊ケース)というように、用語と実装が他社と異なるので注意が必要です。

🧑‍🎓

「3次Ogdenモデル」という言葉がAbaqusとCOMSOLで意味が違うってことですか?

🎓

はい、混乱しやすい点です。Abaqusや学術文献では、「N項Ogdenモデル」のNは級数の項数を指します。一方、COMSOLでは「次数 (Order)」が1, 2, 3とあり、これはひずみエネルギー式の展開次数に関連し、COMSOLの「次数N」は、実質的には「2N-1項Ogdenモデル」に対応します。つまり、COMSOLで「次数3」を選ぶと、内部的には5項のOgdenモデルが使われます。このため、他ソフトや論文の定数値をそのままCOMSOLに入力する時は変換が必要です。マニュアルの式を確認することが必須です。

🧑‍🎓

無料やオープンソースのCAEソフト(CalculiX, Code_Asterなど)でもOgdenモデルは使えますか?

🎓

はい、主要なオープンソースソルバーもサポートしています。CalculiX(プリプロセッサはPREPO7)では、「*HYPERELASTIC」カードの「OGDEN」タイプで定義できます。Code_Asterでは、`DEFI_MATERIAU`コマンド内の`ELAS_HYPER`キーワードで`OGDEN`を指定します。ただし、商用ソフトのようにグラフィカルな材料フィッティングツールは通常付属していないため、定数は外部ツール(例えばMCalibrationの無料版やPythonスクリプト)で求めてから手動で入力する必要があります。また、収束性のチューニング(ソルバー設定)もユーザーが深く理解して行う必要があり、ハードルは少し高いです。

トラブルシューティング

よくあるエラーと対策

🧑‍🎓

Ogden材料を使ったシミュレーションで「負のヤコビアン」や「過度の変形」エラーが出て収束しません。まず何を疑えばいいですか?

🎓

まず疑うべきは「材料定数の妥当性」です。特に、圧縮やせん断を含む多軸変形が生じている領域で、材料モデルが物理的に不安定(非凸)になっている可能性が高いです。対策ステップは:1) 材料定数がフィッティングされた変形モード(通常は引張)以外で評価する。単純な要素試験(片持ち梁の曲げなど)を小変形で実行し、異常な柔らかさが出ないか確認。2) 初期剛性

$$ \mu_0 $$
が常識的な範囲か(ゴムで0.1〜10 MPa程度)確認。3) 要素タイプが適切か(非圧縮性にはハイブリッド要素必須)。4) 最初の増分ステップを極端に小さく(例えば全変位の0.1%)して、安定に初期化できるか試す。

🧑‍🎓

「体積ロック」の具体的な症状と、Ogdenモデルにおける確認方法を教えてください。

🎓

症状は、モデルが異常に硬く変形せず、反力が過大に計算されることです。特に曲げ変形で顕著です。Ogdenモデルは非圧縮性を厳密に扱うため、標準的な変位要素(AbaqusのC3D8など)を使うと確実にロックします。確認方法:Abaqusの場合、メッシュを「ハイブリッド」要素(C3D8H, C3D10H)に変えて再計算し、結果が大幅に柔らかくなるか見ます。また、要素ごとの体積ひずみ(VOLUMETRIC STRAIN)を可視化し、ほとんどゼロになっていることを確認します。ロックしている場合は、体積ひずみに不自然な変動や大きな値が現れます。

🧑‍🎓

フィッティングツールで良い精度が出たのに、実際の部品シミュレーションで応力が実験と合わない場合、材料モデル以外の原因は何が考えられますか?

🎓

材料モデル以外の主要因は3つあります。1) Mullins効果(応力ソフトニング):試験片は初回負荷時のデータですが、部品シミュレーションは往復動作をしている場合、履歴依存の軟化が無視されています。Abaqusでは「Mullins Effect」をOgdenと組み合わせて定義できます。2) ひずみ速度依存性(粘弾性):ゴムは速度が速いほど硬くなります。定常的な超弾性モデルではこれを考慮できません。3) 境界条件と摩擦:部品の接触面での摩擦係数(例えば0.5〜0.7)が現実的か。これらの影響は、単純な試験片フィッティングではキャッチできない、構造シミュレーション固有の問題です。

🧑‍🎓

Ogden定数に温度依存性を持たせたい時はどうすればいいですか?

🎓

直接「Ogdenモデルの定数が温度の関数」として定義できるソフトウェアは限られます。一般的なアプローチは2つです。1) 各温度(例えば-20°C, 23°C, 80°C)で個別に材料試験を行い、それぞれの温度で独立したOgden定数セットをフィッティングします。そして、ソフトウェアの「温度依存テーブル」機能(Abaqusでは`*ELASTIC, DEPENDENCIES=1`と`*EXPANSION`を組み合わせる等)を使って、定数セットを温度に対して補間させます。2) 物理ベースのアプローチとして、Ogdenモデルの

$$ \mu_k $$
が絶対温度Tに比例する(
$$ \mu_k \propto T $$
)と仮定し、基準温度での定数から他の温度の定数をスケーリングする方法もあります。いずれにせよ、広い温度範囲での試験データが不可欠です。

この記事の評価
ご回答ありがとうございます!
参考に
なった
もっと
詳しく
誤りを
報告
参考になった
0
もっと詳しく
0
誤りを報告
0
Written by NovaSolver Contributors
Anonymous Engineers & AI — サイトマップ
プロフィールを見る