An adaptive finite element model for steerable needles.
Ontology highlight
ABSTRACT:
SUBMITTER: Terzano M
PROVIDER: S-EPMC7502456 | biostudies-literature | 2020 Oct
REPOSITORIES: biostudies-literature
Publications
An adaptive finite element model for steerable needles.
Terzano Michele M Dini Daniele D Rodriguez Y Baena Ferdinando F Spagnoli Andrea A Oldfield Matthew M
Biomechanics and modeling in mechanobiology 20200309 5
Penetration of a flexible and steerable needle into a soft target material is a complex problem to be modelled, involving several mechanical challenges. In the present paper, an adaptive finite element algorithm is developed to simulate the penetration of a steerable needle in brain-like gelatine material, where the penetration path is not predetermined. The geometry of the needle tip induces asymmetric tractions along the tool-substrate frictional interfaces, generating a bending action on the ...[more]
Similar Datasets
Project description:We consider the design of an effective and reliable adaptive finite element method (AFEM) for the nonlinear Poisson-Boltzmann equation (PBE). We first examine the two-term regularization technique for the continuous problem recently proposed by Chen, Holst, and Xu based on the removal of the singular electrostatic potential inside biomolecules; this technique made possible the development of the first complete solution and approximation theory for the Poisson-Boltzmann equation, the first provably convergent discretization, and also allowed for the development of a provably convergent AFEM. However, in practical implementation, this two-term regularization exhibits numerical instability. Therefore, we examine a variation of this regularization technique which can be shown to be less susceptible to such instability. We establish a priori estimates and other basic results for the continuous regularized problem, as well as for Galerkin finite element approximations. We show that the new approach produces regularized continuous and discrete problems with the same mathematical advantages of the original regularization. We then design an AFEM scheme for the new regularized problem, and show that the resulting AFEM scheme is accurate and reliable, by proving a contraction result for the error. This result, which is one of the first results of this type for nonlinear elliptic problems, is based on using continuous and discrete a priori L(∞) estimates to establish quasi-orthogonality. To provide a high-quality geometric model as input to the AFEM algorithm, we also describe a class of feature-preserving adaptive mesh generation algorithms designed specifically for constructing meshes of biomolecular structures, based on the intrinsic local structure tensor of the molecular surface. All of the algorithms described in the article are implemented in the Finite Element Toolkit (FETK), developed and maintained at UCSD. The stability advantages of the new regularization scheme are demonstrated with FETK through comparisons with the original regularization approach for a model problem. The convergence and accuracy of the overall AFEM algorithm is also illustrated by numerical approximation of electrostatic solvation energy for an insulin protein.
| S-EPMC3178275 | biostudies-literature
Project description:BACKGROUND: Biological mass transport processes determine the behavior and function of cells, regulate interactions between synthetic agents and recipient targets, and are key elements in the design and use of biosensors. Accurately predicting the outcomes of such processes is crucial to both enhancing our understanding of how these systems function, enabling the design of effective strategies to control their function, and verifying that engineered solutions perform according to plan. METHODS: A Galerkin-based finite element model was developed and implemented to solve a system of two coupled partial differential equations governing biomolecule transport and reaction in live cells. The simulator was coupled, in the framework of an inverse modeling strategy, with an optimization algorithm and an experimental time series, obtained by the Fluorescence Recovery after Photobleaching (FRAP) technique, to estimate biomolecule mass transport and reaction rate parameters. In the inverse algorithm, an adaptive method was implemented to calculate sensitivity matrix. A multi-criteria termination rule was developed to stop the inverse code at the solution. The applicability of the model was illustrated by simulating the mobility and binding of GFP-tagged glucocorticoid receptor in the nucleoplasm of mouse adenocarcinoma. RESULTS: The numerical simulator shows excellent agreement with the analytic solutions and experimental FRAP data. Detailed residual analysis indicates that residuals have zero mean and constant variance and are normally distributed and uncorrelated. Therefore, the necessary and sufficient criteria for least square parameter optimization, which was used in this study, were met. CONCLUSION: The developed strategy is an efficient approach to extract as much physiochemical information from the FRAP protocol as possible. Well-posedness analysis of the inverse problem, however, indicates that the FRAP protocol provides insufficient information for unique simultaneous estimation of diffusion coefficient and binding rate parameters. Care should be exercised in drawing inferences, from FRAP data, regarding concentrations of free and bound proteins, average binding and diffusion times, and protein mobility unless they are confirmed by long-range Markov Chain-Monte Carlo (MCMC) methods and experimental observations.
| S-EPMC1940256 | biostudies-literature
Project description:Nasal tip mechanical stability is important for functional and cosmetic nasal airway surgery. Palpation of the nasal tip provides information on tip strength to the surgeon, though it is a purely subjective assessment. Providing a means to simulate nasal tip deformation with a validated model can offer a more objective approach in understanding the mechanics and nuances of the nasal tip support and eventual nasal mechanics as a whole. Herein we present validation of a finite element (FE) model of the nose using physical measurements recorded using an ABS plastic-silicone nasal phantom. Three-dimensional photogrammetry was used to capture the geometry of the phantom at rest and while under steady state load. The silicone used to make the phantom was mechanically tested and characterized using a linear elastic constitutive model. Surface point clouds of the silicone and FE model were compared for both the loaded and unloaded state. The average Hausdorff distance between actual measurements and FE simulations across the nose were 0.39 ± 1.04 mm and deviated up to 2 mm at the outermost boundaries of the model. FE simulation and measurements were in near complete agreement in the immediate vicinity of the nasal tip with millimeter accuracy. We have demonstrated validation of a two-component nasal FE model, which could be used to model more complex modes of deformation where direct measurement may be challenging. This is the first step in developing a nasal model to simulate nasal mechanics and ultimately the interaction between geometry and airflow.
| S-EPMC5332342 | biostudies-literature
Project description:The functional biomechanics of the lumbar spine have been better understood by finite element method (FEM) simulations. However, there are still areas where the behavior of soft tissues can be better modeled or described in a different way. The purpose of this research is to develop and validate a lumbar spine section intended for biomechanical research. A FE model of the 50th percentile adult male (AM) Total Human Model for Safety (THUMS) v6.1 was used to implement the modifications. The main modifications were to apply orthotropic material properties and nonlinear stress-strain behavior for ligaments, hyperelastic material properties for annulus fibrosus and nucleus pulposus, and the specific content of collagenous fibers in the annulus fibrosus ground substance. Additionally, a separation of the nucleus pulposus from surrounding bones and tissues was implemented. The FE model was subjected to different loading modes, in which intervertebral rotations and disc pressures were calculated. Loading modes contained different forces and moments acting on the lumbar section: axial forces (compression and tension), shear forces, pure moments, and combined loading modes of axial forces and pure moments. The obtained ranges of motion from the modified numerical model agreed with experimental data for all loading modes. Moreover, intradiscal pressure validation for the modified model presented a good agreement with the data available from the literature. This study demonstrated the modifications of the THUMS v6.1 model and validated the obtained numerical results with existing literature in the sub-injurious range. By applying the proposed changes, it is possible to better model the behavior of the human lumbar section under various loads and moments.
| S-EPMC10424670 | biostudies-literature
Project description:Shape plays an important role in determining the biomechanical response of a structure. Specimen-specific finite element (FE) models have been developed to capture the details of the shape of biological structures and predict their biomechanics. Shape, however, can vary considerably across individuals or change due to aging or disease, and analysis of the sensitivity of specimen-specific models to these variations has proven challenging. An alternative to specimen-specific representation has been to develop generic models with simplified geometries whose shape is relatively easy to parameterize, and can therefore be readily used in sensitivity studies. Despite many successful applications, generic models are limited in that they cannot make predictions for individual specimens. We propose that it is possible to harness the detail available in specimen-specific models while leveraging the power of the parameterization techniques common in generic models. In this work we show that this can be accomplished by using morphing techniques to parameterize the geometry of specimen-specific FE models such that the model shape can be varied in a controlled and systematic way suitable for sensitivity analysis. We demonstrate three morphing techniques by using them on a model of the load-bearing tissues of the posterior pole of the eye. We show that using relatively straightforward procedures these morphing techniques can be combined, which allows the study of factor interactions. Finally, we illustrate that the techniques can be used in other systems by applying them to morph a femur. Morphing techniques provide an exciting new possibility for the analysis of the biomechanical role of shape, independently or in interaction with loading and material properties.
| S-EPMC2813388 | biostudies-literature
Project description:Isolating the role(s) of microstructural pathological features in affecting diastolic filling is important in developing targeted treatments for heart diseases. We developed a microstructure-based constitutive model of the myocardium and implemented it in an efficient open-source finite element modeling framework to simulate passive inflation of the left ventricle (LV) in a representative 3D geometry based on experimentally measured muscle fiber architecture. The constitutive model was calibrated using previous tissue-level biaxial mechanical test data derived from the canine heart and validated with independent sets of measurements made at both the isolated constituent and organ level. Using the validated model, we investigated the load taken up by each tissue constituent and their effects on LV passive inflation. The model predicts that the LV compliance is sensitive to the collagen ultrastructure, specifically, the collagen fiber azimuthal angle with respect to the local muscle fiber direction and its waviness. The model also predicts that most of the load in the sub-epicardial and sub-endocardial regions is taken up, respectively, by the muscle fibers and collagen fiber network. This result suggests that normalizing LV passive stiffness by altering the collagen fiber network and myocyte stiffness is most effective when applied to the sub-endocardial and sub-epicardial regions, respectively. This finding may have implication for the development of new pharmaceutical treatments targeting individual cardiac tissue constituents to normalize LV filling function in heart diseases. STATEMENT OF SIGNIFICANCE: Current constitutive models describing the tissue mechanical behavior of the myocardium are largely phenomenological. While able to represent the bulk tissue mechanical behavior, these models cannot distinguish the contribution of the tissue constituents and their ultrastructure to heart function. Although microstructure-based constitutive models can be used to isolate the role of tissue ultrastructure, they have not been implemented in a computational framework that can accommodate realistic 3D organ geometry. The present study addresses these issues by developing and validating a microstructure-based computational modeling framework, which is used to investigate the role of tissue constituents and their ultrastructure in affecting heart function.
| S-EPMC6677579 | biostudies-literature
Project description:Assessment of heart function in zebrafish larvae using electrocardiography (ECG) is a potentially useful tool in developing cardiac treatments and the assessment of drug therapies. In order to better understand how a measured ECG waveform is related to the structure of the heart, its position within the larva and the position of the electrodes, a 3D model of a 3 days post fertilisation (dpf) larval zebrafish was developed to simulate cardiac electrical activity and investigate the voltage distribution throughout the body. The geometry consisted of two main components; the zebrafish body was modelled as a homogeneous volume, while the heart was split into five distinct regions (sinoatrial region, atrial wall, atrioventricular band, ventricular wall and heart chambers). Similarly, the electrical model consisted of two parts with the body described by Laplace's equation and the heart using a bidomain ionic model based upon the Fitzhugh-Nagumo equations. Each region of the heart was differentiated by action potential (AP) parameters and activation wave conduction velocities, which were fitted and scaled based on previously published experimental results. ECG measurements in vivo at different electrode recording positions were then compared to the model results. The model was able to simulate action potentials, wave propagation and all the major features (P wave, R wave, T wave) of the ECG, as well as polarity of the peaks observed at each position. This model was based upon our current understanding of the structure of the normal zebrafish larval heart. Further development would enable us to incorporate features associated with the diseased heart and hence assist in the interpretation of larval zebrafish ECGs in these conditions.
| S-EPMC5100939 | biostudies-literature
Project description:Alar rim retraction is the most common unintended consequence of tissue remodeling that results from overresection of the cephalic lateral crural cartilage; however, the complex tissue remodeling process that produces this shape change is not well understood.To simulate how resection of cephalic trim alters the stress distribution within the human nose in response to tip depression (palpation) and to simulate the internal forces generated after cephalic trim that may lead to alar rim retraction cephalically and upward rotation of the nasal tip.A multicomponent finite element model was derived from maxillofacial computed tomography with 1-mm axial resolution. The 3-dimensional editing function in the medical imaging software was used to trim the cephalic portion of the lower lateral cartilage to emulate that performed in typical rhinoplasty. Three models were created: a control, a conservative trim, and an aggressive trim. Each simulated model was imported to a software program that performs mechanical simulations, and material properties were assigned. First, nasal tip depression (palpation) was simulated, and the resulting stress distribution was calculated for each model. Second, long-term tissue migration was simulated on conservative and aggressive trim models by placing normal and shear force vectors along the caudal and cephalic borders of the tissue defect.The von Mises stress distribution created by a 5-mm tip depression revealed consistent findings among all 3 simulations, with regions of high stress being concentrated to the medial portion of the intermediate crus and the caudal septum. Nasal tip reaction force marginally decreased as more lower lateral cartilage tissue was resected. Conservative and aggressive cephalic trim models produced some degree of alar rim retraction and tip rotation, which increased with the magnitude of the force applied to the region of the tissue defect.Cephalic trim was performed on a computerized composite model of the human nose to simulate conservative and aggressive trims. Internal forces were applied to each model to emulate the tissue migration that results from decades of wound healing. Our simulations reveal that the degree of tip rotation and alar rim retraction is dependent on the amount of cartilage that was resected owing to cephalic trim. Tip reaction force is marginally reduced with increasing tissue volume resection.NA.
| S-EPMC5847269 | biostudies-literature
Project description:The lack of a universal simulation method for triboelectric nanogenerator (TENG) makes the device design and optimization difficult before experiment, which protracts the research and development process and hinders the landing of practical TENG applications. The existing electrostatic induction models for TENGs have limitations in simulating TENGs with complex geometries and their dynamic behaviors under practical movements due to the topology change issues. Here, a dynamic finite element method (FEM) model is proposed. The introduction of air buffer layers and the moving mesh method eliminates the topology change issues during practical movement and allows simulation of dynamic and time-varying behaviors of TENGs with complex 2D/3D geometries. Systematic investigations are carried out to optimize the air buffer thickness and mesh densities, and the optimized results show excellent consistency with the experimental data and results based on other existing methods. It also shows that a 3D disk-type rotating TENG can be simulated using the model, clearly demonstrating the capability and superiority of the dynamic FEM model. Moreover, the dynamic FEM model is used to optimize the shape of the tribo-material, which is used as a preliminary example to demonstrate the possibility of designing a TENG-based sensor.
| S-EPMC7506595 | biostudies-literature
Project description:BackgroundCurrently, there is uncertainty regarding the long-term outcome of medial patellofemoral ligament reconstructions (MPFLr). Our objectives were: (1) to develop a parametric model of the patellofemoral joint (PFJ) enabling us to simulate different surgical techniques for MPFLr; (2) to determine the negative effects on the PFJ associated with each technique, which could be related to long-term deterioration of the PFJ.MethodsA finite element model of the PFJ was created based on CT data from 24 knees with chronic lateral patellar instability. Patella contact pressure and maximum MPFL-graft stress at five angles of knee flexion (0, 30, 60, 90 and 120°) were analysed in three types of MPFLr: anatomic, non-anatomic with physiometric behaviour, and non-anatomic with non-physiometric behaviour.ResultsAn increase in patella contact pressure was observed at 0 and 30° of knee flexion after both anatomic and non-anatomic MPFLr with physiometric behaviour. In both reconstructions, the ligament was tense between 0 and 30° of knee flexion, but at 60, 90 and 120°, it had no tension. In the third reconstruction, the behaviour was completely the opposite.ConclusionA parametric model of the PFJ enables us to evaluate different types of MPFLr throughout the full range of motion of the knee, regarding the effect on the patellofemoral contact pressure, as well as the kinematic behaviour of the MPFL-graft and the maximum MPFL-graft stress.
| S-EPMC6611858 | biostudies-literature