Ogden — CAE用語解説
理論と物理
Ogdenの基本概念
Ogdenモデルって、よく聞くけど、フックの法則やネオ・フックモデルと何が根本的に違うんですか?
根本的な違いは、ひずみエネルギーの表現形式です。フックの法則は微小変形の線形弾性、ネオ・フックは主伸びの2乗和で表される最も単純な超弾性モデルです。一方、Ogdenモデルは主伸び
「主伸びのべき乗」と言われてもイメージが湧きません。具体的に、
はい、実際にあります。例えば、天然ゴムの単軸引張試験データをフィッティングすると、3項Ogdenモデル (
材料定数
数値解法と実装
実装と安定性
Ogdenモデルを有限要素法で実装する時、特に気をつけるべき数値的な問題はありますか?
主な課題は「体積ロック」と「非凸性」です。Ogdenモデルは非圧縮性を仮定することが多いため、混合形式(変位-圧力)の要素(例えばAbaqusのC3D8HやC3D10H)を使わないと、過剰剛性(体積ロック)が発生します。また、
「非凸性」による不安定性をソフトウェア側でどうやって検出・防止しているんですか?
主要なソルバーは、各積分点で接線剛性マトリックスの正定値性をチェックしています。例えば、Abaqus/Standardでは、材料ヤコビアン(接線剛性)が正定値でなくなると、「負の固有値」警告を出し、ステップを縮小したりします。また、実装上の工夫として、主伸び
ひずみエネルギー関数の第2変分、つまり材料ヤコビアン
その通りで、手計算は非常に煩雑です。実用的には、自動微分(Automatic Differentiation, AD)や数値微分を活用しています。例えば、COMSOL Multiphysicsの「非線形構造材料」モジュールでは、ユーザーが入力したひずみエネルギー式をシンボリックに処理し、必要な1次および2次偏微分を自動生成しています。また、Abaqusのような商用コードでは、Ogdenモデル用に最適化された閉じた形のヤコビアン式があらかじめハードコードされています。その式は、主伸びと主方向を使ったスペクトル分解に基づいていて、効率的に計算できるようになっています。
実践ガイド
材料定数の決定
実際にゴム材料のOgden定数を決めたい時、どんな試験データが必要で、どうやってフィッティングするんですか?
最低限、単軸引張試験データは必須です。理想的には、単軸引張、等軸引張(または圧縮)、平面引張(純せん断)の3種類の試験データがあると、材料の挙動を多軸状態でよく規定できます。フィッティングは、Abaqus/CAEに付属する「材料キャリブレーション」機能や、MSC Marcの「Material Fitting」ツール、あるいは専用ソフトであるMCalibration(PolymerFEM社)を使うのが一般的です。これらのツールは、実験データと理論応力-伸び曲線の誤差(通常二乗誤差)を最小化するように、
項数
良い質問です。項数を増やしすぎると「過学習」の問題が起きます。実験データにはノイズが含まれるため、N=6も使うとノイズまで完璧にフィッティングし、見かけ上の精度は上がっても、実際の多軸状態での予測精度が落ち、数値的不安定性も増します。実務的な指針としては、まずN=2または3から始めます。ISO 17200などのゴム材料試験規格に基づくデータであれば、N=3で十分な場合がほとんどです。フィッティング後、未使用の変形モード(例えば平面引張データだけでフィッティングした場合の等軸圧縮挙動)をシミュレーションで予測し、妥当かどうかを確認する「検証」プロセスが不可欠です。
フィッティング結果の「良さ」を判断する定量的な指標はありますか?
主に2つあります。1つは決定係数
ソフトウェア比較
各ソフトウェアでの実装
Ansys、Abaqus、COMSOLでOgdenモデルを定義する時、入力方法やオプションに違いはありますか?
大きな違いがあります。Abaqusでは、「Hyperelastic」材料モデルのサブオプションとして「Ogden」を選び、項数N(最大6)と各項の
「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) 初期剛性
「体積ロック」の具体的な症状と、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モデルの
関連トピック
なった
詳しく
報告