A new method for analyzing the bone-remodeling process using quasi-molecular dynamics was proposed in this study. The effect of pressure due to bone marrow, which could not be considered in previous methods, was also considered in this method. Bone-remodeling response of the 2D vertebral body of lumbar spine to a uniaxial compressive displacement of 1.8564mm. corresponding to approximately 2kN of compressive load, was studied. Converged shape change of the cortical shell and rearrangement of cancellous bone structure matched well with a normal shape of the vertebral body. The calculated responses in the spinal elements also shows closed results compared with experimental results.