Herein, the bending behavior of a hot‐rolled hyper‐sandelin steel subjected to hot‐dip galvanizing (HDG) is numerically simulated. Two coating types are analyzed, that is, a pure zinc coating and a zinc‐based coating with 3% Sn addition by weight. It is experimentally observed that the tin addition influences the thickness of each intermetallic phase and increases the ductility of the intermetallic phases. Therefore, the key of the numerical simulation is to assume suitable constitutive laws to phase layers, by considering that only few mechanical parameters are available in the literature. Such an approach has shown quite satisfactory results and it can be promising as a numerical tool to be used by HDG industry, being quite simple but quite accurate to estimate flexural strength and damage of Zn‐based coated steels.