The process of bone remodeling is regulated by metabolic activities of many bone cells. While osteoclasts and osteoblasts are responsible for bone resorption and formation, respectively, activities of these cells are believed to be controlled by a mechanosensory system of osteocytes embedded in the extracellular bone matrix. Several experimental and theoretical studies have suggested that the strain-derived interstitial fluid flow in lacuno-canalicular porosity serves as the prime mover for bone remodeling. Previously, we constructed a mathematical model for trabecular bone remodeling that interconnects the microscopic cellular activities with the macroscopic morphological changes in trabeculae through the mechanical hierarchy. This model assumes that fluid-induced shear stress acting on osteocyte processes is a driving force for bone remodeling. The validity of this model has been demonstrated with a remodeling simulation using a two-dimensional trabecular model. In this study, to investigate the effects of loading frequency, which is thought to be a significant mechanical factor in bone remodeling, we simulated morphological changes of a three-dimensional single trabecula under cyclic uniaxial loading with various frequencies. The results of the simulation show the trabecula reoriented to the loading direction with the progress of bone remodeling. Furthermore, as the imposed loading frequency increased, the diameter of the trabecula in the equilibrium state was enlarged by remodeling. These results indicate that our simulation model can successfully evaluate the relationship between loading frequency and trabecular bone remodeling. (C) 2011 Elsevier Ltd. All rights reserved.