Orientador: Prof. Dr. Juan Pablo Julca Avila / Dissertação (mestrado) - Universidade Federal do ABC, Programa de Pós-Graduação em Engenharia Mecânica, 2017. / Este trabalho de mestrado tem como principal objetivo o desenvolvimento de um código computacional para análise estática tridimensional de vigas longas submetidas a cargas de torção e tração. Para tal, um elemento de viga espacial de Euler, com pequenas deformações, foi desenvolvido usando a abordagem corrotacional. Um esquema para o tratamento das rotações finitas foi implementado baseando-se em pseudovetores e quaternions. As configurações de equilíbrio foram obtidas com base em dois esquemas incrementais iterativos: o método de Newton Raphson versão total e o método de comprimento de arco cilíndrico com controle de carga e deslocamento. Critérios de convergência baseados em deslocamentos, energia de deformação e força residual foram considerados. Foram simuladas vigas espaciais retas com imperfeições geométricas submetidas a carregamentos combinados de torção e compressão, obtendo-se resultados comparativos satisfatórios. Simulações com imposição de deslocamentos no topo de uma viga curva em catenária também foram realizadas e comparadas satisfatoriamente. O código computacional desenvolvido é capaz de construir / This master's work has as main purpose the development of a computational code for three-dimensional static analysis of slender beams subjected to torsional and tensile loads. For this, an Euler space beam element with small deformations was developed using the corrotational approach. A scheme for the finite rotations treatment was implemented based on pseudovectors and quaternions. The equilibrium configurations were obtained based on two iterative incremental schemes: the Newton Raphson method total version and the cylindrical arc length method with load and displacement control. Convergence criteria based on displacements, strain energy and residual force were taken into account. Direct spatial beams with geometric imperfections subjected to combined loading of torsion and compression were simulated in order to obtain satisfactory comparative results. Simulations with imposition of displacements on top of a curved catenary beam were also performed and compared satisfactorily. This computational code is able to construct slender beam equilibrium configurations until the beginning of the stability loss.

Esse trabalho trata da análise computacional de problemas aeroelásticos. Os casos de interesse envolvem escoamentos compressíveis em torno de estruturas de cascas de materiais compósitos laminados. A solução do problema de interação fluido- estrutura é conduzida por meio do esquema particionado denominado improved serial staggered procedure (ISS), o que permite a utilização de algoritmos independentes para a análise de cada componente do sistema. Um elemento triangular plano com três nós destinado à análise de cascas de materiais compósitos laminados é obtido com base na formulação de dois elementos de alta performance desenvolvidos para a análise de membranas e de placas. A flexibilidade ao corte transversal, a ausência de travamento por razão de aspecto em problemas com flexão no plano da casca são características do elemento. A análise dinâmica não linear é realizada através do procedimento corrotacional de conservação aproximada de energia (AECCP) que foi implementado com base na formulação corrotacional independente de elemento (element independent CR formulation – EICR). Essa abordagem permite que problemas transientes com grandes rotações e forte não linearidade geométrica sejam integrados por longos períodos de tempo, mantendo a estabilidade, a precisão e a eficiência da solução. A simulação do escoamento é realizada através do método explícito de Taylor-Galerkin de dois passos utilizando elementos finitos tetraédricos em malhas não estruturadas para a discretização espacial. As equações governantes do escoamento são dadas segundo a formulação lagrangeana-euleriana arbitrária (ALE), permitindo que malhas móveis sejam empregadas na análise dos problemas aeroelásticos. Um método simples e eficiente é adotado para a movimentação da malha do fluido com base na distância do nó às superfícies de contorno do escoamento. A utilização de malhas não coincidentes na discretização do fluido e da estrutura é possibilitada pelo emprego do esquema de projeção nodal para a transferência de informações entre os meios na região de interface. A apresentação, implementação e verificação de cada componente da ferramenta de análise é realizada independentemente graças à natureza do esquema de solução particionado. Por fim são analisados problemas de interação fluido-estrutura, onde as respostas obtidas são comparadas com resultados experimentais e numéricos de outros autores. Também são propostos e analisados novos problemas envolvendo estruturas de compósitos laminados. / This work deals with computational analysis of aeroelastic problems, specifically those related to compressible flows around laminated composite shell structures. The fluid- structure problem is solved by a partitioned scheme called improved serial staggered procedure (ISS), allowing to use independent algorithms for the analysis of each system component. A 3-node triangular flat element for laminated composite shells is obtained based on the formulation of two high performance plate and membrane elements. Shear flexibility and no aspect ratio locking on in-plane bending problems are some of the element characteristics. The nonlinear dynamic analysis is performed using the approximately energy conserving corotational procedure (AECCP), which is implemented based on the element independent corotational formulation (EICR). This approach enables long time transient problems with strong nonlinearities and large rotations to be efficiently solved while keeping stability and accuracy. The flow simulation is performed by the explicit two step Taylor-Galerkin method employing tetrahedral finite elements in unstructured meshes for the space discrezation. The flow governing equations are given in arbitrary lagrangian-eulerian formulation, enabling aeroelastic simulations with moving meshes. A simple and efficient method is adopted to move the fluid mesh based on the distance of each node to the boundary surfaces. Transference of information between non-matching fluid and structure meshes are performed by the node-projection scheme. Each component of the aeroelastic analysis tool is independently presented, implemented and verified due to the partitioned nature of the scheme adopted for the solution of the coupled system. Some aeroelastic problems are analyzed and results are compared to other experimental and numerical works. New problems of aeroelastic analysis of laminated composite structures are proposed and solved.

Conselho Nacional de Desenvolvimento CientÃfico e TecnolÃgico / A utilizaÃÃo de vigas laminadas de parede fina nas Engenharias AeronÃutica, Civil, MecÃnica e Naval tem crescido bastante nos Ãltimos anos, devido a suas elevadas relaÃÃes rigidez/peso e resistÃncia/peso. Esses elementos estruturais normalmente apresentam paredes finas devido à alta resistÃncia dos materiais compÃsitos. Outra caracterÃstica importante à que, mesmo sem apresentar grandes deformaÃÃes e sem que o limite elÃstico do material seja ultrapassado, estas vigas apresentam comportamento nÃo linear geomÃtrico devido à sua elevada esbeltez, acarretando em grandes deslocamentos e rotaÃÃes. Dependendo da laminaÃÃo utilizada, as vigas de materiais compÃsitos podem apresentar diversos acoplamentos entre esforÃos e deformaÃÃes, tornando a sua anÃlise bem mais complexa do que a anÃlise de vigas de material isotrÃpico. Neste trabalho, foram desenvolvidos dois elementos finitos de pÃrtico espacial para anÃlise nÃo linear geomÃtrica de vigas laminadas de parede fina. As propriedades seccionais da viga sÃo avaliadas atravÃs de teorias de vigas laminadas de parede fina apropriadas, em que desprezam-se os efeitos do empenamento e do cisalhamento transversal. Tais teorias de vigas laminadas conduzem a uma matriz constitutiva 4x4, onde diferentes nÃveis de acoplamento entre esforÃos e deformaÃÃes de viga sÃo considerados, por exemplo, pode-se obter uma matriz constitutiva diagonal ou cheia. A abordagem corrotacional independente do elemento à utilizada para o tratamento de grandes deslocamentos e rotaÃÃes de corpo rÃgido no espaÃo. No Ãmbito local, sÃo utilizados dois elementos, um baseado na teoria linear e outro na descriÃÃo Lagrangeana Total. O tratamento matemÃtico das grandes rotaÃÃes no espaÃo à realizado por meio do tensor das rotaÃÃes (fÃrmula de Rodrigues), juntamente com o conceito do pseudovetor. As implementaÃÃes dos elementos finitos propostos neste trabalho sÃo realizadas no software de cÃdigo aberto FAST. A metodologia de trabalho segue o roteiro clÃssico de mÃtodos computacionais, incluindo formulaÃÃo, implementaÃÃo, verificaÃÃo e validaÃÃo dos resultados. A verificaÃÃo à realizada atravÃs dos modelos tridimensionais de elementos finitos de casca e sÃlido desenvolvidos no software comercial ABAQUS. A validaÃÃo à realizada por meio da comparaÃÃo com resultados de ensaios experimentais encontrados na literatura. No que diz respeito à avaliaÃÃo das propriedades seccionais, pode-se verificar uma Ãtima concordÃncia entre as teorias de vigas laminadas adotadas neste trabalho e os resultados numÃricos e de ensaios experimentais, para todas as laminaÃÃes e carregamentos considerados. No caso dos elementos de pÃrtico espacial, verificou-se uma Ãtima concordÃncia entre os resultados dos elementos finitos propostos neste trabalho e os resultados analÃticos e computacionais disponÃveis na literatura. Observa-se tambÃm que o elemento baseado na descriÃÃo Lagrangeana à mais eficiente do que o elemento baseado na teoria linear no que tange à capacidade de apresentar uma resposta satisfatÃria com uma malha menos refinada. / The use of thin walled laminate beams in Aeronautical, Civil, Mechanical and Naval Enginee- ring is increasing in the last years. This is due to their high stiffness/weight and strength/weight ratios. Composite beams and other structural elements tend to have thin walls due to the elevated strength of the material. Other important aspect is that, even without reaching large strains and without overcoming the elastic limit of the material, such beams present geometric nonlinear behavior due to high their slenderness, leading to large displacements and rotations. Depen- ding on the composite layup, the beams of composite materials can present several couplings between generalized stresses and strains, requiring a more complex analysis procedure when compared to isotropic beams. In this work, two three-dimensional space frame finite elements that can be used to analyze composite thin-walled beams subjected to geometric non-linearity were developed. The cross-section properties of the beams are evaluated through suitable thin walled beam theories, where the effects of the warping and transverse shear are neglected. Such theories yield a 4x4 constitutive matrix for the laminate and different levels of coupling between generalized stresses and strains can be considered. Depending of such couplings, the constitu- tive matrix can either be full or diagonal. The element independent corotational approach was used in order to consider large displacaments and rigid body rotations in space. In the local coordinate system, two elements are used, one based on the linear strain theory and the other on the Total Lagrangian formulation. The mathematical treatment of the large rotations in the space is performed by means of the rotation tensor (Rodriguesâs formula) in conjunction with the concept of the pseudovector. The computational implementations of the two finite elements proposed in this work were done in the open source software FAST ( Finite Element Analysis Tool ). The methodology used follows the classical steps used in computational methods, in- cluding formulation, implementation, verification and validation of results. Such verification is accomplished through shell and solid three-dimensional finite element models developed in the ABAQUS commercial software. The validation is performed by means of comparison with the experimental results found in literature. Regarding the evaluation of cross-sectional properties, one can observe a good agreement between the laminated beam theories adopted in this work and numerical and experimental results for all composite layups and load conditions conside- red. In the case of space frame elements, a good agreement is obtained between the results of finite elements proposed in this work and the analytical and computational results available in the literature. It is also observed that the element based on the Lagrangian formulation is more efficient than the element based on the linear theory regarding the ability to provide a satisfatory response with a less refined mesh

No presente trabalho são avaliados os principais parâmetros de projeto na análise da estabilidade global de estruturas em concreto pré-moldado de múltiplos pavimentos, considerando a não-linearidade física (NLF), não-linearidade geométrica (NLG) e ligações semi-rígidas. Os sistemas estruturais analisados são constituídos por pórticos com ligação viga-pilar semi-rígida e pilares engastados na fundação. As ligações viga-pilar têm sua tipologia definida por chumbadores retos e capa de concreto moldado in loco com armadura de continuidade. A NLF é avaliada segundo a construção dos diagramas M x N x 1/r, onde foi considerado de forma consistente o efeito da força normal, da armadura passiva, da armadura ativa e a reologia do concreto. São propostas funções e definidos coeficientes redutores de rigidez que foram comparados com valores apontados em expressões normativas que contemplam de forma aproximada a NLF. A NLG é avaliada com o auxílio do programa ANSYS e na forma aproximada segundo o coeficiente z e o método da carga lateral fictícia (P-). Apresentam-se modelos analíticos de caracterização de rigidez e resistência de ligações viga-pilar ao momento fletor negativo e positivo. Na análise numérica é feito um estudo de caso completo de uma estrutura típica de múltiplos pavimentos em concreto pré-moldado com o auxílio do programa ANSYS, avaliando-se diferentes formas de consideração da NLF e da NLG. É analisada a distribuição de esforços na estrutura frente às combinações de ações utilizadas e os modelos de comportamento das ligações, para algumas variações de geometria e carregamento. Dentre as conclusões, pode ser citado que os coeficientes redutores obtidos segundo o diagrama M x N x 1/r divergem das indicações normativas para consideração simplificada de NLF. Quanto à consideração simplificada de análise NLG, o coeficiente z apresentou resultados melhores na previsão dos esforços de segunda ordem em relação aos obtidos por 0,95.z. / In this research it was evaluated the main design parameters in the analysis of global stability of multi-storey precast concrete framed structures considering the physical nonlinearity (PNL), geometric nonlinearity (GNL) and semi-rigid connections behavior. The structural systems considered are those consisting of frames with beam-to-column semi-rigid connection and columns inset in the foundation. The beam-to-column connections are defined by two dowels and cast-in-place concrete cap with longitudinal reinforced bars. The PNL is evaluated according the M x N x 1/r diagrams, where are considered the effect of normal force, the reinforcing steel, prestressing steel and creep. Functions are proposed and reducing stiffness coefficients are defined, with were compared with codes that include the approximated PNL. The GNL is evaluated according to ANSYS software by no approximate analisys and its approximate according to the coefficient z and (P-) method. This research presents an analytical model characterization of stiffness and strength of beam-to-column connections to the negative and positive bending moment. In the numerical analysis is made a full case study of a typical multi-storey precast concrete framed structure with the aid of the ANSYS software, evaluating different forms of consideration of the PNL and the GNL. It analysed the bending moment distribution according the combinations of actions used and the behavioral models of connections, for some variations at geometry and loading. The reduction coefficients of flexural rigidity obtained according the M x N x 1/r diagrams differ from normative indications for simplified account of PNL. In GNL simplified analysis, the coefficient z showed better results in predicting the second-order effects with respect to those obtained by the coefficient 0,95. z.

Esse trabalho trata da análise computacional de problemas aeroelásticos. Os casos de interesse envolvem escoamentos compressíveis em torno de estruturas de cascas de materiais compósitos laminados. A solução do problema de interação fluido- estrutura é conduzida por meio do esquema particionado denominado improved serial staggered procedure (ISS), o que permite a utilização de algoritmos independentes para a análise de cada componente do sistema. Um elemento triangular plano com três nós destinado à análise de cascas de materiais compósitos laminados é obtido com base na formulação de dois elementos de alta performance desenvolvidos para a análise de membranas e de placas. A flexibilidade ao corte transversal, a ausência de travamento por razão de aspecto em problemas com flexão no plano da casca são características do elemento. A análise dinâmica não linear é realizada através do procedimento corrotacional de conservação aproximada de energia (AECCP) que foi implementado com base na formulação corrotacional independente de elemento (element independent CR formulation – EICR). Essa abordagem permite que problemas transientes com grandes rotações e forte não linearidade geométrica sejam integrados por longos períodos de tempo, mantendo a estabilidade, a precisão e a eficiência da solução. A simulação do escoamento é realizada através do método explícito de Taylor-Galerkin de dois passos utilizando elementos finitos tetraédricos em malhas não estruturadas para a discretização espacial. As equações governantes do escoamento são dadas segundo a formulação lagrangeana-euleriana arbitrária (ALE), permitindo que malhas móveis sejam empregadas na análise dos problemas aeroelásticos. Um método simples e eficiente é adotado para a movimentação da malha do fluido com base na distância do nó às superfícies de contorno do escoamento. A utilização de malhas não coincidentes na discretização do fluido e da estrutura é possibilitada pelo emprego do esquema de projeção nodal para a transferência de informações entre os meios na região de interface. A apresentação, implementação e verificação de cada componente da ferramenta de análise é realizada independentemente graças à natureza do esquema de solução particionado. Por fim são analisados problemas de interação fluido-estrutura, onde as respostas obtidas são comparadas com resultados experimentais e numéricos de outros autores. Também são propostos e analisados novos problemas envolvendo estruturas de compósitos laminados. / This work deals with computational analysis of aeroelastic problems, specifically those related to compressible flows around laminated composite shell structures. The fluid- structure problem is solved by a partitioned scheme called improved serial staggered procedure (ISS), allowing to use independent algorithms for the analysis of each system component. A 3-node triangular flat element for laminated composite shells is obtained based on the formulation of two high performance plate and membrane elements. Shear flexibility and no aspect ratio locking on in-plane bending problems are some of the element characteristics. The nonlinear dynamic analysis is performed using the approximately energy conserving corotational procedure (AECCP), which is implemented based on the element independent corotational formulation (EICR). This approach enables long time transient problems with strong nonlinearities and large rotations to be efficiently solved while keeping stability and accuracy. The flow simulation is performed by the explicit two step Taylor-Galerkin method employing tetrahedral finite elements in unstructured meshes for the space discrezation. The flow governing equations are given in arbitrary lagrangian-eulerian formulation, enabling aeroelastic simulations with moving meshes. A simple and efficient method is adopted to move the fluid mesh based on the distance of each node to the boundary surfaces. Transference of information between non-matching fluid and structure meshes are performed by the node-projection scheme. Each component of the aeroelastic analysis tool is independently presented, implemented and verified due to the partitioned nature of the scheme adopted for the solution of the coupled system. Some aeroelastic problems are analyzed and results are compared to other experimental and numerical works. New problems of aeroelastic analysis of laminated composite structures are proposed and solved.

A análise de compósitos laminados apresenta grandes desafios, pois, diferentemente dos materiais isotrópicos homogêneos, os compósitos laminados são constituídos de materiais heterogêneos e anisotrópicos. Além disso, as distribuições de tensões interlaminares obtidas com as formulações convencionais são descontínuas e imprecisas. Sua melhoria, portanto, é imprescindível para buscar e modelar critérios de falha relacionados às estruturas formadas por compósitos laminados. Diante disso, este trabalho se concentrou no desenvolvimento e implementação computacional de um elemento finito posicional de pórtico plano laminado cuja cinemática é descrita ao longo da espessura do laminado de acordo com a teoria Layerwise. A formulação do elemento considera a não linearidade geométrica, originada pela ocorrência de grandes deslocamentos e rotações, e admite deformações moderadas, em função da lei constitutiva de Saint-Venant-Kirchhoff. O desenvolvimento deste trabalho se iniciou com uma preparação teórica sobre mecânica dos sólidos deformáveis e métodos numéricos para que fossem adquiridos os subsídios teóricos necessários ao desenvolvimento de códigos computacionais, à interpretação dos resultados e à tomada de decisões quando das análises numéricas. A formulação desenvolvida é Lagrangiana total com emprego do método dos elementos finitos baseado em posições. Inicialmente o elemento finito posicional de pórtico plano homogêneo é proposto, uma vez que sua cinemática possibilita uma expansão natural para o caso laminado. Os graus de liberdade são compostos por posições nodais e por vetores generalizados que representam o giro e a variação na altura da seção transversal. A eficiência do elemento é constatada através de análises realizadas em problemas de pórtico sujeitos a grandes deslocamentos e rotações. Os resultados obtidos apresentaram excelente concordância com soluções numéricas e analíticas disponíveis na literatura. Uma expansão natural da cinemática é empregada na formulação do elemento laminado. Os graus de liberdade do elemento são as posições nodais e as componentes de vetores generalizados associados às seções transversais de cada lâmina. Dessa forma, as lâminas têm liberdade para variação de espessura e giro independente das demais, mas com as posições compatibilizadas nas interfaces. Os resultados de análises numéricas realizadas em vários exemplos demonstram a eficiência da formulação proposta, pois as distribuições de deslocamentos e tensões ao longo da espessura do laminado apresentaram excelente concordância com as obtidas a partir de análises numéricas utilizando um elemento finito bidimensional em uma discretização bastante refinada. Os exemplos analisados contemplam problemas com seção laminada fina ou espessa. / The analysis of laminated composites presents challenges because, unlike homogeneous isotropic materials, the laminated composites are made up of heterogeneous and anisotropic materials. Moreover, the distribution of interlaminar stresses obtained with conventional formulations are discontinuous and inaccurate. His improvement is therefore essential to check and modeling failure criteria related to structures formed by laminates. Thus, this work focused on developing and computational implementation of a positional finite element of laminated plane frame whose kinematics is described throughout the thickness of the laminate according to Layerwise theory. The formulation element considers the geometric nonlinearity, caused by the occurrence of large displacements and rotations, and admits moderate deformation, in the constitutive law function of Saint-Venant-Kirchhoff. The development of this work began with a theoretical preparation on mechanics of deformable solids and numerical methods for the acquired of the theoretical support needed for the development of computational codes, interpretation of results and decision-making when of the numerical analyzes. The developed formulation is total Lagrangian with use of the finite element method based on positions. Initially the positional finite element of homogeneous plane frame is proposed, since their kinematic enables a natural expansion for the laminate case. The degrees of freedom are composed of nodal positions and generalized vectors representing the spin and the variation in the height of the cross section. The efficiency of the element is verified through analyzes performed in frame problems subject to large displacements and rotations. The results showed excellent agreement with numerical and analytical solutions available in the literature. A natural expansion of the kinematics is used in the formulation of the laminate element. The degrees of freedom of the element are the nodal positions and components of the generalized vectors associated to cross-sections of each lamina. Thus, the laminas are free for the thickness variation and for independent spin, but with the positions matched in the interfaces. The results of numerical analysis performed in various examples show the effectiveness of the proposed formulation, since the distributions of displacements and stresses through the thickness of the laminate agreed well with those obtained from numerical analysis using a discretization with two-dimensional finite elements in a very refined. The examples discussed include problems with thin or thick laminated section.

The subject of this thesis is analysis of the plastic bearability of plates. The first part is dedicated to a description of pertinent theories. It is primarily a theory of material nonlinearity and plasticity, including a description of plasticity models, which are commonly used in engineering analysis. Next is the hinted theory of geometric nonlinearity and discusses the key principles of limit plastic bearability of plates. This thesis includes ultimate strength calculations of plates of selected shapes, their solution by both, a manual calculation and using computational software. The obtained values are subsequently evaluated. The appendices contain the procedures for modelling in the computational programs RFEM and ANSYS.

Thesis deals with flexural-lateral buckling of asymmetric cross sections, especially assessment according to ČSN EN 1993-1-1. Models are based on bar, shell and solid finite element method. Assessment of software Scia Engineer and RFEM are analysed in bar model. Plate and volume model uses shell and solid elements and initial imperfections. Computing is geometrically nonlinear. Correspondence of bar and plate model is studied. Optimal amplitude is suggested. Simple beam loaded with continuous line loading is taken as an example.

[pt] O design otimizado de estruturas considerando não-linearidades tem sido amplamente pesquisado nas décadas recentes. A análise de elementos finitos aplicada à otimização topológica é prejudicada pela deformação excessiva de elementos de baixa densidade sob alta compressão, o que impede o processo de encontrar uma solução ótima. Dois métodos, o esquema Interpolação de Energia e a técnica de Hiperelasticidade Aditiva, são implementados para superar essa dificuldade no problema de minimização da flexibilidade, e modelos de materiais hiperelásticos são usados para investigar suas influências na topologia otimizada. O Método das Assíntotas Móveis é usado para atualizar as variáves de projeto cujas sensibilidades foram calculadas pelo método adjunto. A equação de estado é resolvida através do método de Newton-Raphson com um incremento de carga ajustável para reduzir o custo computacional. Resultados de dois problemas de referência são comparado com aqueles já estabelecidos na literatura. O uso de diferentes modelos hiperelásticos apresentou pouca influência no design final da estrutura. O método de Interpolação de Energia foi capaz de convergir para cargas muito maiores que o método padrão, enquanto a Hiperelasticidade Aditiva apresentou dificuldades de convergência em estado plano de deformação. / [en] The optimized design of structures considering nonlinearities has been widely researched in the recent decades. The finite element analysis applied to topology optimization is jeopardized by the excessive deformation of low-density elements under high compression, which hinders the process of finding an optimal solution. Two methods, the Energy Interpolation scheme and the Additive Hyperelasticity technique, are implemented to overcome this difficulty in the minimum compliance problem, and hyperelastic material models are used to investigate their influence on the optimized topology. The Method of Moving Asymptotes is used to update the design variables whose sensitivities were calculated from the adjoint method. The state equation is solved through the Newton-Raphson method with an adjusting load step to reduce computational cost. Results for two benchmark problems are compared with those already established in the literature. The use of different hyperelastic models presented little influence on the final design of the structure. The Energy Interpolation method was able to converge for much higher loads than the default method, while the Additive Hyperelasticity presented convergence difficulties in plane strain.

[pt] O principal objetivo deste trabalho é o desenvolvimento de um programa computacional para a analise não-linear estática e dinâmica de pórticos planos com ligações flexíveis (semi-rígidas). Inicialmente é apresentada a metodologia de solução não-linear e as formulações dos elementos finitos adotados na base computacional implementada. Em seguida, são estudados vários exemplos de sistemas estruturais estáticos com caminhos de equilíbrio fortemente não-lineares com a finalidade de testar os programas implementados. Então é apresentada a formulação do problema dinâmico com a definição das equações diferenciais ordinárias de movimento e as expressões das matrizes de massa e amortecimento. A solução desse sistema de equações diferenciais ordinárias é obtida por métodos de integração numérica implícitos ou explícitos. Alguns destes métodos são apresentados neste trabalho e incorporados ao programa computacional em conjunto com estratégias adaptativas de incremento automático do intervalo de tempo de integração (delta)t. Por fim, o sistema computacional desenvolvido é utilizado na modelagem e obtenção da resposta estrutural estática e dinâmica de alguns sistemas estruturais planos com comportamento eminentemente não-linear. Através destes resultados são analisados alguns fenômenos importantes de instabilidade estática e dinâmica, bem como possíveis mecanismos de colapso e a influência de parâmetros físicos e geométricos no comportamento estrutural. / [en] The main objective of this thesis is to develop a numerical methodology for the nonlinear static and dynamics analysis of plane frames with semi-rigid connections. Initially, the formulations of the adopted finite elements are presented and implemented together with numerical methodologies for the solution of the non-linear equilibrium equations. Then, some examples of strongly nonlinear structural systems under static loads are studied to check the methodology. Subsequently, the ordinary differential equations of motion are derived and the corresponding damping and mass matrices are presented. The solution of this system of ordinary differential equations is obtained by implicit or explicit numerical integration methods. Some of these methods are presented in this work and incorporated into the computational program together with adaptive strategies for the automatic increment of the time step (delta)t. Finally, the computational system here developed is used to study the static and dynamic response of some plain structural systems with an inherent nonlinear behavior. A detailed parametric study is carried out to identify the influence of physical and geometric parameters on the structural behavior. This enables the analysis of some important static and dynamic instability phenomena and identification of possible mechanisms of collapse.