CN118009878A - 一种大尺寸多边坐标测量系统的不确定度评价方法、系统、设备及介质 - Google Patents
一种大尺寸多边坐标测量系统的不确定度评价方法、系统、设备及介质 Download PDFInfo
- Publication number
- CN118009878A CN118009878A CN202410079102.0A CN202410079102A CN118009878A CN 118009878 A CN118009878 A CN 118009878A CN 202410079102 A CN202410079102 A CN 202410079102A CN 118009878 A CN118009878 A CN 118009878A
- Authority
- CN
- China
- Prior art keywords
- uncertainty
- coordinate
- measuring
- measuring station
- coordinates
- 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
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
-
- 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/002—Measuring arrangements characterised by the use of optical techniques for measuring two or more coordinates
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明属于外形尺寸测量技术领域,具体涉及一种大尺寸多边坐标测量系统的不确定度评价方法、系统、设备及介质。所述方法,将获取的目标点坐标输入误差传递模型获得坐标值对应的三个坐标不确定度分量,再结合反射镜光学中心偏移引入的不确定度,计算多边坐标测量系统的不确定度系数,进行不确定度评价。本发明从多边坐标测量系统基本原理入手,简化系统坐标系从而缩减参数,可减小计算量,并给出具体数学解析式,有利于测量现场的快速评估,同时便于计量人员学习与运用。
Description
技术领域
本发明属于外形尺寸测量技术领域,具体涉及一种大尺寸多边坐标测量系统的不确定度评价方法、系统、设备及介质。
背景技术
航空航天、轨道交通、船舶以及风能等几米至上百米的大尺寸装备制造的生产质量需要通过高精度大尺寸测量技术进行保证。为满足各类结构大尺寸零部件测量精度需求,目前已发展了多种大尺寸坐标测量系统,包括3坐标机、多边坐标测量系统、激光跟踪绝对测距仪、近景摄影测量系统、室内定位系统(indoor Global Positioning System:iGPS)、室内空间测量定位系统(workspace Measuring and Positioning System:wMPS)等,其中多边坐标测量系统具有最高的理论精度且具有柔性,大范围的优点。
多边坐标测量系统的测量不确定度是指系统对某个被测位置的测量分散程度。非线性、多输入量的系统不确定度评价方法有不确定度传播定律法和蒙特卡洛法两类。不确定度传播定律法是通过求解测量不确定度系数,从而评估各项间输入量对结果的影响,具有方便快捷的优势。现有技术中,如孙威等人在《多边法坐标测量系统中解算方式对测量精度的影响研究》中公开:与同步解算方式相比较,预先准确标定系统参数方式能有效提升测量精度。但孙威等人在技术文献中提及的方法仅考虑了系统参数标定过程一个变量,并未对整个测量系统的不确定度评定方法提出方案,而大尺寸多边坐标测量系统的输出结果是多维变量,因此其测量不确定度评价往往通过耗时较长的蒙特卡洛法来评估。然而,蒙特卡洛法虽然能有效评估不同误差模型的测量方法,但需要数据量非常大且耗时长、成本高。
发明内容
针对现有技术中采用蒙特卡洛法对大尺寸多边坐标测量系统的不确定度进行评价时存在耗时长、成本高的问题,本发明的目的在于提供一种大尺寸多边坐标测量系统的不确定度评价方法、系统、设备及介质,从多边坐标测量系统基本原理入手,简化系统坐标系从而缩减参数,减小计算量,有利于测量现场的快速评估,同时便于计量人员学习与运用。
本发明主要通过以下技术方案实现:一种大尺寸多边坐标测量系统的不确定度评价方法,对由4个台激光跟踪干涉仪及一个反射镜组成的大尺寸多边坐标测量系统进行不确定度评价,一台激光跟踪干涉仪作为一个测站,用于构建三维的坐标系;所述方法将获取的目标点坐标输入误差传递模型获得坐标值对应的三个坐标不确定度分量,再结合反射镜光学中心偏移引入的不确定度,计算多边坐标测量系统的不确定度系数,进行不确定度评价。
为了更好地实现本发明,进一步地,所述方法包括以下步骤:
步骤S1、建立坐标系;
具体是指,先通过四个测站建立三维的坐标系,并记录在该坐标系中各测站的坐标;
步骤S2、获取目标点坐标;
具体是指,使用各测站对应的激光跟踪干涉仪跟踪反射镜,将反射镜移动到目标位置,多边坐标测量系统输出的空间坐标,即为目标点坐标;
步骤S3、获取三个坐标值对应的三个坐标不确定度分量;
具体是指,将目标点坐标输入误差传递模型,计算得到三个坐标值对应的三个坐标不确定度分量;
步骤S4、根据三个坐标值对应的三个坐标不确定度分量、反射镜光学中心偏移引入的不确定度计算多边坐标测量系统的不确定度系数。
为了更好地实现本发明,进一步地,所述步骤S3中误差传递模型根据各测站的坐标、各测站与目标点的实际距离、各测站的测距不确定度计算三个坐标值对应的三个坐标不确定度分量。
为了更好地实现本发明,进一步地,所述步骤S3中各测站与目标点的实际距离通过以下方法获得:先获取反射镜在初始点时各测站与初始点之间的距离;然后,移动反射镜到N个离散的不同的测量点,各测站跟踪反光镜从初始点移动到测量点,并测得各测站与初始点、测量点二者长度之差,记为长度增量dlpk;最后,各测站与初始点的距离与长度增量dlpk;之和即为各测站与目标点的实际距离;所述N大于10。
为了更好地实现本发明,进一步地,所述步骤S3中各测站的测距不确定度根据激光跟踪干涉仪自身的测量不确定度、环境补偿不全所引入的不确定度进行计算。
为了更好地实现本发明,进一步地,所述步骤S4中多边坐标测量系统的不确定度系数通过坐标系中X、Y、Z三轴对应的三个不确定度分量ux、uy、uz以及反射镜光学中心偏移引入的不确定度ub四个参数分别平方后求和再开根号。
为了更好地实现本发明,进一步地,所述反射镜为猫眼逆反射镜。
本发明提供的一种大尺寸多边坐标测量系统的不确定度评价方法,对由4个测站及一个反射镜组成的大尺寸多边坐标测量系统进行坐标测量不确定度评价;所述方法,先通过四个测站建立坐标系并记录各个测站在该坐标系中的坐标;接着,以反射镜的光学中心作为被测点,将反射镜从初始点分别移动到N个不同位置的控制点,通过各测站对应的激光跟踪干涉仪获取测站分别到初始点、N个控制点的相对长度;然后,根据激光跟踪干涉仪自身的测量不确定度、环境补偿不全所引入的不确定度,计算各测站的测距不确定度和反射镜光学中心偏移引入的不确定度;最后,将反射镜放置在待测量位置,获取此时被测点的坐标,由被测点的坐标、各测站在坐标系中的坐标、各测站的测距不确定度和反射镜光学中心偏移引入的不确定度计算多边坐标测量系统的不确定度。
本发明还提供了一种大尺寸多边坐标测量系统的不确定度评价系统,包括误差传递模型,用于在上述的方法中实现获得坐标值对应的三个坐标不确定度分量的作用。
本发明还提供了一种电子设备,包括存储器和处理器;所述存储器上存储有计算机程序;当所述计算机程序在所述处理器上执行时,实现上述的不确定度评价方法。
一种计算机可读存储介质,所述计算机可读存储介质上存储有计算机指令;当所述计算机指令在上述的电子设备上执行时,实现上述的不确定度评价方法。
本发明的有益效果如下:
本发明将获取的目标点坐标输入误差传递模型获得坐标值对应的三个坐标不确定度分量,再结合反射镜光学中心偏移引入的不确定度,计算多边坐标测量系统的不确定度系数,从多边坐标测量系统基本原理入手,简化系统坐标系从而缩减参数,减小计算量,有利于测量现场的快速评估,同时便于计量人员学习与运用。
附图说明
图1为四个测站构建的测站坐标系。
图2为本发明主要步骤流程图。
具体实施方式
实施例1:
本实施例提供了一种大尺寸多边坐标测量系统的不确定度评价方法,对由4个台激光跟踪干涉仪及一个反射镜组成的大尺寸多边坐标测量系统进行不确定度评价,一台激光跟踪干涉仪作为一个测站,用于构建三维的坐标系;所述方法将获取的目标点坐标输入误差传递模型获得坐标值对应的三个坐标不确定度分量,再结合反射镜光学中心偏移引入的不确定度,计算多边坐标测量系统的不确定度系数,进行不确定度评价。
本实施例采用不同于现有技术的不确定度评价方法,从多边坐标测量系统基本原理入手,简化系统坐标系从而缩减参数,减小计算量,有利于测量现场的快速评估,同时便于计量人员学习与运用。
实施例2:
在实施例1的基础上,本实施例中以猫眼逆反射镜为例说明一种大尺寸多边坐标测量系统的不确定度评价方法的详细内容。
所述方法包括以下步骤:
步骤S1、建立坐标系;
具体是指,先通过四个测站建立三维的坐标系,并记录在该坐标系中各测站的坐标;
步骤S2、获取目标点坐标;
具体是指,使用各测站对应的激光跟踪干涉仪跟踪猫眼逆反射镜,将猫眼逆反射镜移动到目标位置,多边坐标测量系统输出的空间坐标,即为目标点坐标;
步骤S3、获取三个坐标值对应的三个坐标不确定度分量;
具体是指,将目标点坐标输入误差传递模型,计算得到三个坐标值对应的三个坐标不确定度分量;
步骤S4、根据三个坐标值对应的三个坐标不确定度分量、猫眼逆反射镜光学中心偏移引入的不确定度计算多边坐标测量系统的不确定度系数。
进一步地,所述步骤S1具体包括:
步骤S11、将4台激光跟踪干涉仪分别固定于不同的4个位置,并将激光跟踪干涉仪按顺序设为测站A、测站B、测站C、测站D;
步骤S12、建立测站对应的三维坐标系,设测站A的测量中心为坐标系原点,由测站A的测量中心指向测站B的方向为X轴的正方向;设测站C的测量中心位于XOY平面上,测站A、测站B、测站C、测站D的坐标分别为:(0,0,0)、(xB,0,0)、(xC,yC,0)、(xD,yD,zD)。
进一步地,所述步骤S3中误差传递模型根据各测站的坐标、各测站与目标点的实际距离、各测站的测距不确定度计算三个坐标值对应的三个坐标不确定度分量。
进一步地,所述步骤S3中各测站与目标点的实际距离通过以下方法获得:先获取反射镜在初始点时各测站与初始点之间的距离;然后,移动反射镜到N个离散的不同的测量点,各测站跟踪反光镜从初始点移动到测量点,并测得各测站与初始点、测量点二者长度之差,记为长度增量dlpk;最后,各测站与初始点的距离与长度增量dlpk;之和即为各测站与目标点的实际距离;所述N大于10。
进一步地,所述步骤S3中各测站的测距不确定度根据激光跟踪干涉仪自身的测量不确定度、环境补偿不全所引入的不确定度进行计算。
进一步地,所述步骤S4中多边坐标测量系统的不确定度系数通过坐标系中X、Y、Z三轴对应的三个不确定度分量ux、uy、uz以及反射镜光学中心偏移引入的不确定度ub四个参数分别平方后求和再开根号。
本实施例所述的不确定度评价方法,首先,通过测站构建坐标系,并获得各测站在坐标系中的坐标;接着,以反射镜的光学中心作为被测点,将反射镜从初始点分别移动到N个不同位置的控制点,通过各测站对应的激光跟踪干涉仪获取测站到初始点及N个控制点的相对长度;然后,按照激光跟踪干涉仪自身的测量不确定度、环境补偿不全所引入的不确定度,计算各测站的测距不确定度和反射镜光学中心偏移引入的不确定度;最后,由各测站在坐标系中的坐标、各测站的测距不确定度和反射镜光学中心偏移引入的不确定度计算多边坐标测量系统的不确定度。
实施例3:
本实施例基于实施例1或实施例2,提供一种具体使用场景。
本实施例采用的多边坐标测量系统由4台激光跟踪干涉仪及猫眼逆反射镜组成,大尺寸多边坐标测量系统的不确定度评价方法包括以下操作:
首先,将4激光跟踪干涉仪分别固定于不同的4个位置,并按任意顺序将激光跟踪干涉仪设为测站A、B、C、D;
其次,建立测站坐标系,设测站A的测量中心为坐标系原点,由测站A的测量中心指向B的方向为X轴的正方向;设测站C的测量中心位于XOY平面上,A、B、C、D的坐标分别为:(0,0,0)、(xB,0,0)、(xC,yC,0)、(xD,yD,zD);
然后,使多边坐标测量系统的激光跟踪干涉仪跟踪猫眼逆反射镜,将猫眼逆反射镜放置待测量位置,多边坐标测量系统的输出为测站坐标系下某个位置的空间坐标x、y、z,按照以下误差传递公式来评估x、y、z标准不确定度ux、uy、uz:
最后,多边坐标测量系统的不确定度系数uG计算方法如式(4)(4)所示:
式(1)(1)(2)、(1)(2)(2)、(1)(3)(2)中所用到的参数有:xB、xC、yC、uA、uB、uC、lPA、lPB、lPC共9个独立参数。
其中:xB为测站B的X轴坐标值,xC为测站C的X轴坐标值,yC为测站C的Y轴坐标值;uA、uB、uC分别为测站A、测站B、测站C的测距不确定度;lPA、lPB、lPC分别为测站A、测站B、测站C到猫眼逆反射镜光学中心的长度;ub为猫眼逆反射镜光学中心偏移引入的不确定度。
在另一具体实施方式中,提供了本实施例中不确定度评估公式推导过程。
现有3台激光跟踪干涉仪:测站A、测站B、测站C;并按照实施例2的步骤S1建立测站坐标系,如图1所示,其中P为目标点,即猫眼逆反射镜在被测位置时猫眼逆反射镜光学中心的坐标。
若已知xB、xC、yC,根据多边坐标测量原理,被测位置对应目标点P的坐标x、y、z的求解公式如下:
式(7)中“±”号可按测量要求取正负。
函数随机误差合成方法如下:设被测量Y的估计值y由n个相互独立的测量值x1,x2,…,xn的函数求得,即:y=f(x1,x2,…,xn);y的函数随机误差合成为uG的计算公式为:
其中,δx1,δx2,…,δxn分别为测量值x1,x2,…,xn的误差,函数f对变量xi的偏微分。
式(8)(5)(7)中的x、y、z分别与独立的测量值lPA、lPB、lPC成函数关系,以下根据多边法测量结果中x、y、z方向的坐标误差极限δx、δy、δz来表征不确定度ux,uy,uz:
其中:δA、δB、δC分别为测站A、测站B、测站C的激光测距误差。
通过极限误差来表征不确定度,则有式(1)(12)(3)-(14)。
在另一具体实施方式中,提供了所用到xB、xC、yC、lPA、lPB、lPC、x、y、z、uA、uB、uC共各个独立参数的测量计算方法。
一、获取xB、xC、yC。
xB、xC、yC是测站B的X轴坐标值、测站C的X轴坐标值、测站C的Y轴坐标值,直接读取即可。
二、获取lPA、lPB、lPC。
激光跟踪干涉仪直接得到的测量值是相对距离,一般按以下方式来表示:定义一个初始点P0,设lA0、lB0、lC0、lD0分别为测站A、测站B、测站C、测站D至初始点P0的距离。当某测站k跟踪猫眼逆反射镜从初始点P0移动到测量点P时,可测得测站k分别与P0、P连线之间的长度之差dlPk(k=A,B,C,D);k=A,B,C,D。dlPk本质就是增量,所以dlPk与lk0的和为测站k到测量点P的距离。
所以lPA、lPB、lPC及lPD一般表示为:
其中,dlPk是激光跟踪干涉仪直接测量值。
三、获取目标点P的坐标x、y、z。
计算被测位置坐标的推导过程中需额外引入测站D到目标点的距离lPD及坐标xD、yD、zD,但实际计算时仅需使用到xB、xC、yC、lPA、lPB、lPC 6个参数即可计算被测位置坐标。具体地,将获得的xB、xC、yC、lPA、lPB、lPC 6个参数带入式(5)、式(6)、式(7),获得目标点P的坐标x、y、z。
四、解算控制点最小值。
需要说明的时,xB、xC、yC、xD、yD、zD、lA0、lB0、lC0、lD0需要通过冗余信息构建方程,再进行解算;解算方法如下。
(1)将4台激光跟踪干涉仪跟踪同一猫眼逆反射镜,并将猫眼反射镜放置到P0进行测量,并将猫眼反射镜分别移动至相对分散的N个位置、测得N个控制点相对于初始点的相对长度,其中N不小于20,第i个控制点对应的方程如下:
其中,fi(δ)为第i个控制点到测站D距离的估算值与测量值的差。
xi、yi分别为第i个控制点的X轴、Y轴坐标,且 δ为(xB,xC,yC,xD,yD,zD,lA0,lB0,lC0,lD0)。
(2)通过卷尺等长度测量器具测量4台激光跟踪干涉仪的位置间距以及4台激光跟踪干涉仪分别到初始点P0的距离从而获得待求参数的迭代初始值。
(3)通过Matlab软件,将N个控制点的方程的平方和最小化进行求解,令:
函数φ(δ)的极小点δ*就是δ的解,即:
为求δ*由多元函数φ(δ)取极值的必要条件,令φ(δ)的梯度函数g(δ)=▽φ(δ)=0即:
g(δ)=▽φ(δ)=0.5×Dfi T(x)×fi(x)=0 (19)
其中,记为J(δ)。
δ有解时,N至少大于10,即控制点的数量需大于10。
五、获取uA、uB、uC。
测站A对应的测长不确定度uA、测站B对应的测长不确定度uB以及测站C对应的测长不确定度uC均可以按以下方法进行评估。用测站k表示测站A、测站B、测站C中任意一个测站。
(1)本发明采用Etalon的激光跟踪干涉仪,正态分布95%置信区间,测距不确定度U=0.2μm+0.3μm/m,激光跟踪干涉仪测距的标准不确定度为:
uIFM=0.1+0.15×10-3/l (20)
(2)由于温度、湿度及气压等随机变化的环境因素引起的空气折射率变化,从而造成激光跟踪干涉仪测距误差。按照空气温度梯度为0.5℃/m,环境补偿不全所引入的标准不确定度:
(3)测距的标准不确定度合成:
其中,l为测站到被测点的距离(单位:m)。
六、获取ub。
本实施例中采用猫眼逆反射镜的光学中心误差为±4μm,假设服从均匀分布,猫眼逆反射镜光学中心偏移引入的标准不确定度计算公式如下:
已知猫眼逆反射镜的光学中心误差的情况下,ub可以直接解算。
将获得的获取各测站的测长不确定度uA、uB、uC以及猫眼逆反射镜光学中心偏移引入的标准不确定度ub,即可获得多边坐标测量系统的不确定度系数uG。
实施例4:
本实施例在实施例3的基础上,提供一种验证方法。
本实施例通过估算坐标误差极限以作为本发明的验算复核方法,具体方法以下进行说明及解释。
设目标点P分别与测站A、测站B、测站C的理论距离为lPA、lPB、lPC,并设距离测量值分别为lPA+δlPAi、lPB+δlPBi、lPC+δlPCi,其中δlPAi、δlPBi、δlPCi为各测站对P点激光测距的误差;i=1,2,…,n。
设目标点P在测站坐标系下的理论坐标为x0、y0、z0,实际测量值为xi、yi、zi。
设空间坐标理论值x0、y0、z0分别关于激光测距理论值的函数关系为:
其中l=(lPA,lPB,lPC)T。
设空间坐标测量值xi、yi、zi分别关于激光测距理论值的函数关系为:
其中δl=(δlPAi,δlPBi,δlPCi)T。
令(δlPAi,δlPBi,δlPCi)分别取表1中的值,然后计算核验值M核验=((x0-xi)2+(y0-yi)2+(z0-zi)2)0.5,不考虑ub的前提下,计算结果中的最大值约等于不确定度uG。由于uG是根据误差极限来计算的所以,uG大于M核验,但不超出1μm(0.001mm)。
表1测距误差取值
具体参数举例说明。
步骤a.先按照本发明的方法对某点坐标的不确定度估算。
设测站A、B、C以及被测点P的坐标值见表2。
表2A、B、C以及被测点P的坐标值(mm)
然后可求得目标点P分别至测站A、B、C的理论距离值lPA、lPB、lPC见表3。
表3P分别至A、B、C的理论距离值(mm)
根据激光干涉测距的不确定度U=0.1μm+0.15μm/m,有uA、uB、uC,见表4。
表4 3测站对P点的测距不确定度(mm)
根据实施例3中式(1)、(2)、(3),进而估算出ux、uy、uz,见表5。
表5x、y、z方向的不确定度
| ux | 0.0020 |
| uy | 0.0030 |
| uz | 0.0049 |
不计入ub时坐标不确定度uG估算为:步骤b.然后根据前文的验核方法进行验核。
首先计算x0、y0、z0以及xi、yi、zi。
计算得到(x0,y0,z0)=(2000mm,4000mm,2000mm)。
按照表1设置的测距误差有:
表6模拟测距误差值组合(mm)
然后计算xi、yi、zi。
令:
最后按照8中测距误差组合计算核验值M核验=((x0-xi)2+(y0-yi)2+(z0-zi)2)0.5。
表7 8中测距误差组合计算核验值(mm)
M核验最大值0.0059mm与步骤a.中计算的不确定度估算值uG=0.0061mm仅相差0.0002mm低于0.001mm(即1μm)。
本实施例中所述验证方法的理由说明如下:
因为x0-xi,y0-yi,z0-zi有以下关系:
令δx=x0-xi,δy=y0-yi,δz=z0-zi,并且令δx=(δx,δy,δz),则有:
根据一阶泰勒展开公式,有:
δx=M×δl (9)
其中
对式(6)取2范数,有:
||δx||2=||M×δl||2 (10)根据2范数的相容性这一特点,有不等式:
||δx||2=||M×δl||2≤||M||2×||δl||2 (11)即有:
||δx||2=(δx2+δy2+δz2)0.5≤||M||2×||δl||2 (12)
即P的理论坐标与测量坐标点的间距所表征的误差最大不超出||M||2×||δl||2。其中M测站位置和测量点相对测站的位置确定了;而||δl||2取最大值(即测距不确定度数值uA、uB、uC)时可估算出坐标误差的最大值,可逼近极限误差。
实施例5:
本实施例在上述实施例1-实施例3任一项的基础上,提出一种大尺寸多边坐标测量系统的不确定度评价系统,包括误差传递模型,用于在上述的方法中实现获得坐标值对应的三个坐标不确定度分量的作用。
实施例:6:
本发明还提供了一种电子设备,包括存储器和处理器;所述存储器上存储有计算机程序;当所述计算机程序在所述处理器上执行时,实现上述的不确定度评价方法。
实施例7:
一种计算机可读存储介质,所述计算机可读存储介质上存储有计算机指令;当所述计算机指令在上述的电子设备上执行时,实现上述的不确定度评价方法。
以上所述,仅是本发明的较佳实施例,并非对本发明做任何形式上的限制,凡是依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化,均落入本发明的保护范围之内。
Claims (10)
1.一种大尺寸多边坐标测量系统的不确定度评价方法,对由4个台激光跟踪干涉仪及一个反射镜组成的大尺寸多边坐标测量系统进行不确定度评价,一台激光跟踪干涉仪作为一个测站,用于构建三维的坐标系;其特征在于:所述方法将获取的目标点坐标输入误差传递模型获得坐标值对应的三个坐标不确定度分量,再结合反射镜光学中心偏移引入的不确定度,计算多边坐标测量系统的不确定度系数,进行不确定度评价。
2.根据权利要求1所述的一种线缆屏蔽效能测量方法,其特征在于,所述方法包括以下步骤:
步骤S1、建立坐标系;
具体是指,先通过四个测站建立三维的坐标系,并记录在该坐标系中各测站的坐标;
步骤S2、获取目标点坐标;
具体是指,使用各测站对应的激光跟踪干涉仪跟踪反射镜,将反射镜移动到目标位置,多边坐标测量系统输出的空间坐标,即为目标点坐标;
步骤S3、获取三个坐标值对应的三个坐标不确定度分量;
具体是指,将目标点坐标输入误差传递模型,计算得到三个坐标值对应的三个坐标不确定度分量;
步骤S4、根据三个坐标值对应的三个坐标不确定度分量、反射镜光学中心偏移引入的不确定度计算多边坐标测量系统的不确定度系数。
3.根据权利要求2所述的一种线缆屏蔽效能测量方法,其特征在于,所述步骤S3中误差传递模型根据各测站的坐标、各测站与目标点的实际距离、各测站的测距不确定度计算三个坐标值对应的三个坐标不确定度分量。
4.根据权利要求3所述的一种线缆屏蔽效能测量方法,其特征在于,所述步骤S3中各测站与目标点的实际距离通过以下方法获得:先获取反射镜在初始点时各测站与初始点之间的距离;然后,移动反射镜到N个离散的不同的测量点,各测站跟踪反光镜从初始点移动到测量点,并测得各测站与初始点、测量点二者长度之差,记为长度增量dlpk;最后,各测站与初始点的距离与长度增量dlpk;之和即为各测站与目标点的实际距离;所述N大于10。
5.根据权利要求3所述的一种线缆屏蔽效能测量方法,其特征在于,所述步骤S3中各测站的测距不确定度根据激光跟踪干涉仪自身的测量不确定度、环境补偿不全所引入的不确定度进行计算。
6.根据权利要求3所述的一种线缆屏蔽效能测量方法,其特征在于,所述步骤S4中多边坐标测量系统的不确定度系数通过坐标系中X、Y、Z三轴对应的三个不确定度分量ux、uy、uz以及反射镜光学中心偏移引入的不确定度ub四个参数分别平方后求和再开根号。
7.根据权利要求3-6任一项所述的一种线缆屏蔽效能测量方法,其特征在于,所述反射镜为猫眼逆反射镜。
8.一种大尺寸多边坐标测量系统的不确定度评价系统,其特征在于,包括误差传递模型,用于在如权利要求1-7任一项所述的方法中实现获得坐标值对应的三个坐标不确定度分量的作用。
9.一种电子设备,其特征在于,包括存储器和处理器;所述存储器上存储有计算机程序;当所述计算机程序在所述处理器上执行时,实现如权利要求1-7任一项所述的不确定度评价方法。
10.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质上存储有计算机指令;当所述计算机指令在如权利要求9所述的电子设备上执行时,实现如权利要求1-7任一项所述的不确定度评价方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410079102.0A CN118009878B (zh) | 2024-01-19 | 2024-01-19 | 一种大尺寸多边坐标测量系统的不确定度评价方法、系统、设备及介质 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410079102.0A CN118009878B (zh) | 2024-01-19 | 2024-01-19 | 一种大尺寸多边坐标测量系统的不确定度评价方法、系统、设备及介质 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN118009878A true CN118009878A (zh) | 2024-05-10 |
| CN118009878B CN118009878B (zh) | 2025-09-09 |
Family
ID=90955112
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202410079102.0A Active CN118009878B (zh) | 2024-01-19 | 2024-01-19 | 一种大尺寸多边坐标测量系统的不确定度评价方法、系统、设备及介质 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN118009878B (zh) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4884889A (en) * | 1987-11-19 | 1989-12-05 | Brown & Sharpe Manufacturing Company | Calibration system for coordinate measuring machine |
| CN104374317A (zh) * | 2014-11-06 | 2015-02-25 | 北京工业大学 | 一种基于激光跟踪仪多点测量技术的机床误差标定方法 |
| CN108180831A (zh) * | 2017-12-30 | 2018-06-19 | 北京工业大学 | 基于lt多站位测量的cmm坐标误差修正系统不确定度分析方法 |
-
2024
- 2024-01-19 CN CN202410079102.0A patent/CN118009878B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4884889A (en) * | 1987-11-19 | 1989-12-05 | Brown & Sharpe Manufacturing Company | Calibration system for coordinate measuring machine |
| CN104374317A (zh) * | 2014-11-06 | 2015-02-25 | 北京工业大学 | 一种基于激光跟踪仪多点测量技术的机床误差标定方法 |
| CN108180831A (zh) * | 2017-12-30 | 2018-06-19 | 北京工业大学 | 基于lt多站位测量的cmm坐标误差修正系统不确定度分析方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN118009878B (zh) | 2025-09-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111238454A (zh) | 一种空间三维坐标测量仪系统及其三维控制网测量方法 | |
| CN104596420B (zh) | 激光跟踪仪测量基准立方镜中心点位置的精测方法 | |
| CN114239170B (zh) | 一种飞机装配的误差补偿方法 | |
| Zha et al. | Volumetric error compensation of machine tool using laser tracer and machining verification | |
| CN109304730A (zh) | 一种基于激光测距仪的机器人运动学参数标定方法 | |
| CN102589571B (zh) | 一种空间三维视觉计算验证方法 | |
| CN110516350B (zh) | 一种基于各向异性加权的ers点误差修正方法 | |
| WO2021098808A1 (zh) | 激光跟踪仪站位确定方法、系统、电子设备及介质 | |
| Acero et al. | Verification of an articulated arm coordinate measuring machine using a laser tracker as reference equipment and an indexed metrology platform | |
| Li et al. | A novel orientation and position measuring system for large & medium scale precision assembly | |
| CN106052556A (zh) | 一种基于激光追踪仪多站位测量的三坐标测量机空域坐标修正方法 | |
| CN110926366A (zh) | 一种基于激光跟踪仪多站布局的曲面轮廓测量方法 | |
| CN113587819A (zh) | 一种基于三路激光跟踪的大尺度空间位姿动态测量方法与测量精度验证方法 | |
| CN110779503B (zh) | 一种三维精密控制网测量方法 | |
| CN114485392A (zh) | 基于激光跟踪仪的航天器大尺寸机械基准建立方法和系统 | |
| CN103737433A (zh) | 大尺寸构件复合式精度测量方法 | |
| CN114754672B (zh) | 一种机器人与测量仪器的坐标转换关系确定方法 | |
| CN107588929B (zh) | 球幕投影/跟踪系统标定方法及标定器 | |
| CN116021341B (zh) | 激光跟踪干涉仪测量场在机床坐标系中自标定方法及系统 | |
| Ding et al. | A common points selection method based on the uniformity of FLTMMS for pre-alignment of particle accelerators | |
| CN118009878A (zh) | 一种大尺寸多边坐标测量系统的不确定度评价方法、系统、设备及介质 | |
| CN115682982A (zh) | 基于可变大尺寸多面体校准三维测量结果的系统及方法 | |
| GB2639760A (en) | Digital pre-assembly fabrication control method for long-span arch bridge | |
| CN114034244B (zh) | 一种基于Makima算法获得CMM测量范围内任意点几何误差的方法 | |
| CN116877347B (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 |