After anterior interbody fusion in the lumbar spine, the accelerated degeneration of the disc adjacent to the fusion levels was clinically observed. To understand the stress distribution of the adjacent disc, this study created a finite element model of the lumbar spine from L1-L5 vertebral body. The fusion model modified from the intact model was used to simulate the anterior interbody fusion. Various loading conditions, which included flexion, extension, lateral bending, and torsion, were applied to the finite element model to study the corresponding stress distribution. From the finite element model calculation, at a lower fusion site or more fusion levels, the stress of the disc adjacent to interbody fusion increased more than upper fusion site or single fusion level under flexion, torsion and lateral bending. Larger stress increase was estimated at the upper disc adjacent to interbody fusion than the lower disc adjacent to interbody fusion. In stress distribution, the upper disc adjacent to interbody fusion had a little alteration under torsion.