A micromechanically motivated model is proposed to capture nonlinear effects and switching phenomena present in ferroelectric polycrystalline materials. The changing remnant state of the ferroelectric crystal is accounted for by means of so-called back fields—such as back stresses—to resist or assist further switching processes in the crystal depending on the local loading history. To model intergranular effects present in ferroelectric polycrystals, the computational model elaborated is embedded into a mixed polygonal finite element approach, whereby an individual ferroelectric grain is represented by one single irregular polygonal finite element. This computationally efficient coupled simulation framework is shown to reproduce the specific characteristics of the responses of ferroelectric polycrystals under complex electromechanical loading conditions in good agreement with experimental observations.