The present study aims to model the optical response of plasmonic metasurfaces made of a periodic arrangement of metallic particles with arbitrary shape and subwavelength dimensions. By combining homogenization with quasistatic plasmonic eigenmode expansion, the metasurface is replaced by a zero-thickness interface associated with frequency-dependent effective susceptibilities. The resulting discontinuities of the fields are responsible for strong interaction with the incoming light at the resonances when the complex permittivity of the metal passes close to the real permittivity of an eigenmode. Our modeling provides a physical picture of resonances in plasmonic metasurfaces, and it allows for a huge decrease in the numerical cost of their computations. In addition, comparisons with direct numerics in two dimensions evidence its predictive force at any incidence, particle shape, and arrangement.