## Abstract

The insertion of the needle is difficult because the deformation and displacement of the organs are the key elements in the surgical act. Liver and tumor modeling are essential in the development of the needle insertion model. The role of the needle is to deliver into the tumor an active chemotherapeutic agent. We describe in this chapter the deformation of the needle during its insertion into the human liver in the context of surgery simulation of the high- robotic-assisted intraoperative treatment of liver tumors based on the integrated imaging-molecular diagnosis. The needle is a bee barbed type modeled as a flexible thread within the framework of the Cosserat (micropolar) elasticity theory.

### Keywords

- bee needle
- human liver
- Cosserat elasticity theory

## 1. Introduction

The flexible bee needles are useful tools to transport drugs into the liver tumors [1, 2]. The insertion trajectory of the needle must avoid the ribs, blood vessels, and other organs to protect the liver [3, 4, 5, 6] (Figure 1a). The bee needle assures reduced insertion forces and small tissue deformations because of the tip deflections. The furthermost current publications on the surgical needle navigation into the liver can be demonstrated in [7, 8, 9]. The bee needle is shown in Figure 1b. The front angle has 157 deg., the back angle, 110 deg., the height is 0.5 mm, and the tip thickness 0.15 mm.

A number of scientific researches have been carried out on the collision free trajectory of the needle to the target. The surgical event requires experience in imaging the tumor location based on the liver structure and the microstructural interaction between the needle and the liver. Several studies have revealed that the needle flexibility is essential to achieve a good precision in the handling.

The strain and stress fields and the topological changes of the liver are not to be neglected during the needle navigation towards the tumor [10, 11, 12, 13]. Details of the forces during needle insertion into the liver are find in [14], the real time collision detection for virtual surgery in [15] and the minimal hierarchical collision detection in [16]. Optimization is required to modify the needle trajectory in order to protect the liver [17, 18], to manage the tumor risk [19], and to change the robot architecture [20, 21, 22]. The inverse sonification problem for capturing hardly detectable details in a medical image is treated in [23], and the control in [24, 25, 26, 27]. Microscopic investigation of the human liver offers details on its microanatomy with emphases to the granular, fibrillar components and irregular solid–fluid interfaces [28, 29, 30]. The basic unit of the liver is the hepatic lobule which is a hexagonal element with comprised the portal triad -portal vein, hepatic artery and the bile duct [31, 32]. Lobuli form two layers membranes with internal space of 100A and the cellular elements with twisted, spiraling fibers braided into the helical and screw-shaped gaps (pores) of 40–100

In this chapter we try to answer a few questions such as how is the deformation of the needle and how the free-collision trajectories are determined.

## 2. Deformation of the needle

The needle is a bee barbed needle and it is modeled as a flexible thread within the framework of the Cosserat (micropolar) elasticity theory [37, 38, 39, 40, 41, 42]. The Cosserat elasticity is applied to describe the interaction between the needle and the human liver. Let us consider a serial surgical robot composed of a revolute joint and a flexible needle. A Lagrange frame

The Euler axes are oriented with respect to the Lagrange axes by Euler angles

The functions

where * s* which is the coordinate along the central line of the needle. The functions

The function

So, the needle is rigid along the tangential direction and the total length of the needle

The link between the position vector

We introduce the inertia of the needle characterized by

where _{0} the area of the cross section,

The equations which describe the deformation are

where * A* and

*are the bending stiffness and respectively the torsional stiffness of the needle, related to the Lame constants*C

*is the radius of the cross section of the needle, and*a

The Eqs. (7–11) are solved by using the cnoidal method [43].

In short, this method is reducible to a generalization of the Fourier series with the cnoidal functions as the fundamental basis function. This is because the cnoidal functions are much richer than the trigonometric or hyperbolic functions, that is, the modulus

To understand the cnoidal method, consider now a nonlinear system of equations that govern the motion of a dynamical system

with

where

This system of equations can be reduced to Weierstrass equations of the type

We introduce the function transformation

where the theta function

and

Further, we write the solution (15) under the form

for

In (20) we recognize the expression [43].

with

The second term

If

As a result, the cnoidal method yields to solutions consisting of a linear superposition and a nonlinear superposition of cnoidal waves.

Therefore, by applying the cnoidal method, the closed form solutions of the Euler angles

where

with

Our objective is to determine the functions which measure the bending of the needle (

Figures 5 and 6 show that the deformation of the needle for both routes. We see that the deformation is small with no tendency to chaos. The strains are described by localized solitons which propagate for a long time without changes. The soliton is a localized wave with an infinite number of degrees of freedom. This wave conserves its properties even after interaction with another wave. In short, this wave acts somewhat like particles [43]. The system of Eqs. (7–11) has unique properties. These properties are locally preserved such as an infinite number of exact solutions expressed in terms of the Jacobi elliptic functions or the hyperbolic functions, and the simple formulae for nonlinear superposition of explicit solutions.

## 3. Determination of the free-collision trajectories

Let us present a model for determining the collision-free trajectories by using the Fibonacci sequence. The trajectories are determined from the restrictions of avoiding the collisions with blood vessels, ribs and surrounding tissues, and also the interference of needles with each other. We consider that more needles are planned to be inserted into the liver [45].

The trajectory of each needle

The sequence terms in the Fibonacci problem are 1, 1, 2, 3, 5, 8, 13, 21, 34, 55, 89, 144. It is clear that each term is a sum of two preceding sequence term. In fact, the sequence can be recursively defined in the form

The Fibonacci sequence was highlighted in nature for example in the arrangement of the flower petals, in seeds and in the spiral arrangement of pine cones and pineapple. The keys on a piano are divided into Fibonacci numbers, and numerous classical compositions implement the golden section. Such an example is found in the Alleluia Choir in Handel’s Messiah and in many of Chopin’s preludes [50].

Development of a target drug delivery technique usually consisted of three steps. In step 1, a 2D ultrasound image of the tumor is obtained. The size and position of the tumor are analyzed by the surgeon who decides the number of required needles. In step 2, the insertion points on the skin of the needles and the positions of the targets in the tumor are established. In step 3, once the needles are placed at the target points on the skin, they are guided according to a precise surgical planning based on the optimization of the collision-free trajectories to avoid the ribs, blood vessels and tissues in the abdominal area.

The base coordinate frame is established on the first needle with the vertical axis

The first control parameter is the length of the

and

where

The choice of the scaling number

The possible collision point between the needle and the tissue is analyzed by an identifier to check the minimum distance between needle and the surrounding tissue [45]. The minimum distance is expressed as

with

where

As application, the case of a tumor with a difficult location in the vicinity of the portal tree of the vascular territory in the liver, is considered (Figure 8a). The tumor image seen on the microscope is shown in Figure 8b. White and gray denote forbidden areas while the shade of purple are safe regions. The tumor is drawn in red. The Fibonacci algorithm is applied to three-needles with restrictions to avoid the collision with the tissues, blood vessels, ribs and previously inserted needles.

The task of our simulation is to determine the boundaries of each needle as a collision-free surface which represents the feasible insertion area based on given constraints. Then, the optimal trajectory of each needle can be chosen in this surface automatically.

Once all needles are placed at the predetermined epidermis and the ordering of entry is chosen to be 1, 2 and 3, following the first needle’s insertion, the operation is repeated for the needles 2 and 3.

Simulation of the collision-free trajectories for three needles is presented in Figure 9. The insertion scheme is determined by the Fibonacci spirals. A set of free-collision trajectories (red) in the immediate vicinity of the epidermis, is suggested.

From these possible collision-free trajectories (red) the green paths corresponding to the Fibonacci spirals (black) are chosen. These trajectories avoid the blood vessels (purple) and the coasts (brown) in all directions until the tumor. Optimal solution with 3 collision-free trajectories to the target is displayed in Figure 10. The Fibonacci spirals with the centers in ribs (brown) are displayed for needles 1 and 3. For the needle 2, the Fibonacci.

Three locally optimal collision-free trajectories for the surgical needle corresponding to three different entry points into the skin A, B and C are displayed in Figure 11.

## 4. Conclusions

The study investigates the navigation of a flexible needle into the human liver. The role of the needle is to deliver into the tumor an active chemotherapeutic agent. The deformation of the needle during its insertion into the human liver is describe in this chapter in the context of intraoperative treatment of liver tumors based on the integrated imaging-molecular diagnosis. The needle is a bee barbed type modeled as a flexible thread within the framework of the Cosserat (micropolar) elasticity theory. The Cosserat elasticity describes the interaction between the needle and the human liver by incorporating the local rotation of points and the couple stress as well as the force stress representing the chiral properties of the human liver.

## Acknowledgments

This work was supported by a grant of the Romanian Ministry of Research and Innovation project PN-IIIP2-2.1-PED-2019-0085 CONTRACT 447PED/2020 (Acronim POSEIDON).

## Conflict of interest

The authors of this paper certify that they have no affiliations with or involvement in any organization or entity with any financial or nonfinancial interest in the subject matter or materials discussed in this manuscript.

## References

- 1.
Chiroiu, V., Nedelcu, N., Munteanu, L., Rugina, C., Ionescu, M., Dragne, C., Modeling the flexible needle insertion into the human liver, Proceedings of the Romanian Academy, series A: Mathematics, Physics, Technical Sciences, Information Science (2020) - 2.
Sahlabadi, M., Hutapea, P., Tissue defomation and insertion force of bee-stinger inspired surgical needles , Journal of medical device, vol.12, 034501, 1–3 (2018) - 3.
Dionigi, R., Medical Intelligence Unit :Recent Advances in Liver Surgery , Landes Bioscience Austin, Texas (2009) - 4.
Dremin, V., Potapova, E., Zherebtsov, E., Kandurova, K., Shupletsov, V., Alekseyev, A., Mamoshin, A., Dunaev A ., Optical percutaneous needle biopsy of the liver: a pilot animal and clinical study , Scientific Reports, volume 10, Article number: 14200 (2020) - 5.
Abolhassani, N., Patel, R., Moallem, M., Control of soft tissue deformation during robotic needle insertion , Minimally Invasive Therapy, 15(3), 165–176 (2006) - 6.
Abolhassani, N., Patel, R., Moallem, M., Needle insertion into soft tissue: A survey , Medical Engineerings & Physics, 29, 413–431 (2007) - 7.
Torzilli, G., Minagawa, M., Takayama, T., Accurate preoperative evaluation of liver mass lesions without fine-needle biopsy, Hepatology, 30(4), 889–893 (1999) - 8.
Baker, N.E., Emerging mechanisms of cell competition , Nature Reviews Genetics (2020) - 9.
Petrowsky, H., Fritsch, R., Guckenberger, M., De Oliveira, M. L., Dutkowski, P., Clavien, P.A., Modern therapeutic approaches for the treatment of malignant liver tumours , Nature Reviews Gastroenterology & Hepatology (2020) - 10.
Gino van den Bergen, Collision detection in interactive 3D environments , Elsevier (2004) - 11.
Okamura, A.M., Simone, C., O’Leary, M.D., Force modeling for needle insertion into soft tissue , IEEE Trans Biomed Eng., 51,1707–16 (2004) - 12.
DiMaio, S.P., Salcudean, S.E., Needle insertion modeling and simulation , IEEE Trans Robot Automat., 19, 864–75 (2003) - 13.
DiMaio, S.P., Salcudean, S.E., Simulated interactive needle insertion , Proceedings of the 10th IEEE Symposium on Haptic Interfaces for Virtual Environment & Teleoperator Systems, 344–51 (2002) - 14.
Maurin, B., Barbe, L., Bayle, B., Zanne, P., In vivo study of forces during needle insertions , Scientific Workshop on Medical Robotics, Navigation and Visualization (MRNV04), Germany, Remagen, 415–422 (2004) - 15.
Lombardo, J-C, Cani, M-P, Neyret, F., Real time collision detection for virtual surgery , Computer Animation (CA’99), May 1999, Geneva, Switzerland, pp.82–90 (1999) - 16.
Zachmann, G., Minimal hierarchical collision detection, Proc. ACM Symposium on Virtual Reality Software and Technology (VRST), Hong Kong, China, 121–128 (2002) - 17.
Brișan, C., Boantă, C., Chiroiu, V., Introduction in optimisation of industrial robots. Theory and applications, Editura Academiei, Bucharest (2019) - 18.
Kataoka, H., Washio, T., Audette, M., Mizuhara, K., A model for relations between needle deflection, force, and thickness on needle insertion , Proceedings of the Medical Image Computing and Computer-Assisted Intervention Conference, 966–974 (2001) - 19.
Pisla D., Vaida, C., Birlescu I., Nadim, A.H., Gherman, B.,Corina Radu, Plitea N., Risk Management for the Reliability of Robotic Assisted Treatment of Non-resectable Liver Tumors , Appl. Sci., 10(1), 52 (2020) - 20.
Birlescu I., Manfred, H., Vaida C., Plitea N., Nayak A., Pisla, D., Complete Geometric Analysis Using the Study SE(3) Parameters for a Novel, Minimally Invasive Robot Used in Liver Cancer Treatment , Symmetry, 11(12), 1491 (2019) - 21.
Vaida, C., Plitea, N., Pisla, D., Gherman, B., Orientation module for surgical instruments - a systematical approach , Meccanica, 48(1), 145–158 (2013) - 22.
Munteanu, L., Rugină, C., Dragne, C., Chiroiu, V., On the robotic control based on interactive activities of subjects, Proceedings of the Romanian Academy, series A : Mathematics, Physics, Technical Sciences, Information Science, 21(1), (2020) - 23.
Chiroiu, V., Munteanu, L., Ioan, R., Dragne, C., Majercsik, L., Using the Sonification for Hardly Detectable Details in Medical Images, Scientific Reports, 9, article number 17711 (2019) - 24.
Korayem, M.H., Nikoobin, A., Azimirad, V., Trajectory optimization of flexible link manipulators in point-to point motion , Robotica, 27, 825–840 (2009) - 25.
Chiroiu, V., Munteanu, L., Ioan, R., Mosneguțu, V., Girip, I., On the dL algorithm for controlling the hybrid systems , Acta Electronica, Special Issue Proceedings of the XXIXth SISOM, 60(1–2), 58–65, Mediamira Science Publishing (2019) - 26.
Chiroiu, V., Munteanu, L., Dragne, C., Știrbu, C., On the diferential dynamic logic model for hybrid systems , Acta Technica Napocensis - series: Applied Mathematics, Mechanics, and Engineering, 61(4), 2018 - 27.
Majercsik, L., On the dL control appliedto a Stewart platform with flexible joints, Romanian Journal of Mechanics, 4(1), 27–38 (2019) - 28.
Motta, PM, The three-dimensional fine structure of the liver as revealed by scanning electron microscopy , Int. Rev. Cytol.,Suppl 6, 347–399 (1977) - 29.
Ma, M.H., Biempica, L., The normal human liver cell , Am. J. Pathol., 62, 353–370 (1971) - 30.
Jackson RL, Morrisett JD, Gotto AM., Lipoprotein structure and metabolism, Physiol. Rev. 56, 259–316 (1976) - 31.
Saxena, R., Theise, N.D., Crawford, J.M., Microanatomy of the human liver-exploring the hidden interfaces , Hepatology, 30(6), 1339–46 (1999) - 32.
Si-Tayeb, K., Lemaigre, F.P., Duncan, S.A., Organogenesis and development of the liver , Dev. Cell.,18(2), 175–89 (2010) - 33.
Yi–Je Lim, Dhanannjay Deo, Tejinder P. Singh, Daniel B. Jones, Suvranu De, In Situ Measurement and Modeling of Biomechanical Response of Human Cadaveric Soft Tissues for Physics-Based Surgical Simulation, Surg Endosc., 23(6), 1298–1307 (2009) - 34.
Narayan, K.S., Steele, W.J., Busch, H., Evidence that the granular and fibrillar components of nucleoli contain 28 and 65 RNA, respectively , Exp. Cell. Res., 43, 483–492 (1966) - 35.
Habenschus, M.D., Nardini, V., Dias, L.G., Rocha, B.A., Barbosa Jr. F. , Moraes de Oliveira, A.R., In vitro enantioselective study of the toxicokinetic effects of chiral fungicide tebu , Ecotoxicol. Environ. Saf., 181, 96–105 (2019)c onazole in human liver microsomes - 36.
Takahashi, M., Takani, D., Haba, M., Hosokawa, M., Investigation of the chiral recognition ability of human carboxylesterase 1 using indomethacin esters , Chirality, 32(1), 73–80 (2020) - 37.
Cosserat, E., Cosserat, F., Theorie des corps deformable , Herman et Fils, Paris, 1909 - 38.
Eringen, A.C., Microcontinuum Field Theory , volume II. Fluent Media, Springer, New York (2001) - 39.
Eringen, A.C., Kafadar, C.B., Polar field theories. In A. C. Eringen,editor, Continuum Physics, volume IV, pages 1–75. Academic Press, NewYork (1976) - 40.
Iesan, D., Existence theorems in the theory of micropolar elasticity . International Journal of Engineering Science, 8, 777–791 (1970) - 41.
Iesan, D., Existence theorems in micropolar elastostatics , International Journal of Engineering Science, 9, 59–78 (1971) - 42.
Forest, S., Sievert, R., Nonlinear microstrain theories , International Journal of Solids and Structures, 43(24), 7224–7245 (2006) - 43.
Munteanu, L., Donescu, St., Introduction to Soliton Theory :Applications to Mechanics, Book Series Fundamental Theories of Physics, vol.143, Kluwer Academic Publishers, Dordrecht, Boston (Springer Netherlands) (2004) - 44.
Chiroiu, V., Munteanu, L., Gliozzi, A.S., Application of the Cosserat theory for modeling the reinforcement carbon nanotube beams , CMC: Computers, Materials & Continua, 19(1), 1–16 (2010) - 45.
Dragne, C., Chiroiu, V., Munteanu, L., Brişan, C., Rugină, C., Ioan, R., Stănescu, N-D, Stan, A.F., On the collision free-trajectories of a multiple-needle robot based on the Fibonacci sequence, New Trends in Mechanism and Machine Science, Volume 89 of the Mechanisms and Machine Science series, Chapter 20 pp.1–12 Springer (2019) - 46.
Lai, A.C., Loreti, P., Velluci, P.: A Model for Robotic Hand Based on Fibonacci Sequence. In: Proceedings of the 11 ^{th}International Conference on Informatics in Control, Automation and Robotics (ICINCO-2014), pp.577–584 (2014) - 47.
Aghili, F., Parsa, K.: Design of a reconfigurable space robot with lockable telescopic joints . In Conference IEEE/RSJ, International Conference on Intelligent Robots and Systems (2006) - 48.
Burton, D. M.: Elementary number theory (5th ed.). New York: McGraw-Hill (2002) - 49.
Fox, W.P.: Fibonacci Search in Optimization of Unimodal Functions. Department of Mathematic Francis Marion University, Florence, SC 29501 (2002) - 50.
Silverman, J. H.: A friendly introduction to number theory . (3rd ed.) Upper Saddle River, NJ: Pearson Education (2006)