Dans le cadre des systèmes rate-independent, un modèle de plasticité avec endommagement, visant à la description des processus de rupture ductile, est proposé et étudié par une formulation variationnelle. Une fissure cohésive, ou ductile, se produit lorsque le champ de déplacement subit une discontinuité, tout en étant encore associé à une contrainte de traction non nulle. Prévoir et décrire efficacement les phénomènes de rupture ductile est une tâche cruciale pour de nombreux matériaux d'ingénierie (métaux, polymères, ...), comme en témoigne le grand intérêt de la communauté scientifique sur le sujet. Modèles d'endommagement à gradient ont été fructueusement utilisé pour la description des ruptures fragiles: dans ce cas, une fois que le niveau d'endommagement atteint sa valeur maximale, une fissure est créée lorsque la traction entre les deux faces opposées tombe immédiatement à zéro. D'autre part, le modèle de plasticité parfaite pourrait décrire la formation de la glisse plastique au niveau constant de contrainte. Par conséquent, afin de décrire les effets typiques d'une rupture cohésive, l'idée principale consiste à coupler, par une approche variationnelle, le modèle de plasticité parfaite et un modèle d'endommagement à gradient. L'utilisation d'une approche variationnelle se traduit par une formulation faible et sans dérivées, fournit des moyens efficaces pour traiter les notions de bifurcation et de stabilité, est intrinsèquement discret et indique une manière naturelle et rationnelle pour définir des algorithmes numériques efficaces. L'incorporation des effets d'endommagement dans un modèle de plasticité n'est pas une idée nouvelle. Néanmoins, le modèle proposé présente de nombreux aspects originaux comme le couplage entre la plasticité et l'endommagement et la façon avec laquelle l'évolutions des variables se trouvent. L'approche variationnelle s'appuie simplement sur ​​trois concepts: une condition d'irréversibilité, une condition de stabilité globale, locale ou différentielle et le bilan énergétique. Le modèle résultant possède une grande flexibilité dans les réponses possibles couplées, en fonction des paramètres constitutifs. Ces diverses réponses sont d'abord examinées avec un test d'une barre unidimensionnelle en traction quasi-statique en assumant une évolution homogène qui met en évidence les principales caractéristiques du modèle. La discussion sur la stabilité des solutions homogènes conduit à l'existence d'une longueur de la barre critique qui à son tour dépend de la longueur interne caractéristique du matériel. En considérant des barres plus longues par rapport à cette valeur critique, on démontre que la réponse homogène devient instable. Par conséquence une localisation doit apparaitre dans la barre. Une construction de localisation est ensuite proposée, qui prend explicitement en compte la condition d'irréversibilité sur le champ d'endommagement. Ceci permet d'étudier l'évolution non homogène et la réponse globale. Il s'avère que, en général, une fissure cohésive apparaît au centre de la zone d'endommagement avant la rupture. A ce stade, la déformation plastique se localise comme une mesure de Dirac qui devient responsable de cette fissure cohésive. On obtient la loi cohésive associée en termes de paramètres du modèle et retrouve la loi de fracture cohésive postulée par Barenblatt. Enfin, un schéma de résolution numérique est proposé, qui est basé sur un algorithme de minimisation alternée, et mis en œuvre par une librairie d'éléments finis uniquement pour le test de barre en traction. Même si l'espace d'éléments finis adoptés ne peut pas incorporer les discontinuités, les résultats numériques s'accordent parfaitement avec les solutions analytiques. Néanmoins, les développements futurs visent à étendre les simulations dans un cadre à deux / trois dimensions et de tester une méthode d'éléments finis généralisée.

Cette étude fait partie d'un programme de recherche concernant la réparation de structures composites par collage de patchs externes. Les objectifs principaux de ce programme sont d'une part l'identification de l'ensemble des facteurs susceptibles d'influencer les performances à long terme de ce type de réparation, et d'autre part de déterminer dans quelle mesure l'utilisation de tels assemblages peuvent s'avérer une solution optimale. La conception d'un tel système passe obligatoirement par le développement d'un outil de simulation et de prédiction robuste du fait des divers mécanismes d'endommagement pouvant intervenir de fa?on très complexe et de la rupture finale du système résultant d'une propagation des zones endommagées. Cette étude compose d'une et d'autre l'aspect de la modélisation numérique, et l'aspect expérimental. Le dialogue entre les résultats numériques et expérimentaux permet, d'une part de comprendre les mécanismes d'endommagement et l'évolution de ce dernier dans le système réparé, d'autre part de valider le modèle numérique. En particulier, nos efforts ont été concentrés, en utilisant le logiciel LS-dyna, sur l'application des modèles de zone cohésive (MZC). Le comportement au délaminage d'un composite carbone/époxyde et de l'adhésif sont étudié avec les essais en mode I, mode II et mode mixte. Une étude paramétrique de MZC est effectuée. Le modèle de zone cohésive validé est utilisé pour modéliser le comportement en traction des composites réparés par collage de patches externes

Le collage structural est une technique d’assemblage de plus en plus demandée aujourd’hui dans beaucoup de domaines comme l’automobile, l’aéronautique, l’aérospatial et dans d’autres domaines comme la construction, le sport et les loisirs. Cette technique très avantageuse, permet l’assemblage de matériaux semblables ou différents à l’aide d’un adhésif, la réduction importante du poids et la répartition uniforme des charges sur l’assemblage. Malgré ses avantages, le collage souffre encore de quelques inconvénients liés à l’existence de défauts dans les joints de colle. Parmi eux, il existe des défauts qui sont situés à l’interface colle/substrat comme un « kissing bond » ou un mauvais état de surface, qui restent indétectables ou difficilement détectables utilisant les techniques de contrôle non destructives. Donc, afin de prendre en compte l’existence des défauts d’adhésion dans les assemblages collés lors de la phase de conception, il est nécessaire de fournir un modèle analytique capable de prédire la propagation de fissure. Dans cette thèse, un modèle analytique qui prédit la propagation de fissure et qui évalue la résistance effective d’un assemblage collé contenant des défauts d’adhésion a été développé. Un défaut a généralement une géométrie complexe, et une étude générique est difficilement réalisable ce qui nous amène à considérer des géométries de défauts idéales. Le modèle a été vérifié par des expériences réalisées sur des éprouvettes DCB. Des simulations numériques utilisant la méthode de zone cohésive ont été réalisées également pour décrire plus complètement le processus de décohésion et simuler les essais expérimentaux. La dernière partie de ce travail a été dédiée à l’étude de la fissuration des éprouvettes en alliage de titane. Profitant de la collaboration avec Safran et Alphanov, les substrats ont subi un traitement de surface laser en laissant des zones non traitées. Le but de cette partie était de vérifier le modèle analytique proposé avec des configurations plus complexes. / Structural adhesive bonding has known an increasing use in many fields like aeronautics, aerospace and automotive and other fields like construction and sports. This very advantageous technique allows the assembly of similar or different materials using an adhesive, the significant reduction in weight and a uniform distribution of loads on the assembly. Despite its advantages, the bonding still suffers from some disadvantages related to the existence of defects in the bonded joints. Among them, there are defects that are located at the interface glue / substrate as "kissing bond" or poor surface due to bad surface treatment, which remain undetectable or hardly detectable using non-destructive control techniques. Therefore, in order to take into account the existence of adhesion defects in bonded assemblies during the design phase, it is necessary to provide an analytical model capable of predicting crack propagation and estimate the criticality of a defect. In this thesis, an analytical model that predicts crack propagation and evaluates the effective strength of a bonded assembly containing adhesion defects has been developed. A defect usually has a complex geometry, and a generic study is difficult to achieve, which leads us to consider ideal defect geometries. The model was verified by experiments performed on DCB specimens. Numerical simulations using the cohesive zone method were also performed to more fully describe the decohesion process and to simulate the experimental tests. The last part of this work was devoted to the study of titanium alloy assembly containing patterns. Taking advantage of the collaboration with Safran and Alphanov, the substrates underwent a laser surface treatment leaving untreated areas. The purpose of this part was to check the proposed analytical model with more complex configurations.

Dans le cadre du projet ASTRID, le procédé de fabrication du combustible nucléaire par métallurgie des poudres, pour les RNR, est revisité dans le but d’être simplifié. En particulier, nous cherchons à supprimer l’étape de granulation mécanique de la poudre qui lui confère un bon comportement à l’écoulement. La problématique, est d’arriver à faire s’écouler spontanément et rapidement une poudre, à travers un orifice par lequel elle ne s’écoule pas naturellement. Par ailleurs, la poudre alterne entre des phases d’écoulement et de non-écoulement.Pour garantir l’écoulement de la poudre à travers le moule de presse, nous appliquons des vibrations horizontales. Les vibrations permettent de faire s’écouler la poudre et d’atteindre le débit requis par les cadences de fabrication. Toutefois, elles ont également le défaut, lors des phases de non-écoulement, de compacter la poudre et de retarder les futurs écoulements. Tout l’art du remplissage des moules de presse assisté par vibrations consiste à maîtriser le caractère ambivalent des vibrations.L’évolution singulière de la compacité des poudres d’actinides, lors de la phase de non écoulement, nous a amené à définir un modèle stochastique unidimensionnel de compaction. La confrontation des résultats de simulation, avec les lois empiriques de compaction, a permis d’identifier le sens physique des paramètres d’ajustement des lois empiriques. Nous avons aussi proposé une nouvelle loi d’évolution de la compacité à deux exponentielles étirées. Cette nouvelle loi rend compte, non seulement de la cinétique de compaction des poudres d’actinides mais également de l’ensemble des résultats que nous avons trouvé dans la littérature. / In the framework of the ASTRID project, the nuclear fuel production process by powder metallurgy, for Fast Neutron Reactors, is revisited in order to be simplified. In particular, we seek to remove the mechanical granulation step of the powder which gives a good flow behavior during the filling of press molds. The aim is to reach a spontaneous and quick powder flow through a hole in which the powder does not flow without external energy. Furthermore, the powder alternates between flow phases during the filling of press molds, and non-flow phases during the compaction and ejection of the pellet.We hence apply horizontal vibrations to ensure the flow of the powder through the press mold. The vibrations help the powder to flow and increase the production rates. However, they have the disadvantage to compact the powder and delay the future flows, during the non-flow phases. The art of filling the press mold assisted by vibrations is to master / control the ambivalent nature of the vibrations.The remarkable packing fraction evolution of actinides powders, during the non-flow phases, allows us to define a simple 1D stochastic model to understand the compaction kinetics. The comparison of the stochastic model with the empirical compaction laws found in the literature helps us to identify the physical meaning of fitting parameters proposed by the empirical models. Furthermore, we have also proposed a new compaction law with two-stretched exponentials. This new law not only reflects the compaction kinetics of actinides powders, but also of all the compaction data we found in the literature.

Les composites stratifiés de plis unidirectionnels en carbone/époxy sont fortement utilisés pour alléger les structures aéronautiques tout en conservant de bonnes propriétés structurales. Toutefois, les avantages de ce type de matériau ne sont pas encore pleinement exploités de par le manque de confiance accordée aux modèles de prévision de l’endommagement, dont notamment ceux concernant le délaminage. C’est pourquoi l’objectif de cette thèse était de développer une stratégie de modélisation du délaminage adaptée aux structures composites stratifiées. Cette stratégie s’est appuyée sur le développement d’un modèle de zone cohésive prenant en compte les ingrédients nécessaires à la bonne description de l’amorçage et de la propagation de la fissure : (i) un critère d’amorçage avec un renforcement en compression/cisaillement hors-plan, (ii) une loi de propagation décrivant l’évolution de la ténacité en fonction de la mixité de mode et (iii) la prise en compte du couplage inter/intralaminaire. Pour identifier ce nouveau modèle, une procédure d’identification efficace, s’appuyant sur un essai de traction sur plaque rainurée, a été mise en place. Cette procédure d’identification a permis de démontrer que la ténacité semble indépendante (i) de l’orientation des plis adjacents à l’interface et (ii) de l’empilement étudié. De même, pour décrire l’évolution de la ténacité, une nouvelle loi de propagation adaptée au matériau carbone/époxy a été proposée. Pour finir, la stratégie de modélisation, complétée par une stratégie de calcul, a été appliquée sur différents cas structuraux pour mettre en avant ses apports et ses premières limites. / The carbon/epoxy laminated composites of unidirectional plies are strongly used in order to reduce the weight of aeronautical structures while at the same time proposing good structural properties. However, the advantages of this kind of material are not fully exploited due to the lack of confidence in damage models, like ones concerning delamination. Then, the purpose of this work was the development of a strategy to model delamination in laminated composite structures. This strategy was based on the development of a cohesive zone model taking into account the ingredients necessary to the well description of the onset of delamination and the crack growth: (i) an onset criterion with an out-of-plan compression/shearing reinforcement, (ii) a propagation law describing the evolution of the fracture toughness as a function of mixed mode ratio and (iii) the inter/intralaminar coupling. To identify this new model, an efficient identification procedure, basing on a tensile test on notched specimen, has been proposed. This identification procedure has demonstrated that the fracture toughness seems to be independent of (i) the orientation of plies closed to the interface and (ii) the stacking sequence. Furthermore, to describe the evolution of the fracture toughness, a new propagation law adapted to carbon/epoxy material has been proposed. Finally, the strategy to model delamination, completed by a calculation strategy, has been applied on several structural cases to prove its contributions and its first limitations.

Thèse numérisée par la Direction des bibliothèques de l'Université de Montréal.

Le travail présenté a pour objectif la prédiction numérique de la résistance résiduellede grandes structures vis-à-vis d’évènements accidentels, tels que ceux rencontrés p.ex. dans le cas de la collision de navires ou d’impact d’oiseaux en aéronautique. Cesévènements peuvent dans certain cas conduire à la rupture, qui est ici considéréeductile. La difficulté de cette étude, consiste à reproduire dans une méthodologieunifiée basée sur la méthode des éléments finis les étapes successives menant àla ruine ultime de la structure. Ces étapes sont : l’endommagement ductile, lalocalisation de la déformation et la propagation de la fissure. Un élément essentiel pour la conception d’un modèle de fissuration ductile prédictif est le traitement numérique de la phase transitoire critique de localisation de la déformation induite par l’endommagement dans une bande de matière étroite.A cet effet, trois points de vue différents en termes de champ de déplacement àtravers la bande de localisation sont proposés. Ces trois approches se distinguentpar le type de discontinuité considérée : forte, faible et régularisée (expression nonlinéaire). Un cadre variationnel consistant est élaboré pour chacune des trois approches.Ainsi la cinématique enrichie est incorporée dans la formulation de l’élément fini enutilisant la méthode des éléments finis enrichis (X-FEM). Puis, la performance deces méthodes est évaluée vis-à-vis de leur capacité à modéliser la phase transitoireentre endommagement diffus (mécanique des milieux continus) et propagation defissure (mécanique de la rupture). Ces travaux sont réalisés dans le contexte dematériaux ductiles. D’après les analyses réalisées, la combinaison du modèle de ’discontinuité fortecohésive’ et la X-FEM semble être la plus prometteuse des trois approches étudiéespour allier physique et numérique. Le développement d’un tel modèle est discutéen détail. Enfin, deux critères supplémentaires sont définis : le premier pour lepassage de l’endommagement diffus au modèle de bande cohésive et un deuxièmepour le passage du modèle de bande cohésive à la rupture. / The present work aims at numerically predicting the current residual strengthof large engineering structures made of ductile metals regarding accidental events,e.g. ships collision or bird strike in aviation, which may potentially lead to failure.With this aim in view, the challenge consists in reproducing within a unified finiteelement (FE)-based methodology the successive steps of micro-voiding-induceddamage, strain localization and crack propagation, if any.A key ingredient for a predictive ductile fracture model is the proper numericaltreatment of the critical transition phase of damage-induced strain localizationinside a narrow band. For this purpose, three different viewpoints in terms ofdisplacement field across the localization band are proposed involving a strong,weak and (non-linearly) regularized discontinuity, respectively.A consistent variational framework is elaborated for each of the three methods,whereby the enriched kinematics is embedded into the FE formulation using theeXtended FEM. Then, within a comparative procedure, the performance of thesemethods is assessed regarding their ability of modeling the transition phase betweendiffuse damage (continuum mechanics framework) and crack propagation (fracturemechanics framework), always in the context of ductile materials.According to the aforementioned analyses, the combination of the strong discontinuitycohesive model and the X-FEM appears to be the most promising of thethree studied approaches to bring together physics and numerics. The developmentof such a model is discussed in detail. Finally, two supplementary criteria aredefined: the first one for the passage from diffuse damage to the cohesive bandmodel and the second one for the passage from the cohesive band model to thecrack.

Du point de son comportement mécanique, la maçonnerie apparaît comme un géomatériau quasi-fragile composite, caractérisé par la coexistence de trois échelles d'études (macroscopique, mésocopique et microscopique) et deux comportements intimement couplés, l'endommagement fragile et la plasticité non-associée. Suite à une revue bibliographique étendue des nombreuses modélisations déjà proposées (approche macroscopique phénoménologique, approche micromécanique, analyse multi-échelle), il semble que deux points importants restent encore difficiles à traiter : la prise en compte de la structure de l'appareil et la localisation de la déformation liée au caractère adoucissant et non-associé des comportements. <br />Afin de répondre à ces problèmes, une modélisation micromécanique discrète dans un cadre dynamique non-régulier a été proposée. La maçonnerie est vue comme un milieu divisé composé d'une collection de corps déformables reliés par des lois non-régulières. Le modèle de base pour les joints vifs, composé du contact unilatéral et du frottement de Coulomb est enrichi pour les joints de mortier par des modèles de zone cohésive, fragile ou à endommagement progressif. Le cadre numérique de cette modélisation s'appuie sur la méthode ``Non-Smooth Contact Dynamics'' apte à intégrer la dynamique non-régulière en présence de contraintes unilatérales et de frottement sec. Le comportement des corps est traité par des méthodes aux éléments finis adaptées au caractère discret des structures modélisées. Les lois non-régulières d'interfaces sont résolues quant à elles par des méthodes itératives dédiées aux problèmes de complémentarité de grande taille.<br /> Enfin, une collaboration interdisciplinaire a été conduite avec les architectes et les archéologues du bâti pour mener à bien des études d'édifices monumentaux en utilisant la stéréophotogrammétrie numérique.

Cette thèse est une contribution à l'amélioration de la connaissance des comportements mécaniques des tissus biologiques, en particulier dans le contexte accidentologique. L'objectif de cette étude biomécanique est de faire le lien entre les notions de lésions utilisées en clinique et celles d'endommagement et de rupture utilisées en mécanique. Elle porte sur la modélisation du comportement des ligaments dans l'articulation du genou humain, et s'intéresse plus particulièrement à l'insertion du ligament dans l'os. Les lésions qui touchent cette structure se produisent soit dans le ligament, soit dans une région proche de la zone de transition entre l'os et le ligament. Ce problème est original sur les plans anatomique et mécanique. Sur le plan anatomique, il n'existe pas de description fine (c'est-à-dire à l'échelle microscopique) de l'architecture de cette transition, qui paraît brutale à l'oeil nu et qui est un lieu possible de lésion ligamentaire. Sur le plan mécanique, il s'agit de l'étude d'un tissu biologique, qui comprend de plus la transition entre un tissu dur minéralisé (l'os) et un tissu mou non minéralisé à fibres longues (le ligament). La méthodologie a donc été de comprendre le fonctionnement de l'articulation du genou et les lésions et traumatismes dont elle est le lieu, de connaître la composition et l'organisation des tissus qui constituent cette transition, au moyen d'un étude histologique, de décrire son comportement mécanique dans le fonctionnement du ligament, par le développement d'un protocole expérimental, et enfin de développer un modèle de comportement de l'insertion ligamentaire. L'étude porte sur les ligaments croisé postérieur et latéral externe, ligaments dont les mécanismes lésionnels dans le cas des accidents de la route ont été plus particulièrement étudiés. L'architecture de la zone de transition en microscopies optique et électronique a été décrite et a montré qu'elle est la superposition d'un front de minéralisation et d'un changement de la structure du tissu (présence de cartilage fibreux). Cette transition se fait sur une longueur d'environ 300 micromètres. Les essais expérimentaux de traction ont été réalisés sur la structure insertion ligamentaire - ligament - insertion ligamentaire prélevée sur sujets d'anatomie. La sollicitation a été appliquée soit dans la direction des fibres ligamentaires, soit une direction réaliste du point de vue physiologique. Les protocoles développés permettent de faire des essais cycliques sans endommager les tissus et des essais à rupture sous sollicitations quasi-statiques (1 mm/s et 20 mm/s) et dynamiques (0.5 m/s et 1 m/s). Pour les essais dans la configuration physiologique, les genoux ont été testés en extension complète (station érigée) et à 120 degrés de flexion (position de conduite). Les résultats expérimentaux montrent que, pour les vitesses et amplitudes testées, les ligaments dissipent de l'énergie par frottement interne et que leur comportement est très dépendant de l'angle de flexion du genou. La rupture se produit toujours par arrachement osseux au niveau de l'os cortical en quasi-statique. En dynamique, pour la vitesse testée, il y a toujours décohésion entre fibres dans le ligament, mais la rupture se produit le plus souvent par arrachement de l'os spongieux en profondeur. L'étude du comportement du ligament (influence de l'orientation de l'insertion ligamentaire) et de la micro-structure de l'insertion nous ont menés à choisir un modèle d'interface pour décrire son comportement. Deux modèles de zones cohésives (couplant frottement, adhésion et éventuellement endommagement) ont été développés pour prédire les lésions par arrachement osseux et par décohésion entre les fibres ligamentaires.

La perméabilité des roches est fortement influencée par la présence de fractures car ces dernières constituent un chemin préférentiel pour l’écoulement des fluides. Ainsi la présence de fractures naturelles est un facteur déterminant pour la productivité d’un réservoir. Dans le cas de roches à faible conductivité, des techniques de stimulation telle que la fracturation hydraulique sont utilisées pour en augmenter la perméabilité et rendre le réservoir exploitable d’un point de vue économique. A l’inverse, dans le cas du stockage géologique, la présence de fractures dans la roche représente un danger dans la mesure où elle facilite le transport et la migration des espèces disséminées dans la roche. Pour le stockage de CO2, les fuites par les fractures présentes dans les couvertures du réservoir et la réactivation des failles constituent un risque majeur. Et en ce qui concerne le stockage géologique de déchets radioactifs, la circulation de fluide dans des réseaux de fractures nouvelles ou réactivées au voisinage de la zone de stockage peut aboutir à la migration de matériaux nocifs. Il est donc important de prévoir les effets de la présence de fractures dans un réservoir. Le but de cette thèse est le développement d’un outil numérique pour la simulation d’un réseau de fractures et de son évolution sous sollicitation hydro-mécanique. Grâce à sa commodité, la méthode des éléments finis étendue (XFEM) sera retenue et associée à un modèle de zone cohésive. La méthode XFEM permet en effet l’introduction de fissures dans le modèle sans nécessairement remailler en cas de propagation des fissures. L’écoulement du fluide dans la fissure et les échanges de fluide entre la fissure et le milieu poreux seront pris en compte via un couplage hydro-mécanique. Le modèle est validé avec une solution analytique asymptotique pour la propagation d’une fracture hydraulique plane dans un milieu poroélastique en 2D comme en 3D. Puis, nous étudions la propagation de fractures hydrauliques sur trajets inconnus. Les fissures sont initialement introduites comme des surfaces de fissuration potentielles étendues. Le modèle de zone cohésive sépare naturellement les domaines adhérents et ouverts. Les surfaces potentielles de fissuration sont alors actualisées de manière implicite par un post-traitement de l’état cohésif. Divers exemples de réorientation de fissures hydrauliques et de compétition entre fissures voisines sont analysés. Enfin, nous présentons l’extension du modèle aux jonctions de fractures hydrauliques / The permeability of rocks is widely affected by the presence of fractures as it establishes prevailing paths for the fluid flow. Natural cracks are then a critical factor for a reservoir productiveness. For low permeability rocks, stimulation techniques such as hydrofracturing have been experienced to enhance the permeability, so that the reservoir becomes profitable. In the opposite, when it comes to geological storage, the presence of cracks constitutes a major issue since it encourages the leak and migration of the material spread in the rock. In the case of CO2 storage, the scenario of leakage across the reservoir seal through cracks or revived faults is a matter of great concern. And as for nuclear waste storage, the fluid circulation in a fracture network around the storage cavity can obviously lead to the migration of toxic materials. It is then crucial to predict the effects of the presence of cracks in a reservoir. The main purpose of this work is the design of a numerical tool to simulate a crack network and its evolution under hydromechanical loading. To achieve this goal we chose the eXtended Finite Element Method (XFEM) for its convenience, and a cohesive zone model to handle the crack tip area. The XFEM is a meshfree method that allows us to introduce cracks in the model without necessarily remeshing in case of crack propagation. The fluid flow in the crack as well as the exchanges between the porous rock and the crack are accounted for through an hydro-mechanical coupling. The model is validated with an analytical asymptotic solution for the propagation of a plane hydraulic fracture in a poroelastic media, in 2D as well as in 3D. Then we study the propagation of hydraulic fractures on non predefined paths. The cracks are initially introduced as large potential crack surfaces so that the cohesive law will naturally separate adherent and debonding zones. The potential crack surfaces are then updated based on a directional criterion appealing to cohesive integrals only. Several examples of crack reorientation and competition between nearby cracks are presented. Finally, we extend our model to account for the presence of fracture junctions