Compressive strength of heterogeneous masonry walls containing blends of brick types

The study presents a systematic approach for the evaluation of the compression strength of masonry walls composed of heterogeneous mixes of different types of blocks. First of all, the mechanics of a compressed heterogeneous masonry stack is investigated through a series of experimental tests and Finite Element models, then it is reviewed and discussed. Then, the problem of deriving the necessary material parameters entering the Hilsdorf formula is addressed. Solutions for the correct evaluation of the lacking data are presented based on the existing literature data. Finally, the well-known Hilsdorf formula is extended to the field of block blends with different mechanical properties. A deep experimental investigation on stacks and wallets made with fired clay, limestone and sandstone blocks is introduced for the first time. The comparison of the experimental data with the proposed theory points out the very good predictive capability of the extended Hilsdorf formula derived herein.


Introduction
The knowledge of the mechanical properties of an existing building is at the base of its structural assessment. Especially for masonry structures, several models have been developed over the years trying to predict their global capacity. The technological improvement of the last twenty years led to the development of sophisticated numerical models that needs a large computational time to provide the problem solution. Besides, sometimes it happens that even no convergence or mindless results could be achieved. Both numerical and analytical models are often supported by experimental tests, which can prove the accuracy of the proposed model. Lofti et al. [1] and Lourenco et al. [2] reproduced numerically the behaviour of a masonry panel considering the exact geometry of both components by introducing the concept of the interface model. However, the method is nowadays far from practical use when a large number of units is involved. A few years later, Zucchini and Lourenco [3] discussed several existing results concerning methods for masonry homogenization and validated experimentally their micro-mechanical model based on Representative Volume Elements.
On the other hand, formulas and relations could be even found experimentally from the observation of the physical behaviour and regression analysis of the data sets. Recently, Ferretti [4] calibrated a dimensionless power model for predicting the compressive strength of solid-clay-brick masonry walls by collecting a large base of data from the available literature in which both stacks and wallets were included. The calibration of the model involving dimensionless parameters extended the formula proposed by the EC6 [5] to more general instances made of units of only one material. However, it is to stress that, unlike some non-European studies, the largest effort has been put in masonry systems in which the mortar is weaker than units, and therefore in the present study cases with strong mortar and weak units will be disregarded, although the proposed theory applies in principle even to those sets.
Concerning the compressive strength of a layered stack composed of two alternating materials, it has been widely studied since the beginning of the twentieth century (for a list of contributions see Proske and Van Gelders [6]). The solution introduced in 1969 by Hilsdorf [7] and partially described earlier in a work of Haller in 1958 [8], based on the study of the mutual horizontal stress under vertical compressive loading, answers the question: ''Why does a stack of two materials fail under compression at a stress larger than the minimum compressive strength of the two materials?'' A detailed discussion of the compressive strength of brick and mortar stacks can be found in Aprile et al. [9] and Como [10]. The equilibrium and compatibility equations are: where the subscript v and h stand for the vertical and horizontal directions, b is for bricks, m is for mortar, and t is the thickness of the material layers. In the case of masonry wallets (Fig. 1A), the horizontal direction is affected by the presence of the head joints. These weak elements increase the tensile stress required on bricks for fulfilling the equilibrium with the compressed mortar layers, compared to the case of masonry stacks (Fig. 1B). The resulting stress increment is roughly proportional to the number of mortar layers that interact with each brick. Therefore an influence factor k (k = 2 for wallets and k = 1 in the case of stacks) is added in the formula (1) stating the horizontal equilibrium. The solution of the problem gives the ratio of the horizontal and vertical stresses: In which m is the Poisson ratio, q is the ratio of the elastic moduli E m /E b and g is the ratio of the thicknesses t m /t b . Finally, the masonry compressive strength can be computed through the Mohr-Coulomb criterion: where the subscript M stands for masonry and f cb and f tb are the compressive and tensile strengths of bricks. The collapse is produced by vertical cracks running on the blocks whilst the mortar is in a triaxial compression state, which maintains it in an elastic range well beyond the uniaxial compressive strength [11]. Several studies [12][13][14][15][16] have demonstrated that a proportional loading with confinement stresses of 30-40% of the vertical compressive stress, can allow cement mortars to raise stresses over 10 times the uniaxial strength. However, real ancient masonry panels include normally several types of bricks, due to the reuse of demolition parts of other buildings. Stonemasonry walls include often layers of fired clay bricks [10]. Even adobe brick walls show inclusions of fired bricks, mainly on the lower storeys or around the openings. It is thus necessary to understand what is the compressive capacity of a wall containing brick blends.
A first experimental and theoretical study was presented in Benedetti [17], in which tests conducted on model wallets made with mixed adobe and fired bricks were presented and analysed. Afterwards, several experiments were performed at the University of Firenze [18,19] by mixing three different types of materials (namely fired bricks, limestone and sandstone blocks), but the relevant interpretation theory was not set out.
In what follows the experiments are reconsidered and a general theory able to interpret the behaviour of walls with general brick admixtures is presented and verified against the experimental data.

Outline of the study, objective and original contribution
The research is based on the following steps: (1) Definition of the material parameters and the type of analysis relevant for the problem under study, (2) Presentation of the experimental results of two wide investigations on the compressive strength of masonry stacks and wallets including brick admixtures, (3) Set up of the Hilsdorf Extended Theory (HET) and its verification against the presented experimental data.
This work wants to provide an analytical tool for a rapid and direct evaluation of the strength of heterogeneous masonry walls. Assorted brick materials and their random position in the same wall are typical of old masonry buildings for which the manufacturing followed the local availability of the materials, especially after exceptional events, like earthquakes or floods, when the reconstruction needs to be fast for restoring the normal life quickly. The local behaviour of blended masonry observed in some specimens can be extended to existing blended walls explaining local failures or estimating their mechanical properties having the wall pattern.
Most of the experimental studies presented in the literature use only a reduced set of mechanical parameters for mortars and bricks. Therefore the first contribution of this study is aimed at introducing rules for the derivation of the missing data and the selection of the structural problem type.
Then, the second contribution is the extension of the Hilsdorf formula to general compositions of bricks with different properties and dimensions. In most of the numerical benchmark and case history analyses present in the literature, the masonry properties are assumed constant over the whole structure, and sometimes extracted by similarity from other studies. The proposed approach fills the gap of finding homogenized materials' properties to be assigned to local elements in most of the mesoscale and large scale numerical models built for a close structural assessment of real heterogeneous structures. However, to maintain the study relevant for practical applications, in this first derivation, the effect of the wall spatial organization and brick clustering is introduced through the scalar influence factor k introduced in Eq. (2), which involves only the number of mortar courses interacting with each bricklayer. 2 Parameters influencing the brick-mortar interaction problem

Definition of the elastic interaction problem
Masonry walls can be organised with many different bonding patterns, which in principle could produce varying levels of mortar confinement, but in any case, a set of parallel stacks can be identified in every wall. Therefore, even if the problem is intrinsically 3D, it can be studied as a 2D problem either for in-plane stress or in-plane strain condition. A detailed elastic study has been worked out, assuming perfect bonding between brick and mortar layers, for a 2D infinite masonry stack and a 3D crossing bond column of the same size. Since the problem is assumed linear elastic, no interface constitutive laws were considered for the connection zones of the different material layers. The used material data are: In Fig. 2 the Finite Element Models used in this study are presented, in which a quarter of a 3D column with a section of 250Á250 mm 2 has been modelled, together with the equivalent 2D stack of the same size. Concerning the 3D model, 8-nodes solid brick elements were used either for mortar or masonry units, while 4-nodes elements were chosen in the case of the 2D FE model, leading to the linear description of the stresses as the best compromise between the model accuracy and the computational cost. In detail, 3 and 12 solid bricks layers are used for modelling the thickness of mortar and brick courses in the 3D model, while 24 and 5 plane elements layers are used in the 2D model. Both models were created by using the Italian release of the STRAND Finite Element software. Figure 3 shows the comparison among the horizontal stresses obtained in 3D, 2D plane stress and 2D plane strain models under the same unitary load pressure equal to 1 MPa at the top and the bottom of the stacks. Since the two models are self-balanced along with the vertical direction, only lateral restraints have been considered. In particular, the 2D mesh is a subset of an infinite stack with a whole model width of 250 mm, and a thickness of 250 mm. The restraints simply eliminate the possible rigid body motions of the 2D displacement field. The 3D model is split by the two symmetry planes and includes, in addition to the restraints eliminating the rigid body motions, the horizontal restraints forcing a zero displacement across those planes.
As is fairly evident, the 3D and 2D plane strain analyses match very well, while 2D plane stress has a large error in the mortar confinement estimation. This result is well documented in the literature [20]. Therefore, in the study of the interaction in brick blends a 2D plane strain analysis will be carried out. The definition of compressive strength for stacks of brick blends requires not only a precise knowledge of the compressive and tensile strength of the involved bricks but also representative values of mortar and brick Poisson's ratios, very elusive and difficult to be measured experimentally.
Concerning the tensile strength of the bricks, it is to cite that at least three different values can be obtained from laboratory tests, namely flexural, splitting, or direct tensile strengths. Thus, conversion rules are necessary, since the flexural strength is larger than the uniaxial tensile strength, while the splitting strength is smaller.
The flexural strength depends on the shape of the specimen, but generally, its value is 50% above the uniaxial tensile strength. Regarding the splitting strength f t,sp , by assuming a ratio z of the vertical to horizontal stresses in the cylinders subject to splitting, a Mohr-Coulomb criterion can allow defining the strength conversion [21]: According to Li [22], the factor z can be computed in the centre of the cylinder as a function of the loading arc opening h, normally resulting in values very near to 3.0: However, in many studies, the tensile strength of bricks is not even taken into consideration, though it is fundamental for a Mohr-Coulomb analysis. Many studies can be found in the literature dealing with the ratio of tensile to compressive uniaxial strength of mortars, bricks and stones [23,24]. A wide and detailed analysis of experimental data led to defining a very typical functional form that links the two strengths: Usually, the coefficient a seems to vary with the material type (lime mortar, concrete, brick, stone, etc.), the coefficient b is almost constant among materials and holds 2/3 [25].
As an example, in the following Table 1 the specimens' geometry and the results of an experimental campaign, done in Forlì, Italy [26] on XVII century bricks, are presented. In particular, two specimens were collected for each brick of interest: a solid prism was cut from the original brick and then tested in compression evaluating its compressive capacity f c , while a brick core was drilled from the same unit estimating its tensile strength through the splitting test f sp . The compressive and splitting strength are used in defining the pure tensile strength of the brick as in formula (4). Then this value is compared with formula (6). The agreement between the two values is noticeable.
The evaluation of the Poisson ratio is problematic by either using experimental data or looking at the results of mechanical theories. A deep and detailed discussion of the importance and the criteria for the evaluation of this parameter can be found in rock mechanics contributions [16,27]. Zhang presented in Chinese [28] a detailed discussion of the relationship of the Coulomb friction angle /, the rigidity R and the Poisson ratio m, which was later reconsidered by Logo [29]. It is to point out that the rigidity ratio has values fairly equal to the Hoeck-Brown constant m i [30]: Anson and Newman [31] discussed the difference between Poisson ratio values obtained through static and dynamic experiments and the effect of the porosity on the detected values. In fact, tensile and compressive strengths, elastic moduli and Poisson ratio are all approximately inversely proportional to the porosity of the material.
It is well known that in uniaxial compression tests the Poisson ratio (measured as the ratio of the horizontal to the vertical strains), is weakly decreasing in the elastic range, and then is rapidly increasing with the onset of the microcracking, finally even overpassing the limit value 0.5 near the collapse [21,32].
However, as discussed before, the masonry mortar courses are mainly in a triaxial compression state in which microcracking is hindered so that the Poisson ratio has only marginal variations with the confinement ratio [15,33].
Therefore, since the mortar-brick interaction is effective until the cracking of bricks or the crushing of mortar, it is natural to assume that the Poisson ratio is holding constant in this pre-failure range, and its value is linked (through porosity) to the strength ratio R.
Equation (8b, c) give values that better represent the behaviour detected with the experiments.
Therefore, even when only the compressive strength is known, it is possible to gather the complete set of brick and mortar parameters by the following procedure: (a) The tension strength is evaluated from splitting tests or by using Eq. (6); (b) The Coulomb friction angle is determined from the ratio R; (c) The Poisson ratio is obtained from Eq. (8b).
Finally, once the U factor is computed, the compressive strength of masonry is obtained via the Hilsdorf theory. The value is computed with formula (4) by using z = 3 b The value is computed by using formula (6) with a = 0.25

Elastic analysis of heterogeneous masonry stacks
A masonry stack composed of a blend of different bricks is neither a random material nor a regularly layered composite. It is thus significant to examine horizontal stress values arising from the interaction of mortar with blends of different bricks, which can lay either as staggered clusters or in interspersed sequences.
In what follows, according to the previous results of the finite element investigations, several plane strain stacks composed of ten courses of three typical brick materials will be investigated. The mechanical properties are chosen, on average, according to Augenti et al. [37]: The selected thickness of the layers produces two g values, namely 0.364 (blocks 55 mm, mortar 20 mm), and 0.25 (blocks 60 mm, mortar 15 mm). In the following Fig. 5, the stacks are presented.
The nomenclature is as follows: A stands for alternate units and B for clustered sets (sequence of clusters of the same brick type), C for clay units, L for limestone units and S for sandstone units. The number is the percentage of every material composing the stack, 2.75 and 4.00 are the unit/mortar thickness ratios used for the simulation.
Even some compositions of three materials were examined. In Fig. 5d the three studied sample block sequences are listed.
Since the bonding pattern could be different according to the content of different bricks and their position, the analysis of any possible combination is not achievable. Therefore, considering alternated and clustered brick distributions as the two limit bounds of all the possible configurations, other layouts are surely included in the study.
The analysis was carried out as plane strain on finite element models with unitary vertical stress. The models were composed of 4000-7000 square fournode finite elements with a side of 5 mm belonging to the whole model side of 120 mm, and then the symmetry conditions were set at the two end faces.
To compare the horizontal stress in mortar courses and different brick materials, the stress values at the symmetry plane of the stacks were extracted from the models. However, due to the local combinations of block types, the horizontal stresses are mildly varying among different mortar and block positions, and even in different layers of the same course.
Therefore average stress values were obtained by summing up the stress resultants for every material and dividing it by the given material area crossing the symmetry plane.   The effectiveness of the averaging procedure has been checked by computing the equilibrium on the symmetry plane of the stress resultants of the different material groups. A fairly good error of less than 0.1% of the tension and compression resultants has been obtained.
The following Fig. 6a, b present the results for the mortar horizontal stress and the clay brick horizontal stress, as it is computed in the various models by varying the sandostone or limestone content and the ratio h = t b /t m . Given the stresses from the FE model, the adimensional stress ratios are then obtained by dividing by 1 MPa as a reference stress value.
The numerical results clearly point out that the mortar horizontal confinement stress r hm is approximately independent of the bond layout of the stacks, while the average stiffness of the blocks, given as the homogenized elastic modulus of the components through their volumetric ratios, is resulting in a linear influence on it. Thus, since clay is stiffer than limestone but softer than sandstone, the two diagrams of Fig. 6a show respectively ascending and descending paths with the increase in brick content. On the other hand, the brick tensile stress r hb is depending strongly on both the bonding pattern and the thickness ratio h, even if the stress ratios are contained in a narrower range.
Since the collapse of masonry is mainly linked to the combined compressive-tensile biaxial stress of the blocks, obtaining a well-approximated evaluation formula for the block horizontal stress is an important task.
The performed numerical experiments point out that the mortar stress ratio U is strongly dependent on the homogenized properties of the block blend, while the block distribution has almost no influence. Therefore it is a realistic assumption that the collapse of a heterogeneous panel will be driven by the averaged equilibrium of the horizontal stresses in the elastic range until a diffused cracking of the blocks will affect the local stability of the separated stacks, leading to failure.
Therefore, in the following, the extended Haller-Hilsdorf theory (HET) will be discussed and verified avoiding any complex non-linear numerical simulation and considering only an elastic-brittle behaviour for the bricks. The proposed extension of the theory and the procedure for the statement of the missing parameters will be verified against experimental results of masonry panels realized with both homogeneous and heterogeneous blended block systems.

Experimental compressive strength of wallets with blended blocks
In what follows, the results of a large experimental investigation carried out at the University of Firenze is discussed [18,19]. The results concern mainly the compressive strength of scaled-down stacks and wallets (Figs. 7 and 8) composed of fired clay bricks (C) and units sawn from lime-(L) and sandstone (S) blocks. The mortar used for building each stack is a weak lime mortar according to what is typically present in historical constructions.  The stacks were composed of 6 blocks of 30Á30Á14 mm 3 sawn bricks, and 7 layers of mortar of 2.5 mm thickness, obtaining thus scaled down 1:4 models of half brick stacks.
The wallets had dimensions of 225Á63 mm 2 and a height of 264 mm and were built in a running bond pattern, with 16 layers of bricks and 17 mortar courses of 2.5 mm thickness. Displacement-controlled compressive tests were performed up to the panels' failure collecting the load, the four corner vertical displacements, and six vertical strains on both faces. The strains were measured by placing X strain gauges on each sample according to the setup illustrated in Fig. 8. Table 2 lists the mechanical data of the used materials.
In blended panels, since the block distribution can affect the results, similar panels were erected by interchanging the block type (CS-SC, CL-LC, SL-LS).
Two batches of three stacks (sets I and III, samples 1-3) and two sets of single or two wallets (sets II and IV, samples 1-2) were built and then tested in compression. The two sets III and IV aimed at completing and checking the findings of the previous series with some more combinations. The following Table 3 lists the results of the compression tests.
In what follows the theory extending the Hilsdorf formula to wallets with every blend of blocks will be presented in detail, and the comparison with the experiments will allow checking its effectiveness in predicting the strength of real walls.
Several techniques exist concerning the assessment of existing masonry walls, and usually, minor destructive tests are preferred to other methods due to their low invasiveness. For instance, Pelà et al. [38,39] and Segura et al. [40] inferred the strength of existing walls or masonry panels by testing under compression masonry cores and stacks even considering samples with different shapes and sizes. Concerning the experimental campaign discussed in this work, while the effects of a scaling procedure on the mechanical properties can be considered negligible for pure compressive tests [41], a linear relation was found among the strengths of stacks and wallets by Thamboo et al. [42], showing a ratio near to 4/3 which can be explained in term of the influence factor k already introduced and discussed. Thus, direct comparisons among the results of different test setups should be done with great caution.

The extended version of Hilsdorf masonry compressive strength
Several studies have compared experimental values of masonry wall compressive strength with available prediction formulas. Drougkas [43] presented a list of literature data and tested several formulas, including the Hilsdorf one, but many of the lacking input data were simulated by analogy. Boffill [44] considered historic brickworks with low strength mortar. Even some multiple leaf wallets were tested and discussed [45,46]. To point out the importance of the data entering in the Hilsdorf formula, the analysis carried out by Drougkas [43] has been considered again, by deriving the missing data with the formulas proposed in this study (Online Source 1) and plotting the computed strength versus the experimental one (Fig. 9). In the cases of Running, Flemish and English bond, the ratio k = 2 of two mortar layers for every block has been set out according to Aprile et al. [9], due to the presence of the head joints.
The new analysis is worthy of only 34% average error (coefficient of variation 1.27), while the numerical analysis of Drougkas [43] is as high as 43% average error (coefficient of variation of 1.77). Both the analyses tend to overestimate (positive error) the experimental strengths of the masonry stacks, but, except for a few samples with very uncommon material combinations, the use of the proposed formulas could certainly replace any complex FE analysis with reasonable accuracy. In particular, from Fig. 9 it is apparent that the proposed formula does overestimate the resistance of weak stacks, while it is certainly on the safe side for strong stacks and wallets. The reason for this overestimation can be found in the premature collapse by crushing of the mortar layers, which is not considered in the proposed formula.

The failure of the mortar courses
In respect of the previous discussion, the mortar plasticization and crushing in biaxial compression rarely happens in the collapse of the masonry structure. This is consistent with the experiments of Lumantarna [47], in which several combinations of bricks and mortars have been tested on-site and in reconstructed specimens in the laboratory. The average strength of masonry stacks tested is only 65% of the uniaxial strength of blocks and 243% of the uniaxial mortar strength. A rough analytical justification can be obtained by examining the multi-axial compression failure of the mortar in light of the Coulomb criterion [48,49].
By using formula (7) and recognizing that: The following result for the failure compression is obtained: In general, U is large when R is small, and the reverse. In any case, the mortar will have a large compressive resistance due to confinement, since the denominator is small.

A new formula for brick blends
The starting point for the HET formulation is again the compatibility and equilibrium equations of the stack. The solution is developed in a plane strain condition, where mortar and blocks are defined by their thicknesses t m and t b , their elastic moduli E m and E b,i , and their Poisson's ratios m m and m b,i , respectively, with a given combination of volumetric fractions q i . In particular, concerning a general masonry wall containing randomly blended units, the volumetric fraction of block type i is computed by considering the number of blocks n b,i in the panel with their corresponding thicknesses t b,i over the whole composition of the wall.
The equilibrium equation in the horizontal and vertical directions are easily written referring to the mortar and block courses present in the wall: A non-dimensional formula is obtained by dividing the horizontal stress by the vertical stress, and each term by the total thickness of the block courses: where s represents the ratio r h /r v , and kg is the ratio of the mortar thickness to the total block thickness. In practice, k represents the ratio of the number of mortar courses over the number of blocks, while g is the ratio of the mortar course average thickness, over the block average thickness.
The horizontal strain of all the materials is defined by the following formula, in which compression and shortening were considered positive: The compatibility condition states that the horizontal strain of the mortar and the linked bricks must be the same. By denoting the constrained modulus as D, and introducing the non-dimensional stresses s hb,i and s hm , the compatibility equations hold: The Eq. (16) encodes the compatibility equation of the horizontal flexibilities e h of the two materials, and can be solved by extracting s hb,i as a function of s hm . By substituting all the block nondimensional stresses in Eq. (13), the stress ratio s hm in mortar is computed: Since mortar is considered not involved in the failure process, it is necessary to introduce hypotheses capable to split the horizontal stress resultant of the mortar in the horizontal stresses acting on the different blocks composing the blend. Failure can thus be set either at the lower vertical stress, by producing collapse in the weakest block type, or with a sort of average stress which will bring all the blocks to the limit cracked state at the same time. The first situation can be thought of as a limit elastic brittle one, while the second is a rigid plastic condition. In what follows the two solutions are presented and discussed.

Limit elastic brittle condition for a brick blend
In this case, the block stresses are derived from the mortar stresses by using the relative stiffness H b,i in the horizontal direction of the blocks composing the blend. Since: The mortar stress ratio s hm is split using the different block type stiffness ratios K b,i computed in terms of the material rigidities H b,i of the blocks and their volumetric fractions q i : Then, the failure stress f M,EL is computed based on the Coulomb formula (3) for all the present block types, and the minimum strength is considered the overall failure stress: Thus, for a stack of blended blocks, the forecast elastic stress ratios are s hm for the mortar, and s hm k g K b,i for the different block types. The effectiveness of these formulas can be checked against the results of the FEM simulations presented above. In the following Fig. 10 the main results for the CL and CS combinations are presented, by referring to the average of interspersed and clustered patterns.
As is evident from the figures, the stress ratios are in very good agreement over the whole range of volumetric ratios of the studied blends. Indeed, the coefficient of determination R 2 corresponding to the prediction capability of the proposed formulas shows values higher than 0.95 for quite all the numerical investigations. However, lower values are found dealing with the stresses associated either with limestone blocks or thicker clay bricks in the bricksandstone blends.

Limit rigid plastic condition for a brick blend
If the collapse is postponed till cracks span all over the vertical masonry length, the block blend components are contributing to the final strength based on their respective volumetric ratios, as a sort of strength average. Therefore, the combination formula is obtained from the strength of the various local single block type strengths combined with their shares, so that the wallet strength becomes equal to each component's strength if the volumetric ratio of this component is unity: As will be evident in the comparison with the experiments, the two hypotheses lead to similar results, if the ductility of the system due to the heterogeneity is limited.

Comparison of HET with experimental results
The proposed formulas were used in the interpretation of the presented experiments. In addition to these formulas, since in the literature the fractional power   formula proposed by EC6 [5] is widely used in checking results (see for example Ferretti [4]), an extended EC6 formula never proposed before in the literature, is derived for the case of blended blocks: Suitable values for the parameters are k = 0.6 7 0.7, a = 0.25 7 0.30, b = 0.75 7 0.70.
The formula (22) is based on the consideration that in a stack of random elastic perfectly plastic blocks which yield approximately at the same strain, the strength can be computed as the average modulus multiplied by the yielding strain, and therefore by the weighted harmonic mean of the strengths of the various components: (Online Source 2) shows the comparison with the experiments, while Fig. 11 collects all the strength value pairs comparing the performances of the analytical formulas against the experimental outcomes. Although the experiments are very complex, with materials with very different and scattered properties, the agreement of the proposed formulas is considerable, and even the very simple Euro Code 6 power formula proves to be a good evaluation tool. However, it is to mention that a correct evaluation of the data like f t , /, m is of paramount importance in deriving the mechanical parameters to be used in the proposed formulas. Obviously, when only one type of block is present, plastic and elastic masonry strengths are exactly equal.

Conclusions
In this study, a systematic procedure, capable of evaluating the compressive strength of masonry wallets with blended blocks, has been presented and checked against available experimental investigations.
Although the procedure is considering the collapse as a limit elastic process ending with the onset of a diffuse cracking in the block blend forming the masonry, the result is a precise strength evaluation, if the mechanical parameters of the materials composing the walls are defined with care.
The mortar-unit interaction problem typical of Haller Hilsdorf theory has been solved as a plane strain equilibrium problem based on a wide numerical investigation.
The importance of material parameters often lacking in experimental studies, namely the tensile strength and the Poisson ratio, has been stressed. Some useful formulas able to fill the gap, allowing the estimation of these values from the uniaxial compressive strength, were reviewed.
From the analysis, it is clear that the heterogeneity due to a dispersed distribution of blocks showing different properties is producing a sort of apparent ductility which comes out in a harmonic averaging of the block strengths. Moreover, the horizontal compressive stress in the mortar layers is mostly depending on the mortar Poisson ratio only, while the horizontal tensile stress in the blocks is even a function of the stiffness ratio of each block type in turn.
A new analytical tool, namely the Hilsdorf Extended Theory (HET) was derived in both the elastic and plastic ranges. A wide experimental investigation on 1:4 scale specimens carried out at the University of Firenze is discussed in detail.
The comparison of the theoretical results with the outcomes of the experimental investigation and the corresponding R 2 values point out the good predictive capability of both the derived analytical formulas.
A further future step of the theoretical derivation will concern the introduction of the effect of the bonding pattern and the analysis of blended block walls with multiple leaves. It is finally to cite that one of the most important and disregarded components of the real walls is the presence of voids, which is rarely investigated although it constitutes a phase with a significant strengthless volume fraction, introducing in this way a consistent capacity reduction.
Acknowledgements The authors want to acknowledge the deep influence of prof. Silvia Briccoli Bati in diffusing the culture of experimental testing for masonry structures. Acting as the Tutor of Daniela Fenocchi and Sara Rosi, she designed and directed more than 20 years ago the experiments, which are used now in checking the new Extended Hilsdorf Theory.
Funding This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors.

Declarations
Conflict of interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.