Search for Optimal Parameters in the Control Structure of a Surgical System for Soft Tissue Operations Based on In Vitro Experiments on Cardiovascular Tissue

Abstract

The surgical robots currently used in cardiac surgery are equipped with a remote center of motion (RCM) mechanism that enables the required spherical workspace. The dynamics model of the surgical robot’s RCM mechanism presented in this work includes a direct current (DC) motor, an optimal proportional–integral–derivative (PID) controller, and a LuGre friction model that takes into account the Stribeck effect and surface deformation. A finite element method (FEM) analysis of transients was carried out using the energy hypothesis of von Mises with an optimal input signal from the mechatronic system with a PID controller obtained using the Runge–Kutta differentiation method in the Dormand–Prince ordinary differential equations variant (ODE45). Five criteria were adopted for the objective function: the safety factor related to the stress function in the time-varying strength problem, the first natural frequency related to stiffness and the resonance phenomenon, the buckling coefficient in the statics problem related to stability, the static factor of safety, and the displacement of the operating tip. The force inputs to the dynamics model were derived from in vitro force measurements on cardiovascular tissue using a force sensor. The normality of the statistical distribution of the experimental data was confirmed using the Kolmogorov–Smirnov statistical test. The problem of multi-criteria optimization was solved using the non-sorter genetic algorithm (NSGA-II), the finite element method, and the von Mises distortion energy hypothesis. Velocity input signals for the transient dynamics model were obtained from a second in vitro experiment on cardiovascular tissue using the minimally robotic invasive surgery (MIRS) technique. An experienced cardiac surgeon conducted the experiment in a modern method using the Robin Heart Vision surgical robot, and a system of four complementary metal–oxide–semiconductor (CMOS) optical sensors and ariel performance analysis system (APAS-XP 2002) software were used to obtain the endoscopic tool trajectory signal. The trajectory signal was accurate to ±2 [mm] in relation to the adopted standard, and it was smoothed using the Savitzky–Golay (SG) polynomial smoothing, whose parameters were optimally selected using the Durbin–Watson (DW) statistical test.