[go: up one dir, main page]

CN109760051B - Rope length change determination method for rope-driven super-redundancy degree of freedom robot - Google Patents

Rope length change determination method for rope-driven super-redundancy degree of freedom robot Download PDF

Info

Publication number
CN109760051B
CN109760051B CN201910041486.6A CN201910041486A CN109760051B CN 109760051 B CN109760051 B CN 109760051B CN 201910041486 A CN201910041486 A CN 201910041486A CN 109760051 B CN109760051 B CN 109760051B
Authority
CN
China
Prior art keywords
rope
segment
force
vertebral
vector
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201910041486.6A
Other languages
Chinese (zh)
Other versions
CN109760051A (en
Inventor
谷海宇
魏承
王萍萍
赵阳
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Harbin Institute of Technology Shenzhen
Original Assignee
Harbin Institute of Technology Shenzhen
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 Harbin Institute of Technology Shenzhen filed Critical Harbin Institute of Technology Shenzhen
Priority to CN201910041486.6A priority Critical patent/CN109760051B/en
Publication of CN109760051A publication Critical patent/CN109760051A/en
Application granted granted Critical
Publication of CN109760051B publication Critical patent/CN109760051B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Manipulator (AREA)

Abstract

本发明提供一种绳驱超冗余自由度机器人的绳索长度变化确定方法,属于机器人控制技术领域。本发明首先在绳驱超冗余自由度机器人的绳驱串并联机构上建立坐标系;定义绳段矢量并由绳段矢量求得绳索力作用在绳索导向板上的力;然后利用绳索力作用在绳索导向板上的力得到绳索力的等效旋量;通过绳索力与等效力矩之间映射关系得到雅可比矩阵;再结合转动自由度的角速度与雅可比矩阵得到绳长变化速率,最后将得到的绳长变化速率进行积分得到绳索长度的变化。本发明解决了现有绳驱超冗余自由度机器人的控制技术误差较大的问题。本发明可用于绳驱超冗余自由度机器人的控制技术。

Figure 201910041486

The invention provides a method for determining the change of rope length of a rope-driven super-redundant degree-of-freedom robot, which belongs to the technical field of robot control. The invention first establishes a coordinate system on the rope-driven series-parallel mechanism of the rope-driven super-redundant degree-of-freedom robot; defines the rope segment vector and obtains the force acting on the rope guide plate by the rope force from the rope segment vector; and then uses the rope force to act on the rope guide plate. The equivalent screw of the rope force is obtained from the force on the rope guide plate; the Jacobian matrix is obtained through the mapping relationship between the rope force and the equivalent moment; the rate of change of the rope length is obtained by combining the angular velocity of the rotational degree of freedom and the Jacobian matrix, and finally Integrate the resulting rate of change in rope length to obtain the change in rope length. The invention solves the problem that the control technology error of the existing rope-driven super-redundant degree-of-freedom robot is relatively large. The invention can be used for the control technology of the rope-driven super-redundant degree-of-freedom robot.

Figure 201910041486

Description

一种绳驱超冗余自由度机器人的绳索长度变化确定方法A Rope Length Variation Determination Method for Rope-Driven Super-Redundant Degree-of-Freedom Robots

技术领域technical field

本发明涉及一种绳驱超冗余自由度机器人的绳索长度变化确定方法,属于机器人控制技术领域。The invention relates to a method for determining the change of rope length of a rope-driven super-redundant degree-of-freedom robot, which belongs to the technical field of robot control.

背景技术Background technique

绳驱超冗余自由度机器人:Rope-driven super-redundant degrees of freedom robot:

绳驱超冗余自由度机器人是使用绳索作为传动介质的串并联机构。此类机器人由密集的运动关节串联而成,并使用绳索并联驱动,自由度数目众多,运动灵活,在狭窄受限的环境中具有极强的运动能力。由于绳索只能承受单向力,对机构只能起到单方向约束能力,因此绳索数目一般大于自由度数目。不同于传统机器人的关节运动控制方式,此类机器人需要对绳索长度和速度进行控制。而对此进行控制的前提是求得完成期望运动时,绳索需要满足运动状态(绳长和绳速)。此处所述的超冗余自由度机器人为由两自由度万向节串联而成的机器人。The rope-driven super-redundant degree-of-freedom robot is a series-parallel mechanism using rope as the transmission medium. This kind of robot is composed of dense kinematic joints connected in series and driven by ropes in parallel, with a large number of degrees of freedom, flexible movement, and strong movement ability in a narrow and restricted environment. Since the rope can only bear one-way force and can only restrain the mechanism in one direction, the number of ropes is generally greater than the number of degrees of freedom. Different from the joint motion control method of traditional robots, such robots need to control the length and speed of the rope. The premise of this control is to find that the rope needs to meet the motion state (rope length and rope speed) to complete the desired motion. The super-redundant degree-of-freedom robot described here is a robot composed of two-degree-of-freedom universal joints connected in series.

绳驱机器人逆运动学:Inverse kinematics of a rope-driven robot:

绳驱机器人逆运动学包含两部分:由操作空间运动状态(末端运动)求解关节空间运动状态和由关节运动状态求解绳索运动状态。其中,由操作空间运动状态求解关节空间运动状态的方法与传统冗余机器人的逆运动学求解方法相同,即通过速度空间的雅可比矩阵映射求解或通过正运动学方程求数值解。The inverse kinematics of the rope-driven robot consists of two parts: solving the joint space motion state from the operation space motion state (end motion) and solving the rope motion state from the joint motion state. Among them, the method of solving the joint space motion state from the operation space motion state is the same as the inverse kinematics solution method of the traditional redundant robot, that is, the solution by the Jacobian matrix mapping of the velocity space or the numerical solution by the forward kinematics equation.

而由关节运动状态求解绳索运动状态的部分,传统的求解方法为绳索多段累积方法,建立单关节处的绳索长度与关节角度的映射关系,利用数值方法求解,然后对所有绳段进行叠加,求得绳索的长度,该方法的优点是建模简单,但求解复杂。其中涉及到求解非线性方程组,存在多解现象,以及多段叠加过程中的误差累积。For the part of solving the rope motion state from the joint motion state, the traditional solution method is the rope multi-segment accumulation method. The mapping relationship between the rope length and the joint angle at a single joint is established, and the numerical method is used to solve it. Then all rope segments are superimposed to find The length of the rope is obtained. The advantage of this method is that the modeling is simple, but the solution is complicated. It involves solving nonlinear equations, the phenomenon of multiple solutions, and the accumulation of errors in the multi-segment stacking process.

发明内容SUMMARY OF THE INVENTION

本发明为解决现有绳驱超冗余自由度机器人的控制技术误差较大的问题,提供了一种绳驱超冗余自由度机器人的绳索长度变化确定方法。In order to solve the problem that the control technology error of the existing rope-driven super-redundant degree-of-freedom robot is relatively large, the invention provides a method for determining the change of the rope length of the rope-driven super-redundant degree-of-freedom robot.

本发明所述一种绳驱超冗余自由度机器人的绳索长度变化确定方法,通过以下技术方案实现:The method for determining the change of the rope length of a rope-driven super-redundant degree-of-freedom robot according to the present invention is realized by the following technical solutions:

步骤一、在绳驱超冗余自由度机器人的绳驱串并联机构上建立坐标系;定义绳段矢量并由绳段矢量求得绳索力作用在绳索导向板上的力;Step 1. Establish a coordinate system on the rope-driven series-parallel mechanism of the rope-driven super-redundant degree-of-freedom robot; define the rope segment vector and obtain the force of the rope force acting on the rope guide plate from the rope segment vector;

步骤二、利用绳索力作用在绳索导向板上的力得到绳索力的等效旋量;Step 2: Obtain the equivalent screw of the rope force by using the force acting on the rope guide plate by the rope force;

步骤三、通过绳索力与关节等效力矩之间映射关系得到雅可比矩阵;Step 3: Obtain the Jacobian matrix through the mapping relationship between the rope force and the joint equivalent moment;

步骤四、结合转动自由度的角速度与雅可比矩阵得到绳长变化速率;Step 4: Combining the angular velocity of the rotational degrees of freedom and the Jacobian matrix to obtain the rate of change of the rope length;

步骤五、将步骤四中得到的绳长变化速率进行积分得到绳索长度的变化。Step 5: Integrate the rate of change of the rope length obtained in the step 4 to obtain the change of the rope length.

本发明最为突出的特点和显著的有益效果是:The most prominent feature and significant beneficial effect of the present invention are:

本发明所涉及的一种绳驱超冗余自由度机器人的绳索长度变化确定方法,从速度空间进行逆运动学求解,求得绳索长度与转动自由度(关节)速度对应的雅可比矩阵,进而积分得到绳索长度的变化量,进而可以对所述绳驱超冗余自由度机器人进行控制。本发明方法运算简单,求得的绳索长度的变化值唯一并且误差很小,因此也能够使得绳驱超冗余自由度机器人的控制精度提高,相比传统方法,本发明方法能够有效提高绳驱超冗余自由度机器人的控制精度约20%。The invention relates to a method for determining the change of rope length of a rope-driven super-redundant degree-of-freedom robot. Inverse kinematics is solved from the speed space to obtain the Jacobian matrix corresponding to the rope length and the rotational degree of freedom (joint) speed, and then The integral obtains the variation of the rope length, and then the rope-driven super-redundant degree-of-freedom robot can be controlled. The method of the invention is simple in operation, the obtained rope length has a unique change value and the error is small, so the control accuracy of the rope-driven super-redundant degree-of-freedom robot can also be improved. Compared with the traditional method, the method of the invention can effectively improve the rope drive The control accuracy of the hyper-redundant degree-of-freedom robot is about 20%.

附图说明Description of drawings

图1为绳驱超冗余自由度机器人的绳驱串并联机构结构示意图;Figure 1 is a schematic structural diagram of a rope-driven series-parallel mechanism of a rope-driven super-redundant degree-of-freedom robot;

图2为绳驱超冗余自由度机器人的绳驱串并联机构建立的各坐标系示意图;Figure 2 is a schematic diagram of each coordinate system established by the rope-driven series-parallel mechanism of the rope-driven super-redundant degree-of-freedom robot;

图3为本发明中椎节结构示意图;3 is a schematic diagram of a vertebral segment structure in the present invention;

图4为本发明流程图;Fig. 4 is the flow chart of the present invention;

1.基座,2、绳驱机械臂,21.椎节,211.万向节,212.上板,213.下板,22.绳索。1. Base, 2. Rope-driven robotic arm, 21. Vertebral segment, 211. Universal joint, 212. Upper plate, 213. Lower plate, 22. Rope.

具体实施方式Detailed ways

具体实施方式一:结合图1、图4对本实施方式进行说明,本实施方式给出的一种绳驱超冗余自由度机器人的绳索长度变化确定方法,具体包括以下步骤:Embodiment 1: This embodiment will be described with reference to FIG. 1 and FIG. 4. A method for determining the change of rope length of a rope-driven super-redundant degree of freedom robot provided in this embodiment specifically includes the following steps:

步骤一、在绳驱超冗余自由度机器人的绳驱串并联机构上建立坐标系;定义绳段矢量并由绳段矢量求得绳索力作用在绳索导向板上的力;Step 1. Establish a coordinate system on the rope-driven series-parallel mechanism of the rope-driven super-redundant degree-of-freedom robot; define the rope segment vector and obtain the force of the rope force acting on the rope guide plate from the rope segment vector;

步骤二、利用绳索力作用在绳索导向板上的力得到绳索力的等效旋量;Step 2: Obtain the equivalent screw of the rope force by using the force acting on the rope guide plate by the rope force;

步骤三、通过绳索力与关节等效力矩之间映射关系得到雅可比矩阵;Step 3: Obtain the Jacobian matrix through the mapping relationship between the rope force and the joint equivalent moment;

步骤四、结合转动自由度(关节)的角速度与雅可比矩阵得到绳长变化速率;Step 4. Combine the angular velocity of the rotational degrees of freedom (joints) and the Jacobian matrix to obtain the rope length change rate;

步骤五、将步骤四中得到的绳长变化速率进行积分得到绳索长度的变化。Step 5: Integrate the rate of change of the rope length obtained in the step 4 to obtain the change of the rope length.

具体实施方式二:本实施方式与具体实施方式一不同的是,所述步骤一具体包括以下过程:Embodiment 2: The difference between this embodiment and Embodiment 1 is that the step 1 specifically includes the following processes:

如图1、图2所示.,绳驱超冗余自由度机器人的绳驱串并联机构包括:基座1以及设置在基座1上的绳驱机械臂2;将所述绳驱机械臂2上任意一个万向节211以及该万向节211至下一万向节211之间的部分定义为一个椎节21;任意椎节21中,定义靠近万向节211的绳索导向板为下板213,远离万向节211的绳索导向板为上板212;上板212和下板213的边缘均匀设置有与绳索22数量相等的绳孔,每一绳索均依次穿过所有椎节21的下板213和上板212。As shown in Figures 1 and 2, the rope-driven series-parallel mechanism of the rope-driven super-redundant degree-of-freedom robot includes: a base 1 and a rope-driven mechanical arm 2 arranged on the base 1; 2. Any one of the above universal joints 211 and the part between the universal joint 211 and the next universal joint 211 are defined as a vertebral segment 21; in any vertebral segment 21, the rope guide plate close to the universal joint 211 is defined as the lower one. Plate 213, the rope guide plate away from the universal joint 211 is the upper plate 212; the edges of the upper plate 212 and the lower plate 213 are evenly provided with rope holes equal to the number of the ropes 22, and each rope passes through all the vertebral segments 21 in turn. Lower plate 213 and upper plate 212 .

如图2、图3所示,以椎节的下板中心Ci为原点建立椎节坐标系{Ci-xiyizi},i表示椎节的序号,i=1,2,...,I;I表示椎节的总数目;定义由下板中心Ci指向下板上第1个绳孔的方向为xi轴方向、垂直于下板的方向为zi轴方向、yi轴同时垂直于xi轴和zi轴;在万向节中心处建立D-H(Denavit德纳维和Hartenberg哈登伯格在1955年提出一种通用的方法)坐标系{On-xnynzn},n表示转动自由度的序号,n=1,2,...,N;N表示转动自由度的总数目,N=2I;也就是说在每一个万向节中心均具有两个方向的转动自由度;位于椎节i的上板上的绳孔在惯性坐标系{O-XYZ}中的矢量记为

Figure GDA0002242193010000031
其在椎节坐标系{Ci-xiyizi}中的矢量为
Figure GDA0002242193010000032
j表示绳索的序号,椎节上j绳索通过的绳孔的序号也为j,j=1,2,...,M;M为绳索的总数,N=3I;位于椎节i的下板上的绳孔在惯性坐标系{O-XYZ}中的矢量记为
Figure GDA0002242193010000033
其在椎节坐标系{Ci-xiyizi}中的矢量为
Figure GDA0002242193010000034
As shown in Figure 2 and Figure 3, the vertebral segment coordinate system {C i -x i y i z i } is established with the center C i of the lower plate of the vertebral segment as the origin, i represents the serial number of the vertebral segment, i=1,2, ...,I; I represents the total number of vertebral segments; the direction from the center C i of the lower plate to the first rope hole in the lower plate is defined as the xi axis direction, the direction perpendicular to the lower plate is the z i axis direction, The y i axis is perpendicular to both the x i axis and the z i axis; the DH (a general method proposed by Denavit and Hartenberg Hardenberg in 1955) is established at the center of the universal joint. The coordinate system {O n -x n y n z n }, n represents the serial number of rotational degrees of freedom, n=1,2,...,N; N represents the total number of rotational degrees of freedom, N=2I; has rotational degrees of freedom in two directions; the vector of the rope hole located on the upper plate of vertebral segment i in the inertial coordinate system {O-XYZ} is denoted as
Figure GDA0002242193010000031
Its vector in the vertebral coordinate system {C i -x i y i z i } is
Figure GDA0002242193010000032
j represents the serial number of the rope, and the serial number of the rope hole through which the j rope passes through the vertebral segment is also j, j=1,2,...,M; M is the total number of ropes, N=3I; located on the lower plate of the vertebral segment i The vector of the rope hole on the inertial coordinate system {O-XYZ} is denoted as
Figure GDA0002242193010000033
Its vector in the vertebral coordinate system {C i -x i y i z i } is
Figure GDA0002242193010000034

从椎节i+1的下板绳孔j到椎节i的上板绳孔j的矢量为绳段矢量

Figure GDA0002242193010000035
其计算公式为:The vector from the lower plate rope hole j of vertebral segment i+1 to the upper plate rope hole j of vertebral segment i is the rope segment vector
Figure GDA0002242193010000035
Its calculation formula is:

Figure GDA0002242193010000036
Figure GDA0002242193010000036

从椎节i的下板绳孔j到椎节i-1的上板绳孔j的矢量为绳段矢量

Figure GDA0002242193010000037
其计算公式为:The vector from the lower plate rope hole j of vertebral segment i to the upper plate rope hole j of vertebral segment i-1 is the rope segment vector
Figure GDA0002242193010000037
Its calculation formula is:

Figure GDA0002242193010000038
Figure GDA0002242193010000038

则作用在椎节i的下板和上板上的绳索作用力为:Then the rope force acting on the lower plate and the upper plate of the vertebral segment i is:

Figure GDA0002242193010000039
Figure GDA0002242193010000039

Figure GDA00022421930100000310
Figure GDA00022421930100000310

其中,

Figure GDA00022421930100000311
为绳索j作用在椎节i的下板上的作用力;为绳索j作用在椎节i的上板上的作用力;
Figure GDA00022421930100000313
为绳段矢量
Figure GDA00022421930100000314
的单位矢量;
Figure GDA00022421930100000315
为绳段矢量
Figure GDA00022421930100000316
的单位矢量;fj表示绳索j的拉力大小。in,
Figure GDA00022421930100000311
is the force of the rope j acting on the lower plate of the vertebral segment i; is the force of the rope j acting on the upper plate of the vertebral segment i;
Figure GDA00022421930100000313
for rope segments vector
Figure GDA00022421930100000314
the unit vector of ;
Figure GDA00022421930100000315
for rope segments vector
Figure GDA00022421930100000316
The unit vector of ; f j represents the pulling force of the rope j.

其他步骤及参数与具体实施方式一相同。Other steps and parameters are the same as in the first embodiment.

具体实施方式三:本实施方式与具体实施方式二不同的是,绳段矢量的单位矢量

Figure GDA0002242193010000042
绳段矢量的单位矢量
Figure GDA0002242193010000044
Embodiment 3: The difference between this embodiment and Embodiment 2 is that the rope segment vector unit vector of
Figure GDA0002242193010000042
rope segment vector unit vector of
Figure GDA0002242193010000044

其他步骤及参数与具体实施方式一或二相同。Other steps and parameters are the same as in the first or second embodiment.

具体实施方式四:本实施方式与具体实施方式二不同的是,步骤二中所述绳索力的等效旋量具体为:Embodiment 4: The difference between this embodiment and Embodiment 2 is that the equivalent screw of the rope force described in step 2 is specifically:

其中,Si为椎节i上绳索产生的绳索力的等效旋量,

Figure GDA0002242193010000046
为绳索j作用在椎节i上的力在椎节i的椎节坐标系原点处产生的力旋量;f=[f1 f2 … fM]T
Figure GDA0002242193010000047
表示由绳索力向第i个椎节上力旋量映射的雅克比矩阵。Among them, S i is the equivalent screw of the rope force generated by the rope on the vertebral segment i,
Figure GDA0002242193010000046
is the force screw produced by the force of the rope j acting on the vertebral segment i at the origin of the vertebral segment coordinate system of the vertebral segment i; f=[f 1 f 2 ... f M ] T ;
Figure GDA0002242193010000047
Represents the Jacobian matrix that maps from the rope force to the force screw on the ith vertebral segment.

其他步骤及参数与具体实施方式二相同。Other steps and parameters are the same as in the second embodiment.

具体实施方式五:本实施方式与具体实施方式四不同的是,所述

Figure GDA00022421930100000411
的具体计算过程包括:Embodiment 5: This embodiment differs from Embodiment 4 in that the
Figure GDA00022421930100000411
The specific calculation process includes:

将椎节所受的绳索力分为两类:The rope forces on the vertebral segments are divided into two categories:

A、驱动椎节i运动的绳索,其编号为:j=i,I+i,2I+i;A. The ropes that drive the movement of the vertebral segment i are numbered as: j=i, I+i, 2I+i;

B、干扰椎节i+1,i+2,...,I运动的绳索编号为j=i+1,i+2,…I,I+i+1,I+i+2,…2I,2I+i+1,2I+i+2,…3I;B. The ropes that interfere with the movement of vertebral segments i+1, i+2,...,I are numbered j=i+1,i+2,...I,I+i+1,I+i+2,...2I ,2I+i+1,2I+i+2,…3I;

当绳索j为椎节i的驱动绳索时,其作用在椎节i上的力只有其作用在椎节i的下板上的作用力,即为:When the rope j is the driving rope of the vertebral segment i, the force acting on the vertebral segment i is only the force acting on the lower plate of the vertebral segment i, namely:

Figure GDA0002242193010000048
Figure GDA0002242193010000048

此时,

Figure GDA0002242193010000049
在椎节i的椎节坐标系原点处产生的力旋量为:at this time,
Figure GDA0002242193010000049
The force screw generated at the origin of the vertebral segment coordinate system of vertebral segment i is:

Figure GDA00022421930100000410
Figure GDA00022421930100000410

当绳索j为椎节i的干扰力时,其作用在椎节i上的力为其作用在椎节i的下板和下板上的作用力之的矢量和:When the rope j is the interference force of the vertebral segment i, the force acting on the vertebral segment i is the vector sum of the forces acting on the lower plate and the lower plate of the vertebral segment i:

Figure GDA0002242193010000051
Figure GDA0002242193010000051

此时,

Figure GDA0002242193010000052
在椎节i的椎节坐标系原点处产生的力旋量为:at this time,
Figure GDA0002242193010000052
The force screw generated at the origin of the vertebral segment coordinate system of vertebral segment i is:

Figure GDA0002242193010000053
Figure GDA0002242193010000053

其他步骤及参数与具体实施方式一、二或三相同。Other steps and parameters are the same as in the first, second or third embodiment.

具体实施方式六:本实施方式与具体实施方式五不同的是,步骤三中所述获得雅可比矩阵的过程为:Embodiment 6: The difference between this embodiment and Embodiment 5 is that the process of obtaining the Jacobian matrix described in step 3 is:

由式(7)、(9)得到为:Obtained from formulas (7) and (9) for:

Si转化为τ的计算方法为:The calculation method for converting Si to τ is:

Figure GDA0002242193010000056
Figure GDA0002242193010000056

其中,τi为Si引起的关节等效力矩;Ji为椎节i上的力矢量向关节等效力矩映射的雅克比矩阵,其组成为:Among them, τ i is the joint equivalent moment caused by Si; J i is the Jacobian matrix that maps the force vector on the vertebral segment i to the joint equivalent moment, and its composition is:

Figure GDA0002242193010000057
Figure GDA0002242193010000057

其中,z0,z1,…z2i-1为D-H坐标系的Z轴在惯性坐标系中的单位向量,p0,p1,…p2i为D-H坐标系的原点在惯性坐标系中的矢量坐标;Among them, z 0 , z 1 ,...z 2i-1 is the unit vector of the Z axis of the DH coordinate system in the inertial coordinate system, and p 0 , p 1 ,...p 2i is the origin of the DH coordinate system in the inertial coordinate system. vector coordinates;

全部椎节所受的力旋量转化为关节等效力矩为:The torque of all vertebral segments converted into the equivalent moment of the joint is:

Figure GDA0002242193010000058
Figure GDA0002242193010000058

结合式(14):Combining formula (14):

τ=Jtf (14)τ=J t f (14)

得到Jt为:Get Jt as:

Figure GDA0002242193010000061
Figure GDA0002242193010000061

其中,Jt为绳索力向关节等效力矩映射的雅可比矩阵。Among them, J t is the Jacobian matrix that maps the rope force to the joint equivalent moment.

其他步骤及参数与具体实施方式一、二、三、四或五相同。Other steps and parameters are the same as in the first, second, third, fourth or fifth embodiment.

具体实施方式七:本实施方式与具体实施方式二、三、四、五或六不同的是,步骤四中所述得到绳长变化速率的具体过程包括:Embodiment 7: The difference between this embodiment and Embodiments 2, 3, 4, 5 or 6 is that the specific process of obtaining the rate of change of the rope length described in step 4 includes:

将绳索用力代替,利用虚功原理能够得到:The rope is replaced by force, and the principle of virtual work can be used to obtain:

由式(4)得

Figure GDA0002242193010000063
并将其与式(3)带入式(16)得到:From formula (4), we get
Figure GDA0002242193010000063
And put it and formula (3) into formula (16) to get:

Figure GDA0002242193010000064
Figure GDA0002242193010000064

δ为变分运算符,τn表示由绳索导致的第n个转动自由度(关节)的等效力矩的大小;qn表示第n个转动自由度的角速度;δ is a variational operator, τ n represents the magnitude of the equivalent moment of the n-th rotational degree of freedom (joint) caused by the rope; q n represents the angular velocity of the n-th rotational degree of freedom;

由式(2)以及

Figure GDA0002242193010000065
得到
Figure GDA0002242193010000066
则式(17)能够改写为:By formula (2) and
Figure GDA0002242193010000065
get
Figure GDA0002242193010000066
Then formula (17) can be rewritten as:

Figure GDA0002242193010000067
Figure GDA0002242193010000067

由于

Figure GDA0002242193010000068
lj为绳索j的长度,则式(18)能够转化为:because
Figure GDA0002242193010000068
l j is the length of rope j, then equation (18) can be transformed into:

fTδl=τTδq (19)f T δl=τ T δq (19)

其中,δl=[δl1 δl2 … δlM]T;τ=[τ1 τ2 … τN]T;δq=[δq1 δq2 … δqN]TWherein, δl=[δl 1 δl 2 … δl M ] T ; τ=[τ 1 τ 2 … τ N ] T ; δq=[δq 1 δq 2 … δq N ] T ;

将式(14)带入(19)中得到:Taking equation (14) into (19), we get:

fTδl=fTJt Tδq (20)f T δl=f T J t T δq (20)

由式(20)得:From formula (20), we get:

δl=Jt Tδq (21)δl=J t T δq (21)

由于所有约束均为定常约束,因此绳长变化速率

Figure GDA0002242193010000069
与转动自由度的角速度之间的关系为:Since all constraints are steady constraints, the rate of change of rope length
Figure GDA0002242193010000069
The relationship with the angular velocity of the rotational degrees of freedom is:

Figure GDA00022421930100000610
Figure GDA00022421930100000610

其中,上标“·”表示对时间求导。in, The superscript "·" indicates derivation with respect to time.

求解逆运动学关系时(22)的关键是求取雅可比矩阵Jt,由式(14)可知Jt能够通过绳索力与等效力矩之间映射关系得到。When solving the inverse kinematics relationship (22), the key is to obtain the Jacobian matrix J t . From equation (14), it can be known that J t can be obtained through the mapping relationship between the rope force and the equivalent moment.

其他步骤及参数与具体实施方式一至六相同。Other steps and parameters are the same as those in the first to sixth embodiments.

本发明还可有其它多种实施例,在不背离本发明精神及其实质的情况下,本领域技术人员当可根据本发明作出各种相应的改变和变形,但这些相应的改变和变形都应属于本发明所附的权利要求的保护范围。The present invention can also have other various embodiments. Without departing from the spirit and essence of the present invention, those skilled in the art can make various corresponding changes and deformations according to the present invention, but these corresponding changes and deformations are all It should belong to the protection scope of the appended claims of the present invention.

Claims (1)

1. A method for determining the length change of a rope-driven super-redundancy freedom robot is characterized by comprising the following steps:
firstly, establishing a coordinate system on a rope drive serial-parallel mechanism of a rope drive super-redundancy freedom degree robot; defining a rope segment vector and obtaining the force of the rope force acting on the rope guide plate by the rope segment vector, wherein the specific process comprises the following steps:
defining any universal joint and the part between the universal joint and the next universal joint as a vertebral joint; in any vertebral segment, the rope guide plate close to the universal joint is a lower plate, and the rope guide plate far away from the universal joint is an upper plate;
with the center C of the lower plate of the vertebral segmentiEstablishing a vertebral level coordinate system for the origin { Ci-xiyiziI represents the serial number of the vertebral segment, I is 1, 2. I represents the total number of vertebral segments; is defined by the center C of the lower plateiThe direction pointing to the 1 st rope hole on the lower plate is xiThe axial direction, the direction perpendicular to the lower plate, is ziAxial direction, yiThe axis being perpendicular to xiAxis and ziA shaft; establishing a D-H coordinate system { O } at the center of the gimbaln-xnynznN represents a serial number of a rotational degree of freedom, and N is 1, 2. N represents the total number of rotational degrees of freedom, N ═ 2I; the vector of the rope hole on the upper plate of the vertebral segment i in the inertial coordinate system { O-XYZ } is recorded as
Figure FDA0002242192000000011
It is in the vertebral level coordinate system { Ci-xiyiziThe vector in (f) isj represents the serial number of the rope, and the serial number of a rope hole for the j rope to pass through on the vertebral segment is also j, j is 1, 2. M is the total number of ropes, and M is 3I; the vector of the rope hole on the lower plate of the vertebral segment i in the inertial coordinate system (O-XYZ) is recorded as
Figure FDA0002242192000000013
It is in the vertebral level coordinate system { Ci-xiyiziThe vector in (f) is
Figure FDA0002242192000000014
The vector from the lower plate rope hole j of the vertebral segment i +1 to the upper plate rope hole j of the vertebral segment i is a rope segment vector
Figure FDA0002242192000000015
The calculation formula is as follows:
Figure FDA0002242192000000016
the vector from the lower plate rope hole j of the vertebral segment i to the upper plate rope hole j of the vertebral segment i-1 is a rope segment vector
Figure FDA0002242192000000017
The calculation formula is as follows:
Figure FDA0002242192000000018
the rope forces acting on the lower and upper plates of the vertebral segment i are then:
Figure FDA0002242192000000019
Figure FDA00022421920000000110
wherein,
Figure FDA00022421920000000111
the force of cord j on the lower plate of vertebra segment i;
Figure FDA00022421920000000112
the force of the cord j on the upper plate of the vertebral segment i;
Figure FDA00022421920000000113
is a vector of rope segment
Figure FDA00022421920000000114
A unit vector of (a);
Figure FDA00022421920000000115
is a vector of rope segment
Figure FDA00022421920000000116
Unit vector of, rope segment vector
Figure FDA00022421920000000117
Unit vector ofRope segment vectorUnit vector of
Figure FDA00022421920000000120
fjRepresenting the magnitude of the tension of the rope j;
step two, utilize the rope force to act on the equivalent momentum of the rope force of the rope deflector, the concrete process includes:
wherein S isiThe equivalent amount of rotation of the rope force generated by the rope on the vertebral segment i,
Figure FDA0002242192000000022
a force momentum generated at the origin of the vertebral level coordinate system of the vertebral level i for the force of the cord j acting on the vertebral level i; f ═ f1f2… fM]TA Jacobian matrix representing a mapping of the force momentum from the cable force to the ith vertebral level;
the above-mentioned
Figure FDA0002242192000000024
The specific calculation process comprises the following steps:
the cable force on the vertebral level is divided into two categories:
A. rope for driving vertebral segment i to move is numbered as follows: j ═ I, I + I,2I + I;
B. the number of the ropes disturbing the movement of the vertebral segments I +1, I +2, … I, I +1, I +2, … 2I,2I +1,2I +2, … 3I;
when cord j is the drive cord for vertebra i, the forces acting on vertebra i are:
at this time, the process of the present invention,
Figure FDA0002242192000000026
the amount of torque generated at the origin of the vertebral level coordinate system for vertebra i is:
Figure FDA0002242192000000027
when the cord j is the distracting force of the vertebra segment i, the forces acting on the vertebra segment i are:
Figure FDA0002242192000000028
at this time, the process of the present invention,
Figure FDA0002242192000000029
the amount of torque generated at the origin of the vertebral level coordinate system for vertebra i is:
Figure FDA00022421920000000210
step three, obtaining a Jacobian matrix through a mapping relation between the rope force and the joint equivalent moment, wherein the specific process comprises the following steps:
obtained from the formulae (7) and (9)Comprises the following steps:
Figure FDA0002242192000000031
Sithe calculation method for the conversion to tau is as follows:
Figure FDA0002242192000000032
wherein, tauiIs SiInduced joint equivalent moment; j. the design is a squareiThe Jacobian matrix which is mapped from the force vector on the vertebral level i to the joint equivalent moment comprises the following components:
Figure FDA0002242192000000033
wherein z is0,z1,…z2i-1Z in D-H coordinate systemUnit vector of axis in inertial frame, p0,p1,...p2iVector coordinates of an origin of the D-H coordinate system in an inertial coordinate system;
the force rotation quantity borne by all vertebral segments is converted into joint equivalent torque as follows:
Figure FDA0002242192000000034
combined formula (14):
τ=Jtf (14)
to obtain JtComprises the following steps:
wherein, JtA Jacobian matrix for mapping the rope force to the joint equivalent moment;
step four, combining the angular speed of the rotational freedom degree and the Jacobian matrix to obtain the rope length change rate, and the concrete process comprises the following steps:
the rope force is replaced, and the principle of virtual work is utilized to obtain the following result:
Figure FDA0002242192000000036
delta is a variation operator, taunThe magnitude of the equivalent moment representing the nth rotational degree of freedom caused by the rope; q. q.snAn angular velocity representing the nth rotational degree of freedom;
is represented by the formula (2) andto obtain
Figure FDA0002242192000000042
Equation (17) can be rewritten as:
Figure FDA0002242192000000043
due to the fact that
Figure FDA0002242192000000044
ljFor the length of rope j, equation (18) can be converted into:
fTδl=τTδq (19)
wherein δ l ═ δ l [ δ l ═ δ l1δl2… δlM]T;τ=[τ1τ2… τN]T;δq=[δq1δq2… δqN]T(ii) a Substituting formula (14) into (19) yields:
δl=Jt Tδq (21)
since all constraints are constant constraints, the rate of change of the rope length
Figure FDA0002242192000000045
The relationship with the angular velocity of the rotational degree of freedom is:
Figure FDA0002242192000000046
wherein,
Figure FDA0002242192000000047
superscript "·" denotes derivation over time;
and step five, integrating the rope length change rate obtained in the step four to obtain the change of the rope length.
CN201910041486.6A 2019-01-16 2019-01-16 Rope length change determination method for rope-driven super-redundancy degree of freedom robot Active CN109760051B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910041486.6A CN109760051B (en) 2019-01-16 2019-01-16 Rope length change determination method for rope-driven super-redundancy degree of freedom robot

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910041486.6A CN109760051B (en) 2019-01-16 2019-01-16 Rope length change determination method for rope-driven super-redundancy degree of freedom robot

Publications (2)

Publication Number Publication Date
CN109760051A CN109760051A (en) 2019-05-17
CN109760051B true CN109760051B (en) 2020-02-07

Family

ID=66452493

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910041486.6A Active CN109760051B (en) 2019-01-16 2019-01-16 Rope length change determination method for rope-driven super-redundancy degree of freedom robot

Country Status (1)

Country Link
CN (1) CN109760051B (en)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111906762A (en) * 2020-06-10 2020-11-10 哈尔滨工业大学 A method for determining the joint angle of a snake-shaped manipulator
CN112249367B (en) * 2020-10-13 2022-04-05 哈尔滨工业大学 Minor planet probe maneuvering inspection device
CN112318492B (en) * 2020-10-13 2022-04-08 哈尔滨工业大学 A rope-driven snake-shaped manipulator and its control method in the state of rope failure
CN112959310B (en) * 2021-02-04 2022-04-19 清华大学深圳国际研究生院 Method for evaluating operating performance of rope-driven flexible mechanical arm
CN112975925B (en) * 2021-02-08 2021-10-26 西安电子科技大学 Rope-driven snakelike mechanical arm motion data processing method containing rope hole gaps
CN114714342B (en) * 2022-04-24 2024-01-30 哈尔滨工业大学(深圳) Rope-driven flexible arm driving rope hysteresis deformation measuring device and compensation control method thereof
CN115416046A (en) * 2022-08-31 2022-12-02 北京精密机电控制设备研究所 Joint module for super-redundant mechanical arm
CN115890673B (en) * 2022-11-23 2024-09-06 哈尔滨工业大学(深圳) A method, device and storage medium for sensing the whole arm force of a rope-driven flexible mechanical arm
CN117697765B (en) * 2024-02-05 2024-04-16 泓浒(苏州)半导体科技有限公司 Method and system for improving conveying accuracy of wafer mechanical arm based on sensing feedback

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107479564A (en) * 2017-07-13 2017-12-15 西北工业大学 The method that super redundant space robot carries out mission planning using kinematic solution
CN107545126B (en) * 2017-09-28 2019-11-26 大连理工大学 A Dynamic Response Analysis Method of Aggregated Tensioned Integral Structure Based on Sliding Cable Elements of Multibody System
CN108555914B (en) * 2018-07-09 2021-07-09 南京邮电大学 A DNN neural network adaptive control method based on tendon-driven dexterous hand
CN108908332B (en) * 2018-07-13 2021-07-30 哈尔滨工业大学(深圳) Control method and system of super-redundant flexible robot, and computer storage medium
CN109176494B (en) * 2018-09-28 2022-03-29 哈尔滨工业大学(深圳) Self-calibration method and system for rope-driven multi-joint flexible robot and storage medium

Also Published As

Publication number Publication date
CN109760051A (en) 2019-05-17

Similar Documents

Publication Publication Date Title
CN109760051B (en) Rope length change determination method for rope-driven super-redundancy degree of freedom robot
CN110561425B (en) Rope-driven flexible robot force and position hybrid control method and system
Vähä et al. Dynamic model of excavator
CN109176494B (en) Self-calibration method and system for rope-driven multi-joint flexible robot and storage medium
CN108466289B (en) Parallel robot dynamics modeling method considering joint friction
CN113305809B (en) Force sense simulation control method of fully-constrained space rope driving parallel mechanism
CN105196284B (en) Three-degree-of-freedom tandem type self-gravity-balance passive mechanical arm
CN109746913B (en) Method and system for teaching robot posture keeping dragging
CN103558079A (en) Multi-degree-of-freedom loading method based on parallel mechanism driving force closed loop
CN104477420A (en) Nine-rope driving robot device for simulating zero gravity and low gravity
CN113465861B (en) Interference force compensation method of two-degree-of-freedom electro-hydraulic vibration table
CN107270061B (en) A control method for a six-degree-of-freedom hydraulic motion platform with connecting rods
CN108515518B (en) A workspace solution method for flexible support industrial robot
CN113465859A (en) Interference force compensation method of six-degree-of-freedom electro-hydraulic vibration table
CN112936287B (en) A flexible robot control method and device based on dynamic iterative learning
CN101844349A (en) Redundant parallel robot system
CN107121930B (en) Six Degree-of-Freedom Parallel Platform motion control method based on dummy model
CN109657282A (en) A kind of H-type motion platform modeling method based on lagrangian dynamics
CN116141319A (en) Control method, device, equipment and storage medium of three-rope parallel robot
KR102484428B1 (en) Control system for flexible joint robots
CN112231888B (en) Ideal dynamics model establishing method based on laser tracking measurement system mechanical structure
CN114211503B (en) Rope-driven flexible robot track control method and system based on visual feedback
CN114859698B (en) A trajectory tracking control method and device for a 6-3-PUS parallel mechanism
CN110298129B (en) Modeling method for telescopic vibration characteristic of straight-arm type aerial work platform arm support
CN116126003A (en) Wave compensation system modeling and pose control method based on Stewart platform

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant