Bone is a dynamic living tissue that undergoes continuous adaptation of its mass and structure in response to mechanical and biological environment demands. Studies of bone adaptation have focused on metabolic or mechanical stimulus, but mathematical models of bone adaptation considering both, are not available by now. In this paper, we propose a mathematical model of bone adaptation during a remodeling cycle due to mechanical stimulus with the introduction of osteocytes as mechanotransducers. The model captures qualitatively very well the bone adaptation and cell interactions during the bone remodeling.