Fracture Toughness of 3Y-TZP Dental Ceramics by Using Vickers Indentation Fracture and SELNB Methods

Article information

J. Korean Ceram. Soc.. 2019;56(1):37-48
Publication date (electronic) : 2018 December 19
doi : https://doi.org/10.4191/kcers.2019.56.1.01
*Art & Architecture Faculty, Yadegar-e Imam Khomeini (RAH) Shahre-Rey Branch, Islamic Azad University, Tehran, Iran
**Faculty of Materials & Manufacturing Processes, Malek-Ashtar University of Technology, Tehran, Iran
Corresponding author: Alireza Moradkhani, E-mail: moradkhani.a@srbiau.ac.ir, Tel: +98-9125907184, Fax: +98-2122935341
Received 2018 May 17; Revised 2018 August 02; Revised 2018 August 14; Accepted 2018 August 17.

Abstract

The objective of this research is to analyze the fracture toughness of pure and silica co-doped yttria-stabilized tetragonal zirconia polycrystal (3Y-TZP) bioceramics containing 0.1 and 0.2 wt.% of alumina, and sintered at a temperature of 1500 °C. Because of the relatively easy preparation of the test specimens and the high speed of testing, the Vickers indentation fracture (VIF) technique is more frequently used to evaluate the fracture toughness of biomaterials and hard biological tissues. The Young’s modulus and hardness values were obtained by means of nanoindentation and indentation methods. The fracture toughness values of 3Y-TZP bioceramics were calculated and analyzed using 15 equations related to the VIF technique, and loadings of 49.03 and 196.13 N with a Vickers diamond. For validation, the results were compared with fracture toughness values obtained by the single-edge laser-notch beam (SELNB) method with an almost atomically sharp laser-machined initial notch.

1. Introduction

Stabilized zirconia, a polycrystalline ceramic material, is an important structural ceramic that is used as a biomaterial for dental and bone implants.1,2) Yttria-stabilized tetragonal zirconia polycrystal (3Y-TZP) is known to be an important structural ceramic with excellent mechanical properties, such as high fracture toughness, strength, and hardness.35) In clinical applications, the performance of dental ceramics is influenced by several parameters, including fracture toughness.6) 3Y-TZP with 3 mol% Y2O3.711) as a phase stabilizer is well-known for dental application by machining pre-sintered blocks.7)

The ability to preserve the tetragonal phase in PSZ at room temperature leads to favorable mechanical properties, especially a high KIC value.11) When stress in Y-TZP ceramic causes a crack, (t) → (m) phase transformation occurs at the crack tip; as a result, an expansion in the ceramic causes compressive stresses around the crack and prevents its propagation, and thus, KIC increases. This reinforcing mechanism is known as the transformation toughening mechanism, and it makes zirconia much tougher than other ceramic materials.1214) For example, the fracture toughness of Y-TZP, which is used as a core material in three- or four-unit fixed dental prostheses (FDPs), has been reported up to 10 MPa.m1/2.7,15) Another important factor that should be noted regarding Y-TZP ceramics is their susceptibility to aging or low-temperature degradation (LTD). The LTD process occurs by itself and is time-dependent. It also affects the (t) → (m) phase transformation mechanism16) and may also reduce the physical and mechanical properties, such as density and fracture toughness.17,18)

The fracture resistance of brittle materials such as bioceramics is evaluated by their KIC value. In fact, KIC is the critical value of stress intensity factor in the opening mode where there is the occurrence of fracture initiation and unstable crack propagation.14,19,20)

There are different techniques for measuring the KIC values of materials, which can be divided into two groups. One group includes conventional fracture mechanics and allows the determination of KIC values by using notches and induced precracks. This group of techniques, which has been employed for many years to determine the KIC values of different parameters, includes methods such as Single Edge Notch Beam (SENB), Single Edge Precracked Beam (SEPB), and Chevron Notched Beam (CNB). Some of these methods have also been used for obtaining the KIC values of biomaterials.21,22) The second group, which is known as indentation fracture (IF) methods, is based on the use of a sharp indenter such as a Vickers/Knoop indenter. In general, depending on the type of material and the required degree of accuracy, each one of the abovementioned methods is used differently.21,23) However, according to the existing reports, IF techniques are more commonly employed in research works for determining the fracture toughness of ceramic materials, biomaterials, hard biological tissues, etc..2429)

The IF method for biomaterials and hard biological tissues can be classified into various groups, such as Vickers indentation fracture (VIF) test, cube corner indentation fracture (CCIF) test, Vickers crack opening displacement (VCOD) test, and interface indentation fracture (IIF) test.29) Some of the advantages of the IF method include the requirement of less equipment and raw materials, and the easier and faster preparation of test specimens.24,28,30) In the VIF technique, the KIC value of the specimens is often determined by directly measuring the length of cracks created by the Vickers diamond tip around the indented zone of a sample, and using the relevant governing equations.29) A variety of equations exists in this method; in some papers the high accuracy of this technique for some ceramics31) has been reported and in some other papers the existence of large errors, even up to 48%26) has been reported.

In the present research, nanoindentation and indentation tests have been used to determine the hardness and Young’s modulus values of various specimens of pure and silica-doped TZ-3YB, TZ-PX-242A, and TZ-3YSB-E bioceramics containing 0.1 and 0.2 wt.% of α-alumina and sintered at 1500°C. The KIC values of samples have been evaluated and analyzed by using 15 equations from the VIF method. Moreover, the values obtained by the VIF method have also been compared with those obtained by the single-edge laser-notch beam (SELNB) technique. Although better results can be obtained by adding other additives or changing other parameters (e.g. using the fabrication method or sintering temperature for mechanical properties of 3Y-TZP), the specimens in this research have been investigated because of their extensive applications as dental materials and because of other favorable characteristics, such as improved phase composition (because the tetragonal phases are suitably obtained in the 3Y composition) and aging properties and the positive effect of alumina and silica addition on LTD.32)

2. Experimental Procedure

2.1. Samples preparation

In order to compute and analyze the values of KIC, three types of commercially-available 3Y-TZP powders (made by Tosoh Co., Japan), including TZ-3YB, TZ-PX-242A and TZ-3YSB-E, were utilized to prepare the test specimens. The chemical compositions of these powders and their abbreviated designations used in the present research are listed in Table 1. All these powders contain 3–4 wt.% of acrylic binder; although, to a larger extent, their chemical compositions are slightly different. The average size of powder grains is about 29 nm.32) The specific surface areas for standard-grade (pure) TZ-3YB, slightly higher Y2O3-containing biomedical grade TZ-PX-242A along with small quantities of α-Alumina (~ 0.1 wt.%), and also for slightly lower Y2O3-containing biomedical grade TZ-3YSB-E along with large quantities of α-Alumina (~ 0.2 wt.%) are about 16, 10, and 7 m2/g, respectively. Because of its higher translucency, TZ-PX-242A grade is used nowadays for the preparation of porcelain-free, full-contour fixed partial dentures (FPDs), while TZ-3YSB-E grade needs veneering because of its considerable opaqueness. The powders were dry-pressed for 1 min at a pressure of 150 MPa and then pressed using a CIP machine (ABRA Fluid CO.) at a pressure of 250 MPa. The pressed samples in stainless steel mold of 3.34″ size (in the form of 25 × 5 × 2 mm3 plates) were pre-sintered in air for 2 h. The pre-sintering temperature was 1000°C for 0.1A-TZ and TZ specimens and about 1100°C for the 0.2A-TZ specimen.

Chemical Composition of the 3Y-TZP Powders (wt.%)

The pre-sintered samples were divided into two groups. One group was left untreated and used as the control in the analysis of the KIC variations. The samples of the second group were infiltrated with a silica sol, synthesized in situ using the sol–gel method through the hydrolysis of a sylan-based solution (Dynasylan® 6490 made by Evonik, Germany). In this procedure, the specimens were immersed in a combination of ethanol (C2H5OH) and sylan-based solution. Hydrolysis was performed at room temperature by adding aqueous ammonia (20 wt.%) to the samples drop-by-drop. The samples were infiltrated for 30 min. They were then washed with ethanol, dried at room temperature for 24 h, and pyrolyzed at 700°C for 2 h, at a heating rate of 100°C/h. By weighing the specimens before and after the infiltration-drying-pyrolysis cycle, the amount of silica introduced during the infiltration process was determined. Finally, all the specimens of both groups were sintered for 5 h at 1500°C, and the heating/cooling rate was about 5°C/min.

2.2. Mechanical properties measurement

The densities of the samples were obtained by the Archimedes method and based on the European Standard EN 1389:2003.33) These densities were obtained by considering the theoretical density of ρT = 6.08 g/cm3,32) and by averaging five measurements, the relevant SD was also determined.

For analyzing the diagrams obtained from the nanoindentation test, the graphs of force (P) versus penetration depth (h), during sample loading/unloading, were plotted using Equations (1), (2), and (3).34)

(1) Hv=PmaxAc
(2) E=π2A×dpdh
(3) 1Er=(1-v2)E+(1-vi2)Ei

In these equations, Pmax, Ac and Hv indicate the maximum loading, contact area, and Vickers hardness, respectively. Also, dp/dh represents the material stiffness, and Er denotes the calculated modulus. A is the area of all formed micro-cracks around the indented zone, v is Poisson’s ratio, and index i specifies the indenter.

Young’s modulus and Poisson’s ratio of the diamond indenter are 1141 and 0.07 GPa, respectively; the same values are used in the nanohardness measurement device.34) The maximum load applied during the nanohardness test was 800 mN. The nanoindentation tests were performed by using a nanoindentation XP instrument (Nanoinstruments Innovation Center, MTS systems, TN, USA) and a Berkovich diamond indenter with a tip radius of 50 nm. The instrument was operated by applying the continuous stiffness measurement methodology.35) The maximum depth of penetration into the material was 1300 nm. Before performing the nanoindentation tests, the system was calibrated by using the fused silica standard; all the tests were performed at room temperature. Based on the graphs obtained from the nanoindentation tests, the hardness and Young’s modulus values of the samples were determined as a function of the penetration depth. In addition, the Vickers indentation test was used to measure and compare the hardness values. For this purpose, the samples were polished to 1 μm, and the Vickers hardness tests were carried out by applying a Vickers diamond loading of 100 N for a duration of 15 s. The average hardness values were obtained through 7 measurements. The Vickers microhardness test was carried out under laboratory conditions and by using a Tukon 2100B hardness tester.

The VIF test method were used to measure the KIC values of samples polished to 1 μm. Loadings of 29.42, 49.03, 147.09, 196.13, and 245.15 N were applied in order to find the c/a ratios. After establishing the ranges of c/a < 2.5 and c/a ≥ 2.5, for using the equations of Palmqvist and median crack models, two loads of 49.03 N and 196.13 N were respectively applied for 10 s. The KIC values were computed by using the equations of the curve fitting technique and the loading of 196.13 N. By using an optical microscope (Olympus Imaging Corp., Tokyo, Japan), the diagonal and the length of cracks were measured immediately after unloading. To determine the KIC value of each specimen, the test was repeated 10 times and the average value was obtained from the measurements.

In order to compare the KIC values of the samples obtained by the VIF method, these values are also obtained through other approaches, such as the single-edge laser-notch beam (SELNB) method.3639) In the SELNB method, an ultra-sharp V-notch is made on a specimen by laser pulses, and the specimen is then subjected to a three-point bending (TPB) test. The radius of the notch tip can affect the resulting KIC values.40,41) In order to minimize the effect of the notch radius on the KIC, its sharpness must be less than 2 μm.21) Because of the difficulties involved in machining a notch of the desired length, and to avoid unexpected cracks in brittle materials, infrared ultra-short laser pulses are used to create ultra-sharp V-notches on test samples. The notch depth in all the specimens was between 90 and 130 μm. The most important drawback of the SELNB technique is that, due to the saturation effect of the plasma, longer notches cannot be created by this process. Fig. 1 shows a laser notching process for producing a notch of proper geometry and for avoiding unwanted major damage around the created notch. After notching, the laser-notched samples were tested at room temperature by a TPB machine with an opening of 10 mm and constant displacement velocity of 0.15 mm/min. Fig. 2 shows a schematic of a laser-notched specimen during testing. The test was repeated 6 times for each group of samples, and the average KIC values of these samples were obtained. Equations (4) and (5) were used to calculate the KIC values.42)

Fig. 1

View of a laser-produced V-notch in 0.2A-TZ/SiO2 specimen sintered for 5 h at 1500°C.

Fig. 2

Schematic of a laser-notched sample in the SELNB test.

(4) KIC=3LsPmax2BW32·α12·Y*
(5) Y*=(1.989-0.356β)-(1.217+0.315β)α+(3.212+0.705β)α2-(3.222+0.02β)α3+(1.226-0.015β)α4(1-2α)(1-α)32

Pmax is the applied load at fracture, B is the width of the specimens bar, W is the thickness of the specimen bar, α is the notched depth ( aW), Y* is the stress intensity shape factor, and β=WLs.

3. Results and Discussion

3.1. Relative density and microstructure

Table 2 shows the relative density percentages of various samples (TZ, 0.1A-TZ and 0.2A-TZ) under different conditions. As is observed, the relative density of the 0.1A-TZ specimen in all three pre-sintered, sintered (at 1500ºC for 5 h), and silica-doped states is higher than that of the other two specimens. The gravimetrically determined amounts of silica in the infiltrated series of specimens are 0.15, 0.10, and 0.15 wt.% for the TZ, 0.1A-TZ, and 0.2A-TZ samples, respectively. The smallest amount of infiltrated silica is found in the 0.1A-TZ sample, which can be attributed to its higher density relative to the density of other specimens. Fig. 3 shows the SEM images taken from the thermally etched and polished surfaces of the pure and silica-doped 0.1A-TZ samples. As is observed, the size of the grains in both ceramic grades, under similar sintering conditions, is almost the same, and there is no distinct secondary phase. In Fig. 3(b), which shows the grain pattern of the silica-doped 0.1A-TZ specimen, the shapes of grains have changed slightly and their jagged edges have become rounder; this can be attributed to the silica enrichment of the sample. Silica also plays a major role in reducing the amount of strain within the grains.43) Other research works have reported the formation of amorphous grains in alumina/3Y-TZP and silica-enriched materials through TEM analysis and energy-dispersive X-ray spectroscopy (EDS).32,44)

Relative Densities Values in Percentages of Presintered, Sintered (1500°C, 5 h) Pure and Silica-doped Specimens (mean ± SD)

Fig. 3

SEM images taken from thermally-etched and polished surface of (a) pure 0.1A-TZ sample and (b) silica-doped 0.1A-TZ sample.

3.2. Nanoindentation and indentation test

Figure 4 shows the load–displacement (p–h) curves obtained from the nanoindentation test performed on the specimens. A lack of discontinuity in these graphs indicates that fracture has not occurred during the loading/unloading of test samples and signifies the fact that the mechanical integrity of the specimens has been maintained throughout the nanoindentation test. Thus, the obtained empirical results can be reliably used for the evaluation of the samples’ mechanical properties. Also, the elastoplastic behavior of specimens, which have slightly entered the plastic zone due to the application of mN size loads, is significant. This phenomenon may indicate that samples are neither fractured due to the applied milli-Newton loads, nor are they deformed in a ductile manner. Figs. 5 and 6 respectively show the hardness (H) and Young’s modulus (E) values of pure and silica-doped TZ, 0.1A-TZ, and 0.2A-TZ specimens versus the penetration depth. At a penetration depth of more than 50 nm, the E and H values of all the specimens are almost independent of the penetration depth. As is observed, the E and H values of the 0.1A-TZ samples are greater than those of the other specimens, and the E and H values of all silica-doped samples are higher than those of pure samples. The increase or reduction in the hardness can be attributed to the increase or reduction in the density. Fan et al.45) obtained the hardness and Young’s modulus values of a sample made by slip casting technique and consisting of 60 vol.% alumina and 40 vol.% Y-TZP as 15.6 and 340 GPa, respectively. Samodurova et al.32) determined the hardness values of pure and silica-doped 3Y-TZP and 0.05/0.25 wt.% Al2O3-3Y-TZP specimens in the range of 13.5–15.7 GPa. By examining and comparing the Young’s modulus values of samples, it can be realized that Young’s modulus is a function of the relative density. Since Young’s modulus has a direct relationship to density according to the ASTM C769 standard46) and Eq. (6), by increasing the density of samples, which reduces their porosity and the speed of sound in the samples, Young’s modulus increases as well, and vice versa.

Fig. 4

Load–displacement (P–h) curves obtained from nanoindentation tests for (a) pure TZ, 0.1A-TZ, and 0.2A-TZ samples and (b) silica-doped TZ, 0.1A-TZ, and 0.2A-TZ samples.

Fig. 5

Hardness–displacement (H–h) curves obtained from nanoindentation tests for (a) pure TZ, 0.1A-TZ, and 0.2A-TZ samples and (b) silica-doped TZ, 0.1A-TZ. and 0.2A-TZ samples.

Fig. 6

Young’s modulus–displacement (E–h) curves obtained from nanoindentation tests for (a) pure TZ, 0.1A-TZ, and 0.2A-TZ samples and (b) silica-doped TZ, 0.1A-TZ, and 0.2A-TZ samples.

(6) E=ρv2

In the above equation, E is a sample’s Young’s modulus, ρ is its density, and v is the velocity of sound passing through the sample. The average E and H values obtained from 20 nanoindentation tests for each specimen are listed in Table 3.

Obtained Hardness and Young’s Modulus of TZ, 0.1A-TZ, and 0.2A-TZ: Average Elastic Modulus (E) and Hardness (H) Obtained by Nanoindentation and Vickers Hardness (Hv)

In order to compare nanoindentation hardness test results, the Vickers microhardness values of samples are listed in Table 3. As is observed, the Vickers microhardness of each sample is less than the amount of hardness obtained from the nanoindentation test. This discrepancy is justified considering the different ways of measuring hardness in the nanoindentation and Vickers indentation tests, because contrary to the nanoindentation hardness test values, the Vickers microhardness values are measured with the help of residual imprint upon unloading, which results in errors that are mostly caused by elastic relaxation occurring during or after unloading. This phenomenon is important in ceramic materials, and it can explain the existing difference between the results of the two tests. However, the Vickers microhardness of the samples were determined immediately after their unloading so that the difference between the nanoindentation hardness test and the Vickers microhardness test is somewhat reduced. Fig. 7 shows an image of the trace left from the Vickers microhardness in pure 0.1A-TZ samples in loading of 100 N at 0.01 s after unloading.

Fig. 7

Image of the trace left from Vickers microhardness in pure 0.1A-TZ samples in loading of 100 N at 0.01 s after unloading.

3.3. Fracture toughness assessment

3.3.1. Overview of models for the determination of fracture toughness by VIF method

In general, two models of Palmqvist and median/half penny cracks are considered in the VIF approach.47) In the Palmqvist crack model, the length of the formed cracks is assumed to extend only from the tip end of the indented zone,48) while in the median/half penny crack model, the crack length is assumed to extend radially from the center of the indented zone.49) Also, the Palmqvist model and the median crack model mostly occur at low and high loads, respectively. A simple way to make a distinction between these two models is that in the Palmqvist crack model, after polishing the surface of a tested specimen, the formed crack is usually seen separate from the diamond-shaped surface trace left by the Vickers indenter; however, in the median crack model, the crack is completely attached to the surface trace left by the indenter.50,51) Another way of recognizing these two models is that if c/a < 2.5, the model will be assumed as the Palmqvist crack model, and if c/a ≥ 2.5, the model will be considered as the median crack model.52,53) Fig. 8 shows these two models and their differences.

Fig. 8

Crack traces left in the sample by the Vickers indenter: (a) Palmqvist crack model and (b) median/half penny crack model.

Lawn et al.54) analyzed the elastoplastic behavior of the indented zone in various specimens. They assumed the Palmqvist and median crack models due to the presence of tensile stresses during loading, and obtained Eq. (7).

(7) KIC=αEHvPC3/2

In the above equation, KIC is the fracture toughness of the material, P is the amount of loading by Vickers diamond, C is the length of the surface trace of the median crack measured from the center of the indent, Hv is the Vickers hardness of the sample, E is the Young’s modulus, and α is a numerical constant, which is obtained empirically.

After Evans and Charles, other researchers presented numerous equations for calculating the KIC value. More than 30 of these equations have been presented,55) and some of them have been listed in Table 4. Several equations in this table use the curve fitting technique for the acquired data. This group of equations can be used for both the Palmqvist and median crack models. However, the formulation of this group of equations is also based on the two mentioned models, which still make the calculation of KIC dependent on crack morphology.29)

Equations Used in Calculation of Fracture Toughness of Samples

In Table 4, KIC denotes fracture toughness (MPa.m1/2), P is the amount of loading by Vickers diamond (N), c is the average crack length (m), a is the diagonal half-length of Vickers impression (m), l = ca (m), Hv is the Vickers hardness (GPa), and E is the Young’s modulus of the specimen (GPa), tAve is the average thickness of microcracks formed around the indented zone (mm), and A is the total area of crack traces formed around the indent (mm2). F and y in Equations (11) and (13) from Table 4 are determined as.62)

(8) F=-1.59-0.34x-2.02x2+11.24x3-24.97x4+16.32x5
(9) y=-1.59-0.34x-2.02x2+11.24x3-24.97x4+15.32x5

The value of x in Eqs. (8) and (9) is obtained from:

(10) x=Log(ca)

3.3.2. Fracture toughness of samples

Figure 9 shows the indented surface views of pure and silica-doped TZ, 0.1A-TZ, and 0.2 A-TZ specimens, with crack traces at the corners of the indented zone left by the Vickers indenter. Fig. 10 illustrates the c/a ratios with respect to the loadings applied on the test samples. As is observed, for loads of 29.42 and 49.03 N, this ratio indicates a Palmqvist crack system, and for loads of 147.09, 196.13, and 245.15 N, it signifies a median crack system. Hence, by using a loading of 49.03 N in the Palmqvist crack system and a loading of 196.13 N in the median crack system, more reliable KIC results can be obtained. A loading of 196.13 N has also been considered in calculating the KIC values by the curve fitting technique. Although, the Palmqvist crack pattern is normally created in a hard material, the KIC values of many brittle materials can be computed by employing both models, and the choice of which of these two models to use depends on the amount of loading applied by Vickers indenter.50,51) In many brittle materials, Palmqvist cracks only occur at low magnitudes of loading.65) Also, at high loading magnitudes, Palmqvist cracks only form in fairly tough materials such as WC-Co composites containing more than 6 wt.% CO.53)

Fig. 9

Indented surface views of 3Y-TZP composites: (a) Pure TZ at 49.03 N, (b) Pure 0.1A-TZ at 49.03 N, (c) Pure 0.2A-TZ at 49.03 N, (d) Silica-doped TZ at 196.13 N, (e) Silica-doped 0.1A-TZ at 196.13 N, and (f) Silica-doped 0.2A-TZ at 196.13 N.

Fig. 10

c/a ratios versus applied loads in the Vickers indentation test for various pure and silica-doped TZ, 0.1A-TZ, and 0.2A-TZ specimens.

Figure 11 shows the fracture toughness values obtained by the VIF method (using the 15 equations listed in Table 4) and also by the SELNB approach. It can be seen in Fig. 11 that the variations in KIC obtained from the VIF equations for each specimen are close to each other in many cases, and the difference between them is less than 5%; this difference, in comparison with the results of the SELNB method, is less than 10%. Other research works have maintained that, for brittle materials, the KIC results obtained by the SENB approach are more reliable than those obtained by the VIF method,66) and that the values reported by the SELNB method are more accurate than the values obtained by the SENB approach.67,68) Therefore, in the present work, the results of the SELNB method are considered as the KIC values of samples. Eq. (15) in Table 4 is the most recent of these equations, and by comparing the results obtained from this and other equations and also by the SELNB approach, it is found that Eq. (15) yields fairly accurate results. The accuracy of this equation in determining the KIC values of some other brittle materials have also verified the findings of this research.6974) By analyzing the obtained results, it is realized that by adding 0.1 wt.% alumina to pure TZ, the KIC values increase by about 10%, and that the addition of 0.2 wt.% alumina diminishes the KIC values by about 6%. This trend is also seen in silica-doped samples, with an increase and reduction of KIC values by 12 and 4%, respectively, as the mentioned weight fractions of alumina are added. The maximum and minimum KIC values (6.36 and 5.2 MPa.m1/2) are observed in silica co-doped 0.1A-TZ and pure 0.2A-TZ samples, respectively. These results somewhat match the results of other research works on pure and silica-doped 3Y-TZP and 0.05/0.25 wt.% alumina-3Y-TZP samples.32) According to Eq. (11), presented by Griffith, the KIC and Young’s modulus values have a direct relationship.

Fig. 11

Fracture toughness values obtained from 15 different equations in Table 4 (pertaining to VIF method) and from SELNB test for (a) pure TZ, 0.1A-TZ, and 0.2A-TZ samples and (b) silica-doped TZ, 0.1A-TZ, and 0.2A-TZ samples.

(11) KIC=2Eγ

In the above equation, E is the Young’s modulus and Γ denotes the surface energy. Also, by examining the Young’s modulus and KIC values of specimens, it is observed that the trends of these parameters match and overlap each other. Therefore, the increase in KIC can be attributed to the increase in density and, consequently, the increase in Young’s modulus.

4. Conclusions

The relative density, hardness, and Young’s modulus values of pure and silica-doped TZ, 0.1A-TZ, and 0.2A-TZ biomaterials have been investigated in this research through nanoindentation and indentation tests. The fracture toughness values of these materials have been computed and analyzed by means of 15 equations related to the VIF method and compared with those obtained by the SELNB approach. The findings of this research can be summarized as follows:

  • The relative density, hardness, Young’s modulus, and fracture toughness of 3Y-TZP biomaterials are improved with the addition of 0.1 wt.% of alumina. These properties diminish by adding an extra amount of alumina (up to 0.2 wt.%). Also, the mentioned properties are more improved in silica-doped samples than in pure undoped samples. The SEM images of some specimens indicate that there is no distinct secondary phase of alumina present and that the grains of silica-doped samples are often rounder, which can be attributed to the silica enrichment of the samples.

  • The maximum relative density, hardness, Young’s modulus, and fracture toughness values of the samples are 99.8%, 17.3 GPa, 282 GPa, and 6.36 MPa.m1/2, respectively, which have been obtained for the silica-doped 0.1A-TZ. Also, by adding 0.1 wt.% of alumina, the hardness and the Young’s modulus of pure 3Y-TZP specimen increased by 6% and 26%, respectively. Moreover, those of the silica-doped 3Y-TZP specimen increased by 12% and 13%, respectively, at the sintering temperature of 1500°C.

  • In computing the fracture toughness values of all the pure and silica-doped Alumina/3Y-TZP samples by the VIF method, c/a < 2.5 for loadings of 29.42 and 49.03 N. Therefore, it is more appropriate to use the Palmqvist crack model equations to calculate the fracture toughness. For loadings of 147.09, 196.13, and 245.15 N, c/a ≥ 2.5, and therefore, it is more appropriate to use the median crack model equations to obtain the fracture toughness values.

  • By comparing the fracture toughness results obtained using the SELNB and VIF methods, it is found that the values obtained by VIF method, by taking proper loadings into consideration, have a discrepancy of less than 10% relative to the results using SELNB approach. Also, the difference between the results of various VIF equations for each sample is less than 5%.

  • By adding 0.1 wt.% alumina to pure 3Y-TZP material, its fracture toughness is increased by about 10%. Also, the addition of 0.2 wt.% alumina reduced the fracture toughness by about 6%. This trend is also observed more prominently in silica-doped samples, which showed the increase and reduction of fracture toughness by 12 and 4%, respectively, as the mentioned weight fractions of alumina were added.

References

1. Elshazly ES, El-Hout SM, El-Sayed Ali M. Yttria Tetragonal Zirconia Biomaterials: Kinetic Investigation. J Mater Sci Technol 27(4):332–37. 2011;
2. Guazzato M, Albakry M, Ringer SP, Swain MV. Strength, Fracture Toughness and Microstructure of a Selection of All-Ceramic Materials: Part II. Zirconia-Based Dental Ceramics. Dent Mater 20(5):4449–56. 2004;
3. Matsui K, Yamakawa T, Uehara M, Enomoto N, Hojo J. Mechanism of Alumina-Enhanced Sintering of Fine Zirconia Powder: Influence of Alumina Concentration on the Initial Stage Sintering. J Am Ceram Soc 91(6):1888–97. 2008;
4. Park MJ, Yang SK, Kang JB. Effects of Composition and Additives on the Mechanical Characteristics of 3Y-TZP. J Korean Ceram Soc 43(10):640–45. 2006;
5. Yang SK, Bae KM, Cho BR, Kang JB. Effect on Mechanical Properties of 3Y-TZP; (I) Addition of Monoclinic Zirconia. J Korean Ceram Soc 42(6):411–16. 2005;
6. Yilmaz H, Aydin C, Gul BE. Flexural Strength and Fracture Toughness of Dental Core Ceramics. J Prosthet Dent 98(2):120–28. 2007;
7. Denry I, Kelly JR. State of the Art of Zirconia for Dental Applications. Dent Mater 24(3):299–307. 2008;
8. Egilmez F, Ergun G, Cekic-Nagas I, Vallittu PK, Lassila LV. Factors Affecting the Mechanical Behavior of Y-TZP. J Mech Behave Biomed Mater 37:78–87. 2014;
9. Tekeli S. Fracture Toughness (KIC ), Hardness, Sintering and Grain Growth Behavior of 8YSCZ/Al2O3 Composites Produced by Colloidal Processing. J Alloys Compd 391(1–2):217–24. 2005;
10. Piconi C, Maccauro G. Zirconia as a Ceramic Biomaterial. Biomaterials 20(1):1–25. 1999;
11. Chevalier J, Gremillard L, Deville S. Low-Temperature Degradation of Zirconia and Implications for Biomedical Implants. Annu Rev Mater Res 37:1–32. 2007;
12. Camposilvan E, Garcia Marro F, Mestra A, Anglada M. Enhanced Reliability of Yttria-Stabilized Zirconia for Dental Applications. Acta Biomater 17:36–46. 2015;
13. Majić Renjo M, Ćurković L, Štefančić S, Ćorić D. Indentation Size Effect of Y-TZP Dental Ceramics. Dent Mater 30(12):371–76. 2014;
14. Harada K, Shinya A, Yokoyama D, Shinya A. Effect of Loading Conditions on the Fracture Toughness of Zirconia. J Prosthodont Res 57(2):82–87. 2013;
15. Sailer I, Gottnerb J, Kanelb S, Hammerle CH. Randomized Controlled Clinical Trial of Zirconia-Ceramic and Metal-Ceramic Posterior Fixed Dental Prostheses: a 3-year Follow-Up. Int J Prosthodont 22(6):553–60. 2009;
16. Kobayashi K, Kuwajima H, Masaki T. Phase Change and Mechanical Properties of ZrO2–Y2O3 Solid Electrolyte after Ageing. Solid State Ionics 3–4:489–95. 1981;
17. Kim HT, Han JS, Yang JH, Lee JB, Kim SH. The Effect of Low Temperature Aging on the Mechanical Property & Phase Stability of YTZP Ceramics. J Adv Prosthodont 1(3):113–17. 2009;
18. Pereira GKR, Venturini AB, Silvestri T, Dapieve KS, Montagner AF, Soares FZM, Valandro LF. Lowtemperature Degradation of Y-TZP Ceramics: A Systematic Review and Meta-Analysis. J Mech Behav Biomed Mater 55:151–63. 2015;
19. Freimann SW. Brittle Fracture Behavior of Ceramics. Am Ceram Soc Bull 67(2):392–402. 1988;
20. Pabst RF, Kromp K, Popp G. Fracture Toughness-Measurement and Interpretation. Proc Br Ceram Soc 32:89–94. 1982;
21. Fischer H, Waindich A, Telle R. Influence of Preparation of Ceramic SEVNB Specimens on Fracture Toughness Testing Results. Dent Mater 24:618–22. 2008;
22. Kailer A, Marc S. On the Feasibility of the Chevron Notch Beam Method to Measure Fracture Toughness of Fine-Grained Zirconia Ceramics. Dent Mater 32(10):1256–62. 2016;
23. Gogotsi GA. Fracture Toughness of Ceramics and Ceramic Composites. Ceram Int 29:777–84. 2003;
24. Denry IL, Holloway JA. Elastic Constants, Vickers Hardness, and Fracture Toughness of Fluorrichterite-Based Glass–Ceramics. Dent Mater 20:213–19. 2004;
25. Garrido MA, Giráldez I, Ceballos L, Rodríguez J. On the Possibility of Estimating the Fracture Toughness of Enamel. Dent Mater 30:1224–33. 2014;
26. Fischer H, Marx R. Fracture Toughness of Dental Ceramics: Comparison of Bending and Indentation Method. Dent Mater 18:12–19. 2002;
27. Mahoney E, Holt A, Swain M, Kilpatrick N. The Hardness and Modulus of Elasticity of Primary Molar Teeth: an Ultra-Microindentation Study. J Dent 28:589–94. 2000;
28. Şakar-Deliormanli A, Güden M. Microhardness and Fracture Toughness of Dental Materials by Indentation Method. J Biomed Mater Res, Part B 76(2):257–64. 2006;
29. Kruzic JJ, Kim DK, Koester KJ, Ritchie RO. Indentation Techniques for Evaluating the Fracture Toughness of Biomaterials and Hard Tissues. J Mech Behav Biomed Mater 2(4):384–95. 2009;
30. Mullins LP, Bruzzi MS, McHugh PE. Measurement of the Microstructural Fracture Toughness of Cortical Bone Using Indentation Fracture. J Biomech 40(14):3285–88. 2007;
31. Moradkhani A, Baharvandi H, Tajdari M, Latifi H, Martikainen J. Determination of Fracture Toughness Using the Area of Micro-Crack Tracks Left in Brittle Materials by Vickers Indentation Test. J Adv Ceram 2:87–102. 2013;
32. Samodurova A, Kocjan A, Swain MV, Kosmač T. The Combined Effect of Alumina and Silica Co-Doping on the Ageing Resistance of 3Y-TZP Bioceramics. Acta Biomater 11:477–87. 2015;
33. European Standard DIN EN 1389:2003. Advanced Technical Ceramics, Ceramic Composites, Physical Properties-Determination of Density and Apparent Porosity (BSI 81.060.30 Publication, 1994; https://shop.bsigroup.com/en/ProductDetail/?pid=000000000000345631&_ga=2.10794420.1058334025.1539768367-258163731.1539768367).
34. Oliver WC, Pharr GM. An Improved Technique for Determining Hardness and Elastic Modulus Using Load and Displacement Sensing Indentation Experiment. J Mater Res 7(6):1564–83. 1992;
35. Pharr G. Measurement of Mechanical Properties by Ultra-Low Load Indentation. Mater Sci Eng, A 253(1–2):151–59. 1998;
36. Carlton HD, Elmer JW, Freeman DC, Schaeffer RD, Derkach O, Gallegos GF. Laser Notching Ceramics for Reliable Fracture Toughness Testing. J Eur Ceram Soc 36(1):227–34. 2016;
37. Pastor JY. How to Measure the Real Fracture Toughness in Brittle Materials? Past, Present and Future Techniques. In : Proceeding of the International Conference on Experimental Mechanics (ICEM15). Porto, Portugal; 2012;
38. Palacios T, Pastor JY. Influence of the Notch Root Radius on the Fracture Toughness of Brittle Metals: Nanostructure Tungsten Alloy, a Case Study. Int J Refract Met Hard Mater 52:44–9. 2015;
39. Pastor JY, LLorca J, Martín A, Peña JI, Oliete PB. Fracture Toughness and Strength of Al2O3–Y3Al5O12 and Al2O3–Y3Al5O12–ZrO2 Directionally Solidified Eutectic Oxides up to 1900K. J Eur Ceram Soc 28(12):2345–51. 2008;
40. Damani R, Gstrein R, Danzer R. Critical Notch-Root Radius Effect in SENB-S Fracture Toughness Testing. J Eur Ceram Soc 16(7):695–702. 1996;
41. Nishida T, Hanaki Y, Pezzotti G. Effect of Notch-Root Radius on the Fracture Toughness of a Fine-Grained Alumina. J Am Ceram Soc 77(2):606–8. 1994;
42. Guinea GV, Pastor JY, Planas J, Elices M. Stress Intensity Factor, Compliance and CMOD for a General Three-Point-Bend Beam. Int J Fract 89(2):103–16. 1998;
43. Gremillard L, Chevalier J, Epicier T, Fantozzi G. Improving the Durability of a Biomedical-Grade Zirconia Ceramic by the Addition of Silica. J Am Ceram Soc 85(2):401–7. 2002;
44. Gremillard L, Epicier T, Chevalier J, Fantozzi G. Microstructural Study of Silica-Doped Zirconia Ceramics. Acta Mater 48(18–19):4647–52. 2000;
45. Fan K, Pastor JY, Ruiz-Hervias J, Gurauskis J, Baudin C. Determination of Mechanical Properties of Al2O3/Y-TZP Ceramic Composites: Influence of Testing Method and Residual Stresses. Ceram Int 42(16):18700–10. 2016;
46. ASTM C769-98. Standard Test Method for Sonic Velocity in Manufactured Carbon and Graphite Materials for Use in Obtaining an Approximate Young’s Modulus ASTM International. West Conshohocken, PA: 1998.
47. Nastic A, Merati A, Bielawski M, Bolduc M, Fakolujo O, Nganbe M. Instrumented and Vickers Indentation for the Characterization of Stiffness, Hardness and Toughness of Zirconia Toughened Al2O3 and SiC Armor. J Mater Sci Technol 31(8):773–83. 2015;
48. Palmqvist S. The Work for the Formation of a Crack during Vickers Indentation as a Measure of the Toughness of Hard Metals. Arch Eisenhuettenwes 33:629–34. 1962;
49. Anstis GR, Chantikul P, Lawn BR, Marshall DB. A Critical Evaluation of Indentation Techniques for Measuring Fracture Toughness: I, Direct Crack Measurements. J Am Ceram Soc 64(9):533–38. 1981;
50. Barsoum MW. Fundamentals of Ceramics p. 368–69. Taylor & Francis Group. New York: 2003.
51. Lin JD, Duh JG. Fracture Toughness and Hardness of Ceria and Yttria-Doped Tetragonal Zirconia Ceramics. Mater Chem Phys 78(1):253–61. 2002;
52. Bamzai KK, Kotru PN, Wanklyn BM. Fracture Mechanics, Crack Propagation and Microhardness Studies on Flux Grown ErAlO3 Single Crystals. J Mater Sci Technol 16(04):405–10. 2000;
53. Bhat M, Kaur B, Kumar R, Bamzai KK, Kotru PN, Wanklyn BM. Effect of Ion Irradiation on Dielectric and Mechanical Characteristics of ErFeO3 Single Crystals. Nucl Instrum Methods Phys Res, Sect B 234(4):494–508. 2005;
54. Lawn BR, Evans AG, Marshall DB. Elastic/Plastic Indentation Damage in Ceramics: The Median/Radial Crack System. J Am Ceram Soc 63(9–10):574–81. 1980;
55. Quinn GD, Bradt RC. On the Vickers Indentation Fracture Toughness Test. J Am Ceram Soc 90(3):673–80. 2007;
56. Niihara K. A Fracture Mechanics Analysis of Indentation-Induced Palmqvist Crack in Ceramic. J Mater Sci Lett 2(5):221–23. 1983;
57. Lawn HR, Fuller ER. Equilibrium Penny-Like Cracks in Indentation Fracture. J Mater Sci 10(12):2016–24. 1975;
58. Evans AG, Wilshaw TR. Quasi-Static Solid Particle Damage in Brittle Solid-I. Observation Analysis and Implications. Acta Metall 24(10):939–56. 1976;
59. Laugier MT. New formula for indentation toughness in ceramics. J Mater Sci Lett 6:355–6. 1987;
60. Evans AG, Charles EA. Fracture Toughness Determinations by Indentation. J Am Ceram Soc 59(7–8):371–72. 1976;
61. Shetty DK, Wright IG, Mincer PN, Cluar AH. Indentation Fracture of WC-Co Cermets. J Mater Sci 20(5):1873–82. 1985;
62. Evans AG. Fracture Toughness: the Role of Indentation Techniques ASTM Special Technical Publication; 1979.
63. JIS R. 1607. Testing methods for fracture toughness of high performance ceramics Japanese Standard Association; 1990.
64. Lankford J. Indentation Microfracture in the Palmqvist Crack Regime: Implications for Fracture Toughness Evaluation by the Indentation Method. J Mater Sci Lett 1(11):493–95. 1982;
65. Niihara K, Morena R, Hasselman DPH. Evaluation of KIC of Brittle Solids by the Indentation Method with Low Crack-to-Indent Ratios. J Mater Sci Lett 1(1):13–6. 1982;
66. Xu J, Tang D, Lee KJ, Lim HB, Park K-S, Cho W. Comparison of Fracture Toughness Evaluating Methods in 3Y-TZP Ceramics Reinforced with Al2O3 Particles. J Ceram Process Res 13(6):83–7. 2012;
67. Tuan WH, Chen RZ, Wang TC, Cheng CH, Kuo PS. Mechanical Properties of Al2O3/ZrO2 Composites. J Eur Ceram Soc 22(16):2827–33. 2002;
68. Szutkowska M. Fracture Resistance Behavior of Alumina–Zirconia Composites. J Mater Process Technol 153(1):868–74. 2004;
69. Moradkhani A, Baharvandi H, Samani MMM. Mechanical Properties and Microstructure of B4C–Nano-TiB2–Fe/Ni Composites under Different Sintering Temperatures. Mater Sci Eng A 665:141–53. 2016;
70. Latifi H, Moradkhani A, Baharvandi H, Martikainen J. Fracture Toughness Determination and Microstructure Investigation of a B4C–NanoTiB2 Composite with Various Volume Percent of Fe and Ni Additives. Mater Des 62:392–400. 2014;
71. Moradkhani A, Baharvandi H. Analyzing the Microstructures of W-ZrC Composites Fabricated through Reaction Sintering and Determining their Fracture Toughness Values by Using the SENB and VIF Methods. Eng Fract Mech 189:501–13. 2018;
72. Moradkhani A, Baharvandi H. Effects of Additive Amount, Testing Method, Fabrication Process and Sintering Temperature on the Mechanical Properties of Al2O3/3Y-TZP Composites. Eng Fract Mech 191:446–60. 2018;
73. Moradkhani A, Baharvandi H. Mechanical Properties and Fracture Behavior of B4C-Nano/Micro SiC Composites Produced by Pressureless Sintering. Int J Refract Met Hard Mater 70:107–15. 2018;
74. Moradkhani A, Baharvandi H. Determining the Fracture Resistance of B4C-NanoSiB6 Nanocomposite by Vickers Indentation Method and Exploring its Mechanical Properties. Int J Refract Met Hard Mater 68:159–65. 2017;

Article information Continued

Fig. 1

View of a laser-produced V-notch in 0.2A-TZ/SiO2 specimen sintered for 5 h at 1500°C.

Fig. 2

Schematic of a laser-notched sample in the SELNB test.

Fig. 3

SEM images taken from thermally-etched and polished surface of (a) pure 0.1A-TZ sample and (b) silica-doped 0.1A-TZ sample.

Fig. 4

Load–displacement (P–h) curves obtained from nanoindentation tests for (a) pure TZ, 0.1A-TZ, and 0.2A-TZ samples and (b) silica-doped TZ, 0.1A-TZ, and 0.2A-TZ samples.

Fig. 5

Hardness–displacement (H–h) curves obtained from nanoindentation tests for (a) pure TZ, 0.1A-TZ, and 0.2A-TZ samples and (b) silica-doped TZ, 0.1A-TZ. and 0.2A-TZ samples.

Fig. 6

Young’s modulus–displacement (E–h) curves obtained from nanoindentation tests for (a) pure TZ, 0.1A-TZ, and 0.2A-TZ samples and (b) silica-doped TZ, 0.1A-TZ, and 0.2A-TZ samples.

Fig. 7

Image of the trace left from Vickers microhardness in pure 0.1A-TZ samples in loading of 100 N at 0.01 s after unloading.

Fig. 8

Crack traces left in the sample by the Vickers indenter: (a) Palmqvist crack model and (b) median/half penny crack model.

Fig. 9

Indented surface views of 3Y-TZP composites: (a) Pure TZ at 49.03 N, (b) Pure 0.1A-TZ at 49.03 N, (c) Pure 0.2A-TZ at 49.03 N, (d) Silica-doped TZ at 196.13 N, (e) Silica-doped 0.1A-TZ at 196.13 N, and (f) Silica-doped 0.2A-TZ at 196.13 N.

Fig. 10

c/a ratios versus applied loads in the Vickers indentation test for various pure and silica-doped TZ, 0.1A-TZ, and 0.2A-TZ specimens.

Fig. 11

Fracture toughness values obtained from 15 different equations in Table 4 (pertaining to VIF method) and from SELNB test for (a) pure TZ, 0.1A-TZ, and 0.2A-TZ samples and (b) silica-doped TZ, 0.1A-TZ, and 0.2A-TZ samples.

Table 1

Chemical Composition of the 3Y-TZP Powders (wt.%)

Powder grade Abbreviations Y2O3 α-Al2O3 SiO2 Fe2O3 Na2O Specific surface area (m2/g)
TZ-3YB TZ 5.27 < 0.01 < 0.005 < 0.005 < 0.02 ~ 16
TZ-PX-242A 0.1A-TZ 5.34 0.1 < 0.005 < 0.005 < 0.01 ~ 10
TZ-3YSB-E 0.2A-TZ 5.18 0.2 < 0.005 < 0.005 < 0.01 ~ 7

Table 2

Relative Densities Values in Percentages of Presintered, Sintered (1500°C, 5 h) Pure and Silica-doped Specimens (mean ± SD)

TZ 0.1A-TZ 0.2A-TZ
Pre-sintered 52.3±0.3 55.6±0.4 53.1±0.5
As sintered 99.3±0.2 99.6±0.1 99.0±0.2
Silica-doped/sintered 99.4±0.2 99.8±0.1 98.9±0.1

Table 3

Obtained Hardness and Young’s Modulus of TZ, 0.1A-TZ, and 0.2A-TZ: Average Elastic Modulus (E) and Hardness (H) Obtained by Nanoindentation and Vickers Hardness (Hv)

Materials Nanoindentation Vickers indentation

H (GPa) E (GPa) Hv (GPa)
Pure TZ 14.3±0.6 215±7 11.7±0.1
0.1A-TZ 15.2±0.4 272±6 12.9±0.2
0.2A-TZ 13.4±0.4 193±4 11.5±0.3

Silica-doped TZ 15.4±0.7 250±5 12.5±0.3
0.1A-TZ 17.3±0.6 282±5 14.1±0.3
0.2A-TZ 14.5±0.5 235±8 12 ±0.2

Table 4

Equations Used in Calculation of Fracture Toughness of Samples

No. Equation Crack type Author Ref.
1 KIC=0.018Hva1/2(EHv)0.4(ca-l)-.05 Palmqvist Niihara et al. [56]
2 KIC=0.0515(Pc3/2) Palmqvist Lawn and Fuller [57]
3 KIC=0.079(Pa3/2)Log(4.5ac) Palmqvist Evans and Wilshaw [58]
4 KIC=0.015(EHv)2/3(Pc3/2)(la)-0.5 Palmqvist Laugier [59]
5 KIC=0.016(EHv)1/2PC3/2 Median Anstis et al. [49]
6 KIC=0.0752Pc3/2 Median Evans and Charles [60]
7 KIC=0.0726Pc3/2 Median Lawn and Fuller [57]
8 KIC=0.014(EHv)1/2PC3/2 Median Lawn and Evans [54]
9 KIC=0.16Hva1/2(ca)3/2 Median Evans and Charles [60]
10 KIC=0.0889(Hv.Pi=14ci)1/2 Curve fitting technique Shetty et al. [61]
11 KIC=0.4636(EHv)2/5Pa3/210F Curve fitting technique Evans [62]
12 KIC=0.018(EHv)1/2PC3/2 Curve fitting technique Japanese Standards Association [63]
13 KIC=(Hva1/2)(EHv)2/510y Curve fitting technique Evans [62]
14 KIC=0.0782(Hva1/2)(EHv)25(Ca)-1.56 Curve fitting technique Lankford [64]
15 KIC=0.00366(EHv)12tAve3/2PA3/2 Curve fitting technique Moradkhani et al. [31]