Objective To establish a three-dimensional(3D) finite element model of cervical vertebrae (C1-7), and study its biomechanical properties under muscle force by cervical traction, so as to provide references for clinical treatment. Methods On the basis of nonlinear finite element model of normal cervical vertebrae and combined with clinical traction methods, cervical traction at the extension angle of 0°, 10°, 20°, 30°, 40° under the same traction weight, was simulated by finite element analysis (FEA) software to obtain and select the joint force and muscle force that were appropriate for FEA on the model. Results In the process of cervical extension by traction, under the muscle force, the average maximum equivalent stress of cervical vertebrae, intervertebral disc and uncovertebral joints increased by 4.86, 1.79, 0.69 MPa, respectively, and the average maximum relative displacement of cervical vertebrae in sagittal and vertical axis direction increased by 11.1, 1.26 mm, respectively. The biomechanical properties of cervical traction were similar to the FEA results reported in the literature. Conclusions Neck muscles play an active role in promoting the stress and displacement of cervical vertebrae, intervertebral discs and uncovertebral joints and it should be taken into consideration when performing cervical traction in clinic. In addition, the traction angle should not be too large: 0°-20° is generally recommended as a relatively safe angle range at the initial stage.