Background. Vertebral column decancellation (VCD) is a new spinal osteotomy technique to correct thoracolumbar kyphotic deformity (TLKD). Relevant biomechanical research is needed to evaluate the safety of the technique and the fixation system. We aimed to develop an accurate finite element (FE) model of the spine with TLKD following VCD and to provide a reliable model for further biomechanical analysis. Methods. A male TLKD patient who had been treated with VCD on L2 and instrumented from T10 to L4 was a volunteer for this study. The CT scanning images of the postoperative spine were used for model development. The FE model, simulating the spine from T1 to the sacrum, includes vertebrae, intervertebral discs, spinal ligaments, pedicle screws, and rods. The model consists of 509580 nodes and 445722 hexahedrons. The ranges of motion (ROM) under different loading conditions were calculated for validation. The stresses acting on rods, screws, and vertebrae were calculated. Results. The movement trend, peak stress, and ROM calculated by the current FE model are consistent with previous studies. The FE model in this study is able to simulate the mechanical response of the spine during different motions with different loading conditions. Under axial compression, the rod was the part bearing the peak stress. During flexion, the stress was concentrated on proximal pedicle screws. Under extension and lateral bending, an osteotomized L1 vertebra bore the greatest stress on the model. During tests, ligament disruption and unit deletion were not found, indicating an absence of fracture and fixation breakage. Discussion. A subject-specific FE model of the spine following VCD is developed and validated. It can provide a reliable and accurate digital platform for biomechanical analysis and surgical planning.