[go: up one dir, main page]

CN111687834B - Inverse Priority Impedance Control System and Control Method for Redundant Manipulator Arm of Mobile Manipulator - Google Patents

Inverse Priority Impedance Control System and Control Method for Redundant Manipulator Arm of Mobile Manipulator Download PDF

Info

Publication number
CN111687834B
CN111687834B CN202010369760.5A CN202010369760A CN111687834B CN 111687834 B CN111687834 B CN 111687834B CN 202010369760 A CN202010369760 A CN 202010369760A CN 111687834 B CN111687834 B CN 111687834B
Authority
CN
China
Prior art keywords
manipulator
control
task
inverse
priority
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
Application number
CN202010369760.5A
Other languages
Chinese (zh)
Other versions
CN111687834A (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.)
Shenzhen Wanzhida Technology Co ltd
Original Assignee
Guangxi University of Science and Technology
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 Guangxi University of Science and Technology filed Critical Guangxi University of Science and Technology
Priority to CN202010369760.5A priority Critical patent/CN111687834B/en
Publication of CN111687834A publication Critical patent/CN111687834A/en
Application granted granted Critical
Publication of CN111687834B publication Critical patent/CN111687834B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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/1628Programme controls characterised by the control loop
    • B25J9/1633Programme controls characterised by the control loop compliant, force, torque control, e.g. combined with position control
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B25HAND TOOLS; PORTABLE POWER-DRIVEN TOOLS; MANIPULATORS
    • B25JMANIPULATORS; CHAMBERS PROVIDED WITH MANIPULATION DEVICES
    • B25J9/00Programme-controlled manipulators
    • B25J9/16Programme controls
    • B25J9/1679Programme controls characterised by the tasks executed
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02PCLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P90/00Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
    • Y02P90/02Total factory control, e.g. smart factories, flexible manufacturing systems [FMS] or integrated manufacturing systems [IMS]

Landscapes

  • Engineering & Computer Science (AREA)
  • Robotics (AREA)
  • Mechanical Engineering (AREA)
  • Manipulator (AREA)

Abstract

The invention discloses a mobile manipulator redundant mechanical arm reverse priority impedance control system and a control method. The mechanical arm is convenient to move and comprises a movable moving platform, a mechanical arm and a control console for controlling the mechanical arm; the manipulator comprises a mechanical arm, a mounting seat, a vertical column, an output gripper and a vertical cylinder; the mechanical arm comprises a vertical lifting moving platform, an arm section I, an arm section II, an arm section III and an arm section IV; a vertical rail is arranged on the left surface of the vertical column, and the vertical lifting moving table is vertically arranged on the vertical rail in a sliding manner; the lower end of the vertical column is fixedly connected to the upper surface of the mounting seat, the cylinder seat of the vertical cylinder is fixedly connected to the upper surface of the mounting seat positioned at the left side of the vertical track, and the upper end of the telescopic rod of the vertical cylinder is fixedly connected to the lower surface of the vertical lifting moving table; the mount pad is fixed on moving platform.

Description

移动机械手冗余机械臂逆优先级阻抗控制系统及控制方法Reverse priority impedance control system and control method for redundant mechanical arms of mobile manipulator

技术领域Technical Field

本发明涉及机械手冗余机械臂逆优先级阻抗控制技术领域,尤其涉及移动机械手冗余机械臂逆优先级阻抗控制系统及控制方法。The present invention relates to the technical field of reverse priority impedance control of redundant manipulator arms of manipulators, and in particular to a reverse priority impedance control system and a control method for redundant manipulator arms of mobile manipulators.

背景技术Background Art

目前工业机器人采用的控制方法是把机械手上每一个关节都当作一个单独的伺服机构,即每个轴对应一个伺服器,每个伺服器通过总线控制,由控制器统一控制并协调工作;The control method currently used by industrial robots is to treat each joint on the manipulator as a separate servo mechanism, that is, each axis corresponds to a servo, and each servo is controlled by a bus and uniformly controlled and coordinated by a controller;

六个自由度的机械臂是具有完成空间定位能力最小自由度数的机械臂,多于六个自由度的机械臂统一称为冗余机械臂;A six-degree-of-freedom robot is a robot with the minimum number of degrees of freedom to complete spatial positioning. Robots with more than six degrees of freedom are collectively referred to as redundant robots.

现在机械手的机械臂阻抗控制方法在不同层次结构下不能实现期望的阻抗控制任务,因此,设计一种能在不同层次结构下可使机械手的冗余机械臂实现期望的阻抗控制任务的方法显得非常必要。The current impedance control method of the manipulator arm cannot achieve the desired impedance control task under different hierarchical structures. Therefore, it is very necessary to design a method that can enable the redundant manipulator arm of the manipulator to achieve the desired impedance control task under different hierarchical structures.

发明内容Summary of the invention

本发明是为了解决现有机械臂阻抗控制方法在不同层次结构下不能实现期望的阻抗控制任务的不足,提供一种一是能控制机械手的平衡,机械手便于移动;二是能在不同层次结构下可使机械手的冗余机械臂实现期望的阻抗控制任务的移动机械手冗余机械臂逆优先级阻抗控制系统及控制方法。The present invention aims to solve the deficiency that the existing robot arm impedance control method cannot achieve the desired impedance control task under different hierarchical structures, and provides a mobile robot arm redundant robot arm reverse priority impedance control system and control method which can control the balance of the robot arm to facilitate the movement of the robot arm and enable the redundant robot arms of the robot arm to achieve the desired impedance control task under different hierarchical structures.

以上技术问题是通过下列技术方案解决的:The above technical problems are solved by the following technical solutions:

移动机械手冗余机械臂逆优先级阻抗控制系统,包括机械手和控制机械手的控制台;还包括可移动的移动平台;所述机械手包括机械臂、安装座、竖直柱、输出抓手和竖直气缸;安装座固定在移动平台上;A mobile manipulator redundant manipulator reverse priority impedance control system includes a manipulator and a console for controlling the manipulator; and also includes a movable mobile platform; the manipulator includes a manipulator arm, a mounting seat, a vertical column, an output gripper and a vertical cylinder; the mounting seat is fixed on the mobile platform;

机械臂包括竖直升降移动台、臂段一、臂段二、臂段三和臂段四;The mechanical arm includes a vertical lifting and moving platform, arm segment one, arm segment two, arm segment three and arm segment four;

在竖直柱的左表面上设有竖直轨道,竖直升降移动台上下竖直滑动设置在竖直轨道上;竖直柱的下端固定连接在安装座的上表面上,竖直气缸的缸座固定连接在位于竖直轨道左方的安装座的上表面上,竖直气缸的伸缩杆竖直朝上布置,竖直气缸的伸缩杆上端固定连接在竖直升降移动台的下表面上;竖直升降移动台能在竖直气缸的伸缩杆的带动下可沿着竖直轨道上下运动,形成第一自由度;A vertical track is provided on the left surface of the vertical column, and the vertical lifting and moving platform is vertically slidably arranged on the vertical track; the lower end of the vertical column is fixedly connected to the upper surface of the mounting seat, the cylinder seat of the vertical cylinder is fixedly connected to the upper surface of the mounting seat located on the left side of the vertical track, the telescopic rod of the vertical cylinder is arranged vertically upward, and the upper end of the telescopic rod of the vertical cylinder is fixedly connected to the lower surface of the vertical lifting and moving platform; the vertical lifting and moving platform can move up and down along the vertical track driven by the telescopic rod of the vertical cylinder, forming a first degree of freedom;

臂段一包括A1段管和伸缩连接在A1段管的左管口内的A2段管,在A1段管内的右端固定设有伸缩杆水平朝左布置的一号气缸,一号气缸的伸缩杆固定连接在A2段管的右端;The arm section 1 includes an A1 section pipe and an A2 section pipe telescopically connected to the left pipe opening of the A1 section pipe. A No. 1 cylinder with a telescopic rod horizontally arranged to the left is fixedly provided at the right end of the A1 section pipe. The telescopic rod of the No. 1 cylinder is fixedly connected to the right end of the A2 section pipe.

臂段二包括B1段管和伸缩连接在B1段管的左管口内的B2段管,在B1段管内的右端固定设有伸缩杆水平朝左布置的二号气缸,二号气缸的伸缩杆固定连接在B2段管的右端;The arm section 2 includes a B1 section pipe and a B2 section pipe telescopically connected to the left pipe opening of the B1 section pipe. A No. 2 cylinder with a telescopic rod horizontally arranged to the left is fixedly provided at the right end of the B1 section pipe. The telescopic rod of the No. 2 cylinder is fixedly connected to the right end of the B2 section pipe.

在竖直升降移动台的左端设有由一台一号减速电机驱动的一号水平旋转轴,A1段管的右端固定连接在一号水平旋转轴上,使臂段一能水平转动形成第二自由度;在一号水平旋转轴上还设有能控制一号水平旋转轴转动的一号电磁刹车;A No. 1 horizontal rotating shaft driven by a No. 1 reduction motor is provided at the left end of the vertical lifting and moving platform, and the right end of the A1 section pipe is fixedly connected to the No. 1 horizontal rotating shaft, so that the arm section 1 can rotate horizontally to form a second degree of freedom; a No. 1 electromagnetic brake that can control the rotation of the No. 1 horizontal rotating shaft is also provided on the No. 1 horizontal rotating shaft;

在A2段管的左端设有由一台二号减速电机驱动的二号水平旋转轴,B1段管的右端固定连接在二号水平旋转轴上,使臂段二能水平转动形成第三自由度;在二号水平旋转轴上还设有能控制二号水平旋转轴转动的二号电磁刹车;A second horizontal rotation axis driven by a second reduction motor is provided at the left end of the A2 section tube, and the right end of the B1 section tube is fixedly connected to the second horizontal rotation axis, so that the arm section 2 can rotate horizontally to form a third degree of freedom; a second electromagnetic brake is also provided on the second horizontal rotation axis to control the rotation of the second horizontal rotation axis;

在B2段管的左端设有由一台三号减速电机驱动的三号水平旋转轴,臂段三的右端固定连接在三号水平旋转轴上,使臂段三能水平转动形成第四自由度;在三号水平旋转轴上还设有能控制三号水平旋转轴转动的三号电磁刹车;A No. 3 horizontal rotation axis driven by a No. 3 reduction motor is provided at the left end of the B2 section tube, and the right end of the arm section 3 is fixedly connected to the No. 3 horizontal rotation axis, so that the arm section 3 can rotate horizontally to form a fourth degree of freedom; a No. 3 electromagnetic brake is also provided on the No. 3 horizontal rotation axis to control the rotation of the No. 3 horizontal rotation axis;

在臂段三的左端设有由一台四号减速电机驱动的能在左右竖直面上转动的一号横向竖直旋转轴,臂段四的右端固定连接在一号横向竖直旋转轴上,使臂段四能在左右竖直面上竖直转动形成第五自由度;在一号横向竖直旋转轴上还设有能控制一号横向竖直旋转轴转动的四号电磁刹车;At the left end of arm segment 3, there is a No. 1 horizontal vertical rotation axis driven by a No. 4 reduction motor that can rotate on the left and right vertical planes. The right end of arm segment 4 is fixedly connected to the No. 1 horizontal vertical rotation axis, so that arm segment 4 can rotate vertically on the left and right vertical planes to form the fifth degree of freedom; on the No. 1 horizontal vertical rotation axis, there is also a No. 4 electromagnetic brake that can control the rotation of the No. 1 horizontal vertical rotation axis;

在臂段四的左端设有由一台五号减速电机驱动的能在前后竖直面上转动的一号纵向竖直旋转轴,输出抓手的右端固定连接在一号纵向竖直旋转轴上,使输出抓手的右端能在前后竖直面上竖直转动形成第六自由度;在一号纵向竖直旋转轴上还设有能控制一号纵向竖直旋转轴转动的五号电磁刹车;At the left end of arm section 4, there is a No. 1 longitudinal vertical rotation axis driven by a No. 5 reduction motor and capable of rotating on the front and rear vertical planes. The right end of the output gripper is fixedly connected to the No. 1 longitudinal vertical rotation axis, so that the right end of the output gripper can rotate vertically on the front and rear vertical planes to form the sixth degree of freedom; the No. 1 longitudinal vertical rotation axis is also provided with a No. 5 electromagnetic brake capable of controlling the rotation of the No. 1 longitudinal vertical rotation axis;

A2段管能在一号气缸的伸缩杆的带动下在A1段管内左右伸缩移动形成第七自由度;Driven by the telescopic rod of the No. 1 cylinder, the A2 section of pipe can move left and right in the A1 section of pipe to form the seventh degree of freedom;

B2段管能在二号气缸的伸缩杆的带动下在B1段管内左右伸缩移动形成第八自由度;Driven by the telescopic rod of the No. 2 cylinder, the B2 section of pipe can telescope left and right inside the B1 section of pipe to form the eighth degree of freedom;

一根一号水平管的左端水平固定连接在竖直柱的右表面上,在一号水平管内左右滑动设有平衡调节块,在一号水平管内的左端固定连接有伸缩杆水平朝右的平衡调节气缸,平衡调节气缸的伸缩杆的右端固定连接在平衡调节块上;The left end of a No. 1 horizontal tube is horizontally fixedly connected to the right surface of the vertical column. A balance adjustment block is provided in the No. 1 horizontal tube for sliding left and right. A balance adjustment cylinder with a telescopic rod facing horizontally to the right is fixedly connected to the left end of the No. 1 horizontal tube. The right end of the telescopic rod of the balance adjustment cylinder is fixedly connected to the balance adjustment block.

一号电磁刹车的控制端、二号电磁刹车的控制端、三号电磁刹车的控制端、四号电磁刹车的控制端、五号电磁刹车的控制端、一号减速电机的控制端、二号减速电机的控制端、三号减速电机的控制端、四号减速电机的控制端、五号减速电机的控制端、一号气缸的控制端、二号气缸的控制端、平衡调节气缸的控制端和竖直气缸的控制端分别控制连接在控制台上。The control end of electromagnetic brake No. 1, the control end of electromagnetic brake No. 2, the control end of electromagnetic brake No. 3, the control end of electromagnetic brake No. 4, the control end of electromagnetic brake No. 5, the control end of reduction motor No. 1, the control end of reduction motor No. 2, the control end of reduction motor No. 3, the control end of reduction motor No. 4, the control end of reduction motor No. 5, the control end of cylinder No. 1, the control end of cylinder No. 2, the control end of balance adjustment cylinder and the control end of vertical cylinder are respectively controlled and connected to the control console.

平衡调节气缸对平衡调节块进行左右移动控制就能控制竖直柱的平衡,进而控制机械手的平衡。便于机械手移动。由于本发明有的机械手有八个自由度,灵活性好,可靠性高,易于完成控制任务。The balance adjustment cylinder controls the balance adjustment block to move left and right to control the balance of the vertical column, and then controls the balance of the manipulator. It is convenient for the manipulator to move. Since the manipulator of the present invention has eight degrees of freedom, it has good flexibility, high reliability, and is easy to complete the control task.

一种移动机械手冗余机械臂逆优先级阻抗控制方法,包括如下步骤:A method for reverse priority impedance control of redundant mechanical arms of a mobile manipulator comprises the following steps:

步骤1、建立冗余机械臂运动学模型,并给出冗余机械臂零空间向量的梯度方向策略;Step 1: Establish a kinematic model of the redundant manipulator and provide a gradient direction strategy of the redundant manipulator null space vector;

步骤2、建立通过奇异鲁棒解得出消除奇异性算法的任务优先级解决策略;Step 2: Establish a task priority solution strategy for eliminating singularity through a singular robust solution;

步骤3、建立奇异鲁棒解逆运动学分析模型;Step 3, establishing a singular robust solution inverse kinematics analysis model;

步骤4、建立多任务冗余机械臂的逆优先控制策略;Step 4: Establish an inverse priority control strategy for the multi-task redundant manipulator;

步骤5、对具有主要任务和次要任务的冗余机械臂的反向控制方程进行简化处理;Step 5, simplifying the inverse control equation of the redundant manipulator with the primary task and the secondary task;

步骤6、建立机械手的逆优先力控制策略;Step 6: Establish the inverse priority force control strategy of the manipulator;

步骤7、采用关节速度来解决机械手逆优先阻抗控制中的外力与关节加速度之间的关系,从而得到机械手的逆优先阻抗控制保证;Step 7: Use joint velocity to solve the relationship between external force and joint acceleration in the inverse priority impedance control of the manipulator, so as to obtain the inverse priority impedance control guarantee of the manipulator;

步骤8、将位置控制空间的逆优先计算扩展到力控制空间的逆优先计算,从而获得机械手速度级逆优先阻抗控制的总体框架。Step 8: Extend the inverse priority calculation of the position control space to the inverse priority calculation of the force control space, thereby obtaining the overall framework of the manipulator speed-level inverse priority impedance control.

本方案冗余机械臂在关节空间中的运动是按照相反的顺序导出的;然后将笛卡尔阻抗控制与逆优先阻抗控制相结合,解决了逆分层阻抗控制问题,将笛卡尔阻抗控制行为分为高优先阻抗控制和低优先阻抗控制。其中高优先阻抗控制任务不会干扰低优先阻抗控制任务,关节空间中的运动按相反的顺序会受到影响,要在相应的投影算子中工作;最后实现了高优先阻抗控制任务,避免了低优先阻抗控制任务中可能出现的奇异性引起的变形。因此,所提出的逆优先阻抗控制方法可使冗余机械臂在适当的层次结构下实现期望的阻抗控制任务。In this scheme, the motion of the redundant manipulator in the joint space is derived in reverse order; then the Cartesian impedance control is combined with the inverse priority impedance control to solve the inverse hierarchical impedance control problem, and the Cartesian impedance control behavior is divided into high priority impedance control and low priority impedance control. Among them, the high priority impedance control task will not interfere with the low priority impedance control task, and the motion in the joint space will be affected in the opposite order, and it is necessary to work in the corresponding projection operator; finally, the high priority impedance control task is realized, avoiding the deformation caused by the singularity that may appear in the low priority impedance control task. Therefore, the proposed inverse priority impedance control method can enable the redundant manipulator to achieve the desired impedance control task under the appropriate hierarchy.

作为优选,建立冗余机械臂运动学模型,并给出冗余机械臂零空间向量的梯度方向策略实现过程如下:As a preferred method, a kinematic model of a redundant manipulator is established, and the gradient direction strategy implementation process of the redundant manipulator null space vector is given as follows:

定义末端执行器在笛卡尔空间的位姿、速度分别为x、

Figure GDA0004090627790000031
其关节空间的角位置、角速度分别为q、
Figure GDA0004090627790000032
J为n自由度机器人的雅可比矩阵,其中x∈Rn
Figure GDA0004090627790000033
Figure GDA0004090627790000034
J∈Rm×n;冗余自由度机械臂的正运动学方程可用下式描述:Define the position and velocity of the end effector in Cartesian space as x,
Figure GDA0004090627790000031
The angular position and angular velocity of its joint space are q,
Figure GDA0004090627790000032
J is the Jacobian matrix of the n-DOF robot, where x∈R n ,
Figure GDA0004090627790000033
Figure GDA0004090627790000034
J∈R m×n ; The forward kinematics equation of the redundant degree of freedom manipulator can be described by the following formula:

Figure GDA0004090627790000035
Figure GDA0004090627790000035

式(1)也被称为机械臂运动学速度模型;Formula (1) is also called the robot arm kinematic velocity model;

考虑到最小二乘法的解,最优问题可列为:Taking into account the solution of the least squares method, the optimal problem can be stated as:

Figure GDA0004090627790000036
Figure GDA0004090627790000036

式(1)的解可通过寻找最佳

Figure GDA0004090627790000037
来解决;The solution of formula (1) can be obtained by finding the optimal
Figure GDA0004090627790000037
to solve;

Figure GDA0004090627790000038
Figure GDA0004090627790000038

因此,式(1)的伪逆解可表示为:Therefore, the pseudo-inverse solution of formula (1) can be expressed as:

Figure GDA0004090627790000039
Figure GDA0004090627790000039

式中J+——雅可比矩阵的伪逆Where J + ——pseudo-inverse of Jacobian matrix

I——单位矩阵I——Identity Matrix

Figure GDA00040906277900000310
——任意零空间矢量
Figure GDA00040906277900000310
——Arbitrary null space vector

Figure GDA00040906277900000311
——最小范数解,定义了手的动作
Figure GDA00040906277900000311
——Minimum norm solution, defining the hand movements

Figure GDA00040906277900000312
——齐次解,在末端不产生动作
Figure GDA00040906277900000312
——Homogeneous solution, no action is produced at the end

式(4)代表了末端执行器的位置和姿态控制;在式(4)中加入任意残差,可得到包含零空间的一般表达式;利用上述方程可在零向量上实现多任务优化;Formula (4) represents the position and attitude control of the end effector. By adding any residual to Formula (4), a general expression including the null space can be obtained. The above equation can be used to realize multi-task optimization on the null vector.

然而,上述方程忽略了雅可比矩阵的病态;正则化方程可通过添加额外的正则化值来修改,However, the above equation ignores the ill-conditioning of the Jacobian matrix; the regularization equation can be modified by adding additional regularization values,

Figure GDA0004090627790000041
Figure GDA0004090627790000041

其中λ≥0是加权矩阵,

Figure GDA0004090627790000042
为加权系数,并且满足where λ≥0 is the weighting matrix,
Figure GDA0004090627790000042
is the weighting coefficient and satisfies

Figure GDA0004090627790000043
Figure GDA0004090627790000043

上述方程的解就可表示为:The solution to the above equation can be expressed as:

Figure GDA0004090627790000044
Figure GDA0004090627790000044

式(7)也被称为冗余机械臂运动学模型;Formula (7) is also called the redundant manipulator kinematic model;

冗余机械臂零空间向量的位置相关标量指数的关节限制梯度方向的关节限制函数为:The joint constraint function of the joint constraint gradient direction of the position-dependent scalar index of the redundant manipulator null space vector is:

Figure GDA0004090627790000045
Figure GDA0004090627790000045

作为优选,建立通过奇异鲁棒解得出消除奇异性算法的任务优先级解决策略如下:As a preferred method, the task priority solution strategy of eliminating singularity algorithm through singular robust solution is established as follows:

在雅可比矩阵的冗余机械臂求解中,优化任务是在主任务的零空间中实现的;反向任务运动学是在正向任务运动学基础上建立的:In the redundant manipulator solution of the Jacobian matrix, the optimization task is implemented in the null space of the main task; the inverse task kinematics is established based on the forward task kinematics:

Figure GDA0004090627790000046
Figure GDA0004090627790000046

其中

Figure GDA0004090627790000047
Figure GDA0004090627790000048
表示task1和task2in
Figure GDA0004090627790000047
and
Figure GDA0004090627790000048
Represents task1 and task2

从表达式(9)中得出冗余机械臂的逆运动学方程为:From expression (9), the inverse kinematics equation of the redundant manipulator is obtained as follows:

Figure GDA0004090627790000049
Figure GDA0004090627790000049

Task1作为主任务,task2作为辅助任务;也就是说,task2

Figure GDA00040906277900000418
是在task1
Figure GDA00040906277900000417
的零空间中实现的;冗余机械臂的最终逆运动学表达式如下所示:Task1 is the main task and task2 is the auxiliary task; that is, task2
Figure GDA00040906277900000418
It is in task1
Figure GDA00040906277900000417
The final inverse kinematics expression of the redundant manipulator is as follows:

Figure GDA00040906277900000410
Figure GDA00040906277900000410

其中

Figure GDA00040906277900000411
in
Figure GDA00040906277900000411

Figure GDA00040906277900000412
是投影矩阵,它给出了次要任务到主要任务的适用范围;
Figure GDA00040906277900000413
Figure GDA00040906277900000414
是所需的命令速度;
Figure GDA00040906277900000415
是主要任务,
Figure GDA00040906277900000416
是次要任务;
Figure GDA00040906277900000412
is the projection matrix, which gives the applicability of the secondary task to the primary task;
Figure GDA00040906277900000413
and
Figure GDA00040906277900000414
is the desired command speed;
Figure GDA00040906277900000415
is the main task,
Figure GDA00040906277900000416
It is a secondary task;

如果两个相关任务是相互依赖的,则相应的雅可比矩阵是奇异的;如果任务雅可比矩阵是奇异的,则相应的任务是不满足的;在这种情况下,雅可比相关矩阵将是奇点,定义为算法奇点;If two related tasks are interdependent, the corresponding Jacobian matrix is singular; if the task Jacobian matrix is singular, the corresponding tasks are unsatisfiable; in this case, the Jacobian correlation matrix will be a singular point, defined as an algorithmic singularity;

也就是说,如果That is, if

Figure GDA0004090627790000051
Figure GDA0004090627790000051

其中ρ(·)是矩阵的秩;where ρ(·) is the rank of the matrix;

很明显,算法的奇异性是由次要任务和主要任务之间的任务冲突引起的;此外,基于任务优先级的冗余机械臂逆运动学旨在提供更好的控制主要任务的有效性;It is obvious that the singularity of the algorithm is caused by the task conflict between the secondary and primary tasks; moreover, the redundant manipulator inverse kinematics based on task priority is designed to provide better effectiveness in controlling the primary task;

因此,让位置控制方向作为主要任务,从而使位置保证控制方向任务的准确性;然后建立通过奇异鲁棒解得出消除奇异性算法的任务优先级解决策略方程:Therefore, let the position control direction as the main task, so that the position guarantees the accuracy of the control direction task; then establish the task priority solution strategy equation of the algorithm to eliminate the singularity through the singular robust solution:

Figure GDA0004090627790000052
Figure GDA0004090627790000052

作为优选,建立奇异鲁棒解逆运动学分析模型如下:As a preferred method, a singular robust solution inverse kinematics analysis model is established as follows:

基于雅可比伪逆解会出现运动学奇异性,这是由次矩阵引起的;对于运动奇异性问题,应该还要给出DLS(阻尼最小二乘)解;Based on the Jacobi pseudo-inverse solution, kinematic singularities will occur, which are caused by the secondary matrix; for the kinematic singularity problem, the DLS (damped least squares) solution should also be given;

对DLS解的成本函数可修改为:The cost function of the DLS solution can be modified as follows:

Figure GDA0004090627790000053
Figure GDA0004090627790000053

因此,上述方程的奇异鲁棒伪逆解可表示为:Therefore, the singular robust pseudo-inverse solution of the above equation can be expressed as:

Figure GDA0004090627790000054
Figure GDA0004090627790000054

式(15)即为奇异鲁棒解逆运动学分析模型,设λ=ηΔI,上述DLS解等价于附加正则化解,标量值η平衡了任务精度和奇异性;Formula (15) is the inverse kinematics analysis model of the singular robust solution. Let λ = η Δ I. The above DLS solution is equivalent to the additional regularization solution. The scalar value η balances the task accuracy and singularity.

对于雅可比矩阵伪逆解的计算,可给出雅可比矩阵的奇异值SVD分解形式For the calculation of the pseudo-inverse solution of the Jacobian matrix, the singular value SVD decomposition form of the Jacobian matrix can be given

J=UΣVT (16)J=UΣV T (16)

其中U∈Rm×m,V∈Rn×n,∑∈Rm×n,U是由列向量ui组成的一元矩阵,V是由列向量vi组成的一元矩阵,∑是m×n对角矩阵的块矩阵,该对角矩阵包含J的奇异值σi≥0以递减顺序包含n-m个零列向量;where U∈R m×m , V∈R n×n , ∑∈R m×n , U is a one-element matrix consisting of column vectors u i , V is a one-element matrix consisting of column vectors vi , ∑ is a block matrix of an m×n diagonal matrix containing the singular values σ i ≥ 0 of J in decreasing order and containing nm zero column vectors;

Figure GDA0004090627790000055
Figure GDA0004090627790000055

其中,r≤m是矩阵J的秩;Where r≤m is the rank of the matrix J;

对于运动奇异性,参考计算伪逆解时需要的奇异值分解SVD,大的产生的关节速度是由于最小的奇异值迅速接近0,如下式:For motion singularity, refer to the singular value decomposition SVD required to calculate the pseudo-inverse solution. The large joint velocity is due to the small singular value rapidly approaching 0, as shown in the following formula:

Figure GDA0004090627790000056
Figure GDA0004090627790000056

因子λ0将影响奇异性,λ0值越高,阻尼越大,联合速度就越接近奇异点;此外,定义可变阻尼因子的策略也不同;我们可得到The factor λ 0 will affect the singularity. The higher the value of λ 0 , the greater the damping and the closer the joint velocity is to the singularity point. In addition, the strategy for defining the variable damping factor is also different; we can get

Figure GDA0004090627790000061
Figure GDA0004090627790000061

从上式中,我们可看出参数δ>0监视最小的奇异值。From the above formula, we can see that the parameter δ>0 monitors the smallest singular value.

作为优选,建立多任务冗余机械臂的逆优先控制策略如下;As a preferred method, the inverse priority control strategy of the multi-task redundant manipulator is established as follows;

引入逆优先投影矩阵

Figure GDA0004090627790000062
该矩阵包括了不依赖于第k个任务的最低优先级l-k-1任务的相应元素的零空间,所以得出Introducing the inverse priority projection matrix
Figure GDA0004090627790000062
This matrix includes the null space of the corresponding elements of the lowest priority lk-1 tasks that do not depend on the kth task, so we have

Figure GDA0004090627790000063
Figure GDA0004090627790000063

Figure GDA0004090627790000064
Figure GDA0004090627790000064

其中Ji|j是与j-th任务线性无关的i-th任务的所有组件相关联的雅可比矩阵;where Ji |j is the Jacobian matrix associated with all components of the i-th task that are linearly independent of the j-th task;

所以,优先级推导公式如下:Therefore, the priority derivation formula is as follows:

Figure GDA0004090627790000065
Figure GDA0004090627790000065

在上述推导中,k=1,l-1…,1,其中1为正整数;初始值

Figure GDA0004090627790000066
In the above derivation, k = 1, l-1 , ..., 1, where 1 is a positive integer; the initial value
Figure GDA0004090627790000066

为了给出线性无关雅可比矩阵J的一般计算形式,定义逆增广Jacobian矩阵为:In order to give the general calculation form of the linearly independent Jacobian matrix J, the inverse augmented Jacobian matrix is defined as:

Figure GDA0004090627790000067
Figure GDA0004090627790000067

就有了There you have it

Figure GDA0004090627790000068
Figure GDA0004090627790000068

其中

Figure GDA0004090627790000069
表示
Figure GDA00040906277900000610
的行;in
Figure GDA0004090627790000069
express
Figure GDA00040906277900000610
's line;

于是乎,

Figure GDA00040906277900000611
的伪逆解就可表示为:So,
Figure GDA00040906277900000611
The pseudo-inverse solution of can be expressed as:

Figure GDA00040906277900000612
Figure GDA00040906277900000612

Figure GDA00040906277900000613
and
Figure GDA00040906277900000613

其中,Tk表示矩阵

Figure GDA00040906277900000614
的扩充;Where T k represents the matrix
Figure GDA00040906277900000614
expansion;

最终的反向优先级投影可写为:The final reverse priority projection can be written as:

Figure GDA00040906277900000615
Figure GDA00040906277900000615

于是,我们就可得出伪逆解的表达式:Therefore, we can derive the expression of the pseudo-inverse solution:

Figure GDA0004090627790000071
Figure GDA0004090627790000071

建立多任务冗余机械臂的逆优先控制策略方程如下:The inverse priority control strategy equation for the multi-task redundant manipulator is established as follows:

Figure GDA0004090627790000072
Figure GDA0004090627790000072

作为优选,对具有主要任务和次要任务的冗余机械臂的反向控制方程进行简化处理如下:Preferably, the inverse control equation of the redundant manipulator with primary and secondary tasks is simplified as follows:

对于六自由度或七自由度冗余机械臂来说,没有足够的自由度来完成多个层次的任务;实行双任务优先控制是必要的;也就是说,机械手的运动控制是主要任务和次要任务;For a six- or seven-degree-of-freedom redundant manipulator, there are not enough degrees of freedom to complete multiple levels of tasks; it is necessary to implement dual-task priority control; that is, the motion control of the manipulator is the primary task and the secondary task;

具有主要任务和次要任务的冗余机械臂的反向控制方程如下The inverse control equation for a redundant manipulator with primary and secondary tasks is as follows

Figure GDA0004090627790000073
Figure GDA0004090627790000073

上面的公式与前面的表达式(11)有很大的不同,但是算法框架是相似的;在上面的方程中,

Figure GDA0004090627790000074
是次要任务,
Figure GDA0004090627790000075
是主任务;主任务在主任务的指定零空间中实现;反向优先级的核心点是投影矩阵
Figure GDA0004090627790000076
的计算;
Figure GDA0004090627790000077
的表达式如公式(30):The above formula is quite different from the previous expression (11), but the algorithm framework is similar; in the above equation,
Figure GDA0004090627790000074
It's a secondary task.
Figure GDA0004090627790000075
is the main task; the main task is implemented in the specified null space of the main task; the core point of the reverse priority is the projection matrix
Figure GDA0004090627790000076
Calculation of
Figure GDA0004090627790000077
The expression of is as follows:

Figure GDA0004090627790000078
Figure GDA0004090627790000078

利用前述公式(22)-(28)中的导子,可得到简化后具有主要任务和次要任务的冗余机械臂的反向控制方程:Using the derivations in the above formulas (22)-(28), the simplified inverse control equations of the redundant manipulator with primary and secondary tasks can be obtained:

Figure GDA0004090627790000079
Figure GDA0004090627790000079

作为优选,建立机械手的逆优先力控制策略如下:As a preferred method, the inverse priority force control strategy of the manipulator is established as follows:

力控制空间中机械手的动力学可写成:The dynamics of the manipulator in the force control space can be written as:

Figure GDA00040906277900000710
Figure GDA00040906277900000710

其中X是笛卡尔空间中的位置,M(X)是惯性矩阵,H(X,X)是非线性力,F是输入控制力,Fe是接触力;Where X is the position in Cartesian space, M(X) is the inertia matrix, H(X, X) is the nonlinear force, F is the input control force, and Fe is the contact force;

此外,基于雅可比矩阵的变换可得到输入关节力矩In addition, the input joint torque can be obtained based on the transformation of the Jacobian matrix

τ=JT(q)F (33)τ=J T (q)F (33)

操纵器在力控制空间中的期望运动方程可定义如下:The desired equation of motion of the manipulator in the force control space can be defined as follows:

Figure GDA00040906277900000711
Figure GDA00040906277900000711

其中Md和Bd是惯性和阻尼矩阵;Fd是指令力,Fe是接触力;Where Md and Bd are the inertia and damping matrices; Fd is the command force, and Fe is the contact force;

因此,环境和操纵器响应之间的关系可写成Therefore, the relationship between the environment and the manipulator response can be written as

Figure GDA00040906277900000712
Figure GDA00040906277900000712

上述两个方程的组合如下The combination of the above two equations is as follows

Figure GDA00040906277900000713
Figure GDA00040906277900000713

从上面的方程可看出,如果Me、Be和Ke已知,则Md和Bd的调整将影响系统响应;From the above equations, it can be seen that if Me , Be , and Ke are known, then adjustments to Md and Bd will affect the system response;

力控制使机械手能够与环境或人类相互作用;另外,在某些情况下,没有必要实现全方位的力控制,也没有必要保证全方位的力控制,也就是说,有时我们只是想保证某个方向的力跟踪控制精度;Force control enables the manipulator to interact with the environment or humans; in addition, in some cases, it is not necessary to achieve full-range force control, nor is it necessary to ensure full-range force control, that is, sometimes we just want to ensure the force tracking control accuracy in a certain direction;

因此有必要对机械手进行分级力控制;也就是说,有必要给出一个新的层次力控制框架;从上面的方程我们可得到期望的层次力控制关系如下Therefore, it is necessary to perform hierarchical force control on the manipulator; that is, it is necessary to give a new hierarchical force control framework; from the above equations we can get the expected hierarchical force control relationship as follows

Figure GDA0004090627790000081
Figure GDA0004090627790000081

Figure GDA0004090627790000082
Figure GDA0004090627790000082

所以这两个方程的积分公式可写成So the integral formula of these two equations can be written as

Figure GDA0004090627790000083
Figure GDA0004090627790000083

Figure GDA0004090627790000084
Figure GDA0004090627790000084

如果机械手末端执行器能够跟踪期望的笛卡尔速度为

Figure GDA0004090627790000085
Figure GDA0004090627790000086
则可实现机械手的精确力控制;笛卡尔速度与关节速度的关系应借鉴逆优先控制;因此,可得到机械手的逆优先力控制策略的方程:If the robot end effector can track the desired Cartesian velocity
Figure GDA0004090627790000085
and
Figure GDA0004090627790000086
Then the precise force control of the manipulator can be achieved; the relationship between Cartesian speed and joint speed should refer to the inverse priority control; therefore, the equation of the inverse priority force control strategy of the manipulator can be obtained:

Figure GDA0004090627790000087
Figure GDA0004090627790000087

上述方程所要求的关节速度将保证机械手的力控制;值得一提的是,上述力控制律只是速度级控制律,它依赖于内速度环控制;如果内位置控制效果良好,则可实现精确的力控制;由于内速度环控制可实现低频位置跟踪,所以外力环可实现低频力跟踪。The joint speed required by the above equation will ensure the force control of the manipulator; it is worth mentioning that the above force control law is only a speed-level control law, which relies on the inner speed loop control; if the inner position control effect is good, accurate force control can be achieved; since the inner speed loop control can achieve low-frequency position tracking, the outer force loop can achieve low-frequency force tracking.

作为优选,采用关节速度来解决机械手逆优先阻抗控制中的外力与关节加速度之间的关系,从而得到机械手的逆优先阻抗控制保证的实现方式如下:As a preferred method, the joint velocity is used to solve the relationship between the external force and the joint acceleration in the reverse priority impedance control of the manipulator, so that the reverse priority impedance control of the manipulator is guaranteed to be implemented as follows:

当机械手实施力控制时,机械手在一定程度上起到了发起者的作用,也就是说,机械手已经做好了响应外部环境的准备;当机械臂作为阻抗控制模型工作时,机械臂会被动地响应外力;When the manipulator implements force control, the manipulator plays the role of an initiator to a certain extent, that is, the manipulator is ready to respond to the external environment; when the manipulator works as an impedance control model, the manipulator responds to the external force passively;

外力与关节加速度的对应阻抗关系可表示为The corresponding impedance relationship between external force and joint acceleration can be expressed as

Figure GDA0004090627790000088
Figure GDA0004090627790000088

Figure GDA0004090627790000089
Figure GDA0004090627790000089

参考速度可表示为The reference speed can be expressed as

Figure GDA00040906277900000810
Figure GDA00040906277900000810

Figure GDA0004090627790000091
Figure GDA0004090627790000091

因此,机械手的逆优先阻抗控制保证的表达式为:Therefore, the expression guaranteed by the reverse priority impedance control of the manipulator is:

Figure GDA0004090627790000092
Figure GDA0004090627790000092

作为优选,将位置控制空间的逆优先计算扩展到力控制空间的逆优先计算,从而获得机械手速度级逆优先阻抗控制的总体框架实现方式如下:As a preferred method, the inverse priority calculation of the position control space is extended to the inverse priority calculation of the force control space, so as to obtain the overall framework of the manipulator speed level inverse priority impedance control as follows:

混合阻抗应用就是上述两种策略的结合,即笛卡尔任务可分为两种情况:第一种是位置控制子空间,阻抗控制是在该子空间中实现的;第二个是力控制子空间,力控制在该子空间中实现;Hybrid impedance application is a combination of the above two strategies, i.e., the Cartesian task can be divided into two cases: the first is the position control subspace, in which impedance control is implemented; the second is the force control subspace, in which force control is implemented;

因此选择一个选择矩阵;外力与位置响应的关系如下Therefore, a selection matrix is selected; the relationship between external force and position response is as follows

Figure GDA0004090627790000093
Figure GDA0004090627790000093

Figure GDA0004090627790000094
Figure GDA0004090627790000094

所以期望速度的简化形式可表示为So the simplified form of the expected speed can be expressed as

Figure GDA0004090627790000095
Figure GDA0004090627790000095

Figure GDA0004090627790000096
Figure GDA0004090627790000096

然后我们得到了基于反向优先级的解决方案Then we get the solution based on reverse priority

Figure GDA0004090627790000097
Figure GDA0004090627790000097

考虑到n层任务,相应的阻抗控制任务也属于n层框架,因此,机械手速度级逆优先阻抗控制的总体框架表达式如下Considering the n-layer task, the corresponding impedance control task also belongs to the n-layer framework. Therefore, the overall framework expression of the manipulator speed level reverse priority impedance control is as follows:

Figure GDA0004090627790000098
Figure GDA0004090627790000098

表达式(52)解决了将位置控制空间的逆优先计算扩展到力控制空间的逆优先计算的机械手逆优先混合阻抗控制中去,能在不同层次结构下可使机械手的冗余机械臂实现期望的阻抗控制任务。Expression (52) solves the problem of extending the inverse priority calculation of the position control space to the inverse priority calculation of the force control space in the manipulator inverse priority mixed impedance control, which can enable the redundant manipulator arms of the manipulator to achieve the desired impedance control task under different hierarchical structures.

本发明能够达到如下效果:The present invention can achieve the following effects:

本发明能控制机械手的平衡,机械手便于移动,能在不同层次结构下可使机械手的冗余机械臂实现期望的阻抗控制任务。The present invention can control the balance of the manipulator, the manipulator is easy to move, and can enable the redundant manipulator arms of the manipulator to achieve the desired impedance control task under different hierarchical structures.

附图说明BRIEF DESCRIPTION OF THE DRAWINGS

图1为本发明力控制的动力学方案示意图。FIG. 1 is a schematic diagram of the kinetic scheme of force control of the present invention.

图2为本发明阻抗控制的动力学方案示意图。FIG. 2 is a schematic diagram of the kinetic scheme of impedance control of the present invention.

图3为本发明混合阻抗控制的动力学方案示意图。FIG3 is a schematic diagram of a kinetic scheme of hybrid impedance control according to the present invention.

图4为本发明实施例的一种七自由度机械手连接结构示意图。FIG. 4 is a schematic diagram of a connection structure of a seven-degree-of-freedom manipulator according to an embodiment of the present invention.

图5为本发明实施例的一种电路原理连接结构示意框图。FIG. 5 is a schematic block diagram of a circuit principle connection structure according to an embodiment of the present invention.

具体实施方式DETAILED DESCRIPTION

下面结合附图与实施例对本发明作进一步的说明。The present invention will be further described below in conjunction with the accompanying drawings and embodiments.

实施例,移动机械手冗余机械臂逆优先级阻抗控制系统,参见图4、图5所示,包括机械手和控制机械手的控制台S31;还包括可移动的移动平台S41;所述机械手包括机械臂、安装座S1、竖直柱S2、输出抓手S12和竖直气缸S23;安装座固定在移动平台上;Embodiment, mobile manipulator redundant manipulator reverse priority impedance control system, as shown in FIG4 and FIG5, includes a manipulator and a control console S31 for controlling the manipulator; also includes a movable mobile platform S41; the manipulator includes a manipulator arm, a mounting seat S1, a vertical column S2, an output gripper S12 and a vertical cylinder S23; the mounting seat is fixed on the mobile platform;

机械臂包括竖直升降移动台S3、臂段一S6、臂段二S7、臂段三S8和臂段四S10;The mechanical arm includes a vertical lifting and moving platform S3, arm segment one S6, arm segment two S7, arm segment three S8 and arm segment four S10;

在竖直柱的左表面上设有竖直轨道S24,竖直升降移动台上下竖直滑动设置在竖直轨道上;竖直柱的下端固定连接在安装座的上表面上,竖直气缸的缸座S21固定连接在位于竖直轨道左方的安装座的上表面上,竖直气缸的伸缩杆S22竖直朝上布置,竖直气缸的伸缩杆上端固定连接在竖直升降移动台的下表面上;竖直升降移动台能在竖直气缸的伸缩杆的带动下可沿着竖直轨道上下运动,形成第一自由度;A vertical track S24 is provided on the left surface of the vertical column, and the vertical lifting and moving platform is vertically slidably arranged up and down on the vertical track; the lower end of the vertical column is fixedly connected to the upper surface of the mounting seat, the cylinder seat S21 of the vertical cylinder is fixedly connected to the upper surface of the mounting seat located on the left side of the vertical track, the telescopic rod S22 of the vertical cylinder is arranged vertically upward, and the upper end of the telescopic rod of the vertical cylinder is fixedly connected to the lower surface of the vertical lifting and moving platform; the vertical lifting and moving platform can move up and down along the vertical track driven by the telescopic rod of the vertical cylinder, forming a first degree of freedom;

臂段一包括A1段管S13和伸缩连接在A1段管的左管口内的A2段管S14,在A1段管内的右端固定设有伸缩杆水平朝左布置的一号气缸S25,一号气缸的伸缩杆固定连接在A2段管的右端;The arm section 1 includes an A1 section pipe S13 and an A2 section pipe S14 which is telescopically connected to the left pipe opening of the A1 section pipe. A No. 1 cylinder S25 with a telescopic rod horizontally arranged to the left is fixedly provided at the right end of the A1 section pipe. The telescopic rod of the No. 1 cylinder is fixedly connected to the right end of the A2 section pipe.

臂段二包括B1段管S16和伸缩连接在B1段管的左管口内的B2段管S17,在B1段管内的右端固定设有伸缩杆水平朝左布置的二号气缸37,二号气缸的伸缩杆固定连接在B2段管的右端;The arm section 2 includes a B1 section pipe S16 and a B2 section pipe S17 telescopically connected to the left pipe opening of the B1 section pipe. A second cylinder 37 with a telescopic rod horizontally arranged to the left is fixedly provided at the right end of the B1 section pipe. The telescopic rod of the second cylinder is fixedly connected to the right end of the B2 section pipe.

在竖直升降移动台的左端设有由一台一号减速电机S26驱动的一号水平旋转轴S4,A1段管的右端固定连接在一号水平旋转轴上,使臂段一能水平转动形成第二自由度;在一号水平旋转轴上还设有能控制一号水平旋转轴转动的一号电磁刹车S32;A No. 1 horizontal rotation axis S4 driven by a No. 1 reduction motor S26 is provided at the left end of the vertical lifting and moving platform. The right end of the A1 section pipe is fixedly connected to the No. 1 horizontal rotation axis, so that the arm section 1 can rotate horizontally to form a second degree of freedom; a No. 1 electromagnetic brake S32 is also provided on the No. 1 horizontal rotation axis to control the rotation of the No. 1 horizontal rotation axis;

在A2段管的左端设有由一台二号减速电机S27驱动的二号水平旋转轴S15,B1段管的右端固定连接在二号水平旋转轴上,使臂段二能水平转动形成第三自由度;在二号水平旋转轴上还设有能控制二号水平旋转轴转动的二号电磁刹车S33;A second horizontal rotation axis S15 driven by a second reduction motor S27 is provided at the left end of the A2 section pipe, and the right end of the B1 section pipe is fixedly connected to the second horizontal rotation axis, so that the arm section 2 can rotate horizontally to form a third degree of freedom; a second electromagnetic brake S33 capable of controlling the rotation of the second horizontal rotation axis is also provided on the second horizontal rotation axis;

在B2段管的左端设有由一台三号减速电机S28驱动的三号水平旋转轴S18,臂段三的右端固定连接在三号水平旋转轴上,使臂段三能水平转动形成第四自由度;在三号水平旋转轴上还设有能控制三号水平旋转轴转动的三号电磁刹车S34;A No. 3 horizontal rotation axis S18 driven by a No. 3 reduction motor S28 is provided at the left end of the B2 section tube, and the right end of the arm section 3 is fixedly connected to the No. 3 horizontal rotation axis, so that the arm section 3 can rotate horizontally to form a fourth degree of freedom; a No. 3 electromagnetic brake S34 is also provided on the No. 3 horizontal rotation axis to control the rotation of the No. 3 horizontal rotation axis;

在臂段三的左端设有由一台四号减速电机S29驱动的能在左右竖直面上转动的一号横向竖直旋转轴S9,臂段四的右端固定连接在一号横向竖直旋转轴上,使臂段四能在左右竖直面上竖直转动形成第五自由度;在一号横向竖直旋转轴上还设有能控制一号横向竖直旋转轴转动的四号电磁刹车S35;At the left end of arm segment 3, there is a No. 1 horizontal vertical rotation axis S9 driven by a No. 4 reduction motor S29 and capable of rotating on the left and right vertical planes. The right end of arm segment 4 is fixedly connected to the No. 1 horizontal vertical rotation axis, so that arm segment 4 can rotate vertically on the left and right vertical planes to form the fifth degree of freedom; on the No. 1 horizontal vertical rotation axis, there is also a No. 4 electromagnetic brake S35 capable of controlling the rotation of the No. 1 horizontal vertical rotation axis;

在臂段四的左端设有由一台五号减速电机S30驱动的能在前后竖直面上转动的一号纵向竖直旋转轴S11,输出抓手的右端固定连接在一号纵向竖直旋转轴上,使输出抓手的右端能在前后竖直面上竖直转动形成第六自由度;在一号纵向竖直旋转轴上还设有能控制一号纵向竖直旋转轴转动的五号电磁刹车S36;At the left end of arm section 4, there is a No. 1 longitudinal vertical rotation axis S11 driven by a No. 5 reduction motor S30 and capable of rotating on the front and rear vertical planes. The right end of the output gripper is fixedly connected to the No. 1 longitudinal vertical rotation axis, so that the right end of the output gripper can rotate vertically on the front and rear vertical planes to form the sixth degree of freedom; the No. 1 longitudinal vertical rotation axis is also provided with a No. 5 electromagnetic brake S36 capable of controlling the rotation of the No. 1 longitudinal vertical rotation axis;

A2段管能在一号气缸的伸缩杆的带动下在A1段管内左右伸缩移动形成第七自由度;Driven by the telescopic rod of the No. 1 cylinder, the A2 section of pipe can move left and right in the A1 section of pipe to form the seventh degree of freedom;

B2段管能在二号气缸的伸缩杆的带动下在B1段管内左右伸缩移动形成第八自由度;Driven by the telescopic rod of the No. 2 cylinder, the B2 section of pipe can telescope left and right inside the B1 section of pipe to form the eighth degree of freedom;

一根一号水平管S39的左端水平固定连接在竖直柱的右表面上,在一号水平管内左右滑动设有平衡调节块S40,在一号水平管内的左端固定连接有伸缩杆水平朝右的平衡调节气缸S38,平衡调节气缸的伸缩杆的右端固定连接在平衡调节块上;The left end of a No. 1 horizontal tube S39 is horizontally fixedly connected to the right surface of the vertical column. A balance adjustment block S40 is provided in the No. 1 horizontal tube for sliding left and right. A balance adjustment cylinder S38 with a telescopic rod facing horizontally to the right is fixedly connected to the left end of the No. 1 horizontal tube. The right end of the telescopic rod of the balance adjustment cylinder is fixedly connected to the balance adjustment block.

一号电磁刹车的控制端、二号电磁刹车的控制端、三号电磁刹车的控制端、四号电磁刹车的控制端、五号电磁刹车的控制端、一号减速电机的控制端、二号减速电机的控制端、三号减速电机的控制端、四号减速电机的控制端、五号减速电机的控制端、一号气缸的控制端、二号气缸的控制端、平衡调节气缸的控制端和竖直气缸的控制端分别控制连接在控制台上。第七自由度能让机械臂自由伸长,第八自由度能让机械臂自由伸长,大大增加了作业范围和灵活性。The control end of electromagnetic brake No. 1, electromagnetic brake No. 2, electromagnetic brake No. 3, electromagnetic brake No. 4, electromagnetic brake No. 5, reduction motor No. 1, reduction motor No. 2, reduction motor No. 3, reduction motor No. 4, reduction motor No. 5, cylinder No. 1, cylinder No. 2, balance cylinder and vertical cylinder are respectively connected to the control console. The seventh degree of freedom allows the robot arm to extend freely, and the eighth degree of freedom allows the robot arm to extend freely, greatly increasing the range of operation and flexibility.

平衡调节气缸对平衡调节块进行左右移动控制就能控制竖直柱的平衡,进而控制机械手的平衡。机械手便于移动。移动平台包括汽车。The balance adjustment cylinder controls the balance adjustment block to move left and right to control the balance of the vertical column, thereby controlling the balance of the manipulator. The manipulator is easy to move. The mobile platform includes a car.

由于本发明有的机械手有八个自由度,灵活性好,可靠性高,易于完成控制任务。Since the manipulator of the present invention has eight degrees of freedom, it has good flexibility, high reliability, and is easy to complete control tasks.

移动机械手冗余机械臂逆优先级阻抗控制方法,参见图1-3所示。包括如下步骤:The reverse priority impedance control method of the redundant manipulator arm of the mobile manipulator is shown in Figure 1-3. It includes the following steps:

步骤1、建立冗余机械臂运动学模型,并给出冗余机械臂零空间向量的梯度方向策略实现过程如下:Step 1: Establish a kinematic model of the redundant manipulator and give the gradient direction strategy of the redundant manipulator null space vector. The implementation process is as follows:

定义末端执行器在笛卡尔空间的位姿、速度分别为x、

Figure GDA0004090627790000111
其关节空间的角位置、角速度分别为q、
Figure GDA0004090627790000112
J为n自由度机器人的雅可比矩阵,其中x∈Rn
Figure GDA0004090627790000113
Figure GDA0004090627790000114
J∈Rm×n;冗余自由度机械臂的正运动学方程可用下式描述:Define the position and velocity of the end effector in Cartesian space as x,
Figure GDA0004090627790000111
The angular position and angular velocity of its joint space are q,
Figure GDA0004090627790000112
J is the Jacobian matrix of the n-DOF robot, where x∈R n ,
Figure GDA0004090627790000113
Figure GDA0004090627790000114
J∈R m×n ; The forward kinematics equation of the redundant degree of freedom manipulator can be described by the following formula:

Figure GDA0004090627790000115
Figure GDA0004090627790000115

式(1)也被称为机械臂运动学速度模型;Formula (1) is also called the robot arm kinematic velocity model;

考虑到最小二乘法的解,最优问题可列为:Taking into account the solution of the least squares method, the optimal problem can be stated as:

Figure GDA0004090627790000116
Figure GDA0004090627790000116

式(1)的解可通过寻找最佳

Figure GDA0004090627790000121
来解决;The solution of formula (1) can be obtained by finding the optimal
Figure GDA0004090627790000121
to solve;

Figure GDA0004090627790000122
Figure GDA0004090627790000122

因此,式(1)的伪逆解可表示为:Therefore, the pseudo-inverse solution of formula (1) can be expressed as:

Figure GDA0004090627790000123
Figure GDA0004090627790000123

式中J+——雅可比矩阵的伪逆Where J + ——pseudo-inverse of Jacobian matrix

I——单位矩阵I——Identity Matrix

Figure GDA0004090627790000124
——任意零空间矢量
Figure GDA0004090627790000124
——Arbitrary null space vector

Figure GDA0004090627790000125
——最小范数解,定义了手的动作
Figure GDA0004090627790000125
——Minimum norm solution, defining the hand movements

Figure GDA0004090627790000126
一一齐次解,在末端不产生动作
Figure GDA0004090627790000126
A homogeneous solution, no action is produced at the end

式(4)代表了末端执行器的位置和姿态控制;在式(4)中加入任意残差,可得到包含零空间的一般表达式;利用上述方程可在零向量上实现多任务优化;Formula (4) represents the position and attitude control of the end effector. By adding any residual to Formula (4), a general expression including the null space can be obtained. The above equation can be used to realize multi-task optimization on the null vector.

然而,上述方程忽略了雅可比矩阵的病态;正则化方程可通过添加额外的正则化值来修改,However, the above equation ignores the ill-conditioning of the Jacobian matrix; the regularization equation can be modified by adding additional regularization values,

Figure GDA0004090627790000127
Figure GDA0004090627790000127

其中λ≥0是加权矩阵,

Figure GDA0004090627790000128
为加权系数,并且满足where λ≥0 is the weighting matrix,
Figure GDA0004090627790000128
is the weighting coefficient and satisfies

Figure GDA0004090627790000129
Figure GDA0004090627790000129

上述方程的解就可表示为:The solution to the above equation can be expressed as:

Figure GDA00040906277900001210
Figure GDA00040906277900001210

式(7)也被称为冗余机械臂运动学模型;Formula (7) is also called the redundant manipulator kinematic model;

冗余机械臂零空间向量的位置相关标量指数的关节限制梯度方向的关节限制函数为:The joint constraint function of the joint constraint gradient direction of the position-dependent scalar index of the redundant manipulator null space vector is:

Figure GDA00040906277900001211
Figure GDA00040906277900001211

步骤2、建立通过奇异鲁棒解得出消除奇异性算法的任务优先级解决策略如下:Step 2: Establish the task priority solution strategy of eliminating singularity algorithm through singular robust solution as follows:

在雅可比矩阵的冗余机械臂求解中,优化任务是在主任务的零空间中实现的;反向任务运动学是在正向任务运动学基础上建立的:In the redundant manipulator solution of the Jacobian matrix, the optimization task is implemented in the null space of the main task; the inverse task kinematics is established based on the forward task kinematics:

Figure GDA00040906277900001212
Figure GDA00040906277900001212

其中

Figure GDA00040906277900001213
Figure GDA00040906277900001214
表示task1和task2in
Figure GDA00040906277900001213
and
Figure GDA00040906277900001214
Represents task1 and task2

从表达式(9)中得出冗余机械臂的逆运动学方程为:From expression (9), the inverse kinematics equation of the redundant manipulator is obtained as follows:

Figure GDA00040906277900001215
Figure GDA00040906277900001215

Task1作为主任务,task2作为辅助任务;也就是说,

Figure GDA00040906277900001216
是在
Figure GDA00040906277900001217
的零空间中实现的;冗余机械臂的最终逆运动学表达式如下所示:Task1 is the main task and task2 is the auxiliary task; that is,
Figure GDA00040906277900001216
is
Figure GDA00040906277900001217
The final inverse kinematics expression of the redundant manipulator is as follows:

Figure GDA0004090627790000131
Figure GDA0004090627790000131

其中

Figure GDA0004090627790000132
in
Figure GDA0004090627790000132

Figure GDA0004090627790000133
是投影矩阵,它给出了次要任务到主要任务的适用范围;
Figure GDA0004090627790000134
Figure GDA0004090627790000135
是所需的命令速度;
Figure GDA0004090627790000136
是主要任务,
Figure GDA0004090627790000137
是次要任务;
Figure GDA0004090627790000133
is the projection matrix, which gives the applicability of the secondary task to the primary task;
Figure GDA0004090627790000134
and
Figure GDA0004090627790000135
is the desired command speed;
Figure GDA0004090627790000136
is the main task,
Figure GDA0004090627790000137
It is a secondary task;

如果两个相关任务是相互依赖的,则相应的雅可比矩阵是奇异的;如果任务雅可比矩阵是奇异的,则相应的任务是不满足的;在这种情况下,雅可比相关矩阵将是奇点,定义为算法奇点;If two related tasks are interdependent, the corresponding Jacobian matrix is singular; if the task Jacobian matrix is singular, the corresponding tasks are unsatisfiable; in this case, the Jacobian correlation matrix will be a singular point, defined as an algorithmic singularity;

也就是说,如果That is, if

Figure GDA0004090627790000138
Figure GDA0004090627790000138

其中ρ(·)是矩阵的秩;where ρ(·) is the rank of the matrix;

很明显,算法的奇异性是由次要任务和主要任务之间的任务冲突引起的;此外,基于任务优先级的冗余机械臂逆运动学旨在提供更好的控制主要任务的有效性;It is obvious that the singularity of the algorithm is caused by the task conflict between the secondary and primary tasks; moreover, the redundant manipulator inverse kinematics based on task priority is designed to provide better effectiveness in controlling the primary task;

因此,让位置控制方向作为主要任务,从而使位置保证控制方向任务的准确性;然后建立通过奇异鲁棒解得出消除奇异性算法的任务优先级解决策略方程:Therefore, let the position control direction as the main task, so that the position guarantees the accuracy of the control direction task; then establish the task priority solution strategy equation of the algorithm to eliminate the singularity through the singular robust solution:

Figure GDA0004090627790000139
Figure GDA0004090627790000139

步骤3、建立奇异鲁棒解逆运动学分析模型如下:Step 3: Establish the singular robust solution inverse kinematics analysis model as follows:

基于雅可比伪逆解会出现运动学奇异性,这是由次矩阵引起的;对于运动奇异性问题,应该还要给出DLS(阻尼最小二乘)解;Based on the Jacobi pseudo-inverse solution, kinematic singularities will occur, which are caused by the secondary matrix; for the kinematic singularity problem, the DLS (damped least squares) solution should also be given;

对DLS解的成本函数可修改为:The cost function of the DLS solution can be modified as follows:

Figure GDA00040906277900001310
Figure GDA00040906277900001310

因此,上述方程的奇异鲁棒伪逆解可表示为:Therefore, the singular robust pseudo-inverse solution of the above equation can be expressed as:

Figure GDA00040906277900001311
Figure GDA00040906277900001311

式(15)即为奇异鲁棒解逆运动学分析模型,设λ=η2I,上述DLS解等价于附加正则化解,标量值η平衡了任务精度和奇异性;Formula (15) is the inverse kinematics analysis model of the singular robust solution. Let λ = η 2 I. The above DLS solution is equivalent to the additional regularization solution. The scalar value η balances the task accuracy and singularity.

对于雅可比矩阵伪逆解的计算,可给出雅可比矩阵的奇异值SVD分解形式For the calculation of the pseudo-inverse solution of the Jacobian matrix, the singular value SVD decomposition form of the Jacobian matrix can be given

J=UΣVT (16)J=UΣV T (16)

其中U∈∈Bm×m,V∈Rn×n,∑∈Rm×n,U是由列向量ui组成的一元矩阵,V是由列向量vi组成的一元矩阵,∑是m×n对角矩阵的块矩阵,该对角矩阵包含J的奇异值σi≥0以递减顺序包含n-m个零列向量;where U∈∈B m×m , V∈R n×n , ∑∈R m×n , U is a one-element matrix consisting of column vectors u i , V is a one-element matrix consisting of column vectors vi , ∑ is a block matrix of an m×n diagonal matrix containing the singular values σ i ≥ 0 of J in decreasing order and containing nm zero column vectors;

Figure GDA00040906277900001312
Figure GDA00040906277900001312

其中,r≤n是矩阵J的秩;Where r≤n is the rank of the matrix J;

对于运动奇异性,参考计算伪逆解时需要的奇异值分解SVD,大的产生的关节速度是由于最小的奇异值迅速接近0,如下式:For motion singularity, refer to the singular value decomposition SVD required to calculate the pseudo-inverse solution. The large joint velocity is due to the small singular value rapidly approaching 0, as shown in the following formula:

Figure GDA0004090627790000141
Figure GDA0004090627790000141

因子λ0将影响奇异性,λ0值越高,阻尼越大,联合速度就越接近奇异点;此外,定义可变阻尼因子的策略也不同;我们可得到The factor λ 0 will affect the singularity. The higher the value of λ 0 , the greater the damping and the closer the joint velocity is to the singularity point. In addition, the strategy for defining the variable damping factor is also different; we can get

Figure GDA0004090627790000142
Figure GDA0004090627790000142

从上式中,我们可看出参数δ>0监视最小的奇异值。From the above formula, we can see that the parameter δ>0 monitors the smallest singular value.

步骤4、建立多任务冗余机械臂的逆优先控制策略如下;Step 4: Establish the inverse priority control strategy of the multi-task redundant manipulator as follows;

引入逆优先投影矩阵

Figure GDA0004090627790000143
该矩阵包括了不依赖于第k个任务的最低优先级l-k-1任务的相应元素的零空间,所以得出Introducing the inverse priority projection matrix
Figure GDA0004090627790000143
This matrix includes the null space of the corresponding elements of the lowest priority lk-1 tasks that do not depend on the kth task, so we have

Figure GDA0004090627790000144
Figure GDA0004090627790000144

Figure GDA0004090627790000145
Figure GDA0004090627790000145

其中Ji|j是与j-th任务线性无关的i-th任务的所有组件相关联的雅可比矩阵;where Ji |j is the Jacobian matrix associated with all components of the i-th task that are linearly independent of the j-th task;

所以,优先级推导公式如下:Therefore, the priority derivation formula is as follows:

Figure GDA0004090627790000146
Figure GDA0004090627790000146

在上述推导中,k=1,l-1,…,1,其中l为正整数;初始值

Figure GDA0004090627790000147
In the above derivation, k = 1, l-1, ..., 1, where l is a positive integer; the initial value
Figure GDA0004090627790000147

为了给出线性无关雅可比矩阵J的一般计算形式,定义逆增广Jacobian矩阵为:In order to give the general calculation form of the linearly independent Jacobian matrix J, the inverse augmented Jacobian matrix is defined as:

Figure GDA0004090627790000148
Figure GDA0004090627790000148

就有了There you have it

Figure GDA0004090627790000149
Figure GDA0004090627790000149

其中

Figure GDA00040906277900001410
表示
Figure GDA00040906277900001411
的行;in
Figure GDA00040906277900001410
express
Figure GDA00040906277900001411
's line;

于是乎,

Figure GDA00040906277900001412
的伪逆解就可表示为:So,
Figure GDA00040906277900001412
The pseudo-inverse solution of can be expressed as:

Figure GDA00040906277900001413
Figure GDA00040906277900001413

Figure GDA00040906277900001414
and
Figure GDA00040906277900001414

其中,TR表示矩阵

Figure GDA0004090627790000151
的扩充;Among them, TR represents the matrix
Figure GDA0004090627790000151
expansion;

最终的反向优先级投影可写为:The final reverse priority projection can be written as:

Figure GDA0004090627790000152
Figure GDA0004090627790000152

于是,我们就可得出伪逆解的表达式:Therefore, we can derive the expression of the pseudo-inverse solution:

Figure GDA0004090627790000153
Figure GDA0004090627790000153

建立多任务冗余机械臂的逆优先控制策略方程如下:The inverse priority control strategy equation for the multi-task redundant manipulator is established as follows:

Figure GDA0004090627790000154
Figure GDA0004090627790000154

步骤5、对具有主要任务和次要任务的冗余机械臂的反向控制方程进行简化处理如下:Step 5: Simplify the inverse control equation of the redundant manipulator with primary and secondary tasks as follows:

对于六自由度或七自由度冗余机械臂来说,没有足够的自由度来完成多个层次的任务;实行双任务优先控制是必要的;也就是说,机械手的运动控制是主要任务和次要任务;For a six- or seven-degree-of-freedom redundant manipulator, there are not enough degrees of freedom to complete multiple levels of tasks; it is necessary to implement dual-task priority control; that is, the motion control of the manipulator is the primary task and the secondary task;

具有主要任务和次要任务的冗余机械臂的反向控制方程如下The inverse control equation for a redundant manipulator with primary and secondary tasks is as follows

Figure GDA0004090627790000155
Figure GDA0004090627790000155

上面的公式与前面的表达式(11)有很大的不同,但是算法框架是相似的;在上面的方程中,

Figure GDA0004090627790000156
是次要任务,
Figure GDA0004090627790000157
是主任务;主任务在主任务的指定零空间中实现;反向优先级的核心点是投影矩阵
Figure GDA0004090627790000158
的计算;
Figure GDA0004090627790000159
的表达式如公式(30):The above formula is quite different from the previous expression (11), but the algorithm framework is similar; in the above equation,
Figure GDA0004090627790000156
It's a secondary task.
Figure GDA0004090627790000157
is the main task; the main task is implemented in the specified null space of the main task; the core point of the reverse priority is the projection matrix
Figure GDA0004090627790000158
Calculation of
Figure GDA0004090627790000159
The expression of is as follows:

Figure GDA00040906277900001510
Figure GDA00040906277900001510

利用前述公式(22)-(28)中的导子,可得到简化后具有主要任务和次要任务的冗余机械臂的反向控制方程:Using the derivations in the above formulas (22)-(28), the simplified inverse control equations of the redundant manipulator with primary and secondary tasks can be obtained:

Figure GDA00040906277900001511
Figure GDA00040906277900001511

步骤6、建立机械手的逆优先力控制策略如下:Step 6: Establish the reverse priority force control strategy of the manipulator as follows:

力控制空间中机械手的动力学可写成:The dynamics of the manipulator in the force control space can be written as:

Figure GDA00040906277900001512
Figure GDA00040906277900001512

其中X是笛卡尔空间中的位置,M(X是惯性矩阵,

Figure GDA00040906277900001513
是非线性力,F是输入控制力,Fe是接触力;where X is the position in Cartesian space, M (X is the inertia matrix,
Figure GDA00040906277900001513
is the nonlinear force, F is the input control force, and Fe is the contact force;

此外,基于雅可比矩阵的变换可得到输入关节力矩In addition, the input joint torque can be obtained based on the transformation of the Jacobian matrix

τ=JT(q)F (33)τ=J T (q)F (33)

操纵器在力控制空间中的期望运动方程可定义如下:The desired equation of motion of the manipulator in the force control space can be defined as follows:

Figure GDA00040906277900001514
Figure GDA00040906277900001514

其中Md和Bd是惯性和阻尼矩阵;Fd是指令力,Fe是接触力;Where Md and Bd are the inertia and damping matrices; Fd is the command force, and Fe is the contact force;

力控制的动力学方案如图1所示;The kinetic scheme of force control is shown in Figure 1 ;

因此,环境和操纵器响应之间的关系可写成Therefore, the relationship between the environment and the manipulator response can be written as

Figure GDA0004090627790000161
Figure GDA0004090627790000161

上述两个方程的组合如下The combination of the above two equations is as follows

Figure GDA0004090627790000162
Figure GDA0004090627790000162

从上面的方程可看出,如果Me、Be和Ke已知,则Md和Bd的调整将影响系统响应;From the above equations, it can be seen that if Me , Be , and Ke are known, then adjustments to Md and Bd will affect the system response;

力控制使机械手能够与环境或人类相互作用;另外,在某些情况下,没有必要实现全方位的力控制,也没有必要保证全方位的力控制,也就是说,有时我们只是想保证某个方向的力跟踪控制精度;Force control enables the manipulator to interact with the environment or humans; in addition, in some cases, it is not necessary to achieve full-range force control, nor is it necessary to ensure full-range force control, that is, sometimes we just want to ensure the force tracking control accuracy in a certain direction;

例如,当机械手与刨床相互作用时,只需在垂直方向保持精确的力跟踪控制,而另一个方向不需要精确的力跟踪控制;在其他情况下,位置方向力控制比姿态方向力控制更为重要;For example, when the manipulator interacts with a planer, it only needs to maintain precise force tracking control in the vertical direction, while the other direction does not need precise force tracking control; in other cases, position direction force control is more important than attitude direction force control;

因此有必要对机械手进行分级力控制;也就是说,有必要给出一个新的层次力控制框架;从上面的方程我们可得到期望的层次力控制关系如下Therefore, it is necessary to perform hierarchical force control on the manipulator; that is, it is necessary to give a new hierarchical force control framework; from the above equations we can get the expected hierarchical force control relationship as follows

Figure GDA0004090627790000163
Figure GDA0004090627790000163

Figure GDA0004090627790000164
Figure GDA0004090627790000164

所以这两个方程的积分公式可写成So the integral formula of these two equations can be written as

Figure GDA0004090627790000165
Figure GDA0004090627790000165

Figure GDA0004090627790000166
Figure GDA0004090627790000166

如果机械手末端执行器能够跟踪期望的笛卡尔速度为

Figure GDA0004090627790000167
Figure GDA0004090627790000168
则可实现机械手的精确力控制;笛卡尔速度与关节速度的关系应借鉴逆优先控制;因此,可得到机械手的逆优先力控制策略的方程:If the robot end effector can track the desired Cartesian velocity
Figure GDA0004090627790000167
and
Figure GDA0004090627790000168
Then the precise force control of the manipulator can be achieved; the relationship between Cartesian speed and joint speed should refer to the inverse priority control; therefore, the equation of the inverse priority force control strategy of the manipulator can be obtained:

Figure GDA0004090627790000169
Figure GDA0004090627790000169

上述方程所要求的关节速度将保证机械手的力控制;值得一提的是,上述力控制律只是速度级控制律,它依赖于内速度环控制;如果内位置控制效果良好,则可实现精确的力控制;由于内速度环控制可实现低频位置跟踪,所以外力环可实现低频力跟踪。The joint speed required by the above equation will ensure the force control of the manipulator; it is worth mentioning that the above force control law is only a speed-level control law, which relies on the inner speed loop control; if the inner position control effect is good, accurate force control can be achieved; since the inner speed loop control can achieve low-frequency position tracking, the outer force loop can achieve low-frequency force tracking.

步骤7、采用关节速度来解决机械手逆优先阻抗控制中的外力与关节加速度之间的关系,从而得到机械手的逆优先阻抗控制保证的实现方式如下:Step 7: Use joint velocity to solve the relationship between external force and joint acceleration in the inverse priority impedance control of the manipulator, so as to obtain the implementation method of the inverse priority impedance control guarantee of the manipulator as follows:

当机械手实施力控制时,机械手在一定程度上起到了发起者的作用,也就是说,机械手已经做好了响应外部环境的准备;当机械臂When the manipulator implements force control, the manipulator plays the role of an initiator to a certain extent, that is, the manipulator is ready to respond to the external environment; when the manipulator

作为阻抗控制模型工作时,机械臂会被动地响应外力;阻抗控制的动力学方案如图2所示;When working as an impedance-controlled model, the robot passively responds to external forces; the dynamic scheme of impedance control is shown in Figure 2;

外力与关节加速度的对应阻抗关系可表示为The corresponding impedance relationship between external force and joint acceleration can be expressed as

Figure GDA0004090627790000171
Figure GDA0004090627790000171

Figure GDA0004090627790000172
Figure GDA0004090627790000172

参考速度可表示为The reference speed can be expressed as

Figure GDA0004090627790000173
Figure GDA0004090627790000173

Figure GDA0004090627790000174
Figure GDA0004090627790000174

因此,机械手的逆优先阻抗控制保证的表达式为:Therefore, the expression guaranteed by the reverse priority impedance control of the manipulator is:

Figure GDA0004090627790000175
Figure GDA0004090627790000175

步骤8、将位置控制空间的逆优先计算扩展到力控制空间的逆优先计算,从而获得机械手速度级逆优先阻抗控制的总体框架实现方式如下:Step 8: Extend the inverse priority calculation of the position control space to the inverse priority calculation of the force control space, so as to obtain the overall framework of the manipulator speed level inverse priority impedance control. The implementation method is as follows:

混合阻抗应用就是上述两种策略的结合,即笛卡尔任务可分为两种情况:第一种是位置控制子空间,阻抗控制是在该子空间中实现的;第二个是力控制子空间,力控制在该子空间中实现;Hybrid impedance application is a combination of the above two strategies, i.e., the Cartesian task can be divided into two cases: the first is the position control subspace, in which impedance control is implemented; the second is the force control subspace, in which force control is implemented;

因此选择一个选择矩阵;外力与位置响应的关系如下Therefore, a selection matrix is selected; the relationship between external force and position response is as follows

Figure GDA0004090627790000176
Figure GDA0004090627790000176

Figure GDA0004090627790000177
Figure GDA0004090627790000177

所以期望速度的简化形式可表示为So the simplified form of the expected speed can be expressed as

Figure GDA0004090627790000178
Figure GDA0004090627790000178

Figure GDA0004090627790000179
Figure GDA0004090627790000179

然后我们得到了基于反向优先级的解决方案Then we get the solution based on reverse priority

Figure GDA00040906277900001710
Figure GDA00040906277900001710

混合阻抗控制的动力学方案图3所示;The kinetic scheme of hybrid impedance control is shown in Figure 3;

考虑到n层任务,相应的阻抗控制任务也属于n层框架,因此,机械手速度级逆优先阻抗控制的总体框架表达式如下Considering the n-layer task, the corresponding impedance control task also belongs to the n-layer framework. Therefore, the overall framework expression of the manipulator speed level reverse priority impedance control is as follows:

Figure GDA0004090627790000181
Figure GDA0004090627790000181

表达式(52)解决了将位置控制空间的逆优先计算扩展到力控制空间的逆优先计算的机械手逆优先混合阻抗控制中去;能在不同层次结构下可使机械手的冗余机械臂实现期望的阻抗控制任务。Expression (52) solves the problem of extending the inverse priority calculation of the position control space to the inverse priority calculation of the force control space in the manipulator inverse priority mixed impedance control; it can enable the redundant manipulator arms of the manipulator to achieve the desired impedance control task under different hierarchical structures.

Claims (1)

1.一种移动机械手冗余机械臂逆优先级阻抗控制方法,其特征在于,包括如下步骤:1. A method for reverse priority impedance control of redundant mechanical arms of a mobile manipulator, characterized in that it comprises the following steps: 步骤1、建立冗余机械臂运动学模型,并给出冗余机械臂零空间向量的梯度方向策略;Step 1: Establish a kinematic model of the redundant manipulator and provide a gradient direction strategy of the redundant manipulator null space vector; 步骤2、建立通过奇异鲁棒解得出消除奇异性算法的任务优先级解决策略;Step 2: Establish a task priority solution strategy for eliminating singularity through a singular robust solution; 步骤3、建立奇异鲁棒解逆运动学分析模型;Step 3, establishing a singular robust solution inverse kinematics analysis model; 步骤4、建立多任务冗余机械臂的逆优先控制策略;Step 4: Establish an inverse priority control strategy for the multi-task redundant manipulator; 步骤5、对具有主要任务和次要任务的冗余机械臂的反向控制方程进行简化处理;Step 5, simplifying the inverse control equation of the redundant manipulator with the primary task and the secondary task; 步骤6、建立机械手的逆优先力控制策略;Step 6: Establish the inverse priority force control strategy of the manipulator; 步骤7、采用关节速度来解决机械手逆优先阻抗控制中的外力与关节加速度之间的关系,从而得到机械手的逆优先阻抗控制保证;Step 7: Use joint velocity to solve the relationship between external force and joint acceleration in the inverse priority impedance control of the manipulator, so as to obtain the inverse priority impedance control guarantee of the manipulator; 步骤8、将位置控制空间的逆优先计算扩展到力控制空间的逆优先计算,从而获得机械手速度级逆优先阻抗控制的总体框架;Step 8: Extend the inverse priority calculation of the position control space to the inverse priority calculation of the force control space, thereby obtaining the overall framework of the manipulator speed level inverse priority impedance control; 建立冗余机械臂运动学模型,并给出冗余机械臂零空间向量的梯度方向策略实现过程如下:The kinematic model of the redundant manipulator is established, and the gradient direction strategy of the redundant manipulator null space vector is given as follows: 定义末端执行器在笛卡尔空间的位姿、速度分别为x、
Figure FDA0004090627770000011
其关节空间的角位置、角速度分别为q、
Figure FDA0004090627770000012
J为n自由度机器人的雅可比矩阵,其中x∈Rn
Figure FDA0004090627770000013
J∈Rmn;冗余自由度机械臂的正运动学方程可用下式描述:
Define the position and velocity of the end effector in Cartesian space as x,
Figure FDA0004090627770000011
The angular position and angular velocity of its joint space are q,
Figure FDA0004090627770000012
J is the Jacobian matrix of the n-DOF robot, where x∈R n ,
Figure FDA0004090627770000013
J∈R mn ; The forward kinematics equation of the redundant degree of freedom manipulator can be described by the following formula:
Figure FDA0004090627770000014
Figure FDA0004090627770000014
式(1)也被称为机械臂运动学速度模型;Formula (1) is also called the robot arm kinematic velocity model; 考虑到最小二乘法的解,最优问题可列为:Taking into account the solution of the least squares method, the optimal problem can be stated as:
Figure FDA0004090627770000021
Figure FDA0004090627770000021
式(1)的解可通过寻找最佳
Figure FDA0004090627770000022
来解决;
The solution of formula (1) can be obtained by finding the optimal
Figure FDA0004090627770000022
to solve;
Figure FDA0004090627770000023
Figure FDA0004090627770000023
因此,式(1)的伪逆解可表示为:Therefore, the pseudo-inverse solution of formula (1) can be expressed as:
Figure FDA0004090627770000024
Figure FDA0004090627770000024
式中J+——雅可比矩阵的伪逆Where J + ——pseudo-inverse of Jacobian matrix I——单位矩阵I——Identity Matrix
Figure FDA0004090627770000025
——任意零空间矢量
Figure FDA0004090627770000025
——Arbitrary null space vector
Figure FDA0004090627770000026
——最小范数解,定义了手的动作
Figure FDA0004090627770000026
——Minimum norm solution, defining the hand movements
Figure FDA0004090627770000027
——齐次解,在末端不产生动作
Figure FDA0004090627770000027
——Homogeneous solution, no action is produced at the end
式(4)代表了末端执行器的位置和姿态控制;在式(4)中加入任意残差,可得到包含零空间的一般表达式;利用上述方程可在零向量上实现多任务优化;Formula (4) represents the position and attitude control of the end effector. By adding any residual to Formula (4), a general expression including the null space can be obtained. The above equation can be used to realize multi-task optimization on the null vector. 然而,上述方程忽略了雅可比矩阵的病态;正则化方程可通过添加额外的正则化值来修改,However, the above equation ignores the ill-conditioning of the Jacobian matrix; the regularization equation can be modified by adding additional regularization values,
Figure FDA0004090627770000028
Figure FDA0004090627770000028
其中λ≥0是加权矩阵,
Figure FDA0004090627770000029
为加权系数,并且满足
where λ≥0 is the weighting matrix,
Figure FDA0004090627770000029
is the weighting coefficient and satisfies
Figure FDA00040906277700000210
Figure FDA00040906277700000210
上述方程的解就可表示为:The solution to the above equation can be expressed as:
Figure FDA00040906277700000211
Figure FDA00040906277700000211
式(7)也被称为冗余机械臂运动学模型;Formula (7) is also called the redundant manipulator kinematic model; 冗余机械臂零空间向量的位置相关标量指数的关节限制梯度方向的关节限制函数为:The joint constraint function of the joint constraint gradient direction of the position-dependent scalar index of the redundant manipulator null space vector is:
Figure FDA0004090627770000031
Figure FDA0004090627770000031
建立通过奇异鲁棒解得出消除奇异性算法的任务优先级解决策略如下:The task priority solution strategy for establishing the algorithm for eliminating singularities through singular robust solutions is as follows: 在雅可比矩阵的冗余机械臂求解中,优化任务是在主任务的零空间中实现的;反向任务运动学是在正向任务运动学基础上建立的:In the redundant manipulator solution of the Jacobian matrix, the optimization task is implemented in the null space of the main task; the inverse task kinematics is established based on the forward task kinematics:
Figure FDA0004090627770000032
Figure FDA0004090627770000032
其中
Figure FDA0004090627770000033
Figure FDA0004090627770000034
表示task1和task2
in
Figure FDA0004090627770000033
and
Figure FDA0004090627770000034
Represents task1 and task2
从表达式(9)中得出冗余机械臂的逆运动学方程为:From expression (9), the inverse kinematics equation of the redundant manipulator is obtained as follows:
Figure FDA0004090627770000035
Figure FDA0004090627770000035
Task1作为主任务,task2作为辅助任务;也就是说,task2
Figure FDA0004090627770000036
是在task1
Figure FDA0004090627770000037
的零空间中实现的;冗余机械臂的最终逆运动学表达式如下所示:
Task1 is the main task and task2 is the auxiliary task; that is, task2
Figure FDA0004090627770000036
It is in task1
Figure FDA0004090627770000037
The final inverse kinematics expression of the redundant manipulator is as follows:
Figure FDA0004090627770000038
Figure FDA0004090627770000038
其中
Figure FDA0004090627770000039
in
Figure FDA0004090627770000039
Figure FDA00040906277700000310
是投影矩阵,它给出了次要任务到主要任务的适用范围;
Figure FDA00040906277700000311
Figure FDA00040906277700000312
是所需的命令速度;
Figure FDA00040906277700000313
是主要任务,
Figure FDA00040906277700000314
是次要任务;
Figure FDA00040906277700000310
is the projection matrix, which gives the applicability of the secondary task to the primary task;
Figure FDA00040906277700000311
and
Figure FDA00040906277700000312
is the desired command speed;
Figure FDA00040906277700000313
is the main task,
Figure FDA00040906277700000314
It is a secondary task;
如果两个相关任务是相互依赖的,则相应的雅可比矩阵是奇异的;如果任务雅可比矩阵是奇异的,则相应的任务是不满足的;在这种情况下,雅可比相关矩阵将是奇点,定义为算法奇点;If two related tasks are interdependent, the corresponding Jacobian matrix is singular; if the task Jacobian matrix is singular, the corresponding tasks are unsatisfiable; in this case, the Jacobian correlation matrix will be a singular point, defined as an algorithmic singularity; 也就是说,如果That is, if
Figure FDA0004090627770000041
Figure FDA0004090627770000041
其中p(·)是矩阵的秩;where p(·) is the rank of the matrix; 很明显,算法的奇异性是由次要任务和主要任务之间的任务冲突引起的;此外,基于任务优先级的冗余机械臂逆运动学旨在提供更好的控制主要任务的有效性;It is obvious that the singularity of the algorithm is caused by the task conflict between the secondary and primary tasks; moreover, the redundant manipulator inverse kinematics based on task priority is designed to provide better effectiveness in controlling the primary task; 因此,让位置控制方向作为主要任务,从而使位置保证控制方向任务的准确性;然后建立通过奇异鲁棒解得出消除奇异性算法的任务优先级解决策略方程:Therefore, let the position control direction as the main task, so that the position guarantees the accuracy of the control direction task; then establish the task priority solution strategy equation of the algorithm to eliminate the singularity through the singular robust solution:
Figure FDA0004090627770000042
Figure FDA0004090627770000042
建立奇异鲁棒解逆运动学分析模型如下:The singular robust solution inverse kinematics analysis model is established as follows: 基于雅可比伪逆解会出现运动学奇异性,这是由次矩阵引起的;对于运动奇异性问题,应该还要给出DLS(阻尼最小二乘)解;Based on the Jacobi pseudo-inverse solution, kinematic singularities will occur, which are caused by the secondary matrix; for the kinematic singularity problem, the DLS (damped least squares) solution should also be given; 对DLS解的成本函数可修改为:The cost function of the DLS solution can be modified as follows:
Figure FDA0004090627770000043
Figure FDA0004090627770000043
因此,上述方程的奇异鲁棒伪逆解可表示为:Therefore, the singular robust pseudo-inverse solution of the above equation can be expressed as:
Figure FDA0004090627770000044
Figure FDA0004090627770000044
式(15)即为奇异鲁棒解逆运动学分析模型,设λ=η2I,上述DLS解等价于附加正则化解,标量值η平衡了任务精度和奇异性;Formula (15) is the inverse kinematics analysis model of the singular robust solution. Let λ = η 2 I. The above DLS solution is equivalent to the additional regularization solution. The scalar value η balances the task accuracy and singularity. 对于雅可比矩阵伪逆解的计算,可给出雅可比矩阵的奇异值SVD分解形式For the calculation of the pseudo-inverse solution of the Jacobian matrix, the singular value SVD decomposition form of the Jacobian matrix can be given J=UΣVT (16)J=UΣV T (16) 其中U∈Rm×n,V∈Rn×n,∑∈Rm×n,U是由列向量ui组成的一元矩阵,V是由列向量vi组成的一元矩阵,∑是m×n对角矩阵的块矩阵,该对角矩阵包含J的奇异值σi≥0以递减顺序包含n-m个零列向量;where U∈R m×n , V∈R n×n , ∑∈R m×n , U is a one-element matrix consisting of column vectors u i , V is a one-element matrix consisting of column vectors vi , ∑ is a block matrix of an m×n diagonal matrix containing the singular values σ i ≥ 0 of J in decreasing order and containing nm zero column vectors;
Figure FDA0004090627770000051
Figure FDA0004090627770000051
其中,r≤m是矩阵J的秩;Where r≤m is the rank of the matrix J; 对于运动奇异性,参考计算伪逆解时需要的奇异值分解SVD,大的产生的关节速度是由于最小的奇异值迅速接近0,如下式:For motion singularity, refer to the singular value decomposition SVD required to calculate the pseudo-inverse solution. The large joint velocity is due to the small singular value rapidly approaching 0, as shown in the following formula:
Figure FDA0004090627770000052
Figure FDA0004090627770000052
因子λ0将影响奇异性,λ0值越高,阻尼越大,联合速度就越接近奇异点;此外,定义可变阻尼因子的策略也不同;我们可得到The factor λ 0 will affect the singularity. The higher the value of λ 0 , the greater the damping and the closer the joint velocity is to the singularity point. In addition, the strategy for defining the variable damping factor is also different; we can get
Figure FDA0004090627770000053
Figure FDA0004090627770000053
从上式中,我们可看出参数δ>0监视最小的奇异值;From the above formula, we can see that the parameter δ>0 monitors the smallest singular value; 建立多任务冗余机械臂的逆优先控制策略如下;The inverse priority control strategy of the multi-task redundant manipulator is established as follows; 引入逆优先投影矩阵
Figure FDA0004090627770000054
该矩阵包括了不依赖于第k个任务的最低优先级l-k-1任务的相应元素的零空间,所以得出
Introducing the inverse priority projection matrix
Figure FDA0004090627770000054
This matrix includes the null space of the corresponding elements of the lowest priority lk-1 tasks that do not depend on the kth task, so we have
Figure FDA0004090627770000055
Figure FDA0004090627770000055
Figure FDA0004090627770000056
Figure FDA0004090627770000056
其中Ji|j是与j-th任务线性无关的i-th任务的所有组件相关联的雅可比矩阵;where Ji |j is the Jacobian matrix associated with all components of the i-th task that are linearly independent of the j-th task; 所以,优先级推导公式如下:Therefore, the priority derivation formula is as follows:
Figure FDA0004090627770000061
Figure FDA0004090627770000061
在上述推导中,k=l,l-1,…,1,其中1为正整数;初始值
Figure FDA0004090627770000062
In the above derivation, k = l, l-1, ..., 1, where 1 is a positive integer; the initial value
Figure FDA0004090627770000062
为了给出线性无关雅可比矩阵J的一般计算形式,定义逆增广Jacobian矩阵为:In order to give the general calculation form of the linearly independent Jacobian matrix J, the inverse augmented Jacobian matrix is defined as:
Figure FDA0004090627770000063
Figure FDA0004090627770000063
就有了There you have it
Figure FDA0004090627770000064
Figure FDA0004090627770000064
其中
Figure FDA0004090627770000065
表示
Figure FDA0004090627770000066
的行;
in
Figure FDA0004090627770000065
express
Figure FDA0004090627770000066
's line;
于是乎,
Figure FDA0004090627770000067
的伪逆解就可表示为:
So,
Figure FDA0004090627770000067
The pseudo-inverse solution of can be expressed as:
Figure FDA0004090627770000068
Figure FDA0004090627770000068
Figure FDA0004090627770000069
and
Figure FDA0004090627770000069
其中,Tk表示矩阵
Figure FDA00040906277700000610
的扩充;
Where T k represents the matrix
Figure FDA00040906277700000610
expansion;
最终的反向优先级投影可写为:The final reverse priority projection can be written as:
Figure FDA00040906277700000611
Figure FDA00040906277700000611
于是,我们就可得出伪逆解的表达式:Therefore, we can derive the expression of the pseudo-inverse solution:
Figure FDA0004090627770000071
Figure FDA0004090627770000071
建立多任务冗余机械臂的逆优先控制策略方程如下:The inverse priority control strategy equation for the multi-task redundant manipulator is established as follows:
Figure FDA0004090627770000072
Figure FDA0004090627770000072
对具有主要任务和次要任务的冗余机械臂的反向控制方程进行简化处理如下:The inverse control equations of a redundant manipulator with primary and secondary tasks are simplified as follows: 对于六自由度或七自由度冗余机械臂来说,没有足够的自由度来完成多个层次的任务;实行双任务优先控制是必要的;也就是说,机械手的运动控制是主要任务和次要任务;For a six- or seven-degree-of-freedom redundant manipulator, there are not enough degrees of freedom to complete multiple levels of tasks; it is necessary to implement dual-task priority control; that is, the motion control of the manipulator is the primary task and the secondary task; 具有主要任务和次要任务的冗余机械臂的反向控制方程如下The inverse control equation for a redundant manipulator with primary and secondary tasks is as follows
Figure FDA0004090627770000073
Figure FDA0004090627770000073
上面的公式与前面的表达式(11)有很大的不同,但是算法框架是相似的;在上面的方程中,
Figure FDA0004090627770000074
是次要任务,
Figure FDA0004090627770000075
是主任务;主任务在主任务的指定零空间中实现;反向优先级的核心点是投影矩阵
Figure FDA0004090627770000076
的计算;
Figure FDA0004090627770000077
的表达式如公式(30):
The above formula is quite different from the previous expression (11), but the algorithm framework is similar; in the above equation,
Figure FDA0004090627770000074
It's a secondary task.
Figure FDA0004090627770000075
is the main task; the main task is implemented in the specified null space of the main task; the core point of the reverse priority is the projection matrix
Figure FDA0004090627770000076
Calculation of
Figure FDA0004090627770000077
The expression of is as follows:
Figure FDA0004090627770000078
Figure FDA0004090627770000078
利用前述公式(22)-(28)中的导子,可得到简化后具有主要任务和次要任务的冗余机械臂的反向控制方程:Using the derivations in the above formulas (22)-(28), the simplified inverse control equations of the redundant manipulator with primary and secondary tasks can be obtained:
Figure FDA0004090627770000079
Figure FDA0004090627770000079
建立机械手的逆优先力控制策略如下:The inverse priority force control strategy of the manipulator is established as follows: 力控制空间中机械手的动力学可写成:The dynamics of the manipulator in the force control space can be written as:
Figure FDA0004090627770000081
Figure FDA0004090627770000081
其中X是笛卡尔空间中的位置,M(X)是惯性矩阵,
Figure FDA0004090627770000082
是非线性力,F是输入控制力,Fe是接触力;
where X is the position in Cartesian space, M(X) is the inertia matrix,
Figure FDA0004090627770000082
is the nonlinear force, F is the input control force, and Fe is the contact force;
此外,基于雅可比矩阵的变换可得到输入关节力矩In addition, the input joint torque can be obtained based on the transformation of the Jacobian matrix τ=JT(q)F (33)τ=J T (q)F (33) 操纵器在力控制空间中的期望运动方程可定义如下:The desired equation of motion of the manipulator in the force control space can be defined as follows:
Figure FDA0004090627770000083
Figure FDA0004090627770000083
其中Md和Bd是惯性和阻尼矩阵;Fd是指令力,Fe是接触力;Where Md and Bd are the inertia and damping matrices; Fd is the command force, and Fe is the contact force; 因此,环境和操纵器响应之间的关系可写成Therefore, the relationship between the environment and the manipulator response can be written as
Figure FDA0004090627770000084
Figure FDA0004090627770000084
上述两个方程的组合如下The combination of the above two equations is as follows
Figure FDA0004090627770000085
Figure FDA0004090627770000085
从上面的方程可看出,如果Me、Be和Ke已知,则Md和Bd的调整将影响系统响应;From the above equations, it can be seen that if Me , Be , and Ke are known, then adjustments to Md and Bd will affect the system response; 力控制使机械手能够与环境或人类相互作用;另外,在某些情况下,没有必要实现全方位的力控制,也没有必要保证全方位的力控制,也就是说,有时我们只是想保证某个方向的力跟踪控制精度;Force control enables the manipulator to interact with the environment or humans; in addition, in some cases, it is not necessary to achieve full-range force control, nor is it necessary to ensure full-range force control, that is, sometimes we just want to ensure the force tracking control accuracy in a certain direction; 因此有必要对机械手进行分级力控制;也就是说,有必要给出一个新的层次力控制框架;从上面的方程我们可得到期望的层次力控制关系如下Therefore, it is necessary to perform hierarchical force control on the manipulator; that is, it is necessary to give a new hierarchical force control framework; from the above equations we can get the expected hierarchical force control relationship as follows
Figure FDA0004090627770000086
Figure FDA0004090627770000086
Figure FDA0004090627770000087
Figure FDA0004090627770000087
所以这两个方程的积分公式可写成So the integral formula of these two equations can be written as
Figure FDA0004090627770000091
Figure FDA0004090627770000091
Figure FDA0004090627770000092
Figure FDA0004090627770000092
如果机械手末端执行器能够跟踪期望的笛卡尔速度为
Figure FDA0004090627770000093
Figure FDA0004090627770000094
则可实现机械手的精确力控制;笛卡尔速度与关节速度的关系应借鉴逆优先控制;因此,可得到机械手的逆优先力控制策略的方程:
If the robot end effector can track the desired Cartesian velocity
Figure FDA0004090627770000093
and
Figure FDA0004090627770000094
Then the precise force control of the manipulator can be achieved; the relationship between Cartesian speed and joint speed should refer to the inverse priority control; therefore, the equation of the inverse priority force control strategy of the manipulator can be obtained:
Figure FDA0004090627770000095
Figure FDA0004090627770000095
上述方程所要求的关节速度将保证机械手的力控制;值得一提的是,上述力控制律只是速度级控制律,它依赖于内速度环控制;如果内位置控制效果良好,则可实现精确的力控制;由于内速度环控制可实现低频位置跟踪,所以外力环可实现低频力跟踪;The joint speed required by the above equation will ensure the force control of the manipulator; it is worth mentioning that the above force control law is only a speed-level control law, which depends on the inner speed loop control; if the inner position control effect is good, accurate force control can be achieved; since the inner speed loop control can achieve low-frequency position tracking, the outer force loop can achieve low-frequency force tracking; 采用关节速度来解决机械手逆优先阻抗控制中的外力与关节加速度之间的关系,从而得到机械手的逆优先阻抗控制保证的实现方式如下:The joint velocity is used to solve the relationship between the external force and the joint acceleration in the inverse priority impedance control of the manipulator, so that the inverse priority impedance control of the manipulator is guaranteed to be implemented as follows: 当机械手实施力控制时,机械手在一定程度上起到了发起者的作用,也就是说,机械手已经做好了响应外部环境的准备;当机械臂作为阻抗控制模型工作时,机械臂会被动地响应外力;When the manipulator implements force control, the manipulator plays the role of an initiator to a certain extent, that is, the manipulator is ready to respond to the external environment; when the manipulator works as an impedance control model, the manipulator responds to the external force passively; 外力与关节加速度的对应阻抗关系可表示为The corresponding impedance relationship between external force and joint acceleration can be expressed as
Figure FDA0004090627770000096
Figure FDA0004090627770000096
Figure FDA0004090627770000097
Figure FDA0004090627770000097
参考速度可表示为The reference speed can be expressed as
Figure FDA0004090627770000101
Figure FDA0004090627770000101
Figure FDA0004090627770000102
Figure FDA0004090627770000102
因此,机械手的逆优先阻抗控制保证的表达式为:Therefore, the expression guaranteed by the reverse priority impedance control of the manipulator is:
Figure FDA0004090627770000103
Figure FDA0004090627770000103
将位置控制空间的逆优先计算扩展到力控制空间的逆优先计算,从而获得机械手速度级逆优先阻抗控制的总体框架实现方式如下:The inverse priority calculation of the position control space is extended to the inverse priority calculation of the force control space, so as to obtain the overall framework of the manipulator speed level inverse priority impedance control as follows: 混合阻抗应用就是上述两种策略的结合,即笛卡尔任务可分为两种情况:第一种是位置控制子空间,阻抗控制是在该子空间中实现的;第二个是力控制子空间,力控制在该子空间中实现;Hybrid impedance application is a combination of the above two strategies, i.e., the Cartesian task can be divided into two cases: the first is the position control subspace, in which impedance control is implemented; the second is the force control subspace, in which force control is implemented; 因此选择一个选择矩阵;外力与位置响应的关系如下Therefore, a selection matrix is selected; the relationship between external force and position response is as follows
Figure FDA0004090627770000104
Figure FDA0004090627770000104
Figure FDA0004090627770000105
Figure FDA0004090627770000105
所以期望速度的简化形式可表示为So the simplified form of the expected speed can be expressed as
Figure FDA0004090627770000106
Figure FDA0004090627770000106
Figure FDA0004090627770000107
Figure FDA0004090627770000107
然后我们得到了基于反向优先级的解决方案Then we get the solution based on reverse priority
Figure FDA0004090627770000108
Figure FDA0004090627770000108
考虑到n层任务,相应的阻抗控制任务也属于n层框架,因此,机械手速度级逆优先阻抗控制的总体框架表达式如下Considering the n-layer task, the corresponding impedance control task also belongs to the n-layer framework. Therefore, the overall framework expression of the manipulator speed level reverse priority impedance control is as follows:
Figure FDA0004090627770000111
Figure FDA0004090627770000111
表达式(52)解决了将位置控制空间的逆优先计算扩展到力控制空间的逆优先计算的机械手逆优先混合阻抗控制中去;能在不同层次结构下可使机械手的冗余机械臂实现期望的阻抗控制任务;Expression (52) solves the problem of extending the inverse priority calculation of the position control space to the inverse priority calculation of the force control space in the manipulator inverse priority mixed impedance control; it can enable the redundant manipulator arms of the manipulator to achieve the desired impedance control task under different hierarchical structures; 适用于移动机械手冗余机械臂逆优先级阻抗控制方法的移动机械手冗余机械臂逆优先级阻抗控制系统,所述移动机械手冗余机械臂逆优先级阻抗控制系统包括机械手和控制机械手的控制台(S31);还包括可移动的移动平台(S41);所述机械手包括机械臂、安装座(S1)、竖直柱(S2)、输出抓手(S12)和竖直气缸(S23);安装座固定在移动平台上;A mobile manipulator redundant manipulator arm reverse priority impedance control system applicable to a mobile manipulator redundant manipulator arm reverse priority impedance control method, the mobile manipulator redundant manipulator arm reverse priority impedance control system comprising a manipulator and a control console (S31) for controlling the manipulator; and also comprising a movable mobile platform (S41); the manipulator comprises a manipulator arm, a mounting seat (S1), a vertical column (S2), an output gripper (S12) and a vertical cylinder (S23); the mounting seat is fixed on the mobile platform; 机械臂包括竖直升降移动台(S3)、臂段一(S6)、臂段二(S7)、臂段三(S8)和臂段四(S10);The mechanical arm comprises a vertical lifting and moving platform (S3), arm segment one (S6), arm segment two (S7), arm segment three (S8) and arm segment four (S10); 在竖直柱的左表面上设有竖直轨道(S24),竖直升降移动台上下竖直滑动设置在竖直轨道上;竖直柱的下端固定连接在安装座的上表面上,竖直气缸的缸座(S21)固定连接在位于竖直轨道左方的安装座的上表面上,竖直气缸的伸缩杆(S22)竖直朝上布置,竖直气缸的伸缩杆上端固定连接在竖直升降移动台的下表面上;竖直升降移动台能在竖直气缸的伸缩杆的带动下可沿着竖直轨道上下运动,形成第一自由度;A vertical track (S24) is provided on the left surface of the vertical column, and the vertical lifting and moving platform is vertically slidably arranged on the vertical track; the lower end of the vertical column is fixedly connected to the upper surface of the mounting seat, the cylinder seat (S21) of the vertical cylinder is fixedly connected to the upper surface of the mounting seat located on the left side of the vertical track, the telescopic rod (S22) of the vertical cylinder is arranged vertically upward, and the upper end of the telescopic rod of the vertical cylinder is fixedly connected to the lower surface of the vertical lifting and moving platform; the vertical lifting and moving platform can move up and down along the vertical track driven by the telescopic rod of the vertical cylinder, forming a first degree of freedom; 臂段一包括A1段管(S13)和伸缩连接在A1段管的左管口内的A2段管(S14),在A1段管内的右端固定设有伸缩杆水平朝左布置的一号气缸(S25),一号气缸的伸缩杆固定连接在A2段管的右端;The arm section 1 includes an A1 section pipe (S13) and an A2 section pipe (S14) telescopically connected to the left pipe opening of the A1 section pipe, and a No. 1 cylinder (S25) with a telescopic rod horizontally arranged to the left is fixedly provided at the right end of the A1 section pipe, and the telescopic rod of the No. 1 cylinder is fixedly connected to the right end of the A2 section pipe; 臂段二包括B1段管(S16)和伸缩连接在B1段管的左管口内的B2段管(S17),在B1段管内的右端固定设有伸缩杆水平朝左布置的二号气缸(37),二号气缸的伸缩杆固定连接在B2段管的右端;The arm section 2 comprises a B1 section pipe (S16) and a B2 section pipe (S17) telescopically connected to the left pipe opening of the B1 section pipe, a No. 2 air cylinder (37) with a telescopic rod horizontally arranged to the left is fixedly provided at the right end of the B1 section pipe, and the telescopic rod of the No. 2 air cylinder is fixedly connected to the right end of the B2 section pipe; 在竖直升降移动台的左端设有由一台一号减速电机(S26)驱动的一号水平旋转轴(S4),A1段管的右端固定连接在一号水平旋转轴上,使臂段一能水平转动形成第二自由度;在一号水平旋转轴上还设有能控制一号水平旋转轴转动的一号电磁刹车(S32);A No. 1 horizontal rotating shaft (S4) driven by a No. 1 reduction motor (S26) is provided at the left end of the vertical lifting and moving platform, and the right end of the A1 section pipe is fixedly connected to the No. 1 horizontal rotating shaft, so that the arm section 1 can rotate horizontally to form a second degree of freedom; and a No. 1 electromagnetic brake (S32) capable of controlling the rotation of the No. 1 horizontal rotating shaft is also provided on the No. 1 horizontal rotating shaft; 在A2段管的左端设有由一台二号减速电机(S27)驱动的二号水平旋转轴(S15),B1段管的右端固定连接在二号水平旋转轴上,使臂段二能水平转动形成第三自由度;在二号水平旋转轴上还设有能控制二号水平旋转轴转动的二号电磁刹车(S33);A second horizontal rotation axis (S15) driven by a second reduction motor (S27) is provided at the left end of the A2 section pipe, and the right end of the B1 section pipe is fixedly connected to the second horizontal rotation axis, so that the arm section 2 can rotate horizontally to form a third degree of freedom; a second electromagnetic brake (S33) capable of controlling the rotation of the second horizontal rotation axis is also provided on the second horizontal rotation axis; 在B2段管的左端设有由一台三号减速电机(S28)驱动的三号水平旋转轴(S18),臂段三的右端固定连接在三号水平旋转轴上,使臂段三能水平转动形成第四自由度;在三号水平旋转轴上还设有能控制三号水平旋转轴转动的三号电磁刹车(S34);A No. 3 horizontal rotation axis (S18) driven by a No. 3 reduction motor (S28) is provided at the left end of the B2 section tube, and the right end of the arm section 3 is fixedly connected to the No. 3 horizontal rotation axis, so that the arm section 3 can rotate horizontally to form a fourth degree of freedom; a No. 3 electromagnetic brake (S34) capable of controlling the rotation of the No. 3 horizontal rotation axis is also provided on the No. 3 horizontal rotation axis; 在臂段三的左端设有由一台四号减速电机(S29)驱动的能在左右竖直面上转动的一号横向竖直旋转轴(S9),臂段四的右端固定连接在一号横向竖直旋转轴上,使臂段四能在左右竖直面上竖直转动形成第五自由度;在一号横向竖直旋转轴上还设有能控制一号横向竖直旋转轴转动的四号电磁刹车(S35);A first horizontal vertical rotation axis (S9) driven by a fourth reduction motor (S29) and capable of rotating on a left and right vertical plane is provided at the left end of the arm segment 3, and the right end of the arm segment 4 is fixedly connected to the first horizontal vertical rotation axis, so that the arm segment 4 can rotate vertically on the left and right vertical planes to form a fifth degree of freedom; a fourth electromagnetic brake (S35) capable of controlling the rotation of the first horizontal vertical rotation axis is also provided on the first horizontal vertical rotation axis; 在臂段四的左端设有由一台五号减速电机(S30)驱动的能在前后竖直面上转动的一号纵向竖直旋转轴(S11),输出抓手的右端固定连接在一号纵向竖直旋转轴上,使输出抓手的右端能在前后竖直面上竖直转动形成第六自由度;在一号纵向竖直旋转轴上还设有能控制一号纵向竖直旋转轴转动的五号电磁刹车(S36);A No. 1 longitudinal vertical rotation axis (S11) driven by a No. 5 reduction motor (S30) and capable of rotating on the front and rear vertical planes is provided at the left end of the arm section 4. The right end of the output gripper is fixedly connected to the No. 1 longitudinal vertical rotation axis, so that the right end of the output gripper can rotate vertically on the front and rear vertical planes to form a sixth degree of freedom; a No. 5 electromagnetic brake (S36) capable of controlling the rotation of the No. 1 longitudinal vertical rotation axis is also provided on the No. 1 longitudinal vertical rotation axis; A2段管能在一号气缸的伸缩杆的带动下在A1段管内左右伸缩移动形成第七自由度;Driven by the telescopic rod of the No. 1 cylinder, the A2 section of pipe can move left and right in the A1 section of pipe to form the seventh degree of freedom; B2段管能在二号气缸的伸缩杆的带动下在B1段管内左右伸缩移动形成第八自由度;Driven by the telescopic rod of the No. 2 cylinder, the B2 section of pipe can move left and right in the B1 section of pipe to form the eighth degree of freedom; 一根一号水平管(S39)的左端水平固定连接在竖直柱的右表面上,在一号水平管内左右滑动设有平衡调节块(S40),在一号水平管内的左端固定连接有伸缩杆水平朝右的平衡调节气缸(S38),平衡调节气缸的伸缩杆的右端固定连接在平衡调节块上;The left end of a No. 1 horizontal tube (S39) is horizontally fixedly connected to the right surface of the vertical column, and a balance adjustment block (S40) is provided in the No. 1 horizontal tube for sliding left and right. A balance adjustment cylinder (S38) with a telescopic rod facing horizontally to the right is fixedly connected to the left end of the No. 1 horizontal tube, and the right end of the telescopic rod of the balance adjustment cylinder is fixedly connected to the balance adjustment block; 一号电磁刹车的控制端、二号电磁刹车的控制端、三号电磁刹车的控制端、四号电磁刹车的控制端、五号电磁刹车的控制端、一号减速电机的控制端、二号减速电机的控制端、三号减速电机的控制端、四号减速电机的控制端、五号减速电机的控制端、一号气缸的控制端、二号气缸的控制端、平衡调节气缸的控制端和竖直气缸的控制端分别控制连接在控制台上。The control end of electromagnetic brake No. 1, the control end of electromagnetic brake No. 2, the control end of electromagnetic brake No. 3, the control end of electromagnetic brake No. 4, the control end of electromagnetic brake No. 5, the control end of reduction motor No. 1, the control end of reduction motor No. 2, the control end of reduction motor No. 3, the control end of reduction motor No. 4, the control end of reduction motor No. 5, the control end of cylinder No. 1, the control end of cylinder No. 2, the control end of balance adjustment cylinder and the control end of vertical cylinder are respectively controlled and connected to the control console.
CN202010369760.5A 2020-04-30 2020-04-30 Inverse Priority Impedance Control System and Control Method for Redundant Manipulator Arm of Mobile Manipulator Active CN111687834B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010369760.5A CN111687834B (en) 2020-04-30 2020-04-30 Inverse Priority Impedance Control System and Control Method for Redundant Manipulator Arm of Mobile Manipulator

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010369760.5A CN111687834B (en) 2020-04-30 2020-04-30 Inverse Priority Impedance Control System and Control Method for Redundant Manipulator Arm of Mobile Manipulator

Publications (2)

Publication Number Publication Date
CN111687834A CN111687834A (en) 2020-09-22
CN111687834B true CN111687834B (en) 2023-06-02

Family

ID=72476984

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010369760.5A Active CN111687834B (en) 2020-04-30 2020-04-30 Inverse Priority Impedance Control System and Control Method for Redundant Manipulator Arm of Mobile Manipulator

Country Status (1)

Country Link
CN (1) CN111687834B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114800527A (en) * 2022-06-06 2022-07-29 山东大学 Force application control method and system for tail end of mobile operation mechanical arm

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4999553A (en) * 1989-12-28 1991-03-12 The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration Method and apparatus for configuration control of redundant robots
CN108247631A (en) * 2017-12-06 2018-07-06 西北工业大学 A kind of autonomous robust of mechanical arm for improving path trace performance keeps away unusual method
CN108839025A (en) * 2018-07-12 2018-11-20 杭州电子科技大学 A kind of motion planning method and device of mobile mechanical arm
CN110104216A (en) * 2019-01-28 2019-08-09 西北工业大学深圳研究院 A kind of collaboration paths planning method for kinematic redundancy dual-arm space robot

Family Cites Families (23)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS6416389A (en) * 1987-07-11 1989-01-19 Agency Ind Science Techn Control system of multi-joint type arm robot having redundancy
EP1854425A1 (en) * 2006-05-11 2007-11-14 BrainLAB AG Position determination for medical devices with redundant position measurement and weighting to prioritise measurements
JP4508164B2 (en) * 2006-06-26 2010-07-21 トヨタ自動車株式会社 Articulated robot and its control program
JP2009066685A (en) * 2007-09-11 2009-04-02 Sony Corp Robot device, and control method for robot device
US8170718B2 (en) * 2008-12-18 2012-05-01 GM Global Technology Operations LLC Multiple priority operational space impedance control
US8364314B2 (en) * 2009-04-30 2013-01-29 GM Global Technology Operations LLC Method and apparatus for automatic control of a humanoid robot
JP2011183527A (en) * 2010-03-10 2011-09-22 Toyota Motor Corp Joint target value determining method of redundant robot and control device of redundant robot
CN205363898U (en) * 2015-12-31 2016-07-06 中国科学院沈阳自动化研究所 Open -ended modularity arm
CN105676636B (en) * 2016-01-11 2018-12-07 北京邮电大学 A kind of redundancy space manipulator Multipurpose Optimal Method based on NSGA-II algorithm
CN107650120B (en) * 2016-07-26 2022-04-19 深圳华清精密科技有限公司 Method for determining all singular configurations of 9-degree-of-freedom mechanical arm
DE102016222255B3 (en) * 2016-11-14 2018-04-12 Kuka Roboter Gmbh Robotic arm, mobile robot and logistics system
CN106708069B (en) * 2017-01-19 2020-01-10 中国科学院自动化研究所 Coordinated planning and control method of underwater movement robot
CN107290956B (en) * 2017-08-01 2019-08-20 浙江大学 A Position Control Method of Single Joint Flexible Manipulator Based on State Feedback
CN107351085B (en) * 2017-08-21 2019-11-08 西北工业大学 A Space Robot Collision Avoidance Method Based on Multiple Control Points
CN108153153B (en) * 2017-12-19 2020-09-11 哈尔滨工程大学 A learning variable impedance control system and control method
CN108098786B (en) * 2017-12-19 2020-12-04 扬州大学 Endoscopic robotic arm for fusion reactors
CN108241339B (en) * 2017-12-27 2020-09-04 北京航空航天大学 Motion solving and configuration control method of humanoid mechanical arm
CN108621163A (en) * 2018-05-10 2018-10-09 同济大学 A kind of redundancy tow-armed robot cooperation control method towards remittance tenon technique
CN109079780B (en) * 2018-08-08 2020-11-10 北京理工大学 Task Hierarchical Optimal Control Method of Distributed Mobile Manipulator Based on Generalized Coordinates
CN209631457U (en) * 2018-12-17 2019-11-15 清研同创机器人(天津)有限公司 A kind of Omni-mobile formula robot spraying system
CN110497378A (en) * 2019-09-04 2019-11-26 上海海事大学 An automatic loading and unloading robot for parts processing
CN110524522B (en) * 2019-09-19 2024-03-26 浙江工业大学 Multi-degree-of-freedom redundant serial-parallel anthropomorphic mechanical arm
CN111538234B (en) * 2020-07-08 2020-10-09 深圳市优必选科技股份有限公司 Task hierarchical control method and device, robot and readable storage medium

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4999553A (en) * 1989-12-28 1991-03-12 The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration Method and apparatus for configuration control of redundant robots
CN108247631A (en) * 2017-12-06 2018-07-06 西北工业大学 A kind of autonomous robust of mechanical arm for improving path trace performance keeps away unusual method
CN108839025A (en) * 2018-07-12 2018-11-20 杭州电子科技大学 A kind of motion planning method and device of mobile mechanical arm
CN110104216A (en) * 2019-01-28 2019-08-09 西北工业大学深圳研究院 A kind of collaboration paths planning method for kinematic redundancy dual-arm space robot

Also Published As

Publication number Publication date
CN111687834A (en) 2020-09-22

Similar Documents

Publication Publication Date Title
CN110039542B (en) Visual servo tracking control method with speed and direction control function and robot system
Zhang et al. Integrated design of mechanical structure and control algorithm for a programmable four-bar linkage
Cheah et al. Feedback control for robotic manipulator with uncertain kinematics and dynamics
US20100152899A1 (en) Systems and methods of coordination control for robot manipulation
CN111687827B (en) Control method and control system for coordinating and operating weak rigid member by two robots
CN106493735B (en) There are the flexible mechanical arm disturbance observation control methods of external disturbance
Dong et al. A new image-based visual servoing method with velocity direction control
Colbaugh et al. Global regulation of uncertain manipulators using bounded controls
CN111687835B (en) System and method for controlling reverse priority impedance of redundant mechanical arm of underwater mechanical arm
DE102020120116B4 (en) Method for impedance-based multi-tasking tracking control, impedance-based multi-tasking tracking controller and force and/or torque controlled robot
Chung et al. Torque optimizing control with singularity-robustness for kinematically redundant robots
CN111687834B (en) Inverse Priority Impedance Control System and Control Method for Redundant Manipulator Arm of Mobile Manipulator
Liu et al. Research of driving force coordination mechanism in parallel manipulator with actuation redundancy and its performance evaluation
CN118024250A (en) Humanoid robot control method based on high-order full-drive theory
You et al. Semi-autonomous bilateral teleoperation of hexapod robot based on haptic force feedback
CN117532616A (en) Master-slave heterogeneous similarity mapping control method and system for redundant hydraulic mechanical arm
CN111687832B (en) Space Manipulator Redundant Manipulator Inverse Priority Impedance Control System and Control Method
CN111702753A (en) Redundant manipulator inverse priority impedance control system and control method
Li et al. Improved sliding mode control for mobile manipulators based on an adaptive neural network
CN111687833B (en) Manipulator anti-priority impedance control system and control method
Subudhi et al. On the singular perturbation approach to trajectory control of a multilink manipulator with flexible links and joints
CN117032281A (en) Six-foot robot foot-end nonsingular terminal and generalized impedance hybrid control method
Zhong et al. Theoretical and experimental study on remote dynamic balance control for a suspended wheeled mobile manipulator
Zhao et al. Nonlinear composite bilateral control framework for n-DOF teleoperation systems with disturbances
CN111687822B (en) Double-robot guiding track multi-space self-adaptive interpolation method

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
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20250609

Address after: 518000 1002, Building A, Zhiyun Industrial Park, No. 13, Huaxing Road, Henglang Community, Longhua District, Shenzhen, Guangdong Province

Patentee after: Shenzhen Wanzhida Technology Co.,Ltd.

Country or region after: China

Address before: 545006 Liuzhou City East District Road, central city, No. 268, No.

Patentee before: GUANGXI University OF SCIENCE AND TECHNOLOGY

Country or region before: China