結果評価 — CAE用語解説

根本的に異なる物理量です。ひずみは「変形の度合い」、応力は「内部に生じる力の強さ」です。例えば、ゴムを引っ張ると大きく伸びますが（ひずみ大）、引っ張る力は小さい（応力小）です。一方、鋼材は少ししか伸びませんが（ひずみ小）、大きな力に耐えます（応力大）。降伏応力

用途が全く違います。ミーゼス応力（相当応力）は、延性材料（例えばS45C鋼）の降伏判定に使います。3次元の応力状態を1つのスカラー値に換算するための指標です。一方、主応力は最大・中間・最小の3つの値で、脆性材料（鋳鉄やガラス）の破壊や、コンクリートの設計で使われます。JIS G 3101の一般構造用圧延鋼材（SS400）の評価では、ミーゼス応力が基準になります。

絶対値だけでなく、相対変位が重要です。例えば、自動車のドアとボディの間の隙間は、ボディのたわみとドアのたわみの「差」で決まります。また、回転機械では、シャフトの「たわみ角」がベアリング寿命に直結します。Ansysの結果では、Directional Deformation（特定方向の変位）とTotal Deformation（全変位）を分けて表示するので、設計要求に合わせて見るものを選びます。

構造物に蓄えられた弾性エネルギーの総量で、単位はジュール[J]です。主に2つの用途があります。1つは「衝撃吸収」の評価。自動車のバンパーや衝撃吸収材は、このエネルギーを大きくすることで乗員保護を図ります。もう1つは「最適化」の目的関数。形状最適化で「ひずみエネルギーを最小化」すると、剛性が最大化された形状が得られます。式で書くと、体積積分で

「メッシュ収束性」を確認します。具体的には、最大応力値や注目点の応力をプロットし、メッシュサイズを半分にするごとの変化率が、例えば2%以下になった時点で収束と判断します。Abaqus/CAEには「メッシュ収束スタディ」を自動実行するスクリプトがあります。ただし、応力特異点（鋭い角）では理論上無限大に発散するので、収束しません。その場合は幾何学的なフィレットを追加するなど、モデル自体の見直しが必要です。

FEMでは、まず積分点で応力・ひずみが計算され、その後「平滑化」されて節点値が得られます。積分点の値が最も正確です。特に塑性変形や破壊解析では、積分点の履歴を追う必要があります。一方、節点値は可視化やレポート作成に便利で、特に「節点平均化」された値は滑らかなコンターになりますが、要素境界で不連続が生じる場合があります。実務では、まず積分点の結果で現象を理解し、報告書には節点平均化した結果を使うことが多いです。

結果の信頼性を直接示す指標です。非線形解析（接触や塑性）では、Newton-Raphson法の各増分ステップで平衡反復を行います。反復回数が多すぎる（例えば50回以上）か、残差が大きいまま（デフォルト許容値の10倍以上）で終了した場合、結果は物理的に不正確な可能性が高いです。Ansys Mechanicalのソルバー出力で「Force Convergence」のグラフが滑らかに0に収束しているかを必ず確認します。発散している場合は、荷重条件や接触設定を見直す必要があります。

まず「反力のチェック」です。例えば、全体に1000Nの荷重をかけたなら、固定部の反力の合計も1000N（許容誤差1%以内）になっているか確認します。次に「質量のチェック」。CAD質量とFEMモデルの質量が大きく異なっていないか。これらが合わなければ、境界条件や材料密度の設定ミスです。その後、「変形形状の確認」。常識的な変形をしているか（押したら凹むなど）。これらを通過して初めて、応力値の詳細を見る意味があります。

ほとんどの場合、そのまま使ってはいけません。特に鋭い角や穴縁の応力は、メッシュ依存性が強く、現実の部品では局所的な降伏で緩和されます。実務では「構造応力法」や「ノッチ応力法」に基づいて評価します。例えば溶接構造の疲労評価では、IIW（国際溶接学会）勧告に従い、溶接 toe から一定距離（板厚の0.4t, 1.0tなど）離れた点の応力を読み、線形補外して評価応力を算出します。単純にピーク値を見るのは危険です。

それだけではありません。まず「どの強度を基準とするか」が重要です。静荷重に対する延性材料では、降伏応力に対する安全率が一般的です。式は

最低限、以下の4点を含めます。1. **変形図**: 変形量を誇張表示した全体図。2. **重要断面の応力コンター図**: 評価対象部をクローズアップ。3. **荷重-変位曲線** (非線形解析の場合): 構造の剛性や最大耐力が分かる。4. **要約表**: 最大応力・変位の値、発生位置、使用材料の強度、計算された安全率を一覧に。背景色は白、コンターはViridisやPlasmaなどの知覚均一カラーマップを使うと、色覚特性に関わらず情報が伝わりやすくなります。

基本的な物理量は同じですが、後処理の機能性と自動化に差があります。Ansys Mechanicalは「プローブ」機能が強力で、点や線、面上の結果をリアルタイムで数値表示でき、グラフ作成も直感的です。一方、Abaqus/CAEは「Field Output」と「History Output」の概念が明確で、特定節点の時系列データの取得に優れています。また、AbaqusのPythonスクリプトによる結果抽出・処理の自動化は、大量のデータを扱う際に圧倒的に効率的です。

カップリングされた物理量を同時に評価する必要が出てきます。例えば、圧電素子の解析では、印加電圧（電気）→変形（構造）→発生応力（構造）という流れだけでなく、逆に変形から生じる分極（電気）も評価対象になります。COMSOLでは「派生値」を定義して、異なる物理場の変数を組み合わせた独自の評価指標（例えば、熱応力と機械応力の合成）を簡単に計算できます。この柔軟性が最大の特徴です。

ソルバーエンジン（CalculiXはAbaqusと同様の要素・解法を採用）の理論精度は遜色ありません。しかし、「結果評価の作業効率」と「検証機能」に大きな差があります。有償ソフトは、疲労評価ツール（Ansys nCode、FE-Safe）、自動レポート生成、規格に基づく評価（例えば船舶構造用鋼材の共通構造規格）がパッケージ化されています。無償ソフトでは、基本的な応力/変位の抽出はできますが、高度な評価はユーザー自身がスクリプトを書いて実装する必要があります。

「飛び値」と呼ばれる現象で、主に3つの原因が考えられます。1. **ひずみの要素**: 非常に細長い要素（アスペクト比>20）や内角が極端に小さい要素があると、数値誤差で異常値が発生します。2. **剛体運動の残留**: 拘束不足で微小な剛体運動が残っており、ソルバーが無理やり平衡を達成した結果。3. **接触の不安定**: 接触条件が急激に切り替わるポイント。まずは要素品質をチェックし、次に反力のバランスと変形モードを確認します。

まず「変位の拡大率」を1.0（実際の変位）に設定して確認します。可視化のため通常は100倍などに拡大しているので、食い込みはその見かけ上の可能性があります。実際の食い込みが発生する場合は、接触設定の「ペナルティ法」の剛性が低すぎる、または「ラグランジュ乗数法」でギャップ/貫入の許容値が緩すぎるのが原因です。固定面が動くのは、選択ミスやローカル座標系の設定ミスがほとんどです。節点番号を表示して、該当節点の自由度が実際に拘束されているかを一つ一つ確認します。

構造に「面外変形」などの幾何学的非線形性が顕在化している可能性が高いです。例えば薄板や細長い梁では、大きな変形により「膜応力」が発生し、荷重と応力・変位の関係が線形から外れます。これを検証するには、小さな荷重（例えば1N）と大きな荷重（1000N）をかけた時の「荷重-変位曲線」が直線かどうかを確認します。線形解析の前提が崩れているので、この場合は「大変形解析（Geometric Nonlinearity）」をONにして再計算する必要があります。

ソルバーのアルゴリズム（例えば疎行列ソルバーの前処理法）やデフォルトの収束許容値、要素公式、積分点の数が異なると、数値誤差の範囲で結果が変わります。一般的に、最大応力や最大変位で2〜5%以内の差であれば、実務上は「一致している」とみなします。ただし、疲労寿命のように対数スケールで評価する量では、小さな応力の差が寿命に与える影響は桁単位で変わるため、より厳密な一致が要求されます。ベンチマーク用の標準問題（NAFEMSなど）でツールの挙動を把握しておくことが重要です。