Étude du comportement post-coulée de réfractaires électrofondus à Très Haute Teneure en Zircone (THTZ)

Étude du comportement post-coulée de réfractaires électrofondus à Très Haute Teneure en Zircone (THTZ)

Modèles numériques de plasticité de transformation

Dans ce paragraphe, nous allons nous intéresser aux travaux réalisés sur la plasticité de transformation induite par la transformation martensitique des aciers. Dans ses travaux de thèse, N. T. TRINH (TRI 08) a fait un travail de synthèse des modèles numériques de plasticité de transformation dans les aciers.

Mécanismes

Le phénomène de plasticité de transformation est généralement expliqué par deux mécanismes : le mécanisme de Greenwood-Johnson (GRE 65) et le mécanisme de Magee (MAG 66). • Mécanisme de Greenwood-Johnson. Lors d’une transformation de phase, les deux phases en présence (phase mère et phase fille) n’ont pas la même compacité. Cette différence de compacité entre les deux phases engendre, dans la phase ayant la limite d’élasticité la plus faible, des déformations plastiques à l’échelle microscopique. Sans influence d’une contrainte extérieure, cette déformation plastique va s’orienter de façon aléatoire et ne laisser apparaître, à l’échelle macroscopique, que la déformation isotrope induite par le changement de volume. En revanche, si une contrainte extérieure est appliquée au cours de la transformation de phase, la déformation plastique sera alors  orientée, entraînant ainsi à l’échelle macroscopique une déformation plastique irréversible. • Mécanisme de Magee. Le mécanisme de Magee est, quant à lui, décrit dans le cas de la transformation martensitique. Les plaquettes formées au cours de la transformation vont s’orienter sous l’effet d’une contrainte externe. De ce fait, la variation de volume qui résulte de la transformation martensitique ne sera pas isotrope et entraînera une déformation irréversible à l’échelle macroscopique. La prépondérance de l’un ou l’autre des deux mécanismes dépend du type de transformation considéré. Il est d’usage de prendre en compte le mécanisme de GreenwoodJohnson dans le cas de transformation par diffusion, et de prendre en compte celui de Magee dans le cas de la transformation martensitique. Dans le cas des matériaux THTZ, la transformation de phase s’apparentant à une transformation martensitique, on devrait donc plutôt considérer l’effet Magee. La modélisation des effets mécaniques des transformations martensitiques est encore un sujet ouvert. La majorité des travaux réalisés sur des aciers et sur la transformation martensitique ne prennent en compte que l’effet Greenwood-Johnson et négligent l’effet Magee. Par exemple, ne considérer que l’effet Greenwood-Johnson peut conduire à une modélisation satisfaisante d’un acier 16MND5 (TAL 01). La figure 1-29 représente la comparaison entre les résultats expérimentaux et un modèle de simulation numérique de type Leblond (LEB 89I, LEB 89II). Figure 1-29 : plasticité de transformation en fonction de la température, comparaison entre simulation numérique (modèle où l’effet Magee est négligé) et résultats expérimentaux (TAL 01) Dans le cas de la transformation martensitique, ces travaux démontrent que la description de la plasticité de transformation par le seul effet Greenwood-Johnson peut être réaliste, dans le cas où la contrainte appliquée est dite faible, c’est-à-dire très inférieure à la limite d’élasticité de la phase la plus faible en présence pendant la transformation. Nous allons donc nous intéresser par la suite plus spécialement à ce type de modèles. 

Modèle de Leblond

Dans la littérature, la déformation totale ε d’un matériau est communément décomposée en plusieurs termes : cp cp tp T e d ε ε ε ε ε ε σ = + + + + 1.45 Avec ε e la déformation élastique, ε d la déformation due à la dilatation thermique, εT cp et εσ cp les déformations dues à la plasticité pure dépendante de la température (T) et de la contrainte (σ), ε tp le terme de déformation plastique de transformation. Le fait de se placer dans le cas de contraintes faibles a pour conséquence de minimiser les termes de déformation dus à la plasticité pure devant celui dû à la plasticité de transformation. Dans les modèles prenant en compte uniquement l’effet Greenwood-Johnson, les termes εT cp T et εσ cp σ sont négligés et la déformation totale s’écrit : e d tp ε = ε + ε + ε 1.46 Le modèle de Leblond (LEB 89I, LEB 89II) définit la vitesse de plasticité de transformation comme une fonction dépendante de la proportion de phase transformée (ϕ(z) ) et de la contrainte appliquée (S étant la partie déviatorique du tenseur des contraintes) : K z z S tp 2 3 * (‘ ) * * • • ε = ϕ 1.47 K est alors un paramètre de plasticité de transformation à définir, constant, homogène à l’inverse d’une contrainte. La fonction ) ϕ(z est normalisée en z avec φ(0)=0 et φ(1)=1. Différentes études conduisent à des valeurs de K et ) ϕ(z différentes. Dans le cas du modèle de Leblond (LEB 89), lorsque les contraintes appliquées sont faibles, l’auteur propose : ϕ(z) = z 1( − ln z) et 1 12 2 y th K σ ∆ε = avec th 12 ∆ε la différence de compacité entre les deux phases. J. C. Videau (VID 96) a étudié expérimentalement la variation du taux de transformation et a proposé deux fonctions qu’il a comparées aux résultats expérimentaux : φ(z) = z(1-ln(z)) et φ(z) = z(2-z). Les résultats sont présentés sur la figure 1- 30 pour trois charges appliquées différentes :  Figure 1-30 : comparaison simulation – expérience de l’évolution de phase transformée en fonction de la contrainte appliquée (VID 96) Une des particularités du modèle de Leblond est de faire intervenir cette plasticité de transformation uniquement au moment de la transformation de phase. Ce modèle a ensuite été amélioré (LEB 89) par la prise en compte de nouvelles hypothèses. Des écrouissages isotropes, puis cinématiques, peuvent par exemple être pris en compte. Le mécanisme de Magee peut aussi être rajouté au modèle de Leblond, comme dans le cas du modèle de Fischer (FIS 00). 1.2.3.3 Modèle de Videau Videau et al. (VID 94) se sont aussi intéressés à la modélisation mécanique multiaxiale de la plasticité induite par la transformation. Dans cette étude, les auteurs décrivent une formulation viscoplastique. Ils admettent une additivité du mode de déformation par plasticité de transformation par le mécanisme de Greenwood-Johnson et de la déformation viscoplastique « classique ». La déformation inélastique du matériau s’écrit alors : ine vp pt ε = ε + ε 1.48 On définit alors deux domaines d’élasticité avec deux fonctions de charges indépendantes, une pour la viscoplasticité classique et l’autre pour la plasticité de CHAPITRE 1 : Etude bibliographique 48 transformation. Un terme d’écrouissage cinématique X et un terme d’écrouissage isotrope R sont ajoutés au modèle. Les fonctions de charge du domaine de viscoplasticité classique ( vp f ) et de plasticité de transformation pt f sont définies de la façon suivante : pt pt pt vp vp vp d vp f X R f X R = − − = − ( ) _ φ σ σ 1.49 Classiquement les lois d’écoulement vont s’écrire à l’aide des lois de normalité : pt pt vp vp p n n • • • • = = ε ε ν 1.50 avec vp n et pt n les directions d’écoulement. σ σ ∂ ∂ = ∂ ∂ = pt pt vp vp f n f n 1.51 Les équations du modèle s’écrivent : pt pt n vp vp p K f K f = = • • ν 1.52 Les paramètres n, Kvp et Kpt sont à déterminer par méthode inverse, d’après des essais menés sur le matériau de l’étude. En ce qui concerne l’écrouissage, le couplage des variables cinématiques et isotropes mène à la définition d’une matrice d’écrouissage :                 =         pt vp pt vp pt vp C C C C X X α α et                         =         pt vp pt vp vp pt vp r r q q R q q R R 0 0 1.53 Avec les évolutions des variables d’écrouissage : ( ) ( ) pt pt pt pt vp vp vp vp p n d X n d X = − = − • • • • α α ν et 1( ) 1( ( )) 0 pt pt pt vp vp vp vp r p b R r b R R = − = − − • • • • ν 1.54 où C, Cvp, Cpt,dvp,dpt,bvp et bpt sont des caractéristiques du matériau. CHAPITRE 1 : Etude bibliographique 49 Ce modèle a été implémenté dans le code de calcul ZéBuLon. La figure 1-31 représente les résultats fournis par le calcul pour la simulation d’une « trempe à induction » d’une éprouvette entaillée axisymétrique en acier : Figure 1-31 : Evolution de la température et des contraintes résiduelles dans la section médiane (VID 94)

Résumé

Dans ce chapitre nous avons mené une étude bibliographique de la simulation numérique des procédés de fonderie. Nous avons vu que la simulation numérique est découplée en deux parties : -Une première partie thermique : celle-ci s’appuie sur un essai où le moule de coulée est instrumenté par des thermocouples qui permettent de suivre l’évolution de la température au cours du refroidissement en différents points de la structure de coulée. Grâce à ce suivi des températures, on peut définir un paramètre d’échange thermique de chaleur à l’interface matériau/moule, qui permet de construire le modèle numérique thermique et de reconstruire l’évolution de la température au sein du matériau, inaccessible dans notre cas en raison de la très haute température de coulée. -Une deuxième partie mécanique : celle-ci part des résultats thermiques qui en définissent le chargement. Il convient pour l’étude mécanique de définir les lois de comportement qui régissent le comportement mécanique du matériau en fonction de la température. Dans notre cas, nous avons mis l’accent sur la transformation de phase de la zircone et sur le phénomène de plasticité de transformation induit. Dans son étude sur les matériaux THTZ, L. Massard (MAS 05) a mis en évidence la déformation importante d’une éprouvette soumise à une contrainte de traction ou de compression, dans l’intervalle de température où la zircone subit son changement de phase de tétragonale à monoclinique. Par analogie avec le cas de la transformation martensitique des aciers, ce comportement est typique d’une plasticité de transformation : écoulement plastique anormal provoqué par l’application d’une contrainte au cours de sa transformation (VID 94). Ce comportement particulier est pris en compte par les modèles de simulation numérique grâce à l’introduction d’un terme de plasticité de transformation. Nous proposons un modèle simple inspiré du modèle proposé par Leblond (LEB 89I, LEB 89II). 

Table des matières

INTRODUCTION GENERALE
Introduction du chapitre 1
1.1. Modélisation numérique thermique
1.1.1 Formulation mathématique des échanges thermiques
1.1.1.1 Premier et second principe de la thermodynamique
1.1.1.2 Equations de fermeture
1.1.1.3 Formulation faible pour la méthode des éléments finis de type Galerkin
1.1.2 Description du phénomène d’échange thermique à l’interface
1.1.3 Mise en évidence expérimentale des résistances thermiques d’interfaces
1.1.4 Détermination du paramètre h par méthode inverse
1.1.5 De la thermique à la mécanique
1.1.6 Résumé
1.2. Plasticité de transformation (TRIP : Transformation Induced Plasticity)
1.2.1 Transformation de phase de la zircone
1.2.1.1 Transformation de type martensitique
1.2.1.2 Mécanismes de transformation
1.2.2 Influence de la contrainte
1.2.2.1 Cas des THTZ
1.2.2.2 Cas des métaux
1.2.3 Modèles numériques de plasticité de transformation
1.2.3.1 Mécanismes
1.2.3.2 Modèle de Leblond
1.2.3.3 Modèle de Videau
1.3. Résumé
Introduction du chapitre 2
2.1 Fluage : état de l’art
2.1.1 Description générale d’un essai de fluage
2.1.2 Mécanismes de fluage
2.1.3 Dissymétrie de comportement
2.1.4 Récapitulaté
2.2 Caractérisation expérimentale
2.2.1 Démarche expérimentale
2.2.1.1 Mode de sollicitation
2.2.1.2 Hypothèses faites pour l’analyse des résultats
2.2.1.4 Essais réalisés
2.2.2 Dispersion expérimental
2.2.3 Résultats expérimentaux : Matériau ZS
2.2.4 Résultats expérimentaux : Matériau ZB 6
2.2.5 Résultats expérimentaux : Matériau ZBY1
2.2.6 Récapitulatif et comparaison du comportement en fluage des trois matériaux
2.3 Détermination d’une loi de comportement
2.3.1 Décomposition de la déformation
2.3.2 Fonction de charge et écrouissage
2.3.3 Loi de comportement viscoplastique
2.3.4 Méthode d’optimisation
2.3.5 Identification de la loi de comportement
2.3.5.1 Matériau ZS
2.3.5.2 Matériau ZB
2.3.5.3 Matériau ZBY1
2.3.6 Résumé
2.4 Mesure du module d’élasticité en fonction de la température
2.4.1 Module d’Young
2.4.1.1 Méthode de mesure du module d’Young mascroscopique
2.4.1.2 Résultats des mesures
2.4.1.4 Conclusion
2.4.2 Coefficient de dilatation thermique
2.5 Modélisation des déformations irréversibles induites par la transformation de la zircone
2.5.1 Application du modèle de Leblond au cas des matériaux THTZ
2.5.2 Mise en évidence expérimentale
2.5.3 Loi de comportement
2.6 Récapitulatif
Introduction du chapitre
3.1 Instrumentation d’une caisse de coulée et présentation du maillage associé
3.1.1 Description d’une caisse de coulée
3.1.2 Géométrie de la caisse de coulée
3.1.3 Positionnement des thermocouples dans la caisse de coulée
3.1.4 Description du maillage utilisé pour les calculs
3.2 Analyses des données expérimentales et mise en données de la simulation
3.2.1 Analyse des résultats expérimentaux par rapport à la distance de la dalle pour les deux faces
3.2.2 Analyse des résultats expérimentaux dans la hauteur pour chacune des deux faces
3.2.3 Mise en données des calculs thermiques
3.2.3.1 Propriétés thermiques des différents constituants de la caisse de coulée
3.2.3.2 Echanges de chaleur dans la caisse de coulée
3.2.4 Récapitulatif
3.3 Simulation numérique de la partie thermique
3.3.1 Etude de l’échange de chaleur entre la dalle et le moule en graphite au cours du refroidissement
3.3.2 Etude du paramètre de résistance thermique d’interface
3.3.2.1 Dans l’épaisseur de graphite (MM)
3.3.2.2 Dans l’épaisseur d’agent de calage (MA)
3.3.3 Etude des paramètres thermique du sable
3.3.3.1 Méthode d’optimisation et résultats
3.3.3.2 Propriétés thermo physiques du mélange sable alumine fournit par Saint Gobain
3.3.3.3 Comparaison des valeurs optimisées par simulation numérique avec les données fournies par Saint Gobain et la littérature
3.4 Résumé
Introduction du chapitre 4
4.1 Refroidissement d’une dalle de THTZ dans un moule en graphite
4.1.1 Points représentatifs du refroidissement au sein de la dalle
4.1.2 Détermination de la durée maximale de refroidissement à simuler
4.1.3 Analyse des courbes de refroidissement
4.1.3.1 Zone pâteuse
4.1.3.2 Zone de transformation de phase de la zircone
4.1.5 Résumé
4.2 Paramètres influençant la cinétique de refroidissement
4.2.1 Influence de l’épaisseur de graphite
4.2.2 Influence du type d’agent de calage
4.2.3 Résumé
Introduction du chapitre 5
5.1 Etude de l’intervalle de température où le comportement est viscoplastique
5.1.1 Mise en données
5.1.1.1 Variables d’élasticité
5.1.1.2 Lois de comportement viscoplastique
5.1.2 Module d’Young constant en fonction de la température
5.1.2.1 Cas de l’élasticité
5.1.2.2 Cas viscoplastique
5.1.3 Module d’Young variable en fonction de la température
5.1.4 Influence de la loi de viscoplasticité
5.1.5 Résumé
5.2 Etude de la transformation de phase de la zircone
5.2.1 Expansion volumique
5.2.1.1 Mise en données
5.2.1.2 Résultats
5.2.1.3 Conclusion
5.2.2 Déformations irréversibles induites par la transformation de phase de la zircone
5.2.2.1 Première étude : Conséquence de la plasticité de transformation sur les contraintes
5.2.2.2 Comparaison des deux modèles de plasticité identifiés expérimentalement
5.2.2.3 Deuxième étude : influence de la plasticité de transformation sans fluage actif
5.2.3 Résumé
5.3 Etude du comportement à basse température
5.3.1 Mise en données
5.3.2 Résultats
5.4 Application du modèle thermomécanique au cas du moulage en moule céramique
5.4.1 Etude thermique
5.4.1.1 Géométrie de la caisse de coulée
5.4.1.2 Coulée instrumentée en usine
5.4.1.3 Maillage et données thermiques
5.4.1.4 Ajustement du modèle thermique
5.4.1.5 Comparaison des deux procédés de fabrication
5.4.2 Mécanique
5.5 Résumé
Résumé et Perspectives
ANNEXE A
A.1 Méthodes de mesures
A.1.1 Flexion trois points, cyclage par palier (F3Ppalier)
A.1.2 Flexion trois points, cyclage constant (F3Pcst)
A.1.3 Traction-compression alternée (TCalt)
A.1.4 Récapitulatif
A.2 Résultats : Matériau ZS
A.2.1 Méthode F3Ppalier
A.2.2 Méthode F3Pcst et F3Pcst800
A.2.3 Méthode TCalt (GEMH)
A.2.4 Récapitulatif des résultats du matériau ZS et conclusions
A.3 Résultat : Matériau ZB
A.3.1 Méthode F3Ppalier
A.3.2 Méthode F3Pcst et F3Pcst800
A.3.3 Méthode TCalt (GEMH) et ultrasonore
A.3.4 Récapitulatif des résultats du matériau ZB et conclusions
A.4 Comparaison matériau ZS et matériau ZB
ANNEXE B
B-1 Matériau ZS
B.1.1 Essais réalisés à 1400°C
B.1.1.1 Sous 0.5 MPa
B.1.1.2 Sous 1 MPa
B.1.1.3 Sous 2 MPa
B.1.2 Essais réalisés à 1350°C
B.1.3 Essais réalisés à 1300°C
B-2 Matériau ZB
B.2.1 Essais réalisés à 1400°C
B.2.1.1 Sous 0.5 MPa
B.2.1.2 Sous 1 MPa
B.2.1.3 Sous 2 MPa
B.2.2 Essais réalisés à 1350°C
B.2.3 Essais réalisés à 1300°C
B-3 Matériau ZY
B.3.1 Essais réalisés à 1400°C
B.3.2 Essais réalisés à 1350°C et 1300°C
ANNEXE C
C-1 Agent de calage : mélange céramique
C.1.1 Conductivité thermique
C.1.2 Chaleur spécifique
C-2 Sable furanique

projet fin d'etudeTélécharger le document complet

Télécharger aussi :

Laisser un commentaire

Votre adresse e-mail ne sera pas publiée. Les champs obligatoires sont indiqués avec *