[足式机器人]Part3机构运动微分几何学分析与综合Ch01-3 平面运动微分几何学——【读书笔记】

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

---

1.2.2 瞬心线-2

（3）平面连杆机构（二自由度开链串联机构）的瞬心线
见例1-4：铰链点B在固定坐标下中的轨迹为圆，连杆平面$({\vec{i}}_m,{\vec{j}}_m)$绕B点旋转，此时轨迹曲线（圆）为曲线${\Gamma }_{Om}$，其弧长$s=a_1\phi$，$ds=a_{1}d\phi$，其曲率为：$k_{Om}=k_a=\frac{1}{a_1}$，为常数（AB连杆长的倒数），所以曲率$k_{Om}$的各阶倒数均为零。

• 圆曲线上Frenet标架的标矢$\vec{\alpha }$在连杆运动坐标系 { B : i ⃗ m , j ⃗ m } \{B:{{{\vec{i}}}_{m}},{{{\vec{j}}}_{m}}\} {B:i m​,j ​m​}中的方向角为$\theta =\frac{\pi }{2}+\phi -\gamma$，从而有$k_a-\dot{\theta }=\dot{\gamma }$，结合前述内容，可得连杆上瞬心的位置为：
  $\{\begin{array}{c}{\theta }_{p0}=n\pi +(\phi -\gamma ),n=0,1\\ r_{p0}=\mp \frac{1}{\dot{\gamma }}\end{array}$
  其中：$\theta =\angle ({\vec{i}}_m,\vec{\alpha })$，$\dot{\gamma }=\frac{d\gamma }{ds}=\frac{d\gamma }{dt}\cdot \frac{dt}{ds}=\frac{d\gamma }{dt}\cdot \frac{dt}{l_{1}d\phi }=\frac{{\omega }_2}{l_1{\omega }_1}$，${\omega }_1,{\omega }_2$为连杆AB和连杆BC的角速度。
• 则将上述确定的连杆点$P_0({\theta }_{p0},r_{p0})$在连杆平面（动坐标系）与机架平面（定坐标系）上分别绘制出动瞬心线${\pi }_m$与定瞬心线${\pi }_f$，则动瞬心线${\pi }_m$在动坐标系 { B : i ⃗ m , j ⃗ m } \{B:{{{\vec{i}}}_{m}},{{{\vec{j}}}_{m}}\} {B:i m​,j ​m​}，以及定瞬心线${\pi }_f$在定坐标系 { A : i ⃗ f , j ⃗ f } \{A:{{{\vec{i}}}_{f}},{{{\vec{j}}}_{f}}\} {A:i f​,j ​f​}中的矢量方程分别为：
  $\{\begin{array}{c}{\pi }_m:{\vec{R}}_m=\frac{1}{\dot{\gamma }}\vec{\beta }=-\frac{\cos (\phi -\gamma )}{\dot{\gamma }}{\vec{i}}_m-\frac{\sin (\phi -\gamma )}{\dot{\gamma }}{\vec{j}}_m\\ {\pi }_f:{\vec{R}}_f=(\frac{1}{\dot{\gamma }}-a_1)\vec{\beta }=(a_1-\frac{1}{\dot{\gamma }})(\cos \phi {\vec{i}}_f+\sin \phi {\vec{j}}_f)\end{array}$
  其中：连杆坐标系与机架坐标系之间的转换关系为：

例1-5： 平面曲柄摇杆机构的动瞬心线${\pi }_m$与定瞬心线${\pi }_f$
四杆机构如下图所示，机架AD上建立机架坐标系 { A ; i ⃗ f , j ⃗ f } \{A;{{{\vec{i}}}_{f}},{{{\vec{j}}}_{f}}\} {A;i f​,j ​f​}，连杆BC上建立构件坐标系 { B ; i ⃗ m , j ⃗ m } \{B;{{{\vec{i}}}_{m}},{{{\vec{j}}}_{m}}\} {B;i m​,j ​m​}，为描述连杆平面的运动，需要首先确定出铰链点B的位移，与构件AB的转角，或者说需要对该四杆机构进行位移求解——建立闭环矢量方程：（本质上是通过约定坐标系，对机构进行运动学分析，其意义是为了在指定坐标系下——机架坐标系——描述连杆机构各可动构件/运动副/构件点的平面运动）

（此处各角度约定与向量约定依据原文，实际运动学分析约定可参考文章：待补充）

建立闭环矢量方程：${\vec{a}}_1+{\vec{a}}_2={\vec{a}}_3+{\vec{a}}_4$
将矢量方程分别投影到坐标轴${\vec{i}}_f$与${\vec{j}}_f$上，消去矢量${\vec{a}}_3$的方位角可得：
${a_1}^2+{a_2}^2+{a_4}^2-{a_3}^2+2a_1{a}_2\cos (\phi -\gamma )-2a_1{a}_4\cos \phi -2a_2{a}_4\cos \gamma =0$

• 上式求得$\gamma \leftrightarrow \phi$的关系式，而在求解过程中，角度的关系式及变化范围不尽相同，往往通过Grashof准则来对平面铰链四杆机构进行分类：对于Grashof运动链，取与最短杆相邻的杆作机架的曲柄摇杆机构作为基本机构；对于非Grashof运动链，取以最长杆作机架的双摇杆机构作为基本机构。对所有四杆机构，均可由基本机构转化得到，但需用相应的变换数学模型及装配模式。

将上式对铰链点B的轨迹曲线弧长$s$求导，简化后可得：
$\dot{\gamma }=\frac{d\gamma }{ds}=\frac{a_2\sin (\phi -\gamma )-a_4\sin \phi }{a_2{a}_4\sin \gamma +a_1{a}_2\sin (\phi -\gamma )}$

• 由于位移方程是非线性的，随不同运动位置求解容易出现奇异，建议按四杆机构的基本形式建立位移方程求解，然后变换机架得到所使用的四杆机构类型位移解，肖大准教授早有论述，在此不赘述。

（1）曲柄摇杆机构
给定一个平面曲柄摇杆机构，其杆长分别为：$a_1=1,a_2=3,a_3=3.5,a_4=5$，其瞬心线如下所示：

（2）双曲柄机构
将上图的四杆机构进行变换：AB杆为机架，BC杆为原动件，则四杆机构变换称双曲柄机构，此时四杆机构变换称双曲柄机构。

例1-6： 平面曲柄滑块机构的动瞬心线${\pi }_m$与定瞬心线${\pi }_f$

其闭环矢量方程为：
${\vec{a}}_1+{\vec{a}}_2=\vec{S}+\vec{E}$
将矢量方程向固定坐标系 { A : i ⃗ f , j ⃗ f } \{A:{{{\vec{i}}}_{f}},{{{\vec{j}}}_{f}}\} {A:i f​,j ​f​}的坐标轴${\vec{j}}_f$上投影，可得：
$a_1\sin \phi +a_2\sin \gamma -e=0$

• 上式确定了平面曲柄滑块机构连杆BC的倾角与输入杆AB的转角之间的关系

对铰链点B的轨迹曲线弧长$s$求导，化简后可得：
$\dot{\gamma }=\frac{d\gamma }{ds}=-\frac{\cos \phi }{a_2\cos \gamma }$
令：$a_1=1,a_2=2.5,e=1.2$，其连杆相对机架的转角函数为：$\gamma =\gamma (\phi )$——

1.2.3-1 点轨迹的Euler-Savary公式

如上节所述，运动刚体$\sum \ast$相对于固定刚体$\sum$的平面运动，可以通过运动刚体$\sum \ast$相对固定刚体$\sum$的线位移及转角（即：线/角位移）的三种方式来表示 （齐次旋转矩阵、 { R o m : i ⃗ m , j ⃗ m } \{{{R}_{om}}:{{{\vec{i}}}_{m}},{{{\vec{j}}}_{m}}\} {Rom​:i m​,j ​m​}、 { R o m : α ⃗ , β ⃗ } \{{{R}_{om}}:\vec{\alpha },\vec{\beta }\} {Rom​:α ,β ​}），每种方式都由三个独立参数（自由度）确定刚体”的运动。然而，这三种方式中却没有直观表达刚体的运动本质，或体现刚体的运动学与轨迹几何学之间的内在联系。 由相伴运动方法推导出刚体运动的动瞬心线${\pi }_m$与定瞬心线${\pi }_f$，这两条瞬心线在瞬心点处活动标架重合、微弧长相等，并随刚体运动而相切地纯滚动。 显然，瞬心线本身隐含着刚体相对运动的内在联系信息，以其为出发点研究刚体$\sum \ast$上图形，如点或者直线在固定刚体上轨迹的几何性质，无疑是既顺畅自然又简洁直观的方式，已有的研究结果表明相伴运动方法将大大简化研究问题，并使得其曲率公式形式优美。

本节仍沿用上节的相伴运动方法，分别以运动刚体平面上的动瞬心线${\pi }_m$和固定平面上的定瞬心线${\pi }_f$作为原曲线，将运动刚体上的点分别看作与动、定瞬心线作相伴运动，可得到瞬时在瞬心线活动标架中描述的点轨迹的曲率公式以及高阶曲率特征。

任意瞬时平面运动刚体$\sum \ast$上的一点$P$总是与瞬心点相互对应，表明：点$P$的轨迹曲线和瞬心线可以分别视为彼此的相伴曲线——在固定坐标系 { O f : i ⃗ f , j ⃗ f } \{{{O}_{f}}:{{{\vec{i}}}_{f}},{{{\vec{j}}}_{f}}\} {Of​:i f​,j ​f​}中，将点$P$的轨迹曲线${\Gamma }_P$看成定瞬心线${\pi }_f$的相伴曲线——以定瞬心线${\pi }_f$为原曲线来描述点$P$在固定平面上的轨迹${\Gamma }_P$：

则有：${\Gamma }_P:{\vec{R}}_p={\vec{R}}_f+v_1{\vec{E}}_{1f}+v_2{\vec{E}}_{2f}$

• 其中：$(v_1,v_2)$为运动刚体$\sum \ast$上点$P$在定瞬心线活动标架 { R f : E ⃗ 1 f , E ⃗ 2 f } \{{{R}_{f}}:{{{\vec{E}}}_{1f}},{{{\vec{E}}}_{2f}}\} {Rf​:E 1f​,E 2f​}中的相对坐标，若以极坐标$(r,\alpha )$表示，则有：$\{\begin{array}{c}{v}_1=r\cos \alpha \\ v_2=r\sin \alpha \end{array}$

将上式对动、定瞬心线的弧长$\sigma$求导（定瞬心线弧长），可得：$\frac{d{\vec{R}}_P}{d\sigma }=(1+\frac{d{v}_1}{d{\sigma }_f}-k_f{v}_2){\vec{E}}_{1f}+(\frac{d{v}_2}{d{\sigma }_f}+k_f{v}_1){\vec{E}}_{2f}$

考察点$P$在动瞬心线${\pi }_m$上的Frenet标架 { R m : E ⃗ 1 m , E ⃗ 2 m } \{{{R}_{m}}:{{{\vec{E}}}_{1m}},{{{\vec{E}}}_{2m}}\} {Rm​:E 1m​,E 2m​}上的运动，则有：
${\vec{R}}_{Pm}={\vec{R}}_m+u_1{\vec{E}}_{1m}+u_2{\vec{E}}_{2m}$

• 其中：$(u_1,u_2)$为点$P$在动瞬心线${\pi }_m$上是Frenet标架 { R m : E ⃗ 1 m , E ⃗ 2 m } \{{{R}_{m}}:{{{\vec{E}}}_{1m}},{{{\vec{E}}}_{2m}}\} {Rm​:E 1m​,E 2m​}中的坐标参数

即从运动刚体$\sum \ast$上固接的运动坐标系 { O m : i ⃗ m , j ⃗ m } \{{{O}_{m}}:{{{\vec{i}}}_{m}},{{{\vec{j}}}_{m}}\} {Om​:i m​,j ​m​}原点引矢量${\vec{R}}_m$到瞬心点（动瞬心线${\pi }_m$上），再由瞬心点（动瞬心线上Frenet标架）考察点$P(u_1,u_2)$的运动，即把点$P$的运动看成与动瞬心线${\pi }_m$的相伴运动，

将上式对动瞬心线弧长${\sigma }_m$进行求导可得：$\frac{d{\vec{R}}_{Pm}}{d{\sigma }_m}=(1+\frac{d{u}_1}{d{\sigma }_m}-k_m{u}_2){\vec{E}}_{1m}+(\frac{d{u}_2}{d{\sigma }_m}+k_m{u}_1){\vec{E}}_{2m}$

由于点$P$是运动刚体$\sum \ast$上的固定点，即在运动刚体坐标系 { O m : i ⃗ m , j ⃗ m } \{{{O}_{m}}:{{{\vec{i}}}_{m}},{{{\vec{j}}}_{m}}\} {Om​:i m​,j ​m​}中，点$P$的绝对运动变换率为零，即$\frac{d{\vec{R}}_{Pm}}{d{\sigma }_m}=0$，从而点$P$的运动满足Cesaro不动点条件：
$\{\begin{array}{c}1+\frac{d{u}_1}{d{\sigma }_m}-k_m{u}_2=0\\ \frac{d{u}_2}{d{\sigma }_m}+k_m{u}_1=0\end{array}\Rightarrow \{\begin{array}{c}\frac{d{u}_1}{d{\sigma }_m}=k_m{u}_2-1\\ \frac{d{u}_2}{d{\sigma }_m}=-k_m{u}_1\end{array}$
上式表明了在活动标架沿动瞬心线${\pi }_m$运动时，观察到点$P$的变化（上式是点$P$为瞬心的特款）：犹如小船在小河中前进，艄公所观察到小河岸边景物的变化一样，该变化的性质由小河的弯曲程度体现，那么$(\frac{d{u}_1}{d{\sigma }_m},\frac{d{u}_2}{d{\sigma }_m})$是运动刚体$\sum \ast$上点$P$在动瞬心线${\pi }_m$的Frenet标架 { R ⃗ m : E ⃗ 1 m , E ⃗ 2 m } \{{{{\vec{R}}}_{m}}:{{{\vec{E}}}_{1m}},{{{\vec{E}}}_{2m}}\} {R m​:E 1m​,E 2m​}的相对运动分量也就是动瞬心线${\pi }_m$本身几何性质的体现。

且知，动瞬心线与定瞬心线的Frenet 标架重合，在瞬心点处相切地纯滚动，即同一点$P$在瞬时重合的活动标架 { R ⃗ f : E ⃗ 1 f , E ⃗ 2 f } \{{{{\vec{R}}}_{f}}:{{{\vec{E}}}_{1f}},{{{\vec{E}}}_{2f}}\} {R f​:E 1f​,E 2f​}与 { R ⃗ m : E ⃗ 1 m , E ⃗ 2 m } \{{{{\vec{R}}}_{m}}:{{{\vec{E}}}_{1m}},{{{\vec{E}}}_{2m}}\} {R m​:E 1m​,E 2m​}中具有相同的投影坐标及其变化率，即∶
$u_1=v_1,u_2=v_2,\frac{d{u}_1}{d{\sigma }_m}=\frac{d{v}_1}{d{\sigma }_f},\frac{d{u}_2}{d{\sigma }_m}=\frac{d{v}_2}{d{\sigma }_f}$
进而可得：$\frac{d{\vec{R}}_P}{d\sigma }=(k_f-k_m)(-v_2{\vec{E}}_{1f}+v_1{\vec{E}}_{2f})=rk\ast (-\sin \alpha {\vec{E}}_{1f}+\cos \alpha {\vec{E}}_{2f})$
进而可得：$d{\sigma }_P=\left|\frac{d{\vec{R}}_P}{d\sigma }\right|d\sigma =rk\ast d\sigma$

将$\{\begin{array}{c}{v}_1=r\cos \alpha \\ v_2=r\sin \alpha \end{array}$对$\sigma$求导，并利用不动点条件化简，得：
$\frac{dr}{d\sigma }=-\cos \alpha ,\frac{d\alpha }{d\sigma }=-k_m+\frac{\sin \alpha }{r}$

详细推导：
$\{\begin{array}{c}{v}_1=r\cos \alpha \\ v_2=r\sin \alpha \end{array}\Rightarrow \{\begin{array}{c}\frac{d{v}_1}{d\sigma }=\frac{dr}{d\sigma }\cos \alpha -r\sin \alpha \frac{d\alpha }{d\sigma }=k_m{u}_2-11\\ \frac{d{v}_2}{d\sigma }=\frac{dr}{d\sigma }\sin \alpha +r\cos \alpha \frac{d\alpha }{d\sigma }=-k_m{u}_{1}2\end{array}$
$\{\begin{array}{c}1\cdot \sin \alpha +2\cdot \cos \alpha \\ v_1=r\cos \alpha ,v_2=r\sin \alpha \end{array}\Rightarrow \frac{dr}{d\sigma }=(k_m{u}_2-1)\cos \alpha -k_m{u}_1\sin \alpha =-\cos \alpha$
将$\frac{dr}{d\sigma }=-\cos \alpha$带入$2$中，可得：
$-\cos \alpha \sin \alpha +r\cos \alpha \frac{d\alpha }{d\sigma }=-k_{m}r\cos \alpha \Rightarrow \frac{d\alpha }{d\sigma }=-k_m+\frac{\sin \alpha }{r}$

进而建立P点轨迹曲线${\Gamma }_P$的Frenet标架 { R ⃗ P : E ⃗ 1 ∗ , E ⃗ 2 ∗ } \{{{{\vec{R}}}_{P}}:{{{\vec{E}}}_{1}}*,{{{\vec{E}}}_{2}}*\} {R P​:E 1​∗,E 2​∗}为：
$\{\begin{array}{c}{\vec{E}}_1\ast =\frac{{\vec{R}}_P}{d{\sigma }_P}=\frac{{\vec{R}}_P}{d{\sigma }_f}/\left|\frac{{\vec{R}}_P}{d{\sigma }_f}\right|=-\sin \alpha {\vec{E}}_{1f}+\cos \alpha {\vec{E}}_{2f}\\ {\vec{E}}_2\ast =\vec{k}\times {\vec{E}}_1\ast =-\cos \alpha {\vec{E}}_{1f}-\sin \alpha {\vec{E}}_{2f}\end{array}$
进而可得：
$k_{\Gamma }=\frac{d{\vec{E}}_1\ast }{d{\sigma }_P}\cdot {\vec{E}}_2\ast =\frac{1}{r}+\frac{\sin \alpha }{r^{2}k\ast }=\frac{r+\sin \alpha /k\ast }{r^2}$

详细推导：
$\frac{d{\vec{E}}_1\ast }{d{\sigma }_P}=\frac{1}{r^{2}k\ast }(-\cos \alpha \frac{d\alpha }{ds}{\vec{E}}_{1f}-\sin \alpha k_f{\vec{E}}_{2f}-\sin \alpha \frac{d\alpha }{ds}{\vec{E}}_{2f}-\cos \alpha k_f{\vec{E}}_{1f})$
$\Rightarrow \frac{d{\vec{E}}_1\ast }{d{\sigma }_P}=\frac{-\cos \alpha }{r^{2}k\ast }(\frac{d\alpha }{ds}+k_f){\vec{E}}_{1f}+\frac{-\sin \alpha }{r^{2}k\ast }(\frac{d\alpha }{ds}+k_f){\vec{E}}_{2f}$
$k_{\Gamma }=\frac{d{\vec{E}}_1\ast }{d{\sigma }_P}\cdot {\vec{E}}_2\ast =(\frac{-\cos \alpha }{r^{2}k\ast }(\frac{d\alpha }{ds}+k_f){\vec{E}}_{1f}+\frac{-\sin \alpha }{r^{2}k\ast }(\frac{d\alpha }{ds}+k_f){\vec{E}}_{2f})\cdot (-\cos \alpha {\vec{E}}_{1f}-\sin \alpha {\vec{E}}_{2f})=\frac{\frac{d\alpha }{ds}+k_f}{r^{2}k\ast }=\frac{k\ast +\frac{\sin \alpha }{r}}{r^{2}k\ast }$

由此得到运动刚体运动刚体$\sum \ast$上点$P$瞬时在固定刚体$\sum$上的轨迹曲线的曲率公式，也就是著名的Euler-Savary公式。若将曲率$k$，用曲率半径$\rho$的倒数代替，则得到另一种表达形式∶${\rho }_{\Gamma }(r+\frac{\sin \alpha }{k\ast })=r^2$
即平面Euler-Savary公式描述了平面运动刚体上点的位置（在Frenet标架 { R ⃗ f : E ⃗ 1 f , E ⃗ 2 f } \{{{{\vec{R}}}_{f}}:{{{\vec{E}}}_{1f}},{{{\vec{E}}}_{2f}}\} {R f​:E 1f​,E 2f​}中的极坐标$r,\alpha$）、轨迹的曲率半径（中心）${\rho }_{\Gamma }$以及动、定瞬心线诱导曲率$k\ast$之间的关系。

将点P转换到运动坐标系 { O m : i ⃗ m , j ⃗ m } \{{{O}_{m}}:{{{\vec{i}}}_{m}},{{{\vec{j}}}_{m}}\} {Om​:i m​,j ​m​}中，以直角坐标$(x_{pm},y_{pm})$表示，则可得到点$P$轨迹曲线${\Gamma }_P$曲率的另一表达形式∶
$\{\begin{array}{c}{k}_{\Gamma }=\frac{F}{G}\\ F={(x_{pm}-a)}^2+{(y_{pm}-b)}^2-\frac{D^2}{4}\\ G={[{(x_{pm}+\frac{\sin \theta }{k_{Om}-\dot{\theta }})}^2+{(y_{pm}-\frac{\cos \theta }{k_{Om}-\dot{\theta }})}^2]}^{\frac{3}{2}}\\ a=-\frac{({\dot{k}}_{Om}-\ddot{\theta })\cos \theta +(2k_{Om}-\dot{\theta })(k_{Om}-\dot{\theta })\sin \theta }{2{(k_{Om}-\dot{\theta })}^3}\\ b=-\frac{({\dot{k}}_{Om}-\ddot{\theta })\sin \theta +(2k_{Om}-\dot{\theta })(k_{Om}-\dot{\theta })\cos \theta }{2{(k_{Om}-\dot{\theta })}^3}\end{array}$

详细推导：
设：$d{\sigma }_m={(\frac{{\dot{\theta }}^2}{a^2}+\frac{b^2}{a^4})}^{\frac{1}{2}}ds\Rightarrow \frac{ds}{d{\sigma }_m}=\frac{a^2}{{({\dot{\theta }}^2{a}^2+b^2)}^{\frac{1}{2}}}=A$
则将动瞬心线${\pi }_m$的Frenet标架基矢量改写为：$\{\begin{array}{c}{\vec{R}}_m=\frac{1}{a}(-\sin \theta {\vec{i}}_m+\cos \theta {\vec{j}}_m)\\ {\vec{E}}_{1m}=[\frac{-\dot{\theta }a\cos \theta +b\sin \theta }{a^2}{\vec{i}}_m-\frac{\dot{\theta }a\sin \theta +b\cos \theta }{a^2}{\vec{j}}_m]A\\ {\vec{E}}_{2m}=[\frac{\dot{\theta }a\sin \theta +b\cos \theta }{a^2}{\vec{i}}_m+\frac{-\dot{\theta }a\cos \theta +b\sin \theta }{a^2}{\vec{j}}_m]A\end{array}$
且已知：$\{\begin{array}{c}{u}_1=r\cos \alpha \\ u_2=r\sin \alpha \end{array},k_{\Gamma }=\frac{r+\sin \alpha /k\ast }{r^2},k\ast =\frac{a^3}{{({\dot{\theta }}^2{a}^2+b^2)}^{\frac{1}{2}}}$
${\vec{R}}_{Pm}={\vec{R}}_m+u_1{\vec{E}}_{1m}+u_2{\vec{E}}_{2m}$，其坐标参数可表示为$(x_{Pm},y_{Pm})$，可得：
$\{\begin{array}{c}{x}_{Pm}=-\frac{\sin \theta }{a}+\frac{-\dot{\theta }a\cos \theta +b\sin \theta }{{({\dot{\theta }}^2{a}^2+b^2)}^{\frac{1}{2}}}r\cos \alpha +\frac{\dot{\theta }a\sin \theta +b\cos \theta }{{({\dot{\theta }}^2{a}^2+b^2)}^{\frac{1}{2}}}r\sin \alpha \frac{}{}\frac{}{}1\\ y_{Pm}=\frac{\cos \theta }{a}-\frac{\dot{\theta }a\sin \theta +b\cos \theta }{{({\dot{\theta }}^2{a}^2+b^2)}^{\frac{1}{2}}}r\cos \alpha +\frac{-\dot{\theta }a\cos \theta +b\sin \theta }{{({\dot{\theta }}^2{a}^2+b^2)}^{\frac{1}{2}}}r\sin \alpha \frac{}{}\frac{}{}2\end{array}$
则对上式进行处理：
$\{\begin{array}{c}1\frac{}{}\sin \theta \text{-}2\frac{}{}\cos \theta :x_{Pm}\sin \theta -y_{Pm}\cos \theta =-\frac{1}{a}+\frac{b}{{({\dot{\theta }}^2{a}^2+b^2)}^{\frac{1}{2}}}r\cos \alpha +\frac{\dot{\theta }a}{{({\dot{\theta }}^2{a}^2+b^2)}^{\frac{1}{2}}}r\sin \alpha \frac{}{}\frac{}{}3\\ 1\frac{}{}\cos \theta \text{+}2\frac{}{}\sin \theta :x_{Pm}\cos \theta +y_{Pm}\sin \theta =\frac{\dot{\theta }a}{{({\dot{\theta }}^2{a}^2+b^2)}^{\frac{1}{2}}}r\cos \alpha +\frac{b}{{({\dot{\theta }}^2{a}^2+b^2)}^{\frac{1}{2}}}r\sin \alpha \frac{}{}\frac{}{}4\end{array}$
其中根据$3\frac{}{}\frac{\dot{\theta }a}{b}+4$：
$(x_{Pm}\sin \theta -y_{Pm}\cos \theta )\frac{\dot{\theta }a}{b}+(x_{Pm}\cos \theta +y_{Pm}\sin \theta )=-\frac{\dot{\theta }}{b}+\frac{\frac{{\dot{\theta }}^2{a}^2}{b}+b}{{({\dot{\theta }}^2{a}^2+b^2)}^{\frac{1}{2}}}r\sin \alpha$
进而求解$\sin \alpha$：
$\sin \alpha =\frac{(x_{Pm}\sin \theta -y_{Pm}\cos \theta )\dot{\theta }a+(x_{Pm}\cos \theta +y_{Pm}\sin \theta )b+\dot{\theta }}{r{({\dot{\theta }}^2{a}^2+b^2)}^{\frac{1}{2}}}$
且知：$r^2=|\overrightarrow{P{P}_O}{|}^2=|{\vec{R}}_{Pm}-{\vec{R}}_m{|}^2={(x_{Pm}+\frac{\sin \theta }{a})}^2+{(y_{Pm}-\frac{\cos \theta }{a})}^2$
解得：$k_{\Gamma }=\frac{r+\sin \alpha /k\ast }{r^2}=\frac{r^2+\frac{x_{Pm}(\dot{\theta }a\sin \theta +b\cos \theta )+y_{Pm}(-\dot{\theta }a\cos \theta +b\sin \theta )+\dot{\theta }}{a^3}}{r^3}$
令：$G=r^3$，整理分子后可得：$k_{\Gamma }=\frac{{(x_{Pm}+B)}^2+{(y_{Pm}+C)}^2+\frac{1}{4{k^{\ast }}^2}}{G},B=\frac{b\cos \theta +(2a+\dot{\theta })a\sin \theta }{2a^3},C=\frac{b\sin \theta -(2a+\dot{\theta })a\cos \theta }{2a^3}$

若在某一瞬时点$P$轨迹曲线${\Gamma }_P$的曲率$k_{\Gamma }$为零，即该点三个无限接近位置在一直线上，或者说轨迹曲线${\Gamma }_P$与直线二阶接触（可参考吴大任《微分几何》3.6＆3.7），称其为该轨迹曲线${\Gamma }_P$上的拐点。令$\frac{r+\sin \alpha /k\ast }{r^2}$为零，可得∶
$$r+D\sin \alpha =0$$
其中：$D=\frac{1}{k\ast }$。这是一个圆的方程，表示在该瞬时轨迹产生拐点的运动刚体上点都汇集在该圆上，称该圆为运动刚体上的拐点圆。

其直径为$\left|1/k\ast \right|$且圆心位于瞬心线法线上，并且与瞬心线相切。随着刚体平面运动的不同瞬时（或位置），拐点圆在刚体上的位置和大小随之变化。虽然拐点圆过瞬心点$P_0$，但属于拐点圆上奇点——其轨迹线的曲率为无穷大，即轨迹线上出现尖点。

由Euler-Savary公式可知，在瞬心线的Frenet 标架参考系下，可求得点$P$到瞬心$P_0$的矢量$\overrightarrow{P{P}_0}$、点$P$到其轨迹曲线${\Gamma }_P$，曲率中心$O_P$的矢量$\overrightarrow{P{O}_P}$（曲率半径${\rho }_{\Gamma }$），点$P$到拐点圆上相应点$P_I$的矢量$\overrightarrow{P{P}_I}$，并且三个矢量存在关系式∶$\overrightarrow{P{P}_I}\cdot \overrightarrow{P{O}_P}={(\overrightarrow{P{P}_0})}^2$，由于该式右端恒为正值，从而曲率中心$O_P$与拐点$P_I$必位于点$P$的同侧。若某一瞬时，运动刚体上的点恰为拐点，即$\overrightarrow{P{P}_I}=0$，其轨迹曲线的曲率中心应在无穷远处。

补充：《高等机构学》白师贤 第六章 平面无线接近位置运动几何学（待补充）

平面运动刚体上点的轨迹线的曲率仅仅反映了该点无限接近三位置的几何特征，为揭示更大范围上该点的运动几何特征，需研究平面运动点的高阶曲率特征。当平面运动点的轨迹曲线的（高阶）曲率为特定值或者为零时，其轨迹曲线具有特定的几何意义，研究这些点在运动刚体上的位置，对机构的分析与综合具有重要的指导意义。