SUMMARYDue to its ability to take into account discontinuities, the discontinuous Galerkin (DG) method presents some advantages for modeling crack initiations and propagations. This concept has been recently applied to 3D simulations and to elastic thin bodies. In this last case, the assumption of small elastic deformations before crack initiations or propagations reduces drastically the applicability of the framework to a reduced number of materials. To remove this limitation, a full-DG formulation of non-linear Kirchhoff-Love shells is presented and is used in combination with an elasto-plastic finite deformations model. The results obtained by this new formulation are in agreement with other continuum elasto-plastic shell formulations. Then this full-DG formulation of Kirchhoff-Love shells is coupled with the cohesive zone model to perform thin body fracture simulations. As this method allows considering elasto-plastic constitutive laws in combination with the cohesive model, accurate results compared to the experiments are found. In particular, the crack path and propagation rate of a blasted cylinder are shown to match experimental results. One of the main advantages of this framework is its ability to run in parallel with a high speed-up factor, allowing the simulation of ultra fine meshes.