摘要Background Thoracolumbar burst fracture is a common clinical injury,and the fracture mechanism is still controversial.The aim of this research was to study the formation of intracanal fracture fragments in thoracolumbar burst fractures and to provide information for the prevention of thoracolumbar bursts fractures and reduction of damage to the nervous system.Methods A nonlinear three-dimensional finite element model of T11-L3 segments was established,and the injury processes of thoracolumbar bursts were simulated.The intact finite element model and the finite element model after the superior articular were impacted by 100 J of energy in different directions.The distribution and variation of stress in the superior posterior region of the L1 vertebral body were analyzed.Abaqus 6.9 explicit dynamic solver was used as finite element software in calculations.Results A three-dimensional nonlinear finite element model of the thoracolumbar spine was created.In the intact model,stress was concentrated in the superior posterior region of the L1 vertebral body.The stress peak was a maximum for the extension impact load and a minimum for the flexion impact load.The stress peak and contact force in the facet joint had close correlation with time.The stress peak disappeared after excision of the superior articular process.Conclusions The three-dimensional nonlinear finite element model was suitable for dynamic analysis.The contact force in the facet joint,which can be transferred to the superior posterior vertebral body,may explain the spinal canal fragment in thoracolumbar burst fractures.