Herein, the bending behavior of hot‐dip galvanizing hot‐rolled hypersandelin plates, together with the corresponding crack pattern, is numerically investigated by using a hybrid model developed in Ansys LS‐DYNA environment, by combining the discrete element method (DEM) with finite element method (FEM). The experimental bending tests here simulated and available in the literature are performed by considering two types of bath, that is, a pure zinc bath and a technological bath, consisting in a zinc bath with 3% Sn addition by weight. The results related to the bending behavior are compared with both experimental data and other FE numerical results, previously obtained by some of the present authors. Moreover, the numerical crack patterns are compared with the experimental observations of the longitudinal section of the plates at the end of testing, by means of a light optical microscope. A quite satisfactory agreement is observed.