CN108972550B - Inverse kinematics solving method of concentric tube robot - Google Patents

Inverse kinematics solving method of concentric tube robot Download PDF

Info

Publication number
CN108972550B
CN108972550B CN201810752589.9A CN201810752589A CN108972550B CN 108972550 B CN108972550 B CN 108972550B CN 201810752589 A CN201810752589 A CN 201810752589A CN 108972550 B CN108972550 B CN 108972550B
Authority
CN
China
Prior art keywords
tube
robot
length
inner tube
concentric
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
CN201810752589.9A
Other languages
Chinese (zh)
Other versions
CN108972550A (en
Inventor
宋霜
王杰
张迪
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shenzhen Graduate School Harbin Institute of Technology
Original Assignee
Shenzhen Graduate School Harbin Institute of Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shenzhen Graduate School Harbin Institute of Technology filed Critical Shenzhen Graduate School Harbin Institute of Technology
Priority to CN201810752589.9A priority Critical patent/CN108972550B/en
Publication of CN108972550A publication Critical patent/CN108972550A/en
Application granted granted Critical
Publication of CN108972550B publication Critical patent/CN108972550B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • BPERFORMING OPERATIONS; TRANSPORTING
    • B25HAND TOOLS; PORTABLE POWER-DRIVEN TOOLS; MANIPULATORS
    • B25JMANIPULATORS; CHAMBERS PROVIDED WITH MANIPULATION DEVICES
    • B25J9/00Programme-controlled manipulators
    • B25J9/16Programme controls
    • B25J9/1602Programme controls characterised by the control system, structure, architecture
    • B25J9/1607Calculation of inertia, jacobian matrixes and inverses

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Automation & Control Theory (AREA)
  • Robotics (AREA)
  • Mechanical Engineering (AREA)
  • Manipulator (AREA)
  • Numerical Control (AREA)

Abstract

The invention provides a concentric tube robot inverse kinematics solving method, which comprises the following steps: s1, generating a DDA table; s2, use of DDA table; step S2 includes the following substeps; s21, determining according to the coordinate information of the target pointDetermining a corresponding column in a DDA table; s22, determining the length l of the outer tube1And length l of the inner tube2(ii) a S23, determining the rotation angle phi of the inner tube relative to the world coordinate system2(ii) a S24, applying the inverse kinematics solution result of the concentric tube robot to the motion control of the concentric tube robot consisting of a straight tube and a bent tube, wherein the straight tube is an outer tube, and the bent tube is an inner tube. The invention has the beneficial effects that: the precision of the motion control of the concentric tube robot can be improved.

Description

Inverse kinematics solving method of concentric tube robot
Technical Field
The invention relates to a concentric tube robot, in particular to an inverse kinematics solving method of the concentric tube robot.
Background
In recent years, a concentric tube robot has played an increasingly important role as a novel continuous robot in the medical field. Due to the characteristics of small size and flexible bending, the concentric tube robot can better show the advantages in the complex and multi-obstacle narrow space. However, it is difficult to solve the inverse kinematics of a concentric tube robot relative to a robot having a definite rigid joint. So far, when the concentric tube robot is composed of variable-curvature concentric tubes or the number of the concentric tubes exceeds 2, it is extremely difficult or even impossible to solve the inverse kinematics.
There are currently some generally accepted approaches to solving the forward and reverse kinematics of concentric tube robots. First, assuming that the concentric tubes are piecewise constant curvature, there is only a bending force on the tube-to-tube interaction. Based on the assumption, a kinematic model of the concentric tube robot is established, and positive kinematics is solved by utilizing the Euler-Bernoulli beam principle. The forward kinematics was analyzed step by step in the reverse direction to solve the inverse kinematics of the concentric tube robot. Another is the exponential-product of-exponential method, which uses the related knowledge of lie group algebra and rotation theory to solve the forward and inverse kinematics.
Although there are currently some feasible methods in both modeling of the kinematics of a concentric tube robot and solution of the kinematic equations, there are some problems with these methods. Both of these problems reduce the accuracy of the concentric tube robot motion control. The related problems are now expressed as follows:
1. when the kinematic model is built, the concentric tubes are idealized to be of a piecewise constant curvature, which is difficult to realize in real-world tube processing.
2. The most common concentric tube material today is a superelastic nitinol tube, but over time the tube undergoes irreversible plastic deformation, indirectly resulting in less and less accurate kinematic solutions.
3. Almost all models do not consider the problems of tube-to-tube shearing and axial deformation.
Disclosure of Invention
In order to solve the problems in the prior art, the invention provides an inverse kinematics solving method of a concentric tube robot, which can improve the motion control precision of the concentric tube robot.
The invention provides a concentric tube robot inverse kinematics solving method, which comprises the following steps:
s1, generating a DDA table;
s2, use of DDA table;
step S2 includes the following substeps;
s21, determining a corresponding column in the DDA table according to the coordinate information of the target point;
s22, determining the length l of the outer tube1And length l of the inner tube2
S23, determining the rotation angle phi of the inner tube relative to the world coordinate system2
S24, applying the inverse kinematics solution result of the concentric tube robot to the motion control of the concentric tube robot consisting of a straight tube and a bent tube, wherein the straight tube is an outer tube, and the bent tube is an inner tube.
As a further improvement of the present invention, step S1 includes: first, the outer tube length of the concentric tube robot remains at a maximum, i.e./1=l1maxAngle of rotation phi of the inner tube coordinate system relative to the world coordinate system2Remains at 0; the inner tube is then moved axially relative to the outer tube to a length l of the inner tube2From l1maxMovement changes to2maxObtaining a DDA table with 6 rows and n columns, wherein the first three rows of the DDA table are respectively the length l of the outer tube1Inner tubeLength l of2And angle of rotation phi of the inner tube2The last three rows of the table are the components of the concentric tube robot tip position in the x, y, z axes for the same column input.
As a further development of the invention, the inner tube length l2From l1maxChange to2maxThis course of movement represents the course of movement of the inner tube relative to the outer tube; in the initial state, the inner pipe and the outer pipe are completely overlapped, namely the inner pipe is completely nested in the outer pipe; in the movement process, the inner pipe moves axially relative to the outer pipe, namely the inner pipe gradually extends out of the outer pipe; assuming that the arc length of the inner tube extending out of the outer tube is denoted by s, the length l of the inner tube2From l1maxMovement changes to2maxIs the length l of the inner tube2Gradually increasing in certain increment, the inner pipe gradually extends out of the outer pipe from s-0 to s-l21In which l21Showing the arc length of the 1 st bend of the inner tube distal from the mounting end.
As a further improvement of the present invention, in step S2, the coordinates (x) of any point in the feasible space of the concentric tube robot are knowni,yi,zi) The data obtained from the DDA table are the length l of the outer tube1Length of the inner tube l2And the rotation angle phi of the inner tube relative to the world coordinate system2
As a further improvement of the present invention, step S21 includes: when the inner tube of the concentric tube robot rotates for one circle relative to a world coordinate system, a locus circle is drawn in space by a point at any arc length s on the inner tube, and the size of the locus circle is determined by the extension arc length of the inner tube relative to the tail end of the outer tube; therefore, the radius r of the track circle is used to select the corresponding column in the DDA table, and the rotating radius r corresponding to the target point is matched with the rotating radius r corresponding to the end point in the DDA tableiAnd comparing one by one, and selecting the column with the minimum difference value to prepare for solving the input quantity of the following robot.
As a further improvement of the present invention, step S22 includes: in step S21, one of the columns in the DDA table has been locked, that is, the specific value (l) of the input amount when the robot reaches the target position1i,l2i) (ii) a In the generation of DDA tablesDDA Table is generated on the premise that the outer tube length reaches a maximum value l1=l1max(ii) a Therefore, when the target point does not satisfy this condition,/1iAnd l2iThe value of (c) needs to be corrected; assuming that the inner tube rotates to obtain a first trajectory circle when the length of the outer tube is maximum and the end point of the concentric tube robot rotates to obtain a second trajectory circle when the target position is actually reached, the correct value is represented as (l)1,l2),
l1=l1i+e
l2=l2i+e
Wherein e>0 is a correction value, i.e. the center O of the first locus circle1And a center O of a second locus circle2The distance of (c).
As a further improvement of the present invention, step S23 includes: assuming T is the target location, all that needs to be determined is the vector O1T and vector O1The angle of rotation between G; suppose the vector O1T and vector O1The angle of rotation between G is theta2iThen there is
Figure GDA0002946633490000041
θ2iThe angle being in the range of [0,180 °];
Assuming that the target point is rotated within 0,180 deg. counterclockwise, the angle of rotation is positive, and vice versa, there is,
Figure GDA0002946633490000042
alpha represents the rotation angle theta capable of distinguishing the two side spaces2iThe coordinate value of the sign is determined.
The invention has the beneficial effects that: through the scheme, the motion control precision of the concentric tube robot can be improved.
Drawings
FIG. 1 is a flow chart of an inverse kinematics solution method for a concentric tube robot according to the present invention.
FIG. 2 is a schematic diagram of the rotation of the inner tube of the inverse kinematics solution method of the concentric tube robot according to the present invention.
FIG. 3 is a schematic diagram of length correction of an inverse kinematics solution method of a concentric tube robot according to the present invention.
Detailed Description
The invention is further described with reference to the following description and embodiments in conjunction with the accompanying drawings.
As shown in fig. 1 to 3, an inverse kinematics solution method for a concentric tube robot includes the following steps:
s1, generating a DDA table, wherein DDA is an abbreviation of English Date-drive Approach and indicates that the method is a method based on data acquisition and processing;
s2, use of DDA table;
step S2 includes the following substeps;
s21, determining a corresponding column in the DDA table according to the coordinate information of the target point;
s22, determining the length l of the outer tube 11And the length l of the inner tube 22
S23, determining the rotation angle phi of the inner tube 2 relative to the world coordinate system2
S24, applying the inverse kinematics solution result of the concentric tube robot to the motion control of the concentric tube robot consisting of a straight tube and a bent tube, wherein the straight tube is the outer tube 1, and the bent tube is the inner tube 2.
The invention provides an inverse kinematics solving method of a concentric tube robot, which is a method based on data acquisition and table lookup solving, and is applied to the concentric tube robot consisting of a straight tube and a bent tube, and in order to express the use of the inverse kinematics solving method clearly, the detailed description is as follows:
1. DDA Table Generation
The generation of the DDA table is obtained based on the manner of data acquisition. As shown in FIG. 2, first the length of the outer tube 1 of the concentric tube robot is kept at a maximum, i.e./1=l1maxThe rotation angle Φ of the coordinate system of the inner tube 2 with respect to the world coordinate system is kept at 0. Then the length l of the inner tube 22From l1maxExercise of sportsTo l2maxAcquiring coordinate data of the tail end of the concentric tube robot and corresponding robot motion input quantity in the motion process, namely the length l of the inner tube 22Recording the coordinate data (x) of the tail end of the concentric tube robot in each extension incrementi,yi,zi) And the motion input l of the concentric tube robot1i,l2i,Φ2i. The size of the motion increment of the inner tube 2 finally determines the calculation precision of inverse kinematics by using the method. Finally, a DDA table of 6 rows and n columns is obtained, wherein the first three rows of the table are respectively the length l of the outer tube 11Length l of inner tube 22And the angle of rotation phi of the inner tube 22The last three rows of the table are the components of the concentric tube robot tip position in the x, y, z axes for the same column input.
Length l of inner tube 22From l1maxChange to2maxThis course of movement represents the course of movement of the inner tube 2 relative to the outer tube 1; in an initial state, the inner tube 2 and the outer tube 1 are completely overlapped, namely, the inner tube 2 is completely nested in the outer tube 1; in the movement process, the inner tube 2 moves axially relative to the outer tube 1, namely, the inner tube 2 gradually extends out of the outer tube 1; assuming that the arc length of the inner tube 2 protruding from the outer tube 1 is denoted by s, the length l of the inner tube 22From l1maxMovement changes to2maxIs the length l of the inner tube 22The inner tube 2 is gradually extended out of the outer tube 1 in increments from s to l21In which l21Showing the arc length of the curved portion of the inner pipe 2 at the 1 st segment away from the mounting end.
2. The use of DDA tables;
inverse kinematics solution is the process of finding the robot input quantity required by the robot end to reach the target position of the known robot. That is, the coordinate (x) of any point in the feasible space of the known concentric tubei,yi,zi) The data obtained from the DDA table are the length l of the outer tube 11Length l of inner tube 22And a rotation angle phi2(since the outer tube 1 is a straight tube, the angle of rotation phi here1No determination is required). The solving process can be divided into three steps, which are described separately below.
1) Determining corresponding columns in the DDA table according to the coordinate information of the target point
When the inner tube 2 of the concentric tube robot rotates one circle relative to the world coordinate system, a point at any arc length s on the inner tube 2 draws a circle in space, as shown in fig. 2. It will be readily apparent from the figures that the length of the arc of protrusion of the inner tube 2 relative to the end of the outer tube 1 determines the size of this circle. The radius r of the circle is used here to select the corresponding column in the DDA table. The rotating radius r corresponding to the target point is compared with the rotating radius r corresponding to the end point in the DDA tableiAnd comparing one by one, and selecting the column with the minimum difference value to prepare for solving the input quantity of the following robot.
2) Determining the length l of the outer and inner tubes1,l2
In the last step, one of the columns in the DDA table has been locked, that is, the specific value (l) of the input amount when the robot reaches the target position1i,l2i). But from the analysis of the feasible space of a dual-tube concentric tube robot, it can be seen that this set of values is not necessarily the final correct solution. In the generation of the front table, the DDA table is generated on the premise that the length of the outer tube 1 reaches the maximum value l1=l1max. Therefore, when the target point does not satisfy this condition,/1iAnd l2iThe value of (c) needs to be corrected. As shown in fig. 3, it is assumed that the dotted circle is obtained by rotating the inner tube 2 when the length of the outer tube 1 is maximum, and the solid circle is obtained by rotating the end point of the concentric tube robot when the target position is actually reached. Then, the correct value should be expressed as (l)1,l2),
l1=l1i+e
l2=l2i+e
Wherein e>0 is a correction number, i.e. two centers O in FIG. 31And O2The distance of (c).
3) Determining the rotation angle phi of the inner tube 2 relative to the world coordinate system2
Assuming that T is the target location, as shown in FIG. 3, all that needs to be determined is the vector O1T and vector O1G. Suppose that the rotation angle is θ2iThen there is
Figure GDA0002946633490000071
Where theta is2iThe angle being in the range of [0,180 °]It is clear that when the concentric tube robot rotates forward and backward, the results are not the same, and therefore the sign of the rotation angle needs to be determined. Assume a counterclockwise rotation of the target point of [0,180 °]Within the scope of which the angle of rotation is positive and vice versa, there are,
Figure GDA0002946633490000072
where α represents a rotation angle θ capable of distinguishing spaces on both sides2iThe coordinate value of the sign is determined. If the inner tube coordinate system F2(0) Is shown in fig. 3, then α is the y coordinate value of the target point, otherwise it is the projection value of the target point on the y axis.
The inverse kinematics problem of the double-tube concentric tube robot can be calculated by the method.
The invention provides an inverse kinematics solving method of a concentric tube robot, which is based on inverse kinematics solving of data acquisition, and the method does not need to consider the specific curvature value of a concentric tube, whether the concentric tube is a constant curvature tube or not, and the influence of the shearing, bending, friction and other acting forces between the tubes on the precision of a conventional inverse kinematics calculating method. Since these factors, which are practically difficult to control, are already included in the DDA table. At the same time, this method also solves the problem of irreversible plastic deformation of the concentric tubes over time, since the change in curvature of the tubes is also included in the DDA table. Therefore, the inverse kinematics solution method is simple to operate and high in precision, and can effectively improve the precision of the motion control of the concentric tube robot.
The foregoing is a more detailed description of the invention in connection with specific preferred embodiments and it is not intended that the invention be limited to these specific details. For those skilled in the art to which the invention pertains, several simple deductions or substitutions can be made without departing from the spirit of the invention, and all shall be considered as belonging to the protection scope of the invention.

Claims (7)

1. The inverse kinematics solving method of the concentric tube robot is characterized by comprising the following steps of:
s1, generating a DDA table;
s2, use of DDA table;
step S2 includes the following substeps;
s21, determining a corresponding column in the DDA table according to the coordinate information of the target point;
s22, determining the length l of the outer tube1And length l of the inner tube2
S23, determining the rotation angle phi of the inner tube relative to the world coordinate system2
S24, applying the inverse kinematics solution result of the concentric tube robot to the motion control of the concentric tube robot consisting of a straight tube and a bent tube, wherein the straight tube is an outer tube, and the bent tube is an inner tube.
2. The concentric tube robot inverse kinematics solution method according to claim 1, wherein step S1 comprises: first, the outer tube length of the concentric tube robot remains at a maximum, i.e./1=l1maxAngle of rotation phi of the inner tube coordinate system relative to the world coordinate system2Remains at 0; the inner tube is then moved axially relative to the outer tube to a length l of the inner tube2From l1maxMovement changes to2maxObtaining a DDA table with 6 rows and n columns, wherein the first three rows of the DDA table are respectively the length l of the outer tube1Length of the inner tube l2And angle of rotation phi of the inner tube2The last three rows of the table are the components of the concentric tube robot tip position in the x, y, z axes for the same column input.
3. The concentric tube robot inverse kinematics solution method according to claim 2, wherein: length l of inner tube2From l1maxChange to2maxThe motion process represents the inner tube relativeThe movement process of the outer pipe; in the initial state, the inner pipe and the outer pipe are completely overlapped, namely the inner pipe is completely nested in the outer pipe; in the movement process, the inner pipe moves axially relative to the outer pipe, namely the inner pipe gradually extends out of the outer pipe; assuming that the arc length of the inner tube extending out of the outer tube is denoted by s, the length l of the inner tube2From l1maxMovement changes to2maxIs the length l of the inner tube2Gradually increasing in certain increment, the inner pipe gradually extends out of the outer pipe from s-0 to s-l21In which l21Showing the arc length of the 1 st bend of the inner tube distal from the mounting end.
4. The concentric tube robot inverse kinematics solution method according to claim 3, wherein: in step S2, the coordinates (x) of any point in the feasible space of the concentric tube robot are knowni,yi,zi) The data obtained from the DDA table are the length l of the outer tube1Length of the inner tube l2And the rotation angle phi of the inner tube relative to the world coordinate system2
5. The concentric tube robot inverse kinematics solution according to claim 4, comprising, at step S21: when the inner tube of the concentric tube robot rotates for one circle relative to a world coordinate system, a locus circle is drawn in space by a point at any arc length s on the inner tube, and the size of the locus circle is determined by the extension arc length of the inner tube relative to the tail end of the outer tube; therefore, the radius r of the track circle is used to select the corresponding column in the DDA table, and the rotating radius r corresponding to the target point is matched with the rotating radius r corresponding to the end point in the DDA tableiAnd comparing one by one, and selecting the column with the minimum difference value to prepare for solving the input quantity of the following robot.
6. The concentric tube robot inverse kinematics solution method according to claim 5, wherein step S22 comprises: in step S21, one of the columns in the DDA table has been locked, that is, the specific value (l) of the input amount when the robot reaches the target position1i,l2i) (ii) a In the generation of DDA tables, the DDA tables are generatedThe premise is that the length of the outer tube reaches a maximum value l1=l1max(ii) a Therefore, when the target point does not satisfy this condition,/1iAnd l2iThe value of (c) needs to be corrected; assuming that the inner tube rotates to obtain a first trajectory circle when the length of the outer tube is maximum and the end point of the concentric tube robot rotates to obtain a second trajectory circle when the target position is actually reached, the correct value is represented as (l)1,l2),
l1=l1i+e
l2=l2i+e
Wherein e>0 is a correction value, i.e. the center O of the first locus circle1And a center O of a second locus circle2The distance of (c).
7. The concentric tube robot inverse kinematics solution method according to claim 6, wherein step S23 comprises: assuming T is the target location, all that needs to be determined is the vector O1T and vector O1The angle of rotation between G; suppose the vector O1T and vector O1The angle of rotation between G is theta2iThen there is
Figure FDA0002946633480000031
θ2iThe angle being in the range of [0,180 °];
Assuming that the target point is rotated within 0,180 deg. counterclockwise, the angle of rotation is positive, and vice versa, there is,
Figure FDA0002946633480000032
alpha represents the rotation angle theta capable of distinguishing the two side spaces2iThe coordinate value of the sign is determined.
CN201810752589.9A 2018-07-10 2018-07-10 Inverse kinematics solving method of concentric tube robot Active CN108972550B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810752589.9A CN108972550B (en) 2018-07-10 2018-07-10 Inverse kinematics solving method of concentric tube robot

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810752589.9A CN108972550B (en) 2018-07-10 2018-07-10 Inverse kinematics solving method of concentric tube robot

Publications (2)

Publication Number Publication Date
CN108972550A CN108972550A (en) 2018-12-11
CN108972550B true CN108972550B (en) 2021-05-04

Family

ID=64537692

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810752589.9A Active CN108972550B (en) 2018-07-10 2018-07-10 Inverse kinematics solving method of concentric tube robot

Country Status (1)

Country Link
CN (1) CN108972550B (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112828893B (en) * 2021-01-26 2022-05-13 南京理工大学 Multi-segment linear drive continuum mechanical arm dynamics modeling method
CN113246121B (en) * 2021-04-08 2022-04-08 浙江工业大学 Concentric tube robot control method based on Actor-Critic deep reinforcement learning

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2019509103A (en) * 2016-03-04 2019-04-04 コヴィディエン リミテッド パートナーシップ Inverse kinematics control system for robotic surgical system
CN108068108B (en) * 2016-11-16 2021-02-02 沈阳高精数控智能技术股份有限公司 Method for controlling motion of in-plane swinging telescopic robot
CN107045552B (en) * 2017-04-28 2020-04-24 齐鲁工业大学 Motion control method and control device of quadruped robot based on sine diagonal gait and rapid table look-up method
CN107097231A (en) * 2017-07-06 2017-08-29 哈尔滨工业大学深圳研究生院 A kind of concentric tube robot precise motion control method of view-based access control model servo
CN107943056B (en) * 2017-12-25 2020-09-25 成都信息工程大学 Incomplete constraint wheeled robot track tracking control method based on table lookup method

Also Published As

Publication number Publication date
CN108972550A (en) 2018-12-11

Similar Documents

Publication Publication Date Title
US10682758B2 (en) Numerical method for obtaining the inverse kinematics of six-degree-of-freedom serial robot with an offset wrist
US11667035B2 (en) Path-modifying control system managing robot singularities
CN108972550B (en) Inverse kinematics solving method of concentric tube robot
CN110861088A (en) Motion optimization method of redundant degree of freedom robot
CN111168719B (en) Robot calibration method and system based on positioning tool
CN113632135A (en) System and method for low latency, high performance pose fusion
CN105382835A (en) Robot path planning method for passing through wrist singular point
CN109176526A (en) A kind of three axis Cartesian robot space circular arc interpolation methods
CN108062071B (en) Real-time measuring method for parameter curve track servo contour error
WO2022193668A1 (en) Inverse solution method and apparatus for arm angle interval of mechanical arm, and terminal device
CN104842371A (en) Robot hand-eye calibration method based on non-minimized and non-optimized algorithm
CN110345936A (en) The track data processing method and its processing system of telecontrol equipment
Luo et al. Deep motion network for freehand 3d ultrasound reconstruction
Milenkovic Continuous path control for optimal wrist singularity avoidance in a serial robot
CN112356032B (en) Posture smooth transition method and system
CN113334381A (en) Movable decoupling continuum robot control method
Wang et al. Learning-based visual-strain fusion for eye-in-hand continuum robot pose estimation and control
CN103412337A (en) Dual-independent rotating magnetic bar electromagnetic tracking-based position tracking method and electromagnetic tracking system
CN107505838B (en) A kind of method of accurate control machine fish C-shaped turning motion
Yuan et al. Static modeling and analysis of continuum surgical robots
CN115424701A (en) Bone surface follow-up technology about optimal path planning
Du et al. Continuum robots: Developing dexterity evaluation algorithms using efficient inverse kinematics
CN111618853B (en) Continuous robot deviation correction kinematics equivalent method and control device
CN105629880A (en) Inverse kinematics numerical solution for series manipulator with triangular apocenter mechanism
JP2021122894A (en) Inverse kinematics calculation device and inverse kinematics calculation method

Legal Events

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