Nonlinear Elastic Material Property Estimation of Lower Extremity Residual Limb Tissues

The interface stresses between the residual limb and prosthetic socket have been studied to investigate prosthetic fit. Finite-element models of the residual limb-prosthetic socket interface facilitate investigation of the mechanical interface and may serve as a potential tool for future prosthetic socket design. However, the success of such residual limb models to date has been limited, in large part due to inadequate material formulations used to approximate the mechanical behavior of residual limb soft tissues. Nonlinear finite-element analysis was used to simulate force-displacement data obtained during in vivo rate-controlled (1, 5, and 10 mm/s) cyclic indentation of the residual limb soft tissues of seven individuals with transtibial amputation. The finite-element models facilitated determination of an appropriate set of nonlinear elastic material coefficients for bulk soft tissue at discrete clinically relevant test locations. Axisymmetric finite-element models of the residual limb bulk soft tissue in the vicinity of the test location, the socket wall and the indentor tip were developed incorporating contact analysis, large displacement, and large strain, and the James-Green-Simpson nonlinear elastic material formulation. Model dimensions were based on medical imaging studies of the residual limbs. The material coefficients were selected such that the normalized sum of square error (NSSE) between the experimental and finite-element model indentor tip reaction force was minimized. A total of 95% of the experimental data were simulated using the James-Green-Simpson material formulation with an NSSE less than 5%. The respective James-Green-Simpson material coefficients varied with subject, test location, and indentation rate. Therefore, these coefficients cannot be readily extrapolated to other sites or individuals, or to the same site and individual some time after testing

Controlling the Error on Target Motion through Real-time Mesh Adaptation: Applications to Deep Brain Stimulation

We present an error-controlled mesh refinement procedure for needle insertion simulation and apply it to the simulation of electrode implantation for deep brain stimulation, including brain shift. Our approach enables to control the error in the computation of the displacement and stress fields around the needle tip and needle shaft by suitably refining the mesh, whilst maintaining a coarser mesh in other parts of the domain. We demonstrate through academic and practical examples that our approach increases the accuracy of the displacement and stress fields around the needle without increasing the computational expense. This enables real-time simulations. The proposed methodology has direct implications to increase the accuracy and control the computational expense of the simulation of percutaneous procedures such as biopsy, brachytherapy, regional anesthesia, or cryotherapy and can be essential to the development of robotic guidance.Comment: 21 pages, 14 figure

Finite element model creation and stability considerations of complex biological articulation : the human wrist joint

The finite element method has been used with considerable success to simulate the behaviour of various joints such as the hip, knee and shoulder. It has had less impact on more complicated joints such as the wrist and the ankle. Previously published finite element studies on these multi bone joints have needed to introduce un-physiological boundary conditions in order to establish numerical convergence of the model simulation. That is necessary since the stabilising soft tissue mechanism of these joints is usually too elaborate in order to be fully included both anatomically and with regards to material properties. This paper looks at the methodology of creating a finite element model of such a joint focussing on the wrist and the effects additional constraining has on the solution of the model. The study shows that by investigating the effects each of the constraints, a better understanding on the nature of the stabilizing mechanisms of these joints can be achieved

A three-dimensional finite element model of maximal grip loading in the human wrist

The aim of this work was to create an anatomically accurate three-dimensional finite element model of the wrist, applying subject-specific loading and quantifying the internal load transfer through the joint during maximal grip. For three subjects, representing the anatomical variation at the wrist, loading on each digit was measured during a maximal grip strength test with simultaneous motion capture. The internal metacarpophalangeal joint load was calculated using a biomechanical model. High-resolution magnetic resonance scans were acquired to quantify bone geometry. Finite element analysis was performed, with ligaments and tendons added, to calculate the internal load distribution. It was found that for the maximal grip the thumb carried the highest load, an average of 72.2 ¡ 20.1 N in the neutral position. Results from the finite element model suggested that the highest regions of stress were located at the radial aspect of the carpus. Most of the load was transmitted through the radius, 87.5 per cent, as opposed to 12.5 per cent through the ulna with the wrist in a neutral position. A fully three-dimensional finite element analysis of the wrist using subject-specific anatomy and loading conditions was performed. The study emphasizes the importance of modelling a large ensemble of subjects in order to capture the spectrum of the load transfer through the wrist due to anatomical variation

Full modelling of high-intensity focused ultrasound and thermal heating in the kidney using realistic patient models

Objective: High-intensity focused ultrasound (HIFU) therapy can be used for non-invasive treatment of kidney (renal) cancer, but the clinical outcomes have been variable. In this study, the efficacy of renal HIFU therapy was studied using nonlinear acoustic and thermal simulations in three patients. Methods: The acoustic simulations were conducted with and without refraction in order to investigate its effect on the shape, size and pressure distribution at the focus. The values for the attenuation, sound speed, perfusion and thermal conductivity of the kidney were varied over the reported ranges to determine the effect of variability on heating. Furthermore, the phase aberration was studied in order to quantify the underlying phase shifts using a second order polynomial function. Results: The ultrasound field intensity was found to drop on average 11.1 dB with refraction and 6.4 dB without refraction. Reflection at tissue interfaces was found to result in a loss less than 0.1 dB. Focal point splitting due to refraction significantly reduced the heating efficacy. Perfusion did not have a large effect on heating during short sonication durations. Small changes in temperature were seen with varying attenuation and thermal conductivity, but no visible changes were present with sound speed variations. The aberration study revealed an underlying trend in the spatial distribution of the phase shifts. Conclusion: The results show that the efficacy of HIFU therapy in the kidney could be improved with aberration correction. Significance: A method is proposed by which patient specific pre-treatment calculations could be used to overcome the aberration and therefore make ultrasound treatment possible.Comment: Journal paper, IEEE Transactions on Biomedical Engineering (2018

Intimal and medial contributions to the hydraulic resistance of the arterial wall at different pressures: a combined computational and experimental study

The hydraulic resistances of the intima and media determine water flux and the advection of macromolecules into and across the arterial wall. Despite several experimental and computational studies, these transport processes and their dependence on transmural pressure remain incompletely understood. Here, we use a combination of experimental and computational methods to ascertain how the hydraulic permeability of the rat abdominal aorta depends on these two layers and how it is affected by structural rearrangement of the media under pressure. Ex vivo experiments determined the conductance of the whole wall, the thickness of the media and the geometry of medial smooth muscle cells (SMCs) and extracellular matrix (ECM). Numerical methods were used to compute water flux through the media. Intimal values were obtained by subtraction. A mechanism was identified that modulates pressure-induced changes in medial transport properties: compaction of the ECM leading to spatial reorganization of SMCs. This is summarized in an empirical constitutive law for permeability and volumetric strain. It led to the physiologically interesting observation that, as a consequence of the changes in medial microstructure, the relative contributions of the intima and media to the hydraulic resistance of the wall depend on the applied pressure; medial resistance dominated at pressures above approximately 93 mmHg in this vessel

Implementing vertex dynamics models of cell populations in biology within a consistent computational framework

The dynamic behaviour of epithelial cell sheets plays a central role during development, growth, disease and wound healing. These processes occur as a result of cell adhesion, migration, division, differentiation and death, and involve multiple processes acting at the cellular and molecular level. Computational models offer a useful means by which to investigate and test hypotheses about these processes, and have played a key role in the study of cell–cell interactions. However, the necessarily complex nature of such models means that it is difficult to make accurate comparison between different models, since it is often impossible to distinguish between differences in behaviour that are due to the underlying model assumptions, and those due to differences in the in silico implementation of the model. In this work, an approach is described for the implementation of vertex dynamics models, a discrete approach that represents each cell by a polygon (or polyhedron) whose vertices may move in response to forces. The implementation is undertaken in a consistent manner within a single open source computational framework, Chaste, which comprises fully tested, industrial-grade software that has been developed using an agile approach. This framework allows one to easily change assumptions regarding force generation and cell rearrangement processes within these models. The versatility and generality of this framework is illustrated using a number of biological examples. In each case we provide full details of all technical aspects of our model implementations, and in some cases provide extensions to make the models more generally applicable

Computational modelling of void growth in swelled hydrogels

The nature and the large notable distinguishing features of polymeric gels explain their pervasive use as biomaterials in both regenerative medicine and tissue engineering. With regard to their biocompatibility, their ability to withstand large deformation and their significant capacity of solvent absorption, these biomaterials are often selected owing to their versatile mechanical properties and especially the closeness to soft biological tissues, amongst others. A finite-strain theory for the study of the overall behaviour of a porous polymeric gel where microvoids are present is presented. The swollen polymeric gel is modeled as a two-component body composed of two incompressible components, namely, an elastic porous polymer imbibed with a solvant. The chemical equilibrium is assumed to be preponderate at the interface between the porous polymer and the environment where the chemical potential of the solvent is fixed. The initially dry porous polymer undergoes large deformation induced by absorption of a solvent from the environment and mechanical loading. In this paper an attempt is made towards obtaining an estimation of the macroscopic responses of the swollen porous polymer to prescribed proportional loadings. To this end, a two-level representation of the material at hand for which the Representative Volume Element (RVE) imbibed with a solvent is a simple axisymmetric cylinder composed of a homogeneous matrix surrounding a spherical void, is considered. The computational study addresses the situation where the RVE is subjected to prescribed axial and lateral overall stresses under conditions of constant overall stress triaxiality. For fixed values of the Flory-Huggins parameter and the nominal concentration of the solvent, the overall stress-strain behaviour of the RVE model, the influence of the initial porosity, and the prescribed stress triaxiality ratio have been outlined