Orientador: Osvaldo Luís Manzoli / Resumo: Este trabalho apresenta uma proposta de análise numérica do concreto com agregados reciclados (CAR). A ferramenta desenvolvida facilita a avaliação do comportamento mecânico do material, verificando a viabilidade de sua utilização para fins estruturais. Considerando que no Brasil o emprego atual do material ainda é bastante limitado em comparação com outros países, o CAR pode vir a ser utilizado mais ampla e adequadamente, trazendo vantagens de caráter estrutural, econômico, e principalmente ambiental. Em termos de modelagem do CAR, a ocorrência de particularidades em nível mesoscópico exige um modelo em escala mais refinada para representar as propriedades alteradas do agregado reciclado (AR), de modo que foi desenvolvido um gerador de agregados reciclados para representar sua composição variada. Os mecanismos de falha que ocorrem nas interfaces características do agregado reciclado e da matriz de concreto introduzem não linearidade ao problema mecânico. Para que haja a representação dessa não linearidade, a estratégia de modelagem proposta recorre ao emprego da técnica de fragmentação de malha de elementos finitos juntamente com um modelo constitutivo de dano. Essa técnica de fragmentação baseia-se no emprego de elementos finitos com alta razão de aspecto, os quais poderão ser utilizados para representar as fases adicionais desse material compósito, correspondentes às diferentes zonas de transição interfaciais (ZTIs) intrínsecas do CAR em mesoescala, e consequentemente, cam... (Resumo completo, clicar acesso eletrônico abaixo) / Abstract: This research presents a numerical analysis proposal for recycled aggregate concrete (RAC), whose recycled aggregates are obtained from the crushing of concrete waste. The developed tool facilitates the mechanical behavior evaluation of the RAC, verifying the feasibility of its use for structural purposes. Considering that in Brazil the current use of the material is still quite limited in comparison with other countries, the RAC could be used more widely and properly, bringing structural, economic and mainly environmental advantages. Due to the occurrence of mesoscopic specificities, the RAC modeling requires a more refined scale model to represent the altered properties of recycled aggregate, so that a recycled aggregate generator was developed to represent its variable composition. The failure mechanisms that occur at the several interfaces of the recycled aggregate (RA) and concrete matrix introduce nonlinearity to the mechanical problem. In order to represent this nonlinearity, the proposed modeling strategy uses the finite element mesh fragmentation technique together with robust and stable damage constitutive model. This fragmentation technique is based on the use of finite elements with high aspect ratio, which can be used to represent the additional phases of this composite material, corresponding to the different mesoscale RAC intrinsic interfacial transition zones (ITZs), as well as potential pathways for the propagation of fractures. A constitutive model based on ... (Complete abstract click electronic access below) / Mestre

[pt] O método híbrido de elementos finitos, proposto por Pian com base no potencial de Hellinger-Reissner, provou ser um avanço conceitual entre as formulações de discretização, tendo sido explorado extensivamente desde então por códigos acadêmicos e comerciais, também levando em conta uma série independente dos mais recentes desenvolvimentos chamados métodos de Trefftz. O método híbrido de elementos de contorno é uma generalização bem sucedida da formulação original de Pian, em que funções de Green são usadas como funções de interpolação no domínio, possibilitando assim a modelagem robusta e precisa de formas arbitrárias submetidas a vários tipos de ações.Mais recentemente, uma proposição de Przemieniecki - para a análise geral de vibração livre de elementos de treliça e viga - foi incorporada à formulação de elementos híbridos de contorno e estendida para a análise de problemas dependentes do tempo fazendo uso de um processo de superposição modal avançada que leva em conta condições iniciais gerais assim como ações de corpo gerais, além de efeitos inerciais. A presente contribuição pretende trazer para elementos finitos os melhoramentos conceituais obtidos no contexto do método híbrido de elementos de contorno. Uma grande família de macro elementos finitos híbridos é introduzida para o tratamento unificado em 2D e 3D, de problemas estáticos e transientes de elasticidade e potencial com base nas soluções fundamentais não-singulares. É também mostrado que materiais nãohomogêneos, como os novos materiais com gradação funcional, podem ser tratados consistentemente, pelo menos para problemas de potencial. Alguns exemplos numéricos simples são apresentados como ilustração dos desenvolvimentos teóricos. / [en] The hybrid finite element method, proposed by Pian on the basis of the Hellinger-Reissner potential, has proved itself a conceptual breakthrough among the discretization formulations, and has been extensively explored both academically and in commercial codes also taking into account an independent series of more recent developments called Trefftz methods. The hybrid boundary element method is a successful generalization of Pian´s original formulation, in which Green´s functions are taken as interpolation functions in the domain, thus enabling the robust and accurate modeling of arbitrarily shaped bodies submitted to several types of actions. More recently, a proposition by Przemieniecki - for the generalized free vibration analysis of truss and beam elements - was incorporated into the hybrid boundary element formulation and extended to the analysis of time-dependent problems by making use of an advanced mode superposition procedure that takes into account general initial conditions as well as general body actions, besides the inertial effect. The present contribution aims to bring to finite elements the conceptual improvements obtained in the frame of the hybrid boundary element method. A large family of hybrid, macro finite elements is introduced for the unified treatment of 2D and 3D, static and transient problems of elasticity and potential on the basis of nonsingular fundamental solutions. It is also shown that nonhomogeneous materials, as the novel functionally graded materials, may be dealt with consistently, at least for potential problems. Some simple numerical examples are shown to illustrate the theoretical developments.

Orientador: Fernando de Azevedo Silva / Banca: Ricardo Marques Barreiros / Banca: Fábio Minoru Yamaji / Resumo: Este trabalho tem o objetivo principal de compreender a metodologia de especificação e análise de linha de transmissão de torque através de cardans / juntas universais e avaliar o quanto o resultado oriundo desta metodologia influenciou na falha do suporte do mancal do cardan em uma aplicação real. Foram realizadas três linhas de trabalhos distintas, primeiro foi compreendido toda a metodologia e realizado os cálculos analíticos de acelerações e cargas de mancais nos cardans / juntas universais, em segundo foi realizado cálculos estruturais via MEF - Método dos Elementos Finitos com vários tipos de carregamentos, incluindo os obtidos através dos cálculos analíticos, e como resposta foram gerados valores de tensões normais máximas e mínimas em função de cada carregamento. Os resultados de tensões via MEF foram comparados com as tensões resultantes da terceira linha de trabalho obtidas através da instrumentação de um suporte de mancal de cardan em campo durante a aplicação real do caminhão em análise. Os resultados geram a conclusão de que os carregamentos oriundos dos cálculos analíticos possuem muito pouca influência na falha do suporte do mancal do cardan e converge para as conclusões das linhas de pesquisas já iniciadas praticamente pela totalidade dos autores do estado da arte. Uma planilha de cálculo paramétrica foi desenvolvida no intuito de efetuar a análise com base na teoria / metodologia suportanto o projetista diminuindo o tempo de aprovação, mas é necessário a continuidade de linhas de pesquisas neste assunto no intuito de prever com maior precisão o comportamento dinâmico da aplicação de linhas de transmissão de torque através de cardans e juntas universais / Abstract: Understand the methodology of specification and analysis of torque transmission line through cardan and universal joint and evaluate how much these results are influent in the gimbal bearing bracket failure in a real case is the most important objective of this job. Three distinct task lines were applied, first it was understood the entire methodology and performed all analytical calculations about accelerations and gimbal / universal joints bearing loads. In second and parallel, it was performed structural calculations through the FEA - Finite Elements Method considering several loads including the load result acquired in previous step, and the results were explicit in stress. The results of maximum and minimum stresses for all loads of the second step were compared with results of field measures stress acquired by instrumentation of a new gimbal bearing bracket with extensometers sensors during a real application of a truck under analysis. All results generate the conclusion that loads from the analytical calculations have little influence in the gimbal bearing bracket failure and converge to the conclusions of the research lines already begun for almost authors all commented in the art state. A spreadsheet was created based on the current theory to support designers during the Project, minimizing the approval time, but is necessary the studies continuity in this assumption to predict with more precision the dynamic behavior of the application of torque transmission line through gimbals and universal joints / Mestre

Ingeniero Civil / Los muros esbeltos con discontinuidades en su base, tales como muros bandera o muros con abertura tipo puerta, se dan mucho en la práctica constructiva debido a necesidades arquitectónicas en fachadas, primeros pisos o en subterráneos. Sin embargo y aunque existen datos experimentales que muestran el comportamiento de estos elementos sometidos a cargas cíclicas, estos no incluyen análisis que consideren el tamaño de la abertura como parámetro, ni recomendaciones o guías de diseño en la normativa vigente que consideren los efectos de las discontinuidades en la base. Los efectos del corte en las deformaciones unitarias en la base de un muro esbelto no son considerados en los modelos de diseño actuales, es por esto que se hace necesario estimar la magnitud de estos efectos para ser incorporados a los modelos de flexión tradicionales. En este trabajo se analizan los efectos del corte en las deformaciones de compresión máximas y tracción máximas en la base de un muro rectangular, encontrándose que la relación de aspecto tiene gran importancia en ambos casos, y la deriva de techo toma importancia en la amplificación de compresión. Se calibra un modelo matemático a través de un análisis de elementos finitos no lineal, capaz de estimar la amplificación de la compresión y reducción de la tracción en la base mediante parámetros globales del muro: largo de confinamiento, relación de aspecto, deriva plástica de techo, cuantía longitudinal de borde y el largo de promedio de deformación. También se incluye un modelo simplificado que considera la deriva plástica y la relación de aspecto como variables para determinar la amplificación de compresión máxima, y la relación de aspecto para la reducción de tracción máxima. Los modelos simplificados tienen una desviación estándar de 20% y un promedio de 1.03 en la precisión (valor obtenido sobre valor esperado) en el caso de compresión, y una desviación del 9% y promedio 0.99 en el caso de la tracción. A través de análisis de flexión no lineales, se propone un modelo que estima el desplazamiento de fluencia de muros con aberturas de borde o muros bandera, que depende de las dimensiones de la abertura y la cuantía longitudinal de borde. El modelo propuesto tiene una desviación estándar de 13% y un promedio de precisión de 1.03. Además, se proponen rangos de validez para 3 modelos que estiman la curvatura última en el rango plástico de este tipo de muros. El modelo propuesto con mayor rango de validez tiene un promedio (valor propuesto sobre valor esperado) de 1 y una desviación estándar de 29%. Finalmente se estudian 3 casos de muros con aberturas basales en el centro, describiéndose el comportamiento global y local en secciones críticas para dar luces a futura investigación y al diseño en la práctica común.

El abatimiento y control de las aguas subterráneas son actividades implementadas y monitoreadas de forma continua a lo largo de la vida de proyectos de minería, ya sea del tipo a cielo abierto o subterránea. Una implementación apropiada y eficiente de estas actividades dependen de los estudios hidrogeológicos de grande porte, que permiten evaluar los sistemas de control de agua más adecuados, necesarios en las actividades mineras. La búsqueda de la eficiencia técnica-económica de estos procesos demandan análisis numéricos de flujo tridimensionales de toda la región en estudio, caracterizada por profundas y complejos estratos de materiales permeables abajo del nivel freático, como normalmente abarcan proyectos de minería, donde las soluciones analíticas ya no pueden ser más aplicadas. Esta tesis busca contribuir a una mejor comprensión de soluciones numéricas que representan el comportamiento de flujo subterráneo a través de un estudio de caso de una mina subterránea que presenta problemas de infiltración en su interior. En este estudio fueron incorporados elementos discretos 1D como representación de caminos cársticos preferenciales, en um modelo 3D de elementos finitos realizado con el programa FEFLOW. Fueron discutidas las ventajas de incorporar tales caminos preferenciales de forma explícita, cuantificando el flujo que pasa por medio de éstas, ya que estas estructuras alimentan a la mina a través de conexiones directas con un río adyacente. Estos elementos discretos permiten representar de forma más realista el medio hidrogeológico y al mismo tiempo evaluar de mejor forma los efectos en el comportamiento de flujo subterráneo debido a la impermeabilización superficial del río como una solución más coherente para este tipo de problemas de infiltración.

Memoria para optar al título de Ingeniero Civil / En este trabajo de título se presenta el estudio numérico en base a elementos finitos de una conexión a momento tipo End-Plate , la cual consiste en una placa de acero soldada a la sección transversal de la viga, donde la unión se produce al apernar la placa a la columna solo que en este caso la placa se conectará a una plancha unida a placas de continuidad instaladas en el alma de la columna, a diferencia de una unión por eje fuerte en que la placa End-Plate es apernada a las alas y no por la zona del alma. Para llevar a cabo el diseño de dicha conexión, en el programa de análisis estructural ETABS 2016 se modeló un edificio de oficinas de 4 pisos de altura compuesto por marcos a momento ubicado en Santiago, Chile, y su diseño fue hecho de acuerdo a lo establecido por la norma NCh433 Diseño sísmico de edificios y por las especificaciones AISC 341 y AISC 360. Con los perfiles vigas y columnas obtenidos, se diseña la conexión de acuerdo a la norma AISC 358-16 Prequalified Connections for Special and Intermediate Steel Moment Frames for Seismic Applications debido a que la conexión End-Plate está precalificada para su uso y la unión propuesta pretende ser análoga a una unión por eje fuerte. Finalizado el diseño, en el programa de elementos finitos ANSYS v17.2 se lleva a cabo el estudio numérico y consiste en el análisis de 4 configuraciones de la conexión, bajo diferentes niveles de carga axial, sometidos a desplazamientos controlados cuya aplicación se hace de acuerdo al protocolo de carga establecido en AISC 341-16 con el objetivo de identificar zonas de falla y de verificar que se cumpla la filosofía de diseño empleada en el diseño de la conexión. Al analizar los resultados obtenidos es posible observar que la conexión propuesta para uniones por eje débil presenta un comportamiento semi-rígido según lo definido en la norma AISC 341-16 y además provoca una disminución en la rigidez de aquellas vigas que llegan a la columna por el eje fuerte debido a plastificaciones locales observadas en las alas y alma de la columna a medida que se incrementa el número de vigas que llega al nodo de conexión y el nivel de carga axial. Sin embargo, a pesar de que se cumple con el requisito de resistencia mínima para que la conexión pueda ser utilizada en marcos especiales a momento, se observan plastificaciones locales en los alrededores de la zona panel de la columna cuando se conecta más de una viga en la zona de unión, lo que lleva a proponer para futuras investigaciones el análisis de implementación de otros métodos de rigidización de dicha zona para evitar que la columna incursione en rango inelástico o bien llevar a cabo estudios de desempeño sísmico de otras conexiones pre-calificadas cuando son implementadas por eje débil para ampliar la gama de conexiones disponibles para este tipo de uniones y establecer con mayor claridad cuando es factible su utilización.

[pt] Em operações de perfuração, é importante ter uma estimativa do gradiente de fratura de determinada formação, como o objetivo de poder estimar o máximo peso da lama necessária para a próxima seção da formação rochosa, antes de colocar o revestimento. Por esse motivo, são realizados ensaios de Leak-Off para determinar a máxima pressão que a nova seção do poço pode sustentar sem fraturar ou perder fluido, sendo realizados durante a fase de perfuração do poço, nas formações imediatamente abaixo de cada sapata de revestimento. A finalidade deste estudo é realizar a simulação numérica deste ensaio usando um programa comercial de elementos finitos para calcular a pressão de Leak-Off, considerando a taxa de bombeamento como um dado de entrada do problema e a pressão na parede do poço como uma resposta. O trabalho abordará dois cenários: uma formação rochosa permeável e uma formação rochosa impermeável constituída de rocha de sal. Deste modo, para um estudo de uma análise real de um ensaio de Leak-Off em rocha permeável e rocha impermeável, torna-se necessário a utilização de modelos que considerem poroelasticidade e fluência, respectivamente. Para a caracterização do comportamento da rocha de sal foi usada uma lei de fluência de duplo mecanismo de deformação, que considera apenas os mecanismos de fluência: planar e indefinido. Os modelos estudados foram bidimensionais e analisados sobre a hipótese de deformação plana. A hipótese, de que a pressão de Leak-Off é atingida quando a tensão tangencial efetiva é igual a zero (calculada graficamente), foi considerada. Finalmente, foram simulados alguns casos para verificar os possíveis efeitos da alteração da permeabilidade e influência do fluido de pressurização no cálculo da pressão de Leak-Off em rocha permeável. / [en] In drilling operations, it is important to have an estimate of the fracture gradient of certain formation, in order to being able to estimate the maximum mud weight required for the next section of the rock formation, before casing is placed. For this reason, Leak-Off tests are conducted to determine the maximum pressure that the new section of the well can support without fracturing or losing fluid, being made during the drilling phase of the well in formations immediately below each casing shoe. The purpose this study is to perform the numerical simulation of this test using a commercial finite element program to calculate the Leak-Off pressure, considering the pumping rate as an input data of the problem and pressure at the borehole wall as the answer. The work will address two cases: a permeable rock formation and an impermeable rock formation composed of salt rock. Thus, for studying a real analysis of a Leak-Off test in permeable and impermeable rock, it becomes necessary to use models that consider poroelasticity and creep, respectively. For the characterization of salt rock it was used a creep law of dual mechanism of deformation, which considers only the dislocation creep mechanisms: planar and undefined. The two-dimensional models were studied and analyzed under the assumption of plane strain. The hypothesis that the Leak-Off pressure is reached when the effective tangential stress is zero (estimated graphically), was considered. Finally, some cases were simulated to check the possible effects of changing the permeability and the influence of pressurized fluid in the calculation of Leak-Off pressure in the permeable rock.

[pt] As lajes nervuradas quando comparadas às lajes maciças apresentam redução do volume de concreto, não comprometem a eficiência da estrutura e permitem a construção de lajes mais econômicas e com maiores vãos. Este trabalho tem como objetivo verificar o comportamento de três tipos de lajes nervuradas: a laje nervurada tradicional, a laje nervurada rotacionada e a laje nervurada tridirecional. Na primeira fase deste trabalho foi realizado um estudo comparativo de vários modelos usando diferentes tipos de elementos, a fim de determinar aquele que melhor representa o comportamento das lajes nervuradas. O modelo selecionado utiliza elementos de casca para representar a capa de concreto e elementos de viga para representar as nervuras, ambos com seis graus de liberdade por nó. Elementos de ligação rígida foram usados para conectar os elementos de casca e os de viga a fim de capturar a posição relativa entre a capa e as nervuras. Uma vez selecionado o modelo, foi desenvolvido um estudo dos sistemas bidirecional, rotacionado e tridirecional, no regime elástico-linear. Os resultados encontrados através do programa de elementos finitos Robot permitiram comparar os três tipos de lajes em termos dos deslocamentos obtidos no Estado Limite de Serviço, da quantidade de aço determinada através do dimensionamento no Estado Limite Último, e do volume de concreto. As lajes nervuradas tradicionais apresentaram um comportamento estrutural melhor, com lajes mais rígidas e mais econômicas sob o ponto de vista da quantidade de materiais utilizados. / [en] When compared to flat slabs, the waffle slabs reduce the volume of concrete, do not compromise the efficiency of the structure and also allow the construction of more economical slabs with longer spans. This study intends to investigate the behavior of three types of waffle slabs: a traditional waffle slab, a rotated ribbed slab and three-way slab. In the first part of this study, it was carried out a comparative analysis of various models using different types of finite element techniques in order to determine which one best represents the behavior of waffle slabs. The selected model uses shell elements to represent the concrete cover and beam elements to represent the ribs, both with six degrees of freedom per node. Rigid link elements were used to connect the shell elements and the beam elements in order to capture the relative position between the concrete cover and the ribs. Once the model was selected, a study of two-way, three-way and rotational systems were developed in the linear elastic regime. The results found using the finite element program Robot allowed the comparison between the three types of slabs in terms of, displacements obtained from the Service Limit State, amount of steel determined by designing in the Ultimate Limit State and volume of concrete. The traditional ribbed slabs had an overall better structural behavior, resulting in slabs that are more rigid and more economical from the point of view of the amount of materials used.

Orientador: Jozué Vieira Filho / Resumo: As técnicas de detecção de danos baseadas na Impedância Eletromecânica (EMI) baseiam-se na capacidade dos materiais piezoeléctricos em atuar como sensores e atuadores e contribuem para o desenvolvimento de sistemas de Structural Health Monitoring (SHM). As técnicas clássicas baseadas na EMI utilizam um transdutor Pb-Lead Zirconate Titanite (PZT ) ligado à estrutura monitorada e medem a assinatura de impedância do PZT. No entanto, as técnicas baseadas na EMI dependem de diferentes fatores, como faixa de frequência, número de PZT, temperatura ambiente, tipo de estrutura, entre outros. Assim, para demonstrar a eficácia dos métodos baseados na EMI, faz-se necessário realizar experimentos práticos, o que não é uma tarefa trivial, considerando tais fatores. Portanto, neste trabalho são estudados e propostos procedimentos para criar modelos numéricos, usando elementos finitos (FE), de técnicas baseadas na EMI usando o software PZFlex®. Além disso, os modelos desenvolvidos são usados para propor uma técnica inovadora de compensação de temperatura em sistemas baseados na EMI. Foram modeladas, simuladas e analisadas algumas estruturas clássicas como placa-de-alumínio/PZT e tubo-de-aço/PZT. Os resultados das simulações foram comparados com os equivalentes obtidos com modelos experimentais reais e mostraram-se fortemente correlacionados, indicando que o modelo proposto pode ser uma ferramenta poderosa para o desenvolvimento de técnicas de SHM baseadas na EMI. Foram realizadas simulações ... (Resumo completo, clicar acesso eletrônico abaixo) / Abstract: The damage detection techniques based on the Electromechanical Impedance (EMI) rely on the ability of piezoelectric materials in acting as sensors and actuators and contribute to the development of Structural Health Monitoring (SHM) systems. The classical EMI-based techniques use a Pb-Lead Zirconate Titanite (PZT) transducer bonded to the monitored structure and measure the impedance signature of the PZT. However, the techniques based on EMI depend on different factors, such as frequency range, number of PZTs, environmental temperature, type of structure, among others. Thus, to demonstrate the effectiveness of EMI-based methods, it is necessary to carry out practical experiments, which is not a trivial task considering such factors. Therefore, in this work, it is studied and proposed procedures to create numerical models of techniques based on the EMI, using finite elements (FE) and the PZFlex® software. In addition, the developed models are used to propose an innovative temperature compensation technique for EMI-based systems. Some classical structures were modelled, simulated and analyzed, such as aluminum plate/PZT and steel pipe/PZT. The results of the simulations were compared with the equivalents obtained by experimental models and showed to be strongly correlated, indicating that the proposed model can be a powerful tool for the development of EMI-based SHM techniques. Simulations were performed to analyze the behavior of the signatures under the effect of temperature,... (Complete abstract click electronic access below) / Doutor

[pt] Um método para a simulação dinâmica bi-dimensional de bóias cilíndricas de superfície, sujeitas à ação de ondas e ancorada ao fundo marinho através de uma linha flexível é apresentado neste trabalho. O modelo de elementos finitos adotado por Lustosa (2002) é empregado na descrição do movimento do cabo de ancoragem, onde elementos de vigas de Euler-Bernoulli de dois nós são empregados, com a cinemática do movimento descrita através de grandezas referidas ao sistema local co-rotacionado. No modelo numérico, são consideradas as não-linearidades geométricas resultantes dos grandes deslocamentos da linha e que causam o acoplamento das rigidezes axial e flexional. Na caracterização do movimento da linha consideram-se os seguintes efeitos: peso próprio, empuxo, carregamento hidrodinâmico da correnteza marinha, deslocamentos impostos pela unidade flutuante, forças de inércia e, se presente, a ação de flutuadores. Os deslocamentos do cabo são obtidos da solução do sistema algébrico de equações não-lineares, resultante da integração temporal de Newmark das equações diferenciais temporais do movimento, cuja convergência é garantida através do método iterativo de Newton-Raphson. Para a bóia, considerada para efeito de análise como indeformável, três graus-de-liberdade são considerados: dois de deslocamentos lineares e um de deslocamento angular. As ondas de superfície, monocromáticas e bi-dimensionais, são representadas no modelo através da teoria linear de Airy associada à modificação empírica de Wheeler (Wheeler, 1969). Forças resultantes da ação das ondas sobre a bóia de superfície são obtidas através da integração numérica da equação proposta por Chitrapu et al. (1998) que compõe-se de um termo viscoso - resultante da parcela de arrasto da equação de Morison - e de um termo não-viscoso - considerado como resultante da soma da força de Froude-Krylov (Chakrabarti, 1987) à parcela de inércia da equação de Morison -. A integração numérica desta equação sobre a superfície “molhada” instantânea do cilindro é realizada, a cada incremento de tempo da análise, empregando-se o método da quadratura adaptativa de Simpson. Finalmente, na integração temporal das equações de movimento da bóia de superfície utiliza-se o método de Runge-Kutta de quinta ordem, com o controle adaptável do passo temporal. O procedimento descrito foi implementado em um programa de computador e os resultados obtidos de algumas análises numéricas comparados com os fornecidos por outras simulações independentes, de forma a verificar a adequabilidade da formulação proposta no estudo da dinâmica de bóias flutuantes. / [en] In this work a method for two-dimensional dynamic analysis of a cylindrical buoy under surfaced waves and moored to the seafloor by a flexible riser is presented. The finite element model used to describe the dynamics of the maritime line was developed by Lustosa (2002). This is a two node beam element based on the Euler-Bernoulli theory, with the element kinematics referred to a co-rotacional coordinate system attached to the element. The model considers the representation of geometric non-linearities due to the line motion large displacements resulting in axial-bending couplings in the element total stiffness representation. The element displacements in the line are due to the following external loadings: self weight, buoyancy force, hydrodynamic forces due to maritime currents, surface buoy imposed displacements, the forces of inertia, and, if present, the action of floaters. The buoy two displacements and one rotation degrees-of-freedom are measured at its center of gravity and are obtained from the solution of the non-linear algebraic system resulting from the Newmark temporal integration of the differential equations of motion, which solution convergence is granted by a Newton- Rapshon iterative method. The linear wave theory - Airy Theory -, associated to the Wheeler empirical modification, is used in this work to represent the monochromatic bi-dimensional surface waves. To obtain the wave forces acting on the surface buoy while in motion, the numerical integration of the equation presented by Chitrapu et al. (1998) is employed. This is the result of a two part equation: the first is evaluated as the Morison viscous drag force and, in a second part, is the sum of the - Froude-Krylov - and the Morison inertia forces. The numerical integration of the resulting equation is evaluated over the instantaneous wet surface of the buoy, in each time-step, using the adaptive recursive Simpson´s method. Finally the temporal integration of the buoy equations of motion are evaluated using a fifth-order Runge-Kutta method, with an adaptive control of the temporal time step. This procedure has been implemented in a computer program and the results obtained from some numerical analysis are compared to the ones obtained from independent simulations, so as to verify the ability of this proposed method to describe the dynamic of floating buoys.