Discrete-continuum approach to assess 3D failure modes of masonry arch bridges

There are two main objectives of this research. First, a full masonry arch bridge with all structural components are considered. The failure mechanism of spandrel wall and backfill-masonry interaction are successfully simulated using a 3D discrete-continuum model as validated by previously published experimental data. Moreover, the influence of the frictional resistance between soil and masonry components is discussed. Second, two different skew arches, with different bond patterns, are analysed to understand the influence of construction method (helicoidal and false) on the damage pattern and capacity. The results of the analysis demonstrated that discrete and mixed discrete-continuum approaches can predict complex 3D collapse mechanisms of masonry arch bridges and provides detailed information about their damage progression.(undefined

Discrete element modeling of masonry structures: Validation and application

The failure mechanism and maximum collapse load of masonry structures may change significantly under static and dynamic excitations depending on their internal arrangement and material properties. Hence, it is important to understand correctly the nonlinear behavior of masonry structures in order to adequately assess their safety and propose efficient strengthening measures, especially for historical constructions. The discrete element method (DEM) can play an important role in these studies. This paper discusses possible collapse mechanisms and provides a set of parametric analyses by considering the influence of material properties and cross section morphologies on the out of plane strength of masonry walls. Detailed modeling of masonry structures may affect their mechanical strength and displacement capacity. In particular, the structural behavior of stacked and rubble masonry walls, portal frames, simple combinations of masonry piers and arches, and a real structure is discussed using DEM. It is further demonstrated that this structural analysis tool allows obtaining excellent results in the description of the nonlinear behavior of masonry structures

Static and impact response of a single-span stone masonry arch

Unreinforced masonry structures are susceptible to man-made hazards such as impact and blast loading. However, the literature on this subject mainly focuses on masonry wall behavior, and there is a knowledge gap about the behavior of masonry arches under high-strain loading. In this context, this research aims to investigate both quasistatic and impact response of a dry-joint stone masonry arch using the discrete element method. Rigid blocks with noncohesive joint models are adopted to simulate dry-joint assemblages. First, the employed modeling strategy is validated utilizing the available experimental findings, and then sensitivity analyses are performed for both static and impact loading, considering the effect of joint friction angle, contact stiffness, and damping parameters. The outcomes of this research strengthen the existing knowledge in the literature regarding the computational modeling of masonry structures that are subjected to usual and extreme loading conditions. The results highlight that applied discontinuum-based numerical models are more sensitive to stiffness parameters in high-strain loading than static analysis

Simulation of uniaxial tensile behavior of quasi-brittle materials using softening contact models in DEM

This study proposes new contact models to be incorporated into discrete element method (DEM) to more accurately simulate the tensile softening in quasi-brittle materials, such as plain concrete and masonry with emphasis on fracture mechanism and post-peak response. For this purpose, a plain concrete specimen (double notched) and stack bonded masonry prism under direct tensile test are modeled. Furthermore, mixed mode crack propagation is investigated in concrete and brickwork assemblages. Two modeling approaches are proposed, the simplified and detailed meso modeling, both based on DEM. In the simplified meso-model, a smooth contact surface is considered between two separate blocks, whereas the internal structure of the material is explicitly represented as a tessellation into random polyhedral blocks in the detailed meso-model. Furthermore, recently developed tensile softening contact constitutive models implemented into a commercial discrete element code (3DEC) are used to simulate the softening behavior of concrete and masonry. As an important novel contribution, it is indicated that the proposed computational models successfully capture the complete (pre- and post-peak) material behavior and realistically replicate the cracking mechanism. Additionally, a sensitivity analysis demonstrates the influence of the various micro-contact parameters on the overall response of the examined materials.Authors would like to express their gratitude to Itasca Educational Partnership Program (IEP) for their kind support and providing 3DEC softwar

In-plane static response of dry-joint masonry arch-pier structures

The majority of historical masonry structures include arches and vaults, constructed with or without (dry-joint) any mortar. This paper focuses on dry-joint masonry, because it is common all around the world among architectural heritage. Furthermore, even if there was a mortar in the original construction, it typically suffers from deterioration over its lifetime, often causing total loss of mortar in many of the joints. Due to large horizontal thrust that can be produced, depending on their geometry, arches are typically supported by heavy buttresses. These structures tend to be difficult to model due to their nonlinear nature and inherent discontinuity, which makes it challenging to evaluate their stability. In that context, it is necessary to have realistic numerical models to better diagnose their structural behaviour in a seismic event and, ultimately, to perform only necessary and beneficial interventions. The main goal of this paper is to assess the seismic performance of various dry-joint arch forms with different masonry pier types (i.e. monolithic and regularly coursed) subjected to incrementally increasing lateral loads proportional to the mass (pushover). To achieve this goal, a parametric study is performed on arch curvature and pier morphology. Moreover, the influence of steel tie-rod reinforcement is also examined on the proposed masonry models. These complex masonry arch systems can be simulated with discrete element modeling (DEM) approach. In this research, a commercial three-dimensional discrete element code, 3DEC, is used; in which masonry units are modeled as distinct blocks with zero tensile strength at their joints. The results reveal that pointed arches provide better seismic resistance than the circular arch form. Furthermore, implemented steel tie-rods yield significant increase in stability for the arch-pier structures, which is quantified on different arch curvatures.- (undefined

Lateral capacity of URM walls: a parametric study using macro and micro limit analysis predictions

This research investigates the texture influence of masonry walls’ lateral capacity by comparing analytical predictions performed via macro and micro limit analysis. In particular, the effect of regular and quasi-periodic bond types, namely Running, Flemish, and English, is investigated. A full factorial dataset involving 81 combinations is generated by varying geometrical (panel and block aspect ratio, bond type) and mechanical (friction coefficient) parameters. Analysis of variance (ANOVA) approach is used to investigate one-way and two-way factor interactions for each parameter in order to assess how it affects the horizontal load multiplier. Macro and micro limit analysis predictions are compared, and the differences in terms of mass-proportional horizontal load multiplier and failure mechanism are critically discussed. Macro and micro limit analysis provide close results, demonstrating the reliability of such approaches. Furthermore, results underline how the panel and block aspect ratio had the most significant effect on both the mean values and scatter of results, while no significant effect could be attributed to the bond types.This work was partly financed by FCT/MCTES through national funds (PIDDAC) under the R&D Unit ISISE under reference UIDB/04029/2020. This study has been partly funded by the STAND4HERITAGE project that has received funding from the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation program (Grant agreement No. 833123), as an Advanced Grant

Analysis and prediction of masonry wallette strength under combined compression-bending via stochastic computational modeling

The out-of-plane flexural bending capacity of masonry is a fundamental property for understanding the behavior of masonry structures. This study investigates the behavior of unreinforced masonry wallettes subjected to combined compression-flexural loading using the discrete element method (DEM), and provides a novel framework to estimate the masonry strength. A simplified micro-modeling strategy is utilized to analyze a masonry wallette, including the variation of the mechanical properties in masonry units and joints. Stochastic DEM analyses are performed to simulate brickwork assemblages, assuming a strong unit-weak joint material model typical of most masonry buildings, including historical ones. Once the proposed computational approach is validated against the experimental findings, the effect of spatial and non-spatial variation of mechanical prop-erties is explored. Two failure types are identified: joint failure and brick failure. For each failure mechanism, the variability of the response and the effects of the modeling parameters on the load-carrying capacity is quantified. Afterward, Lasso regression is employed to determine predictive equations in terms of the material properties and vertical pressure on the wallette. The results show that the most important parameters changing the response are the joint tensile strength and the amount of vertical stress for joint failure, whereas the unit tensile strength dominates the response for brick failure. Overall, this research proposes a novel framework adopting validated advanced computational models that feed on simple test results to generate data that is further utilized for training response prediction models for complex structures.- (undefined

Comparison of in-plane and out-of-plane failure modes of masonry arch bridges using discontinuum analysis

This research aims to provide a better understanding of the structural behavior of masonry arch bridges using advanced modeling strategies. Two main contributions are achieved in this article; first, triggering mechanisms for the out of plane failure of spandrel walls are established; second, the influence of soil backfill on the behavior and strength of the bridges is presented through a comprehensive parametric study. Here, masonry arch bridges are modeled using a discontinuum approach, composed of discrete blocks, including also a continuum mesh to replicate infill material, adopting a framework of discrete element modeling. The equations of motion for each block are solved by an explicit finite-difference method, using the commercial software 3DEC. The results of the preliminary analyses are compared with analytical solutions and limit state analysis for validation purposes. Different arch bridge models, representing common geometrical properties in the northwest Iberian Peninsula are analyzed. Transverse effects, damage patterns and collapse mechanisms are discussed under different types of loading. The analysis demonstrated the severe capacity reduction due to spandrel wall failures and the importance of soil backfill in results, only possible by taking advantage of the performed numerical modeling strategy

Tensile fracture mechanism of masonry wallettes parallel to bed joints: A stochastic discontinuum analysis

Nonhomogeneous material characteristics of masonry lead to complex fracture mechanisms, which require substantial analysis regarding the influence of masonry constituents. In this context, this study presents a discontinuum modeling strategy, based on the discrete element method, developed to investigate the tensile fracture mechanism of masonry wallettes parallel to the bed joints considering the inherent variation in the material properties. The applied numerical approach utilizes polyhedral blocks to represent masonry and integrate the equations of motion explicitly to compute nodal velocities for each block in the system. The mechanical interaction between the adjacent blocks is computed at the active contact points, where the contact stresses are calculated and updated based on the implemented contact constitutive models. In this research, different fracture mechanisms of masonry wallettes under tension are explored developing at the unit–mortar interface and/or within the units. The contact properties are determined based on certain statistical variations. Emphasis is given to the influence of the material properties on the fracture mechanism and capacity of the masonry assemblages. The results of the analysis reveal and quantify the importance of the contact properties for unit and unit–mortar interfaces (e.g., tensile strength, cohesion, and friction coefficient) in terms of capacity and corresponding fracture mechanism for masonry wallettes.This research received no external funding

Macro vs Micro Limit Analysis models for the seismic assessment of in-plane masonry walls made with quasi-periodic bond types

Masonry bond patterns can considerably affect the seismic performance of in-plane walls. Although several numerical and experimental works addressed this topic, few attempts tried to investigate such an issue using analytical formulations. This paper aims to compare macro and micro limit analysis models investigating masonry walls arranged with different bond types, namely Running, Flemish and English. A dataset involving 81 combinations is generated by varying geometrical (panel aspect ratio, block aspect ratio, bond type) and mechanical (friction coefficient) parameters. Finally, one-way and two-way factor interactions are used to evaluate how each parameter affects the horizontal load multiplier and assess matching among the two adopted formulations.This work was partly financed by FCT/MCTES through national funds (PIDDAC) under the R&D Unit ISISE under reference UIDB/04029/2020. This study has been partly funded by the STAND4HERITAGE project that has received funding from the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation programme (Grant agreement No. 833123), as an Advanced Gran