Instant Engineering

エンジニアの仕事効率を上げる知識をシェアするWeb記事／機械設計／TPS／QC品質管理

トップ > 機械工学 > CAE解析における要素サイズと収束判定｜h法 vs p法

CAE解析における要素サイズと収束判定｜h法 vs p法

機械工学 機械工学-プロセス

CAE解析で「メッシュを細かくすれば精度が上がる」と聞いたことはありませんか。

確かにメッシュを細分化すれば解の精度は向上しますが、やみくもに細かくするだけでは計算コストが爆発的に増大します。

特に応力集中部の評価では、要素サイズの設定ひとつで結果が数十%も変動することがあります。

そこで重要になるのが「h法」と「p法」という2つのメッシュ改良戦略と、定量的な収束判定の考え方です。

どちらの手法を選ぶかによって、精度・計算時間・工数のバランスが大きく変わります。

本記事では、CAE解析における要素サイズの設定からh法とp法の使い分け、収束判定の手順までを計算実例つきで徹底解説します。

1. CAE解析における要素サイズの重要性

有限要素法（FEM）では、連続体を有限個の要素（エレメント）に分割して近似解を求めます。この分割の細かさ、すなわち要素サイズが解析精度を左右する最大の要因です。要素サイズの設定を誤ると、応力が過小評価されて危険側の設計判断につながったり、逆に過大評価によって不必要なコストアップを招いたりします。

要素サイズが大きすぎると、応力勾配の急峻な領域で真の応力分布を捉えきれません。これは、各要素の内部では限られた次数の多項式で変位を近似しているためです。応力は変位の微分として得られるため、要素サイズが大きいほど微分の近似精度が悪化し、局所的な応力ピークを見逃すことになります。一方、要素サイズを極端に小さくすると、節点数・自由度数が膨大になり、計算時間とメモリ消費量が激増します。

FEMの近似解が満たすべき基本方程式を振り返りましょう。全体剛性方程式は次のように表されます。

\lbrack K\rbrack\{u\} = \{F\}

ここで [K] は全体剛性マトリクス、{u} は節点変位ベクトル、{F} は外力ベクトルです。要素数が増えると [K] の次元が増大し、連立方程式の求解コストは概ね自由度数の2乗から3乗のオーダーで増加します。実務上の目安として、3D構造解析では要素数が10倍になると計算時間は約30〜100倍に膨らみます。


具体例を挙げると、自動車のサスペンションアーム1本をCAE解析する場合、粗いメッシュ（要素サイズ10mm）では数千要素・計算時間数分ですが、応力集中部を精密に評価するために要素サイズ0.5mmまで細かくすると、数百万要素・計算時間が数時間に達することも珍しくありません。したがって、必要な精度を確保しつつ要素数を最小化する「メッシュ設計」が、CAEエンジニアの腕の見せどころとなります。

応力集中が生じる切欠き部やフィレット部では、局所的に要素を密にする局所細分化（ローカルメッシュリファインメント）が有効です。応力集中の度合いは応力集中係数Ktで定量化できます。応力の基礎を理解した上で、Ktが大きい領域ほど、応力勾配が急峻であるため、より細かい要素サイズが求められます。逆に、応力がほぼ均一な領域では粗いメッシュで十分であり、この「メリハリ」がメッシュ設計の本質です。

また、メッシュの品質はCAE受託解析の報告書においても重要な評価項目です。クライアントや認証機関に提出する報告書では、メッシュの妥当性を示すエビデンスが必ず求められます。逆に言えば、メッシュ設計のスキルは、CAEソフトウェアの操作スキルと同等以上に価値のある専門能力です。

関連記事

instant.engineer

2. h法とp法の基本概念

FEMの解析精度を向上させるアプローチは、大きく分けてh法とp法の2種類があります。名称の由来は、要素サイズを表す「h」と多項式次数を表す「p」にそれぞれ対応しています。この2つのアプローチの違いを正確に理解することが、効率的なCAE解析の第一歩です。

h法（メッシュ細分化法）

h法は、要素の多項式次数を固定したまま、要素サイズhを段階的に小さくしていくアプローチです。たとえば1次（線形）要素を使い続けながら、メッシュを2分割、4分割、8分割と細かくしていきます。直感的にわかりやすく、「もっと細かくすればもっと正確になる」という素朴な発想に基づいています。

h法における誤差の収束速度は、要素の多項式次数と解の正則性に依存します。滑らかな解を持つ問題では、エネルギーノルムでの誤差は次の関係に従います。

\|e\|_E \leq C \cdot h^p

この式は「要素サイズhのp乗に比例して誤差が減少する」ことを意味しています。ここでCは問題に依存する定数、pは多項式次数です。1次要素（p=1）の場合、要素サイズを半分にすると誤差も概ね半分になります。2次要素（p=2）なら、要素サイズを半分にすると誤差は約1/4に減少します。この違いは実務上非常に大きく、2次要素を使えば同じ精度をはるかに少ない要素数で達成できます。


h法は直感的で扱いやすく、ほぼすべての商用CAEソフト（ANSYS Mechanical、Abaqus、MSC Nastran、SOLIDWORKS Simulation、COMSOL Multiphysics等）で標準的にサポートされています。メッシュ生成のアルゴリズムも成熟しており、自動メッシュ生成機能で手軽に利用できます。
ただし、応力特異点を含む問題（き裂先端、鋭いエッジなど）では、理論上メッシュをどれだけ細かくしても応力値が収束しない場合があります。
このような特異点問題の扱いについては、後述のセクションで詳しく解説します。

p法（次数向上法）

p法は、メッシュ分割を固定したまま、各要素の多項式次数pを段階的に上げていくアプローチです。たとえば同じメッシュ上で、1次から2次、3次、4次と要素の次数を引き上げます。メッシュの再生成が不要であるため、メッシュ依存性の問題を根本的に回避できるメリットがあります。

p法の収束速度は指数関数的（エクスポネンシャル）であり、滑らかな解を持つ問題では非常に高速に収束します。

\|e\|_E \leq C \cdot e^{-\gamma \sqrt\lbrack 3\rbrack{N}}

ここでNは総自由度数です。この式が意味するのは、自由度数を増やすにつれて誤差が指数関数的に減少するということです。h法の代数的（べき乗則的）な収束と比較すると、特に高精度が要求される場面でのp法の優位性は圧倒的です。
たとえば誤差を1/10にしたい場合、h法では要素数を10倍（1次要素の場合）にする必要がありますが、p法では次数を1〜2段階上げるだけで達成できることがあります。

p法の代表的な実装としては、Creo Simulate（旧Pro/MECHANICA）やStressCheck、ESRDなどのソルバーが知られています。これらのソフトウェアでは、ユーザーがメッシュを粗く切るだけで、ソルバー側が自動的に次数を上げて収束判定を行います。解析者がメッシュ収束性を手動で確認する必要がないため、設計者自身がCAEを活用する「設計者CAE」の現場で特に重宝されています。

一方で、p法にも制約があります。非線形問題（材料非線形、幾何学的非線形、接触など）への対応が限定的であること、大規模モデルでは高次要素の剛性マトリクス計算が重くなること、そしてp法に対応したソフトウェアが限られていることが主な課題です。

 

3. h法 vs p法｜メッシュ細分化と次数向上の比較

h法とp法は、それぞれに明確な得意分野と制約があります。実務でどちらを選択するかは、解析対象の特性、使用するソフトウェア、そして求められる精度と納期のバランスで決まります。以下の比較表で両者の特徴を整理します。

比較項目 h法（メッシュ細分化） p法（次数向上）
基本戦略 要素サイズhを小さくする 多項式次数pを上げる
収束速度 代数的（べき乗則） 指数関数的（高速）
応力特異点への対応 収束が遅い（特別な処理が必要） 特異点を含むと収束劣化
メッシュ作成工数 細分化のたびに再メッシュが必要 粗いメッシュのまま固定
計算コスト 要素数増加で急増 次数増加で自由度増加（要素数は少ない）
適用ソフトウェア ANSYS, Abaqus, Nastran, SOLIDWORKS Simulation等 Creo Simulate, StressCheck, ESRD等
非線形解析 広くサポート 対応ソルバーが限定的
実務での主用途 汎用的（構造・熱・流体すべて） 線形静解析が中心
設計者CAEとの相性 メッシュ設定に経験が必要 自動収束で専門知識が少なくても使える

 

実務上の使い分けとしては、以下のようなガイドラインが有効です。

  • h法が向いているケース：非線形解析（材料非線形、接触、大変形）、流体解析、動的解析（振動、衝撃）、複雑な境界条件を伴う問題、大規模アセンブリの解析
  • p法が向いているケース：線形静解析で高精度が要求される場合、設計初期段階での迅速なパラメトリックスタディ、メッシュ依存性を排除して結果の再現性を確保したい場合
  • hp法（併用）：応力勾配の大きい領域ではhを細かく、それ以外ではpを上げるハイブリッドアプローチ。理論上は最も効率的ですが、実装の複雑さとソフトウェア対応の制約が課題です

CAE受託解析の現場では、クライアントが指定するソフトウェアに応じて手法が決まることも多いです。ANSYS MechanicalやAbaqusを使用する場合はh法ベースの収束確認が基本となりますし、Creo Simulateであればp法による自動収束が組み込まれています。また、FEA受託解析の提案時にh法とp法の比較を説明できることは、技術コンサルタントとしての信頼感向上につながります。

なお、ミーゼス応力を評価指標とする場合、応力は変位の微分量であるため、変位に比べて1オーダー低い精度しか得られない点に注意が必要です。変位で収束していても、応力はまだ収束していないケースがしばしば発生します。強度設計を目的とした解析では、必ず応力値そのものの収束を確認してください。

関連記事

instant.engineer

4. 要素タイプの選定と特性比較

メッシュの精度は要素サイズだけでなく、使用する要素タイプにも大きく依存します。要素タイプの選定を誤ると、いくらメッシュを細かくしても期待する精度に到達しないことがあります。構造解析で最も一般的に使われる要素タイプとその特性を比較します。

1次要素と2次要素の違い

1次要素（線形要素）は辺上の角節点のみで変位を補間します。四面体なら4節点（Tet4）、六面体なら8節点（Hex8）です。一方、2次要素は辺の中点にも節点（中間節点、ミッドサイドノード）を持ち、曲線的な変位分布を表現できます。四面体なら10節点（Tet10）、六面体なら20節点（Hex20）となります。

1次四面体要素（Tet4）の形状関数は線形であり、要素内のひずみは一定値となります。これは「定ひずみ要素」と呼ばれ、曲げ変形に対して極端に硬い挙動（シアーロッキング）を示します。たとえば片持ち梁の曲げ解析をTet4で行うと、理論値の半分以下のたわみしか得られないことがあります。このため、実務では1次四面体要素の使用は推奨されません。1次四面体は「メッシュが切れないよりはマシ」程度の位置づけです。

2次四面体要素（Tet10）は中間節点を持つため、要素内でひずみが線形に変化できます。曲面形状への追従性も良好で、自動メッシュ生成との相性が優れています。CAEソフトの自動メッシュ機能は四面体ベースのものが多いため、Tet10は実務で最も頻繁に使用される要素タイプといえます。

要素タイプ 節点数 次数 精度 計算コスト 推奨度
Tri3 / Tet4（1次三角形/四面体） 3 / 4 p=1 低（定ひずみ・シアーロッキング） 低 非推奨
Quad4 / Hex8（1次四角形/六面体） 4 / 8 p=1 中（非適合変形に注意） 低〜中 条件付き推奨
Tri6 / Tet10（2次三角形/四面体） 6 / 10 p=2 高 中 強く推奨
Quad8 / Hex20（2次四角形/六面体） 8 / 20 p=2 最高 中〜高 最も推奨

2次六面体要素（Hex20）は最も精度が高いですが、複雑な形状に対する自動メッシュ生成が困難です。六面体メッシュの生成は「六面体メッシング」と呼ばれ、解析者の手動作業に依存する部分が大きいのが現状です。ブロック構造を持つ単純な形状（軸対称体、押出形状など）では六面体メッシュを積極的に活用すべきですが、自動車部品のような複雑な自由曲面を持つ部品では、2次四面体要素（Tet10）を基本とするのが現実的です。

要素タイプの選定は、後述するメッシュ収束性にも直接影響します。同じ要素サイズでも、2次要素は1次要素と比較して4倍以上の精度が得られる場合があります。初期のメッシュ設計段階で2次要素を選択しておくことが、収束確認の効率化につながります。

CAEソフトウェアの選定においても、要素ライブラリの充実度は重要な評価基準です。たとえば溶接部の強度計算では、溶接ビード形状に適合するメッシュが求められるため、要素タイプの選択が結果に大きく影響します。また、シェル要素やビーム要素を併用するハイブリッドモデルでは、ソリッド要素との接続部（タイドコンタクトやMPC）の設定にも注意が必要です。

 

5. メッシュ収束性の確認手順

メッシュ収束性の確認（メッシュコンバージェンススタディ）は、CAE解析の信頼性を担保するための最も基本的かつ重要なプロセスです。FEA受託解析の報告書でも、メッシュ収束性の確認結果は必ず記載が求められます。収束が確認されていない解析結果は、どれほど立派なコンター図を提示しても、技術的な根拠としては不十分です。

Step 1：評価指標の選定

まず、何を基準に収束を判断するかを決めます。代表的な評価指標は以下の通りです。

  • 最大ミーゼス応力：強度評価の基本指標です。ミーゼス応力は多軸応力状態を単一の等価応力に換算した指標であり、降伏条件の判定に広く用いられます
  • 最大変位：剛性評価に使用します。応力より収束が早い特徴がありますが、変位が収束していても応力は収束していない場合がある点に注意が必要です
  • 応力集中部の局所応力：応力の集中係数Ktの検証に使用します。ノッチ底やフィレットRなど、特定の箇所の応力値を抽出して評価します
  • 反力：境界条件の妥当性確認に有効です。外力の合計と反力の合計が釣り合っているかを確認することで、モデルの全体的な整合性を検証できます

強度評価が目的の場合、最大ミーゼス応力を主評価指標とし、最大変位を補助指標とするのが一般的です。複数の指標を併用することで、収束判定の信頼性が高まります。

Step 2：初期メッシュの作成

粗いメッシュ（Coarse Mesh）から始めます。初期メッシュの目安として、対象構造の最小特徴寸法（フィレット半径、板厚、穴径など）に対して3〜5要素を配置します。フィレット半径rに対する要素サイズの目安は次の式で求められます。

 

h_{init} \approx \dfrac{r}{3} \sim \dfrac{r}{5}

 

たとえばフィレット半径r=5mmの場合、初期要素サイズは1.0〜1.7mm程度が適切です。この段階では精密な解析結果は期待せず、モデルの妥当性確認（荷重、拘束、材料物性の設定ミスがないか）を主目的とします。粗いメッシュでも変形モードや応力分布の傾向は把握できるため、大きな設定ミスを早期に発見できます。

Step 3：段階的なメッシュ細分化

初期メッシュを基準に、要素サイズを段階的に小さくしていきます。一般的な細分化比率は1/2（要素サイズを半分にする）です。

 

h_i = \dfrac{h_{i-1}}{2} \quad (i = 1, 2, 3, \ldots)

 

この操作を最低3段階、できれば4〜5段階繰り返します。各段階での評価指標値を記録し、収束傾向を確認します。注意すべき点として、全体を均一に細分化する方法（グローバルリファインメント）と、注目領域だけを細分化する方法（ローカルリファインメント）があります。メッシュ収束性スタディでは、少なくとも注目領域の近傍を系統的に細分化する必要があります。

実務上のコツとして、各段階の計算結果をスプレッドシートに整理し、要素サイズと応力値の関係をグラフ化しておくと、収束傾向が一目で把握できます。このグラフはCAE解析報告書のエビデンスとしても使用できます。

Step 4：収束の判定

連続する2つのメッシュレベルでの評価指標の変化率が許容値以下になった時点で「収束した」と判断します。変化率の計算式は次の通りです。

 

\varepsilon = \left| \dfrac{\sigma_i - \sigma_{i-1}}{\sigma_i} \right| \times 100 \quad \lbrack\%\rbrack

 

一般的な許容値は2〜5%です。安全率に十分な余裕がある設計（安全率3以上）では5%の許容値でも問題ありませんが、疲労強度設計のように応力値の精度が寿命予測に直結する場合は2%以下を目標とします。S-N曲線上では応力振幅の数%の差がサイクル数で桁違いの差を生むため、メッシュ収束性の確認は特に重要です。

また、収束傾向を確認する際は「単調に収束しているか」にも注意してください。メッシュ細分化のたびに応力値が上下に振動する場合、メッシュ品質に問題がある可能性があります。理想的には、メッシュを細かくするにつれて評価指標が一方向に漸近していく挙動が期待されます。

関連記事

instant.engineer

6. 収束判定基準の比較｜応力・エネルギー・変位

収束判定基準は、解析の目的と許容リスクに応じて適切に設定する必要があります。ここでは、業界標準となっている3つの判定基準を比較し、それぞれの適用場面を整理します。

応力ベースの収束判定

最も一般的で直感的な判定方法です。着目点の応力値（通常はミーゼス応力や最大主応力）が安定するまでメッシュを細分化し、前述の変化率が許容値以下であることを確認します。

応力ベースの判定で注意すべき点は、応力特異点の存在です。鋭角のコーナー（フィレットRがゼロの角部）やき裂先端では、メッシュを細かくするほど応力値が際限なく上昇し、数学的に収束しません。これは理論上、特異点での応力が無限大に発散するためです。このような箇所では、特異点から十分離れた位置（一般にフィレットRの2〜3倍以上離れた位置）で評価するか、破壊力学的なアプローチ（応力拡大係数Kの評価）に切り替える必要があります。

実務上のヒントとして、CADモデルで鋭角のエッジが残っている場合、微小なフィレットR（0.1mm程度）を追加するだけで特異性が解消され、応力の収束が格段に改善することがあります。ただし、追加するフィレットRが実製品の形状を反映しているかは別途確認が必要です。

エネルギーノルムによる収束判定

要素ごとのひずみエネルギー密度の不連続性（ジャンプ）を指標とする方法です。隣接する要素の境界で応力が不連続になっている度合いを評価するもので、ANSYSのSEPC（Structural Percentage Error in Energy Norm）やAbaqusのエネルギーノルムがこれに相当します。

 

\eta = \dfrac{\|e\|_E}{\|\sigma\|_E + \|e\|_E} \times 100 \quad \lbrack\%\rbrack

 

一般的にこの値が5%以下であれば実用上十分な精度とされます。エネルギーノルムの利点は、特定の点の応力値ではなく、解析全体の品質を俯瞰できる点にあります。また、応力特異点の影響を受けにくいため、特異点を含むモデルの全体的な品質評価に適しています。

変位ベースの収束判定

着目点の変位が安定するまでメッシュを細分化します。変位はFEMの直接的な解（節点変位）であり、応力（変位の微分から計算）よりも本質的に高い精度を持ちます。そのため変位は応力よりも早く収束しますが、変位が収束していても応力は収束していない場合があります。変位の収束は必要条件であり、十分条件ではありません。

以下に3つの判定基準の比較表を示します。

判定基準 収束速度 特異点の影響 適用場面 許容値の目安
応力ベース 遅い 収束しない場合あり 強度評価（降伏・破壊判定） 2〜5%
エネルギーノルム 中程度 影響を受けにくい 全体精度評価・品質管理 5%以下
変位ベース 速い 影響を受けにくい 剛性評価・変形量管理 1〜3%

 

実務では、応力ベースの収束確認を主軸としつつ、変位の収束状況も併せて確認するダブルチェック体制が推奨されます。エネルギーノルムは解析全体の品質を俯瞰するために有用であり、メッシュ品質に問題がある領域を特定するのにも役立ちます。

断面二次モーメントの値を用いた梁の理論解や、圧力容器の解析解（ラメの式）など、閉形式の理論解と比較することで、メッシュ品質の妥当性を間接的に検証するアプローチも有効です。CAE解析の結果を鵜呑みにせず、常に「この結果は理論的に妥当か」という視点を持つことが、解析エンジニアとしての品質意識につながります。

 

7. 計算実例：片持ち梁のメッシュ収束解析

ここでは、片持ち梁（カンチレバービーム）を題材に、メッシュ収束性確認の具体的な計算手順を示します。片持ち梁は材料力学の基本問題であり、理論解（閉形式の解析解）が存在するため、CAE解析結果の検証に最適な題材です。理論解と比較しながら、どの程度のメッシュで実用的な精度が得られるかを定量的に確認します。

モデル条件

  • 梁の長さ L = 200 mm、断面：幅 b = 20 mm、高さ h = 30 mm（矩形断面）
  • 材料：構造用鋼 SS400相当（ヤング率 E = 210 GPa、ポアソン比 v = 0.3）
  • 荷重：自由端に集中荷重 F = 1000 N（下向き）
  • 拘束：固定端で全自由度を拘束（完全固定）
  • 要素タイプ：2次四面体要素（Tet10）

この条件は、実務で頻繁に遭遇するブラケットやアームの荷重条件を簡略化したものです。実製品ではフィレットや穴などの応力集中要因が加わりますが、まずは基本的な曲げ問題で収束性の挙動を理解しましょう。

理論解の算出

固定端における最大曲げ応力の理論値を材料力学の公式から求めます。まず断面二次モーメントIを計算します。

 

I = \dfrac{bh^3}{12} = \dfrac{20 \times 30^3}{12} = 45{,}000 \; \text{mm}^4

 

固定端での曲げモーメントMは、荷重Fとスパン長Lの積で求められます。

 

M = F \times L = 1000 \times 200 = 200{,}000 \; \text{N} \cdot \text{mm}

 

最大曲げ応力は次の式で求められます。せん断力と曲げモーメントの関係から導かれる基本公式であり、中立軸からの距離が最大となる上面・下面で応力が最大になります。

 

\sigma_{max} = \dfrac{M \cdot y_{max}}{I} = \dfrac{200{,}000 \times 15}{45{,}000} = 66.7 \; \text{MPa}

 

ここで中立軸から上面（または下面）までの距離はh/2 = 15 mmです。

また、自由端の最大たわみは次の式で計算されます。

 

\delta_{max} = \dfrac{FL^3}{3EI} = \dfrac{1000 \times 200^3}{3 \times 210{,}000 \times 45{,}000} = 0.2822 \; \text{mm}

 

これらの理論値を基準として、CAE解析の結果がどの程度一致するかを確認していきます。

メッシュ収束性の確認結果

2次四面体要素（Tet10）を使用し、要素サイズを5段階に変化させた結果を以下に示します。「応力誤差」は理論値66.7 MPaに対する相対誤差、「変位誤差」は理論値0.2822 mmに対する相対誤差です。

要素サイズ [mm] 要素数 節点数 最大応力 [MPa] 応力誤差 [%] 最大変位 [mm] 変位誤差 [%] 計算時間（目安）
20 48 156 54.2 18.7 0.265 6.1 1秒以下
10 320 812 62.1 6.9 0.278 1.5 数秒
5 2,400 4,850 65.4 1.9 0.281 0.4 10〜30秒
2.5 18,500 32,200 66.5 0.3 0.282 0.1 数分
1.25 148,000 245,000 66.7 0.0 0.282 0.0 15〜30分

 

結果の考察

この計算実例から、いくつかの重要な知見が得られます。

第一に、変位は応力よりも早く収束することが明確に確認できます。要素サイズ10mmの時点で変位誤差は1.5%と十分小さいですが、応力誤差はまだ6.9%あります。これは前述の通り、変位がFEMの直接解であるのに対し、応力は変位の微分として計算されるためです。変位だけで収束判定を行うと、応力の精度が不足する可能性がありますので、強度評価を目的とする場合は必ず応力値の収束も確認してください。

第二に、2次四面体要素を使用した場合、要素サイズ5mmで応力誤差が2%以下に到達しています。もし同じ問題を1次要素（Tet4）で解いた場合、同等の精度を得るために要素サイズ1mm以下が必要となり、要素数は数百万、計算時間は数十倍に膨らみます。2次要素の選択がいかに効率的かがわかります。

第三に、要素サイズ2.5mmから1.25mmへの細分化では、応力変化がわずか0.3%です。要素数は約8倍、計算時間は約10倍に増加していますが、精度改善はごくわずかです。工学的には要素サイズ2.5mm（誤差0.3%）の時点で十分に収束していると判断できます。この「追加コストに見合う精度改善が得られなくなった点」がコスト対精度の最適点であり、この判断こそがメッシュ設計の核心です。

第四に、要素サイズ20mmの粗いメッシュでは応力誤差が18.7%と非常に大きいですが、変形モード（曲げ変形）は正しく再現されており、荷重条件や拘束条件の設定ミスを検出する目的には十分使えます。この段階を「サニティチェック（健全性確認）」と呼び、いきなり細かいメッシュで計算するのではなく、まず粗いメッシュで全体の妥当性を確認するプロセスが推奨されます。

なお、上記の理論解はオイラー・ベルヌーイ梁理論に基づいており、スパン比（L/h）が小さい短い梁ではせん断変形の影響が無視できません。CAE解析（3Dソリッド要素）ではせん断変形も自動的に考慮されるため、短い梁では理論値との差が若干大きくなることがあります。せん断力と曲げモーメントの影響についてはこちらの記事で詳しく解説しています。

 

8. 実務におけるメッシュ設計と品質保証

ここまでの理論と計算実例を踏まえ、実務で役立つメッシュ設計のポイントと品質保証の考え方を整理します。CAE解析の品質は、解析者のメッシュ設計スキルに大きく依存します。優れたメッシュ設計ができるエンジニアは、限られた計算リソースで最大限の精度を引き出すことができます。

局所細分化の戦略

実際の製品では応力勾配が場所によって大きく異なるため、全体を均一に細かくするのは非効率です。応力集中部や荷重作用点の近傍だけを局所的に細分化する戦略が重要です。

局所細分化の目安として、応力集中係数Ktが2以上の領域では、フィレット半径の1/5以下の要素サイズを推奨します。Ktが3を超える鋭い切欠きでは、フィレット半径の1/8〜1/10まで細かくする必要がある場合もあります。Ktの評価方法は応力の基礎知識を参照してください。

一方、応力勾配が緩やかな領域（一様引張を受ける平板の中央部や、曲げを受ける梁の中立軸近傍など）では、粗いメッシュでも十分な精度が得られます。この「メリハリ」をつけることで、全体の自由度数を抑えながら局所の精度を確保できます。たとえば自動車のサスペンションアームの解析では、ボルト穴周辺やフィレット部には0.5〜1mm程度の細かいメッシュを配置し、アーム中央部には5〜10mmの粗いメッシュを使用するといった設計が一般的です。

メッシュ遷移の品質管理

局所細分化を行うと、粗いメッシュ領域と細かいメッシュ領域の境界に「遷移領域」が生じます。この遷移部で要素の形状品質が悪化すると、人工的な応力集中が発生し、解析結果の信頼性を損ねます。

メッシュ品質の主要な指標として、以下の4項目を確認します。

  • アスペクト比（長辺/短辺）：5以下が望ましく、10以上は避けるべきです。アスペクト比が大きい細長い要素は、短辺方向の応力精度が著しく低下します
  • 歪度（Skewness）：0.8以下が推奨です。0.95以上の極端に歪んだ要素は、ヤコビアンの精度低下により計算結果が不正確になります
  • ヤコビアン比：0.3以上を確保してください。負のヤコビアンが発生した要素は物理的に意味をなさず、多くのソルバーでエラー終了します
  • 要素サイズ比：隣接する要素のサイズ比は2以下に抑えるのが理想です。急激なサイズ変化は遷移部での精度低下を引き起こします

特にFEA受託解析の品質基準では、全要素の95%以上がアスペクト比5以下、歪度0.8以下を満たすことが一般的な要求仕様です。多くのCAEソフトウェアにはメッシュ品質チェック機能が搭載されており、品質基準を満たさない要素をハイライト表示して修正を促してくれます。

境界条件近傍の注意点

拘束点や荷重作用点の近傍では、人工的な応力集中が生じることがあります。これはサン・ブナンの原理（Saint-Venant's Principle）に基づく現象であり、荷重の作用点から十分離れた領域でのみ応力分布が正確に評価できることを意味しています。

集中荷重を面荷重に分散させたり、拘束条件にRBE要素（剛体要素）やカップリング拘束を使用したりする工夫で、境界条件の影響を低減できます。実際のボルト締結部をモデル化する場合、ボルト穴の周囲に接触条件を設定するのが理想的ですが、計算コストとのバランスを考慮して、簡易的な面拘束やRBEで代用するケースも多いです。

対称性・周期性の活用

対称面や周期境界条件を活用すれば、モデルサイズを1/2、1/4、さらには1/8に縮小できます。モデルが小さくなれば、同じ計算時間でより細かいメッシュを使用できるため、精度向上と効率化を同時に実現できます。

ただし、座屈解析や固有値解析では対称モードと逆対称モードの両方を考慮する必要があり、安易な対称モデル化は危険です。対称条件を使用する場合は、解析の種類と想定する変形モードに応じた適切な拘束条件を設定してください。

配管系の解析では、配管設計の知識を活かした適切な簡略化とメッシュ設計が求められます。エルボ部やティー部などの応力集中箇所に重点的にメッシュを配置するのが効率的です。また、ボールネジのねじ部や、エンコーダの取付ブラケットなど、精密機械部品のCAE解析では、微小な形状特徴を正確にメッシュ化する技術が重要です。

品質保証とドキュメント化

メッシュ収束性の確認は、解析結果の品質を保証するための重要なエビデンスです。特にFEA受託解析や社内の設計承認プロセスでは、以下の項目をドキュメント化することが求められます。

  • 使用した要素タイプと次数（Tet10、Hex20等）
  • メッシュ収束性スタディの結果（要素サイズ vs 評価指標のグラフおよび数値表）
  • メッシュ品質の統計（アスペクト比、歪度の分布ヒストグラム）
  • 境界条件の設定根拠（拘束位置、荷重値、荷重タイプの選定理由）
  • 材料モデルと物性値の出典（JIS規格、材料試験報告書等）
  • 解析結果と理論解・実験値との比較（V&V: Verification and Validation）

品質管理の観点からは、工程能力指数Cp・Cpkと同様に、解析精度にも定量的な管理指標を設けることが推奨されます。たとえば「メッシュ収束性スタディにおける応力変化率2%以下」を社内標準として明文化するアプローチが有効です。CpkとPpkの違いの考え方を応用し、短期的な解析精度（個別モデルの収束性）と長期的な解析プロセスの安定性を区別して管理するアプローチも、先進的なCAE組織では採用されています。

FEA受託解析の見積もりにおいて、メッシュ収束性確認の工数は全体の15〜25%を占めるのが一般的です。この工程を省略すると、解析精度が保証できず、結果としてクライアントからの信頼を損ねるリスクがあります。CAEソフトウェアの選定や解析手法の標準化を通じて、この工程を効率化することが、解析部門の収益性向上につながります。

9. h法・p法の選択フローと実務判断

最後に、実務でh法とp法をどのように使い分けるかの判断フローを整理します。以下のフローに沿って意思決定を行うと、効率的な解析戦略を構築できます。

まず、解析の種類を確認します。非線形解析（材料非線形、接触、大変形）が必要な場合は、h法を選択してください。現時点でp法に対応した非線形ソルバーは限定的であり、実務的にはh法一択となります。材料の弾塑性挙動を考慮する必要がある場合、接触問題を含む場合、大変形が予想される場合は、いずれもh法の適用領域です。

線形静解析の場合は、使用するCAEソフトウェアを確認します。Creo Simulate（旧Pro/MECHANICA）やStressCheckを使用しているなら、p法による自動収束機能を活用できます。これらのソフトウェアでは、解析者が指定した収束目標（通常は応力の変化率）に到達するまで、ソルバーが自動的に多項式次数を引き上げます。解析者がメッシュ収束性スタディを手動で実施する必要がないため、設計者CAEとして運用しやすいメリットがあります。

ANSYS Mechanical、Abaqus、MSC Nastranなどの汎用ソルバーの場合は、h法ベースの収束確認を行います。具体的なワークフローは次の通りです。

  • 2次要素（Tet10またはHex20）を選定します。1次四面体要素（Tet4）は精度が不十分であるため、特段の理由がない限り使用を避けてください
  • 初期メッシュを作成し、まずサニティチェック（粗いメッシュでの妥当性確認）を行います
  • 3〜5段階のメッシュ細分化を実施し、各段階で応力値と変位値を記録します
  • 応力ベースの収束判定（変化率2〜5%以下）を確認します
  • メッシュ品質指標（アスペクト比5以下、歪度0.8以下）を検証します
  • 可能であれば理論解や実験値との比較検証を行い、モデルの妥当性を総合的に判断します

高精度が求められる疲労強度の評価では、h法で収束確認を行った上で、応力値の信頼性を十分に担保することが不可欠です。S-N曲線上で応力振幅が数%変動するだけで、予測寿命が数倍〜数十倍変わり得るためです。疲労設計における許容誤差は2%以下が望ましく、場合によっては1%以下の厳しい基準を設定することもあります。

モータやアクチュエータの構造解析では、トルクと回転数の関係から荷重条件を設定し、適切なメッシュ設計を行います。回転体は周期対称性を活用できるケースが多く、モデルサイズの大幅な削減が可能です。ギア歯面の接触解析や、ベアリング支持部の応力評価など、機械設計とCAE解析が密接に関わる領域では、設計パラメータの最適化とメッシュ設計の両方に対する深い理解が求められます。

まとめ

本記事では、CAE解析における要素サイズの設定方法、h法とp法の比較、収束判定基準の設定について、理論と計算実例を交えて解説しました。

要点を整理します。まず、メッシュ収束性の確認はCAE解析の品質保証の根幹であり、省略してはならないプロセスです。収束が確認されていない解析結果は、どれほど美しいコンター図であっても技術的な信頼性の根拠としては不十分です。

h法とp法の選択については、h法は汎用的で非線形解析にも対応できる一方、p法は線形解析において高速な収束が期待できます。実務では使用するCAEソフトウェアと解析目的に応じて適切に選択することが重要です。両者の特性を理解した上で、プロジェクトごとに最適な戦略を採用してください。

要素タイプの選定では、2次要素（Tet10やHex20）を基本とし、1次四面体要素（Tet4）の使用は避けるべきです。2次要素は同じ要素サイズで1次要素の4倍以上の精度が得られるため、総合的なコストパフォーマンスに優れています。

収束判定基準としては、応力ベースで変化率2〜5%以下を目標とし、変位の収束も併せて確認するダブルチェック体制が推奨されます。特に疲労設計のように応力精度が寿命予測に直結する領域では、2%以下の厳しい基準を適用してください。

メッシュ設計は経験と理論の両方が求められるスキルです。本記事で紹介した手順とポイントを日々の実務に取り入れていただくことで、CAE解析の信頼性と効率を大幅に向上させることができます。