Considerations for reporting finite element analysis studies in biomechanics

Simulation-based medicine and the development of complex computer models of biological structures is becoming ubiquitous for advancing biomedical engineering and clinical research. Finite element analysis (FEA) has been widely used in the last few decades to understand and predict biomechanical phenomena. Modeling and simulation approaches in biomechanics are highly interdisciplinary, involving novice and skilled developers in all areas of biomedical engineering and biology. While recent advances in model development and simulation platforms offer a wide range of tools to investigators, the decision making process during modeling and simulation has become more opaque. Hence, reliability of such models used for medical decision making and for driving multiscale analysis comes into question. Establishing guidelines for model development and dissemination is a daunting task, particularly with the complex and convoluted models used in FEA. Nonetheless, if better reporting can be established, researchers will have a better understanding of a model's value and the potential for reusability through sharing will be bolstered. Thus, the goal of this document is to identify resources and considerate reporting parameters for FEA studies in biomechanics. These entail various levels of reporting parameters for model identification, model structure, simulation structure, verification, validation, and availability. While we recognize that it may not be possible to provide and detail all of the reporting considerations presented, it is possible to establish a level of confidence with selective use of these parameters. More detailed reporting, however, can establish an explicit outline of the decision-making process in simulation-based analysis for enhanced reproducibility, reusability, and sharing.     

Biomechanical measurements of torsion-tension coupling in human cadaveric femurs

The mechanical behavior of human femurs has been described in the literature with regard to torsion and tension but only as independent measurements. However, in this study, human femurs were subjected to torsion to determine if a simultaneous axial tensile load was generated. Fresh frozen human femurs (n=25) were harvested and stripped of soft tissue. Each femur was mounted rigidly in a specially designed test jig and remained at a fixed axial length during all experiments. Femurs were subjected to external and internal rotation applied at a constant angulation rate of 0.1 deg/s to a maximum torque of 12?N?m. Applied torque and generated axial tension were monitored simultaneously. Outcome measurements were extracted from torsion-versus-tension graphs. There was a strong relationship between applied torsion and the resulting tension for external rotation tests (torsion/tension ratio=551.7±283.8?mm, R(2)=0.83±0.20, n=25), internal rotation tests (torsion/tension ratio=495.3±233.1?mm, R(2)=0.87±0.17, n=24), left femurs (torsion/tension ratio=542.2±262.4?mm, R(2)=0.88±0.13, n=24), and right femurs (torsion/tension ratio=506.7±260.0?mm, R(2)=0.82±0.22, n=25). No statistically significant differences were found for external versus internal rotation groups or for left versus right femurs when comparing torsion/tension ratios (p=0.85) or R(2) values (p=0.54). A strongly coupled linear relationship between torsion and tension for human femurs was exhibited. This suggests an interplay between these two factors during activities of daily living and injury processes.     

The effect of boundary condition on the biomechanics of a human pelvic joint under an axial compressive load: a three-dimensional finite element model

The finite element (FE) model of the pelvic joint is helpful for clinical diagnosis and treatment of pelvic injuries. However, the effect of an FE model boundary condition on the biomechanical behavior of a pelvic joint has not been well studied. The objective of this study was to study the effect of boundary condition on the pelvic biomechanics predictions. A 3D FE model of a pelvis using subject-specific estimates of intact bone structures, main ligaments and bone material anisotropy by computed tomography (CT) gray value was developed and validated by bone surface strains obtained from rosette strain gauges in an in vitro pelvic experiment. Then three FE pelvic models were constructed to analyze the effect of boundary condition, corresponding to an intact pelvic joint, a pelvic joint without sacroiliac ligaments and a pelvic joint without proximal femurs, respectively. Vertical load was applied to the same pelvis with a fixed prosthetic femoral stem and the same load was simulated in the FE model. A strong correlation coefficient (R(2)=0.9657) was calculated, which indicated a strong correlation between the FE analysis and experimental results. The effect of boundary condition changes on the biomechanical response depended on the anatomical location and structure of the pelvic joint. It was found that acetabulum fixed in all directions with the femur removed can increase the stress distribution on the acetabular inner plate (approximately double the original values) and decrease that on the superior of pubis (from 7 MPa to 0.6 MPa). Taking sacrum and ilium as a whole, instead of sacroiliac and iliolumber ligaments, can influence the stress distribution on ilium and pubis bone vastly. These findings suggest pelvic biomechanics is very dependent on the boundary condition in the FE model.     

Experimental validation of a finite element model of a human cadaveric tibia

Finite element (FE) models of long bones are widely used to analyze implant designs. Experimental validation has been used to examine the accuracy of FE models of cadaveric femurs; however, although convergence tests have been carried out, no FE models of an intact and implanted human cadaveric tibia have been validated using a range of experimental loading conditions. The aim of the current study was to create FE models of a human cadaveric tibia, both intact and implanted with a unicompartmental knee replacement, and to validate the models against results obtained from a comprehensive set of experiments. Seventeen strain rosettes were attached to a human cadaveric tibia. Surface strains and displacements were measured under 17 loading conditions, which consisted of axial, torsional, and bending loads. The tibia was tested both before and after implantation of the knee replacement. FE models were created based on computed tomography (CT) scans of the cadaveric tibia. The models consisted of ten-node tetrahedral elements and used 600 material properties derived from the CT scans. The experiments were simulated on the models and the results compared to experimental results. Experimental strain measurements were highly repeatable and the measured stiffnesses compared well to published results. For the intact tibia under axial loading, the regression line through a plot of strains predicted by the FE model versus experimentally measured strains had a slope of 1.15, an intercept of 5.5 microstrain, and an R(2) value of 0.98. For the implanted tibia, the comparable regression line had a slope of 1.25, an intercept of 12.3 microstrain, and an R(2) value of 0.97. The root mean square errors were 6.0% and 8.8% for the intact and implanted models under axial loads, respectively. The model produced by the current study provides a tool for simulating mechanical test conditions on a human tibia. This has considerable value in reducing the costs of physical testing by pre-selecting the most appropriate test conditions or most favorable prosthetic designs for final mechanical testing. It can also be used to gain insight into the results of physical testing, by allowing the prediction of those variables difficult or impossible to measure directly.     

Deformation and stress distribution of the human foot after plantar ligaments release: A cadaveric study and finite element analysis

The majority of foot deformities are related to arch collapse or instability, especially the longitudinal arch. Although the relationship between the plantar fascia and arch height has been previously investigated, the stress distribution remains unclear. The aim of this study was to explore the role of the plantar ligaments in foot arch biomechanics. We constructed a geometrical detailed three-dimensional (3-D) finite element (FE) model of the human foot and ankle from computer tomography images. The model comprised the majority of joints in the foot as well as bone segments, major ligaments, and plantar soft tissue. Release of the plantar fascia and other ligaments was simulated to evaluate the corresponding biomechanical effects on load distribution of the bony and ligamentous structures. These intrinsic ligaments of the foot arch were sectioned to simulate different pathologic situations of injury to the plantar ligaments, and to explore bone segment displacement and stress distribution. The validity of the 3-D FE model was verified by comparing results with experimentally measured data via the displacement and von Mise stress of each bone segment. Plantar fascia release decreased arch height, but did not cause total collapse of the foot arch. The longitudinal foot arch was lost when all the four major plantar ligaments were sectioned simultaneously. Plantar fascia release was compromised by increased strain applied to the plantar ligaments and intensified stress in the midfoot and metatarsal bones. Load redistribution among the centralized metatarsal bones and focal stress relief at the calcaneal insertion were predicted. The 3-D FE model indicated that plantar fascia release may provide relief of focal stress and associated heel pain. However, these operative procedures may pose a risk to arch stability and clinically may produce dorsolateral midfoot pain. The initial strategy for treating plantar fasciitis should be non-operative.     

Development of a finite element model of the tibia for short-duration high-force axial impact loading

Finite element (FE) models can allow computer simulations of impact loading, providing a useful companion to cadaveric testing. These models allow injury evaluations to be conducted under a variety of conditions, but must be validated against experimental data. An FE model of a cadaveric tibia was developed using geometry from CT scans, and the quality of the mesh was evaluated. Loading and boundary conditions from experimental tests were simulated, and the model was optimised to best represent the response of natural bone to impacts. The model was shown to have good agreement for impact force, duration, impulse and strain during simulation of three non-injurious and one injurious axial impact when compared with experimental test data for the specimen. Failure criteria were evaluated for their ability to predict fracture. This model of the tibia can be used for future injury prediction assessment studies.     

Subject-specific finite element modelling of the human foot complex during walking: sensitivity analysis of material properties,boundary and loading conditions

The objective of this study was to develop and validate a subject-specific framework for modelling the human foot. This was achieved by integrating medical image-based finite element modelling, individualised multi-body musculoskeletal modelling and 3D gait measurements. A 3D ankle–foot finite element model comprising all major foot structures was constructed based on MRI of one individual. A multi-body musculoskeletal model and 3D gait measurements for the same subject were used to define loading and boundary conditions. Sensitivity analyses were used to investigate the effects of key modelling parameters on model predictions. Prediction errors of average and peak plantar pressures were below 10% in all ten plantar regions at five key gait events with only one exception (lateral heel, in early stance, error of 14.44%). The sensitivity analyses results suggest that predictions of peak plantar pressures are moderately sensitive to material properties, ground reaction forces and muscle forces, and significantly sensitive to foot orientation. The maximum region-specific percentage change ratios (peak stress percentage change over parameter percentage change) were 1.935–2.258 for ground reaction forces, 1.528–2.727 for plantar flexor muscles and 4.84–11.37 for foot orientations. This strongly suggests that loading and boundary conditions need to be very carefully defined based on personalised measurement data.     

Quantitative comparison of ligament formulation and pre-strain in finite element analysis of the human lumbar spine

Data has been published that quantifies the nonlinear, anisotropic material behaviour and pre-strain behaviour of the anterior longitudinal, supraspinous (SSL), and interspinous ligaments of the human lumbar spine. Additionally, data has been published on localized material properties of the SSL. These results have been incrementally incorporated into a previously validated finite element model of the human lumbar spine. Results suggest that the effects of increased ligament model fidelity on bone strain energy were moderate and the effects on disc pressure were slight, and do not justify a change in modelling strategy for most clinical applications. There were significant effects on the ligament stresses of the ligaments that were directly modified, suggesting that these phenomena should be included in FE models where ligament stresses are the desired metric.     

Poroelastic finite element analysis of a bone specimen under cyclic loading

It had been suggested that the fluid embodied in bone lacunar-canalicular porosity may play an important role in bone remodelling [Weinbaum et al., 1994. Journal of Biomechanics 27, 339-360]. In this paper a finite element model of a poroelastic prismatic solid of rectangular cross-section is considered to simulate bone behaviour, precisely as in the previous work by Zhang and Cowin [Zhang and Cowin, 1994. Journal of Mechanical Physics of Solids 42, 1575-1599]. This solid is subject to combined cyclic axial and bending loads at its end. The objectives of the study are: (1) to verify the accuracy of the simplifying hypotheses underlying the analytical solutions established by the above authors; (2) to provide further insight into the behaviour of that solid; (3) to test the advantages in generality and versatility and the computing costs of general-purpose finite element codes in poroelastic analysis. The study is parametric with respect to the fluid leakage coefficient, to the ratio of the bending moment and axial load, and to the ratio of the characteristic relaxation time of the pore pressure over the excitation period. Results show that, for all the cases considered, the pore pressure distribution along the section height of the poroelastic beam exhibits a very good matching with previous analytical results. Stresses transversal with respect to the beam axis (assumed as constant or zero in previous analytical solutions) are evaluated. The analysis pointed out that: (1) the effects due to end-loads with zero resultants practically extinguish within a distance from the beam end almost equal to a typical length of the cross-section; (2) cross-sections remain plane above that distance; (3) the transversal total stresses are three orders of magnitude lower than axial stress.     

Modeling the biomechanics of the mandible: a three-dimensional finite element study.

Three-dimensional finite element models of a partially edentulated human mandible were generated to calculate the mechanical response to simulated isometric biting and mastication loads. The level of mesh refinement was established via a convergence test and showed that a model with over 30,000 degrees of freedom was required to obtain analysis accuracy. The functional loading cases included muscle loading based on an algorithm that assigns muscle forces in accordance with muscle cross-sectional area, while maintaining static equilibrium. Results were found for isometric application of unilateral and bilateral bite and mastication loading, and two different sets of displacement boundary conditions were imposed at the condyles. The mechanical response is shown in terms of displacements, principal strains, and a new measure called the 'mechanical intensity scalar'. For each load case studied, there was substantial bending in the molar region of the corpus and high tensile strains in the anterior portion of the ramus.     

A dynamic finite element analysis of impulsive loading of the extension-splinted rabbit knee

A dynamic nonlinear finite element model was developed to study juxtarticular stresses in the splinted rabbit knee, an established laboratory model for creating osteoarthrosis due to impulsive loading. Plane strain finite element results were validated by comparison with corresponding experimental data. Parametric effects studied included the input tibial displacement speed, the local bone density distribution, and the modulus of cartilage and subchondral bone. While the computed resultant contact force magnitude was sensitive to a number of model parameters, the stress patterns, when normalized to a given resultant force magnitude, were not. Despite comparable force peaks, the finite element results showed approximately six-fold higher effective strain rate levels for a severely impulsive loading protocol known to induce rapid osteoarthrosis, versus those for a mildly impulsive loading protocol not usually associated with cartilage damage. A propensity for elevated shear in the deep cartilage layer near the contact periphery, observed in nearly all computed stress distributions, is consistent with previous experimental findings of fissuring at that level in the impulsively loaded rabbit knee.     

Experimental and finite element analysis of the rat ulnar loading model-correlations between strain and bone formation following fatigue loading

The rat forelimb compression model has been used widely to study bone response to mechanical loading. We used strain gages to assess load sharing between the ulna and radius in the forelimb of adult Fisher rats. We used histology and peripheral quantitative computed tomography (pQCT) to quantify ulnar bone formation 12 days after in vivo fatigue loading. Lastly, we developed a finite element model of the ulna to predict the pattern of surface strains during compression. Our findings indicate that at the mid-shaft the ulna carries 65% of the applied compressive force on the forelimb. We observed large variations in fatigue-induced bone formation over the circumference and length of the ulna. Bone formation was greatest 1-2 mm distal to the mid-shaft. At the mid-shaft, we observed woven bone formation that was greatest medially. Finite element analysis indicated a strain pattern consistent with a compression-bending loading mode, with the greatest strains occurring in compression on the medial surface and lesser tensile strains occurring laterally. A peak strain of -5190 microepsilon (for 13.3N forelimb compression) occurred 1-2 mm distal to the mid-shaft. The pattern of bone formation in the longitudinal direction was highly correlated to the predicted peak compressive axial strains at seven cross-sections (r2 = 0.89, p = 0.014). The in-plane pattern of bone formation was poorly correlated to the predicted magnitude of axial strain at 51 periosteal locations (r2 = 0.21, p < 0.001), because the least bone formation was observed where tensile strains were highest. These findings indicate that the magnitude of bone formation after fatigue loading is greatest in regions of high compressive strain.     

Effect of partial-thickness tear on loading capacities of the supraspinatus tendon: a finite element analysis

Partial-thickness tears of the supraspinatus tendon frequently occur at its insertion on the greater tubercule of the humerus, causing pain and reduced strength and range of motion. The goal of this work was to quantify the loss of loading capacity due to tendon tears at the insertion area. A finite element model of the supraspinatus tendon was developed using in vivo magnetic resonance images data. The tendon was represented by an anisotropic hyperelastic constitutive law identified with experimental measurements. A failure criterion was proposed and calibrated with experimental data. A partial-thickness tear was gradually increased, starting from the deep articular-sided fibres. For different values of tendon tear thickness, the tendon was mechanically loaded up to failure. The numerical model predicted a loss in loading capacity of the tendon as the tear thickness progressed. Tendon failure was more likely when the tendon tear exceeded 20%. The predictions of the model were consistent with experimental studies. Partial-thickness tears below 40% tear are sufficiently stable to persist physiotherapeutic exercises. Above 60% tear surgery should be considered to restore shoulder strength.     

Mixed-mode loading of the cement-bone interface: a finite element study

While including the cement-bone interface of complete cemented hip reconstructions is crucial to correctly capture their response, its modelling is often overly simplified. In this study, the mechanical mixed-mode response of the cement-bone interface is investigated, taking into account the effects of the well-defined microstructure that characterises the interface. Computed tomography-based plain strain finite element analyses models of the cement-bone interface are built and loaded in multiple directions. Periodic boundaries are considered and the failure of the cement and bone fractions by cracking of the bulk components are included. The results compare favourably with experimental observations. Surprisingly, the analyses reveal that under shear loading no failure occurs and considerable normal compression is generated to prevent interface dilation. Reaction forces, crack patterns and stress fields provide more insight into the mixed-mode failure process. Moreover, the cement-bone interface analyses provide details which can serve as a basis for the development of a cohesive law.     

How plate positioning impacts the biomechanics of the open wedge tibial osteotomy; A finite element analysis

A numerical model of the medial open wedge tibial osteotomy based on the finite element method was developed. Two plate positions were tested numerically. In a configuration, (a), the plate was fixed in a medial position and (b) in an anteromedial position. The simulation took into account soft tissues preload, muscular tonus and maximal gait load.

The maximal stresses observed in the four structural elements (bone, plate, wedge, screws) of an osteotomy with plate in medial position were substantially higher (1.13–2.8 times more) than those observed in osteotomy with an anteromedial plate configuration. An important increase (1.71 times more) of the relative micromotions between the wedge and the bone was also observed. In order to avoid formation of fibrous tissue at the bone wedge interface, the osteotomy should be loaded under 18.8% (～50 kg) of the normal gait load until the osteotomy interfaces union is achieved.     

How plate positioning impacts the biomechanics of the open wedge tibial osteotomy; a finite element analysis

A numerical model of the medial open wedge tibial osteotomy based on the finite element method was developed. Two plate positions were tested numerically. In a configuration, (a), the plate was fixed in a medial position and (b) in an anteromedial position. The simulation took into account soft tissues preload, muscular tonus and maximal gait load.The maximal stresses observed in the four structural elements (bone, plate, wedge, screws) of an osteotomy with plate in medial position were substantially higher (1.13-2.8 times more) than those observed in osteotomy with an anteromedial plate configuration. An important increase (1.71 times more) of the relative micromotions between the wedge and the bone was also observed. In order to avoid formation of fibrous tissue at the bone wedge interface, the osteotomy should be loaded under 18.8% (approximately 50 kg) of the normal gait load until the osteotomy interfaces union is achieved.     

Effect of mechanical loading on heterotopic ossification in cervical total disc replacement: a three-dimensional finite element analysis

The development of heterotopic ossification (HO) is considered one of the major complications following cervical total disc replacement (TDR). Even though previous studies have identified clinical and biomechanical conditions that may stimulate HO, the mechanism of HO formation has not been fully elucidated. The objective of this study is to investigate whether mechanical loading is a biomechanical condition that plays a substantial role to decide the HO formation. A finite element model of TDR on the C5–C6 was developed, and HO formation was predicted by simulating a bone adaptation process under various physiological mechanical loadings. The distributions of strain energy on vertebrae were assessed after HO formation. For the compressive force, most of the HO formation occurred on the vertebral endplates uncovered by the implant footplate which was similar to the Type 1 HO. For the anteriorly directed shear force, the HO was predominantly formed in the anterior parts of both the upper and lower vertebrae as the Type 2 HO. For both the flexion and extension moments, the HO shapes were similar to those for the shear force. The total strain energy was reduced after HO formation for all loading conditions. Two distinct types of HO were predicted based on mechanically induced bone adaptation processes, and our findings were consistent with those of previous clinical studies. HO formation might have a role in compensating for the non-uniform strain energy distribution which is one of the mechanical parameters related to the bone remodeling after cervical TDR.     

Parameter identification of the human lower limb under dynamic, transient torsional loading

The response of the lower limb to dynamic, transient torsional loading applied at the foot has been measured for a male test subject. The dynamic loading was provided by a computer controlled pneumatic system which applied single haversine (i.e. half cycle of a sine wave) axial moment pulses of variable amplitude (0-100 Nm) and duration (50-600 ms). Potentiometers measured the absolute rotations of the three leg segments. Test variables included rotation direction, weight bearing and joint flexion. Two approaches were explored for specifying parameters (i.e. inertia, damping, stiffness) of a three degree-of-freedom dynamic system model which best duplicated the measured response. One approach involved identification of linear parameters by means of optimization while the other approach entailed estimation. Parameter estimates, which included non-linear, asymmetric stiffness functions, were derived from the literature. The optimization was undertaken so as to identify parameter dependence on test variables. Results indicate that parameter values are influenced by test variables. Results also indicate that the non-linear, estimated model better approximates the experimental data than the linear, identified model. In addition to identifying parameters of a three degree-of-freedom model, parameters were also identified for a single degree-of-freedom model where the motion variable was intended to indicate the rotation of the in vivo knee. It is concluded that the simpler model offers good accuracy in predicting both magnitude and time of occurrence of peak knee axial rotations. Model motion fails to track the measured knee rotation subsequent to the peak, however.