A patient-specific finite element methodology to predict damage accumulation in vertebral bodies under axial compression, sagittal flexion and combined loads

Due to the inherent limitations of DXA, assessment of the biomechanical properties of vertebral bodies relies increasingly on CT-based finite element (FE) models, but these often use simplistic material behaviour and/or single loading cases. In this study, we applied a novel constitutive law for bone elasticity, plasticity and damage to FE models created from coarsened pQCT images of human vertebrae, and compared vertebral stiffness, strength and damage accumulation for axial compression, anterior flexion and a combination of these two cases. FE axial stiffness and strength correlated with experiments and were linearly related to flexion properties. In all loading modes, damage localised preferentially in the trabecular compartment. Damage for the combined loading was higher than cumulated damage produced by individual compression and flexion. In conclusion, this FE method predicts stiffness and strength of vertebral bodies from CT images with clinical resolution and provides insight into damage accumulation in various loading modes.     

Biomechanical comparison of conventional and optimised locking plates for the fixation of intraarticular calcaneal fractures: a finite element analysis

Intraarticular calcaneal fractures can result in poor prognosis. Although operative fixation can improve the functional outcomes in most cases, surgical complications such as loss of reduction and wound healing problems may increase the risk of reoperation. Hence, this study aimed to design calcaneal locking plate with a lower profile and better biomechanical performance and to compare the redesigned plate with the traditional calcaneal plate via the finite element method. A Sanders’ type II-C intraarticular calcaneal fracture was simulated. Two fixation models utilising the branch-like calcaneal locking plate and the full plate were constructed. Topology optimisation was conducted to generate a new calcaneal plate design. A biomechanical comparison among the three groups of plates was performed using the finite element method. For the fracture simulated in this study, the optimised plate was superior to the traditional plate in terms of fixation stability and safety but was reduced in volume by approximately 12.34%. In addition, more rational stress distributions were observed in the redesigned plate, underscoring the superiority of this new design in terms of fatigue strength. These results demonstrate that the topology optimisation can be used to design a new implant with a minimised profile and no loss of fixation stability.     

Off-axis loads cause failure of the distal radius at lower magnitudes than axial loads: a finite element analysis

Distal radius (Colles') fractures are a common fall-related injury in older adults and frequently result in long-term pain and reduced ability to perform activities of daily living. Because the occurrence of a fracture during a fall depends on both the strength of the bone and upon the kinematics and kinetics of the impact itself, we sought to understand how changes in bone mineral density (BMD) and loading direction affect the fracture strength and fracture initiation location in the distal radius. A three-dimensional finite element model of the radius, scaphoid, and lunate was used to examine changes of +/-2% and +/-4% BMD, and both axial and physiologically relevant off-axis loads on the radius. Changes in BMD resulted in similar percent changes in fracture strength. However, modifying the applied load to include dorsal and lateral components (assuming a dorsal view of the wrist, rather than an anatomic view) resulted in a 47% decrease in fracture strength (axial failure load: 2752N, off-axis: 1448N). Loading direction also influenced the fracture initiation site. Axially loaded radii failed on the medial surface immediately proximal to the styloid process. In contrast, off-axis loads, containing dorsal and lateral components, caused failure on the dorsal-lateral surface. Because the radius appears to be very sensitive to loading direction, the results suggest that much of the variability in fracture strength seen in cadaver studies may be attributed to varying boundary conditions. The results further suggest that interventions focused on reducing the incidence of Colles' fractures when falls onto the upper extremities are unavoidable may benefit from increasing the extent to which the radius is loaded along its axis.     

Experimental validation of intact and implanted distal femur finite element models

Four finite element (FE) models of intact and distal femur of knee replacements were validated relative to measured bone strains. FE models of linear tetrahedrons were used. Femoral replacements with cemented stemless, cemented and noncemented femoral stems of the PFC Sigma Modular Knee System were analyzed. Bone strains were recorded at ten locations on the cortex. The magnitude of the FE bone strains corresponded to the mean measured strains, with an overall agreement of 10%. Linear regression between the FE and mean experimental strains produced slopes between 0.94 and 1.06 and R(2) values between 0.92 and 0.99. RSME values were less than 12%. The FE models were able to adequately replicate the mechanical behavior of distal femur reconstructions.     

Characterization of human passive muscles for impact loads using genetic algorithm and inverse finite element methods

The objective of this study is to identify the dynamic material properties of human passive muscle tissues for the strain rates relevant to automobile crashes. A novel methodology involving genetic algorithm (GA) and finite element method is implemented to estimate the material parameters by inverse mapping the impact test data. Isolated unconfined impact tests for average strain rates ranging from 136 s⁻¹ to 262 s⁻¹ are performed on muscle tissues. Passive muscle tissues are modelled as isotropic, linear and viscoelastic material using three-element Zener model available in PAMCRASH^TM explicit finite element software. In the GA based identification process, fitness values are calculated by comparing the estimated finite element forces with the measured experimental forces. Linear viscoelastic material parameters (bulk modulus, short term shear modulus and long term shear modulus) are thus identified at strain rates 136 s⁻¹, 183 s⁻¹ and 262 s⁻¹ for modelling muscles. Extracted optimal parameters from this study are comparable with reported parameters in literature. Bulk modulus and short term shear modulus are found to be more influential in predicting the stress-strain response than long term shear modulus for the considered strain rates. Variations within the set of parameters identified at different strain rates indicate the need for new or improved material model, which is capable of capturing the strain rate dependency of passive muscle response with single set of material parameters for wide range of strain rates.     

Cortical bone viscoelasticity and fixation strength of press-fit femoral stems: finite element model

Many cementless implant designs rely upon a diaphyseal press-fit in conjunction with a porous coated implant surface to achieve primary or short term fixation, thereby constraining interface micromotion to such a level that bone ingrowth and consequent secondary or long-term fixation, i.e., osseointegration, can occur. Bone viscoelasticity, however, has been found to affect stem primary stability by reducing push-out load. In this investigation, an axisymmetric finite element model of a cylindrical stem and diaphyseal cortical bone section was created in order to parametrically evaluate the effect of bone viscoelasticity on stem push-out while controlling coefficient of friction (mu = 0.15, 0.40, and 1.00) and stem-bone diametral interference (delta = 0.01, 0.05, 0.10, and 0.50 mm). Based on results from a previous study, it was hypothesized that stem-bone interference (i.e., press-fit) would elicit a bone viscoelastic response which would reduce the initial fixation of the stem as measured by push-out load. Results indicate that for all examined combinations of mu and delta, bone viscoelastic behavior reduced the push-out load by a range of 2.6-82.6% due to stress relaxation of the bone. It was found that the push-out load increased with mu for each value of delta, but minimal increases in the push-out load (2.9-4.9%) were observed as delta was increased beyond 0.10 mm. Within the range of variables reported for this study, it was concluded that bone viscoelastic behavior, namely stress relaxation, has an asymptotic affect on stem contact pressure, which reduces stem push-out load. It was also found that higher levels of coefficient of friction are beneficial to primary fixation, and that an interference "threshold" exists beyond which no additional gains in push-out load are achieved.     

Determination of cellular strains by combined atomic force microscopy and finite element modeling

Many organs adapt to their mechanical environment as a result of physiological change or disease. Cells are both the detectors and effectors of this process. Though many studies have been performed in vitro to investigate the mechanisms of detection and adaptation to mechanical strains, the cellular strains remain unknown and results from different stimulation techniques cannot be compared. By combining experimental determination of cell profiles and elasticities by atomic force microscopy with finite element modeling and computational fluid dynamics, we report the cellular strain distributions exerted by common whole-cell straining techniques and from micromanipulation techniques, hence enabling their comparison. Using data from our own analyses and experiments performed by others, we examine the threshold of activation for different signal transduction processes and the strain components that they may detect. We show that modulating cell elasticity, by increasing the F-actin content of the cytoskeleton, or cellular Poisson ratio are good strategies to resist fluid shear or hydrostatic pressure. We report that stray fluid flow in some substrate-stretch systems elicits significant cellular strains. In conclusion, this technique shows promise in furthering our understanding of the interplay among mechanical forces, strain detection, gene expression, and cellular adaptation in physiology and disease.     

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.     

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.     

Porcine models in spinal research: calibration and comparative finite element analysis of various configurations during flexion-extension

This study was conducted to develop and calibrate a detailed 3-dimensional finite element model of the porcine lumbar spine and to compare this model with various configurations in flexion and extension. Computed tomography scans obtained from the L4-L5 lumbar segment of a Landrace x Large White pig were used to generate a solid volume. The various passive components were characterized by using a step-by-step calibration procedure in which the material properties of the anatomic structures were modified to match the corresponding in vitro data set-points retrieved from the literature. The range of motion of the totally assembled intact model was assessed under a 10-Nm flexion-extension moment and compared with data from a bilateral complete and hemifacetectomy configuration. In addition, the results from our porcine model were compared with published data regarding range of motion in a human finite element model in order to predict the configuration of the porcine model that most closely represented the human spine. Both the intact and hemifacetectomy configurations of the porcine model were comparable to the human spine. However, qualitative analysis of the instantaneous axis of rotation revealed a dissimilarity between the intact porcine model and human spine behavior, indicating the hemifacetectomy configuration of the porcine model as the most appropriate for spinal instrumentation studies. The present 3-dimensional finite element porcine model offers an additional tool to improve understanding of the biomechanics of the porcine spine and to decrease the expense of spinal research.     

Explicit finite element modelling of heel pad mechanics in running: inclusion of body dynamics and application of physiological impact loads

Many heel pathologies including plantar heel pain may result from micro tears/trauma in the subcutaneous tissues, in which internal tissue deformation/stresses within the heel pad play an important role. Previously, many finite element models have been proposed to evaluate stresses inside the heel pad, but the majority of these models only focus on static loading boundary conditions. This study explored a dynamics modelling approach to the heel pad subjected to realistic impact loads during running. In this model, the inertial property and action of the body are described by a lumped parameter model, while the heel/shoe interactions are modelled using a viscoelastic heel pad model with contact properties. The impact force pattern, dynamic heel pad deformation and stress states predicted by the model were compared with published experimental data. Further parametrical studies revealed the model responses, in terms of internal stresses in the skin and fatty tissue, change nonlinearly when body dynamics changes. A reduction in foot's touchdown velocity resulted in a less severe impact landing and stress relief inside the heel pad, for example peak von-Mises stress in fatty tissue, was reduced by 11.3%. Applications of the model may be extendable to perform iterative analyses to further understand the complex relationships between body dynamics and stress distributions in the soft tissue of heel pad during running. This may open new opportunities to study the mechanical aetiology of plantar heel pain in runners.     

Biomechanical analysis and comparison of 12 dental implant systems using 3D finite element study

Finite element analysis plays an important role in dental implant design. The objective of this study was to show the effect of the overall geometry of dental implants on their biomechanics after implantation. In this study, 12 dental implants, with the same length, diameter and screw design, were simulated from different implant systems. Numerical model of right mandibular incisor bone segment was generated from CT data. The von-Mises stress distributions and the total deformation distributions under vertical/lateral load were compared for each implant by scores ranking method. The implants with cylindrical shapes had highest scores. Results indicated that cylindrical shape represented better geometry over taper implant. This study is helpful in choosing the optimal dental implant for clinical application and also contributes to individual implant design. Our study could also provide reference for choice and modification of dental implant in any other insertion sites and bone qualities.     

A finite element model technique to determine the mechanical response of a lumbar spine segment under complex loads

This study presents a CT-based finite element model of the lumbar spine taking into account all function-related boundary conditions, such as anisotropy of mechanical properties, ligaments, contact elements, mesh size, etc. Through advanced mesh generation and employment of compound elements, the developed model is capable of assessing the mechanical response of the examined spine segment for complex loading conditions, thus providing valuable insight on stress development within the model and allowing the prediction of critical loading scenarios. The model was validated through a comparison of the calculated force-induced inclination/deformation and a correlation of these data to experimental values. The mechanical response of the examined functional spine segment was evaluated, and the effect of the loading scenario determined for both vertebral bodies as well as the connecting intervertebral disc.     

An efficient and accurate prediction of the stability of percutaneous fixation of acetabular fractures with finite element simulation

Posterior wall fracture is one of the most common fracture types of the acetabulum and a conventional approach is to perform open reduction and internal fixation with a plate and screws. Percutaneous screw fixations, on the other hand, have recently gained attention due to their benefits such as less exposure and minimization of blood loss. However their biomechanical stability, especially in terms interfragmentary movement, has not been investigated thoroughly. The aims of this study are twofold: (1) to measure the interfragmentary movements in the conventional open approach with plate fixations and the percutaneous screw fixations in the acetabular fractures and compare them; and (2) to develop and validate a fast and efficient way of predicting the interfragmentary movement in percutaneous fixation of posterior wall fractures of the acetabulum using a 3D finite element (FE) model of the pelvis. Our results indicate that in single fragment fractures of the posterior wall of the acetabulum, plate fixations give superior stability to screw fixations. However screw fixations also give reasonable stability as the average gap between fragment and the bone remained less than 1 mm when the maximum load was applied. Our finite element model predicted the stability of screw fixation with good accuracy. Moreover, when the screw positions were optimized, the stability predicted by our FE model was comparable to the stability obtained by plate fixations. Our study has shown that FE modeling can be useful in examining biomechanical stability of osteosynthesis and can potentially be used in surgical planning of osteosynthesis.     

Configuration of anchorage holes affects fixation of the acetabular component in cemented total hip replacement--a finite element study

Our survey of current practice among UK orthopaedic surgeons shows wide variations in fixation techniques. The aim of this study, is to investigate the effect of drilling different configurations of anchorage holes in the acetabulum on implant stability. To avoid variables that could incur during in vitro testing, we used commercially available COSMOS finite element analysis package to investigate the stress distributions, deformations, and strains on the cement mantle when drilling three large anchorage holes and six smaller ones, with straight and rounded cement pegs. The results, which are in line with our in vitro studies on simulated reconstructed acetabulae, indicate better stability of the acetabular component when three larger holes than six smaller holes are drilled and when the necks of the anchorage holes are rounded. The longevity of total hip replacements could be improved by drilling three large anchorage holes, rather than many smaller ones, as initially proposed by Charnley.     

Statistical methods in finite element analysis

Finite element analysis (FEA) is a commonly used tool within many areas of engineering and can provide useful information in structural analysis of mechanical systems. However, most analyses within the field of biomechanics usually take no account either of the wide variation in material properties and geometry that may occur in natural tissues or manufacturing imperfections in synthetic materials. This paper discusses two different methods of incorporating uncertainty in FE models. The first, Taguchi's robust parameter design, uses orthogonal matrices to determine how to vary the parameters in a series of FE models, and provides information on the sensitivity of a model to input parameters. The second, probabilistic analysis, enables the distribution of a response variable to be determined from the distributions of the input variables. The methods are demonstrated using a simple example of an FE model of a beam that is assigned material properties and geometry over a range similar to an orthopaedic fixation plate. In addition to showing how each method may be used on its own, we also show how computational effort may be minimised by first identifying the most important input variables before determining the effects of imprecision.