CN105737762B - 一种航空发动机叶片型面测量方法 - Google Patents
一种航空发动机叶片型面测量方法 Download PDFInfo
- Publication number
- CN105737762B CN105737762B CN201610300993.3A CN201610300993A CN105737762B CN 105737762 B CN105737762 B CN 105737762B CN 201610300993 A CN201610300993 A CN 201610300993A CN 105737762 B CN105737762 B CN 105737762B
- Authority
- CN
- China
- Prior art keywords
- msub
- mrow
- plane
- formula
- blade
- 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.)
- Active
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/24—Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures
- G01B11/25—Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures by projecting a pattern, e.g. one or more lines, moiré fringes on the object
- G01B11/2518—Projection by scanning of the object
- G01B11/2522—Projection by scanning of the object the position of the object changing and being recorded
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/24—Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures
- G01B11/2433—Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures for measuring outlines by shadow casting
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明公开了一种航空发动机叶片型面测量方法,所述测量方法基于一个测量系统,所述测量系统包括计算机、双目工业相机、投影仪、转台和定位平面;所述测量方法是:先对整个系统进行标定,再以定位平面为基准对转台进行定位,然后将叶片固定在转台上,由投影仪投影蓝色结构光到被测叶片上,双目工业相机分别抓拍变形条纹,然后进行相位解算及匹配,得到单次测量的叶片点云数据,计算机上的程序自动控制转动转台,重复投影抓拍解算过程,得到叶片的360°多视角测量数据,利用定位的四元数及旋转中心对测量数据进行计算,自动完成叶片型面的拼接,获得叶片型面的完整融合数据。本发明具有测量精度高和测量效率高的双重优势。
Description
技术领域
本发明涉及叶片型面的测量技术,具体是一种航空发动机叶片型面的测量方法。
背景技术
叶片是航空发动机的核心部件,其型面轮廓的测量技术是逆向工程、质量检测、叶片设计以及维修指导等作业的重要步骤,因此,关于航空发动机叶片型面的测量技术一直以来都是研究的热点。
目前,对于航空发动机叶片的测量技术通常分为两类。其一,基于单点的高精度测量方法,如三坐标法,一些研究采用改进辅助测具与测量路径设计的方法来进行测量,此类测量方法的测量精度高,但其存在的问题是:测量效率低、测量成本高、且无法获取全部的面型信息。其二,采用面结构光进行三维测量的方法,典型的是相位测量轮廓方法,此类测量方法单次测量可以获取视场内全部的面型信息,相较三坐标法而言,其测量精度较低,然而,它具有测量效率高、可获取整体面型信息等优点,是目前行业内的研究热点。
相位测量方法要想获取全部结构的面型信息,必须要进行多次测量数据的融合。目前,行业内通常利用商业软件进行手工融合,基于人为因素的影响,存在极大的融合误差,进而影响融合精度。在数据的自动融合方面,可以采取标记点的方式进行一定量的高精度的数据。然而,叶片是有特殊要求的,标记点的融合虽然可以获得一定量的高精度的数据,由于需要在被测物体表面贴标记点,这势必对被测量表面的数据有一定的影响,因而,此技术不适合叶片的测量。
发明内容
本发明的发明目的在于:针对上述航空发动机叶片的特殊性以及现有测量技术的不足,提供一种具有精度和效率双重优势的航空发动机叶片型面测量方法。
本发明所采用的技术方案是,一种航空发动机叶片型面测量方法,所述测量方法基于一个测量系统,该测量系统包括计算机、双目工业相机、投影仪、转台和定位平面,其中,所述计算机用于控制转台、投影仪、双目工业相机和处理测量数据,所述投影仪用于投影结构光到被测物体上,所述双目工业相机用于抓拍被测物体的变形条纹图,所述定位平面用于转台的定位,作为参照对比的基准;所述测量方法是:先对整个系统进行标定,再以定位平面为基准对转台进行定位,然后将叶片固定在转台上,由投影仪投影蓝色结构光到被测叶片上,双目工业相机分别抓拍变形条纹,然后进行相位解算及匹配,得到单次测量的叶片点云数据,计算机上的处理程序自动控制转动转台,重复投影抓拍解算过程,得到叶片的360°多视角测量数据,利用定位的四元数及旋转中心对测量数据进行计算,自动完成叶片型面的拼接,获得叶片型面的完整融合数据。
所述利用定位的四元数及旋转中心对测量数据进行计算的具体步骤是:
1).平面拟合
设转台转动两次,每次转动角度为α,共有三个测量平面,定义转台定位后的初始平面为P1,第一次转动后的平面为P2,第二次转动后的平面为P3,经过三次抓拍采集,得到三个平面的点云数据集Pt(xk,yk,zk),其中t=1,2,3;k=1,2,3……n;
对平面数据进行拟合,设平面模型为
atxk+btyk+ctzk+dt=0 (1)
式中,at,bt,ct,dt为待拟合平面Pt(t=1,2,3)的平面参数;Pt(at,bt,ct)为平面法向量;
按下式用最小二乘法拟合平面方程,使采集的点云数据坐标值Pt(xk,yk,zk)到平面模型的距离和最小;
L=∑D2 t=∑(at*xk+bt*yk+ct*zk+dt)2/(at 2+bt 2+ct 2) (2)
式中,D2 t是点云数据坐标值Pt(xk,yk,zk)到平面模型的距离;L是点云数据坐标值Pt(xk,yk,zk)到平面模型的距离和;
2).计算虚拟旋转角度
按下式分别计算出平面P2和P1、P3的平面夹角θ21和θ23,
式中,θij为平面Pi和Pj的夹角,其中i=2,j=1,3;
再按下式计算虚拟旋转角w
w=(θ21+θ23)/2 (4)
3).计算虚拟旋转轴
设O21和O23分别为平面P1和P2、P2和P3的交线,平面P21为过O21的角θ21的平分面,平面P23为过O23的θ23的角平分面,平面P21和P23的交线即为转轴;
通过计算平面向量叉积,分别得到交线O21和O23的法向量n21和n23,算式如下
nij=(bi*cj-ci*bj,ci*aj-ai*cj,ai*bj-bi*aj) (5)
分别计算平面P2和P1、P3的法线夹角角平分线的向量m21和m23,算式如下
mij=((ai+aj)/2,(bi+bj)/2,(ci+cj)/2) (6)
则m21和n21的叉积为平面P21的法向量r21,m23和n23的叉积则为平面P23的法向量r23,算式如下
rij=(mij[2]*nij[3]-mij[3]*nij[2],mij[3]*nij[1]-mij[1]*nij[3],mij[1]*nij[2]-mij[2]*nij[1]) (7)
设角平分面Pij的方程为aijx+bijy+cijz+dij=0,则aij=rij[1],bij=rij[2],cij=rij[3],系统dij的计算公式为
dij=di*T1+dj*T2 (8)
式中,T1和T2的定义分别如公式(9)和公式(10)
T1=(aj*bij-aij*bj)/(aj*bi-ai*bj) (9)
T2=(aij*bi-ai*bij)/(aj*bi-ai*bj) (10)
计算平面P21和P23的法线叉积,就可以得到转轴的方向矢量v,公式为
v=(b21*c23-c21*b23,c21*a23-a21*c23,a21*b23-b21*a23) (11)
4).计算旋转中心点
得到旋转轴的方向矢量后,再计算旋转轴中心点(Xo,Yo,Zo),就可以得旋转轴,计算公式如下
5).四元数计算
利用上面计算得到的旋转轴矢量和角度,可以进一步转化为四元数quat,计算公式如下
quat=(cosw,v[1]*sinw,v[2]*sinw,v[3]*sinw) (13)
6).测量数据定位融合
计算得到四元数和旋转中心后,转台每次旋转角度β,则每次测量的数据G(xs,ys,zs)(s=1,2,3…n)相对中心点按照公式(14)做平移运算后生成数据G、(xs、,ys、,zs、),之后按照公式(15)运算,最终计算得到的点云数据G、、(xs、、,ys、、,zs、、),直接完成数据的拼接;
G、(xs、,ys、,zs、)=G(xs-Xo,ys-Yo,zs-Zo),(s=1,2,3...n) (14)
其中quat、,tn0、tn1、tn2以及tn3的定义如公式(16)、公式(17)所示
quat、=quat*α/β (16)
在测量数据定位融合之后,还利用ICP算法,对转台机械误差进行补偿。
本发明的有益效果是:
1.本发明通过对虚拟转轴及旋转中心的高精度定位,进而实现了航空发动机叶片型面的多视角点云数据的融合,从而很好的解决了结构光测量叶片轮廓的难题;从实际实验和测量结果分析来看,本发明具有自动测量、测量精度高、测量成本低、无需辅助设备及二次操作等优点;
2.本发明结合ICP(iterative closest point)算法,能够有效补偿测量系统转台的机械误差带来的拟合误差,进而可靠地实现了航空发动机叶片复杂型面的高精度360°扫描且自动融合,具有精度和效率的双重优势。
附图说明
下面结合附图对本发明作进一步的说明。
图1是本发明所采用测量系统的结构示意图。
图2是本发明的测量流程图。
图3是本发明的转台进行平面拟合的示意图。
图4是本发明的定位原理示意图。
具体实施方式
本发明为航空发动机叶片型面的测量方法,首先,该测量方法基于一个测量系统(已公开的现有技术),参见图1所示,该测量系统包括计算机A、双目工业相机(B1、B2)、投影仪C、转台D和定位平面E。其中,计算机A用于控制转台D、投影仪C和双目工业相机(B1、B2)动作,并且,该计算机A上设有对应的数据处理程序,通过处理程序对扫描采集到的测量数据进行处理。投影仪C用于投影结构光到转台D上的被测物体-叶片F(即航空发动机叶片)上。双目工业相机(B1、B2)用于抓拍投影仪C投影到被测物体-叶片F上以后所产生的变形条纹图。定位平面E用于对转台D进行定位,作为参照对比的基准。
参见图2所示,本发明的测量方法是:先对整个测量系统进行标定,再以定位平面E为基准对转台D进行定位,然后将叶片F固定在转台D上,由投影仪C投影蓝色结构光到被测叶片F上,双目工业相机(B1、B2)分别抓拍变形条纹,然后进行相位解算及匹配,得到单次测量的叶片点云数据,计算机A上的处理程序自动控制转动转台D,重复投影抓拍解算过程,得到叶片F的360°多视角测量数据,利用定位的四元数及旋转中心对测量数据进行计算,自动完成叶片型面的拼接,获得叶片型面的完整融合数据;最后,利用ICP算法,对转台D的机械误差进行补偿,获取高精度的被测物体-叶片F的完整融合数据。
上述利用定位的四元数及旋转中心对测量数据进行计算的具体步骤是:
1).平面拟合
在测量系统中,平面的精度及拟合精度直接影响到转轴及角度的定位精度,因而对于测量平面,选择高精度的陶瓷平晶作为基准平面,该基准平面的精度在微米级;平面拟合的具体措施是(参见图3所示),设转台转动两次,每次转动角度为α,共有三个测量平面,定义转台定位后的初始平面为P1,第一次转动后的平面为P2,第二次转动后的平面为P3,经过三次抓拍采集,得到三个平面的点云数据集Pt(xk,yk,zk),其中t=1,2,3;k=1,2,3……n;
对平面数据进行拟合,设平面模型为
atxk+btyk+ctzk+dt=0 (1)
式中,at,bt,ct,dt为待拟合平面Pt(t=1,2,3)的平面参数;Pt(at,bt,ct)为平面法向量;
对于测量后的平面点云数据的拟合,采用最小二乘法实现,具体是,按下式用最小二乘法拟合平面方程,使采集的点云数据坐标值Pt(xk,yk,zk)到平面模型的距离和最小;
L=∑D2 t=∑(at*xk+bt*yk+ct*zk+dt)2/(at 2+bt 2+ct 2) (2)
式中,D2 t是点云数据坐标值Pt(xk,yk,zk)到平面模型的距离;L是点云数据坐标值Pt(xk,yk,zk)到平面模型的距离和;
2).计算虚拟旋转角度
参见图4所示,旋转定位平面P2和P1、P3之间分别具有夹角,按下式分别计算出平面P2和P1、P3的平面夹角θ21和θ23,
式中,θij为平面Pi和Pj的夹角,其中i=2,j=1,3;;
再按下式计算虚拟旋转角w
w=(θ21+θ23)/2 (4)
3).计算虚拟旋转轴
参见图4所示,设O21和O23分别为平面P1和P2、P2和P3的交线,平面P21为过O21的角θ21的平分面,平面P23为过O23的θ23的角平分面,平面P21和P23的交线即为转轴O;
通过计算平面向量叉积,分别得到交线O21和O23的法向量n21和n23,算式如下
nij=(bi*cj-ci*bj,ci*aj-ai*cj,ai*bj-bi*aj) (5)
分别计算平面P2和P1、P3的法线夹角角平分线的向量m21和m23,算式如下
mij=((ai+aj)/2,(bi+bj)/2,(ci+cj)/2) (6)
则m21和n21的叉积为平面P21的法向量r21,m23和n23的叉积则为平面P23的法向量r23,算式如下
rij=(mij[2]*nij[3]-mij[3]*nij[2],mij[3]*nij[1]-mij[1]*nij[3],mij[1]*nij[2]-mij[2]*nij[1]) (7)
设角平分面Pij的方程为aijx+bijy+cijz+dij=0,则aij=rij[1],bij=rij[2],cij=rij[3],系统dij的计算公式为
dij=di*T1+dj*T2 (8)
式中,T1和T2的定义分别如公式(9)和公式(10)
T1=(aj*bij-aij*bj)/(aj*bi-ai*bj) (9)
T2=(aij*bi-ai*bij)/(aj*bi-ai*bj) (10)
计算平面P21和P23的法线叉积,就可以得到转轴的方向矢量v,公式为
v=(b21*c23-c21*b23,c21*a23-a21*c23,a21*b23-b21*a23) (11)
4).计算旋转中心点
得到旋转轴的方向矢量后,再计算旋转轴中心点(Xo,Yo,Zo),就可以得旋转轴,计算公式如下
5).四元数计算
利用上面计算得到的旋转轴矢量和角度,可以进一步转化为四元数quat,计算公式如下
quat=(cosw,v[1]*sinw,v[2]*sinw,v[3]*sinw) (13)
6).测量数据定位融合
计算得到四元数和旋转中心后,转台每次旋转角度β,则每次测量的数据G xs,ys,zs)(s=1,2,3…n)相对中心点按照公式(14)做平移运算后生成数据G、(xs、,ys、,zs、),之后按照公式(15)运算,最终计算得到的点云数据G、、(xi、、,yi、、,zi、、),直接完成数据的拼接;
G、(xs、,ys、,zs、)=G(xs-Xo,ys-Yo,zs-Zo),(s=1,2,3...n) (14)
其中quat、,tn0、tn1、tn2以及tn3的定义如公式(16)、公式(17)所示
quat、=quat*α/β (16)
为了验证上述算法的有效性,对定位平面P1、P2和P3的数据代入上述公式模拟旋转过程,以P3为基准,模拟P1旋转两次,P2旋转一次,利用求取的四元数和旋转中心进行处理后,平面P1、P2和P3完全重合,通过标准偏差分析,P3和P2的标准偏差为0.021mm,P3和P1的标准偏差为0.026mm,由于扫描仪的精度为0.02mm,因此可以看出上述算法具有优异的精度。
定位融合后,由于测量物体-叶片要旋转360°,基于转台旋转角度的机械误差,每次旋转的角度有误差,这对于大尺寸的测量,该误差对实际的测量结果融合影响很小,但对于高精度的小尺寸叶片测量而言,随着旋转次数的增加,误差会出现累加,进而该误差会足够影响到实际的融合结果,因此,为了消除该影响,在后处理中采用了ICP算法(已公开的现有技术)进行误差消除,以此对测量系统的转台的机械误差进行有效、可靠地补偿,获取高精度的叶片的完整融合数据。
上述具体技术方案仅用以说明本发明,而非对其限制。尽管参照上述具体技术方案对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对上述具体技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明技术方案的精神和范围。
Claims (2)
1.一种航空发动机叶片型面测量方法,所述测量方法基于一个测量系统,该测量系统包括计算机、双目工业相机、投影仪、转台和定位平面,其中,所述计算机用于控制转台、投影仪、双目工业相机和处理测量数据,所述投影仪用于投影结构光到被测物体上,所述双目工业相机用于抓拍被测物体的变形条纹图,所述定位平面用于转台的定位,作为参照对比的基准;其特征在于,所述测量方法是:先对整个系统进行标定,再以定位平面为基准对转台进行定位,然后将叶片固定在转台上,由投影仪投影蓝色结构光到被测叶片上,双目工业相机分别抓拍变形条纹,然后进行相位解算及匹配,得到单次测量的叶片点云数据,计算机上的处理程序自动控制转动转台,重复投影抓拍解算过程,得到叶片的360°多视角测量数据,利用定位的四元数及旋转中心对测量数据进行计算,自动完成叶片型面的拼接,获得叶片型面的完整融合数据;
所述利用定位的四元数及旋转中心对测量数据进行计算的具体步骤是:
1).平面拟合
设转台转动两次,每次转动角度为α,共有三个测量平面,定义转台定位后的初始平面为P1,第一次转动后的平面为P2,第二次转动后的平面为P3,经过三次抓拍采集,得到三个平面的点云数据集Pt(xk,yk,zk),其中t=1,2,3;k=1,2,3……n;
对平面数据进行拟合,设平面模型为
atxk+btyk+ctzk+dt=0 (1)
式中,at,bt,ct,dt为待拟合平面Pt(t=1,2,3)的平面参数;Pt(at,bt,ct)为平面法向量;
按下式用最小二乘法拟合平面方程,使采集的点云数据坐标值Pt(xk,yk,zk)到平面模型的距离和最小;
L=∑D2 t=∑(at*xk+bt*yk+ct*zk+dt)2/(at 2+bt 2+ct 2) (2)
式中,D2 t是点云数据坐标值Pt(xk,yk,zk)到平面模型的距离;L是点云数据坐标值Pt(xk,yk,zk)到平面模型的距离和;
2).计算虚拟旋转角度
按下式分别计算出平面P2和P1、P3的平面夹角θ21和θ23,
<mrow>
<msub>
<mi>&theta;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>=</mo>
<mi>arccos</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>a</mi>
<mi>i</mi>
</msub>
<msub>
<mi>a</mi>
<mi>j</mi>
</msub>
<mo>+</mo>
<msub>
<mi>b</mi>
<mi>i</mi>
</msub>
<msub>
<mi>b</mi>
<mi>j</mi>
</msub>
<mo>+</mo>
<msub>
<mi>c</mi>
<mi>i</mi>
</msub>
<msub>
<mi>c</mi>
<mi>j</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>/</mo>
<msqrt>
<mrow>
<msup>
<msub>
<mi>a</mi>
<mi>i</mi>
</msub>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<msub>
<mi>b</mi>
<mi>i</mi>
</msub>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<msub>
<mi>c</mi>
<mi>i</mi>
</msub>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mo>/</mo>
<msqrt>
<mrow>
<msup>
<msub>
<mi>a</mi>
<mi>j</mi>
</msub>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<msub>
<mi>b</mi>
<mi>j</mi>
</msub>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<msub>
<mi>c</mi>
<mi>j</mi>
</msub>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
式中,θij为平面Pi和Pj的夹角,其中i=2,j=1,3;
再按下式计算虚拟旋转角w
w=(θ21+θ23)/2 (4)
3).计算虚拟旋转轴
设O21和O23分别为平面P1和P2、P2和P3的交线,平面P21为过O21的角θ21的平分面,平面P23为过O23的θ23的角平分面,平面P21和P23的交线即为转轴;
通过计算平面向量叉积,分别得到交线O21和O23的法向量n21和n23,算式如下
nij=(bi*cj-ci*bj,ci*aj-ai*cj,ai*bj-bi*aj) (5)
分别计算平面P2和P1、P3的法线夹角角平分线的向量m21和m23,算式如下
mij=((ai+aj)/2,(bi+bj)/2,(ci+cj)/2) (6)
则m21和n21的叉积为平面P21的法向量r21,m23和n23的叉积则为平面P23的法向量r23,算式如下
rij=(mij[2]*nij[3]-mij[3]*nij[2],mij[3]*nij[1]-mij[1]*nij[3],mij[1]*nij[2]-mij[2]*nij[1]) (7)
设角平分面Pij的方程为aijx+bijy+cijz+dij=0,则aij=rij[1],bij=rij[2],cij=rij[3],系统dij的计算公式为
dij=di*T1+dj*T2 (8)
式中,T1和T2的定义分别如公式(9)和公式(10)
T1=(aj*bij-aij*bj)/(aj*bi-ai*bj) (9)
T2=(aij*bi-ai*bij)/(aj*bi-ai*bj) (10)
计算平面P21和P23的法线叉积,就可以得到转轴的方向矢量v,公式为
v=(b21*c23-c21*b23,c21*a23-a21*c23,a21*b23-b21*a23) (11)
4).计算旋转中心点
得到旋转轴的方向矢量后,再计算旋转轴中心点(Xo,Yo,Zo),就可以得旋转轴,计算公式如下
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>X</mi>
<mi>o</mi>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mo>-</mo>
<msub>
<mi>b</mi>
<mn>23</mn>
</msub>
<mo>*</mo>
<msub>
<mi>d</mi>
<mn>21</mn>
</msub>
<mo>+</mo>
<msub>
<mi>b</mi>
<mn>21</mn>
</msub>
<mo>*</mo>
<msub>
<mi>d</mi>
<mn>23</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>a</mi>
<mn>21</mn>
</msub>
<mo>*</mo>
<msub>
<mi>b</mi>
<mn>23</mn>
</msub>
<mo>-</mo>
<msub>
<mi>a</mi>
<mn>23</mn>
</msub>
<mo>*</mo>
<msub>
<mi>b</mi>
<mn>21</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>Y</mi>
<mi>o</mi>
<mo>=</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>a</mi>
<mn>21</mn>
</msub>
<mo>*</mo>
<msub>
<mi>d</mi>
<mn>23</mn>
</msub>
<mo>-</mo>
<msub>
<mi>a</mi>
<mn>23</mn>
</msub>
<mo>*</mo>
<msub>
<mi>d</mi>
<mn>21</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>a</mi>
<mn>21</mn>
</msub>
<mo>*</mo>
<msub>
<mi>b</mi>
<mn>23</mn>
</msub>
<mo>-</mo>
<msub>
<mi>a</mi>
<mn>23</mn>
</msub>
<mo>*</mo>
<msub>
<mi>b</mi>
<mn>21</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>Z</mi>
<mi>o</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>12</mn>
<mo>)</mo>
</mrow>
</mrow>
5).四元数计算
利用上面计算得到的旋转轴矢量和角度,可以进一步转化为四元数quat,计算公式如下
quat=(cosw,v[1]*sinw,v[2]*sinw,v[3]*sinw) (13)
6).测量数据定位融合
计算得到四元数和旋转中心后,转台每次旋转角度β,则每次测量的数据G(xs,ys,zs)(s=1,2,3…n)相对中心点按照公式(14)做平移运算后生成数据G`(xs`,ys`,zs`),之后按照公式(15)运算,最终计算得到的点云数据G``(xs``,ys``,zs``),直接完成数据的拼接;
G`(xs`,ys`,zs`)=G(xs-Xo,ys-Yo,zs-Zo),(s=1,2,3...n) (14)
其中quat`,tn0、tn1、tn2以及tn3的定义如公式(16)、公式(17)所示
quat`=quat*α/β (16)
2.根据权利要求1所述航空发动机叶片型面测量方法,其特征在于,在测量数据定位融合之后,还利用ICP算法,对转台机械误差进行补偿。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610300993.3A CN105737762B (zh) | 2016-05-09 | 2016-05-09 | 一种航空发动机叶片型面测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610300993.3A CN105737762B (zh) | 2016-05-09 | 2016-05-09 | 一种航空发动机叶片型面测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105737762A CN105737762A (zh) | 2016-07-06 |
CN105737762B true CN105737762B (zh) | 2018-03-16 |
Family
ID=56288366
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610300993.3A Active CN105737762B (zh) | 2016-05-09 | 2016-05-09 | 一种航空发动机叶片型面测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105737762B (zh) |
Families Citing this family (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106762762B (zh) * | 2016-12-08 | 2018-05-08 | 沈阳黎明航空发动机(集团)有限责任公司 | 航空发动机导向叶片模拟工作状况的检查系统及方法 |
CN106510099A (zh) * | 2016-12-30 | 2017-03-22 | 哈尔滨福特威尔科技有限公司 | 全自动完整鞋楦三维数据测量装置及方法 |
CN106840026A (zh) * | 2017-01-11 | 2017-06-13 | 江苏科技大学 | 一种基于红外投线仪的三维测量系统及方法 |
CN107270820B (zh) * | 2017-06-05 | 2019-05-17 | 上海交通大学 | 一种大型薄壁构件壁厚在位测量系统和方法 |
CN108267097B (zh) * | 2017-07-17 | 2020-07-10 | 先临三维科技股份有限公司 | 基于双目三维扫描系统的三维重构方法和装置 |
CN107806836A (zh) * | 2017-09-28 | 2018-03-16 | 中国航发动力股份有限公司 | 一种基于五轴联动点光源的精锻叶片型面检测方法 |
CN108151668B (zh) * | 2017-12-15 | 2020-05-22 | 西安交通大学 | 一种叶片型面全数据测量拼接方法及装置 |
CN108917687B (zh) * | 2018-04-26 | 2020-08-04 | 太原理工大学 | 一种航空发动机叶片前后缘微小弧面检测方法 |
CN108844459B (zh) * | 2018-05-03 | 2020-07-03 | 华中科技大学无锡研究院 | 一种叶片数字化样板检测系统的标定方法及装置 |
CN108458659A (zh) * | 2018-05-03 | 2018-08-28 | 华中科技大学无锡研究院 | 一种叶片非接触式检测装置及方法 |
CN109489580B (zh) * | 2018-12-10 | 2021-09-24 | 华东理工大学 | 一种航空发动机叶片表面加工的在机点云检测及补偿方法 |
CN110749290B (zh) * | 2019-10-30 | 2021-06-01 | 易思维(杭州)科技有限公司 | 基于三维投影的特征信息快速定位方法 |
CN111504228A (zh) * | 2020-04-21 | 2020-08-07 | 东方电气集团东方锅炉股份有限公司 | 基于结构光的空气预热器转子三维测量与可视化方法 |
CN111612768B (zh) * | 2020-05-21 | 2023-06-06 | 重庆大学 | 一种采用结构光空间定位和二维工业ct检测叶片方法 |
CN111692997B (zh) * | 2020-06-09 | 2021-08-13 | 西安交通大学 | 一种数据驱动的矢量尾喷口面积在位测量方法 |
CN112013787B (zh) * | 2020-10-21 | 2021-01-26 | 四川大学 | 基于叶片自特征的叶片三维轮廓重建方法 |
CN113375594B (zh) * | 2021-06-08 | 2022-05-06 | 四川大学青岛研究院 | 飞机蒙皮型面数字化检测方法 |
CN113375598A (zh) * | 2021-08-10 | 2021-09-10 | 四川大学 | 基于自基准面的叶片三维型面轮廓高精度匹配方法 |
CN114660259B (zh) * | 2022-03-17 | 2023-04-14 | 中国民用航空飞行学院 | 一种飞机机翼疲劳裂纹检测设备 |
CN114993204B (zh) * | 2022-06-07 | 2023-03-14 | 上海交通大学 | 大尺寸构件轮廓动态测量与辅助特征定位方法及系统 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101373135A (zh) * | 2008-07-01 | 2009-02-25 | 南京航空航天大学 | 基于瞬时随机光照的三维曲面测量装置及方法 |
CN201340256Y (zh) * | 2008-12-09 | 2009-11-04 | 武汉理工大学 | 用于双目立体视觉测量的试验台 |
CN103453849A (zh) * | 2013-07-18 | 2013-12-18 | 黑龙江科技大学 | 多光学传感器协同的复杂曲面零件三维测量方法与系统 |
CN103630085A (zh) * | 2013-10-15 | 2014-03-12 | 广东工业大学 | 多角度叶片测量装置及方法 |
CN104457616A (zh) * | 2014-12-31 | 2015-03-25 | 苏州江奥光电科技有限公司 | 一种360度三维成像的投影装置 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3781438B2 (ja) * | 1993-02-24 | 2006-05-31 | 與語 照明 | 3次元表面形状測定装置 |
-
2016
- 2016-05-09 CN CN201610300993.3A patent/CN105737762B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101373135A (zh) * | 2008-07-01 | 2009-02-25 | 南京航空航天大学 | 基于瞬时随机光照的三维曲面测量装置及方法 |
CN201340256Y (zh) * | 2008-12-09 | 2009-11-04 | 武汉理工大学 | 用于双目立体视觉测量的试验台 |
CN103453849A (zh) * | 2013-07-18 | 2013-12-18 | 黑龙江科技大学 | 多光学传感器协同的复杂曲面零件三维测量方法与系统 |
CN103630085A (zh) * | 2013-10-15 | 2014-03-12 | 广东工业大学 | 多角度叶片测量装置及方法 |
CN104457616A (zh) * | 2014-12-31 | 2015-03-25 | 苏州江奥光电科技有限公司 | 一种360度三维成像的投影装置 |
Also Published As
Publication number | Publication date |
---|---|
CN105737762A (zh) | 2016-07-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105737762B (zh) | 一种航空发动机叶片型面测量方法 | |
CN108921901B (zh) | 一种基于精密二轴转台和激光跟踪仪的大视场相机标定方法 | |
CN108253906B (zh) | 一种桥壳圆度圆柱度检测装置工件轴线定位误差补偿方法 | |
CN103615980B (zh) | 一种板件上圆孔参数的测量方法及系统 | |
CN103729841B (zh) | 一种基于方靶模型和透视投影的相机畸变校正方法 | |
CN106767540B (zh) | 一种交会测量相机光轴与反射镜夹角误差标定方法 | |
JP2004037270A (ja) | キャリブレーション用データ測定装置、測定方法及び測定プログラム、並びにコンピュータ読取可能な記録媒体、画像データ処理装置 | |
CN104034354B (zh) | 一种imu定位定向系统装调工艺 | |
CN106468544B (zh) | 基于光电自准直仪的卫星高精度测角方法 | |
CN103697811B (zh) | 一种相机与结构光源结合获取物体轮廓三维坐标的方法 | |
CN104880215B (zh) | 多传感器平台相机激光一体化检校的方法 | |
CN107121079A (zh) | 一种基于单目视觉的曲面高度信息测量装置及方法 | |
CN106289317B (zh) | 一种单镜头数字航测相机的单机检校方法及装置 | |
CN106949835A (zh) | 医用放疗设备等中心三维坐标精准测量方法 | |
CN102538707B (zh) | 一种对工件进行三维定位的装置及方法 | |
Wu et al. | A calibration method for spatial pose of a laser beam | |
JP4052489B2 (ja) | 平面形状測定解析方法 | |
CN109406525A (zh) | 桥梁表观病害的检测系统及其检测方法 | |
CN109002567A (zh) | 基于实际光学系统装调检测光路的计算机辅助装调方法 | |
CN105931292A (zh) | 一种基于仿射标定的多方向莫尔层析方法 | |
CN115560698A (zh) | 一种基于定位标记的计算机可视化子孔径拼接方法 | |
CN111238422B (zh) | 一种基于小偏摆的三平面参考镜平面度绝对测量方法 | |
CN107085854A (zh) | 一种基于图像处理技术的铁路轨道二维地图的创建方法 | |
Naletto et al. | Innovative optical setup for testing a stereo camera for space applications | |
CN110118530A (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |