Pour se renseigner et pour acheter des produits les consommateurs vont sur les pages de marques sur Facebook. Ces utilisateurs font du bouche à oreille au sein de communautés de marques sur Facebook. L’objectif de cette thèse est de proposer une analyse des déterminants au bouche à oreille au sein des communautés de marques sur Facebook. Les principaux déterminants au bouche à oreille que nous présentons dans le cadre de cette thèse sont l’engagement à la marque, les motivations au partage, et la disposition générale au partage. / EWord of Mouth (eWOM) influences consumers who belong to brand communities on Facebook. The aim of this thesis is to analyze the antecedents for eWOM on brand communities on Facebook. The thesis proposes an analysis of eWOM’s antecedents: brand engagement, identification to brand communities, general tendency to share on Facebook.

La thèse a pour corpus le Canon des Poèmes (Shijing 詩經). S’appuyant sur l’étude de diverses interprétations des Poèmes, elle propose une application numérique pour interpréter le procédé incitatif (Xing 興) qui, jusque-là, n’a pas été clairement défini au plan rythmique. On y considère le sinogramme à la manière d’une note musicale, ou plus précisément d’une sino-syllabe (hanzi -yinjie 漢字-音節). La thèse peut alors se baser sur trois principes d’analyse :1. le contraste d’effectifs des sino-syllabes ;2. le calcul combinatoire des sino-syllabes ;3. l’échelle de perception des sino-syllabes.En suivant ces trois principes d’analyse, la thèse expose une méthode de lecture expérimentale qui distingue deux aspects dans les sino-syllabes aperçues pendant la lecture : les composants fondamentaux (qui créent le rythme dominant) et les idées principales (qui constituent le sujet thématique). La collaboration harmonieuse du rythme dominant avec le sujet thématique, révèle un système de figuration qui transforme les propos en langage poétique sur le plan rythmique. Ce système sert de référence interne au corpus pour objectiver le rapport rythmique entre le motif initial (= premier groupe de vers) et le motif final (= deuxième groupe de vers), dans la strophe créée par le procédé de Xing. On propose de nommer ce système « figure rythmique » (jielü xiuci 節律修辭). / The present thesis has chosen to center its attention on The Book of Odes. After the study of diverse interpretations of these poems, the thesis is brought to suggest using numeric humanities in the analysis of what Xing calls “the incentive process” whose rhythmic figure has not yet been clearly defined. At the condition that a sinogram could be considered a musical note, it would be more accurate to call it a “sino-syllable” (hanzi-yinjie 漢字-音節). The results are relayed by these three principles: 1. The contrast of occurrences between sino-syllables 2. The sino-syllables’ combinatorial calculation 3. The scale of perception in comprehending sino-syllablesApplying these principles of analysis, the thesis then reveals an experimental method of reading which sets forth two aspects: the basic components that create the dominant rhythm and the mean ideas that constitute the thematic subject.The harmonious relationship of the dominant rhythm with the thematic subject reveals a system that rhythmically transforms words to poetic language. This metaphorical system, called “rhythmic figure” (jiulü xiuci 節律修辭) for the purpose of this research, serves as an internal reference to the corpus to manifest the rhythmic relationship between the first and final patterns in verses produced with Xing’s method.

Lorsque l’extension de la croûte sous les parties hautes des chaînes de montagnes est colinéaire à la direction de convergence, il est traditionnellement admis que le moteur est l’effondrement gravitaire. Pourtant, des études récentes remettent en cause ce paradigme en montrant que l’érosion induit un soulèvement et de l’extension dans la partie centrale des chaînes de montagne à faible taux de convergence. L’objectif de notre étude est d’étudier l’impact de la dénudation de la topographie sur le régime sismo-tectonique des chaînes de montagnes.La première partie de ce travail présente une compilation de données dans les chaînes de montagnes afin de dégager des relations entre régime sismo-tectonique et érosion. Sur la base de ces observations, un modèle cinématique simple permettant de prévoir le régime de la chaîne est proposé. Ainsi, pour les chaînes à faible taux de convergence et d’élévation moyenne, ce modèle prédit de l’extension lorsque le taux de dénudation est 15% plus élevé que le taux de convergence.La deuxième partie est consacrée au développement d’un modèle thermo-mécanique 2D en éléments finis pour étudier l’impact des processus de surface sur la déformation des Pyrénées. Les résultats montrent que la réponse isostatique à l’érosion permet de réactiver des structures pré-existantes. La cinématique d’un plan de faille hérité peut être prédite grâce au gradient du profil des vitesses de surface horizontales. Ainsi, un plan situé dans la zone d’érosion est réactivé en faille normale alors qu’en bordure de cette zone une faille est réactivée en régime inverse. Ces résultats suggèrent que la déformation actuelle des Pyrénées pourrait être la conséquence d’un processus d’érosion.Compte tenu du faible nombre d’études quantifiant les taux d’érosion dans les Pyrénées, les modèles développés dans la deuxième partie souffrent d’une forte incertitude. Pour y remédier, nous avons cherché à les quantifier dans les Pyrénées Centrales grâce à une étude qui combine deux types de données : taux de dénudation des bassins versants à partir des isotopes cosmogéniques, et vitesses d’incision à partir des sédiments piégés dans les karsts). Ces résultats sont présentés dans la troisième partie. Les profils de dénudation obtenus sont compatibles avec un rejeu en faille normale d’un plan situé dans la Zone Nord Pyrénéenne des Pyrénées Centrales. Dans les Alpes, une bonne corrélation apparaît entre la valeur du taux d’érosion et la vitesse verticale géodésique, ce qui pose la question de l’impact de la déglaciation tardi-Wurmienne dans les Alpes sur la déformation actuelle. Un modèle numérique détaillant cette relation est présenté dans le quatrième chapitre. Les résultats montrent que la déglaciation des Alpes occidentales est contrôlée par l’hétérogénéité rhéologique de la croute. Certains de nos modèles prédisent des vitesses de surrection compatibles avec celles mises en évidence par la géodésie. / When mountain ranges upper parts express crustal extension direction collinear to the convergence direction, it is traditionally accepted that the extensive motor is gravitational collapse. However, recent studies challenge this paradigm by showing that erosion induces uplift and extension in the central part of the low convergent mountain ranges. Our goal is to investigate the impact of the denudation on the seismotectonic regime of mountain ranges.In order to identify a relationship between seismotectonic regime and erosion, the first part of this work presents a compilation of data in the mountain ranges. Based on these observations, a simple kinematic model is proposed to predict the seismotectonic regime of the chain. Thus, for low convergence rate chains with a moderate mean elevation, this model predicts an extension regime when the denudation rate is 15% higher than the convergence rate.The second part is devoted to the development of thermomechanical 2D finite element model to study the impact of surface processes on the deformation of the Pyrenees. The results show that the isostatic response to erosion reactivates pre-existing structures. The kinematics of an inherited fault plane can be predicted due to the gradient of the horizontal surface velocities profile. Thus, a plane located in the eroded zone is reactivated in normal fault when in a border area of this same plane is reactivated in reverse fault. These results suggest that the current deformation in the North Pyrenean Zone could be the result of surface processes.Given the small number of studies quantifying erosion rates in the Pyrenees, the models developed in the second part suffer from high uncertainty. To remedy this, we sought to quantify it in the central Pyrenees through a study that combines two types of data: watershed denudation rates from cosmogenic isotopes concentration, and incision rates from sediments buried in the karst. These results are presented in chapter 3. Denudation profiles obtained are consistent with a replay of a normal fault plane located in the North Zone of the Central Pyrenees.In the Alps, a good correlation appears between the value of the rate of erosion and geodetic vertical velocities, which raises the question of the impact of the late-würmian deglaciation in the Alps on the present deformation. A numerical model detailing this relationship is presented in the fourth chapter. The results show that deglaciation of the western Alps is controlled by the rheological heterogeneity of the crust. Some of our models predict uplift rates consistent with those highlighted by geodesy.

Dans un disque axisymétrique en équilibre, les vitesses galactocentriques radiales et verticales sont théoriquement nulles. Pourtant, de grands relevés spectroscopiques ont révélé que les étoiles du disque de la Voie Lactée sont animées de vitesses non nulles dans les directions radiale et verticale. Les structures en vitesse radiale sont généralement associées aux composantes non-axisymétriques du potentiel. Celles en vitesse verticale non nulle sont souvent associées à des excitations externes. Nous avons montré que la réponse stellaire à une perturbation spirale produit un déplacement radial et des mouvements verticaux non nuls. La structure du champ moyen de vitesse obtenue est cohérente avec les observations. De plus un modèle simple reposant sur une linéarisation des équations d'Euler reproduit naturellement ce résultat. Nous concluons que ces structures observées pourraient aussi être engendrées par des perturbations internes non-axisymétriques. / In an equilibrium axisymmetric galactic disc, the mean galactocentric radial and vertical velocities are expected to be zero everywhere. Recent spectroscopic surveys have however shown that stars of the Milky Way disc exhibit non-zero mean velocities outside of the Galactic plane in both the radial and vertical velocity components. While radial velocity structures have already often been assumed to be linked with non-axisymmetric components of the potential, non-zero vertical velocity structures are usually rather attributed to excitations by external sources. We show that the stellar response to a spiral perturbation induces both a radial velocity flow and non-zero vertical motions. The resulting structure of the mean velocity field is qualitatively similar to the observations. Such a pattern also emerges from an analytic toy model based on linearized Euler equations. In conclusion, non-axisymmetric internal perturbations can also be the source of the observed mean velocity patterns.

Il est difficile de quantifier les dissipations d'énergie se produisant dans les ouvrages lors d'un séisme de manière réaliste, en particulier dans le cas de structures en béton armé. Ceci s'explique par la diversité des sources de dissipations. Classiquement, elles sont introduites dans la modélisation des structures sous la forme d'un amortissement global peu physique. Le problème qui se pose est de savoir comment prendre en compte l'amortissement de manière plus physique. Ces travaux visent à apporter des éléments de réponse à cette problématique. Pour cela, deux objectifs ont été poursuivis : le premier consiste à qualifier et quantifier expérimentalement les sources d'amortissement dues au béton, le second vise à développer une méthode de calcul, peu couteuse et permettant de modéliser le comportement global mais aussi l'amortissement de manière réaliste. Une campagne d'essais de flexion 3 points alternée a été réalisée afin de quantifier les mécanismes contribuant à l'amortissement. Cette campagne est effectuée sur des poutres saines mais également sur des poutres pré-endommagées. L'analyse des expériences a permis d’identifier différentes grandeurs physiques influençant l’amortissement. Cette campagne a été modélisée à l'aide de trois lois de comportements, une loi s'est révélée pertinente. La loi de comportement ainsi identifiée a permis de développer une modélisation simplifiée, en vue d'études probabilistes. Cette dernière est fondée sur une loi de comportement simple couplée à une actualisation de l'amortissement visqueux. Trois actualisations de l'amortissement ont été développées et confrontées à une campagne d'essais sismiques mono axiaux sur un poteau. / It is hard to determine the energy dissipation that occurs during an earthquake, especially where reinforced concrete structures are concerned. The reason for this is the many different causes of energy dissipation, these dissipations typically creep into the essential pattern of the structures as a uniform, slight damping, and which is heavily quantify. The challenge is therefore to ascertain how to carry out damping in a way that relies more on the laws of physics themselves. This study aims at bringing some clarifications to this problem. In order to achieve this, two objectives were targeted during the case study: the first consisted in experimentally qualifying and quantifying the sources of damping in concrete, the second aims at developing a method which model both the overall behaviour and the damping in a realistic way with low computational costs. Reverse 3-point bending tests were carried out to determine and quantify the mechanisms responsible for damping. The tests were carried out on sound beams, and also on pre-damaged beams. It was possible to relate the damping to the damage, the intensity of the load and the erosion of the crack surfaces. These tests were modeled using three physical constitutive laws, one proved to be relevant. The behaviour law thus identified allowed us to develop a simplified model to be used during probabilistic assessments. This model is based on a simple constitutive law, coupled with the updating of viscous damping. It is carried out according to the evolutions of the mechanical properties of the structure and the load. Three updates of the damping were developed and subjected to a serie of mono axial seismic tests on a column.

Les astéroïdes Cérès et Vesta ont motivé la mission spatiale Dawn parce qu'ils représentent deux embryons planétaires différents restés relativement intacts depuis leur formation. Vesta est large- ment considéré comme le corps parent des météorites HED témoins d'une activité magmatique probablement due à la présence de l'isotope radioactif 26Al qui était suffisamment abondant pour permettre la fusion interne des corps rocheux primitifs. La composition d'une surface planétaire peut être mesurée grâce à l'analyse des rayons gammas qu'elle produit. Pour la sonde Dawn cela est rendu possible par l'instrument GRaND et la scintillation d'un cristal de BGO. Cette thèse présente l'analyse des spectres gammas de Vesta par deux outils de séparation aveugle de source: l'analyse en composantes indépendantes (ICA) et la factorisation en matrice non-négative (NMF). Ces méthodes sont aussi appliquées à un jeu de données lunaire comparable et déjà bien interprété. Des spectres synthétiques lunaires permettent de tester ICA et NMF. La séparation de spectres élémentaires s'avère délicate même si on peut distinguer les éléments K, Th et Fe en raison des propriétés statistiques de leur signaux sources plus favorables. On mesure la sensibilité d'ICA-NMF à la variabilité chimique de la surface pour des Lunes artificielles, ce qui permet d'expliquer l'absence de séparation d'un signal élémentaire clair dans le cas de Vesta. Malgré les observations de la sonde Dawn et le nombre important d'informations fournies par les HED, il n'y a pas de consensus sur la formation des HED. On met souvent en avant l'existence d'un océan magmatique global sur Vesta, alors que la migration de la principale source chaleur, contenue dans le premier minéral fondu, le plagioclase, ne permet pas la fusion totale. On met en oeuvre un modèle de migration des magmas, basé sur les équations de la compaction. On adapte ce modèle en utilisant un diagramme d'équilibre de phase olivine-anorthite-quartz. Cela permet de calculer l'évolution de la minéralogie en fonction du temps et de la profondeur. Les résultats montrent que les eucrites et les diogénites pourraient être une caractéristique commune des gros corps accrétés tôt dans l'histoire du système solaire. / Asteroids Vesta and Ceres motivated the space mission Dawn because they represent two different planetary embryos that remained relatively intact since their formation. Vesta is broadly considered as the parent body of the HED meteorites suite that are witnesses of a magmatic activity probably due to the presence of the radioactive isotope 26 Al which was present in significant amount to cause internal melting of primitive rocky bodies. The composition of a planetary surface can be quantified through the analysis of the gamma rays it produces. This is made possible for the Dawn spacecraft by the instrument GRaND and the scintillation of a BGO crystal. This thesis presents the analysis of gamma ray spectra from Vesta by two blind source separation methods: the independent component analysis and the non negative matrix factorization. These methods are also applied to an equivalent lunar dataset already well interpreted. Lunar synthetic spectra are used to test ICA and NMF. The separation of elementary spectra is delicate although K, Th and Fe can be discriminated due to the more favorable statistical properties of their source signals. The sensitivity of separation to the chemical variability is assessed based on artificial lunar spectra, which allows to explain the lack of separation of a clear elemental signal in the case of Vesta. Despite the observations of Dawn and the important collection of HED data, there is no consensus on the conditions of the vestan magmatism. A global magma ocean is often put forward, whereas the migration of the heat source, contained in the easiest mineral to melt, plagioclase, does not allow it. A model of melt migration is implemented, based on two-phase flow equations. This model is combined with the olivine-anorthite-quartz equilibrium phase diagram. This allows to predict the mineralogy as a function of depth and time. Results obtained show that eucrites and diogenites may be a common feature of large bodies accreted early in solar system history.

La thèse contient 6 chapitres. Le premier chapitre contient une introduction à la régression linéaire et aux problèmes Lasso et Lasso bayésien. Le chapitre 2 rappelle les algorithmes d’optimisation convexe et présente l’algorithme FISTA pour calculer l’estimateur Lasso. La statistique de la convergence de cet algorithme est aussi donnée dans ce chapitre en utilisant l’entropie et l’estimateur de Pitman-Yor. Le chapitre 3 est consacré à la comparaison des méthodes quasi-Monte Carlo et Monte Carlo dans les calculs numériques du Lasso bayésien. Il sort de cette comparaison que les points de Hammersely donne les meilleurs résultats. Le chapitre 4 donne une interprétation géométrique de la fonction de partition du Lasso bayésien et l’exprime en fonction de la fonction Gamma incomplète. Ceci nous a permis de donner un critère de convergence pour l’algorithme de Metropolis Hastings. Le chapitre 5 présente l’estimateur bayésien comme la loi limite d’une équation différentielle stochastique multivariée. Ceci nous a permis de calculer le Lasso bayésien en utilisant les schémas numériques semi implicite et explicite d’Euler et les méthodes de Monte Carlo, Monte Carlo à plusieurs couches (MLMC) et l’algorithme de Metropolis Hastings. La comparaison des coûts de calcul montre que le couple (schéma semi-implicite d’Euler, MLMC) gagne contre les autres couples (schéma, méthode). Finalement dans le chapitre 6 nous avons trouvé la vitesse de convergence du Lasso bayésien vers le Lasso lorsque le rapport signal/bruit est constant et le bruit tend vers 0. Ceci nous a permis de donner de nouveaux critères pour la convergence de l’algorithme de Metropolis Hastings. / The thesis contains 6 chapters. The first chapter contains an introduction to linear regression, the Lasso and the Bayesian Lasso problems. Chapter 2 recalls the convex optimization algorithms and presents the Fista algorithm for calculating the Lasso estimator. The properties of the convergence of this algorithm is also given in this chapter using the entropy estimator and Pitman-Yor estimator. Chapter 3 is devoted to comparison of Monte Carlo and quasi-Monte Carlo methods in numerical calculations of Bayesian Lasso. It comes out of this comparison that the Hammersely points give the best results. Chapter 4 gives a geometric interpretation of the partition function of the Bayesian lasso expressed as a function of the incomplete Gamma function. This allowed us to give a convergence criterion for the Metropolis Hastings algorithm. Chapter 5 presents the Bayesian estimator as the law limit a multivariate stochastic differential equation. This allowed us to calculate the Bayesian Lasso using numerical schemes semi-implicit and explicit Euler and methods of Monte Carlo, Monte Carlo multilevel (MLMC) and Metropolis Hastings algorithm. Comparing the calculation costs shows the couple (semi-implicit Euler scheme, MLMC) wins against the other couples (scheme method). Finally in chapter 6 we found the Lasso convergence rate of the Bayesian Lasso when the signal / noise ratio is constant and when the noise tends to 0. This allowed us to provide a new criteria for the convergence of the Metropolis algorithm Hastings.

Les suspensions colloïdales se trouvent un peu partout autour de nous, dans les matériauxde constructions, en cosmétique, dans l'alimentation, en biologie. Elles sont composésde particules nanométriques ou micrométriques dispersés dans un gaz, un liquide ou unsolide.Cette thèse porte sur les suspensions colloïdales dans des solutions ioniques,où les colloïdes portent une charge électriques, par exemple des particules de silicedans une solution aqueuse de chlorure de sodium, à un pH basique. Les colloïdes,ici approximés par des sphères, peuvent varier significativement en taille,ce qui peut avoir un effet important sur le comportement de ces systèmes.Cette étude vise à améliorer la compréhension de ces suspensions colloïdales chargéespar des modèles théoriques résolus par des simulations numériques.Un des défis de ces simulations est le grand nombre de degrés de libertés. Pour chaque(micro-)ion il y a des centaines de molécules de solvant, et pour chaque colloïdedes centaines voire des milliers d'ions. Pour s'en sortir, nous avons calculéles interactions effectives à l'échelle colloïdale. Nous avons repris et développéplusieurs approches, présentant chacune un compromis en terme de temps de calcul etprécision.La variation en taille des colloïdes peut introduire des effets intéressants,observés expérimentalement, notamment le fractionnement des suspensions en plusieursphases cristallines quand on augmente la concentration en colloïdes.Des techniques de simulations Monte-Carlo simples associées aux interactions inter-colloïdescalculées précédemment ont permis d'obtenir des résultats en bon accord avec l'expérience. / Colloidal suspensions are found a bit everywhere around us, in construction materials,in cosmetics, in food, in biology. They are composed of nanometric or micrometric particlesdispersed in a gas, a liquid or sometimes a solid.This thesis is about colloidal suspensions in ionic solutions, where colloids bear anelectric charge, for example silica particles in an aqueous solution of sodium chloride,at a basic pH. The colloids, here approximated by spheres, can vary significantly in size,which can have an important effect on the behavior of these systems.This study aims at improving the understanding of these charged colloidal suspensionsby theoretical models solved by numerical simulations.of these charged colloidal suspensionsby theoretical models solved by numerical simulations.One of the challenge of theses simulations is the huge number of degrees of freedom.For each (micro-)ion there is hundreds of solvent molecules, and for each colloidthere is hundreds if not thousands of ions. To get away with it, we calculated theeffective interactions at the colloidal scale. We took and developed several approaches,each showing a trade-off in terms of computational time and accuracy.The size variation of colloids may introduce interesting effects, experimentallyobserved, notably the fractionation of suspensions in several crystalline phaseswhen the colloidal concentration is increased. Some simple Monte-Carlo simulationtechniques in combination with the inter-colloid interactions computed previouslyallowed us to obtain results in good agreement with experiments.

Les pièces de structure aéronautique utilisées dans l’aéronautique sont fabriquées en plusieurs étapes. Des étapes, comme le traitement thermique, génèrent des contraintes résiduelles. Les enlèvements de matière réalisés par les opérations d’usinage peuvent alors conduire à la réorganisation des contraintes résiduelles dans la matière et ainsi à la déformation de la pièce. La mesure in-situ de ces déformations devient nécessaire lorsque ce phénomène est étudié. Ces travaux de thèse abordent cette problématique dans le cadre du projet SIMP-Aero. L’objectif de ces travaux de thèse est d’adapter la méthode de corrélation d’images numériques à la mesure de la déformation de pièce pendant l’usinage, c’est-à-dire dans un centre d’usinage. Pour cela, plusieurs améliorations sont apportées à la méthode. Premièrement, les mouvements du système optique sont pris en compte afin que ceux-ci n’altère pas la qualité des mesures. Ensuite, les copeaux présents sur les images sont détectés et filtrés par un algorithme. Au final, la méthode développée permet de mesurer des champs de déplacement durant toute la séquence d’usinage, sans devoir l’interrompre, avec une incertitude de mesure de l’ordre du centième de millimètre. / Structural aluminum alloy parts used in aeronautics are manufactured in several steps, from forming processes and heat treatments to final machining. Some of the process steps induce residual stresses. The material removal during machining release these residual stresses and thus, leads to the part deformation. The in-situ measurement of these deformations becomes necessary when this phenomenon is studied. The present work address this problematic in the context of the ANR SIMP-Aero Project. It aims to define a reliable experimental technique dedicated to the measurement of part deformations during machining of large aeronautical parts. The backbone of the technique relies on Digital Image Correlation (DIC). Mainly as a consequence of the harsh constraints environment of machining, the customization of DIC is required. First, movements of the optical system are quantified and compensated for the proper measurement of the workpiece displacement. Then, the metal chips that fly between the observed surface and the acquisition system are detected and filtered by the algorithm. Finally, the developed method enables the measurement of displacement fields throughout the whole machining sequence, without interrupting it, and a measurement uncertainty of around one hundredth of a millimeter is ensured.

L’accent mis actuellement par l’industrie pétrolière sur l’augmentation de l’efficacité des réservoirs, ainsi que sur l’intérêt grandissant pour l’exploitation d’autres sources d’énergie enfouies profondément sous terre a suscité un regain d’intérêt pour la mécanique de la fracturation des roches en général et la fracturation hydraulique en particulier. La fracturation hydraulique, appelée de manière informelle “fracturation”, est un processus qui consiste généralement à injecter de l’eau, sous haute pression dans une formation rocheuse via le puits. Ce processus vise à créer de nouvelles fractures dans la roche et à augmenter la taille, l’étendue et la connectivité des fractures existantes. Des avancées récentes dans la modélisation et la simulation de fractures hydrauliques ont eu lieu, au confluent de facteurs qui incluent une activité accrue, une tendance vers une complexité accrue et une compréhension approfondie du modèle mathématique sous-jacent et de ses défis intrinsèques. Cependant, certaines des caractéristiques très importantes de ce processus ont été négligées. Parmi les caractéristiques négligées, on peut citer l’incapacité de la grande majorité des modèles existants de s’attaquer à la fois à la propagation de fractures hydrauliques dans la roche intacte, à l’inititation de nouvelles fractures ainsi qu’à la réactivation des fractures existantes. Une autre caractéristique qui a été ignorée est sa dimension intrinsèque en trois dimensions, négligée par la plupart des modèles actuallement proposés. Parmi tous les différents types de méthodes numériques développées pour évaluer le mécanisme du phénomène de fracturation, très peu sont capables de représenter la totalité des mecanismes mis en jeu. Dans la présente thèse, l’initiation et la propagation de fissures induites par les fluides dans des roches isotropes transversales sont simulées à l’aide d’un modèle hydromécanique (HM) couplé basé la méthode XFEM (eXtended Finite Element Method) et un modèle de zones cohésives. Le HM-XFEM développé dans cette thèse est une extension des modèles précédemment développés dans l’équipe hydro-géomécanique multi échelle de GeoRessources. L’accent a été porté plus particulièrement sur la prise en compte de l’anisotropie du milieu et sur son influence sur le chemin de propagation. Ce dernier est défini à partir du le concept d’angle de bifurcation introduit auparavant dans la littérature. En complément des développements réalisés dans le modèle HM-XFEM, effort a été fait pour mieux comprendre l’initiation de la fissure en utilisant la méthode des éléments discrets (DEM) à l’aide du logiciel open source YADE Open DEM. La nature différente des deux méthodes, DEM étant une méthode discontinue et XFEM, une méthode continue, révèle les potentiels des deux méthodes et permet de comparer correctement la méthode qui convient le mieux au problème à résoudre, compte tenu des objectifs de la conception / Current emphasis in petroleom industry toward increasing the reservoirs efficiency, along with the interest in exploitation of other sources of energy buried deep underground created a renewed interest in rock fracture mechanics in general and hydraulic fracturing specifically. Hydraulic fracturing, informally referred to as “fracking,” is an oil and gas well development process that typically involves injecting water, under high pressure into a bedrock formation via the well. This process is intended to create new fractures in the rock as well as increase the size, extent, and connectivity of existing fractures. However some of the very important features of this process have been overlooked. Among these neglected features one can name of inability of the vast majority of existing models to tackle at once the propagation of hydraulic fractures in fractured rocks-masses where a competing dipole mechanism exists between fracturing of the intact rock and re-activation of exiting fracture networks. Another feature that has been ignored is its intrinsically three dimensionality which is neglected by most models. Among all different types of numerical methods that have been developed in order to assess the mechanism of fracturing phenomenon very few, if any, can handle the entire complexity of such process. In the present thesis, fluid-driven crack initiation and propagation in transverse isotropic rocks is simulated using a coupled model comprising of eXtended Finite Element Method (XFEM) and cohesive zone models. The HM XFEM developed in this thesis is an extension to previous models developed introduced in multiscale hydrogeomechanics team of GeoRessources. An emphasis is put on considering the anisotropic nature of the medium and on studying its influence on the propagation path. This latter is investigated by the concept of bifurcation angle previously introduced in literature. In complementary efforts was made to have a better understanding of crack initiation in transversely isotropic media, we also used the discrete element method (DEM) in order to gain insights into the mechanisms at stake. Both methods exhibit their advantages and disadvantages in modeling fracturing phenomenon. The different nature of two methods, DEM being a discontinuous and XFEM being a continuous method, reveals potentials of both methods and renders a good comparison of which method suits the problem in hand the best, considering the the objectives of the design