Numerical solution of the eXtended Pom-Pom model for viscoelastic free surface flows

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

Evaluation of a bounded high order upwind scheme for 3D incompressible free surface flow computations

fit the context of normalized variable formulation (NVF) of Leonard and total variation diminishing (TVD) constraints of Harten. this paper presents an extension of it previous work by the authors for solving unsteady incompressible flow problems. The main contributions of the paper are threefold. First, it presents the results of the development and implementation of a bounded high order upwind adaptative QUICKEST scheme in the 3D robust code (Freeflow), for the numerical solution of the full incompressible Navier-Stokes equations. Second, it reports numerical simulation results for 1D hock tube problem, 2D impinging jet and 2D/3D broken clam flows. Furthermore, these results are compared with existing analytical and experimental data. and third, it presents the application of the numerical method for solving 3D free surface flow problems. (C) 2007 IMACS. Published by Elsevier B.V. All rights reserved,

Determination of the numerical parameters of a continuous damage model for the structural analysis of clay brick masonry

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

A planar flexible robotic manipulator

Presents the dynamic modelling of a flexible robotic manipulator with two flexible links and two revolute joints, which rotates in the horizontal plane. The dynamic equations are derived using the Newton-Euler formulation and the finite element method, based on elementary beam theory, which is used to discretize the displacements such that the small motion is represented in terms of nodal displacements. Computer simulation results are presented to illustrate this study. The dynamic model becomes necessary for use in future design and control applications.

HIGH STRENGTH CONCRETE: THE REINFORCED CONCRETE OF THE CENTURY 21(ST)

This work is about the 21st century reinforced concrete analysis under the point of view of its constituent materials. First of all it is described the theoretical approach of the bending elements calculated based on the Norms BAEL 91 standarts. After that, numerical load-displacement are presented from reinforced concrete beams and plates validated by experimental data. The numerical modellings has been carried on in the program CASTEM 2000. In this program a elastoplastic model of Drucker-Prager defines the rupture surface of the concrete in non associative plasticity. The crack is smeared on the Gauss points of the finite elements with formation criterion starting from the definition of the rupture surface in the branch traction-traction of the Rankine model. The reinforcements were modeled in a discrete approach with perfect bond. Finally, a comparative analysis is made between the numerical results and calculated criteria showing the future of high performance reinforced concrete in this beginning of 21st century.

NUMERICAL SIMULATION of THE BEHAVIOR of THE BEHAVIOR of REINFORCED CONCRETE BEAMS and COLUMNS BY FEM CAASTEM 2000 PROGRAM

The objective of this paper is the numerical study of the behavior of reinforced concrete beams and columns by non-linear numerical simulations. The numerical analysis is based on the finite element method implemented in CASTEM 2000. This program uses the constitutive elastoplastic perfect model for the steel, the Drucker-Prager model for the concrete and the Newton-Raphson for the solution of non-linear systems. This work concentrates on the determination of equilibrium curves to the beams and force-strain curves to the columns. The numeric responses are confronted with experimental results found in the literature in order to check there liability of the numerical analyses.

Levenberg-Marquardt application to two-phase nonlinear parameter estimation for finned-tube coil evaporators

A procedure for calculation of refrigerant mass flow rate is implemented in the distributed numerical model to simulate the flow in finned-tube coil dry-expansion evaporators, usually found in refrigeration and air-conditioning systems. Two-phase refrigerant flow inside the tubes is assumed to be one-dimensional, unsteady, and homogeneous. In themodel the effects of refrigerant pressure drop and the moisture condensation from the air flowing over the external surface of the tubes are considered. The results obtained are the distributions of refrigerant velocity, temperature and void fraction, tube-wall temperature, air temperature, and absolute humidity. The finite volume method is used to discretize the governing equations. Additionally, given the operation conditions and the geometric parameters, the model allows the calculation of the refrigerant mass flow rate. The value of mass flow rate is computed using the process of parameter estimation with the minimization method of Levenberg-Marquardt minimization. In order to validate the developed model, the obtained results using HFC-134a as a refrigerant are compared with available data from the literature.

Vibration Attenuation in Rotating Machines Using Smart Spring Mechanism

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

Experimental Active Vibration Control in Truss Structures Considering Uncertainties in System Parameters

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

Modelagem Matemática e Computacional de fenômenos eletrocinéticos em meios porosos carregados eletricamente

In this work we present a mathematical and computational modeling of electrokinetic phenomena in electrically charged porous medium. We consider the porous medium composed of three different scales (nanoscopic, microscopic and macroscopic). On the microscopic scale the domain is composed by a porous matrix and a solid phase. The pores are filled with an aqueous phase consisting of ionic solutes fully diluted, and the solid matrix consists of electrically charged particles. Initially we present the mathematical model that governs the electrical double layer in order to quantify the electric potential, electric charge density, ion adsorption and chemical adsorption in nanoscopic scale. Then, we derive the microscopic model, where the adsorption of ions due to the electric double layer and the reactions of protonation/ deprotanaç~ao and zeta potential obtained in modeling nanoscopic arise in microscopic scale through interface conditions in the problem of Stokes and Nerst-Planck equations respectively governing the movement of the aqueous solution and transport of ions. We developed the process of upscaling the problem nano/microscopic using the homogenization technique of periodic structures by deducing the macroscopic model with their respectives cell problems for effective parameters of the macroscopic equations. Considering a clayey porous medium consisting of kaolinite clay plates distributed parallel, we rewrite the macroscopic model in a one-dimensional version. Finally, using a sequential algorithm, we discretize the macroscopic model via the finite element method, along with the interactive method of Picard for the nonlinear terms. Numerical simulations on transient regime with variable pH in one-dimensional case are obtained, aiming computational modeling of the electroremediation process of clay soils contaminated

Influence of Ridge Inclination and Implant Localization on the Association of Mandibular Kennedy Class I Removable Partial Denture

The aim of this study was to evaluate the tendency of displacement of the supporting structures of the distal extension removable partial denture (DERPD) associated to the implant with different inclinations of alveolar ridge and implant localizations through a two-dimensional finite-element method. Sixteen mandibular models were fabricated, presenting horizontal, distally descending, distally ascending, or descending-ascending ridges. All models presented the left canine and were rehabilitated with conventional DERPD or implant-retained prosthesis with the ERA system. The models were obtained by the AutoCAD software and transferred to the finite-element software ANSYS 9.0 for analysis. A force of 50 N was applied on the cusp tips of the teeth, with 5 points of loading of 10 N. The results were visualized by displacement maps. For all ridge inclinations, the assembly of the DERPD with distal plate retained by an anterior implant exhibited the lowest requisition of the supporting structures. The highest tendency of displacement occurred in the model with distally ascending ridge with incisal rest. It was concluded that the association of the implant decreased the displacement of the DERPD, and the anterior positioning of the implant associated to the DERPD with the distal plate preserved the supporting structures for all ridges.

Influence of customized composite resin fibreglass posts on the mechanics of restored treated teeth

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

Influence of Implant Inclination Associated With Mandibular Class I Removable Partial Denture

The aim of this study was to use two-dimensional finite element method to evaluate the displacement and stress distribution transmitted by a distal extension removable partial denture (DERPD) associated with an implant placed at different inclinations (0, 5, 15, and 30 degrees) in the second molar region of the edentulous mandible ridge. Six hemimandibular models were created: model A, only with the presence of the natural tooth 33; model B, similar to model A, with the presence of a conventional DERPD replacing the missing teeth; model C, similar to the previous model, with a straight implant (0 degrees) in the distal region of the ridge, under the denture base; model D, similar to model C, with the implant angled at 5 degrees in the mesial direction; model E, similar to model C, with the implant angled at 15 degrees in the mesial direction; and model F, similar to ME, with the implant angled at 30 degrees in the mesial direction. The models were created with the use of the AutoCAD 2000 program (Autodesk, Inc, San Rafael, CA) and processed for finite element analysis by the ANSYS 8.0 program (Swanson Analysis Systems, Houston, PA). The force applied was vertical of 50 N on each cusp tip. The results showed that the introduction of the RPD overloaded the supporting structures of the RPD and that the introduction of the implant helped to relieve the stresses of the mucosa alveolar, cortical bone, and trabecular bone. The best stress distribution occurred in model D with the implant angled at 5 degrees. The use of an implant as a support decreased the displacement of alveolar mucosa for all inclinations simulated. The stress distribution transmitted by the DERPD to the supporting structures was improved by the use of straight or slightly inclined implants. According to the displacement analysis and von Mises stress, it could be expected that straight or slightly inclined implants do not represent biomechanical risks to use.

Evaluation of the influence of location of osseointegrated implants associated with mandibular removable partial dentures

Purpose: The aim of this research was to assess, by means of, the bi-dimensional finite element method, the best implant location in the alveolar edge, through stress distribution and support structure displacement of a distal extension removable partial denture associated with an osseointegrated implant of 10.0 x .75 mm, acting as abutment for the denture base.Methods and Materials: Five models in sagittal cut were used to represent: model A-hemi arch containing natural tooth 33 and the distal alveolar edge; model B-similar to model A, but with a conventional removable partial denture to replace the absent teeth; model C (MC)-similar to the previous one, with an implant in the distal region of the edge under the denture base; model D-similar to MC, with the implant in the central region of the edge; model E-similar to MC, with an implant in the mesial region of the edge. With the aid of the finite element program ANSYS 8.0, the models were loaded with strictly vertical forces of 50 N on each cusp tip. Displacement and von Mises Maps were plotted for visualization of results.Results: The introduction of implant diminished the tendency of intrusion of the removable partial denture in all situations. The maximum stress was observed on implant in all situations. Approximating implant in direction of support teeth was benefit for stress distribution.Conclusion: Model D presented the lowest value for maximum tendency to displacement when compared with those found in the other models; model E demonstrated better relief with regard to demand from the abutment tooth; locating the implant near of the abutment tooth influenced positively the distribution of stresses on the analyzed structures.

Evaluation of Different Retention Systems on a Distal Extension Removable Partial Denture Associated With an Osseointegrated Implant

The aim of this study was to evaluate the biomechanical behavior of a mandibular distal extension removable partial denture (DERPD) associated with an implant and different retention system, by bidimensional finite element method. Five hemimandible models with a canine and external hexagon implant at second molar region associated with DERPD were simulated: model A, hemimandible with a canine and a DERPD; model B, hemimandible with a canine and implant with a healing abutment associated to a DERPD; model C, hemimandible with a canine and implant with an ERA attachment associated to a DERPD; model D, hemimandible with a canine and implant with an O'ring attachment associated to a DERPD; and model E, hemimandible with a canine and implant-supported prosthesis associated to a DERPD. Cusp tips were loaded with 50 N of axial or oblique force (45 degrees). Finite element analysis was performed in ANSYS 9.0. model E showed the higher displacement and overload in the supporting tissues; the patterns of stress distribution around the dental apex of models B, C, and D were similar. The association between a DERPD and an osseointegrated implant using the ERA or O'ring systems shows lower stress values. Oblique forces showed higher stress values and displacement. Oblique forces increased the displacement and stress levels in all models; model C displayed the best stress distribution in the supporting structures; healing abutment, ERA, and O'ring systems were viable with RPD, but DERPD association with a single implant-supported prosthesis was nonviable.