フルネ-セレの公式の導出|空間曲線の曲率と捩率の定義と公式

ベクトル解析
ベクトル解析

3次元ユークリッド空間$\R^3$上の滑らかな曲線$C:\m{r}=\m{r}(t)$に対して

  • 「進む向き」を表す単位接ベクトル$\m{v}_1$
  • 「曲がる向き」を表す主法線ベクトル$\m{v}_2$
  • 「ねじれる向き」を表す従法線ベクトル$\m{v}_3$

を考えることができます.

これらのベクトルは正規直交基底をなしており,${\m{v}_1}’$, ${\m{v}_2}’$, ${\m{v}_3}’$を$\m{v}_1$, $\m{v}_2$, $\m{v}_3$の線形結合で表す公式をフルネ-セレ(Frenet-Serret)の公式(フレネ-セレの公式)といいます.

「フレネ-セレの公式」と言われることも多いですが,”Frenet”はフランス人名なので「フルネ」の方が本来の発音に近いです.なお,”Serret”もフランス人名で,こちらは「セレ」で本来の発音に近いです.

この記事では

  • 空間曲線の基礎知識
  • フルネ-セレの公式と曲率・捩率
  • フルネ-セレの公式の証明

を順に解説します.

空間曲線の基礎知識

以下,$\alpha<\beta$とします.

基本的に閉区間$[\alpha,\beta]$上で定義された連続関数を考えますが,導関数が関わる場合は開区間$(\alpha,\beta)$上で考えます.

空間曲線の考え方と定義

曲線とは直観的には「ちぎれていない曲線」のことですが,もう少しきちんと定義しましょう.

時刻$t=\alpha$から時刻$t=\beta$まで連続的に点が動く状況を考え,時刻$t$での点の位置を$\m{r}(t)$と表すことにすると,この$\m{r}$は閉区間$[\alpha,\beta]$で定義された連続写像と考えることができますね.

閉区間[α,β]上の連続写像r
連続写像$[\alpha,\beta]\to\R^3$は$\R^3$上の曲線と捉えられる

このことを踏まえて,次のように空間曲線を定義します.

1変数ベクトル値関数$\m{r}:[\alpha,\beta]\to\R^3$を空間曲線という.

この定義で$[\alpha,\beta]$は1次元ユークリッド空間$\R$の部分位相空間,$\R^3$は3次元ユークリッド空間です.つまり,$[\alpha,\beta]$, $\R^3$はどちらもいわゆる「直線距離」を備えた距離空間です.

空間曲線の導関数

空間曲線$\m{r}:(\alpha,\beta)\to\R^3$の各成分が微分可能であるとき,$\m{r}$は微分可能であるといいます.

空間曲線$\m{r}:[\alpha,\beta]\to\R^3$を

\begin{align*}\m{r}=\bmat{r_1\\r_2\\r_3},\quad
r_i:[\alpha,\beta]\to\R\quad(i=1,2,3)\end{align*}

と表すとき,$\m{r}$が(開区間$(\alpha,\beta)$上)微分可能であるとは,$r_1,r_2,r_3$がいずれも(開区間$(\alpha,\beta)$上)微分可能であることをいう.また,このとき$\m{r}$の導関数$\m{r}’$を

\begin{align*}\m{r}’=\bmat{r’_1\\r’_2\\r’_3}\end{align*}

で定める.

空間曲線の内積の微分公式

[内積の微分公式]微分可能な空間曲線$\m{r},\m{s}:[\alpha,\beta]\to\R^3$を考える.任意の$t\in(\alpha,\beta)$に対して,

\begin{align*}&\frac{d\anb{\m{r},\m{s}}}{dt}(t)=\anb{\m{r}'(t),\m{s}(t)}+\anb{\m{r}(t),\m{s}'(t)}\end{align*}

が成り立つ.ただし,$\anb{\cdot,\cdot}$は通常の内積(標準内積)である.

通常の関数の積の微分公式$(fg)’=f’g+fg’$と同様の形をしているので覚えやすいですね.

$\m{r}=\bmat{r_1\\r_2\\r_3}$, $\m{s}=\bmat{s_1\\s_2\\s_3}$とすると,内積の定義より

\begin{align*}\anb{\m{r},\m{s}}=r_1s_1+r_2s_2+r_3s_3=\sum_{i=1}^{3}r_is_i\end{align*}

なので,積の微分公式を用いて

\begin{align*}\od{\anb{\m{r},\m{s}}}{t}
&=\sum_{i=1}^{3}\bra{r’_is_i+r_is’_i}
\\&=\sum_{i=1}^{3}r’_is_i+\sum_{i=1}^{3}r_is’_i
\\&=\anb{\m{r}’,\m{s}}+\anb{\m{r},\m{s}’}\end{align*}

が成り立つ.

空間曲線の外積の微分公式

[外積の微分公式]微分可能な空間曲線$\m{r},\m{s}:[\alpha,\beta]\to\R^3$を考える.任意の$t\in(\alpha,\beta)$に対して,

\begin{align*}\frac{d(\m{r}\times \m{s})}{dt}(t)=\m{r}'(t)\times \m{s}(t)+\m{r}(t)\times \m{s}'(t)\end{align*}

が成り立つ.ただし,$\times$は$\R^3$上の外積である.

内積のときと同じく,外積でも通常の関数の積の微分公式$(fg)’=f’g+fg’$と同様の形をしているので覚えやすいですね.

$\m{r}=\bmat{r_1\\r_2\\r_3}$, $\m{s}=\bmat{s_1\\s_2\\s_3}$とすると,外積の定義より

\begin{align*}\m{r}\times \m{s}=\bmat{r_2s_3-r_3s_2\\r_3s_1-r_1s_3\\r_1s_2-r_2s_1}\end{align*}

なので,各成分で積の微分公式を用いて

\begin{align*}\od{(\m{r}\times \m{s})}{t}
&=\bmat{(r’_2s_3+r_2s’_3)-(r’_3s_2+r_3s’_2)\\(r’_3s_1+r_3s’_1)-(r’_1s_3+r_1s’_3)\\(r’_1s_2+r_1s’_2)-(r’_2s_1+r_2s’_1)}
\\&=\bmat{r’_2s_3-r’_3s_2\\r’_3s_1-r’_1s_3\\r’_1s_2-r’_2s_1}+\bmat{r_2s’_3-r_3s’_2\\r_3s’_1-r_1s’_3\\r_1s’_2-r_2s’_1}
\\&=\m{r}’\times \m{s}+\m{r}\times \m{s}’\end{align*}

が従う.

補題(ノルムが一定のベクトル)

次の補題はベクトル解析・曲線曲面論ではよく用いられるので,ここで準備しておきます.

微分可能な空間曲線$\m{r}:[\alpha,\beta]\to\R^3$のノルム$\|\m{r}\|$が一定なら,任意の$t\in(\alpha,\beta)$に対して

\begin{align*}\anb{\m{r}(t),\m{r}'(t)}=0\end{align*}

が成り立つ.

$\|\m{r}\|$が一定なら,内積の微分公式と実内積の対称性より

\begin{align*}0&\equiv(\|\m{r}\|^2)’=\anb{\m{r},\m{r}}’
\\&=\anb{\m{r}’,\m{r}}+\anb{\m{r},\m{r}’}=2\anb{\m{r},\m{r}’}\end{align*}

が分かる.よって,両辺を2で割って$0\equiv\anb{\m{r},\m{r}’}$が従う.





フルネ-セレの公式と曲率・捩率

以下,空間曲線$\m{r}:[\alpha,\beta]\to\R^3$は

  1. $\m{r}$は$C^\infty$級($\m{r}$は何回でも微分可能)
  2. $\m{r}’\neq\m{0}$($\m{r}$が動く速さは0にならない)
  3. $\m{r}^{\prime\prime}\times \m{r}’\neq\m{0}$(曲線の軌跡は曲がっている)

を満たすとします.

単位接ベクトル$\m{v}_1$・主法線ベクトル$\m{v}_2$・従法線ベクトル$\m{v}_3$の定義

空間曲線$C:\m{r}=\m{r}(t)$に対し,

  • 単位接ベクトル$\m{v}_1$:曲線$C$の「進む向き」
  • 主法線ベクトル$\m{v}_2$:曲線$C$の「曲がる向き」
  • 従法線ベクトル$\m{v}_3$:曲線$C$の「ねじれる向き」

を次で定義します.

空間曲線$C:\m{r}=\m{r}(t)$に対し,$\m{v}_1,\m{v}_2,\m{v}_3:(\alpha,\beta)\to\R^3$を

\begin{align*}\m{v}_1:=\frac{\m{r}’}{\|\m{r}’\|},\quad
\m{v}_2:=\frac{{\m{v}_1}’}{\|{\m{v}_1}’\|},\quad
\m{v}_3:=\m{v}_1\times \m{v}_2\end{align*}

で定義し,$\m{v}_1$を$C$の単位接ベクトル,$\m{v}_2$を$C$の主法線ベクトル,$\m{v}_3$を$C$の従法線ベクトルという.

$\m{v}_1$, $\m{v}_2$, $\m{v}_3$の正規直交性

これらのベクトル$\m{v}_1,\m{v}_2,\m{v}_3$については,次を当たり前にしておきましょう.

空間曲線$C:\m{r}=\m{r}(t)$の単位接ベクトル$\m{v}_1$,主法線ベクトル$\m{v}_2$,従法線ベクトル$\m{v}_3$は$\R^3$の正規直交基底をなす.

$\m{v}_1$, $\m{v}_2$の定義より,

\begin{align*}&\|\m{v}_1\|=\nor{\frac{\m{r}’}{\|\m{r}’\|}}=\frac{1}{\|\m{r}’\|}\|\m{r}’\|=1
\\&\|\m{v}_2\|=\nor{\frac{\m{v}’_1}{\|\m{v}’_1\|}}=\frac{1}{\|\m{v}’_1\|}\|\m{v}’_1\|=1\end{align*}

である.また,$\|\m{v}_1\|=1$だから,補題より

\begin{align*}\anb{\m{v}_1,\m{v}_2}=\anb{\m{v}_1,\frac{\m{v}’_1}{\|\m{v}’_1\|}}=\frac{1}{\|\m{v}’_1\|}\anb{\m{v}_1,\m{v}’_1}=0\end{align*}

である.さらに,外積の性質より$\m{v}_3=\m{v}_1\times \m{v}_2$は$\m{v}_1$, $\m{v}_2$と直交し,

\begin{align*}\|\m{v}_3\|=\abs{\|\m{v}_1\|\|\m{v}_2\|\sin{\bra{\pm\frac{\pi}{2}}}}=1\end{align*}

である.

フルネ-セレの公式

単位接ベクトル$\m{v}_1$,主法線ベクトル$\m{v}_2$,従法線ベクトル$\m{v}_3$が$\R^3$の(正規直交)基底をなすことから,これら3本のベクトル値関数の導関数$\m{v}’_1$, $\m{v}’_2$, $\m{v}’_3$を$\m{v}_1$, $\m{v}_2$, $\m{v}_3$の線形結合で一意に表すことができます.

この公式は

  • 1847年にジャン・フレデリック・フルネ(Jean Frédéric Frenet)
  • 1851年にジョセフ・アルフレッド・セレ(Joseph Alfred Serret)

によって独立に発見されましたのでフルネ-セレの公式といいます.

[フルネ-セレの公式]空間曲線$C:\m{r}=\m{r}(t)$は$\|\m{r}'(t)\|\equiv1$を満たすとする.このとき,ある非負値関数$\kappa:(\alpha,\beta)\to\R_{\ge0}$と実数値関数$\tau:(\alpha,\beta)\to\R$が存在して,任意の$t\in(\alpha,\beta)$に対して

\begin{align*}\begin{cases}\m{v}’_1(t)=\kappa(t)\m{v}_2(t)\\\m{v}’_2(t)=-\kappa(t)\m{v}_1(t)+\tau(t)\m{v}_3(t)\\\m{v}’_3(t)=-\tau(t)\m{v}_2(t)\end{cases}\end{align*}

が成り立つ.

$\|\m{r}'(t)\|\equiv1$を満たすとき,曲線上を点が進む速さが1ということなので,パラメータ$t$は弧長パラメータであるといいます.

フルネ-セレの公式は

\begin{align*}[\m{v}’_1(t),\m{v}’_2(t),\m{v}’_3(t)]=[\m{v}_1(t),\m{v}_2(t),\m{v}_3(t)]\bmat{0&\kappa(t)&0\\-\kappa(t)&0&\tau(t)\\0&-\tau(t)&0}\end{align*}

と表すこともできますね.

曲率・捩率の定義

上のフルネ-セレの公式の

  • 非負値関数$\kappa:(\alpha,\beta)\to[0,\infty)$を空間曲線$C:\m{r}=\m{r}(t)$の曲率
  • 実数値関数$\tau:(\alpha,\beta)\to\R$を空間曲線$C:\m{r}=\m{r}(t)$の捩率(れいりつ)

という.

第1式の両辺でノルムをとれば

\begin{align*}\|\m{v}’_1\|=|\kappa|\|\m{v}_2\|=\kappa\end{align*}

となるので,曲率は単位接ベクトルの導関数$\m{v}’_1$のノルムとして定義することもできます.





フルネ-セレの公式の証明

それでは,フルネ-セレの公式を証明しましょう.

[フルネ-セレの公式(再掲)]空間曲線$C:\m{r}=\m{r}(t)$は$\|\m{r}'(t)\|\equiv1$を満たすとする.このとき,ある非負値関数$\kappa:(\alpha,\beta)\to\R_{\ge0}$と実数値関数$\tau:(\alpha,\beta)\to\R$が存在して,任意の$t\in(\alpha,\beta)$に対して

\begin{align*}\begin{cases}\m{v}’_1(t)=\kappa(t)\m{v}_2(t)\\\m{v}’_2(t)=-\kappa(t)\m{v}_1(t)+\tau(t)\m{v}_3(t)\\\m{v}’_3(t)=-\tau(t)\m{v}_2(t)\end{cases}\end{align*}

が成り立つ.

任意の$t\in(\alpha,\beta)$に対して$\m{v}_1(t)$, $\m{v}_2(t)$, $\m{v}_3(t)$は$\R^3$の基底をなすことに注意する.

第1式$\m{v}’_1=\kappa\m{v}_2$の証明

$\m{v}_2$の定義の両辺に$\|\m{v}’_1\|$をかけて

\begin{align*}\m{v}_2=\frac{\m{v}’_1}{\|\m{v}’_1\|}\iff \m{v}’_1=\|\m{v}’_1\|\m{v}_2\end{align*}

が成り立つので,$\kappa=\|\m{v}’_1\|$とおけば$\kappa$は非負値で第1式が成り立つ.

第3式$\m{v}’_3=-\tau\m{v}_2$の証明

$\m{v}_1$, $\m{v}_2$, $\m{v}_3$は$\R^3$の基底をなすから,

\begin{align*}\m{v}’_3=a_{1}\m{v}_1+a_{2}\m{v}_2+a_{3}\m{v}_3\end{align*}

で$a_{i}:(\alpha,\beta)\to\R$($i=1,2,3$)を定めることができる.$a_{1}=a_{3}\equiv0$を示せばよい.

$\|\m{v}_3\|\equiv1$だから,補題より$\anb{\m{v}_3,\m{v}’_3}\equiv0$が成り立つので

\begin{align*}0\equiv\anb{\m{v}_3,a_{1}\m{v}_1+a_{2}\m{v}_2+a_{3}\m{v}_3}=a_{3}\end{align*}

が成り立つ.また,内積の微分公式と,上で示した第1式を併せて

\begin{align*}a_{1}&=\anb{a_{1}\m{v}_1+a_{2}\m{v}_2,\m{v}_1}=\anb{\m{v}’_3,\m{v}_1}
\\&=\anb{\m{v}_3,\m{v}_1}’-\anb{\m{v}_3,\m{v}’_1}
\\&=0′-\anb{\m{v}_3,\kappa\m{v}_2}=-\kappa\anb{\m{v}_3,\m{v}_2}=0\end{align*}

が成り立つ.よって,$\tau=-a_2$とおけば第3式が成り立つ.

第2式$\m{v}’_2=-\kappa\m{v}_1+\tau\m{v}_3$の証明

$\m{v}_1$, $\m{v}_2$, $\m{v}_3$は$\R^3$の基底をなすから,

\begin{align*}\m{v}’_2=b_{1}\m{v}_1+b_{2}\m{v}_2+b_{3}\m{v}_3\end{align*}

で$b_{i}:(\alpha,\beta)\to\R$($i=1,2,3$)を定めることができる.$b_{2}\equiv0$, $b_{1}=-\kappa$, $b_{3}=\tau$を示せばよい.

$\|\m{v}_2\|=1$だから,補題より$\anb{\m{v}_2,\m{v}’_2}\equiv0$が成り立つので

\begin{align*}0\equiv\anb{\m{v}_2,b_{1}\m{v}_1+b_{2}\m{v}_2+b_{3}\m{v}_3}=b_{2}\end{align*}

が成り立つ.また,内積の微分公式と,上で示した第1式を併せて

\begin{align*}\kappa&=\kappa\anb{\m{v}_2,\m{v}_2}=\anb{\m{v}_2,\kappa\m{v}_2}=\anb{\m{v}_2,\m{v}’_1}
\\&=\anb{\m{v}_2,\m{v}_1}’-\anb{\m{v}’_2,\m{v}_1}
\\&=0′-\anb{b_{1}\m{v}_1+b_{3}\m{v}_3,\m{v}_1}=-b_{1}\end{align*}

だから$b_{1}=-\kappa$が成り立つ.さらに,内積の微分公式と,上で示した第3式を併せて

\begin{align*}\tau&=\tau\anb{\m{v}_2,\m{v}_2}=\anb{\m{v}_2,\tau\m{v}_2}=-\anb{\m{v}_2,\m{v}’_3}
\\&=-\{\anb{\m{v}_2,\m{v}_3}’-\anb{\m{v}’_2,\m{v}_3}\}
\\&=0′-\anb{-\kappa\m{v}_1+b_{3}\m{v}_3,\m{v}_3}=b_{3}\end{align*}

だから$b_{3}=\tau$が成り立つ.よって,第2式が成り立つ.

管理人

プロフィール

山本やまもと 拓人たくと

元予備校講師.講師として駆け出しの頃から予備校の生徒アンケートで抜群の成績を残し,通常の8倍の報酬アップを提示されるなど頭角を表す.

飛び級・首席合格で大学院に入学しそのまま首席修了するなど数学の深い知識をもち,本質をふまえた分かりやすい授業に定評がある.

現在はオンライン家庭教師,社会人向け数学教室での講師としての教育活動とともに,京都大学で数学の研究も行っている.専門は非線形偏微分方程式論.大学数学系YouTuberとしても活動中.

趣味は数学,ピアノ,甘いもの食べ歩き.公式LINEを友達登録で【限定プレゼント】配布中.

Twitter・大学数学YouTube・公式LINEを見てみる

コメント