An Evaluation Method of Line Profile Errors Based on Least Condition Principle
LU Xincheng1 HUANG Meifa1 TANG Zhemin1 ZHONG Yanru2
1.School of Mechanical and Electrical Engineering，Guilin University of Electronic Technology，
Guilin，Guangxi，541004
2. School of Computer Science and Information Security，Guilin University of Electronic Technology，
Guilin，Guangxi，541004
Abstract： A measuring method of simulating actual measuring tool processes was proposed to evalu⁃ ate the line profile errors. The method followed the least condition principle of shape tolerances in the na⁃ tional standard. The position relationship between the measuring points and the corresponding envelope boundary was analyzed，the set of measurement points was regarded as a rigid body. The relative mo⁃ tions were calculated between the rigid body and the boundary which might occur in the processes of the boundary shrinking to the minimum region，and finally all the measured points were located in the mini⁃ mum containment region. The results show that evaluation processes are more global and directivity than that of the large ⁃ scale search of common optimization algorithm，which may effectively avoid the local solutions of the search results caused by the algorithm defects. The proposed method is suitable for the error evaluation of line profiles.
Key words： line profile；least condition principle；error evaluation；virtual measurement
0 引言在大型船舶、航空和航天等领域，常常需要评定自由曲线的加工质量，因此，线轮廓度误差的准确评定有着重要的意义。由于自由曲线的特殊性，使其有时难以通过简单的几何方程来表示，此外，根据国家标准GB/T 1958⁃2004，形状公差的评定需要满足最小条件。在曲线轮廓度的评定中存在两大难点：一是求算测量点到理想曲线的最
短距离，二是调整测量点与理想曲线的相对位置以满足最小条件。在现有的曲线轮廓度误差评定方法中，大多是将曲线切分为一系列可近似为直线的曲线段，以点到直线的距离来代替点到曲线的距离进行计算 。
［］ 1⁃3求算最短距离的方法主要有分割逼近 、粒子
［］ 4群搜索 、构造近似曲线 等。调整测量点位姿的
［］ 5 ［］ 6优化算法主要有遗传算法 、拟牛顿算法 、禁忌
［］ 3 ［］ 7搜索算法 等。但现有方法仍然存在两个问题： ①
［］ 8
搜索算法效率有待进一步提高； ②搜索算法的全局性难以保证，容易陷入局部解的情况。为解决
上述问题，本文提出一种模拟实际量具测量过程，满足最小条件的曲线轮廓度误差的评定方法。该方法通过分析全局测点特性，快速寻找可能满足最小条件的特征点，然后模拟给定轮廓边界的收缩过程，使所有测量点最终都处于最小包容区域内，最后求算轮廓边界的间距来进行曲线轮廓度误差的评定。
1 线轮廓度公差带的定义与分析
线轮廓度公差用于控制平面曲线轮廓的形状公差，是实际轮廓线所允许的变动全量，线轮廓度公差带是包络一系列圆（直径为公差值t）的两包络曲线之间的区域，各圆圆心应位于理想轮廓曲线上 见图1，其中被测实际轮廓线上的离散测量
［］ 9点必须位于包络一系列圆的两包络曲线之间。
图1 线轮廓度公差带的定义
Fig.1 Defination of line proflie tolerance zone
根据定义，在线轮廓度误差的评定过程中，保持所有测量点之间相对位置不变，即把实际轮廓作刚化处理，刚化的轮廓相对于上下两包络线或包容边界进行旋转和平移，每一步的搜索都要使离散测量点与包络线接触，并使两包络线之间的距离逐渐减小，直至两者间的距离为最小，此时满足了最小条件。
根据刚体运动学原理，对最小区域的搜索实际上是一个有明确指向性的逼近过程，不会陷入局部优化解的情况，理论上具有更高的计算效率。2 理想曲线轮廓数学模型
线轮廓度误差的评定需要构建理想轮廓作为后续评定的参照，理想轮廓可以通过插值拟合来实现。目前，常用的曲线拟合方法有Bezier 法、B
Σ
样条和非均匀有理 B 样条（ NURBS）。由于
i NURBS在IGES、STEP和PHIGS等标准中广泛使用 ，因此本文中的理想轮廓选用NURBS曲
［］ 10
线来建立。
2.1 曲线数学模型的构建
对于理想曲面数学模型的描述通常是通过插值拟合理想点来实现。一条p次的NURBS曲线可以表示为一分段有理多项式矢函数 ：
［］ 10 其中， Pi为曲线上的控制点，这些控制点形成控制多边形来约束曲线的形状，见图2； ωi为权因子，首末权因子ω0、ωn > 0，其余ωi ≥ 0，且按顺序排列的p个权因子不同时为0，以防止分母为零、保留凸包性及曲线不因权因子而退化为一点 ； Ni, ( u ) ［］ 11
p表示定义在非周期且非均匀节点矢量U上的p次B样条基函数，函数形式及计算方法依据de Boor⁃ Cox给出的基函数递推式［］ 10 ：
图2 控制点约束曲线形状（ p = 3）
Fig.2 Control point constraint curve shape（ p = 3）
计算过程中，上式可能会出现0/0，规定0/0= 0。节点矢量的表达式如下：
U ={ a, … ,, a up 1, up 2, …, um 1, b, … , b }
+ + -- p p +1 p +1
式中, m为控制点数。
通常将两端节点的重复度设为p+ 1。参数u的取值范围为［ a，b ］，在多数实际应用过程中设定a= 0， b= 1。
令 2.2 确定NURBS曲线的阶次
NURBS曲线由一系列基函数组合而成，为了保证曲线各连接处具有足够的光滑度，参数曲线需要满足一定的连续性，以符号Gn表示。一般情况下，拟合得到的曲线至少要满足G2连续性［］ 11 ，即
曲线连接处存在关于弧长的二阶导矢，使得该点二阶可微，根据n=p- 1可知NURBS曲线为3次曲线。
2.3 NURBS曲线的局部修改性
考虑到NURBS曲线存在局部修改性，即移动控制点Pi或改变权因子ωi ，仅会影响到对应区间u∈ ［ ui，ui+ ］内的曲线形状。为简便起见，可令
p+ 1
对所有的i，均满足ωi = 1，则对所有的i，均满足Ri， p（ u）=Ni， u），即完全通过控制点来实现曲线
p（
形状的调整。当节点矢量处于不同区间时，有理基函数Ri，（u）
p有如下形式： 其中， Bi， p（ u）表示p次的伯恩斯坦基函数，其二项式递推式如下： 3 模拟接触法计算最小距离
在公差建模领域中，常用到小位移旋量（ small displacement torsor，SDT）理论的观点，该理论存在两个前提假设：刚体假设和小变动假设。假设零件的刚性足够大，其表面形变远小于公差值，可忽略不计，且零件的公差值相对于名义尺寸是极小量 ，因此在评定过程中，测点只会在极小
［］ 12的范围内变动，该范围内的两包络线段的曲率半径极大，曲率值极小，可看作相互平行的两条直线段，见图3。此外，对于平行的包络线段，在向理想曲线收缩的过程中，可能产生的弯曲变形可忽略不计，认为其始终保持平行并同时平行于对应的理想曲线段。
图3 测点运动范围内近似直线的包络线
Fig.3 The envelopes of approximate straight lines in
the range of motion of measuring points
假设测量区域内的三点Qi、Qi+ 1和 Qi+ 2，虚线框表示测点附近的微小变动区域，则区域内的包络边界Li、Li+ 1和Li+ 2可看做直线（图3），上下边界逐渐向理想曲线收缩，上边界与Qi+ 2接触，此时接触点在边界的推动下沿边界Li+ 2的法向矢量Cn （ ui+ 2）移动。由于当前仅有一个接触点，可令其余测点均沿C （ ui+ 2）平动。
当边界收缩一段距离后使得Qi+ 1与下边界接触（图4），此时的上下包络边界的间距小于图3的上下包络边界的间距。接触点Qi+ 1在下边界的推动下沿 Li+ 1的法向矢量 Cn（ ui+ 1）运动，由于 Cn （ ui+ 1）和 Cn（ ui+ 2）的方向不同，因此所有测量点将以 Qi+ 1和 Qi+ 2连线上的一点为旋转中心进行转动。
图4 测点在边界推动下的转动
Fig.4 The rotation of measuring points driven
by the boundary
当旋转至Qi与上边界接触时（图5），3个测量点之间的位置关系满足国家标准GB / T 1958 ⁃ 2004中描述的终止条件，测量点相对于包络边界无法进一步运动，包容区域处于最小区域，两包络线的间距即为线轮廓度误差。
图5 包容边界达到最小区域
Fig.5 The containment boundary reaches the
minimum area
4 最小距离求算过程
通过式（ 1）可求得理想曲线的表达式，对式（ 3）求一阶导数可得到节点u处的曲线切向矢量： 根据直线垂直的特性，可求出对应节点的法向矢量Cn：
Cn ( u· ) C' ( u )=0 （ 5）基于上述前提与说明，现构造测量点附近平行包络线段，通过寻找线段间最短距离来满足最小区域条件并进行线轮廓度误差求算，具体步骤如下：
（ 1）假设给定的公差值为a，设计曲线的节点矢量为U，测量点的坐标为Qi（xi ， yi），通过式（ 1）求出设计曲线上点的坐标C( ( x( , y( )， k表示设
k ) k ) k ) i i i计曲线上的点。对节点矢量U中每一个参数ui计算点距：
d = (- xi x( ) +( yi - y( )
( k ) k ) 2 k ) 2 i i i
若d( ∈ ［ a，λa ］，则继续下一点；反之则沿曲线方向
k ) i
改变( x( , y( )直到满足di∈ ［ a，λa ］，见图6。其中
k ) k ) i i λ的取值范围可根据经验在（ 1，1.5 ］范围内选取。
图6 调整理想点与测量点间距
Fig.6 Adjusting the distance between ideal points and
measuring points
（ 2）结合式（ 4）计算各等分点处切向矢量，根据矢量C' ( u )与坐标C( k ) ( x( k ) , y( k ) )计算切线方程， i i i再由点到直线距离方程，计算各测点Qi与对应切线的距离Di：
式中， ABC、、 为切线方程的参数。
找出其中的最大值Di， 。
max
通过式（ 5）求出对应法向矢量C （ ui），曲线上各点C( k )沿各自法向矢量方向计算距离为Di， 的i max新坐标点，且过新坐标点作平行于切线的直线，见图7。此时， Di， 对应的测点已在包络线段上，包
max
络边界L1接触到第一个测点Q1。
图7 构造平行线
Fig.7 Tectonic parallel line
（ 3）当接触第一个测点后， L1沿理想点的反向法向矢量方向向理想曲线靠近， Q1在L1的推动下沿着相同的矢量方向移动。为了保持所有测点的相对位置不变，令所有测点均沿Q1的移动方向改变位置。对于除L1外的其余包络线段，沿各自反向法向矢量方向向理想曲线靠近。在靠近的过程中，所有包络线段的移动速度v均相同。
由于除Q1外其余测点与各自包络线段的移动方向不同，故会出现图8中的三种情况。由于移动速度v相同，因此可能出现离测点较远的包络线段比离测点近的包络线段先接触测点的情况，故此处不能依据距离关系寻找第二接触点Q2。
测点运动方向的方向矢量为Cn （ u2），包络线段运动方向的方向矢量为Cn （ ui），通过向量夹角公式
（ a） θ< 90° (b) θ= 90° （ c） θ> 90°
图8 向量夹角的三种情况
Fig.8 Three cases of vector angle
① θ∈ ［ 0o，90o），计算测点到包络线段的距离S，则两者接触时所需时间
② θ= 90o，两者接触所需时间
③ θ∈（ 90o，180o ］，两者接触所需时间
计算出所有时间T后， Tmin对应的测点将作为第二
接触点Q2。 n
（ 4）找到Q2后，计算Q2接触L2所需的移动距离。以θ∈（ 90o，180o ］为例，构建Q2与L2的几何关系见图9。Q2向A移动， L2向Q2移动， x为Q2和L2的移动距离， h为Q2到A的距离， S为Q2到L2的距离，由相关几何条件可知：
即可求得移动距离x。
图9 Q2向L2平移
Fig.9 Q2 translation to L2
（ 5）当2个测点均接触到对应包络线段后，第三接触点Q3的搜索是平移和旋转的复合运动。根据国家标准GB/T 1958⁃2004中关于形状公差及其评定的描述，在最小包容区域的判定中有3种准则，基于此可得到间距最小的两平行平面。同理，在线轮廓度误差评定中可以用此判断最小区域。
Q3的复合运动可将其拆分为平移和旋转2个独立运动依次进行。在Q2接触到L2后， Q1、L1和Q2、L2存在如下2种位置情况（图10）：
①当处于图10a的位置关系时，仍可按照搜索Q2的方式搜索Q3，即只存在平移运动。可以证明，
随着L1和L2沿矢量方向平移相同距离后， Q1和Q2能在保持相对间距不变的情况下，平移至处于新位置的两包络线段上，且两点连线与原位置上两点连线平行。
图10 包络线段与测点的2种位置情况
Fig.10 Two positions of envelopes and
measuring points
②当处于图10b的位置关系时，则会产生旋转运动，旋转中心位于Q1和Q2连线上，见图11。可以证明在Q1和Q2连线上存在一点O使得连线绕O点旋转一定角度α后，其两端点B1、B2与包络线段的间距相等。实际情况下，为满足上述条件， O点会随Q1Q2的旋转而不断改变位置。
图11 旋转运动
Fig.11 Rotational motion
设线段Q1Q2长度为L，旋转中心O点将线段Q1Q2分成两部分LQ1O和LQ2O，对应比例关系为
将线段 Q1Q2长度代入式（ 12）中，可得线段LQ1O和LQ2O的长度
根据图11中几何关系，可以证明△ OQ1B1∽ △ OQ2B2。结合式（ 13），可得线段LQ1 和LQ2 长度：
B1 B2
考虑到△ OQ1B1和△ OQ2B2均为等腰三角形，二者的两底角相等，均为（ 180° - α）/ 2，因此两端点B1、B2到包络线段的距离s1、s2分别为
式（ 16）表明，对于任意α，与 β1、β2存在比例p 使s1=s2，随着α不断变化，旋转中心O点同时会沿线段 Q1Q2移动，当确定一个角度α后可得到唯一对应的p。
此外，在旋转的过程中， β1和 β2也会随α发生改变，其与旋转角α的关系为
式中， βp1、βp2分别为转动后β1和β2对应的夹角。
设（ B x1，y1），B （ x2，y2），则转动中心O点（ x0， y0）的坐标分量为
（ 6）线段Q1Q2旋转方向由L1和L2的平移矢量方向决定。分别计算旋转前与旋转后测点到各自包络线段的距离并更新旋转中心点位置。当距离的最小值为0时， Q3接触对应包络线段L3，此时若满足GB/T 11337⁃2004中的判定要求，则完成最小区域的搜索；若不满足判定要求，则继续旋转直到满足要求。
5 实例分析
本实例利用Hexagon青岛公司生产的三坐标测量机，以接触式测量方法获取零件表面部分区域的测点数据。测量部分为某阶梯轴零件的曲面表面，通过对比仪器自带评定软件与本算法的结果，分析本算法的可行性。表1给出了实际测量点数据。
根据测量数据构造NURBS曲线，设曲线次数p= 3，节点矢量U=（ 0，0，0，0，0.340 5，0.436 7， 0.534 6，0.640 6，0.749 6，0.870 0，1，1，1，1，1），可单独通过改变控制点或权因子来调整曲线形状，因此可设ωi= 1。
（ 1）调整理想点与测点间距。计算表1中的测点与理想点的间距，若理想点与测点的间距大于1.5倍公差值，则在理想曲线上选取新的理想点，
最终使理想点与测点的间距小于1.5倍公差值。完成调整后构造理想点处曲线的切线，并计算各测点到对应切线的距离。
（ 2）计算切线方程并搜索最大距离。由理想点的坐标计算得到理想点处的切线方程，再结合测量点数据搜索测点到各自切线的距离Di，其中最大距离Di， max= 0.586 8 mm，对应测点序号为18，见表2。 （ 3）包络边界接触第一接触点。对所有切线沿各自法向矢量方向，在距离测点0.586 8 mm处构造包络线段。此时第一接触点与对应包络线段接触，且矢量模长均为0.586 8 mm。
（ 4）搜索第二接触点。根据式（ 7）计算各测点法向矢量方向与18号测点法向矢量方向的夹角θ，见表3。设所有包络线段沿各自法向矢量方向的移动速度v= 1 mm/s，根据θ计算各测点在18 号测点法向矢量方向上的分速度，并通过该速度计算测点移动至与对应包络线段接触的过程中，所需时间T的最短时间来确定第二接触点。
（ 5）搜索第三接触点。由式（ 11）计算出移动距离x，将测点沿第一接触点法向矢量Cn （ u18）方向平移距离D18，包络线段沿各自法向矢量方向平移
相同距离，完成第二接触点的接触。
此时两个接触点满足图10b中的情况，需通过旋转来搜索第三接触点。反复循环式（ 12）~式（ 17）来调整测点位置，并设定s1（或s2）旋转变换前后的误差ε< 0.01时，停止计算。
使用测量仪器自带的线轮廓度误差评定方法所给出的曲线轮廓度误差评定结果为0.484 mm，本文给出算法的评定结果为0.438 2 mm，这表明本算法是有效的，也是正确的。
6 结论
（ 1）根据轮廓度公差的定义，分析了轮廓度误差的评定过程及相关特征，并建立了模拟实际测量的虚拟量具模型。在遵循物理规律的情况下，模拟实际量具的评定过程，使评定过程中各阶段得到的误差值均比前一阶段误差值小，不存在上下浮动情况，且评定过程中始终考虑所有测量点的整体运动，不存在陷入传统优化算法中局部解的情况，本文方法具有全局性和单一性。
（ 2）以测量仪器分析数据为对比进行曲线轮廓度误差的评定并给出评定结果，结果表明本算法具备有效性及正确性。
将此方法拓展到复杂曲面的轮廓度误差评定中是今后的工作重点。
参考文献：
［ 1 ］ 廖平. 分割逼近法快速求解点到复杂平面曲线最小距离［］ J. 计算机工程与应用， 2009，45（10）：163⁃164. LIAO Ping. Fast Calculating Minimum Distance Be⁃ tween Point and Complex Curve with subdivision Approximating Algorithm ［］ J . Computer Engineer⁃ ing and Applications，2009，45（10）：163⁃164.
［ 2 ］ 温秀兰，赵艺兵，王东霞，等. 改进遗传算法与拟随机序列结合评定自由曲线轮廓度误差［ J. ］ 光学精密工程， 2012，20（4）：836⁃842.
WEN Xiulan，ZHAO Yibing，WANG Dongxia，et al. Evaluating Freeform Curve Profile Error Based
on improved Genetic Algorithm and Quasi Random Sequence ［］ J . Optics and Precision Engineering， 2012，20（4）：836⁃842.
［ 3 ］ 廖平. 基于遗传算法和分割逼近法精确计算复杂曲面轮廓度误差［］机械工程学报， J. 2010，46（10）：2⁃7. LIAO Ping. Calculating of Complex Surface Profile Error Based on Subdivision Approach Algorithm and Genetic Algorithm ［］ J . Journal of Mechanical Engineering，2010，46（10）：2⁃7.
［ 4 ］ 王宇春，孙和义，唐文彦，等. 最小条件下一般二次曲面轮廓度误差的评定［ J. ］ 仪器仪表学报， 2014， 35（8）：1804⁃1809.
WANG Yuchun，SUN Heyi，TANG Wenyan，et al. Evaluating General Quadric Profile Error Based on Least Condition Principle ［］ J . Chinese Journal of Scientific Instrument，2014，35（8）：1804⁃1809.
［ 5 ］ 王东霞，温秀兰，赵艺兵. 基于CAD模型引导测量的自由曲面定位及轮廓度误差评定［ J. ］ 光学精密工程， 2012，20（12）：2721⁃2727.
WANG Dongxia ， WEN Xiulan ， ZHAO Yibing. Localization and Profile Error Evaluation of Free⁃ form Surface Based on CAD Model ⁃ directed Mea⁃ surment ［］ J . Optics and Precision Engineering， 2012，20（12）：2721⁃2727.
［ 6 ］ 郭慧，潘家祯. 基于微粒群算法的叶片曲面形状误差评定［ J. ］东华大学学报（自然科学版），2008，34 （ 5）：769⁃772.
GUO Hui，PAN Jiazhen. Form Error Evaluation of Cared Blade Surface Based on Particle Swarm Opti⁃ mization ［］ J . Journal of Donghua University（Natural Science），2008，34（5）：769⁃772.
［ 7 ］ 刘玉君，朱秀莉，纪卓尚，等. 船体外板曲面形状误差评定方法分析［ J. ］ 哈尔滨工程大学学报， 2006， 27（5）：636⁃638.
LIU Yujun，ZHU Xiuli，JI Zhuoshang，et al. Detec⁃ tion of Profile Error of Ship Hull Plate Processing Surface ［］ J . Journal of Harbin Engineering Universi⁃ ty，2006，27（5）：636⁃638.
［］ 8 BYRD R H，LU P，NOCEDAL J，et al. A Limited Memory Algorithm for Bound Constrained Optimi⁃ zation ［］ J . Siam Journal on Scientific Computing， 1995，16（5）：1190⁃1208.
［ 9 ］ 李柱，徐振高，蒋向前.互换性与测量技术——几何产品技术规范与认证GPS ［ M ］ . 北京：高等教育出版社， 2015：170⁃171.
LI Zhu，XU Zhengao，JIANG Xiangqian. Geometri⁃ cal Product Specifications and Verification ［］ M .
Beijing：Higher Education Press，2015：170⁃171. ［］ 10 PIEGL L，TILLER W. 非均匀有理B样条［ M ］ .2版 .赵罡，穆国旺，王拉柱，译. 北京：清华大学出版社， 2010：46⁃47.
PIEGL L， TILLER W. The NURBS Book ［］ M . 2nd ed. ZHAO Gang， MU Guowang， WANG Lazhu，Trans. Beijing：Tsinghua Universi⁃ ty Press，2010：46⁃47
［ 11 ］ 施法中.计算机辅助几何设计与非均匀有理B样条［ M. ］北京：高等教育出版社， 2001：310⁃312.
SHI Fazhong. CAGD & NURBS ［］ M . Beijing： Higher Education Press，2001：310⁃312.
［ 12 ］ 茅健，张恒，曹衍龙. 基于SDT的平面尺寸公差模型研究及应用［ J. ］ 工程设计学报， 2010，17（5）： 322⁃324.
MAO Jian，ZHANG Heng，CAO Yanlong. Re⁃ search and Application on Model of Size Tolerance for Plane Based on SDT ［］ J . Journal of Engineer⁃ ing Design，2010，17（5）：322⁃324.
(编辑 胡佳慧)
作者简介：陆辛成，男， 1991年生，硕士研究生。研究方向为新一代GPS数字化评定。E ⁃ mail：1248895609@qq. com。黄美发（通信作者），男， 1962年生，教授、博士研究生导师。研究方向为机电系统精度设计和智能测量方法。E⁃mail：hmhmf@guet.edu.cn.
（上接第2319页）
［］ 18 JOE H. Multivariate Models and Multivariate De⁃ pendence Concepts ［］ M . Boca Raton：CRC Press， 1997.
［］ 19 LI Dianqing，TANG Xiaosong，PHOON K K，et al. Bivariate Simulation Using Copula and Its Ap⁃ plication to Probabilistic Pile Settlement Analysis ［］ J . International Journal for Numerical and Ana⁃ lytical Methods in Geomechanics, 2013, 37（6）： 597⁃617.
［］ 20 TU J，CHOI K K, PARK Y H. A New Study on Reliability ⁃ based Design Optimization ［］ J . Journal of Mechanical Design，1999，121（4）：557⁃564.
（编辑 袁兴玲）
作者简介：王倩蓉，女， 1994年生，硕士研究生。研究方向为可靠性优化设计。姜潮（通信作者），男， 1978年生，教授、博士研究生导师。研究方向为复杂装备不确定性分析与可靠性设计理论和汽车CAE技术。获国家科技进步二等奖1项和省部级科技进步一等奖2项，教育部自然科学一等奖1项。E⁃mail：jiangc@hnu.edu.cn。