L'approche D.N.L.R. (Distribution of Non Linear Relaxation) repose sur la théorie des fluctuations et l'analyse modale de la dissipation. Elle permet la description d'un grand nombre de comportements physiques.<br /> L'objectif de ce travail consiste à valider le modèle mécanique sur des chargements cycliques<br /> complexes et à intégrer les équations constitutives dans un code de calcul par éléments finis (Msc-Marc). Les schémas numériques habituels utilisant les surfaces de charge ne peuvent être utilisés dans ce cadre. La formulation thermodynamique sous-jacente traite la dissipation à l'aide de réorganisation interne, dont la cinétique est fortement non linéaire. Nous avons mis en place dans le schéma d'intégration des procédures numériques qui accélèrent le calcul. Il s'agit de la technique de pas variable qui s'avère performante lorsque qu'on la compare au modèle N.L.K. (Non Linear Kinematics), à priori plus favorable sur ce plan.<br />L'utilisation de la méthode à grand incrément de temps (MAGIT) prévue pour le calcul des cycles, n'a pas permis de gain de temps par rapport au schéma avec accélération du pas d'intégration.

Les études d'accidents montrent que l'abdomen est devenu le segment prioritaire à protéger pour un occupant ceinturé lors d'un choc automobile frontal, les lésions étant imputables au phénomène de sous marinage. Par ailleurs, un mauvais positionnement de la ceinture de sécurité peut également être à l'origine de lésions lors du déclenchement des prétensionneurs. Ces deux phénomènes entrainent un chargement dynamique de l'abdomen par la ceinture, les vitesses pouvant aller de 4m/s à 8m/s. Le but de cette étude est d'accroitre notre compréhension du comportement de l'abdomen et des mécanismes lésionnels afin d'élaborer un modèle prédictif du risque lésionnel utilisable pour le développement et l'amélioration des systèmes de protection de l'occupant. Pour étudier la réponse de l'abdomen face à un chargement dynamique par une ceinture huit essais sur huit Sujets Humains Post Mortem (SHPM) ont été réalisés. Quatre essais simulaient le phénomène de sous marinage et quatre celui de la prétension de la ceinture. Des bilans d'efforts complets, des mesures de déplacement, de déformées externes de la paroi abdominale, des radiographies et des mesures de pressions ont été réalisés. Les essais disponibles dans la littérature et les résultats expérimentaux ont conduit à développer un modèle masse-ressort-amortisseur reproduisant les principes pilotant le comportement de l'abdomen. Ce modèle souligne le caractère viscoélastique, non linéaire et inertiel de l'abdomen. L'exploitation des mesures réalisées associées aux bilans lésionnels a permis d'apporter des méthodes de validation des modèles numériques par éléments finis et d'améliorer la formulation du modèle être humain HUMOS2

Cette étude vise à modéliser les mécanismes d'endommagement sous chargements multiaxiaux en Fatigue à Grand Nombre de Cycles (FGNC). Le rôle des hétérogénéités microstructurales introduites en particulier par le procédé d'Usinage Assisté Laser fait l'objet d'une attention particulière. On propose de tenir compte de l'aspect dispersé des données en fatigue à l'aide d'une approche de type plan critique et d'un modèle à deux échelles. Plus précisément, une approche probabiliste de type Weibull est d'abord introduite pour prendre en compte le rôle des hétérogénéités microstructurales sur le comportement en fatigue à l'échelle macroscopique. Une deuxième partie de l'étude concerne l'utilisation d'un modèle de couplage plasticité - endommagement à l'échelle mésoscopique afin de refléter les effets de cumul d'endommagement observés en FGNC. Les prédictions des deux approches sont confrontées à des données expérimentales relatives à des chargements complexes (traction biaxiale, traction-torsion hors phase, asynchrone ou encore chargement d'amplitude variable). On montre un très bon comportement d'ensemble ainsi que la possibilité d'associer les deux modélisations pour obtenir une loi d'endommagement probabiliste dont une application directe est l'obtention de courbes P-S-N (Probabilité-Contrainte-Nombre de cycles) pour n'importe quel type de chargement. La dernière partie du travail est consacrée à l'étude de l'interaction entre un procédé d'usinage par assistance Laser (UAL) et la tenue en fatigue des composants métalliques obtenus (acier et alliage de titane). Le rôle des caractéristiques d'intégrité de surface (la microstructure, l'écrouissage, les contraintes résiduelles, l'état de surface

La première étape d'un calcul de durée de vie en fatigue oligocyclique consiste à obtenir l'historique des contraintes et des déformations sur des points jugés critiques. Pour éviter des calculs complets par éléments finis (force brute) particulièrement longs, les corrections élasto-plastiques locales évaluent ces historiques en un temps très court. Les approches existantes manquant de précision, une nouvelle méthode de calcul qui s'inspire des modèles d'homogénéisation a été proposée. Elle a été implémentée dans le code de calcul ZéBuLoN et validée sur des cas de charge quelconques en multiaxial non proportionnel aléatoire, ce qui n'est le cas d'aucune autre méthode de la littérature. Ces historiques servent d'entrée aux méthodes de prévision de durée de vie. Un nouvel algorithme de comptage de cycles a été utilisé pour extraire une série de cycles d'un chargement multiaxial aléatoire. Il s'agit d'une technique de rainow qui permet de conserver toutes les composantes du chargement, et de définir un "cycle" dans l'espace des contraintes déviatoriques. Des calculs de durée de vie ont été réalisés pour évaluer la pertinence de l'ensemble de la chaîne de calcul ainsi constituée. On montre pour finir des applications à une pièce industrielle (bras de châssis de grue).

Les aciers inoxydables austénitiques de type 304L et 316L sont largement employés dans le circuit primaire des centrales nucléaires à Réacteurs à Eau sous Pression (REP). Le retour d'expérience indique la présence d'un nombre limité de fissures intergranulaires dues à la corrosion sous contrainte (CSC) sur des composants en acier inoxydable écroui. Il a été démontré qu'une pré-déformation importante associée à un chargement cyclique favorise la propagation des fissures de CSC. L'objectif de l'étude est d'améliorer la compréhension du rôle de la pré-déformation par traction ou par laminage sur les mécanismes d'amorçage de la CSC pour les aciers inoxydables austénitiques. Le comportement mécanique des matériaux écrouis a été caractérisé et des essais d'amorçage en milieu primaire simulé ont été réalisés sur des éprouvettes entaillées. L'ensemble des essais d'amorçage réalisés a confirmé un fort effet du trajet de chargement sur la sensibilité à l'amorçage des matériaux étudiés, quel que soit le niveau de pré-déformation. Un critère global a été proposé pour réunir les deux aspects de l'amorçage de la fissuration que sont la densité de fissures et leur profondeur. Ce critère est utile pour caractériser l'amorçage tandis que la profondeur maximale de fissure est le paramètre pertinent pour définir la transition entre propagation lente et propagation rapide. Des cartes de sensibilité à l'amorçage vrai ont été établies. Une profondeur critique de fissure de 10 à 20 µm a été déterminée pour les aciers 316L A et B pré-déformés par laminage. Elle est comprise entre 20 µm et 50 µm pour les matériaux pré-déformés par traction. Une ébauche de modèle d'ingénierie applicable aux aciers inoxydables austénitiques a été proposée : l'effet de la température est négligeable dans la gamme 290°C-360°C et l'impact de la contrainte sur le temps pour obtenir la transition varie comme (max/Rp0.2, T°C)11,5. L'effet du trajet de chargement ainsi que de l'écrouissage de surface dû à l'usinage sont intégrés indirectement à l'indice de contrainte, à ce stade du modèle. L'effet " matériau " observé dans cette étude tient principalement à l'effet du trajet de déformation. La puissance élevée de la dépendance à la contrainte traduit l'intégration de différents paramètres favorisant la localisation de la déformation. Pour cette raison, il sera nécessaire de définir le champ des contraintes locales pour parvenir à une modélisation plus physique.

De nos jours, le transport de marchandise occupe une place importante dans la vie économique des sociétés modernes. Le problème de collecte et de distribution avec fenêtres de temps à un seul véhicule est un des problèmes les plus rencontrés. Ayant un ensemble de demandes à satisfaire, le véhicule doit transporter des biens de fournisseurs à leurs clients respectifs en respectant les fenêtres de temps et sa capacité. Dans ce travail, nous présentons un état de l'art du 1-PDPTW et nous proposons plusieurs approches évolutionnistes pour traiter ses deux cas : statique et dynamique. Nos approches utilisent principalement des algorithmes évolutionnistes basés sur l'utilisation d'opérateurs génétiques spéciaux conçus dans le but d'améliorer la qualité des solutions et de diminuer le temps de calcul. Elles sont aussi basées sur la Pareto optimalité pour fournir un ensemble de solutions viables. Quelques unes de nos approches utilisent des bornes inférieures de distance et de retard dans le but d'évaluer les résultats obtenus. Une recherche Tabou, constituant un étage d'hybridation, peut être appliquée pour l'amélioration des solutions obtenues par les algorithmes évolutionnistes. Enfin nous présentons quelques simulations et résultats élaborés à partir de benchmarks spécialement conçus pour le 1-PDPTW ainsi que d'autres provenant de la littérature.

Cette thèse s'ajoute aux différents travaux de recherche qui traitent des alliages d'aluminium fortement utilisés dans l'industrie aéronautique et contribue fortement à comprendre le comportement élastoplastique en chargement cyclique à contrainte imposée du 2017A. L'apport essentiel de ce travail est l'étude de l'anisotropie propre du matériau utilisé à travers le suivi de l'évolution des différents paramètres caractérisant la plasticité cyclique de notre matériau. En effet, nous avons caractérisé cette anisotropie en comparant le comportement du matériau en traction-compression avec celui de la torsion alternée selon l'évolution cyclique de la réponse contrainte-déformation, l'évolution de l'état stabilisé, l'évolution des variables d'écrouissages cinématique et isotrope ainsi que l'anisotropie selon le comportement en fatigue et endommagement. Pour mieux affiner la partie expérimentale de ce travail, des investigations microstructurales des faciès de rupture de toutes les éprouvettes utilisées ont été effectuées afin de mieux comprendre les mécanismes d'endommagement cyclique dans notre matériau. Dans la partie numérique de cette thèse, nous avons réalisé des simulations numériques en utilisant la dernière version du modèle multimécanismes qui tient compte de l'anisotropie du matériau. Les résultats de ces simulations, réalisées en considérant les mêmes conditions de nos essais expérimentaux, confirment les capacités de cette nouvelle version à estimer le comportement élastoplastique d'un matériau anisotrope.

Ce travail de thèse vise à une meilleure compréhension et prise en compte des mécanismes d'endommagement ductile au cours des procédés de mise en forme à froid. Le tréfilage, le laminage de fils et le laminage à pas de pèlerin sont pris comme exemples d'application. Une grande attention est portée également à la méthodologie d'identification des paramètres des modèles d'endommagement implémentés. Trois approches de l'endommagement ductile ont été étudiées : des critères de rupture phénoménologique, des modèles phénoménologiques couplés et des modèles micromécaniques. Ces modèles ont été implémentés dans Forge®, et les algorithmes ont dû être adaptés à sa formulation mixte vitesse-pression et à son élément fini (P1+/P1). Parallèlement aux développements numériques, différents essais mécaniques ont été effectués sur trois matériaux différents (acier haut carbone, acier inoxydable et alliage de zirconium) pour identifier les paramètres des modèles d'écrouissage et d'endommagement. Des essais de traction in situ sous micro-tomographie aux rayons X ont été exploités pour l'identification des mécanismes d'endommagement ductile (germination, croissance et coalescence), ainsi que l'identification du modèle micromécanique. Enfin, la dernière partie de ce travail a consisté à réaliser des études comparatives de ces différents modèles d'endommagement pour les différents procédés de mise en forme et différents matériaux. Concernant le tréfilage et le laminage de l'acier inoxydable, un bon accord entre les simulations numériques et les résultats expérimentaux a été obtenu. Pour l'acier à haute teneur en carbone (perlitique), le modèle micro-mécanique GTN a donné le meilleur résultat à la fois qualitativement et quantitativement pour le tréfilage ultime. Par ailleurs, la comparaison entre les différents modèles sur différents procédés (laminage du fil avec l'acier haut carbone, laminage à pas de pèlerin avec l'alliage de zirconium) a aidé d'une part à mettre en évidence le rôle important du troisième invariant de la contrainte déviatorique sur la localisation de l'endommagement dans les procédés de mise en forme où le cisaillement est dominant. D'autre part, elle montre que le processus d'identification lui-même devrait se fonder sur la mesure de la microstructure afin de garantir un résultat précis lors de l'application aux procédés.

Devant le défi grandissant de diminuer la libération dans l'atmosphère des gaz à effet de serre plusieurs solutions ont été envisagées. Une des plus étudiée est le stockage géologique du CO2 par injection dans des aquifères salins profonds. Ces structures géologiques sont présentes dans les bassins sédimentaires et se caractérisent par leur épaisseur (plusieurs centaines de mètres) et leur étendue géographique (l'ensemble du bassin parisien, par exemple, pour la couche du Dogger). Le second interêt de ces structures est que les roches carbonatées qui les composent sont généralement très poreuses (10 à 30% de porosité) et présentent une perméabilité suffisamment élevée pour faciliter l'injection de fluides. Cependant, le comportement et l'évolution à long terme de ces réservoirs, après l'injection de CO2 dans la saumure en place, est incertain. Aussi, la résistance mécanique et les propriétés d'écoulement de ces roches sont étudiées pour réaliser des modèles prédictifs d'évolution, afin de garantir l'intégrité et la pérennité des sites de stockage. Le travail de cette thèse repose sur un protocole d'acidification retardée, développé et utilisé à l'IFPEN, permettant une altération homogène d'échantillons de carbonates. Notre objectif a été de déterminer les effets de ce type d'altération sur le comportement hydromécanique de la roche calcaire de Lavoux. La caractérisation des propriétés d'écoulement de la roche a été combinée à une étude microstructurale. Les paramètres étudiés sont la porosité et la perméabilité avant et après altération. Les résultats d'évolution de la perméabilité ont été analysés en s'appuyant sur une méthode de soustraction d'images de microtomographie RX, qui permet de visualiser la mobilité de particules fines au sein du réseau poreux. Ces dernières sont aussi révélées à postériori dans les effluents par diffraction laser. L'étude microstructurale de la roche après altération a permis la mise en évidence des mécanismes de dissolution et d'altération. Le comportement mécanique de notre roche à différents stades d'altération a été étudié à l'aide d'essais triaxiaux classiques réalisés à différents états de confinement. Ces essais ont montré que l'altération induit une transition du domaine de rupture fragile vers le domaine de rupture ductile. Afin de mieux comprendre ce comportement, nous avons mené des essais mécaniques uniaxiaux avec observation et acquisition d'images numériques " in-situ " : i) sur une presse électromécanique classique équipée de caméras CCD et de microscopie optique, ii) à l'aide d'un dispositif de compression miniaturisé adapté à l'enceinte d'un microscope électronique à balayage. La technique de corrélation d'images numériques a pu être appliquée tout au long du chargement mécanique à différentes échelles, depuis celle de l'échantillon, jusqu'à celle de la microstructure. Les mesures de champs de déformation ainsi réalisées ont permis de mettre en évidence les modes d'endommagement et de rupture selon l'état d'altération. Une cellule triaxiale miniaturisée pour être montée dans l'enceinte du microtomographe a spécialement été développée à l'IFPEN pour cette étude. Cette cellule nous a permis de mener des essais à 5 et 10 MPa de confinement avec observation volumique " in-situ " par tomographie. Nous avons pu identifier pour trois états d'altération, trois types de comportements à la rupture, avec l'observation i) d'une bande de cisaillement en régime de rupture fragile, ii) d'une étroite bande de cisaillement en compaction, iii) d'une large bande de compaction en cisaillement. L'ensemble des résultats obtenus sur les propriétés mécaniques et sur les propriétés d'écoulement (notamment la variation de perméabilité et la présence de fines dans le réseau poreux) tend à montrer que l'altération affecte fortement la microporosité du calcaire, avec des conséquences quant au comportement cisaillant compactant de notre roche.

Le système de failles de Haiyuan qui borde le plateau du Tibet au Nord-Est est un système majeur sénestre. Au cours du dernier siècle, deux grands séismes (M~8) ont rompu ce système de failles: le séisme de Haiyuan en 1920 et le séisme de Gulang en 1927. A l'aide d'interférométrie radar à synthèse d'ouverture, nous analysons les variations spatiales et temporelles de la déformation intersismique au travers de la faille de Haiyuan, dans une zone étendue (150x150 km2) qui couvre l'extrémité Ouest de la rupture de 1920 et la lacune sismique de Tianzhu. Avec une approche dite en Small Baseline, nous traitons cinq séries temporelles d'images SAR, acquises par le satellite Envisat le long de tracks descendantes et ascendantes pendant la période allant de 2003 à 2009. Les cartes de vitesse moyenne de déformation dans la ligne de visée du satellite ainsi obtenues sont cohérentes avec un mouvement sénestre au travers de la faille et montrent des variations latérales du gradient de vitesse dans la zone de faille. Nous inversons ces cartes de vitesse moyenne en LOS pour obtenir le taux de chargement à court terme en profondeur et la distribution du glissement dans la partie sismogène le long du plan de faille. Le taux de chargement en profondeur est d'environ 5mm/an. Les sections de faille ayant rompu en 1920 et une grande partie de la lacune sismique de Tianzhu sont bloquées en surface. Entre ces deux sections, un segment de 35 km de long, qui montre une forte activité micro-sismique, glisse de manière asismique avec un taux de glissement horizontal qui atteint presque 5 mm/an. Cependant, le taux de glissement asismique le long de la partie sismogène varie le long du plan de faille et atteint localement des taux supérieurs au chargement tectonique, suggérant des variations temporelles du glissement asismique. La comparaison de profils moyens de vitesse parallèle à la faille issus de données InSAR sur les périodes 1993-1998 (données ERS) et 2003-2009 suggèrent une migration vers la surface du glissement asismique sur une période de 20~ans. Une analyse en séries temporelles des données Envisat, en appliquant un lissage temporel, montrent une accélération du taux de glissement asismique pendant l'année 2007. Cette accélération est précédée et a probablement été déclenchée par un séisme de magnitude 4.7 au sein même du glissement asismique. Enfin, nous étudions la relation entre l'évolution spatio-temporelle du glissement asismique en surface et la rugosité de la trace de la faille à l'aide d'une analyse multi-échelle. Nous montrons que les propriétés élastiques de la croûte cassante contrôlent la rugosité de la faille, qui exerce à son tour un contrôle sur la distribution de glissement asismique en surface. Le glissement asismique est fait de spasmes qui interagissent les uns avec les autres en suivant une loi d'échelle similaire à la loi de Gutenberg-Richter pour les séismes.