FEM-Based Compression Fracture Risk Assessment in Osteoporotic Lumbar Vertebra L1

Investigation of the Role of Osteoporotic Vertebra Degeneration on the Stability of the Lumbar Spine: In Silico Modelling under Compressive Loading

An evaluation of the impact of osteoporosis on loss of spinal stability, with or without intervertebral disc degeneration, using computational analysis is presented. The research also investigates the correlation between osteoporosis and intervertebral disc degeneration. Three-dimensional finite element models of human lumbar spine segments were used to assess the influence of osteoporosis on spinal stability. Five different models of age-related degeneration were created using various material properties for trabecular bone and intervertebral discs. Calculation results indicate that in a spine with osteoporosis, the deformation of the intervertebral discs can increase by more than 30% when compared to a healthy spine. Thus, intervertebral disc deformation depends not only on the degree of degeneration of the discs themselves, but their deformation is also influenced by the degree of osteoporosis of the vertebrae. Additionally, the load-bearing capacity of the spine can decrease by up to 30% with osteoporosis, regardless of the degree of intervertebral disc deformation. In conclusion, osteoporosis can contribute to intervertebral disc degeneration.

Modeling the Effect of Annulus Fibrosus Stiffness on the Stressed State of a Vertebral L1 Body and Nucleus Pulposus

The investigation examines the transference of stiffness from intervertebral discs (IVDs) to the lumbar body of the L1 vertebra and the interactions among adjacent tissues. A computational model of the vertebra was developed, considering parameters such as cortical bone thickness, trabecular bone elasticity, and the nonlinear response of the nucleus pulposus to external loading. A nonlinear dynamic analysis was performed, revealing certain trends: a heightened stiffness of the annulus fibrosus correlates with a significant reduction in the vertebral body's ability to withstand external loading. At a supplied displacement of 6 mm, the vertebra with a degenerative disc reached its yielding point, whereas the vertebrae with a healthy annulus fibrosus exhibited a strength capacity exceeding 20%. The obtained findings and proposed methodology are potentially useful for biomedical engineers and clinical specialists in evaluating the condition of the annulus fibrosus and predicting its influence on the bone components of the spinal system.

Bone Density Micro-CT Assessment during Embedding of the Innovative Multi-Spiked Connecting Scaffold in Periarticular Bone to Elaborate a Validated Numerical Model for Designing Biomimetic Fixation of Resurfacing Endoprostheses

Our team has been working for some time on designing a new kind of biomimetic fixation of resurfacing endoprostheses, in which the innovative multi-spiked connecting scaffold (MSC-Scaffold) that mimics the natural interface between articular cartilage and periarticular trabecular bone in human joints is the crucial element. This work aimed to develop a numerical model enabling the design of the considered joint replacement implant that would reflect the mechanics of interacting biomaterials. Thus, quantitative micro-CT analysis of density distribution in bone material during the embedding of MSC-Scaffold in periarticular bone was applied. The performed numerical studies and corresponding mechanical tests revealed, under the embedded MSC-Scaffold, the bone material densification affecting its mechanical properties. On the basis of these findings, the built numerical model was modified by applying a simulated insert of densified bone material. This modification led to a strong correlation between the re-simulation and experimental results (FVU = 0.02). The biomimetism of the MSC-Scaffold prototype that provided physiological load transfer from implant to bone was confirmed based on the Huber–von Mises–Hencky (HMH) stress maps obtained with the validated finite element (FE) model of the problem. The micro-CT bone density assessment performed during the embedding of the MSC-Scaffold prototype in periarticular bone provides insight into the mechanical behaviour of the investigated implant-bone system and validates the numerical model that can be used for the design of material and geometric features of a new kind of resurfacing endoprostheses fixation.

Effect of pedicle screw angles on the fracture risk of the human vertebra: A patient-specific computational model

The assessment of a human vertebra's stability after a screws fixation procedure and its fracture risk is still an open clinical problem. The accurate evaluation of fracture risk requires that all fracture mechanical determinants such as geometry, constitutive behavior, loading modes, and screws angulation are accounted for, which requires biomechanics-based analyses. As such, in the present work we investigate the effect of pedicle screws angulation on a patient-specific model of non osteoporotic lumbar vertebra, derived from clinical CT images. We propose a novel computational approach of fracture analysis and compare the effects of fixation stability in the lumbar spine. We considered a CT-based three-dimensional FE model of bilaterally instrumented L4 vertebra virtually implanting pedicle screws according to clinical guidelines. Nine screws trajectories were selected combining three craniocaudal and mediolateral angles, thus investigated through a parametric computational analysis. Bone was modeled as an elastic material with element-wise inhomogeneous properties fine-tuned on CT data. We implemented a custom algorithm to identify the thin cortical layer correctly from CT images ensuring reliable material properties in the computational model. Physiological motion (i.e., flexion, extension, axial rotation, lateral bending) was further accomplished by simultaneously loading the vertebra and the implant. We simulated local progressive damage of the bone by using a quasi-static force-driven incremental approach and considering a stress-based fracture criterion. Ductile-like and brittle-like fractures were found. Statistical analyses show significant differences comparing screws trajectories and averaging the results among six loading modes. In particular, we identified the caudomedial trajectory as the least critical case, thus safer from a clinical perspective. Instead, medial and craniolaterally oriented screws entailed higher peak and average stresses, though no statistical evidence classified such loads as the most critical scenarios. A quantitative validation procedure will be required in the future to translate our findings into clinical practice. Besides, to apply the results to the target osteoporotic population, new studies will be needed, including a specimen from an osteoporotic patient and the effect of osteoporosis on the constitutive model of bone.