後処理 — CAE用語解説

それは大きな誤解だ。後処理の本質は、膨大な数値データから「意味のある情報」を抽出し、「意思決定に必要な根拠」を可視化することだ。例えば、単に応力分布を見せるのではなく、降伏応力245MPaに対して最大ミーゼス応力がどこで320MPaに達しているのか、その領域の体積が全体の何%を占めるのかを定量的に示すことが重要になる。

分野によって変わる。構造解析なら変位、ひずみ、安全率、接触圧力。熱解析なら温度分布、熱流束、熱伝達率。流体解析（CFD）なら速度ベクトル、圧力分布、流線、渦度、乱流エネルギーだ。自動車のホイール空力解析では、回転に伴う複雑な渦構造をQクリテリオンやλ₂クリテリオンといった渦識別手法で可視化し、空気抵抗（Cd値）への寄与を評価する。

材料と破壊基準によって計算式が異なる。最もシンプルなのは、延性材料に対するミーゼス応力ベースの安全率だ。降伏応力をσ_y、相当応力をσ_eqとすると、安全率nは

それは「補間」または「平滑化」と呼ばれる処理だ。最も一般的なのは、要素内の形状関数を用いた補間だ。1次四面体要素（テトラ）の場合、要素内の任意の点での物理量φは、節点値φ_iと形状関数N_iを使って

「クリッピング」や「プローブ」機能を使う。例えば、Ansysの「Section Plane」では、ユーザーが定義した平面（例えば、x=0.1m）とメッシュ要素の交線を計算する。交線上の点では、所属する要素の形状関数から物理量を補間して値を求める。パスに沿ったグラフ（「チャート」）では、パスを細かい点群で離散化し、各点で同様の補間計算を行い、線グラフとしてプロットする。

基本的にはその通りだ。変位ベクトルuにスケール係数αを乗じて表示変位u_dispを計算する：

ダメだ。まずは「モデルが意図通りに動いたか」の基本チェックから始める。具体的には、1) 反力の確認：例えば片持ち梁なら、先端に10Nの荷重をかけた時、固定端の反力の合計が10N（ベクトル和）になっているか。2) 質量の確認：材料密度と体積から計算されるモデル総質量がCADデータと大きく乖離していないか（例：CADで1.5kgなのにFEMで0.01kgはメッシュ欠損の疑い）。これらを飛ばすと、後の結果は全て無意味になる。

ほとんどのソフトに「シンボル表示」機能がある。Ansysなら、固定支持は三角形のマーク、強制変位は矢印、分布圧力は面を覆う均一な矢印群として表示される。このシンボルの向きと大きさ（スケール）を必ず確認せよ。特に回転拘束やモーメント荷重は視認しづらいので、変位結果（微小変形でもスケールアップ表示）で挙動を確認する。例えばモーメントをかけたら、ねじりや曲げ変形が起きるはずだ。

異なるメッシュサイズ（粗い、標準、細かい）で計算した、注目する最大応力値や最大変位値を表や折れ線グラフにまとめる。横軸を要素数や節点数、縦軸を応力値とする。要素を細かくしても結果の変化が例えば2%以内に収まれば、収束したと判断できる。Abaqus/CAEの「XY Data」機能でこのグラフを作成できる。収束していない場合、応力集中部のメッシュをさらに細かくする必要がある。

最低限これだけは守れ：1) カラーバーの最大/最小値は、材料の許容値（降伏応力など）や規格の基準値と比較できるように設定せよ。例えばS45Cの降伏応力は365MPaなので、最大値を400MPaに設定するなど。2) 単位を必ず記載（MPa, mm, °C）。3) 変形スケール係数（Deformation Scale Factor: 200など）を図中に明記。これを怠ると、変形の程度を誤解される。JIS Z 8316（製図－表示の一般原則）も参考にすると良い。

哲学の違いが大きい。Ansys Mechanicalは「エンジニア向け」で、安全率、疲労寿命（nCode連携）、自動レポート生成など、設計判断に直結する機能が充実。Abaqus/CAEは「研究者・開発者向け」で、ユーザー定義フィールド出力（UVARM）やPythonスクリプトによる高度なカスタム後処理が強み。COMSOLは「マルチフィジックス統合」が売りで、異なる物理場（構造、電磁、流体）の結果を同じ画面で重ね表示し、相互の影響を直感的に評価できる。

根本的な違いは「計算機能の有無」だ。ビューア（Ensight, FieldView, ParaView）はあくまで可視化専門。計算済みのデータを読み込んで綺麗に表示するのは得意だが、新しい物理量を計算できない。一方、有償ソフトの後処理モジュール（Ansysの「Mechanical」や「CFD-Post」）は、結果データからさらに計算を行う「派生結果」機能が強い。例えば、応力テンソルから主応力やミーゼス応力をその場で計算したり、時間歴データからフーリエ変換して周波数応答を求めたりできる。

ParaViewはオープンソースでカスタマイズ性が高く、Pythonスクリプトでバッチ処理や特殊な可視化（粒子追跡の複雑なフィルタリングなど）を自動化できる。研究や新手法の開発に向く。STAR-CCM+のビューアは、自動車業界などで標準的な評価指標（例えば、ダウンフォース係数Clの時系列グラフ、ヨー角ごとの空力特性曲線）をワンクリックで出力するテンプレートが豊富で、業務効率が圧倒的に高い。定常的な業務レポート作成なら商用ビューアが有利だ。

まずは「数値的特異点」を疑え。鋭い角に集中荷重がかかっている、またはリブの付け根など鋭いエッジ部分でメッシュが粗い場合、理論上応力は無限大に発散する。この見かけ上の高応力はメッシュ依存性が強く、要素を細かくするほど値が上がり続ける。対策は、1) 現実的な面圧荷重に変更する、2) 角にフィレットを追加する、3) 結果を見る際に「スパイクを除外した平均値」や「ある距離離れた場所の値」を評価する。Ansysでは「Probe」機能で特異点から数mm離れた位置の応力を拾う。

「接触設定」の問題である可能性が高い。特に初歩的なのは、接触ペアを定義していない、または定義しても「接触アルゴリズム」が適切でない場合だ。Ansysでは「Bonded」（離れない）と設定すべきところを「Frictionless」（摩擦無しすべり）にすると貫通する。対策は、1) 接触領域を再確認し、接触・目標面を正しく設定する、2) ピンチングを防ぐため「Adjust to Touch」で初期接触を調整する、3) 非線形接触では、荷重を小さなステップ（サブステップ）で適用し、各ステップで接触状態が正しく検出されるようにする。

計算誤差ではなく「メッシュの粗さ」が原因であることが多い。壁面近傍では速度勾配が急峻なので、メッシュが粗いと速度ベクトル自体が離散的にしか求まらず、流線を積分すると不自然な軌道になる。対策は二つ。第一に、後処理側で「流線の開始点の数を増やす」「積分ステップを細かくする」設定を変えてみる。根本的解決は第二で、境界層メッシュを適切に設定する。例えば、乱流モデルを使う場合、無次元壁面距離y+を1前後にするように、最初の層の厚さを0.01mm程度に設定し、層数を15〜20層にする。

これは表示スケールの問題だ。変位や応力のコンター表示では、カラーバーの最大/最小値が「モデル全体」に基づいて自動設定される。大きな部品の変位が10mm、小さなボルトの変位が0.1mmの場合、ボルトはほぼ一色（最小値の色）で塗りつぶされ、詳細が見えなくなる。対策：1) その部品だけを選択し、「表示対象」を「選択したボディ」に変更して個別にプロットする。2) カラーバーの範囲を手動で設定する（例えば0mmから0.5mm）。Abaqusでは「オーダー制限」を解除する設定もある。