Elsevier

Engineering Fracture Mechanics

Volume 167, November 2016, Pages 56-67
Engineering Fracture Mechanics

A failure scenario of ceramic laminates with strong interfaces

https://doi.org/10.1016/j.engfracmech.2016.04.022Get rights and content

Highlights

  • Paper explains in detail scenario of damage mechanism of layered ceramic composite.

  • The reason for crack bifurcation in the compressive layer is explained.

  • The paper contributes to the better understanding of damage of ceramic laminates.

Abstract

Over the last years many researchers have put a lot of effort into designing layered structures combining different materials in order to improve low fracture toughness and mechanical reliability of ceramics. It has been proven that an effective way is to create layered ceramics with strongly bonded interfaces. Significant internal residual stresses are developed within the composite layers after the cooling process from the sintering temperature, due to the different coefficients of thermal expansion of individual composite constituents. Residual stresses can significantly change the crack behaviour. Suitable choice of material of layers and ratio between layer thicknesses can lead to higher value of the so-called apparent fracture toughness, i.e. higher resistance of the ceramic laminate to the crack propagation.

The paper deals with a description of the specific crack behaviour in the layered alumina–zirconia ceramics. Attention is devoted to the differences in the stress field description in the vicinity of the crack front. Two-dimensional and three-dimensional numerical models are developed for this purpose. 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 is investigated on the basis of linear elastic fracture mechanics. Fracture parameters are computed numerically and by routines of authors. The sharp change of the crack propagation direction is estimated using the Sih’s criterion based on the strain energy density factor, and conditions for crack bifurcation are determined. Estimated crack behaviour is qualitatively in a good agreement with experimental observations.

Introduction

Ceramic materials are widely used in many engineering applications for their temperature stability, hardness, mechanical strength or chemical inertness. However, ceramics are very brittle materials; therefore they are mainly utilized under compression loading. Over the last years many researchers have been working on improving of mechanical properties of ceramics by fabricating ceramic composites materials. It was shown that material interfaces play an important role in mechanical behaviour and material properties of composite, which can cause improving of fracture toughness of the ceramic laminate. The presence of multiple constituents influences the stress distribution in the composite. Therefore, the influence of different materials and geometry combination and weak and strong material interfaces on the resulting material behaviour were studied both experimentally and theoretically, see e.g. [1], [2], [3], [4], [5], [6], [7], [8], [9] for details.

Improvement of resistance against crack propagation through a ceramic composite body can be caused by several factors; (a) presence of residual stresses developed during the fabrication of ceramics due to different coefficients of thermal expansion of constituents and cooling down from high temperatures; (b) material interfaces; (c) consuming of more fracture energy during crack deflection or/and bifurcation [10], [11], [12], [13], [14], [15], see Fig. 1. Mentioned mechanisms were proved experimentally for different layered ceramic materials with strongly bonded interfaces, see [16], [17], [18], [19], [20], [21], [22], [23], [24], [25], [26], [27], [28], [29]. In [30] a procedure for estimation of apparent fracture toughness of ceramic laminate based on generalized linear elastic fracture mechanics was published. The specific crack behaviour in the compressive layer of layered ceramics composites was the subject of several studies [31], [32], [33], [34], where procedures, derived in the frame of linear elastic fracture mechanics, for estimation of crack deflection or bifurcation were published. However, comprehensive interpretation of failure of layered ceramic composite and explanation of the role of residual stresses is still missing in the literature.

As is mentioned above, the presence of residual stresses developed during sintering process in individual layers can significantly influence the fracture behaviour of the whole composite. Compressive residual stresses can cause closure of crack, its deflection from the original direction of crack propagation, or crack bifurcation. Furthermore a crack can stay arrested at the bimaterial interface under specific conditions. These mechanisms lead to consumption of more fracture energy, and contribute to increase the (apparent) fracture toughness of the composite [2], [3], [4], [7], [8], [9], [10], [11], [12], [13], [14], [15].

The motivation of this work is to investigate and describe the specific crack behaviour and failure mechanism of ceramic laminate with special focus on crack behaviour in the compressive layer of multilayer composites. For this purpose, numerical modelling and linear elastic fracture mechanics procedures are used.

A material system based on alumina and zirconia is herein chosen. These two constituents can be taken as typical examples of convenient materials for layered ceramics. Behaviour of such kinds of laminates was experimentally studied e.g. in [20], [21], [22], [23]. Considered ceramic laminate is created by 9 layers; 5 are made of Al2O3/5vol.%t-ZrO2 (alumina with tetragonal zirconia, referred as ATZ) and 4 are made of Al2O3/30vol.%m-ZrO2 (alumina with monoclinic zirconia, referred as AMZ). Material characteristics are taken from the literature [20], [21] and summarized in Table 1. Table 2 shows grain sizes of the laminate components. The laminate is subjected to four-point bending (4PB) test, see Fig. 2.

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 the following expressions [35]:σres,ATZ=TsfT0(αAMZ-αATZ)dT1EATZ+1EAMZ·ρ·N+1N-1,σres,AMZ=-σres,ATZ·ρ·N+1N-1,where ρ=tATZtAMZ, EAMZ=EAMZ1-νAMZ, EATZ=EATZ1-νATZ, Tsf [°C] is a stress free temperature, T0 [°C] room temperature, α [10−6 K−1] the coefficient of thermal expansion, N [–] the number of layers, t [mm] the thickness of layer, E [MPa] Young’s modulus and ν [–] Poisson’s ratio. Calculated values of residual stresses for the studied configurations are shown in Table 3 and are in good agreement with experimental observation [36]. Stress free temperature Tsf = 1200 °C is used in the calculations.

From Table 3, it can be noted that the compressive residual stresses in the AMZ layers are about six hundred megapascals and higher for all considered cases. These stresses are responsible for typical stepwise crack propagation, as is shown in the scheme in Fig. 2.

Section snippets

Numerical model of crack propagation in tensile layer

In the previous studies [30], [34]D numerical models were used to determine the crack behaviour in the kind of laminate being examined. In the following, the influence of real crack front shape on crack behaviour is analyzed and 3D numerical models are used. The aim is to determine differences in the stress fields between results obtained by the use of 2D and 3D numerical models and to investigate whether the modelling of natural crack front curvature can significantly influence results

Crack touching the interface between two layers

When the crack front touches the interface between two different materials, the stress singularity exponent changes from 0.5 (characteristic value for homogeneous materials) to a generic value ranging in the interval (0;1). This change depends on the elastic properties of both materials. This well-known fact has been deeply studied in the past, see e.g. [4], [5], [7], [41], [42], [43], [44], [45] for details. Therefore, authors focus only on the specific differences between 2D and 3D solutions.

Crack propagation in compressive layer

In the following, a considerable effort is made to apply the procedure described in the section devoted to the crack propagation in the tensile layer to the case of crack propagation just behind the material interface. However, semi-elliptical stable crack front shape is not found by this procedure. The last stable crack front shape is found in the case of the crack touching the material interface, see configuration in Fig. 4c). It means that, for crack propagation in the compressive (AMZ)

Results: Description of damage mechanism of ceramic laminate

On the basis of the results obtained in previous sections, the crack behaviour in the compressive layer can be divided into five stages, see Fig. 11 for explanation:

  • 1.

    The crack is retarded after passing through the ATZ/AMZ interface and perpendicular to the material interface (the stress intensity factor and the strain energy density factors decrease). An external load is necessary for further crack propagation. The compressive stresses do not allow further direct crack propagation under mode I

Conclusions

The present paper focuses on the crack behaviour in the multilayered ceramic composite subjected to the bending load. Crack behaviour ahead of ATZ/AMZ interface, at the interface and behind the interface is studied. Special attention is devoted to the differences in the stress fields obtained from 2D and 3D FE numerical models and to the crack propagation in the compressive layer. The crack propagation in the compressive layer responsible for higher apparent fracture toughness of the composite

Acknowledgements

The work has been supported through Grant No. 15-09347S of the Czech Science Foundation.

The equipment and the base of research infrastructure IPMINFRA and CEITEC – Central European Institute of Technology (LQ 1602) have been used during the research activities.

References (52)

  • Z. Chlup et al.

    Fracture behaviour of alumina and zirconia thin layered laminate

    J Eur Ceram Soc

    (2012)
  • R. Bermejo et al.

    High failure resistance layered ceramics using crack bifurcation and interface delamination as reinforcement mechanisms

    Engng Fract Mech

    (2010)
  • R. Bermejo et al.

    Threshold strength evaluation on an Al2O3–ZrO2 multilayered system

    J Eur Ceram Soc

    (2007)
  • R. Bermejo et al.

    Processing optimisation and fracture behaviour of layered ceramic composites with highly compressive layers

    Compos Sci Technol

    (2007)
  • F.J. Gotor et al.

    Processing and characterisation of cermet/hardmetal laminates with strong interfaces

    Mater Des

    (2014)
  • D.H. Kuo et al.

    Fracture of multilayer oxide composites

    Mater Sci Engng A

    (1998)
  • P. Zhou et al.

    Laminated ZrB2–SiC ceramic with improved strength and toughness

    Scripta Mater

    (2011)
  • R. Bermejo et al.

    High-temperature mechanical behaviour of flaw tolerant alumina–zirconia multilayered ceramics

    Acta Mater

    (2007)
  • L. Náhlík et al.

    Estimation of apparent fracture toughness of ceramic laminates

    Comput Mater Sci

    (2009)
  • O. Ševeček et al.

    Prediction of the crack bifurcation in layered ceramics with high residual stresses

    Engng Fract Mech

    (2013)
  • M. Kotoul et al.

    Crack growth in ceramic laminates with strong interfaces and large compressive residual stresses

    Theor Appl Fract Mech

    (2012)
  • C.R. Chen et al.

    Numerical analysis on special cracking phenomena of residual compressive inter-layers in ceramic laminates

    Engng Fract Mech

    (2010)
  • L. Náhlík et al.

    Prediction of crack propagation in layered ceramics with strong interfaces

    Engng Fract Mech

    (2010)
  • M. Ševčík et al.

    Estimation of the critical configuration of a crack arrested at the interface between two materials

    Comput Mater Sci

    (2012)
  • P. Hutař et al.

    Multilayer polymer pipes failure assessment based on a fracture mechanics approach

    Engng Fail Anal

    (2013)
  • M. Ševčík et al.

    Numerical estimation of the fatigue crack front shape for a specimen with finite thickness

    Int J Fatigue

    (2012)
  • Cited by (7)

    • Effect of residual stresses to the crack path in alumina/zirconia laminates

      2020, Journal of the European Ceramic Society
      Citation Excerpt :

      In the case of indentation cracks, the direction of propagation is from the nearly stress-free surface to the interior of laminate where residual stresses are high, i.e. only elastic mismatch is controlling the trajectory at the beginning. In the case of notched bending bars, the crack front is usually curved as well-known from the literature [48–50]. In the central part, plane strain conditions increase stress triaxiality which enhances brittle fracture and thus crack propagates faster in the interior than at the surface with plane stress region [51].

    • Effects of micro and macro features on the thermal shock resistance of laminated Al<inf>2</inf>O<inf>3</inf>–MgO composites

      2020, Ceramics International
      Citation Excerpt :

      However, the length of stage I for the ATMC composite was obviously shorter than that of the LAMC-1 and LAMC-2 composites, and it can almost be ignored. It was easily observed that LAMC-1 and LAMC-2 composites exhibited a slight ductile characteristic in Fig. 5a, which is similar to results in other studies [16,21,22]. Rapid transformation of the displacement occurred as the applied stress exceeded a critical value, with a brittle characteristic of the rupture process appearing (stage II).

    • Fabrication, mechanical properties and fracture behaviors of the laminated Al<inf>2</inf>O<inf>3</inf>–ZrB<inf>2</inf>–MgO / Al<inf>2</inf>O<inf>3</inf>–TiN–MgO ceramic composite

      2020, Ceramics International
      Citation Excerpt :

      The relation between the residual flexural strength, up to σA in formula (5), of ATM 4 sample and indentation force in the front and side directions are shown in Fig. 6. ψ and ξ were estimated to be 1.24 and 0.016 according to the literature [29]. The flexural strengths under different Vickers indentations were lower than the original flexural strengths in the front and side directions.

    View all citing articles on Scopus
    View full text