CN106864776B - 一种基于对接环的捕获目标卫星的方法与系统 - Google Patents
一种基于对接环的捕获目标卫星的方法与系统 Download PDFInfo
- Publication number
- CN106864776B CN106864776B CN201710040120.8A CN201710040120A CN106864776B CN 106864776 B CN106864776 B CN 106864776B CN 201710040120 A CN201710040120 A CN 201710040120A CN 106864776 B CN106864776 B CN 106864776B
- Authority
- CN
- China
- Prior art keywords
- angular velocity
- mechanical arm
- docking ring
- point
- relative position
- 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
- 238000003032 molecular docking Methods 0.000 abstract description 58
- 239000013598 vector Substances 0.000 abstract description 40
- 238000000034 method Methods 0.000 abstract description 19
- 210000000078 claw Anatomy 0.000 abstract description 5
- 239000011159 matrix material Substances 0.000 description 23
- 238000004364 calculation method Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 6
- 238000005516 engineering process Methods 0.000 description 5
- 230000006872 improvement Effects 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 230000010354 integration Effects 0.000 description 3
- 238000005259 measurement Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 230000009466 transformation Effects 0.000 description 2
- 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 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000007667 floating Methods 0.000 description 1
- 239000000446 fuel Substances 0.000 description 1
- 230000014509 gene expression Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 230000008439 repair process Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64G—COSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
- B64G4/00—Tools specially adapted for use in space
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64G—COSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
- B64G4/00—Tools specially adapted for use in space
- B64G2004/005—Robotic manipulator systems for use in space
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Aviation & Aerospace Engineering (AREA)
- Manipulator (AREA)
Abstract
本发明公开了一种基于对接环的捕获目标卫星的方法,包括:系统控制手眼相机测量得到空间机器人机械臂末端点与对接环中心点的位姿偏差参数;系统求得所述空间机器人机械臂末端点与对接环中心点的相对位置矢量在对接环面的投影,得出距离机械臂末端最近的抓捕点,并得到相对位姿偏差。本发明还公开了一种捕获目标卫星的系统,系统判断所述相对位姿偏差与设定的临界值大小,以及当不满足临界值时通过规划机械臂末端的运动,进而控制机械臂末端的手爪闭合,进而实现目标的捕获,且达到时时接近的都是对接环上离得最近的一个点,即“动态最近点”,实现快速抓捕的目的,广泛应用于空间合作/非合作卫星抓捕的技术领域。
Description
技术领域
本发明涉及空间合作/非合作卫星抓捕的技术领域,具体为一种基于对接环的捕获目标卫星的方法与系统。
背景技术
自主在轨服务是利用安装在服务卫星上的机械臂捕获轨道上的卫星进行操作。比如,辅助未成功进入预定轨道的卫星入轨;对发生故障的卫星进行维修、更换故障元器件、辅助将未展开机构展开;对燃料耗尽但其他系统工作正常的卫星进行燃料加注,延长使用寿命;将废弃卫星和空间碎片进行回收或者送入坟墓轨道。对于大多数自主在轨服务任务而言,首先要解决操作目标的抓捕问题。目前,在轨捕获技术仍然是航天领域具有挑战性的研究方向,内容涉及空间交会对接技术、目标测量技术、高性能计算机软硬件技术、先进GNC技术等,具有跨学科、跨领域的特点。以往的目标捕获研究,如日本的ETS-VII项目、美国的“轨道快车项目”等主要针对的是合作目标卫星的捕获。而对于非合作目标,尤其是运动复杂的非合作目标的捕获,研究较少。非合作目标抓捕面临许多难题,如:
(1)非合作目标运动形式较为复杂;
(2)非合作目标抓捕的过程中,目标卫星与服务的空间机器人之间的相对位置测量较为困难;
(3)非合作目标上未安装用于目标捕获的手柄装置。
因此,有必要进行改进。
发明内容
为了解决上述技术问题,本发明的目的是提供一种可快速抓捕目标卫星的方法与系统。
本发明所采用的技术方案是:
本发明提供一种基于对接环的捕获目标卫星的方法,包括以下步骤:
S1、系统控制手眼相机测量得到空间机器人机械臂末端点与对接环中心点的相对位置矢量Δr和相对姿态θ;
S2、系统根据对接环面的法向量zA,求得所述空间机器人机械臂末端点与对接环中心点的相对位置矢量Δr在对接环面的投影bn,得出距离机械臂末端最近的抓捕点,并计算得到机械臂末端点到对接环“动态最近点”的相对位置偏差Δp和相对姿态偏差ΔO;
S3、系统判断所述相对位置偏差模||Δp||和相对姿态偏差模||ΔO||是否小于设定的临界值εp和εo,若二者均小于临界值,则机械臂末端的手爪闭合,实现目标的捕获;反之,则进入S4;
S4、系统根据所得的相对位置偏差Δp和相对姿态偏差ΔO,根据以下公式规划机械臂末端的线速度与角速度,
其中,Kp、Kv分别为规划的比例、微分参数;分别表示目标卫星的线速度与角速度;
S5、系统根据S4中得到的机械臂末端规划的线速度和角速度,通过雅可比矩阵反解,求得机械臂各关节的角速度;
S6、系统根据所得的关节角速度,驱动机械臂各关节运动,直至机械臂末端与“动态最近点”的相对位置偏差模和相对姿态偏差模小于设定的临界值εp和εo,手抓闭合并完成目标卫星的捕获。
作为该技术方案的改进,所述机械臂末端点到对接环“动态最近点”的相对位置偏差Δp的计算公式为其中,R为对接环的半径。
作为该技术方案的改进,所述机械臂末端的线速度满足
其中,ts为设定的平滑启动时间,ved、vem分别表示规划的线速度以及设定的最大线速度值。
进一步地,机械臂末端的角速度满足
其中,ωed、ωem分别表示规划的角速度以及设定的最大角速度值。
进一步地,所述临界值εp小于等于50毫米,εo小于等于5°。
进一步地,所述关节角的角速度求解公式为:
其中,表示关节角速度,表示在末端坐标系中,关于机械臂参数的雅克比矩阵,分别表示在末端坐标系中,机械臂末端相对于基座的速度与角速度。
另一方面,本发明还提供一种基于对接环的捕获目标卫星的系统,包括:
第一模块,用于执行步骤S1、系统控制手眼相机测量得到空间机器人机械臂末端点与对接环中心点的相对位置矢量Δr和相对姿态θ;
第二模块,用于执行步骤S2、系统根据对接环面的法向量zA,求得所述空间机器人机械臂末端点与对接环中心点的相对位置矢量Δr在对接环面的投影bn,得出距离机械臂末端最近的抓捕点,并计算得到机械臂末端点到对接环“动态最近点”的相对位置偏差Δp和相对姿态偏差ΔO;
第三模块,用于执行步骤S3、系统判断所述相对位置偏差模||Δp||和相对姿态偏差模||ΔO||是否小于设定的临界值εp和εo,若二者均小于临界值,则机械臂末端的手爪闭合,实现目标的捕获;反之,则进入S4;
第四模块,用于执行步骤S4、系统根据所得的相对位置偏差Δp和相对姿态偏差ΔO,根据以下公式规划机械臂末端的线速度与角速度,
其中,Kp、Kv分别为规划的比例、微分参数;分别表示目标卫星的线速度与角速度;
第五模块,用于执行步骤S5、系统根据S4中得到的机械臂末端规划的线速度和角速度,通过雅可比矩阵反解,求得机械臂各关节的角速度;
第六模块,用于执行步骤S6、系统根据所得的关节角速度,驱动机械臂各关节运动,直至机械臂末端与“动态最近点”的相对位置偏差模和相对姿态偏差模小于设定的临界值εp和εo,手抓闭合并完成目标卫星的捕获。
本发明的有益效果是:
本发明提供的一种可快速抓捕目标卫星的方法与系统,通过利用空间机器人机械臂末端点与对接环中心点的相对位置偏差Δr和二者z轴之间的夹角θ,通过运算处理得到机械臂末端点到对接环“动态最近点”的相对位置偏差Δp和相对姿态偏差ΔO,并将二者的模||Δp||和||ΔO||与设定的临界值εp和εo的对比判断,以及当不满足临界值时通过规划机械臂末端的运动,进而控制机械臂末端的手爪闭合,进而实现目标的捕获,且达到时时接近的都是对接环上离得最近的一个点,即“动态最近点”,实现快速抓捕的目的。
附图说明
下面结合附图对本发明的具体实施方式作进一步说明:
图1是本发明一实施例的机械臂末端与目标星对接环相对位姿示意图;
图2是本发明一实施例的服务卫星捕获目标卫星结构模型示意图;
图3是本发明一实施例的基于对接环“动态最近点”捕获的自主轨迹规划流程示意图;
图4是本发明一实施例的服务卫星捕获目标卫星简化模型示意图。
具体实施方式
需要说明的是,在不冲突的情况下,本申请中的实施例及实施例中的特征可以相互组合。
本发明提供一种基于对接环的捕获目标卫星的方法,包括以下步骤:
S1、系统控制手眼相机测量得到空间机器人机械臂末端点与对接环中心点的相对位置矢量Δr和相对姿态θ;
S2、系统根据对接环面的法向量zA,求得所述空间机器人机械臂末端点与对接环中心点的相对位置矢量Δr在对接环面的投影bn,得出距离机械臂末端最近的抓捕点,并计算得到机械臂末端点到对接环“动态最近点”的相对位置偏差Δp和相对姿态偏差ΔO;
S3、系统判断所述相对位置偏差模||Δp||和相对姿态偏差模||ΔO||是否小于设定的临界值εp和εo,若二者均小于临界值,则机械臂末端的手爪闭合,实现目标的捕获;反之,则进入S4;
S4、系统根据所得的相对位置偏差Δp和相对姿态偏差ΔO,根据以下公式规划机械臂末端的线速度与角速度,
其中,Kp、Kv分别为规划的比例、微分参数;分别表示目标卫星的线速度与角速度;
S5、系统根据S4中得到的机械臂末端规划的线速度和角速度,通过雅可比矩阵反解,求得机械臂各关节的角速度;
S6、系统根据所得的关节角速度,驱动机械臂各关节运动,直至机械臂末端与“动态最近点”的相对位置偏差模和相对姿态偏差模小于设定的临界值εp和εo,手抓闭合并完成目标卫星的捕获。
作为该技术方案的改进,所述机械臂末端点到对接环“动态最近点”的相对位置偏差Δp的计算公式为其中,R为对接环的半径。
作为该技术方案的改进,所述机械臂末端的线速度满足
其中,ts为设定的平滑启动时间,ved、vem分别表示规划的线速度以及设定的最大线速度值。
进一步地,机械臂末端的角速度满足
其中,ωed、ωem分别表示规划的角速度以及设定的最大角速度值。
进一步地,所述临界值εp小于等于50毫米,εo小于等于5°。
进一步地,所述关节角的角速度求解公式为:
其中,表示关节角速度,表示在末端坐标系中,关于机械臂参数的雅克比矩阵,分别表示在末端坐标系中,机械臂末端相对于基座的速度与角速度。
另一方面,本发明还提供一种基于对接环的捕获目标卫星的系统,包括:
第一模块,用于执行步骤S1、系统控制手眼相机测量得到空间机器人机械臂末端点与对接环中心点的相对位置矢量Δr和相对姿态θ;
第二模块,用于执行步骤S2、系统根据对接环面的法向量zA,求得所述空间机器人机械臂末端点与对接环中心点的相对位置矢量Δr在对接环面的投影bn,得出距离机械臂末端最近的抓捕点,并计算得到机械臂末端点到对接环“动态最近点”的相对位置偏差Δp和相对姿态偏差ΔO;
第三模块,用于执行步骤S3、系统判断所述相对位置偏差模||Δp||和相对姿态偏差模||ΔO||是否小于设定的临界值εp和εo,若二者均小于临界值,则机械臂末端的手爪闭合,实现目标的捕获;反之,则进入S4;
第四模块,用于执行步骤S4、系统根据所得的相对位置偏差Δp和相对姿态偏差ΔO,根据以下公式规划机械臂末端的线速度与角速度,
其中,Kp、Kv分别为规划的比例、微分参数;分别表示目标卫星的线速度与角速度;
第五模块,用于执行步骤S5、系统根据S4中得到的机械臂末端规划的线速度和角速度,通过雅可比矩阵反解,求得机械臂各关节的角速度;
第六模块,用于执行步骤S6、系统根据所得的关节角速度,驱动机械臂各关节运动,直至机械臂末端与“动态最近点”的相对位置偏差模和相对姿态偏差模小于设定的临界值εp和εo,手抓闭合并完成目标卫星的捕获。
参照图1,是本发明一实施例的机械臂末端与目标星对接环相对位姿示意图。
用于空间合作/非合作卫星对接环“动态最近点”捕获的方法,其主要步骤包括:
S1:根据安装在机械臂上的手眼相机测量,得到末端点与对接环中心点的相对位置矢量Δr和二者z轴之间的夹角(相对姿态)θ。
S2:根据对接环面的法相量zA,求出相对位置矢量Δr在对接环面的投影bn。由于对接环的半径R的大小已知,可以得出距离机械臂末端最近的抓捕点T,从而得到机械臂末端点到对接环“动态最近点”的相对位置偏差Δp,并通过计算求得机械臂末端点到对接环“动态最近点”的相对姿态偏差ΔO。
S3:判断相对位置偏差模||Δp||和姿态偏差模||ΔO||是否小于初始时设定的临界值εp和εo,若二者小于临界值,则机械臂末端的手爪闭合,实现目标的捕获;反之,进行下一步。
S4:由得到的相对位置偏差Δp和相对姿态偏差ΔO,根据公式(1)规划机械臂末端的线速度与角速度,
S5:根据S4中得到的机械臂末端规划的线速度和角速度,通过雅可比矩阵反解,求得机械臂各个关节角的角速度。
S6:由S5中得到的关节角速度,驱动机械臂各个关节运动,直至机械臂末端与“动态最近点”的相对位姿小于初始时设定的临界值εp和εo,手抓闭合,完成目标卫星的捕获。
参照图2,是本发明一实施例的服务卫星捕获目标卫星结构模型示意图。该服务卫星由基座和搭载基座上的空间机械臂组成,机械臂是由连杆和基座组成,在机械臂末端安装有测量目标卫星位姿信息的手眼相机。机械臂第一关节J1连接空间机器人基座和连杆B1,第二关节J2连接连杆B1和连杆B2,以此类推。
参照图3,是本发明的用于空间合作/非合作卫星对接环“动态最近点”捕获的方法的一个实施例的流程示意图。其包括步骤:
S1:根据安装在机械臂上的手眼相机测量,得到末端点与对接环中心点的相对位置矢量Δr和二者z轴之间的夹角θ。
S2:根据对接环面的法相量zA,求出相对位置矢量Δr在对接环面的投影bn。由于对接环的半径R的大小已知,可以得出距离机械臂末端最近的抓捕点T,从而得到机械臂末端点到对接环“动态最近点”的相对位置偏差Δp,通过计算求得机械臂末端点到对接环“动态最近点”的相对姿态偏差ΔO。
S3:判断相对位置偏差模||Δp||和姿态偏差模||ΔO||是否小于初始时设定的临界值εp和εo。若二者小于临界值,则机械臂末端的手爪闭合,实现目标的捕获。反之,进行下一步。所述临界值εp小于等于50毫米,εo小于等于5°。
S4:由得到的相对位置偏差Δp和相对姿态偏差ΔO,根据公式(1)规划机械臂末端的线速度与角速度。
S5:根据S4中得到的机械臂末端规划的线速度和角速度,通过雅可比矩阵反解,求得机械臂末端各个关节角的角速度。
S6:由S5中得到的关节角速度,驱动机械臂各个关节运动,直至机械臂末端与“动态最近点”的相对位姿模小于初始时设定的临界值εp和εo,手抓闭合,完成目标卫星的捕获。
参照图4,作为本发明的一个实施例的物理模型示意图,其中:
B0:刚体0,即空间机器人基座;
Bi(i=1,……,n):刚体i,机械臂的第i个连杆;
Ji:连接Bi-1和Bi的关节;
Ci:Bi的质心;
OI:惯性系原点;
Og:整个系统的质心;
∑I、∑E:分别表示惯性系、机械臂末端坐标系;
∑i:Bi的固连坐标系,Zi正向为Ji的旋转方向;
iAj∈R3×3:∑j相对于∑i的姿态变换矩阵,当i缺省时,表示∑j相对于惯性系的姿态变换矩阵,iAj=[inj,ioj,iaj];
E,O:分别为单位矩阵和零矩阵;
mi,M:mi是Bi的质量,且
iIj∈R3×3:在坐标系∑i表示下,刚体Bj相对于其质心的惯量矩阵,当i缺省时,表示在惯性系中的表示;
ki∈R3:表示Ji旋转方向的单位矢量;
ri∈R3:Ci的位置矢量;
rg∈R3:系统质心的位置矢量;
pi∈R3:Ji的位置矢量;
pe∈R3:机械臂末端的位置矢量;
ai,bi∈R3:在∑i下,分别为从Ji指向Ci,Ci指向Ji+1的位置矢量;
li∈R3:从∑i原点指向∑i+1原点的矢量,li=ai+bi;
分别表示∑i相对于∑j的线速度和角速度,在∑k中的表示,如果∑j或者∑k为惯性系,则可以省去相应的符号j或者k;
v0,ω0∈R3:基座的线速度和角速度;
vi,ωi∈R3:Bi的线速度和角速度;
θi:杆件i旋转的角度;
Θ∈Rn:关节角速度,即Θ=[θ1,…,θn];
Ψb∈R3:基座的姿态角,用z-y-x欧拉角表示,即Ψb=[αb,βb,γb]T
Ψ0:基座姿态欧拉角,按照x-y-z欧拉角的形式表示;
Ψe:空间机器人机械臂末端的姿态欧拉角,按照x-y-z欧拉角的形式表示。
对于矢量的表示方法,作如下规定:iv表示在固连于第i个刚体的坐标系∑i中的速度矢量,而当v没有左上标i的标识时表示在惯性坐标系中的速度矢量。为简化两个矢量叉乘计算,定义如下叉乘算子:若p=[x,y,z]T,则
因此,可以把两个向量的叉乘计算转化为一个矩阵乘上一个向量的运算,即
可以得出,空间机器人机械臂中的各个刚体质心的位置矢量如公式(3)所示。同理,空间机器人机械臂末端的位置矢量如公式(4)所示。
对公式(3)和公式(4)两边进行求导,由此可以得出空间机器人机械臂中的各个刚体质心的线速度以及空间机器人机械臂末端的线速度为,
根据空间机器人机械臂中各个刚体的运动规律,空间机器人机械臂中各个刚体质心的角速度和机械臂末端的角速度可以表示为,
将式(6)和式(8)写成矩阵的形式:
在公式(9)中,Jb代表和空间机器人基座的运动状态相关的雅克比矩阵。
在公式(10)中,块矩阵Jbv和块矩阵Jbw是Jb的分块矩阵,二者分别与空间机器人基座的线速度和角速度相关。
Jm代表与空间机器人机械臂运动状态相关的雅可比矩阵,其结果和与之对应的基座固定的机器人的雅可比矩阵一致,即:
类似地,Jmv、Jmw是Jm的分块矩阵。
当空间机器人系统处于自由漂浮状态且该系统不受任何外力/力矩作用时,系统的动量守恒。假设空间机器人系统在初始时线动量和角动量为零,那么整个空间机器人系统角动量不但符合公式(12)完整约束方程,而且还符合公式(13)的非完整约束方程:
将式(5)和式(7)带入到(13)可得
其中,
JRi=[k1,k2,...,ki,0,...,0]∈R3xn (17)
r0i=ri-r0 (18)
参照图3所示的规划流程,空间机械臂末端期望的线速度和期望的角速度按照公式(1)进行规划。
在公式(1)中,Kp、Kv分别为规划的比例、微分参数(对应于分解运动控制的PD参数);分别表示通过手眼相机测量出的目标的运动速度;Δp表示械臂末端相对于目标星对接环“动态最近点”的位置误差;ΔO表示机械臂末端相对于目标抓捕点的姿态指向误差。以下,详细说明Δp和ΔO的求算方法。
根据安装在机械臂末端的相机可以得出机械臂末端相对于对接环中心点的相对位置Δr和在空间机器人机械臂末端坐标系下表示的目标卫星对接环平面的法向量zA。机械臂末端坐标轴zE与对接环法向量zA轴之间的夹角为,
θ=arc cos(zE·zA) (19)
相对位置Δr在对接环平面法相向量zA的投影an如下式所示。
an=||Δr||zA cosθ (20)
相对位置Δr在对接环平面的投影bn如下式所示。
bn=Δr-an (21)
根据公式(21)可以得出,空间机器人机械臂末端点和对接环上捕获点的相对位置矢量。
该方法中,仅期望机械臂末端坐标轴z轴与捕获点坐标轴z轴重合,空间机器人末端点和对接环上捕获点的相对姿态矢量如公式(23)所示。
ΔOe=en·θ (23)
在公式(23)中,en为旋转轴,且en=zE×zA,en=[ex ey ez]T。
根据公式(22)和公式(23),可以得出了空间机器人末端点和目标星对接环“动态最近点”的相对位置偏差与相对姿态偏差。
在实际应用中,机械臂末端的运动速度不能超过设计的最大值,且需要保证机械臂在启动过程的平稳性。因而对空间机器人机械臂末端运动速度做下述处理。
对于机械臂末端的线速度,
对于机械臂末端的角速度,
其中,ts为设定的平滑启动时间,如ts=2s。这样,可以确保机械臂在0-ts内平滑启动到合成速度最大值为vem(方向矢量仍然是vem);而当启动后的抓捕过程中,机械臂末端速度≥vem时,将其限制在vem(方向矢量仍然是vem),这样的处理在工程中非常重要。
将公式(12)代入公式(9)中,可以将空间机械臂的末端线速度和角速度如下表示:
当机械臂末端线速度和角速度在机械臂末端坐标系中表示时,可以得出,
在公式(29)中,
在公式(33)中,表示与空间机器人机械臂相关的雅可比矩阵,在末端坐标系表示。当采用x-y-z欧拉角表示空间机器人基座的姿态时,空间机器人基座的欧拉角的速率与基座姿态角速度的关系表达式如下。
对上式进行数值积分,可以得出基座姿态的欧拉角表示。另外,空间机器人系统还满足角动量守恒方程,即
其中,
JTi=[k1×(ri-p1),k2×(ri-p2),...,ki×(ri-pi),0,...,0]∈R3xn (41)
r0g=rg-r0 (42)
由于Hs是非奇异的,因此:
对于公式(43)中的各个矢量,当同时在空间机器人机械臂的末端坐标系表示时,最终空间机器人基座的角速度如下。
通过搭载在基座上的姿态传感器,我们可以测量出空间机器人基座的姿态及姿态角速度,由此可得出如下关系。
公式(45)的左侧部分表示在末端坐标系下,空间机器人末端相对于空间机器人基座的相对速度,由此可以将公式(45)写为:
在公式(46)中,是在末端坐标系下,机械臂末端相对于空间机器人基座的线速度。在末端坐标系下,机械臂末端相对于空间机器人基座的角速度。而且二者同时在空间机器人末端坐标系中表示。
则针对动力学奇异回避问题转换为针对方程(47)的运动学奇异回避问题,而雅可比矩阵只包含空间机器人运动学参数,具有简单的表达式。当关节角速度确定后,可采用数值积分的方法得到关节角,其中一种最简单的方法为一步积分法。
通过公式(48)得到的关节角速度,驱动机械臂各个关节运动,直至机械臂末端与“动态最近点”的相对位姿小于初始时设定的临界值εp和εo,手爪闭合,完成目标卫星的捕获。
本发明提供的一种可快速抓捕目标卫星的方法与系统,通过利用空间机器人机械臂末端点与“动态最近点”的相对位置偏差Δp和相对姿态偏差ΔO与设定的临界值εp和εo的对比判断,以及当不满足临界值时通过规划机械臂末端的运动,进而控制机械臂末端的手爪闭合,进而实现目标的捕获,且达到时时接近的都是对接环上离得最近的一个点,即“动态最近点”,实现快速抓捕的目的。
以上是对本发明的较佳实施进行了具体说明,但本发明创造并不限于所述实施例,熟悉本领域的技术人员在不违背本发明精神的前提下还可做出种种的等同变形或替换,这些等同的变形或替换均包含在本申请权利要求所限定的范围内。
Claims (7)
1.一种基于对接环的捕获目标卫星的方法,其特征在于,包括以下步骤:
S1、系统控制手眼相机测量得到空间机器人机械臂末端点与对接环中心点的相对位置矢量Δr和相对姿态θ;
S2、系统根据对接环面的法向量zA,求得所述空间机器人机械臂末端点与对接环中心点的相对位置矢量Δr在对接环面的投影bn,得出距离机械臂末端最近的抓捕点,并计算得到机械臂末端点到对接环“动态最近点”的相对位置偏差Δp和相对姿态偏差ΔO;
S3、系统判断所述相对位置偏差模||Δp||和相对姿态偏差模||ΔO||是否小于设定的临界值εp和εo,若二者均小于临界值,则机械臂末端的手爪闭合,实现目标的捕获;反之,则进入S4;
S4、系统根据所得的相对位置偏差Δp和相对姿态偏差ΔO,根据以下公式规划机械臂末端的线速度与角速度,
其中,Kp、Kv分别为规划的比例、微分参数;分别表示目标卫星的线速度与角速度;
S5、系统根据S4中得到的机械臂末端规划的线速度和角速度,通过雅可比矩阵反解,求得机械臂各关节的角速度;
S6、系统根据所得的关节角速度,驱动机械臂各关节运动,直至机械臂末端与“动态最近点”的相对位置偏差模和相对姿态偏差模小于设定的临界值εp和εo,手抓闭合并完成目标卫星的捕获。
2.根据权利要求1所述的基于对接环的捕获目标卫星的方法,其特征在于:所述机械臂末端点到对接环“动态最近点”的相对位置偏差Δp的计算公式为其中,R为对接环的半径。
3.根据权利要求2所述的基于对接环的捕获目标卫星的方法,其特征在于:所述机械臂末端的线速度满足
其中,ts为设定的平滑启动时间,ved、vem分别表示规划的线速度以及设定的最大线速度值,vex、vey、vez分别表示规划的机械臂末端的线速度在机械臂末端坐标系x-y-z三个方向的分量,t是时间,表示机械臂末端线速度,即机械臂末端实际的线速度。
4.根据权利要求3所述的基于对接环的捕获目标卫星的方法,其特征在于:机械臂末端的角速度满足
其中,ωed、ωem分别表示规划的角速度以及设定的最大角速度值,ωex、ωey、ωez分别表示规划的机械臂末端的角速度在机械臂末端坐标系x-y-z三个方向的分量,t是时间,表示机械臂末端角速度,即机械臂末端实际的角速度。
5.根据权利要求4所述的基于对接环的捕获目标卫星的方法,其特征在于:所述临界值εp小于等于50毫米,εo小于等于5°。
6.根据权利要求5所述的基于对接环的捕获目标卫星的方法,其特征在于,所述关节角的角速度求解公式为:
其中,表示关节角速度,表示在末端坐标系中,关于机械臂参数的雅克比矩阵,分别表示在末端坐标系中,机械臂末端相对于基座的速度与角速度。
7.一种基于对接环的捕获目标卫星的系统,其特征在于,包括:
第一模块,用于执行步骤S1、系统控制手眼相机测量得到空间机器人机械臂末端点与对接环中心点的相对位置矢量Δr和相对姿态θ;
第二模块,用于执行步骤S2、系统根据对接环面的法向量zA,求得所述空间机器人机械臂末端点与对接环中心点的相对位置矢量Δr在对接环面的投影bn,得出距离机械臂末端最近的抓捕点,并计算得到机械臂末端点到对接环“动态最近点”的相对位置偏差Δp和相对姿态偏差ΔO;
第三模块,用于执行步骤S3、系统判断所述相对位置偏差模||Δp||和相对姿态偏差模||ΔO||是否小于设定的临界值εp和εo,若二者均小于临界值,则机械臂末端的手爪闭合,实现目标的捕获;反之,则进入S4;
第四模块,用于执行步骤S4、系统根据所得的相对位置偏差Δp和相对姿态偏差ΔO,根据以下公式规划机械臂末端的线速度与角速度,
其中,Kp、Kv分别为规划的比例、微分参数;分别表示目标卫星的线速度与角速度;
第五模块,用于执行步骤S5、系统根据S4中得到的机械臂末端规划的线速度和角速度,通过雅可比矩阵反解,求得机械臂各关节的角速度;
第六模块,用于执行步骤S6、系统根据所得的关节角速度,驱动机械臂各关节运动,直至机械臂末端与“动态最近点”的相对位置偏差模和相对姿态偏差模小于设定的临界值εp和εo,手抓闭合并完成目标卫星的捕获。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710040120.8A CN106864776B (zh) | 2017-01-18 | 2017-01-18 | 一种基于对接环的捕获目标卫星的方法与系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710040120.8A CN106864776B (zh) | 2017-01-18 | 2017-01-18 | 一种基于对接环的捕获目标卫星的方法与系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106864776A CN106864776A (zh) | 2017-06-20 |
CN106864776B true CN106864776B (zh) | 2019-04-12 |
Family
ID=59158409
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710040120.8A Active CN106864776B (zh) | 2017-01-18 | 2017-01-18 | 一种基于对接环的捕获目标卫星的方法与系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106864776B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108492333B (zh) * | 2018-03-30 | 2021-12-03 | 哈尔滨工业大学 | 基于星箭对接环图像信息的航天器姿态估计方法 |
CN108594757B (zh) * | 2018-05-15 | 2021-01-01 | 南京旭上数控技术有限公司 | 一种基于位置和姿态约束的机器人小线段前瞻规划方法 |
CN109284768B (zh) * | 2018-07-25 | 2022-04-01 | 西北工业大学 | 一种空间抓捕过程不确定性重构和预测方法 |
CN111591474B (zh) * | 2020-02-28 | 2021-12-07 | 上海航天控制技术研究所 | 一种航天器在轨操作系统对准式手眼标定方法 |
CN111891411B (zh) * | 2020-07-16 | 2022-04-12 | 北京控制工程研究所 | 一种三爪式通用对接机构抓捕控制方法 |
CN111942621B (zh) * | 2020-07-17 | 2022-03-04 | 北京控制工程研究所 | 一种基于多任务学习的在轨自主加注控制方法及系统 |
CN112060088A (zh) * | 2020-08-28 | 2020-12-11 | 华中科技大学 | 一种变时延条件下非合作目标精确抓捕遥操作方法 |
CN115416030B (zh) * | 2022-02-24 | 2023-06-20 | 哈尔滨工业大学 | 一种机械臂远距离跟踪翻滚卫星对接环的运动规划方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA2831692C (en) * | 2007-03-09 | 2015-01-27 | Macdonald, Dettwiler And Associates Inc. | Satellite refuelling system and method |
CN101733746A (zh) * | 2009-12-22 | 2010-06-16 | 哈尔滨工业大学 | 空间机器人非合作目标自主识别与捕获方法 |
CN103331759B (zh) * | 2013-06-28 | 2015-06-10 | 哈尔滨工业大学 | 用于空间大型机械臂末端执行器的大容差捕获机构 |
US9302793B2 (en) * | 2014-03-21 | 2016-04-05 | The Boeing Company | Spacecraft docking system |
CN104097791B (zh) * | 2014-06-24 | 2016-06-15 | 上海微小卫星工程中心 | 一种基于磁敏感器和星敏感器的全姿态捕获方法及其装置 |
-
2017
- 2017-01-18 CN CN201710040120.8A patent/CN106864776B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN106864776A (zh) | 2017-06-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106864776B (zh) | 一种基于对接环的捕获目标卫星的方法与系统 | |
CN109048890B (zh) | 基于机器人的协调轨迹控制方法、系统、设备及存储介质 | |
Zhang et al. | Robust control of an aerial manipulator based on a variable inertia parameters model | |
CN106166749B (zh) | 一种空间移动多臂机器人的移动轨迹规划方法 | |
Yashin et al. | Aerovr: Virtual reality-based teleoperation with tactile feedback for aerial manipulation | |
CN103009389B (zh) | 一种冗余空间机械臂在轨抓捕的轨迹规划方法 | |
Dong et al. | Position-based visual servo control of autonomous robotic manipulators | |
CN105116910B (zh) | 一种对地面点凝视成像的卫星姿态控制方法 | |
Yan et al. | Virtual-base modeling and coordinated control of a dual-arm space robot for target capturing and manipulation | |
CN106892137A (zh) | 一种空间非合作翻滚目标的大容差捕获方法及系统 | |
CN101733746A (zh) | 空间机器人非合作目标自主识别与捕获方法 | |
CN104015191B (zh) | 基于基座卫星角速度的空间机械臂工具坐标下的运动补偿方法 | |
CN109426147B (zh) | 捕获卫星后组合航天器的自适应增益调整控制方法 | |
CN108381553A (zh) | 一种用于空间非合作目标捕获的相对导航近距离跟踪方法及系统 | |
CN105539890A (zh) | 一种模拟空间机械臂捕获目标卫星的地面三维空间微重力的装置与方法 | |
CN105382843B (zh) | 一种抓捕末阶段机械臂与操作平台协调控制方法 | |
Xu et al. | A pose measurement method of a non-cooperative GEO spacecraft based on stereo vision | |
CN107450309B (zh) | 一种绳系组合体姿态机动协同控制方法 | |
CN106055901A (zh) | 一种自由漂浮空间机器人捕获翻滚目标的时机确定方法 | |
CN108436913A (zh) | 一种力协调的多臂机器人柔顺控制方法 | |
Ma et al. | Hand-eye servo and impedance control for manipulator arm to capture target satellite safely | |
Xu et al. | Modeling and planning of a space robot for capturing tumbling target by approaching the dynamic closest point | |
Zhao et al. | Minimum base disturbance control of free-floating space robot during visual servoing pre-capturing process | |
CN107696033A (zh) | 一种基于视觉测量的空间机械臂轨迹滚动规划方法 | |
CN110470297A (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 |