Research Article  Open Access
Daniele Baraldi, Claudia Brito De Carvalho Bello, Antonella Cecchi, "Refined Rigid Block Model for InPlane Loaded Masonry", Advances in Civil Engineering, vol. 2020, Article ID 8844759, 13 pages, 2020. https://doi.org/10.1155/2020/8844759
Refined Rigid Block Model for InPlane Loaded Masonry
Abstract
In this work, a refined rigid block model is proposed for studying the inplane behavior of regular masonry. The rigid block model is based on an existing discrete/rigid model with rigid blocks and elastoplastic interfaces that already proven its effectiveness in representing masonry behavior in linear and nonlinear fields. In this case, the proposed model is improved by assuming rigid quadrilateral elements connected by onedimensional nonlinear interfaces, which are adopted both to represent mortar (or dry) joints between the blocks and also to represent inner potential cracks into the blocks. Furthermore, the softening behavior of interfaces in tension and shear is taken into account. Several numerical tests are performed by considering masonry panels with regular texture subjected to compression and shear. Particular attention is given to the collapse mechanisms and the pushover curves obtained numerically and compared with existing numerical and laboratory results. Furthermore, the numerical tests aim to evaluate the applicability limits of the proposed model with respect to existing results.
1. Introduction
Masonry is a heterogeneous structural material made of natural or artificial blocks connected by dry or mortar joints. Due to the large amount of masonry historical buildings that may be found in Italy and in other south European countries and considering the frequent seismic events in such areas, the assessment of the behavior of masonry and the development of numerical models for masonry structural elements or buildings is an active field of research in Architecture and Civil Engineering.
Among the different numerical strategies that are frequently adopted for studying masonry and that may be found in literature [1–3], this work mainly focuses on a rigid block model, together with a finite element (FE) model based on a total rotating strain crack approach (FETRSCM) assumed as a comparison for some numerical tests. The first model may be considered a micro or local model, while the latter may be considered a macro or diffused model. A comparative analysis between the numerical models and a comparison with experimental results allow to perform a multiscale analysis of masonry behavior and to validate the new proposed model.
Rigid block models are based on several hypotheses that are typical of the wide field of discrete approaches. In this field, the discrete element (DE) model or method is the most common tool adopted in solid and structural mechanics. In general, DE models are characterized by the definition of elements and contacts between the elements. Following the original definition of the method [4–6], a DE model focuses on specimens composed by independent elements, which may be subjected to large displacements, and with a continuously updated determination of the contacts between the elements during a numerical test. DE approaches may be successfully adopted for studying masonry, thanks to its heterogeneity, and several models based on the DE method, originally introduced for soil and rock mechanics, were extended to masonry mechanics [7].
In literature, numerous discrete models based on elastic [8–11] or shear deformable [12] blocks, connected by zerothickness nonlinear interfaces are present, also considering cyclic loads [11]. Furthermore, in the field of discrete modeling for masonry, comparisons between such models and the continuous ones are frequently done [13, 14], also by performing multiscale approaches [15].
The rigid block model adopted here is particularly suitable for studying historical masonry, since it considers blocks as rigid bodies and dry or mortar joints as nonlinear interfaces. It is based on the original numerical model introduced by Cecchi and Sab [16], which follows small displacements hypothesis. Hence, it considers independent elements in contact, but it does not update the contact topology during numerical tests, given that small displacements are assumed. This contribution was extended to the field of material nonlinearity adopting a tension strength criterion and a Mohr–Coulomb yield criterion for restraining interface actions [17]. Such updated model was still based on small displacement hypothesis, which can be considered reasonable for masonry even in case of nonlinear behavior. Furthermore, no contact detection algorithms were proposed, allowing to obtain fast numerical solutions by means of a stiffness matrix approach, which turned out to be an advantage with respect to more general DE models extended to masonry. The resulting nonlinear rigid model turned out to be similar to the rigid block models adopted for performing limit analysis of masonry [18–22]. However, the original nonlinear rigid block model was able to perform pushover tests in load control, but it was not able to reproduce the softening behavior typical of masonry, given that it considers an elasticperfectly plastic interface behavior, which does not respect the actual response in tension and shear typical of mortar joints. However, this preliminary nonlinear rigid block model was calibrated for modeling inplane loaded masonry panels with dry joints [23–25], showing a good agreement between more refined DE models, namely the combined finite elementdiscrete element (FEDE) model [26], and laboratory results.
In this contribution, the new aspect introduced for improving the original rigid block model regards the adoption of a softening law in tension and shear for representing the postelastic behavior of mortar joints and the potential cracking of the blocks. In this latter case, a new inner block joint type is introduced, by subdividing each block into two half blocks, in order to simulate the potential cracking of the resisting elements due to tension and shear. This last aspect was already considered by other authors for better modeling masonry behavior in case of elastic blocks both in and outofplane loaded [8, 10, 11] and also in case of rigid blocks [27] by performing limit analyses.
Focusing then on FE models for masonry, among the different models that may be found in literature, two main approaches may be considered: micro and macromodeling [28]. The first approach envisages the use of standard heterogeneous FE models, with joints and blocks modeled by different materials and specific FEs, or with joints often modeled by means of zerothickness interface elements [8, 29]. The second approach envisages the definition of a continuous material equivalent to the heterogeneous one, modeled by a unique type of FE. In this case, homogenization approaches may be adopted by identifying a periodic cell typical of the masonry material and by adopting local linear and nonlinear mechanical parameters [30–33], also by performing multiscale analyses [34]. In particular, multiscale approaches are focused on Cauchy or micropolar/nonlocal continuum models [35–39]. Another possibility in the field of continuous modeling is an approach based on the evaluation of the overall properties of a masonry assemblage by fitting experimental tests. Orthotropic and also isotropic damage models are often adopted for this purpose [40]. Focusing on damage models, discrete or smeared crack approaches may be adopted by extending existing models originally introduced for modeling concrete material. The continuous total strain crack model (TSCM) is a proper tool for modeling strengthened masonry structures [41–43] and homogenous isotropic materials as rammed earth masonry [44, 45]. Moreover, it has also been adopted for modeling unreinforced masonry panels under inplane loads [40, 46, 47]. An interesting comparison between several numerical finite element models for the assessment of the behavior of unreinforced masonry has been performed by Bartoli et al. [48], also considering the performance of the total strain crack model.
This work is dedicated to the assessment of the nonlinear behavior of masonry panels with regular texture subjected to vertical and horizontal inplane actions. The aim of this work is to validate the proposed improvements of an existing rigid block model by performing several numerical tests and comparing the results of the updated model with respect to existing laboratory results and other accurate numerical results. The calibration already performed by the authors for rigid, FE/DE, and FETRSCM models [49] is extended to the case of the updated rigid block model. Before performing tests on masonry walls, the nonlinear behavior of mortar joints is evaluated for first by performing some numerical tests on two blocks specimens. Then, the experimental campaign by Raijmakers and Vermeltfoort [50] on almost squareshaped masonry panels is taken as reference, together with the work of Lourenço and Rots [8], that studied the same panels by means of FE models with interface elements. Such case studies have been recently considered by the authors in a comparison between continuous and discrete models focusing on nonlinear behavior [51] and by other authors in a similar comparison by means of a multiscale/multimodel approach [15]. Here, the pioneering work of Page [29] is also taken into consideration as a further case study.
The paper is organized as follows: the numerical models adopted are introduced, focusing on the updated rigid block model; then, the case studies considered are described, and the results of the numerical tests performed with the adopted models are showed and discussed, accounting also for laboratory results. Some comments and possible future developments are presented in the final part of the work.
2. Numerical Models
2.1. Refined Rigid Block Model–Geometry
A oneleaf masonry panel with regular texture and following the socalled “running bond” pattern is considered (Figure 1(a)). Block dimensions are width b, height a, and thickness s, whereas mortar joint thickness is e, which is assumed to be equal for both head (vertical) and bed (horizontal) joints. A refined rigid block model is defined by introducing rigid elements coincident with one half of each actual resisting element, hence having width b/2 and height a (Figure 1(b)). The resulting discretization is characterized by half blocks aligned both horizontally and vertically and in contact with four neighboring half blocks by means of four interfaces coincident with the four half block edges. Horizontal interfaces are always represented by mortar interfaces, whereas vertical interfaces can be, alternatively, mortar interfaces or inner block interfaces; the latter interface case may represent a potential vertical crack for the original block, which can be obtained in case of a compressive action on the masonry structure without or with a contemporary shear action.
(a)
(b)
Assuming a twodimensional coordinate system Oy_{1}y_{2}, a plane stress state, and considering the hypothesis of rigid elements and small displacements, the inplane displacement of each half block B_{i,j} is a rigid body motion defined by half block center inplane translations u_{ij} = {u_{1}u_{2}}^{T} and half block inplane rotation ω_{3,ij} with respect to its center:where Ω_{ij} = skw (ω_{3,ij}). Nodal degrees of freedom of a generic element can be collected in a vector .
The interactions between the half blocks through the interfaces are represented by the resultants of the corresponding contact stresses, which are given by a normal force F_{n}, a shear force F_{s} and a bending moment M: f_{loc} = {F_{n} F_{s} M}^{T}, which are defined considering the local orientation of the interface. The corresponding interface actions in the global coordinate system are f = {F_{1} F_{2} M}^{T} = R f_{loc}, where R is the rotation matrix of the interface depending on its horizontal or vertical orientation. Such interactions depend on the relative displacements and rotations between adjacent half blocks, defined as d = {d_{1} d_{2} δ}^{T}, which can be written in the form of local relative displacements d_{loc} = {d_{n} d_{s} δ}^{T} depending on interface orientation, d = {d_{1} d_{2} δ}^{T} = R d_{loc}. These relative displacements can also be written as a function of half block global displacements in vector form by introducing a compatibility matrix H detailed for each block and collecting the distances between each half block center and interface center:withwhere k_{1} = 1, k_{2} = 0 for a horizontal interface, and k_{1} = 0, k_{2} = 1 for a vertical one. It is worth noting that the compatibility matrix H assembled over the entire masonry specimen considered turns out to be coincident with the matrix commonly adopted for performing limit analysis of masonry [20].
2.2. Refined Rigid Block Model–Linear and Nonlinear Behavior
In the elastic field, the constitutive relation that defines the interface actions between adjacent elements is f_{loc} = K_{loc}d_{loc}, where K_{loc} is the local stiffness matrix of the interface collecting the stiffness parameters depending on interface type and geometric dimensions (area A and moment of inertia I). For instance, K_{loc} = {k_{n} k_{s} k_{r}}^{T} = {K_{n}A K_{s}A K_{n}I}^{T}, where K_{n} and K_{s} are, respectively, interface normal and shear stiffness, that in case of mortar joints are defined as the function of mortar elastic modulus E_{m} and Poisson’s ratio ν_{m}, together with mortar joint thickness e:
With respect to the original rigid block model proposed by authors [16, 17], this work introduces the elastic stiffness parameters of inner block joints, which are determined in the same manner of mortar joints parameters, but they are assumed to be several times larger than those of mortar joints in order to simulate the extremely limited deformability of a block.
It is worth noting that the proposed rigid block model allows to define the elastic stiffness of a masonry assemblage, by assembling the local stiffness matrices over the entire structural element considered. The elastic stiffness matrix of the entire assemblage is not updated during the following numerical tests, namely by adopting a modified Newton–Raphson method based on initial stiffness for studying the nonlinear behavior of masonry specimens, in order to avoid, at each step of the test and for each interface, the evaluation of the secant stiffness or the tangential stiffness, which may lead to a local or global singular stiffness matrix.
The nonlinear behavior of the model is taken into account by adopting a nonlinear part of the constitutive laws after the initial elastic behavior. A tensile strength and a Mohr–Coulomb yield criterion are assumed for defining the elastic limits of tension and shear actions, respectively; whereas a maximum eccentricity criterion is assumed for the bending moment by introducing l as one half of the interface length:
Differently, with respect to the work of Lourenço and Rots [8], in this contribution, an elasticperfectly plastic behavior in compression is assumed for both mortar and inner block interfaces, by introducing a limit in compression F_{c} = −f_{c}A.
These criterions hold for both mortar joints and inner block joints and are based on the corresponding material properties represented by tensile strength f_{t}, compressive strength f_{c}, cohesion c, and friction ratio μ.
Furthermore, in this work, the softening behavior in case of tensile and shear failure, typical of brittle materials, is assumed by introducing the fracture energy G_{I} for the first (tension) and G_{II} for the second (sliding or shear) modes of cracking.
Starting from the elastic limits of Equation (5), the corresponding relative displacements are given byand the softening behavior is defined as follows if the relative normal and shear displacements are larger than the corresponding elastic limits:
The above expressions are characterized by an exponential coefficient that tends to zero for increasing relative displacement in normal and tangential direction. For simplicity, the softening behavior in case of bending moment is not directly introduced as a function of relative rotation, given that the limit in terms of bending moment already accounts for the updated normal force of the interface.
In the following numerical tests, the elastic and inelastic parameters for mortar joints and inner block joints will be assumed coincident with those adopted by other numerical tests assumed as reference.
2.3. FETSRCM
The FETRSCM is a macromodeling approach originally introduced for concrete structures. The model can be extended to masonry by assuming a heterogeneous material as a continuous one, with bricks and mortar considered together as a single homogenous material through a total strain rotating crack model approach [52]. Recently [49], the FETRSCM has been adopted by authors for modeling inplane loaded masonry and compared with a simple rigid block model and a combined FEDE model. Such a model has been also compared by other authors with other numerical tools for assessing the seismic behavior of unreinforced masonry [48]. The continuous material follows a stressstrain law in compression and tension, where the mechanical properties are assigned in agreement with the information from the experimental test, and they are not determined by means of a proper homogenization procedure. An elastic modulus E and a Poisson’s ratio ν are assumed for the continuous material, and a damaged plasticity model is adopted for representing the material nonlinear behavior, based on the incremental plasticity theory and on the concept of isotropic damage elasticity to describe the irreversible damage that occurs during the fracturing process. In particular, incremental plasticity relates plastic stress increments with plastic strain increments by means of a constitutive matrix that adopts a softening postcrack parameter and a shear retention factor. Two failure mechanisms are assumed: cracking in tension and crushing in compression. An exponential softening law and a parabolic law are assumed, respectively, in tension and compression. Both constitutive laws are defined by limit stresses (f_{t} in tension and f_{c} in compression) and by the corresponding fracture energies (G_{t} and G_{c}), which are considered in terms of dimensionless values with respect to finite element size h [53]: G_{t} = G_{I}/h and G_{c} = G_{C}/h. In particular, G_{C} is obtained by integrating the ellipsoid law for the cap model in compression adopted in the study by Lourenco and Rots [8], and h is defined as a function of the area (or the volume) of the finite element adopted in the numerical model. In the following numerical tests, a commercial computer code [54] is adopted, and quadrilateral finite elements in the plane state are assumed, with a regular mesh of almost squareshaped elements is defined; hence, h is assumed to be coincident to the height or the width of the generic finite element.
3. Numerical Tests
3.1. Introduction to the Tests
Several numerical tests are performed in order to evaluate the effectiveness of the proposed model for correctly simulating the behavior of masonry panels subject to compression and shear.
3.2. Calibration of Mortar Interface Nonlinear Parameters
The nonlinear behavior of the interfaces in the rigid block model is calibrated for first by performing a set of simple numerical tests involving a few blocks connected by mortar interfaces and subjected to simple tensile and shear actions.
For this purpose, the simple tension and shear tests on specimens made of two blocks (Figure 2(a) and 2(b), respectively), having dimensions b = 250 mm, a = 50 mm, and s = 120 mm, similar but not coincident to those of the work assumed as reference [8], are taken into consideration and reproduced by the proposed refined rigid block model. For simplicity, the mechanical parameters adopted for the shear tests, namely c = 0.87 MPa, μ = 0.73, and G_{II} = 0.058–0.13 σ N/mm, are extended to the tension test by setting f_{t} = c/1.4 = 0.62 MPa and G_{I} = 0.012 N/mm. Such mechanical parameters have been determined by Van der Pluijm [55, 56] by means of laboratory experimental tests in tension and shear.
Results in terms of normal stresses given by interface normal force F_{n} over interface area A for increasing interface relative normal displacement d_{n} are shown in Figure 3(a), with a softening branch after reaching the tensile strength set as input of the model. Similarly, results in terms of interface shear stresses given by interface shear force F_{s} over interface area A for increasing interface relative tangential displacement d_{s} are shown in Figure 3(b) for three different levels of normal compression. Such a compression level increases the shear strength of the mortar interface and turns out to be close to the residual shear strength attained at the end of each test, which is equal to the contribution given by friction.
(a)
(b)
3.3. Panels Subjected to Compression and Shear
The case studies considered for the numerical tests are given by the laboratory tests on rectangular masonry panels performed by Raijmakers and Vermeltfoort [50]. The geometric characteristics of the masonry panels are as follows: blocks dimensions are width b = 204 mm, height a = 50 mm, and thickness s = 98 mm. Mortar joint thickness is e = 10 mm. The overall dimensions of the specimen are length L = 990 mm, height H = 1000 mm, and thickness s = 98 mm. These dimensions are obtained with 4.5 blocks in horizontal direction and 18 blocks in vertical direction (Figure 4). A steel beam at panel base is used for fixing block displacements and rotation, and a similar steel beam at panel top is used for avoiding relative vertical block displacements and for applying different levels of vertical compression (p = −0.30; −1.21;−2.12 MPa) and an increasing horizontal force F.
The main mechanical parameters adopted for mortar interfaces and inner block interfaces are listed in Table 1 and were determined by performing laboratory tests on blocks, mortar, and mortar joints. In particular, mortar interface stiffness parameters are obtained starting with a material having E_{m} = 800 MPa, E_{b} = 4672 MPa, and ν_{m} = 0.14, and their nonlinear behavior in shear is characterized by cohesion c = 1.4 f_{t} and friction ratio μ = 0.75, whereas inner block interfaces nonlinear behavior is considered only for tensile and bending failure. The FETRSCM adopts E = 3960 MPa, ν = 0.14, f_{c} = 8.8 MPa, and G_{C} = 2 N/mm.

The numerical results obtained with the proposed refined rigid block model turn out to be in quite good agreement with experimental results. Focusing on the loaddisplacement curves (Figure 5), the initial stiffness of the panels is correctly reproduced by the numerical model in the three different compression cases. The level of initial failure of the panels is correctly determined with p = 1.21 MPa, whereas it is slightly underestimated with p = 0.30 MPa and overestimated with p = 2.12 MPa. The residual strength of the panels is correctly determined only with p = 0.30 MPa, probably thanks to the small level of compression of the case study.
(a)
(b)
(c)
Figure 6 shows the deformed configurations (left column, Figures 6(a)–6(c)) and the interface damage (right column, Figures 6(d)–6(f)) obtained at d = 2 mm during the pushover tests for the three compression cases considered. Deformed configurations are characterized by the opening of vertical mortar joints (cyan color in Figures 6(d)–6(f)) and by the sliding of horizontal joints (red color in Figures 6(d)–6(f)) along the diagonal direction of the panel from panel top left corner to panel bottom right corner. Increasing the level of compression, the interface damage increases and involves also vertical inner block interfaces, which are mainly characterized by bending failure (yellow color in Figures 6(d)–6(f)). The deformed configurations reached at d = 2 mm with p = 1.21 and 2.12 MPa (Figures 6(b) and 6(c)) are similar, but a slightly large vertical deformation can be observed with 2.12 MPa. In these numerical simulations, the compressive strength of mortar joints is never reached; hence, the corresponding representation is not highlighted in the legend of Figure 6.
(a)
(b)
(c)
(d)
(e)
(f)
Figure 7 collects the deformed configurations and the tangential stresses obtained with the FETSRCM at d = 2 mm during the pushover tests where the damage along the diagonal of the panel is more evident in the case with low compression level and less evident in the case of high compression level.
(a)
(b)
(c)
The comparison between the numerical results of the proposed refined rigid block model with respect to laboratory tests and numerical results obtained with FETSRCM shows that the behavior of mortar interfaces subjected to compression has to be taken into account into the rigid model, given that the hypothesis of elastic behavior and unlimited strength in compression assumed in the rigid model does not allow to correctly reproduce the softening branch of the pushover tests, together with the compressive failure of some mortar joints close to the top left and bottom right corner of the panels. However, even if loaddisplacement curves obtained with the FETSRCM are very close to the actual behavior of the panels, the deformed configurations are not able to highlight the local damage of mortar joints or blocks.
A further comparison between the results given by the proposed refined rigid block model with respect to the original one is briefly described in appendix in terms of loaddisplacement curves. For further details on the results given by the standard rigid block model, the recent contribution by authors [49] can be taken into consideration.
3.4. Masonry Deep Beam Subjected to a Vertical Force
In order to further evaluate the proposed refined rigid block model, the “Page test” [29] is simulated. Such test is one of the most used in the calibration of numerical models for masonry panels [1, 37, 57], and it was used by authors for an initial calibration of the original nonlinear discrete/rigid block model [17]. The specimen is made of pressed clay blocks having dimensions b = 122 mm, a = 34 mm, and s = 37 mm, with mortar joint thickness e = 5 mm; overall panel dimensions are L = 757 mm and H = 457 mm. The panel is supported at each side of the base, and it is loaded by a compressive symmetric force P applied by a stiff steel beam long 381 mm. Load and restraint conditions are shown in Figure 8. Brick and mortar linear and nonlinear mechanical parameters are presented in Table 2 where several data are taken from the work of Lourenço [1]. The FETRSCM adopts E = 4125 MPa, ν = 0.15, f_{c} = 13 MPa, and G_{C} = 1 N/mm.

Figure 9 shows the loaddisplacement curve obtained with the rigid block and continuous numerical models. The refined rigid block model reaches an ultimate force close to 122.2 kN that turns out to be in quite good agreement with experimental result equal to 109.2 kN, whereas the FETSRCM is in excellent agreement with such results, since it reaches an ultimate force close to 108.2 kN. In this case, the continuous model shows a larger deformability with respect to the rigid one, and the pushover curves are characterized by the same initial stiffness up to 20 kN; then the FETSRCM starts to exhibit damage in compression, which causes a stiffness reduction of the model.
Figures 10 and 11(a) show, respectively, the deformed configuration of the deep beam and the damage map at P = 100 kN, which correspond to the last load increment considered by Page. The deformed shape of the panel is in agreement with the one proposed by Lourenço [1], and the damage of the panel is characterized by the shear, bending, and tension failure of several vertical joints, in quite good agreement with the numerical results of Page [29] and Lourenço [1]. Figure 11(b) shows the damage map at P = 120 kN, which is characterized by the compression failure of several bed joints (blue color) and the shear failure of two inner block joints (red color, dashed line) close to the upper steel beam, together with the further tension and shear failure of some vertical mortar joints.
(a)
(b)
Figure 12 shows the deformed configuration and the map of normal stresses in y_{2} direction obtained with the FETSRCM at P = 100 kN. The deformed configuration is in quite good agreement with that obtained with the rigid model, even if the continuous model shows a more diffused deformability along the entire specimen. The map of normal stresses clearly shows the load transfer from the upper steel beam and the lower supports.
4. Conclusions
In this work, a refined rigid block model has been proposed for studying the inplane behavior of masonry with regular texture. The refined rigid block model is based on an existing model with rigid blocks and elastoplastic interfaces that already demonstrated its effectiveness in representing masonry behavior in linear and nonlinear fields, and it was compared with other commercial DE models and laboratory results. However, the original elasticperfectly plastic law adopted for restraining interface actions of dry or mortar joints was not sufficiently in agreement with actual material behavior; moreover, the model was not able to represent the possible block cracking. For these reasons, the proposed model improves the original one by assuming rigid quadrilateral elements connected by onedimensional nonlinear interfaces, which are adopted both to represent mortar (or dry) joints between the blocks and also to represent inner potential cracks into the blocks. Furthermore, the softening behavior of interfaces in tension and shear is taken into account. The computational effort required for performing numerical tests with the updated rigid block model is still smaller than that typical of commercial DE and FE models, since the degrees of freedom of the rigid block model are lumped at block centers.
Several numerical tests are performed for first by simulating the tension and shear cracking of two blocks connected by a mortar joints and then by considering masonry panels with regular texture subjected to compression and shear. Particular attention is given to the collapse mechanisms and the pushover curves obtained numerically and compared with existing numerical and laboratory results. The numerical simulation of the three shear tests with different compression levels originally performed by Raijmakers and Vermeltfoort [50] turned out to be in quite good agreement with the laboratory test and other numerical results obtained with a FETRSCM approach. Furthermore, the numerical simulation of the deep beam test originally performed by Page [29] showed a good agreement with laboratory results both in terms of ultimate load and crack pattern.
Further developments of this work will focus on a more accurate simulation of the behavior of joints in compression, by taking into account, for instance, of the softening behavior typical of the compressed mortar joints, in order to better reproduce the postpeak branch of the pushover tests. The proposed rigid block model will be further validated and calibrated by simulating other laboratory tests and more complex case studies, such as masonry walls and façades with openings.
Appendix
The refined rigid block model proposed in this work can be used as the original one by simply avoiding the deformability and potential damage of the inner block interfaces. However, the original rigid block model is characterized by a smaller number of degrees of freedom involved in the numerical tests, namely one half of the overall number of degrees of freedom, with respect to the refined rigid block model. A comparison between the loaddisplacement curves and the damage maps obtained with the refined and standard rigid block models is showed in the following figure by performing the first three pushover tests considered in Section 3.3.
As it has already been highlighted in the introduction and in Section 2.2, the original rigid block model is not able to reproduce the softening behavior typical of mortar joints in tension and shear; hence, the original rigid block model is not able to represent a possible reduction of shear capacity after the damage has occurred on the masonry panel. This aspect is less evident in the shear test with the lower compression level (Figure 13(a)), where also the refined rigid block model is not able to correctly represent the larger shear resistance of the panel and the subsequent softening and residual resistance. The difference between the refined and standard rigid block model is more evident in the other two cases with a larger level of compression on the panels (Figures 13(b) and 13(c)).
(a)
(b)
(c)
Data Availability
The data used to support this study are available from the corresponding authors upon request.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
This work was financed by the Research Project PRIN 2017 (grant 2017HFPKZY_002), project “Modelling of constitutive laws for traditional and innovative building materials.”
References
 P. B. Lourenço, “Computational strategies for masonry structures,” Delft University, Delft, Netherlands, 1996, Ph.D. thesis. View at: Google Scholar
 P. Roca, M. Cervera, G. Gariup, and L. Pelà, “Structural analysis of masonry historical constructions. classical and advanced approaches,” Archives of Computational Methods in Engineering, vol. 17, no. 3, pp. 299–325, 2010. View at: Publisher Site  Google Scholar
 H. Smoljanovic, N. Zivaljic, and Z. Nikolic, “Overview of the methods for the modelling of historical masonry structures,” Gradevinar, vol. 65, no. 7, pp. 603–618, 2013. View at: Google Scholar
 P. A. Cundall, “A computer model for simulating progressive large scale movements in blocky rock systems,” in Proceedings of the Symposium of International Society of Rock Mechanics, Nancy, France, 1971. View at: Google Scholar
 P. A. Cundall and R. D. Hart, “Development of generalized 2d and 3d distinct element programs for modeling jointed rock,” ITASCA Consulting Group, Minneapolis, MN, USA, 1985, Final report. View at: Google Scholar
 P. A. Cundall and R. D. Hart, “Numerical modelling of discontinua,” Engineering Computations, vol. 9, no. 2, pp. 101–113, 1992. View at: Publisher Site  Google Scholar
 J. V. Lemos, “Discrete element modeling of masonry structures,” International Journal of Architectural Heritage, vol. 1, no. 2, pp. 190–213, 2007. View at: Publisher Site  Google Scholar
 P. B. Lourenço and J. G. Rots, “Multisurface interface model for analysis of masonry structures,” Journal of Engineering Mechanics, vol. 123, no. 7, pp. 660–668, 1997. View at: Publisher Site  Google Scholar
 A. M. D’Altri, S. de Miranda, G. Castellazzi, and V. Sarhosis, “A 3D detailed micromodel for the inplane and outofplane numerical analysis of masonry panels,” Computers & Structures, vol. 206, pp. 18–30, 2018. View at: Publisher Site  Google Scholar
 L. Macorini and B. A. Izzuddin, “A nonlinear interface element for 3D mesoscale analysis of brickmasonry structures,” International Journal for Numerical Methods in Engineering, vol. 85, no. 12, pp. 1584–1608, 2011. View at: Publisher Site  Google Scholar
 E. Minga, L. Macorini, and B. A. Izzuddin, “A 3D mesoscale damageplasticity approach for masonry structures under cyclic loading,” Meccanica, vol. 53, no. 7, pp. 1591–1611, 2018. View at: Publisher Site  Google Scholar
 F. Cannizzaro, B. Pantò, S. Caddemi, and I. Caliò, “A Discrete MacroElement Method (DMEM) for the nonlinear structural assessment of masonry arches,” Engineering Structures, vol. 168, pp. 243–256, 2018. View at: Publisher Site  Google Scholar
 M. Pepe, M. Sangirardi, E. Reccia, M. Pingaro, P. Trovalusci, and G. de Felice, “Discrete and continuous approaches for the failure analysis of masonry structures subjected to settlements,” Frontiers in Built Environment, vol. 6, no. 43, pp. 1–14, 2020. View at: Publisher Site  Google Scholar
 M. Pepe, M. Pingaro, P. Trovalusci, E. Reccia, and L. Leonetti, “Micromodels for the inplane failure analysis of masonry walls: limit analysis, FEM and FEM/DEM approaches,” Frattura Ed Integrita Strutturale, vol. 14, no. 51, pp. 504–516, 2020. View at: Publisher Site  Google Scholar
 E. Reccia, L. Leonetti, P. Trovalusci, and A. Cecchi, “A multiscale/multidomain model for the failure analysis of masonry walls: a validation with a combined FEM/DEM approach,” International Journal for Multiscale Computational Engineering, vol. 16, no. 4, pp. 325–343, 2018. View at: Publisher Site  Google Scholar
 A. Cecchi and K. Sab, “A comparison between a 3D discrete model and two homogenised plate models for periodic elastic brickwork,” International Journal of Solids and Structures, vol. 41, no. 910, pp. 2259–2276, 2004. View at: Publisher Site  Google Scholar
 D. Baraldi and A. Cecchi, “Discrete approaches for the nonlinear analysis of in plane loaded masonry walls: molecular dynamic and static algorithm solutions,” European Journal of Mechanics—A/Solids, vol. 57, pp. 165–177, 2016. View at: Publisher Site  Google Scholar
 R. K. Livesley, “Limit analysis of structures formed from rigid blocks,” International Journal for Numerical Methods in Engineering, vol. 12, no. 12, pp. 1853–1871, 1978. View at: Publisher Site  Google Scholar
 C. Baggio and P. Trovalusci, “Limit analysis for notension and frictional threedimensional discrete systems,” Mechanics of Structures and Machines, vol. 26, no. 3, pp. 287–304, 1998. View at: Publisher Site  Google Scholar
 M. C. Ferris and F. TinLoi, “Limit analysis of frictional block assemblies as a mathematical program with complementarity constraints,” International Journal of Mechanical Sciences, vol. 43, no. 1, pp. 209–224, 2001. View at: Publisher Site  Google Scholar
 M. Gilbert, C. Casapulla, and H. M. Ahmed, “Limit analysis of masonry block structures with nonassociative frictional joints using linear programming,” Computers & Structures, vol. 83, no. 1314, pp. 873–887, 2006. View at: Publisher Site  Google Scholar
 F. Portioli, C. Casapulla, L. Cascini, M. D’Aniello, and R. Landolfo, “Limit analysis by linear programming of 3D masonry structures with associative friction laws and torsion interaction effects,” Archive of Applied Mechanics, vol. 83, no. 10, pp. 1415–1438, 2013. View at: Publisher Site  Google Scholar
 D. Baraldi, E. Reccia, A. Cazzani, and A. Cecchi, “Comparative analysis of numerical discrete and finite element models: the case of inplane loaded periodic brickwork,” Composites: Mechanics, Computations, Applications, An International Journal, vol. 4, no. 4, pp. 319–344, 2013. View at: Publisher Site  Google Scholar
 D. Baraldi, E. Reccia, and A. Cecchi, “DEM & FEM/DEM models for laterally loaded masonry walls,” in Proceedings of the 5th ECCOMAS Thematic Conference on Computational Methods in Structural Dynamics and Earthquake Engineering—COMPDYN 2015, pp. 2144–2157, Crete Island, Greece, 2015. View at: Google Scholar
 D. Baraldi, E. Reccia, and A. Cecchi, “In plane loaded masonry walls: DEM and FEM/DEM models. A critical review,” Meccanica, vol. 53, no. 7, pp. 1613–1628, 2018. View at: Publisher Site  Google Scholar
 A. Munjiza, The Finite/Discrete Element Method, John Wiley & Sons, Chicester, UK, 2004.
 F. Portioli, L. Cascini, C. Casapulla, and M. D’Aniello, “Limit analysis of masonry walls by rigid block modelling with cracking units and cohesive joints using linear programming,” Engineering Structures, vol. 57, pp. 232–247, 2013. View at: Publisher Site  Google Scholar
 P. B. Lourenço, J. G. Rots, and J. Blaauwendraad, “Two approaches for the analysis of masonry structures: micro and macromodeling,” Heron, vol. 40, 4 pages, 1995. View at: Google Scholar
 A. W. Page, “Finite element model for masonry,” Journal of the Structural Division, vol. 104, pp. 1267–1285, 1978. View at: Google Scholar
 P. De Buhan and G. de Felice, “A homogenization approach to the ultimate strength of brick masonry,” Journal of the Mechanics and Physics of Solids, vol. 45, no. 7, pp. 1085–1104, 1997. View at: Publisher Site  Google Scholar
 R. Luciano and E. Sacco, “Homogenization technique and damage model for old masonry material,” International Journal of Solids and Structures, vol. 34, no. 24, pp. 3191–3208, 1997. View at: Publisher Site  Google Scholar
 G. Milani, “Simple homogenization model for the nonlinear analysis of inplane loaded masonry walls,” Computers & Structures, vol. 89, no. 1718, pp. 1586–1601, 2011. View at: Publisher Site  Google Scholar
 E. Bertolesi, G. Milani, and P. B. Lourenço, “Implementation and validation of a total displacement nonlinear homogenization approach for inplane loaded masonry,” Computers & Structures, vol. 176, pp. 13–33, 2016. View at: Publisher Site  Google Scholar
 D. Addessi and E. Sacco, “A multiscale enriched model for the analysis of masonry panels,” International Journal of Solids and Structures, vol. 49, no. 6, pp. 865–880, 2012. View at: Publisher Site  Google Scholar
 I. Stefanou, J. Sulem, and I. Vardoulakis, “Threedimensional Cosserat homogenization of masonry structures: elasticity,” Acta Geotechnica, vol. 3, no. 1, pp. 71–83, 2008. View at: Publisher Site  Google Scholar
 G. Salerno and G. De Felice, “Continuum modeling of periodic brickwork,” International Journal of Solids and Structures, vol. 46, no. 5, pp. 1251–1267, 2009. View at: Publisher Site  Google Scholar
 M. L. De Bellis and D. Addessi, “A cosserat based multiscale model for masonry structures,” International Journal for Multiscale Computational Engineering, vol. 9, no. 5, pp. 543–563, 2011. View at: Publisher Site  Google Scholar
 A. Bacigalupo and L. Gambarotta, “Nonlocal computational homogenization of periodic masonry,” International Journal for Multiscale Computational Engineering, vol. 9, no. 5, pp. 565–578, 2011. View at: Publisher Site  Google Scholar
 D. Baraldi and A. Cecchi, “Discrete element model for inplane loaded viscoelastic masonry,” International Journal for Multiscale Computational Engineering, vol. 12, no. 2, pp. 155–175, 2014. View at: Publisher Site  Google Scholar
 F. da Porto, G. Guidi, E. Garbin, and C. Modena, “Inplane behavior of clay masonry walls: experimental testing and finiteelement modeling,” Journal of Structural Engineering, vol. 136, no. 11, pp. 1379–1392, 2010. View at: Publisher Site  Google Scholar
 B. Ghiassi, D. V. Oliveira, P. B. Lourenço, and G. Marcari, “Numerical study of the role of mortar joints in the bond behavior of FRPstrengthened masonry,” Composites Part B: Engineering, vol. 46, pp. 21–30, 2013. View at: Publisher Site  Google Scholar
 V. Gattulli, G. Lampis, G. Marcari, and A. Paolone, “Simulations of FRP reinforcement in masonry panels and application to a historic facade,” Engineering Structures, vol. 75, pp. 604–618, 2014. View at: Publisher Site  Google Scholar
 X. Wang, B. Ghiassi, D. V. Oliveira, and C. C. Lam, “Modelling the nonlinear behaviour of masonry walls strengthened with textile reinforced mortars,” Engineering Structures, vol. 134, pp. 11–24, 2017. View at: Publisher Site  Google Scholar
 V. G. Haach, G. Vasconcelos, and P. B. Lourenço, “Study of the behaviour of reinforced masonry wallets subjected to diagonal compression through numerical modelling,” in Proceedings of the 9th International Masonry Conference, Guimaraes, Portugal, 2014. View at: Google Scholar
 L. Miccoli, D. V. Oliveira, R. A. Silva, U. Müller, and L. Schueremans, “Static behaviour of rammed earth: experimental testing and finite element modelling,” Materials and Structures, vol. 48, no. 10, pp. 3443–3456, 2015. View at: Publisher Site  Google Scholar
 C. B. de Carvalho Bello, A. Cecchi, E. Meroi, and D. V. Oliveira, “Experimental and numerical investigations on the behaviour of masonry walls reinforced with an innovative sisal FRCM system,” Key Engineering Materials, vol. 747, pp. 190–195, 2017. View at: Publisher Site  Google Scholar
 C. B. de Carvalho Bello, G. Boscato, E. Meroi, and A. Cecchi, “Nonlinear continuous model for three leaf masonry walls,” Construction and Building Materials, vol. 244, p. 118356, 2020. View at: Publisher Site  Google Scholar
 G. Bartoli, M. Betti, P. Biagini et al., “Epistemic uncertainties in structural modeling: a blind benchmark for seismic assessment of slender masonry towers,” Journal of Performance of Constructed Facilities, vol. 31, no. 5, Article ID 04017067, 2017. View at: Publisher Site  Google Scholar
 D. Baraldi, C. B. de Carvalho Bello, A. Cecchi, E. Meroi, and E. Reccia, “Nonlinear behavior of masonry walls: FE, DE, and FE/DE models,” Composites: Mechanics, Computations, Applications: An International Journal, vol. 10, no. 3, pp. 253–272, 2019. View at: Publisher Site  Google Scholar
 T. M. Raijmakers and A. T. Vermeltfoort, “Deformation Controlled meso shear tests on masonry piers,” TNOBOUW/TU Eindhoven, Build. and Constr. Res., Eindhoven, Netherlands, 1992, Report B921156. View at: Google Scholar
 M. Petracca, L. Pelà, R. Rossi, S. Zaghi, G. Camata, and E. Spacone, “Microscale continuous and discrete numerical models for nonlinear analysis of masonry shear walls,” Construction and Building Materials, vol. 149, pp. 296–314, 2017. View at: Publisher Site  Google Scholar
 J. G. Rots and J. Blaauwendraad, “Crack models for concrete: discrete or smeared, fixed, multi‐directional or rotating,” Heron, vol. 34, no. 1, pp. 1–59, 1989. View at: Google Scholar
 J. G. Rots, “Computational modeling of concrete fracture,” Delft University of Technology, Delft, Netherlands, 1988, Ph.D. thesis. View at: Google Scholar
 TNO DIANA, DIsplacement Method ANAlyser, Release 9.4, User’s Manual, 2010, https://dianafea.com/manuals/d94/Diana.html.
 R. van der Pluijm, “Material properties and its components under tension and shear,” in Proceedings of the 6th Canadian Masonry Symposium, Saskatoon, Canada, June 1992. View at: Google Scholar
 R. van der Pluijm, “Shear behaviour of bed joints,” in Proceedings of the 6th North American Masonry Conference, pp. 125–136, Philadelphia, PA, USA, 1993. View at: Google Scholar
 J. Lopez, S. Oller, E. Oñate, and J. Lubliner, “A homogeneous constitutive model for masonry,” International Journal for Numerical Methods in Engineering, vol. 46, no. 10, pp. 1651–1671, 1999. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2020 Daniele Baraldi et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.