In this study, a comprehensive 3D non-linear finite element modelling procedure is presented for the simulation of reinforced concrete flat plates, retrofitted with externally bonded carbon fibre-reinforced polymer (CFRP) strips and subjected to punching shear loads. The presented modelling procedure may be used to conduct parametric numerical analyses on the response of reinforced concrete elements similar to considered in this study. The modelling procedure is applied on: (i) a solid, un-retrofitted flat plate and (ii) eight flat plates with openings and retrofitted with externally bonded CFRP strips. In the procedure, debonding of CFRP strips from concrete surfaces is accounted for by defining appropriate bilinear shear stress-slip and normal stress-gap relationships. In the model, the corner lift-offs observed at the corners of the flat plates are also considered by defining appropriate restrains and frictional support conditions. Results obtained from the proposed modelling procedure are verified using experimental data available in the literature. From the verification, accuracy of the proposed modelling procedure is confirmed in terms of structural responses such as load displacement relationships and crack propagations of the test specimens.