Three-dimensional finite element analysis of MHD duct flow by the penalty function formulation

A numerical scheme based on the Finite Element Method (FEM) is presented to calculate the full solution of a three-dimensional steady magnetohydrodynamic (MHD) flow with moderately high Hartmann numbers and interaction parameters. An incompressible, viscous and electrically conducting liquid-metal is considered. Assuming a low magnetic Reynolds number, the solution method solves the coupled Navier-Stokes and Maxwell's equations through the use of a penalty function method. Results are presented for Hartmann numbers in the range 10(2)-10(3).

Finite element analysis of anisotropic optical waveguide with arbitrary index profile

This work presents the application of a scalar finite element formulation for Ex (TE-like) modes in anisotropic planar and channel waveguides with diagonal permittivity tensor, diffused in both transversal directions. This extended formulation considers explicitly both the variations of the refractive index and their spatial derivates inside of each finite element. Dispersion curves for Ex modes in planar and channel waveguides are shown, and the results compared with solutions obtained by other formulations.

Investigations into the use of the finite element method and artificial neural networks in the non-destructive analysis of metallic tubes

This work presents an investigation into the use of the finite element method and artificial neural networks in the identification of defects in industrial plants metallic tubes, due to the aggressive actions of the fluids contained by them, and/or atmospheric agents. The methodology used in this study consists of simulating a very large number of defects in a metallic tube, using the finite element method. Both variations in width and height of the defects are considered. Then, the obtained results are used to generate a set of vectors for the training of a perceptron multilayer artificial neural network. Finally, the obtained neural network is used to classify a group of new defects, simulated by the finite element method, but that do not belong to the original dataset. The reached results demonstrate the efficiency of the proposed approach, and encourage future works on this subject.

Electroosmotic pumping in rectangular microchannels: A numerical treatment by the finite element method

In this work, the analysis of electroosmotic pumping mechanisms in microchannels is performed through the solution of Poisson-Boltzmann and Navier Stokes equations by the Finite Element Method. This approach is combined with a Newton-Raphson iterative scheme, allowing a full treatment of the non-linear Poisson-Boltzmann source term which is normally approximated by linearizations in other methods.

Simulations of incompressible fluid flows by a least squares finite element method

In this work simulations of incompressible fluid flows have been done by a Least Squares Finite Element Method (LSFEM) using velocity-pressure-vorticity and velocity-pressure-stress formulations, named u-p-ω) and u-p-τ formulations respectively. These formulations are preferred because the resulting equations are partial differential equations of first order, which is convenient for implementation by LSFEM. The main purposes of this work are the numerical computation of laminar, transitional and turbulent fluid flows through the application of large eddy simulation (LES) methodology using the LSFEM. The Navier-Stokes equations in u-p-ω and u-p-τ formulations are filtered and the eddy viscosity model of Smagorinsky is used for modeling the sub-grid-scale stresses. Some benchmark problems are solved for validate the numerical code and the preliminary results are presented and compared with available results from the literature. Copyright © 2005 by ABCM.

Finite element analysis of warpage in laminated aluminium alloy plates for machining of primary aeronautic parts

The aim of this paper consists in presenting a method of simulating the warpage in 7xxx series aluminium alloy plates. To perform this simulation finite element software MSC.Patran and MSC.Marc were used. Another result of this analysis will be the influence on material residual stresses induced on the raw material during the rolling process upon the warpage of primary aeronautic parts, fabricated through machining (milling) at Embraer. The method used to determinate the aluminium plate residual stress was Layer Removal Test. The numerical algorithm Modified Flavenot Method was used to convert layer removal and beam deflection in stress level. With such information about the level and profile of residual stresses become possible, during the step that anticipate the manufacturing to incorporate these values in the finite-element approach for modelling warpage parts. Based on that warpage parameter surely the products are manufactured with low relative vulnerability propitiating competitiveness and price. © 2007 American Institute of Physics.

Simulation of needle-type corona electrodes by the finite element method

This paper describes a software tool, called LEVSOFT, suitable for the electric field simulations of corona electrodes by the Finite Element Method (FEM). Special attention was paid to the user friendly construction of geometries with corners and sharp points, and to the fast generation of highly refined triangular meshes and field maps. The execution of self-adaptive meshes was also implemented. These customized features make the code attractive for the simulation of needle-type corona electrodes. Some case examples involving needle type electrodes are presented.

Study of the phenomena of fatigue using the finite element method

About 99% of mechanical failures are consequence of the phenomena of fatigue, which consists on the progressive weakening of the resistant section of a mechanical component due to the growing of cracks caused by fluctuating loadings. A broad diversity of factors influences the fatigue life of a mechanical component, like the surface finishing, scale factors, among others, but none is as significantly as the presence of geometric severities. Stress concentrators are places where fatigue cracks have a greater probability to occur, and so on, the intuit of this work is to develop a consistent and trustfully methodology to determine the theoretical stress concentration factor of mechanical components. Copyright © 2007 SAE International.

Stress distribution in ceramic restorations over natural tooth using finite element analysis. lithium disilicate x aluminum oxide material

Background: Data on stress distribution in tooth-restoration interface with different ceramic restorative materials are limited. The aim of this chapter was to assess the stress distribution in the interface of ceramic restorations with laminate veneer or full-coverage crown with two different materials (lithium dissilicate and densely sintered aluminum oxide) under different loading areas through finite element analysis. Materials and Methods: Six two-dimensional finite element models were fabricated with different restorations on natural tooth: laminate veneer (IPS Empress, IPS Empress Esthetic and Procera AllCeram) or full-coverage crown (IPS e.max Press and Procera AllCeram). Two different loading areas (L) (50N) were also determined: palatal surface at 45° in relation to the long axis of tooth (L1) and perpendicular to the incisal edge (L2). A model with higid natural tooth was used as control. von Mises equivalent stress (σ vM) and maximum principal stress (σ max) were obtained on Ansys software. Results: The presence of ceramic restoration increased σ vM and σ max in the adhesive interface, mainly for the aluminum oxide (Procera AllCeram system) restorations. The full-coverage crowns generated higher stress in the adhesive interface under L1 while the same result was observed for the laminate veneers under L2. Conclusions: Lithium dissilicate and densely sintered aluminum oxide restorations exhibit different behavior due to different mechanical properties and loading conditions. © 2011 Nova Science Publishers, Inc.

Micromechanics of dentin/adhesive interface in function of dentin depth: 3d finite element analysis

Objectives: The aim of this study was to analyze the stress distribution on dentin/adhesive interface (d/a) through a 3-D finite element analysis (FEA) varying the number and diameter of the dentin tubules orifice according to dentin depth, keeping hybrid layer (HL) thickness and TAǴs length constant. Materials and Methods: 3 models were built through the SolidWorks software: SD - specimen simulating superficial dentin (41 x 41 x 82 μm), with a 3 μm thick HL, a 17 μm length Tag, and 8 tubules with a 0.9 μm diameter restored with composite resin. MD - similar to M1 with 12 tubules with a 1.2 μm diameter, simulating medium dentin. DD - similar to M1 with 16 tubules with a 2.5 μm diameter, simulating deep dentin. Other two models were built in order to keep the diameter constant in 2.5 μm: MS - similar to SD with 8 tubules; and MM - similar to MD with 12 tubules. The boundary condition was applied to the base surface of each specimen. Tensile load (0.03N) was performed on the composite resin top surface. Stress field (maximum principal stress in tension - σMAX) was performed using Ansys Wokbench 10.0. Results: The peak of σMAX (MPa) were similar between SD (110) and MD (106), and higher for DD (134). The stress distribution pathway was similar for all models, starting from peritubular dentin to adhesive layer, intertubular dentin and hybrid layer. The peak of σMAX (MPa) for those structures was, respectively: 134 (DD), 56.9 (SD), 45.5 (DD), and 36.7 (MD). Conclusions: The number of dentin tubules had no influence in the σMAX at the dentin/adhesive interface. Peritubular and intertubular dentin showed higher stress with the bigger dentin tubules orifice condition. The σMAX in the hybrid layer and adhesive layer were going down from superficial dentin to deeper dentin. In a failure scenario, the hybrid layer in contact with peritubular dentin and adhesive layer is the first region for breaking the adhesion. © 2011 Nova Science Publishers, Inc.

Influence of buccal cusp reduction when using porcelain laminate veneers in premolars. A comparative study using 3-D finite element analysis

Objectives: Based on a maxillary premolar restored with laminate veneer and using the 3-D finite element analysis (FEA) and mCT data, the aim of this study was to evaluate the influence of different types of buccal cusp reduction on the stress distribution in the porcelain laminate veneer and in the resin luting cement layer. Methods: Two 3-D FEA models (M) of a maxillary premolar were built from mCT data. The buccal cusp reduction followed two configurations: Mt-buccal cusp completely covered by porcelain laminate veneer; and Mp-buccal cusp partially covered by porcelain laminate veneer. The loading (150 N in 458) was performed on the top of the buccal cusp. The finite element software (Ansys Workbench 10.0) was used to obtain the maximum shear stress (σmax) and maximum principal stress (σmax). Results: The Mp showed reduced the stress (σmax) in porcelain laminate veneer (from-2.3 to 24.5 MPa) in comparison with Mt (from-5.3 to 27.4 MPa). The difference between the peak and lower stress values of σmax in Mp (-6.8 to 26.7 MPa) and Mt (-5.3 to 27.4 MPa) was similar for the resin luting cement layer. The structures not exceeded the ultimate tensile strength or the shear bond strength. Conclusions: Cusp reduction did not affect significant increase in σmax and τmax. The Mt showed better stress distribution (τmax) than Mp. © 2011 Published by Elsevier Ireland on behalf of Japan Prosthodontic Society.

Comparison of different designs of implant-retained overdentures and fixed full-arch implant-supported prosthesis on stress distribution in edentulous mandible - A computed tomography-based three-dimensional finite element analysis

A finite element analysis was used to compare the effect of different designs of implant-retained overdentures and fixed full-arch implant-supported prosthesis on stress distribution in edentulous mandible. Four models of an human mandible were constructed. In the OR (O'ring) group, the mandible was restored with an overdenture retained by four unsplinted implants with O'ring attachment; in the BC (bar-clip) -C and BC groups, the mandibles were restored with overdentures retained by four splinted implants with bar-clip anchor associated or not with two distally placed cantilevers, respectively; in the FD (fixed denture) group, the mandible was restored with a fixed full-arch four-implant-supported prosthesis. Models were supported by the masticatory muscles and temporomandibular joints. A 100-N oblique load was applied on the left first molar. Von Mises (σvM), maximum (σmax) and minimum (σmin) principal stresses (in MPa) analyses were obtained. BC-C group exhibited the highest stress values (σvM=398.8, σmax=580.5 and σmin=-455.2) while FD group showed the lowest one (σvM=128.9, σmax=185.9 and σmin=-172.1). Within overdenture groups, the use of unsplinted implants reduced the stress level in the implant/prosthetic components (59.4% for σvM, 66.2% for σmax and 57.7% for σmin versus BC-C group) and supporting tissues (maximum stress reduction of 72% and 79.5% for σmax, and 15.7% and 85.7% for σmin on the cortical and trabecular bones, respectively). Cortical bone exhibited greater stress concentration than the trabecular bone for all groups. The use of fixed implant dentures and removable dentures retained by unsplinted implants to rehabilitate edentulous mandible reduced the stresses in the periimplant bone tissue, mucosa and implant/prosthetic components. © 2013 Elsevier Ltd.

Free tools and strategies for the generation of 3D finite element meshes: Modeling of the cardiac structures

The Finite Element Method is a well-known technique, being extensively applied in different areas. Studies using the Finite Element Method (FEM) are targeted to improve cardiac ablation procedures. For such simulations, the finite element meshes should consider the size and histological features of the target structures. However, it is possible to verify that some methods or tools used to generate meshes of human body structures are still limited, due to nondetailed models, nontrivial preprocessing, or mainly limitation in the use condition. In this paper, alternatives are demonstrated to solid modeling and automatic generation of highly refined tetrahedral meshes, with quality compatible with other studies focused on mesh generation. The innovations presented here are strategies to integrate Open Source Software (OSS). The chosen techniques and strategies are presented and discussed, considering cardiac structures as a first application context. © 2013 E. Pavarino et al.

A finite element approach for predicting the ultimate rotation capacity of RC beams

This paper presents a numerical approach to model the complex failure mechanisms that define the ultimate rotational capacity of reinforced concrete beams. The behavior in tension and compression is described by a constitutive damage model derived from a combination of two specific damage models [1]. The nonlinear behavior of the compressed region is treated by the compressive damage model based on the Drucker-Prager criterion written in terms of the effective stresses. The tensile damage model employs a failure criterion based on the strain energy associated with the positive part the effective stress tensor. This model is used to describe the behavior of very thin bands of strain localization, which are embedded in finite elements to represent multiple cracks that occur in the tensioned region [2]. The softening law establishes dissipation energy compatible with the fracture energy of the concrete. The reinforcing steel bars are modeled by truss elements with elastic-perfect plastic behavior. It is shown that the resulting approach is able to predict the different stages of the collapse mechanism of beams with distinct sizes and reinforcement ratios. The tensile damage model and the finite element embedded crack approach are able to describe the stiffness reduction due to concrete cracking in the tensile zone. The truss elements are able to reproduce the effects of steel yielding and, finally, the compressive damage model is able to describe the non-linear behavior of the compressive zone until the complete collapse of the beam due to crushing of concrete. The proposed approach is able to predict well the plastic rotation capacity of tested beams [3], including size-scale effects.

Influence of the implant diameter with different sizes of hexagon: Analysis by 3-dimensional finite element method

The aim of this study was to evaluate the stress distribution in implants of regular platforms and of wide diameter with different sizes of hexagon by the 3-dimensional finite element method. We used simulated 3-dimensional models with the aid of Solidworks 2006 and Rhinoceros 4.0 software for the design of the implant and abutment and the InVesalius software for the design of the bone. Each model represented a block of bone from the mandibular molar region with an implant 10 mm in length and different diameters. Model A was an implant 3.75 mm/regular hexagon, model B was an implant 5.00 mm/regular hexagon, and model C was an implant 5.00 mm/ expanded hexagon. A load of 200 N was applied in the axial, lateral, and oblique directions. At implant, applying the load (axial, lateral, and oblique), the 3 models presented stress concentration at the threads in the cervical and middle regions, and the stress was higher for model A. At the abutment, models A and B showed a similar stress distribution, concentrated at the cervical and middle third; model C showed the highest stresses. On the cortical bone, the stress was concentrated at the cervical region for the 3 models and was higher for model A. In the trabecular bone, the stresses were less intense and concentrated around the implant body, and were more intense for model A. Among the models of wide diameter (models B and C), model B (implant 5.00 mm/regular hexagon) was more favorable with regard to distribution of stresses. Model A (implant 3.75 mm/regular hexagon) showed the largest areas and the most intense stress, and model B (implant 5.00 mm/regular hexagon) showed a more favorable stress distribution. The highest stresses were observed in the application of lateral load.