Constitutive modeling of alumina sintering: grain-size effect on dominant densification mechanism

Constitutive modeling of alumina sintering: grain-size effect on dominant densification mechanism

Computational Materials Science 32 (2005) 196–202 www.elsevier.com/locate/commatsci Constitutive modeling of alumina sintering: grain-size effect on d...

317KB Sizes 0 Downloads 53 Views

Computational Materials Science 32 (2005) 196–202 www.elsevier.com/locate/commatsci

Constitutive modeling of alumina sintering: grain-size effect on dominant densification mechanism Zeming He *, J. Ma School of Materials Engineering, Nanyang Technological University, Nanyang Avenue, Singapore 639798, Singapore Received 25 March 2004; received in revised form 3 August 2004; accepted 25 August 2004

Abstract In the present work, alumina powders with the initial grain sizes of 0.9 and 7.0 lm, respectively, were sintered at different temperatures. Constitutive laws for densification were employed to model the sintering process of alumina ceramics. Based on the constitutive laws employed and the experimental results obtained, the dominant densification mechanism was identified and the effect of grain size on dominant densification mechanism was discussed. The activation energy for densification was also evaluated. In the investigated sintering temperature range, interface reaction was identified as the controlling process in sintering of alumina powders with the initial grain size of 0.9 lm, while grainboundary diffusion was identified as the dominant process in sintering of alumina powders with the initial grain size of 7.0 lm. The activation energies for densification of the finer and coarser grain size alumina ceramics were determined as 342 and 384 kJ mol1, respectively, which provided a strong support on the densification mechanism investigation. Ó 2004 Elsevier B.V. All rights reserved. PACS: 65.70.+y; 66.10.Cb Keywords: Constitutive modeling; Densification; Grain size; Grain-boundary diffusion; Interface reaction; Activation energy

1. Introduction Ceramic products are generally fabricated from the powder metallurgy route. Sintering is one of

*

Corresponding author. Tel.: +65 6790 4590; fax: +65 6790 9081. E-mail address: [email protected] (Z. He).

the key procedures in the process for achieving materials with good performance. If residual stress or density inhomogeneity is present during densification, the formed defects could result in the property deterioration of the prepared components. Therefore, good sintering process control is necessary. It requires a full understanding of the constitutive laws which govern the densification behavior.

0927-0256/$ - see front matter Ó 2004 Elsevier B.V. All rights reserved. doi:10.1016/j.commatsci.2004.08.006

Z. He, J. Ma / Computational Materials Science 32 (2005) 196–202

A sintering constitutive law is a mathematic model used to describe materials densification under a specific mechanism. Several mechanisms are involved in the densification process, and the exact mechanism in operation, however, depends on the sintering conditions such as temperature, pressure, and material parameters such as grain size and porosity. It is well accepted that diffusional mechanisms dominate the sintering process of ceramic component at high temperature and low pressure. It is also noted that the dominant path for diffusion is usually along the grain boundary [1]. The classical grain-boundary diffusion mechanism assumes that grain boundaries act as perfect sources and sinks for the diffusing atoms during the diffusion process and the energy provided is all available to drive the diffusional flux along the grain boundaries. In practice, however, some energy is expended for materials to be added to, or removed from, a grain boundary, that is, an interface reaction occurs [2,3]. If the process of interface reaction is taken into account, the rate of densification can then be controlled either by the rate of diffusional transport between the sources and sinks, or the rate at which the sources and sinks can provide, or accept, material for the diffusional process. Furthermore, it is noted that the grain size has the important correlation with the dominance of the two processes [4]. To investigate the effect of grain size on the dominant densification mechanism, alumina powders with the initial grain size of 0.9 and 7.0 lm were pressurelessly sintered at different temperatures in the present work. Constitutive laws for densification [5–7] were employed to model the sintering process of alumina ceramics. Based on the constitutive laws employed and the experimental results obtained, materials dominant densification mechanism, grain-boundary diffusion or interface reaction, is identified and the effect of grain size on dominant densification mechanism is discussed. The activation energies for densification of the prepared alumina ceramics were also evaluated. The sintering stage investigated in the present work is limited to Stage 1, where the pores form interconnected channels when the relative density of the component is less than 0.9 [8].

197

2. Constitutive laws When grain-boundary diffusion is the controlling process, the volumetric strain rate, e_ v , is expressed by the model [5] as  3 e_ 0b L0 e_ v ¼ 9 fb ðqÞðrm  rs Þ ð1Þ r0 L for Stage 1 sintering, where the subscript ÔbÕ indicates grain-boundary diffusion controlled sintering. e_ 0b is the uniaxial strain rate under an applied stress r0 for a material of the initial grain size L0. L is grain size of the material, fb(q) is the dimensionless function of the relative density q, rm is the mean stress and rs is the sintering stress. Following the model [5], the volumetric strain rate in Stage 1 sintering for interface reaction controlled sintering is expressed as "  32  3 #   L0 re j3ðrm  rs Þj 2 þ fr ðqÞ fr ðqÞ cr ðqÞ e_ v ¼ 3_e0r r0 L r0  rffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi r0 3ðrm  rs Þ ð2Þ  j3ðrm  rs Þj r0

where the subscript ÔrÕ indicates interface reaction controlled sintering and re is effective stress. fr(q) and cr(q) are dimensionless functions of the relative density. Eqs. (1) and (2) are general expressions to describe material densification behavior when grainboundary diffusion and interface reaction are the dominant mechanisms respectively. However, specific forms of the variables shown in these equations are required for modeling the sintering process. The dimensionless density functions are mechanism dependent. For the grain-boundary diffusion controlled sintering, Du and Cocks [6] proposed the expression for fb(q) during Stage 1 sintering fb ðqÞ ¼

0:54ð1  q0 Þ qðq  q0 Þ

2

2

ð3Þ

where q0 is the green density of the compact. For the interface reaction controlled sintering, the dimensionless density functions fr(q) and cr(q) for Stage 1 sintering proposed by Du and Cocks [7] are fr ðqÞ ¼

1 ð1  q0 Þ 3 q32 q12 ðq  q Þ 0 0

ð4Þ

198

Z. He, J. Ma / Computational Materials Science 32 (2005) 196–202

and cr ðqÞ ¼

2 ð1  q0 Þ 3 q32 q12 ðq  q Þ 0 0

ð5Þ

It is considered that the sintering stress is the function of relative density and grain size, but independent of densification mechanism. In Stage 1, the expression of the sintering stress is to be [6]   3ð2cs  cb Þ 2 2q  q0 q rs ¼ ð6Þ L 1  q0 where cs and cb are surface energy per unit area and grain-boundary energy per unit area respectively. The relationship between the rate constant for densification, e_ 0 (_e0b , or e_ 0r ), and the absolute sintering temperature can be expressed as   Qd e_ 00 _e0 ¼ exp ð7Þ RT RT where e_ 00 is a material constant independent of temperature, and Qd is the activation energy for densification. T and R are absolute sintering temperature and the gas constant, respectively.

3. Experimental procedure The raw materials used in the experiment were 99.9% a-alumina powders (Baikowski, USA) with initial grain sizes of 0.9 and 7.0 lm, respectively. An Enerpac hydraulic press was used to form the powders into pellets under a uniaxial pressure of 200 MPa. After dry pressing, the green compact was subjected to sintering in a high-temperature furnace (Carbolite, UK). The sintering profile started with a heating rate of 2 K min1 up to 600 °C, and then progressed at 15 K min1 until the holding sintering temperature. The sintering temperatures chosen were 1400, 1450, and 1500 °C for 0.9 lm powders, and 1600, 1650, and 1700 °C for 7.0 lm powders. At a given sintering temperature, different holding time was set for achieving the samples with a certain density range. Finally, the compacts were cooled at a rate of 20 K min1 to room temperature.

The as-sintered sample was weighed using a Scientech SA 510 analytical balance and its diameter and the thickness were measured by a vernier calliper and a digital micrometer. The relative density was calculated using the measured density divided by the theoretical density of alumina (3.98 g cm3). The thermally-etched sample was gold coated and imaged using a JSM-5310 scanning electron microscope (SEM) (JEOL, Japan). The grain size was determined by counting the number of grains traversed by different straight lines of known length. The final mean grain size was taken as 1.56 times the average intercept length [9]. The densification rate of the component at a given sintering temperature was obtained from the gradient of the curve relative density against holding time. Software Mathcad (Mathsoft, USA) was used to do the modeling.

4. Results and discussion As an illustration, the SEM image of the sample sintered at 1450 °C for 300 min is shown in Fig. 1. As stated, the grain sizes of the samples were measured from the SEM images. Fig. 2 shows the grain size-relative density trajectories of the samples sintered at different temperatures. For both the finer and the coarser ceramics, it is noted that the grain sizes obtained at different sintering temperatures but at a fixed relative density are approximately the same. This phenomenon has been previously

Fig. 1. SEM image of sample with the initial grain size of 0.9 lm sintered at 1450 °C for 300 min.

Z. He, J. Ma / Computational Materials Science 32 (2005) 196–202 16 F1400 F1450 F1500 C1600 C1650 C1700

14 12 10

L (µm) 8 6 4 2 0 0.5

0.6

0.7

ρ

0.8

0.9

1

Fig. 2. Grain size-relative density trajectories of samples sintering at different temperatures (ÔFÕ refers to the sample with the initial grain size of 0.9 lm and ÔCÕ refer to the sample with the initial grain size of 7.0 lm).

described by Gupta [10], and it also occurred in our earlier work of sintering of PZT ceramics [11]. In Eqs. (1) and (2), re and rm are zero for pres_ when sureless sintering. The densification rate, q, grain-boundary diffusion and interface reaction dominant are derived as  3 L0 q_ ¼ Ab q fb ðqÞrs ð8Þ L and

  L0 2 q_ ¼ Ar q f ðqÞr2s L r

ð9Þ

respectively, where Ab and Ar are constants. Software Mathcad was used to perform the modeling. A simple Euler scheme was adopted, where the change of the relative density after a time step, Dt, can be expressed as qtþDt ¼ qt þ q_ t Dt

ð10Þ

In Eq. (10), the specific forms of densification rate are expressed by Eqs. (8) and (9), because both the grain-boundary diffusion and the interface reaction mechanisms are assumed and later one of them will be identified. During modeling, the green densities were taken as 0.5083 for finerand 0.5500 for coarser-grained ceramics. The values of cs and cb were taken as 1 and 0.5 J/m2, respectively [12]. To avoid the added complication of the selection of the grain-growth law, the grain

199

size was converted to the function of the relative density, with respect to Fig. 2. The constant Ab and Ar, however, cannot be chosen directly. They can be obtained from the modeling process, i.e. the appropriate rate constant can be determined by least square fitting of the experimental data to the constitutive model [11]. Figs. 3 and 4, respectively, show the modeling curves for the densification behavior (relative density versus holding time) of the finer and the coarser powders sintered at different temperatures. For comparison, the experimental points are also shown in the corresponding figures. For the densification behavior of 0.9 lm alumina powder in the investigated sintering temperature range, it is found that the grain-boundary diffusion mechanism always overestimates the densification rate as holding time increases, while the interface reaction mechanism provides a good agreement between modeling and experiment (Fig. 3). On the contrary, for the alumina samples with the initial grain size of 7.0 lm, a reverse trend is shown. That is, as holding time increases, interface reaction mechanism always overestimates the densification rate, while grain-boundary diffusion gives good prediction to the experimental results (Fig. 4). From the modeling results, it is considered that interface reaction is a more appropriately dominant densification mechanism for the sintering of fine-grained (with the initial grain size of 0.9 lm) alumina in the sintering temperature range of 1400–1500 °C. On the other hand, grain-boundary diffusion is a more appropriately dominant densification mechanism for the sintering of coarsegrained (with the initial grain size of 7.0 lm) alumina in the sintering temperature range of 1600–1700 °C. Grain-boundary diffusion and interface reaction are two sequential processes in the densification of materials. Therefore, the slower mechanism would dominate the sintering process. The grain-size effect on the dominant densification mechanism could be explained by considering the distance of the diffusion path from matter source to sink. The smaller the grain size is, the shorter the transport path and hence the faster the grainboundary diffusion. Therefore, the interface reaction becomes the controlling process for sintering

200

Z. He, J. Ma / Computational Materials Science 32 (2005) 196–202

1

1

0.9

0.9 0.8

0.8

ρ

ρ

0.7

0.7 Experiment

0.6

Experiment

0.6

Model GB

Model GB

Model IR

Model IR

0.5 0

500

1000

1500

t (min)

(a)

0.5

0

500

(b)

1000

1500

t (min)

1 0.9 0.8

ρ 0.7 Experiment

0.6

Model GB Model IR

0.5

0

200

400

600

t (min)

(c)

Fig. 3. Modeling results of densification process of 0.9 lm initial grain size samples sintered at different temperatures: (a) 1400 °C (b) 1450 °C and (c) 1500 °C (ÔGDÕ and ÔIRÕ refer to grain-boundary diffusion and interface reaction respectively).

of finer grain size. As grain size increases, that is as the distance that the materials have to travel along grain boundary increases, the grain-boundary diffusion becomes more important. With respect to Fig. 2, it is noted that the difference in grain sizes at different sintering temperatures but at a given relative density is minimal. If the grain-size difference is ignored, the relationship between the densification rate at a fixed density and the activation energy for densification is expressed from Eq. (7) as _ Þ¼ lnðqT

Qd þC RT

where C is a constant.

ð11Þ

Fig. 5 shows the variations of the logarithm of _ with the reciprocal of the absolute sintering qT temperature of the two investigated alumina ceramic series. The activation energies for densification were evaluated from the slopes in the figure. The values determined are 342 kJ mol1 for alumina with 0.9 lm initial grain size and 384 kJ mol1 for alumina with 7.0 lm initial grain size. The activation energy for densification of the fine-grained ceramics is close to that measured by Ma [13] with 325 kJ mol1, when interface reaction is the dominant mechanism. The energy value of the coarse-grained ceramics is close to that reported by Frost and Ashby [14] with 419 kJ mol1, when grain-boundary diffusion is the dominant mechanism.

Z. He, J. Ma / Computational Materials Science 32 (2005) 196–202

0.8

0.8

0.7

0.7

ρ

201

ρ 0.6

0.6 Experiment

Experiment

Model GB

Model GB

Model IR

0.5

0

150

Model IR

300

450

t (min)

(a)

0.5

0

150

300

450

t (min)

(b)

0.8

0.7

ρ 0.6 Experiment Model GB Model IR

0.5

0

100

200

300

t (min)

(c)

Fig. 4. Modeling results of densification process of 7.0 lm initial grain size samples sintered at different temperatures: (a) 1600 °C (b) 1650 °C and (c) 1700 °C (ÔGDÕ and ÔIRÕ refer to grain-boundary diffusion and interface reaction respectively).

2

1

In (ρT) 0 -1 -1 (min K ) •

Finer grain size Coarser grain size Linear fitting (Coarser grain size) Linear fitting (Finer grain size)

-1

-2 4.8

5.2

5.6

6 -1

From the values of the activation energies determined in the present work, it is noted that the results provide a strong support on the dominant densification mechanism investigation. Interface reaction is the appropriate mechanism for describing the densification process of fine-grained alumina, and grain-boundary diffusion is the appropriate mechanism for describing the densification process of coarse-grained alumina.

6.4

-1

10000×T (K )

_ with reciprocal of Fig. 5. Variations of logarithm of qT absolute sintering temperature at relative density of 0.8 in finer-grained samples and of 0.7 in coarser-grain samples.

5. Conclusions Interface reaction was identified as the controlling process in sintering of alumina powders with

202

Z. He, J. Ma / Computational Materials Science 32 (2005) 196–202

the initial grain size of 0.9 lm in the temperature range of 1400–1500 °C. Grain-boundary diffusion was identified as the dominant mechanism in sintering of alumina powders with the initial grain size of 7.0 lm in the temperature range of 1600– 1700 °C. It was noted that grain size was an indication of the distance of the diffusion path and it had important effect on dominant densification mechanism. The mechanism could transit from interface reaction to grain-boundary diffusion as grain size increased. The activation energies for densification of the investigated fine- and coarse-grained alumina ceramics were determined as 342 and 384 kJ mol1, respectively, which provided the strong supporting data for dominant densification mechanism investigation.

References [1] W.D. Kingery, H.K. Bowen, D.R. Uhlmann, Introduction to Ceramics, John Wiley and Sons, New York, 1975. [2] M.F. Ashby, On interface-reaction control of Narbarro– Herring creep and sintering, Scripta Metall. 3 (1969) 837– 842. [3] B. Burton, Interface reaction controlled diffusional creep: a consideration of grain boundary dislocation climb sources, Mater. Sci. Eng. 10 (1972) 9–13.

[4] J. Ma, A.C.F. Cocks, A constitutive model for the sintering of fine grain alumina, in: N.A. Fleck, A.C.F. Cocks (Eds.), IUTAM Symposium on Mechanics of Granular and Porous Materials, Cambridge, UK, 1996, pp. 129–138. [5] A.C.F. Cocks, The structure of constitutive laws for the sintering of fine grained materials, Acta Metall. Mater. 42 (1994) 2191–2210. [6] Z.Z. Du, A.C.F. Cocks, Constitutive models for the sintering of ceramic components: I. Material models, Acta Metall. Mater. 40 (1992) 1969–1979. [7] Z.Z. Du, A.C.F. Cocks, Sintering of fine-grained materials by interface reaction controlled grain boundary diffusion, Int. J. Solids Struct. 31 (1994) 1429–1445. [8] M.F. Ashby, Sintering and isostatic pressing diagrams, Technical report, University of Cambridge, Cambridge, 1990. [9] M.N. Rahaman, L.C. De Jonghe, R.J. Brook, Effect of shear stress on sintering, J. Am. Ceram. Soc. 69 (1986) 53– 58. [10] T.K. Gupta, Possible correlation between density and grain size during sintering, J. Am. Ceram. Soc. 55 (1972) 276–277. [11] Z. He, J. Ma, Constitutive modeling of the densification of PZT ceramics, J. Phys. Chem. Solids 64 (2003) 177–183. [12] J. Zhao, M.P. Harmer, Effect of pore distribution on microstructure development: I. Matrix pores, J. Am. Ceram. Soc. 71 (1988) 113–120. [13] J. Ma, Constitutive modeling of the densification of porous ceramic components, Ph.D. thesis, University of Cambridge, Cambridge, 1997. [14] H.J. Forst, M.F. Ashby, Deformation Mechanism Maps, Pergamon, Oxford, 1982.