シミュレーション — CAE用語解説

それは大きな誤解だ。メッシュは計算の「基礎グリッド」であり、品質が結果を左右する。単に細かいだけではダメで、例えば応力集中部では細かく、変形の少ない部分では粗くする「適応メッシュ」が基本だ。Ansysのマニュアルでは、曲率や狭い隙間に自動で細かいメッシュを生成する「Curvature」や「Proximity」機能が推奨されている。

分野によって異なる。構造解析なら、弾性体の変位と応力の関係を表す

収束とは、反復計算を繰り返すうちに解の変化量が許容誤差（例えば、残差ノルムが1e-6以下）に収まることを指す。終わることではない。非線形解析（接触や塑性）では、荷重を小さなステップに分割して解く「増分解析」を行うが、あるステップで反復計算が50回を超えても収束しない場合、「収束が悪い」と判断する。これは接触条件の急激な変化や材料の不安定な挙動が原因だ。

その通り、現実の拘束や環境をモデル化するものだ。例えば「対称条件」は、モデルが対称であればその面の法線方向変位を拘束し、面内の回転を拘束する。これで計算コストを1/2や1/4に削減できる。また、「遠方境界条件」はCFDでよく使い、流入速度や圧力出口を設定する。Ansys Fluentでは「velocity-inlet」や「pressure-outlet」として指定する。間違えると非物理的な流れが発生するので注意が必要だ。

絶対に気にする必要がある。直接法（例：PARDISO, MUMPS）はメモリを大量に消費するが、数値的に頑健で、接触や不連続を含む非線形問題の各ステップで有効だ。反復法（例：CG法、GMRES）は大規模問題でメモリ効率が良いが、条件数が悪い（剛性マトリクスが悪条件）と収束しない。Abaqus/Standardでは、デフォルトで直接法を使うが、大規模線形弾性問題には反復法ソルバーをオプションで選択できる。メモリ32GBのPCで100万節点の問題を解く場合、直接法ではメモリ不足になる可能性が高い。

一次要素（線形要素）は角点のみに節点があり、二次要素は辺の中点にも節点がある。曲げや応力集中を表現するには二次要素の方が精度が圧倒的に高い。しかし、計算コストも高い。目安として、同じ精度を得るために必要な一次要素の節点数は、二次要素の約8倍（3次元の場合）と言われる。ただし、接触解析では二次要素は接触面の処理が複雑になり、収束性が悪化することがある。そのため、AbaqusのExplicit（動的解析）では、一次要素がデフォルトで推奨されることが多い。

離散化誤差は、連続体を有限個の要素で近似することによる本質的な誤差。メッシュを細かくすれば減らせる。丸め誤差は、コンピュータが浮動小数点数を有限桁で表現することによる誤差だ。通常、離散化誤差の方が桁違いに大きい。例えば、片持ち梁の先端たわみを一次要素で計算すると、理論値に対して10%以上の誤差が出ることもある。一方、丸め誤差は倍精度計算では約1e-16のオーダーだ。ただし、剛性マトリクスの条件数が極端に悪い（例えば1e12以上）と、丸め誤差が増幅されて解が意味をなさなくなる「数値的不安定」が起きる。

時間刻みは重要な違いだが、本質は解法の安定性だ。陽解法（Abaqus/Explicit, LS-DYNA）は条件付き安定で、時間刻み

最初から全域を均一に細かくするのは非効率だ。実務では、まず「コースメッシュ」で全体の挙動を把握し、応力の高い領域や関心領域を特定する。次に、その領域のみメッシュを細かくした「ローカルリファインメント」を行う。Ansys Workbenchの「Convergence」ツールのように、関心出力（最大応力など）の変化が例えば2%以内に収まるまで、要素サイズを自動で1/2ずつ細かくする機能を使うと良い。最終的には、中程度のメッシュと最も細かいメッシュの結果を比較し、差が許容範囲内であることを確認する。これで2〜3回の計算で済む。

実験データをそのまま使うと、数値的な振動を引き起こすことがある。ソルバーはデータポイント間を線形補間するため、入力ポイントが多すぎると勾配が不連続に変化し、収束性が悪化する。特に降伏点直後はデータを多めに、その後はなめらかに変化するように間引くのがコツだ。また、ひずみ硬化則（例：Swift則

完全一致は現実的に不可能で、目的に応じた許容誤差を設定する。自動車業界の構造耐久性解析では、応力集中部の最大応力で±10%、全体の変形モードや固有振動数で±5%以内を目安とするケースが多い。ただし、これは材料定数、境界条件、メッシュを適切に設定した上での話だ。実験自体にも誤差（ひずみゲージの貼り付け誤差、荷重セルの精度など）がある。重要なのは、トレンド（パラメータを変えた時の変化傾向）が一致しているかどうかだ。規格として、ASME V&V 20やJSME S 012では、検証と妥当性確認のガイドラインが定められている。

最低限、以下の5項目は必須だ：1. 使用ソフトウェアとバージョン（例：Ansys Mechanical 2023 R2）、2. 要素タイプと総節点数/要素数（例：二次四面体要素（SOLID187）、節点数約50万）、3. 材料モデルと定数（例：弾性：E=210GPa, ν=0.3、塑性：双線形硬化、降伏応力355MPa）、4. 境界条件と荷重条件の詳細（拘束位置、荷重値と方向）、5. 収束判定基準（例：力の残差ノルム 1.0e-5）。これがないと、第三者が結果を再現・評価できない。自動車業界では、さらにFMEA（故障モード影響解析）の観点から、想定外の荷重ケースについての考察も求められる。

コアとなる数値解法の哲学が異なる。Abaqusは非線形解析に非常に強く、特に「接触」と「材料非線形」の組み合わせが得意だ。そのため自動車のゴムブッシュやボルト締結部の解析でシェアが高い。一方、Ansys Mechanicalは「ワークフローの統合性」と「マルチフィジックス連成」に強みがある。Ansys Workbench上で流体（Fluent）と構造（Mechanical）を容易に連成させ、熱応力や流体構造連成（FSI）を解ける。また、HPC（高性能計算）ライセンスの柔軟性もAnsysが優れていることが多い。予算が許せば、用途によって使い分けるか、両方使える環境が理想だ。

COMSOLの核は「弱形式」でユーザーが偏微分方程式（PDE）を直接入力できる点だ。AnsysやAbaqusが「物理インターフェース」（構造、流体、電磁気など）を提供するのに対し、COMSOLはそれらを「PDEの組み合わせ」として実装している。これにより、学術的な新しい物理モデルや、複雑に結合した現象（例えば、電界で誘起される化学反応とそれによる熱変形）を一つの環境でモデリングできる。ただし、大規模な3D構造解析や複雑な接触問題では、Abaqusのような特化型ソルバーに比べて計算効率が落ちる傾向がある。メッシュ生成機能も、AnsysやAltair HyperMeshに比べると素直ではないという評価だ。

特定の分野では実用レベルに達している。CalculiXはAbaqusの入力ファイル形式（.inp）と互換性があり、線形・非線形構造解析が可能だ。ただし、GUIが貧弱で、前処理・後処理には別のツール（PrePoMax, Gmsh, Paraview）が必要になる。Code_Asterはフランス電力公社（EDF）が開発したもので、破壊力学や地盤解析に強い機能を持つ。両者とも有償ソフトのような高度なコンタクトアルゴリズムやソルバーの最適化は不足しており、計算時間が長くなりがちだ。しかし、教育目的や、自社でカスタムソルバーを開発する際のベースとしては極めて有用。商用ソフトのライセンス費用が年間数百万円かかることを考えれば、コストパフォーマンスは高い。

特定の問題に対して、アルゴリズムが極限までチューニングされている点だ。MSC Nastranは航空宇宙業界で長年使われており、線形静解析・固有値解析（振動・座屈）の速度と精度、特に大規模モデルを扱う能力に定評がある。その入力形式（Bulk Data Deck）は業界のデファクトスタンダードだ。LS-DYNAは陽解法のエキスパートで、衝突・爆発・プレス成形などの極限の非線形動的問題において、接触アルゴリズムと材料モデルの豊富さで群を抜いている。これらの専用ソルバーは、汎用ソルバーが苦手とする「悪条件問題」を解くための特殊な数値的処方箋を多く内蔵している。ただし、ライセンスは高額で、習得にも時間がかかる。

まず、ログファイルで「どの自由度（U1, U2, U3, ROTなど）が収束していないか」を確認せよ。変位が収束しないなら境界条件が不適切（剛体移動が残っている）かもしれない。次に、最初の増分ステップで収束しない場合、初期状態に問題がある。接触条件が「突然」適用されていないか確認し、接触面には微小な初期ギャップまたはオーバーラップを設ける。Abaqusでは「スタビライゼーション」や「自動ステップ」をオンにし、Ansysでは「弱いばね」を導入して初期の数値的不安定性を抑える。荷重が大きすぎる場合は、荷重を1/10にして試し、段階的に増やす「リダクション」も有効だ。

要素が極端に変形し、数学的にマッピングが成立しなくなった状態だ。主な原因は3つ：1. メッシュ品質が悪い（スキュー角が大きい、アスペクト比が大きい）、2. 材料モデルが不適切（例えば、大きなひずみまで弾性のみで塑性を考慮していない）、3. 荷重/変位条件が現実的でない。対策は、まずメッシュを改善する