The present paper concentrates on the numerical modelling of masonry vaults, adopting a type of analysis first developed at the University of Parma for applied mechanics, based on the use of non-smooth dynamics software, through a Differential Variational Inequalities (DVI) formulation specifically developed for the 3D discrete elements method. It allows to follow large displacements and the opening and closure of cracks in dynamic field, typical of the masonry vaulted structures. Once the modelling instrument was calibrated, thanks to the comparison with the recurrent damage mechanisms previously analysed, it was also applied to foresee the behavior of the same structure with different actions and with different types of strengtheni...