A Tomografia por Impedância elétrica (TIE) tem como objetivo gerar imagens da distribuição de resistividade dentro de um domínio. A TIE injeta correntes em eletrodos alocados na fronteira do domínio e mede potenciais elétricos através dos mesmos eletrodos. A TIE é considerada um problema inverso, não-linear e mal posto. Atualmente, para gerar uma solução do problema inverso, existem duas classes de algoritmos para estimar a distribuição de resistividade no interior do domínio, os que estimam variações da distribuição de resistividade do domínio e os absolutos, que estimam a distribuição de resistividade. Variações da distribuição de resistividade são o resultado da solução de um sistema linear do tipo Ax = b. O objetivo do presente trabalho é avaliar o desempenho da Programação Linear (PL) na solução do sistema linear, avaliar o algoritmo quanto a propaga- ção de erros numéricos e avaliar os efeitos de restringir o espaço solução através de restrições de PL. Os efeitos do uso de Programação Linear é avaliado tanto em métodos que geram imagens de diferenças, como o Matriz de Sensibilidade, como em métodos absolutos, como o Gauss-Newton. Mostra-se neste trabalho que o uso da PL diminui o erro numérico propagado quando comparado ao uso do algoritmo LU Decomposition. Resulta também que reduzir o espaço solução, diretamente através de restrições de PL, melhora a resolução em resistividade e a resolução espacial da imagem quando comparado com o uso de LU Decomposition. / Electrical impedance tomography (EIT) generates images of the resistivity distribution of a domain. The EIT method inject currents through electrodes placed on the boundary of the domain and measures electric potentials through the same electrodes. EIT is considered an inverse problem, non-linear and ill-conditioned. There are two classes of algorithms to estimate the resistivity distribution inside the domain, difference images algorithms, which estimate resistivity distribution variations, and absolute images algorithms, which estimate the resistivity distribution. Resistivity distribution variations are the solution of a linear system, say Ax = b. In this work, the main objective is to evaluate the performance of Linear Programming (LP) solving an EIT linear system from the point of view of the numerical error propagation and the ability to constrain the solution space. The impact of using LP to solve an EIT linear system is evaluated on a difference image algorithm and on an absolute algorithm. This work shows that the use of LP diminishes the numerical error propagation compared to LU Decomposition. It is also shown that constraining the solution space through LP improves the resistivity resolution and the spatial resolution of the images when compared to LU Decomposition.

Risers em configuração de catenária podem apresentar nível de tração muito baixo próximos ao fundo do mar. Isso pode ocorrer em algumas condições de lançamento, em que sua estrutura pode se configurar de forma quase vertical. Quando se trata de tubos flexíveis ou cabos umbilicais, a composição interna do riser contém armaduras helicoidais de tração. Essas podem induzir a ocorrência de giro axial quando o tubo é solicitado à tração. Se esse movimento não for permitido, surgirá um momento de torção na estrutura. O baixo nível de tração da configuração de catenária combinado com o momento de torção surgido durante o lançamento do riser pode levar a uma forma de instabilidade estrutural que culmina na formação de um laço. Isso é indesejável uma vez que, se existe o laço, dependendo dos esforços submetidos à linha, é possível que o laço se transforme em uma dobra, danificando a estrutura. O presente trabalho analisa as condições de formação de laços em configurações de catenária. Para isso, foram utilizados critérios de estabilidade aplicados a um modelo de elementos finitos, que leva em conta as não linearidades geométrica e de contato entre o riser e o solo. Foi utilizada uma formulação cinematicamente exata de elemento de viga através de uma descrição lagrangiana atualizada, que pode tratar de forma correta as grandes rotações que são impostas ao riser para induzir o surgimento do momento de torção. É mostrado que uma expressão analítica baseada na Fórmula de Greenhill pode prever o fenômeno com boa concordância com os resultados numéricos, mesmo considerando-se fenômenos como contato unilateral com atrito e correntezas marítimas. Além disso, foi feita uma análise paramétrica para prever a formação do laço para diversas geometrias de catenária, procurando generalizar as conclusões obtidas. / Catenary risers can show a very low tension level close to the seabed. This can occur in some launching conditions, in which the structure can be almost vertical. When dealing with flexible pipes or umbilical cables, their internal composition includes the helical tension armors. These can induce an axial rotation to the riser when it is tensioned. If this movement is constrained, a torsion moment will act upon the structure. The low tension level present in the catenary configuration combined with the torsion moment that occurs during the riser launching can lead to a structural instability, resulting in a loop formation. This is undesirable since, if the loop is present, depending on the loads applied to the riser it is possible that it turns into a kink, causing damage. The present work analyzes the loop formation conditions in catenary riser. For that, stability criteria were applied to a finite element model that takes into account the geometric nonlinearities and also the contact between the riser and the seabed. A kinematically-exact beam formulation was adopted, dealing correctly with the large rotations that are imposed to the riser, in order to induce the torsion moment. It is showed that an analytical expression based on Greenhills formula can predict the phenomenon with a good agreement with numerical results, even including phenomena as frictional unilateral contact and sea currents. Furthermore, a parametric analysis was done to predict the loop formation to some catenary geometries, in order to generalize the obtained conclusions.

Estruturas que estão sujeitas a altas temperaturas absolutas, à convecção natural, ou ainda, estruturas que trocam calor na ausência de um meio físico, apresentam relevante transferência de calor por radiação térmica. Este fenômeno é importante para diversas aplicações e processos, como, por exemplo, no funcionamento de coletores solares, satélites, fornos industriais, motores a combustão e usinas nucleares. O presente trabalho de mestrado apresenta a aplicação do método da otimização topológica (MOT) no projeto de estruturas que trocam calor substancialmente por radiação térmica no interior de confinamentos, através da distribuição de material refletor ou de aquecedores. Por meio do MOT, cuja principal característica é a liberdade de distribuição do material dentro de um domínio inicial, é possível adicionar ou remover material de uma determinada região do domínio, criando ou desfazendo fronteiras, de forma livre, visando à obtenção de um projeto otimizado. O algoritmo de otimização é baseado no Método das Assíntotas Móveis (MMA) e é complementado pelo Método dos Elementos Finitos (MEF), para a análise do fenômeno de radiação em confinamentos. Ambos são implementados através do software Matlab. Os casos considerados são o da distribuição de material refletor de radiação térmica ou de aquecedores, sujeitos a uma eventual restrição nas quantidades destes materiais, sobre uma superfície plana, de forma a extremizar-se a irradiação ou a minimizar-se a temperatura em determinada área específica do domínio de projeto. Este problema depende, dentre outros fatores, da geometria da superfície e dos ângulos dos raios incidentes sobre ela. / Structures subjected to high absolute temperatures or to natural convection, as well structures that exchange heat in the absence of a physical medium present significant heat transfer through thermal radiation. This phenomenon is important for several applications and processes, such as in the operation of solar collectors, satellites, industrial furnaces, combustion engines and nuclear plants. The present work shows the application of topology optimization to the design of structures that exchange heat substantially by thermal radiation within an enclosure, through the distribution of reflective material or heaters. However, the design of such radiant enclosures is not trivial and it is necessary to use robust and systematic design tools, such as optimization techniques. Topology optimization is a numerical method which allows finding the layout, or topology, of a structure such that a prescribed objective is maximized or minimized subjected to design constraints. The optimization algorithm, based on the method of moving asymptotes (MMA), and the finite element method for analysis of the phenomenon of radiation in enclosures, are implemented using $Matlab^\\circledR$. The cases considered are the distribution of thermal radiation reflective material or heaters, subjected to a volume fraction constraint of these materials on a flat surface, in order to extremize the irradiation or to minimize the temperature in a specified region of the design domain. This problem depends, among other factors, on the geometry of the surfaces that exchange heat through thermal radiation.

Flexible pipes are used in the offshore oil production to transport fluid and gas from the sea bead to the floating stations, and vice versa. These pipes have several concentric layers, of different materials, geometries and structural functions, since they are exposed to adverse operating environments, subjected to high internal and external pressures, high axial stresses and a series of dynamic loads. The local analysis is an important stage of a flexible pipe design and it consists on determining the stresses and strains distributions along the layers of the pipe. Multipurpose finite element packages are commonly used in the local analysis of flexible pipes, but they possess many limitations due to its generic nature, varying from the absence of specific tools for model creation to heavy restrictions of the number of degrees-of-freedom to make computational processing feasible. At the Polytechnic School of the University of São Paulo, within a research line in progress, several finite macroelements were formulated specifically for structural analysis of flexible pipes, taking into account their particularities, such as geometric patterns and layers assemblage. However, the numerical tools that implement these elements present very high memory and processing consumptions, limiting its usage for large-scale models. Therefore, this work has been motivated by memory and processing limitations of finite element structural analysis of flexible pipes for offshore applications. In this context, the Element-by-Element method, which does not require the global stiffness matrix, was chosen for its potential in memory reduction and processing capabilities, given its scalability and ease of parallelization. After an extensive literature review on numerical methods regarding the EBE method, it was chosen the Element-by-Element Diagonal Preconditioned Conjugate Gradient Method (EBE-PCG) algorithm. Aiming higher computational performance, the finite macroelements formulated by (PROVASI, 2013) were converted to the C++ language, implemented and parallelized in a new analysis tool, named as PipeFEM. The diagonal preconditioned EBE-PCG algorithm was implemented and parallelized with OpenMP. The scalability of the PCG algorithm is directly influenced by the efficiency of the matrix-vector product, an operation that, in the element-by-element method, is computed in a local basis with the blocks that comprise the model, and that requires synchronization techniques when performed in parallel. Four different synchronization strategies were developed, being the one based on geometric- and mesh- based mappings the most efficient of them. Numerical experiments showed a reduction of almost 92% in the EBE-PCG solution time of the parallelized version in comparison to the sequential one. In order to compare the efficiency of PipeFEM with the well-established finite element package ANSYS, a simplified flexible pipe was modeled in both software. PipeFEM was approximately 82 times faster than ANSYS to solve the problem, spending 24.27 seconds against 33 minutes and 18 seconds. In addition to this, PipeFEM required much less memory, 61.8MB against 6.8GB in ANSYS. In comparison to the dense version of MacroFEM, a reduction of more than three orders of magnitude was achieved in memory consumption. Despite the low the rate of convergence presented by the diagonal preconditioner, the implementation is very efficient in computational terms. Therefore, the objectives of this work were fulfilled with the development and application of the EBE method, allowing a reduction of memory and simulation costs. / Tubos flexíveis são utilizados na produção offshore de petróleo para o transporte de fluidos e gás natural das estruturas submersas até as estações flutuantes, e vice-versa. Estes tubos possuem diversas camadas concêntricas, de diferentes materiais, geometrias e funções estruturais, pois são expostos a ambientes adversos de operação, nos quais são submetidos à elevadas pressões internas e externas, elevados carregamentos e tensões axiais, além de uma série de carregamentos dinâmicos. A análise local é uma etapa importante do dimensionamento de um tubo flexível e consiste em determinar as distribuições de tensões e deformações ao longo das camadas do tubo. Pacotes multiuso de elementos finitos são comumente utilizados na análise local de tubos flexíveis, mas, devido as suas naturezas genéricas, possuem limitações que variam desde a ausência de ferramentas específicas para a criação de modelos até restrições pesadas no número total de graus de liberdade para tornar exequível o processo computacional. Na Escola Politécnica da Universidade de São Paulo, dentro de uma linha de pesquisa em andamento, diversos macroelementos finitos foram formulados especificamente para a análise estrutural de tubos flexíveis, levando em consideração as suas particularidades, como por exemplo padrões de geometrias e de montagem de camadas. Entretanto, a ferramenta numérica que implementa esses elementos apresenta elevado consumo de memória e de processamento, o que limita o seu uso para modelos de grande escala. Portanto, este trabalho foi motivado por limitações de memória e processamento em análises estruturais com o método dos elementos finitos para tubos flexíveis de aplicações offshore. Neste contexto, o método elemento-a-elemento, caracterizado pela eliminação da matriz global de rigidez, foi escolhido devido ao seu potencial de redução de consumo de memória e às suas capacidades de processamento, dada a sua escalabilidade e facilidade de paralelização. Após uma extensa revisão bibliográfica em métodos numéricos a respeito do método EBE, foi escolhido a versão diagonalmente precondicionada do método do gradiente conjugado (EBE-PCG). Com o intuito de se obter maior performance computacional, os macroelementos finitos formulados por (PROVASI, 2013) foram convertidos para a linguagem C++, paralelizados e implementado em uma nova ferramenta de análise chamada de PipeFEM, totalmente escrita em C++ e que explora paralelismo em todos as etapas. O algoritmo EBE-PCG foi implementado e paralelizado com OpenMP. A escalabilidade do algoritmo PCG é diretamente influenciada pela eficiência do produto entre matriz e vetor, uma operação que no método elemento-a-elemento é calculada na base local com os blocos que compõem o modelo, o que requer técnicas de sincronização quando realizada de modo paralelo. Quatro diferentes estratégias de sincronização foram desenvolvidas, sendo a mais eficiente delas a que utilizada mapeamentos baseados em características da geometria e malha. Experimentos numéricos mostraram uma redução de quase 92% no tempo de simulação do algoritmo PCG da versão paralelizada em relação à sequencial. De modo a comparar a eficiência do PipeFEM com o pacote bem estabelecido de elementos finitos, ANSYS, um tubo simplificado foi modelado em ambos os programas. PipeFEM foi aproximadamente 82 vezes mais rápido do que o ANSYS, gastando 24.27 segundos contra 33 minutos e 18 segundos. Além disso, PipeFEM consumiu muito menos memória, 61.8MB contra 6.8GB in ANSYS. Em comparação com a versão densa do MacroFEM, uma redução superior a três ordens de grandeza no consum e de memória foi obtida. Assim, apesar da baixa taxa de convergência apresentada pelo precondicionador diagonal, a implementação está muito eficiente em termos computacionais. Portanto, os objetivos deste trabalho foram alcançados com o desenvolvimento e aplicação do método EBE, o que permitiu uma redução considerável dos custos de simulação e memória.

Este trabalho é uma contribuição importante para o desenvolvimento de um modelo numérico de arênquima pulmonar capaz de simular manobras de ventilação mecânica. O pulmão é um órgão estruturalmente complexo e hierarquizado. Por isso uma revisão bibliográfica multidisciplinar foi realizada. A revisão apresenta as propriedades mecânicas do parênquima, sua morfologia e os efeitos da tensão superficial para a contração e adesão dos septos interalveolares. Um aspecto importante para o modelo de alvéolo é um modelo de contato com adesão causada pela tensão superficial. Um modelo de contato adesivo simplificado foi desenvolvido e simulado em uma estrutura com parâmetros de mesma ordem de grandeza de um alvéolo real. A simulação foi realizada com o método dos elementos finitos não-linear e foi necessário empregar o método da corda para evitar divergência em pontos limites. Os resultados numéricos se aproximaram de resultados experimentais globai no pulmão com pressões de mesma ordem de grandeza. / This work is an important step on the development of a computational model of lung parenchyma capable to simulate mechanical ventilation maneuveres. The lung has a complex and hierarchized structure. Therefore a multidisciplinary literature review was held. The review presents the mechanical properties of the parenchyma, its morphology and the effects of surface tension to septums contraction and adhesion. An important aspect for the alveolus model is a contact model which includes the adhesion caused by surface tension. Thus a simplified model was developed and then simulated in a structure with properties of the same order of magnitude of a real alveolus. The simulation was performed with the nonlinear finite element method. The implementation of the arc-length method was also necessary in order to prevent diversion at limit points. The numerical results were close to whole lung experimental results with pressure levels of the same order of magnitude.

O emprego das estruturas de membrana é cada dia mais frequente em edificações de relevância civil e arquitetônica, em especial para a cobertura de grandes vãos. Sua aplicabilidade, contudo, vai além da construção civil, sendo igualmente importante nas indústrias das engenharias mecânica, naval, oceânica, aeroespacial e biomédica: aeronaves, satélites, paraquedas, airbags, velas de embarcações, moinhos de vento e até aplicações biomecânicas com tecidos humanos ou artificiais utilizam a tecnologia das estruturas de membrana. O comportamento mecânico de grande parte das membranas estruturais pode ser idealizado como uma casca isótropa de pequena espessura reforçada por uma membrana ortótropa. O objetivo desta pesquisa de mestrado é dar continuidade aos estudos referentes às teorias de casca geometricamente exatas desenvolvidas por [1], [2], [3] e [4], e sua generalização para o âmbito das membranas iniciada em [5]. Pretende-se contribuir, principalmente, para o desenvolvimento de equações constitutivas ortótropas consistentes para grandes deformações, apresentar uma metodologia para calibração dos parâmetros materiais destas equações constitutivas e para a análise de estabilidade com vistas ao estudo do fenômeno do enrugamento. / The use of membrane structures is becoming increasingly common in buildings of architectural and civil engineering importance, especially to cover large spans. Its applicability, however, goes beyond the construction industry being equally important in the industries of mechanical, naval, ocean, aerospace and biomedical engineering: aircraft, satellites, parachutes, airbags, sails of boats, windmills and even biomechanical applications with human and artificial tissues use the technology of membrane structures. The mechanical behavior of most structural membranes can be idealized as an isotropic thin shell reinforced by an orthotropic membrane. The objective of this work is to continue the studies on the geometrically exact shell theories developed by [1], [2] , [3] and [4] and its generalization to the scope of the membranes initially studied by [5]. It aims intended to contribute mainly to the development of consistent orthotropic constitutive equations for large deformations, to present a methodology for the calibration of the material parameters in those constitutive relations and for the stability analysis in order to study the phenomenon of wrinkling.

Este trabalho utiliza uma teoria não linear tridimensional de barras de [1], que considera o empenamento da seção transversal como uma das variáveis do problema, para a análise de perfis de aço por meio do Método dos Elementos Finitos. A equação constitutiva apresentada em [2] é aqui expressa de forma completa, com os termos de ordem cúbica e superiores antes desprezados. A influência da inclusão desses termos sobre a estabilidade da barra pode assim ser estudada com maior profundidade. O modelo constitutivo obtido foi implementado em um programa computacional de elementos finitos, e diversos exemplos numéricos são apresentados para validar a teoria aqui desenvolvida. Foi possível obter a instabilidade por flexo-torção com a equação constitutiva completa aqui desenvolvida, onde antes não era possível em [2]. Melhorando o modelo constitutivo. / This work presents a nonlinear tridimensional beam theory of [1], which encompasses cross-sectional warping as one of the model variables, and its application to the analysis of steel profiles with the Finite Element Method. The constitutive equation presented in [2] is fully developed, including the third and higher order terms neglected therein. The influence of the inclusion of such terms on the beam stability could then be studied more deeply. The obtained constitutive model was implemented in a computational Finite Element program and several numerical examples are presented to validate the theory here developed. The complete constitutive equation made possible to determine instability in composite flexural torsional states, a clear contribution to the model developed in [2].

Este trabalho é uma contribuição para o estudo de modelos de representação da não linearidade física de estruturas de concreto. Nele estudou-se o Modelo de Campo Distorcido de Tensões, ou Disturbed Stress Field Model - DSFM, através da implementação no programa MATLAB®. O programa desenvolvido implementa também o Modified Compression Field Theory - MCFT, que é a base para o DSFM. Neste trabalho também foi utilizado o programa VecTor2 desenvolvido na Universidade de Toronto. Neste documento é apresentada a base teórica das formulações estudadas, o programa desenvolvido é descrito brevemente e são apresentados resultados obtidos. / This research consists on the study of nonlinear behavioral models for reinforced concrete. The main focus of this study is the Disturbed Stress Field Model - DSFM. In order to study this model, it was implemented in the MATLAB® platform. The basis for the DSFM is the Modified Compression Field Theory MCFT, this model was also implemented using MATLAB®. The VecTor2 program developed at the University of Toronto was also used. This document presents the formulation for these analytical models, a description of the developed tools and the results obtained through them.

In this work, we present the formulation and implementation of two elastoplastic constitutive equations for kinematically exact thin-walled rod models. The first uses the fact that first order strains due to cross sectional shear stresses and warping are considered to formulate a small strains three-dimensional elastoplastic constitutive model. Given the kinematical hypothesis of non-deformability of the cross section in the projection of its plane, we may also assume that plastic deformations may occur due only to the cross sectional normal stresses, thereby allowing us to formulate a second, simple one-dimensional framework. Our approach adopts a standard additive decomposition of the strains together with a linear elastic relation for the elastic part of the deformation. Both ideal plasticity and plasticity with (linear) isotropic hardening are considered. The models have a computational implementation within a finite element thin-walled rod model and, following the kinematics adopted, we implement this equation on models with consideration of the warping of the cross sections, having 7 degrees of freedom. The formulation and implementation presented is validated by the analysis of problems known in the literature and comparison of the results. We believe that simple elastoplastic models combined with robust thin-walled rod finite element may be a useful tool for the analysis of thin-walled rod structures, such as, e.g., steel structures. / Neste trabalho, apresentamos a formulação e implementação de duas equações constitutivas elastoplásticas simples para modelos de barra de parede fina cinematicamente exatos. O primeiro usa o fato de deformações de primeira ordem devido a esforço cortante na seção transversal e empenamento serem considerados para formular um modelo constitutivo elastoplástico tridimensional para pequenas deformações. Dada a hipótese cinemática de não deformabilidade da seção transversal da barra na projeção de seu plano, podemos também assumir que deformações plásticas ocorrem devido apenas às tensões normais à seção transversal, nos permitindo formular um segundo modelo unidimensional simples. Nossa abordagem adota uma decomposição aditiva padrão das deformações com uma relação elástica linear para a parte elástica das deformações. Tanto plasticidade ideal quanto plasticidade com encruamento isótropo (linear) são considerados. Os modelos resultantes têm uma implementação computacional com elementos finitos de barras e, de acordo com a cinemática adotada, implementamos esta equação com consideração do empenamento das seções transversais, possuindo 7 graus de liberdade. A formulação e implementação apresentadas são validadas pela análise de problemas conhecidos na literatura e comparação dos resultados. Acreditamos que modelos elastoplásticos simples combinados com um elemento finito de barras robusto podem ser uma ferramenta útil para a análise de estruturas reticuladas como, por exemplo, estruturas de aço.

The pipelines used to transport oil and gas from the wellheads to the distribution and refining sites can be subjected to high levels of pressure and temperature. Under such conditions, the pipelines tend to expand, but, if the expansion is inhibited, a significant compressive axial force can arise, leading to their buckling, which can occur in the horizontal or vertical plane. In this context, the objective of the present work is to analyze the upheaval buckling of pipelines, considering the internal pressure to which they are subjected during the transportation of oil and gas as its only triggering. Using the concept of effective axial force, it aims at discussing two different approaches for considering the internal pressure in buckling problems: distributed loads dependent on pipeline curvature and equivalent compressive axial forces with follower and non-follower characteristics. It also discusses the influence of using static or dynamic analysis for such approaches. Concerning the upheaval buckling itself, the work intends to analyze and compare the influence of the soil imperfection amplitudes to the influence of the friction between the pipeline and the ground in the critical loads and in the post-buckling configurations of the pipeline. Besides theoretical research, the objectives are achieved through the development of various numerical models, since geometrically-simple models, without the consideration of the interaction between the pipeline and the ground, until more complex models, with the use of contact models to detect the ground and its imperfections. The models are developed in Giraffe (Generic Interface Readily Accessible for Finite Elements) using geometrically-exact finite element models of beams, undergoing large displacements and finite rotations. Through the research, it is concluded that there is an equivalence between the application of the internal pressure as a distributed load dependent on pipeline curvature and the application of the internal pressure as a follower compressive axial force. Besides this, it is demonstrated that the type of the analysis (static or dynamic) depends on the nature of the physical system analyzed. With the aid of results presented in terms of internal pressure, classical results about the influence of the imperfection amplitudes and of the friction between the pipeline and the ground in buckling are confirmed. It is also showed that the imperfection amplitudes analyzed play a more important role in the post-buckling configurations of the pipeline than the friction. / Os dutos utilizados para transportar petróleo e gás natural das reservas até os locais de distribuição e refino podem estar submetidos a elevados níveis de pressão e temperatura. Sob tais condições, os dutos tendem a se expandir, porém, se a expansão é inibida, uma força axial de compressão significativa pode surgir nos dutos, ocasionando a flambagem lateral ou vertical dos mesmos. Dentro desse contexto, o objetivo do presente trabalho é analisar a flambagem vertical de dutos, considerando a pressão interna à qual eles estão submetidos durante o transporte de petróleo e gás natural como o único parâmetro desencadeador da flambagem. Usando o conceito de força axial efetiva, o trabalho objetiva discutir duas abordagens diferentes para considerar a pressão interna nos problemas de flambagem: carregamentos distribuídos dependentes da curvatura do duto e forças axiais de compressão equivalentes à pressão com caráter seguidor e não seguidor. O trabalho também discute a influência de usar a análise estática ou dinâmica para analisar essas abordagens de carregamento. Com relação à flambagem vertical propriamente dita, o trabalho pretende analisar e comparar a influência das amplitudes das imperfeições presentes no solo com a influência do atrito entre o duto e o solo nas cargas críticas e nas configuração pós-críticas do duto. Além de pesquisa teórica, os objetivos são atingidos através do desenvolvimento de vários modelos numéricos, desde modelos geometricamente simples, sem a consideração da interação entre o duto e o solo, até modelos mais complexos, com o uso de modelos de contato para detectar o solo e suas imperfeições. Os modelos são desenvolvidos no Giraffe (Generic Interface Readily Accessible for Finite Elements) usando elementos finitos geometricamente exatos de viga, sujeitos a grandes deslocamentos e rotações finitas. Através da pesquisa, conclui-se que existe uma equivalência entre a aplicação da pressão interna como um carregamento distribuído dependente da curvatura do duto e a aplicação da pressão interna como uma força axial de compressão seguidora. Além disso, demonstra-se que o tipo de análise (estática e dinâmica) depende da natureza do sistema físico analisado. Com a ajuda de resultados apresentados em termos de pressão interna, os resultados clássicos sobre a influência das amplitudes das imperfeições e do atrito entre o duto e o solo são confirmados. Também é mostrado que as amplitudes das imperfeições analisadas desempenham uma maior influência nas configurações pós-críticas do duto do que o atrito.