EFFECT – Efficient finite element code

The theme of this dissertation is the finite element method applied to mechanical structures. A new finite element program is developed that, besides executing different types of structural analysis, also allows the calculation of the derivatives of structural performances using the continuum method of design sensitivities analysis, with the purpose of allowing, in combination with the mathematical programming algorithms found in the commercial software MATLAB, to solve structural optimization problems. The program is called EFFECT – Efficient Finite Element Code. The object-oriented programming paradigm and specifically the C ++ programming language are used for program development. The main objective of this dissertation is to design EFFECT so that it can constitute, in this stage of development, the foundation for a program with analysis capacities similar to other open source finite element programs. In this first stage, 6 elements are implemented for linear analysis: 2-dimensional truss (Truss2D), 3-dimensional truss (Truss3D), 2-dimensional beam (Beam2D), 3-dimensional beam (Beam3D), triangular shell element (Shell3Node) and quadrilateral shell element (Shell4Node). The shell elements combine two distinct elements, one for simulating the membrane behavior and the other to simulate the plate bending behavior. The non-linear analysis capability is also developed, combining the corotational formulation with the Newton-Raphson iterative method, but at this stage is only avaiable to solve problems modeled with Beam2D elements subject to large displacements and rotations, called nonlinear geometric problems. The design sensitivity analysis capability is implemented in two elements, Truss2D and Beam2D, where are included the procedures and the analytic expressions for calculating derivatives of displacements, stress and volume performances with respect to 5 different design variables types. Finally, a set of test examples were created to validate the accuracy and consistency of the result obtained from EFFECT, by comparing them with results published in the literature or obtained with the ANSYS commercial finite element code.

POST-PROCESSING TECHNIQUES SUITABILITY FOR MESOLEVEL FREE BOUNDARY FLOWS

Reliable flow simulation software is inevitable to determine an optimal injection strategy in Liquid Composite Molding processes. Several methodologies can be implemented into standard software in order to reduce CPU time. Post-processing techniques might be one of them. Post-processing a finite element solution is a well-known procedure, which consists in a recalculation of the originally obtained quantities such that the rate of convergence increases without the need for expensive remeshing techniques. Post-processing is especially effective in problems where better accuracy is required for derivatives of nodal variables in regions where Dirichlet essential boundary condition is imposed strongly. In previous works influence of smoothness of non-homogeneous Dirichlet condition, imposed on smooth front was examined. However, usually quite a non-smooth boundary is obtained at each time step of the infiltration process due to discretization. Then direct application of post-processing techniques does not improve final results as expected. The new contribution of this paper lies in improvement of the standard methodology. Improved results clearly show that the recalculated flow front is closer to the ”exact” one, is smoother that the previous one and it improves local disturbances of the “exact” solution.

Post-processing Techniques in Free Boundary Flows during Liquid Composite Moulding Processes

Post-processing a finite element solution is a well-known technique, which consists in a recalculation of the originally obtained quantities such that the rate of convergence increases without the need for expensive remeshing techniques. Postprocessing is especially effective in problems where better accuracy is required for derivatives of nodal variables in regions where Dirichlet essential boundary condition is imposed strongly. Consequently such an approach can be exceptionally good in modelling of resin infiltration under quasi steady-state assumption by remeshing techniques and with explicit time integration, because only the free-front normal velocities are necessary to advance the resin front to the next position. The new contribution is the post-processing analysis and implementation of the freeboundary velocities of mesolevel infiltration analysis. Such implementation ensures better accuracy on even coarser meshes, which in consequence reduces the computational time also by the possibility of employing larger time steps.

Critical Velocity obtained using Simplified Models of the Railway Track: Viability and Applicability

Increased demands on the capacity of the railway network gave rise to new issues related to the dynamic response of railway tracks subjected to moving vehicles. Thus, it becomes important to evaluate the applicability of traditionally used simplified models which have a closed form solution. Regarding simplified models, transversal vibrations of a beam on a visco-elastic foundation subjected to a moving load are considered. Governing equations are obtained by Hamilton’s principle. Shear distortion, rotary inertia and effect of axial force are accounted for. The load is introduced as a time varying force moving at a constant velocity. Transversal vibrations induced by the load are solved by the normal-mode analysis. Reflected waves at the extremities of the full beam are avoided by introduction of semi-infinite elements. Firstly, the critical velocity obtained from this model is compared with results of an undamped Euler- Bernoulli formulation with zero axial force. Secondly, a finite element model in ABAQUS is examined. The new contribution lies in the introduction of semi- infinite elements and in the first step to a systematic comparison, which have not been published so fa

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

The influence of geometric factors on the wall shear stress distribution in realistic human coronary arteries

Dissertação apresentada na Faculdade de Ciências e Tecnologia da Universidade Nova de Lisboa para a obtenção do grau de Mestre em Engenharia Biomédica. A presente dissertação foi desenvolvida no Erasmus Medical Center em Roterdão, Holanda

The Influence of Damping on Vibration Induced by High-Speed Trains

Passage of high-speed trains may induce high ground and track vibrations, which, besides increasing wheel, rail and track deterioration, may have a negative impact on the vehicle stability and on the passengers comfort. In this paper two distinct analyses are presented. The first one is dedicated to efficient decoupling of rail and soil vibrations by suggesting new interface materials in rail-sleeper fixing system, i.e. in the part where damping efficiency can be directly controlled and tested. The second analysis concerns with an adequate model of soils damping. Proper understanding and correct numerical simulation of this behaviour can help in suggesting soil improvement techniques.

Modelling a cable structure for a new wind energy production device

Dissertação apresentada na Faculdade de Ciências e Tecnologia da Universidade Nova de Lisboa para obtenção do grau de Mestre em Engenharia Mecânica Especialização em Concepção e Produção

Robustness of corroded reinforced concrete structures. A structural performance approach

Structure and Infrastructure Engineering, 1-17

Fatigue life assessment of exhaust systems for naval gas turbines

Dissertação para obtenção do Grau de Mestre em Engenharia Mecânica

An application of reliability-based robustness assessment of steel moment resisting frame structures under post-mainshock cascading events

The paper presented herein proposes a reliability-based framework for quantifying the structural robustness considering the occurrence of a major earthquake (mainshock) and subsequent cascading hazard events, such as aftershocks that are triggered by the mainshock. These events can significantly increase the probability of failure of buildings, especially for structures that are damaged during the mainshock. The application of the proposed framework is exemplified through three numerical case studies. The case studies correspond to three SAC steel moment frame buildings of 3-, 9-, and 20- stories, which were designed to pre-Northridge codes and standards. Twodimensional nonlinear finite element models of the buildings are developed using the Open System for Earthquake Engineering Simulation framework (OpenSees), using a finite-length plastic hinge beam model and a bilinear constitutive law with deterioration, and are subjected to multiple mainshock-aftershock seismic sequences. For the three buildings analyzed herein, it is shown that the structural reliability under a single seismic event can be significantly different from that under a sequence of seismic events. The reliability-based robustness indicator used shows that the structural robustness is influenced by the extent by which a structure can distribute damage.

Identificação modal e actualização de modelos de elementos finitos

Dissertação para obtenção do Grau de Mestre em Engenharia Civil - Perfil Estruturas

Damage propagation in composite Materials meso-mechanical models

Composite materials have a complex behavior, which is difficult to predict under different types of loads. In the course of this dissertation a methodology was developed to predict failure and damage propagation of composite material specimens. This methodology uses finite element numerical models created with Ansys and Matlab softwares. The methodology is able to perform an incremental-iterative analysis, which increases, gradually, the load applied to the specimen. Several structural failure phenomena are considered, such as fiber and/or matrix failure, delamination or shear plasticity. Failure criteria based on element stresses were implemented and a procedure to reduce the stiffness of the failed elements was prepared. The material used in this dissertation consist of a spread tow carbon fabric with a 0°/90° arrangement and the main numerical model analyzed is a 26-plies specimen under compression loads. Numerical results were compared with the results of specimens tested experimentally, whose mechanical properties are unknown, knowing only the geometry of the specimen. The material properties of the numerical model were adjusted in the course of this dissertation, in order to find the lowest difference between the numerical and experimental results with an error lower than 5% (it was performed the numerical model identification based on the experimental results).

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.