Harmonic quasiconformal mappings of Riemannian manifolds