Microsoft Word - numero_34_art_12 L. Náhlík et alii, Frattura ed Integrità Strutturale, 34 (2015) 116-124; DOI: 10.3221/IGF-ESIS.34.12 116 Focussed on Crack Paths Estimation of stepwise crack propagation in ceramic laminates with strong interfaces L. Náhlík, K. Štegnerová Institute of Physics of Materials, Academy of Sciences of the Czech Republic, Zizkova 22, 616 62 Brno, Czech Republic and Faculty of Mechanical Engineering, Brno University of Technology, Technicka 2, 616 69 Brno, Czech Republic nahlik@ipm.cz, stegnerova@ipm.cz P. Hutař Institute of Physics of Materials, Academy of Sciences of the Czech Republic, Zizkova 22, 616 62 Brno, Czech Republic hutar@ipm.cz ABSTRACT. During the last years many researchers put so much effort to design layered structures combining different materials in order to improve low fracture toughness and mechanical reliability of the ceramics. It has been proven, that an effective way is to create layered ceramics with strongly bonded interfaces. After the cooling process from the sintering temperature, due to the different coefficients of thermal expansion of individual constituents of the composite, significant internal residual stresses are developed within the layers. These stresses can change the crack behaviour. This results to the higher value of so-called apparent fracture toughness, i.e. higher resistance of the ceramic laminate to the crack propagation. The contribution deals with a description of the specific crack behaviour in the layered alumina-zirconia ceramic laminate. The main aim is to clarify crack behaviour in the compressive layer and provide computational tools for estimation of crack behaviour in the field of strong residual stresses. The crack propagation was investigated on the basis of linear elastic fracture mechanics. Fracture parameters were computed numerically and by author’s routines. Finite element models were developed in order to obtain a stress distribution in the laminate containing a crack and to simulate crack propagation. The sharp change of the crack propagation direction was estimated using Sih’s criterion based on the strain energy density factor. Estimated crack behaviour is qualitatively in a good agreement with experimental observations. Presented approach contributes to the better understanding of the toughening mechanism of ceramic laminates and can be advantageously used for design of new layered ceramic composites and for better prediction of their failure. KEYWORDS. Ceramic laminates; Crack behaviour; Residual stresses; Strain energy density factor; Crack propagation direction. INTRODUCTION aterial interfaces play an important role in material (composite) behaviour and its properties. The presence of regions with different material and mechanical properties and the existence of an interface between them have a pronounced influence on the stress distribution of composite bodies. The influence of material interfaces was M L. Náhlík et alii, Frattura ed Integrità Strutturale, 34 (2015) 116-124; DOI: 10.3221/IGF-ESIS.34.12 117 intensively studied during last years for different material combinations, from both theoretical and experimental sides, see e.g. [1-8] for detailed information. The specific stress distribution of layered composites can be used in design of new composite materials with improved mechanical properties, e.g. with higher flexural strength or higher resistance to crack propagation, etc. The existence of interfaces and their toughening effect are, with the advantage, used in the manufacturing of “flaw tolerant” multilayer ceramic composites [9]. There is a new class of ceramic materials under development, the ceramic laminates with controlled crack resistance and crack propagation [10-12]. The paper is focused on ceramic laminates with strongly bonded layers. The toughening mechanism of strongly bonded laminates is based on deflection and/or bifurcation of crack path caused by internal residual stresses inside the layers [13,14]. Mechanism of internal stresses evolution is based on unequal shrinkage of physically bonded materials during sintering [15]. This mechanism was experimentally proved for different kinds of ceramic laminates with strongly bonded interfaces [16-27], however only few works exists dealing with the prediction of laminate properties [28] or explaining specific crack behaviour in the ceramic laminates [29-32]. This paper focuses on the crack behaviour in the compressive layer of the laminate and describes by use of linear elastic fracture mechanics specific crack propagation (see Fig. 1) in this layer. Figure 1: Crack bifurcation in AMZ layer. With courtesy of R. Bermejo [23]. The basic principle of the toughening of the layered ceramic composite is based on selection of two ceramics with different values of coefficients of thermal expansion. The ceramics are layered by system A-B-A symmetrically and parallelly to the axis of symmetry of the layered composite body. Cooling down from the sintering temperature leads to the development of residual stresses, which are usually of tensile nature in the outside layers. Description of the studied layered ceramic composite For the study material system based on alumina and zirconia was chosen. These two constituents can be taken as typical examples of convenient materials for layered ceramics. Behaviour of such kind of laminates were studied e.g. in [20-23]. Considered ceramic laminate was created by 9 layers. 5 were made of Al2O3/5vol.%t-ZrO2 (alumina with tetragonal zirconia, reffered as ATZ) and 4 were made of Al2O3/30vol.%m-ZrO2 (alumina with monoclinic zirconia, reffered as AMZ). Material characteristics were taken from literature [20,21] and are summarized in Tab. 1. The laminate was subjected to four-point bending (4PB) test, see Fig. 2. Property ATZ AMZ Young´s modulus E [GPa] 390 280 Poisson´s ratio ν [-] 0.22 0.22 Fracture toughness KIC [MPa.m0.5] 3.2 2.6 Coefficient of thermal expansion t [10 -6 K -1 ] 9.8 8.0 Strength f [MPa] 422 90 Table 1: Material properties of studied laminate [20,21]. L. Náhlík et alii, Frattura ed Integrità Strutturale, 34 (2015) 116-124; DOI: 10.3221/IGF-ESIS.34.12 118 Figure 2: Scheme of studied ceramic laminate with typical crack path (dimensions in [mm]) and orientation of residual stresses in composite layers. The magnitude of residual stresses in the case of crack absence in the layered composite can be estimated from forces balance in individual layers. This leads to following expressions [33]:  0 , 1 1 1 ' ' 1 sf T AMZ ATZ T res ATZ ATZ AMZ dT N E E N             , , , 1 1 res AMZ res ATZ N N          , (1) where: ATZ AMZ t t   , ' 1 AMZ AMZ AMZ E E    , ' 1 ATZ ATZ ATZ E E    , Tsf [°C] is stress free temperature, T0 [°C] room temperature, α [10-6K-1] coefficient of thermal expansion, N [-] number of layers, t [mm] thickness of layer, E [MPa] Young´s modulus and ν [-] Poisson´s ratio. Calculated values of residual stresses for studied configurations are shown in the Tab. 2. tATZ : tAMZ tATZ [mm] tAMZ [mm] σres,ATZ [MPa] σres,AMZ [MPa] 2 : 1 0.4288 0.2140 236.6 -592.6 5 : 1 0.5170 0.1038 109.8 -683.6 7 : 1 0.5384 0.0770 80.6 -704.6 10 : 1 0.5556 0.0556 57.7 -721.0 Table 2: Residual stresses calculated for different ratio of layer thicknesses (tATZ : tAMZ) by Eq. 1. It is evident from the Tab. 2 that the compressive residual stresses in the AMZ layers are circa six hundred megapascals and higher for all considered cases. These stresses are responsible for typical stepwise crack propagation like is shown in the scheme in Fig. 2. MODEL OF CRACK PROPAGATION Estimation of crack propagation direction ih´s criterion [34,35] based on strain energy density factor S was used for the determination of crack propagation direction due to brittle nature of the composite. The criterion was implemented in finite element system Ansys [36] to perform calculations with small crack increments automatically: 2 211 12 222I I II IIS a K a K K a K   , (2) where KI and KII are stress intensity factors corresponding to mode I and II respectively, S L. Náhlík et alii, Frattura ed Integrità Strutturale, 34 (2015) 116-124; DOI: 10.3221/IGF-ESIS.34.12 119   11 1 1 cos cos 16 a k        ,  12 1 sin 2 cos 1 16 a k        ,      22 1 1 1 cos 1 cos 3cos 1 16 a k             , constant 3 4k   for plane strain or 3 1 k      for plane stress,  2 1 E     is shear modulus,  is polar coordinate originating from the crack tip. The criterion postulates that the crack will propagate in the direction where factor S is the minimum. The angle of next crack propagation direction 0 can be determined from the conditions: 0 0 2 2 0  0 S S                      (3) The crack will propagate under angle 0 if the value of S will reach its critical value Scr, i.e. for S = Scr. In special case (pure mode I) can be the critical value Scr related to fracture toughness KIc:   21 2 4 IC cr K S     for plane strain condition. (4) Fracture toughness of AMZ layer is 2.6 MPa.m0.5. Corresponding value of Scr = 825.10-8 MPa.m. Figure 3: Detail of fine mesh around the crack tip. Numerical model For the numerical modelling finite element (FE) method was chosen. Numerical models were developed in commercial system Ansys. Models contained circa 100 000 elements PLANE 183. Boundary conditions corresponding to 4PB, i.e. vertical displacements on the lower side in the locations of rigid supports were equal to zero and the loading forces P/2 on the upper side of the models were applied, see Fig. 2 for details. Thermal load by change of temperature from 1200°C to 20°C was applied together with mechanical loading to generate residual stresses in the layers. Material characteristics written in the Tab. 1 were used. All calculations were performed under condition of plane strain. Very fine mesh was used around the crack tip to well describe the stress distribution there (see Fig. 3), which is of crucial importance for determination of crack behavior. The size of smallest elements was circa 1.10-5 mm. On the base of former experience [28] the crack propagation was modelled from first ATZ/AMZ interface through the AMZ layer. The initial crack L. Náhlík et alii, Frattura ed Integrità Strutturale, 34 (2015) 116-124; DOI: 10.3221/IGF-ESIS.34.12 120 increment was chosen 1 m to describe crack path precisely and model was remeshed after each step. Stress intensity factors were calculated by author’s routine from displacements close to the crack tip in each step. New crack propagation direction was determined on the base of author´s Ansys macros based on Eqs. 2 and 3 after each step. Hundreds of calculations were performed to obtained realistic crack path in compressive AMZ layer during each simulation. PC based workstations were used for extensive numerical calculations. RESULTS AND DISCUSSION ifferent material layers create regions with different material properties. It was shown in the former work of authors [7] that suitable selection of materials with different elastic properties can lead to the higher values of applied load for crack propagation through interface. This value of applied load can be higher than the one for crack propagation in the individual composite constituents. The effect of crack retardation is stronger in the studied case due to acting of residual stresses developed in the layers during sintering process. These residual stresses have pronounced influence on the crack behaviour in the composite, damage mechanism of the composite and consequently on the value of apparent fracture toughness of the composite body. In the numerical FE studies four different ratios of layer thicknesses were considered (see Tab. 2). The calculations performed were focused on the description of mechanism leading to the higher apparent fracture toughness of ceramic laminate with strongly bonded interfaces. Sih’s criterion (Eq. 3) was applied in each step of the simulations. The crack started to strongly deflect in certain depth a´ under the first ATZ/AMZ interface due to acting of compressive stresses in all simulations. In this moment an additional external load was necessary for the crack propagation. The depth of deflection and values of external force P are shown in the Tab. 3. tATZ : tAMZ a’ [mm] Pmax [N] 2 :1 0.046 16.2 5 : 1 0.020 24.8 7 : 1 0.017 34.6 10 : 1 0.024 48.0 Table 3: The depth a´ in AMZ layer under the first ATZ/AMZ interface, where the crack changes mode of propagation and values of the force P acting in the moment of crack deflection. Hundreds of calculations (steps) were performed in each simulation to obtain crack path in the AMZ layer. Results of simulations are shown in the Fig. 4. On the base of results obtained the crack behaviour in the compressive layer can be divided to four stages, see Fig. 5 for the explanation: 1) After passing perpendicularly through the ATZ/AMZ interface the crack is retarded (the stress intensity factor and the strain energy density factors decrease). An external load is necessary for further crack propagation. The compressive stresses don’t allow further direct crack propagation under mode I like in tensile loaded ATZ layer. 2) In the depth a´ under the interface the character (mode) of crack propagation changes from mode I to mode II. The stress distribution around the crack tip in the depth a´ enables crack deflection or bifurcation. Further crack propagation is controlled by compressive stresses. The crack propagates parallel or nearly parallel to the material interfaces. 3) When the crack tip is close to the next AMZ/ATZ interface the presence of tensile stresses in the ATZ layer growths in importance and the controlling mode of crack propagation starts to be again mode I and the crack deflect to the direction (nearly) perpendicular to AMZ/ATZ interface. 4) The crack passing through (nearly) perpendicularly the AMZ/ATZ interface. The resistance to the crack propagation is the highest between stages 3 and 4. Behaviour of the crack at the stage 4 was studied in [28]. D L. Náhlík et alii, Frattura ed Integrità Strutturale, 34 (2015) 116-124; DOI: 10.3221/IGF-ESIS.34.12 121 a) Result for tATZ : tAMZ = 2. Crack path in the AMZ layer after 270 steps (left), detail of crack deflection after 33 steps. b) Result for tATZ : tAMZ = 5. Crack path in the AMZ layer after 491 steps (left), detail of crack deflection. c) Result for tATZ : tAMZ = 7. Crack path in the AMZ layer after 315 steps (left), detail of crack deflection. d) Result for tATZ : tAMZ = 10. Crack path in the AMZ layer after 306 steps (left), detail of crack deflection. e) Detail of stress distribution in the laminates: tATZ : tAMZ = 2 (left), tATZ : tAMZ = 10. Figure 4: Simulation of crack propagation in the compressive layer for different ratio of layer thicknesses. Magnitude of longitudinal stress component is displayed. L. Náhlík et alii, Frattura ed Integrità Strutturale, 34 (2015) 116-124; DOI: 10.3221/IGF-ESIS.34.12 122 Figure 5: Scheme of stages of crack propagation in the compressive layer of ceramic laminate with strongly bonded interfaces. The crack can relatively easy propagate in the outside layers from some flaws, scratches or initial notches under external bending load of the composite. In the tensile layers the crack propagation corresponds to the mode I and cracks propagate perpendicularly to the material interface. Due to this fact it is necessary to prepare the composite with low tensile stresses, i.e. with wide tensile layer (ATZ in studied case) and thin compressive layer (AMZ) with high compressive stresses. The resistance to the crack propagation in the tensile layer is lower due to tensile residual stresses than in the case of crack propagation in the homogeneous ceramics, see e.g. [28]. When the crack passed through the interface to the compressive layer (AMZ) acting compressive residual stresses change the mode of crack propagation from mode I to shear mode II. This causes strong crack deflection or can lead to bifurcation of the propagating crack. Crack propagation parallel to the material interface in the compressive layer is evident from experiments published in [37]. This phenomenon is supported by results obtained as well. The parallel crack propagation is controlled by specific stress distribution in the compressive AMZ layer and influenced by the vicinity of tensile layers. After some propagation the bending nature of loading, existence of material imperfections (flaws, pores) and the vicinity of material interface lead to the other change of crack propagation direction close to AMZ/ATZ interface [32]. These effects contribute to the characteristic stepwise crack propagation in the strongly bonded ceramic laminates. CONCLUSIONS he paper presented focuses on the crack behaviour in the multilayered ceramic composite subjected to the bending load. Crack propagation in the compressive layer responsible for higher apparent fracture toughness of the composite was investigated by means of finite element method. Sih’s criterion based on strain energy density factor was used for estimation of the crack behaviour in the layers and for the estimation of crack path. It was shown that the crack after passing through the first material interface between tensile and compressive layer retards due to acting of strong residual stresses and the crack is sharply deflected. An additional external load is necessary for further crack propagation at this moment. The depth of crack deflection was determined for different thicknesses of the composite layers. The stress distribution allows in this moment bifurcation of the crack as well, like was shown e.g. in the works [29,31]. Further crack propagation in the compressive layer is more or less parallel to the material interface. The stages leading to the stepwise crack propagation were described in detail. Finite element method implemented in commercial system Ansys with author’s routines was used for numerical simulations. The paper contributes to the better understanding of damage of strongly bonded ceramic laminates and their toughening mechanism. ACKNOWLEDGEMENT his work was supported through the Grant No. 15-09347S of the Czech Science Foundation and the specific academic research grant (K. Štegnerová) No. FSI-S-14-2311 provided to Brno University of Technology, Faculty of Mechanical Engineering. T T L. Náhlík et alii, Frattura ed Integrità Strutturale, 34 (2015) 116-124; DOI: 10.3221/IGF-ESIS.34.12 123 REFERENCES [1] Kolednik, A., Predan, J., Gubeljak, N., Fischer, D.F., Modeling fatigue crack growth in a bimaterial specimen with the configurational forces concept, Mat. Sci. Eng. A, 519 (2009) 172–183. DOI: 10.1016/j.msea.2009.04.059. [2] Marsavina, L., Sadowski, T., Knec, M., Crack propagation paths in four point bend Aluminium–PMMA Specimens, Eng. Frac. Mech., 108 (2013) 139-151. DOI: 10.1016/j.engfracmech.2013.02.029. [3] Leguillon, D., Martin, E., The strengthening effect caused by an elastic contrast—part I: the bimaterial case, Int. J. Fract., 179 (2013) 157–167. DOI: 10.1007/s10704-012-9787-y. [4] He, M.Y., Hutchinson, J.W., Crack deflection at an interface between dissimilar elastic materials, Int. J. Solids Struct., 25 (1989) 1053–1067. [5] He, M.Y., Evans, A.G., Hutchinson, J.W., Crack deflection at an interface between dissimilar elastic materials: Role of residual stresses, Int. J. Solids Struct., 31 (1994) 3443-3455. DOI: 10.1016/0020-7683(94)90025-6. [6] Sistaninia, M., Kolednik, O., Effect of a single soft interlayer on the crack driving force, Eng. Frac. Mech., 130 (2014) 21–41. DOI: 10.1016/j.engfracmech.2014.02.026. [7] Knésl, Z., Náhlík, L., Radon, J., Influence of interface on fatigue threshold values in elastic biomaterials, Comp. Mater. Sci., 28 (2003) 620–627. DOI:10.1016/j.commatsci.2003.08.018. [8] Wang, B., Siegmund, T., Simulation of fatigue crack growth at plastically mismatched bi-material interfaces, Int. J. Plasticity, 22 (2006) 1586-1609. DOI: 10.1016/j.ijplas.2005.09.004. [9] Minatto, F.D., Milak, P., De Noni Jr., A., Hotza, D., Montedo, O. R. K., Multilayered ceramic composites – a review, Adv. Appl. Ceram., 114 (2015) 127-138. DOI: 10.1179/1743676114Y.0000000215. [10] Šestáková, L., Bermejo, R., Chlup, Z., Danzer, R., Strategies for fracture toughness, strength and reliability optimisation of ceramic–ceramic laminates, Int. J. Mat. Res., 102 (2011) 1-14. DOI: 10.3139/146.110523. [11] Bermejo, R., Chlup, Z., Šestáková, L., Ševeček, O., Danzer, R., Strategies to optimize the strength and fracture resistence of ceramic laminate, Mechanical Properties and Performance of Engineering Ceramics and Composites VII: Ceramic Engineering and Science Proceedings, 33 (2012) 163-174. [12] She, J., Inoue, T., Ueno, K., Multilayer Al2O3/SiC ceramics with improved mechanical behaviour, J. Eur. Ceram. Soc., 20 (2000) 1771-1775. DOI:10.1016/S0955-2219(00)00048-0. [13] Hadraba, H., Drdlík, D., Chlup, Z., Máca, K., Dlouhý, I., Cihlář, J., Layered ceramic composites via control of electrophoretic deposition kinetics, J. Eur. Ceram. Soc., 33 (2013) 2305–2312. DOI:10.1016/j.jeurceramsoc.2013.01.026. [14] Hadraba, H., Klimeš, J., Maca, K., Crack propagation in layered Al2O3/ZrO2 composites prepared by electrophoretic deposition, J. Mater. Sci., 42 (2007) 6404–6411. DOI: 10.1007/s10853-006-1197-y. [15] Oechsner, M., Hillman, C., Lange, F. F., Crack bifurcation in laminar ceramic composites, J. Am. Ceram. Soc., 79 (1996) 1834-1838. [16] Lugovy, M., Slyunyayev, V., Orlovskaya, N., Blugan, G., Kuebler, J., Lewis, M., Apparent fracture toughness of Si3N4-based laminates with residual compressive or tensile stresses in surface layers, Acta Mater., 53 (2005) 289–296. DOI:10.1016/j.actamat.2004.09.022. [17] Wei, P., Chen, L., Okubo, A., Hirai, T., Tough multilayered α–β SI3N4 ceramics prepared by spark plasma sintering, Mater. Lett., 49 (2001) 239–243. DOI: 10.1016/S0167-577X(00)00377-3. [18] Tomaszewski, H., Weglarz,H., Wajler, A., Boniecki, M., Kalinski, D., Multilayer ceramic composites with high failure resistence, J. Eur. Ceram. Soc., 27 (2007) 1373–1377. DOI:10.1016/j.jeurceramsoc.2006.04.030. [19] Chlup, Z., Hadraba, H., Slabáková, L., Drdlík, D., Dlouhý, I., Fracture behaviour of alumina and zirconia thin layered laminate, J. Eur. Ceram. Soc., 32 (2012) 2057–2061. DOI:10.1016/j.jeurceramsoc.2011.09.006. [20] Bermejo, R., Danzer, R., High failure resistance layered ceramics using crack bifurcation and interface delamination as reinforcement mechanisms, Eng. Frac. Mech., 77 (2010) 2126–2135. DOI: 10.1016/j.engfracmech.2010.02.020. [21] Bermejo, R., Llanes, L., Anglada, M., Supancic, P., Lube, T., Thermal Shock Behavior of an Al2O3/ZrO2 Multilayered Ceramic with Residual Stresses due to Phase Transformation, Key Engineering Materials, 290 (2005) 191-198. DOI: 10.4028/www.scientific.net/KEM.290.191. [22] Bermejo, R., Torres, Y., Baudín, C., Sánchez-Herencia, A.J., Pascual, J., Anglada, M., Llanes, L., Threshold strength evaluation on an Al2O3-ZrO2 multilayered system, J. Eur. Ceram. Soc., 27 (2007) 1443–1448. DOI:10.1016/j.jeurceramsoc.2006.05.037. [23] Bermejo, R., Baudín, C., Moreno, R., Llanes, L., Sánchez-Herencia, A.J., Processing optimisation and fracture behaviour of layered ceramic composites with highly compressive layers, Compos.Sci. Technol., 67 (2007) 1930–1938. DOI:10.1016/j.compscitech.2006.10.010. L. Náhlík et alii, Frattura ed Integrità Strutturale, 34 (2015) 116-124; DOI: 10.3221/IGF-ESIS.34.12 124 [24] Gotor, F.J., Bermejo, R., Córdoba, J.M., Chicardi, E., Medri, V., Dalle Fabbriche, D., Torres, Y., Processing and characterisation of cermet/hardmetal laminates with strong interfaces, Materials and Design, 58 (2014) 226–233. DOI: 10.1016/j.matdes.2014.01.076. [25] Kuo, D.H., Kriven, W.M., Fracture of multilayer oxide composites, Mat.Sci. Eng. A, 241 (1998) 241–250. [26] Zhou, P., Hu, P., Zhang, X., Han, W., Laminated ZrB2–SiC ceramic with improved strength and toughness, Scripta Mater., 64 (2011) 276–279. DOI:10.1016/j.scriptamat.2010.10.005. [27] Bermejo, R., Sánchez-Herencia, A.J., Llanes, L., Baudín, C., High-temperature mechanical behaviour of flaw tolerant alumina–zirconia multilayered ceramics, Acta Mater., 55 (2007) 4891–4901. DOI:10.1016/j.actamat.2007.05.005. [28] Náhlík, L., Šestáková, L., Hutař, P., Estimation of apparent fracture toughness of ceramic laminates, Comp. Mater. Sci., 46 (2009) 614–620. DOI:10.1016/j.commatsci.2009.04.005. [29] Ševeček, O., Bermejo, R., Kotoul, M., Prediction of the crack bifurcation in layered ceramics with high residual stresses, Eng. Frac. Mech., 108 (2013) 120–138. DOI: 10.1016/j.engfracmech.2013.03.013. [30] Kotoul, M., Ševeček, O., Vysloužil, T., Crack growth in ceramic laminates with strong interfaces and large compressive residual stresses, Theor. Appl. Frac.Mech., 61 (2012) 40–50. DOI: 10.1016/j.tafmec.2012.08.005. [31] Chen, C.R., Bermejo, R., Kolednik, O., Numerical analysis on special cracking phenomena of residual compressive inter-layers in ceramic laminates, Eng. Frac. Mech., 77 (2010) 2567–2576. DOI:10.1016/j.engfracmech.2010.06.020. [32] Náhlík, L., Šestáková, L., Hutař, P., Bermejo, R., Prediction of crack propagation in layered ceramics with strong interfaces, Eng. Frac. Mech., 77 (2010) 2192–2199. DOI:10.1016/j.engfracmech.2010.02.023. [33] Oel, H.J., Fréchette, V.D., Stress Distribution in Multiphase Systems: I, Composites with Planar Interfaces. J. Am. Ceram. Soc., 50 (1967) 542–549. DOI: 10.1111/j.1151-2916.1967.tb14992.x. [34] Sih, G.C., in: G.C. Sih (Ed.), A Special Theory of Crack Propagation: Methods of Analysis and Solutions of Crack Problems, Mechanics of Fracture, Noordhoff, International Publishing, Leyden, (1973) 21–45. [35] Berto, F., Lazzarin, P., Recent developments in brittle and quasi-brittle failure assessment of engineering materials by means of local approaches, Mat. Sci. Eng. R, 75 (2014) 1–48. DOI: 10.1016/j.mser.2013.11.001. [36] Ansys v14.5, Ansys, Inc., (2012). [37] Tarlazzi, A., Roncari, E., Pinasco, P., Guicciardi, S., Melandri, C., de Portu, G., Tribological behaviour of Al2O3/ZrO2-ZrO2 laminated composites, Wear, 244 (2000) 29–40. DOI:10.1016/S0043-1648(00)00429-4.