Cyclic constitutive model for masonry joint damage and energy dissipation using the distinct element method

This paper presents a cyclic joint constitutive model within a Distinct Element Method framework to simulate the in-plane response of unreinforced masonry structures. The model combines multi-surface failure criteria, including tensile cut-off, Coulomb friction, and an elliptical compression cap. It incorporates exponential softening, a unified damage scalar for stiffness degradation, and a hardening–softening law for compression. Shear-induced dilatancy is captured via an uplift-correction mechanism with an exponential dilatancy-decay law, while stiffness degradation governs energy dissipation. The model is validated at both material and structural scales. Material-level simulations of cyclic compression and shear tests show close agreement with experimental data. Structural-scale validation on full-height calcium-silicate walls under combined compression and cyclic lateral loading demonstrates the ability to reproduce rocking-dominated, shear-dominated, and hybrid failure mechanisms. The model successfully replicated global hysteretic force–drift loops, capturing stiffness decay and energy dissipation, as well as local failures like cracking, sliding, and toe crushing. The model also reproduced the drift-dependent transition from rocking to friction-controlled sliding, a key mechanism for earthquake assessment. By integrating these features into a single, efficient framework, the proposed constitutive model provides a robust tool for evaluating seismic performance and conserving heritage.