Research Article  Open Access
Beata Śmielak, Leszek Klimek, Jacek Świniarski, "The Use of the Finite Elements Method (FEM) to Determine the Optimal Angle of Force Application in Relation to Grooves Notched into a Zirconia Coping with the Aim of Reducing Load on a Connection with Veneering Ceramic", BioMed Research International, vol. 2019, Article ID 7485409, 9 pages, 2019. https://doi.org/10.1155/2019/7485409
The Use of the Finite Elements Method (FEM) to Determine the Optimal Angle of Force Application in Relation to Grooves Notched into a Zirconia Coping with the Aim of Reducing Load on a Connection with Veneering Ceramic
Abstract
Objective of Study. To investigate, using the FEM, the influence of different notching angles on a zirconium dioxide coping with the aim of establishing the optimal connection conditions with veneering ceramic. Materials and Methods. To calculate the stresses in the connection between zirconia coping and veneering ceramic, a model comprising grooves cut perpendicular was adopted. Such a notch profile was used to design the shape and spacing of the grooves on an FEM model simulating a zirconium dioxide coping. For discretization purposes we used twentynode solid BRICK elements featuring intermediate nodes with three degrees of freedom in each node. The model was divided into 117 745 finished elements and 439 131 nodes. The problem was solved with a GLU type contact. The same load F = 1N divided by the number of nodes on the external surface was applied to each node of the outer surface of the base. In subsequent computing variants the F load changed the orientation by angle from 0° to 45° every 15°. Results. The highest level of material strain occurs at angle = 0° =309 MPa and the lowest at angle = 45° =220 MPa. The highest positive stress pressure occurs at angle = 0° =251 MPa, =354 MPa and the lowest at angle = 15°, =171 MPa, =186 MPa. In the case of tangential stresses on the copingveneering ceramic connection, the highest values were noted at angle = 15° =44,4 MPa and the lowest at angle = 45° =32,7 MPa. Conclusions. To reduce the load on the zirconiaveneering ceramic connection, the notches should be made at an angle of = 45°.
1. Introduction
Despite their many advantages, zirconiabased restorations have poor connection parameters with veneering ceramic. More often than is the case with complex metalporcelain structures veneering porcelain is at risk of chipping (1562%), cracking (2550%), delamination (> 10.7%), and major fractures (333%) [1–5].
The majority of authors are of the opinion that the quality of the connection between ceramic and zirconium dioxide depends on a suitably developed coping surface and the creation of microattachments between uneven areas of the base and the liquid ceramic [6, 7]. However, because of its considerable hardness, zirconium dioxide is more resistant to conventional machining than traditional ceramic and requires more aggressive action. The effects of acid etching are measurable on the nanoscale [8, 9]. In turn, abrasive blasting is a highly controversial approach [10, 11]. Fischer et al. report that abrasive blasting is not necessary at all [12]. Many studies reject the hypothesis that it significantly improves adhesion, regardless of the size of the aluminium oxide grains [13, 14]. Other studies point to the negative influence of sandblasting on the structure of zirconium [15].
On one hand, other studies have demonstrated promising results for conditioning with the following lasers Nd:YAG, Er:YAG, CO2, XeCl, KrF, and ArF [16–19]. Various studies have compared the influence of different wavelengths, frequency, laser power, continuous wave mode, and pulse mode on roughness parameters and surface quality in SEM images. No one, on the other hand, has focused on the orientation of notches in relation to the load and its impact on connection strain. According to one study by Fischer et al., connection quality is dependent on the mechanical anchoring, the type and concentration of surface defects, and the level of stress on the ceramic layer [12].
The finite element method (FEM) allows us to establish the ideal conditions for surface development using a laser and assess connection quality after applying appropriate force. The FEM is a modern technique that is also currently widely used in scientific studies in the field of biomechanics. Taking into account the studies mentioned above, the question asked in the present paper is whether notches made with a diode laser should take the form of peripheral lines around the coping as well as lines running perpendicular to those already formed.
The aim of the present study was to find the optimal angle of force application in relation to the grooves notched into a zirconia coping, so as to reduce to a maximum the load on the connection between a coping and veneering ceramic.
2. Material and Methods
Ten rectangular 3YTZP Ceramill Zi sintered zirconia samples (Amann Girrbach AG, Koblach, Austria) with dimensions 10 x 10 x 5 mm underwent laser etching (Figure 1). One randomly chosen sample was selected for further testing. The choice of groove pattern was based on the results of previous studies in which the optimal values of the characteristic dimensions were determined [20]. To calculate the stresses in the connection between zirconia coping and veneering ceramic, we adopted a model comprising grooves cut perpendicular to one another with a width of 0.075 mm, a depth of 0.017 mm, and spacing of 0.115 mm. Such a notch profile was used to design the shape and spacing of the grooves on an FEM model simulating a zirconium dioxide coping.
The mechanical properties of the materials adopted for the FEM calculations are shown in Table 1.

The mechanical properties of materials used for FEM calculation purposes are presented in Table 1.
Table 1 is reproduced from [20] [under the Creative Commons Attribution License/public domain].
Laser scanning microscopy showed that, at the point where the notches intersected, a crater had formed with an average depth of 0.01 mm as a result of the crisscrossing of laser beams during machining (this site can be said to have been machined twice, hence the recess). This crater was also taken into account in the adopted model Figures 2(b) and 3(b) as well as Figure 3(c). This crater is also taken into account in the adopted model. For discretization purposes we used twentynode solid BRICK elements featuring intermediate nodes with three degrees of freedom in each node. The FEM model is presented in Figure 2. The zirconium dioxide coping is marked blue and the ceramic is marked green. The zirconium coping was supported in each direction in accordance with the symmetry conditions.
(a)
(b)
(c)
The model was divided into 117 745 solid elements as well as 439 131 nodes. The FEM mesh was chosen in such a way that a twofold reduction in the size of the mesh would not cause a change in stress of greater than 2% between the zirconia and ceramic layers. This problem was solved by means of a BONDED type contact, i.e., with no possibility of separation, between the zirconia and the ceramic layer. The axial load was applied in the base shown in Figure 4(a) using RBE (RIGID BODY ELEMENT) elements in such a way that the same load F = 1N divided by the number of nodes on the external surface is applied to each node of the outer surface of the base. This type of load discretization corresponds to the even load of the hold along the entire circumference. In further calculation variants the F load changed the orientation by angle from 0° to 45° every 15°. The location of the angle and load are shown in Figure 4(b).
3. Results
The results are presented in the form of colour maps of normal stresses, i.e., stresses perpendicular to the surface of the connection between the zirconia coping and the veneering ceramic, tangential stresses, i.e., the stresses lying on this surface, and reduced according to the Mises hypothesis. The results show the level of strain on the zirconia coping and the ceramic. They should be treated in a qualitative rather than a quantitative sense, as they show which material experiences more strain.
The colour code, ranging from navy blue to red in the computer printout legend, illustrates the increase in stress values. Identical colouring in a given area of the mathematical model indicates approximately the same stress value at a given moment.
3.1. Results of Calculations
The results of the calculations for reduced stresses are compared with each other in the form of a map of reduced stresses according to the Huber (Mises) hypothesis and presented in Figure 5.
(a)
(b)
(c)
(d)
The results of these calculations are presented in Figure 6, so as to facilitate interpretation of the calculation results. It shows the maximum reduced strain as a function of the angle of load rotation.
The connection tightness may be disturbed by the ceramic breaking off from the base. This means that the adhesive bonds have become detached. The possibility of such a situation arising is illustrated by the distributions of pressure between the zirconia coping and the ceramic. The concept of normal pressure means that this value is always oriented perpendicularly to the contact surfaces of both materials, regardless of the orientation of the adjacent surfaces towards each other. In this case, our primary concern will not be the highest positive values, because here the ceramic presses onto the coping, but rather negative values (occurring on the opposite surface). In the latter scenario the ceramic breaks off from the coping and local layers of the connection are at risk of separating, which can lead to delamination.
Figure 7 shows the pressure distribution on the contact surface of the ceramic. A positive value denotes pressure on the connection while a negative value indicates that the connection has been broken. The results of the calculations are presented in Figure 7 in the form of a graph showing both pressure on the connection (a positive stress value) and detachment of the adhesive connection (a negative stress value).
(a)
(b)
(c)
(d)
As Figure 8 shows, the greatest connection strain caused by surface pressures can be observed at load angle = 0°. The lowest pressures were observed at angle = 15°. The tensile stress break forces are lower at angle = 45° than at angle = 0° but greater than at angle = 15°.
In addition to the reduced stresses, corresponding in a straight line to the first signs of internal cracks in the ceramic, as well as the normal stresses responsible for the detachment of adhesive bonds, it is also worth paying attention to the shear stresses occurring in the connection surface between the ceramic and the zirconia coping. This type of stress is particularly dangerous because in addition to thermal deformations it can initiate delamination between the layers, which when combined with normal stresses can quickly result in the ceramic breaking off from the coping. Figure 9 shows the distributions of tangential stresses on the surfaces of the connection between the zirconia coping and the veneering ceramic, whereas Figure 10 shows the distribution of these stresses as a function of the angle of load application in relation to the arrangement of the grooves.
(a)
(b)
(c)
(d)
As can be seen in this case, shear stresses were highest where minimum normal stresses were observed ( = 15°). From Figure 10 it can be seen that, just as with the reduced stresses, the least strained connections resulting from shearing forces were observed at angle = 45°.
4. Discussion
The FEM allows us to study structures with complex shapes and composed of various materials and subjected to any loads, without the need for their physical construction or loading. The tests are carried out on computer models of the examined objects. The accuracy of the calculations depends on the number of elements that the model is divided into. To achieve credible results the model should be divided into a large number of elements, which at the same time increases the degrees of freedom and, as a consequence, significantly extends the computation time [21]. Besides this, the boundary conditions, especially in terms of loading, should be very carefully analysed.
In the presented FEM tests, it was assumed that there was an ideal connection between the different materials and that it remained undamaged despite the increase in loading. In clinical reality, no ideal connection exists. This is due to inaccuracies in laboratorymade connections [22].
The purpose of the calculations was to determine reduced stresses in ceramics, the pressure between a zirconia coping layer and the veneering ceramic, and shear stresses occurring at the boundary between the coping and the veneering layer.
The notches made on the coping surface and the placement of the veneering material in these places were aimed at reducing the load on the adhesive connections between the ceramic and the coping to facilitate mechanical load transfer (mechanical lock). This in turn was designed to ensure the combination of both materials was not purely adhesive in character but rather mechanicaladhesive.
To ensure the complete and correct transfer of the load via mechanical phenomena such as shearing of the coping or the veneering ceramic, the grooves on the zirconia coping should have walls set perpendicular to the contact surface of both materials, which would offer the best solution for the connection. Then, instead of shearing the adhesive layer the lateral forces would be transferred through the volume of material present in the created groove. However, the use of laser technology to make grooves makes it impossible to construct groove walls for the coping and veneering ceramic that would run perpendicular to the contact surface. There are also limitations regarding their size and depth.
If a clean mechanical connection cannot be achieved, the inclined groove walls are loaded with surface pressures resulting from the shape of the connection as well as with shearing forces of the ceramiccoping connection. It is well known that a flat connection only has the strength of an adhesive connection. A mixed adhesivemechanical connection partially eliminates the risk of adhesive connection detachment by transferring loads via microattachments, i.e., “locks” formed by the notches in the coping surface. However, there is a tendency for stress to accumulate in the corners of grooves as a result of transferred mechanical loads. The accumulation or otherwise stress concentrations or material strain only affects small areas. In the case of an actual nonidealized surface, we can assume that these concentrations will be blurred by mechanical microlocks (irregularities on the processed surface).
It is important to stress that the calculation results presented here are qualitative and not quantitative in character. They show which elements are subject to most strain. For computation purposes a loading force of 1N was adopted. The application of a different force will of course change the stress values, but it will not alter the relationship between them. Not all aspects can be modelled according to the FEM calculation, because there is no table with critical stress values for the adhesion and for the shearing of the bond.
The calculations show that when it comes to reduced stresses, the most advantageous situation is one where the loading force is directed at an angle of 45° in relation to the groove system (Figure 5(d)) and this arrangement is the most effective in terms of the material load. However, most often damage to the copingveneering ceramic connection occurs within the interface between these two materials and therefore an analysis of contact pressures and shear stresses is more important. As far as pressure is concerned, the most favourable conditions arise when the load is applied at an angle of 15° to the groove system (Figures 7(b) and 8). However, an analysis of the tangential stresses responsible for the delamination of the connection shows that these forces are highest at an angle of 15° and lowest at an angle of 45°. Summing up, it appears that the most favourable load system is one where the loading force is applied at an angle of close to 45° in relation to the groove system.
5. Conclusions
(1)A change in the orientation of the notches in relation to the load reduces connection strain and improves connection stability.(2)It has been demonstrated that the notches should be made at an angle = 45° so as to reduce stress concentrations and thus minimise the stress occurring in the corners of the connection.(3)It is recommended that cuts be made in two mutually perpendicular directions so as to make orientation independent of the load.
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare no conflicts of interest.
Authors’ Contributions
Beata Śmielak and Leszek Klimek contributed to conceptualization; Jacek Świniarski contributed to methodology; Jacek Świniarski contributed to software; Beata Śmielak, Leszek Klimek, and Jacek Świniarski contributed to validation; Jacek Świniarski contributed to formal analysis; Beata Śmielak contributed to investigation; Beata Śmielak contributed to resources; Leszek Klimek contributed to data curation; Beata Śmielak contributed to writingoriginal draft preparation; Beata Śmielak contributed to writingreview and editing; Beata Śmielak contributed to visualization; Leszek Klimek contributed to supervision; Beata Śmielak contributed to project administration; Beata Śmielak contributed to funding acquisition.
References
 F. Komine, M. B. Blatz, and H. Matsumura, “Current status of zirconiabased fixed restorations,” Journal of oral science, vol. 52, no. 4, pp. 531–539, 2010. View at: Publisher Site  Google Scholar
 M. Schmitter, D. Mueller, and S. Rues, “Chipping behaviour of allceramic crowns with zirconia framework and CAD/CAM manufactured veneer,” Journal of Dentistry, vol. 40, no. 2, pp. 154–162, 2012. View at: Publisher Site  Google Scholar
 P. Vigolo and S. Mutinelli, “Evaluation of zirconiumoxidebased ceramic singleunit posterior fixed dental prostheses (FDPs) generated with two CAD/CAM systems compared to porcelainfusedtometal singleunit posterior FDPs: a 5year clinical prospective study,” Journal of Prosthodontics, vol. 21, no. 4, pp. 265–269, 2012. View at: Publisher Site  Google Scholar
 P. Benetti, F. Pelogia, L. F. Valandro, M. A. Bottino, and A. D. Bona, “The effect of porcelain thickness and surface liner application on the fracture behavior of a ceramic system,” Dental Materials, vol. 27, no. 9, pp. 948–953, 2011. View at: Publisher Site  Google Scholar
 I. Sailer, A. Feher, F. Filser, L. J. Gauckler, H. Luthy, and C. H. Hammerle, “Fiveyear clinical results of zirconia frameworks for posterior fixed partial dentures,” The International Journal of Prosthodontics, vol. 20, pp. 383–386, 2007. View at: Google Scholar
 A. Eichler, “Tetragonal Ydoped zirconia: structure and ion conductivity,” Physical Review B, vol. 64, pp. 174103–174108, 2001. View at: Google Scholar
 M. Inokoshi, K. Yoshihara, N. Nagaoka et al., “Structural and chemical analysis of the zirconia veneering ceramic interface,” Journal of Dental Research, vol. 95, no. 1, pp. 102–109, 2015. View at: Publisher Site  Google Scholar
 B. Smielak and L. Klimek, “Effect of hydrofluoric acid concentration and etching duration on select surface roughness parameters for zirconia,” The Journal of Prosthetic Dentistry, vol. 113, no. 6, pp. 596–602, 2015. View at: Publisher Site  Google Scholar
 A. Casucci, E. Osorio, R. Osorio et al., “Influence of different surface treatments on surface zirconia frameworks,” Journal of Dentistry, vol. 37, no. 11, pp. 891–897, 2009. View at: Publisher Site  Google Scholar
 J. Teng, H. Wang, Y. Liao, and X. Liang, “Evaluation of a conditioning method to improve coreveneer bond strength of zirconia restorations,” Journal of Prosthetic Dentistry, vol. 107, pp. 380–387, 2009. View at: Google Scholar
 H. Kim, H. Lim, Y. Park, and M. Vang, “Effect of zirconia surface treatments on the shear bond strength of veneering ceramic,” The Journal of Prosthetic Dentistry, vol. 105, no. 5, pp. 315–322, 2011. View at: Publisher Site  Google Scholar
 J. Fischer, P. Grohmann, and B. Stawarczyk, “Effect of zirconia surface treatments on the shear strength of zirconia/veneering ceramic composites,” Dental Materials, vol. 27, no. 3, pp. 448–454, 2008. View at: Publisher Site  Google Scholar
 V. Lughi and V. Sergo, “Low temperature degradation aging of zirconia: a critical review of the relevant aspects in dentistry,” Dental Materials, vol. 26, no. 8, pp. 807–820, 2010. View at: Publisher Site  Google Scholar
 K. Tada, T. Sato, and M. Yoshinari, “Influence of surface treatment on bond strength of veneering ceramics fused to zirconia,” Dental Materials Journal, vol. 31, no. 2, pp. 287–296, 2012. View at: Publisher Site  Google Scholar
 B. Stawarczyk, M. Özcan, M. Roos et al., “Fracture load and failure analysis of zirconia single crowns veneered withpressed and layered ceramics after chewing simulation,” Dental Materials Journal, vol. 30, no. 4, pp. 554–562, 2011. View at: Publisher Site  Google Scholar
 S. Armani, M. H. Tabatabase, S. F. Namdar, and N. Chiniforush, “Effect od different lasers and particle abrasion on surface characteristics of zirconia ceramics,” Journal of Dentistry, vol. 11, pp. 233–241, 2014. View at: Google Scholar
 B. Ersu, B. Yuzugullu, A. Ruya Yazici, and S. Canay, “Surface roughness and bond strength of glass infiltrated aluminaceramics prepared using various surface treatments,” Journal of Dentistry, vol. 37, no. 11, pp. 848–856, 2009. View at: Publisher Site  Google Scholar
 L. Hao and J. Lawrence, “Effects of Nd:YAG laser treatment on the wettability characteristics of a zirconiabased bioceramic,” Optics and Lasers in Engineering, vol. 44, no. 8, pp. 803–814, 2006. View at: Publisher Site  Google Scholar
 G. Nicolas, M. Autric, W. Marine, and G. A. Shafeev, “Laser induced surface modifications on ZrO2 ceramics,” in Proceedings of the 1996 EMRS Spring Conference, pp. 289–292, June 1996. View at: Google Scholar
 B. Śmielak, L. Klimek, and J. Świniarski, “The use of the FEM to identify the optimal groove dimensions ensuring the least stressed connection between a zirconia coping and veneering ceramic,” Materials, vol. 11, no. 12, article no 2360, 2018. View at: Publisher Site  Google Scholar
 J. Orr and C. Mitchell, “The evaluation of postretained crowns using a custom designed fatigue machine,” Journal of Oral Rehabilitation, vol. 23, no. 3, pp. 179–185, 1996. View at: Publisher Site  Google Scholar
 C. Mitchell, J. Orr, and J. Kennedy, “A semiempirical model for prediction of how postretained crowns will fail under compressive loading,” Journal of Dental Research, vol. 71, no. 9, pp. 1613–1618, 2016. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2019 Beata Śmielak et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.