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.

Estudo comparativo da resistência à tração de juntas adesivas de sobreposição simples e dupla

O método de união com ligações adesivas está cada vez mais a ser utilizado na conceção de estruturas mecânicas, por causa das vantagens significativas desta técnica em comparação com as ligações tradicionais. De facto, as juntas com ligação adesiva estão sob investigação intensa há bastante tempo. Entre as vantagens, destaca-se a redução de peso e possibilidade de unir diferentes materiais, incluindo compósitos, sem danificar as estruturas a ligar. Os adesivos comerciais variam desde resistentes e frágeis (por exemplo, Araldite® AV138) a menos resistentes e dúcteis (por exemplo, Araldite® 2015). Uma nova família de adesivos de poliuretano combina elevada resistência e ductilidade (por exemplo, Sikaforce® 7888). Este trabalho compara o desempenho à tração dos três adesivos supracitados, em juntas de alumínio (Al6082-T651) de sobreposição simples e dupla, com variação dos valores de comprimento de sobreposição (LO). A análise numérica de modelos de dano coesivo (MDC) foi realizada para analisar as tensões de arrancamento (oy) e as de corte (txy) na camada adesiva, para estudar a variável de dano do MDC durante o processo de rotura e para avaliar a capacidade MDC na previsão da resistência da junta. A análise da resistência e da variável de dano ajudou na compreensão das diferenças entre os adesivos no que se refere ao processo de rotura e resistência da junta. Observou-se que as juntas de sobreposição dupla apresentam uma distribuição de tensões bastante mais favorável relativamente às juntas de sobreposição simples, principalmente devido à eliminação da flexão do substrato interior. Como resultado, a resistência destas juntas foi tipicamente superior ao dobro da observada para as juntas de sobreposição simples, com exceção de algumas configurações de junta em que houve plastificação extensa ou mesmo rotura dos substratos por tração. O trabalho proposto permitiu também concluir que as previsões MDC são tipicamente precisas, e qual a família de adesivos é mais adequada para cada configuração de junta, com a clara vantagem para o Sikaforce® 7888. Como resultado deste trabalho, foram propostas diretrizes de conceção para juntas adesivas.

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.

Cohesive law estimation of adhesive joints in mode II condition

The adhesive bonding technique enables both weight and complexity reduction in structures that require some joining technique to be used on account of fabrication/component shape issues. Because of this, adhesive bonding is also one of the main repair methods for metal and composite structures by the strap and scarf configurations. The availability of strength prediction techniques for adhesive joints is essential for their generalized application and it can rely on different approaches, such as mechanics of materials, conventional fracture mechanics or damage mechanics. These two last techniques depend on the measurement of the fracture toughness (GC) of materials. Within the framework of damage mechanics, a valid option is the use of Cohesive Zone Modelling (CZM) coupled with Finite Element (FE) analyses. In this work, CZM laws for adhesive joints considering three adhesives with varying ductility were estimated. The End-Notched Flexure (ENF) test geometry was selected based on overall test simplicity and results accuracy. The adhesives Araldite® AV138, Araldite® 2015 and Sikaforce® 7752 were studied between high-strength aluminium adherends. Estimation of the CZM laws was carried out by an inverse methodology based on a curve fitting procedure, which enabled a precise estimation of the adhesive joints’ behaviour. The work allowed to conclude that a unique set of shear fracture toughness (GIIC) and shear cohesive strength (ts0) exists for each specimen that accurately reproduces the adhesive layer’ behaviour. With this information, the accurate strength prediction of adhesive joints in shear is made possible by CZM.

Using neural networks and support vector regression to relate marchetti dilatometer test parameters and maximum shear modulus

In the last two decades, small strain shear modulus became one of the most important geotechnical parameters to characterize soil stiffness. Finite element analysis have shown that in-situ stiffness of soils and rocks is much higher than what was previously thought and that stress-strain behaviour of these materials is non-linear in most cases with small strain levels, especially in the ground around retaining walls, foundations and tunnels, typically in the order of 10−2 to 10−4 of strain. Although the best approach to estimate shear modulus seems to be based in measuring seismic wave velocities, deriving the parameter through correlations with in-situ tests is usually considered very useful for design practice.The use of Neural Networks for modeling systems has been widespread, in particular within areas where the great amount of available data and the complexity of the systems keeps the problem very unfriendly to treat following traditional data analysis methodologies. In this work, the use of Neural Networks and Support Vector Regression is proposed to estimate small strain shear modulus for sedimentary soils from the basic or intermediate parameters derived from Marchetti Dilatometer Test. The results are discussed and compared with some of the most common available methodologies for this evaluation.

Using neural networks and support vector regression to relate marchetti dilatometer test parameters and maximum shear modulus

In the last two decades, small strain shear modulus became one of the most important geotechnical parameters to characterize soil stiffness. Finite element analysis have shown that in-situ stiffness of soils and rocks is much higher than what was previously thought and that stress-strain behaviour of these materials is non-linear in most cases with small strain levels, especially in the ground around retaining walls, foundations and tunnels, typically in the order of 10−2 to 10−4 of strain. Although the best approach to estimate shear modulus seems to be based in measuring seismic wave velocities, deriving the parameter through correlations with in-situ tests is usually considered very useful for design practice.The use of Neural Networks for modeling systems has been widespread, in particular within areas where the great amount of available data and the complexity of the systems keeps the problem very unfriendly to treat following traditional data analysis methodologies. In this work, the use of Neural Networks and Support Vector Regression is proposed to estimate small strain shear modulus for sedimentary soils from the basic or intermediate parameters derived from Marchetti Dilatometer Test. The results are discussed and compared with some of the most common available methodologies for this evaluation.

Estudo de Ferramentas FEA Comerciais na Simulação Numérica de Processos de Conformação de Chapas Metálicas. Aplicação à Indústria Automóvel.

Os componentes obtidos através da conformação plástica de chapas têm uma grande importância, tanto na etapa de concepção do produto como na etapa de produção na indústria automóvel. Isto comprova-se pelo facto de, em média, cada automóvel integrar cerca de 500 componentes estampados para construir o chassis e a carroçaria [Alves 2003]. Deste total de componentes, 50 são de grandes dimensões (portas, tejadilho, painéis inferior e laterais, entre outros) e necessitam, em média, de cinco ferramentas para o seu fabrico, sendo o custo estimado para cada ferramenta de 230 000 € [Col 2000, Alves 2003]. Para além da indústria automóvel, a conformação plástica de chapas metálicas é um processo tecnológico presente nas indústrias relativas à aeroespacial, petrolífera, decoração, alimentar, entre outras. Do ponto de vista do enquadramento económico, cerca de 20% do custo total de um automóvel novo é devido à incorporação de componentes metálicos estampados. [Alves 2003]. A pressão do “Mercado Global” faz com que os custos relativos à matéria-prima, energia e mão-de-obra sejam uma constante em termos de redução do seu impacte no orçamento das empresas fornecedoras destes produtos. É neste contexto que surge a necessidade da realização deste estudo de Benchmark de Softwares, tornando-se bastante importante, quer ao nível da competitividade industrial, quer ao nível da inovação para novos produtos. A análise por elementos finitos desempenha um papel primordial no tryout virtual e otimização das ferramentas e processos de conformação plástica. Os objetivos principais deste estudo de simulação numérica são a identificação e comparação dos resultados obtidos pelo AUTOFORM e pelo PAMSTAMP, para cada uma das variáveis identificadas como as mais influentes na robustez dos processos de estampagem de chapa metálica. Estas variáveis identificadas são: consumo de material (Draw-in) após conformação; forças de conformação; valores de variação de espessura e dos valores de extensão e resultados de Springback. Os resultados obtidos são comparados com os resultados experimentais e, desta forma, avalia-se a capacidade inovadora e a eficácia de cada um dos softwares, obtendo-se assim, uma orientação mais real para o software mais indicado aos objetivos impostos pela indústria automóvel. Para este efeito, a indústria automóvel, como maior impulsionador e motor da investigação na área da simulação numérica aplicada aos processos de estampagem, tem aderido em peso ao Benchmarking. Um exemplo disto, é o que acontece nas conferências Numisheet. O Benchmark #2 da conferência Numisheet 2008 é analisado pormenorizadamente e os resultados numéricos e experimentais são comparados e apresentados. Dois materiais distintos (aço HC260LAD e liga de alumínio AC170), assim como três modelos com geometrias diferentes (com e sem freios) são apresentados neste relatório. Com vista à redução dos ciclos tentativa-erro, tem-se adotado ciclos virtuais ou numéricos e tem-se incrementado a interatividade entre as fases de concepção e projeto, num conceito muito próprio, mas cada vez mais abrangente, denominado “produção virtual”. É nesta filosofia que se insere a simulação numérica dos processos de conformação de chapa.

Optimisation and design of post-tensioning anchorage corner blisters in concrete box girder bidges

The design of anchorage blisters of internal continuity post-tensioning tendons of bridges built by the cantilever method, presents some peculiarities, not only because they are intermediate anchorages but also because these anchorages are located in blisters, so the prestressing force has to be transferred from the blister the bottom slab and web of the girder. The high density of steel reinforcement in anchorage blisters is the most common reason for problems with concrete cast in situ, resulting in zones with low concrete compacity, leading to concrete crushing failures under the anchor plates. A solution may involve improving the concrete compression and tensile strength. To meet these requirements a high-performance fibre reinforced self-compacting mix- ture (HPFRC) was used in anchorage corner blisters of post-tensioning tendons, reducing the concrete cross-section and decreasing the reinforcement needed. To assess the ultimate capacity and the adequate serviceability of the local anchorage zone after reducing the minimum concrete cross-section and the confining reinforcement, specified by the anchorage device supplier for the particular tendon, load transfer tests were performed. To investigate the behaviour of anchorage blisters regarding the transmission of stresses to the web and the bottom slab of the girder, and the feasibility of using high performance concrete only in the blister, two half scale models of the inferior corner of a box girder existing bridge were studied: a reference specimen of ordinary reinforced concrete and a HPFRC blister specimen. The design of the reinforcement was based in the tensile forces obtained on strut-and-tie models. An experimental program was carried out to assess the models used in design and to study the feasibility of using high performance concrete only in the blister, either with casting in situ, or with precast solutions. A non-linear finite element analysis of the tested specimens was also performed and the results compared.

Electro-Optical/Infrared sensor turret integration on an aircraft - structural im-pact on LOCKHEED MARTIN C-130H and design methodology

The growing need to patrol and survey large maritime and terrestrial areas increased the need to integrate external sensors on aircraft in order to accomplish those patrols at increasingly higher altitudes, longer range and not depending upon vehicle type. The main focus of this work is to elaborate a practical, simple, effective and efficient methodology for the aircraft modification procedure resulting from the integration of an Elec-tro-Optical/Infra-Red (EO/IR) turret through a support structure. The importance of the devel-opment of a good methodology relies on the correct management of project variables as time, available resources and project complexity. The key is to deliver a proper tool for a project de-sign team that will be used to create a solution that fulfils all technical, non-technical and certi-fication requirements present in this field of transportation. The created methodology is inde-pendent of two main inputs: sensor model and aircraft model definition, and therefore it is in-tended to deliver the results for different projects besides the one that was presented in this work as a case study. This particular case study presents the development of a structure support for FLIR STAR SAPHIRE III turret integration on the front lower fuselage bulkhead (radome) of the LOCKHEED MARTIN C-130 H. Development of the case study focuses on the study of local structural analysis through the use of Finite Element Method (FEM). Development of this Dissertation resulted in a cooperation between Faculty of Science and Technology - Universidade Nova de Lisboa and the company OGMA - Indústria Aeronáutica de Portugal

Fatigue life assessment of endodontic instruments subjected to multi planar curvature

This present study aimed to investigate the fatigue life of unused (new) endodontic instruments made of NiTi with control memory by Coltene™ and subjected to the multi curvature of a mandibular first molar root canal. Additionally, the instrument‟s structural behaviour was analysed through non-linear finite element analysis (FEA). The fatigue life of twelve Hyflex™ CM files was assessed while were forced to adopt a stance with multiple radius of curvature, similar to the ones usually found in a mandibular first molar root canal; nine of them were subjected to Pecking motion, a relative movement of axial type. To achieve this, it was designed an experimental setup with the aim of timing the instruments until fracture while worked inside a stainless steel mandibular first molar model with relative axial motion to simulate the pecking motion. Additionally, the model‟s root canal multi-curvature was confirmed by radiography. The non-linear finite element analysis was conducted using the computer aided design software package SolidWorks™ Simulation, in order to define the imposed displacement required by the FEA, it was necessary to model an endodontic instrument with simplified geometry using SolidWorks™ and subsequently analyse the geometry of the root canal CAD model. The experimental results shown that the instruments subjected to pecking motion displayed higher fatigue life values and higher lengths of fractured tips than those with only rotational relative movement. The finite element non-linear analyses shown, for identical conditions, maximum values for the first principal stress lower than the yield strength of the material and those were located in similar positions to the instrument‟s fracture location determined by the experimental testing results.

Development of a new U 200 pultrusion die to produce GRP composite profiles using different thermosetting resins

In this work, a steel heated pultrusion die was designed, developed and manufactured to produce U200 glass fibre reinforced thermosetting matrix (GRP) profiles. The finite element analysis (FEA) was used to predict and optimise the developed die heating by using cylindrical electrical powered cartridges. To assess the new die performance it was mounted in the 120 kN pultrusion line of the Portuguese company Vidropol SA and used to produce continuously U200 profiles able to meet all requirements specified for the E23 grade accordingly to the European Standard EN 13706: 2002. After setting up the type, orientation and sequence of layers in laminate, orthophthalic, isophthalic and bisphenolic unsaturated polyester as well as vinylester resins were used to produce glass fibre reinforced U 200 composite profiles. An appropriated catalyst system was selected and the processing variables optimised for each case, namely, pultrusion pull-speed and die temperature. Finally, the produced U200 profiles were submitted to visual inspection, calcination and mechanical tests, namely, flexural, tensional and interlaminar shear strength (ILSS) tests, to assess their accomplishment with the EN 13706 requirements.

Using a new developed die and different types of thermosetting resins to produce U200 profiles by pultrusion

In this work, a new steel heated pultrusion die was designed, developed and manufactured to produce U200 glass fibre reinforced thermosetting matrix (GRP) profiles. The finite element analysis (FEA) was used to predict and optimise the developed die heating by using cylindrical electrical powered cartridges. To assess the new die performance it was mounted in the 120 kN pultrusion line of the Portuguese company Vidropol SA and used to produce continuously U200 profiles able to meet all requirements specified for the E23 grade accordingly to the European Standard EN 13706: 2002. After setting up the type, orientation and sequence of layers in the U 200 laminate, different types of thermosetting resins were used in its production. Orthophthalic, isophthalic and bisphenolic unsaturated polyester as well as vinylester resins were used to produce glass fibre reinforced U 200 composite profiles. All applied resins were submitted to SPI gel tests in order to select the more appropriated catalyst system and optimise the processing variables to be used in each case, namely, pultrusion pull-speed and die temperature. The best pultrusion operational conditions were selected by varying and monitoring the pull-speed and die temperature and, at the same time, measuring the temperature on the manufactured U 200 profile during processing. Finally, the produced U200 profiles were submitted to visual inspection, calcination and mechanical tests, namely, flexural, tensional and interlaminar shear strength (ILSS) tests, to assess their accomplishment with the EN 13706 requirements.

Evaluation of temperature influence on embedded kinematic system

Dissertação de mestrado integrado em Engenharia Mecânica

Development of procedures for the design, optimization and manufacturing of customized orthopaedic and trauma implants: Geometrical/anatomical modelling from 3D medical imaging

Tese de Doutoramento (Programa Doutoral em Engenharia Biomédica)

Optimising the design and manufacture of a customised, prescription fit, Maxillo-facial prosthesis

This thesis is a continuation of the Enterprise-Ireland Research Innovation Fund (RIF) Project entitled’ "Design and Manufacturing of Customised Maxillo-Facial Prostheses" The primary objective of this Internal Research Development Program (IRDP) project was to investigate two fundamental design changes 1 To incorporate the over-denture abutments directly into the implant. 2 To remove the restraining wings by the addition of screws, which affix the. implant to the dense material of the jawbone. The prosthetic was redesigned using the ANSYS Finite Element Analysis software program and analysed to* • Reduce the internal von Mises stress distribution The new prosthetic had a -63.63 % lower von Mises stress distribution when compared with the original prosthetic. • Examine the screw preload effects. A maximum relative displacement of 22 6 * lO^mm between the bone and screw was determined, which is well below the critical threshold of micromotion which prevents osseointegration • Investigate the prosthetic-bone contact interface. Three models of the screw, prosthesis, and bone, were studied. (Axisymmetnc, quarter volume, and full volume), a recommended preload torque of 0 32 Nm was applied to the prosthetic and a maximum von Mises stress of 1.988 MPa was predicted • Study the overdenture removal forces. This analysis could not be completed because the correct plastic multilinear properties of the denture material could not be established The redesigned prosthetic was successfully manufactured on a 3-axis milling machine with an indexing system The prosthetic was examined for dimensional quality and strength The research established the feasibility of the new design and associated manufacturing method.