Modelação da circulação oceânica no Arquipélago de Cabo Verde

O conhecimento dos oceanos é limitado pela disponibilidade de dados disponíveis. Este estudo mostra como modelos numéricos podem aumentar o conhecimento com pouca informação local, e de forma económica. Neste trabalho é efetuado a implementação do modelo numérico MOHID, bidimensional (2D) para simular a propagação da maré e tridimensional (3D) para estudar o escoamento baroclíncio no Arquipélago de Cabo Verde. O modelo 2D é baseado num modelo de escala regional com passo espacial de 6 km (nível 1), com dois modelos encaixados de passo espacial de 3 km (nível 2). Também foi utilizado para impor a maré (nível 0) no modelo 3D (nível 1, com passo de 6 km). O modelo tridimensional tem 50 camadas na vertical e fornece informação sobre a circulação geral. A fronteira lateral aberta foi forçada pelo modelo global de maré FES 2004 e pelos resultados do projecto MyOcean fornecendo condições de temperatura, salinidade e níveis médios diários para o modelo 3D. O forçamento meteorológico foi fornecido pelos resultados do modelo meteorológico GFS (Global Forecasting System). Os resultados dos modelos foram analisados com as informações conhecidas como as principais correntes da região, as imagens de satélite de temperatura da superfície do mar, perfis verticais obtidos pelas bóias argo e os resultados do projeto MyOcean. O modelo reproduz o padrão de circulação conhecido para a região e os resultados apresentam boa concordância com a informação disponível. O modelo mostra que o vento tem influência sobre o escoamento à superfície, mas tem pouca influência sobre os níveis.

Caracterização hidrodinâmica e hidrodispersiva de um latossolo amarelo na microrregião do brejo paraibano

Para modelar o fluxo de água e o transporte de soluto em solos são necessários parâmetros hidrodinâmicos e hidrodispersivos como dados de entrada nos modelos numéricos. O experimento de infiltração axissimétrica utilizando um infiltrômetro de anel simples com um traçador conservativo (KCl) é um método simples e eficiente para determinar esses parâmetros no campo. O objetivo deste trabalho consistiu na caracterização hidrodinâmica e hidrodispersiva de um Latossolo Amarelo situado na microrregião do brejo paraibano, no município de Areia, PB, em nível de campo, por meio da infiltrometria de anel simples. Os experimentos foram realizados em uma área cultivada com feijão em uma malha de 50 x 50 m. A condutividade hidráulica saturada, Ks, e a sorvidade, S, foram estimadas pelas análises de infiltração tridimensional axissimétrica para tempos curto e longo. O soluto inerte (Cl-) foi utilizado para o cálculo da fração de água móvel, Φ, por meio da medição da concentração de soluto na camada do solo (0-15 cm) ao final do experimento de infiltração. Um modelo numérico de transferência de soluto com base no conceito de convecção-dispersão (CD) foi utilizado para definir os parâmetros (D e R) de transporte. A infiltrometria de anel simples demonstrou ser um método simples e eficiente para obtenção dos parâmetros hidrodinâmicos e hidrodispersivos. Observou-se boa concordância entre os valores medidos e calculados (v, D e R) com o modelo CD observado pelos coeficientes de determinação (r²).

Sistema de visualización SIG del Servei Meteorològic de Catalunya

El sistema de visualización SIG del Servei Meteorològic de Catalunya (SMC) consiste en un sistema que visualiza la información meteorológica del SMC mediante un visor SIG. La herramienta permitirá a los técnicos del SMC poder consultar productos meteorológicos, como modelos numéricos, productos de Meteosat, radar, descargas eléctricas, etc. El sistema se compone de un servidor MapServer, utilizando MapCache como sistema de caché, y de un visor SIG web implementado en JavaScript, utilizando las librerías Dojo Toolkit y OpenLayers.

Emprego de ferramentas numéricas na avaliação do módulo de elasticidade em vigas roliças de madeira

Vigas são elementos estruturais encontrados na maioria das construções civis. Dentre os materiais de engenharia, destaca-se a madeira, por ter resistência mecânica satisfatória aliada a baixa densidade. A madeira roliça apresenta-se como boa solução na confecção de vigas, uma vez que não precisa ser processada, como é o caso da madeira serrada. O projeto de elementos estruturais de madeira requer o conhecimento de suas propriedades físicas e mecânicas, obtidas segundo as premissas de documentos normativos. Em se tratando da madeira roliça, os documentos normativos nacionais que tratam da determinação das propriedades de resistência e rigidez estão vigentes há mais de vinte anos sem revisão técnica. De forma geral, tanto as normas nacionais como as internacionais idealizam geometria troncocônica para as peças roliças de madeira, implicando equações simplificadas incapazes de prever a influência das irregularidades da forma na determinação do módulo de elasticidade longitudinal. Este trabalho objetiva avaliar a influência das irregularidades da geometria em peças roliças de madeira Corymbia citriodora e Pinus caribaea no cálculo do módulo de elasticidade longitudinal. Para tanto, utilizou-se do ensaio de flexão estática a três pontos, considerando também um modelo matemático simplificado, assumindo seção circular constante para a forma do elemento. As irregularidades das peças são consideradas nos modelos numéricos, constituídos de elementos finitos de barra e tridimensionais. Os resultados encontrados revelam equivalência estatística entre os módulos de elasticidade para ambas as formas de cálculo, indicando ser plausível a consideração de seção circular constante para as peças de madeira aqui avaliadas.

Modelamento geomecânico para o dimensionamento de pilares de carvão

O dimensionamento de pilares em mineração de carvão, no Brasil, sempre foi considerado conservativo em relação às metodologias e características das camadas que estava sendo lavrada. Torna-se necessário reverter o status-quo existente nesse tema para que a atividade de mineração em carvão se torne mais lucrativa e segura ao mesmo tempo. Problemas geomecânicos, que causaram acidentes fatais, diminuição das reservas de carvão, alto custo de extração e baixa recuperação na lavra são decorrentes ao baixo conhecimento do comportamento geomecânico do maciço rochoso e das metodologias de modelamento geomecânicos. Esse trabalho propõe uma nova técnologia de mapeamento geomecânico e caracterização geomecânica para determinar parâmetros geomecânicos das rochas que formam o sistema piso-pilarteto de uma camada de carvão. O dimensionamento de pilares, então, deve ser baseado nessas informações que alimentam um modelo geomecânico que utiliza métodos numéricos para solucionálos. Dentro do processo de aquisição de informações geomecânicas são propostos novos métodos de mapeamento geomecânico por janela amostral e o cálculo dos índices espaçamento médio e freqüência de descontinuidades. Baseado nas características quantitativas das descontinuidades é proposto o cálculo do RQD* (Rock Quality Designation teórico) para a camada de carvão com base na função de distribuição que se ajusta ao histograma de distribuição do espaçamento de descontinuidades. Esse procedimento refletiu melhor a qualidade do maciço rochoso, baseado no índice RMR (Rock Mass Rating). Considerando-se a qualidade do maciço rochoso e utilizando critérios de rupturas empíricos, foram estimados os parâmetros geomecânicos das rochas que alimentaram os modelos geomecânicos.(Continua)_ Foram gerados modelos numéricos para simular situações reais encontradas na camada Bonito em Santa Catarina e propor um novo paradigma para dimensionamento de pilares, baseado em simulação numérica, considerando todo o entorno do pilar, compreendendo o sistema piso-pilarteto. Os resultados desses modelos são comparados com o monitoramento das deformações dos pilares para confrontar as simulações numéricas com o comportamento geomecânico do maciço insitu. Os resultados apresentaram grande qualidade e concordância com o comportamento dos pilares apresentados no monitoramento das deformações, provando que a metodologia proposta aqui pode ser utilizada em futuros projetos de dimensionamento de pilares.

Simulação do movimento da água no solo utilizando modelo numérico

Um dos principais desafios para o século XXI é o uso sustentável dos recursos hídricos. A agricultura é a maior usuária desse recurso, sendo responsável por aproximadamente 70% dos gastos globais de água. Para o uso racional da água é necessário o conhecimento das propriedades hidráulicas do solo na condição de não-saturação. A simulação do movimento da água no solo através da modelagem é uma ferramenta importante no entendimento do sistema solo-água-planta, pois permite a previsão e estudo do comportamento da redistribuição da água no solo de maneira eficiente. Porém a forte dependência entre os parâmetros hidráulicos do solo (teor de água, potencial matricial e condutividade hidráulica) fazem da modelagem um assunto complicado, necessitando de modelos numéricos que utilizem intervalos discretos de tempo e espaço. Esses intervalos devem ser suficientemente pequenos para que dentro de um intervalo as variações dos parâmetros hidráulicos sejam insignificantes. Nesta tese é proposto um algoritmo para a descrição do movimento vertical da água no solo, definindo o intervalo de tempo como função de uma variação máxima admissível do teor de água. O algoritmo foi testado para alguns conjuntos de dados hidráulicos de solos, utilizando as condições de contorno de um experimento de perfil instantâneo. A eficiência do algoritmo foi verificada em situações práticas, como na previsão do teor de água na Capacidade de Campo e na predição do tempo para ser atingida determinada condição de umidade no solo. Utilizando o algoritmo, também foram estudados o comportamento da retenção e dos parâmetros de condutividade em relação aos parâmetros de drenagem. Concluiu-se que o algoritmo descreve adequadamente o processo de redistribuição de água no solo utilizando intervalos suficientemente pequenos de tempo e espaço; a forma de discretização do tempo não comprometeu a acurácia das medidas; a previsão do teor de água na Capacidade de Campo e do tempo para ser atingida determinada condição de umidade no solo foram satisfatórias; e o algoritmo permite o estudo das relações entre os parâmetros hidráulicos do solo.

Influência de parâmetros físicos no dimensionamento de chaminés de equilíbrio simples de usinas hidrelétricas

Tendo em vista que as chaminés de equilíbrio são estruturas importantes na proteção dos circuitos hidráulicos das usinas hidrelétricas contra os fenômenos transientes que ocorrem rotineiramente durante sua operação, realizou-se uma investigação sobre a influência dos parâmetros físicos no desempenho e conseqüentemente no dimensionamento de chaminés de equilíbrio simples localizadas a jusante das turbinas, caso aplicado a Usina Hidrelétrica de Serra da Mesa. Os parâmetros estudados foram: área da chaminé, comprimento e área do túnel de fuga, perda de carga, vazão, nível no canal de fuga e tempo de manobra. Foram coletados dados na Usina Hidrelétrica de Serra da Mesa e em uma estrutura laboratorial existente no LAHE, ambas pertencentes a FURNAS Centrais Elétricas S.A, referentes a situações de regimes permanentes (operação normal) e de regimes transientes. Estes dados foram utilizados na determinação da perda de carga linear nos túneis e das características das oscilações em suas chaminés de equilíbrio. Para o estudo da perda de carga linear em túneis, também foram feitas medições, em regime permanente, na Usina Hidrelétrica de Salto Forqueta – CERTEL - Cooperativa de Eletrificação Teutônia, e utilizados os dados da Usina Hidrelétrica Bugres fornecidos pela CEEE – Companhia Estadual de Energia Elétrica –RS. Com estes dados foram estabelecidas relações entre o coeficiente de atrito e o número de Reynolds para cada estrutura e realizada a sua comparação com valores teóricos Verificou-se, através da comparação com dados experimentais que os quatro modelos numéricos disponíveis estavam aptos à aplicação na seqüência da pesquisa, porém, elegeu-se o modelo CHAMINE-IPH devido a sua maior versatilidade e simplicidade de uso. Com a análise dos resultados das simulações numéricas realizadas foi possível a identificação da influência dos parâmetros físicos no comportamento do nível d’água em chaminés de equilíbrio simples durante regimes transientes. Constatou-se que para a reprodução numérica de fenômenos transitórios em chaminés de equilíbrio a jusante das turbinas, é fundamental a consideração da oscilação do nível no canal de fuga e dos coeficientes de perda de carga adequados.

Análise numérica e experimental da solidificação das ligas Al 3,5% Cu e AA36

A formação interna da estrutura de uma peça fundida deve determinar as características mecânicas dessa peça, quando em serviço. Dessa forma, torna-se fundamental conhecer e parametrizar os elementos que influem na formação da estrutura do componente sólido, para o desenvolvimento adequado de sistemas de solidificação que permitam conduzir a melhoria da qualidade de uma peça no processo de solidificação. Modelos numéricos podem ser empregados para a determinação da formação estrutural, entretanto estes devem ser aferidos adequadamente para garantir sua aplicabilidade. No trabalho em questão, desenvolveu-se uma metodologia que permitiu a aferição de um modelo numérico para a previsão das condições de resfriamento de ligas de alumínio. Elaborou-se um conjunto de coquilhas e moldes que permitiram o acoplamento e verificação da estrutura formada na solidificação em diferentes situações de resfriamento, observando principalmente a direção de extração de calor e a taxa de resfriamento associada. O modelo numérico aplicado respondeu coerentemente aos resultados obtidos permitindo a previsão da estrutura solidificada pela observação dos resultados de evolução das isotermas solidus e liquidus no interior dos sistemas.

Simulação computacional do comportamento de estruturas de aço sob incêndio

A análise do comportamento estrutural sob incêndio constitui uma parte importante da engenharia de proteção contra incêndio, especialmente no caso de estruturas de aço, por sua alta condutividade térmica e relativa esbeltez das seções. As dificuldades econômicas e práticas, associadas à avaliação do comportamento estrutural, por meio de ensaios em escala real, têm estimulado o desenvolvimento e uso de métodos de simulação numérica. Esta tese trata da simulação numérica do comportamento de estruturas de aço sob condições de incêndio e se divide em três partes. As duas primeiras partes foram desenvolvidas na Universidade de Liége, na Bélgica, usando-se o programa SAFIR como ferramenta numérica. A terceira parte foi desenvolvida de forma independente, na Universidade Federal do Rio Grande do Sul. Na primeira parte, é feito um estudo comparativo entre o uso de elementos finitos de viga e de casca, na modelagem de vigas simplesmente apoiadas, sujeitas a flambagem lateral por torção. Os esforços de torção, presentes no caso de flambagem lateral, podem levar a uma plastificação da seção transversal e à conseqüente redução da rigidez à torção da seção. Da mesma forma, a degradação das propriedades do material, com o aumento da temperatura, leva à redução da rigidez à torção Havia dúvidas se o modelo com elementos de viga, com uma rigidez à torção constante, poderia fornecer uma resposta aceitável. O estudo mostrou que uma resposta com boa precisão pode ser conseguida, usando-se elementos de viga, desde que o módulo de elasticidade transversal seja ajustado para refletir seu valor correspondente à temperatura de falha. Isso implica um processo iterativo, uma vez que a temperatura de falha não é previamente conhecida. Por outro lado, a degradação da rigidez à torção, por efeitos de plastificação, pode ser ignorada. Na segunda parte, é feita a comparação entre as modelagens bidimensional e tridimensional, de galpões industriais de um andar, sob incêndio. Comumente, a estrutura de galpões industriais é composta por pórticos-tipo, dispostos em paralelo. A análise desses galpões é comumente feita pela simulação no plano do pórtico de aço ou, simplesmente, da treliça da cobertura Na análise bidimensional, importantes efeitos fora do plano são ignorados, como a redistribuição de esforços devido à degradação do material ou à expansão térmica, ou instabilidade lateral dos elementos. A importância desses efeitos e a adequabilidade do modelo 2D para representar o comportamento real são discutidas. Na terceira parte, um modelo numérico para a simulação tridimensional do comportamento de estruturas de aço sob incêndio é apresentado. O modelo é baseado no conceito de rótulas plásticas generalizadas, com modificações para melhor representar a formação e expansão da plastificação no elemento. A descrição cinemática adotada permite obter bons resultados, mesmo com o uso de poucos elementos. A determinação do vetor de esforços internos no elemento, incluindo os efeitos da temperatura, é detalhada. O procedimento foi validado por comparação com ensaios e com modelos numéricos mais sofisticados, como o programa SAFIR. Os resultados demonstram que o procedimento proposto pode ser usado como uma forma alternativa de análise 3D de estruturas sob incêndio, com precisão razoável e baixo esforço computacional.

Estudo teórico e experimental das ligações diretas contra-fiadas entre paredes de blocos de concreto em escala real e reduzida 1:4

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

Análise teórica e experimental vibro-acústica utilizando a técnica de matrizes compactas

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

Análise dinâmica do amortecimento em estruturas compostas com material viscoelástico

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

Controle térmico de componentes eletrônicos com tubos de calor

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

Avaliação do impacto do reservatório de Três Irmãos sobre a superfície potenciométrica do aquífero livre na cidade de Pereira Barreto (SP): uma abordagem numérica e geoestatística

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

Programa computacional com interface gráfica para identificação estocástica de parâmetros modais de estruturas civis: aplicação em pontes e torres de linha de transmissão

O interesse no comportamento dinâmico de estruturas metálicas vem crescendo nas últimas décadas no Brasil, em decorrência de acidentes com colapso total de algumas estruturas devido às vibrações ambientes em diversas regiões do país. Na região amazônica, por exemplo, onde esse tipo de estrutura deve vencer obstáculos como florestas e rios de grande largura, casos de colapso total de estruturas metálicas também são relatados. O foco principal dessa dissertação é o estudo do comportamento modal de estruturas metálicas submetidas às vibrações ambientes cuja magnitude das forças de excitação é desconhecida. Dois estudos de caso são apresentados: no primeiro deles, o comportamento modal de uma torre de linha de transmissão de energia elétrica é investigado; e no segundo caso, tanto o comportamento modal como os níveis de desconforto de uma ponte são estudados. Os estudos realizados neste último caso visam avaliar os níveis de desconforto da ponte quando submetida às excitações ambientes como rajadas de vento e o tráfego de veículo de acordo a norma brasileira NBR 8800 (1986). Em ambos os estudos de caso foram realizadas análises experimentais e computacionais. Na etapa experimental, ambas as estruturas foram monitoradas com emprego de um conjunto de acelerômetros de baixa freqüência e também de um sistema de aquisição apropriados para ensaios de vibração de estruturas civis. Como é muito difícil medir a magnitude das forças de excitação ambientes, foram utilizados os métodos de identificação estocásticos SSI-DATA e SSI-COV para extração de parâmetros modais de estruturas civis a partir somente dos dados de resposta coletados nos ensaios de vibração. Entre as atividades desenvolvidas nessa etapa, destaca-se a criação de um programa computacional com recursos do Graphical User Interface (GUI) da plataforma Matlab®, destinado à identificação modal de estruturas civis com o emprego dos referidos métodos estocásticos. Esse programa é constituído de três módulos: o primeiro é destinado ao processamento e tratamento dos sinais coletados nos ensaios de vibração; o segundo é utilizado para adicionar as informações do posicionamento dos acelerômetros utilizados nos arquivos dos sinais de resposta; e o terceiro e último módulo é destinado à identificação a partir dos arquivos de dados de resposta processados nos dois primeiros módulos. Na etapa das análises teóricas, foram criados modelos numéricos utilizando o método dos elementos finitos para simular o comportamento dinâmico das estruturas analisadas. Comparando os resultados obtidos em ambas as etapas de análise, verifica-se que resultados experimentais e teóricos apresentaram parâmetros bastante próximos entre si nos primeiros modos de vibração. Os resultados experimentais mostraram que ambos os métodos estocásticos foram muito eficientes na identificação das estruturas ensaiadas.