Este estudo estende a abordagem micromecânica baseada em células computacionais e a abordagem baseada em deformação utilizando o critério de CTOA para descrever extensão dúctil de trincas longitudinais em dutos de aços de elevada tenacidade sob alta pressão. Ensaios laboratoriais realizados em aços API 5L X60 e X70 utilizando corpos-de-prova padronizados de fratura com trinca profunda, produzem os resultados necessários para medir as curvas de resistência à propagação de trincas e calibrar os parâmetros dos modelos de Gurson e CTOA para estes materiais. O foco central deste trabalho é a aplicação da metodologia de células computacionais e o critério de CTOA para prever as pressões de colapso medidas experimentalmente em dutos de paredes finas contendo trincas longitudinais. O programa experimental inclui espécimes pré-trincados com 508 mm (20") e 219 mm (8 5/8") com diferentes profundidades de trincas. Foram conduzidas análises de elementos finitos, sob estado plano de deformações e no contexto 3D pleno utilizando modelos detalhados dos dutos para descrever a extensão de trinca com o incremento de pressão. As simulações numéricas demonstram a eficiência e limitações de ambas as abordagens para descrever a resposta de crescimento de dúctil e prever a pressão de colapso para os dutos testados. Enquanto o critério de CTOA ainda parece ter limitada aplicabilidade para prever o comportamento dúctil dos dutos, a previsão da resposta dúctil do dutos pré-trincados utilizando o modelo de células computacionais mostrou boa concordância com as pressões de colapso medidas experimentalmente. / This study extends a micromechanics approach based upon the computational cell methodology incorporating the Gurson model and a deformation-based approach using the CTOA criterion to describe ductile crack extension of longitudinal crack-like defects in high pressure pipeline steels. Laboratory testing of API 5L X60 and X70 steels at room temperature using standard, deeply cracked fracture specimens provides the data needed to measure the crack growth resistance curve and to calibrate the Gurson and the CTOA parameters for these materials. A central focus of the paper is the application of the cell methodology and the CTOA criterion to predict experimentally measured burst pressures for thin-walled gas pipeline containing longitudinal cracks. The experimental program includes precracked pipe specimens with 508 mm (20 inches) O.D. and 219 mm (8 5/8 inches) O.D. with varying crack depth to thickness ratios (at). Plane-strain and full 3D computations are conducted on detailed finite element models for the pipe specimens to describe crack extension with increased pressure. The numerical simulations demonstrate the effectiveness and limitations of both approaches to describe crack growth response and to predict the burst pressure for the tested pipes. While the CTOA criterion still appears to have limited applicability to predict ductile cracking behavior for the pipe specimens, the cell model predictions of the ductile response for the precracked pipes show good agreement with experimentally measured burst pressures.

O presente estudo mostra que o conceito de materiais com gradação funcional (MGF) pode ser aplicado no desenvolvimento de componentes de fibrocimento mais econômicos e com desempenho mecânico melhorado. Primeiramente, este trabalho estabelece o conceito de fibrocimento com gradação funcional. Depois, modelam-se, em elementos finitos, as tensões em telhas onduladas sob flexão estática e cargas de vento e mostra-se que, para esses tipos de carregamento, regiões específicas das telhas são submetidas a tensões de baixa intensidade justificando a aplicação do conceito de MGF para a otimização da distribuição da resistência em tais componentes. A parte experimental deste trabalho é constituída de quatro estudos: a) produção de fibrocimentos com o teor de fibras variando ao longo da espessura; b) desenvolvimento de um método de escolha de formulações para fibrocimentos com gradação funcional; c) desenvolvimento de compósitos cimentícios para extrusão em altas velocidades, destinados a gradação de fibrocimentos; e d) aplicação pré-industrial para avaliar a aplicabilidade de técnica de modificação local das propriedades de telhas com gradação funcional. Mostra-se, experimentalmente, que fibrocimentos contendo 1,0% de fibras de PVA, em massa, estrategicamente distribuídas ao longo da espessura, apresentam módulo de ruptura (MOR) similar ao dos fibrocimentos homogêneos com 1,8% de fibras, comprovando a eficácia deste tipo de gradação na redução do custo de placas de fibrocimento. O segundo trabalho experimental avalia formulações de fibrocimento constituídas de seis matérias-primas empregando-se a estratégia screening, uma metodologia promissora para gerar regras de mistura, otimizar custos e desempenho e facilitar a escolha de formulações locais para fibrocimentos com gradação funcional. Compósitos cimentícios com fibras de PVA ou vidro álcali resistente (AR), destinadas à extrusão em altas velocidades, foram desenvolvidos para aplicação pré-industrial. Aqueles contendo 4,0% de fibra de vidro, em volume, apresentaram resistência à tração média igual a 12,0 MPa, enquanto compósitos com 3,0% de fibras de PVA apresentaram resistência igual a 7,5 MPa. A aplicação destas misturas entre as camadas de telhas onduladas de fibrocimento, durante o experimento pré-industrial, resultou em melhoras substanciais do desempenho local. / The present study applies the functionally graded materials (FGM) concept on the development of fiber cements as an alternative to improve the mechanical performance and reduce cost of production of asbestos-free corrugated sheets. Primarily, this work establishes the concept of functionally graded fiber cement. Then, finite element modeling (FEM) is applied in the evaluation of the stresses distributions in corrugated sheets under bending and static wind loads. The simulations show that for these load cases, some regions of corrugated sheets are submitted to low stress intensities what justify the application of FGM concept for optimization of the strength distribution in such components. The experimental part of this work consists of four different studies: a) production of functionally graded fiber cements with gradation through thickness; b) development of a method for choosing formulations in functionally graded fiber cements; c) development of fiber reinforced cementitious composites for high speed extrusion and; d) pre-industrial application to evaluate the suitability of functionally graded fiber cements production. The first experimental study shows that functionally graded cementitious composites with 1.0% of PVA fibers, in mass, strategically distributed through thickness, presented similar modulus of ruptures (MORs) to that homogeneous ones with 1.8% of PVA fibers what proves that gradation of fiber content through thickness is a good approach to reduce cost of fiber cement pads without affecting the mechanical performance of the composite. The second experimental work examines sixcomponent fiber cement formulations employing the screening strategy which is a promising methodology to generate coherent mixture rules, optimize cost and performance and even facilitate the choosing of formulations for functionally graded fiber cements. PVA fibers and glass fiber reinforced cementitious composites for high speed extrusion were developed forpre-industrial application. Composites with 4.0% of alkali-resistant (AR) glass fiber, in volume, presented average tensile strength of 12.0 MPa while composites with 3.0% of PVA fibers presented 7.5 MPa. The application of these mixtures between layers of fiber cement corrugated sheets results in improvement of the performance of the product.

Le présent travail concerne l’étude de l’usinage de l’alliage de titane Ti-6Al-4V, matériau très apprécié par les industries aéronautique, biomédicale et de l’énergie. Les qualités des alliages de titane sont nombreuses : haute résistance aux températures élevées et à la corrosion, haute résistance mécanique, biocompatibilité, etc. Cependant, certaines propriétés physiques de ces matériaux, comme leur faible conductivité thermique, conduisent à des difficultés lors de leur mise en forme par usinage. Des études ont été et sont toujours conduites afin de comprendre le comportement de ces matériaux lors de leur mise en forme. Peu de travaux portent sur la prise en compte de la microstructure dans le comportement des alliages de titane lors du procédé d’usinage. Cette dimension constitue l’une des originalités de ce travail de thèse. Les phénomènes microstructuraux sont caractérisés à travers une étude expérimentale en coupe orthogonale de l’alliage Ti-6Al-4V. Les efforts, la température, la morphologie des copeaux et la microstructure sont analysés et interprétés. Une étude numérique du processus de coupe par simulation éléments finis est employée pour comprendre le rôle de l’endommagement et de la recristallisation. A partir des conclusions de ces différentes études, la construction d’un nouveau modèle de comportement est proposée. Ce modèle est appliqué à une modélisation élément fini pour différentes conditions de coupe afin d’étudier l’influence des paramètres d’usinage. Le modèle est validé par comparaison aux résultats expérimentaux. Il est ensuite exploité afin de proposer une analyse du processus de la coupe et notamment de la formation du copeau. / The present work concerns the study of the machining of titanium alloy Ti-6Al-4V. This material is much appreciated by the aerospace, biomedical and energy industries for its advantageous properties: high resistance to high temperatures and corrosion, high mechanical strength, biocompatibility, etc. However, certain physical properties of these materials, such as their low thermal conductivity, lead to difficulties during the machining process. Studies have been and are still conducted to understand the behavior of these materials during their shaping. Few studies consider the influence of microstructure on the behavior of titanium alloys during the machining process. This is one of the originalities of the present work. The microstructural phenomena are characterized through an experimental study of orthogonal cutting of the Ti-6Al-4V alloy. Machining forces, temperature, chip morphology and microstructure are analyzed and discussed. A numerical study of the finite element simulation process is used to understand the role of damage and recrystallization. From the conclusions of these different studies, the construction of a new model of behavior is proposed. This model is applied to finite element modeling for different cutting conditions to study the influence of machining parameters. The model is validated by comparison with the experimental results. It is then used to propose an analysis of the microstructural phenomena during the cutting process and the formation of the chip.

La caractérisation du comportement des matériaux bitumineux est primordiale pour pouvoir dimensionner les structures de chaussées, mais aussi pour prédire de manière plus précise leur durée de vie. En effet, ces matériaux sont soumis à des phénomènes complexes, mécaniques, thermiques, physiques et chimiques qui apparaissent souvent de manière couplée. Devant la complexité des problèmes observés et avec l’apparition des nouveaux matériaux et nouvelles techniques de fabrication, des études avancées de laboratoire et des modélisations rhéologiques sont nécessaires. Les méthodes existantes permettent de caractériser les matériaux en laboratoire grâce à l’utilisation de presses hydrauliques très couteuses. De plus, la plupart de ces essais ne sont pas adaptés à une utilisation directe sur chaussées. Dans cette thèse, un essai dynamique qui repose sur le phénomène de propagation d’onde a été développé pour caractériser les enrobés bitumineux dans le domaine viscoélastique linéaire. Il s’agit d’un essai non destructif, bien moins couteux que les essais existants et qui pourrait à terme s’adapter à une utilisation in-situ. Des méthodes complexes d’analyse inverse de propagation d’ondes en milieu viscoélastique ont été mises au point afin d’obtenir les propriétés thermomécaniques des enrobés bitumineux à partir des signaux mesurés expérimentalement. Les propriétés viscoélastiques linéaires de différents enrobés bitumineux déterminées à partir de l’essai dynamique développé ont été comparées aux résultats obtenus avec des essais cycliques plus conventionnels de traction-compression. / Characterizing the behaviour of bituminous materials is essential for the design of pavement structures and to predict their service life more accurately. Indeed, these materials are subjected to complex phenomenon, mechanical, thermal, physical and chemical that are often coupled. Due to the complexity of the phenomenon observed and with the emergence of new materials and new production techniques, advances laboratories studies and rheological modelling are necessary. With the existing methods, it is possible to characterize the materials in laboratories using expensive hydraulic presses. Moreover, most of these tests cannot be used directly on pavement structures. In this thesis, a dynamic test based on wave propagation was developed to characterize bituminous mixtures in the linear viscoelastic domain. It is a nondestructive test, far less expensive than the existing tests and that could be adapted for in situ measurements. Complex methods of inverse analysis of wave propagation in a viscoelastic medium were developed to obtain the thermomechanical properties of bituminous mixtures from the signals measured experimentally. The linear viscoelastic properties of different mixtures obtained from the developed dynamic test were compared to the results of more conventional cyclic tension-compression tests.

Cette thèse traite de la modélisation du remodelage osseux. Nous présentons tout d'abord un modèle général continu tenant compte de la réponse cellulaire à un stimulus mécanique. Ce modèle est appliqué à des géométries 2D et 3D macroscopiques afin de se rapprocher des problématiques réelles, ainsi que sur des géométries mésoscopiques d'os trabécullaires en 2D. Cependant la complexité du remodelage osseux ne permet pas d'avoir une approche unique de modélisation. Ainsi, dans un second temps, le cas particulier du remodelage osseux orthodontique est étudié. Un nouveau modèle spécifique est développé tenant compte de l'influence du ligament parodontal sur le remodelage osseux, et intégrant l'influence du taux d'oxygène qui contrôle les évolutions de densités cellulaires. Des données expérimentales in vitro sont extraites de la littérature et servent de données d'entrées du modèle développé afin d'obtenir l'évolution de la densité osseuse alentours d'une racine dentaire cylindrique en 3D. / This work deals with modelization of bane remodeling. We present first a madel thal accounts for the cellular res panse to a mechanical stimulus in a general case at a continuous scale. This madel is applied to 2D and 3D geometries at macroscopic scale to mimic real cases, as weil as 2D trabecular-type geometries at mesoscopic scale. However, the complexity of bane remodeling does not allow a unique approach. Th us, the thesis work is focused on the particular case of orthodontie bane re mode ling. A new specifie madel is developed accounting for the influence of the periodontal ligament on orthodontie bane remodeling by integrating the oxygen concentration effect controling the evolutions of cellular densities. The cellular experimental data in vitro are extracted from the literature, and serve as input data of the developed madel in arder to ablain the evolution of bane density around the root of a 3D cylindrical tooth.

Le système ostéo-articulaire a un rôle de soutien essentiel pour l’ensemble du corps humain. Il supporte les forces de gravité et les efforts produits par les activités quotidiennes. Ainsi l’os optimise et adapte sa masse et sa géométrie à travers le processus de remodelage osseux. Cette régulation spatio-temporelle peut subir des déséquilibres métaboliques comme l’ostéoporose qui conduisent à la survenue de fracture de l’extrémité supérieure du fémur lors d’une chute latérale ou de traumatismes divers. En effet l’ostéoporose se traduit par une diminution au niveau de la densité osseuse et des détériorations au niveau de la microarchitecture du tissu osseux augmentant ainsi le risque de fractures.Dans le cadre de cette thèse une modélisation mécanique du comportement de l’os ainsi qu’une modélisation biologique des activités cellulaires ont été proposées dans un premier temps. Dans un second temps, la pharmacocinétique de quatre médicaments a été modélisé soient l’Alendronate, le Denosumab,le Romosozumab et l’Odanacatib. Finalement, le couplage de ces modèles a permis d’avoir un modèle mécano-biologique couplé aux effets des médicaments contre l’ostéoporose qui permet de prédire l’évolution de la densité minérale osseuse et celle de l’endommagement par fatigue permettant ainsi d’analyser l’évolution de la qualité osseuse. Ce modèle a été implémenté au code de calcul par éléments finis ABAQUS/standard à travers sa routine utilisateur UMAT. Le modèle a été appliqué pour simuler différents scénarii de remodelage sur des fémurs humains (2D et 3D). Différents facteurs ont été analysés tels que l’amplitude des activités physiques, les doses de médicaments injectées, la durée des traitements, etc. les résultats obtenus sont cohérents (qualitativement) avec les études cliniques existantes. En conclusion, le modèle mécano-biologique couplé aux effets des médicaments proposés contribue à l’analyse fine du comportement de l’os et l’application des algorithmes a permis d’effectuer des essais virtuels permettant d’analyser les effets combinés de nombreux facteurs pluridisciplinaires caractérisant la qualité osseuse. / The osteo-articular system plays the role of crucial support for the whole human body. It supports the gravity forces as well as the efforts generated by daily activities. Thus the bone optimizes and adapts its mass and its geometry through the process of bone remodeling. This spatio temporal regulation can undergo metabolic imbalances such as osteoporosis which lead to the occurrence of the upper end of the femur fracture during a fall side or various traumas. In fact, osteoporosis is reflected into a decrease in bone density and level of damage at the level of the microarchitecture of bone tissue increasing therefore the risk of fractures. In the framework of this thesis, a mechanical modelling of the behaviour of the os as well as a biological modelling of cellular activities were firstly proposed. Secondly, the pharmacokinetics of four drugs were modeled which are Alendronate, Denosumab, the Romosozumaband the Odanacatib. Finally, the coupling of these models allowed us to have a mechanic-biological model coupled with the effects of drugs against osteoporosis that can predict the evolution of bone mineral density and of damage by fatigue allowing to analyze the evolution of the bone quality. This model has been implemented by finite elements ABAQUS/standard through its user routine UMAT.The model has been applied to simulate different scenarios of remodeling on human femurs (2D and3D). Different factors were analysed such as the range of physical activities, the doses of injected drugs,the duration of treatment, etc. Obtained results are consistent (qualitatively) with existing clinical studies.As a conclusion, the mechanic-biological model coupled with the effects of the proposed drugs contributes to the fine analysis of the bone behaviour and the application of algorithms allowed to conduct virtual tests in order to analyze the combined effects of many multidisciplinary factors characterizing the bone quality.

Les dépôts effectués sur substrats de silicium sont très courants notamment dans le domaine de la microélectronique. Les propriétés physiques recherchées pour ce type de structures dépendent fortement de celles de la couche. Il apparaît donc essentiel de connaître les paramètres élastiques ainsi que l’épaisseur des films considérés. De plus, la détection de certains défauts concernant la couche est souvent recherchée. L’objectif de ce travail a été de contribuer à la caractérisation de structures du type couche sur substrat. Pour cela, les ultrasons-lasers présentent de nombreux avantages puisqu’ils autorisent entre autres leur contrôle non destructif sans contact. Les ondes acoustiques de surface dans une gamme de fréquence s’étendant jusqu’à 45 MHz ont été utilisées. Nous avons développé différents modèles analytiques et les résultats expérimentaux ont aussi été comparés à certaines simulations par éléments finis. Plus particulièrement, nous avons montré qu’il était possible d’obtenir l’ensemble des paramètres élastiques du substrat et de la couche ainsi que l’épaisseur de cette dernière. Par ailleurs, nous nous sommes aussi intéressés à la détection de certains défauts en régime impulsionnel mais aussi quasi-monochromatique. Des résultats originaux concernant l’effet d’une absence de couche de forme déterminée sur le premier mode de Rayleigh ou bien encore de problèmes d’adhésion ont été présentés. Sur ce dernier point, une méthode innovante permettant de distinguer un fort niveau d’adhésion d’un faible a aussi été introduite. / Thin films deposited on silicon substrates are very common especially in microelectronic applications. The physical properties expected for these types of structures depend on the properties of the layer. Therefore, it is essential to know elastic parameters and thickness of the films considered. Furthermore, some layer defects detection is often required. The aim of this work was to contribute to the characterization of structures composed of a single layer deposited on a substrate. For this, laser ultrasonics offers many advantages since it is nondestructive and non-contact method. Surface acoustic waves excited in a frequency range up to 45 MHz are used. Analytical models have been developed and experimental results have also been compared with some finite element simulations. More specifically, we have shown that it was possible to obtain the thickness of the layer and all elastic parameters of both substrate and layer. In addition, we have been interested in the detection of various defects with a broadband and quasi-monochromatic excitation. Original results concerning the effect of a lack of layer with a specific geometry on the first Rayleigh mode or even adhesion problems were presented. Concerning this last point, an innovative method to tell difference between high and low adhesion was introduced.

Un nouveau capteur acoustique sphérique autonome est proposé. Il est destiné à être immergé dans un milieu liquide ou pâteux pour mesurer certaines propriétés physiques du milieu et récupérer l'énergie vibratoire ambiante pour assurer son autonomie. Le capteur est composé de deux coquilles hémisphériques en plexiglas et d'une bague piézoélectrique en PZ26 fixée entre les deux coquilles. Cette structure peut être utilisée aussi bien en excitateur que capteur. Un modèle de simulation de la récupération d'énergie vibratoire a été développé en considérant seulement deux modes de vibration: mode épaisseur et mode radial. Pour chaque mode, le comportement de l’anneau est décrit par un circuit électromécanique équivalent reliant les ports mécaniques (forces et vitesses) au port électrique (tension et courant). Ce choix est guidé par la possibilité de combiner la partie électromécanique avec l'électronique qui traite l'énergie directement dans un simulateur basé sur Spice. Pour valider cette approche, une simulation par éléments finis a été réalisée et comparée aux résultats produits par le circuit électromécanique. Les fréquences de résonance ont également été vérifiées expérimentalement avec un analyseur d'impédance. Toutes ces vérifications donnent des résultats en très bon accord avec le modèle électromécanique proposé en termes de fréquences de résonance, de tension et de puissance collectées. Enfin, plusieurs validations expérimentales sont présentées avec un prototype de capteur sphérique. Ces validations montrent l’adéquation des prédictions avec les résultats expérimentaux. Finalement, un test du circuit de récupération est effectué en situation réelle. / A new spherical autonomous acoustic sensor is proposed. It is intended to be immersed in a liquid or pasty medium to measure some physical properties of the medium and should harvest ambient energy to ensure its autonomy. The sensor is composed of two Plexiglas half-spherical shells and a PZ26 piezoelectric ring clamped between the two shells. This structure can be used as well as in exciter or sensor. A simulation model of vibrational energy harvesting has been developed considering only two modes of vibration: thickness and radial modes. For each mode, the ring behavior is described by an equivalent electromechanical circuit connecting the mechanical ports (forces and velocities) to the electrical port (voltage and current). This choice is guided by the possibility to combine the electromechanical part with the electronics that process the energy directly in a Spice based simulator. To validate this approach, a finite elements simulation was realized and compared to the electromechanical circuit results. Resonance frequencies were also verified experimentally with an impedance analyzer. All these verifications give results in very good agreement with the proposed electromechanical model, as well as in terms of resonant frequencies, harvested voltage and power. Finally several experimental investigations are presented with a prototype of spherical sensor. These validations show the adequacy of the predictions with the experimental results. Finally, a test of the harvesting circuit is done in real situation.

Les structures mécaniques utilisées de nos jours ne cessent d’évoluer en utilisant des matériaux composites ou à gradient fonctionnel afin de répondre aux enjeux de résistance accrue, allégement de la structure et amélioration des performances. Ceux-ci nécessitent un contrôle adéquat de leur état de santé afin de s’assurer de l’intégrité de la structure. L’utilisation des ondes guidées ultrasonores fournit un moyen efficace et rapide d’inspection sur de longues distances. Néanmoins, ces ondes présentent certaines caractéristiques complexes qui rendent la tâche très difficile. L’utilisation d’outils d’analyse tels que les modèles numériques constitue un grand atout pour ce type d’application. Dans ce contexte, l’objectif de cette de thèse est le développement d’un outil de modélisation performant, permettant d’étudier la propagation des ondes guidées ultrasonores avec une grande précision et une faible consommation de ressources et de temps de calculs. De ce fait, l’intérêt est porté sur des méthodes numériques d’ordres élevés dont les propriétés de convergence sont beaucoup améliorées que les méthodes classiques. En particulier, la méthode semi-analytique éléments finis de Galerkin discontinue pour la détermination des courbes de dispersion des ondes guidées est développée. La méthode est applicable aux structures planes et cylindriques fabriquées de matériaux isotropes, anisotropes et hétérogènes (à gradient fonctionnel de propriétés). Une étude comparative sur l’analyse des performances de ces méthodes est effectuée. Celle-ci a démontré la capacité de la méthode à modéliser la propagation des ondes guidées ultrasonores dans des guides d’ondes à section arbitraire avec des performances prometteuses par rapport à la méthode des éléments finis classique. / The mechanical structures used today are constantly evolving using composite or functionally gradient materials to meet the challenges of increased strength, lightening the structure and improving performance. These require adequate control of their state of health to ensure the integrity of the structure. The use of Ultrasound Guided Waves (UGW) provides an efficient and fast way of inspection over long distances. Nevertheless, these waves have some complex features that make the task very difficult. The use of analysis tools such as numerical models is a great asset for this type of application. In this context, the objective of this thesis is the development of a powerful modeling tool, allowing to study the propagation of UGWs with a great precision, less computational time and consumption of resources. Accordingly, we are interested in higher order numerical methods whose convergence properties are much improved than the classical methods. In particular, a semi-analytical discontinuous Galerkin finite element method (SADG-FE) is developped for the determination of the dispersion properties of guided waves in arbitrary cross-section waveguides. The method is applicable to plates and cylindrical structures made of isotropic, anisotropic heterogeneous (functionally graded) materials. The performance analysis of these methods and their comparisons are performed with respect to the models based on the classical finite element method. The results demonstrated the ability of the proposed method to model the propagation of ultrasounic guided waves in arbitrary section waveguides with promising performance over the conventional finite element method.

Orientador: Osvaldo Luís Manzoli / Abstract: This work presents a numerical method based on Discrete Fracture Model (DFM) and the Finite Element Method (FEM), where the fractures are approximated by a reduced model. The flow along and across the fracture is described by a simplified set of equations considering both conductive fractures and barriers. The coupled hydromechanical model is composed of a linear poroelastic Biot medium and a nonlinear model based on damage mechanics for the fractures, which captures the nonlinear normal deformation and shear dilation according to the Barton-Bandis model. Both flow and geomechanical models are approximated using the finite element model. Fractures are explicitly represented by three-node standard finite elements with high aspect ratio (i.e. ratio between the largest and the smallest element dimensions) and appropriate constitutive laws. These interface high aspect ratio elements represent a regularization method which continuously approximate the discontinuous pressure and displacement fields on a narrow material band around the fracture. The complete mathematical formulation is presented together with the algorithm suggested for its numerical implementation. The efficiency of the proposed method is demonstrated through numerical examples, as well as the effects of fractures in the hydraulic properties of porous rocks and its dependency of the stress state. / Resumo: Este trabalho apresenta um método numérico baseado no Modelo de Fratura Discreta (MFD) e no Método dos Elementos Finitos (MEF), onde as fraturas são aproximadas por um modelo reduzido. O fluxo ao longo e através da fratura é descrito por um conjunto simplificado de equações, considerando tanto fraturas condutoras quanto barreiras. O modelo hidromecânico acoplado é composto por um meio poroelástico linear e um modelo não linear para fraturas, baseado na mecânica do dano e que captura a deformação normal não linear e a dilatância ao cisalhamento de acordo com o modelo de Barton-Bandis. Os modelos de fluxo e geomecânico são aproximados usando o método dos elementos finitos. As fraturas são explicitamente representadas por elementos finitos triangulares de três nós com elevada razão de aspecto (isto é, a razão entre a maior e a menor dimensão do elemento) e leis constitutivas apropriadas. Esses elementos de elevada razão de aspecto representam um método de regularização que aproxima de forma contínua os campos de pressão e deslocamento descontínuos em uma estreita faixa material ao redor da fratura. A formulação matemática completa é apresentada juntamente com o algoritmo sugerido para sua implementação numérica. A eficiência do método proposto é demonstrada através de exemplos numéricos, bem como os efeitos de fraturas nas propriedades hidráulicas de rochas porosas e sua dependência do estado de tensão. / Mestre