发明名称 一种空间绳系机器人目标抓捕后复合体姿态协调控制方法
摘要 本发明涉及一种空间绳系机器人目标抓捕后复合体姿态协调控制方法,引进了系绳拉力,设计了系绳拉力和推力器的协调姿态控制方法,首先设计了利用系绳拉力和推力器进行协调控制方法;然后考虑到利用系绳进行姿态控制的局限性,设计了仅仅利用推力器进行姿态控制协调控制方法,最后设计了控制器的切换控制方法,进行协调姿态控制。
申请公布号 CN103135552B 申请公布日期 2015.11.04
申请号 CN201310018221.7 申请日期 2013.01.18
申请人 西北工业大学 发明人 黄攀峰;王东科;孟中杰;刘正雄
分类号 G05D1/08(2006.01)I;G05B13/04(2006.01)I 主分类号 G05D1/08(2006.01)I
代理机构 西北工业大学专利中心 61204 代理人 王鲜凯
主权项 一种空间绳系机器人目标抓捕后复合体姿态协调控制方法,其特征在于步骤如下:步骤1:针对空间绳系机器人目标抓捕后复合体姿态协调控制,设计两种控制器:系绳和推力器协调控制器:系绳控制力矩为<img file="FDA0000706466320000011.GIF" wi="705" he="176" />对俯仰角和偏航角进行控制;其中:T<sub>l</sub>为系绳控制力矩,R为平台本体坐标系到目标本体坐标系之间的转换矩阵,d为抓捕点在复合体本体坐标系下位置矢量,×为叉乘算子;x、y和z为复合体质心三个位置分量;所述<img file="FDA0000706466320000012.GIF" wi="853" he="133" />其中:σ修正罗德里格参数,I<sub>3</sub>为3×3单位向量,σ<sup>×</sup>为矢量叉乘运算的反对称矩阵<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msup><mi>&sigma;</mi><mo>&times;</mo></msup><mo>=</mo><mfenced open='(' close=')'><mtable><mtr><mtd><mn>0</mn></mtd><mtd><mo>-</mo><msub><mi>&sigma;</mi><mn>3</mn></msub></mtd><mtd><msub><mi>&sigma;</mi><mn>2</mn></msub></mtd></mtr><mtr><mtd><msub><mi>&sigma;</mi><mn>3</mn></msub></mtd><mtd><mn>0</mn></mtd><mtd><mo>-</mo><msub><mi>&sigma;</mi><mn>1</mn></msub></mtd></mtr><mtr><mtd><mo>-</mo><msub><mi>&sigma;</mi><mn>2</mn></msub></mtd><mtd><msub><mi>&sigma;</mi><mn>1</mn></msub></mtd><mtd><mn>0</mn></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000706466320000013.GIF" wi="536" he="235" /></maths>σ<sup>T</sup>为σ的转置向量;所述x、y和z为复合体质心三个位置分量由复合体轨道动力学方程<maths num="0002" id="cmaths0002"><math><![CDATA[<mfenced open='{' close=''><mtable><mtr><mtd><mover><mi>x</mi><mrow><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mrow></mover><mo>-</mo><mn>2</mn><mi>n</mi><mover><mi>z</mi><mo>&CenterDot;</mo></mover><mo>=</mo><msub><mi>a</mi><mi>x</mi></msub><mo>+</mo><msub><mi>a</mi><mi>tx</mi></msub></mtd></mtr><mtr><mtd><mover><mi>y</mi><mrow><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mrow></mover><mo>+</mo><msup><mi>n</mi><mn>2</mn></msup><mi>y</mi><mo>=</mo><msub><mi>a</mi><mi>y</mi></msub><mo>+</mo><msub><mi>a</mi><mi>ty</mi></msub></mtd></mtr><mtr><mtd><mover><mi>z</mi><mrow><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mrow></mover><mo>+</mo><mn>2</mn><mi>n</mi><mover><mi>x</mi><mo>&CenterDot;</mo></mover><mo>-</mo><mn>3</mn><msup><mi>n</mi><mn>2</mn></msup><mi>z</mi><mo>=</mo><msub><mi>a</mi><mi>z</mi></msub><mo>+</mo><msub><mi>a</mi><mi>tz</mi></msub></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000706466320000014.GIF" wi="516" he="273" /></maths>得到,其中:a<sub>x</sub>、a<sub>y</sub>和a<sub>z</sub>为空间绳系机器人自身控制力产生的加速度,a<sub>tx</sub>、a<sub>ty</sub>和a<sub>tz</sub>为系绳拉力产生的加速度,n为轨道平均角速度;系绳拉力F<sub>l</sub>=(k<sub>pl</sub>+k<sub>dl</sub>s)η,其中:k<sub>pl</sub>、k<sub>dl</sub>分别为比例微分系数,s为复变量,η为d和系绳拉力矢量F之间的夹角;滚转角的推力器控制力矩T<sub>c</sub>(1)=k<sub>p</sub>φ<sub>e</sub>+k<sub>d</sub>ω<sub>ex</sub>,k<sub>p</sub>和k<sub>d</sub>分别为比例和微分参数,φ<sub>e</sub>为滚转角,ω<sub>ex</sub>为角速度误差ω<sub>e</sub>的第一个分量利用全推力器进行姿态控制的控制器:控制器为<img file="FDA0000706466320000015.GIF" wi="853" he="77" />其中:T<sub>c</sub>为空间绳系机器人自身姿态控制机构的控制力矩,ω为空间绳系机器人相对于惯性坐标系空间角速度在本体坐标系下的分量,ω由动力学/运动学方程<maths num="0003" id="cmaths0003"><math><![CDATA[<mfenced open='{' close=''><mtable><mtr><mtd><mi>J</mi><mover><mi>&omega;</mi><mo>&CenterDot;</mo></mover><mo>=</mo><mo>-</mo><msup><mi>&omega;</mi><mo>&times;</mo></msup><mi>J&omega;</mi><mo>+</mo><msub><mi>T</mi><mi>c</mi></msub><mo>+</mo><msub><mi>T</mi><mi>l</mi></msub></mtd></mtr><mtr><mtd><mover><mi>&sigma;</mi><mo>&CenterDot;</mo></mover><mo>=</mo><mi>G</mi><mrow><mo>(</mo><mi>&sigma;</mi><mo>)</mo></mrow><mi>&omega;</mi></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000706466320000021.GIF" wi="473" he="164" /></maths>得到,<img file="FDA0000706466320000022.GIF" wi="57" he="59" />和<img file="FDA0000706466320000023.GIF" wi="52" he="58" />分别是ω和修正罗德里格参数σ的对时间的导数,J为空间绳系机器人的正定对称的转动惯量矩阵,ω<sup>×</sup>为矢量叉乘运算的反对称矩阵;所述<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msup><mi>&omega;</mi><mo>&times;</mo></msup><mo>=</mo><mfenced open='(' close=')'><mtable><mtr><mtd><mn>0</mn></mtd><mtd><mo>-</mo><msub><mi>&omega;</mi><mn>3</mn></msub></mtd><mtd><msub><mi>&omega;</mi><mn>2</mn></msub></mtd></mtr><mtr><mtd><msub><mi>&omega;</mi><mn>3</mn></msub></mtd><mtd><mn>0</mn></mtd><mtd><mo>-</mo><msub><mi>&omega;</mi><mn>1</mn></msub></mtd></mtr><mtr><mtd><mo>-</mo><msub><mi>&omega;</mi><mn>2</mn></msub></mtd><mtd><msub><mi>&omega;</mi><mn>1</mn></msub></mtd><mtd><mn>0</mn></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000706466320000024.GIF" wi="530" he="239" /></maths>ω<sub>1</sub>、ω<sub>2</sub>和ω<sub>3</sub>为ω的三个分量,所述T<sub>c</sub>有不等式约束:|T<sub>c</sub>|≤|T<sub>cmax</sub>|,其中:T<sub>l</sub>为系绳产生的控制力矩,G(σ)定义如下:<img file="FDA0000706466320000025.GIF" wi="776" he="130" />I为3×3单位向量;ω<sub>d</sub>为本体系下期望角速度;姿态跟踪误差σ<sub>e</sub>和角速度误差ω<sub>e</sub>由误差动力学方程<maths num="0005" id="cmaths0005"><math><![CDATA[<mfenced open='{' close=''><mtable><mtr><mtd><mi>J</mi><msub><mover><mi>&omega;</mi><mo>&CenterDot;</mo></mover><mi>e</mi></msub><mo>+</mo><mi>J</mi><msub><mover><mi>&omega;</mi><mo>&CenterDot;</mo></mover><mi>d</mi></msub><mo>+</mo><msup><mi>&omega;</mi><mo>&times;</mo></msup><mi>J&omega;</mi><mo>=</mo><msub><mi>T</mi><mi>c</mi></msub></mtd></mtr><mtr><mtd><msub><mover><mi>&sigma;</mi><mo>&CenterDot;</mo></mover><mi>e</mi></msub><mo>=</mo><mi>G</mi><mrow><mo>(</mo><msub><mi>&sigma;</mi><mi>e</mi></msub><mo>)</mo></mrow><msub><mi>&omega;</mi><mi>e</mi></msub></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000706466320000026.GIF" wi="506" he="164" /></maths>得到,<img file="FDA0000706466320000027.GIF" wi="172" he="69" />和<img file="FDA0000706466320000028.GIF" wi="60" he="70" />分别为ω<sub>e</sub>、σ<sub>e</sub>和ω<sub>d</sub>对时间的导数;s=ω<sub>e</sub>+cσ<sub>e</sub>为滑动超平面,参数c>0;步骤2:当两种控制器的使用切换条件:<maths num="0006" id="cmaths0006"><math><![CDATA[<mfenced open='{' close=''><mtable><mtr><mtd><mo>|</mo><mi>&psi;</mi><mo>|</mo><mo>&le;</mo><msub><mi>&psi;</mi><mi>s</mi></msub></mtd></mtr><mtr><mtd><mo>|</mo><mi>&theta;</mi><mo>|</mo><mo>&le;</mo><msub><mi>&theta;</mi><mi>s</mi></msub></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000706466320000029.GIF" wi="195" he="174" /></maths>同时满足时,控制器从系绳和推力器协调控制器切换到仅利用全推力器进行姿态控制的控制器,ψ为俯仰角角度,θ为偏航角角度,ψ<sub>s</sub>、θ<sub>s</sub>为设置的控制器切换时的俯仰角和偏航角角度;实现了对抓捕后复合体的姿态控制。
地址 710072 陕西省西安市友谊西路127号