第 6 章 曲面上的曲线:测地线与协变变化
在地球仪上,从北京到巴黎画一条“尽量短”的路线,为什么它通常不是地图上的水平直线?把曲线当作空间曲线,它的加速度有一部分指向球心;把曲线当作球面上的路径,又要区分它是在曲面内拐弯,还是因为球面本身弯曲才产生了这部分加速度。
这一章把这两个角度分开。法曲率描述曲面在所选方向上的外在弯曲,测地曲率描述曲线相对于曲面内部的转弯。协变导数只保留切向变化,Christoffel 符号则把这种切向变化写成坐标方程。测地线就是测地曲率为零的曲面曲线;在球面上,它们正是大圆的弧段。
曲面上的曲线与加速度分解
设 r(u,v) 是 C2 正则参数化,曲面上的 C2 曲线写成
α(t)=r(u(t),v(t)).
只要 α′(t)=0,可以用弧长 s 作为局部参数。速度为
dtdα=ru
单位切向量定义为
T=dsdα=∣α
沿曲线选定单位法向 N,再定义切平面内与 T 正交的单位向量
J(T)=N×T.
这里的 J(T) 依赖曲面的定向;换成 −N 时,J(T) 也换成相反方向。
因为 ⟨T,T⟩=1,对弧长求导可得
0=dsd⟨T,T⟩=2⟨ds
所以 dT/ds 垂直于 T。在三维空间中,垂直于 T 的平面由 J(T) 与 N 张成,于是存在唯一的两个数 k,使得
dsdT=kgJ(T)+k
定义
- kg 为测地曲率,它是加速度在曲面切平面内的分量;
- kn 为法曲率,它是加速度沿曲面法向的分量。

图示说明:曲面上的曲线有单位切向 T,切平面内的转弯方向 J(T),以及曲面法向 N。
若使用任意参数 t,不必先求出弧长函数。因为
dsdT=∣α′∣1
而 N 与 α′ 垂直,所以
kn=⟨dsdT
对 kg 同样处理:
kg=⟨dsd
这里的分母不能漏掉。直接用 ⟨α′′,N⟩ 只在曲线已经用单位速度参数化时才等于法曲率。

图示说明:空间曲线的加速度分解为切平面内的测地分量和法向的曲面弯曲分量。
协变导数:只看切向变化
设 X(t) 是沿 α(t) 的切向量场。它可以写成
X(t)=X1(t)ru+X2(t)
把它当作 R3 中的向量场求普通导数,得到
dtdX=dtdX
其中采用重复指标求和约定,r1=ru,r2=r。这个导数可能包含法向分量;曲面内部的观察者无法把法向分量当作沿曲面的变化。于是定义协变导数
∇α˙X=(dtdX
即普通导数在切平面上的正交投影。
若 α 用弧长参数化,X=T,则
∇TT=kgJ(T).
因为 dT/ds=kgJ(T)+knN,投影去掉的恰好是 。所以
∇TT=0⟺kg=0.
这给出了测地线的几何定义:曲线沿自身前进时,不在曲面的切平面内转弯。

图示说明:普通导数先在空间中计算,再投影回曲面切平面,剩下的就是协变导数。
Christoffel 符号的来源
为了把协变导数写成坐标公式,设 q1=u,q2=v,并记
gij=⟨ri,r
由于 rij 是空间向量,它可以唯一分解为切向部分和法向部分:
rij=Γijkrk
其中
b11=e,b12=b21
系数 Γijk 称为 Christoffel 符号。它们不是张量的分量,而是坐标基底 ri 随位置变化时出现的连接系数。
对度量系数求导:
∂igjℓ=∂
写出另外两式:
∂jgiℓ=⟨rj
∂ℓgij=⟨rℓi
C2 条件给出混合偏导相等:rij=rji、、。把前两式相加再减去第三式,交叉项两两抵消,剩下
∂igjℓ+∂
将分解式代入右端。因为 N⊥rℓ,
⟨rij,rℓ⟩=Γ
所以
2Γijkgkℓ=∂
第一基本形式正定,矩阵 (gij) 有逆矩阵 (gij)。两边乘以 g,利用 ,得到
Γijm=2
这个推导同时说明了三个事实:Γijm=Γjim 来自混合偏导相等;它们只由第一基本形式及其一阶导数决定;正定性保证求逆这一步合法。

图示说明:坐标基向量 ru,rv 沿曲面移动时会发生切向和法向变化,Christoffel 符号记录其中的切向部分。
协变导数的坐标表达
沿 α(t) 的切向量场写成 X=Xjrj。由上面的分解式,普通导数的切向部分为
(dtdX)⊤=(
因此
(∇α˙X)k
如果 X=α˙=q˙jrj,则
(∇α˙α˙)k=
这就是测地线方程的左侧。注意这里的 qi 是曲面坐标,不是空间中的笛卡尔坐标。
测地线方程与参数条件
若参数 t 是仿射参数,例如弧长参数或常速参数,测地线满足
q¨k+Γij
为什么需要“仿射参数”这个条件?同一条几何曲线换成任意非线性参数后,速度会变化,普通的二阶方程会多出沿切向量的项。几何上,测地线的正确条件是
∇TT=0
或 kg=0;若使用任意参数,应允许
∇α˙α˙=λ(t)α˙
其中 λ 只反映参数速度的改变。选择仿射参数可以令 λ=0,于是得到上面的标准方程。
由 ∇TT=0 和加速度分解还可见,测地线的空间加速度只有法向分量:
dsdT=knN.
它通常不是 R3 中的直线。平面上的测地线碰巧是直线,是因为平面的法向不随位置改变;在球面上,大圆弧是测地线,但它在空间中是圆而不是直线。

图示说明:测地线方程抵消坐标基底变化造成的假加速度,使协变加速度为零。
完整例题:球面上的纬线与大圆
考虑单位球面 S2 的参数化
r(θ,φ)=(sinφcosθ,sinφsinθ,cosφ),
其中 0<φ<π。求固定余纬度 φ=φ0 的纬线的法曲率与测地曲率,并判断它什么时候是测地线。
方法选择
题目要求区分法曲率和测地曲率,直接把纬线写成空间曲线并使用加速度分解最清楚。若只解 Christoffel 方程,可以判断测地性,却不容易看见曲线的空间加速度怎样分成两个方向。算出曲率后,再用坐标方程验证大圆是同一个结论。
固定 0<φ0<π,记 s0=sinφ、,纬线为
α(t)=(s0cost,s0sint,c
求速度与加速度:
α′(t)=(−s
用 Christoffel 方程复核赤道
球面的第一基本形式来自上一章:
gθθ=sin2φ,gθφ
由 Christoffel 公式可以算出与本题相关的非零系数
Γθφθ=Γφθθ=
赤道取 θ=t,φ=π/2,于是 θ˙=1,φ˙。
- θ 方程为 θ¨+2cotφθ˙φ˙=,代入后是 。
两条方程都成立,和 kg=0 的几何判断一致。这次复核还说明:赤道的测地性不是因为它在空间中“看起来最平”,而是因为球面度量下的协变加速度为零。

图示说明:球面上一条大圆由过球心的平面截出;其切向加速度指向球面法向,因此没有切平面内的转弯。

图示说明:赤道是大圆并且 kg=0,一般纬线是小圆且 ∣kg∣=∣cotφ。
交互可以让曲线在球面上从赤道移动到不同纬度,并实时分解加速度。操作前先预测:球面半径固定时 kn 是否改变,以及哪一个量会在赤道处变成零;观察时分别核对法向分量和切平面分量,而不是只看总加速度大小。
交互可以从球面上一点改变初始切向方向,观察由测地线方程产生的轨迹。尝试让初始方向指向经线方向、赤道方向和一般斜方向,预测哪些轨迹会闭合为大圆,并用“轨迹所在平面是否经过球心”检查结果。
练习
练习 1:测地曲率的含义
下列哪一项等价于一条正则曲面曲线是测地线?
练习 2:分解一条单位速曲线的加速度
某曲面曲线已经用弧长参数化,在某点满足 kg=3,kn=4。求该点的 ∣dT/ds∣,并说明 是否一定沿曲面法向。
答案:因为 J(T) 与 N 都是单位向量且互相正交,
dsdT
练习 3:Christoffel 符号的读取
在某坐标点,度量矩阵为 gij=diag(4,9),且该点附近所有 gij 都是常数。求该点的 Christoffel 符号,并写出测地线方程。
答案:因为所有 ∂igjℓ=0,Christoffel 公式给出所有 Γij。因此测地线方程为
练习 4:球面上的小圆
单位球面上固定余纬度 φ0=π/4 的纬线,求其法曲率和测地曲率的绝对值,并判断它是否是测地线。
答案:由纬线例题,
∣kn∣=1,∣kg∣=∣cot(π/4
练习 5:迁移到任意参数
一条曲面曲线用参数 t 表示,速度为 ∣α′(t)∣=2,且 ⟨α′′(t。求法曲率,并说明若误把 当作法曲率会错在哪里。
答案:任意参数下
kn=∣α′∣2

“测地线是曲面上的直线”只能作为直觉,不能当定义。可靠的判定是 kg=0,或在仿射参数下检查 q。球面大圆在空间中有圆形加速度,却仍然是球面上的测地线;普通纬线则相反,空间上也是圆,但在球面内部有非零的测地转弯。