Esta dissertação tem como objetivo apresentar e verificar uma série de métodos de análise numérica não convencional em relação ao método de elementos finitos (MEF) utilizado comercialmente. Inicialmente é realizada uma revisão bibliográfica para posterior construção do referencial teórico. Dentre os itens abordados neste referencial teórico estão uma breve descrição e histórico de vasos de pressão, método MEFG (Elementos Finitos Degenerados), o método H-FEM (Elementos Finitos de Hermite), método Stable GFEM (Metódo de Elementos Finitos Generalizados com Estabilização) e o método FEM B-Spline (Método de Elementos Finitos B-Spline). No referencial teórico é construído todo equacionamento de cada método numérico para posterior programação de cada um no programa Matlab®. Os resultados gerados abordam os estudos dos efeitos de travamento devido ao esforço cisalhante sobre a primeira frequência natural em placas circulares, análise dos efeitos de regularidade dos espaços de aproximação obtidos segundo os métodos citados anteriormente, estudo do erro relativo das frequências naturais em relação à solução de referência para um intervalo definido por 50% dos primeiros modos aproximados numericamente e estudos dos efeitos de vibração forçada produzido por forças de impulso aplicadas no plano de simetria radial de cascas de revolução cilíndricas e esféricas. / This dissertation has as objective presents and checks a series of non-conventional methods for numerical analysis in comparison with the finite element method (FEM) applied commercially. Initially is performed a bibliographic revision for a posterior building of theoretical reference. Among the approached items of this theoretical reference are a quick description and history about pressure vessels, GFEM method, H-FEM method, Stable GFEM method and FEM B-Spline method. On the theoretical reference is built all equation basis of each method for posterior programming of each one in the software Matlab®. The generated results approaches the studies of locking due shear forces in circular plates under the first natural frequency, analysis of approximation spaces regularity obtained according to the methods cited before, study of relative error of natural frequencies than the reference solution for an interval of 50% of first modes numerically approximated and studies of effects due forced vibration generated by impulse forces applied in the radial symmetric plan of cylindrical shells and spheres.

Neste trabalho são descritas análises estática e dinâmica em regime elástico de pórticos espaciais com incerteza no carregamento utilizando o método de elementos finitos (MEF). A abordagem do problema dinâmico estrutural de vibrações forçadas será feita dentro do regime elástico linear com amortecimento viscoso. O pórtico é abordado com teoria de viga fina onde são considerados os acoplamentos devido às forças normais e a torção sobre o modelo de flexão em viga fina de Bernoulli-Euler. A formulação forte do problema dinâmico da viga com carga axial é obtida a partir das equações de equilíbrio de Euler-Lagrange, decorrentes do princípio variacional de Hamilton e, a formulação forte do problema dinâmico do eixo longo sujeito a torção é determinada utilizando o modelo de torção de Saint-Venan. Para a determinação da formulação fraca para o elemento de viga tridimensional foi aplicado o método dos resíduos ponderados Galerkin. Também é descrita neste trabalho, a análise da confiabilidade estrutural considerando a incerteza no carregamento e, em algumas propriedades mecânicas dos materiais, utilizando o método de Monte Carlo (MMC). A série de Neumann foi utilizada como alternativa para reduzir o tempo de processamento do problema dinâmico. A análise simultânea das diversas variáveis foi abordada utilizando a estatística multivariada. Os resultados da analise estática e de vibrações livres dos exemplos numéricos são apresentados com o intuito de validar os métodos contidos neste trabalho comparando-os com resultados obtidos utilizando um software comercial de análise estrutural. / In this thesis, static and dynamic analysis of elastic space frames with regime uncertainty in loading using the finite element method (FEM) are described. A structural approach to the dynamic problem of forced vibrations will be made within the linear elastic regime with viscous damping. The structure is approached with thin beam theory which the couplings are considered normal forces due to twisting and bending about the model thin Bernoulli-Euler beam. The strong formulation of the dynamic problem of the beam with axial load is obtained from the equilibrium equations of Euler-Lagrange equations arising from the variational principle of Hamilton and the strong formulation of the dynamic problem of the long shaft under torsion is determined using the model of twist of Saint-Venan. To determine the weak formulation for three-dimensional beam element applies the Galerkin method of weighted residues. Is also described in this study, the analysis of structural reliability considering the uncertainty in loading and some mechanical properties of materials using the Monte Carlo method (MMC). Neumann series will be used as an alternative to reduce the processing time of the dynamic problem. The simultaneous analysis of several variables is addressed using multivariate statistics. Some results of numerical examples are presented in order to validate the methods contained in this work compared with the results obtained from structural analysis of commercial software.

Orientador: Kátia Lucchesi Cavalca Dedini / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecânica / Made available in DSpace on 2018-08-16T19:43:14Z (GMT). No. of bitstreams: 1 Tuckmantel_FelipeWenzeldaSilva_M.pdf: 5340309 bytes, checksum: a62c87427ed691bb455fd23f46cde6f0 (MD5) Previous issue date: 2010 / Resumo: O estudo de máquinas rotativas ocupa uma posição destacada no contexto de máquinas e estruturas, tendo em vista a grande aplicação das máquinas rotativas e a importância de conhecer o seu funcionamento com fins de otimização e de projeto para diminuição de custos, melhora de qualidade e aumento da vida útil dos componentes. Estas máquinas, e seus componentes, possuem características peculiares em seu funcionamento, as quais se refletem em seu comportamento dinâmico e, conseqüentemente, nos problemas mecânicos a estes associados. A análise do comportamento das máquinas rotativas envolve considerável complexidade, já que são constituídas, geralmente, por um sistema de equipamentos que abrangem, além do rotor, eixos e mancais, fundação e outros equipamentos auxiliares, como, por exemplo, selos de fluxo. Portanto, para se estudar o comportamento dinâmico desses sistemas, é necessário determinar a interação de todos os componentes para o entendimento dos fenômenos envolvidos. O sistema rotativo foi matematicamente modelado através da modelagem do eixo por elementos finitos, através de um conjunto de discos rígidos e segmentos de eixo com massa e elasticidade distribuídas, e dos mancais hidrodinâmicos, cuja representação ocorre a partir de seus coeficientes equivalentes de rigidez e amortecimento. A fundação também teve sua influência contemplada pelas representações por Parâmetros Concentrados, Impedância Mecânica e Coordenadas Mistas. Portanto, neste trabalho, desenvolveu-se a configuração do pacote computacional, a padronização da entrada e saída dos arquivos, a implementação dos elementos de viga e de disco rígido, a excitação por desbalanceamento, a determinação dos carregamentos nos mancais e a resposta dinâmica do sistema no domínio da freqüência, tendo como objetivo robustez, flexibilidade, confiabilidade, além de usabilidade e interface amigável do pacote / Abstract: The Rotordynamics occupies a prominent position in the design of machines and structures, due to the wide application of rotating machinery and the importance of the knowledgment of its operation for purposes of optimization and design for cost reduction, quality improvement and increased component life cycle. These machines and their components have peculiar characteristics in its operation, which will reflect in its dynamic behavior and, consequently, on the mechanical problems associated with these. The behavior of rotating machinery involves considerable complexity, as they are composed, usually by a system of equipment which includes, besides the rotor, shafts and bearings, foundation and other auxiliary equipment, for example, shaft seals. Therefore, to study the dynamic of these systems, it is necessary to determine the interaction of all components for understanding the phenomena involved. The rotation system was mathematically modeled by the shaft modeling through the Finite Element Method, through a set of disks and shaft segments with elasticity and mass distribution, and hydrodynamic bearings, whose representation occurs through the coefficients of equivalent stiffness and damping. The foundation also had its influence covered by the formulation of Concentrated Parameters, Mechanical Impedance and Mixed Coordinates. Therefore, in this work, it is developed the computational package configuration, standardization of input and output files, the implementation of beam elements and disks, the unbalance excitation, the determination of loads on the bearings and the frequency dynamic response of the system, aiming at robustness, flexibility, reliability, usability and user-friendly of the software / Mestrado / Mecanica dos Sólidos e Projeto Mecanico / Mestre em Engenharia Mecânica

Orientadores: José Roberto Nunhez, Nicolas Spogis / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Química / Made available in DSpace on 2018-08-17T20:12:44Z (GMT). No. of bitstreams: 1 Nasato_DanielSchiochet_M.pdf: 8429075 bytes, checksum: 08b5ea34c892e0a82506be022a1613ee (MD5) Previous issue date: 2011 / Resumo: Este trabalho propõe um acoplamento entre o método dos elementos discretos (DEM), utilizado para simular o manuseio de materiais particulados, e o método dos elementos finitos (FEM), utilizado para estimar respostas estruturais a uma dada condição de carregamento. O acoplamento é de 1 via, onde as forças na estrutura calculadas pelo software de DEM são passadas ao software de FEM. Foi proposta uma validação numérica a fim de avaliar a precisão do acoplamento em passar as forças do software de DEM para o software de FEM. Apenas efeitos estáticos são avaliados neste primeiro teste. Uma validação experimental foi realizada utilizando-se medidas em um tanque de mistura, permitindo assim a captura de efeitos dinâmicos das partículas. Com o modelo desenvolvido e validado, foram realizadas simulações em um equipamento industrial utilizando os softwares de DEM e FEM e o acoplamento desenvolvido, como parte de um projeto realizado com uma empresa do setor de mineração / Abstract: This work proposes a coupling between the discrete element method (DEM), used to simulate bulk materials handling, and the finite element method (FEM), used to simulate structural response to a certain loading condition. It is a one-way coupling, where the forces in the structure calculated by the DEM software are sent to the FEM software. It was proposed a numerical validation in order to evaluate the precision of the coupling to pass the forces from the DEM to the FEM software. Only static effects are evaluated in this initial test. An experimental validation was done using data measured from a mixing tank, allowing dynamic effects to be considered. With the model developed and validated simulations were done in industrial equipment using DEM and FEM software and the developed coupling, as part of a project with a mining sector company / Mestrado / Desenvolvimento de Processos Químicos / Mestre em Engenharia Química

Orientador: Philippe Remy Bernard Devloo / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Civil, Arquitetura e Urbanismo / Made available in DSpace on 2018-08-18T01:44:23Z (GMT). No. of bitstreams: 1 Cecilio_DiogoLira_M.pdf: 5163879 bytes, checksum: 102fa941f977c8a36107d3d37f05c5af (MD5) Previous issue date: 2011 / Resumo: Este estudo revisa os fundamentos da teoria da plasticidade incremental em uma formulação elastoplástica do comportamento no ponto material e sua ligação com o método dos elementos finitos. Documenta a implementação na linguagem de programação C++ dos modelos elastoplásticos de Drucker-Prager, Mohr-Coulomb e von Mises. É realizado uma série de simulações de carregamento uniaxial para os modelos implementados e finalmente solucionados três problemas elastoplásticos em elementos finitos. Em cada um dos problemas resolvidos foi utilizado o modelo elastoplástico adequado ao material do problema em questão. Os problemas solucionados foram todos em estado plano de deformação e consistiram em: i) um ensaio de compressão diametral; ii) um problema de estabilidade de taludes; iii) e um tubo de aço pressurizado internamente. A implementação da formulação elastoplástica foi validada comparando os resultados das simulações com os da literatura / Abstract: This work reviews the fundamentals of the theory of plasticity in an incremental formulation of the elastoplastic behavior at a material point and its connection with the finite element method. It documents the implementation of the elastoplastic models of Drucker-Prager, Mohr-Coulomb and von Mises in the programming language C++. A series of simulations of uniaxial loading models are carried out and three problems in elastoplastic finite element method are solved. In each of the modeled cases the most appropriate elastoplastic model among those implemented is employed. The problems were solved in plane strain state and comprise: i) a diametral compression test of concrete or Brazilian Disc Test; ii) a problem of slope stability; iii) and an internally pressurized vessel. The implementation of the elastoplastic formulation was validated by comparing the simulation results with the results from the literature / Mestrado / Estruturas / Mestre em Engenharia Civil

Orientador: Renato Pavanello / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecânica / Made available in DSpace on 2018-08-18T00:16:08Z (GMT). No. of bitstreams: 1 Bispo_RafaelSantana_M.pdf: 3924997 bytes, checksum: eedf2e79312457d2dfe0c4ce2418dfe4 (MD5) Previous issue date: 2011 / Resumo: Atualmente, existe uma grande tendência no incremento da produção de energia elétrica através de fontes renováveis. Em especial, a geração de energia elétrica produzida através de parques eólicos tem sido bastante adotada. O projeto desses equipamentos envolve a modelagem dinâmica acoplada solo-fluido-estrutura que pode ser estudada usando-se a formulação particionada, onde o problema da interação entre os meios é tratado de maneira iterativa. Nesse tipo de técnica, é possível que modelos fisicamente heterogêneos, chamados de partições, possam utilizar diferentes técnicas de discretização, como por exemplo o domínio do fluido ser baseado em uma formulação de Elementos de Contorno e o domínio estrutura baseado em uma formulação em Elementos Finitos. Neste trabalho, é realizado um estudo dinâmico de turbinas eólicas, utilizando tratamento particionado e Multiplicadores de Lagrange afim de se obter as frequências características e as curvas de resposta em frequência do sistema em análise. A discretização do problema é realizada através do Método dos Elementos Finitos (FEM) utilizando elemento de pórtico e quadrilateral de Wilson. Desta forma, a resolução de problemas de interação, utilizando a formulação particionada, é estudada com a finalidade de avaliar a convergência e a viabilidade da técnica em problemas harmônicos estruturais / Abstract: Currently, there is a great tendency in increasing the production of electricity through renewable sources. In this context, the generation of electric energy produced by wind farms has been widely adopted. The design of these devices involves the dynamic modeling of coupling fluid-structure-soil that can be studied using the partitioned formulation, where the problem of interaction between the parties is iterative manner. In this type of technique, it is possible that physically heterogeneous models, called partitions, can use different discretization techniques, such as the domain of fluid is based on a formulation of boundary element or based on a finite element formulation. In this paper a harmonic simulation of wind turbines, using partitioned treatment and Lagrange multipliers is studied in order to obtain the characteristic frequencies and frequency response function of the system under analysis. The discretization of the problem is performed using the Finite Element Method (FEM), as well as beam elements and the quadrilateral Wilson element. Thus, the resolution of elastodynamics problems, using the partitioned formulation is studied with the aim of assessing the feasibility and convergence of this technique, applied to dynamic harmonic analysis / Mestrado / Mecanica dos Sólidos e Projeto Mecanico / Mestre em Engenharia Mecânica

Orientador: Hugo Enrique Hernández Figueroa / Tese (doutorado) - Universidade Estadual de Campinas, Faculdade de Engenharia Elétrica e de Computação / Made available in DSpace on 2018-08-18T02:45:36Z (GMT). No. of bitstreams: 1 MoresJunior_JoseAmilton_D.pdf: 15207023 bytes, checksum: 112938cdefab40754d08619dd4e4bb3a (MD5) Previous issue date: 2011 / Resumo: Uma estratégia numérica eficiente para calcular os parâmetros de dispersão de alta ordem em guias ópticos e apresentada. No que e de nosso conhecimento, pela primeira vez na literatura um estudo sistemático dos erros envolvidos no processo de calculo numérico de tais parâmetros e feito, mostrando que a presente estratégia e capaz de modelá?los com precisão superior aquela que pode ser fornecida pelo estado da arte dos processos de fabricação. Tal abordagem combina um simulador modal completamente vetorial baseado no Método dos Elementos Finitos e um algoritmo apropriado de diferenciação por diferenças finitas. O seu desempenho foi cuidadosamente validado por meio da analise de varias geometrias?chaves. Adicionalmente, a otimização dos parâmetros de dispersão de alta ordem pode ser feita acoplando?se ao esquema um algoritmo genético, como mostrado aqui por meio do projeto de Fibras de Cristais Fotônicos apropriadas para aplicações em amplificação paramétrica. Por fim, e proposto um novo Modelo Vetorial de Quatro Ondas para Guias com Alto Contraste de Índice de Refração que e deixado para servir de base em trabalhos futuros / Abstract: An efficient numerical strategy to compute the higher-order dispersion parameters of optical waveguides is presented. To the best of our knowledge, for the first time in the literature a systematic study of the errors involved in the higher-order dispersions' numerical calculation process is made, showing that the present strategy can model those parameters with greater accuracy than the one that can be provided by optical waveguides' state of the art manufacturing processes. Such strategy combines a full-vectorial finite element modal solver and a proper finite difference differentiation algorithm. Its performance has been carefully assessed through the analysis of several key geometries. In addition, the optimization of those higher-order dispersion parameters can also be carried out by coupling to the present scheme a genetic algorithm, as shown here through the design of Photonic Crystal Fibers suitable for parametric amplification applications. Finally, a new Four-Wave Vectorial Model for High Refractive-Index-Contrast Optical Waveguides is proposed and left to be used as a basis for future works / Doutorado / Telecomunicações e Telemática / Doutor em Engenharia Elétrica

Orientador: Renato Pavanello / Tese (doutorado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecanica / Made available in DSpace on 2018-08-05T04:07:44Z (GMT). No. of bitstreams: 1 Souto_Carlosd'Andrade_D.pdf: 10416734 bytes, checksum: 002d9b3eec5201986dcb26153a06af7f (MD5) Previous issue date: 2005 / Resumo: Neste trabalho é desenvolvido um método modal de resolução para modelos de elementos finitos aplicados em problemas dinâmicos elasto - acústicos cíclico simétricos. O método combina a teoria de propagação de ondas em sistemas cíclico - simétricos com o método de subestruturaçãocom interfaces- fixas de Craig Bampton e foi desenvolvido originalmente para estruturas e sistemas hidroelásticos, sendo extendido aqui para o caso elasto - acústico.Um programa de elementos finitos escrito em C++ utilizando a filosofia da programação orientada a objetos desenvolvido no Departamento de Mecânica Computacional foi utilizado para realizar as implementações computacionais. Como pré-processador foi utilizado o ANSYS também utilizado no pós-processamento, juntamente com o MATLAB. Foram adotadas como variáveis deslocamentos para o domínio estrutural e potencial de deslocamentos para o domínio fluido. Foi implementada a resolução de problemas multidomínios (estrutura, fluido e interface) onde a interface é um domínio puramente geométrico. Também são analisados métodos de redução modal e subestruturação para sistemas fluido - estrutura sendo proposto um método de redução modal simétrica para problemas elasto - acústicos. Diversos exemplos numéricos são analisados e os resultados são comparados com soluções numéricas dos modelos completos obtidas por meio do ANSYS / Abstract: ln this work a modal method to perform eigenanalysis of elastoacoustic cyclic - symmetric systems by the finite element method is developed. The method combines component mode synthesis and wave propagation in cyclic symmetric domains was originally developed to structural and hydroelastic systems and is now extended to the elastoacoustic problem. A C++ object oriented finite element code developed by the Computational Mechanics Department was used on the computational implementation. ANSYS was used on pre-processing work and MATLAB and ANSYS were used to post-process the results. Displacements and potential of displacements were choosen as variables for the structure and the fluid respectively, leading to a symmetric coupled system. Using a multi-domain approach, the structure and fluid domains are necessarily coupled by the interface, a geometric domain. Modal reduction and substructuring methods for elastoacoustic problems were also investigated and a symmetric modal reduction is proposed. Some numerical examples were analyzed and the results are compared to analytical solutions andANSYS results / Doutorado / Mecanica dos Sólidos e Projeto Mecanico / Doutor em Engenharia Mecânica

Orientador: João Frederico da Costa Azevedo Meyer / Dissertação (mestrado) - Universidade Estadual de Campinas, Instituto de Matematica, Estatistica e Computação Cientifica / Made available in DSpace on 2018-08-05T12:58:22Z (GMT). No. of bitstreams: 1 Salvatierra_MarcosMarreiro_M.pdf: 1600898 bytes, checksum: a8beabb556b24c734508735602125989 (MD5) Previous issue date: 2005 / Resumo: A proposta deste trabalho é criar um modelo para descrever computacionalmente o convívio entre duas espécies competidoras com características de migração na presença de um material impactante tóxico. As equações a serem utilizadas deverão incluir os fenômenos de dispersão populacional, processos migratórios, dinâmicas populacionais densidade-dependentes e efeitos tóxicos de um material impactante evoluindo no meio, provocando um decaimento proporcional. Recorrendo a um instrumental consagrado, embora com desenvolvimento relativamente recente, será usado um sistema clássico do tipo Lotka-Volterra (conseqüentemente não-linear) combinado a Equações Diferenciais Parciais de Dispersão-Migração. O primeiro passo é a formulação variacional discretizada deste sistema visando o uso de Elementos Finitos combinados a um método de Crank-Nicolson. Em segundo lugar, virá a formulação de um algoritmo (conjuntamente com sua programação em ambiente MATLAB) que aproxima as soluções discretas relativas a cada população em cada ponto e ao longo do intervalo de tempo considerado nas simulações. Por fim, serão obtidas saídas gráficas úteis dos pontos de vista quantitativo e qualitativo para uso em conjunto com especialistas de áreas de ecologia e meio ambiente na avaliação e na calibração de modelos e programas, bem como no estudo de estratégias de preservação, impacto e recuperação de ambientes / Abstract: The purpose of this work is to create a model to computationally describe the coexistence of two competing species with migration features in the presence of a toxic impactant material. The equations must include the phenomena of populational dispersion, migratory processes, density-dependent populational dynamics and toxic effects of the evolutive presence of an impactant material developing in the environment, generating a proportional decrease in both populations. Resorting to well-established, although relatively recent, mathematical instruments a Lotka - Volterra type (and consequently nonlinear) system, including characteristics of a Migration-Dispersion PDE. The first step is the discrete variational formulation of this system aiming for the use of the Finite Element Method toghether with a Crank-Nicolson Method. Second, the formulation of an algorithm (together with a programme in MATLAB environment) that approximates the relative discrete solutions to each population in each point and along of the time interval considered in the simulations. Lastly, useful graphics will be obtained of the quantitative and qualitative viewpoints for use with specialists of the fields of ecology and environment and in the evaluation and calibration of models and programmes / Mestrado / Mestre em Matemática Aplicada

Orientador: Alberto Luiz Serpa / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecanica / Made available in DSpace on 2018-08-05T19:36:12Z (GMT). No. of bitstreams: 1 Paulino_LeonardoPereira_M.pdf: 2501868 bytes, checksum: 87eef8b309aa247d3e5c6c4f572b1323 (MD5) Previous issue date: 2004 / Resumo: Em muitos problemas não lineares na mecânica estrutural é necessário o emprego de métodos de solução que controlam os deslocamentos e o nível de carregamento simultaneamente. Um dos métodos mais conhecidos para estas situações é o método do arc-length. Este método consiste na introdução de uma restrição (" esfera generalizada", por exemplo) à equação de equilíbrio não linear do problema. A solução do problema de equilíbrio juntamente com a restrição" esférica" mantém o deslocamento limitado, evitando a divergência durante o processo de resolução. O método clássico do arc-length apresenta certa complexidade por envolver duas fases de cálculo, predição e correção da solução, além do fato de que mais de uma solução satisfaz as equações envolvidas, sendo necessário o emprego de alguns critérios ainda não suficientemente robustos para a escolha da solução adequada. Propõe-se neste trabalho o estudo de uma formulação alternativa para o método do arc-length através do emprego dos conceitos de programação matemática, onde um problema de minimização com restrições é formulado, minimizando-se uma função objetivo que caracterize o problema de equilíbrio e esteja sujeita a restrições que garantam o controle dos deslocamentos. Além disso, uma restrição adicional é imposta ao problema para buscar a unicidade da solução. A verificação do desempenho da metodologia proposta é feita através da solução de problemas de estruturas reticuladas envolvendo efeitos de grandes deslocamentos. A implementação foi feita em Matlab pela facilidade e disponibilidade de ferramentas computacionais tais como as funções do "toolbox" de otimização / Abstract: Non-linear problems in structural mechanics in general require the use of solution methods that control the displacements and the load leveI simultaneously. One of the most popular method used in these cases is the arc-length method. The arc-length method introduces one additional constraint equation to the non-linear equilibrium equations of the problem. The solution of the equilibrium equations with the additional constraints should keep a limited displacement to avoid numerical divergence during the solution procedure. The classical arc length method involves two solution phases, prediction and correction. A difficulty of this method is that it can re-compute a solution already determined, requiring the use of some cri teria that are not robust enough to choose the adequated solution. This work proposes an alternative formulation for the arc-length method using the concepts of mathematical programming, where a constrained minimization problem is formulated. The objective function is established in terms of the equilibrium residue and the arc-length constraint should ensure displacements and load leveI control. Besides this, an additional constraint equation is employed to guarantee the uniqueness of the solution, improving the solution path compared to the classical arc-length method. To verify the proposed methodology computational performance some problems involving the effects of large displacements in truss structures are presented. The implementation was done in the software Matlab due to its facility and the computational tools that are available, such as the optimization toolbox functions / Mestrado / Mecanica dos Sólidos e Projeto Mecanico / Mestre em Engenharia Mecânica