Biomechanical differences of Coflex-F and pedicle screw fixation combined with TLIF or ALIF--a finite element study

Lumbar interbody fusion is a common procedure for treating lower back pain related to degenerative disc diseases. The Coflex-F is a recently developed interspinous spacer, the makers of which claim that it can provide stabilisation similar to pedicle screw fixation. Therefore, this study compares the biomechanical behaviour of the Coflex-F device and pedicle screw fixation with transforaminal lumbar interbody fusion (TLIF) or anterior lumbar interbody fusion (ALIF) surgeries by using finite element analysis. The results show that the Coflex-F device combined with ALIF surgery can provide stability similar to the pedicle screw fixation combined with TLIF or ALIF surgery. Also, the posterior instrumentations (Coflex-F and pedicle screw fixation) combined with TLIF surgery had lower stability than when combined with ALIF surgery.     

Estimation of trunk muscle forces using the finite element method and in vivo loads measured by telemeterized internal spinal fixation devices.

Direct quantitative measurement of muscle forces is not possible. Forces in the trunk muscles were estimated for standing and flexion of the upper body using three-dimensional, nonlinear finite element models of the lumbar spine with and without an internal spinal fixation device. Muscle forces assumed were two pairs dorsally and one ventrally, each representing several muscles. Muscle forces in the model with internal fixators were varied in discrete steps until the implant loads calculated closely corresponded to those measured in a patient with an instrumented implant. The calculated angles between adjacent lumbar vertebrae were compared with corresponding values measured on X-ray films of a patient as well as with literature values and served as a second criterion for predicting muscle forces. For the model without an implant, the muscle forces of the first model were slightly varied until the lumbar spine shape and the intradiscal pressure were physiological. The abdomen was shown to have a considerable supporting function for flexion.     

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.     

Effects of lumbo-pelvic rhythm on trunk muscle forces and disc loads during forward flexion: A combined musculoskeletal and finite element simulation study

Previous in-vivo studies suggest that the ratio of total lumbar rotation over pelvic rotation (lumbo-pelvic rhythm) during trunk sagittal movement is essential to evaluate spinal loads and discriminate between low back pain and asymptomatic population. Similarly, there is also evidence that the lumbo-pelvic rhythm is key for evaluation of realistic muscle and joint reaction forces and moments predicted by various computational musculoskeletal models. This study investigated the effects of three lumbo-pelvic rhythms defined based on in-vivo measurements on the spinal response during moderate forward flexion (60°) using a combined approach of musculoskeletal modeling of the upper body and finite element model of the lumbosacral spine. The muscle forces and joint loads predicted by the musculoskeletal model, together with the gravitational forces, were applied to the finite element model to compute the disc force and moment, intradiscal pressure, annular fibers strain, and load-sharing. The results revealed that a rhythm with high pelvic rotation and low lumbar flexion involves more global muscles and increases the role of the disc in resisting spinal loads, while its counterpart, with low pelvic rotation, recruits more local muscles and engages the ligaments to lower the disc loads. On the other hand, a normal rhythm that has balanced pelvic and lumbar rotations yields almost equal disc and ligament load-sharing and results in more balanced synergy between global and local muscles. The lumbo-pelvic rhythm has less effect on the intradiscal pressure and annular fibers strain. This work demonstrated that the spinal response during forward flexion is highly dependent on the lumbo-pelvic rhythm. It is therefore, essential to adapt this parameter instead of using the default values in musculoskeletal models for accurate prediction of muscle forces and joint reaction forces and moments. The findings provided by this work are expected to improve knowledge of spinal response during forward flexion, and are clinically relevant towards low back pain treatment and disc injury prevention.     

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.     

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.     

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.     

Validation performance comparison for finite element models of the human brain

The objective of this study was to compare the performance of six validated brain finite element (FE) models to localized brain motion validation data in five experimental configurations. Model performance was measured using the objective metric CORA (CORrelation and Analysis), where higher ratings represent better correlation. The KTH model achieved the highest average CORA rating, and the ABM received the highest average rating among models robustly validated against five cadaver impacts in three directions. This technique can be more frequently employed to build better models and, when associated limitations are well understood, to compare inter-model performance under similar conditions.     

Experimental validation of finite element predicted bone strain in the human metatarsal

Metatarsal stress fracture is a common injury observed in athletes and military personnel. Mechanical fatigue is believed to play an important role in the etiology of stress fracture, which is highly dependent on the resulting bone strain from the applied load. The purpose of this study was to validate a subject-specific finite element (FE) modeling routine for bone strain prediction in the human metatarsal. Strain gauge measurements were performed on 33 metatarsals from seven human cadaveric feet subject to cantilever bending, and subject-specific FE models were generated from computed tomography images. Material properties for the FE models were assigned using a published density-modulus relationship as well as density-modulus relationships developed from optimization techniques. The optimized relationships were developed with a ‘training set’ of metatarsals (n = 17) and cross-validated with a ‘test set’ (n = 16). The published and optimized density elasticity equations provided FE-predicted strains that were highly correlated with experimental measurements for both the training (r² ≥ 0.95) and test (r² ≥ 0.94) sets; however, the optimized equations reduced the maximum error by 10% to 20% relative to the published equation, and resulted in an X = Y type of relationship between experimental measurements and FE predictions. Using a separate optimized density-modulus equation for trabecular and cortical bone did not improve strain predictions when compared to a single equation that spanned the entire bone density range. We believe that the FE models with optimized material property assignment have a level of accuracy necessary to investigate potential interventions to minimize metatarsal strain in an effort to prevent the occurrence of stress fracture.     

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 and finite element analysis for prediction of kidney injury under blunt impact

Kidneys are third most injured organs in abdominal trauma after liver and spleen; this study therefore is an attempt to understand the behaviour of kidneys under blunt trauma. Dynamic impact tests were performed on 20 fresh porcine kidneys to study the injury propagation in the organ, and the acceleration of the impactor was measured. A kidney model was developed with structural details like capsule and cortex. The kidney cortex was modelled with solid hexahedral elements and the capsule was modelled with quadratic shell elements. The material models for the capsule and cortex were used from the experimental data reported in our previous study. The developed model was calibrated using previous and current experimental results to reproduce the injuries of the organ in terms of acceleration of the impactor, and the injuries sustained by the organ during the experiments. The developed kidney model is observed to be robust and can be integrated with the available human body finite element models to simulate accidents and to predict or simulate injuries.     

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.     

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.     

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.     

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.