Vibrações naturais em um sistema de interação viga-água sem o efeito de onda de superfície livre

Neste trabalho é apresentado um estudo sobre o comportamento dinâmico de sistemas flexíveis com interação viga-água. Especificamente, é escolhida uma barragem freqüentemente encontrada em problemas de engenharia de irrigação e represas. O modelo matemático utilizado para descrever tal fenômeno tem como variáveis principais a pressão hidrodinâmica e os deslocamentos horizontais da barragem. Uma vez formuladas as equações governantes de tal sistema, que resultam em uma equação de onda para a pressão e uma equação de Euler-Bernoulli para os deslocamentos, é utilizada a técnica de separação de variáveis para proceder à sua solução. Por simplicidade, assumiu-se que no modelo não havia o efeito de onda da superfície livre. Primeiro, foi resolvida a equação para a pressão, sendo calculadas as suas componentes temporais e espaciais. Depois, foi resolvida a equação de Euler-Bernoulli junto com a condição de interface. Foram simulações para o material específico, e apresentados os resultados conseguidos. No final do trabalho são enunciadas as conclusões pertinentes.

Análise de estruturas de concreto armado sujeitas à corrosão de armadura por cloretos através do método dos elementos finitos

A corrosão da annadura é um dos principais problemas para um grande número de es1ruturas de concreto annado por causa da redução substancial da área de aço, da perda de aderência, da fonnação de fissuras e lascamento do concreto. O modelo tridimensional proposto por este trabalho é baseado no processo de corrosão iniciado por ataque de cloretos e propagado a partir de micropilhas de corrosão, utilizando o método dos elementos finitos. A influência da temperatura, da umidade relativa do concreto, da ligação dos cloretos com compostos do concreto, da fonnação de fissuras provenientes do carregamento e da resistência ôhmica do concreto são levados em consideração. O modelo é dividido em 4 módulos que interagem entre si, sendo eles: análise estrutural, transferência de calor, difusão de cloretos e corrosão eletroquímica. O modelo computacional para a análise estrutural utiliza dois procedimentos. No primeiro busca-se a resposta da estrutura para um carregamento instantâneo, considerando-se o material com comportamento elasto-viscoplástico, associado com modelos de annadura incorporada e fissuras distribuídas.No segundo, detennina-se a resposta da estrutura ao longo do tempo com a teoria da solidificação proposta por Bazant e Prasannan (1989), utilizando as cadeias Kelvin, apresentando ótimos resultados com uma grande economia computacional. O processo de transferência de calor é baseado na fonnulação tridimensional da equação de Helmholtz, com condições de contorno de fluxos provenientes da radiação solar, convecção e radiação, apresentando bons resultados em relação a fonnulações analíticas. Para análise de vigas com a variação da temperatura média diária, constatou-se que as mesmas se encontram em equilíbrio ténnico com o ambiente A análise de difusão de cloretos implementada apresentou um comportamento idêntico ao da segunda lei de Fick para casos unidimensionais. Para análises de vigas carregadas se constatou que as fissuras provenientes de cargas de serviço não influenciam o processo de difusão e que o coeficiente de difusão referencial é o principal parâmetro na análise. O modelo de corrosão eletroquímica é baseado na fonnulação proposta por Liu (1996), obtendo-se como resultado a taxa de corrosão em função do tempo.

Análise experimental e numérica da solidificação direcional de materiais de mudança de fase

O presente trabalho apresenta uma análise da solidificação unidirecional a partir de um estudo experimental e numérico. A fim de possibilitar a visualização do processo são evitadas as altas temperaturas. Desta forma, materiais com baixo ponto de fusão, abaixo de 50°C, translúcidos e transparentes, são fundidos e solidificados em um solidificador direcional de forma cilíndrica. Visando um processo de solidificação controlado por difusão pura, este solidificador é resfriado pela sua base. São obtidos no experimento os campos de temperatura, além da posição da interface sólido-líquido ao longo do tempo. Na simulação numérica do processo é resolvida a equação da energia, formulada em termos da entalpia usando uma malha fixa. Os resultados são comparados e analisados de modo a validar a simulação, que passa assim a servir para a investigação de outras situações de interesse no âmbito da fusão-solidificação.

Cálculo da complexidade exata de algoritmos do tipo divisão-e-conquista através das equações características

A equação de complexidade de um algoritmo pode ser expressa em termos de uma equação de recorrência. A partir destas equações obtém-se uma expressão assintótica para a complexidade, provada por indução. Neste trabalho, propõem-se um esquema de solução de equações de recorrência usando equações características que são resolvidas através de um "software" de computação simbólica, resultando em uma expressão algébrica exata para a complexidade. O objetivo é obter uma forma geral de calcular a complexidade de um algoritmo desenvolvido pelo método Divisão-e-Conquista.

Impactos comerciais da Área de Livre Comércio das Américas : uma aplicação do modelo gravitacional

Este trabalho analisa os efeitos que a formação da Área de Livre Comércio das Américas (ALCA) poderá exercer sobre o fluxo de comércio de um grupo de países que farão parte deste bloco, a partir da utilização de um modelo gravitacional. Inicialmente, apresentamos a fundamentação teórica sobre a formação de blocos econômicos, enfatizando os principais conceitos como criação de desvio de comércio. Um breve relato sobre o andamento das negociações para a formação da ALCA é apresentado a seguir, com ênfase no seu desenvolvimento histórico e temas ainda pendentes. A terceira parte do trabalho destina-se comentar alguns estudos sobre o mesmo tema que utilizaram metodologias diferentes para analisar o mesmo tema. Na quarta parte, são estimados os efeitos sobre o comércio de três importantes blocos já em vigor na região (Pacto Andino, MERCOSUL e NAFTA) através do modelo gravitacional. Este modelo permite identificar os efeitos sobre o comércio que podem ser atribuídos exclusivamente a formação dos blocos, separando-os dos demais fatores que influenciam o padrão de comércio. Os resultados obtidos por meio da equação gravitacional servem de base para a criação de cenários prospectivos da ALCA. No cenário mais otimista, o comércio do Brasil com os principais parceiros da ALCA cresceria significativamente, assim como os fluxos comerciais entre os demais países analisados.

Simulação numérica de um sistema de armazenamento de calor em estufas de plasticultura

Apresentar um modelo para simular um sistema de armazenamento de calor no solo em estufas para plasticultura é o objetivo do presente trabalho. O sistema consiste num feixe de tubos enterrados no solo. A convecção forçada de ar no seu interior realiza a troca térmica necessária para manter as estufas sob faixas desejadas de temperatura. O objetivo do modelo é investigar os efeitos no calor armazenado e a influência das variáveis, tais como diâmetro, comprimento, espaçamento entre os tubos e a velocidade de ar no canal provocam no sistema. O solo é tratado como um meio difusivo e avalia-se a contribuição do termo de condensação e evaporação da água contida no ar em escoamento nos tubos. A equação da energia é resolvida para o solo e para o ar. Os tubos de seção transversal circular são modelados como tubos de seção transversal quadrada com o objetivo de que as simulações possam ser processadas em coordenadas cartesianas. O programa resolve situações tridimensionais, transientes e emprega o Método dos Volumes Finitos para integrar as equações diferenciais governantes. O modelo original é baseado no modelo de Gauthier et al., 1997, tendo sido os resultados do mesmo foram usados para a validação do presente estudo. Um circuito de água quente é também projetado e apresentado para o aquecimento das estufas. A água circula através de mangueiras sobre o solo e é aquecida por um sistema de queimadores a gás liqüefeito de petróleo ou óleo combustível, transferindo assim calor para o interior da mesma. O projeto de aquecimento foi realizado através de um programa de parceria entre a Ufrgs, Sebrae, Fapergs e a Agropecuária Clarice.

Interação de espirais em 2D : redução da dinâmica à interação de defeitos e exploração de novas estruturas espaço-temporais

O presente trabalho apresenta um anova proposta de tratamento de estruturas espirais em meios contínuos oscilatórios na vizinhança de bifurcações de Hopf supercríticas. Tais estruturas são normalmente descritas pela Equação de Cinzburg-Landau Complexa a qual usa um campo complexo associado a essas oscilações. A proposta apresentada reduz a dinâmica de espirais à interação entre os centros das mesmas. Inicialmente, comparamos numericamente as duas descrições e com os ganhos computacionais decorrentes da abordagem reduzida caracterizamos finamente as estruturas espaço-temporais formadas nesses sistemas: em vez dos estados congelados mencionados anteriormente na literatura encontrou-se uma dinâmica espaço-temporal intermitente. Esse regime ocorre em duas fases distintas: Líquido de Vórtices e Vidros de Vórtices. Esta última evolui em escalas de tempo ultralentas como fenômenos semelhantes encontrados na Mecânica Estatística, apesar de sua origem puramente determinista.

Comparação entre resultados experimentais e computacionais do comportamento térmico de um ambiente

O presente trabalho apresenta uma comparação entre resultados computacionais, obtidos através do "software" EnergyPlus e experimentais do comportamento térmico de um ambiente condicionado e não condicionado. Para tanto, monitoraram-se os dados climáticos de radiação, velocidade do vento e temperatura, no período de 11 a 20 de janeiro de 2002 e produziu-se um arquivo climático. Simultaneamente, fez-se a aquisição das temperaturas de uma sala-teste, localizada no terceiro pavimento de um prédio na cidade de Porto Alegre, bem como das salas adjacentes. As temperaturas do ar de insuflamento e de retorno dos condicionadores de ar, localizados na sala-teste, foram medidas durante o dia, em seis dias do período de monitoramento. Mediu-se também a velocidade do ar de retorno e determinou-se a potência sensível de refrigeração. Os ganhos de calor interno da sala também foram medidos e utilizaramse as formulações apresentadas pela ASHRAE, 2001, para determiná-los em relação às pessoas e à infiltração. Tais dados foram declarados ao programa como variáveis de entrada. As simulações do comportamento térmico da sala-teste foram implementadas informando-se ao programa a temperatura das salas ou os coeficientes de uma equação representativa das mesmas. Por considerar que a primeira representava melhor as condições térmicas das salas contíguas, utilizou-se esta modalidade para análise As simulações foram conduzidas, alterando-se opções fornecidas pelo programa: modelo de céu isotrópico e anisotrópico, coeficiente de convecção simples e detalhado. Os resultados da carga térmica e temperatura ambiente da sala-teste obtidos nas simulações foram comparados com os dados experimentais do período de monitoramento. A melhor concordância foi obtida para o modelo anisotrópico, coeficiente de convecção detalhado. Constatou-se uma grande discrepância entre as cargas térmicas, para o modelo de convecção simples. Assim, conclui-se que o "software" EnergyPlus representa bem o comportamento térmico de uma edificação "termicamente pesada" para coeficiente de convecção detalhado, necessitando pesquisa para as demais edificações.

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.

Atividade de predação de adultos de Cosmoclopius nigroannulatus Stal (Hemiptera: Reduviidae) sobre ninfas de Spartocera dentiventris (Berg) (Hemiptera: Coreidae)

O comportamento alimentar de insetos predadores determina o impacto que causam às populações com as quais interagem. Este trabalho teve como objetivo avaliar a resposta funcional e a extração de alimento por Cosmoclopius nigroannulatus Stal (Hemiptera: Reduviidae) em diferentes densidades de ninfas de 1° ínstar de Spartocera dentiventris (Berg) (Hemiptera: Coreidae). Os insetos foram obtidos de uma criação mantida em uma lavoura experimental de fumo em Porto Alegre, RS, Brasil. O experimento foi em laboratório (27 ± 1°C; 80 ± 5% UR; fotofase de 12h), sendo utilizados 10 adultos, recém-emergidos, de cada sexo de C. nigroannulatus em cada uma das cinco densidades (N) de ninfas (5, 15, 25, 35 e 45) de S. dentiventris. Os predadores foram observados individualmente por cinco dias, a cada 24h (T), registrando-se o número e o peso de ninfas mortas e/ou consumidas (Na), o peso do predador (Pb) e o das ninfas remanescentes e o tempo gasto para a ingestão do alimento (Ti). Foram ainda estimados o tempo de manuseio (Tm), tempo de busca (Tb), eficiência de busca (E) e taxa de ataque (a’). Em relação à extração de alimento foram estimadas a quantidade e o percentual de alimento extraído de cada ninfa (Pe), a quantidade extraída por minuto (Pe/min), as sobras (S) e a taxa de consumo relativa (RCR). Para a avaliação da extração de alimento cada tratamento foi pesado, com os seguintes resultados: 5,1 14,7, 29,8, 30,6 e 44,8mg de ninfas (oferecidas aos machos) e 5,5, 14,7, 31,1, 37,4 e 48,5mg de ninfas (oferecidas às fêmeas). fêmeas). Foi encontrada correlação positiva entre Na e a densidade (N), sendo que as fêmeas ingeriram mais ninfas que os machos. Tm foi maior nos machos (3,0h) que nas fêmeas (1,9h), tendo o tempo de manuseio total (Tm x Na) aumentado com a densidade. Tb, E e a' apresentaram correlação negativa em relação à densidade. Os componentes da resposta funcional foram ajustados ao modelo randômico da equação dos discos de Holling (Na= N {1-exp[-a’(T-TmNa)]}), evidenciando o tipo II (curvilinear) de resposta funcional como característico da espécie. Com o aumento da quantidade de ninfas oferecidas evidenciou-se aumento no consumo total de ninfas ingeridas em mg (Na), bem como em Pe, Pb, S e RCR. Já em relação ao Na e Pe, em termos percentuais, e Pe/min verificou-se uma diminuição, sendo que as fêmeas consumiram e ganharam mais peso que os machos. Com os machos, a estabilização de Na e Pe ocorreu em torno da quantidade 29,8 mg de ninfas oferecidas. Quanto às fêmeas, a estabilização do consumo foi por volta de 31,1 mg de ninfas oferecidas. Não se encontrou correlação entre o tempo de ingestão (Ti) e a quantidade de ninfas, sendo que as fêmeas ingerem mais rapidamente (19,9 ± 2,26min) que os machos (23,9 ± 1,29min). Os resultados sugerem que o comportamento alimentar de C. nigroannulatus sobre ninfas de 1° ínstar de S. dentiventris está diretamente relacionado à quantidade de presas disponíveis, podendo ser considerado um agente potencial para o controle biológico.

Modelo numérico para análise à flexão de elementos estruturais com protensão aderente e não aderente

Na protensão não aderente, a armadura permanece livre para se movimentar ao longo de seu perfil em todas as seções com exceção das de ancoragem. Não há aderência entre concreto e armadura, e a hipótese da compatibilidade de deformações entre o aço e concreto não é aplicável, tornando inviável o desenvolvimento de uma solução analítica. Visando colaborar para a maior compreensão do comportamento à flexão das estruturas com protensão não aderente e para o desenvolvimento de critérios nacionais de projeto, um modelo numérico foi implementado utilizando o elemento finito do tipo híbrido para pórticos planos. Nesta formulação, a equação para as solicitações ao longo do elemento é a função de interpolação, e as forças as variáveis interpoladas. Como esta função resulta das condições de equilíbrio, sem hipóteses arbitrárias, o método é considerado exato para forças e curvaturas. Elementos longos são possíveis, de maneira que um único elemento finito pode ser utilizado para um vão de viga ou pilar, reduzindo o esforço computacional. O caráter exato da formulação contribui para a boa modelagem dos cabos não aderentes, já que a tensão nestas armaduras depende das curvaturas de todas as seções do elemento. O modelo numérico proposto prevê a não linearidade geométrica, carregamentos cíclicos e a construção composta. Relações constitutivas já consolidadas na literatura são empregadas para os materiais. A cadeia de Maxwell é utilizada para representar o comportamento reológico do concreto e do aço de protensão, respeitando as características de envelhecimento de cada material. Inúmeros exemplos são apresentados, permitindo a comparação entre resultados numéricos e experimentais Uma análise paramétrica foi realizada, caracterizando o desempenho do modelo numérico frente a variações nos parâmetros de entrada. Discutem-se, ainda, dois critérios de ruptura para a utilização do modelo numérico no estudo do comportamento à flexão das estruturas com protensão não aderente. Os resultados indicam o bom desempenho do modelo numérico e a sua adequação para a realização de pesquisas sobre o assunto.

Derivação e solução de equações modelo da dinâmica de gases rarefeitos

Neste trabalho, duas equações modedlo na área da dinâmica de gases rarefeitos, são derivadas a partir de algumas soluções exatas da equação linearizada de Boltzmann homogênea e não homogênea. Em adição, uma versão analítca do método de ordenadas discretas é usado para resolver problemas clássicos nesta área, descritos pelo "Modelo S". Resultados numéricos são apresentados para os problemas de fluxo de Couette, fluxo de Poiseuille, "Creep" Térmico, Deslizamento Térmico e problema de Kramers.

Desenvolvimento de um método de volumes finitos de alta ordem para a simulação de escoamentos de fluidos viscoelásticos

O uso da mecânica de fluidos computacional no estudo de processos envolvendo o escoamento de fluidos poliméricos está cada vez mais presente nas indústrias de transformação de polímeros. Um código computacional voltado a esta função, para que possa ser aplicado com sucesso, deve levar a predições mais próximas possível da realidade (modelagem), de uma forma relativamente rápida e eficiente (simulação). Em relação à etapa de modelagem, o ponto chave é a seleção de uma equação constitutiva que represente bem as características reológicas do fluido, dentre as diversas opções existentes. Para a etapa de simulação, ou seja, a resolução numérica das equações do modelo, existem diversas metodologias encontradas na literatura, cada qual com suas vantagens e desvantagens. Neste tópico se enquadra o trabalho em questão, que propõe uma nova metodologia para a resolução das equações governantes do escoamento de fluidos viscoelásticos. Esta se baseia no método dos volumes finitos, usando o arranjo co-localizado para as variáveis do problema, e na utilização de aproximações de alta ordem para os fluxos médios lineares e não-lineares e para outros termos não lineares que surgem da discretização das equações constitutivas. Nesta metodologia, trabalha-se com os valores médios das variáveis nos volumes durante todo o processo de resolução, sendo que os valores pontuais são obtidos ao final do procedimento via deconvolução. A solução do sistema de equações não lineares, resultante da discretização das equações, é feita de forma simultânea, usando o método de Newton São mostrados então, resultados da aplicação da metodologia proposta em problemas envolvendo escoamentos de fluidos newtonianos e fluidos viscoelásticos. Para descrever o comportamento reológico destes últimos, são usadas duas equações constitutivas, que são o modelo de Oldroyd-B e o modelo de Phan-Thien-Tanner Simplificado. Por estes resultados pode-se ver que a metodologia é muito promissora, apresentando algumas vantagens frente às metodologias convencionais em volumes finitos. A implementação atual da metodologia desenvolvida está restrita a malhas uniformes e, consequentemente, soluções para problemas com geometrias complexas, que necessitam de refinamento localizado da malha, foram obtidas somente para baixos números de Weissenberg, devido a limitação do custo computacional. Esta restrição pode ser contornada, tornando o seu uso competitivo.

Vibrações livres e forçadas no modelo de Timoshenko

O objetivo deste trabalho é o cálculo modal da resposta impulso distribu ída de uma viga descrita pela equação de Timoshenko e das vibrações forçadas, devidas a influência de cargas externas. Os modos vibratórios foram obtidos com o uso da base dinâmica, gerada por uma resposta livre e suas derivadas. Esta resposta é caracterizada por condições iniciais impulsivas. Simulações foram realizadas para os modos, a resposta impulso distribuída e vibrações forçadas em vigas apoiadas em uma extremidade e na outra livre, fixa, deslizante ou apoiada, sujeitas a cargas oscilatórias espacialmente concentradas ou distribuídas através de pulsos.

Sistemas concentrados e distribuídos através da análise modal adjunta

Este estudo visou a utilização direta da análise modal adjunta em sistemas não-clásicos de natureza concentrada e distribuída, explorando o conceitos de biortogonalidade modal e da resposta impulso evolutiva e estacionária. Discutiu-se a natureza do espectro para diversas classes de sistemas, em particular para sistemas com efeitos giroscópicos e com atrito interno. A teoria foi desenvolvida de maneira paralela para sistemas de primeira ordem e de segunda ordem, porém é apresentada de maneira independente, ou seja, para os sistemas de segunda ordem não é utilizada a formulação de espaço de estado, a qual reduz o sistema para primeira ordem. Assim, as relações de biortogonalidade para sistemas de segunda ordem são apresentadas e utilizadas de maneira direta na expansão modal da resposta dos sistemas. A forma dos modos de vibração em variadas aplicações é determinada de maneira exata com o uso da base dinâmica gerada pela resposta impulso espacial. No cálculo dos autovalores foi introduzida uma aproximação polinomial para a resposta impulso espacial Os coeficientes dessa aproximação foram obtidos por recursão, a partir de uma equação em diferenças associada à equação característica do problema modal. Simulações numéricas foram realizadas para obter a resposta impulso evolutiva, respostas forçadas e modos de vibração de sistemas não-clássicos concentrados, formulados através de modelos ou aproximações, e sistemas distribuídos, formulados através de modelos ou incluindo distúrbios e acoplamento através das condições decontorno. Os resultados deste estudo permitiram concluir, através das simulações numéricas realizadas, a importância da base dinâmica no sentido de simplificar os cálculos para obtenção dos autovalores, dos modos de vibração e, consequentemente, da resposta do sistema, seja concentrado, distribuído com ou sem acoplamento das condições de contorno. A utilização da análise modal adjunta, desde que ocorra uma ordenação adequada dos autovalores e modos, mostrou-se um método eficiente na obtenção direta da resposta de sistemas não-clássicos de segunda ordem, ou seja, sem redução ao espaço de estado.