This study presents a finite element model (FEM) to reproduce the punching shear failure of fibre-reinforced concrete (FRC) slabs. The mechanical response of the material is characterised through a cohesive crack model that uses a trilinear softening law, which allows for an accurate representation of the fibre bridging effect during the fracture process. The cohesive model using a trilinear diagram as softening function has been successfully applied in the past to symmetric three-point bending tests [1], involving Mode I failure, and to non-symmetric three-point bending tests [2], where Mode I-dominant mixed-mode fracture takes place. In this study, the model is validated for punching shear, where Mode II becomes predominant. The experimental work by Su et al. [3] is used as a benchmark, where circular flat slabs of an Engineered Cementitious Composite (ECC) are tested using different fiber proportions and span-depth ratios of the punching shear test. The FEM models are developed with OOFEM, a free finite element code, and use a smeared crack formulation for the material model [4]. The results show that this simulation strategy provides good agreement in terms of the load-deflection curve and also identifies the principal cracking patterns observed experimentally.