发明名称 一种空间多臂复杂连接联合体自适应控制方法
摘要 本发明公开了一种空间多臂复杂连接联合体自适应控制方法,本发明所提出的空间多臂复杂连接联合体的高精度强适应控制方法不依赖联合体的具体参数值,无需在轨对参数进行在线辨识,只需输出可测,控制器对参数不确定性具有强鲁棒性。本发明提出的控制方法,将简单自适应控制方法推广到相对阶数为高阶非近似严格正实性的被控对象。本发明处理机械臂闭链控制问题时,将机械臂系统闭链运动约束条件引入到轨迹规划中,通过轨迹规划给出各机械臂满足运动约束条件的理想轨迹,从而将机械臂闭链控制问题简化为各机械臂系统跟踪满足运动学约束的理想轨迹的跟踪问题。
申请公布号 CN103624784B 申请公布日期 2015.10.21
申请号 CN201310547914.5 申请日期 2013.11.06
申请人 北京控制工程研究所 发明人 王大轶;汤亮;何英姿;刘敏;魏春岭
分类号 B25J9/16(2006.01)I 主分类号 B25J9/16(2006.01)I
代理机构 中国航天科技专利中心 11009 代理人 臧春喜
主权项 一种空间多臂复杂连接联合体自适应控制方法,其特征在于步骤如下:(1)将空间多臂复杂连接联合体的自适应控制分解成联合体的姿态控制和机械臂的轨迹跟踪控制,联合体的姿态控制采用步骤(2)所示方法完成,机械臂的轨迹跟踪控制采用步骤(3)所示方法完成;(2)按照退步控制方法,将联合体的姿态控制模型按动力学方程与运动学方程进行退步简单自适应姿态控制器设计,具体实现过程如下:(2.1)首先根据联合体姿态控制要求设计参考动力学方程和参考运动学方程;参考动力学方程为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mfenced open='' close=''><mtable><mtr><mtd><msub><mover><mi>x</mi><mo>.</mo></mover><mi>m</mi></msub><mo>=</mo><msub><mi>A</mi><mi>m</mi></msub><msub><mi>x</mi><mi>m</mi></msub><mo>+</mo><msub><mi>B</mi><mi>m</mi></msub><msub><mi>u</mi><mi>m</mi></msub></mtd></mtr><mtr><mtd><msub><mi>y</mi><mi>m</mi></msub><mo>=</mo><msub><mi>C</mi><mi>m</mi></msub><msub><mi>x</mi><mi>m</mi></msub></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000703159880000011.GIF" wi="420" he="150" /></maths>其中状态变量x<sub>m</sub>=y<sub>m</sub>=ω<sub>d</sub>,<img file="FDA0000703159880000012.GIF" wi="167" he="72" />A<sub>m</sub>=0 B<sub>m</sub>=C<sub>m</sub>=I<sub>3</sub>,I<sub>3</sub>为3阶单位阵,ω<sub>d</sub>为期望姿态角速度;参考运动学方程为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mover><mi>q</mi><mo>.</mo></mover><mrow><mn>0</mn><mi>d</mi></mrow></msub></mtd></mtr><mtr><mtd><msub><mover><mi>q</mi><mo>.</mo></mover><mi>d</mi></msub></mtd></mtr></mtable></mfenced><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mfenced open='[' close=']'><mtable><mtr><mtd><mo>-</mo><msup><msub><mi>q</mi><mi>d</mi></msub><mi>T</mi></msup></mtd></mtr><mtr><mtd><msub><mi>q</mi><mrow><mn>0</mn><mi>d</mi></mrow></msub><msub><mi>I</mi><mn>3</mn></msub><mo>+</mo><msub><mover><mi>q</mi><mo>~</mo></mover><mi>d</mi></msub></mtd></mtr></mtable></mfenced><msub><mi>&omega;</mi><mi>d</mi></msub><mo>,</mo></mrow>]]></math><img file="FDA0000703159880000013.GIF" wi="564" he="174" /></maths>式中q<sub>0d</sub>、q<sub>d</sub>分别为期望四元数的标部和矢量部分,<img file="FDA0000703159880000014.GIF" wi="57" he="78" />为q<sub>d</sub>的反对称阵,(2.2)设计联合体姿态控制中间控制律ω<sub>m</sub>:ω<sub>m</sub>=A<sub>bd</sub>ω<sub>d</sub>‑kq<sub>e</sub>,A<sub>bd</sub>为卫星期望姿态角体坐标系到卫星实际体坐标系的转换矩阵,k为正定的控制系数矩阵;q<sub>e</sub>由联合体姿态四元数误差方程<maths num="0003" id="cmaths0003"><math><![CDATA[<mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>q</mi><mrow><mi>e</mi><mn>0</mn></mrow></msub><mo>=</mo><msub><mi>q</mi><mn>0</mn></msub><msub><mi>q</mi><mrow><mn>0</mn><mi>d</mi></mrow></msub><mo>+</mo><msup><mi>q</mi><mi>T</mi></msup><msub><mi>q</mi><mi>d</mi></msub></mtd></mtr><mtr><mtd><msub><mi>q</mi><mi>e</mi></msub><mo>=</mo><msub><mi>q</mi><mrow><mn>0</mn><mi>d</mi></mrow></msub><mi>q</mi><mo>-</mo><msub><mi>q</mi><mn>0</mn></msub><msub><mi>q</mi><mi>d</mi></msub><mo>+</mo><mover><mi>q</mi><mo>~</mo></mover><msub><mi>q</mi><mi>d</mi></msub></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000703159880000015.GIF" wi="479" he="179" /></maths>确定,式中q<sub>e0</sub>为误差四元数的标部,q<sub>0</sub>、q分别为实际四元数的标部和矢量部分,<img file="FDA0000703159880000016.GIF" wi="41" he="76" />为q的反对称阵;(2.3)设计联合体的退步简单自适应姿态控制律T:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>T</mi><mo>=</mo><mi>Kr</mi><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>K</mi><mi>x</mi></msub></mtd><mtd><msub><mi>K</mi><mi>u</mi></msub></mtd><mtd><msub><mi>K</mi><mi>e</mi></msub></mtd></mtr></mtable></mfenced><msup><mfenced open='[' close=']'><mtable><mtr><mtd><msup><msub><mi>&omega;</mi><mi>d</mi></msub><mi>T</mi></msup></mtd><mtd><msup><msub><mover><mi>&omega;</mi><mo>.</mo></mover><mi>d</mi></msub><mi>T</mi></msup></mtd><mtd><msup><mi>e</mi><mi>T</mi></msup></mtd></mtr></mtable></mfenced><mi>T</mi></msup><mo>,</mo></mrow>]]></math><img file="FDA0000703159880000017.GIF" wi="809" he="81" /></maths>式中K=[K<sub>x</sub> K<sub>u</sub> K<sub>e</sub>]为自适应控制系数矩阵,由参数自适应律确定,K<sub>x</sub> K<sub>u</sub> K<sub>e</sub>分别为组成自适应控制矩阵K的子矩阵;r为变量,<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mi>r</mi><mo>=</mo><msup><mfenced open='[' close=']'><mtable><mtr><mtd><msup><msub><mi>&omega;</mi><mi>d</mi></msub><mi>T</mi></msup></mtd><mtd><msup><msub><mover><mi>&omega;</mi><mo>.</mo></mover><mi>d</mi></msub><mi>T</mi></msup></mtd><mtd><msup><mi>e</mi><mi>T</mi></msup></mtd></mtr></mtable></mfenced><mi>T</mi></msup><mo>,</mo></mrow>]]></math><img file="FDA0000703159880000018.GIF" wi="423" he="79" /></maths>e为联合体动力学方程输出误差,e=ω‑ω<sub>m</sub>,ω为联合体姿态角速度输出,参数自适应律<img file="FDA0000703159880000019.GIF" wi="56" he="74" />为<img file="FDA00007031598800000110.GIF" wi="387" he="69" />K(0)=K<sub>0</sub>,K<sub>0</sub>为自适应控制参数矩阵的初始值,Γ为正定的自适应控制矩阵,σ为正定的控制系数;(2.4)利用执行机构按照步骤(2.3)设计的联合体退步简单自适应姿态控制律T对联合体进行姿态控制;(3)按照退步控制方法,将机械臂轨迹跟踪控制系统按动力学方程与运动学方程进行机械臂的退步简单自适应轨迹跟踪控制器设计,具体实现过程如下:(3.1)首先根据机械臂轨迹规划设计参考动力学方程和参考运动学方程;参考动力学方程为:<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><mfenced open='' close=''><mtable><mtr><mtd><msub><mover><mi>x</mi><mo>.</mo></mover><mrow><mn>2</mn><mi>m</mi></mrow></msub><mo>=</mo><msub><mi>A</mi><mrow><mn>2</mn><mi>m</mi></mrow></msub><msub><mi>x</mi><mrow><mn>2</mn><mi>m</mi></mrow></msub><mo>+</mo><msub><mi>B</mi><mrow><mn>2</mn><mi>m</mi></mrow></msub><msub><mi>u</mi><mrow><mn>2</mn><mi>m</mi></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>y</mi><mrow><mn>2</mn><mi>m</mi></mrow></msub><mo>=</mo><msub><mi>C</mi><mrow><mn>2</mn><mi>m</mi></mrow></msub><msub><mi>x</mi><mrow><mn>2</mn><mi>m</mi></mrow></msub></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000703159880000021.GIF" wi="493" he="148" /></maths>状态变量x<sub>2m</sub>=y<sub>2m</sub>=ω<sub>ar</sub>,<img file="FDA0000703159880000022.GIF" wi="201" he="76" />A<sub>2m</sub>=0,B<sub>2m</sub>=C<sub>2m</sub>=I<sub>7</sub>,I<sub>7</sub>为单位阵,ω<sub>ar</sub>为各机械臂期望的关节角速度;参考运动学方程为直接的积分关系:<img file="FDA0000703159880000023.GIF" wi="221" he="85" />θ<sub>ar</sub>为机械臂各关节角;(3.2)设计机械臂轨迹跟踪控制中间控制律ω<sub>arm</sub>:<img file="FDA0000703159880000024.GIF" wi="417" he="92" /><img file="FDA00007031598800000210.GIF" wi="73" he="90" />为运动学输出误差,<img file="FDA0000703159880000025.GIF" wi="332" he="80" />θ<sub>a</sub>为机械臂的实际关节角,k<sub>a</sub>为正定控制系数矩阵;(3.3)设计机械臂的退步简单自适应轨迹跟踪控制T<sub>a</sub>:<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><msub><mi>T</mi><mi>a</mi></msub><mo>=</mo><msub><mi>K</mi><mi>a</mi></msub><msub><mi>r</mi><mi>a</mi></msub><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>K</mi><mi>xa</mi></msub></mtd><mtd><msub><mi>K</mi><mi>ua</mi></msub></mtd><mtd><msub><mi>K</mi><mi>ea</mi></msub></mtd></mtr></mtable></mfenced><msup><mfenced open='[' close=']'><mtable><mtr><mtd><msup><msub><mover><mi>&theta;</mi><mo>.</mo></mover><mi>ar</mi></msub><mi>T</mi></msup></mtd><mtd><msup><msub><mover><mi>&theta;</mi><mrow><mo>.</mo><mo>.</mo></mrow></mover><mi>ar</mi></msub><mi>T</mi></msup></mtd><mtd><msup><msub><mi>e</mi><mn>2</mn></msub><mi>T</mi></msup></mtd></mtr></mtable></mfenced><mi>T</mi></msup><mo>,</mo></mrow>]]></math><img file="FDA0000703159880000026.GIF" wi="912" he="86" /></maths>式中K<sub>a</sub>=[K<sub>xa</sub> K<sub>ua</sub> K<sub>ea</sub>]为自适应控制系数矩阵,由自适应律确定,K<sub>xa</sub> K<sub>ua</sub> K<sub>ea</sub>为组成自适应控制矩阵K<sub>a</sub>的子矩阵,<maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><msub><mi>r</mi><mi>a</mi></msub><mo>=</mo><msup><mfenced open='[' close=']'><mtable><mtr><mtd><msup><msub><mover><mi>&theta;</mi><mo>&CenterDot;</mo></mover><mi>ar</mi></msub><mi>T</mi></msup></mtd><mtd><msup><msub><mover><mi>&theta;</mi><mrow><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mrow></mover><mi>ar</mi></msub><mi>T</mi></msup></mtd><mtd><msup><msub><mi>e</mi><mn>2</mn></msub><mi>T</mi></msup></mtd></mtr></mtable></mfenced><mi>T</mi></msup><mo>,</mo></mrow>]]></math><img file="FDA0000703159880000027.GIF" wi="460" he="86" /></maths>e<sub>2</sub>为机械臂系统动力学方程输出误差,e<sub>2</sub>=ω<sub>a</sub>‑ω<sub>arm</sub>,ω<sub>a</sub>为机械臂各关键角输出,自适应律<img file="FDA0000703159880000028.GIF" wi="68" he="82" />为:<img file="FDA0000703159880000029.GIF" wi="766" he="82" />K<sub>a0</sub>为自适应控制参数矩阵的初始值,Γ<sub>a</sub>为正定的自适应控制矩阵,σ<sub>2</sub>为正定的控制系数;(3.4)利用机械臂各关节上的执行机构按照步骤(3.3)设计的机械臂退步简单自适应轨迹跟踪控制律T<sub>a</sub>对机械臂进行轨迹跟踪控制。
地址 100080 北京市海淀区北京2729信箱