A method for the calculation of translationally invariant wave functions for systems of identical fermions with arbitrary potential of pair interaction is developed. It is based on the well-known result that the essential dynamic part of the Hamiltonian for the system of identical particles is the Reduced Hamiltonian operator describing relative movement of two particles inside the system. The eigenfunctions of this operator take into account all correlations caused by interaction. These eigenfunctions are basic for the construction of the components (i.e. the functions with a lower degree of antisymmetry) of the system wave functions. The main problem of this approach appears to be antisymmetrization of the components. The developed universal algorithm for antisymmetrization gives a possibility to perform this operation in a simple way and keep numerical approximations under control.