Figure 1.
Characterization of the remote center of motion (RCM). The motion of the surgical tool is constrained by the trocar point, resulting in an RCM error, which is defined as the minimum distance between the tool axis and the trocar point .
Figure 1.
Characterization of the remote center of motion (RCM). The motion of the surgical tool is constrained by the trocar point, resulting in an RCM error, which is defined as the minimum distance between the tool axis and the trocar point .
Figure 2.
Three kinematic chains were used for numerical validation in simulation. (A) A 6-DOF endoscope holder. (B) A 7-DOF manipulator holding a 3-DOF robotic surgical tool. (C) A 7-DOF manipulator holding a 5-DOF robotic surgical tool.
Figure 2.
Three kinematic chains were used for numerical validation in simulation. (A) A 6-DOF endoscope holder. (B) A 7-DOF manipulator holding a 3-DOF robotic surgical tool. (C) A 7-DOF manipulator holding a 5-DOF robotic surgical tool.
Figure 3.
Average solving time (ms) of successful IK solves with unconstrained random targets for each kinematic chain and IK implementation. Considering all evaluated kinematic chains, INVJ and INVJ+HQP showed the best performance, with solving times of less than 0.2 ms.
Figure 3.
Average solving time (ms) of successful IK solves with unconstrained random targets for each kinematic chain and IK implementation. Considering all evaluated kinematic chains, INVJ and INVJ+HQP showed the best performance, with solving times of less than 0.2 ms.
Figure 4.
Simulation environment for KC-1 (endoscope holder) tracking a 10 cm circular path. The trajectory followed by the endoscope tip is visualized in blue, and snapshots (a–c), shown on the right, depict the robot executing the tracking task at different trajectory steps.
Figure 4.
Simulation environment for KC-1 (endoscope holder) tracking a 10 cm circular path. The trajectory followed by the endoscope tip is visualized in blue, and snapshots (a–c), shown on the right, depict the robot executing the tracking task at different trajectory steps.
Figure 5.
Average runtime (ms) of successful IK solves for a constrained tracking task with KC-1 (endoscope holder) kinematic chain. The best time performance is achieved by the INVJ-based IK solvers with average solving times of 0.07 ms.
Figure 5.
Average runtime (ms) of successful IK solves for a constrained tracking task with KC-1 (endoscope holder) kinematic chain. The best time performance is achieved by the INVJ-based IK solvers with average solving times of 0.07 ms.
Figure 6.
Simulation environment for KC-2 Kinematic chain (OpenRST) tracking a 10 cm circular path. In the close-up view, the tool tip target orientation is shown. The trajectory followed by the surgical tool tip is visualized in blue, and snapshots (a–c), shown on the right, depict the robot executing the tracking task at different trajectory steps.
Figure 6.
Simulation environment for KC-2 Kinematic chain (OpenRST) tracking a 10 cm circular path. In the close-up view, the tool tip target orientation is shown. The trajectory followed by the surgical tool tip is visualized in blue, and snapshots (a–c), shown on the right, depict the robot executing the tracking task at different trajectory steps.
Figure 7.
Average runtime (ms) of successful IK solves for a constrained tracking task with KC-2 (OpenRST) kinematic chain. The concurrent solvers demonstrated slightly better performance than their respective optimization-based single-method IK solvers, with INVJ showing the worst time performance due to its low solve rate.
Figure 7.
Average runtime (ms) of successful IK solves for a constrained tracking task with KC-2 (OpenRST) kinematic chain. The concurrent solvers demonstrated slightly better performance than their respective optimization-based single-method IK solvers, with INVJ showing the worst time performance due to its low solve rate.
Figure 8.
Simulation environment for KC-3 Kinematic chain (Hyper-redundant RST) tracking a 10 cm circular path. In the close-up view, the tool tip target orientation. The trajectory followed by the surgical tool tip is visualized in blue, and snapshots (a–c), shown on the right, depict the robot executing the tracking task at different trajectory steps.
Figure 8.
Simulation environment for KC-3 Kinematic chain (Hyper-redundant RST) tracking a 10 cm circular path. In the close-up view, the tool tip target orientation. The trajectory followed by the surgical tool tip is visualized in blue, and snapshots (a–c), shown on the right, depict the robot executing the tracking task at different trajectory steps.
Figure 9.
Average runtime (ms) of successful IK solves for a constrained tracking task with KC-3 (Hyper-redundant RST) kinematic chain. The concurrent solvers demonstrated slightly better performance than their respective optimization-based single-method IK solvers, with INVJ+HQP showing the best average solving time performance.
Figure 9.
Average runtime (ms) of successful IK solves for a constrained tracking task with KC-3 (Hyper-redundant RST) kinematic chain. The concurrent solvers demonstrated slightly better performance than their respective optimization-based single-method IK solvers, with INVJ+HQP showing the best average solving time performance.
Figure 10.
Endoscope positioning experiment (A) Experimental setup. A 6-DOF manipulator (VS 050, Denso Corp.) with a rigid endoscope mounted. (B) The RCM constraint and the helix trajectory followed for each IK solver.
Figure 10.
Endoscope positioning experiment (A) Experimental setup. A 6-DOF manipulator (VS 050, Denso Corp.) with a rigid endoscope mounted. (B) The RCM constraint and the helix trajectory followed for each IK solver.
Figure 11.
(A) Average IK solving time for the endoscope positioning task. (B) RCM error for each IK solver obtained from motion capture for each step of the helix trajectory.
Figure 11.
(A) Average IK solving time for the endoscope positioning task. (B) RCM error for each IK solver obtained from motion capture for each step of the helix trajectory.
Figure 12.
Snapshots of a robotic endoscope holder tracking a 4-DOF helix path subjected to an RCM constraint. The labels (a–h) indicate the sequence of the robot’s movement along the path.
Figure 12.
Snapshots of a robotic endoscope holder tracking a 4-DOF helix path subjected to an RCM constraint. The labels (a–h) indicate the sequence of the robot’s movement along the path.
Figure 13.
Endoscope positioning experiment (A) A 7-DOF manipulator (Gen3, Kinova) with a 3-DOF robotic surgical tool attached to the end effector. (B) The RCM constraint and the Lissajous trajectory followed for each IK solver.
Figure 13.
Endoscope positioning experiment (A) A 7-DOF manipulator (Gen3, Kinova) with a 3-DOF robotic surgical tool attached to the end effector. (B) The RCM constraint and the Lissajous trajectory followed for each IK solver.
Figure 14.
(A) Average IK solving time for the surgical tool pose control task. (B) RCM error for each IK solver obtained from motion capture for each step of the Lissajous trajectory.
Figure 14.
(A) Average IK solving time for the surgical tool pose control task. (B) RCM error for each IK solver obtained from motion capture for each step of the Lissajous trajectory.
Figure 15.
Snapshots of the robotic surgical tool tracking a 6-DOF Lissajous path subjected to an RCM constraint. The labels (a–h) indicate the sequence of the robot’s movement along the path.
Figure 15.
Snapshots of the robotic surgical tool tracking a 6-DOF Lissajous path subjected to an RCM constraint. The labels (a–h) indicate the sequence of the robot’s movement along the path.
Table 1.
Parameters used for the unconstrained inverse kinematics benchmark.
Table 1.
Parameters used for the unconstrained inverse kinematics benchmark.
Parameter | Value |
---|
EE task error type | log6 |
Max. EE task error | 1 × 10 |
Max. computation time (s) | 10 × 10 |
INVJ | 1.0 |
INVJ | 1.0 |
NLO EE task weight () | 20 |
NLO RCM task weight () | 100 |
NLO weight () | 1 × 10 |
HQP | 1.0 |
HQP | 1.0 |
Table 2.
A comparison of solve rate and solving time for a total of nine IK solvers: five single-method and four concurrent. The targets are randomly selected from reachable poses. Evaluation is performed for three kinematic chains commonly found in robot-assisted surgical applications: KC-1 (Endoscope Holder), KC-2 (Manipulator + OpenRST), and KC-3 (Manipulator + Hyper-redundant RST). The highest values in each column are highlighted in bold.
Table 2.
A comparison of solve rate and solving time for a total of nine IK solvers: five single-method and four concurrent. The targets are randomly selected from reachable poses. Evaluation is performed for three kinematic chains commonly found in robot-assisted surgical applications: KC-1 (Endoscope Holder), KC-2 (Manipulator + OpenRST), and KC-3 (Manipulator + Hyper-redundant RST). The highest values in each column are highlighted in bold.
| KC-1: Endoscope | KC-2: OpenRST | KC-3: Hyper RST | All |
---|
IK Algorithm | Solve Rate (%) | Avg. Time (ms) | Solve Rate (%) | Avg. Time (ms) | Solve rate (%) | Avg. Time (ms) | Avg. Solve rate (%) | Avg. Time (ms) |
---|
KDL | 35.5 | 0.243 | 81.2 | 0.717 | 68.1 | 0.941 | 61.7 | 0.708 |
TRAC-IK | 99.9 | 0.379 | 100 | 0.359 | 99.7 | 0.649 | 99.8 | 0.462 |
INVJ | 99.7 | 0.267 | 99.9 | 0.144 | 100 | 0.150 | 99.8 | 0.187 |
NLO | 47.3 | 1.913 | 98.8 | 3.745 | 93.3 | 5.477 | 79.8 | 4.058 |
HQP | 43.4 | 0.941 | 94.4 | 1.648 | 93.3 | 2.082 | 77.0 | 1.690 |
INVJ+NLO | 98.5 | 1.191 | 99.6 | 0.566 | 99.9 | 0.729 | 99.4 | 0.826 |
INVJ+HQP | 99.9 | 0.265 | 100 | 0.148 | 100 | 0.161 | 99.9 | 0.190 |
HQP+NLO | 57.7 | 1.696 | 96.8 | 2.872 | 98.2 | 2.924 | 84.2 | 1.831 |
INVJ+NLO+HQP | 98.9 | 4.038 | 99.5 | 5.063 | 100 | 5.584 | 99.5 | 4.685 |
Table 3.
Parameters used for the constrained inverse kinematics benchmark.
Table 3.
Parameters used for the constrained inverse kinematics benchmark.
Parameter | Value |
---|
EE task error type | log6 1 |
Max. EE task error | 1 × 10 |
Max. RCM error | 1 × 10 |
Max. computation time (s) | 10 × 10 |
INVJ | 1.0 |
INVJ | 1.0 |
NLO EE task weight () | 20 |
NLO RCM task weight () | 100 |
NLO weight () | 1 × 10 |
HQP | 1.0 |
HQP | 1.0 |
Table 4.
A comparison of the solve rate and the solving time for single-method and concurrent IK algorithms. The target trajectories are two circular paths with diameters of 3 cm and 10 cm. Evaluation is performed for the KC-1 kinematic chain (endoscope holder).
Table 4.
A comparison of the solve rate and the solving time for single-method and concurrent IK algorithms. The target trajectories are two circular paths with diameters of 3 cm and 10 cm. Evaluation is performed for the KC-1 kinematic chain (endoscope holder).
| KC-1: 6-DOF Endoscope Holder |
---|
| Circumference D = 3 cm | Circumference D = 10 cm |
---|
IK Algorithm | Solve Rate (%) | Avg. Time (ms) | Solve Rate (%) | Avg. Time (ms) |
---|
INVJ | 100 | 0.044 | 100 | 0.072 |
NLO | 100 | 0.668 | 52 | 0.829 |
HQP | 100 | 0.449 | 100 | 0.578 |
INVJ+NLO | 100 | 0.055 | 100 | 0.072 |
INVJ+HQP | 100 | 0.045 | 100 | 0.068 |
Table 5.
A comparison of the solve rate and the solving time for single-method and concurrent IK algorithms. The target trajectories are two circular paths with diameters of 3 cm and 10 cm. Evaluation is performed for the KC-2 (Manipulator + OpenRST) kinematic chain.
Table 5.
A comparison of the solve rate and the solving time for single-method and concurrent IK algorithms. The target trajectories are two circular paths with diameters of 3 cm and 10 cm. Evaluation is performed for the KC-2 (Manipulator + OpenRST) kinematic chain.
| KC-2: 7-DOF Manipulator + 3-DOF RST |
---|
| Circumference D = 3 cm | Circumference D = 10 cm |
---|
IK Algorithm | Solve Rate (%) | Avg. Time (ms) | Solve Rate (%) | Avg. Time (ms) |
---|
INVJ | 43 | 6.729 | 25 | 4.843 |
NLO | 100 | 1.455 | 99 | 2.523 |
HQP | 100 | 1.983 | 100 | 2.965 |
INVJ+NLO | 100 | 1.264 | 100 | 2.390 |
INVJ+HQP | 100 | 1.776 | 100 | 3.037 |
Table 6.
A comparison of the solve rate and the solving time for single-method and concurrent IK algorithms. The target trajectories are two circular paths with diameters of 3 cm and 10 cm. Evaluation is performed for the KC-3 kinematic chain (Manipulator + Hyper-redundant RST).
Table 6.
A comparison of the solve rate and the solving time for single-method and concurrent IK algorithms. The target trajectories are two circular paths with diameters of 3 cm and 10 cm. Evaluation is performed for the KC-3 kinematic chain (Manipulator + Hyper-redundant RST).
| KC-3: 7-DOF Manipulator + 5-DOF RST |
---|
| Circumference D = 3 cm | Circumference D = 10 cm |
---|
IK Algorithm | Solve Rate (%) | Avg. Time (ms) | Solve Rate (%) | Avg. Time (ms) |
---|
INVJ | 96 | 1.614 | 69 | 4.136 |
NLO | 100 | 1.553 | 100 | 2.015 |
HQP | 99 | 0.746 | 100 | 2.478 |
INVJ+NLO | 100 | 1.344 | 100 | 1.848 |
INVJ+HQP | 99 | 0.708 | 100 | 2.366 |
Table 7.
A comparison of solve rate and time performance for single and concurrent IK algorithms considering all target paths.
Table 7.
A comparison of solve rate and time performance for single and concurrent IK algorithms considering all target paths.
IK Algorithm | Solve Rate (%) | Avg. Time (ms) | RCM Error (m) | EE Error (m) |
---|
INVJ | 72 | 1.991 |
|
|
NLO | 91 | 1.564 |
|
|
HQP | 99 | 1.534 |
|
|
INVJ+NLO | 99 | 1.162 |
|
|
INVJ+HQP | 99 | 1.334 |
|
|
Table 8.
Solve rate and solving time performance for a 6-DOF robotic endoscope holder tracking a 4D helix path (D = 5 cm).
Table 8.
Solve rate and solving time performance for a 6-DOF robotic endoscope holder tracking a 4D helix path (D = 5 cm).
4D Helix Path Tracking (D = 5 cm) |
---|
IK Algorithm | Solve Rate (%) | Avg. Time (ms) | RCM Error (mm) |
---|
NLO | 8.0 | 3.362 | - |
HQP | 100 | 1.340 | 0.44 |
INVJ+NLO | 99.6 | 0.180 | 0.43 |
INVJ+HQP | 100 | 0.189 | 0.42 |
Table 9.
Solve rate and solving time performance for a 7-DOF robotic manipulator + 3 DOF RST tracking a 6D Lissajous path.
Table 9.
Solve rate and solving time performance for a 7-DOF robotic manipulator + 3 DOF RST tracking a 6D Lissajous path.
6D Lissajous Path Tracking |
---|
IK Algorithm | Solve Rate (%) | Avg. Time (ms) | Avg. RCM Error (mm) |
---|
NLO | 100 | 4.414 | 0.99 |
HQP | 100 | 3.774 | 0.95 |
INVJ+NLO | 100 | 3.321 | 0.97 |
INVJ+HQP | 99.9 | 2.758 | 0.94 |