This paper outlines a computational model for the analysis of the piezoelectric behaviour of the vertebral body remodelling process. Particular attention is paid to the algorithms for the simulation of the stress energy density for each point of the geometry and the distribution of the density in the bone. In addition, the model takes into account the piezoelectric effect and the anisotropy (transversal isotropy) of the bone. A model for internal anisotropic piezoelectric bone remodelling of a human vertebra is discussed in detail. The model consists of the implementation of an algorithm which includes the elastic and electric variables in a single equation using boundary element method. The presented results show a good agreement with biological data and the model does not include any electric additional charge.