[pt] O método de elementos finitos é uma ferramenta amplamente utilizada na análise e simulação de fenômenos físicos. Uma etapa crucial desse método consiste na discretização de um domínio contínuo em uma malha de elementos finitos. A precisão da solução obtida na análise está diretamente relacionada á qualidade da malha utilizada, a qual é determinada pela forma e consistência de seus elementos. Grande parte dos trabalhos existentes na literatura é baseada em métodos clássicos, como Avanço de Fronteira, Triangulação de Delaunay com Restrição e Decomposição de Domínios. Utilizando uma abordagem diferente, esta dissertação propõe um método de geração a posteriori, em que restrições do domínio são inseridas iterativamente na malha, a qual é deformada localmente a fim de satisfazer cada nova restrição. Para validar o método proposto, os resultados obtidos são analisados através de critérios de qualidade comumente utilizados na literatura. / [en] The finite element method (FEM) is a widely used tool in the analysis and simulation of physical phenomena. The discretization of a continuous domain into a finite element mesh is an essential FEM step. The accuracy of the solution is highly dependent on the quality of the input mesh, which is mainly evaluated through the shape of the elements. Most of the related works is based on traditional methods, such as Advancing Front, Constrained Delaunay Triangulation and Domain Decomposition. Unlike these methods, this work proposes an iterative mesh generation method with deferred constraint insertion, in which an initially regular mesh is locally deformed in order to satisfy each new domain constraint. In addition, in order to validate the proposed method, this work evaluates each output mesh according to quality criteria commonly used in the literature.

[pt] Trata-se de um estudo sobre métodos de previsão de vida à fadiga de dutos com mossas submetidos a carregamento de pressão interna. Um método é proposto e com ele realiza-se previsões utilizando o método SN, sendo a influência das tensões médias calculadas segundo a parábola de Gerber. Adicionalmente são analisados outros métodos disponíveis na literatura, sendo abordados seus embasamentos teóricos, coeficientes de segurança e validações experimentais. A análise de tensões dos diferentes métodos foi realizada utilizando o Método de Elementos Finitos com o pacote computacional Abaqus 6.10, salvo nos casos em que foram fornecidas equações vindas de regressões de resultados numéricos. Todos os resultados avaliados foram comparados com resultados experimentais disponíveis na literatura, sendo um total de oito curvas avaliadas em seis casos em que se alteraram a espessura do tubo, o diâmetro do punção que gerou a mossa, a pressão inicial e o material do tubo. Os resultados indicam que as diferenças conceituais entre os diferentes métodos se refletem nas previsões geradas e que o método proposto faz previsões próximas aos resultados experimentais, sem o uso de fatores de segurança ou de confiabilidade. / [en] This study aims to analyze different methods to predict the fatigue life in pipelines with dents loaded with internal pressure. A new method is proposed which makes previsions based on the SN method; the mean pressure is taken into account according to the parabola of Gerber. Some other methods available in the literature are analyzed in terms of their safety factors and experimental validation. The stress analysis has been performed for the different methods using Abaqus 6.10, a Finite Element analysis software, except in those cases where a regression curve, based on numerical results, are provided. All the results found have been compared against experimental results available in the literature. A total of eight methods have been evaluated in six cases with different linepipe wall thicknesses, indentor diameters, materials and initial pressures. The results indicate that that different methods produce different life predictions, and the new method presented predict lives with a good match to the experimental data without using safety or reliability factors.

[pt] Este trabalho apresenta um estudo comparativo de métodos numéricos para solução das equações da elasticidade em problemas de otimização topológica. Um sistema computacional é desenvolvido em MATLAB para solução de problemas de otimização topológica usando malhas poligonais não estruturadas em domínios bidimensionais arbitrários. Dois métodos numéricos são implementados e comparados com o método dos elementos finitos (FEM) em relação à precisão e à eficiência computacional: o recém proposto Método dos Elementos Virtuais (VEM) e o Método dos Elementos Finitos Suavizados (SFEM). A principal característica que distingue estes métodos do FEM é que as funções de base canônicas não são obtidas de forma explícita. A utilização de projetores locais apropriados permite a extração do componente linear das deformações dos elementos e, por conseguinte, o cálculo da matriz de rigidez se reduz a avaliações de quantidades puramente geométricas. Exemplos numéricos representativos, usando malhas convexas e não convexas, para minimização da flexibilidade são apresentados para ilustrar as potencialidades dos métodos estudados. / [en] This work presents a comparative study of numerical methods for solving the elasticity equations in topology optimization problems. A computational framework is developed in MATLAB for solving topology optimization problems using unstructured polygonal meshes in arbitrary two-dimensional domains. Two numerical methods are implemented and compared with the finite element method (FEM) with respect to accuracy and computational efficiency: the recentlyproposed Virtual Element Method (VEM) and the Smoothed Finite Element Method (SFEM). The key characteristic that distinguish these methods from the FEM is that the canonical basis functions are not computed explicitly. The use of appropriate local projection maps allows the extraction of the linear component of the element deformations and, therefore, the computation of the stiffness matrix is reduced to the evaluation of purely geometric quantities. Representative numerical examples, using convex and non-convex meshes, for compliance minimization are presented to illustrate the capabilities of the methods studied.

[pt] O processo de revestimento de substratos é uma das principais etapas na manufatura de diversos produtos, tais como eletrônicos transparentes, fitas magnéticas e adesivos. Em sua grande maioria, os estudos disponíveis na literatura consideram que o líquido de revestimento é um fluido Newtoniano. No entanto, em grande parte das aplicações, este fluido é mais complexo, sendo composto por partículas dispersas em uma fase contínua. Nesse sentido, efeitos reológicos relacionados à presença das partículas no escoamento ainda não são bem estabelecidos. Por exemplo, variações nos valores locais de taxa de cisalhamento e concentração induzem à diferentes mecanismos de migração de partículas que, por sua vez, podem mudar os padrões do escoamento e a estrutura final do filme. Além disso, se o tamanho das partículas é pequeno o suficiente de modo que forças coloidais são relevantes, fenômenos locais como redução ou aumento de viscosidade podem ocorrer dependendo das contribuições macroscópicas do escoamento. Neste trabalho, as equações de conservação de massa, momento e transporte de partículas são consideradas juntas com um modelo de viscosidade sensível a variações locais tanto de concentração quanto da taxa de cisalhamento. Nesse sentido, o Modelo de Fluxo Difusivo é empregado com o intuito de capturar fluxos difusivos de partículas em casos nos quais interações irreversíveis entre partículas são relevantes. As equaçõoes de conservação são acopladas a um sistema de equações diferencias que descreve o mapeamento numérico entre o domínio físico e outro de referência, e são discretizadas usando o Método de Elementos Finitos com a formulação fraca de Galerkin. O sistema não-linear resultante é resolvido por meio do método de Newton associado à um robusto método LU de decomposição matricial. Os resultados mostram que o comportamento reológico de suspensões coloidais afeta a distribuição final de partículas no filme depositado e os limites de operação de um processo típico de revestimento de substrato por extrusão. / [en] Coating is the main step in the manufacturing process of many common and new products such as paper, adhesive and magnetic tapes, and flexible and transparent electronics. Most of the available analyses consider the coating liquid as Newtonian. However, in many applications the coating liquid is a complex fluid composed by particles suspended in a viscous medium. The effect of suspended particles in the coating liquid on the coating flow remains unknown. Local variation of particle concentration associated with different mechanisms of particle migration may lead to viscosity gradients, which can change the flow pattern and final particle structure in the coated film. Moreover, if the particles are small enough, colloidal forces are important and hence non-Newtonian behavior such as shear-thinning and shear-thickening may occur. In this study, the conservation equations of mass, momentum and particle transport, are considered in which the viscosity is assumed to vary not only with local particle concentration but also with local shear rate (i.e. shear sensitive liquids). With this in mind, we employ the Diffusive Flux Model (DFM) to capture particle diffusion where irreversible particle-particle interactions are dominant. The conservation equations together with the set of differential equations that describes the mapping between the unknown flow domain to a fixed reference domain are solved using the Galerkin/Finite Element Method (GFEM). The resulting non-linear system of algebraic equations is solved using Newton s method coupled with a robust LU frontal solver. The results show how the rheological behavior of colloidal suspensions affects the particle distribution on the deposited liquid film and the process limits in a typical slot coating flow.

[pt] Após a década de 80 tornou-se frequente a utilização de poços direcionais, tanto poços horizontais quanto os de longo alcance, nas operações de desenvolvimentos de campos petrolíferos, os quais permitiram grande avanço na exploração. Este cenário, porém, traz grandes desafios para operações de perfuração relacionadas à limpeza de poços. A limpeza de poços consiste na remoção dos cascalhos de dentro do anular do poço através da circulação do fluido de perfuração. Devido à inclinação do poço, os cascalhos sofrem a ação da gravidade sendo empurrados para o fundo do canal, já que possuem densidade maior que a da fase líquida, assim, transportados a uma velocidade menor que a velocidade do escoamento. Ao se depositarem no fundo do canal, formam um leito estacionário, obstruindo parte do anular e diminuindo a vazão, o que gera problemas, como por exemplo, redução da taxa de penetração, desgaste prematuro da broca, elevação do torque e arraste, aprisionamento da coluna, perda de circulação, dentre outros, podendo gerar a perda do poço. Desta forma, o perfeito entendimento do processo de sedimentação e transporte de partículas sólidas suspensas em fluido é fundamental para a otimização do processo de perfuração de poços. Este trabalho analisa o escoamento laminar bidimensional de suspensões de partículas sólidas devido a um gradiente de pressão entre duas placas paralelas, representando uma descrição simplificada do escoamento que ocorre em um anular de poço durante o processo de perfuração. O perfil do leito de partículas ao longo do canal e a relação vazão-diferença de pressão para diferentes condições de escoamento são determinados pela solução numérica das equações que descrevem o problema. A formulação matemática leva a um sistema acoplado de três equações diferenciais: conservação de massa e de quantidade de movimento e a equação de transporte, que engloba os efeitos de difusão de partículas devido à frequência de interação entre as partículas, ao gradiente de viscosidade e à diferença de densidade entre o líquido e as partículas. O sistema é resolvido pela técnica de Elementos Finitos, através do método de Galerkin. Os resultados obtidos serão de extrema importância no desenvolvimento de modelos mais precisos que descrevam o processo de transporte de sólidos em anulares de poços. / [en] After the 80 s, the use of directional wells, both horizontal and long range wells, became frequent in the development of oil fields, which allowed great progress in exploration. This scenario, however, brings great challenges to operations related to wellbore cleaning. Wellbore cleaning consists in the removal of cuttings from within the annular through the circulation of drilling fluid. Due to the inclination of the well, the cuttings undergo the action of gravity and more to the bottom of the channel, as they have higher density than the liquid phase. They may be transported with a speed less than the speed of the liquid flow. When deposited on the channel, the cuttings form a stationary bed, blocking part of the annular and decreasing the flow rate that causes problems, such as reducing the rate of penetration, premature wear of the bit, high torque and drag, trapping column, loss of circulation and others, which may cause the loss of the well. Therefore, the fundamental understanding of particle sedimentation and transport in a suspending flowing liquid is necessary for drilling operation design and optimization. This research studies the laminar two-dimensional flow of solid particles suspended in a liquid due to a pressure gradient between two parallel plates, representing a simplified description of the flow that occurs in an annular during the drilling process. The profile of the particle bed along the channel and the flow rate pressure difference relationship for different flow conditions are determined. The mathematical formulation leads to a coupled system of three differential equations: mass and momentum conservation and transport equation, which includes the effects of particle diffusion due to the frequency of interaction between the particles, the gradient of viscosity and the difference density between the liquid and the particles. The system is solved by the finite element method Galerkin. The results will be of extreme importance in the development of more accurate models that describe the solids transport process in annular space.

[pt] As zonas de falha são compostas por um núcleo, onde a maior parte da deformação é acomodada, e uma zona de dano, com deformação menos intensa. A zona de dano pode atuar como caminho de fluxo preferencial devido à presença de fraturas, ou como barreira devido às bandas de deformação. Portanto, sua caracterização é essencial para a adoção de estratégias de produção adequadas em campos de petróleo. Os métodos geofísicos geralmente utilizados, porém, dificilmente permitem a identificação das zonas de dano devido à baixa resolução sísmica. Como alternativa, empregam-se observações em afloramentos superficiais. Contudo, há uma grande dispersão de dados, que pode estar relacionada a uma variedade de fatores, como as propriedades da rocha protólita e os mecanismos de deformação atuantes. Neste sentido, este trabalho apresenta duas metodologias baseadas no método dos elementos finitos (MEF) para analisar a formação e evolução das zonas de dano em escala de reservatório. Na primeira abordagem, a zona de falha é totalmente representada através de um meio contínuo, enquanto que na segunda, a falha é representada como um plano por meio de uma descontinuidade. Em ambas aproximações, a zona de dano é estabelecida através das regiões plastificadas. Os resultados numéricos obtidos se aproximaram das observações de campo e possibilitaram a identificação das vantagens e limitações das duas abordagens baseadas no MEF. Por fim, os resultados também permitiram identificar os principais parâmetros geomecânicos que influenciam o desenvolvimento das zonas de dano, bem como os diferentes mecanismos de deformação que ocorrem ao longo da zona de dano. / [en] Fault zones are composed of two structural domains: the fault core, which accommodates most of the deformation, and a damage zone, with less intense deformation. The damage zone may act as a preferential flow path due to the presence of fractures, or as a barrier due to deformation bands. Therefore, the characterization of geological fault zones is essential for the adoption of adequate production strategies in oil fields. Generally, geophysical methods are used to characterize geological faults in the field. However, they hardly allow the identification of damage zones due to low seismic resolution. Alternatively, damage zones are analyzed through surface outcrops. Nonetheless, there is a wide dispersion of data in this type of study, which may be related to various factors, such as the properties of the host rock and the acting deformation mechanisms. Therefore, it is interesting to carry out this type of study in conjunction with numerical modeling to understand better the damage zone formation process. In this study, we present two methodologies based on the finite element method (FEM) to analyze the formation and evolution of damage zones at a reservoir scale. In the first methodology, the entire fault zone is represented through a continuum medium, while in the second methodology, the fault core is represented as a plane through a discontinuity. In both approaches, the damage zone is defined through the regions where plastic deformations were triggered. The numerical results obtained were close to field observations. They enabled the identification of the advantages and limitations of the two approaches based on the MEF. Finally, the results also allowed to identify the main parameters that influence the development of the damage zones and the different deformation mechanisms that occur along the damage zone.

[pt] Compósitos cimentícios estão ganhando cada vez mais relevância na indústria da construção civil. No entanto, as diretrizes para o projeto do material compósito e dos seus elementos estruturais são ainda incipientes, pois mecanismos de ponte de transferência de forças providos pelas fibras ainda estão sob investigação. Este trabalho apresenta uma estratégia de modelagem de elementos finitos que leva em consideração a estrutura de nível mesoestrutural do material cimentício reforçado com fibras. Desta forma, quatro fases do material são consideradas no modelo numérico: agregados graúdos, argamassa, zona de transição interfacial (ZTI) e fibras. A argamassa e os agregados são modelados usando elementos contínuos triangulares com comportamento linear-elástico. As fibras são incluídas usando elementos de treliça unidimensionais acopladas a elementos bidimensionais contínuos. Uma técnica de fragmentação de malha é usada para introduzir elementos de interface nas arestas dos elementos de argamassa e na interface entre agregados e argamassa para representar a ZTI. O método Take-and-Place, proposto por Wriggers e Moftah (2006), foi adotado neste estudo para incluir agregados no modelo. Primeiro, os agregados são gerados seguindo uma curva de Fuller, que define um empacotamento entre os agregados perfeitos. Na segunda fase, os agregados são introduzidos no modelo garantindo a não sobreposição entre eles. Finalmente, as fibras são adicionadas. Para validar a metodologia proposta, testes experimentais foram simulados com sucesso em um framework de simulação numérica – GeMA. Por fim, o trabalho explora a influência do empacotamento fibra-agregado na resposta mecânica e nos padrões de fraturamento de compósitos cimentícios fibrosos. / [en] Fiber Reinforced Concrete (FRC) materials are gaining more relevance in the construction industry. However, the guidelines for the design of the composite material and of structural elements thereof are incipient and the stress bridging mechanisms are still under investigation. This work presents a finite element modelling strategy that takes into account the material meso-level structure. Four phases of the FRC material are considered in the model: coarse aggregates, mortar, interfacial transition zone (ITZ), and fibers. The mortar and aggregates are modelled using triangular linear elements with linear–elastic behavior. Fibers are included using one-dimensional truss elements which are coupled to the matrix through the technique proposed by Congro (2021). Zero-thickness interface elements are introduced at the interface between mortar elements, and at the interface between aggregates and mortar to represent the ITZ. The Take-and-Place method, obtained from Wriggers and Moftah (2006), was adopted in this study to include aggregates in the model. First, the aggregates are generated following a Fuller s curve that means a perfect aggregate package. In the second phase, the aggregates are placed in the model without overlapping. Finally, fibers were added. A mesh fragmentation technique is used to introduced zero-thickness interface elements at the interface between mortar elements, and at the interface between aggregates and mortar to represent the ITZ. To validate the proposed methodology, direct tensile test models were successfully reproduced in finite element analyses performed in an in-house framework – GeMA. Based on the obtained results, the authors could explore the influence of the fibers-aggregate packing in the mechanical response of the composite material.

[pt] A interação entre o cimento e outros constituintes desempenha um papel importante em várias aplicações de Engenharia, como nas indústrias de construçãocivil e de óleo e gás (OeG). Na indústria da construção, os compósitos cimentícios reforçados com fibras (CRF) ganharam grande destaque por suas excelentes propriedades mecânicas. As fibras podem aumentar a resistência crítica à fissuração do compósito, melhorando a durabilidade do concreto convencional e controlando a propagação de fissuras na matriz cimentícia. Além disso, as fibras desenvolvem um mecanismo de ponte de transferência de tensões na interface, alterando o comportamento pós-pico do compósito. Por outro lado, na indústria de OeG, cimento e aço são elementos estruturais essenciais que devem garantir a integridade de poços e fornecer isolamento para a passagem de fluidos, especialmente em cenários de abandono. Esse mecanismo na interface é considerado crítico, uma vez que uma interação não eficaz pode permitir a formação de caminhos de vazamento no microanular ao longo da interface cimento-aço, gerando a formação de fissuras. Neste sentido, um estudo abrangente dos mecanismos de dano desenvolvidos na interface do cimento é essencial em ambas as aplicações para entender o comportamento mecânico do material. Portanto, faz-se necessário o desenvolvimento de modelos de elementos finitos que considerem os mecanismos de pullout (descolamento, adesão e atrito) e os parâmetros de interface que governam o comportamento mecânico local do cimento. Embora existam numerosos estudos experimentais e modelos numéricos na literatura, o estado-da-arte atual carece de formulações que investiguem os mecanismos de mapeamento de dano e as interações de transferência de tensão na interface do cimento, especialmente considerando diferentes tipos de matriz de cimento e geometrias de fibra de aço.Esta tese aborda uma lacuna crítica na literatura ao propor a modelagem numérica do descolamento interfacial e mecanismos de evolução de dano para materiais cimentícios avançados e em aplicações de integridade de poços. Modelos de elementos finitos elastoplásticos, incorporando formulações coesivas baseadas em superfícies de contato, são empregados para simular o comportamento da interface do cimento. Além disso, ensaios experimentais de caracterização mecânica e análises de microtomografia são realizados para validar e apoiar os resultados do modelo numérico, avaliando a resistência ao cisalhamento e a propagação de dano na interface do cimento. Assim sendo, esta pesquisa pode oferecer contribuições para engenheiros de diferentes áreas aprimorarem o desempenho mecânico e prototipar novos materiais avançados por meio da investigação da evolução do dano. Os modelos de elementos finitos desenvolvidos emergem como ferramentas valiosas para avaliações de desempenho do cimento de maneira eficaz, simulando confiavelmente o comportamento de pullout/pushout. / [en] The interaction between cement and other constituents plays an important role in several engineering applications, such as in the construction and oil and gas (OandG) industries. In the construction industry, fiber-reinforced cementitious composites (FRC) have gained wide prominence for their excellent mechanical properties. Fibers can increase the post-cracking strength of the composite, improving concrete durability and controlling crack propagation in the cement matrix. Moreover, they perform a bridging mechanism at the interface, changing the material post-peak behavior. On the other hand, in the OandG industry, cement and steel are essential structural elements that should ensure well integrity and provide zonal isolation. This interaction is considered critical since a strong bond may prevent the generation of microannulus leakage paths along the cement and steel interface, which also can lead to crack propagation. In this sense, a comprehensive study of the damage mechanisms developed at the cement interface is essential in both applications to understand the material mechanical behavior. Therefore, it is possible to develop finite element models that consider the pullout mechanisms (debonding, adhesion, and friction) and the interface parameters that govern the local mechanical behavior of cement. While numerous experimental studies and numerical models exist, the current state-of-the-art lacks formulations investigating damage mapping and stress transfer interactions at the cement interface, particularly considering different cement matrix types and steel fiber geometries. This thesis addresses a critical gap in the literature by proposing the numerical modeling of interfacial debonding and damage evolution mechanisms for cement advanced materials and well integrity applications. Elastoplastic finite element models, incorporating surface-based cohesive formulations with contact, are employed to simulate cement interface behavior. Additionally, mechanical characterization tests and microCT analyses are conducted to validate and support the numerical model results, assessing shear strength and damage propagation at the cement interface. Therefore, this research can offer insights for engineers across disciplines to enhance mechanical performance and prototype new advanced materials by damage evolution investigation. The developed finite element models emerge as valuable tools for cost-effective evaluations of cement performance through reliably simulating pullout/pushout behavior.

[pt] O dimensionamento de estrutura de contenção flexíveis atirantadas geralmente representa um desafio para os engenheiros civis, devido a complexidade do tema e as severas consequências que uma eventual ruptura pode causar. Visando melhorar a compreensão do comportamento dessas estruturas, bem como o seu dimensionamento, o presente trabalho descreve os mecanismos que governam o comportamento dessas contenções, como por exemplo: o arqueamento das tensões; a redistribuição das tensões; concentração de tensões nas ancoragens, etc. Embora estes mecanismos tenha uma forte influência sobre o comportamento de cortinas atirantadas, eles geralmente não são considerados em projetos, já que na prática brasileira, os dimensionamentos dessas contenções são realizados predominante-mente a partir de métodos de cálculo simples, como por exemplo o cálculo por equilíbrio limite, que não leva em conta esses efeitos. Portanto, no presente trabalho são discutidas as limitações dessa abordagem, e em contrapartida são apresentadas, através de um estudo de caso, os benefícios da aplicação do método dos elementos finitos em análises de escavações. / [en] The design of anchored retaining wall usually represents a challenge for civil engineers, due to the complexity of the issue and the severe consequences that a possible rupture can cause. In order to improve understanding of the behavior of these structures, as well as its design, the present work describes the mechanisms that govern the behavior of these structures, for example: the arching of the stress, the redistribution of stress, concentration of stresses in the anchorages, etc. Although these mechanisms have a strong influence on behavior of the anchored retaining wall, they are usually not considered in design, because the Brazilian practice, the projects of these structures are made predominantly from simple calculation methods, for example by calculating the equilibrium limit, that does not account for these effects. Therefore, in the present work the limitations of this approach are discussed, and on the other hand are presented, via a case study, the benefits of the finite element method applied to the excavation analysis.

[pt] Apresentamos um método computacional na GPU que lida com eventos de fragmentação dinâmica, simulados por meio de zona coesiva. Implementamos uma estrutura de dados topológica simples e especializada para malhas com triângulos ou tetraedros, projetada para rodar eficientemente e minimizar ocupação de memória na GPU. Apresentamos um código dinâmico paralelo, adaptativo e distribuído que implementa a formulação de modelo zona coesiva extrínsica (CZM), onde elementos são inseridos adaptativamente, onde e quando necessários. O principal objetivo na implementação deste framework computacional reside na habilidade de adaptar a malha de forma dinâmica e consistente, inserindo elementos coesivos nas facetas fraturadas e inserindo e removendo elementos e nós no caso da malha adaptativa. Apresentamos estratégias para refinar e simplificar a malha para lidar com simulações dinâmicas de malhas adaptativas na GPU. Utilizamos uma versão de escala reduzida do nosso modelo para demonstrar o impacto da variação de operações de ponto flutuante no padrão final de fratura. Uma nova estratégia de duplicar nós conhecidos como ghosts também é apresentado quando distribuindo a simulação em diversas partições de um cluster. Deste modo, resultados das simulações paralelas apresentam um ganho de desempenho ao adotar estratégias como distribuir trabalhos entre threads para o mesmo elemento e lançar vários threads por elemento. Para evitar concorrência ao acessar entidades compartilhadas, aplicamos a coloração de grafo para malhas não-adaptativas e percorrimento nodal no caso adaptativo. Experimentos demonstram que a eficiência da GPU aumenta com o número de nós e elementos da malha. / [en] A GPU-based computational framework is presented to deal with dynamic failure events simulated by means of cohesive zone elements. We employ a novel and simplified topological data structure relative to CPU implementation and specialized for meshes with triangles or tetrahedra, designed to run efficiently and minimize memory requirements on the GPU. We present a parallel, adaptive and distributed explicit dynamics code that implements an extrinsic cohesive zone formulation where the elements are inserted on-the-fly, when needed and where needed. The main challenge for implementing a GPU-based computational framework using an extrinsic cohesive zone formulation resides on being able to dynamically adapt the mesh, in a consistent way, by inserting cohesive elements on fractured facets and inserting or removing bulk elements and nodes in the adaptive mesh modification case. We present a strategy to refine and coarsen the mesh to handle dynamic mesh modification simulations on the GPU. We use a reduced scale version of the experimental specimen in the adaptive fracture simulations to demonstrate the impact of variation in floating point operations on the final fracture pattern. A novel strategy to duplicate ghost nodes when distributing the simulation in different compute nodes containing one GPU each is also presented. Results from parallel simulations show an increase in performance when adopting strategies such as distributing different jobs amongst threads for the same element and launching many threads per element. To avoid concurrency on accessing shared entities, we employ graph coloring for non-adaptive meshes and node traversal for the adaptive case. Experiments show that GPU efficiency increases with the number of nodes and bulk elements.