CN110487210A - 蜂窝芯表面轮廓测量方法 - Google Patents

蜂窝芯表面轮廓测量方法 Download PDF

Info

Publication number
CN110487210A
CN110487210A CN201910829610.5A CN201910829610A CN110487210A CN 110487210 A CN110487210 A CN 110487210A CN 201910829610 A CN201910829610 A CN 201910829610A CN 110487210 A CN110487210 A CN 110487210A
Authority
CN
China
Prior art keywords
honeycomb core
axis
line laser
measurement
data
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201910829610.5A
Other languages
English (en)
Other versions
CN110487210B (zh
Inventor
康仁科
董志刚
秦炎
鲍岩
朱祥龙
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Dalian University of Technology
Original Assignee
Dalian University of Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Dalian University of Technology filed Critical Dalian University of Technology
Priority to CN201910829610.5A priority Critical patent/CN110487210B/zh
Publication of CN110487210A publication Critical patent/CN110487210A/zh
Application granted granted Critical
Publication of CN110487210B publication Critical patent/CN110487210B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B11/00Measuring arrangements characterised by the use of optical techniques
    • G01B11/24Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Length Measuring Devices By Optical Means (AREA)

Abstract

本发明公开了一种蜂窝芯表面轮廓测量方法,包括如下步骤:进行规定路径扫描蜂窝芯;利用扫描后的轮廓数据,分别计算出线激光测头绕Z轴的安装误差角αZ、绕Y轴的安装误差角αy和绕X轴的安装误差角αx;由数控机床带动传感器按照一定路径扫描待测蜂窝芯表面;根据线激光测头的安装误差角和测头的运行坐标,将测量的在线激光坐标系下的测量数据经坐标变换统一到全局坐标系;将各个路径下测量的蜂窝芯数据进行拼接,获得整个待测蜂窝芯的表面数据。该方法有效解决了线激光安装误差角难以确定的问题,方法操作方便、精度高、适用性强、成本低,能够满足实际加工后的蜂窝芯表面的测量要求。

Description

蜂窝芯表面轮廓测量方法
技术领域
本发明涉及线激光测量技术领域,尤其涉及一种蜂窝芯表面轮廓测量方法。
背景技术
蜂窝芯是一种低密度蜂窝状材料,因其较高的比强度、比刚度等优质特性而广泛应用于航空航天、高速列车等领域。但由于其薄壁、弱刚度的特性,其加工面形精度的高精度测量缺乏有效的手段。
申请号为201610585321.1、201610585419.7的发明专利公开了一种用于蜂窝芯面形测量方法,在待测工件表面覆上一层薄膜,以一定的真空度将薄膜吸附压在蜂窝芯表面,通过激光位移传感器测量薄膜表面,再进一步提取出蜂窝芯面形。申请号为201711002057.5,名称为“蜂窝芯面形测量数据毛刺去除方法”、申请号为201711002057.5,名称为“蜂窝芯面形测量数据毛刺去除方法”和201810379213.8,名称为“蜂窝芯面形测量数据毛刺去除方法”的发明专利中均公开了基于激光位移传感器的测量数据的处理方法。而对于实际的蜂窝芯需要测量的尺寸通常较大,如何能够利用激光位移传感器高效精确地测量出蜂窝芯的表面同样是亟待解决的问题。
线激光位移传感器在机器视觉、逆向工程、工业检测、生物医学等领域被广泛应用,其基于激光三角法并能单次测量一条线的轮廓,具有测量精度高、效率高的优点。线激光直接测量的数据是两个维度的,其通过沿另外一个维度的扫描运动补充第三维度信息,实现三维测量。
为了实现大尺寸物体的测量,只需将传感器安装在机床上,由机床的运动带动其实现对物体的扫描。但由于安装误差的存在,传感器坐标系的方向与全局坐标系不能完全平行,线激光直接测量的数据结果与理想结果有一定误差。要想最终获得被测对象的整体面形数据,需要准确获得每次传感器测量数据的全局坐标,才能最终拼接成被测对象的整体面形。测量数据由传感器坐标系变换到全局坐标系,需要确定两个坐标系之间变换所需的平移矩阵和旋转矩阵。平移矩阵能够可以由传感器的运动位置确定,根据机床反馈的坐标很容易获得。而旋转矩阵反应的是传感器坐标系三个轴的方向相对于全局坐标系不平行的偏角误差,根据三个轴的偏角即可确定该旋转矩阵。
申请号为201810269833.6,名称为“采用二维线激光扫描仪获取三维精确数据的方法”的发明专利公开了一种用于根据线激光绕x、y和z轴与理想位置的角度偏差,对线激光测量的数据进行校准的方法。其虽然给出了由这三个误差角确定坐标变换中旋转矩阵的方法,但其并未给出这三个误差角的确定方式。而如何高效高精度地确定出这三个偏角的具体值才是制约线激光实现大尺寸扫描的关键所在。在实际应用时,尤其在一些在线测量的场合,由于现场有切削液或者切屑等会对线激光造成污染,在加工时往往需要拆下测头,在测量时再重新安装。每次线激光重新安装后,其三个误差角需要重新确定。因此,一种易于操作的、高精度的、快速的线激光安装误差角的确定非常关键,是线激光能否通过数控机床实现便携高精度测量、甚至在位测量蜂窝芯面形的关键。
发明内容
本发明针对以上问题的提出,而公开了一种蜂窝芯表面轮廓测量方法。
本发明基于对蜂窝芯本身的测量,确定出线激光绕三个坐标轴的旋转误差角,再通过机械调节或者坐标补偿的方式对测量结果进行校准,获得线激光在三维测量中准确的三维数据。并经过合适的扫描路径对整个蜂窝芯表面进行扫描,经过到全局坐标变换后能实现最终蜂窝芯面形高精度的数据拼接
为达到上述目的,本发明所采取的方法包括如下步骤:
一种蜂窝芯表面轮廓测量方法,包括如下步骤:
步骤1:测量平台搭建
在三轴数控机床的工作台上搭建蜂窝芯的测量平台,其中O-XYZ坐标系为机床运动坐标系;
将蜂窝芯固定于工作台,根据待测蜂窝芯表面对面的面形形状,设计对应支撑结构,使三个方向的蜂窝壁都沿着Z方向;
线激光传感器安装在三轴数控机床,使线激光的入射方向沿着Z方向;线激光沿Y轴或X轴做直线扫描运动,做固定步距扫描测量;线激光测量时记录每条激光线对应的测头运动坐标;
步骤2:确定线激光测头的安装误差角
将线激光测头安装的倾斜误差依次分解为绕Z轴旋转αZ角度,绕Y轴旋转αy角度,绕X轴旋转αx角度;线激光测量后的数据坐标系与机床坐标系存在错位,另建立线激光测量数据坐标系为ol-xlylzl,其中xl轴和zl轴为线激光测量数据本身的坐标,xl轴为激光线方向,zl轴为测量高度值,yl轴为沿机床运动方向为线激光补充的坐标轴;
为确定线激光测头绕Z轴的安装误差角αZ,绕Y轴的安装误差角αy和绕X轴的安装误差角αx,进行规定路径扫描蜂窝芯;
利用扫描后的数据计算出线激光测头绕Z轴的安装误差角αZ、绕Y轴的安装误差角αy和绕X轴的安装误差角αx
步骤3:扫描待测蜂窝芯
线激光测头校准后,由数控机床带动传感器按照一定路径扫描,确保整个蜂窝芯都被扫描到;
步骤4:被测蜂窝芯表面重建
根据线激光测头的安装误差角和所记录的机床的运行坐标,确定出坐标变换矩阵,将测量的在线激光坐标系下的测量数据经坐标变换统一到全局坐标系;坐标变换关系为:
P=Rl×Pl+T (9)
其中P(X,Y,Z)为测量后的数据在全局坐标系下对应的坐标,Pl(xl,0,zl)为线激光测量的数据,旋转变换矩阵Rl为:
Rl=Rxx)×Ryy)×RzZ) (10)
Rx,Ry,Rz分别为αx,αy,αZ三个误差角所对应的旋转变换矩阵;T=[x0,y0,z0]T为测量数据对应的传感器测头在机床中的坐标;
对全局坐标系下的测量数据进行蜂窝芯表面的重建。
进一步地,所述线激光沿Y轴做直线运动,做固定步距扫描测量的实现方法为:
利用Y轴电机的编码器脉冲信号或光栅尺反馈信号作为线激光测量的触发信号,编码器脉冲的间隔即为线激光测量的步距Δy
将机床的三个轴的编码器脉冲信号或者是光栅尺反馈信号,接入脉冲计数处理***,在Y轴电机的编码器脉冲信号或光栅尺反馈信号触发线激光测量时,同时触发该脉冲计数处理***记录下三个轴的脉冲数,根据脉冲数与实际运动的距离的关系,能够确定测头运动坐标[x0,y0,z0]。
进一步地,所述为确定线激光测头的3个安装误差角,进行扫描蜂窝芯的规定路径为:
线激光三次扫描蜂窝芯:
线激光第一次扫描蜂窝芯,线激光沿Y轴运动,测量蜂窝芯的轮廓数据;
线激光第二次扫描蜂窝芯,线激光沿X轴正向平移距离w,以与第一次扫描时相同的Y向路径再次扫描蜂窝芯,测量蜂窝芯的轮廓数据,其中,w小于线激光的扫描宽度w0
线激光第三次扫描蜂窝芯,线激光沿Z轴正向抬高距离h,以与第二次扫描时相同的Y向路径再次扫描蜂窝芯,测量蜂窝芯的轮廓数据,其中,h小于线激光的量程;
所述蜂窝芯的形状包括如相互连接的六边形,圆形,多边形的规则形状或非规则形状。
进一步地,所述线激光三次扫描蜂窝芯时:
第二次相对第一次扫描测头沿X轴正向平移距离w,沿Y向扫描路径相同其中,平移距离w满足:
w=β×w0 (1)
β取值为0.5~0.75;
两次测量的数据中沿X轴方向各有(1-β)w0宽度范围内的测量数据,对应相同的蜂窝芯轮廓;
扫描前,调节线激光到蜂窝芯的距离,使三次扫描时蜂窝芯都在线激光测量量程范围内。
进一步地,所述计算线激光测头绕Z轴的安装误差角αZ的具体过程为:
在xlolyl平面内比较第一次和第二次测量的相同蜂窝芯轮廓数据,第二次测量相对于第一次测量,蜂窝芯轮廓沿yl轴方向的偏移量为lz,线激光测头绕Z轴的安装误差角满足:
tanαZ=lz/w (2)
αZ=arctan(lz/w) (3)
其中lz确定的具体过程为:
将第二次测量的数据沿yl轴正方向和负方向依次平移i行(i=-nmax,-nmax+1,…,-1,0,1,…nmax-1,nmax),其中nmax对应移动的数据行数范围,取值应大于预估的偏移量;计算平移后两次测量重合区域测量数据的重合数据量(或重合的面积),找到最大重合数据量(或重合的面积)对应的平移行数n1
偏移量为lz为:
lz=n1×Δy (4)。
进一步地,所述计算线激光测头绕Y轴的安装误差角αy的具体过程为:
在xlolzl平面内比较第一次和第二次测量的相同蜂窝芯轮廓数据,第二次测量相对于第一次测量,蜂窝芯轮廓沿zl轴方向的偏移量ly,线激光测头绕Y轴的安装误差角满足:
tanαy=ly/w (5)
αy=arctan(ly/w) (6)
其中ly确定的具体过程为:
将第二次测量的数据平移w个距离后,重复测量部分蜂窝芯的轮廓重合,计算重合的数据第二次测量相对于第一次测量的差值,ly取为所有差值的平均值。
进一步地,所述计算线激光测头绕X轴的安装误差角αx的具体过程为:
比较第二次和第三次测量的相同蜂窝芯轮廓数据,第三次测量相对于第二次测量,蜂窝芯轮廓沿zl轴方向的偏移量h',线激光测头绕Z轴的安装误差角αx满足:
cosαx=h/h' (7)
αx=arccos(h/h') (8)
其中h'确定的具体过程为:
第三次相对于第二次测量的在xlolyl平面内蜂窝芯的轮廓基本重合,计算重合的数据第三次测量相对于第二次测量的差值,h'取为所有差值的平均值。
进一步地,所述扫描待测蜂窝芯的具体过程为:
由数控机床带动传感器在YOZ平面内做扫描运动,其在Y方向分量用于补充线激光缺失的维度,同时传感器在Z方向根据蜂窝芯高度的变化做相应的运动,使待测面在线激光的量程范围内;每完成一次扫描后,向X轴正方向或负方向运动距离wx,wx小于线激光扫描宽度w0,传感器继续在YOZ平面内做扫描运动,以此类推进行往复多次扫描,直到完成整个蜂窝芯的测量。
进一步地,所述被测蜂窝芯表面重建时,线激光测头的误差角αZ、αy和αx,可以跟据各个角度对测量数据造成的误差,直接补偿到测量数据,其所对应的旋转矩阵分别为:
则其所对应的旋转矩阵Rl为:
进一步地,所述被测蜂窝芯表面重建时,线激光测头的误差角αZ、αy和αx,可以跟据各个角度的具体值,通过机械结构进行精确调节。
与现有技术相比,本发明具有以下有益效果:
1.测量精度高。本发明通过精确的计算,所得到的误差角精度高,进而精确地测量出蜂窝芯的表面。
2.操作方便。本发明所需要的测量步骤简单容易实现,使用方便。
3.适用性强。本发明能够适用于传感器需要不断安装拆卸的场合,可以实现在位测量。
4.成本低。不需要额外的标准件,只对待测的蜂窝芯进行测量,即可实现对线激光安装误差角的校准。
附图说明
为了更清楚的说明本发明的实施例或现有技术的技术方案,下面将对实施例或现有技术描述中所需要使用的附图做一简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是本发明的蜂窝芯表面轮廓测量方法流程图。
图2是本发明实施例中的线激光测量示意图。
图3是本发明实施例中的待测量蜂窝芯表面示意图。
图4是本发明实施例中线激光测头绕Z轴的安装误差角αZ的计算示意图。
图5是本发明实施例中线激光测头绕Y轴的安装误差角αy的计算示意图。
图6是本发明实施例中线激光测头绕X轴的安装误差角αx的计算示意图。
图7是本发明实施例一中线激光测量的蜂窝芯样件。
图8是本发明实施例一中线激光测量蜂窝芯样件的扫描路径。
图9是本发明实施例一中线激光测量蜂窝芯样件的测量数据结果。
图10是本发明实施例一中蜂窝芯样件的测量数据的其中一个拼接位置的俯视图,其中图中标注的直线为不同扫描位置的分界线。
图11是本发明实施例一中蜂窝芯样件的测量数据的其中一个拼接位置的侧视图,其中图中标注的直线为不同扫描位置的分界线。
具体实施方式
为使本发明的实施例的目的、技术方案和优点更加清楚,下面结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚完整的描述:
如图1所示:本发明采用的技术方案具体实现步骤为:
1.测量平台搭建
在三轴数控机床的工作台上搭建蜂窝芯的测量平台,其中O-XYZ坐标系为机床运动坐标系;将蜂窝芯固定于工作台,根据待测蜂窝芯表面对面的面形形状,设计对应支撑结构,使三个方向的蜂窝壁都沿着Z方向;
线激光传感器安装在三轴数控机床,使线激光的入射方向沿着Z方向;线激光沿Y轴做直线扫描运动,做固定步距扫描测量,本实施例中线激光沿Y轴做直线运动,沿X轴或其他可行的实施方案中,与沿Y轴同理,;利用Y轴电机的编码器脉冲信号或光栅尺反馈信号作为线激光测量的触发信号,编码器脉冲的间隔即为线激光测量的步距Δy;扫描时的状态如图2所示。
线激光测量时记录每条激光线对应的测头运动坐标;将机床的三个轴的编码器脉冲信号或者是光栅尺反馈信号,接入脉冲计数处理***,在Y轴电机的编码器脉冲信号或光栅尺反馈信号触发线激光测量时,同时触发该脉冲计数处理***记录下三个轴的脉冲数,根据脉冲数与实际运动的距离的关系,能够确定测头运动坐标[x0,y0,z0];
所述的蜂窝芯的形状不受限于某个具体的蜂窝芯形状;如相互连接的六边形,圆形,多边形,等规则与非规则的蜂窝芯均适用;示例如图3所示。
2.线激光测头的安装误差校准
线激光测头安装的倾斜误差可以依次分解为绕Z轴旋转αZ角度,绕Y轴旋转αy角度,绕X轴旋转αx角度;线激光测量后的数据坐标系与机床坐标系存在错位,另建立线激光测量数据坐标系为ol-xlylzl,其中xl轴和zl轴线激光测量数据本身的坐标,xl轴为激光线方向,zl轴为测量高度值,yl轴为通过机床运动为线激光补充的坐标轴,如图2所示。
1)蜂窝芯扫描
线激光三次扫描蜂窝芯:
线激光第一次扫描蜂窝芯,线激光沿Y轴运动,测量蜂窝芯的轮廓数据;
线激光第二次扫描蜂窝芯,线激光沿X轴正向平移距离w(w小于线激光的扫描宽度w0),以与第一次扫描时相同的Y向路径再次扫描蜂窝芯,测量蜂窝芯的轮廓数据;
第二次相对第一次扫描测头沿X轴正向平移距离w,沿Y向扫描路径相同其中,平移距离w满足:
w=β×w0 (1)
β取值为0~1;
两次测量的数据中沿X轴方向各有(1-β)w0宽度范围内的测量数据,对应相同的蜂窝芯轮廓;
线激光第三次扫描蜂窝芯,线激光沿Z轴正向抬高距离h(h小于线激光的量程),以与第二次扫描时相同的Y向路径再次扫描蜂窝芯,测量蜂窝芯的轮廓数据;h小于线激光的量程;
扫描前,调节线激光到蜂窝芯的距离,使三次扫描时蜂窝芯都在线激光测量量程范围内。
2)线激光测头绕Z轴的安装误差角αZ的计算与校准
如图4所示,在xlolyl平面内比较第一次和第二次测量的相同蜂窝芯轮廓数据,第二次测量相对于第一次测量,蜂窝芯轮廓沿yl轴方向的偏移量为lz,线激光测头绕Z轴的安装误差角满足:
tanαZ=lz/w (2)
αZ=arctan(lz/w) (3)
计算第一次和第二次测量蜂窝芯轮廓沿Y轴方向的偏移量为lz的具体过程为:
将第二次测量的数据沿yl轴正方向和负方向依次平移i行(i=-nmax,-nmax+1,…,-1,0,1,…nmax-1,nmax),其中nmax对应移动的数据行数范围,取值应大于预估的偏移量;计算平移后两次测量重合区域测量数据的重合数据量(或重合的面积),找到最大重合数据量(或重合的面积)对应的平移行数n1
偏移量为lz为:
lz=n1×Δy (4)
3)线激光测头绕Y轴的安装误差角αy的计算与校准
如图5所示,在xlolzl平面内比较第一次和第二次测量的相同蜂窝芯轮廓数据,第二次测量相对于第一次测量,蜂窝芯轮廓沿zl轴方向的偏移量ly,线激光测头绕Y轴的安装误差角满足:
tanαy=ly/w (5)
αy=arctan(ly/w) (6)
计算第一次和第二次测量蜂窝芯轮廓沿Z坐标方向的偏移量为ly的具体过程为:
将第二次测量的数据平移w个距离后,重复测量部分蜂窝芯的轮廓重合,计算重合的数据第二次测量相对于第一次测量的差值,ly取为所有差值的平均值。
4)线激光测头绕X轴的安装误差角αx的计算与校准
如图6所示,比较第二次和第三次测量的相同蜂窝芯轮廓数据,第三次测量相对于第二次测量,蜂窝芯轮廓沿zl轴方向的偏移量h',线激光测头绕Z轴的安装误差角αx满足:
cosαx=h/h' (7)
αx=arccos(h/h') (8)
第二次和第三次测量蜂窝芯轮廓沿Z轴方向的偏移量为h'的具体过程为:
第三次相对于第二次测量的在xlolyl平面内蜂窝芯的轮廓基本重合,计算重合的数据第三次测量相对于第二次测量的差值,h'取为所有差值的平均值。
3.扫描待测蜂窝芯
线激光测头校准后,由数控机床带动传感器按照一定路径扫描,确保整个蜂窝芯都被扫描到;
由数控机床带动传感器在YOZ平面内做扫描运动,其在Y方向分量用于补充线激光缺失的维度,同时传感器在Z方向根据蜂窝芯高度的变化做相应的运动,使待测面在线激光的量程范围内;每完成一次扫描后,向X轴正方向(或负方向)运动距离wx,wx小于线激光扫描宽度w0,传感器继续在YOZ平面内做扫描运动,以此类推进行往复多次扫描,直到完成整个蜂窝芯的测量。
4.被测蜂窝芯表面重建
根据线激光测头的安装误差角和所记录的机床的运行坐标,确定出坐标变换矩阵,将测量的在线激光坐标系下的测量数据经坐标变换统一到全局坐标系;坐标变换关系为:
P=Rl×Pl+T (9)
其中P(X,Y,Z)为测量后的数据在全局坐标系下对应的坐标,Pl(xl,0,zl)为线激光测量的数据,旋转变换矩阵Rl为:
Rl=Rxx)×Ryy)×RzZ) (10)
Rx,Ry,Rz分别为αx,αy,αZ三个误差角所对应的旋转变换矩阵;T=[x0,y0,z0]T为测量数据对应的传感器测头在机床中的坐标;
对全局坐标系下的测量数据进行蜂窝芯表面的重建,由于传感器沿X轴方向的进给小于激光线的宽度,相邻的两次扫描的数据有重叠的部分,对于重叠区域的数据,以重叠中心为界,分别选取两侧数据。
其特征还在于所述被测蜂窝芯表面重建时,线激光测头的误差角αZ、αy和αx,可以跟据各个角度对测量数据造成的误差,直接补偿到测量数据,其所对应的旋转矩阵分别为:
则其所对应的旋转矩阵Rl为:
另外,所述被测蜂窝芯表面重建时,线激光测头的误差角αZ、αy和αx,可以跟据各个角度的具体值,通过机械结构进行精确调节。
实施例1
1.测量平台搭建
本实施例所选用的激光位移传感器为基恩士公司的超高速轮廓测量仪(LJ-V7060)。该测量仪,采用蓝色半导体激光,可以实现稳定和超高速测量,采样间隔可达16μs,其Z轴方向重复测量精度可达0.4μm,测量轮廓数据间隔20μm,测量的激光线长度为15mm。测量时,将轮廓测量仪测头安装在三轴数控机床上,测量的标准件固定在机床工作台,由机床控制轮廓测量仪测头的运动,对标准件的表面进行扫描。
测量的蜂窝芯为具有六边形结构单元的蜂窝芯。
2.线激光测头的安装误差校准
1)蜂窝芯扫描
线激光三次扫描蜂窝芯:
线激光第一次扫描蜂窝芯,线激光沿Y轴运动,测量蜂窝芯的轮廓数据;
线激光第二次扫描蜂窝芯,线激光沿X轴正向平移距离10mm(线激光的扫描宽度15mm),以与第一次扫描时相同的Y向路径再次扫描蜂窝芯,测量蜂窝芯的轮廓数据;
线激光第三次扫描蜂窝芯,线激光沿Z轴正向抬高距离5mm(线激光的量程16mm),以与第二次扫描时相同的Y向路径再次扫描蜂窝芯,测量蜂窝芯的轮廓数据;
扫描前,调节线激光到蜂窝芯的距离,使三次扫描时蜂窝芯都在线激光测量量程范围内。
2)线激光测头绕Z轴的安装误差角αZ的计算与校准
如图4所示,在xlolyl平面内比较第一次和第二次测量的相同蜂窝芯轮廓数据,第二次测量相对于第一次测量,蜂窝芯轮廓沿yl轴方向的偏移量为lz,线激光测头绕Z轴的安装误差角满足:
tanαZ=lz/w(2)
αZ=arctan(lz/w) (3)
计算第一次和第二次测量蜂窝芯轮廓沿Y轴方向的偏移量为lz的具体过程为:
将第二次测量的数据沿yl轴正方向和负方向依次平移i行(i=-nmax,-nmax+1,…,-1,0,1,…nmax-1,nmax),其中nmax对应移动的数据行数范围,取值应大于预估的偏移量;计算平移后两次测量重合区域测量数据的重合数据量(或重合的面积),找到最大重合数据量(或重合的面积)对应的平移行数n1
偏移量为lz为:
lz=n1×Δy (4)
3)线激光测头绕Y轴的安装误差角αy的计算与校准
如图5所示,在xlolzl平面内比较第一次和第二次测量的相同蜂窝芯轮廓数据,第二次测量相对于第一次测量,蜂窝芯轮廓沿zl轴方向的偏移量ly,线激光测头绕Y轴的安装误差角满足:
tanαy=ly/w (5)
αy=arctan(ly/w) (6)
计算第一次和第二次测量蜂窝芯轮廓沿Z坐标方向的偏移量为ly的具体过程为:
将第二次测量的数据平移w个距离后,重复测量部分蜂窝芯的轮廓重合,计算重合的数据第二次测量相对于第一次测量的差值,ly取为所有差值的平均值。
4)线激光测头绕X轴的安装误差角αx的计算与校准
如图6所示,比较第二次和第三次测量的相同蜂窝芯轮廓数据,第三次测量相对于第二次测量,蜂窝芯轮廓沿zl轴方向的偏移量h',线激光测头绕Z轴的安装误差角αx满足:
cosαx=h/h' (7)
αx=arccos(h/h') (8)
第二次和第三次测量蜂窝芯轮廓沿Z轴方向的偏移量为h'的具体过程为:
第三次相对于第二次测量的在xlolyl平面内蜂窝芯的轮廓基本重合,计算重合的数据第三次测量相对于第二次测量的差值,h'取为所有差值的平均值。
3.扫描待测蜂窝芯
线激光测头校准后,由数控机床带动传感器按照一定路径扫描,确保整个蜂窝芯都被扫描到;
由数控机床带动传感器在YOZ平面内做扫描运动,其在Y方向分量用于补充线激光缺失的维度,同时传感器在Z方向根据蜂窝芯高度的变化做相应的运动,使待测面在线激光的量程范围内;每完成一次扫描后,向X轴正方向(或负方向)运动距离14mm,传感器继续在YOZ平面内做扫描运动,以此类推进行往复多次扫描,直到完成整个蜂窝芯的测量。
4.被测蜂窝芯表面重建
本实施例中,通过Matlab软件对被测蜂窝芯表面重建,根据线激光测头的安装误差角和所记录的机床的运行坐标,确定出坐标变换矩阵,将测量的在线激光坐标系下的测量数据经坐标变换统一到全局坐标系;坐标变换关系为:
P=Rl×Pl+T (9)
其中P(X,Y,Z)为测量后的数据在全局坐标系下对应的坐标,Pl(xl,0,zl)为线激光测量的数据,旋转变换矩阵Rl为:
Rl=Rxx)×Ryy)×RzZ) (10)
其中,Rx,Ry,Rz分别为αx,αy,αZ三个误差角所对应的旋转变换矩阵;T=[x0,y0,z0]T为测量数据对应的传感器测头在机床中的坐标;
对全局坐标系下的测量数据进行蜂窝芯表面的重建,由于传感器沿X轴方向的进给小于激光线的宽度,相邻的两次扫描的数据有重叠的部分,对于重叠区域的数据,以重叠中心为界,分别选取两侧数据。
其特征还在于所述被测蜂窝芯表面重建时,线激光测头的误差角αZ、αy和αx,可以跟据各个角度对测量数据造成的误差,直接补偿到测量数据,其所对应的旋转矩阵分别为:
则其所对应的旋转矩阵Rl为:
5.被测蜂窝芯结果观测
所扫描的蜂窝芯如图7所示,扫描时的路径如图8所示。利用本发明所提出的方法所构建的蜂窝芯表面轮廓整体图如图9所示,其中图10是所测量蜂窝芯样件的测量数据的其中一个拼接位置的俯视图,其中图中标注的直线为不同扫描位置的分界线。图11是所测量蜂窝芯样件的测量数据的其中一个拼接位置的侧视图,其中图中标注的直线为不同扫描位置的分界线。从分界线可以看出在拼接位置无明显错位情况,验证了本发明所提出的测量方法的可靠性。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,根据本发明的技术方案及其发明构思加以等同替换或改变,都应涵盖在本发明的保护范围之内。

Claims (10)

1.一种蜂窝芯表面轮廓测量方法,其特征在于包括如下步骤:
步骤1:测量平台搭建
在三轴数控机床的工作台上搭建蜂窝芯的测量平台,其中O-XYZ坐标系为机床运动坐标系;
将蜂窝芯固定于工作台,根据待测蜂窝芯表面对面的面形形状,设计对应支撑结构,使三个方向的蜂窝壁都沿着Z方向;
线激光传感器安装在三轴数控机床,使线激光的入射方向沿着Z方向;线激光沿Y轴或X轴做直线扫描运动,做固定步距扫描测量;线激光测量时记录每条激光线对应的测头运动坐标;
步骤2:确定线激光测头的安装误差角
将线激光测头安装的倾斜误差依次分解为绕Z轴旋转αZ角度,绕Y轴旋转αy角度,绕X轴旋转αx角度;线激光测量后的数据坐标系与机床坐标系存在错位,另建立线激光测量数据坐标系为ol-xlylzl,其中xl轴和zl轴为线激光测量数据本身的坐标,xl轴为激光线方向,zl轴为测量高度值,yl轴为沿机床运动方向为线激光补充的坐标轴;
为确定线激光测头绕Z轴的安装误差角αZ,绕Y轴的安装误差角αy和绕X轴的安装误差角αx,进行规定路径扫描蜂窝芯;
利用扫描后的数据计算出线激光测头绕Z轴的安装误差角αZ、绕Y轴的安装误差角αy和绕X轴的安装误差角αx
步骤3:扫描待测蜂窝芯
线激光测头校准后,由数控机床带动传感器按照一定路径扫描,确保整个蜂窝芯都被扫描到;
步骤4:被测蜂窝芯表面重建
根据线激光测头的安装误差角和所记录的机床的运行坐标,确定出坐标变换矩阵,将测量的在线激光坐标系下的测量数据经坐标变换统一到全局坐标系;坐标变换关系为:
P=Rl×Pl+T (9)
其中P(X,Y,Z)为测量后的数据在全局坐标系下对应的坐标,Pl(xl,0,zl)为线激光测量的数据,旋转变换矩阵Rl为:
Rl=Rxx)×Ryy)×RzZ) (10)
Rx,Ry,Rz分别为αx,αy,αZ三个误差角所对应的旋转变换矩阵;T=[x0,y0,z0]T为测量数据对应的传感器测头在机床中的坐标;
对全局坐标系下的测量数据进行蜂窝芯表面的重建。
2.根据权利要求1所述的蜂窝芯表面轮廓测量方法,其特征还在于所述线激光沿Y轴做直线运动,做固定步距扫描测量的实现方法为:
利用Y轴电机的编码器脉冲信号或光栅尺反馈信号作为线激光测量的触发信号,编码器脉冲的间隔即为线激光测量的步距Δy
将机床的三个轴的编码器脉冲信号或者是光栅尺反馈信号,接入脉冲计数处理***,在Y轴电机的编码器脉冲信号或光栅尺反馈信号触发线激光测量时,同时触发该脉冲计数处理***记录下三个轴的脉冲数,根据脉冲数与实际运动的距离的关系,能够确定测头运动坐标[x0,y0,z0]。
3.根据权利要求1所述的蜂窝芯表面轮廓测量方法,其特征还在于所述为确定线激光测头的3个安装误差角,进行扫描蜂窝芯的规定路径为:
线激光三次扫描蜂窝芯:
线激光第一次扫描蜂窝芯,线激光沿Y轴运动,测量蜂窝芯的轮廓数据;
线激光第二次扫描蜂窝芯,线激光沿X轴正向平移距离w,以与第一次扫描时相同的Y向路径再次扫描蜂窝芯,测量蜂窝芯的轮廓数据,其中,w小于线激光的扫描宽度w0
线激光第三次扫描蜂窝芯,线激光沿Z轴正向抬高距离h,以与第二次扫描时相同的Y向路径再次扫描蜂窝芯,测量蜂窝芯的轮廓数据,其中,h小于线激光的量程;
所述蜂窝芯的形状包括如相互连接的六边形,圆形,多边形的规则形状或非规则形状。
4.根据权利要求3所述的蜂窝芯表面轮廓测量方法,其特征还在于所述线激光三次扫描蜂窝芯时:
第二次相对第一次扫描测头沿X轴正向平移距离w,沿Y向扫描路径相同其中,平移距离w满足:
w=β×w0 (1)
β取值为0.5~0.75;
两次测量的数据中沿X轴方向各有(1-β)w0宽度范围内的测量数据,对应相同的蜂窝芯轮廓;
扫描前,调节线激光到蜂窝芯的距离,使三次扫描时蜂窝芯都在线激光测量量程范围内。
5.根据权利要求1~4任一项所述的蜂窝芯表面轮廓测量方法,其特征还在于所述计算线激光测头绕Z轴的安装误差角αZ的具体过程为:
在xlolyl平面内比较第一次和第二次测量的相同蜂窝芯轮廓数据,第二次测量相对于第一次测量,蜂窝芯轮廓沿yl轴方向的偏移量为lz,线激光测头绕Z轴的安装误差角满足:
tanαZ=lz/w (2)
αZ=arctan(lz/w) (3)
其中lz确定的具体过程为:
将第二次测量的数据沿yl轴正方向和负方向依次平移i行(i=-nmax,-nmax+1,…,-1,0,1,…nmax-1,nmax),其中nmax对应移动的数据行数范围,取值应大于预估的偏移量;计算平移后两次测量重合区域测量数据的重合数据量或重合的面积,找到最大重合数据量或重合的面积对应的平移行数n1
偏移量为lz为:
lz=n1×Δy (4)。
6.根据权利要求1~4任一项所述的蜂窝芯表面轮廓测量方法,其特征还在于所述计算线激光测头绕Y轴的安装误差角αy的具体过程为:
在xlolzl平面内比较第一次和第二次测量的相同蜂窝芯轮廓数据,第二次测量相对于第一次测量,蜂窝芯轮廓沿zl轴方向的偏移量ly,线激光测头绕Y轴的安装误差角满足:
tanαy=ly/w (5)
αy=arctan(ly/w) (6)
其中ly确定的具体过程为:
将第二次测量的数据平移w个距离后,重复测量部分蜂窝芯的轮廓重合,计算重合的数据第二次测量相对于第一次测量的差值,ly取为所有差值的平均值。
7.根据权利要求1~4任一项所述的蜂窝芯表面轮廓测量方法,其特征还在于所述计算线激光测头绕X轴的安装误差角αx的具体过程为:
比较第二次和第三次测量的相同蜂窝芯轮廓数据,第三次测量相对于第二次测量,蜂窝芯轮廓沿zl轴方向的偏移量h',线激光测头绕Z轴的安装误差角αx满足:
cosαx=h/h' (7)
αx=arccos(h/h') (8)
其中h'确定的具体过程为:
第三次相对于第二次测量的在xlolyl平面内蜂窝芯的轮廓基本重合,计算重合的数据第三次测量相对于第二次测量的差值,h'取为所有差值的平均值。
8.根据权利要求1所述的蜂窝芯表面轮廓测量方法,其特征还在于所述扫描待测蜂窝芯的具体过程为:
由数控机床带动传感器在YOZ平面内做扫描运动,其在Y方向分量用于补充线激光缺失的维度,同时传感器在Z方向根据蜂窝芯高度的变化做相应的运动,使待测面在线激光的量程范围内;每完成一次扫描后,向X轴正方向或负方向运动距离wx,wx小于线激光扫描宽度w0,传感器继续在YOZ平面内做扫描运动,以此类推进行往复多次扫描,直到完成整个蜂窝芯的测量。
9.根据权利要求1所述的蜂窝芯表面轮廓测量方法,其特征还在于所述被测蜂窝芯表面重建时,线激光测头的误差角αZ、αy和αx,可以跟据各个角度对测量数据造成的误差,直接补偿到测量数据,其所对应的旋转矩阵分别为:
则其所对应的旋转矩阵Rl为:
10.根据权利要求1所述的蜂窝芯表面轮廓测量方法,其特征还在于所述被测蜂窝芯表面重建时,线激光测头的误差角αZ、αy和αx,可以跟据各个角度的具体值,通过机械结构进行精确调节。
CN201910829610.5A 2019-09-03 2019-09-03 蜂窝芯表面轮廓测量方法 Active CN110487210B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910829610.5A CN110487210B (zh) 2019-09-03 2019-09-03 蜂窝芯表面轮廓测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910829610.5A CN110487210B (zh) 2019-09-03 2019-09-03 蜂窝芯表面轮廓测量方法

Publications (2)

Publication Number Publication Date
CN110487210A true CN110487210A (zh) 2019-11-22
CN110487210B CN110487210B (zh) 2020-12-11

Family

ID=68556259

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910829610.5A Active CN110487210B (zh) 2019-09-03 2019-09-03 蜂窝芯表面轮廓测量方法

Country Status (1)

Country Link
CN (1) CN110487210B (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111811425A (zh) * 2020-06-19 2020-10-23 成都飞机工业(集团)有限责任公司 一种航空复合材料蜂窝壁板蒙皮浅孔检测装置和检测方法
CN112763486A (zh) * 2020-11-30 2021-05-07 成都飞机工业(集团)有限责任公司 一种基于线激光扫描的复材壁板阵列孔检测方法
CN113865508A (zh) * 2021-09-28 2021-12-31 南京航空航天大学 一种蜂窝夹芯复合材料声衬的通孔率自动化检测装置与方法
CN114739290A (zh) * 2022-04-01 2022-07-12 上海拓璞数控科技股份有限公司 线激光扫描化铣胶刻线的路径规划方法及***
CN117647210A (zh) * 2024-01-29 2024-03-05 成都飞机工业(集团)有限责任公司 一种蜂窝芯构件复杂型面轮廓的检测方法

Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1078550A (zh) * 1992-04-07 1993-11-17 本田技研工业株式会社 非接触位置检测设备
CN102243067A (zh) * 2011-04-15 2011-11-16 中国科学院长春光学精密机械与物理研究所 太阳能聚光镜面形检测装置
EP2679362A2 (en) * 2012-06-27 2014-01-01 NGK Insulators, Ltd. Method for manufacturing ceramic honeycomb structure
CN104019772A (zh) * 2014-05-20 2014-09-03 泰安至诚公路工程监理咨询有限公司 一种桥梁结构砼表面蜂窝麻面面积测定仪及其检测方法
CN204613125U (zh) * 2015-05-18 2015-09-02 无锡惠科电工高新技术有限公司 一种线扫描式蜂窝陶瓷检测装置
CN105423881A (zh) * 2015-10-15 2016-03-23 哈尔滨工业大学 大型高速回转装备蜂窝面圆柱形测头测量及数据处理方法
EP3075718A1 (en) * 2015-03-31 2016-10-05 NGK Insulators, Ltd. Method for manufacturing ceramic formed body, and apparatus for manufacturing ceramic formed body
CN106017352A (zh) * 2016-07-22 2016-10-12 大连理工大学 一种蜂窝芯面形的测量方法
CN106017353A (zh) * 2016-07-22 2016-10-12 大连理工大学 一种蜂窝芯面形测量装置
US20180238683A1 (en) * 2015-10-28 2018-08-23 Ngk Insulators, Ltd. Method for inspecting end face and device for inspecting end face, of honeycomb structure
CN108648153A (zh) * 2018-04-25 2018-10-12 大连理工大学 蜂窝芯面形测量数据毛刺去除方法

Patent Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1078550A (zh) * 1992-04-07 1993-11-17 本田技研工业株式会社 非接触位置检测设备
CN102243067A (zh) * 2011-04-15 2011-11-16 中国科学院长春光学精密机械与物理研究所 太阳能聚光镜面形检测装置
EP2679362A2 (en) * 2012-06-27 2014-01-01 NGK Insulators, Ltd. Method for manufacturing ceramic honeycomb structure
CN104019772A (zh) * 2014-05-20 2014-09-03 泰安至诚公路工程监理咨询有限公司 一种桥梁结构砼表面蜂窝麻面面积测定仪及其检测方法
EP3075718A1 (en) * 2015-03-31 2016-10-05 NGK Insulators, Ltd. Method for manufacturing ceramic formed body, and apparatus for manufacturing ceramic formed body
CN204613125U (zh) * 2015-05-18 2015-09-02 无锡惠科电工高新技术有限公司 一种线扫描式蜂窝陶瓷检测装置
CN105423881A (zh) * 2015-10-15 2016-03-23 哈尔滨工业大学 大型高速回转装备蜂窝面圆柱形测头测量及数据处理方法
US20180238683A1 (en) * 2015-10-28 2018-08-23 Ngk Insulators, Ltd. Method for inspecting end face and device for inspecting end face, of honeycomb structure
CN106017352A (zh) * 2016-07-22 2016-10-12 大连理工大学 一种蜂窝芯面形的测量方法
CN106017353A (zh) * 2016-07-22 2016-10-12 大连理工大学 一种蜂窝芯面形测量装置
CN108648153A (zh) * 2018-04-25 2018-10-12 大连理工大学 蜂窝芯面形测量数据毛刺去除方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
YAN QIN等: "Detection of honeycomb cell walls from measurement data based on Harris corner detection algorithm", 《MEASUREMENT SCIENCE AND TECHNOLOGY》 *
刘蕾等: "一种针对蜂窝芯外形曲面特征的测量方法", 《中国测试》 *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111811425A (zh) * 2020-06-19 2020-10-23 成都飞机工业(集团)有限责任公司 一种航空复合材料蜂窝壁板蒙皮浅孔检测装置和检测方法
CN112763486A (zh) * 2020-11-30 2021-05-07 成都飞机工业(集团)有限责任公司 一种基于线激光扫描的复材壁板阵列孔检测方法
CN112763486B (zh) * 2020-11-30 2022-05-10 成都飞机工业(集团)有限责任公司 一种基于线激光扫描的复材壁板阵列孔检测方法
CN113865508A (zh) * 2021-09-28 2021-12-31 南京航空航天大学 一种蜂窝夹芯复合材料声衬的通孔率自动化检测装置与方法
CN114739290A (zh) * 2022-04-01 2022-07-12 上海拓璞数控科技股份有限公司 线激光扫描化铣胶刻线的路径规划方法及***
CN114739290B (zh) * 2022-04-01 2024-01-26 上海拓璞数控科技股份有限公司 线激光扫描化铣胶刻线的路径规划方法及***
CN117647210A (zh) * 2024-01-29 2024-03-05 成都飞机工业(集团)有限责任公司 一种蜂窝芯构件复杂型面轮廓的检测方法
CN117647210B (zh) * 2024-01-29 2024-05-14 成都飞机工业(集团)有限责任公司 一种蜂窝芯构件复杂型面轮廓的检测方法

Also Published As

Publication number Publication date
CN110487210B (zh) 2020-12-11

Similar Documents

Publication Publication Date Title
CN110487210A (zh) 蜂窝芯表面轮廓测量方法
JP2764103B2 (ja) アナログ測定プローブの使用方法および位置決め装置
US10209107B2 (en) Geometric error identification method of multi-axis machine tool and multi-axis machine tool
CN105404238B (zh) 一种在机激光测量中测头位姿的线性化标定方法
US4982504A (en) Method for determining positional errors and for compensating for such errors, and apparatus for carrying out the method
CN110530296A (zh) 一种线激光安装误差角确定方法
CN102607457B (zh) 基于惯性导航技术的大尺寸三维形貌测量装置及方法
US8786243B2 (en) Method and device for preparing error map and numerically controlled machine tool having error map preparation function
CN109794938A (zh) 一种适用于曲面结构的机器人制孔误差补偿装置及其方法
CN110500978A (zh) 点激光传感器的光束方向矢量和零点位置在线标定方法
CN107588742A (zh) 一种基于线结构光的圆柱齿轮齿廓偏差测量方法
CN102554710B (zh) 一种机床回转工作台转角误差的测量装置和方法
CN101334267A (zh) 数字影像测头矢量坐标变换标定与误差修正方法及其装置
CN103307984A (zh) 一种用于可调桨叶片的激光测量装置、***及方法
CN108827187A (zh) 一种用于对工件进行三维轮廓测量的测量***
CN106989670B (zh) 一种机器人协同的非接触式高精度大型工件跟踪测量方法
CN109520420A (zh) 一种转台回转中心的空间坐标确定方法
CN109163658A (zh) 一种可提供位置和角度基准的光学基准件的标定方法
CN103712572A (zh) 结构光源与相机结合的物体轮廓三维坐标测量装置
JPS63292005A (ja) 走り誤差補正をなした移動量検出装置
CN101788265B (zh) 发动机缸体结合面孔组快速测量全局统一标定方法
Iwasawa et al. Development of a measuring method for several types of programmed tool paths for NC machine tools using a laser displacement interferometer and a rotary encoder
TW201711796A (zh) 工具機旋轉軸定位精度檢測裝置與方法
CN109520417A (zh) 机床几何误差及旋转台转角定位误差检定装置和方法
CN111536877A (zh) 一种三坐标测量机上的线激光传感器姿态标定方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant