L'utilisation des moyens de calcul de type éléments finis couplés aux progrès techniques en matière de mesures de champs expérimentaux, tant microstructuraux que mécaniques, ouvrent la voie à une richesse de dialogue expérience/calcul qui ne se limite plus à la seule réponse macroscopique obtenue sur un volume élémentaire représentatif. L'objectif de ce travail concerne l'exposé de la méthodologie adoptée et des outils développés pour l'identification de paramètres pertinents de lois de comportement de matériaux hétérogènes à partir de la définition de leur microstructure et des mécanismes responsables, au sein du volume élémentaire représentatif, de leur comportement élasto-plastique voire de leur endommagement. Nous nous intéresserons tout particulièrement à préciser la notion d'hétérogénéités, qui est fortement dépendante de la base de mesure choisie, que se soit pour l'expérimentation ou pour les simulations numériques. Nous aborderons d'autre part, les questions relatives à la pertinence de mesures de surface pour décrire le comportement d'un volume de matière puis nous nous focaliserons sur l'intérêt d'utiliser des mesures cinématiques, obtenues en cours d'essais mécaniques in situ comme conditions aux limites des simulations numériques. Nous discuterons de la mise en place de cette méthodologie à partir de résultats obtenus à partir d'alliages métalliques et de matériaux virtuels.

Nous avons caractérisé de manière détaillée les  fortes hétérogénéités de texture et de microstructure (macrozones) par EBSD dans des zones déformées différemment d’un disque en IMI834 de structure bimodale αP/αS, utilisé en motorisation aéronautique. Des échantillons spécifiques ont permis de montrer que l’amplitude des signaux ultrasonores rétrodiffusés dans ces échantillons dépendait de la présence et de l’organisation des macrozones. Partant de ces observations, une analyse ultrasonore des signaux rétrodiffusés dans les parties du disque soumis à de fortes déformations lors du forgeage a permis de sélectionner un jeu d’éprouvettes ayant des caractéristiques microstructurales différentes et destinées à des essais en fatigue-dwell. Une forte dispersion des durées de vie liée aux microstructures présentes a été observée. L’analyse fractographique des surfaces de rupture et leurs études par EBSD ont permis de localiser les sites d’amorçage et les chemins de propagation des fissures ainsi que leurs caractéristiques cristallographiques. Les résultats des essais en fatigue-dwell montrent que la réduction de la tenue en service est fortement influencée par la présence de macrozones formées d’agglomérats de grains αP allongés associés à des lamelles αS plus épaisses et organisées en colonies et d’autre part présentant majoritairement des axes c à moins de 30° de l’axe de sollicitation / We have characterized in detail the strong heterogeneity of texture and microstructure (macrozones) by EBSD in areas, differently deformed, of a disk in IMI834 with a bimodal structure αP/αS, used in aircraft engines. Specific samples have shown that the amplitude of back-scattered ultrasonic signals depended on the presence and organization of macrozones. Based on these observations, an analysis of back-scattered ultrasonic signals in parts of the disc subject to strong deformation during forging led to select a set of specimens with different microstructural features for fatigue-dwell tests. A high dispersion of lifetimes associated with these microstructures was observed. The fractographic analysis of fracture surfaces and their EBSD studies were used to locate crack initiation sites and propagation paths and their crystallographic characteristics. The results of fatigue-dwell tests showed that the lifetime reduction is strongly influenced by the presence of macrozones formed by agglomerates of grains elongated αP associated with thicker αS lamellae, organized in colonies and with c-axis predominantly at less than 30 ° of the solicitation axis

Il est crucial de garantir l'intégrité des cuves de réacteurs à eau pressurisée (REP) en cas de fonctionnement accidentel : dans ce contexte, la compréhension et la modélisation des mécanismes de rupture fragile des aciers constituent des éléments décisifs de l'évaluation complexe des durées de vie des cuves.<br /><br />Les modèles d'approche locale de la rupture par clivage constituent l'un des principaux outils de prédiction de la ténacité des aciers faiblement alliés. La dispersion des contraintes à rupture est interprétée comme un effet de la distribution des défauts dans la microstructure, mais l'effet des hétérogénéités mécaniques n'est pas pris en compte. Or, en dessous d'une température de transition de comportement Ta (de l'ordre de 25°C), les mécanismes de déformation sont grandement affectés par la température et la vitesse de déformation.<br /><br />Notre approche consiste à prendre en compte l'effet des hétérogénéités de contraintes dans un critère local d'amorçage du clivage. Les résultats de calculs de microstructure sont utilisés pour proposer une description statistique de l'évolution des distributions de contraintes locales. Cette approche statistique permet de proposer un modèle d'approche locale de la rupture dépendant à la fois des hétérogénéités mécaniques et des distributions de tailles de défauts.<br /><br />Le comportement du matériau et son évolution sont caractérisés aux échelles microscopique et macroscopique dans le domaine de température [25°C,-196°C]. Des essais de traction simple, de sauts de vitesse et de température, et de ténacité sont réalisés.<br /><br />Nous proposons un modèle de comportement micromécanique décrivant le comportement plastique en dessous de la température de transition Ta. La loi de comportement est basée sur les mécanismes de déformation décrits dans la bibliographie et identifiées par méthode inverse à partir des essais mécaniques. Les observations au MET et la caractérisation du comportement activé thermiquement permettent de fixer plusieurs paramètres du modèle.<br /><br />Des simulations sont réalisées afin de modéliser les distributions de contrainte principale Σ1 dans deux microstructures bainitiques correspondant au volume élémentaire de l'approche locale de la rupture. L'effet de la température et de la triaxialité sur l'évolution des hétérogénéités est caractérisé. Nous proposons une fonction de distribution décrivant la distribution des valeurs locales de Σ en fonction des contraintes principales et équivalente <Σ> et <Σmises> moyennes dans la microstructure.<br /><br />Cette fonction est utilisée pour formuler un modèle d'approche locale de la rupture intégrant la distribution des tailles de défauts critiques et les distributions de Σ1.. On montre que dans certains cas, la dispersion des contraintes locales suffit à expliquer les dispersions des contraintes à rupture à l'échelle du volume élémentaire. Les dispersions de contraintes à rupture sont en accord avec celles prédites par le modèle de Beremin. La prise en compte des hétérogénéités mécaniques permet d'introduire une dépendance de la probabilité de rupture en fonction de la température, de la déformation et de la triaxialité. Il reste à appliquer le modèle d'approche locale au calcul d'éprouvettes CT et de comparer les dispersions de ténacités simulées à celles mesurées expérimentalement.

Ce travail est consacré à l’étude par simulation de dynamique moléculaire l’effet d’une perturbation contrôlée, utilisant des molécules photoisomérisables sur les effets de taille finie et de surface dans les matériaux amorphes. Le présent modèle utilise le potentiel de Lennard-Jones et l’algorithme de Verlet pour résoudre les équations du mouvement. Nous avons réalisé plusieurs boîtes de simulation de tailles différentes dans le but d’étudier l’effet de la taille de la boîte sur les propriétés du matériau. Nous avons utilisé ces boîtes pour simuler le matériau à différentes températures, le matériau était avec un chromophore dilué à l’intérieur, le chromophore s’isomérisant ou non,et avec une période courte ou longue. L’étude utilise des simulations à l’échelle atomique par dynamique moléculaire où l’énergie produite par le processus de photo-isomérisation est absorbée par un thermostat,après sa dégradation thermique à l’intérieur de la matrice hôte. Par une modification continue de la forme des photochromes trans vers cis et inversement, une photo isomérisation est introduite de manière périodique. Dans ces conditions, nous avons montré que ce modèle reproduit bien les propriétés statiques et dynamiques d’un matériau amorphe. Nous avons trouvé que des hétérogénéités dynamiques apparaissent dans un liquide surfondu pour une température qui égale 120K. Nous avons constaté que le temps de relaxation diminue en même temps que la force des mouvements coopératifs lorsque la taille du système est devient plus petite. Ces résultats suggèrent la présence de deux mécanismes physiques concurrents différents avec différentes échelles de longueur. / This work is devoted to the study by molecular dynamics simulation the effect of a controlled disturbance, using photo-isomerizable molecules on the finite size effects and surface in amorphous materials. This model uses the potential of Leonard Jones and the Verlet algorithm to solve the equations of motion.We realized several simulation boxes of different sizes in order to study the effect of the size of the box on the material properties. We used these boxes to simulate the material at different temperatures, the material was diluted with a chromophore inside, the chromophore isomerized or not, and with a short or long period. The study uses simulations at the atomic scale molecular dynamics where the energy produced by the photoisomerization is absorbed by a thermostat after thermal degradation within the host matrix. By a continuous change in the form of photochromes trans tocis and vice versa, a photo isomerization is introduced periodically. Under these conditions, we have shown that this model reproduces well the static and dynamic properties of amorphous material. We found that dynamic heterogeneities appear in a supercooled liquid to a temperature equal to 120 K. We found that the relaxation time decreases along with the strength of cooperative movements when the system size is becoming smaller. These results suggest the presence of two different physical mechanisms competitors with different length scales

La présente étude porte sur les interactions entre deux types de matériaux susceptibles d'être mis en jeu au sein d'un stockage géologique de déchets radioactifs : les matériaux métalliques, constituant notamment le surconteneur de déchets de haute activité, et les matériaux argileux, telle que la roche hôte argileuse. Les interactions entre ces deux matériaux en conditions de stockage pourraient en effet modifier leurs propriétés de confinement initiales. Les travaux présentés ont visé à déterminer l'influence d'hétérogénéités (vides technologiques et fractures) et d'activités bactériennes sur ces interactions, notamment en termes d'évolution de propriétés chimiques et hydrauliques de l'argile. Dans cet objectif, deux expériences intégrées en cellules de percolation ont été mises en œuvre, à 60 °C, pendant 13 mois : la première en présence de deux souches bactériennes (BSR, BFR), la seconde constituant un témoin abiotique. Ces expérimentations ont permis la circulation d'une eau synthétique dont la composition est représentative de celle de l'eau porale de Tournemire au travers d'un compact de fer pulvérulent puis d'une carotte d'argilite du Toarcien de Tournemire, artificiellement fissurée. L'une des deux demi-carottes d'argile contenait également un cylindre de fer massif. Les caractérisations post-mortem (MEB, MEB/EDS, Raman, DRX, tomographie aux rayons X) ont permis d'étudier deux interfaces : l'interface fer pulvérulent/argilite et l'interface fer massif/argilite.Concernant le déroulement des expérimentations, il convient de noter en premier lieu que la fissure a probablement été circulante durant la totalité de l'essai, ce que tendent à confirmer les modélisations couplées chimie/transport avec le code HYTEC. Toutefois, aucune phase néoformée n'a été mise en évidence à son niveau. Par ailleurs, la survie bactérienne dans la cellule biotique a été confirmée au cours de l'expérimentation par le suivi régulier de la population et par une analyse de la diversité génétique à la fin des essais. Une nette diminution de la concentration en sulfates en sortie de cellule confirme l'activité des BSR.Les caractérisations solides en fin d'essais ont mis en évidence dans la zone de fer pulvérulent, avec et sans bactéries, de la magnétite et de la chukanovite, cette dernière étant localisée majoritairement à l'interface avec l'argilite. Une zone enrichie en fer (10 µm) a été identifié à l'interface dans l'argilite. La vitesse moyenne de corrosion du fer pulvérulent a été estimée à 0,2 µm/an (valeur basse). A l'interface fer massif/argilite, deux faciès ont été mis en évidence. Le premier faciès, identifié dans les deux dispositifs, est composé de chukanovite en couche interne et de sidérite en couche externe. L'extension de la perturbation de l'argilite est de l'ordre de 30 µm. Le second faciès, mis en évidence seulement en présence de bactéries, montre la présence de sulfure de fer (mackinawite) et d'avancées locales de corrosion. Enfin, les simulations HYTEC réalisées ont permis de mieux comprendre les mécanismes biogéochimiques observés expérimentalement, notamment l'effet du pH, et de mieux quantifier certains paramètres cinétiques clefs. / This study focuses on the interactions between two materials that may be introduced in a geological disposal of radioactive waste: metallic materials such as the high-level waste overpack, and clay materials such as the clay host rock. Indeed, the interactions between these two materials in such conditions could induce a change of their initial confinement properties. This work aimed at determining the influence of heterogeneities (technological gaps and fractures) and bacterial activities on these interactions, in terms of evolution of chemical and hydraulic properties of clayey materials. To this end, two percolation cells have been conducted during 13 months: the first one with two bacteria (SRB, IRB), the second one without bacteria. These experiments, carried out at 60°C, involved circulating synthetic water representative of the Tournemire pore water through iron powder and through Toarcian artificially cracked argillite from Tournemire. An iron rod was also placed into the argillite. Thus, solid characterizations (SEM, SEM/EDS, Raman, XRD, X-ray tomography) allowed the study of both interfaces: the iron powder/argillite interface and the iron rod/argillite interface.The water probably circulated into the crack during the entire test, which was confirmed by reactive transport modeling with the HYTEC reactive transport code. However, no secondary phase was identified in the crack. In addition, bacteria survival in the biotic cell was confirmed during the experiment by monitoring their population and by analyzing their genetic diversity at the end of the experiment. A strong decrease in sulfate concentration was measured in the output, which confirms the SRB activity.Solid characterization conducted at the end of the experiments have highlighted, with and without bacteria, the occurrence of magnetite and chukanovite in the iron powder, the latter being mainly located close to the argillite interface. In the argillite, a Fe-enriched zone (10 µm) was identified. The mean corrosion rate was estimated at 0.2 µm/y (lower bound). At the iron rod/argillite interface, two corrosion facies were observed. The first, identified in both cells, is mainly constituted of chukanovite in the inner layer and siderite in the outer layer. Extent of the argillite perturbation reaches about 30 µm. The second, only observed with bacteria, highlights the presence of iron sulfide precipitation (mackinawite) and localized corrosion patterns. Finally, HYTEC simulations have enabled us to better understand the observed biogeochemical processes, such as the pH effect, and to better quantify some key kinetic parameters.

La caractérisation des fissures dans les roches est devenu un thème de recherche majeur pour de nombreuses applications géophysiques dans le domaine de la subsurface ou des réservoirs (entre autres l'exploration sismique, l'étude continue de l'évolution d'une zone de production en pétrole, gaz ou eau, le stockage de déchets toxiques, la géothermie ... ). Cet intérêt pour les fissures provient de ce qu'elles se comportent à la fois comme des réserves de fluide et comme des chemins d'écoulement potentiels. Au cours de ce travail, nous nous sommes consacrés à l'étude des vitesses d'ondes se propageant dans les roches fissurées. Notre objectif est de modéliser, à partir d'une méthode différentielle auto-cohérente, l'influence de la fissuration et de la saturation en fluide sur les propriétés élastiques effectives et, par conséquent, sur les vitesses des ondes. Le modèle construit repose sur une hypothèse de géométrie simplifiée pour les inclusions poreuses: la porosité dite rigide est représentée par des sphères tandis que la porosité de type compressible (les fissures) est décrite par des ellipsoïdes de révolution. Des distributions isotropes et anisotropes de. fissures ont été considérées. Ce modèle permet * d'accéder aux modules élastiques effectifs d'une roche sèche ou . d'une roche saturée à l'intérieur de laquelle le fluide est non relaxé. Aussi les vitesses calculées correspondent-elles à des fréquences élevées, de l'ordre de celles mises en jeu en laboratoire (MHz). Une seconde étape de cette étude permet de prédire la dispersion entre les vitesses à haute (MHz) et basse (Hz-kHz) fréquences. La dispersion de vitesse est attribuée à un écoulement local du fluide contrôlé par la ,variation de compressibilité des inclusions poreuses. Dans cette optique, le modèle différentiel auto-cohérent et la limite basse fréquence de la théorie de Biot ont été combinés. Nos résultats soulignent qu'une dispersion importante est attendue pour les rêches fissurées, plus particulièrement dans le cas des ondes de cisaillement. Lorsque deux fluides coexistent dans l'espace poreux et que la saturation n'e st pas uniforme, la modélisation des vitesses nécessite de prendre en compte une seconde famille assimilées à des poches correspondant à des zones de saturations différentes. significatives d'une échelle beaucoup plus grande que celle des pores/fissures et so second mécanisme de dispersion appelé mécanisme de l'écoulement de poche. élaboré qui permet de décrire les vitesses en fonction de la saturation en te hétérogénéités tant à l'échelle microscopique qu'à l'échelle des poches. Nous avons cherché à modéliser le comportement des vitesses mesurées en laboratoire sur des échanti llons pendant des cycles d'imbibition/drainage. L'hystérésis de vitesse observé expérimentalement est simulé en considérant des distributions de fluide différentes suivant le processus de saturation impliqué. Les vitesses prédites s'accordent avec les données expérimentales. A l'aide d'une démarche analogue, on a montré que les anomalies de vitesses détectées avant certains séismes peuvent aussi être significatives d'une saturation non uniforme. Un dernier volet de ce travail, motivé par les recherches menées sur le stockage des déchets nucléaires, a été consacré à l'étude de la fissuration thermique dans des roches ignées. Deux processus de fissuration méritent d'être examinés. Le premier dépend de ta dilatation thermique des minéraux et a été envisagé d'un point de vue expérimental. Le second est lié à la saturation en fluide de la roche. Un modèle a été réalisé afin d'analyser, pour des roches fissurées et saturées, l'effet de la fissuration thermique sur la perméabilité. Cette dernière s'avère contrôlée par deux effets antagonistes: la fermeture des fissures et l'amélioration de la connectivité. L'augmentation de perméabilité attendue pour une élévation de température atteignant 300°C reste modérée (1 à 2 ordres de grandeur).

Nous utilisons des simulations par dynamique moléculaire hors d'équilibre pour étudier le vieillissement des verres de silice après une trempe rapide. Le présent modèle utilise la version modifiée du potentiel de van Beest et l'algorithme de Verlet pour résoudre les équations du mouvement. Ce potentiel reproduit relativement bien la structure et les propriétés dynamiques de la silice. Nous avons effectué des simulations dans un système comportant jusqu'à 9000 atomes pour ne pas avoir des effets liés à la taille de la boite. Dans les conditions que nous utilisons ici, nous avons montré que ce modèle reproduit bien les propriétés statiques et dynamiques de la silice amorphe. Nous avons trouvé que des hétérogénéités dynamiques apparaissent dans la silice surfondue lorsque la température est inférieure à 4000 K. Ces hétérogénéités dynamiques augmentent pendant le processus de vieillissement. Les atomes d'oxygène et de silicium montrent des similarités du point de vue des hétérogénéités dynamiques mais ces dernières évoluent différemment en fonction du temps. Nous avons étudié de manière plus approfondie le comportement des hétérogénéités dynamiques dans la silice surfondue et leur vieillissement. Nous avons aussi observé dans cette gamme de températures, l'apparition des mouvements en chaînes pour les hétérogénéités dynamiques. Nous avons étudié le vieillissement de différentes grandeurs physiques après une trempe rapide. Ces chaînes se comportent d'une manière différente pour les atomes d'oxygène et les atomes de silicium. Nous observons aussi différents types d'agrégation d'atomes. En étudiant le vieillissement de ces différents agrégats, nous trouvons un vieillissement différent dans les diverses agrégations dynamiques.

Les déformations et ruptures gravitaires sont des processus d'évolution communs à tout massif rocheux, et se produisant à différentes échelles. Il est largement admis que ces processus sont influencés par différents facteurs parmi lesquels les hétérogénéités structurales préexistantes (hérités de l'histoire géologique du massif) et les processus d'altération physico-chimique sont de premier ordre. Cependant, l'influence respective de ces paramètres reste floue. Pour avancer dans ce domaine, nous avons utilisé une approche couplée de modélisation physique et numérique en deux et trois dimensions, en nous basant sur deux exemples naturels bien documentés que sont le glissement de La Clapière (Alpes-Maritimes, France) et les écroulements de Randa (canton du Valais, Suisse). Le premier objectif de cette thèse a été l'amélioration de la technique de modélisation physique existante afin de pouvoir d'une part introduire un grand nombre de discontinuités (fractures et failles) dans des modèles physiques, et d'autre part afin de produire des modèles en trois dimensions présentant une topographie réaliste issue de modèles numériques de terrain. Un protocole expérimental de modélisation numérique a ensuite été développé afin d'étudier l'effet de l'altération physico-chimique sur la perte de stabilité des versants rocheux et ce, à grande échelle. Les résultats ont montré, au travers de l'application de la méthode couplée de modélisation physique et numérique (en deux et trois dimensions) que les failles et fractures affectant un massif rocheux (et héritées de son histoire géologique) sont des structures influençant fortement la cinématique des déstabilisations gravitaires, notamment au travers des paramètres propres de ces structures (géométrie, persistance et densité). Il a également été montré que la topographie tridimensionnelle (de premier et de second ordre) ainsi que les propriétés mécaniques à grande échelle (incluant l'altération) ont une grande influence sur le volume mobilisé lors des déstabilisations de versant.

La présente étude porte sur les interactions entre deux types de matériaux susceptibles d'être mis en jeu au sein d'un stockage géologique de déchets radioactifs : les matériaux métalliques, constituant notamment le surconteneur de déchets de haute activité, et les matériaux argileux, telle que la roche hôte argileuse. Les interactions entre ces deux matériaux en conditions de stockage pourraient en effet modifier leurs propriétés de confinement initiales. Les travaux présentés ont visé à déterminer l'influence d'hétérogénéités (vides technologiques et fractures) et d'activités bactériennes sur ces interactions, notamment en termes d'évolution de propriétés chimiques et hydrauliques de l'argile. Dans cet objectif, deux expériences intégrées en cellules de percolation ont été mises en œuvre, à 60 °C, pendant 13 mois : la première en présence de deux souches bactériennes (BSR, BFR), la seconde constituant un témoin abiotique. Ces expérimentations ont permis la circulation d'une eau synthétique dont la composition est représentative de celle de l'eau porale de Tournemire au travers d'un compact de fer pulvérulent puis d'une carotte d'argilite du Toarcien de Tournemire, artificiellement fissurée. L'une des deux demi-carottes d'argile contenait également un cylindre de fer massif. Les caractérisations post-mortem (MEB, MEB/EDS, Raman, DRX, tomographie aux rayons X) ont permis d'étudier deux interfaces : l'interface fer pulvérulent/argilite et l'interface fer massif/argilite.Concernant le déroulement des expérimentations, il convient de noter en premier lieu que la fissure a probablement été circulante durant la totalité de l'essai, ce que tendent à confirmer les modélisations couplées chimie/transport avec le code HYTEC. Toutefois, aucune phase néoformée n'a été mise en évidence à son niveau. Par ailleurs, la survie bactérienne dans la cellule biotique a été confirmée au cours de l'expérimentation par le suivi régulier de la population et par une analyse de la diversité génétique à la fin des essais. Une nette diminution de la concentration en sulfates en sortie de cellule confirme l'activité des BSR.Les caractérisations solides en fin d'essais ont mis en évidence dans la zone de fer pulvérulent, avec et sans bactéries, de la magnétite et de la chukanovite, cette dernière étant localisée majoritairement à l'interface avec l'argilite. Une zone enrichie en fer (10 µm) a été identifié à l'interface dans l'argilite. La vitesse moyenne de corrosion du fer pulvérulent a été estimée à 0,2 µm/an (valeur basse). A l'interface fer massif/argilite, deux faciès ont été mis en évidence. Le premier faciès, identifié dans les deux dispositifs, est composé de chukanovite en couche interne et de sidérite en couche externe. L'extension de la perturbation de l'argilite est de l'ordre de 30 µm. Le second faciès, mis en évidence seulement en présence de bactéries, montre la présence de sulfure de fer (mackinawite) et d'avancées locales de corrosion. Enfin, les simulations HYTEC réalisées ont permis de mieux comprendre les mécanismes biogéochimiques observés expérimentalement, notamment l'effet du pH, et de mieux quantifier certains paramètres cinétiques clefs.

Le calcul de la distribution de dose en Radiothérapie externe se fait en routine clinique à l'aide de Systèmes de Planification de Traitement (TPS) commerciaux. Les algorithmes de calcul de ces TPS ont énormément progressé ces dernières années. Cependant ils sont basés sur des approximations qui restent acceptables pour la plupart des conditions cliniques mais qui montrent leurs limites dans certains cas notamment avec des petites tailles de champ d'irradiation et/ou des faibles densités massiques dans un milieu. Or ces deux conditions sont pourtant réunies dans le cadre de la radiothérapie stéréotaxique des tumeurs bronchiques. Si quelques études ont été réalisées pour des densités massiques classiques de poumon, aucune n'a été réalisée pour des densités pulmonaires très faibles comme par exemple lorsque le patient est traité en inspiration profonde (" Deep Inspiration BreathHold ", i.e. DIBH). Mes travaux de recherche de thèse proposent une étude du calcul de dose pour différentes densités massiques et différentes tailles de champ en se basant sur un modèle Monte-Carlo (MC). La première étape modélise un accélérateur de type TrueBeam(r) (Varian, Palo Alto, CA) en utilisant les données du constructeur. Le modèle est construit à l'aide de la plateforme GATE basée sur la librairie Geant4. Les éléments principaux de la tête de l'appareil sont modélisés. Les espaces de phases (fichiers de particules) fournis par le constructeur au format " .IAEAphsp " sont situés en amont des mâchoires. Pour valider ce modèle, une série de champs simples (3x3 à 20x20 cm2) dans un fantôme d'eau sont implémentés pour des faisceaux de photons de 6X FF (" Flattening Filter "), 6X FFF, 10X FF et 10X FFF (" Flattening Filter Free "). Les résultats (profils, rendements de dose) sont comparés à des mesures de référence obtenues dans une cuve d'eau : respectivement 99% et 97% des points de dose des rendements et des profils respectent les critères de gamma index de 2%-2mm. Une fois le modèle validé, nous avons réalisé une série de simulations pour des champs de petites tailles (3x3 à 8x8 cm2) avec des fantômes hétérogènes de formes simples, pour lesquels la mesure reste accessible. Pour cette dernière, ont été insérés des films radio-chromiques dans des fantômes composés de plaques de PMMA et de deux types de liège de densité 0,12 et 0,24 correspondant respectivement aux poumons en DIBH et en respiration libre. Les résultats du modèle MC pour les quatre énergies ont été confrontés aux mesures expérimentales et aux algorithmes AAA et Acuros (Varian). De façon générale l'algorithme AAA surestime la dose au sein de l'hétérogénéité pulmonaire pour les petites tailles de champ et les faibles densités massiques. Par exemple, pour un champ de 3x3 cm2 et une densité de 0,12 au sein de l'hétérogénéité, une surestimation de la dose absorbée dans le poumon de 16% est mise en évidence pour l'algorithme AAA. Enfin, le modèle est utilisé pour trois cas non mesurables : un objet-test numérique cylindrique hétérogène, des données tomodensitométriques d'un patient en DIBH pour un champ fixe et en arc-thérapie en condition de stéréotaxie pulmonaire. Les résultats ont démontré respectivement pour les études sur TDM une surestimation de la dose dans la tumeur de 7% et 5,4% et dans le poumon de 14% et 9,6% par AAA. D'un point de vue clinique, cela se traduit par un sous-dosage du patient et donc un risque de récidive. / For clinical routine in external Radiotherapy, dose computation is achieved using commercial Treatment Planning Systems (TPS). Since ten years, TPS algorithms have been improved. However they include approximations that are acceptable in most of the clinical cases but they show their limits in some particular conditions for example in presence of small fields and/or low mass y media. And these two conditions are found in the context of stereotactic body radiation therapy (SBRT) of lung tumor. Some studies were published for standard lung densities but none for very low y like in lung during Deep Inspiration Breath Hold (DIBH). This work is a study of dose computation based on a Monte Carlo (MC) model, for different field sizes and mass densities. The first step was to model a TrueBeam(r) linac (Varian, Palo Alto, CA) using data furnished by the manufacturer. This model is built using the Geant4-based GATE platform. The main compounds of the linac head are modeled. Space phase files (i.e. particles files) are furnished by Varian in "IAEAphsp" format and are integrated to the model above the main jaws. To validate this model, a set of simple fields (from 3x3 to 20x20 cm2) in a water phantom is implemented for different photon energies: 6FF, 6FFF, 10FF and 10FFF (FFF = "Flattening Filter Free"). Percentage depth dose (PDD) and lateral profiles are compared to reference measurement in a water tank: respectively 99% and 97% of all the points of these curves passed the Gamma Index test (2% 2mm). Once this validation was completed, a set of simulation was achieved with small field sizes (3x3 to 8x8 cm2) for simple heterogeneous phantoms for which the measurement was achievable. For this purpose, radiochromic films were inserted in phantoms made of PMMA slabs and two types of cork. Cork densities were 0.12 and 0.24 that correspond respectively to lungs during DIBH and free breathing. Results of the MC model for four energies are compared to experimental measurements and to AAA and Acuros Varian's algorithms. AAA algorithm overestimates the dose inside the lung heterogeneity for small field sizes and low density. As an example in the case a 3x3 cm2 field, inside the heterogeneity of density 0.12 an over estimation of 16% in the lung is observed for AAA. The model is finally used for three non-measurable cases: a cylindrical digital reference object and computerized tomography data of a patient during DIBH with a static and stereotactic arc field. Results showed respectively for CT studies an overestimation of dose in the tumor of 7% and 5.4% and in the lungs of 14% and 9.6% by AAA. From a clinical point of view, this means under-dosing the patient and thus a risk of recurrence.