FANUC三点圆分中宏程序：原理、代码与实操避坑

干机加工的老师傅应该都有这种经历：法兰盘、轴承座、圆盘类零件一到手，找圆心就成了最磨人的一道工序。以前车间里最原始的办法是拉表，磁力表座吸在主轴上，百分表顶住内孔壁，摇手轮一点点看表针摆动，碰上直径大一点的孔，来回张罗十几分钟很正常。后来同事在FANUC系统上写了一套三点圆分中宏程序，用三个点的坐标直接算出圆心写进G54，在厂里用了好多年，可靠性是拿批量生产验证过的。今天就把这套宏程序从原理到代码再到实操细节完整整理出来，给经常在法兰克加工中心上干圆类活的朋友做个参考。

这套宏程序解决的痛点很明确：圆孔、外圆、法兰端面这类零件，要的不是“分中”，而是“找圆心”。矩形件用寻边器碰两条边取中点就行，圆的中心靠两条边是碰不出来的。你当然可以拉表，但拉表对操机经验要求高，碰上没加工过的毛坯孔，表针乱摆，很多时候纯粹在碰运气。三点圆分中宏程序的好处是，把找圆心的几何问题转化成三个点的坐标采集问题，只要点取得到、取点位置合理，圆心就是确定的，跟操作者手感没太大关系。

这个内容适合谁看？我总结三类人。一类是在FANUC加工中心上经常干圆类零件的操机师傅，拿过去就能用；一类是做工艺编排的编程员，需要给现场配一套可靠的找正方案；还有一类是在学宏程序的CNC学徒，想找一个真实、能落地、见过生产检验的例子，而不是教材里那种永远跑不起来的演示程序。

1. 为什么需要三点圆分中：从拉表到宏程序的思路转变

1.1 圆类零件找正才是真正的磨人活

很多人刚接触加工中心时觉得分中很简单：装好寻边器，主轴一转，碰左边记X，碰右边记X，取中间，Y方向再来一次。这是矩形工件的玩法。一旦工件是个圆盘、圆环，或者上面有个大孔需要加工，这套思路就卡住了——你找不到“两条边”在哪里。

拉表是圆类零件最传统的找正方法。具体操作是把磁力表座吸在主轴上，百分表头压到孔壁或者外圆表面上，手动转主轴，看表针的跳动量。表针摆动越大，说明主轴回转中心和圆的中心偏离越大。你一边转一边用手轮调整X、Y坐标，直到表针在任意角度位置跳动都在你接受的范围内。这套方法有效，但是有几个问题很现实。

第一是慢。大一点的孔，从粗调到精调，快也得十分钟，慢起来没谱。批量件一天干几十次，人耗在里面，谁干谁知道。第二是挑手感。好多老师傅拉表是“靠手指尖那点感觉”的，新手照着学，按表力度不一样，读出来的数也不一样。第三是毛坯面、半精加工面上拉表，表头打在不规则表面上，跳动的规律很乱，很难判断中心在哪。

三点圆分中宏程序绕开了这些问题。它的逻辑很简单：在圆的任意三个位置分别取一个点的坐标，三点定圆的几何关系是确定的，把坐标喂给程序，圆心就出来了。你不需要让主轴回转中心去对正圆心，也不需要看懂表针的细微摆动规律。对你来说，难度从“手感经验”转移到了“把点取准”，而把点取准这件事，机器本身的坐标精度就能保证大半。

1.2 三点定圆的数学原理与宏程序适配性

三点确定一个圆的几何常识，初中几何就讲过：平面内不在同一直线上的三个点，可以唯一确定一个圆。问题是怎么把这个“唯一确定”变成可计算的公式，并且让人脑能轻松转换成宏程序。

设三个点的坐标分别为P1(X1,Y1)、P2(X2,Y2)、P3(X3,Y3)，圆心为(Xc,Yc)。圆的方程为：

(X - Xc)² + (Y - Yc)² = R²

把三个点的坐标代入，得到三个方程。三个方程两两相减，二次项X²和Y²正好消掉，变成关于Xc和Yc的线性方程组。解这个方程组，整理后就是下面这套公式：

D = X1(Y2-Y3) + X2(Y3-Y1) + X3(Y1-Y2) Xc = [ (X1²+Y1²)(Y2-Y3) + (X2²+Y2²)(Y3-Y1) + (X3²+Y3²)(Y1-Y2) ] / (2D) Yc = [ (X1²+Y1²)(X3-X2) + (X2²+Y2²)(X1-X3) + (X3²+Y3²)(X2-X1) ] / (2D)

之所以把这套公式做成宏程序，是因为它的运算结构非常适合老系统。整套公式只用到了加减乘除、平方和开方，没有三角函数、没有反三角函数、没有坐标旋转，FANUC系统自带的运算功能完全够用，哪怕是十多年前的系统，跑起来也没有任何压力。

这里有个关键点：分母D的几何含义。当三个点在同一条直线上时，D等于0，公式没有意义。换句话说，三个点如果共线，就绝对不能确定一个圆。这个就是宏程序里首个必须做的保护逻辑。实际使用中还有一种更隐蔽的情况，三个点不严格共线，但分布非常差，让D趋近于零，这时候圆心会被放大到很远的地方，程序要能识别出这种异常而不是傻乎乎地输出一个错得离谱的坐标。

1.3 几种分中方式的对比

用表格把这几种方法放到一起看，优点和边界条件会很清楚：

拉表精度的上限其实不低，否则也不会在模具车间存活这么多年。但它的问题是慢和依赖人。当工件尺寸大、批量大、工艺节拍紧张的时候，时间成本就是最大的成本。寻边器四点分中很实用，但只解决矩形或者对称结构的找中心问题。三点圆分中宏程序，恰恰填上了“圆的找正”这个空档，尤其适合没有测头、只有普通寻边器的绝大多数车间。

自动测头的精度和速度无可挑剔，但并不是所有车间都配得起测头，多数中小企业的主力设备还是“三轴+寻边器”的配置。即使有测头的车间，这套三点计算宏程序也有价值——测头采集到坐标之后，最终还是要靠一段计算程序算出圆心，那套计算逻辑和这套手动版的内核完全一致。

2. 宏程序核心代码逐行拆解

2.1 完整程序与变量约定

这套宏程序在车间里验证多年后，固定成了O9010号程序。完整代码如下：

O9010 (3 POINT CIRCLE CENTER) (INPUT: #1=X1 #2=Y1 #3=X2 #4=Y2 #5=X3 #6=Y3) (CALL: G65 P9010 A_ B_ C_ I_ J_ K_) (OUTPUT: #110=CX #111=CY #112=CR) #11=#1*#1+#2*#2 #12=#3*#3+#4*#4 #13=#5*#5+#6*#6 #14=#1*[#4-#6]+#3*[#6-#2]+#5*[#2-#4] IF[ABS[#14]LT0.001]GOTO900 #110=[#11*[#4-#6]+#12*[#6-#2]+#13*[#2-#4]]/[2*#14] #111=[#11*[#5-#3]+#12*[#1-#5]+#13*[#3-#1]]/[2*#14] #112=SQRT[[#1-#110]*[#1-#110]+[#2-#111]*[#2-#111]] #5221=#110 #5222=#111 M99 N900 #3000=1(3P LINEAR) M99

先看输入变量。三个点的坐标通过G65调用传入，使用地址A、B、C、I、J、K。FANUC宏程序的自变量映射关系是固定的：A对应#1，B对应#2，C对应#3，I对应#4，J对应#5，K对应#6。所以在调用的时候要特别注意顺序：

G65 P9010 A[X1] B[Y1] C[X2] I[Y2] J[X3] K[Y3]

也就是A放第一个点的X，B放第一个点的Y，C放第二个点的X，I放第二个点的Y。很多第一次用的人会在这步卡住，习惯性以为I、J、K是X、Y、Z方向上的三个循环量，结果把第二点的Y放错位置，算出来的圆心自然不对。

程序里#11、#12、#13是三个点到原点距离的平方，也就是(X²+Y²)项，在公式里要反复用，提前算出来既节省NC语句条数，也让后面的圆心公式看起来清爽。D值放在#14里，用ABS函数判断是否小于0.001来判断三点是否过于接近共线。

2.2 计算模块背后的逻辑

从#11到#14这四行，实际上就是把前面那套展开公式分步落地。中间计算过程全部使用局部变量#11~#14，避免在宏程序里反复写一长串表达式，万一按错一个括号，排查起来会让人想拆机床。

#110和#111是算出来的圆心X、Y坐标，放在公共变量区里。使用公共变量而不是局部变量是有讲究的：宏程序执行完M99返回后，局部变量#1~#33会被清空，如果你把圆心结果藏在#100以内的某个局部变量里，调用一结束数据就没了。公共变量#100~#199在程序结束后仍然保留，可以供后续程序或者操作者查看。

#112算的是半径。用的第一点到圆心的距离，事实上任意一个点到圆心距离都相等。之所以额外算这个半径，不是为了写坐标系，而是给操作者一个直观的校验参考。手动取点时，如果寻边器碰的是内孔，算出来的#112会比你实际用卡尺量到的孔半径小一个寻边器半径；如果碰的是外圆，会大一个寻边器半径。这个差异不是程序错了，而是测量基准不同，后面实操部分会详细说。

整套计算没有用任何坐标转换指令。三个点在同一个坐标系下，公式直接算出的是该坐标系下的圆心坐标。这意味着只要你取点时读的是机械坐标，算出来的圆心就是机械坐标系里的圆心。反过来，如果你读的是工件坐标，算出来的圆心就是工件坐标系里的相对位置，直接拿去写G54就会出错。

2.3 坐标写入与报警保护

程序末尾用了两个系统变量：

#5221=#110 #5222=#111

在FANUC系统上，#5221和#5222分别对应G54工件坐标系的X偏置和Y偏置。直接把计算出的圆心坐标赋值进去，G54就被设定好了。这条路径在绝大多数FANUC系统上都好用，只要你机床没有改过系统变量的映射关系。

如果你不想赌系统变量地址，或者想把圆心写到G55、G56，可以用G10 L2指令替代。G10 L2是FANUC专门用来写入工件偏置的指令，P代码对应坐标系：

G10 L2 P1 X#110 Y#111 (G54) G10 L2 P2 X#110 Y#111 (G55) G10 L2 P6 X#110 Y#111 (G59)

P1对应G54，P2对应G55，依次到P6对应G59。用G10 L2的好处是语义清楚，不依赖系统变量地址，而且可以一次把X和Y同时写入。缺点是部分老旧系统的G10格式支持程度需要现场确认。我的建议是：车间里统一用哪一种写进去，就在机床说明书上翻到“系统变量”那一章确认一下，两种方式各验证一次，以后就固定下来，别换来换去。

报警保护使用#3000是宏程序的标准做法：

IF[ABS[#14]LT0.001]GOTO900 ... N900 #3000=1(3P LINEAR)

当#14接近0时，程序跳到N900，触发3001号报警，屏幕上显示信息“3P LINEAR”，提示操作者三个点有问题。为什么用#3000报警而不是简单地输出一个提示行？因为宏程序一旦判断数据异常，最优策略就是立刻中止并报警。如果只是打印一段提示文字然后继续往下执行，机床会拿着算出来的离谱坐标写进G54，后面加工直接就废了。报警会强制中断程序，逼着操作者去检查三个点的数据，这个保护逻辑才是整套宏程序能稳定用多年的关键。

3. 实操指南：手动取点、自动测头与精度验证

3.1 手动取点：一套安全高效的流程

手动模式是这套宏程序最经典的使用方式。不需要测头，不需要额外的硬件，一台普通的FANUC加工中心加一个寻边器就能跑起来。完整流程分七步。

第一步，准备工作。把工件压到工作台上，先用手轮把主轴大致移到孔中心附近。装好寻边器，主轴转速调到500~600转/分。转速太低，寻边器接触不稳定，转速太高，碰偏了容易伤到偏心球头。

第二步，取第一个点。用X轴或Y轴方向缓缓接近孔壁，让寻边器的偏心球头轻触内壁。当球头刚发生明显偏摆时，停下来。这时候打开位置显示界面，必须切换到“机械坐标”页面，记下当前的X坐标和Y坐标。这是P1点。

第三步，退开，移动机床，从另一个方向再靠近孔壁，取得P2点。如此再操作一次，取得P3点。三个点的顺序没有强制要求，但覆盖角度最好拉开。

第四步，调用宏程序。在MDI方式下输入：

G65 P9010 A[X1] B[Y1] C[X2] I[Y2] J[X3] K[Y3]

启动运行。程序瞬间算完，G54自动写入圆心。

第五步，检查G54。打开工件坐标系界面，确认G54的X、Y已经是刚刚算出来的数值。这一步别省，程序跑完到确认数据之间，才是最容易被忽略的环节。

第六步，把主轴移到G54的X0Y0，用百分表或者直接用寻边器验证孔中心。如果三个点取得没问题，这个位置应该就是孔的真正中心。

第七步，确认无误后开始加工，或者在批量生产里把G54锁定，后续工件直接沿用。

这套流程最核心的难点其实在于第二步的“记机械坐标”。很多操机新手习惯看屏幕上最大最显眼的那个坐标数值，但那个往往是相对坐标或者工件坐标。三个点必须统一使用机械坐标。如果你把工件坐标当成机械坐标抄进宏程序，而当前G54又不是零偏置，算出来的圆心整个就是错的，而且错得很隐蔽，光看数值很难发现。

3.2 与触发式测头配合的自动取点思路

车间里如果有触发式测头，手动取点这套流程可以大幅升级。自动测头的核心机制是FANUC的G31跳过指令。当测头在移动中碰到工件表面时，跳过信号触发，当前运动被打断，当前位置被记录下来。在宏程序里，当前机械坐标可以通过#5021、#5022、#5023读取，分别对应X、Y、Z轴。

思路是把三点采集也做成自动化。操作者只需要把主轴快速移动到孔内三个方向的预停位置，然后启动测量宏程序，测头会从预停位置慢速进给，依次碰三个点，自动把坐标存入宏程序需要的变量，最后调用同一套圆心计算模块。核心示意如下：

(探测第一个点) G31 G91 X-15. F60 #1=#5021 #2=#5022 (探测第二个点) G31 G91 Y-15. F60 #3=#5021 #4=#5022 (探测第三个点) G31 G91 X15. F60 #5=#5021 #6=#5022

G31进给速度F60是安全起见，用慢速探测能减少测头碰撞瞬间的过冲量。这段代码只是演示采集思路，实际使用时各家测头系统都有自己封装好的探测循环，有的还需要测头回退、信号确认等附加动作。无论测头程序多复杂，最终都要产出一组三个点的XY坐标，给到这套三点圆心计算程序。测量前端可以换，计算内核不需要动，这也是拆开写宏程序的好处。

3.3 取点过程中的误差控制细节

自动测头解决了取点一致性问题，但手动取点才是更多车间的常态，所以误差控制尤其重要。

第一个细节是三个点的分布角度。理想情况是三个点大致均匀分布在圆周上，也就是相互间隔接近120度。间隔越均匀，三点构成的几何约束越强，即使某个点有微小误差，对圆心位置的影响也会被摊薄。如果三个点全部挤在一小段圆弧里，D值变小，同样的单点误差会被放大几倍甚至几十倍。用我车间里师傅的话说，“三点扒在一块儿，圆心都飘到车间外面去了”。

第二个细节是寻边器半径的处理。很多人刚拿到这套宏程序时习惯问：取点之后要不要把寻边器半径减掉？答案是圆心不用减。原因是寻边器球头碰孔壁时，球头中心并不是落在孔壁上，而是落在距离孔壁一个半径的轨迹圆上。这个轨迹圆和实际孔圆是同心的，只是半径小了寻边器半径。三个点整体向圆心缩了一个常数，圆的中心位置不变。程序算出来的#112会接近“孔半径减去寻边器半径”，但这个半径只是参考，不影响G54写入。如果你确实需要程序给出实际孔半径，可以在#112基础上加一个寻边器半径常量，属于简单偏移，但圆心计算不受影响。

第三个细节是工件表面状态。内孔边缘如果有毛刺，或者有切削液带起的铁屑附着，寻边器碰到的是毛刺顶端，坐标自然偏了。手动取点前用油石把孔口毛刺蹭掉，用气枪吹干净，表面擦干净，这三个动作能让分中精度至少稳定一个台阶。

第四个细节是寻边器转速和接触力度。寻边器球头刚接触孔壁时会产生一个肉眼可见的偏摆，这时候就应该停止进给。如果还继续往里推，球头被压死，读数就会偏。手动取点最理想的状态是“刚刚碰偏就停”，这个手感多练几次就能掌握。

3.4 分中完成后如何快速验证

程序把G54写好后，先别急着直接开始铣孔。花一分钟做一次验证，能避免大批量报废。

最简单的验证方法：把主轴手动移动到G54原点，也就是XO和YO处，装上一把百分表，让表头压到内孔壁上，手动慢慢转动主轴，看表针的摆动幅度。如果三点取点质量正常，表针跳动应该在两三个丝以内。如果表针摆动很大，说明某个点取错了或者取点位置太差，回到宏程序的输入数据重新排查。

另一种验证方法是用另一组三个点做一次独立计算。取三个跟之前不同位置的坐标，再调用一次宏程序，算出的圆心如果和第一次结果吻合，说明取点和计算都正常。如果两次算出的圆心相差很大，基本可以断定至少有一组点的取点质量有问题。这种方法不依赖百分表，适合在批量生产中快速抽检。

批量生产的环境里，我更推荐把刀具路径跑一遍首件，用实际加工尺寸反推圆心精度。因为首件加工出来后，孔的坐标位置和理论位置的偏差直接反映了分中结果。首件合格，整套基准就锁定了；后面只要工件装夹方式和毛坯状态不变，G54坐标可以直接沿用。

4. 多年使用磨出来的常见问题与避坑技巧

4.1 三点分布不当导致的计算失真

一个典型的故障现象：宏程序跑完，G54写入的圆心坐标和实际孔中心差了很远，甚至数值大得离谱。第一次遇到这种情况，大多数人会怀疑宏程序公式错了，反复检查代码，最后才发现问题出在三个点的位置选择上。

三点接近共线或分布角度过小时，由于取点本身的微小误差，计算出的圆心会在垂直方向上被急剧放大。可以想象一根很短的圆弧，稍微改变端点位置，拟合出的圆心就会大幅漂移。这种情况在使用手动寻边器时尤其容易发生，因为人对“均匀分布”容易没有感念，随手在相邻的位置取三点，对于一个直径一百多的孔，三个点之间的弧长可能只有十几毫米，算出来的圆心能偏出几毫米。

解决方式很朴实。取第一点之后，把角度大致分成三等分再取后两个点，心里默念一下“隔开点”，或者用分度盘辅助看角度。如果实际工件条件只能取到一段小圆弧，比如孔被其他结构遮挡了大半，三点圆分中就不适合了，这时候老老实实拉表反而更靠谱。这套宏程序的边界条件，一定要在操作规程里写清楚。

4.2 坐标写入失败与系统变量差异

另一种常见情况是：宏程序运行正常，M99正常返回，但是打开坐标系界面一看，G54里面的数值根本没变化。很多人遇到这种情况第一反应是程序没执行成功，反复重新跑，结果还是一样。

这大概率是系统变量地址不匹配。FANUC的标准设置有#5221对应G54的X偏置、#5222对应G54的Y偏置，但不同时期、不同系列的系统，变量分配的细节并不完全一致。部分系统还支持扩展工件坐标系G54.1，其中P1的变量地址又不一样。老系统、新系统、带特殊选项的系统，都可能出现差异。

排查方法很简单。在MDI方式下手动执行一行赋值指令，比如写#5221=50，然后去坐标系界面看G54的X是不是变成了50。如果变了，说明标准变量有效；如果没变，就说明这台机器的变量映射被改过，需要查随机说明书里“系统变量一览表”。更通用的做法是改用G10 L2 P1 X#110 Y#111写入，不依赖系统变量地址，多数情况下能绕开这个坑。

4.3 取点重复性差导致圆心飘

同一个孔，今天分中得出的圆心和明天分中得出的圆心不一样，偏移量有时候超过半毫米。这种情况在手动取点模式下偶尔会发生，而且容易让人误以为宏程序不稳定。

问题几乎都出在取点环节。寻边器接触孔壁时，如果转速不同，球头压紧状态不同，或者孔壁表面附着切削液，每次触发的位置就会产生偏差。还有个容易被忽略的点：手动取点过程中，人移动X轴和Y轴的顺序不同，丝杠反向间隙造成的位置差异也会混进去。

针对反向间隙，有一个小技巧：取每个点时，最后一步进给尽量保持同一个方向接近孔壁。比如三个点都用X轴正向进给的方式去碰，尽量避免有时候正向、有时候反向。这样丝杠间隙对三点的影响是一致的，在后续计算中部分误差会被抵消。当然，最彻底的办法还是上测头，触发式测头的触发力和触发方向一致，重复性好得多。如果车间有测头条件，手动取点主要用于毛坯或者单件，精度要求较高的批量件一律走测头。

4.4 常见报警与故障速查表

把多年使用中遇到的问题集中整理成速查表，现场排查照着这个表来，效率能高很多。

报警3001的处理核心是检查三点坐标输入，尤其要看一下是不是有两个点的坐标几乎一样。车间里出现过把第二点和第三点抄重复的情况，C和J两个地址输入成了同一套数值，程序判定三点共线，直接报警。这类问题纯属操作失误，格式校验和人工核对输入数据可以避免大部分低级错误。

关于G10 L2写入方式和重点。很多操作者第一次见到G10 L2没概念，其实它就是FANUC系统里用程序写入工件坐标系的官方通道。G10 L2 P1的P代码对应坐标系，X和Y地址后面的数值就是要写入的偏置值。用G10写完以后，G54界面立刻就能看到更新的数值，非常直观。我个人的习惯是，正式加工前用G10 L2写完后，再手动去坐标系界面看一眼，确认无误才开跑。

最后聊点个人感受。这套宏程序能稳定用这么多年，靠的不只是那段计算公式，更重要的是每次取点时候的认真劲。宏程序不挑机床，挑的是数据来源，三点坐标有一组不干净，算出来的圆心就一定会偏。新接触这套程序的朋友，建议先在废料上取三点、分中、再拉表验证一轮，把机械坐标读数这个习惯练扎实了再上正经活。程序本身几分钟就能跑熟，真正的功夫，都在取点的那几分钟里。