曲率・曲率半径の感覚的な意味と求め方

  • 曲率半径とは,曲線を「局所的に円の弧」とみなしたときの円の半径。

  • 曲率とは,曲率半径の逆数。

曲率・曲率半径について解説します。曲率半径を題材とした入試問題もときどき出題されます。

曲率について

  • 二階微分可能な曲線 y=f(a)y=f(a)(a,f(a))(a,f(a)) 付近で円に近似できます。その円を曲率円,半径を曲率半径と言います。曲率半径が大きいほどカーブはゆるいです。 曲率

  • 曲率は曲率半径の逆数です。曲線の(局所的な)曲がり具合を表します。曲率が大きいほどカーブは急です。

  • 直線の場合(曲がっていない場合)曲率は 00,曲率半径は \infty とみなせます。

  • このページでは,二階微分可能な曲線を考えます。y=xy=|x|x=0x=0 など,とがった点では円弧で近似できないので曲率は定義されません。

  • 曲率がどの点でも一定な曲線は円です。曲率の変化率が一定であるような曲線はクロソイド曲線と呼ばれるものです。→クロソイド曲線の性質とその証明

曲率半径を求める公式

曲率半径を求める公式1

y=f(x)y=f(x) の点 A(a,f(a))A(a,f(a)) における曲率半径は,

R=(1+f(a)2)32f(a)R=\dfrac{(1+f'(a)^2)^{\frac{3}{2}}}{|f''(a)|}

ただし,f(a)=0f''(a)=0 のときは曲率半径は \infty とみなす。

この公式は覚える必要は無いですが,導出の考え方は覚えておくとよいです。他にも導出方法はありますが,入試で題材にしやすい&曲率円の中心も求まる導出方法を解説します。

方針は比較的簡単で面白いですが,計算はけっこう大変です。

AA における曲率円の中心を,AA における法線と AA に近い点 BB における法線の交点,の極限」と解釈します。

公式の証明の概略
  1. まずは A(a,f(a))A(a,f(a)) における法線の方程式を求める(→ 法線ベクトルの3通りの求め方と応用):
    (yf(a))f(a)=ax(y-f(a))f'(a)=a-x 曲率の公式の証明

  2. 同様に B(b,f(b))B(b,f(b)) における法線の方程式も求まるので,二つの法線の交点の座標 C(xb,yb)C(x_b,\:y_b) を求める(詳細は省略)

  3. bbaa に近づけていくと (xb,yb)(x_b,y_b) は「AA 付近で f(x)f(x) を円弧とみなしたときの円」の中心となる。
    つまり,曲率円の中心は (limbaxb,limbayb)(\displaystyle\lim_{b\to a}x_b,\:\lim_{b\to a}y_b)
    実際に極限を計算すると,
    limbayb=1f(a)+f(a)+f(a)2f(a)\displaystyle\lim_{b\to a}y_b=\dfrac{1}{f''(a)}+f(a)+\dfrac{f'(a)^2}{f''(a)}
    limbaxb=af(a)f(a)f(a)3f(a)\displaystyle\lim_{b\to a}x_b=a-\dfrac{f'(a)}{f''(a)}-\dfrac{f'(a)^3}{f''(a)}

  4. 曲率円の中心と AA の距離が曲率半径である。
    R2=(limbaxba)2+(limbaybf(a))2=(1+f(a)2)3f(a)2R^2=(\displaystyle\lim_{b\to a}x_b-a)^2+(\lim_{b\to a}y_b-f(a))^2\\=\dfrac{(1+f'(a)^2)^3}{f''(a)^2}

注:なお,この記事では「関数を局所的に円弧で近似したときの円」が曲率円だと定義していますが,これは感覚的な表現で数学的に厳密ではありません。

納得出来ない人は,上記の1〜3で導出した円こそが曲率円の定義で,上記の4が曲率半径の定義だ! と理解してもOKです。

例題

具体的な計算例として,放物線の曲率半径を求めてみます。

例題

f(x)=x2f(x)=x^2(a,a2)(a,a^2) における曲率半径を求めよ。

解答

f(a)=2a,f(a)=2f'(a)=2a,\:f''(a)=2 より,曲率半径は

R=(1+4a2)322R=\dfrac{(1+4a^2)^{\frac{3}{2}}}{2}

つまり曲率半径は a2a^2 が大きいほど大きくなる。放物線は原点付近が最も急カーブで,離れるに連れてカーブが緩やかになる。

媒介変数表示で表された曲線の曲率半径

曲率半径を求める公式2

曲線 (x(t),y(t))(x(t),y(t)) 上の点 (x(t0),y(t0))(x(t_0),y(t_0)) における曲率半径は,

R=(x2+y2)32xyyxR=\dfrac{(x'^2+y'^2)^{\frac{3}{2}}}{|x'y''-y'x''|}

ただし,x=dx(t)dtt=t0x'=\left.\dfrac{dx(t)}{dt}\right|_{t=t_0}y=dy(t)dtt=t0y'=\left.\dfrac{dy(t)}{dt}\right|_{t=t_0}x=d2x(t)dt2t=t0x''=\left.\dfrac{d^2x(t)}{dt^2}\right|_{t=t_0}y=d2y(t)dt2t=t0y''=\left.\dfrac{d^2y(t)}{dt^2}\right|_{t=t_0} です。

y=f(x)y=f(x) 版の公式から導出します。(x0x'\neq 0 という条件がつきますが)計算が楽です。

導出

y=f(x)y=f(x) 版の公式:
R={1+(dydx)2}32d2ydx2R=\dfrac{\left\{1+\left(\frac{dy}{dx}\right)^2\right\}^{\frac{3}{2}}}{\left|\frac{d^2y}{dx^2}\right|}

に,dydx=dtdxdydt=yx\dfrac{dy}{dx}=\dfrac{dt}{dx}\dfrac{dy}{dt}=\dfrac{y'}{x'} および d2ydx2=ddx(yx)=dtdxddt(yx)=1xxyyxx2\dfrac{d^2y}{dx^2}=\dfrac{d}{dx}\left(\dfrac{y'}{x'}\right)=\dfrac{dt}{dx}\dfrac{d}{dt}\left(\dfrac{y'}{x'}\right)=\dfrac{1}{x'}\dfrac{x'y''-y'x''}{x'^2} を代入すると得る。

曲率と加速度ベクトル

重要な性質

曲率は,曲線上を速度1で進む物体の加速度ベクトルの大きさ(つまり,速度ベクトルを微分したものの大きさ)

大雑把な導出

時刻 ss における速度ベクトルを v(s)v(s) とする。

下図を見ると,Δs\Delta s が小さい正の数のとき,

v(s+Δs)v(s)ΔθΔsR|v(s+\Delta s)-v(s)|\fallingdotseq\Delta\theta \fallingdotseq\dfrac{\Delta s}{R}

曲率と加速度

よって,limΔs0v(s+Δs)v(s)Δs=1R\displaystyle\lim_{\Delta s\to 0}\dfrac{|v(s+\Delta s)-v(s)|}{\Delta s}=\dfrac{1}{R}

加速度の定義より,左辺は加速度ベクトルの大きさと一致する。

なお,上記の「重要な性質」を曲率の定義とすることも多いです。

上記の「重要な性質」から曲率半径を求める公式2を導出することもできます。合成関数の微分公式を使ってひたすら微分するだけです。

曲率半径を求める公式2の導出
  • 曲線の媒介変数表示を (x(t),y(t))(x(t),y(t)) とする
  • ss を弧長パラメタとする(つまり,曲線上を速度 11 で進む物体の時刻 ss における位置が (x(s),y(s))(x(s),y(s)) と表せる)
  • dsdt=(dxdt)2+(dydt)2=x2+y2\dfrac{ds}{dt}=\sqrt{\left(\dfrac{dx}{dt}\right)^2+\left(\dfrac{dy}{dt}\right)^2}=\sqrt{x'^2+y'^2} に注意する。

速度ベクトルの xx 成分は,dxds=dtdsdxdt=xx2+y2\dfrac{dx}{ds}=\dfrac{dt}{ds}\dfrac{dx}{dt}=\dfrac{x'}{\sqrt{x'^2+y'^2}}

加速度ベクトルの xx 成分は,
d2xds2=dtdsddt(xx2+y2)=1x2+y2xx2+y2x2xx+2yy2x2+y2x2+y2=xxx2+yyxx2+y2x2+y2=xy2xyy(x2+y2)2=y(yxxy)(x2+y2)2\dfrac{d^2x}{ds^2}=\dfrac{dt}{ds}\dfrac{d}{dt}\left(\dfrac{x'}{\sqrt{x'^2+y'^2}}\right)\\ =\dfrac{1}{\sqrt{x'^2+y'^2}}\dfrac{x''\sqrt{x'^2+y'^2}-x'\frac{2x'x''+2y'y''}{2\sqrt{x'^2+y'^2}}}{x'^2+y'^2}\\ =\dfrac{x''-\frac{x''x'^2+y''y'x'}{x'^2+y'^2}}{x'^2+y'^2}\\ =\dfrac{x''y'^2-x'y'y''}{(x'^2+y'^2)^2}\\ =\dfrac{y'(y'x''-x'y'')}{(x'^2+y'^2)^2}

対称性より yy 成分は,
d2yds2=x(xyyx)(x2+y2)2\dfrac{d^2y}{ds^2}=\dfrac{x'(x'y''-y'x'')}{(x'^2+y'^2)^2}

よって,重要な性質より

1R=(d2xds2)2+(d2yds2)2=xyyx(x2+y2)32\dfrac{1}{R}=\sqrt{\left(\dfrac{d^2x}{ds^2}\right)^2+\left(\dfrac{d^2y}{ds^2}\right)^2}\\ =\dfrac{|x'y''-y'x''|}{(x'^2+y'^2)^{\frac{3}{2}}}

線路は曲率が急激に変化しないようになっているらしいです。