Método de Otimização Topológica em Estruturas Tridimensionais

The topology optimization problem characterize and determine the optimum distribution of material into the domain. In other words, after the definition of the boundary conditions in a pre-established domain, the problem is how to distribute the material to solve the minimization problem. The objective of this work is to propose a competitive formulation for optimum structural topologies determination in 3D problems and able to provide high-resolution layouts. The procedure combines the Galerkin Finite Elements Method with the optimization method, looking for the best material distribution along the fixed domain of project. The layout topology optimization method is based on the material approach, proposed by Bendsoe & Kikuchi (1988), and considers a homogenized constitutive equation that depends only on the relative density of the material. The finite element used for the approach is a four nodes tetrahedron with a selective integration scheme, which interpolate not only the components of the displacement field but also the relative density field. The proposed procedure consists in the solution of a sequence of layout optimization problems applied to compliance minimization problems and mass minimization problems under local stress constraint. The microstructure used in this procedure was the SIMP (Solid Isotropic Material with Penalty). The approach reduces considerably the computational cost, showing to be efficient and robust. The results provided a well defined structural layout, with a sharpness distribution of the material and a boundary condition definition. The layout quality was proporcional to the medium size of the element and a considerable reduction of the project variables was observed due to the tetrahedrycal element

Modelagem de máquina pneumática recolhedora de frutos de café em terreiro usando análise por elementos finitos

Teve-se o objetivo de desenvolver um modelo matemático por meio de análise de elementos finitos, utilizando o programa computacional ANSYS®, versão 5.7, para otimizar o projeto de máquina recolhedora de frutos de café no terreiro. A modelagem da máquina foi realizada com base no levantamento das características aerodinâmicas dos frutos de café e da vazão de ar necessária para o transporte pneumático dos frutos. Foram obtidas, experimentalmente, as pressões estáticas nos dutos da máquina, sendo esses valores comparados com os resultados determinados pelo programa ANSYS, no intuito de validar o modelo. Com base nos resultados numéricos obtidos, concluiu-se que a modelagem desenvolvida apresentou resultados próximos aos determinados experimentalmente, obtendo erro relativo médio nos valores simulados de pressão de 9,2%. Por meio da modelagem, identificaram-se faixas de pressão que dificultariam o transporte pneumático dos frutos de café em alguns pontos da máquina. Esses problemas foram corrigidos e, com isso, o fluxo de ar proporcionado pelo ventilador foi suficiente para succionar os frutos de café no terreiro e transportá-los para dentro do reservatório da máquina. A modelagem desenvolvida atendeu às necessidades propostas no trabalho para o recolhimento dos frutos de café utilizando transporte pneumático eficientemente.

Simulação do comportamento não linear de materiais quase-frágeis via elementos finitos com fissura embebida

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

Análise da execução de um algoritmo numérico cbs para simulação de escoamento utilizando elementos finitos

This work intends to analyze the application and execution time of a numerical algorithm that simulates incompressible and isothermal flows. It was used the explicit scheme of the Characteristic Based Split (CBS) algorithm and the Artificial Compressibility (AC) scheme for coupling pressure-velocity equations. The discretization was done with the finite element method using a bilinear elements grid. The free software GNU-Octave was used for implementation and execution of routines. The results were analyzed using the classic lid-driven cavity problem. This work shows results for tests with several Reynolds' number. The results for these tests show a good agreement when compared with previous ones obtained from bibliography. The code runtime's analysis shows yet that the matrix's assembly is the part of greater consumption time in the implementation.

Método de detecção interferométrica de fase, com baixa profundidade de modulação, aplicado à medição de deslocamentos nanométricos em atuadores e mini-manipuladores piezoelétricos

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

Estudo de campo elétrico em linha de transmissão utilizando o método dos elementos de contorno

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

Aplicação da técnica de elementos finitos na estratificação de armazenadores térmicos com barreiras de meios porosos

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

Inclusão de pilares numa formulação do método dos elementos de contorno para análise linear de pavimentos de edifícios

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

Implementação do método das características na modelagem de problemas de convecção natural em cavidades cilíndricas

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

Análise da estampagem de chapas metálicas utilizando elementos finitos

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

Análise da forma da partícula na aplicação do Método dos Elementos Discretos (DEM)

O presente trabalho analisa o efeito da forma da partícula por meio do método dos Elementos Discretos (DEM). São analisadas quatro partículas com formas e esfericidades diferentes – uma esférica e três não-esféricas. O ângulo de repouso é o principal instrumento de avaliação dos escoamentos. A análise inicia com a calibração da partícula esférica usando dados da literatura, em seguida, estes parâmetros calibrados são empregados para a simulação dos escoamentos das partículas não-esféricas. São feitas comparações dos esforços computacionais e, estas informações, são usadas para verificar as vantagens que a partícula esférica oferece sobre as outras três formas. Nesse cenário, procedimentos foram desenvolvidos para ajudar no processo de calibração dos ângulos de repousos, tendo como base, o conhecimento da sensibilidade de alguns parâmetros DEM. Os resultados mostram a influência das partículas não-esféricas e, principalmente, que é possível obter com a partícula esférica, escoamentos similares aos das partículas não-esféricas, com ganho computacional.

Modelagem 1D e 2,5D de dados do método CSEM marinho em meios com anisotropia transversal inclinada

Neste trabalho apresentamos a solução do campo eletromagnético gerado por um dipolo elétrico horizontal em meios transversalmente isotrópicos com eixo de simetria vertical (TIV) e com eixo de simetria inclinado (TII). Para modelos unidimensionais, o campo eletromagnético foi obtido por duas metodologias distintas: (1) solução semi-analítica das equações de Maxwell com auxílio de potenciais vetores no caso TIV e (2) em modelos com anisotropia transversal inclinada o campo eletromagnético foi separado em primário e secundário, e então, o campo secundário foi calculado pelo método de elementos finitos no domínio (k_x, k_y, z) da transformada de Fourier. Para estruturas bidimensionais, foi aplicada a mesma metodologia usado nos modelos TII unidimensionais, onde o campo secundário foi calculado pelo método de elementos finitos no domínio (x, k_y, z), da transformada de Fourier, com a utilização de malhas não estruturadas para discretização dos modelos. Estas respostas foram usados para avaliar os efeitos da anisotropia elétrica nos dados CSEM marinho 1D e 2,5D.

Modelagem 2,5D dos campos usados no Método Eletromagnético a Multi-Frequência - EMMF

Esta tese mostra a modelagem 2,5D de dados sintéticos do Método Eletromagnético a Multi-frequência (EMMF). O trabalho é apresentado em duas partes: a primeira apresenta os detalhes dos métodos usados nos cálculos dos campos gerados por uma bobina horizontal de corrente colocada sobre a superfície de modelos bidimensionais; e a segunda, usa os resultados obtidos para simular os dados medidos no método EMMF, que são as partes real e imaginária da componente radial do campo magnético gerado pela bobina. Nesta segunda parte, observamos o comportamento do campo calculado em diversos modelos, incluindo variações nas propriedades físicas e na geometria dos mesmos, com o intuito de verificar a sensibilidade do campo observado com relação às estruturas presentes em uma bacia sedimentar. Com esta modelagem, podemos observar as características dos dados e como as duas partes, real e imaginária, contribuem com informações distintas e complementares. Os resultados mostram que os dados da componente radial do campo magnético apresentam muito boa resolução lateral, mesmo estando a fonte fixa em uma única posição. A capacidade desses dados em distinguir e resolver estruturas alvo será fundamental para o trabalho futuro de inversão, bem como para a construção de seções de resistividade aparente.

Ensaios sobre a contribuição do ar nos resultados da modelagem do modo transverso-elétrico do método magnetotelúrico

Apresentamos a solução numérica do modelo da falha infinita para o modo TE utilizando o método dos elementos finitos para investigar o efeito do ar na solução TE. A comparação de nossa solução para o campo elétrico normalizado com a solução de Weaver evidenciou uma discrepância esperada devido ao efeito do ar, que foi intencionalmente negligenciado por Weaver, com o intento de facilitar a solução do problema. O problema analítico, levando em conta a presença do ar, é substancialmente mais difícil. Uma solução do problema real, parcialmente analítica e parcialmente numérica, foi apresentada por Sampaio. Devido a uma discrepância inesperada entre a nossa solução por elementos finitos e a solução exata de Sampaio, recalculamos a sua solução. Iterando a parte numérica da solução de Sampaio até ordens mais altas pudemos verificar que esta apresenta comportamento divergente. Para examinar o efeito do ar na resistividade aparente, comparamos três perfis de resistividade aparente obtidos a partir de nossa solução por elementos finitos com os mesmos perfis obtidos pela solução de Weaver. A discrepância encontrada entre estas soluções foi muito pequena, situando-se na faixa de erro instrumental de campo.

Método dos elementos de contorno na resolução do problema de segunda ordem em placas delgadas

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