The liquid crystals anchoring theory is considered on the molecular level. Earlier calculations of mesogenic molecules interaction energy with graphite and PE surface give results being one-two orders higher than experimentally obtained values. To explain the experimental values of anchoring energy we proposed its dependence on the interaction between liquid crystal layers. In our calculations we have simulated the mesogen molecules by virtual rod-like molecules. The molecules orientation was given by polar and azimuthal angles towards the directors of neighboring liquid crystal layers. Obtained values of polar and azimuthal anchoring energy as well as their dependences on the order parameter are in good agreement with experimental data.