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.

Prestress losses in reinforced concrete beams strengthened with NSM-CFRP laminates

Using prestressed near surface mounted fibre reinforced polymers (NSM-FRP) is nowadays regaining the attention from the scientific community for the strengthening of existing reinforced concrete (RC) structures. The application of prestressed internal FRP bars and externally bonded prestressed FRPs has already been deeply investigated and revealed considerable benefits when compared to the corresponding passive solutions. A certain amount of prestress provides benefits mainly associated to structural integrity and material durability. Immediately after prestress transference, it is possible to close some of the existing cracks, decreasing the susceptibility of the element to corrosion and, a certain amount of deflection can be recovered due to the creation of a negative curvature. However, very few studies have been carried out to properly assess the preservation of prestress over time. In this context, several reinforced concrete beams strengthened with prestressed NSM carbon FRP (CFRP) laminates were prestressed and monitored for about 40 days. The data obtained from these experimental programs is in this paper presented and analysed. The observed prestress losses were later modelled using finite elements analysis and, although this topic is not addressed in this paper, the obtained results revealed considerable precision. The largest strain losses in the CFRP laminate were found to be mainly located in the extremities of the bonded length, while in the central zone most of the applied pre-strain was retained over time. The highest CFRP strain losses were observed in the first 6 to 12 days after prestress transfer, suggesting that the application of prestressed NSM-FRP will be very effective over time.

Deteção de colisões com guardas de segurança

Dissertação de mestrado integrado em Engenharia Eletrónica Industrial e Computadores

Caracteritació de models funcionals d'àmfores romanes mitjançant Finite Element Analysis

Projecte de recerca elaborat a partir d’una estada al Laboratory of Archaeometry del National Centre of Scientific Research “Demokritos” d’Atenes, Grècia, entre juny i setembre 2006. Aquest estudi s’emmarca dins d’un context més ampli d’estudi del canvi tecnològic que es documenta en la producció d’àmfores de tipologia romana durant els segles I aC i I dC en els territoris costaners de Catalunya. Una part d’aquest estudi contempla el càlcul de les propietats mecàniques d’aquestes àmfores i la seva avaluació en funció de la tipologia amforal, a partir de l’Anàlisi d’Elements Finits (AEF). L’AEF és una aproximació numèrica que té el seu origen en les ciències d’enginyeria i que ha estat emprada per estimar el comportament mecànic d’un model en termes, per exemple, de deformació i estrès. Així, un objecte, o millor dit el seu model, es dividit en sub-dominis anomenats elements finits, als quals se’ls atribueixen les propietats mecàniques del material en estudi. Aquests elements finits estan connectats formant una xarxa amb constriccions que pot ser definida. En el cas d’aplicar una força determinada a un model, el comportament de l’objecte pot ser estimat mitjançant el conjunt d’equacions lineals que defineixen el rendiment dels elements finits, proporcionant una bona aproximació per a la descripció de la deformació estructural. Així, aquesta simulació per ordinador suposa una important eina per entendre la funcionalitat de ceràmiques arqueològiques. Aquest procediment representa un model quantitatiu per predir el trencament de l’objecte ceràmic quan aquest és sotmès a diferents condicions de pressió. Aquest model ha estat aplicat a diferents tipologies amforals. Els resultats preliminars mostren diferències significatives entre la tipologia pre-romana i les tipologies romanes, així com entre els mateixos dissenys amforals romans, d’importants implicacions arqueològiques.

Estrategies per la inversió tridimensional de dades magnetotel.lúriques:

Projecte de recerca elaborat a partir d’una estada a la Dublin Institute for Advanced Studies, Irlanda, entre setembre i desembre del 2009.En els últims anys s’ha realitzat un important avanç en la modelització tridimensional en magnetotel•lúrica (MT) gracies a l'augment d’algorismes d’inversió tridimensional disponibles. Aquests codis utilitzen diferents formulacions del problema (diferències finites, elements finits o equacions integrals), diverses orientacions del sistema de coordenades i, o bé en el conveni de signe, més o menys, en la dependència temporal. Tanmateix, les impedàncies resultants per a tots els valors d'aquests codis han de ser les mateixes una vegada que es converteixen a un conveni de signe comú i al mateix sistema de coordenades. Per comparar els resultats dels diferents codis hem dissenyat models diferents de resistivitats amb estructures tridimensional incrustades en un subsòl homogeni. Un requisit fonamental d’aquests models és que generin impedàncies amb valors importants en els elements de la diagonal, que no són menyspreables. A diferència dels casos del modelització de dades magnetotel.lúriques unidimensionals i bidimensionals, pel al cas tridimensional aquests elements de les diagonals del tensor d'impedància porten informació sobre l'estructura de la resistivitat. Un dels models de terreny s'utilitza per comparar els diferents algoritmes que és la base per posterior inversió dels diferents codis. Aquesta comparació va ser seguida de la inversió per recuperar el conjunt de dades d'una estructura coneguda.

Influence of the implanted pulse generator as reference electrode in finite element model of monopolar deep brain stimulation.

Electrical deep brain stimulation (DBS) is an efficient method to treat movement disorders. Many models of DBS, based mostly on finite elements, have recently been proposed to better understand the interaction between the electrical stimulation and the brain tissues. In monopolar DBS, clinically widely used, the implanted pulse generator (IPG) is used as reference electrode (RE). In this paper, the influence of the RE model of monopolar DBS is investigated. For that purpose, a finite element model of the full electric loop including the head, the neck and the superior chest is used. Head, neck and superior chest are made of simple structures such as parallelepipeds and cylinders. The tissues surrounding the electrode are accurately modelled from data provided by the diffusion tensor magnetic resonance imaging (DT-MRI). Three different configurations of RE are compared with a commonly used model of reduced size. The electrical impedance seen by the DBS system and the potential distribution are computed for each model. Moreover, axons are modelled to compute the area of tissue activated by stimulation. Results show that these indicators are influenced by the surface and position of the RE. The use of a RE model corresponding to the implanted device rather than the usually simplified model leads to an increase of the system impedance (+48%) and a reduction of the area of activated tissue (-15%).

A temperature predicting model for manufacturing processes requiring coiling

A model for predicting temperature evolution for automatic controling systems in manufacturing processes requiring the coiling of bars in the transfer table is presented. Although the method is of a general nature, the presentation in this work refers to the manufacturing of steel plates in hot rolling mills. The predicting strategy is based on a mathematical model of the evolution of temperature in a coiling and uncoiling bar and is presented in the form of a parabolic partial differential equation for a shape changing domain. The mathematical model is solved numerically by a space discretization via geometrically adaptive finite elements which accomodate the change in shape of the domain, using a computationally novel treatment of the resulting thermal contact problem due to coiling. Time is discretized according to a Crank-Nicolson scheme. Since the actual physical process takes less time than the time required by the process controlling computer to solve the full mathematical model, a special predictive device was developed, in the form of a set of least squares polynomials, based on the off-line numerical solution of the mathematical model.

Numerical study of Petrov-Galerkin formulations for the shallow water wave equations

The behavior of Petrov-Galerkin formulations for shallow water wave equations is evaluated numerically considering typical one-dimensional propagation problems. The formulations considered here use stabilizing operators to improve classical Galerkin approaches. Their advantages and disadvantages are pointed out according to the intrinsic time scale (free parameter) which has a particular importance in this kind of problem. The influence of the Courant number and the performance of the formulation in dealing with spurious oscillations are adressed.

Numerical Investigations on Soil Structure Interaction of Multistorey Frames

Frames are the most widely used structural system for multistorey buildings. A building frame is a three dimensional discrete structure consisting of a number of high rise bays in two directions at right angles to each other in the vertical plane. Multistorey frames are a three dimensional lattice structure which are statically indeterminate. Frames sustain gravity loads and resist lateral forces acting on it. India lies at the north westem end of the Indo-Australian tectonic plate and is identified as an active tectonic area. Under horizontal shaking of the ground, horizontal inertial forces are generated at the floor levels of a multistorey frame. These lateral inertia forces are transferred by the floor slab to the beams, subsequently to the columns and finally to the soil through the foundation system. There are many parameters that affect the response of a structure to ground excitations such as, shape, size and geometry of the structure, type of foundation, soil characteristics etc. The Soil Structure Interaction (SS1) effects refer to the influence of the supporting soil medium on the behavior of the structure when it is subjected to different types of loads. Interaction between the structure and its supporting foundation and soil, which is a complete system, has been modeled with finite elements. Numerical investigations have been carried out on a four bay, twelve storeyed regular multistorey frame considering depth of fixity at ground level, at characteristic depth of pile and at full depth. Soil structure interaction effects have been studied by considering two models for soil viz., discrete and continuum. Linear static analysis has been conducted to study the interaction effects under static load. Free vibration analysis and further shock spectrum analysis has been conducted to study the interaction effects under time dependent loads. The study has been extended to four types of soil viz., laterite, sand, alluvium and layered.The structural responses evaluated in the finite element analysis are bending moment, shear force and axial force for columns, and bending moment and shear force for beams. These responses increase with increase in the founding depth; however these responses show minimal increase beyond the characteristic length of pile. When the soil structure interaction effects are incorporated in the analysis, the aforesaid responses of the frame increases upto the characteristic depth and decreases when the frame has been analysed for the full depth. It has been observed that shock spectrum analysis gives wide variation of responses in the frame compared to linear elastic analysis. Both increase and decrease in responses have been observed in the interior storeys. The good congruence shown by the two finite element models viz., discrete and continuum in linear static analysis has been absent in shock spectrum analysis.

A finite element method for second order nonvariational elliptic problems

We propose a numerical method to approximate the solution of second order elliptic problems in nonvariational form. The method is of Galerkin type using conforming finite elements and applied directly to the nonvariational (nondivergence) form of a second order linear elliptic problem. The key tools are an appropriate concept of “finite element Hessian” and a Schur complement approach to solving the resulting linear algebra problem. The method is illustrated with computational experiments on three linear and one quasi-linear PDE, all in nonvariational form.

An improved finite element space for discontinuous pressures

We consider incompressible Stokes flow with an internal interface at which the pressure is discontinuous, as happens for example in problems involving surface tension. We assume that the mesh does not follow the interface, which makes classical interpolation spaces to yield suboptimal convergence rates (typically, the interpolation error in the L(2)(Omega)-norm is of order h(1/2)). We propose a modification of the P(1)-conforming space that accommodates discontinuities at the interface without introducing additional degrees of freedom or modifying the sparsity pattern of the linear system. The unknowns are the pressure values at the vertices of the mesh and the basis functions are computed locally at each element, so that the implementation of the proposed space into existing codes is straightforward. With this modification, numerical tests show that the interpolation order improves to O(h(3/2)). The new pressure space is implemented for the stable P(1)(+)/P(1) mini-element discretization, and for the stabilized equal-order P(1)/P(1) discretization. Assessment is carried out for Poiseuille flow with a forcing surface and for a static bubble. In all cases the proposed pressure space leads to improved convergence orders and to more accurate results than the standard P(1) space. In addition, two Navier-Stokes simulations with moving interfaces (Rayleigh-Taylor instability and merging bubbles) are reported to show that the proposed space is robust enough to carry out realistic simulations. (c) 2009 Elsevier B.V. All rights reserved.

Um modelo numerico para analise de estruturas de materiais compostos considerando efeitos viscoelasticos e falhas progressivas

Neste trabalho apresentam-se pocedimentos para análise não linear de estruturas de materiais compostos laminados reforçados por fibras. A formulação é baseada em uma descrição cinemática incremental Lagrangeana Total, que permite o tratamento de deslocamentos arbitrariamente grandes com pequenas deformações, utilizando elementos finitos tridimensionais degenerados deduzidos para a análise de cascas. As estruturas são consideradas como submetidas a cargas mecânicas e a ações de temperatura e de umidade. O material é suposto elástico linear com propriedades dependentes, ou não, dos valores da temperatura e da concentração de umidade, ou viscoelástico linear com uma relação constitutiva em integral hereditária , e com comportamento higrotermo-reologicamente simples. As lâminas são consideradas como sujeitas a falhas, as quais são detectadas através de critérios macroscópicos, baseados em tensões ou em deformações. As equações não lineares de equilíbrio são resolvidas através de procedimentos iterativos e as deformações dependentes do tempo são avaliadas pelo método das variáveis de estado. Diversos exemplos numéricos de estruturas submetidas à flexão, flambagem elástica e viscoelástica e falhas são apresentados.

Finite Element Modeling for Development and Optimization of a Bone Plate for Mandibular Fracture in Dogs

This study aimed to develop a plate to treat fractures of the mandibular body in dogs and to validate the project using finite elements and biomechanical essays. Mandible prototypes were produced with 10 oblique ventrorostral fractures (favorable) and 10 oblique ventrocaudal fractures (unfavorable). Three groups were established for each fracture type. Osteosynthesis with a pure titanium plate of double-arch geometry and blocked monocortical screws offree angulanon were used. The mechanical resistance of the prototype with unfavorable fracture was lower than that of the fcworable fracture. In both fractures, the deflection increased and the relative stiffness decreased proportionally to the diminishing screw number The finite element analysis validated this plate study, since the maximum tension concentration observed on the plate was lower than the resistance limit tension admitted by the titanium. In conclusion, the double-arch geometry plate fixed with blocked monocortical screws has sufficient resistance to stabilize oblique,fractures, without compromising mandibular dental or neurovascular structures. J Vet Dent 24 (7); 212 - 221, 2010