This paper presents the development of a three dimensional computational model, based on the Discrete Element Method (DEM), which was used to investigate the effect of the angle of skew on the load carrying capacity of twenty-eight different in geometry single span stone masonry arches. Each stone of the arch was represented as a distinct block. Mortar joints were modelled as zero thickness interfaces which can open and close depending on the magnitude and direction of the stresses applied to them. The variables investigated were the arch span, the span: rise ratio and the skew angle. At each arch, a full width vertical line load was applied incrementally to the extrados at quarter span until collapse. At each load increment, the crack deve...