Understanding the biomechanical effects of cervical disc degeneration (CDD) on the cervical spine is fundamental for understanding the mechanisms of spinal disorders and improving clinical treatment. While the biomechanical effects of CDD on segmental flexibility and the posterior facets have been reported, a clear understanding of the effect of the motion loading method on facet joint forces after CDD is still lacking. Therefore, the objective of this study was to determine the effect of the motion loading method on facet joint forces after CDD. A three-dimensional nonlinear finite element (FE) model of the cervical spine (C3-C7) was developed and validated to represent normal conditions. This normal model was modified to create six degenerative models simulating mild, moderate, and severe grades of disc degeneration at C5-C6. While under a follower compressive preload (73.6 N), a 1-Nm moment was applied to all models to determine range of motion (ROM). A displacement load was applied to all degenerative models under the same follower load, making the C5-C6 degeneration segment motion same to the ROM of C5-C6 in normal model, and facet joint forces were computed. Compared with the normal model, ROM of the C5-C6 degenerative segments dramatically declined in all postures with increasing degenerative pathologies in the disc. The ROM in the adjacent normal segments of the degenerative segments also declined, with the exception of C4-C5 during extension. Under a 1-Nm moment load, there were not obvious changes in facet joint forces in the C5-C6 degenerative segment with increasing grades of degeneration, but facet joint forces in the adjacent normal segments did increase. Under a displacement load, the facet joint forces of the C5-C6 degenerative segment increased with increasing grades of degeneration. Facet joint forces were positively correlated with the ROM of the degenerative segment, demonstrating that the motion loading method had a significant effect on facet joint forces after CDD. Loading conditions must be strictly controlled in future finite element analysis studies to improve the comparability between models built by different units. Copyright © 2019 Elsevier Ltd. All rights reserved.