Research Papers

Mechanics-Based Interactive Modeling for Medical Flexible Needle Insertion in Consideration of Nonlinear Factors

[+] Author and Article Information
Shan Jiang

Centre for Advanced Mechanisms and Robotics,
School of Mechanical Engineering,
Tianjin University,
Tianjin 300072, China
e-mail: shanjmri@tju.edu.cn

Xingji Wang

Centre for Advanced Mechanisms and Robotics,
School of Mechanical Engineering,
Tianjin University,
Tianjin 300072, China

1Corresponding author.

Contributed by the Design Engineering Division of ASME for publication in the JOURNAL OF COMPUTATIONAL AND NONLINEAR DYNAMICS. Manuscript received November 18, 2013; final manuscript received May 30, 2015; published online June 30, 2015. Assoc. Editor: Carmen M. Lilley.

J. Comput. Nonlinear Dynam 11(1), 011004 (Jan 01, 2016) (11 pages) Paper No: CND-13-1290; doi: 10.1115/1.4030747 History: Received November 18, 2013; Revised May 30, 2015; Online June 30, 2015

A mechanics-based model of flexible needle insertion into soft tissue is presented in this paper. Different from the existing kinematic model, a new model has been established based on the quasi-static principle, which also incorporates the dynamics of needle motions. In order to increase the accuracy of the model, nonlinear characteristics of the flexible needle and the soft tissue are both taken into account. The nonlinear Winkler foundation model and the modified Euler–Bernoulli theory are applied in this study, providing a theoretical framework to study insertion and deformation of needles. Galerkin method and iteration cycle analysis are applied in solving a series of deformation control equations to obtain the needle deflection. The parameters used in the mechanics-based model are obtained from the needle force and needle insertion experiment. Sensitivity studies show that the model can respond reasonably to changes in response to variations in different parameters. A 50 mm needle insertion simulation and a 50 mm corresponding needle insertion experiment are conducted to prove the validity of the model. At last, a study on different needle tip bevel demonstrates that the mechanics-based model can precisely predict the needle deflection when more than one parameter is changed. The solution can also be used in optimizing trajectory of the needle tip, enabling the needle to reach the target without touching important physiological structures such as blood vessels with the help of dynamic trajectory planning.

Copyright © 2016 by ASME
Your Session has timed out. Please sign back in to continue.


Abolhassani, N., Patel, R., and Moallem, M., 2007, “Needle Insertion Into Soft Tissue: A Survey,” Med. Eng. Phys., 29(4), pp. 413–431. [CrossRef] [PubMed]
Park, W., Wang, Y., and Chirikjian, G. S., 2009, “Path Planning for Flexible Needles Using Second Order Error Propagation,” International Workshop on Algorithmic Foundations of Robotics, Guanajuato, Mexico, pp. 583–595. [CrossRef]
Kallem, V., and Cowan, N. J., 2009, “Image Guidance of Flexible Tip-Steerable Needles,” IEEE Trans. Robot., 25(1), pp. 191–196. [CrossRef] [PubMed]
Park, Y. L., Elayaperumal, S., Daniel, B., Ryu, S. C., Shin, M., Savall, J., Black, R. J., Moslehi, B., and Cutkosky, M. R., 2010, “Real-Time Estimation of 3-D Needle Shape and Deflection for MRI-Guided Interventions,” IEEE ASME Trans. Mech., 15(6), pp. 906–915. [CrossRef]
Majewicz, A., Marra, S., van Vledder, M., Lin, M., Choti, M., Song, D., and Okamura, A., 2011, “Behavior of Tip-Steerable Needles in Ex Vivo and In Vivo Tissue,” IEEE Trans. Biomed. Eng., 59(10), pp. 2705–2715. [CrossRef]
Misra, S., Reed, K. B., Schafer, B. W., Ramesh, K. T., and Okamura, A. M., 2010, “Mechanics of Flexible Needles Robotically Steered Through Soft Tissue,” Int. J Robot. Res., 29(13), pp. 1640–1660. [CrossRef]
Rucker, D. C., Jones, B. A., and Webster, R. J., 2010, “A Geometrically Exact Model for Externally Loaded Concentric-Tube Continuum Robots,” IEEE Trans. Robot., 26(5), pp. 769–780. [CrossRef] [PubMed]
Ahmadian, M. T., and Janabi-Sharifi, F., 2010, “Modeling, Simulation, and Optimal Initiation Planning for Needle Insertion Into the Liver,” ASME J. Biomech. Eng., 132(4), p. 041001. [CrossRef]
Avramidis, I. E., and Morfidis, K., 2006, “Bending of Beams on Three-Parameter Elastic Foundation,” Int J. Solids. Struct., 43(2), pp. 357–375. [CrossRef]
Taylor, R. L., Filippou, F. C., Saritas, A., and Auricchio, F., 2003, “A Mixed Finite Element Method for Beam and Frame Problems,” Comput Mech., 31(1–2), pp. 192–203. [CrossRef]
Naidu, N. R., and Rao, G. V., 1995, “Stability Behaviour of Uniform Columns on a Class of Two Parameter Elastic Foundation,” Comput Struct., 57(3), pp. 551–553. [CrossRef]
Reddy, J. N., 2004, Nonlinear Finite Element Analysis, Oxford University, New York.
Goksel, O., Dehghan, E., and Salcudean, S. E., 2009, “Modeling and Simulation of Flexible Needles,” Med. Eng. Phys., 31(9), pp. 1069–1078. [CrossRef] [PubMed]
Omar, M. A., and Shabana, A. A., 2001, “A Two-Dimensional Shear Deformable Beam for Large Rotation and Deformation Problems,” J. Sound Vib., 243(3), pp. 565–576. [CrossRef]
Patel, B. P., Ganapathi, M., and Touratier, M., 1999, “Nonlinear Free Flexural Vibrations/Post-Buckling Analysis of Laminated Orthotropic Beams/Columns on a Two Parameter Elastic Foundation,” Compos. Struct., 46(2), pp. 189–196. [CrossRef]
Ansari, M., Esmailzadeh, E., and Younesian, D., 2010, “Internal-External Resonance of Beams on Non-Linear Viscoelastic Foundation Traversed by Moving Load,” Nonlinear Dyn., 61(1–2), pp. 163–182. [CrossRef]
Ong, R. E., Herrell, S. D., Miga, M. I., and Galloway, R. L., 2008, “A Kidney Deformation Model for Use in Non-Rigid Registration During Image-Guided Surgery,” Proc. SPIE, 6918, pp. W9180–W9180. [CrossRef]
Jiang, S., Liu, S., and Feng, W., 2011, “PVA Hydrogel Properties for Biomedical Application,” J. Mech. Behav. Biomed., 4(7), pp. 1228–1233. [CrossRef]
van Gerwen, D. J., Dankelman, J., and van den Dobbelsteen, J. J., 2012, “Needle–Tissue Interaction Forces—A Survey of Experimental Data,” Med. Eng. Phys., 34(6), pp. 665–680. [CrossRef] [PubMed]
Yankelevsky, D. Z., Eisenberger, M., and Adin, M. A., 1989, “Analysis of Beams on Nonlinear Winkler Foundation,” Comput. Struct., 31(2), pp. 287–292. [CrossRef]
Mozer, P. C., Partin, A. W., and Stoianovici, D., 2009, “Robotic Image-Guided Needle Interventions of the Prostate,” Rev. Urol., 11(1), pp. 7–15. [PubMed]
Jiang, S., Li, P., Yu, Y., Liu, J., and Yang, Z., 2014, “Experimental Study of Needle-Tissue Interaction Forces: Effect of Needle Geometries, Insertion Methods and Tissue Characteristics,” ASME J. Biomech. Eng., 47(13), pp. 3344–3353. [CrossRef]


Grahic Jump Location
Fig. 1

The biomechanical model of the insertion system with forces acting on the whole needle as it interacts with the elastic medium, where linput is the applied needle base displacement working as the input of the system. The two dimensional model incorporates needle insertion forces (cutting force Ft, distribution forces fyfz), soft tissue Winkler foundation parameters k1k3 and flexible needle parameters (flexural rigidity EI, bevel angle α).

Grahic Jump Location
Fig. 2

The needle movement during two steering phases, dotted line represents the initial needle position (a) the insertion phase and (b) the adjustment phase

Grahic Jump Location
Fig. 3

Needle insertion using a robotic assisted needle insertion system to measure the insertion force

Grahic Jump Location
Fig. 4

Insertion forces profile into PVA phantom at six different insertion velocities. The experiment data and simulation input used in this study are acquired at a fixed velocity. The whole insertion process is continuous without interruption.

Grahic Jump Location
Fig. 5

The relationship between insertion force and insertion depth during puncture process. The averaged maximum puncture force (subtract friction force from total insertion force) is 0.1304 N.

Grahic Jump Location
Fig. 6

The horizontal force along the needle measured by a six-axis force sensor and the simulation result obtained by the mechanics model

Grahic Jump Location
Fig. 7

Depiction of the mechanics-based needle insertion model and methods to solve the reasonable numerical solution of needle deflection

Grahic Jump Location
Fig. 8

Five sequential movements with 10 mm depth of insertion phase. Each movement corresponds to a deflection curve obtained by the simulation result.

Grahic Jump Location
Fig. 9

Needle tip trajectories of needle–tissue interaction model with five movements. During each 10 mm insertion, only the tip position is marked out, so the horizontal axis contains the needle length.

Grahic Jump Location
Fig. 10

Five insertion phases and one adjustment phase of the whole needle beam. The base of the needle is forced to insert into the soft tissue. It can be seen that the displacement of needle tip is reduced comparing to Fig. 8.

Grahic Jump Location
Fig. 11

The needle deflection curves of different slenderness ratio λ in 200 mm insertion simulation. Other simulation parameters are the same as the needle insertion simulation shown in Fig. 8.

Grahic Jump Location
Fig. 12

The needle deflection curves of different friction coefficient η in 200 mm insertion simulation. Other simulation parameters are the same as the needle insertion simulation shown in Fig. 8.

Grahic Jump Location
Fig. 13

The needle deflection curves of different puncture force Ft, caused by different tissue stiffness in 30 mm insertion

Grahic Jump Location
Fig. 14

Robotic assisted needle insertion system and image acquisition system setup, used for needle deformation measurement

Grahic Jump Location
Fig. 15

Experimental data showing needle trajectories. In each 10 mm insertion, the needle tip position is measured to form a single curve. Then, the experimental data and simulation results for needle deflection are compared in one figure. The experimental image and dotted box indicate the needle insertion depth in the PVA phantom.

Grahic Jump Location
Fig. 16

Needle insertion simulation error between the simulated model and the experiment result

Grahic Jump Location
Fig. 17

Insertion force measured during the needle insertion into the PVA phantom. Except different bevel angles, other dimension parameters of the surgical needles are the same.

Grahic Jump Location
Fig. 18

Needle insertion experiment result, needle insertion simulation result, and simulation error of two surgical needles with different bevel angles. The insertion depth is set as 50 mm.




Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In