MATHEMATICA notebook for computing tetrahedral finite element shape functions and matrices for the Helmholtz equation

A MATHEMATICA notebook to compute the elements of the matrices which arise in the solution of the Helmholtz equation by the finite element method (nodal approximation) for tetrahedral elements of any approximation order is presented. The results of the notebook enable a fast computational implementation of finite element codes for high order simplex 3D elements reducing the overheads due to implementation and test of the complex mathematical expressions obtained from the analytical integrations. These matrices can be used in a large number of applications related to physical phenomena described by the Poisson, Laplace and Schrodinger equations with anisotropic physical properties.

Finite element analysis and fracture resistance testing of a new intraradicular post

Objectives: The objective of the present study was to evaluate a prefabricated intraradicular threaded pure titanium post, designed and developed at the Sao Jose dos Campos School of Dentistry - UNESP, Brazil. This new post was designed to minimize stresses observed with prefabricated post systems and to improve cost-benefits. Materials and and methods: Fracture resistance testing of the post/core/root complex, fracture analysis by microscopy and stress analysis by the finite element method were used for post evaluation. The following four prefabricated metal post systems were analyzed: group 1, experimental post; group 2, modification of the experimental post; group 3, Flexi Post, and group 4, Para Post. For the analysis of fracture resistance, 40 bovine teeth were randomly assigned to the four groups (n=10) and used for the fabrication of test specimens simulating the situation in the mouth. The test specimens were subjected to compressive strength testing until fracture in an EMIC universal testing machine. After fracture of the test specimens, their roots were sectioned and analyzed by microscopy. For the finite element method, specimens of the fracture resistance test were simulated by computer modeling to determine the stress distribution pattern in the post systems studied. Results: The fracture test presented the following averages and standard deviation: G1 (45.63 +/- 8.77), G2 (49.98 +/- 7.08), G3 (43.84 +/- 5.52), G4 (47.61 +/- 7.23). Stress was homogenously distributed along the body of the intraradicular post in group 1, whereas high stress concentrations in certain regions were observed in the other groups. These stress concentrations in the body of the post induced the same stress concentration in root dentin. Conclusions: The experimental post (original and modified versions) presented similar fracture resistance and better results in the stress analysis when compared with the commercial post systems tested (08/2008PA/CEP).

The role of chain scission and chain branching in high density polyethylene during thermo-mechanical degradation

The mechanical and thermo-oxidative degradation of high density polyethylene (HDPE) was measured in a twin-screw extruder using various processing conditions. Two types of HDPE, Phillips and Ziegler-Natta, having different levels of terminal vinyl unsaturation were analysed. Mild screw profiles, having mainly conveying elements, have short mean residence times then profiles with kneading discs and left hand elements. Carbonyl and traps-vinylene group concentrations increased, whereas vinyl group concentration decreased with number of extrusions. Higher temperature profiles intensified these effects. The thermo-mechanical degradation mechanism begins with chain scission in the longer chains due to their higher probability of entanglements. These macroradicals then react with the vinyl terminal unsaturations of other chains producing chain branching. Shorter chains are more mobile, not suffering scission but instead are used for grafting the macroradicals, increasing the molecular weight. Increase in the levels of extrusion temperature, shear and vinyl end groups content facilitates the thermo-mechanical degradation reducing the amount of both, longer chains via chain scission and shorter chains via chain branching, narrowing the polydispersity. Phillips HDPE produces a higher level of chain branching than does the Ziegler-Natta type. (C) 2004 Elsevier Ltd. All rights reserved.

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.

Study of the thermo-mechanical and electrical properties of conducting composites containing natural rubber and carbon black

This work shows the preparation and characterization of composites obtained by mixing natural rubber (NR) and carbon black (CB) in different percentages aiming suitable mechanical properties, processability and electrical conductivity for future applications as transducers in pressure sensors. The composites NR/CB are characterized through dc conductivity, thermal analysis using differential scanning calorimetry (DSC), dynamic mechanical thermal analysis (DMA), thermogravimetry (TGA) and stress-strain test. The electrical conductivity changed from 10-9 to 10 Sm-1 depending on the percentage of CB in the composite. Besides, it was found a linear (and reversible) dependence of the conductivity on the applied pressure in the range from 0 to 1.6 MPa for the sample 80/20 (NR/CB wt%).

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.

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.