CN106541419A - 一种机器人轨迹误差的测量方法 - Google Patents
一种机器人轨迹误差的测量方法 Download PDFInfo
- Publication number
- CN106541419A CN106541419A CN201610891618.0A CN201610891618A CN106541419A CN 106541419 A CN106541419 A CN 106541419A CN 201610891618 A CN201610891618 A CN 201610891618A CN 106541419 A CN106541419 A CN 106541419A
- Authority
- CN
- China
- Prior art keywords
- robot
- rotating
- point
- camera
- motion
- 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
- 238000000034 method Methods 0.000 title claims abstract description 48
- 230000033001 locomotion Effects 0.000 claims abstract description 106
- 238000004364 calculation method Methods 0.000 claims description 5
- 239000003054 catalyst Substances 0.000 claims description 2
- 238000005259 measurement Methods 0.000 abstract description 15
- 238000003384 imaging method Methods 0.000 abstract description 4
- 238000010586 diagram Methods 0.000 description 6
- 238000004422 calculation algorithm Methods 0.000 description 5
- 238000004088 simulation Methods 0.000 description 4
- 238000001514 detection method Methods 0.000 description 2
- 238000000691 measurement method Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000012636 effector Substances 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000003550 marker Substances 0.000 description 1
- 230000009347 mechanical transmission Effects 0.000 description 1
- 238000010408 sweeping Methods 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B25—HAND TOOLS; PORTABLE POWER-DRIVEN TOOLS; MANIPULATORS
- B25J—MANIPULATORS; CHAMBERS PROVIDED WITH MANIPULATION DEVICES
- B25J19/00—Accessories fitted to manipulators, e.g. for monitoring, for viewing; Safety devices combined with or specially adapted for use in connection with manipulators
- B25J19/02—Sensing devices
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B25—HAND TOOLS; PORTABLE POWER-DRIVEN TOOLS; MANIPULATORS
- B25J—MANIPULATORS; CHAMBERS PROVIDED WITH MANIPULATION DEVICES
- B25J9/00—Programme-controlled manipulators
- B25J9/16—Programme controls
- B25J9/1628—Programme controls characterised by the control loop
- B25J9/1653—Programme controls characterised by the control loop parameters identification, estimation, stiffness, accuracy, error analysis
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B25—HAND TOOLS; PORTABLE POWER-DRIVEN TOOLS; MANIPULATORS
- B25J—MANIPULATORS; CHAMBERS PROVIDED WITH MANIPULATION DEVICES
- B25J9/00—Programme-controlled manipulators
- B25J9/16—Programme controls
- B25J9/1674—Programme controls characterised by safety, monitoring, diagnostic
Landscapes
- Engineering & Computer Science (AREA)
- Robotics (AREA)
- Mechanical Engineering (AREA)
- Manipulator (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明涉及一种机器人轨迹误差的循迹测量方法,本发明可以测量机器人全程或局部最大运动误差:若需要测量机器人某一段时间内或全程最大运动误差,在机器人运动过程中,通过视觉成像系统全程实时采集机器人末端实际运动轨迹与轨迹发生器产生的理论轨迹,通过比较两种轨迹,可得到机器人某一段时间内或全程最大运动误差。本发明可以测量机器人运动过程中关键离散点运动误差:在机器人运动过程中,控制机器人在某一关键点处停止运动,通过视觉成像系统实时采集此时机器人末端实际空间点坐标与轨迹发生器产生的理论空间点坐标,计算两点之间距离,可得到关键离散点处的机器人运动误差。本发明采用循迹测量的方法,在机器人运动过程中,对机器人末端的运动轨迹进行实时误差测量。
Description
技术领域
本发明涉及机器人运动轨迹误差测量方法,具体涉及一种轨迹发生器和视觉测量系统相结合的机器人轨迹误差的测量方法。
背景技术
在先技术(机械科学与技术,2011年2期,20卷,第252页,钱瑞明《挠性机器人动态误差集成激光测量与补偿研究》)提出一种基于3个激光发生器和3个激光位置检测器(PSD)的挠性构件动态误差测量方法,该法可以测量除杆长方向外的5个变形分量,且光路和测量模型简单;建立了PSD上光点位置与构件各误差分量及机器人末端执行器动态误差之间的关系,并给出其补偿控制方法。但其控制过程复杂,计算繁琐,容易出现因控制误差而导致测量结果不准确的问题。
在先技术(机械传动,2013年5期,6卷,第50页,王良文《用于物体捕捉的四足步行机器人的运动误差模型》)提出在四足机器人的机体上,安装图像捕捉系统,用于引导机器人完成对目标物的抓取。并在对抓取状态逆运动学分析的基础上,建立精确的图像捕捉系统误差与机器人的工作臂参数误差之间的关系,给出了详细的计算公式,得到机器人运动过程中的运动误差并进行补偿。此机器人运动误差的获取是通过采集外界环境的改变的图像来逆向计算出机器人工作臂运动误差,测量精度受后期图像处理精度限制。
在先技术(浙江大学,专利号:201010552545.5)提到一种基于扫频激光干涉的圆轨迹运动误差快速测量系统:通过总分光镜将来自扫频激光器的光信号分成X向光信号和Y向光信号,并打到安装在机床导轨的靶镜上。最后,通过X向检测机构和Y向检测机获取机床导轨的圆轨迹运动误差。此发明装置只能检测平面上机器人轨迹运动误差,应用领域较窄。
发明内容
本发明的目的在于提供一种机器人轨迹误差的测量方法,本发明可以实现机器人全程或局部最大运动误差测量:在机器人运动过程中,通过视觉成像系统全程实时采集机器人末端实际运动轨迹,同时轨迹发生器产生机器人运动的理论轨迹,通过比较两种轨迹,可得到机器人某一段时间内或全程最大运动误差。
本发明提出的机器人轨迹误差的测量方法,所述方法通过用于轨迹发生器的旋转双棱镜系统和用于轨迹图像采集的双目视觉测量系统实现,所述方法用于对机器人全程或局部最大运动误差进行测量,其中:
所述旋转双棱镜系统包括第一旋转双棱镜1和第二旋转双棱镜2,第一旋转双棱镜1和第二旋转双棱镜2同轴布置;所述双目视觉测量系统包括第一相机3、第二相机4、第一支杆5、第二支杆6和底板7,第一相机3与第一支杆5一端相连,第一支杆5另一端固定于底板7上;第二相机4与第二支杆6一端相连,第二支杆6另一端固定于底板7上。第一相机3和第二相机4同时对机器人末端9进行拍摄,并设置同样的拍摄间隔;对拍摄的照片进行特征匹配,匹配的内容为机器人末端标记点10和激光点11;具体步骤如下:
(1)在关节机器人末端粘贴一个标记点10。
(2)根据关节机器人各关节的运动角度以及角速度,计算出标记点10随时间变化的理论运动位置曲线T1(X,Y,Z);
(3)基于旋转双棱镜逆向方法,根据步骤(2)所述的理论运动位置曲线T1(X,Y,Z)可以计算出旋转双棱镜1和旋转双棱2的旋转角度(θ1(t),θ2(t))和旋转速度(v1(t),v2(t));其中:
(4)控制机器人和旋转双棱镜系统同时动作,入射激光从第一旋转双棱镜1入射面中心垂直入射,并根据已知的旋转角度(θ1(t),θ2(t))和旋转速度(v1(t),v2(t)),控制第一旋转双棱镜1和第二旋转双棱镜2旋转;使得出射激光打在机器人末端9上标记点10的理论运动位置曲线T1(X,Y,Z)上,作为机器人的理论轨迹发生器;
(5)第一相机3和第二相机4对机器人末端9进行拍照;
(6)对同一时刻第一相机3和第二相机4拍摄的照片,进行机器人末端标记点的图像匹配,得到机器人末端标记点在第一相机3内的图像坐标系坐标(x1,y1),在第二相机4内的图像坐标系坐标(x2,y2);
(7)对同一时刻第一相机3和第二相机4拍摄的照片,进行机器人末端上激光点的图像匹配,得到激光点在第一相机3内的图像坐标系坐标(x3,y3),在第二相机4内的图像坐标系坐标(x4,y4);
(8)基于双目视觉测量方法,通过步骤(6)得到的标记点10在第一相机3内的图像坐标系坐标(x1,y1),与在第二相机4内的图像坐标系坐标(x2,y2),可计算出机器人运动时末端上标记点10的三维坐标(X1,Y1,Z2);由步骤(7)得到的激光点11在第一相机3内的图像坐标系坐标(x3,y3),与在第二相机4内的图像坐标系坐标(x4,y4),可计算出机器人运动时末端上激光点11的三维坐标(X2,Y2,Z2);
(9)根据得到的一系列激光点11的三维坐标,进行机器人实际运动轨迹拟合,得到机器人实际运动轨迹曲线为T2(X,Y,Z);
(10)通过比较实际运动轨迹和理论运动轨迹两条曲线间距离,可得到机器人全程或局部最大运动误差。
本发明中,步骤(3)中所述旋转双棱镜逆向方法,具体步骤如下:
(1)第一旋转双棱镜1和第二旋转双棱镜2的楔角为α,入射光从第一旋转双棱镜1垂直入射,并打在第一旋转双棱镜1入射面中心。已知目标点P(Xp,Yp,Zp),并设置出射点位置精度为Δ0;
(2)假设出射光在第二旋转双棱镜2上的出射点为第二旋转双棱镜2出射面中心N0(0,0),以为出射光向量,得到此时出射光俯仰角ρ1和方位角
ρ1=arccos(zP),
(3)令第一旋转双棱镜1旋转角度θ1=0°且静止不动,仅第二旋转双棱镜2旋转。基于旋转双棱镜矢量折射定理可得:ρ=arccos(cosδ1cosδ2-sinδ1sinδ2cos△θr);
式中:△θr=θ1-θ2为两棱镜夹角,δ1为棱镜1的偏向角,即光束通过棱镜1后出射光偏离入射光的角度。可求得δ1=arcsin(n·sinα)-α。δ2为棱镜2的偏向角,即光束通过棱镜2后出射光偏离入射光的角度。可求得其中,为棱镜2的等效折射率,其值为γr=arctan(tanδ1·cos△θr),βr=arccos(sinδ1sin△θr)。将步骤(2)得到的出射光俯仰角ρ1带入公式,求出双棱镜夹角为Δθr,即第二旋转双棱镜2旋转角度为θ2=-Δθr。此时,基于矢量折射定理,根据第一旋转双棱镜1和第二旋转双棱镜2的旋转角度(0°,-Δθr),计算此时出射光的方位角为
(4)根据步骤(2)得到的方位角和步骤三得到方位角将棱镜1和棱镜2的旋转角度同时增大后,即可求得第一旋转双棱镜1和第二旋转双棱镜2旋转角度(θ1,θ2)为
(5)根据步骤(4)得到的两棱镜转角(θ1,θ2),基于旋转双棱镜矢量折射定理,求解实际打在机器人末端点P1(X1 p,Y1 p,Z1 p)和棱镜2出射点位置N1(x1 n,y1 n);
(6)计算偏差判断是否满足精度要求,即是否满足Δ﹤Δ0,Δ0为给定的目标点精度;
(7)若Δ≥δ,回到步骤(2)中,计算以为出射光向量时的俯仰角ρ1和方位角重复步骤(2)、(3)、(4)、(5)和(6),若Δ﹤Δ0,则结束,得出旋转双棱镜转角解(θ1,θ2)。
本发明中,步骤(10)中所述全程或局部最大运动误差计算方法,具体步骤如下:
(1)根据得到的理论运动轨迹曲线T1(X,Y,Z)和实际运动轨迹曲线T2(X,Y,Z)。将两条运动轨迹曲线放置在同一空间坐标下。
(2)将理论运动轨迹曲线T1(X,Y,Z)等分成n1个点,对于其中的任一点Pn(Xn,Yn,Zn),可以计算出此点处的切线方程,方程斜率为kn;
(3)将实际运动轨迹曲线T2(X,Y,Z)等分成n2个点,对于其中的任一点Pm(Xm,Ym,Zm),可以计算出此点处的切线方程,方程斜率为km;
(4)进行对应点的匹配,设置精度为δ0,对于实际轨迹曲线上的任一点Pm(Xm,Ym,Zm),切线方程的斜率为km,在理论轨迹曲线上找一点Pn(Xn,Yn,Zn),同时满足|km-kn|﹤δ0、|km-1-kn-1|﹤δ0和|km+1-kn+1|﹤δ0。则认为实际曲线上的Pm点对应理论曲线上的Pn点。则机器人运动到Pm处的运动误差为
(5)重复步骤(4)计算出机器人实际运动曲线上任一点的运动误差,可得到机器人全程或局部最大运动误差。
本发明的有益效果在于:
1.本方法属于非接触测量机器人运动轨迹误差方法,不需要直接接触机器人。保证了在不破坏机器人的情况下,能精确得到机器人运动轨迹误差。
2.易于得到精准的机器人运动轨迹误差,本发明采用双目视觉测量机器人运动轨迹误差的方法,通过两个相机对机器人末端标记点和激光扫描点拍照,利用双目视觉测量原理,可以精确的得到机器人末端标记点和激光扫描点三维坐标值,通过计算标记点和扫描点的距离得到机器人运动轨迹误差,精度最高可达到微米级别。
3.易于实时得到机器人运动轨迹误差,本发明采用循迹测量的方法,在机器人运动过程中,对机器人末端的运动轨迹进行实时误差测量。
4.控制方便。该装置采用独立控制方式,只需要控制旋转双棱镜旋转达到指定的角度(θ1,θ2),保证通过旋转双棱镜的激光持续打在机器人末端的理论轨迹上;对于第一相机3和第二相机4,只需使用PC机控制第一相机3和第二相机4同时拍摄照片,并设置相同的拍摄时间间隔,控制过程简单,易于实现。
附图说明
图1是本发明一种机器人的轨迹误差的循迹测量方法的总示意图。由三个部分组成,分别为关节机器人8、双目视觉测量系统和旋转双棱镜系统。
图2是机器人末端9标记点和激光点11示意图。
图3为轨迹发生器产生机器人末端理论轨迹的流程图。
图4为旋转双棱镜逆向算法示意图。
图5为得到机器人末端9实际运动轨迹曲线和理论运动轨迹曲线流程图。
图6为全程或局部最大运动误差计算方法流程图。
图7为机器人局部或全程运动轨迹最大误差测量流程图。
图8机器人理论运动轨迹仿真图。
图9旋转双棱镜1和旋转双棱2的旋转角度(θ1(t),θ2(t))仿真图。
图10旋转双棱镜1和旋转双棱2的旋转速度(v1(t),v2(t))仿真图。
图11机器人理论运动轨迹和实际运动轨迹仿真图。
图中标号:1为第一旋转双棱镜,2为第二旋转双棱镜,3为第一相机,4为第二相机,5为第一支杆,6为第二支杆,7为底板,8为机器人,9为机器人末端,10为标记点,11为激光点。
具体实施方式
下面通过实施例结合附图进一步说明本发明。
实施例1:
本发明提供一种机器人轨迹误差的测量方法,可以实现机器人全程或局部最大运动误差测量:在机器人运动过程中,通过视觉成像系统全程实时采集机器人末端实际运动轨迹,同时轨迹发生器产生机器人运动的理论轨迹,通过比较两种轨迹,可得到机器人某一段时间内或全程最大运动误差。
2.本发明的目的通过下述几个部分来实现:包括用于轨迹发生器的旋转双棱镜系统和轨迹图像采集的双目视觉测量系统。根据机器人理论运动轨迹,旋转双棱镜产生高精度的光束扫描轨迹照射在机器人末端上,同时保证该光束扫描轨迹的跟踪时间和跟踪坐标与机器人运动实际的同步性,从而高精度模拟机器人末端运动理论轨迹。后续文字我们称该扫描轨迹为机器人理论轨迹。双目视觉测量系统用于采集机器人末端上标记点的实际运动轨迹和激光点的理论运动轨迹,最后根据得到的理论运动轨迹点坐标和实际运动轨迹点坐标计算出理论运动论轨迹与实际运动轨迹的误差,得到机器人末端运动时的运动轨迹误差。
3.结合图1、图2和图3说明旋转双棱镜产生机器人末端理论轨迹的步骤:S1:如图2所示,在关节机器人末端9中心粘贴一标记点10。S2:根据关节机器人各关节的运动角度以及角速度等信息可以计算出标记点10随时间变化的理论运动位置方程T1(X,Y,Z,t)。S3:基于旋转双棱镜逆向算法,根据标记点空间位置T1(X,Y,Z,t)可以计算出旋转双棱镜1和旋转双棱2的旋转角度(θ1(t),θ2(t))。S4:根据第一旋转双棱镜1和第二旋转双棱2的旋转角度(θ1(t),θ2(t)),可以计算得出旋转速度S4:激光从第一旋转双棱镜1入射面中心垂直入射,第二旋转双棱镜2出射,并根据已知的(θ1(t),θ2(t))和(v1(t),v2(t))控制第一旋转双棱镜1和第二旋转双棱镜2角旋转。使得出射激光打在机器人末端9上标记点10的理论位置,作为机器人的理论轨迹发生器。
4.结合图4和图5说明旋转双棱镜逆向算法的过程,图4中第一旋转双棱镜1和第二旋转双棱镜2的楔角为α,入射光从第一旋转双棱镜1垂直入射,并打在第一旋转双棱镜1入射面中心。已知目标点P(Xp,Yp,Zp),并设置出射点位置精度为Δ0。根据逆向算法算可求出第一旋转双棱镜转角(θ1,θ2),其步骤如下:
S1:假设出射光在第二旋转双棱镜2上的出射点为第二旋转双棱镜2出射面中心N0(0,0),以为出射光向量,求出出射光向量的俯仰角ρ和方位角φ。
S2:令第一旋转双棱镜1不动,仅第二旋转双棱镜2旋转,根据出射光俯仰角ρ求出第二双棱镜夹角Δθr,此时的出射光方位角为φ0。根据出射光方位角φ,第一旋转双棱镜1和第二旋转双棱镜2同时旋转φ-φ0,即得到第一旋转双棱镜1和第二旋转双棱镜2转角(θ1,θ2)。或者令第二旋转双棱镜2不动,仅第一旋转双棱镜1旋转,根据出射光俯仰角ρ求出第一旋转双棱镜1和第二旋转双棱镜2夹角Δθr,此时的出射光方位角为φ0。根据出射光方位角φ,第一旋转双棱镜1和第二旋转双棱镜2同时旋转φ-φ0,即得到第一旋转双棱镜1和第二旋转双棱镜2转角(θ1,θ2)。
S3:根据第一旋转双棱镜1和第二旋转双棱镜2转角(θ1,θ2),求解实际打在机器人末端点P1(X1 p,Y1 p,Z1 p)和第二旋转双棱镜2出射点位置N1(x1 n,y1 n)。
S4:计算偏差判断是否满足精度要求,即是否满足Δ﹤Δ0,Δ0为给定的目标点精度。
S5:若Δ≥δ,将作为出射光向量,重复步骤(2)、(3)、(4);若Δ﹤Δ0,则结束,得出第一旋转双棱镜1和第二旋转双棱镜2转角解(θ1,θ2)。
5.机器人全程或局部最大运动误差计算原理如下:在得到机器人实际运动轨迹和理论运动轨迹的基础上,通过比较两条曲线间距离,可得到机器人全程或局部最大运动误差。
根据下列步骤完成全程或局部机器人运动轨迹最大误差的测量:
S1:如附图2所示,在关节机器人末端粘贴一个标记点10。
S2:根据关节机器人各关节的运动角度以及角速度等信息可以计算出标记点随时间变化的理论运动位置方程T1(X,Y,Z,t)为(单位cm),参数方程为(单位cm),其运动曲线如附图8所示,为一空间椭圆曲线。
S3:基于旋转双棱镜逆向算法,根据S2计算所得标记点空间位置T1(X,Y,Z,t)可以计算出旋转双棱镜1和旋转双棱2的旋转角度(θ1(t),θ2(t))和旋转速度(v1(t),v2(t)),如附图9和附图10所示。
S4:同时控制机器人和轨迹发生器同时动作。入射激光从旋转双棱镜1入射面中心垂直入射,并根据已知的(θ1(t),θ2(t))和(v1(t),v2(t))控制旋转双棱镜1和旋转双棱镜2旋转。使得出射激光打在机器人末端9上标记点10的理论位置,作为机器人的理论轨迹发生器。
S5:如图1所示,第一相机3和第二相机4对机器人末端9进行拍照。
S6:对同一时刻第一相机3和第二相机4拍摄的照片,进行机器人末端标记点的图像匹配,得到机器人末端标记点在第一相机3内的图像坐标系坐标(x1,y1)和第二相机4内的图像坐标系坐标(x2,y2)。
S7:针对同一时刻第一相机3和第二相机4拍摄的照片,进行机器人末端上激光点的图像匹配,得到激光点在第一相机3内的图像坐标系坐标(x3,y3)和第二相机4内的图像坐标系坐标(x4,y4)。
S8:基于双目视觉测量原理,由(x1,y1)和(x2,y2)可计算出机器人运动时末端上标记点10的实际三维坐标(X1,Y1,Z2)。由(x3,y3)和(x4,y4)可计算出机器人运动时末端上标记点的理论三维坐标(X2,Y2,Z2)即激光点三维坐标。
S9:根据得到的所有标记点10和激光点11的三维坐标,进行机器人实际运动轨迹和理论运动轨迹拟合,如附图11所示,实线为得到的机器人理论运动轨迹,虚线为得到的机器人实际运动轨迹。
S10:根据全程或局部最大运动误差计算方法,如附图11所示,在机器人实际运动轨迹上的取一点Pm(42.43,42.43,28.28)(单位cm),可以在机器人理论运动轨迹上找到对应的点Pn(37.98,46.45,30.97)(单位cm),得到A点的机器人运动误差为△=6.57cm。S11:重复步骤S10,计算出机器人实际运动曲线上任一点的运动误差,可得到机器人全程或局部最大运动误差为8.98cm。
Claims (3)
1.一种机器人轨迹误差的测量方法,其特征在于所述方法通过用于轨迹发生器的旋转双棱镜系统和用于轨迹图像采集的双目视觉测量系统实现,所述方法用于对机器人全程或局部最大运动误差进行测量,其中:
所述旋转双棱镜系统包括第一旋转双棱镜(1)和第二旋转双棱镜(2),第一旋转双棱镜(1)和第二旋转双棱镜(2)同轴布置;所述双目视觉测量系统包括第一相机(3)、第二相机(4)、第一支杆(5)、第二支杆(6)和底板(7),第一相机(3)与第一支杆(5)一端相连,第一支杆(5)另一端固定于底板(7)上;第二相机(4)与第二支杆(6)一端相连,第二支杆(6)另一端固定于底板(7)上;第一相机(3)和第二相机(4)同时对机器人末端(9)进行拍摄,并设置同样的拍摄间隔;对拍摄的照片进行特征匹配,匹配的内容为机器人末端标记点(10)和激光点(11);具体步骤如下:
(1)在关节机器人末端粘贴一个标记点(10);
(2)根据关节机器人各关节的运动角度以及角速度,计算出标记点(10)随时间变化的理论运动位置曲线T1(X,Y,Z);
(3)基于旋转双棱镜逆向方法,根据步骤(2)所述的理论运动位置曲线T1(X,Y,Z)可以计算出旋转双棱镜(1)和旋转双棱(2)的旋转角度(θ1(t),θ(2)(t))和旋转速度(v1(t),v(2)(t));其中:
(4)控制机器人和旋转双棱镜系统同时动作,入射激光从第一旋转双棱镜(1)入射面中心垂直入射,并根据已知的旋转角度(θ1(t),θ(2)(t))和旋转速度(v1(t),v(2)(t)),控制第一旋转双棱镜(1)和第二旋转双棱镜(2)旋转;使得出射激光打在机器人末端(9)上标记点(10)的理论运动位置曲线T1(X,Y,Z)上,作为机器人的理论轨迹发生器;
(5)第一相机(3)和第二相机(4)对机器人末端(9)进行拍照;
(6)对同一时刻第一相机(3)和第二相机(4)拍摄的照片,进行机器人末端标记点的图像匹配,得到机器人末端标记点在第一相机(3)内的图像坐标系坐标(x1,y1),在第二相机(4)内的图像坐标系坐标(x(2),y(2));
(7)对同一时刻第一相机(3)和第二相机(4)拍摄的照片,进行机器人末端上激光点的图像匹配,得到激光点在第一相机(3)内的图像坐标系坐标(x(3),y(3)),在第二相机(4)内的图像坐标系坐标(x(4),y(4));
(8)基于双目视觉测量方法,通过步骤(6)得到的标记点(10)在第一相机(3)内的图像坐标系坐标(x1,y1),与在第二相机(4)内的图像坐标系坐标(x(2),y(2)),可计算出机器人运动时末端上标记点(10)的三维坐标(X1,Y1,Z(2));由步骤(7)得到的激光点(11)在第一相机(3)内的图像坐标系坐标(x(3),y(3)),与在第二相机(4)内的图像坐标系坐标(x(4),y(4)),可计算出机器人运动时末端上激光点(11)的三维坐标(X(2),Y(2),Z(2));
(9)根据得到的一系列激光点(11)的三维坐标,进行机器人实际运动轨迹拟合,得到机器人实际运动轨迹曲线为T2(X,Y,Z);
(10)通过比较实际运动轨迹和理论运动轨迹两条曲线间距离,可得到机器人全程或局部最大运动误差。
2.根据权利要求1所述的方法,其特征在于步骤(3)中所述旋转双棱镜逆向方法,具体步骤如下:
(1)第一旋转双棱镜(1)和第二旋转双棱镜(2)的楔角为α,入射光从第一旋转双棱镜(1)垂直入射,并打在第一旋转双棱镜(1)入射面中心;已知目标点P(Xp,Yp,Zp),并设置出射点位置精度为Δ0;
(2)假设出射光在第二旋转双棱镜(2)上的出射点为第二旋转双棱镜(2)出射面中心N0(0,0),以为出射光向量,得到此时出射光俯仰角ρ1和方位角
ρ1=arccos(zP),
(3)令第一旋转双棱镜(1)旋转角度θ1=0°且静止不动,仅第二旋转双棱镜(2)旋转;基于旋转双棱镜矢量折射定理可得:ρ=arccos(cosδ1cosδ2-sinδ1sinδ2cos△θr);
式中:△θr=θ1-θ2为两棱镜夹角,δ1为棱镜1的偏向角,即光束通过棱镜(1)后出射光偏离入射光的角度;可求得δ1=arcsin(n·sinα)-α;δ2为棱镜(2)的偏向角,即光束通过棱镜(2)后出射光偏离入射光的角度;可求得其中,为棱镜(2)的等效折射率,其值为γr=arctan(tanδ1·cos△θr),βr=arccos(sinδ1sin△θr);将步骤(2)得到的出射光俯仰角ρ1带入公式,求出双棱镜夹角为Δθr,即第二旋转双棱镜(2)旋转角度为θ2=-Δθr;此时,基于矢量折射定理,根据第一旋转双棱镜(1)和第二旋转双棱镜(2)的旋转角度(0°,-Δθr),计算此时出射光的方位角为
(4)根据步骤(2)得到的方位角和步骤三得到方位角将棱镜(1)和棱镜(2)的旋转角度同时增大后,即可求得第一旋转双棱镜(1)和第二旋转双棱镜(2)旋转角度(θ1,θ2)为
(5)根据步骤(4)得到的两棱镜转角(θ1,θ2),基于旋转双棱镜矢量折射定理,求解实际打在机器人末端点P1(X1p,Y1p,Z1p)和棱镜(2)出射点位置N1(x1n,y1n);
(6)计算偏差判断是否满足精度要求,即是否满足Δ﹤Δ0,Δ0为给定的目标点精度;
(7)若Δ≥δ,回到步骤(2)中,计算以为出射光向量时的俯仰角ρ1和方位角重复步骤(2)、(3)、(4)、(5)和(6),若Δ﹤Δ0,则结束,得出旋转双棱镜转角解(θ1,θ2)。
3.根据权利要求1所述的方法,其特征在于步骤(10)中所述全程或局部最大运动误差计算方法,具体步骤如下:
(1)根据得到的理论运动轨迹曲线T1(X,Y,Z)和实际运动轨迹曲线T2(X,Y,Z)。将两条运动轨迹曲线放置在同一空间坐标下。
(2)将理论运动轨迹曲线T1(X,Y,Z)等分成n1个点,对于其中的任一点Pn(Xn,Yn,Zn),可以计算出此点处的切线方程,方程斜率为kn;
(3)将实际运动轨迹曲线T2(X,Y,Z)等分成n2个点,对于其中的任一点Pm(Xm,Ym,Zm),可以计算出此点处的切线方程,方程斜率为km;
(4)进行对应点的匹配,设置精度为δ0,对于实际轨迹曲线上的任一点Pm(Xm,Ym,Zm),切线方程的斜率为km,在理论轨迹曲线上找一点Pn(Xn,Yn,Zn),同时满足|km-kn|﹤δ0、|km-1-kn-1|﹤δ0和|km+1-kn+1|﹤δ0。则认为实际曲线上的Pm点对应理论曲线上的Pn点。则机器人运动到Pm处的运动误差为
(5)重复步骤(4)计算出机器人实际运动曲线上任一点的运动误差,可得到机器人全程或局部最大运动误差。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610891618.0A CN106541419B (zh) | 2016-10-13 | 2016-10-13 | 一种机器人轨迹误差的测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610891618.0A CN106541419B (zh) | 2016-10-13 | 2016-10-13 | 一种机器人轨迹误差的测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106541419A true CN106541419A (zh) | 2017-03-29 |
CN106541419B CN106541419B (zh) | 2019-01-25 |
Family
ID=58368663
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610891618.0A Expired - Fee Related CN106541419B (zh) | 2016-10-13 | 2016-10-13 | 一种机器人轨迹误差的测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106541419B (zh) |
Cited By (23)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107053176A (zh) * | 2017-04-09 | 2017-08-18 | 北京工业大学 | 一种六自由度机器人末端空间曲线轨迹的误差建模方法 |
CN107272015A (zh) * | 2017-07-05 | 2017-10-20 | 同济大学 | 高精度视觉导引激光跟踪方法 |
CN107391340A (zh) * | 2017-07-21 | 2017-11-24 | 郑州云海信息技术有限公司 | 一种整机柜服务器节点热拔插系统和控制方法 |
CN107803855A (zh) * | 2017-12-13 | 2018-03-16 | 大连四达高技术发展有限公司 | 一种机器人自动定位机构 |
CN108858218A (zh) * | 2018-06-28 | 2018-11-23 | 北京航星机器制造有限公司 | 一种适应自动物流的机械臂制孔装置及制孔方法 |
CN109186969A (zh) * | 2018-07-28 | 2019-01-11 | 西安交通大学 | 一种伺服进给运动动态性能视觉检测方法 |
CN109296024A (zh) * | 2018-11-30 | 2019-02-01 | 徐州市产品质量监督检验中心 | 一种无人挖掘机采装位姿精度检测方法 |
CN109514557A (zh) * | 2018-12-13 | 2019-03-26 | 北京炎凌嘉业机电设备有限公司 | 一种3d视觉机器人轨迹修正系统 |
CN109596125A (zh) * | 2017-09-30 | 2019-04-09 | 北京柏惠维康科技有限公司 | 一种机器人的空间坐标系转换关系的确定方法和装置 |
CN109884590A (zh) * | 2019-03-28 | 2019-06-14 | 河北工业大学 | 一种工业机器人轨迹精度检测装置及方法 |
CN111176304A (zh) * | 2020-03-18 | 2020-05-19 | 常州市贝叶斯智能科技有限公司 | 机器人运动底盘质检方法、装置、智能设备以及介质 |
CN111428626A (zh) * | 2020-03-23 | 2020-07-17 | 北京明略软件系统有限公司 | 一种移动物体的识别方法、装置、及存储介质 |
CN111805531A (zh) * | 2020-06-30 | 2020-10-23 | 同济大学 | 一种管道内窥机器人 |
CN111975780A (zh) * | 2020-08-25 | 2020-11-24 | 厦门众合天元科技有限公司 | 一种工业机器人运动轨迹设定装置及其使用方法 |
CN113977558A (zh) * | 2021-11-29 | 2022-01-28 | 湖南交通职业技术学院 | 一种可视化动态显示并联机器人末端轨迹的装置及方法 |
CN114454216A (zh) * | 2022-03-07 | 2022-05-10 | 云鲸智能(深圳)有限公司 | 机器人的精度检测方法、装置、机器人及存储介质 |
CN114454177A (zh) * | 2022-03-15 | 2022-05-10 | 浙江工业大学 | 一种基于双目立体视觉的机器人末端位置补偿方法 |
CN114563982A (zh) * | 2022-01-24 | 2022-05-31 | 中铁九桥工程有限公司 | 一种移动设备在圆管上运动轨迹的控制方法 |
CN115752321A (zh) * | 2022-11-09 | 2023-03-07 | 中山大学 | 医疗机器人运动轨迹测量比对方法及计算机可读存储介质 |
CN115847427A (zh) * | 2023-02-07 | 2023-03-28 | 成都秦川物联网科技股份有限公司 | 双重识别的协作机器人监测工业物联网系统及其控制方法 |
CN116197918A (zh) * | 2023-05-05 | 2023-06-02 | 北京华晟经世信息技术股份有限公司 | 基于动作记录分析的机械手控制系统 |
CN118123849A (zh) * | 2024-05-08 | 2024-06-04 | 鹏城实验室 | 机器人轨迹控制方法、装置、设备及存储介质 |
CN118238154A (zh) * | 2024-05-28 | 2024-06-25 | 深圳市欣茂鑫实业有限公司 | 一种用于自动上料的机械臂控制方法及系统 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6057967A (en) * | 1993-04-16 | 2000-05-02 | Nippon Telegraph And Telephone Corporation | Apparatus for extracting pattern features |
US20070257216A1 (en) * | 2003-08-27 | 2007-11-08 | Withers Michael J | Method and Apparatus for Investigating a Non-Planar Sample |
CN102135776A (zh) * | 2011-01-25 | 2011-07-27 | 解则晓 | 基于视觉定位的工业机器人控制系统及其控制方法 |
CN103231375A (zh) * | 2013-04-28 | 2013-08-07 | 苏州大学 | 基于距离误差模型的工业机器人标定方法 |
CN203509345U (zh) * | 2013-08-15 | 2014-04-02 | 中国电子科技集团公司第四十八研究所 | 一种焊接轨迹自动修正系统 |
CN104793334A (zh) * | 2015-04-02 | 2015-07-22 | 同济大学 | 一种级联粗精耦合光学扫描装置 |
CN104820400A (zh) * | 2015-04-18 | 2015-08-05 | 桂林鸿程机电设备有限公司 | 一种三维焊接机器人混合控制方法 |
CN105180834A (zh) * | 2015-05-28 | 2015-12-23 | 华中科技大学 | 一种叶片进排气边三维非接触式测量装置 |
DE102014014968A1 (de) * | 2014-10-14 | 2016-04-14 | Rwth Aachen | Optisches Messverfahren und Vorrichtung zur Bestimmung von Position und Orientierung von Werkstücken und/oder Maschinen |
CN105583825A (zh) * | 2016-03-14 | 2016-05-18 | 陈杨 | 一种工业机器人轨迹检测装置 |
-
2016
- 2016-10-13 CN CN201610891618.0A patent/CN106541419B/zh not_active Expired - Fee Related
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6057967A (en) * | 1993-04-16 | 2000-05-02 | Nippon Telegraph And Telephone Corporation | Apparatus for extracting pattern features |
US20070257216A1 (en) * | 2003-08-27 | 2007-11-08 | Withers Michael J | Method and Apparatus for Investigating a Non-Planar Sample |
CN102135776A (zh) * | 2011-01-25 | 2011-07-27 | 解则晓 | 基于视觉定位的工业机器人控制系统及其控制方法 |
CN103231375A (zh) * | 2013-04-28 | 2013-08-07 | 苏州大学 | 基于距离误差模型的工业机器人标定方法 |
CN203509345U (zh) * | 2013-08-15 | 2014-04-02 | 中国电子科技集团公司第四十八研究所 | 一种焊接轨迹自动修正系统 |
DE102014014968A1 (de) * | 2014-10-14 | 2016-04-14 | Rwth Aachen | Optisches Messverfahren und Vorrichtung zur Bestimmung von Position und Orientierung von Werkstücken und/oder Maschinen |
CN104793334A (zh) * | 2015-04-02 | 2015-07-22 | 同济大学 | 一种级联粗精耦合光学扫描装置 |
CN104820400A (zh) * | 2015-04-18 | 2015-08-05 | 桂林鸿程机电设备有限公司 | 一种三维焊接机器人混合控制方法 |
CN105180834A (zh) * | 2015-05-28 | 2015-12-23 | 华中科技大学 | 一种叶片进排气边三维非接触式测量装置 |
CN105583825A (zh) * | 2016-03-14 | 2016-05-18 | 陈杨 | 一种工业机器人轨迹检测装置 |
Cited By (35)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107053176A (zh) * | 2017-04-09 | 2017-08-18 | 北京工业大学 | 一种六自由度机器人末端空间曲线轨迹的误差建模方法 |
CN107272015A (zh) * | 2017-07-05 | 2017-10-20 | 同济大学 | 高精度视觉导引激光跟踪方法 |
CN107391340A (zh) * | 2017-07-21 | 2017-11-24 | 郑州云海信息技术有限公司 | 一种整机柜服务器节点热拔插系统和控制方法 |
CN107391340B (zh) * | 2017-07-21 | 2020-10-20 | 苏州浪潮智能科技有限公司 | 一种整机柜服务器节点热拔插系统和控制方法 |
CN109596125A (zh) * | 2017-09-30 | 2019-04-09 | 北京柏惠维康科技有限公司 | 一种机器人的空间坐标系转换关系的确定方法和装置 |
CN109596125B (zh) * | 2017-09-30 | 2022-03-11 | 北京柏惠维康科技有限公司 | 一种机器人的空间坐标系转换关系的确定方法和装置 |
CN107803855A (zh) * | 2017-12-13 | 2018-03-16 | 大连四达高技术发展有限公司 | 一种机器人自动定位机构 |
CN108858218B (zh) * | 2018-06-28 | 2021-10-12 | 北京航星机器制造有限公司 | 一种适应自动物流的机械臂制孔装置及制孔方法 |
CN108858218A (zh) * | 2018-06-28 | 2018-11-23 | 北京航星机器制造有限公司 | 一种适应自动物流的机械臂制孔装置及制孔方法 |
CN109186969A (zh) * | 2018-07-28 | 2019-01-11 | 西安交通大学 | 一种伺服进给运动动态性能视觉检测方法 |
CN109296024B (zh) * | 2018-11-30 | 2023-04-07 | 徐州市产品质量监督检验中心 | 一种无人挖掘机采装位姿精度检测方法 |
CN109296024A (zh) * | 2018-11-30 | 2019-02-01 | 徐州市产品质量监督检验中心 | 一种无人挖掘机采装位姿精度检测方法 |
CN109514557A (zh) * | 2018-12-13 | 2019-03-26 | 北京炎凌嘉业机电设备有限公司 | 一种3d视觉机器人轨迹修正系统 |
CN109884590A (zh) * | 2019-03-28 | 2019-06-14 | 河北工业大学 | 一种工业机器人轨迹精度检测装置及方法 |
CN109884590B (zh) * | 2019-03-28 | 2024-05-24 | 湖南第一师范学院 | 一种工业机器人轨迹精度检测装置及方法 |
CN111176304A (zh) * | 2020-03-18 | 2020-05-19 | 常州市贝叶斯智能科技有限公司 | 机器人运动底盘质检方法、装置、智能设备以及介质 |
CN111428626A (zh) * | 2020-03-23 | 2020-07-17 | 北京明略软件系统有限公司 | 一种移动物体的识别方法、装置、及存储介质 |
CN111428626B (zh) * | 2020-03-23 | 2023-05-23 | 北京明略软件系统有限公司 | 一种移动物体的识别方法、装置、及存储介质 |
CN111805531A (zh) * | 2020-06-30 | 2020-10-23 | 同济大学 | 一种管道内窥机器人 |
CN111975780A (zh) * | 2020-08-25 | 2020-11-24 | 厦门众合天元科技有限公司 | 一种工业机器人运动轨迹设定装置及其使用方法 |
CN111975780B (zh) * | 2020-08-25 | 2021-08-17 | 厦门众合天元科技有限公司 | 一种工业机器人运动轨迹设定装置及其使用方法 |
CN113977558A (zh) * | 2021-11-29 | 2022-01-28 | 湖南交通职业技术学院 | 一种可视化动态显示并联机器人末端轨迹的装置及方法 |
CN114563982B (zh) * | 2022-01-24 | 2023-05-09 | 中铁九桥工程有限公司 | 一种移动设备在圆管上运动轨迹的控制方法 |
CN114563982A (zh) * | 2022-01-24 | 2022-05-31 | 中铁九桥工程有限公司 | 一种移动设备在圆管上运动轨迹的控制方法 |
CN114454216B (zh) * | 2022-03-07 | 2023-10-10 | 云鲸智能(深圳)有限公司 | 机器人的精度检测方法、装置、机器人及存储介质 |
CN114454216A (zh) * | 2022-03-07 | 2022-05-10 | 云鲸智能(深圳)有限公司 | 机器人的精度检测方法、装置、机器人及存储介质 |
CN114454177A (zh) * | 2022-03-15 | 2022-05-10 | 浙江工业大学 | 一种基于双目立体视觉的机器人末端位置补偿方法 |
CN115752321A (zh) * | 2022-11-09 | 2023-03-07 | 中山大学 | 医疗机器人运动轨迹测量比对方法及计算机可读存储介质 |
CN115847427A (zh) * | 2023-02-07 | 2023-03-28 | 成都秦川物联网科技股份有限公司 | 双重识别的协作机器人监测工业物联网系统及其控制方法 |
CN115847427B (zh) * | 2023-02-07 | 2024-07-16 | 成都秦川物联网科技股份有限公司 | 双重识别的协作机器人监测工业物联网系统及其控制方法 |
CN116197918A (zh) * | 2023-05-05 | 2023-06-02 | 北京华晟经世信息技术股份有限公司 | 基于动作记录分析的机械手控制系统 |
CN118123849A (zh) * | 2024-05-08 | 2024-06-04 | 鹏城实验室 | 机器人轨迹控制方法、装置、设备及存储介质 |
CN118123849B (zh) * | 2024-05-08 | 2024-07-02 | 鹏城实验室 | 机器人轨迹控制方法、装置、设备及存储介质 |
CN118238154A (zh) * | 2024-05-28 | 2024-06-25 | 深圳市欣茂鑫实业有限公司 | 一种用于自动上料的机械臂控制方法及系统 |
CN118238154B (zh) * | 2024-05-28 | 2024-08-02 | 深圳市欣茂鑫实业有限公司 | 一种用于自动上料的机械臂控制方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN106541419B (zh) | 2019-01-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106541419B (zh) | 一种机器人轨迹误差的测量方法 | |
CN106546170B (zh) | 一种机器人运动轨迹关键点误差测量方法 | |
CN102607457B (zh) | 基于惯性导航技术的大尺寸三维形貌测量装置及方法 | |
CN102937418B (zh) | 一种扫描式物体表面三维形貌测量方法及装置 | |
US7145647B2 (en) | Measurement of spatial coordinates | |
CN102927018B (zh) | 一种用于离心泵piv相机测准及调节装置及其方法 | |
CN111811395B (zh) | 基于单目视觉的平面位姿动态测量方法 | |
CN103900489A (zh) | 一种线激光扫描三维轮廓测量方法及装置 | |
CN107272015A (zh) | 高精度视觉导引激光跟踪方法 | |
CN104296655B (zh) | 一种激光跟踪仪像旋公式初始角的标定方法 | |
KR101216706B1 (ko) | 예인수조용 입자영상유속계 시스템 | |
US20180283851A1 (en) | Motion detection device and three-dimensional shape measurement device using same | |
CN110017852B (zh) | 一种导航定位误差测量方法 | |
CN108253939A (zh) | 可变视轴单目立体视觉测量方法 | |
CN102944188A (zh) | 一种点扫描三维形貌测量系统标定方法 | |
CN106403900A (zh) | 飞行物追踪定位系统及方法 | |
CN208818162U (zh) | 定位机器人 | |
CN105423954A (zh) | 一种基于视觉测量的柔壁喷管的测量方法 | |
CN206339207U (zh) | 一种轨迹准确度重复性测量仪 | |
CN110686595A (zh) | 非正交轴系激光全站仪的激光束空间位姿标定方法 | |
CN105716547A (zh) | 一种机械工件平面度快速测量装置及方法 | |
CN108413865A (zh) | 基于三维测量和坐标系转换的二次反射镜面型检测方法 | |
CN202304767U (zh) | 一种基于多传感器的三维轮廓测量装置 | |
CN102155918B (zh) | 一种柔性喷管五自由度运动参数的测量装置和方法 | |
CN110030940B (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190125 |