This study presents a robust contact constitutive model in the distinct element method (DEM) framework for simulating the mechanical behavior of masonry structures. The model is developed within the block-based modeling strategy, where the masonry unit is modeled as deformable blocks with potential crack surfaces in the middle of the bricks, while the mortar joints are defined as zero-thickness interfaces. The modeling strategy implements multi-surface plasticity with damage mechanics, including a tension cut-off, Coulomb failure criterion, and an elliptical compressive cap for the damage in tension, shear, and compression, respectively. Two new features are introduced in this contact model: a piecewise linear softening function for strength degradation in tension and shear and a hardening/softening function to phenomenologically define the compressive damage of masonry composite into the unit-mortar interface. The constitutive model is implemented in commercial DEM software using the small displacement configuration and validated against material and experimental tests on masonry walls subjected to constant pre-compression and monotonically increasing in-plane load. The experimental and numerical results regarding the force-displacement relationship and damage pattern produced by the proposed constitutive model are compared and critically discussed, demonstrating the capability of DEM coupled with the suitable constitutive law in simulating the behavior of masonry structures.

@en

The assessment of the seismic performance of unreinforced masonry cross-vaults is still a challenge in numerical analysis, due to complex curved geometries and bond patterns, and uncertainties related to the selection of adequate modeling strategies, including but not limited to that of material properties, damping scheme, and unit/joint idealization. This paper presents the results of a collaborative effort to validate, against the shake table test of both unstrengthened and strengthened masonry cross-vault specimens as part of the SERA Project Blind Prediction and Post-diction Competition, various discontinuum-based numerical approaches. First, the geometry of the cross-vault is created using a Python-based computational framework to accurately represent the brick arrangement and the shape of the vault. Then, the geometry is converted into an assemblage of deformable blocks and analyzed using the Distinct Element Method (DEM). An elasto-softening contact model based on fracture energy is implemented in the masonry joints to simulate crushing, tensile, and shear failures. The performance of the proposed strategy, conceived for the unstrengthened configuration of the tested vault specimen and then adapted to include the presence of cementitious repairs, shows satisfactory agreement with both qualitative and quantitative experimental responses, also revealing critical insights and lessons learned through the blind/post-prediction exercise.@en