To investigate the bone remodeling phenomenon around implant device, 3-D mathematical simulation model was developed. Strain energy density from the finite element method was chosen for the indicator for remodeling process. Recursive calculations continued until converged results between FEM and mathematical model. For a osteo-integration example, bone-remodeling process in a implanted tibia of beagle was adapted. Calculated results indicated that the bone densities around screw pitch were increased which indicates firm fixations between the bone and implant. Screw design parameters have an influence on initial stability of the implant rather than remodeling process.