Objective To study dynamic characteristics of human lumbar spine using three-dimensional finite element method. Methods Finite element model of lumbar spine (L1~5) was developed and validated based on CT images, and the modal analysis was also conducted. Results A total of top 30-order modal parameters were extracted to obtain dynamic characteristics of the lumbar spine under free boundary conditions. Resonance frequencies of the model were concentrately distributed, but the amplitude of each order varied greatly. Amplitude near L5 segment was much larger, indicating L5 was easily to be injured. This lumbar modal analysis could provide a basis for its further dynamic analysis. Parameters such as natural frequency, modal shape and vibration amplitude of the lumbar spine would be helpful for both lumbar dynamic analysis and optimal design of man-machine interface mechanical equipment.