In this study, a biological process of bone remodeling was considered as a closed loop feedback control system, which enables bone to optimize and renew itself over a lifetime. A novel idea of combining strain-adaptive and damage-induced remodeling algorithms at Basic Multicellular Unit (BMU) level was introduced. In order to make the outcomes get closer to clinical observation, the stochastic occurrence of microdamage was involved and a hypothesis that remodeling activation probability is related to the value of damage rate was assumed. Integrated with Finite Element Analysis (FEA), the changes of trabecular bone in morphology and material properties were simulated in the course of five years. The results suggest that deterioration and anisotropy of trabecluar bone are inevitable with natural aging, and that compression rather than tension can be applied to strengthen the ability of resistance to fracture. This investigation helps to gain more insight the mechanism of bone loss and identify improved treatment and prevention for osteoporosis or stress fracture.