CN113310445B - 一种多仪器组合测量***的校准方法 - Google Patents
一种多仪器组合测量***的校准方法 Download PDFInfo
- Publication number
- CN113310445B CN113310445B CN202110576171.9A CN202110576171A CN113310445B CN 113310445 B CN113310445 B CN 113310445B CN 202110576171 A CN202110576171 A CN 202110576171A CN 113310445 B CN113310445 B CN 113310445B
- Authority
- CN
- China
- Prior art keywords
- measuring
- instrument
- length
- formula
- instrument combined
- 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
- 238000005259 measurement Methods 0.000 title claims abstract description 96
- 238000000034 method Methods 0.000 title claims abstract description 21
- 238000013519 translation Methods 0.000 claims abstract description 23
- 239000013598 vector Substances 0.000 claims description 15
- 239000011159 matrix material Substances 0.000 claims description 11
- CLOMYZFHNHFSIQ-UHFFFAOYSA-N clonixin Chemical compound CC1=C(Cl)C=CC=C1NC1=NC=CC=C1C(O)=O CLOMYZFHNHFSIQ-UHFFFAOYSA-N 0.000 claims description 6
- 238000006243 chemical reaction Methods 0.000 claims description 4
- 230000004927 fusion Effects 0.000 claims description 4
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 2
- BJQHLKABXJIVAM-UHFFFAOYSA-N bis(2-ethylhexyl) phthalate Chemical compound CCCCC(CC)COC(=O)C1=CC=CC=C1C(=O)OCC(CC)CCCC BJQHLKABXJIVAM-UHFFFAOYSA-N 0.000 abstract description 2
- 238000011161 development Methods 0.000 abstract description 2
- 238000011160 research Methods 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 4
- 238000012360 testing method Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B21/00—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B21/00—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant
- G01B21/02—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant for measuring length, width, or thickness
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Length Measuring Devices With Unspecified Measuring Means (AREA)
Abstract
本发明涉及一种多仪器组合测量***的校准方法,采用多齿分度台作为标准器校准多仪器组合测量***的坐标系旋转误差;采用带有激光干涉仪的直线导轨作为标准器校准多仪器组合测量***的坐标系平移误差;采用长度标尺作为标准器校准多仪器组合测量***的公共测量点和非公共测量点的空间距离长度误差。本发明弥补了现有单台测量仪器的校准方法不能直接应用于多仪器组合测量***的校准的问题,保障了多仪器组合测量***的量值统一、准确可靠,助力多仪器组合测量***的研发升级和现场应用。
Description
技术领域
本发明涉及一种多仪器组合测量***的校准方法,是一种针对多仪器组合测量***的外参数及坐标测量的准确度的校准方法,属于计量测试领域。
背景技术
面对航空航天、船舶等大型装备制造中大空间、高精度、多信息、低成本以及高效可靠的测量需求,由多台或多种测量仪器组成的多仪器组合测量***成为重要测量手段,得到广泛应用。如图1所示,多仪器组合测量***中,各测量仪器在各自的仪器坐标系下进行测量,多台测量仪器通过一组测量公共测量点计算其外参数(坐标系旋转、平移转换参数)以将各仪器坐标***一至***坐标系;公共测量点由两台及以上测量仪器测得,其坐标由这些仪器的测量数据融合计算获得;非公共测量点仅由一台测量仪器测量测得,其坐标由该仪器的测量数据经该仪器的外参转换至***坐标系下获得。综上,多仪器组合测量***中不可避免地涉及各测量仪器的坐标***一、公共测量点的数据融合,以及非公共测量点的数据转换的问题,为确保多仪器组合测量***的量值统一、准确可靠,其外参数及坐标测量的测试不能直接套用单台测量仪器的校准方法(如《JJF 1242激光跟踪三维坐标测量***校准规范》、《JJF 1408关节臂式坐标测量机校准规范》等),目前现有技术中多仪器组合测量***的整体校准问题尚未解决。
发明内容
本发明为弥补现有单台测量仪器的校准方法不能直接应用于多仪器组合测量***的校准的问题,提供一种多仪器组合测量***的校准方法,通过分别测试多仪器组合测量***的外参数、公共测量点坐标及非公共测量点坐标,确保多仪器组合测量***的量值统一、准确可靠。
本发明采取以下技术方案:
一种多仪器组合测量***的校准方法,采用多齿分度台作为标准器校准多仪器组合测量***的坐标系旋转误差;采用带有激光干涉仪的直线导轨作为标准器校准多仪器组合测量***的坐标系平移误差;采用长度标尺作为标准器校准多仪器组合测量***的公共测量点和非公共测量点的空间距离长度误差。
优选的,多仪器组合测量***的坐标系旋转误差校准包括以下步骤:
步骤一、由N台测量仪器组成的多仪器组合测量***,将其中1台测量仪器固定在多齿分度台上,其余N-1台测量仪器固定在多齿分度台旁,N台测量仪器的整体布局与其现场应用时一致;
步骤二、多齿分度台刻度线归零,标定***中N台测量仪器的外参数,将多齿分度台上的测量仪器的坐标系旋转矩阵记作R0;
步骤三、多齿分度台旋转角度θ,记录此时多齿分度台示值θr1;
步骤四、标定***中N台测量仪器的外参数,重复三遍,将多齿分度台上的测量仪器的坐标系旋转矩阵记作R11、R12、R13;
步骤五、分别计算坐标系旋转矩阵R11、R12、R13与R0之间的变化量ΔR11、ΔR12、ΔR13,如式(1)所示;结合罗德里格斯公式计算对应的旋转角θm11、θm12、θm13,如式(2)所示;并计算与θr1的差值Δθ11、Δθ12、Δθ13,如式(3)所示;取最大值记作该角度处的坐标系旋转误差Δθ1,如式(4)所示;
其中,r1js表示矩阵ΔR1j中第s个元素,j=1,2,3;
θm1j=arcsin(u2+v2+w2) (2);
Δθ1j=|θm1j-θr1| (3);
其中,j=1,2,3;
Δθ1=max{Δθ11,Δθ12,Δθ13} (4);
步骤六、按照步进θ,依次旋转多齿分度台至示值θr2、θr3、…、θrM处,重复步骤四、步骤五,测得各处的坐标系旋转误差Δθ2、Δθ3、…、ΔθM,并取最大值作为本次校准的坐标系旋转误差Δθ,如式(5)所示;
Δθ=max{Δθ1,Δθ2,...,ΔθM} (5)。
优选的,多仪器组合测量***的坐标系平移误差校准包括以下步骤:
步骤一、由N台测量仪器组成的多仪器组合测量***,将其中1台测量仪器固定在直线导轨的移动工作台上,其余N-1台测量仪器固定在直线导轨旁,N台测量仪器的整体布局与其现场应用时一致;
步骤二、激光干涉仪复位,标定***中N台测量仪器的外参数,将移动工作台上的测量仪器的坐标系平移向量记作T0;
步骤三、移动工作台沿直线导轨移动位移d,记录此时激光干涉仪器的示值dr1;
步骤四、标定***中N台测量仪器的外参数,重复三遍,将移动工作台上的测量仪器的坐标系平移向量记作T11、T12、T13;
步骤五、分别计算坐标系平移向量T11、T12、T13与T0之间的变化量ΔT11、ΔT12、ΔT13,如式(6)所示,结合向量模的计算公式计算对应的平移量dm11、dm12、dm13,如式(7)所示,并计算与dr1的差值Δd11、Δd12、Δd13,如式(8)所示,取最大值记作该位置处的坐标系平移误差Δd1,如式(9)所示;
其中,tx1j、ty1j、tz1j表示向量ΔT1j中3个元素,j=1,2,3;
其中,j=1,2,3;
Δd1j=|dm1j-dr1| (8);
其中,j=1,2,3;
Δd1=max{Δd11,Δd12,Δd13} (9);
步骤六、按照步进d,依次沿直线导轨将移动工作台移至示值dr2、dr3、…、drM处,重复步骤四、步骤五,测得各处的坐标系平移误差Δd2、Δd3、…、ΔdM,并取最大值作为本次校准的坐标系平移误差Δd,如式(10)所示;
Δd=max{Δd1,Δd2,...,ΔdM} (10)。
优选的,多仪器组合测量***的公共测量点的空间长度示值误差校准包括以下步骤:
步骤一、由N台测量仪器组成的多仪器组合测量***,按照与现场应用时一致的整体布局放置,在多仪器组合测量***的测量范围内布置M根长度标尺,每根长度标尺的两端各安装有一个靶标座,可适配多仪器组合测量***中的不同测量仪器的靶标,且保持测量点不随靶标类型改变而变动,长度标尺的方向包括水平、垂直、倾斜等不同方向;
步骤二、保证M根长度标尺两端的测量点均可被多仪器组合测量***中2台及以上的测量仪器测得,将长度标尺两端的测量点作为公共测量点,配合其他公共测量点,标定多仪器组合测量***的外参数并通过数据融合计算公共测量点的坐标;重复三遍,将第j遍中第i根长度标尺两端的公共测量点坐标分别记作(xi1j,yi1j,zi1j)和(xi2j,yi2j,zi2j),i=1,2,…,M,j=1,2,3;
步骤三、利用坐标(xi1j,yi1j,zi1j)、(xi2j,yi2j,zi2j),计算第1根长度标尺的长度lm1j,如式(11)所示,并计算其与该长度标尺的参考值lr1j的差值Δl1j,如式(12)所示;取Δl11、Δl12、Δl13中最大值作为第1根长度标尺处的空间长度示值误差Δl1,如式(13)所示;
其中,j=1,2,3;
Δl1j=|lm1j-lr1| (12);
其中,j=1,2,3;
Δl1=max{Δl11,Δl12,Δl13} (13);
步骤四、重复步骤三,依次计算剩余M-1根长度标尺处的空间长度示值误差Δli,i=2,3,…,M;取M根长度标尺处的空间长度示值误差的最大值作为本次校准的公共测量点的空间长度示值误差Δl,如式(14)所示;
Δl=max{Δl1,Δl2,...,ΔlM} (14)。
进一步的,多仪器组合测量***的非公共测量点的空间长度示值误差校准包括以下步骤:
步骤一、由N台测量仪器组成的多仪器组合测量***,按照与现场应用时一致的整体布局放置,标定多仪器组合测量***的外参数;
步骤二、在多仪器组合测量***的测量范围内布置M根长度标尺,每根长度标尺的两端各安装有一个靶标座,可适配多仪器组合测量***中的不同测量仪器的靶标,且保持测量点不随靶标类型改变而变动,长度标尺的方向包括水平、垂直、倾斜方向;
步骤三、多仪器组合测量***对M根长度标尺进行测量,保证每根长度标尺的测量点均仅被一台测量仪器测量,且同一根长度标尺两端的测量点不被同一台测量仪器测量,即M根长度标尺上的测量点均为非公共测量点;根据多仪器组合测量***的外参数完成数据转换,获得非公共测量点的坐标;重复三遍;将第j遍中第i根长度标尺两端的公共测量点坐标分别记作(xi1j,yi1j,zi1j)和(xi2j,yi2j,zi2j),i=1,2,…,M,j=1,2,3;
步骤四、利用坐标(xi1j,yi1j,zi1j)、(xi2j,yi2j,zi2j),计算第1根长度标尺的长度lm1j,如式(11)所示,并计算其与该长度标尺的参考值lr1j的差值Δl1j,如式(12)所示;取Δl11、Δl12、Δl13中最大值作为第1根长度标尺处的空间长度示值误差Δl1,如式(13)所示;
步骤五、重复步骤四,依次计算剩余M-1根长度标尺处的空间长度示值误差Δli,i=2,3,…,M。取M根长度标尺处的空间长度示值误差的最大值作为本次校准的非公共测量点的空间长度示值误差Δl,如式(14)所示。
本发明的有益效果在于:针对外参数、公共测量点及非公共测量点的不同测试需求,提供一种***的、可溯源的多仪器组合测量***校准方法,弥补了现有单台测量仪器的校准方法不能直接应用于多仪器组合测量***的校准的问题,保障了多仪器组合测量***的量值统一、准确可靠,助力多仪器组合测量***的研发升级和现场应用。
附图说明
图1是多仪器组合测量***示意图。
图2是多仪器组合测量***的坐标系旋转误差校准示意图。
图3是多仪器组合测量***的坐标系平移误差校准示意图。
图4是多仪器组合测量***的公共测量点空间长度示值误差校准示意图。
图5是多仪器组合测量***的非公共测量点空间长度示值误差校准示意图。
具体实施方式
下面结合实施例和附图对本发明的优选实施方式进行详细说明。
实施例一:
以由2台测量仪器组成的多仪器组合测量***为例,多仪器组合测量***的坐标系旋转误差校准如图2所示,包括以下步骤:
步骤一、将第一测量仪器固定在多齿分度台上,第二测量仪器固定在多齿分度台旁,2台测量仪器的整体布局与其现场应用时一致。
步骤二、多齿分度台刻度线归零,按照多仪器组合测量***的操作手册标定***中2台测量仪器的外参数,将多齿分度台上的测量仪器的坐标系旋转矩阵记作R0。
步骤三、多齿分度台旋转角度θ,记录此时多齿分度台示值θr1。
步骤四、按照多仪器组合测量***的操作手册标定***中2台测量仪器的外参数,重复三遍,将多齿分度台上的测量仪器的坐标系旋转矩阵记作R11、R12、R13。
步骤五、分别计算坐标系旋转矩阵R11、R12、R13与R0之间的变化量ΔR11、ΔR12、ΔR13,如式(15)所示,结合罗德里格斯公式计算对应的旋转角θm11、
θm12、θm13,如式(16)所示,并计算与θr1的差值Δθ11、Δθ12、Δθ13,如式(17)所示,取最大值记作该角度处的坐标系旋转误差Δθ1,如式(18)所示。
其中,r1js表示矩阵ΔR1j中第s个元素,j=1,2,3。
θm1j=arcsin(u2+v2+w2) (16)
Δθ1j=|θm1j-θr1| (17)
其中,j=1,2,3。
Δθ1=max{Δθ11,Δθ12,Δθ13} (18)
步骤六、按照步进θ,依次旋转多齿分度台至示值θr2、θr3、…、θrM处,重复步骤四、步骤五,测得各处的坐标系旋转误差Δθr2、Δθr3、…、ΔθrM,并取最大值作为本次校准的坐标系旋转误差Δθ,如式(19)所示。
Δθ=max{Δθ1,Δθ2,...,ΔθM} (19)。
实施例二:
以由2台测量仪器组成的多仪器组合测量***为例,多仪器组合测量***的坐标系平移误差校准如图3所示,包括以下步骤:
步骤一、将第一测量仪器固定在直线导轨的移动工作台上,第二测量仪器固定在直线导轨旁,2台测量仪器的整体布局与其现场应用时一致。
步骤二、激光干涉仪复位,按照多仪器组合测量***的操作手册标定***中2台测量仪器的外参数,将移动工作台上的测量仪器的坐标系平移向量记作T0。
步骤三、移动工作台沿直线导轨移动位移d,记录此时激光干涉仪的示值dr1。
步骤四、按照多仪器组合测量***的操作手册标定***中2台测量仪器的外参数,重复三遍,将移动工作台上的测量仪器的坐标系平移向量记作T11、T12、T13。
步骤五、分别计算坐标系平移向量T11、T12、T13与T0之间的变化量ΔT11、ΔT12、ΔT13,如式(20)所示,结合向量模的计算公式计算对应的平移量dm11、dm12、dm13,如式(21)所示,并计算与dr1的差值Δd11、Δd12、Δd13,如式(22)所示,取最大值记作该位置处的坐标系平移误差Δd1,如式(23)所示。
其中,tx1j、ty1j、tz1j表示向量ΔT1j中3个元素,j=1,2,3。
其中,j=1,2,3。
Δd1j=|dm1j-dr1| (22)
其中,j=1,2,3。
Δd1=max{Δd11,Δd12,Δd13} (23)
步骤六、按照步进d,依次沿直线导轨将移动工作台移至示值dr2、dr3、…、drM处,重复步骤四、步骤五,测得各处的坐标系平移误差Δdr2、Δdr3、…、ΔdrM,并取最大值作为本次校准的坐标系平移误差Δd,如式(24)所示。
Δd=max{Δd1,Δd2,...,ΔdM} (24)
实施例三:
以由2台测量仪器组成的多仪器组合测量***配合3根长度标尺为例,多仪器组合测量***的公共测量点的空间长度示值误差校准如图4所示,包括以下步骤:
步骤一、由2台测量仪器组成的多仪器组合测量***,按照与现场应用时一致的整体布局放置,在多仪器组合测量***的测量范围内布置3根长度标尺,每根长度标尺的两端各安装有一个靶标座,可适配多仪器组合测量***中的不同测量仪器的靶标,且保持测量点不随靶标类型改变而变动,长度标尺的方向包括水平、垂直、倾斜等不同方向。
步骤二、保证3根长度标尺两端的测量点均可被多仪器组合测量***中2台的测量仪器测得,将长度标尺两端的测量点作为公共测量点,配合其他公共测量点,按照多仪器组合测量***的操作手册标定多仪器组合测量***的外参数并通过数据融合计算公共测量点的坐标。重复三遍。将第j遍中第i根长度标尺两端的公共测量点坐标分别记作(xi1j,yi1j,zi1j)和(xi2j,yi2j,zi2j),i=1,2,3,j=1,2,3。
步骤三、利用坐标(x11j,y11j,z11j)、(x12j,y12j,z12j),计算第1根长度标尺的长度lm1j,如式(25)所示,并计算其与该长度标尺的参考值lr1j的差值Δl1j,如式(26)所示。取Δl11、Δl12、Δl13中最大值作为第1根长度标尺处的空间长度示值误差Δl1,如式(27)所示。
其中,j=1,2,3。
Δl1j=|lm1j-lr1| (26)
其中,j=1,2,3。
Δl1=max{Δl11,Δl12,Δl13} (27)
步骤四、重复步骤三,依次计算剩余2根长度标尺处的空间长度示值误差Δli,i=2,3。取3根长度标尺处的空间长度示值误差的最大值作为本次校准的公共测量点的空间长度示值误差Δl,如式(28)所示。
Δl=max{Δl1,Δl2,Δl3} (28)。
实施例四:
以由2台测量仪器组成的多仪器组合测量***配合3根长度标尺为例,多仪器组合测量***的非公共测量点的空间长度示值误差校准如图5所示,包括以下步骤:
步骤一、由2台测量仪器组成的多仪器组合测量***,按照与现场应用时一致的整体布局放置,按照多仪器组合测量***的操作手册标定多仪器组合测量***的外参数。
步骤二、在多仪器组合测量***的测量范围内布置3根长度标尺,每根长度标尺的两端各安装有一个靶标座,可适配多仪器组合测量***中的不同测量仪器的靶标,且保持测量点不随靶标类型改变而变动,长度标尺的方向包括水平、垂直、倾斜等不同方向。
步骤三、多仪器组合测量***对3根长度标尺进行测量,保证每根长度标尺的测量点均仅被一台测量仪器测量,且同一根长度标尺两端的测量点不被同一台测量仪器测量,即3根长度标尺上的测量点均为非公共测量点。根据多仪器组合测量***的外参数完成数据转换,获得非公共测量点的坐标。重复三遍。将第j遍中第i根长度标尺两端的公共测量点坐标分别记作(xi1j,yi1j,zi1j)和(xi2j,yi2j,zi2j),i=1,2,3,j=1,2,3。
步骤四、利用坐标(x11j,y11j,z11j)、(x12j,y12j,z12j),计算第1根长度标尺的长度lm1j,如式(25)所示,并计算其与该长度标尺的参考值lr1j的差值Δl1j,如式(26)所示。取Δl11、Δl12、Δl13中最大值作为第1根长度标尺处的空间长度示值误差Δl1,如式(27)所示。
步骤五、重复步骤四,依次计算剩余2根长度标尺处的空间长度示值误差Δli,i=2,3。取3根长度标尺处的空间长度示值误差的最大值作为本次校准的非公共测量点的空间长度示值误差Δl,如式(28)所示。
附图1-5中的附图标记如下:
1第一测量仪器,
2第二测量仪器,
3第三测量仪器,
4第四测量仪器,
5第N-1测量仪器,
6第N测量仪器,
7被测工件,
8公共测量点,
9激光干涉仪,
10直线导轨,
11第一长度标尺,
12第二长度标尺,
13第三长度标尺,
14非公共测量点,
15多齿分度台,
16移动工作台。
Claims (1)
1.一种多仪器组合测量***的校准方法,其特征在于:
采用多齿分度台作为标准器校准多仪器组合测量***的坐标系旋转误差;
采用带有激光干涉仪的直线导轨作为标准器校准多仪器组合测量***的坐标系平移误差;
采用长度标尺作为标准器校准多仪器组合测量***的公共测量点和非公共测量点的空间距离长度误差;
多仪器组合测量***的坐标系旋转误差校准包括以下步骤:
步骤1.1、由N台测量仪器组成的多仪器组合测量***,将其中1台测量仪器固定在多齿分度台上,其余N-1台测量仪器固定在多齿分度台旁,N台测量仪器的整体布局与其现场应用时一致;
步骤1.2、多齿分度台刻度线归零,标定***中N台测量仪器的外参数,将多齿分度台上的测量仪器的坐标系旋转矩阵记作R0;
步骤1.3、多齿分度台旋转角度θ,记录此时多齿分度台示值θr1;
步骤1.4、标定***中N台测量仪器的外参数,重复三遍,将多齿分度台上的测量仪器的坐标系旋转矩阵记作R11、R12、R13;
步骤1.5、分别计算坐标系旋转矩阵R11、R12、R13与R0之间的变化量ΔR11、ΔR12、ΔR13,如式(1)所示;结合罗德里格斯公式计算对应的旋转角θm11、θm12、θm13,如式(2)所示;并计算与θr1的差值Δθ11、Δθ12、Δθ13,如式(3)所示;取最大值记作该角度处的坐标系旋转误差Δθ1,如式(4)所示;
其中,r1js表示矩阵ΔR1j中第s个元素,j=1,2,3;
θm1j=arcsin(u2+v2+w2) (2);
Δθ1j=|θm1j-θr1| (3);
其中,j=1,2,3;
Δθ1=max{Δθ11,Δθ12,Δθ13} (4);
步骤1.6、按照步进θ,依次旋转多齿分度台至示值θr2、θr3、…、θrM处,重复步骤1.4、步骤1.5,测得各处的坐标系旋转误差Δθ2、Δθ3、…、ΔθM,并取最大值作为本次校准的坐标系旋转误差Δθ,如式(5)所示;
Δθ=max{Δθ1,Δθ2,...,ΔθM} (5);
多仪器组合测量***的坐标系平移误差校准包括以下步骤:
步骤2.1、由N台测量仪器组成的多仪器组合测量***,将其中1台测量仪器固定在直线导轨的移动工作台上,其余N-1台测量仪器固定在直线导轨旁,N台测量仪器的整体布局与其现场应用时一致;
步骤2.2、激光干涉仪复位,标定***中N台测量仪器的外参数,将移动工作台上的测量仪器的坐标系平移向量记作T0;
步骤2.3、移动工作台沿直线导轨移动位移d,记录此时激光干涉仪器的示值dr1;
步骤2.4、标定***中N台测量仪器的外参数,重复三遍,将移动工作台上的测量仪器的坐标系平移向量记作T11、T12、T13;
步骤2.5、分别计算坐标系平移向量T11、T12、T13与T0之间的变化量ΔT11、ΔT12、ΔT13,如式(6)所示,结合向量模的计算公式计算对应的平移量dm11、dm12、dm13,如式(7)所示,并计算与dr1的差值Δd11、Δd12、Δd13,如式(8)所示,取最大值记作该位置处的坐标系平移误差Δd1,如式(9)所示;
其中,tx1j、ty1j、tz1j表示向量ΔT1j中3个元素,j=1,2,3;
其中,j=1,2,3;
Δd1j=|dm1j-dr1| (8);
其中,j=1,2,3;
Δd1=max{Δd11,Δd12,Δd13} (9);
步骤2.6、按照步进d,依次沿直线导轨将移动工作台移至示值dr2、dr3、…、drM处,重复步骤2.4、步骤2.5,测得各处的坐标系平移误差Δd2、Δd3、…、ΔdM,并取最大值作为本次校准的坐标系平移误差Δd,如式(10)所示;
Δd=max{Δd1,Δd2,...,ΔdM} (10);
多仪器组合测量***的公共测量点的空间长度示值误差校准包括以下步骤:
步骤3.1、由N台测量仪器组成的多仪器组合测量***,按照与现场应用时一致的整体布局放置,在多仪器组合测量***的测量范围内布置M根长度标尺,每根长度标尺的两端各安装有一个靶标座,可适配多仪器组合测量***中的不同测量仪器的靶标,且保持测量点不随靶标类型改变而变动,长度标尺的方向包括水平、垂直、倾斜方向;
步骤3.2、保证M根长度标尺两端的测量点均可被多仪器组合测量***中2台及以上的测量仪器测得,将长度标尺两端的测量点作为公共测量点,配合其他公共测量点,标定多仪器组合测量***的外参数并通过数据融合计算公共测量点的坐标;重复三遍,将第j遍中第i根长度标尺两端的公共测量点坐标分别记作(xi1j,yi1j,zi1j)和(xi2j,yi2j,zi2j),i=1,2,…,M,j=1,2,3;
步骤3.3、利用坐标(x11j,y11j,z11j)、(x12j,y12j,z12j),计算第1根长度标尺的长度lm1j,如式(11)所示,并计算其与该长度标尺的参考值lr1j的差值Δl1j,如式(12)所示;取Δl11、Δl12、Δl13中最大值作为第1根长度标尺处的空间长度示值误差Δl1,如式(13)所示;
其中,j=1,2,3;
Δl1j=|lm1j-lr1| (12);
其中,j=1,2,3;
Δl1=max{Δl11,Δl12,Δl13} (13);
步骤3.4、重复步骤3.3,依次计算剩余M-1根长度标尺处的空间长度示值误差Δli,i=2,3,…,M;取M根长度标尺处的空间长度示值误差的最大值作为本次校准的公共测量点的空间长度示值误差Δl,如式(14)所示;
Δl=max{Δl1,Δl2,...,ΔlM} (14);
多仪器组合测量***的非公共测量点的空间长度示值误差校准包括以下步骤:
步骤4.1、由N台测量仪器组成的多仪器组合测量***,按照与现场应用时一致的整体布局放置,标定多仪器组合测量***的外参数;
步骤4.2、在多仪器组合测量***的测量范围内布置M根长度标尺,每根长度标尺的两端各安装有一个靶标座,可适配多仪器组合测量***中的不同测量仪器的靶标,且保持测量点不随靶标类型改变而变动,长度标尺的方向包括水平、垂直、倾斜方向;
步骤4.3、多仪器组合测量***对M根长度标尺进行测量,保证每根长度标尺的测量点均仅被一台测量仪器测量,且同一根长度标尺两端的测量点不被同一台测量仪器测量,即M根长度标尺上的测量点均为非公共测量点;根据多仪器组合测量***的外参数完成数据转换,获得非公共测量点的坐标;重复三遍;将第j遍中第i根长度标尺两端的公共测量点坐标分别记作(xi1j,yi1j,zi1j)和(xi2j,yi2j,zi2j),i=1,2,…,M,j=1,2,3;
步骤4.4、利用坐标(xi1j,yi1j,zi1j)、(xi2j,yi2j,zi2j),计算第1根长度标尺的长度lm1j,如式(11)所示,并计算其与该长度标尺的参考值lr1j的差值Δl1j,如式(12)所示;取Δl11、Δl12、Δl13中最大值作为第1根长度标尺处的空间长度示值误差Δl1,如式(13)所示;
步骤4.5、重复步骤4.4,依次计算剩余M-1根长度标尺处的空间长度示值误差Δli,i=2,3,…,M;取M根长度标尺处的空间长度示值误差的最大值作为本次校准的非公共测量点的空间长度示值误差Δl,如式(14)所示;
N和M为正整数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110576171.9A CN113310445B (zh) | 2021-05-26 | 2021-05-26 | 一种多仪器组合测量***的校准方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110576171.9A CN113310445B (zh) | 2021-05-26 | 2021-05-26 | 一种多仪器组合测量***的校准方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113310445A CN113310445A (zh) | 2021-08-27 |
CN113310445B true CN113310445B (zh) | 2023-06-06 |
Family
ID=77374962
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110576171.9A Active CN113310445B (zh) | 2021-05-26 | 2021-05-26 | 一种多仪器组合测量***的校准方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113310445B (zh) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102589430A (zh) * | 2012-02-01 | 2012-07-18 | 哈尔滨工业大学 | 多仪器坐标统一化装置的校准方法 |
CN106524905A (zh) * | 2016-10-11 | 2017-03-22 | 北京工业大学 | 一种基于激光追踪仪多站位测量的四轴机床标定方法 |
CN107265221A (zh) * | 2016-04-08 | 2017-10-20 | 奥的斯电梯公司 | 用于多个3d传感器校准的方法和*** |
CN111238372A (zh) * | 2020-02-24 | 2020-06-05 | 重庆市计量质量检测研究院 | 双复合式坐标测量***的联合位置误差的同步检测方法 |
CN111678471A (zh) * | 2020-06-09 | 2020-09-18 | 无锡身为度信息技术有限公司 | 圆柱坐标测量机转台误差辨识与补偿方法 |
CN112509067A (zh) * | 2021-02-02 | 2021-03-16 | 中智行科技有限公司 | 多传感器联合标定方法、装置、电子设备和存储介质 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7840372B2 (en) * | 2007-07-06 | 2010-11-23 | The University Of British Columbia | Self-calibration method and apparatus for on-axis rotary encoders |
-
2021
- 2021-05-26 CN CN202110576171.9A patent/CN113310445B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102589430A (zh) * | 2012-02-01 | 2012-07-18 | 哈尔滨工业大学 | 多仪器坐标统一化装置的校准方法 |
CN107265221A (zh) * | 2016-04-08 | 2017-10-20 | 奥的斯电梯公司 | 用于多个3d传感器校准的方法和*** |
CN106524905A (zh) * | 2016-10-11 | 2017-03-22 | 北京工业大学 | 一种基于激光追踪仪多站位测量的四轴机床标定方法 |
CN111238372A (zh) * | 2020-02-24 | 2020-06-05 | 重庆市计量质量检测研究院 | 双复合式坐标测量***的联合位置误差的同步检测方法 |
CN111678471A (zh) * | 2020-06-09 | 2020-09-18 | 无锡身为度信息技术有限公司 | 圆柱坐标测量机转台误差辨识与补偿方法 |
CN112509067A (zh) * | 2021-02-02 | 2021-03-16 | 中智行科技有限公司 | 多传感器联合标定方法、装置、电子设备和存储介质 |
Non-Patent Citations (1)
Title |
---|
何俊 ; 张福民 ; 张画迪 ; 曲兴华 ; .基于球心拟合的多边激光跟踪***自标定方法.《红外与激光工程.2020,49(08),1-7. * |
Also Published As
Publication number | Publication date |
---|---|
CN113310445A (zh) | 2021-08-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP2764103B2 (ja) | アナログ測定プローブの使用方法および位置決め装置 | |
CN107042528A (zh) | 一种工业机器人的运动学标定***及方法 | |
JPH02284216A (ja) | アナログプローブ較正方法と装置 | |
CN110487182B (zh) | 一种基于多维动态标准器的坐标转换方法 | |
CN111595238B (zh) | 一种基于多站法的激光跟踪仪精度现场评价*** | |
CN109163658B (zh) | 一种可提供位置和角度基准的光学基准件的标定方法 | |
CN102721393A (zh) | 一种精密旋转工作台测量***误差在位自标定方法 | |
US5649368A (en) | Method for calibrating a coordinate measuring apparatus having two pivot axes | |
CN111854587B (zh) | 一种导轨五自由度运动误差在线测量装置及方法 | |
CN110530296A (zh) | 一种线激光安装误差角确定方法 | |
CN110220454A (zh) | 一种三坐标定位机构的位姿标定方法 | |
CN112229321B (zh) | 基于lasso算法求解三坐标测量机21项几何误差的方法 | |
JP2561861B2 (ja) | 組合わされたスケールと干渉計 | |
JP2020046301A (ja) | 工作機械の計測誤差評価方法及びプログラム | |
CN113514017B (zh) | 一种并联驱动机构动平台位姿测量方法 | |
CN112945109B (zh) | 基于水平位移台的激光位移计阵列***参数标定方法 | |
CN113310445B (zh) | 一种多仪器组合测量***的校准方法 | |
CN112697074B (zh) | 动态待测物角度测量仪及测量方法 | |
JP2013171043A (ja) | 加工品に対して加工工具を位置決めするための構造および方法 | |
CN113513986B (zh) | 几何公差测量设备及其测量方法 | |
CN114136357A (zh) | 一种适用于面结构光传感器的测试方法及测试*** | |
Klobucar et al. | Automatic high resolution measurement set-up for calibrating precise line scales | |
CN112050731A (zh) | 基于大型虚拟计量器的测量数据配准精度评价方法 | |
CN117091514B (zh) | 双层浮动读数头的光栅位移测量装置、方法、介质及设备 | |
CN111505540B (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 |