Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / Os perfis tubulares sem costura são largamente utilizados em diversos países, porém, no Brasil, o uso desses perfis era bastante limitado, restringindo-se praticamente a coberturas espaciais. A situação do mercado brasileiro começou a se alterar em razão da maior oferta desses perfis proporcionada pelo início da produção pela V&M do BRASIL no ano de 2000, sendo assim diante da novidade da tecnologia, impõe-se a necessidade de divulgação e implementação do uso desse tipo de perfil, além de uma ampliação do número de trabalhos de pesquisa para melhor compreensão de seu comportamento estrutural. O presente trabalho apresenta uma análise paramétrica de ligações tipo KT e T, através da determinação da resistência última da ligação obtida pelo critério de deformação limite proposto por diversos pesquisadores, comparando os resultados com as prescrições do Eurocode 3, do Projeto de Norma Brasileira e do CIDECT. Adicionalmente procedeu-se na comparação do comportamento numérico global do nó de uma treliça, com a análise do comportamento numérico do nó isolado. As ligações objeto deste estudo são compostas por perfil tubular retangular (RHS) no banzo e perfis circulares (CHS) nas diagonais e montantes. Os modelos em elementos finitos, foram desenvolvidos no programa Ansys 12.0, as não-linearidades físicas e geométricas foram incorporadas aos modelos, a fim de se mobilizar totalmente a capacidade resistente da ligação. Concluiu-se entre outras coisas que. para o nó da treliça em estudo neste trabalho, o carregamento de projeto máximo permitido através do dimensionamento pelas normas fica a favor da segurança quando comparado com o critério de deformação limite. / The hollow section seamless are widely used in many countries, however, in Brazil, the use of these section was rather limited, being restricted to just the spatial roof. The situation of the Brazilian market began to change due to the increased supply of these section provided by the start of production by V&M do BRASIL in 2000, thus before the novelty of the technology, imposes the need to disseminate and implement the use of section type, and an increased number of research studies to better understand its structural behavior. This study presents a parametric analysis of joints type "KT" and "T", by determining the ultimate strength of the joint obtained by the criterion of deformation limit proposed by several researchers, comparing the results with the requirements of Eurocode 3, Project Brazilian Standard and CIDECT Additionally we proceeded to compare the behavior of the global node number of a truss, with the numerical analysis of the behavior of the isolated node. The joint object of this study are composed of rectangular hollow section (RHS) in the chord and circular hollow section (CHS) on the brace member. The finite element models were developed in the program Ansys 12.0. The non-linearities physical and geometric have been incorporated in the models, to mobilize the full load capacity of the joint. It was concluded among other things for the node of the truss studied in this work, the maximum allowable loading project through the design standards is to promote security when compared with the criterion limit of deformation.

Diversos pesquisadores têm estudado o comportamento e o emprego de aduelas de concreto, que constituem as vigas segmentadas em sistemas estruturais, de maneira especial pontes e viadutos. Por esta razão, inúmeros trabalhos têm sido publicados nos últimos anos respaldados por testes experimentais e análises numéricas. O comportamento destas vigas contrasta com as clássicas vigas monolíticas em diversos aspectos, pois, a estrutura é composta de partes de elementos de concreto pré-moldado que, após serem posicionados no local definitivo, são protendidos. A protensão pode ser aderente ou não aderente. A principal vantagem deste sistema de construção é a rapidez e o alto controle de qualidade, por isso é largamente utilizado, havendo uma demanda de estudo de previsão do seu real comportamento No presente trabalho apresenta-se uma modelagem numérica via elementos finitos, para simular o comportamento de vigas compostas por aduelas justapostas sem material ligante entre as juntas. A protensão aplicada é aderente e a análise considera a não linearidade da região da junta. Assim sendo, o objetivo desta investigação é dar uma contribuição ao estudo do comportamento estrutural estático de vigas segmentadas, atentando para o comportamento das juntas, utilizando um programa comercial. Para o modelo são empregadas técnicas usuais de discretização, via método dos elementos finitos (MEF), por meio do programa de elementos finitos SAP2000[93]. O modelo proposto é constituído de elementos de placa próprios para concreto para representar a viga, a protensão é introduzida por meio de barras bidimensionais que transferem as tensões ao longo de seu comprimento e as juntas são implementadas utilizando elementos de contato. A analise é bidimensional e considera os efeitos das perdas de protensão. Este trabalho de pesquisa objetiva também o estudo de elementos de contato especialmente as características de deformação para esta ferramenta computacional. A definição dos parâmetros para o modelo foi feita com base em dados experimentais disponíveis na literatura. O modelo numérico foi calibrado e confrontado com resultados experimentais obtidos em ensaios de laboratório. / Several researchers have studied the behavior and the use of concrete staves, which are segmented beams in structural systems, especially bridges and overpasses, in civil engineering. For this reason, numerous studies have been published in recent years supported by experimental tests and numerical analyzes. The behavior of these beams in contrast to the conventional monolithic beams in different ways because the structure is composed of portions of pre-cast concrete elements which after being placed in situ are prestressed. The prestressing can be adhered or otherwise attached. The main advantage of this system of construction is speed and high quality control, so it is widely used; there is a demand forecast study their actual behavior. In this paper, we present a finite element model to simulate the behavior of composite beams prestressed by adherent with staves. Therefore, the objective of this research is to make a contribution to the study of static structural behavior of beams targeted observing the behavior of joints using a commercial program. For the model are employed standard techniques of discretization via the finite element method (FEM), using the finite element program SAP2000 [90]. The model is made of elements suitable for plaque to represent the actual beam, the prestressing is introduced by means of two-dimensional bar which transfer the stresses along its length and the joints are implemented using contact elements proposed by the program, called Link. The analysis is two dimensional and considers the effects of prestressing losses. This research also aims to study the contact elements especially the characteristics of deformation to this computational tool. The definition of parameters for the model was based on experimental data available in literature. The numerical model was calibrated and compared with experimental results obtained in laboratory trials.

Este trabalho apresenta uma formulação tensorial genérica para parametrização das rotações do tipo vetorial destinada ao estudo de grandes rotações no espaço tridimensional. Esta formulação é compatível com as parametrizações de Euler e de Rodrigues. É dada ênfase aos que aqui se denominou parâmetros generalizados de Rodrigues, que fornecem expressões simples, computacionalmente mais eficientes que a parametrização clássica de Euler. A formulação apresentada é adequada para uso em métodos numéricos baseados nas projeções de Galerkin, como o método dos elementos finitos, podendo ser implementada com facilidade em programas já existentes de elementos finitos. Apresentam-se aqui expressões para o tensor das rotações e suas derivadas, bem como os tensores necessários à análise incremental. As formas fracas são construídas tanto com projeção ortogonal como não-ortogonal, correspondentes à aplicação do Teorema dos Trabalhos Virtuais e Teorema das Potências Virtuais, respectivamente. Os modelos propostos foram aplicados em um programa de elementos finitos utilizando formulações cinemáticas Lagrangiana total e Lagrangiana atualizada e foram resolvidos vários exemplos, dentre eles alguns clássicos da literatura, de forma a avaliar sua validade e aplicação. / This work presents a generic formulation of vector-type for the parameterization of large rotations in three-dimensional space. This formulation adapts to the Euler and the Rodrigues parameterization. Special distinction is made to the here named generalized Rodrigues parameters which result in very simple and computationally efficient expressions. The attained formulation is convenient for numerical procedures employing Galerkin projection like the finite element method and can be readily implemented in a FE code. The expressions of rotation tensor and its derivatives, which lead to a consistent linearization, are herein derived. The necessary tensor quantities employed in the derivation of the tangent bilinear weak form of incremental analysis are obtained too. The weak forms are constructed here with both orthogonal and non-orthogonal projections, corresponding to the application of the virtual work theorem or virtual power theorem respectively. The formulation is implemented within a finite element code in total Lagrangian and updated Lagrangian framework and assessment of the scheme is made by means of several numerical simulations.

Neste trabalho foi desenvolvida uma plataforma computacional para análise via método dos elementos finitos de estruturas de concreto armado convencional e reforçado com fibras de aço. A ferramenta numérica desenvolvida foi obtida por meio do acoplamento do programa FEMOOP, denominado solver com o pré e pós-processador GiD. Esse acoplamento foi possibilitado por meio da programação de um conjunto de arquivos denominados arquivos de customização, responsáveis pelo trabalho conjunto dos programas. Utiliza-se uma única interface gráfica com caixas de diálogo vinculadas ao código do solver, responsáveis por aplicar as condições de contorno do problema, tipo de análise, e aplicação dos materiais nos seus respectivos elementos finitos. Para a representação do concreto, foram implementados elementos finitos planos isoparamétricos quadrilaterais e triangulares e para as armaduras elementos finitos isoparamétricos unifilares lineares e quadráticos representados por meio do modelo discreto. Para o comportamento do concreto, foi considerado um modelo elástico não-linear com comportamento isotrópico até o limite de ruptura, acoplado a um modelo de amolecimento linear na tração. As fissuras são representadas pelo modelo de fissuração distribuída do tipo rotacional. Como critério de resistência para o concreto podem-se usar o modelo de Ottosen ou o modelo de Willam e Warnke de cinco parâmetros implementados na plataforma. Especificamente para considerar a presença de fibras de aço descontínuas na matriz de concreto, é utilizado o critério de ruptura proposto por SEOW e SWADDIWUDHIPONG (2005), que é uma alteração no meridiano de tração do critério proposto por Willam e Warnke. Para o concreto reforçado com fibras de aço fissurado considera-se o trecho pós-fissuração proposto no modelo de tração de LIM et. al (1987). As armaduras têm seu comportamento descrito através de um modelo elasto-plástico bilinear. A interação entre as armaduras e o concreto foi considerada como de aderência perfeita. Como se trata da modelagem de um material com comportamento não-linear, foi implementado para resolução das equações de equilíbrio o método de Newton-Raphson. Por fim, a plataforma final obtida foi avaliada por meio da simulação de vigas de concreto armado convencional e reforçado com fibras de aço disponíveis na literatura, que confirmaram a eficiência das implementações efetuadas. / In this work a computational platform for the analysis of reinforced concrete structures reinforced or not with steel fibers has been developed. This tool is based on the finite element method and has been obtained by the coupling of FEMOOP, denominated solver, with the pre and post-processor program GiD. The coupling has been possible by programming a set of customization files responsible for the communication between the two base programs. A single graphical interface with particular dialog boxes which are linked to the solver facilities is used to apply the boundary conditions, type of analysis, and material properties in the finite element model. For the geometrical representation of concrete elements, plane isoparametric quadrilateral and triangular finite elements have been implemented, while for the steel reinforcement bars, discrete isoparametric truss finite elements with linear end quadratic interpolation have been used. In order to model the mechanical behavior of concrete materials, a nonlinear isotropic elastic model together with a tension softening linear model has been adapted. Cracks are represented through a rotational smeared crack model. Both Ottosen and 5 parameters Willam-Warnke models can be used as the strength criterion of concrete. A failure model proposed by SEOW and SWADDIWUDHIPONG (2005), based on an adaptation of the Willam-Warnke model where a modification of the tension meridian is introduced, is used to consider the discontinuous steel fibers dispersed into the concrete mass. The post-cracking behavior of the steel fiber reinforced concrete considers the tension model proposed by LIM et. al (1987). The steel rebars have their behavior described by a bilinear elastoplastic model. A perfect bond between concrete and the reinforcing bars is assumed. For the solution of the nonlinear equations the Newton-Raphson method is used. The developed computational platform has been evaluated through a set of numerical simulations of tests performed in conventionally reinforced and steel fiber reinforced concrete beams available on the literature. These simulations confirm the efficiency of the current implementation.

Mais recentemente a exploração offshore de petróleo e gás tem se intensificado na costa brasileira. Uma das características da exploração offshore no Brasil está ligada ao fato de os hidrocarbonetos se situarem a grandes profundidades no mar, frequentemente abaixo de uma espessa camada de sal, exigindo o desenvolvimento de novas tecnologias para vencer os desafios para sua prospecção e extração. A profundidade dos poços pode variar entre dois mil e sete mil metros abaixo da superfície do mar. Um dos grandes desafios é garantir uma boa conexão dos equipamentos de extração com a unidade de produção na superfície, o que se faz com risers. Há distintas configurações de risers, entre as quais os verticais, em catenária, lazy waves ou steep waves, sendo que o riser em catenária e o vertical, em particular, serão estudados aqui. Este trabalho visa a estudar os efeitos dinâmicos não lineares no riser, que, por ser extremamente esbelto, exige a consideração da não linearidade geométrica. Entre eles, destaca-se a possibilidade de ocorrência da ressonância paramétrica. Os carregamentos dinâmicos podem ser provenientes da movimentação da unidade flutuante, correntezas ou escoamento interno. Neste trabalho serão elaborados modelos em elementos finitos, utilizando o software Abaqus, que possui mais recursos que a maioria dos softwares comerciais para as análises que se pretende realizar. Esses modelos simularão inicialmente a configuração de equilíbrio estático do riser e posteriormente serão adicionados os carregamentos dinâmicos. Na representação do encontro do riser com o solo no fundo do mar, serão utilizados elementos de contato elásticos. Pretende-se discutir a viabilidade e as limitações no uso de programas generalistas de análise estrutural pelo método dos elementos finitos, no confronto com programas dedicados, e ainda comparando os resultados obtidos com aqueles que decorrem de soluções analíticas, ou experimentais, quando disponíveis. / More recently the offshore oil and gas exploitation has being intensified on the Brazilian coast. One of the characteristics of offshore exploitation in Brazil is linked to the fact that hydrocarbons are located at great depths under the sea bed, often below a thick layer of salt, requiring the development of new technologies to meet these challenges. The depth of the wells can vary between 2,000 and 7,000 meters below the sea surface. A major challenge is to ensure a good connection between the extraction system and the production unit at the sea surface, which is provided by risers. There are different configurations of risers, including the vertical, catenary, lazy-wave or steep-wave; the catenary and vertical risers will be studied in this work. This work aims at studying the nonlinear dynamic effects on the riser, which, being extremely slender, requires consideration of the geometric nonlinearity. Among them, there is the possibility of occurrence of parametric resonance. The dynamic loads may come from the motion of the floating unit and external or internal fluid flow. This work will be developed using finite-element models in Abaqus, which has more resources than most commercial softwares to carry out the type of analyses of interest. Initially these models simulate the static equilibrium configuration of the catenary riser, which will be later disturbed by dynamic loads. The modeling of the so called touch-down zone (TDZ) will use elastic contact elements. It is intended to discuss the feasibility and limitations when using generalist codes for structural analysis by the finite-element method, in comparison with dedicated codes, and also to correlate these results with those from analytical and experimental studies, whenever available.

Nesta tese discute-se a necessidade de se desenvolver novas ferramentas para auxiliar o projeto e análise de estruturas de cabos e membranas. Esse tipo de estrutura, essencialmente não linear, é geralmente analisada por meio do Método dos Elementos Finitos, combinado com o Método de Newton-Raphson, para a resolução do sistema de equações não lineares resultante. Porém, a ausência de um campo de tensão de tração sobre toda estrutura composta por elementos finitos de cabos e membranas pode gerar uma matriz de rigidez tangente indeterminada, levando à divergência da solução pelo Método de Newton-Raphson. O Método da Relaxação Dinâmica é uma alternativa interessante para resolver problemas não lineares complicados de equilíbrio estático, na qual o problema do equilíbrio estático é resolvido por uma análise dinâmica, com integração no tempo. A resposta transiente é fictícia e não tem significado físico, entretanto a parte estacionária é a solução do problema de equilíbrio estático. Nesta tese, apresenta-se uma contextualização histórica sobre o Método da Relaxação Dinâmica, apontando as contribuições mais relevantes já desenvolvidas por outros autores. Propõe-se um procedimento de sintonia da massa dos elementos, capaz de uniformizar as condições impostas ao incremento de tempo, para se obter estabilidade do processo de integração numérica. Implementam-se as formulações dos elementos finitos adotados, bem como um algoritmo de enrugamento para os elementos de membrana e diversas rotinas de pós-processamento, no programa de elementos finitos SATS (A System for the Analysis of Taut Structures), desenvolvido pelo autor desta tese, em colaboração com seu orientador. A implementação desenvolvida é aplicada a uma série de exemplos relativos ao projeto e análise de estruturas de cabos e membranas, permitindo verificar a eficiência dos procedimentos de amortecimento e cinético e de sintonia de massa propostos. / This thesis discusses the need to develop new tools to assist the design and analysis of cables and membrane structures. This type of structures, essentially non-linear is generally analyzed using the Finite Element Method, where in most cases the solution is obtained by the Newton-Raphson Method. However, the absence of a tension stress field over the entire structure composed only with cable and membrane finite element can generate a non-positive definite tangent stiffness matrix, leading to the divergence of Newton-Raphson iterations. The Method of Dynamic Relaxation is an interesting alternative to solve complicated nonlinear problems of static equilibrium, replaced by an equivalent dynamic analysis. The transient solution is fictitious and without physical meaning, and the stationary phase provides the static equilibrium solution. This thesis presents a historical contextualization of the Dynamic Relaxation Method, highlighting the most relevant contributions already developed by other authors. A procedure for the tuning of the element masses is proposed, which is capable of making uniform the restrictions imposed to the time steps in order to preserve the stability of the numerical integration. Some adopted finite element formulations are implemented, as well as an algorithm for representing the wrinkling of membrane elements and several post-processing routines, in the SATS (A System for the Analysis of Taut Structures) finite element program, developed by the author of this thesis, in collaboration with his advisor. The developed implementation is applied to a series of examples on the design and analysis of cables and membrane structures, allowing verification of the efficiency of the procedures proposed for kinetic damping and mass tuning.

O objetivo deste trabalho é aplicar uma formulação viscoelástica bidimensional sob não linearidade geométrica para estado plano de deformação baseada em grandezas energicamente conjugadas a partir do gradiente da transformação e do primeiro tensor das tensões de Piola- Kirchhoff. Esta formulação viscoelástica é composta pela equação constitutiva hiperelástica de Kircchoff- Saint Venant e a função de fluência do modelo viscoelástico de Kelvin-Voigt. A formulação é aplicada em análises numéricas com método dos elementos finitos. Para tal aplicação, foi proposto o uso de elementos triangulares T6 que não apresentaram o fenômeno de travamento e são adequados para a utilização em malhas não estruturadas. Esta implementação foi aplicada para prever deslocamentos no concreto projetado de túneis. / The objective of this work is to apply a geometrically nonlinear two-dimensional viscoelastic formulation for plane strain where deformation gradient is energetically conjugated with first Piola-Kircchoff stress tensor. This formulation is from constitutive equation St. Venant Kirchhoff and creep compliance Kelvin-Voigt rheological model. This is applied in numerical analysis with finite element method. Use of triangular elements was proposed. Locking is not observed in numerical examples and this kind of element is appropriate for use in unstructured mesh. Displacement prediction of shotcrete support in a tunnel was made with this formulation and excellent results were found.

O objetivo precípuo deste estudo é a implementação computacional de modelos constitutivos elásticos e elastoplásticos para materiais com gradação funcional em regime de grandes deslocamentos e elevadas deformações. Para simular numericamente um problema estrutural, são empregados aqui elementos finitos sólidos (tetraédrico e hexaédrico) com ordem de aproximação polinomial qualquer. Grandezas da Mecânica Não Linear do Contínuo, como deformação e tensão, são utilizadas na formulação deste estudo. Para reproduzir os grandes deslocamentos, é empregada a análise não linear geométrica. A descrição adotada aqui é a Lagrangiana total, e o equilíbrio da estrutura é expresso pelo Princípio da Mínima Energia Potencial Total. Com relação à resposta elástica do material, são usadas leis constitutivas hiperelásticas, nas quais a relação tensão-deformação é obtida a partir de um potencial escalar. O comportamento elastoplástico do material é definido pela decomposição da deformação nas parcelas elástica e plástica, pelo critério de plastificação de von-Mises, pela lei de fluxo associativa, pelas condições de consistência e de complementaridade, pelo parâmetro de encruamento isotrópico e pelo tensor das tensões inversas, relacionado ao encruamento cinemático. Duas formulações elastoplásticas são utilizadas aqui: a de Green-Naghdi, na qual a deformação é decomposta de forma aditiva; e a hiperelastoplástica, em que o gradiente é decomposto de forma multiplicativa. É empregado também o conceito de material com gradação funcional (GF), a qual é definida como a variação gradual (contínua e suave) das propriedades constitutivas do material. A solução numérica do equilíbrio de forças é feita via método iterativo de Newton-Raphson. Para satisfazer o critério de plastificação, são utilizadas as estratégias de previsão elástica, e de correção plástica via algoritmos de retorno. Basicamente foram desenvolvidos cinco programas computacionais: o gerador automático das funções de forma; o gerador de malhas de elementos finitos sólidos; o código para análise de materiais em regime elástico; o código para análise de materiais em regime elastoplástico; e o programa de pós-processamento. Além desses, o aluno teve contato com os programas EPIM3D e DD3IMP ao longo de seu estágio de doutorado na Universidade de Coimbra (Portugal). Os programas EPIM3D e DD3IMP são empregados para analisar, respectivamente, materiais em regime elastoplástico, e processos de conformação de metais. Para o problema da barra sob tração uniaxial uniforme, são descritas equações e soluções analíticas para materiais homogêneos e com GF em regime elastoplástico. Para reduzir o tempo de simulação, foi empregada a programação em paralelo. De acordo com os resultados das simulações numéricas, as principais conclusões são: o refinamento da malha de elementos finitos melhora a precisão dos resultados para materiais em regimes elástico e elastoplástico; as formulações elastoplásticas de Green-Naghdi e hiperelastoplástica parecem ser equivalentes para pequenas deformações; a formulação hiperelastoplástica é equivalente ao modelo mecânico dos programas EPIM3D e DD3IMP para materiais em regime de pequenas deformações elásticas; foram constatados ganhos significativos, em termos de tempo de simulação, com a paralelização dos códigos computacionais de análise estrutural; e os programas desenvolvidos são capazes de simular - com precisão - problemas complexos, como a membrana de Cook e o cilindro fino transversalmente tracionado. / The main objective of this study is the computational implementation of elastic and elastoplastic constitutive models for functionally graded materials in large deformation regime. In order to numerically simulate a structural problem, the finite elements used are solids (tetrahedric and hexahedric) of any order of approximation. Entities from Nonlinear Continnum Mechanics, as strain and stress, are used in the present formulation. To reproduce the finite displacements, the geometrically nonlinear analysis is employed. The description adopted here is the total Lagrangian, and the structural equilibrium is expressed by means of the Principal of Minimum Total Potential Energy. Regarding the elastic material response, hyperelastic constitutive laws are used, in which the stress-strain relation is obtained from a scalar potential. The elastoplastic material behavior is defined by the strain decomposition in the elastic and plastic parts, by the von-Mises yield criterion, by the associative flow law, by the consistency and complementarity conditions, by the isotropic hardening parameter, and by the backstress tensor, related to the kinematic hardening. Two elastoplastic formulations are used here: the Green-Naghdi one, in which the strain is additively decomposed; and the hyperelastoplasticiy, in which the gradient is multiplicatively decomposed. The concept of functionally graded (FG) material, in which the constitutive properties vary gradually (continuous and smoothly), is also used. The numerical solution of the forces equilibrium is obtained via Newton-Raphson iterative procedure. In order to satisfy the yield criterion, the strategies of elastic prediction and plastic correction (via return algorithms) are used. Basically, five computer codes have been developed: the automatic shape functions generator; the solid mesh generator; the code for analysis of materials in the elastic regime; the code for analysis of materials in the elastoplastic regime; and the post-processor. Besides these, the student had contact with the programs EPIM3D and DD3IMP during his doctoral stage in the University of Coimbra (Portugal). The programs EPIM3D and DD3IMP are employed to analyze, respectively, materials in the elastoplastic regime, and sheet-metal forming processes. For the problem of the bar under uniform uniaxial tension, equations and analytical solutions are described for homogeneous and FG materials. To reduce the simulation time, the parallel programming has been employed. According to the numerical simulation results, the main conclusions are: the results accuracy is improved with mesh refinement for materials in the elastic and elastoplastic regimes; the Green-Naghdi elastoplastic formulation and the hyperelastoplasticity appear to be equivalent for small strains; the hyperelastoplastic formulation is equivalent to the mechanical model of the programs EPIM3D and DD3IMP for materials the small elastic strains regime; simulation time reduction has been obtained with the parallelization of the computer codes for structural analysis; the developed programs are capable of simulating, precisely, complex problems, such as the Cook\'s membrane and the pulled thin cylinder.

O trabalho trata da formulação e implementação numérica de modelos matemáticos do comportamento de estruturas considerando-se as não-linearidades física e geométrica. O equilíbrio na posição deslocada é formulado via Princípio dos Trabalhos Virtuais, empregando-se o método dos elementos finitos para a discretização espacial das estruturas e busca de soluções aproximadas. Inicialmente destaca-se com base no caso de treliças planas o emprego de medidas de deformação e tensão conjugadas energeticamente. Particularizando-se a formulação geral do equilíbrio para os pórticos planos apresenta-se uma análise crítica das formulações lagrangiana total e atualizada. Em seguida, tendo-se em vista aplicações às estruturas em concreto armado, aborda-se o comportamento não-linear físico pela mecânica do dano em meios contínuos, empregando-se os modelos de dano para o concreto propostos por Mazars e La Borderie. Estendem-se os estudos do comportamento estrutural não-linear físico (dano) e geométrico incorporando-se a análise dinâmica. Utiliza-se para integração no domínio do tempo o método implícito de Newmark combinado com o procedimento incremental e iterativo de Newton-Raphson. O amortecimento é levado em conta por meio da regra de Rayleigh. Exemplos consistindo de análises não lineares estática e dinâmica de estruturas reticulares planas compõem numéricas. Os resultados obtidos ilustram o desempenho e as potencialidades das formulações empregadas. / The formulation and numerical implementation of mathematical models of the structural behavior of plane frames considering material and geometrical nonlinearities are treated in this work. The Principle of Virtual Work is presented in order to characterize the equilibrium in the displaced position. The structure is discretized by the finite element method. Plane trusses analysis is performed in order to show how important it is to take into account the strain and stress tensors energetically conjugated. Being particularized, the general equilibrium formulation for the plane frames a critical analysis of the total and updated lagrangian formulations is presented. The material nonlinear behavior is modeled by continuum damage mechanics by using both Mazars and La Boderie\'s damage models. The nonlinear analysis is extended including the dynamic response of the reinforced concrete plane frames. The usual iterative Newton-Raphson technique is used combined with implicit Newmark method in order to carry out the integration at time. The damping is introduced by means of the Rayleigh\'s rule. Numerical analysis by a suitable computer program show the theoretical results considering static and dynamic response of plane frames.

Sabe-se que materiais compostos laminados são, hoje em dia, geralmente usados nas indústrias aeronáutica, aeroespacial, naval e outras, principalmente por causa de suas atrativas propriedades se comparadas aos materiais isotrópicos, como alta rigidez/peso, alta resistência, alto amortecimento e boas propriedades relacionadas ao isolamento térmico e acústico, entre outras. Porém, o comportamento de estruturas feitas de materiais compostos pode ser aperfeiçoado através da utilização de materiais inteligentes. Dentre os diferentes tipos comercialmente disponíveis de materiais inteligentes, os materiais piezelétricos são amplamente usados como sensores e atuadores para o monitoramento e controle de estruturas. O efeito piezelétrico direto define que uma deformação mecânica aplicada ao material é convertida em uma carga elétrica. Por outro lado, o efeito piezelétrico inverso define que um potencial elétrico aplicado ao material é convertido em deformação mecânica. Estes efeitos governam a interação eletromecânica nos materiais piezelétricos. O Método dos Elementos Finitos, uma ferramenta amplamente reconhecida e poderosa para a análise de estruturas complexas, é capaz de realizar a integração dos componentes inteligentes e das partes estruturais clássicas. Sendo assim, o comportamento estático e dinâmico, linear e geometricamente não-linear, de estruturas compostas laminadas delgadas com lâminas piezelétricas incorporadas é analisado neste trabalho usando o Método dos Elementos Finitos (MEF). Elementos triangulares, chamados GPL-T9, com três nós e seis graus de liberdade por nó (três componentes de deslocamento e três de rotação) e um grau de liberdade por camada piezelétrica (potencial elétrico) são usados. Para a análise estática não-linear as equações de equilíbrio são solucionadas usando o Método do Controle de Deslocamentos Generalizados (MCDG) enquanto a solução dinâmica é obtida usando o Método de Newmark com Formulação Lagrangeana Atualizada (FLA). O sistema de equações é resolvido usando o Método dos Gradientes Conjugados (MGC) e nos casos não-lineares um esquema iterativo-incremental é empregado. Diversos exemplos numéricos são apresentados e comparados com resultados obtidos por outros autores com diferentes tipos de elementos e diferentes formulações. A concordância entre estes resultados demonstra a validade e a eficácia dos modelos desenvolvidos. / It is well known that laminate composite materials are nowadays commonly used in the aeronautical, aerospace, naval and other industries mainly because their attractive properties as compared to isotropic materials, such as higher stiffness/weight, higher strength, higher damping and good properties related to thermal or acoustic isolation, among others. However, the behavior of structures made of composite materials can be improved using smart materials. Among several kinds of commercially available smart materials, the piezoelectric materials are widely used as sensors and actuators for the monitoring and control of structures. The direct piezoelectric effect states that a mechanical strain applied to the material is converted to an electric charge. On the other hand, the converse piezoelectric effect states that an electric potential applied to the material is converted to mechanical strain. These effects govern the electromechanical interaction in piezoelectric materials. The finite element method, a widely accepted and powerful tool for analyzing complex structures, is capable of dealing with the integration of smart components and classic structural parts. So, linear and geometrically nonlinear static and dynamic behavior of thin laminate composite structures embedded with piezoelectric layers are analyzed in this work using the Finite Element Method (FEM). Triangular elements, called GPL-T9, with three nodes and six degrees of freedom per node (three displacement and three rotation components) and one degree of freedom per piezoelectric layer (electrical potential) are used. For static analysis the nonlinear equilibrium equations are solved using the Generalized Displacement Control Method (GDCM) while the dynamic solution is performed using the classical Newmark Method with an Updated Lagrangean Formulation (ULF). The system of equations is solved using the Gradient Cojugate Method (GCM) and in nonlinear cases an iterative-incremental scheme is employed. Several numerical examples are presented and compared with results obtained by other authors with different kind of elements and different schemes. The agreement among these results demonstrates the validity and effectiveness of the developed models.