Optimisation of Intervertebral Disc Mechanical Properties and the Impact of Vertebral Alignment in Subject‐Specific Finite Element Models

• Published in: International journal for numerical methods in biomedical engineering Vol. 41; no. 6; pp. e70052 - n/a
• Main Authors: Kelly, Emily S., Javadi, Akbar A., Holsgrove, Timothy P., Ward, Michael, Williams, David, Williams, Jenny, Holt, Cathy, Meakin, Judith R.
• Format: Journal Article
• Language: English
• Published: Hoboken, USA John Wiley & Sons, Inc 01.06.2025; Wiley Subscription Services, Inc
• Subjects: Alignment; Animals; Axial compression; Basic Research; Biomechanical Phenomena; Boundary conditions; Finite Element Analysis; Finite element method; finite element modelling; Humans; In vitro methods and tests; Intervertebral Disc - physiology; intervertebral disc properties; Intervertebral discs; lumbar spine; Lumbar Vertebrae - physiology; Mathematical models; Mechanical properties; Models, Biological; Optimization; Segments; Spine; Spine (lumbar); Stiffness matrix; Stress, Mechanical; subject‐specific modelling; Swine; Vertebrae; Weight-Bearing; finite element modelling; intervertebral disc properties; lumbar spine; subject‐specific modelling
• ISSN: 2040-7939; 2040-7947; 2040-7947
• DOI: 10.1002/cnm.70052

ABSTRACT Subject‐specific finite element models could improve understanding of how spinal loading varies between people, based on differences in morphology and tissue properties. However, determining accurate subject‐specific intervertebral disc (IVD) properties can be difficult due to the spine's complex behaviour, in six degrees of freedom. Previous studies optimising IVD properties have utilised axial compression alone or range of motion data in three axes. This study aimed to optimise IVD properties using 6‐axis force‐moment data, and compare the resultant model's accuracy against a model optimised using IVD pressure data. Additionally, model vertebral alignment was assessed to determine if differences between imaged specimen alignment and in vitro 6‐axis test alignment affected the optimisation process. A finite element model of a porcine lumbar motion segment was developed, with generic IVD properties. The model loading and boundary conditions replicated in vitro 6‐axis stiffness matrix testing of the same specimen. The model was then optimised twice, once using experimental IVD pressures and once using forces and moments. A second model with geometry based on the specimen's vertebral alignment from the 6‐axis testing was also developed and optimised. The 6‐axis force‐moment optimised model had more accurate overall 6‐axis load‐displacement behaviour, but less accurate IVD pressures than the pressure optimised model. Neither optimised model fully captured spinal behaviours, due to model and optimisation process limitations. The 6‐axis vertebral alignment model had lower error and different optimised IVD properties than the imaged vertebral alignment model. Thus, vertebral alignment affected segment stiffness, so should be considered when developing spine models. Optimisation methods were compared for a porcine lumbar intervertebral disc model. Neither optimising with 6‐axis forces and moments nor optimising with disc pressures fully captured spinal behaviour. Addressing model limitations could improve the process. Flexion‐extension optimisation was more accurate with geometry based on experimental rather than imaged vertebral alignment.