CN115476356A - 一种无需力传感器的空间机械臂自适应混合阻抗控制方法 - Google Patents

一种无需力传感器的空间机械臂自适应混合阻抗控制方法 Download PDF

Info

Publication number
CN115476356A
CN115476356A CN202211015538.0A CN202211015538A CN115476356A CN 115476356 A CN115476356 A CN 115476356A CN 202211015538 A CN202211015538 A CN 202211015538A CN 115476356 A CN115476356 A CN 115476356A
Authority
CN
China
Prior art keywords
adaptive
matrix
force sensor
mechanical arm
impedance control
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.)
Pending
Application number
CN202211015538.0A
Other languages
English (en)
Inventor
李剑飞
张晓东
张沛
梁常春
赵志军
辛鹏飞
田�健
许哲
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Institute of Spacecraft System Engineering
Original Assignee
Beijing Institute of Spacecraft System Engineering
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Beijing Institute of Spacecraft System Engineering filed Critical Beijing Institute of Spacecraft System Engineering
Priority to CN202211015538.0A priority Critical patent/CN115476356A/zh
Publication of CN115476356A publication Critical patent/CN115476356A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • BPERFORMING OPERATIONS; TRANSPORTING
    • B25HAND TOOLS; PORTABLE POWER-DRIVEN TOOLS; MANIPULATORS
    • B25JMANIPULATORS; CHAMBERS PROVIDED WITH MANIPULATION DEVICES
    • B25J9/00Programme-controlled manipulators
    • B25J9/16Programme controls
    • B25J9/1656Programme controls characterised by programming, planning systems for manipulators
    • B25J9/1669Programme controls characterised by programming, planning systems for manipulators characterised by special application, e.g. multi-arm co-operation, assembly, grasping

Landscapes

  • Engineering & Computer Science (AREA)
  • Robotics (AREA)
  • Mechanical Engineering (AREA)
  • Feedback Control In General (AREA)
  • Manipulator (AREA)

Abstract

本发明提供了一种无需力传感器的空间机械臂自适应混合阻抗控制方法,包括:建立解耦递推自适应动力学控制器,对空间机械臂惯性参数进行自适应处理;建立无需力传感器的外力估计器,实现无需力传感器对机械臂任意位置受力进行估计;根据解耦递推自适应动力学控制器和无需力传感器的外力估计器,建立自适应混合阻抗控制,得到无需力传感器的自适应混合阻抗控制策略,在根据需求划分的任务子空间实现位置控制或阻抗控制。本发明的方法可显著提高空间机械臂执行接触任务时对自身动力学参数以及环境参数不确定性的自适应能力。

Description

一种无需力传感器的空间机械臂自适应混合阻抗控制方法
技术领域
本发明涉及一种空间机械臂自适应混合阻抗控制方法,属于空间机械臂控制技术领域。
背景技术
空间机械臂的高精度力控制一般需要末端六维力传感器或关节力传感器测量外部力的大小,但对于空间机械臂而言,末端六维力传感器或关节力传感器大大增加了系统的重量和成本,并且难以满足空间环境要求,以此提出一种无需力传感器空间机械臂自适应混合阻抗控制方法,基于复合自适应动力学自适应控制算法和高精度的外力估计,实现一种无需力传感器的复合自适应混合阻抗控制方法,对空间机器人的在轨应用具有重要的意义。
基于动力学模型的前馈设计是空间机器人控制算法的重要组成部分,可以大大提高空间机器人控制系统的精度、响应速度和稳定性,但在实际的系统中,通常很难精确地获得所有的惯性参数,并且各个多体系统间存在动力学特性耦合,对动力学建模的准确性和计算的实时性提出很大挑战。虽然自适应控制方案可以在参数不确定性较大的情况下保证机械臂的跟踪精度,但基于动力学模型的自适应控制算法参数多且相互耦合,难以实现参数的自适应化。
发明内容
本发明所要解决的技术问题是:克服现有技术的不足,本发明提供了一种无需力传感器的空间机械臂自适应混合阻抗控制方法,提出自适应混合阻抗控制策略,显著提高空间机械臂执行接触任务时对自身动力学参数以及环境参数不确定性的自适应能力。
为了达到上述目的,本发明的主要技术方案为:一种无需力传感器的空间机械臂自适应混合阻抗控制方法,包括:
建立解耦递推自适应动力学控制器,对空间机械臂惯性参数进行自适应处理;
建立无需力传感器的外力估计器,实现无需力传感器对机械臂任意位置受力进行估计;
根据解耦递推自适应动力学控制器和无需力传感器的外力估计器,建立自适应混合阻抗控制,得到无需力传感器的自适应混合阻抗控制策略,在根据需求划分的任务子空间实现位置控制或阻抗控制。
进一步的,所述建立解耦递推自适应动力学控制器,对空间机械臂惯性参数进行自适应处理,包括:
建立复合解耦控制律;
建立解耦复合自适应律,对不确定动力学参数进行处理;
采用所述复合解耦控制律和解耦复合自适应律,令机械臂在自由空间按照非规则连续轨迹运行,直到惯性参数收敛性满足要求。
进一步的,所述建立复合解耦控制律,包括:
根据如下控制律方程递推计算关节输出控制力矩:
Figure BDA0003811031340000021
Figure BDA0003811031340000022
Figure BDA0003811031340000023
Figure BDA0003811031340000024
Figure BDA0003811031340000025
τi=z6 TFri
其中,第i个连杆的广义速度
Figure BDA0003811031340000026
n为正整数;
vi为空间机械臂第i个连杆的线速度在其本体坐标系∑i中的表示,ωi为空间机械臂第i个连杆的线速度在其本体坐标系∑i中的表示;
iTi-1为坐标系∑i-1到坐标系∑i的广义变换矩阵:
Figure BDA0003811031340000027
iCi-1是旋转矩阵,
Figure BDA0003811031340000031
为坐标系∑i原点在坐标系∑i-1中的表示;
τi为关节i的输出控制力矩,
Figure BDA0003811031340000032
为关节角速度
Figure BDA0003811031340000033
的第i个分量,q为关节角度,
Figure BDA0003811031340000034
为关节参考速度
Figure BDA0003811031340000035
轨迹的第i个分量;z6为列向量,z6=[0,0,0,0,0,1]T;Vri为连杆i的参考速度;Vr0为Vri的初值;Fri为连杆i的参考力;
Mi是连杆i的惯性矩阵,Ci是连杆i的科里奥利和离心力矩阵,Gi是连杆i的重力矩阵;Yri是动力学方程的线性回归矩阵,θi是相应的空间机械臂惯性参数向量;
Figure BDA0003811031340000036
分别是Mi,Ci,Gi,θi估计值,KDi是正定对称矩阵;
Figure BDA0003811031340000037
Figure BDA0003811031340000038
Figure BDA0003811031340000039
其中,()×表示叉乘矩阵;mi为连杆i的质量;I3×3为单位矩阵;iIi为连杆i的转动惯量阵,
Figure BDA00038110313400000310
为连杆i的质心位置,iCci为本体坐标系∑i相对惯性系的旋转矩阵。
进一步的,所述建立解耦复合自适应律,对不确定动力学参数进行处理,包括:
利用机械臂的机械能构建与动力学参数的线性关系:
Ej=Uj(Vjj
其中,Ej为连杆j的机械能,Uj为机械臂的机械能与动力学参数的线性关系的回归矩阵;
机械臂的功率为:
Figure BDA00038110313400000311
其中,τ为以关节i的输出控制力矩τi作为元素构成的列向量;
利用一阶滤波器得出应用于复合自适应算法的动力学的线性关系式:
Figure BDA0003811031340000041
Figure BDA0003811031340000042
其中,λ为正常数,p为微分算子,Wj为机械臂的机械能与动力学参数的线性关系的关于惯性参数θj的回归矩阵,j=1,2,...,n;
每个关节的自适应律由如下公式求解:
Figure BDA0003811031340000043
其中,
Figure BDA0003811031340000044
滑模函数si=Vi-Vri,Pi为正定矩阵。
进一步的,所述采用所述复合解耦控制律和解耦复合自适应律,令机械臂在自由空间按照非规则连续轨迹运行,直到惯性参数收敛性满足要求,包括:
所述非规则连续轨迹的表达式为:
qd=qd(t)
Figure BDA0003811031340000045
其中,qd为关节角的期望轨迹,
Figure BDA0003811031340000046
为关节角速度的期望轨迹;
关节参考速度
Figure BDA0003811031340000047
其中,Λ为正定矩阵;
将关节参考速度表达式
Figure BDA0003811031340000048
联合所述复合解耦控制律和解耦复合自适应律,对机械臂惯性参数θi进行自适应控制。
进一步的,所述建立无需力传感器的外力估计器,实现无需力传感器对机械臂任意位置受力进行估计,包括:
机械臂动力学模型的拉格朗日的表示为:
Figure BDA0003811031340000049
其中,H为拉格朗日形式方程的惯性矩阵;C为拉格朗日形式方程的离心科氏力矩阵;G为拉格朗日形式方程的重力矩阵;f是关节摩擦力,τe是机械臂受到的外力等效到关节空间的值;
机械臂的动量为
Figure BDA00038110313400000511
则按照如下公式估计机械臂受到的外力等效到关节空间的值τe
Figure BDA0003811031340000051
其中,Ke为正定常矩阵,
Figure BDA0003811031340000052
为通过库仑摩擦模型拟合的摩擦力f的估计值,p(0)为设定的p的初值,t为时间;
Fe为末端笛卡尔空间机械臂受到的外力和外力矩的组合,如果无力传感器,则由关节外力估计值
Figure BDA0003811031340000053
通过下式计算得到估计值
Figure BDA0003811031340000054
Figure BDA0003811031340000055
J为机械臂的雅各比矩阵。
进一步的,所述自适应混合阻抗控制为双环架构,无需力传感器的自适应混合阻抗控制策略采用内环和外环控制。
进一步的,所述自适应混合阻抗控制策略的外环,根据期望的混合阻抗方程获得参考控制量,包括:
建立期望的混合阻抗方程:
Figure BDA0003811031340000056
其中,S为由0和1组成的对角矩阵,用于划分力控和位置环控制子空间;Md、Bd、Kd分别为期望的机械臂惯性矩阵、阻尼矩阵、刚度矩阵;X为机械臂末端位姿,Xd为期望末端位姿轨迹;
任务空间的参考速度为:
Figure BDA0003811031340000057
关节空间的参考速度
Figure BDA0003811031340000058
进一步的,所述自适应混合阻抗控制策略的内环,采用所述控制律方程的解耦递推自适应控制律,并将
Figure BDA0003811031340000059
代入
Figure BDA00038110313400000510
一种无需力传感器的空间机械臂自适应混合阻抗控制系统,包括:
第一模块,用于建立解耦递推自适应动力学控制器,对空间机械臂惯性参数进行自适应处理;
第二模块,用于建立无需力传感器的外力估计器,实现无需力传感器对机械臂任意位置受力进行估计;
第三模块,用于根据解耦递推自适应动力学控制器和无需力传感器的外力估计器,建立自适应混合阻抗控制,得到无需力传感器的自适应混合阻抗控制策略,在根据需求划分的任务子空间实现位置控制或阻抗控制。
本发明与现有技术相比的优点在于:
(1)本发明所提出的方法在空间机械臂存在不确定或因空间环境变化的动力学参数时可保证对不确定参数的适应以及控制的收敛性,并无需力传感器实现空间机械臂任意部位的外力感知,且实现空间机械臂自适应混合阻抗控制,实现对机械臂任务空间的位置和力的协调控制,显著提高空间机械臂执行接触任务时对自身动力学参数以及环境参数不确定性的自适应能力。
(2)本发明提出一种由跟踪误差和预测误差共同驱动参数自适应控制器,将复杂的多自由度动力学系统分解为多个相对简单的、具有相同特性的动力学子系统模块,针对各子系统模块实现控制参数的自适应化,从而提高控制系统的抗干扰能力和系统稳定性,并且在此基础上,根据系统动力学模型估计外部作用力,实现无需力传感器的自适应混合阻抗控制。
(3)本发明的方法通过将力位混合控制和阻抗控制的优点结合,较好地处理精度和稳定性的矛盾;基于解耦动力学模型设计控制算法,并结合自适应律处理建模的不确定性,得到高精度动力学模型,在此基础上,通过估计外部作用力,建立双环控制策略,该策略基于内环和外环控制,其中外环根据期望的混合阻抗方程获得参考控制量,内环采用上述递推复合自适应控制算法,使系统的状态量跟踪参考控制量,进而实现了期望的混合阻抗方程。
附图说明
为了更清楚地说明本发明实施例的技术方案,下面将对实施例中所需要使用的附图作简单的介绍,显而易见的,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性和劳动性的前提下,还可以根据这些附图获得其它附图。
图1是本发明的方法中坐标系关系图;
图2内环/外环控制策略结构图示意图;
图3是七自由度空间机械臂初始构型示意图;
图4是关节位置误差图;
图5是参数收敛情况图;
图6是机械臂末端在水平面上的轨迹(单位:m)图;
图7是末端位姿误差仿真图;
图8是z方向接触力仿真图;
图9是关节转矩仿真图。
具体实施方式
为了更好的理解本发明的技术方案,下面结合附图对本发明实施例进行详细描述。
应当明确,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例,都属于本发明保护的范围。
一种无需力传感器的空间机械臂自适应混合阻抗控制方法,包括步骤如下:
步骤1,建立解耦递推自适应动力学控制器,对空间机械臂惯性参数进行自适应处理。具体步骤为:
a)建立复合解耦控制律。
自由度为n的空间机械臂坐标系见图1,记空间机械臂第i个连杆的线速度和角速度在其本体坐标系∑i中分别表示为vi3×1和ωi3×1,定义第i个连杆的广义速度Vi为线速度与角速度的组合:
Figure BDA0003811031340000081
i=1,2,...,n;n为正整数;
定义iTi-1为坐标系∑i-1到坐标系∑i的广义变换矩阵:
Figure BDA0003811031340000082
其中,iCi-1是旋转矩阵,
Figure BDA0003811031340000083
为坐标系∑i原点在坐标系∑i-1中的表示。
根据以下控制律方程递推计算关节输出控制力矩:
Figure BDA0003811031340000084
Figure BDA0003811031340000085
Figure BDA0003811031340000086
Figure BDA0003811031340000087
τi=z6 TFri (7)
其中,τi为关节i的输出控制力矩,
Figure BDA0003811031340000088
为关节角速度
Figure BDA0003811031340000089
(关节角度记为q)的第i个分量,
Figure BDA00038110313400000810
为关节参考速度
Figure BDA00038110313400000811
轨迹的第i个分量,将在c)中和步骤3中进行定义;z6为列向量,z6=[0,0,0,0,0,1]T,Vri为连杆i的参考速度;Vr0为Vri的初值;Fri为连杆i的参考力;Mi,Ci,Gi分别是连杆i的惯性矩阵、科里奥利和离心力矩阵、重力矩阵,Yri是动力学方程的线性回归矩阵,θi是相应的惯性参数向量;
Figure BDA00038110313400000812
是Mi,Ci,Gi,θi估计值,KDi6×6是正定对称矩阵。
Mi,Ci,Gi的具体形式如下所示,式中iIi为连杆i的转动惯量阵,
Figure BDA00038110313400000813
为连杆i的质心位置,iCci为本体坐标系∑i相对惯性系的旋转矩阵。
Figure BDA00038110313400000814
Figure BDA0003811031340000091
Figure BDA0003811031340000092
其中,()×表示叉乘矩阵;mi为连杆i的质量;I3×3为单位矩阵;g为重力加速度;
b)建立解耦复合自适应律,对不确定动力学参数进行处理。
由于动力学参数未知,上述控制律使用的是惯性参数的估计值,为了建立递推的复合自适应律,利用机械臂的机械能构建与动力学参数的线性关系:
Figure BDA0003811031340000093
Ej为连杆j的机械能,Ej(0)为常数,Uj为机械臂的机械能与动力学参数的线性关系的回归矩阵。
机械能的导数即是机械臂的功率,也就是关节做功的功率:
Figure BDA0003811031340000094
τ为以关节i的输出控制力矩τi作为元素构成的列向量;
可见上述功率是可计算量,为了避免使用关节角加速度,利用一阶滤波器,得出可应用于复合自适应算法的动力学的线性关系:
Figure BDA0003811031340000095
λ为正常数,p为微分算子,Wj为机械臂的机械能与动力学参数的线性关系的关于惯性参数θj的回归矩阵,j=1,2,...,n,每个关节的自适应律可由下式求解:
Figure BDA0003811031340000101
其中,
Figure BDA0003811031340000102
滑模函数si=Vi-Vri,Pi为正定矩阵。
c)采用以上复合解耦控制律和解耦复合自适应律,令机械臂在自由空间按照非规则连续轨迹运行足够时间,直到惯性参数收敛性满足要求。
记非规则连续轨迹的表达式为:
Figure BDA0003811031340000103
qd为关节角的期望轨迹,
Figure BDA0003811031340000104
为关节角速度的期望轨迹;
则令式(5)中用到的
Figure BDA0003811031340000105
按如下公式计算
Figure BDA0003811031340000106
其中,Λ为正定矩阵,联合以上控制律和自适应律,得到完整的机械臂惯性参数θi辨识过程控制算法。以上自适应算法以跟踪误差和预测误差两类驱动信号为驱动源进行自适应控制,参数收敛速度和精度高于一般模型参考自适应算法,且计算量为O(n),低于一般的复合自适应算法1个阶次。
步骤2,建立无需力传感器的外力估计器,实现无需力传感器对机械臂任意位置受力进行估计。
针对无关节力传感器的机械臂关节,估计关节输出力矩和机械臂在空间受到的六维力。根据以上步骤得到准确的机械臂动力学模型,其拉格朗日的表示为:
Figure BDA0003811031340000107
H,C,G分别为拉格朗日形式方程的惯性矩阵、离心科氏力矩阵、重力矩阵。f是关节摩擦力,τe是机械臂受到的外力等效到关节空间的值。
机械臂的动量为
Figure BDA0003811031340000108
则按照如下公式估计机械臂受到的外力等效到关节空间的值τe
Figure BDA0003811031340000111
其中,Ke为正定常矩阵;
Figure BDA0003811031340000112
为通过库仑摩擦模型拟合的摩擦力f的估计值,可以用关节角速度来计算;输出力矩τ是控制算法计算的结果,为已知量;p(0)为设定的p的初值,t为时间。
自适应混合阻抗控制中需用到的Fe为末端笛卡尔空间机械臂受到的外力和外力矩的组合,如果无力传感器,则可由关节外力估计值
Figure BDA0003811031340000113
通过下式计算:
Figure BDA0003811031340000114
J为机械臂的雅各比矩阵;
步骤3,根据步骤1中的解耦递推自适应动力学控制器和步骤2中的无需力传感器的外力估计器,建立自适应混合阻抗控制,得到无需力传感器的自适应混合阻抗控制策略,在根据需求划分的任务子空间实现位置控制或阻抗控制。
建立自适应混合阻抗控制架构为双环架构,如图2所示,该策略基于内环和外环控制。
(1)外环根据期望的混合阻抗方程获得参考控制量。
首先建立期望的混合阻抗方程:
Figure BDA0003811031340000115
其中,S为由0和1组成的对角矩阵,用于划分力控和位置环控制子空间;Md、Bd、Kd分别为期望的机械臂惯性矩阵、阻尼矩阵、刚度矩阵,X为机械臂末端位姿,可由关节角根据运动学方程计算;Xd为期望末端位姿轨迹。
定义任务空间的参考速度为:
Figure BDA0003811031340000116
关节空间的参考速度为:
Figure BDA0003811031340000117
(2)内环采用公式(1)~(7)的解耦递推自适应控制律,并将公式(23)代入公式(5)。
以上共同构成无需力传感器的空间机械臂自适应混合阻抗方法。
一种无需力传感器的空间机械臂自适应混合阻抗控制系统,包括:
第一模块,用于建立解耦递推自适应动力学控制器,对空间机械臂惯性参数进行自适应处理;
第二模块,用于建立无需力传感器的外力估计器,实现无需力传感器对机械臂任意位置受力进行估计;
第三模块,用于根据解耦递推自适应动力学控制器和无需力传感器的外力估计器,建立自适应混合阻抗控制,得到无需力传感器的自适应混合阻抗控制策略,在根据需求划分的任务子空间实现位置控制或阻抗控制。
实施例:
以中国空间站大型机械臂为仿真模型,如图3,采用解耦递推自适应自适应控制算法对机械臂进行轨迹跟踪控制,得到的关节位置误差如图4,参数收敛情况如图5。从图5中可以看出,算法在轨迹跟踪误差的收敛精度和收敛时间都有优势,在参数也可以收敛到真值。
在此基础上,开展自适应混合阻抗控制仿真,机械臂在竖直方向为力控制,其他方向为位置控制。运动约束为保持末端姿态角不变,水平方向跟踪一个心形线,其轨迹方程为:
Figure BDA0003811031340000121
t为时间,单位为秒。竖直方向给定向下Fcz=10N的命令。
图6给出了机械臂末端在XY平面内的轨迹,可以看出在位控子空间很好地完成了轨迹跟踪任务。图8是z方向接触力的变化曲线,仿真初始时,机械臂末端未与地面接触,末端在偏置力命令的作用下以所设阻抗向下运动,经过约12秒与地面接触,有略微的冲击。图7给出了位控方向的跟踪误差,在保证竖直方向的接触力收敛的基础上,水平方向仍然具有很高的控制精度,位置控制误差小于0.01mm,角度控制误差优于3×10-5rad,最大误差出现在接触的一瞬间。z方向属于力控制子空间,图8中其偏差为与初始位置的偏差。图9是关节的输出转矩,转矩大小满足电机输出的要求,并且比较光滑,无高频震荡。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明保护的范围之内。
本发明说明书中未作详细描述的内容属本领域技术人员的公知技术。

Claims (10)

1.一种无需力传感器的空间机械臂自适应混合阻抗控制方法,其特征在于,包括:
建立解耦递推自适应动力学控制器,对空间机械臂惯性参数进行自适应处理;
建立无需力传感器的外力估计器,实现无需力传感器对机械臂任意位置受力进行估计;
根据解耦递推自适应动力学控制器和无需力传感器的外力估计器,建立自适应混合阻抗控制,得到无需力传感器的自适应混合阻抗控制策略,在根据需求划分的任务子空间实现位置控制或阻抗控制。
2.根据权利要求1所述的一种无需力传感器的空间机械臂自适应混合阻抗控制方法,其特征在于:所述建立解耦递推自适应动力学控制器,对空间机械臂惯性参数进行自适应处理,包括:
建立复合解耦控制律;
建立解耦复合自适应律,对不确定动力学参数进行处理;
采用所述复合解耦控制律和解耦复合自适应律,令机械臂在自由空间按照非规则连续轨迹运行,直到惯性参数收敛性满足要求。
3.根据权利要求2所述的一种无需力传感器的空间机械臂自适应混合阻抗控制方法,其特征在于:所述建立复合解耦控制律,包括:
根据如下控制律方程递推计算关节输出控制力矩:
Figure FDA0003811031330000011
Figure FDA0003811031330000012
Figure FDA0003811031330000013
Vr0=O6×1
Figure FDA0003811031330000014
Figure FDA0003811031330000015
τi=z6 TFri
其中,第i个连杆的广义速度
Figure FDA0003811031330000021
n为正整数;
vi为空间机械臂第i个连杆的线速度在其本体坐标系∑i中的表示,ωi为空间机械臂第i个连杆的线速度在其本体坐标系∑i中的表示;
iTi-1为坐标系∑i-1到坐标系∑i的广义变换矩阵:
Figure FDA0003811031330000022
iCi-1是旋转矩阵,
Figure FDA0003811031330000023
为坐标系∑i原点在坐标系∑i-1中的表示;
τi为关节i的输出控制力矩,
Figure FDA0003811031330000024
为关节角速度
Figure FDA0003811031330000025
的第i个分量,q为关节角度,
Figure FDA0003811031330000026
为关节参考速度
Figure FDA0003811031330000027
轨迹的第i个分量;z6为列向量,z6=[0,0,0,0,0,1]T;Vri为连杆i的参考速度;Vr0为Vri的初值;Fri为连杆i的参考力;
Mi是连杆i的惯性矩阵,Ci是连杆i的科里奥利和离心力矩阵,Gi是连杆i的重力矩阵;Yri是动力学方程的线性回归矩阵,θi是相应的空间机械臂惯性参数向量;
Figure FDA0003811031330000028
分别是Mi,Ci,Gi,θi估计值,KDi是正定对称矩阵;
Figure FDA0003811031330000029
Figure FDA00038110313300000210
Figure FDA00038110313300000211
其中,mi为连杆i的质量;I3×3为单位矩阵;iIi为连杆i的转动惯量阵,
Figure FDA00038110313300000212
为连杆i的质心位置,iCci为本体坐标系∑i相对惯性系的旋转矩阵。
4.根据权利要求3所述的一种无需力传感器的空间机械臂自适应混合阻抗控制方法,其特征在于:所述建立解耦复合自适应律,对不确定动力学参数进行处理,包括:
利用机械臂的机械能构建与动力学参数的线性关系:
Ej=Uj(Vjj
其中,Ej为连杆j的机械能,Uj为机械臂的机械能与动力学参数的线性关系的回归矩阵;
机械臂的功率为:
Figure FDA0003811031330000031
其中,τ为以关节i的输出控制力矩τi作为元素构成的列向量;
利用一阶滤波器得出应用于复合自适应算法的动力学的线性关系式:
Figure FDA0003811031330000032
Figure FDA0003811031330000033
其中,λ为正常数,p为微分算子,Wj为机械臂的机械能与动力学参数的线性关系的关于惯性参数θj的回归矩阵,j=1,2,...,n;
每个关节的自适应律由如下公式求解:
Figure FDA0003811031330000034
其中,
Figure FDA0003811031330000035
滑模函数si=Vi-Vri,Pi为正定矩阵。
5.根据权利要求4所述的一种无需力传感器的空间机械臂自适应混合阻抗控制方法,其特征在于:所述采用所述复合解耦控制律和解耦复合自适应律,令机械臂在自由空间按照非规则连续轨迹运行,直到惯性参数收敛性满足要求,包括:
所述非规则连续轨迹的表达式为:
qd=qd(t)
Figure FDA0003811031330000036
其中,qd为关节角的期望轨迹,
Figure FDA0003811031330000037
为关节角速度的期望轨迹;
关节参考速度
Figure FDA0003811031330000038
其中,Λ为正定矩阵;
将关节参考速度表达式
Figure FDA0003811031330000039
联合所述复合解耦控制律和解耦复合自适应律,对机械臂惯性参数θi进行自适应控制。
6.根据权利要求5所述的一种无需力传感器的空间机械臂自适应混合阻抗控制方法,其特征在于:所述建立无需力传感器的外力估计器,实现无需力传感器对机械臂任意位置受力进行估计,包括:
机械臂动力学模型的拉格朗日的表示为:
Figure FDA0003811031330000041
其中,H为拉格朗日形式方程的惯性矩阵;C为拉格朗日形式方程的离心科氏力矩阵;G为拉格朗日形式方程的重力矩阵;f是关节摩擦力,τe是机械臂受到的外力等效到关节空间的值;
机械臂的动量为
Figure FDA0003811031330000042
则按照如下公式估计机械臂受到的外力等效到关节空间的值τe
Figure FDA0003811031330000043
其中,Ke为正定常矩阵,
Figure FDA0003811031330000044
为通过库仑摩擦模型拟合的摩擦力f的估计值,p(0)为设定的p的初值,t为时间;
Fe为末端笛卡尔空间机械臂受到的外力和外力矩的组合,如果无力传感器,则由关节外力估计值
Figure FDA0003811031330000045
通过下式计算得到估计值
Figure FDA0003811031330000046
Figure FDA0003811031330000047
J为机械臂的雅各比矩阵。
7.根据权利要求6所述的一种无需力传感器的空间机械臂自适应混合阻抗控制方法,其特征在于:所述自适应混合阻抗控制为双环架构,无需力传感器的自适应混合阻抗控制策略采用内环和外环控制。
8.根据权利要求7所述的一种无需力传感器的空间机械臂自适应混合阻抗控制方法,其特征在于:所述自适应混合阻抗控制策略的外环,根据期望的混合阻抗方程获得参考控制量,包括:
建立期望的混合阻抗方程:
Figure FDA0003811031330000051
其中,S为由0和1组成的对角矩阵,用于划分力控和位置环控制子空间;Md、Bd、Kd分别为期望的机械臂惯性矩阵、阻尼矩阵、刚度矩阵;X为机械臂末端位姿,Xd为期望末端位姿轨迹;
任务空间的参考速度为:
Figure FDA0003811031330000052
关节空间的参考速度
Figure FDA0003811031330000053
9.根据权利要求8所述的一种无需力传感器的空间机械臂自适应混合阻抗控制方法,其特征在于:所述自适应混合阻抗控制策略的内环,采用所述控制律方程的解耦递推自适应控制律,并将
Figure FDA0003811031330000054
代入
Figure FDA0003811031330000055
Vr0=O6×1
Figure FDA0003811031330000056
10.一种无需力传感器的空间机械臂自适应混合阻抗控制系统,其特征在于,包括:
第一模块,用于建立解耦递推自适应动力学控制器,对空间机械臂惯性参数进行自适应处理;
第二模块,用于建立无需力传感器的外力估计器,实现无需力传感器对机械臂任意位置受力进行估计;
第三模块,用于根据解耦递推自适应动力学控制器和无需力传感器的外力估计器,建立自适应混合阻抗控制,得到无需力传感器的自适应混合阻抗控制策略,在根据需求划分的任务子空间实现位置控制或阻抗控制。
CN202211015538.0A 2022-08-23 2022-08-23 一种无需力传感器的空间机械臂自适应混合阻抗控制方法 Pending CN115476356A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211015538.0A CN115476356A (zh) 2022-08-23 2022-08-23 一种无需力传感器的空间机械臂自适应混合阻抗控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211015538.0A CN115476356A (zh) 2022-08-23 2022-08-23 一种无需力传感器的空间机械臂自适应混合阻抗控制方法

Publications (1)

Publication Number Publication Date
CN115476356A true CN115476356A (zh) 2022-12-16

Family

ID=84422301

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211015538.0A Pending CN115476356A (zh) 2022-08-23 2022-08-23 一种无需力传感器的空间机械臂自适应混合阻抗控制方法

Country Status (1)

Country Link
CN (1) CN115476356A (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117226825A (zh) * 2023-08-21 2023-12-15 北京航空航天大学杭州创新研究院 一种基于环境刚度在线辨识的机器人交互力估计方法
CN117733829A (zh) * 2024-01-11 2024-03-22 北方工业大学 多自由度工业机器人的关节力位协调自适应控制方法及系统

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117226825A (zh) * 2023-08-21 2023-12-15 北京航空航天大学杭州创新研究院 一种基于环境刚度在线辨识的机器人交互力估计方法
CN117226825B (zh) * 2023-08-21 2024-03-22 北京航空航天大学杭州创新研究院 一种基于环境刚度在线辨识的机器人交互力估计方法
CN117733829A (zh) * 2024-01-11 2024-03-22 北方工业大学 多自由度工业机器人的关节力位协调自适应控制方法及系统
CN117733829B (zh) * 2024-01-11 2024-08-27 北方工业大学 多自由度工业机器人的关节力位协调自适应控制方法及系统

Similar Documents

Publication Publication Date Title
CN106483964B (zh) 一种基于接触力观测器的机器人柔顺控制方法
CN115476356A (zh) 一种无需力传感器的空间机械臂自适应混合阻抗控制方法
Wen et al. Elman fuzzy adaptive control for obstacle avoidance of mobile robots using hybrid force/position incorporation
CN108621163A (zh) 一种面向汇榫工艺的冗余双臂机器人协作控制方法
Polverini et al. Sensorless and constraint based peg-in-hole task execution with a dual-arm robot
CN110427043B (zh) 基于作业飞行机器人重心偏移的位姿控制器设计方法
CN112276954B (zh) 基于有限时间输出状态受限的多关节机械臂阻抗控制方法
CN112558621A (zh) 一种基于解耦控制的飞行机械臂系统
CN108638068B (zh) 一种携带冗余度机械臂的飞行机器人控制系统设计方法
CN114942593B (zh) 一种基于干扰观测器补偿的机械臂自适应滑模控制方法
CN112809666B (zh) 一种基于神经网络的5-dof机械臂力位跟踪算法
CN109426147B (zh) 捕获卫星后组合航天器的自适应增益调整控制方法
CN111506095B (zh) 一种双刚体特征点间饱和固定时间相对位姿跟踪控制方法
CN109656132B (zh) 一种空间机器人有限时间协调控制方法
CN115625711A (zh) 一种考虑末端力的双臂机器人协同控制方法
CN116069050A (zh) 一种抗扰动的旋翼无人机吊载系统减摇控制方法及装置
CN114055467A (zh) 基于五自由度机器人的空间位姿在线仿真系统
CN117301064A (zh) 一种基于固定时间收敛的双臂机器人安全协同控制方法
WO2023165174A1 (zh) 构建机器人的控制器的方法、机器人的运动控制方法、装置以及机器人
CN116872197A (zh) 单杆机械臂的自适应神经网络反演控制方法及其控制系统
CN117055593A (zh) 一种旋翼无人机避障滑模容错控制方法
Anuar et al. Sliding mode control with adaptive gain of quadrotor with rigid manipulator
CN116141329A (zh) 一种基于自校正控制的空间闭环双臂机器人运动控制方法
CN113618743B (zh) 一种针对多源干扰的无人机机械臂末端位姿控制方法
Liu et al. On adaptive force/motion control of constrained robots

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