[足式机器人]Part3机构运动微分几何学分析与综合Ch03-1 空间约束曲线与约束曲面微分几何学——【读书笔记】

本文仅供学习使用
本文参考：
《机构运动微分几何学分析与综合》-王德伦、汪伟
《微分几何》吴大任

Ch01-4 平面运动微分几何学

• [3.1 空间曲线微分几何学概述](#3.1 空间曲线微分几何学概述)
• • [3.1.1 矢量表示](#3.1.1 矢量表示)
  • [3.1.2 Frenet标架](#3.1.2 Frenet标架)

---

连杆机构中的连杆与连架杆构成运动副，该运动副元素的特征点或特征线在机架坐标系 中的运动轨迹曲线或曲面 称为约束曲线或约束曲面，是联系刚体运动与机构运动综合的桥梁，其几何性质是机构运动综合的理论基础，既是曲线与曲面的几何学研究内容，也是连杆机构运动几何学分析与综合的课题。然而，研究曲线与曲面的几何学，微分几何学方法无疑是自然而然的选择，将其与机构运动学结合，形成以点与线的运动方式 研究约束曲线与曲面几何性质，为机构运动几何学分析与综合提供理论依据。

为方便阅读后续内容，在第3.1和第3.2节简单概述微分几何学基本知识；采用微分几何量方法研究连杆机构中典型而又重要的约束曲线与约束曲面，称为空间约束曲线与约束曲面微分几何学。

3.1 空间曲线微分几何学概述

3.1.1 矢量表示

在直角坐标中表达一条空间曲线 Γ \Gamma Γ时，有：
{ x = x ( t ) y = y ( t ) z = z ( t ) \left\{ \begin{matrix} x=x(t) \\ y=y(t) \\ z=z(t) \\ \end{matrix} \right. ⎩ ⎨ ⎧x=x(t)y=y(t)z=z(t)

式中， t t t为曲线的参数，若置换自变量或者消去参数 t t t,则可写成：
{ y = y ( x ) z = z ( x ) \left\{ \begin{matrix} y=y(x) \\ z=z(x) \\ \end{matrix} \right. {y=y(x)z=z(x)

或者写成隐函数形式：
{ F 1 ( x , y , z ) = 0 F 2 ( x , y , z ) = 0 \left\{ \begin{matrix} {{F}{1}}(x,y,z)=0 \\ {{F}{2}}(x,y,z)=0 \\ \end{matrix} \right. {F1(x,y,z)=0F2(x,y,z)=0

若将上述 x , y , z x,y,z x,y,z置于空间固定坐标系 { O : i , j , k } \{O:i,j,k\} {O:i,j,k}中，则曲线 Γ \Gamma Γ以参数 t t t表示的矢量方程为：
Γ : R = x ( t ) i + y ( t ) j + z ( t ) k \Gamma :R=x(t)i+y(t)j+z(t)k Γ:R=x(t)i+y(t)j+z(t)k

可以将其简化为:
R = R ( t ) R=R(t) R=R(t)
式(3.4)和式(3.5)为空间曲线 Γ \Gamma Γ的矢量表达式， t t t为曲线 Γ \Gamma Γ的一般参数。在 第1章平面曲线的微分几何学 中引入了圆矢量函数用来描述曲线的矢量方程，使得形式简洁并便于计算。因此对于空间曲线 Γ \Gamma Γ的矢量方程式(3.4),可以选择任意两个坐标轴上的分量用圆矢量函数 进行描述。例如，将曲线 Γ \Gamma Γ上任意点的矢径在坐标平面 O − i j O-ij O−ij上的投影矢量用圆矢量函数描述，如下图所示。

则其矢量方程可以写出另一种形式：
Γ : R = r ( φ ) e I ( φ ) + z ( φ ) k \Gamma :R=r(\varphi ){{e}_{I(\varphi )}}+z(\varphi )k Γ:R=r(φ)eI(φ)+z(φ)k

对于空间曲线 Γ \Gamma Γ，弧长参数 s s s为其自然参数，且与一般参数 t t t的关系为：
s = ∫ t a t b ∣ d R d t ∣ d t , d s = ∣ d R ∣ = ( d x d t ) 2 + ( d y d t ) 2 + ( d z d t ) 2 d t s=\int_{{{t}{a}}}^{{{t}{b}}}{\left| \frac{dR}{dt} \right|dt,ds=\left| dR \right|}=\sqrt{{{(\frac{dx}{dt})}^{2}}+{{(\frac{dy}{dt})}^{2}}+{{(\frac{dz}{dt})}^{2}}}dt s=∫tatb dtdR dt,ds=∣dR∣=(dtdx)2+(dtdy)2+(dtdz)2 dt

空间曲线 Γ \Gamma Γ的矢量方程用弧长参数 s s s表示为：
Γ : R = R ( s ) , s a ≤ s ≤ s b \Gamma :R=R(s),{{s}{a}}\le s\le {{s}{b}} Γ:R=R(s),sa≤s≤sb

注 ：书中为 Γ : R = R ( s ) , s a ≤ a ≤ s b \Gamma :R=R(s),{{s}{a}}\le a\le {{s}{b}} Γ:R=R(s),sa≤a≤sb
【例3-1】 球面曲线如下图所示：

对于球面曲线 Γ \Gamma Γ,习惯于将直角坐标系 { O : i , j , k } \{O:i,j,k\} {O:i,j,k}原点置于球心，则用直角坐标表示为：
{ x = x ( t ) , y = y ( t ) , z = z ( t ) x 2 + y 2 + z 2 = R 2 \left\{ \begin{matrix} x=x(t),y=y(t),z=z(t) \\ {{x}^{2}}+{{y}^{2}}+{{z}^{2}}={{R}^{2}} \\ \end{matrix} \right. {x=x(t),y=y(t),z=z(t)x2+y2+z2=R2

式中， R R R为球面半径， t t t为球面曲线的参数，若置换自变量或者消去参数 t t t,可写成：
{ z = z ( x , y ) x 2 + y 2 + z 2 = R 2 \left\{ \begin{matrix} z=z(x,y) \\ {{x}^{2}}+{{y}^{2}}+{{z}^{2}}={{R}^{2}} \\ \end{matrix} \right. {z=z(x,y)x2+y2+z2=R2

由于球面曲线上的点始终分布在一球面上，因此往往用球面坐标表示曲线为：
δ = δ ( t ) , φ = φ ( t ) , r = R \delta =\delta (t),\varphi =\varphi (t),r=R δ=δ(t),φ=φ(t),r=R

式中， δ \delta δ是由原点0到曲线上点 P P P的有向线段 O P OP OP与 k k k的夹角； φ \varphi φ是 O P OP OP在 O − i j O-ij O−ij面上的投影与i的夹角， δ \delta δ和 φ \varphi φ的取值范围分别为 [ 0 , π ] [0,\pi ] [0,π]和 [ 0 , 2 π ] [0,2\pi ] [0,2π]。点 P P P在坐标系 { O : i , j , k } \{O:i,j,k\} {O:i,j,k}中的球面坐标与直角坐标之间具有如下转换关系：
x = R sin ⁡ δ cos ⁡ φ , y = R sin ⁡ δ sin ⁡ φ , z = R cos ⁡ δ x=R\sin \delta \cos \varphi ,y=R\sin \delta \sin \varphi ,z=R\cos \delta x=Rsinδcosφ,y=Rsinδsinφ,z=Rcosδ

将上述 x , y , z x,y,z x,y,z置于坐标系 { O : i , j , k } \{O:i,j,k\} {O:i,j,k}中，则球面曲线以参数 t t t表示的矢量方程为：
Γ : R = R ( t ) = x ( t ) i + y ( t ) j + z ( t ) k \Gamma :R=R(t)=x(t)i+y(t)j+z(t)k Γ:R=R(t)=x(t)i+y(t)j+z(t)k

若通过圆矢量函数表示球面曲线的矢量方程，则为：
R = R sin ⁡ δ ( φ ) e I ( φ ) + R cos ⁡ δ ( φ ) k R=R\sin \delta (\varphi ){{e}_{I(\varphi )}}+R\cos \delta (\varphi )k R=Rsinδ(φ)eI(φ)+Rcosδ(φ)k

比较式(E3-1.1)、式(E3-1.4)与式(E3-1.6)可知，采用矢量表示的球面曲线比其他方式表达要简单的多。
【例3-2】 圆柱面曲线如下图所示：

圆柱面曲线在直角坐标系 { O : i , j , k } \{O:i,j,k\} {O:i,j,k}中的方程为：
{ x = r 0 cos ⁡ φ y = r 0 sin ⁡ φ z = z ( φ ) \left\{ \begin{matrix} x={{r}{0}}\cos \varphi \\ y={{r}{0}}\sin \varphi \\ z=z(\varphi ) \\ \end{matrix} \right. ⎩ ⎨ ⎧x=r0cosφy=r0sinφz=z(φ)

式中， r 0 {{r}{0}} r0为圆柱面半径。若通过圆矢量函数表示圆柱面曲线的矢量方程，则为：
R = r 0 e I ( φ ) + z ( φ ) k R={{r}{0}}{{e}_{I(\varphi )}}+z(\varphi )k R=r0eI(φ)+z(φ)k

3.1.2 Frenet标架

空间曲线 Γ ⃗ : R ⃗ = R ⃗ ( s ) \vec{\varGamma}:\vec{R}=\vec{R}\left( s \right) Γ :R =R (s) 在任意点 P P P处有两个无限接近位置的点连线组成切线，其单位切矢 α ⃗ ( s ) = d R ⃗ ( s ) d s \vec{\alpha}\left( s \right) =\frac{\mathrm{d}\vec{R}\left( s \right)}{\mathrm{d}s} α (s)=dsdR (s)始终指向曲线弧长增加的方向，将切矢 α ⃗ ( s ) \vec{\alpha}\left( s \right) α (s) 对弧长参数求导，可得：

d α ⃗ ( s ) d s = k ( s ) β ⃗ ( s ) \frac{\mathrm{d}\vec{\alpha}\left( s \right)}{\mathrm{d}s}=k\left( s \right) \vec{\beta}\left( s \right) dsdα (s)=k(s)β (s)

其中， k ( s ) k\left( s \right) k(s) 称为曲线 Γ ⃗ \vec{\varGamma} Γ 在点 P P P 处的曲率 ，即三个无限接近位置点构成空间曲线在该点处的密切平面 ，曲率是空间曲线在密切平面内的弯曲程度，体现了曲线的切矢的倾斜角对弧长参数的变化率。与平面曲线曲率不同，空间曲线的曲率非负 。 β ⃗ ( s ) \vec{\beta}\left( s \right) β (s) 称为曲线 Γ ⃗ \vec{\varGamma} Γ 在点 P P P 处的主法矢 ，指向了曲线在该点的曲率中心 。当 k ( s ) ≠ 0 k\left( s \right) \ne 0 k(s)=0 时，其倒数 ρ ( s ) = 1 / k ( s ) \rho \left( s \right) =1/k\left( s \right) ρ(s)=1/k(s) 称为曲线 Γ ⃗ \vec{\varGamma} Γ 的曲率半径 ，则曲线 Γ ⃗ \vec{\varGamma} Γ 曲率中心 C C C的矢量为：

R ⃗ C = R ⃗ P + ρ ⋅ β ⃗ \vec{R}{\mathrm{C}}=\vec{R}{\mathrm{P}}+\rho \cdot \vec{\beta} R C=R P+ρ⋅β

由空间曲线 Γ ⃗ \vec{\varGamma} Γ 在点 P P P 处的切矢 α ⃗ ( s ) \vec{\alpha}\left( s \right) α (s) 和主法矢 β ⃗ ( s ) \vec{\beta}\left( s \right) β (s) 可以构建矢量 γ ⃗ ( s ) = α ⃗ ( s ) × β ⃗ ( s ) \vec{\gamma}\left( s \right) =\vec{\alpha}\left( s \right) \times \vec{\beta}\left( s \right) γ (s)=α (s)×β (s) ，称之为曲线的副法矢 ，从而在空间曲线 Γ ⃗ \vec{\varGamma} Γ 上构造了单位右手系正交标架 { R ⃗ ( s ) ; α ⃗ ( s ) , β ⃗ ( s ) , γ ⃗ ( s ) } \left\{ \vec{R}\left( s \right) ;\vec{\alpha}\left( s \right) ,\vec{\beta}\left( s \right) ,\vec{\gamma}\left( s \right) \right\} {R (s);α (s),β (s),γ (s)} ，称为曲线 Γ ⃗ \vec{\varGamma} Γ 在点 P P P 的Frenet标架。

对于空间曲线 Γ ⃗ \vec{\varGamma} Γ 在 P P P 点的Frenet标架 { R ⃗ ( s ) ; α ⃗ ( s ) , β ⃗ ( s ) , γ ⃗ ( s ) } \left\{ \vec{R}\left( s \right) ;\vec{\alpha}\left( s \right) ,\vec{\beta}\left( s \right) ,\vec{\gamma}\left( s \right) \right\} {R (s);α (s),β (s),γ (s)} ，其中标矢 α ⃗ ( s ) \vec{\alpha}\left( s \right) α (s) 和 β ⃗ ( s ) \vec{\beta}\left( s \right) β (s) 确定了密切平面 ， β ⃗ ( s ) \vec{\beta}\left( s \right) β (s) 和 γ ⃗ ( s ) \vec{\gamma}\left( s \right) γ (s) 确定的平面称为法平面 ，而 α ⃗ ( s ) \vec{\alpha}\left( s \right) α (s) 和 γ ⃗ ( s ) \vec{\gamma}\left( s \right) γ (s) 确定的平面称为从切平面。可见Frenet标架由三个同空间曲线紧密联系的向量所组成，其微分运算公式为：

{ d R ⃗ ( s ) d s = α ⃗ ( s ) d α ⃗ ( s ) d s = k ( s ) β ⃗ ( s ) d β ⃗ ( s ) d s = − k ( s ) α ⃗ ( s ) + τ ( s ) γ ⃗ ( s ) d γ ⃗ ( s ) d s = − τ ( s ) β ⃗ ( s ) \begin{cases} \begin{array}{c} \frac{\mathrm{d}\vec{R}\left( s \right)}{\mathrm{d}s}=\vec{\alpha}\left( s \right)\\ \frac{\mathrm{d}\vec{\alpha}\left( s \right)}{\mathrm{d}s}=k\left( s \right) \vec{\beta}\left( s \right)\\ \end{array}\\ \begin{array}{c} \frac{\mathrm{d}\vec{\beta}\left( s \right)}{\mathrm{d}s}=-k\left( s \right) \vec{\alpha}\left( s \right) +\tau \left( s \right) \vec{\gamma}\left( s \right)\\ \frac{\mathrm{d}\vec{\gamma}\left( s \right)}{\mathrm{d}s}=-\tau \left( s \right) \vec{\beta}\left( s \right)\\ \end{array}\\ \end{cases} ⎩ ⎨ ⎧dsdR (s)=α (s)dsdα (s)=k(s)β (s)dsdβ (s)=−k(s)α (s)+τ(s)γ (s)dsdγ (s)=−τ(s)β (s)

其中， τ ( s ) \tau \left( s \right) τ(s) 称为空间曲线 Γ ⃗ \vec{\varGamma} Γ 在点 P P P 处的挠率 ，它衡量了曲线在点 $P$ 的(密切平面)副法矢 $\vec{\gamma}\left( s \right)$ 倾斜角对弧长的变化率，从而描述了曲线在该点偏离密切平面的程度。上式也称为空间曲线的Frenet公式。

由Frenet公式可以得到空间曲线 Γ ⃗ \vec{\varGamma} Γ 曲率 k k k 和挠率 τ \tau τ 的表达式为：

k = ∣ d 2 R ⃗ ( s ) d s 2 ∣ , τ = ( d R ⃗ ( s ) d s , d 2 R ⃗ ( s ) d s 2 , d 3 R ⃗ ( s ) d s 3 ) / ∣ d 2 R ⃗ ( s ) d s 2 ∣ 2 k=\left| \frac{\mathrm{d}^2\vec{R}\left( s \right)}{\mathrm{d}s^2} \right|,\tau =\left( \frac{\mathrm{d}\vec{R}\left( s \right)}{\mathrm{d}s},\frac{\mathrm{d}^2\vec{R}\left( s \right)}{\mathrm{d}s^2},\frac{\mathrm{d}^3\vec{R}\left( s \right)}{\mathrm{d}s^3} \right) /\left| \frac{\mathrm{d}^2\vec{R}\left( s \right)}{\mathrm{d}s^2} \right|^2 k= ds2d2R (s) ,τ=(dsdR (s),ds2d2R (s),ds3d3R (s))/ ds2d2R (s) 2

若空间曲线 Γ ⃗ \vec{\varGamma} Γ 是以一般参数 t t t 进行描述的，则其曲率 k ( s ) k(s) k(s) 和挠率 τ ( s ) \tau \left( s \right) τ(s) 的表达式为：
k = ∣ d R ⃗ d t × d 2 R ⃗ d t 2 ∣ / ∣ d R ⃗ d t ∣ 3 , τ = ( d R ⃗ d t , d 2 R ⃗ d t 2 , d 3 R ⃗ d t 3 ) / ( d R ⃗ d t × d 2 R ⃗ d t 2 ) 2 k=\left| \frac{\mathrm{d}\vec{R}}{\mathrm{d}t}\times \frac{\mathrm{d}^2\vec{R}}{\mathrm{d}t^2} \right|/\left| \frac{\mathrm{d}\vec{R}}{\mathrm{d}t} \right|^3,\tau =\left( \frac{\mathrm{d}\vec{R}}{\mathrm{d}t},\frac{\mathrm{d}^2\vec{R}}{\mathrm{d}t^2},\frac{\mathrm{d}^3\vec{R}}{\mathrm{d}t^3} \right) /\left( \frac{\mathrm{d}\vec{R}}{\mathrm{d}t}\times \frac{\mathrm{d}^2\vec{R}}{\mathrm{d}t^2} \right) ^2 k= dtdR ×dt2d2R / dtdR 3,τ=(dtdR ,dt2d2R ,dt3d3R )/(dtdR ×dt2d2R )2

对于空间曲线来说，曲率 k ( s ) k(s) k(s) 和挠率 τ ( s ) \tau \left( s \right) τ(s) 不依赖于坐标系的选定。是空间曲线的不变量 ，能够唯一地确定空间曲线，可以将 k = k ( s ) , τ = τ ( s ) k=k\left( s \right) ,\tau =\tau \left( s \right) k=k(s),τ=τ(s) 称为空间曲线的自然方程。于是有：

定理3.1 ：在区间 0 ⩽ s ⩽ l 0\leqslant s\leqslant l 0⩽s⩽l上任意给定连续可微函数 $k\left( s \right) >0 $ 和连续函数 τ ( s ) \tau \left( s \right) τ(s) 以及初始右手系正交标架 { R ⃗ 0 ; α ⃗ 0 , β ⃗ 0 , γ ⃗ 0 } \left\{ \vec{R}_0;\vec{\alpha}_0,\vec{\beta}_0,\vec{\gamma}_0 \right\} {R 0;α 0,β 0,γ 0} ,则一定有且仅有一条以 s s s 为弧长、以 k ( s ) k\left( s \right) k(s) 为曲率、 τ ( s ) \tau \left( s \right) τ(s) 为挠率的空间有向曲线。

建立了空间曲线 Γ ⃗ \vec{\varGamma} Γ 上 P P P 点处的Frenet标架 { R ⃗ ( s ) ; α ⃗ ( s ) , β ⃗ ( s ) , γ ⃗ ( s ) } \left\{ \vec{R}\left( s \right) ;\vec{\alpha}\left( s \right) ,\vec{\beta}\left( s \right) ,\vec{\gamma}\left( s \right) \right\} {R (s);α (s),β (s),γ (s)} ,可将曲线 Γ ⃗ \vec{\varGamma} Γ 在点 P P P的邻域内按照泰勒公式展开。假定曲线 Γ ⃗ \vec{\varGamma} Γ 在点 P P P 处的弧长为 s s s,则有：

R ⃗ ( s + Δ s ) = R ⃗ ( s ) + d R ⃗ ( s ) d s Δ s + 1 2 ! d 2 R ⃗ ( s ) d s 2 ( Δ s ) 2 + ⋯ + 1 n ! d n R ⃗ ( s ) d s n ( Δ s ) n + ε n ( s , Δ s ) ( Δ s ) n \vec{R}\left( s+\varDelta s \right) =\vec{R}\left( s \right) +\frac{\mathrm{d}\vec{R}\left( s \right)}{\mathrm{d}s}\varDelta s+\frac{1}{2!}\frac{\mathrm{d}^2\vec{R}\left( s \right)}{\mathrm{d}s^2}\left( \varDelta s \right) ^2+\cdots +\frac{1}{n!}\frac{\mathrm{d}^n\vec{R}\left( s \right)}{\mathrm{d}s^n}\left( \varDelta s \right) ^n+\varepsilon _{\mathrm{n}}\left( s,\varDelta s \right) \left( \varDelta s \right) ^n R (s+Δs)=R (s)+dsdR (s)Δs+2!1ds2d2R (s)(Δs)2+⋯+n!1dsndnR (s)(Δs)n+εn(s,Δs)(Δs)n

式中， lim ⁡ Δ s → 0 ε n ( s , Δ s ) = 0 , d R ⃗ d s = α ⃗ , d 2 R ⃗ d s 2 = k β ⃗ , d 3 R ⃗ d s 3 = − k 2 α ⃗ + d k d s β ⃗ + k τ γ ⃗ \lim_{\varDelta s\rightarrow 0} \varepsilon _{\mathrm{n}}\left( s,\varDelta s \right) =0,\frac{\mathrm{d}\vec{R}}{\mathrm{d}s}=\vec{\alpha},\frac{\mathrm{d}^2\vec{R}}{\mathrm{d}s^2}=k\vec{\beta},\frac{\mathrm{d}^3\vec{R}}{\mathrm{d}s^3}=-k^2\vec{\alpha}+\frac{\mathrm{d}k}{\mathrm{d}s}\vec{\beta}+k\tau \vec{\gamma} limΔs→0εn(s,Δs)=0,dsdR =α ,ds2d2R =kβ ,ds3d3R =−k2α +dsdkβ +kτγ ，并以此可以得到矢径 R ⃗ ( s ) \vec{R}\left( s \right) R (s) 关于弧长参数的各阶导数。