Otimização topológica de estruturas bidimensionais contínuas submetidas a restrições de flexibilidade e tensão

Neste trabalho é resolvido o problema da minimização do volume de estruturas bidimensionais contínuas submetidas a restrições sobre a flexibilidade (trabalho das forças externas) e sobre as tensões, utilizando a técnica chamada otimização topológica, que visa encontrar a melhor distribuição de material dentro de um domínio de projeto pré-estabelecido. As equações de equilíbrio são resolvidas através do método dos elementos finitos, discretizando a geometria e aproximando o campo de deslocamentos. Dessa forma, essas equações diferenciais são transformadas em um sistema de equações lineares, obtendo como resposta os deslocamentos nodais de cada elemento. A distribuição de material é discretizada como uma densidade fictícia constante por elemento finito. Esta densidade define um material isotrópico poroso de uma seqüência pré-estabelecida (SIMP). A otimização é feita através da Programação Linear Seqüencial. Para tal, a função objetivo e as restrições são sucessivamente linearizadas por expansão em Série de Taylor. A análise de sensibilidade para a restrição de flexibilidade é resolvida utilizando o cálculo da sensibilidade analítico adaptado para elementos finitos de elasticidade plana. Quando as restrições consideradas são as tensões, o problema torna-se mais complexo. Diferente da flexibilidade, que é uma restrição global, cada elemento finito deve ter sua tensão controlada. A tensão de Von Mises é o critério de falha considerado, cuja sensibilidade foi calculada de acordo com a metodologia empregada por Duysinx e Bendsøe [Duysinx e Bendsøe, 1998] Problemas como a instabilidade de tabuleiro e dependência da malha sempre aparecem na otimização topológica de estruturas contínuas. A fim de minimizar seus efeitos, um filtro de vizinhança foi implementado, restringindo a variação da densidade entre elementos adjacentes. Restrições sobre as tensões causam um problema adicional, conhecido como singularidade das tensões, fazendo com que os algoritmos não convirjam para o mínimo global. Para contornar essa situação, é empregada uma técnica matemática de perturbação visando modificar o espaço onde se encontra a solução, de forma que o mínimo global possa ser encontrado. Esse método desenvolvido por Cheng e Guo [Cheng e Guo, 1997] é conhecido por relaxação-ε e foi implementado nesse trabalho.

Simulação numérica da interação de escoamentos tridimensionais de fluidos compressíveis e incompressíveis e estruturas deformáveis usando o método de elementos finitos

Neste trabalho apresenta-se um algoritmo para a simulação de problemas tridimensionais de interação fluido-estrutura utilizando a técnica de elementos finitos. Um esquema de Taylor-Galerkin de dois passos e elementos tetraédricos lineares são empregados para o fluido, que pode ser compressível ou incompressível. É adotada uma formulação lagrangeana-euleriana arbitrária (ALE), compatível com o movimento da interface fluidoestrutura. Um método ftacionado de correção de velocidade é utilizado para os fluidos incompressíveis. A estrutura é analisada usando elementos triangulares com três nós e seis graus de liberdade por nó (três componentes de deslocamentos e três componentes de rotação). Os efeitos da não-linearidade geométrica são incluídos. O método de Newmark é empregado para integrar no tempo as equações dinâmicas de equilíbrio, usando-se uma descrição lagrangeana atualizada. O sistema de equações alge'bricas é solucionado através do método dos gradientes conjugados e o sistema não-linear, resultante de deslocamentos e rotacões finitas da estrutura, é solucionado com um esquema incremental-iterativo. O código é otimizado para aproveitar as vantagens do processamento vetorial.

Otimização em dois níveis da orientação e da topologia de cascas laminadas

Este trabalho é uma contribuição para o conhecimento de metodologias de projeto de estruturas de material composto, aplicando métodos de otimização estrutural a cascas laminadas e apresentando uma estratégia em dois níveis. No primeiro nível é realizada a minimização da flexibilidade da estrutura, tendo como variável de projeto a orientação de cada lâmina da estrutura. Utiliza-se Programação Linear Seqüencial (SLP) e direção de tensão principal para otimização da orientação. No segundo nível minimiza-se o volume de cada lâmina, usando a flexibilidade total da estrutura como restrição e a densidade relativa como variável de projeto, também através de SLP. Para evitar aparecimento de áreas com densidades intermediárias, utiliza-se um Método de Continuação, dividindo o nível de otimização topológica em duas ou mais etapas. As formulações desenvolvidas permitem a solução de problemas com múltiplos casos de carregamento. Para a solução da equação de equilíbrio de casca laminada, utiliza-se um elemento finito de casca degenerado de oito nós com integração explícita na direção da espessura. A implementação desse elemento é feita de modo a facilitar a obtenção das derivadas da matriz de rigidez, necessárias na linearização das funções objetivo e restrições. Evita-se assim o uso de derivadas numéricas. Resultados para vários tipos de estrutura são apresentados, incluindo comparações entre diferentes carregamentos, condições de contorno, número de lâminas, espessuras, etc. As soluções obtidas, formas de análise e possíveis aplicações são discutidas.

Comportamento mecânico e hidráulico de misturas de solo-cimento-bentonita para a aplicação em barreiras verticais de contenção de contaminantes

O presente trabalho tem como objetivo investigar o comportamento hidráulico e mecânico de misturas de solo-bentonita (SB) e de solo-cimento-bentonita (SCB), bem como analisar a influência da variação do teor de bentonita e cimento com a finalidade de propiciar resultados capazes de auxiliar para um projeto construtivo de barreiras verticais de contaminantes. O programa experimental consistiu na realização de ensaios de compressão não-confinada, ensaios triaxiais convencionais adensados não-drenados (CIU) e ensaios de condutividade hidráulica a fim de se estudar o comportamento das misturas de SB e SCB em termos de tensão-deformação-resistência e permeabilidade. Foram analisados diferentes fatores água-cimento, teores de umidade, tensões confinantes. Também foi avaliado o comportamento desses materiais através da definição de três valores distintos de teor de bentonita nas amostras de SCB. Já para as amostras de SB foram avaliados dois teores de bentonita. Em um primeiro instante, os resultados dos ensaios realizados para cada tipo de mistura foram analisados separadamente. A análise global dos resultados permitiu identificar as alterações provocadas na resistência, deformabilidade e condutividade hidráulica pela variação dos parâmetros. Os resultados dos ensaios de compressão não-confinada demonstraram que a variação do teor de umidade não implicou na variação de resistência ao cisalhamento de pico para amostras com mesmo fator águacimento. Os ensaios CIU indicaram que o aumento do teor de bentonita nas amostras acarretou em um aumento da resistência ao cisalhamento de pico e da rigidez do material. A condutividade hidráulica das misturas de SB apresentaram valores inferiores em relação às misturas de SCB.

Projeto, construção e validação de transdutores para medição de esforços em semeadoras-adubadoras

A mecanização agrícola tem sido um importante fator que vem contribuindo para atender o aumento na demanda por alimentos. Porém, o desenvolvimento de máquinas e implementos com maior desempenho operacional ainda carecem de uma base de informações tecnológicas sistematizadas para a implementação em novos projetos de produtos. Desta forma, o principal objetivo deste trabalho foi o de implementar um sistema capaz de gerar informações para auxiliar no desenvolvimento de novos projetos de produtos, ou mesmo na melhoria das máquinas e implementos já existentes. Para tanto, foi construído um conjunto de transdutores de força – dinamômetros do tipo anel octogonal estendido, capazes de medir solicitações em duas direções ortogonais e o momento associado a estes dois esforços. Associados a estes transdutores foram utilizados outros recursos para a formação de uma instrumentação embarcada, a fim de possibilitar os registros, em tempo real, dos dados gerados para serem armazenados em bancos de dados para posterior análises e avaliações. Com o sistema implementado foram realizados diversos conjuntos de testes de campo, onde se ensaiou uma semeadora adubadora comercial com seus elementos ativos e/ou seus componentes isolados da máquina. No primeiro teste, foi avaliada uma linha completa de uma semeadora-adubadora, onde foram medidos os esforços no disco de corte da palha, na haste sulcadora para deposição do adubo e na unidade de semeadura. Nos testes seguintes, foram realizadas avaliações em dois modelos de hastes sulcadoras, trabalhando com diversos parâmetros distintos, a fim de comparar os respectivos desempenhos entre os dois modelos estudados. Também, verificaram-se avaliações das solicitações dos esforços de tração para sulcadores de adubo tipo disco duplo e, para tracionar somente os discos de corte. Finalmente, mediram-se os esforços para o arraste do chassi durante o deslocamento em campo. Os resultados dos sensores de força combinados com a instrumentação embarcada foram considerados satisfatórios durante as medições nos testes de campo, visto a sua repetibilidade tendo como referência um sistema de calibrações prévias. Desta forma, confirmouse que os registros dos testes de campo podem servir de referência para a análise de projeto de máquinas semeadoras-adubadoras, uma vez que os mesmos foram obtidos em condições reais durante as operações de campo.

Estudo da instabilidade de um talude rodoviário em São Jerônimo/RS

Esta dissertação de mestrado apresenta um a investigação das causas que levara mum talude rodoviário à ruptura. Esse talude localiza-se no município de São Jerônimo, no Km 7+500 da rodovia RS/470. A ruptura foi progressiva, tendo início logo após a abertura do corte para a execução da rodovia e envolveu um a extensão demais de 100 metros. Foram levantados aspectos geológicos e geomorfológicos da região, indicando tratar-se de uma área localizada na Depressão Periférica, com relevo suave ondulado. Os dois solos que compõe o talude foram descritos e amostrados. Foram realizados ensaios de caracterização e ensaios de cisalhamento para as diversas condições de deformação dos materiais. Foram feitos diversos ensaios de cisalhamento direto: a) com amostras indeformadas, buscando-se determinar parâmetros de resistência para a condição de pico e, b ) ensaios especiais, utilizando a técnica de reversões múltiplas para determinar parâmetros de resistência na condição de grandes deformações. Foram feitos também ensaios de cisalhamento por torção (ring shear), buscando-se parâmetros de resistência para a condição residual. A estabilidade do talude foi verificada para duas geometrias, a inicial e a pós-ruptura, utilizando-se os parâmetros de resistência ao cisalhamento obtidos nos ensaios realizados. A retro-análise da situação inicial indicou que as rupturas iniciais mobilizaram um valor de resistência ao cisalhamento correspondente a um a condição de deformação intermediária entre a condição de pico e a de grandes deformações. A retro-análise da condição atua lindicou que os solos ainda não atingiram a condição residual, sendo que a resistência mobilizada é melhor representada pelos parâmetros na condição de grandes deformações.

Métodos multigrid paralelos em malhas não estruturadas aplicados à simulação de problemas de dinâmica de fluidos computacional e transferência de calor

Fenômenos naturais, tecnológicos e industriais podem, em geral, ser modelados de modo acurado através de equações diferenciais parciais, definidas sobre domínios contínuos que necessitam ser discretizados para serem resolvidos. Dependendo do esquema de discretização utilizado, pode-se gerar sistemas de equações lineares. Esses sistemas são, de modo geral, esparsos e de grande porte, onde as incógnitas podem ser da ordem de milhares, ou até mesmo de milhões. Levando em consideração essas características, o emprego de métodos iterativos é o mais apropriado para a resolução dos sistemas gerados, devido principalmente a sua potencialidade quanto à otimização de armazenamento e eficiência computacional. Uma forma de incrementar o desempenho dos métodos iterativos é empregar uma técnica multigrid. Multigrid são uma classe de métodos que resolvem eficientemente um grande conjunto de equações algébricas através da aceleração da convergência de métodos iterativos. Considerando que a resolução de sistemas de equações de problemas realísticos pode requerer grande capacidade de processamento e de armazenamento, torna-se imprescindível o uso de ambientes computacionais de alto desempenho. Uma das abordagens encontradas na literatura técnica para a resolução de sistemas de equações em paralelo é aquela que emprega métodos de decomposição de domínio (MDDs). Os MDDs são baseados no particionamento do domínio computacional em subdomínios, de modo que a solução global do problema é obtida pela combinação apropriada das soluções obtidas em cada um dos subdomínios Assim, neste trabalho são disponibilizados diferentes métodos de resolução paralela baseado em decomposição de domínio, utilizando técnicas multigrid para a aceleração da solução de sistemas de equações lineares. Para cada método, são apresentados dois estudos de caso visando a validação das implementações. Os estudos de caso abordados são o problema da difusão de calor e o modelo de hidrodinâmica do modelo UnHIDRA. Os métodos implementados mostraram-se altamente paralelizáveis, apresentando bons ganhos de desempenho. Os métodos multigrid mostraram-se eficiente na aceleração dos métodos iterativos, já que métodos que utilizaram esta técnica apresentaram desempenho superior aos métodos que não utilizaram nenhum método de aceleração.

Estudo crítico da metodologia de cálculo para silos prismáticos multicelulares

A armazenagem de grãos é atualmente um importante tópico a ser explorado na região sul do Brasil, e vem se desenvolvendo consideravelmente nos últimos anos. Atualmente a armazenagem é feita através de silos cilíndricos que podem ser de grandes diâmetros ou elevados. Na área da armazenagem é importante implantar novas tecnologias como o caso de silos multicelulares prismáticos que são bastante difundidos nos Estados Unidos. Para estudar estes silos propõe-se analisar um caso real sob a ótica morfológica e dimensional. O estudo baseia-se na comparação entre métodos de análise estrutural analítica convencional com análise numérica através do cálculo por elementos finitos, além de testar as análises de cálculo mais avançadas. Para estudar um caso real buscaram-se informações detalhadas através de desenhos fornecidos por empresas que fabricam estes tipos de silos. Adotando um caso prático como padrão, na seqüência fez-se a análise da interação entre o número de módulos do silo com o objetivo de verificar a influência do número de módulos no dimensionamento dos componentes do silo. Analisou-se o referido projeto utilizando procedimentos clássicos de análise de dimensionamento (métodos analíticos) e compararam-se os mesmos com a análise utilizando o método dos elementos finitos. Adicionalmente executaram-se análises mais avançadas utilizando o método dos elementos finitos como não linearidade geométrica e física, flambagem e análise de freqüências e modos.Após as referidas análises pôde-se concluir que analisando dimensionalmente somente um módulo do silo, pode-se obter resultados que representam todo o conjunto de módulos. Além disso, verificou-se que a análise por elementos finitos é mais apropriada, pois permite menos simplificações de cálculo, apresentando assim resultados mais reais. Quanto às análises adicionais, verificou-se que a carga crítica de flambagem considerando a análise não linear geométrica por elementos finitos é maior do que a carga crítica obtida por análises convencionais.

Comportamento de um solo colapsível artificialmente cimentado

Esta tese apresenta duas técnicas diferentes de preparação de amostras para produzir corpos-de-prova “ideais”, i.e. amostras artificialmente moldadas nas condições desejadas. O desafio fundamental da moldagem de amostras altamente colapsíveis cimentadas artificialmente é conservar a estrutura porosa com elevado índice de vazios enquanto ocorre a cura do cimento. Na primeira técnica isto é alcançado acrescentandose partículas de poliestireno expandido à mistura de solo. A segunda técnica utiliza duas etapas de preparo de forma que a cura parcial do cimento da primeira etapa de mistura já aconteceu quando a segunda etapa for realizada. Ambas as técnicas são utilizadas para produzir amostras colapsíveis que imitam o comportamento de solos altamente colapsíveis naturais. Ao mesmo tempo, as técnicas permitem a variação controlada dos parâmetros fundamentais que podem influenciar comportamento de colapso do solo. O trabalho apresentado nesta tese analisou os parâmetros que caracterizam e controlam o comportamento mecânico de solos colapsíveis. O programa experimental foi executado em amostras colapsíveis cimentadas artificialmente que possuem as características físicas mais relevantes de solos colapsíveis naturais: elevado índice de vazios, baixo nível de cimentação e condições não saturadas antes da inundação. O programa de laboratório foi fundamentado em ensaios de compressão não confinada, ensaios de condutividade hidráulica, testes para obtenção da curva característica do solo, análises de microscopia ótica, ensaios de cisalhamento direto, ensaios triaxiais, ensaios edométricos convencionais e com sucção controlada e ensaio de compressão isotrópica com controle de sucção Com base nos resultados experimentais, foi possível observar a importante influência do índice de vazios inicial e do teor de adição de cimento no processo de degradação da cimentação (desestruturação do solo) e no fenômeno de colapso por inundação. Isto é importante para se compreender o comportamento mecânico fundamental de solos colapsíveis cimentados e permitir o desenvolvimento de modelos constitutivos apropriados. É apresentado um modelo constitutivo (Wheeler, 2004) baseado em solos não saturados incluindo o efeito de cimentação (bonding) e o processo de desestruturação do solo. O modelo sugerido foi analisado qualitativamente com base nos resultados obtidos do programa experimental.

Determinação do comportamento carga-recalque de sapatas em solos residuais a partir de ensaios SPT

O comportamento carga-recalque de sapatas assentes em solos residuais é perfeitamente representado a partir da interpretação de provas de carga. No entanto, a execução de provas de carga em fundações reais ou em placa é indesejada pela maioria dos construtores e contratantes por razões de cronograma e custo. Por outro lado, sondagens de simples reconhecimento (SPT) são de rápida execução, exigidas pela normalização vigente de projeto e execução de fundações (ASSOCIAÇÃO BRASILEIRA DE NORMAS TÉCNICAS, 1996) e largamente utilizados na geotecnia nacional, sendo tal ensaio normalizado pela da norma NBR 6.484/2001. Deste modo, nesta pesquisa, pretende-se estabelecer uma metodologia semi-empírica para estimar o comportamento de sapatas assentes em solos residuais a partir de correlação estatística de resultados de sondagens de simples reconhecimento (NSPT) com a tensão admissível e os recalques sob tensão de trabalho, obtidos a partir de resultados de provas de carga em placa ou fundações reais assentes em solos residuais. Deste modo, a partir de retro-análise de um conjunto de resultados de provas de carga executados em diferentes solos residuais, encontrados na literatura, e o uso de conceitos de estatística, obteve-se diversas correlações que possibilitam estimar o comportamento tensão-recalque deste tipo de solo através dos resultados médios dos ensaios de SPT Como resultados destas correlações foram obtidas equações para a previsão de recalques e tensão admissível de fundações assentes sobre solos residuais, bem como do módulo de elasticidade de tais solos. As equações de estimativa de recalques são dependentes, além dos resultados de SPT, da geometria da fundação e da tensão aplicada. Já as equações para a estimativa da tensão admissível e módulo de elasticidade são dependentes somente dos valores médios do ensaio de SPT. E finalmente, para testar as equações desenvolvidas, foi realizada uma análise “Classe A” no comportamento carga-recalque de três provas de carga realizadas em locais diferentes, mostrando a eficiência da metodologia desenvolvida para a determinação da tensão admissível e dos recalques em sapatas sob a aplicação de carga de trabalho assentes em solos residuais.

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.

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.

Investigação geoambiental de uma área contaminada por resíduos industriais

O presente estudo tem por objetivo apresentar e discutir os resultados de uma investigação geoambiental realizada em uma área contaminada por resíduos industriais, na qual foi depositada borra oleosa ácida. Sondagens de simples reconhecimento (SPT) foram empregadas para a identificação do perfil estratigráfico da área, e ensaios laboratoriais (granulometria e limites de Atterberg) foram utilizados para caracterização das amostras de solo. Os resultados mostraram que o subsolo da área é constituído por um perfil de solo residual composto essencialmente pelas frações argila e silte. Uma série de sondagens de cone resistivo (RCPT) foi realizada na área com o objetivo de detectar a presença de contaminantes no subsolo, através de medidas da resistividade elétrica. As verticais resistivas foram comparadas com os resultados de duas campanhas de geofísica superficial por eletro resistividade, que foram obtidas em trabalhos anteriores de avaliação da contaminação da área. Com o propósito de avaliar a eficiência do cone resistivo na identificação das zonas contaminadas, os ensaios de RCPT foram realizados próximos a poços de monitoramento, onde foram coletadas amostras da água subterrâneas para análises quantitativas da contaminação. O principal efeito da contaminação do meio físico é redução dos parâmetros físico-químicos, pH e resistividade elétrica da água subterrânea, em razão da elevada quantidade de ácido sulfúrico e metais que são liberados pelo resíduo borra ácida depositado na área. Através das verticais resistivas foi possível identificar uma zona altamente contaminada até 5 m de profundidade, caracterizada por valores muito baixos de resistividade. Os dados fornecidos pelo cone resistivo apresentaram concordância satisfatória com as análises químicas e físico-químicas da água subterrânea. Nas zonas superficiais do subsolo (< 5 m de profundidade) os valores de resistividade obtidos pelo cone resistivo foram inferiores àqueles obtidos através da geofísica de superfície. A compatibilidade entre os resultados foi observada para as cotas mais profundas do subsolo.

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.

Resistência ao cisalhamento de solos não saturados : análise experimental e teórica

Um equipamento triaxial convencional originalmente utilizado para ensaios em amostras de solos saturados saturada foi modificado para que fosse possível a realização de ensaios triaxiais em solos não saturados. As principais modificações foram a instalação de uma pedra cerâmica de alto valor de entrada de ar (300 kPa) na nova base da câmara triaxial e uma nova linha de aplicação de pressão para aplicação de pressão de ar na amostra. A técnica de translação de eixos é utilizada para evitar cavitação no sistema de aplicação de pressão. O solo estudado nesta pesquisa consiste de um colúvio de arenito localizado nas encostas da Formação Serra Geral, entre os municípios de Timbé do Sul (SC) e São José dos Ausentes (RS). Ensaios triaxiais convencionais foram realizados nas condições drenado e não drenado em amostras indeformadas para tensões de confinamento iniciais variando de 50 até 500 kPa; ensaios triaxiais convencionais foram realizados na condição drenada em amostras remoldadas para as tensões de confinamento iniciais de 50, 100 e 200kPa. Ensaios triaxiais com sucção controlada foram realizados na condição drenado em amostras indeformadas para as tensões normais líquidas de 50, 100, 150 e 200 kPa em níveis de sucção variando de 25 até 150 kPa. A partir dos resultados dos ensaios triaxiais determinou-se os parâmetros de resistência ao cisalhamento: c’ – intercepto coesivo efetivo; f´ - ângulo de atrito interno efetivo; fb – ângulo de atrito interno que quantifica a contribuição da sucção na resistência ao cisalhamento. O comportamento encontrado foi de uma envoltória de resistência ao cisalhamento não saturada bilinear com acréscimo na resistência ao cisalhamento para baixos níveis de sucção até aproximadamente 75 kPa, e após um decréscimo na resistência ao cisalhamento. A tensão normal líquida não influencia no valor de fb e a sucção não provoca alteração no valor de f´. O valor de fb apresenta valores superiores ao valor de f´, quando se utiliza o critério de ruptura clássico de resistência ao cisalhamento. Alguns critérios de determinação dos valores máximos de ruptura foram aplicados para um melhor entendimento do comportamento resistente deste material.