A section of a Riemannian fibration is called harmonic if it is harmonic as a map into the total space. The Euler-Lagrange equation for the such sections is formulated. This formula is applied to distributions which are sections of Grassmannian bundle and examples of harmonic distributions on Riemannian manifolds are constructed.