Vane viscometers are often used to investigate the low shear rate properties of plastic fluids. The shear stress is determined by assuming that the material is held in the space between the vane blades so that it behaves like a rigid cylinder. Experimental evidence supports this assumption and the aim of the present study is to model numerically the yield process in a vane rheometer using viscoelastic and plastic fluids. The finite element method has been used to model the behaviour of Herschel-Bulkley (Bingham), Casson and viscoelastic (Maxwell type) fluids. The penalty function approach for the pressure approximation and a rotating reference frame are used together with fine meshes containing more than 1300 elements. The results show that for Herschel-Bulkley (Bingham), and Casson fluids a rotating rigid cylinder of fluid is trapped inside the periphery of the vane, the shear stress is uniformly distributed over the surface of the cylinder. Finally a modified second order fluid is used to simulate the viscoelastic behaviour, anticipated to be an intermediate between the elastic deformation and the plastic flow, to provide a more realistic simulation of the yield process about a vane. In this case, contrast with the concentration of the elastic strain rate at the blade tips, a nearly uniform distribution of the plastic shear rate is still found. This implies that the plastic shear always distributes uniformly during the entire yielding process. Evidently the assumption of uniform shear on a rotating cylinder of material occluded in the blades of a vane is a valid and useful model for many types of fluid possessing a yield stress.