Modelling of fracture behaviour in biomaterials☆

Residual dentin thickness and biomechanical performance of post-and-core-restored mandibular premolars: A finite element analysis study

STATEMENT OF PROBLEM

Endodontically treated teeth often require post-and-core restorations for structural support because of extensive hard tissue loss. Assessment of the effect of the residual dentin thickness on the biomechanical performance of these restorations is lacking.

PURPOSE

The purpose of this study was to evaluate the residual dentin thickness in mandibular premolars after post-and-core restorations using cone beam computed tomography (CBCT) and to analyze the stress distribution with finite element analysis (FEA).

MATERIAL AND METHODS

The CBCT data from 236 mandibular premolars having undergone post-and-core restorations were examined. An imaging software program (NNT; NewTom) was used to measure the buccolingual and mesiodistal root diameters in cross-sections 5 to 11 mm from the radiologic apex. The CBCT derived measurements were subsequently integrated into an FEA model. A 3-dimensional (3D) mandibular premolar model reflecting the residual dentin thickness was created with a computer-aided design software program (Hypermesh 9.0; Altair Engineering). A static force of 100 N was applied directly to the buccal cusp tip at 45, 60, 75, and 90 degrees to the long axis of the tooth, and the stress distribution of dentin was analyzed by using an FEA software program (ANSYS APDL 18.0; Ansys Inc).

RESULTS

CBCT analysis showed that the buccolingual root diameter was wider than the mesiodistal diameter and that the residual dentin thickness of the buccal aspect was approximately 0.3 mm thinner compared with the lingual aspect along the root. The proportions of residual dentin thickness values in the buccolingual direction of the mandibular premolar teeth no less than 1 mm exceeded 96.2% at 5 to 11 mm from the apex. The proportions of residual dentin thickness values in the mesiodistal direction of the mandibular premolar teeth >1 mm were 92.1% and 88.2% at 11 mm from the apex after post space preparation and decreased further to 70.8% and 58.9% at 5 mm from the apex. The von Mises stresses of the mandibular premolar model with residual dentin thickness were mainly localized to the cervical area (region C, cervix) and the post apex (region A, apex) in the buccolingual direction. Tensile and compressive stress were concentrated on regions C and A on the buccal and lingual sides, respectively. The actual residual dentin thickness model demonstrated higher maximum tensile stress compared with the 1-mm residual dentin thickness model under various loading conditions.

CONCLUSIONS

During the process of post space preparation in mandibular premolars, sufficient dentin thickness should be retained in the apical region. The tensile stresses of mandibular premolars after clinical post-and-core restoration were mainly concentrated in the cervical area and the post apex, and the maximum tensile stress value was higher than the ideal 1-mm residual dentin thickness model.

Study on Optimization of the Structural Mechanical Properties of Personalized Porous Implant Prosthesis

Porous implant prostheses can effectively reduce the stress shielding effect. Still, the single elastic modulus prosthesis cannot adapt to the individual skeletal variability, so it is necessary to optimize the structural parameters of the prosthesis to overcome the individual variability. In this regard, this study analyzes the law of structural parameters and mechanical properties after selecting the type of porous structure (diamond structure). It proposes the optimization method of the structural parameters on this basis. First, the functional relationship equations between the unit mass of the porous implant prosthesis, the elastic modulus of the porous implant prosthesis, and the structural parameters were established respectively. Second, the support rod length and radius of the porous implant prosthesis are optimized by a genetic algorithm to form the optimization design method of the porous implant prosthesis. Finally, the feasibility and effectiveness of the optimized design of the porosity implanted prosthesis were verified by animal experiments, and the optimized implanted prosthesis with optimized structural parameters increased bone growth by 20-30% compared to the control group in the animal body. The proposed method provides a theoretical basis and technical support for the rehabilitation of patients and the production of prostheses by physicians.

On the fracture behavior of molar teeth with MOD cavity preparation

Mesial-occlusal-distal (MOD) cavity preparations are commonly used to restore damaged teeth. While numerous in vitro cavity designs have been devised and tested, no analytical frameworks for assessing their resistance to fracture seem to exist. This concern is addressed here by resorting to a 2D slice specimen cut from restored molar teeth with rectangular-base MOD cavity. The evolution of damage due to axial cylindrical indentation is followed in situ. The failure begins with a rapid debonding along the tooth/filler interface and continues with unstable cracking from the cavity corner. The debonding load q_d is fairly fixed while the failure load q_f is insensitive to the presence of filler, increasing with cavity wall thickness h and reducing with cavity depth D. The growth of the corner crack is studied using a 2D fracture analysis in conjunction with the FEM technique. The ratio h = h/D emerges as a viable system parameter. A simple expression for q_f given in terms of h and dentin toughness K_C is developed that predicts well the test data. In vitro studies on full-fledged molar teeth with MOD cavity preparation show that the fracture resistance of filled cavities often exceeds by a large margin that of unfilled ones. Indications are that this may reflect load sharing with the filler. Thus, the fracture resistance of the unfilled cavity provides a lower bound to a compromised MOD filling after long-term aging in the mouth. This bound is well predicted by the slice model. Finally, it is recommended that MOD cavities be prepared, if applicable, such that h > D regardless of the tooth size.

On interaction between fatigue of reconstruction plate and time-dependent bone remodeling

BACKGROUND AND OBJECTIVE

The fibula free flap (FFF) has been extensively used to repair large segmental bone defects in the maxillofacial region. The reconstruction plate plays a key role in maintaining stability and load-sharing while the fibula unites with adjacent bone in the course of healing and remodeling. However, not all fibula flaps would fully unite, and fatigue of prosthetic devices has been recognized as one major concern for long-term load-bearing applications. This study aims to develop a numerical approach for predicting the fatigue life of the reconstruction plate by taking into account the effect of ongoing bone remodeling.

METHODS

The patient-specific mandible reconstruction with a prosthetic system is studied in this work. The 3D finite element model with heterogeneous material properties obtained from clinical computerized tomography (CT) data is developed for bone, and eXtended Finite Element Method (XFEM) is adopted for the fatigue analysis of the plate. During the remodeling process, the changing apparent density and Young's modulus of bone are simulated in a step-wise fashion on the basis of Wolff's law, which is correlated with the specific clinical follow-up. The maximum biting forces were considered as the driving force on the bone remodeling, which are measured clinically at different time points (4, 16 and 28 months) after reconstruction surgery.

RESULTS

Under various occlusal loadings, the interaction between fatigue crack growth and bone remodeling is investigated to gain new insights for the future design of prosthetic devices. The simulation results reveal that appropriate remodeling of grafted bone could extend the fatigue life of fixation plates in a positive way. On the other hand, the rising occlusal load associated with healing and remodeling could lead to fatigue fracture of fixation plate and potentially cause severe bone resorption.

CONCLUSION

This study proposes an effective approach for more realistically predicting fatigue life of prosthetic devices subject to a tissue remodeling condition in-silico. It is anticipated to provide a guideline for deriving an optimal design of patient-specific prosthetic devices to better ensure longevity.

Mechanical failure of posterior teeth due to caries and occlusal wear- A modelling study

OBJECTIVES

The purpose of the present work is to explore the effect of occlusal wear and different types and degrees of caries on the mechanical performance and structural integrity of posterior teeth.

METHODS

Three-dimensional (3D) computational models with different combinations of caries parameters (caries location, caries size and caries induced pulp shrinkage) and occlusal wear factors (enamel thickness, marginal ridge height and cuspal slope) were developed and analyzed using the extended finite element method (XFEM) to identify the stress distribution, crack initiation load and ultimate fracture load values. The effect of a non-drilling conservative treatment using resin infiltration on the recovery of mechanical properties of carious molar teeth was also investigated.

RESULTS

Presence of fissural caries, worn proximal marginal ridge and decreased enamel thickness due to occlusal wear, imparted a significant negative effect on the crack initiation load value of the lower molar models. Accordingly, models with intact and strong proximal marginal ridge, generally exhibited higher crack initiation loading, regardless of caries size and location. Presence of fissure caries drastically decreased (55%-70%) the crack initiation load compared to sound teeth. The depth of the fissural lesion and the presence of proximal caries did not have a major effect on crack initiation load values. However, increasing the caries size resulted in lower final fracture load values in most of the cases. Accordingly, the groups with combined and connected large fissural and proximal lesions experienced the largest drop in the fracture load values compared to sound tooth models. The worst condition consisted of two connected large proximal and fissural caries with no proximal marginal ridge, in which the fracture load dramatically decreased to only 25% of that for sound teeth with intact marginal ridge. On the other hand, decreased cuspal slope due to occlusal wear and shrinkage of the pulp due to caries appeared to have a protective role and a direct relation with the fracture resistance of the tooth. Following the application of resin infiltration on the carious models, the crack initiation load and the fracture load could recover up to 75% and 90% of the values for the corresponding sound tooth models, respectively.

SIGNIFICANCE

Presence of fissural caries, if not treated (either with remineralization, resin infiltration or restoration), can be a major risk factor in the initiation of tooth fracture. When combined with decreased enamel thickness and loss of proximal marginal ridge due to mechanical or chemical wear, the weakening effect of the caries will be amplified specially in teeth with steep cuspal slopes. The application of a conservative treatment with resin infiltration can be an effective approach in prevention of further mechanical failure of demineralized enamel. The findings of this study emphasize the importance of early interventions in the management of caries for the prevention of future cuspal or tooth fracture especially in subjects with higher risk factors for tooth fracture such as caries, wear and bruxism.

Influence of Cavity Designs on Fracture Behavior of a Mandibular First Premolar with a Severely Curved h-shaped Canal

INTRODUCTION

This study aimed to investigate the fracture behavior of a mandibular first premolar with a severely curved h-shaped canal using the extended finite element method.

METHODS

Following the micro-computed tomographic data, models of the intact tooth, the conservative endodontic cavity (CEC), the modified conservative endodontic cavity (MCEC), and the traditional endodontic cavity (TEC) were created. All models were subjected with a total load of 600 N perpendicularly to the contact areas. The distributions of maximum principal stress were recorded. The evolution of cracks in the enamel and dentin was simulated with the extended finite element method.

RESULTS

The intact tooth showed the highest crack initiation load and the smallest stress concentration area. The CEC and MCEC showed higher crack initiation loads and smaller stress concentration areas compared with the TEC. On the occlusal surface, tensile stress was centralized around the distal fossa and the distal margins of cavities. In the root, the stress was concentrated at the mesiolingual side for all models and at the internal surface on the bifurcation section for the CEC. Cracks originated at buccal side of the distal fossa and microcracks were formed confined to the enamel, and then cracks occurred in the dentin below the bone level.

CONCLUSIONS

For the mandibular first premolar with a severely curved h-shaped canal, the MCEC preserved the fracture resistance equally as well as the CEC and reduced the stress concentration on the bifurcation section. The fracture initiated in the enamel, forming microcracks on the buccal side of the distal fossa and then occurred as an irreparable fracture in the dentin.

Relationship between Changes in Chemical Composition of Enamel Subsurface Lesions and the Emitted Nonlinear Optical Signals: An in vitro Study

The development of new diagnostic technologies based on the light scattering and autofluorescence properties of dental tissues is required to improve the diagnostic ability of initial caries lesions earlier than previously done and promoting the potential of treatment without surgical intervention. The aim of this study is to correlate fluorescence-based results provided by multiphoton microscopy (MPM) with confocal Raman microscopy records using phosphate level at 960 cm-1 and the organic matrix at ∼2,931 cm-1 in healthy and demineralized human enamel. Measurements on 14 teeth were made using two incident lights of different wavelengths, released by confocal Raman microscopy and MPM. Raman phosphate peak intensity at 960 cm-1 along with organic to mineral ratio at (2,931/430 cm-1) and nonlinear optical signals (second harmonic generation [SHG] and intrinsic two-photon excited fluorescence [I2PEF]) were recorded from the demineralized and healthy enamel sites. Raman spectral maps showed that the higher the organic/mineral ratio in the demineralized enamel, the lower the intensity of mineral component in the same zone. MPM revealed new optical indicators of carious lesion as shown by the presence of a red-shifted fluorescence peak in the 650- to 750-nm area of the fluorescence spectrum of demineralized enamel. Moreover, on sample regions with insignificant autofluorescence, the emergence of the SHG signal could be noted. By comparing I2PEF images with the structural motifs observed by the confocal Raman imaging system, the morphological similarity of the acquired images was quite evident. Any change in the I2PEF spectra reflects alterations in the chemical composition of enamel. These findings may provide an important basis for potentially valuable applications of photonic tools in the clinical diagnosis of tooth pathological conditions, besides exposing the fundamental role of organic matrix in enamel integrity and reparation.

Modelling of stress distribution and fracture in dental occlusal fissures

The aim of this study was to investigate the fracture behaviour of fissural dental enamel under simulated occlusal load in relation to various interacting factors including fissure morphology, cuspal angle and the underlying material properties of enamel. Extended finite element method (XFEM) was adopted here to analyse the fracture load and crack length in tooth models with different cusp angles (ranging from 50° to 70° in 2.5° intervals), fissural morphologies (namely U shape, V shape, IK shape, I shape and Inverted-Y shape) and enamel material properties (constant versus graded). The analysis results showed that fissures with larger curved morphology, such as U shape and IK shape, exhibit higher resistance to fracture under simulated occlusal load irrespective of cusp angle and enamel properties. Increased cusp angle (i.e. lower cusp steepness), also significantly enhanced the fracture resistance of fissural enamel, particularly for the IK and Inverted-Y shape fissures. Overall, the outcomes of this study explain how the interplay of compositional and structural features of enamel in the fissural area contribute to the resistance of the human tooth against masticatory forces. These findings may provide significant indicators for clinicians and technicians in designing/fabricating extra-coronal dental restorations and correcting the cuspal inclinations and contacts during clinical occlusal adjustment.

Numerical fatigue analysis of premolars restored by CAD/CAM ceramic crowns

OBJECTIVES

The purpose of this study was to estimate the fatigue life of premolars restored with two dental ceramics, lithium disilicate (LD) and polymer infiltrated ceramic (PIC) using the numerical method and compare it with the published in vitro data.

METHODS

A premolar restored with full-coverage crown was digitized. The volumetric shape of tooth tissues and crowns were created in Mimics^®. They were transferred to IA-FEMesh for mesh generation and the model was analyzed with Abaqus. By combining the stress distribution results with fatigue stress-life (S-N) approach, the lifetime of restored premolars was predicted.

RESULTS

The predicted lifetime was 1,231,318 cycles for LD with fatigue load of 1400N, while the one for PIC was 475,063 cycles with the load of 870N. The peak value of maximum principal stress occurred at the contact area (LD: 172MPa and PIC: 96MPa) and central fossa (LD: 100MPa and PIC: 64MPa) for both ceramics which were the most seen failure areas in the experiment. In the adhesive layer, the maximum shear stress was observed at the shoulder area (LD: 53.6MPa and PIC: 29MPa).

SIGNIFICANCE

The fatigue life and failure modes of all-ceramic crown determined by the numerical method seem to correlate well with the previous experimental study.

Identification of dynamic load for prosthetic structures

Dynamic load exists in numerous biomechanical systems, and its identification signifies a critical issue for characterizing dynamic behaviors and studying biomechanical consequence of the systems. This study aims to identify dynamic load in the dental prosthetic structures, namely, 3-unit implant-supported fixed partial denture (I-FPD) and teeth-supported fixed partial denture. The 3-dimensional finite element models were constructed through specific patient's computerized tomography images. A forward algorithm and regularization technique were developed for identifying dynamic load. To verify the effectiveness of the identification method proposed, the I-FPD and teeth-supported fixed partial denture structures were investigated to determine the dynamic loads. For validating the results of inverse identification, an experimental force-measuring system was developed by using a 3-dimensional piezoelectric transducer to measure the dynamic load in the I-FPD structure in vivo. The computationally identified loads were presented with different noise levels to determine their influence on the identification accuracy. The errors between the measured load and identified counterpart were calculated for evaluating the practical applicability of the proposed procedure in biomechanical engineering. This study is expected to serve as a demonstrative role in identifying dynamic loading in biomedical systems, where a direct in vivo measurement may be rather demanding in some areas of interest clinically.

Topological design of all-ceramic dental bridges for enhancing fracture resistance

Layered all-ceramic systems have been increasingly adopted in major dental prostheses. However, ceramics are inherently brittle, and they often subject to premature failure under high occlusion forces especially in the posterior region. This study aimed to develop mechanically sound novel topological designs for all-ceramic dental bridges by minimizing the fracture incidence under given loading conditions. A bi-directional evolutionary structural optimization (BESO) technique is implemented within the extended finite element method (XFEM) framework. Extended finite element method allows modeling crack initiation and propagation inside all-ceramic restoration systems. Following this, BESO searches the optimum distribution of two different ceramic materials, namely porcelain and zirconia, for minimizing fracture incidence. A performance index, as per a ratio of peak tensile stress to material strength, is used as a design objective. In this study, the novel XFEM based BESO topology optimization significantly improved structural strength by minimizing performance index for suppressing fracture incidence in the structures. As expected, the fracture resistance and factor of safety of fixed partial dentures structure increased upon redistributing zirconia and porcelain in the optimal topological configuration. Dental CAD/CAM systems and the emerging 3D printing technology were commercially available to facilitate implementation of such a computational design, exhibiting considerable potential for clinical application in the future. Copyright © 2015 John Wiley & Sons, Ltd.

Modeling of damage driven fracture failure of fiber post-restored teeth

Mechanical failure of biomaterials, which can be initiated by either violent force, or progressive stress fatigue, is a serious issue. Great efforts have been made to improve the mechanical performances of dental restorations. Virtual simulation is a promising approach for biomechanical investigations, which presents significant advantages in improving efficiency than traditional in vivo/in vitro studies. Over the past few decades, a number of virtual studies have been conducted to investigate the biomechanical issues concerning dental biomaterials, but only with limited incorporation of brittle failure phenomena. Motivated by the contradictory findings between several finite element analyses and common clinical observations on the fracture resistance of post-restored teeth, this study aimed to provide an approach using numerical simulations for investigating the fracture failure process through a non-linear fracture mechanics model. The ability of this approach to predict fracture initiation and propagation in a complex biomechanical status based on the intrinsic material properties was investigated. Results of the virtual simulations matched the findings of experimental tests, in terms of the ultimate fracture failure strengths and predictive areas under risk of clinical failure. This study revealed that the failure of dental post-restored restorations is a typical damage-driven continuum-to-discrete process. This approach is anticipated to have ramifications not only for modeling fracture events, but also for the design and optimization of the mechanical properties of biomaterials for specific clinically determined requirements.

Mechanics analysis of molar tooth splitting

A model for the splitting of teeth from wedge loading of molar cusps from a round indenting object is presented. The model is developed in two parts: first, a simple 2D fracture mechanics configuration with the wedged tooth simulated by a compact tension specimen; second, a full 3D numerical analysis using extended finite element modeling (XFEM) with an embedded crack. The result is an explicit equation for splitting load in terms of indenter radius and key tooth dimensions. Fracture experiments on extracted human molars loaded axially with metal spheres are used to quantify the splitting forces and thence to validate the model. The XFEM calculations enable the complex crack propagation, initially in the enamel coat and subsequently in the interior dentin, to be followed incrementally with increasing load. The fracture evolution is shown to be stable prior to failure, so that dentin toughness, not strength, is the controlling material parameter. Critical conditions under which tooth splitting in biological and dental settings are likely to be met, however rare, are considered.

Numerical Simulation of Biomechanical Behaviours in Novel Dental Restorations

Besides the prevention strategies against early stage dental caries, restoration is a preferable way to prevent decayed tooth from further deterioration. This study aimed to compare the mechanical strengths of carious tooth, traditionally restored tooth, and novel conservatively restored teeth under occlusal loading. The two-dimensional (2D) finite element method (FEM) was applied to quantify and compare maximum tensile stresses thereby predicting the initiation of crack. Taking into consideration of peak tensile stresses, it was found that the conservative (minimal intervention) restorations exhibited better fracture resistance than traditional restoration.

Thermally induced fracture for core-veneered dental ceramic structures

Effective and reliable clinical uses of dental ceramics necessitate an insightful analysis of the fracture behaviour under critical conditions. To better understand failure characteristics of porcelain veneered to zirconia core ceramic structures, thermally induced cracking during the cooling phase of fabrication is studied here by using the extended finite element method (XFEM). In this study, a transient thermal analysis of cooling is conducted first to determine the temperature distributions. The time-dependent temperature field is then imported to the XFEM model for viscoelastic thermomechanical analysis, which predicts thermally induced damage and cracking at different time steps. Temperature-dependent material properties are used in both transient thermal and thermomechanical analyses. Three typical ceramic structures are considered in this paper, namely bi-layered spheres, squat cylinders and dental crowns with thickness ratios of either 1:2 or 1:1. The XFEM fracture patterns exhibit good agreement with clinical observation and the in vitro experimental results obtained from scanning electron microscopy characterization. The study reveals that fast cooling can lead to thermal fracture of these different bi-layered ceramic structures, and cooling rate (in terms of heat transfer coefficient) plays a critical role in crack initiation and propagation. By exploring different cooling rates, the heat transfer coefficient thresholds of fracture are determined for different structures, which are of clear clinical implication.

An investigation of dentinal fluid flow in dental pulp during food mastication: simulation of fluid-structure interaction

This study uses fluid-structure interaction (FSI) simulation to investigate the relationship between the dentinal fluid flow in the dental pulp of a tooth and the elastic modulus of masticated food particles and to investigate the effects of chewing rate on fluid flow in the dental pulp. Three-dimensional simulation models of a premolar tooth (enamel, dentine, pulp, periodontal ligament, cortical bone, and cancellous bone) and food particle were created. Food particles with elastic modulus of 2,000 and 10,000 MPa were used, respectively. The external displacement loading (5 μm) was gradually directed to the food particle surface for 1 and 0.1 s, respectively, to simulate the chewing of food particles. The displacement and stress on tooth structure and fluid flow in the dental pulp were selected as evaluation indices. The results show that masticating food with a high elastic modulus results in high stress and deformation in the tooth structure, causing faster dentinal fluid flow in the pulp in comparison with that obtained with soft food. In addition, fast chewing of hard food particles can induce faster fluid flow in the pulp, which may result in dental pain. FSI analysis is shown to be a useful tool for investigating dental biomechanics during food mastication. FSI simulation can be used to predict intrapulpal fluid flow in dental pulp; this information may provide the clinician with important concept in dental biomechanics during food mastication.

Stress Analysis of a Class II MO-Restored Tooth Using a 3D CT-Based Finite Element Model

A computational method has been developed for stress analysis of a restored tooth so that experimental effort can be minimized. The objectives of this study include (i) developing a method to create a 3D FE assembly model for a restored tooth based on CT images and (ii) conducting stress analysis of the restored tooth using the 3D FE model established. To build up a solid computational model of a tooth, a method has been proposed to construct a 3D model from 2D CT-scanned images. Facilitated with CAD tools, the 3D tooth model has been virtually incorporated with a Class II MO restoration. The tooth model is triphasic, including the enamel, dentin, and pulp phases. To mimic the natural constraint on the movement of the tooth model, its corresponding mandible model has also been generated. The relative high maximum principal stress values were computed at the surface under loading and in the marginal region of the interface between the restoration and the tooth phases.

Computational Fracture Modelling in Bioceramic Structures

Bioceramics have rapidly emerged as one of major biomaterials in modern biomedical applications because of its outstanding biocompatibility. However, one drawback is its low tensile strength and fracture toughness due to brittleness and inherent microstructural defects, which to a certain extent prevents the ceramics from fully replacing metals used as load-bearing prostheses. This paper aims to model the crack initiation and propagation in ceramic fixed partial denture, namely dental bridge, by using two recently developed methods namely continuum-to-discrete element method (CDEM) in ELFEN and extended finite element methods (XFEM) in ABAQUS. Unlike most existing studies that typically required prescriptions of initial cracks, these two new approaches will model crack initiation and propagation automatically. They are applied to a typical prosthodontic example, thereby demonstrating their applicability and effectiveness in biomedical applications.

Mechanics of longitudinal cracks in tooth enamel

A study is made of longitudinal "channel" cracking in tooth enamel from axial compressive loading. The cracks simulate those generated in the molar and premolar teeth of humans and animals by natural tooth function. Contact loading tests are made on extracted human molars with hard and soft indenting plates to determine the evolution of such cracks with increasing load. Fracture is largely stable, with initial slow growth followed by acceleration as the cracks approach completion around an enamel side wall. A simple power law relation expresses the critical load for full fracture in terms of characteristic tooth dimensions-base radius and enamel thickness-as well as enamel toughness. Extended three-dimensional finite element modeling with provision for growth of embedded cracks is used to validate this relation. The cracks leave "fingerprints" that offer valuable clues to dietary habits, and provide a basis for a priori prediction of bite forces for different animals from measured tooth dimensions.

3D micro-crack propagation simulation at enamel/adhesive interface using FE submodeling and element death techniques

This study investigates micro-crack propagation at the enamel/adhesive interface using finite element (FE) submodeling and element death techniques. A three-dimensional (3D) FE macro-model of the enamel/adhesive/ceramic subjected to shear bond testing was generated and analyzed. A 3D micro-model with interfacial bonding structure was constructed at the upper enamel/adhesive interface where the stress concentration was found from the macro-model results. The morphology of this interfacial bonding structure (i.e., resin tag) was assigned based on resin tag geometry and enamel rod arrangement from a scanning electron microscopy micrograph. The boundary conditions for the micro-model were determined from the macro-model results. A custom iterative code combined with the element death technique was used to calculate the micro-crack propagation. Parallel experiments were performed to validate this FE simulation. The stress concentration within the adhesive occurred mainly at the upper corner near the enamel/adhesive interface and the resin tag base. A simulated fracture path was found at the resin tag base along the enamel/adhesive interface. A morphological observation of the fracture patterns obtained from in vitro testing corresponded with the simulation results. This study shows that the FE submodeling and element death techniques could be used to simulate the 3D micro-stress pattern and the crack propagation noted at the enamel/adhesive interface.

Reductions in the effects of damping on stress concentration in premolars by post-endodontic restorations: a non-linear finite element study

The aim of this study was to measure the structural damping constants of premolars after treatment with a cast Co-Cr post-core system or permanent root filling, and to evaluate the stress damping effects of these restored premolars. Both the damping ratio and the natural frequency (NF) of the cast Co-Cr post-core restored premolars and the permanent root-filled premolars were detected by in-vitro NF testing experiments. Unprepared premolars served as the control. The damping constants beta of the samples were calculated from the measured damping ratios and natural frequencies. The measured damping constants beta of the test premolars were then used for dynamic finite element (FE) analyses. Stress contours and damping effects of stresses in each treated type of premolar were computed and compared using ANSYS. The measured damping constants beta were 0.75 x 10(-5) for the unprepared premolars, 0.69 x 10(-5) for the root-filled premolars with coronal restoration, and 0.72 x 10(-5) for the cast Co-Cr post-core restored premolars. The unprepared intact premolars demonstrated the highest stress dissipation effects with a ratio of 29.3 per cent at the middle root opposite to the loading side. However, no stress dissipation effects were found in the premolars that had been restored with the cast Co-Cr post-core system. The FE analysis showed that metallic post treatment attenuated the damping properties of the premolar. The effects of damping on stress concentration were significantly lower in restored premolars than in untreated vital premolars. These findings suggest that future research on post material should take the damping property into consideration.

Predicting failure in mammalian enamel

Dentition is a vital element of human and animal function, yet there is little fundamental knowledge about how tooth enamel endures under stringent oral conditions. This paper describes a novel approach to the issue. Model glass dome specimens fabricated from glass and back-filled with polymer resin are used as representative of the basic enamel/dentine shell structure. Contact loading is used to deform the dome structures to failure, in simulation of occlusal loading with opposing dentition or food bolus. To investigate the role of enamel microstructure, additional contact tests are conducted on two-phase materials that capture the essence of the mineralized-rod/organic-sheath structure of dental enamel. These materials include dental glass-ceramics and biomimicked composites fabricated from glass fibers infiltrated with epoxy. The tests indicate how enamel is likely to deform and fracture along easy sliding and fracture paths within the binding phase between the rods. Analytical relations describing the critical loads for each damage mode are presented in terms of material properties (hardness, modulus, toughness) and tooth geometry variables (enamel thickness, cusp radius). Implications in dentistry and evolutionary biology are discussed.