3角形のケース: 各節点p0, p1, p2 での値が既知 (10, 20, 30) の場合に、要素内の点 p での値はいくつになるのだろう?これを求めるのが要素内補間である。 要素内補間(ようそないほかん)とは、数値解析において要素の各節点の既知量から、要素内の値を補間して求めることをいう。この要素内補間は、内挿とも呼ばれることがある。要素内補間は、例えば地図の等高線、CAD、CAE、CGなど、要素が使用される図形処理において、要素内の任意の位置の値を計算する際にも使用される。
与えられた節点情報のみ(要素情報は使用しない)から補間する手法もあるが、これらは本説明に含まれない。
要素には、線分(2節点)、3角形(3節点)、4面体(4節点)などがある。
線分要素内の補間
線分要素のケース: p0, p1 の各節点値が既知(10, 20)の場合に、線分内のp の値を求める。 全体座標系と局所座標系の関係
線分要素内の点p は、節点p0, p1 により
![{\displaystyle {\boldsymbol {p}}={\boldsymbol {p}}_{0}+({\boldsymbol {p}}_{1}-{\boldsymbol {p}}_{0})u}](https://wikimedia.org/api/rest_v1/media/math/render/svg/95df9b1073dfa433f41eaf6d46648991c322b0b7)
と表せる。ここで局所座標u は 0 < u < 1 を満たす実数で、
![{\displaystyle u={\frac {|p-p_{0}|}{|p_{1}-p_{0}|}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3e28574bde8a545337348241258c98e9392a28e7)
である。直感的には、基点を p0 として、そこから p1 までの距離の比率 u で線分内の座標値 p を求めたことになる。
- u = 0 の場合には点 p は p0 を示し、u = 1 の場合には p1 を示す。
- 前述の表現は、直線のパラメトリックまたは、媒介変数u による定義と呼ばれることもある。
線形補間
全体座標系と局所座標系の関係と同様に、節点 p0, p1 での各値を C0, C1 とし線形補間すると、点 p での値 C は
![{\displaystyle C=C_{0}+(C_{1}-C_{0})u}](https://wikimedia.org/api/rest_v1/media/math/render/svg/43f354610bb14c443c38b7256499c8f7bdaba587)
と表せる。
計算例
節点p0 (5, 5, 0), p1 (10, 10, 0) で各節点の値がそれぞれC0 = 10, C1 = 20 の場合に、線分の中点p (15/2, 15/2, 0) での値C は、
![{\displaystyle u={\frac {|{\boldsymbol {p}}-{\boldsymbol {p}}_{0}|}{|{\boldsymbol {p}}_{1}-{\boldsymbol {p}}_{0}|}}=0.5}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4713f39c6daf775738144259ed2c95ce1760b96a)
より
![{\displaystyle C=C_{0}+(C_{1}-C_{0})u=15}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8970a6515c263813e5bd5d6f2453ada9606a07ec)
になる。
3角形要素内の補間
面要素のケース: 接点 p0, p1, p3 での値が既知でそれぞれ (10, 20, 30) の場合に、面内の点 p の値を求める。
局所座標系上での位置座標を求めれば点 p での値が求まる。 局所座標系と全体座標系の変換
点 p0 を原点とする、3角形で構成される局所座標系は、基底ベクトル
![{\displaystyle {\begin{aligned}{\boldsymbol {e}}_{u}&={\boldsymbol {p}}_{1}-{\boldsymbol {p}}_{0}\\{\boldsymbol {e}}_{v}&={\boldsymbol {p}}_{2}-{\boldsymbol {p}}_{0}\\\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/542f5667ac86df9f4f11cce92378f154b4e6ed6c)
から得られ、全体座標系への変換行列は、点 p の局所座標系の座標を up, vp 、全体座標系の座標を xp, yp, zp とすると以下の通りとなる。
![{\displaystyle {\boldsymbol {p}}={\boldsymbol {p}}_{0}+u_{p}{\boldsymbol {e}}_{u}+v_{p}{\boldsymbol {e}}_{v}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/83897ea0bbe500b21e09ff6f27b021cd0fe8fe4a)
または成分で表せば、
![{\displaystyle \left({\begin{array}{c}x_{p}\\y_{p}\\z_{p}\\\end{array}}\right)=\left({\begin{array}{c}x_{p0}\\y_{p0}\\z_{p0}\end{array}}\right)+\left({\begin{array}{ccc}x_{eu}&x_{ev}\\y_{eu}&y_{ev}\\z_{eu}&z_{ev}\end{array}}\right)\left({\begin{array}{c}u_{p}\\v_{p}\end{array}}\right)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/acc7876f7264d713f1d984834d6440c3bc7e74c8)
線形補間
節点p0, p1, p2 での各値をC0, C1, C2とすると、要素内の点p での値C は
![{\displaystyle C=C_{0}+u_{p}(C_{1}-C_{0})+v_{p}(C_{2}-C_{0})}](https://wikimedia.org/api/rest_v1/media/math/render/svg/19fde36a2d4619dd8409b76d5a1e2fb3b705a331)
と表せる。ここでup, vp は点p の局所座標系での座標である。
- up = 0, vp = 0 の場合には、C = C0 (p0 の値)を示す。
- up = 1, vp = 0 の場合には、C = C1 (p1 の値)を示す。
- up = 0, vp = 1 の場合には、C = C2 (p2 の値)を示す。
- up, vp ≥ 0 かつ up + vp ≤ 1 の場合には、点p は要素の内部に存在する。
計算例
重心位置と p1, p2 の中点位置(赤色は既知量) 節点座標がp0 (0, 0, 0), p1 (5, 0, 0), p2 (5, 5, 0) とし、各節点の既知量は、それぞれC0 = 10, C1 = 20,C2 = 30とすると、
- 重心位置での値
- 重心位置pG の座標は (10/3, 5/3, 0) で、全体座標系から局所座標系での座標を求めると、(uG, vG ) = (1/3, 1/3) となる。これを補間式にあてはめると C G = 20 となる。重心位置のため、平均値 (C0 + C1 + C2 ) / 3と同じである。
- p1 とp2 の中点位置座標での値
- p1 とp2 の中点位置p12 の全体座標系での座標は (5, 5/2, 0) で、全体座標系から局所座標系での座標 (u12, v12 ) を求めると、(1/2, 1/2) となる。これを補間式にあてはめると C12 = 25 となる。
4面体要素内の補間
4面体要素のケース: 節点p0, p1, p3, p4 での値が既知 (10, 20, 30, 40) の場合に、要素内の点p の値を求める。 局所座標系と全体座標系の変換
4面体で構成される局所座標系は、基底ベクトル (eu, ev, ew ) から得られ、全体座標系への変換行列は、局所座標系の座標をup, vp, wp と置くと以下の通りとなる。
![{\displaystyle {\boldsymbol {p}}={\boldsymbol {p}}_{0}+u_{p}{\boldsymbol {e}}_{u}+v_{p}{\boldsymbol {e}}_{v}+w_{p}{\boldsymbol {e}}_{w}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/7a71e54889283595932d9703d03b72c99f310be6)
または成分で表せば、
![{\displaystyle \left({\begin{array}{c}x_{p}\\y_{p}\\z_{p}\\\end{array}}\right)=\left({\begin{array}{c}x_{p0}\\y_{p0}\\z_{p0}\end{array}}\right)+\left({\begin{array}{ccc}x_{eu}&x_{ev}&x_{ew}\\y_{eu}&y_{ev}&y_{ew}\\z_{eu}&z_{ev}&z_{ew}\end{array}}\right)\left({\begin{array}{c}u_{p}\\v_{p}\\w_{p}\end{array}}\right)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e1cd322a384ae90349dc28638d1fc75804c1cd9b)
ただし
![{\displaystyle {\begin{aligned}&{\boldsymbol {e}}_{u}={\boldsymbol {p}}_{1}-{\boldsymbol {p}}_{0}\\&{\boldsymbol {e}}_{v}={\boldsymbol {p}}_{2}-{\boldsymbol {p}}_{0}\\&{\boldsymbol {e}}_{w}={\boldsymbol {p}}_{3}-{\boldsymbol {p}}_{0}\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8f04bd8a999b7333140f72d955b393fdf0b1845e)
である。
全体座標系から局所座標系への変換は、局所座標系から全体座標系への変換行列の逆行列を求めることで得られる。
![{\displaystyle \left({\begin{array}{c}u_{p}\\v_{p}\\w_{p}\end{array}}\right)=\left({\begin{array}{ccc}x_{eu}&x_{ev}&x_{ew}\\y_{eu}&y_{ev}&y_{ew}\\z_{eu}&z_{ev}&z_{ew}\end{array}}\right)^{-1}\left(\left({\begin{array}{c}x_{p}\\y_{p}\\z_{p}\\\end{array}}\right)-\left({\begin{array}{c}x_{p0}\\y_{p0}\\z_{p0}\end{array}}\right)\right)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8483d56295f7b18d594bbaed0bd991aec1d54ba6)
線形補間
節点p0, p1, p2, p3 での各値をC0, C1, C2, C3 とすると、点p の値C は
![{\displaystyle C=C_{0}+u_{p}(C_{1}-C_{0})+v_{p}(C_{2}-C_{0})+w_{p}(C_{3}-C_{0})}](https://wikimedia.org/api/rest_v1/media/math/render/svg/9970de8dd409b374a165b6376ad266f0c9c96826)
と表せる。ここでup, vp, wp は点p の局所座標系での座標である。
- u = 0, v = 0, w = 0 の場合には、C = C0 (p0 の値)を示す。
- u = 1, v = 0, w = 0 の場合には、C = C1 (p1 の値)を示す。
- u = 0, v = 1, w = 0 の場合には、C = C2 (p2 の値)を示す。
- u = 0, v = 0, w = 1 の場合には、C = C3 (p3 の値)を示す。
- u, v, w ≥ 0 かつ u + v + w ≤ 1 の場合には、点p は要素の内部に存在する。
計算例
節点座標がp0 (0, 0, 0), p1 (5, 0, 0), p2 (5, 5, 0), p3 (0, 0, 5)、各節点の既知量はそれぞれC0 = 10, C1 = 20, C2 = 30, C3 = 40 とする。
このとき、重心位置pG での座標は、(5/2, 5/4, 5/4) で、全体座標系から局所座標系での座標を求めると、 (uG, vG, wG ) = (1/4, 1/4, 1/4)となる。これを補間式にあてはめるとC = 25 となる。重心位置のため、平均値 (C0 + C1 + C2 + C3) / 4 と同じになる。
関連項目