paper proposes a new method of transferable calculation of elastic constants of the second, third and fourth orders. The method is based on the assumption that the potential energy of bulk crystals may be presented as the decomposition into irreducible energies of two, three, and four atomic clusters. To justify this decomposition, we proved that the irreducible energy of atomic clusters obeys additional symmetry O3×NCl, which depends on the number of atoms constituting clusters (NCl), and not depends on cluster structure. Additional symmetry allows defining the bases invariants, which depend on interatomic distances, and irreducible energies of clusters are depended on bases invariants only. In turn, these dependences made it possible to calculate the elastic modulus of complexes that exist in different crystal modifications as functions of the same phenomenological parameters of the model of interatomic interactions. Then we presented an example of the calculation of elastic constants BCC, HCP and FCC of the metal phases of Co polymorphs. The results confirm that the proposed method is transferable and makes it possible to calculate the elastic moduli of the second, third and fourth orders with the accuracy not worse than was achieved using the quantum chemistry models.