アイソパラメトリック写像シミュレーター 戻る
FEMシミュレーター

アイソパラメトリック写像シミュレーター — 4節点四辺形要素

基準正方形 (ξ,η)∈[-1,1]² から物理座標 (x,y) への写像とヤコビ行列式を可視化。コーナー節点を動かして、要素歪みと反転(detJ≤0)が起こる仕組みを学べます。

パラメータ設定
Δx₁
Δy₁
ξ
η

既定の標準四辺形は 4 隅 (0,0),(1,0),(1,1),(0,1)。Corner 1 のみを (Δx₁, Δy₁) だけ摂動します。

計算結果
0.550
物理座標 x
0.525
物理座標 y
0.2125
ヤコビ行列式 det J
85.0 %
歪み比 detJ/0.25
基準空間 (ξ,η) と物理空間 (x,y)

左=基準正方形 [-1,1]²/右=摂動した物理四辺形。赤丸=評価点 (ξ,η) と (x,y)。det J ≤ 0 で「要素反転」警告。

理論・主要公式

4 節点四辺形(Q4)要素の形状関数(双線形)。コーナー番号は 1=(-1,-1)、2=(+1,-1)、3=(+1,+1)、4=(-1,+1) に対応:

$$N_1=\tfrac{1}{4}(1-\xi)(1-\eta),\;N_2=\tfrac{1}{4}(1+\xi)(1-\eta),\;N_3=\tfrac{1}{4}(1+\xi)(1+\eta),\;N_4=\tfrac{1}{4}(1-\xi)(1+\eta)$$

物理座標への写像(座標と未知量に同じ形状関数を使う=アイソパラメトリック):

$$x=\sum_{i=1}^{4} N_i(\xi,\eta)\,x_i,\qquad y=\sum_{i=1}^{4} N_i(\xi,\eta)\,y_i$$

ヤコビ行列とヤコビ行列式(要素歪みインジケータ):

$$J=\begin{bmatrix}\partial x/\partial\xi & \partial y/\partial\xi\\ \partial x/\partial\eta & \partial y/\partial\eta\end{bmatrix},\qquad dx\,dy=|\det J|\,d\xi\,d\eta$$

標準正方形(Δx=Δy=0)では det J = 0.25。det J ≤ 0 では要素が反転しており、剛性行列の数値積分が破綻します。

アイソパラメトリック写像シミュレーターとは

🙋
FEM の教科書で「アイソパラメトリック要素」って出てきたんですけど、何が「同じ(iso)」なんですか?
🎓
ざっくり言うと、「形を補間する関数」と「変位を補間する関数」に同じものを使う、っていう意味だ。基準になる正方形 $(\xi,\eta)\in[-1,1]^2$ を用意しておいて、形状関数 $N_i$ で実際の四辺形に写す。同時に節点の変位も同じ $N_i$ で要素内部に補間する。1 つの基準要素で任意形状の要素を作れるのが嬉しい点だよ。上のシミュレーターで Corner 1 を動かしてみて、右側の物理四辺形が変形するのを見てごらん。
🙋
「ヤコビ行列式」のカードがあるんですが、これは何を表してるんですか?
🎓
ヤコビ行列式は「基準空間の小さな面積が、物理空間でどれだけ拡大されるか」の比率だ。式で言うと $dx\,dy = |\det J|\,d\xi\,d\eta$。剛性行列を計算するときの面積積分で必ず出てくる。Δx=Δy=0 の理想的な単位正方形だと det J = 0.25 ぴったり。Corner 1 を動かしていくと、この値が変動する。歪み比カードで標準値 0.25 の何倍かが分かるよ。
🙋
Δx₁ や Δy₁ をマイナス方向にぐっと動かすと「要素反転」って赤字が出ました!これは何が起こってるんですか?
🎓
それが det J ≤ 0 になった状態だ。物理空間で四辺形が「裏返って」しまっている。例えば Corner 1 を Corner 3 の向こう側まで動かすと、辺が交差して蝶ネクタイ型になる。こうなると形状関数による写像が一対一でなくなり、剛性行列の積分が破綻する。実務の大変形解析で「ソルバーが落ちました」というときの最頻原因の一つで、リメッシュや ALE で対処するんだ。
🙋
背景に格子模様が見えますが、あれは何ですか?
🎓
基準空間で $\xi$ と $\eta$ を 0.05 刻みで動かしたときの曲線群を、物理空間に写したものだよ。標準四辺形だと真っ直ぐな格子だけど、Corner 1 を動かすと格子線が湾曲する。これがヤコビ行列が場所ごとに違う値を取るということを視覚的に表している。FEM では数値積分(ガウス積分)でこの効果を取り込むんだ。

よくある質問

座標補間に使う関数の次数と、変位補間に使う関数の次数の関係で 3 つに分類されます。両者が同じ次数のときが「アイソパラメトリック」(標準的な FEM 要素はこれ)、座標補間の方が低次なら「サブパラメトリック」、高次なら「スーパーパラメトリック」です。アイソパラメトリックは数値計算上の都合がよく、剛体並進・回転の表現が自然に成立するため、商用 FEM ソルバーの大半が採用しています。
Q4 要素では基準座標で $(\pm 1/\sqrt{3},\pm 1/\sqrt{3})$ の 4 点に置く 2×2 ガウス積分が標準です(完全積分)。曲げを正確に評価できる反面、せん断ロッキングが発生しやすい問題があります。せん断ロッキング対策として、せん断項のみ 1×1 で評価する「縮約積分」や、選択的縮約積分、$\bar{B}$ 法などが用いられます。要素の歪みが大きい領域では、積分点の物理座標に対する写像が不正確になるため、メッシュ品質の管理が重要です。
商用 FEM ソルバーでは、要素歪み比(Jacobian Ratio)として、要素内の最大 det J と最小 det J の比を使うことが多く、一般に 1〜10 程度なら良好、10〜40 程度で警告、40 以上または負値(要素反転)でエラー停止という扱いです。アスペクト比(最大辺/最小辺)も別途チェックされ、線形要素なら 5〜10 以内、二次要素なら 20 程度までが目安です。歪んだ要素は応力勾配が大きい領域で精度を大きく落とすため、応力集中部では特に注意が必要です。
はい、まったく同じ枠組みで 8 節点六面体要素(H8)が定義できます。基準立方体 $(\xi,\eta,\zeta)\in[-1,1]^3$ から物理空間への写像を、トリリニア形状関数 $N_i = \tfrac{1}{8}(1\pm\xi)(1\pm\eta)(1\pm\zeta)$ で行い、ヤコビ行列は 3×3 になります。ヤコビ行列式が体積比を表す点も同じです。三角形・四面体要素では面積座標・体積座標を使った別形式になりますが、考え方は共通です。

実世界での応用

商用 FEM ソルバーの基盤:Abaqus、ANSYS、Nastran、LS-DYNA など主要な商用 FEM ソルバーの大半が、アイソパラメトリック要素を採用しています。線形 4 節点四辺形(Q4/CPS4・S4 等)から二次 8 節点(Q8/CPS8 等)、3 次元の六面体・テトラまで、要素の中身はすべてこの理論で組み立てられています。任意形状の構造をメッシュ化できる柔軟性は、現代の CAE が成立する前提となっています。

大変形解析・成形シミュレーション:金属プレス成形・ゴムシール圧縮・衝突解析など、要素自体が大きく変形する解析では、要素歪み比(Jacobian Ratio)と det J の正値性を逐次チェックすることが必須です。det J が小さくなりすぎたら適応的にメッシュを切り直す「リメッシュ」、または基準格子と物理格子を分離する ALE(Arbitrary Lagrangian-Eulerian)法が用いられます。

等幾何解析(IGA)への発展:アイソパラメトリックの考え方を発展させたのが、CAD で使われる NURBS 基底をそのまま要素として使う等幾何解析(Isogeometric Analysis, IGA)です。CAD と CAE のメッシュ変換ロスを排除でき、薄板・シェル解析で高い精度が得られます。基本となる「同じ関数で形と未知量を補間する」という哲学は、Q4 要素から IGA まで一貫しています。

教育・FEM 入門の標準題材:形状関数・ヤコビ行列・数値積分という FEM の核心概念を、目に見える形で理解できる最小単位が Q4 アイソパラメトリック要素です。多くの教科書で最初に登場し、コードを実装する演習でも最初の課題として扱われます。本シミュレーターはまさにこの教育題材を可視化したものです。

よくある誤解と注意点

最も多い誤解は、「四辺形要素は常に正方形に近い形でなければならない」と考えてしまうことです。実際には台形・斜方形・任意の凸四辺形でも、det J が正で十分な大きさを保てば計算は成立します。むしろ複雑な実形状をメッシュで埋めるには、ある程度の歪みを許容しないと不可能です。シミュレーターで Δx₁=0.3, Δy₁=0.2 程度に設定すると、台形に近い形でも det J の歪み比が 60〜80 % 程度を保ち、十分実用的な要素であることが分かります。

次に多いのが、「det J が大きいほど良い要素だと思い込む」ことです。det J は単に基準空間と物理空間の面積比であり、要素のサイズに比例して増えるだけです。重要なのは det J の絶対値ではなく、要素内の積分点間での det J の比(Jacobian Ratio)と、最小 det J が正であることです。本シミュレーターの「歪み比」カードは標準正方形 0.25 を基準にした相対値で、100 % 付近が理想・40 % 以下で警告・0 % 以下で反転、と読み取れます。

最後に、「形状関数の総和 ΣN_i が 1 になる条件を見落とす」点に注意してください。ΣN_i = 1 は剛体並進を正しく表現するための必須条件で、これが破れると要素は剛体運動でも内部応力を発生させてしまいます。$N_1+N_2+N_3+N_4 = \tfrac{1}{4}[(1-\xi)(1-\eta)+(1+\xi)(1-\eta)+(1+\xi)(1+\eta)+(1-\xi)(1+\eta)] = 1$ が任意の $(\xi,\eta)$ で成立することを式変形で確認できます。新しい要素タイプを実装する際は、この条件と「節点位置で $N_i$ がクロネッカーのデルタになる」条件の両方を必ず検証してください。