基于逐次逼近的垂直度误差评定软件.doc

基于逐次逼近的垂直度误差评定软件林翔（福建商业高等专科学校福建福州350012）【摘要】关于空间直线对基准平面的垂直度误差及评定准则，国际GB/T1958-2004中没有相应的计算方法。本文独辟蹊径，全局搜索包容小圆柱，逐次将包容柱的直径下降，使垂直度误差向“最小区域”逼近。【关键词】垂直度误差，最小区域，高精度SoftwareofPerpendicularityErrorEvaluationBasedontheSuccessiveApproximationLinXiang(FujianCommercialCollege,Fuzhou,Fujian,350012)Abstract:InGB/T1958-2004,thereisnotthecalculationmethodontheperpendicularityerrorofstraightlinetotheplaneandtheevaluationcriterion.Thispaperhastheinnovation,searchesforthesmallouter-containedcylinder,declinesitsdiameterstepbystepandmakesthecalculationoftheperpendicularityerrorapproachtothesmallestarea.Keywords:perpendicularityerror,thesmallestarea,high-precision一、关于垂直度误差按《产品几何量技术规范(GPS)形状和位置公差检测规定》（GB/T1958-2004）规定，垂直度度误差属于位置公差中的定向误差范畴，指实际被测直线相对于给定基准平面在垂直方向上的变动量。直观地说，就是求取一个垂直于基准平面的小圆柱，该小圆柱包容被测直线上的检测点，如果小圆柱的直径达到最小，则其直径就是这条实测直线对于基准平面的垂直度误差值。在精密检测行业中，三维空间直线相对于基准平面的垂直度问题，是计算位置度误差的典型问题。垂直度误差问题的解决，对于解决其他位置度误差计算，具有基础性意义。在精密测量与评定的实践中，往往是当通过三坐标测量机检测获得空间某直线上一系列测量点位置坐标，然后采用满足最小区域判定原则的计算方法，求取这些离散点与基准平面间的垂直度误差。从《GB/T1958-2004》中空间直线对于基准平面的垂直度误差概念的描述可知，如果把直线上各样本点投影到基准平面，此时求取垂直度误差，其实际就是求基准平面各投影点的最小外包圆（或称外接圆），如此，空间的问题就转化为二维平面的问题。亦即对于被测的三维空间“直线”，用三坐标测量机测得一系列点Pk（xk，yk，zk）（k=1~n），则这条姑且称之为“直线”的轮廓线，就以点集Pk（k=1~n）来描述；那么评定“直线”与给定基准平面π0（设方程为：a0*x+b0*y+c0*z+d=0）的垂直度误差f，就转化为评定Pk（k=1~n）相对于π0的垂直度误差。文献[2]规定，用“最小区域法”求出的垂直度误差值，满足“最小条件”原则，计算结果最佳。而对于如何求得符合“最小条件”的垂直度误差，众人见仁见智，推出了各自的算法。经笔者观察分析，提出如下新的算法，其计算过程与结果具有真正的“最小区域”意义，满足“最小条件”的原则。二、新算法设从被测直线上测量得到的样本点Pk（xk，yk，zk）（k=1~n），基准平面π0方程为：a0*x+b0*y+c0*z+d=0，=（）为π0的法向矢量，P’k（x’k，y’k，z’k）（k=1~n）为Pk（k=1~n）在π0上的投影，包容Pk（k=1~n）的小圆柱在平面π0上的投影是一个小圆C。显然，如果小圆柱包容了Pk（xk，yk，zk）（k=1~n），则P’k（k=1~n）必然全部落在小圆C之内。设圆C的圆心坐标为Pq（xq，yq，zq）,那么显然小包容柱的中心线方程即为l1：，δ为圆柱的半径，2*δ即Pk（k=1~n）对于平面π0的空间垂直度误差值f，它是符合文献[2]关于“最小区域”条件下垂直度误差定义的。显然，求Pk（k=1~n）相对于平面π0的垂直度误差的三维空间问题，可以转化成为二维平面的问题，即求包容图一基准平面π0与小圆C关系图P’k（k=1~n）的小圆C的最小直径，如图一中所示。可见，求垂直度误差值的核心问题是找到小圆C，唯一的硬指标是小圆C必须是P’k（k=1~n）的最小外包圆。要