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.

Methods Used for Assessing Stresses in Buccomaxillary Prostheses: Photoelasticity, Finite Element Technique, and Extensometry

The authors describe a literature revision on assessing stresses in buccomaxillary prostheses photoelasticity, finite element technique, and extensometry. They describe the techniques and the importance for use of each method in buccomaxillary prostheses with implants and the need of accomplishing more studies in this scarce literary area.

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.

Interpolating EFGM for computing continuous and discontinuous electromagnetic fields

Purpose - This paper proposes an interpolating approach of the element-free Galerkin method (EFGM) coupled with a modified truncation scheme for solving Poisson's boundary value problems in domains involving material non-homogeneities. The suitability and efficiency of the proposed implementation are evaluated for a given set of test cases of electrostatic field in domains involving different material interfaces.Design/methodology/approach - the authors combined an interpolating approximation with a modified domain truncation scheme, which avoids additional techniques for enforcing the Dirichlet boundary conditions and for dealing with material interfaces usually employed in meshfree formulations.Findings - the local electric potential and field distributions were correctly described as well as the global quantities like the total potency and resistance. Since, the treatment of the material interfaces becomes practically the same for both the finite element method (FEM) and the proposed EFGM, FEM-oriented programs can, thus, be easily extended to provide EFGM approximations.Research limitations/implications - the robustness of the proposed formulation became evident from the error analyses of the local and global variables, including in the case of high-material discontinuity.Practical implications - the proposed approach has shown to be as robust as linear FEM. Thus, it becomes an attractive alternative, also because it avoids the use of additional techniques to deal with boundary/interface conditions commonly employed in meshfree formulations.Originality/value - This paper reintroduces the domain truncation in the EFGM context, but by using a set of interpolating shape functions the authors avoided the use of Lagrange multipliers as well Mathematics in Engineering high-material discontinuity.

Automation and cybernetics: Control of a flexible one-link manipulator

States that control is of the essence in cybernetics. Summarizes the dynamic equations for a flexible one-link manipulator moving in the horizontal plane. Employs the finite element method, based on elementary beam theory, during the process of formulation. Develops and instruments a one-link flexible manipulator in order to control its vibration modes. Uses a simple second-order vibration model which permits vibrations on the rod to be estimated using the hub angle. The validation of the dynamic model and the structural analysis of the flexible manipulator is reached using proper infrared cameras and active light sources for determining actual positions of objects in space. Shows that the performance of the control is satisfactory, even under perturbation action.

Application of a model updating technique for health condition monitoring and damage detection of flexible structures

This paper discusses the application of a damage detection methodology to monitor the location and extent of partial structural damage. The methodology combines, in an iterative way, the model updating technique based on frequency response functions (FRF) with monitoring data aiming at identifying the damage area of the structure. After the updating procedure reaches a good correlation between the models, it compares the parameters of the damage structure with those of the undamaged one to find the deteriorated area. The influence of the FEM mesh size on the evaluation of the extent of the damage has also been discussed. The methodology is applied using real experimental data from a spatial frame structure.

Optimal design of smart structures using bonded piezoelectrics for vibration control

Smart material technology has become an area of increasing interest for the development of lighter and stronger structures which are able to incorporate actuator and sensor capabilities for collocated control. In the design of actively controlled structures, the determination of the actuator locations and the controller gains, is a very important issue. For that purpose, smart material modelling, modal analysis methods, control and optimization techniques are the most important ingredients to be taken into account. The optimization problem to be solved in this context presents two interdependent aspects. The first one is related to the discrete optimal actuator location selection problem, which is solved in this paper using genetic algorithms. The second is represented by a continuous variable optimization problem, through which the control gains are determined using classical techniques. A cantilever Euler-Bernoulli beam is used to illustrate the presented methodology.

Influence of convergence angle of tooth preparation on the fracture resistance of Y-TZP-based all-ceramic restorations

Objective: To investigate the influence of the convergence angle of tooth preparation on the fracture load of Y-TZP-based ceramic (YZ-Vita YZ) substructure (SB) veneered with a feldspathic porcelain (VM9-Vita VM9). Methods: Finite element stress analysis (FEA) was performed to examine the stress distribution of the system. Eighty YZ SB were fabricated using a CAD-CAM system and divided into four groups (n = 20), according to the total occlusal convergence (TOC) angle: G6-6° TOC; G12-12° TOC; G20-20° TOC; and G20MOD-20° TOC with modified SB. All SB were veneered with VM9, cemented in a fiber reinforced epoxy resin die, and loaded to failure. Half of the specimens from each group (n = 10) were cyclic fatigued (106 cycles) before testing. Failure analysis was performed to determine the fracture origin. Data were statistically analyzed using Anova and Tukey's tests (α = 0.05). Results: The greatest mean load to fracture value was found for the G20MOD, which was predicted by the FEA. Cyclic fatigue did not significantly affect the load of fracture. Catastrophic failure originating from the internal occlusal surface of the SB was the predominant failure mode, except for G20MOD. Significance: The YZ-VM9 restorations resisted greater compression load than the usual physiological occlusal load, regardless of the TOC angle of preparations. Yet, the G20MOD design produced the best performance among the experimental conditions evaluated. © 2013 Academy of Dental Materials. Published by Elsevier Ltd. All rights reserved.

Strain rate of the South American lithospheric plate by SIRGAS-CON geodetic observations

A multiyear solution of the SIRGAS-CON network was used to estimate the strain rates of the earth surface from the changing directions of the velocity vectors of 140 geodetic points located in the South American plate. The strain rate was determined by the finite element method using Delaunay triangulation points that formed sub-networks; each sub-network was considered a solid and homogeneous body. The results showed that strain rates vary along the South American plate and are more significant on the western portion of the plate, as expected, since this region is close to the subduction zone of the Nazca plate beneath the South American plate. After using Euler vectors to infer Nazca plate movement and to orient the velocity vectors of the South American plate, it was possible to estimate the convergence and accommodation rates of the Nazca and South American plates, respectively. Strain rate estimates permitted determination of predominant contraction and/or extension regions and to establish that contraction regions coincide with locations with most of the high magnitude seismic events. Some areas with extension and contraction strains were found to the east within the stable South American plate, which may result from different stresses associated with different geological characteristics. These results suggest that major movements detected on the surface near the Nazca plate occur in regions with more heterogeneous geological structures and multiple rupture events. Most seismic events in the South American plate are concentrated in areas with predominant contraction strain rates oriented northeast-southwest; significant amounts of elastic strain can be accumulated on geological structures away from the plate boundary faults; and, behavior of contractions and extensions is similar to what has been found in seismological studies. © 2013 Elsevier Ltd.

Desenvolvimento de pino intra-radicular, resistência à fratura e análise da distribuição de tensões através do método dos elementos finitos

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Avaliação da resistência á flexão e determinação das propriedades mecânicas de resinas compostas e sua aplicação no método dos elementos finitos

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)