Parametric convergence sensitivity and validation of a finite element model of the human lumbar spine

The primary objective of this study was to generate a finite element model of the human lumbar spine (L1–L5), verify mesh convergence for each tissue constituent and perform an extensive validation using both kinematic/kinetic and stress/strain data. Mesh refinement was accomplished via convergence of strain energy density (SED) predictions for each spinal tissue. The converged model was validated based on range of motion, intradiscal pressure, facet force transmission, anterolateral cortical bone strain and anterior longitudinal ligament deformation predictions. Changes in mesh resolution had the biggest impact on SED predictions under axial rotation loading. Nonlinearity of the moment-rotation curves was accurately simulated and the model predictions on the aforementioned parameters were in good agreement with experimental data. The validated and converged model will be utilised to study the effects of degeneration on the lumbar spine biomechanics, as well as to investigate the mechanical underpinning of the contemporary treatment strategies.     

Comparison of eight published static finite element models of the intact lumbar spine: Predictive power of models improves when combined together

Finite element (FE) model studies have made important contributions to our understanding of functional biomechanics of the lumbar spine. However, if a model is used to answer clinical and biomechanical questions over a certain population, their inherently large inter-subject variability has to be considered. Current FE model studies, however, generally account only for a single distinct spinal geometry with one set of material properties. This raises questions concerning their predictive power, their range of results and on their agreement with in vitro and in vivo values.     

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.     

Pedicle screw-based posterior dynamic stabilisation of the lumbar spine: in vitro cadaver investigation and a finite element study

Pedicle screw-based dynamic constructs either benefit from a dynamic (flexible) interconnecting rod or a dynamic (hinged) screw. Both types of systems have been reported in the literature. However, reports where the dynamic system is composed of two dynamic components, i.e. a dynamic (hinged) screw and a dynamic rod, are sparse. In this study, the biomechanical characteristics of a novel pedicle screw-based dynamic stabilisation system were investigated and compared with equivalent rigid and semi-rigid systems using in vitro testing and finite element modelling analysis. All stabilisation systems restored stability after decompression. A significant decrease in the range of motion was observed for the rigid system in all loadings. In the semi-rigid construct the range of motion was significantly less than the intact in extension, lateral bending and axial rotation loadings. There were no significant differences in motion between the intact spine and the spine treated with the dynamic system (P>0.05). The peak stress in screws was decreased when the stabilisation construct was equipped with dynamic rod and/or dynamic screws.     

The influence of different magnitudes and methods of applying preload on fusion and disc replacement constructs in the lumbar spine: a finite element analysis

In a finite element (FE) analysis of the lumbar spine, different preload application methods that are used in biomechanical studies may yield diverging results. To investigate how the biomechanical behaviour of a spinal implant is affected by the method of applying the preload, hybrid-controlled FE analysis was used to evaluate the biomechanical behaviour of the lumbar spine under different preload application methods. The FE models of anterior lumbar interbody fusion (ALIF) and artificial disc replacement (ADR) were tested under three different loading conditions: a 150 N pressure preload (PP) and 150 and 400 N follower loads (FLs). This study analysed the resulting range of motion (ROM), facet contact force (FCF), inlay contact pressure (ICP) and stress distribution of adjacent discs. The FE results indicated that the ROM of both surgical constructs was related to the preload application method and magnitude; differences in the ROM were within 7% for the ALIF model and 32% for the ADR model. Following the application of the FL and after increasing the FL magnitude, the FCF of the ADR model gradually increased, reaching 45% at the implanted level in torsion. The maximum ICP gradually decreased by 34.1% in torsion and 28.4% in lateral bending. This study concluded that the preload magnitude and application method affect the biomechanical behaviour of the lumbar spine. For the ADR, remarkable alteration was observed while increasing the FL magnitude, particularly in the ROM, FCF and ICP. However, for the ALIF, PP and FL methods had no remarkable alteration in terms of ROM and adjacent disc stress.     

Simple finite element models for use in the design of therapeutic footwear

Therapeutic footwear is frequently prescribed in cases of rheumatoid arthritis and diabetes to relieve or redistribute high plantar pressures in the region of the metatarsal heads. Few guidelines exist as to how these interventions should be designed and what effect such interventions actually have on the plantar pressure distribution. Finite element analysis has the potential to assist in the design process by refining a given intervention or identifying an optimal intervention without having to actually build and test each condition. However, complete and detailed foot models based on medical image segmentation have proven time consuming to build and computationally expensive to solve, hindering their utility in practice. Therefore, the goal of the current work was to determine if a simplified patient-specific model could be used to assist in the design of foot orthoses to reduce the plantar pressure in the metatarsal head region. The approach is illustrated by a case study of a diabetic patient experiencing high pressures and pain over the fifth metatarsal head. The simple foot model was initially calibrated by adjusting the individual loads on the metatarsals to approximate measured peak plantar pressure distributions in the barefoot condition to within 3%. This loading was used in various shod conditions to identify an effective orthosis. Model results for metatarsal pads were considerably higher than measured values but predictions for uniform surfaces were generally within 16% of measured values. The approach enabled virtual prototyping of the orthoses, identifying the most favorable approach to redistribute the patient’s plantar pressures.     

Three-dimensional finite element stress analysis of the dentate human mandible.

The biomechanical events which accompany functional loading of the human mandible are not fully understood. The techniques normally used to record them are highly invasive. Computer modelling offers a promising alternative approach in this regard, with the additional ability to predict regional stresses and strains in inaccessible locations. In this study, we built two three-dimensional finite element (FE) models of a human mandible reconstructed from tomographs of a dry dentate jaw. The first model was used for a complete mechanical characterization of physical events. It also provided comparative data for the second model, which had an increased vertical corpus depth. In both cases, boundary conditions included rigid restraints at the first right molar and endosteal cortical surfaces of the articular eminences of temporal bones. Groups of parallel multiple vectors simulated individual masticatory muscle loads. The models were solved for displacements, stresses, strains, and forces. The simulated muscle loads in the first model deformed the mandible helically upward and toward its right (working) side. The highest principal stresses occurred at the bite point, anterior aspects of the coronoid processes, symphyseal region, and right and left sides of the mandibular corpus. In general, the observed principal stresses and strains were highest on the periosteal cortical surface and alveolar bone. At the symphyseal region, maximum principal stresses and strains were highest on the lower lingual mandibular aspect, whereas minimum principal stresses and strains were highest on its upper labial side. Subcondylar principal strains and condylar forces were higher on the left (balancing or nonbiting) side than on the right mandibular side, with condylar forces more concentrated on the anteromedial aspect of the working-side condyle and on the central and lateral aspects of the left. When compared with in vivo strain data from macaques during comparable biting events, the predictive strain values from the first model were qualitatively similar. In the second model, the reduced tensile stress on the working-side, and decreased shear stress bilaterally, confirmed that lower stresses occurred on the lower mandibular border with increased jaw depth. Our results suggested that although the mandible behaved in a beam-like manner, its corpus acted more like a combination of open and closed cross sections due to the presence of tooth sockets, at least for the task modelled.(ABSTRACT TRUNCATED AT 400 WORDS)     

Biomechanical implications of lumbar spinal ligament transection

Many lumbar spine surgeries either intentionally or inadvertently damage or transect spinal ligaments. The purpose of this work was to quantify the previously unknown biomechanical consequences of isolated spinal ligament transection on the remaining spinal ligaments (stress transfer), vertebrae (bone remodelling stimulus) and intervertebral discs (disc pressure) of the lumbar spine. A finite element model of the full lumbar spine was developed and validated against experimental data and tested in the primary modes of spinal motion in the intact condition. Once a ligament was removed, stress increased in the remaining spinal ligaments and changes occurred in vertebral strain energy, but disc pressure remained similar. All major biomechanical changes occurred at the same spinal level as the transected ligament, with minor changes at adjacent levels. This work demonstrates that iatrogenic damage to spinal ligaments disturbs the load sharing within the spinal ligament network and may induce significant clinically relevant changes in the spinal motion segment.     

Uncertainty analysis of material properties and morphology parameters in numerical models regarding the motion of lumbar vertebral segments

Abstract

The kinematics of a spinal motion segment is determined by the material properties of the soft-tissue and the morphology. The material properties can vary within subjects and between vertebral levels, leading to a wide possible range of motion of a spinal segment independently on its morphology. The goal of this numerical study was to identify the most influential material parameters concerning the kinematics of a spinal motion segment and their plausible ranges. Then, a method was tested to deduce the material properties automatically, based on a given ROM and morphology. A fully parametric finite element model of the morphology and material properties of a lumbar spinal motion segment was developed. The impact of uncertainty of twelve spinal material parameters, as well as the size of the gap between the articular surfaces of the facet joints was examined. The simulation results were compared to our own in vitro data. The flexibility of a lumbar segment was especially influenced by the properties of the anterior annulus region, the facet gap size and the interspinous ligament. The high degree of uncertainty in the material properties and facet gap size published in the literature can lead to a wide scatter in the motion of a spinal segment, with a range of 6°-17° in the intact condition in flexion/extension, from 5°-22° in lateral bending and from 3°-14° in axial rotation. Statistical analysis of the variability might help to estimate the sensitivity and total uncertainty propagated through biomechanical simulations, affecting the reliability of the predictions.     

Parametric convergence sensitivity and validation of a finite element model of the human lumbar spine

The primary objective of this study was to generate a finite element model of the human lumbar spine (L1-L5), verify mesh convergence for each tissue constituent and perform an extensive validation using both kinematic/kinetic and stress/strain data. Mesh refinement was accomplished via convergence of strain energy density (SED) predictions for each spinal tissue. The converged model was validated based on range of motion, intradiscal pressure, facet force transmission, anterolateral cortical bone strain and anterior longitudinal ligament deformation predictions. Changes in mesh resolution had the biggest impact on SED predictions under axial rotation loading. Nonlinearity of the moment-rotation curves was accurately simulated and the model predictions on the aforementioned parameters were in good agreement with experimental data. The validated and converged model will be utilised to study the effects of degeneration on the lumbar spine biomechanics, as well as to investigate the mechanical underpinning of the contemporary treatment strategies.     

High-resolution finite element models with tissue strength asymmetry accurately predict failure of trabecular bone

The ability to predict trabecular failure using microstructure-based computational models would greatly facilitate study of trabecular structure–function relations, multiaxial strength, and tissue remodeling. We hypothesized that high-resolution finite element models of trabecular bone that include cortical-like strength asymmetry at the tissue level, could predict apparent level failure of trabecular bone for multiple loading modes. A bilinear constitutive model with asymmetric tissue yield strains in tension and compression was applied to simulate failure in high-resolution finite element models of seven bovine tibial specimens. Tissue modulus was reduced by 95% when tissue principal strains exceeded the tissue yield strains. Linear models were first calibrated for effective tissue modulus against specimen-specific experimental measures of apparent modulus, producing effective tissue moduli of (mean±S.D.) 18.7±3.4 GPa. Next, a parameter study was performed on a single specimen to estimate the tissue level tensile and compressive yield strains. These values, 0.60% strain in tension and 1.01% strain in compression, were then used in non-linear analyses of all seven specimens to predict failure for apparent tensile, compressive, and shear loading. When compared to apparent yield properties previously measured for the same type of bone, the model predictions of both the stresses and strains at failure were not statistically different for any loading case (p>0.15). Use of symmetric tissue strengths could not match the experimental data. These findings establish that, once effective tissue modulus is calibrated and uniform but asymmetric tissue failure strains are used, the resulting models can capture the apparent strength behavior to an outstanding level of accuracy. As such, these computational models have reached a level of fidelity that qualifies them as surrogates for destructive mechanical testing of real specimens.     

Density-based load estimation using two-dimensional finite element models: a parametric study

A parametric investigation was conducted to determine the effects on the load estimation method of varying: (1) the thickness of back-plates used in the two-dimensional finite element models of long bones, (2) the number of columns of nodes in the outer medial and lateral sections of the diaphysis to which the back-plate multipoint constraints are applied and (3) the region of bone used in the optimization procedure of the density-based load estimation technique. The study is performed using two-dimensional finite element models of the proximal femora of a chimpanzee, gorilla, lion and grizzly bear. It is shown that the density-based load estimation can be made more efficient and accurate by restricting the stimulus optimization region to the metaphysis/epiphysis. In addition, a simple method, based on the variation of diaphyseal cortical thickness, is developed for assigning the thickness to the back-plate. It is also shown that the number of columns of nodes used as multipoint constraints does not have a significant effect on the method.     

Effect of lumbar fasciae on the stability of the lower lumbar spine

The biomechanical effect of tensioning the lumbar fasciae (LF) on the stability of the spine during sagittal plane motion was analysed using a validated finite element model of the normal lumbosacral spine (L4-S1). To apply the tension in the LF along the direction of the fibres, a local coordinate was allocated using dummy rigid beam elements that originated from the spinous process. Up to 10 Nm of flexion and 7.5 Nm of extension moment was applied with and without 20 N of lateral tension in the LF. A follower load of 400 N was additionally applied along the curvature of the spine. To identify how the magnitude of LF tension related to the stability of the spine, the tensioning on the fasciae was increased up to 40 N with an interval of 10 N under 7.5 Nm of flexion/extension moment. A fascial tension of 20 N produced a 59% decrease in angular motion at 2.5 Nm of flexion moment while there was a 12.3% decrease at 10 Nm in the L5-S1 segment. Its decrement was 53 and 9.6% at 2.5 Nm and 10 Nm, respectively, in the L4-L5 segment. Anterior translation was reduced by 12.1 and 39.0% at the L4-L5 and L5-S1 segments under 10 Nm of flexion moment, respectively. The flexion stiffness shows an almost linear increment with the increase in fascial tension. The results of this study showed that the effect of the LF on the stability of the spine is significant.     

An upper bound computational model for investigation of fusion effects on adjacent segment biomechanics of the lumbar spine

Abstract

Prediction of the biomechanical effects of fusion surgery on adjacent segments is a challenge in computational biomechanics of the spine. In this study, a two-segment L3-L4-L5 computational model was developed to simulate the effects of spinal fusion on adjacent segment biomechanical responses under a follower load condition. The interaction between the degenerative segment (L4-5) and the adjacent segment (L3-4) was simulated using an equivalent follower spring. The spring stiffness was calibrated using a rigid fusion of a completely degenerated disc model at the L4-5 level, resulting in an upper bound response at the adjacent (L3-4) segment. The obtained upper bound equivalent follower spring was used to simulate the upper bound biomechanical responses of fusion of the disc with different degeneration grades. It was predicted that as the disc degeneration grade at the degenerative segment decreased, the effect on the adjacent segment responses decreased accordingly after fusion. The data indicated that the upper bound computational model can be a useful computational tool for evaluation of the interaction between segments and for investigation of the biomechanical mechanisms of adjacent segment degeneration after fusion.     

A multibody modelling approach to determine load sharing between passive elements of the lumbar spine

The human spinal segment is an inherently complex structure, a combination of flexible and semi-rigid articulating elements stabilised by seven principal ligaments. An understanding of how mechanical loading is shared among these passive elements of the segment is required to estimate tissue failure stresses. A 3D rigid body model of the complete lumbar spine has been developed to facilitate the prediction of load sharing across the passive elements. In contrast to previous multibody models, this model includes a non-linear, six degrees of freedom intervertebral disc, facet bony articulations and all spinal ligaments. Predictions of segmental kinematics and facet joint forces, in response to pure moment loading (flexion–extension), were compared to published in vitro data. On inclusion of detailed representation of the disc and facets, the multibody model fully captures the non-linear flexibility response of the spinal segment, i.e. coupled motions and a mobile instantaneous centre of rotation. Predicted facet joint forces corresponded well with reported values. For the loading case considered, the model predicted that the ligaments are the main stabilising elements within the physiological motion range; however, the disc resists a greater proportion of the applied load as the spine is fully flexed. In extension, the facets and capsular ligaments provide the principal resistance. Overall patterns of load distribution to the spinal ligaments are in agreement with previous predictions; however, the current model highlights the important role of the intraspinous ligament in flexion and the potentially high risk of failure. Several important refinements to the multibody modelling of the passive elements of the spine have been described, and such an enhanced passive model can be easily integrated into a full musculoskeletal model for the prediction of spinal loading for a variety of daily activities.     

A finite element study relating to the rapid filling phase of the human ventricles

During the rapid diastolic filling phase at rest, the ventricles of the human heart double approximately in volume. In order to investigate whether the ventricular filling pressures measured under physiological conditions can give rise to such an extensive augmentation in ventricular volumes, a finite element model of the human right and left ventricles has been developed, taking into account the nonlinear mechanical behavior and effective compressibility of the myocardial tissue. The results were compared with the filling phase of the human left ventricle as extrapolated from measurements documented in the literature. We arrived at the conclusion that the ventricular pressures measured during the rapid filling phase cannot be the sole cause of the rise of the observed ventricular volumes. We rather advocate the assumption that further dilating mechanisms might be part of ventricular activity thus heralding a multiple function of the ventricular muscle body. A further result indicates that under normal conditions the influence of the viscoelasticity of the tissue should not be disregarded in ventricular mechanics.     

Using 3D finite element models verified the importance of callus material and microstructure in biomechanics restoration during bone defect repair

Background: There is lack of further observations on the microstructure and material property of callus during bone defect healing and the relationships between callus properties and the mechanical strength. Methods: Femur bone defect model was created in rabbits and harvested CT data to reconstruct finite element models at 1 and 2 months. Three types of assumed finite element models were compared to study the callus properties, which assumed the material elastic property as heterogeneous (R-model), homogenous (H-model) or did not change from 1 to 2 months (U-model). Results: The apparent elastic moduli increased at 2 months (from 355.58 ± 132.67 to 1139.30 ± 967.43 MPa) in R-models. But there was no significant difference in apparent elastic moduli between R-models (355.58 ± 132.67 and 1139.30 ± 967.43 MPa) and H-models (344.79 ± 138.73 and 1001.52 ± 692.12 MPa) in 1 and 2 months. A significant difference of apparent elastic moduli was found between the R-model (1139.30 ± 967.43 MPa) and U-model group (207.15 ± 64.60 MPa) in 2 months. Conclusions: This study showed that the callus structure stability remodeled overtime to achieve a more effective structure, while the material quality of callus tissue is a very important factor for callus strength. At the meantime, this study showed an evidence that the material heterogeneity maybe not as important as it is in bone fracture model.     

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

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

Discretization error when using finite element models: Analysis and evaluation of an underestimated problem

Mesh convergence tests are often insufficiently performed in finite element analyses. There are many parameters which may have an effect on the mesh convergence behavior. The aim of this study was to identify the influence of different parameters on the mesh convergence behavior.For this purpose we used a simplified axis-symmetrical model of a single pedicle screw flank with surrounding bone to simulate a pull-out test. In parameter studies, the flank radii and the contact conditions at the bone–screw interface were varied. These parameter studies were carried out using an implicit and explicit solver. Thereby, the convergence criteria and the number of the substeps for the implicit nonlinear iteration process as well as the velocity and the material density for the explicit approach were considered.The mesh convergence behavior was influenced by varying the flank radii and the contact conditions. The implicit calculations led to a reaction force, which converged rapidly to a certain value with increasing mesh density, whereas the maximum von-Mises stress showed substantial convergence problems. The number of substeps and the convergence criteria of the iteration process strongly influenced the implicit solutions. In contrast, the maximum von-Mises stresses resulting from explicit calculations converged to a certain value after only a few refinement steps. Different pull-out velocities substantially affected the mesh convergence behavior, while the material density showed only a negligible influence.The results indicated the need to perform an appropriate mesh convergence test when using finite element methods. We were able to show that different parameters strongly influence the mesh convergence behavior and we demonstrated that convergence tests do not always lead to a satisfactory or acceptable solution.