Microsoft Word - numero_33_art_11 A.Spagnoli et alii, Frattura ed Integrità Strutturale, 33 (2015) 80-88; DOI: 10.3221/IGF-ESIS.33.11 80 Focussed on characterization of crack tip fields Interpreting experimental fracture toughness results of quasi-brittle natural materials through multi-parameter approaches Andrea Spagnoli, Andrea Carpinteri, Sabrina Vantadori Dipartimento di Ingegneria Civile, dell’Ambiente del Territorio e Architettura, Università di Parma, Italy spagnoli@unipr.it ABSTRACT. Natural stones like marbles are often employed as façade panels to externally cover buildings. These natural materials tend to exhibit a quasi-brittle nonlinear fracture behaviour which, conversely to concrete counterpart, has much less been studied in the literature. In the present paper, an experimental campaign on the so-called red Verona marble is carried out, and the results are discussed together with some previously published results on the white Carrara marble. The analysis of the two marbles as a whole allows us to discuss size effect and to point out the need for additional parameters in order to describe their fracture behaviour. The study focuses on a two-parameter model which accounts for a characteristic material length. A quantitative correlation between material microstructure features, obtained from thin sections technique, and the characteristic material length is proposed. KEYWORDS. Fracture process zone; Marble; Quasi-brittle behaviour; Size effect; Two-parameter model. INTRODUCTION uasi-brittle materials are typically characterized by a stable nonlinear load-deflection response followed by an unstable fracture. Such a behaviour is generally attributed to the development of a finite-size process zone ahead of the crack tip. The quasi-brittle behaviour has largely been examined for concrete both experimentally and theoretically [1]. When fracture toughness in quasi-brittle materials has to be determined, an effective crack length (spanning a traction-free segment and a stressed segment embedded in the process zone) is needed to be defined, with respect to which the critical stress intensity factor has to be computed. Accordingly, in the realm of concrete material testing, Shah and co-workers [2] proposed a two-parameter method to calculate the effective crack length in the verge of unstable cracking. Other theories are available in the literature to describe experimental size-dependent fracture data. These theories (e.g. theory of critical distances [3,4], finite fracture mechanics [5], cohesive crack models [6], size effect laws [7-10], non-local stress failure criteria [11], and crack extension resistance curves [12]) are characterized by a fil rouge, that is to say, the near-crack tip stress field and/or the fracture toughness is described by means of multi-parameter approaches. A size-dependence discussion of the RILEM work-of-fracture method can be found in Ref. [13]. Natural materials, like calcite marbles used for building claddings [14, 15], often exhibit a quasi-brittle mechanical behaviour [16, 17] which, conversely to concrete counterparts, has much less been investigated despite the fact that, under different deteriorating in-service actions (temperature cycles, wind and seismic loads, chemical actions, humidity changes, etc.), fracture behaviour might significantly affect structural failure. Q http://www.gruppofrattura.it/pdf/rivista/numero33/Audioslides/Spagnoli1/Sp1.html A.Spagnoli et alii, Frattura ed Integrità Strutturale, 33 (2015) 80-88; DOI: 10.3221/IGF-ESIS.33.11 81 In the present paper, an experimental campaign on the so-called red Verona marble is performed, and the results are examined with previously published results on the white Carrara marble [18]. The analysis of the two marbles allows us to discuss size effect and to point out the need for additional parameters in order to describe their fracture behaviour. The study focuses on a two-parameter model which accounts for a characteristic material length. A tentative quantitative correlation between material microstructure features, obtained from thin sections technique, and the characteristic material length is presented. MATERIAL MICROSTRUCTURAL CHARACTERIZATION AND EXPERIMENTAL MECHANICAL TESTS he results of some fracture experimental tests on two different types of natural stones, used as cladding elements in building façades, are presented. In particular, white Carrara marble, which is extracted from Carrara area in the Alpi Apuane metamorphic complex, and red Verona marble, which is a variety of limestone sedimentary rock taking its name from the city of Verona in Northern Italy, are analysed. For the former material, the results published in Refs [18,19] are herein elaborated, while the experimental results for the latter are presented in this paper for the first time. A picture of a thin section of Carrara marble is shown in Fig. 1. From a qualitative viewpoint, the picture indicates common features of calcite grains with granoblastic texture and anhedral grain shape. The grain size distribution is characterized by the prevalence of large grains (about 200 m on average), followed by smaller grains (30–50 m). Grain boundaries are typically ranging from polygonal shape to interlobate one in both large and small crystals. A third group of calcite grains is detected as fine grained aggregates (<10 m) commonly rimming large porphyroblasts. Calcite crystals are characterized by various degree of plastic deformation, as is demonstrated by the occurrence of grains with extinction features from homogeneous to undulated; often twin lamellae occur. Twin planes or cleavage planes are well developed in most of the large grains and occur with various orientations. Two sets of twin lamellae may also occur in the same grains, testifying the activation of both the twinning planes of calcite. Locally, cleavage planes may become fracture planes never crossing through adjacent grains. No anisotropic features are present in the orthogonal sections in terms of either grain size variations, shape preferred orientation or compositional layers. Figure 1: Carrara marble microstructure observed by a thin sections technique: (a) parallel nicol; (b) crossed nicol. Fig. 2 shows a thin section of Verona marble. The typical structure of a sedimentar calcareous rock is evident, with pale pink nodules set in a darker pink matrix. Nodules show a wackestone texture, with skeletal grains (thin-shelled bivalves, foraminifera, gastropods, ostracodes) set in an abundant (>70%) micritic matrix (microcrystalline calcite with average size lower than 5mm). The matrix shows a packstone texture, with compacted thin-shelled bivalves (>70%) and interstitial micritic matrix. Porosity is very low (< 1%, visive estimate in thin section) with holes, variable in shape and size, empty or cemented by a mosaic of sparite (calcite crystals ranging from 10 μm to 500 μm in size). The tensile strengths of the analysed marbles are determined from bending tests on prismatic specimens. Ten 4-point bend tests were performed on Carrara marble [18]. The nominal sizes of the specimens are 180mm in span, 60mm in depth and 30mm in width. Further, twenty 3-point bend tests are performed on Verona marble, where the nominal sizes of the specimens are 150mm in span, 30mm in depth and 60mm in width. A summary of the tensile strength is reported in Tab. 1. It can be noted that the mean value of strength is similar for the two marbles, with the Verona marble showing T A.Spagnoli et alii, Frattura ed Integrità Strutturale, 33 (2015) 80-88; DOI: 10.3221/IGF-ESIS.33.11 82 a comparatively larger scatter in the results. This trend can be attributed to the markedly heterogeneous nature of Verona marble, which presents weak phases of clay inclusions. Figure 2: Verona marble microstructure observed by a thin sections technique: (a) parallel nicol; (b) crossed nicol. Mean value (MPa) Stand. Dev. (MPa) Coeff. Variation White Carrara marble 14.91 2.97 0.20 Red Verona marble 12.27 3.55 0.29 Table 1: Tensile strength of the two marbles, obtained from bending tests Fracture tests are then considered for the two marbles. Three-point bend tests on single edge notched beams have been performed (Fig. 3). The nominal sizes of Carrara marble specimens are [18]: depth W = 60mm, width B = W/2, span S = 3W, initial notch length a0 = W/10. The nominal sizes of Verona marble specimens are: W = 30mm, B = W/2, S = 4W, initial notch length a0 = W/3. The relative crack length is defined as /a W  . Figure 3: Three-point bend notched specimen. During the tests, the Crack Mouth Opening Displacement (CMOD) is recorded via a clip gauge mounted across the edge notch of the beam, using two metallic gauge holders. The tests are performed through a servo-mechanical Instron 8862 machine, under a CMOD controlled procedure (CMOD rate in the range of 0.50-0.80 m/min) in order to be able to follow the softening branch beyond the peak load. In some tests, unloading-reloading cycles in the softening phase are performed. The experimental curves of nominal stress N (equal to  26 4PS BW ) against CMOD are reported in Fig. 4 for Carrara marble and Verona marble. S W a0 P a CMOD CTOD FPZ a0 x A.Spagnoli et alii, Frattura ed Integrità Strutturale, 33 (2015) 80-88; DOI: 10.3221/IGF-ESIS.33.11 83 Figure 4: Experimental N-CMOD curves for: (a) Carrara marble (6 tests); (b) Verona marble (7 tests). In Tab. 2 and 3, fracture test results obtained from the elaboration of the curves plotted in Fig. 4 are illustrated. In particular, the peak nominal stress and the initial crack length are reported in the 2nd and 3rd column of both tables (the other quantities reported in Tab. 2 and 3 are described in the following Sections). Using the classical Linear Elastic Fracture Mechanics (LEFM), the nominal critical stress intensity factor ICK can be worked out on the basis of the initial notch length (relative length, 0  ). The following expression is used [20]:  I NK aF   (1) where         2 1.5 1.99 1 2.15 3.93 2.71 1 2 1 F                (2) It can be noted that the mean values of the nominal critical stress intensity factors for the two marbles are very similar. Specimen No. N (MPa) a0 (mm) ICK (MPa√m) E (GPa) a-a0 (mm) s ICK (MPa√m) CTODc (m) Q (mm) 1 7.06 6.07 0.98 44.6 5.35 1.32 6.9 54.98 2 6.87 6.20 0.96 28.8 5.82 1.32 11.3 60.50 3 6.95 6.40 0.99 31.1 5.29 1.31 9.9 54.86 4 7.04 6.80 1.03 50.1 3.28 1.23 4.4 32.68 5 7.38 6.57 1.06 42.6 0.88 1.12 2.4 8.33 6 6.84 6.87 1.01 85.3 3.83 1.24 4.2 38.97 Mean 7.02 6.49 1.01 42.6 4.08 1.26 6.5 41.72 Stand. Dev. 0.19 0.32 0.04 11.3 1.85 0.08 3.5 19.52 Coeff. Var. 0.03 0.05 0.04 0.26 0.45 0.06 0.5 0.47 Table 2: Results of fracture tests for white Carrara marble [18] (a) (b) A.Spagnoli et alii, Frattura ed Integrità Strutturale, 33 (2015) 80-88; DOI: 10.3221/IGF-ESIS.33.11 84 Specimen No. N (MPa) a0 (mm) ICK (MPa√m) E (GPa) a-a0 (mm) s ICK (MPa√m) CTODc (m) Q (mm) 1 7.53 8.74 1.30 53.3 0.92 1.41 2.7 10.49 2 5.08 9.74 0.96 52.9 3.16 1.29 5.1 44.81 3 7.43 8.93 1.30 45.9 0.69 1.38 2.7 7.82 4 3.66 9.74 0.69 42.6 1.25 0.77 2.2 15.09 5 5.07 9.08 0.90 60.3 1.34 1.02 2.1 15.95 6 3.17 10.16 0.62 28.9 0.72 0.66 2.1 8.44 7 7.95 9.59 1.48 50.9 0.68 1.57 2.7 7.76 Mean 5.70 9.43 1.04 47.8 1.25 1.16 2.8 15.77 Stand. Dev. 1.95 0.52 0.33 10.1 0.89 0.35 1.1 13.25 Coeff. Var. 0.34 0.05 0.32 0.21 0.71 0.30 0.4 0.84 Table 3: Results of fracture tests for red Verona marble. PITFALL OF SIMILITUDE AND SIZE EFFECT he results of strength and toughness for the two marbles seem to indicate that they have similar mechanical characteristics. This might be used instrumentally in the following to tentatively treat the two marbles as a whole, i.e. as a single material with the same strength and toughness but different material microstructure features and different specimen and crack sizes. In particular, the nominal peak stress of the notched specimens multiplied by the dimensionless factor in Eq. 2 is plotted against the initial notch length in Fig. 5. Such an adjustment of the nominal stress is for accounting the different relative initial notch lengths in Carrara marble and Verona marble specimens under consideration. In this way, the dependence on 1 a as the crack length is made to vary can be isolated. From the observation of the results in Fig. 5, it can be noted that there is a crack length dependence of the strength, but this deviates from the -0.50 exponent of LEFM. As a matter of fact, considering the mean values of nominal peak stress and the nominal initial notch length for the two marbles, a power law fitting gives us an exponent equal to -0.29 which is in absolute value smaller than that of LEFM. This might be an indication of a deviation from LEFM and, hence, of a pitfall of similitude. In other words, nonlinear phenomena are likely to take place and, therefore, some additional parameters are needed to describe the fracture behaviour of such materials. Figure 5: Distribution of the adjusted nominal peak stress against the initial notch length for the two marbles examined. The power law fitting is based on the mean values of strength (star symbols) at nominal initial notch length (fitting slope equal to -0.29). Mean values of tensile strength (see the two symbols at a0 = 0) are reported as a reference. T A.Spagnoli et alii, Frattura ed Integrità Strutturale, 33 (2015) 80-88; DOI: 10.3221/IGF-ESIS.33.11 85 MULTI-PARAMETER APPROACHES TO INTERPRET QUASI-BRITTLE BEHAVIOUR imilar to concrete, which is the archetypal quasi-brittle structural material, natural stones are also characterized by a quasi-brittle fracture behaviour. Therefore, the research findings of the concrete community might be borrowed when natural stones are investigated. As is well known, the quasi-brittle behaviour is due to, different from the ductile behaviour of metals, the development of a sizable nonlinear zone ahead of the fracture front. Such a zone is almost entirely filled by a Fracture Process Zone (FPZ, see Fig. 3), whose size is equal or proportional to a characteristic (material) length l. Depending on the size of the structural component, the FPZ may affect different fractions of the resistant cross-section, or it might be even larger than the cross-section itself. For instance (see Ref. [1]), l can encompass order of magnitudes, e.g. from 10-100nm for a silicone waver to about 3m for a dam concrete up to 10m for a grouted soil mass. If D/l (D = characteristic size of the resistant cross-section) is, say, greater than 100, the FPZ size is negligible with respect to that of the structure, and LEFM holds. On the contrary, when D/l is in the range of about 5-100, the nonlinear behaviour taking place in the FPZ cannot be disregarded and, hence, size effects intervene and multi-parameter approaches are required when fracture energy and toughness have to be determined. The simplest conceptual model to characterize the nonlinear behaviour of the finite FPZ is the cohesive model by Hillerborg [21]. Accordingly, the FPZ for Mode I fracture can be described by a fictitious line crack that transmits normal stress  which is a monotonically decreasing function of the crack opening displacement w, namely  f w  . By definition,  0f is equal to the tensile strength tf of the material and, at a certain critical value cw of the crack opening displacement, the transmitted force across the crack becomes null, i.e.   0cf w  . The area underneath the curve  f w  is equal to the total energy ( FG ) dissipated by fracture per unit area of crack plane, as the crack faces are completely detached at a given point. The quantity FG has been proved to be equal to the energy dissipated during the crack extension (the latter being the energy which has to be equated, according to Griffith, to the energy release rate of the structure) [22]. Several shapes of the cohesive law  f w  are available in the literature, since the decaying exponential law of Hillerborg [21]. The most popular law for concrete is probably a bilinear law presented in Ref. [23]. It is important to stress at this point is that, regardless of the shape of the  f w  law, fracture resistance is described by (at least) two parameters ( tf and FG ) due to the sizable FPZ of quasi-brittle materials. In addition, as was pointed out by Planas et al. [24], the initial negative slope in the softening stress-separation law (sometimes it is useful to isolate a fracture energy fG defined by such an initial slope) controls the peak load of a structure and hence size effect, so that three parameters are needed to identify fracture. Quasi-brittle behaviour is well captured by cohesive crack models which are nonlinear in nature. This implies that any attempt to directly apply LEFM to quasi-brittle materials needs some adaptation. Therefore, in addition to fracture toughness related to fracture energy through the Irwin relationship ( IC FK E G , with E E  for plane stress and  21E E    for plane strain), at least one material parameter is needed. The easiest way of adapting LEFM is to define an equivalent crack (regarded as a sharp traction free crack) whose tip is located at some distance ahead the real tip in the FPZ ( 0a a a   ), where a is a material parameter. This directly implies that some non-zero crack tip opening displacement (crack blunting) occurs at the end of the length 0a , that is to say, at the beginning of the FPZ. One possible definition of the material parameter a is related to the Irwin-like inelastic zone size 2 1 IC t K f       (e.g. see Ref. [3]). Alternatively, the additional material parameter can be based on the attainment of a critical value of CTOD ( cw ) as is suggested by the Jenq-Shah model [2] (that is shown in the next Section). S A.Spagnoli et alii, Frattura ed Integrità Strutturale, 33 (2015) 80-88; DOI: 10.3221/IGF-ESIS.33.11 86 APPLICATION OF THE TWO-PARAMETER MODEL BY JENQ AND SHAH his model [2], reminiscent of that by Wells-Cottrell for metals [25], is based on the observation that the initial crack length a0 in quasi-brittle materials grows slowly well before the peak load is attained. This nonlinear stable stage terminates when the CTOD reaches a critical value (in other words, when the crack propagates to a critical extent) and IK attains a value s ICK which differs from the nominal ICK (calculated on the basis of a0). If the geometric and loading conditions are such that, the stress intensity factor is monotonically increasing with the crack length (being the load constant), as the critical condition explained before takes place at the peak load in the case of a 3-point bend beam with an edge crack. From the value of a, using the LEFM formula, we have:    4 0 N a CMOD w x V E     (3) where     2 3 2 0.66 0.76 2.28 3.87 2.04 1 V            (4) The opening displacement along the crack can be expressed as follows [2]:     2 2 1 1.081 1.149 x x x COD w x CMOD a a a                      (5) Considering the compliance parameter C, defined as C CMOD P , we can write from Eq. 3 (obviously P is equal to  24 6N BW S ):  2 6Sa C V EBW  (6) Two compliances measurements, at the initial crack length and at the unknown crack length at peak load, are needed to determine the required parameters. From the former, the Young modulus E is calculated and, from the latter, the effective crack length a at failure is determined. In this way the equivalent crack length (containing the additional parameter 0a a a   ) is estimated and, using the adapted LEFM at the peak load (Eq. 1), the effective toughness s ICK is worked out. This effective toughness is larger than the nominal one based on the initial crack length a0 and the peak load, and is dependent on the material parameter a . In Fig. 6 for Verona marble, two examples of calculation of initial compliance iC and failure compliance uC are shown. In particular, one example is based on the secant compliance at peak load, whereas the other example on the tangent compliance in the unloading-reloading cycle after peak load. Posing 0x a , we can calculate CTOD from Eq. 6 and, considering the peak load in Eq. 3, the critical value CTODc can be obtained. The additional parameter implicitly introduced in the Jenq-Shah model can be expressed by the length parameter Q, defined as follows [2]: 2 c s IC E CTOD Q K        (7) The results of the two marbles (see Tab. 2 and 3) show that, in Carrara marble, a significant increase of the initial crack length can be observed prior to reaching the critical condition at peak load ( 0a a =4.08mm on average), which is about 20 times the mean calcite grain size (about 200 m on average), regarded as a characteristic material length. In addition, the length parameter Q is on average equal to 41.72mm, that is, a Q/W = 0.70. On the other hand, Verona marble exhibits a comparatively smaller crack growth with respect to the initial length, with an extent of 1.25mm on average, which is about 250 times the mean microstructure size of the material (about 5 m on T A.Spagnoli et alii, Frattura ed Integrità Strutturale, 33 (2015) 80-88; DOI: 10.3221/IGF-ESIS.33.11 87 average; note that, in this case, the identification of a single material characteristic length from thin section is not straightforward). The length parameter Q is on average equal to 15.77mm, which corresponds to Q/W = 0.53. Figure 6: Calculation of the actual crack length for two tests related to Verona marble: (a) compliance at peak load (specimen 1); (b) compliance during unloading-loading cycle in the post-peak phase (specimen 2). The trend of results shows that there is an increase of both stable crack growth and length parameter Q with the increasing of the characteristic length (of the two marbles) estimated from thin sections technique. This trend is in line with the fact that the size of FPZ is proportional to microstructure characteristic length of the material. In addition, the values of Q/W, together with the relative initial notch length (equal to 0.10 and 0.33 for Carrara marble and Verona marble specimens, respectively), indicate a roughly similar size effect for the two marble specimens [2], as is confirmed by the obtained ratio of nominal against effective critical stress intensity factors (on average equal to 0.80 and 0.89 for Carrara marble and Verona marble specimens, respectively). Instead, the Irwin-like inelastic zone size does not seem to be in line with the trend of microstructure characteristic length of the material (on average equal to 1.45mm and 2.27mm for Carrara marble and Verona marble specimens, respectively). CONCLUSIONS ifferent from the concrete behaviour, the quasi-brittle fracture behaviour of natural rocks has attracted much less attention. In the present paper, the experimental results of the so-called Verona marble are presented and discussed with some previously published results on Carrara marble. Despite their roughly similar nominal tensile strength and fracture toughness, the two marbles are characterized by very different characteristic length of material microstructure, spanning about two order of magnitudes (with Carrara marble having the larger length). A two- parameter model proposed by Shah and coworkers seems to be suitable to describe the experimental results under discussion. In particular, it is shown that the additional length parameter, represented by the effective crack length, required to describe the nonlinear fracture behaviour, appears to be consistent with the characteristic size of the material microstructure. REFERENCES [1] Bazant, Z.P., Concrete fracture models: testing and practice. Engng Frac. Mech., 69 (2002) 165-205. [2] Jenq, Y., Shah, S.P., Two Parameter Fracture Model for Concrete. J. Eng. Mech., 111 (1985) 1227-1241. [3] Susmel, L., Taylor, D., The Theory of Critical Distances as an alternative experimental strategy for the determination of KIc and Kth. Engng Frac. Mech., 77 (2010) 1492-1501. [4] Askes, H., Susmel, L., Understanding cracked materials: is Linear Elastic Fracture Mechanics obsolete?, Fatigue Fract. Engng Mater. Struct., (2014). Available on-line. doi: 10.1111/ffe.12183. D (a) (b) A.Spagnoli et alii, Frattura ed Integrità Strutturale, 33 (2015) 80-88; DOI: 10.3221/IGF-ESIS.33.11 88 [5] Taylor, D., Cornetti, P., Pugno, N., The fracture mechanics of finite crack extension. Engng Frac. Mech., 72 (2005) 1021-1038. [6] Li, V.C., Maalej, M., Toughening in cement based composites. Part I: Cement, mortar, and concrete, 18 (1996) 223- 237. [7] Bazant, Z.P., Size effect in blunt fracture: concrete, rock, metal. J. Engng Mech. ASCE, 110 (1984) 518-535. [8] Carpinteri, Al., Scaling laws and renormalization groups for strength and toughness of disordered materials. Int. J. Solids Structures, 31 (1994) 291-302. [9] Carpinteri, An., Spagnoli, A., A fractal analysis of size effect on fatigue crack growth, Int. J. Fatigue, 26(2) (2004) 125- 133. [10] Carpinteri, An., Spagnoli, A., Vantatori, S., A multifractal analysis of fatigue crack growth and its application to concrete, Engng Frac. Mech., 77(6) (2010) 974-984. [11] Gomez, F.J., Elices, M., Berto, F., Lazzarin, P., Local strain energy to assess the static failure of U-notches in plates under mixed mode loading, Int. J. Fract., 145 (2007) 29-45. [12] Xu, S., Reinhardt, H.W., Crack Extension Resistance and Fracture Properties of Quasi-Brittle Softening Materials like Concrete Based on the Complete Process of Fracture, Int. J. Fat., 92 (1998) 71-99. [13] Bazant, Z.P., Analysis of work-of-fracture method for measuring fracture energy of concrete, J. Engng Mech. ASCE, 122 (1996) 138-144. [14] Ferrero, A.M., Migliazza, M., Spagnoli, A., Theoretical modelling of bowing in cracked marble slabs under cyclic thermal loading. Construction and Building Materials, 23 (2009) 2151-2159. [15] Spagnoli, A., Ferrero, A.M., Migliazza, M., A micromechanical model to describe thermal fatigue and bowing of marble. Int. J. Solids Struct., 48 (2011) 2557-2564. [16] Alber, M., Hauptfleisch, U., Generation and visualization of microfractures in Carrara marble for estimating fracture toughness, fracture shear and fracture normal stiffness, Int. J. Rock Mechanics & Mining Sci., 36 (1999) 1065-1071. [17] Aliha, M.R.M., Ayatollahi, M.R., Smith, D.J., Pavier, M.J., Geometry and size effects on fracture trajectory in a limestone rock under mixed mode loading, Engng Fract. Mech., 77 (2010) 2200-2212 [18] Migliazza, M., Ferrero, A.M., Spagnoli, A., Experimental investigation on crack propagation in Carrara marble subjected to cyclic loads. Int. J. Rock Mechanics & Mining Sci., 48 (2011) 1038-1044. [19] Ferrero, A.M., Migliazza, M., Spagnoli, A., Zucali, M. Micromechanics of intergranular cracking due to anisotropic thermal expansion in calcite marbles, Engineering Fracture Mechanics, 130 (2014) 42-52. [20] Tada, H., Paris, P.C., Irwin, G.R., The Stress Analysis of Crack Handbook. Del Research Corporation, St.Louis (1985). [21] Hillerborg, A., Modéer, M., Petersson, P.E., Analysis of crack formation and crack growth in concrete by means of fracture mechanics and finite lements. Cement Concrete Res, 6 (1976) 773-782. [22] Rice, J.R., Mathematical analysis in the mechanics of fracture. In: Liebowitz H, editor. Fracture – an advanced treatise, New York: Academic Press., 2 (1968) 191-308. [23] Petersson, P.E., Crack growth and development of fracture zones in plain concrete and similar materials. Report TVBM-1006. Division of Building Materials, Lund Institute of Technology, Lund, Sweden, (1981). [24] Planas, J., Elices, M., Guinea, G.V., Measurement of the fracture energy using three-point bend tests: Part 2 – Influence of bulk energy dissipation. Mater. Struct., 25 (1992) 305-312. [25] Wells A.A. Unstable crack propagation in metals-clevage and fast fracture. Symposium on Crack Propagation, Cransfield., 1 (1961) 210-230.