Trabecular bone’s capability to adapt its architecture to the surrounding environment has been a subject of research for more than a century. The purpose of this research is to simulate and improve a semimechanistic bone remodeling theory. A computer code was developed using finite element method, and also nonuniformly distributed osteocytes scheme was introduced into the model. The regulatory mechanism proposed in the semi-mechanistic model was implemented and studied in a MATLAB code. Also, in addition to simulating the osteocyte mechanosensory role and separation of osteoclast and osteoblast activities as proposed by others, effect of a nonuniform distribution of osteocytes on bone remodelling process was investigated in this study. Trabecular-like structures for different loading conditions and directionalities similar to actual human bone have been resulted for a square plate using this new approach. Because the function of osteocyte cells has been of great interest as a part of the signal transduction system in the bone remodeling, the computation domain was divided into three parts resembling real morphology of spongy bone. It is assumed that number of osteocytes decreases linearly three to one for each element by moving away from the matrix-porosity interface. Results of this research showed that the investigated semi-mechanistic model, despite its simplicity, has a good scientific merit and including varying osteocytes distribution in it will help to predict more realistic morphologies in time-dependent simulations. Ultimate goal of this research is finding a relation between bone diseases and number, also distribution of osteocytes in the bone.