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 PDFInfo
- 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
Links
Images
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B25—HAND TOOLS; PORTABLE POWER-DRIVEN TOOLS; MANIPULATORS
- B25J—MANIPULATORS; CHAMBERS PROVIDED WITH MANIPULATION DEVICES
- B25J9/00—Programme-controlled manipulators
- B25J9/16—Programme controls
- B25J9/1628—Programme controls characterised by the control loop
- B25J9/1633—Programme controls characterised by the control loop compliant, force, torque control, e.g. combined with position control
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B25—HAND TOOLS; PORTABLE POWER-DRIVEN TOOLS; MANIPULATORS
- B25J—MANIPULATORS; CHAMBERS PROVIDED WITH MANIPULATION DEVICES
- B25J9/00—Programme-controlled manipulators
- B25J9/16—Programme controls
- B25J9/1679—Programme controls characterised by the tasks executed
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/02—Total 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
Description
技术领域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
在竖直升降移动台的左端设有由一台一号减速电机驱动的一号水平旋转轴,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
在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
在臂段三的左端设有由一台四号减速电机驱动的能在左右竖直面上转动的一号横向竖直旋转轴,臂段四的右端固定连接在一号横向竖直旋转轴上,使臂段四能在左右竖直面上竖直转动形成第五自由度;在一号横向竖直旋转轴上还设有能控制一号横向竖直旋转轴转动的四号电磁刹车;At the left end of
在臂段四的左端设有由一台五号减速电机驱动的能在前后竖直面上转动的一号纵向竖直旋转轴,输出抓手的右端固定连接在一号纵向竖直旋转轴上,使输出抓手的右端能在前后竖直面上竖直转动形成第六自由度;在一号纵向竖直旋转轴上还设有能控制一号纵向竖直旋转轴转动的五号电磁刹车;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、建立奇异鲁棒解逆运动学分析模型;
步骤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、其关节空间的角位置、角速度分别为q、J为n自由度机器人的雅可比矩阵,其中x∈Rn, J∈Rm×n;冗余自由度机械臂的正运动学方程可用下式描述:Define the position and velocity of the end effector in Cartesian space as x, The angular position and angular velocity of its joint space are q, J is the Jacobian matrix of the n-DOF robot, where x∈R n , J∈R m×n ; The forward kinematics equation of the redundant degree of freedom manipulator can be described by the following formula:
式(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:
式(1)的解可通过寻找最佳来解决;The solution of formula (1) can be obtained by finding the optimal to solve;
因此,式(1)的伪逆解可表示为:Therefore, the pseudo-inverse solution of formula (1) can be expressed as:
式中J+——雅可比矩阵的伪逆Where J + ——pseudo-inverse of Jacobian matrix
I——单位矩阵I——Identity Matrix
——任意零空间矢量 ——Arbitrary null space vector
——最小范数解,定义了手的动作 ——Minimum norm solution, defining the hand movements
——齐次解,在末端不产生动作 ——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,
其中λ≥0是加权矩阵,为加权系数,并且满足where λ≥0 is the weighting matrix, is the weighting coefficient and satisfies
上述方程的解就可表示为:The solution to the above equation can be expressed as:
式(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:
作为优选,建立通过奇异鲁棒解得出消除奇异性算法的任务优先级解决策略如下: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:
其中和表示task1和task2in and Represents task1 and task2
从表达式(9)中得出冗余机械臂的逆运动学方程为:From expression (9), the inverse kinematics equation of the redundant manipulator is obtained as follows:
Task1作为主任务,task2作为辅助任务;也就是说,task2是在task1的零空间中实现的;冗余机械臂的最终逆运动学表达式如下所示:Task1 is the main task and task2 is the auxiliary task; that is, task2 It is in task1 The final inverse kinematics expression of the redundant manipulator is as follows:
其中 in
是投影矩阵,它给出了次要任务到主要任务的适用范围;和是所需的命令速度;是主要任务,是次要任务; is the projection matrix, which gives the applicability of the secondary task to the primary task; and is the desired command speed; is the main task, 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
其中ρ(·)是矩阵的秩;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:
作为优选,建立奇异鲁棒解逆运动学分析模型如下: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:
因此,上述方程的奇异鲁棒伪逆解可表示为:Therefore, the singular robust pseudo-inverse solution of the above equation can be expressed as:
式(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;
其中,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:
因子λ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
从上式中,我们可看出参数δ>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;
引入逆优先投影矩阵该矩阵包括了不依赖于第k个任务的最低优先级l-k-1任务的相应元素的零空间,所以得出Introducing the inverse priority projection matrix 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
其中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:
在上述推导中,k=1,l-1,…,1,其中1为正整数;初始值 In the above derivation, k = 1, l-1 , ..., 1, where 1 is a positive integer; the initial value
为了给出线性无关雅可比矩阵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:
就有了There you have it
其中表示的行;in express 's line;
于是乎,的伪逆解就可表示为:So, The pseudo-inverse solution of can be expressed as:
和 and
其中,Tk表示矩阵的扩充;Where T k represents the matrix expansion;
最终的反向优先级投影可写为:The final reverse priority projection can be written as:
于是,我们就可得出伪逆解的表达式:Therefore, we can derive the expression of the pseudo-inverse solution:
建立多任务冗余机械臂的逆优先控制策略方程如下:The inverse priority control strategy equation for the multi-task redundant manipulator is established as follows:
作为优选,对具有主要任务和次要任务的冗余机械臂的反向控制方程进行简化处理如下: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
上面的公式与前面的表达式(11)有很大的不同,但是算法框架是相似的;在上面的方程中,是次要任务,是主任务;主任务在主任务的指定零空间中实现;反向优先级的核心点是投影矩阵的计算;的表达式如公式(30):The above formula is quite different from the previous expression (11), but the algorithm framework is similar; in the above equation, It's a secondary task. 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 Calculation of The expression of is as follows:
利用前述公式(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:
作为优选,建立机械手的逆优先力控制策略如下: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:
其中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:
其中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
上述两个方程的组合如下The combination of the above two equations is as follows
从上面的方程可看出,如果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
所以这两个方程的积分公式可写成So the integral formula of these two equations can be written as
如果机械手末端执行器能够跟踪期望的笛卡尔速度为和则可实现机械手的精确力控制;笛卡尔速度与关节速度的关系应借鉴逆优先控制;因此,可得到机械手的逆优先力控制策略的方程:If the robot end effector can track the desired Cartesian velocity and 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:
上述方程所要求的关节速度将保证机械手的力控制;值得一提的是,上述力控制律只是速度级控制律,它依赖于内速度环控制;如果内位置控制效果良好,则可实现精确的力控制;由于内速度环控制可实现低频位置跟踪,所以外力环可实现低频力跟踪。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
参考速度可表示为The reference speed can be expressed as
因此,机械手的逆优先阻抗控制保证的表达式为:Therefore, the expression guaranteed by the reverse priority impedance control of the manipulator is:
作为优选,将位置控制空间的逆优先计算扩展到力控制空间的逆优先计算,从而获得机械手速度级逆优先阻抗控制的总体框架实现方式如下: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
所以期望速度的简化形式可表示为So the simplified form of the expected speed can be expressed as
然后我们得到了基于反向优先级的解决方案Then we get the solution based on reverse priority
考虑到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:
表达式(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
在竖直升降移动台的左端设有由一台一号减速电机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
在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
在臂段三的左端设有由一台四号减速电机S29驱动的能在左右竖直面上转动的一号横向竖直旋转轴S9,臂段四的右端固定连接在一号横向竖直旋转轴上,使臂段四能在左右竖直面上竖直转动形成第五自由度;在一号横向竖直旋转轴上还设有能控制一号横向竖直旋转轴转动的四号电磁刹车S35;At the left end of
在臂段四的左端设有由一台五号减速电机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、其关节空间的角位置、角速度分别为q、J为n自由度机器人的雅可比矩阵,其中x∈Rn, J∈Rm×n;冗余自由度机械臂的正运动学方程可用下式描述:Define the position and velocity of the end effector in Cartesian space as x, The angular position and angular velocity of its joint space are q, J is the Jacobian matrix of the n-DOF robot, where x∈R n , J∈R m×n ; The forward kinematics equation of the redundant degree of freedom manipulator can be described by the following formula:
式(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:
式(1)的解可通过寻找最佳来解决;The solution of formula (1) can be obtained by finding the optimal to solve;
因此,式(1)的伪逆解可表示为:Therefore, the pseudo-inverse solution of formula (1) can be expressed as:
式中J+——雅可比矩阵的伪逆Where J + ——pseudo-inverse of Jacobian matrix
I——单位矩阵I——Identity Matrix
——任意零空间矢量 ——Arbitrary null space vector
——最小范数解,定义了手的动作 ——Minimum norm solution, defining the hand movements
一一齐次解,在末端不产生动作 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,
其中λ≥0是加权矩阵,为加权系数,并且满足where λ≥0 is the weighting matrix, is the weighting coefficient and satisfies
上述方程的解就可表示为:The solution to the above equation can be expressed as:
式(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:
步骤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:
其中和表示task1和task2in and Represents task1 and task2
从表达式(9)中得出冗余机械臂的逆运动学方程为:From expression (9), the inverse kinematics equation of the redundant manipulator is obtained as follows:
Task1作为主任务,task2作为辅助任务;也就是说,是在的零空间中实现的;冗余机械臂的最终逆运动学表达式如下所示:Task1 is the main task and task2 is the auxiliary task; that is, is The final inverse kinematics expression of the redundant manipulator is as follows:
其中 in
是投影矩阵,它给出了次要任务到主要任务的适用范围;和是所需的命令速度;是主要任务,是次要任务; is the projection matrix, which gives the applicability of the secondary task to the primary task; and is the desired command speed; is the main task, 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
其中ρ(·)是矩阵的秩;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:
步骤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:
因此,上述方程的奇异鲁棒伪逆解可表示为:Therefore, the singular robust pseudo-inverse solution of the above equation can be expressed as:
式(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;
其中,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:
因子λ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
从上式中,我们可看出参数δ>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;
引入逆优先投影矩阵该矩阵包括了不依赖于第k个任务的最低优先级l-k-1任务的相应元素的零空间,所以得出Introducing the inverse priority projection matrix 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
其中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:
在上述推导中,k=1,l-1,…,1,其中l为正整数;初始值 In the above derivation, k = 1, l-1, ..., 1, where l is a positive integer; the initial value
为了给出线性无关雅可比矩阵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:
就有了There you have it
其中表示的行;in express 's line;
于是乎,的伪逆解就可表示为:So, The pseudo-inverse solution of can be expressed as:
和 and
其中,TR表示矩阵的扩充;Among them, TR represents the matrix expansion;
最终的反向优先级投影可写为:The final reverse priority projection can be written as:
于是,我们就可得出伪逆解的表达式:Therefore, we can derive the expression of the pseudo-inverse solution:
建立多任务冗余机械臂的逆优先控制策略方程如下:The inverse priority control strategy equation for the multi-task redundant manipulator is established as follows:
步骤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
上面的公式与前面的表达式(11)有很大的不同,但是算法框架是相似的;在上面的方程中,是次要任务,是主任务;主任务在主任务的指定零空间中实现;反向优先级的核心点是投影矩阵的计算;的表达式如公式(30):The above formula is quite different from the previous expression (11), but the algorithm framework is similar; in the above equation, It's a secondary task. 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 Calculation of The expression of is as follows:
利用前述公式(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:
步骤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:
其中X是笛卡尔空间中的位置,M(X是惯性矩阵,是非线性力,F是输入控制力,Fe是接触力;where X is the position in Cartesian space, M (X is the inertia matrix, 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:
其中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
上述两个方程的组合如下The combination of the above two equations is as follows
从上面的方程可看出,如果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
所以这两个方程的积分公式可写成So the integral formula of these two equations can be written as
如果机械手末端执行器能够跟踪期望的笛卡尔速度为和则可实现机械手的精确力控制;笛卡尔速度与关节速度的关系应借鉴逆优先控制;因此,可得到机械手的逆优先力控制策略的方程:If the robot end effector can track the desired Cartesian velocity and 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:
上述方程所要求的关节速度将保证机械手的力控制;值得一提的是,上述力控制律只是速度级控制律,它依赖于内速度环控制;如果内位置控制效果良好,则可实现精确的力控制;由于内速度环控制可实现低频位置跟踪,所以外力环可实现低频力跟踪。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
参考速度可表示为The reference speed can be expressed as
因此,机械手的逆优先阻抗控制保证的表达式为:Therefore, the expression guaranteed by the reverse priority impedance control of the manipulator is:
步骤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
所以期望速度的简化形式可表示为So the simplified form of the expected speed can be expressed as
然后我们得到了基于反向优先级的解决方案Then we get the solution based on reverse priority
混合阻抗控制的动力学方案图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:
表达式(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)
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)
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)
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)
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 |
-
2020
- 2020-04-30 CN CN202010369760.5A patent/CN111687834B/en active Active
Patent Citations (4)
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 |