プリズム要素 — CAE用語解説

その名の通り、底面が三角形で、それが高さ方向に押し出されてできる六面体の要素です。正確には、2つの三角形の面と3つの四角形の面で構成されます。例えば、底面の三角形の3節点と、それらを高さ方向に移動させた3節点、合計6節点で定義されるのが一次プリズム要素です。二次要素だと、辺の中間に節点が追加されて合計15節点になります。

良い質問だ。理由は主に2つある。第一に、複雑な幾何形状、特に曲面を持つ構造物をメッシュ分割する際、三角形要素は曲面への追従性が高い。しかし、三角形要素はせん断ロッキングの問題が顕著で、特に板・シェル構造の曲げ解析には向かない。そこで、三角形メッシュの面内方向の柔軟性と、プリズム要素の厚さ方向（高さ方向）の直線性を組み合わせることで、薄肉構造の曲げ挙動を精度良く捉えられるんだ。

板やシェルでは、厚さ方向のひずみ分布は線形と仮定される（Mindlin板理論など）。プリズム要素は、高さ方向に直線の辺を持っているため、この厚さ方向の線形変形モードを自然に表現できる。一方、四面体要素で板をモデル化すると、厚さ方向の節点配置が不規則になり、この線形分布を表現するのが難しく、過剰剛性（ロッキング）を引き起こしやすい。プリズム要素の形状関数は、面内方向（ξ, η）と厚さ方向（ζ）を分離して記述できる利点があるんだ。

標準的な設定では、一次要素（線形）で面内方向に1点、厚さ方向に2点のガウス積分点を使うことが多い。これは「減積分」と呼ばれ、せん断ロッキングを軽減する効果がある。AnsysのSolid186（20節点二次要素）のデフォルト設定は、完全積分で3x3x3=27点だが、シェル状の構造には「URI（Uniform Reduced Integration）」オプションを選び、面内2x2、厚さ方向2点の合計8点に減らすことが推奨される場合もある。過剰な減積分はゼロエネルギー変形（アワーミグラス）を生むリスクがあるから注意が必要だ。

はい、代表的なのが「Assumed Natural Strain（ANS）」法や「Enhanced Assumed Strain（EAS）」法だ。AbaqusのC3D6（一次プリズム）やC3D15（二次プリズム）要素は、デフォルトでせん断ロッキング対策が施されている。特に複合シェルをソリッド要素でモデル化する場合、Abaqusでは「非適合モード」を有効にした要素タイプを選ぶことで、曲げ変形の精度を大幅に向上させることができる。この場合、要素内部に追加的な変形自由度が導入される。

一般的に、同じ節点数であれば、プリズムメッシュの方が四面体メッシュより接続性が良く（各節点が接続する要素数が少なく）、結果として剛性マトリックスのバンド幅が小さくなる傾向がある。これは直接法ソルバーの計算時間短縮に寄与する。しかし、プリズムメッシュは要素のアスペクト比（厚さ方向の長さと面内の辺長の比）が極端に大きくなりやすい。アスペクト比が1000を超えるような薄い要素では、マトリックスの条件数が悪化し、反復法ソルバーの収束性が低下する。MSC Nastranでは、この問題を緩和するためにスケーリングオプションが用意されている。

航空機の複合材パネルや自動車のスポット溶接周辺の応力集中解析が良い例だ。まず、CADでの中面形状をインポートし、シェルメッシュ（主に三角形）を作成する。次に、このシェルメッシュを「押し出し」または「ソリッドメッシュ生成」機能を使って、厚さ方向に数層（例えば複合材なら1層積層ごとに1層）のプリズム要素を生成する。この時、厚さ方向の要素数は少なくとも2層、できれば3層以上が望ましい。1層だとせん断ひずみが一定と仮定されることになり、精度が落ちる。

これは非常に重要なポイントだ。一般的な経験則として、アスペクト比（面内辺長/要素厚さ）は10:1から20:1以内に収めるのが安全だ。航空宇宙分野の薄肉構造では100:1以上になることもあるが、その場合は要素公式や積分法の選択が結果に大きく影響する。Ansysのマニュアルでは、静解析では1000:1まで可能だが、固有値解析や動解析では100:1以下を推奨している。常にメッシュ感度解析を行い、応力や変位がアスペクト比に敏感でないことを確認する必要がある。

接続部で注意が必要だ。例えば、薄肉部をプリズム要素で、リブやブラケットなどの肉厚部を四面体要素でモデル化する場合、両者の界面で節点が一致していれば直接接続できる。しかし、プリズム要素の三角形面と四面体要素の三角形面を接続する場合は問題ないが、プリズムの四角形面と四面体を接続する場合は、四角形面を2つの三角形に分割する「Tria Connectivity」などのオプションをソフトウェアが内部で適用する。この時、局部剛性が変化する可能性があるので、応力集中が生じる接合部では、可能な限り均一な要素タイプ（全て二次要素など）を使うのがベストプラクティスだ。

各社で呼び名と実装が異なる。Ansys Mechanicalでは、一次要素が「SOLID185」、二次要素が「SOLID186」だ。Abaqus/Standardでは、一次完全積分が「C3D6」、一次減積分が「C3D6R」、二次要素が「C3D15」となる。COMSOL Multiphysicsでは、「プリズム」要素カテゴリがあり、一次から三次までのラグランジュ要素を選択できる。特にCOMSOLは、電磁界解析や流体解析でもプリズム境界層メッシュを多用する点が特徴的だ。

層ごとに材料と繊維配向角を定義する「複合ソリッド」機能では、各社が力を入れている。Ansys Composite PrepPost (ACP) は、プリプリポスト処理に特化しており、プリズム要素による層別メッシュ生成と破壊基準（Tsai-Wu, Puckなど）の評価が直感的に行える。Abaqus/CAEにも「Composite Layup」機能があり、C3D6RやC3D15要素を用いた解析が可能だ。MSC Nastranの「PCOMP」プロパティカードも業界標準の一つで、SOL 400などの非線形ソルバーと連携して複合材の損傷進展を追跡できる。

CFDでは「プリズム境界層メッシュ」が必須技術だ。壁面近くの流速勾配が急峻な領域（境界層）を、壁面に直交するように何層も重ねたプリズム要素で捕捉する。例えば、Star-CCM+の「プリズムレイヤーメッシャー」や、Fluent Meshingの「Inflation」機能がこれに当たる。第一層の厚さは無次元壁面距離y+の値（例えば1以下）から逆算して決定する。これにより、乱流モデル（SST k-ωなど）の精度が大幅に向上する。一方、流域内部は四面体やポリヘドロン要素で埋め、境界層メッシュと接続するハイブリッドメッシュが一般的だ。

プリズム要素で多いのは、極端に大きなアスペクト比と、押し出し方向が曲面の法線方向から大きくずれている場合だ。後者は、中面から押し出しメッシュを生成する際、複雑な曲面では隣接する要素の押し出し方向が不連続になり、要素が「捩じれた」形状になることが原因だ。対策としては、1) 面内メッシュを細かくして形状追従性を高める、2) 押し出し方向を「平均曲面法線」に統一する、3) どうしてもダメな局部は、二次要素（C3D15など）に切り替えて形状表現力を上げる、などがある。

これは古典的な問題だ。層間のせん断応力を一次プリズム要素（C3D6）の1層で評価すると、せん断ひずみが一定と仮定されるため、応力が過小評価される。少なくとも厚さ方向に3層以上の要素が必要で、二次要素（C3D15）なら2層でも改善される。さらに、Abaqusの「連続体シェル要素」（SC6R, SC8R）のように、厚さ方向の積分点を多く持つ特殊要素を使う方法もある。また、はく離自体をモデル化するには、層間に「コヒーシブ要素」（COH3D6など）を挿入するのが現在のベストプラクティスだ。この場合、コヒーシブ要素の厚さはゼロに近い仮想的な値とし、そのトラクション-分離則を定義する。

まず確認すべきは「メッシュ品質」だ。プリズムレイヤーの成長率（隣接層の厚さの比）が1.2を超えすぎていないか。急激に厚さが増えると数値拡散が大きくなる。次に、境界層メッシュと内部のコアメッシュ（四面体など）の接続部で、要素サイズが急激に変化していないか。サイズ比が5倍以上あると、数値的な「反射」が起きて収束を乱す。最後に、y+の値が使用する乱流壁面モデル（壁関数か低レイノルズ数モデルか）に適しているか。例えばSST k-ωモデルで低レイノルズ数条件を解くなら、y+は1以下が理想だが、これが達成できていないと物理的に誤った結果になる。