[pt] O presente trabalho tem por objetivo desenvolver uma formulação e certas estratégias que permitam a análise da perda de estabilidade de colunas esbeltas submetidas a carregamento axial periódico, fenômeno este conhecido como ressonância paramétrica. Uma excitação é dita paramétrica quando aparece nas equações de movimento do sistema na forma de coeficientes variáveis com o tempo - geralmente periódicos - e não como uma não homogeneidade. A coluna é descrita pela formulação clássica de Navier. O presente trabalho trata a coluna considerando-se um e três graus de liberdade com ou sem não-linearidades. As equações de movimento são obtidas utilizando-se o princípio de Hamilton através do método de Ritz. A equação linear (equação de Mathieu) e a equação de Duffing com pequeno amortecimento, são resolvidas de forma aproximada pelo método das múltiplas escalas, revelando a possibilidade de instabilização da posição de equilíbrio em diversas regiões do espaço definido pelos parâmetros de controle. A mesma conclusão é mostrada utilizando-se procedimentos computacionais para a resolução dos sistemas de equações lineares e nãolineares, com ou sem imperfeição geométrica inicial, podendo-se obter assim, a resposta do sistema, planos fase, seções de Poincaré e diagramas de bifurcação. Mostra-se a partir dos resultados numéricos, que a coluna submetida a cargas axiais harmônicas, pode tanto apresentar soluções com o mesmo período da força excitadora, quanto oscilações subarmônicas e superarmônicas de diversas ordens, além de movimentos caóticos. / [en] The main aim of the present work is to develop a formulation and some strategies for the instability analysis of slender columns under an axial harmonic force this phenomenon is known as parametric ressonance. An excitation is said to be parametric if it appears as timedependent - often periodic - coefficients in the equations governing the motion of the system,and not as an inhomogeneous term.The column is described by Navier classical formulation. The present work consider the column with one or three degrees of freedom with or without nonlinearities. The equations governing the motion are obtained by the Ritz method.The linear equation (Mathieu equation) and the Duffing equation with small damping are solved in an approximate way using multiple scales techniques, revealing the possibility of destabilizing the static equilibrium position in certain regions of the control space. A similar conclusion is obtained by employing numerical methods for the solution of linear and nonlinear equation systems with or without initial geometrical imperfections.This enables one to obtain time response, phase space, projections Poincaré sections and bifurcation diagrams. These numerical results show that the column with nonlinearities and loaded by a periodic longitudinal force can present various solutions with the same period as the forcing and subharmonic e superharmonic oscillations, as well as chaotic motions.

[pt] Estruturas offshore têm demandado, em função do aumento da profundidade da lâminha de água e da severidade do ambiente, análises de vibração cada vez mais confiáveis. Em face de oscilações com grandes deslocamentos, torna-se imprescindível uma análise não linear dessas estruturas. Métodos numéricos como os elementos finitos constituem-se numa tarefa computacionalmente custosa, uma vez que os acoplamentos modais tornam necessários modelos com muitos graus de liberdade. Isso dificulta as análises paramétricas e prolonga os ciclos de projeto para estruturas offshore. Uma alternativa a esses problemas é o uso de modelos de ordem reduzida. Os modos normais não lineares têm-se mostrado uma ferramenta eficiente na derivação de modelos de ordem reduzida para análises de vibrações não lineares. Isso ocorre porque um número menor de modos não lineares, em relação aos modelos com modos lineares, é necessário para se obter o mesmo nível de precisão num modelo reduzido. Esse trabalho utiliza modelos de ordem reduzida, obtidos por meio de análise modal não linear, para o estudo de vibração de modelos simplificados de estruturas offshore. Três exemplos de aplicação são utilizados: pêndulo invertido, torre articulada e plataforma spar. Além dos métodos baseado no procedimento de Galerkin e o assintótico, um procedimento numérico alternativo é proposto para obtenção dos modos, podendo ser utilizado para construção dos modos essencialmente não lineares. As vibrações livres e forçadas são estudadas. A estabilidade das soluções é analisada utilizando-se a teoria de Floquet, diagramas de bifurcação e de Mathieu e seções de Poincaré. As seções de Poincaré são também utilizadas para identificar a multiplicidade dos modos não lineares e a existência de multimodos. Os resultados são comparados com a solução obtida da integração numérica do sistema original de equações, mostrando uma boa precisão dos modelos reduzidos. / [en] The increasing water depth and the ocean adverse environment demand more accurate vibration analysis of offshore structures. Due to large amplitude oscillations, a nonlinear vibration analysis becomes necessary. Numerical methods such as finite element constitute a computationally expensive task when applied to these problems, since the occurrence of modal coupling demands a high number of degrees-of-freedom. A feasible possibility to overcome these difficulties is the use of low order models. The nonlinear normal modes have been shown to be an effective tool in the derivation of reduced order models in nonlinear dynamics. In the use of nonlinear modal analysis fewer modes are required to achieve a given level of accuracy in comparison to the use of linear modes. This work uses the nonlinear normal modes to derive low dimensional models to study the vibration of simplified models of offshore structures. Three examples are considered: an inverted pendulum, an articulated tower and a spar platform. Both free and forced vibrations are studied. The asymptotic and Galerkin-based methods are used to derive the normal modes. In addition, an alternative numerical procedure to construct such modes is proposed, which can be used to derive coupled modes. The solution stability is determined by the use of the Floquet theory, bifurcation and Mathieu diagrams, and Poincaré sections. The Poincaré sections are also used to investigate the multiplicity of modes and multimodes. The results obtained from the numerical integration of the original system are favourably compared with those of the reduced order models, showing the accuracy of the reduced models.

[pt] Neste trabalho estuda-se uma estratégia para o controle ativo não-linear de estruturas flexíveis submetidas a carregamentos dinâmicos. O algoritmo de controle é deduzido com base na teoria de controle ótimo não-linear com realimentação de estado, utilizando uma representação tensorial. Desenvolve-se as equações polinomiais de controle para diferentes ordens, partindo-se do controle linear clássico até um controle não-linear de quinta ordem. A estratégia é particularizada para aplicação em sistemas com um grau de liberdade que apresentem não-linearidades quadráticas e cúbicas, que podem representar, de forma aproximada, a maioria dos elementos estruturais encontrados nas Engenharias Civil e Mecânica, tais como vigas, arcos, placas e cascas. Determina-se analiticamente os coeficientes de ganhos até a terceira ordem, e utiliza-se os mesmos para estudar o feito das forças de controle sobre a não- linearidade e estabilidade do sistema. Vários exemplos numéricos de aplicação são apresentados, utilizando-se diferentes tipos de excitação. Uma atenção especial é dedicada a sistemas caracterizados pela coexistência de dois vales potenciais, um deles correspondente a uma posição de equilíbrio pós- flambagem. A influência do sistema de controle sobre a carga de escape é estudada. O efeito do retardo na aplicação das forças de controle é analisado tanto numericamente como analiticamente, utilizando-se o método das múltiplas escalas para desenvolver expressões que permitem encontrar situações críticas de retardo. Como exemplo de aplicação prática, estuda-se o problema de uma viga flambada submetida a um carregamento dinâmico lateral. / [en] The present thesis studies a strategy for the active non- linear control of dynamically loaded flexible structures. The control method is based on the non-linear optimal control theory using state feedback and the solution of the non-linear optimal control problem is obtained by representing system non-linearities and performance indices by power series with the help of algebraic tensor theory. General polynomial representations of the non-linear control law are obtained up to the fifth order. This methodology is applied to systems with quadratic and cubic nonlinearities, capable of representing most of the elements usually used in civil and mechanical engineering structures, such as beams, plates, shells and arcs. Control gains up to the third order are analytically derived and the effect of the control forces on the system is studied. Special emphasis is placed on systems susceptible to chaotic vibrations, escape from a potential well and dynamic jumps. Several examples are provided to illustrate the control approach. Strongly nonlinear systems subjected to free vibration, simple harmonic excitations, impact and ground acceleration are tested. The variation of the dynamic buckling load with the degree of the control algorithms is studied for the problem of structures with two potential wells, one of them corresponding to a post-buckling equilibrium position. The effect of time delay on controlled systems is studied analytically and numerically. The studied methodology is also applied to control the oscillations of simply supported buckled beams, in order to mitigate the effects of dynamic loading on the vibration amplitudes and prevent dangerous instability phenomena. / [es] En este trabajo estudia se una estrategia para el control activo no-lineal de extructuras flexibles sometidas la cargas dinámicas. EL algoritmo de control es deduzido con base en la teoría de control ótimo no lineal con realimentación de estado, utilizando una representación tensorial. Se dearrollan las ecuaciones polinomiales de control para diferentes órdenes, desde el control lineal clásico hasta el control no lineal de quinta orden. Se particulariza la estrategia para la aplicación en sistemas con un grado de liberdad que presenten no linealidades cuadráticas y cúbicas, que pueden representar, de forma aproximada, la mayoría de los elementos extructurales encontrados en las Ingenierías Civil y Mecánica, tales como vigas, arcos y placas. Se determinan analiticamente los coeficientes de ganancias hasta tercer orden, y se utilizan los mismos para estudiar las fuerzas de control sobre la no linealidad y estabilidad del sistema. Varios ejemplos numéricos de aplicación son presentados, utilizando diferentes tipos de excitación. Atención especial se le dedica a los sistemas caracterizados por la coexistencia de dos vales potenciales, uno de ellos correspondiendo a una posición de equilibrio posflameado. Se estudia también la influencia del sistema de control sobre la carga de escape. EL efecto de retardo en la aplicación de las fuerzas de control se anlaiza tanto numérica como analíticamente, utilizando el método de las múltiples escalas para dearrollar expresiones que permiten encontrar situaciones críticas de retardo. Como ejemplo de aplicación práctica, se estudia el problema de una viga flameada sometida a una carga dinámica lateral.

[pt] Com base nas equações de Donnell para cascas abatidas, estudam-se as vibrações não-lineares e a instabilidade dinâmica de cascas cilíndricas carregadas axialmente, considerando o efeito simultâneo de cargas estáticas e harmônicas. Para isto, o problema é primeiro reduzido a um sistema finito de graus de liberdade usando o método de Galerkin. O sistema de equações diferenciais de movimento não-lineares é resolvido através do método de Runge-Kutta. Na análise do fenômeno de acoplamento modal foi usado um modelo com dois graus de liberdade, que reflete de maneira consistente o comportamento pós-crítico inicial da casca. Os resultados obtidos com esse modelo simplificado foram comparados com diversas modelagens encontradas na literatura, permitindo assim aferir o modelo utilizado. Para o estudo da interação modal, escolheram-se três modelos distintos com freqüências e cargas críticas próximas ou semelhantes. Para estudar o comportamento não- linear da casca, diversas estratégias numéricas foram usadas para se obter os mapas de Poincaré, expoentes de Lyapunov, pontos fixos estáveis e instáveis, diagramas de bifurcação e bacias de atração. Particular atenção foi dada a dois fenômenos de instabilidade dinâmica que podem ocorrer sob estas condições de carregamento, a saber: excitação paramétrica dos modos de flexão e escape da bacia de atração pré-flambagem. Os cálculos foram realizados nas regiões principal e secundária de instabilidade paramétrica associadas com a menor freqüência natural da casca. Mostra-se de forma detalhada a determinação dos limites de instabilidade no espaço de controle e a identificação dos mecanismos de escape relacionados com estes limites. Os resultados mostram a importância do acoplamento e da interação modal na solução pós-crítica e no comportamento dinâmico não-linear de cascas cilíndricas. / [en] Based on Donnell shallow shell equations, the nonlinear vibrations and dynamic instability of axially loaded circular cylindrical shells under both static and harmonic forces are analyzed. For this, the problem is first reduced to that of a finite degree-of-freedom system by the Galerkin method. The resulting set of coupled non-linear ordinary differential equations of motion is, in turn, solved by the Runge-Kutta method. For the study of modal coupling, a 2 DOF model was used that describes consistently the initial post-buckling behavior of the shell. This model was compared favorably with others found in literature. For the analysis of modal interaction three different models were used considering shells with close or equal frequencies and critical loads. To study the non-linear behavior of the shell several numerical strategies were used to obtain Poincaré maps, Lyapunov exponents, stable and unstable fixed points, bifurcation diagrams and basins of attraction. Particular attention is paid to two dynamic instability phenomena that may arise under these loading conditions:parametric excitation of flexural modes and escape from the pre-buckling potential well.Calculations are carried out for the principal and secondary instability regions associated with the lowest natural frequency of the shell. Special attention is given to the determination of the instability boundaries in control space and the identification of the bifurcational events connected with these boundaries. The results also clarify the importance of modal coupling and modal interaction to the post-buckling solution and non-linear dynamic behavior of cylindrical shells.

[pt] Nesta tese são estudadas as vibrações não-lineares de membranas circulares inicialmente tracionadas sujeitas a deformações finitas. O material da membrana é modelado como um material hiperelástico neo-Hookeano, isotrópico e incompressível. Baseada na teoria de deformações finitas para membranas hiperelásticas, uma formulação variacional é desenvolvida. Primeiro a solução da membrana sob tração radial uniforme é obtida e então as equações de movimento da membrana são obtidas pelo princípio de Hamilton. A partir das equações linearizadas, as freqüências e os modos de vibração da membrana são obtidos analiticamente. Os modos naturais são usados para aproximar o campo de deformações não-linear usando o método de Galerkin e modelos de ordem reduzida são deduzidos através do método de Karhunen-Loève e de métodos analíticos. Além disso, estuda-se a influência da variação da massa específica e da espessura ao longo da direção radial da membrana nas vibrações. A seguir a mesma metodologia é utilizada para uma membrana anular. Por fim, estudam-se as vibrações não-lineares da membrana anular acoplada a uma inclusão rígida que insere tensões de tração na membrana, pois, devido ao seu peso próprio, provoca deslocamentos estáticos transversais e axissimétricos na membrana. Os mesmos problemas são analisados por elementos finitos utilizando o programa comercial Abaqus. / [en] This work presents an analysis of the nonlinear vibration response of a prestretched hyperelastic circular membrane subjected to finite deformations. The membrane material is assumed to be isotropic, homogeneous and neo-Hookean. Based on the theory of finite deformations for hyperelastic membranes, a variational formulation is developed. First the exact solution of the membrane under a uniform radial stretch is obtained and then the equations of motion of the pre-stretched membrane are derived using the Hamilton’s principle. From the linearized equations of motion, the natural frequencies and mode shapes of the membrane are obtained analytically. Then the natural modes are used to approximate the nonlinear deformation field using the Galerkin method. Several reduced order models are tested using the Karhunen-Loève method and analytical methods. Besides, the influence of the variation of the membrane thickness and material density along the radial direction of the membrane on the vibrations is investigated. The same methodology it is used for the annular membrane. Finally, the non-linear vibrations of the annular membrane coupled to a rigid inclusion are studied. The rigid inclusion inserts traction forces in the membrane and its own weight causes static transverse and radial displacements in the membrane. The same problems are analyzed by finite elements using the commercial program Abaqus®.

[pt] Arcos abatidos são usados com frequência para vencer grandes vãos. Exemplos incluem pontes em arco e coberturas de grandes espaços como galpões industriais e estádios. Em muitos casos empregam-se arcos atirantados ou apoiados em estruturas flexíveis, fazendo com que os apoios se movam quando o arco é carregado. Isto aumenta a flexibilidade do sistema e a probabilidade de perda de estabilidade na presença de cargas estáticas e dinâmicas. Em muitos casos estas estruturas podem ser modeladas como arcos com apoios elásticos. No presente trabalho resolve-se o problema de estabilidade estática de forma analítica e através de uma aproximação usando o método de Ritz, servindo a solução analítica para aferir a precisão do modelo numérico. A seguir, com base neste estudo, desenvolve-se, usando o método de Ritz, a formulação para análise das vibrações não lineares do arco com apoios elásticos, assunto inédito na literatura. Os resultados mostram a grande influência dos apoios nas vibrações não lineares e na estabilidade do arco sob cargas estáticas e dinâmicas. / [en] Shallow arches are often used to overcome large spans, for example, arch bridges or steel roofs to cover large spaces such as industrial sheds and stadiums. In many cases the arches are tied or are supported by a flexibility structure, causing that supports to move when the arch has been loaded. This increases the flexibility of the system and the probability of loss of stability in the presence of static and dynamic loads. In many cases, these structures can be modeled as arches with elastic supports. In the present work the static stability has been solved analytically and through the Ritz method, serving the analytical solution to assess the accuracy of the numerical model. Then, based on this study, the analysis of nonlinear vibrations of shallow arches with elastic supports is developed, using the Ritz method, a subject not yet studied in the literature. The results show the noticeable influence of the supports on the nonlinear vibration and stability of shallow arches under static and dynamic loads.

[pt] Um dos resultados fundamentais na mecânica clássica é que, para sistemas lineares com n graus de liberdade, existem n modos de vibração ortogonais e que as freqüências naturais são independentes da amplitude de vibração. Além disso, qualquer movimento da estrutura pode ser obtido como uma combinação linear desses modos. No caso de sistemas não-lineares, isto não mais se verifica e a relação entre freqüência, amplitude e os modos de vibração precisa ser determinada. A obtenção dessas informações para estruturas se dá em geral pelo uso de programas de análise não-linear baseados em uma formulação em elementos finitos. Contudo, isto é um procedimento custoso computacionalmente. Uma abordagem mais viável é o uso de modelos discretos compatíveis de baixa dimensão, por meio dos quais as freqüências e os modos não- lineares são obtidos. Neste trabalho é proposto um procedimento para a derivação de modelos de redução de dimensão para vigas e pórticos planos esbeltos. As equações diferenciais de movimento são obtidas a partir da aplicação das técnicas variacionais a um funcional não-linear de energia. A obtenção do modelo se dá através do emprego dos métodos de Ritz ou Galerkin para a redução espacial e do balanço harmônico para redução no tempo. Os modos lineares são utilizados como uma primeira aproximação para os modos não-lineares. As relações freqüência-amplitude são satisfatoriamente obtidas para vibrações livre e forçada (não-amortecida e amortecida). Entretanto, essas curvas apresentam, em geral, no regime não-linear, pontos limites, sendo obtidas, portanto, com uso do método do controle de comprimento de arco. Uma correção para o modo- linear é obtida com uso dos métodos dos elementos finitos e da perturbação. Um estudo paramétrico e das condições de contorno é apresentado para vigas. O comportamento não-linear de pórticos em L é também analisado. Para esses pórticos é estudada a influência de cargas axiais e da geometria. Os resultados são comparados com soluções analíticas encontradas na literatura. / [en] One of the fundamental results in classical mechanics is that linear systems with n degrees of freedom have n orthogonal vibration modes and n natural frequencies which are independent of the vibration amplitude. Any motion of the system can be obtained as a linear combination of these modes. This does not hold for nonlinear systems in which case amplitude dependent vibrations modes and frequencies must be obtained. One way of obtaining these informations for arbitrary structures is to use a nonlinear finite element software. However, this is a cumbersome and time consuming procedure. A better approach is to derive a consistent low dimensional model from which the nonlinear frequencies and mode shapes can be derived. In this work a procedure for the derivation of low dimensional models for slender beams and portal frames is proposed. The differential equations of motion are derived from the application of variational techniques to a nonlinear energy functional. The linear vibration modes are used as a first approximation for the nonlinear modes. The Galerkin and Ritz methods are used in the model for the spatial reduction and the harmonic balance method for the reduction in time domain. This allows the analysis of the free and forced (damped or undamped) vibrations of the structure in non- linear regime. However nonlinear resonance curves usually presents limit points. To obtain these curves, a methodology for the solution of non-linear equations based on an arc-length procedure is derived. Based on the finite element methods and using the basic ideas of the perturbation theory, a correction for the nonlinear vibration modes is derived. The influence of boundary conditions, geometric, and force parameters on the beam response is analyzed. The behavior of L frames is studied. For this kind of frame, the influence of axial loading and geometric parameters on the response is studied. The results are compared with analytical solutions found in the literature.

[pt] Treliças espaciais de forma piramidal são um componente básico de diversas estruturas, incluindo desde nanoestruturas de carbono até domos geodésicos para cobertura de grandes espaços. Estas estruturas, tal como a treliça plana de von Mises, apresentam uma resposta altamente não linear na presença de cargas estáticas e dinâmicas. A não linearidade é particularmente significativa, mesmo para baixos níveis de carregamento, quando estas estruturas são abatidas. Neste trabalho apresenta-se uma formulação não linear exata para uma treliça piramidal composta de n barras e, a partir desta formulação, analisa-se a perda de estabilidade e vibrações não lineares destas estruturas sob cargas estáticas e dinâmicas. Para compreensão do comportamento não linear são usadas respostas no tempo e planos de fase, diagramas de bifurcação, perfis de energia e bacias de atração. Os resultados mostram a importância da não linearidade na dinâmica e estabilidade da estrutura. / [en] Spatial trusses in the form of a regular pyramid are a basic configuration of many structures, from carbon nanostructures to geodesic domes used as roof of large spaces. These structures, as the planar von Mises truss, present a highly nonlinear response when submitted to static and dynamic loads. The nonlinearity is particularly significant when these structures are shallow, even when low load levels are applied. This paper presents a exact nonlinear formulation for a pyramidal truss made of n equal bars and analyze the loss of stability and nonlinear vibrations of the structure under static and dynamic loads. To understand to nonlinear behavior of the pyramidal truss, time responses, phase planes, bifurcation diagrams, energy profiles and basins of attraction are used. The results show the importance of the nonlinearity in the dynamics and stability of the structure.

[pt] Em virtude de sua eficiência, elementos estruturais de paredes finas com seções abertas são comuns em estruturas de aço, sendo secção em I, L, C e T usuais na prática de engenharia. A maior parte das vigas de parede fina tem uma boa resistência à flexão em relação ao eixo principal de inércia, mas uma baixa rigidez à flexão em relação ao eixo de menor inércia e uma baixa rigidez em torção. É por isso que estes elementos apresentam em geral uma instabilidade que leva a um acoplamento de flexo-torção. Muitas destas estruturas trabalham em um regime não linear e uma formulação não linear que leve em conta grandes deslocamentos e os acoplamentos inerentes é necessária. Neste trabalho um modelo não linear para vigas de seção aberta e paredes finas, considerando grandes deslocamentos, os efeitos de encurtamento e acoplamentos em flexão e torção é adotado. Inicialmente um estudo das frequências naturais, das cargas críticas e da relação frequência-carga axial é apresentado para diversos perfis. Com base nestes resultados, faz-se um estudo detalhado do comportamento dinâmico não linear destes perfis destacando o efeito do acoplamento não linear na região de ressonância e sua influência na estabilidade dinâmica da estrutura. Para isto são usadas diversas ferramentas de dinâmica não linear, tais como diagramas de bifurcação, respostas no tempo e plano de fase e bacias de atração. Os resultados mostram que a consideração dos acoplamentos não lineares é essencial para se avaliar o nível de segurança destas estruturas. / [en] Due to its efficiency, thin-walled structural elements with open sections, such as I, L, C and Z profiles, are common in steel structures, being usual in engineering applications. Most thin-walled beams have a good flexural strength around of the principal axis of inertia, but a low one around the axis of lower inertia as well as low torsional stiffness. That is why these elements, generally, show instabilities that lead to flexural torsional coupling. Many of these structures do not work in a linear range and a non-linear formulation that takes into account large displacements and associated couplings is required. This dissertation presents a nonlinear model for extensional beams with thin-walled open section, considering large displacements, and flexural-torsional couplings. Initially a study of the natural frequencies, critical load and axial load vs. frequency relation is presented for different profile kinds. Based on these results, a detailed study of the dynamic behavior of non-linear profiles is made, highlighting the effect of non-linear coupling in the resonance region and its influence on the dynamic stability of the structure. For this, various tools of nonlinear dynamics are used, such as bifurcation diagrams, time histories and phase-space portraits and basins of attraction. The results show that the consideration of non-linear couplings is essential to availed the safety level of these structures.

[pt] Cascas cilíndricas são usadas em muitas aplicações de engenharia e, devido a sua forma e capacidade de transporte de carga, são bastante usadas na indústria aeroespacial e em estruturas civis. Elas minimizam a quantidade de material do qual são fabricadas, tornando-se assim estruturas muito leves e esbeltas. Em décadas recentes tem se procurado criar novos materiais que conjuguem múltiplas propriedades como maior resistência, melhor proteção térmica, proteção contra corrosão e adequado nível de amortecimento, dentre outras. Uma classe de materiais que podem atender simultaneamente várias destas exigências é o chamado material com gradação funcional, onde as propriedades do material variam de forma contínua em uma ou mais direções. Materiais com gradação funcional são particularmente indicados para a construção de cascas. Como a maioria destas estruturas estão sujeitas a cargas dinâmicas, torna-se importante o estudo do comportamento dinâmico de cascas fabricadas com materiais com gradação funcional. O objetivo deste trabalho é estudar as vibrações não lineares de cascas cilíndricas esbeltas com gradação funcional. Para isto utiliza-se a teoria não linear de cascas de Sanders, considerada uma das teorias mais precisas para a análise de cascas esbeltas. Inicialmente, derivam-se as equações de movimento considerando um estado de tensões iniciais. Usando as equações linearizadas, obtêm-se às frequências naturais e as cargas críticas, sendo estes resultados comparados favoravelmente com resultados encontrados na literatura para materiais homogêneos e com gradação funcional. A seguir, usando uma expansão modal que atende as condições de contorno e continuidade, além de expressar os acoplamentos modais característicos de cascas cilíndricas no regime não linear, as equações de movimento são discretizadas usando-se o método de Galerkin. As equações algébricas resultantes são resolvidas pelo método de Newton-Raphson, sendo assim obtida a relação não linear frequência-amplitude. Finalmente, realiza-se uma análise paramétrica para estudar a influência da geometria da casca, da gradação do material funcional e dos modos de vibração no grau e tipo de não linearidade da casca cilíndrica, sendo esta a principal contribuição deste trabalho de pesquisa. / [en] Cylindrical shells are used in many engineering applications and, due to its shape and load carrying capacity, are frequently used in aerospace and civil structures. They minimize the amount of material from which they are manufactured, thus making it a very lightweight and slender structure. In recent decades, there has been a search for new materials that combine multiple properties such as increased strength, better thermal protection, corrosion protection and appropriate damping level, among others. A material that can meet several of these requirements simultaneously is the so called functionally graded material, where the material properties vary continuously in one or more directions. Functionally graded materials are particularly suitable for the construction of shells. As most of these structures are subjected to dynamic loads, it is important to study the dynamic behavior of shells made of functionally graded materials. The objective of this work is to study the nonlinear vibrations of slender functionally graded cylindrical shells. For this, the Sanders non-linear shell theory, which is considered one of the most precise theories for the analysis of slender shells, is adopted. Initially, the equations of motion are derived considering an initial stress state. Using the linearized equations of motion, the natural frequencies and critical loads are obtained. These results compare favorably with results reported in the literature for homogeneous and functionally graded shells. Then, using a modal expansion that satisfies the boundary and continuity conditions and expresses the modal couplings characteristic of cylindrical shells in the nonlinear regime, the equations of motion are discretized using the Galerkin method. The resulting algebraic equations are solved by the Newton-Raphson method, thus obtaining the nonlinear frequency-amplitude relation. Finally, a parametric analysis is conducted to study the influence of the geometry of the shell, the gradient of the functional material and vibration modes on the degree and type of nonlinearity of the cylindrical shell, which is the main contribution of this research work.