Application of hte finite element method in Dentistry

Introduction: The finite element method (FEM) involves a series of computational procedures to calculate the stress in each element, which performs a model solution. Such a structural analysis allows the determination of stress resulting from external force, pressure, thermal change, and other factors. This method is extremely useful for indicating mechanical aspects of biomaterials and human tissues that can hardly be measured in vivo. The results obtained can then be studied using visualization software within the FEM environment to view a variety of parameters, and to fully identify implications of the analysis. Objective: An overview to show application of FEM in dentistry was undertaken. Literature review: This paper shows the basic concept, advances, advantages, limitations and applications of finite element method (FEM) in dentistry. Conclusion: It is extremely important to verify what the purpose of the study is in order to correctly apply FEM.

Scientific use of the finite element method in Orthodontics

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Stress distribution in delayed replanted teeth splinted with different orthodontic wires: a three-dimensional finite element analysis

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

Finite element analysis of the influence of geometry and design of zirconia crowns on stress distribution

Purpose: To evaluate the influence of the geometry and design of prosthetic crown preparations on stress distribution in compression tests, using finite element analysis (FEA). Materials and Methods: Six combinations of 3D drawings of all-ceramic crowns (yttria-stabilized zirconia framework and porcelain veneer) were evaluated: F, flat preparation and simplified crown; FC, flat preparation and crown with contact point; FCM, flat preparation and modified crown; A, anatomical preparation and simplified anatomical crown framework; AC, anatomical preparation and crown with contact point; and ACM, anatomical preparation and modified crown. Bonded contact types at all interfaces with the mesh were assigned, and the material properties used were according to the literature. A 200 N vertical load was applied at the center of each model. The maximum principal stresses were quantitatively and qualitatively analyzed. Results: The highest values of tensile stress were observed at the interface between the ceramics in the region under the load application for the simplified models (F and A). Reductions in stress values were observed for the model with the anatomical preparation and modified infrastructure (ACM). The stress distribution in the flat models was similar to that of their respective anatomical models. Conclusions: The modified design of the zirconia coping reduces the stress concentration at the interface with the veneer ceramic, and the simplified preparation can exert a stress distribution similar to that of the anatomical preparation at and near the load point, when load is applied to the center of the crown.

Stress distribution around osseointegrated implants with different internal-cone connections: photoelastic and finite element analysis

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Fatigue surviving, fracture resistance, shear stress and finite element analysis of glass fiber posts with different diameters

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

A coupling technique for non-matching finite element meshes

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

Computer graphics applied to anatomy: a study of two bio-cad modeling methods on finite element analysis of human edentulous hemi-mandible

Modeling is a step to perform a finite element analysis. Different methods of model construction are reported in literature, as the Bio-CAD modeling. The purpose of this study was to perform a model evaluation and application using two methods of Bio-CAD modeling from human edentulous hemi-mandible on the finite element analysis. From CT scans of dried human skull was reconstructed a stereolithographic model. Two methods of modeling were performed: STL conversion approach (Model 1) associated to STL simplification and reverse engineering approach (Model 2). For finite element analysis was used the action of lateral pterygoid muscle as loading condition to assess total displacement (D), equivalent von-Mises stress (VM) and maximum principal stress (MP). Two models presented differences on the geometry regarding surface number (1834 (model 1); 282 (model 2)). Were observed differences in finite element mesh regarding element number (30428 nodes/16683 elements (model 1); 15801 nodes/8410 elements (model 2). D, VM and MP stress areas presented similar distribution in two models. The values were different regarding maximum and minimum values of D (ranging 0-0.511 mm (model 1) and 0-0.544 mm (model 2), VM stress (6.36E-04-11.4 MPa (model 1) and 2.15E-04-14.7 MPa (model 2) and MP stress (-1.43-9.14 MPa (model 1) and -1.2-11.6 MPa (model 2). From two methods of Bio-CAD modeling, the reverse engineering presented better anatomical representation compared to the STL conversion approach. The models presented differences in the finite element mesh, total displacement and stress distribution.

Effect of crown-to-implant ratio on peri-implant stress: a finite element analysis

The aim of this study was to evaluate stress distribution in the fixation screws and bone tissue around implants in single-implant supported prostheses with crowns of different heights (10,12.5, 15 mm crown-to-implant ratio 1:1, 1.25:1, 1.5:1, respectively). It was designed using three 3-Dmodels. Each model was developed with a mandibular segment of bone block including an internal hexagon implant supporting a screw-retained, single metalceramic crown. The crown height was set at 10, 12.5, and 15 mm with crown-to-implant ratio of 1:1, 1.25:1, 1.5:1, respectively. The applied forces were 200 N (axial) and 100 N (oblique). The increase of crown height showed differences with the oblique load in some situations. By von Mises'criterion, a high stress area was concentrated at the implant/fixation screw and abutment/implant interfaces at crown-to-implant ratio of 1:1, 1.25:1, 1.5:1, respectively. Using the maxiinum principal criteria, the buccal regions showed higher traction stress intensity, whereas the distal regions showed the largest compressive stress in all models. The increase of C/I ratio must be carefully evaluated by the dentist since the increase of this C/I ratio is proportional to the increase of average stress for both screw fixation (C/I 1:1 to 1:1.25 ratio = 30.1% and C/I 1:1 to 1 :1.5 ratio = 46.3%) and bone tissue (C/I 1:1 to 1:1.25 ratio = 30% and C/I 1:1 to 1:1.5 ratio = 51.5%). (C) 2014 Elsevier B.V. All rights reserved.

Dynamic analysis for different finite element models of beams with viscoelastic damping layer

Many new viscoelastic materials have been developed recently to help improve noise and vibration levels in mechanical structures for applications in automobile and aeronautical industry. The viscoelastic layer treatment applied to solid metal structures modifies two main properties which are related to the mass distribution and the damping mechanism. The other property controlling the dynamics of a mechanical system is the stiffness that does not change much with the viscoelastic material. The model of such system is usually complex, because the viscoelastic material can exhibit nonlinear behavior, in contrast with the many available tools for linear dynamics. In this work, the dynamic behavior of sandwich beam is modeled by finite element method using different element types which are then compared with experimental results developed in the laboratory for various beams with different viscoelastic layer materials. The finite element model is them updated to help understand the effects in the damping for various natural frequencies and the trade-off between attenuation and the mass add to the structure.

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)

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.