The objective of this study was to investigate the process of mandibular bone remodeling induced by implant-supported overdentures. computed tomography (CT) images were collected from edentulous patients to reconstruct the geometry of the mandibular bone and overdentures supported by implants. Based on the theory of strain energy density (SED), bone remodeling models were established using the user material subroutine (UMAT) in abaqus. The stress distribution in the mandible and bone density change was investigated to determine the effect of implant number on the remodeling of the mandibular bone. The results indicated that the areas where high Mises stress values were observed were mainly situated around the implants. The stress was concentrated in the distal neck region of the distal-most implants. With an increased number of implants, the biting force applied on the dentures was almost all taken up by implants. The stress and bone density in peri-implant bone increased. When the stress reached the threshold of remodeling, the bone density began to decrease. In the posterior mandible area, the stress was well distributed but increased with decreased implant numbers. Changes in bone density were not observed in this area. The computational results were consistent with the clinical data. The results demonstrate that the risk of bone resorption around the distal-most implants increases with increased numbers of implants and that the occlusal force applied to overdentures should be adjusted to be distributed more in the distal areas of the mandible.