Análise numérica de uma torre de resfriamento de grande porte

Este trabalho tem como objetivo principal determinar a troca de calor em uma torre de resfriamento de grande porte, do tipo seca. A força motriz deste sistema é o empuxo resultante da diferença de temperatura entre o ar interior e exterior da torre. A partir de uma certa velocidade o vento externo também desempenha um papel importante. A convecção natural ocorre em regime turbulento, devido às grandes dimensões. A modelagem numérica é feita pelo Método dos Volumes Finitos com malhas criadas sobre Elementos Finitos, através dos softwares comerciais Fluent e CFX. Os resultados são obtidos para diferentes potências dissipadas no trocador de calor. A partir dos resultados são realizados estudos sobre a viabilidade de mudanças na geometria da torre, a fim de se obter melhorias no escoamento, bem como uma troca de calor mais eficiente. Faz-se ainda uma regressão a fim de obter correlações entre a potência dissipada no trocador de calor e a vazão mássica de ar necessária para a troca. Estas correlaçõoes são de utilidade, para aplicações posteriores, na implantação de um simulador para o processo de troca de calor nesta torre. O trocador de calor também é enfocado a fim de se determinar o coeficiente global de troca térmica, a fim de disponibilizá-lo para aplicação em outros trabalhos

Uso do método das interfaces coesivas na análise da simulação do comportamento de fratura em materiais frágeis

Este trabalho apresenta uma abordagem numérica da fratura baseada no Método dos Elementos Finitos na qual são considerados materiais frágeis como objeto de estudo. As simulações realizadas são possíveis através de um algoritmo computacional que consiste na implementação de molas entre as faces dos elementos finitos. Esta composição constitui o método denominado Interfaces Coesivas. As molas fornecem o efeito de coesão e é representado por leis constitutivas do tipo bi-linear. Sua descrição é apresentada, bem como sua implementação na análise. Tanto o modo I como o modo II de ruptura são considerados. Com o objetivo de verificar o presente modelo com outros resultados e soluções, são analisados exemplos numéricos e experimentais extraídos da literatura. Tais exemplos tratam da nucleação, propagação e eventuais paradas de trincas em peças estruturais submetidas aos mais diversos carregamentos e condições de contorno. É ainda discutido como os parâmetros de fratura do material interferem no processo de ruptura. Uma tentativa é feita de determinar estes parâmetros, em especial no caso de concreto.

Análise dos fatores influentes na tensão última de protensão em cabos não aderentes

A protensão não aderente é um sistema de pós-tensão caracterizado pela liberdade de deslizamento da armadura em relação ao concreto, ao longo de todo o perfil do cabo, com exceção das ancoragens. Devido à falta de aderência entre aço e concreto, a compatibilidade de deformações dos materiais na seção transversal não existe. O que se verifica é a compatibilidade de deslocamentos, com a equivalência entre os alongamentos do cabo e das fibras de concreto adjacentes ao mesmo. Isto acarreta complexidade no projeto de tais estruturas. No que se refere ao dimensionamento no Estado Limite Último, a tensão atuante nas armaduras não aderentes, também denominada tensão última de protensão, não é, a princípio, conhecida. Para sua obtenção, é fundamental a precisa determinação das curvaturas ao longo do elemento, tornando o desenvolvimento analítico de uma solução, praticamente inviável. Geralmente, são utilizados critérios empíricos para a previsão da tensão última de protensão em armaduras não aderentes. Estes levam em conta na sua formulação, parâmetros considerados como de grande influência no valor da tensão última de protensão. A fim de avaliar a significância de alguns destes fatores no valor da tensão última de protensão, um estudo paramétrico foi realizado. Parâmetros como a taxa de armadura, a tensão efetiva inicial de protensão, a relação entre a altura do elemento e o seu vão e, também, o tipo de carregamento, foram investigados. Para tal, um protótipo foi idealizado por meio de um modelo numérico Este modelo utiliza a formulação de elemento finito do tipo híbrido para pórticos planos, que é caracterizado pela precisa obtenção das curvaturas, e mostrou-se adequado para utilização na presente pesquisa. A análise dos resultados obtidos permitiu a identificação e quantificação da influência dos parâmetros estudados, no valor da tensão última de protensão. A relevância dos valores de tensão obtidos, em função dos parâmetros adotados, na capacidade portante dos elementos, foi também avaliada.

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.

Ferramentas para análise dinâmica e estrutural de um reboque de linha leve

Neste trabalho é desenvolvida uma metodologia de projeto para identificar as regiões críticas da estrutura de um reboque de linha leve sendo tracionado em pavimentos do tipo rodovia de baixa qualidade e estrada secundária de muito baixa qualidade. Para tanto, são levantados alguns dados experimentais da estrutura, necessários para a aproximação e simulação dinâmica de um modelo simplificado. A excitação da base é realizada por atuadores que simulam as oscilações verticais de um perfil de estrada, a qual é definida de acordo com os estudos realizados por Dodds e Robson (1973). Isto permite a determinação de um histórico de carregamentos das regiões da estrutura do chassi sob a ação das molas da suspensão. Em seguida, é gerado um modelo estrutural simplificado do reboque em elementos finitos, chamado de global, no qual são determinadas as regiões sob ação das maiores tensões. Tendo identificada a região mais crítica da estrutura, é criado um modelo local desta parte, onde se pode observar a distribuição de tensões com mais detalhe, permitindo a identificação dos pontos de concentração de tensões. Desta forma, com a aplicação do método de análise global-local é possível a obtenção de resultados detalhados quanto aos esforços da estrutura com um menor custo computacional.

O modelo de Gurson para dano dúctil : estratégia computacional e aplicações

Este trabalho trata do emprego do modelo de dano dúctil de Gurson e de alguns aspectos relativos a sua implementação computacional. Emprega-se como ponto de partida o programa Metafor, código de elementos finitos voltado para a simulação de problemas de conformação mecânica, contemplando portanto grandes deformações plásticas. Estudam-se alguns casos simples procurando identificar a influência de cada parâmetro do modelo de Gurson na resposta. É discutida a aplicação dos modelos de nucleação de vazios usualmente empregados em conjunto com a superfície de escoamento de Gurson em situações onde há reversão de solicitação. São propostas alterações nos modelos usuais de nucleação de forma a minorar resultados incoerentes verificados. Apresenta-se um algoritmo alternativo ao método de Newton-Raphson para a solução do sistema de equações associado à forma fraca das equações de equilíbrio, em situações onde a rigidez fica muito pequena. Tal algoritmo é implementado e testado. Um algoritmo de integração de tensões é apresentado e implementado, mostrando-se sua vantagem do ponto de vista de robustez, ou seja, obtenção da resposta correta a despeito do tamanho de incremento empregado. Também é discutido um procedimento para desconsiderar elementos que apresentem ruptura local Estudam-se problemas envolvendo geração de calor por deformação plástica em conjunto com a formulação do modelo de Gurson, mostrando que a consideração de acoplamento termo-mecânico com dano traz vantagens em algumas simulações. É discutida a implementação do modelo de Gurson sobre a versão tridimensional do Metafor, que contempla a formulação lagrangiana-euleriana arbitrária (LEA). Ao final, é feita uma simulação inicial do comportamento de espumas metálicas.

Estudo e desenvolvimento de um transdutor de posição linear magnética diferencial

Este trabalho tem como objetivo o estudo e o desenvolvimento de um transdutor magnético de posição, cujo princípio de funcionamento é baseado na relação entre a distribuição do fluxo magnético em sua estrutura, provido por uma fonte de fluxo magnético constante, c.c. ou c.a., e a posição linear da parte móvel do transdutor ao longo de sua direção de deslocamento ou excursão. O transdutor pode ter entreferros de medição para medir a indução magnética nas culatras laterais, com fluxo magnético de excitação constante, c.c. ou c.a., ou não ter entreferros de medição, caso em que apenas fluxo magnético de excitação c.a. é empregado. Em ambos os casos, com fluxo magnético de excitação c.a., bobinas exploradoras foram usadas para medir o fluxo magnético estabelecido nas culatras de acordo com a posição da parte móvel do transdutor. Para o primeiro caso foi desenvolvido um protótipo no qual o material das culatras é aço maciço. Para o segundo foram desenvolvidos dois protótipos, um com barra de aço maciça, e o outro com núcleo de lâminas de Fe-Si. Além de se desenvolver uma modelagem analítica para estudar o transdutor, empregou-se, como ferramenta auxiliar de simulação, uma modelagem eletromagnética por elementos finitos. Os resultados teóricos, simulados e experimentais são esclarecedores e demonstram que o transdutor confeccionado com lâminas de Fe-Si e excitação através de fluxo magnético c.a. apresenta uma saída linear. É demonstrado também que a metodologia proposta e utilizada na modelagem e na concepção do transdutor são adequadas e indicam o domínio da teoria e da tecnologia propostas. As configurações sugeridas são extensivamente analisadas por este trabalho No caso do transdutor com fluxo magnético c.a. proposto e resultante desse estudo, sua topologia pode ser considerada inovadora, não sofrendo um efeito de histerese no seu sinal de saída devido à histerese magnética de suas culatras ferromagnéticas como acontece no transdutor com fluxo magnético c.c.. Isso, entre outras razões, ensejou o depósito de sua patente. Adicionalmente, é robusto, de fácil confecção, com enrolamentos que não requerem uma montagem de precisão como em outros transdutores eletromagnéticos de posição, e de baixo custo.

Aproximações de Galerkin clássico e Galerkin/mínimos quadrados de escoamentos viscosos incompressíveis : aplicação à medição de fluxos

Desde a antigüidade a medição do escoamento dos fluidos tem sido uma marca de nossa civilização, ajudando a predizer a fertilidade das terras e o consumo d’água em fontes e aquedutos. Nos nossos dias, a área de medição de fluxo está bem estabelecida e ainda desperta grande interesse nas linhas de pesquisa da mecânica dos fluidos experimental e computacional. Em particular, o estudo da medição de fluxo com elementos intrusivos, tais como placas de orifício, é de grande interesse dado o preço baixo do medidor, e sua boa precisão quando comparada à sua simplicidade de instalação e manutenção. Esta dissertação tem como objetivo o estudo da aplicação de elementos finitos aos escoamentos de fluidos viscosos - via aproximação clássica de Galerkin e Galerkin/mínimos-quadrados (GLS) – com particular ênfase na aproximação das equações de Navier-Stokes incompressível no escoamento newtoniano através de um canal obstruído por uma placa de orifício. Inicialmente, são apresentadas as dificuldades do método de Galerkin clássico na aproximação de escoamentos incompressíveis; ou seja, através da simulação de escoamentos viscosos bem conhecidos - como o escoamento no interior de uma cavidade e através de uma expansão súbita - fica evidenciada a restrição imposta pela condição de Babuška-Brezzi quando da escolha dos subespaços aproximantes do campo de velocidade e pressão Como alternativa às patologias do método de Galerkin clássico, esta dissertação emprega a metodologia de Galerkin/mínimos-quadrados na simulação acima mencionada da placa de orifício, a qual permite o uso de elementos de igual-ordem para aproximar velocidade e pressão e capturar de maneira estável escoamentos sujeitos a altos números de Reynolds. Os testes computacionais realizados se apresentaram fisicamente realistas quando comparados com a literatura e dados experimentais, sendo todos desenvolvidos no Laboratório de Mecânica dos Fluidos Aplicada e Computacional (LAMAC) do Departamento de Engenharia Mecânica da Universidade Federal do Rio Grande do Sul.

Adaptação via movimento de malhas em escoamentos compressíveis

Resumo não disponível.

Projeto, construção e teste de equipamento para ensaios de fadiga de risers flexíveis

Apesar da vida em fadiga de risers flexíveis ser de grande importância na explotação “off-shore”, os dados disponíveis são bastante restritos. A interação entre as suas várias camadas, as inúmeras combinações realizadas entre elas para se obter diferentes características e os recentes avanços em termos de materiais, nos obrigam a estudá-los constantemente. Este trabalho tinha como objetivo principal a construção de um equipamento para possibilitar o estudo da vida em fadiga de risers, através de movimentos de tração combinada com flexão, referentes ao teste “a” da classe II da norma API 17B, bem como todos os outros testes da classe I da referida norma. Foram projetadas a estrutura e todos os demais componentes, que posteriormente foram redimensionados pelos resultados obtidos em simulações numéricas (elementos finitos). Posteriormente, de posse de todos os componentes, realizou-se a montagem do equipamento, seguida de um teste final, que foi acompanhado por uma análise extensométrica dos principais pontos da estrutura. A conclusão básica deste trabalho foi obtida através da comparação das análises por elementos finitos e extensométrica, onde foi constatado que o equipamento estava apto a entrar em operação para a realização de ensaios estáticos ou dinâmicos, observando-se apenas algumas restrições comentadas.

Programa computacional para análise dinâmica de estruturas incluindo amortecimento viscoelástico

Neste trabalho apresenta-se a implementação da matriz de amortecimento viscoelástica para um programa computacional de análise de cascas laminadas de materiais compósitos. A formulação apresentada permite realizar análises dinâmicas de estruturas laminadas com a consideração do efeito do amortecimento para dois modelos diferentes: Kelvin e Zener. A matriz de amortecimento foi implementada de duas formas: proporcional à massa ou proporcional à rigidez. A equação do movimento do sistema dinâmico foi resolvida utilizando-se o método de Newmark para integração direta. Para o modelo Zener foi desenvolvida uma análise para um elemento com 1 grau de liberdade. Apresentam-se exemplos de aplicações da formulação para modelos viscosos, implementadas no programa de elementos finitos, submetidos a diferentes tipos de carregamentos, como carga distribuída e cargas de impacto com diferentes tipos de excitações. Comparações entre o comportamento dos modelos Kelvin e Zener foram realizadas para validar os resultados obtidos.

Análise de escoamentos incompressíveis utilizando simulação de grandes escalas e adaptação de malhas

No presente estudo, são apresentadas soluções numéricas de problemas de Engenharia, na área de Dinâmica dos Fluidos Computacional, envolvendo fluidos viscosos, em escoamentos incompressíveis, isotérmicos e não isotérmicos, em regime laminar e turbulento, podendo envolver transporte de massa. Os principais objetivos deste trabalho são a formulação e a aplicação de uma estratégia de adaptação automática de malhas e a inclusão de modelos de viscosidade turbulenta, integrados com um algoritmo utilizado para simular escoamentos de fluidos viscosos bi e tridimensionais, no contexto de malhas não estruturadas. O estudo é dirigido no sentido de aumentar o conhecimento a respeito das estruturas de escoamentos turbulentos e de estudar os efeitos físicos no transporte de quantidades escalares propiciando, através de técnicas de adaptação automática de malhas, a obtenção de soluções numéricas precisas a um custo computacional otimizado. As equações de conservação de massa, de balanço de quantidade de movimento e de quantidade escalar filtradas são utilizadas para simular as grandes escalas de escoamentos turbulentos e, para representar as escalas submalha, são utilizados dois modelos de viscosidade turbulenta: o modelo de Smagorinsky clássico e o modelo dinâmico. Para obter soluções numéricas com precisão, é desenvolvida e implementada uma estratégia de adaptação automática de malhas, a qual é realizada simultaneamente e interativamente com a obtenção da solução. O estudo do comportamento da solução numérica é fundamentado em indicadores de erro, com o propósito de mapear as regiões onde certos fenômenos físicos do escoamento ocorrem com maior intensidade e de aplicar nestas regiões um esquema de adaptação de malhas. A adaptação é constituída por processos de refinamento/desrefinamento e por um processo de suavização laplaciana. Os procedimentos para a implementação dos modelos de viscosidade turbulenta e a estratégia de adaptação automática de malhas são incorporados ao código computacional de elementos finitos tridimensionais, o qual utiliza elementos tetraédricos lineares. Aplicações de escoamentos de fluidos viscosos, incompressíveis, isotérmicos e não isotérmicos em regime laminar e turbulento são simuladas e os resultados são apresentados e comparados com os obtidos numérica ou experimentalmente por outros autores.

Túneis rodoviários antigos da Madeira

Este trabalho pretende contribuir para um melhor conhecimento sobre os túneis rodoviários antigos existentes na ilha da Madeira. Para tal foi efetuado um levantamento dos túneis construídos entre as décadas de 40 e 80, num total de 28 túneis. Os primeiros túneis rodoviários construídos na década de 40 caracterizam-se por uma geometria retangular ou quadrada. Com passar das décadas e com a evolução ao nível das tecnologias utilizadas, os túneis evoluíram para secções em arco de maior dimensão. Constata-se que estes foram escavados em maciços rochosos com boas características geotécnicas, pelo que eram construídos sem qualquer revestimento, recorrendo-se à execução de revestimentos apenas em zonas pontuais. Para efetuar uma análise paramétrica bidimensional foi utilizado o método dos elementos finitos, onde foram adotadas secções de escavação e materiais com parâmetros de resistência considerados representativos dos túneis em estudo. Este software permitiu analisar o comportamento do maciço com o faseamento construtivo utilizado na escavação destes túneis. Através da caracterização realizada verifica-se que os túneis antigos apresentam-se, em geral, em bom estado de conservação, pelo que a grande maioria deles encontram-se em funcionamento e atualmente inseridos em vias alternativas utilizadas essencialmente por turistas.

Modelação e estudo paramétrico do comportamento de vigas de betão armado com Varões de GFRP

A introdução de um novo material para a execução de estruturas, implica um amplo estudo sobre o mesmo, para que seja aceite como um material de utilização, ou solução, corrente na engenharia civil. Os varões plásticos reforçados com fibra de vidro, GFRP,apresentam-se com diversas vantagens para uma solução de substituição das tradicionais armaduras de aço no betão armado. Um Estado da Arte foi apurado para determinar o ponto de situação do conhecimento e dos estudos realizados por vários autores, desde o início das pesquisas até aos mais recentes estudos disponíveis. Neste seguimento, o presente estudo incide na análise do comportamento à flexão de vigas contínuas de betão armado com varões de GFRP, submetidas a um carregamento distribuído. O trabalho realiza-se com base em 18 modelos de elementos finitos para uma análise paramétrica, sendo os parâmetros estudados, o vão, a percentagem de armadura de flexão e a esbelteza da secção. Para validação do modelo foram, modeladas 7 vigas reais ensaiadas no âmbito de trabalhos de investigação de dois autores Os resultados das modelações das análises paramétricas foram comparados com uma análise elástica e uma análise plástica, utilizando a metodologia do ACI 440.1R-06, para a carga última de dimensionamento. Outras comparações foram feitas através dos limites de deformação do EN 1992-1-1, 2004. Analisa-se a distribuição de esforços, não elásticos, devido à perda de rigidez da viga. Os modelos das vigas mais esbeltas e com menor percentagem de armadura mostraram maior concordância com a análise realizada. As deformações comparadas com o limite L/250 mostram concordância de resultados.

Pré-dimensionamento de vigas de betão armado de acordo com critérios de ductilidade

O dimensionamento de uma estrutura de betão armado depende em grande parte do tipo de metodologia de cálculo utilizado, sendo a análise não linear a que melhor reflete o comportamento da estrutura. Neste tipo de análise é preciso ter em conta alguns aspetos relativos à ductilidade como, a abertura de fendas, a formação de rótulas plásticas e as relações não lineares do material. Entende-se por ductilidade a capacidade de deformação da estrutura sem perda significativa de resistência após atingida a sua capacidade elástica. Esta pode ser quantificada através de diversos indicadores. No caso das vigas é corrente o uso do índice de ductilidade de flecha, que pode ser útil no pré-dimensionamento estrutural pois pode permitir aos projetistas pré-dimensionar vigas segundo um critério de ductilidade. Assim realizou-se uma pesquisa de maneira a identificar quais os fatores mais importantes na ductilidade de uma viga e através do programa cálculo de elementos finitos DIANA 9.4 simularam-se 162 vigas variando a percentagem de armadura mecânica total, o comprimento do vão, a altura útil da secção, a redistribuição de momentos e a classe de betão. Observou-se que o aumento da percentagem de armadura mecânica total origina uma diminuição no índice de ductilidade de cerca de 25% a cada incremento de 0.1 de percentagem de armadura mecânica total, a esbelteza provoca um aumento de 5% a cada 5 unidades no betão C30 e de 2 a 3% no betão C70 e a redistribuição de momentos provocou um aumento na ductilidade. Os resultados mostram que os fatores que mais influenciam na ductilidade das vigas são a percentagem mecânica de armadura total e a classe de resistência do betão, no entanto existem outras que deverão sofrer maior aprofundamento em trabalhos futuros como a utilização de armadura à compressão e os níveis de confinamento do betão.