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

One of the characteristic features of the city of Utrecht is its extensive system of canals and wharf cellars, whose constructions date back as early as the 1200s, and which are now considered as one of the historical properties of the city. A typical wharf cellar in Utrecht comprises a masonry barrel vault with multi-layered rings for the cellar interior, masonry piers which are interconnected to the other wharf cellars, and spandrel walls for the façades. Due to increased traffic volume and urbanization which caused the increase of dead load and traffic load, it is important to assess the structural safety and state of maintenance of these historical structures. In this paper, a novel safety assessment framework for these structures is presented and applied to the analysis of a typical masonry wharf cellar in central Utrecht. The geometry of the cellar is first parametrically generated, which is then used to create a block-based numerical model for analysis using the Distinct Element Method (DEM), where bricks units are modelled as discrete blocks separated by zero thickness interfaces. Traffic loads in accordance with the Dutch Standard traffic model for regular vehicles and emergency service vehicles are calculated and the dispersion through the filling soil is modelled. The ultimate load due to these load configurations is then assessed. The analysis results can be used to identify the critical load cases and the failure mechanisms of the wharf cellar, while also providing general insights into the safety and stability of the cellars, thus aiding engineers in their efforts to extend the lifespan of these historical structures.@en

The city of Utrecht is famously known for the system of canals and the wharf cellars integrated to the heart of the city, whose construction dates back to the 1300s. Due to increased traffic volume which caused the increase in dead load and traffic load, it is important to assess the safety and state of maintenance of these historical structures. In this paper, a safety assessment framework for wharf cellars is introduced and the application to a wharf cellar as a case study in central Utrecht is provided. The geometry of the wharf cellar is parametrically generated and used for the numerical analysis using the distinct element method (DEM), where arch units and piers are modeled as discrete blocks separated by zero-thickness interfaces. Traffic load models in accordance with the Dutch guideline for emergency vehicles are calculated. Unlike traditional approaches, the three-dimensional load distribution through the soil is modeled. The structure’s compliance with this load is assessed, and the failure load and mechanism are observed. The analysis result can be used to help engineers on providing insights into the safety and stability of the cellars in an effort to extend the lifespan of the historical 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