Three-dimensional electrical impedance, tomography: A topology optimization approach

Electrical impedance tomography is a technique to estimate the impedance distribution within a domain, based on measurements on its boundary. In other words, given the mathematical model of the domain, its geometry and boundary conditions, a nonlinear inverse problem of estimating the electric impedance distribution can be solved. Several impedance estimation algorithms have been proposed to solve this problem. In this paper, we present a three-dimensional algorithm, based on the topology optimization method, as an alternative. A sequence of linear programming problems, allowing for constraints, is solved utilizing this method. In each iteration, the finite element method provides the electric potential field within the model of the domain. An electrode model is also proposed (thus, increasing the accuracy of the finite element results). The algorithm is tested using numerically simulated data and also experimental data, and absolute resistivity values are obtained. These results, corresponding to phantoms with two different conductive materials, exhibit relatively well-defined boundaries between them, and show that this is a practical and potentially useful technique to be applied to monitor lung aeration, including the possibility of imaging a pneumothorax.

A numerical study of effect of grain boundaries on elastic and plastic properties in nanocomposite materials

Nanocomposite materials have received considerable attention in recent years due to their novel properties. Grain boundaries are considered to play an important role in nanostructured materials. This work focuses on the finite element analysis of the effect of grain boundaries on the overall mechanical properties of aluminium/alumina composites. A grain boundary is incorporated into the commonly used unit cell model to investigate its effect on material properties. By combining the unit cell model with an indentation model, coupled with experimental indentation measurements, the ''effective'' plastic property of the grain boundary is estimated. In addition, the strengthening mechanism is also discussed based on the Estrin-Mecking model.

The subcycled Newmark algorithm

The popular Newmark algorithm, used for implicit direct integration of structural dynamics, is extended by means of a nodal partition to permit use of different timesteps in different regions of a structural model. The algorithm developed has as a special case an explicit-explicit subcycling algorithm previously reported by Belytschko, Yen and Mullen. That algorithm has been shown, in the absence of damping or other energy dissipation, to exhibit instability over narrow timestep ranges that become narrower as the number of degrees of freedom increases, making them unlikely to be encountered in practice. The present algorithm avoids such instabilities in the case of a one to two timestep ratio (two subcycles), achieving unconditional stability in an exponential sense for a linear problem. However, with three or more subcycles, the trapezoidal rule exhibits stability that becomes conditional, falling towards that of the central difference method as the number of subcycles increases. Instabilities over narrow timestep ranges, that become narrower as the model size increases, also appear with three or more subcycles. However by moving the partition between timesteps one row of elements into the region suitable for integration with the larger timestep these the unstable timestep ranges become extremely narrow, even in simple systems with a few degrees of freedom. As well, accuracy is improved. Use of a version of the Newmark algorithm that dissipates high frequencies minimises or eliminates these narrow bands of instability. Viscous damping is also shown to remove these instabilities, at the expense of having more effect on the low frequency response.

The suitability of different FEA models for studying root fractures caused by wedge effect

Finite element analysis (FEA) utilizing models with different levels of complexity are found in the literature to study the tendency to vertical root fracture caused by post intrusion (""wedge effect""). The objective of this investigation was to verify if some simplifications used in bi-dimensional FEA models are acceptable regarding the analysis of stresses caused by wedge effect. Three plane strain (PS) and two axisymmtric (Axi) models were studied. One PS model represented the apical third of the root entirely, in dentin (PS-nG). The other models included gutta-percha in the apical third, and differed regarding dentin-post relationship: bonded (PS-B and Axi-B) or nonbonded (PS-nB and Axi-nB). Mesh discretization and material properties were similar for all cases. Maximum principal stress (sigma(max)) was analyzed as a response to a 165 N longitudinal load. Stress magnitude and orientation varied widely (PS-nG: 10.3 MPa; PS-B: 0.8 MPa; PS-nB: 10.4 MPa; Axi-13: 0.2 MPa, Axi-nB: 10.8 MPa). Axi-nB was the only model where all (sigma(max) vectors at the apical third were perpendicular to the model plane. Therefore, it is adequate to demonstrate the tendency to vertical root fractures caused by wedge effect. Axi-13 showed only part of the (sigma(max) perpendicular to the model plane while PS models showed sigma(max) on the model plane. In these models, sigma(max) orientation did not represent a situation where vertical root fracture would occur due to wedge effect. Adhesion between post and dentin significantly reduced (c) 2007 Wiley Periodicals, Inc.

Earth pressure on cantilever walls at design retained heights

There are many methods for the analysis and design of embedded cantilever retaining walls. They involve various different simplifications of the pressure distribution to allow calculation of the limiting equilibrium retained height and the bending moment when the retained height is less than the limiting equilibrium value, i.e. the serviceability case. Recently, a new method for determining the serviceability earth pressure and bending moment has been proposed. This method makes an assumption defining the point of zero net pressure. This assumption implies that the passive pressure is not fully mobilised immediately below the excavation level. The finite element analyses presented in this paper examine the net pressure distribution on walls in which the retained height is less, than the limiting equilibrium value. The study shows that for all practical walls, the earth pressure distributions on the front and back of the wall are at their limit values, Kp and K-a respectively, when the lumped factor of safety F-r is less than or equal to2.0. A rectilinear net pressure distribution is proposed that is intuitively logical. It produces good predictions of the complete bending moment diagram for walls in the service configuration and the proposed method gives results that have excellent agreement with centrifuge model tests. The study shows that the method for determining the serviceability bending moment suggested by Padfield and Mair(1) in the CIRIA Report 104 gives excellent predictions of the maximum bending moment in practical cantilever walls. It provides the missing data that have been needed to verify and justify the CIRIA 104 method.

Knowledge-based system on optimum design of liquid retaining structures with genetic algorithms

This paper delineates the development of a prototype hybrid knowledge-based system for the optimum design of liquid retaining structures by coupling the blackboard architecture, an expert system shell VISUAL RULE STUDIO and genetic algorithm (GA). Through custom-built interactive graphical user interfaces under a user-friendly environment, the user is directed throughout the design process, which includes preliminary design, load specification, model generation, finite element analysis, code compliance checking, and member sizing optimization. For structural optimization, GA is applied to the minimum cost design of structural systems with discrete reinforced concrete sections. The design of a typical example of the liquid retaining structure is illustrated. The results demonstrate extraordinarily converging speed as near-optimal solutions are acquired after merely exploration of a small portion of the search space. This system can act as a consultant to assist novice designers in the design of liquid retaining structures.

A coupled knowledge-based expert system for design of liquid-retaining structures

This paper describes a coupled knowledge-based system (KBS) for the design of liquid-retaining structures, which can handle both the symbolic knowledge processing based on engineering heuristics in the preliminary synthesis stage and the extensive numerical crunching involved in the detailed analysis stage. The prototype system is developed by employing blackboard architecture and a commercial shell VISUAL RULE STUDIO. Its present scope covers design of three types of liquid-retaining structures, namely, a rectangular shape with one compartment, a rectangular shape with two compartments and a circular shape. Through custom-built interactive graphical user interfaces, the user is directed throughout the design process, which includes preliminary design, load specification, model generation, finite element analysis, code compliance checking and member sizing optimization. It is also integrated with various relational databases that provide the system with sectional properties, moment and shear coefficients and final member details. This system can act as a consultant to assist novice designers in the design of liquid-retaining structures with increase in efficiency and optimization of design output and automated record keeping. The design of a typical example of the liquid-retaining structure is also illustrated. (C) 2003 Elsevier B.V All rights reserved.

On the induced electric field gradients in the human body for magnetic stimulation by gradient coils in MRI

Prior theoretical studies indicate that the negative spatial derivative of the electric field induced by magnetic stimulation may he one of the main factors contributing to depolarization of the nerve fiber. This paper studies this parameter for peripheral nerve stimulation (PNS) induced by time.-varying gradient fields during MRI scans. The numerical calculations are based on an efficient, quasi-static, finite-difference scheme and an anatomically realistic human, full-body model. Whole-body cylindrical and planar gradient sets in MRI systems and various input signals have been explored. The spatial distributions of the induced electric field and their gradients are calculated and attempts are made to correlate these areas with reported experimental stimulation data. The induced electrical field pattern is similar for both the planar coils and cylindrical coils. This study provides some insight into the spatial characteristics of the induced field gradients for PNS in MRI, which may be used to further evaluate the sites where magnetic stimulation is likely to occur and to optimize gradient coil design.

Determinação do momento crítico de flambagem lateral com distorção em vigas mistas contínuas de aço e concreto com perfis de alma senoidal

As vigas mistas de aço e concreto estão sendo largamente utilizadas em construções de edifícios e pontes. Ao se combinar o aço com o concreto obtêm-se estruturas mais econômicas, uma vez que se tira proveito das melhores características de cada material. Nas regiões de momento negativo de uma viga mista contínua, a mesa inferior e parte da alma estão comprimidas, se a alma do perfil não tiver rigidez suficiente para evitar a flexão lateral, ela distorcerá gerando um deslocamento lateral e um giro na mesa comprimida, caracterizando um modo de flambagem denominado flambagem lateral com distorção (FLD). O procedimento de verificação à FLD da EN 1994-1-1:2004 originou o método de cálculo da ABNT NBR 8800:2008, entretanto a EN 1994-1-1:2004 não fornece expressão para o cálculo do momento crítico elástico, enquanto a ABNT NBR 8800:2008 prescreve uma formulação proposta por Roik, Hanswille e Kina (1990) desenvolvida para vigas mistas com perfis de alma plana. Embora as normas prescrevam um método de verificação à FLD para vigas mistas com perfis de alma plana, poucos estudos têm sido feitos sobre esse estado-limite. Além disso, tanto a ABNT NBR 8800:2008 quanto as normas internacionais não abordam perfis de alma senoidal. Neste trabalho, foram implementadas análises de flambagem elástica, com auxílio do software ANSYS 14.0 (2011), em modelos de elementos finitos que retratem o comportamento à FLD de vigas mistas de aço e concreto com perfis de alma plana e senoidal. Os modelos numéricos foram constituídos pelo perfil de aço, por uma mola rotacional que restringe parcialmente o giro da mesa superior e uma restrição ao deslocamento lateral, ao longo de todo o comprimento da viga. Os resultados numéricos são comparados com os obtidos pelas formulações de Roik, Hanswille e Kina (1990) e de Hanswille (2002), adaptadas para levar em consideração a corrugação da alma do perfil de aço. Para avaliação das formulações supracitadas e da consistência da modelagem numérica adotada, o momento crítico elástico foi determinado para vigas mistas com perfis de aço de alma plana. Como resultado, um método para o cálculo do momento crítico elástico de vigas mistas de alma senoidal é proposto.

Análise numérico-experimental do comportamento da aderência aço-concreto

Dissertação (mestrado)—Universidade de Brasília, Faculdade de Tecnologia, Departamento de Engenharia Civil e Ambiental, 2016.

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.

Simulação numérica do comportamento dinâmico de uma laje

Nesta dissertação pretende-se simular o comportamento dinâmico de uma laje de betão armado aplicando o Método de Elementos Finitos através da sua implementação no programa FreeFEM++. Este programa permite-nos a análise do modelo matemático tridimensional da Teoria da Elasticidade Linear, englobando a Equação de Equilíbrio, Equação de Compatibilidade e Relações Constitutivas. Tratando-se de um problema dinâmico é necessário recorrer a métodos numéricos de Integração Directa de modo a obter a resposta em termos de deslocamento ao longo do tempo. Para este trabalho escolhemos o Método de Newmark e o Método de Euler para a discretização temporal, um pela sua popularidade e o outro pela sua simplicidade de implementação. Os resultados obtidos pelo FreeFEM++ são validados através da comparação com resultados adquiridos a partir do SAP2000 e de Soluções Teóricas, quando possível.

Response of a multi-domain continental margin to compression: Study from seismic reflection-refraction and numerical modelling in the Tagus Abyssal Plain

The effects of the Miocene through Present compression in the Tagus Abyssal Plain are mapped using the most up to date available to scientific community multi-channel seismic reflection and refraction data. Correlation of the rift basin fault pattern with the deep crustal structure is presented along seismic line IAM-5. Four structural domains were recognized. In the oceanic realm mild deformation concentrates in Domain I adjacent to the Tore-Madeira Rise. Domain 2 is characterized by the absence of shortening structures, except near the ocean-continent transition (OCT), implying that Miocene deformation did not propagate into the Abyssal Plain, In Domain 3 we distinguish three sub-domains: Sub-domain 3A which coincides with the OCT, Sub-domain 3B which is a highly deformed adjacent continental segment, and Sub-domain 3C. The Miocene tectonic inversion is mainly accommodated in Domain 3 by oceanwards directed thrusting at the ocean-continent transition and continentwards on the continental slope. Domain 4 corresponds to the non-rifted continental margin where only minor extensional and shortening deformation structures are observed. Finite element numerical models address the response of the various domains to the Miocene compression, emphasizing the long-wavelength differential vertical movements and the role of possible rheologic contrasts. The concentration of the Miocene deformation in the transitional zone (TC), which is the addition of Sub-domain 3A and part of 3B, is a result of two main factors: (1) focusing of compression in an already stressed region due to plate curvature and sediment loading; and (2) theological weakening. We estimate that the frictional strength in the TC is reduced in 30% relative to the surrounding regions. A model of compressive deformation propagation by means of horizontal impingement of the middle continental crust rift wedge and horizontal shearing on serpentinized mantle in the oceanic realm is presented. This model is consistent with both the geological interpretation of seismic data and the results of numerical modelling.

Projecto de fundações e estrutura de um edifício destinado a pavilhão gimnodesportivo

Este trabalho tem como objectivo a elaboração do projecto de estruturas de um edifício destinado a pavilhão gimnodesportivo, caracterizando as suas diferentes fases de execução, desde a etapa inicial de concepção até à fase final de dimensionamento. Trata-se de um projecto complexo de uma estrutura com elementos estruturais em betão armado e pré-esforçado, e com muros de contenção. Na concepção do edifício foram utilizados os critérios gerais de dimensionamento presentes na regulamentação Europeia (Eurocódigos), uma vez que estes elementos representam o futuro da regulamentação de estruturas em termos Europeus, vindo substituir a nível nacional o “Regulamento de Segurança e Acções para Estruturas de Betão Armado (RSA)” e o “Regulamento para Estruturas de Betão Armado e Pré- Esforçado (REBAP)”. A adopção das normas europeias representam assim um elevado desafio devido ao aumento da complexidade na concepção e dimensionamento de estruturas que estes regulamentos traduzem, principalmente o Eurocódigo 8, que define de um modo mais detalhado e complexo a análise sísmica, relativamente à regulamentação actual em vigor. Devido à elevada complexidade que os projectos de estruturas apresentam, utilizam-se actualmente ferramentas de cálculo automático. No dimensionamento deste edifício foi utilizado um programa tridimensional de elementos finitos para a modelação da estrutura. Pretende-se com a escolha deste projecto e dos métodos de dimensionamento presentes nos Eurocódigos, o desenvolvimento de um trabalho detalhado e correcto, permitindo assim adquirir conhecimentos importantes relativamente às futuras normas, e pôr em prática as competências e os conhecimentos obtidos ao longo curso.

Modelização de um dispositivo de travagem automático de via, através do método dos elementos finitos

Apresenta-se a modelização, através do Método dos Elementos Finitos, de um sistema de proteção automática de comboios, denominado DTAV – Dispositivo de Travagem Automático de Via. A modelização é realizada com o auxílio do software Ansys 13.0. Descreve-se a evolução histórica dos sistemas de proteção automática de comboios, desde a época vitoriana até à atualidade. São categorizados os vários tipos de sistemas de proteção existentes, consoante as suas funcionalidades e os seus princípios físicos de funcionamento. É apresentada uma breve descrição da normalização em vigor, aplicada aos sistemas de proteção automáticos de comboios. Descreve-se em pormenor o sistema DTAV, nomeadamente a sua funcionalidade, o seu princípio físico de funcionamento e o conjunto de equipamentos de que é constituído. Apresenta-se uma breve introdução sobre o método dos elementos finitos, enquadrando o modelo criado com os princípios físicos em que se fundamenta, nomeadamente através da descrição das equações de eletromagnetismo. Desenvolve-se modelo do sistema DTAV, evidenciado as etapas da construção e os resultados obtidos, validados por um conjunto de especificações e ensaios prévios realizados em laboratório.