Periprosthetic bone loss following total hip arthroplasty (THA) is a serious concern leading to the premature failure of prosthetic implant. Therefore, investigating bone remodeling in response to hip arthroplasty is of paramount for the purpose of designing long lasting prostheses. In this study, a thermodynamic-based theory, which considers the coupling between the mechanical loading and biochemical affinity as stimulus for bone formation and resorption, was used to simulate the femoral density change in response to THA. The results of the numerical simulations using 3D finite element analysis revealed that in Gruen zone 7, after remarkable postoperative bone loss, the bone density started recovering and got stabilized after 9% increase. The most significant periprosthetic bone loss was found in Gruen zone 7 (−17.93%) followed by zone 1 (−13.77%). Conversely, in zone 4, bone densification was observed (+4.63%). The results have also shown that the bone density loss in the posterior region of the proximal metaphysis was greater than that in the anterior side. This study provided a quantitative figure for monitoring the distribution variation of density throughout the femoral bone. The predicted bone density distribution before and after THA agree well with the bone morphology and previous results from the literature.