CHAPTER 7
Investigation of Elastic-Plastic Transitional Stresses in Zirconia-Based Ceramic Dental Implants under Uniaxial Compression
1Research Scholar, Department of Mathematics, Punjabi University, India, E-mail: shivdevshahi93@gmail.com
2Professor, Department of Mathematics, Punjabi University, Patiala, India, E-mail: sbsingh69@yahoo.com
3Professor Emeritus, Canadian Research and Development Center of Sciences and Cultures, Montreal, Canada
7.1 INTRODUCTION
Zirconia (ZrO2) has been considered as a wear-resistant and osteoconductive ceramic, highly suitable for making stress-bearing crowns of dental implants. Various properties, namely, less inflammatory response, compatibility of zirconia in contact with gum tissues, lower plaque retention, and highly aesthetic tooth-like coloration has made zirconia a suitable alternative to titanium implants [2, 6, 19]. The production of ceramic dental implants have improved in past few years [7, 20]. The yttria-stabilized tetragonal zirconia polycrystal, popularly known as Y-TZP has been a preference for some years now, having the highest fracture strength [10, 12, 15, 27]. Muhammad et al. [11] confirmed in their chapter that zirconia is highly flexible in all directions of the lattice plane. The Poisson ratio obtained, i.e., from 0.16 to 0.31 also indicate high anisotropy. The stiffness constants so obtained will also be used in this chapter for the determination of elastic-plastic transitional stresses.
This chapter is concerned with the investigation of elastic-plastic transition stress in Zirconia-based crowns of ceramic dental implants, further comparing the results with stresses in titanium-based implants. The chapter also checks the similarities between the implants and human tooth enamel using Seth’s transition theory. The crowns will be modeled as spherical shells exhibiting transversely isotropic macrostructural symmetry. Equations for modeling spherical shells made of isotropic materials are available in most standard textbooks [3, 9, 13, 14, 17, 29]. Following the classical methods of stress and strain determination, Miller [16] evaluated solutions for stresses and displacements in a thick spherical shell subjected to internal and external pressure loads. You et al. [30] presented a highly precise model to carry out elastic analysis of thick-walled spherical pressure vessels. The authors have studied the behavior of shells particularly when some assumptions, such as (i) incompressibility of material used, (ii) creep strain law derived by Norton, (iii) yield condition of Tresca, and (iv) associated flow rules; were made. The need of utilization of these specially appointed semi-experimental laws in elastic-plastic transition depends on approach that the transition is a linear phenomenon which is unrealistic. Deformation fields related with the irreversible phenomenon, such as elastic-plastic disfigurements, creep relaxation, fatigue, and crack, etc. are non-linear in character. The traditional measures of deformation are not adequate to manage transitions. The concept of generalized strain measures and transition theory given by Seth [21] has been applied to find elastic-plastic stresses in various problems by solving the non-linear differential equations at the transition points. Shahi and Singh [24] successfully calculated elastic-plastic transitional stresses in human tooth enamel and dentine using this theory. The results for hydroxyapatite (HAP)-Ca10(PO4)6(OH)2, which structures 95% of enamel and 50% of dentine by weight, were also obtained for comparative analysis. The theory has been used to solve various problems of stress and strain determination in structures modeled in the form of discs and shells [25, 28]. All these problems based on the recognition of the transition state as separate state necessitates showing the existence of the used constitutive equation for that state.
7.2 GOVERNING EQUATIONS
We consider a spherical shell of constant thickness with internal and external radius a and b respectively under external pressure p. The external pressure will act radially to simulate the state of axial compression (Figure 7.1).

FIGURE 7.1 Structure of the crown part of the implant.
1. Displacement Coordinates: The components of displacement in spherical coordinates (r, θ, ø) are taken as:
![]() | (7.1) |
where: ß is position function depending on r.
The generalized components of strain are given by Seth [21, 22] as:>
![]() | (7.2) |
where: n is the measure and rαβ’ = βP; P is a function of ß and ß is a function of r.
2. Stress-Strain Relation: The stress-strain relations for isotropic material are given by Sokolinokoff [26]:
![]() |
where Tij and ekl are the stress and strain tensors respectively. These nine equations contain a total of 81 coefficients eijkl, but not all the coefficients are independent. The symmetry of Tij and eij reduces the number of independent coefficients to 36. For Transversely isotropic materials which have a plane of elastic symmetry, these independent coefficients reduce to 5. The constitutive equations for transversely isotropic media are given by Altenbach et al. [1]:
![]() | (7.3) |
Substituting (Eqn. 2) in (Eqn. 3), we get:
![]() | (7.4) |
![]() |
3. Equation of Equilibrium: The equations of equilibrium are:
![]() |
![]() | (7.5) |
Substituting (Eqn. 4) in (Eqn. 5), we see that the equations of equilibrium are all satisfied except:
![]() | (7.6) |
![]() | (7.7) |
From (Eqn. 7), one may also say that:
![]() | (7.8) |
(Eqn. 8) is satisfied by Tθθ and Tϕϕ as given by (Eqn. 2). If c21 = c31, c22– c33– c32– c23 the equation of equilibrium from (Eqn. 6) becomes:
![]() | (7.9) |
4. Critical Points or Turning Points: By substituting (Eqn. 4) into (Eqn. 9), we get a non-linear differential equation in terms of β:
![]() | (7.10) |
where;
![]() |
and
![]() |
where P is function of ß and ß is function of r only.
5. Transition Points: The transition points of ß in (Eqn. 10) are P = 0, P → –1 and P → ± ∞
To solve the Elastoplastic stress problems we consider the case of P → ± ∞
6. Boundary Condition: The boundary conditions of the problem are given by:
![]() | (7.11) |
7.3 PROBLEM SOLUTION
For finding the elastic-plastic stresses, the transition function is taken through the principal stresses at the transition point P → ± ∞, we define the transition function ζ as:
![]() | (7.12) |
where ζ be the transition function unction of r only. Taking the logarithmic differentiation of (Eqn. 12), with respect to r and using (Eqn. 10), we get:
![]() | (7.13) |
Taking the asymptotic value of (Eqn. 13) as P → ± ∞ and integrating, we get:
![]() | (7.14) |
where A is a constant of integration and C1 = (c33 – c13)/c33. From (Eqn. 12) and (7.14), we have:
![]() | (7.15) |
Using boundary condition from (Eqn. 11) in (Eqn. 15), we get:
![]() | (7.16) |
Substituting (Eqn. 16) in to (Eqn. 15) and using (Eqn. 16) in equation of equilibrium, we get:
![]() | (7.17) |
1. Initial Yielding: From (Eqn. 17), it is seen that |Tθθ–Trr| is maximum at the outer surface (that is at r = b), therefore yielding of the shell will take place at the external surface of the shell:
![]() | (7.18) |
Using (Eqn. 18) in Eqns. (7.15)–(7.17), we get the transitional stresses as in non-dimensional components as:
![]() | (7.19) |
where;
![]() |
and
![]() |
2. Fully-Plastic State: For fully-plastic case [23], C1 → 0; therefore stresses and pressure from (Eqn. 19) becomes:
![]() | (7.20) |
where,
![]() |
and
![]() |
7.4 NUMERICAL RESULTS AND DISCUSSION
The above investigations elaborate the initial yielding and fully plastic state of a crown made of zirconia and titanium modeled in the form of spherical shell subjected to external pressure, to analyze uniaxial compression. The elastic constants for the same are taken from the literature [11] which has been obtained by ultrasonic resonance spectroscopy, a nondestructive measure to obtain the stiffness constants. The results obtained for both types of crowns are compared with enamel. Enamel is made up of HAP mineral, 95% by vol. All these materials exhibit transversely isotropic macrostructural symmetry.
In Figure 7.2, the curves are plotted for pressure at initial yielding at various radius ratios. It is observed that the intensity of pressure at initial yielding increases with an increase in the thickness of the shell. It has been observed that titanium had the lowest yield strength and yielded at lower levels of stress as compared to enamel and Zirconia. Figures 7.3 and 7.4 show the trends of radial and circumferential stresses at initial yielding. Maximum stresses were observed at the external surface of the shell. In Figure 7.5, the curves are plotted for the pressure required at a fully plastic state for various radius ratios. It has been observed that shells exhibited high plasticity when the thickness of the shell was between ratios 1 < R0 < 3, particularly Zirconia. A significant drop in the levels of plasticity is observed when the thickness increases. The plasticity of Zirconia is inferred to be greater than that of titanium. Figures 7.6 and 7.7 represent the trends of radial and circumferential stresses at a fully plastic state. The observations infer to the fact that the principal stress differences were maximum at the external surface of the crowns.

FIGURE 7.2 Pressure in the Shell at initial yielding for radius ratio R0.

FIGURE 7.3 Radial stresses at initial yielding.

FIGURE 7.4 Circumferential stresses at initial yielding.

FIGURE 7.5 Pressure in the shell at fully plastic state for radius ratio R0.

FIGURE 7.6 Radial stresses at fully plastic state.

FIGURE 7.7 Circumferential stresses at fully plastic state.
7.5 CONCLUSIONS
The findings allow us to conclude that zirconia has a greater resemblance with enamel with the necessary elastic and plastic limit, which demonstrates considerable ability to suppress a crack growth. Varying values of pressure required for initial yielding and fully plastic state were calculated for various radius ratios depending on the geometry of the crown sample. Trends of the graphs were similar for enamel and HAP due to enamel’s composition. The significant difference between stress buildup at the inner and outer layer of the implant crown is observed by varying the radii ratios.
KEYWORDS
• boundary condition
• critical points
• equilibrium
• hydroxyapatite
• initial yielding
• stress-strain relation
REFERENCES
1. Altenbach, H., Altenbach, J., & Kissing, W., (2004). Mechanics of Composite Structural Elements. Springer-Verlag.
2. Bianchi, A. E., Bosetti, M., Dolci, G., Sberna, M. T., Sanfilippo, F., & Cannas, M. (2004). In vitro and in vivo follow-up of titanium transmucosal implants with zirconia collar. J. Appl. Biomater. Biomech., 2, 143–150.
3. Boyle, J. T., & Spence, J., (1983). Stress Analysis for Creep. Butterworth’s Coy. Ltd. London.
4. Cuy, J. L., Mann, A. B., Livi, K. J., Teaford, M. F., & Weihs, T. P., (2002). Nanoindentation mapping of the mechanical properties of human molar tooth enamel. Arch. Oral. Biol., 47, 281–291.
5. Zaytsev, D., Grigoriev, S., & Panfilov, P., (2012). Deformation behavior of human dentin under uniaxial compression. Int. J. Biomater, 8. Article ID: 854539.
6. Degidi, M. Artese, L. Scarano, A. Perrotti, V. Gehrke, P. & Piattelli, A. (2006). Inflammatory infiltrate, microvessel density, nitric oxide synthase expression, and proliferative activity in peri-implant soft tissues around titanium oxide healing caps. J. Periodontol, 77, 73–80.
7. Depprich, R., Ommerborn, M., Zipprich, H., Naujoks, C., Handschel, J., & Wiesmann, H. P., (2008). Behavior of osteoblastic cells cultured on titanium and structured zirconia surfaces. Head Face Med., 4, 29.
8. Menéndez-Proupina, E., Cervantes-Rodríguezb, S., Osorio-Pulgara, R., Franco-Cisternaa, M., Camacho-Montesc, H., & Fuentesb, M. E., (2011). Computer simulation of elastic constants of hydroxyapatite and fluorapatite. Journal of the Mechanical Behavior of Biomedical Materials, 4, 1011–1020.
9. Fung, Y. C., (1965). Foundations of Solid Mechanics. Englewood Cliffs, N.J. Prentice-Hall.
10. Hallmann, L., Mehl, A., Ulmer, P., Reusser, E., Stadler, J., & Zenobi, R. (2012). The influence of grain size on low-temperature degradation of dental zirconia. J. Biomed. Mater. Res. B Appl. Biomat., 100, 447–456.
11. Muhammad, I. D., et al., (2014). Modeling the elastic constants of cubic zirconia using molecular dynamics simulations. Advanced Materials Research, 845, 387–391.
12. Inokoshi, M., Zhang, F., De Munck, J., Minakuchi, S., Naert, I., & Vleugels, J. (2014). Influence of sintering conditions on low-temperature degradation of dental zirconi. Dent. Mater., 30, 669–678.
13. Kraus, H., (1980). Creep Analysis. Wiley, New York, Usa.
14. Lubhan, D., & Felger, R. P., (1961). Plasticity and Creep of Metals. Wiley, New York, USA.
15. Lughi, V., & Sergo, V. (2010). Low temperature degradation-aging- of zirconia: A critical review of the relevant aspects in dentistry. Dent. Mater., 26, 807–820.
16. Miller, G. K., (1995). Stresses in a spherical pressure vessel undergoing creep and dimensional changes. International Journal of Solids and Structures, 32, 2077–2093.
17. Parkus, H., (1976). Thermo-Elasticity. Springer-Verlag, Wien, New York.
18. Craig, R. G., Peyton, F. A., & Johnson, D. W., (1961). Compressive properties of enamel, dental cements, and gold. J. Dent. Res., 40, 936–945.
19. Rimondini, L., Cerroni, L., Carrasi, A., & Torricelli, P. (2002). Bacterial colonization of zirconia ceramic surfaces: An in vitro and in vivo study. Int. J. Oral Maxillofac Implants, 17, 793–798.
20. Scarano, A., Piatelli, M., Caputi, S., Favero, G. A., & Piattelli, A. (2004). Bacterial adhesion on commercially pure titanium and zirconium oxide disks: An in vivo human study. J. Periodontol., 75, 292–296.
21. Seth, B. R., (1962). Transition theory of elastic-plastic deformation, creep, and relaxation. Nature, 195, 896–897.
22. Seth, B. R., (1966). Measure concept in mechanics. International Journal of Non-Linear Mechanics, 1(1), 35–40.
23. Seth, B. R., (1963). Elastic-plastic transition in shells and tubes under pressure. Z. Angew. Math. Mech., 43, 345–351.
24. Shahi, S. & Singh, S. B. (2020). Elastic plastic transitional stress analysis of human tooth enamel and dentine under external pressure using Seth’s transition theory, materials physics, and chemistry. Applied Mathematics and Chemo-Mechanical Analysis. AAP, Ch: 3 (in press).
25. Shahi, S., Singh, S. B., & Thakur, P., (2019). Modeling creep parameter in rotating discs with rigid shaft exhibiting transversely isotropic and isotropic material behavior. Journal of Emerging Technologies and Innovative Research, 6(1), 387–395.
26. Sokolinokoff, I. S., (1956). Mathematical Theory of Elasticity (2nd edn.), McGraw-Hill, New York.
27. Stawarczyk, B., Ozcan, M., Hallmann, L., Ender, A., Mehl, A., & Hämmerlet, C. H., (2013). The effect of zirconia sintering temperature on flexural strength, grain size, and contrast ratio. Clin. Oral. Investig., 17, 269–274.
28. Thakur, P., Shahi, S., Gupta, N., & Singh, S. B., (2017). Effect of mechanical load and thickness profile on creep in a rotating disc by using Seth’s transition theory. AIP Conf. Proc., Amer. Inst. of Physics, USA, 1859(1), 020024. doi.org/10.1063/1.49 90177.
29. Timoshenko, S. P., & Woinowsky-Krieger, S., (1959). Theory of Plates and Shells (2nd edn.). New York: McGraw-Hill.
30. You, L. H., Zhang, J. J., & You, X. Y., (2005). Elastic analysis of internally pressurized thick-walled spherical pressure vessels of functionally graded materials. International Journal of Pressure Vessels and Piping, 82, 347–354.