Availableonlineatwww.sciencedirect.com SCIENCE DIRECT E噩≈S ELSEVIER Journal of the European Ceramic Society 26(2006)343-349 www.elsevier.com/locate/jeurceran Prediction of crack deflection in porous/dense ceramic laminates D. Leguillona*,S.Tariolleb, E. Martin, T. Chartier, J.L. Besson a LMM, CNRS UMR7607, Universite P et M. Curie, case 162, 4 place Jussieu. 75252 Paris Cedex 05, france b CES CNRS UMR 5146. Ecole des mines de saint-Etienne. St Etienne. france c LCTS. CNRS UMR 5801 Universite bordeaux Pessac france d sPCTS, CNRS UMR 6638, E. N.S. de Ceramique Industrielles, Limoges, france Received 24 June 2004; received in revised form 4 November 2004; accepted 12 November 2004 Abstract The arrangement of ceramic layers in laminated structures is an interesting way to enhance the flaw tolerance of brittle ceramic materials The interfaces are expected to deflect cracks, increasing the fracture energy of the laminate compared to a monolithic material and thus raising Laminates have been fabricated with alternating dense and porous layers of the same material, i.e. SiC or B. C, in order to obtain a good chemical compatibility between the laminas and almost no thermal residual stresses. Porosity, in the porous layers, is achieved by incorporating organic particles which are removed during the debinding step. In this context, the target of this study is to predict the volume fraction of pores, in the porous layer, required to cause crack deflection. The proposed criterion derives from an energy balance. It relies on a two-scale analysis taking into account the laminated structure of the material. It can be written in terms of two relevant material parameters: the ratio of Youngs moduli of the dense and porous materials and toughness ratios. a unique function depending on the volume fraction of pores can be used to express the two above-mentioned ratios Assuming a cubic lattice of spherical voids, the parameters of the porous ceramic depend linearly on the porosity and vanish at the point of percolation of pores. As a consequence, the criterion can be rewritten in term of a single parameter: the porosity. Crack deflection is permitted only for very high values of the porosity. Predicted values agree satisfactorily with experiments on Sic and B.C. The comparison with the He and hutchinson criterion shows that this latter underestimates the correct value C 2004 Elsevier Ltd. All rights reserved Keywords: Fracture; Toughness and toughening: Laminates; Porosity 1. Introduction ally fabricated by stacking layers of different compositions in a suitable sequence. Alternating dense and porous layer Laminated structures with weak interfaces or interphases of the same material offers the best chemical compatibility constitute one strategy to improve the flaw tolerance of brit- and almost no thermal residual stresses. Crack deflection ha tle ceramic materials. These structures have proved to be ef- been observed within such systems as the result of the pres- ficient in increasing fracture energy by promoting crack de- ence of porous interlayers. - The question is thus to predict flection mechanisms at interfaces between layers: cracks that the volume fraction of pores, in the porous layers, required form in one layer are deflected along the interface with adja- to cause crack deflection nt layers which increases the fracture energy of the laminate Two kinds of ceramic laminates are analysed silicon car- compared to that of a monolithic material and thus raises the bide(Sic)2-7 and boron carbide(B4C).7, 8 The specimens apparent toughness. Laminar ceramic composites are gener- are fabricated by stacking layers obtained by tape-casting and then by lamination, debinding and sintering. Porosity of Corresponding author. Tel. +33 144 275 322; fax: +33 144 275 259 porous layers is introduced by adding pore-forming agents E-mail address: dol(@ccr, jussieu. fr(D. Leguillon) such as corn starch or polymer particles. These organic par- 0955-2219/S-see front matter c 2004 Elsevier Ltd. All rights reserved doi: 10.1016/j. jeurceramsoc 2004.1 1.003
Journal of the European Ceramic Society 26 (2006) 343–349 Prediction of crack deflection in porous/dense ceramic laminates D. Leguillona,∗, S. Tariolleb, E. Martinc, T. Chartierd, J.L. Bessond a LMM, CNRS UMR7607, Universit´e P. et M. Curie, case 162, 4 place Jussieu, 75252 Paris Cedex 05, France b CES, CNRS UMR 5146, Ecole des Mines de Saint-Etienne, St Etienne, France c LCTS, CNRS UMR 5801, Universit´e Bordeaux 1, Pessac, France d SPCTS, CNRS UMR 6638, E.N.S. de C´eramiques Industrielles, Limoges, France Received 24 June 2004; received in revised form 4 November 2004; accepted 12 November 2004 Available online 19 January 2005 Abstract The arrangement of ceramic layers in laminated structures is an interesting way to enhance the flaw tolerance of brittle ceramic materials. The interfaces are expected to deflect cracks, increasing the fracture energy of the laminate compared to a monolithic material and thus raising the toughness. Laminates have been fabricated with alternating dense and porous layers of the same material, i.e. SiC or B4C, in order to obtain a good chemical compatibility between the laminas and almost no thermal residual stresses. Porosity, in the porous layers, is achieved by incorporating organic particles which are removed during the debinding step. In this context, the target of this study is to predict the volume fraction of pores, in the porous layer, required to cause crack deflection. The proposed criterion derives from an energy balance. It relies on a two-scale analysis taking into account the laminated structure of the material. It can be written in terms of two relevant material parameters: the ratio of Young’s moduli of the dense and porous materials and toughness ratios. A unique function depending on the volume fraction of pores can be used to express the two above-mentioned ratios. Assuming a cubic lattice of spherical voids, the parameters of the porous ceramic depend linearly on the porosity and vanish at the point of percolation of pores. As a consequence, the criterion can be rewritten in term of a single parameter: the porosity. Crack deflection is permitted only for very high values of the porosity. Predicted values agree satisfactorily with experiments on SiC and B4C. The comparison with the He and Hutchinson criterion shows that this latter underestimates the correct value. © 2004 Elsevier Ltd. All rights reserved. Keywords: Fracture; Toughness and toughening; Laminates; Porosity 1. Introduction Laminated structures with weak interfaces or interphases constitute one strategy to improve the flaw tolerance of brittle ceramic materials. These structures have proved to be ef- ficient in increasing fracture energy by promoting crack de- flection mechanisms at interfaces between layers: cracks that form in one layer are deflected along the interface with adjacent layers which increases the fracture energy of the laminate compared to that of a monolithic material and thus raises the apparent toughness. Laminar ceramic composites are gener- ∗ Corresponding author. Tel.: +33 144 275 322; fax: +33 144 275 259. E-mail address: dol@ccr.jussieu.fr (D. Leguillon). ally fabricated by stacking layers of different compositions in a suitable sequence.1 Alternating dense and porous layers of the same material offers the best chemical compatibility and almost no thermal residual stresses. Crack deflection has been observed within such systems as the result of the presence of porous interlayers.2–9 The question is thus to predict the volume fraction of pores, in the porous layers, required to cause crack deflection. Two kinds of ceramic laminates are analysed: silicon carbide (SiC)2–7 and boron carbide (B4C).7,8 The specimens are fabricated by stacking layers obtained by tape-casting, and then by lamination, debinding and sintering. Porosity of porous layers is introduced by adding pore-forming agents such as corn starch or polymer particles. These organic par- 0955-2219/$ – see front matter © 2004 Elsevier Ltd. All rights reserved. doi:10.1016/j.jeurceramsoc.2004.11.003
guillon et al. Journal of the European Ceramic Sociery 26(2006)343-349 Fig. 1. Macrostructure of boron carbide composite with interlayers obtained by the use of corn starch(left), microstructure of porous boron carbide obtained with 20 vol. of corn starch(right) ticles are burned out during the debinding step, prior to sin- The present model relies on a two-scale analysis where ering Dense and porous layers have the same thickness af- the laminated structure is taken into account as explained in ter sintering(e a 100 um). Laminated specimens made of 20 Section 2. The criterion proposed in Section 3 derives from an layers have been tested under 3-point flexure loading. De- energy balance and can take into account a complementary tails on the fabrication of specimens and the measurement stress criterion. This makes it possible to avoid the above- of elastic and failure parameters can be found in Reynaud mentioned drawback, since virtual crack extension lengths and co-workers,6 and in Tariolle et al. 7,8 Fig. 1 shows mi- are known. The criterion is expressed in terms of two relevant shapes within the porous interlayer, and depicts the pore parameters: the Young s moduli and toughness ratios of the crographs of the layered B4C material dense and porous layers In Section 4, experiments show that The analyses of crack deflections by interfaces are gen- these ratios can be written quite simply as functions of the erally based on two models due to He and Hutchinson porosity. Sections 5 and 6 are devoted to the study of the crack Both are carried out in an unbounded domain made of two deflection at the dense/porous and porous/dense interfaces. A elastic materials. In the first one, the primary crack lies in short Section 7 is dedicated to the analysis of the role of the one material and impinges on the interface. Two virtual crack relative thickness of porous layers. Section 8 deals with the extensions are considered, either deflected along the interface minor influence of Poissons ratios (Fig. 2b)or straight in the adjacent layer. The energy release rates at the tip of these two extensions are compared. The 2. The asymptotics of the problem drawback of this approach is the arbitrary choice of the two increment lengths. In the second He and Hutchinson model, I The model is based on a two-scale analysis, the small pa- the primary crack is along the interface and the ability of the rameter being the layers thickness. At the macro scale the crack to kink out of the interface is studied. The two criteria laminated micro-structure is ignored, as a first approxima- involve the toughness of the materials and of the interface. tion, the whole laminate is treated as a homogeneous mate- Curiously, it is often this second paper that is used to inter- rial. It is homogenized using a rule of mixture for simplicity pret the experimental results of cracks deflection in ceramic since more sophisticated homogenization processes do not laminates, although the main assumption, a long primary in- bring significant differences in the final results. 12 There is a terface crack, is not fulfilled. Moreover, in any case, it is clear primary crack for which the tip undergoes the classical mode from Fig. 2b that the laminated environment of the crack tip I singular field. The antisymmetric mode II is inhibited due is ignored in these approaches to the symmetries Within this framework. the dis placement solution(so- called far field) prior to any crack growth can be written in ED °(x1,x2)={(0,0)+kr√F1(6)+ (1) Here xi and x2 stand for the Cartesian coordinates and rand 8 for the polar ones. The coefficient ky is the stress intensity factor and the angular shape function is denoted u(0 Considering now a small crack extension e, the perturbed solution is expressed as a correction brought to the initial term(Eq ( 1)) Fig 2 Schematic view of the inner(unbounded)stretched domain for a U((x1, x2)=U(1, x2)+small correction deflection at a porous/dense interface: (a)the present analysis(P: porous, D dense, H: homogenized); (b) He and Hutchinson approach The small correction is assumed to vanish as e-0
344 D. Leguillon et al. / Journal of the European Ceramic Society 26 (2006) 343–349 Fig. 1. Macrostructure of boron carbide composite with interlayers obtained by the use of corn starch (left), microstructure of porous boron carbide obtained with 20 vol.% of corn starch (right). ticles are burned out during the debinding step, prior to sintering. Dense and porous layers have the same thickness after sintering (e ≈ 100m). Laminated specimens made of 20 layers have been tested under 3-point flexure loading. Details on the fabrication of specimens and the measurement of elastic and failure parameters can be found in Reynaud and co-workers5,6 and in Tariolle et al.7,8 Fig. 1 shows micrographs of the layered B4C materials and depicts the pore shapes within the porous interlayers.8 The analyses of crack deflections by interfaces are generally based on two models due to He and Hutchinson.10,11 Both are carried out in an unbounded domain made of two elastic materials. In the first one,10 the primary crack lies in one material and impinges on the interface. Two virtual crack extensions are considered, either deflected along the interface (Fig. 2b) or straight in the adjacent layer. The energy release rates at the tip of these two extensions are compared. The drawback of this approach is the arbitrary choice of the two increment lengths. In the second He and Hutchinson model,11 the primary crack is along the interface and the ability of the crack to kink out of the interface is studied. The two criteria involve the toughness of the materials and of the interface. Curiously, it is often this second paper that is used to interpret the experimental results of cracks deflection in ceramic laminates, although the main assumption, a long primary interface crack, is not fulfilled. Moreover, in any case, it is clear from Fig. 2b that the laminated environment of the crack tip is ignored in these approaches. Fig. 2. Schematic view of the inner (unbounded) stretched domain for a deflection at a porous/dense interface: (a) the present analysis (P: porous, D: dense, H: homogenized); (b) He and Hutchinson approach.10 The present model relies on a two-scale analysis where the laminated structure is taken into account as explained in Section 2. The criterion proposed in Section 3 derives from an energy balance and can take into account a complementary stress criterion. This makes it possible to avoid the abovementioned drawback, since virtual crack extension lengths are known. The criterion is expressed in terms of two relevant parameters: the Young’s moduli and toughness ratios of the dense and porous layers. In Section 4, experiments show that these ratios can be written quite simply as functions of the porosity. Sections 5 and 6 are devoted to the study of the crack deflection at the dense/porous and porous/dense interfaces. A short Section 7 is dedicated to the analysis of the role of the relative thickness of porous layers. Section 8 deals with the minor influence of Poisson’s ratios. 2. The asymptotics of the problem The model is based on a two-scale analysis, the small parameter being the layers thickness. At the macro scale the laminated micro-structure is ignored, as a first approximation, the whole laminate is treated as a homogeneous material. It is homogenized using a rule of mixture for simplicity, since more sophisticated homogenization processes do not bring significant differences in the final results.12 There is a primary crack for which the tip undergoes the classical mode I singular field. The antisymmetric mode II is inhibited due to the symmetries. Within this framework, the displacement solution (socalled far field) prior to any crack growth can be written in plane elasticity: U- 0(x1, x2) = U- 0(0, 0) + kI √r u- I (θ) +··· (1) Here x1 and x2 stand for the Cartesian coordinates and r and θ for the polar ones. The coefficient kI is the stress intensity factor and the angular shape function is denoted u- I (θ). Considering now a small crack extension , the perturbed solution is expressed as a correction brought to the initial term (Eq. (1)): U- (x1, x2) = U- 0(x1, x2) + small correction (2) The small correction is assumed to vanish as → 0.
D. Leguillon et al / Journal of the European Ceramic Society 26(2006)343-349 and porous layers and in the homogenized remaining part. The third equation expresses that the crack faces are free of traction. Finally, the last one is the matching condition with he mode I term involved in the far field(Eq (1)) Similarly, a crack extension e either in the next layer or along the interface(a deflection at the interface porous/dense is illustrated in Fig. 2a)leads to the following expansion =U0,0)+k√eW(y1,y2,)+ where u= e/e is the dimensionless crack extension length Here w must fulfil the same system of equations(Eq (4)), the traction free condition(Eq. (43)) being extended to the faces of the new extension of the crack 3. The deflection criterion Within this framework, the leading term of the change in Fig 3. The inner stretched domain and its artificial outer boundary r. The potential energy between the two states(prior to and follow- origin is located at the tip of the primary crack (i.e. before branching) ing a crack extension)is written: 13 The micro scale is obtained by stretching the domain 8W=ki/A/Ep E around the primary crack tip by 1/e, e being the layer thick- ness. Considering the limit e-0, the problem is now set- where d stands for the specimen depth(plane elasticity ). Ed tled in an unbounded domain, so-called inner domain. In or- and Ep are respectively the Youngs moduli of the dense and der to have tractable computations, this domain is artificially the porous ceramics. Poisson's ratios play a minor role as bounded at a large distance(roo>>1, where roo is the radius shown below in Fig 12. The function A is numerically derived of the artificial boundary and I the dimensionless stretched from the displacement field W using a contour integral: 13, thickness of the layers, say r=200 for instance, see Fig 3) from the primary crack tip. Moreover, only few dense(D v(W(y1,y,p),√pu(6) ig. 2a)and porous(P in Fig. 2a) layers, say 3 or 4, are kept being replaced by the homogenized material (H in Fig et with The validity of this simplification will be discussed in Section v(y, 5)=2/(o(U)xy-o(y)ny)ds Using the change of variable y =xi le(p=rle), Uo can be I is any contour in the inner domain surrounding and located far from the crack tip and its extension, n is its normal pointing U(1, x2)=y(ey, ey2) toward the crack tip. For practical reasons, the outer artificial boundary ro is selected( see Fig 3). The integral y in Ec =(0,0)+kr√ew(y1,y2,0)+ (3) (7)is contour independent for any and y fulfilling th equilibrium n(Eq. 41) where the 0 in W recalls that in a first step there is no crack a necessary condition for the crack growth is a conse extension. The function w is the solution to the follo quence of an energy balance 8W≥c→kf(p Vy·a=0 g=C: V,w g. n=along the crack faces E A(Ep/Ed, p)-A(Ep/Ed, O) The first equation is the balance of momentum(equilibrium). Here G is the toughness in the direction of fracture and ed The symbol nabla Vy holds for derivatives with respect to is the newly created crack surface. This expression(Eq (8)) yI and y2. The second equation is the constitutive law, C must be considered twice, once for a deflection(index def in is the elastic operator, it takes different values in the dense e following) and once for a penetration in the next layer
D. Leguillon et al. / Journal of the European Ceramic Society 26 (2006) 343–349 345 Fig. 3. The inner stretched domain and its artificial outer boundary Γ ∞. The origin is located at the tip of the primary crack (i.e. before branching). The micro scale is obtained by stretching the domain around the primary crack tip by 1/e, e being the layer thickness. Considering the limit e→0, the problem is now settled in an unbounded domain, so-called inner domain. In order to have tractable computations, this domain is artificially bounded at a large distance (r∞ 1, where r∞ is the radius of the artificial boundary and 1 the dimensionless stretched thickness of the layers, say r∞ = 200 for instance, see Fig. 3) from the primary crack tip. Moreover, only few dense (D in Fig. 2a) and porous (P in Fig. 2a) layers, say 3 or 4, are kept in the vicinity of the primary crack tip, the remaining part being replaced by the homogenized material (H in Fig. 2a). The validity of this simplification will be discussed in Section 3. Using the change of variable yi = xi/e (ρ = r/e), U- 0 can be expanded as (near field): U- 0(x1, x2) = U- 0(ey1, ey2) = U- 0(0, 0) + kI √e W- (y1, y2, 0) +··· (3) where the 0 in W- recalls that in a first step there is no crack extension. The function W- is the solution to the following problem: −∇y · σ = 0 σ = C : ∇yW- σ · n- = 0 along the crack faces W- behaves like √ρ u- I (θ) at infinity (4) The first equation is the balance of momentum (equilibrium). The symbol nabla y holds for derivatives with respect to y1 and y2. The second equation is the constitutive law, C is the elastic operator, it takes different values in the dense and porous layers and in the homogenized remaining part. The third equation expresses that the crack faces are free of traction. Finally, the last one is the matching condition with the mode I term involved in the far field (Eq. (1)). Similarly, a crack extension either in the next layer or along the interface (a deflection at the interface porous/dense is illustrated in Fig. 2a) leads to the following expansion: U- (x1, x2) = Ueµ(ey1, ey2) = U- 0(0, 0) + kI √e W- (y1, y2, µ) +··· (5) where µ = /e is the dimensionless crack extension length. Here W- must fulfil the same system of equations (Eq. (4)), the traction free condition (Eq. (43)) being extended to the faces of the new extension of the crack. 3. The deflection criterion Within this framework, the leading term of the change in potential energy between the two states (prior to and following a crack extension) is written:13 W = k2 I A Ep Ed , µ − A Ep Ed , 0 ed (6) where d stands for the specimen depth (plane elasticity). Ed and Ep are respectively the Young’s moduli of the dense and the porous ceramics. Poisson’s ratios play a minor role as shown below in Fig. 12. The functionAis numerically derived from the displacement field W- using a contour integral:13,14 A Ep Ed , µ = ψ(W- (y1, y2, µ), √ρ u- I (θ)) (7) with ψ(U- , V- ) = 1 2 Γ (σ(U- )n-V- − σ(V- )n-U- ) dS Γ is any contour in the inner domain surrounding and located far from the crack tip and its extension, n- is its normal pointing toward the crack tip. For practical reasons, the outer artificial boundary Γ ∞ is selected (see Fig. 3). The integral ψ in Eq. (7) is contour independent for any U- and V- fulfilling the equilibrium equation (Eq. 41). A necessary condition for the crack growth is a consequence of an energy balance: W ≥ Gc d ⇒ k2 I f Ep Ed , µ ≥ Gc (8) with f Ep Ed , µ = A(Ep/Ed, µ) − A(Ep/Ed, 0) µ Here Gc is the toughness in the direction of fracture and d is the newly created crack surface. This expression (Eq. (8)) must be considered twice, once for a deflection (index def in the following) and once for a penetration in the next layer
D. Leguillon et al. /Journal of the European Ceramic Sociery 26(2006)343-349 (index pen). Deflect promoted if the above inequality holds true for deflection but is wrong for penetration, it lead Ep/Ed fpen(Ep/Ed 0.6 where Gpen and Gder are the toughness of the next layer(pen- 04 etration mechanism)and of the interface(deflection mecha nism). It will be assumed in the following that the toughness of the interface is that of the porous material. Indeed, if the interface was stronger then the crack would grow within the porous medium at a short distance from the interface. Such a choice is also suggested in Fujita et al. I5 03 Clearly the dimensionless crack increment lengths ud Fig. 5. Youngs moduli ratio vs pore volume fraction V: SiC with Polyamide and Apen play a role in the above relation. Now we make the particles -7(diamonds), SiC with com starch particles, 7, (squares),BaC following reasonable additional assumption: if the crack pen with corn starch particles(triangles). Shear modulus ratio vs pore volume Apen=1. The deflected extension length remains to be deter mined. It could be done using a maximum stress criterion. 3 For simplicity, we assume he toughness parameters. They can be expressed in terms of the porosity as proposed in the next section hH11mt么u udef=Ppen=I (10) A simplification of the geometry of the since it refers to a characteristic length of the microstructure part of the material being replaced by a homogenized(aver (the layer thickness) that does not exist in their approach aged )one. Fig. 4 compares the function g defined in Eq(12) (Fig 2a and b). Nevertheless, complete computations relying when 3 and 4 layers are kept around the crack tip. Obviously, on a stress criterion have been carried out and it has been the simplified geometry retained in this model seems to be but that the final results are not deeply modified l orosity observed that the deflection length increases with the sufficient to our purpose Finally the criterion takes the simplified form: Ep 4. Elastic and fracture parameters of the porous with The two Figs. 5 and 6 show that a unique function depend E Adet(ep/Ed, D)-A(Ep/E ng on the volume fraction of pores V can be used to express (12) the elastic and fracture parameters of the porous material. In Ed/ Apen( Ep/Ed, 1)-A(Ep/ the first case(Eq (13)the parameters depend linearly on the It is clear from this expression that the crucial point is the volume fraction of pores V, whereas they depend linearly on knowledge of two relevant data: the ratios of the elastic and 02 g 4 layers 0.2 PreD Fig. 6. Toughness ratio vs. pore volume fraction V: Fig 4. The function g(Eq(12)(solid lines)vs the Young's moduli ratio particles -7(diamonds ) SiC with com starch particle quare), B4C Ep/Ed for two models of geometry(3 and 4 layers)at the porous/de with corn starch particles7, (circles), SiC with PTFE (triangles). H(V (Eq (13))dashed line, K(V(Eq (14) solid line
346 D. Leguillon et al. / Journal of the European Ceramic Society 26 (2006) 343–349 (index pen). Deflection is promoted if the above inequality holds true for deflection but is wrong for penetration, it leads to: fdef(Ep/Ed, µdef) fpen(Ep/Ed, µpen) ≥ Gc def Gc pen (9) where Gc pen and Gc def are the toughness of the next layer (penetration mechanism) and of the interface (deflection mechanism). It will be assumed in the following that the toughness of the interface is that of the porous material. Indeed, if the interface was stronger then the crack would grow within the porous medium at a short distance from the interface. Such a choice is also suggested in Fujita et al.15 Clearly the dimensionless crack increment lengths µdef and µpen play a role in the above relation. Now we make the following reasonable additional assumption: if the crack penetrates the next layer then it breaks it completely: pen = e ⇒ µpen = 1. The deflected extension length remains to be determined. It could be done using a maximum stress criterion.13 For simplicity, we assume here: µdef = µpen = 1 (10) This choice is less arbitrary than that of He and Hutchinson11 since it refers to a characteristic length of the microstructure (the layer thickness) that does not exist in their approach (Fig. 2a and b). Nevertheless, complete computations relying on a stress criterion have been carried out and it has been observed that the deflection length increases with the porosity but that the final results are not deeply modified.12,16 Finally the criterion takes the simplified form: g Ep Ed ≥ Gc def Gc pen (11) with g Ep Ed = Adef(Ep/Ed, 1) − A(Ep/Ed, 0) Apen(Ep/Ed, 1) − A(Ep/Ed, 0) (12) It is clear from this expression that the crucial point is the knowledge of two relevant data: the ratios of the elastic and Fig. 4. The function g (Eq. (12)) (solid lines) vs. the Young’s moduli ratio Ep/Ed for two models of geometry (3 and 4 layers) at the porous/dense interface. Fig. 5. Young’s moduli ratio vs. pore volume fraction V: SiC with Polyamide particles5–7 (diamonds), SiC with corn starch particles3,7,8 (squares), B4C with corn starch particles7,8 (triangles). Shear modulus ratio vs. pore volume fraction V: B4C with corn starch particles8 (circles). The dashed line is the function H(V) (Eq. (13)), the solid line is the function K(V) (Eq. (14)). toughness parameters. They can be expressed in terms of the porosity as proposed in the next section. A simplification of the geometry of the inner domain (Figs. 2a and 3) has been used. Only three or four layers have been considered around the crack tip, the remaining part of the material being replaced by a homogenized (averaged) one. Fig. 4 compares the function g defined in Eq. (12) when 3 and 4 layers are kept around the crack tip. Obviously, the simplified geometry retained in this model seems to be sufficient to our purpose. 4. Elastic and fracture parameters of the porous ceramic The two Figs. 5 and 6 show that a unique function depending on the volume fraction of pores V can be used to express the elastic and fracture parameters of the porous material. In the first case (Eq. (13)) the parameters depend linearly on the volume fraction of pores V, whereas they depend linearly on Fig. 6. Toughness ratio vs. pore volume fraction V: SiC with polyamide particles5–7 (diamonds), SiC with corn starch particles3,7,8 (squares), B4C with corn starch particles7,8 (circles), SiC with PTFE particles3 (triangles). H(V) (Eq. (13)) dashed line, K(V) (Eq. (14)) solid line.
D. Leguillon et al / Journal of the European Ceramic Society 26(2006)343-349 the largest surface fraction of pores S in the second case(Eq (14) Ep= H(Ed; Gp=H(V)Gd with (13) Ep=k(Ed; Gp=k(vGa with 4S K(V)=1-= (14) where Gp and Ga are the toughness respectively of the porous and the dense ceramics Fig. 7. The function g(Eq (12))(solid line)vs the Youngs moduli ratie In both cases a cubic lattice of spherical voids is assumed Ep/Ed and the He and Hutchinson approach (HH, dotted line)compare colation condition of pores(V=T/6=0.52). Such a choice is consistent with the process used to create porosity by addition (see Section 3) of spherical pyrolysable particles of constant diameters. 2-8 Higher volume fraction of pores can be obtained using parti- GAf=G cles of different sizes It must be pointed out that, throughout this paper, the term The function g(Eq (12)(solid line)and the toughness ra- "dense"is related to the stiffer material. The sintering leads tio(dashed line)are plotted versus the Youngs moduli ratio to a residual close porosity in the bulk(2.5% in SiC and 6% in (Fig. 7). As expected, this figure shows clearly that no de B4C) independent of that obtained by addition of pore form- flection can occur at such an interface, since the inequality ing agents. This residual porosity is formed of pores much (Eq (II)) never holds true. The He and Hutchinson(HH) smaller than those resulting from the addition of particles. It approach 0 is also plotted(dotted line)and leads to the same can be ignored in the present analysis but has to be taken into conclusion account in the measures V=V-Vo (15) 6. Deflection at the porous/dense interface where V is the actual porosity and Vo the initial(residual) one if no pore forming agent is added The so-called next layer is now a dense one. The function As a consequence of Eqs. (13)and(14), the deflection g(Eq (12))and the toughness ratio are again plotted versus criterion(Eq. (11)) can be rewritten the Youngs moduli ratio( Fig 8)at a porous/dense interface It is only assumed that the two ratios(Youngs moduli and h(V)≥H(v) with h(V)=g (16 toughness)follow the same rule(whatever this rule, i.e. Eqs (16)and(17)or any other ). The toughness ratio in Eq (I1) k(V)≥K(V) with k(V)=g (17) 0.8 The two following Figs. 5 and 6 exhibit experimental mea- res of the Youngs moduli and the toughness ratios for the Gp/Gd,/ two types of laminates and different additional pore form ing agents: corn starch, polyamide and PTFE. Data are taken from Blanks et al. Reynaud and co-workers 5-7and Tariolle 0.4 et al. 7, 8 In both cases the better matching is obtained using the surface fraction of pores(Eq(17)(dashed lines) Ep/Ed 5. Deflection at the dense/porous interface Fig 8. The function g(Eq (12))(solid line)vs the In this case, the so-called next layer is a porous one Ep ea and the He and Hutchinson approach o(HH, dotted line)compare and the toughness ratio in Eq.(9)or Eq(11)equals 1 5 to the toughness ratio Gs/G&(dashed line)at the porous/dense interface
D. Leguillon et al. / Journal of the European Ceramic Society 26 (2006) 343–349 347 the largest surface fraction of pores S in the second case (Eq. (14)): Ep = H(V)Ed; Gc p = H(V)Gc d with H(V) = 1 − 6V π (13) or Ep = K(V)Ed; Gc p = K(V)Gc d with K(V) = 1 − 4S π = 1 − 6V π 2/3 (14) where Gc p and Gc d are the toughness respectively of the porous and the dense ceramics. In both cases a cubic lattice of spherical voids is assumed and the parameters of the porous ceramic vanish at the percolation condition of pores (V = π/6 = 0.52). Such a choice is consistent with the process used to create porosity by addition of spherical pyrolysable particles of constant diameters.2–8 Higher volume fraction of pores can be obtained using particles of different sizes.9 It must be pointed out that, throughout this paper, the term “dense” is related to the stiffer material. The sintering leads to a residual close porosity in the bulk (2.5% in SiC and 6% in B4C) independent of that obtained by addition of pore forming agents. This residual porosity is formed of pores much smaller than those resulting from the addition of particles. It can be ignored in the present analysis but has to be taken into account in the measures: V = V˜ − V0 (15) where V˜ is the actual porosity and V0 the initial (residual) one if no pore forming agent is added. As a consequence of Eqs. (13) and (14), the deflection criterion (Eq. (11)) can be rewritten: h(V) ≥ H(V) with h(V) = g 1 − 6V π (16) or k(V) ≥ K(V) with k(V) = g 1 − 6V π 2/3 (17) The two following Figs. 5 and 6 exhibit experimental measures of the Young’s moduli and the toughness ratios for the two types of laminates and different additional pore forming agents: corn starch, polyamide and PTFE. Data are taken from Blanks et al.,3 Reynaud and co-workers5–7 and Tariolle et al.7,8 In both cases the better matching is obtained using the surface fraction of pores (Eq. (17)) (dashed lines). 5. Deflection at the dense/porous interface In this case, the so-called next layer is a porous one and the toughness ratio in Eq. (9) or Eq. (11) equals 115 Fig. 7. The function g (Eq. (12)) (solid line) vs. the Young’s moduli ratio Ep/Ed and the He and Hutchinson approach10 (HH, dotted line) compared to the toughness ratio Gc p/Gc d (dashed line) at the dense/porous interface. (see Section 3): Gc def = Gc pen = Gc p (18) The function g (Eq. (12)) (solid line) and the toughness ratio (dashed line) are plotted versus the Young’s moduli ratio (Fig. 7). As expected, this figure shows clearly that no de- flection can occur at such an interface, since the inequality (Eq. (11)) never holds true. The He and Hutchinson (HH) approach10 is also plotted (dotted line) and leads to the same conclusion. 6. Deflection at the porous/dense interface The so-called next layer is now a dense one. The function g (Eq. (12)) and the toughness ratio are again plotted versus the Young’s moduli ratio (Fig. 8) at a porous/dense interface. It is only assumed that the two ratios (Young’s moduli and toughness) follow the same rule (whatever this rule, i.e. Eqs. (16) and (17) or any other). The toughness ratio in Eq. (11) Fig. 8. The function g (Eq. (12)) (solid line) vs. the Young’s moduli ratio Ep/Ed and the He and Hutchinson approach10 (HH, dotted line) compared to the toughness ratio Gc p/Gc d (dashed line) at the porous/dense interface
D. Leguillon et al. /Journal of the European Ceramic Sociery 26(2006)343-349 0.3 、H 0.6 025 Fig9. The function h(Eq (16)(solid line)vs the volume fraction of pores Fig. Il. The function g(Eq (12))(solid line)vs the Young s moduli ratio and the He and Hutchinson approach(HH, dotted line)compared to the Ep/Ea for different thicknesses ratio ep/ea(0. 25, 0.5, 1, 2)compared to the toughness ratio H(V)=GS/Gs(Eq (13))(dashed line) toughness ratio Gp/Ga(dashed line)at the porous/dense interface the present analysis). Between the two values, the primary crack changes direction but kinks out immediately. In this (19) latter case, the work of fracture is not strongly increased, the The answers brought by the present analysis and the He and Hutchinson one 0 differ now significantly, the agreement 7. The influence of the porous layers thickness with the experiments will be discussed below The criterion can be also plotted versus the volume fraction Throughout this paper, the porous material is considered as of pore V. It is illustrated in the two following figures derived homogenous, this implies that the pores size is much smaller from the single Fig8. In the first one(Fig. 9), the elastic and than the layers thickness. This property must not be forgotten fracture parameters depend linearly on the volume fraction of especially in this section. The porous layers thickness can pores(Eq (16),whereas in the second (Fig. 10)they depend be diminished provided it does not interfere with the pores linearly on the surface fraction of pores(Eq. (17). Clearly diameter. Typically, at least one decade must separate these the predicted porosity that causes crack deflection(arrows in two characteristic lengths, the ratio between the pore diameter Figs. 9 and 10)is above 40% in both cases. If one does not and the layer thickness must not exceed 0.1 forget the residual porosity in the sintered parts(Eq. (15)) Numerical results show that thin porous layers tend to pro- these results are in a good agreement with the experiments mote crack deflection In Fig. 11(azoom of Fig 8 in the range of Reynaud and co-workers",and Tariolle",, while the He 0-0.3 for Ep/Ed), the function g (Eq (12)is plotted versus and Hutchinson approach 0 underestimates it. Blanks et al. the Young s moduli ratio Ep/Ed for different values of the found a wide range of values: between 34% and 44%. below thicknesses ratio pled, where ep and ed are respectively the 34%no deflection was observed. above 44% an extensive porous and dense layer thicknesses. Results are summarized deflection was obtained (this last value is in agreement with in Table 1. The model of porosity relies on function H defined in Eq (12). The porosity required to promote deflection de- creases with the porous layers thickness. Nevertheless, this effect remains small, the reduction is only about 10% when he relative porous layers thickness is divided by 8. 8. The influence of the poisson 's ratios The Poissons ratios of the components have been omitted in the above discussion. The next figure shows that they play 041 Table 1 Youngs moduli ratio Ep/Ed and pore volume fraction V promoting crack Fig 10. The function k(Eq (17)(solid line)vs the surface fraction of pores and the He and Hutchinson approach o(HH, dotted line)compared to the toughness ratio(Eq (14))K(V)= GS/GA(dashed line). 0.445 0.416
348 D. Leguillon et al. / Journal of the European Ceramic Society 26 (2006) 343–349 Fig. 9. The function h (Eq. (16)) (solid line) vs. the volume fraction of pores and the He and Hutchinson approach10 (HH, dotted line) compared to the toughness ratio H(V) = Gc p/Gc d (Eq. (13)) (dashed line). is: Gc def Gc pen = Gc p Gc d (19) The answers brought by the present analysis and the He and Hutchinson one10 differ now significantly, the agreement with the experiments will be discussed below. The criterion can be also plotted versus the volume fraction of pore V. It is illustrated in the two following figures derived from the single Fig. 8. In the first one (Fig. 9), the elastic and fracture parameters depend linearly on the volume fraction of pores (Eq. (16)), whereas in the second (Fig. 10) they depend linearly on the surface fraction of pores (Eq. (17)). Clearly the predicted porosity that causes crack deflection (arrows in Figs. 9 and 10) is above 40% in both cases. If one does not forget the residual porosity in the sintered parts (Eq. (15)), these results are in a good agreement with the experiments of Reynaud and co-workers7,8 and Tariolle8,9 while the He and Hutchinson approach10 underestimates it. Blanks et al.3 found a wide range of values: between 34% and 44%. Below 34% no deflection was observed, above 44% an extensive deflection was obtained (this last value is in agreement with Fig. 10. The function k (Eq. (17)) (solid line) vs. the surface fraction of pores and the He and Hutchinson approach10 (HH, dotted line) compared to the toughness ratio (Eq. (14)) K(V) = Gc p/Gc d (dashed line). Fig. 11. The function g (Eq. (12)) (solid line) vs. the Young’s moduli ratio Ep/Ed for different thicknesses ratio ep/ed (0.25, 0.5, 1, 2) compared to the toughness ratio Gc p/Gc d (dashed line) at the porous/dense interface. the present analysis). Between the two values, the primary crack changes direction but kinks out immediately. In this latter case, the work of fracture is not strongly increased, the design goal is not attained. 7. The influence of the porous layers thickness Throughout this paper, the porous material is considered as homogenous, this implies that the pores size is much smaller than the layers thickness. This property must not be forgotten, especially in this section. The porous layers thickness can be diminished provided it does not interfere with the pores diameter. Typically, at least one decade must separate these two characteristic lengths, the ratio between the pore diameter and the layer thickness must not exceed 0.1. Numerical results show that thin porous layers tend to promote crack deflection. In Fig. 11 (a zoom of Fig. 8 in the range 0–0.3 for Ep/Ed), the function g (Eq. (12)) is plotted versus the Young’s moduli ratio Ep/Ed for different values of the thicknesses ratio ep/ed, where ep and ed are respectively the porous and dense layer thicknesses. Results are summarized in Table 1. The model of porosity relies on function H defined in Eq. (12). The porosity required to promote deflection decreases with the porous layers thickness. Nevertheless, this effect remains small, the reduction is only about 10% when the relative porous layers thickness is divided by 8. 8. The influence of the Poisson’s ratios The Poisson’s ratios of the components have been omitted in the above discussion. The next figure shows that they play Table 1 Young’s moduli ratio Ep/Ed and pore volume fraction V promoting crack deflection for various porous/dense layers thickness ratios ep/ed ep/ed 2 1 0.5 0.25 Ep/Ed 0.120 0.150 0.175 0.205 V 0.461 0.445 0.432 0.416
D. Leguillon et al. Journal of the European Ceramic Society 26(2006)343-349 0.3 cause crack deflection,for which the assumption seems to be 0.2 Acknowledgement v=0.16 This work was supported by the French Ministry 0.1 Research in the ACI program "Surfaces and Interfa 2001-2003 02 03 References nctiong(Eq (12))(solid lines)vs the Youngs moduli ra- 1. Chartier, T, Merle, D and Besson, J. L, Laminar ceramic composites wo different Poissons ratio compared to the toughness ratio J. Eur Ceram Soc. 1995. 15. 101-10 Gp/Ga(dashed line) at the porous/dense interface 2. Clegg, w.J. Blanks, K.S., Davis, J. B. and Lanckmans, F, Porous es as crack deflecting interlayers a minor role. It is assumed that both dense and porous ceram- Eng. Mater,1997,132-136,1866-1869 ics have the same Poissons ratio. Two values are compared Blanks, K. S, Kristofferson, A, Carlstrom, E. and Clegg, W. J in Fig. 12: v=0.16(B4C)(note that v=0. 17 for SiC)and Crack deflection in ceramic laminates using porous interlayers. J. EE. Ceram.Soc.,1998,18,1945-1951 a realistic value met in many materials v=0. 3. Clearly the 4. Davis. J. B. Kristoffersson. A. Carlstrom. E. and Clegg. W. J. Fab- deviation in Youngs moduli ratio causing deflection is weak (6%)(Fig 12). Using the porosity function H(Eq (13),it leads to about l% deviation in the pore volume fraction V. It llaires monolithiques et composites en Carbure de silicin. Ph. D. thesis no. 282TD, Ecole des mines, St- is obviously a negligible effect. 6. Reynaud, properties and mechanical behaviour of Sic dense-porous laminates 9. Conclusion J.Eur. Ceram.Soe,2005,25,589597 7. Tariolle, S, Reynaud, C, Thevenot, F, Chartier, T. and Besson, J The first conclusion to draw is that deflection is very L, Preparation and mechanical properties of SiC-SiC and B:C-B laminates. J. Solid State Chem. 2004. 177. 487-492. difficult to promote by porous layers obtained by the ad 8. Tariolle, S, Carbure de bore monolithique et preux er composites dition of spherical pore forming particles. This prediction correlates well with experimental results: Reynaud and co- 328TD. Ecole des Mines. St-Etienne. France. 2004 orkers'-did not observe any extensive deflection in Sic 9. Ma, J, Wang, H, Weng, L. and Tan, G. E. B, Effec below v=42%(i. e. VN40%)and Tariolle et al. 7, found a interlayers on crack deflection in ceramic laminates. J. ignificant deflection for a rather high value V=52%(1.e Soc.,2004,24,825-831 M. Y and Va46%)in B4 C. It is in a good agreement with the present between dissimilar elastic materials. Int. J. Solids struct. 153-167 model, while the one based on the second He and Hutchin- 11 He. M. Y and Hutchinson, J. w Kinking of a crack out of an son approach, I as proposed by Clegg et al. 2-4 tends to sig- nificantly underestimate the experimental values. The other 12 terface. J. App..,1989,111,270-278 herti Tazi, O, Comportement a la rupture d'un assemblage de well-known He and Hutchinson approach neglects the lam materiaux fragiles, Ph D. thesis, University P. and M. Curie, Paris inated micro-structure of the material and gives an erroneous france. 2005 low value of the porosity causing crack deflection 13. Leguillon, D, Strength or toughness? A criterion for crack onset at notch. Eur J. Mech. A/Solids. 2002. 21.61-72 Finally, emphasis must be put on Figs. 7 and 8 that are 14. Leguillon, D, Sanchez-Palencia, E, Computation of Singular Solu- in a way intrinsic, porosity does not occur explicitly. They tions in Elliptic Problems and Elasticity. Masson, Paris, John Wiley play the role of 'master curves' and can be used whatever New York, 1987. the dependence of the elastic and fracture properties on the 15. Fujita, H, Jefferson, G. McMeeking, R.M. and Zok,F.W,Mul- porosity. The only assumption is that the Youngs modulus lite/alumina mixtures for use as porous m oxide fiber com. and the toughness of the porous material follow the same rule posites. J. Am. Ceram. Soc., 2004, 87(2), This may be wrong for small porosity values. But the present 16. Martin, E. and Leguillon, D, Energetic conditions for interfacial fa ure in the vicinity of a matrix crack in brittle matrix composites. Int. analysis deals only with large values of the porosity that can J. Solids struct.,2004,41,6937-6948
D. Leguillon et al. / Journal of the European Ceramic Society 26 (2006) 343–349 349 Fig. 12. The function g (Eq. (12)) (solid lines) vs. the Young’s moduli ratio Ep/Ed for two different Poisson’s ratio compared to the toughness ratio Gc p/Gc d (dashed line) at the porous/dense interface. a minor role. It is assumed that both dense and porous ceramics have the same Poisson’s ratio. Two values are compared in Fig. 12: ν = 0.16 (B4C) (note that ν = 0.17 for SiC) and a realistic value met in many materials ν = 0.3. Clearly the deviation in Young’s moduli ratio causing deflection is weak (<6%) (Fig. 12). Using the porosity function H (Eq. (13)), it leads to about 1% deviation in the pore volume fraction V. It is obviously a negligible effect. 9. Conclusion The first conclusion to draw is that deflection is very difficult to promote by porous layers obtained by the addition of spherical pore forming particles. This prediction correlates well with experimental results: Reynaud and coworkers5–7 did not observe any extensive deflection in SiC below V˜ = 42% (i.e. V ≈ 40%) and Tariolle et al.7,8 found a significant deflection for a rather high value V˜ = 52% (i.e. V ≈ 46%) in B4C. It is in a good agreement with the present model, while the one based on the second He and Hutchinson approach,11 as proposed by Clegg et al.,2–4 tends to significantly underestimate the experimental values. The other well-known He and Hutchinson approach10 neglects the laminated micro-structure of the material and gives an erroneous low value of the porosity causing crack deflection. Finally, emphasis must be put on Figs. 7 and 8 that are in a way intrinsic, porosity does not occur explicitly. They play the role of ‘master curves’ and can be used whatever the dependence of the elastic and fracture properties on the porosity. The only assumption is that the Young’s modulus and the toughness of the porous material follow the same rule. This may be wrong for small porosity values. But the present analysis deals only with large values of the porosity that can cause crack deflection, for which the assumption seems to be valid. Acknowledgement This work was supported by the French Ministry of Research in the ACI program “Surfaces and Interfaces 2001–2003”. References 1. Chartier, T., Merle, D. and Besson, J. L., Laminar ceramic composites. J. Eur. Ceram. Soc., 1995, 15, 101–107. 2. Clegg, W. J., Blanks, K. S., Davis, J. B. and Lanckmans, F., Porous interfaces as crack deflecting interlayers in ceramic laminates. Key Eng. Mater., 1997, 132–136, 1866–1869. 3. Blanks, K. S., Kristofferson, A., Carlstrom, E. and Clegg, W. J., ¨ Crack deflection in ceramic laminates using porous interlayers. J. Eur. Ceram. Soc., 1998, 18, 1945–1951. 4. Davis, J. B., Kristoffersson, A., Carlstrom, E. and Clegg, W. J., Fab- ¨ rication and crack deflection in ceramic laminates with porous interlayers. J. Am. Ceram. Soc., 2000, 83(10), 2369–2374. 5. Reynaud, C., C´eramiques lamellaires monolithiques et composites en Carbure de Silicium. Ph.D. thesis no. 282TD, Ecole des Mines, StEtienne, France, 2002. 6. Reynaud, C., Thevenot, F., Chartier, T. and Besson, J. L., Mechanical ´ properties and mechanical behaviour of SiC dense-porous laminates. J. Eur. Ceram. Soc., 2005, 25, 589–597. 7. Tariolle, S., Reynaud, C., Thevenot, F., Chartier, T. and Besson, J. ´ L., Preparation and mechanical properties of SiC–SiC and B4C–B4C laminates. J. Solid State Chem., 2004, 177, 487–492. 8. Tariolle, S., Carbure de Bore monolithique et poreux et composites lamellaires, ´elaboration, propri´et´es, renforcement. Ph.D. thesis no. 328TD, Ecole des Mines, St-Etienne, France, 2004. 9. Ma, J., Wang, H., Weng, L. and Tan, G. E. B., Effect of porous interlayers on crack deflection in ceramic laminates. J. Eur. Ceram. Soc., 2004, 24, 825–831. 10. He, M. Y. and Hutchinson, J. W., Crack deflection at an interface between dissimilar elastic materials. Int. J. Solids Struct., 1989, 25(9), 153–167. 11. He, M. Y. and Hutchinson, J. W., Kinking of a crack out of an interface. J. Appl. Mech., 1989, 111, 270–278. 12. Cherti Tazi, O., Comportement a la rupture d’un assemblage de ` materiaux fragiles, Ph.D. thesis, University P. and M. Curie, Paris, ´ France, 2005. 13. Leguillon, D., Strength or toughness? A criterion for crack onset at a notch. Eur. J. Mech. A/Solids, 2002, 21, 61–72. 14. Leguillon, D., Sanchez-Palencia, E., Computation of Singular Solutions in Elliptic Problems and Elasticity. Masson, Paris, John Wiley, New York, 1987. 15. Fujita, H., Jefferson, G., McMeeking, R. M. and Zok, F. W., Mullite/alumina mixtures for use as porous matrices in oxide fiber composites. J. Am. Ceram. Soc., 2004, 87(2), 261–267. 16. Martin, E. and Leguillon, D., Energetic conditions for interfacial failure in the vicinity of a matrix crack in brittle matrix composites. Int. J. Solids Struct., 2004, 41, 6937–6948