Análise estática e dinâmica de placas e cascas de materiais compósitos laminados usando elementos finitos hexaédricos de oito nós com integração reduzida

Neste trabalho implementou-se o elemento hexaédrico com um ponto de integração para análise estática e dinâmica de placas e cascas de materiais compósitos laminados com ou sem enrijecedores. O elemento está livre de travamento volumétrico e travamento de cisalhamento, não apresentando modos espúrios. São também incluídos problemas com não-linearidade geométrica A matriz de rigidez e de massa são dadas de forma explícita, reduzindo o tempo computacional, especialmente na análise não-linear. Para evitar o travamento de cisalhamento as componentes de deformações são referidas a um sistema co-rotacional. O travamento volumétrico é também eliminado, já que a parte dilatacional (esférica) da matriz gradiente é avaliada apenas no ponto central do elemento. Para a solução das equações de equilíbrio na análise estática, utilizam-se um método direto baseado na eliminação de Gauss ou um método iterativo de gradientes conjugados com precondicionamento executado através da eliminação incompleta de Choleski. Para a análise dinâmica, as equações de equilíbrio são integradas através do método explícito ou implícito de Taylor-Galerkin ou do método implícito de Newmark. Para análise não-linear utiliza-se o Método Generalizado de Controle dos Deslocamentos. Através de exemplos numéricos demonstra-se a eficiência e o potencial do elemento tridimensional na análise linear e não-linear de placas e cascas de materiais laminados. Os resultados são comparados com trabalhos que utilizam diferentes elementos de placas e cascas.

Análise numérica de problemas de fratura elástica e elastoplástica

Este trabalho apresenta o método dos elementos finitos em conjunto com métodos numéricos especificos para a solução de problemas de fratura. Esta é uma poderosa ferramenta para a análise de fraturas e soluções confiáveis são obtidas para problemas complexos de Engenharia tanto no campo linear como no não-linear. O elemento finito. implementado é do tipo isoparamétrico quadrâtico da família Serendipity. Com dois graus de liberdade por nó, permite discretizar em estado plano de tensão ou deformação estruturas com geometrias bastante variadas. Para a análise linear são implementadas quatro técnicas consagradas para a avaliação do fator de intensidade de tensão no modo I de fratura: extrapolação de doslocamentos (usando malha somente com elementos convencionais e malha mesclada com elementos especiais), taxa de liberação de energia de defermação, extensão virtual da trinca e o método da integral J, descartando-se neste caso a hipótese de descarregamento. A linguagem de programação adotada é o FORTRAN 77. A entrada de dados é feita por intermédio de arquivos previamente preparados. Os resultados obtidos são confrontados com resultados experimentais e computacionais fornecidos por outros programas. Analisam-se placas, estruturas de uso na indústria e simulam-se ensaios como o corpo de prova de flexão em três pontos e o corpo de prova de tensão. compacto.

Análise estática e dinâmica, linear e não-linear geométrica, através de elementos hexaédricos de oito nós com um ponto de integração

Para a análise estática e dinâmica, linear e não-linear de placas, cascas e vigas, implementa-se neste trabalho o elemento hexaédrico com integração reduzida, livre de travamento volumétrico e travamento de cisalhamento e que não apresenta modos espúrios. Na formulação do elemento, utiliza-se apenas um ponto de integração. Desta forma, a matriz de rigidez é dada de forma explícita e o tempo computacional é significativamente reduzido, especialmente em análise não-linear. Os modos espúrios são suprimidos através de um procedimento de estabilização que não exige parâmetros especificados pelo usuário. Para evitar o travamento de cisalhamento, desenvolve-se o vetor de deformações num sistema co-rotacional e remove-se certos termos não constantes nas componentes de deformações de cisalhamento. O travamento volumétrico é resolvido fazendo-se com que a parte dilatacional (esférica) da matriz gradiente seja avaliada apenas no ponto central do elemento. Como a eliminação do travamento de cisalhamento depende de uma abordagem no sistema local, emprega-se um procedimento co-rotacional para obter o incremento de deformação no sistema local e atualizar os vetores de tensões e forças internas na análise não-linear Para a solução das equações de equilíbrio na análise estática, utilizam-se métodos diretos baseados na eliminação de Gauss ou métodos iterativos de Gradientes Conjugados Precondicionado elemento-por-elemento (EBE). Para a análise dinâmica, as equações de equilíbrio são integradas através do método explícito de Taylor-Galerkin ou do método implícito de Newmark. Através de exemplos numéricos demonstra-se a eficiência e o potencial do elemento tridimensional na análise de casca, placas e vigas submetidas a grandes deslocamentos e grande rotações. Os resultados são comparados com trabalhos que utilizam elementos clássicos de placa e casca.

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.

Identificação de danos em estruturas bi-dimensionais via matriz de flexibilidade baseado em um modelo de dano contínuo

No presente trabalho, o modelo de identificação de danos apresentado por Stutz et al. (2005) é utilizado. A contribuição do presente trabalho consiste em avaliar alguns pontos da identificação de danos em vigas e, em seguida, expandir o modelo para identificar danos estruturais em placas. Uma avaliação do comportamento das frequências naturais e da matriz de flexibilidade para uma viga de Euler-Bernoulli simplesmente apoiada, na presença de danos simulados pelo parâmetro de coesão é realizada. Essa análise, permite também o conhecimento das regiões onde há maior sensibilidade ao dano, ajudando a traçar estratégias para melhorar a identificação de danos em regiões que sofrem poucas alterações na presença de falhas estruturais. Comparou-se o comportamento dos dois primeiros modos de vibração da viga simplesmente apoiada na presença de um dano estrutural, com os dois primeiros modos de vibração da estrutura intacta e corrompidos por ruído. Diversos métodos de localização de danos e de otimização são avaliados na tentativa de identificar os danos simulados através do campo de danos proposto por Stutz et al. (2005) na presença de dados ruidosos. Após a apresentação de resultados da identificação de danos obtidos para uma viga de Euler-Bernoulli simplesmente apoiada, uma análise do comportamento das frequências naturais e da matriz de flexibilidade de uma viga de Euler- Bernoulli engastada na presença de danos é apresentada, assim como os resultados de identificação de danos considerando-se diversos cenários e níveis de ruído. Uma importante contribuição do presente trabalho consiste em propor um método de identificação de danos via matriz de flexibilidade onde o campo de defeitos para a placa de Kirchoff é modelado via MEF. Uma análise do comportamento da matriz de flexibilidade devido à presença de danos na placa é apresentada, assim como os resultados numéricos da identificação de danos estruturais com e sem a presença de dados ruidosos. Com a finalidade de reduzir o custo computacional na identificação de danos em estruturas complexas, uma hibridização entre o método de otimização por enxame de particulas (PSO, do inglês, Particle Swarm Optimization) e o método de otimização Levenberg-Marquardt é proposta. Resultados numéricos da hibridização para uma estrutura do tipo placa são apresentados.

Controle activo de estruturas adaptativas com sensores e actuadores piezoeléctricos

Nos últimos anos o estudo de estruturas inteligentes tem atraído vários investigadores devido às suas potenciais vantagens numa larga gama de aplicações, tais como controle de forma, supressão de vibrações, atenuação de ruído e detecção de dano. O uso de materiais “inteligentes” tal como os materiais piezoeléctricos na forma de lâminas ou “patches”, embebidas ou coladas na superfície de estruturas construídas de materiais compósitos, permite assim obter estruturas que por um lado são adaptativas e por outro revelam excelentes propriedades mecânicas, aumentando assim bastante o desempenho e a fiabilidade de sistemas estruturais. Os materiais piezoeléctricos têm a propriedade de gerar uma carga eléctrica sob a acção duma carga mecânica e o reverso, isto é, aplicando um campo eléctrico nos elementos piezoeléctricos da estrutura, esta deforma-se. Neste trabalho, é apresentado um modelo de elementos finitos, baseado na teoria clássica de placas, desenvolvido para a análise do controle activo em estática e dinâmica lineares de estruturas integrando sensores e actuadores piezoeléctrico na forma de lâminas, os quais introduzem um grau de liberdade referente ao potencial eléctrico, por cada camada piezoeléctrica do elemento finito. É utilizado método de Newmark para a solução iterativa das equações de equilíbrio. Apresentam-se os resultados obtidos em três exemplos ilustrativos.

Optimização de estruturas laminadas adaptativas com controlo activo

Neste trabalho apresenta-se um modelo de elementos finitos, baseados na teoria clássica das placas, desenvolvido para a análise de controlo activo em dinâmica de estruturas de tipo placa/casca integrando sensores e actuadores piezoeléctricos. O controlo é iniciado através de uma optimização prévia do núcleo laminado de modo a diminuir a amplitude da vibração. É usado um algoritmo de controlo baseado na ligação entre as lâminas piezoeléctricas sensoras e actuadoras para obter um mecanismo de controlo da resposta dinâmica da estrutura. A resolução por elementos finitos usa um elemento placa/casca triangular plano de 3 nós, e em cuja formulação se introduz um grau de liberdade referente ao potencial eléctrico, por cada camada piezoeléctrica do elemento finito. Apresentam-se os resultados obtidos em dois exemplos ilustrativos.

Análise geometricamente não-linear de estruturas adaptativas piezolaminadas

Neste trabalho apresenta-se um modelo de elementos finitos, baseado na teoria clássica de placas, para a análise linear e não-linear de estruturas do tipo placa/casca integrando sensores e actuadores piezoeléctricos. É usado um simples e eficiente elemento placa/casca triangular plano de 3 nós, e em cuja formulação se introduz um grau de liberdade referente ao potencial eléctrico, por cada camada piezoeléctrica do elemento finito. É utilizada a formulação Lagrangeana actualizada associada à tecnica de Newton - Raphson para a solução iterativa das equações de equilibrio .O modelo pode ser aplicado a cascas piezolaminadas com geometria e carregamento arbitrários. Apresentam-se vários exemplos ilustrativos cujos resultados mostram a eficiencia do modelo proposto.

Placas elasto-electro-magnéticas

Neste trabalho apresenta-se um modelo de elementos finitos, baseado na teoria clássica de placas, para a análise linear e não-linear de estruturas do tipo placa/casca integrando sensores e actuadores piezoeléctricos. É usado um simples e eficiente elemento placa/casca triangular plano de 3 nós, e em cuja formulação se introduz um grau de liberdade referente ao potencial eléctrico, por cada camada piezoeléctrica do elemento finito. É utilizada a formulação Lagrangeana actualizada associada à tecnica de Newton - Raphson para a solução iterativa das equações de equilibrio .O modelo pode ser aplicado a cascas piezolaminadas com geometria e carregamento arbitrários. Apresentam-se vários exemplos ilustrativos cujos resultados mostram a eficiencia do modelo proposto.

Otimização de placas compósitas laminadas utilizando redes neuronais e enxames de partículas

Trabalho Final de mestrado para obtenção do grau de Mestre em engenharia Mecância

Criação de programa informático em ambiente Matlab® para estudo de estruturas planas pelo Método de Elementos Finitos.

Os modelos a ser analisados pelo Método de Elementos Finitos são cada vez mais complexos e, nos tempos que correm, seria impensável realizar tais análises sem um apoio computorizado. Existe para esta finalidade uma vasta gama de programas que permitem realizar tarefas que passam pelo desenho de estruturas, análise estática de cargas, análise dinâmica e vibrações, visualização do comportamento físico (deformações) em tempo real, que permitem a otimização da estrutura. Sob o pretexto de permitir a qualquer utilizador uma análise de estruturas simples com o Método dos Elementos Finitos, surge esta tese, onde se irá criar de raiz um programa com interface gráfica no ambiente MATLAB® para análise de estruturas simples com dois tipos de elemento finito, triangular de deformação constante e quadrangular de deformação linear. O software desenvolvido, verificado por comparação com um software comercial dedicado para o efeito, efetua malhagem com elementos bidimensionais triangulares e quadriláteros e resolve modelos arbitrados pelo Método de Elementos Finitos, representando estes resultados visualmente e em formato tabular.

Uma implementação geral de elementos de contorno e sua aplicação a problemas de flexão de placas

o presente trabalho aborda a aplicação do método dos elementos de contorno (MEC) para solução de problemas de flexão linear e geometricamente não-linear de placas semiespessas. Os modelos de placa empregados consideraram a influência do cisalhamento através de teorias de primeira ordem, especificamente as de Mindlin e Reissner. Uma formulação integral unificada dos modelos de placa utilizados é desenvolvida para o operador de Navier do problema, onde foram mantidos alguns termos de ordem superior no tensor deformação de Green. A formulação integral do problema de membrana acoplado ao de flexão é igualment desenvolvida, levando a um sistema de equações integrais não-lineares que descreve completamente problemas de placas que envolvem grandes deslocamentos. Estas equações podem ser particularizadas para problemas de flexão linear e estabilidade elástica. Tendo em vista a necessidade de se considerar derivadas dos deslocamentos translacionais, as equações integrais correspondentes ao gradiente dos deslocamentos também foram deduzidas, caracterizando uma formulação hipersingular. o método empregado para solução numérica do sistema de equações integrais foi o método direto dos elementos de contorno. Um tratamento das integrais fortemente singulares presentes nas equações foi realizado, baseado em expansões assint6ticas dos núcleos. Deste procedimento resulta uma abordagem regularizada que emprega apenas quadraturas padrão de Gauss-Legendre.

Placas maxilo-dentárias em Rincossauros Hyperodapedon Huxley, 1859. Histologia e morfogênese maxilo-dental: nova abordagem.

O estudo das estruturas da placas maxilo-dentárias dos rincossauros Hyperodapedon Huxley, 1859, do Triássico da Formação Santa Maria, em nova abordagem histológica e ontogênica resultou na identificação da natureza do esmalte aprismático verdadeiro, de variados elementos histológicos dentinários e osteológicos, e de um centro de ossificação periosteal primário. Também foram encontradas evidências histológicas dos mecanismos de fusionamento maxilo-dentinário e amelo-maxilar. Com estes elementos, inferimos os modelos de organogênese dental, da ontogênese maxilar e dos mecanismos de fusionamento maxilodental. Encontrou-se uma singular e raríssima coroa dental, imatura e ainda não erupcionada, na região posterior da placa dental e assim evidenciou-se a correta posição da margem odontogenicamente ativa. Adicionalmente inferiu-se a localização da posição da lâmina dentária embrionária. Constatou-se a não formação de alvéolos dentários, de cemento radicular e do espaço necessário à formação do ligamento periodontal e, assim, se deduziu a não formação do folículo dental embriônico. As presenças de especiais elementos anatômicos e histológicos nos tecidos ósseos periapicais evidenciam o crescimento radicular contínuo, enquanto a forma e o fusionamento radicular imediato depõe a favor de uma função dentária fisiológica diferenciada para as baterias dentárias maxilares dos Rincossauros do gênero Hyperodapedon. Os mecanismos que possibilitaram o controle embriônico para a deposição das lamelas de tecido ósseo coronal e seu preciso fusionamento sobre o esmalte dentário, declinam por modificações nas funções tardias do órgão reduzido do esmalte e pela presença de uma membrana oral com funções osteogênicas e também protetivas, situada nas porções posteriores da placa maxilo-dentária em desenvolvimento. Mudanças heterocrônicas no tempo de diferenciação das células da crista neural embriônica e em seus derivados, como a lâmina dentária e órgãos dentários embrionários ou correlacionadas com a organogênese das placas maxilo-dentárias e seus anexos periodontais, todos como condições plesiomórficas para Diápsidas Triássicos, poderiam ser as causas responsáveis pela origem e evolução deste estranho aparelho estomatognático nos clados de Hyperodapedon sp..

Proteção passiva de estruturas e sistemas elétricos contra incêndio aplicada a uma unidade de refino de petróleo

Nesta dissertação visa-se estudar e propor alternativas de solução para a proteção de estruturas e sistemas elétricos contra fogo numa unidade de craqueamento catalítico de uma refinaria de petróleo, por meio de proteção passiva. A proteção passiva tem por finalidade garantir a integridade das estruturas sujeitas a incêndio, durante um determinado período de tempo, para possibilitar, no caso da refinaria, a realização de procedimentos de parada da unidade de forma segura e controlar o incêndio a fim de diminuir a possibilidade de propagação do fogo para outras áreas. Com base em técnicas de análise de riscos fez-se a identificação de zonas potencialmente sujeitas a cenários de acidente envolvendo jato de fogo e/ou incêndio em poça. A delimitação das áreas onde haveria necessidade de proteção passiva foi realizada com base em modelos para jatos de fogo e incêndio em poça já estabelecidos na literatura. O dimensionamento da proteção passiva de estruturas e sistemas elétricos com o uso de diversos materiais usados comercialmente para este fim foi estimado com base em equações empíricas desenvolvidas por Jeanes, 1980, Stanzak, 1973 e PABCO, 1984, e, para alguns casos particulares foi feita uma verificação por solução numérica da equação da condução do calor em meio sólido.Assim, foram determinados quais os materiais mais adequados em cada caso de aplicação e qual a espessura em que deve ser aplicado para que a temperatura no elemento estrutural ou no sistema elétrico não atinja a sua determinada temperatura crítica em um período de tempo pré-determinado. Para os casos de elementos estruturais como colunas de sustentação da unidade de seção cilíndrica, o principal material para proteção passiva é a argamassa projetada e para perfil I, é o emprego de placas de gesso. Já para o caso de sistemas elétricos, podem ser utilizadas tanto tintas intumescentes quanto as mantas reforçadas com fibras minerais, esta escolha depende da geometria do sistema em que será empregado. Da comparação entre estes dois métodos pode-se concluir que o dimensionamento da proteção passiva fazendo o uso das correlações empíricas é menos conservativo que para o caso do uso da equação da difusão do calor resolvida por método numérico. Porém, os resultados diferem dentro de um limite considerado aceitável (em torno de 15%) levando-se em consideração os erros embutidos em cada método de cálculo. É importante mencionar que as correlações empíricas são de mais simples aplicação por possuir apenas operações matemáticas básicas. Usando as correlações empíricas para os perfis cilíndricos de aço (diâmetro de 0,1524 m e espessura de parede de 0,0254 m), a espessura de revestimento estimada com o uso das correlações empíricas necessária para garantir que a temperatura na interface entre os dois materiais não atinja 550°C em duas horas seria de 13,5 mm para argamassa projetada, 19,7 mm para vermiculita com silicato de sódio e 34,5 mm para recobrimento com concreto com proteção do tipo contorno. Fazendo o mesmo cálculo pelo método numérico proposto, os resultados foram de 15,53 mm para argamassa projetada, 22,06 mm para vermiculita com silicato de sódio e 38,98 mm para recobrimento com concreto com proteção do tipo contorno. Fazendo o mesmo cálculo pelo método numérico proposto, os resultados foram de 15,53 mm para argamassa projetada, 22,06 mm para vermiculita com silicato de sódio e 38,98 mm para recobrimento com concreto com proteção do tipo contorno. Cabe ressaltar que com a realização desta dissertação busca-se uma integração entre o mestrado acadêmico e o meio empresarial com o desenvolvimento de trabalhos de natureza acadêmica que tenham aplicação direta na prática. Espera-se assim permitir que Universidade dê retorno à sociedade que a mantém e propiciar que setores da sociedade possam usufruir da capacidade disponível na academia.

Análise dinâmica de placas delgadas utilizando elementos finitos triangulares e retangulares

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