The application of wave induced forces to a two-dimensional finite element long wave hydrodynamic model

This paper describes the modification of a two-dimensional finite element long wave hydrodynamic model in order to predict the net current and water levels attributable to the influences of waves. Tests examine the effects of the application of wave induced forces, including comparisons to a physical experiment. An example of a real river system is presented with comparisons to measured data, which demonstrate the importance of simulating the combined effects of tides and waves upon hydrodynamic behavior. (C) 2002 Elsevier Science Ltd. All rights reserved.

Comparação de desempenho entre a formulação direta do Método dos Elementos de Contorno com Funções Radiais e o Método dos Elementos Finitos em problemas de Poisson e Helmholtz

O presente trabalho objetiva avaliar o desempenho do MECID (Método dos Elementos de Contorno com Interpolação Direta) para resolver o termo integral referente à inércia na Equação de Helmholtz e, deste modo, permitir a modelagem do Problema de Autovalor assim como calcular as frequências naturais, comparando-o com os resultados obtidos pelo MEF (Método dos Elementos Finitos), gerado pela Formulação Clássica de Galerkin. Em primeira instância, serão abordados alguns problemas governados pela equação de Poisson, possibilitando iniciar a comparação de desempenho entre os métodos numéricos aqui abordados. Os problemas resolvidos se aplicam em diferentes e importantes áreas da engenharia, como na transmissão de calor, no eletromagnetismo e em problemas elásticos particulares. Em termos numéricos, sabe-se das dificuldades existentes na aproximação precisa de distribuições mais complexas de cargas, fontes ou sorvedouros no interior do domínio para qualquer técnica de contorno. No entanto, este trabalho mostra que, apesar de tais dificuldades, o desempenho do Método dos Elementos de Contorno é superior, tanto no cálculo da variável básica, quanto na sua derivada. Para tanto, são resolvidos problemas bidimensionais referentes a membranas elásticas, esforços em barras devido ao peso próprio e problemas de determinação de frequências naturais em problemas acústicos em domínios fechados, dentre outros apresentados, utilizando malhas com diferentes graus de refinamento, além de elementos lineares com funções de bases radiais para o MECID e funções base de interpolação polinomial de grau (um) para o MEF. São geradas curvas de desempenho através do cálculo do erro médio percentual para cada malha, demonstrando a convergência e a precisão de cada método. Os resultados também são comparados com as soluções analíticas, quando disponíveis, para cada exemplo resolvido neste trabalho.

Simulação numérica do comportamento dinâmico de uma laje

Nesta dissertação pretende-se simular o comportamento dinâmico de uma laje de betão armado aplicando o Método de Elementos Finitos através da sua implementação no programa FreeFEM++. Este programa permite-nos a análise do modelo matemático tridimensional da Teoria da Elasticidade Linear, englobando a Equação de Equilíbrio, Equação de Compatibilidade e Relações Constitutivas. Tratando-se de um problema dinâmico é necessário recorrer a métodos numéricos de Integração Directa de modo a obter a resposta em termos de deslocamento ao longo do tempo. Para este trabalho escolhemos o Método de Newmark e o Método de Euler para a discretização temporal, um pela sua popularidade e o outro pela sua simplicidade de implementação. Os resultados obtidos pelo FreeFEM++ são validados através da comparação com resultados adquiridos a partir do SAP2000 e de Soluções Teóricas, quando possível.

Study of the one dimensional and transient bioheat transfer equation: Multi-layer solution development and applications

: In this work we derive an analytical solution given by Bessel series to the transient and one-dimensional (1D) bioheat transfer equation in a multi-layer region with spatially dependent heat sources. Each region represents an independent biological tissue characterized by temperature-invariant physiological parameters and a linearly temperature dependent metabolic heat generation. Moreover, 1D Cartesian, cylindrical or spherical coordinates are used to define the geometry and temperature boundary conditions of first, second and third kinds are assumed at the inner and outer surfaces. We present two examples of clinical applications for the developed solution. In the first one, we investigate two different heat source terms to simulate the heating in a tumor and its surrounding tissue, induced during a magnetic fluid hyperthermia technique used for cancer treatment. To obtain an accurate analytical solution, we determine the error associated with the truncated Bessel series that defines the transient solution. In the second application, we explore the potential of this model to study the effect of different environmental conditions in a multi-layered human head model (brain, bone and scalp). The convective heat transfer effect of a large blood vessel located inside the brain is also investigated. The results are further compared with a numerical solution obtained by the Finite Element Method and computed with COMSOL Multi-physics v4.1 (c). (c) 2013 Elsevier Ltd. All rights reserved.

Interlaminar fracture characterization of a carbon-epoxy composite in pure mode II

The interlaminar fracture toughness in pure mode II (GIIc) of a Carbon-Fibre Reinforced Plastic (CFRP) composite is characterized experimentally and numerically in this work, using the End-Notched Flexure (ENF) fracture characterization test. The value of GIIc was extracted by a new data reduction scheme avoiding the crack length measurement, named Compliance-Based Beam Method (CBBM). This method eliminates the crack measurement errors, which can be non-negligible, and reflect on the accuracy of the fracture energy calculations. Moreover, it accounts for the Fracture Process Zone (FPZ) effects. A numerical study using the Finite Element Method (FEM) and a triangular cohesive damage model, implemented within interface finite elements and based on the indirect use of Fracture Mechanics, was performed to evaluate the suitability of the CBBM to obtain GIIc. This was performed comparing the input values of GIIc in the numerical models with the ones resulting from the application of the CBBM to the numerical load-displacement (P-) curve. In this numerical study, the Compliance Calibration Method (CCM) was also used to extract GIIc, for comparison purposes.

Single-lap joints of similar and dissimilar adherends bonded with an acrylic adhesive

In this study, the tensile strength of single-lap joints (SLJs) between similar and dissimilar adherends bonded with an acrylic adhesive was evaluated experimentally and numerically. The adherend materials included polyethylene (PE), polypropylene (PP), carbon-epoxy (CFRP), and glass-polyester (GFRP) composites. The following adherend combinations were tested: PE/PE, PE/PP, PE/CFRP, PE/GFRP, PP/PP, CFRP/CFRP, and GFRP/GFRP. One of the objectives of this work was to assess the influence of the adherends stiffness on the strength of the joints since it significantly affects the peel stresses magnitude in the adhesive layer. The experimental results were also used to validate a new mixed-mode cohesive damage model developed to simulate the adhesive layer. Thus, the experimental results were compared with numerical simulations performed in ABAQUS®, including a developed mixed-mode (I+II) cohesive damage model, based on the indirect use of fracture mechanics and implemented within interface finite elements. The cohesive laws present a trapezoidal shape with an increasing stress plateau, to reproduce the behaviour of the ductile adhesive used. A good agreement was found between the experimental and numerical results.

Ligação Adesiva de Polímeros Particularmente Difíceis de Colar

As ligações adesivas têm sido cada vez mais utilizadas nos últimos anos em detrimento de outros métodos tais como a soldadura, ligações aparafusadas e ligações rebitadas. Os plásticos de Engenharia têm um papel cada vez mais preponderante na indústria, devido às suas excelentes propriedades. Neste trabalho foram considerados três polímeros diferentes, o Policloreto de Vinilo (PVC) e o Polipropileno (PP) dado o seu baixo custo e peso e a superfície quimicamente inerte e o Politetrafluoretileno (PTFE) devido às suas boas propriedades químicas e excelentes propriedades de deslizamento. No entanto, estes materiais possuem uma baixa energia de superfície e, por isso, são muito difíceis de colar com mais relevância para o PTFE. Assim, após um estudo preliminar foi escolhido, para realizar as colagens necessárias, um adesivo da Tamarron Technology “Tam Tech Adhesive”, próprio para este tipo de substratos difíceis de colar. Posteriormente foi efetuada a sua caraterização através de ensaios de provetes maciços à tração. O principal objetivo deste trabalho foi estudar juntas de sobreposição simples de materiais poliméricos difíceis de colar tais como o PTFE, PP e PVC com recurso a um adesivo que não necessitasse de preparação de superfície. Foram fabricadas juntas de sobreposição simples (JSS) segundo os métodos Lap Shear (LS) e Block Shear (BS) dos três materiais referidos anteriormente e realizados os respetivos ensaios para avaliar o comportamento mecânico das ligações adesivas. Os materiais utilizados como substratos foram também submetidos a ensaios de tração com a finalidade de obter o módulo de elasticidade e as suas propriedades de resistência. Os substratos envolvidos nas juntas adesivas não sofreram qualquer preparação especial das superfícies. Na maioria dos casos consistiu apenas numa limpeza das superfícies com álcool etílico. Contudo, para o PTFE também se experimentou a preparação por abrasão com lixa e por chama. Foi também efetuado um trabalho de simulação numérica por elementos finitos utilizando um modelo de dano coesivo triangular. As resistências ao corte obtidas são superiores em BS comparativamente a LS, exceção feita aos substratos de PTFE aonde os resultados são similares. O tratamento por chama melhorou a resistência mecânica das juntas. Verificou-se também que o modelo numérico simulou adequadamente o comportamento das juntas principalmente das LS.

Projeto de um centro de roda para veículo de competição

A presente tese surgiu da cooperação entre a empresa Semog Racing e o Instituto Superior de Engenharia do Porto, no âmbito da unidade curricular de DPEST, englobado no 2º ano do Mestrado de Engenharia Mecânica, ramo de Construções Mecânicas. Esta dissertação tem como objetivo principal, o projeto de um novo centro de roda para substituir o existente num veículo de competição da Semog Racing. Este centro deverá ser dimensionado para suportar o peso do veículo e os esforços transmitidos à roda. Pretende-se igualmente que o novo modelo apresente um design apelativo e um baixo custo de produção para poder vir ser a comercializado. Na primeira fase, será estudado o centro de roda original fornecido pela Semog Racing. Este estudo engloba a modelação do componente através do software SolidWorks® e uma fase de simulações para diferentes condições de carregamento. O estudo será complementado com a realização de ensaios experimentais para validação do modelo numérico. A segunda fase é dedicada ao desenvolvimento do novo modelo de centro de roda, focando as características mecânicas e o design. Este tem como base o estudo numérico realizado para a roda original, tendo sempre como objetivo final garantir que o novo centro de roda cumpra todos os requisitos. O caminho seguido no processo de otimização é suportado em simulações numéricas pelo método dos elementos finitos, o qual permite aferir quais os pontos críticos a corrigir. No final, será apresentado um novo modelo de centro de roda capaz de suportar as cargas previstas de serviço, que apresente um baixo custo de fabrico e design apelativo.

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.

Avaliação numérica e experimental da carga admissível de vigas compostas para monocarris de sistemas de elevação

Este trabalho foi realizado com o apoio, e em colaboração da empresa Flexcrane, que constrói e comercializa soluções flexíveis de pontes rolantes para cargas ligeiras. Com cargas máximas de 2 toneladas estas pontes e seus caminhos de rolamento são construídas essencialmente enformados a frio. O objetivo deste trabalho é a verificação numérica e experimental das cargas que este tipo de estrutura suporta, nas suas diferentes configurações. Numa aproximação por engenharia inversa, com a realização deste trabalho, foi possível determinar a carga admissível das vigas compostas utilizadas, assim como o seu comportamento quando em carga. Para obtenção destes resultados foram realizadas simulações computacionais recorrendo ao método de elementos finitos, utilizando para isso o software Ansys Workbench e ainda ensaios experimentais que permitiram validar os resultados obtidos numericamente e consolidar o conhecimento sobre o comportamento desta solução comercial. Por ser de construção modular, são possíveis diferentes combinações para o produto final, por esta razão foram estudados numericamente diferentes configurações possíveis, com e sem reforços. Foi realizada a validação experimental para a configuração mais usual, fazendo recurso a um dispositivo para ensaio de estruturas e componentes, pertencente ao Isep. Foi ainda utilizada a técnica extensométrica para avaliar as tensões em pontos selecionados e compara-los com os obtidos numericamente. Concluiu-se pelo ótimo desempenho dos produtos Flexcrane.

Shear strengthening of damaged reinforced concrete beams with hybrid composite plates

This paper aims to evaluate experimentally the potentialities of Hybrid Composite Plates (HCPs) technique for the shear strengthening of reinforced concrete (RC) beams that were previously subjected to intense damage in shear. HCP is a thin plate of Strain Hardening Cementitious Composite (SHCC) reinforced with Carbon Fiber Reinforced Polymer (CFRP) laminates. For this purpose, an experimental program composed of two series of beams (rectangular and T cross section) was executed to assess the strengthening efficiency of this technique. In the first step of this experimental program, the control beams, without steel stirrups, were loaded up to their shear failure, and fully unloaded. Then, these pre-damaged beams were shear strengthened by applying HCPs to their lateral faces by using a combination of epoxy adhesive and mechanical anchors. The bolts were applied with a certain torque in order to increase the concrete confinement. The obtained results showed that the increase of load carrying capacity of the damaged strengthened beams when HCPs were applied with epoxy adhesive and mechanical anchors was 2 and 2.5 times of the load carrying capacity of the corresponding reference beams (without HCPs) for the rectangular and T cross section beam series, respectively. To further explore the potentialities of the HCPs technique for the shear strengthening, the experimental tests were simulated using an advanced numerical model by a FEM-based computer program. After demonstration the good predictive performance of the numerical model, a parametric study was executed to highlight the influence of SHCC as an alternative for mortar, as well as the influence of torque level applied to the mechanical anchors, on the load carrying capacity of beams strengthened with the proposed technique.

A sixth-order finite volume method for the 1D biharmonic operator: application to intramedullary nail simulation

A new very high-order finite volume method to solve problems with harmonic and biharmonic operators for one- dimensional geometries is proposed. The main ingredient is polynomial reconstruction based on local interpolations of mean values providing accurate approximations of the solution up to the sixth-order accuracy. First developed with the harmonic operator, an extension for the biharmonic operator is obtained, which allows designing a very high-order finite volume scheme where the solution is obtained by solving a matrix-free problem. An application in elasticity coupling the two operators is presented. We consider a beam subject to a combination of tensile and bending loads, where the main goal is the stress critical point determination for an intramedullary nail.

Mortar-Techniken zur Behandlung von Grenzschichtproblemen

Elliptic differential equations, finite element method, mortar element method, streamline diffusion FEM, upwind method, numerical method, error estimate, interpolation operator, grid generation, adaptive refinement

Stabilized finite element methods applied to transient convection-diffusion-reaction and population balance equations

Magdeburg, Univ., Fak. für Mathematik, Diss., 2011

Development of a new isogeometric finite element and its application for Lamb wave based structural health monitoring

Magdeburg, Univ., Fak. für Maschinenbau, Diss., 2012