MIATERIALS ENE S ENGINEERING ELSEVIER Materials Science and Engineering A272(1999)398-409 www.elsevier.com/locate/msea Finite element analysis of stresses associated with transformations in magnesia partially stabilized zirconia C R. Chen a, * S X. Li, Q. Zhang b State Key Laboratory for Fatigue and Fracture of Materials, Institute of Metal Research, Chinese Academy of sciences Shenyang 110015, People s Republic of b Department of Materials Science, Tshinghua University, Beijing 100084, People's Republic of China Received 10 February 1999: received in revised form 8 July 1999 Abstract Anisotropy finite element analysis was performed to study stresses associated with the tetragonal precipitates and the tetragonal monoclinic transformation in Mgo partially stabilized zirconia. Stresses were assumed to be caused by the lattice misfits between the product phase and the parent phase. In the finite element calculations, the tetragonal and monoclinic precipitates were assumed to be lenticular in shape and the anisotropic elastic constants of the cubic, tetragonal and monoclinic phases were considered. The purpose of this paper is to obtain some knowledge about how stresses respond when the microstructure is changed, and how stress fields affect microstructure development. The finite element results show that increasing a/as and creasing cy/ae can create a favorable stress field to reduce the growth rate of tetragonal precipitates during the heat treatment. Stresses associated with the single tetragonal precipitate in an infinite cubic matrix are distinctly larger than stresses associated ith the homogeneously distributed tetragonal precipitates. The twin orientations within the monoclinic particle greatly affect stresses associated with the tetragonal to monoclinic transformation. o 1999 Elsevier Science S.A. All rights reserved Keywords: Finite element analysis; Transformations; Magnesia partially stabilized zirconia 1. Introduction crack induces a transformation and how the transfor- affects crack growth, it is also necessary The knowledge about stresses associated with a solid the stresses associated with the transformation transformation is becoming important for developing Stresses and strains associated with transformations advanced ceramics and shape memory alloys. Consider- can be investigated by theoretical, experimental and able experimental evidence has shown that the tetrago- numerical methods. The Eshelby theory and the nal zirconia (t-ZrO,) transformation to the monoclinic Khachaturyan theory are two well known theoretical zirconia(m-ZrO )in the region of a propagating crack methods which are widely used for studying solid trans- can be utilized to increase the fracture toughness of formations. For example, Porter [6] applied Eshelby's partially stabilized ZrO2 and composites containing transformed-inclusion method to study transformation ZrO2 particles [1-4]. For understanding how the mi- toughening in magnesia partially stabilized zirconia crostructures develop and how to control the mi- (Mg-PSZ), and Lanteri [7] applied the theory of crostructure development, it is necessary to study the Khachaturyan to explain the habit plane and precipi- stresses that are produced as the transformation pro- tate morphology observed in Mg-PSZ. Detailed ceeds. because the microstructure that results from the about the theoretical methods can be found elsewhere. solid transformation is largely determined by the stress e.g. in Refs. [8-ll]. For the experimental study of that develops as the product phase precipitates grow strain fields associated with solid transformations, high and interact [5]. For understanding how a propagating resolution electron microscopy [12] and convergent beam electron diffraction [13] have bee However, these experimental methods can only measure chrch(@imr ac cn (C.R. Chen) strains in thin samples. due to the surface stress relax 0921-5093/99/- see front matter c 1999 Elsevier Science S.A. All rights reserved. PI:s0921-5093099)00507-9
Materials Science and Engineering A272 (1999) 398–409 Finite element analysis of stresses associated with transformations in magnesia partially stabilized zirconia C.R. Chen a,*, S.X. Li a , Q. Zhang b a State Key Laboratory for Fatigue and Fracture of Materials, Institute of Metal Research, Chinese Academy of Sciences, Shenyang 110015, People’s Republic of China b Department of Materials Science, Tshinghua Uni6ersity, Beijing 100084, People’s Republic of China Received 10 February 1999; received in revised form 8 July 1999 Abstract Anisotropy finite element analysis was performed to study stresses associated with the tetragonal precipitates and the tetragonal to monoclinic transformation in MgO partially stabilized zirconia. Stresses were assumed to be caused by the lattice misfits between the product phase and the parent phase. In the finite element calculations, the tetragonal and monoclinic precipitates were assumed to be lenticular in shape, and the anisotropic elastic constants of the cubic, tetragonal and monoclinic phases were considered. The purpose of this paper is to obtain some knowledge about how stresses respond when the microstructure is changed, and how stress fields affect microstructure development. The finite element results show that increasing at /ac and decreasing ct /ac can create a favorable stress field to reduce the growth rate of tetragonal precipitates during the heat treatment. Stresses associated with the single tetragonal precipitate in an infinite cubic matrix are distinctly larger than stresses associated with the homogeneously distributed tetragonal precipitates. The twin orientations within the monoclinic particle greatly affect stresses associated with the tetragonal to monoclinic transformation. © 1999 Elsevier Science S.A. All rights reserved. Keywords: Finite element analysis; Transformations; Magnesia partially stabilized zirconia www.elsevier.com/locate/msea 1. Introduction The knowledge about stresses associated with a solid transformation is becoming important for developing advanced ceramics and shape memory alloys. Considerable experimental evidence has shown that the tetragonal zirconia (t-ZrO2) transformation to the monoclinic zirconia (m-ZrO2) in the region of a propagating crack can be utilized to increase the fracture toughness of partially stabilized ZrO2 and composites containing ZrO2 particles [1–4]. For understanding how the microstructures develop and how to control the microstructure development, it is necessary to study the stresses that are produced as the transformation proceeds, because the microstructure that results from the solid transformation is largely determined by the stress that develops as the product phase precipitates grow and interact [5]. For understanding how a propagating crack induces a transformation and how the transformation affects crack growth, it is also necessary to study the stresses associated with the transformation. Stresses and strains associated with transformations can be investigated by theoretical, experimental and numerical methods. The Eshelby theory and the Khachaturyan theory are two well known theoretical methods which are widely used for studying solid transformations. For example, Porter [6] applied Eshelby’s transformed-inclusion method to study transformationtoughening in magnesia partially stabilized zirconia (Mg-PSZ), and Lanteri [7] applied the theory of Khachaturyan to explain the habit plane and precipitate morphology observed in Mg-PSZ. Detailed reviews about the theoretical methods can be found elsewhere, e.g. in Refs. [8–11]. For the experimental study of strain fields associated with solid transformations, high resolution electron microscopy [12] and convergent beam electron diffraction [13] have been employed. However, these experimental methods can only measure strains in thin samples. Due to the surface stress relax- * Corresponding author. E-mail address: chrch@imr.ac.cn (C.R. Chen) 0921-5093/99/$ - see front matter © 1999 Elsevier Science S.A. All rights reserved. PII: S0921-5093(99)00507-9
C.R. Chen et al. Materials Science and Engineering 4272(1999)398-409 ation, strains in thin samples are different from that in tribute perpendicularly to each bulk materials. For numerical study of stresses and plate, however, the crowded precipitates are strains associated with a transformation, the finite ele- parallel to each other, and the I precipitates are ment(FE)method is a powerful tool. The FE method perpendicular to the Mgo plates at most locations and has been applied to study mechanical problems associ- are parallel to the Mgo plates at some locations. Such ed with the transformations in nickel base supera distribution pattern is due to the thermal residual loys [14], shape memory alloys [15], Zro2 toughened stress field near the Mgo plates [30]. The t-m trans AlO3[16-18], partially stabilized zirconia [18, 19], and formation involves the largest strain energy, and it tetragonal zirconia polycrystals [20]. Although the Fe difficult for an entire t-ZrO2 precipitate to transform modeling of transformation is rather simple at present, as a single unit. Thus the twin-related thin plates are this method has great potential for improvement in formed in the m-ZrO2 particle to reduce the strain the analysis of a transformation energy resulting from the martensitic transformation In the partially stabilized zirconia(PSZ), a disper- [31]. The shear strain in each plate is large, however, on of metastable t-ZrO, in cubic zirconia (c-ZrO, long-range shear strains of the m-Zro2 particle are gives rise to a powerful toughening mechanism [21]. eliminated because of the twin- relation [1, 24, 27, 32 Because the toughening in PSz is achieved by the The m-ZrO2 particle takes the shape of the t-ZrO careful control of the microstructure and phase distri- precipitate. Cancellation of shear strains with pairs of bution [22] the knowledge of stresses associated with twin-related plates ensures that a m-ZrO2 particle fits the c-t transformation is helpful for conducting the into the space available for it in the matrix (i.e. the heat treatment of PSZ. The aim of the heat treatment space it occupied before the t-m transformation) is first, to develop the t-ZrO, precipitates to a condi- with the minimum possible strain [27]. The twin orien- ion of metastability so that they will transform to the tations in the m-ZI ffer from pre m-ZrO2 precipitates in the presence of a propagating precipitate, and are dependent on the local stress state crack, and second, to maximize the volume of criti- [33]. In some m-ZrO2 particles, the twinned plates are cally metastable t-ZrO2 precipitates [23]. Knowledge parallel to the short dimension of the precipitate, while bout stresses associated with the t-m transforma- in some other particles the twinned plates are perpen- tion is important for heat treatment and the applica dicular to the short dimension of the precipitates tion of PSZ. The t-m transformation in the process [26, 27 of heat treatment degrades the toughness of PSz and In this paper, the FE method is applied to stuc thus, should be avoided. Under the applied load, the stresses associated with a t-precipitate in an infinite strength of PSz is limited by the stress at which the c-ZrO2 matrix;(b) interactions between t-precipitates t-m transformation occurs [24], and this transforma-(c)stresses associated with the t-m transformation of tion is widely believed to be triggered by shear stresses a precipitate in an infinite c-ZrO2 matrix. There are [25]. The t-m transformation induced by a propagat- two features in the FE calculations: first, the elastic ing crack can greatly affect the growth rate and path anisotropy of c-Zro2, t-ZrO2 and m-ZrO2 are consid- of the crack ered; second, the t-ZrO2 and m-zrO2 pre In Mg-PSZ, the t-ZrO, precipitates are lenticular in cipitates are assumed to be lenticular in shape. The shape. The tetragonal c-axis is parallel to(100) of the purpose of this paper is to obtain some knowledge C-ZrO2 matrix and is the rotation axis of the particle, about how microstresses of a transformation depend thus three possible variants may occur [3, 26, 27]. The on microstructures and how microstresses affect the typical aspect ratio of particles is about 4: 1 [] or 5: 1 formation of the microstructure during the transfor [1, 7]. The equilibrium particle shape is controlled by mation the balance between the interfacial energy and strain computer simulation base the inzburg-Landau phenomenological theory was per- 2. Material data about Mg-PSZ formed to study the shape evolution of a t-ZrO, pre- cipitate and the microstructure development of t-ZrO2 2. 1. Geometry and crystallography orientations of precipitates. The simulation results show that a lense- t-ZrO2 and m-ZrO2 in Mg-PSZ like shape appears during growth of a t-ZrO2 particle, and upon further coarsening, the shape relaxes into a The t-ZrO2 precipitates in Mg-PSZ form as lenticu rhombus bounded by facets [28, 29]. During the heat lar particles, and the tetragonal c-axis is parallel to the treatment of Mg-PSZ, the arrangement of t-ZrO2 pre- short dimension of the particles. The cross-section of cipitates is influenced by the local stress distribution. the particle is shown as Fig. la. When viewed along In the complex PSZ stabilized by high content Mg the short axis, the shape of t-ZrO, is circular in the and low content Y,O, together, at the region far from ideal case, shown as Fig. lb; in real materials, due to the Mgo precipitate plates, the t-ZrO2 precipitates dis- the interaction between precipitates, when viewed
C.R. Chen et al. / Materials Science and Engineering A272 (1999) 398–409 399 ation, strains in thin samples are different from that in bulk materials. For numerical study of stresses and strains associated with a transformation, the finite element (FE) method is a powerful tool. The FE method has been applied to study mechanical problems associated with the transformations in nickel base superalloys [14], shape memory alloys [15], ZrO2 toughened Al2O3 [16–18], partially stabilized zirconia [18,19], and tetragonal zirconia polycrystals [20]. Although the FE modeling of transformation is rather simple at present, this method has great potential for improvement in the analysis of a transformation. In the partially stabilized zirconia (PSZ), a dispersion of metastable t-ZrO2 in cubic zirconia (c-ZrO2) gives rise to a powerful toughening mechanism [21]. Because the toughening in PSZ is achieved by the careful control of the microstructure and phase distribution [22], the knowledge of stresses associated with the ct transformation is helpful for conducting the heat treatment of PSZ. The aim of the heat treatment is first, to develop the t-ZrO2 precipitates to a condition of metastability so that they will transform to the m-ZrO2 precipitates in the presence of a propagating crack, and second, to maximize the volume of critically metastable t-ZrO2 precipitates [23]. Knowledge about stresses associated with the tm transformation is important for heat treatment and the application of PSZ. The tm transformation in the process of heat treatment degrades the toughness of PSZ and, thus, should be avoided. Under the applied load, the strength of PSZ is limited by the stress at which the tm transformation occurs [24], and this transformation is widely believed to be triggered by shear stresses [25]. The tm transformation induced by a propagating crack can greatly affect the growth rate and path of the crack. In Mg-PSZ, the t-ZrO2 precipitates are lenticular in shape. The tetragonal c-axis is parallel to 100 of the c-ZrO2 matrix and is the rotation axis of the particle, thus three possible variants may occur [3,26,27]. The typical aspect ratio of particles is about 4:1 [3] or 5:1 [1,7]. The equilibrium particle shape is controlled by the balance between the interfacial energy and strain energy. A computer simulation based on the Ginzburg–Landau phenomenological theory was performed to study the shape evolution of a t-ZrO2 precipitate and the microstructure development of t-ZrO2 precipitates. The simulation results show that a lenselike shape appears during growth of a t-ZrO2 particle, and upon further coarsening, the shape relaxes into a rhombus bounded by facets [28,29]. During the heat treatment of Mg-PSZ, the arrangement of t-ZrO2 precipitates is influenced by the local stress distribution. In the complex PSZ stabilized by high content MgO and low content Y2O3 together, at the region far from the MgO precipitate plates, the t-ZrO2 precipitates distribute perpendicularly to each other; near the MgO plate, however, the crowded t-ZrO2 precipitates are parallel to each other, and the parallel precipitates are perpendicular to the MgO plates at most locations and are parallel to the MgO plates at some locations. Such distribution pattern is due to the thermal residual stress field near the MgO plates [30]. The tm transformation involves the largest strain energy, and it is difficult for an entire t-ZrO2 precipitate to transform as a single unit. Thus the twin-related thin plates are formed in the m-ZrO2 particle to reduce the strain energy resulting from the martensitic transformation [31]. The shear strain in each plate is large, however, long-range shear strains of the m-ZrO2 particle are eliminated because of the twin-relation [1,24,27,32]. The m-ZrO2 particle takes the shape of the t-ZrO2 precipitate. Cancellation of shear strains with pairs of twin-related plates ensures that a m-ZrO2 particle fits into the space available for it in the matrix (i.e. the space it occupied before the tm transformation) with the minimum possible strain [27]. The twin orientations in the m-ZrO2 can differ from precipitate to precipitate, and are dependent on the local stress state [33]. In some m-ZrO2 particles, the twinned plates are parallel to the short dimension of the precipitate, while in some other particles the twinned plates are perpendicular to the short dimension of the precipitates [26,27]. In this paper, the FE method is applied to study: (a) stresses associated with a t-precipitate in an infinite c-ZrO2 matrix; (b) interactions between t-precipitates; (c) stresses associated with the tm transformation of a precipitate in an infinite c-ZrO2 matrix. There are two features in the FE calculations: first, the elastic anisotropy of c-ZrO2, t-ZrO2 and m-ZrO2 are considered; second, the t-ZrO2 and m-ZrO2 precipitates are assumed to be lenticular in shape. The purpose of this paper is to obtain some knowledge about how microstresses of a transformation depend on microstructures and how microstresses affect the formation of the microstructure during the transformation. 2. Material data about Mg-PSZ 2.1. Geometry and crystallography orientations of t-ZrO2 and m-ZrO2 in Mg-PSZ The t-ZrO2 precipitates in Mg-PSZ form as lenticular particles, and the tetragonal c-axis is parallel to the short dimension of the particles. The cross-section of the particle is shown as Fig. 1a. When viewed along the short axis, the shape of t-ZrO2 is circular in the ideal case, shown as Fig. 1b; in real materials, due to the interaction between precipitates, when viewed
C.R. Chen et al. Materials Science and Engineering 4272(1999)398-409 along the short axis, the shape of t-ZrO2 can be nearly cles the short dimension lies in the plane of the plates, rectangular. Because the tetragonal c-axis can be along and there is an angle between am and xe shown in Fig anyone of the (100> cube directions, there are three 2a; in some m-ZrO2 particles, the normal to the plates kinds of crystallography relation between t-ZrO2 pre- is along the short dimension of the particle, and there is cipitates and c-ZrO, matrix:[001]/[001]e, [001,/[100] an angle between cm and E, shown as Fig 2b. The b, and [001M/[010]. Typical aspect ratio of t-ZrO, precip- axis is normal to the plane of Fig. 2a and b itates is about 4 or 5. In the present FE calculations, the aspect ratio of the particles is assumed to be 4.5. 2 2. Lattice parameters The t-m transformation is a martensitic transfor- mation. After the t-ZrO2 particle is transformed to the For calculating stresses and strains caused by the m-zrO2 particle, the m-zrO2 particle takes the shape of misfits of lattice parameters between the product phase the t-ZrO2 particle, but the twin-related thin plates ar and the parent phase in Mg-PSZ, the lattice parameters formed in the m-ZrO2 particle. In some m-ZrO2 parti- of c-ZrO2, t-ZrO2 and m-zrO2 without constraints should be known. The lattice parameters of Mg-PSZ can be found in many publications, such as in Refs [6, 7, 26, 27, 34]. We selected the lattice parameters of Mg-PSZ from [26] as follows: c-ZrO2,a=0.50778nm t-ZrO2,a1=0.50803nm,c=0.51903nm. m-ZrO,: am=0.51172 nm, b 0.51770mm.c C[o01]t 0.53031nm,B=9891° Y【ol For understanding the stresses associated with the t-Zr0, precipitate in an(Mg, Y)-PSZ designed in Ref. 30]. the FE calculation is also performed for this material. This material is based on the Mg-PSZ with Dh=4.5 the molar ratio MgO/Zro2=14/86, and Y2O3 is added at the molar ratio Y,O3/ZrO,= 2/98, and some o- MgO-Al2O3 is also added. From Ref [30], the lattice parameters of c-ZrO2 and t-ZrO2 are as follows: c-ZrO,: a=0.5102 nm Fig. 1. Morphology of the t-ZrO, precipitate in Mg-PSZ:(a) the t-ZrO2:a1=0.5122nm,c1=0.5170nm cross-section, and(b) viewed along the short dimension 3. Elastic Under the local coordinate system xv=, shown in Figs. I and 2, the constitutive relation of each phase can be expressed as toi=[Celi where to) and fe) are stresses and strains, and the sequence of stress- es and strains into, and ai is [xr yy z p F, 2Ap xy]. Here [Cl, is the elastic constant matrix, and [C] is different for c-Zr02, t-ZrO2, and m- For c-Zro C12 C11 C1200 0 C12C12C1000 000C4400 where C=401. C=96 and C=56 GPa. The Fig. 2. The cross-section of the m-ZrO2 particle in Mg-PSZ:(a) constants are selected from Refs. [7, 35, 36] l/cu, and(b)aml/a, For t-ZrO
400 C.R. Chen et al. / Materials Science and Engineering A272 (1999) 398–409 along the short axis, the shape of t-ZrO2 can be nearly rectangular. Because the tetragonal c-axis can be along anyone of the 100 cube directions, there are three kinds of crystallography relation between t-ZrO2 precipitates and c-ZrO2 matrix: [001]t //[001]c, [001]t //[100]c and [001]t //[010]c. Typical aspect ratio of t-ZrO2 precipitates is about 4 or 5. In the present FE calculations, the aspect ratio of the particles is assumed to be 4.5. The tm transformation is a martensitic transformation. After the t-ZrO2 particle is transformed to the m-ZrO2 particle, the m-ZrO2 particle takes the shape of the t-ZrO2 particle, but the twin-related thin plates are formed in the m-ZrO2 particle. In some m-ZrO2 particles the short dimension lies in the plane of the plates, and there is an angle between am and xi , shown in Fig. 2a; in some m-ZrO2 particles, the normal to the plates is along the short dimension of the particle, and there is an angle between cm and zi , shown as Fig. 2b. The bm axis is normal to the plane of Fig. 2a and b. 2.2. Lattice parameters For calculating stresses and strains caused by the misfits of lattice parameters between the product phase and the parent phase in Mg-PSZ, the lattice parameters of c-ZrO2, t-ZrO2 and m-ZrO2 without constraints should be known. The lattice parameters of Mg-PSZ can be found in many publications, such as in Refs. [6,7,26,27,34]. We selected the lattice parameters of Mg-PSZ from [26] as follows: c-ZrO2, ac=0.50778 nm. t-ZrO2, at=0.50803 nm, ct=0.51903 nm. m-ZrO2: am=0.51172 nm, bm=0.51770 nm, cm= 0.53031 nm, b=98.91°. For understanding the stresses associated with the t-ZrO2 precipitate in an (Mg, Y)-PSZ designed in Ref. [30], the FE calculation is also performed for this material. This material is based on the Mg-PSZ with the molar ratio MgO/ZrO2=14/86, and Y2O3 is added at the molar ratio Y2O3/ZrO2=2/98, and some aMgO–Al2O3 is also added. From Ref. [30], the lattice parameters of c-ZrO2 and t-ZrO2 are as follows: c-ZrO2: ac=0.5102 nm. t-ZrO2: at=0.5122 nm, ct=0.5170 nm. 2.3. Elastic constants Under the local coordinate system xi yi zi shown in Figs. 1 and 2, the constitutive relation of each phase can be expressed as {s}i=[C]i{o}i , where {s}i and {o}i are stresses and strains, and the sequence of stresses and strains in{s}i and {o}i is [xi xi , yi yi , zi zi , yi zi , zi xi , xi yi] T. Here [C]i is the elastic constant matrix, and [C]i is different for c-ZrO2, t-ZrO2, and mZrO2. For c-ZrO2: [C]i=Ã Ã Ã Ã Ã Æ È C11 C12 C12 0 0 0 C12 C11 C12 0 0 0 C12 C12 C11 0 0 0 0 0 0 C44 0 0 0 0 0 0 C44 0 0 0 0 0 0 C44 Ã Ã Ã Ã Ã Ç É where C11=401, C12=96, and C44=56 GPa. These constants are selected from Refs. [7,35,36]. For t-ZrO2: Fig. 1. Morphology of the t-ZrO2 precipitate in Mg-PSZ: (a) the cross-section, and (b) viewed along the short dimension. Fig. 2. The cross-section of the m-ZrO2 particle in Mg-PSZ: (a) cm//ct , and (b) am//at .
C.R. Chen et al. Materials Science and Engineering 4272(1999)398-409 Cll C12 C1300 0 xi, y and =i directions. For modeling the stress field C12 C11 C1300 0 associated with a t-ZrO2 precipitate, a lenticular shape C13C13C33000 region with the elastic constants and thermal expansion 000C4400 coefficients of the t-phase is located within a region with the elastic constants and thermal expansion coeffi 00000C66 cients of the c-phase. At a given temperature To,we assume that there is no strain in both the c-phase and where C1=327,C3=264,C12=100,C13=62,C4= t-phase. When a temperature change AT= l K is uni- 59, and C66=64 GPa. These constants are selected formly applied to the whole FE model, the shape of the from Ref [35]. By comparing th c t-phase is changed due to the anisotropic thermal ex- the t-ZrO2 is softer along the cr-axis. By comparing the ansion coefficients of the t-phase defined above. Be- Cl of t-ZrO2 with [Cl; of c-ZrO2, we know the t-ZrO, cause the t-phase region and the c-phase region are is softer than the c-zro connected at the interface ange o For m-Zro t-phase region is constrained by the c-phase region, thus the interaction between the t-phase and the c- C11 C12 C130 C15 0 phase creates a stress field within and around the t-phase. Such a stress field is the same as the stress field C12C22C230C250 associated with a t-precipitate within the c-phase CI3C23C330C350 matrix 000C440C46 After a t-ZrO2 particle is transformed to a m-ZrO2 C15C25C350C550 particle, a large shear strain is created in each thin plate 000C460C66 within the m-ZrO2 particle because the lattice parame- ter B of m-ZrO2 is not 90. However, the effect of the where C1=358, C2=426, 33 =240, C12=44, due to the twin relationship between the plates, thus the C shear strain of these plates can only affect the matrix GPa. These constants are selected from Refs. [35,37, 38]. region closely near the interface. Thus we do not con- sider shear stains in the twin-related plates of the m-ZrO, particle when calculating stresses in the matrix 3. Method for realizing lattice misfits between phases As an approximate calculation, only the anisotropy of in the fe calculation volume expansion of the particle needs to consider When the microstructure of m-ZrO, is in the form as in The stresses created by a transformation are assumed Fig. 2a, i.e. cmllct, for calculating stresses created by the to be due to the lattice misfits between the produc t-m transformation, the thermal expansion coeffi phase and the parent phase. In the FE modeling, we use cients of m-ZrO2 particles are assumed to be a thermal expansion method to realize the effect of Afx=(am cos(B-900)-aDla lattice misfits. The phases are assumed to be coherently △ connected at the interfaces. The thermal expansion coefficients of c-ZrO2 matrix are assumed to be zero, fz=(cm-a)/a fz=0, where x, y and z are axes of the When the microstructure of m-ZrO, is in the form as global coordinate system xyz. The thermal expansion in Fig. 2b, i.e. am//an, the thermal expansion coefficients coefficients of t-Zro are assumed to be. of m-ZrO, particles are assumed to be △/x=( fx=/y=(a1-a-)/a △/y=(bm-a1) where xi y, and =, are axes of the local coordinate system xy of the t-ZrO2 precipitate. These thermal △fz=( C cOS(B-90°)-c1)/a expansion coefficients are not the real thermal expan In this way, the unconstrained lattice misfits between sion coefficients of the materials. The only purpose of phases are realized when a temperature change AT defining these thermal expansion coefficients in this K is applied to the whole Fe model, and the stresses paper is to realize the shape changes of the transformed created by the t-m transformation can be obtained region within the FE model. Under the free state, if the We take the axes [100], [010) and [ool] of the c-Zro temperature T is changed AT=l K, the shape of the matrix as the global coordinate system xyz. Befo C-phase will not change, while the relative dimension performing the FE calculation, the elastic constant changes of the t-phase will be x, fr, and fz, along the matrix [C], and the thermal expansion coefficients (fJ
C.R. Chen et al. / Materials Science and Engineering A272 (1999) 398–409 401 [C]i=Ã Ã Ã Ã Ã Æ È C11 C12 C13 0 0 0 C12 C11 C13 0 0 0 C13 C13 C33 0 0 0 0 0 0 C44 0 0 0 0 0 0 C44 0 0 0 0 0 0 C66 Ã Ã Ã Ã Ã Ç É where C11=327, C33=264, C12=100, C13=62, C44= 59, and C66=64 GPa. These constants are selected from Ref. [35]. By comparing C33 with C11, we know the t-ZrO2 is softer along the ct -axis. By comparing the [C]i of t-ZrO2 with [C]i of c-ZrO2, we know the t-ZrO2 is softer than the c-ZrO2. For m-ZrO2: [C]i=Ã Ã Ã Ã Ã Æ È C11 C12 C13 0 C15 0 C12 C22 C23 0 C25 0 C13 C23 C33 0 C35 0 0 0 0 C44 0 C46 C15 C25 C35 0 C55 0 0 0 0 C46 0 C66 Ã Ã Ã Ã Ã Ç É where C11=358, C22=426, C33=240, C12=144, C13=67, C23=127, C44=99, C55=79, C66=130, C15= −26, C25=38, C35= −23, and C46= −39 GPa. These constants are selected from Refs. [35,37,38]. 3. Method for realizing lattice misfits between phases in the FE calculation The stresses created by a transformation are assumed to be due to the lattice misfits between the product phase and the parent phase. In the FE modeling, we use a thermal expansion method to realize the effect of lattice misfits. The phases are assumed to be coherently connected at the interfaces. The thermal expansion coefficients of c-ZrO2 matrix are assumed to be zero, i.e. fX=fY=fZ=0, where x, y and z are axes of the global coordinate system xyz. The thermal expansion coefficients of t-ZrO2 are assumed to be: fZi =(ct−ac)/ac fXi =fYi =(at−ac)/ac (1) where xi , yi and zi are axes of the local coordinate system xi yi zi of the t-ZrO2 precipitate. These thermal expansion coefficients are not the real thermal expansion coefficients of the materials. The only purpose of defining these thermal expansion coefficients in this paper is to realize the shape changes of the transformed region within the FE model. Under the free state, if the temperature T is changed DT=1 K, the shape of the c-phase will not change, while the relative dimension changes of the t-phase will be fXi , fYi and fZi along the xi , yi and zi directions. For modeling the stress field associated with a t-ZrO2 precipitate, a lenticular shape region with the elastic constants and thermal expansion coefficients of the t-phase is located within a region with the elastic constants and thermal expansion coeffi- cients of the c-phase. At a given temperature T0, we assume that there is no strain in both the c-phase and t-phase. When a temperature change DT=1 K is uniformly applied to the whole FE model, the shape of the t-phase is changed due to the anisotropic thermal expansion coefficients of the t-phase defined above. Because the t-phase region and the c-phase region are connected at the interface, the shape change of the t-phase region is constrained by the c-phase region, thus the interaction between the t-phase and the cphase creates a stress field within and around the t-phase. Such a stress field is the same as the stress field associated with a t-precipitate within the c-phase matrix. After a t-ZrO2 particle is transformed to a m-ZrO2 particle, a large shear strain is created in each thin plate within the m-ZrO2 particle because the lattice parameter b of m-ZrO2 is not 90°. However, the effect of the shear strain of these plates on the matrix is eliminated due to the twin relationship between the plates, thus the shear strain of these plates can only affect the matrix region closely near the interface. Thus we do not consider shear stains in the twin-related plates of the m-ZrO2 particle when calculating stresses in the matrix. As an approximate calculation, only the anisotropy of volume expansion of the particle needs to consider. When the microstructure of m-ZrO2 is in the form as in Fig. 2a, i.e. cm//ct , for calculating stresses created by the tm transformation, the thermal expansion coeffi- cients of m-ZrO2 particles are assumed to be: DfXi =(am cos (b−90°)−at )/ac DfYi =(bm−at )/ac DfZi =(cm−at )/ac (2) When the microstructure of m-ZrO2 is in the form as in Fig. 2b, i.e. am//at , the thermal expansion coefficients of m-ZrO2 particles are assumed to be: DfXi =(am−at )/ac DfYi =(bm−at )/ac DfZi =(cm cos (b−90°)−ct )/ac (3) In this way, the unconstrained lattice misfits between phases are realized when a temperature change DT=1 K is applied to the whole FE model, and the stresses created by the tm transformation can be obtained. We take the axes [100], [010] and [001] of the c-ZrO2 matrix as the global coordinate system xyz. Before performing the FE calculation, the elastic constant matrix [C]i and the thermal expansion coefficients { f }i
02 C.R. Chen et al. Materials Science and Engineering 4272(1999)398-409 t-zrO2 are according to the Mg-PSZ [26 as: a 0.50778,a1=0.50803, and c=0.51903nm. Thus the lattice misfits are x=fy=(a-ae/ae=0.049%, z (c-ad/ac=2.216%. A three-dimensional model is alse established. A lenticular shape precipitate is assumed to be in an infinite large spherical c-ZrO, matrix. Due to the symmetry, only one-eighth of the precipitate needs -Zr 02 to be analyzed The model is shown as Fig. 3b. The FE results show that under plane strain assumption, al though the tensile o. of the matrix near the edge of precipitate is overestimated, characteristics of stresses of the model in Fig. 3a are quite similar to that of the model in Fig. 3b Fig. 4a shows the distribution of stress a. in the c-ZrO2 matrix near the t-zrO2 precipitate under plane for calculating the stresses associated with the single tate in an infinite c-ZrO, matrix: (a) the model und ssumption, and(b)the three-dimensional model. Ur,fr, z] of particles under the local coordinate system xv=, must be transformed to be under the global coordinate system xyz. The sequence of stresses and strains in o, and (e) is [xx, yy, zz, yz, zx, xy? However, in some FE software packages, such as AN OXx (MPa SYS, the sequence of stresses and strains is [xx, yy, zz, xy, yz, zx]. Therefore, when using these FE software packages, the sequence of elements in elastic constant matrix [C] should be adjusted. The detailed descriptions of the anisotropy FE analysis can be seen in our previous work [39, 40 4. Results and discussion 4.1. The single I-ZrO2 precipitate in an infinite c-ZrO 4.1.1. Characteristics of stresses associated with a t-ZrO2 precipitate in the c-ZrO2 matrix For investigating stresses associated with the single zro, precipitate in the infinite c-zro2 matrix, a Fe model is established as in Fig 3a. Constraints: u,=0 at G=-5 0 at 0 at y=0: u,=0 at y=8(8 is the thickness of the model). The 3-dimensional an isotropy element is used. There are only one layer of Fig. 4. The stresses in the c-ZrO, matrix near the single t-Zro elements along the y-direction, thus this is a plane precipitate. the lattice parameters are according to Mg-PSZ:(a)a2, strain model. The lattice parameters of c-ZrO2 and b)orr, and (c)
402 C.R. Chen et al. / Materials Science and Engineering A272 (1999) 398–409 Fig. 3. Models for calculating the stresses associated with the single t-ZrO2 precipitate in an infinite c-ZrO2 matrix: (a) the model under plane strain assumption, and (b) the three-dimensional model. t-ZrO2 are according to the Mg-PSZ [26] as: ac= 0.50778, at=0.50803, and ct=0.51903 nm. Thus the lattice misfits are fXi =fYi =(at−ac)/ac=0.049%, fZi = (ct−ac)/ac=2.216%. A three-dimensional model is also established. A lenticular shape precipitate is assumed to be in an infinite large spherical c-ZrO2 matrix. Due to the symmetry, only one-eighth of the precipitate needs to be analyzed. The model is shown as Fig. 3b. The FE results show that under plane strain assumption, although the tensile szz of the matrix near the edge of precipitate is overestimated, characteristics of stresses of the model in Fig. 3a are quite similar to that of the model in Fig. 3b. Fig. 4a shows the distribution of stress szz in the c-ZrO2 matrix near the t-ZrO2 precipitate under plane Fig. 4. The stresses in the c-ZrO2 matrix near the single t-ZrO2 precipitate, the lattice parameters are according to Mg-PSZ: (a) szz, (b) sxx, and (c) sxz. =[ fXi , fYi , fZi ] T of particles under the local coordinate system xi yi zi must be transformed to be under the global coordinate system xyz. The sequence of stresses and strains in {s} and {o} is [xx, yy, zz, yz, zx, xy] T. However, in some FE software packages, such as ANSYS, the sequence of stresses and strains is [xx, yy, zz, xy, yz, zx] T. Therefore, when using these FE software packages, the sequence of elements in elastic constant matrix [C] should be adjusted. The detailed descriptions of the anisotropy FE analysis can be seen in our previous work [39,40]. 4. Results and discussion 4.1. The single t-ZrO2 precipitate in an infinite c-ZrO2 matrix 4.1.1. Characteristics of stresses associated with a t-ZrO2 precipitate in the c-ZrO2 matrix For investigating stresses associated with the single t-ZrO2 precipitate in the infinite c-ZrO2 matrix, a FE model is established as in Fig. 3a. Constraints: uz=0 at z=0; ux=0 at x=0; uy=0 at y=0; uy=0 at y=d (d is the thickness of the model). The 3-dimensional anisotropy element is used. There are only one layer of elements along the y-direction, thus this is a plane strain model. The lattice parameters of c-ZrO2 and
C.R. Chen et al. Materials Science and Engineering 4272(1999)398-409 (A) (B) [oo1]c -tD/h= 00 千 TTTTn D/h D/h=3 3000 102 2500 -Oxx D/h of U/o with D/h, here Uo is the per volume elastic strain energy of precipitate when D/h=l; and(b) variation of stresses with D/h. B The dependence of the elastic strain energy on the orientation of the t-precipitate: (a) the orientation relation -precipitate and the where U is the elastic strain energy of the t-precipitate, and U is the elastic strain energy when 0=0(Uo changes with the aspect ratio v nd e C-matrix, where [01o, is assumed to coincide with [o1ols, and B is the angle between [ool, and [ool]e; and(b) the between U/Uo strain assumption. Near the side face of the precipitate, large, the corresponding shear strain ex: is large. The the o. in the matrix is compressive and the largest region around N is a stress highly concentrated region the precipitate, the o in the matrix is tensile, and the parts of the precipitate, however, it becomes large the distribution of stress o in the matrix near the ing N(the largest oz is about 500 MPa); the shear stress o is so small that it can be ignored. In the much smaller than the z-direction lattice misfit f, the C-matrix, the largest a: is above 200 MPa, and O of the matrix is distinctly smaller than o .. The a is tensile at the location N, and compressive at the MPa, and is located near Ni, but not at N. Because a location N2. Fig. 4c shows the shear stress ox of the of the t-precipitate near N, becomes a large tensile stress,and because the tensile o of the matrix at N, is matrix near the precipitate. The largest ox: is near the significantly decreased if the 3-dimensional model is location Nr. The shear stress ox: is smaller than the applied, thus we do not think that the sharp edge of the stress o. However, com the shear strain e with precipitate is a stress singular site he strain E, the larges of the matrix is several The average stresses of the t-precipitate are: o2 times the largest E. The reason is that the elastic 576 and ox=-488 MPa (Under the plane strain modulus C44 is much smaller than Cul in the matrix assumption, the compressive a is underestimated. If (C4/Cu=0.14). Thus even the shear stress ox is not the model in Fig. 3b is used, then 0=-816 MPa)
C.R. Chen et al. / Materials Science and Engineering A272 (1999) 398–409 403 Fig. 5. A. The effects of the aspect ratio D/h on the per volume elastic strain energy U and the average stresses of the t-ZrO2 precipitate: (a) varies of U/U0 with D/h, here U0 is the per volume elastic strain energy of precipitate when D/h=1; and (b) variation of the average stresses with D/h. B. The dependence of the elastic strain energy on the orientation of the t-precipitate: (a) the orientation relation between the t-precipitate and the c-matrix, where [010]t is assumed to coincide with [010]c, and u is the angle between [001]t and [001]c; and (b) the relationship between U/U0 and u, where U is the elastic strain energy of the t-precipitate, and U0 is the elastic strain energy when u=0 (U0 changes with the aspect ratio D/h). large, the corresponding shear strain oxz is large. The region around N1 is a stress highly concentrated region. Within the t-precipitate, szz is compressive in most parts of the precipitate, however, it becomes large tensile near N1, and the tensile szz increases approaching N1(the largest szz is about 500 MPa); the shear stress sxz is so small that it can be ignored. In the C-matrix, the largest szz is above 1200 MPa, and is located at N1; The largest shear stress sxz is about 390 MPa, and is located near N1, but not at N1. Because szz of the t-precipitate near N1 becomes a large tensile stress, and because the tensile szz of the matrix at N1 is significantly decreased if the 3-dimensional model is applied, thus we do not think that the sharp edge of the precipitate is a stress singular site. The average stresses of the t-precipitate are: szz= − 576 and sxx= −488 MPa. (Under the plane strain assumption, the compressive szz is underestimated. If the model in Fig. 3b is used, then szz= −816 MPa). strain assumption. Near the side face of the precipitate, the szz in the matrix is compressive, and the largest compressive szz is near the location N2; near the edge of the precipitate, the szz in the matrix is tensile, and the largest tensile szz is at the location N1. Fig. 4b shows the distribution of stress sxx in the matrix near the precipitate. Because the x-direction lattice misfit fx is much smaller than the z-direction lattice misfit fz, the sxx of the matrix is distinctly smaller than szz. The sxx is tensile at the location N1, and compressive at the location N2. Fig. 4c shows the shear stress sxz of the matrix near the precipitate. The largest sxz is near the location N1. The shear stress sxz is smaller than the stress szz. However, compare the shear strain oxz with the strain ozz, the largest oxz of the matrix is several times the largest ozz. The reason is that the elastic modulus C44 is much smaller than C11 in the matrix (C44/C11=0.14). Thus even the shear stress sxz is not
C.R. Chen et al. Materials Science and Engineering 4272(1999)398-409 4.1.2. Effect of the aspect ratio D/h of the t-zro tends to be along one of the (100> directions of the precipitate For understanding how the aspect ratio D/h of the 4.1.4. Effect of the edge thickness of the I-precipitate t-precipitate affects the stresses and elastic strain energy, on the local strain we performed FE calculations for the cases of D/h=l In our FE models, the edge thickness of the t-precip 2, 3, 4,45 and 5 by using the model in Fig. 3b. The Fe itate is assumed to be zero. From Fig. 4a we see that the results show that the larger the aspect ratio D/h is, the concentration of stress a. exists at the location n,. for smaller the stress o= and the per volume elastic strain checking if the concentration of o. at N, is due to the energy U of the precipitate will be. Fig 5Aa shows how fact that the edge of the t-precipitate is assumed too he per volume elastic strain energy U of the t-preci sharp, we blunted the sharp edge by a small circular tate varies with the aspect ratio D/. With increasing boundary. The FE results show that the tensile stress a D/, the U of the t-precipitate decreases drastically of the matrix near the blunted edge of the t-zro when D/h is small, and decreases smoothly when D/h large. The elastic strain energy provides the resistance to Precipitate is increased, thus we know that the concen tration of stress o at the location n, is not due to the transformation, thus the increasing of aspect ratio D/h sharp edge assumption of the t- precipitate in our FE is favorable to the transformation [41]. For a given models. Experimental observation [3] has shown that volume of the t-ZrO2 precipitate, with increasing aspect ratio D/h, the interfacial energy increases, while the elastic strain energy decreases. Fig 5Ab shows how the average stresses o2 and Oxr of the t-ZrO2 precipitate ary with the aspect ratio D/h. Both o, and oxxare compressive, and orr is equal to oxx When D/h is small, O2 is much larger than xr. With increasing D/h,a decreases drastically, while oxx decreases smoothly. Thus the increase of D/h can effectively reduce stress o. Why the per volume elastic strain energy decreases with the increasing aspect ratio is mainly due to the variation of stress a. and strain a 4.1.3. The dependence of elastic strain energy on orientation of the I-precipitate Experimental evidence has shown that the c-axis of the t-precipitate with a lenticular shape always tends to XX(MPa) be along one of the (100> directions of the c-matrix. For understanding the reason we assume that an infi- nitely large c-matrix contains a t-precipitate, where the [010] axis of the t-precipitate is assumed to coincide with he [010] axis of the c-matrix. As shown as in Fig. 5Ba G=180 we rotate the t-precipitate about the [010] axis in the H=300 model, and calculate how the elastic strain energy of the precipitate varies with the angle 0. The plane strain condition is assumed. Because the symmetry condition (b)0 can not be used when the [ool] axis of the t-precipitate is not along the (100> directions of the c-matrix, the whole t-precipitate must be included in the FE model Ox2(MPa) The results are shown in Fig. 5Bb. Although the effect of t-precipitate orientation on the elastic strain energy of the t-precipitate is not strong, the dependence of elastic strain energy on the orientation of the t-precipi clear. When 0=00 and 0=90, the elastic strain energy is the lowest; when 0=45, the elastic strain energy is the highest. Thus we know for a lenticular cipitate, when the c-axis of the t-precipitate is along one of thedirections of the c-matrix, the elastic strain Fig. 6. The stresses in the c-ZrO2 matrix near the single t-zrO2 energy of the t-precipitate is the smallest. This result precipitate, the lattice parameters are according to a(Mg, Y)-PSZ:(a) may partly explain why the c-axis of the t-precipitate
404 C.R. Chen et al. / Materials Science and Engineering A272 (1999) 398–409 4.1.2. Effect of the aspect ratio D/h of the t-ZrO2 precipitate For understanding how the aspect ratio D/h of the t-precipitate affects the stresses and elastic strain energy, we performed FE calculations for the cases of D/h=1, 2, 3, 4, 4.5 and 5 by using the model in Fig. 3b. The FE results show that the larger the aspect ratio D/h is, the smaller the stress szz and the per volume elastic strain energy U of the precipitate will be. Fig. 5Aa shows how the per volume elastic strain energy U of the t-precipitate varies with the aspect ratio D/h. With increasing D/h, the U of the t-precipitate decreases drastically when D/h is small, and decreases smoothly when D/h is large. The elastic strain energy provides the resistance to transformation, thus the increasing of aspect ratio D/h is favorable to the transformation [41]. For a given volume of the t-ZrO2 precipitate, with increasing aspect ratio D/h, the interfacial energy increases, while the elastic strain energy decreases. Fig. 5Ab shows how the average stresses szz and sxx of the t-ZrO2 precipitate vary with the aspect ratio D/h. Both szz and sxx are compressive, and syy is equal to sxx. When D/h is small, szz is much larger than sxx. With increasing D/h, szz decreases drastically, while sxx decreases smoothly. Thus the increase of D/h can effectively reduce stress szz. Why the per volume elastic strain energy decreases with the increasing aspect ratio is mainly due to the variation of stress szz and strain ozz. 4.1.3. The dependence of elastic strain energy on the orientation of the t-precipitate Experimental evidence has shown that the c-axis of the t-precipitate with a lenticular shape always tends to be along one of the 100 directions of the c-matrix. For understanding the reason, we assume that an infi- nitely large c-matrix contains a t-precipitate, where the [010] axis of the t-precipitate is assumed to coincide with the [010] axis of the c-matrix. As shown as in Fig. 5Ba, we rotate the t-precipitate about the [010] axis in the model, and calculate how the elastic strain energy of the t-precipitate varies with the angle u. The plane strain condition is assumed. Because the symmetry condition can not be used when the [001] axis of the t-precipitate is not along the 100 directions of the c-matrix, the whole t-precipitate must be included in the FE model. The results are shown in Fig. 5Bb. Although the effect of t-precipitate orientation on the elastic strain energy of the t-precipitate is not strong, the dependence of elastic strain energy on the orientation of the t-precipitate is clear. When u=0° and u=90°, the elastic strain energy is the lowest; when u=45°, the elastic strain energy is the highest. Thus we know for a lenticular shape t-precipitate, when the c-axis of the t-precipitate is along one of the 100 directions of the c-matrix, the elastic strain energy of the t-precipitate is the smallest. This result may partly explain why the c-axis of the t-precipitate tends to be along one of the 100 directions of the c-matrix. 4.1.4. Effect of the edge thickness of the t-precipitate on the local strain In our FE models, the edge thickness of the t-precipitate is assumed to be zero. From Fig. 4a we see that the concentration of stress szz exists at the location N1. For checking if the concentration of szz at N1 is due to the fact that the edge of the t-precipitate is assumed too sharp, we blunted the sharp edge by a small circular boundary. The FE results show that the tensile stress szz of the matrix near the blunted edge of the t-ZrO2 precipitate is increased, thus we know that the concentration of stress szz at the location N1 is not due to the sharp edge assumption of the t-precipitate in our FE models. Experimental observation [3] has shown that Fig. 6. The stresses in the c-ZrO2 matrix near the single t-ZrO2 precipitate, the lattice parameters are according to a (Mg, Y)-PSZ: (a) szz, (b) sxx, and (c) sxz
C.R. Chen et al. Materials Science and Engineering 4272(1999)398-409 4.1.5. The response of hen the lattice misfits between the t-ZrO2 prec and the c-ZrO, matrix are One disadvantage of Mg-PSZ during the heat treat- ment is that the t-ZrO2 precipitates grow too fast. Thi makes some t-ZrO2 precipitates easy to contact each other perpendicularly, and induces the t-m transfor- mation. As a result, some t-zrO2 precipitates are trans- formed to m-ZrO2 particles during the heat treatment and thus the toughness of Mg-PSZ is degraded. Exper- Fig. 7. Three t-zrO2 precipitates exist in an infinite c-ZrO2 matrix, for studying the interaction between particle I and particle 2. imental results [30] show that the growth rate of t-ZrO precipitates in(Mg, Y)-PSZ is slower than in Mg-PSZ, the edge of the precipitate is very sharp, and an ellip- and thus is easy for controlling the heat treatment. As soid can not suitably represent the sharp edge of the a result, contact between the t-zrO2 precipitates during precipitate. Thus in our following FE calculations, edge thickness of the t-ZrO, precipitate is assumed to t-ZrO2 precipitates become more stable. After long time be zero as a first approximation heat treatment, the toughness and strength of (b)0 Oxz(MP A=-300 D=-120 Fig. 8. The stresses in the c-ZrO2 matrix near the perpendicularly neighboring t-ZrO2 precipitates, where the lattice parameters are according to Mg-PSZ: (a)o.,(b)a and (c)a
C.R. Chen et al. / Materials Science and Engineering A272 (1999) 398–409 405 Fig. 7. Three t-ZrO2 precipitates exist in an infinite c-ZrO2 matrix, for studying the interaction between particle 1 and particle 2. 4.1.5. The response of stresses when the lattice misfits between the t-ZrO2 precipitate and the c-ZrO2 matrix are changed One disadvantage of Mg-PSZ during the heat treatment is that the t-ZrO2 precipitates grow too fast. This makes some t-ZrO2 precipitates easy to contact each other perpendicularly, and induces the tm transformation. As a result, some t-ZrO2 precipitates are transformed to m-ZrO2 particles during the heat treatment, and thus the toughness of Mg-PSZ is degraded. Experimental results [30] show that the growth rate of t-ZrO2 precipitates in (Mg, Y)-PSZ is slower than in Mg-PSZ, and thus is easy for controlling the heat treatment. As a result, contact between the t-ZrO2 precipitates during the heat treatment can be avoided, and the t-ZrO2 precipitates become more stable. After long time heat treatment, the toughness and strength of the edge of the precipitate is very sharp, and an ellipsoid can not suitably represent the sharp edge of the precipitate. Thus in our following FE calculations, the edge thickness of the t-ZrO2 precipitate is assumed to be zero as a first approximation. Fig. 8. The stresses in the c-ZrO2 matrix near the perpendicularly neighboring t-ZrO2 precipitates, where the lattice parameters are according to Mg-PSZ: (a) szz, (b) sxx, and (c) sxz.
C.R. Chen et al. Materials Science and Engineering 4272(1999)398-409 this (Mg, Y)-PSZ are distinctly better than Mg-PSZ. trolled by the No. 2 precipitate. Compare Fig. 8 with One important reason that the t-Zro2 precipitates in Fig ig. 4, near the No. I precipitate, o only changes a this(Mg, Y)-PSZ grow slower is due to the changes of little, Oxx near the No. I precipitate changes greatly, lattice misfits between the t-ZrO, precipitate and the and becomes a large compressive stress, thus the No. 2 c-ZrO, matrix by adding some Y,O, in Mg-PSZ precipitate resists the growth of the No. I precipitate The lattice parameters of this(Mg, Y)-PSZ are se- along the x-direction; shear stress between the No. I lected from [30] as follows: e=0.5102 nm, a,=0.5122 and No. 2 precipitates is reduced nm and C,=0.5170 nm. The lattice misfits between Stresses associated with two parallel t-ZrO2 precipi- t-ZrO2 precipitate and c-ZrO2 matrix are: fa=(a,-ad/ tates are also calculated. The c-axis stress is distinctly ac=0.392% and fc=(c-ad/ac=1.333%. By compar- enhanced, while the shear stress in the region between ing the lattice misfits of this(Mg, Y)-PSZ with that of the precipitates is reduced (Mg, Y-PSZ and Mg-PSZ; however, fa in this(Mg, 4.2.2. Stresses associated with the homogeneously dis- Y)-PSZ is distinctly larger than that in Mg-PSZ, and fe tributed I-precipitates in this(Mg, Y)-PSz is smaller than that in Mg-PSZ An idealized Mg- PSZ is shown in Fig. 9a: the distri When the lattice parameters ae, a, and c, are assumed bution of the t-zro2 precipitates in the c-ZrO2 matrix is to be of this(Mg, Y)-PSZ, stresses z, Orr and ox in assumed to be periodical, homogeneous and symmetri C-ZrO2 matrix near a single t-ZrO2 precipitate are cal, and the size of the precipitate is assumed to be shown as Fig. 6a-c. Compare Fig. 6 with Fig. 4, uniform. The volume fraction of t-ZrO2 is 28.5%. The characteristics of stresses o. and o. have no distinct FE model is established as in Fig. 9b. Constraints change, while the oxr distribution in Fig. 6b is much u,=0 at z=0; ux=0 at x=0: u,=0 at y=0; all nodes different from that in Fig. 4b. In Fig. 4b oxx is com- at z=w have a common unknown us; all nodes at pressive at N2, and tensile at Ni; in Fig. 6b oxr is tensile x= w have a common unknown ux; all nodes at y= 8 at N2, and compressive at N. Thus in this(Mg, Y)- have a common unknown u,(8 is the thickness of the PSZ, the compressive Oxx near the edge of the precipi- model). The 3-dimensional anisotropy element is used tate reduces the growth rate of the precipitate during the transformation. In this(Mg, Y-PSZ, the average stresses of t-precipitate are: 02=-437 MPa and 1272 MPa. Thus we know the compressive stress o of the precipitate of this (Mg, Y)-PSz is distinctly larger than that of Mg-PSZ (ox=-488 MPa) 4.2. Interaction between t-precipitates 4.2.1. Interaction between two perpendicularly neighbor ing t-ZrO, precipitate For studying the interaction between perpendicular neighboring t-ZrO2 particles, three t-ZrO2 precipitates are assumed to exist in an infinite c-ZrO, matrix as shown as Fig. 7. Because of symmetry, only one-fourth 2 of Fig. 7 needs to be analyzed. The constraints are the same as the single precipitate case in Fig 3a. The lattice parameters are according to Mg-PSZ The FE results show that in the No. I precipitate, the average stresses are:a. 515 MPa and oxx 1 125 MPa. Compared with the single precipitate case the compressive Oxx distinctly becomes larger, while the WD=145 compressive az is reduced a little. In the No. 2 precip. itate, the average stresses are o=-274 MPa and 668 MPa. Compared with the single precipitate case, the a-axis compressive stress(o) of the No. 2 precipitate is reduced. Fig. 8a-c show the stresses a in the c-Zro. matrix Fig. 9.(a) An idealized Mg-PSZ, the t-ZrO, precipitates with uniform precipitates. The stress o. is mainly controlled by the size are periodically, homogeneously and symmetrically distributed in the c-ZrO2 matrix.(b) The model for calculating the stresses in the No. I precipitate, while the stress xx is mainly con- idealized Mg-PSZ in(a)
406 C.R. Chen et al. / Materials Science and Engineering A272 (1999) 398–409 this (Mg, Y)-PSZ are distinctly better than Mg-PSZ. One important reason that the t-ZrO2 precipitates in this (Mg, Y)-PSZ grow slower is due to the changes of lattice misfits between the t-ZrO2 precipitate and the c-ZrO2 matrix by adding some Y2O3 in Mg-PSZ. The lattice parameters of this (Mg, Y)-PSZ are selected from [30] as follows: ac=0.5102 nm, at=0.5122 nm and ct=0.5170 nm. The lattice misfits between t-ZrO2 precipitate and c-ZrO2 matrix are: fa=(at−ac)/ ac=0.392% and fc=(ct−ac)/ac=1.333%. By comparing the lattice misfits of this (Mg, Y)-PSZ with that of Mg-PSZ, we know fc is much larger than fa in both this (Mg, Y)-PSZ and Mg-PSZ; however, fa in this (Mg, Y)-PSZ is distinctly larger than that in Mg-PSZ, and fc in this (Mg, Y)-PSZ is smaller than that in Mg-PSZ. When the lattice parameters ac, at and ct are assumed to be of this (Mg, Y)-PSZ, stresses szz, sxx and sxz in c-ZrO2 matrix near a single t-ZrO2 precipitate are shown as Fig. 6a–c. Compare Fig. 6 with Fig. 4, characteristics of stresses szz and sxz have no distinct change, while the sxx distribution in Fig. 6b is much different from that in Fig. 4b. In Fig. 4b sxx is compressive at N2, and tensile at N1; in Fig. 6b sxx is tensile at N2, and compressive at N1. Thus in this (Mg, Y)- PSZ, the compressive sxx near the edge of the precipitate reduces the growth rate of the precipitate during the transformation. In this (Mg, Y)-PSZ, the average stresses of t-precipitate are: szz= −437 MPa and sxx= −1272 MPa. Thus we know the compressive stress sxx of the precipitate of this (Mg, Y)-PSZ is distinctly larger than that of Mg-PSZ (sxx= −488 MPa). 4.2. Interaction between t-precipitates 4.2.1. Interaction between two perpendicularly neighboring t-ZrO2 precipitates For studying the interaction between perpendicularly neighboring t-ZrO2 particles, three t-ZrO2 precipitates are assumed to exist in an infinite c-ZrO2 matrix as shown as Fig. 7. Because of symmetry, only one-fourth of Fig. 7 needs to be analyzed. The constraints are the same as the single precipitate case in Fig. 3a. The lattice parameters are according to Mg-PSZ. The FE results show that in the No. 1 precipitate, the average stresses are: szz= −515 MPa and sxx= − 1125 MPa. Compared with the single precipitate case, the compressive sxx distinctly becomes larger, while the compressive szz is reduced a little. In the No. 2 precipitate, the average stresses are szz= −274 MPa and sxx= −668 MPa. Compared with the single precipitate case, the a-axis compressive stress (szz ) of the No. 2 precipitate is reduced. Fig. 8a–c show the stresses szz, sxx and sxz in the c-ZrO2 matrix near the t-ZrO2 precipitates. The stress szz is mainly controlled by the No. 1 precipitate, while the stress sxx is mainly controlled by the No. 2 precipitate. Compare Fig. 8 with Fig. 4, near the No. 1 precipitate, szz only changes a little; sxx near the No. 1 precipitate changes greatly, and becomes a large compressive stress, thus the No. 2 precipitate resists the growth of the No. 1 precipitate along the x-direction; shear stress between the No. 1 and No. 2 precipitates is reduced. Stresses associated with two parallel t-ZrO2 precipitates are also calculated. The c-axis stress is distinctly enhanced, while the shear stress in the region between the precipitates is reduced. 4.2.2. Stresses associated with the homogeneously distributed t-precipitates An idealized Mg-PSZ is shown in Fig. 9a: the distribution of the t-ZrO2 precipitates in the c-ZrO2 matrix is assumed to be periodical, homogeneous and symmetrical, and the size of the precipitate is assumed to be uniform. The volume fraction of t-ZrO2 is 28.5%. The FE model is established as in Fig. 9b. Constraints: uz=0 at z=0; ux=0 at x=0; uy=0 at y=0; all nodes at z=w have a common unknown uz ; all nodes at x=w have a common unknown ux; all nodes at y=d have a common unknown uy (d is the thickness of the model). The 3-dimensional anisotropy element is used. Fig. 9. (a) An idealized Mg-PSZ, the t-ZrO2 precipitates with uniform size are periodically, homogeneously and symmetrically distributed in the c-ZrO2 matrix. (b) The model for calculating the stresses in the idealized Mg-PSZ in (a)
C.R. Chen et al. Materials Science and Engineering 4272(1999)398-409 G=700 2 △ XX(MPa) 3Oxz(MPa) A=-280 E=70 Fig. 10. The stresses in the c-Zro, matrix of the idealized Mg-PSZ (a)a, and(b) There are only one layer of elements along the y-direc tion. The lattice misfits between the t-ZrO2 precipitate and the c-ZrO2 matrix are according to Mg-PSZ The FE results show that the average stresses in the No I precipitate are: 02=-62 MPa and Oxr=-377 Fig. Il. The stresses created by m transformation when the twin structure of mZrO, is in the c. The stresses are shown MPa. The oz is much smaller than in the single precip- at the z-x plane: (a)Ao(b)Ac (c) △a itate case. Fig. 10a and b show the distributions of o. and o in the matrix. The a distribution near the no I precipitate is the same as the a distribution near the misfits between the precipitate and the matrix. The No. 2 precipitate, thus the ox distribution does not lattice parameters are according to Mg-PSZ. If the twin need to be shown. Compare Fig. 10 with Fig. 4, at orientation of the m-ZrO2 particle is in the form of Fig location NI, O is still a high tensile stress, while the 2a, then the changes of lattice misfit between precipitate compressive o2 near the side face of the No. I precipi and matrix during the t-m transformation are: Afz tate is greatly reduced. A large shear stress x exists (cm -CD/ac= 2.221%, Afy=(bm -as/ac=1.904% and near the edge of the precipitate. From this calculation Afx=(am cos(B-90)-adlac=-0489%. Because we know that when t-ZrO, precipitates distribute ho- Afr is large, we use the three-dimensional model in Fig nogeneously, stresses will be significantly smaller than 3b. The changes of the average stresses of precipitate the single precipitate case during the t→ m transformation are:△a=-1202 MPa and Aorr=-612 MPa. Fig. lla-c show the 4.3. Stresses created by the t-m transformation stresses△,△σxand△ a in the c-zro2 matrix at the z-x plane created by the t-m transformation. Com- For understanding the characteristics about the paring Fig. ll with Fig. 4, the characteristics of the stresses created by the t-m transformation, an m- stress distributions are similar ZrO, particle is assumed to be in an infinite c-ZrO If the twin orientation of the m-ZrO, particle is in the matrix. The stresses created by the t-m transforma- form of Fig. 2b, then the changes of the lattice misfit tion are assumed to be due to the change of lattice between the precipitate and the matrix during the t-m
C.R. Chen et al. / Materials Science and Engineering A272 (1999) 398–409 407 Fig. 10. The stresses in the c-ZrO2 matrix of the idealized Mg-PSZ: (a) szz, and (b) sxz. Fig. 11. The stresses created by the tm transformation when the twin structure of m-ZrO2 is in the form cm//ct . The stresses are shown at the z–x plane: (a) Dszz, (b) Dsxx, and (c) Dsxz. There are only one layer of elements along the y-direction. The lattice misfits between the t-ZrO2 precipitate and the c-ZrO2 matrix are according to Mg-PSZ. The FE results show that the average stresses in the No. 1 precipitate are: szz= −62 MPa and sxx= −377 MPa. The szz is much smaller than in the single precipitate case. Fig. 10a and b show the distributions of szz and sxx in the matrix. The sxx distribution near the No. 1 precipitate is the same as the szz distribution near the No. 2 precipitate, thus the sxx distribution does not need to be shown. Compare Fig. 10 with Fig. 4, at location N1, szz is still a high tensile stress, while the compressive szz near the side face of the No. 1 precipitate is greatly reduced. A large shear stress sxz exists near the edge of the precipitate. From this calculation we know that when t-ZrO2 precipitates distribute homogeneously, stresses will be significantly smaller than the single precipitate case. 4.3. Stresses created by the tm transformation For understanding the characteristics about the stresses created by the tm transformation, an mZrO2 particle is assumed to be in an infinite c-ZrO2 matrix. The stresses created by the tm transformation are assumed to be due to the change of lattice misfits between the precipitate and the matrix. The lattice parameters are according to Mg-PSZ. If the twin orientation of the m-ZrO2 particle is in the form of Fig. 2a, then the changes of lattice misfit between precipitate and matrix during the tm transformation are: DfZ= (cm−ct )/ac=2.221%, DfY=(bm−at )/ac=1.904% and DfX=(am cos (b−90°)−at )/ac= −0.489%. Because DfY is large, we use the three-dimensional model in Fig. 3b. The changes of the average stresses of precipitate during the tm transformation are: Dszz= −1202 MPa and Dsxx= −612 MPa. Fig. 11a–c show the stresses Dszz, Dsxx and Dsxz in the c-ZrO2 matrix at the z–x plane created by the tm transformation. Comparing Fig. 11 with Fig. 4, the characteristics of the stress distributions are similar. If the twin orientation of the m-ZrO2 particle is in the form of Fig. 2b, then the changes of the lattice misfit between the precipitate and the matrix during the tm