This paper focusses on computational aspects related to a recently proposed anisotropic viscoelastic model for soft biological tissues at large strains [1]. A key aspect of this model is the generalisation of micromechanically motivated one‐dimensional constitutive equations to three dimensions by numerical integration over the unit sphere. A strong effect of this procedure on the accuracy and in particular on the material symmetry of the model is observed. Finally a finite element example of an artery subject to normotensive blood pressure is presented. (© 2009 Wiley‐VCH Verlag GmbH & Co. KGaA, Weinheim)