Computational efficiency of numerical approximations of tangent moduli for finite element implementation of a fiber-reinforced hyperelastic material model

In this study, we evaluated computational efficiency of finite element (FE) simulations when a numerical approximation method was used to obtain the tangent moduli. A fiber-reinforced hyperelastic material model for nearly incompressible soft tissues was implemented for 3D solid elements using both the approximation method and the closed-form analytical method, and validated by comparing the components of the tangent modulus tensor (also referred to as the material Jacobian) between the two methods. The computational efficiency of the approximation method was evaluated with different perturbation parameters and approximation schemes, and quantified by the number of iteration steps and CPU time required to complete these simulations. From the simulation results, it can be seen that the overall accuracy of the approximation method is improved by adopting the central difference approximation scheme compared to the forward Euler approximation scheme. For small-scale simulations with about 10,000 DOFs, the approximation schemes could reduce the CPU time substantially compared to the closed-form solution, due to the fact that fewer calculation steps are needed at each integration point. However, for a large-scale simulation with about 300,000 DOFs, the advantages of the approximation schemes diminish because the factorization of the stiffness matrix will dominate the solution time. Overall, as it is material model independent, the approximation method simplifies the FE implementation of a complex constitutive model with comparable accuracy and computational efficiency to the closed-form solution, which makes it attractive in FE simulations with complex material models.     

Comparison of micro-level and continuum-level voxel models of the proximal femur

Continuum-level finite element (FE) models became standard computational tools for the evaluation of bone mechanical behavior from in vivo computed tomography scans. Such scans do not account for the anisotropy of the bone. Instead, local mechanical properties in the continuum-level FE models are assumed isotropic and are derived from bone density, using statistical relationships. Micro-FE models, on the other hand, incorporate the anisotropic structure in detail. This study aimed to quantify the effects of assumed isotropy, by comparing continuum-level voxel models of a healthy and a severely osteoporotic proximal femur with recently analyzed micro-FE models of the same bones. The micro-model element size was coarsened to generate continuum FE models with two different element sizes (0.64 and 3.04 mm) and two different density–modulus relationships found in the literature for wet and ash density. All FE models were subjected to the same boundary conditions that simulated a fall to the side, and the stress and strain distributions, model stiffness and yield load were compared. The results indicated that the stress and strain distributions could be reproduced well with the continuum models. The smallest differences between the continuum-level model and micro-level model predictions of the stiffness and yield load were obtained with the coarsest element size. Better results were obtained for both continuum-element sizes when isotropic moduli were based on ash density rather than wet density.     

Combined finite element and multibody musculoskeletal investigation of a fractured clavicle with reconstruction plate

This paper addresses the evaluation of clavicle fixation devices, by means of computational models. The aim was to develop a method for comparison of stress distribution in various fixation devices, to determine whether the use of multibody musculoskeletal input in such model is applicable and to report the approach. The focus was on realistic loading and the motivation for the work is that the treatment can be enhanced by a better understanding of the loading of the clavicle and fixation device. The method can be used to confirm the strength of customised plates, for optimisation of new plates and to complement experimental studies. A finite element (FE) mesh of the clavicle geometry was created from computed tomography data and imported into the FE solver where the model was subjected to muscle forces and other boundary conditions from a multibody musculoskeletal model performing a typical activity of daily life. A reconstruction plate and screws were also imported into the model. The combination models returned stresses and displacements of plausible magnitudes in all included parts and the result, upon further development and validation, may serve as a design guideline for improved clavicle fixation.     

An efficient numerical method for distributed-loop models of the urine concentrating mechanism

In this study we describe an efficient numerical method, based on the semi-Lagrangian (SL) semi-implicit (SI) method and Newton's method, for obtaining steady-state (SS) solutions of equations arising in distributed-loop models of the urine concentrating mechanism. Dynamic formulations of these models contain large systems of coupled hyperbolic partial differential equations (PDEs). The SL method advances the solutions of these PDEs in time by integrating backward along flow trajectories, thus allowing large time steps while maintaining stability. The SI approach controls stiffness arising from transtubular transport terms by averaging these terms in time along flow trajectories. An approximate SS solution of a dynamic formulation obtained via the SLSI method can be used as an initial guess for a Newton-type solver, which rapidly converges to a highly accurate numerical approximation to the solution of the ordinary differential equations that arise in the corresponding SS model formulation. In general, it is difficult to specify a priori for a Newton-type solver an initial guess that falls within the radius of convergence; however, the initial guess generated by solving the dynamic formulation via the SLSI method can be made sufficiently close to the SS solution to avoid numerical instability. The combination of the SLSI method and the Newton-type solver generates stable and accurate solutions with substantially reduced computation times, when compared to previously applied dynamic methods.     

Accounting for spatial variation of trabecular anisotropy with subject-specific finite element modeling moderately improves predictions of local subchondral bone stiffness at the proximal tibia

IntroductionPreviously, a finite element (FE) model of the proximal tibia was developed and validated against experimentally measured local subchondral stiffness. This model indicated modest predictions of stiffness (R² = 0.77, normalized root mean squared error (RMSE%) = 16.6%). Trabecular bone though was modeled with isotropic material properties despite its orthotropic anisotropy. The objective of this study was to identify the anisotropic FE modeling approach which best predicted (with largest explained variance and least amount of error) local subchondral bone stiffness at the proximal tibia.MethodsLocal stiffness was measured at the subchondral surface of 13 medial/lateral tibial compartments using in situ macro indentation testing. An FE model of each specimen was generated assuming uniform anisotropy with 14 different combinations of cortical- and tibial-specific density-modulus relationships taken from the literature. Two FE models of each specimen were also generated which accounted for the spatial variation of trabecular bone anisotropy directly from clinical CT images using grey-level structure tensor and Cowin’s fabric-elasticity equations. Stiffness was calculated using FE and compared to measured stiffness in terms of R² and RMSE%.ResultsThe uniform anisotropic FE model explained 53–74% of the measured stiffness variance, with RMSE% ranging from 12.4 to 245.3%. The models which accounted for spatial variation of trabecular bone anisotropy predicted 76–79% of the variance in stiffness with RMSE% being 11.2–11.5%.ConclusionsOf the 16 evaluated finite element models in this study, the combination of Synder and Schneider (for cortical bone) and Cowin’s fabric-elasticity equations (for trabecular bone) best predicted local subchondral bone stiffness.     

A coupled sharp-interface immersed boundary-finite-element method for flow-structure interaction with application to human phonation

A new flow-structure interaction method is presented, which couples a sharp-interface immersed boundary method flow solver with a finite-element method based solid dynamics solver. The coupled method provides robust and high-fidelity solution for complex flow-structure interaction (FSI) problems such as those involving three-dimensional flow and viscoelastic solids. The FSI solver is used to simulate flow-induced vibrations of the vocal folds during phonation. Both two- and three-dimensional models have been examined and qualitative, as well as quantitative comparisons, have been made with established results in order to validate the solver. The solver is used to study the onset of phonation in a two-dimensional laryngeal model and the dynamics of the glottal jet in a three-dimensional model and results from these studies are also presented.     

Finite element modelling approaches for well-ordered porous metallic materials for orthopaedic applications: cost effectiveness and geometrical considerations

The mechanical properties of well-ordered porous materials are related to their geometrical parameters at the mesoscale. Finite element (FE) analysis is a powerful tool to design well-ordered porous materials by analysing the mechanical behaviour. However, FE models are often computationally expensive. This article aims to develop a cost-effective FE model to simulate well-ordered porous metallic materials for orthopaedic applications. Solid and beam FE modelling approaches are compared, using finite size and infinite media models considering cubic unit cell geometry. The model is then applied to compare two unit cell geometries: cubic and diamond. Models having finite size provide similar results than the infinite media model approach for large sample sizes. In addition, these finite size models also capture the influence of the boundary conditions on the mechanical response for small sample sizes. The beam FE modelling approach showed little computational cost and similar results to the solid FE modelling approach. Diamond unit cell geometry appeared to be more suitable for orthopaedic applications than the cubic unit cell geometry.     

Numerical Models for the Simulation of Flexible Artificial Heart Valves: Part I - Computational Methods

Abstract

A numerical model of the coupled motion of a flexing surface in a high Reynolds number flow is presented for the simulation of flexible polyurethane heart valves in the aortic position. This is achieved by matching a Lagrangian dynamic leaflet model with a panel method based flow solver. The two models are coupled via the time-dependent pressure field using the unsteady Bernoulli equation.

Incorporation of sub-cycling in the dynamic model equations and fast pre conditioning techniques in the panel method solver yields efficient convergence and near real-time simulations of valve motion. The generality of dynamic model allows different material properties and/or geometries to be studied easily and interactively. This interactivity is realized by embedding the models within a design environment created using the software IRIS Explorer ^TM.

Two flow domains are developed, an infinite domain and an internal domain using conformal mapping theory. In addition bending stress on the valve is computed using a simple stress model based on spline and circle equation techniques.     

HR-pQCT-based homogenised finite element models provide quantitative predictions of experimental vertebral body stiffness and strength with the same accuracy as μFE models

This study validated two different high-resolution peripheral quantitative computer tomography (HR-pQCT)-based finite element (FE) approaches, enhanced homogenised continuum-level (hFE) and micro-finite element (μFE) models, by comparing them with compression test results of vertebral body sections. Thirty-five vertebral body sections were prepared by removing endplates and posterior elements, scanned with HR-pQCT and tested in compression up to failure. Linear hFE and μFE models were created from segmented and grey-level CT images, and apparent model stiffness values were compared with experimental stiffness as well as strength results. Experimental and numerical apparent elastic properties based on grey-level/segmented CT images (N=35) correlated well for μFE (r2=0.748/0.842) and hFE models (r2=0.741/0.864). Vertebral section stiffness values from the linear μFE/hFE models estimated experimental ultimate apparent strength very well (r2=0.920/0.927). Calibrated hFE models were able to predict quantitatively apparent stiffness with the same accuracy as μFE models. However, hFE models needed no back-calculation of a tissue modulus or any kind of fitting and were computationally much cheaper.     

Study of age-related changes in Middle ear transfer function

Abstract

Osteoporosis (OP) is common with advancing age. Several studies have shown a strong correlation between OP and otosclerosis. However, no studies have investigated OP of the malleus, incus or stapes in the human middle ear, its effect on middle ear transfer function. Here, we investigate whether these three ossicles develop OP, and how this affects middle ear transfer function. The effect of OP on middle ear transfer function was investigated in simulations based on a finite element (FE) method. First, the FE model used in our previous study was refined, and optimized by introducing viscoelastic properties to selected soft tissues of the middle ear. Then, the FE model was used to simulate OP of the three ossicles and assess its influence on middle ear transfer function. Other possible age-related changes, such as stiffness of the joints or ligaments in the middle ear, were also investigated. The results indicated that OP of the ossicles could increase the high frequency displacement of both the umbo and stapes footplate (FP). However, the stiffness of the middle ear soft tissue can lead to the decrease of middle ear gain at lower frequencies. Furthermore, loosening of these joints or ligaments could increase displacement of the umbo and stapes FP. In conclusion, although age-related hearing loss is most commonly conceived of as sensorineural hearing loss (SNHL), we found that age-related changes may also include OP and changes in joint stiffness, but these will have little effect on middle ear transfer function in elderly people.     

Validation of density–elasticity relationships for finite element modeling of human pelvic bone by modal analysis

In total hip arthroplasty and particularly in revision surgery, computer assisted pre-operative prediction of the best possible anchorage strategy for implant fixation would be a great help to the surgeon. Computer simulation relies on validated numerical models. In the current study, three density–elasticity relationships (No. 1–3) from the literature for inhomogeneous material parameter assignment from CT data in automated finite element (FE) modeling of long bones were evaluated for their suitability for FE modeling of human pelvic bone. Numerical modal analysis was conducted on 10 FE models of hemipelvic bone specimens and compared to the gold standard provided by experimental modal analysis results from a previous in-vitro study on the same specimens. Overall, calculated resonance frequencies came out lower than measured values. Magnitude of mean relative deviation of numerical resonance frequencies with regard to measured values is lowest for the density–elasticity relationship No. 3 (−15.9%) and considerably higher for both density–elasticity relationships No. 1 (−41.1%) and No. 2 (−45.0%). Mean MAC values over all specimens amount to 77.8% (No. 1), 78.5% (No. 2), and 83.0% (No. 3). MAC results show, that mode shapes are only slightly influenced by material distribution. Calculated resonance frequencies are generally lower than measured values, which indicates, that numerical models lack stiffness. Even when using the best suited (No. 3) out of three investigated density–elasticity relationships, in FE modeling of pelvic bone a considerable underestimation of model stiffness has to be taken into account.     

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.     

Numerical models for the simulation of flexible artificial heart valves: part I--computational methods

A numerical model of the coupled motion of a flexing surface in a high Reynolds number flow is presented for the simulation of flexible polyurethane heart valves in the aortic position. This is achieved by matching a Lagrangian dynamic leaflet model with a panel method based flow solver. The two models are coupled via the time-dependent pressure field using the unsteady Bernoulli equation. Incorporation of sub-cycling in the dynamic model equations and fast pre conditioning techniques in the panel method solver yields efficient convergence and near real-time simulations of valve motion. The generality of dynamic model allows different material properties and/or geometries to be studied easily and interactively. This interactivity is realized by embedding the models within a design environment created using the software IRIS Explorer. Two flow domains are developed, an infinite domain and an internal domain using conformal mapping theory. In addition bending stress on the valve is computed using a simple stress model based on spline and circle equation techniques.     

Effects of stiffness and volume on the transit time of an erythrocyte through a slit

By using a fully coupled fluid–cell interaction model, we numerically simulate the dynamic process of a red blood cell passing through a slit driven by an incoming flow. The model is achieved by combining a multiscale model of the composite cell membrane with a boundary element fluid dynamics model based on the Stokes flow assumption. Our concentration is on the correlation between the transit time (the time it takes to finish the whole translocation process) and different conditions (flow speed, cell orientation, cell stiffness, cell volume, etc.) that are involved. According to the numerical prediction (with some exceptions), the transit time rises as the cell is stiffened. It is also highly sensitive to volume increase inside the cell. In general, even slightly swollen cells (i.e., the internal volume is increased while the surface area of the cell kept unchanged) travel dramatically slower through the slit. For these cells, there is also an increased chance of blockage.     

Estimation of conditions evoking fracture in finger bones under pinch loading based on finite element analysis

A finger finite element (FE) model was created from CT images of a Japanese male in order to obtain a shape-biofidelic model. Material properties and articulation characteristics of the model were taken from the literature. To predict bone fracture and realistically represent the fracture pattern under various loading conditions, the ESI-Wilkins-Kamoulakos rupture model in PAM-CRASH (ESI Group S.A., Paris, France) was utilized in this study with parameter values of the rupture model determined by compression testing and simulation of porcine fibula. A finger pinch simulation was then conducted to validate the finger FE model. The force-displacement curve and fracture load from the pinch simulation was compared to the result of finger pinch test using cadavers. Simulation results are coincident with the test result, indicating that the finger FE model can be used in an analysis of finger bone fracture during pinch accident. With this model, several pinch simulations were conducted with different pinching object’s stiffness and pinching energy. Conditions for evoking finger bone fracture under pinch loading were then estimated based on these results. This study offers a novel method to predict possible hazards of manufactured goods during the design process, thus finger injury due to pinch loading can be avoided.     

Deformation of articular cartilage during static loading of a knee joint – Experimental and finite element analysis

Novel conical beam CT-scanners offer high resolution imaging of knee structures with i.a. contrast media, even under weight bearing. With this new technology, we aimed to determine cartilage strains and meniscal movement in a human knee at 0, 1, 5, and 30 min of standing and compare them to the subject-specific 3D finite element (FE) model. The FE model of the volunteer?s knee, based on the geometry obtained from magnetic resonance images, was created to simulate the creep. The effects of collagen fibril network stiffness, nonfibrillar matrix modulus, permeability and fluid flow boundary conditions on the creep response in cartilage were investigated. In the experiment, 80% of the maximum strain in cartilage developed immediately, after which the cartilage continued to deform slowly until the 30 min time point. Cartilage strains and meniscus movement obtained from the FE model matched adequately with the experimentally measured values. Reducing the fibril network stiffness increased the mean strains substantially, while the creep rate was primarily influenced by an increase in the nonfibrillar matrix modulus. Changing the initial permeability and preventing fluid flow through noncontacting surfaces had a negligible effect on cartilage strains. The present results improve understanding of the mechanisms controlling articular cartilage strains and meniscal movements in a knee joint under physiological static loading. Ultimately a validated model could be used as a noninvasive diagnostic tool to locate cartilage areas at risk for degeneration.     

Transmembrane voltage induced on altered erythrocyte shapes exposed to RF fields

In this article, the transmembrane voltage induced on erythrocyte, codocyte, ovalocyte and spherocyte cell models exposed to a linearly polarised electromagnetic plane wave of frequency 1800 MHz is calculated. For this purpose, a finite element (FE) numerical technique with adaptive meshing is used. The results show that the value of the induced voltage on the original erythrocyte shape is higher than the one observed on the rest of the altered cell geometries studied. The erythrocyte shape and the membrane electric permittivity are shown to play a fundamental role on the values of the induced transmembrane voltage.     

Intervertebral disc creep behavior assessment through an open source finite element solver

Degenerative Disc Disease (DDD) is one of the largest health problems faced worldwide, based on lost working time and associated costs. By means of this motivation, this work aims to evaluate a biomimetic Finite Element (FE) model of the Intervertebral Disc (IVD). Recent studies have emphasized the importance of an accurate biomechanical modeling of the IVD, as it is a highly complex multiphasic medium. Poroelastic models of the disc are mostly implemented in commercial finite element packages with limited access to the algorithms. Therefore, a novel poroelastic formulation implemented on a home-developed open source FE solver is briefly addressed throughout this paper. The combination of this formulation with biphasic osmotic swelling behavior is also taken into account. Numerical simulations were devoted to the analysis of the non-degenerated human lumbar IVD time-dependent behavior. The results of the tests performed for creep assessment were inside the scope of the experimental data, with a remarkable improvement of the numerical accuracy when compared with previously published results obtained with ABAQUS^®. In brief, this in-development open-source FE solver was validated with literature experimental data and aims to be a valuable tool to study the IVD biomechanics and DDD mechanisms.     

Evaluating a suitable level of model complexity for finite element analysis of the intact acetabulum

To enable large-scale multi-factorial finite element (FE) studies, the FE models used must be as computationally efficient as is feasible, while maintaining a suitable level of definition. The present study seeks to find an optimum level of model complexity for use in such large-scale studies by investigating which model attributes are most influential over the chosen model outputs of principal stress and strain in the intact acetabulum. A multi-factorial sensitivity study was carried out using 128 FE models, representing combinations of the following variables: bone stiffness distribution, imposed muscle loading, boundary condition location, hip joint contact conditions and patient's bone anatomy. The relative sensitivity of each input factor was analysed, and it was concluded that the optimum level of model definition must include CT-dependent trabecular bone properties and a sliding interface at the hip joint. It was found that it was not essential to describe the ligamentous sacroiliac and pubic symphysis joints; these could be rigidly fixed in space; and for the normal walking load case, muscle forces may be neglected. It was also concluded that a variety of bone anatomies should be included in a multi-factorial analysis if results are to be inferred for a wider population.     

An anisotropic elastic-viscoplastic damage model for bone tissue

A new anisotropic elastic-viscoplastic damage constitutive model for bone is proposed using an eccentric elliptical yield criterion and nonlinear isotropic hardening. A micromechanics-based multiscale homogenization scheme proposed by Reisinger et al. is used to obtain the effective elastic properties of lamellar bone. The dissipative process in bone is modeled as viscoplastic deformation coupled to damage. The model is based on an orthotropic ecuntric elliptical criterion in stress space. In order to simplify material identification, an eccentric elliptical isotropic yield surface was defined in strain space, which is transformed to a stress-based criterion by means of the damaged compliance tensor. Viscoplasticity is implemented by means of the continuous Perzyna formulation. Damage is modeled by a scalar function of the accumulated plastic strain ${D(\kappa)}$ , reducing all element s of the stiffness matrix. A polynomial flow rule is proposed in order to capture the rate-dependent post-yield behavior of lamellar bone. A numerical algorithm to perform the back projection on the rate-dependent yield surface has been developed and implemented in the commercial finite element solver Abaqus/Standard as a user subroutine UMAT. A consistent tangent operator has been derived and implemented in order to ensure quadratic convergence. Correct implementation of the algorithm, convergence, and accuracy of the tangent operator was tested by means of strain- and stress-based single element tests. A finite element simulation of nano- indentation in lamellar bone was finally performed in order to show the abilities of the newly developed constitutive model.