Nesta tese desenvolvemos um método numérico para diagonalizar o Hamiltoniano de Heisenberg iterativamente. O método consiste basicamente em diagonalizar cadeias de spins, cada vez maiores, em que cada passo da diagonalização corresponde à adição de um novo spin à cadeia. A base de vetores para calcular o Hamiltoniano de uma cadeia de N spins, HN, é construída por meio do produto direto dos autovetores do Hamiltoniano Hn-1 da rede diagonalizada no passo anterior, pelos autoestados correspondentes ao N-ésimo spin adicionado. Além de usar a comutação do Hamiltoniano com a componente azimutal do spin total, Sz, prática comum em outros métodos, usufruímos da conservação com o quadrado do spin total, S2. Para uma classe específica de redes também implementamos a simetria de reflexão. Obtemos o espectro completo de energia de cadeias de spins 1/2 com até 20 sítios, para as quais mostramos resultados da dependência com a temperatura da susceptibilidade magnética e do calor específico, para redes com impurezas tipo spin substitucionais, com defeitos nas ligações ou com efeitos de bordas, isto é, para sistemas sem invariância translacional. Usualmente essa restrição impõe enormes dificuldades em métodos tradicionais. Para diagonalizar cadeias com um número maior de sítios, implementamos um procedimento que seleciona os estados de mais baixa energia para serem usados na base de vetores do passo seguinte. Com esse tipo de truncamento de estados, fomos capazes de obter o estado fundamental e alguns estados de baixa energia de cadeias com mais de uma centena de sítios, com precisão de até cinco algarismos significativos. Nossos resultados reproduzem os da literatura para os casos conhecidos, em geral sistemas homogêneos. As aproximações desenvolvidas recentemente no contexto da Teoria do Funcional da Densidade, aplicada ao modelo de Heisenberg, e que também se aplicam a sistemas inomogêneos, estão em conformidade com nossos resultados numericamente exatos. Generalizamos o método para diagonalizar escadas de spins 1/2. Calculamos o estado fundamental e o gap de energia desse sistema, onde variamos a razão entre os acoplamentos ao longo das pernas da escada e ao longo dos degraus da mesma; nossos resultados são comparados com os da literatura. Apresentamos também a implementação do método iterativo no modelo de Hubbard, que descreve um sistema de spins itinerantes. Sabe-se que no regime de alta repulsão Coulombiana entre os spins e densidade um (número de spins igual ao número de sítios da cadeia), esse modelo é mapeado no modelo de Heisenberg, resultado que é verificado numericamente em nosso procedimento por meio do cálculo de energias de ambos os modelos em um regime paramétrico apropriado. / In this Thesis we develop a numerical method to diagonalize the Heisenberg model iteratively. In essence, we diagonalize spin chains in steps, each one corresponding to an addition of a spin to a smaller chain. The basis vectors to calculate the Hamiltonian of a N-spin chain, HN, is built by means of the direct product of the eigenvectors of the (N-1)-spin Hamiltonian, diagonalized on the previous step, by the eigenstates of the N-th added spin. Besides the common use of the conservation of the z-component of the total spin, Sz, we also exploit the conservation of the squared total spin, S2. For a specific class of spin systems we also implemented the reflection symmetry. We obtain the entire energy spectrum of spin-1/2 chains up to 20 sites, for which we show the temperature dependence of the magnetic susceptibility and specific heat, for systems with substitutional impurity spins, bond defects, border effects, i.e., for systems without translational invariance. This normally imposes enormous restrictions in many traditional methods. In order to diagonalize chains with a larger number of sites we implemented a procedure that selects lower energy states to be used in the basis vector on the next step. Using this truncation scheme, we are able to obtain low-lying energy states for chains with more than a hundred sites, up to five significant figures of accuracy. Our results reproduce those of the literature for the known cases, in general homogeneous systems. The approaches recently developed in the context of Density Functional Theory to the Heisenberg model, which also apply to inhomogeneous systems, are consistent with our numerical results. We generalize the method to diagonalize spin-1/2 ladders. We calculate the ground-state and the energy gap of this system, for arbitrary ratio of the couplings along the lags or over the rungs of the ladder. We also present the implementation of our iterative method to the Hubbard model, which describes a system of itinerant spins. It is known that in the regime of high Coulomb repulsion between the spins and unitary density (number of spins equal to the number of sites in the chain), this model is mapped onto Heisenberg one, a result which is verified numerically in our procedure by calculating the energy spectrum of both models in na appropriated parametric regime.

Transformadores de potência, como produtos não seriados, estão sujeitos a uma série de possíveis desvios de processo. As divergências que transformadores podem apresentar em relação ao seu projeto original podem ter consequências em sua capacidade de suportar carregamentos mecânicos que lhe serão aplicados. Seu transporte da fábrica às subestações elétricas é um dos principais eventos durante o qual suas estruturas são carregadas. Este trabalho foca em descrever e analisar um transformador de potência genérico em um modelo bidimensional, sob o ponto de vista da Mecânica Analítica e da análise de variabilidades com o objetivo de identificar parâmetros importantes à integridade estrutural do equipamento, especialmente aqueles relacionados ao modo de falha de destacamento das chapas do núcleo nos step laps. O modelo numérico do transformador consiste na representação de corpos rígidos conectados entre si por elementos de rigidez e amortecimento, tendo o conjunto do núcleo ferromagnético particularmente mais detalhado de modo a descrever o atrito existente em suas junções. A avaliação dos esforços nesses locais é feita por meio do método dos multiplicadores de Lagrange, ao considerar o atrito como uma restrição ao movimento, sendo que as forças e momentos limites para que ocorra o deslizamento planar das chapas do núcleo são determinados utilizando-se Newton-Raphson ou outras técnicas de otimização. / Power transformers, as non serialized products, are subject to a number of process deviations. The divergences that transformers present from their original design may have consequences on their ability to withstand mechanical loads which will be applied on them. Their transportation from workshop to power plants is one of the main events their structures are loaded. This work focuses on describing and analyzing a generic power transformer in a two-dimensional model, under the sight of Analytical Mechanics and Structural Reliability in order to identify important parameters related to their structural integrity, especially those concerning the failure mode of dismantling of the ferromagnetic core in the step-lap joints. The numerical model of the transformer consists of a representation of rigid bodies, connected to each other by stiness and damping elements, having the assembly of the ferromagnetic core more detailed in order to describe the friction present on their step-lap joint. The evaluation of friction in these locations is made using the Lagrange Multipliers Method, by considering it as a movement restriction, and the limit values of forces and moments to which the planar sliding of the core sheets occurs, damaging the behavior of the electric machine, are determined using Newton-Raphson and Optimization Methods.

Wave-based models of sound propagation can be used to predict and synthesize sounds as they would be heard naturally in room acoustic environments. The numerical simulation of such models with traditional time-stepping grid-based methods can be an expensive process, due to the sheer size of listening environments (e.g., auditoriums and concert halls) and due to the temporal resolution required by audio rates that resolve frequencies up to the limit of human hearing. Finite difference methods comprise a simple starting point for such simulations, but they are known to suffer from approximation errors that may necessitate expensive grid refinements in order to achieve sufficient levels of accuracy. As such, a significant amount of research has gone into designing finite difference methods that are highly accurate while remaining computationally efficient. The problem of designing and using accurate finite difference schemes is compounded by the fact that room acoustics models require complex boundary conditions to model frequency-dependent wall impedances over non-trivial geometries. The implementation of such boundary conditions in a numerically stable manner has been a challenge for some time. Stable boundary conditions for finite difference room acoustics simulations have been formulated in the past, but generally they have only been useful in modelling trivial geometries (e.g., idealised shoebox halls). Finite volume methods have recently been shown to be a viable solution to the problem of complex boundary conditions over non-trivial geometries, and they also allow for the use of energy methods for numerical stability analyses. Finite volume methods lend themselves naturally to fully unstructured grids and they can simplify to the types of grids typically used in finite difference methods. This allows for room acoustics simulation models that balance the simplicity of finite difference methods for wave propagation in air with the detail of finite volume methods for the modelling of complex boundaries. This thesis is an exploration of these two distinct, yet related, approaches to wave-based room acoustic simulations. The overarching theme in this investigation is the balance between accuracy, computational efficiency, and numerical stability. Higher-order and optimised schemes in two and three spatial dimensions are derived and compared, towards the goal of finding accurate and efficient finite difference schemes. Numerical stability is analysed using frequency-domain analyses, as well as energy techniques whenever possible, allowing for stable and frequency-dependent boundary conditions appropriate for room acoustics modelling. Along the way, the use of non-Cartesian grids is investigated, geometric relationships between certain finite difference and finite volume schemes are explored, and some problems associated to staircasing effects at boundaries are considered. Also, models of sound absorption in air are incorporated into these numerical schemes, using physical parameters that are appropriate for room acoustic scenarios.

Neste trabalho um método de partículas para simular a dinâmica de sólidos elásticos e interação fluido estrutura e implementado. O método e baseado no Moving Particle Semi-implicit (MPS), originalmente desenvolvido para escoamentos incompressíveis com superfície livre. A estratégia principal do MPS e substituir os operadores diferenciais das equações governantes por operadores diferenciais discretos em uma distribuição de nos irregulares, derivados de um modelo de interação entre partículas. Inicialmente são apresentados os detalhes da formulação do método e modelos constitutivos utilizados. Uma condição simplificada de fragmentação e proposta, assim como um algoritmo de detecção de contato, permitindo a fragmentação entre vários sólidos. No caso da interação fluido-estrutura, as partículas de superfície do solido são tratadas como partículas de fluido e as pressões destas partículas são calculadas pela resolução da equação de Poisson para a pressão, tal como as partículas de fluido. Desta forma, o acoplamento entre solido e fluido e realizado utilizando o deslocamento e velocidade do solido elástico, como condições de contorno do fluido, e a pressão na interface, obtida pela resolução do movimento do fluido, e aplicada ao movimento do solido elástico. São apresentados e detalhados os algoritmos de solido elástico, fragmentação, colisão e acoplamento fluido-estrutura. Validações qualitativas e quantitativas do método são realizadas para casos estáticos e dinâmicos sujeitos a diferentes condições de contorno, comparando os resultados numéricos obtidos pelo MPS, outros métodos numéricos, soluções analíticas e medições experimentais presentes na literatura. / In this work a particle method to simulate the dynamics of elastic solids and fluid-structure interaction is implemented. It is based on the Moving Particle Semi-implicit Method (MPS), which was originally developed for incompressible flows with free surface. The main strategy of the MPS is to replace the differential operators of the governing equations by discrete differential operators on irregular nodes, which are derived from a model of interaction between particles. Initially details of the method and constitutive equations are shown. A simplified condition of fragmentation and collision between solids are proposed to allow the investigation of fragmentation amount multiple solids. In case of fluid-structure interaction, the solid\'s surface particles are treated as a fluid particle and the pressures of the surface particles are computed by solving Poisson equation for the pressure, just as the fluid particles. Therefore, the coupling between solid and fluid is done by using the displacement and velocity of elastic solid as the boundary conditions of the fluid, and the pressure at the interface, which is obtained when solving the fluid motion, is used to calculate the motion of the elastic solid. The algorithms for elastic solid, fragmentation, collision and fluid-structure interaction are presented and detailed. The qualitative and quantitative validations of the method are carried out herein considering static and dynamic cases subjected to deferent boundary conditions by comparing the numerical results from MPS with other numerical, analytical and experimental results available in the literature.

Esta dissertação discute o desenvolvimento de modelos analíticos e numéricos para as características elétricas de transistores de fios quânticos. Sendo assim, realizou-se um estudo implementando uma sequência de formalismos e ferramentas computacionais para solução auto-consistente das equações de Schrödinger e Poisson para poços e fios quânticos. Com a utilização deste método numérico pode-se determinar os auto-estados os níveis de energias e as densidades eletrônicas de portadores livres, dentre outros parâmetros relevantes para dispositivos de fio quântico. Adicionalmente, realizou-se um estudo analítico das heteroestruturas semicondutoras de interesse para a área de dispositivos de dimensionalidade reduzida. Este estudo levou a obtenção de resultados referentes ao desenvolvimento de modelos teóricos para as características elétricas de dispositivos baseados no mecanismo de tunelamento ressonante. Os resultados obtidos para a característica corrente-tensão (I-V) nas heteroestruturas investigadas foram contrastados satisfatoriamente com os encontrados na literatura. Este ferramental analítico foi então aplicado para computar o coeficiente de transmissão eletrônico de um diodo de fio quântico com tunelamento ressonante. / This dissertation discusses the development of analytical and numerical models for the electrical characteristics of quantum wire transistors. A study is carried out, implementing a sequence of formalisms and computational tools for the self-consistent solution of the equations of Schrödinger and Poisson in quantum wells and quantum wires. By using this numerical formulation it is possible to determine the eigenstates, energy levels and free-carrier electronic density, among other relevant parameters for quantum wire devices. In addition, we also conducted an analytical study concerning semiconductor heetrostrucures of interest for reduced dimensionality devices applications. This study led to results regarding the development of theoretical models for the electrical characteristics of devices based on the resonant tunneling mechanism. The results obtained for the current-voltage (I-V) characteristics in the investigated heterostructures were satisfactorily compared to those available at the published literature and this analytical tool was then used to compute the electronic transmission coefficient in a resonant tunneling quantum wire diode.

O Tanque de Provas Numérico (TPN) é um laboratório pioneiro em hidrodinâmica aplicada e fruto de uma colaboração entre a indústria brasileira de petróleo (PETROBRAS S.A.) e as principais instituições de pesquisa do país. Seu principal objetivo é atuar como parceiro da indústria offshore e de petróleo, colaborando para a obtenção da autossuficiência da produção nacional de petróleo como uma poderosa ferramenta para projeto e análise de sistemas flutuantes de produção de óleo e gás. O coração do TPN é um cluster de computadores SMP que é hoje um dos maiores agrupamentos do Brasil para fins de pesquisa. Um dos focos de atenção do TPN é a aplicação do Moving Particles Semi-implicit Method (MPS) na exploração de soluções para muitos problemas de Engenharia. Por trabalhar sem a necessidade do uso de malhas (método tradicional Euleriano), o método tem diversas aplicações na simulação de corpos flutuantes e na hidrodinâmica aplicada, sendo atualmente utilizado para realização de estudos sobre a influência do movimento de ondas em navios; simulações de fenômenos que envolvem fragmentações; superfícies livres; grandes deformações; dinâmica de fluidos em condições extremas, como é o caso em processos prospecção do petróleo onde muitas vezes é difícil e economicamente inviável fazer ensaios físicos. Devido ao grande número de partículas utilizadas na simulação de sistemas complexos pelo método MPS, é necessário aproveitar de forma eficiente os recursos computacionais disponíveis para a análise de modelos com o refinamento adequado às aplicações práticas. Com tera-FLOPS disponíveis na rede cluster do TPN para modelagem computacional, há uma grande necessidade de uma solução computacional paralela altamente escalonável que, além disto, seja fácil de manutenção e extensibilidade. Dentro desta linha de pesquisa, foi desenvolvida uma solução com essas características através do emprego de modernas técnicas de engenharia de software. / The Numerical Offshore Tank (TPN) is a pioneer laboratory in applied hydrodynamics and result of collaboration between the Brazilian oil (Petrobras SA) and the major research institutions in the country. Its main purpose is to act as a partner of industry and offshore oil, contributing to the achievement of self-sufficiency of domestic oil production as a powerful tool for design and analysis of floating production systems for oil and gas. The heart of TPN is a cluster of SMP computers that is now one of the largest groupings of Brazil for research purposes. One focus of attention of TPN is the application of Moving Particles Semi-implicit Method (MPS) in exploring solutions to many engineering problems. By working without the use of mesh (Eulerian traditional method), the method has several applications in the simulation of floating bodies and applied hydrodynamics, currently being used for studies on the influence of the movement of ships in waves; simulations of phenomena involving fragmentation; free surfaces, large deformations; fluid dynamics in extreme conditions, as is the case in processes where petroleum exploration is often difficult and uneconomical to do physical tests. Due to the high number of particles used in the simulation of complex systems by the MPS method, it is necessary to efficiently take advantage of the computational resources available for the analysis of models with the refinement suitable for practical applications. With tera-FLOPS available in the TPN network cluster for computational modeling, there is a great need for a parallel highly scalable solution which, moreover, must be easy maintenance and extensibility. Within this line of research, we developed a solution with these characteristics through the use of modern software engineering techniques.

Esta dissertação apresenta uma revisão histórica dos trabalhos acerca de superfícies minimais, ressaltando a pertinência da analogia entre a busca de superfícies de mínima área e a busca de formas de membranas estruturais sujeitas a um estado de tensões superficiais, homogêneo e isótropo. São colocados alguns conceitos geométricos das superfícies parametrizáveis, com base na geometria diferencial, a fim de realizar o equilíbrio diferencial de membranas e determinar as suas equações de equilíbrio. Além disso, é apresentada uma metodologia puramente geométrica para a determinação de superfícies minimais, baseada na minimização do funcional da área, dado pela soma das áreas das facetas triangulares nas quais a superfície é discretizada. O trabalho discute a formulação matemática do problema e apresenta resultados obtidos tanto por meio das rotinas implementadas no software MATLAB quanto por meio daquelas da biblioteca de otimização deste mesmo software. Finalmente, são realizados alguns exemplos e um teste de convergência, comparando as superfícies resultantes dos métodos numéricos com suas respectivas respostas analíticas. A geometria final de um dos exemplos é verificada por meio da analogia dos filmes de sabão, realizando-se uma análise não-linear de equilíbrio através do software Ansys. As soluções foram bastante satisfatórias, resultando em formas muito próximas das analíticas e com pequenos erros relativos das áreas. O teste de convergência também comprovou que o refinamento da discretização leva a uma solução mais próxima da desejada. Portanto, os procedimentos apresentados podem ser empregados no processo de busca da forma de membranas estruturais. / This dissertation presents a historical review on the theoretical developments on minimal surfaces, highlighting the important analogy between the problems of finding minimal area surfaces and finding membrane surfaces with homogeneous and isotropic stress fields. Some geometric concepts of the parametric surfaces are placed, on the basis of differential geometry, in order to do the differential equilibrium of membranes and to achieve its equilibrium equations. Moreover, a purely geometric methodology for the determination of minimal surfaces is presented, based on the minimization of the area functional, which is computed by the simple addition of a finite number of triangular facet areas in which the surface is divided. It discusses the mathematical formulation of the problem as well as some results obtained with the algorithms implemented in MATLAB and others obtained with the aid of MATLAB optimization routines. Finally, some examples and a convergence test are produced, comparing their analytical and numerical results. The final geometry of one of examples is verified by means of the soap film analogy, with a nonlinear equilibrium analysis through Ansys. The solutions have been sufficiently satisfactory, resulting forms very close to the analytical ones and with small areas relative errors. Convergence test also confirm that the method lead to numerical solutions as close to the analytical one as required, as long as the triangular facets mesh is refined. Therefore, the presented procedures can be used in structural membranes form finding.

Este trabalho estuda a técnica de reconstrução de imagens conhecido como tomografia por impedância elétrica em um domínio bidimensional. Esta técnica consiste na alocação de eletrodos na fronteira do volume e uma fonte injeta padrões de corrente através dos eletrodos e medem-se as voltagens resultantes na fronteira. Com estes dados estima-se a condutividade (ou resistividade) do interior do domínio criando-se uma imagem do mesmo. A tomografia por impedância elétrica é um problema inverso e mal posto no sentido de Hadamard. Estudam-se diversos métodos de solução para resolver o problema direto usando métodos numéricos como diferenças finitas e volumes finitos. Proporemos os métodos numéricos a serem aplicados na solução do problema direto os quais serão testados com problemas onde a solução analítica é conhecida. Posteriormente aplicaremos os métodos propostos ao problema especifico. Uma questão importante na reconstrução de imagens é propor a maneira como aproximar o Jacobiano (ou matriz de sensibilidade) do problema, assim desenvolvemos uma técnica para a aproximação do mesmo usando os dados fornecidos pelo problema direto. / In this work is studied the technique of reconstruction of images known as electrical impedance tomography for a two-dimensional domain. This technique consists in the allocation of electrodes on the border of the volume and a source injects patterns of current through the electrodes and then measuring voltages through the other electrodes. With these data it is estimated the conductivity (or resistivity) on the interior of the domain and an image is create of it. The electrical impedance tomography is an inverse and ill conditioned problem in the Hadamard sense. In this work, is studying some numerical methods to solve the direct problem and are applied numerical methods such as the finite difference method and the finite volume method. It is proposed some numerical methods to solve the direct problem which will be tested with analytical problems where the solution is known. Later, apply the methods proposed to the specific issue. An important issue in the reconstruction problems is about the Jacobian (or sensitivity matrix) aproximation, thus proposing a technique for the calculation of even using the data provided by the direct problem. Keywords:

Este trabalho objetivou a ampliação dos conhecimentos sobre a dinâmica de chamas de difusão turbulenta. A abordagem escolhida foi o método numérico da simulação de grandes escalas turbulentas (LES) acoplado ao modelo de combustão do tipo folha de chama ou \"Flamesheet\". A modelagem dos efeitos de dissipação das pequenas escalas turbulentas, foi realizada pelo modelo de Smagorinsky (1963). Foi feita a verificação e validação do código implementado pelo método das soluções manufaturadas, descrito por Roache (2002), e pelas soluções de engenharia do tipo tampa móvel em cavidade laminar, escoamento de Poiseuille e jato laminar. Neste trabalho foi estudada a chama de difusão turbulenta do tipo D, padronizada pelo \"SANDIA National Laboratories\". Foi verificado que, quando estudado o jato isotérmico, este obedece ao critério de jato auto-sustentável, definido por Hussein et al. (1994). Já os resultados obtidos para o jato reativo, quando comparados aos dados experimentais obtidos por Masri et al. (1996) e Barlow e Frank (1998), apresentaram boa concordância para os campos de velocidades e razoável correspondência para os campos de fração de mistura e temperatura, indicando que maiores estudos são necessários. / The mean goal of this work was to better understand the dynamics of turbulent diffusion flames. The turbulent flow was resolved appling the methodology of Large Eddy Simulation, coupled with the Flamesheet model for reactive systems. The effects of the turbulent subgrid dissipation were accounted with the Smagorinsky (1963) subgrid model. The numerical code verification and validation were accomplished with the manufactured solutions method, described by Roache (2002), and with the method of engineering solutions of the type of moveable wall in laminar cavity flow, Poiseuille flow and laminar jet flow. In this work was studied the turbulent flame of the type D flow, standardized by the SANDIA National Laboratories. It was verified that, when the isotermal jet flow was studied, the criterium of a self-sustained jet, defined by Hussein et al. (1994), was achieved. The results for the reactive flow, when compared against experimental data obtained by Masri et al. (1996) e Barlow e Frank (1998), presented resonable agreement for the velocity field and satisfactory correspondence for mixture fraction and temperature fields, indicating that further studies are necessary.

A utilização de perfis formados a frio na construção civil tem sido motivada pela elevada eficiência estrutural, expressa pela relação entre capacidade resistente e peso, e pela facilidade de fabricação, caracterizada pela possibilidade de produção de elementos com diferentes seções transversais. Devido à alta esbeltez dos elementos que os constituem, o projeto desses perfis, à temperatura ambiente ou em situação de incêndio, é governado pelos fenômenos de instabilidade local, distorcional e global. O objetivo dessa Tese é o desenvolvimento de ferramentas computacionais que possibilitem a avaliação do comportamento estrutural de pilares de aço formados a frio em situação de incêndio. Para esse fim são desenvolvidos dois programas de computador. O primeiro programa, denominado de ATERM, tem o objetivo de determinar o campo de temperaturas em regime transiente, de estruturas bidimensionais formadas por qualquer material, submetidas a qualquer curva temperatura-tempo de incêndio, com base no método dos elementos finitos. A esse programa acopla-se o programa ATERM-DIM para o dimensionamento em regime plástico de vigas de aço continuamente travadas, em situação de incêndio. Resultados obtidos do ATERM são comparados aos obtidos do programa sueco de computador Super Tempcalc. O segundo programa, denominado de INSTAB, realiza análises lineares e não lineares de estabilidade de perfis formados a frio, considerando as instabilidades local, global e distorcional, empregando o método das faixas finitas splines, para material elastofrágil, considerando a redução das propriedades mecânicas devido à temperatura. Valores de esforços resistentes obtidos pelo programa INSTAB são comparados aos determinados por meio do programa comercial de elementos finitos ANSYS, o qual considera o comportamento elastoplástico do material, e a resultados obtidos por meio de um método simplificado de dimensionamento de perfis formados a frio em situação de incêndio, proposto pelo autor para fins de normatização brasileira. Essa comparação permite avaliar o efeito da elastoplasticidade nos perfis axialmente comprimidos. Pretende-se, com base neste estudo, fornecer subsídios para o desenvolvimento de procedimentos para a verificação estrutural de perfis formados a frio em situação de incêndio. / The use of cold-formed steel profile in construction has increased because of its high efficiency, expressed as the ratio between load capacity and weight, and ease of manufacturing, characterized by the possibility of production of elements with different cross sections. Due the high slenderness ratio of the sections elements, the design of these profiles, either at room temperature or in case of fire, is determined by local, distortional and global buckling phenomena. The aim of this Thesis is to develop computational tools that allow the assessment of the structural behavior of cold-formed steel columns in fire. For this purpose, two softwares are developed. The first one, called ATERM, allows determining the temperature field under transient analysis of two-dimensional structures formed by any material, subjected to any time-temperature fire curve, and is based on the finite element method. This software interacts with another program, ATERM-DIM, used for plastic design of lateral restrained beams steel in fire. Results of ATERM are compared to those obtained from the Swedish software Super Tempcalc. The second program, named INSTAB, can perform linear and nonlinear stability studies of cold formed profiles, taking into account the local, distortional and global buckling, by means of the splines finite strip method for elastofragile material, considering the reduction of mechanical properties caused by the increase in temperature. Resistant values obtained by the INSTAB software are compared with results from the commercial finite element program ANSYS, which considers the plastic behavior of the material and also with results obtained by means of a simplified method for design of cold-formed profiles in fire, proposed by the author for the Brazilian fire standard. This comparison allows analyzing the effect of elastoplasticity in columns of cold-formed steel. Another objective of this study is to provide background for the development of procedures for the structural analysis of cold formed profiles in fire.