A 3-D finite element study of the influence of crown-implant ratio on stress distribution

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Root canal filling: fracture strength of fiber-reinforced composite-restored roots and finite element analysis

The aims of this study were to evaluate the effect of root canal filling techniques on root fracture resistance and to analyze, by finite element analysis (FEA), the expansion of the endodontic sealer in two different root canal techniques. Thirty single-rooted human teeth were instrumented with rotary files to a standardized working length of 14 mm. The specimens were embedded in acrylic resin using plastic cylinders as molds, and allocated into 3 groups (n=10): G(lateral) - lateral condensation; G(single-cone) - single cone; G(tagger) - Tagger's hybrid technique. The root canals were prepared to a length of 11 mm with the #3 preparation bur of a tapered glass fiber-reinforced composite post system. All roots received glass fiber posts, which were adhesively cemented and a composite resin core was built. All groups were subjected to a fracture strength test (1 mm/min, 45°). Data were analyzed statistically by one-way ANOVA with a significance level of 5%. FEA was performed using two models: one simulated lateral condensation and Tagger's hybrid technique, and the other one simulated the single-cone technique. The second model was designed with an amount of gutta-percha two times smaller and a sealer layer two times thicker than the first model. The results were analyzed using von Mises stress criteria. One-way ANOVA indicated that the root canal filling technique affected the fracture strength (p=0.004). The G(lateral) and G(tagger) produced similar fracture strength values, while G(single-cone) showed the lowest values. The FEA showed that the single-cone model generated higher stress in the root canal walls. Sealer thickness seems to influence the fracture strength of restored endodontically treated teeth.

Zirconia-based dental crown to support a removable partial denture: a three-dimensional finite element analysis using contact elements and micro-CT data

Veneer fracture is the most common complication in zirconia-based restorations. The aim of this study was to evaluate the mechanical behavior of a zirconia-based crown in a lower canine tooth supporting removable partial denture (RPD) prosthesis, varying the bond quality of the veneer/coping interface. Microtomography (μCT) data of an extracted left lower canine were used to build the finite element model (M) varying the core material (gold core - MAu; zirconia core - MZi) and the quality of the veneer/core interface (complete bonded - MZi; incomplete bonded - MZi-NL). The incomplete bonding condition was only applied for zirconia coping by using contact elements (Target/Contact) with 0.3 frictional coefficients. Stress fields were obtained using Ansys Workbench 10.0. The loading condition (L = 1 N) was vertically applied at the base of the RPD prosthesis metallic support towards the dental apex. Maximum principal (σmax) and von Mises equivalent (σvM) stresses were obtained. The σmax (MPa) for the bonded condition was similar between gold and zirconia cores (MAu, 0.42; MZi, 0.40). The incomplete bonded condition (MZi-NL) raised σmax in the veneer up to 800% (3.23 MPa) in contrast to the bonded condition. The peak of σvM increased up to 270% in the MZi-NL. The incomplete bond condition increasing the stress in the veneer/zirconia interface.

Three-dimensional finite element analysis of anterior single implant-supported prostheses with different bone anchorages

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Straight and angulated abutments in platform switching: influence of loading on bone stress by three-dimensional finite element analysis

PURPOSE: In view of reports in the literature on the benefits achieved with the use of platform switching, described as the use of an implant with a larger diameter than the abutment diameter, the goal being to prevent the (previously) normal bone loss down to the first thread that occurs around most implants, thus enhancing soft tissue aesthetics and stability and the need for implant inclination due to bone anatomy in some cases, the aim of this study was to evaluate bone stress distribution on peri-implant bone, by using three-dimensional finite element analysis to simulate the influence of implants with different abutment angulations (0 and 15 degrees) in platform switching. METHODS: Four mathematical models of an implant-supported central incisor were created with varying abutment angulations: straight abutment (S1 and S2) and angulated abutment at 15 degrees (A1 and A2), submitted to 2 loading conditions (100 N): S1 and A1-oblique loading (45 degrees) and S2 and A2-axial loading, parallel to the long axis of the implant. Maximum (σmax) and minimum (σmin) principal stress values were obtained for cortical and trabecular bone. RESULTS: Models S1 and A1 showed higher σmax in cortical and trabecular bone when compared with S2 and A2. The highest σmax values (in MPa) in the cortical bone were found in S1 (28.5), followed by A1 (25.7), S2 (11.6), and A2 (5.15). For the trabecular bone, the highest σmax values were found in S1 (7.53), followed by A1 (2.87), S2 (2.85), and A2 (1.47). CONCLUSIONS: Implants with straight abutments generated the highest stress values in bone. In addition, this effect was potentiated when the load was applied obliquely.

Complex modal analysis of a vertical rotor through finite element method

Natural frequencies were analyzed (axial, torsional and flexural) and frequency response of a vertical rotor with a hard disk at the edge through the classical and complex modal analysis. The mathematical modeling was based on the theory of Euler-Bernoulli beam. The equation that rules the movement was obtained through the Lagrangian formulation. The model considered the effects of bending, torsion and axial deformation of the shaft, besides the gravitational and gyroscopic effects. The finite element method was used to discretize the structure into hollow cylindrical elements with 12 degrees of freedom. Mass, stiffness and gyroscopic matrices were explained consistently. This type of tool, based on the use of complex coordinates to describe the dynamic behavior of rotating shaft, allows the decomposition of the system in two submodes, backward and forward. Thus, it is possible to clearly visualize that the orbit and direction of the precessional motion around the line of the rotating shaft is not deformed. A finite element program was developed using Matlab ®, and numerical simulations were performed to validate this model.

Finite element analysis of bonded model Class I 'restorations' after shrinkage

Objectives. The C-Factor has been used widely to rationalize the changes in shrinkage stress occurring at the tooth/resin-composite interfaces. Experimentally, such stresses have been measured in a uniaxial direction between opposed parallel walls. The situation of adjoining cavity walls has been neglected. The aim was to investigate the hypothesis that: within stylized model rectangular cavities of constant volume and wall thickness, the interfacial shrinkage-stress at the adjoining cavity walls increases steadily as the C-Factor increases. Methods. Eight 3D-FEM restored Class I 'rectangular cavity' models were created by MSC.PATRAN/MSC.Marc, r2-2005 and subjected to 1% of shrinkage, while maintaining constant both the volume (20 mm(3)) and the wall thickness (2 mm), but varying the C-Factor (1.9-13.5). An adhesive contact between the composite and the teeth was incorporated. Polymerization shrinkage was simulated by analogy with thermal contraction. Principal stresses and strains were calculated. Peak values of maximum principal (MP) and maximum shear (MS) stresses from the different walls were displayed graphically as a function of C-Factor. The stress-peak association with C-Factor was evaluated by the Pearson correlation between the stress peak and the C-Factor. Results. The hypothesis was rejected: there was no clear increase of stress-peaks with C-Factor. The stress-peaks particularly expressed as MP and MS varied only slightly with increasing C-Factor. Lower stress-peaks were present at the pulpal floor in comparison to the stress at the axial walls. In general, MP and MS were similar when the axial wall dimensions were similar. The Pearson coefficient only expressed associations for the maximum principal stress at the ZX wall and the Z axis. Significance. Increase of the C-Factor did not lead to increase of the calculated stress-peaks in model rectangular Class I cavity walls. (C) 2011 Academy of Dental Materials. Published by Elsevier Ltd. All rights reserved.

Three-Dimensional Finite Element Analysis of Stress Distribution on Different Bony Ridges With Different Lengths of Morse Taper Implants and Prosthesis Dimensions

This finite element analysis (FEA) compared stress distribution on different bony ridges rehabilitated with different lengths of morse taper implants, varying dimensions of metal-ceramic crowns to maintain the occlusal alignment. Three-dimensional FE models were designed representing a posterior left side segment of the mandible: group control, 3 implants of 11 mm length; group 1, implants of 13 mm, 11 mm and 5 mm length; group 2, 1 implant of 11 mm and 2 implants of 5 mm length; and group 3, 3 implants of 5 mm length. The abutments heights were 3.5 mm for 13- and 11-mm implants (regular), and 0.8 mm for 5-mm implants (short). Evaluation was performed on Ansys software, oblique loads of 365N for molars and 200N for premolars. There was 50% higher stress on cortical bone for the short implants than regular implants. There was 80% higher stress on trabecular bone for the short implants than regular implants. There was higher stress concentration on the bone region of the short implants neck. However, these implants were capable of dissipating the stress to the bones, given the applied loads, but achieving near the threshold between elastic and plastic deformation to the trabecular bone. Distal implants and/or with biggest occlusal table generated greatest stress regions on the surrounding bone. It was concluded that patients requiring short implants associated with increased proportions implant prostheses need careful evaluation and occlusal adjustment, as a possible overload in these short implants, and even in regular ones, can generate stress beyond the physiological threshold of the surrounding bone, compromising the whole system.

Assessment of nose protector for sport activities: finite element analysis

There has been a significant increase in the number of facial fractures stemming from sport activities in recent years, with the nasal bone one of the most affected structures. Researchers recommend the use of a nose protector, but there is no standardization regarding the material employed. Clinical experience has demonstrated that a combination of a flexible and rigid layer of ethylene vinyl acetate (EVA) offers both comfort and safety to practitioners of sports. The aim of the present study was the investigation into the stresses generated by the impact of a rigid body on the nasal bone on models with and without an EVA protector. For such, finite element analysis was employed. A craniofacial model was constructed from images obtained through computed tomography. The nose protector was modeled with two layers of EVA (1 mm of rigid EVA over 2 mm of flexible EVA), following the geometry of the soft tissue. Finite element analysis was performed using the LS Dyna program. The bone and rigid EVA were represented as elastic linear material, whereas the soft tissues and flexible EVA were represented as hyperelastic material. The impact from a rigid sphere on the frontal region of the face was simulated with a constant velocity of 20 m s-1 for 9.1 mu s. The model without the protector served as the control. The distribution of maximal stress of the facial bones was recorded. The maximal stress on the nasal bone surpassed the breaking limit of 0.130.34 MPa on the model without a protector, while remaining below this limit on the model with the protector. Thus, the nose protector made from both flexible and rigid EVA proved effective at protecting the nasal bones under high-impact conditions.

FINITE ELEMENT ANALYSIS OF MAXILLARY BONE STRESS CAUSED BY ARAMANY CLASS IV OBTURATOR PROSTHESES

Statement of problem. The retention of an Aramany Class IV removable partial dental prosthesis can be compromised by a lack of support. The biomechanics of this obturator prosthesis result in an unusual stress distribution on the residual maxillary bone. Purpose. This study evaluated the biomechanics of an Aramany Class IV obturator prosthesis with finite element analysis and a digital 3-dimensional (3-D) model developed from a computed tomography scan; bone stress was evaluated according to the load placed on the prosthesis. Material and methods. A 3-D model of an Aramany Class IV maxillary resection and prosthesis was constructed. This model was used to develop a finite element mesh. A 120 N load was applied to the occlusal and incisal platforms corresponding to the prosthetic teeth. Qualitative analysis was based on the scale of maximum principal stress; values obtained through quantitative analysis were expressed in MPa. Results. Under posterior load, tensile and compressive stresses were observed; the tensile stress was greater than the compressive stress, regardless of the bone region, and the greatest compressive stress was observed on the anterior palate near the midline. Under an anterior load, tensile stress was observed in all of the evaluated bone regions; the tensile stress was greater than the compressive stress, regardless of the bone region. Conclusions. The Aramany Class IV obturator prosthesis tended to rotate toward the surgical resection when subjected to posterior or anterior loads. The amount of tensile and compressive stress caused by the Aramany Class IV obturator prosthesis did not exceed the physiological limits of the maxillary bone tissue. (J Prosthet Dent 2012;107:336-342)

GENERALIZED FINITE ELEMENT METHOD ON NONCONVENTIONAL HYBRID-MIXED FORMULATION

The generalized finite element method (GFEM) is applied to a nonconventional hybrid-mixed stress formulation (HMSF) for plane analysis. In the HMSF, three approximation fields are involved: stresses and displacements in the domain and displacement fields on the static boundary. The GFEM-HMSF shape functions are then generated by the product of a partition of unity associated to each field and the polynomials enrichment functions. In principle, the enrichment can be conducted independently over each of the HMSF approximation fields. However, stability and convergence features of the resulting numerical method can be affected mainly by spurious modes generated when enrichment is arbitrarily applied to the displacement fields. With the aim to efficiently explore the enrichment possibilities, an extension to GFEM-HMSF of the conventional Zienkiewicz-Patch-Test is proposed as a necessary condition to ensure numerical stability. Finally, once the extended Patch-Test is satisfied, some numerical analyses focusing on the selective enrichment over distorted meshes formed by bilinear quadrilateral finite elements are presented, thus showing the performance of the GFEM-HMSF combination.

Interpolation estimate for a finite-element space with embedded discontinuities

We consider a recently proposed finite-element space that consists of piecewise affine functions with discontinuities across a smooth given interface Γ (a curve in two dimensions, a surface in three dimensions). Contrary to existing extended finite element methodologies, the space is a variant of the standard conforming Formula space that can be implemented element by element. Further, it neither introduces new unknowns nor deteriorates the sparsity structure. It is proved that, for u arbitrary in Formula, the interpolant Formula defined by this new space satisfies Graphic where h is the mesh size, Formula is the domain, Formula, Formula, Formula and standard notation has been adopted for the function spaces. This result proves the good approximation properties of the finite-element space as compared to any space consisting of functions that are continuous across Γ, which would yield an error in the Formula-norm of order Graphic. These properties make this space especially attractive for approximating the pressure in problems with surface tension or other immersed interfaces that lead to discontinuities in the pressure field. Furthermore, the result still holds for interfaces that end within the domain, as happens for example in cracked domains.

Finite element analysis of stress distribution in anchor teeth in surgically assisted rapid palatal expansion

The treatment of a transverse maxillary deficiency in skeletally mature individuals should include surgically assisted rapid palatal expansion. This study evaluated the distribution of stresses that affect the expander's anchor teeth using finite element analysis when the osteotomy is varied. Five virtual models were built and the surgically assisted rapid palatal expansion was simulated. Results showed tension on the lingual face of the teeth and alveolar bone, and compression on the buccal side of the alveolar bone. The subtotal Le Fort I osteotomy combined with intermaxillary suture osteotomy seemed to reduce the dissipation of tensions. Therefore, subtotal Le Fort I osteotomy without a step in the zygomaticomaxillary buttress, combined with intermaxillary suture osteotomy and pterygomaxillary disjunction may be the osteotomy of choice to reduce tensions on anchor teeth, which tend to move mesiobuccally (premolar) and distobuccally (molar)

Morse taper implants at different bone levels: a finite element analysis of stress distribution

AIM: To explore the biomechanical effects of the different implantation bone levels of Morse taper implants, employing a finite element analysis (FEA). METHODS: Dental implants (TitamaxCM) with 4x13 mm and 4x11 mm, and their respective abutments with 3.5 mm height, simulating a screwed premolar metal-ceramic crown, had their design performed using the software AnsysWorkbench 10.0. They were positioned in bone blocks, covered by 2.5 mm thickness of mucosa. The cortical bone was designed with 1.5 mm thickness and the trabecular bone completed the bone block. Four groups were formed: group 11CBL (11 mm implant length on cortical bone level), group 11TBL (11 mm implant length on trabecular bone level), group 13CBL (13mm implant length on cortical bone level) and group 13TBL (13 mm implant length on trabecular bone level). Oblique 200 N loads were applied. Von Mises equivalent stresses in cortical and trabecular bones were evaluated with the same design program. RESULTS: The results were shown qualitatively and quantitatively by standard scales for each type of bone. By the results obtained, it can be suggested that positioning the implant completely in trabecular bone brings harm with respect to the generated stresses. Its implantation in the cortical bone has advantages with respect to better anchoring and locking, reflecting a better dissipation of the stresses along the implant/bone interfaces. In addition, the search for anchoring the implant in its apical region in cortical bone is of great value to improve stabilization and consequently better stress distribution. CONCLUSIONS: The implant position slightly below the bone in relation to the bone crest brings advantages as the best long-term predictability with respect to the expected neck bone loss.