圆螺旋线求曲率半径

引言

看论文的时候看到给出了螺旋线的曲率半径求解公式：

便想着这个公式是怎么求出来的。

上述螺旋线我们知道h(转一圈，又称为导程)，螺旋的半径是r。

我们可以用参数方程来描述螺旋线：

$$\{\begin{array}{rl}x& =rcos\theta \\ y& =rsin\theta \\ z& =\frac{h}{2\pi }\theta \end{array}$$

为什么z的表达式中系数是$\frac{h}{2\pi }$呢？想一下，$\theta$角从0转到$2\pi$，这样整个旋转了一圈，把螺旋线拉开截面方向就从0到$2\pi r$，纵向方向从0到h，两个都是线性增加，很明显
$$\{\begin{array}{rl}z(0)& =0\\ z(2\pi )& =h\\ z(\theta )& =k\theta \end{array}$$

则有

$$k=\frac{z(2\pi )-z(0)}{2\pi -0}=\frac{h}{2\pi }$$

这就解释了z的表达式中系数k是$\frac{h}{2\pi }$。
写出了参数方程，参考维基百科[1]参数方程的曲率公式（曲率半径是曲率的倒数）：

对于一个以参数化形式给出的空间曲线$\mathbf{r}(t)=(x(t),y(t),z(t))\mathbf{r}(t)=(x(t),y(t),z(t))$其曲率为

列出需要的表达式：
$$\{\begin{array}{rl}{x}^{\prime }(\theta )=& -rsin\theta \\ x^{\prime \prime }(\theta )=& -rcos\theta \\ y^{\prime }(\theta )=& rcos\theta \\ y^{\prime \prime }(\theta )=& -rsin\theta \\ z^{\prime }(\theta )=& \frac{h}{2\pi }\\ z^{\prime \prime }(\theta )=& 0\end{array}$$

代入计算
$$\begin{array}{rl}\kappa & =\frac{\sqrt{(z^{\prime \prime }{y}^{\prime }-y^{\prime \prime }{z}^{\prime }{)}^2+(x^{\prime \prime }{z}^{\prime }-z^{\prime \prime }{x}^{\prime }{)}^2+(y^{\prime \prime }{x}^{\prime }-x^{\prime \prime }{y}^{\prime }{)}^2}}{(x^{\prime 2}+y^{\prime 2}+z^{\prime 2}{)}^{\frac{3}{2}}}\\ & =\frac{\sqrt{[0\cdot rcos\theta -(-rsin\theta )\cdot \frac{h}{2\pi }{]}^2+[-rcos\theta \cdot \frac{h}{2\pi }-0\cdot (-rsin\theta ){]}^2+[(-rsin\theta )\cdot (-rsin\theta )-(-rcos\theta )\cdot rcos\theta {]}^2}}{[(-rsin\theta {)}^2+(rcos\theta {)}^2+(\frac{h}{2\pi }{)}^2{]}^{\frac{3}{2}}}\\ & =\frac{\sqrt{(\frac{h}{2\pi }{)}^2{r}^2+r^4}}{((\frac{h}{2\pi }{)}^2+r^2{)}^{\frac{3}{2}}}\\ & =\frac{r\sqrt{(\frac{h}{2\pi }{)}^2+r^2}}{((\frac{h}{2\pi }{)}^2+r^2{)}^{\frac{3}{2}}}\\ & =\frac{r}{(\frac{h}{2\pi }{)}^2+r^2}\end{array}$$

曲率求倒数，所以圆螺旋线曲率半径为

$$\rho =\frac{(\frac{h}{2\pi }{)}^2+r^2}{r}$$

和论文所述一致

参考

维基百科-曲率