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.

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.

Dynamic stationary response of reinforced plates by the boundary element method

A direct version of the boundary element method (BEM) is developed to model the stationary dynamic response of reinforced plate structures, such as reinforced panels in buildings, automobiles, and airplanes. The dynamic stationary fundamental solutions of thin plates and plane stress state are used to transform the governing partial differential equations into boundary integral equations (BIEs). Two sets of uncoupled BIEs are formulated, respectively, for the in-plane state ( membrane) and for the out-of-plane state ( bending). These uncoupled systems are joined to formamacro-element, in which membrane and bending effects are present. The association of these macro-elements is able to simulate thin-walled structures, including reinforced plate structures. In the present formulation, the BIE is discretized by continuous and/or discontinuous linear elements. Four displacement integral equations are written for every boundary node. Modal data, that is, natural frequencies and the corresponding mode shapes of reinforced plates, are obtained from information contained in the frequency response functions (FRFs). A specific example is presented to illustrate the versatility of the proposed methodology. Different configurations of the reinforcements are used to simulate simply supported and clamped boundary conditions for the plate structures. The procedure is validated by comparison with results determined by the finite element method (FEM).

Analysis of mechanical behavior in bending of polymer composites using the finite elements method

The use of composite materials has increased in the recent decades, mainly in the aeronautics and automotives industries. In the present study is elaborated a computational simulation program of the bending test using the finite elements method, in the commercial software ANSYS. This simulation has the objective of analyze the mechanical behavior in bending of two composites with polymeric matrix reinforced with carbon fibers. Also are realized bending tests of the 3 points to obtain the resistances of the materials. Data from simulation and tests are used to make a comparison between two failures criteria, Tsai-Wu and Hashin criterion. Copyright © 2009 SAE International.

Complex modal analysis of a slender vertical rotor by finite elements method

In this paper, 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 modal and complex analysis. 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. The classical modal analysis, usually applied to stationary structures, does not consider an important characteristic of rotating machinery which are the methods of forward and backward whirl. Initially, through the traditional modal analysis, axial and torsional natural frequencies were obtained in a static shaft, since they do not suffer the influence of gyroscopic effects. Later research was performed by complex modal analysis. 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 (TM) and numerical simulations were performed to validate this model. Natural frequencies and directional frequency forced response (dFRF) were obtained using the complex modal analysis for a simple vertical rotor and also for a typical drill string used in the construction of oil wells.

Wavelet-Galerkin method for one-dimensional elastoplasticity and damage problems: Constitutive modeling and computational aspects

This work presents an analysis of the wavelet-Galerkin method for one-dimensional elastoplastic-damage problems. Time-stepping algorithm for non-linear dynamics is presented. Numerical treatment of the constitutive models is developed by the use of return-mapping algorithm. For spacial discretization we can use wavelet-Galerkin method instead of standard finite element method. This approach allows to locate singularities. The discrete formulation developed can be applied to the simulation of one-dimensional problems for elastic-plastic-damage models. (C) 2007 Elsevier B.V. All rights reserved.

An application of the Adomian's decomposition method to the analysis of MHD duct flows

A new approach is proposed in this work for the treatment of boundary value problems through the Adomian's decomposition method. Although frequently claimed as accurate and having fast convergence rates, the original formulation of Adomian's method does not allow the treatment of homogeneous boundary conditions along closed boundaries. The technique here presented overcomes this difficulty, and is applied to the analysis of magnetohydrodynamic duct flows. Results are in good agreement with finite element method calculations and analytical solutions for square ducts. Therefore, new possibilities appear for the application of Adomian's method in electromagnetics.

Discrete generator coordinate: Symmetry restoration in finite-dimensional spaces

Group theoretical-based techniques and fundamental results from number theory are used in order to allow for the construction of exact projectors in finite-dimensional spaces. These operators are shown to make use only of discrete variables, which play the role of discrete generator coordinates, and their application in the number symmetry restoration is carried out in a nuclear BCS wave function which explicitly violates that symmetry. © 1999 Published by Elsevier Science B.V. All rights reserved.

Induced variational method from supersymmetric quantum mechanics and the screened Coulomb potential

The formalism of supersymmetric quantum mechanics supplies a trial wave function to be used in the variational method. The screened Coulomb potential is analyzed within this approach. Numerical and exact results for energy eigenvalues are compared.

A general technique to embed non-uniform discontinuities into standard solid finite elements

A general technique to embed non-uniform displacement discontinuities into standard solid finite elements is presented. The technique is based on the decomposition of the kinematic fields into a component related to the deformation of the solid portion of the element and one related to the rigid-body motion due to a displacement discontinuity. This decomposition simplifies the incorporation of discontinuity interfaces and provides a suitable framework to account for non-uniform discontinuity modes. The present publication addresses two families of finite element formulations: displacement-based and stress hybrid finite element. © 2005 Elsevier Ltd. All rights reserved.