Twinning is a major deformation mechanism in various materials, especially when few dislocation slip systems are operative. It is the case of zinc-rich coatings in galvanised steel sheets, made of pancake grains on a substrate and where the slip systems with a non-vanishing component along the c-axis present high critical resolved shear stress values. In addition, the abrupt lattice orientation change associated to twinning, the stress relaxation during its propagation and the localised nature of its early stages make it difficult to reproduce this deformation mechanism by using classical crystal plasticity models conceived for dislocation slip. In this sense, this work proposes a hierarchy of three twinning models in combination with a dislocation slip crystal plasticity model, for the case of a ZnAlMg coating. These three models are: a relaxed-Taylor model applied to individual crystal orientations of the coating, a “pseudo-slip” model for twinning and a localised twinning model. The latter incorporates a linear softening in the material law accounting for the unstable twinning initiation and enforces twinning lattice reorientation. A microstructure portion extracted from an in-situ SEM tensile experiment on galvanised steel is used to perform 3D full-field finite element simulations within a finite strain formulation. SEM observations and EBSD acquisitions are used to compare simulation and experimental results during the different steps of the in-situ SEM test, regarding the deformation and damage modes of the zinc-rich coating. The focus is set on twinning evolution inside some individual grains, and the pros and the cons of the three models are finally discussed.