要素内補間

要素内補間（ようそないほかん）とは、数値解析において要素の各節点の既知量から、要素内の値を補間して求めることをいう。この要素内補間は、内挿とも呼ばれることがある。要素内補間は、例えば地図の等高線、CAD、CAE、CGなど、要素が使用される図形処理において、要素内の任意の位置の値を計算する際にも使用される。
与えられた節点情報のみ（要素情報は使用しない）から補間する手法もあるが、これらは本説明に含まれない。
要素には、線分（2節点）、3角形（3節点）、4面体（4節点）などがある。
線分要素内の補間.
全体座標系と局所座標系の関係.
線分要素内の点p は、節点p0, p1 により
と表せる。ここで局所座標"u" は 0 < "u" < 1 を満たす実数で、
である。直感的には、基点を "p"0 として、そこから "p"1 までの距離の比率 "u" で線分内の座標値 "p" を求めたことになる。
線形補間.
全体座標系と局所座標系の関係と同様に、節点 p0, p1 での各値を "C"0, "C"1 とし線形補間すると、点 "p" での値 "C" は
と表せる。
計算例.
節点p0 (5, 5, 0), p1 (10, 10, 0) で各節点の値がそれぞれ"C"0 = 10, "C"1 = 20 の場合に、線分の中点p (15/2, 15/2, 0) での値"C" は、
より
になる。
3角形要素内の補間.
局所座標系と全体座標系の変換.
点 p0 を原点とする、3角形で構成される局所座標系は、基底ベクトル
から得られ、全体座標系への変換行列は、点 p の局所座標系の座標を "up", "vp" 、全体座標系の座標を "xp", "yp", "zp" とすると以下の通りとなる。
または成分で表せば、
線形補間.
節点p0, p1, p2 での各値を"C"0, "C"1, "C"2とすると、要素内の点p での値"C" は
と表せる。ここで"up", "vp" は点p の局所座標系での座標である。
計算例.
節点座標がp0 (0, 0, 0), p1 (5, 0, 0), p2 (5, 5, 0) とし、各節点の既知量は、それぞれ"C"0 = 10, "C"1 = 20,"C"2 = 30とすると、
4面体要素内の補間.
局所座標系と全体座標系の変換.
4面体で構成される局所座標系は、基底ベクトル (eu", ev", "ew" ) から得られ、全体座標系への変換行列は、局所座標系の座標を"up", "vp", "wp" と置くと以下の通りとなる。
または成分で表せば、
ただし
である。
全体座標系から局所座標系への変換は、局所座標系から全体座標系への変換行列の逆行列を求めることで得られる。
線形補間.
節点p0, p1, p2, p3 での各値を"C"0, "C"1, "C"2, "C"3 とすると、点p の値"C" は
と表せる。ここで"up", "vp", "wp" は点p の局所座標系での座標である。
計算例.
節点座標がp0 (0, 0, 0), p1 (5, 0, 0), p2 (5, 5, 0), p3 (0, 0, 5)、各節点の既知量はそれぞれ"C"0 = 10, "C"1 = 20, "C"2 = 30, "C"3 = 40 とする。
このとき、重心位置pG での座標は、(5/2, 5/4, 5/4) で、全体座標系から局所座標系での座標を求めると、 ("u"G, "v"G, "w"G ) = (1/4, 1/4, 1/4)となる。これを補間式にあてはめると"C" = 25 となる。重心位置のため、平均値 ("C"0 + "C"1 + "C"2 + "C"3) / 4 と同じになる。