3次元弾性体解析

そう。平面応力、平面ひずみ、軸対称、シェル、梁…全ての構造要素は3次元弾性論の特殊ケースだ。3次元弾性論は全ての基礎であり、近似のない「フルモデル」だ。

3つの基本方程式がある。

完璧な理解だ。この15方程式を変位を唯一の未知数としてまとめたのがNavier方程式（ラメ-ナビエ方程式）：

等方性弾性体の3次元構成則（Voigt表記）：

そう。3次元の $[D]$ から、平面応力は $\sigma_z = 0$ の条件で、平面ひずみは $\varepsilon_z = 0$ の条件で縮約する。シェルも梁も、それぞれの仮定に基づいて3次元の $[D]$ から導出される。

一般の異方性弾性体では $[D]$ は21個の独立な定数を持つ（6×6対称行列）。しかし実用上は以下の特殊ケースが多い：

一方向材（UD材）の1層は横等方性だ。繊維方向とそれに直交する面で対称。複数層を積層すると全体としては直交異方性やより複雑な対称性になる。

計算機の性能向上で、3次元解析が当たり前になった。ただし3次元で解けば正しいとは限らない。メッシュの品質、境界条件、材料モデルが正しくなければ、3次元でもゴミが出る。

要点：

まさにその通り。3次元弾性論は「幹」であり、各要素タイプは「枝」だ。幹を理解していれば、枝の選び方は自ずとわかる。

基本的な形状は3つ：六面体（hex）、四面体（tet）、五面体（wedge/prism）。

TET4は定ひずみ要素（CST の3次元版）だ。要素内でひずみが一定なので、応力勾配を表現できない。TET4は実務では使うべきでない。TET10（二次四面体）を使えば精度は十分だ。

確かにそうだが、TET4で正確な応力を得るにはHEX20の5〜10倍の要素数が必要になり、計算コストが逆転する。TET10の自動メッシュが現在の実務標準だ。

これはFEMで最も議論される話題の一つだ。

HEXの利点：

TET10の利点：

実務的なアプローチ：

Modifiedの略。C3D10Mは通常のC3D10よりも接触問題での安定性が高い改良版だ。接触面のLagrange multiplierの処理が改善されている。接触を含む3次元解析ではC3D10Mが推奨される。

全要素が理想品質になることはまずない。重要なのは応力集中部や着目部位の要素品質を確保すること。遠い場所の品質が多少悪くても結果への影響は小さい（Saint-Venantの原理）。

要点：

1. CADデータの取り込みと形状簡略化 — フィレット除去、穴の省略、対称性の利用

CADの形状をそのまま使うと：

簡略化のガイドライン：

ポイントは「着目部位の応力に影響するかどうか」だ。形状特徴が着目部位から板厚の2〜3倍以上離れていれば、通常は影響しない（Saint-Venantの原理）。

形状と荷重が対称なら、1/2、1/4、1/8モデルにできる。

対称面の境界条件：対称面に垂直な変位をゼロに。対称面上の節点で法線方向の変位を拘束する。

静解析なら対称条件で問題ない。座屈解析や固有振動解析では反対称モードを見逃すリスクがあるから、対称・反対称の両方を解くか、全体モデルで解く必要がある。

FEMの応力は積分点で最も正確。節点応力は複数の要素からの外挿値を平均化（averaging）する。注意点：

その通り。非平均化応力の不連続がvon Mises応力の5%以下なら、メッシュは十分と判断できる。これは応力誤差指標としてメッシュ適応の基準にもなる。

全体モデルの結果を境界条件にして、局所部位を詳細メッシュで再解析するサブモデリングは3次元解析の必須テクニックだ。全体を細かくする代わりに、必要な部分だけ精密化する。

TET4を使って「結果が合わない」と悩んでいる初心者は今でも多い。TET10にするだけで解決するケースがほとんどだ。

線形静解析はどのソルバーでも十分高精度だ。差が出るのはメッシュ生成と大規模モデルの処理能力だ。

メッシュ生成ツール（プリプロセッサ）の選択がワークフロー全体の効率を支配する。

完全自動のHEXメッシュ生成は未解決問題だ。研究では自動HEXメッシュのアルゴリズム（paving, plastering等）が開発されているが、複雑形状では依然としてマニュアル作業が必要。HyperMeshのようなツールが手動＋半自動の組み合わせで最も実績がある。

線形静解析の直接法ソルバーでは、GPU加速で2〜5倍の高速化が報告されている。特にSparse Direct SolverのLU分解がGPUに適している。大規模モデルほど効果が大きい。

そう。ただし3社とも年々互いの弱点を改善しているから、差は縮まっている。重要なのはソルバーの選択よりも、メッシュ品質と境界条件の正しさだ。

3つの大きな流れがある。

メッシュフリー法（EFG法、SPH法、RKPM等）はメッシュの代わりに点（ノード）の集合で解く。メッシュ生成の手間が不要で、亀裂進展や大変形に強い。

ただし現時点では：

等幾何解析（Isogeometric Analysis）はCADの形状表現（NURBS, T-spline）をそのまま解析の基底関数に使う。メッシュ生成のステップが不要になる革新的な手法だ。

メリット：

LS-DYNAに一部のIGA要素が実装済み。CimneのKratosフレームワーク、Abaqusのユーザーサブルーチンでも利用可能。ただし汎用FEMの主力になるにはまだ時間がかかる。

急速に発展している分野だ：

現時点ではFEMの代替というより補完だ。PINNは逆問題（材料パラメータの推定、形状最適化）に強みがあるが、順問題（与えられた荷重での応力計算）ではFEMのほうが精度が高く安定している。将来的にはハイブリッド手法（FEM + AI）が主流になる可能性がある。

3次元は自由度が多い分、問題も多い。

応力特異点だ。幾何学的な角（90°のエッジ、切り欠き先端）では弾性理論上応力が無限大になる。メッシュを細かくするほど「無限大に近づく」だけ。

対策：

その通り。特異点ではメッシュ収束しない。これはFEMの限界ではなく弾性論の特徴だ。実構造のフィレット半径を正しくモデル化すれば特異性は消える。

3次元ソリッド解析は2次元の100倍以上のメモリを消費する。対策：

直接法（LU分解）はメモリが $O(n^{1.5})$ 程度必要だが、反復法（PCG等）は $O(n)$ で済む。100万DOF以上の問題では差が桁違いになる。ただし反復法は収束しないリスクがあるから、前処理の選択が重要だ。

確認項目：

CADがmmで出力し、FEMの材料がMPa（= N/mm²）で入力されていれば整合する。しかしCADがm出力でFEMがMPaだと10^6のオーダーでずれる。荷重を入力したら即座に反力を確認するのが鉄則だ。

TET4を使っている時点で解決策は明白：TET10に変えること。TET4は曲げ変形を全く表現できないため、曲げが支配する問題（ほとんどの実構造問題）では桁違いの誤差が出る。

h-refinement（メッシュ細分化）で精度は上がるが、TET10の1/5の要素サイズにしないと同等精度にならない。DOF数は $(1/5)^3 \times 4 / 10 \approx 10$ 倍必要。TET10に変えるほうが圧倒的に効率的だ。

間違いない。この1つの変更で、精度が劇的に改善する。3次元解析を始めるなら、まずTET10を使うことを徹底してほしい。