サブモデリング — CAE用語解説

単なる分割ではありません。大規模な全体モデル（グローバルモデル）の解析結果を境界条件として利用し、関心領域だけを詳細に再解析する「局所詳細化」技術です。例えば、橋梁全体の応力分布を求めた後、ボルト穴周りの応力集中を詳細に調べたい時などに使います。

主に「強制変位」です。全体モデルで求めた、サブモデル境界上の節点変位を、サブモデル解析の強制変位条件として適用します。これがサン・ベナンの原理に基づく基本的な方法です。応力や反力も使えますが、変位が最も一般的で安定しています。

「物体の一部に加えられた力系が、それと静力学的に等価な他の力系で置き換えられても、置き換えられた場所以外の部分の応力分布には影響が少ない」という原理です。つまり、遠方からの影響は大局的な変形として捉えればよく、局所詳細部の近傍で境界条件を正確に与えれば、内部の応力は正確に求まる、という理屈です。

理論的には大きくても構いませんが、実務的には段階的に細かくすることを推奨します。例えば全体モデルの要素サイズが10mmなら、第一段階のサブモデルで2mm、さらに詳細な部位で0.5mmといった具合です。いきなり100倍細かくすると、境界での変位補間誤差が大きくなり、不自然な応力が発生するリスクがあります。

ソフトウェアが自動的に「空間補間」を行います。サブモデル境界上の各節点について、全体モデル上の近傍の節点の変位データから、形状関数を使って補間計算します。例えば、Abaqusでは「*BOUNDARY」条件に「SUBMODEL」と指定するだけで、このプロセスがバックグラウンドで実行されます。

基本的な考え方は同じです。全体モデルの要素（例えば4節点四面体要素）内にサブモデルの節点が位置すると仮定します。その要素の節点変位

その場合は「ドライビングノード」の概念が重要になります。シェル要素の中面の変位と回転を、ソリッド要素の対応する面上の節点にどのように分配するか、という問題です。Ansysでは「Shell-to-Solid Submodeling」という専用機能があり、シェルの面内変形と曲げに基づいてソリッド節点の3次元変位を決定します。手動で設定するのは非常に煩雑です。

支配方程式そのもの

標準的なワークフローは以下の通りです。
1. 全体モデル（コースメッシュ）を作成・解析し、応力集中が見られる領域を特定する。
2. サブモデルの切り出し領域を決定。経験則として、関心領域から最低でも要素2〜3層分外側に境界を設ける。
3. サブモデルを作成（詳細メッシュ、局所的な形状修正や異種材料の追加も可能）。
4. ソフトウェア機能を用いて、全体モデル結果からサブモデル境界条件を生成・適用。
5. サブモデルを解析。
6. 検証：サブモデル境界付近の応力を全体モデル結果と比較し、大きな乖離がないか確認する。

業界のベストプラクティスでは、サブモデル境界から少し内側（例えば詳細領域の深さの10%の位置）での応力（ミーゼス応力など）を比較し、差が5%以内であることが望ましいとされます。10%を超える場合は、サブモデルの切り出し領域が小さすぎる、または境界近くに大きな応力勾配がある可能性があります。その場合は切り出し領域を大きくする必要があります。

その通りです。サブモデリングの真価は、詳細化だけでなく「設計変更の検討」にあります。具体例を挙げましょう。
・全体モデル：溶接部をビード形状なしの単純な接合としてモデル化。
・サブモデル：実際の溶接ビード形状（JIS Z 3312の規格に基づく）を追加し、溶接金属と熱影響部（HAZ）に異なる材料特性（例えば、SM490A母材に対し溶接金属はヤング率を1%変化させる）を割り当て、残留応力の初期条件を導入する。
これにより、全体モデルでは評価できない局所的な疲労強度を詳細に評価できます。

通常、ソフトウェアのサブモデリング機能を利用する場合は、全体モデルの結果ファイル（Abaqusなら.odb, Ansysなら.rst）とモデルデータベースファイル（.cae, .db）の両方が必要です。ソフトウェアが境界節点の位置関係を特定し、自動補間するために、ジオメトリとメッシュの情報も参照するからです。単なる数値データのリストでは不十分です。

大きな違いは「プロセス管理」にあります。Abaqusでは、全体モデルのジョブ（Job-1）を実行した後、同じモデルデータベース（.cae）内でサブモデル部品（Part）を作成し、「Model > Edit Attributes > Submodel」で設定します。一つのCAEファイルで完結します。一方、Ansys Workbenchでは、「Submodeling」システムという専用のテンプレートを使い、全体モデルシステムとサブモデルシステムをリンクさせてプロジェクトツリーを構成します。データ連携がグラフィカルで分かりやすいです。

CalculiXでは、標準機能としての「自動サブモデリング」はありませんが、手動で実現する方法はあります。全体モデルの境界節点の変位を計算し、それをサブモデルの入力ファイル（.inp）の「*BOUNDARY」条件として記述すれば原理的には可能です。Code_Asterでは「MACR_RECAL」というオペレーターがあり、これを用いてメッシュの一部を再計算する「再局所化」が可能で、サブモデリングに近いことができます。ただし、商用ソフトのような高度な自動補間やGUIサポートは乏しいです。

可能ですが、注意が必要です。AnsysやAbaqusはいずれも非線形解析からのサブモデリングをサポートしています。しかし、サブモデル解析自体は線形解析であることがほとんどです。これは、全体モデルから受け取る「変位」が、ある特定の荷重状態（最終状態または特定の増分ステップ）での値だからです。サブモデル内で新たに接触や大きな変形などの強い非線形性が生じる場合は、その取り扱いが複雑になり、単純なサブモデリングでは不適切な場合があります。

COMSOLでは「サブモデリング」という専用機能名ではなく、「研究」の機能を組み合わせて実現します。具体的には、全体モデルを解く「研究1」を作成し、その解を「解」として保存します。次に、詳細なジオメトリで新たな「コンポーネント」を作り、「研究2」を追加します。この研究2の初期条件や境界条件として、「研究1の解」をインポートする「線形押し出し」や「一般投影」などの機能を使って変位場を適用します。柔軟ですが、設定に慣れが必要です。

最も多い原因は2つです。1つ目は、全体モデルの結果ファイル（.odb等）のパスが間違っている、またはサブモデル解析時に読み込めていないこと。2つ目は、サブモデルの境界面（カットバウンダリ）が、全体モデルのジオメトリ/メッシュからはみ出していることです。ソフトウェアは、サブモデルの各境界節点が、全体モデルの「どの要素内にあるか」を検索します。境界が全体モデルの外にあると、補間のための「ドライビングノード」データが見つからず、このエラーになります。まずはサブモデルの切り出し領域を確認してください。

これは典型的な「境界効果」です。主な原因は3つ考えられます。
1. サン・ベナン原理の適用限界：境界が関心領域に近すぎる。対策は切り出し領域を大きくする。
2. 補間誤差：全体モデルのメッシュが粗すぎて、詳細な変形モードを捉えられていない。全体モデルのメッシュをある程度細かくするか、境界から離れた位置で結果を評価する。
3. 不連続な境界条件：Abaqusなどで「*SUBMODEL, INTERPOLATE=ELCENTRIC」など、節点ではなく要素中心のデータを補間に使う設定を試すと平滑化される場合があります。

動解析からのサブモデリングは「時間歴サブモデリング」と呼ばれ、より高度な設定が必要です。単に変位だけでなく、速度や加速度も境界条件として考慮する必要がある場合があります。Ansysでは「Submodeling」システムの詳細設定で「Load step and time」を指定します。また、全体モデルが動的である場合、サブモデルも動的解析として設定しなければ、慣性効果を見逃す可能性があります。この場合、サブモデルは「モーダル過渡応答解析」などを行い、全体モデルから得た変位時歴を強制変位として与える形になります。ソフトウェアのマニュアルで「transient submodeling」を確認してください。

これは重要なポイントです。サブモデリングの基本形（変位境界条件）では、境界は「強制変位」で固定されているので、内部に新たな荷重を加えても境界は動きません。これは、サブモデルが全体モデルの一部として「拘束されている」状態を表しています。もし、局部荷重によって境界の反力が大きく変化し、それが全体モデルの変形に影響を与えるほどであれば、サブモデリングの前提が崩れています。その場合は、局部荷重を全体モデルに反映させて再解析するか、またはサブモデリングではなく「子モデル法」など、境界の反力を全体モデルに戻して連成させる高度な手法を検討する必要があります。