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

A method for predicting needle insertion deflection in soft tissue based on cutting force identification

The deflection modeling during the insertion of bevel-tipped flexible needles into soft tissues is crucial for robot-assisted flexible needle insertion into specific target locations within the human body during percutaneous biopsy surgery. This paper proposes a mechanical model based on cutting force identification to predict the deflection of flexible needles in soft tissues. Unlike other models, this method does not require measuring Young's modulus (E ) and Poisson's ratio (ν ) of tissues, which require complex hardware to obtain. In the model, the needle puncture process is discretized into a series of uniform-depth puncture steps. The needle is simplified as a cantilever beam supported by a series of virtual springs, and the influence of tissue stiffness on needle deformation is represented by the spring stiffness coefficient of the virtual spring. By theoretical modeling and experimental parameter identification of cutting force, the spring stiffness coefficients are obtained, thereby modeling the deflection of the needle. To verify the accuracy of the proposed model, the predicted model results were compared with the deflection of the puncture experiment in polyvinyl alcohol (PVA) gel samples, and the average maximum error range predicted by the model was between 0.606 ± 0.167 mm and 1.005 ± 0.174 mm, which showed that the model can successfully predict the deflection of the needle. This work will contribute to the design of automatic control strategies for needles.

Finite-element modelling of interactions of needle with tympanic membrane and middle ear

The tympanic membrane (TM) is one of the most common routes to access the middle ear and inner ear for the treatment of hearing and balance pathologies. Since the TM is a soft thin biological tissue with small dimensions, using needles seems to be among the most practical interventional approaches. In this study, we proposed a finite-element (FE) analysis of needle-TM interactions that combines a 3D model of the TM and other main middle-ear structures in gerbil, and a 2D model of needle insertion into the TM based on the cohesive zone method (CZM). The TM was modelled using a 1st-order Ogden hyperelastic material and its properties were obtained by fitting to the experimental force-displacement plots of large deformation in the TM under needle indentation. The cohesive parameters were also acquired by calibrating the puncture force against the experimental data of needle insertion into the TM. These FE models were then used to obtain the deformation behaviour of the TM and other middle-ear structures due to the insertion force applied at different locations on the TM. Moreover, we investigated the effect of the TM thickness, the geometry of the needle (i.e., diameter and tip angle), and needle material on the insertion of needles into the TM. We also studied the penetration success of deformable needles.

Experimental Study of Needle Insertion into Gerbil Tympanic Membrane

The perforation characteristics and fracture-related mechanical properties of the tympanic membrane (TM) greatly affect surgical procedures like myringotomy and tympanostomy performed on the middle ear. We analyzed the most important features of the gerbil TM perforation using an experimental approach that was based on force measurement during a 2-cycle needle insertion/extraction process. Fracture energy, friction energy, strain energy, and hysteresis loss were taken into consideration for the analysis of the different stages of needle insertion and extraction. The results demonstrated that (1) although the TM shows viscoelastic behavior, the contribution of hysteresis loss was negligible compared to other irreversible dissipated energy components (i.e., fracture energy and friction energy). (2) The TM puncture force did not substantially change during the first hours after animal death, but interestingly, it increased after 1 week due to the drying effects of soft tissue. (3) The needle geometry affected the crack length and the most important features of the force-displacement plot for the needle insertion process (puncture force, puncture displacement, and jump-in force) increased with increasing needle diameter, whereas the insertion velocity only changed the puncture and jump-in forces (both increased with increasing insertion velocity) and did not have a noticeable effect on the puncture displacement. (4) The fracture toughness of the gerbil TM was almost independent of the needle geometry and was found to be around 0.33 ± 0.10 kJ/m2.

Optimizing autoinjector devices using physics-based simulations and Gaussian processes

Autoinjectors are becoming a primary drug delivery option to the subcutaneous space. These devices need to work robustly and autonomously to maximize drug bio-availability. However, current designs ignore the coupling between autoinjector dynamics and tissue biomechanics. Here we present a Bayesian framework for optimization of autoinjector devices that can account for the coupled autoinjector-tissue biomechanics and uncertainty in tissue mechanical behavior. The framework relies on replacing the high fidelity model of tissue insertion with a Gaussian process (GP). The GP model is accurate yet computationally affordable, enabling a thorough sensitivity analysis that identified tissue properties, which are not part of the autoinjector design space, as important variables for the injection process. Higher fracture toughness decreases the crack depth, while tissue shear modulus has the opposite effect. The sensitivity analysis also shows that drug viscosity and spring force, which are part of the design space, affect the location and timing of drug delivery. Low viscosity could lead to premature delivery, but can be prevented with smaller spring forces, while higher viscosity could prevent premature delivery while demanding larger spring forces and increasing the time of injection. Increasing the spring force guarantees penetration to the desired depth, but it can result in undesirably high accelerations. The Bayesian optimization framework tackles the challenge of designing devices with performance metrics coupled to uncertain tissue properties. This work is important for the design of other medical devices for which optimization in the presence of material behavior uncertainty is needed.

Modelling biological puncture: a mathematical framework for determining the energetics and scaling

Biological puncture systems use a diversity of morphological tools (stingers, teeth, spines etc.) to penetrate target tissues for a variety of functions (prey capture, defence, reproduction). These systems are united by a set of underlying physical rules which dictate their mechanics. While previous studies have illustrated form-function relationships in individual systems, these underlying rules have not been formalized. We present a mathematical model for biological puncture events based on energy balance that allows for the derivation of analytical scaling relations between energy expenditure and shape, size and material response. The model identifies three necessary energy contributions during puncture: fracture creation, elastic deformation of the material and overcoming friction during penetration. The theoretical predictions are verified using finite-element analyses and experimental tests. Comparison between different scaling relationships leads to a ratio of released fracture energy and deformation energy contributions acting as a measure of puncture efficiency for a system that incorporates both tool shape and material response. The model represents a framework for exploring the diversity of biological puncture systems in a rigorous fashion and allows future work to examine how fundamental physical laws influence the evolution of these systems.