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

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

Avaliação de ligações soldadas de perfis tubulares em estruturas de aço através do método dos elementos finitos.

A necessidade de se obter soluções para os diversos problemas estruturais na engenharia, associado ao amplo crescimento e difusão da construção metálica fez surgir, por volta dos anos 60, os perfis tubulares. As ligações soldadas entre perfis tubulares são as ligações de maior difusão nas construções com este tipo de perfil, sendo necessário portanto, estudos mais detalhados que tornem a utilização de ligações soldadas mais seguras e otimizadas. O Eurocode 3 possui um item específico para as ligações entre perfis tubulares, em especial as ligações soldadas devido a suas particularidades. Infelizmente, a norma brasileira NBR8800 não aborda as ligações envolvendo perfis tubulares. Para a execução da análise numérica no presente trabalho, faz-se necessário uma revisão bibliográfica. Com estes resultados, foram feitas modelagens de ligações tubulares com o método dos elementos finitos de forma a otimizar os modelos a serem utilizados em uma análise paramétrica futura. Desta forma, foram desenvolvidos dois modelos numéricos, um considerando ligação tipo T entre perfis tubulares quadrados e outro para uma ligação tipo K entre perfis circulares. Estes modelos foram caracterizados através de elementos de casca com seis graus de liberdade por nó considerando-se adicionalmente o efeito de membrana. A análise não-linear realizada considerou a não-linearidade do material através do critério de plastificação de Von Mises através de uma lei constitutiva tensão versus deformação bi-linear e a não-linearidade geométrica foi atribuída através da formulação de Lagrange atualizado. Dentre as principais conclusões obtidas no presente trabalho, pode-se citar que os resultados para as ligações tipo T, o Eurocode 3, fornece resultados que precisam ser observados com cautela. Todavia, para as ligações do tipo K, os resultados numéricos mostraram-se sempre inferiores aos valores através do Eurocode 3, representando um dimensionamento a favor da segurança.

Influência do comportamento semi-rígido de placas de base e de ligações viga-coluna na resposta dinâmica de pórticos de aço.

Tradicionalmente, na análise e dimensionamento de estruturas de aço, assume-se que as ligações viga-coluna são rígidas ou flexíveis (rotuladas). Por outro lado, é de conhecimento geral que a grande maioria das ligações viga-coluna apresenta um comportamento intermediário, ou seja, semi-rígido. Inúmeros trabalhos de pesquisa têm sido desenvolvidos nos últimos vinte e cinco anos, de forma a estudar o comportamento desse tipo de ligação. Um dos principais objetivos desta investigação é o de propor uma metodologia de análise que represente de forma apropriada a influência do comportamento semi-rígido de placas de base e de ligações viga-coluna, sobre a resposta dinâmica (linear e não-linear) de estruturas de aço. Outra contribuição desta dissertação diz respeito à investigação do comportamento dinâmico (linear e não-linear) de pórticos de aço, a partir da consideração de ligações viga-coluna simétricas e não-simétricas e especialmente das placas de base. A análise estrutural é desenvolvida com base no emprego do programa de elementos finitos ANSYS [27]. Nos modelos em elementos finitos foram considerados os efeitos de não-linearidade geométrica (efeitos de segunda ordem), o comportamento não-linear das placas de base e das ligações viga-coluna e, bem como, o efeito de histerese que ocorre quando a estrutura é submetida a cargas cíclicas. Os resultados alcançados indicaram que o fenômeno físico da ressonância não ocorre no que se refere à resposta dinâmica dos modelos semi-rígidos não-lineares. A ressonância não ocorre na resposta dos modelos devido ao fato de que, na análise dinâmica não-linear, o efeito de histerese presente nas ligações (placas de base e viga-coluna), essencialmente com comportamento não-linear, provoca um amortecimento na resposta dinâmica da estrutura.

Estudo da quantificação de incertezas para o problema de contaminação de meios porosos heterogêneos

As técnicas de injeção de traçadores têm sido amplamente utilizadas na investigação de escoamentos em meios porosos, principalmente em problemas envolvendo a simulação numérica de escoamentos miscíveis em reservatórios de petróleo e o transporte de contaminantes em aquíferos. Reservatórios subterrâneos são em geral heterogêneos e podem apresentar variações significativas das suas propriedades em várias escalas de comprimento. Estas variações espaciais são incorporadas às equações que governam o escoamento no interior do meio poroso por meio de campos aleatórios. Estes campos podem prover uma descrição das heterogeneidades da formação subterrânea nos casos onde o conhecimento geológico não fornece o detalhamento necessário para a predição determinística do escoamento através do meio poroso. Nesta tese é empregado um modelo lognormal para o campo de permeabilidades a fim de reproduzir-se a distribuição de permeabilidades do meio real, e a geração numérica destes campos aleatórios é feita pelo método da Soma Sucessiva de Campos Gaussianos Independentes (SSCGI). O objetivo principal deste trabalho é o estudo da quantificação de incertezas para o problema inverso do transporte de um traçador em um meio poroso heterogêneo empregando uma abordagem Bayesiana para a atualização dos campos de permeabilidades, baseada na medição dos valores da concentração espacial do traçador em tempos específicos. Um método do tipo Markov Chain Monte Carlo a dois estágios é utilizado na amostragem da distribuição de probabilidade a posteriori e a cadeia de Markov é construída a partir da reconstrução aleatória dos campos de permeabilidades. Na resolução do problema de pressão-velocidade que governa o escoamento empregase um método do tipo Elementos Finitos Mistos adequado para o cálculo acurado dos fluxos em campos de permeabilidades heterogêneos e uma abordagem Lagrangiana, o método Forward Integral Tracking (FIT), é utilizada na simulação numérica do problema do transporte do traçador. Resultados numéricos são obtidos e apresentados para um conjunto de realizações amostrais dos campos de permeabilidades.

Avaliação do efeito da interação solo estrutura sobre o comportamento estrutural de edificações em aço e mistas (aço concreto).

O conhecimento do tema da interação solo-estrutura permite que as edificações sejam projetadas de maneira mais realista com o comportamento físico. Há décadas atrás seria inviável um dimensionamento considerando a deformabilidade do solo de fundação, e as estruturas eram dimensionadas adotando-se as fundações como indeslocáveis, sob uma base rígida. Essa consideração conduz a respostas estruturais inadequadas, por vezes comprometendo a segurança e a estabilidade do conjunto estrutural. Atualmente, o avanço tecnológico permite a obtenção de resultados de milhões de cálculos matemáticos em questões de segundos, podendo-se obter soluções mais arrojadas e dinâmicas, facilitando o dimensionamento estrutural de forma a atender ao novo padrão de construção e arquitetura. A relevância de tal assunto motivou a análise numérica de um edifício de 4 pavimentos de estrutura mista (aço-concreto), considerando o efeito da interação solo-estrutura. As análises foram conduzidas com o programa ANSYS, de elementos finitos, substituindo-se os apoios indeslocáveis por molas discretas, lineares e elásticas, com rigidez equivalente ao solo, conforme hipótese de Winkler. Os parâmetros dos solos de fundação foram adquiridos a partir de correlações empíricas existentes na literatura e da utilização do programa computacional PLAXIS para a determinação das constantes elásticas das molas. Neste trabalho, foram comparados os resultados de reações verticais, esforços normais, momentos fletores e deslocamentos da estrutura obtidos pelo modelo clássico de projeto, que considera apoios indeslocáveis, e pelo modelo de Winkler, que considera a interação solo-estrutura. As análises foram executadas para seis diferentes tipos de solos argilosos, siltosos e arenosos. Os resultados obtidos mostraram claramente a redistribuição dos momentos fletores, esforços normais e reações verticais nos pilares com diferenças significativas para os pilares de canto e periféricos. Observou-se uma tendência de alívio dos esforços nos pilares mais solicitados, adotando a estrutura assentada em uma base rígida. As análises ressaltaram a relevância da interação solo-estrutura, com efeitos provenientes do rearranjo do solo de fundação nos elementos estruturais das edificações.

Influência da Velocidade, do Espaçamento e do Número de Veículos sobre a Resposta Dinâmica de Pontes Rodoviárias de Concreto Armado.

Nesta dissertação são avaliados os efeitos dinâmicos provenientes da travessia de comboios de veículos sobre o tabuleiro irregular de obras de arte rodoviárias de concreto armado. O modelo matemático empregado para simular o comportamento do sistema veículo-ponte considera a participação da massa e da rigidez das viaturas na definição das freqüências do sistema e, conseqüentemente, a força de interação entre os veículos e a ponte é afetada pela flexibilidade desta. A ponte é modelada a partir do emprego de elementos finitos de barra unidimensionais e discretizado com massas concentradas e flexibilidade distribuída. O modelo de veículo empregado baseia-se no veículo TB-12 preconizado pela norma brasileira NBR 7188. Este veículo é simulado por sistemas de massas, molas e amortecedores sendo descrito por graus de liberdade à translação e rotação no plano. As irregularidades da pista são definidas por um modelo não-determinístico com base na densidade espectral do pavimento. O carregamento sobre a ponte é constituído por sucessões de veículos deslocando-se com velocidade constante sobre a obra. Devido à própria natureza das irregularidades da pista e do comboio de veículos, atenção especial é concentrada na fase permanente da resposta do sistema. São estudadas as respostas de dois modelos estruturais existentes, com base em tabuleiros isostáticos, em concreto armado, com e sem balanços, em seção do tipo T e duplo T, respectivamente, em termos de deslocamentos e esforços nas seções onde ocorrem os efeitos máximos. As conclusões do trabalho versam sobre a influência da velocidade, espaçamento e do número de veículos, referentes a situações distintas de carregamento, no que tange a resposta dinâmica das pontes rodoviárias de concreto armado. A magnitude dos efeitos dinâmicos associados à interação dos veículos com o pavimento irregular também é investigada.

Métodos sem malha: aplicações do Método de Galerkin sem elementos e do Método de Interpolação de Ponto em casos estruturais.

Apesar de serem intensamente estudados em muitos países que caminham na vanguarda do conhecimento, os métodos sem malha ainda são pouco explorados pelas universidades brasileiras. De modo a gerar uma maior difusão ou, para a maioria, fazer sua introdução, esta dissertação objetiva efetuar o entendimento dos métodos sem malha baseando-se em aplicações atinentes à mecânica dos sólidos. Para tanto, são apresentados os conceitos primários dos métodos sem malha e o seu desenvolvimento histórico desde sua origem no método smooth particle hydrodynamic até o método da partição da unidade, sua forma mais abrangente. Dentro deste contexto, foi investigada detalhadamente a forma mais tradicional dos métodos sem malha: o método de Galerkin sem elementos, e também um método diferenciado: o método de interpolação de ponto. Assim, por meio de aplicações em análises de barras e chapas em estado plano de tensão, são apresentadas as características, virtudes e deficiências desses métodos em comparação aos métodos tradicionais, como o método dos elementos finitos. É realizado ainda um estudo em uma importante área de aplicação dos métodos sem malha, a mecânica da fratura, buscando compreender como é efetuada a representação computacional da trinca, com especialidade, por meio dos critérios de visibilidade e de difração. Utilizando-se esses critérios e os conceitos da mecânica da fratura, é calculado o fator de intensidade de tensão através do conceito da integral J.

Elementos finitos em fluidos dominados pelo fenômeno de advecção: um método semi-Lagrangeano.

Os escoamentos altamente convectivos representam um desafio na simulação pelo método de elementos finitos. Com a solução de elementos finitos de Galerkin para escoamentos incompressíveis, a matriz associada ao termo convectivo é não simétrica, e portanto, a propiedade de aproximação ótima é perdida. Na prática as soluções apresentam oscilações espúrias. Muitos métodos foram desenvolvidos com o fim de resolver esse problema. Neste trabalho apresentamos um método semi- Lagrangeano, o qual é implicitamente um método do tipo upwind, que portanto resolve o problema anterior, e comparamos o desempenho do método na solução das equações de convecção-difusão e Navier-Stokes incompressível com o Streamline Upwind Petrov Galerkin (SUPG), um método estabilizador de reconhecido desempenho. No SUPG, as funções de forma e de teste são tomadas em espaços diferentes, criando um efeito tal que as oscilações espúrias são drasticamente atenuadas. O método semi-Lagrangeano é um método de fator de integração, no qual o fator é um operador de convecção que se desloca para um sistema de coordenadas móveis no fluido, mas restabelece o sistema de coordenadas Lagrangeanas depois de cada passo de tempo. Isto prevê estabilidade e a possibilidade de utilizar passos de tempo maiores.Existem muitos trabalhos na literatura analisando métodos estabilizadores, mas não assim com o método semi-Lagrangeano, o que representa a contribuição principal deste trabalho: reconhecer as virtudes e as fraquezas do método semi-Lagrangeano em escoamentos dominados pelo fenômeno de convecção.

Análise de critérios de dimensionamento de carga em dutos de concreto instalados em vala na situação de recobrimento mínimo.

Perante diversas situações da engenharia são utilizadas formulações empíricas de dimensionamento baseadas em dados de campo e experiência profissional que definem muito o caráter subjetivo da metodologia padrão de projeto. O presente trabalho de pesquisa aborda os diversos métodos de obtenção dos esforços gerados em dutos enterrados submetidos a cargas dinâmicas e estáticas e sua posterior reavaliação através de modelagem numérica com o programa Plaxis 3D. Os métodos analíticos não convencionais foram comparados com o método padrão de cálculo sendo que o mesmo demonstrou ter uma boa precisão mesmo sem considerar outros fatores importantes como a parcela de resistência devida à coesão do solo e sua deformabilidade. A modelagem numérica demonstrou o conservadorismo do método de Marston e o subdmensionamento do espraiamento em prisma devido aos efeitos locais ocasionados pela adoção do recobrimento mínimo e sobrecarga dinâmica elevada. Também se observou, através da modelagem 3D, que a utilização dos dois métodos clássicos favorecem a obtenção de resultados dentro da razoabilidade.Verificou-se também, como resultado desta pesquisa, que a proposta de um método clássico modificado permite uma melhor aproximação da carga que atinge o duto.

Estabilidade de talude com efeito sísmico a partir dos métodos de equilíbrio limite e de elementos finitos.

A Costa Peruana apresenta alta atividade sísmica, tornando imprescindível a execução de análises de estabilidade que considerem os eventos sísmicos. Com o desenvolvimento de novas ferramentas numéricas, a análise dinâmica está se tornando cada vez mais importante e usual na fase de projeto, não se justificando mais a execução de análises estáticas em locais tão vulneráveis a sismos. A presente dissertação tem como objetivo realizar a análise da estabilidade dos Penhascos de Lima, Peru, considerando quatro taludes distintos, em uma região bastante afetada por abalos sísmicos. As análises foram executadas pelos métodos de equilíbrio limite, elementos finitos e pseudo-estático, buscando-se uma comparação entre os diferentes métodos. O trabalho apresenta uma descrição completa dos taludes em estudo, uma avaliação da condição sismológica da região e finalmente a estabilidade dos Penhascos de Lima, fazendo uso dos programas computacionais Slide (método do equilíbrio limite) e Plaxis (método dos elementos finitos) em 2D. Os resultados mostraram que os três métodos adotados forneceram fatores de segurança compatíveis, principalmente quando se considera perfis menos estratificados. Para perfis homogêneos, as diferenças obtidas foram da ordem de 0,5 a 1,0 %. As análises ressaltaram a importância de considerar a condição dinâmica, e mostraram-se bastante sensíveis aos valores de carga sísmica adotado.

Análise das deformações dos aterros que compõem a obra do Arco Metropolitano do Rio de Janeiro.

O presente trabalho tem como objetivo a análise numérica do comportamento dos aterros instrumentados que compõem a obra do Arco Metropolitano do Rio de Janeiro. Os resultados da instrumentação de campo serão retroanalisados, juntamente com dados disponíveis na literatura, buscando-se a definição de parâmetros confiáveis, representativos do comportamento da argila compressível da região em estudo. O Arco Metropolitano do Rio de Janeiro é uma rodovia projetada que servirá como ligação entre as principais rodovias que cortam o município do Rio de Janeiro. Devido a presença de grandes espessuras de solo compressível em alguns trechos da região, cortados pelo traçado da rodovia, instrumentos de monitoramento, como placas de recalque e inclinômetros, foram utilizados para avaliar os deslocamentos verticais e horizontais dos aterros durante o processo construtivo. Para este trabalho foram selecionados trechos de aterros representativos, devido à magnitude dos recalques, qualidade do resultado da instrumentação e diferentes métodos construtivos. A partir da análise dos parâmetros de projeto e dos parâmetros encontrados na literatura, procede-se à simulação numérica do processo construtivo dos aterros selecionados com o programa PLAXIS, de elementos finitos, através de modelagem bidimensional. Os resultados numéricos são confrontados com a instrumentação de campo (placas de recalque e inclinômetros) e com os resultados de previsões teóricas (teoria de adensamento unidimensional). Os resultados comprovaram que a modelagem numérica mostrou-se uma ferramenta adequada para a previsão dos recalques totais, tempos de adensamento e ganho de resistência ao longo do tempo. A retroanálise do comportamento de aterros sobre solos moles permite a reavaliação das premissas de projeto, uma vez que as limitações das teorias de análise e a dificuldade na seleção de parâmetros, muitas vezes acarretam em estimativas de recalque incoerentes com as observações de campo.

Morfologia e topografia de superfícies de aço inoxidável duplex UNS S32205, submetidas à erosão por impactos de partículas de alumina em fluxo de ar.

Os aços inoxidáveis duplex vêm sendo cada vez mais utilizados em aplicações que demandam maior resistência à corrosão que os aços inoxidáveis tradicionais. Suas aplicações em estruturas ao ar livre, ou submersas, motivam estudos sobre seu desempenho em variadas situações, como no caso de exposição ao desgaste erosivo.O presente trabalho estuda os efeitos da erosão por impacto de partículas sólidas em superfícies de amostras de um aço inoxidável duplex UNS S32205, tanto do ponto de vista morfológico quanto do topográfico. Superfícies de amostras polidas e posteriormente submetidas a ataques erosivos com partículas de alumina em fluxo de ar, foram examinadas por microscopia eletrônica de varredura, para verificação das formas de desgaste obtidas, e submetidas a ensaios de rugosidade, para determinação do perfil e avaliação da evolução do processo erosivo. Foi feito um exame mais detalhado de formatos de crateras, para comparação com mecanismos propostos em modelos descritos na literatura e com resultados de simulações pelo método dos elementos finitos. Nesse aspecto, o método utilizado mostrou-se tão eficaz quanto aos métodos tradicionais de estudos balísticos com uma só partícula, além de mais abrangente e econômico. As medidas de rugosidade superficial e as informações estatísticas derivadas sinalizaram a possibilidade de utilização dessa ferramenta para avaliar as condições de efetiva perda de massa, que caracterizam a erosão.

Aplicação da equação de Poisson-Boltzmann modificada em sistemas biológicos: análise da partição iônica em um eritrócito

Neste trabalho, a partição iônica e o potencial de membrana em um eritrócito são analisados via equação de Poisson-Boltzmann modificada, considerando as interações não eletrostáticas presentes entre os íons e macromoléculas, assim como, o potencial β. Este potencial é atribuído à diferença de potencial químico de referência entre os meios intracelular e extracelular e ao transporte ativo de íons. O potencial de Gibbs-Donnan via equação de Poisson-Boltzmann na presença de carga fixa em um sistema contendo uma membrana semipermeável também é estudado. O método de aproximação paraboloide em elementos finitos em um sistema estacionário e unidimensionalé aplicado para resolver a equação de Poisson-Boltzmann em coordenadas cartesianas e esféricas. O parâmetro de dispersão relativo às interações não eletrostáticas écalculado via teoria de Lifshitz. Os resultados em relação ao potencial de Gibbs-Donnan mostram-se adequados, podendo ser calculado pela equação de Poisson-Boltzmann. No sistema contendo um eritrócito, quando o potencial β é considerado igual a zero, não se verifica a diferença iônica observada experimentalmente entre os meios intracelular e extracelular. Dessa forma, os potenciais não eletrostáticos calculados via teoria de Lifshitz têm apenas uma pequena influência no que se refere à alta concentração de íon K+ no meio intracelular em relação ao íon Na+

A new Lagrangian-Eulerian shell-fluid coupling algorithm based on level sets

We propose a computational method for the coupled simulation of a compressible flow interacting with a thin-shell structure undergoing large deformations. An Eulerian finite volume formulation is adopted for the fluid and a Lagrangian formulation based on subdivision finite elements is adopted for the shell response. The coupling between the fluid and the solid response is achieved via a novel approach based on level sets. The basic approach furnishes a general algorithm for coupling Lagrangian shell solvers with Cartesian grid based Eulerian fluid solvers. The efficiency and robustness of the proposed approach is demonstrated with a airbag deployment simulation. It bears emphasis that in the proposed approach the solid and the fluid components as well as their coupled interaction are considered in full detail and modeled with an equivalent level of fidelity without any oversimplifying assumptions or bias towards a particular physical aspect of the problem.