平面度误差的测量及数据处理资料

1、学习资料保密保密学校代码：学 号:课程设计说明书题 目:平面度误差的测量及数据处理 学生姓名：保密学院：保密班级：保密指导教师：保密仅供学习与参考摘要平面度误差是将被测实际表面与理想平面进行比较， 两者之间的线值距离即为平 面度误差值；或通过测量实际表面上若干点的相对高度差，再换算以线值表示的平面 度误差值。本文就平面度误差的数学模型与按最小二乘法建立理想平面（评定基准） 的数学模型展开分析讨论；并结合案列分析，得出比较客观地评定平面度误差或者测 量较大平面的平面度误差，最小二乘法是最佳方法 关键词：最小二乘法；平面度误差；最佳方法AbstractFlatness error is measu

2、red by actual surface with an ideal plane are compared, the line between the two values of distance, which is the flatness error values, or by measuring the actual surface on several points of relative height difference, conversion to line value representation of flatness error value .This paper stu

3、dies the mathenatical model of flatness error and ideal plane made by least square method .With illust ration of practical cases, the author reaches the conclusion that least sauare method is the best one in judging and measuring larger plane' s flatness error.Keywords : Leastaquare method ; Fla

4、tness error； Best method目录第一章 平面度的测量方法 11.1 引言 11.2 平面度误差的测量 1第二章 平面度的评定 321 最小区域法 32.2 最小二乘法测量平面度误差的原理 32.21 建立被测实际表面的数学模型 4第三章 用MATL/实现的过程 63.1 软件编程 63.2 平面度误差的最小二乘法评定及其评定结果的不确定度 93.3 小结 11第一章 平面度的测量方法1.1 引言平面度误差是指被测实际表面对其理想平面的变动量。 理想平面是评定平面度误 差的评定基准， 而评定基准的方位不同求得的平面度误差值也不同。 若用水平仪、 自 准直仪按节距法测量实际表面

5、上各点相对于测量基准的平面度误差时， 确定评定基准 的方法可用：简便法、最小二乘法、和最小包容区域法 2 。本文着重分析介绍最小二 乘法来确定评定基准， 从而求得平面度误差值。 最小二乘法能准确而充分地利用全部 的原始观测数据提供的信息，比较客观地评定出不需要经过多次试算的平面度误差， 而且可直接运用于电子计算机运算，列入 MATLAB件，是的平面度误差的计算达到 迅速、准确、可靠 7 。1.2 平面度误差的测量平面度误差的测量方法很多，常用的有如下所列的方法： 光隙法： 将被测直线和测量基线间形成的光隙， 与标准光隙比较， 测量不同方 向的若干个截面上的直线度误差， 取其中最大值， 作为平面

6、度误差近似值的方法。 该 方法适用于磨削或研磨加工的小平面的平面度误差测量。 指示器法： 将被测零件支撑在平板上， 平板工作面为测量基准， 按一定的方式 布点，如图 3 所示，用指示器对被测面上各点进行测量并记录所测数据，然后，按一 定的方法评定其误差值。该方法适用于中、小平面的平面度误差测量。 光轴法： 以几何光轴建立测量基面， 测出被测面相对测量基面的偏离量， 进而 评定平面度误差值的方法，该方法适用于一般精度大平面的平面度误差测量。 干涉法：利用光波干涉原理，根据干涉条纹形状、条数，来确定平面度误差值 的方法，该方法适用于精研表面的平面度误差测量。 三坐标测量机 ：三坐标测量机是综合利用

7、精密机械、微电子、光栅和激光干涉仪等先进技术的测量仪器，其原理是在三个相互垂直的方向上，有导向机构、 测长元件、数显装置，有一个能够放置工件的工作台，测头可以以手动或机动方式， 轻快地移动到被测点上， 由读数设备和数显装置， 把被测点的坐标值显示出来的一种 测量设备 4 。除上面介绍的几种方法外， 还有液面法、自准直仪法等， 娄底工贸职业中专学校 实验室在教学中，常采用光学合像水平仪测定平板的平面度误差。第二章 平面度的评定平面度误差是被测实际表面相对理想平面的变动量。根据GBT1 1337-2004 平面度误差检测 对形状误差的定义， 理想平面的位置应符合最小条件， 平面度误差大 小，等于包

8、容实际表面且距离为两平行平面之间的宽度。平面度误差的评定方法有： 最小包容区域法，对角线法，三远点法和最小二乘法 5 。2.1 最小区域法它是两平行包容平面与实际被测要素的接触状态符合平面度误差判别法中某一 准则时， 此时两平行平面之问的距离， 为平面度误差。 用最小条件评定平面度误差有 3 种判断准则：三角形准则，交叉准则，直线准则。平面度误差的最小区域判别方法 是：由两个平行平面包容被测实际表面时， 至少有 3点或 4点相接触，接触点的高低 分布，有下列 3 种形式之一者，即符合最小区域。 三角形准则：如果被测实际表面上有3个最高（低）点及一个最低（高）点分别与 两个包容平面相接触，并且最

9、高 （低）点能投影到 3个最低（高）点之间，形成最小包容， 称为三角形准则。 交叉准则： 两个高点与两个低点某些实际表面具有鞍形的形状特征， 其与上下 包容面各有两个接触点， 若两高点的投影位于两低点连线之两侧， 形成最小包容， 称 为交叉准则。 直线准则：如果被测表面上的同一截面内有两个最高（低）点与一个低（高）点分 别和两个平行的包容面相接触，并且有一个最低 （高）点的投影要落在两高 （或两低） 点连线之上，此时也形成最小包容，称为直线准则。 对角线法：评定基准平面是通过实际平面的一条对角线和平行于另一条对角线 的平面，各测点对此平面的偏差中最大正值与最大负值的绝对值之和为被测实际平面 的

10、平面度误差值。 三点法：评定基准平面是通过被测实际平面上相距较远且不在一条直线上的三 点建立的平面。偏离此平面的最大值和最小值的绝对值之和为平面度误差这三种方法虽然都是针对同一被测平面， 但由于评定基准平面选取不同， 得到的 误差值也不同 6 。2.2 最小二乘法测量平面度误差的原理最小二乘平面是个理想的平面，它使从实际被测轮廓上各点到该平面的距离的平 方和为最小。因此，最小二乘法的目标函数，即各测点到最小二乘基准平面的距离的 平方和nF(A,B,C)=i 1d.2IJ(2-1)其满足最小化时，F (A,B,C)的(A,B,C)即最小二乘基准平面的法向量,f dijmaxdumin(2-2)即

11、为最小二乘法平面度误差值2.2.1建立被测实际表面的数学模型平面度误差是指被测实际表面不平的程度，而平面在空间直角坐标系中，它的一 般方程为;AX+BY+CZ+D=0假设测得的点的坐标是(冷,乙)(1=1,2，,n),其中A,B,C为待定参数。根据最小二乘法原理得到最小二乘平面的待定系数A,B,C为B=nA= xiYiZiX| Z|YiXiYi/nxiX|2X|YiXiYi22Y|Z|yIy|XiYiYi(2-3)Z|X|YinX|Yi2/2XiX|X|Z|XIX|xi yi2YiXiYiYi ziyiXi YiYi(2-4)Z|X|YinX|Yi22Xi Z|X|XiYi/XiX|xi Yi

12、22YiZixi YiyiYIXiYiYiC=(2-5)由于现在知道A,B,C的系数，即确定了最小二乘平面。接下来就是就能知道改平面 的上方的点和下方的点到该平面的距离，由于在该平面上下方都有点，所以点的本身就带有矢量性。点到平面的距离公式为dizi Axi By i Ci(1 A2(2-6)根据最小二乘法的原理有:n2dimin由于观测到的数据位于最小二乘平面（基准平面）的上下两侧, 离有最大值和最小值的两个点！即平面度为djmaxdjmin（2-7）所以有点到平面的距（2-8）第三章 用MATL实现的过程3.1 软件编程A,B,C，将最小根据观测到本文是通过MATLAB件求得最小二乘平面（

13、基准平面）的系数二乘平面（基准平面） 拟合在三维坐标中和求出基准平面的平面度误差 的数据如下：x（ mm） y（ mm） z（ mm）110.13243.982-0.003219.64144.092-0.001330.44444.217-0.002430.61829.359-0.002517.74129.210-0.00165.45029.068-0.00774.91246.453-0.012815.88746.580-0.004929.74546.740-0.0021030.10415.7460.0031118.72115.614-0.004125.73815.464-0.0101317.6

14、3047.032-0.0101429.96647.174-0.007155.98944.945-0.0111616.80845.071-0.0051728.76145.209-0.0021828.90133.1250.0001917.10832.988-0.007203.97532.836-0.008214.18814.508-0.0082217.69014.6640.0032330.90914.8170.0042430.14032.118-0.0032518.37831.981-0.007264.37131.820-0.0042718.42131.982-0.0082830.88932.12

15、7-0.0042931.11113.1760.0013017.42213.018-0.005现在将 MATLAB码输入软件中：function testisqnonlin111 clc;clear all;close all;%读入数据data = xlsread(' 数据.xls');% 3列对应x,y,zxdata = data(:,2);ydata = data(:,3);zdata = data(:,4);%画出原始数据点，红色figure;plot3(xdata,ydata,zdata,'r*');hold on;%初始值abc 都为0jyy #11A

16、EcDE1x (nil)y ( mm)E (fflilTi)z110. 13Z43. 35Z一Q QQ33219, 64144b Q9E一Q DOI43 諷2L7-o. 0023430, 618姐353一 Cu 0Q26 1517. T4L29. 210-0. ooi7ft吕 i'RAg hg-h me注解：将数据放入到 excel表格中如下图所示注解：将保存好的数据文件与MATLAB的代码文件放到同一个文件夹下。然后就可以启动 MATLAB软件打开代码文件运行一下，就可以出 来结果图啦。祝你好运！x0 = 0 0 0'%调用最小二乘法非线性拟合函数 x =ls qnon li

17、 n( (x) myfun( x,xdata,ydata,zdata),x0)%求平面度误差%平面方程的系数 a = x(1);b = x(2);c = x( 3);%根据公式求didi = (zdata - a*xdata - b*ydata - c)/sqrt(aA2+bA2+1);dmax = max(di);dmin = mi n(di);f = dmax - dmi n %就是要求的平面度误差%画拟合的平面xx = lin space( min( xdata),max(xdata),200);yy = lin space( min( ydata),max(ydata),200);X,

18、Y = meshgrid(xx,yy);Z = a*X+b*Y+c;surf(X,Y,Z)；colormap;hold on ;grid on;%按照这个函数去拟合，使得式(2)成立function F = myfun(x,xdata,ydata,zdata)F = (zdata - x(1)*xdata - x(2) *ydata - x( 3)/sqrt(1 + x(1)A2+x(2F2);得到结果为：0. 0003-0. 0001-0. OCIoSf 二0. 0097图1最小二乘平面系数与平面度误差111 0图2 MATLAB最小二乘法平面方程三维图3.2平面度误差的最小二乘法评定及其评

19、定结果的不确定度测量不确定度是指测量结果变化的不肯定，是表征被测量的真值在某个量值范围 的一个估计，是测量结果含有的一个参数，用以表示被测量值的分散性。这种测量不 确定度的定义表明，一个完整的测量结果应该包含被测量值的估计与分散性参数两部 分9。根据测量不确定度定义，在测量实践中如何对测量不确定度进行合理的评定，这 是必须解决的基本问题。对于一个实际测量过程，影响测量结果的精度有多方面的因 素，因此测量不确定一般包含若干个分量， 各确定度分量不论性质如何，皆可用两类 方法进行评定，即A类评定与B类评定。其中一些分量由一系列观测数据的统计分 析来评定，称为A类评定；另一些分量不是用一系列观测数据

20、的统计分析法，而是基于经验或其他信息所认定的概率分布来评定，称为B类评定。所有的不确定度分量均 用标准差表征，它们或是由随机误差而引起的，或是由系统误差而引起的，都对测量 结果的分散性产生相应的影响10。现在将平面方程的系数A,B,C的表达式分别展开yiXiYi2yi(3-1)nZi令 V=xixi2yiXi yiP», qQ* 3 = 丫兀心乃莎整理后得到777(3-2)V = nu<i)+ 2ptq - uq - wp nrC = （血-坷+ （妙-卩曲龙巧可+ （Pt"羿）»问力=-呻迟岭+仗少-迟絹+（Pt尬迟炳/V|(pt - uqZi + (pq

21、 -昭 + (nu -卅)工焙/卩设按最小二乘法评定公式得到最大偏离值点和最小偏离值的坐标分别为，带入公式minf dijmax |dj得到平窗L5 = 令 一(f + 点耳 + B%) - zt(c + Axl + Byj 令珀-闯,亠二*-刃,卩二-幼再令可=(仇-叫)亠+ O 7gMj/lz占2 = (w-护)亠 + (pq - nty/v£3=(pq-nt) + (nu - p2)Aj/K得到澤砂7 -叼& -辽呐-吃"F 展开f 平窗圧=AZ- £1 + Wo + 勺可 + 5 + 础R + £2X1 + 勺y 1 + £1

22、+ 臥 + 勺儿转化为矩阵形式，令 pi-" ! 匚-1_ （巧+ E汽+ £讥），7=-（5 +尤1十£"1）4 V-（巧+£的严勺儿）得到 :'_：一 j 一 二一* 厂 丫 一：卜 r因被测点的x,丫坐标值在平面度的非敏感方向上，其不确定度可以忽略，因此上式为 坐标值z的线性函数，若被测点的测量结果间互相独立，且具有相同的标准差，即有相同的不确定度则可以用不确定度的公式，得到用最小二乘法评定的平面度的不确定度为 若各点被测点直接测量结果间有关，且各自的确定度不同，贝U要计入协方差的影响, 为将表达式简化起见，将矩阵增广，即当时，令

23、，当时，这时候最小二乘法评定的平Y1+ 叼可 + f3)2，n面度误差的不确定度为® +切*)十1 + 十叽+训)+ §(巧由于测量了 30组数据，所以n=30。由excel表格就可求出，所以P=571.190、q=965.116 u=13762.216 w=35543.224 t=18360.108 v=3.86x。由于最小二乘平面的方程系数已经算出来，即A=0.0003, B=-0.0001，C=-0.0058,在利用点到excel 计算出最小二乘平面的距离dizi Axi Byi C.(1 A2 B2)(3-3)所以能够得到偏离最小二乘平面的最大值和最小值所对应的坐标(17.690,14.664,0.003)此为最大值，(29.966,47.174, -0.007)为最小值，将数据带 入到公式中得到的结果不确定度为 气=1创m。3.3小结:最小二乘法评定平面度误差从数理统计分析观点来看，能够准确地和充分地处理 原始观测数据所提供的信息，比较客观地评定平面度误差，有利于采用电子计算机进 行运算，尤其适用于测量较大平面的平面度误差，测量较大平面时，因其布点较多， 若不用最小二乘法，则往往需要经过许多次试算，使得计算变得繁琐。而采用最