Design of a 4-DoF (degree of freedom) hybrid-haptic device for laparoscopic surgery

This paper presents a novel kinematics architecture with 4 DoFs (degrees of freedom) intended to be used as a haptic interface for laparoscopic surgery. The proposed architecture is a result of an association of serial and parallel kinematics chains, with each one handling a part of the whole device DoF. The serial chain allows one to handle the translation and self-rotation and the parallel chain handles the two tilt motions, and this in a disjoint way as the natural gesture of the surgeon. The proposed hybrid-haptic device (HH device) benefits from the split DoF to ensure a good kinematic performance, large workspace, as well as gravity compensation. The kinematics study of the HH device is presented and followed by the optimal dimensional synthesis and the gravity compensation model.


Introduction
Applications requiring interactions between humans and robots have grown in recent years. Among these applications, there may be mentioned robotic surgery and rehabilitation in the medical field. However, robots are powerful and rigid systems. This raises questions about the safety of using such robots with humans. Haptic can be an answer to these questions. In fact, for robotic surgery for instance, haptic is an additional feature giving more immersion to the surgeon. Some commercially available haptic devices are used to control surgical robots. For instance, Sigma.7 from Force Dimension ® is used in Tobergte et al. (2011) to control the MiroSurge robotic system designed by the German Aerospace Center (DLR). This device is hybrid since it is composed of a parallel part, based on the delta robot, and a spherical serial part. Also, a 3-DOF PHANTOM haptic device is used in Tavakoli et al. (2006) to control a surgery robot. The Phantom haptic device is augmented with 2 additional degrees of freedom to form a surgeon-robot interface with full force reflection capabilities. In Medellin-Castillo et al. (2020), the Omni Phantom from Sensable ® and the Falcon from Novit ® are used to evaluate the use of haptics and virtual reality technologies as an orthognathic surgery training tool. However, versatile and commercially haptic devices have a predefined workspace and degree of freedom, which are not optimized for the intended surgical operation. Therefore, several haptic devices developed specially for roboticassisted surgery are presented in the literature (van den Bedem et al., 2009;Saafi et al., 2017Saafi et al., , 2018Preault et al., 2018). The design of haptic devices is a challenging task since these devices have to be transparent and offer the best kinesthetic sense to the surgeon.
Minimally invasive surgery (MIS) requires mechanisms with at least 4 DoFs (degrees of freedom -three rotations around a fixed center and one translation) (Nisar et al., 2017). Therefore, spherical mechanisms are the first candidates for developing MIS haptic devices. van den Bedem et al. (2009) designed a 4-DoF haptic device with spherical serial architecture. However, all actuators were placed on the joint axes. This increases the actuator-required torques for compensating the gravity and also increases the weight of the end- effector. Saafi et al. (2017) used a spherical parallel manipulator as a structure for the haptic device. However, as parallel manipulators, this device suffers from the presence of the parallel singularity inside its workspace. Saafi et al. (2015) presented solutions to overcome the singularity issues such as actuator redundancy, developing and optimizing a new spherical parallel manipulator in Saafi et al. (2018). However, limitations have been reported in  regarding the limited motion of the self-rotation. Another kinematic base of the delta structure was studied in Preault et al. (2018). This kinematic solves the limitation of the self-rotation; however, the center of rotation (CoR) is at the bottom of the parallel mechanism. This increases the gravity compensation torques and the size of the required actuators.
In this paper, a new hybrid design of a haptic device is presented. This design connects a serial part and a parallel part around a fixed CoR. Hence, the gravity compensation torque is decreased since the weights of the two parts can be balanced. In addition, the self-rotation is decoupled from other DoFs and has the required range. This paper is organized as follows: Sect. 2 presents the minimally invasive surgery and the importance of teleoperation robotics. The new design of a 4-DoF hybrid-haptic device (HH device) is presented in Sect. 3. In Sect. 4, the geometric parameters of the 3-RRR (revolute-revolute-revolute) planar parallel manipulator (PPM), which is the parallel part of the HH device, are optimized. The haptic model is studied in Sect. 5. Section 6 determines the position of the 3-RRR PPM end-effector using sensors placed on the universal joint. The gravity compensation problem is modeled in Sect. 7. Finally, Sect. 8 concludes the paper.

Minimally invasive surgery
MIS is a laparoscopic surgical procedure where tools enter the patient's abdominal cavity through tiny incisions, as shown in Fig. 1. This procedure offers the patients the following advantages: shorter hospital stays, quicker recovery times, less pain and discomfort, less chance of infection  and bleeding, and much smaller scars. Nevertheless, MIS presents some drawbacks listed below: longer duration due to the complexity of the procedure and the physical tiredness of the surgeon; limited motions due to the limited workspace (abdominal cavity) and tool types; and need for in-depth experience due to the complexity and the 2D vision feedback.
To improve minimally invasive surgery and cope with the previous drawbacks, the use of teleoperation robotic systems has emerged. The tools are held by slave robots and the surgeon manipulates these robots using master devices as presented in Fig. 2. This improvement makes the surgeon more comfortable and decreases movement limitations as well as the duration of the surgical operation.
MIS limits the motion of tools inside the patient's body, as illustrated in Fig. 3. Only 4 DoFs are maintained and lim- ited to rotations around the incision point and one translation within the tool axis. This limitation makes the design of master and slave robots for MIS teleoperation systems challenging. In fact, the development and control of teleoperation systems for MIS are still up to date (Boabang et al., 2020;Iijima et al., 2020;Saracino et al., 2020). Since there are always possible improvements, many studies focus on the identification of the required workspace for MIS (Thakre et al., 2008;Pisla et al., 2008;Cavusoglu et al., 2001). In fact, the workspace size depends on the surgical procedure (Konietschke et al., 2003). Here, we adopt the workspace identified in Laribi et al. (2012) through surgeon gesture records using a motion capture system. The workspace is described by a cone with a half-apex angle of 26 • . The translation range is around 112 mm. The self-rotation depends on the dexterity level of the surgeon as well as on the used technique. Nevertheless, it is considered here to be between ± 90 • . The kinematic of the master device must offer a workspace larger than the required one to guarantee the efficiency of the proposed device. As highlighted in the introduction, the expert surgeon who tested the first prototype of the haptic device reported the issue of the self-rotation limitation. This has been taken into account in the new design. Peirs et al. (2004) measured the required force for MIS. The measurements were made in two directions: radial and axial. The radial force range is estimated at 1.7 N and the axial force at 2.5 N.
The next section deals with a new design of the hybridhaptic device for robotized MIS.

New hybrid spherical kinematics for the master device
As introduced in the previous section, the MIS-required motions are rotations around a fixed point, called the CoR, and a translation. A new kinematics offering these movements is presented in Fig. 4. The proposed architecture is obtained by associating serial and parallel kinematics chains, where each one handles a part of the whole device DoF. The serial chain, composed of a prismatic joint and a revolute joint, allows us to handle the translation as well as the self-rotation. The parallel chain handles the two tilt motions. The required DoFs are generated in a disjointed way; therefore, the two chains are connected through a universal joint. This association helps to cope with the drawbacks of serial architecture. This latter is not suitable for developing haptic devices since the actuators are placed on joint axes. In this case, the required gravity compensation torques are increased and the dynamic behavior of the device is not negligible. Therefore, a parallel chain is linked to the serial chain through the universal joint to ensure the force control for the two tilt motions. Here, a 3-RRR planar parallel mechanism with 2 DoFs is selected as a parallel structure for the following reasons: reducing the required actuator torques, having a singular-free workspace, and having the most compact workspace after optimization (Saafi and Lamine, 2020). The self-rotation is decoupled from the parallel part because it requires a wide range of motion. The computer-aided design (CAD) is presented in Fig. 5 to illustrate the appearance of the new architecture.
The haptic device measures the motion using sensors. Therefore, all joints of the serial part will be equipped with sensors (Fig. 6). This choice facilitates the resolution of the forward kinematic model of the master device. This model is required to measure the motion of the surgeon and then to control the surgical slave robot. The position of the endeffector can be expressed using the active angles of the new device as follows: where λ is the distance between the end-effector and the CoR (located at point O). β 1 and β 2 are the rotation angles of the universal joint in the X and Y directions, respectively. β 3 is the self-rotation angle. The orientations of the tool axis, Z E , of the haptic device are expressed using ZY Z Euler angles (ψ, θ , and ϕ): ψ and θ for the tilt and ϕ for the self-rotation.
In order to illustrate the kinematic performance of the new proposed design, the dexterity index is measured using the inverse conditioning number of the orientation Jacobian matrix, J s , of the serial part of the HH device. Its expression is as follows: The mechanism composed of the universal joint and the revolute joint can be considered a spherical wrist capable only of rotations around the center O. In this case, the Jacobian matrix J s relating the angular velocity of the moving platform and the active joint velocities is composed of the axes of the joints (Nouaille et al., 2010). The Jacobian matrix expression is as follows: where Y 1 is the second axis of the universal joint, as illustrated in Fig. 6. Figure 6. The serial chain geometric parameters. Figure 7 shows the dexterity distribution of the new haptic device in the plane (ψ, θ ). The new device presents good kinematic performances, since the dexterity is allowed to be higher than 0.6. The serial part does not require optimization since all joints are decoupled. However, the haptic device must be compact. It must fit in a box with a height equal to 500 mm and width equal to 200 mm. Only the parallel part requires optimization. The next section deals with the optimization of the 3-RRR planar parallel manipulator.

Optimization of the 3-RRR redundant planar parallel manipulator
The redundant 3-RRR PPM is composed of three identical branches, as illustrated in Fig. 8. Each branch is composed of two links: one connected to the base and the other connected to the end-effector. The end-effector position is defined by point E and its (x, y) coordinates. The three active joints are placed on the vertices of an equilateral triangle such that l 0 is the edge length. The base reference frame (O , X Y ) is placed at the center of the equilateral triangle. l 1 and l 2 are the sizes of the first link and the second link, respectively. θ 1 , θ 2 , and θ 3 are the active joint angles. This section focuses on the geometric optimization of the redundant PPM to fit the required MIS workspace.
As shown in Fig. 9, the projection of the required workspace in the PPM plane is a circle. The radius of this circle depends on the distance, denoted h, between the center of rotation (center of the universal joint, Fig. 4) and the PPM plane. The angle α is considered equal to 26°. The two distances h and l 0 are chosen manually to have a compact structure. The two geometric constants are equal to (h = 150 mm, l 0 = 200 mm) .
The radius of the projected circle is equal to the following expression: The border of the required workspace is shown in Fig. 10. The circle is discretized to 12 points defined as follows: The optimization design vector is  where x L1 and x L2 are the length of the first link and the second link, respectively. The optimization problem is formulated as follows: subject to P i ∈ W S, i = 1, . . ., 12.
The point P i is within the workspace of the PPM if and only if the following condition is met: where (x i , y i ) are the coordinates of P i (i = 1, · · ·, 12) and (x j , y j ) are the coordinates of the points (A, B, and C). The proposed optimization approach is based on the minimization of an objective function f (x L ) which is the quadratic sum of the robot links (x L1 , x L2 ). The optimization is subject to one constraint which involves the workspace (WS) and aims to guarantee that the haptic device WS fits the prescribed one. A genetic algorithm is used to minimize the objective function, where the constraints are nonlinear and the lower and upper bounds of x Li are [50,200] in millimeters. Furthermore, the initial guess is x init = [200,200] in millimeters. As a result, the computed optimal design vector is 99.8827, 88.7472) [mm].
(10) Figure 11 presents the workspace of the optimized PPM. Nevertheless, to move the serial singularity away from the border of the prescribed workspace, each component of x is increased by 5 mm. The following design vector was selected: (11) Figure 12 presents the workspace of the selected PPM with the prescribed area.
In this section, the geometrical parameters of the 3-RRR redundant PPM are optimized to obtain a compact structure. The next section focuses on the study of the haptic model.

Haptic model
Haptic devices apply forces and torques to the user's hand using actuators. As described before, for the proposed design, the self-rotation and the translation will each be equipped by a specific actuator. However, the universal joint will not be actuated. Instead, the optimized 3-RRR redundant PPM will be equipped with actuators. This section will focus on the development of the haptic model for the 3-RRR redundant PPM.
The kinematic model of the 3-RRR PPM is as follows: where J x is a 3 × 2 matrix called the parallel part of the Jacobian matrix, and J θ is a 3 × 3 diagonal matrix called the serial part of the Jacobian matrix.
where Table 1. Working mode signs.

Working modes Sign
For the non-redundant PPM, the relation between the force vector at the end-effector, F , and the actuated joint torques, τ , is as follows: where J is the Jacobian matrix. The gravity effects are neglected because of the horizontal placement of the PPM. For the redundant PPM, only the force vector can be expressed using the actuator joint torques as follows: where G is the inverse of the Jacobian matrix. The pseudo-inverse is used to express the torque vector τ using the force vector F as follows: where (G T ) + is the pseudo-inverse of G T , I is the 3 × 3 unit matrix, and V is an arbitrary vector. (I − (G T ) + G)V is the null-space torque which can change the actuator torques τ without affecting the force vector F . To get the optimal torque distribution, the following optimization approach is employed: subject to τ i < τ max , i = 1, 2, 3.
As shown in the haptic model, the active angles (θ 1 , θ 2 , and θ 3 ) are required to determine the Jacobian matrix. Since the 3-RRR PPM will not be equipped with sensors, the inverse kinematic model (IKM) will be used to determine the required angles. The next section studies the determination of the active angles using the sensors placed on the universal joint.

Determination of the 3-RRR PPM active angles
The IKM gives the active angles (θ 1 , θ 2 , θ 3 ) as a function of the 3-RRR PPM end-effector position (x, y). The expression of the active angles is as follows: where Each angle θ i (i = 1, 2, 3) has two possible solutions depending on the sign (±) in Eq. (20). The combination of all the solutions gives eight working modes (WMs) for the 3-RRR PPM enumerated in Table 1. For the redundant PPM shown in Fig. 8, WM1 is considered with the sign (+) for the three angles' expressions.
The IKM is not enough to determine the active angles since the position of the PPM end-effector is not known. Figure 13 shows the approach used to determine the position of the PPM end-effector, E, defined by the coordinates (x, y). The sensors measuring the tilt angles are placed on the universal joint at the center of rotation of the haptic device.
The PPM end-effector position can be expressed using ψ and θ as follows: x = −h tan θ cos ψ, y = −h tan θ sin ψ.  As detailed in the previous section, the gravity effect of the redundant 3-RRR PPM is neglected since it will be placed horizontally. However, the haptic device requires gravity compensation. The next section deals with the gravity compensation for the new haptic device.

Gravity compensation
The hybrid-haptic device can be modeled as illustrated in Fig. 14. The upper part is defined by its center of gravity G 1 and the mass m 1 and the lower part by its center of gravity G 2 and the mass m 2 . The upper part is composed of the gripper handle as well as prismatic and revolute joints, while the main component of the lower part is the passive prismatic joint. The upper part is expected to be heavier than the lower part that will compensate the gravity partially by the counterbalance effect.
The remaining part is compensated by the redundant 3-RRR PPM by generating the force F g .
The gravity compensation problem can be formulated as follows: Figure 16. Gravity compensation-required torques. where P 1 and P 2 are the gravity vectors of the upper and lower blocks, respectively. To reduce the required force F g to compensate the gravity, the two blocks must be balanced as much as possible. This issue will be handled when designing the prototype of the master device.
A numerical simulation of the gravity compensation model is carried out to calculate the required actuator torques. Table 2 presented the gravity compensation model parameters.
The compensation force F g is given by the following expression: It is obvious that the compensation force is maximum when the tilt angle θ is maximum. Therefore, the simulation is made with θ equalling θ max = 26 • . The angle ψ is varied between 0 • and 360 • . The compensation force F g , calculated using Eq. (24), is shown in Fig. 15. The maximal gravity compensation force is equal to 2.6 • N. The required actuator torques for gravity compensation are illustrated in Fig. 16. The maximum required actuated torque for compensating the gravity is equal to 0.2 Nm.

Conclusions
In the present work, a new hybrid spherical manipulator with 4 DoFs was presented and studied. The proposed manipulator, named the hybrid-haptic device (HH device), is intended to be used as a haptic device for laparoscopic surgery. The HH device is an association of two kinematic chains, serial and parallel. As the natural gesture of the surgeon, the HH device allowed them to handle disjointly the mobilities two by two, the two tilt motions and the translation and selfrotation. The serial chain is built out with installed sensors on its joints to compute the orientation and the position of the end-effector. However, the parallel chain, which is a redundant 3-RRR PPM with 2 DoFs, is used to ensure the force feedback and to compensate the gravity issue. The geometric parameters of the PPM are optimized to fit the prescribed workspace. The haptic and compensation models were discussed. The CAD prototype is presented and used to validate the proposed kinematic. Future works will focus on sensors and actuator selections, the development of a first prototype of the HH device, as well as the possible use of the same kinematics for the slave surgical robot.
Data availability. All the data used in this paper can be obtained from the corresponding author upon request.
Author contributions. HS designed the new master device. HS and MAL co-wrote the paper. The research work was supervised by SZ.