Authors
Abstract
Background: About 80% of the population will experience back pain in their lifetime; however, although many patients have low back pain associated with disc degeneration, the exact course of degeneration is still unclear. The disc degeneration disorder has affected one-third of the world's young population. During degeneration, the disc undergoes morphological and biochemical changes, which in turn alter the tissue hydration, permeability, and ultimately the load-bearing capacity of the disc. Therefore, the finite element model, designed to study the relationship between frequent loading and disk degeneration, must be able to analyze the complex loading in the in-vivo conditions. The aim of this study was to construct and update models of finite element components with the prolastic properties so that different quasi-static loads could be investigated by presenting a personalized model in different individuals and applied in clinical studies to simulate the daily biomechanical behavior for accurate diagnosis and treatment.
Method: This study simulated three different modes of finite element modeling, including axial symmetry method, parametric model and precision model with poroelastic mechanical properties and its results were compared with experimental in-vivo experiments.
Results: To validate the constructed models, the results of three different quasi-static creep experiments were performed, including short-term creep, long-term creep and creep under regular daily activities, the results of which predicted changes. The results predicted height changes, axial displacement of the spine and the intradiscal pressure of the nucleus.
Conclusion: All the proposed results indicated that the models presented in quasi-statistic behavior predicted acceptable results and have sufficient validity to be examined in other quasi-statistic experiments. Therefore, it is possible to take a step forward in examining the results of clinical activities in determining the process of intervertebral disc degeneration.
- Chagnon A, Aubin CE, Villemure I. Biomechanical influence of disk properties on the load transfer of healthy and degenerated disks using a poroelastic finite element model. J Biomech Eng. 2010;132(11):111006.
- Nikkhoo M, Wang JL, Parnianpour M, El-Rich M, Khalaf K. Biomechanical response of intact, degenerated and repaired intervertebral discs under impact loading - Ex-vivo and In-Silico investigation. J Biomech. 2018;70:26-32.
- Williams JR, Natarajan RN, Andersson GB. Inclusion of regional poroelastic material properties better predicts biomechanical behavior of lumbar discs subjected to dynamic loading. J Biomech. 2007;40(9):1981-7.
- Natarajan RN, Williams JR, Lavender SA, Andersson GBJ. Poro-elastic finite element model to predict the failure progression in a lumbar disc due to cyclic loading. Comput Struct. 2007;85(11-14):1142-51.
- Galbusera F, Schmidt H, Noailly J, Malandrino A, Lacroix D, Wilke HJ, et al. Comparison of four methods to simulate swelling in poroelastic finite element models of intervertebral discs. J Mech Behav Biomed Mater. 2011;4(7):1234-41.
- Jacobs NT, Cortes DH, Peloquin JM, Vresilovic EJ, Elliott DM. Validation and application of an intervertebral disc finite element model utilizing independently constructed tissue-level constitutive formulations that are nonlinear, anisotropic, and time-dependent. J Biomech. 2014;47(11):2540-6.
- Abaqus Unified FEA. Dassault Systèmes; 2019 [Available from: https://www.3ds.com/products-services/simulia/products/abaqus/.
- Panjabi MM, Goel V, Oxland T, Takata K, Duranceau J, Krag M, et al. Human lumbar vertebrae. Quantitative three-dimensional anatomy. Spine (Phila Pa 1976). 1992;17(3):299-306.
- Nikkhoo M, Haghpanahi M, Parnianpour M, Wang J-L. Dynamic responses of intervertebral disc during static creep and dynamic cyclic loading: A parametric poroelastic finite element analysis. Biomedical Engineering: Applications, Basis and Communications. 2013;25(01):1350013.
- Ferguson SJ, Ito K, Nolte LP. Fluid flow and convective transport of solutes within the intervertebral disc. J Biomech. 2004;37(2):213-21.
- Argoubi M, Shirazi-Adl A. Poroelastic creep response analysis of a lumbar motion segment in compression. J Biomech. 1996;29(10):1331-9.
- Lee KK, Teo EC. Poroelastic analysis of lumbar spinal stability in combined compression and anterior shear. J Spinal Disord Tech. 2004;17(5):429-38.
- Castro AP, Wilson W, Huyghe JM, Ito K, Alves JL. Intervertebral disc creep behavior assessment through an open source finite element solver. J Biomech. 2014;47(1):297-301.
- Schmidt H, Kettler A, Heuer F, Simon U, Claes L, Wilke HJ. Intradiscal pressure, shear strain, and fiber strain in the intervertebral disc under combined loading. Spine (Phila Pa 1976). 2007;32(7):748-55.
- Shih SL, Liu CL, Huang LY, Huang CH, Chen CS. Effects of cord pretension and stiffness of the Dynesys system spacer on the biomechanics of spinal decompression- a finite element study. BMC Musculoskelet Disord. 2013;14:191.
- Khoz Z, Nikkhoo M, Cheng Ch. Parametric patient-specific finite element modeling of lumbar spine based on anatomical parameters. Iranian Journal of Orthopaedic Surgery. 2018;16:195-202.
- Naserkhaki S, Jaremko JL, Adeeb S, El-Rich M. On the load-sharing along the ligamentous lumbosacral spine in flexed and extended postures: Finite element study. J Biomech. 2016;49(6):974-82.
- Catia. Dassault Systèmes; 2017 [Available from: https://www.3ds.com/products-services/catia/.
- Altair HyperMesh. Michigan Altair Engineering, Inc; 2017 [Available from: https://altairhyperworks.com/product/hypermesh.
- Schmidt H, Bashkuev M, Dreischarf M, Rohlmann A, Duda G, Wilke HJ, et al. Computational biomechanics of a lumbar motion segment in pure and combined shear loads. J Biomech. 2013;46(14):2513-21.
- Shirazi-Adl A, Ahmed AM, Shrivastava SC. A finite element study of a lumbar motion segment subjected to pure sagittal plane moments. J Biomech. 1986;19(4):331-50.
- Whyne CM, Hu SS, Lotz JC. Burst fracture in the metastatically involved spine: development, validation, and parametric analysis of a three-dimensional poroelastic finite-element model. Spine (Phila Pa 1976). 2003;28(7):652-60.
- Schmidt H, Galbusera F, Wilke HJ, Shirazi-Adl A. Remedy for fictive negative pressures in biphasic finite element models of the intervertebral disc during unloading. Comput Methods Biomech Biomed Engin. 2011;14(3):293-303.
- Tyrrell AR, Reilly T, Troup JD. Circadian variation in stature and the effects of spinal loading. Spine (Phila Pa 1976). 1985;10(2):161-4.
- Schmidt H, Shirazi-Adl A, Galbusera F, Wilke HJ. Response analysis of the lumbar spine during regular daily activities--a finite element analysis. J Biomech. 2010;43(10):1849-56.
- Heuer F, Schmidt H, Klezl Z, Claes L, Wilke HJ. Stepwise reduction of functional spinal structures increase range of motion and change lordosis angle. J Biomech. 2007;40(2):271-80.
- Reilly T, Tyrrell A, Troup JD. Circadian variation in human stature. Chronobiol Int. 1984;1(2):121-6.
- McMillan DW, Garbutt G, Adams MA. Effect of sustained loading on the water content of intervertebral discs: implications for disc metabolism. Ann Rheum Dis. 1996;55(12):880-7.
- Adams MA, McMillan DW, Green TP, Dolan P. Sustained loading generates stress concentrations in lumbar intervertebral discs. Spine (Phila Pa 1976). 1996;21(4):434-8.
)