5.7 空间曲线的曲率与挠率

本节将应用多元函数微分学的知识进一步对空间曲线的形态作更深入的研究。为此，先介绍 Frenet 标架，然后研究空间曲线的曲率与挠率。

5.7 空间曲线的曲率与挠率

7.1 Frenet 标架

在曲线的研究中，用曲线上建立活动的坐标系来取代原有的固定的直角坐标系会带来许多方便。例如，船在大海中航行，对于时刻 $t$ 在海平面上或空中所发现的目标，往往以此时刻船的位置为坐标原点所建立的某一坐标系来测量目标的位置。当 $t$ 变化时，这一坐标系将作为一刚体随船的航线而变化。

回顾空间直角坐标系，我们是用三个相互正交的平面的交线来构造坐标系的。要在曲线上建立活动的坐标架，自然也希望寻找三个由曲线特征所确定的相互正交的平面，用它们的交线来构造坐标架。以下设曲线 $\Gamma$ 的自然参数方程为 $\mathbf{r}=\mathbf{r}(s)$，而 $\Gamma$ 的一般参数方程 $\mathbf{r}=\mathbf{r}(t)$ 表示。为了区分它们的导数，我们用 ${\mathbf{r}}^{\prime },{\mathbf{r}}^{\prime \prime },\dots$ 分别表示 $\frac{\mathrm{d}\mathbf{r}}{\mathrm{d}s},\frac{{\mathrm{d}}^2\mathbf{r}}{\mathrm{d}{s}^2},\dots$，用 $\dot{\mathbf{r}},\ddot{\mathbf{r}},\dots$ 分别表示 $\frac{\mathrm{d}\mathbf{r}}{\mathrm{d}t},\frac{{\mathrm{d}}^2\mathbf{r}}{\mathrm{d}{t}^2},\dots$，并且假定所出现的函数 $\mathbf{r}$ 的导数都是连续的。

1. 法平面与切线

设 ${\mathbf{r}}^{\prime }(s)\ne 0$。在上一节中已经指出，${\mathbf{r}}^{\prime }(s_0)$ 是 $\Gamma$ 在点 $\mathbf{r}(s_0)$ 处其正向与 $\Gamma$ 正向一致的单位切向量，记作 $\mathbf{T}(s_0)$。即

T(s₀) = r'(s₀), (7.1)$$\mathbf{T}(s_0)={\mathbf{r}}^{\prime }(s_0),\text{\hspace{1em}(7.1)}$$

于是 $\Gamma$ 在点 $\mathbf{r}(s_0)$ 的法平面方程为

(ρ - r(s₀)) · r'(s₀) = 0,$$(\mathbf{\rho }-\mathbf{r}(s_0))\cdot {\mathbf{r}}^{\prime }(s_0)=0,$$

其中 $\mathbf{\rho }$ 是法平面上动点的径向。$\Gamma$ 在点 $\mathbf{r}(s_0)$ 的切线方程为

ρ = r(s₀) + λ r'(s₀),$$\mathbf{\rho }=\mathbf{r}(s_0)+\lambda {\mathbf{r}}^{\prime }(s_0),$$

其中 $\mathbf{\rho }$ 是切线上动点的径向，$\lambda \in \mathbb{R}$ 为参数。

2. 密切平面与次法线

通过 $\Gamma$ 上点 $\mathbf{r}(s_0)$ 的切线的平面称为 $\Gamma$ 在点 $\mathbf{r}(s_0)$ 的切平面，显然，这种切平面有无穷多个，其中哪一个与曲线 $\Gamma$ 最贴近呢？如果 $\Gamma$ 是平面曲线，那么显然 $\Gamma$ 所在的平面就是与 $\Gamma$ 最贴近的切平面。当 $\Gamma$ 是空间曲线时，与它最贴近的切平面可以用如下的方法得到：将 $\Gamma$ 在点 $\mathbf{r}(s_0)$ 的切线与 $\Gamma$ 上 $\mathbf{r}(s_0)$ 邻近的点 $\mathbf{r}(s_0+\Delta s)$ 所确定的平面记作 ${\pi }^{\prime }$，当 $\Delta s\to 0$ 时点 $\mathbf{r}(s_0+\Delta s)$ 将沿 $\Gamma$ 趋于点 $\mathbf{r}(s_0)$，如果此时 ${\pi }^{\prime }$ 有极限位置 $\pi$，那么就认为 $\pi$ 是与曲线 $\Gamma$ 最贴近的切平面，称之为曲线 $\Gamma$ 在点 $\mathbf{r}(s_0)$ 处的密切平面；而把密切平面在点 $\mathbf{r}(s_0)$ 的法线，称为曲线 $\Gamma$ 在点 $\mathbf{r}(s_0)$ 的次法线。下面来求它们的方程。

由平面 ${\pi }^{\prime }$ 的构造可知，向量 ${\mathbf{r}}^{\prime }(s_0)$ 与 $\mathbf{r}(s_0+\Delta s)-\mathbf{r}(s_0)$ 均在 ${\pi }^{\prime }$ 上，从而 ${\pi }^{\prime }$ 的法向量可取为

r'(s₀) × (r(s₀ + Δ s) - r(s₀)).$${\mathbf{r}}^{\prime }(s_0)\times (\mathbf{r}(s_0+\Delta s)-\mathbf{r}(s_0)).$$

对向量 $\mathbf{r}(s_0+\Delta s)-\mathbf{r}(s_0)$ 的各分量分别应用一元函数的 Taylor 公式，可得

r(s₀ + Δ s) - r(s₀) = r'(s₀) Δ s + (1)/(2!) (r''(s₀) + ε) Δ s²,$$\mathbf{r}(s_0+\Delta s)-\mathbf{r}(s_0)={\mathbf{r}}^{\prime }(s_0)\Delta s+\frac{1}{2!}({\mathbf{r}}^{\prime \prime }(s_0)+\mathbf{\epsilon })\Delta s^2,$$

其中 $\Delta s\to 0$ 时 $\mathbf{\epsilon }\to 0$。注意到 ${\mathbf{r}}^{\prime }(s_0)\times {\mathbf{r}}^{\prime }(s_0)=0$，从而

r'(s₀) × (r(s₀ + Δ s) - r(s₀)) = (Δ s²)/(2) r'(s₀) × (r''(s₀) + ε),$${\mathbf{r}}^{\prime }(s_0)\times (\mathbf{r}(s_0+\Delta s)-\mathbf{r}(s_0))=\frac{\Delta s^2}{2}{\mathbf{r}}^{\prime }(s_0)\times ({\mathbf{r}}^{\prime \prime }(s_0)+\mathbf{\epsilon }),$$

所以，${\mathbf{r}}^{\prime }(s_0)\times ({\mathbf{r}}^{\prime \prime }(s_0)+\mathbf{\epsilon })$ 也是 ${\pi }^{\prime }$ 的法向量。令 $\Delta s\to 0$，得到这个向量积的极限为 ${\mathbf{r}}^{\prime }(s_0)\times {\mathbf{r}}^{\prime \prime }(s_0)$。因此，当 ${\mathbf{r}}^{\prime }(s_0)\times {\mathbf{r}}^{\prime \prime }(s_0)\ne 0$ 时，可得 $\Gamma$ 在点 $\mathbf{r}(s_0)$ 次法线的方向向量（简称次法向量）为

r'(s₀) × r''(s₀).$${\mathbf{r}}^{\prime }(s_0)\times {\mathbf{r}}^{\prime \prime }(s_0).$$

由于 $\Vert {\mathbf{r}}^{\prime }(s)\Vert =1$，即向量 ${\mathbf{r}}^{\prime }(s)$ 的模不变，由本章习题5.5 ($A$) 第 4 题可知

r'(s₀) perp r''(s₀), (7.2)$${\mathbf{r}}^{\prime }(s_0)\perp {\mathbf{r}}^{\prime \prime }(s_0),\text{\hspace{1em}(7.2)}$$

于是由向量积的定义可知，$\Gamma$ 在点 $\mathbf{r}(s_0)$ 处的单位次法向量（记为 $\mathbf{B}(s_0)$）为

B(s₀) = (r'(s₀) × r''(s₀))/(|r'(s₀) × r''(s₀)|) = (r'(s₀) × r''(s₀))/(|r''(s₀)|) (7.3)$$\mathbf{B}(s_0)=\frac{{\mathbf{r}}^{\prime }(s_0)\times {\mathbf{r}}^{\prime \prime }(s_0)}{\Vert {\mathbf{r}}^{\prime }(s_0)\times {\mathbf{r}}^{\prime \prime }(s_0)\Vert }=\frac{{\mathbf{r}}^{\prime }(s_0)\times {\mathbf{r}}^{\prime \prime }(s_0)}{\Vert {\mathbf{r}}^{\prime \prime }(s_0)\Vert }\text{\hspace{1em}(7.3)}$$

所以，$\Gamma$ 在点 $\mathbf{r}(s_0)$ 的密切平面方程为

B(s₀) · (ρ - r(s₀)) = 0,$$\mathbf{B}(s_0)\cdot (\mathbf{\rho }-\mathbf{r}(s_0))=0,$$

或

(r'(s₀) × r''(s₀)) · (ρ - r(s₀)) = 0.$$({\mathbf{r}}^{\prime }(s_0)\times {\mathbf{r}}^{\prime \prime }(s_0))\cdot (\mathbf{\rho }-\mathbf{r}(s_0))=0.$$

$\Gamma$ 在 $\mathbf{r}(s_0)$ 处的次法线方程为

ρ = r(s₀) + λ B(s₀),$$\mathbf{\rho }=\mathbf{r}(s_0)+\lambda \mathbf{B}(s_0),$$

或

ρ = r(s₀) + λ (r'(s₀) × r''(s₀)).$$\mathbf{\rho }=\mathbf{r}(s_0)+\lambda ({\mathbf{r}}^{\prime }(s_0)\times {\mathbf{r}}^{\prime \prime }(s_0)).$$

3. 从切平面与主法线

由 $\Gamma$ 在点 $\mathbf{r}(s_0)$ 即点 $P$ 的切向量 $\mathbf{T}(s_0)={\mathbf{r}}^{\prime }(s_0)$ 与次法向量 $\mathbf{B}(s_0)$ 所确定的平面，称为曲线 $\Gamma$ 在点 $\mathbf{r}(s_0)$ 的从切平面（图 5.34），它显然既垂直于法平面也垂直于密切平面。从切平面在点 $\mathbf{r}(s_0)$ 的法线称为曲线 $\Gamma$ 在点 $\mathbf{r}(s_0)$ 的主法线，并取

N(s₀) = B(s₀) × T(s₀) (7.4)$$\mathbf{N}(s_0)=\mathbf{B}(s_0)\times \mathbf{T}(s_0)\text{\hspace{1em}(7.4)}$$

作为 $\Gamma$ 在点 $\mathbf{r}(s_0)$ 的主法线的单位方向向量（简称单位主法向量）。于是 $\Gamma$ 过点 $\mathbf{r}(s_0)$ 的从切平面与主法线的方程分别为

(ρ - r(s₀)) · N(s₀) = 0, ρ = r(s₀) + λ N(s₀).$$(\mathbf{\rho }-\mathbf{r}(s_0))\cdot \mathbf{N}(s_0)=0,\mathbf{\rho }=\mathbf{r}(s_0)+\lambda \mathbf{N}(s_0).$$

主法向量的表达式 (7.4) 还可以简化。把表达式 (7.3) 与 (7.1) 代入 (7.4) 得

N(s₀) = (1)/(|r''(s₀)|) (r'(s₀) × r''(s₀)) × r'(s₀),$$\mathbf{N}(s_0)=\frac{1}{\Vert {\mathbf{r}}^{\prime \prime }(s_0)\Vert }({\mathbf{r}}^{\prime }(s_0)\times {\mathbf{r}}^{\prime \prime }(s_0))\times {\mathbf{r}}^{\prime }(s_0),$$

再由 (7.2) 式可知，$\mathbf{N}(s_0)$ 与 ${\mathbf{r}}^{\prime \prime }(s_0)$ 平行且同向。又由于 $\Vert \mathbf{N}(s_0)\Vert =1$，所以

N(s₀) = (r''(s₀))/(|r''(s₀)|). (7.5)$$\mathbf{N}(s_0)=\frac{{\mathbf{r}}^{\prime \prime }(s_0)}{\Vert {\mathbf{r}}^{\prime \prime }(s_0)\Vert }.\text{\hspace{1em}(7.5)}$$

容易看出，曲线 $\Gamma$ 在点 $\mathbf{r}(s_0)$ 的从切平面、密切平面和法平面相互正交，它们的三条交线分别是 $\Gamma$ 的切线、主法线和次法线（图 5.34）。而相互正交的三个单位向量 $\mathbf{T},\mathbf{N},\mathbf{B}$ 构成以点 $P$ 为原点的空间（右手）直角坐标系，称为曲线 $\Gamma$ 在点 $P=\mathbf{r}(s)$ 处的 Frenet 标架。当点 $P$ 沿曲线 $\Gamma$ 移动时，这个标架贴附在 $\Gamma$ 上作为一个刚体随之而移动，因此又称其为活动标架。

以上，标架向量 $\mathbf{T},\mathbf{N},\mathbf{B}$ 是在曲线 $\Gamma$ 的方程由自然参数方程 $\mathbf{r}=\mathbf{r}(s)$ 给出的情形下推导出的。如果 $\Gamma$ 的方程由一般参数方程 $\mathbf{r}=\mathbf{r}(t)$ 给出，则可以证明，当 $\dot{\mathbf{r}}(t_0)\ne 0,\ddot{\mathbf{r}}(t_0)\ne 0$ 时，$\Gamma$ 上点 $\mathbf{r}(t_0)$ 处的单位切向量 $\mathbf{T}(t_0)$ 及单位次法向量 $\mathbf{B}(t_0)$ 的计算公式分别为

T(t₀) = fracdotr(t₀)|dotr(t₀)| 与 B(t₀) = fracdotr(t₀) × ddotr(t₀)|dotr(t₀) × ddotr(t₀)| (7.6)$$\mathbf{T}(t_0)=\frac{\dot{\mathbf{r}}(t_0)}{\Vert \dot{\mathbf{r}}(t_0)\Vert }\text{与}\mathbf{B}(t_0)=\frac{\dot{\mathbf{r}}(t_0)\times \ddot{\mathbf{r}}(t_0)}{\Vert \dot{\mathbf{r}}(t_0)\times \ddot{\mathbf{r}}(t_0)\Vert }\text{\hspace{1em}(7.6)}$$

而按照 (7.6) 式算出 $\mathbf{T}(t_0)$ 及 $\mathbf{B}(t_0)$ 时，则 $\mathbf{r}(t_0)$ 处的单位主法向量 $\mathbf{N}(t_0)$ 便可由公式 $\mathbf{N}(t_0)=\mathbf{B}(t_0)\times \mathbf{T}(t_0)$ 算出。上述向量 $\mathbf{T}(t_0),\mathbf{N}(t_0),\mathbf{B}(t_0)$ 即构成 $\Gamma$ 上 $\mathbf{r}(t_0)$ 处的 Frenet 标架向量。

例 7.1 螺旋线的 Frenet 标架与相关平面方程

求螺旋线 $\mathbf{r}=(a\cos t,a\sin t,kt)$ 的 Frenet 标架、密切平面以及从切平面的方程。

解

由于 $\dot{\mathbf{r}}=(-\text{a}\sin t,\text{a}\cos t,k)$, $\ddot{\mathbf{r}}=(-\text{a}\cos t,-\text{a}\sin t,0)$, $\dot{\mathbf{r}}\times \ddot{\mathbf{r}}=\text{a}(k\sin t,-k\cos t,\text{a})$,

所以

|dotr| = √(a² + k²), |dotr × ddotr| = a√(a² + k²).$$\Vert \dot{\mathbf{r}}\Vert =\sqrt{{\text{a}}^2+k^2},\Vert \dot{\mathbf{r}}\times \ddot{\mathbf{r}}\Vert =\text{a}\sqrt{{\text{a}}^2+k^2}.$$

于是由公式 (7.6) 即得

T = fracdotr|dotr| = (1)/(√(a² + k²))(-asin t, acos t, k),$$\mathbf{T}=\frac{\dot{\mathbf{r}}}{\Vert \dot{\mathbf{r}}\Vert }=\frac{1}{\sqrt{{\text{a}}^2+k^2}}(-\text{a}\sin t,\text{a}\cos t,k),$$B = fracdotr × ddotr|dotr × ddotr| = (1)/(√(a² + k²))(ksin t, -kcos t, a),$$\mathbf{B}=\frac{\dot{\mathbf{r}}\times \ddot{\mathbf{r}}}{\Vert \dot{\mathbf{r}}\times \ddot{\mathbf{r}}\Vert }=\frac{1}{\sqrt{{\text{a}}^2+k^2}}(k\sin t,-k\cos t,\text{a}),$$

从而得

N = B × T = (-cos t, -sin t, 0).$$\mathbf{N}=\mathbf{B}\times \mathbf{T}=(-\cos t,-\sin t,0).$$

上述三向量 $\mathbf{T},\mathbf{N},\mathbf{B}$ 所构成的坐标系即为所求 Frenet 标架。

密切平面的方程为

B(t) · (ρ - r(t)) = 0,$$\mathbf{B}(t)\cdot (\mathbf{\rho }-\mathbf{r}(t))=0,$$

即

ksin t(x - acos t) - kcos t(y - asin t) + a(z - kt) = 0.$$k\sin t(x-\text{a}\cos t)-k\cos t(y-\text{a}\sin t)+\text{a}(z-kt)=0.$$

从切平面的方程为

N(t) · (ρ - r(t)) = 0,$$\mathbf{N}(t)\cdot (\mathbf{\rho }-\mathbf{r}(t))=0,$$

即

-cos t(x - acos t) - sin t(y - asin t) = 0.$$-\cos t(x-\text{a}\cos t)-\sin t(y-\text{a}\sin t)=0.$$

7.2 曲率

1. 曲率的定义与计算

通俗地讲，所谓曲率就是曲线上各点的弯曲程度。在工程技术、生产实践和自然科学中，不少问题与曲线的弯曲程度有关。例如，在铁路转弯处，铁道的外轨要比内轨高，而且转弯愈急，“弯曲”愈大处，内、外轨的高度差也就愈大，这是由于做转动的列车在“弯曲”较大处所受的向心力愈大的缘故。再如，车床的主轴由于自重总会产生弯曲变形，如果弯曲过大，就会影响车床的精度和正常运行。为了从数量上刻画曲线的弯曲程度，需要引入一个度量概念——曲率。

如果所考察的平面曲线 $\Gamma$ 是半径为 $R$ 的圆周（图 5.35），曲线上各点的弯曲程度是一样的。这时，我们可以用单位弧段上，曲线切线所转过的角度 $\frac{\Delta \theta }{\Delta s}$ 来度量曲线的弯曲程度，其中 $\Delta \theta$ 是在弧段 $\Delta s$ 上切线所转过的角度。显然对圆周来说，

(Δ θ)/(Δ s) = C (常数)$$\frac{\Delta \theta }{\Delta s}=C\text{ (常数)}$$

即，曲线上切线的转角 $\theta$ 随弧长 $s$ 均匀变化。所以，我们就用 $\frac{\Delta \theta }{\Delta s}$ 来作为圆周弯曲程度的度量，称为圆周的曲率。

对于一般的曲线来说，各点处的弯曲程度不尽相同（图 5.36），换句话说，曲线上切线的转角 $\theta$ 随弧长 $s$ 的变化是非均匀的。这时 $\frac{\Delta \theta }{\Delta s}$ 并非常数，它只能表示从 $s_0$ 到 $s_0+\Delta s$ 弧段上的平均弯曲程度，称为平均曲率。为了精确地刻画曲线在一点的弯曲程度，需要令 $\Delta s\to 0$，对平均曲率取极限。容易看出，这就要利用导数的概念。为此，我们引入以下定义：

定义 7.1 (曲率)

设空间光滑曲线 $\Gamma$ 的方程为 $\mathbf{r}=\mathbf{r}(s)(a⩽s⩽b)$，其中 $s$ 为自然参数，$\Gamma$ 上点 $M$ 对应的参数为 $s$，$\Gamma$ 上点 $M$ 附近点 $N$ 对应的参数为 $s+\Delta s$，$\Gamma$ 在点 $M$ 处的切线向量 ${\mathbf{r}}^{\prime }(s)$ 与点 $N$ 处的切线向量 ${\mathbf{r}}^{\prime }(s+\Delta s)$ 的夹角为 $\Delta \theta$，称极限值

lim_(Δ s → 0) | (Δ θ)/(Δ s) |$$\underset{\Delta s\to 0}{\lim }\left|\frac{\Delta \theta }{\Delta s}\right|$$

为 $\Gamma$ 在点 $M$ 处的曲率，记为 $\kappa$ (读作 kappa)，即

κ = lim_(Δ s → 0) | (Δ θ)/(Δ s) | (7.1)$$\kappa =\underset{\Delta s\to 0}{\lim }\left|\frac{\Delta \theta }{\Delta s}\right|\text{\hspace{1em}(7.1)}$$

根据定义，曲率 $\kappa$ 显然与切线向量 ${\mathbf{r}}^{\prime }(s)$ 有关，那么它们的关系究竟如何？曲率又如何计算？下面就来讨论这些问题。为此，先来证明下述简单命题。

引理

设 $\mathbf{e}(t)$ 为定义在 $U(t_0)\subseteq \mathbb{R}$ 上的单位向量值函数，$\Delta \mathbf{e}=\mathbf{e}(t_0+\Delta t)-\mathbf{e}(t_0)$，$\Delta \theta$ 为 $\mathbf{e}(t_0)$ 与 $\mathbf{e}(t_0+\Delta t)$ 的夹角，则

lim_(Δ t → 0) (|Δ e|)/(Δ θ) = 1 (7.7)$$\underset{\Delta t\to 0}{\lim }\frac{\Vert \Delta \mathbf{e}\Vert }{\Delta \theta }=1\text{\hspace{1em}(7.7)}$$

证明

把向量 $\mathbf{e}(t)$ 的起点固定在空间一点（例如坐标原点），则当 $t$ 在 $U(t_0)$ 中变化时，$\mathbf{e}(t)$ 的终点轨迹便是单位球面上的一段曲线。设 $\mathbf{e}(t_0)$ 与 $\mathbf{e}(t_0+\Delta t)$ 的终点分别为 $Q_0$ 与 $Q$，则 $\Delta \mathbf{e}=\overrightarrow{Q_{0}Q}$（图 5.37）。于是，对于任意固定的 $Q$，在 $\Delta ◯Q_{0}Q$ 内有

|Δ e| = |overrightarrowQ₀ Q| = 2 |e(t₀)| sin (Δ θ)/(2) = 2 sin (Δ θ)/(2)$$\Vert \Delta \mathbf{e}\Vert =\Vert \overrightarrow{Q_{0}Q}\Vert =2\Vert \mathbf{e}(t_0)\Vert \sin \frac{\Delta \theta }{2}=2\sin \frac{\Delta \theta }{2}$$

从而

lim_(Δ t → 0) (|Δ e|)/(Δ θ) = lim_(Δ θ → 0) (2 sin (Δ θ)/(2))/(Δ θ) = 1$$\underset{\Delta t\to 0}{\lim }\frac{\Vert \Delta \mathbf{e}\Vert }{\Delta \theta }=\underset{\Delta \theta \to 0}{\lim }\frac{2\sin \frac{\Delta \theta }{2}}{\Delta \theta }=1$$

图 5.37 单位向量在球面上的运动

由上述引理立即可以导出曲率 $\kappa$ 的计算公式。

定理 7.1

设空间曲线 $\Gamma$ 的方程为 $\mathbf{r}=\mathbf{r}(s)$，$s$ 为自然参数，$\mathbf{r}(s)$ 二阶可导，则 $\Gamma$ 上点 $\mathbf{r}(s)$ 处的曲率为

κ(s) = |r''(s)| (7.8)$$\kappa (s)=\Vert {\mathbf{r}}^{\prime \prime }(s)\Vert \text{\hspace{1em}(7.8)}$$

证明

我们知道 $\mathbf{T}(s)={\mathbf{r}}^{\prime }(s)$ 是 $\Gamma$ 在点 $\mathbf{r}(s)$ 处与 $\Gamma$ 正向一致的单位切向量。设 $\mathbf{T}(s)$ 与 $\mathbf{T}(s+\Delta s)$ 的夹角为 $\Delta \theta$，由曲率的定义及 (7.7) 式，注意到 $\Delta \mathbf{T}=\Delta \mathbf{e}$，可得

κ(s) = lim_(Δ s → 0) | (Δ θ)/(Δ s) | = lim_(Δ s → 0) | (Δ T)/(Δ θ) | · lim_(Δ s → 0) | (Δ θ)/(Δ s) |$$\kappa (s)=\underset{\Delta s\to 0}{\lim }\left|\frac{\Delta \theta }{\Delta s}\right|=\underset{\Delta s\to 0}{\lim }\left|\frac{\Delta \mathbf{T}}{\Delta \theta }\right|\cdot \underset{\Delta s\to 0}{\lim }\left|\frac{\Delta \theta }{\Delta s}\right|$$= lim_(Δ s → 0) ( | (Δ T)/(Δ θ) | (Δ θ)/(Δ s) ) = lim_(Δ s → 0) | (Δ T)/(Δ s) | = lim_(Δ s → 0) | (Δ T)/(Δ s) |$$=\underset{\Delta s\to 0}{\lim }\left(\left|\frac{\Delta \mathbf{T}}{\Delta \theta }\right|\frac{\Delta \theta }{\Delta s}\right)=\underset{\Delta s\to 0}{\lim }\left|\frac{\Delta \mathbf{T}}{\Delta s}\right|=\underset{\Delta s\to 0}{\lim }\left|\frac{\Delta \mathbf{T}}{\Delta s}\right|$$= |T'(s)| = |r''(s)|$$=\Vert {\mathbf{T}}^{\prime }(s)\Vert =\Vert {\mathbf{r}}^{\prime \prime }(s)\Vert$$

几何解释

注意 $\kappa$ 等于 $\mathbf{r}(s)$ 的切向量 $\mathbf{T}(s)$ 的导数的模。由于 $\mathbf{T}(s)$ 刻画切线的方向，因此曲率 $\kappa$ 可看作是切线方向对弧长的转动率，转动越“快”则曲率越大。另一方面，由于切向量 $\mathbf{T}$ 在密切平面上，因此曲率 $\kappa$ 实际上反映了曲线各点在密切平面上的弯曲程度。

推论

设空间曲线 $\Gamma$ 由一般的参数方程 $\mathbf{r}=\mathbf{r}(t)$ 给定，$\mathbf{r}(t)$ 二阶可导且 $\dot{\mathbf{r}}(t)\ne 0$，则 $\Gamma$ 在点 $\mathbf{r}(t)$ 的曲率为

κ(t) = frac|dotr(t) × ddotr(t)||dotr(t)|³ (7.9)$$\kappa (t)=\frac{\Vert \dot{\mathbf{r}}(t)\times \ddot{\mathbf{r}}(t)\Vert }{\Vert \dot{\mathbf{r}}(t){\Vert }^3}\text{\hspace{1em}(7.9)}$$

证明

由链式法则可知

r'(s) = (dr)/(ds) = (dr)/(dt) · (dt)/(ds), r''(s) = (d²r)/(ds²) = (d²r)/(dt²) · ((dt)/(ds))² + (dr)/(dt) · (d²t)/(ds²) (7.10)$${\mathbf{r}}^{\prime }(s)=\frac{\mathrm{d}\mathbf{r}}{\mathrm{d}s}=\frac{\mathrm{d}\mathbf{r}}{\mathrm{d}t}\cdot \frac{\mathrm{d}t}{\mathrm{d}s},{\mathbf{r}}^{\prime \prime }(s)=\frac{{\mathrm{d}}^2\mathbf{r}}{\mathrm{d}{s}^2}=\frac{{\mathrm{d}}^2\mathbf{r}}{\mathrm{d}{t}^2}\cdot {\left(\frac{\mathrm{d}t}{\mathrm{d}s}\right)}^2+\frac{\mathrm{d}\mathbf{r}}{\mathrm{d}t}\cdot \frac{{\mathrm{d}}^{2}t}{\mathrm{d}{s}^2}\text{\hspace{1em}(7.10)}$$

注意到 $\Vert {\mathbf{r}}^{\prime }(s)\Vert =1$，利用向量积的定义与 (7.2) 式，由 (7.8) 式可得

κ(s) = |r''(s)| = |r''(s) × r'(s)|$$\kappa (s)=\Vert {\mathbf{r}}^{\prime \prime }(s)\Vert =\Vert {\mathbf{r}}^{\prime \prime }(s)\times {\mathbf{r}}^{\prime }(s)\Vert$$

将 (7.10) 式代入上式右端并注意到 $\dot{s}(t)=\Vert \dot{\mathbf{r}}(t)\Vert$，得

κ(t) = | ((dt)/(ds))³ ddotr × dotr + (d²t)/(ds²) (dt)/(ds) dotr × dotr | = frac|dotr × ddotr||dotr|³$$\kappa (t)=\Vert {\left(\frac{\mathrm{d}t}{\mathrm{d}s}\right)}^3\ddot{\mathbf{r}}\times \dot{\mathbf{r}}+\frac{{\mathrm{d}}^{2}t}{\mathrm{d}{s}^2}\frac{\mathrm{d}t}{\mathrm{d}s}\dot{\mathbf{r}}\times \dot{\mathbf{r}}\Vert =\frac{\Vert \dot{\mathbf{r}}\times \ddot{\mathbf{r}}\Vert }{\Vert \dot{\mathbf{r}}{\Vert }^3}$$

对于平面曲线 $\mathbf{r}=(x(t),y(t),0)$，由 (7.9) 式易得

κ = frac|xddoty - yddotx|[(dotx)² + (doty)²]³/2 (7.11)$$\kappa =\frac{|x\ddot{y}-y\ddot{x}|}{[(\dot{x}{)}^2+(\dot{y}{)}^2{]}^{3/2}}\text{\hspace{1em}(7.11)}$$

对于平面曲线 $y=y(x)$，即 $\mathbf{r}=(x,y(x),0)$，若 $y(x)$ 二阶可导，由 (7.11) 式易得

κ = frac|y''|[1 + (y')²]³/2 (7.12)$$\kappa =\frac{|y^{\prime \prime }|}{[1+(y^{\prime }{)}^2{]}^{3/2}}\text{\hspace{1em}(7.12)}$$

例 7.2 证明：曲线 $L$ 为直线的充分必要条件是其曲率处处为零

证 必要性 设 $L$ 为直线，取其长度 $s$ 为参数，得其方程为 $\mathbf{r}={\mathbf{r}}_0+s{\mathbf{e}}_1$（${\mathbf{e}}_1$ 为直线 $L$ 的单位方向向量），于是 ${\mathbf{r}}^{\prime }={\mathbf{e}}_1,{\mathbf{r}}^{\prime \prime }=0$，故由 (7.8) 式知 $\kappa (s)\equiv 0$。

充分性 设 $L$ 的曲率 $\kappa \equiv 0$，即 $\Vert {\mathbf{r}}^{\prime \prime }(s)\Vert \equiv 0$，从而 ${\mathbf{r}}^{\prime \prime }(s)=0$。令 $\mathbf{r}(s)=(x(s),y(s),z(s))$，便有 $x^{\prime \prime }(s)=0,y^{\prime \prime }(s)=0,z^{\prime \prime }(s)=0$。积分两次得 $x(s)=c_{1}s+c_1^{\prime },y(s)=c_{2}s+c_2^{\prime },z(s)=c_{3}s+c_3^{\prime }$。这就是以 $s$ 为参数的空间直线的参数方程。

例 7.3 求半径为 $R$ 的圆的曲率

解法一 取圆心为坐标原点 $O$，建立平面直角坐标系。圆的参数方程为

x = Rcos t y = Rsin t (0 ≤slant t ≤slant 2π)$$\{\begin{array}{l}x=R\cos t\\ y=R\sin t\end{array}(0⩽t⩽2\pi )$$

故 $\dot{x}=-R\sin t,\ddot{x}=-R\cos t,\dot{y}=R\cos t,\ddot{y}=-R\sin t$。代入公式 (7.11) 中可算得 $\kappa =\frac{1}{R}$。

解法二 在圆周上任取两点 $P_1,P_2$，圆在点 $P_1$ 与 $P_2$ 切线的夹角 $\Delta \theta$ 等于 $\angle P_{1}O{P}_2$ 所夹的圆心角，其中 $O$ 为圆心。于是弧长 $\Delta s=\text{arc }{P}_1{P}_2=R\Delta \theta$。由曲率定义直接可得

κ = lim_(Δ s → 0) | (Δ θ)/(Δ s) | = 1/R$$\kappa =\underset{\Delta s\to 0}{\lim }\left|\frac{\Delta \theta }{\Delta s}\right|=\frac{1}{R}$$

由此可见，圆周上各点的弯曲程度相同，且曲率等于半径的倒数，半径愈小，曲率愈大。

例 7.4 求螺旋线 $\boldsymbol{r}=(a\cos t, a\sin t, kt)$ 的曲率

解 由于 $\dot{\mathbf{r}}(t)=(-a\sin t,a\cos t,k),\ddot{\mathbf{r}}(t)=(-a\cos t,-a\sin t,0)$，从而

dotr × ddotr = (kasin t, -kacos t, a²)$$\dot{\mathbf{r}}\times \ddot{\mathbf{r}}=(ka\sin t,-ka\cos t,a^2)$$

$\Vert \dot{\mathbf{r}}\times \ddot{\mathbf{r}}\Vert =a\sqrt{a^2+k^2},\Vert \dot{\mathbf{r}}\Vert =\sqrt{a^2+k^2}$。于是由公式 (7.9) 得

κ = (a)/(a² + k²)$$\kappa =\frac{a}{a^2+k^2}$$

2. 曲率半径与曲率圆

如上所述，曲率是曲线在各点弯曲程度的度量。曲线在不同的点处其弯曲程度可能不同，为了更形象地表示曲线在一点处的弯曲程度，下面引入曲率圆的概念。

曲率半径与曲率圆

设 $P$ 为平面曲线 $\Gamma$ 上一点，过点 $P$ 在 $\Gamma$ 所在平面上作此曲线的法线，在法线一侧的法线上取点 $Q$ 使 $\Vert \overrightarrow{PQ}\Vert =\frac{1}{\kappa }$，其中 $\kappa$ 是 $\Gamma$ 在点 $P$ 的曲率。以 $Q$ 为圆心、$\Vert \overrightarrow{PQ}\Vert$ 为半径的圆称为平面曲线 $\Gamma$ 在点 $P$ 处的曲率圆。

如果 $\Gamma$ 是一条空间曲线，那么在点 $P$ 与 $\Gamma$ 最贴近的切平面是点 $P$ 处的密切平面 $\pi$。在点 $P$ 处曲线的主法线上取点 $Q$，使 $\overrightarrow{PQ}$ 的正向与 $\mathbf{N}(s)$ 的正向相同，且使 $\Vert \overrightarrow{PQ}\Vert =\frac{1}{\kappa }$。以 $Q$ 为圆心、$\frac{1}{\kappa }$ 为半径且在平面 $\pi$ 上的圆称为 $\Gamma$ 在点 $P$ 处的曲率圆或密切圆（图 5.38）。此曲率圆的圆心 $Q$ 与半径 $R$ 分别称为曲线 $\Gamma$ 在点 $P$ 处的曲率中心与曲率半径。若取弧长作参数，则曲率中心的向径 $r_0$ 与曲率半径 $R$ 分别为

r₀ = rₚ(s) + R(s)N(s) (7.13)$$r_0={\mathbf{r}}_p(s)+R(s)\mathbf{N}(s)\text{\hspace{1em}(7.13)}$$R = (1)/(κ) (7.14)$$R=\frac{1}{\kappa }\text{\hspace{1em}(7.14)}$$

其中 ${\mathbf{r}}_p$ 为点 $P$ 的向径。

图 5.38 空间曲线的曲率圆

当 $\Gamma$ 为平面曲线时，密切平面就是 $\Gamma$ 所在的平面，而主法线就是 $\Gamma$ 的法线。这时 (7.13) 与 (7.14) 就分别是平面曲线 $\Gamma$ 在点 $P$ 处的曲率中心的向径和曲率半径的表达式。

例 7.5 铁路转弯时的摇晃问题

解 列车在拐弯时将受到向心力的作用，如果向心力的变化不连续，则会产生摇晃。由力学知识可知，轨道上一点处向心力的大小为 $\frac{m{v}^2}{R}$，其中 $m$ 是物体的质量，$v$ 是运动的速度，$R$ 是轨道在该点处的曲率半径。如果列车由直道直接进入圆弧形轨道（图 5.39），尽管轨道是光滑连接的，但由于直线的曲率半径为无穷大，因此在直道末端 $B$ 处，向心力的大小会发生跳跃，这就会导致摇晃震动。为了减小摇晃，必须让轨道的曲率半径 $R$ 连续变化。

图 5.39 轨道连接处的向心力跳跃

容易求得立方抛物线 $y=a{x}^3(a>0)$ 在任一点的曲率半径为

R = frac(1 + 9a²x⁴)³/26a|x|$$R=\frac{(1+9a^2{x}^4{)}^{3/2}}{6a|x|}$$

当 $x\ne 0$ 时，$R$ 随 $x$ 连续变化，且当 $x\to 0$ 时 $R$ 趋于无穷大。因此，如果在修建铁路时，先在直道末端接上一段立方抛物线 $BC$（图 5.40），再在立方抛物线上选取适当的点 $C$ 处与圆弧形轨道 $CD$ 相接，使此立方抛物线在点 $C$ 处的曲率半径近似于圆弧 $CD$ 的半径。这样，在从直轨 $AB$ 转弯轨 $BCD$ 时，曲率半径 $R$ 将接近连续变化，从而减小列车的摇晃震动。

例 7.6 砂轮直径的选择问题

解 为了在磨削时不使砂轮磨削到工件里面去，即多磨掉不应该磨去的部分，砂轮的半径应小于工件截线上各点曲率半径的最小值。为此先求抛物线 $y=0.4x^2$ 上任一点 $(x,y)$ 处的曲率半径。由于 $y^{\prime }=0.8x,y^{\prime \prime }=0.8$，所以曲率半径 $R$ 为

R = (1)/(κ) = frac(1 + (y')²)³/2|y''| = frac(1 + 0.64x²)³/20.8$$R=\frac{1}{\kappa }=\frac{(1+(y^{\prime }{)}^2{)}^{3/2}}{|y^{\prime \prime }|}=\frac{(1+0.64x^2{)}^{3/2}}{0.8}$$

容易求出抛物线的顶点 $(0,0)$ 处曲率半径 $R$ 取到最小值

min R(x) = (1)/(0.8) = 1.25$$\min R(x)=\frac{1}{0.8}=1.25$$

所以选用砂轮的直径应略小于 $2.5\text{ cm}$。

图 5.40 抛物线过渡段设计

图 5.41 抛物线曲率半径分布

3. 渐伸线与渐屈线

若一条曲线 $C$ 上任一点的切线均为另一条曲线 $\Gamma$ 相应点的法线，则称 $C$ 为 $\Gamma$ 的渐伸线，称 $\Gamma$ 为 $C$ 的渐屈线（图 5.42）。

图 5.42 渐伸线与渐屈线

渐伸线与渐屈线

下面我们来求已知曲线 $C$ 的渐伸线 $\Gamma$ 的方程。设 $C$ 的方程为 $\mathbf{r}=\mathbf{r}(s)$，所求渐伸线的方程为 $\mathbf{\rho }=\mathbf{\rho }(s)$。由渐伸线的定义可知

ρ(s) = r(s) + α(s)T(s) (7.15)$$\mathbf{\rho }(s)=\mathbf{r}(s)+\alpha (s)\mathbf{T}(s)\text{\hspace{1em}(7.15)}$$

其中 $\mathbf{T}(s)$ 为曲线 $C$ 的单位切向量，$\alpha (s)$ 为待定数量值函数。为求 $\alpha (s)$，将 (7.15) 式两端对 $s$ 求导，并注意到 $\mathbf{T}(s)={\mathbf{r}}^{\prime }(s)$，得

ρ'(s) = [1 + α'(s)]r'(s) + α(s)r''(s)$${\mathbf{\rho }}^{\prime }(s)=[1+{\alpha }^{\prime }(s)]{\mathbf{r}}^{\prime }(s)+\alpha (s){\mathbf{r}}^{\prime \prime }(s)$$

两端点乘 ${\mathbf{r}}^{\prime }(s)$，由于 ${\mathbf{r}}^{\prime }(s)\perp {\mathbf{r}}^{\prime \prime }(s)$，可得

ρ'(s) · r'(s) = 1 + α'(s)$${\mathbf{\rho }}^{\prime }(s)\cdot {\mathbf{r}}^{\prime }(s)=1+{\alpha }^{\prime }(s)$$

注意到 ${\mathbf{\rho }}^{\prime }(s)$ 是渐伸线 $\Gamma$ 的切向量，由渐伸线的定义知 ${\mathbf{\rho }}^{\prime }(s)\cdot {\mathbf{r}}^{\prime }(s)=0$，从而有

1 + α'(s) = 0$$1+{\alpha }^{\prime }(s)=0$$

解得 $\alpha (s)=a-s$，其中 $a$ 为任意常数。代入 (7.15) 式便得渐伸线 $\Gamma$ 的方程为

ρ(s) = r(s) + (a - s)T(s)$$\mathbf{\rho }(s)=\mathbf{r}(s)+(a-s)\mathbf{T}(s)$$

对于渐屈线的方程本书不作一般讨论，仅指出如下简单情形：平面曲线的曲率中心轨迹就是此曲线的渐伸线（证明略）。因此，平面曲线 $\mathbf{r}=\mathbf{r}(s)$ 的渐伸线方程就是 (7.13) 式。

例 7.7 求椭圆 $\frac{x^2}{a^2} + \frac{y^2}{b^2} = 1$ 的渐伸线方程

解 若利用公式 (7.13) 需要引入弧长参数比较麻烦。我们直接根据曲率中心的定义来求此渐伸线方程。设椭圆上点的向径为 $\mathbf{r}(t)$，所求渐伸线上点的向径为 $\mathbf{\rho }(t)$。将所给椭圆看成一空间曲线，化成参数方程为

r(t) = (acos t, bsin t, 0), t ∈ [0, 2π]$$\mathbf{r}(t)=(a\cos t,b\sin t,0),t\in [0,2\pi ]$$

其单位切向量为

eₜ(t) = fracdotr(t)|dotr(t)| = (1)/(√(a²sin² t + b²cos² t)) (-asin t, bcos t, 0)$${\mathbf{e}}_t(t)=\frac{\dot{\mathbf{r}}(t)}{\Vert \dot{\mathbf{r}}(t)\Vert }=\frac{1}{\sqrt{a^2{\sin }^{2}t+b^2{\cos }^{2}t}}(-a\sin t,b\cos t,0)$$

从而椭圆上点 $\mathbf{r}(t)$（凹一侧）的单位法向量为

eₙ(t) = (0, 0, 1) × eₜ(t) = (1)/(√(a²sin² t + b²cos² t)) (-bcos t, -asin t, 0)$${\mathbf{e}}_n(t)=(0,0,1)\times {\mathbf{e}}_t(t)=\frac{1}{\sqrt{a^2{\sin }^{2}t+b^2{\cos }^{2}t}}(-b\cos t,-a\sin t,0)$$

于是由曲率中心的定义有

ρ(t) = r(t) + R(t)eₙ(t) (7.16)$$\mathbf{\rho }(t)=\mathbf{r}(t)+R(t){\mathbf{e}}_n(t)\text{\hspace{1em}(7.16)}$$

由 (7.14) 式与 (7.11) 式可得

R(t) = (1)/(κ(t)) = frac(a²sin² t + b²cos² t)³/2ab$$R(t)=\frac{1}{\kappa (t)}=\frac{(a^2{\sin }^{2}t+b^2{\cos }^{2}t{)}^{3/2}}{ab}$$

把 $\mathbf{r}(t),R(t),{\mathbf{e}}_n(t)$ 代入 (7.16) 式后化简即得所求渐伸线的方程为

ρ(t) = ( 1/a(a² - b²)cos³ t, -1/b(a² - b²)sin³ t, 0 ), t ∈ [0, 2π]$$\mathbf{\rho }(t)=\left(\frac{1}{a}(a^2-b^2){\cos }^{3}t,-\frac{1}{b}(a^2-b^2){\sin }^{3}t,0\right),t\in [0,2\pi ]$$

或化成直角坐标方程为

(ax)²/3 + (by)²/3 = (a² - b²)²/3$$(ax{)}^{2/3}+(by{)}^{2/3}=(a^2-b^2{)}^{2/3}$$

7.3 挠率

如上所述，曲率是从曲线 $\Gamma$ 在点 $P$（即 $\mathbf{r}(s)$）处的密切平面上考察 $\Gamma$ 的弯曲程度。对于空间曲线 $\Gamma$ 来说，它可能同时还向点 $P$ 的密切平面之外弯曲，我们把这种弯曲叫做扭曲。现在，我们来考察 $\Gamma$ 在点 $P$ 的扭曲情况，换句话说，就是研究 $\Gamma$ 在点 $P$ 偏离密切平面的情况。

当曲线 $\Gamma$ 扭曲时，随着点在曲线上移动，相应的密切平面的位置将随着之改变，从而它的次法向量也将随之转动。在本节 7.2 段已经指出，曲率 $\kappa (s)=\Vert {\mathbf{T}}^{\prime }(s)\Vert$ 反映了曲线切线方向的转动快慢程度，同理，$\Vert {\mathbf{B}}^{\prime }(s)\Vert$ 反映了曲线次法线方向的转动快慢程度，因而它刻画了曲线偏离密切平面的程度，即曲线扭曲的程度。现在我们来考察向量 ${\mathbf{B}}^{\prime }(s)$。

因为 $\mathbf{B}(s)$ 是单位向量，故 ${\mathbf{B}}^{\prime }(s)\perp \mathbf{B}(s)$。又由 $\mathbf{B}(s)=\mathbf{T}(s)\times \mathbf{N}(s)$ 求导得

B' = T' × N + T × N' (7.17)$${\mathbf{B}}^{\prime }={\mathbf{T}}^{\prime }\times \mathbf{N}+\mathbf{T}\times {\mathbf{N}}^{\prime }\text{\hspace{1em}(7.17)}$$

注意到 $\mathbf{N}=\frac{{\mathbf{r}}^{\prime \prime }}{\Vert {\mathbf{r}}^{\prime \prime }\Vert }=\frac{{\mathbf{T}}^{\prime }}{\Vert {\mathbf{T}}^{\prime }\Vert }$，从而 ${\mathbf{T}}^{\prime }\times \mathbf{N}=0$。代入上式得

B' = T × N' (7.18)$${\mathbf{B}}^{\prime }=\mathbf{T}\times {\mathbf{N}}^{\prime }\text{\hspace{1em}(7.18)}$$

这表明 ${\mathbf{B}}^{\prime }\perp \mathbf{T}$，前已说明 ${\mathbf{B}}^{\prime }\perp \mathbf{B}$，故 ${\mathbf{B}}^{\prime }$ 必与 $\mathbf{B}\times \mathbf{T}$（即 $\mathbf{N}$）是共线的，不妨设

B' = -τ N (7.19)$${\mathbf{B}}^{\prime }=-\tau \mathbf{N}\text{\hspace{1em}(7.19)}$$

其中 $\tau$ 为一数量，为确定 $\tau$，用 $\mathbf{N}$ 与上式两端作内积，得

τ = -B' · N (7.20)$$\tau =-{\mathbf{B}}^{\prime }\cdot \mathbf{N}\text{\hspace{1em}(7.20)}$$

而由 (7.19) 式知 $\Vert {\mathbf{B}}^{\prime }\Vert =|\tau |$。

定义 7.2 (挠率)

为曲线 $\mathbf{r}=\mathbf{r}(s)$ 在点 $\mathbf{r}(s)$ 处的挠率 $\tau (s)$ 定义为

τ(s) = -B'(s) · N(s) (7.20)$$\tau (s)=-{\mathbf{B}}^{\prime }(s)\cdot \mathbf{N}(s)\text{\hspace{1em}(7.20)}$$

根据前面的分析，挠率 $\tau (s)$ 的绝对值 $|\tau (s)|=\Vert {\mathbf{B}}^{\prime }(s)\Vert$，反映了曲线在点 $\mathbf{r}(s)$ 处的扭曲程度。

由定义7.2 可以导出挠率的另一计算公式（证略）。

τ(s) = ([r'(s) r''(s) r'''(s)])/(|r''(s)|²) (7.21)$$\tau (s)=\frac{[{\mathbf{r}}^{\prime }(s){\mathbf{r}}^{\prime \prime }(s){\mathbf{r}}^{\prime \prime \prime }(s)]}{\Vert {\mathbf{r}}^{\prime \prime }(s){\Vert }^2}\text{\hspace{1em}(7.21)}$$

如果曲线的方程为

r = r(t) = (x(t), y(t), z(t)), t₁ ≤slant t ≤slant t₂$$\mathbf{r}=\mathbf{r}(t)=(x(t),y(t),z(t)),t_1⩽t⩽t_2$$

那么，当 $\dot{\mathbf{r}}(t)\times \ddot{\mathbf{r}}(t)\ne 0$ 时，由变换 $s=s(t)$ 不难将曲线的挠率公式 (7.21) 转换（习题5.7($A$) 第 13 题）为

τ(t) = frac[dotr(t) ddotr(t) dddotr(t)]|dotr(t) × ddotr(t)|² (7.22)$$\tau (t)=\frac{[\dot{\mathbf{r}}(t)\ddot{\mathbf{r}}(t)\stackrel{...}{r}(t)]}{\Vert \dot{\mathbf{r}}(t)\times \ddot{\mathbf{r}}(t){\Vert }^2}\text{\hspace{1em}(7.22)}$$

例 7.8 证明：曲线 $\boldsymbol{r}=\boldsymbol{r}(s)$ 是平面曲线的充要条件是其任一点的挠率为零

证 必要性 设 $\mathbf{r}=\mathbf{r}(s)$ 是平面曲线，从而其任一点的密切平面都是同一平面，即曲线 $\mathbf{r}=\mathbf{r}(s)$ 所在的平面，于是次法向量 $\mathbf{B}(s)$ 是一常向量。由 (7.20) 式知 $\tau (s)=-{\mathbf{B}}^{\prime }\cdot \mathbf{N}=0$。

充分性 设 $\tau (s)\equiv 0$，则由 (7.19) 式得 ${\mathbf{B}}^{\prime }\equiv 0$，由其坐标式易知 $\mathbf{B}$ 为一常向量，记为 ${\mathbf{B}}_0$。由于 $\mathbf{T}\perp \mathbf{B}$，因而 $\mathbf{T}\cdot {\mathbf{B}}_0=0$，即 ${\mathbf{r}}^{\prime }(s)\cdot {\mathbf{B}}_0=0$，从而

(r(s) · B₀)' = 0$$(\mathbf{r}(s)\cdot {\mathbf{B}}_0{)}^{\prime }=0$$

将上式两端积分得

r(s) · B₀ = C (常数)$$\mathbf{r}(s)\cdot {\mathbf{B}}_0=C\text{ (常数)}$$

令 $s=0$，得 $C=\mathbf{r}(0)\cdot {\mathbf{B}}_0$，代入上式得

(r(s) - r(0)) · B₀ = 0$$(\mathbf{r}(s)-\mathbf{r}(0))\cdot {\mathbf{B}}_0=0$$

故曲线上任一点 $\mathbf{r}(s)$ 均在过点 $\mathbf{r}(0)$ 且以 ${\mathbf{B}}_0$ 为法向量的平面上，故为平面曲线。

例 7.9 求螺旋线 $\boldsymbol{r}(t)=(a\cos t, a\sin t, kt)$ 的挠率

解 由于 $\dot{\mathbf{r}}(t)=(-a\sin t,a\cos t,k),\ddot{\mathbf{r}}(t)=(-a\cos t,-a\sin t,0),\stackrel{...}{r}(t)=(a\sin t,-a\cos t,0)$，从而

[dotr(t) ddotr(t) dddotr(t)] = -asin t & acos t & k -acos t & -asin t & 0 asin t & -acos t & 0 = ka²$$[\dot{\mathbf{r}}(t)\ddot{\mathbf{r}}(t)\stackrel{...}{r}(t)]=\left|\begin{array}{ccc}-a\sin t& a\cos t& k\\ -a\cos t& -a\sin t& 0\\ a\sin t& -a\cos t& 0\end{array}\right|=k{a}^2$$

$\Vert \dot{\mathbf{r}}\times \ddot{\mathbf{r}}{\Vert }^2=a^2(a^2+k^2)$，所以由公式 (7.22) 得

τ(t) = (ka²)/(a²(a² + k²)) = (k)/(a² + k²)$$\tau (t)=\frac{k{a}^2}{a^2(a^2+k^2)}=\frac{k}{a^2+k^2}$$

5.7 空间曲线的曲率与挠率