Elementos de contorno adaptables

Como es bien sabido, en el método de los elementos finitos se suele hablar de dos tipos de convergencia. La primera, o convergencia h, se refiere a la mejora del resultado que se obtiene refinando la malla. Debido a la correspondencia elemento-variables nodales-funciones de interpolación, ello implica un ajuste progresivo de los resultados en aquellas zonas donde se produce el refinamiento. Se trata del método más usado cuando, de forma pragmática, se desea tener una idea de la convergencia de los resultados. Su principal inconveniente radica en el hecho que cada refinamiento exige el cálculo de matrices de rigidez diferentes de las anteriores, de modo que la información debe ser rehecha en cada caso y, por tanto, los costes son elevados. El segundo método analiza la convergencia p, o refinamiento de la aproximación mediante el incremento del grado del polinomio definido sobre cada elemento. Se trata de abandonar la idea de asociar a cada nodo el valor físico de la variable correspondiente en la aproximación típica: u ~ a1Ø1 + a2Ø2 + a3Ø3+ … + anØn; donde las funciones Ø son unidad en el nodo correspondiente y cero en el resto. Por el contrario, se vuelve a la idea original de Ritz, semejante al de un desarrollo en la serie de Fourier, donde las funciones Ø están definidas globalmente y los coeficientes de ponderación no tienen por qué presentar un significado físico concreto. Evidentemente la vuelta no es total; se siguen manteniendo elementos y dentro de cada uno de ellos se establece una jerarquía de funciones Øi. Con esta situación intermedia entre la globalidad absoluta de Ritz y la correspondencia absoluta de la discretización con las variables se consigue, por un lado, mantener una versatilidad suficiente para el ajuste por trozos y, por otro, refinar la aproximación de forma inteligente ya que, al igual que sucede en una serie de Fourier, cada término que se añade produce un efecto menor, lo que posibilita el truncamiento cuando se alcanza un determinado nivel de precisión. Además, puesto que cada Ø tiene un soporte perfectamente definido desde un principio, cada etapa del refinamiento aprovecha todos los cálculos anteriores y sólo se necesita evaluar los nuevos términos de la matriz de rigidez. La primera idea fue propuesta por Zienckiewicz et al.(1970), y posteriormente han desarrollado el método Szabo et al.(1978), Babuska (1975,1978), Peano (1978)etc. El proceso operativo incluye así: a)Establecimiento de una malla amplia sobre el dominio a analizar; b)Definición de una jerarquía de funciones de interpolación dentro de cada elemento; c)Establecimiento de un "indicador" de las zonas que precisen la adición de nuevas funciones jerarquizadas; d)Establecimiento de un "estimador a posteriori" que evalúe el error cometido y precise el momento en que pueda ser detenido el proceso. Un método que sigue los pasos anteriores se denomina autoadaptable y, como se puede comprender, resulta interesantísimo para problemas no triviales. En este artículo, se contempla la posibilidad de extender las ideas anteriores al método de los elementos de contorno.

Elementos de contorno adaptables

Como es bien sabido, en el método de los elementos finitos se suele hablar de dos tipos de convergencia. La primera, o convergencia h, se refiere a la mejora del resultado que se obtiene refinando la malla. Debido a la correspondencia elemento-variables nodales-funciones de interpolación, ello implica un ajuste progresivo de los resultados en aquellas zonas donde se produce el refinamiento. Se trata del método más usado cuando, de forma pragmática, se desea tener una idea de la convergencia de los resultados. Su principal inconveniente radica en el hecho que cada refinamiento exige el cálculo de matrices de rigidez diferentes de las anteriores, de modo que la información debe ser rehecha en cada caso y, por tanto, los costes son elevados. El segundo método analiza la convergencia p, o refinamiento de la aproximación mediante el incremento del grado del polinomio definido sobre cada elemento. Se trata de abandonar la idea de asociar a cada nodo el valor físico de la variable correspondiente en la aproximación típica: u ~ a1Ø1 + a2Ø2 + a3Ø3+ … + anØn; donde las funciones Ø son unidad en el nodo correspondiente y cero en el resto. Por el contrario, se vuelve a la idea original de Ritz, semejante al de un desarrollo en la serie de Fourier, donde las funciones Ø están definidas globalmente y los coeficientes de ponderación no tienen por qué presentar un significado físico concreto. Evidentemente la vuelta no es total; se siguen manteniendo elementos y dentro de cada uno de ellos se establece una jerarquía de funciones Øi. Con esta situación intermedia entre la globalidad absoluta de Ritz y la correspondencia absoluta de la discretización con las variables se consigue, por un lado, mantener una versatilidad suficiente para el ajuste por trozos y, por otro, refinar la aproximación de forma inteligente ya que, al igual que sucede en una serie de Fourier, cada término que se añade produce un efecto menor, lo que posibilita el truncamiento cuando se alcanza un determinado nivel de precisión. Además, puesto que cada Ø tiene un soporte perfectamente definido desde un principio, cada etapa del refinamiento aprovecha todos los cálculos anteriores y sólo se necesita evaluar los nuevos términos de la matriz de rigidez. La primera idea fue propuesta por Zienckiewicz et al.(1970), y posteriormente han desarrollado el método Szabo et al.(1978), Babuska (1975,1978), Peano (1978)etc. El proceso operativo incluye así: a)Establecimiento de una malla amplia sobre el dominio a analizar; b)Definición de una jerarquía de funciones de interpolación dentro de cada elemento; c)Establecimiento de un "indicador" de las zonas que precisen la adición de nuevas funciones jerarquizadas; d)Establecimiento de un "estimador a posteriori" que evalúe el error cometido y precise el momento en que pueda ser detenido el proceso. Un método que sigue los pasos anteriores se denomina autoadaptable y, como se puede comprender, resulta interesantísimo para problemas no triviales. En este artículo, se contempla la posibilidad de extender las ideas anteriores al método de los elementos de contorno.

Método multigrid algébrico: reutilização das estruturas multigrid no transporte de contaminantes

A necessidade de obter solução de grandes sistemas lineares resultantes de processos de discretização de equações diferenciais parciais provenientes da modelagem de diferentes fenômenos físicos conduz à busca de técnicas numéricas escaláveis. Métodos multigrid são classificados como algoritmos escaláveis.Um estimador de erros deve estar associado à solução numérica do problema discreto de modo a propiciar a adequada avaliação da solução obtida pelo processo de aproximação. Nesse contexto, a presente tese caracteriza-se pela proposta de reutilização das estruturas matriciais hierárquicas de operadores de transferência e restrição dos métodos multigrid algébricos para acelerar o tempo de solução dos sistemas lineares associados à equação do transporte de contaminantes em meio poroso saturado. Adicionalmente, caracteriza-se pela implementação das estimativas residuais para os problemas que envolvem dados constantes ou não constantes, os regimes de pequena ou grande advecção e pela proposta de utilização das estimativas residuais associadas ao termo de fonte e à condição inicial para construir procedimentos adaptativos para os dados do problema. O desenvolvimento dos códigos do método de elementos finitos, do estimador residual e dos procedimentos adaptativos foram baseados no projeto FEniCS, utilizando a linguagem de programação PYTHONR e desenvolvidos na plataforma Eclipse. A implementação dos métodos multigrid algébricos com reutilização considera a biblioteca PyAMG. Baseado na reutilização das estruturas hierárquicas, os métodos multigrid com reutilização com parâmetro fixo e automática são propostos, e esses conceitos são estendidos para os métodos iterativos não-estacionários tais como GMRES e BICGSTAB. Os resultados numéricos mostraram que o estimador residual captura o comportamento do erro real da solução numérica, e fornece algoritmos adaptativos para os dados cuja malha retornada produz uma solução numérica similar à uma malha uniforme com mais elementos. Adicionalmente, os métodos com reutilização são mais rápidos que os métodos que não empregam o processo de reutilização de estruturas. Além disso, a eficiência dos métodos com reutilização também pode ser observada na solução do problema auxiliar, o qual é necessário para obtenção das estimativas residuais para o regime de grande advecção. Esses resultados englobam tanto os métodos multigrid algébricos do tipo SA quanto os métodos pré-condicionados por métodos multigrid algébrico SA, e envolvem o transporte de contaminantes em regime de pequena e grande advecção, malhas estruturadas e não estruturadas, problemas bidimensionais, problemas tridimensionais e domínios com diferentes escalas.

Simulação e modelagem computacional de dados de espalhamento à baixos ângulos - enfoque em estruturas de alta simetria

Esta tese apresenta uma abordagem para a criação rápida de modelos em diferentes geometrias (complexas ou de alta simetria) com objetivo de calcular a correspondente intensidade espalhada, podendo esta ser utilizada na descrição de experimentos de es- palhamento à baixos ângulos. A modelagem pode ser realizada com mais de 100 geome- trias catalogadas em um Banco de Dados, além da possibilidade de construir estruturas a partir de posições aleatórias distribuídas na superfície de uma esfera. Em todos os casos os modelos são gerados por meio do método de elementos finitos compondo uma única geometria, ou ainda, compondo diferentes geometrias, combinadas entre si a partir de um número baixo de parâmetros. Para realizar essa tarefa foi desenvolvido um programa em Fortran, chamado de Polygen, que permite modelar geometrias convexas em diferentes formas, como sólidos, cascas, ou ainda com esferas ou estruturas do tipo DNA nas arestas, além de usar esses modelos para simular a curva de intensidade espalhada para sistemas orientados e aleatoriamente orientados. A curva de intensidade de espalhamento é calculada por meio da equação de Debye e os parâmetros que compõe cada um dos modelos, podem ser otimizados pelo ajuste contra dados experimentais, por meio de métodos de minimização baseados em simulated annealing, Levenberg-Marquardt e algorítmicos genéticos. A minimização permite ajustar os parâmetros do modelo (ou composição de modelos) como tamanho, densidade eletrônica, raio das subunidades, entre outros, contribuindo para fornecer uma nova ferramenta para modelagem e análise de dados de espalhamento. Em outra etapa desta tese, é apresentado o design de modelos atomísticos e a sua respectiva simulação por Dinâmica Molecular. A geometria de dois sistemas auto-organizado de DNA na forma de octaedro truncado, um com linkers de 7 Adeninas e outro com linkers de ATATATA, foram escolhidas para realizar a modelagem atomística e a simulação por Dinâmica Molecular. Para este sistema são apresentados os resultados de Root Mean Square Deviations (RMSD), Root Mean Square Fluctuations (RMSF), raio de giro, torção das hélices duplas de DNA além da avaliação das ligações de Hidrogênio, todos obtidos por meio da análise de uma trajetória de 50 ns.

Um modelo multiescala concorrente para representar o processo de fissuração do concreto.

Este trabalho propõe uma técnica de modelagem multiescala concorrente do concreto considerando duas escalas distintas: a mesoescala, onde o concreto é modelado como um material heterogêneo, e a macroescala, na qual o concreto é tratado como um material homogêneo. A heterogeneidade da estrutura mesoscópica do concreto é idealizada considerando três fases distintas, compostas pelos agregados graúdos e argamassa (matriz), estes considerados materiais homogêneos, e zona de transição interfacial (ZTI), tratada como a parte mais fraca entre as três fases. O agregado graúdo é gerado a partir de uma curva granulométrica e posicionado na matriz de forma aleatória. Seu comportamento mecânico é descrito por um modelo constitutivo elástico-linear, devido a sua maior resistência quando comparado com as outras duas fases do concreto. Elementos finitos contínuos com alta relação de aspecto em conjunto com um modelo constitutivo de dano são usados para representar o comportamento não linear do concreto, decorrente da iniciação de fissuras na ZTI e posterior propagação para a matriz, dando lugar à formação de macrofissuras. Os elementos finitos de interface com alta relação de aspecto são inseridos entre todos os elementos regulares da matriz e entre os da matriz e agregados, representando a ZTI, tornando-se potenciais caminhos de propagação de fissuras. No estado limite, quando a espessura do elemento de interface tende a zero (h ?0) e, consequentemente, a relação de aspecto tende a infinito, estes elementos apresentam a mesma cinemática da aproximação contínua de descontinuidades fortes (ACDF), sendo apropriados para representar a formação de descontinuidades associados a fissuras, similar aos modelos coesivos. Um modelo de dano à tração é proposto para representar o comportamento mecânico não linear das interfaces, associado à formação de fissuras, ou até mesmo ao eventual fechamento destas. A fim de contornar os problemas causados pela malha de elementos finitos de transição entre as malhas da macro e da mesoescala, que, em geral, apresentam diferenças expressivas 5 de refinamento, utiliza-se uma técnica recente de acoplamento de malhas não conformes. Esta técnica é baseada na definição de elementos finitos de acoplamento (EFAs), os quais são capazes de estabelecer a continuidade de deslocamento entre malhas geradas de forma completamente independentes, sem aumentar a quantidade total de graus de liberdade do problema, podendo ser utilizados tanto para acoplar malhas não sobrepostas quanto sobrepostas. Para tornar possível a análise em multiescala em casos nos quais a região de localização de deformações não pode ser definida a priori, propõe-se uma técnica multiescala adaptativa. Nesta abordagem, usa-se a distribuição de tensões da escala macroscópica como um indicador para alterar a modelagem das regiões críticas, substituindo-se a macroescala pela mesoescala durante a análise. Consequentemente, a malha macroscópica é automaticamente substituída por uma malha mesoscópica, onde o comportamento não linear está na iminência de ocorrer. Testes numéricos são desenvolvidos para mostrar a capacidade do modelo proposto de representar o processo de iniciação e propagação de fissuras na região tracionada do concreto. Os resultados numéricos são comparados com os resultados experimentais ou com aqueles obtidos através da simulação direta em mesoescala (SDM).

Elastografia em imagens de ultrassom utilizando elementos de contorno.

Este trabalho apresenta uma nova metodologia para elastografia virtual em imagens simuladas de ultrassom utilizando métodos numéricos e métodos de visão computacional. O objetivo é estimar o módulo de elasticidade de diferentes tecidos tendo como entrada duas imagens da mesma seção transversal obtidas em instantes de tempo e pressões aplicadas diferentes. Esta metodologia consiste em calcular um campo de deslocamento das imagens com um método de fluxo óptico e aplicar um método iterativo para estimar os módulos de elasticidade (análise inversa) utilizando métodos numéricos. Para o cálculo dos deslocamentos, duas formulações são utilizadas para fluxo óptico: Lucas-Kanade e Brox. A análise inversa é realizada utilizando duas técnicas numéricas distintas: o Método dos Elementos Finitos (MEF) e o Método dos Elementos de Contorno (MEC), sendo ambos implementados em Unidades de Processamento Gráfico de uso geral, GpGPUs ( \"General Purpose Graphics Units\" ). Considerando uma quantidade qualquer de materiais a serem determinados, para a implementação do Método dos Elementos de Contorno é empregada a técnica de sub-regiões para acoplar as matrizes de diferentes estruturas identificadas na imagem. O processo de otimização utilizado para determinar as constantes elásticas é realizado de forma semi-analítica utilizando cálculo por variáveis complexas. A metodologia é testada em três etapas distintas, com simulações sem ruído, simulações com adição de ruído branco gaussiano e phantoms matemáticos utilizando rastreamento de ruído speckle. Os resultados das simulações apontam o uso do MEF como mais preciso, porém computacionalmente mais caro, enquanto o MEC apresenta erros toleráveis e maior velocidade no tempo de processamento.

Desenvolvimento teórico e experimental de FSS com elementos fractais de Gosper em estruturas de multicamadas

The fractal self-similarity property is studied to develop frequency selective surfaces (FSS) with several rejection bands. Particularly, Gosper fractal curves are used to define the shapes of the FSS elements. Due to the difficulty of making the FSS element details, the analysis is developed for elements with up to three fractal levels. The simulation was carried out using Ansoft Designer software. For results validation, several FSS prototypes with fractal elements were fabricated. In the fabrication process, fractals elements were designed using computer aided design (CAD) tools. The prototypes were measured using a network analyzer (N3250A model, Agilent Technologies). Matlab software was used to generate compare measured and simulated results. The use of fractal elements in the FSS structures showed that the use of high fractal levels can reduce the size of the elements, at the same time as decreases the bandwidth. We also investigated the effect produced by cascading FSS structures. The considered cascaded structures are composed of two FSSs separated by a dielectric layer, which distance is varied to determine the effect produced on the bandwidth of the coupled geometry. Particularly, two FSS structures were coupled through dielectric layers of air and fiberglass. For comparison of results, we designed, fabricated and measured several prototypes of FSS on isolated and coupled structures. Agreement was observed between simulated and measured results. It was also observed that the use of cascaded FSS structures increases the FSSs bandwidths and, in particular cases, the number of resonant frequencies, in the considered frequency range. In future works, we will investigate the effects of using different types of fractal elements, in isolated, multilayer and coupled FSS structures for applications on planar filters, high-gain microstrip antennas and microwave absorbers

Desenvolvimento de estruturas integradas de filtros e antenas para aplicações em sistemas de comunicações sem fio

This paper presents a study of the integration of filters and microstrip antennas, yielding devices named as filtennas for applications in wireless communications systems. The design of these structures is given from the observation of filtennas based integration between horn antennas and frequency selective surfaces (FSS), used in the band X. The choice of microstrip line structures for the development of a new configuration filtennas justifies the wide application of these transmission lines, in recent decades, always resulting in the production of circuit structures with planar light-weight, compact size, low cost, easy to construct and particularly easy to integrate with other microwave circuits. In addition, the antenna structure considered for the composition of filtennas consists of a planar monopole microstrip to microstrip filters integrated in the feed line of the antenna. In particular, are considered elliptical monopole microstrip (operating in UWB UWB) microstrip filters and (in structures with associated sections in series and / or coupled). In addition, the monopole microstrip has a proper bandwidth and omnidirectional radiation pattern, such that its integration with microstrip filters results in decreased bandwidth, but with slight changes in the radiation pattern. The methods used in the analysis of monopoles, and filters were filtennas finite elements and moments by using commercial software Ansoft Designer and HFSS Ansoft, respectively. Specifically, we analyze the main characteristics of filtennas, such as radiation pattern, gain and bandwidth. Were designed, constructed and measures, several structures filtennas, for validation of the simulated results. Were also used computational tools (CAD) in the process of building prototypes of planar monopoles, filters and filtennas. The prototypes were constructed on substrates of glass-fiber (FR4). Measurements were performed at the Laboratory for Telecommunications UFRN. Comparisons were made between simulated and measured, and found good agreement in the cases considered

Modelagem computacional para avaliação da interação entre painéis de alvenaria e estrutura de suporte em concreto armado

The developed study proposes a new computer modeling efficient and easy to apply in usual project situations to evaluate the interaction between masonry panels and support structure. The proposed model simulates the behavior of the wall exclusively using frame finite elements, thus compounding an equivalent frame. The validation was performed in two ways: firstly, through the analysis of various panels of generic plans, comparing the results obtained from equivalent frame model with the ones from a reference model, which uses shell finite elements in discretization of the walls; and in a second step, comparing with the results of the experimental model of Rosenhaupt. The analyzes considered the linear elastic behavior for materials and consisted basically in the evaluation of vertical displacements and efforts in support beams, and tensions at the base of walls. Was also evaluated, from flat and threedimensional modeling of some walls from a real project, important aspects of the wall-beam interaction, e.g.: the presence of openings of doors and windows, arranged in any position; conditions of support and linking of beams; interference of moorings between walls; and consideration of wind action. The analysis of the achieved results demonstrated the efficiency of the proposed modeling, since they have very similar aspects in the distribution of stresses and efforts, always with intensities slightly larger than those of the reference and experimental models.

Análise experimental e numérica do comportamento de estruturas de solo reforçado com geossintéticos : solo granular versus fino

As estruturas de solo reforçado com geossintéticos são normalmente constituídas por solos granulares com boas propriedades físicas e mecânicas. O uso de apenas este tipo de solos pode proporcionar o aumento, por vezes insustentável, do custo da execução das estruturas e o aumento do seu impacto ambiental. Deste modo, as estruturas de solo reforçado perdem a sua vantagem competitiva em relação a outros tipos de estruturas (muros de betão, muros de gravidade, muros de gabiões, etc.). Para resolver este problema podem ser utilizados outros tipos de solos (solos locais, finos, com propriedades físicas e mecânicas piores mas, no entanto, mais baratos) para a execução deste tipo de estruturas. De forma geral, com este estudo pretendeu-se contribuir para o incremento do conhecimento sobre a utilização de solos finos para a construção de estruturas de solo reforçado (muros e taludes). Para tal avaliaram-se as diferenças no comportamento mecânico dos materiais compósitos (solo granular reforçado versus solo fino reforçado) e das estruturas de solo reforçado constituídas com os dois tipos de solos. Assim, os objetivos deste estudo foram avaliar: a influência de vários parâmetros nas propriedades mecânicas e na capacidade de carga dos solos reforçados com geossintéticos; a influência de vários parâmetros no dimensionamento das estruturas de solo reforçado; e o comportamento das estruturas dimensionadas (incluindo a estabilidade global e a influência do processo construtivo) recorrendo a uma ferramenta numérica (PLAXIS). Para cumprir os objetivos propostos foram realizadas análises experimentais em laboratório (análise do comportamento do solo reforçado através de ensaios triaxiais e de California Bearing Ratio) e análises numéricas (dimensionamento de estruturas de solo reforçado; modelação numérica do comportamento através de uma ferramenta numérica comercial com o método dos elementos finitos). Os resultados dos ensaios experimentais mostraram que o comportamento mecânico e a capacidade de carga do solo foram incrementados com a inclusão das camadas de geossintético. Este efeito variou com os diversos parâmetros analisados mas, de forma geral, foi mais importante no solo fino (solo com propriedades mecânicas piores). As análises numéricas mostraram que as estruturas de solo fino precisaram de maior densidade de reforços para serem estáveis. Além disso, as estruturas de solo fino foram mais deformáveis e o efeito do seu processo construtivo foi mais importante (principalmente para estruturas de solo fino saturado).

Projeto de fundações e estruturas de um edifício destinado a uma igreja

Trabalho de Projeto para obtenção do grau de Mestre em Engenharia Civil na Área de Especialização em Estruturas

Análise estática e dinâmica de painéis corrugados construídos em materiais compósitos

Os painéis corrugados têm cada vez maior aplicabilidade em engenharia. São frequentemente utilizados em estruturas sandwich, em telhados e têm muitas outras aplicações nas indústrias civil, mecânica, aeronáutica e naval. Esta crescente aplicação deve-se ao facto de as suas corrugações conferirem maior rigidez ao painel devido ao aumento do rácio resistência/peso, evitando-se assim o recurso a reforços estruturais. Este facto, associado às características dos materiais compósitos laminados, que comparativamente com os materiais “tradicionais”, evidenciam uma relação rigidez/peso, bem como um comportamento mecânico; no mesmo sentido, contribui para uma potencial utilização deste tipo de solução, em variadíssimas situações na área das engenharias. Este tipo de painéis é particularmente adequado em situações de carga em que se verificam esforços de compressão e torção. Com o presente trabalho, pretende-se efectuar um estudo de análise do comportamento mecânico, quer em termos de análise estática linear, quer de vibrações livres de um painel corrugado em material compósito, sendo analisadas várias variáveis, nomeadamente: o tipo de material, a geometria da corrugação, a sequência de empilhamento das camadas de fibra, entre outras. Para este efeito, e após o enquadramento do tipo de problema que se possui, é utilizado o software comercial de análise por elementos finitos, Ansys®.

Estudo das tensões interlaminares e rotura de elementos curvos em materiais compósitos

Este trabalho aborda o comportamento mecânico de elementos curvos em materiais compósitos aquando sujeitos a forças e momentos nas extremidades. Os modos de falha mais comuns são estudados em detalhe: falha transversal na matriz devido a tensões de flexão (tensões circunferenciais), ou falha por delaminagem devido as tensões interlaminares de tracção verificadas ao longo da espessura na região curva (tensões radiais). Inicialmente, faz-se uma abordagem sobre a influência de Z-Pins na resistência interlaminar, uma abordagem sobre os efeitos de bordo, e sobre tensões interlaminares com base em teorias Layerwise de ordem superior. Salienta-se, ainda, a importância do estudo feito com o World-Wide Failure Exercise (WWFE) (Hinton et al. 2004) na comunidade científica relacionada com o tema dos materiais compósitos. O presente trabalho, é essencialmente dividido em duas partes: modelos analíticos e modelos numéricos. São estudados provetes curvos, em forma de C (semicircular), geometria analisada por Ko e Jackson (1989); e em forma de L (cantoneira de abas iguais), analisados por Sun e Kelly (1988). Começa-se por fazer uma análise teórica, através da Teoria Multicamada, ao trabalho Ko e Jackson, autores da mesma. Analisou-se também, através desta teoria, a parte em quarto de círculo da viga estudada por Sun e Kelly. Após ter a solução analítica bem compreendida, essa teoria é aplicada ao principal estudo deste trabalho, um provete semicircular, usando os empilhamentos de Sun e Kelly. É ainda estudado, e aplicado a este último, o critério de falha em três dimensões (3D) proposto por Hashin (1980). Na segunda parte, é feita uma validação numérica dos três casos anteriores, com recurso ao software comercial de elementos finitos (EF), ANSYS. Estes estudos de EF têm a diferença, relativamente aos originais, de serem em 3D, podendo assim ser retiradas todas as tensões necessárias para se aplicar o critério de falha. No final, é feita uma comparação geral entre os resultados analíticos e numéricos.

Dimensionamento de um silo em betão armado para stockagem de cimento com fundações indirectas

O presente trabalho é o corolário do curso de engenharia civil, para a obtenção do grau de Mestre, tendo-se optado pela elaboração do dimensionamento de um Silo em betão armado. O projecto inclui os elementos necessários para uma correcta execução desta estrutura. Os silos são construções, que em função da sua finalidade de utilização estão sujeitos a diversas acções peculiares. No presente trabalho foram consideradas as acções, indicadas nos Eurocódigos para este tipo de estruturas, nomeadamente a acção gerada pela presença do material ensilado, que conduz a um conjunto de pressões horizontais e verticais nas paredes do Silo, bem como a acção causada pelo gradiente de temperatura, devido ao material ensilado (120ºC). O estudo realizado para o dimensionamento do Silo, foi feito considerando a discretização do mesmo através de um programa de elementos finitos, tendo a introdução da estrutura sido realizada em computador com recurso ao “SAP2000”. Foram estudados 2 modelos: O primeiro partindo de uma configuração unifilar em que se consideraram apenas elementos ditos frames (barras), o qual serviu basicamente para o pré dimensionamento da estrutura. O segundo mais completo em que a estrutura foi analisada a partir de um conjunto que incluía a modelação de todos os elementos previamente dimensionados, com recurso a outros softwares (Robot), folhas de cálculo, tabelas e ábacos. Reuniram-se também elementos relativos a estruturas deste tipo tendo-se consultado bibliografia da especialidade, de modo a adquirir um conhecimento mais profundo necessário à realização deste projecto.

Aspectos das contenções tipo "Berlim" provisórias

O presente relatório insere-se no âmbito do trabalho final de Mestrado do curso de Engenharia Civil, área de especialização de Estruturas, do Instituto Superior de Engenharia de Lisboa. O estágio realizou-se na área de projecto do núcleo de engenharia de estruturas na empresa Mota-Engil, Engenharia e Construção, S.A. (MEEC) e teve como objectivo o desenvolvimento de competências no âmbito das contenções tipo “Berlim” provisórias utilizada nas escavações, nomeadamente, os seus aspectos de projecto e sua execução. Foi feita uma descrição geral deste tipo de estruturas de contenção onde se discutem as suas características, vantagens e inconvenientes. Descreveu-se e analisou-se um caso de estudo, onde se desenvolveram os aspectos de cálculo com recurso aos métodos clássicos da mecânica dos solos e a um programa de cálculo automático. O programa de cálculo automático utilizado foi o programa PLAXIS, no qual foram criados modelos simples de estruturas de contenção em madeira, e dos quais se obtiveram resultados relativos ao comportamento das estruturas. Procedeu-se ainda a uma breve introdução ao Eurocódigo 7 (EN 1997), referindo-se os estados limites últimos aplicados ao mesmo caso de estudo e destacando-se as principais diferenças que este veio introduzir no cálculo, procedendo-se de seguida a uma breve análise dos resultados obtidos.