Numerical study of lateral bending influence on lumbar intervertebral disc
Sandra Mikuckytė^{1} , Vytautas Ostaševičius^{2}
^{1, 2}Institute of Mechatronics, Kaunas University of Technology, Kaunas, Lithuania
^{1}Corresponding author
Vibroengineering PROCEDIA, Vol. 15, 2017, p. 7176.
https://doi.org/10.21595/vp.2017.19401
Received 4 November 2017; accepted 14 November 2017; published 1 December 2017
The objective of this study is to investigate the influence of in the frontal plane constrained and free lateral bending moves on lumbar spine intervertebral disc and to decide whether lateral bending could be beneficial for lumbar spine rehabilitation. A three dimensional poroelastic finite element model of L4L5 intervertebral disc was used to investigate the influence of lateral bending moves on healthy and degenerated lumbar intervertebral discs. Results of finite element modeling show that stress values in moderately degenerated annulus fibrosus increase more when lateral bending is combined with flexion and extension, also, additional flexion, extension and axial rotation movements reduce fluid velocity in nucleus pulposus compared with pure (in the frontal plane constrained) lateral bending. These results suggest that pure lateral bending exercise could be more beneficial than other moves in order to achieve better nutrition of intervertebral disc, and that degeneration of intervertebral discs should be evaluated before defining an intensity of lateral bending exercise.
Keywords: Intervertebral disc, lateral bending, finite element analysis, poroelasticity.
1. Introduction
While flexion and extension of lumbar spine are performed every day multiple times, lateral bending is significantly less frequent during daily activities and therefore less researched, although lateral bending has notable influence on spine because it is primarily accomplished by movement of the spine [1].
Pure lateral bending moment is often applied to the spinal finite element models as one of the loads (together with flexion, extensions or axial rotation) in various intervertebral disc degeneration or finite element sensitivity analysis [24] but the influence of lateral bending on lumbar spine is not the main concern of these studies. Their results related with lateral bending show that maximum stress values in annulus fibrosus caused by lateral bending are higher compared with stress values caused by flexion, extension or axial rotation and increase even more when intervertebral disc is degenerated, on the other hand, lateral bending causes the lowest facet joint forces. Also, increased intervertebral disc degeneration reduces range of motion of lateral bending and intradiscal pressure even obtains negative values, however, failure of ligaments does not have significant influence on range of motion during lateral bending. Modeling assumptions, chosen in mentioned studies, do not allow to evaluate fluid flow paths and pressure distribution within intervertebral disc and to predict the effects of various loads on intervertebral disc nutrition.
Fluid flow within the intervertebral disc is usually evaluated by using poroelastic models of intervertebral disc, but current studies are mostly concentrated on investigation of axial compression load and diurnal cycle influence on intervertebral disc [56] and effects of other loads are not sufficiently researched.
The purpose of this study is to investigate the influence of lateral bending on stress values, pore pressure and fluid flow velocity of healthy and degenerated lumbar spine intervertebral discs in order to explore the possibilities to use lateral bending exercises for lumbar spine rehabilitation.
2. Methods
A three dimensional poroelastic finite element model of 35.9 mm wide and 51.9 mm long lumbar spine intervertebral disc L4L5 with a total height of 11.4 mm was created with Comsol MultiPhysics (COMSOL Inc.) software. Model consisted of nucleus pulposus (NP), which occupied 43 % of disc volume, annulus fibrosus (AF) and 0.5 mm thick cartilaginous endplates (CEP). Composite behaviour of annulus fibrosus was represented by annulus fibrosus matrix and seven layers of annulus fibers which were modeled as tensiononly elements with an inclination angle of ±24° (Fig. 1).
Fig. 1. L4L5 intervertebral disc: a) disc geometry with distinct regions of nucleus pulposus, annulus fibrosus and cartilage endplates, b) finite element model
a)
b)
Intervertebral disc was considered as fluid saturated porous media and calculation of fluid flow was based on Darcy’s generalized law. Tissues permeability was considered isotropic and straindependant [7]:
where ${k}_{0}$ – initial permeability; ${e}_{0}$ – initial voids ratio. Voids ratio $e$ is defined as:
where $\varphi $ – tissue porosity, which depends on volumetric strain ${\epsilon}_{vol}$:
Annulus fibers were defined with elastic modulus of 100 MPa and Poisson ratio of 0.33. Other material properties representing healthy (1st grade) and moderately degenerated (approximately 3rd grade) discs were taken from literature and are summarized in Table 1.
Intervertebral disc swelling was simulated by applying osmotic pressure gradient:
where ${\varphi}_{i}$ and ${\varphi}_{e}$ – internal and external osmotic coefficients, $R$ – universal gas constant, $T$ – absolute temperature, ${c}_{F}$ – fixed charge density, ${\gamma}_{i}$ and ${\gamma}_{e}$ – activity coefficients, ${c}_{e}$ – external salt concentration.
Fixed charge density depends on volumetric deformations and is defined as:
where ${c}_{F0}$ – initial fixed charge density, ${n}_{0}$ – initial fluid volume fraction, $J$ – determinant of the deformation gradient tensor.
An inferior surface of intervertebral disc was constrained in displacement. Model was validated by applying pure lateral bending moments of 110 Nm and checking simulation results with in vivo results given in literature [11].
In order to investigate the influence of in frontal plane constrained lateral bending move, 10° lateral bending angle was prescribed to the disk. This angle was chosen as close to the maximum failurefree lateral bending angle. During time transient analysis maximum bending value was achieved in one second time. In addition to lateral bending in frontal plane, small rotations of 2° and 5° in sagittal and transverse planes were applied to simulate situations when lateral bending move is performed unconstrained (not strictly in frontal plane) due to the physical condition of the subject.
Table 1. Material properties of healthy and degenerated discs used in finite element analysis (Material properties are set according to Maladrino et al., 2009; Natarajan et al., 2006; Ferguson et al., 2004 [810]; or assumed)
Spinal component

Elastic modulus (MPa)

Poisson ratio

Initial porosity

Initial permeability (mm^{4}/Ns)

M


Healthy
disc

Annulus fibrosus

2.5

0.17

0.75

7.5·10^{16}

8.5

Nucleus Pulposus

1

0.17

0.83

7.5·10^{16}

8.5


Cartilage endplates

5

0.17

0.80

7.5·10^{15}

8.5


Degenerated
disc

Annulus fibrosus

12

0.15

0.60

8.5·10^{16}

8.5

Nucleus Pulposus

1.6

0.15

0.78

8.5·10^{16}

8.5


Cartilage endplates

4

0.15

0.80

5·10^{15}

8.5

3. Results
Fig. 2 shows pure 10° LB simulation results in frontal section view. Stress values in annulus fibrosus is significantly higher than in nucleus pulposus. Within healthy intervertebral disc the highest von Misses stress zones occur at the lateral posterior annulus rim and the maximum value is 3.15 MPa. Stress values significantly increase in annulus fibrosus of degenerated intervertebral disc and reach value of 9.12 MPa. However, there are no significant changes of stress values within nucleus pulposus of healthy and degenerated discs. Disc degeneration also increases absolute values of pore pressure and pressure distribution between nucleus pulposus and annulus fibrosus becomes less smooth. Calculations show that fluid flow velocity in degenerated disc also increases, but higher fluid velocity value in degenerated disc does not mean that disc degeneration has any beneficial effects. Due to the reason that exact lateral bending angle was prescribed in this study instead of lateral bending moments, it could be assumed that reactions and strain values of degenerated disc are higher, as degeneration tends to increase a stiffness of intervertebral disc, so that the same bending angle could induce critical stress values and increase risk of failure.
The results of combined loading modelling (Fig. 3) reveal that additional 2° flexion, extension or axial rotation moves lead to only minor changes in stress values of healthy annulus, while additional 5° moves change maximal stress values more noticeably. It was calculated that 2° additional flexion and extension moves cause a minor increase in annulus fibrosus stress values, while 2° axial rotations reduce the stress level of pure lateral bending. Additional loads of 5° reduce maximal stress values in all four combined loading cases. Both 2° and 5° additional flexion and extension moves increase stress values in moderately degenerated annulus, while 2° and 5° axial rotations reduce stress values compared with pure lateral bending.
Fig. 2. Results of pure 10° lateral bending in a) healthy; b) degenerated intervertebral discs. Frontal plane section view
a) Health disc
b) Degenerated disc
Fig. 3. Von Misses stress values induced by pure lateral bending and additional loads of 2° and 5° in healthy (1st grade) and degenerated (3rd grade) discs: PBL – pure 10° lateral bending, $E$ – lateral bending + extension; $F$ – lateral bending + flexion; LAR and RAR – lateral bending + axial rotation to the left or right
Combined loading and disc degeneration also have the influence on intervertebral disc pore pressure (Fig. 4). Nucleus pulposus pore pressure is almost 2.5 times higher and annulus fibrosus pore pressure is almost 4.5 times higher in degenerated disc. While changes in the nucleus pulposus pore pressure due to additional rotations are insignificant in both – healthy and degenerated discs, changes in annulus fibrosus are more noticeable and their patterns are different for healthy and degenerated discs. Left axial rotation has the most influence on healthy disc annulus fibrosus pore pressure, while additional extension affects pore pressure of degenerated annulus fibrosus the most.
Combined loading also influences fluid flow velocity through the intervertebral disc (Fig. 5). Fluid velocity within healthy nucleus pulposus is reduced by all 2° and 5° additional moves compared with pure lateral bending values. 2° additional extension move has the most significant influence of all 2° moves while 5° additional rotations reduce nucleus pulposus fluid velocity the most. The highest fluid flow velocity is at the boundaries of an intervertebral disc and at the circumferential area between nucleus pulposus and annulus fibrosus.
Fig. 4. Pore pressure comparison between loading cases: a) healthy intervertebral disc; b) degenerated intervertebral disc (PBL – pure 10° lateral bending, $E$ – lateral bending + 2° extension; $F$ – lateral bending + 2° flexion; LAR and RAR – lateral bending + 2° axial rotation to the left or right)
a)
b)
Fig. 5. Fluid velocity in healthy nucleus pulposus: PBL – pure 10° lateral bending, $E$ – lateral bending + extension; $F$ – lateral bending + flexion; LAR and RAR – lateral bending + axial rotation to the left or right
4. Conclusions
The results of lateral bending of poroelastic intervertebral disc show that stress values in intervertebral disc increase when more degenerated intervertebral disc is bent, and this agrees with results of the study where nucleus pulposus was modeled as fluid filled cavity and annulus fibrosus matrix was assumed hyperelastic [2], so various intervertebral disc modeling methods present similar results of lateral bending induced stress values of intervertebral disc.
Disc degeneration also increases pore pressure of intervertebral disc, especially within annulus fibrosus, and this could lead to intervertebral disc disorders such as disk herniation.
Due to significant disc degeneration influence, the condition of patient’s spine should be evaluated before defining an intensity of lateral bending exercise. While combined loading simulations do not show any significant changes in annulus stress values and more comprehensive study is necessary to define exact relationships between combination of different loads and their induced stress, obtained results could imply that additional flexion and extension loads are more dangerous than axial rotations and could induce higher annulus fibrosus stresses than pure lateral bending.
Fluid flow velocity reduced during combined loading could suggests that in frontal plane constrained lateral bending exercise is more beneficial than other moves in order to achieve better nutrition of intervertebral disc.
Acknowledgements
This research was funded by a Grant (No. SEN10/15) from the Research Council of Lithuania. Project acronym: “CaSpine”.
References
 Lee R. Y., Wong T. K. Relationship between the movements of the lumbar spine and hip. Human Movement Science, Vol. 21, Issue 4, 2002, p. 481494. [Search CrossRef]
 Rohlmann A., Zander Th, Schmidt H., Wilke H. J., Bergmann G. Analysis of the influence of disc degeneration on the mechanical behaviour of a lumbar motion segment using a finite element method. Journal of Biomechanics, Vol. 39, 2006, p. 24842490. [Search CrossRef]
 Ellingson A. M., Shaw M. N., Giambini H., An K. Comparative role of disc degeneration and ligament failure on functional mechanics of the lumbar spine. Computer Methods in Biomechanics and Biomedical Engineering, Vol. 19, Issue 9, 2015, p. 10091018. [Search CrossRef]
 Zander Th., Dreischarf M., Timm A., Baumann W. W., Schmidt H. Impact of material and morphological parameters on the mechanical response of the lumbar spine – a finite element sensitivity study. Journal of Biomechanics, Vol. 53, 2017, p. 185190. [Search CrossRef]
 Massey Ch. J., van Donkelaar C. C., Vresilovic E., Zavaliangos A., Marcolongo M. Effects of aging and degeneration on the human intervertebral disc during the diurnal cycle: a finite element study. Journal of Orthopaedic Research, Vol. 30, Issue 1, 2012, p. 122128. [Search CrossRef]
 Malandrino A., Noailly J., Lacroix D. Numerical exploration of the combined effect of nutrient supply, tissue condition and deformation in the intervertebral disc. Journal of Biomechanics, Vol. 47, 2014, p. 15201525. [Search CrossRef]
 Argoubi M., Shirazi Adl A. Poroelastic creep response analysis of a lumbar motion segment in compression. Journal of Biomechanics, Vol. 29, Issue 10, 1996, p. 13311339. [Search CrossRef]
 Maladrino A., Planell J. A., Lacroix D. Statistical factorial analysis on the poroelastic material properties sensitivity of the lumbar intervertebral disc under compression, flexion and axial rotation. Journal of Biomechanics, Vol. 42, 2009, p. 27802788. [Search CrossRef]
 Natarajan R. N., Williams J. R., Andersson G. B. J. Modeling changes in intervertebral disc mechanics with degeneration. The Journal of Bone and Joint Surgery, Vol. 88, Issue 2, 2006, p. 3640. [Search CrossRef]
 Ferguson S. J., Ito K., Nolte L. P. Fluid flow and convective transport of solutes within the intervertebral disc. Journal of Biomechanics, Vol. 37, 2004, p. 213221. [Search CrossRef]
 Schmidt H., Kettler A., Heuer F., Simon U., Claes L., Wilke H. J. Intradiscal pressure, shear strain, and fiber strain in the intervertebral disc under combined loading. Spine, Vol. 32, Issue 7, 2007, p. 748755. [Search CrossRef]
Cited By
Computer Methods in Biomechanics and Biomedical Engineering
M. A. Chetoui, D. Ambard, P. Royer, S. Le Floc’h

2019
