Simulating damage phenomena in masonry structures remains a significant challenge because of the intricate and heterogeneous nature of this material. An accurate evaluation of fracture behavior is essential for assessing the bearing capacity of these structures, thereby mitigating dramatic failures. This paper proposes an innovative adaptive concurrent multiscale model for evaluating the bearing capacity of in-plane masonry structures under in-plane loadings. Developed within a Finite Element (FE) set, the proposed model employs a domain decomposition scheme to solve a combination of fine- and coarse-scale sub-models concurrently. In regions requiring less detail, the masonry is represented by homogeneous linear elastic macro-elements. The material properties for these macro-elements are derived through a first-order computational homogenization strategy. Conversely, in areas with higher resolution needs, the masonry is modeled by accurately depicting individual brick units and mortar joints. To capture strain localization effectively in these finer regions, a Phase Field Cohesive Zone Model (PF-CZM) formulation is employed as the fracture model. The adaptive nature derives from the fact that at the beginning of the analysis, the model is entirely composed of coarse regions. As nonlinear phenomena develop, these regions are progressively deactivated and replaced by finer regions. An activation criterion identifies damage-prone regions of the domain, thereby triggering the transition from macro to micro scales. The proposed model’s validity was assessed through multiscale numerical simulations applied to a targeted case study, with the results compared to those from a direct numerical simulation. The results confirm the effectiveness and accuracy of this innovative approach for analyzing masonry failure.
Greco et al. (2025) studied this question.