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

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

2、idealplanearecompared,thelinebetweenthetwovaluesofdistance,wliichistlieflatnesserrorvalues,orbymeasuringtlieactualsurfaceonseveralpointsofrelativeheightdifference,conversiontolinevaluerepresentationofflatnessenorvalue.Tliispaperstudiesthemathenaticalmodelofflatnessenorandidealplanemadebyleastsquarem

3、ethod.Withilhistrationofpracticalcases,theauthorreachestheconclusionthatleastsauaremethodisthebestoneinjudgingandmeasuringlargerplane'sflatnessenor.Keywords:Leastaqiiaremethod;Flatnesserror：Bestmethod目录第一章平面度的测量方法11.1引言11. 2平面度误差的测量1第二章平面度的评定32. 1最小区域法32. 2最小二乘法测量平面度误差的原理32. 2.1建立被测实际表面的数学模型4第三章

4、用MATLA实现的过程63.1软件编程63. 2平面度误差的最小二乘法评定及其评定结果的不确定度93. 3小结11第一章平面度的测量方法1.1 引言平面度误差是指被测实际表面对其理想平面的变动量。理想平面是评定平面度误差的评定基准，而评定基准的方位不同求得的平面度误差值也不同。若用水平仪、U准直仪按节距法测量实际表面上各点相对于测量基准的平面度误差时，确定评定基准的方法可用：简便法、最小二乘法、和最小包容区域法闺。本文着重分析介绍最小二乘法来确定评定基准，从而求得平面度误差值。最小二乘法能准确而充分地利用全部的原始观测数据提供的信息，比较客观地评定出不需要经过多次试算的平面度误差,而且可直接运

5、用于电子计算机运算，列入MATLAB软件，是的平面度误差的计算达到迅速、准确、可靠"L1.2 平面度误差的测量平面度误差的测量方法很多，常用的有如下所列的方法：光隙法：将被测直线和测量基线间形成的光隙，与标准光隙比较，测量:不同方向的若干个截面上的直线度误差，取其中最大值，作为平面度误差近似值的方法。该方法适用于磨削或研磨加工的小平面的平面度误差测量。指示器法：将被测零件支撑在平板上，平板工作面为测曷基准，按一定的方式布点，如图3所示，用指示器对被测面上各点进行测量并记录所测数据，然后，按一定的方法评定其误差值。该方法适用于中、小平面的平面度误差测量。光轴法：以几何光轴建立测量基面，

6、测出被测面相对测量基面的偏离量，进而评定平面度误差值的方法，该方法适用于一般精度大平面的平面度误差测量。干涉法：利用光波干涉原理，根据干涉条纹形状、条数，来确定平面度误差值的方法，该方法适用于精研表面的平面度误差测量。三坐标测量机：三坐标测量机是综合利用精密机械、微电子、光栅和激光干涉仪等先进技术的测量仪器，其原理是在三个相互垂直的方向上，有导向机构、测长元件、数显装置，有一个能够放置工件的工作台，测头可以以手动或机动方式,轻快地移动到被测点上，由读数设备和数显装置，把被测点的坐标值显示出来的一种测量设备。除上面介绍的几种方法外，还有液面法、自准宜仪法等，娄底匚贸职业中专学校实验室在教学中，常

7、采用光学合像水平仪测定平板的平面度误差。第二章平面度的评定平面度误差是被测实际表面相对理想平面的变动量。根据GB/T11337-2004平面度误差检测对形状误差的定义，理想平面的位置应符合最小条件，平面度误差大小，等于包容实际表面且距离为两平行平面之间的宽度。平面度误差的评定方法有:最小包容区域法，对角线法，三远点法和最小二乘法。2.1 最小区域法它是两平行包容平面与实际被测要素的接触状态符合平面度误差判别法中某一准则时，此时两平行平面之间的距离，为平面度误差。用最小条件评定平面度误差有3种判断准则：三角形准则，交叉准则，直线准则。平面度误差的最小区域判别方法是：由两个平行平面包容被测实际表面

8、时，至少有3点或4点相接触，接触点的高低分布，有下列3种形式之一者，即符合最小区域。三角形准则：如果被测实际表面上有3个最高（低）点及一个最低（高）点分别与两个包容平而相接触，并且最高（低）点能投影到3个最低（高）点之间，形成最小包容,称为三角形准则。交义准则：两个高点与两个低点某些实际表面具有鞍形的形状特征，其与上下包容面各有两个接触点，若两高点的投影位于两低点连线之两侧，形成最小包容，称为交叉准则。直线准则：如果被测表面上的同一截面内有两个最高（低）点与一个低（高）点分别和两个平行的包容面相接触，并且有一个最低（高）点的投影要落在两高（或两低）点连线之上，此时也形成最小包容，称为直线准则。

9、对角线法：评定基准平面是通过实际平面的一条对角线和平行于另一条对角线的平面，各测点对此平面的偏差中最大正值与最大负值的绝对值之和为被测实际平面的平面度误差值。三点法：评定基准平面是通过被测实际平面上相距较远且不在一条直线上的三点建立的平面。偏离此平面的最大值和最小值的绝对值之和为平面度误差这三种方法虽然都是针对同一被测平面，但由于评定基准平面选取不同，得到的误差值也不同同。2.2 最小二乘法测量平面度误差的原理最小二乘平面是个理想的平面，它便从实际被测轮廓.上各点到该平面的距离的平方和为最小。因此，最小二乘法的目标函数，即各测点到最小二乘基准平面的距离的平方和7(2-1)F(A,B,C)=Xd

10、f1-1其满足最小化时,F（A,B,C）的（A,B,C）即最小二乘基准平面的法向量,maxEin(2-2)即为最小二乘法平面度误差值叫2.2.1 建立被测实际表面的数学模型平而度误差是指被测实际表面不平的程度，而平面在空间直角坐标系中，它的一般方程为;AX+BY+CZ+D=O;假设测得的点的坐标是（0%4）3=1,2,，n）,其中A,B,C为待定参数。根据最小二乘法原理得到最小二乘平面的待定系数A,B,C为n2%z%Zy送ZyJ匕yl2ry-xiysee-1yl2%"%EES5JyleEEnZE/*zlzl技xyEyEEEZ=B力丝力ZSE*2yl>%”eEE%ylnEz/.y

11、ll空力><"ZEZqhrzlZAylZEZc=由于现在知道A,B,C的系数，即确定了最小二乘平面。接下来就是就能知道改平面的上方的点和下方的点到该平面的距离，由于在该平面上下方都有点，所以点的本身就带有矢量性。点到平面的距离公式为(2-6)尸By二C7(1+A2+B2)根据最小二乘法的原理有:=min(2-7)由于观测到的数据位于最小二乘平面（基准平面）的上下两侧，所以有点到平面的距离有最大值和最小值的两个点！即平面度为(2-8)f"KLx-KjL第三章用MATLA实现的过程3.1 软件编程本文是通过MATLAB软件求得最小二乘平面（基准平面）的系数A,B,C

12、,将最小二乘平面（基准平面）拟合在三维坐标中和求出基准平面的平面度误差。根据观测到的数据如下：x(mm)y(mm)z(mm)10.13243.982-0.00319.64144.092一0.00130.44444.217-0.00230.61829.359-0.00217.74129.210-0.0015.45029.068-0.0074.91246.453-0.01215.88746.580-0.00429.74546.740一0.0021030.10415.7460.0031118.72115.614-0.004125.73815.464-0.0101317.63047.032-0.010

13、1429.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.127-0.0042931.Il

14、l13.1760.0013017.42213.018-0.005现在将MATLAB代码输入软件中：functiontestIsqnonlinlllclearall;closeall;%读入数据data=xlsread(，数据.xls，);%3列对应x,y,zxdata=data(：,2);ydata=data(:,3);zdata=data(:,4);%画出原始数据点，红色figure;plot3(xdata,ydata,zdata,'r*');holdon;%初始值abc都为0注解：将数据放入到excel表格中如卜.图所示juy/人/ABCDE1x(nun>y(mn)z

15、(nm)2110.13243.982-0.0033219.64144.092-0.0014330.44444.217-0.0025430.61829.359-0.0026517.74129.210-0.0017A57QA6R-nnn7注解：将保存好的数据文件与MATLAB的代码文件放到同一个文件夹下。然后就可以启动MATLAB软件打开代码文件运行一下，就可以出来结果图啦。祝你好运！x0=000'%调用最小二乘法非线性拟合函数x=lsqnonlin(x)myfun(x,xdata,ydata,zdata),xO)%求平面度误差%平面方程的系数a=x(l);b=x(2);c=x(3);%根

16、据公式求didi=(zdata-a*xdata-b*ydata-c)/sqrt(a2+b2+1);dmax=max(di);dmin=min(di);f=dmax-dmin%就是要求的平面度误差%画拟合的平面xx=1inspace(min(xdata),max(xdata),200);yy=1inspace(min(ydata),max(ydata),200);X,Y=meshgrid(xx,yy);Z=a*X+b*Y+c;surf(X,Y,Z);colormap;hoidon;gridon;%按照这个函数去拟合，使得式(2)成立functionF=myfunfx,xdata,ydata,zd

17、ata)F=(zdata-x(l)*xdata-x(2)*ydata-x(3)/sqrt(1+x(l)2+x(2)2);得到结果为:0.0003-0.0001-0.00580.0097图1最小二乘平面系数与平面度误差x10图2MATLAB最小二乘法平面方程三维图3.2 平面度误差的最小二乘法评定及其评定结果的不确定度测量不确定度是指测量结果变化的不肯定，是表征被测量的真值在某个量值范围的一个估计，是测量结果含有的一个参数，用以表示被测量值的分散性。这种测量不确定度的定义表明，一个完整的测量结果应该包含被测量值的估计与分散性参数两部分也根据测最不确定度定义，在测最实践中如何对测最不确定度进行合理

18、的评定，这是必须解决的基本问题。对于一个实际测量过程，影响测量结果的精度有多方面的因素，因此测量不确定一般包含若干个分量，各确定度分量不论性质如何，皆可用两类方法进行评定，即A类评定与B类评定。其中一些分量由一系列观测数据的统计分析来评定，称为A类评定；另一些分量不是用一系列观测数据的统计分析法，而是基于经验或其他信息所认定的概率分布来评定，称为B类评定。所有的不确定度分量均用标准差表征，它们或是由随机误差而引起的，或是由系统误差而引起的，都对测量结果的分散性产生相应的影响口叫现在将平面方程的系数A,B,C的表达式分别展开nZ4Z%令（3-1）p=X。q=EH，u=X靖,3=2£,t

19、=整理后得到V=nuo)+2ptquq2o)p2nt2(3-2)C=(ua)-t2)SZi+(qt-po)SXiZj+(pt-qu)£力zJ/VA=(qt-p(o)SZi+(nw-q2)SxjZi+(pt-nt)工力zJ/VB=(pt-uq)EZi+(pq-nt)SxiZj4-(nu-p2)EyizJ/V设按最小二乘法评定公式得到最大偏离值点和最小偏离值的坐标分别为，带入公式f=NL1TdiL得到f平面LS=【ZR(C+AxR+ByR)HzL(C+AxL+ByL)令Ax=xr，l,Ay=yR-yL，Az=zr-zL再令J=(qt-p<*)Ax+(pt-uq)Ay/V£

20、2=(几3-q2)»x+(pq-nt)Ay/V£3=(PQ-nt)Ax+(nu-p2)A),/V得至旷平而LS=Az一£1£ZjwXXiZt-£3S%4展开f平面LS=Az-£l+£2yoXo+£3Zi+/o+Zn岛+£2Xl+£3yi+tl+£2Xn+£3yn转化为矩阵形式，令Z=(ZiZ2Z3Zn)-(E1+E2Xo+E3yo)'y=-(E1+E2X1+E3Y1).-(E14-£2Xn+E3yn)-得到f乎就5=口一(J+QXr+£3yR)ZR

21、-1+(£1+QXl+叼丫1吐+3因被测点的X,Y坐标值在平面度的非敏感方向上，其不确定度可以忽略，因此上式为坐标值Z的线性函数，若被测点的测量结果间互相独立，且具有相同的标准差，即有相同的不确定度则可以用不确定度的公式，得到用最小二乘法评定的平面度的不确定度为若各点被测点直接测量结果间有关，且各自的确定度不同，则要计入协方差的影响,为将表达式简化起见，将矩阵增广，即当时，令，当时，这时候最小二乘法评定的平面度误差的不确定度为U平面LS=l-(ErE2XR+E3yR)2+l+(£l+E2XL+E3yL)2+SiLo(£l+E2Xi+E3yi)2)/2U()由于测量

22、了30组数据，所以n=30。由excel表格就可求出，所以P=57L190、q=965.116u=13762.216w=35543.224t=18360.108v=3.86x。由于最小二乘平面的方程系数已经算出来，KPA=0.0003,B=-0.0001,0-0.0058,在利用点到excel计算出最小二乘平面的距离(3-3)A-Ax,-Byx-C7(1+A2+B2)所以能够得到偏离最小二乘平面的最大值和最小值所对应的坐标(17.690,14.664,0.003)此为最大值，(29.966,47.174,-0.007)为最小值，将数据带入到公式中得到的结果不确定度为3.3小结:最小二乘法评定平面度误差从数理统计分析观点来看，能够准确地和充分地处理原始观测数据所提供的信息，比较客观地评定平面度误差，有利于采用电子计算机进行运算，尤其适用于测量较大平面的平面度误差，测量较大平面时，因其布点较多,若不用最