Development of numerical methodologies for parameter identification and shape optimization in metal forming simulations

Por parte da indústria de estampagem tem-se verificado um interesse crescente em simulações numéricas de processos de conformação de chapa, incluindo também métodos de engenharia inversa. Este facto ocorre principalmente porque as técnicas de tentativa-erro, muito usadas no passado, não são mais competitivas a nível económico. O uso de códigos de simulação é, atualmente, uma prática corrente em ambiente industrial, pois os resultados tipicamente obtidos através de códigos com base no Método dos Elementos Finitos (MEF) são bem aceites pelas comunidades industriais e científicas Na tentativa de obter campos de tensão e de deformação precisos, uma análise eficiente com o MEF necessita de dados de entrada corretos, como geometrias, malhas, leis de comportamento não-lineares, carregamentos, leis de atrito, etc.. Com o objetivo de ultrapassar estas dificuldades podem ser considerados os problemas inversos. No trabalho apresentado, os seguintes problemas inversos, em Mecânica computacional, são apresentados e analisados: (i) problemas de identificação de parâmetros, que se referem à determinação de parâmetros de entrada que serão posteriormente usados em modelos constitutivos nas simulações numéricas e (ii) problemas de definição geométrica inicial de chapas e ferramentas, nos quais o objetivo é determinar a forma inicial de uma chapa ou de uma ferramenta tendo em vista a obtenção de uma determinada geometria após um processo de conformação. São introduzidas e implementadas novas estratégias de otimização, as quais conduzem a parâmetros de modelos constitutivos mais precisos. O objetivo destas estratégias é tirar vantagem das potencialidades de cada algoritmo e melhorar a eficiência geral dos métodos clássicos de otimização, os quais são baseados em processos de apenas um estágio. Algoritmos determinísticos, algoritmos inspirados em processos evolucionários ou mesmo a combinação destes dois são usados nas estratégias propostas. Estratégias de cascata, paralelas e híbridas são apresentadas em detalhe, sendo que as estratégias híbridas consistem na combinação de estratégias em cascata e paralelas. São apresentados e analisados dois métodos distintos para a avaliação da função objetivo em processos de identificação de parâmetros. Os métodos considerados são uma análise com um ponto único ou uma análise com elementos finitos. A avaliação com base num único ponto caracteriza uma quantidade infinitesimal de material sujeito a uma determinada história de deformação. Por outro lado, na análise através de elementos finitos, o modelo constitutivo é implementado e considerado para cada ponto de integração. Problemas inversos são apresentados e descritos, como por exemplo, a definição geométrica de chapas e ferramentas. Considerando o caso da otimização da forma inicial de uma chapa metálica a definição da forma inicial de uma chapa para a conformação de um elemento de cárter é considerado como problema em estudo. Ainda neste âmbito, um estudo sobre a influência da definição geométrica inicial da chapa no processo de otimização é efetuado. Este estudo é realizado considerando a formulação de NURBS na definição da face superior da chapa metálica, face cuja geometria será alterada durante o processo de conformação plástica. No caso dos processos de otimização de ferramentas, um processo de forjamento a dois estágios é apresentado. Com o objetivo de obter um cilindro perfeito após o forjamento, dois métodos distintos são considerados. No primeiro, a forma inicial do cilindro é otimizada e no outro a forma da ferramenta do primeiro estágio de conformação é otimizada. Para parametrizar a superfície livre do cilindro são utilizados diferentes métodos. Para a definição da ferramenta são também utilizados diferentes parametrizações. As estratégias de otimização propostas neste trabalho resolvem eficientemente problemas de otimização para a indústria de conformação metálica.

A importância da matemática e da modelação matemática no mestrado integrado em ciências farmacêuticas

A Matemática e as Ciências Farmacêuticas encontram-se relacionadas desde há muito, no entanto, foi a partir do séc. XVII, período de notável agitação cultural e científico que os métodos experimentais foram sustentados com cálculos matemáticos. Esta ciência e as técnicas de modelagem matemática tornaram-se numa ferramenta amplamente utilizada, de tal modo, que nos dias de hoje são consideradas como fundamentais na generalidade das profissões e em especial nas Ciências Farmacêuticas. Contudo, para muitos ainda não é vista como fundamental e essencial para a formação de futuros farmacêuticos. Deste modo, pretende-se demonstrar como a Matemática e as técnicas de modelagem se tornaram ao longo dos anos nesta poderosa ferramenta. Quer pelos instrumentos, quer pelas competências que nos proporcionam. Pretende-se também, com recurso aos conteúdos programáticos desta unidade curricular, avaliar se os conhecimentos, sistemas de avaliação e distribuição da carga horária são efetuados de forma homogénea pelas diferentes instituições portuguesas, públicas ou privadas que lecionam o Mestrado Integrado em Ciências Farmacêuticas. Verificou-se que a Matemática é uma ciência plena de capacidades e recursos e que estabelece uma relação interdisciplinar com as Ciências Farmacêuticas. Quer pela componente utilitária, quer pela componente formativa que proporciona. A análise dos conteúdos programáticos demonstra que apesar de serem transversais, as Universidades que não lecionam Sistemas de Equações Lineares e Equações diferenciais deveriam faze-lo e também realizarem um melhor controlo da carga horária por temática.

Análise dinâmica de vigas assimétricas em curva com parede fina e aberta submetidas a ação de carga móveis

Fundamentalmente, o presente trabalho faz uma análise elástica linear de pontes ou vigas curvas assimétricas de seção transversal aberta e de parede fina, com propriedades físicas, geométricas e raio de curvatura constantes ao longo do eixo baricêntrico. Para tanto, utilizaram-se as equações diferenciais de VLASOV considerando o acoplamento entre as deformações nas direções vertical, transversal, axial de torcão nal. Na solução do sistema de quatro equações com derivadas parciais foi utilizado um apropriado método numérico de integração (Diferenças Finitas Centrais). A análise divide-se, basicamente, em dois tipos: análise DINÂMICA e ESTATICA. Ambas são utilizadas também na determinação do coeficiente de impacto (C.M.D.). A primeira refere-se tanto na determinação das características dinâmicas básicas (frequências naturais e respectivos modos de vibração), como também na determinação da resposta dinâmica da viga, em tensões e deformações, para cargas móveis arbitrárias. Vigas com qualquer combinação das condições de contorno, incluindo bordos rotulados e engastados nas três direções de flexão e na torção, são consideradas. 0s resultados da análise teórica, obtidos pela aplicação de programas computacionais implementados em microcomputador (análise estática) e no computador B-6700 (análise dinâmica), são comparados tanto com os da bibliografia técnica como também com resultados experimentais, apresentando boa correlação.

Otimização de geometrias aerodinâmicas utilizando métodos inversos

O objetivo deste trabalho é a obtenção de uma técnica para a modelagem otimizada de corpos submetidos a fluxos de alta velocidade, como aerofólios em escoamentos transônicos e outras geometrias aerodinâmicas. A técnica é desenvolvida através de expansões em séries de Fourier para um conjunto de equações diferenciais com interrelação com as condições de contorno, sendo uma equação para a parte superior e outra para a parte inferior do aerofólio. O método de integração temporal empregado baseia-se no esquema explícito de Runge-Kutta de 5 estágios para as equações da quantidade de movimento e na relação de estado para a pressão. Para a aproximação espacial adota-se um esquema em volumes finitos no arranjo co-localizado em diferenças centrais. Utiliza-se dissipação artificial para amortecer as frequências de alta ordem do erro na solução das equações linearizadas. A obra apresenta a solução de escoamentos bi e tridimensionais de fluidos compressíveis transônicos em torno de perfis aerodinâmicos. Os testes num´ericos são realizados para as geometrias do NACA 0012 e 0009 e asas tridimensionais usando as equações de Euler, para número de Mach igual a 0.8 e ® = 0o. Os resultados encontrados comparam favoravelmente com os dados experimentais e numéricos disponíveis na literatura.

Uma implementação geral de elementos de contorno e sua aplicação a problemas de flexão de placas

o presente trabalho aborda a aplicação do método dos elementos de contorno (MEC) para solução de problemas de flexão linear e geometricamente não-linear de placas semiespessas. Os modelos de placa empregados consideraram a influência do cisalhamento através de teorias de primeira ordem, especificamente as de Mindlin e Reissner. Uma formulação integral unificada dos modelos de placa utilizados é desenvolvida para o operador de Navier do problema, onde foram mantidos alguns termos de ordem superior no tensor deformação de Green. A formulação integral do problema de membrana acoplado ao de flexão é igualment desenvolvida, levando a um sistema de equações integrais não-lineares que descreve completamente problemas de placas que envolvem grandes deslocamentos. Estas equações podem ser particularizadas para problemas de flexão linear e estabilidade elástica. Tendo em vista a necessidade de se considerar derivadas dos deslocamentos translacionais, as equações integrais correspondentes ao gradiente dos deslocamentos também foram deduzidas, caracterizando uma formulação hipersingular. o método empregado para solução numérica do sistema de equações integrais foi o método direto dos elementos de contorno. Um tratamento das integrais fortemente singulares presentes nas equações foi realizado, baseado em expansões assint6ticas dos núcleos. Deste procedimento resulta uma abordagem regularizada que emprega apenas quadraturas padrão de Gauss-Legendre.

Modelagem estocástica de sistemas hierárquicos de estoques

Esta dissertação apresenta a Modelagem Estocástica de Sistemas Hierárquicos de Estoques (SHEs) como um instrumento de uso gerencial, orientado para a maximização do nível de serviço oferecido aos consumidores finais e a minimização do nível global de estoques em uma cadeia de abastecimento. Tal modelagem envolve otimização de sistemas não-lineares com integrais indefinidas, normalmente tratados através de simulação, dada a impossibilidade de estabelecer uma solução analítica para o problema. Neste trabalho, diferentemente, a solução é obtida através da otimização do sistema, a partir de uma modelagem que estabelece os níveis de estoque-alvo, os pontos de reposição de estoques, os níveis de serviço oferecidos aos consumidores finais e o nível ótimo global de estoques. A sua operacionalização utiliza um aplicativo computacional, projetado e desenvolvido especificamente para este fim. Estruturas hierárquicas, quando otimizadas localmente, em cada um de seus níveis, podem apresentar resultados subótimos em nível global. Assim, não é suficiente identificar os ótimos locais para compor uma solução otimizadora para o sistema. Para resolver esse problema, a maior parte dos pesquisadores utiliza simulação. A qualidade das soluções assim obtidas geralmente é inferior à qualidade da solução obtida mediante otimização global do sistema A presente modelagem tem como objetivo gerar soluções mais qualificadas para o problema da otimização de SHEs. Nesta dissertação, as principais abordagens para otimização de SHEs são apresentadas e analisadas, com destaque para o trabalho de Rosenbaum (1981), cujos resultados, obtidos através de simulação, são comparados com os resultados deste trabalho.

ONNIS-GI: uma rede neural oscilatória para segmentação de imagens implementada em arquitetura maciçamente paralela

A presente tese apresenta a concepção de uma rede neural oscilatória e sua realização em arquitetura maciçamente paralela, a qual é adequada à implementação de chips de visão digitais para segmentação de imagens. A rede proposta, em sua versão final, foi denominada ONNIS-GI (Oscillatory Neural Network for Image Segmentation with Global Inhibition) e foi inspirada em uma rede denominada LEGION (Locally Excitatory Globally Inhibitory Oscillator Network), também de concepção recente. Inicialmente, é apresentada uma introdução aos procedimentos de segmentação de imagens, cujo objetivo é o de situar e enfatizar a importância do tema abordado dentro de um contexto abrangente, o qual inclui aplicações de visão artificial em geral. Outro aspecto abordado diz respeito à utilização de redes neurais artificiais em segmentação de imagens, enfatizando as denominadas redes neurais oscilatórias, as quais têm apresentado resultados estimulantes nesta área. A implementação de chips de visão, integrando sensores de imagens e redes maciçamente paralelas de processadores, é também abordada no texto, ressaltando o objetivo prático da nova rede neural proposta. No estudo da rede LEGION, são apresentados resultados de aplicações originais desenvolvidas em segmentação de imagens, nos quais é verificada sua propriedade de separação temporal dos segmentos. A versão contínua da rede, um arranjo paralelo de neurônios baseados em equações diferenciais, apresenta elevada complexidade computacional para implementação em hardware digital e muitos parâmetros, com procedimento de ajuste pouco prático. Por outro lado, sua arquitetura maciçamente paralela apresenta-se particularmente adequada à implementação de chips de visão analógicos com capacidade de segmentação de imagens. Com base nos bons resultados obtidos nas aplicações desenvolvidas, é proposta uma nova rede neural, em duas versões, ONNIS e ONNIS-GI, as quais suplantam a rede LEGION em diversos aspectos relativos à implementação prática. A estrutura dos elementos de processamento das duas versões da rede, sua implementação em arquitetura maciçamente paralela e resultados de simulações e implementações em FPGA são apresentados, demonstrando a viabilidade da proposta. Como resultado final, conclui-se que a rede ONNIS-GI apresenta maior apelo de ordem prática, sendo uma abordagem inovadora e promissora na solução de problemas de segmentação de imagens, possuindo capacidade para separar temporalmente os segmentos encontrados e facilitando a posterior identificação dos mesmos. Sob o ponto de vista prático, a nova rede pode ser utilizada para implementar chips de visão digitais com arquitetura maciçamente paralela, explorando a velocidade de tais topologias e apresentando também flexibilidade para implementação de procedimentos de segmentação de imagens mais sofisticados.

Modelo de transmissão da dengue com competição larval uniforme

A dengue representa um sério problema de saúde pública no Brasil, que apresenta condições climáticas favoráveis ao desenvolvimento e proliferação do Aedes aegypti, vetor transmissor da doeça. O mosquito Aedes aegypti passa por diferentes estágios de desenvolvi- mento com características distintas, logo, para uma descrição mais aproximada da história de vida desta espécie e, consequentemente, do comportamento da epidemia, considera-se neste trabalho um modelo com distribuição etária, fundamental para a determinação das propriedades de estabilidade de populações que têm fases distintas de desenvolvimento. O modelo com distribuição etária para a população de mosquitos é descrito por um conjunto de equações diferenciais ordinárias com retardo de difícil análise e implementação. Inicialmente, apresenta-se o modelo SEIR com dinâmica vital, onde a população de mosquitos estabiliza rapidamente e, após, incorpora-se a ele um modelo com competição larval uniforme para população de insetos, com distribuição etária, o que provoca um período de instabilidade nas populações de larvas e adultos, estágios de desenvolvimento considerados para o vetor. A análise realizada investiga as consequências que este período de instabilidade provoca na evolução da epidemia.

Decomposição de Calderon e suas aplicações na teoria da regularidade em equações elípticas

Este trabalho tem por objetivo estudar a regularidade de soluções de Equações Diferenciais Parciais Elípticas da forma Lu = f, para f 2 Lp(­), onde p > 1. Para isto, usamos a Decomposição de Calderon-Zygmund e um resultado que é consequência deste, o Teorema da Interpolação de Marcinkiewicz. Além disso, usando quocientes-diferença provamos a regularidade das soluções para o caso p = 2 e L = ¡¢ de uma forma alternativa.

Solução da equação de difusão unidimensional transiente para o estudo da dispersão de poluentes na camada limite planetária

Neste trabalho apresenta-se uma solução analítica para a dispersão vertical turbulenta em uma Camada Limite Convectiva e em uma Camada Limite Estável. A equação analisada considera a difusão com velocidades finitas, o que representa o transporte turbulento fisicamente correto. Considerando o caráter não-local, adicionam-se na equação que representa uma fonte área instantânea, termos como: o tempo de relaxação, a assimetria, a escala de tempo Lagrangeana e a velocidade turbulenta vertical. A solução é obtida utilizando-se a técnica da Transformada de Laplace. Os parâmetros que encerram a turbulência são derivados da teoria de difusão estatística de Taylor combinada com a teoria de similaridade. Foram utilizados coeficientes de difusão especáficos para cada uma das camadas. A transformada inversa é obtida através do esquema numérico de quadratura Gaussiana. São apresentadas várias simulações para diferentes alturas de fonte área e obtém-se o valor da concentração para alturas próximas ao solo e próximas ao topo da Camada Limite Planetária. A inserção do termo de contra-gradiente na equação resultou em uma pequena influência na concentração de poluentes, observada de forma mais expressiva na Camada Limite Convectiva.

Otimização topológica de estruturas contínuas submetidas a restrições de flexibilidade, volume e frequência natural

Este trabalho trata dos problemas de otimização de minimização de volume com restrição de flexibilidade e freqüência natural e minimização de flexibilidade com restrição de volume. Os problemas são resolvidos para estruturas bidimensionais e tridimensionais. As equações diferenciais de equilíbrio são solucionadas de forma aproximada através do método dos elementos finitos, em um subespaço de dimensão finita. O método utilizado no estudo é o da otimização topológica, o qual consiste em encontrar dentro de um domínio pré-existente uma distribuição ideal de material. São avaliadas técnicas como programação linear e critério de ótimo. Em ambos os casos são utilizadas sensibilidades calculadas analiticamente. Para a otimização com restrição modal, problemas característicos como autovalores repetidos e normalização do autovetor são tratados. Ferramentas usadas na otimização topológica, como método da continuação, penalização e filtragem são discutidos. São abordados também problemas e características inerentes ao processo de otimização topológica, tais como instabilidades de tabuleiros, dependência de malha e sensibilidade da topologia a diferentes condições de contorno. Os resultados obtidos permitem avaliações referentes à otimização topológica (geometrias, ou seja, topologias resultantes) sob diferentes condições, utilizando-se as ferramentas discutidas nesse trabalho.

Avaliação de projetos na presença de volatilidade: estudo de caso de empresas no setor de telecomunicações

Quando as empresas decidem se devem ou não investir em determinado projeto de investimentos a longo prazo (horizonte de 5 a 10 anos), algumas metodologias alternativas ao Fluxo de Caixa Descontado (FCD) podem se tornar úteis tanto para confirmar a viabilidade do negócio como para indicar o melhor momento para iniciar o Empreendimento. As análises que levam em conta a incerteza dos fluxos de caixa futuros e flexibilidade na data de início do projeto podem ser construídos com a abordagem estocástica, usando metodologias como a solução de equações diferenciais que descrevem o movimento browniano. Sob determinadas condições, as oportunidades de investimentos em projetos podem ser tratados como se fossem opções reais de compra, sem data de vencimento, como no modelo proposto por McDonald-Siegel (1986), para a tomada de decisões e momento ótimo para o investimento. Este trabalho analisa a viabilidade de investimentos no mercado de telecomunicações usando modelos não determinísticos, onde a variável mais relevante é a dispersão dos retornos, ou seja, que a variância representa o risco associado a determinado empreendimento.

Uma solução híbrida para a equação de Fokker-Planck dependente da energia usando a técnica da transformada de Laplace e diferenças finitas

Neste trabalho é obtida uma solução híbrida para a equação de Fokker-Planck dependente da energia, muito utilizada em problemas de implantação iônica. A idéia consiste na aplicação da transformada de Laplace na variável de energia e aplicação de um esquema de diferenças finitas nas variáveis espacial e angular desta equação. Tal procedimento gera um problema matricial simbólico para a energia transformada. Para resolver este sistema, procede-se a inversão de Laplace da matriz (sI+A), onde s é um parâmetro complexo, I a matriz identidade e A uma matriz quadrada gerada pela discretização das variáveis espacial e angular. A matriz A não é diagonalizável, desta forma, contorna-se este problema decompondo esta matriz na soma de outras duas, onde uma delas é diagonalizável. É gerado então um método iterativo de inversão, semelhante ao método da fonte fixa associado ao método de diagonalização, do qual o resultado fornecido são os valores para o fluxo de partículas do sistema. A partir disto pode-se determinar a energia depositada no sistema eletrônico e nuclear do alvo. Para validar os resultados obtidos faz-se a simulação de implantação de íons de B em Si numa faixa energética de 1keV a 50MeV, comparam-se os resultados com simulação gerada numericamente pelo software SRIM2003.

Rafaella : um esquema para um novo paradigma de criptografia de chave pública

O presente trabalho apresenta um novo esquema de criptografia de chave pública baseado no emprego de funções para representar as mensagens original e cifrada. No esquema proposto – denominado Rafaella -, o processo de cifração consiste na aplicação de um deslocamento no argumento da função que representa a mensagem, de modo que se f(x) descreve a mensagem original, então f(x+z) representa a respectiva mensagem cifrada. O deslocamento z representa um número complexo que, no esquema proposto, representa a forma das chaves privadas dos participantes. A dificuldade da resolução do problema inversos concentra-se na obtenção das partes real e imaginária do deslocamento z, que pode ser efetuada através de método de força bruta, ou da resolução de um problema de contorno. A segunda alternativa envolve a resolução de equações diferenciais. Dentre os métodos disponíveis para a resolução de equações diferenciais, o emprego dos chamados grupos de Lie constitui, via de regra, a estratégia mais apropriada para a obtenção de soluções analíticas, que demandam menor tempo de processamento do que as formulações numéricas. Mesmo assim, a solução obtida através da utilização dos grupos de Lie requer elevado número de operações simbólicas.

Uma resenha sobre modelos de apreçamento de derivativos

Apresento aqui uma abordagem que unifica a literatura sobre os vários modelos de apreçamento de derivativos que consiste em obter por argumentos intuitivos de não arbitragem uma Equação Diferencial Parcial(EDP) e através do método de Feynman-Kac uma solução que é representada por uma esperança condicional de um processo markoviano do preço do derivativo descontado pela taxa livre de risco. Por este resultado, temos que a esperança deve ser tomada com relação a processos que crescem à taxa livre de risco e por este motivo dizemos que a esperança é tomada em um mundo neutro ao risco(ou medida neutra ao risco). Apresento ainda como realizar uma mudança de medida pertinente que conecta o mundo real ao mundo neutro ao risco e que o elemento chave para essa mudança de medida é o preço de mercado dos fatores de risco. No caso de mercado completo o preço de mercado do fator de risco é único e no caso de mercados incompletos existe uma variedade de preços aceitáveis para os fatores de risco pelo argumento de não arbitragem. Neste último caso, os preços de mercado são geralmente escolhidos de forma a calibrar o modelo com os dados de mercado.