宏程序铣削圆求

%
O1000
#1=10 (球半径)
#13=20 (毛坯料长度的一半大小)
#2=5 (铣刀半径)
#3=0 (Z方向初抬进刀点)
#4=10 (最终加工深度,与球半径相等)
#10=1.5*#2 (切削步距)
#17=1 (Z方向每层切削量)
#14=0.15 (开粗的切削余量)
G91G28Z0.
G40G80G49
M6T1
S2000M3
G0G90G54X0Y0
G43H01Z[#1+10.]M08
WHILE[#3LE#4] DO1
#11=[#3*[2*#1-#3]]
#5=SQRT[#11] (计算当前位置的球体半径)
#6=#5+#2+#14 (将铣刀半径和余量算进去,得到实际走刀半径)
X#6Y0.
Z[-#3+1.]
G1Z-#3F500.
G2I-#6F1000
#7=#13-#6-#2 (计算同一切削层内XY方向剩余余量)
#8=FIX[#7/#10] (计算剩余走刀次数)
WHILE[#8GE0]DO2
#9=#7-#8*#10
#12=#6+#9 (计算剩余毛料每偏置一步长所得的半径)
/IF[#8EQ0] THEN#12=[#13-#2] (此处加条件是为了让刀具不铣出去,只在毛料内部走刀)
G1X#12F1000.
G2I-#12
#8=#8-1
END2
G0Z[#1+10.]
#3=#3+#17
END1
G91G28Z0.
G28Y0.
M30
%
说明:1、此程序的问题点在于半径是变化的，如何在知道高度的情况下求解得到半径。
2、#11=[#3*[2*#1-#3]] 是由相交弦定理推导过来。（圆内两条相交弦，被分成的两线段之积相等），通过相交弦定理，问题即迎刃而解。
3、当[#13-#2]LT#6时，可跳过该程序段，否则会过切，一般毛料大小的一半与球半径之差至少要大于一个刀具直径才用到。
4、此程序是以最顶面为零点来编程，采用从中心下刀往外开粗加工，有兴趣的朋友可以将其改为从外往里面进刀。
5、对宏程序感兴趣的朋友可以试着用三角函数编写球面精加工。

6、假定毛坯料为40X40X13,用虎钳装夹,留3mm夹位。