p適応 — CAE用語解説

有限要素法における要素の「多項式次数(p)」を適応的に変化させる手法です。h適応がメッシュサイズ(h)を細かくするのに対し、p適応は同じメッシュ形状のまま、要素内の近似関数の次数を上げて精度を高めます。例えば、線形要素(p=1)から二次要素(p=2)、三次要素(p=3)へと上げていきます。

支配方程式の近似解の滑らかさが向上します。例えば、一次元の変位u(x)を線形要素で近似すると折れ線ですが、二次要素では放物線、三次要素では三次曲線で近似できます。特に応力やひずみは変位の微分なので、

はい。誤差の収束率は、メッシュサイズをh、多項式次数をpとすると、エネルギー・ノルムで

主に要素剛性マトリックスの構成と数値積分（ガウス積分）が変わります。高次要素では、要素内の節点（または自由度）の数が増えるため、形状関数が複雑になり、剛性マトリックスのサイズとバンド幅が大きくなります。例えば、二次四面体要素（10節点）は線形四面体要素（4節点）に比べ、一要素あたりの自由度が2.5倍になります。

特性が異なります。h適応では未知数の数（メッシュ数）が爆発的に増えますが、マトリックスはスパースで条件数は悪化しません。p適応では全体の未知数は緩やかに増加しますが、一要素あたりの自由度が増え、マトリックスのバンド幅が広がり、条件数が悪化します。特にpを5以上にすると、反復ソルバーの収束性が極端に悪くなることがあり、前処理技術が重要になります。

【向いている】幾何形状が単純で材料が均質、解が滑らかな問題。例えば、弾性体の遠方応力場、規則的な熱伝導、高周波電磁界解析などです。【向かない】不連続や特異性がある問題。具体的にはき裂先端、鋭い角での応力集中、衝撃波（CFD）、異材界面などです。これらの領域ではh適応で局所的にメッシュを細かくする方が効率的です。

誤差指標に基づいて自動判定するのが一般的です。例えば、Zienkiewicz-Zhuの誤差推定法で各要素の誤差を評価し、誤差が大きい要素の次数を上げます。Ansys Mechanicalでは「Adaptive p-Method」として実装されており、目標精度（例えばエネルギー・ノルム誤差2%）を設定すると、自動的に次数分布を決定します。最初は全域p=2から始め、最大p=8まで上げることもあります。

必須です。特に初めてのケースでは、p適応とh適応（メッシュ収束）の両方で同じ物理量（最大応力、固有値など）が収束することを確認します。また、pを上げても結果が変わらなくなった点（p収束）を見極めます。実務では、ASME V&V 20規格に沿って、計算結果の不確かさを評価する際、異なる離散化手法による検証を推奨しています。

主に「Ansys Mechanical」の静解析/モーダル解析で利用可能です。「Analysis Settings」内で「Adaptive p-Method」を有効にし、目標精度を設定します。最大多項式次数はデフォルトで8まで設定可能です。ただし、非線形解析や接触問題には対応していない場合が多いので、マニュアルで確認が必要です。

Abaqus/Standardには「p適応」という名称の自動機能はあまりありません。ただし、ユーザーが手動で要素ライブラリから高次要素（例えば二次要素C3D20、三次要素は標準ではない）を選択してメッシュを切ることはできます。自動誤差推定に基づく次数適応は、AnsysやCOMSOLに比べると前面に出ていない印象です。

COMSOL Multiphysicsは「メッシュ適応化」の機能の中で、h適応とp適応を組み合わせた「hp適応」を強力にサポートしています。特に波動伝搬や固有値問題で有効です。物理場インターフェースごとに離散化の次数を個別に設定でき（例えば構造は二次、電磁気は三次）、「メッシュ」研究ステップで誤差指標に基づき自動で次数を調整して収束させられます。

Code_Asterでは、高次要素（二次、三次）の利用は可能ですが、解析中に自動で次数を変化させる「適応」機能は限定的です。CalculiXも同様で、主に手動での要素選択になります。オープンソースで本格的なhp-FEMを実装しているのは、例えば「NGSolve」のような学術志向の強いコードです。産業用途では商用ソフトの機能が先行しています。

それは「ギブス現象」の可能性が高いです。解に急峻な勾配や不連続がある領域で、高次多項式で無理に近似しようとすると、過剰振動が発生します。対策は、その領域ではpを上げずにh適応でメッシュを細かくする（hp適応の考え方）、または「フィルタリング」や「投影法」で結果を平滑化します。Ansysでは、応力結果を「ノード平均」ではなく「要素平均」で表示すると見かけ上はましになります。

剛性マトリックスの条件数が極端に悪化しているためです。高次形状関数は互いに直交性が低く、数値的にも線形従属に近くなり、マトリックスが特異に近づきます。対策としては、1) ソルバーを直接法（スパース直接ソルバー）に切り替える、2) 反復ソルバーを使う場合は、前処理（例えば不完全LU分解）の許容値を厳しくする、3) 積分点数を増やして数値積分精度を上げる、などがあります。

誤差指標が局所的に非常に大きいことを示しています。まず、その要素の幾何形状（アスペクト比が極端に悪いなど）や材料設定に問題がないか確認してください。また、その場所に特異性（鋭い角、集中荷重点）がある場合は、pを上げても誤差は減らないので、設計上その形状をフィレットするか、その周辺のみh適応でメッシュを細分化するのが現実的な対策です。ソフトのバグではなく、問題の性質に起因することが多いです。

主に二つです。1) 要素剛性マトリックスの生成：高次要素では数値積分のコストが膨大になります。2) ソルバー：反復ソルバーの収束回数が増大します。プロファイリングツールで確認し、1)が支配的なら積分点数を必要最小限に減らす（ただし精度低下に注意）、2)が支配的なら直接法ソルバーに切り替えることを検討します。メモリ消費も激増するので、計算リソースの見積もりは必須です。