Ch01-4 平面运动微分几何学
连杆机构中的连杆与连架杆构成运动副，该运动副元素的特征点或特征线在机架坐标系中的运动轨迹曲线或曲面称为约束曲线或约束曲面，是联系刚体运动与机构运动综合的桥梁，其几何性质是机构运动综合的理论基础，既是曲线与曲面的几何学研究内容，也是连杆机构运动几何学分析与综合的课题。然而，研究曲线与曲面的几何学，微分几何学方法无疑是自然而然的选择，将其与机构运动学结合，形成以点与线的运动方式研究约束曲线与曲面几何性质，为机构运动几何学分析与综合提供理论依据。
为方便阅读后续内容，在第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
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) 称为空间曲线的自然方程。于是有：
建立了空间曲线 Γ ⃗ \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) 关于弧长参数的各阶导数。