前処理 — CAE用語解説

いい質問だ。CAEにおける「前処理」は、ソルバーが計算を実行できるように幾何形状を離散化し、物理条件を定義するまでの一連の作業全般を指す。具体的には、幾何モデルの修復、メッシュ生成、材料物性の割り当て、境界条件・荷重条件の設定、接触定義、解析タイプの選択などが含まれる。例えば、Ansys Mechanicalで言えば、GeometryからModelまでのワークフローがほぼ前処理に相当する。

実務ではそれが一番の落とし穴だ。設計用のCADデータは、製造や外観を重視しており、解析に不要な微小な面（例えばR0.1mmのフィレット）や、完全な接合ではない「接触状態」の部品が含まれる。ソルバーはこれらの情報をそのまま解釈できない。例えば、厚さ0.5mmのリブと本体が0.01mmの隙間でモデリングされていたら、剛性が過小評価され、実体とは全く異なる振動モードが計算されてしまう。前処理は、こうした「設計モデル」を「解析モデル」に変換する工程なんだ。

経験則として、応力集中を評価しない部位では、応力評価対象部の板厚の10%以下のフィレットや穴は除去することが多い。また、固有値解析では質量や剛性への寄与が小さい（全体の1%未満）部品は除外する。重要なのは、その判断が解析目的に依存することだ。疲労寿命を評価する場合は、微小な切り欠きも全て再現する必要がある。Abaqus/CAEの「Geometry Edit」ツールや、ANSYS SpaceClaimの「Pull」機能は、こうした形状の簡略化（デフェーチャリング）に特化している。

主要なアルゴリズムは2つある。1つは「Advancing Front法」で、境界から内部に向かって三角形や四角形を順次生成していく方法だ。もう1つは「Delaunay三角形分割」で、まず点群を生成し、それらをDelaunay条件（三角形の外接円内に他の点が含まれない）を満たすように接続する。実用的なソルバーはこれらを組み合わせている。例えば、Altair HyperMeshの2DメッシャーはAdvancing Frontが基幹だが、メッシュ品質を改善するためにDelaunay再分割も部分的に利用する。

数学的な本質は、任意の3次元領域を「位相的に直方体」にマッピングできるかどうか、だ。複雑な形状は、このマッピングが非常に難しくなる。そのため、六面体メッシュ生成では「マルチブロック法」がよく使われる。これは、形状をいくつかの六面体ブロックに分割し、各ブロック内で構造化メッシュを生成する方法だ。Siemens NX NastranやMSC Apexの一部機能はこれを自動化しようとしているが、依然としてユーザーによるブロック分割（デコンポジション）が必要な場合が多い。一方、四面体は任意の形状を埋め尽くせる単体なので、アルゴリズムが比較的容易なんだ。

数値的な誤差、特に「剛性マトリクスの条件数」を悪化させる。極端なアスペクト比（例えば100:1）の要素があると、要素剛性マトリクスの最大固有値と最小固有値の比が大きくなり、連立一次方程式のソルバーが解を求める際に数値的不安定性を引き起こす。具体的には、反復ソルバー（PCG法など）の収束が遅くなったり、直接ソルバー（スパースLU分解）でも丸め誤差が増幅される。許容基準は解析タイプによるが、静的弾性解析ではアスペクト比10以下、スキュー角60度以下を目標とするのが一般的だ。Ansysのメッシュメトリクスチェックはこれらの値を自動で評価してくれる。

確立されたベストプラクティスがある。まず「解析目的の明確化」だ。何を、どの精度で知りたいのか。次に「幾何モデルの評価と修復」。CADからのインポートエラー（隙間、重複面、微小エッジ）をチェックし、解析に不要な特徴を除去する。その後「シミュレーション方針の決定」：対称性の利用、線形/非線形の選択、接触の定義方針だ。これが済んでから初めて「メッシュ生成」に移る。メッシュを切ってから「あ、この部品は剛体で良かった」と気付くと、全てやり直しになるからな。

これは物理現象と計算コストのトレードオフだ。「ボンデッド」は最も単純で、接触面を完全に接着（変位連続）とみなす。はめあい部の大まかな応力を見る初期検討向けだ。「フリクション」は接触・離脱・滑りを許し、摩擦を考慮する。ボルト締結部の荷重伝達や、ブレーキパッドの解析に必須だが、計算は非線形になり収束性の問題が生じる。「フリーフリクション」は摩擦を無視した接触で、離脱と滑りのみを考慮する。多くの非線形静解析でバランスの良い選択肢だ。Ansysでは「Frictional」の摩擦係数を0に設定すれば「Frictionless」と同じになる。

それは最も危険な考え方だ。無闇に細かくすると、要素数が膨大になり計算時間が現実的でなくなるだけでなく、丸め誤差の蓄積でかえって精度が落ちる場合さえある。実務では「収束性検討」を行う。例えば、応力集中部の最大応力値に注目し、メッシュサイズを段階的に細かく（例えば2mm, 1mm, 0.5mm）して計算する。結果の変化が許容範囲内（例えば5%以内）になった時点のメッシュサイズを採用する。このプロセスは、Abaqusの「Mesh Convergence Study」機能や、自前のPythonスクリプトで自動化できる。

用途による。統合環境（Ansys Workbench, SIMULIA Abaqus/CAE, COMSOL Multiphysics）は、ワークフローの連携がスムーズで学習コストが比較的低い。特にパラメータスタディや最適化と連携させる場合に強い。一方、専用ツール（Altair HyperMesh, ANSA）は、超複雑なアセンブリモデル（自動車のフルボディなど）や、特定の業界規格（NASTRAN, LS-DYNAの入力デッキ）に特化した高度なメッシュ制御機能を持つ。航空宇宙業界では、数百万要素のモデルを扱うため、HyperMeshのバッチメッシング機能が重宝される。

設計思想が全く異なる。CATIAやSolidWorksは「履歴ベースのパラメトリックCAD」で、フィーチャーの作成順序や寸法拘束が木構造で管理されている。これは設計変更には強いが、他ソフトからのインポートモデルを編集するには不向きだ。一方、SpaceClaimは「ダイレクトモデリング」が中心で、面やエッジを直接押し引き（Pull）して形状を変更する。履歴に縛られないため、どんなCAD由来のモデルでも素早く修正（デフェーチャリング、中実化→中空化など）できる。CAE前処理に特化して「形状をいじる」ことに焦点を当てたツールと言える。

可能だが、対象が限られる。Gmshはスクリプト（.geoファイル）によるメッシュ生成に優れ、研究用途や一定の形状パターンが反復するモデルでは非常に強力だ。しかし、GUI操作による複雑なアセンブリモデルのインタラクティブな修復機能は、有償ツールに比べて貧弱だ。Salomeは統合プラットフォームを目指しており、幾何・メッシュ・後処理まで一貫しているが、商用ソフトのような高度な接触自動定義やメッシュ品質のバッチチェック機能は期待できない。コストを抑えつつ、特定の反復作業を自動化したい場合や、教育用途には十分有用だ。

ほぼ100%メッシュの品質問題だ。特に、六面体要素や2次要素（中間節点あり）で発生しやすい。「負のヤコビアン」は、要素の形状が極端に歪んでおり、等パラメータ座標から実座標への写像が一意でなくなっている（要素が自分自身と交差している）状態を示す。具体的には、凹んだ六面体要素や、中間節点が辺から大きく飛び出した「バウンド」要素が原因だ。対策は、メッシュ品質チェックで「ワープ角」や「スキュー」を確認し、該当要素を再メッシュする。Abaqusでは「メッシュ→検証」ツールで、ヤコビアンチェックの閾値を緩和（「許容範囲」を増やす）する暫定対策もあるが、根本解決にはならない。

まず「初期接触状態」を疑え。部品同士が最初からめり込んでいないか？ 微小な隙間が意図せず空いていないか？ 例えば、Ansys Mechanicalでは「Contact Tool」で初期接触状態を可視化できる。めり込みは「剛体運動」を引き起こし、ソルバーが発散する。隙間は接触が突然発生する「不連続」を生み、収束を難しくする。対策は、CAD段階で適切なクリアランスを設定するか、解析ソフトの「接触面オフセット」機能で調整する。次に「メッシュ不一致」だ。接触面同士のメッシュサイズが違いすぎると、接触力の伝達が不安定になる。接触面のメッシュサイズは同程度に揃えるのが原則だ。

64ビット版ソフトウェアを使用していることが大前提だ。その上で、ワークフローを分割する。例えば、HyperMeshでは「フェザーライト」モードを利用し、表示していないコンポーネントのデータをメモリからアンロードする。幾何修復とメッシュ生成を別セッションで行うのも手だ。また、アセンブリ全体を一度にメッシュせず、部品ごと（サブシステムごと）にメッシュを生成し、後に結合する「サブストラクチャリング」手法を検討する。物理メモリが足りない場合は、SSDを搭載したマシンで仮想メモリのページファイルをSSD上に設定し、容量を確保する（速度は落ちる）。根本的には、ワークステーションのRAMは可能な限り大容量（現在では128GB以上が望ましい）にすべきだ。