A comprehensive, geometrically accurate, nonlinear C0-T1 three-dimensional finite element (FE) model was developed for the biomechanical study of human cervical spine and related disorders. The model was developed with anatomic detail from the computed tomography (CT) images of a 46-year old female healthy volunteer, and applied the finite element model processing softwares such as MIMICS13.1, Hypermesh11.0, Abaqus 6.12-1, etc., for developing, preprocessing, calculating and analysing sequentially. The stress concentration region and the range of motion (ROM) of each vertebral level under axial rotation, flexion, extension, and lateral bending under physiologic static loadings were observed and recorded. The model was proven reliable, which was validated with the range of motion in previous published literatures. The model predicted the front and side parts of the foramen magnum and contralateral pedicle and facet was the stress concentration region under physiological loads of the upper spine and the lower spine, respectively. The development of this comprehensive, geometrically accurate, nonlinear cervical spine FE model could provide an ideal platform for theoretical biomechanical study of human cervical spine and related disorders.