Evidence-Based Complementary and Alternative Medicine
Volume 2018 (2018), Article ID 2798396, 8 pages
https://doi.org/10.1155/2018/2798396

Research Article

Biomechanical Effects of Lateral Bending Position on Performing Cervical Spinal Manipulation for Cervical Disc Herniation: A Three-Dimensional Finite Element Analysis

Xuecheng Huang,1,2 Linqiang Ye,² Zixian Wu,1,2 Lichang Liang,³ Qianli Wang,² Weibo Yu,⁴ De Liang,⁵ and Xiaobing Jiang²

¹First School of Clinical Medicine, Guangzhou University of Chinese Medicine, Guangzhou 510405, China
²Department of Spinal Surgery, First Affiliated Hospital of Guangzhou University of Chinese Medicine, Guangzhou 510405, China
³Shenzhen Hospital of Guangzhou University of Chinese Medicine, Shenzhen 518034, China
⁴Department of Orthopaedic and Traumatology, Tongde Hospital of Zhejiang Province, Hangzhou 310012, China
⁵Laboratory Affiliated to National Key Discipline of Orthopaedic and Traumatology of Chinese Medicine, Guangzhou University of Chinese Medicine, Guangzhou 510405, China

Correspondence should be addressed to De Liang; 453223076@qq.com and Xiaobing Jiang; spinedrjxb@sina.com

Received 31 January 2018; Accepted 2 May 2018; Published 11 June 2018

Academic Editor: Morry Silberstein

Copyright © 2018 Xuecheng Huang 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.

Abstract.

Background. Most studies report that the common position of cervical spinal manipulation (CSM) for treating symptomatic cervical disc herniation (CDH) is lateral bending to the herniated side. However, the rationality of lateral bending position on performing CSM for CDH is still unclear. Objective. The purpose of this study is to investigate the biomechanical effects of lateral bending position on performing CSM for CDH. Methods. A finite element (FE) model of CDH (herniated on the left side) was generated in C5-6 segment based on the normal FE model. The FE model performed CSM in left lateral bending position, neutral position, and right lateral bending position, respectively. Cervical disc displacement, annulus fiber stress, and facet joint stress were observed during the simulation of CSM. Results. The cervical disc displacement on herniated side moved forward during CSM, and the maximum forward displacements were 0.23, 0.36, and 0.45 mm in left lateral bending position, neutral position, and right lateral bending position, respectively. As the same trend of cervical disc displacement, the annulus fiber stresses on herniated side from small to large were 7.40, 16.39, and 22.75 MPa in left lateral bending position, neutral position, and right lateral bending position, respectively. However, the maximum facet stresses at left superior cartilage of C6 in left lateral bending position, neutral position, and right lateral bending position were 6.88, 3.60, and 0.12 MPa, respectively. Conclusion. Compared with neutral position and right lateral bending position, though the forward displacement of cervical disc on herniated side was smaller in left lateral bending position, the annulus fiber stress on herniated side was declined by sharing load on the left facet joint. The results suggested that lateral bending to the herniated side on performing CSM tends to protect the cervical disc on herniated side. Future clinical studies are needed to verify that.

1. Introduction

Cervical disc herniation (CDH) is a common cause of cervical radiculopathy which occurs in approximately 85.4 of every 100000 persons [1, 2]. Symptoms of CDH which include pain and disability innervated by the nerve root can arise from the nerve root compression, inflammation, or both [2]. Patients with symptomatic CDH are initially treated conservatively, with surgery reserved for those cases that remain unresponsive to conservative treatment. Up to 90% of patients with symptomatic CDH will have significant improvement in symptoms with conservative treatment [2].

Cervical spinal manipulation (CSM) is one of the important conservative treatments for symptomatic CDH [3–5]. Studies reported that CSM reduces pain in patients with symptomatic CDH [6–10]. The manipulation commonly included two key procedures of rotation to the healthy side and lateral bending to the affected (herniated) side [6–10]. Wu et al. suggested that compression of the nerve root is relieved by a small displacement between it and the herniated disk during CSM with rotation to the healthy side [11, 12]. However, the rationality of lateral bending to the herniated side on performing CSM for CDH is still unclear [13].

The purpose of the present study is to compare the biochemical effects of different lateral bending positions on performing CSM through three-dimensional finite element analysis, so as to evaluate the rationality of lateral bending to the herniated side on performing CSM for CDH.

2. Materials and Methods

2.1. A Finite Element Model of Intact Cervical Spine (C3–7)

Following the methods of Zhong Jun Mo et al. 2014 [14], a normal three-dimensional finite element (FE) model of intact cervical spine was built using digitized image data of a C3-7 motion segment. The image data of C3-7 was obtained by a 64-detector CT scanner (Discovery CT750 HD, GE Healthcare, USA) at 0.625-mm interval from a healthy female volunteer (25 years old, 55 kg, and 165 cm) without any radiographic evidence of degenerative sign.

The slice images were imported into medical image processing software (Mimics 10.1, Materialise Inc., Belgium) to reconstruct the vertebrae geometry volume of C3–7. The geometry of other structures (the annulus fibrosis, nucleus pulpous, facet cartilage), which were difficult to separate from the CT images, was modeled using the solid modeling software, SolidWorks 2014 (SolidWorks Corp, Dassault Systèmes, Concord, M A).

Finite element modeling software (ABAQUS 2016, Simulia Inc., USA) was used to build and mesh the cervical spine components. The vertebrae were made up of a solid volume (cancellous bone) and a layer of shell (cortical bone and endplate) with a thickness of 0.4 mm [14, 15]. The intervertebral disc was constructed as a continuum structure partitioned into nucleus pulposus and annulus fibrosus. The nucleus pulposus was 43% of the total disc volume and located slightly posterior to the center of the disc [16]. The annulus fibrosus was modeled as a composite structure: the annulus ground substance reinforced by inclined fibers acting at approximately 30° from the transverse plane [17]. Seven intervertebral ligaments were incorporated, including anterior longitudinal ligament (ALL), posterior longitudinal ligament (PLL), capsular ligament (CL), flaval ligament (FL), interspinous ligament (ISL), supraspinous ligament (SSL), and transverse ligament (TL) with the suggested insertion site [18]. The cancellous bone was meshed into tetrahedron elements (C3D4), while the cortical element was meshed into triangle shell elements (S3) [14]. The annulus ground substance, the nucleus pulposus, and anterior plate were meshed into hexahedron elements (C3D8R) [14]. All the ligaments were modeled as tension-only axial connector, and the annulus fiber was meshed as tension-only truss elements (T3D2) [14]. In the FE model of intact cervical spine, the total number of nodes and elements was 216287 and 966930, respectively (Figure 1). Convergence within 1% was achieved in the intact model to ensure that the results were not relevant to the mesh density [14, 15].

Figure 1: The normal FE model of C3-7.

2.2. Validation of the Normal FE Model of C3-7 and C5-6

The assigned material properties (Table 1) selected from various sources in the literature were assumed to be linear, homogeneous, and isotropic [19–21]. Tied contact interfaces were used to ensure that the disc and ligament were attached to the vertebra, preventing any relative movement during the simulations. Surface-based, finite-sliding contact with a friction coefficient 0 was defined for facet joints [22]. The validation of the normal model was conducted according to the published FE model and human cadaveric cervical spines. For model validation of C3-7, the inferior endplate of C7 was fixed at six degrees of freedom in the same way as in vitro experiments [23, 24]. A follower preload of 50N was applied to the superior endplate center of C3 to simulate the head weight in the normal FE model of C3-7 [15]. In addition to the follower preload, a moment of 1.0Nm applied to the superior endplate center of C3 to simulate flexion, extension, lateral bending, and axial rotation motions in this model [24]. Furthermore, the FE model of C5-6 was also validated for using the same experimental and simulated loading protocols, except the follower preload of 73.6N and the moment of 1.8Nm [23].

Table 1: Material properties of the C3-7 finite element model.

Component
Young’s modulus (MPa)
Poisson’s
ratio
Cross-sectional area (mm²)

Cortical bone
12000
0.29
-

Cancellous bone
100
0.29
-

Endplate
1200
0.29
-

Annulus ground substance
3.4
0.45
-

Annulus fiber
450
0.40
-

Nucleus pulposus
1
0.49
-

Cartilage
10.4
0.40
-

ALL
10
0.30
12

PLL
10
0.30
45

CL
10
0.30
14

FL
1.5
0.30
5

ISL
1.5
0.30
13

SSL
1.5
0.30
13

TL
17
0.30
10

ALL: anterior Longitudinal ligament; PLL: posterior longitudinal ligament; CL: capsular ligament; FL: flaval ligament; ISL: interspinous ligament; SSL: supraspinous ligament; TL: transverse ligament.

2.3. Simulation of Cervical Disc Herniation in C5–6 Spinal Motion Segment

Similar to simulation methods reported by Hussain et al. [25], the FE model of the degenerative cervical spine was generated based on the developed FE model of the normal cervical spine. The cervical disc herniation (herniated on the left side) was simulated in C5-6 spinal motion segment. The changes in geometry and material properties (Table 2) used to simulate the degeneration were adapted from the clinical classification of degeneration of the cervical spine and the results of the previously published literature [18, 25, 26]. For the herniated disc of C5-6, the left posterolateral annulus was weakened as mid-then-outer annulus fibers tear [27], allowing herniation of nuclear material into the outer annular structure as a contained protrusion, and part of the disc pass into the vertebral canal space as an extrusion, as shown in Figure 2.

Table 2: Elastic tissue material properties of the AF and NP.

Description
Element type
Young’s modulus (MPa)
Poisson’s ratio

Normal disc

    AF
3D solid (4 node)
2.50
0.40

    NP
3D solid (4 node)
1.50
0.49

Moderately degenerated disc

    AF
3D solid (4 node)
2.50
0.40

    NP
3D solid (4 node)
2.00
0.49

Severely degenerated disc

    AF
3D solid (4 node)
5.00
0.20

    NP
3D solid (4 node)
4.00
0.25

AF Young modulus represents the modulus of lateral AF.

AF: annulus fibrosus; NP: nucleus pulposus.

Figure 2: The cervical disc herniation model in C5-6.

The Simulation and Loading of CSM in C5–6 with CDH

In this research, we used the material properties of moderately degenerated disc to investigate the biomechanical effects of lateral bending position on performing CSM for CDH. The mechanical loading steps in sequence to simulate CSM were as follows [12]:

(1)The inferior endplate of C6 was fixed.(2)The C5 vertebra, along with the entire model, was rotated 2°, 0, -2° around y-axis to simulate the left lateral bending, neutral, and right lateral bending positions, respectively.(3)In the left lateral bending, neutral, and right lateral bending positions, the cervical FE model was rotated 2° to the right side (the opposite side of CDH) around z-axis (vertical axis) to simulate rotation to the right side, respectively.(4)The cervical FE model continued to rotate a further 0.5° to the right side within 0.15 seconds so as to simulate the high-velocity, low-amplitude CSM.

2.4. Analysis

Cervical disc displacement, annulus fiber stress (von Mises stress) and facet joint stress (von Mises stress) were observed during the simulation of CSM.

3. Results

3.1. Validation of the Normal FE Model

The range of motion (ROM) of each functional spinal unit in the FE model of C3-7 is shown in Figure 3(a), which is consistent with in vitro experimental results [24]. Under 1.0 Nm moment and 50N follower load, the overall ROM of C3-7 was 17.53° in flexion, 10.01° in extension, 8.96° in axial rotation, and 11.05° in lateral bending, respectively. The results of C5-6 FE model validation shown in Figure 3(b) were also similar to the previous FE model (Ganbat et al.) [19] and in vitro experimental data (Moroney et al.) [23].

(a)

(b)

Figure 3: Validation of the normal FE model of C3-7 and C5-6.

3.2. Cervical Disc Displacement

When performing CSM, the left posterolateral (the herniated side) cervical disc moved left, forward, and up in three positions: left lateral bending, neutral, and right lateral bending (Figure 4). The maximum forward displacements of left posterolateral (the herniated side) cervical disc were 0.23, 0.36, and 0.45 mm in the left lateral bending, neutral, and right lateral bending positions, respectively.

Figure 4: Cervical disc displacement in the left lateral bending, neutral, and right lateral bending positions during CSM.

3.3. Annulus Fiber Stress

When performing CSM, the maximum annulus fiber stresses were 22.52, 17.61, and 23.60 MPa in left lateral bending position, neutral position, and right lateral bending position, respectively. The distribution of annulus fibers stress in left lateral bending position was concentrated on the right posterolateral (the heathy side) of cervical disc. However, the distributions of annulus fibers stress in neutral position and right lateral bending position were both concentrated on the left posterolateral (the herniated side) of cervical disc (Figure 5). And the annulus fiber stress of the same node on herniated side was 7.40 MPa in left lateral bending position, which was smaller than those in neutral position (16.39 MPa) and in right lateral bending position (22.75 MPa) while performing CSM (Figure 5, black arrows).

Figure 5: Annulus fiber stresses in the left lateral bending, neutral, and right lateral bending positions during CSM.

3.4. Facet Joint Stress

When performing CSM, the maximum facet stresses at left superior cartilage of C6 in left lateral bending position, neutral position, and right lateral bending position were 6.88, 3.60, and 0.12 MPa, respectively (Figure 6).

Figure 6: Facet joint stresses in the left lateral bending, neutral, and right lateral positions bending during CSM.

4. Discussion

Cervical disc herniation is most common at the C5-6 level [28, 29]. Patients with symptomatic CDH have good outcome associated with performing CSM at the level of CDH with a proper procedure [30]. However, some reports showed that the aggravating cervical disc rupture can occur during a process of CSM with preexisting CDH under the wrong manipulative technique [31–33]. Rotation to the healthy side and lateral bending to the herniated side are the common procedures on performing CSM for CDH. However, the rationality of lateral bending position is still unclear. In the present study, we first established a three-dimensional FE model of CDH (herniated on the left side) in C5-6 segment and then simulated CSM with rotation to the right side (healthy side) in left lateral bending position, neutral position, and right lateral bending position, to evaluate the biomechanical effects of lateral bending to the herniated side during CSM.

Interestingly, we discovered that the cervical disc on herniated side moved forward during CSM in three positions: left lateral bending, neutral, and right lateral bending. The maximum forward displacements were 0.23, 0.36, and 0.45 mm in left lateral bending position, neutral position, and right lateral bending position, respectively. It reconfirmed that a small forward displacement was generated on the left side (herniated side) of cervical disk during CSM with rotation to the right side (healthy side) [12]. And it suggested that the release of the compressed nerve root from small to large was left lateral bending position, neutral position, and right lateral bending position, respectively.

As the same trend of cervical disc displacement, the annulus fiber stress on herniated side from small to large was 7.40, 16.39, and 22.75 MPa in left lateral bending position, neutral position, and right lateral bending position, respectively. The distribution of annulus fibers stress in left lateral bending position was concentrated on the right posterolateral (the healthy side) of cervical disc. However, the distributions of annulus fibers stress in neutral position and right lateral bending position were both concentrated on the left posterolateral (the herniated side) of cervical disc. It suggested that the cervical disc on herniated side may be damaged with high annulus fibers stress during CSM [34]. Compared with left lateral bending position, the left posterolateral (the herniated side) of cervical disc is acted upon by larger annulus fibers stresses in neutral position and right lateral bending position, which could lead to more torsion on the herniated side and, thereby, may cause injury to the annulus fibrosus and therefore aggravation of disc herniation [35–37].

The facet joint stress showed the opposite trend with the cervical disc displacement and annulus fiber stress. The maximum facet stresses at left superior cartilage of C6 were 6.88, 3.60, and 0.12 MPa in left lateral bending position, neutral position, and right lateral bending position during CSM. As we known, facet joints play a role in mechanical function which can contribute to spinal stability and load sharing between spinal columns [38, 39]. The left facet joint contact force increased in left lateral bending position, rotation to right side, or combination of both [40]. In the present study, facet stress at left superior cartilage of C6 was largest in left bending position during CSM with rotation to right side. In other words, the left facet joint sustained more mechanical loading and restricted segmental motion in left bending position with rotation to right side.

Compared with neutral position and right lateral bending position, though the forward displacement of cervical disc on herniated side was smaller in left lateral bending position, the annulus fiber stress on herniated side was declined by sharing load on the left facet joint. In summary, lateral bending to the herniated side on performing CSM for treating CDH tends to protect the cervical disc on herniated side during such manipulation. Future clinical studies are needed to verify the biomechanical effects of lateral bending position on performing CSM for treating CDH.

5. Limitations

The present study has certain shortcomings. First, the current results were based on the FE model of a healthy female who did not have cervical disc herniation. The simulation of cervical disc herniation in C5-6 segment was made by the model adjustment, but other parameters and geometry did not change from the normal FE model. The interpretation of the results should be cautious because the results were drawn from the adjusted model. Second, the loading conditions were highly idealized and could not represent the complicated condition of cervical spinal manipulation. Therefore, it should be kept in mind that the present results were driven by these assumptions. However, the results of FE analysis represent trends rather than precise values because of the necessary simplifications and assumptions concerning the geometry, material properties of the different tissues, contact behavior, and applied loads [41]. Third, the musculature’s effect on the stability of the cervical spinal was not considered in the present study. It remains to be determined how much the current findings would vary if this limitation is appropriately addressed in our future FE models.

6. Conclusions

Compared with neutral position and right lateral bending position, though the forward displacement of cervical disc on herniated side was smaller in left lateral bending position, the annulus fiber stress on herniated side was declined by sharing load on the left facet joint. In summary, lateral bending to the herniated side on performing CSM for treating CDH tends to protect the cervical disc on herniated side. Future clinical studies are needed to verify that.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Disclosure

Xuecheng Huang is the first author. Linqiang Ye, Zixian Wu, and Lichang Liang are the co-first authors.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This research was supported by the National Natural Science Foundation of China (Grant no. 81273871), the Application and Promotion Project of Characteristic Diagnosis and Treatment Technique of South of the Five Ridges of First Affiliated Hospital of Guangzhou University of Chinese Medicine (Grant no. 2016JY03), and the Science and Technology Planning Project of Guangdong Province (2015A020214013). The authors thank Jingmin Huang for her help in language.

References

J. J. Wong, P. Côté, J. J. Quesnele, P. J. Stern, and S. A. Mior, “The course and prognostic factors of symptomatic cervical disc herniation with radiculopathy: A systematic review of the literature,” The Spine Journal, vol. 14, no. 8, pp. 1781–1789, 2014.

B. I. Woods and A. S. Hilibrand, “Cervical radiculopathy: epidemiology, etiology, diagnosis, and treatment,” Journal of Spinal Disorders & Techniques, vol. 28, no. 5, pp. E251–E259, 2015.

G. Bronfort, M. Haas, R. L. Evans, and L. M. Bouter, “Efficacy of spinal manipulation and mobilization for low back pain and neck pain: A systematic review and best evidence synthesis,” The Spine Journal, vol. 4, no. 3, pp. 335–356, 2004.

J. Herzog, “Use of cervical spine manipulation under anesthesia for management of cervical disk herniation, cervical radiculopathy, and associated cervicogenic headache syndrome,” Journal of Manipulative and Physiological Therapeutics, vol. 22, no. 3, pp. 166–170, 1999.

X. Wei, S. Wang, J. Li et al., “Complementary and alternative medicine for the management of cervical radiculopathy: an overview of systematic reviews,” Evidence-Based Complementary and Alternative Medicine, vol. 2015, Article ID 793649, 10 pages, 2015.

J. Fernández-Carnero, C. Fernández-de-las-Peñas, and J. A. Cleland, “Immediate hypoalgesic and motor effects after a single cervical spine manipulation in subjects with lateral epicondylalgia,” Journal of Manipulative and Physiological Therapeutics, vol. 31, no. 9, pp. 675–681, 2008.

C. K. Peterson, C. W. A. Pfirrmann, J. Hodler et al., “Symptomatic, Magnetic Resonance Imaging-Confirmed Cervical Disk Herniation Patients: A Comparative-Effectiveness Prospective Observational Study of 2 Age- and Sex-Matched Cohorts Treated with Either Imaging-Guided Indirect Cervical Nerve Root Injections or Spinal Manipulative Therapy,” Journal of Manipulative and Physiological Therapeutics, vol. 39, no. 3, pp. 210–217, 2016.

C. K. Peterson, C. Schmid, S. Leemann, B. Anklin, and B. K. Humphreys, “Outcomes from magnetic resonance imaging-confirmed symptomatic cervical disk herniation patients treated with high-velocity, low-amplitude spinal manipulative therapy: A prospective cohort study with 3-month follow-up,” Journal of Manipulative and Physiological Therapeutics, vol. 36, no. 8, pp. 461–467, 2013.

R. Van Schalkwyk and G. F. Parkin-Smith, “A clinical trial investigating the possible effect of the supine cervical rotatory manipulation and the supine lateral break manipulation in the treatment of mechanical neck pain: A Pilot Study,” Journal of Manipulative and Physiological Therapeutics, vol. 23, no. 5, pp. 324–331, 2000.

L. G. Zhu, X. Wei, and S. Q. Wang, “Does cervical spine manipulation reduce pain in people with degenerative cervical radiculopathy? A systematic review of the evidence, and a meta-analysis,” Clinical Rehabilitation, 2015.

D. W. Evans, “Mechanisms and effects of spinal high-velocity, low-amplitude thrust manipulation: previous theories,” Journal of Manipulative and Physiological Therapeutics, vol. 25, no. 4, pp. 251–262, 2002.

L.-P. Wu, Y.-Q. Huang, D. Manas, Y.-Y. Chen, J.-H. Fan, and H.-G. Mo, “Real-time monitoring of stresses and displacements in cervical nuclei pulposi during cervical Spine manipulation: A finite element model analysis,” Journal of Manipulative and Physiological Therapeutics, vol. 37, no. 8, pp. 561–568, 2014.

J. G. Pickar, “Neurophysiological effects of spinal manipulation,” The Spine Journal, vol. 2, no. 5, pp. 357–371, 2002.

Z. J. Mo, Y. B. Zhao, L. Z. Wang, Y. Sun, M. Zhang, and Y. B. Fan, “Biomechanical effects of cervical arthroplasty with U-shaped disc implant on segmental range of motion and loading of surrounding soft tissue,” European Spine Journal, vol. 23, no. 3, pp. 613–621, 2014.

Z. Mo, Y. Zhao, C. Du, Y. Sun, M. Zhang, and Y. Fan, “Does Location of Rotation Center in Artificial Disc Affect Cervical Biomechanics?” The Spine Journal, vol. 40, no. 8, pp. E469–E475, 2015.

G. Denozière and D. N. Ku, “Biomechanical comparison between fusion of two vertebrae and implantation of an artificial intervertebral disc,” Journal of Biomechanics, vol. 39, no. 4, pp. 766–775, 2006.

H.-J. Kim, H.-J. Chun, H.-M. Lee et al., “The biomechanical influence of the facet joint orientation and the facet tropism in the lumbar spine,” The Spine Journal, vol. 13, no. 10, pp. 1301–1308, 2013.

H.-G. Du, S.-H. Liao, Z. Jiang et al., “Biomechanical analysis of press-extension technique on degenerative lumbar with disc herniation and staggered facet joint,” Saudi Pharmaceutical Journal, vol. 24, no. 3, pp. 305–311, 2016.

D. Ganbat, Y. H. Kim, K. Kim, Y. J. Jin, and W. M. Park, “Effect of mechanical loading on heterotopic ossification in cervical total disc replacement: a three-dimensional finite element analysis,” Biomechanics and Modeling in Mechanobiology, vol. 15, no. 5, pp. 1191–1199, 2016.

S. Lee, Y. Im, K. Kim, Y. Kim, W. Park, and K. Kim, “Comparison of Cervical Spine Biomechanics After Fixed- and Mobile-Core Artificial Disc Replacement,” The Spine Journal, vol. 36, no. 9, pp. 700–708, 2011.

D. Liang, L.-Q. Ye, X.-B. Jiang et al., “Biomechanical effects of cement distribution in the fractured area on osteoporotic vertebral compression fractures: A three-dimensional finite element analysis,” Journal of Surgical Research, vol. 195, no. 1, pp. 246–256, 2015.

T. Xie, J. Qian, Y. Lu, B. Chen, Y. Jiang, and C. Luo, “Biomechanical comparison of laminectomy, hemilaminectomy and a new minimally invasive approach in the surgical treatment of multilevel cervical intradural tumour: A finite element analysis,” European Spine Journal, vol. 22, no. 12, pp. 2719–2730, 2013.

S. P. Moroney, A. B. Schultz, J. A. A. Miller, and G. B. J. Andersson, “Load-displacement properties of lower cervical spine motion segments,” Journal of Biomechanics, vol. 21, no. 9, pp. 769–779, 1988.

M. M. Panjabi, J. J. Crisco, A. Vasavada et al., “Mechanical properties of the human cervical spine as shown by three-dimensional load-displacement curves,” The Spine Journal, vol. 26, no. 24, pp. 2692–2700, 2001.

M. Hussain, R. N. Natarajan, G. Chaudhary, H. S. An, and G. B. J. Andersson, “Posterior facet load changes in adjacent segments due to moderate and severe degeneration in C5-C6 disc: A poroelastic C3-T1 finite element model study,” Journal of Spinal Disorders & Techniques, vol. 25, no. 4, pp. 218–225, 2012.

M. Hussain, R. N. Natarajan, H. S. An, and G. B. J. Andersson, “Patterns of height changes in anterior and posterior cervical disc regions affects the contact loading at posterior facets during moderate and severe disc degeneration: A poroelastic C5-C6 finite element model study,” The Spine Journal, vol. 35, no. 18, pp. E873–E881, 2010.

K. R. Wade, P. A. Robertson, A. Thambyah, and N. D. Broom, “How healthy discs herniate: A biomechanical and microstructural study investigating the combined effects of compression rate and flexion,” The Spine Journal, vol. 39, no. 13, pp. 1018–1028, 2014.

E. Okada, M. Matsumoto, D. Ichihara et al., “Aging of the cervical spine in healthy volunteers: A 10-year longitudinal magnetic resonance imaging study,” The Spine Journal, vol. 34, no. 7, pp. 706–712, 2009.

M. Teraguchi, N. Yoshimura, H. Hashizume et al., “Prevalence and distribution of intervertebral disc degeneration over the entire spine in a population-based cohort: The Wakayama Spine Study,” Osteoarthritis and Cartilage, vol. 22, no. 1, pp. 104–110, 2014.

M. J. Hubka, S. P. Phelan, P. M. Delaney, and V. L. Robertson, “Rotary manipulation for cervical radiculopathy: Observations on the importance of the direction of the thrust,” Journal of Manipulative and Physiological Therapeutics, vol. 20, no. 9, pp. 622–627, 1997.

D. G. Malone, N. G. Baldwin, F. J. Tomecek et al., “Complications of cervical spine manipulation therapy: 5-year retrospective study in a single-group practice.,” Neurosurgical Focus, vol. 13, no. 6, p. ecp1, 2002.

L. Padua, R. Padua, M. LoMonaco, and P. A. Tonali, “Radiculomedullary complications of cervical spinal manipulation,” Spinal Cord, vol. 34, no. 8, pp. 488–492, 1996.

S.-H. Tseng, S.-M. Lin, Y. Chen, and C.-H. Wang, “Ruptured cervical disc after spinal manipulation therapy: report of two cases.,” The Spine Journal, vol. 27, no. 3, pp. E80–82, 2002.

K. Wang, H. Wang, Z. Deng, Z. Li, H. Zhan, and W. Niu, “Cervical traction therapy with and without neck support: A finite element analysis,” Musculoskeletal Science and Practice, vol. 28, pp. 1–9, 2017.

P. J. Fazey, S. Song, Å. Mønsås et al., “An MRI investigation of intervertebral disc deformation in response to torsion,” Clinical Biomechanics, vol. 21, no. 5, pp. 538–542, 2006.

J. D. M. Drake, C. D. Aultman, S. M. McGill, and J. P. Callaghan, “The influence of static axial torque in combined loading on intervertebral joint failure mechanics using a porcine model,” Clinical Biomechanics, vol. 20, no. 10, pp. 1038–1045, 2005.

S. P. Veres, P. A. Robertson, and N. D. Broom, “The influence of torsion on disc herniation when combined with flexion,” European Spine Journal, vol. 19, no. 9, pp. 1468–1478, 2010.

X. Rong, Z. Liu, B. Wang, H. Chen, and H. Liu, “The facet orientation of the subaxial cervical spine and the implications for cervical movements and clinical conditions,” The Spine Journal, vol. 42, no. 6, pp. E320–E325, 2017.

X. Rong, B. Wang, C. Ding et al., “The biomechanical impact of facet tropism on the intervertebral disc and facet joints in the cervical spine,” The Spine Journal, 2017.

X. Liu, Z. Huang, R. Zhou et al., “The effects of orientation of lumbar facet joints on the facet joint contact forces: an in vitro biomechanical study,” The Spine Journal, p. 1, 2017.

A. Rohlmann, T. Zander, and G. Bergmann, “Comparison of the biomechanical effects of posterior and anterior spine-stabilizing implants,” European Spine Journal, vol. 14, no. 5, pp. 445–453, 2005.