Même si le transport par pipelines de bruts conventionnels reste relativement aisé, il devient plus problématique si le brut est composé d'un grand nombre de molécules hydrocarbonées comme les paraffines. Ces bruts paraffiniques ont la particularité d'avoir une température de prise en gel élevée, souvent supérieure à la température de l'environnement dans lequel il s'écoule. Lorsque l'écoulement dans la conduite est arrêté, la température diminue et entraîne la formation de bouchons de gel de plusieurs centaines de mètres. Après un long arrêt, il devient difficile de redémarrer l'écoulement. Le manque de compréhension des phénomènes impliqués lors du transport d'un brut paraffinique conduit à une surestimation des pressions de redémarrage et donc à un surcoût pour les compagnies pétrolières. Ainsi, ces travaux de thèse s'attachent à améliorer la compréhension du comportement d'un brut paraffinique et à enrichir sa modélisation afin que les simulations numériques prévoient des pressions et des temps de redémarrage plus précis. En faisant référence aux travaux théoriques, expérimentaux et numériques, il est décrit le comportement viscoplastique, thixotrope et compressible d'un brut paraffinique. Les propriétés thermiques et mécaniques liées à l'histoire du brut jouent un rôle essentiel sur le comportement du brut au redémarrage. Ainsi, le transport d'un brut paraffinique avant l'arrêt de la conduite correspond à l'écoulement stationnaire non-isotherme dans une conduite 2D axisymétrique d'un fluide viscoplastique dont les propriétés rhéologiques dépendent de la température. Le comportement d'un brut paraffinique au redémarrage est décrit par le déplacement d'un fluide viscoplastique, thixotrope et compressible par un autre fluide. Le caractère thixotrope est pris en compte par le modèle de Houska et la compressibilité est introduite dans le modèle de redémarrage. La loi de comportement est implémentée grâce à une technique de Lagrangien augmenté associée à un algorithme d'Uzawa afin d'appréhender sa non-dérivabilité. Les équations du problème sont discrétisées grâce à une méthode de Volumes Finis et un schéma TVD (Total Variation Diminishing) est utilisé pour traiter les termes de transport. La méthode de résolution habituellement utilisée pour les fluides viscoplastiques incompressibles est adaptée à des situations compressibles. Différents tests numériques sont menés et montrent les bonnes propriétés de convergence des algorithmes de résolution. Les résultats non-isothermes montrent la sensibilité des zones cisaillées/non-cisaillées aux variations de la température. L'influence des effets combinés de la thixotropie et de la compressibilité sur la capacité de redémarrage est également observée. Un fluide thixotrope incompressible ne pouvant être redémarré, peut être remis en écoulement si une compressibilité lui est ajoutée. Enfin, une étude comparative entre les résultats numériques et expérimentaux permet de valider le code numérique.

Le travail de la thèse porte sur la modélisation fine des lois de comportement thermomécaniques et métallurgiques des aciers afin de déterminer l'état mécanique résiduel des structures lors d'une opération de soudage. Il s'agit de définir les cinétiques de transformation de phases à l'état solide ainsi que le comportement mécanique du matériau multiphasé. Dans un cadre théorique original-les matériaux standards généralisés avec des liaisons internes-nous avons proposé des modèles de comportement pour un matériau multiphasé pour lesquels le comportement mécanique propre à chaque phase est supposé connu. Le comportement macroscopique du mélange est ensuite obtenu via l'hypothèse de Reuss (modèle en série). Soulignons que chaque phase peut présenter un comportement spécifique (plastique ou viscoplastique) et l'effet de la plasticité de transformation (TRIP¬-TRansformation Induced Plasticity) peut être prise en compte, ce qui constitue l'originalité de la thèse. Les modèles ont été implantés dans COMSOL-Multiphysics et deux simulations ont été effectuées dans le but de reproduire des essais représentatifs des phénomènes rencontrés lors d'une opération de soudage ou de traitement thermique des aciers. Les résultats fournis par nos modèles ont été comparés aux résultats expérimentaux issus de la littérature. On a pu ainsi montré que les évolutions thermiques, métallurgiques et mécaniques corrèlent de manière satisfaisante avec l'expérience. Cela confirme la performance et la souplesse de notre approche dans la simulation numérique des structures métalliques sous chargements thermomécaniques anisothermes avec la présence de transformations microstructurales.

La demande en termes de miniaturisation des objets manufacturés est sans cesse croissante. La fabrication en grande série des pièces submillimétriques est rendue difficile par un certain nombre de phénomènes liés à l'échelle. L'objectif de cette thèse est de développer une approche complète permettant de prédire la réponse mécanique d'une tôle d'épaisseur submillimétrique lors d'une opération de mise en forme. Dans un premier temps, des essais expérimentaux en traction uniaxiale sur des tôles minces de cuivre, d'aluminium et d'acier inoxydable sont réalisés. Ils permettent, entre autres, de mettre en évidence l'influence de la surface libre et de la distribution de taille de grains sur la réponse mécanique des tôles. Dans un deuxième temps, une loi de comportement élasto-viscoplastique monocristalline, formulée en grandes déformations et prenant en compte les mécanismes physiques de déformation plastique, est développée et implantée dans le code éléments fins ABAQUS®. Cette loi de comportement utilise les densités de dislocations comme variables internes. Pour être fonctionnelle, cette approche nécessite la description de la microstructure du matériau simulé. Un traitement statistique de mesures réelles est développé afin de générer un modèle éléments finis représentant de façon précise la morphologie des grains, leurs tailles et leurs orientations cristallographiques. Dans un troisième temps, la mise en oeuvre de l'approche complète développée permet de reproduire les phénomènes que nous observons expérimentalement, à savoir : · un adoucissement du comportement global du matériau, · l'apparition d'un état de surface appelé peau d'orange, · une forte localisation de la déformation dans l'épaisseur de l'éprouvette, orientée perpendiculairement à la direction de traction. De plus, grâce à ce nouvel outil numérique, nous vérifions l'influence de la surface libre et la prédominance des gros grains sur la réponse globale des structures calculées. Finalement, une machine d'essais in situ est conçue et réalisée dans le but de valider l'approche numérique, dans des conditions proches de celles observées en micro-formage.

L'objectif de cette étude est de caractériser et de modéliser le comportement mécanique et la rupture en situation accidentelle d'injection de réactivité de gaines de crayons combustibles en Zircaloy-4 détendu vierges, hydrurées ou irradiées en réacteurs nucléaires à eau pressurisée. Un modèle est proposé pour décrire le comportement viscoplastique anisotrope du matériau en fonction de la température (de 20°C à 1100°C), la vitesse de déformation (de 3.10-4 s-1 à 5 s-1), la fluence (de 0 à 1026 n.m-2) et des conditions d'irradiation. Le comportement plastique anisotrope et la rupture du matériau non irradié hydruré jusqu'à 1200 ppm est étudié à l'aide d'essais de traction axiale, traction circonférentielle, expansion due à la compression et traction plane circonférentielle réalisés à 25°C, 350°C et 480°C. La résistance mécanique et l'écrouissage du<br />matériau dépendent de la température et des teneurs en hydrogène en solution solide et en hydrures précipités. L'anisotropie plastique du matériau est peu modifiée par l'hydrogène. A température ambiante le matériau est fragilisé par les hydrures, qui rompent pour des déformations plastiques d'autant plus faibles que la teneur en hydrogène est élevée. La ductilité du matériau, croissante en fonction de la température, n'est pas affectée par l'hydrogène à 350°C et 480°C. Les modes de rupture macroscopiques et les mécanismes d'endommagement diffèrent selon la géométrie des éprouvettes, la température et la teneur en hydrogène. Un modèle de type Gurson est finalement proposé pour représenter le comportement viscoplastique anisotrope et la rupture ductile du matériau en fonction de la température et de la teneur en hydrogène.

L'objectif principal de cette étude est de proposer des lois de comportement macroscopiques permettant de décrire la viscosité de matériaux métalliques sur une large gamme de température et de vitesse de déformation. En effet, dans un contexte industriel fortement concurrentiel, le dimensionnement de structures comme le collecteur d'échappement des moteurs Diesel doit être rapide et prédictif. Les lois de comportement des matériaux jouent un rôle capital dans ce processus de dimensionnement. En effet, la multiaxialité et le caractère anisotherme des chargements appliqués aux collecteurs d'échappement induisent un comportement élastoviscoplastique du matériau et in fine de la structure. Bien qu'elles répondent aux principales exigences de coût, de masse et d'usinabilité, les fontes à graphite sphéroïdal, très souvent choisies comme matériau constitutif de ces pièces, ont un comportement complexe marqué par une très forte viscosité à hautes températures. Ainsi, les structures comme le collecteur d'échappement sont soumises principalement à deux contraintes lors de leur dimensionnement : - La prévention de tout risque d'amorçage de fissure lors de leur durée de vie nominale. On doit ainsi développer des méthodes de dimensionnement à la fatigue thermomécanique oligocyclique car ce sont les cycles de démarrage-arrêt du moteur qui sont endommageant pour la structure. - La prévention de tout risque de fuite et de dérives dimensionnelles. En effet, la succession de cycles thermomécaniques sévères peut induire une accumulation de déformation inélastique dans certaines zones de la structure et alors engendrer de fortes déformées globales. Ces "dérives" peuvent alors être à l'origine de défaut d'étanchéité et de perte de gaz d'échappement. Si une méthodologie efficace pour dimensionner les structures à l'initiation de fissure existait d'ores et déjà dans les bureaux d'études de PSA Peugeot Citroën, la prise en compte des phénomènes de dérives dimensionnelles restait un sujet peu exploré mais néanmoins crucial. En effet, les structures fonctionnent de plus en plus souvent à la limite de leurs capacités et une modélisation efficace de leur comportement permettrait des calculs de dimensionnement prédictifs tout en diminuant les temps de conception. De nombreuses solutions sont envisagées dans la littérature afin de modéliser le comportement élastoviscoplastique de matériau comme la fonte à graphite sphéroïdal. Leur capacité à bien décrire l'évolution d'un matériau soumis à des chargements complexes n'est que rarement avérée pour une large gamme de vitesses de déformation mécanique et inélastique. Par ailleurs, elle dépend fortement de la base expérimentale grâce à laquelle les paramètres de la loi ont été identifiés. Nous nous sommes donc attachés, au cours de cette étude, à constituer des bases expérimentales originales, composées d'essais (écrouissage cyclique, déformation progressive à vitesse constante ou Slow-Fast) qui ont permis de mieux comprendre le comportement de la fonte pour des températures variées (de l'ambiante à 700°C) et des niveaux de déformation multiples. En sollicitant le matériau sur de larges gammes de vitesses de déformation notamment inélastique, ils ont contribué à constituer une base de données souvent plus riche que celles que l'on trouve dans la littérature pour un nombre d'essais et donc un temps de développement réduit. Deux modélisations ont alors été proposées pour décrire le comportement élastoviscoplastique d'une fonte à graphite sphéroïdal particulière : la fonte SiMo. La première d'entre elles se base sur un potentiel visqueux innovant décrivant l'évolution de l'écoulement à l'échelle macroscopique en s'inspirant de la physique à l'échelle des dislocations et de la plasticité cristalline. En combinant différents termes valant pour différents régimes de températures, elle a permis de bien décrire la majorité des essais d'écrouissage cyclique et notamment les cinétiques de relaxation. Constituée de 7 paramètres par température identifiée (soit un de plus que les lois utilisées jusqu'alors par PSA Peugeot Citroën), dont l'évolution avec la température est physiquement ! réaliste, elle ne permet toutefois pas de bien décrire les essais réalisés à vitesse lente ou pour de faibles niveaux de déformation. La seconde loi a donc été développée à partir du même potentiel visqueux enrichi pour réduire la viscoplasticité pour les faibles déformations. Couplée à un terme d'écrouissage isotrope, cette nouvelle loi d'écoulement a permis de corriger les principaux défauts de la première modélisation proposée ainsi que de bien décrire l'ensemble des essais de la base expérimentale. Composée de 11 paramètres, obtenus grâce à une stratégie d'identification adaptée et rapide, elle permet une très bonne description des phénomènes visqueux et est capable de reproduire des essais anisothermes représentatifs des sollicitations sur collecteur avec bien plus d'acuité que les lois de référence utilisées par PSA Peugeot Citroën. Les deux lois de comportement, dont l'admissibilité thermodynamique a été vérifiée, ont été intégrées dans un schéma de calcul de structure. Leur implémentation numérique est fondée sur un algorithme implicite de type retour radial permettant d'obtenir des temps de calcul compatibles avec les exigences des bureaux d'études. Ces modélisations ont ensuite été validées sur deux architectures de collecteur d'échappement et la faisabilité de calculs rapides de collecteur d'échappement sous chargement thermique transitoire a alors été démontrée. Différents calculs de type éléments finis ont permis de démontrer l'amélioration certaine apportée par les lois proposées par rapport aux lois de référence pour déterminer la déformée globale de la structure et pour prévenir d'éventuelles fuites de gaz d'échappement. Ces mêmes calculs ont également permis de montrer que l'utilisation de ces deux nouvelles lois était compatible avec le critère de dimensionnement à la fatigue oligocyclique utilisé chez PSA. On peut donc maintenant prendre efficacement en compte, grâce aux modélisations proposées, les deux contraintes de dimensionnement des collecteurs d'échappement dès les calculs en bureaux d'études et ainsi diminuer les temps et les coûts de développement de telles structures.

Les matériaux polymères sont largement utilisés pour des applications structurelles dans le secteur automobile et leurs comportements complexes nécessitent des modèles précis pour la simulation éléments finis. Les polymères possèdent un comportement dépendant du temps et de la vitesse. La dépendance à la vitesse peut être observée par un accroissement de la rigidité et de la limite élastique en fonction de la vitesse de déformation. Le long temps nécessaire pour retrouver des contraintes nulles après sollicitation du matériau met en évidence la dépendance du temps sur le comportement. De plus, particulièrement pour les polymères chargés, le phénomène de cavitation se traduisant par la création et la croissance de micro-cavités et de microfissures conduit à un changement de volume durant la déformation. Dans ce travail, un modèle de comportement est développé pour un polymère semi-cristallin chargé de talc utilisé dans l'industrie automobile. Un modèle constitutif viscoélastique-viscoplastique non-associatif avec endommagement non-local est proposé dans le but de simuler les phénomènes observés expérimentalement. Dans le modèle développé, une surface de charge non symétrique est utilisée pour prendre en compte la pression hydrostatique. La viscoplasticité non-associative couplée avec l'endommagement conduit aux déformations viscoplastiques non-isochoriques caractérisées expérimentalement. Les paramètres du modèle proviennent d'essais expérimentaux réalisés sous différentes conditions et 'a différentes vitesses de déformation. Pour ces essais, plusieurs techniques de mesure, telles que la corrélation d'images et l'extensommetrie optique sont utilisées pour les mesures de champs de déplacements. La bonne corrélation entre les données expérimentales et les simulations numériques mettent en évidence la précision du modèle développé afin de modéliser le comportement des matériaux polymères semi-cristallins.

Pour la réalisation de pièces de structure automobiles (poutre d'absorption de chocs, quartavants, doublures d'aile...), les matériaux composites à matrice polypropylène renforcée de fibres de verre, continues ou coupées, sont pressentis. L'objet de cette étude est de proposer un modèle de comportement adapté à ces matériaux, capable de répondre de manière cohérente aux sollicitations statiques mais aussi dynamiques subies par ces pièces de structure. Deux matériaux différents ont été étudiés: un tissu et une matrice renforcée de fibres coupées.<br />Afin d'identifier les phénomènes dissipatifs au sein de ces matériaux et de mieux comprendre leur évolution avec la vitesse de sollicitation, une étude expérimentale a été menée incluant des observations microscopiques et des essais mécaniques à différentes vitesses de déformation comprises entre 10-4 s-1 et 100 s-1. Les essais réalisés à vitesse quasi-statique (<10-1 s-1) sont des essais de traction, de compression, de charge-décharge et de traction-compression. L'étude en dynamique (≥10-1 s-1) consiste en la réalisation d'essais de traction. Une étude approfondie de cet essai complexe a par ailleurs été réalisée afin d'améliorer l'analyse des résultats.<br />Un modèle de comportement unifié, valable pour des sollicitations statiques et dynamiques, est ensuite proposé. Il prend en compte les phénomènes d'endommagement, de viscoélasticité et de viscoplasticité observés lors de l'étude expérimentale. Son écriture très générale permet de modéliser le comportement de nombreuses familles de matériaux composites. Il a été implémenté dans le code de calcul ZéBuLoN et identifié pour les deux matériaux étudiés.<br />Enfin, plusieurs essais de validation sont présentés, sur structures simples (flexion à différentes vitesses, choc multiaxial) mais aussi sur pièces industrielles (choc sur poutre d'absorption...).

Le découpage de tôles minces en alliages de cuivre est un procédé de mise en forme très utilisé pour la fabrication de composants électroniques. Deux aspects permettent la caractérisation du procédé : d'une part, le profil de découpe représente la qualité du produit et d'autre part, l'effort de découpage conduit à la définition de l'outillage (matériaux du poinçon et de la matrice, la vitesse et la puissance des presses, etc.). La modélisation numérique permet d'obtenir des réponses rapides sur les questions concernant ces deux aspects. Au sein de l'équipe modélisation et mise en forme des matériaux du Laboratoire de Mécanique Appliquée R. Chaléat, un code de calcul par éléments finis dédié à la simulation bidimensionnelle du découpage a été développé depuis quelques années (Blankform©).<br />Afin d'obtenir une meilleure prédiction de l'effort de découpage par simulation numérique, notre démarche commence par la caractérisation expérimentale du comportement thermo-élasto-viscoplastique de la nuance de cuivre CuNiP par des essais de traction à températures et à vitesses de déformation imposées. Suite à ces essais, une loi de comportement particulière permettant la prise en compte de la sensibilité du matériau à la vitesse et à la température a été proposée. Après l'implémentation de cette loi dans le code de calcul, l'application du modèle au découpage d'une tôle mince en CuNiP, dans un cadre adiabatique, montre que l'effort maximal de découpage est correctement prédit.<br />Lors d'une opération de découpage, certaines nuances de cuivre, comme le Cua1, présentent un phénomène de recristallisation dynamique. Les essais de traction sur cette nuance nous ont permis de caractériser le comportement mécanique du matériau. Sur cette base spécifique, un modèle de comportement thermo-élasto-viscoplastique avec prise en compte de la recristallisation dynamique a été proposé dans le cadre de la thermodynamique des processus irréversibles. L'application de ce modèle au découpage montre qu'il permet de prédire qualitativement les phénomènes observés expérimentalement pour cette nuance.<br />La modélisation bidimensionnelle du découpage est souvent utilisée pour des empreintes de formes relativement simples. Cependant, pour des formes complexes, il est nécessaire de résoudre le problème par une analyse tridimensionnelle. A partir d'un cas d'étude industriel, une modélisation tridimensionnelle du découpage a été réalisée afin d'étudier l'influence des paramètres technologiques sur la distribution des contraintes et des déformations. Les résultats numériques obtenus corroborent bien les constatations expérimentales en termes de profil de découpe.

Tissue characterization is a major step in tissue mechanobiological studies. By characterization methods, tissue quality i.e. the combination of tissue structural, compositional and mechanical properties, is determined. This research focuses on mechanical characterization methods. Among all mechanical characterization methods, we propose those ones which are: 1) Non-destructive, (i.e. that reserves the capability of doing other characterization tests at the end of mechanical test; and, 2) In-line, (that enables tissue progression observation during experiment, and without transferring the specimen from one apparatus to another). However, in-line characterization raises the question of whether conducting tissue observation methods during experimentation modifies tissue progression over time. Therefore, the purpose of this study was to deepen our knowledge about the parameters which could affect tissue quality during mechanical testing. This requires a better understanding of viscoelasticity and viscoplasticity, two key behaviors of tissue, affecting the impact of these parameters (e.g. tissue quality, stimulation parameters) on the response of live tissue to biophysical stimuli. Thus, the objectives of this study were: 1. To review the literature to find information about two mechanical behaviors of tissue i.e. viscoelasticity and viscoplasticity, and the way they affect tissue properties 2. To investigate whether diagnostic tests, as mechanical characterization tests to observe tissue properties, affect tissue progression We explain that viscoelasticity and viscoplasticity of tissue originate from structure and components of the extracellular matrix. We also describe the way they affect tissue dynamic competition between repair, enzymatic degradation and mechanical degradation of the extracellular matrix. Moreover, we specify some tissue stimulation parameters, such as stimulation control type or stimulus history, which could affect tissue progression in response to biophysical stimuli because of viscoelasticity and viscoplasticity. Moreover, by conducting a series of 3-day experiments on frshly extracted tendons, we investigated whether applying "stress relaxation" tests at physiological amplitudes affects tissue response. We divided the tendons into two groups based on the characterization protocol (24 and 0 stress relaxation tests each day), and compared the progression of these groups over time. The stress relaxation tests at physiological amplitude modified tissue response to mechanical stimuli in vitro . In general, the modulus increased for 0 stress relaxation tests, while it first decreased and then increased slightly for 24 stress relaxation tests each day. The difference of mechanical properties between the two groups was significant. Therefore, applying stress relaxation tests at physiological amplitude during the rest periods between mechanical stimuli can affect live tissue progression over time. Therefore, it is essential to take into account the viscoelasticity and viscoplasticity of tissue while developing a stimulation protocol for bioreactor studies or clinical applications.

Les différents travaux de recherche menés dans le cadre de la présente thèse, ont pour objectifs de :(i) comprendre les mécanismes et phénomènes prenant part dans la plasticité des aciers en présence d'une transformation de phases diffusive ou martensitique. (ii) développer des outils de prédiction de TRIP capables d'une part de reproduire correctement cette déformation macroscopique pour des cas de chargements complexes et d'autre part de nous offrir des enseignements sur les interactions élasto-visco-plastiques locales entre phase naissante et phase parente. Pour ce faire, de nouvelles campagnes d'essais TRIP à chargements uni- et bi-axiaux constants ou variables ont été menées sur l'acier 35NCD16 pour une transformation austénitemartensite et l'acier 100C6 pour une transformation austénite-perlite. En complément, des essais de traction-compression et relaxation ont été réalisés pour caractériser les propriétés élastoviscoplastique ssur chacune des 2 phases de l'acier 100C6...