Desenvolvimento de programa em linguagem Matlab® para cálculo pelo método de elementos finitos

Para o projeto de qualquer estrutura existente (edifícios, pontes, veículos, máquinas, etc.) é necessário conhecer as condições de carga, geometria e comportamento de todas as suas partes, assim como respeitar as normativas em vigor nos países nos quais a estrutura será aplicada. A primeira parte de qualquer projeto nesta área passa pela fase da análise estrutural, onde são calculadas todas as interações e efeitos de cargas sobre as estruturas físicas e os seus componentes de maneira a verificar a aptidão da estrutura para o seu uso. Inicialmente parte-se de uma estrutura de geometria simplificada, pondo de parte os elementos físicos irrelevantes (elementos de fixação, revestimentos, etc.) de maneira a simplificar o cálculo de estruturas complexas e, em função dos resultados obtidos da análise estrutural, melhorar a estrutura se necessário. A análise por elementos finitos é a ferramenta principal durante esta primeira fase do projeto. E atualmente, devido às exigências do mercado, é imprescindível o suporte computorizado de maneira a agilizar esta fase do projeto. 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 em análise. Porém, estes programas demostram uma certa complexidade durante a introdução dos parâmetros, levando muitas vezes a resultados errados. Assim sendo, é essencial para o projetista ter uma ferramenta fiável e simples de usar que possa ser usada para fins de projeto de estruturas e otimização. Sobre esta base nasce este projeto tese onde se elaborou um programa com interface gráfica no ambiente Matlab® para a análise de estruturas por elementos finitos, com elementos do tipo Barra e Viga, quer em 2D ou 3D. Este programa permite definir a estrutura por meio de coordenadas, introdução de forma rápida e clara, propriedades mecânicas dos elementos, condições fronteira e cargas a aplicar. Como resultados devolve ao utilizador as reações, deformações e distribuição de tensões nos elementos quer em forma tabular quer em representação gráfica sobre a estrutura em análise. Existe ainda a possibilidade de importação de dados e exportação dos resultados em ficheiros XLS e XLSX, de maneira a facilitar a gestão de informação. Foram realizados diferentes testes e análises de estruturas de forma a validar os resultados do programa e a sua integridade. Os resultados foram todos satisfatórios e convergem para os resultados de outros programas, publicados em livros, e para cálculo a mão feitos pelo autor.

Avaliação do desempenho de ligações tubulares em estruturas de veículos pesados de passageiros por simulação numérica

Nesta Dissertação ir-se-á avaliar o desempenho de algumas ligações existentes nos veículos pesados de passageiros, através do Regulamento Eurocódigo 3. Nos últimos anos ocorreram diversos acidentes envolvendo este tipos de veículos, em que os mesmos causaram vítimas mortais e feridos graves. Serão testadas por simulação numérica algumas ligações pertencentes a elementos constituintes da superestrutura, em que esta é normalmente afectada com a ocorrência de acidentes. Assim sendo, o estudo de nós de ligação tem uma importância fulcral para que uma superestrutura suporte situações extremas e que resista a solicitações externas aplicadas. Iniciou-se esta Dissertação com o estudo da sinistralidade e de acidentes que envolvem veículos pesados de passageiros. No capítulo 2 abordou-se um programa que promove simulações numéricas de acidentes e estudo do comportamento dos passageiros em caso de acidente, sendo referido o Regulamento que homologa os veículos pesados de passageiros e os seus principais métodos. Abordaram-se os principais constituintes da estrutura de um veículo pesado de passageiros. No capítulo 3, é referido o Eurocódigo 3 em termos do estudo das ligações tubulares usadas neste tipo de veículos. No capítulo 4, fez-se o estudo e selecção de elementos a utilizar para a simulação numérica de casos preconizados pelo Eurocódigo 3 e estudaram-se três tipos de ligações que são usadas na construção da superestrutura deste tipo de veículos, tendo-se retirado conclusões deste estudo.

Aplicação de métodos numéricos avançados para a previsão de resistência de ligações adesivas

Com a necessidade de encontrar uma forma de ligar componentes de forma mais vantajosa, surgiram as ligações adesivas. Nos últimos anos, a utilização de juntas adesivas em aplicações industriais tem vindo a aumentar, substituindo alguns métodos de ligação tradicionais, por apresentarem vantagens tais como, redução de concentração de tensões, reduzido peso e facilidade de processamento/fabrico. O seu estudo permite prever a sua resistência e durabilidade. Este trabalho refere-se ao estudo de juntas de sobreposição simples (JSS), nas quais são aplicados os adesivos comerciais que variam desde frágeis e rígidos, como o caso do Araldite® AV138, até adesivos mais dúcteis, como o Araldite® 2015 e o Sikaforce® 7888. Estes são aplicados em substratos de alumínio (AL6082-T651) em juntas com diferentes geometrias e diferentes comprimentos de sobreposição (L), sendo sujeitos a esforços de tracção. Foi feita uma análise dos valores experimentais fornecidos e uma posterior comparação destes com diferentes métodos numéricos baseados em Elementos Finitos (EF). A comparação foi feita por uma análise de Modelos de Dano Coesivo (MDC) e segundo os critérios baseados em tensões e deformações do Método de Elementos Finitos Extendido (MEFE). A utilização destes métodos numéricos capazes de simular o comportamento das juntas poderá levar a uma poupança de recursos e de tempo. A análise por MDC revelou que este método é bastante preciso, excepto para os adesivos que sejam bastante dúcteis. A aplicação de uma outra lei coesiva pode solucionar esse problema. Por sua vez a análise por MEFE demonstrou que esta técnica não é particularmente adequada para o crescimento de dano em modo misto e, comparativamente com o MDC, a sua precisão é bastante inferior.

Adhesive selection for single lap bonded joints: Experimentation and advanced techniques for strength prediction

The integrity of multi-component structures is usually determined by their unions. Adhesive-bonding is often used over traditional methods because of the reduction of stress concentrations, reduced weight penalty, and easy manufacturing. Commercial adhesives range from strong and brittle (e.g., Araldite® AV138) to less strong and ductile (e.g., Araldite® 2015). A new family of polyurethane adhesives combines high strength and ductility (e.g., Sikaforce® 7888). In this work, the performance of the three above-mentioned adhesives was tested in single lap joints with varying values of overlap length (LO). The experimental work carried out is accompanied by a detailed numerical analysis by finite elements, either based on cohesive zone models (CZM) or the extended finite element method (XFEM). This procedure enabled detailing the performance of these predictive techniques applied to bonded joints. Moreover, it was possible to evaluate which family of adhesives is more suited for each joint geometry. CZM revealed to be highly accurate, except for largely ductile adhesives, although this could be circumvented with a different cohesive law. XFEM is not the most suited technique for mixed-mode damage growth, but a rough prediction was achieved.

Projeto da estrutura e sistema de movimentação robotizado de linha automática de corte de chapa

Os robots de movimentação de chapa são bastantes úteis para as empresas de metalomecânica. De facto, cada vez mais existem máquinas de corte por jato de água, laser ou outros processos, nos quais os robots apresentam um papel importante na carga e descarga do material. O trabalho realizado apresenta novas soluções aos sistemas de movimentação existentes no mercado, e permite reduzir os custos na movimentação do material. Este projeto serve essencialmente para chapas em trajetória retilínea, e efetuar o seu levantamento do equipamento e deposição em estrutura de suporte (ou viceversa). A vantagem a ter em conta é a diminuição dos custos de movimentação do material. Neste trabalho apresentou-se a metodologia de dimensionamento de um robot automatizado que transporta chapas com um peso máximo de 3500 kg, tendo por base as normas do EC3-P1.8 e o Método de Elementos Finitos (MEF). No decurso do projeto foram abordadas os seguintes temas: Abordagem inicial da geometria através do Solidworks; Dimensionamento da estrutura por software de Elementos Finitos, o Solidworks; Dimensionamento das correntes, carretos/discos ou coroas e rolamentos; Dimensionamento e seleção dos moto-redutores, bomba de vácuo e ventosas; Cálculo das solicitações em cada membro da estrutura por software de análise estrutural, o Multiframe3D, e respetivo dimensionamento das ligações aparafusadas e soldadas; Elaboração dos desenhos de projeto finais, processos de fabrico e custos; Dimensionamento do acionamento, MG e disposição dos dispositivos no quadro elétrico. Como conclusão refere-se que se conseguiu realizar o projeto e obter uma solução final otimizada, com a ajuda de ferramentas importantes, como sejam o MEF, resultando num equipamento cujas solicitações para a estrutura e sistema de movimentação foram otimizadas, resultando num equipamento eficiente, robusto, seguro e de custo reduzido.

Projecto de equipamento de rotação de moldes de grande porte com recurso ao método de elementos finitos

A fabricação de moldes é uma actividade que se desenvolve em Portugal há várias décadas, com qualidade reconhecida mundialmente. A indústria automóvel é um dos principais destinos de grande parte dos moldes produzidos em Portugal. Algumas das peças produzidas nesses moldes possuem dimensões apreciáveis, tendo em conta o peso e dimensão média das peças injectadas em plástico. São exemplo disso os pára-choques dos automóveis e o tabliê. A produção dos moldes para estas peças é complexa, implicando um longo tempo de maquinagem e a manipulação dos moldes em diferentes fases, com vista ao acesso a todas as faces do molde. A IGM – Indústria Global de Moldes, SA. é a empresa responsável pela produção dos moldes para peças de média dimensão dentro do Grupo SIMOLDES. Atendendo à necessidade permanente de rodar a 90o moldes de grande porte, que podem apresentar pesos na ordem das 10 a 30 toneladas, e não existindo no mercado qualquer solução adaptável a esta necessidade, a empresa entendeu por bem levar a efeito o projecto desse equipamento, atendendo ao compromisso custo-benefício que torne viável a realização prática do mesmo. Após os esboços iniciais e uma discussão interactiva e iterativa com a empresa, foram analisadas as diferentes soluções entendidas como viáveis, tendo sido escolhido um dos anteprojectos realizados. Foram ainda discutidas as diversas alternativas de accionamento. Com base nesse anteprojecto, a estrutura foi optimizada e verificada através do Método de Elementos Finitos, tendo sido elaborado o projecto final, com o grau de detalhe necessário à sua fabricação e implementação na empresa.

Estudo e otimização de juntas do tipo Tpeel soldadas, adesivas e híbridas

A necessidade de utilizar métodos de ligação que melhor satisfaçam as necessidades de projeto tem causado a crescente utilização das juntas adesivas, em detrimento dos métodos tradicionais tais como a soldadura, ligações aparafusadas e rebitadas. A sua utilização em diversas aplicações industriais justifica-se pela redução de peso, redução de concentrações de tensões, isolamento acústico e melhor resistência à corrosão. Contudo, também apresentam desvantagens, como a necessidade de preparação das juntas, a fraca resistência a esforços de arrancamento e a complexidade da previsão da sua resistência. As juntas híbridas são obtidas por combinação de uma técnica tradicional com uma ligação adesiva. As juntas híbridas adesivas-soldadas obtêm-se através da combinação da ligação adesiva com a ligação soldada, sendo a soldadura de resistência por pontos a técnica de soldadura mais usada no fabrico deste tipo de juntas. A sinergia entre ligação adesiva e soldadura por pontos oferece vantagens competitivas em relação às ligações adesivas, tais como superior resistência e rigidez, e maior resistência ao arrancamento e à fadiga. No presente trabalho é apresentado um estudo experimental e numérico de juntas T-peel soldadas, adesivas e híbridas (adesivas-soldadas) solicitadas ao arrancamento. Considerouse o adesivo frágil Araldite® AV138 e os adesivos dúcteis Araldite® 2015 e Sikaforce® 7752 e aderentes de aço (C45E). Foi realizada uma análise dos valores experimentais e efetuada uma comparação destes valores com os resultados obtidos pelo Método de Elementos Finitos (MEF) no software ABAQUS®, que incluiu uma análise de tensões na camada de adesivo e previsão do comportamento das juntas por MDC. Observou-se que, dos três adesivos em estudo, o adesivo Sikaforce® 7752 é o que apresenta o melhor desempenho na ligação de juntas T-peel. A boa concordância entre os resultados experimentais e numéricos permitiu validar a utilização de MDC para previsão da resistência de juntas T-peel adesivas e híbridas. Assim, o presente trabalho representa uma base para posterior aplicação no projeto deste tipo de ligação, com as vantagens decorrentes na redução do tempo de projeto e maior facilidade de otimização.

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.

Comparação de técnicas analíticas e numéricas para previsão de resistência de juntas adesivas de sobreposição simples

As juntas adesivas têm vindo a ser usadas em diversas áreas e contam com inúmeras aplicações práticas. Devido ao fácil e rápido fabrico, as juntas de sobreposição simples (JSS) são um tipo de configuração bastante comum. O aumento da resistência, a redução de peso e a resistência à corrosão são algumas das vantagens que este tipo de junta oferece relativamente aos processos de ligação tradicionais. Contudo, a concentração de tensões nas extremidades do comprimento da ligação é uma das principais desvantagens. Existem poucas técnicas de dimensionamento precisas para a diversidade de ligações que podem ser encontradas em situações reais, o que constitui um obstáculo à utilização de juntas adesivas em aplicações estruturais. O presente trabalho visa comparar diferentes métodos analíticos e numéricos na previsão da resistência de JSS com diferentes comprimentos de sobreposição (LO). O objectivo fundamental é avaliar qual o melhor método para prever a resistência das JSS. Foram produzidas juntas adesivas entre substratos de alumínio utilizando um adesivo époxido frágil (Araldite® AV138), um adesivo epóxido moderadamente dúctil (Araldite® 2015), e um adesivo poliuretano dúctil (SikaForce® 7888). Consideraram-se diferentes métodos analíticos e dois métodos numéricos: os Modelos de Dano Coesivo (MDC) e o Método de Elementos Finitos Extendido (MEFE), permitindo a análise comparativa. O estudo possibilitou uma percepção crítica das capacidades de cada método consoante as características do adesivo utilizado. Os métodos analíticos funcionam apenas relativamente bem em condições muito específicas. A análise por MDC com lei triangular revelou ser um método bastante preciso, com excepção de adesivos que sejam bastante dúcteis. Por outro lado, a análise por MEFE demonstrou ser uma técnica pouco adequada, especialmente para o crescimento de dano em modo misto.

Aplicação do modelo escoras e tirantes a vigas parede de acordo com o Eurocódigo 2:2010

O Modelo Escoras e Tirantes surgiu no início do Século XX quando Ritter e Mörsch estabeleceram a analogia entre uma treliça clássica e uma viga de betão armado. Desde então pesquisadores têm estudado esse modelo como método de dimensionamento. A partir dos anos 90 várias normas apresentaram a utilização do modelo escoras e tirantes como relevante no dimensionamento de elementos de betão armado. Os critérios de segurança do Modelo Escoras e Tirantes são neste trabalho explicados de acordo com o Eurocódigo 2:2010, e comparadas com as normas NBR 6118:2014 e ACI 318:2011. É unanime dizer que a utilização do método torna-se mais vantajosa para regiões de descontinuidade. Em todos os elementos de betão armado, o Método Escoras e Tirantes são representações dos campos de tensão idealizados por elementos comprimidos e tracionados. Para a definição destes elementos é proposto o processo de “caminho de carga” em que conhecidas as tensões elásticas e suas direções principais utilizando o Método dos Elementos Finitos o modelo das escoras e tirantes é de fácil concepção. É também possível a definição deste a partir de modelos padrão já concebidos para determinados tipos de elementos estruturais de betão armado. Para o elemento descontinuo viga-parede estudado foram apresentados cinco modelos de cálculo até otimizar a solução validando as tensões com o Método dos Elementos Finitos. Em todos os modelos foram analisadas a definição do Modelo Escoras e Tirantes, a resistência das escoras, dos tirantes e dos nós até chegar à solução construtiva da viga-parede. Concluiu-se que a definição do modelo é a chave do dimensionamento.

Numerical simulation of galvanized rebars pullout

The usage of rebars in construction is the most common method for reinforcing plain concrete and thus bridging the tensile stresses along the concrete crack surfaces. Usually design codes for modelling the bond behaviour of rebars and concrete suggest a local bond stress – slip relationship that comprises distinct reinforcement mechanisms, such as adhesion, friction and mechanical anchorage. In this work, numerical simulations of pullout tests were performed using the finite element method framework. The interaction between rebar and concrete was modelled using cohesive elements. Distinct local bond laws were used and compared with ones proposed by the Model Code 2010. Finally an attempt was made to model the geometry of the rebar ribs in conjunction with a material damaged plasticity model for concrete.

Recycled steel fibre reinforced concrete failing in bending and in shear

Recent research is showing that the addition of Recycled Steel Fibres (RSF) from wasted tyres can decrease significantly the brittle behaviour of cement based materials, by improving its toughness and post-cracking resistance. In this sense, Recycled Steel Fibre Reinforced Concrete (RSFRC) seems to have the potential to constitute a sustainable material for structural and non-structural applications. To assess this potential, experimental and numerical research was performed on the use of RSFRC in elements failing in bending and in beams failing in shear. The values of the fracture mode I parameters of the developed RSFRC were determined by performing inverse analysis with test results obtained in three point notched beam bending tests. To assess the possibility of using RSF as shear reinforcement in Reinforced Concrete (RC) beams, three point bending tests were executed with three series of RSFRC beams flexurally reinforced with a relatively high reinforcement ratio of longitudinal steel bars in order to assure shear failure for all the tested beams. By performing material nonlinear simulations with a computer program based on the finite element method (FEM), the applicability of the fracture mode I crack constitutive law derived from the inverse analysis is assessed for the prediction of the behaviour of these beams. The performance of the formulation proposed by RILEM TC 162 TDF and CEB-FIP 2010 for the prediction of the shear resistance of fibre reinforced concrete elements was also evaluated.

Finite integration methods for isospectral flows

In this paper we consider the approximate computation of isospectral flows based on finite integration methods( FIM) with radial basis functions( RBF) interpolation,a new algorithm is developed. Our method ensures the symmetry of the solutions. Numerical experiments demonstrate that the solutions have higher accuracy by our algorithm than by the second order Runge- Kutta( RK2) method.

Estudo numérico sobre o reforço de pilares de betão armado com cantoneiras e presilhas metálicas: análise paramétrica e modelo de cálculo

Dissertação de mestrado integrado em Engenharia Civil

Adaptive Verfahren zur numerischen Berechnung von Reaktions-Diffusions-Systemen

Reaction-diffusion problems, finite elements, unstructured grid, grid adaption, W-method, stiffness, local partitioning, excitable medium, spiral wave drift