The piezoelectricity of bone is known to play a crucial role in bone adaptation and remodeling. The application of an external stimulus such as mechanical strain or electric field has the potential to enhance bone formation and implant osseointegration. Therefore, in the present study, the objective is to investigate bone remodeling under electromechanical stimulation as a step towards establishing therapeutic strategies. For the first time, piezoelectric bone remodeling in the human proximal tibia under electro-mechanical loads was analyzed using the finite element method in an open-source framework. The predicted bone density distributions were qualitatively and quantitatively assessed by comparing with the computed tomography (CT) scan and the bone mineral density (BMD) calculated from the CT, respectively. The effect of model parameters such as uniform initial bone density and reference stimulus on the final density distribution was investigated. Results of the parametric study showed that for different values of initial bone density the model predicted similar but not identical final density distribution. It was also shown that higher reference stimulus value yielded lower average bone density at the final time. The present study demonstrates an increase in bone density as a result of electrical stimulation. Thus, to minimize bone loss, for example, due to physical impairment or osteoporosis, mechanical loads during daily physical activities could be partially replaced by therapeutic electrical stimulation.