WO2018107851A1 - Method and device for controlling redundant robot arm - Google Patents

Method and device for controlling redundant robot arm Download PDF

Info

Publication number
WO2018107851A1
WO2018107851A1 PCT/CN2017/103286 CN2017103286W WO2018107851A1 WO 2018107851 A1 WO2018107851 A1 WO 2018107851A1 CN 2017103286 W CN2017103286 W CN 2017103286W WO 2018107851 A1 WO2018107851 A1 WO 2018107851A1
Authority
WO
WIPO (PCT)
Prior art keywords
redundant
trajectory
robot arm
mechanical arm
joint
Prior art date
Application number
PCT/CN2017/103286
Other languages
French (fr)
Chinese (zh)
Inventor
阳方平
Original Assignee
广州视源电子科技股份有限公司
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 广州视源电子科技股份有限公司 filed Critical 广州视源电子科技股份有限公司
Publication of WO2018107851A1 publication Critical patent/WO2018107851A1/en

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
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B25HAND TOOLS; PORTABLE POWER-DRIVEN TOOLS; MANIPULATORS
    • B25JMANIPULATORS; CHAMBERS PROVIDED WITH MANIPULATION DEVICES
    • B25J9/00Programme-controlled manipulators
    • B25J9/16Programme controls
    • B25J9/1656Programme controls characterised by programming, planning systems for manipulators
    • B25J9/1664Programme controls characterised by programming, planning systems for manipulators characterised by motion, path, trajectory planning

Definitions

  • the present invention relates to the field of robot technology, and in particular, to a control method and apparatus for a redundant mechanical arm.
  • the redundant manipulator When the redundant manipulator performs a Cartesian space-specific task, there is an infinite group solution in the joint space.
  • a set of optimal solutions can be selected according to the optimization index such as joint speed, joint torque or obstacle distance to determine the redundant machine. The movement of the arm.
  • the methods commonly used for multi-objective optimization can be classified into a weighted function method and a Pareto frontier method.
  • the weighting function method can use the weighting coefficient to weight the optimization target to transform the multi-objective optimization problem into a single-objective optimization problem, it is difficult to judge the influence of the weighting coefficient change on the optimization result.
  • an embodiment of the present invention provides a control method and apparatus for a redundant mechanical arm to solve
  • the multi-objective solving method in the prior art calculates a complicated technical problem with complicated steps.
  • an embodiment of the present invention provides a method for controlling a redundant mechanical arm, including:
  • the equation corresponding to the trajectory function is solved according to the target approaching method, and the position and velocity of each joint in the redundant mechanical arm corresponding to the motion trajectory are obtained.
  • an embodiment of the present invention further provides a control device for a redundant mechanical arm, including:
  • An information acquiring unit configured to acquire current point information and target point information where the redundant robot arm is located;
  • a trajectory unit configured to be connected to the information acquiring unit, configured to determine, according to the current point information and the target point information, a trajectory function corresponding to a motion trajectory of the redundant mechanical arm moving from a current point to a target point;
  • a redundancy processing unit connected to the trajectory unit, configured to establish an equation corresponding to the trajectory function by using a redundancy space vector as an independent variable; and solving an equation corresponding to the trajectory function according to a target approach method, to obtain the motion The position and velocity of each joint in the redundant robot arm corresponding to the trajectory.
  • the technical solution for controlling the redundant mechanical arm obtained by the embodiment of the present invention obtains current point information and target point information of the redundant mechanical arm, and determines that the redundant mechanical arm moves from the current point to the current point information and the target point information.
  • the trajectory function corresponding to the motion trajectory of the target point is constructed by using the redundancy space vector as the independent variable, and the equation is solved by the target approach method to obtain the position of each joint in the redundant mechanical arm corresponding to the motion trajectory. speed.
  • FIG. 1 is a schematic flow chart of a method for controlling a redundant mechanical arm according to Embodiment 1 of the present invention
  • FIG. 2 is a schematic flow chart of a method for controlling a redundant mechanical arm according to Embodiment 2 of the present invention
  • FIG. 3 is a schematic diagram of a coordinate system of a nine-degree-of-freedom redundant robot arm according to Embodiment 2 of the present invention.
  • FIG. 4 is a structural block diagram of a control device for a redundant mechanical arm according to Embodiment 3 of the present invention.
  • Embodiment 1 of the present invention provides a control method for a redundant mechanical arm.
  • the method can be performed by a control device of a redundant robot arm, wherein the device can be implemented by hardware and/or software, and can generally be integrated in a control module for controlling the robot.
  • 1 is a schematic flow chart of a method for controlling a redundant mechanical arm according to Embodiment 1 of the present invention. As shown in FIG. 1, the method includes:
  • the current point information and the target point information of the redundant robot arm may be redundant mechanical arm ends.
  • the current point information and the target point information of the terminal may also be the current point information and the target point information of the joints of the redundant robot arm, which are not limited herein.
  • the current point information of the redundant robot arm may include current point information of the end of the arm and current point information of each joint of the redundant manipulator, target point information Target point information at the end of the robot arm can be included.
  • the current point information (target point information) may include information such as position coordinates, angular coordinates, linear velocity, angular velocity, linear acceleration, and/or angular acceleration of the end of the arm or the joint of the robot arm at the current (target) position.
  • the current point information of the redundant manipulator can be obtained directly by the sensor or by using the calculation result of the last moment as the current point information of the redundant manipulator;
  • the target point information of the redundant manipulator can be calculated according to the current point information of the end of the redundant manipulator and the motion track information of the end of the redundant manipulator.
  • the current point information as the angular coordinate, angular velocity and angular acceleration as an example
  • the angular information and angular velocity information of a joint of the redundant mechanical arm at the current position can be obtained by an encoder installed at the joint position, and the angular acceleration information can be differentiated. The operation is obtained; or, the target point information of each joint of the robot arm calculated at the previous time can be directly used as the current point information of each joint of the robot arm at the current time.
  • the trajectory function corresponding to the motion trajectory of the redundant mechanical arm moving from the current point to the target point may include a trajectory corresponding to each joint motion trajectory of the redundant mechanical arm when the end of the redundant mechanical arm moves from the current point to the target point.
  • the trajectory function of each joint of the redundant manipulator can be one or more, that is, the trajectory function corresponding to each joint motion trajectory of the redundant manipulator can be described as a whole by a total trajectory function, or can be used for each joint
  • the motion trajectory is described in the form of one or several trajectory functions.
  • the trajectory function corresponding to the motion trajectory of each joint of the redundant manipulator may include an equation and/or an inequality.
  • the trajectory function of a joint of the redundant manipulator may include the position coordinates of the joint
  • the correspondence between speed and/or acceleration and other variables may also include coordinates of the joint position, moving speed, and/or moving angular velocity when the redundant manipulator moves normally.
  • the range of values, etc., is not limited here.
  • the trajectory function of the end of the redundant mechanical arm may be first determined according to the current point information and the target point information at the end of the redundant mechanical arm. Then, according to the connection relationship between the end of the redundant manipulator and each joint, the relationship between the motion path of each joint and the motion track of the end of the redundant manipulator is determined, and then the redundant robot arm is determined according to the trajectory function of the end of the redundant manipulator. Trajectory function.
  • the correspondence between the redundancy vector and the coordinates of each joint position, speed, or acceleration may be first determined, and then the position vector, speed, or position in the trajectory function corresponding to each joint is replaced by the redundancy space vector according to the correspondence. Acceleration and other independent variables, and then establish the equation corresponding to the trajectory function of each joint of the redundant manipulator with the redundancy space vector as the independent variable.
  • the redundancy space vector may be any vector of the redundancy vector space, and the redundancy vector space may be an existing vector space or a custom space, which is not limited herein.
  • the acceleration of each joint of the redundant mechanical arm can be further determined as needed.
  • all solutions of the equation can be obtained, and then a set of solutions of the equations can be obtained according to actual needs or according to the set selection rules.
  • the target position of the joint at the next moment At the same time, the moving speed and acceleration of the joint when moving from the current position to the target position are determined.
  • a set of solutions of the equation can be arbitrarily obtained and the target position of the joint at the next moment is obtained by the set of solutions and moved from the current position to The moving speed and acceleration of the target position; first, the constraint condition of the equation solution can be set first, then the solution of the equation that satisfies the constraint condition is obtained, and the movement process of the joint from the current time to the next moment is guided by the set of solutions, There are no restrictions.
  • the control method of the redundant mechanical arm acquires current point information and target point information of the redundant mechanical arm, and determines that the redundant mechanical arm moves from the current point to the current point information and the target point information.
  • the trajectory function corresponding to the motion trajectory of the target point is constructed by using the redundancy space vector as the independent variable, and the equation is solved by the target approach method to obtain the position of each joint in the redundant mechanical arm corresponding to the motion trajectory. speed.
  • the present embodiment can reduce the search space of the multi-objective optimization problem of the redundant manipulator under the premise of satisfying the global optimal solution of the design target, and avoid the occurrence of the dimensional explosion problem in the multi-objective solution process, simplifying The amount of calculation required during the redundant manipulator control process increases the reaction speed of the redundant manipulator.
  • the acquiring the current point information and the target point information where the redundant mechanical arm is located includes: acquiring the speed of the end of the redundant mechanical arm at the target point; According to the speed of the end of the redundant manipulator at the target current point, the joint of the redundant manipulator The angular, joint angular velocity and joint angular acceleration are mapped to the redundancy space.
  • the determining, according to the current point information and the target point information, a trajectory function corresponding to a motion trajectory of the redundant robot arm moving from a current point to a target point comprises: determining the redundancy according to an inverse kinematics equation The optimization objective function and constraints corresponding to the motion trajectory of the remaining manipulator moving from the current point to the target point.
  • the solving the equation corresponding to the trajectory function according to the target approaching method, and obtaining the position and velocity of each joint in the redundant mechanical arm corresponding to the motion trajectory includes: in the redundancy space, according to the optimization target a function, a minimum target value corresponding to the optimization target, a corresponding constraint condition, and a preset weighting coefficient, and an equation corresponding to the trajectory function is solved by an auxiliary vector and a single target optimization algorithm to obtain the redundancy corresponding to the motion trajectory The position and speed of the joints in the robot arm.
  • control method of the redundant mechanical arm includes:
  • the obtaining the speed of the end of the redundant mechanical arm at the current point comprises: according to the posture of the redundant mechanical arm end at the current point and the movement of the redundant mechanical arm end
  • the pose of the next sample point in the trajectory determines the velocity of the end of the redundant robot arm at the next sample point.
  • the position of the end of the redundant manipulator refers to the position and posture of the redundant manipulator.
  • the position can be represented by position coordinates, and the posture can be expressed by angle.
  • information such as the angular velocity and/or angular acceleration of the redundant manipulator at the current point can be obtained from an encoder installed at the end of the redundant manipulator, or directly based on the calculation result of the redundant manipulator at the previous moment.
  • the relationship between the end speed of the redundant manipulator and the joint angular velocity of the redundant manipulator ie, the trajectory function of each joint of the redundant manipulator with angular velocity as an independent variable
  • N the degree of freedom of the redundant manipulator
  • N 6+r, r>1
  • J ⁇ R 6 ⁇ N is the Jacobian matrix.
  • equation (1) can be rewritten as:
  • Redundant manipulator joint angular acceleration The expression for the redundancy space vector k is:
  • ⁇ 0 is the joint angle of the previous moment
  • ⁇ t is the sampling step size
  • the optimization objective function and the constraint condition corresponding to the motion trajectory of the redundant mechanical arm moving from the current point to the target point can be flexibly set according to requirements, for example, the range of the redundant mechanical arm movement can be restricted to avoid Obstacle, or limit the bending angle of the redundant manipulator to avoid the joints reaching their motion limits.
  • the optimization objective function includes a minimum joint motion amplitude of the redundant mechanical arm or the redundant mechanical arm is minimally affected by a corresponding gear gap;
  • the constraint condition includes : constraining the joint speed of the redundant mechanical arm within a preset speed range or constraining the joint angle of the redundant mechanical arm within a preset angular range.
  • the optimal objective function for the minimum joint motion of the redundant manipulator can be:
  • the optimal objective function for the redundant manipulator to be minimally affected by the corresponding gear clearance can be:
  • the constraint that constrains the joint speed of the redundant manipulator to a preset speed range may be:
  • the constraint that constrains the joint angle of the redundant manipulator within a predetermined range of angles may be: ⁇ min ⁇ ⁇ ⁇ ⁇ max .
  • the optimization objective function and the constraint equation of the redundant robot arm trajectory can be established based on the redundancy space, and the optimization objective function and the constraint equation of the redundant mechanical arm can be expressed as joint angle, joint angular velocity and/or
  • the joint angular acceleration is an expression of the independent variable. Therefore, by substituting the above equations (7), (8), and (9) into the expression, the optimal objective function and the constraint equation with the redundant space vector k as the independent variable can be obtained.
  • an optimization objective function with minimal joint motion amplitude can be further expressed as an equation for the redundancy space vector k:
  • the optimization objective function with respect to the minimum influence of the mechanical arm on the corresponding gear gap can be further expressed as the equation of the redundancy space vector k:
  • the constraint on the joint angular velocity limit can be further expressed as an equation with the redundant space vector k as an independent variable:
  • the constraint on the joint angle limit can be further expressed as an equation with the redundancy space vector k as an independent variable:
  • the single-objective optimization algorithm for solving the equation corresponding to the trajectory function can be flexibly determined according to needs, and is not limited herein.
  • the single target optimization algorithm comprises: a Newton-Eulerian algorithm, a Nelder-Mead simplex algorithm, an interior point algorithm, a genetic algorithm and/or a Pattern Search algorithm.
  • the Pareto frontier is a solution set for multi-objective optimization problems, which is defined as follows: On the multi-objective feasible domain ⁇ , there is a point x * ⁇ ⁇ , if there is no point x ⁇ ⁇ , so that f(x) ⁇ f(x * ) is established, and there is at least a point x' ⁇ ⁇ such that the strict inequality f(x') > f(x * ) holds, then x * is a Pareto optimal solution.
  • the equation corresponding to the trajectory function is as follows:
  • A is a matrix constraint condition
  • b is a vector constraint condition
  • is a weighting coefficient vector
  • is an auxiliary vector.
  • the weighting coefficient vector ⁇ can be set by the user or the developer as needed, and the auxiliary vector ⁇ is used to convert the original multi-objective optimization into a single-objective optimization problem.
  • the inner point method can be used to solve the equation (12) to obtain the optimal solution of the equation, and then the joints of the redundant mechanical arm and the optimal solution are obtained by the equations (7), (8) and (9). Corresponding angles, angular velocities, and angular angular velocities, resulting in the position and velocity of the joints of the redundant manipulator at the next moment.
  • the wrist coordinate system can be defined according to the structural characteristics of the redundant mechanical arm.
  • the specific coordinates and direction of the wrist coordinate system can be flexibly set according to requirements, which is not limited herein.
  • the wrist coordinate system can be defined as a coordinate system parallel to the coordinate end system of the redundant manipulator and the origin coincides with the origin of the coordinate system of the connecting rod 8.
  • 3 ⁇ w is the angular velocity in the wrist coordinate system with reference to the coordinate system of the connecting rod 3;
  • 3 ⁇ h is the angular velocity in the coordinate system of the redundant manipulator referenced by the coordinate system of the connecting rod 3;
  • 3 v w The linear velocity in the wrist coordinate system with reference to the coordinate system of the connecting rod 3;
  • 3 v h is the linear velocity in the coordinate system of the redundant manipulator referenced by the coordinate system of the connecting rod 3;
  • 3 x h is the connection
  • the rod 3 coordinate system is the value of the unit vector of the Z direction in the redundant robot arm end coordinate system of the reference:
  • c i represents cos( ⁇ i )
  • s i represents sin( ⁇ i )
  • c ij represents cos( ⁇ i + ⁇ j )
  • s ij represents sin( ⁇ i + ⁇ j ).
  • J w is the joint angular velocity Speed with wrist
  • the Jacques matrix has:
  • the Jacobian matrix can be found by the vector product method as follows:
  • any three columns in the equation (19) may be taken as ⁇ , and J * is the remaining columns other than the three columns constituting ⁇ .
  • J 5 , J 6 and J 7 are taken as ⁇ , and the remaining columns (J 1 , J 2 , J 3 , J 4 , J 8 and J 9 ) constitute J * , and can be obtained by the formula (5):
  • J 21 a 2 s 3 - d 6 c 5 - d 4 - a 5 s 5 - d 7 c 5 c 7 + d 7 s 5 c 6 s 7 ,
  • J 23 a 2 c 3 + a 4 c 4 + a 5 c 4 c 5 - d 6 c 4 s 5 - d 7 c 4 s 5 c 7 + d 7 s 4 s 6 s 7 -d 7 c 4 c 5 c 6 s 7 ,
  • J 31 -d 4 -d 7 (c 5 c 7 -s 5 c 6 s 7 )-d 6 c 5 -a 5 s 5 ,
  • J 33 a 4 c 4 + d 7 (s 7 (s 4 s 6 -c 4 c 5 c 6 )-c 4 s 5 c 7 )+a 5 c 4 c 5 -d 6 c 4 s 3 ,
  • J 41 d 7 (s 7 (c 4 s 6 +s 4 c 5 c 6 )+s 4 s 5 c 7 )-a 4 s 4 -a 5 s 4 c 5 +d 6 s 4 s 5 ,
  • J 42 d 7 (s 7 (s 4 s 6 -c 4 c 5 c 6 )-c 4 s 5 c 7 )+a 4 c 4 +a 5 c 4 c 5 -d 6 c 4 s 5 ,
  • J 84 s 7 (s 4 s 6 -c 4 c 5 c 6) -c 4 s 5 c 7,
  • J 85 -s 7 (c 4 s 6 + s 4 c 5 c 6) -s 4 s 5 c 7,
  • J 94 c 8 (s 4 c 6 +c 4 c 5 s 6 )-s 8 (c 7 (s 4 s 6 -c 4 c 5 c 6 ))+c 4 s 5 s 7 ,
  • J 95 s 8 (c 7 (s 4 s 6 +s 4 c 5 c 6 ))-s 4 s 5 s 7 -c 8 (c 4 c 6 -s 4 c 5 s 6 ),
  • J 96 s 8 (c 5 s 7 - s 5 c 6 c 7 ) + c 8 s 5 s 6 .
  • Matrix The ith element.
  • the joint angle ⁇ can be obtained as a function of the redundancy vector (k 5 , k 6 , k 7 ) as an independent variable:
  • ⁇ 0 is the current joint angle and ⁇ t is the sampling step size.
  • H i ( ⁇ (k 5 , k 6 , k 7 )) is the element corresponding to the second row and the fourth column of the link transformation matrix 0 T i ( ⁇ (k 5 , k 6 , k 7 )):
  • d i , ⁇ i and a i are Denavit-Hartenberg parameters, and specific values thereof are shown in Table 1.
  • the control method of the redundant mechanical arm obtaineds the speed of the end of the redundant mechanical arm at the target point, and the joint angle, the joint angular velocity of the redundant mechanical arm and the speed of the end of the redundant mechanical arm at the target point
  • the joint angular acceleration is mapped to the redundancy space.
  • the optimal objective function and constraints corresponding to the motion trajectory of the redundant manipulator moving from the current point to the target point are determined, and the redundant space vector is used as the independent variable to establish the motion.
  • the equation corresponding to the trajectory is solved by the auxiliary vector and the single-objective optimization algorithm to solve the equation corresponding to the trajectory function, and the position and velocity of each joint in the redundant mechanical arm corresponding to the motion trajectory are obtained.
  • a mathematical module for multi-objective constrained optimization problem is established, which can provide a theoretical basis for solving a multi-objective optimization problem, and reduce redundant manipulators under the premise of satisfying the global optimal solution of the design target.
  • the search space of the target optimization problem avoids the occurrence of dimensional explosion problems in the multi-objective solution process, simplifies the calculation amount required in the redundant robot arm control process, and improves the reaction speed of the redundant mechanical arm.
  • Embodiment 3 of the present invention provides a control device for a redundant mechanical arm.
  • the device can be implemented by hardware and/or software, generally integrated in a control module for controlling the robot, and the control of the redundant manipulator can be realized by performing a control method of the redundant manipulator.
  • 4 is a structural block diagram of a control device for a redundant mechanical arm according to the embodiment. As shown in FIG. 4, the device includes:
  • the information acquiring unit 410 is configured to acquire current point information and target point information where the redundant robot arm is located;
  • the trajectory unit 420 is connected to the information acquiring unit, and configured to determine, according to the current point information and the target point information, a trajectory function corresponding to a motion trajectory of the redundant mechanical arm moving from a current point to a target point;
  • the redundancy processing unit 430 is connected to the trajectory unit, and is configured to establish an equation corresponding to the trajectory function by using a redundancy space vector as an independent variable; and solving an equation corresponding to the trajectory function according to a target proximity method, to obtain the The position and velocity of each joint in the redundant robot arm corresponding to the motion trajectory.
  • the control device of the redundant mechanical arm obtained by the third embodiment of the present invention obtains the current point information and the target point information of the redundant mechanical arm through the information acquiring unit, and determines the redundancy according to the acquired current point information and the target point information by the trajectory unit.
  • the trajectory function corresponding to the motion trajectory of the robot arm moving from the current point to the target point is established by the redundancy processing unit with the redundancy space vector as the independent variable, and the equation is solved by the target approach method to obtain the motion.
  • the position and velocity of each joint in the redundant robot arm corresponding to the trajectory.
  • the embodiment of the present invention can reduce the search space of the multi-objective optimization problem of the redundant manipulator under the premise of obtaining the global optimal solution satisfying the design goal, and avoid the occurrence of the dimensional explosion problem in the multi-objective solution process. Simplify the amount of calculations required during redundant manipulator control and increase the response speed of redundant manipulators.
  • the information acquiring unit 410 is specifically configured to: acquire a speed of the end of the redundant mechanical arm at the target point; and the redundant machine according to a speed of the redundant mechanical arm end at the target point
  • the joint angle of the arm, the joint angular velocity, and the joint angular acceleration are mapped to the redundancy space.
  • the acquiring the speed of the end of the redundant mechanical arm at the current point comprises: according to the posture of the redundant mechanical arm end at the current point and the end of the redundant mechanical arm The pose of the next sample point in the trajectory determines the velocity of the end of the redundant robot arm at the next sample point.
  • trajectory unit 420 is specifically configured to determine an optimization objective function and a constraint condition corresponding to a motion trajectory of the redundant mechanical arm moving from a current point to a target point according to an inverse kinematics equation.
  • the optimization objective function includes a minimum joint motion amplitude of the redundant mechanical arm or the redundant mechanical arm is minimally affected by a corresponding gear gap;
  • the constraint condition includes: a joint of the redundant mechanical arm The speed constraint is within a preset speed range or the joint angle of the redundant robot arm is constrained within a preset angle range.
  • the solving the equation corresponding to the trajectory function according to the target approaching method, and obtaining the position and velocity of each joint in the redundant mechanical arm corresponding to the motion trajectory includes: in the redundancy space, according to the optimization target a function, a minimum target value corresponding to the optimization target, a corresponding constraint condition, and a preset weighting coefficient, and an equation corresponding to the trajectory function is solved by an auxiliary vector and a single target optimization algorithm to obtain the redundancy corresponding to the motion trajectory The position and speed of the joints in the robot arm.
  • the single target optimization algorithm includes: a Newton-Eulerian algorithm, a Nelder-Mead simplex algorithm, an interior point algorithm, a genetic algorithm, and/or a Pattern Search algorithm.
  • the control device of the redundant mechanical arm provided by this embodiment can execute the control method of the redundant mechanical arm provided by any embodiment of the present invention, and has the corresponding functional modules and beneficial effects of the control method for executing the redundant mechanical arm.
  • control method of the redundant robot arm provided by any embodiment of the present invention.

Landscapes

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

Abstract

A method and device for controlling a redundant robot arm. The method comprises: acquiring information about a current point where a redundant robot arm is located and information about a target point; according to the current point information and the target point information, determining a trajectory function corresponding to a motion trajectory of the redundant robot arm moving from the current point to the target point; establishing an equation corresponding to the trajectory function by using a redundancy spatial vector as an independent variable; and solving the equation corresponding to the trajectory function according to an objective approaching method, so as to obtain the position and speed of each joint in the redundant robot arm corresponding to the motion trajectory. By means of the technical solution, on the premise of obtaining an optimal global solution satisfying a design goal, the search space for a multi-objective optimization problem of a redundant robot arm can be reduced, the occurrence of dimension explosion in a multi-objective solution process can be avoided, the computation required in a redundant robot arm control process can be simplified, and the reaction speed of the redundant robot arm can be improved.

Description

冗余机械臂的控制方法及装置Redundant mechanical arm control method and device 技术领域Technical field
本发明涉及机器人技术领域,尤其涉及一种冗余机械臂的控制方法及装置。The present invention relates to the field of robot technology, and in particular, to a control method and apparatus for a redundant mechanical arm.
背景技术Background technique
近年来,随着人工智能技术和机械控制技术的提高以及机器人所执行任务复杂程度的提高,冗余机械臂(关节数大于6个的机械臂)越来越多的被应用到了机器人中来完成各种较为复杂任务。In recent years, with the improvement of artificial intelligence technology and mechanical control technology and the complexity of tasks performed by robots, redundant mechanical arms (manipulators with more than 6 joints) have been applied to robots to complete. A variety of more complex tasks.
冗余机械臂执行笛卡尔空间特定任务时在关节空间有无穷组解,在具体应用时,一般可以根据关节速度、关节力矩或障碍物距离等优化指标选择一组最优解来确定冗余机械臂的移动过程。在机器人执行较复杂的任务时,通常需要同时对冗余机械臂的多个目标进行优化,目前,常用来进行多目标优化的方法可归类为加权函数法和帕累托前沿法。但是,加权函数法虽然可以利用加权系数对优化目标进行加权从而将多目标优化问题转化为单目标优化问题,但是,此类方法难以判断加权系数变化对优化结果造成的影响,每次只能获得一个解,需不断改变加权系数的值,计算步骤繁琐且存在一定的局限性;帕累托前沿类方法的计算复杂程度随机械臂自由度的增加而急剧增加,当冗余机械臂关节数较多时,该类方法也存在计算复杂、步骤繁琐的问题。When the redundant manipulator performs a Cartesian space-specific task, there is an infinite group solution in the joint space. In the specific application, a set of optimal solutions can be selected according to the optimization index such as joint speed, joint torque or obstacle distance to determine the redundant machine. The movement of the arm. When the robot performs more complex tasks, it is usually necessary to optimize multiple targets of the redundant manipulator at the same time. At present, the methods commonly used for multi-objective optimization can be classified into a weighted function method and a Pareto frontier method. However, although the weighting function method can use the weighting coefficient to weight the optimization target to transform the multi-objective optimization problem into a single-objective optimization problem, it is difficult to judge the influence of the weighting coefficient change on the optimization result. A solution needs to constantly change the value of the weighting coefficient. The calculation steps are cumbersome and have certain limitations. The computational complexity of the Pareto frontier method increases sharply with the increase of the degree of freedom of the robot arm. For a long time, this kind of method also has the problems of complicated calculation and complicated steps.
发明内容Summary of the invention
有鉴于此,本发明实施例提供一种冗余机械臂的控制方法及装置,以解决 现有技术中多目标求解方法计算复杂、步骤繁琐的技术问题。In view of this, an embodiment of the present invention provides a control method and apparatus for a redundant mechanical arm to solve The multi-objective solving method in the prior art calculates a complicated technical problem with complicated steps.
第一方面,本发明实施例提供了一种冗余机械臂的控制方法,包括:In a first aspect, an embodiment of the present invention provides a method for controlling a redundant mechanical arm, including:
获取冗余机械臂所在的当前点信息和目标点信息;Obtaining current point information and target point information where the redundant robot arm is located;
根据所述当前点信息和所述目标点信息确定所述冗余机械臂从当前点移动到目标点的运动轨迹对应的轨迹函数;Determining, according to the current point information and the target point information, a trajectory function corresponding to a motion trajectory of the redundant robot arm moving from a current point to a target point;
以冗余度空间向量为自变量建立所述轨迹函数对应的方程;Establishing an equation corresponding to the trajectory function by using a redundancy space vector as an independent variable;
根据目标接近法求解所述轨迹函数对应的方程,得到所述运动轨迹对应的所述冗余机械臂中各关节的位置和速度。The equation corresponding to the trajectory function is solved according to the target approaching method, and the position and velocity of each joint in the redundant mechanical arm corresponding to the motion trajectory are obtained.
第二方面,本发明实施例还提供了一种冗余机械臂的控制装置,包括:In a second aspect, an embodiment of the present invention further provides a control device for a redundant mechanical arm, including:
信息获取单元,用于获取冗余机械臂所在的当前点信息和目标点信息;An information acquiring unit, configured to acquire current point information and target point information where the redundant robot arm is located;
轨迹单元,与所述信息获取单元相连,用于根据所述当前点信息和所述目标点信息确定所述冗余机械臂从当前点移动到目标点的运动轨迹对应的轨迹函数;a trajectory unit, configured to be connected to the information acquiring unit, configured to determine, according to the current point information and the target point information, a trajectory function corresponding to a motion trajectory of the redundant mechanical arm moving from a current point to a target point;
冗余度处理单元,与所述轨迹单元相连,用于以冗余度空间向量为自变量建立所述轨迹函数对应的方程;根据目标接近法求解所述轨迹函数对应的方程,得到所述运动轨迹对应的所述冗余机械臂中各关节的位置和速度。a redundancy processing unit, connected to the trajectory unit, configured to establish an equation corresponding to the trajectory function by using a redundancy space vector as an independent variable; and solving an equation corresponding to the trajectory function according to a target approach method, to obtain the motion The position and velocity of each joint in the redundant robot arm corresponding to the trajectory.
本发明实施例提供的控制冗余机械臂的技术方案,获取冗余机械臂的当前点信息和目标点信息,根据所获取的当前点信息和目标点信息确定冗余机械臂从当前点移动到目标点的运动轨迹对应的轨迹函数,以冗余度空间向量为自变量建立该轨迹函数对应的方程,通过目标接近法求解该方程以得到运动轨迹对应的冗余机械臂中各关节的位置和速度。本发明实施例通过采用上述技术方案,可以在得到满足设计目标的全局最优解的前提下,减少冗余机械臂多目标优化问题的搜索空间,避免多目标求解过程中维度爆炸问题的发生,简化冗余机械 臂控制过程中所需的计算量,提高冗余机械臂的反应速度。The technical solution for controlling the redundant mechanical arm provided by the embodiment of the present invention obtains current point information and target point information of the redundant mechanical arm, and determines that the redundant mechanical arm moves from the current point to the current point information and the target point information. The trajectory function corresponding to the motion trajectory of the target point is constructed by using the redundancy space vector as the independent variable, and the equation is solved by the target approach method to obtain the position of each joint in the redundant mechanical arm corresponding to the motion trajectory. speed. By adopting the above technical solution, the embodiment of the present invention can reduce the search space of the multi-objective optimization problem of the redundant manipulator under the premise of obtaining the global optimal solution satisfying the design goal, and avoid the occurrence of the dimensional explosion problem in the multi-objective solution process. Simplify redundant machinery The amount of calculation required during the arm control process increases the reaction speed of the redundant manipulator.
附图说明DRAWINGS
通过阅读参照以下附图所作的对非限制性实施例所作的详细描述,本发明的其它特征、目的和优点将会变得更明显:Other features, objects, and advantages of the present invention will become more apparent from the Detailed Description of Description
图1为本发明实施例一提供的一种冗余机械臂的控制方法的流程示意图;1 is a schematic flow chart of a method for controlling a redundant mechanical arm according to Embodiment 1 of the present invention;
图2为本发明实施例二提供的一种冗余机械臂的控制方法的流程示意图;2 is a schematic flow chart of a method for controlling a redundant mechanical arm according to Embodiment 2 of the present invention;
图3为本发明实施例二提供的一种九自由度冗余机械臂的坐标系示意图FIG. 3 is a schematic diagram of a coordinate system of a nine-degree-of-freedom redundant robot arm according to Embodiment 2 of the present invention; FIG.
图4为本发明实施例三提供的一种冗余机械臂的控制装置的结构框图。4 is a structural block diagram of a control device for a redundant mechanical arm according to Embodiment 3 of the present invention.
具体实施方式detailed description
下面结合附图和实施例对本发明作进一步的详细说明。可以理解的是,此处所描述的具体实施例仅仅用于解释本发明,而非对本发明的限定。另外还需要说明的是,为了便于描述,附图中仅示出了与本发明相关的部分而非全部内容。The present invention will be further described in detail below with reference to the accompanying drawings and embodiments. It is understood that the specific embodiments described herein are merely illustrative of the invention and are not intended to limit the invention. It should also be noted that, for ease of description, only some, but not all, of the present invention are shown in the drawings.
实施例一Embodiment 1
本发明实施例一提供一种冗余机械臂的控制方法。该方法可由冗余机械臂的控制装置执行,其中,该装置可通过硬件和/或软件实现,一般可集成在用于控制机器人的控制模块中。图1是本发明实施例一提供的冗余机械臂的控制方法的流程示意图,如图1所示,该方法包括:Embodiment 1 of the present invention provides a control method for a redundant mechanical arm. The method can be performed by a control device of a redundant robot arm, wherein the device can be implemented by hardware and/or software, and can generally be integrated in a control module for controlling the robot. 1 is a schematic flow chart of a method for controlling a redundant mechanical arm according to Embodiment 1 of the present invention. As shown in FIG. 1, the method includes:
S110、获取冗余机械臂所在的当前点信息和目标点信息。S110. Acquire current point information and target point information where the redundant robot arm is located.
本实施例中,冗余机械臂的当前点信息和目标点信息可以是冗余机械臂末 端的当前点信息和目标点信息,也可以是冗余机械臂各关节的当前点信息和目标点信息,此处不作限制。考虑到当前点信息和目标点信息的实用性,可选的,冗余机械臂所在的当前点信息可以包括机械臂末端的当前点信息和冗余机械臂各关节的当前点信息,目标点信息可以包括机械臂末端的目标点信息。其中,当前点信息(目标点信息)可以包括机械臂末端或机械臂各关节在当前(目标)位置的位置坐标、角度坐标、线速度、角速度、线加速度和/或角加速度等信息等。In this embodiment, the current point information and the target point information of the redundant robot arm may be redundant mechanical arm ends. The current point information and the target point information of the terminal may also be the current point information and the target point information of the joints of the redundant robot arm, which are not limited herein. Considering the practicability of the current point information and the target point information, optionally, the current point information of the redundant robot arm may include current point information of the end of the arm and current point information of each joint of the redundant manipulator, target point information Target point information at the end of the robot arm can be included. The current point information (target point information) may include information such as position coordinates, angular coordinates, linear velocity, angular velocity, linear acceleration, and/or angular acceleration of the end of the arm or the joint of the robot arm at the current (target) position.
示例性的,在获取冗余机械臂的当前点信息时,可以直接通过传感器获取冗余机械臂的当前点信息或者通过采用上一时刻的计算结果作为冗余机械臂的当前点信息;在获取冗余机械臂的目标点信息时,可以根据冗余机械臂末端的当前点信息和冗余机械臂末端的运动轨迹信息计算得到冗余机械臂末端需要移动到的目标点信息。以当前点信息为角度坐标、角速度和角加速度为例,冗余机械臂某一关节在当前位置处的角度信息和角速度信息可以通过安装在该关节位置的编码器获得,其角加速度信息可以通过差分运算获得;或者,也可以直接调用上一时刻计算得到的机械臂各关节的目标点信息作为机械臂各关节在当前时刻的当前点信息。Exemplarily, when acquiring the current point information of the redundant manipulator, the current point information of the redundant manipulator can be obtained directly by the sensor or by using the calculation result of the last moment as the current point information of the redundant manipulator; When the target point information of the redundant manipulator is used, the target point information to which the end of the redundant manipulator needs to be moved can be calculated according to the current point information of the end of the redundant manipulator and the motion track information of the end of the redundant manipulator. Taking the current point information as the angular coordinate, angular velocity and angular acceleration as an example, the angular information and angular velocity information of a joint of the redundant mechanical arm at the current position can be obtained by an encoder installed at the joint position, and the angular acceleration information can be differentiated. The operation is obtained; or, the target point information of each joint of the robot arm calculated at the previous time can be directly used as the current point information of each joint of the robot arm at the current time.
S120、根据所述当前点信息和所述目标点信息确定所述冗余机械臂从当前点移动到目标点的运动轨迹对应的轨迹函数。S120. Determine, according to the current point information and the target point information, a trajectory function corresponding to a motion trajectory of the redundant robot arm moving from a current point to a target point.
本实施例中,冗余机械臂从当前点移动到目标点的运动轨迹对应的轨迹函数可以包含在冗余机械臂末端从当前点移动到目标点时冗余机械臂各关节运动轨迹对应的轨迹函数,冗余机械臂各关节的轨迹函数可以是一个或多个,即,冗余机械臂各关节运动轨迹对应的轨迹函数可以以一个总的轨迹函数进行整体描述,也可以以每个关节的运动轨迹对应一个或几个轨迹函数的形式进行描述。 在此,需要指出的是,冗余机械臂各关节的运动轨迹对应的轨迹函数可以包含等式和/或不等式,例如,冗余机械臂某一关节的轨迹函数可以包含该关节的位置坐标、速度和/或加速度与其它变量(如机械臂末端的位置坐标、速度和/或加速度)的对应关系,还可以包含冗余机械臂正常移动时该关节位置坐标、移动速度和/或移动角速度的取值范围等,此处不作限制。In this embodiment, the trajectory function corresponding to the motion trajectory of the redundant mechanical arm moving from the current point to the target point may include a trajectory corresponding to each joint motion trajectory of the redundant mechanical arm when the end of the redundant mechanical arm moves from the current point to the target point. Function, the trajectory function of each joint of the redundant manipulator can be one or more, that is, the trajectory function corresponding to each joint motion trajectory of the redundant manipulator can be described as a whole by a total trajectory function, or can be used for each joint The motion trajectory is described in the form of one or several trajectory functions. Here, it should be pointed out that the trajectory function corresponding to the motion trajectory of each joint of the redundant manipulator may include an equation and/or an inequality. For example, the trajectory function of a joint of the redundant manipulator may include the position coordinates of the joint, The correspondence between speed and/or acceleration and other variables (such as position coordinates, velocity and/or acceleration at the end of the arm) may also include coordinates of the joint position, moving speed, and/or moving angular velocity when the redundant manipulator moves normally. The range of values, etc., is not limited here.
具体的,在确定冗余机械臂从当前点移动到目标点的运动轨迹对应的轨迹函数时,可以首先根据冗余机械臂末端的当前点信息和目标点信息确定冗余机械臂末端的轨迹函数,然后根据冗余机械臂末端与各关节的连接关系确定各关节运动轨迹与冗余机械臂末端运动轨迹之间的关系,进而根据冗余机械臂末端的轨迹函数确定冗余机械臂各关机的轨迹函数。Specifically, when determining a trajectory function corresponding to the motion trajectory of the redundant mechanical arm moving from the current point to the target point, the trajectory function of the end of the redundant mechanical arm may be first determined according to the current point information and the target point information at the end of the redundant mechanical arm. Then, according to the connection relationship between the end of the redundant manipulator and each joint, the relationship between the motion path of each joint and the motion track of the end of the redundant manipulator is determined, and then the redundant robot arm is determined according to the trajectory function of the end of the redundant manipulator. Trajectory function.
S130、以冗余度空间向量为自变量建立所述轨迹函数对应的方程。S130. Establish an equation corresponding to the trajectory function by using a redundancy space vector as an independent variable.
具体的,可以首先确定冗余度向量与各关节位置坐标、速度或加速度等变量的对应关系,然后根据该对应关系采用冗余度空间向量替换各关节对应的轨迹函数中的位置坐标、速度或加速度等自变量,进而建立冗余机械臂各关节以冗余度空间向量为自变量的轨迹函数对应的方程。其中,冗余度空间向量可以是冗余度向量空间的任一向量,冗余度向量空间可以是已有的向量空间,也可以是自定义的空间,此处不作限制。Specifically, the correspondence between the redundancy vector and the coordinates of each joint position, speed, or acceleration may be first determined, and then the position vector, speed, or position in the trajectory function corresponding to each joint is replaced by the redundancy space vector according to the correspondence. Acceleration and other independent variables, and then establish the equation corresponding to the trajectory function of each joint of the redundant manipulator with the redundancy space vector as the independent variable. The redundancy space vector may be any vector of the redundancy vector space, and the redundancy vector space may be an existing vector space or a custom space, which is not limited herein.
S140、根据目标接近法求解所述轨迹函数对应的方程,得到所述运动轨迹对应的所述冗余机械臂中各关节的位置和速度。S140. Solve an equation corresponding to the trajectory function according to a target approaching method, and obtain a position and a velocity of each joint in the redundant mechanical arm corresponding to the motion trajectory.
本实施例中,在得到冗余机械臂各关节的位置和速度后,可以根据需要进一步求得冗余机械臂各关节的加速度。示例性的,当冗余机械臂某一关节的轨迹函数对应的方程存在有限组解时,可以求得方程的所有解,然后根据实际需要或根据设定的选取规则选取方程的一组解得到该关节在下一时刻的目标位置 并同时确定该关节从当前位置移动到目标位置时的移动速度和加速度。当冗余机械臂某一关节的轨迹函数对应的方程存在无穷组解时,可以任意求得该方程的一组解并通过该组解得到该关节在下一时刻的目标位置以及从当前位置移动到目标位置的移动速度和加速度;也可以首先设定方程解的约束条件,然后求符合该约束条件的方程的解,并以该组解指导该关节从当前时刻到下一时刻的移动过程,此处不作限制。为保证冗余机械臂各关节都可以高效地移动,减少各关节的输出力矩和运动、规避奇异与避免关节运动的极限,优选的,可以根据符合约束条件的方程的解得到运动轨迹对应的冗余机械臂中各关节的位置、速度和加速度。In this embodiment, after obtaining the position and speed of each joint of the redundant mechanical arm, the acceleration of each joint of the redundant mechanical arm can be further determined as needed. Illustratively, when there is a finite set of solutions for the equation corresponding to the trajectory function of a joint of a redundant manipulator, all solutions of the equation can be obtained, and then a set of solutions of the equations can be obtained according to actual needs or according to the set selection rules. The target position of the joint at the next moment At the same time, the moving speed and acceleration of the joint when moving from the current position to the target position are determined. When there is an infinite group solution for the equation corresponding to the trajectory function of a joint of a redundant manipulator, a set of solutions of the equation can be arbitrarily obtained and the target position of the joint at the next moment is obtained by the set of solutions and moved from the current position to The moving speed and acceleration of the target position; first, the constraint condition of the equation solution can be set first, then the solution of the equation that satisfies the constraint condition is obtained, and the movement process of the joint from the current time to the next moment is guided by the set of solutions, There are no restrictions. In order to ensure that the joints of the redundant manipulator can move efficiently, reducing the output torque and motion of each joint, avoiding the singularity and avoiding the limit of joint motion, it is preferable to obtain the redundancy corresponding to the motion trajectory according to the solution of the equation satisfying the constraint condition. The position, velocity and acceleration of the joints in the remaining manipulators.
本发明实施例一提供的冗余机械臂的控制方法,获取冗余机械臂的当前点信息和目标点信息,根据所获取的当前点信息和目标点信息确定冗余机械臂从当前点移动到目标点的运动轨迹对应的轨迹函数,以冗余度空间向量为自变量建立该轨迹函数对应的方程,通过目标接近法求解该方程以得到运动轨迹对应的冗余机械臂中各关节的位置和速度。本实施例通过采用上述技术方案,可以在得到满足设计目标的全局最优解的前提下,减少冗余机械臂多目标优化问题的搜索空间,避免多目标求解过程中维度爆炸问题的发生,简化冗余机械臂控制过程中所需的计算量,提高冗余机械臂的反应速度。The control method of the redundant mechanical arm provided by the first embodiment of the present invention acquires current point information and target point information of the redundant mechanical arm, and determines that the redundant mechanical arm moves from the current point to the current point information and the target point information. The trajectory function corresponding to the motion trajectory of the target point is constructed by using the redundancy space vector as the independent variable, and the equation is solved by the target approach method to obtain the position of each joint in the redundant mechanical arm corresponding to the motion trajectory. speed. By adopting the above technical solution, the present embodiment can reduce the search space of the multi-objective optimization problem of the redundant manipulator under the premise of satisfying the global optimal solution of the design target, and avoid the occurrence of the dimensional explosion problem in the multi-objective solution process, simplifying The amount of calculation required during the redundant manipulator control process increases the reaction speed of the redundant manipulator.
实施例二Embodiment 2
图2为本发明实施例二提供的一种冗余机械臂的控制方法的流程示意图。本实施例在上述实施例的基础上进行优化,进一步地,所述获取冗余机械臂所在的当前点信息和目标点信息包括:获取所述冗余机械臂末端在所述目标点的速度;根据所述冗余机械臂末端在目标当前点的速度将所述冗余机械臂的关节 角、关节角速度和关节角加速度映射到冗余度空间。2 is a schematic flow chart of a method for controlling a redundant mechanical arm according to Embodiment 2 of the present invention. The embodiment is optimized on the basis of the foregoing embodiment. Further, the acquiring the current point information and the target point information where the redundant mechanical arm is located includes: acquiring the speed of the end of the redundant mechanical arm at the target point; According to the speed of the end of the redundant manipulator at the target current point, the joint of the redundant manipulator The angular, joint angular velocity and joint angular acceleration are mapped to the redundancy space.
进一步地,所述根据所述当前点信息和所述目标点信息确定所述冗余机械臂从当前点移动到目标点的运动轨迹对应的轨迹函数包括:根据逆运动学方程,确定所述冗余机械臂从当前点移动到目标点的运动轨迹对应的优化目标函数和约束条件。Further, the determining, according to the current point information and the target point information, a trajectory function corresponding to a motion trajectory of the redundant robot arm moving from a current point to a target point comprises: determining the redundancy according to an inverse kinematics equation The optimization objective function and constraints corresponding to the motion trajectory of the remaining manipulator moving from the current point to the target point.
进一步地,所述根据目标接近法求解所述轨迹函数对应的方程,得到所述运动轨迹对应的所述冗余机械臂中各关节的位置和速度包括:在冗余度空间中,根据优化目标函数、与优化目标对应的最小目标值、对应的约束条件和预设的加权系数,通过辅助向量和单目标优化算法求解所述轨迹函数对应的方程,得到所述运动轨迹对应的所述冗余机械臂中各关节的位置和速度。Further, the solving the equation corresponding to the trajectory function according to the target approaching method, and obtaining the position and velocity of each joint in the redundant mechanical arm corresponding to the motion trajectory includes: in the redundancy space, according to the optimization target a function, a minimum target value corresponding to the optimization target, a corresponding constraint condition, and a preset weighting coefficient, and an equation corresponding to the trajectory function is solved by an auxiliary vector and a single target optimization algorithm to obtain the redundancy corresponding to the motion trajectory The position and speed of the joints in the robot arm.
相应的,如图2所示,本实施例提供的冗余机械臂的控制方法包括:Correspondingly, as shown in FIG. 2, the control method of the redundant mechanical arm provided by this embodiment includes:
S210、获取所述冗余机械臂末端在所述目标点的速度。S210. Acquire a speed of the end of the redundant mechanical arm at the target point.
优选的,所述获取所述冗余机械臂末端在所述当前点的速度包括:根据所述冗余机械臂末端在所述当前点的位姿和所述冗余机械臂末端在所述运动轨迹中的下一采样点的位姿确定所述冗余机械臂末端在所述下一采样点的速度。其中,冗余机械臂末端的位姿指的是冗余机械臂的位置和姿态,位置可以用位置坐标来表示,姿态可以用角度来表示。示例性的,可以根据安装在冗余机械臂末端的编码器获取冗余机械臂在当前点的角速度和/或角加速度等信息,或者直接根据冗余机械臂在上一时刻的计算结果得到冗余机械臂末端在当前点的角速度和/或角加速度等信息,然后根据冗余机械臂在当前点的角度和下一采样的角度得到冗余机械臂在下一采样点的角速度,进而得到冗余机械臂在下一采样点的速度。Preferably, the obtaining the speed of the end of the redundant mechanical arm at the current point comprises: according to the posture of the redundant mechanical arm end at the current point and the movement of the redundant mechanical arm end The pose of the next sample point in the trajectory determines the velocity of the end of the redundant robot arm at the next sample point. The position of the end of the redundant manipulator refers to the position and posture of the redundant manipulator. The position can be represented by position coordinates, and the posture can be expressed by angle. Illustratively, information such as the angular velocity and/or angular acceleration of the redundant manipulator at the current point can be obtained from an encoder installed at the end of the redundant manipulator, or directly based on the calculation result of the redundant manipulator at the previous moment. The angular velocity and/or angular acceleration of the end of the mechanical arm at the current point, and then the angular velocity of the redundant manipulator at the next sampling point according to the angle of the redundant manipulator at the current point and the angle of the next sampling, thereby obtaining redundancy The speed of the arm at the next sampling point.
S220、根据所述冗余机械臂末端在所述目标点的速度将所述冗余机械臂的 关节角、关节角速度和关节角加速度映射到冗余度空间。S220, according to the speed of the redundant mechanical arm end at the target point, the redundant mechanical arm Joint angle, joint angular velocity, and joint angular acceleration are mapped to the redundancy space.
示例性的,冗余机械臂末端速度与冗余机械臂的关节角速度的关系(即冗余机械臂各关节以角速度为自变量的轨迹函数)可以表示为:Exemplarily, the relationship between the end speed of the redundant manipulator and the joint angular velocity of the redundant manipulator (ie, the trajectory function of each joint of the redundant manipulator with angular velocity as an independent variable) can be expressed as:
Figure PCTCN2017103286-appb-000001
Figure PCTCN2017103286-appb-000001
其中,
Figure PCTCN2017103286-appb-000002
为机械臂末端的速度,
Figure PCTCN2017103286-appb-000003
vi(i=x,y,z)表示线速度,ωi(i=θ,Ψ,Φ)表示角速度;
Figure PCTCN2017103286-appb-000004
为机械臂的关节角速度,
Figure PCTCN2017103286-appb-000005
N为冗余机械臂的自由度,N=6+r,r>1;J∈R6×N为雅克比矩阵。
among them,
Figure PCTCN2017103286-appb-000002
For the speed of the end of the arm,
Figure PCTCN2017103286-appb-000003
v i (i=x, y, z) represents a linear velocity, and ω i (i = θ, Ψ, Φ) represents an angular velocity;
Figure PCTCN2017103286-appb-000004
For the joint angular velocity of the arm,
Figure PCTCN2017103286-appb-000005
N is the degree of freedom of the redundant manipulator, N=6+r, r>1; J∈R 6×N is the Jacobian matrix.
假设α是雅可比矩阵J中的任意r(r=n-6>1)列,使得剩下的6列构成非奇异矩阵J*,则式(1)可以改写为:Assuming that α is an arbitrary r (r=n-6>1) column in the Jacobian matrix J, so that the remaining 6 columns constitute a non-singular matrix J * , then equation (1) can be rewritten as:
Figure PCTCN2017103286-appb-000006
Figure PCTCN2017103286-appb-000006
从而,任何满足式(2)的关节角速度
Figure PCTCN2017103286-appb-000007
可以进一步改写为:
Thus, any joint angular velocity that satisfies equation (2)
Figure PCTCN2017103286-appb-000007
Can be further rewritten as:
Figure PCTCN2017103286-appb-000008
Figure PCTCN2017103286-appb-000008
其中,
Figure PCTCN2017103286-appb-000009
是满足式(3)的一个特解,
Figure PCTCN2017103286-appb-000010
是式(3)的齐次解,且
among them,
Figure PCTCN2017103286-appb-000009
Is a special solution that satisfies equation (3).
Figure PCTCN2017103286-appb-000010
Is a homogeneous solution of equation (3), and
Figure PCTCN2017103286-appb-000011
Figure PCTCN2017103286-appb-000011
Figure PCTCN2017103286-appb-000012
其中,i=1,2,...,r,αj是矩阵α的第j列;
Figure PCTCN2017103286-appb-000013
为列矩阵,且其前r(r=N-6)行为0;
Figure PCTCN2017103286-appb-000014
为列矩阵,且其第i行为1,前r行中除第i行外的其他行为0。定义冗余度空间向量k=(k1,…,kr)T为冗余度向量空间的任一向量,则有:
make
Figure PCTCN2017103286-appb-000012
Where i=1, 2, . . . , r, α j is the jth column of the matrix α;
Figure PCTCN2017103286-appb-000013
Is a column matrix, and its pre-r(r=N-6) behaves as 0;
Figure PCTCN2017103286-appb-000014
It is a column matrix, and its ith behavior is 1, and the other behaviors of the first r rows except the ith row are 0. Defining the redundancy space vector k=(k 1 ,...,k r ) T is any vector of the redundancy vector space, then:
Figure PCTCN2017103286-appb-000015
Figure PCTCN2017103286-appb-000015
假设角速度
Figure PCTCN2017103286-appb-000016
以S(k)进行表示,角度θ以P(k)进行表示,角加速度
Figure PCTCN2017103286-appb-000017
以A(k)进行表示,则冗余机械臂关节角速度
Figure PCTCN2017103286-appb-000018
可表示为冗余度空间向量k的函数:
Assumed angular velocity
Figure PCTCN2017103286-appb-000016
Expressed as S(k), the angle θ is represented by P(k), angular acceleration
Figure PCTCN2017103286-appb-000017
Represented by A(k), the angular velocity of the redundant manipulator joint
Figure PCTCN2017103286-appb-000018
Can be expressed as a function of the redundancy space vector k:
Figure PCTCN2017103286-appb-000019
Figure PCTCN2017103286-appb-000019
冗余机械臂关节角度θ关于冗余度空间向量k的表达式为:The expression of the redundant manipulator joint angle θ with respect to the redundancy space vector k is:
Figure PCTCN2017103286-appb-000020
Figure PCTCN2017103286-appb-000020
冗余机械臂关节角加速度关于冗余度空间向量k的表达式为:Redundant manipulator joint angular acceleration The expression for the redundancy space vector k is:
Figure PCTCN2017103286-appb-000022
Figure PCTCN2017103286-appb-000022
其中,θ0为上一时刻的关节角度,
Figure PCTCN2017103286-appb-000023
为上一时刻的关节角速度,Δt为采样步长。
Where θ 0 is the joint angle of the previous moment,
Figure PCTCN2017103286-appb-000023
For the joint angular velocity at the previous moment, Δt is the sampling step size.
S230、根据逆运动学方程,确定所述冗余机械臂从当前点移动到目标点的运动轨迹对应的优化目标函数和约束条件。S230. Determine, according to an inverse kinematics equation, an optimization objective function and a constraint condition corresponding to a motion trajectory of the redundant mechanical arm moving from a current point to a target point.
本实施例中,冗余机械臂从当前点移动到目标点的运动轨迹对应的优化目标函数和约束条件可以根据需要灵活进行设定,例如,可以对冗余机械臂移动的范围进行限制以躲避障碍物,或者,对冗余机械臂的弯曲角度进行限制以避免各关节达到其运动极限等。考虑到各约束条件的实用性,可选的,所述优化目标函数包括所述冗余机械臂的关节运动幅度最小或者所述冗余机械臂受对应的齿轮间隙影响最小;所述约束条件包括:将所述冗余机械臂的关节速度约束在预设的速度范围内或者将所述冗余机械臂的关节角度约束在预设的角度范围内。示例性的,冗余机械臂关节运动幅度最小的优化目标函数可以为:
Figure PCTCN2017103286-appb-000024
冗余机械臂受对应的齿轮间隙影响最小的优化目标函数可以为:
Figure PCTCN2017103286-appb-000025
将冗余机械臂的关节速度约束在预设的速度范围内的约束条件可以为:
Figure PCTCN2017103286-appb-000026
将冗余机械臂的关节角度约束在预设的角度范围内 的约束条件可以为:θmin≤θ≤θmax
In this embodiment, the optimization objective function and the constraint condition corresponding to the motion trajectory of the redundant mechanical arm moving from the current point to the target point can be flexibly set according to requirements, for example, the range of the redundant mechanical arm movement can be restricted to avoid Obstacle, or limit the bending angle of the redundant manipulator to avoid the joints reaching their motion limits. In view of the practicability of each constraint, optionally, the optimization objective function includes a minimum joint motion amplitude of the redundant mechanical arm or the redundant mechanical arm is minimally affected by a corresponding gear gap; the constraint condition includes : constraining the joint speed of the redundant mechanical arm within a preset speed range or constraining the joint angle of the redundant mechanical arm within a preset angular range. Illustratively, the optimal objective function for the minimum joint motion of the redundant manipulator can be:
Figure PCTCN2017103286-appb-000024
The optimal objective function for the redundant manipulator to be minimally affected by the corresponding gear clearance can be:
Figure PCTCN2017103286-appb-000025
The constraint that constrains the joint speed of the redundant manipulator to a preset speed range may be:
Figure PCTCN2017103286-appb-000026
The constraint that constrains the joint angle of the redundant manipulator within a predetermined range of angles may be: θ min ≤ θ ≤ θ max .
S240、以冗余度空间向量为自变量建立所述轨迹函数对应的方程。S240. The equation corresponding to the trajectory function is established by using a redundancy space vector as an independent variable.
本实施例中,可以基于冗余度空间建立冗余机械臂轨迹的优化目标函数和约束方程,由于冗余机械臂的优化目标函数和约束方程都可以表示为以关节角度、关节角速度和/或关节角加速度为自变量的表达式,因此,将上式(7)、(8)、(9)代入表达式中即可得到以冗余度空间向量k为自变量的优化目标函数和约束方程。例如,关于关节运动幅度最小的优化目标函数可以进一步表示为冗余度空间向量k的方程:In this embodiment, the optimization objective function and the constraint equation of the redundant robot arm trajectory can be established based on the redundancy space, and the optimization objective function and the constraint equation of the redundant mechanical arm can be expressed as joint angle, joint angular velocity and/or The joint angular acceleration is an expression of the independent variable. Therefore, by substituting the above equations (7), (8), and (9) into the expression, the optimal objective function and the constraint equation with the redundant space vector k as the independent variable can be obtained. . For example, an optimization objective function with minimal joint motion amplitude can be further expressed as an equation for the redundancy space vector k:
Figure PCTCN2017103286-appb-000027
Figure PCTCN2017103286-appb-000027
关于机械臂受对应的齿轮间隙影响最小的优化目标函数可以进一步表示为冗余度空间向量k的方程:The optimization objective function with respect to the minimum influence of the mechanical arm on the corresponding gear gap can be further expressed as the equation of the redundancy space vector k:
Figure PCTCN2017103286-appb-000028
Figure PCTCN2017103286-appb-000028
关于关节角速度限制的约束条件可以进一步表示为以冗余度空间向量k为自变量的方程:The constraint on the joint angular velocity limit can be further expressed as an equation with the redundant space vector k as an independent variable:
Figure PCTCN2017103286-appb-000029
Figure PCTCN2017103286-appb-000029
关于关节角度限制的约束条件可以进一步表示为以冗余度空间向量k为自变量的方程:The constraint on the joint angle limit can be further expressed as an equation with the redundancy space vector k as an independent variable:
Figure PCTCN2017103286-appb-000030
Figure PCTCN2017103286-appb-000030
S250、在冗余度空间中,根据优化目标函数、与优化目标对应的最小目标值、对应的约束条件和预设的加权系数,通过辅助向量和单目标优化算法求解所述轨迹函数对应的方程,得到所述运动轨迹对应的所述冗余机械臂中各关节 的位置和速度。S250. In the redundancy space, according to the optimization objective function, the minimum target value corresponding to the optimization target, the corresponding constraint condition, and the preset weighting coefficient, the equation corresponding to the trajectory function is solved by the auxiliary vector and the single target optimization algorithm. Obtaining each joint in the redundant mechanical arm corresponding to the motion trajectory The location and speed.
示例性的,求解轨迹函数对应的方程的单目标优化算法可以根据需要灵活确定,此处不作限制。优选的,所述单目标优化算法包括:牛顿-欧拉算法、Nelder-Mead单纯形算法、内点算法、基因算法和/或Pattern Search算法。Exemplarily, the single-objective optimization algorithm for solving the equation corresponding to the trajectory function can be flexibly determined according to needs, and is not limited herein. Preferably, the single target optimization algorithm comprises: a Newton-Eulerian algorithm, a Nelder-Mead simplex algorithm, an interior point algorithm, a genetic algorithm and/or a Pattern Search algorithm.
帕累托前沿是多目标优化问题的解集,其定义如下:在多目标可行域Ω上,存在一点x*∈Ω,若不存在点x∈Ω,使得f(x)≤f(x*)成立,且至少存在一点x′∈Ω使严格不等式f(x′)>f(x*)成立,则x*为一个帕累托最优解。以双目标优化问题为例,其轨迹函数对应的方程如下:The Pareto frontier is a solution set for multi-objective optimization problems, which is defined as follows: On the multi-objective feasible domain Ω, there is a point x * ∈ Ω, if there is no point x ∈ Ω, so that f(x) ≤ f(x * ) is established, and there is at least a point x' ∈ Ω such that the strict inequality f(x') > f(x * ) holds, then x * is a Pareto optimal solution. Taking the double-objective optimization problem as an example, the equation corresponding to the trajectory function is as follows:
Figure PCTCN2017103286-appb-000031
Figure PCTCN2017103286-appb-000031
上式中,F(x)={f1(x) f2(x)}为优化目标函数,
Figure PCTCN2017103286-appb-000032
为与优化目标函数对应的最小目标值,A为矩阵约束条件,b为向量约束条件,ω为加权系数向量,γ为辅助向量。其中,加权系数向量ω可由用户或开发商根据需要进行设置,辅助向量γ用于将原来的多目标优化为题转换为单目标优化问题。
In the above formula, F(x)={f 1 (x) f 2 (x)} is the optimization objective function,
Figure PCTCN2017103286-appb-000032
For the minimum target value corresponding to the optimization objective function, A is a matrix constraint condition, b is a vector constraint condition, ω is a weighting coefficient vector, and γ is an auxiliary vector. The weighting coefficient vector ω can be set by the user or the developer as needed, and the auxiliary vector γ is used to convert the original multi-objective optimization into a single-objective optimization problem.
本实施例中,可试用内点法求解式(12)得到方程的最优解,进而通过式(7)、式(8)和式(9)得到冗余机械臂各关节与该最优解对应的角度、角速度和角角速度,从而得到冗余机械臂各关节在下一时刻的位置和速度。In this embodiment, the inner point method can be used to solve the equation (12) to obtain the optimal solution of the equation, and then the joints of the redundant mechanical arm and the optimal solution are obtained by the equations (7), (8) and (9). Corresponding angles, angular velocities, and angular angular velocities, resulting in the position and velocity of the joints of the redundant manipulator at the next moment.
以九自由度的冗余机械臂为例(假设该九自由度冗余机械臂的D-H坐标如图3所示,该九自由度冗余机械臂连杆D-H参数表如表1所示),求解冗余机械臂各关节位置和速度的过程可以为:Taking a redundant manipulator with nine degrees of freedom as an example (assuming that the DH coordinates of the nine-degree-of-freedom redundant manipulator are as shown in Fig. 3, the parameter table of the nine-degree-of-freedom redundant manipulator link DH is shown in Table 1) The process of solving the joint position and speed of the redundant manipulator can be:
表1Table 1
θi θ i di(mm)d i (mm) αi α i ai(mm)a i (mm) 范围range
00 d1 d 1 90°90° 00 1550≤d1≤82761550 ≤ d 1 ≤ 8276
θ2 θ 2 00 00 a2=2200a 2 =2200 -90≤θ2≤90-90≤θ 2 ≤90
θ3 θ 3 00 -90°-90° 00 -90≤θ3≤90-90≤θ 3 ≤90
θ4 θ 4 d4=2800d 4 = 2800 90°90° α4=-375α 4 =-375 -180≤θ4≤180-180≤θ 4 ≤180
θ5 θ 5 00 -90°-90° α5=-a4=375α 5 =-a 4 =375 -180≤θ5≤180-180≤θ 5 ≤180
θ6 θ 6 d4=2800d 4 = 2800 90°90° 00 -180≤θ6≤180-180≤θ 6 ≤180
θ7 θ 7 00 -90°-90° 00 -90≤θ7≤90-90≤θ 7 ≤90
θ8 θ 8 d8=1650d 8 =1650 90°90° 00 -90≤θ8≤90-90≤θ 8 ≤90
θ9 θ 9 00 00 ah=350a h =350 -90≤θ9≤90-90≤θ 9 ≤90
首先,把冗余机械臂末端的速度转换为以连杆3的坐标系为参考:First, the speed of the end of the redundant manipulator is converted to the reference of the coordinate system of the link 3:
Figure PCTCN2017103286-appb-000033
Figure PCTCN2017103286-appb-000033
Figure PCTCN2017103286-appb-000034
Figure PCTCN2017103286-appb-000034
然后,定义腕部坐标系。考虑到计算的简便性,可以根据冗余机械臂的结构特点定义腕部坐标系,其中,腕部坐标系的具体坐标与方向可以根据需要灵活设置,此处不作限定。示例性的,可以定义腕部坐标系为与冗余机械臂末端坐标系平行、原点与连杆8坐标系原点重合的坐标系,此时,当腕部和冗余机械臂末端的速度以连杆3坐标系为参考时,二者的角速度和线速度存在如下关系:Then, define the wrist coordinate system. Considering the simplicity of the calculation, the wrist coordinate system can be defined according to the structural characteristics of the redundant mechanical arm. The specific coordinates and direction of the wrist coordinate system can be flexibly set according to requirements, which is not limited herein. Exemplarily, the wrist coordinate system can be defined as a coordinate system parallel to the coordinate end system of the redundant manipulator and the origin coincides with the origin of the coordinate system of the connecting rod 8. At this time, when the speed of the end of the wrist and the redundant manipulator is connected When the rod 3 coordinate system is a reference, the angular velocity and the linear velocity of the two have the following relationship:
3ωw3ωn    (15) 3 ω w = 3 ω n (15)
3vw3vh-3Rw(wωh×wPh)=3vh-3ωh×3Ph3vh-3ωh×a h 3xh    (16) 3 v w = 3 v h - 3 R w ( w ω h × w P h )= 3 v h - 3 ω h × 3 P h = 3 v h - 3 ω h × a h 3 x h (16)
其中,3ωw是以连杆3坐标系为参考的腕部坐标系中的角速度;3ωh是以连杆3坐标系为参考的冗余机械臂末端坐标系中的角速度;3vw是以连杆3坐标系为参考的腕部坐标系中的线速度;3vh是以连杆3坐标系为参考的冗余机械臂末端坐标系中的线速度;3xh是以连杆3坐标系为参考的冗余机械臂末端坐标系中Z方向的单位矢量的值: Where 3 ω w is the angular velocity in the wrist coordinate system with reference to the coordinate system of the connecting rod 3; 3 ω h is the angular velocity in the coordinate system of the redundant manipulator referenced by the coordinate system of the connecting rod 3; 3 v w The linear velocity in the wrist coordinate system with reference to the coordinate system of the connecting rod 3; 3 v h is the linear velocity in the coordinate system of the redundant manipulator referenced by the coordinate system of the connecting rod 3; 3 x h is the connection The rod 3 coordinate system is the value of the unit vector of the Z direction in the redundant robot arm end coordinate system of the reference:
Figure PCTCN2017103286-appb-000035
Figure PCTCN2017103286-appb-000035
上式中,ci表示cos(θi),si表示sin(θi),cij表示cos(θij),sij表示sin(θij)。In the above formula, c i represents cos(θ i ), s i represents sin(θ i ), c ij represents cos(θ ij ), and s ij represents sin(θ ij ).
定义3Jw为联系关节角速度
Figure PCTCN2017103286-appb-000036
与腕部速度
Figure PCTCN2017103286-appb-000037
的雅克比矩阵,则有:
Definition 3 J w is the joint angular velocity
Figure PCTCN2017103286-appb-000036
Speed with wrist
Figure PCTCN2017103286-appb-000037
The Jacques matrix has:
Figure PCTCN2017103286-appb-000038
Figure PCTCN2017103286-appb-000038
用矢量积法可求出该雅克比矩阵如下:The Jacobian matrix can be found by the vector product method as follows:
3Jw=[J1 J2 J3 J4 J5 J6 J7 J8 J9]    (19) 3 J w =[J 1 J 2 J 3 J 4 J 5 J 6 J 7 J 8 J 9 ] (19)
其中,among them,
Figure PCTCN2017103286-appb-000039
Figure PCTCN2017103286-appb-000039
Figure PCTCN2017103286-appb-000040
Figure PCTCN2017103286-appb-000040
Figure PCTCN2017103286-appb-000041
Figure PCTCN2017103286-appb-000041
Figure PCTCN2017103286-appb-000042
Figure PCTCN2017103286-appb-000042
Figure PCTCN2017103286-appb-000043
Figure PCTCN2017103286-appb-000043
Figure PCTCN2017103286-appb-000044
Figure PCTCN2017103286-appb-000044
本实施例中,可在式(19)中取任意三列作为α,J*为除组成α的三列之外的其余列。示例性的,取J5、J6和J7作为α,其余列(J1、J2、J3、J4、J8和J9)组成J*,通过式(5)可以得到:In the present embodiment, any three columns in the equation (19) may be taken as α, and J * is the remaining columns other than the three columns constituting α. Illustratively, J 5 , J 6 and J 7 are taken as α, and the remaining columns (J 1 , J 2 , J 3 , J 4 , J 8 and J 9 ) constitute J * , and can be obtained by the formula (5):
Figure PCTCN2017103286-appb-000045
Figure PCTCN2017103286-appb-000045
其中,among them,
J21=a2s3-d6c5-d4-a5s5-d7c5c7+d7s5c6s7J 21 = a 2 s 3 - d 6 c 5 - d 4 - a 5 s 5 - d 7 c 5 c 7 + d 7 s 5 c 6 s 7 ,
J23=a2c3+a4c4+a5c4c5-d6c4s5-d7c4s5c7+d7s4s6s7-d7c4c5c6s7J 23 = a 2 c 3 + a 4 c 4 + a 5 c 4 c 5 - d 6 c 4 s 5 - d 7 c 4 s 5 c 7 + d 7 s 4 s 6 s 7 -d 7 c 4 c 5 c 6 s 7 ,
J31=-d4-d7(c5c7-s5c6s7)-d6c5-a5s5J 31 =-d 4 -d 7 (c 5 c 7 -s 5 c 6 s 7 )-d 6 c 5 -a 5 s 5 ,
J33=a4c4+d7(s7(s4s6-c4c5c6)-c4s5c7)+a5c4c5-d6c4s3J 33 = a 4 c 4 + d 7 (s 7 (s 4 s 6 -c 4 c 5 c 6 )-c 4 s 5 c 7 )+a 5 c 4 c 5 -d 6 c 4 s 3 ,
J41=d7(s7(c4s6+s4c5c6)+s4s5c7)-a4s4-a5s4c5+d6s4s5J 41 =d 7 (s 7 (c 4 s 6 +s 4 c 5 c 6 )+s 4 s 5 c 7 )-a 4 s 4 -a 5 s 4 c 5 +d 6 s 4 s 5 ,
J42=d7(s7(s4s6-c4c5c6)-c4s5c7)+a4c4+a5c4c5-d6c4s5J 42 =d 7 (s 7 (s 4 s 6 -c 4 c 5 c 6 )-c 4 s 5 c 7 )+a 4 c 4 +a 5 c 4 c 5 -d 6 c 4 s 5 ,
J84=s7(s4s6-c4c5c6)-c4s5c7 J 84 = s 7 (s 4 s 6 -c 4 c 5 c 6) -c 4 s 5 c 7,
J85=-s7(c4s6+s4c5c6)-s4s5c7 J 85 = -s 7 (c 4 s 6 + s 4 c 5 c 6) -s 4 s 5 c 7,
J86=c5c7-s5c6s7J 86 =c 5 c 7 -s 5 c 6 s 7 ,
J94=c8(s4c6+c4c5s6)-s8(c7(s4s6-c4c5c6))+c4s5s7J 94 =c 8 (s 4 c 6 +c 4 c 5 s 6 )-s 8 (c 7 (s 4 s 6 -c 4 c 5 c 6 ))+c 4 s 5 s 7 ,
J95=s8(c7(s4s6+s4c5c6))-s4s5s7-c8(c4c6-s4c5s6),J 95 =s 8 (c 7 (s 4 s 6 +s 4 c 5 c 6 ))-s 4 s 5 s 7 -c 8 (c 4 c 6 -s 4 c 5 s 6 ),
J96=s8(c5s7-s5c6c7)+c8s5s6J 96 = s 8 (c 5 s 7 - s 5 c 6 c 7 ) + c 8 s 5 s 6 .
Figure PCTCN2017103286-appb-000046
Figure PCTCN2017103286-appb-000047
由式(20)可得:
make
Figure PCTCN2017103286-appb-000046
with
Figure PCTCN2017103286-appb-000047
From equation (20):
Figure PCTCN2017103286-appb-000048
Figure PCTCN2017103286-appb-000048
Figure PCTCN2017103286-appb-000049
Figure PCTCN2017103286-appb-000049
Figure PCTCN2017103286-appb-000050
Figure PCTCN2017103286-appb-000050
Figure PCTCN2017103286-appb-000051
Figure PCTCN2017103286-appb-000051
Figure PCTCN2017103286-appb-000052
Figure PCTCN2017103286-appb-000052
Figure PCTCN2017103286-appb-000053
Figure PCTCN2017103286-appb-000053
其中,
Figure PCTCN2017103286-appb-000054
为矩阵
Figure PCTCN2017103286-appb-000055
的第i个元素。
among them,
Figure PCTCN2017103286-appb-000054
Matrix
Figure PCTCN2017103286-appb-000055
The ith element.
由式(6)可得:Available from formula (6):
Figure PCTCN2017103286-appb-000056
Figure PCTCN2017103286-appb-000056
其中,α=(J5 J6 J7)。Where α = (J 5 J 6 J 7 ).
因此,将式(21)至(26)中的
Figure PCTCN2017103286-appb-000057
(j=1,2,3,4,8和9)用
Figure PCTCN2017103286-appb-000058
(i=5,6和7,j=1,2,3,4,8和9)代替,将式(21)至(26)中的
Figure PCTCN2017103286-appb-000059
(k=1,2,3,4,5和6)用-Jik(i=5,6和7,k=1,2,3,4,5和6)代替即可得到
Figure PCTCN2017103286-appb-000060
Figure PCTCN2017103286-appb-000061
的值。以
Figure PCTCN2017103286-appb-000062
为例,将式(21)至(26)中的
Figure PCTCN2017103286-appb-000063
(j=1,2,3,4,8和9)用
Figure PCTCN2017103286-appb-000064
(j=1,2,3,4,8和9)代替,可以得到:
Therefore, in the equations (21) to (26)
Figure PCTCN2017103286-appb-000057
(j=1, 2, 3, 4, 8 and 9)
Figure PCTCN2017103286-appb-000058
(i=5, 6 and 7, j=1, 2, 3, 4, 8 and 9) instead, in the equations (21) to (26)
Figure PCTCN2017103286-appb-000059
(k=1, 2, 3, 4, 5, and 6) can be obtained by replacing -J ik (i=5,6 and 7, k=1, 2, 3, 4, 5, and 6)
Figure PCTCN2017103286-appb-000060
with
Figure PCTCN2017103286-appb-000061
Value. Take
Figure PCTCN2017103286-appb-000062
As an example, the equations (21) to (26)
Figure PCTCN2017103286-appb-000063
(j=1, 2, 3, 4, 8 and 9)
Figure PCTCN2017103286-appb-000064
(j=1, 2, 3, 4, 8 and 9) instead, you can get:
Figure PCTCN2017103286-appb-000065
Figure PCTCN2017103286-appb-000065
将式(21)至(26)中的
Figure PCTCN2017103286-appb-000066
(k=1,2,3,4,5和6)用-J5k(k=1,2,3,4,5和6)代替,可得到:
In the formulas (21) to (26)
Figure PCTCN2017103286-appb-000066
(k=1, 2, 3, 4, 5, and 6) Substituting -J 5k (k = 1, 2, 3, 4, 5, and 6) gives:
Figure PCTCN2017103286-appb-000067
Figure PCTCN2017103286-appb-000067
由式(19)可知
Figure PCTCN2017103286-appb-000068
因此可得:
Known by equation (19)
Figure PCTCN2017103286-appb-000068
So available:
Figure PCTCN2017103286-appb-000069
Figure PCTCN2017103286-appb-000069
Figure PCTCN2017103286-appb-000070
Figure PCTCN2017103286-appb-000070
Figure PCTCN2017103286-appb-000071
Figure PCTCN2017103286-appb-000071
Figure PCTCN2017103286-appb-000072
Figure PCTCN2017103286-appb-000072
Figure PCTCN2017103286-appb-000073
Figure PCTCN2017103286-appb-000073
Figure PCTCN2017103286-appb-000074
Figure PCTCN2017103286-appb-000074
其中,
Figure PCTCN2017103286-appb-000075
为矩阵
Figure PCTCN2017103286-appb-000076
的第i个元素。同理可以得出:
among them,
Figure PCTCN2017103286-appb-000075
Matrix
Figure PCTCN2017103286-appb-000076
The ith element. The same can be concluded:
Figure PCTCN2017103286-appb-000077
Figure PCTCN2017103286-appb-000077
Figure PCTCN2017103286-appb-000078
Figure PCTCN2017103286-appb-000078
Figure PCTCN2017103286-appb-000079
Figure PCTCN2017103286-appb-000079
Figure PCTCN2017103286-appb-000080
Figure PCTCN2017103286-appb-000080
Figure PCTCN2017103286-appb-000081
Figure PCTCN2017103286-appb-000081
Figure PCTCN2017103286-appb-000082
Figure PCTCN2017103286-appb-000082
Figure PCTCN2017103286-appb-000083
Figure PCTCN2017103286-appb-000083
Figure PCTCN2017103286-appb-000084
Figure PCTCN2017103286-appb-000084
Figure PCTCN2017103286-appb-000085
Figure PCTCN2017103286-appb-000085
Figure PCTCN2017103286-appb-000086
Figure PCTCN2017103286-appb-000086
Figure PCTCN2017103286-appb-000087
Figure PCTCN2017103286-appb-000087
Figure PCTCN2017103286-appb-000088
Figure PCTCN2017103286-appb-000088
其中,among them,
Figure PCTCN2017103286-appb-000089
Figure PCTCN2017103286-appb-000089
Figure PCTCN2017103286-appb-000090
Figure PCTCN2017103286-appb-000090
Figure PCTCN2017103286-appb-000091
Figure PCTCN2017103286-appb-000091
将关节角速度
Figure PCTCN2017103286-appb-000092
表示为冗余度向量(k5,k6,k7)的函数:
Joint angular velocity
Figure PCTCN2017103286-appb-000092
A function expressed as a redundancy vector (k 5 , k 6 , k 7 ):
Figure PCTCN2017103286-appb-000093
Figure PCTCN2017103286-appb-000093
将式(28)代入式(8)中,可以得到关节角度θ以冗余度向量(k5,k6,k7)为自变量的函数:Substituting equation (28) into equation (8), the joint angle θ can be obtained as a function of the redundancy vector (k 5 , k 6 , k 7 ) as an independent variable:
Figure PCTCN2017103286-appb-000094
Figure PCTCN2017103286-appb-000094
其中,θ0是当前关节角度,Δt是采样步长。Where θ 0 is the current joint angle and Δt is the sampling step size.
将式(28)代入式(10-1)、将式(29)代入式(10-2)可以得到关节运动幅度最小以及机械臂受对应的齿轮间隙影响最小的方程分别为: Substituting equation (28) into equation (10-1) and equation (29) into equation (10-2) can obtain the equations with the smallest joint motion and the minimum influence of the mechanical arm on the corresponding gear gap:
Figure PCTCN2017103286-appb-000095
Figure PCTCN2017103286-appb-000095
Figure PCTCN2017103286-appb-000096
Figure PCTCN2017103286-appb-000096
其中,Hi(θ(k5,k6,k7))是连杆变换矩阵0Ti(θ(k5,k6,k7))中第二行第四列对应的元素:Where H i (θ(k 5 , k 6 , k 7 )) is the element corresponding to the second row and the fourth column of the link transformation matrix 0 T i (θ(k 5 , k 6 , k 7 )):
0Ti(θ(k5,k6,k7))=0T1(d11T22)…i-1Tii),(i=1,…,9)    (32) 0 T i (θ(k 5 , k 6 , k 7 ))= 0 T 1 (d 11 T 22 )... i-1 T ii ), (i=1,..., 9) (32)
其中,i-1Tii)是Denavit-Hartenberg矩阵:Where i-1 T ii ) is the Denavit-Hartenberg matrix:
Figure PCTCN2017103286-appb-000097
Figure PCTCN2017103286-appb-000097
其中,di、αi和ai是Denavit-Hartenberg参数,其具体数值如表1所示。Wherein, d i , α i and a i are Denavit-Hartenberg parameters, and specific values thereof are shown in Table 1.
将式(28)代入式(11-1)、将式(29)代入式(11-2)可以得到关节角速度限制以及关节角度限制的约束条件为:Substituting equation (28) into equation (11-1) and substituting equation (29) into equation (11-2) can obtain joint angular velocity limits and joint angle constraints.
Figure PCTCN2017103286-appb-000098
Figure PCTCN2017103286-appb-000098
由式(33)可得:From formula (33):
Figure PCTCN2017103286-appb-000099
Figure PCTCN2017103286-appb-000099
由式(34)可以进一步得到式(12)中的矩阵A和向量b如下:From the equation (34), the matrix A and the vector b in the equation (12) can be further obtained as follows:
Figure PCTCN2017103286-appb-000100
Figure PCTCN2017103286-appb-000100
Figure PCTCN2017103286-appb-000101
Figure PCTCN2017103286-appb-000101
从而,多目标优化问题可以转化为:Thus, multi-objective optimization problems can be translated into:
Figure PCTCN2017103286-appb-000102
Figure PCTCN2017103286-appb-000102
由式(30)、(31)和(35)可以求得冗余度向量(k5,k6,k7)的最优值,将冗余度向量(k5,k6,k7)的最优值代入式(28)和式(29)即可得到冗余机械臂各关节在下一采样点的角速度和角度,从而得到冗余机械臂各关节的位置和速度。By the formula (30), (31) and (35) can be determined redundancy vector (k 5, k 6, k 7) the optimum value, the redundancy vector (k 5, k 6, k 7) The optimal value is substituted into equations (28) and (29) to obtain the angular velocity and angle of the joints of the redundant manipulator at the next sampling point, thereby obtaining the position and velocity of each joint of the redundant manipulator.
本发明实施例二提供的冗余机械臂的控制方法,获取冗余机械臂末端在目标点的速度,根据冗余机械臂末端在目标点的速度将冗余机械臂的关节角、关节角速度和关节角加速度映射到冗余度空间,根据逆运动学方程确定冗余机械臂从当前点移动到目标点的运动轨迹对应的优化目标函数和约束条件,以冗余度空间向量为自变量建立运动轨迹对应的方程,通过辅助向量和单目标优化算法求解轨迹函数对应的方程,得到运动轨迹对应的冗余机械臂中各关节的位置和速度。本实施例通过采用上述技术方案,建立多目标约束优化问题的数学模块,可以为求解多目标优化问题提供理论基础,在得到满足设计目标的全局最优解的前提下,减少冗余机械臂多目标优化问题的搜索空间,避免多目标求解过程中维度爆炸问题的发生,简化冗余机械臂控制过程中所需的计算量,提高冗余机械臂的反应速度。The control method of the redundant mechanical arm provided by the second embodiment of the present invention obtains the speed of the end of the redundant mechanical arm at the target point, and the joint angle, the joint angular velocity of the redundant mechanical arm and the speed of the end of the redundant mechanical arm at the target point The joint angular acceleration is mapped to the redundancy space. According to the inverse kinematics equation, the optimal objective function and constraints corresponding to the motion trajectory of the redundant manipulator moving from the current point to the target point are determined, and the redundant space vector is used as the independent variable to establish the motion. The equation corresponding to the trajectory is solved by the auxiliary vector and the single-objective optimization algorithm to solve the equation corresponding to the trajectory function, and the position and velocity of each joint in the redundant mechanical arm corresponding to the motion trajectory are obtained. In this embodiment, by using the above technical solution, a mathematical module for multi-objective constrained optimization problem is established, which can provide a theoretical basis for solving a multi-objective optimization problem, and reduce redundant manipulators under the premise of satisfying the global optimal solution of the design target. The search space of the target optimization problem avoids the occurrence of dimensional explosion problems in the multi-objective solution process, simplifies the calculation amount required in the redundant robot arm control process, and improves the reaction speed of the redundant mechanical arm.
实施例三 Embodiment 3
本发明实施例三提供一种冗余机械臂的控制装置。该装置可由硬件和/或软件实现,一般集成在用于控制机器人的控制模块中,可通过执行冗余机械臂的控制方法实现对冗余机械臂的控制。图4所示为本实施例提供的冗余机械臂的控制装置的结构框图,如图4所示,该装置包括:Embodiment 3 of the present invention provides a control device for a redundant mechanical arm. The device can be implemented by hardware and/or software, generally integrated in a control module for controlling the robot, and the control of the redundant manipulator can be realized by performing a control method of the redundant manipulator. 4 is a structural block diagram of a control device for a redundant mechanical arm according to the embodiment. As shown in FIG. 4, the device includes:
信息获取单元410,用于获取冗余机械臂所在的当前点信息和目标点信息;The information acquiring unit 410 is configured to acquire current point information and target point information where the redundant robot arm is located;
轨迹单元420,与所述信息获取单元相连,用于根据所述当前点信息和所述目标点信息确定所述冗余机械臂从当前点移动到目标点的运动轨迹对应的轨迹函数;The trajectory unit 420 is connected to the information acquiring unit, and configured to determine, according to the current point information and the target point information, a trajectory function corresponding to a motion trajectory of the redundant mechanical arm moving from a current point to a target point;
冗余度处理单元430,与所述轨迹单元相连,用于以冗余度空间向量为自变量建立所述轨迹函数对应的方程;根据目标接近法求解所述轨迹函数对应的方程,得到所述运动轨迹对应的所述冗余机械臂中各关节的位置和速度。The redundancy processing unit 430 is connected to the trajectory unit, and is configured to establish an equation corresponding to the trajectory function by using a redundancy space vector as an independent variable; and solving an equation corresponding to the trajectory function according to a target proximity method, to obtain the The position and velocity of each joint in the redundant robot arm corresponding to the motion trajectory.
本发明实施例三提供的冗余机械臂的控制装置,通过信息获取单元获取冗余机械臂的当前点信息和目标点信息,通过轨迹单元根据所获取的当前点信息和目标点信息确定冗余机械臂从当前点移动到目标点的运动轨迹对应的轨迹函数,通过冗余度处理单元以冗余度空间向量为自变量建立该轨迹函数对应的方程,通过目标接近法求解该方程以得到运动轨迹对应的冗余机械臂中各关节的位置和速度。本发明实施例通过采用上述技术方案,可以在得到满足设计目标的全局最优解的前提下,减少冗余机械臂多目标优化问题的搜索空间,避免多目标求解过程中维度爆炸问题的发生,简化冗余机械臂控制过程中所需的计算量,提高冗余机械臂的反应速度。The control device of the redundant mechanical arm provided by the third embodiment of the present invention obtains the current point information and the target point information of the redundant mechanical arm through the information acquiring unit, and determines the redundancy according to the acquired current point information and the target point information by the trajectory unit. The trajectory function corresponding to the motion trajectory of the robot arm moving from the current point to the target point is established by the redundancy processing unit with the redundancy space vector as the independent variable, and the equation is solved by the target approach method to obtain the motion. The position and velocity of each joint in the redundant robot arm corresponding to the trajectory. By adopting the above technical solution, the embodiment of the present invention can reduce the search space of the multi-objective optimization problem of the redundant manipulator under the premise of obtaining the global optimal solution satisfying the design goal, and avoid the occurrence of the dimensional explosion problem in the multi-objective solution process. Simplify the amount of calculations required during redundant manipulator control and increase the response speed of redundant manipulators.
进一步地,所述信息获取单元410具体用于:获取所述冗余机械臂末端在所述目标点的速度;根据所述冗余机械臂末端在所述目标点的速度将所述冗余机械臂的关节角、关节角速度和关节角加速度映射到冗余度空间。 Further, the information acquiring unit 410 is specifically configured to: acquire a speed of the end of the redundant mechanical arm at the target point; and the redundant machine according to a speed of the redundant mechanical arm end at the target point The joint angle of the arm, the joint angular velocity, and the joint angular acceleration are mapped to the redundancy space.
进一步地,所述获取所述冗余机械臂末端在所述当前点的速度包括:根据所述冗余机械臂末端在所述当前点的位姿和所述冗余机械臂末端在所述运动轨迹中的下一采样点的位姿确定所述冗余机械臂末端在所述下一采样点的速度。Further, the acquiring the speed of the end of the redundant mechanical arm at the current point comprises: according to the posture of the redundant mechanical arm end at the current point and the end of the redundant mechanical arm The pose of the next sample point in the trajectory determines the velocity of the end of the redundant robot arm at the next sample point.
进一步地,所述轨迹单元420具体用于根据逆运动学方程,确定所述冗余机械臂从当前点移动到目标点的运动轨迹对应的优化目标函数和约束条件。Further, the trajectory unit 420 is specifically configured to determine an optimization objective function and a constraint condition corresponding to a motion trajectory of the redundant mechanical arm moving from a current point to a target point according to an inverse kinematics equation.
进一步地,所述优化目标函数包括所述冗余机械臂的关节运动幅度最小或者所述冗余机械臂受对应的齿轮间隙影响最小;所述约束条件包括:将所述冗余机械臂的关节速度约束在预设的速度范围内或者将所述冗余机械臂的关节角度约束在预设的角度范围内。Further, the optimization objective function includes a minimum joint motion amplitude of the redundant mechanical arm or the redundant mechanical arm is minimally affected by a corresponding gear gap; the constraint condition includes: a joint of the redundant mechanical arm The speed constraint is within a preset speed range or the joint angle of the redundant robot arm is constrained within a preset angle range.
进一步地,所述根据目标接近法求解所述轨迹函数对应的方程,得到所述运动轨迹对应的所述冗余机械臂中各关节的位置和速度包括:在冗余度空间中,根据优化目标函数、与优化目标对应的最小目标值、对应的约束条件和预设的加权系数,通过辅助向量和单目标优化算法求解所述轨迹函数对应的方程,得到所述运动轨迹对应的所述冗余机械臂中各关节的位置和速度。Further, the solving the equation corresponding to the trajectory function according to the target approaching method, and obtaining the position and velocity of each joint in the redundant mechanical arm corresponding to the motion trajectory includes: in the redundancy space, according to the optimization target a function, a minimum target value corresponding to the optimization target, a corresponding constraint condition, and a preset weighting coefficient, and an equation corresponding to the trajectory function is solved by an auxiliary vector and a single target optimization algorithm to obtain the redundancy corresponding to the motion trajectory The position and speed of the joints in the robot arm.
进一步地,所述单目标优化算法包括:牛顿-欧拉算法、Nelder-Mead单纯形算法、内点算法、基因算法和/或Pattern Search算法。Further, the single target optimization algorithm includes: a Newton-Eulerian algorithm, a Nelder-Mead simplex algorithm, an interior point algorithm, a genetic algorithm, and/or a Pattern Search algorithm.
本实施例提供的冗余机械臂的控制装置可执行本发明任意实施例提供的冗余机械臂的控制方法,具备执行冗余机械臂的控制方法相应的功能模块和有益效果。未在本实施例中详尽描述的技术细节,可参见本发明任意实施例所提供的冗余机械臂的控制方法。The control device of the redundant mechanical arm provided by this embodiment can execute the control method of the redundant mechanical arm provided by any embodiment of the present invention, and has the corresponding functional modules and beneficial effects of the control method for executing the redundant mechanical arm. For technical details not fully described in this embodiment, reference may be made to the control method of the redundant robot arm provided by any embodiment of the present invention.
注意,上述仅为本发明的较佳实施例及所运用技术原理。本领域技术人员会理解,本发明不限于这里所述的特定实施例,对本领域技术人员来说能够进 行各种明显的变化、重新调整和替代而不会脱离本发明的保护范围。因此,虽然通过以上实施例对本发明进行了较为详细的说明,但是本发明不仅仅限于以上实施例,在不脱离本发明构思的情况下,还可以包括更多其他等效实施例,而本发明的范围由所附的权利要求范围决定。 Note that the above are only the preferred embodiments of the present invention and the technical principles applied thereto. Those skilled in the art will appreciate that the present invention is not limited to the specific embodiments described herein and can be made by those skilled in the art Various obvious changes, modifications, and substitutions are made without departing from the scope of the invention. Therefore, the present invention has been described in detail by the above embodiments, but the present invention is not limited to the above embodiments, and other equivalent embodiments may be included without departing from the inventive concept. The scope is determined by the scope of the appended claims.

Claims (10)

  1. 一种冗余机械臂的控制方法,其特征在于,包括:A method for controlling a redundant mechanical arm, comprising:
    获取冗余机械臂所在的当前点信息和目标点信息;Obtaining current point information and target point information where the redundant robot arm is located;
    根据所述当前点信息和所述目标点信息确定所述冗余机械臂从当前点移动到目标点的运动轨迹对应的轨迹函数;Determining, according to the current point information and the target point information, a trajectory function corresponding to a motion trajectory of the redundant robot arm moving from a current point to a target point;
    以冗余度空间向量为自变量建立所述轨迹函数对应的方程;Establishing an equation corresponding to the trajectory function by using a redundancy space vector as an independent variable;
    根据目标接近法求解所述轨迹函数对应的方程,得到所述运动轨迹对应的所述冗余机械臂中各关节的位置和速度。The equation corresponding to the trajectory function is solved according to the target approaching method, and the position and velocity of each joint in the redundant mechanical arm corresponding to the motion trajectory are obtained.
  2. 根据权利要求1所述的冗余机械臂的控制方法,其特征在于,所述获取冗余机械臂所在的当前点信息和目标点信息包括:The method for controlling a redundant robot arm according to claim 1, wherein the acquiring current point information and target point information of the redundant robot arm comprises:
    获取所述冗余机械臂末端在所述目标点的速度;Obtaining a velocity of the end of the redundant mechanical arm at the target point;
    根据所述冗余机械臂末端在所述目标点的速度将所述冗余机械臂的关节角、关节角速度和关节角加速度映射到冗余度空间。The joint angle, joint angular velocity, and joint angular acceleration of the redundant robot arm are mapped to the redundancy space according to the speed of the redundant robot arm end at the target point.
  3. 根据权利要求2所述的冗余机械臂的控制方法,其特征在于,所述获取所述冗余机械臂末端在所述当前点的速度包括:The control method of the redundant robot arm according to claim 2, wherein the obtaining the speed of the end of the redundant mechanical arm at the current point comprises:
    根据所述冗余机械臂末端在所述当前点的位姿和所述冗余机械臂末端在所述运动轨迹中的下一采样点的位姿确定所述冗余机械臂末端在所述下一采样点的速度。Determining the end of the redundant manipulator according to the pose of the end of the redundant manipulator at the current point and the pose of the next sample point of the end of the redundant manipulator in the motion trajectory The speed of a sample point.
  4. 根据权利要求1或2所述的冗余机械臂的控制方法,其特征在于,所述根据所述当前点信息和所述目标点信息确定所述冗余机械臂从当前点移动到目标点的运动轨迹对应的轨迹函数包括:The control method of the redundant robot arm according to claim 1 or 2, wherein the determining, by the current point information and the target point information, that the redundant robot arm moves from a current point to a target point The trajectory function corresponding to the motion trajectory includes:
    根据逆运动学方程,确定所述冗余机械臂从当前点移动到目标点的运动轨迹对应的优化目标函数和约束条件。And determining, according to the inverse kinematics equation, an optimization objective function and a constraint condition corresponding to the motion trajectory of the redundant mechanical arm moving from the current point to the target point.
  5. 根据权利要求4所述的冗余机械臂的控制方法,其特征在于,所述优化 目标函数包括所述冗余机械臂的关节运动幅度最小或者所述冗余机械臂受对应的齿轮间隙影响最小;所述约束条件包括:将所述冗余机械臂的关节速度约束在预设的速度范围内或者将所述冗余机械臂的关节角度约束在预设的角度范围内。A method of controlling a redundant robot arm according to claim 4, wherein said optimizing The objective function includes a minimum degree of joint motion of the redundant mechanical arm or the redundant mechanical arm is minimally affected by a corresponding gear gap; the constraint condition includes constraining the joint speed of the redundant mechanical arm to a preset The joint angle of the redundant manipulator is constrained within a predetermined range of angles within the speed range.
  6. 根据权利要求5所述的冗余机械臂的控制方法,其特征在于,所述根据目标接近法求解所述轨迹函数对应的方程,得到所述运动轨迹对应的所述冗余机械臂中各关节的位置和速度包括:The control method of the redundant robot arm according to claim 5, wherein the equation corresponding to the trajectory function is solved according to a target approach method, and each joint in the redundant mechanical arm corresponding to the motion trajectory is obtained. The location and speed include:
    在冗余度空间中,根据优化目标函数、与优化目标对应的最小目标值、对应的约束条件和预设的加权系数,通过辅助向量和单目标优化算法求解所述轨迹函数对应的方程,得到所述运动轨迹对应的所述冗余机械臂中各关节的位置和速度。In the redundancy space, according to the optimization objective function, the minimum target value corresponding to the optimization target, the corresponding constraint condition and the preset weighting coefficient, the equation corresponding to the trajectory function is solved by the auxiliary vector and the single target optimization algorithm, and obtained The motion trajectory corresponds to the position and velocity of each joint in the redundant robot arm.
  7. 根据权利要求6所述的冗余机械臂的控制方法,其特征在于,所述单目标优化算法包括:牛顿-欧拉算法、Nelder-Mead单纯形算法、内点算法、基因算法和/或Pattern Search算法。The control method of a redundant robot arm according to claim 6, wherein the single target optimization algorithm comprises: a Newton-Eulerian algorithm, a Nelder-Mead simplex algorithm, an interior point algorithm, a genetic algorithm, and/or a Pattern. Search algorithm.
  8. 一种冗余机械臂的控制装置,其特征在于,包括:A control device for a redundant mechanical arm, comprising:
    信息获取单元,用于获取冗余机械臂所在的当前点信息和目标点信息;An information acquiring unit, configured to acquire current point information and target point information where the redundant robot arm is located;
    轨迹单元,与所述信息获取单元相连,用于根据所述当前点信息和所述目标点信息确定所述冗余机械臂从当前点移动到目标点的运动轨迹对应的轨迹函数;a trajectory unit, configured to be connected to the information acquiring unit, configured to determine, according to the current point information and the target point information, a trajectory function corresponding to a motion trajectory of the redundant mechanical arm moving from a current point to a target point;
    冗余度处理单元,与所述轨迹单元相连,用于以冗余度空间向量为自变量建立所述轨迹函数对应的方程;根据目标接近法求解所述轨迹函数对应的方程,得到所述运动轨迹对应的所述冗余机械臂中各关节的位置和速度。a redundancy processing unit, connected to the trajectory unit, configured to establish an equation corresponding to the trajectory function by using a redundancy space vector as an independent variable; and solving an equation corresponding to the trajectory function according to a target approach method, to obtain the motion The position and velocity of each joint in the redundant robot arm corresponding to the trajectory.
  9. 根据权利要求8所述的冗余机械臂的控制装置,其特征在于,所述信息 获取单元具体用于:获取所述冗余机械臂末端在所述当前点的速度;根据所述冗余机械臂末端在所述当前点的速度将所述冗余机械臂的关节角、关节角速度和关节角加速度映射到冗余度空间。A control device for a redundant robot arm according to claim 8, wherein said information The acquiring unit is specifically configured to: acquire a velocity of the end of the redundant mechanical arm at the current point; and combine an articulated angle and an articulated angular velocity of the redundant mechanical arm according to a speed of the redundant mechanical arm end at the current point And joint angular acceleration is mapped to the redundancy space.
  10. 根据权利要求8或9所述的冗余机械臂的控制装置,其特征在于,所述轨迹单元具体用于根据逆运动学方程,确定所述冗余机械臂从当前点移动到目标点的运动轨迹对应的优化目标函数和约束条件。 The control device of the redundant robot arm according to claim 8 or 9, wherein the trajectory unit is specifically configured to determine the movement of the redundant mechanical arm from the current point to the target point according to an inverse kinematics equation. The optimization objective function and constraints corresponding to the trajectory.
PCT/CN2017/103286 2016-12-16 2017-09-25 Method and device for controlling redundant robot arm WO2018107851A1 (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN201611169869.4 2016-12-16
CN201611169869.4A CN106625666B (en) 2016-12-16 2016-12-16 The control method and device of redundant mechanical arm

Publications (1)

Publication Number Publication Date
WO2018107851A1 true WO2018107851A1 (en) 2018-06-21

Family

ID=58822194

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2017/103286 WO2018107851A1 (en) 2016-12-16 2017-09-25 Method and device for controlling redundant robot arm

Country Status (2)

Country Link
CN (1) CN106625666B (en)
WO (1) WO2018107851A1 (en)

Cited By (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111300414A (en) * 2020-03-06 2020-06-19 陕西理工大学 Dual-criterion redundant mechanical arm self-movement planning method
CN111399514A (en) * 2020-03-30 2020-07-10 浙江钱江机器人有限公司 Robot time optimal trajectory planning method
CN112706163A (en) * 2020-12-10 2021-04-27 中山大学 Mechanical arm motion control method, device, equipment and medium
CN112914601A (en) * 2021-01-19 2021-06-08 深圳市德力凯医疗设备股份有限公司 Obstacle avoidance method and device for mechanical arm, storage medium and ultrasonic equipment
CN113084803A (en) * 2021-03-31 2021-07-09 信阳师范学院 Multi-task control method for redundant manipulator based on multilayer structure
CN113276121A (en) * 2021-05-31 2021-08-20 华南理工大学 Redundant manipulator moving obstacle avoidance method based on quadratic programming
CN113843793A (en) * 2021-09-21 2021-12-28 兰州大学 Mobile redundant mechanical arm model prediction control method with obstacle avoidance function
CN113878578A (en) * 2021-09-30 2022-01-04 上海景吾智能科技有限公司 Dynamic self-adaptive positioning method and system suitable for composite robot
CN114211500A (en) * 2021-12-31 2022-03-22 华南理工大学 Self-adaptive fuzzy neural network planning method
CN114211492A (en) * 2021-12-22 2022-03-22 武汉鼎元同立科技有限公司 Optimal trajectory planning method of multi-degree-of-freedom mechanical arm based on model
CN114227687A (en) * 2021-12-28 2022-03-25 深圳市优必选科技股份有限公司 Robot control method and device, terminal equipment and storage medium
CN114329986A (en) * 2021-12-30 2022-04-12 华中科技大学 Redundant mechanical arm anthropomorphic kinematics inverse solution solving method based on human muscle bone model
CN114670190A (en) * 2022-03-08 2022-06-28 西北工业大学 Redundant mechanical arm inverse kinematics method based on analytical numerical value mixing method
CN114734441A (en) * 2022-04-15 2022-07-12 北京邮电大学 Method for optimizing motion capability of mechanical arm in failure fault space of joint part
CN115462908A (en) * 2022-09-16 2022-12-13 哈尔滨工业大学 Main manipulator structure of minimally invasive surgery robot
WO2023024278A1 (en) * 2021-08-24 2023-03-02 深圳市优必选科技股份有限公司 Robot joint pose optimization method, robot control method, and robot
CN116985147A (en) * 2023-09-27 2023-11-03 睿尔曼智能科技(北京)有限公司 Mechanical arm inverse kinematics solving method and device

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106625666B (en) * 2016-12-16 2019-03-01 广州视源电子科技股份有限公司 The control method and device of redundant mechanical arm
JP6931457B2 (en) 2017-07-14 2021-09-08 オムロン株式会社 Motion generation method, motion generator, system and computer program
CN107703756B (en) * 2017-11-03 2021-03-02 广州视源电子科技股份有限公司 Kinetic model parameter identification method and device, computer equipment and storage medium
CN108537404B (en) * 2018-03-06 2021-10-22 中国人民解放军63920部队 Extraterrestrial celestial body detection sampling area collectability assessment method, medium and equipment
CN108638055B (en) * 2018-04-11 2020-07-14 北京控制工程研究所 Autonomous obstacle avoidance planning method for seven-degree-of-freedom space manipulator
CN111345894B (en) * 2018-12-21 2022-08-02 上海微创医疗机器人(集团)股份有限公司 Mechanical arm and surgical robot
CN110125942B (en) * 2019-06-21 2022-06-10 上海工程技术大学 Plane trajectory tracking method for mobile snakelike mechanical arm
CN113001537B (en) * 2019-12-20 2022-08-02 深圳市优必选科技股份有限公司 Mechanical arm control method, mechanical arm control device and terminal equipment
CN113378349B (en) * 2021-03-25 2022-05-20 北京航空航天大学 Numerical stabilization algorithm for S-R-S structure seven-degree-of-freedom mechanical arm inverse kinematics analytic solution
CN114800477B (en) * 2022-05-25 2023-03-14 华东交通大学 Minimum flow-based redundant hydraulic mechanical arm motion planning method
CN117301038A (en) * 2022-06-22 2023-12-29 瑞龙诺赋(上海)医疗科技有限公司 Mechanical arm adjusting method and device, electronic equipment and storage medium

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS61271508A (en) * 1985-05-27 1986-12-01 Daikin Ind Ltd Action controller for multi-joint robot
JP2012130999A (en) * 2010-12-22 2012-07-12 Kawasaki Heavy Ind Ltd Control device and method for controlling robot arm
CN103147577A (en) * 2013-02-27 2013-06-12 中联重科股份有限公司 Control method, equipment, system and construction machinery for multi-joint mechanical arm support
CN103984230A (en) * 2014-05-09 2014-08-13 大连大学 Zero-disturbance optimization control method for base of space manipulator
US20140358283A1 (en) * 2011-09-21 2014-12-04 Seiko Epson Corporation Robot and robot control method
CN104526695A (en) * 2014-12-01 2015-04-22 北京邮电大学 Space manipulator track planning method for minimizing base seat collision disturbance
CN106625666A (en) * 2016-12-16 2017-05-10 广州视源电子科技股份有限公司 Method and device for controlling redundant mechanical arm

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2569412B2 (en) * 1991-04-30 1997-01-08 工業技術院長 Compliance Control Method for Redundant Robot Manipulator
JP3379540B2 (en) * 1991-10-31 2003-02-24 株式会社安川電機 Arm control method for redundant arm robot
CN103492133B (en) * 2011-04-19 2016-04-13 Abb研究有限公司 There is the industrial robot of motion redundancy arm and the method for controlling this robot
CN103009389B (en) * 2012-11-30 2015-07-08 北京控制工程研究所 Track planning method of redundant space mechanical arm for on-track catching
CN104331547B (en) * 2014-10-23 2017-05-10 北京控制工程研究所 Space mechanical arm structure parameter optimization method based on operability
CN105183009B (en) * 2015-10-15 2017-11-21 哈尔滨工程大学 A kind of redundant mechanical arm method for controlling trajectory

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS61271508A (en) * 1985-05-27 1986-12-01 Daikin Ind Ltd Action controller for multi-joint robot
JP2012130999A (en) * 2010-12-22 2012-07-12 Kawasaki Heavy Ind Ltd Control device and method for controlling robot arm
US20140358283A1 (en) * 2011-09-21 2014-12-04 Seiko Epson Corporation Robot and robot control method
CN103147577A (en) * 2013-02-27 2013-06-12 中联重科股份有限公司 Control method, equipment, system and construction machinery for multi-joint mechanical arm support
CN103984230A (en) * 2014-05-09 2014-08-13 大连大学 Zero-disturbance optimization control method for base of space manipulator
CN104526695A (en) * 2014-12-01 2015-04-22 北京邮电大学 Space manipulator track planning method for minimizing base seat collision disturbance
CN106625666A (en) * 2016-12-16 2017-05-10 广州视源电子科技股份有限公司 Method and device for controlling redundant mechanical arm

Cited By (28)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111300414A (en) * 2020-03-06 2020-06-19 陕西理工大学 Dual-criterion redundant mechanical arm self-movement planning method
CN111300414B (en) * 2020-03-06 2022-07-15 陕西理工大学 Dual-criterion redundant mechanical arm self-movement planning method
CN111399514A (en) * 2020-03-30 2020-07-10 浙江钱江机器人有限公司 Robot time optimal trajectory planning method
CN112706163B (en) * 2020-12-10 2022-03-04 中山大学 Mechanical arm motion control method, device, equipment and medium
CN112706163A (en) * 2020-12-10 2021-04-27 中山大学 Mechanical arm motion control method, device, equipment and medium
CN112914601B (en) * 2021-01-19 2024-04-02 深圳市德力凯医疗设备股份有限公司 Obstacle avoidance method and device for mechanical arm, storage medium and ultrasonic equipment
CN112914601A (en) * 2021-01-19 2021-06-08 深圳市德力凯医疗设备股份有限公司 Obstacle avoidance method and device for mechanical arm, storage medium and ultrasonic equipment
CN113084803A (en) * 2021-03-31 2021-07-09 信阳师范学院 Multi-task control method for redundant manipulator based on multilayer structure
CN113276121A (en) * 2021-05-31 2021-08-20 华南理工大学 Redundant manipulator moving obstacle avoidance method based on quadratic programming
CN113276121B (en) * 2021-05-31 2022-08-09 华南理工大学 Redundant manipulator moving obstacle avoiding method based on quadratic programming
WO2023024278A1 (en) * 2021-08-24 2023-03-02 深圳市优必选科技股份有限公司 Robot joint pose optimization method, robot control method, and robot
CN113843793A (en) * 2021-09-21 2021-12-28 兰州大学 Mobile redundant mechanical arm model prediction control method with obstacle avoidance function
CN113878578A (en) * 2021-09-30 2022-01-04 上海景吾智能科技有限公司 Dynamic self-adaptive positioning method and system suitable for composite robot
CN113878578B (en) * 2021-09-30 2024-01-16 上海景吾智能科技有限公司 Dynamic self-adaptive positioning method and system suitable for composite robot
CN114211492A (en) * 2021-12-22 2022-03-22 武汉鼎元同立科技有限公司 Optimal trajectory planning method of multi-degree-of-freedom mechanical arm based on model
CN114211492B (en) * 2021-12-22 2023-07-07 武汉鼎元同立科技有限公司 Optimal track planning method of multi-degree-of-freedom mechanical arm based on model
CN114227687B (en) * 2021-12-28 2023-08-15 深圳市优必选科技股份有限公司 Robot control method and device, terminal equipment and storage medium
CN114227687A (en) * 2021-12-28 2022-03-25 深圳市优必选科技股份有限公司 Robot control method and device, terminal equipment and storage medium
CN114329986A (en) * 2021-12-30 2022-04-12 华中科技大学 Redundant mechanical arm anthropomorphic kinematics inverse solution solving method based on human muscle bone model
CN114211500B (en) * 2021-12-31 2023-05-30 华南理工大学 Planning method of self-adaptive fuzzy neural network
CN114211500A (en) * 2021-12-31 2022-03-22 华南理工大学 Self-adaptive fuzzy neural network planning method
CN114670190A (en) * 2022-03-08 2022-06-28 西北工业大学 Redundant mechanical arm inverse kinematics method based on analytical numerical value mixing method
CN114670190B (en) * 2022-03-08 2023-10-24 西北工业大学 Redundant mechanical arm inverse kinematics method based on analysis numerical mixing method
CN114734441A (en) * 2022-04-15 2022-07-12 北京邮电大学 Method for optimizing motion capability of mechanical arm in failure fault space of joint part
CN114734441B (en) * 2022-04-15 2023-11-24 北京邮电大学 Joint part failure fault space mechanical arm movement capacity optimization method
CN115462908A (en) * 2022-09-16 2022-12-13 哈尔滨工业大学 Main manipulator structure of minimally invasive surgery robot
CN116985147A (en) * 2023-09-27 2023-11-03 睿尔曼智能科技(北京)有限公司 Mechanical arm inverse kinematics solving method and device
CN116985147B (en) * 2023-09-27 2023-12-22 睿尔曼智能科技(北京)有限公司 Mechanical arm inverse kinematics solving method and device

Also Published As

Publication number Publication date
CN106625666B (en) 2019-03-01
CN106625666A (en) 2017-05-10

Similar Documents

Publication Publication Date Title
WO2018107851A1 (en) Method and device for controlling redundant robot arm
US11845186B2 (en) Inverse kinematics solving method for redundant robot and redundant robot and computer readable storage medium using the same
CN110275436B (en) RBF neural network self-adaptive control method of multi-single-arm manipulator
Ananthanarayanan et al. Real-time Inverse Kinematics of (2n+ 1) DOF hyper-redundant manipulator arm via a combined numerical and analytical approach
JP5261495B2 (en) Real-time self-collision and obstacle avoidance using weight matrix
CN110682286B (en) Real-time obstacle avoidance method for cooperative robot
Hasan et al. Artificial neural network-based kinematics Jacobian solution for serial manipulator passing through singular configurations
US8406989B1 (en) Method for adaptive obstacle avoidance for articulated redundant robot arm
CN110370256B (en) Robot and path planning method, device and controller thereof
Tang et al. Swarm robots search for multiple targets based on an improved grouping strategy
JPH06314103A (en) Controller and passive sensing device
WO2017132905A1 (en) Method and apparatus for controlling motion system
Kouabon et al. A Learning Framework to inverse kinematics of high DOF redundant manipulators
CN111975771A (en) Mechanical arm motion planning method based on deviation redefinition neural network
Navarro-Alarcon et al. A Lyapunov-stable adaptive method to approximate sensorimotor models for sensor-based control
Waltersson et al. Planning and control for cable-routing with dual-arm robot
Chembuly et al. An efficient approach for inverse kinematics and redundancy resolution of spatial redundant robots for cluttered environment
Zhao et al. A learning-based multiscale modelling approach to real-time serial manipulator kinematics simulation
Benzaoui et al. Redundant robot manipulator control with obstacle avoidance using extended jacobian method
Banga et al. Modeling and simulation of robotic arm movement using soft computing
Wang et al. Whole body control of a dual-arm mobile robot using a virtual kinematic chain
Kawaharazuka et al. Dynamic task control method of a flexible manipulator using a deep recurrent neural network
Al-Faiz et al. Human arm inverse kinematic solution based geometric relations and optimization algorithm
CN115533920A (en) Collaborative planning method and system for solving inverse kinematics of rope-driven mechanical arm and computer storage medium
Crenganis et al. Inverse kinematics of a 7 DOF manipulator using adaptive neuro-fuzzy inference systems

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 17881861

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

32PN Ep: public notification in the ep bulletin as address of the adressee cannot be established

Free format text: NOTING OF LOSS OF RIGHTS PURSUANT TO RULE 112(1) EPC (EPO FORM 1205A DATED 22.10.2019)

122 Ep: pct application non-entry in european phase

Ref document number: 17881861

Country of ref document: EP

Kind code of ref document: A1