Numerical modelling of experimental uniaxial and compression-shear tests on traditional stone masonry specimens, using 2d particle models

Stone masonry walls are the main structural elements of many historic buildings. Their restoration and preservation are a major concern given the increasing interest in the rehabilitation of built historical heritage and the implementation of preventive measures to mitigate seismic risk. The accurate structural assessment of the existing rubble stone masonry is a very complex and difficult task, due its composite and complex nature. The heterogeneity and uncertainty in material properties of its constituents, mortar and stone, the variability of the stone units positioning and geometry, among others, make its experimental characterization and accurate numerical modelling still nowadays a challenging task. In this context, the research presented aims to contribute to a better understanding of the in-plane shear behaviour of unreinforced two leaf rubble stone masonry walls, typical stone masonry of ancient buildings that are representative in Portugal. Rubble stone masonry specimens, built with traditional Portuguese construction techniques, tested experimentally under monotonic compression-shear loading conditions are here numerically evaluated, adopting a micro-modelling approach using a 2D particle model (2D-PM). In the 2D-PM model the stone and mortar elements are represented as particle assemblies that interact with each other, thus capable of representing their inherent physical and material heterogeneity. The numerical model is generated through a mapping process of the stone units and mortar joints. The experimental campaign conducted on rubble stone masonry wall specimens allowed to collect the data necessary for the calibration of the PM model parameters, namely the stone-stone and mortarmortar elastic and strength contact properties. The validation of the 2D-PM models of the rubble stone masonry specimens under compression and combined compression-shear loading conditions is performed using the Parmac2D software. The presented results show that 2D-PM models can predict the crack propagation, the final failure modes, the maximum shear strength, and the wall ductility observed experimentally. Parametric studies that allow a better agreement between the numerical predictions and the experimental response, showing the relevance of the stone-mortar interface strength properties in the overall macroscopic behaviour, are also presented.