Fracture behaviour of human skin in deep needle insertion can be captured using validated cohesive zone finite-element method

Medical needles have shown an appreciable contribution to the development of novel medical devices and surgical technologies. A better understanding of needle-skin interactions can advance the design of medical needles, modern surgical robots, and haptic devices. This study employed finite element (FE) modelling to explore the effect of different mechanical and geometrical parameters on the needle's force-displacement relationship, the required force for the skin puncture, and generated mechanical stress around the cutting zone. To this end, we established a cohesive FE model, and identified its parameters by a three-stage parameter identification algorithm to closely replicate the experimental data of needle insertion into the human skin available in the literature. We showed that a bilinear cohesive model with initial stiffness of 5000 MPa/mm, failure traction of 2 MPa, and separation length of 1.6 mm can lead to a model that can closely replicate experimental results. The FE results indicated that while the coefficient of friction between the needle and skin substantially changes the needle reaction force, the insertion velocity does not have a noticeable effect on the reaction force. Regarding the geometrical parameters, needle cutting angle is the prominent factor in terms of stress fields generated in the skin tissue. However, the needle diameter is more influential on the needle reaction force. We also presented an energy study on the frictional dissipation, damage dissipation, and strain energy throughout the insertion process.