Finite element aided design evolution of composite leaf spring

This paper shows the process of the virtual production development of the mechanical connection between the top leaf of a dual composite leaf spring system to a shackle using finite element methods. The commercial FEA package MSC/MARC has been used for the analysis. In the original design the joint was based on a closed eye-end. Full scale testing results showed that this configuration achieved the vertical proof load of 150 kN and 1 million cycles of fatigue load. However, a problem with delamination occurred at the interface between the fibres going around the eye and the main leaf body. To overcome this problem, a second design was tried using transverse bandages of woven glass fibre reinforced tape to wrap the section that is prone to delaminate. In this case, the maximum interlaminar shear stress was reduced by a certain amount but it was still higher than the material’s shear strength. Based on the fact that, even with delamination, the top leaf spring still sustained the maximum static and fatigue loads required, the third design was proposed with an open eye-end, eliminating altogether the interface where the maximum shear stress occurs. The maximum shear stress predicted by FEA is reduced significantly and a safety factor of around 2 has been obtained. Thus, a successful and safe design has been achieved.

An Exponentially Convergent Nonpolynomial Finite Element Method for Time-Harmonic Scattering from Polygons

In recent years nonpolynomial finite element methods have received increasing attention for the efficient solution of wave problems. As with their close cousin the method of particular solutions, high efficiency comes from using solutions to the Helmholtz equation as basis functions. We present and analyze such a method for the scattering of two-dimensional scalar waves from a polygonal domain that achieves exponential convergence purely by increasing the number of basis functions in each element. Key ingredients are the use of basis functions that capture the singularities at corners and the representation of the scattered field towards infinity by a combination of fundamental solutions. The solution is obtained by minimizing a least-squares functional, which we discretize in such a way that a matrix least-squares problem is obtained. We give computable exponential bounds on the rate of convergence of the least-squares functional that are in very good agreement with the observed numerical convergence. Challenging numerical examples, including a nonconvex polygon with several corner singularities, and a cavity domain, are solved to around 10 digits of accuracy with a few seconds of CPU time. The examples are implemented concisely with MPSpack, a MATLAB toolbox for wave computations with nonpolynomial basis functions, developed by the authors. A code example is included.

First order k-th moment finite element analysis of nonlinear operator equations with stochastic data

We develop and analyze a class of efficient Galerkin approximation methods for uncertainty quantification of nonlinear operator equations. The algorithms are based on sparse Galerkin discretizations of tensorized linearizations at nominal parameters. Specifically, we consider abstract, nonlinear, parametric operator equations J(\alpha ,u)=0 for random input \alpha (\omega ) with almost sure realizations in a neighborhood of a nominal input parameter \alpha _0. Under some structural assumptions on the parameter dependence, we prove existence and uniqueness of a random solution, u(\omega ) = S(\alpha (\omega )). We derive a multilinear, tensorized operator equation for the deterministic computation of k-th order statistical moments of the random solution's fluctuations u(\omega ) - S(\alpha _0). We introduce and analyse sparse tensor Galerkin discretization schemes for the efficient, deterministic computation of the k-th statistical moment equation. We prove a shift theorem for the k-point correlation equation in anisotropic smoothness scales and deduce that sparse tensor Galerkin discretizations of this equation converge in accuracy vs. complexity which equals, up to logarithmic terms, that of the Galerkin discretization of a single instance of the mean field problem. We illustrate the abstract theory for nonstationary diffusion problems in random domains.

An improved finite element space for discontinuous pressures

We consider incompressible Stokes flow with an internal interface at which the pressure is discontinuous, as happens for example in problems involving surface tension. We assume that the mesh does not follow the interface, which makes classical interpolation spaces to yield suboptimal convergence rates (typically, the interpolation error in the L(2)(Omega)-norm is of order h(1/2)). We propose a modification of the P(1)-conforming space that accommodates discontinuities at the interface without introducing additional degrees of freedom or modifying the sparsity pattern of the linear system. The unknowns are the pressure values at the vertices of the mesh and the basis functions are computed locally at each element, so that the implementation of the proposed space into existing codes is straightforward. With this modification, numerical tests show that the interpolation order improves to O(h(3/2)). The new pressure space is implemented for the stable P(1)(+)/P(1) mini-element discretization, and for the stabilized equal-order P(1)/P(1) discretization. Assessment is carried out for Poiseuille flow with a forcing surface and for a static bubble. In all cases the proposed pressure space leads to improved convergence orders and to more accurate results than the standard P(1) space. In addition, two Navier-Stokes simulations with moving interfaces (Rayleigh-Taylor instability and merging bubbles) are reported to show that the proposed space is robust enough to carry out realistic simulations. (c) 2009 Elsevier B.V. All rights reserved.

Numeric simulation of pollutant dispersion by a control-volume based on finite element method

The dispersion of pollutants in the environment is an issue of great interest as it directly affects air quality, mainly in large cities. Experimental and numerical tools have been used to predict the behavior of pollutant species dispersion in the atmosphere. A software has been developed based on the control-volume based on the finite element method in order to obtain two-dimensional simulations of Navier-Stokes equations and heat or mass transportation in regions with obstacles, varying position of the pollutant source. Numeric results of some applications were obtained and, whenever possible, compared with literature results showing satisfactory accordance. Copyright (C) 2010 John Wiley & Sons, Ltd.

Influence of Cusp Inclination on Stress Distribution in Implant-Supported Prostheses. A Three-Dimensional Finite Element Analysis

Purpose: The aim of this study was to assess the influence of cusp inclination on stress distribution in implant-supported prostheses by 3D finite element method.Materials and Methods: Three-dimensional models were created to simulate a mandibular bone section with an implant (3.75 mm diameter x 10 mm length) and crown by means of a 3D scanner and 3D CAD software. A screw-retained single crown was simulated using three cusp inclinations (10 degrees, 20 degrees, 30 degrees). The 3D models (model 10d, model 20d, and model 30d) were transferred to the finite element program NeiNastran 9.0 to generate a mesh and perform the stress analysis. An oblique load of 200 N was applied on the internal vestibular face of the metal ceramic crown.Results: The results were visualized by means of von Mises stress maps. Maximum stress concentration was located at the point of application. The implant showed higher stress values in model 30d (160.68 MPa). Cortical bone showed higher stress values in model 10d (28.23 MPa).Conclusion: Stresses on the implant and implant/abutment interface increased with increasing cusp inclination, and stresses on the cortical bone decreased with increasing cusp inclination.

Platform Switching: Biomechanical Evaluation Using Three-Dimensional Finite Element Analysis

Purpose: The objective of this study was to evaluate, using three-dimensional finite element analysis (3D FEA), the stress distribution in peri-implant bone tissue, implants, and prosthetic components of implant-supported single crowns with the use of the platform-switching concept. Materials and Methods: Three 3D finite element models were created to replicate an external-hexagonal implant system with peri-implant bone tissue in which three different implant-abutment configurations were represented. In the regular platform (RP) group, a regular 4.1-mm-diameter abutment (UCLA) was connected to regular 4.1-mm-diameter implant. The platform-switching (PS) group was simulated by the connection of a wide implant (5.0 mm diameter) to a regular 4.1-mm-diameter UCLA abutment. In the wide-platform (WP) group, a 5.0-mm-diameter UCLA abutment was connected to a 5.0-mm-diameter implant. An occlusal load of 100 N was applied either axially or obliquely on the models using ANSYS software. Results: Both the increase in implant diameter and the use of platform switching played roles in stress reduction. The PS group presented lower stress values than the RP and WP groups for bone and implant. In the peri-implant area, cortical bone exhibited a higher stress concentration than the trabecular bone in all models and both loading situations. Under oblique loading, higher intensity and greater distribution of stress were observed than under axial loading. Platform switching reduced von Mises (17.5% and 9.3% for axial and oblique loads, respectively), minimum (compressive) (19.4% for axial load and 21.9% for oblique load), and maximum (tensile) principal stress values (46.6% for axial load and 26.7% for oblique load) in the peri-implant bone tissue. Conclusion: Platform switching led to improved biomechanical stress distribution in peri-implant bone tissue. Oblique loads resulted in higher stress concentrations than axial loads for all models. Wide-diameter implants had a large influence in reducing stress values in the implant system. INT J ORAL MAXILLOFAC IMPLANTS 2011;26:482-491

Effect of Metal-Ceramic or All-Ceramic Superstructure Materials on Stress Distribution in a Single Implant-Supported Prosthesis: Three-Dimensional Finite Element Analysis

Purpose: This three-dimensional finite element analysis study evaluated the effect of different material combinations on stress distribution within metal-ceramic and all-ceramic single implant-supported prostheses. Materials and Methods: Three-dimensional finite element models reproducing a segment of the maxilla with a missing left first premolar were created. Five groups were established to represent different superstructure materials: GP, porcelain fused to gold alloy; GR, modified composite resin fused to gold alloy; TP, porcelain fused to titanium; TR, modified composite resin fused to titanium; and ZP, porcelain fused to zirconia. A 100-N vertical force was applied to the contact points of the crowns. All models were fixed in the superior region of bone tissue and in the mesial and distal faces of the maxilla section. Stress maps were generated by processing with finite element software. Results: Stress distribution and stress values of supporting bone were similar for the GP, GR, TP, and ZP models (1,574.3 MPa, 1,574.3 MPa, 1,574.3 MPa, and 1,574.2 MPa, respectively) and different for the TR model (1,838.3 MPa). The ZP model transferred less stress to the retention screw (785 MPa) than the other groups (939 MPa for GP, 961 MPa for GR, 1,010 MPa for TP, and 1,037 MPa for TR). Conclusion: The use of different materials to fabricate a superstructure for a single implant-supported prosthesis did not affect the stress distribution in the supporting bone. The retention screw received less stress when a combination of porcelain and zirconia was used. Int J Oral Maxillofac Implants 2011;26:1202-1209

Mechanical behavior of ceramic veneer in zirconia-based restorations: a 3-dimensional finite element analysis using microcomputed tomography data

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

Three-Dimensional Finite Element Analysis of Vertical and Angular Misfit in Implant-Supported Fixed Prostheses

Purpose: Three-dimensional finite element analysis was used to evaluate the effect of vertical and angular misfit in three-piece implant-supported screw-retained fixed prostheses on the biomechanical response in the peri-implant bone, implants, and prosthetic components. Materials and Methods: Four three-dimensional models were fabricated to represent a right posterior mandibular section with one implant in the region of the second premolar (2PM) and another in the region of the second molar (2M). The implants were splinted by a three-piece implant-supported metal-ceramic prosthesis and differed according to the type of misfit, as represented by four different models: Control = prosthesis with complete fit to the implants; UAM (unilateral angular misfit) = prosthesis presenting unilateral angular misfit of 100 pm in the mesial region of the 2M; UVM (unilateral vertical misfit) = prosthesis presenting unilateral vertical misfit of 100 pm in the mesial region of the 2M; and TVM (total vertical misfit) = prosthesis presenting total vertical misfit of 100 pm in the platform of the framework in the 2M. A vertical load of 400 N was distributed and applied on 12 centric points by the software Ansys, ie, a vertical load of 150 N was applied to each molar in the prosthesis and a vertical load of 100 N was applied at the 2PM. Results: The stress values and distribution in peri-implant bone tissue were similar for all groups. The models with misfit exhibited different distribution patterns and increased stress magnitude in comparison to the control. The highest stress values in group UAM were observed in the implant body and retention screw. The groups UVM and TVM exhibited high stress values in the platform of the framework and the implant hexagon, respectively. Conclusions: The three types of misfit influenced the magnitude and distribution of stresses. The influence of misfit on peri-implant bone tissue was modest. Each type of misfit increased the stress values in different regions of the system. INT J ORAL MAXILLOFAC IMPLANTS 2011;26:788-796

Biomechanical evaluation of internal and external hexagon platform switched implant-abutment connections: An in vitro laboratory and three-dimensional finite element analysis

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

Cortical Bone Stress Distribution in Mandibles with Different Configurations Restored with Prefabricated Bar-Prosthesis Protocol: A Three-Dimensional Finite-Element Analysis

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

Evaluation of Bone Insertion Level of Support Teeth in Class I Mandibular Removable Partial Denture Associated With an Osseointegrated Implant: A Study Using Finite Element Analysis

Purpose: This study evaluated the influence of distal extension removable partial denture associated with implant in cases of different bone level of abutment tooth, using 2D finite element analysis.Materials and Methods: Eight hemiarch models were simulated: model A-presenting tooth 33 and distal extension removable partial denture replacing others teeth, using distal rest connection and no bone lost; model B-similar to model A but presenting distal guide plate connection; model C-similar to model A but presenting osseointegrated implant with ERA retention system associated under prosthetic base; model D-similar to model B but presenting osseointegrated implant as described in model C; models E, F, G, and H were similar to models A, B, C, and D but presenting reduced periodontal support around tooth 33. Using ANSYS 9.0 software, the models were loaded vertically with 50 N on each cusp tip. For results, von Mises Stress Maps were plotted.Results: Maximum stress value was encountered in model G (201.023 MPa). Stress distribution was concentrated on implant and retention system. The implant/removable partial denture association decreases stress levels on alveolar mucosa for all models.Conclusions: Use of implant and ERA system decreased stress concentrations on supporting structures in all models. Use of distal guide plate decreased stress levels on abutment tooth and cortical and trabecular bone. Tooth apex of models with reduced periodontal support presented increased stress when using distal rest. (Implant Dent 2011;20:192-201)

Finite Element Stress Analysis of Edentulous Mandibles with Different Bone Types Supporting Multiple-Implant Superstructures

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