上海专业网站建设公司站霸网络,萝岗区网站建设推广,网站建设交什么税,免费网站seo排名优化本文仅供学习使用 本文参考#xff1a; 《机构运动微分几何学分析与综合》-王德伦、汪伟 《微分几何》吴大任 Ch01-4 平面运动微分几何学 3.1 空间曲线微分几何学概述3.1.1 矢量表示3.1.2 Frenet标架 连杆机构中的连杆与连架杆构成运动副#xff0c;该运动副元素的特征点或特… 本文仅供学习使用 本文参考 《机构运动微分几何学分析与综合》-王德伦、汪伟 《微分几何》吴大任 Ch01-4 平面运动微分几何学 3.1 空间曲线微分几何学概述3.1.1 矢量表示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} xx(t) \\ yy(t) \\ zz(t) \\ \end{matrix} \right. ⎩ ⎨ ⎧xx(t)yy(t)zz(t) 式中 t t t为曲线的参数若置换自变量或者消去参数 t t t,则可写成 { y y ( x ) z z ( x ) \left\{ \begin{matrix} yy(x) \\ zz(x) \\ \end{matrix} \right. {yy(x)zz(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 :Rx(t)iy(t)jz(t)k Γ:Rx(t)iy(t)jz(t)k 可以将其简化为: R R ( t ) RR(t) RR(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 :Rr(\varphi ){{e}_{I(\varphi )}}z(\varphi )k Γ:Rr(φ)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 :RR(s),{{s}_{a}}\le s\le {{s}_{b}} Γ:RR(s),sa≤s≤sb 注书中为 Γ : R R ( s ) , s a ≤ a ≤ s b \Gamma :RR(s),{{s}_{a}}\le a\le {{s}_{b}} Γ:RR(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} xx(t),yy(t),zz(t) \\ {{x}^{2}}{{y}^{2}}{{z}^{2}}{{R}^{2}} \\ \end{matrix} \right. {xx(t),yy(t),zz(t)x2y2z2R2 式中 R R R为球面半径 t t t为球面曲线的参数若置换自变量或者消去参数 t t t,可写成 { z z ( x , y ) x 2 y 2 z 2 R 2 \left\{ \begin{matrix} zz(x,y) \\ {{x}^{2}}{{y}^{2}}{{z}^{2}}{{R}^{2}} \\ \end{matrix} \right. {zz(x,y)x2y2z2R2 由于球面曲线上的点始终分布在一球面上因此往往用球面坐标表示曲线为 δ δ ( t ) , φ φ ( t ) , r R \delta \delta (t),\varphi \varphi (t),rR δδ(t),φφ(t),rR 式中 δ \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 δ xR\sin \delta \cos \varphi ,yR\sin \delta \sin \varphi ,zR\cos \delta xRsinδcosφ,yRsinδsinφ,zRcosδ 将上述 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 :RR(t)x(t)iy(t)jz(t)k Γ:RR(t)x(t)iy(t)jz(t)k 若通过圆矢量函数表示球面曲线的矢量方程则为 R R sin δ ( φ ) e I ( φ ) R cos δ ( φ ) k RR\sin \delta (\varphi ){{e}_{I(\varphi )}}R\cos \delta (\varphi )k RRsinδ(φ)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 \\ zz(\varphi ) \\ \end{matrix} \right. ⎩ ⎨ ⎧xr0cosφyr0sinφzz(φ) 式中 r 0 {{r}_{0}} r0为圆柱面半径。若通过圆矢量函数表示圆柱面曲线的矢量方程则为 R r 0 e I ( φ ) z ( φ ) k R{{r}_{0}}{{e}_{I(\varphi )}}z(\varphi )k Rr0eI(φ)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 CR 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 ) kk\left( s \right) ,\tau \tau \left( s \right) kk(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)Δs2!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) 关于弧长参数的各阶导数。
