基于三坐标测量机的四次EB样条曲线曲面的应用

摘要　文章讨论了四次EB样条的构造、算法及性质，论述了三坐标测量机检测数据的四次EB样条曲线曲面的拟合，运用快速迭代算法建立了被测零件轮廓面的误差评定模型。并给出了应用实例。
关键词　三坐标测量机　EB样条　误差评定

Application of quartic BE-spline based on CMM

Chen Hua, et al

　　在复杂机械零件（如弧面分度凸轮）、模具等带有复杂三维空间曲面零件制造、测量及CAD的造型中，常常需要对三维空间曲面进行描述与控制，从而在计算机屏幕上生成图形或满足NC加工的需要、或将测量的数据拟合成样条曲面（曲线）以满足误差评定的需要。本文讨论的四次EB样条曲线曲面通过引入几何意义十分明显的可调参数，大大提高了曲线曲面的拟合程度。

1　四次EB样条

1.1　四次EB样条曲线
　　四次EB样条曲线Pi(t)是由基本曲线Gi(t)经调配曲线Fi(t)调配而成的。即

Pi（t）=Gi（t）+kFi（t）
(0≤t≤1;i=1,2…，n-3) (1)

　　引入k值可以通过改变它对曲线的形状进行调整，使其与特征多边形的差异可大可小。当k=0时，Pi（t）=Gi（t），四次EB样条曲线就是基本曲线，成为特征顶点的插值曲线。基本曲线段的起始点和终止点通过相应的特征多边形顶点ri+1,ri+2;如果按某种方式拓宽两个控制顶点r0和rn+1,则可使曲线｛Gi（t）｝通过所有的给定顶点r1,r2，…，rn。如果这n个给定点是某条曲线上的型值点，则曲线就是这n个型值点的插值曲线［1］。
　　四次EB样条曲线的矢值方程式为

1.2　四次EB样条曲面
　　四次EB样条曲面可以看作是U、V两个方向EB样条曲线的直积。给定m×n个空间位置向量ri,j(i=1,2,…,n;j=1,2,…,m)构造曲面，其矢值方程式为

1.3　四次EB样条的特性［2］
　　（1）由于在EB样条方程中引入了一个几何意义比较明显的参量K，使设计者能直观地对样条曲线（曲面）的形状进行调整，从而使特征多边形的差异可大可小。
　　（2）四次EB样条可应用给定的一些特殊的特征顶点来生成含有直线和尖点的曲线。
　　（3）四次EB样条曲线段端点处的切矢仅与基本曲线段有关，而与调配曲线段无关。
　　（4）用四次EB样条曲线（曲面）作插值曲线（曲面）时，其设计方法和计算方法均不变，仅需将参量K取为零，这对计算程序的编制极为有用。
　　(5)当K=1/6时，四次EB样条曲线（曲面）又成了B样条曲线（曲面）。由此可见，四次EB样条曲线（曲面）是B样条曲线（曲面）的拓广。

2　测量数据拟合曲面的误差评定

　　在评定三坐标测量机的拟合曲面与实测数据间的误差时，需计算空间实测点到拟合曲面的距离。本文运用了一种快速迭代收敛算法［3，4］求空间任一点到拟合曲面的距离。
　　设P（x，y，z）为曲面外空间一点，求P点到参数曲面∑s（u，w）的距离。其中Qi（ui,wi）为曲面上任意一点，作为P点在曲面上的映射点，是迭代初始点。Qi+1（ui+1,wi+1）为P点在曲面上的垂足。∑s(u,w)为Qi点的切平面，S为PQi+1在切平面上的交点。当Qi与Qi+1重合时,PQi+1为所求点到曲面的距离。Qui,Qwi为Qi在u,w方向上的两个切向量。因此有

　　若PQi+1为点到曲面的距离，则必然S与Qi+1重合，S亦为该点在曲面上的垂足。根据微分几何点到平面的距离为最小距离的条件得

　　搜索方法是，求出（11）式中的Δu,Δw。当（Δu、Δw）＞ε时（ε是任意给定小的正数），用Qi+1代替Qi进行下一轮计算。当（Δu、Δw）＜ε时迭代停止，此时

PQi+1=P-Qi+1=P-Qi=P-S=d

由于规定曲面法矢指向实体之外，所以P在实体之外时d＞0,P在实体之内时d＜0。

3　应用实例

　　弧面分度凸轮机构（见图2）由于结构紧凑,传递扭矩大，特别适用于高速、高精度分度。因其分度运动曲线，可设计成适于各种工作场合的最优运动曲线，从而在运动平稳性、降噪、延长机构使用寿命等方面都具有良好性能。图3是由四次EB样条生成的弧面分度凸轮张量积曲面的轮廓面光照图，其形状为一空间螺旋带。

　　在秦川机床厂和庆安公司的协助下，分别用LK及UKM三坐标测量机，以径向、旋向两种方法对弧面分度凸轮的样件进行了测量。根据以上的测量和误差评定方法，用Visual C++开发了弧面分度凸轮的误差评定软件。取样件设计参数见附表。图4所示为样件在R=57mm的螺旋线上的误差曲线图。该零件设计要求：曲线误差± 0.05mm，轮廓度公差带为0.1mm的双向等距配置。对该零件的评定结果为：凸轮半径57mm的螺旋线轮廓度误差0.147mm，面轮廓度误差 0.165mm。

测量样件设计参数表

旋向曲线类型转位角停歇角分度角中心距分度盘回转半径滚子高度滚子直径

右旋变形正弦 120° 240° 45° 108±0.012mm 49.2mm 19mm ?25-0.005mm

图4　旋向扫描螺旋线误差曲线