发明名称 一种自由漂浮空间机器人捕获翻滚目标的时机确定方法
摘要 本发明公开了一种自由漂浮空间机器人捕获翻滚目标的时机确定方法,包括建立了翻滚目标的运动方程;分析了自由漂浮空间机器人的工作空间;提出了确定最佳抓捕时机的三个准则;最后以实例验证了本发明提出的方法的有效性。本发明提出的确定时机的三个准则考虑了避免机械臂与翻滚目标碰撞,使末端执行器具有更好的操作灵活性和操作能力,以及快速地完成捕获,有益于未来空间机器人安全、可靠、快速地执行捕获翻滚目标任务。
申请公布号 CN106055901A 申请公布日期 2016.10.26
申请号 CN201610388972.1 申请日期 2016.06.02
申请人 西北工业大学 发明人 罗建军;宗立军;王明明;袁建平;朱战霞
分类号 G06F19/00(2011.01)I;G05D1/12(2006.01)I;G05D1/10(2006.01)I 主分类号 G06F19/00(2011.01)I
代理机构 西安通大专利代理有限责任公司 61200 代理人 陆万寿
主权项 一种自由漂浮空间机器人捕获翻滚目标的时机确定方法,其特征在于,包括以下步骤:步骤一、建立翻滚目标的运动方程;翻滚目标模型简化为长方刚体,坐标系ox<sub>t</sub>y<sub>t</sub>z<sub>t</sub>为定义在目标上的本体坐标系,其中,设o为目标质心,各坐标轴与通过质心的惯量主轴重合;定义刚体绕各轴的主惯量分别是Ι<sub>x</sub>,Ι<sub>y</sub>,Ι<sub>z</sub>;假设初始时刻本体坐标系和惯性坐标系重合,随后刚体以角速度ω旋转,其在本体坐标系三个轴向的分量分别记为ω<sub>x</sub>,ω<sub>y</sub>,ω<sub>z</sub>;假设目标在空间不受任何外力,则其姿态动力学方程表示为:<maths num="0001"><math><![CDATA[<mrow><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mi>I</mi><mi>x</mi></msub><msub><mover><mi>&omega;</mi><mo>&CenterDot;</mo></mover><mi>x</mi></msub><mo>-</mo><mrow><mo>(</mo><msub><mi>I</mi><mi>y</mi></msub><mo>-</mo><msub><mi>I</mi><mi>z</mi></msub><mo>)</mo></mrow><msub><mi>&omega;</mi><mi>y</mi></msub><msub><mi>&omega;</mi><mi>z</mi></msub><mo>=</mo><mn>0</mn></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>I</mi><mi>y</mi></msub><msub><mover><mi>&omega;</mi><mo>&CenterDot;</mo></mover><mi>y</mi></msub><mo>-</mo><mrow><mo>(</mo><msub><mi>I</mi><mi>z</mi></msub><mo>-</mo><msub><mi>I</mi><mi>x</mi></msub><mo>)</mo></mrow><msub><mi>&omega;</mi><mi>z</mi></msub><msub><mi>&omega;</mi><mi>x</mi></msub><mo>=</mo><mn>0</mn></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>I</mi><mi>z</mi></msub><msub><mover><mi>&omega;</mi><mo>&CenterDot;</mo></mover><mi>z</mi></msub><mo>-</mo><mrow><mo>(</mo><msub><mi>I</mi><mi>x</mi></msub><mo>-</mo><msub><mi>I</mi><mi>y</mi></msub><mo>)</mo></mrow><msub><mi>&omega;</mi><mi>x</mi></msub><msub><mi>&omega;</mi><mi>y</mi></msub><mo>=</mo><mn>0</mn></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001008198130000011.GIF" wi="1099" he="230" /></maths>使用四元数描述刚体的姿态变换矩阵,得到本体坐标系到惯性坐标系的姿态变换矩阵A(q)如式(2)所示:<maths num="0002"><math><![CDATA[<mrow><mi>A</mi><mrow><mo>(</mo><mi>q</mi><mo>)</mo></mrow><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><mrow><msubsup><mi>q</mi><mn>4</mn><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>q</mi><mn>1</mn><mn>2</mn></msubsup><mo>-</mo><msubsup><mi>q</mi><mn>2</mn><mn>2</mn></msubsup><mo>-</mo><msubsup><mi>q</mi><mn>3</mn><mn>2</mn></msubsup></mrow></mtd><mtd><mrow><mn>2</mn><msub><mi>q</mi><mn>4</mn></msub><msub><mi>q</mi><mn>3</mn></msub><mo>+</mo><mn>2</mn><msub><mi>q</mi><mn>1</mn></msub><msub><mi>q</mi><mn>2</mn></msub></mrow></mtd><mtd><mrow><mo>-</mo><mn>2</mn><msub><mi>q</mi><mn>4</mn></msub><msub><mi>q</mi><mn>2</mn></msub><mo>+</mo><mn>2</mn><msub><mi>q</mi><mn>1</mn></msub><msub><mi>q</mi><mn>3</mn></msub></mrow></mtd></mtr><mtr><mtd><mrow><mo>-</mo><mn>2</mn><msub><mi>q</mi><mn>4</mn></msub><msub><mi>q</mi><mn>3</mn></msub><mo>+</mo><mn>2</mn><msub><mi>q</mi><mn>1</mn></msub><msub><mi>q</mi><mn>2</mn></msub></mrow></mtd><mtd><mrow><msubsup><mi>q</mi><mn>4</mn><mn>2</mn></msubsup><mo>-</mo><msubsup><mi>q</mi><mn>1</mn><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>q</mi><mn>2</mn><mn>2</mn></msubsup><mo>-</mo><msubsup><mi>q</mi><mn>3</mn><mn>2</mn></msubsup></mrow></mtd><mtd><mrow><mn>2</mn><msub><mi>q</mi><mn>4</mn></msub><msub><mi>q</mi><mn>1</mn></msub><mo>+</mo><mn>2</mn><msub><mi>q</mi><mn>2</mn></msub><msub><mi>q</mi><mn>3</mn></msub></mrow></mtd></mtr><mtr><mtd><mrow><mn>2</mn><msub><mi>q</mi><mn>4</mn></msub><msub><mi>q</mi><mn>2</mn></msub><mo>+</mo><mn>2</mn><msub><mi>q</mi><mn>1</mn></msub><msub><mi>q</mi><mn>3</mn></msub></mrow></mtd><mtd><mrow><mo>-</mo><mn>2</mn><msub><mi>q</mi><mn>4</mn></msub><msub><mi>q</mi><mn>1</mn></msub><mo>+</mo><mn>2</mn><msub><mi>q</mi><mn>2</mn></msub><msub><mi>q</mi><mn>3</mn></msub></mrow></mtd><mtd><mrow><msubsup><mi>q</mi><mn>4</mn><mn>2</mn></msubsup><mo>-</mo><msubsup><mi>q</mi><mn>1</mn><mn>2</mn></msubsup><mo>-</mo><msubsup><mi>q</mi><mn>2</mn><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>q</mi><mn>3</mn><mn>2</mn></msubsup></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001008198130000012.GIF" wi="1453" he="214" /></maths>其中,<img file="FDA0001008198130000013.GIF" wi="470" he="63" />为表示姿态的单位四元数,前三个参数代表欧拉转轴的方向,第四个参数代表欧拉转角的大小,四元数各元素和本体坐标系下角速度的分量满足式(3)所示的姿态运动学方程:<maths num="0003"><math><![CDATA[<mrow><mfenced open = "[" close = "]"><mtable><mtr><mtd><msub><mover><mi>q</mi><mo>&CenterDot;</mo></mover><mn>1</mn></msub></mtd></mtr><mtr><mtd><msub><mover><mi>q</mi><mo>&CenterDot;</mo></mover><mn>2</mn></msub></mtd></mtr><mtr><mtd><msub><mover><mi>q</mi><mo>&CenterDot;</mo></mover><mn>3</mn></msub></mtd></mtr><mtr><mtd><msub><mover><mi>q</mi><mo>&CenterDot;</mo></mover><mn>4</mn></msub></mtd></mtr></mtable></mfenced><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mfenced open = "[" close = "]"><mtable><mtr><mtd><mn>0</mn></mtd><mtd><msub><mi>&omega;</mi><mi>z</mi></msub></mtd><mtd><mrow><mo>-</mo><msub><mi>&omega;</mi><mi>y</mi></msub></mrow></mtd><mtd><msub><mi>&omega;</mi><mi>x</mi></msub></mtd></mtr><mtr><mtd><mrow><mo>-</mo><msub><mi>&omega;</mi><mi>z</mi></msub></mrow></mtd><mtd><mn>0</mn></mtd><mtd><msub><mi>&omega;</mi><mi>x</mi></msub></mtd><mtd><msub><mi>&omega;</mi><mi>y</mi></msub></mtd></mtr><mtr><mtd><msub><mi>&omega;</mi><mi>y</mi></msub></mtd><mtd><mrow><mo>-</mo><msub><mi>&omega;</mi><mi>x</mi></msub></mrow></mtd><mtd><mn>0</mn></mtd><mtd><msub><mi>&omega;</mi><mi>z</mi></msub></mtd></mtr><mtr><mtd><mrow><mo>-</mo><msub><mi>&omega;</mi><mi>x</mi></msub></mrow></mtd><mtd><mrow><mo>-</mo><msub><mi>&omega;</mi><mi>y</mi></msub></mrow></mtd><mtd><mrow><mo>-</mo><msub><mi>&omega;</mi><mi>z</mi></msub></mrow></mtd><mtd><mn>0</mn></mtd></mtr></mtable></mfenced><mfenced open = "[" close = "]"><mtable><mtr><mtd><msub><mi>q</mi><mn>1</mn></msub></mtd></mtr><mtr><mtd><msub><mi>q</mi><mn>2</mn></msub></mtd></mtr><mtr><mtd><msub><mi>q</mi><mn>3</mn></msub></mtd></mtr><mtr><mtd><msub><mi>q</mi><mn>4</mn></msub></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001008198130000014.GIF" wi="1222" he="278" /></maths>假设抓捕目标上只有唯一的抓捕点,抓捕点在本体坐标系下的位置向量记为<img file="FDA0001008198130000015.GIF" wi="62" he="62" />则在惯性坐标系下,抓捕点的位置矢量表示为:<maths num="0004"><math><![CDATA[<mrow><msubsup><mi>r</mi><mi>a</mi><mi>I</mi></msubsup><mo>=</mo><mi>i</mi><mi>n</mi><mi>v</mi><mrow><mo>(</mo><mi>A</mi><mo>(</mo><mi>q</mi><mo>)</mo><mo>)</mo></mrow><msubsup><mi>r</mi><mi>a</mi><mi>t</mi></msubsup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001008198130000016.GIF" wi="1019" he="70" /></maths>其中,inv表示对矩阵求逆;通过辨识等方法得到目标的惯性参数、初始角速度以及抓捕点在目标本体系下的位置坐标,结合方程(1)~(4),给出抓捕点在惯性空间的运动轨迹;步骤二、计算自由漂浮空间机器人的工作空间;在卫星基座本体坐标系下,自由漂浮空间机器人的运动学方程如式(5)所示:<maths num="0005"><math><![CDATA[<mrow><mfenced open = "[" close = "]"><mtable><mtr><mtd><mrow><msub><mmultiscripts><mi>v</mi><mprescripts/><none/><mi>o</mi></mmultiscripts><mi>e</mi></msub></mrow></mtd></mtr><mtr><mtd><mrow><msub><mmultiscripts><mi>&omega;</mi><mprescripts/><none/><mi>o</mi></mmultiscripts><mi>e</mi></msub></mrow></mtd></mtr></mtable></mfenced><mo>=</mo><msub><mmultiscripts><mi>J</mi><mprescripts/><none/><mi>o</mi></mmultiscripts><mi>g</mi></msub><mrow><mo>(</mo><mi>&theta;</mi><mo>,</mo><msub><mi>m</mi><mi>i</mi></msub><mo>,</mo><msub><mi>I</mi><mi>i</mi></msub><mo>)</mo></mrow><mover><mi>&theta;</mi><mo>&CenterDot;</mo></mover><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001008198130000021.GIF" wi="1078" he="142" /></maths>其中,<sup>o</sup>v<sub>e</sub>,<sup>o</sup>ω<sub>e</sub>分别为末端执行器的线速度和角速度,<sup>o</sup>J<sub>g</sub>为空间机器人的广义雅克比矩阵,上标“o”表示在基座本体坐标系下的表示;求解式(6)所示的方程,得到对应矩阵<sup>o</sup>J<sub>g</sub>奇异自由漂浮空间机器人的奇异臂型:det[<sup>o</sup>J<sub>g</sub>]=0    (6)方程(6)的解对应机械臂关节空间下的一组超曲面Q<sub>s,i</sub>(i=1,2,...),这些超曲面为动力学奇异臂型相应关节转角q<sub>s</sub>的集合;采用自由漂浮空间机器人的虚拟机械臂模型,该模型下,机械臂末端的位置矢量表示为:<maths num="0006"><math><![CDATA[<mrow><msub><mi>p</mi><mi>e</mi></msub><mo>=</mo><msub><mi>r</mi><mi>g</mi></msub><mo>+</mo><msub><mmultiscripts><mi>A</mi><mprescripts/><none/><mi>I</mi></mmultiscripts><mi>C</mi></msub><mrow><mo>(</mo><msub><mmultiscripts><mover><mi>b</mi><mo>^</mo></mover><mprescripts/><none/><mi>C</mi></mmultiscripts><mn>0</mn></msub><mo>+</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mmultiscripts><mi>A</mi><mprescripts/><none/><mi>C</mi></mmultiscripts><mi>i</mi></msub><mo>(</mo><mrow><msub><mmultiscripts><mover><mi>a</mi><mo>^</mo></mover><mprescripts/><none/><mi>i</mi></mmultiscripts><mi>i</mi></msub><mo>+</mo><msub><mmultiscripts><mover><mi>b</mi><mo>^</mo></mover><mprescripts/><none/><mi>i</mi></mmultiscripts><mi>i</mi></msub></mrow><mo>)</mo><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001008198130000022.GIF" wi="1229" he="135" /></maths>其中,r<sub>g</sub>为系统质心的位置矢量,<sup>I</sup>A<sub>C</sub>和<sup>C</sup>A<sub>i</sub>分别为系统质心坐标系到惯性系和连杆本体坐标系到系统质心坐标系的变换矩阵,<img file="FDA0001008198130000023.GIF" wi="38" he="62" />和<img file="FDA0001008198130000024.GIF" wi="43" he="71" />称为虚拟连杆矢量;系统质心到机械臂末端的距离表示为:<maths num="0007"><math><![CDATA[<mrow><mi>R</mi><mrow><mo>(</mo><mi>q</mi><mo>)</mo></mrow><mo>=</mo><mo>|</mo><mo>|</mo><msub><mi>p</mi><mi>e</mi></msub><mo>-</mo><msub><mi>r</mi><mi>g</mi></msub><mo>|</mo><mo>|</mo><mo>=</mo><mo>|</mo><mo>|</mo><msub><mmultiscripts><mover><mi>b</mi><mo>^</mo></mover><mprescripts/><none/><mi>C</mi></mmultiscripts><mn>0</mn></msub><mo>+</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mmultiscripts><mi>A</mi><mprescripts/><none/><mi>C</mi></mmultiscripts><mi>i</mi></msub><mrow><mo>(</mo><msub><mmultiscripts><mover><mi>a</mi><mo>^</mo></mover><mprescripts/><none/><mi>i</mi></mmultiscripts><mi>i</mi></msub><mo>+</mo><msub><mmultiscripts><mover><mi>b</mi><mo>^</mo></mover><mprescripts/><none/><mi>i</mi></mmultiscripts><mi>i</mi></msub><mo>)</mo></mrow><mo>|</mo><mo>|</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001008198130000025.GIF" wi="1259" he="134" /></maths>符号“||·||”表示向量的长度,式(8)定义了以系统质心为球心、半径为R(q)的球面;每组使雅克比矩阵奇异的关节转角q<sub>s</sub>映射到惯性空间都对应式(8)定义的球面,基于同样的考虑,每个超曲面Q<sub>s,i</sub>映射到惯性空间将对应以系统质心为球心,由以下两个半径定义的球环:<maths num="0008"><math><![CDATA[<mrow><msub><mi>R</mi><mrow><mi>m</mi><mi>i</mi><mi>n</mi><mo>,</mo><mi>i</mi></mrow></msub><mo>=</mo><munder><mrow><mi>m</mi><mi>i</mi><mi>n</mi></mrow><mrow><mi>q</mi><mo>&Element;</mo><msub><mi>Q</mi><mrow><mi>s</mi><mo>,</mo><mi>i</mi></mrow></msub></mrow></munder><mi>R</mi><mrow><mo>(</mo><mi>q</mi><mo>)</mo></mrow><mo>,</mo><msub><mi>R</mi><mrow><mi>m</mi><mi>a</mi><mi>x</mi><mo>,</mo><mi>i</mi></mrow></msub><mo>=</mo><munder><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow><mrow><mi>q</mi><mo>&Element;</mo><msub><mi>Q</mi><mrow><mi>s</mi><mo>,</mo><mi>i</mi></mrow></msub></mrow></munder><mi>R</mi><mrow><mo>(</mo><mi>q</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>9</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001008198130000031.GIF" wi="1190" he="87" /></maths>记W<sub>i</sub>={R(q):R<sub>min,i</sub>≤R(q)≤R<sub>max,i</sub>},则超曲面Q<sub>s,i</sub>映射到惯性空间对应集合W<sub>i</sub>;如果方程(6)的解对应多个超曲面,则惯性空间有对应的多个集合W<sub>1</sub>,W<sub>2</sub>,...,W<sub>i</sub>,...,自由漂浮空间机器人的路径相关工作空间定义为<img file="FDA0001008198130000036.GIF" wi="582" he="62" />符号“∪”表示取集合的并集,路径无关工作空间由W<sub>PIW</sub>=W<sub>reach</sub>/W<sub>PDW</sub>计算得到,W<sub>reach</sub>表示空间机器人的可达工作空间,符号“A/B”表示集合A与B的差集;步骤三、确定捕获翻滚目标的最佳时机;设定三个准则用于确定空间机器人捕获翻滚目标的最佳时机;准则1:选择沿捕获方向,当抓捕点是翻滚目标上距离空间机器人系统最近的点时进行抓捕;翻滚目标模型简化为长方刚体,长方体一个顶点A代表抓捕点,其它的顶点B,C,D分别与点A形成翻滚目标长方体模型的三个相邻的棱边;假设空间机器人沿惯性系x轴正方向位于目标的前方,并沿惯性系x轴方向捕获翻滚目标,则确定捕获时机的第一个准则表述为:<img file="FDA0001008198130000032.GIF" wi="134" he="62" />且<img file="FDA0001008198130000033.GIF" wi="138" he="63" />且<img file="FDA0001008198130000034.GIF" wi="139" he="63" />其中,<img file="FDA0001008198130000035.GIF" wi="234" he="63" />分别表示翻滚目标上点A,B,C,D在惯性系下x方向的坐标;准则2:捕获时抓捕点位于空间机器人的路径无关工作空间,步骤二中用符号W<sub>PIW</sub>表示自由漂浮空间机器人的路径无关工作空间,如果将抓捕点到空间机器人系统质心的距离记为r<sub>dis</sub>,则确定抓捕时机的第二个准则表述为:r<sub>dis</sub>∈W<sub>PIW</sub>准则3:翻滚目标运动时,抓捕点的位置可能多次符合准则1和准则2,选择准则1和准则2第一次符合的时刻作为最佳抓捕时机,设定空间机器人的最小反应时间t<sub>react</sub>,并认为只有当确定的时机大于该时间时,空间机器人才可能完成对目标的捕获;因此,确定唯一的最佳抓捕时机:在t≥t<sub>react</sub>的前提下,符合准则1和准则2的最近时刻。
地址 710072 陕西省西安市碑林区友谊西路127号