任意位置弦齿厚计算程序（非固定弦）

发一个以前写的任意位置弦齿厚计算程序（非固定弦）。
这是和他人探讨问题时做的，界面十分简易，大家多包涵。
也许有人用得着。
  

  
M法向模数
alpha n法向压力角
beta螺旋角
z齿数
xn法向变位系数
da外径
hc要计算的位置的齿高，也就是要计算的位置到外径的距离

sc要计算的位置的弦齿厚

可以更改模数啊，只不过更改之后必须输入有效的外径，否则就有可能出现计算错误。这是个简易程序，没有作任何错误检查，所以需要用户保证输入有效的数据。

如果你算出了任意位置（如径=Dx）的圆弧齿厚，那么你算出该状态的弦齿厚时，
注意，你实际算出的并不是齿高=（外径-Dx）/2位置的弦齿厚，而是比这个位置靠下一点。

我的程序的算法是这样的，在齿高处做一平面，该平面将和左右两侧齿面各有一条交线，我计算的就是这两条曲线的最小距离。其实我的算法是适用于斜齿轮也适用于任意齿形的任意位置弦齿厚的。主楼的程序只是渐开线情况下的特例。

下了，感觉不错，直齿螺旋角可以给定一个不为零的数，比如0.00000001。外圆尺寸给定是对的，可以认为实测外圆。很好

计算斜齿可以，结果正确。计算直齿，好像只能取0.0....?替代啦，螺旋角β(beta)取0，会出错！

太好了~~太谢谢了

根据齿轮手册90版上的特例进行了对比，发现和我计算的在4位以后就出现了不同，欢迎探讨一下

这是很久以前我在http://bbs.cmiw.cn/viewthforum.php?mod=viewthread&tid=12207与阿松讨论弦齿厚时所做的。
最终我和他的结果在小数点后14位才开始有不同。
他的方法是根据立体几何推导出了公式。
我的方法是基于解析几何，先写出渐开面和平面方程，求出交线，再求导数找到最短距离。

继续用上面的例子
齿数15，模数1，压力角20，螺旋角0，变位系数为0，齿顶圆直径为17
则图示中的弦齿厚为d*sin(π/2*z)=15*sin(π/30)=1.567926949014

看看是不是上面这个原因。