Abstract
Structure and mechanical characteristics of dense ceramic composites synthesised by reactive hot pressing of TiC–B4C powder mixtures at 1800–1950°C under 30 MPa were investigated by X-ray diffraction, scanning electron microscopy and energy dispersive X-ray spectroscopy (SEM and EDX). The results show that during hot pressing solid-phase chemical reaction 2TiC + B4C = 2TiB2 + 3C has occurred with final products like TiB2–TiC–C, TiB2–C or TiB2–B4C–C hetero-modulus composite formation with around one micrometer size carbon precipitates. The fracture toughness depends on the amount of graphite precipitation and has a distinct maximum K1C = 10 MPa m1/2 at nearly 7 vol.-% of carbon precipitate. The fracture toughness behaviour is explained by the developed model of crack propagation. Within the model, it is shown that pores (voids) and low-modulus carbon inclusions blunt the cracks and can increase ceramic toughness in some cases.
Introduction
Carbide-based ceramics are known and used for more than 100 years as they possess excellent properties such as high melting temperature, high hardness and high electrical conductivity (see for instance review works edited by Matkovich [1] and references within.). Exceptional properties make it possible to use them as structural materials under extreme conditions, at ultrahigh temperatures [2-4] and as a lightweight armour [5]. However, the main issue for all materials based on covalently bond structures is their brittleness. Classic way to enhance ceramic fracture resistance, e.g. toughness, is by building layered systems, for instance, adding metals [6] or other materials [7]. In the case of combination with metals or polymers the material heat resistance decreases significantly. Composites containing purely refractory phases do not possess such flow. In this case the toughening mechanisms can be classified under three main groups [8]. First in an intrinsic fracture toughness enhancement based on frontal process zone expansion which increases the fracture energy consumed in the damaged zone ahead of a crack tip. This includes microcracking due to localised internal stresses around the second phase inclusions [9] as well as so called transformation toughening [8]. To second group belong cases when extrinsic fracture toughness enhancement is provided by means of crack surface bridging with second phase, for instance fibers [10] or high aspect ratio grains [11]. The third group contains materials with the crack deflection mechanisms [8,11].
Lately the new way of brittle material toughening was developed for so called hetero-modulus ceramics (HMC). It is provided by combination of high Young's modulus (400–700 GPa) matrix with low Young's modulus (15–50 GPa) graphite or graphite-like boron nitride inclusions [12-14]. An additional advantage of HMC, as is mentioned in Shabalin et al. [15], is their machinability with the conventional tools. Low-modulus phases of graphite and hexagonal boron nitride are not tough materials on their own account [16]. Thus, the improvement of super hard ceramic mechanical characteristics with those phase addition should be explained with the interphase structure phenomena. In some respect a soft inclusion works in the same manner as an empty void e.g. pore.
It is well known that an increase in porosity leads to a dramatic strength decrease in most cases. For instance, Bris et al. [17] showed the effect in steels while Samborski and Sadowski [18] observed the same for alumina and magnesia. Li and Aubertin [19] summarised similar results for a wide range of non-metal materials. However, there are other works[20,21] claiming that in some cases porosity can have a negligible effect on the fracture energy and toughness or even positively influence mechanical properties [20,22]. It is worth mentioning that Rice [20] presented results for alumina fracture toughness improvement with small (up to 15 vol.-%) amount of ∼1 µm pores and rather similar data were obtained for porous MgO. Gnesin [22] has also mentioned the possibility of alumina fracture toughness increasing by the little ‘intergranular’ pores.
Theoretical analysis of crack front – void interaction was mentioned earlier by Broek [23] who considered the possibility of crack capture in a metal plate by the round hole and presented experimental investigation of such capture in aluminium-based alloys. However it was also shown that all the beneficial effects are being downturned by the weakening of the material between the hole and a crack and by the increase of the crack length as soon as the void became a part of the total defect.
One of the early works [24] has evaluated porous ceramic fracture toughness where spherical pores were shown as being able to improve matrix mechanical characteristics by blunting of the crack tip. Unfortunately, the authors did not take into account a case of the pores affecting the crack front (the case was later analysed in Martin [25]). The developed model predictions were experimentally confirmed with limited types of ceramic foams only. Yoshida et al. [26] evaluated fracture energy of porous alumina by accounting the crack branching and demonstrated by the experiments the positive effect of pores on the material fracture energy. It should be noted though that the model was only applicable to the high-porosity (>35%) materials which mechanical characteristics were not claimed being higher than that of the dense material.
In summary, although there is some theoretical and practical understanding of porosity influence upon ceramic toughness [21,23,24,26] there is still no universal explanation for all experimental results [18,20]. The existing models do not explain the behaviour of porous alumina and can hardly be used to explain the observed ceramic matrix reinforcement with low-modulus inclusions. Hence, the correct description of crack – void interaction in brittle materials as well as the link between fracture toughness and the strength of ceramics remains as an unresolved problem.
As has been shown previously [27], sintering of 2ТіС–В4С powder at 2100°C leads to TiB2 phase creation, with extra carbon segregated as graphite inclusions in diboride matrix. Such inclusions are emerging as a result of chemical reaction; they are almost spherical and thus ideal objects for an investigation.
The paper presents an investigation of the almost spherical low-modulus inclusion influence on fracture toughness of reactive hot-pressed TiC–TiB2–B4C–C HMC and demonstrates, with developed crack propagation model how ‘soft’ inclusions can toughen ceramics.
Experimental
Initial powder composition and sintering temperature.
The bulk densities of obtained materials were measured using the Archimedes method and the crystalline phases in the sintered specimens were determined by X-ray diffractometry (XRD) (DRON–4M, St. Petersburg, Russia). Microstructural observations and local compositional analysis measurements were taken with scanning electron microscopy and energy dispersive X-ray spectroscopy (SEM and EDX) (Carl Zeiss, Germany and EDAX, USA). For further investigation, some specimen surfaces were polished with diamond abrasive media. Hardness HV measurements with a Vickers indenter were performed with two different loads of 4.9 N and 147 N for 20 s holding time on the polished surfaces. Fracture toughness was estimated by measuring the crack lengths generated by Vickers indenter with a load of 147 N [28]. The heat effect and adiabatic temperature were calculated using thermochemistry data from the NIST Chemistry WebBook [29].
Results and discussion
The XRD spectra of hot-pressed ceramics are presented in Figure 1. It is evident that the sintering process causes chemical reaction between TiC and B4C components with production of TiB2 and graphite formation. This process can be described, in line with [27], by the following chemical reaction:
XRD spectra of sintered samples.

Composition* and density of the samples after hot pressing.
*Content of different phases was estimated basing on the initial powder compositions and considering reaction (1) occurring during the sintering process.
Typical structure of cleaved (not polished) stoichiometric composite (sample number 6, Tables 1 and 2) is shown on Figure 2. It should be noted that the presented images were taken by the backscatter detector in composition mode; this implies that the level of grey is roughly proportional to the average Z at the imaging point, for instance carbon inclusions will appear almost black. The first image (Figure 2) can be described as from the material constructed of large grey appearing grains (1) alternated with non-uniform mixture based on some darker appearing phase (2). The large areas of bright submicron particles (3) as well as many of darker spherical inclusions with sizes between 0.5 and 3 µm (4) should be also noted. Considering XRD results (Figure 1) the areas (1) should be titanium diboride, areas (2) – mixture of carbon and titanium diboride, (4) could be carbon inclusions and white particles – titanium carbide remnants.
SEM micrograph (fracture surface, sample 6, backscattered electrons): 1 – TiB2 grains, 2 – C-TiB2 mixture, 3 – TiC remnants, 4 – graphite clusters.
The distribution of C and Ti atoms in a polished sample 8 section was investigated with EDX (Figure 3) and XRD analysis (Table 2) proved light grey and grey areas to be TiC and TiB2 as well as black ones to be the graphite inclusions. Typical ultrathin eutectic structure of TiC–TiB2 matrix is also presented in Figure 3. TiC content in the initial powder increasing produces further composite structure changes, as is shown in Figure 4(b) where TiC-based matrix with distributed carbon inclusions is presented. As it could be expected, samples 1–5 with greater boron carbide amount have similar structure with B4C–TiB2 matrix and carbon inclusions (Figure 4(a)).
SEM backscatter image and EDX line profiles (polished section, sample 8). SEM micrographs (fracture surfaces, backscattered electron image): (a) – sample 1, (b) – sample 9.

It should be noted that most of the composites were fully densified at temperatures not more than 1850°C during 16 minutes while as is shown in Wang et al. [2] titanium diboride can be densified at similar conditions during 60 minutes and TiC–C composites were densified by Shabalin et al. [12] at temperatures more than 2200°C. Charge shrinkage in presented case can be intensified because of new phase formation during the hot pressing. Reaction (1) heat effect increases almost linearly with temperature from ΔH ≈ −170 kJ mole−1 at 1000°C to ΔH ≈ −220 kJ mole−1 at 1800°C [29]. Besides, boron carbide formation enthalpy turns positive at 1000°C thus, B4C phase becomes metastable at higher temperatures, and may begin to decompose in contact with TiC grains. Boron atoms diffuse into titanium carbide displacing carbon and forming TiB2 clusters. Carbon atoms remaining after boron carbide decomposition, as well as those of having been displaced form graphite interlayer between B4C and newly appearing TiB2 particles. Similar interlayer can appear between TiC and TiB2 as well. These soft graphite areas emerging under external pressure provoke intergranular slips and lead to the charge consolidation. It should be also noted that 16-minute densification temperature values (minimum hot-pressing temperature at which maximum possible density can be achieved during 16 minutes) are lower (1800°C) near the middle of composition interval (samples 3–6, Tables 1 and 2) where reaction mass content is higher. Latter proved that the densification process of sintered materials depends strongly on the reaction (1).
Sintered material microhardness/composition dependence (Figure 5(a)) has a minimum for stoichiometric sample 6. In general, the hardness is lower than observed in a similar ceramics created by the laser processing [30]. The same minimum can be seen in a corresponding fracture toughness curve though alongside with two distinct maximums (Figure 5(b)). It should be noted (Table 3) that the maximums are in the same graphite content. It is clear that sintered material fracture toughness depends not on the type of matrix but predominately on carbon inclusion content. So, as it was shown in [15], graphite being low-modulus phase should block cracks in the material.
Mechanical properties of sintered ceramics.
Theoretical description
To confirm the possibility of a ceramic material strengthening with the carbon clusters (inclusions) a model for porous material fracture toughness estimation has been developed. The model is based on the experimental fact of crack front bending between reinforcing particles [25] before final destruction of the material. In [25], matrix strengthening was connected with a crack front critical linear strain TC (the crack front elastic energy per length unit which should be accumulated for crack propagation) and elastic energy increasing because of front bending, but neither TC value nor reinforcing particle destruction conditions were represented.
Let us consider a loaded crack with stress intensity coefficient K1 in a brittle material. The expression for linear strain T can be obtained with integration of elastic energy density considering Irvine relations[33] for stress distribution near the crack tip:

integration,
– crack tip radius estimated as

<< C/2 expression (2) transforms into:

Then the critical linear strain TC0 can be estimated as

Here

Now we can consider uniqueness of crack propagation in a matrix containing inclusions. As it was shown in Martin [25] crack front in such structure is sagging until the particles are crashed. In this case crack propagation work (γ) includes not only energy for material cleavage (γ0) but also crack front elastic energy increasing:

Let the voids be spherical with the average radius Rp and concentration η, evenly distributed in ceramic matrix with fracture energy and Young's modulus γ0 and E0 respectively. Reduction of load-carrying intersection of the sample causes fracture energy decreasing:

Critical linear strain (5) alters in a similar way:

On the other hand, as is shown in Cherepanov [36], maximum stress σ in a crack tip with circumferential (Rp) void can be estimated as:

Thus the void leads to the significant decreasing of stress concentration, so linear strain is reduced to



Hence, if Rp > r0 a spherical pore can not only weaken the material, but also detain crack in it with effectiveness increasing with the pore radius.
Such confinement is shown in Martin [25] as leading to crack front sagging in a way presented in Figure 6 which causes matrix destruction energy increasing (7) with the elongation
The crack front is sagging between two spherical pores.

Thus matrix fracture energy is

Here R1 is the crack front curvature (see Figure 7) which can be calculated as
Crack front geometrical characteristics.




So now we have expressions to estimate effective fracture energies of ‘pore’ (13) and matrix (15) based on matrix and porosity properties as well as on crack front bending parameters (ϕ and θ). The latter will change during front bending process reaching their critical values ϕc and θc when over-pore material will at last be crashed.
To find the critical crack front curvature parameters we should consider that the crack sagging leads to stress redistribution between incurved and domed areas (see
and
, Figure 7). For the crack propagation in incurved area there should be accumulated critical linear strain TC1 (9). Linear strain unit is N, so it is natural to suppose domed/incurved area interaction to be force-like. The similar approach for dislocation movement was developed in Foreman and Makin [37]. So, taking (11) into account, we can write the expression for linear strain (T2) in
area (Figure 7):

In principle the crack could grow in a way shown in Figure 6 until T2 reaches TC1 crashing over-pore material. But such crack front bending leads to R1 decreasing (16) hence causes increasing of matrix effective fracture energy (15) and at some R1 value Figure 6 ceases representing the most energy-optimal way of crack propagation. In this case, concaved
The crack front is ‘cutting’ over-pore material.
front area begins to move as it is shown in Figure 8, when simultaneously with
lengthening
shortening occurs resulting in φ increasing. This is the way of ‘cutting’ over-pore material instead of ‘crashing’ it. Hence formula (7) for the effective crack propagation work goes to:



New (Figure 8) front elongation value can be estimated as follows:

The crack grows in a second (Figure 8) way when the front propagation energy (7) is more than that of (21) or, considering (14) and (23) the corresponding condition can be presented as

Such movement can be described as φ magnification at a constant θ and in fact is front curvature relaxation and logically leads to R1 increasing (16). Thus at some φ value the condition (24) cannot be fulfilled. Then the crack front once more can propagate in a way presented in Figure 6 which in its turn can be described as θ magnification at a constant φ. Analysis shows that during the crack growth both ways (Figs. 6 and 8) can occur swapping each other. But this process diminishes
region – the main cause of crack arresting and simultaneously increases T2 (20) and after reaching critical values of φc and θc over-pore material should be overcome in a way of crashing


To estimate these critical parameters we developed an analytical programme which begins with a straight crack front (θ = 0 and φ = 0). Then step 1 begins with θ increasing while condition (24) checking and when it is fulfilled step 2 begins with φ increasing. When at certain φ value (24) ceases being fulfilled, the front propagation way changes once more and step 1 continues. With this step by step θ and φ magnification critical front curvature φc and θc values of maximum material fracture resistance can be achieved. Using these parameters one can obtain fracture energy and fracture toughness of ceramic material with spherical pores as follows:


– Young's modulus of porous material [38].
To analyse material properties vs. void concentration dependences of the ceramic matrix materials we have used mechanical characteristics of a model material. Matrix properties were considered to be the same as those of TiC: E = 500 GPa, K1C = 5 MPa m1/2. Fracture toughness on porosity dependences for different pore radii d is presented in Figure 9(a) and clearly indicates that spherical low-modulus clusters can improve ceramic mechanical characteristics significantly. As it is shown in Figure 9(b), there is the optimal pore diameter value which is supposed to be 0.3–1 µm for the investigated matrix. Larger pores are easier to be ‘cut’ by a crack front while smaller ones are not so efficient in crack tip blunting.
Dependence of ceramic fracture toughness on voids concentration with different void diameters (a) and on void diameter at the concentration of 12% (b).
Additional experimental evidence of the conclusion is presented by Jiang et al. [14] who have showed that h-BN inclusion influence on B4C matrix and the effect is dictated by the inclusion size: 10% of submicron BN particles could strengthen the material while 2–3 µm ones just weakened it.
Further calculations were performed considering complex matrix of the investigated materials. Fracture toughness value for TiC–TiB2 matrix (6.5 MPa m1/2; samples 6–10) was used based on published data, [39] and that for B4C–TiB2 matrix (5 MPa m1/2; samples 1–6) using [40] data. Theoretical curves show fair correspondence to the experimental characteristics of sintered ceramics (Figure 10). The authors are fully aware that the presented here indentation toughness has its limitations when applied to a porous composite system but here we concentrate mainly on comparative fracture toughness of different compositions rather than on absolute fracture toughness values. On this basis the authors would like to stress the validity of the main statement: low-modulus inclusions and voids (pores) may strengthen ceramic matrix. The exact influence will depend on the ceramic properties but one most likely should have submicron pores to improve fracture toughness.
Dependence of ceramic fracture toughness on graphite cluster content: (a) – TiC–TiB2–C (samples 6–10); (b) – B4C–TiB2–C (samples 1–6).
Conclusions
A reactive hot-pressing method for TiB2–TiC–C, TiB2–C and TiB2–B4C–C HMC composite production was developed using B4C and TiС precursors.
High-speed densification at relatively low temperatures (not more than 1950°C during 16 minutes) was a result of solid-phase chemical reactions between starting components during the charge consolidation.
A model for HMC fracture toughness estimation has been developed. The model corresponded to the obtained experimental data and proved that the crack tip blunting in graphite inclusions was the basic cause of sintered composite high properties (K1C ≈ 10 MPa m1/2).
Footnotes
Acknowledgement
The authors acknowledge Prof. I. Shabalin for his helpful suggestions concerning the idea to publish this paper.
