This approach has been tested on some real systems, receiving interesting confirmation to its viability [@Disc1; @Disc2]. However, being derived via a phenomenological approach, the Lichnerowicz energy momentum tensor must be completed by the specification of an ansatz linking the fluid index to the thermodynamical variables of the system, so closing the dynamical problem. In what follows, we implement the Lichnerowicz treatment to the viscous Bianchi I cosmology, by pursuing two different tasks. On one hand, we study the solution with matter creation, fixing the fluid index via the request of incompressibility. Results show that the Universe evolves trough an intrinsic shear-driven anisotropic solution, meaning that also in the Lichnerowicz scenario the solution with matter creation exists and, actually, such a phenomenon is enhanced, being therefore not related to non-physical effects of the Eckart formulation. On the other hand, we analyze the isotropic limit near the singularity, by reducing the three scale factors of the Bianchi I model to be equal. The latter is known in literature as the flat Robertson-Walker Universe, for which only bulk viscosity may be relevant due to the homogeneity and isotropy of the model, preventing shear among different layers. In this specific study we see 