This study presents a detailed failure analysis of seven masonry churches, which were severely damaged during the 2012 Emilia-Romagna (Italy) seismic sequence, by means of 3D FE non-linear dynamic simulations. The main aims of the numerical investigations are: (1) to carry out a sufficiently wide sensitivity study on different specific case studies to have an insight into the role played by the geometry – which is always unique for churches – and by the irregularities; (2) to validate or address the limits of applicability of the most widespread conventional static approaches recommended by Italian Code, i.e. kinematic limit and pushover analyses. Non-linear dynamic analyses are carried out assuming that masonry behaves as a non-lin...