Three-Dimensional Finite Element Analysis of L4-5 Degenerative Lumbar Disc Traction under Different Pushing Heights

Ding, Huaili; Liao, Lijun; Yan, Peichun; Zhao, Xiaolin; Li, Min

doi:https://doi.org/10.1155/2021/1322397

Journal of Healthcare Engineering

On this page

Abstract Introduction Materials and Methods Results Conclusions Discussion Data Availability Consent Conflicts of Interest Authors’ Contributions Acknowledgments References Copyright Related Articles

Special Issue

Data and Intelligence in Cyber Health

View this Special Issue

Research Article | Open Access

Volume 2021 | Article ID 1322397 | https://doi.org/10.1155/2021/1322397

Three-Dimensional Finite Element Analysis of L4-5 Degenerative Lumbar Disc Traction under Different Pushing Heights

Huaili Ding,¹Lijun Liao,¹Peichun Yan,¹Xiaolin Zhao,¹and Min Li¹

Academic Editor: Khin Wee Lai

Received14 May 2021

Accepted02 Jul 2021

Published22 Jul 2021

Abstract

Objective. To study and analyze the changes of intervertebral foramen height and area of the degenerative L4-5 intervertebral disc under different pushing heights by the finite element method. Methods. CT and MRI images of T12-S1 segments were obtained from a healthy volunteer who met the inclusion criteria. A DR machine was used to capture images of the lumbar lateral section before and after simultaneous pushing of the L4 and L5 spinous processes by manipulation called Daogaijinbei, and the measurement showed that the displacement changes of L4 and L5 were both approximately 10 cm, so the pushing height was set at 0–10 cm. A three-dimensional finite element model of the entire normal lumbar spine was established using Mimics 16.0, Geomagic Studio 2014, Hypermesh 13.0, MSC.Patran 2012, and so on. The disc height and nucleus area of the lumbar disc of the normal entire lumbar disc model were adjusted to establish models of the L4-5 disc with mild, moderate, and severe degeneration. Changes of disc height and area of the L4-5 degenerative intervertebral disc under different pushing heights were calculated. Results. The size of the L4-5 intervertebral foramen was analyzed from the height and area of the intervertebral foramen, and the results showed the following: (1) as for the normal lumbar disc and a lumbar of the L4-5 disc with mild and moderate degeneration, the height of the L4-5 intervertebral foramen and its area both increased during pushing between 0 and 8 cm. After the pushing height reached 8 cm, the height and area of the L4-5 intervertebral foramen gradually became stable; (2) as for the L4-5 disc with severe degeneration, during the process of pushing, the height and area of the L4-5 intervertebral foramen increased slightly, but this change was not obvious. Conclusions. After the spinal manipulation, the sizes of the L4-5 intervertebral foramen of the L4-5 disc with mild and moderate degeneration were significantly larger than those before pushing; in contrast, the size of L4-5 intervertebral foramen of the L4-5 disc with severe lumbar degeneration was not significantly changed.

1. Introduction

Lumbar intervertebral disc degeneration can cause biomechanical changes of the lumbar vertebrae and destroy their mechanical environment, thus resulting in a series of degenerative diseases of the lumbar intervertebral discs. Lumbar disc degeneration is a complex process in which the lumbar spine is affected by a variety of complex factors and changes in stroma composition or stroma state, and it is confluence by lumbar disc aging, biomechanical factors [1–3], nutritional factors [4–6], and so on. The main clinical manifestations of degenerative disease of the lumbar intervertebral discs are lumbago and leg pain, which are mainly caused by extrusion of a nerve root by protrusion after degeneration of the lumbar intervertebral disc. As lots of studies [7, 8] have showed, traditional Chinese massage therapy can effectively treat lumbar intervertebral disc degeneration. In this study, the relationship between the degree of L4-5 intervertebral disc degeneration and the corresponding intervertebral height and area was studied by the finite element analysis under different pushing heights of the traditional Chinese manipulation called Daogaijinbei.

2. Materials and Methods

2.1. Research Object

Creation of a three-dimensional finite element model of the whole normal lumbar and the models of the L4-5 disc with mild, moderate, and severe degeneration established based on CT and MRI images of the lumbar vertebrae of a healthy male volunteer, aged 26 years, with height 171 cm, weight 60.5 kg, no lumbar spinal stenosis, lumbar spondylolisthesis, or lumbar instability, and no ossification of intervertebral space, posterior longitudinal ligament, ligamentum flavum, vertebral anterior margin hyperplasia, and complex disease is not suitable for the study. The finite element method can simulate the material properties, morphological structure, boundary conditions, and load conditions of lumbar vertebrae in the digital form, observe any parameter that can change the influence on the mechanics of the entire lumbar vertebrae, and can be repeated loading and arbitrary testing. In addition, the research object referred to most literatures was 1 healthy adult, so the research object in this study was 1 healthy adult male volunteer.

2.2. Establishment of a Normal Full Lumbar Finite Element Model

In the imaging department of the affiliated rehabilitation hospital of Fujian University of Traditional Chinese Medicine, a healthy male volunteer was examined by CT imaging. Imaging examination and clinical confirmation showed that the volunteer was in good health, and his lumbar vertebrae were normal. The man voluntarily agreed to participate in the research subject and signed the informed consent form. The study passed the requirements of the hospital ethics committee. Images of the lumbar spine were acquired by CT scanning of T12-S1 segments using a Toshiba Aquilion 16 slice spiral CT, and MRI images were acquired using a Siemens 3.0 T MRI machine to continuously scan the lumbar spine. The scanning data were directly saved in a DICOM (digital imaging and communications in medicine) format. After completion of the lumbar CT and MRI image acquisition, Mimics 16.0, Geomagic Studio 2014, and Hypermesh 13.0 were used to complete the creation of a three-dimensional finite element model of the lumbar spine. After the model was established, MSC.Patran 2012 was used to calculate the size of the L4-5 intervertebral foramen under different pushing heights. Due to the complex range of motion of intervertebral joints, facet joints were not taken into account in this study, and the study was limited to the changes of intervertebral foramen.

The preprocessed CT images were imported into Mimics software, and the mask between the vertebrae was divided by edit masks. Then, the region growing tool was used to establish the vertebral lamina, and the vertebral model was generated from these laminas. Next, the model was filled with cavities and deburred using edit masks to obtain a model of the lumbar vertebrae with a smooth surface. Finally, a 3D model of the entire lumbar spine was generated using Calculate 3D, and the 3D model was postprocessed and meshed. Through the above steps, a finite element mesh model of the whole lumbar region was obtained as shown in Figure 1, and the mesh model of the intervertebral disc is shown in Figures 2 and 3.

Referring to the published literature [9–13], the parameter settings of L1-L5 lumbar vertebrae structural materials are given in Table 1.

3. Establishment of L4-5 Degenerative Disc Models

At present, there is no uniform standard for the specific grade of intervertebral disc degeneration. In this study, referring to the MRI grading standard of the Pfirrmann classification of intervertebral disc degeneration [14] (Table 2) and the Thompson five-level system [15], the process of intervertebral disc degeneration was considered from two aspects: disc height and nucleus pulposus area.

Classes I and II were categorized as a normal lumbar disc, III as mild degeneration, and IV and V as severe degeneration. On the basis of the effective normal three-dimensional finite element model, the degeneration of a lumbar intervertebral disc was simulated by reducing the height of the intervertebral disc and the area of the nucleus pulposus. The L4-5 intervertebral disc mild degeneration model was generated from the normal lumbar disc model by reducing the height by 15% while keeping the nucleus pulposus area unchanged, the moderate L4-5 intervertebral disc degeneration model was obtained by a 33% height reduction and a 67% nucleus pulposus area reduction, and the severe L4-5 intervertebral disc degeneration model was obtained by a 70% height reduction and a 67% nucleus pulposus area reduction.

Based on the above L4-5 intervertebral disc degeneration classification method, the parameters of the normal lumbar model were adjusted to establish whole lumbar models of mild, moderate, and severe degeneration of the L4-5 spine and of the L4-5 intervertebral disc, as shown in Figures 4–12.

3.1. Loading Mechanical Parameters

Lumbar vertebrae account for approximately 13.9% of human bodyweight [16]. The weight of the volunteer used for this model was 60.5 kg, making the weight of lumbar vertebrae approximately 8.41 kg, which is 84.1 N. Meanwhile, considering the pressure between other tissues, this study measured the waist pressure range of the doctor’s palm on the subject’s waist during the normal performance of lumbar extension with the help of a pressure sensor, and the maximum pressure was 99 N [17]. Thus, for convenience of calculation, the force applied to the L4 and L5 spinous processes was defined as 100 N in the finite element analysis.

The corresponding nodes on the upper surface of vertebral body L1 and the lower surface of vertebral body L5 are restrained by sliding support, and the displacement changes in flexion and extension directions are allowed, while different degrees of displacement are released. Considering the changes in the normal physiological curvature of the lumbar spine, the stress changes of the lumbar spine under different jacking heights were studied.The force is applied perpendicular to the spinous process nodes, and some nodes of L4 and L5 are constrained by releasing different displacement boundaries on the surfaces of L4 and L5 anterior longitudinal ligaments. The lateral sections of the lumbar spine before and after the action of the spinous process of L4 and L5 with the reverse cover were taken with a DR camera, respectively. The safe movement range of the spinous process of L4 and L5 was measured to be 10 cm. After the validity verification of the lumbar model, the mechanical parameters were loaded. The safe range of L4 and L5 lumbar vertebrae was within 10 cm, and the dynamic loading was carried out with 1 cm as the step length.The stress and distribution of each lumbar structure were observed. During the simulation calculation, the specific boundary constraint conditions of the spine manipulation bed simulating the inverted cover method on the whole lumbar spine are shown in Figure 13.

3.2. Calculation Method of Intervertebral Foramen Size

Research [18] shows that the maximum height and area of an intervertebral foramen decrease with the aggravation of intervertebral disc degeneration, while the maximum width of an intervertebral foramen is not affected by intervertebral disc degeneration, and that intervertebral disc degeneration is closely related to the height and area of the intervertebral foramen, and the calculation method used in the model established in this study is shown in Figure 14.

Using the method shown in the figure, we analyzed the changes in the lumbar structure when subjected to the stress of a normal L4-5 intervertebral disc and degenerative L4-5 discs with different degrees of degeneration, under different pushing heights of 0–10 cm. The size of the intervertebral foramen was measured and analyzed in the normal group, mild degeneration group, moderate degeneration group, and severe degeneration group.

4. Results

In the normal, mild degeneration, moderate degeneration, and severe degeneration groups, the L4 and L5 spinous processes were pushed up to a height of 0–10 cm, and the finite element experiment was conducted every 1 cm step. Considering that the L4-L5 structure is relatively symmetrical, the height and area of the intervertebral foramen on the one side were measured. As the results showed, the L4-5 intervertebral foramen heights of a normal lumbar disc and a lumbar of the L4-5 disc with mild and moderate degeneration were 20.1–25.9 mm, 18.4–24.9 mm, 14.5–20.2 mm, respectively, and the areas of the L4-5 intervertebral foramen were 235.2–341.9 mm², 199.9–313.3 mm², and 166.8–253.6 mm², respectively. As for the L4-5 disc with severe degeneration, during the pushing process, the height of the L4-5 intervertebral foramen was 9.1–10.2 mm, and the area was 121.6–154.1 mm². The changes during the pushing with different pushing heights are shown in Figure 15.

(a)

(b)

From Figure 15, the heights of the L4-5 intervertebral foramen and its areas of the normal lumbar disc and a lumbar of the L4-5 disc with mild and moderate degeneration both increased when the pushing height is between 0 and 8 cm, and the height and area of the L4-5 intervertebral foramen gradually became stable when the pushing height reached 8 cm or larger; however, the height and area of the L4-5 intervertebral foramen of the L4-5 disc with severe degeneration increased slightly, but this change was not obvious during the pushing.

5. Conclusions

After the spinal manipulation, the sizes of the L4-5 intervertebral foramen of the L4-5 disc with mild and moderate degeneration were significantly larger than those before pushing; in contrast, the size of L4-5 intervertebral foramen of the L4-5 disc with severe lumbar degeneration was not significantly changed.

6. Discussion

The size of the intervertebral foramen is closely related to the normal physiological morphology of spinal nerves and blood vessels, and some studies have shown that intervertebral disc degeneration is closely related to the height and area of the intervertebral foramen [19–22]. In this study, we analyzed the stress changes of the degenerative L4-5 intervertebral disc during pushing to a height of 0–10 cm and also analyzed the height and area of intervertebral foramen in these biomechanical environments. The height and area of the L4-5 foramen in mild, moderate, and severe degeneration of the L4-5 intervertebral disc were smaller than those of the normal lumbar spine when the pushing height was 0 cm, i.e., with no pushing. The height and area of the L4-5 intervertebral foramen decreased with increased L4-5 intervertebral disc degeneration. For the lumbar vertebrae with mild and moderate degeneration of the L4-5 intervertebral disc, the height and area of the L4-5 intervertebral foramen showed a continuously increasing trend when the pushing height was 0–8 cm. When the pushing height was 8–10 cm, the height and area of the L4-5 intervertebral foramen gradually became stable, which indicated that the height and area of the L4-5 intervertebral foramen of mild and moderate degeneration of the L4-5 intervertebral disc could be improved under spinal manipulation, and the height and area of the L4-5 intervertebral foramen tended to remain stable after the pushing height reached 8 cm. For the L4-5 intervertebral disc with severe degeneration, the height and area of L4-5 intervertebral foramen did not increase significantly; therefore, the pushing effect of spinal manipulation did not have a significant effect on the size of the L4-5 intervertebral foramen [23, 24].

Data Availability

The data used to support the findings of this study cannot be shared as no datasets were generated or analyzed during the current study.

Informed consent was obtained from all individual participants included in the study references.

Conflicts of Interest

The authors declare that there are no conflicts of interest.

Authors’ Contributions

Huaili Ding and Lijun Liao contributed equally to the article.

Acknowledgments

The authors acknowledge the Research and Development of Spinal Draft Bed, Science and Technology Program of Fujian Province (2018Y4006), 2018.4-2022.4, Biomechanical Mechanism Analysis of the Treatment of Lumbar Intervertebral Disc Herniation by “DaoGaiJinBei” Manipulation in South Shaolin, Project of Fujian Maternity and Child Health Hospital (YCXB20-03), 2020.12.01-2023.11.30, and the Mechanism of “Liu’s Acupuncture” in the Treatment of Autism Model Rats Based on mTOR Signaling Pathway, Project of Fujian Maternity and Child Health Hospital (YCXM20-17), 2020.12.01-2023.11.30.

References

F. M. K. Williams and P. N. Sambrook, “Neck and back pain and intervertebral disc degeneration: R,” Best Practice & Research Clinical Rheumatology, vol. 25, no. 1, pp. 69–79, 2011.
View at: Publisher Site | Google Scholar
C. A. Ranson, R. W. Kerslake, A. F. Burnett, M. E. Batt, and S. Abdi, “Magnetic resonance imaging of the lumbar spine in asymptomatic professional fast bowlers in cricket,” Journal of Bone and Joint Surgery British Volume, vol. 87-B, no. 8, pp. 1111–1116, 2005.
View at: Publisher Site | Google Scholar
G. Cinotti, C. D. Rocca, S. Romeo, F. Vittur, R. Toffanin, and G. Trasimeni, “Degenerative changes of porcine intervertebral disc induced by vertebral endplate injuries,” Spine, vol. 30, no. 2, pp. 174–180, 2005.
View at: Publisher Site | Google Scholar
J. Niinimäki, A. Korkiakoski, O. Parviainen et al., “Association of lumbar artery narrowing, degenerative changes in disc and endplate and apparent diffusion in disc on postcontrast enhancement of lumbar intervertebral disc,” Magnetic Resonance Materials in Physics, Biology and Medicine, vol. 22, no. 2, pp. 101–109, 2009.
View at: Publisher Site | Google Scholar
A. Korkiakoski, J. Niinimäki, J. Karppinen et al., “Association of lumbar arterial stenosis with low back symptoms: a cross-sectional study using two-dimensional time-of-flight magnetic resonance angiography,” Acta Radiologica, vol. 50, no. 1, pp. 48–54, 2009.
View at: Publisher Site | Google Scholar
H. E. Gruber, N. Ashraf, J. Kilburn et al., “Vertebral endplate architecture and vascularization: application of micro-computerized tomography, a vascular tracer, and immunocytochemistry in analyses of disc degeneration in the aging sand rat,” Spine, vol. 30, no. 23, pp. 2593–2600, 2005.
View at: Publisher Site | Google Scholar
H. Xue, “Clinical study on the treatment of lumbar disc herniation with double positioning oblique pulling method,” Acta Sinica Sinica, vol. 32, no. 3, pp. 477–479, 2017.
View at: Google Scholar
J.-F. Chen, X.-N. Shu, S.-J. Tang, Y. Wu, and Y.-P. Zhang, “Influence of lumbar disc degeneration on the efficacy of lumbar fixed-point rotation manipulation in sitting position: a finite element study,” Journal of Acupuncture and Tuina Science, vol. 14, no. 4, pp. 295–299, 2016.
View at: Publisher Site | Google Scholar
S. Jin, Establishment and Biomechanical Analysis of Lumbar Finite Element Model, Dalian Medical University, Dalian, China, 2016.
H. Deng, Three-dimensional Finite Element Modeling of Lumbosacral Vertebrae and Biomechanical Analysis of Three Internal Fixation Methods, China Medical University, Taichung, Taiwan, 2011.
Y. Shen, Finite Element Analysis of Lumbar Spondylolisthesis with a New Posterior Dynamic Internal Fixation System, China Medical University, Taichung, Taiwan, 2012.
C.-C. Lo, K.-J. Tsai, Z.-C. Zhong, S.-H. Chen, and C. Hung, “Biomechanical differences of Coflex-F and pedicle screw fixation combined with TLIF or ALIF—a finite element study,” Computer Methods in Biomechanics and Biomedical Engineering, vol. 14, no. 11, pp. 947–956, 2011.
View at: Publisher Site | Google Scholar
Y. Wang, X.-D. Yi, and C.-D. Li, “The influence of artificial nucleus pulposus replacement on stress distribution in the cartilaginous endplate in a 3-dimensional finite element model of the lumbar intervertebral disc,” Medicine, vol. 96, no. 50, Article ID e9149, 2017.
View at: Publisher Site | Google Scholar
C. W. A. Pfirrmann, A. Metzdorf, M. Zanetti, J. Hodler, and N. Boos, “Magnetic resonance classification of lumbar intervertebral disc degeneration,” Spine, vol. 26, no. 17, pp. 1873–1878, 2001.
View at: Publisher Site | Google Scholar
J. P. Thompson, R. H. Pearce, M. T. Schechter, M. E. Adams, I. K. Y. Tsang, and P. B. Bishop, “Preliminary evaluation of a scheme for grading the gross morphology of the human intervertebral disc,” Spine, vol. 15, no. 5, pp. 411–415, 1990.
View at: Publisher Site | Google Scholar
H.-J. Kim, H.-J. Chun, K.-T. Kang et al., “The biomechanical effect of pedicle Screwsʼ insertion angle and position on the superior adjacent segment in 1 segment lumbar fusion,” Spine, vol. 37, no. 19, pp. 1637–1644, 2012.
View at: Publisher Site | Google Scholar
Y. J. Jin, Y. E. Kim, J. H. Seo, H. W. Choi, and T.-A. Jahng, “Effects of rod stiffness and fusion mass on the adjacent segments after floating mono-segmental fusion: a study using finite element analysis,” European Spine Journal, vol. 22, no. 5, pp. 1066–1077, 2013.
View at: Publisher Site | Google Scholar
M. Sheng, L. fan, N. Wang et al., “Imaging changes and clinical significance of intervertebral foramen in different degrees of intervertebral disc degeneration,” Chinese Journal of Orthopedics, vol. 22, no. 10, pp. 730–732, 2009.
View at: Google Scholar
S. Naserkhaki, J. L. Jaremko, M. El-Rich et al., “Effects of inter-individual lumbar spine geometry variation on load-sharing: geometrically personalized finite element study,” Journal of Biomechanics, vol. 49, 2016.
View at: Publisher Site | Google Scholar
J. Huang, H. Yan, F. Jian, X. Wang, and H. Li, “Numerical analysis of the influence of nucleus pulposus removal on the biomechanical behavior of a lumbar motion segment,” Computer Methods in Biomechanics and Biomedical Engineering, vol. 18, no. 14, pp. 1516–1524, 2015.
View at: Publisher Site | Google Scholar
C. Song, X. F. Li, Z. D. Liu et al., “Biomechanical assessment of a novel L4/5 level interspinous implant using three-dimensional finite element analysis,” European Review for Medical and Pharmacological Sciences, vol. 18, no. 1, pp. 86–94, 2014.
View at: Google Scholar
D. V. Ambati, E. K. Wright, R. A. Lehman, D. G. Kang, S. C. Wagner, and A. E. Dmitriev, “Bilateral pedicle screw fixation provides superior biomechanical stability in transforaminal lumbar interbody fusion: a finite element study,” The Spine Journal, vol. 15, no. 8, pp. 1812–1822, 2015.
View at: Publisher Site | Google Scholar
G. Li-Xin and F. Wei, “Impact of material properties of intervertebral disc on dynamic response of the human lumbar spine to vertical vibration: a finite element sensitivity study,” Medical & Biological Engineering & Computing, vol. 57, 2018.
View at: Publisher Site | Google Scholar
L. Li, T. Shen, and Y.-K. Li, “A finite element analysis of stress distribution and disk displacement in response to lumbar rotation manipulation in the sitting and side-lying positions,” Journal of Manipulative and Physiological Therapeutics, vol. 40, no. 8, pp. 580–586, 2017.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2021 Huaili Ding et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

709

Downloads

679

Citations