Finite Element Lumbar Spine Facet Contact Parameter Predictions are Affected by the Cartilage Thickness Distribution and Initial Joint Gap Size

Abstract

Current finite element modeling techniques utilize geometrically inaccurate cartilage distribution representations in the lumbar spine. We hypothesize that this shortcoming severely limits the predictive fidelity of these simulations. Specifically, it is unclear how these anatomically inaccurate cartilage representations alter range of motion and facet contact predictions. In the current study, cadaveric vertebrae were serially sectioned, and images were taken of each slice in order to identify the osteochondral interface and the articulating surface. A series of custom-written algorithms were utilized in order to quantify each facet joint's three-dimensional cartilage distribution using a previously developed methodology. These vertebrae-dependent thickness cartilage distributions were implemented on an L1 through L5 lumbar spine finite element model. Moments were applied in three principal planes of motion, and range of motion and facet contact predictions from the variable thickness and constant thickness distribution models were determined. Initial facet gap thickness dimensions were also parameterized. The data indicate that the mean and maximum cartilage thickness increased inferiorly from L1 to L5, with an overall mean thickness value of 0.57 mm. Cartilage distribution and initial facet joint gap thickness had little influence on the lumbar range of motion in any direction, whereas the mean contact pressure, total contact force, and total contact area predictions were altered considerably. The data indicate that range of motion predictions alone are insufficient to establish model validation intended to predict mechanical contact parameters. These data also emphasize the need for the careful consideration of the initial facet joint gap thickness with respect to the spinal condition being studied.

In vivo loads in the lumbar L3–4 disc during a weight lifting extension

“…Great efforts have been made to improve FE modeling by including various ligaments, muscle activities and realistic facet interactions (Shirazi-Adl, 2006) (Shirazi-Adl et al, 2005) (Natarajan et al, 2006) (Rohlmann et al, 2006) (Goel et al, 1993) (Woldtvedt et al, 2011) (Shirazi-Adl, 1994b). In our study, the kinematics of the each disc endplate was obtained from the net balance of all external and internal forces and moments, such as from the ligaments, muscles and facet joints.…”

Experimental and Computational Approach Investigating Burst Fracture Augmentation Using PMMA and Calcium Phosphate Cements

“…The concept of using finite element (FE) analysis for modeling the spine is relatively new, but in recent years an increasing number of studies have adopted the technique to evaluate specific aspects of the complex biomechanics of the spine 30,38,50,65,66. In addition, there are examples of the application of FE modeling for the assessment of the vertebroplasty procedure on osteoporotic vertebrae using PMMA cement 2,11,25,47,53,59,61,62.…”

A Computational Model to Describe the Regional Interlamellar Shear of the Annulus Fibrosus

“…Finite element (FE) modeling is a useful tool for predicting physiological loading environments [13][14][15][16][17]. Predictions of relative deformations and local stresses within the lamellar and interlamellar regions resulting from physiologically relevant loading conditions could be used to guide experiments that link biological response to local deformation (mechanobiology).…”