1117 words
6 minutes
工科数学分析-4.2

Ex.1 环面的参数方程:#

考虑yOzyOz平面上的圆周C:(ya)2+z2=b2C:(y-a)^2+z^2=b^2,其中a>b>0a>b>0
改曲线绕zz轴一周形成环面.求曲面方程
z=bsinu,x=rcosv=(a+bcosu)cosv,y=rsinv=(a+bcosu)sinvz=b\sin u,x=r\cos v=(a+b\cos u)\cos v,y=r\sin v=(a+b\cos u)\sin v

Ex.2 莫比乌斯带#

矩形的一组对边反向粘合
是两个旋转的结合,一个转速是另一个的两倍
AB=2l,AA=2πa,a>lAB=2l,AA'=2\pi a,a>l
开始的时候细杆占据的位置 {(a,0,u)lu}\{(a,0,u)|l\ge|u|\}
细杆中点CC沿xOyxOy平面上的圆周x2+y2=a2x^2+y^2=a^2逆时针旋转
细杆CC自转,使LLzz轴所在的平面与xOyxOy平面垂直,且当CCθ\theta角时,细杆绕CCθ2\frac{\theta}{2}
z=ucosθ2z=ucos\frac{\theta}{2}
x=OPcosθ=(a+usinθ2)cosθx=OP'\cos\theta=(a+u\sin\frac{\theta}{2})\cos\theta
x=OPsinθ=(a+usinθ2)sinθx=OP'\sin\theta=(a+u\sin\frac{\theta}{2})\sin\theta

5.7 空间曲线的曲率与挠率#

5.7.1 FrenetFrenet标架#

设空间曲线Γ\Gamma的自然参数方程为r=r(s)r=r(s)
一般的参数方程r=r(t)r=r(t)
求导:
drds=r,d2rds2=r,drdt=r˙,d2rdt2=r¨\frac{dr}{ds}=r',\frac{d^2r}{ds^2}=r'',\frac{dr}{dt}=\dot{r},\frac{d^2r}{dt^2}=\ddot{r}

1.切线与法平面#

与曲线Γ:r=r(s)\Gamma:r=r(s)正向一致的单位切向量
T(s)=r(s)T(s)=r'(s)
dr=(dx,dy,dz),(ds)2=(dx)2+(dy)2+(dz)2=dr2dr=(dx,dy,dz),(ds)^2=(dx)^2+(dy)^2+(dz)^2=\|dr\|^2
drds=1\|\frac{dr}{ds}\|=1
于是Γ\Gamma在点r(S0)r(S_0)的切线方程为ρ=r(S0)=λr(S0)\rho=r(S_0)=\lambda r'(S_0)
法平面方程
(ρr(S0))r(S0)=0(\rho-r(S_0))r'(S_0)=0

2.主法线与次法线#

因为 r(S)2=r(S)r(S)\|r'(S)\|^2=r'(S)\cdot r'(S)
SS求导
0=2r(S)r(S)0=2r'(S)r''(S)
说明
r(S)r(S),r(S)×r(S)0r'(S)\perp r''(S),r'(S)\times r''(S)\neq0

N(S0)=r(S0)r(S0)N(S_0)=\frac{r''(S_0)}{\|r''(S_0)\|}
为曲线Γ\Gamma主法向量,对应的法线称为主法线.

B(S0)=T(S0)×N(S0)B(S_0)=T(S_0)\times N(S_0)
为曲线Γ\Gamma在点r(S0)r(S_0)次法向量,对应的法线称为次法线.
B(S)B(S)N(S)N(S)张成法平面

3.密切平面#

过曲线Γ\Gamma上的一点r(S0)r(S_0)的平面有无穷多个,哪一个与Γ\Gamma最贴近?
Γ\Gamma是一个空间曲线时,考虑r(S0)r(S_0)临近的点,r(S0+ΔS)r(S_0+\Delta S)到平面的距离,设Σ\Sigma是过点r(S0)r(S_0)平面的单位法向量
r(S0+ΔS)r(S_0+\Delta S)到平面Σ\Sigma的距离
d=λ=n(r(S0+ΔS)r(S0))d=\lambda=n\cdot (r(S_0+\Delta S)-r(S_0))
r(S0+ΔS)r(S0)=λn+nr(S_0+\Delta S)-r(S_0)=\lambda n+n^{\perp}
TaylorTaylor公式
r(S0+ΔS)r(S0)=r(S0)ΔS+12r(S0)(ΔS)2+o((ΔS)2)r(S_0+\Delta S)-r(S_0)=r'(S_0)\Delta S+\frac{1}{2}r''(S_0)(\Delta S)^2+o((\Delta S)^2)
于是
d=n(r(S0)ΔS+12r(S0)(ΔS)2+o((ΔS)2))d=n\cdot(r'(S_0)\Delta S+\frac{1}{2}r''(S_0)(\Delta S)^2+o((\Delta S)^2))

(1)(1)Σ\Sigma不含切线,则nr(S0)0n\cdot r'(S_0)\neq0,此时r(S0+ΔS)r(S0)r(S_0+\Delta S) \approx r(S_0)
偏差dd是一阶无穷小
(2)(2)Σ\Sigma含切线但不含主法线,则
nr(S0)=0,nr(S0)0n\cdot r'(S_0)=0,n\cdot r''(S_0)\neq0
偏差dd是二阶无穷小
(3)(3)Σ\Sigma含切线与主法线,则
nr(S0)=0,nr(S0)=0n\cdot r'(S_0)=0,n\cdot r''(S_0)=0
偏差d=o((Δ)2),r(S0+Δ)r(S0)d=o((\Delta)^2),r(S_0+\Delta)\approx r(S_0)最佳逼近
因此,含曲线Γ\Gamma在点r(S0)r(S_0)处的切线和主法线的平面称为密切平面,方程为
B(S0)(ρr(S0))=0B(S_0)\cdot (\rho-r(S_0))=0
注:密切平面是平面到曲线上的点的距离为三阶无穷小的平面

4.从切平面#

切向量T(S0)T(S_0)与次法向量B(S0)B(S_0)确定的平面为从切平面
三个单位向量T,N,BT,N,B互相正交,它们构成以r(S0)r(S_0)为原点的右手直角坐标系,称为FrenetFrenet标架
Γ\Gamma的方程由一般的参数方程r=r(t)r=r(t)给出
r˙(t0)0,r¨(t0)0\dot{r}(t_0)\neq0,\ddot{r}(t_0)\neq0时,单位切向量,主法向量,次法向量分别为
T(t0)=r˙(t0)r˙(t0),N(t0)=r¨(t0)r¨(t0),B(t0)=T(t0)×N(t0)=r˙(t0)×r¨(t0)r˙(t0)×r¨(t0)T(t_0)=\frac{\dot{r}(t_0)}{\|\dot{r}(t_0)\|},N(t_0)=\frac{\ddot{r}(t_0)}{\|\ddot{r}(t_0)\|},B(t_0)=T(t_0)\times N(t_0)=\frac{\dot{r}(t_0)\times\ddot{r}(t_0)}{\|\dot{r}(t_0)\times\ddot{r}(t_0)\|}

例 7.1#

求曲线r=(acost,asint,kt)r=(a\cos t,a\sin t,kt)FrenetFrenet标架,密切平面,从切平面的方程
解:
r˙(t0)=(asint,acost,k)\dot{r}(t_0)=(-a\sin t,a\cos t,k)
r¨(t0)=(acost,asint,0)\ddot{r}(t_0)=(-a\cos t,-a\sin t,0)
r˙(t0)×r¨(t0)=(aksint,akcost,a2)\dot{r}(t_0)\times\ddot{r}(t_0)=(ak\sin t,-ak\cos t,a^2)
r˙=a2+k2,r¨=a,r˙×r¨=aa2+k2\|\dot{r}\|=\sqrt{a^2+k^2},\|\ddot{r}\|=a,\|\dot{r}\times\ddot{r}\|=a\sqrt{a^2+k^2}
于是
T=1a2+k2(asint,acost,k)T=\frac{1}{\sqrt{a^2+k^2}}(-a\sin t,a\cos t,k)
N=(sint,cost,0)N=(-\sin t,\cos t,0)
T=1a2+k2(sint,cost,a)T=\frac{1}{\sqrt{a^2+k^2}}(-\sin t,\cos t,a)
密切平面的方程为B(t)(ρr(t)=0)B(t)\cdot(\rho-r(t)=0),即
ksint(xacost)kcost(yasint)+a(zkt)=0k\sin t(x-a\cos t)-k\cos t(y-a\sin t)+a(z-kt)=0
(ksint)x(kcost)y+az=akt(k\sin t)x-(k\cos t)y+az=akt
从切平面方程为N(t)(ρr(t))=0N(t)\cdot (\rho-r(t))=0
cost(xacost)+sint(yasint)=0\cos t(x-a\cos t)+\sin t(y-a\sin t)=0
xcost+ysint=ax\cos t+y\sin t=a

5.7.2 曲率#

1.定义与计算#

S0S_0S0+ΔSS_0+\Delta S的弧段的平均弯曲程度,称为平均曲率,记为κ=ΔθΔS\overline{\kappa}=|\frac{\Delta \theta}{\Delta S}|

定义 7.1#

假设空间光滑的曲线Γ\Gamma的方程为r=r(s),(αSβ)r=r(s),(\alpha\le S\le\beta)
切向量r(S),r(S+ΔS)r'(S),r'(S+\Delta S)的夹角为Δθ\Delta \theta若极限κ=limΔS0ΔθΔS\kappa=\lim_{\Delta S\to 0}|\frac{\Delta \theta}{\Delta S}|存在
则称之为Γ\Gamma在点r(S)r(S)的曲率
ΔT=r(S+ΔS)r(S)\Delta T=r'(S+\Delta S)-r'(S)
因为
ΔθΔS=ΔθΔTΔTΔS=ΔTΔS(ΔTΔθ)1|\frac{\Delta \theta}{\Delta S}|=\frac{\Delta\theta}{\|\Delta T\|}\frac{\|\Delta T\|}{\Delta S}=\|\frac{\Delta T}{\Delta S}\|(\|\frac{\Delta T}{\Delta\theta}\|)^{-1}
所以
κ=limΔS0ΔTΔS(ΔTΔθ)1\kappa=\lim_{\Delta S\to0}\|\frac{\Delta T}{\Delta S}\|(\|\frac{\Delta T}{\Delta\theta}\|)^{-1}
又因为

ΔT2=r(S+ΔS)r(S)2=22cosΔθ=4sin2Δθ2\begin{align} \|\Delta T\|^2&=\|r'(S+\Delta S)-r'(S)\|^2\\ &=2-2\cos \Delta\theta\\ &=4\sin^2\frac{\Delta \theta}{2} \end{align}

所以
ΔTΔθ2=4sin2Δθ2(Δθ)2=1\|\frac{\Delta T}{\Delta \theta}\|^2=\frac{4\sin^2\frac{\Delta \theta}{2}}{(\Delta \theta)^2}=1
因此
κ=limΔS0ΔTΔS=limΔS0ΔTΔS=limΔS0r(S+ΔS)r(S)ΔS=r(S)=T(S)\kappa=\lim_{\Delta S\to0}\|\frac{\Delta T}{\Delta S}\|=\|\lim_{\Delta S\to 0}\frac{\Delta T}{\Delta S}\|=\|\lim_{\Delta S\to0}\frac{r'(S+\Delta S)-r'(S)}{\Delta S}\|=\|r''(S)\|=\|T'(S)\|

推论#

设空间曲线Γ\Gamma由一般参数方程r=r(t)r=r(t)给出,r(t)r(t)二阶可导且r˙(t)0\dot{r}(t)\neq0,则Γ\Gammar(t)r(t)的曲率为

证明:
r(s)=drds=drdtdtds,r(s)=r˙(t)dtdsr'(s)=\frac{dr}{ds}=\frac{dr}{dt}\frac{dt}{ds},r'(s)=\dot{r}(t)\frac{dt}{ds}
r(s)=d2rds2=d2rdt2(dtds)2+drdtd2tds2r''(s)=\frac{d^2r}{ds^2}=\frac{d^2r}{dt^2}(\frac{dt}{ds})^2+\frac{dr}{dt}\frac{d^2t}{ds^2}
因为r(s)=1\|r'(s)\|=1,所以
1=r˙dtds,r(s)r(s)=01=\|\dot{r}\||\frac{dt}{ds}|,r'(s)\cdot r''(s)=0
由向量积的定义
r(s)×r(s)=r(s)r(s)sin(r(s),r(s))=r(s)\|r'(s)\times r''(s)\|=\|r'(s)\|\|r''(s)\|\sin(\langle r'(s),r''(s)\rangle)=\|r''(s)\|
于是

κ(t)=r(s)=r(s)×r(S)=r(t)×(r¨(t)(dtds)2+r(s)×r˙(t)d2tds2)=d˙(t)dtds×r¨(t)(dtds)2=r˙(t)×r¨(t)dtds3=r˙(t)×r¨(t)r˙(t)3\begin{align} \kappa(t)&=\|r''(s)\|=\|r'(s)\times r''(S)\|\\ &=\|r'(t)\times (\ddot{r}(t)(\frac{dt}{ds})^2+r'(s)\times \dot{r}(t)\frac{d^2t}{ds^2})\|\\ &=\|\dot{d}(t)\frac{dt}{ds}\times \ddot{r}(t)(\frac{dt}{ds})^2\|\\ &=\|\dot{r}(t)\times \ddot{r}(t)\||\frac{dt}{ds}|^3\\ &=\frac{\|\dot{r}(t)\times \ddot{r}(t)\|}{\|\dot{r}(t)\|^3} \end{align}

对平面曲线r=(x(t),y(t),0)r=(x(t),y(t),0)
r˙=(x˙,y˙,0)\dot{r}=(\dot{x},\dot{y},0)
r¨=(x¨,y¨,0)\ddot{r}=(\ddot{x},\ddot{y},0)
r˙×r¨=(x˙y¨x¨y˙,0,0)\dot{r}\times \ddot{r}=(\dot{x}\ddot{y}-\ddot{x}\dot{y},0,0)
κ=x˙y¨x¨y˙(x˙2+y˙2)32\kappa=\frac{|\dot{x}\ddot{y}-\ddot{x}\dot{y}|}{(\dot{x}^2+\dot{y}^2)^{\frac{3}{2}}}
r=(x,y(x),0)r=(x,y(x),0)
κ=y(x)(1+y(x)2)32\kappa=\frac{|y''(x)|}{(1+y'(x)^2)^{\frac{3}{2}}}

例 7.2#

证明:曲线LL为直线\Leftrightarrow它的曲率处处为零
proof
κ=r(s)=0r(s)=0r(s)=r0,r(s)r(s0)=(ss0)r0\kappa=\|r''(s)\|=0\Leftrightarrow r''(s)=0\Leftrightarrow r'(s)=r_0,\Leftrightarrow r(s)-r(s_0)=(s-s_0)r_0

例 7.3#

求半径为RR的圆的曲率
ΔS=RΔθ\Delta S=R\Delta \theta
κ=limΔS0ΔθΔS1R\kappa=\lim_{\Delta S\to 0}|\frac{\Delta \theta}{\Delta S}|\frac{1}{R}

例 7.4#

r=(acost,asint,kt)r=(a\cos t,a\sin t,kt)的曲率
κ(t)=r˙(t)×r¨(t)r˙(t)3=aa2+k2(a2+k2)32=aa2+k2\kappa(t)=\frac{\|\dot{r}(t)\times \ddot{r}(t)\|}{\|\dot{r}(t)\|^3}=\frac{a\sqrt{a^2+k^2}}{(\sqrt{a^2+k^2})^{\frac{3}{2}}}=\frac{a}{a^2+k^2}

2.曲率半径与曲率圆#

Γ\Gamma为平面曲线,过点P0ΓP_0\in \Gamma在凹向侧的法线,取点QQ,使PQ=1κ|PQ|=\frac{1}{\kappa},以QQ为圆心,PQ|PQ|为半径的圆称为Γ\Gamma在点PP的曲率圆
Γ\Gamma为空间曲线,则在点PPΓ\Gamma最贴近的平面是密切平面,记为Π\Pi
在主法线上取点QQ,使PQ\vec{PQ}NN同向,且PQ=1κ\|PQ\|=\frac{1}{\kappa},以QQ为圆心,1κ\frac{1}{\kappa}为半径的圆,
称为Γ\GammaPP曲率圆密切圆
若以弧长为参数,则曲率中心的向径rQr_Q与曲率半径RR,分别为
rQrP=R(s)N(s),R=1κr_Q-r_P=R(s)N(s),R=\frac{1}{\kappa}