The Effect of Preparation Design on Stress Distribution in Lithium Disilicate Laminate Veneer using 3D Finite Element Analysis.
RESEARCH ARTICLE

The Effect of Preparation Design on Stress Distribution in Lithium Disilicate Laminate Veneer using 3D Finite Element Analysis.

The Open Dentistry Journal 12 Oct 2023 RESEARCH ARTICLE DOI: 10.2174/0118742106243827230921095501

Abstract

Background:

A layer of ceramic bonded to the facial surfaces of prepared teeth to create a ceramic esthetic veneer, also known as a laminate veneer.

Objective:

This study aimed to examine the effect of preparation design on the stress distribution in lithium disilicate veneer at different angulations of masticatory loading.

Methods:

The finite element approach assessed the structure stress distribution. So, Exocad dental CAD was used to obtain the models, which consist of an upper central incisor in three dimensions (3D), and prepare them with the same software for an indirect veneer. The models were imported into software (SolidWorks; Dassault Systems) to construct solid models. The external load used in this study was 100 N at 60 degrees (intercuspal) and 125 degrees (protrusive) in the palatal surface, two mm below the incisal edge. The models were restored using e-max laminate veneers with different tooth preparation designs (palatal chamfer and window preparation designs). Analysis and calculations were performed on the stress distributions and maximum stress values of the tooth structures, cement layer, and laminate veneer.

Results:

With a window preparation design, the maximum stress values in the tooth and laminate veneer were higher. In the palatal chamfer design, stresses were dispersed more consistently through the cement layer and laminate veneer. Higher stress was generated by protrusive movement than by intercuspal movement.

Conclusion:

Compared to window preparation, the palatal chamfer design for lithium disilicate laminate veneer showed a more favourable stress distribution. Regardless of the preparation design, maximum stress values were observed in the cervical and incisor areas of the tooth structure.

Keyword: laminate veneer, Preparation designs, Palatal chamfer design, Window design, Stress distribution, Three-dimensional finite element analysis.

1. INTRODUCTION

A ceramic layer is bound to the facial surfaces of prepared teeth to create a ceramic esthetic veneer, also known as a laminate veneer. Ceramic laminate veneers are produced specifically for each patient and are created in a dental lab. Most ceramic laminate veneers are made by machining or heat-pressing lithium disilicate or leucite-reinforced ceramic [1].

The principle of conservative dentistry with less invasive procedures involved ceramic laminate veneer preparations. Furthermore, this method has been effectively applied to correct flaws, such as diastema, minor fractures, and stains and improve patients' aesthetic conditions [2, 3].

Various tooth preparation geometries for ceramic laminate veneers are described in the literature. Two of them are used in this study because they offer more benefits. These are (1) the “window preparation,” which involves minimal preparation of the labial surface and doesn't involve the incisal edge; and (2) the “palatal chamfer preparation,” which involves labial surface preparation with incisal edge reduction. The durability of ceramic veneers is not primarily determined by the materials used. Previous studies have demonstrated that the preparation design impacts the ultimate result [4, 5]. Even though a rounded shoulder or chamfer has been proven to be the best marginal design for obtaining an obvious finish line and suitable structural resistance, more approval must be given regarding the potential lingual extension of the veneer preparations [6, 7].

Since practically all of the information about the clinical attitude of these restorations was based on personal records, the effect of various tooth preparation geometries on the fracture strength and stress distribution of laminate veneers is still controversial [8, 9].

The several events in the examined models can be better understood by finite element analysis (FEA), a non-destructive test connected to in vitro tests. The finite element method enables a more exact system behaviour to check by better reviewing the stress distribution within the materials or at an interface. Every element's mechanical property can be written as a function of node translation. When such nodes are loaded in specific ways, the model operates similarly to the structure it represents [10-12].

So, the study's objective in vitro is to examine the effects of tooth preparation designs on the stress distributions in lithium disilicate veneer at different angulations of masticatory force loading (60 and 125 degrees). The null hypothesis is that there is no difference in stress distribution between the window preparation design and the palatal chamfer design.

2. MATERIALS AND METHODS

The model consisting of maxillary central incisors was obtained from Exocad dental CAD (Exocad GnbH, Germany) and prepared by the same software (Exocad dental CAD) with the protocol for clinical preparation:

• 1 mm interproximal reduction.

• The cervical border is 1 mm coronal to the cement-enamel junction (CEJ).

• 0.7 mm reduction of the buccal aspect for palatal chamfer preparation (extension two mm toward the palatal aspect to cover the incisal edge with the chamfer finish line).

• 0.7 mm reduction of the buccal aspect for window preparation.

Finally, two digitally prepared tooth models were obtained (Fig. 1).

The ceramic veneer (IPS e-max Press, Ivoclar Vivdent, Schaan, Liechtenstein) with a thickness of 0.5mm and a cement layer of 0.2mm (RelyX™, 3MESPE, USA), according to manufacture instructions, was designed on the digitally prepared tooth using Exocad dental CAD (Exocad GnbH, Germany) [13]. The models consist of four main volumes: the enamel, the dentin, the cement, and the veneer structures. The models were imported into software (SolidWorks; Dassault Systems) to construct solid models. Four models have been obtained, as reported in Table 1. The models didn't include roots, PDL, or bone. The mesh is formed by a mix of tetrahedral, linear triangular, and linear quadrilateral elements (243309 elements and 314365 nodes), and the solids are considered isotropic, linearly elastic, and homogenous. The mechanical properties used in this study are elastic modulus and Poisson's ratio for laminate veneer, cement layer, enamel, and dentin.

Fig. (1A). Palatal chamfer.
Fig. (1B). Window preparation.
Table 1.
3D-FEM model descriptions.
Model No. Preparation Technique Load Applications
Mod. 1 Window 60°
Mod. 2 Palatal chamfer 60°
Mod. 3 Window 125°
Mod. 4 Palatal chamfer 125°
Table 2.
Material's mechanical properties.
Material Elastic Modulus (GPa) Poisson Ratio
Enamel [14] 84.1 0.3
Dentin [14] 18.6 0.31
Resin cement 5 0.35 [14]
Lithium disilicate-based ceramic (IPS.e-max Press) 95 0.3 [15]
Note: *a: Values provided by the manufacturer.
**GPa: Gigapascal.

Regarding the laminate veneer and cement layer, the mechanical properties (elastic modulus and Poisson's ratio) of a disilicate lithium ceramic (IPS e-max press, Ivoclar Vivadent, Schaan, Liechtenstein) and a methacrylate resin-based luting agent (RelyX™ Veneer Cement, 3M ESPE, USA) were taken into account. The mechanical properties of the isotropic materials are described in Table 2. The external load used in this study was 100 N at 60 degrees (intercuspal) and 125 degrees (protrusive) on the palatal surface, two mm below the incisal edge. The load was applied separately to selected fixed nodes for all models.

2.1. Experimental Test

The finite element solution of each model was carried out using a three-dimensional planar strain model and the Abaqus/19 package. A structural linear static analysis assessed the stress distribution in important areas. Maximum principal stress was used to analyze the stress distribution in tooth structure, cement, and veneer.

3. RESULTS

There are two fundamental finite element models. Table 3 reports the maximum principal stress that was recorded. Maximum stress areas are reported in Table 4. The stress distributions in laminate veneer under intercuspal load for Mod. 1 and Mod. 2 are compared (Fig. 2a). The stress in the veneer was uniformly distributed in both models. The maximum stress value in Mod. 1 was 200 MPa in the cervical area, while in Mod. 2, the maximum value reached was 148 MPa in the incisal region. In Fig. (2b), the stress distributions in laminate veneer under protrusive load are marked (Mod. 3 and Mod. 4). The maximum value that occurred in Mod. 3 was 236 MPa. In Mod. 3, the stress concentrated in the cervical area and progressively decreased in the buccal area. In Mod. 4, maximum concentration stress (196 MPa) occurred in the incisal region and gradually decreased in the cervical region.

Table 3.
Maximum stress value.
Model No. Laminate Veneer
(MPa)
Cement Layer
(MPa)
Tooth Structure
(MPa)
Mod. 1 200 25 115
Mod. 2 148 11.5 82
Mod. 3 236 32.8 124
Mod. 4 196 15 94
Note: *MPa: Megapascal.
Table 4.
Maximum stress areas.
Models No. Max. Stress Area
Veneer Cement Tooth Structure
Mod. 1 Cervical area Incisal area Cervical area
Mod. 2 Incisal area Buccal surface Cervical, incisal area
Mod. 3 Cervical area Incisal area Cervical area
Mod. 4 Incisal area Buccal surface Cervical, incisal area
Fig. (2a). Mod.1: stress distribution in laminate veneer under tearing load with window preparation.
Mod.2: stress distribution in laminate veneer under tearing load with palatal chamfer preparation.
Fig. (2b). Mod.3: stress distribution in laminate veneer under protrusive load with window preparation.
Mod.4: stress distribution in laminate veneer under protrusive load with palatal chamfer preparation.

The stress distribution in the cement layer under intercuspal movement in Mod. 1 and Mod. 2 was compared (Fig. 3a). The maximum stress value that occurred in Mod. 1 was 25 MPa. In Mod2, the maximum stress value recorded was 11.5 MPa. In Mod. 2, the stress was focused on the buccal and incisal areas. In Fig. (3b), the stress distribution in the cement layer under the protrusive load was marked (Mod. 3 and Mod. 4). In Mod. 3, the stress reached a maximum value of 32.8 MPa, and the stress was concentrated in the cervical and incisal areas. In Mod. 4, the maximum stress value was 15 MPa. The high stress was concentrated in the buccal and incisal areas in Mod. 4.

Fig. (3a). Mod.1: stress distribution in cement layer under tearing load with window preparation.
Mod.2: stress distribution in cement layer under tearing load with palatal chamfer preparation.
Fig. (3b). Mod.3: stress distribution in cement layer under protrusive load with window preparation.
Mod.4: stress distribution in cement layer under protrusive load with palatal chamfer preparation.
Fig. (4a). Mod.1: stress distribution in tooth structure(lingual surface) under tearing load with window preparation.
Mod.2: stress distribution in tooth structure(lingual surface) under tearing load with palatal chamfer preparation.
Fig. (4b). Mod.3: stress distribution in tooth structure(lingual surface) under protrusive load with window preparation.
Mod.4: stress distribution in tooth structure(lingual surface) under protrusive load with palatal chamfer preparation.

In Fig. (4a), the stress distributions in dental tissue under intercuspal load were marked (Mod. 1 and Mod. 2). In Mod. 1, the stress in the incisal region was decreased. The maximum stress value (115 MPa) was recorded in the cervical area of the palatal aspect. In Mod. 2, the maximum stress value was 82 MPa. Stress was concentrated in the cervical and incisal areas of the palatal aspect in Mod. 2. In Fig. (4b), reports the stress distributions under protrusive movement in tooth structure (Mod. 3 and Mod. 4). In Mod. 3, the higher stress (124 MPa) is recorded in the cervical area of the palatal aspect and reduced in the incisal area of dental tissue. In Mod. 4, the maximum stress value (94 MPa) was recorded in the cervical and incisal areas of the palatal aspect.

4. DISCUSSION

Finite element analysis was used in this study to monitor and analyze stress distribution. It is often regarded as the most accurate method because it allows researchers to quantitatively assess complicated biological systems for which other approaches are impractical or unsolvable [16-18].

A mathematical representation of a real object or phenomenon is a numerical model, such as one that uses FEM technology. The fact that a mathematical model is virtual—that is, it resides in a computer and is entirely programm- able—gives it significant advantages over in-vivo testing. The investigator can easily change the test parameters, model parameters, and geometries. Additionally, the investigator can copy any preferred result and reiterate the test (simulation) anytime. All of these things are not achievable during in vivo studies. Consequently, a well-validated mathematical model gives the investigator an extremely potent tool for analysis [19, 20].

The likelihood of fluctuation in the results due to sampling error is eliminated since numerical values obtained from stress analysis are mathematical calculations without variance and because variables can be changed with computer precision. Any number of times the same FEM analysis was repeated, it would always produce the same results. As a result, it is evident that outcomes are never the result of chance events but rather are the result of variable control. Because of this, a conclusion from a FEM study is typically reached without using formal inferential statistical analysis [21, 22].

The result of this study reveals that the window preparation technique develops higher stress than the palatal chamfer preparation in both 125° (protrusive) and 60° (intercuspal) load conditions (Figs. 2-4), and this is agreed with by Zarone et al. in 2005 [23]. Such evidence appears to be due to the load direction. When a tooth is subjected to therapeutic veneer treatment, the first substrate located by the load direction is enamel, which performs similarly regardless of the ceramic and cement agent, reflecting similar results for dentin between all groups. This result can be clarified because the thickness of the remaining dentin is the same for each group and close enough to the load to prevent variations in behaviour groups. The maximum principal stress result demonstrated that, regardless of the restorative material or occlusal condition, the palatal chamfer design reduces stress in both the tooth structure and the laminate veneer, as seen in other investigations [23, 24]. (Table 3).

Regarding loading circumstances, performing window preparation makes the intercuspal and protrusive movements more stressful (Figs. 2-4). While in palatal chamfer preparation, the opposite situation was noticed. In this study, most of the tooth structure can evenly transport the load from the tooth to the cement layer and then to the veneer because the tooth structure is subject to compressive force.

The e-max veneer is a barrier during functional motions by absorbing most stress and preserving underlying dental tissue according to the stress variance between the veneer and tooth structure (Table 3). All models saw larger stress concentrations at the veneer structure level. This phenomenon can be explained by the therapeutic system's overall mechanical behaviour. The e-max veneer has a higher elastic modulus than the cement and tooth structure. Since the prepared tooth's weakest link is the incisal and cervical regions, high stress is concentrated in these areas. The cement layer was a crucial area for biomechanical behaviour [25]. The outcome demonstrated that the material with a lower modulus of elasticity (less stiffness) produces less concentration stress, as reported by other evidence [26]. The cement layer, however, served as a stress absorber (Fig. 3).

Considering the three layers of the examined restorative system (tooth structure, cement layer, and laminate veneer), when utilizing window preparation, the stress was consistently localized at the cervical level (Figs. 2-4). When employing palatal chamfer preparation, the stress was dispersed from the cervical area of the palatal aspect to the area of load application in the middle region of the incisal area. The stress within the veneer and tooth structure decreases as the preparation depth increases. Prior 3D-FEA research claimed that stress was localized in the incisal area [27, 28]. Nevertheless, this study also found the stress concentrated in the cervical area of the palatal surface, which Tsoukindas A. et al. had previously reported on in 2020 [29].

The volume of restoration would be increased by eliminating contact points proximally, typically done in patients with significant discolouration or deformity. In the current investigation, the window preparation and palatal chamfer designs, intended to create a good finish line on the proximal surface, failed to preserve contact points proximally. Castelnuovo et al. showed that the strength of the palatal chamfer design is lower than that of the butt joint design because of the thickness of the ceramic [8]. Based on the current findings, the palatal chamfer preparation design may be a better alternative for laminate veneer since it has a lower maximum main stress and a more consistent stress distribution in the cement layer. This study shows a significant difference in stress distribution between the window and palatal chamfer designs, so the null hypothesis was rejected.

The FEA, as a computational technique, did not take into account all aspects of the oral environment, such as the anisotropic nature of dentin and the nonlinear viscoelastic properties of the periodontal ligament, which both required a significant amount of experimental data. Other limitations include that the mean masticatory load in the anterior region could vary from 93 N to 206 N depending on age, race, gender, tooth loss, TMJ condition and pain, periodontal support of teeth, weight, and height [30-32].

CONCLUSION

Concerning different loading angles, the window and palatal chamfer designs varied mechanically within the confines of the current investigation. Compared to window preparation, the palatal chamfer design for lithium disilicate laminate veneer showed a more favourable stress distribution. Regardless of the preparation design, maximum stress values were observed in the cervical and incisor areas of the tooth structure.

LIST OF ABBREVIATIONS

% = Percentage
2D = Two Dimensional
3D = Three Dimensional
CAD = Computer Aid Design
e.g = For example
FEA = Finite Element Analysis
FEM = Finite Element Method
Fig = Figure
GPa = Gigapascal
Max = Maximum
MM = Millimeter
MPa = Megapascal
Mod = Model
N = Newton
USA = United State America
Tab = Table

ETHICS APPROVAL AND CONSENT TO PARTICIPATE

Not applicable.

HUMAN AND ANIMAL RIGHTS

No human and animals were used for the studies that are the basis of this research.

CONSENT FOR PUBLICATION

Not applicable.

AVAILABILITY OF DATA AND MATERIAL

All the data and supporting information is provided within the article.

CONFLICT OF INTEREST

The authors declare no conflict of interest, financial or otherwise.

FUNDING

None.

ACKNOWLEDGMENTS

Declared none.

REFERENCES

1
Sakaguchi RL. Craig’s restorative dental materials-e-book. In: Ferracane LJ, Powers JM, Eds. Elsevier Health Sciences. 4th ed. 4th ed.Portland, OR 2011; p. 211.
2
Fradeani M, Redemagni M, Corrado M. Porcelain laminate veneers: 6- to 12-year clinical evaluation-a retrospective study. Int J Periodontics Restorative Dent 2005; 25(1): 9-17.
3
Magne P. Efficient 3D finite element analysis of dental restorative procedures using micro-CT data. Dent Mater 2007; 23(5): 539-48.
4
Adebayo OA, Burrow MF, Tyas MJ. Bond strength test: Role of operator skill. Aust Dent J 2008; 53(2): 145-50.
5
Shetty A, Kaiwar A, Shubhashini N, et al. Survival rates of porcelain laminate restoration based on different incisal preparation designs: An analysis. J Conserv Dent 2011; 14(1): 10-5.
6
Costa VL, Tribst JP, Uemura ES, de Morais DC, Borges AL. Influence of thickness and incisal extension of indirect veneers on the biomechanical behaviour of maxillary canine teeth. JRDE 2018; 43(4)
7
Oladapo BI, Abolfazl Zahedi S, Vahidnia F, Ikumapayi OM, Farooq MU. Three-dimensional finite element analysis of a porcelain crowned tooth. Beni Suef Univ J Basic Appl Sci 2018; 7(4): 461-4.
8
Castelnuovo J, Tjan AHL, Phillips K, Nicholls JI, Kois JC. Fracture load and mode of failure of ceramic veneers with different preparations. J Prosthet Dent 2000; 83(2): 171-80.
9
Hahn P, Gustav M, Hellwig E. An in vitro assessment of the strength of porcelain veneers dependent on tooth preparation. J Oral Rehabil 2000; 27(12): 1024-9.
10
Dal Piva AMO, Tribst JPM, Borges ALS, Souza ROA, Bottino MA. CAD-FEA modeling and analysis of different full crown monolithic restorations. Dent Mater 2018; 34(9): 1342-50.
11
Dartora NR, de Conto Ferreira MB, Moris ICM, et al. Effect of intracoronal depth of teeth restored with endo-crowns on fracture resistance: in vitro and 3-dimensional finite element analysis. J Endod 2018; 44(7): 1179-85.
12
Bragança GF, Mazão JD, Versluis A, Soares CJ. Effect of luting materials, presence of tooth preparation, and functional loading on stress distribution on ceramic laminate veneers: A finite element analysis. J Prosthet Dent 2021; 125(5): 778-87.
13
Bahrami M, Shakeri B, Memarian M. Investigating the stress distribution applied to edentulous ridge from polyamide and cobalt-chrome removable-partial-dentures using three-dimensional finite-element-analysis. J Dentomaxillofac Sci 2020; 9(4): 11-9.
14
Dejak B, Młotkowski A, Langot C. Three-dimensional finite element analysis of molars with thin-walled prosthetic crowns made of various materials. Dent Mater 2012; 28(4): 433-41.
15
Coldea A, Fischer J, Swain MV, Thiel N. Damage tolerance of indirect restorative materials (including PICN) after simulated bur adjustments. Dent Mater 2015; 31(6): 684-94.
16
Cakici EB, Guler MS, Guler C, Cakici F, Sen S. Finite element analysis of thermal stress distribution in different restorative materials used in class V cavities. Niger J Clin Pract 2016; 19(1): 30-4.
17
Silva NR, Aguiar GCR, Rodrigues MP, et al. Effect of resin cement porosity on retention of glass-fibre posts to root dentin: An experimental and finite element analysis. Braz Dent J 2015; 26(6): 630-6.
18
Shukri BM, Al-Hashimi WN. Finite element stress analysis of endodontically treated teeth restored by prefabricated posts. JBCD 2005; 17(3)
19
Dawood SN, Al-Zahawi AR, Sabri LA. Mechanical and thermal stress behavior of a conservative proposed veneer preparation design for restoring misaligned anterior teeth: A 3D finite element analysis. Appl Sci 2020; 10(17): 5814.
20
Menicucci G, Mossolov A, Mozzati M, Lorenzetti M, Preti G. Tooth-implant connection: Some biomechanical aspects based on finite element analyses. Clin Oral Implants Res 2002; 13(3): 334-41.
21
Akça K, Çehreli MC, İplikçioğlu H. Evaluation of the mechanical characteristics of the implant-abutment complex of a reduced-diameter morse-taper implant. Clin Oral Implants Res 2003; 14(4): 444-54.
22
Penteado MM, Mendes Tribst JP, Dal Piva AMO, Archangelo KC, Bottino MA, Souto Borges AL. Influence of different restorative material and cement on the stress distribution of ceramic veneer in upper central incisor. IJDR 2020; 31(2): 236-40.
23
Zarone F, Apicella D, Sorrentino R, Ferro V, Aversa R, Apicella A. Influence of tooth preparation design on the stress distribution in maxillary central incisors restored by means of alumina porcelain veneers: A 3D-finite element analysis. Dent Mater 2005; 21(12): 1178-88.
24
Li Z, Yang Z, Zuo L, Meng Y. A three-dimensional finite element study on anterior laminate veneers with different incisal preparations. J Prosthet Dent 2014; 112(2): 325-33.
25
Peumans M, Voet M, De Munck J, Van Landuyt K, Van Ende A, Van Meerbeek B. Four-year clinical evaluation of a self-adhesive luting agent for ceramic inlays. Clin Oral Investig 2013; 17(3): 739-50.
26
de Andrade GS, Tribst JPM, Dal Piva AMO, et al. A study on stress distribution to cement layer and root dentin for post and cores made of CAD/CAM materials with different elasticity modulus in the absence of ferrule. J Clin Exp Dent 2019; 11(1): 0.
27
Sorrentino R, Apicella D, Riccio C, et al. Nonlinear visco‐elastic finite element analysis of different porcelain veneers configuration. J Biomed Mater Res B Appl Biomater 2009; 91(2): 727.: 736.
28
Walls AWG, Steele JG, Wassell RW. Crowns and other extra-coronal restorations: Porcelain laminate veneers. Br Dent J 2002; 193(2): 73-82, 79-82.
29
Tsouknidas A, Karaoglani E, Michailidis N, Kugiumtzis D, Pissiotis A, Michalakis K. Influence of preparation depth and design on stress distribution in maxillary central incisors restored with ceramic veneers: A 3D finite element analysis. J Prosthodont 2020; 29(2): 151-60.
30
Kaya BM, Ergun G. The effect of post length and core material on root fracture with respect to different post materials. Acta Odontol Scand 2013; 71(5): 1063-70.
31
Koc D, Dogan A, Bek B. Bite force and influential factors on bite force measurements: A literature review. Eur J Dent 2010; 4(2): 223-32.
32
Quiudini PR Jr, Pozza DH, Pinto AS, de Arruda MF, Guimarães AS. Differences in bite force between dolichofacial and brachyfacial individuals: Side of mastication, gender, weight and height. J Prosthodont Res 2017; 61(3): 283-9.