1D有限要素 形状関数シミュレーター

構造解析ソルバーの中核：Ansys、Abaqus、Nastran など、あらゆる商用FEMコードは形状関数を使って要素剛性行列を組み立てます。要素剛性は K = ∫Bᵀ D B dV で計算され、ひずみ-変位行列 B は形状関数の微分そのものです。つまり「どんな形状関数を選ぶか」が、その要素の精度・収束性・適用範囲を直接決めています。形状関数を理解することは、ソルバーがブラックボックスでなくなる第一歩です。

メッシュ品質と要素次数の選択：実務では「線形要素を細かく切るか、二次要素を粗く切るか」という判断が常に発生します。応力を見たい解析では、二次要素のほうが曲げや応力勾配を正確に捉えるため、節点数が同じでも精度が高いことが多いです。一方、接触解析や陽解法の衝突解析では線形要素が好まれます。形状関数の次数の違いを体感しておくと、この使い分けの理由が腑に落ちます。

アイソパラメトリック要素とヤコビアン：自然座標 ξ で書いた形状関数は、座標写像 x(ξ)=ΣN_i·x_i にも使われます。実座標での微分は du/dx = (du/dξ)/(dx/dξ) のようにヤコビアン dx/dξ を介して計算されます。本ツールが表示する du/dξ は、この変換の前段にあたる量です。歪んだ要素でヤコビアンが負やゼロになると要素が破綻するため、メッシュ品質チェックの理論的根拠にもなっています。

有限要素法以外への展開：「節点値を重み関数で内挿する」という発想は、等幾何解析（IGA）のNURBS基底、無メッシュ法のMLS近似、CFDの不連続ガラーキン法など、現代の数値解析の多くに受け継がれています。単位分割という性質は、これらの手法に共通する設計原理です。1次元の単純な形状関数で原理をつかんでおくと、より高度な手法の理解が一気に進みます。

まず多いのが、「形状関数は補間の式であって、解そのものではない」という点の見落としです。形状関数はあくまで「節点値が決まったら要素内をどう埋めるか」を定めるルールであり、節点値 u_i は連立方程式 Ku=f を解いて初めて得られます。形状関数の選び方は解の精度を左右しますが、形状関数だけ眺めても物理的な答えは出ません。本ツールは節点値を手で与えて補間の挙動だけを見るものであり、実際のFEM解析では節点値が未知数である点に注意してください。

次に、「次数を上げれば必ず精度が上がる」と思い込むこと。二次要素は曲がった分布を表せますが、節点値が極端に振れている場合や荷重が局所的な場合、高次の形状関数はかえって振動（オーバーシュート）を起こすことがあります。また二次要素は中間節点を持つため、メッシュの不整合や中間節点のずれに敏感です。次数を上げるよりメッシュを適切に細分化したほうが安定することも多く、要素次数とメッシュ密度はセットで検討すべきです。

最後に、「ΣN=1 や ξ∈[−1,1] は当たり前」と軽視すること。単位分割が成り立たない形状関数を自作すると、剛体並進だけでニセのひずみが生じ、収束しない解になります。また自然座標は要素内で必ず −1〜+1 に収まらねばならず、座標写像が破綻して ξ が範囲外になると形状関数の意味が失われます。歪んだ要素でヤコビアンが負になる不具合は、まさにこの「ξ 空間と実空間の対応が壊れた」状態です。形状関数の素直な性質こそが、有限要素法の信頼性を支えています。