FEM/CFD解析結果の正しい見方 — コンター図の罠から報告書作成まで
なぜ「結果の見方」が重要なのか
コンター図の罠 — 自動スケールの落とし穴
先生、解析結果の見方ってコンター図を見るだけじゃダメなんですか? きれいな色が出てると「できた!」って思っちゃうんですけど…
それ、実務で一番ハマるやつだね。コンター図のカラーバーの最大値が自動スケールだと、応力特異点に引っ張られて全体が真っ青になって、本当に見たい情報が埋もれるんだ。例えばボルト穴のエッジに鋭角があると、そこの応力はメッシュを細かくするたびに際限なく上がる——いわゆる応力特異点。自動スケールだとその特異点の値がカラーバーの上限になるから、本来見たい「フランジ面全体の応力分布」が全部同じ色になって判別できなくなる。カラーバーの最大値がどこで・なぜ出ているかを確認するのが、コンター図を読む第一歩だよ。
正しい確認の順序——変形図→収支→プローブ
いきなりコンター図に飛びつかず、次の順で確認するのが体系的な手順です。
- 変形図(適切な倍率で) — モデルが「物理的にありえる動き」をしているか。拘束面が動いていないか、部品が食い込んでいないか、荷重方向は正しいか。結果解釈の誤りの大半はこの段階で発見できます
- 大域収支 — 反力総和=入力荷重(構造)、質量・エネルギー収支(熱流体)。境界条件検証の反力チェックと同じ規律です
- ポイントプローブ — 理論値・手計算・過去実績と比較できる代表点の数値を確認。公称応力が梁理論と桁で合っているか
- コンター図 — ここで初めて分布を見る。スケールは目的に合わせて手動設定
「見る」とは意思決定の材料を作ること
解析結果の解釈は鑑賞ではなく、「この設計で成立するか」「どこを変えるべきか」という意思決定への変換作業です。したがって評価基準(許容応力・目標圧損・温度上限)を先に明文化し、結果をその基準と比較できる形(数値・マージン・分布)に整えることが「正しい見方」の定義になります。基準なしにコンター図を眺めるのは、ゴールの位置を知らずにボールを蹴っているのと同じです。
数値の裏側——その色は何の値か
応力値ができるまで——積分点・外挿・節点平均
FEMの応力は節点で直接計算されていません。実際のパイプラインは「①要素内の積分点で応力を計算 → ②要素の節点位置へ外挿 → ③複数要素が共有する節点で平均化」です。ポストプロセッサの平均(averaged)表示は③まで、非平均(unaveraged)表示は②までの値です。ここから2つの実務規則が導かれます。
- 平均・非平均の差はメッシュ品質の指標 — 十分細かいメッシュでは両者は近づく。差が大きい(目安5〜10%超)の領域は要素間の応力ジャンプが大きく、メッシュ不足のサイン
- 材料・板厚の境界をまたぐ平均は無効 — 異材界面では応力は本来不連続。境界をまたいだ節点平均は物理に存在しない中間値を作り出すので、部品・材料ごとに平均範囲を区切る
相当応力にはもう1つ罠があります。von Mises応力は成分の非線形関数なので、「成分を平均してからvon Misesを計算」と「各点のvon Misesを平均」は一致しません。ツールがどちらを表示しているかを確認せずにピーク値を議論すると、数%〜十数%の差が「謎の不一致」として現れます。
CFDの表示値の裏側——セル値・補間・壁面量
有限体積のCFDでは変数はセル中心にあり、コンター表示は節点への補間を経ています。特に注意すべきは壁面量(摩擦係数・熱伝達率・y+)で、これらは壁面第1セルの解と壁処理(壁関数か低Re型か)から再構成された量です。同じ流れ場でもツール・設定によって壁面熱流束の算出法が異なるため、ツール間比較では「どの定義の熱伝達率か」(基準温度の取り方)まで揃える必要があります。また断面平均値も「面積平均か質量流量平均か」で値が変わります——圧損評価は通常、質量流量平均の全圧で行います。
変形図・ベクトル・主応力の読み方
変形図の表示倍率は自動で数百倍になっていることがあり、倍率を確認せずに「大変形だ」と誤読する事故が定番です。倍率1で見る癖をつけ、幾何非線形の要否は変形量と板厚・スパンの比で判断します。応力はスカラーの相当応力だけでなく、主応力ベクトル(方向付き)を見ると荷重パスが読めます。圧縮側の座屈リスク、引張主応力と溶接線の直交など、方向情報でしか判断できない設計論点は多く、コンター+主応力ベクトルのセットを標準にする価値があります。
報告書作成の実務
伝わる図の作法
| 項目 | 作法 | 理由 |
|---|---|---|
| カラーバー | 手動で範囲固定。比較図は全図で同一スケール | 自動スケールの図同士は比較不能 |
| 上限の設定 | 許容値や降伏応力に合わせ、超過域を1色で飽和表示 | 「どこがNGか」が一目で伝わる |
| 特異点 | 評価対象外なら注記し、スケールから除外 | 特異値でスケールが潰れるのを防ぐ |
| パスプロット | 評価線に沿った1Dグラフを併記 | コンターの色よりも定量的な議論ができる |
| 変形図 | 倍率を図中に明記 | 誇張表示の誤解を防ぐ |
報告書に載せるべき検証情報
結論の図表に加えて、第三者が信頼性を判断できる検証セクションを設けます。最小構成は、①メッシュ情報と収束性確認の結果(GCIまたは水準比較)、②反力・収支チェックの数値、③収束履歴の要約(全増分収束・最終残差)、④モデル化の主要仮定(境界条件の理想化・材料モデル・省略した形状)とその影響の見積り、の4点です。この4点が揃った報告書は、後日の設計変更・トラブル調査で「再利用できる資産」になります。
よくある誤読の実例
- 特異点のピーク値で疲労寿命を計算 — メッシュ依存の無意味な数値。ホットスポット法・公称応力法など評価法側の規約に従う
- 平均表示で応力ジャンプが隠れたまま「収束した」と判断 — 非平均表示での確認を検証手順に組み込む
- 1ケースの結果で「安全」と断定 — 入力(材料・荷重・境界剛性)のばらつきに対する感度を少数ケースでも確認(スクリーニング感度分析)
- 色の印象で優劣判断 — スケールが違う2図の比較。必ず同一スケール・差分表示で
ツール別の後処理機能
ツール別の要点
| ツール | 平均化の制御 | 実務メモ |
|---|---|---|
| Abaqus/Viewer | Averaging threshold(既定75%):隣接要素の値の差が閾値超なら平均しない | 閾値を100%にすると常時平均・0%で非平均。既定のままだと図ごとに平均範囲が変わる点に注意 |
| Ansys Mechanical | Averaged / Unaveraged / Nodal Difference表示 | Nodal Difference(平均と非平均の差)がメッシュ品質マップとして便利 |
| Nastran系ポスト | 要素応力/節点応力の出力選択(コーナー出力) | CQUAD4のセンター値かコーナー値かで分布の鋭さが変わる |
| ParaView(OpenFOAM等) | Cell Data / Point Data の明示的変換フィルタ | セルデータのまま等値面を切ると階段状に。変換の有無を意識する |
後処理の自動化——スクリプトで再現可能に
報告書の図を手作業のスクリーンショットで作ると、再解析のたびに視点・スケールが揺れて比較不能になります。ParaView(Pythonトレース)、Abaqus(Pythonスクリプト)、Ansys(Mechanicalスクリプト/ACT)はいずれも後処理の完全自動化に対応しており、「視点・スケール・注記を固定した標準図セット」をスクリプト化すると、ケース比較・改訂対応・監査対応のコストが桁で下がります。設計変更のたびに同じ図が同じ規約で再生成される状態が理想です。
先端動向——結果解釈の高度化
1ケース表示から差分・統計の可視化へ
設計比較では「AとBのコンターを並べる」より差分場 \( \Delta\sigma(\mathbf{x}) \) を直接描く方が情報量が多く、UQの文脈では多数ケースの平均場と標準偏差場(どこがばらつくか)の可視化が標準になりつつあります。ばらつきの大きい領域は入力不確かさに敏感な領域であり、設計ロバスト性の議論に直結します。
Web3D・軽量フォーマットによる結果共有
glTF等の軽量3Dフォーマットへ結果をエクスポートし、ブラウザで回転・プローブできる形で共有する運用が広がっています。静止画のコンター図と違い、受け手が自分で視点とスケールを操作できるため、「見せたい絵」ではなく「データそのもの」を渡せるのが本質的な進歩です。設計レビューの合意形成が速くなります。
AIによる異常検知・解釈支援
過去解析のデータベースから「通常と違う結果」を自動フラグする異常検知、結果場から危険部位を提案する学習モデルなど、解釈の初期スクリーニングをAIが担う研究が進んでいます。ただし現状の実務的な位置づけは見落とし防止の補助であり、本記事で述べた収支チェック・平均/非平均確認のような検証規律を置き換えるものではありません。
トラブル対応
症状別の原因と対策
| 症状 | 考えられる原因 | 対策 |
|---|---|---|
| メッシュを細かくすると最大応力が増え続ける | 応力特異点を評価している | 特異点の判別(収束テスト)。評価法をホットスポット/公称応力へ |
| 平均表示と非平均表示で値が大きく違う | メッシュ不足、異材界面をまたぐ平均 | 該当部を局所細分、平均範囲を材料単位に制限 |
| 同じモデルなのにツールで壁面熱伝達率が違う | 基準温度・壁処理・再構成法の定義差 | 定義を文書で確認し、比較は熱流束そのもので行う |
| コンター図が斑・ノイズ状 | セルデータの直接表示、品質不良要素 | ポイントデータ変換、要素品質の確認 |
| 報告後に「前回と比較できない」と指摘される | スケール・視点・平均設定が図ごとに異なる | 後処理スクリプト化で図の規約を固定 |
| 結果の数値が理論オーダーと合わない | 単位系の不整合、出力量の定義違い(全圧/静圧等) | 単位系を入力から一気通貫で確認、量の定義をマニュアルで確認 |
最後の砦は物理感覚
チェック項目が多くて、全部やる自信がないです…。一番大事なものを1つ選ぶなら?
1つ選ぶなら「手計算のオーダー見積りと比べる」だよ。梁の曲げ応力、円管の圧損、集中熱源の温度上昇——単純化した手計算と解析結果が2倍以内で合っていれば、大きな入力ミスはまず無い。逆に桁が違えば、どれだけ図が綺麗でもどこかが壊れている。解析ツールが高度になるほど、この素朴な検算の価値は上がり続けている。コンター図は「答え」じゃなくて「検算に合格した後の詳細情報」——この順序を守るのが、結果解釈のプロの流儀だね。
なった
詳しく
報告