Desenvolvimento de um método numérico para simular escoamentos viscoelásticos axissimétricos com superfícies livres: Modelo PTT

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

Simulação numérica do comportamento dinâmico de uma laje

Nesta dissertação pretende-se simular o comportamento dinâmico de uma laje de betão armado aplicando o Método de Elementos Finitos através da sua implementação no programa FreeFEM++. Este programa permite-nos a análise do modelo matemático tridimensional da Teoria da Elasticidade Linear, englobando a Equação de Equilíbrio, Equação de Compatibilidade e Relações Constitutivas. Tratando-se de um problema dinâmico é necessário recorrer a métodos numéricos de Integração Directa de modo a obter a resposta em termos de deslocamento ao longo do tempo. Para este trabalho escolhemos o Método de Newmark e o Método de Euler para a discretização temporal, um pela sua popularidade e o outro pela sua simplicidade de implementação. Os resultados obtidos pelo FreeFEM++ são validados através da comparação com resultados adquiridos a partir do SAP2000 e de Soluções Teóricas, quando possível.

Métodos numéricos para resolução de equações diferenciais ordinárias lineares baseados em interpolação por spline

In this work we have elaborated a spline-based method of solution of inicial value problems involving ordinary differential equations, with emphasis on linear equations. The method can be seen as an alternative for the traditional solvers such as Runge-Kutta, and avoids root calculations in the linear time invariant case. The method is then applied on a central problem of control theory, namely, the step response problem for linear EDOs with possibly varying coefficients, where root calculations do not apply. We have implemented an efficient algorithm which uses exclusively matrix-vector operations. The working interval (till the settling time) was determined through a calculation of the least stable mode using a modified power method. Several variants of the method have been compared by simulation. For general linear problems with fine grid, the proposed method compares favorably with the Euler method. In the time invariant case, where the alternative is root calculation, we have indications that the proposed method is competitive for equations of sifficiently high order.

Unicidade e discretização para problemas de valor inicial

This paper has two objectives: (i) conducting a literature search on the criteria of uniqueness of solution for initial value problems of ordinary differential equations. (ii) a modification of the method of Euler that seems to be able to converge to a solution of the problem, if the solution is not unique

O chute com o membro dominante e não dominante realizado com a bola parada e em deslocamento no futsal

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

Métodos numéricos para o retoque digital

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

Modelo distribuído aplicado à análise de evaporadores do tipo tubo aletado

Pós-graduação em Engenharia Mecânica - FEIS

Algoritmos de Planificación y Seguimiento de Trayectorias para Robots Agrícolas

La base para realizar cualquier tarea agrícola mediante robots, es la planificación y seguimiento de rutas o trayectorias. Así, el objetivo de esta investigación es desarrollar e implementar algoritmos de seguimiento y planificación (global y local) de trayectorias de robots agrícolas. La planificación global se realizó mediante el algoritmo A* aplicado sobre mapas de cultivo y la planificación local se realizó aplicando A* sobre un mapa 2D obtenido a partir de imágenes 3D de los obstáculos encontrados en el camino. En cuanto el seguimiento de trayectorias, esta se realizó implementando una aproximación numérica de la trayectoria mediante el método de Euler. Los parámetros correspondientes a la dinámica del controlador de la trayectoria del robot fueron obtenidos mediante algoritmos genéticos. El mapa 3D fue generado a partir del sensor Kinect de Microsoft y sus datos procesados usando Matlab 2010b. Los resultados preliminares muestran que es posible implementar estos algoritmos en pequeños robots diseñados para cultivos hilerados. Proveyendo así, una metodología robusta que permite seguir las rutas asignadas con errores inferiores a RMSE=0.1m en trayectorias de 30m.

Difusão do cloreto de sódio no processo de salga de queijos: modelagem matemática com o emprego do método de elementos finitos

A transferência de um soluto (cloreto de sódio), através de uma matriz sólida tridimensional (queijo) foi estudada aplicando-se o método de elementos finitos. A formulação variacional (Galerkin) do problema diferencial (modelo de difusão) teve como base teórica a 2ª lei de Fick. Os procedimentos para integração no tempo foram o de Crank-Nicolson e o de Euler-modificado, que foram escolhidos por apresentarem estabilidade incondicional. O programa computacional desenvolvido mostrou-se versátil para resolver situações de amostragem em condições mais realistas e pode ser aplicado para geometrias complexas. O modelo proposto permitiu uma boa estimativa do ganho de sal no queijo, usando um coeficiente de difusão cujo valor pode ser obtido por extrapolação de dados experimentais. A aplicação do método numérico (MEF), com o esquema de Crank-Nicolson, na simulação da difusão do cloreto de sódio na salga de queijos, mostrou boa aproximação quando os resultados foram comparados com os valores experimentais encontrados na literatura especializada.

A equação de Euler e a IS: estimações para o Brasil

Neste artigo estudamos a relação entre a taxa de juros e o hiato do produto no Brasil através da estimação de modelos Novo-Keynesianos. Para tanto, estimamos os modelos por três métodos: (1) método generalizados dos momentos, (2) máxima verossimilhança utilizando dados de expectativa divulgados pelo Banco Central do Brasil e (3) máxima verossimilhança utilizando variáveis de expectativa estimadas por um modelo VAR. As conclusões são altamente dependentes do método de estimação. Ao utilizar (1), os resultados indicam uma relação espúria entre a taxa de juros real e o hiato. Entretanto, as conclusões originadas de (2) indicam que somente o hiato defasado seria uma variável relevante para a nossa especificação. Ao estimar o modelo por (3), as estimativas corroboram os resultados obtidos com o método generalizados dos momentos.

Aerodinâmica de automóveis baseado no método dos contornos virtuais

Este trabalho desenvolve um método numérico para a solução de escoamentos bidimensionais em torno de geometrias automobilísticas utilizando o método de diferenças finitas. O código computacional resolve as equações de Navier-Stokes e de Euler para uma distribuição adequada dos pontos discretos na malha. O método de integração empregado baseia-se no esquema explícito de Runge-Kutta de 3 estágios para as equações da quantidade de movimento e no de sub-relaxações sucessivas para a pressão na base Gauss-Seidel. Utilizou-se a técnica dos contornos virtuais em coordenadas cartesianas para resolver o escoamento sobre uma geometria simplificada, com a superfície coincidente com a malha computacional, e uma geometria automobilística mais complexa (BMW). Para a certificação da técnica empregada, optou-se pela utilização da teoria do escoamento potencial e pela comparação com dados experimentais encontrados na literatura e outros coletados em túnel de vento em escala reduzida. Houve dificuldade nesta comparação devido à falta de artigos relativos às simulações numéricas de escoamentos sobre automóveis e na aplicação da técnica dos contornos virtuais em geometrias complexas. Os resultados foram satisfatórios, com boas perspectivas para trabalhos futuros, contribuindo assim para o desenvolvimento da área.

Simulação numérica da interação entre uma nuvem incidente de vórtices e uma aerofólio através do método dos vórtices

The study of aerodynamic loading variations has many engineering applications, including helicopter rotor blades, wind turbines and turbo machinery. This work uses a Vortex Method to make a lagrangian description of the a twodimensional airfoil/ incident wake vortex interaction. The flow is incompressible, newtonian, homogeneus and the Reynolds Number is 5x105 .The airfoil is a NACA 0018 placed a angle of attack of the 0° and 5°simulates with the Painel Method with a constant density vorticity panels and a generation poit is near the painel. The protector layer is created does not permit vortex inside the body. The vortex Lamb convection is realized with the Euler Method (first order) and Adans-Bashforth (second order). The Random Walk Method is used to simulate the diffusion. The circular wake has 366 vortex all over positive or negative vorticity located at different heights with respect to the airfoil chord. The Lift was calculated based in the algorithm created by Ricci (2002). This simulation uses a ready algorithm vatidated with single body does not have a incident wake. The results are compared with a experimental work The comparasion concludes that the experimental results has a good agrement with this papper

Simulação numérica de escoamentos hipersônicos sobre corpos rombudos pelo método de elementos finitos

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

Deconvolução de Euler: passado, presente e futuro – um tutorial

Neste tutorial apresentamos uma revisão da deconvolução de Euler que consiste de três partes. Na primeira parte, recordamos o papel da clássica formulação da deconvolução de Euler 2D e 3D como um método para localizar automaticamente fontes de campos potenciais anômalas e apontamos as dificuldades desta formulação: a presença de uma indesejável nuvem de soluções, o critério empírico usado para determinar o índice estrutural (um parâmetro relacionado com a natureza da fonte anômala), a exeqüibilidade da aplicação da deconvolução de Euler a levantamentos magnéticos terrestres, e a determinação do mergulho e do contraste de susceptibilidade magnética de contatos geológicos (ou o produto do contraste de susceptibilidade e a espessura quando aplicado a dique fino). Na segunda parte, apresentamos as recentes melhorias objetivando minimizar algumas dificuldades apresentadas na primeira parte deste tutorial. Entre estas melhorias incluem-se: i) a seleção das soluções essencialmente associadas com observações apresentando alta razão sinal-ruído; ii) o uso da correlação entre a estimativa do nível de base da anomalia e a própria anomalia observada ou a combinação da deconvolução de Euler com o sinal analítico para determinação do índice estrutural; iii) a combinação dos resultados de (i) e (ii), permitindo estimar o índice estrutural independentemente do número de soluções; desta forma, um menor número de observações (tal como em levantamentos terrestres) pode ser usado; iv) a introdução de equações adicionais independentes da equação de Euler que permitem estimar o mergulho e o contraste de susceptibilidade das fontes magnéticas 2D. Na terceira parte apresentaremos um prognóstico sobre futuros desenvolvimentos a curto e médio prazo envolvendo a deconvolução de Euler. As principais perspectivas são: i) novos ataques aos problemas selecionados na segunda parte deste tutorial; ii) desenvolvimento de métodos que permitam considerar interferências de fontes localizadas ao lado ou acima da fonte principal, e iii) uso das estimativas de localização da fonte anômala produzidas pela deconvolução de Euler como vínculos em métodos de inversão para obter a delineação das fontes em um ambiente computacional amigável.

Análise estrutural de Tie-Rod pelo método dos elementos finitos e através de cálculos analíticos

Currently one of the great concerns of the aeronautical industry is in relation to the security and integrity of the aircraft and its equipments / components when under critical flight maneuvers such as during landing / takeoff and emergency maneuvers. The engineers, technicians and scientists are constantly developing new techniques and theories to reduce the design time and testing, ir order to minimize costs. More and more the Finite Element Method is used in the structural analysis of a project as well as theories based on experimental results. This work aimed to estimate the critical load to failure for tensile, compression and buckling of the Tie-Rod, a fixture aircraft widely used on commercial aircrafts. The analysis was performed by finite element method with the assistance of software and by analytical calculations. The results showed that the Finite Element Method provides relative accuracy and convenience in the calculations, indicating critical load values slightly lower than those found analytically for tension and compression. For buckling, the Finite Element Method indicates a critical load very similar to that found analytically following empirical theories, while Euler's theory results in a slightly higher value. The highest risk is to fail by buckling, but the geometric irregularity of Tie-Rod pieces makes difficult the calculations, therefore a practical test must be done before validation of the results