A convenient scheme for coupling a finite element curvilinear mesh to a finite element voxel mesh: application to the heart

Background  

In some cases, it may be necessary to combine distinct finite element meshes into a single system. The present work describes a scheme for coupling a finite element mesh, which may have curvilinear elements, to a voxel based finite element mesh.     

A novel adaptive algorithm for 3D finite element analysis to model extracortical bone growth

Extracortical bone growth with osseointegration of bone onto the shaft of massive bone tumour implants is an important clinical outcome for long-term implant survival. A new computational algorithm combining geometrical shape changes and bone adaptation in 3D Finite Element simulations has been developed, using a soft tissue envelope mesh, a novel concept of osteoconnectivity, and bone remodelling theory. The effects of varying the initial tissue density, spatial influence function and time step were investigated. The methodology demonstrated good correspondence to radiological results for a segmental prosthesis.     

3D finite element model of biofilm detachment using real biofilm structures from CLSM data

In this work, a three‐dimensional model of fluid–structure interactions (FSI) in biofilm systems is developed in order to simulate biofilm detachment as a result of mechanical processes. Therein, fluid flow past the biofilm surface results in a mechanical load on the structure which in turn causes internal stresses in the biofilm matrix. When the strength of the matrix is exceeded parts of the structure are detached. The model is used to investigate the influence of several parameters related to the mechanical strength of the biofilm matrix, Young's modulus, Reynolds number, and biofilm structure on biofilm detachment. Variations in biofilm strength and flow conditions significantly influence the simulation outcome. With respect to structural properties the model is widely independent from a change of Young's modulus. A further result of this work indicates that the change of biofilm structure due to growth or other processes will significantly change the stress distribution in the biofilm and thereby the detachment rate. An increase of the mechanical load by increasing fluid flow results in a flat surface of the remaining biofilm structure. It is concluded that the change of structure during biofilm development is the key determinant in terms of the detachment behavior. Biotechnol. Bioeng. 2009;103: 177–186. © 2008 Wiley Periodicals, Inc.     

A penetration-based finite element method for hyperelastic 3D biphasic tissues in contact. Part II: finite element simulations

The penetration method allows for the efficient finite element simulation of contact between soft hydrated biphasic tissues in diarthrodial joints. Efficiency of the method is achieved by separating the intrinsically nonlinear contact problem into a pair of linked biphasic finite element analyses, in which an approximate, spatially and temporally varying contact traction is applied to each of the contacting tissues. In Part I of this study, we extended the penetration method to contact involving nonlinear biphasic tissue layers, and demonstrated how to derive the approximate contact traction boundary conditions. The traction derivation involves time and space dependent natural boundary conditions, and requires special numerical treatment. This paper (Part II) describes how we obtain an efficient nonlinear finite element procedure to solve for the biphasic response of the individual contacting layers. In particular, alternate linearization of the nonlinear weak form, as well as both velocity-pressure, v-p, and displacement-pressure, u-p, mixed formulations are considered. We conclude that the u-p approach, with linearization of both the material law and the deformation gradients, performs best for the problem at hand. The nonlinear biphasic contact solution will be demonstrated for the motion of the glenohumeral joint of the human shoulder joint.     

3D finite element modeling of pelvic organ prolapse

Objectives: The purpose of this study is to develop a validated 3D finite element model of the pelvic floor system which can offer insights into the mechanics of anterior vaginal wall prolapse and have the ability to assess biomedical device treatment methods. The finite element results should accurately mimic the clinical findings of prolapse due to intra-abdominal pressure (IAP) and soft tissues impairment conditions. Methods: A 3D model of pelvic system was created in Creo Parametric 2.0 based on MRI Images, which included uterus, cervix, vagina, cardinal ligaments, uterosacral ligaments, and a simplified levator plate and rectum. The geometrical model was imported into ANSYS Workbench 14.5. Mechanical properties of soft tissues were based on experimental data of tensile test results from current literature. Studies were conducted for IAP loadings on the vaginal wall and uterus, increasing from lowest to extreme values. Results: Anterior vaginal wall collapse occurred at an IAP value corresponding to maximal valsalva and showed similar collapsed shape as clinical findings. Prolapse conditions exhibited high sensitivity to vaginal wall stiffness, whereas healthy tissues was found to support the vagina against prolapse. Ligament impairment was found to have only a secondary effect on prolapse.     

3D finite element analysis of porous Ti-based alloy prostheses

In this paper, novel designs of porous acetabular cups are created and tested with 3D finite element analysis (FEA). The aim is to develop a porous acetabular cup with low effective radial stiffness of the structure, which will be near to the architectural and mechanical behavior of the natural bone. For the realization of this research, a 3D-scanner technology was used for obtaining a 3D-CAD model of the pelvis bone, a 3D-CAD software for creating a porous acetabular cup, and a 3D-FEA software for virtual testing of a novel design of the porous acetabular cup. The results obtained from this research reveal that a porous acetabular cup from Ti-based alloys with 60 ± 5% porosity has the mechanical behavior and effective radial stiffness (Young’s modulus in radial direction) that meet and exceed the required properties of the natural bone. The virtual testing with 3D-FEA of a novel design with porous structure during the very early stage of the design and the development of orthopedic implants, enables obtaining a new or improved biomedical implant for a relatively short time and reduced price.     

Failure locus of the anterior cruciate ligament: 3D finite element analysis

Anterior cruciate ligament (ACL) disruption is a common injury that is detrimental to an athlete's quality of life. Determining the mechanisms that cause ACL injury is important in order to develop proper interventions. A failure locus defined as various combinations of loadings and movements, internal/external rotation of femur and valgus and varus moments at a 25^o knee flexion angle leading to ACL failure was obtained. The results indicated that varus and valgus movements were more dominant to the ACL injury than femoral rotation. Also, Von Mises stress in the lateral tibial cartilage during the valgus ACL injury mechanism was 83% greater than that of the medial cartilage during the varus mechanism of ACL injury. The results of this study could be used to develop training programmes focused on the avoidance of the described combination of movements which may lead to ACL injury.     

3D finite element model of meniscectomy: changes in joint contact behavior

The goal of this study is to quantify changes in knee joint contact behavior following varying degrees of the medial partial meniscectomy. A previously validated 3D finite element model was used to simulate 11 different meniscectomies. The accompanying changes in the contact pressure on the superior surface of the menisci and tibial plateau were quantified as was the axial strain in the menisci and articular cartilage. The percentage of medial meniscus removed was linearly correlated with maximum contact pressure, mean contact pressure, and contact area. The lateral hemi-joint was minimally affected by the simulated medial meniscectomies. The location of maximum strain and location of maximum contact pressure did not change with varying degrees of partial medial meniscectomy. When 60% of the medial meniscus was removed, contact pressures increased 65% on the remaining medial meniscus and 55% on the medial tibial plateau. These data will be helpful for assessing potential complications with the surgical treatment of meniscal tears. Additionally, these data provide insight into the role of mechanical loading in the etiology of post-meniscectomy osteoarthritis.     

An augmented Lagrangian finite element formulation for 3D contact of biphasic tissues

Biphasic contact analysis is essential to obtain a complete understanding of soft tissue biomechanics, and the importance of physiological structure on the joint biomechanics has long been recognised; however, up to date, there are no successful developments of biphasic finite element contact analysis for three-dimensional (3D) geometries of physiological joints. The aim of this study was to develop a finite element formulation for biphasic contact of 3D physiological joints. The augmented Lagrangian method was used to enforce the continuity of contact traction and fluid pressure across the contact interface. The biphasic contact method was implemented in the commercial software COMSOL Multiphysics 4.2^® (COMSOL, Inc., Burlington, MA). The accuracy of the implementation was verified using 3D biphasic contact problems, including indentation with a flat-ended indenter and contact of glenohumeral cartilage layers. The ability of the method to model multibody biphasic contact of physiological joints was proved by a 3D knee model. The 3D biphasic finite element contact method developed in this study can be used to study the biphasic behaviours of the physiological joints.     

Failure locus of the anterior cruciate ligament: 3D finite element analysis

Anterior cruciate ligament (ACL) disruption is a common injury that is detrimental to an athlete's quality of life. Determining the mechanisms that cause ACL injury is important in order to develop proper interventions. A failure locus defined as various combinations of loadings and movements, internal/external rotation of femur and valgus and varus moments at a 25(o) knee flexion angle leading to ACL failure was obtained. The results indicated that varus and valgus movements were more dominant to the ACL injury than femoral rotation. Also, Von Mises stress in the lateral tibial cartilage during the valgus ACL injury mechanism was 83% greater than that of the medial cartilage during the varus mechanism of ACL injury. The results of this study could be used to develop training programmes focused on the avoidance of the described combination of movements which may lead to ACL injury.     

Development of a 3D finite element model of lens microcirculation

Background

It has been proposed that in the absence of a blood supply, the ocular lens operates an internal microcirculation system. This system delivers nutrients, removes waste products and maintains ionic homeostasis in the lens. The microcirculation is generated by spatial differences in membrane transport properties; and previously has been modelled by an equivalent electrical circuit and solved analytically. While effective, this approach did not fully account for all the anatomical and functional complexities of the lens. To encapsulate these complexities we have created a 3D finite element computer model of the lens.

Methods

Initially, we created an anatomically-correct representative mesh of the lens. We then implemented the Stokes and advective Nernst-Plank equations, in order to model the water and ion fluxes respectively. Next we complemented the model with experimentally-measured surface ionic concentrations as boundary conditions and solved it.

Results

Our model calculated the standing ionic concentrations and electrical potential gradients in the lens. Furthermore, it generated vector maps of intra- and extracellular space ion and water fluxes that are proposed to circulate throughout the lens. These fields have only been measured on the surface of the lens and our calculations are the first 3D representation of their direction and magnitude in the lens.

Conclusion

Values for steady state standing fields for concentration and electrical potential plus ionic and fluid fluxes calculated by our model exhibited broad agreement with observed experimental values. Our model of lens function represents a platform to integrate new experimental data as they emerge and assist us to understand how the integrated structure and function of the lens contributes to the maintenance of its transparency.     

A 3D finite element model of an implanted scapula: importance of a multiparametric validation using experimental data

In order to help to understand the loosening phenomenon around gleno?d prostheses, a 3D finite element model of a previously tested implanted scapula has been developed. The construction of the model was done using CT scans of the tested scapula. Different bone material properties were tested and shell elements or 8 nodes hexaedric elements were used to model the cortical bone. Surface contact elements were introduced on one hand between the bone and the lower part of the plate of the implant, and on the other, between the loading metallic ball and the upper surface of the implant. The results of the model were compared with those issued from in vitro experiments carried out on the same scapula. The evaluation of the model was done for nine cases of loading of 500 N distributed on the implant, in terms of strains (principal strains of six spots around peripheral cortex of the gleno?d) and displacement of four points positioned on the implant. The best configuration of the model presented here, fits with experiments for most of the strains (difference lower than 150microdef) but it seems to be still too stiff (mainly in the lower part). Nevertheless, we want, in this paper, to underline the importance of doing a multiparametric validation for such a model. Indeed, some models can give correct results for one case of loading but bad results for another kind of loading, some others can give good results for one kind of compared parameters (like strains for instance) but bad results for the other one (like displacements).     

3D finite element analysis of uniaxial cell stretching: from image to insight

Mechanical forces play an important role in many microbiological phenomena such as embryogenesis, regeneration, cell proliferation and differentiation. Micromanipulation of cells in a controlled environment is a widely used approach for understanding cellular responses with respect to external mechanical forces. While modern micromanipulation and imaging techniques provide useful optical information about the change of overall cell contours under the impact of external loads, the intrinsic mechanisms of energy and signal propagation throughout the cell structure are usually not accessible by direct observation. This work deals with the computational modelling and simulation of intracellular strain state of uniaxially stretched cells captured in a series of images. A nonlinear elastic finite element method on tetrahedral grids was applied for numerical analysis of inhomogeneous stretching of a rat embryonic fibroblast 52 (REF 52) using a simplified two-component model of a eukaryotic cell consisting of a stiffer nucleus surrounded by a softer cytoplasm. The difference between simulated and experimentally observed cell contours is used as a feedback criterion for iterative estimation of canonical material parameters of the two-component model such as stiffness and compressibility. Analysis of comparative simulations with varying material parameters shows that (i) the ratio between the stiffness of cell nucleus and cytoplasm determines intracellular strain distribution and (ii) large deformations result in increased stiffness and decreased compressibility of the cell cytoplasm. The proposed model is able to reproduce the evolution of the cellular shape over a sequence of observed deformations and provides complementary information for a better understanding of mechanical cell response.     

Smooth surface meshing for automated finite element model generation from 3D image data

Finite element (FE) modelling based on data from three-dimensional high-resolution computed tomography (CT) imaging systems provides a non-invasive method to assess structural mechanics. Automated mesh generation from these voxel based image data can be achieved by direct conversion to hexahedron elements, however these model representations have jagged edges. This paper proposes an automated method to generate smoothed FE meshes from voxel-based image data. Mesh fairing processes are utilized that allow constraints that control the smoothing process, and are computationally efficient. Surfaces of the mesh on the exterior, as well as interfaces between two tissues, can be smoothed by varying fairing parameters and constraint criteria. The method was tested on a variety of real and simulated three-dimensional data sets, resulting in both hexahedron and tetrahedron meshes. It was shown that the fairing process is linearly related to the number of smoothing iterations, and that peak stresses are reduced in FE simulations of the smoothed models. Although developed for micro-CT data sets, this fast and reliable mesh smoothing method could be applied to any three-dimensional image data where node and element connectivity have been defined.     

3D finite element model of the chinchilla ear for characterizing middle ear functions

Chinchilla is a commonly used animal model for research of sound transmission through the ear. Experimental measurements of the middle ear transfer function in chinchillas have shown that the middle ear cavity greatly affects the tympanic membrane (TM) and stapes footplate (FP) displacements. However, there is no finite element (FE) model of the chinchilla ear available in the literature to characterize the middle ear functions with the anatomical features of the chinchilla ear. This paper reports a recently completed 3D FE model of the chinchilla ear based on X-ray micro-computed tomography images of a chinchilla bulla. The model consisted of the ear canal, TM, middle ear ossicles and suspensory ligaments, and the middle ear cavity. Two boundary conditions of the middle ear cavity wall were simulated in the model as the rigid structure and the partially flexible surface, and the acoustic-mechanical coupled analysis was conducted with these two conditions to characterize the middle ear function. The model results were compared with experimental measurements reported in the literature including the TM and FP displacements and the middle ear input admittance in chinchilla ear. An application of this model was presented to identify the acoustic role of the middle ear septa—a unique feature of chinchilla middle ear cavity. This study provides the first 3D FE model of the chinchilla ear for characterizing the middle ear functions through the acoustic-mechanical coupled FE analysis.     

A novel methodology for generating 3D finite element models of the hip from 2D radiographs

Finite element (FE) modelling has been proposed as a tool for estimating fracture risk and patient-specific FE models are commonly based on computed tomography (CT). Here, we present a novel method to automatically create personalised 3D models from standard 2D hip radiographs. A set of geometrical parameters of the femur were determined from seven a–p hip radiographs and compared to the 3D femoral shape obtained from CT as training material; the error in reconstructing the 3D model from the 2D radiographs was assessed. Using the geometry parameters as the input, the 3D shape of another 21 femora was built and meshed, separating a cortical and trabecular compartment. The material properties were derived from the homogeneity index assessed by texture analysis of the radiographs, with focus on the principal tensile and compressive trabecular systems. The ability of these FE models to predict failure load as determined by experimental biomechanical testing was evaluated and compared to the predictive ability of DXA. The average reconstruction error of the 3D models was 1.77 mm (±1.17 mm), with the error being smallest in the femoral head and neck, and greatest in the trochanter. The correlation of the FE predicted failure load with the experimental failure load was r²=64% for the reconstruction FE model, which was significantly better (p<0.05) than that for DXA (r²=24%). This novel method for automatically constructing a patient-specific 3D finite element model from standard 2D radiographs shows encouraging results in estimating patient-specific failure loads.