Background
The reference is made to a Chinese patent with an authorization announcement number of CN109091237B and a name of an auxiliary system of a minimally invasive surgical instrument, a Chinese patent with an authorization announcement number of CN109091238B and a name of an auxiliary system of a split type minimally invasive surgical instrument, a utility model patent with an authorization announcement number of CN210872029U and a name of a doctor operating table, and the minimally invasive surgical robot is designed to be an equipment system which is operated in an operating room by a doctor and controls the surgical instrument to perform an operation under an endoscope. The doctor sits in front of the doctor operation panel, and through watching the 3D image display, operates the doctor arm, and the doctor can utilize the accurate control surgical instruments of patient arm to carry out various operation actions like ordinary operation.
Referring to the patent application publication nos. CN105286999A, CN105286989A, surgical instruments can perform different functions, including clamping, resection, cutting, suturing, stapling, etc. Surgical instruments come in different configurations, including an execution tip, wrist, instrument bar, instrument box, and the like.
In the existing minimally invasive surgery robot system, a doctor outputs an action command by holding a mechanical motion input device with multiple degrees of freedom, the collected action command can be converted into joint motion information of a mechanical arm through a master-slave motion mapping algorithm, the mechanical arm follows the motion of a master mechanical arm, and the current algorithm is poor in real-time performance, accuracy and precision.
Disclosure of Invention
The invention provides a master-slave tracking control method for a minimally invasive surgery robot, which aims to solve the technical problems of poor instantaneity, poor accuracy and low precision of the master-slave tracking control method of the existing minimally invasive surgery robot system and improve the poor instantaneity, poor accuracy and precision.
The invention provides a master-slave tracking control method for a minimally invasive surgery robot, which comprises the following steps:
establishing a D-H parameter table of the main operating arm by using a D-H method, wherein the D-H parameter table comprises the connecting rod lengths a of three joints of M1, M2 and M3iOffset d of the connecting rodiAngle of rotation of connecting rod alphaMiAnd joint angle thetaMiAccording to the principle of space transformation, a homogeneous transformation matrix between two adjacent joints is obtained according to the following formula:
establishing transformation between each adjacent joint according to the above formula, and then solving a homogeneous transformation matrix of the terminal coordinate system of the main operating arm relative to the base coordinate system according to the following formula
Obtaining the position of the tail end of the main operating arm relative to the base coordinate system from the homogeneous transformation matrix;
will be provided with
The following is written:
then, the following steps are carried out:
PMthe position coordinates of the wrist of the main operation arm are obtained;
establishing a corresponding dynamic coordinate system at the tail end position point of the main operating arm, and calculating the rotation angle values theta of the three joints M4, M5 and M6 of the wrist according to the following formulaMiIndicating the posture of thewrist 4;
θMi=AMi/μMi
wherein A isMiRepresents the amount of change, mu, of the joint encoder valueMiThe product of the joint encoder precision, the speed reduction ratio of the speed reducer and the transmission ratio of the mechanical structure;
establishing a kinematic model by a D-H parameter method aiming at the end of the operation arm from a joint S1, a joint S2 and a joint L3 of the operation arm, a joint S4, a joint S5 and a joint S6 of the wrist at the tail end of the surgical instrument, and obtaining a homogeneous transformation matrix of the tail end of the operation arm in a base coordinate system according to a homogeneous transformation matrix calculation formula
The rotation angle values theta of the three joints of the wrist of the surgical instrument are calculated according to the following formula
Si:
θSi=ASi/μSi
Mapping the position base coordinate system of the master operating arm into the base coordinate system described by the slave operating arm, and obtaining a homogeneous transformation matrix of the slave operating arm after mapping as follows:
wherein,
is a transformation matrix from the base coordinate system of the master operation arm to the base coordinate system of the slave operation arm,
a scale factor that is a master-slave mapping;
from the position coordinates of the operating arm
According to a homogeneous transformation matrix of the joints S1, S2 and L3 of the operation arm, which is derived from the configuration of the mechanical arm according to a D-H parameter method:
wherein, the following relation is further established:
psx=xs
psy=ys
psz=zs
according to the three equations, the unknown parameters theta of the joints S1, S2 and L3 can be obtaineds1、θs2、ls3;
The calculation formula of the expected angle values of the wrist three joints at the tail end of the operation arm instrument is as follows:
θs4=θM4
θs5=θM5
θs6=θM6
wherein, thetaM4、θM5、θM6Respectively representing the actual rotation angle values of three joints of the main operation wrist part;
and then the motion controller is used for outputting a corresponding numerical value to the slave operation arm to realize master-slave following.
The beneficial effects of the invention are: the invention discloses a simplified kinematic modeling method combining a heterogeneous mapping mode and a homogeneous mapping mode. The position part of the master-slave manipulator adopts a heterogeneous mapping mode, and the posture part of the master-slave manipulator adopts an isomorphic mapping mode. The complexity of the whole kinematics calculation process is simplified, and the real-time performance, accuracy and precision of the master-slave tracking motion control system are improved.
Further features and aspects of the present invention will become apparent from the following description of the embodiments with reference to the accompanying drawings.
Drawings
FIG. 1 is a schematic structural diagram of a main operating arm of a minimally invasive surgical robot;
FIG. 2 is a schematic diagram of the structure of the wrist shown in FIG. 1;
FIG. 3 is a block diagram of the wrist shown in FIG. 2;
FIG. 4 is a configuration diagram of the main operating arm position adjustment;
FIG. 5 is a schematic view of the slave operating arm in a folded state;
FIG. 6 is a schematic view of the slave arm in a deployed state;
FIG. 7 is a schematic view of the slave arm in a deployed state;
FIG. 8 is a configuration view of the slave arm;
FIG. 9 is a configuration diagram of a surgical instrument;
FIG. 10 is a mechanical schematic of the surgical instrument;
FIG. 11 is a configuration diagram of the wrist.
The symbols in the drawings illustrate that:
1. the device comprises a base, 2, a first arm rod, 3, a second arm rod, 4, a wrist, 4-1, a fixed connecting rod, 4-2, a first L-shaped connecting rod, 4-3, a second L-shaped connecting rod, 4-4, a handle, 4-5, an opening and closing seat, 4-6, a first encoder, 4-7, a second encoder and 4-8, a third encoder, wherein the first arm rod is arranged on the base; 5. the connecting seat is rotated. 201. The base end base, 202, the connecting rod base, 203, the slave end connecting rod I, 204, the slave end connecting rod II, 205, the instrument installation device, 206, the instrument lifting base, 207, the surgical instrument, 208, the motor, 209 and the motor.
Detailed Description
The main operating arm applied to the invention is a doctor mechanical arm in the utility model patent with the name of doctor operating table and the publication number of CN 210872029U. The slave manipulator arm used is the mechanical arm of patent No. 201711314221.6, named the master manipulator arm for minimally invasive surgery.
As shown in fig. 1, the rear end of thefirst arm lever 2 is connected to thebase 1 through a first joint, thefirst arm lever 2 can rotate on the horizontal plane, therotation connecting base 5 is connected to the front end of thefirst arm lever 2 through a second joint, therotation connecting base 5 can rotate on the horizontal plane, thesecond arm lever 3 is connected to therotation connecting base 5 through a third joint, and thesecond arm lever 3 can rotate on the vertical plane.
In the operation process, a doctor holds thewrist 4 with a hand to operate, and the joints of the main operation arm rotate around afirst joint axis 14 of the first joint, asecond joint axis 15 of the second joint and athird joint axis 16 of the third joint respectively, wherein thethird joint axis 16 is vertical to the gravity direction, and the firstjoint axis 14 and the secondjoint axis 15 are vertical to the ground.
As shown in figures 2 and 3, thewrist 4 comprises a fixed connecting rod 4-1, a first L-shaped connecting rod 4-2, a second L-shaped connecting rod 4-3, a handle 4-4, an opening and closing seat 4-5, a first encoder 4-6, a second encoder 4-7 and a third encoder 4-8. Thewrist 4 is provided with a joint M4, a joint M5, and a joint M6.
Fig. 4 shows a main operation arm position adjustment configuration diagram, wherein the three joints M1, M2 and M3 are the first joint, the second joint and the third joint in the configuration diagram shown in fig. 1. As shown in FIGS. 3 and 11, the three joints M4, M5 and M6 are posture-adjustingAnd the 6 joint degrees of freedom respectively determine the position and the posture of the tail end of the main operation arm in a Cartesian coordinate system. According to the configuration of the main operating arm, the kinematic analysis of the main operating arm can be carried out by adopting a D-H method, and a D-H parameter table is established, wherein the D-H parameter table comprises the connecting rod lengths a of three joints of M1, M2 and M3iOffset d of the connecting rodiAngle of rotation of connecting rod alphaMiAnd joint angle thetaMi. For example, in FIG. 1, a1Denotes the distance, a, between the axis of the joint M1 and the axis of the joint M22Denotes the distance, a, between the axis of the joint M2 and the axis of the joint M33Indicating the distance between the axis of the joint M3 and thewrist 4. The link offset is the distance between two adjacent common perpendicular lines along the axial direction, and the link offset diAre all 0. Knowing the parameters of each joint, a homogeneous transformation matrix between two adjacent joints is obtained according to the principle of spatial transformation as follows.
The transformation between each adjacent joint is established according to the above formula, and then the controller can obtain a homogeneous transformation matrix of the coordinate system of the tail end of the main operation arm (namely the coordinate system of the wrist 4) relative to the base coordinate system according to the following formula
And i takes values of 1, 2 and 3, and the position of the tail end of the main operating arm relative to the base coordinate system can be obtained from the homogeneous transformation matrix.
The following can also be written:
then, the following steps are carried out:
PMi.e. the position coordinates of thewrist 4 of the main operating arm.
Establishing a corresponding dynamic coordinate system at the tail end position point of the main operating arm, and calculating the rotation angle values theta of the three joints M4, M5 and M6 of thewrist 4 according to the following formulaMiAnd represents the posture of thewrist 4.
θMi=AMi/μMi
Wherein A isMiRepresents the amount of change, mu, of the joint encoder valueMiThe product of the joint encoder precision, the speed reduction ratio of the speed reducer and the transmission ratio of the mechanical structure. Referring to FIG. 3, θM4For the rotation angle value of the joint M4, the corresponding encoder is the first encoder 4-6, thetaM5For the value of the rotation angle of the joint M5, the corresponding encoder is the second encoder 4-7, thetaM6For the rotational angle value of the joint M6, the corresponding encoder is the third encoder 4-8.
As shown in fig. 5-7, the slave manipulator arm comprises aslave end base 201, alink base 202, a slave end link I203, a slave end link II204, aninstrument lift base 206, and aninstrument mounting device 205. Thelink base 202 is rotatably mounted at one end of theslave base 201 with its axis of rotation aligned with the direction of gravity. Thesurgical instrument 207 may be mounted on theinstrument mounting device 205.
As shown in fig. 8, the configuration of the slave arm is such that the joint S1 is a rotary joint for moving the
link base 202, the slave end link I203, the slave end link II204, and the
instrument lifter base 206 in the entire lateral direction, the joint S2 is a parallelogram-shaped pivoting joint that pivots in the pitch direction about the axis (i.e., extends or folds the slave end link I203, the slave end link II204, and the instrument lifter base 206), the joint L3 is a joint for moving the
instrument attachment device 205 along the
instrument lifter base 206, the joints S1 and S2 define the position of the distal dead point (point P in fig. 7), and the joint L3 that moves in the vertical direction about the axis of the distal dead point defines the position of point P. Such asAs shown in fig. 9 and 10, the surgical instrument adopts a wrist degree of freedom layout of "self-transmission-deflection-end self-transmission", and the joints S4, S5 and S6 of the wrist at the end of the surgical instrument determine the posture of the end of the surgical instrument. Similarly, the slave manipulator end also adopts a D-H parameter method to establish a kinematic model of the slave manipulator, and the homogeneous transformation matrix of the tail end of the slave manipulator in the base coordinate system is obtained according to the homogeneous transformation matrix calculation formula

The position of the surgical instrument tip in the base coordinate system from the manipulator arm can be derived from the homogeneous transformation matrix. The rotation angle values theta of the three joints of the wrist of the surgical instrument can be calculated according to the following formula
SiAre each theta
S4、θ
S5、θ
S6And represents the posture of the wrist from the end of the manipulator instrument.
θSi=ASi/μSi
Wherein A isSiRepresents the variation of the joint encoder value, muSiIs the product of the precision of the joint encoder, the deceleration ratio of the speed reducer and the transmission ratio of the mechanical structure.
According to the designed configurations of the main operating arm and the auxiliary operating arm, the surgical instrument and the wrist three joint axes of the isomorphic main hand at the tail end of the surgical instrument are intersected at one point, and the decoupling of the position and the posture in the kinematic analysis can be realized. By combining the structural characteristics, the invention establishes master-slave mapping with separated pose and isomorphic tail ends. In the position and motion mapping, a position base coordinate system of the master operation arm needs to be mapped into a base coordinate system described by the slave operation arm, and a homogeneous transformation matrix of the slave operation arm position is obtained after mapping as follows:
wherein,
is a transformation matrix from the base coordinate system of the master operation arm to the base coordinate system of the slave operation arm,
the motion of the main manipulator is mapped into the motion of the manipulator after being reduced according to a certain proportion for the scale factor of the master-slave mapping.
The following can also be written:
from the position coordinates of the operating arm
A homogeneous transformation matrix of the joints S1, S2 and L3 of the operation arm is derived according to a D-H parameter method from the configuration of the mechanical arm:
wherein each element in the matrix is a joint 1 to joint 3 parameter (θ)s1、θs2、ls3) Is described in (1). The following relationships are established according to the above description:
psx=xs
psy=ys
psz=zs
according to the three equations, the unknown parameters theta of the joints S1, S2 and L3 can be obtaineds1、θs2、ls3。
And the posture mapping of the instrument and the tail end of the main hand is realized by adopting a one-to-one mapping mode. Therefore, the calculation formula of the expected angle values of the wrist three joints at the tail end of the operation arm instrument is as follows:
θs4=θM4
θs5=θM5
θs6=θM6
wherein, thetaM4、θM5、θM6Respectively representing the actual rotation angle values of the three joints of the main operating wrist. Then the motion controller is used for outputting corresponding numerical values to the slave operation arm, and the purpose of master-slave following is achieved. The invention provides a modeling method for simplifying a kinematics model of a minimally invasive surgery robot according to the concept of master-slave operation arm configuration, namely the concept of separation of postures and isomorphism of tail ends, simplifies the complexity of the whole kinematics calculation process and improves the real-time property of a motion control system.
The above description is only for the purpose of illustrating the preferred embodiments of the present invention and is not to be construed as limiting the present invention, since various modifications and changes will be apparent to those skilled in the art.