Studying the vibration behavior of the entire spine can guide the design of seat comfort and vibration safety. However, due to simplification of traditional biomechanical models, very few studies have analyzed in detail the vibration behavior characteristics of the entire spine inside a seated human body. Therefore, this study aimed to provide guidance and reference for spinal modeling and biomechanical research in ergonomics. A developed finite element model of three-dimensional seated human body was validated and adjusted based on anatomical data of human spine in detail. Static analysis, modal analysis and random response analysis (under vertical excitation between 0 and 20 Hz at 1 m/s2 r.m.s.) were conducted. The range of motion, modal frequencies and the tri-axial transmissibility of the developed models matched well with experimental results. In the vertical resonance mode, the entire spine contained both vertical deformation (60% of the total) and vertical displacement (40%), in addition, the cervical spine, especially the lumbar spine, also contained bending deformation which could alleviate the impact, but led to complex alternating stresses, increasing the risks of the spinal injuries under vertical whole-body vibration. From the bottom to the top of the spine, the frequency distribution of vertical transmissibility became steeper, the peak value increased, and the number of peaks decreased. This study provided new insights into the vibration behavior and frequency response of the entire spine inside a seated human body for improving seat comfort and vibration safety.