Surgical correction of scoliosis: Numerical analysis and optimization of the procedure

A previously developed model is used to numerically simulate real clinical cases of the surgical correction of scoliosis. This model consists of one-dimensional finite elements with spatial deformation in which (i) the column is represented by its axis; (ii) the vertebrae are assumed to be rigid; and (iii) the deformability of the column is concentrated in springs that connect the successive rigid elements. The metallic rods used for the surgical correction are modeled by beam elements with linear elastic behavior. To obtain the forces at the connections between the metallic rods and the vertebrae geometrically, non-linear finite element analyses are performed. The tightening sequence determines the magnitude of the forces applied to the patient column, and it is desirable to keep those forces as small as possible. In this study, a Genetic Algorithm optimization is applied to this model in order to determine the sequence that minimizes the corrective forces applied during the surgery. This amounts to find the optimal permutation of integers 1, ... , n, n being the number of vertebrae involved. As such, we are faced with a combinatorial optimization problem isomorph to the Traveling Salesman Problem. The fitness evaluation requires one computing intensive Finite Element Analysis per candidate solution and, thus, a parallel implementation of the Genetic Algorithm is developed.

Comporta de vigas: dimensionamento e análise estrutural

As comportas de vigas de fundo, objecto de estudo no presente trabalho, são essenciais devido à sua capacidade de controlo do fluxo de água em casos de emergência. Assim este trabalho visou o dimensionamento e análise estrutural de uma comporta dessa natureza. Para o seu dimensionamento foi essencial a norma DIN19704-1:1998, que define todos os padrões que devem ser levados em conta numa construção hidráulica em aço. Deste modo, após uma breve descrição do enquadramento e do estado da arte relativamente a este assunto, foram apresentados ao longo do trabalho, os cálculos dos principais componentes mecânicos que fazem parte dum projecto de uma comporta de vigas. Para comprovar a validade do dimensionamento da comporta realizado através da norma, recorreu-se a um software de análise e simulação por elementos finitos, COSMOSWorks, no sentido de assim prever o comportamento mecânico em análise estática linear, que as solicitações mecânicas em presença têm sobre a comporta. Os resultados da comparação entre o dimensionamento teórico e a análise de tensões através do COSMOSWorks nos elementos mais críticos da estrutura, permitem concluir globalmente que a comporta está bem dimensionada tendo em conta os esforços que tem de suportar.

Influência da modelação de núcleos de betão armado no cálculo automático de estruturas de edifícios

O presente trabalho aborda a questão do comportamento e influência estrutural da modelação de um núcleo de secção aberta de parede fina no modelo tridimensional de elementos finitos de um edifício de betão armado, composto por 5 pisos sujeito à acção sísmica regulamentar. Para o efeito, foram considerados os Eurocódigos estruturais, nomeadamente o Eurocódigo 0, Eurocódigo 1 e Eurocódigo 8. A fim de avaliar a importância do tipo de modelação destes elementos estruturais, elemento fundamental para a caracterização da resposta dinâmica, foi desenvolvido um estudo paramétrico, fazendo-se a análise de quatro formas distintas de modelação do mesmo núcleo. Para os mesmos casos de carregamento, é feita uma análise comparativa entre os resultados obtidos a nível de períodos, percentagem de armadura nos pilares e no próprio núcleo bem como os deslocamentos para os modelos em estudo, com recurso a um programa de cálculo automático. O edifício considerado neste trabalho foi modelado com base no programa comercial de cálculo automático SAP 2000 versão 14.4. Posteriormente apresentam-se e discutem-se os resultados obtidos nos vários modelos em causa e para as várias situações propostas ao longo do trabalho.

Projecto de fundações e estrutura de um edifício de 4 pisos em Lisboa, destinado a habitação

O objectivo deste trabalho consiste na elaboração de um projecto de estruturas de um edifício, desde a fase de concepção estrutural até à fase de pormenorização, passando obviamente pelo dimensionamento. O trabalho foi motivo para aplicação dos conhecimentos adquiridos durante o curso, na prática de projecto de estruturas, nomeadamente na utilização da regulamentação europeia de estruturas que está para entrar em vigor. O projecto foi elaborado a partir do projecto base de arquitectura e diz respeito a um conjunto de 4 edifícios de habitação separados por junta de dilatação, a construir em Lisboa. A estrutura em betão armado constituída por lajes fungiformes maciças apoiadas em pilares, paredes de betão e vigas, foi analisada com recurso a um programa de cálculo automático de elementos finitos SAP 200 V14.2.4. Para a execução deste trabalho recorreu-se à regulamentação europeia, nomeadamente o Eurocódigo 0, o Eurocódigo 1, o Eurocódigo 2 e o Eurocódigo 8, na sua versão final de Norma Portuguesa. Os critérios de dimensionamento bem como as disposições construtivas relativas a armaduras de betão são detalhadamente referidos ao longo do trabalho. Embora a análise estrutural tenha sido efectuada para quatro edifícios, os desenhos de pormenorização de betão armado são apresentados relativamente apenas ao edifício Nº1. Optou-se por esta solução em virtude de os edifícios serem semelhantes.

Determination of stresses around a cylindrical single pile caused by horizontal movements of soils with the finitive element method

This study was carried out with the aim of modeling in 2D, in plain strain, the movement of a soft cohesive soil around a pile, in order to enable the determination of stresses resulting along the pile, per unit length. The problem in study fits into the large deformations problem and can be due to landslide, be close of depth excavations, to be near of zones where big loads are applied in the soil, etc. In this study is used an constitutive Elasto-Plastic model with the failure criterion of Mohr-Coulomb to model the soil behavior. The analysis is developed considering the soil in undrained conditions. To the modeling is used the finite element program PLAXIS, which use the Updated Lagrangian - Finite Element Method (UL-FEM). In this work, special attention is given to the soil-pile interaction, where is presented with some detail the formulation of the interface elements and some studies for a better understand of his behavior. It is developed a 2-D model that simulates the effect of depth allowing the study of his influence in the stress distribution around the pile. The results obtained give an important base about how behaves the movement of the soil around a pile, about how work the finite element program PLAXIS and how is the stress distribution around the pile. The analysis demonstrate that the soil-structure interaction modeled with the UL-FEM and interface elements is more appropriate to small deformations problems.

Lar militar: projecto de execução

Trabalho Final de Mestrado para obtenção do grau de mestre em Engenharia Civil na Área de Especialização em Estruturas

Avaliação experimental do comportamento dinâmico de um edifício do século XIX: aplicação à calibração de modelos numéricos

Trabalho de Dissertação de natureza científica para obtenção do grau de Mestre em Engenharia Civil na Área de Especialização em Estruturas

Determinação de propriedades mecânicas através de ensaios de provetes miniatura, "Small Punch"

Trabalho Final de Mestrado para obtenção do grau de Mestre em Engenharia Mecânica

Colloidal dipolar interactions in 2D smectic-C films

We use a two-dimensional (2D) elastic free energy to calculate the effective interaction between two circular disks immersed in smectic-C films. For strong homeotropic anchoring, the distortion of the director field caused by the disks generates topological defects that induce an effective interaction between the disks. We use finite elements, with adaptive meshing, to minimize the 2D elastic free energy. The method is shown to be accurate and efficient for inhomogeneities on the length scales set by the disks and the defects, that differ by up to 3 orders of magnitude. We compute the effective interaction between two disk-defect pairs in a simple (linear) configuration. For large disk separations, D, the elastic free energy scales as similar to D-2, confirming the dipolar character of the long-range effective interaction. For small D the energy exhibits a pronounced minimum. The lowest energy corresponds to a symmetrical configuration of the disk-defect pairs, with the inner defect at the mid-point between the disks. The disks are separated by a distance that, is twice the distance of the outer defect from the nearest disk. The latter is identical to the equilibrium distance of a defect nucleated by an isolated disk.

Modelação e monitorização do comportamento dinâmico de barragens abóbada: identificação modal automática

Trabalho Final de Mestrado elaborado no Laboratório de Engenharia Civil (LNEC) para obtenção do grau de Mestre em Engenharia Civil pelo Instituto Superior de Engenharia de Lisboa no âmbito do protocolo de cooperação entre o ISEL e o LNEC

Using a cohesive damage model to predict the tensile behaviour of CFRP single-strap repairs

This work addresses both experimental and numerical analyses regarding the tensile behaviour of CFRP single-strap repairs. Two fundamental geometrical parameters were studied: overlap length and patch thickness. The numerical model used ABAQUS® software and a developed cohesive mixed-mode damage model adequate for ductile adhesives, and implemented within interface finite elements. Stress analyses and strength predictions were carried out. Experimental and numerical comparisons were performed on failure modes, failure load and equivalent stiffness of the repair. Good correlation was found between experimental and numerical results, showing that the proposed model can be successfully applied to bonded joints or repairs.

Numerical evaluation of three-dimensional scarf repairs in carbon-epoxy structures

The widespread employment of carbon-epoxy laminates in high responsibility and severely loaded applications introduces an issue regarding their handling after damage. Repair of these structures should be evaluated, instead of their disposal, for cost saving and ecological purposes. Under this perspective, the availability of efficient repair methods is essential to restore the strength of the structure. The development and validation of accurate predictive tools for the repairs behaviour are also extremely important, allowing the reduction of costs and time associated to extensive test programmes. Comparing with strap repairs, scarf repairs have the advantages of a higher efficiency and the absence of aerodynamic disturbance. This work reports on a numerical study of the tensile behaviour of three-dimensional scarf repairs in carbon-epoxy structures, using a ductile adhesive (Araldite® 2015). The finite elements analysis was performed in ABAQUS® and Cohesive Zone Modelling was used for the simulation of damage onset and growth in the adhesive layer. Trapezoidal cohesive laws in each pure mode were used to account for the ductility of the specific adhesive mentioned. A parametric study was performed on the repair width and scarf angle. The use of over-laminating plies covering the repaired region at the outer or both repair surfaces was also tested as an attempt to increase the repairs efficiency. The obtained results allowed the proposal of design principles for repairing composite structures.

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

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

Node adaptation for global collocation with radial basis functions using direct multisearch for multiobjective optimization

Meshless methods are used for their capability of producing excellent solutions without requiring a mesh, avoiding mesh related problems encountered in other numerical methods, such as finite elements. However, node placement is still an open question, specially in strong form collocation meshless methods. The number of used nodes can have a big influence on matrix size and therefore produce ill-conditioned matrices. In order to optimize node position and number, a direct multisearch technique for multiobjective optimization is used to optimize node distribution in the global collocation method using radial basis functions. The optimization method is applied to the bending of isotropic simply supported plates. Using as a starting condition a uniformly distributed grid, results show that the method is capable of reducing the number of nodes in the grid without compromising the accuracy of the solution. (C) 2013 Elsevier Ltd. All rights reserved.

The influence of the boundary conditions on low-velocity impact composite damage

The objective of this work was to study the influence of the boundary conditions on low-velocity impact behaviour of carbon-epoxy composite plates. Experimental work and numerical analysis were performed on [04,904]s laminates. The influence of different boundary conditions on the impacted plates was analysed considering rectangular and square plates. The X-radiography was used as a non-destructive technique to evaluate the internal damage caused by impact loading. A three-dimensional numerical analysis was also performed considering progressive damage modelling. The model includes three-dimensional solid elements and interface finite elements including a cohesive mixed-mode damage model, which allows simulating delamination between different oriented layers. It was verified that plate’s boundary conditions have influence on the delaminated area. Good agreement between experimental and numerical analysis for shape, orientation and size of the delamination was obtained.