O objetivo principal deste trabalho é o desenvolvimento de uma formulação e sua implementação computacional para se analisar, via Método dos Elementos Finitos (MEF),o comportamento dinâmico não linear geométrico de sólidos bidimensionais. Trata-se o comportamento geometricamente não linear através de uma formulação posicional classificada como lagrangeana total com cinemática exata. No estudo do comportamento dinâmico utiliza-se um algoritmo de integração temporal baseado na família de integradores temporais de Newmark. Para a consideração do impacto adota-se uma técnica que utiliza como integrador temporal o algoritmo de Newmark, modificado de forma a garantir sua estabilização, e limita-se a posição de cada nó da estrutura que por ventura sofra impacto. O código computacional desenvolvido é validado através de exemplos tradicionais da literatura científica. Analisam-se exemplos com comportamento apenas não linear geométrico e não linear geométrico dinâmico com ou sem impacto / The main goal of this work is the development of a formulation and its computational implementation, based on the finite element method (FEM), to analyze the dynamic geometrically nonlinear behavior of bidimensional solids. The geometrically nonlinear behavior is treated with a positional formulation classified as total Lagrangean with exact kinematics. In the study of the dynamic behavior, a time integration algorithm based on the family of time integrators of Newmark is applied. In order to consider the impact, a technique based on the time integrator of Newmark, modified to assure its stabilization, is used. This technique limits the position of each node that suffers impact. The developed computational code is validated through benchmarks of scientific literature. Examples with static geometrically nonlinear and dynamic geometrically nonlinear behavior, with or without impact, are analyzed

Tese de mestrado integrado. Engenharia Civil (Estruturas). Universidade do Porto. Faculdade de Engenharia. 2010

O objetivo principal deste trabalho é o desenvolvimento de uma formulação e sua implementação computacional para se analisar, via Método dos Elementos Finitos (MEF),o comportamento dinâmico não linear geométrico de sólidos bidimensionais. Trata-se o comportamento geometricamente não linear através de uma formulação posicional classificada como lagrangeana total com cinemática exata. No estudo do comportamento dinâmico utiliza-se um algoritmo de integração temporal baseado na família de integradores temporais de Newmark. Para a consideração do impacto adota-se uma técnica que utiliza como integrador temporal o algoritmo de Newmark, modificado de forma a garantir sua estabilização, e limita-se a posição de cada nó da estrutura que por ventura sofra impacto. O código computacional desenvolvido é validado através de exemplos tradicionais da literatura científica. Analisam-se exemplos com comportamento apenas não linear geométrico e não linear geométrico dinâmico com ou sem impacto / The main goal of this work is the development of a formulation and its computational implementation, based on the finite element method (FEM), to analyze the dynamic geometrically nonlinear behavior of bidimensional solids. The geometrically nonlinear behavior is treated with a positional formulation classified as total Lagrangean with exact kinematics. In the study of the dynamic behavior, a time integration algorithm based on the family of time integrators of Newmark is applied. In order to consider the impact, a technique based on the time integrator of Newmark, modified to assure its stabilization, is used. This technique limits the position of each node that suffers impact. The developed computational code is validated through benchmarks of scientific literature. Examples with static geometrically nonlinear and dynamic geometrically nonlinear behavior, with or without impact, are analyzed

Nesta tese, apresenta-se o método dos elementos finitos posicional descrito em um referencial Lagrangiano total dedicado à análise de instabilidade de estruturas tridimensionais. Três tipos de elementos finitos são implementados e testados, a saber: os elementos de barra simples, casca e barra geral. A análise de instabilidade para o elemento de barra simples é efetuada determinando-se os pontos críticos ao longo da trajetória de equilíbrio em grandes deslocamentos. Para se determinar essas trajetórias são utilizados os algoritmos de Newton-Raphson e arc-length. Este tipo de análise é particularmente importante na definição de estruturas multi-estáveis de uso crescente na indústria mecânica e aeroespacial. Para o estudo da instabilidade empregando-se os elementos finitos de casca e barra geral realizam-se as análises para pequenos níveis de carga e deslocamentos por meio do cálculo dos autovalores e autovetores da matriz de rigidez da estrutura. Avaliam-se também as trajetórias de equilíbrio em grandes deslocamentos considerando-se pequenas imperfeições na geometria dos elementos estruturais. Quando os elementos de casca são utilizados na modelagem de perfis estruturais esbeltos surgem naturalmente modos de falha locais associados à mudança de forma da seção transversal. Com a finalidade de inserir essas mobilidades no elemento de barra geral propõem-se uma metodologia que considera os aprimoramentos na cinemática da barra. Esses aprimoramentos são tratados como parâmetros nodais generalizados e estão associadas a intensidade da mudança de forma de seção transversal, incluindo os modos de empenamento. Descreve-se originalmente uma metodologia de decomposição da matriz Hessiana usada para o cálculo dos valores e vetores próprios em pequenos deslocamentos. Essa metodologia possui importância adicional pois é utilizada na preparação e avaliação do parâmetro de carga em cinemáticas alternativas da formulação posicional. Utiliza-se o algoritmo de Lanczos na determinação das cargas e modos de falha realizando-se chamadas a biblioteca ARPACK. Os algoritmos são testados em exemplos modelados com os elementos finitos propostos. Próximo aos pontos críticos realiza-se a separação da matriz Hessiana procurando-se possíveis modos de colapso da estrutura. Além dos modos de falha globais é possível se identificar os modos de falha locais e distorcionais. O equilíbrio do sistema mecânico é garantido pelo princípio da estacionariedade da energia potencial total. Nas análises com os elementos de casca e barra geral, a solução do sistema não-linear é obtida empregando-se o método incremental iterativo de Newton-Raphson. Os aprimoramentos sugeridos nesta pesquisa são acoplados ao código computacional utilizado pelo grupo de mecânica computacional do departamento de engenharia de estruturas, onde diversas funcionalidades estão disponíveis, como análise dinâmica e não-linearidade material. Exemplos selecionados são apresentados ao longo da tese para demonstrar a eficiência dos elementos propostos e a aplicabilidade da técnica. Por fim, são realizadas comparações com estratégia de solução já consagradas, como por exemplo: o método das faixas finitas e a teoria generalizadas de vigas. Os resultados obtidos justificam as contribuições originais da presente pesquisa destacando-se a contribuição da formulação posicional ao estudo da instabilidade das estruturas. / This thesis presents the positional finite element method in a total Lagrangian framework dedicated to instability analysis of the three-dimensional structures. Three types of finite elements are implemented and tested, namely: truss, shells and frames. The instability analysis for truss element is computed using equilibrium path in large displacements. The critical points are computed using Newton-Raphson and arc-length algorithm. This analysis is particularly important in the definition of multi-stable and large displacements structures widely used in mechanical and aerospace industry. For shell and frame geometrically non-linear finite elements, the instability phenomenon is studied from the eigenvalues and eigenvectors analysis for small levels of loads and displacements. It is also evaluate the equilibrium trajectories for large displacements, considering small imperfections in the geometry of the structure. When using the shell elements to model the frames structures local failure modes associated with changing of the cross section shape arise. In order to consider the mobility in frame element new improvements are propose in the kinematic. These improvements are treated as generalized nodal parameters and are associated with the intensity of the cross-sectional change, including warping. The originally methodology of decomposition of the Hessian matrix are described and used for calculating eigenvalues and eigenvectors of the stiffness matrix. This methodology has additional importance because it is used in the preparation and evaluation of load parameter in kinematic alternatives of the positional formulation. The Lanczos algorithm is used to determining the loads and failure modes, through calls to ARPACK library for calculating eigenvalues and eigenvectors. The algorithms are tested on the examples modeled by proposed finite elements. Near the critical point takes place the separation of the Hessian matrix for possible identification of the failure modes. In addition to global failure methods, local and distortion failure are captured by this methodology. The balance of the mechanical system is guaranteed by the stationarity of the total potential energy principle. In the analysis using shells and frames elements the solution of the nonlinear system is calculated using the iterative incremental Newton-Raphson method. The improvements suggested in this research are coupled to the computer code used by computational mechanics group of the structures engineering department, where several features are available like dynamic and plasticity analysis. Selected examples are presented throughout the thesis to demonstrate the efficiency of the proposed elements and applicability of the technique. Finally, comparisons are carried out with already established solving strategy such as the finite strip methods and the generalized beam theory. The results justified the original contributions of this research to study of unstable structures.

Nesta tese, apresenta-se o método dos elementos finitos posicional descrito em um referencial Lagrangiano total dedicado à análise de instabilidade de estruturas tridimensionais. Três tipos de elementos finitos são implementados e testados, a saber: os elementos de barra simples, casca e barra geral. A análise de instabilidade para o elemento de barra simples é efetuada determinando-se os pontos críticos ao longo da trajetória de equilíbrio em grandes deslocamentos. Para se determinar essas trajetórias são utilizados os algoritmos de Newton-Raphson e arc-length. Este tipo de análise é particularmente importante na definição de estruturas multi-estáveis de uso crescente na indústria mecânica e aeroespacial. Para o estudo da instabilidade empregando-se os elementos finitos de casca e barra geral realizam-se as análises para pequenos níveis de carga e deslocamentos por meio do cálculo dos autovalores e autovetores da matriz de rigidez da estrutura. Avaliam-se também as trajetórias de equilíbrio em grandes deslocamentos considerando-se pequenas imperfeições na geometria dos elementos estruturais. Quando os elementos de casca são utilizados na modelagem de perfis estruturais esbeltos surgem naturalmente modos de falha locais associados à mudança de forma da seção transversal. Com a finalidade de inserir essas mobilidades no elemento de barra geral propõem-se uma metodologia que considera os aprimoramentos na cinemática da barra. Esses aprimoramentos são tratados como parâmetros nodais generalizados e estão associadas a intensidade da mudança de forma de seção transversal, incluindo os modos de empenamento. Descreve-se originalmente uma metodologia de decomposição da matriz Hessiana usada para o cálculo dos valores e vetores próprios em pequenos deslocamentos. Essa metodologia possui importância adicional pois é utilizada na preparação e avaliação do parâmetro de carga em cinemáticas alternativas da formulação posicional. Utiliza-se o algoritmo de Lanczos na determinação das cargas e modos de falha realizando-se chamadas a biblioteca ARPACK. Os algoritmos são testados em exemplos modelados com os elementos finitos propostos. Próximo aos pontos críticos realiza-se a separação da matriz Hessiana procurando-se possíveis modos de colapso da estrutura. Além dos modos de falha globais é possível se identificar os modos de falha locais e distorcionais. O equilíbrio do sistema mecânico é garantido pelo princípio da estacionariedade da energia potencial total. Nas análises com os elementos de casca e barra geral, a solução do sistema não-linear é obtida empregando-se o método incremental iterativo de Newton-Raphson. Os aprimoramentos sugeridos nesta pesquisa são acoplados ao código computacional utilizado pelo grupo de mecânica computacional do departamento de engenharia de estruturas, onde diversas funcionalidades estão disponíveis, como análise dinâmica e não-linearidade material. Exemplos selecionados são apresentados ao longo da tese para demonstrar a eficiência dos elementos propostos e a aplicabilidade da técnica. Por fim, são realizadas comparações com estratégia de solução já consagradas, como por exemplo: o método das faixas finitas e a teoria generalizadas de vigas. Os resultados obtidos justificam as contribuições originais da presente pesquisa destacando-se a contribuição da formulação posicional ao estudo da instabilidade das estruturas. / This thesis presents the positional finite element method in a total Lagrangian framework dedicated to instability analysis of the three-dimensional structures. Three types of finite elements are implemented and tested, namely: truss, shells and frames. The instability analysis for truss element is computed using equilibrium path in large displacements. The critical points are computed using Newton-Raphson and arc-length algorithm. This analysis is particularly important in the definition of multi-stable and large displacements structures widely used in mechanical and aerospace industry. For shell and frame geometrically non-linear finite elements, the instability phenomenon is studied from the eigenvalues and eigenvectors analysis for small levels of loads and displacements. It is also evaluate the equilibrium trajectories for large displacements, considering small imperfections in the geometry of the structure. When using the shell elements to model the frames structures local failure modes associated with changing of the cross section shape arise. In order to consider the mobility in frame element new improvements are propose in the kinematic. These improvements are treated as generalized nodal parameters and are associated with the intensity of the cross-sectional change, including warping. The originally methodology of decomposition of the Hessian matrix are described and used for calculating eigenvalues and eigenvectors of the stiffness matrix. This methodology has additional importance because it is used in the preparation and evaluation of load parameter in kinematic alternatives of the positional formulation. The Lanczos algorithm is used to determining the loads and failure modes, through calls to ARPACK library for calculating eigenvalues and eigenvectors. The algorithms are tested on the examples modeled by proposed finite elements. Near the critical point takes place the separation of the Hessian matrix for possible identification of the failure modes. In addition to global failure methods, local and distortion failure are captured by this methodology. The balance of the mechanical system is guaranteed by the stationarity of the total potential energy principle. In the analysis using shells and frames elements the solution of the nonlinear system is calculated using the iterative incremental Newton-Raphson method. The improvements suggested in this research are coupled to the computer code used by computational mechanics group of the structures engineering department, where several features are available like dynamic and plasticity analysis. Selected examples are presented throughout the thesis to demonstrate the efficiency of the proposed elements and applicability of the technique. Finally, comparisons are carried out with already established solving strategy such as the finite strip methods and the generalized beam theory. The results justified the original contributions of this research to study of unstable structures.

Made available in DSpace on 2014-07-29T15:03:39Z (GMT). No. of bitstreams: 1 lilian pureza.pdf: 2774999 bytes, checksum: 2a074d04ee02c7e1c87fdbe8c2c68ef6 (MD5) Previous issue date: 2006-10-20 / The aim of this work is to present a formulation and corresponding computational implementation for sizing optimization of plane frames and cable-stayed columns considering geometric non liner behavior. The structural analysis is based on the finite element method using the updated lagrangian approach for plane frame and cable elements, which are represented by plane truss elements. The non linear system is solved by the Newton-Raphson method coupled to load increment strategies such as the arch length method and the generalized displacement parameter method, which allow the algorithm to transpose any critical point that happen to appear along the equilibrium path. In the optimization process the design variables are the heights of the crosssection of the frame elements, the objective function represents the volume of the structure and the constraints impose limits to displacements and critical load. Lateral constraints impose limits to the design variables. The finite difference method is used in the sensitivity analysis of the displacement and critical load constraints. The optimization process is carried out using three different optimization strategies: the sequential quadratic programming algorithm; the interior points algorithm; and the branch and bound method. Some numerical experiments are carried out so as to test the analysis and the sensitivity strategies. Numerical experiments are presented to show the validity of the implementation presented in this dissertation. / O objetivo deste trabalho é a otimização de dimensões de pórticos planos e de colunas estaiadas planas pela minimização do volume da estrutura, considerando os efeitos da não-linearidade geométrica em seu comportamento. A formulação utiliza, para análise das estruturas, elementos finitos de pórtico e de treliça planos e referencial lagrangeano atualizado. O método de Newton-Raphson foi utilizado como estratégia para solução do sistema de equações não lineares. Foram acopladas estratégias especiais para ultrapassagem de pontos críticos que possam existir ao longo da trajetória de equilíbrio, tais como o comprimento de arco cilíndrico e o controle dos deslocamentos generalizados. Na otimização, as variáveis de projeto são as alturas das seções transversais dos elementos, a função objetivo é o volume do material e as restrições dizem respeito a limitações impostas a deslocamentos e à carga limite, além de limitações impostas aos valores das variáveis. A sensibilidade da função objetivo foi obtida por diferenciação direta e a sensibilidade das restrições pelo método das diferenças finitas. Foram utilizados o algoritmo de programação quadrática seqüencial, PQS, o algoritmo de pontos interiores, PI, e o algoritmo de Branch and Bound, B&B. São apresentados exemplos de validação das estratégias de análise não linear e da análise de sensibilidade, além dos exemplos de validação da formulação empregada para a otimização resolvidos pelos métodos implementados.

A interação entre fluidos e estruturas caracteriza um problema multi-físico não linear e está presente numa grande variedade de áreas da engenharia. Este trabalho apresenta o desenvolvi mento de ferramentas computacionais com base no Método dos Elementos Finitos (MEF) para a análise de interação fluido-estrutura (IFE) considerando escoamentos com baixas velocidades. Dada a interdisciplinaridade do tema, se faz necessário o estudo em três diferentes assuntos: a dinâmica das estruturas computacional, a dinâmica dos fluidos computacional, e o problema de acoplamento. No caso da dinâmica das estruturas empregar-se um elemento finito que seja adequado para a simulação de problemas de IFE, que claramente demandam uma análise não linear geométrica, optando-se pelo emprego de uma formulação descrita em posições, a qual evita problemas relativos à aproximação de rotações finitas. Quanto à dinâmica dos fluidos computacional, é empregado um método estável e ao mesmo tempo sensível à movimentação da estrutura, utilizando a descrição Lagrangeana-Euleriana Arbitrária (ALE). Os casos considerados neste trabalho, assim como muitos dos problemas de engenharia, ocorrem com escoamentos em baixas velocidades, implicando na incompressibilidade do fluido, o que demanda, para um método estável, a utilização de elementos que atendam à condição de Ladyzhenskaya-Babuska-Brezzi (LBB). Além disso, é necessário também o emprego de métodos que consigam neutralizar as variações espúrias decorrentes da não-linearidade de possíveis escoamentos com convecção dominante e que surgem com a aplicação do processo clássico de Galerkin. Para superar esse problema, é aplicado o método Streamline-Upwind/Petrov-Galerkin (SUPG), que adiciona difusividade artificial na direção do escoamento, controlando a amplitude dos termos convectivos. No que se refere ao acoplamento fluido-casca, buscam-se modularidade e versatilidade adotando-se o modelo particionado. O modelo de acoplamento implementado garante ainda a utilização de malhas do fluido e da estrutura sem a necessidade de coincidência de nós. / Interaction between fluids and structures characterizes a nonlinear multi-physics problem presente in a wide range of engineering fields. This works presets the development of computational tools based on finite element method (FEM) for fluid-structure interaction (FSI) analysis considering low speed flows (incompressible), as a great part of the engineering problems. Given the topic multidisciplinary nature, it is necessary to study three different subjects: the computational structural dynamics, the computational fluid mechanics and the coupling problem. Regarding structural mechanics, we seek to employ a finite element adequate to FSI simulation, what clearly demands a geometric nonlinear analysis. We chose to employ shell elements with formulation in terms of positions, which avoids problems related to finite rotations approximations. Concerning computational fluid dynamics, we employ a stable method, at same time sensible o structural movements, which is written in the arbitrary Lagrangian-Eulerian (ALE) description. The flow incompressibility demands, for a stable method, the use of elements according to the Ladyzhenskaya-Bbuska-Brezzi (LBB) condition. It is also necessary to employ methods able to neutralize the spurious variations that appears from convection dominated flows when applying the standard Galerking method. In order to overcome this problem, we apply the Streamline-Upwind/Petrov-Galerkin (SUPG) method, which adds artificial diffusivity to the streamline direction, controlling spurious variations. Considering the fluid-shell coupling, we seek modularity and versatility, adopting the partitioned model. The developed coupling model ensure the use of fluid and structure meshes with no need for matching nodes.

A utilização de materiais compósitos tornou-se uma alternativa importante em muitas aplicações dentro de diversas áreas da engenharia, pois seus constituintes podem agregar propriedades mecânicas, térmicas e acústicas ao compósito, garantindo eficiência e baixo custo. Com isso, faz-se necessário um maior conhecimento do comportamento mecânico desses materiais diante das solicitações, principalmente no que diz respeito aos campos de deslocamento, deformações e tensões. O presente trabalho tem por finalidade a análise, em nível macroscópico, de estruturas bidimensionais elásticas constituídas de materiais compósitos particulados, utilizando formulação desenvolvida no contexto do Grupo de Mecânica Computacional (GMEC), do Departamento de Engenharia de Estruturas (SET), da Escola de Engenharia de São Carlos (EESC), da Universidade de São Paulo (USP), no qual se insere a presente pesquisa. A formulação utilizada baseia-se no Método dos Elementos Finitos Posicional (MEFP) e foi desenvolvida em nível mesoscópico por tratar da interação entre matriz e partículas. Tal formulação possibilita a consideração da interação partícula-matriz sem a necessidade de coincidência entre as malhas da matriz e das partículas e sem o aumento do número de graus de liberdade dos problemas, admitindo-se aderência perfeita entre as fases. A formulação considera material isotrópico e comportamento não-linear geométrico das fases. A aplicação da formulação foi aqui proposta com o intuito de avaliar a influência da geometria, tamanho, fração volumétrica, distribuição e propriedades mecânicas das partículas adotadas, no comportamento global da estrutura em nível macroscópico. Foram desenvolvidos e apresentados exemplos de aplicação, com comparação dos resultados numéricos das análises com resultados de ensaios experimentais encontrados na literatura, bem como com resultados de modelos matemáticos de homogeneização e modelos numéricos propostos por outros autores, que utilizaram o método dos elementos finitos e técnicas de homogeneização assintótica. / The use of composite materials has become an important alternative in many applications in different areas of engineering, because their constituents can add mechanical, thermal and acoustic properties to the composite, ensuring efficiency and low cost. Thus, it is necessary a better understanding of the mechanical behavior of these materials, mainly regarding displacement, stress and strain fields. This study aims to analyze, in macroscopic scale, two-dimensional elastic structures made of particulate composite materials, using formulation developed in the context of the Grupo de Mecânica Computacional (GMEC), of Departamento de Engenharia de Estruturas (SET), of Escola de Engenharia de São Carlos (EESC), of Universidade de São Paulo (USP). The formulation is based on the Positional Finite Element Method and was developed in mesoscopic level, considering the matrix-particles interaction and neglecting the interface, by means of kinematic relations used to ensure adherence of the particles to the matrix without introducing new degrees of freedom in the problem. The formulation considers isotropic material and geometric non-linear behavior of the composite phases. The application of the formulation was proposed in this work in order to evaluate the influence of geometry, size, volume fraction, distribution and mechanical properties of the particles adopted in the global behavior of the structure in macroscopic level. Numerical examples were developed and presented in order to compare the numerical results of the analysis with results obtained in experimental studies found in the literature, as well as results of mathematical models and numerical models using finite element method and the asymptotic homogenization technique.

Neste estudo se propõe o desenvolvimento de um modelo numérico para a ligação deslizante entre elementos sólidos bidimensionais, aplicável à simulação de sistemas deslizantes de isolação de base para estruturas. A formulação implementada é baseada no Método dos Elementos Finitos Posicional (MEFP) para análise dinâmica não linear geométrica de estruturas escrita na forma Lagrangeana total. Elementos triangulares planos e isoparamétricos de aproximação cúbica com matriz de massa completa são utilizados principalmente na elaboração da parte sólida dos dispositivos de ligação entre estruturas reticuladas e a base móvel. Esses elementos também poderm ser utilizados na modelagem da estrutura em si, porém, para esse fim, elementos finitos isoparamétricos de barra geral com massa distribuída por unidade de comprimento foram implementados. As equações de movimento são integradas no tempo aplicando o método de Newmark e o problema de deslizamento é resolvido com o algoritmo baseado na técnica dos multiplicadores de Lagrange, onde a restrição das posições de um nó escravo é feita em relação a uma sequência de superfícies mestres. Elementos de barra geral foram usados para simular as superfícies mestres de contato, o que aumenta as possibilidades de aplicações, incluindo mecanismos compostos apenas por barras gerais. Analisam-se exemplos disponíveis na literatura para a validação da formulação proposta e propõem-se aplicações diversas na engenharia das estruturas. / This study proposes the development of a numerical model for the sliding joint between two-dimensional solid elements, applicable to the simulation of sliding base isolation systems. The implemented formulation is based on the Positional Finite Element Method (PFEM) for geometrical nonlinear dynamic analysis of structures written in the total Lagrangian form. Plane and isoparametric triangular cubic approximation elements with full mass matrix are mainly used in the elaboration of the solid part of the devices of joints between reticulated structures and mobile base. These elements can also be used in the modeling of the structure itself, however, for that purpose, isoparametric elements of general bar with mass distributed per unit of length were implemented. The motion equations are integrated in time by applying the Newmark method and the sliding problem is solved with the algorithm based on the technique of Lagrange multipliers, where the constraint of the positions of a slave node is made in relation to a sequence of master surfaces. General bar elements were used to simulate the master contact surface, which increases the possibilities of applications, including mechanisms composed only of general bars. Analyze examples available in the literature for the validation of the proposed formulation and proposed diverse applications in the engineering of the structures.

O objetivo central deste trabalho é implementar modelos constitutivos hiperelásticos não lineares em um código computacional que faz análise não linear geométrica de cascas. São necessários, para este propósito, conceitos sobre álgebras linear e tensorial, cinemática, deformação, tensão, balanços, princípios variacionais, métodos numéricos e hiperelasticidade. Tal programa usa a formulação Lagrangiana posicional, o método dos elementos finitos, o princípio dos trabalhos virtuais e o método iterativo de Newton-Raphson para solução das equações não lineares. O elemento finito de casca possui dez nós, sete parâmetros por nó e variação linear da deformação ao longo da espessura. Para dedução dos novos modelos usou-se a decomposição multiplicativa do gradiente da função mudança de configuração, o tensor deformação de Green-Lagrange e o tensor da tensão de Piola-Kirchhoff de segunda espécie. O código desenvolvido foi usado em simulações de diversos exemplos e apresentou boa precisão na análise mecânica de polímeros naturais altamente deformáveis. A ocorrência do fenômeno travamento não se manifestou nas análises realizadas. A presente pesquisa confirmou outros trabalhos, reforçou a necessidade de se usar modelos hiperelásticos não lineares para simular o comportamento mecânico de polímeros naturais e apresentou resultados condizentes com dados experimentais existentes na literatura científica e às respectivas soluções analíticas. / The main objective of this work is to implement nonlinear hyperelastic constitutive models in a computational code of geometrically nonlinear analysis of shells. For this purpose, concepts of linear and tensor algebras, kinematics, strain, stress, balances, variational principles, numerical methods and hyperelasticity are necessary. Such program uses the positional Lagrangian formulation, the finite element method, the principle of virtual work and the iterative method of Newton-Raphson for the solution of the nonlinear equations. The shell finite element has ten nodes, seven parameters per node and presents linear variation of the strain along the thickness. To achieve the new constitutive models the multiplicative decomposition of the deformation gradient, the Green-Lagrange strain tensor and the second Piola-Kirchhoff stress tensor are used. The developed code is tested for simulations of various examples and presents good accuracy in the mechanical analysis of highly deformable natural rubber. The locking phenomena didn\'t appear in the proposed analysis. The present research confirms other works, corroborates the need of using nonlinear hyperelastic models to simulate the mechanical behavior of natural rubber and presents suitable results when compared to existent experimental data of the scientific literature and to the respective analytical solutions.