Made available in DSpace on 2014-06-11T19:35:03Z (GMT). No. of bitstreams: 0 Previous issue date: 2008-10-31Bitstream added on 2014-06-13T20:06:16Z : No. of bitstreams: 1 archangelo_cm_dr_araca.pdf: 3340208 bytes, checksum: 53b104f834f13f3fac3b5838d3d64527 (MD5) / Objetivo: O objetivo do presente estudo foi avaliar através do método dos elementos finitos tridimensional (3D FEA), o comportamento mecânico de uma faceta laminada de cerâmica e do cimento resionoso em um 1o pré-molar maxilar, variando a redução (parcial ou total) da cúspide vestibular. Material e método: Baseado em imagens micro-tomográficas (μCT) 2 modelos (M) tridimensionais foram construídos no programa Solidworks 2007 variando a redução da cúspide vestibular, sendo Mp o modelo composto pelo dente pré-molar restaurado com uma faceta laminada de cerâmica, com redução parcial da cúspide vestibular e Mt, o modelo semelhante ao Mp, porém com redução total da cúspide vestibular. O programa de elementos finitos ANSYS Workbench 10.0 foi utilizado para realizar as análises. O carregamento (150 N) foi aplicado próximo à ponta de cúspide vestibular, com 45º de inclinação. Foram utilizados os critérios da máxima tensão principal ( max) e da máxima tensão cisalhante ( shear) na faceta laminada e na camada de cimento resinoso. Resultados: O Mt mostrou menor shear em 9 das 12 áreas analisadas quando comparado ao Mp. O Mp mostrou menor variação da max (de -3,13 para 14,4 MPa) que o Mt para a lâmina de cerâmica (de -2,08 para 24,3 MPa) A camada de cimento apresentou variação de max semelhante entre os Mp e Mt. Os critérios de tensões utilizados não evidenciaram valores além dos limites estabelecidos para a falha estrutural do material estudado. Conclusão: O laminado cerâmico e a camada de cimento resinoso apresentaram comportamento mecânico favorável para a restauração do pré-molar frente ao carregamento adotado. O tipo de redução de cúspide não influenciou no aumento expressivo das tensões principal e cisalhante. O Mp apresentou menores valores de max em tração quando comparado ao Mt. O Mp e o Mt apresentaram comportamento favorável com relação as forças cisalhantes, sendo as maiores tensões. / Objectives: Using 3-D finite element analysis (FEA) and μCT data, the aim of this study was to evaluate the behavior of the laminate veneer and the resin luting cement in a maxillary premolar with the buccal cusp partially or completed covered. Material and method: 2 3-D FEA models of a maxillary premolar were built from μCT data. All models were restored with a feldspathic laminate veneer using Solidworks. A similar tooth reduction (0.5mm) was applied for two models. The buccal cusp reduction of the premolar followed two configurations: Model (Mt) - buccal cusp completely covered by veneer; and Model (Mp) - buccal cusp partially covered by veneer. The loading (150N, in 45o) was performed on the buccal cusp of each premolar. Appropriate boundary condition was applied at the periodontal ligament surfaces and ANSYS (10.0) was used to obtain the maximum shear stress ( shear) and maximum principal stress ( max) for the laminate veneer and resin luting cement layer. Results: The Mp showed higher shear than Mt in 9 to 12 analyzed areas. The Mp showed smallest variation of max in laminate veneer (from -3,13 to 14,4 MPa) compared with Mt (from - 2,08 to 24,3 MPa). The resin luting cement showed similar variation of max between Mp and Mt (from -3.13 to 12.2 MPa) and (from -1.61 to 14.1 MPa), respectively. None of the structures exceeded the ultimate tensile strength or the shear bond strength to dentin and enamel. Conclusion: The laminate veneer and the resin luting cement layer showed an adequate behavior when veneer was applied in premolars. Cusp reduction did not affect significant increase in max and shear. In general, the Mp showed better behavior ( max) than Mt. The Mt showed better behavior ( shear) than Mp.

In this dissertation, we study the stationary heat equation in a heterogeneous tridimensional plate, using a "dimension reduction" techinique called hierarchical modelling and we generate model the original problem in a two-dimensional domain. To estimate the error modelling, we develop an asymptotic expansion for the original problem solution and for the aproximate solution. Comparing both solutions with their own asymptotic expansions, we obtain an estimative of the error modelling. We perform some computational experiments, using the Residual Free Bubbles (RFB) Method and the Multiscale Finite Element Method for the diffusion problem and for the diffusion-reaction problem in a two-dimensional domain, with small parameters. Finally, we extend the numeric solutions found the original tridimensional problem. / Neste trabalho, estudamos a equação do calor estacionária em uma placa heterogênea tridimensional. Para a modelagem deste problema, utilizamos uma técnica de redução de dimensão conhecida por Modelagem Hierárquica. Desta forma, geramos um modelo para o problema original em um domínio bidimensional. Com o objetivo de estimar o erro de modelagem, desenvolvemos a expansão assintótica da solução do problema original e da solução aproximada. Comparando as soluções com suas respectivas expansões assintóticas, obtemos uma estimativa para o erro de modelagem. Realizamos alguns experimentos computacionais, desenvolvendo o método Residual Free Bubbles (RFB) e o método de Elementos Finitos Multiescala (MEFM) para o problema de difusão e para o problema de difusão-reação em um domínio bidimensional, com parâmetros pequenos. Com base nestes experimentos, encontramos algumas soluções numéricas para o problema da placa tridimensional.

O presente trabalho tem por finalidade estudar e implementar um modelo numérico de análises cinemáticas de sólidos tridimensionais via método dos elementos finitos posicionais, com consideração de fibras longas ou curtas inseridas de maneira aleatória ou não no domínio da análise. O modelo numérico considera material isotrópico para a matriz e não linearidade geométrica. O domínio do sólido é discretizado por meio de elementos finitos tetraédricos de ordem qualquer, cujos parâmetros nodais são suas posições. A medida de deformação utilizada é a de Green, associada à lei constitutiva de Saint-Venant-Kirchhoff, referenciada pela configuração inicial do corpo, caracterizando o sistema de espaço como Lagrangiano total. O cálculo da posição de equilíbrio é baseado no princípio da mínima energia potencial total. Para a resolução do problema não linear geométrico, adota-se o método iterativo de Newton-Raphson. A inserção das fibras no domínio da análise é feita com a associação das mesmas com elementos finitos unidimensionais curvos de ordem qualquer, cujas posições nodais são dadas em função das posições dos nós dos elementos de sólido. Essa abordagem tem como vantagem o fato de não aumentar o número de graus de liberdade do sistema, ao mesmo tempo em que não limita as posições das fibras dentro do domínio por não ser necessária a coincidência das malhas. Exemplos são apresentados para validação dos desenvolvimentos e implementações realizadas. / This study aims to develop and implement a numerical model of kinematic enrichment, to analyze tridimensional solids based on positional finite element method, considering long and short fibers random distributed inside the domain. The numerical model considers isotropic material and geometric nonlinear behavior for both matrix and fibers. Tetrahedral finite elements with any order of approximation are used to discretize the solid domain, with positions as nodal parameters. Green strain and Saint-Venant-Kirchhoff constitutive law are used, referenced in initial configuration of the body, characterizing the developed formulation as total Lagrangian. The equilibrium is obtained with the application of Total Potential Energy Principle, adopting the Newton-Raphson method to solve the resulting nonlinear system of equations. The fibers are considered in the formulation using curved one-dimensional finite elements with any order of approximation, and the nodal positions of the fibers are related with the nodal positions of the solid elements. The coupling method adopted does not increase the number of degrees of freedom of the system, and does not limit the positions of the fiber nodes to be coincident with solid nodes. Examples are presented in order to validate the developed and implemented formulations.

A análise de compósitos laminados apresenta grandes desafios, pois, diferentemente dos materiais isotrópicos homogêneos, os compósitos laminados são constituídos de materiais heterogêneos e anisotrópicos. Além disso, as distribuições de tensões interlaminares obtidas com as formulações convencionais são descontínuas e imprecisas. Sua melhoria, portanto, é imprescindível para buscar e modelar critérios de falha relacionados às estruturas formadas por compósitos laminados. Diante disso, este trabalho se concentrou no desenvolvimento e implementação computacional de um elemento finito posicional de pórtico plano laminado cuja cinemática é descrita ao longo da espessura do laminado de acordo com a teoria Layerwise. A formulação do elemento considera a não linearidade geométrica, originada pela ocorrência de grandes deslocamentos e rotações, e admite deformações moderadas, em função da lei constitutiva de Saint-Venant-Kirchhoff. O desenvolvimento deste trabalho se iniciou com uma preparação teórica sobre mecânica dos sólidos deformáveis e métodos numéricos para que fossem adquiridos os subsídios teóricos necessários ao desenvolvimento de códigos computacionais, à interpretação dos resultados e à tomada de decisões quando das análises numéricas. A formulação desenvolvida é Lagrangiana total com emprego do método dos elementos finitos baseado em posições. Inicialmente o elemento finito posicional de pórtico plano homogêneo é proposto, uma vez que sua cinemática possibilita uma expansão natural para o caso laminado. Os graus de liberdade são compostos por posições nodais e por vetores generalizados que representam o giro e a variação na altura da seção transversal. A eficiência do elemento é constatada através de análises realizadas em problemas de pórtico sujeitos a grandes deslocamentos e rotações. Os resultados obtidos apresentaram excelente concordância com soluções numéricas e analíticas disponíveis na literatura. Uma expansão natural da cinemática é empregada na formulação do elemento laminado. Os graus de liberdade do elemento são as posições nodais e as componentes de vetores generalizados associados às seções transversais de cada lâmina. Dessa forma, as lâminas têm liberdade para variação de espessura e giro independente das demais, mas com as posições compatibilizadas nas interfaces. Os resultados de análises numéricas realizadas em vários exemplos demonstram a eficiência da formulação proposta, pois as distribuições de deslocamentos e tensões ao longo da espessura do laminado apresentaram excelente concordância com as obtidas a partir de análises numéricas utilizando um elemento finito bidimensional em uma discretização bastante refinada. Os exemplos analisados contemplam problemas com seção laminada fina ou espessa. / The analysis of laminated composites presents challenges because, unlike homogeneous isotropic materials, the laminated composites are made up of heterogeneous and anisotropic materials. Moreover, the distribution of interlaminar stresses obtained with conventional formulations are discontinuous and inaccurate. His improvement is therefore essential to check and modeling failure criteria related to structures formed by laminates. Thus, this work focused on developing and computational implementation of a positional finite element of laminated plane frame whose kinematics is described throughout the thickness of the laminate according to Layerwise theory. The formulation element considers the geometric nonlinearity, caused by the occurrence of large displacements and rotations, and admits moderate deformation, in the constitutive law function of Saint-Venant-Kirchhoff. The development of this work began with a theoretical preparation on mechanics of deformable solids and numerical methods for the acquired of the theoretical support needed for the development of computational codes, interpretation of results and decision-making when of the numerical analyzes. The developed formulation is total Lagrangian with use of the finite element method based on positions. Initially the positional finite element of homogeneous plane frame is proposed, since their kinematic enables a natural expansion for the laminate case. The degrees of freedom are composed of nodal positions and generalized vectors representing the spin and the variation in the height of the cross section. The efficiency of the element is verified through analyzes performed in frame problems subject to large displacements and rotations. The results showed excellent agreement with numerical and analytical solutions available in the literature. A natural expansion of the kinematics is used in the formulation of the laminate element. The degrees of freedom of the element are the nodal positions and components of the generalized vectors associated to cross-sections of each lamina. Thus, the laminas are free for the thickness variation and for independent spin, but with the positions matched in the interfaces. The results of numerical analysis performed in various examples show the effectiveness of the proposed formulation, since the distributions of displacements and stresses through the thickness of the laminate agreed well with those obtained from numerical analysis using a discretization with two-dimensional finite elements in a very refined. The examples discussed include problems with thin or thick laminated section.

Submitted by Renata Lopes (renatasil82@gmail.com) on 2016-07-04T21:05:03Z No. of bitstreams: 1 fernandoluizgoulartcruz.pdf: 3319031 bytes, checksum: b8ba4abba29514e3926b4197c93cdcf2 (MD5) / Approved for entry into archive by Adriana Oliveira (adriana.oliveira@ufjf.edu.br) on 2016-07-13T16:23:44Z (GMT) No. of bitstreams: 1 fernandoluizgoulartcruz.pdf: 3319031 bytes, checksum: b8ba4abba29514e3926b4197c93cdcf2 (MD5) / Made available in DSpace on 2016-07-13T16:23:45Z (GMT). No. of bitstreams: 1 fernandoluizgoulartcruz.pdf: 3319031 bytes, checksum: b8ba4abba29514e3926b4197c93cdcf2 (MD5) Previous issue date: 2012-03-06 / CAPES - Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / O objetivo deste estudo foi avaliar pelo método de elementos finitos (MEF), a influência da orientação do Mini-implante Ortodôntico (MIO), em relação ao carregamento, na distribuição de tensões no tecido ósseo. Por meio da interface gráfica do programa de elementos finitos Abaqus/CAE, foram confeccionados modelos tridimensionais, correspondentes a um MIO, a uma seção do osso maxilar e dentes pré-molar e molar provenientes de cortes de tomografia computadorizada no formato Dicom. Todas as partes foram consideradas como seções sólidas e homogêneas, com propriedades isotrópicas e lineares-elásticas. Foi aplicada no MIO uma força de 2N, decomposta a 45º em duas componentes sobre o plano normal ao seu eixo longitudinal. Para a simulação numérica, o MIO foi instalado em nove orientações proximodistais diferentes. Uma instalação ortogonal à superfície vestibular (angulação de 0º), outras de -10º, -20º, -30º e -40º a favor da força aplicada e 10º, 20º, 30º e 40º contrários à força aplicada. Por meio do critério de von Mises e da análise das tensões principais observou-se que houve uma variação das tensões em função da orientação, sendo que a situação que ocorreu menor concentração para o osso cortical foi a de 0º, enquanto que a -20º e 20º ocorreram maiores concentrações. Conclui-se que houve influência da orientação proximodistal do MIO, em relação ao carregamento, na distribuição de tensões no tecido ósseo. / The aim of this study was to evaluate, by finite element method (FEM), the influence of load direction on stress distribution in bone tissue around orthodontic miniscrew. By means of the graphic interface of Abaqus/CAE finite element program, tridimensional models were made, corresponding to a miniscrew, section of maxillary bone, premolar and molar teeth taken from computerized tomograph cuts in Dicom format. All parts were considered solid and homogeneous sections with isotropic and linear-elastic properties. A 2N force was applied on the miniscrew, decomposed in two components at 45º on the plane normal to its longitudinal axis. For numerical simulation, the miniscrew was placed in nine different mesio-distal directions. A placement orthogonal to the vestibular surface (angulation of 0º), others of -10º, -20º, -30º and -40º in favor of the force applied and 10º, 20º, 30º and 40º contrary to the force applied. By means of von Mises criteria and analysis of main stressed it was observed that there was a variation in stresses as a result of direction, and the most attenuating situation for the cortical bone was that of 0º, whereas at -20º and 20º higher stress concentrations occurred. It was concluded that there was influence of the mesio-distal direction of the miniscrew in relation to load, of stress distribution in bone tissue.

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 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] Este trabalho estuda um modelo acoplado de deformação e fluxo, baseado na teoria de Biot, para solos saturados e não saturados, porém levando-se em consideração a não linearidade física do solo. Na condição não saturada o fluxo do fluido pode ser considerado bifásico ou monofásico. No caso de fluxo bifásico o modelo calcula simultaneamente a pressão de água e do ar, possibilitando assim, a avaliação dos valores da sucção. O grau de saturação é então determinado a partir dos valores obtidos de sucção. Uma vez conhecido o grau de saturação é possível avaliar a permeabilidade relativa do solo para o fluxo de água e ar, tornando o problema não linear em relação à permeabilidade. No problema de fluxo monofásico a pressão de ar é mantida constante, assumindo-se o valor da pressão atmosférica. O comportamento não linear da relação tensão-deformação em solos não saturados é avaliado pelo modelo básico de Barcelona (BBM), baseado na teoria da plasticidade e conceitos de comportamento de solos no estado crítico. Este modelo pode ser compreendido como uma generalização do modelo Cam Clay freqüentemente usado para representação do comportamento mecânico de solos saturados. A influência da sucção na resistência ao cisalhamento de solos não saturados é verificada por meio de algumas aplicações numéricas. Neste trabalho, a não linearidade física de solos secos e saturados também é simulada através do modelo elastoplástico de Lade & Kim. / [en] This dissertation studies a coupled flow and deformation model for saturated and non-saturated soils, based on Biot`s theory but taking into account the nonlinear material behavior of soil. In the non-saturated condition the fluid flow can be considered either as a two-phase or an one-phase problem. In the two-phase case, the model calculates the pressure of air and water simultaneously, thus permitting the evaluation of the pore suction values and the degree of saturation. Once known, the degree of saturation is used for the determination of the relative permeability of the soil considering the flow of water and air, solving a problem that is nonlinear with respect to permeability. In the one-phase flow, the pressure of air is kept constant and equal to the atmospheric pressure. The stress-strain nonlinear behavior of non-saturated soils is modeled through the basic Barcelona model (BBM model) based on plasticity theory and concepts from the critical state behavior of soils. This model may be understood as an extension of the modified Cam Clay model frequently used for representation of the mechanical behavior of saturated soils. The suction influence on the shear resistance of non-saturated soils is examined through some numerical applications. In this work, the material nonlinearity of dry or saturated soils is also simulated through the Lade & Kim elastoplastic model.

[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.

[pt] O processo de fraturamento hidráulico tem sido amplamente usado para aumentar o volume de petróleo e gás extraído na indústria petroleira. Durante a injeção de fluido, uma região determinada do reservatório é fraturada com a finalidade de aumentar a permeabilidade do meio poroso, de tal forma que o fluxo do fluido desde o reservatório para o poço seja favorecido notoriamente. Porém, este processo pode ocasionar danos ambientais tais como contaminação de aguas subterrâneas, vazamentos de gás, fraturas indesejadas nas camadas capeadoras pela injeção de agua e atividade sísmica fazendo primordial um estudo rigoroso do fraturamento hidráulico com a finalidade de reduzir os riscos potenciais associados a esta operação. Umas das metodologias usadas para projetar o fraturamento hidráulico é a simulação computacional. É possível determinar o volume injetado e a potência da bomba de injeção necessária para obter a geometria de fratura (comprimento, abertura e altura) desejada. A modelagem numérica através de elementos coesivos acoplados do processo de fraturamento hidráulico pode ser efetuada considerando o processo transiente ou permanente, tendo geometrias da fratura e curvas de injeção diferentes. Neste trabalho foi simulado numericamente o modelo KGD nos regimes de fluxo transiente e permanente para dois casos de estudo, (1) injeção numa única camada e (2) injeção em três camadas com contraste de tensões e poropressões entre elas. O estudo numérico foi desenvolvido usando o método dos elementos finitos com modelo de zona coesiva no software Abaqus o qual foi comparado com as soluções analíticas do KGD no regime dominado pela rigidez (Vértice-K) para uma camada e de Simonson e Fung para três camadas. / [en] The hydraulic fracturing process has been widely used to improve oil and gas recovery in the petroleum industry. During the fluid injection, the desired section of rock formation is fractured in order to increase the permeability of the medium that can facilitate the flow of oil to a producing well. However, this process can lead to potential environmental risks such as seismic activity, unwanted fractures in cap layers by water injection, water contamination and gas leakage making primordial to develop a rigorous study in order to reduce this environmental risks associated to hydraulic fracturing. One of the studies developed to design the hydraulic fracturing is computational simulation to determine the fluid volume and hydraulic horsepower required in order to produce the wanted fracture geometry (length, opening and height). The numerical modelling of fracturing process by using fully coupled cohesive element hydraulic can be carried out considering either a steady state or a transient analysis, which modify the fracture geometry and injection pressure. In this work, the KGD model is simulated in transient and steady analysis for two cases: (1) injection in a single layer formation and (2) injection in tri-layered formation with stress and porepressure contrast between them. The numerical simulation of a hydraulic fracturing is carried out using the finite element method with the zone cohesive model in Abaqus whose results are compared with analytical solutions of toughness-dominated propagation regime for the one layer formation model and Simonson and Fung analytical solutions for tri-layered formations model.