Adsorption is a retention mechanism of fluid molecules on solid surfaces and presents a wide range of applications, such as refrigeration, separation processes and fuel storage. This work describes the modelling of gas adsorption on porous media and presents an optimisation approach for the design of adsorption systems based on topology optimisation. A comprehensive formulation for the adsorption and desorption phenomena is presented, detailing the derivation of their weak forms and assembling the matrices for the implementation of the Finite Element Method (FEM). The sensitivities of objective functions to topology optimisation design variables are obtained by the adjoint method, which is presented in detail. The numerical implementation is aided by FEniCS and Dolfin Adjoint packages, resulting in a generic software for the optimisation of adsorption systems. The method is applied for the design of Adsorbed Natural Gas (ANG) tanks for on-demand consumption. Optimised tank layouts are proposed and compared with standard tanks, presenting improvement of the design objectives. / Adsorção é um mecanismo de retenção de moléculas de um fluido em superfícies sólidas e apresenta um vasto campo de aplicações, como em refrigeração, processos de separação e armazenamento de combustível. Este trabalho descreve o modelo de adsorção de gás em meios porosos e define uma abordagem de otimização de sistemas de adsorção baseada no Método de Otimização Topológica (MOT). Uma formulação abrangente para os fenômenos de adsorção e dessorção é apresentada, detalhando a dedução de seus formas fracas e montagem das matrizes para a implementação do Método dos Elementos Finitos (MEF). As sensibilidades das funções-objetivo em relação às variáveis de projeto são obtidas pelo método adjunto, que também é apresentado em detalhe. A implementação numérica é auxiliada pelos pacotes FEniCS e Dolfin Adjoint, originando uma ferramenta genérica para a otimização de sistemas de adsorção. Este método é aplicado ao projeto de tanques de Gás Natural Adsorvido (GNA) visando transporte, armazenamento e consumo. Tanques otimizados são propostos e comparados com tanques tradicionais, apresentando melhora dos objetivos de projeto.

Mancais magnéticos com controle uniaxial são mancais híbridos que combinam um mancal radial por imãs permanentes e um mancal axial eletromagnético com controle ativo. Estes mancais, a despeito das vantagens com relação a outras modalidades de mancais magnéticos, têm como principal deficiência a ausência de capacidade de amortecimento de vibrações radiais do rotor. De modo a eliminar tal deficiência, esta tese propõe um novo mancal magnético híbrido que conjuga um mancal com controle uniaxial a um mancal eletrodinâmico. Esta tese tem como objetivo investigar e demonstrar a eficácia desta solução. O objetivo final é apresentar um mancal magnético de elevado desempenho, de arquitetura simples e robusta, com possibilidades de aplicação em máquinas rotativas de alta velocidade. O mancal aqui proposto tem como base um mancal de controle uniaxial cuja arquitetura é definida de modo a apresentar elevada rigidez radial, elevada folga mecânica na porção rotativa, reduzida inércia no rotor, assim como simplicidade na construção e funcionamento. Para tanto, aqui se emprega mancais radiais com imãs disposto em camadas, operando em repulsão. Na direção axial, são usados atuadores eletromagnéticos compatíveis com a elevada rigidez gerada pelos mancais radiais. A estes mancais são associados mancais eletrodinâmicos homopolares de dois tipos, de fluxo magnético axial e de fluxo radial. Definido conceitualmente o novo mancal, protótipos foram desenvolvidos e foram realizados ensaios de validação do conceito do novo mancal, assim como de verificação do desempenho do protótipo. Com relação ao mancal proposto, foi verificada uma levitação estável do rotor e, obtida uma rigidez radial de 44 N/mm mediante uma folga de 2,5 mm. São valores que superam quaisquer resultados encontrados em literatura sobre mancais similares. Em seguida, demonstrou-se, ao menos até uma rotação de 35 Hz (2.100 rpm), a eficácia do mancal eletrodinâmicos de fluxo axial no amortecimento de movimentos assíncronos do rotor. / Single axis controlled magnetic bearings are hybrid bearings that combine a radial bearing based on permanent magnets and an electromagnetic bearing with active control. Such bearings, despite advantages with respect to other modalities of magnetic bearings, present an important problem of absence of vibration damping capability in the radial directions of the rotor. In order to solve this problem, this thesis proposes a new hybrid magnetic bearing that conjugates a single axis controlled bearing with an electrodynamic bearing. This thesis aims investigate and demonstrate the efficiency of the proposed solution. The final goal is to present a magnetic bearing with high capabilities, of simple architecture and robust, with possibilities for applications in high speed rotary machines. The presented bearing has, as the base, a single axis controlled bearing with an architecture suitable for achieving high radial stiffness, large gap in the rotary portion, low rotor inertia, as well as simplicity in its construction and operation. For this, the bearing employs radial bearings based on magnets arranged in layers operating in repulsion mode. In the axial direction, it uses electromagnetic actuators compatible with the high radial stiffness generated by radial bearings. These bearings are associated with electrodynamic homopolar bearings of two types: one of axial magnetic flux and other with radial flux. Once, the new bearing is defined conceptually, tests are conducted to validate the principle of the new bearing, as well, tests to evaluation of the efficiency of the prototype. With respect to the proposed bearing, a stable levitation of the rotor is obtained. Also, a radial stiffness of 44 N/mm is achieved with a gap of 2.5 mm. These values are higher than any results presented in literatures related to similar bearings. Finally, it is demonstrated that at rotations speeds of until 35 Hz (2,100 rpm), the electrodynamic of axial flux is efficient to attenuate asynchronous motions of the rotor.

O objetivo deste trabalho é estudar o fenômeno da força de radiação acústica produzida por ondas progressivas e estacionárias. Neste trabalho o estudo da força produzida por ondas estacionárias é aplicado na análise de um levitador acústico e o estudo da força de radiação acústica por ondas progressivas é feito visando a futura construção de um separador acústico. Neste trabalho é utilizado o método dos elementos finitos para simular o comportamento de um levitador acústico. Primeiramente, é feita a simulação de um levitador acústico que consiste de um transdutor de Langevin com uma face de emissão plana que opera na freqüência de aproximadamente 20 kHz e um refletor plano. O método dos elementos finitos é utilizado para determinar o deslocamento da face do transdutor e o potencial acústico que atua numa esfera pequena. O deslocamento da face do transdutor obtido numericamente é comparado com o medido experimentalmente por um vibrômetro de fibra ótica e o potencial acústico determinado pelo método dos elementos é verificado experimentalmente colocando pequenas esferas de isopor no levitador. Depois de verificar o modelo numérico, o método dos elementos finitos é utilizado na otimização de um levitador acústico composto de um refletor côncavo e um transdutor com face de emissão côncava. Os resultados numéricos mostram que a força de radiação acústica no novo levitador é aumentada em 604 vezes quando comparada com o levitador composto de um transdutor com face plana e refletor plano. Este trabalho também apresenta um modelo numérico para determinar a trajetória de partículas esféricas na presença de uma onda de ultra-som progressiva. O modelo assume que as seguintes forças atuam na partícula: gravidade, empuxo, forças viscosas e força de radiação acústica devido a uma onda progressiva. Com o objetivo de não restringir o tamanho das partículas que podem ser utilizadas no modelo é empregada uma equação empírica do coeficiente de arrasto, válida para uma grande faixa de número de Reynolds. O modelo proposto requer a distribuição de pressão gerada pelo transdutor de ultra-som. A distribuição de pressão é medida experimentalmente utilizando um hidrofone calibrado. A verificação do modelo é feita soltando-se pequenas esferas de vidro (com diâmetros da ordem de 500 m) em frente a um transdutor de ultra-som de 1 MHz e 35 mm de diâmetro. / The objective of this work is to study the acoustic radiation force produced by progressive and standing waves. In this work, the studies related to the acoustic radiation force generated by ultrasonic standing waves are applied in the analysis of an acoustic levitator and the studies involving the acoustic radiation force generated by progressive waves are conducted aiming the design of acoustic separators. In this work, the finite element method is used to simulate an acoustic levitator. First, an acoustic levitator consisting of a 20 kHz Langevin ultrasonic transducer with a plane radiating surface and a plane reflector is simulated by the finite element method. The finite element method is used to determine the transducer face displacement and the acoustic radiation potential that acts on a small sphere. The numerical displacement is compared with that obtained by a fiber-optic vibration sensor and the acoustic radiation potential determined by the finite element method is verified experimentally by placing small Styrofoam spheres in the levitator. After verifying the numerical method, the finite element method was used to optimize an acoustic levitator consisting of a concave-faced transducer and a curved reflector. The numerical results show that the acoustic radiation force in the new levitator is enhanced 604 times compared with the levitator consisting of a plane transducer and a plane reflector. This work also presents a numerical model to determine the trajectory of sphere particles when submitted to ultrasonic progressive waves. This model assumes that the following forces act on the particle: gravity, buoyancy, viscous forces and acoustic radiation force due to the progressive wave. In order not to restrict the model to a small particle size range, the viscous forces that act on the sphere are modeled by an empirical relationship of drag coefficient that is valid for a wide range of Reynolds numbers. The numerical model requires the pressure field radiated by the ultrasonic transducer. The pressure field is obtained experimentally by using a calibrated needle hydrophone. The numerical model validation is done by dropping small glass spheres (on the order of 500 m diameter) in front of a 1-MHz 35-mm diameter ultrasonic transducer.

Computer simulations have become a great tool to assist the medical field. The present work is a study of a non-invasive patient-specific methodology to evaluate the hemodynamic importance of coronary stenosis using computer simulations. The severity of the lesion is evaluated by the Fractional Flow Reserve, calculated by the pressure gradient before and after the lesion. The geometry models are obtained from medical images of Computed Tomography Angiography exams, and the simulations considers the pulsatile flow and the blood as a non-Newtonian fluid. The governing equations of the blood flow are solved by using Finite Element Method applied to the numerical method called Incremental Pressure Correction Scheme, and with the aid of the libraries from the open-source software FEniCS. Computer simulations of three different patients are performed and the results are compared with the invasive FFR measurements. The methodology proposed shows to be feasible in the study and analysis of stenosed coronaries. / Simulações computacionais tornaram-se uma excelente ferramenta para auxiliar a área médica. O presente trabalho é um estudo de uma metodologia não invasiva para avaliar a importância hemodinâmica de artérias coronárias com estenose através do uso de simulações computacionais. A severidade da lesão da artéria é avaliada através do Fractional Flow Reserve, calculado pelo gradiente de pressão antes e depois da lesão. Os modelos geométricos computacionais foram obtidos a partir de imagens médicas de exames de Angiografia por Tomografia Computadorizada e as simulações consideram o fluxo pulsátil e as propriedades não-Newtonianas do sangue. As equações governantes do fluxo de sangue são resolvidas utilizando o Método dos Elementos Finitos aplicado ao método numérico chamado Incremental Pressure Correction Scheme, e com o uso de bibliotecas do programa em código aberto FEniCS. Foram realizados simulações de três pacientes e os resultados foram confrontados com as medidas invasivas do FFR. A metodologia proposta mostrou-se viável para o estudo e análise de coronárias com estenose.

Um material piezocompósito é resultante da combinação de um material piezelétrico com outros materiais não-piezelétricos, oferecendo vantagens substanciais em relação aos materiais piezelétricos convencionais. Diferentes propriedades efetivas podem ser obtidas alterando-se a fração de volume dos constituintes ou a própria topologia da célula unitária do piezocompósito. Materiais com Gradação Funcional (MGF) são materiais compósitos avançados, projetados de tal forma que sua composição varie gradualmente numa direção espacial. A vantagem do conceito MGF é não apresentar interface convencional entre os materiais da inclusão e da matriz, reduzindo assim um problema comum em materiais compósitos laminados, como por exemplo, o surgimento de concentração de tensões mecânicas. O Método de Otimização Topológica (MOT) é uma técnica computacional utilizada para se determinar a distribuição de materiais em uma estrutura ou material de forma sistemática, a fim de se extremizar uma determinada função objetivo. Assim, esse trabalho propõe uma metodologia sistemática e genérica para o projeto de materiais piezocompósitos com gradação funcional (MPGF) utilizando o MOT, tanto para aplicações quasi-estáticas, quanto para aplicações dinâmicas. Dessa forma, divide-se o projeto de materiais piezocompósitos em três grupos. O primeiro grupo consiste em um método de projeto de materiais baseado na combinação do método de homogeneização com o MOT para o projeto de MPGF para aplicações quasi-estáticas, onde o objetivo é projetar materiais piezocompósitos que, de modo geral, maximizem a conversão de energia mecânica em elétrica. A aplicação utilizada como exemplo neste trabalho são materiais empregados em dispositivos de coleta de energia. O segundo grupo visa aplicações dinâmicas de materiais piezocompósitos fonônicos, onde a propriedade de interesse é a possibilidade de se ter faixas de frequência, mais conhecidas por band gaps, nas quais ondas elásticas não se propagam. Assim, neste estudo visa-se o projeto de MPGF fonônicos com largura e posição de band gaps prescritos, empregando estruturas unidimensionais, e a maximização de diversos band gaps, empregando estruturas bidimensionais. O terceiro grupo explora o conceito de gradação geométrica, baseado em repetições de padrão ao longo do domínio de projeto, porém cada repetição tem um ou mais comprimentos modificados, de forma gradual. Dessa forma, suas propriedades alteram-se progressivamente ao longo da estrutura, embora a distribuição de materiais seja discreta, contornando assim possíveis dificuldades de manufatura. Esta abordagem é empregada visando à aplicação na coleta de energia, onde se procura maximizar a potência elétrica gerada em um resistor acoplado aos eletrodos, através da obtenção da topologia otimizada de estruturas piezocompósitas. Exemplos numéricos são apresentados de forma a ilustrar as metodologias de projeto propostas, bem como, analisar a influência dos parâmetros de otimização nos resultados. / Piezocomposite materials result from the combination of a piezoelectric material with other non-piezoelectric materials, offering advantages over conventional piezoelectric materials. Different effective properties can be obtained by changing the volume fraction of constituent materials, the shape of inclusions, or even the topology of the unit cell. Functionally Graded Materials (FGM) are composite materials, which are designed so that its composition varies gradually in space. One of the advantages of FGMs is that there is no conventional interface between the constituent materials, which reduces, for instance, microscopic stress concentration problems in composite materials. Topology Optimization Method (TOM) is a computational technique used to determine the material distribution of a structure or material in a systematic way, in order to maximize a determined objective function. Thus, this study proposes a generic and systematic methodology to design Functionally Graded Piezocomposites Materials (FGPM) using TOM, for quasi-static and dynamic applications. The study is divided into three groups. The first group combines the homogenization method with TOM in order to design FGPM for quasi-static applications, where the goal is to maximize the conversion of mechanical energy into electrical energy. The application used as an example in this study focuses materials used in energy harvesting devices. The second group focuses on dynamic applications of phononic piezocomposite materials, where the property of interest is the possibility of having frequency band gaps, in which elastic waves do not propagate. This study aims to design phononic FGPM with prescribed band gap width using one-dimensional model, and to design phononic FGPM with maximized band gaps using two-dimensional model. The third group investigates the pattern gradation concept, based on pattern repetitions over the design domain, but each pattern has one or more dimensions gradually modified. Thus, properties change gradually along the structure, although the material distribution keeps in the discrete form, thereby circumventing potential manufacturing difficulties. The objective function consists of maximizing the electric power generated in a load resistor. A projection scheme is employed to compute the element densities from design variables and control the length scale of the material density. Numerical examples are presented and discussed using the proposed methods.

Bombas de fluxo são dispositivos importantes em áreas como a Bioengenharia, Medicina, Farmácia, entre outras aplicações clássicas de Engenharia. Princípios para o bombeamento de fluidos baseados em atuadores piezelétricos estão sendo estudados no Departamento de Engenharia Mecatrônica e de Sistemas Mecânicos da Escola Politécnica, que permitem a construção de bombas de fluxo de pequena escala, ou seja, bombas de fluxo de pequena potência para deslocamento de pequenos volumes de fluido com baixo consumo de energia. O presente trabalho estuda bombas de fluxo piezelétricas de diafragma do tipo valve-less para geração de vazão. A bomba de fluxo piezelétrica de diafragma utiliza cerâmica piezelétrica como atuador para mover uma membrana (diafragma) para cima e para baixo como um pistão, que causa uma seqüência de aumento e diminuição do volume da câmara da bomba, forçando a entrada e a saída do fluido na bomba. A direção do fluxo é garantida por válvulas que privilegiam o fluxo em apenas um sentido. O objetivo deste trabalho é o estudo da metodologia de desenvolvimento de uma bomba de fluxo piezelétrica de diafragma de baixo custo do tipo valve-less. Para tanto, será utilizado a modelagem por Método dos Elementos Finitos (MEF) para a realização de análises de sensibilidade dos parâmetros geométricos e construtivos da bomba de fluxo. Serão realizadas simulações de escoamento de fluido pelo Método de Volumes Finitos (MVF) para a realização de análises de sensibilidade dos parâmetros geométricos dos elementos difusor/bocal e o levantamento das curvas características da bomba de fluxo. Por fim, protótipos serão construídos e caracterizados para validação dos resultados computacionais. Serão apresentadas a metodologia empregada e a discussão dos resultados obtidos, de forma a analisar o princípio proposto e os fenômenos físicos em questão. / Flow pumps act as important devices in areas as Bioengineering, Medicine, Pharmacy, among other areas of Engineering. Principles for pumping fluids based on piezoelectric actuators have been studied in the Department of Mechatronic and Mechanical Engineering of Escola Politécnica da Universidade de São Paulo, that allow the construction of small flow pumps, in other words, pumps for displacement of small fluid volumes with low power consumption. The present work studies valveless piezoelectric diaphragm flow pumps for flow generation. The piezoelectric diaphragm flow pump uses a piezoelectric ceramic as actuator to move a membrane (diaphragm) up and down as a piston. Consequently, there is a sequence of increase and decrease in the chamber volume that will force the fluid in and out of the pump. The direction of the flow is guaranteed by valves that privilege the flow in just one pumping direction. The main objective of this work is the study of a methodology to develop a low cost valve-less piezoelectric diaphragm flow pump. A sensitivity analysis is carried using computational simulation through the Finite Element Method (FEM) to study how construction parameters and assembly affect diaphragm flow pump performance. Using CFD simulations through the Finite Volume Method (FVM), a sensitivity analysis is done around nozzle/diffuser element geometrical parameters and pump characteristic curves are obtained. Finally, computational results are validated by prototype construction and characterization. The text presents methodologies employed and discusses the obtained results, analyzing the principle and the related physical phenomena.

O objetivo deste trabalho é estudar sistemas de coxinização, também conhecidos como sistemas de isolamento, desenvolvidos para motores automotivos juntamente com os modelos analíticos para calcular o comportamento dinâmico do motor suspenso, complementados por análises numéricas com o emprego do método dos elementos finitos e estudos experimentais. Estudos para verificar o efeito dos parâmetros de projeto do sistema de coxinização na resposta modal do motor isolado foram realizados utilizando os modelos analíticos desenvolvidos. Foi possível identificar como a variação de tais parâmetros interfere na resposta modal do motor suspenso. Os resultados fornecidos por tais modelos analíticos auxiliaram na compreensão do comportamento modal do sistema, mostrando que são ferramentas úteis para auxiliar o projeto de sistemas de isolamento de motores, especialmente em estágios iniciais do projeto, onde os parâmetros do veículo ainda estão por ser determinados e protótipos ainda não estão disponíveis. / The aim of this work is to study engine mount systems (isolation systems) developed for automotive engines along with the analytical models to calculate the dynamic behavior of the isolated engine, complemented by numerical analysis using the finite element method and experimental analysis. Studies to verify the effect that the isolation system design parameters have on the modal response of the isolated engine were performed using the analytical models developed. It was possible to identify how the variation of such parameters affects the modal response of the isolated engine. The results provided by such analytical models aided in understanding the modal behavior of the system, showing that they are useful tools to assist in the design of engine isolation systems, especially in the early stages of design, where vehicle parameters are still to be determined and prototypes are not yet available.

As ferrovias, muito importantes para o transporte de pessoas e mercadorias, são também grandes fontes de ruído e vibração no meio urbano e rural. Os esforços dinâmicos devido ao contato aço-aço, entre as rodas e trilhos, são transmitidos à fundação da via e geram ondas vibratórias que se propagam pelo solo até atingir um receptor. Devido aos potenciais danos causados por essas excitações, modelos matemáticos capazes de prever e realizar a propagação das vibrações são de grande importância, podendo nortear ações mitigatórias. Neste contexto, este trabalho propõe a criação de um modelo numérico, utilizando-se de dinâmica multicorpos imerso em um ambiente de elementos finitos com não linearidade geométrica, que procure representar a interação dinâmica no sistema veículo-via e que permita analisar o comportamento das forças de contato roda-trilho. Tal modelo é utilizado com o objetivo de analisar e compreender como as forças dinâmicas de contato são influenciadas pelas características do veículo e da via. Alguns estudos sequenciais foram realizados. Iniciou-se com a avaliação do comportamento do veículo ao trafegar em uma superfície plana. Sequencialmente, considerou-se o veículo trafegando sobre modelos de vigas rígidas, representando os trilhos. Neste caso foram estudados cenários de rolamento e deslizamento. Posteriormente, o valor de uma rigidez equivalente da via foi avaliado por meio de um modelo estrutural e incorporado em um estudo para análise das forças de impacto das rodas ao passarem sobre descontinuidades em trilhos. Os resultados deste trabalho podem ser utilizados como dados de entrada para futuros modelos para previsão de propagação de ondas no solo devido às cargas ferroviárias. / Railways are important to freight and people transportation but are also one of the greatest sources of vibration in urban and rural environments. Dynamic loads are generated in steel-steel contact between wheels and rails and are transferred to foundation. They generate oscillatory waves that propagate through the ground, until finding a receptor. Potential damage due to ground borne vibration has increased the importance of predicting and calculating wave propagation with the aid of computational models, serving as guidelines for mitigating actions. On this context, the present work proposes a numerical model, adopting a multibody rigid-flexible dynamics formulation, based on a nonlinear geometrically-exact finite element description. The aim is to represent the dynamic interaction between the vehicle and the railway, permitting one to analyze wheel-rail contact forces. The model represents the dynamic interaction between vehicle and track and allows one to analyze how the dynamic contact forces are influenced by vehicle and railway characteristics. Some sequential studies were performed, starting from the evaluation of a vehicle running on a plane surface. Next, the scenario of a vehicle running on rigid beams representing the rails was simulated. In this case both sliding and sticking scenarios of contact were considered. Then, the value of an equivalent track stiffness was obtained by a structural model. It was taken as input to another study, which analyzed impact forces between wheels and rails when the vehicle experiences the passage by a discontinuity on the track. The results of the present work may be useful as input data for future prediction models of ground borne vibration due to railway dynamic forces.

Este trabalho é definido a partir da elaboração de uma metodologia de cálculo para a determinação da evolução dos parâmetros de anisotropia em chapas metálicas e na definição de uma nova superfície de escoamento anisotrópica. A metodologia proposta é fundamentada nas relações incrementais entre tensão e deformação da teoria clássica da plasticidade utilizando-se do critério de escoamento de Hill e baseado em um estado particular de tensão-deformação comum em processos de conformação mecânica, principalmente em regiões da peça onde ocorram dobras. A nova superfície de escoamento é definida a partir da superfície de Hill original, porém ponderada por um fator de correção o qual contempla uma evolução anisotrópica do encruamento. Este fator de correção é definido em função dos parâmetros evolutivos de anisotropia determinados pela metodologia proposta e um conceito particular é então adotado para a complementação deste novo modelo de superfície de escoamento. A contribuição científica apresentada estabelece o desenvolvimento do conjunto de equações que permite determinar a evolução dos parâmetros de anisotropia com o carregamento plástico e adicionalmente permite avaliar a evolução de nova superfície servindo como indicativo da formabilidade do material. Os resultados teóricos são comparados e discutidos com evoluções obtidas via ensaio normalizado de anisotropia. Além disso é realizada a aplicação do procedimento proposto em conjunto com análises numéricas via método dos elementos finitos para dois casos de processos de conformação. Desta forma tem-se uma contribuição para a representação do fenômeno de anisotropia plástica dado por um modelo evolutivo, cuja aplicação tem interesse em aprimoramento de algoritmos computacionais bem como em aplicações de cunho experimental como validações, predições e aprimoramentos de processos de conformação e outros. / This work is aimed to the elaboration of a new approach for prediction of anisotropy parameters evolution in sheet metal as well as definition of a new anisotropic yield surface. This approach is based on incremental relation between stress and strain from classical plasticity theory by using the Hill-48 yield criteria. Also, it is based on a particular stress-strain state, common to bending regions found in several sheet forming processes. The new yield surface is derived from the original Hill criteria, which was modified by a correcting factor proposed to represent an anisotropic evolution of the continuing plastic loading. This factor is a strict function of the evolutive anisotropy parameters generated by the present methodology. From this point, is adopted a special concept in order to reach a complete definition of the new model for yield surfaces. The scientific contribution presents the development of a set of equations in order to predict the evolution of the anisotropy parameters as function of plastic loading, and additionally, it is taken this new yield surface as an evaluation of material formability degree. The theoretical results are discussed and compared with results obtained from standardized anisotropy tests. Besides, the proposed procedure is applied, aided by numerical analysis using the finite element method to simulate two cases of sheet forming processes. Therefore it is presented a contribution to the representation of plastic anisotropy phenomenon using an evolutive model. The main interests for practical application of this work range from improvement of computational algorithms as well as experimental tests investigations, which are related to validation of this model, prediction and improvement of sheet metal forming processes.

Uniões por prendedores são elementos amplamente utilizados na indústria aeronáutica para a união de partes constituintes da aeronave. Contudo, devido à sua geometria e aos carregamentos sofridos, estes elementos estão frequentemente sujeitos a falhas por fadiga. Assim, para um projeto e dimensionamento bem executado dessas juntas, é necessário conhecer seu comportamento mecânico e o campo de tensões ao qual estão sujeitas. O método dos elementos finitos certamente atende a estas necessidades; porém, o uso de elementos sólidos tridimensionais para a representação destas uniões pode levar a análises demasiadamente demoradas e custosas, sendo desejável o uso de modelos mais simplificados. Nesse trabalho, juntas de topo assimétricas são modeladas pelo método dos elementos finitos, utilizando tanto elementos sólidos tridimensionais quanto elementos de casca, com o objetivo de encontrar um modelo relativamente simples que apresente resultados satisfatórios e requeira um menor tempo de solução. Os resultados numéricos obtidos são comparados com resultados experimentais, que utilizam extensômetros e fotoelasticidade. / Riveted joints are structural elements widely used in the aeronautic industry to connect different parts of the aircraft. However, due to their geometry and working load, they are frequently subjected to fatigue failure. Therefore, to a well executed project of this type of joint, it is necessary to understand the mechanical behavior and the stress field to which they are subjected. The finite element method can certainly answer those needs; however, the use of tridimensional elements to represent this type of joint can lead to time consuming and expensive analysis, being desirable the use of simpler models. In this work, asymmetric butt joints are modeled using the finite element method using both, tridimensional and shell elements, with the objective of finding a model that provides satisfactory results at smaller solution time. The numerical results obtained are compared with experimental results using strain gage and photoelastic procedures.