Next Article in Journal
Model Driven Interoperability for System Engineering
Previous Article in Journal
Modelling—International Open Access Journal of Modelling in Engineering Science
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Tensile Fracture Mechanism of Masonry Wallettes Parallel to Bed Joints: A Stochastic Discontinuum Analysis

1
Architectural Engineering, University of Nebraska-Lincoln, Lincoln, NE 68182, USA
2
Civil Engineering, Bogazici University, Istanbul 34342, Turkey
3
ISISE, Institute of Science and Innovation for Bio-Sustainability, Department of Civil Engineering, University of Minho, 4800-058 Guimarães, Portugal
4
Laboratório Nacional de Engenharia Civil (LNEC), 1700-075 Lisboa, Portugal
5
Itasca Consulting Group, Minneapolis, MN 55401, USA
*
Author to whom correspondence should be addressed.
Submission received: 2 August 2020 / Revised: 14 September 2020 / Accepted: 2 October 2020 / Published: 8 October 2020

Abstract

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

1. Introduction

Masonry is a heterogeneous, nonlinear, and composite construction material consisting of units and mortar. The components within a masonry assemblage may have a considerable difference in terms of stiffness and strength, which makes it challenging to predict the deformation behavior (e.g., stress–displacement curve) of the material at the macro-level. Additionally, the heterogeneous character of mortar, unit, and unit–mortar interface properties inherently leads to large scattering in the experimental results, as noted by several researchers [1,2]. From the computational point of view, most of the numerical models presented in the literature are deterministic and rely on the continuum-based representation of the masonry, described as an equivalent orthotropic continuum associated with either plasticity or damage constitutive laws. In the present research, an alternative approach is employed using the discrete element method (DEM), where the masonry texture is replicated via discrete polyhedral blocks, as shown in Figure 1. Note that generated blocks are grouped to simulate masonry units and unit–mortar interfaces, where each masonry unit in the discrete model is expanded up to the half-thickness of the mortar joints (see Figure 1). Therefore, with the help of the proposed discontinuum-based modeling strategy, cracking in the masonry units and sliding/tension failure modes at the unit–mortar interfaces can be captured. The mechanical interaction among the blocks is represented using nonlinear spring models at the contact points, defined among the adjacent blocks, similarly to work published elsewhere [3,4,5]. Hence, utilizing the discontinuum representation of masonry, this research aimed to better understand the influence of material properties in masonry wallettes on the fracture mechanism and stress–displacement behavior, subjected to tension parallel to bed joints.
In a preceding study by the authors, Pulatsu et al. [6], a deterministic modeling strategy was used to explore the tensile behavior of the masonry wallettes using average contact properties. The contribution of the present study furthers this investigation and explores the effect of statistically varied contact properties. In other words, the emphasis of this work is on the stochastic analysis of masonry wallettes within the discrete element method framework. Therefore, the goals of this research are as follows:
  • To investigate the effect of nonlinear contact properties on the macro behavior of the masonry wallettes subjected to tension parallel to bed joints.
  • To gain a better understanding of the tensile fracture mechanism in masonry, considering strong masonry unit–weak mortar joint and weak masonry unit–strong mortar joint systems.
In the following sections, a brief overview of the applied computational method and the statistical sampling method are presented.

2. Computational Framework

The employed modeling strategy, discrete element method (DEM), is developed by Cundall [7] in the early 1970s to analyze the progressive collapse mechanism of the jointed rock masses. Since then, DEM has been applied to various engineering disciplines ranging from micro- to macro-scale. Fundamentally, DEM was developed to explore the mechanical behavior of the discontinuous systems by considering them as an assembly of rigid and/or deformable blocks, interacting with each other through the contact points. Deformability of the blocks is introduced by constant strain tetrahedral element (CST) discretization of the block, referred to as finite-difference zones, that may be considered as an internal mesh.
In DEM, the mechanical interaction between the adjacent polyhedral blocks is simulated based on the point contact hypothesis [8], as shown in Figure 2. The random generation of the polyhedral blocks is performed using an open-source software package NEPER. The interested readers can access the mathematical background and applications of this software in related sources [9,10,11,12,13]. The elastic response of the contact is controlled by the orthogonal springs both in the normal and shear directions, based on the assigned contact stiffness k n and k s , respectively; whereas the post-peak behavior is governed by the nonlinear parameters, namely tensile strength ( f T ) defined in the normal direction and cohesion ( c ) as well as the friction angle ( ϕ ) in the shear direction (see Figure 2). In the present study, semi-rigid blocks are used considering relatively high stiffness at the finite-difference zones (i.e., 10 times macro elastic stiffness) to lump the linear and nonlinear displacements at the joints as performed in relevant studies [3,6,13]. Therefore, deformations develop at the joints (i.e., along the boundaries of adjacent polyhedral blocks) instead of within the block domain and block deformations have negligible influence on the macro stress–displacement behavior of the system.
The computational procedure of DEM relies on the explicit integration scheme of the equations of motion for the nodal points. The governing equations of motion are solved using the central difference method and turn into a compact form (see Equation (1)), evaluated at the mid-time intervals ( t + = t + Δ t / 2 , t = t Δ t / 2 ; Δ t : time step).
u ˙ i t + = u ˙ i t + [ Σ F i t ( F d ) i ] Δ t m ( F d ) i = λ | F i t | s g n ( u ˙ i t )
where m , u ˙ , and F are the nodal mass, nodal velocity vector, and the resultant nodal force vector, respectively. The resultant nodal force vector is calculated as the sum of the external forces, contact forces, gravity forces, and the contribution to the internal stress in the zones adjacent to the nodes. Furthermore, F d and λ denote the damping force and dimensionless damping constant (default value is 0.8). Note that local type of damping is applied during the analyses, which is proportional to the magnitude of the unbalanced force ( | F i t | ) and applied to the opposite direction of the velocity vector [14]. Thus, the damping force varies from one point to another in the model, depending on the magnitude of the unbalanced force. Once the nodal velocities are computed, nodal displacements ( Δ u i = u ˙ i t + Δ t ) are obtained to update the position of the blocks. Accordingly, the contact stresses, developing among the adjacent blocks, are calculated based on the defined contact stress–displacement laws then multiplied with the contact area to be utilized in Equation (1). It is worth noting that small displacement theory is assumed through the computational procedure; thus, only initially established contacts are taken into consideration and updated during the analyses. At each time step, elastic stress increments ( Δ σ ,   Δ τ ) are computed using the relative contact displacements ( Δ u n ,   Δ u s ) and added to the previous contact stresses ( σ o l d ,   τ o l d ) to obtain the new ones ( σ n e w ,   τ n e w ) as expressed in Equation (2).
σ n e w = σ o l d + Δ σ ;   Δ σ = k n Δ u n τ n e w = τ o l d + Δ σ ;   Δ τ = k s Δ u s
Moreover, the dilatancy is taken into account upon the shear failure at the contact point, considering the dilation angle (ψ) to modify the normal stress increment as follows; Δ σ T o t a l = Δ σ + k n Δ u s t a n ψ . The numerical stability of the explicit integration scheme is ensured by utilizing a sufficiently small time-step during the analysis, calculated as follows.
Δ t = 2 m / k g p
where k g p denotes nodal stiffness obtained by adding zone and contact stiffness [15]. Hence, the computational procedure is executed in a cyclic manner, presented in Figure 3, until the quasi-static equilibrium is reached at the discontinuum system.
Throughout this research, a commercial discrete element code 3DEC, developed by Itasca, is used, where the custom contact constitutive models are adopted into the software (3DEC) via a user-defined constitutive model option [16].

2.1. Contact Constitutive Models

As mentioned earlier in the article, the linear response (elastic deformation) of the discrete models is controlled by the contact stiffnesses (normal and shear springs), which can be computed by Equation (4), given in [17,18].
k n = E / t ; k s = G / t
where t indicates the average thickness of the fracture zone, whereas E and G denote the macro elastic and shear modulus of masonry, respectively. Contact stiffness values have an inverse relationship with the thickness of the fracture zone (or the number of blocks), which means that the contact stiffness increases with a higher number of blocks, as discussed in the literature [13]. In the nonlinear regime, computed contact stresses are updated based on the defined failure criteria. Typically, most discrete element models consider sudden stress drops at the contact points in the normal direction (i.e., zero tensile strength is adopted, or the contact is deleted after the contact stress reaches the ultimate stress) and residual shear strength upon the sliding failure. However, softening behavior (i.e., a gradual decrease of stress with increasing deformation) is a salient feature of the quasi-brittle materials (e.g., concrete, mortar, clay brick, and rock) due to a process of progressive internal crack growth [19,20]. To address this gradual degradation and the corresponding decrease in the load carrying capacity of the material, nonlinear contact constitutive models are implemented in 3DEC to be utilized in tension and shear failure criteria considering exponential softening functions (Figure 4). The mathematical formulations of the exponential functions, given by Lourenço et al. [21], and further details about the contact models can be found in the related study [6]. It is important to note that the Coulomb slip joint model is utilized in the shear direction, where the maximum and residual shear strengths are calculated as written in Equation (5), and the descending branch beyond the peak is defined via softening of the cohesion [22]. Once the contact stress reaches its capacity, the implemented contact models are activated and updated (new) contact stresses are calculated. The contact models (linear and polynomial softening) developed by the authors are available online and free to download for all 3DEC users [23].
τ m a x = c 0 + t a n ϕ 0   σ ;   τ r e s = c r e s + t a n ϕ r e s   σ

2.2. Benchmark Study and Testing Setup

The cracking mechanism of masonry depends on the mechanical properties of its constituents and their interactions. However, most experimental studies are related to the tensile response of masonry components (i.e., bricks) and tension perpendicular to the unit–mortar interface [24,25]. There are a limited number of works available in the literature related to the tensile response of masonry parallel to bed joints [26].
In this research, the well-known test setup and experimental results, presented by Backes [2], are used as the benchmark study, where different types of masonry units and mortar combinations were tested. Throughout the testing program, masonry wallettes, composed of four courses of units with unit dimensions of 240 mm × 115 mm × 113 mm, were subjected to uniaxial tensile forces applied through steel plates that were glued to both sides of the wallette. In Figure 5, an illustration of the test setup and boundary conditions are demonstrated on the discontinuum model, in which the lateral deflections are imposed to the right-hand side of the specimen while the opposite side is restricted. Specifically, the constant nodal velocities are prescribed for all gridpoints at the right-hand side of the specimen, whereas gridpoints at the opposite side are restricted in the horizontal direction. To ensure the numerical stability during the analysis, a relatively low displacement rate (i.e., 1 mm/s) is defined at the gridpoints, allowing to obtain a smooth response from the numerical solutions as discussed in [27]. Reaction forces are extracted from gridpoints, in which the constant velocity is applied, via the implemented subroutine in the software based on FISH functions (an executable programming language in 3DEC).
The experimental findings show that two failure mechanisms develop depending on the type of masonry constituents: (i) tensile and shear debonding at the unit–mortar interfaces (failure mode I) and (ii) tensile failure in the mortar joint and cracking at the units (failure mode II). The former was obtained when the weak mortar joint with strong masonry units was used, whereas the latter was observed when weaker masonry units and relatively strong mortar joint were utilized. Typical failure patterns obtained from the experiment are shown in Figure 6a. Furthermore, these two failures modes led to different stress–displacement (or stress–strain) curves, where the failure mode I resulted in an approximately constant residual stress after the peak due to interlocking effect of the units and friction, and the failure mode II demonstrated a sharp descending branch after the maximum tensile stress was obtained (see Figure 6b), as typical in quasi-brittle materials such as concrete. For further details about the testing, readers are referred to [2].

3. Statistical Representation of the Material Properties

The elastic modulus of the masonry assembly ( E ), the tensile strength of the masonry unit ( f t , u ), the bond (unit–mortar interface) tensile strength ( f t , i ), bond cohesion ( c ), the friction angle ( ϕ ), and the dilation angle ( ψ ) are the random variables considered in the probabilistic analyses. The mean values of variables are either obtained from the experimental data (if applicable) [2] or adopted from similar tests, presented in the literature [22,28,29]. Furthermore, the residual friction angle ( ϕ r e s ) is taken as 20% less than initial friction angle ( ϕ ) for all the joints.
The normal distribution is determined for the contact friction angle and the elastic modulus of masonry following the examples in literature [30,31,32,33], while lognormal distribution is assumed for the other parameters. The distributions of the random variables and their parameters are presented in Table 1 and Table 2 for the strong unit–weak bond (SU-WB) and weak unit–strong bond (WU-SB) systems, respectively. Note that the shear contact parameters for the masonry units (i.e., cohesion ( c u ) and friction angle ( ϕ u )) are kept constant throughout the analyses, where c u = 1.5 f t , u and ϕ u = 45 ° , respectively. Moreover, mode-I and -II fracture energies are computed based on the tensile strength and cohesion parameters, as suggested in [34,35].
Dependency between the random variables is considered via the correlation coefficient, which is a linear dependency ratio. A strong relationship among f t , i , c , and ϕ is assumed and defined by 85% correlation, whereas a moderate relationship between E and f t , u is represented by a 50% correlation. With the predefined statistical distributions and the correlation between the random variables, the Latin Hypercube Sampling (LHS) method is used to derive sample values for the simulations. The LHS method reduces the variance in the Monte Carlo Simulation and, thus, lowers the number of required simulations [36,37]. For each case (SU-WB and WU-SB), 200 simulations—400 simulations in total—were conducted to derive the random parameters. In Figure 7, the graphical illustration of the input variables for WB-SM, together with the corresponding histograms and probability distributions, is presented.
The number of simulations is determined by monitoring the relative difference between the mean values after a number of simulations with the prescribed mean values in Table 1 and Table 2. The percentage of the error is used to express the relative difference. The magnitude of the difference in the mean and standard deviation values after “j” simulations can be tracked in Figure 8 and Figure 9. As can be seen from these figures, the final values of the errors/or differences are smaller than 0.1% for most of the random variables.

4. Results of the Computational Models

In this section, the results of the computational models are presented. A range of values for input parameters ( E , f t , u , f t , i , c , and ϕ ) are established for each simulation described in the previous section, using statistical methods. The results of the analyses clearly show that the macro behavior of the masonry wallettes reveals two separate stress–displacement trends corresponding to SU-WB and WU-SB combinations, presented in Figure 10. As mentioned earlier in the article, similar stress–displacement behaviors were also noted from the experimental observations. Moreover, readers may find a brief discussion regarding the effect of block size and shape provided in the Appendix A.

4.1. Failure Mode I: Strong Unit–Weak Bond (SU-WB) Behavior

When the strength of the units is higher than the mortar joints (i.e., strong unit–weak bond system), experimental work demonstrated that the failure mode is governed by the unit–mortar interfaces. A typical failure (denoted as mode I) indicates a zig-zag failure pattern, including shear deformations (sliding) at the bed joints and the cracks at the head joints. A similar failure mechanism is successfully obtained using the proposed modeling strategy, shown in Figure 11a. Furthermore, performed stochastic analyses provide a wide range of stress–displacement results depending on the defined contact properties (see Figure 11b). Regarding the fracture pattern, no significant difference is observed during the analysis for this particular set of material properties. Results of the analyses show that there is a 94.8% correlation between the bond strength ( f t , i   and   c ) and the overall capacity, and an 81.4% correlation between unit–mortar interface (bond) friction angle ( ϕ ) and the overall capacity (Figure 12). Since the governing mechanism for mode I is concentrated at the bed and head joints, the capacity is strongly affected by bond tensile strength, cohesion, and friction angle, as expected. The strong linear correlation between the tensile capacity and these two modeling parameters even allows capacity prediction without requiring any simulations.

4.2. Failure Mode II: Weak Unit–Strong Bond (WU-SB) Behavior

The second failure mode of masonry wallettes includes more complex fracture patterns compared to mode I. Specifically, when the tensile strength of the masonry units gets close to the bond strength or even smaller than the bond strength, several distinctive fracture patterns are observed via the proposed modeling strategy, categorized as two subgroups (fracture mechanisms) of failure mode II. In the first subgroup, cracks develop through the head joints and units (FM-1 and FM-4 in Figure 13), which was given as a typical fracture pattern in the benchmark study. On the other hand, in the second subgroup, only the units fail under tension, corresponding to the condition where the bond tensile strength is higher than the tensile strength of the masonry units ( f t , i > f t , u ), denoted as FM-2 and FM-3 in Figure 13. In addition, a typical sliding mechanism (same as the mode I) is obtained from the discontinuum models, reflecting the SU-WB combination (Figure 13c). However, this should be expected depending on the determined data set and statistical distribution, in which there is a chance to assign contact properties, where the masonry units may have relatively higher tensile strength than the bond ( f t , i f t , u ). Fracture simulations of FM-1 and FM-2 are shared as video at the supplementary documents.
In the benchmark study, FM-1 was presented as the typical fracture pattern for strong bond masonry wallettes, mentioned in Section 2.1. Our results also show that FM-1 has the highest probability with 43%. In order to compare the stress–displacement curves obtained for FM-1 from the discontinuum analyses, with the experimental result, they are plotted in the same graph in Figure 14. This comparison demonstrates that the determined response envelope for the range of input material properties in DEM is in good agreement with the experimental findings and can be considered representative of possible dispersion of a larger number of testing series. For instance, in Figure 14, the average trend of the numerical solutions resembles the experimental one, although an exact match is not possible due to assumed statistical distributions of the random parameters.
Furthermore, the influence of contact parameters on the system’s overall capacity is investigated considering each failure mode. The results show that there is a 64% correlation between the unit tensile strength and capacity for WU-SB masonry wallettes, where the dispersion of the relation regarding each specific failure mechanisms can be found in Figure 15a. FM-2 and FM-3 indicate a perfect linear correlation (100%) between the tensile strength of the unit and the capacity, while the scattered distribution of the data is related to the mechanisms involving the bond failure (Figure 15a).
A similar correlation factor, 65%, is obtained for the bond tensile strength ( f t , i ) and the ultimate load-carrying capacity. Note that, for bond tensile strength, FM-1 and FM-5 present 92% and 99% correlation, whereas a more scattered data cloud is noted for FM-2 and FM-3 since these failure modes are directly related to the cracking in masonry units (see Figure 15b). Finally, the lowest correlation factor, 54%, is computed for the contact friction angle at the unit–mortar interfaces. Given the fact that only the FM-5 reveals a sliding mechanism at the bed joints, there should be a low correlation expected for the other fracture mechanisms (Figure 15c). Partitioning of the modeling parameters for each fracture mechanism allows deriving further correlation with the load-carrying capacity and corresponding failure.

5. Conclusions

In this research, a stochastic discontinuum analysis is performed within the discrete element method framework to analyze the diverse fracture mechanism of masonry wallettes subjected to tension parallel to the bed joints. According to the results of the analyses, the following conclusions are derived.
  • The proposed modeling strategy allows detailed and realistic (based on comparison with experimental results) fracture patterns of masonry wallettes. This is accomplished by representing the masonry components (masonry units and unit–mortar interfaces) using mechanically interacting discrete polyhedral blocks.
  • The inherent uncertainty in the mechanical properties of masonry constituents yields various possible fracture mechanisms. Through this research, all possible fracture patterns of masonry wallettes under tension parallel to bed joints are explored. There is a great economy in performing such a comprehensive study computationally instead of a very broad and costly experimental campaign.
  • Stochastic analyses enable us to attain the dispersion of the load-carrying capacity concerning the modeling parameters ( E , f t , u , f t , i , c , and ϕ ). They highlight both the direct and indirect influence of these parameters on the fracture patterns and stress–displacement behavior. Direct correlation of the parameters such as the bond tensile strength, friction angle, or the unit tensile strength with the capacity reflects the physically expected phenomena. In the future, results of these determined correlations can be used to estimate the capacity of different arrangements of masonry materials without using numerical simulations and may lead to the derivation of new empirical formulas to predict the tensile capacity of masonry parallel to bed joints.
  • For the given probability distributions of the material parameters, the fracture mechanisms and their likelihood have been found. Moreover, the influence of the modeling parameters on the tensile strength of masonry has been quantified. However, it should be noted that the results presented within this study may be restricted to the statistical distributions and the parameters used.
Currently, the practical applications of this work are restricted to the capacity prediction concerning the change in the modeling parameters for specific fracture mechanisms. Future studies (both experimental testing and stochastic numerical analyses) considering different material properties and masonry unit sizes will help define reference values and empirical equations for tensile strength prediction, which can then be used directly by practitioners. A further step ahead would be using statistical learning methods to predict both the fracture mechanism and the tensile strength in the multi-dimensional parameter space.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2673-3951/1/2/6/s1, Video S1: SU-WM-failure mode I, Video S2: WU-SB (FM-1), Video S3: WU-SB (FM-2).

Author Contributions

Conceptualization, original draft preparation, software, validation, and editing, B.P., S.G.; methodology, supervision, writing—review and editing, E.E.; methodology, supervision, review and editing, P.B.L.; software, review and editing, J.V.L.; software, review and editing, J.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

In this section, the influence of block size and shape are presented regarding capacity and stress–displacement behavior. The former was discussed in a preceding study published by the authors, where no severe block size dependency was observed from the results of discrete element models (see [6]). Note that discontinuum models were grouped according to the size of the polyhedral blocks, denoted as coarse, medium (present study), and fine, that were analyzed using one set of contact properties, presented in [6]. The results are given in Figure A1.
Figure A1. Influence of block-size on the capacity and stress–displacement response, Left: failure mode I; Right: failure mode 2 (taken from [6]).
Figure A1. Influence of block-size on the capacity and stress–displacement response, Left: failure mode I; Right: failure mode 2 (taken from [6]).
Modelling 01 00006 g0a1
On the other hand, the latter (i.e., the effect of polygonal structure (or morphology)) requires more research and only briefly discussed in this section based on the preliminary results for WU-SB. Here, masonry composition is represented by different morphological features, as shown in Figure A1. Each morphology is generated via the NEPER software package, mentioned earlier in the manuscript [9,10,11]. Again, one set of contact properties are utilized that is taken from a preceding study [6]. Although preliminary results do not indicate significant differences in terms of capacity and demonstrate the same fracture mechanism (see Figure A2), more research is needed, especially for different loading conditions and contact properties.
Figure A2. Influence of morphology on the stress–displacement and fracture mechanism (deterministic discontinuum analysis).
Figure A2. Influence of morphology on the stress–displacement and fracture mechanism (deterministic discontinuum analysis).
Modelling 01 00006 g0a2

References

  1. Lourenço, P.B.; Almeida, J.C.; Barros, J.A. Experimental investigation of bricks under uniaxial tensile testing. Mason. Int. 2005, 18, 11–20. [Google Scholar]
  2. Backes, H.P. Tensile strength of masonry. In Proceedings of the 7th International Brick Masonry Conference, Melbourne, Australia, 17–20 February 1985; pp. 779–790. [Google Scholar]
  3. Pulatsu, B.; Erdogmus, E.; Lourenço, P.B.; Quey, R. Simulation of uniaxial tensile behavior of quasi-brittle materials using softening contact models in DEM. Int. J. Fract. 2019, 217, 105–125. [Google Scholar] [CrossRef]
  4. Resende, R.; Lemos, J.V.; Dinis, P.B. Application of a discontinuity model with softening to the analysis of dam foundations using the Discrete Element Method. In Proceedings of the 1st International UDEC/3DEC Symposium: Numerical Modelling of Discrete Materials in Geotechnical Engineering, Civil Engineering and Earth Sciences, Bochum, Germany, 29 September–1 October 2004; pp. 249–255. [Google Scholar]
  5. Sarhosis, V.; Lemos, J.V. A detailed micro-modelling approach for the structural analysis of masonry assemblages. Comput. Struct. 2018, 206, 66–81. [Google Scholar] [CrossRef]
  6. Pulatsu, B.; Erdogmus, E.; Lourenço, P.B.; Lemos, J.V.; Hazzard, J. Discontinuum analysis of the fracture mechanism in masonry prisms and wallettes via discrete element method. Meccanica 2020, 55, 505–523. [Google Scholar] [CrossRef]
  7. Cundall, P.A. A computer model for simulating progressive, large-scale movements in blocky rock systems. In Proceedings of the The International Symposium on Rock Mechanics, Nancy, France, 4–6 October 1971; pp. 47–65. [Google Scholar]
  8. Lemos, J.V. The Basis for Masonry Analysis with UDEC and 3DEC. In Computational Modeling of Masonry Structures Using the Discrete Element Method; ACIE: Hershey PA, USA, 2016; pp. 61–89. [Google Scholar]
  9. Quey, R. Neper Reference Manual. Available online: https://neper.info/docs/neper.pdf (accessed on 17 September 2020).
  10. Quey, R.; Dawson, P.R.; Barbe, F. Large-scale 3D random polycrystals for the finite element method: Generation, meshing and remeshing. Comput. Methods Appl. Mech. Eng. 2011, 200, 1729–1745. [Google Scholar] [CrossRef] [Green Version]
  11. Quey, R.; Renversade, L. Optimal polyhedral description of 3D polycrystals: Method and application to statistical and synchrotron X-ray diffraction data. Comput. Methods Appl. Mech. Eng. 2018, 330, 308–333. [Google Scholar] [CrossRef] [Green Version]
  12. Pulatsu, B.; Kim, S.; Erdogmus, E.; Lourenço, P.B. Advanced Analysis of Masonry Retaining Walls Using Mixed Discrete-Continuum Approach. Available online: https://0-www-icevirtuallibrary-com.brum.beds.ac.uk/doi/abs/10.1680/jgeen.19.00225 (accessed on 2 August 2020).
  13. Pulatsu, B.; Erdogmus, E.; Lourenço, P.B.; Lemos, J.V.; Tuncay, K. Numerical modeling of the tension stiffening in reinforced concrete members via discontinuum models. Comput. Part. Mech. 2020. [Google Scholar] [CrossRef]
  14. Cundall, P.A. Distinct element models of rock and soil structure. In Analytical and Computational Methods in Engineering Rock Mechanics; Brown, E.T., Ed.; George Allen Unwin: London, UK, 1987; pp. 129–163. [Google Scholar]
  15. Lemos, J.V. Block modelling of rock masses. Concepts and application to dam foundations. Rev. Eur. génie Civ. 2008, 12, 915–949. [Google Scholar] [CrossRef]
  16. Itasca Consulting Group Inc. 3DEC Three Dimensional Distinct Element Code. Available online: https://www.itascacg.com/software/3dec (accessed on 2 August 2020).
  17. Kazerani, T. A discontinuum-based model to simulate compressive and tensile failure in sedimentary rock. J. Rock Mech. Geotech. Eng. 2013, 5, 378–388. [Google Scholar] [CrossRef] [Green Version]
  18. Kazerani, T.; Zhao, J. Micromechanical parameters in bonded particle method for modelling of brittle material failure. Int. J. Numer. Anal. Methods Geomech. 2010, 34, 1877–1895. [Google Scholar] [CrossRef]
  19. Lourenço, P.B. Computational Strategies for Masonry Structures; Delft University of Technology: Delft, The Netherlands, 1996. [Google Scholar]
  20. Van Mier, J.G.M.; Man, H.-K. Some Notes on Microcracking, Softening, Localization, and Size Effects. Int. J. Damage Mech. 2009, 18, 283–309. [Google Scholar] [CrossRef]
  21. Lourenço, P.B.; Rots, J.G. Multisurface interface model for analysis of masonry structures. J. Eng. Mech. 1997, 123, 660–668. [Google Scholar] [CrossRef]
  22. Van der Pluijm, R.; Rutten, H.; Ceelen, M. Shear behaviour of bed joints. In Proceedings of the 12th International Brick/Block Masonry Conference, Madrid, Spain, 25–28 June 2000; pp. 1–14. [Google Scholar]
  23. Itasca Consulting Group Inc. UDM Library-Fracture Energy Softening Model. Available online: https://www.itascacg.com/software/udm-library/fracture-energy-softening-model (accessed on 7 September 2020).
  24. van der Pluijm, R. Out-of-Plane Bending of Masonry Behaviour and Strength; Technische Universiteit Eindhoven: Eindhoven, The Netherlands, 1999. [Google Scholar]
  25. Lourenço, P.B.; Barros, J.A.; Almeida, J.C. Characterization of masonry under unixial tension-Report 02-DEC/E-12; Universidade do Minho: Guimarães, Portugal, 2002. [Google Scholar]
  26. Backes, H.P. On the Behavior of Masonry under Tension in the Direction of the Bed Joints. Ph.D. Thesis, Aachen University of Technology, Aachen, Germany, 1985. [Google Scholar]
  27. Pulatsu, B.; Erdogmus, E.; Lourenço, P.B.; Lemos, J.V.; Tuncay, K. Simulation of the in-plane structural behavior of unreinforced masonry walls and buildings using DEM. Structures 2020, 27, 2274–2287. [Google Scholar] [CrossRef]
  28. Atkinson, R.H.; Amadei, B.P.; Saeb, S.; Sture, S. Response of masonry bed joints in direct shear. J. Struct. Eng. 1989, 115, 2276–2296. [Google Scholar] [CrossRef]
  29. Schubert, P. Tensile and flexural strength of masonry - influences, test methods, test results. In Proceedings of the 10th International Brick/Block Masonry Conference, Calgary, AB, Canada, 5–7 July 1994; pp. 895–908. [Google Scholar]
  30. Parisi, F.; Augenti, N. Uncertainty in seismic capacity of masonry buildings. Buildings 2012, 2, 218–230. [Google Scholar] [CrossRef] [Green Version]
  31. Conde, B.; Matos, J.C.; Oliveira, D.V.; Riveiro, B. Probabilistic-based structural assessment of a historic stone arch bridge. Struct. Infrastruct. Eng. 2020, 1–13. [Google Scholar] [CrossRef]
  32. Pina-Henriques, J.; Lourenço, P.B. Masonry compression: A numerical investigation at the meso-level. Eng. Comput. (Swansea, Wales) 2006, 23, 382–407. [Google Scholar] [CrossRef] [Green Version]
  33. Moreira, V.N.; Fernandes, J.; Matos, J.C.; Oliveira, D.V. Reliability-based assessment of existing masonry arch railway bridges. Constr. Build. Mater. 2016, 115, 544–554. [Google Scholar] [CrossRef]
  34. Ghiassi, B.; Vermeltfoort, A.T.; Lourenço, P.B. Masonry mechanical properties. In Numerical Modeling of Masonry and Historical Structures; Ghiassi, B., Milani, G., Eds.; Elsevier: Amsterdam, The Netherlands, 2019; pp. 237–261. ISBN 9780081024393. [Google Scholar]
  35. Lourenço, P.B. Recent advances in masonry structures: Micromodelling and homogenisation. Available online: https://www.worldscientific.com/doi/abs/10.1142/9781848163089_0006 (accessed on 2 August 2020).
  36. Boussabaine, A.; Kirkham, R. Whole Life-Cycle Costing: Risk and Risk Responses; Wiley: Oxford, UK, 2004. [Google Scholar]
  37. Lemaire, M.; Chateauneuf, A.; Mitteau, J.-C. Structural Reliability; ISTE: London, UK, 2009; ISBN 9780470611708. [Google Scholar]
Figure 1. Left: Illustration of a masonry wallette, Right: Discontinuum representation of the masonry.
Figure 1. Left: Illustration of a masonry wallette, Right: Discontinuum representation of the masonry.
Modelling 01 00006 g001
Figure 2. Representation of the point contact.
Figure 2. Representation of the point contact.
Modelling 01 00006 g002
Figure 3. Representation of the point contact.
Figure 3. Representation of the point contact.
Modelling 01 00006 g003
Figure 4. Contact constitutive models. (a) Tension-compression stress–displacement behavior, (b) Shear stress–displacement behavior.
Figure 4. Contact constitutive models. (a) Tension-compression stress–displacement behavior, (b) Shear stress–displacement behavior.
Modelling 01 00006 g004
Figure 5. Boundary conditions; the tensile test of masonry parallel to bed joints.
Figure 5. Boundary conditions; the tensile test of masonry parallel to bed joints.
Modelling 01 00006 g005
Figure 6. Failure modes and macro behavior of masonry wallettes, obtained from the experiment [2]. (a) Left: Failure mode I, Right: Failure mode II, (b) Stress–displacement response of masonry wallettes, corresponding to different failure modes.
Figure 6. Failure modes and macro behavior of masonry wallettes, obtained from the experiment [2]. (a) Left: Failure mode I, Right: Failure mode II, (b) Stress–displacement response of masonry wallettes, corresponding to different failure modes.
Modelling 01 00006 g006
Figure 7. Histograms and probability distributions of the sampled random variables for WU-SB masonry assemblages.
Figure 7. Histograms and probability distributions of the sampled random variables for WU-SB masonry assemblages.
Modelling 01 00006 g007
Figure 8. Differences in the mean values after j number of simulations.
Figure 8. Differences in the mean values after j number of simulations.
Modelling 01 00006 g008
Figure 9. Differences in the standard deviation (SD) values after j number of simulations.
Figure 9. Differences in the standard deviation (SD) values after j number of simulations.
Modelling 01 00006 g009
Figure 10. Stress–displacement behavior of SU-WB and WU-SB masonry wallettes under tension parallel to bed joints with different material properties. (a) SU-WB (failure mode I), (b) WU-SB (failure mode II).
Figure 10. Stress–displacement behavior of SU-WB and WU-SB masonry wallettes under tension parallel to bed joints with different material properties. (a) SU-WB (failure mode I), (b) WU-SB (failure mode II).
Modelling 01 00006 g010
Figure 11. Fracture mechanism and stress–displacement behavior of SU-WB masonry wallettes (experimental result indicated via a thick red dashed line in part b). (a) Typical fracture mechanism of SU-WB, (b) Macro behavior of SU-WB wallettes (DEM and Experiment).
Figure 11. Fracture mechanism and stress–displacement behavior of SU-WB masonry wallettes (experimental result indicated via a thick red dashed line in part b). (a) Typical fracture mechanism of SU-WB, (b) Macro behavior of SU-WB wallettes (DEM and Experiment).
Modelling 01 00006 g011
Figure 12. Differences in the standard deviation values after j number of simulations.
Figure 12. Differences in the standard deviation values after j number of simulations.
Modelling 01 00006 g012
Figure 13. Different fracture mechanisms (FM) of the masonry wallettes subjected to tension parallel to bed joints (WU-SB). (a) Left: FM-1; Right: FM-2, (b) Left: FM-3; Right: FM-4, (c) FM-5.
Figure 13. Different fracture mechanisms (FM) of the masonry wallettes subjected to tension parallel to bed joints (WU-SB). (a) Left: FM-1; Right: FM-2, (b) Left: FM-3; Right: FM-4, (c) FM-5.
Modelling 01 00006 g013
Figure 14. Macro behavior of WU-SB wallettes (DEM results (FM-1) vs. experiment, which is plotted via a thick red dashed line).
Figure 14. Macro behavior of WU-SB wallettes (DEM results (FM-1) vs. experiment, which is plotted via a thick red dashed line).
Modelling 01 00006 g014
Figure 15. Effect of contact properties on the capacity for each fracture mechanism. (a) Influence of unit tensile strength ( f t , u ), (b) Influence of bond tensile strength ( f t , i ) and cohesion ( c ) , (c) Influence of unit–mortar interface friction angle ( ϕ ).
Figure 15. Effect of contact properties on the capacity for each fracture mechanism. (a) Influence of unit tensile strength ( f t , u ), (b) Influence of bond tensile strength ( f t , i ) and cohesion ( c ) , (c) Influence of unit–mortar interface friction angle ( ϕ ).
Modelling 01 00006 g015
Table 1. Random variables and their distributions for the strong unit–weak bond (SU-WB).
Table 1. Random variables and their distributions for the strong unit–weak bond (SU-WB).
Random VariableProb. Distribution Mean   ( μ ) Coeff. of Variation Lognormal   Mean   ( μ l n ) Lognormal   Std   Dev   ( σ l n )
E (GPa)Normal40.30N/AN/A
f t , u   ( MPa ) Lognormal1.00.45−0.09220.4294
f t , i   ( MPa ) Lognormal0.100.45−2.39480.4294
ϕ (degrees)Normal350.30N/AN/A
c   ( MPa ) Lognormal0.150.45−1.98930.4294
ψ (degrees)Lognormal20.300.65010.2936
Table 2. Random variables and their distributions for weak unit–strong bond (WU-SB).
Table 2. Random variables and their distributions for weak unit–strong bond (WU-SB).
Random VariableProb. Distribution Mean   ( μ ) Coeff. of Variation Lognormal   Mean   ( μ l n ) Lognormal   Std   Dev   ( σ l n )
E (GPa)Normal40.30N/AN/A
f t , u   ( MPa ) Lognormal0.250.45−1.47850.4294
f t , i   ( MPa ) Lognormal0.200.45−1.70160.4294
ϕ (degrees)Normal350.30N/AN/A
c   ( MPa ) Lognormal0.300.45−1.29620.4294
ψ (degrees)Lognormal20.300.65010.2936

Share and Cite

MDPI and ACS Style

Pulatsu, B.; Gonen, S.; Erdogmus, E.; Lourenço, P.B.; Lemos, J.V.; Hazzard, J. Tensile Fracture Mechanism of Masonry Wallettes Parallel to Bed Joints: A Stochastic Discontinuum Analysis. Modelling 2020, 1, 78-93. https://0-doi-org.brum.beds.ac.uk/10.3390/modelling1020006

AMA Style

Pulatsu B, Gonen S, Erdogmus E, Lourenço PB, Lemos JV, Hazzard J. Tensile Fracture Mechanism of Masonry Wallettes Parallel to Bed Joints: A Stochastic Discontinuum Analysis. Modelling. 2020; 1(2):78-93. https://0-doi-org.brum.beds.ac.uk/10.3390/modelling1020006

Chicago/Turabian Style

Pulatsu, Bora, Semih Gonen, Ece Erdogmus, Paulo B. Lourenço, Jose V. Lemos, and Jim Hazzard. 2020. "Tensile Fracture Mechanism of Masonry Wallettes Parallel to Bed Joints: A Stochastic Discontinuum Analysis" Modelling 1, no. 2: 78-93. https://0-doi-org.brum.beds.ac.uk/10.3390/modelling1020006

Article Metrics

Back to TopTop