Modelagem e otimização de experimentos para o tratamento térmico de recozimento: um estudo com o algoritmo simulated annealing

Esta dissertação apresenta um estudo da modelagem de experimentos aplicados a um processo industrial de tratamento térmico. A motivação deste trabalho surgiu diante das dificuldades associadas aos processos de recozimento industrial de aços do tipo baixa liga, na tentativa de encontrar temperaturas nas quais as durezas superficiais dos aços atingissem valores suficientemente baixos, adequados para etapas posteriores de fabricação, em especial a usinagem. Inicialmente forem realizados diversos experimentos com diferentes aços, onde a dureza superficial é obtida em função da temperatura de recozimento e dos teores de carbono e silício das amostras utilizadas. Em seguida propôs-se um modelo quadrático para modelar a dureza superficial como função dessas três variáveis. A estimação de parâmetros do modelo proposto foi realizada com o emprego do algoritmo Simulated Annealing, uma meta-heurística para otimização global que procura imitar o processo de recozimento de um material sólido. Finalmente, usando-se o modelo proposto, foi resolvido o chamado problema inverso, o qual consiste na estimação da temperatura de recozimento em função dos teores de carbono e silício e da dureza desejada.

Modelagem e simulação da propagação de ondas em barras não homogêneas envolvendo materiais elásticos não lineares.

O objetivo deste trabalho é tratar da simulação do fenômeno de propagação de ondas em uma haste heterogênea elástico, composta por dois materiais distintos (um linear e um não-linear), cada um deles com a sua própria velocidade de propagação da onda. Na interface entre estes materiais existe uma descontinuidade, um choque estacionário, devido ao salto das propriedades físicas. Empregando uma abordagem na configuração de referência, um sistema não-linear hiperbólico de equações diferenciais parciais, cujas incógnitas são a velocidade e a deformação, descrevendo a resposta dinâmica da haste heterogénea. A solução analítica completa do problema de Riemann associado são apresentados e discutidos.

Prismas do exílio : trajetória intelectual e modelagem do self em Anatol Rosenfeld e Otto maria Carpeaux

Esta tese tem por objetivo investigar as trajetórias intelectuais de Anatol Rosenfeld e Otto Maria Carpeaux. Exilados no Brasil durante a década de 30 do século passado, em função do avanço do nazifascismo na Europa, o alemão Rosenfeld e o austríaco Carpeaux souberam reinventar suas existências na nova terra, transformando-se em dois dos maiores críticos literários do país. Além de examinar os itinerários percorridos por ambos os autores, esta tese terá como foco analisar a relação entre exílio e atividade intelectual, promovendo uma reflexão a respeito do lugar das cidades de Berlim e Viena na construção do arcabouço espiritual de Rosenfeld e Carpeaux. Este estudo também tratará das suas primeiras produções em solo brasileiro, bem como da recepção que obtiveram de seus anfitriões, delineando, assim, os contornos de uma complexa atmosfera intelectual marcada por uma espécie de cordialidade literária. Como último objetivo, esta tese debaterá o processo de modelagem do self por meio do qual ambos os críticos conseguiram aprimorar a própria individualidade, tornando-se, enfim, figuras cuja memória deve ser recuperada.

Modelagem computacional aplicada na estimação de parâmetros cinéticos em processos químicos e cromatográficos

A estimação de parâmetros cinéticos em processos químicos e cromatográficos utilizando técnicas numéricas assistidas por computadores tem conduzido para melhoria da eficiência e o favorecimento da compreensão das fenomenologias envolvidas nos mesmos. Na primeira parte deste trabalho será realizada a modelagem computacional do processo de produção de biodiesel via esterificação, sendo que, o método de otimização estocástica Random Restricted Window (R2W) será correlacionado com os dados experimentais da produção de biodiesel a partir da esterificação do ácido láurico com etanol anidro na presença do catalisador ácido nióbico (Nb2O5). Na segunda parte do mesmo será realizada a modelagem computacional do processo de cromatografia de adsorção (batch process) onde serão correlacionados os dados provenientes dos modelos cinéticos de HASHIM, CHASE e IKM2 com os dados experimentais da adsorção de amoxicilina com quitosana, e também serão correlacionados os dados experimentais da adsorção de Bovine Serum Albumin (BSA) com Streamline DEAE com os dados provenientes de uma nova aplicação do método R2W mediante a implementação de um modelo cinético reversível. Ademais, as constantes cinéticas para cada processo supracitado serão estimadas levando em consideração o valor mínimo da função resíduos quadrados.

Modelagem de bioinvasão do coral-sol (Tubastraea coccinea e T. tagusensis):mecanismos da ocupação e dispersão e identificação de sua potencial distribuição geográfica

Os fatores que explicam a distribuição observada em plantas e animais é uma pergunta que intriga naturalistas, biogeógrafos e ecólogos há mais de um século. Ainda nos primórdios da disciplina de ecologia, as tolerâncias ambientais já haviam sido apontadas como as grandes responsáveis pelo padrão observado da distribuição dos seres vivos, o que mais tarde levou à concepção de nicho ecológico das espécies. Nos últimos anos, o estudo das distribuições dos organismos ganhou grande impulso e destaque na literatura. O motivo foi a maior disponibilidade de catálogos de presença de espécies, o desenvolvimento de bancos de variáveis ambientais de todo o planeta e de ferramentas computacionais capazes de projetar mapas de distribuição potencial de um dado organismo. Estes instrumentos, coletivamente chamados de Modelos de Distribuição de Espécies (MDEs) têm sido desde então amplamente utilizados em estudos de diferentes escopos. Um deles é a avaliação de potenciais áreas suscetíveis à invasão de organismos exóticos. Este estudo tem, portanto, o objetivo de compreender, através de MDEs, os fatores subjacentes à distribuição de duas espécies de corais escleractíneos invasores nativos do Oceano Pacífico e ambas invasoras bem sucedidas de diversas partes do Oceano Atlântico, destacadamente o litoral fluminense. Os resultados mostraram que os modelos preditivos da espécie Tubastraea coccinea (LESSON, 1829), cosmopolita amplamente difundida na sua região nativa pelo Indo- Pacífico demonstraram de maneira satisfatória suas áreas de distribuição nas áreas invadidas do Atlântico. Sua distribuição está basicamente associada a regiões com alta disponibilidade de calcita e baixa produtividade fitoplanctônica. Por outro lado, a aplicação de MDEs foi incapaz de predizer a distribuição de T. tagusensis (WELLS,1982) no Atlântico. Essta espécie, ao contrário de sua congênere, tem distribuição bastante restrita em sua região nativa, o arquipélago de Galápagos. Através de análises posteriores foi possível constatar a mudança no nicho observado durante o processo de invasão. Finalmente, o sucesso preditivo para T. coccinea e o fracasso dos modelos para T. tagusensis levantam importantes questões sobre quais os aspectos ecológicos das espécies são mais favoráveis à aplicação de MDEs. Adicionalmente, lança importantes ressalvas na utilização recentemente tão difundida destas ferramentas como forma de previsão de invasões biológicas e em estudos de efeitos de alterações climáticas sobre a distribuição das espécies.

Modelagem e simulação do escoamento imiscível em meios porosos fractais descritos pela equação de Kozeny-Carman Generalizada

O presente trabalho trata do escoamento bifásico em meios porosos heterogêneos de natureza fractal, onde os fluidos são considerados imiscíveis. Os meios porosos são modelados pela equação de Kozeny-Carman Generalizada (KCG), a qual relaciona a porosidade com a permeabilidade do meio através de uma nova lei de potência. Esta equação proposta por nós é capaz de generalizar diferentes modelos existentes na literatura e, portanto, é de uso mais geral. O simulador numérico desenvolvido aqui emprega métodos de diferenças finitas. A evolução temporal é baseada em um esquema de separação de operadores que segue a estratégia clássica chamada de IMPES. Assim, o campo de pressão é calculado implicitamente, enquanto que a equação da saturação da fase molhante é resolvida explicitamente em cada nível de tempo. O método de otimização denominado de DFSANE é utilizado para resolver a equação da pressão. Enfatizamos que o DFSANE nunca foi usado antes no contexto de simulação de reservatórios. Portanto, o seu uso aqui é sem precedentes. Para minimizar difusões numéricas, a equação da saturação é discretizada por um esquema do tipo "upwind", comumente empregado em simuladores numéricos para a recuperação de petróleo, o qual é resolvido explicitamente pelo método Runge-Kutta de quarta ordem. Os resultados das simulações são bastante satisfatórios. De fato, tais resultados mostram que o modelo KCG é capaz de gerar meios porosos heterogêneos, cujas características permitem a captura de fenômenos físicos que, geralmente, são de difícil acesso para muitos simuladores em diferenças finitas clássicas, como o chamado fenômeno de dedilhamento, que ocorre quando a razão de mobilidade (entre as fases fluidas) assume valores adversos. Em todas as simulações apresentadas aqui, consideramos que o problema imiscível é bidimensional, sendo, portanto, o meio poroso caracterizado por campos de permeabilidade e de porosidade definidos em regiões Euclideanas. No entanto, a teoria abordada neste trabalho não impõe restrições para sua aplicação aos problemas tridimensionais.

Análise da resposta dinâmica experimental de uma passarela tubular mista, aço-concreto, submetida ao caminhar humano.

Este trabalho de pesquisa apresenta como objetivo principal o desenvolvimento de investigação experimental dinâmica sobre estrutura real de uma passarela tubular mista aço-concreto. O sistema estrutural objeto deste trabalho corresponde a uma passarela composta por três vãos (32,5m, 17,5m e 20,0m, respectivamente) e dois balanços (7,50m e 5,0m, respectivamente), com comprimento total de 82,5m. A passarela com estrutura contínua de aço com as ligações soldadas se apoia em quatro pórticos também de aço. Estruturalmente está constituída por duas treliças planas que se interligam através de contraventamentos horizontais fixados na corda superior e inferior da treliça e lajes de concreto, formando um sistema misto com interação completa. A estrutura está submetida correntemente à travessia de pedestres e ciclistas. Testes experimentais foram realizados sobre o sistema estrutural e confrontados com resultados numéricos. Para a modelagem numérica do sistema são empregadas técnicas usuais de discretização, via método dos elementos finitos (MEF), por meio do programa ANSYS. Os resultados experimentais são analisados de acordo com a metodologia desenvolvida, sendo realizada análise modal experimental para a determinação das propriedades dinâmicas: freqüências, modos e taxa de amortecimento, enquanto que os resultados da estrutura, em termos de aceleração de pico, são comparados com os valores limites propostos por diversos autores, normas e recomendações de projeto, para uma avaliação do desempenho da estrutura em relação a vibração quando solicitada pelo caminhar dos pedestres no que diz respeito a critério para conforto humano.

Avaliação de ligações K e T entre perfis estruturais tubulares circulares.

O difundido uso de perfis estruturais tubulares, principalmente devido às vantagens associadas ao comportamento estrutural e estético, levou a uma intensa utilização nos países da Europa, Sudeste Asiático, América do Norte e na Austrália. Países como Canadá, Inglaterra, Alemanha e Holanda fazem uso intensivo dessas estruturas e contam com uma produção corrente, industrializada e contínua com alto grau de desenvolvimento tecnológico. O aumento da oferta deste tipo de perfis, aliado a recentes pesquisas nesta área, leva o Brasil a se inserir neste cenário mundial. Entretanto, as normas brasileiras que regem o dimensionamento destes perfis ainda não se encontram atualizadas, principalmente no que tange ao dimensionamento de ligações de perfis tubulares. Considerando esta perspectiva, este trabalho apresenta uma análise de ligações tipo K e T com perfis tubulares circulares (CHS). Propõe-se estabelecer um quadro comparativo entre as formulações analíticas de dimensionamento proposta pelo Eurocode 3 Parte 1.8, 2 edição do guia de projeto de ligações tubulares do CIDECT, projeto de norma brasileira PN 02.125.03-004 e critérios de deformação limite. Para cada um dos tipos de ligações analisadas, desenvolveu-se um modelo em elementos finitos no programa Ansys, calibrado e validado com resultados experimentais e numéricos existentes na literatura. As não-linearidades física e geométrica foram incorporadas aos modelos, a fim de se mobilizar totalmente a capacidade resistente da ligação. A não-linearidade do material foi considerada através do critério de plastificação de Von Mises através de uma lei constitutiva tensão versus deformação bilinear. A não-linearidade geométrica foi introduzida no modelo através da Formulação de Lagrange Atualizado considerando-se a previsão de grandes deformações de forma a permitir a redistribuição de carregamento na ligação após o escoamento inicial. Esta dissertação propõe ainda, a modelagem de uma treliça planar constituída de perfis tubulares, objetivando estabelecer uma comparação entre um nó isolado e um nó da treliça planar.

Avaliação do comportamento estrutural de subestações de energia elétrica com o uso do aço inoxidável.

A crescente utilização do aço inoxidável como elemento estrutural despertou o interesse de clientes, arquitetos e engenheiros nos últimos anos. Apesar do custo ainda elevado, a sua aplicação na construção civil vem substituindo outros elementos estruturais. Seja por sua alta resistência à corrosão, aumentando a relação custo benefício; sua estética, proporcionando formas cada vez mais ousadas ou; seu apelo ambiental, gerando menos resíduos no meio ambiente. As subestações representam um papel importante no fornecimento de energia. Como possuem grande complexidade para manutenção, foi escolhida a estrutura suporte de seu barramento, para o dimensionamento em aço inoxidável. Desta forma, minimizando as paradas para realização de manutenções das estruturas, possibilitando maior qualidade no fornecimento de energia elétrica. Para fins comparativos foi escolhido o projeto de uma SE existente, cuja estrutura de suporte do barramento, foi construída por treliças formadas por cantoneiras de aço carbono galvanizado. Inicialmente, o dimensionamento foi desenvolvido utilizando perfis H e I funcionando como viga-coluna para os dois tipos de aço. Num segundo momento, a estrutura foi dimensionada como treliças planas. Todos os dimensionamentos foram realizados de acordo com as prescrições normativas do EUROCODE 3. Após realização dos dimensionamentos, foram apresentadas as análises comparativas dos custos envolvidos para os tipos de aço. Abordando o investimento inicial, os gastos com manutenção ao longo da vida e os custos elétricos agregados à redução das paradas para manutenção.

Análise numérica do comportamento de juntas entre aduelas de vigas protendidas.

Diversos pesquisadores têm estudado o comportamento e o emprego de aduelas de concreto, que constituem as vigas segmentadas em sistemas estruturais, de maneira especial pontes e viadutos. Por esta razão, inúmeros trabalhos têm sido publicados nos últimos anos respaldados por testes experimentais e análises numéricas. O comportamento destas vigas contrasta com as clássicas vigas monolíticas em diversos aspectos, pois, a estrutura é composta de partes de elementos de concreto pré-moldado que, após serem posicionados no local definitivo, são protendidos. A protensão pode ser aderente ou não aderente. A principal vantagem deste sistema de construção é a rapidez e o alto controle de qualidade, por isso é largamente utilizado, havendo uma demanda de estudo de previsão do seu real comportamento No presente trabalho apresenta-se uma modelagem numérica via elementos finitos, para simular o comportamento de vigas compostas por aduelas justapostas sem material ligante entre as juntas. A protensão aplicada é aderente e a análise considera a não linearidade da região da junta. Assim sendo, o objetivo desta investigação é dar uma contribuição ao estudo do comportamento estrutural estático de vigas segmentadas, atentando para o comportamento das juntas, utilizando um programa comercial. Para o modelo são empregadas técnicas usuais de discretização, via método dos elementos finitos (MEF), por meio do programa de elementos finitos SAP2000[93]. O modelo proposto é constituído de elementos de placa próprios para concreto para representar a viga, a protensão é introduzida por meio de barras bidimensionais que transferem as tensões ao longo de seu comprimento e as juntas são implementadas utilizando elementos de contato. A analise é bidimensional e considera os efeitos das perdas de protensão. Este trabalho de pesquisa objetiva também o estudo de elementos de contato especialmente as características de deformação para esta ferramenta computacional. A definição dos parâmetros para o modelo foi feita com base em dados experimentais disponíveis na literatura. O modelo numérico foi calibrado e confrontado com resultados experimentais obtidos em ensaios de laboratório.

Análise estática e dinâmica de lajes nervuradas de concreto armado.

Nesta dissertação desenvolve-se um estudo sobre o comportamento estrutural de lajes nervuradas de concreto armado, tendo em mente o emprego destas em sistemas estruturais na engenharia civil, de forma viável sob o ponto de vista técnico e econômico. No cenário atual, diversos tipos de sistemas estruturais têm sido desenvolvidos, dentre estes as lajes nervuradas de concreto armado com base no sistema Reduzcon. Essas lajes nervuradas têm como principais características o baixo consumo de material e o peso próprio bastante reduzido. Este sistema utiliza cubas cilíndricas invertidas metálicas denominadas BRC (Barrote Redutor de Concreto), que suportam o concreto da laje. Assim sendo, este trabalho de pesquisa apresenta um estudo acerca da resposta estática e dinâmica do sistema estrutural em questão, em relação aos deslocamentos translacionais verticais e tensões; e, bem como, da resposta dinâmica do sistema sob o ponto de vista de conforto humano. Ênfase especial é dada às atividades humanas rítmicas, com base no emprego de modelos de carregamento que representem as ações inerentes à ginástica aeróbica. A definição das ações dinâmicas atuantes sobre os modelos estruturais foi feita com base em resultados de testes experimentais que levam em conta os efeitos de multidão. A análise numérica fundamenta-se na modelagem computacional do sistema, através do Método dos Elementos Finitos (MEF), de forma a representar apropriadamente o comportamento estrutural das lajes Reduzcon. São empregadas técnicas usuais de discretização, por meio do programa ANSYS. As vigas de bordo, nervuras e a laje de concreto são simuladas por elementos finitos sólidos tridimensionais.

Estudo do conforto humano em pisos mistos (aço-concreto) submetidos a ações humanas rítmicas.

Este trabalho de pesquisa objetiva o estudo do comportamento dinâmico de pisos mistos (aço-concreto), em edificações de andares múltiplos, sob o ponto de vista de conforto humano, quando essas estruturas encontram-se submetidas às atividades rítmicas provenientes dos seres humanos. A definição das ações dinâmicas atuantes sobre os modelos estruturais foi feita com base em resultados experimentais, com os indivíduos praticando atividades rítmicas e não rítmicas associadas à ginástica aeróbica e saltos à vontade. Os modelos estruturais investigados baseiam-se em edificações mistas de andares múltiplos. O sistema estrutural é do tipo misto (aço-concreto), composto por vigas de aço em seção do tipo I e laje de concreto armado. A análise fundamenta-se na modelagem computacional dos sistemas estruturais, através do Método dos Elementos Finitos (MEF). São empregadas técnicas usuais de discretização, por meio do emprego do programa ANSYS. Uma análise paramétrica foi desenvolvida sobre três modelos estruturais, com dois, três e quatro pavimentos. Os valores das acelerações máximas encontradas na análise são confrontados e comparados com os limites propostos por recomendações internacionais. Os resultados obtidos mostram que os limites recomendados em diversas normas de projeto foram ultrapassados. Esses resultados demonstram que atividades rítmicas oriundas dos seres humanos podem gerar acelerações de pico elevadas, violando critérios de projeto, no que concerne ao conforto humano. Foi observado também que estas ações dinâmicas podem comprometer o conforto humano em pisos adjacentes, próximos do local onde a carga dinâmica está sendo efetivamente aplicada.

Modelagem numérica de ligações viga-coluna em aço sob momento fletor e força normal.

As ligações desempenham um papel fundamental no comportamento global das estruturas de aço. Inúmeros trabalhos de pesquisa têm sido desenvolvidos para entender o comportamento real de uma ligação e sua influência na resistência global dos pórticos. Atualmente, a Norma Brasileira de estruturas de aço de edificações, NBR 8800, considera o comportamento das ligações entre duas situações extremas: rígidas, onde não ocorre nenhuma rotação entre os membros conectados, transferindo momento fletor, força cortante e força normal; ou flexíveis, caracterizadas pela liberdade de rotação entre os membros conectados, impedindo a transmissão de momento fletor. Outras normas de projeto de estruturas de aço, consideram que as ligações apresentam um comportamento intermediário, ou seja, semi-rigídas, que podem estar submetidas a uma combinação de momento fletor e esforço normal. Porém, mesmo com a combinação, estas normas não consideram a presença de esforço normal (tração e/ou compressão). Uma limitação empírica de 5% da resistência plástica da viga é a única condição imposta no Eurocode 3. Para o estudo da ligação semi-rigída será utilizada a filosofia do Método das Componentes, que verifica a resistência da ligação e sua classificação quanto à rigidez rotacional, desenvolvida através de modelos mecânicos (modelos de molas). O objetivo deste trabalho é descrever alguns resultados de caracterização de ligações viga-coluna com placa de extremidade ajustada à altura da viga obtidos através de um modelo de elementos finitos. Para tal, será realizada uma análise não-linear geométrica e de material. Esta análise possibilitará avaliar os principais parâmetros que influenciam no comportamento deste componente no que diz respeito a sua avaliação em termos de distribuição de tensões e deformações no modelo de forma global.

Comportamento Estrutural de Elementos em Aço Inoxidável

Atualmente, a utilização do aço inoxidável em elementos estruturais ainda é por muitos engenheiros e arquitetos, considerada uma solução extravagante para os problemas da engenharia. Todavia, mudanças de atitudes dentro da construção civil e uma transição global para um desenvolvimento sustentável e redução em impactos ambientais tem seguramente provocado um aumento no uso do aço inoxidável. A maioria das normas de projeto de aço inoxidável atuais ainda são baseadas em analogias assumidas com o comportamento de estruturas de aço carbono. Todavia, o aço inoxidável apresenta quatro curvas tensão versus deformação não-lineares sem patamar de escoamento e região de encruamento claramente definidos (tração e compressão, paralela e perpendicular a direção de laminação), modificando assim, seu o comportamento global. Na presente investigação foram utilizados o aço inoxidável austenítico 304, e o aço carbono USI 300, ambos com tensão de escoamento similares, mas com tensões últimas e ductilidades distintos. Em elementos estruturais submetidos a tensões normais de tração, usualmente a ruptura da seção líquida representa um dos estados limites últimos a serem verificados. Com o objetivo de se avaliar a resistência a tração de elementos estruturais aparafusados em aço inoxidável, este trabalho executou um programa experimental inovador envolvendo ligações aparafusadas defasadas sob tração. O programa experimental foi executado em peças de aço carbono e aço inoxidável de forma a comparar as principais semelhanças e diferenças entre estes dois tipos de aços estruturais. O programa experimental possibilitou que conclusões significativas no comportamento a tração destas ligações pudessem ser observadas. Dentre outras variáveis que controlam estes estados limites últimos foi verificada a influência significativa da espessura da placa de aplicação do carregamento, disposição, configuração e número de parafusos da ligação e propriedades do aço inoxidável como ductilidade e razão entre suas tensões de escoamento e de ruptura.

Modelagem matemática da transferência de calor numa placa plana sob o efeito de uma fonte pontual externa de radiação térmica.

Este trabalho apresenta uma modelagem matemática para o processo de aquecimento de um corpo exposto a uma fonte pontual de radiação térmica. O resultado original que permite a solução exata de uma equação diferencial parcial não linear a partir de uma seqüência de problemas lineares também é apresentado. Gráficos gerados com resultados obtidos pelo método de diferenças finitas ilustram a solução do problema proposto.