Objective To establish the finite element model of L1~L5 lumbar motion segments with contact theory by self-compiled software and HyperMesh. Method Based on the spine CT continuous images of L1~L5 segments from the healthy volunteer, the three dimensional finite elements model of L1~L5 vertebrae structure was constructed applying self compiled software and HyperMesh, and added with simulating the endplates, joint cartilage intervertebral disc as well as ligament structure. The contact theory was adopted to solve the process between the high and low articular surface surface contacting. ResultsThe three dimensional finite element model with more structure integrity and good geometric similarity was obtained. The simulating results of the physical loading are matched better with the results of the vitro experiment of biomechanics. Conclusions Selfcompiled software and HyperMesh improved the speed and accuracy for building finite element model of lumbar motion segment, while the introduction of contact theory could enhance its reliability and authenticity.