Objective Aiming at the low accuracy problem of the traditional finite element method ( FEM) in analyzing the piezoelectric femur remodeling after total hip arthroplasty (THA), the edge-based smoothed finite element method ( ES-FEM ) was applied to simulate remodeling process of the piezoelectric femur after prosthesis implantation. Methods According to the adaptive bone remodeling theory, the prosthesis-piezoelectric femur model was established. Based on background mesh of the model, the edge-based smoothing domains were constructed. Then, the smoothed remodeling stimuli were calculated by introducing the gradient smoothing technology. The distribution of apparent density in proximal femur could be obtained. Results After prosthesis implantation, the stress points were transferred from the femoral head to the prosthesis. Hence, the stress shielding phenomenon occurred. The apparent density distributions changed obviously. Compared with the FEM the ES-FEM could soften the numerical model and improve the simulation accuracy. Based on the same mesh, the computation accuracy of the potential and density was improved by about 27% and 30% , respectively. Conclusions The ES-FEM in this study possesses better accuracy in analyzing the bone remodeling after THA, which provides an effective theoretical basis for clinical study of artificial joint replacement.