Bone remodeling process is investigated theoretically and numerically within the framework of the extended adaptive elasticity. The coupling between the internal and surface remodeling is considered for the solution of the bone remodeling process. A semi-analytical solution based on the state-space method was used to analyze remodeling process of inhomogeneous bone materials subjected to multi-field loadings. The results show that, in contrast to the former works without considering the coupling between the internal and the surface remodeling, the magnitude of the electrical field changes is obviously lower and the remodeling time is shorten significantly. This implies that the coupling between the internal and the surface remodeling plays an important role in the bone remodeling process