Abstract

Firstly, an updated numerical model for the numerical analysis of planar unreinforced and confined masonry walls is presented. The model can simulate the main nonlinear material effects of masonry and reinforced concrete. A simplified anisotropic constitutive model for masonry is developed and presented. The criteria for the limit bearing capacity and collapse of masonry are separately defined for normal stresses only and for normal and shear stresses. The presented numerical model is verified and used to analyse the anisotropy effect of masonry on the behaviour of unreinforced and confined two-story anisotropic masonry walls with different coefficients of anisotropy, wall lengths, and quality of masonry under horizontal static forces. The influence of the anisotropy coefficient of the masonry on the response of the walls is discussed in detail, and the main conclusions are given.

1. Introduction

Unreinforced and confined masonry walls are complex structures that are primarily composed of masonry and concrete elements, such as foundations, tie columns, and tie beams. The masonry is formed by masonry units, which are connected by binder material in horizontal and vertical joints. Due to the different properties of vertical and horizontal joints, and especially because the properties of masonry units often differ in the vertical and horizontal directions, the masonry is always anisotropic. If the plaster is included in the bearing capacity of masonry, especially if it is reinforced, the structure of masonry walls becomes even more complicated. Uplifting and slippage of the foundation from the ground often occur under the action of horizontal forces on the walls, which also complicates reliable analysis of these complex structures. The main subject of this paper is a numerical simulation of the anisotropy effect of masonry on the behaviour and bearing capacity of unreinforced and confined masonry walls.

Numerous numerical models for the simulation of the anisotropic properties of masonry have been developed. One of the first nonlinear models for masonry was proposed by Page [13], which considered the masonry to be a two-phase material composed of bricks and mortar. Dhanasekar et al. [4] defined a failure surface of masonry as composed of three elliptical cones. Andreaus [5] carried out a relatively comprehensive study, taking into account different failure modes of masonry using three failure criteria, namely, modified Mohr–Coulomb friction law, maximum tensile strain criterion, and maximum compressive stress criterion. Lourenco et al. [68] defined a limiting surface to describe the orthotropic behaviour of masonry with a material admissible field bounded by a Hill-type yield criterion for compression and a Rankine-type yield criterion for tension. Syrmakezis and Asteris [9] used a third-order polynomial to mathematically describe the failure surface. Berto et al. [10] developed an orthotropic damage model based on the definition of four independent internal damage parameters, two parameters in compression and two parameters in tension. Calderini et al. [11] described the anisotropic damage process and the dissipative behaviour of masonry structures under static incremental and dynamic loads. Pela et al. [12, 13] presented a plane-stress damage macromodel for orthotropic masonry based on the concept of mapped tensors from the orthotropic real space. Lishak et al. [14] developed phenomenological masonry failure criteria for use in numerical calculations based on models where the masonry is considered an orthotropic homogeneous material. Ghiassi et al. [15] presented an orthotropic macromodel for nonlinear static analysis of masonry walls. Grande et al. [16] presented a frame element model for evaluating the nonlinear response of unstrengthened and FRP-strengthened masonry panels subjected to in-plane vertical and lateral loads. Basaran et al. [17] conducted a study whose purpose is the improvement of the behaviors of walls constructed with masonry bricks using reinforced plaster mortars. Calio and Panto [18] used a macroelement modelling approach of infilled frame structures. Validation of macroelement models and cost analysis for unreinforced and confined masonry buildings were considered by Marques and Lourenco [19]. Asteris et al. [20] used a macromodelling approach for the analysis of infilled frame structures. Berardi [21] used a micromechanical procedure in order to evaluate the initiation of damage and failure of masonry with in-plane loads. Ridwan et al. [22] proposed a new simplified formula for the calculation of the equivalent modulus of elasticity of masonry structures. Panto et. al. [23] used a 3D discrete macroelement for modelling the out-of-plane behaviour of infilled frames. Yang et al. [24] proposed an equation for calculating the shear strength of masonry with precast concrete antishear blocks based on the experimental results. Yang et al. [25] proposed a simple and rational stress-strain relationship model applicable to brick masonry under compression.

However, there is still no such numerical model for the analysis of masonry structures that would be appropriate for extensive engineering applications, i.e., which would be simple enough and able to simulate the main nonlinear effects of the masonry structures. Hence, such a model should adequately simulate not only the nonlinear behaviour and anisotropy of masonry but also the real behaviour of reinforced concrete elements and soil-structure coupled systems.

In this paper, a quite simple numerical model for analysis of planar unreinforced and confined masonry structures is briefly described. The model can simulate the main nonlinear effects of masonry, reinforced concrete, and soil, as well as the soil-structure interaction and large displacements of the structure. The model is intended for engineering applications. The emphasis on improving the constitutive masonry model developed by Radnić et al. [26] is the primary purpose of this paper, with the intention of modelling the masonry anisotropy. The directions of the masonry anisotropy correspond to the horizontal and vertical joints of the masonry units. Corresponding stiffness, strength, and limit strains of masonry are modelled for the main direction of anisotropy. Criteria for yielding (limit bearing capacity) and failure (collapse) of masonry are separately defined for normal stresses only and for normal and shear stresses. The yielding of masonry in the compression and cracking of masonry in tension are modelled, as well as the opening and closing of smeared cracks at repeated loading. Constitutive models for reinforced concrete, soil, and contact elements were taken from the study of Radnić et al. [26]. The updated numerical model is verified on the experimental test [27] and numerical results [28].

Using the updated numerical model, the behaviour and bearing capacity of unreinforced and confined two-story masonry walls with different lengths, qualities of masonry, and coefficients of anisotropy are analysed under horizontal static forces. It is assumed that the wall foundations stand on a rigid base with the possibility of uplifting. The walls were loaded with a constant vertical load and various horizontal forces at the floor levels. The force successively increased until wall failure. The effect of masonry anisotropy on the behaviour and limit bearing capacity of the walls is studied. The conclusions of the analysis are presented at the end of this paper.

2. Numerical Model

A detailed description of the basic finite element numerical model can be found in the study of Radnić et al. [26]. In this paper, the constitutive model of masonry was extended to include the effect of anisotropy of masonry and a more adequate model of shear effect on bearing capacity of masonry.

A simplified macromodel for masonry is adopted that is intended for practical applications and able to simulate the main nonlinear effects of this material. At the macrolevel, the masonry is approximated by a representative material whose mechanical properties can adequately describe the actual complex masonry properties. This approach allows the use of large finite elements (rough discretization) and significantly reduces the number of unknown variables in comparison with a micromodel of masonry (i.e., this approach rapidly accelerates the structure analysis).

An orthotropic model for masonry is assumed, which is analogous to a previously developed model by Smilovic et al. [29] for masonry exposed to dynamic loading. The main directions of the orthotropy are horizontal (h) and vertical (), which is the usual case in real masonry walls. For each direction of orthotropy, different compressive strengths of masonry , tensile strengths , moduli of elasticity , limit compressive strains , and limit tensile strains can be defined. Wherein, denote strains at the cracks occurrence , and are the fictitious strains after which the masonry has no tensile stiffness. The coefficient of anisotropy (orthotropy) is defined by

It should be noted thatwhere and are Poisson’s ratios in the horizontal and vertical directions, respectively.

The criteria of limit bearing capacity of masonry were considered separately for normal stresses only and for normal and corresponding shear stresses. Graphical presentations of the adopted orthotropic constitutive masonry model in compression and tension are given in Figure 1, where and are the normal stresses.

As shown in Figure 1, the effect of biaxial compressive stresses on the limit bearing capacity of masonry is disregarded. As the model is intended for practical use, a small number of parameters were assumed that can be easily determined from classical uniaxial tests (strengths, limit strains, and modulus of elasticity).

In biaxial compression, elastic-plastic behaviour with linear unloading is assumed. Elastic behaviour is assumed when the stresses in masonry are smaller than the uniaxial compressive strength. Then completely plastic behaviour is assumed. The crushing of masonry is assumed when the strain of masonry reaches the maximum compressive strain from the uniaxial test.

In biaxial tension, elastic behaviour is adopted until the tensile strength of masonry is reached. Then, it is assumed that smeared cracks occur in the plane perpendicular to the direction of tensile stresses. The tensile stiffness of the cracked masonry is simulated by a gradual decrease in the tensile stress components perpendicular to the crack plane in accordance with the stress-strain relation from the uniaxial tension test.

The opening and closing of cracks are modelled. The transfer of compressive stress across a fully closed crack is modelled as in masonry without cracks. After crack reopening, the tensile stiffness of the cracked masonry perpendicular to the crack plane is taken to be zero. The shear stiffness of the cracked masonry is modelled by reducing the shear modulus depending on the value of the tensile masonry strain perpendicular to the crack plane.

In addition to tension (cracking) and compression (crushing), the collapse of the masonry due to shear stress is modelled. The criterion for the yielding and failure (collapse) of masonry under normal and shear stresses is defined according to Figure 2. The model is simple and includes the relation of normal and shear stresses/strains according to the circular equation. Therefore, the masonry yield criterion is defined in a function of stresses and the failure criterion as a function of strains. A simple relation between normal and shear stresses/strain was adopted for use in practical applications, where basic parameters of masonry are usually not experimentally determined.

The shear modulus of masonry Gm is given by

The adopted isotropic concrete model in Section 4 is described in the study of Radnić et al. [26, 30], which can simulate yielding of concrete in compression and cracking in tension. The adopted constitutive model for contact elements is described in the study of Radnić et al. [26], which can simulate separation and slippage on contact surface. The presented anisotropic masonry model or isotropic concrete model can be used for the soil model, with corresponding material parameters for soil.

3. Validation of the Numerical Model

The numerical model presented in Section 2, previously verified in its basic form [26], is validated based on the results of experimental tests at ETH Zürich by Ganz and Thürlimann [27] and numerical results by Lourenco [28]. The geometry and loading of walls W1 and W2 are presented in Figure 3. The experimentally determined material properties of the masonry are shown in Table 1 [27]. The walls were subjected to a uniform axial load followed by an incremental horizontal force (F) at the top of the slab that caused a horizontal displacement (d). The only difference between walls W1 and W2 is the initial axial force applied to the walls, which resulted in different behaviour and failure modes. The initial vertical load was q = 0.61 MPa (415 kN) for wall W1 and q = 1.91 MPa (1287 kN) for wall W2. Wall W1 showed a ductile response with tensile and shear failure along the diagonal stepped cracks. The behaviour of wall W2 began with relative ductility and was followed by brittle cracking [27]. The adopted parameters of masonry for the presented numerical model are presented in Table 2.

Only some of the numerical results are presented. The experimentally determined force (F) versus displacement (d) is presented in Figure 4. The experimentally determined crack patterns for wall W1 at a displacement of 12.0 mm are shown in Figure 5. The numerical crack patterns for wall W2 at a displacement of 8.0 mm are presented in Figure 6.

The numerical results of the presented numerical model are in good agreement with the experimental results of Ganz and Thürlimann [27] and the numerical results of Lourenco [28].

4. Numerical Tests of the Anisotropy Effects on the Behaviour of Unreinforced and Confined Masonry Walls under Horizontal Static Load

Using the presented numerical model, the effect of masonry anisotropy on behaviour of single two-story unreinforced and confined masonry walls of different lengths under horizontal static force was investigated. The reason for this is the fact that masonry blocks in practice usually have vertical cavities and generally different load bearing capacity and deformability in the vertical and horizontal directions. We believe that this study contributes to the article quality and that it can be useful for understanding the effects of masonry anisotropy in practice.

4.1. Basic Data of the Analysed Walls

Two-story unreinforced and confined masonry walls with different lengths B (3 m, 6 m, and 12 m) were analysed (Figure 7). All walls were 0.24 m thick with 0.6 m wide concrete foundations. The foundations were reinforced with 3Ø10 longitudinal bars in the upper and bottom zones and Ø8 transverse stirrups with 250 mm spacing. The tie beams were reinforced with 4Ø10 bars and tie columns with 4Ø12 bars. The stirrups of the tie beams and tie columns were Ø6 at a 250 mm spacing.

4.2. Finite Element Discretization and Adopted Material Parameters

Finite element discretization of the walls is presented in Figure 8. It is assumed that the wall foundation stands on a rigid base with the possibility of uplifting. For this purpose, contact elements are placed below the foundations. Masonry with high strength and stiffness (stiff masonry) and masonry with low strength and stiffness (soft masonry) were considered. The adopted basic material parameters of the masonry (modulus of elasticity, shear modulus, compressive strength, and tensile strength) for the considered coefficient of anisotropy ca are shown in Table 3. The effect of shear stresses on the bearing capacity of masonry is not included (model according to Figure 2) in the performed analysis of masonry walls for clarity of presentation of the anisotropy effect of masonry under normal stresses only.

The walls were assumed to be in a plane-stress state. A rather rough finite element mesh was used, which should not affect the conclusions of the study. 8-node concrete elements and 3-node reinforcement elements were used. Full compatibility of displacement at the function of concrete and reinforcement was assumed. An elastic-fully plastic model of concrete in pressure was used. A smeared cracks model for concrete in tension was adopted. The mechanism of opening and closing of cracks in concrete and the tensile and shear stiffness of cracked concrete are modelled. Contact elements were taken at the bottom of the foundation.

The adopted basic material parameters of concrete, reinforcement, and contact elements are shown in Table 4. First, the initial displacements and internal stresses for the dead weight and vertical load q of the walls were calculated. Then, the walls were incrementally loaded with horizontal force F at the floor levels until wall failure.

4.3. Numerical Results

Some results of the numerical analysis performed for unreinforced and confined masonry walls with stiff and soft masonry are presented in Figures 915. When evaluating the presented diagrams, it should be borne in mind that the foundation lifting on the tensile side of the wall was modelled and that, in many cases, there was a loss of stability of the wall due to its overturning and not to wall collapse due to excess limit stresses/strains in masonry, concrete, or reinforcement (especially with narrower confined walls).

As shown in the presented figures, limit strengths and displacements (deformability), as well as the stresses in masonry, concrete, and reinforcement, depend on the coefficient of anisotropy ca. As the anisotropy of masonry is higher (smaller ca), the difference in these values between walls with isotropic and anisotropic masonry is also higher. With the same coefficient of anisotropy, the behaviour of the masonry wall depends on the type of the walls (unreinforced and confined), on the height-length ratio of the wall, and on the quality of the masonry (stiff masonry and soft masonry).

Longer walls have a greater anisotropy effect relative to shorter walls. Soft masonry has a greater anisotropy effect relative to stiff masonry. The anisotropy effect is more affected by the horizontal displacement (deformability) of the wall and less affected by the ultimate bearing capacity. Confined masonry walls have a higher anisotropy effect relative to unreinforced masonry walls. For high values of the coefficient ca (lower anisotropy), the difference in behaviour between anisotropic walls and isotropic walls is small.

5. Conclusions

The presented numerical model for nonlinear analysis of planar masonry structures, which is upgrading of the previously developed model [26], can simulate their main nonlinear behaviour. The model is simple and primarily intended for practical application. Further verification of this numerical model is still required.

The ultimate bearing capacity, displacements, and stresses in considered confined and unreinforced masonry walls under horizontal static loads depend on the anisotropy of masonry, but not a lot. The reason is that the reduced masonry strength in the horizontal direction corresponds to significantly lower wall stresses in that direction, compared to vertical direction. As the anisotropy increases, the difference in the above wall values increases. The behaviour of masonry walls with low anisotropy is similar to the behaviour of walls with isotropic masonry. The anisotropy effect also depends on the type of the wall (unreinforced and confined), the quality of the masonry (stiff masonry and soft masonry), and the height-length ratio of the walls, but also not significantly. Confined masonry walls have a greater anisotropy effect than unreinforced masonry walls due to increased horizontal stresses in walls due to tie beams and tie columns. The anisotropy effect depends very little on the quality of the masonry. Longer walls have somewhat greater anisotropy effect than shorter walls. The aforementioned facts are the reason that, due to the tendency to reduce weight, masonry blocks with vertical cavities have long been used in practice. The use of soft masonry, especially those with a higher anisotropy, should be excluded in practice.

Data Availability

All data underlying the findings of the study are presented in this paper.

Conflicts of Interest

The authors declare that they have no conflicts of interest.