[pt] Este trabalho trata da solução numérica das equações de Navier-Stokes, na forma vorticidade-função corrente, via método das Diferenças Finitas e técnicas de aceleração baseadas no uso de malhas múltiplas. Embora outras opções tenham sido consideradas, a que melhor funcionou tratou o problema de forma não acoplada: a solução da equação de vorticidade foi obtida pela uso desta aceleração e a solução da equação de função corrente, uma equação puramente elíptica, foi resolvida via método das relaxações sucessivas. O código desenvolvido foi aplicado a diversos problemas, inclusive ao problema da cavidade com tampa móvel, em diversos números de Reynolds, típico no teste de simuladores em Dinâmica dos Fluidos. Foram testados um método clássico (armazenamento da correção) e o método FAZ (Full Approximation Storage). Os resultados obtidos mostram claramente os ganhos computacionais obtidos na formulação escolhida. Expressando em percentual, valores com 80 por cento de ganho foram obtidos se comparados os resultados do método multigrid com o método iterativo básico utilizado (S.O.R.), indicando o potencial do uso desta técnica para problemas mais complexo incluindo aqueles em coordenadas generalizadas. / [en] This works deals with the numerical solution of the Navier-Stokes equations, written in the stream function-vorticity form, by the finite difference method and acceleration techniques using multiple meshes. Although other solution schemes have been investigated, best results were obtained by treating the problem in a non-coupled form: the solution for the vorticity equation was obtained by the multigrid method and the solution of the streamfunction equation, which is purely elliptic, was solved by the S.O.R. (Successive over relaxation method). The computer code was applied to several problems, including the wall driven problem considering a wide range of Reynolds numbers, which is a typical benchmark problem for testing fluid-dynamic simulations. The classical method (storage of the correction) and the methos FAS (Full Approximation Storage) have been tested. The results obtained clearly show that a very efficient computational scheme has been achieved with the multigrid method. For example, when comparing this method with the basic S.O.R. method, relative gains in the order of 80 per cent have been obtained. This indicates that the present technique has potential use in more complicated fluid dynamics problems including those involving generalized coordinates.

[pt] Este trabalho apresenta o estudo do escoamento transiente em difusores radiais com fronteira móvel alimentados axialmente. O problema é resolvido acoplando-se a modelagem do fluido à dinâmica do movimento da palheta. Admite-se o comportamento Newtoniano para o fluido de trabalho (ar). O escoamento é considerado incompressível, laminar e isotérmico numa geometria axisimétrica. As equações de conservação são resolvidas pelo método dos volumes-finitos utilizando-se um código numérico comercial, o FLUENT, adaptado de modo a permitir que o posicionamento da fronteira móvel faça parte da solução do problema. Analise-se os casos de movimentação súbita imposta à palheta e de fluxo mássico de alimentação periódica. Resultados do comportamento temporal da força axial do fluido sobre a palheta são apresentados. Determina-se também o campo de geração de entropia para o escoamento bem como sua variação para diversos instantes de tempo. / [en] This work presents the study of unsteady flow in a axially fed radial flow between parallel disks on a situation in which the frontal disk can move. The problem is solved coupling the fluid flow modeling and the valve dynamics movement. The fluid (air) is assumed to have a Newtonian behavior. The flow is incompressible, laminar and isothermal in an axisymmetric geometry. The governing equations are solved by finite volume methodology using a commercial software, FLUENT, which have been modified to admit boundary movement as part of the solution. The flow characteristics presented herein relate to situations where: I. the frontal disk is suddenly moved to a new location and ii. prescribed mass flux at the inlet orifice. Results for the behavior of the axial force in the frontal disk can be seen. Another important parameter analyzed is the entropy generation and its time variation.

[pt] Este trabalho apresenta uma análise dinâmica de uma turbina eólica com torre tubular e com rotor de três pás sob carregamentos determinísticos e estocásticos, concentrados e distribuídos, dentre eles os provocados pelo vento e por correntes marinhas. Variações nas condições do solo são permitidas sendo simuladas como molas axiais e torcionais lineares instaladas na base da torre. Os carregamentos aleatórios são devidos à ação do vento sobre as pás da turbina. As simetrias envolvidas permitem que a torre seja modelada como quatro elementos estruturais distintos e superpostos: uma viga com movimento espacial (dois movimentos laterais independentes), um eixo (movimento de torção) e uma barra (movimento axial). Elementos finitos são utilizados para a aproximação da dinâmica da torre, sendo então construído um modelo reduzido, utilizando uma base de dimensão finita com autovetores da estrutura. Os demais componentes da turbina, nacele, eixo e pás são modelados como corpos rígidos. Todos os componentes são agrupados e para modelar seus vínculos são impostas equações de restrição. Um algoritmo baseado no método de Newmark foi desenvolvido para integrar o sistema de equações no tempo. Velocidades, acelerações, forças e torques são obtidos para diversas configurações. / [en] This thesis presents a dynamic analysis of a three-blade wind turbine with tubular cross section under deterministic and stochastic concentrated and distributed loads including those caused by wind and by sea currents. Variations in soil conditions are allowed, being modeled as linear axial and torsional springs installed at the base of the tower. The random loading is due to the action of the wind on the turbine blades. The symmetries involved allow the tower to be modeled as having four separated and superimposed structural behaviors: a beaming spatial (lateral motion), a shaft (torsional) motion and a bar (axial) motion. Finite elements are used for approximating the dynamic of the tower and then a reduced model is constructed using a basis of finite dimension of the structure eigenvectors. The other turbine components, namely, the nacelle, the shaft and the blades are modeled as rigid bodies whose links are imposed by constraints functions. An algorithm based on Newmark method was developed to integrate the system of equations in time. Velocities, accelerations, forces and torques as well as the resulting stresses are obtained for several configurations.

[pt] Neste trabalho estuda-se um sistema de acompanhamento de alvos, também conhecido como sistema de rastreamento de alvos, do tipo pan-tilt atuado por motores de corrente contínua e fixo em um corpo em movimento. Para tanto é montada uma bancada de testes e implementa-se um programa de simulação. A modelagem para o programa é feita no domínio do tempo, permitindo a utilização de equações bastante complexas para representar o sistema, o que não é possível quando se utiliza modelos no domínio da freqüência. Apesar de se modelar o sistema com corpos rígidos, flexibilidades e amortecimentos estruturais devido aos redutores de velocidade são considerados. Erros nos sensores, folgas nos redutores, atritos seco e viscoso, limites de saturação para as correntes e tensões nas armaduras dos motores também são considerados. Um método para a inclusão dos atrasos de tempo para atualização dos sinais de controle e dados obtidos pelos sensores durante a integração numérica das equações de movimento é apresentado. Para controlar o sistema utilizam-se controladores que não requerem o modelo matemático da planta, tanto na bancada de testes como no programa de simulação. Três tipos diferentes de arquitetura de controle são propostas, chamadas neste trabalho de tipo 1, tipo 2 e tipo 3. A complexidade delas aumenta à medida que mais sensores estão disponíveis no sistema. A arquitetura do tipo 1 destina-se a sistemas onde se possui apenas sensores que fornecem os erros angulares de azimute e elevação do alvo. Se, além deste sensor, também houver sensores para medir as posições angulares relativas entre os elos do mecanismo usa-se a arquitetura do tipo 2. Se houver, ainda, sensores de velocidades angulares inerciais pode-se utilizar a arquitetura do tipo 3. Por fim são apresentados resultados de experimentos e simulações onde se compara o desempenho do sistema com cada tipo de arquitetura de controle. / [en] A study on a pan-tilt type target tracking system actuated by permanent magnet DC motors and assembled in a moving body is presented in this work. To achieve such objective, an experimental test bed is constructed and a simulation program is implemented. The mechanical model is derived and simulated in time domain. This approach allows using accurate non-linear equations to represent system behavior, otherwise infeasible in frequency domain. Although the system is modeled with rigid bodies, flexibility and structural damping due to the gearboxes are considered. Sensor errors, backlash in the gearboxes, dry and viscous friction, saturation limits for armature current and tension of the motors are also considered. A method to include the time delays for the control signal updates, as well as time delays due to sensor dynamic response, during the numerical integration of the equations of motion, is presented. Controllers that require no mathematical model of the plant are employed in the experimental test bed and in the simulation program. Three different control architectures are proposed, called in this work type 1, type 2 and type 3. Their complexity increases depending on the number of available sensors. The type 1 is applied to systems with only one sensor that provides the targets angular azimuth and elevation errors. If, besides this sensor, sensors to measure the relative angular positions between the mechanism links are available type 2 architecture is used. In addition, if sensors to measure inertial angular speeds are also available, type 3 architecture can be used. Finally, experimental and numerical results, comparing system performance with each control architecture are presented.

[pt] Suspensões correspondem a uma classe de materiais amplamente utilizada em uma grande variedade de aplicações e indústrias. Devido à sua extrema versatilidade, elas têm sido foco de inúmeros estudos nas últimas décadas. Suspensões também são muito flexíveis e podem apresentar diferentes propriedades reológicas e respostas macroscópicas dependendo da escolha dos parâmetros usados como entrada no sistema. Mais especificamente, a resposta reológica de suspensões está intimamente associada ao arranjo microestrutural das partículas que compõem o meio e a fatores externos, como o quão confinadas elas se encontram e a rigidez das partículas. No presente estudo, o efeito da rigidez, confinamento e vazão na microestrutura de suspensões altamente concentradas é avaliado usando Dinâmica Dissipativa de Partículas com Núcleo Modificado. Precedento este estudo principal, foram necessárias outras duas etapas para garantir um sistema de simulação confiável e representativo, que consistiu, essencialmente, na realização de estudos paramétricos para compreender e estimar os valores adequados para os parâmetros de interacção parede-partícula. O presente trabalho aborda estudos paramétricos realizados para auxiliar na escolha dos parâmetros de entrada para evitar a penetração de partículas em um sistema delimitado por paredes. Inicialmente um sistema mais simples, composto por solvente e paredes é construído e os parâmetros de interação e densidades de parede foram ajustados. Em seguida as interações são definidas para suspensões. Neste último caso, vários parâmetros desempenham um papel na penetração e a maneira tradicional de investigar esses efeitos seria exaustiva e demorada. Por isso, optamos por usar uma abordagem de Machine Learning para realizar este estudo. Uma vez ajustados os parâmetros, o estudo de confinamento pôde ser realizado. O objetivo principal deste estudo foi entender como a microestrutura de suspensões concentradas é afetada pela vazão, rigidez das partículas e confinamento. Verificou-se que partículas muito flexíveis sempre formam um aglomerado gigante independente da razão de confinamento; a diferença está em quão compactadas são as partículas. No caso de partículas rígidas, um confinamento mais forte leva à formação de aglomerados maiores. O estudo final aborda um estudo de aprendizado de máquina realizado para prever a reologia de suspensões não confinadas. Com este trabalho foi possível entender e ajustar parâmetros de simulação e desenvolver um domínio computacional que permite estudar sistematicamente efeitos do confinamento em suspensões. / [en] Suspensions correspond to a class of materials vastly used in a large set of applications and industries. Due to its extreme versatility, they have been the focus of numerous studies over the past decades. Suspensions are also very flexible and can display different rheological properties and macroscopic responses depending on the choice of parameters used as input in the system. More specifically, the rheological response of suspensions is intimately associated to the microstructural arrangement of the particles composing the medium and external factors, such as how strongly they are confined and particle rigidity. In the present study, the effect of particle rigidity, confinement and flow rate on the microstructure of highly concentrated suspensions is studied using CoreModified Dissipative Particle Dynamics. Preceding this main study, two other steps were necessary to guarantee a reliable and realistic simulation system, which consisted, essentially, on performing parametric studies to understand and estimate the appropriate values for wall-particle interaction parameters. The present work address parametric studies performed to assist the input parameters choice to prevent particle penetration in a wall-bounded system. Initially a simpler system, composed of solvent and walls, is built and the interaction parameters and wall densities were adjusted. Following, the interactions are set for suspensions. In the latter case multiple parameters play a role in penetration and the traditional way to investigate these effects would be exhaustive and time consuming. Hence, we choose to use a Machine Learning approach to perform this study. Once the parameters were adjusted, the study of confinement could be carried out. The main goal of this study was to understand how the microstructure of concentrated suspensions is affected by flow rate, particle rigidity and confinement. It was found that very soft particles always form a giant cluster regardless the confinement ratio; the difference being on how packed the particles are. In the rigid case, a stronger confinement leads the formation of larger clusters. The final study addresses a machine learning study carried out to predict the rheology of unconfined suspensions. The main contribution of this work is that it was possible to understand and adjust simulation parameters and develop a computational domain that enables to systematically study confinement effects on suspensions.

[pt] Simulações de dinâmica molecular coarse-grained foram usadas para investigar diferentes polímeros para o encapsulamento de um siRNA para a transposição em modelos de surfactante pulmonar. Os modelos consistiam em uma monocamada de 1,2-dipalmitoil-sn-glicero-3-fosfatidilcolina (DPPC) ou uma mistura 70:30 de DPPC e 1,2-palmitoil-sn-glicero-3-fosfo-(1 -sn-glicerol) (DPPG). Os nanocarreadores escolhidos para encapsular o siRNA anti-TNF foram o polietilenoglicol (PEG) e a polietilenimina (PEI). Acredita-se que o uso do siRNA para o silenciamento genético da citocina TNF pode ser terapêutico para algumas doenças inflamatórias, em especial, pulmonares. As simulações mostraram que as nanopartículas contendo apenas PEG promovem o arraste de fosfolipídeos do modelo de surfactante pulmonar, formando uma coroa lipídica. As nanopartículas contendo apenas PEI ou PEI e PEG causaram certas perturbações à monocamada, porém, nenhum colapso do filme durante a transposição da nanopartícula foi observado. O método de amostragem guarda-chuva foi utilizado para calcular a energia livre de Gibbs da transposição da nanopartícula através do modelo de surfactante pulmonar contendo apenas DPPC. A presença do PEG na nanopartícula causou uma diminuição na energia livre de Gibbs em relação ao siRNA nãoencapsulado, enquanto a presença do PEI não causou mudança. Os resultados implicam que siRNA encapsulado com ambos PEI e PEG melhoram a transposição do siRNA anti-TNF através de modelos de surfactante pulmonar na interface gáslíquido. / [en] Coarse-grained molecular dynamics was used to investigate different polymers to encapsulate the siRNA for its transposition through two lung surfactant models. These models consisted of a monolayer containing either 1,2-dipalmitoylsn-glycero-3-phosphatidylcholine (DPPC) or a 70:30 mixture of DPPC and 1,2- palmitoyl-sn-glycero-3-phospho-(1 -sn-glycerol) (DPPG). The nanocarriers chosen to encapsulate the anti-TNF siRNA were polyethylene glycol (PEG) and polyethyleneimine (PEI). It is believed that the use of siRNA to genetically silence the cytokine TNF might be therapeutic to treat several inflammatory diseases, in particular, pulmonary ones. The simulations showed that the nanoparticles containing PEG promoted the lipid depletion of the lung surfactant model by forming a lipid corona. The nanoparticles containing only PEI, or both PEG and PEI showed some perturbation of the lung surfactant model, however no collapse during the nanoparticle transposition was observed. The umbrella sampling method was used to calculate the Gibbs free energy of transposition through the pure DPPC lung surfactant model. The nanoparticles containing PEG showed a decreased Gibbs free energy as compared with naked siRNA, while PEI nanoparticles have not caused any change. The implication of this finding is that siRNA encapsulated with both PEI and PEG enhances the transposition of anti-TNF siRNA through lung surfactant models at the gas-liquid interface.

[pt] Os últimos leilões do pré-sal para exploração e produção de petróleo e gás no Brasil indicam que as operações de perfuração se tornarão mais intensas nos próximos anos. O processo de perfuração rotativo é amplamente utilizado para alcançar os reservatórios de petróleo e devido à relação diâmetro/comprimento do sistema de perfuração, o modo de vibração torcional está presente em quase todos os processos de perfuração, podendo chegar a um estado crítico indesejável: o fenômeno de stick-slip. Com o intuito de abordar este problema, o modo torcional é isolado e o stick-slip é observado em uma coluna de perfuração em escala reduzida completamente instrumentada. Durante o stick-slip, outro torque pode ser aplicado em uma posição intermediária da bancada de teste. O modelo matemático de parâmetros concentrados é obtido e o modelo é comparado com dados experimentais com o propósito de verificar se o modelo matemático representa o aparato experimental. Uma análise de estabilidade é feita usando o modelo validado com o objetivo de identificar soluções estáveis do sistema. Com isso, observou-se que existe uma faixa do parâmetro de bifurcação na qual soluções de equilíbrio e periódicas estáveis coexistem. Para uma dada situação de stick-slip na faixa de biestabilidade, duas estratégias de mitigação de vibração torcional foram consideradas e consistiram em impor perturbações no sistema por meio do torque na posição intermediária da bancada de teste: (i) torques aplicados apenas contra a direção de movimento do sistema, e (ii) torques aplicados em ambas as direções. As estratégias foram testadas numericamente e apresentaram eficiência de tal modo que o stick-slip foi completamente mitigado: as energias do sistema e o trabalho gerado pelo torque intermediário aplicado foram comparados com o propósito de avaliar a factibilidade e razoabilidade da estratégia. Experimentalmente, o sistema continuou a oscilar, porém apresentou uma significante redução na fase de stick mesmo com limitações de aplicações de torque. / [en] The latter round bids of the pre-salt for exploration and production of oil and natural gas in Brazil indicate the drilling operations will become more intense in coming years. The rotational drilling process is largely used to reach the oil reservoirs and because of diameter-to-length ratio of the drilling system, torsional vibration mode is present in most all drilling processes and may reach an undesired severe stage: the stick-slip phenomenon. In order to address this problem, the torsional vibration mode is isolated and the stick-slip is observed in a fully instrumented drill-string experimental set-up in this work. During this phenomenon, another torque may be applied on an intermediate position of the test bench. The lumped parameter mathematical model is obtained and it is compared to experimental data to validate whether the mathematical model represents the experimental apparatus. A stability analysis is performed using the validated mathematical model in order to identify stable solutions of the system. Therewith, one observed that there is a range of the bifurcation parameter in which stable equilibrium and periodic solutions may coexist. For a given stick-slip situation in bi-stability range, two mitigation strategies of torsional vibration were considered which consisted of imposing perturbations in the system via torques on the intermediate position of the test bench: (i) torques applied only against the direction of motion of the system, and (ii) torques applied in both directions. The strategies were tested numerically and presented eficiency so that the stickslip was completely mitigated: the energies of the system and the work created by the intermediate torque were compared in order evaluate the feasibility and reasonableness of the strategy. Experimentally, the system continued to oscillate, however it presented a significant reduction of stick phase even with limitations of torque applications.

[pt] O presente trabalho avalia a utilização de elementos enriquecidos para obtenção de cargas críticas, freqüências de vibração e seus respectivos modos de peças estruturais bidimensionais (flexão de placas retangulares sujeitas a compressão em seu plano). O método de aproximação empregado foi o de Rayleigh-Ritz voltado para o uso de elementos finitos convencionais enriquecidos com funções de deslocamentos adicionais internas e de contorno. As funções ditas internas são desenvolvidas de forma a não envolver deslocamentos e rotações nodais e no contorno. Já as funções ditas de contorno são concebidas de forma a envolver apenas deslocamentos internos e ao longo de um lado apenas, sem deslocamentos generalizados nodais. Para este estudo foram desenvolvidas duas famílias de funções, uma com termos adicionais trigonométricos e outra com termos adicionais polinomiais. Para o cálculo de cargas críticas e freqüências são utilizadas as matrizes de rigidez elástica, rigidez geométrica e de massa, introduzidas em problemas generalizados de autovalores. Resultados numéricos são obtidos através de procedimentos computacionais utilizando o software Maple. Verifica-se que as funções adicionais trigonométricas, embora mais satisfatórias que as polinomiais quanto à convergência, exigem maior esforço computacional. São comparados resultados de elementos para placas esbeltas (teoria de Kirchhoff), com três e quatro graus de liberdade por nó, onde o quarto grau de liberdade corresponde à derivada mista (torção). Mostra-se que as funções adicionais, não-nodais, requerem o uso do elemento com quatro graus de liberdade por nó, para se ter convergência no cálculo das cargas críticas e freqüências em situações gerais. Outros exemplos abordam preliminarmente a inclusão de efeitos de dano e ortotropia no material, visando a modelagem de lajes comprimidas e pilares com seções retangulares alongadas. Esta modelagem envolvendo combinação de funções adicionais gerais e elementos convencionais representa um passo no desenvolvimento de uma técnica aplicável à combinação de modos globais e localizados de instabilidade / [en] The focus of the present work is to developand evaluate enriched elements used to obtain critical loads, frequencies of vibration and respective modes for two-dimensional structural components (rectangular plates in bending under inplane compressive loading). The Rayleigh-Ritz approximation method has been employed, directed to the use of conventional finite elements enriched by internal and boundary additional displacements functions. The socalled internal functions are do not involve nodal and boundary displacements and rotations. The boundary functions are conceived to include displacements within the element and along one side, without involving any generalized nodal displacements. Two displacement function families were developed, the first with trigonometric additional terms and the second with polynomial additional terms. Critical loads and frequencies, and respective modes, are obtained by the use of elastic stifiness, geometric, and mass matrices, introduced in generalized eigenvalue problems. Numerical results are obtained by computational procedures using Maple software. The trigonometric additional functions, in spite of better convergence properties, demand greater computational effort. The basic elements are classical thin plate elements (Kirchhoff's theory) with three or four degrees of freedom per node, where the fourth degree of freedom corresponds to the mixed derivative (torsion). The results indicate that non- nodal additional functions require the use of elements with four freedom degrees by node to obtain convergence of critical loads and frequencies convergence in general situations. Other examples consist of preliminary approaches to include damage effects, in reinforced orthotropic plates, as modeling columns with wide rectangular sections and compressed slabs. The use of general additional functions combined with conventional elements represents a step on the development of a technique applicable to global and localized instability modes.

[pt] Apresenta-se a aplicação dos algoritmos genéticos para o tratamento do problema inverso em reconstituição de acidentes e análise de colisões com veículos terrestres de estrutura deformável. Define-se como, a partir de restrições impostas, das posições finais, e das deformações encontradas nos veículos após uma colisão, o algoritmo de otimização pode fornecer o conjunto de variáveis e parâmetros que mais provavelmente levam os veículos àquela condição. Todos os procedimentos desenvolvidos foram implementados em imulink/Matlab. Para resolver o problema, foi escolhida a técnica de otimização denominada algoritmo genético, que é indicado para solução de problemas complexos, que envolvem um grande número e variáveis e, conseqüentemente, espaços de soluções de dimensões elevadas. / [en] This work show an applicattion of the genetic algorithm to resolve the inverse problem of accident reconstitution and to analise colisions between vehicles of deformable structure. It is determined how, with imposing of restrictions, final positions and deformations found at vehicles after collision, the optimization algorithmcan give the set of variables and parameters that probably conduct the vehicles to true initial condition. All the developed procedures were implemented at Simulink/Matlab.The optimization technique chose to resolve the inverse problem was the genetic algorithm because it is the most popular to solve complex problems that have a very large number of variables and a elevate dimension space solutions.

[pt] Mudanças das propriedades dos solos devido à ação de carregamentos dinâmicos são responsáveis por danos significativos em geo-estruturas, tais como: barragens, estruturas de concentração, fundações, taludes, etc. A ocorrência do fenômeno da liquefação, em materiais suscetíveis como areias fofas saturadas, representa um tipo de resposta desastrosa de solos. O termo liquefação tem sido empregado para descrever uma variedade de fenômenos no qual tem em comum o desenvolvimento de altas poropressões em materiais saturados sem coesão devido a carregamentos monotônicos , transientes ou ciclios. A previsão da liquefação depende de uma adequada análise do comportamento não-drenado do material, em termos do incremento de poropressões e da perda da rigidez da mistura sólido-fluido, durante e após o período de movimento. O estabelecimento das equações governantes é essencial para elaboração de um modelo matemático realista para descrever o comportamento físico deste fenômeno. As equações a srem consideradas são: equação de movimento da fase sólida, a equação do movimento da mistura sólido-fluido , a equação de continuidade da fase fluida, as equações de acoplamento das fases e as equações constitutivas desses materiais. Nesta tese a resposta dinâmica do solo foi investigada numericamente mediante a técnica dos elementos finitos. A discretização espacial das equações governantes foi feita através de método de Galerkin e a discretização temporal pelo método de Newmark Generalizado. Um modelo constitutivo elasto-plástico foi considerado para descrever o comportamento mecânico da fase sólida, desenvolvido a partir de conceitos da generealização da teoria da plasticidade, que apresenta algumas vantagens em relação aos outros modelos baseados na teoria da plasticidade clássica. A implementação computacional foi escrito em fortran 90. Exemplos numéricos analisados nesta tese comprovam tanto a eficiência do modelo constitutivo na predição do comportamento do solo sobre liquefação como a confiabilidade do programa computacional elaborado nesta pesquisa, em termos da rapidez de processamento e da boa precisão dos resultados, quando comparados com soluções analíticas e outros valores numéricos obtidos por vários autores e diferentes modelos constitutivos. / [en] Changes in soil properties due to the action of dynamic loads are responsible for significant damage of geo-structures such as dams, retaining structures,building foundations, slopes, etc. The occurrence of liquefaction phenomena in susceptible materials, such as loose saturated, represents a type of disastrous response of soil, the term liquefaction has been used to refer to a group of phenomena wich have in common the development of high pore pressures in saturated cohesionless mterial due to monotonic, transient, or cyclic loads. The prediction of soil liquefaction depends of an adequate analysis of the behavior of undrained materials, in terms of increase of pore water pressure and weakening of the solid-fluid mixture, during and after the periodic motion. The establishment of the governing equations is essential to provide a realistic mathematical model to describe the physical behavior of this phenomenon. The system of equations to be considered are: the equilibrium equation of the solid phase, the equilibrium equation of the solid-fluid mixture, the conservation mass of the fluid phase, the coupling equation of phases, and the conservation equations of materials. In this thesis the soil dynamic response was numerically investigated by the finite element method. To obtain the spatial discretization in time was the Generalized Newmark method. An elastic-plastic constitutive model was used to describe the mechanical behavior of the solid phase. This model was developed in the framework of the generalized theory of plasticity, wich has some advantages when compared with other models based on the classical plasticity theory. The computacional implementation was written in fortran 90. Numerical examples considered in this thesis demonstrate the efficiency of the constitutive model to simulated the predicted behavior of soil under liquefaction as well as the reliability of the software developed in this research, in terms of computational effort and good accuracy of the results, when compared with some analytical solutions and other numerical values obtained by various authors and different constitutive models.