Ferroelectric materials exhibit a huge potential for engineering applications -ranging from electrical actuators (inverse piezoelectric effect) to sensor technology (direct piezoelectric effect). To give an example, lead zirconate titanate (PZT) is a typical perovskite ion crystal possessing ferroelectric properties. In this contribution, we are particularly interested in the modelling of microcracking effects in ferroelectric materials.In view of Finite-Element-based simulations, the geometry of a natural grain structure, as observed on the so-called micro-level, is represented by an appropriate mesh. While the response on the grains themselves is approximated by coupled continuum elements, grain boundaries are numerically incorporated via so-called cohesive-type elements. For the sake of simplicity, switching effects in the bulk material will be neglected. The behaviour of the grain boundaries is modelled by means of cohesive-type laws. Identifying grain boundaries as potential failure zones leading to microcracking, cohesive-type elements consequently offer a great potential for numerical simulations. As an advantage, in the case of failure they do not a priori result in ill-conditioned systems of equations as compared with the application of standard continuum elements to localised deformations. Finally, representative constitutive relations for both the bulk material and the grain boundaries, enable two-dimensional studies of low-cycle-fatigue motivated benchmark boundary value problems.