CN101957175B - 基于三点微平面式法向检测方法 - Google Patents

基于三点微平面式法向检测方法 Download PDF

Info

Publication number
CN101957175B
CN101957175B CN 201010281290 CN201010281290A CN101957175B CN 101957175 B CN101957175 B CN 101957175B CN 201010281290 CN201010281290 CN 201010281290 CN 201010281290 A CN201010281290 A CN 201010281290A CN 101957175 B CN101957175 B CN 101957175B
Authority
CN
China
Prior art keywords
coordinate
coordinate system
measuring machine
laser
displacement sensor
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.)
Expired - Fee Related
Application number
CN 201010281290
Other languages
English (en)
Other versions
CN101957175A (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.)
XIAN RUITE RAPID MANUFACTURE ENGINEERING Co Ltd
Xian Jiaotong University
Original Assignee
XIAN RUITE RAPID MANUFACTURE ENGINEERING Co Ltd
Xian Jiaotong University
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 XIAN RUITE RAPID MANUFACTURE ENGINEERING Co Ltd, Xian Jiaotong University filed Critical XIAN RUITE RAPID MANUFACTURE ENGINEERING Co Ltd
Priority to CN 201010281290 priority Critical patent/CN101957175B/zh
Publication of CN101957175A publication Critical patent/CN101957175A/zh
Application granted granted Critical
Publication of CN101957175B publication Critical patent/CN101957175B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Length Measuring Devices By Optical Means (AREA)

Abstract

本发明公开了一种基于三点微平面法向检测方法,基于微分几何思想——微平面近似代替微曲面法,并应用激光位移传感器技术及数据采集技术,通过一定的算法可以测得曲面上待测点的法向量。相关参数的标定方法是:首先建立一个测量空间,然后在该空间中拟定一球面并建立球面方程,将传感器安装好后在球面上采集若干数据点,这些数据点必然满足该球面方程,从而可得出若干个含参数的方程,再应用fminsearch最小优化方法求解方程组,便可解出各参数。本发明不仅能够进行非接触式测量,而且测量精度高、测量方便、快捷,可应用于法向误差校正、法向姿态调整、法向检测等测量领域。

Description

基于三点微平面式法向检测方法
技术领域
本发明涉及测量领域中的一种法向检测技术,具体是一种基于三点微平面法的、采用非接触式测量的曲面法向测量方法的设计及相关参数标定方法。
背景技术
在工业生产中,测量技术发挥着举足轻重的作用,随着技术的进步,特别是激光技术、传感器技术的发展,测量手段更加多样化和先进化,同时生产过程中要求的被测对象也越来越多,很多时候需要测量某一曲面上的法向量,如在飞机制造过程中的钻铆技术,就要求钻铆机的加工轴线和飞机蒙皮上待加工点处的法矢量重合;目前国内关于这方面的技术与文献比较少,因此提出曲面法向测量,并设计一种测量方法及相关参数标定方法,在工业生产中的测量领域具有重要的价值意义。
发明内容
本发明的目的在于提供一种曲面法向测量及相关参数的标定方法。本发明利用了激光传感器技术及数据采集技术并基于三点微平面法的设计思想,可以准确、快速、方便地进行曲面法向的测量。
本发明通过以下技术方案实现,实施步骤如下:
1)建立模型:将3个激光位移传感器环形均布于三坐标测量机的头部,传感器的空间姿态能使激光成锥形射在曲面上,根据精度要求调整锥度来控制3个激光点所形成的微平面大小;
2)建立空间坐标系:首先,用三坐标测量机建立一个世界坐标系1;然后选用三坐标测量机末端关节坐标系作为中间转换坐标系2;另外在每个激光位移传感器上建立一个以激光线射出方向为z轴负向、以激光位移传感器工作参考点为坐标原点的局部坐标系3;
3)采集数据:利用数据采集卡和LabVIEW软件采集激光位移传感器电压信号并转换为位移值Li
4)计算激光点的坐标:利用激光位移传感器测得的位移即可知道激光点在局部坐标系3下的坐标(0,0,wi),wi=-Li,根据欧拉坐标转换公式,分别求得3个激光点在坐标系1下的坐标,坐标转换公式如下:
P i = R 0 * R i * 0 0 w i 1 i=1,2,3
R i = R z 1 * R x 2 + T i
= 1 0 0 0 0 cos ( r z 1 ) sin ( r z 1 ) 0 0 - sin ( r z 1 ) cos ( r z 1 ) 0 0 0 0 1 * cos ( r x 2 ) sin ( r x 2 ) 0 0 - sin ( r x 2 ) cos ( r x 2 ) 0 0 0 0 1 0 0 0 0 1 + 0 0 0 a i 0 0 0 b i 0 0 0 c i 0 0 0 0
式中Pi表示激光点的坐标;R0是三坐标测量机末端位姿,从三坐标测量机软件中获取;Ri是局部坐标系3相对于中间转换坐标系2的齐次变换矩阵,其中rz1,rx2分别是坐标系3绕坐标系2的z轴,x轴旋转角度,(ai,bi,ci)是坐标系3对坐标系2的平移向量,在对rz1、rx2、(ai,bi,ci)进行标定;
5)计算法向量:用步骤4)求得的3个激光点坐标可以解出3个向量,再将其中任意2个向量叉乘即可得三点微平面法向量,公式如下:
P 2 P 1 → = ( P x 1 - P x 2 , P y 1 - P y 2 , P z 1 - P z 2 )
P 2 P 3 → = ( P x 3 - P x 2 , P y 3 - P y 2 , P z 3 - P z 2 )
N → = P 2 P 1 → × P 2 P 3 →
计算得到的法向量应用于法向误差校正、法向姿态调整、法向检测。
所述的激光位移传感器姿态参数标定方法:实施步骤如下:
1)按照(1)中步骤1)将激光位移传感器安装于三维坐标测量机头部并建立如下测量空间:按照(1)中步骤2)使用三坐标测量机建立一个世界坐标系1和建立一个中间转换坐标系2,然后在世界坐标系1中用标准球拟定一个球面方程(x-a)2+(y-b)2+(z-c)2=r2,其中,(a,b,c)表示球心坐标,r表示球半径,球心坐标用三坐标测量机标定,球半径已知;
2)以其中一个传感器操作三坐标测量机在球面上均匀采集30个测点,得到30个位移值;
3)将测得的位移值代入球面方程应用fminsearch最小优化方法求解各参数值;
4)重复步骤2)、步骤3)求解另外2个传感器的姿态参数。
本发明设计原理基于微分几何的思想-微平面近似代替微曲面,采用了激光传感器技术,不仅能够进行非接触式测量,而且测量精度高、测量方便、快捷;另外本发明中的标定算法简单、易实现,且在应用fminsearch最小优化求解时,对初值选取要求不高。
附图说明
图1是本发明的原理说明图;
图2是本发明的模型示意图;
(a)图是世界坐标系1;
(b)图是激光位移传感器布置示意图及中间转换坐标系2;
(c)是激光位移传感器局部坐标系3;
下面结合附图对本发明的内容作进一步详细说明。
具体实施方式
参照图1所示,是基于微分几何思想,在曲面上某点附近选取1个3点式微小平面,用此微小平面的法向量近似代替曲面上该点处的法向量,结果误差受微平面大小的影响。因此,微平面大小根据精度要求选取,1表示传感器。
参照图2所示,1表示传感器,2表示球面。本发明的实施步骤为:建立测量模型——姿态参数标定——曲面法向测量,具体实施方式如下:
(1)将3个激光传感器安装于三坐标测量机上,并使射出的3条激光构成锥形,根据精度要求调节锥度控制微平面的大小。
(2)使用三坐标测量机建立一个世界坐标系1,即附图2中的坐标系1,选用三坐标测量机末端关节处的坐标系为中间转换坐标系2,即附图2中坐标系2,另外分别沿激光线为Z轴,以传感器工作参考点为原点建立3个局部坐标系3。各坐标系的关系如下:坐标系1与坐标系2的变换关系即为三坐标测量机末端位姿,从三坐标测量机软件中获取;坐标系3与坐标系2的变换关系即激光位移传感器姿态参数待标定。
(3)在世界坐标系1中放置一标准球,使用三坐标测量机标定球心坐标,并根据标准球半径建立球面方程:(x-a)2+(y-b)2+(z-c)2=r2,用于参数标定。
(4)以单个传感器为研究对象,控制三坐标测量机在标准球面上采集30个测点,通过数据采集设备采集到各测点的位移值Li,可得到测点在局部坐标系3下的坐标(0,0,wi),wi=-Li
(5)进行坐标变换:将测点坐标(0,0,wi)变换到世界坐标系1中,变换公式为
P i = R 0 * R i * 0 0 w i 1 i=1,2,3
R i = R z 1 * R x 2 + T i
= 1 0 0 0 0 cos ( r z 1 ) sin ( r z 1 ) 0 0 - sin ( r z 1 ) cos ( r z 1 ) 0 0 0 0 1 * cos ( r x 2 ) sin ( r x 2 ) 0 0 - sin ( r x 2 ) cos ( r x 2 ) 0 0 0 0 1 0 0 0 0 1 + 0 0 0 a i 0 0 0 b i 0 0 0 c i 0 0 0 0
式中Pi表示激光点的坐标;R0是三坐标测量机末端位姿,从三坐标测量机软件中获取;Ri是局部坐标系3相对于中间转换坐标系2的齐次变换矩阵,其中rz1,rx2分别是坐标系3绕坐标系2的z轴,x轴旋转角度,(ai,bi,ci)是坐标系3对坐标系2的平移向量,各参数待标定。变换后坐标满足球面方程(xP-a)2+(yP-b)2+(zP-c)2=r2;得到30个球面方程。
(6)在matlab中应用fminsearch最小优化方法求解得到的30个球面方程,可解出坐标变换公式中各个参数。
(7)各参数求解出来后,则可以进行曲面法向测量。实现算法是:通过3个传感器采集数据,在曲面上获得3个点的空间坐标,可构成3个向量,任意2个向量叉乘可求得微平面法向量,即曲面法向量,算法用编程软件实现,得出的向量用于法向误差校正、法向姿态调整等测量领域。
在飞机自动化钻铆技术中,通常要求刀具沿蒙皮加工点处法向进给,以保证孔的加工质量。自动化过程中需要进行自动调姿,必须获取加工点处的法向,然后将法向反馈给控制***进行姿态调整。本发明刚好解决了如何实现法向的获取,而且利用了激光非接触式测量,测量精度高,方便等特点,能够符合实际应用需求。

Claims (2)

1.一种基于三点微平面法向检测方法,其特征在于:按如下步骤进行:
1)建立模型:将3个激光位移传感器环形均布于三坐标测量机的头部,传感器的空间姿态能使激光成锥形射在曲面上,根据精度要求调整锥度来控制3个激光点所形成的微平面大小;
2)建立空间坐标系:首先,用三坐标测量机建立一个世界坐标系1;然后选用三坐标测量机末端关节坐标系作为中间转换坐标系2;另外在每个激光位移传感器上建立一个以激光线射出方向为z轴负向、以激光位移传感器工作参考点为坐标原点的局部坐标系3;
3)采集数据:利用数据采集卡和LabVIEW软件采集激光位移传感器电压信号并转换为位移值Li
4)计算激光点的坐标:利用激光位移传感器测得的位移即可知道激光点在局部坐标系3下的坐标(0,0,wi),wi=-Li,根据欧拉坐标转换公式,分别求得3个激光点在坐标系1下的坐标,坐标转换公式如下:
P i = R 0 * R i * 0 0 w i 1 , i=1,2,3
R i = R z 1 * R x 2 + T i
= 1 0 0 0 0 cos ( r z 1 ) sin ( r z 1 ) 0 0 - sin ( r z 1 ) cos ( r z 1 ) 0 0 0 0 1 * cos ( r x 2 ) sin ( r x 2 ) 0 0 - sin ( r x 2 ) cos ( r x 2 ) 0 0 0 0 1 0 0 0 0 1 + 0 0 0 a i 0 0 0 b i 0 0 0 c i 0 0 0 0
式中Pi表示激光点的坐标;R0是三坐标测量机末端位姿,从三坐标测量机软件中获取;Ri是局部坐标系3相对于中间转换坐标系2的齐次变换矩阵,其中rz1,rx2分别是坐标系3绕坐标系2的z轴,x轴旋转角度,(ai,bi,ci)是坐标系3对坐标系2的平移向量,在对rz1、rx2、(ai,bi,ci)进行标定;
5)计算法向量:用步骤4)求得的3个激光点坐标可以解出3个向量,再将其中任意2个向量叉乘即可得三点微平面法向量,公式如下:
P 2 P 1 → = ( P x 1 - P x 2 , P y 1 - P y 2 , P z 1 - P z 2 )
P 2 P 3 → = ( P x 3 - P x 2 , P y 3 - P y 2 , P z 3 - P z 2 )
N → = P 2 P 1 → × P 2 P 3 →
计算得到的法向量应用于法向误差校正、法向姿态调整、法向检测。
2.根据权利要求1所述的方法,其特征在于,所述的激光位移传感器姿态参数标定方法:实施步骤如下:
(1)按照权利要求1中步骤1)将激光位移传感器安装于三维坐标测量机头部并建立如下测量空间:按照权利要求1中步骤2)使用三坐标测量机建立一个世界坐标系1和建立一个中间转换坐标系2,然后在世界坐标系1中用标准球拟定一个球面方程(x-a)2+(y-b)2+(z-c)2=r2,其中,(a,b,c)表示球心坐标,r表示球半径,球心坐标用三坐标测量机标定,球半径已知;
(2)以其中一个传感器操作三坐标测量机在球面上均匀采集30个测点,得到30个位移值;
(3)将测得的位移值代入球面方程应用fminsearch最小优化方法求解各参数值;
(4)重复步骤(2)、步骤(3)求解另外2个传感器的姿态参数。
CN 201010281290 2010-09-14 2010-09-14 基于三点微平面式法向检测方法 Expired - Fee Related CN101957175B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 201010281290 CN101957175B (zh) 2010-09-14 2010-09-14 基于三点微平面式法向检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 201010281290 CN101957175B (zh) 2010-09-14 2010-09-14 基于三点微平面式法向检测方法

Publications (2)

Publication Number Publication Date
CN101957175A CN101957175A (zh) 2011-01-26
CN101957175B true CN101957175B (zh) 2012-03-28

Family

ID=43484648

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 201010281290 Expired - Fee Related CN101957175B (zh) 2010-09-14 2010-09-14 基于三点微平面式法向检测方法

Country Status (1)

Country Link
CN (1) CN101957175B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108332909A (zh) * 2018-02-07 2018-07-27 哈尔滨工业大学 一种用于质心测量的多点称重支撑结构模型的修正方法

Families Citing this family (26)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102248450B (zh) * 2011-04-20 2013-03-20 上海交通大学 用于大曲率半径曲面法向矢量快速检测方法
CN102284956B (zh) * 2011-05-04 2013-11-20 南京航空航天大学 一种自动钻铆机器人的法向找正方法
CN102322794B (zh) * 2011-05-12 2013-06-05 浙江大学 基于led环形光的实时半径补偿方法及测头装置
CN102636110B (zh) * 2012-03-26 2015-03-25 南京航空航天大学 飞机部件自动钻铆***的基准检测装置及其检测方法
CN102768006B (zh) * 2012-08-02 2015-05-13 西北工业大学 一种用于大曲率曲面零件法向自动测量和调整的方法
CN103419212A (zh) * 2013-05-20 2013-12-04 北京航天金宇电子技术有限公司 一种任意曲面法向找正方法及装置
CN103307968B (zh) * 2013-06-03 2016-11-02 西北工业大学 机器人承载平台姿态的检测方法
CN103447877B (zh) * 2013-09-16 2015-08-19 南京航空航天大学 自主移动机构法矢检测与调姿运动方法
CN103528519B (zh) * 2013-10-21 2016-08-17 中国航空工业集团公司北京航空精密机械研究所 一种点投射光学测头空间位置矢量的测量标定方法
CN104615230B (zh) * 2013-11-01 2019-04-26 腾讯科技(深圳)有限公司 一种操作识别方法和终端
CN103968767B (zh) * 2014-05-14 2016-06-29 哈尔滨工程大学 一种球体识别的激光测量方法
CN105345600B (zh) * 2015-10-20 2020-02-21 上海拓璞数控科技股份有限公司 一种用于航空薄壁件自由曲面钻孔的法向测量与调整方法
CN105783722B (zh) * 2016-04-12 2018-04-06 南京航空航天大学 一种基于四面体构型的距离传感器法矢测量方法
CN106123817B (zh) * 2016-06-14 2018-11-27 昆明理工大学 一种基于任意空间距离的标定***及方法
CN107063086B (zh) * 2017-03-13 2019-03-12 浙江优迈德智能装备有限公司 一种基于双红外传感器的局部曲面法向量确定方法
CN107063089A (zh) * 2017-04-10 2017-08-18 北京航空航天大学 一种曲面法线测量精度验证方法
CN108375352A (zh) * 2018-01-31 2018-08-07 歌尔股份有限公司 基于平行度的两平面位置的调整方法、装置及***
CN109211104B (zh) * 2018-03-22 2020-09-08 中国航空制造技术研究院 一种加工孔位法向在线修正方法
CN109032069B (zh) * 2018-07-19 2020-11-03 西南交通大学 一种采用电涡流位移传感器的非接触式R-test测量仪球心坐标计算方法
CN110842931B (zh) * 2019-07-30 2022-03-22 南京埃斯顿机器人工程有限公司 一种应用于机器人打孔的工具姿态调整方法
CN110470223B (zh) * 2019-09-16 2020-08-25 大连理工大学 一种基于1d位移传感器的三维空间坐标测量方法
CN113074628B (zh) * 2020-01-03 2024-02-27 上海飞机制造有限公司 一种修边型面检验的装置及方法
CN112525131B (zh) * 2020-10-23 2021-09-14 清华大学 非接触式距离传感器阵列测量工件曲率的方法及***
CN112432622A (zh) * 2020-11-05 2021-03-02 中国航空制造技术研究院 基于单测距传感器的曲面法向的测量装置及方法
CN113008132B (zh) * 2021-02-23 2022-05-17 中山大学 一种基于cqp的激光干涉仪与光轴精密定位调装装置及方法
CN114111672A (zh) * 2021-11-26 2022-03-01 南京航空航天大学 一种多位移传感器法向测量的传感器安装位置参数快速标定方法

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3466661B2 (ja) * 1993-06-29 2003-11-17 キヤノン株式会社 画像処理装置及びその方法
KR100457080B1 (ko) * 2002-07-12 2004-11-12 삼성중공업 주식회사 영상에 의한 암반사면 절리의 기하학적 특성조사방법
DE10313038B4 (de) * 2003-03-24 2005-02-17 Klingelnberg Gmbh Vorrichtung zur Erfassung der Lage eines Tastelements in einem Mehrkoordinatenmessgerät
CN100453969C (zh) * 2004-12-31 2009-01-21 上海萨克斯动力总成部件***有限公司 用于测量复杂曲面的三坐标机测量方法及其装置
WO2008014461A2 (en) * 2006-07-28 2008-01-31 Optimet, Optical Metrology Ltd. Double-sided measurement of dental objects using an optical scanner

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108332909A (zh) * 2018-02-07 2018-07-27 哈尔滨工业大学 一种用于质心测量的多点称重支撑结构模型的修正方法

Also Published As

Publication number Publication date
CN101957175A (zh) 2011-01-26

Similar Documents

Publication Publication Date Title
CN101957175B (zh) 基于三点微平面式法向检测方法
CN110370286B (zh) 基于工业机器人和单目相机的定轴运动刚体空间位置识别方法
CN108008177B (zh) 多轴磁阻电流测量方法、装置、设备和***
CN101907448B (zh) 一种基于双目三维视觉的深度测量方法
CN102607459B (zh) Lidar测量数据的拼接方法和装置
CN103886593A (zh) 一种基于三维点云曲面圆孔检测方法
CN102107374A (zh) 轴盘类零件直径尺寸在线检测方法
CN105066915A (zh) 模具曲面加工误差和表面粗糙度在机检测装置及检测方法
CN109124769B (zh) 手术机器人的坐标系标定、控制的方法及***
CN104482911A (zh) 基于误差球的球度误差评定方法
CN103759672A (zh) 雪糕棒平面轮廓尺寸的视觉测量方法
CN104359415B (zh) 一种水火弯板角变形量测量方法及***
CN103673976A (zh) 复合式精度测量坐标系的转换与统一方法及***
CN108317999B (zh) 一种基于惯性传感器的通信铁塔倾角测量方法
CN106871859A (zh) 一种刚体空间定轴转动角度检测方法
CN208780164U (zh) 全深度阵列式测斜仪
CN109115147A (zh) 全深度阵列式测斜仪及深部位移测量方法
CN106568365B (zh) 一种球面孔系复合位置度误差的检测与评定方法
CN107990871A (zh) 一种基于多传感信息融合的相对位姿测量装置与方法
Hou et al. Automatic recognition system of pointer meters based on lightweight CNN and WSNs with on-sensor image processing
CN103075977A (zh) 双目立体视觉***中的点云数据的自动拼合算法
CN110940271A (zh) 基于空间三维测控网的船舶等大型工业制造检测、监测和智能搭载安装方法
CN103040471B (zh) 一种奶牛体型线性鉴定指标获取***及方法
CN109084734A (zh) 基于单目显微视觉的微球姿态测量装置及测量方法
CN204881558U (zh) 模具曲面加工误差和表面粗糙度在机检测装置

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20120328

Termination date: 20150914

EXPY Termination of patent right or utility model