Alliage Homogene-Héterog-Solution Solide PDF

REPUBLIQUE ALGERIEN DEMOCRATIQUE ET POPULAIRE
Ministère de l’enseignement supérieur et de la recherche scientifique

Université El-Hadj Lakhdar-BATNA-
MEMOIRE
Présenté au département des sciences de la matière

Faculté des sciences
Pour obtenir le diplôme de

magistère en physique
Option : physique des matériaux
Par
ZEMMOURI SALIMA
THEME
Etude d’une variété de solutions solides binaires

traitées par la méthode des sous réseaux de Hillert.
Soutenue le : 31/05/2012
Devant le jury :
Président : A. BALGACEM-BOUZIDA Pr U. Batna

Rapporteur : Y. DJABALLAH M.C. A U. Batna
Examinateur : F. BOUHARKAT Pr U. Batna
E. BELBACHA Pr U. Batna
S. BENALIA M.C.A U.Tissimsilt
-0-
PDF created with pdfFactory Pro trial version www.pdffactory.com

Dédicace
« Les oeuvres de l’esprit, de la science et du savoir doivent être au service de l’humanité pour
rendre la vie plus simple et plus belle, et non à l’enrichissement de quelques individus »
Nikola Tesla (1856-1943)
N A mes parents
Ñ A mes sœurs
Ñ A mes frères
Y A toute ma famille
-0-

REMERCIMENT
Tous d’abord je remercie le bon DIEU le tout puissant.
Je tiens à remercie très vivement Monsieur Djaballah Yassine maitre de

conférence à l’université de Batna, je lui suis gré de m’avoir orientée, suivie et
dirigée dans mon travail, tout mon profond respect pour son aide ses précieuses
orientations et ses remarques pertinentes.
Je remercie également Monsieur Belgacem Bouzida Aissa, professeur à

l’université de Batna, qui ma fait l’honneur de s’intéresser à ce travail proposé, et
de présider le jury.
Je n’oublie surtout pas monsieur Belbacha El-Djemai, professeur à l’université

de Batna, pour ses encouragements, son intérêt pour mon sujet et de faire
partie du jury.
Ma profonde reconnaissance va ensuite envers, Madame Bouharkat Fouzia,

Professeur à l’université de Batna, qui m’a aidé et qui m'a souvent encouragé et
conseillé avec beaucoup d'efficacité. Je la remercie encore une fois pour sa
disponibilité, sa gentillesse et sa bonne humeur quotidienne. Qu’il trouve ici
l’expression de ma haute gratitude et mon profond respect. Je lui adresse
également mes vifs remerciements pour m’avoir fait l’honneur d’examiner ce
travail.
Mes remerciements vont également à monsieur Salah Benalia, maître de

conférence au centre universitaire de Tissemsilt, d’avoir accepté de juger ce
travail et de faire partie du jury.
Je tiens cependant à remercier tout particulièrement Amel Naceur professeur

à l’université de Biskra pour son aide inestimable. La bonne maîtrise des
diagrammes de phases n’a été possible que grâce à son soutien, ses remarques
et conseils précieux. Je la remercie également de m’avoir redonné du courage à
chaque fois que j’en avais besoin. Bonne continuation Amel.
Enfin un grand merci à ma famille et à tous mes amis qui, de près ou de loin, ont
contribué au bon déroulement de ce mémoire.
-0-

Introduction Générale ………………………………………………… …………………… .1
Chapitre I
ALLIAGE ET DIAGRAMME DE PHASES

I-1-Alliage …………………………………………………………………………………………….3
I-1-1- Alliage homogène……………………………………………………………………………...3
I-1-2- Alliage hétérogène…………………………………………………………………….4
I-2- La phase ………………………………………………………………………………...4
I-2-1- Homogénéité de la composition ……………………………………………………....4
I-2-2- Structure d’une phase……………………………………………………………….....5
I-2-3- Aspect densité ou masse volumique d’une phase…………………………………..…..5
I-3-Les Solutions Solides ………………………………………………………………..…...6
I-3-1-Les Solutions Solides Primaires……………………………………………………..….6
I-3-2-Les Solutions Solides d'Insertion (ou Interstitielles)……………………………..……..7
I-3-3-Les Solutions Solides de Substitution………………………………………………………….7
I-3-4-Les Solutions solides intermédiaires………………………………………………….....8
I-3-4-1- solutions ordonnées et désordonnées……………………………………….……….9
I-3-4-2-composés intermétalliques ordonnés (Composés Définis)……………………...…..10
I-4- Structures et dénominations des phases métalliques……………………………………10
I-4-1-Groupe d’espace…………………………………………………………………….…..11
I-4-2-Notation Pearson…………………………………………………………………….….12
I-5- Les Diagrammes de Phases binaire………………………………………….……….….12
I-5-1-Règle de phase de Gibbs……………………………………………………………...…13
I-5-2- systèmes binaires isomorphes ……………………………………………………….…14
I-5-3-Systèmes binaires eutectiques…………………………………………..……………….15
I-5-4-Systèmes binaires péritectique…………………………………………………………..16
I -5-5-Systèmes binaires présentant des phases intermédiaires…………………………........17
I -5-5-1-Formation de composés intermédiaires définis Am Bn………………………..….....17
I- 5-5-2-Diagramme binaire avec solution solide secondaire …………………………....…19
I -5-6-Systèmes binaires présentant une miscibilité partielle à l’état liquide….......................19
I -5-7- Ensemble des équilibres triphasés………………………………………21
-0-

Chapitre II
LES BASES THERMODYNAMIQUES DES DIAGRAMMES DE PHASES
II-1- Introduction……………………………………………………………………..…….....22
II-2- Premier principe de la thermodynamique ……………………………………………..22
II-3-Deuxième principe de la thermodynamique…………………………………………….23
II-4-Troisième principe de la thermodynamique ……………………………………………23
II-5-Potentielle chimique……………………………………………………………………...24
II-6-L’activité……………………………………………………………………………….…25
II-7-Théorème d’Euler …………………………………………………….…………………26
II-8-Relation Gibbs- Duhem………………………………………………………………….26
II-9- Relation de Gibbs- Helmholtz…………………………………………………………...27
II-10-L’énergie de Gibbs molaire et le potentiel chimique ………………………………….27
II-11-Critères d’équilibres………………………………………………………………….....28
II-12-Stabilité d’une phase…………………………………………………………………....28
II-13- Construction de la tangente commune et le potentiel chimique…………………..….30
II-14- Détermination des courbes de solidus, liquidus et solvus……………………………..30
II-14-1-La miscibilité totale de A et B………………………………………………….……...31
II-14-2- Gap de miscibilité dans l’état solide………………………………………………….31
II-14-3-Le point eutectique sans isomorphes………………………………………………......32
II-14-4-Le point eutectique dans le cas isomorphe…………………………………………….32
II-15- Présentation de la méthode CALPHAD……………………………………….............35
II-16- Description thermodynamique de phases d’un diagramme d’équilibre……………...37
II-16-1-Les éléments purs……………………………………………………………….……..37
II-16-2- phases stœchiométriques (Composés définis)………………………………………...38
II-16-3-phase étendu (Modèle polynomial).…………………………………………………...38
II-16-4-Phases intermédiaires…………………………………………………………………40
-0-

Chapitre III
LE MODELE DES SOUS-RE SEAUX

III-1-Introduction ……………………………………………………………………………42
III-2- modèle à plusieur sous-réseau………………………………………………………...43
III-2-1-Définition de fraction de site………………………………………………………...43
III-2-2- L’entropie idéal de mélange…………………………………………………………45
III-3- Modèle à deux sous-résau……………………………………………………………..46
III-3-1- L’enthalpie libre molaire ……………………………………………………………….. 46
III-3-1-1- L’enthalpie libre de référence de Gibbs …………………………………..……….47
III-3-1-2- L’enthalpie libre idéal de mélange de Gibbs……………………………………...48
III-3-1-3- L’enthalpie libre d’excès de mélange de Gibbs ……………………………………48
III-3-1- 4-La relation entre l’enthalpie atomique et celle de site………………………….….49
III-3-2-Cas particuliers……………………………………………………………………….49
III-3-2-1-Cas d’une solution …………………………………………………………………..49
III-3-2-2-Cas d’un composé défini……………………………………………………….....…49
III-3-2-3- Modèle à deux sous-réseaux identiques pour un système binaire ……………..…..50
III-3-2-4-modèle à deux sous-réseau non- identique pour un système ternaire…...……..……52
III-3-2-5- modèle à deux sous-réseau non- identique pour un système binaire………………..53
III-3-2-6- modèle à deux sous-réseaux présentant du lacune pour un système binaire………...54
III-4-généralisation du modèle CEF…………………………………………………………..55
III-5-La transformation ordre ↔ désordre ……………………………………………….….56
III-6- Modélisation de la réaction ordre-désordre……………………………………………..57
III-7-Modélisation de l’ordre à courte distance……………………………………….............58
III-8-Application du modèle aux différentes phases…………………………………………..59
III-8-1- Modélisation des phases de structure B2……………………………………………….59
III-8-2- Modélisation de la phase σ et χ ……………………………………………………….61
III-8-3- Modélisation de la phase AB2 (phases de laves)………………………………………..62
III-8-4 -modélisation de la phase CFC ( L12, A1, L10)……………………………………..…..65
III-9 –conclusion………………………………………………………………………………..71
-0-

Chapitre IV
ETUDE DES SYSTEME BINAIRE PAR LA METHODE DES SOUS-RESEAUX
IV-1- introduction ………………………………………………………………………………74

IV-2-Présentation du programme BATNABIN ………………………………………………..74
IV -3- Méthode de calcul………………………………………………………………………...75
IV-4-Présentation du programme PANDAT……………………………………...…………….75
IV-5-Etude du système Cobalt-Antimoine…………………………………………………..…..76
IV-5-1-Les données expérimentales ………………………….………………………….............…76
IV-5-2-Resultats et discussion…………………………………………………………………....78
IV-6-Etude du système Cadmium-Antimoine:………………………………………………….85
IV-6-1-Les données experimentales ………………………….…………………………………….85
IV-6-2- Resultats et discussion ………………………………………………………..................84
Conclusion générale………………….………………………………………………………….90
Bibliographie………………………….………………………………………………………….92
Annexe ……………………………….………………………………………………………..…96
-0-

INTRODUCTION
GENERALE
-0-

Chapitre I Alliages et diagrammes de phases
Introduction générale
Les métaux purs sont très rarement utilisés dans les constructions industrielles car ils
présentent souvent des caractéristiques mécaniques très limitées, et la plupart des matériaux
de construction sont constitués de différents éléments. Lorsque ces éléments ne sont que
partiellement miscibles, le matériau est constitué de plusieurs phase de structures cristallines
et de compositions différentes. Il en résulte une influence très importante sur les
caractéristiques mécaniques du matériau et son comportement.
L’étude de la microstructure du matériau consiste à déterminer la nature des
différentes phases en présence, leurs proportions et leurs constitutions. Cette analyse se fait à
partir des diagrammes d’équilibres de phases qui traduisent les états d’équilibre entre les
différentes phases. Le diagramme de phases d’un mélange est illustration graphique des
conditions de température, de pression et de composition aux quelles les différentes phases
de ce mélange sont thermodynamiquement stables. Il se caractérise par : le nombre de phase,
le nombre de constituants et le degré de liberté.
Les diagrammes de phases présentent souvent des phases liquides, solides primaires,
des composés définis et des phases intermédiaires non-stoechiométriques. Pour calculer un
diagramme de phase il est indispensable de modéliser ces différents types de phases par
différents modèles existant dans la littérature. Notamment, les polynômes de Redlich-Kister
sont souvent utilisés pour la modélisation de l’enthalpie d’excès des solutions solides
désordonnées primaire. Mais les phases intermédiaires ordonnées considérées pendant
longtemps comme stoechiométriques dans un but de simplifier leurs traitements sont
aujourd’hui modélisés par la méthode des sous-réseaus régulières de Hillert et Staffanson
(1997) initialement développé pour les solutions salines par Temekin. Son formalisme très
général permet de prendre également en compte les solutions interstitielles et lacunes et
s’applique pratiquement à tous les types de solutions intermédiaires non-stoechiométriques
[1,2]
L’idée générale du modèle est venue de la représentation des solutions de sels pour
les quels les constituants sont soit des anions, soit des cations qui ne se permutent qu’entre
ion de même signes, c’est le modèle de Temekin. Il a été supposé que dans un cristal ionique
ils existent deux sous réseaux. L’entropie configurationnelle à été décrite séparément pour les
-1-

cations et anions. Plus tard le modèle à été utilisé pour les systèmes métalliques par Hillert et
il a été adapté pour les sous réseaux multiples par Sundman et Agren [3,4].
Dans le cadre des systèmes métalliques, l’hypothèse de plusieurs sous réseaux n’a
aucun sens physique. Les termes plus corrects et utilisés actuellement sont Sublattice
Formalism ou Compound Energy Formalism (CEF). Ce dernier a été appliqué à la
modélisation d’une grande variété de phases et des méthodes ont été développées pour traiter
différentes situations. [5]
Notre travail comporte les parties suivantes :
Ø Dans le premier chapitre on présente un rappel bibliographique sur les alliages
et différents types de diagrammes de phases.
Ø Le deuxième chapitre est consacré à l’étude des bases thermodynamique des
diagrammes de phases ainsi que la méthode d’optimisation et d’établissement
des diagrammes de phases connue sous le nom CALPHAD qui se base sur la
minimisation des enthalpies libre. Les différents modèles utilisés dans le
calcul ont été également présentés.
Ø Le troisième chapitre est une étude détaillée sur le modèle des sous-réseaux.
Après une présentation générale de ce modèle avec les cas particuliers, on
donne les formules des enthalpies libres d’excès pour chaque structure
possible des phases intermédiaires.
Ø Dans le dernier chapitre on a optimisé deux systèmes binaires (Co-Sb) (Cd-Sb)
en utilisant deux code de calcul BATNABIN et PANDAT [6].
-2-

I-1-Alliage
Un alliage est une combinaison d'un métal avec un ou plusieurs autres éléments
chimiques.
Un métal pur a des caractéristiques mécaniques relativement faibles. Le fait d'ajouter
d'autres éléments permet de « durcir » (augmenter les caractéristiques mécaniques).
Également, ces ajouts permettent de modifier les caractéristiques chimiques des métaux (en
particulier leur comportement à la corrosion) ou d'améliorer d'autres caractéristiques (facilité
de mise en œuvre : coulabilité par exemple).
Le métal principal, la plus importante partie du mélange, est appelé « métal de base »
ou « base ». Les éléments ajoutés volontairement sont appelés « éléments d'alliage » (ou
d'addition) et les éléments non désirés sont appelés « impuretés ».
Les éléments d'alliages sont le plus souvent des métaux, mais peuvent également être
d'autres éléments chimiques : le carbone dans l'acier ou la fonte, le silicium dans
l'aluminium,… etc.
I-1-1- Alliage homogène

Un alliage homogène peut être ordonné (les atomes de différentes natures suivent une
alternance stricte) ou désordonné (les différents atomes occupent des places aléatoire
(Fig- I-1)).
a b
Figure-I-1 :
(a) alliage désordonné (b) alliage ordonné
I-1-2- Alliage hétérogène

Lorsque la teneur en élément d'alliage augmente, on peut avoir formation de deux
phases : une phase contenant peu d'éléments d'alliage, et une phase à forte teneur en éléments
d'alliage (Fig-I-2).
Les cristallites à forte teneur sont appelées « précipités ». Les précipités sont souvent
des alliages ordonnés, que l'on appelle « intermétalliques ». Les intermétalliques ainsi formés
-3-

sont parfois par la suite étudiés en tant qu'alliages propres, comme un nouveau matériau, et on
essaie d'en produire en tant que tel et non plus en tant que précipités [7].
Figure- I-2 : Alliage hétérogène
I-2- La phase
Pour définir l’idée actuelle attachée au concept de phases, trois critères doivent être
pris en compte pour définir et identifier une phase à l’intérieur d’un système :
I-2-1- Homogénéité de la composition :

Pour mesurer la concentration locale d’une phase, il faut procéder à la l’analyse
chimique d’un certain prélèvement volumique de la substance. Il est certain qu’il faut préciser
la valeur de ce volume de prélèvement si la définition d’homogénéité doit avoir un sens. En
effet un volume de l’ordre de quelque Angstrom au cube isolerait l’atome, ce qui est une
entité trop petite pour apprécier l’homogénéité d’un mélange. Un volume important
permettrait, au contraire, des erreurs d’interprétation dans le cas de phases dispersées liées
entre elles on peut donc, penser que le micron cube est une échelle intéressante de
prélèvement. Pour les phases solides, il existe précisément un instrument d’analyse chimique
qui analyse les diverses substances à l’intérieur du (μm) 3, c’est la microsonde de casting.
Le principe de cet instrument consiste à exciter, par l’impact d’un faisceau
électronique pointu, les atomes de la matière d’une cible (Fig. I-3). Environ 1 μm 3 de matière
excité participe alors à l’émission d’un spectre de rayon X, qui est analysé dans un
spectrographe. Les longueurs d’ondes émises sont caractéristiques de la nature chimique des
atomes et les intensités sont proportionnelles aux fractions massiques (atomiques). Cet
instrument convient donc relativement bien pour vérifier l’homogénéité des concentrations
chimiques dans les phases solides.
-4-

D’autre part, la thermodynamique phénoménologique est toujours globalisante par

rapport à la réalité du matériau. Il existe dans le volume même d’une phase réputé homogène,
des fluctuations locales de concentrations liées à la diffusion de la matière. Ces fluctuations
peuvent être aléatoire ou bien s’organiser en ondulations de longueur d’onde définie, par
exemple, dans le mécanisme de décomposition spinodale. En général, on ne doit pas
considérer une fluctuation locale de concentration comme un changement de phase. De la
même façon un grain fraîchement cristallisé à partir d’un liquide complexe présente toujours
une hétérogénéité de composition entre le cœur et la surface de ces dendrites. On ne doit pas
pour autant le considérer comme système polyphasé, mais comme une seule phase n’ayant
pas encore totalement atteint son état d’équilibre [8].
Figure - I-3 : Principe de l’analyse par microsonde électronique de Castaing.
I-2-2- Structure d’une phase

Le critère d’analyse n’est pas suffisant pour s’assurer de l’unicité d’une phase dans un
volume donné de matière, même très petit. Or une même phase ne peut posséder qu’une
organisation structurale. On peut donc relier la notion de phase à celle de structure cristalline.
La matière peut avoir une structure de type amorphe ou cristallisé. On doit considérer comme
un changement de phase tout passage amorphe → cristal ou tout changement de groupe de
symétrie cristalline [8].
I-2-3- Aspect densité ou masse volumique d’une phase

A l’intérieur d’une même phase, on ne doit pas non plus observer une discontinuité de
densité local. Ce critère de densité pourrait être rattaché à la notion d’homogénéité de
concentration volumique à la place des fractions molaires ou massiques. Soit ρ la masse
volumique local d’un élément ΔV de matériaux.
-5-

A l’intérieur d’une même phase ρ ne doit pas subir aucune discontinuité. Les
fluctuations de matière par diffusion et les fluctuations locales de volume interdisent de poser
de façon absolue ρ= cte. Mais en pratique, on se rapproche de cette condition lorsque la phase
est dans son état d’équilibre [8].
I-3-Les solutions solides

Les solutions solides sont des agrégats, chimiquement homogènes, résultant de
l'addition d'un ou de plusieurs éléments étrangers (éléments d'alliage) dans un métal pur
lorsque ce dernier est présent en faible teneur dans un alliage. Il s’agit de mélanges intimes à
l’échelle atomique dont les atomes du métalloïde ou du métal B peuvent entrer en solution
dans le métal A (souvent appelé matrice) [9].
Dans certains cas, les métaux en présence sont miscibles en toute proportion. C'est
le cas par exemple de l'or (Au) et de l'argent (Ag) qui forment une solution solide continue.
La plupart du temps, il existe une concentration limite (limite de solubilité) en atomes de
soluté au delà de laquelle la structure cristalline est modifiée. Dans l'intervalle de
concentration borné par cette concentration limite, la solution solide est dite primaire ou
terminale (exemple le laiton α1). Au delà de cette concentration, les deux constituants
peuvent former une solution solide ayant une structure cristalline différente de celle du
métal de base : il s'agit d'une solution solide intermédiaire (exemple le laiton β1).
On observe également des solutions dont l'existence n'est possible que dans un
domaine de concentration limité : on dit alors que les constituants forment un composé
défini.
Les atomes étrangers peuvent entrer en solution avec le métal de base soit en se
substituant à ses atomes, on parle alors de solution solide de substitution ; Soit en se plaçant
aux interstices de son réseau, on parle alors de solution solide d'insertion. On peut
distinguer ces deux types de solutions solides par des mesures combinées de densité et de
paramètre cristallin (par diffraction de rayons X) [10].
I-3-1-Les solutions solides primaires

Sont des solutions solides formées à partir du solvant pur ou du soluté pur. Elles
ont évidemment la structure cristallographique du métal pur dont elles dérivent. Si les deux
métaux purs ont même structure cristallographique, la solution peut être continue. Les
-6-

solutions solides primaires sont définies par la limite de solubilité du ou des éléments
d'alliages (B) dans le métal de base (A). Un certain nombre de règles empiriques, appelées
règles de Hume-Rothery, permet d'évaluer qualitativement cette grandeur [10]:
Ø Effet de taille
Ø Effet de l'électronégativité (polarisation)
Ø Effet de valence :
Ø Effet de concentration électronique :
En résumé, il est sans doute illusoire de vouloir trouver des concentrations
magiques qui régiraient la stabilité de diverses phases. Mais il n’en est pas moins vrai que
la concentration électronique est un paramètre utile pour classer les phases métalliques.
I-3-2-Les solutions solides d'insertion (ou interstitielles)

Lorsque le rayon atomique de l'élément B est suffisamment faible, cet élément peut
occuper les sites interstitiels de la structure A. Seuls les atomes de faible diamètre pourront
conduire à des solutions solides de ce type, par exemple C, N, O, H, B [11]. Il faut que le
rayon de l’atome interstitiel ne dépasse pas 59 % de celui de la matrice A. Les solutions
solides interstitielles sont surtout rencontrées dans les métaux de transition (ou de terres rares,
ou d’actinides). Cela est sans doute dû à la facilité de créer une liaison covalente de forte
énergie entre les orbitales de valence d d’un métal de transition et les orbitales s et p des
atomes légers précités. Dans les métaux alcalins, beaucoup plus électropositifs, il y a tout de
suite tendance à former un composé ionique tel que Li+ H–. Les solutions solides
interstitielles ont toujours un caractère métallique marqué.
Les solutions solides interstitielles peuvent présenter des phénomènes d’ordre: par
exemple, le carbone dissous dans le tantale peut s’ordonner pour donner une phase Ta64C
[10].
I-3-3-Les Solutions Solides de Substitution

La plupart des solutions solides sont des solutions solides de substitution : les
atomes étrangers occupent une fraction des sites réticulaires à la place des atomes de base.
La structure cristalline est en général inchangée mais le paramètre de maille varie avec la
concentration en atomes étrangers. Les atomes de base et de l'élément d'alliage peuvent être
répartis complètement au hasard sur les divers sites du réseau et la solution est dite
-7-

désordonnée, c'est le cas quand les deux éléments constituants sont parfaitement
équivalents (même structure cristalline et dimensions atomiques voisines). Dans d'autres
cas, il y a une tendance plus ou moins marquée à l'acquisition d'un ordre (solution
ordonnée) ou à un rassemblement d'atomes du même type comme cela est indiqué sur la
figure-I-4.
Solution désordonnées Solution parfaitement ordonnées Rassemblement
Figure- I-4 : Type de solutions solides de substitution.
I-3-4-Les Solutions solides intermédiaires

Quand les deux type d’atomes métalliques A et B sont de tailles très différentes, ils
ont tendance à former des composés ordonnés dont le domaine de composition est souvent
très étroit.
Le domaine d'existence de ces phases intermédiaires évoquées précédemment est
compris entre deux valeurs limites de la concentration en élément d'addition. Lorsque la phase
intermédiaire n'existe qu'à l'état ordonné, on parle de composé intermétallique ou composé
défini. Ces composés ne sont généralement stables que sur des domaines de concentration
beaucoup plus étroits que les solutions solides intermédiaires désordonnées. Aussi, peut-on les
qualifier de composés stœchiométriques bien que l'écart à la stœchiométrie soit parfois très
important. Suivant les cas, ces composés ont un caractère métallique plus ou moins marqué. Ce
caractère décroît quand on passe des solutions solides intermédiaires aux composés
intermétalliques à caractère plus ou moins covalent (ou ionique).
Approximativement, on peut classer les phases intermédiaires en plusieurs ensembles :
solutions solides ordonnées, désordonnées (ou à faible énergie de mise en ordre), composé
intermétalliques ordonnés.
-8-

I-3-4-1- Solutions ordonnées et désordonnées

Les phases intermédiaires peuvent être désordonnées ou ordonnées suivant q’un site
du réseau cristallin est occupé indifféremment par un type A ou B, ou au contraire
préférentiellement par un atome de ces atomes [11].
Une solution solide est désordonnée si la probabilité de trouver un atome A en
un site donné ne dépend pas de la nature des atomes occupant les sites voisins. Dans ce cas,
les deux éléments en présence doivent être équivalents. Les solutions sont rarement
idéalement désordonnées. Une solution solide est complètement ordonnée si chaque type
d’atome correspond un type de site donné. Dans ce cas, les liaisons entre les atomes des deux
éléments sont énergétiquement favorisées par rapport aux liaisons entre atomes de même
type. A haute température, l'ordre sera détruit par l'agitation thermique qui provoque un
brassage permanent des atomes sur les divers sites. La structure ordonnée disparaît à partir
d'une température critique Tc (température de transformation ordre-désordre). Les
phénomènes d’ordre se rencontrent aussi bien dans les solutions de substitution que
d’insertion [12]. La figure I-5 donne quelques exemples de structures ordonnées.
a b C
Figure- I-5 : Exemple de structures ordonnées
a- maille CFC(L10) -b- maille CFC(L12) -c- maille CC(L2)
• La structure L 10 (type AuCu)

C'est une structure CFC formée par l'alternance de plans (001) d'atomes de (Cu) et
de plans (001) d'atomes de (Au). Le cube perd ainsi certains de ses éléments de symétrie. La
structure ordonnée devient quadratique avec c/a = 0,93 correspondant à un effet stérique de
tassement de la structure.
• La structure L 12 (type AuCu3)
C'est également une structure CFC dans laquelle les atomes d'une espèce se trouvent
aux sommets du cube et les atomes de l'autre espèce se trouvent aux centres des faces du
-9-

cube, ce qui correspond dans le cas de Au et Cu à la composition stœchiométrique AuCu3. La

symétrie de la maille n'est plus cubique à faces centrées mais cubique simple.
• La structure L 2 (type CsCl ou laiton beta)
Le laiton β a une structure CC. Pour la composition équiatomique, la structure
ordonnée présente une alternance de plans de Cu et de Zn du type (001) ; une espèce
d'atomes occupe le centre du cube et l'autre les sommets du cube. La symétrie cubique
centrée est perdue, elle devient cubique simple.
I-3-4-2-Composés intermétalliques ordonnés (Composés définis)

On distingue généralement pour les composés intermétalliques :
Ø Les composés électroniques (phases de Hume-Rothry) tels que : CuZnβ, Cu3Al ou
Cu5Sn.
Ø Les phases de Laves : du type AB2 ordonné ou le rapport des rayons atomiques
Ra /rB est souvent voisin de 1.2 et qui peuvent se présenter sous trois formes: MgZn2
(hexagonal), MgCu2(cubique), CuAl2. ou MgNi2 (hexagonal) [10].

Ø les phases σ qui apparaissent dans de nombreux alliages de métaux de transition par
exemple autour de la composition équiatomique dans le diagramme (Fe-Cr) sont
souvent considérées comme appartenant à cette catégorie. Elles ont pourtant un
domaine d’existence assez étendu. Leurs structures sont complexes [11].
I-4- Structures et dénominations des phases métalliques

Les métallurgistes ont toujours utilisé des dénominations particulières pour désigner
les structures cristallines des phases solides. Pour les métaux purs, la structure peut être
décrite par un empilement de sphères dures au contact les unes des autres et qui représentent
les atomes ionisés. Ces empilements de sphères se divisent en deux grandes familles, les
empilements compacts ou chaque sphère possède 12 premiers voisins et les empilements non
compacts dont le principal représentant possède 8 premiers voisins autour d’une sphère
centrale. Et pour plus de simplicité des structures métalliques, Ewald et Hermann avaient
imaginé une classification appelée « Strukturbericht » qui a était abandonnée parce qu’elle
présentait des contradictions internes et qu’elle n’était pas exhaustive. Mais pratiquement
l’utilisation de ces notations est encore continue, souvent sans en connaître la signification.
C’est pourquoi nous donnons quelques clefs pour en comprendre la logique [13].
- 10 -

Les corps purs et les solutions solides primaires sont notés An, n représente le numéro
dans la liste des structures des métaux purs et comprend 1 ou 2 chiffres. Les deux principaux
empilements compacts sont ainsi A1, structure cubique à faces centrées et A3, structure
hexagonale compacte. La structure non compacte la plus souvent rencontrée est cubique,
centrée et s’appelle A2. Il est remarquable de constater que 90 % des métaux purs cristallisent
dans l’une des structures A1, A2, ou A3.

Les solutions solides intermédiaires et les composés définis sont notés Bn, Cn et DIn,
selon que leur stoechiométrie est respectivement 1/1 comme AlFe, 1/2 comme AuAl2 ou plus
complexe. Le chiffre I dans DIn désigne la stoechiométrie particulière, par exemple D0n
désigne la stoechiométrie 1/3 comme AlFe3 ; dans tous les cas n est un numéro dans la série
considérée et qui comporte 1 ou 2 chiffres.
Les phases ordonnées qui se rattachent à la structure désordonnée An sont aussi
appelées Lnm le numéro n est identique dans les deux structures ordonnée et désordonnée
tandis que l’indice m est le numéro dans la série engendrée. Par exemple, la structure cubique
à faces centrées A1 engendre fréquemment deux variétés ordonnées : L10 comme AuCu ou
AlTi et L12 comme AuCu3 ou AlNi3. C’est ce jeu de double entrée qui provoque l’incohérence
de cette classification, par exemple les phases ordonnées non compactes CuZn ou AlFe se
rattachent à la structure désordonnée A2 et à ce titre s’appellent L20 mais elles possèdent la
stoechiométrie 1/1 et pour cela s’appellent aussi B2 [14].
I-4-1-Groupe d’espace
Il est vain de vouloir construire une classification cristallographique simplifiée pour
les seules solutions solides métalliques car les variétés cristallines des solutions intermédiaires
aujourd’hui connues couvrent un très large éventail des groupes d’espace. On sait que la
cristallographie distingue 7 groupes cristallins, selon la forme extérieure de la maille, 14
modes de Bravais qui permettent par un jeu de maille multiple de ne pas masquer la plus
haute symétrie d’une structure, 32 groupes d’espace appelés groupes finis, basés sur la
combinaison autour d’un point, les éléments de symétrie ne comportant pas de translation et
230 groupes infinis en incluant la translation dans les éléments de symétrie. On a trouvé des
solutions solides intermédiaires appartenant pratiquement à tous ces groupes d’espace. Alors
que la notation simplifié devient impossible [14].
- 11 -

I-4-2-Notation Pearson
Aujourd’hui environ 30 000 cas. Les notations Pearson sont extrêmement simples :
une lettre minuscule désigne la première lettre du système cristallin de la phase considérée,
suivie d’une lettre majuscule qui est le mode de Bravais, suivie enfin du nombre d’atomes par
maille. Ainsi les structures A1, A2 et A3 sont notées respectivement cF4, cI2, et hP2, (F comme
Faces centrées, I comme Intérieur centré, P comme Primitif). Dans les anciens recueils de
diagrammes de phases le texte descriptif des structures des phases solides se réfère
essentiellement aux notations du Strukturbericht et si nécessaire aux groupes d’espace des
cristallographes. Actuellement les trois désignations de structure sont données comme suit :
Pearson, groupe d’espace et Strukturbericht ainsi que la phase prototype. Par ailleurs, Hume-
Rothery a introduit ces phases sur le tracé du diagramme qui sont souvent désignées par des
lettres grecques. Ces phases sont des analogies structurales (phases appelées composés
électroniques) observées dans les alliages de cuivre, d’argent et d’or pour une même valence
moyenne, puis elle s’est étendue à de très nombreux diagrammes ; pour cette notation on
prend souvent les lettres grecques dans leur ordre alphabétique en composition croissante de
l’élément B.
Exemple : dans les laitons (Cu,Zn), on rencontre successivement les solutions solides
a (cF4, Fm m, A1, Cu), β (cI2, Im m, A2, W), β’ (cP2, Pm m, B2, CsCl), γ (cI52, I 3 m, D82 ,
Cu5Zn8) , δ = β, ε (hP2, P63/mmc, A3, Mg) et η = ε.
Les phases β et β’ sont en relation ordre « désordre : l’utilisation de l’apostrophe pour
désigner la variante ordonnée est fréquente. Parfois aussi l’apostrophe désigne une phase
martensitique qui n’est pas une phase d’équilibre thermodynamique. [14].
I-5- Les Diagrammes de phases binaire

Le diagramme d'équilibre d'un système est une représentation graphique des
domaines d'existence, en fonction des variables d'état (c, T) des différents états sous lesquels
ce système existe à l'équilibre thermodynamique. Un diagramme d'équilibre répond à la
question : la réaction chimique est elle thermodynamiquement possible ? Autrement dit, pour
un système binaire (2 composants), un diagramme d'équilibre permet de représenter les
domaines de stabilité des phases et les conditions d'équilibre entre plusieurs phases en
fonction de la température et de la composition.
- 12 -

Une approche empirique permet de dire qu'un diagramme d'équilibre binaire est
constitué de lignes appelées Frontières qui isolent des surfaces appelées Domaines. Un
domaine peut être constitué d'une phase unique, de deux ou de trois phases au maximum [15].
I-5-1-Règle de phase de Gibbs

La règle des phases établit la relation entre le nombre de degrés de liberté V
(variance), le nombre de composants C et le nombre de phases Φ :
V=C-Φ+n (I-1)
La variance désigne le nombre de facteurs d’équilibre indépendants que l’on peut
faire varier arbitrairement sans que le nombre et la nature des constituants physico-chimiques
varient.
n = nombre de facteurs extérieurs (température, pression).
Si on se limite aux diagrammes de phases métalliques, on prend n=1, parce que la
pression dans ce cas n’influe pas de manière signifiante sur les points de transformation et
l’on obtient alors : [16]
V=C-Φ+1 (I-2)
Comme la variance ne peut être inférieure à zéro, ni à un nombre fractionnaire, on en
déduit que :
C - Φ + 1 > 0 soit Φ < C + 1 (I-3)
Autrement dit, le nombre de phases d’un alliage en équilibre ne peut dépasser le
nombre de composants C + l, cela veut dire qu'un système binaire ne peut avoir plus de trois
phases en équilibre, un système ternaire ne peut avoir plus de quatre, etc.
Ø Le cas d’un diagramme de phase binaire

Dans ce cas, il n'y a que deux composants ; on a donc : v= 3- ϕ .
• Cas d'un domaine monophasé, v=2 (exmp. α ): Dans ce domaine on peut faire
varier la température et la composition sans que rien d'autre ne change.
• Cas d'un domaine biphasé, v=1 (exmp. α + β ): Il n'y a plus qu'un seul degré de
liberté. Dès que l'on fixe la température, on en déduit la composition de chacune des deux
phases et leur proportion. Cette déduction peut se faire à l'aide de la règle de l'horizontale et
de la règle des segments inverses.
- 13 -

• Cas d'un domaine triphasé, v=0 : La variance est nulle, il n'y a aucun degré de
liberté donc tous les paramètres sont fixés (entre autre la température). Le domaine se
reconnaît graphiquement car il s'agit d'un segment horizontal (isotherme) appelé palier
invariant. Ce segment comporte trois points fixes qui représentent les compositions des trois
phases invariantes en équilibre.
I-5-2- systèmes binaires isomorphes

Ce cas est le plus simple des diagrammes de phases dont la solubilité solide mutuelle
des deux composants est complète. Exemple d’un système isomorphe : Ni-Cu (Fig. I-6) [17].
La solubilité complète se produit parce que le Cu et Ni ont la même structure cristalline
(CFC), presque le même rayon, électronégativité et Valence [18].
Ce diagramme est caractérisé par un "fuseau de solidification" délimité par deux
courbes correspondant l'une au début de solidification (liquidus) et l'autre à la fin de
solidification (solidus). Ces deux courbes aboutissent aux températures TA et TB qui
correspondent aux températures de solidification (ou de fusion) des deux métaux purs. Le
fuseau compris entre le liquidus et le solidus correspond à la coexistence de deux phases
constituées l'une par un liquide et l'autre par des cristaux formés par une solution solide de
substitution de structure cristalline identique à celle des métaux A et B.
En fait, lorsque deux métaux cristallisent dans la même structure cristalline et ont des
rayons métalliques voisins, ils donnent en général des alliages monophasés. Dans
l’assemblage de l’un des métaux, tous les atomes peuvent être substitués par des atomes de
l’autre métal sans changement de structure cristalline. Dans ce cas, les deux métaux sont
totalement miscibles à l’état solide et forment une solution solide de substitution, quelle que
soit la composition du mélange [19].
Figure- I-6 : Exemple d’un système binaire isomorphe [17].
- 14 -

I-5-3-Systèmes binaires eutectiques

Lorsque les deux constituants A et B n’ont pas la même structure cristalline et/ou que
les rayons métalliques sont très différents, ils sont en général peu miscibles. L’hypothèse
d’une miscibilité nulle est un cas limite (Fig.I-7).
Aux extrémités de ce diagramme existent deux domaines limités en composition, la
solution solide primaire qui présente la solubilité de B dans A (solution solide α) ayant la
structure de A, et l'autre solution solide primaire qui présente la solubilité de A dans B
(Solution solide β) ayant la structure de B. Si l'on s'intéresse par exemple à un alliage de
composition AxBy à l'état solide, à la température To (Fig.I-7) il se présente sous la forme
d'une seule phase formée par la solution solide α. Lorsque la température diminue pour
atteindre la valeur T1,le maximum de solubilité du métal B dans A est franchi lors du passage
par la température Tα et l'alliage présente alors deux phases. L'une, majoritaire, est formée par
la solution solide α de composition définie par le point Ao et l'autre, riche en B, est formée par
la solution solide β et de composition définie par le point Bo. Cette seconde phase précipite en
général aux joints de grains de la solution solide α. Les points M et N correspondent aux
extrémités d'un palier horizontal dit de transformation eutectique passant par le point d'eutexie
E où se rejoignent les deux liquidus partant des points de fusion des métaux purs A et B.
L'alliage de composition CE dénommé "eutectique" est caractérisé par une température de
solidification constante comme les métaux purs, et présente une microstructure caractéristique
constitué par des lamelles alternées des phases α et β de compositions respectives
correspondant aux points Μ et Ν du palier eutectique. Cette constance du point de
solidification est explicable par la règle des phases ou règle de la variance V qui, à pression
atmosphérique constante, s'exprime par V = C + 1 - ϕ, ou C est le nombre de constituants de
l'alliage et ϕ le nombre de phases distinctes par leur état physique ou leur composition. Pour
l'alliage CE à TE trois phases de compositions CE (liquide), CM et CN solides) coexistent, la
variance est donc nulle et le système n'ayant plus de degré de liberté la température est
parfaitement déterminée [19].
- 15 -

Figure-I-7- systèmes binaire eutectique[19]
I-5-4-Systèmes binaires péritectique

Dans ce cas, la température de fusion de l’eutectique vient s’insérer entre la
température de fusion des deux composés A et B purs : le système est appelé péritectique.
Dans certains alliages la solidification s'effectuant par une transformation invariante
entre trois phases comme dans le cas de la transformation eutectique peut conduire à une
microstructure de l'alliage solidifié complètement différente. Cette transformation, dite
péritectique, peut se définir par :
Cristaux α + liquide → cristaux β.
La solidification débute par la formation des cristaux α. Lorsque la température
décroît, le nombre de ces cristaux augmente et le liquide résiduel s'enrichit en composé B.
Lorsqu'une teneur critique est atteinte, le liquide transforme les cristaux α en cristaux β. C'est
le cas par exemple des aciers dont la teneur en carbone est inférieure à 0,5% pour lesquels la
solidification se termine par une réaction péritectique à 1495°C conduisant à la formation de
ferrite δ (Fig. I-8) [19].
- 16 -

Figure I-8 : diagramme d’équilibre Fer-Carbon métastable à cémentite [19].
I -5-5-Systèmes binaires présentant des phases intermédiaires

I -5-5-1-Formation de composés intermédiaires définis Am Bn
Soit un système binaire A−B donnant lieu à la formation d’un composé intermédiaire
de formule AmBn où m et n sont des nombres entiers. Selon la composition globale, on aura
soit un excès de A soit un excès de B par rapport à la composition AmBn. Le diagramme total
peut être considéré comme deux diagrammes juxtaposés : un diagramme A−AmBn et B−AmBn.
On observera, par exemple, la formation de deux systèmes eutectiques ayant chacun leur
propre température de fusion.
Il faut ajouter que AmBn étant un corps pur, on s’attendrait à ce que le sommet des deux
courbes de liquidus se terminant au point de fusion du composé intermédiaire soit anguleux.
Or ce point est arrondi, avec tangente horizontale commune aux deux courbes. Cet arrondi
indique qu’au point de fusion une partie des molécules AmBn se dissocie. Il s’établit un
équilibre de la forme :
Am + Bn -------------> mA + nB
- 17 -

On dit que le composé intermédiaire a un point de fusion congruent ou encore a un

semitectique. Il existe par contre des cas où le composé intermédiaire se dissocie avant de
parvenir à son point de fusion.
La figure I-8 montre un cas où, en plus de deux transformations eutectiques, il y a
formation de deux composés intermédiaires A2B et AB2. On peut voir que le composé A2B est
stable jusqu’à la température de fusion. On dit que A2B a une température de fusion
congruente.
Figure - I-8 : Formation de composés intermédiaires définis [17]
Il faut mentionner qu’en phase liquide le constituant A2B n’existe pas. À ce point de
fusion, point M, ce composé se dissocie en ses éléments A et B. On montre alors que les
courbes de liquidus en ce point de fusion congruente ont des tangentes horizontales.
Au contraire, l’intermédiaire AB2 se décompose avant d’atteindre la fusion. Parvenu à
la température correspondant au point N, cet intermédiaire laisse place à la formation du
liquide et du composé B pur. On dira que le solide AB2 a un point de fusion incongruente[17].
La courbe de liquidus présente aussi un ou plusieurs maxima correspondant à des
compositions définies, mais contrairement au cas précédent l'écart à la composition
stoechiométrique est très réduit. Ces composés sont appelés intermétalliques si A et B sont des
métaux et ont une liaison de type ionique ou covalente qui définit leur structure cristalline. Par
exemple, beaucoup de composés répondant à la formule AB ont une structure de type NaCl
avec alternance des atomes A et B sur un réseau cubique à faces centrées. C'est le cas des
carbures TiC et NbC qui jouent un rôle important dans les aciers inoxydables.
- 18 -

I- 5-5-2-Diagramme binaire avec solution solide secondaire

La courbe du liquidus présente un ou plusieurs maxima correspondant en général à des
compositions définies telles que AB, A2B, AB2, etc. Autour de ces compositions existent des
intervalles correspondant à des alliages monophasés de structure cristalline en général
différente de celle de A et B, d'ou le nom de solutions solides secondaires ou intermédiaires. Si
par exemple le maximum correspond à la composition A2B, Les solutions solides secondaires
situées autour de ce maximum (zone grisée) de la figure-I-9 auront des compositions
correspondant à gauche à A2B1-x.et à droite à A2-xB, notées AxBy sur le diagramme, et présentant
donc des écarts à la stoechiométrie. Suivant les alliages, le domaine d'existence de la solution
secondaire peut atteindre de un à plusieurs %. De part et d'autre d'un maximum se trouvent de
plus deux domaines de transformation eutectique E1 et E2 [16].
Figure I-9 : solution solide secondaire
I -5-6-Systèmes binaires présentant une miscibilité partielle à l’état liquide

Soit le diagramme apparaissant à la figure -I-10. On remarque que le diagramme est
constitué d’un mélange présentant le phénomène de la cristallisation eutectique. L’une des
branches du liquidus est cependant déformée. Une zone de démixtion du liquide en deux
phases l' + l" apparaît.
- 19 -

Figure- I-10 : Diagramme binaire isobare avec miscibilité partielle en phase liquide [17].
Soit un mélange est représenté par le point L de composition xL. En le refroidissant,

en D, apparaissent les premiers cristaux de A pur en équilibre avec le liquide. À la
température correspondant au point M, le solide A est en équilibre avec le liquide de
composition indiquée par la position du point N. Arrivé en E, le liquide de composition xG est
saturé en B : il ne peut plus en dissoudre. Intervient alors la démixtion du liquide : une
deuxième phase liquide de composition xH apparaît. Le système contient trois phases : il est
invariant tant et aussi longtemps que le liquide de composition xG est présent. Lorsque cette
phase a disparu, la température recommence à diminuer et atteint éventuellement la
température de précipitation de l’eutectique [17].
Cela veut dire que les deux métaux à l’état liquide n’étaient pas toujours solubles en
toute proportion, il existe des mélanges qui se séparent en deux liquides dans certaines limites
de concentration et de température. Le plus souvent, chaque métal peut dissoudre une certaine
quantité de l’autre de sorte que les deux couches sont des solutions diluées de A dans B et de
B dans A. l’insolubilité totale est d’ailleurs théoriquement impossible, mais dans certain cas,
rare, comme celui du système Pb-Al on se rapproche beaucoup.
L’équilibre de deux phases est régi, comme celui de deux phases solides, par deux
courbes conjuguées. La solubilité réciproque augmente avec la température critique, les deux
liquides ont la même composition.
- 20 -

I -5-7- Ensemble des équilibres triphasés

La figure I-11 représente un diagramme hypothétique qui rassemble l’ensemble des
équilibres triphasés solide-liquide susceptible d’être observés [20].
L
Eutectique
Monotectique
Métatectique
Péritectique
Syntectique
Eutectoïde
Péritectoïde
Monotectoïde
T
B pur
A pur
xB
Figure I-11: Ensemble des équilibres triphasés solide-liquide [20].
- 21 -

Chapitre II
LES BASES
THERMODYNAMIQUES
DES DIAGRAMMES
DE PHASES
- 21 -

Chapitre II Les bases thermodynamique des diagramme de phases
II-1- Introduction
La description thermodynamique des équilibres entre phase de certains systèmes
notamment les systèmes multi constitués s’avère irréalisable à partir de méthodes uniquement
expérimentales. Les informations résultantes de ces méthodes sont souvent incomplètes. Le
recours à la modélisation numérique est nécessaire. Elle a été rendue possible par les
développements de l’informatique et par la mise au point de logiciels de calcul
thermodynamique.
En effet différents logiciels de calcul thermodynamique permettent maintenant
d’estimer les diagrammes de phase simples, binaires et systèmes ternaires à partir de la
modélisation de leurs enthalpies libres des phases de ces systèmes. La description de ces
systèmes simples prend en compte les données expérimentales disponibles sur chacun d’eux.
Cette étape peut être assistée par l’utilisation de logiciel d’optimisation des paramètres
thermodynamiques nécessaires.
La méthode CALPHAD (CALculation of PHAse Diagram) qui prévaut dans ce
domaine, est très bien décrite dans de nombreux ouvrages tels que celui de Saunders et
Miodowink [2,21]. L’objectif de ce chapitre n’est donc pas d’établir un état de l’art de cette
méthode mais plutôt d’en donner le principe et de présenter les modèles thermodynamiques
qui sont mise en en œuvre dans la description des différents systèmes qui seront étudiés dans
le chapitre IV. Nous nous s’intéressons, en particulier, à étudier le modèle en sous-reseaux,
qui est généralisé dans les logiciels Calphad. Il sera étudié en détail dans le prochain chapitre.
Avant tout ça on va commencer tout d’abord par l’étude de quelques notions de bases de la
thermodynamique nécessaires pour notre étude.
II-2- Premier principe de la thermodynamique :

Pour un système isolé, c’est à dire pour lequel toute forme d’échange d’énergie avec
l’extérieur a été supprimée. L’énergie contenue à l’intérieur du système se conserve
globalement. Le système isolé ne peut donc évoluer qu’à énergie constante.
Le premier principe s’écrit :
E =W +Q (II-1)
Où E : l’énergie interne, et la chaleur (Q) et le travail (W) sont des modes de
transfert de l’énergie dans le cas des systèmes fermés.
En introduisant une nouvelle fonction thermodynamique l’enthalpie H .
- 22 -

H = E + PV (II-2)
et
dH = dE + PdV + VdP (II-3)
P: la pression, V : le volume
Cette fonction représente la chaleur emmagasinée dans le système sous pression
constante. Pour une phase l’enthalpie est décomposable en chaleur sensible qui représente
l’énergie cinétique globale d’agitation des atomes et chaleur latente qui est la somme de
toutes les énergies de liaison interatomique, ainsi une transition de phase, la variation de
l’énergie globale de liaison des atomes représente la chaleur latente de la transition [8].
II-3-Deuxième principe de la thermodynamique :

La thermodynamique introduit une grandeur non conservative, c’est l’entropie S. Elle
représente la mesure quantitative du désordre globale du système [22,23].
La chaleur reçue par le système lors d’une évolution réversible isobare est reliée tout
à la fois à la variation de l’enthalpie et à celle de l’entropie, c’est la relation de Clausius :
dQ=dH=TdS (II-4)
Ce désordre global d’une phase est la somme de plusieurs contributions :
• L’entropie d’agitation thermique (distribution des phonons);
• L’entropie de configuration (distribution des atomes sur les sites);
• L’entropie électronique (distribution des électrons sur les états quantiques);
• L’entropie magnétique (distribution des spins demi-impairs).
Chacun de ces termes s’exprime par une formule de type Boltzmann [6] :
S = k ln(ϖ eq) (II-5)
avec
ϖ eq : est le nombre de distributions microscopiques à l’équilibre de la population
concernée (phonons, atomes, électrons, spins), qui représentent le même état macroscopique,
et k la constante de Boltzmann.
II-4-Troisième principe de la thermodynamique :

Un troisième principe de la thermodynamique posé par Nernst en 1906, attribue l’entropie
nulle s=0 à toute phase en état d’équilibre stable ou métastable, au zéro de l’échelle de
température thermodynamique.
- 23 -

En introduisant deux fonctions d’état

-l’énergie libre F=E-TS (II-6)
-l’enthalpie libre de Gibbs G=H-TS (II-7)
Ce qui conduit
dF =-P dV—SdT—TdS (II-8)
dG =VdP—SdT—TdS (II-9)
F : est une fonction thermodynamique à V et T constants.
G : est une fonction thermodynamique à P et T constants.
T : la température thermodynamique
Pour les phases condensées (solide et liquide), on peut confondre les deux fonctions F et G
(tant que l’on reste dans les domaines de pression pas trop élevée (lorsque H….> E, G ….>F)
[24].
II-5-Potentielle chimique :
Cette grandeur permet de caractériser la présence d’un constituant dans un milieu. Sa
valeur permet de prévoir comment va s’opérer le déplacement de matière au cours d’une
évolution spontanée : un constituant se déplacera spontanément de la phase où son potentiel
chimique est plus élevé, vers la phase où il est plus faible, pour établir l’équilibre
thermodynamique caractérisé par l’égalité des potentiel chimique entre les deux phases [24].
La composition d’un système physico-chimique peut varier :
• soit par suit d’une transformation physico-chimique (seul possibilité si le système est
fermé )
dE = dW + dQ = -PdV + TdS (II-10)
• soit, en plus, par suite d’un transfère de matière avec l’extérieur si le système est
ouvert
dE = dW + dQ + ∑ µ i dni = − PdV + TdS + ∑ µ i dni (II-11)
Le premier terme représente l’énergie mécanique reçue par le système (-PdV), le deuxième
représente la chaleur reçue par le système (TdS) et le troisième représente l’énergie chimique
reçue par le système lors de l’échange du constituant i ( µi dni ).
Soit un système qui comporte i constituants différents, chaque constituant comporte ni moles
L’enthalpie libre qui est une fonction d’état, est maintenant fonction de T, P et ni.
D’où
- 24 -

dG = ∂G ( )
∂P T , ni
(
dP + ∂G
∂T
)
P , ni
dT + ∑  ∂G  dni
 ∂ni T , P
(II-12)
On définit le potentiel chimique µi par :
µ i =  ∂G 
 (II-13)
 ∂ ni T , P
d’où
dG = VdP − SdT + ∑ µ i dni (II-14)
De la même manière, on peut écrire :

dH=TdS+VdP+ ∑ µi dni (II-15)
dF=-pdV-SdT+ ∑ µi dni (II-16)
Des relations précédentes, découlent quatre expressions égales du potentiel chimique, qui
peuvent chacune lui servir de définition :
µi =  ∂G  =  ∂E  =  ∂H 
 =  ∂F  (II-17)
 ∂ni T , P , n j  ∂ni  S ,V , n j  ∂ni S , P,n j  ∂ni T ,V , n j
II-6-L’activité:
Si deux éléments A et B sont susceptibles de former une solution solide ou liquide,
trois éventualités peuvent être envisagées:
Ø Les atomes de même nature ont des forces d’attraction égales à celle des atomes de
nature différente. Du point de vu chimique les atomes A et B sont indiscernables.
Ø Les atomes de même nature s’attirent plus que les atomes de nature différente. Il y’a
tendance à la ségrégation.
Ø Les atomes de nature différente s’attirent plus que les atomes de même nature, il y’a
tendance à la formation d’un ordre.
Les éléments qui constituent une solution liquide ou solide ont une tension de vapeur qui, à
l’état pur, est fonction de la température appelée P0. Lorsqu’un élément contient un autre
élément en solution sa tension de vapeur n’est plus P0 mais P inférieur à P0. L’activité de
l’élément i dans la solution est définie par la relation:
pi
ai = (II-18)
p i0
Dans le cas des solutions a i = γ i X i ( γ i est le coefficient d’activité dans le système) [25].
Le potentiel chimique s’écrit :
µ i = RTLnai (II-19)
- 25 -

II-7-Théorème d’Euler
Soit une fonction homogène f ( X , Y , Z ) de degré p :
f (λX , λY , λZ ) = λ p f ( X , Y , Z ) (II-20)
df (λX , λY , λZ ) ∂f ∂λX ∂f ∂λY ∂f ∂λZ
− . + . + . = pλ p−1 f ( X , Y , Z ) (II-21)
dλ ∂λX ∂λ ∂λY ∂λ ∂λZ ∂λ
Où
 ∂f   ∂f   ∂f 
 = pλ f ( X , Y , Z )
p −1
X  + Y  + Z (II-22)
 ∂λ X   ∂λ Y   ∂λ Z 
Comme cette relation est varie pour toute valeurs de λ , on peut aussi poser :
 ∂f   ∂f   ∂f 
X  + Y  + Z  = pf ( X , Y , Z ) (II-23)
 ∂λ X   ∂ λ Y   ∂λ Z 
dans le cas d’une fonction d’état extensive, c'est-à-dire homogène de degré 1, on a :
 ∂f   ∂f   ∂f 
X  + Y  + Z  = f ( X ,Y , Z ) (II-24)
 ∂λ X   ∂λ Y   ∂λZ 
Appliqué à la fonction G, la relation précédente (relation d’Euler) donne :
G = ∑ ni G i = ∑ ni µ i (II-25)
i i
G i : est l’enthalpie libre molaire partielle du constant i[26] .
II-8-Relation de Gibbs- Duhem

Pour une transformation élémentaire d’un système ouvert, la variation d’enthalpie libre est :
dG = VdP − SdT + ∑ µ i dni (II-26)
i
A partir de l’expression de l’enthalpie libre G = ∑ ni µ i , on obtient :

i
dG = ∑ µ i dni + ∑ ni dµ i (II-27)
i i
Par identification, il vient :
∑ n dµ
i
i i = VdP − SdT (II-28)
C’est la relation de Gibbs- Duhem. En particulier, pour une transformation à température et

pression constantes :
∑ n dµ
i
i i =0 (II-29)
- 26 -

Pour un système à deux constituants notés avec les indices 1 et 2, la relation s’écrit :
n1 dµ1 + n2 dµ 2 = 0 (II-30)
ni
0ù la fraction molaire X 1 =
∑ ni
i
X 1 dµ1 + X 2 dµ 2 = 0 (II-31)
Cette expression permet en particulier de calculer le potentiel chimique µ 2 lorsque
l’expression µ1 ( X 1 ) est connue. [26].
II-9- Relation de Gibbs- Helmholtz

On va étudier la variation du potentiel chimique avec la température. L’égalité des dérivées
secondes de G par rapport à la température et à la quantité ni donne
 ∂µ i   ∂  ∂G    ∂  ∂G    ∂S 
  =    =    = −  = −S i (II-32)
  
 ∂T  P ,ni  ∂T  ∂ni T , P  P ,ni  ∂ni  ∂T  P ,ni  T , P  ∂ni  T , p,n j ≠i
Où S i est l’enthalpie molaire partielle du constituant i

µi
Il est souvent plus facile d’exprimer les variations de en fonction de la température. On
T
utilise pour cela la relation de Gibbs- Helmholtz, qui se dérive comme suit :
 ∂G 
H = G + TS = G − T   (II-33)
 ∂T  P
 G
∂  H
qui s’écrit aussi  T  = − 2 Relation de Gibbs- Helmholtz
 ∂T  T
 
 P
La dérivation des deux membres de cette relation par rapport à ni donne :
 µi 
∂  Hi
 T  =− 2 (II-34)
 ∂T  T
 
  P , ni
Où H i est l’enthalpie molaire du constituant i [26].
II-10-L’énergie de Gibbs molaire et le potentiel chimique

L’énergie libre de Gibbs d’un système à température et pression constantes est la
somme des potentiels chimiques de ses constituants (formule (II-25)).
- 27 -

L’énergie libre molaire de Gibbs Gm est l’énergie d’une mole de substance, à température et
pression constante s’écrit :
G m = ∑ X i Gi = ∑ X i µ i (II-35)
X i est la fraction molaire de l’élément i
C’est le système contient deux éléments A et B, alors X A + X B =1, l’énergie molaire

libre de Gibbs s’écrit :
Gm = X A G A + X B G B = X A µ A + X B µ B [24] (II-36)
II-11-Critères d’équilibres
Un système fermé, isotherme et isobare est à l’équilibre thermodynamique lorsque sa
fonction enthalpie libre de Gibbs G (plus simplement fonction de Gibbs) est minimale. Cette
fonction de Gibbs est la somme de deux contributions (formule (II-7)).
Cependant l’évolution interne du système obéit à deux impératifs contradictoires la
tendance énergétique voudrait que l’enthalpie s’écoule au maximum vers l’extérieur (dH<O),
tandis que la tendance entropique voudrait que le désordre s’accroisse au maximum (dS>O).
Selon la relation (II-4), les variations de H et S ne sont pas indépendantes et c’est au
minimum de la fonction G que la balance s’équilibre. Toutefois en considérant
l’équation (II-7), on déduit que vers les basses températures l’influence entropique sera faible
et que c’est la tendance enthalpique qui régira l’équilibre, tandis qu’aux très hautes
températures c’est la tendance au désordre qui l’emportera[8].
Pour un système totalement isolé, Clausius a posé comme second principe de la
thermodynamique que l’état d’équilibre thermodynamique est atteint lorsque l’entropie du
système est maximale. Boltzmann a montré ultérieurement que ce maximum d’entropie
représentait le nombre maximal de distributions microscopiques possibles, pour un même état
macroscopique, ce qui équivaut à placer le système dans son état statistiquement le plus
probable [6].
II-12-Stabilité d’une phase

Un extremum se traduit toujours par une différentielle nulle : pour un système fermé,
isotherme et isobare, c’est la fonction de Gibbs G qui présente cet extremum et l’égalité
dG = 0 est appelée condition d’équilibre. Mais un extremum peut être un maximum ou un
minimum. Le premier cas représente un équilibre instable qui ne peut se réaliser, le second un
- 28 -

équilibre stable vers lequel le système tend à s’établir. Les conditions de stabilité concernent
donc les dérivées secondes de la fonction G et traduisent le fait que la forme quadratique
d 2 G est positive lors des échanges internes locaux de chaleur, de volume et de matière. A ces
échanges internes locaux correspondent respectivement les fluctuations locales de
température, de densité et de potentiel chimique. Toutes ces fluctuations doivent rester stables
et réversibles. L’équilibre thermodynamique représente donc une moyenne statistique d’états
microscopiques en fluctuation et qui reste stable lorsqu’elle est prise sur un intervalle fini de
temps.
Pour un système binaire les conditions de stabilité s’expriment par trois inégalités qui
doivent être simultanément satisfaites
∂T ∂P ∂µ i
≥0 ≤0 ≥0 (II-37)
∂S ∂V ∂ni
Avec S l’entropie, T la température, P la pression, V volume, µi . le potentiel

chimique et n le nombre de moles du constituant i dans le système [6].
Les variables qui peuvent être bloquées dans ces dérivées partielles sont assez
variées, par exemple la première inégalité est valable en condition isochore ou en condition
isobare, la seconde en condition isotherme ou en condition isentropique, etc. Ces conditions
de stabilité expriment trois réalités physiques de l’équilibre d’un système : son désordre
augmente avec la température, son volume diminue lorsque la pression s’accroît, le potentiel
chimique d’un de ses constituants augmente lorsque l’on ajoute des moles de ce constituant au
système.
Un état métastable satisfait toutes ces conditions, c’est donc un état d’équilibre
thermodynamiquement stable, mais qui s’établit en l’absence d’une phase dont la présence
engendrerait un équilibre de plus basse enthalpie libre G; la raison de l’absence de cette phase
est qu’elle germe très difficilement dans le système. Le passage spontané de l’équilibre
métastable à l’équilibre stable est possible à chaque instant et imprévisible, le passage inverse
est impossible [6].
Pour un système binaire, si on considère deux phases α et β à équilibre qui sont
associées respectivement aux G α et G β , l’énergie libre total du système s’écrit :
dG = dG α + dG β (II-38)
à température et pression constantes on a les égalité suivantes :
dGα = µ 1α dn1α + µ 2α dn2α
 β (II-39)
dG = µ 1β dn1β + µ 2β dn2β
- 29 -

Comme la quantité totale de A et B est constante on a :
dn1α = − dn1β et dnα2 = − dn2β (II-40)
dG = µ 1α dn1α + µ 2α dn2α − µ1β dn1α − µ 2β dnα2

 (II-41)
( ) ( )
dG = µ 1α − µ 1β dn1α + µ 2α − µ 2β dn2α
Comme dn1α et dn2α sont indépendantes
A l’équilibre, dG = 0 ⇒ µ1α = µ1β et µ 2α = µ 2β
Pour chaque élément, légalité du potentiel chimique dans toutes les phases (pas de transfert
chimique) [24].
II-13- Construction de la tangente commune et le potentiel chimique

La représentation graphique des potentiels chimiques des deux éléments par la
méthode d’intersection est schématisée sur la figure II-1.
L’intersection des deux axes par la tangente de la courbe d’énergie libre de la phase α
à la composition X 2α représente µ1α et µ 2β . Même chose pour la phase β .
Les conditions d’équilibres sont équivalentes à la tangente commune des deux courbes
d’énergie libre :
µ1α = µ1β et µ 2α = µ 2β ⇔ Tangente commune des deux courbes [24].
(a) α et β ne sont pas en équilibre (b) α et β en équilibre
Figure-II-1 : Représentation de la tangente commune
- 30 -

II-14- Détermination des courbes de solidus, liquidus

En considérant un système binaire A-B : α (G α ) et β (G β ) deux phases solides et L(G L ) une
phase liquide. Pour chaque température, la tangente commune entre G α et G L ou entre
G α et G β donne la composition des deux phases en équilibre [24].
II-14-1-La miscibilité totale de A et B

La figure-II-2g correspond à un système simple où les deux constituants sont
miscibles en toutes proportions dans la phase liquide L comme dans la solution solide α .
Deux phases seulement occupent le diagramme d’équilibre. A haute température, la phase L
et, à basse température, la phase α . Les corps purs isomorphes A* et B* fondent
respectivement à TA et TB. Entre ces deux températures, le diagramme présente deux lignes
conjuguées liquide et solide appelées liquidus et solidus formant un fuseau de coexistence des
phases qui se referme aux extrémités du diagramme sur les deux points de fusion des corps
purs.
Les enthalpies libres de mélange calculées à T1, T2, T3,T4,T5 et à T6 sont tracées sur la
figure-II-2. A T1 et T6 °C, les courbes L et α sont totalement disjointes et la phase α et L
sont stables respectivement à toute composition. Au contraire à T3 et T4 °C, les courbes L et
α se coupent au voisinage de la composition X et une tangente commune définit un segment
de conjugaison liquide *<—Solide. Le segment de tangente commune α L de la figure-II-2b
constitue la plus basse enthalpie libre réalisable entre ces deux compositions ; le mélange
biphasé ( α + L) est donc, dans ce domaine, l’équilibre stable. A sa gauche, le domaine α
monophasé et à sa droite le domaine L monophasé. Dans le domaine biphasé et pour une
même température, les potentiels chimiques sont indépendants de la composition (ils se lisent
à l’intersection de la double tangente avec les axes verticaux des deux corps purs (Fig.II-2b).
Aux points de fusion des constituants A et B purs (TA et TB) on remarque: G AL = G αA à T = TA
et G BL = G Bα à T = TB [15].
II-14-2- Gap de miscibilité dans l’état solide

La figure-II-3 illustre la phase de démixtion qui se réalise à T<Tc (température
critique). La figure-II-3b montre les équilibres réalisés à T1 et T2, températures où
apparaissent deux segments de conjugaison séparés, le segment de démixtion solide → solide
et présente le domaine d’existence des deux phases isomorphe alors que à T>Tc il n’ya plus
des deux phases, il existe seulement la nouvelle phase solide α [24].
- 31 -

II-14-3-Le point eutectique sans isomorphes

Les deux solutions solides ne sont pas isomorphes ( α et β ont une structure
cristalline différente). La figure-II-4-d représente les enthalpies libres de mélange des trois
phases α , β et L, à la température invariante eutectique TE. La tangente triple touche une
seule fois chacune des courbes. La figure-II-4,b,c ,e montrent les courbes d’enthalpie libre de
mélange au-dessus et en dessous de la température eutectique TE. Au-dessus de TE, il existe
deux équilibres distincts du liquide soit avec α soit avec β , tandis qu’en dessous de TE le
liquide n’est plus en concurrence pour l’équilibre. La disposition des courbes d’enthalpie libre
pour les trois phases ( α , β , et liquide) aux différentes températures figure-II-3 permet
d’établir le diagramme de phases eutectique [15].
II-14-4-Le point eutectique dans le cas isomorphe

La figure-II-5 représente un équilibre eutectique très semblable à celui de la
figure-II-4, mais ici les deux solutions solides α1 ,α 2 sont isomorphes (même structure
cristalline) et proviennent de l’équilibre de démixtion solide →* solide. Mais la situation
eutectique est plus générale et peut concerner deux solides quelconques. Lors de la
solidification du liquide eutectique (celui de l’équilibre invariant à trois phases), il se forme
simultanément deux phases solides distinctes. La réaction s’opère à température constante. Le
liquide eutectique est celui qui fond à plus basse température. La disposition des courbes
d’enthalpie libre pour les trois phases ( α1 ,α 2 , et liquide) aux différentes températures
(Fig.II-5) permet d’établir le diagramme de phases eutectique [13].
- 32 -

Figure-II-2 : Etablissement du diagramme de phases d’un système présentant une

miscibilitée totale à partir des courbes des enthalpies libres
(a,b,c,d,e,f) Courbes d’enthalpies libres
(g) Diagramme d’équilibre de phases correspondant. [24]
Figure-II-3 : Etablissement du diagramme de phases d’un

système présentant une Gap de miscibilité dans l’état solide
(a) courbe d’enthalpie libre
(b) Diagramme d’équilibre de phases correspondant [24].
solvus
- 33 -

Figure-II-4 : Etablissement du diagramme de phases d’un système non isomorphe

présentant une miscibilité
limitée à partir des courbes des enthalpies libres
(a,b,c,d,e) Courbes d’enthalpies libres (f) Diagramme d’équilibre de phases eutectique [24]
Figure-II-5 : Etablissement du diagramme de phases d’un système isomorphe présentant une

miscibilitée limitée à partir des courbes des enthalpies libres
(a,b,c,d,e) Courbes d’enthalpies libres
(f) Diagramme d’équilibre de phases eutectique.[24].
- 34 -

II-15- Présentation de la méthode CALPHAD

La méthode CALPHAD (Calculation of Phase Diagrams) est basée sur la
minimisation de l’enthalpie libre, et elle regroupe plusieurs modèles thermodynamiques. La
figure II.6 est un organigramme simplifié des différentes étapes de cette méthode pour
calculer les diagrammes de phases. La description thermodynamique du système est basée sur
des données expérimentales (thermodynamiques et diagrammes de phases). Le terme
« description thermodynamique » signifie que les paramètres des modèles de chaque phases
existant dans le système sont obtenus, donc non seulement les propriétés thermodynamiques
des phases qui seront calculés mais aussi le diagramme de phases [6 ].
Figure II.6 : La méthodologie Calphad.[21]
- 35 -

Le calcul des systèmes multiconstitués avec plus de constituants peut être obtenu à partir
d’une description des systèmes avec moins de constituants en utilisant des méthodes
d’extrapolation [27]. La stratégie utilisée dans l’évaluation des systèmes multiconstitués est
représentée sur la figure II.7. En premier, une description thermodynamique des systèmes
binaires est obtenue, ensuite les fonctions thermodynamiques de ces systèmes sont
extrapolées aux systèmes ternaires. Les résultats obtenus sont alors comparés aux données
expérimentales et si nécessaire d’autres fonctions sont ajoutées.
Système Binaire Evaluation (G )

exc
bin
Système ternaire Extrapolation (∑ G ) exc

bin
+ Evaluation (G )
exc
ter
Système quaternaire Extrapolation (∑ G exc

bin + Gter
exc
)
+ Evaluation (G )
exc
qua
.
.
.
.
Figure II.7 : Stratégie utilisée par la méthode Calphad pour les systèmes multiconstitués.
Les informations expérimentales peuvent se classer en deux catégories selon leur

nature. D'une part les informations concernant les grandeurs thermodynamiques et d'autre
part, les informations purement diagrammatiques [21].
Ø Les grandeurs thermodynamiques

Ces grandeurs qu'elles soient partielles ou intégrales, sont directement liées à
l'enthalpie libre ou à sa composante enthalpique et sont donc les informations les plus
pertinentes
- les potentiels chimiques sont obtenus principalement par des techniques d'équilibres
en phase gazeuse (mesure de pression de vapeur, cellule d'effusion de Knudsen…) ou bien par
des techniques potentiométriques (mesures f.é.m de cellules électrolytiques à électrolytes
liquides ou solides) ;
- 36 -

- la calorimétrie permet d'obtenir les chaleurs spécifiques, les enthalpies de réaction et

de formation. Dans ce dernier cas, la mesure peut se faire soit par réaction directe (ex :
calorimétrie de chute) ou alors par réaction indirecte (ex : calorimétrie de dissolution).
Les résultats des méthodes d’estimation par calculs ab-initio qui deviennent de plus en plus
nombreuses et fiables sont également pris en compte dans l’optimisation. Ils sont précieux car
ils permettent d’obtenir des informations dans des situations pour lesquelles l’expérience est
rendue difficile voire impossible comme dans le cas des états métastables [21].
Ø Les données diagrammatiques

Ces données proviennent des diagrammes de phases et traduisent des situations
d'équilibres entre phases. Ce type de données est généralement plus abondant que le
précédent. En effet, elles sont plus facilement accessibles expérimentalement. La plus
fréquemment rencontrées sont :
- les mesures d'analyses thermiques (ATD, DSC …) [21].
II-16- Description thermodynamique de phases d’un diagramme d’équilibre

II-16-1-Les éléments purs
La description de l’enthalpie libre des corps purs nécessite de définir son états de
référence celui couramment utilisé provient de formalisme développé par le SGTE
(Scientifique Groupe Thermo data Europe) qui recense, établit et compiles les descriptions
thermodynamiques des corps purs dans les bases de données. Selon le SGTE, l’état de
référence noté SER (stable élément référence) qui correspond à la structure stable de l’élément
i à 298.15K (sous pression de 105 Pa)
L’évolution avec la température de l’enthalpie libre molaire ainsi référencée d’un
élément i dans son état ϕ (stable ou métastable) est donnée par GHSERi afin de corriger
l’enthalpie libre molaire de i de cette différence de température.
GHSERi = °Giϕ − ° H iSER [28] (II-42)
avec
°Giϕ (T ) = a + bT + cTLnT + dT 2 + eT 3 + fT −1 + hT −9 [28] (II-43)
ϕ
° G i (T ) : L’enthalpie libre molaire de l’élément i dans la structure ϕ à la température T
0
H
SER
i
: L’enthalpie de formation de l’élément i dans son état SER à 298.15K
- 37 -

. Pour représenter ces termes dans un domaine de température étendue il est souvent
nécessaire de définir pour un élément donné plusieurs intervalles de température avec des
valeurs différentes des coefficients a,b,c,dn pour prendre en compte les intervalles de stabilité
des divers phases des corps purs [2].
Notons que les changements d’état de référence (stabilité de réseau) entre deux
phases ϕ et ϕ ' peuvent s’exprimer par une équation du même type mais se résument
généralement à une variation linéaire :
Giϕ − Giϕ = a + bT (Variété allotropique)

'
(II-44)
II-16-2- Phases stœchiométriques (Composés définis)

L’énergie libre molaire de formation d’un alliage binaire imjp s’exprime sous la forme:
°G mm p (T ) = mGHSERi + pGHSERj + ∆ f G m
i j i jp
(II-45)
Où le terme ΔfGimjp représente l’énergie libre molaire de formation du composé

intermétallique
∆fG m =∆fH − T∆ f S
i jp im j p im j p
[28] (II-46)
Une telle expression présente des avantages lorsqu’il n’existe pas de données
expérimentales de capacité calorifique pour le composé considéré. De ce fait, les termes
enthalpiques ΔfHimjp et entropiques ΔfS im jp sont considérés comme indépendants de la
température.
En effet la capacité calorifique du composé est alors simplement égale à la somme
pondéré des capacités calorifiques des éléments purs. Toute fois, la capacité calorifique du
composé étant liée à son enthalpie libre par une relation similaire à la précédente, les valeurs
expérimentales de Cp de ce composé peuvent également être utilisée dans les procédures
d’optimisation en exprimant son enthalpie libre molaire de la manière suivante [19] :
°
Gmm p − mGHSERi − pGHSERi = a + bT + cTLnT .....
i j
(II-47)
II-16-3-phase étendu (Modèle polynomial)

Le modèle présent ci-dessous convient particulièrement à la description des solutions
solide primaires, de certaines solutions intermédiaires présentant un large domaine de non-
stœchiométrie ainsi que la phase liquide rencontrée dans les systèmes intermétalliques.
L’enthalpie libre molaire de ce type de phase s’exprime à partir de trois
contributions :
- 38 -

ϕ ref ϕ id ϕ ex ϕ
G (T )= G + G + G
m
(II-48)
ref ϕ
G est l’enthalpie libre de référence prenant en compte les corps purs pris avant le
mélange.
= ∑ xi G i
ref ϕ ϕ
G (II-49)
i
Le terme de mélange idéal idGϕ est donné par l’entropie de configuration idéale d’un
mélange aléatoire des éléments A et B, appelé mélange idéal :
RT ∑ x Lnx
id ϕ
G = (II-50)
i i
i
ex ϕ
Le terme d’excès G sur lequel porte l’essentiel de la modélisation est décrit dans le
cas d’une solution binaire par le modèle de Redlich-Kister [29]:

ex ϕ ϕ ϕ ϕ
G = x xLi j i, j
(II-51)
Le terme Lϕi, j introduit dans cette équation traduit l’interaction entre les atomes A et B
dans la solution. Dans le cas des solutions dites strictement régulière, il ne varie pas avec la
composition et peut être relié aux énergies des liaisons i-i,j-j et i-j.
Lorsque Lϕi, j est positif, les atomes A et B ont tendance a se repousser et la solution
présente une lacune de miscibilité dans un certain intervalle de T.

Lorsqu’il est négatif, le mélange des atomes A et B stabilise la solution et peut
conduire à une mise en ordre.
Lorsqu’il est faible, on tend vers un mélange idéal.
Dans la plupart des cas que nous étudierons un terme d’interaction indépendant de la
composition ne permet pas de rendre compte de l’ensemble des données expérimentales. Il a
le plus souvent été exprimé sous la forme d’un polynôme de Redlich-Kister :
∑ Lν ,ϕ ( xϕi − xϕj )
υ
ϕ
= (II-52)
L i, j
υ
i, j
Les paramètres
Lν ,ϕ peuvent dépendre du degré de développement ν (entier positif) du
A, B
polynôme, et en général une fonction linéaire de la température, ils s’écrivent :
Lν ,ϕ i , j = aνi , j + biν, j T (II-53)

ex
G ϕ = xi x j ∑ Lνi ,,ϕj ( xi − x j )ν
ν
Où chaque coefficient binaire est partagé entre un terme d’interaction enthalpique et un terme
entropique d’excès [2].
- 39 -

II-16-4-Phases intermédiaires
L’ecart à la stoechiométrique est très possible dans certains cas où la phase ordonnée
devenant étendu, et le modèle utilisé pour décrire ce type de phase est dit modèle des sous-
réseau qui sera étudié en détaille dans le prochain chapitre.
- 40 -

Chapitre III
PRESENTATION
DU MODELE DES
SOUS-RESEAUX
- 41 -

III-1-Introduction
Le modèle de la solution régulière (Hildebrand) est inapte pour décrire les phases qui
présentent un ordre important comme les composés intermétalliques ou ioniques et liquides
ioniques. Pour cela on a choisi un autre modèle pour décrire ce type de phase, appelé modèle
des sous- réseaux.
Une solution ordonnée est en réalité peut être non-stoechiomitrique, elle peut
présenter un domaine de stabilité bien étendu. Les différents éléments constituants ces
solutions adaptent préférentiellement des sites cristallographiques différents comme la montre
la figure-III-1 [30] qui représente l’interpénétration de deux ou plusieurs sous-réseaux. Le
modèle en sous-réseaux sera utilisé pour décrire ces solutions qui peuvent admettre
différentes espèces chimiques (ion, molécule, atomes, lacunes), pour les solutions solides
métalliques chaque type de site cristallographique est appelé sous réseau.
Nous introduirons tout d’abord ce modèle dans le cas général en prenant soin de
définir tous les termes impliqués puis nous étudierons quelques cas particuliers que le modèle
permet de décrire dans les systèmes métalliques.
Figure-III-1: Structure cubique représente l’interpénétration des deux réseaux

cubiques à face centrée.
- 42 -

III-2- Modèle à plusieur sous-réseaux

III-2-1- Définition de fraction de site
Dans cette section les concepts généraux de (CEF) seront définis. Le modèle à deux
sous-réseaux est utilisé dans la section 2 pour illustrer le modèle des sous-réseaux.
La quantité du composant i dans le sous-résau s, est donnée par la fraction de site.
Cette fraction de site est défini comme suit [30] :
nis
y is = (III-1)
Ns
où
nis : est le nombre d’atome de l’élément i dans le sous réseau(s).
et
N s : est le nombre total de site dans ce sous- réseau.
Dans le cas d’un sous-réseau substitutionnel présentant dans une solution solide
interstitielle et qui est occupé par plusieurs genres d’atomes du métal, la fraction de site de ces
atomes du métal est donnée par la relation (III-1).
Par exemple, Pour l’alliage Fe-Cr-C, la fraction de site du fer dans cet alliage est donnée par [31]
n Fe
y Fe = (III-2)
N
L’exposant a été négligé dans l’équation (III-2) car tous les paramètres ont été
considérés pour un sous-réseau substitutionnel. Evidemment, pour un sous-résau
substitutionnel [31] :
n
∑y
i =1
i =1 (III-3)
On peut généralisé le concept de fraction de site en incluant les sites vacants qui sont
importants de les considérer dans le cas des phases interstitielle alors l’équation (III-1)
devient [31] :
nis
yis = N s = nvas + ∑ nis (III-4)
n vas + ∑ nis i
i
Où le nombre totale de site N s dans le sous réseau (s) est égale au nombre totale
d’atome nis plus le nombre totale de site vacants nvas existant dans ce sous-réseau [32] .
- 43 -

On considère que la solution solide interstitielle décrite par un sous réseau

interstitielle et contient seulement du carbone et des lacunes
Exemple : la fraction de site des atomes de carbone est exprimée par [31] :
nc
yc = (III-5)
nva + nc
et
la fraction de sites vacants s’écrit :
nVa
yVa = (III-6)
nVa + nc
Pour un sous-réseau interstitielle : y c + yVa = 1 (III-7)
La fraction molaire xiϕ de l’élément i dans la phase ϕ peut être calculée à partir de fraction de
site comme suit [31] :
∑N y s
i
s
xiϕ =
∑ N (1 − y )
s
s s
(III-8)
va
s
∑y
i
i
s
= 1 − y vas (III-9)
Lorsque l équation (III-8) est appliquée dans le cas des solutions interstitielles, le
rapport(r) entre le nombre de sites interstitiels N i et le nombre de sites substitutionnels
N s est nécessaire de le considérer pour différentes matrices.
Les matrices plus intéressantes sont :cfc, bcc, hcp.
Les différentes valeurs du rapport s’écrivent :
(N N )i
s =1
(N N )
Fcc
i
s =3 (III-10)
bcc
(N N )i
s
hcp
=
1
2
Pour le sous-réseau interstitiel, la relation entre la fraction molaire et la fraction de
site est s’écrit :
xi
yi = (III-11)
r (1 − xi )
ry i
xi = (III-12)
1 + ry i
- 44 -

i
Où r=N (III-13)
Ns
yi et xi sont respectivement la fraction de site et la fraction molaire de l’élément i en
interstice, comme par exemple le carbone, dont l’élément i ne présente que les éléments
d’alliage et non les sites vacants car xi est seulement référée aux éléments en excluant
conventionnellement les sites vacants.
Pour le sous-réseau substitutionnel, la relation entre la fraction molaire et la fraction
de site est s’écrit aussi :
xj
yj = (III-14)
1 − xi
yj
xj = (III-15)
1 + ryi
avec
x j et y j sont respectivement la fraction de site et la fraction molaire des éléments , comme le
fer, chrom et nickel [31] .

L’occupation des sous-réseaux est présenté comme suit (I, J)(K,L) …(M,N), où I,J
sont les éléments qui occupent le premier sous-réseau, K et L pour le deuxième etc …
Dans cette notation les sites vacants sont considérés comme l’un des éléments du
sous-réseau.
GI , J : K : L : désigne l’énergie de Gibbs d’une phase modélisée par trois sous-réseaux,
dont les éléments I,J occupent le 1er sous-réseau , alors que K et L occupent respectivement
le 2ème et le 3ème.
La notation utilisée ici pour indiquer l’occupation des sous-réseaux est la suivante:
les deux points (:) séparent des éléments qui ne sont pas dans le même sous-réseau
alors que la virgule (,) sépare deux éléments interagissant sur le même sous-réseau [30].
III-2-2- L’entropie ideal de mélange

S ideal : est l’entropie idéal qui traduit la répartition aléatoire des éléments sur chacun
de sites. Le nombre d’arrangement possible des espèces pressentant sur le sous-réseau(s)
est :
N (s)
wp = (s)
(III-16)
∏n
i =1
(s)
i
- 45 -

D’où le nombre total d’arrangements possible sur tous les sites s’écrit [1]:
N s!
Wp = ∏ s
(III-17)
s
∏n !
i =1
i
s
ü Pour l’entropie idéale du mélange, on considère que la contribution configurationnelle.

L’enthalpie libre molaire de Gibbs idéal est donnée par:
id
G ϕ = −TS id = RT ∑ N s ∑ y is Lnyis (III-18)
s i
où
yis : représente aussi la contribution des sites vacants.
ü Pour un sous-reseau substitutionnel, l’enthalpie libre de Gibbs idéal est donnée par:
id
G ϕ (substitutionnelle) = −TS id = RT ∑ yi Lny i (III-19)
i
ü Pour un sous-reseau interstitiel (carbon, nitrogène, vacant), l’enthalpie libre de Gibbs

idéal est donnée par :
id
G ϕ (interstiti el) = −TS id = rRT ( yc Lnyc + yVa LnyVa ) (III-20)
Par conséquent, d'après l’équation (III-18), Pour une solution interstitielle l’enthalpie libre
molaire de Gibbs idéal s’écrit [31] :
id
Gϕ =id G ϕ (substitutionnelle)+ id G ϕ (interstiti el )
= −TS id (III-21)
id
G ϕ = RT ∑ yi Lnyi + rRT ( yc Lnyc + yVa LnyVa )
i
III-3- modèle à deux sous-résau

III-3-1-L’enthalpie libre molaire
Le calcul de l’enthalpie libre de Gibbs par le modèle de sous-réseau est illustré ici
pour le cas hypothétique avec deux sous-réseau (A,B)a(C,D)b. Encore l’énergie libre molaire
peut être présentée sous trois formes de contributions, c’est le cas similaire à celui employé
pour les solutions de substitutions :
Gmϕ 2 s = ref G2ϕs + id G2ϕs + exG 2ϕs (III-22)
Gm2s : représente l’enthalpie libre pour deux sous-réseaux
- 46 -

Si l’élément A est substitué par l’élément B sur le premier sous-réseau et même chose pour les
deux éléments (D,C), les fractions de sites sont liées à la composition globale par les relations
suivantes :
a a b b
x ϕA = yA x ϕB = yB xCϕ = yC x ϕD = y D (III-23)
a+b a+b a+b a+b
III-3-1-1- L’enthalpie libre de référence

L’enthalpie libre de référence dans ce cas est définit pour les composés définis. Ces
composés sont obtenus lorsque chaque sous-résau est occupé par un seul élément. La
description d’une tel phase par le modèle (A,B)a(C,D)b, nous permet d’avoir quatre composés
AC, BC, AD et BD [32]. La figure-III-2-a indique une représentation bidimensionnelle des
compositions possibles de cette phase. Chaque point de carré obtenu correspond à une
composition possible pour la phase. Chaque sommet est associé à l’un des composés
stoechiométrique définis par le modèle. La surface de référence d’énergie est représentée sur
la même figure-III-2-b.
La surface de référence est représentée par l’équation suivante :
ref
G 2ϕs = y 1A y c2 ° G Aϕ:C + y 1B yC2 ° GBϕ:C + y 1A y D2 ° G Aϕ:D + y 1B y D2 ° GBϕ:D (III-24)
(a) (b)
Figure-III-2: a- Représentation de l’ensemble des répartitions de site possibles pour
le modèle (A,B)a (C,D)b.
b- Enthalpie libre molaire de référence ref G d’un composé décrit par
(A,B)a (C,D)b en fonction de la répartition sur les sous-réseaux [2].
- 47 -

III-3-1-2- L’enthalpie libre idéal de mélange de Gibbs

L'énergie idéale de Gibbs de mélange peut être déduite à partir de l’équation (III-18) :
id  a
G2ϕs = RT  ( y 1A Lny 1A + y 1B Lny 1B ) +
b
( y C2 Lny C2 + y D2 Lny D2 )

(III-25)
a + b a+b 
III-3-1-3- L’enthalpie libre d’excès de mélange de Gibbs

En considérant que la solution est régulière dans chaque sous-réseaux, Alors
l’énergie d’excès de mélange de Gibbs peut s’écrit comme suit [30] :
ex
G2ϕs = y1A y1B L0A, B:∗ + yC2 yD2 L0∗:C , D (III-26)
où
L0A, B:∗ et L0∗:C , D représente les paramètres d’interaction de la solution régulière dans le sous-
réseaux qui sont indépendants de l’occupation de site des autres sous-réseau. On peut
introduiser le modèle sub-régulier en considérant que les interactions compositionnelles
dépendent de l’occupation de site des autres sous- réseaux.
Si ∗ = C , D, A, B ⇒
ex
G2ϕs = y 1A y 1B yC2 L°A, B:C + y 1A y 1B y D2 L°A, B:D + yC1 y 1D y A2 L°A,C:D + yC1 y1D y B2 L°B ,C:D (III-27)
avec
∗ : Représente le nom de n’importe quel élément existant dans les autres sous-réseaux,
lorsque ces derniers sont occupés par plusieurs éléments.
Finalement quelques fractions de site dépendent à ces paramètres sont ajoutées tel
que :
L°A, B:C = y 1A y 1B yC2 ∑ LνA ,B:C ( y 1A − y 1B )

ν
(III-28)
ν
L°A, B:D = y 1A y 1B y D2 ∑ LνA ,B:D ( y 1A − y 1B )

ν
(III-29)
ν
L°A:C , D = y 1A yC2 y D2 ∑ LνA;C ,D ( yC2 − y D2 )

ν
(III-30)
ν
L°B:C , D = y 1B yC2 y D2 ∑ LνB ;C ,D ( yC2 − y D2 )

ν
(III-31)
ν
Si on considère que les interactions dépendent de la température et de la

composition, un ordre très élevé d’interactions peut être représenté. C’est le cas similaire à
celui des paramètres d’interactions de Redlich-kister [2].
- 48 -

III-3-1- 4-La relation entre l’enthalpie atomique et celle de site

ϕ
La relation entre l’enthalpie libre molaire atomique 0Gma et l’enthalpie libre molaire de site
° ϕ
Gms est la suivante [30] :
°ϕ
° ϕ Gms
=
∑ N s (1 − yVas )
G ma (III-32)
s
ϕ ϕ
Si la description ne fait pas intervenir de lacune (0Gma =0Gms )
III-3-2-Cas particuliers
III-3-2-1-Cas d’une solution

Une solution métallique est décrite par un seul sous-réseau contenant les deux
éléments A et B symbolisé (A,B) . C’est le cas d’une solution de substitution.
Les deux éléments A et B sont mélangés dans le sous-réseau avec des fraction telles
que yA et yB =1. Dans ce cas particulier, les fractions yA et yB sont égales aux compostions xA
et xB . L’enthalpie de Gibbs de la solution est donnée par [30] :
G ϕm =refG ϕ +idG ϕ +exG ϕ (III-33)
où
ref
G ϕ =xAGA ϕ +xBGB ϕ (III-34)
ref
G ϕ : est le terme de référence correspondant à la somme pondérée de l’enthalpie de
Gibbs des éléments pur.
Le terme idéal correspond à une distribution aléatoire des atomes dans la solution et
le terme d’excès décrit à l’aide du polynôme de Redlich-kister, ces deux cas ont été discutés
dans le chapitre précédent.
III-3-2-2-Cas d’un composé défini
Dans le cas d’un composé stoechiométrique, chaque sous-réseaux est occupé par un
seul élément les termes Gid et Gex sont alors nuls et l’enthalpie libre molaire du composé
s’exprime simplement par :
GAaBb=refG+ ∆H f − T∆S f (III-35)
où
ref
G : est la somme pondéré des enthalpies libres des éléments purs dans l’état de référence
choisi. ∆H f et ∆S f sont respectivement l’enthalpie et l’entropie de formation.
- 49 -

III-3-2-3- Modèle à deux sous-réseaux identiques pour un système binaire

Dans le cas où A et B sont admis sur les deux sous-réseaux, une phase ordonnée du
type: (A,B)a (A,B)b est modélisée.
Les composés stoechiométriques définis par le modèle sont alors AaBb, BaAb, AaAb, et
BaBb. Les deux derniers correspondent aux corps purs dans la structure du composé AaBb.
BaAb correspond à l’antistructure de AaBb. Chacun de ces composés est reporté sur la
figure-III-3 [30].
Figure -III-3: Ensemble des réparations de site définies par le modèle (A,B)a (A,B)b[30].
La composition globale de la phase pourra décrire tout le système binaire A-B. Elle est liée à
la répartition sur les sites par les relations suivantes:
ay A' + by A'' ay B' + by B''
x ϕA = x ϕB = (III-36)
a+b a+b
Une composition globale de la phase pourra donc être obtenue par une infinité de répartition
de sites. Les fractions de sites occupées par l’élément A sont liées par la relation suivante:
a+b ϕ a '
y A'' = xA − yA (III-37)
b b
Ainsi, sur la figure-III-3, le lieu de composition constante x ϕA est-il représenté par la droite de
a 1 1
pente − passant par y "A = y 'A = x ϕA . Les droites correspondant aux compostions x ϕA = et
b 3 2
sont représentées sur la figure-III-3 dans le cas a=1 et b=2
La figure-III-4-a représente la nappe d’enthalpie libre molaire en fonction de y 'A' et y 'A pour
des valeurs particulières de l’enthalpie libre molaire des composés stoechiométriques définis
- 50 -

par le modèle lorsque les termes d’interactions sont nuls. Pour chaque composition, cette
nappe présente un minimum. Le lieu de ces minima est reporté sur la figure–III-3 dans le plan
des répartitions de site.
b
b
Figure -III-4: enthalpie libre molaire ref G ϕ + id G ϕ , en kj /mol, à 1000k, en fonction de la

répartition sur les sous-réseaux pour ( A, B )a ( A, B )b avec a=1, b=2
a-représention tridimensionnelle
b-projection des contours iso-enthalpie libre [30].
- 51 -

III-3-2-4-modèle à deux sous-réseau non- identique pour un système ternaire

Soit un composé ternaire dans lequel les éléments A et C partagent le même site alors
que l’élément B occupe seul le sien. Ce type de phase est obtenu lorsque les caractéristiques
des atomes A et C sont proches et différentes de celles de B. C’est le cas pour les composés
III-V ou les carbures M23C6 ou M7C3 par exemple. Une telle phase peut être schématisée de la
manière suivante: ( A, C )a Bb ou ( AyA , C yC )a Bb avec y A + yC = 1
La composition de cette phase est liée aux fractions de site par [30]:
a a b
x ϕA = yA , xCϕ = yC et x ϕB = (III-38)
a+b a+b a+b
L’ensemble des compositions décrit par ce modèle est représenté par

la figure -III-5 dans le cas a=1, b=2.
Figure -III-5: Ensemble des compositions décrit par le modèle (A,C) B.

L’enthalpie libre est la somme des termes suivants [30]. :
ref
G ϕ = y A ° G ϕA:B + y C ° GCϕ:B (III-39)
(III-40)
id
G ϕ = RT
a
[y A Ln( y A ) + y C Ln( y C )]
a +b
ex
G ϕ = y A y C LϕA,C :B (III-41)
L’enthalpie libre de référence est exprimée comme la somme pondérée des enthalpies
libres molaires des composés AaBb et CaBb de même structure que AaBb, respectivement
0
GϕA: B et 0GCϕ : B .
- 52 -

Le paramètre LϕA ,C : B apparaissant dans l’expression de l’enthalpie libre d’excès
représente l’interaction entre les éléments A et C dans le premier sous-réseau lorsque le

second est entièrement occupé par B [30].
III-3-2-5- modèle à deux sous-réseau non- identique pour un système binaire

Le composé AaBb dans lequel A est partiellement substitué par B peut être schématisé
de la manière suivante: ( A, B )a Bb ou (A1− yB , B yB ) Bb
a
C’est un cas particulier du précédent qui décrit un composé binaire. Ce modèle est
principalement utilisé pour décrire les phases de laves
La composition globale de cette phase ϕ est liée aux fractions de site de la manière
suivante:
a a b
x ϕA = yA et x Bϕ = yB + (III-42)
a+b a+b a+b
Les deux composés stoechiométriques définis par ce modèle sont AaBb et BaBb. Ce
dernier correspond à l’élément B pur dans la structure de même symétrie que celle du
composé AaBb. Toutefois, le fait d’introduire les mêmes atomes sur tous les sites entraîne une
dégénérescence de cette structure; ainsi, le composé BaBb correspond à un état de l’élément B
pur, , stable ou métastable.
La figure -III-6 représente l’enthalpie libre molaire d’un composé ainsi décrit référée
aux états α pour A et pour β à la température considérée. Dans le cas représenté par la
figure-III-6, le terme d’interaction LA,B :B est nul. Les différents cas représentés correspondent
à des valeurs différentes de l’enthalpie libre molaire de formation du composé de AaBb,
0
GϕA: B [30].
La surface de référence définie par ref

Gϕ se réduit ici à une simple droite liant 0G ϕA: B
et G Bϕ:B .
Pour le cas (1), cette droite est confondue avec l’axe.
Pour les cas (2) et (3), elle est tracée en pointillés.
La courbe (1) - pour ∆GϕA:B = 0 - représente id Gϕ . Elle traduit simplement le mélange
aléatoire des deux éléments dans le premier sous-réseau.
Les courbes (2) et (3) représentent ref
Gϕ + idGϕ pour 0∆Gϕ = —5 et —7kJ/mol. La
comparaison des trois courbes montre que lorsque l’enthalpie libre molaire de formation du
composé AaBb devient plus négative, l’enthalpie libre molaire totale diminue sur tout le
- 53 -

domaine de composition; c’est l’effet dû au terme de référence. De plus, son minimum se

rapproche de la composition correspondant à AaBb [30].
Figure -III-6: Enthalpie libre molaire d’un composé décrit par (A, B)1B
Sur cette figure, les deux composés stoechiométriques AaBb et BaBb définissent les
bornes du modèle. [30].
III-3-2-6- modèle à deux sous-réseaux présentant de lacune pour un système binaire

Toujours dans le cas d’une phase du système binaire A-B, considérons le cas où le site
de A n’est pas entièrement rempli. Ceci peut être schématisé de la manière suivante:
( A, Va )a Bb ou (AyA , Va yVa ) avec y A + y B = 1
Les composés stoechiométriques définis par ce modèle sont AaBb et VaaBb. Encore
une fois, le composé VaaBb correspond à l’élément B pur dans une structure dégénérée β. La
composition globale de la phase ϕ est reliée à l’occupation du premier sous-réseau par les
relations suivantes:
a b
x ϕA = yA et x ϕB = (III-43)
ay A + b ay A + b
La description faisant intervenir des lacunes, l’enthalpie libre molaire atomique,
0 ϕ ϕ
Gma et l’enthalpie libre molaire de site, 0Gms , sont différentes. Les trois termes d’enthalpie
libres deviennent :
- 54 -

° ϕ
° ϕ G ms (III-44)
G ma =
ay A + b
avec (III-45)
ϕ ϕ ϕ ϕ
G ms = ref
G + G + G
id iex
où
(III-46)
ref
G ϕ = y A ° G ϕA: B + y Va ° G Va
ϕ
:B
id
G ϕ = RT
a
[ y A Ln ( y A ) + y Va Ln ( y Va )] (III-47)
a+b
(III-48)
ex
Gϕ = y A y Va LϕA ,Va : B
III-4-généralisation du modèle CEF

On peut généralisé Les équations (III-24) (III-25) (III-26) pour un sous-résau
multiple. Pour cela, on définit une matrice de fractions de sites comme suit [32] :
 y11 y 12 y31 ... y 1C 

 2 2 2 2
 y1 y 2 y3 ... y C 
y = .  (III-49)
 
. 
 i i i i 
 y1 y 2 y 3 ... y C 
Où Chaque rangée représente un sous-résau et chaque colonne représente un

composant. L’occupation I peut etre utilisé pour définir les composant présentant dans les
différents sous-réseaux. Si il existe un seul composant dans chaque sous-réseau, l’occupation
est de l’ordre zero et noté par I0. I1représente l’occupation du premier ordre. c-à-dire que tous
les sous-réseaux sont occupés par un seul atome sauf un seul qui est occupé par deux
éléments différents. I2 indique que tous les sous-réseaux sont occupés par un seul atome sauf
deux sous-résaux qui sont occupés par deux éléments différents. Les paramètres d'interaction
s'appellent les paramètres réciproques. Le choix plus général de l’occupation I d'un ordre plus
supérieur peut être défini comme suit IZ avec ces occupations ne doit contenir aucun
composant plus d'une fois dans le même sous-réseau. L’enthalpie libre de référence de Gibbs
peut s’écrit comme suit [32] :
ref
Gϕ = ∑ PI (Y ).GI° (III-50)
I
- 55 -

G I° : représente l’enthalpie de Gibbs du composé définit par I
PI (Y ) : représente le produit des fractions de sites de la matrice correspondante.

L’enthalpie idéale de Gibbs est donnée par l’équation (III-18).
L’enthalpie libre d’excès est donnée par :
ex
G ϕ = ∑ PI 1 (Y ).LI 1 (III-51)
I1
Si la solution est considérée régulière dans les différents sous-réseaux. Pour un ordre
d’interaction très élevé l’équation (III-48) deviendra :
ex
G ϕ = ∑∑ PIZ (Y ).LIZ (III-52)
Z >0 IZ
L’enthalpie molaire de Gibbs totale s’écrit [32] :
ϕ
Gms = ∑ PI 0 (Y ).GI°0 + RT ∑ N s ∑ yi Ln ( yi ) + ∑∑ PIZ (Y ).LIZ (III-
I0 s i Z >0 IZ
53)
III-5-La transformation ordre ↔ désordre
La transformation de type α ' ↔ α concerne l’équilibre entre une phase ordonnée à

longue distance à basse température et désordonnée à longue distance à haute température : on
dit qu’une phase est ordonnée à longue distance si on peut distinguer sur son réseau cristallin
deux ou plusieurs «sous-réseaux» chacun préférentiellement occupé par un type d’atomes.
C’est le cas des composés définis, Il peut se produire que l’occupation des deux sous-réseaux
devienne identique en composition vers les hautes températures, lors de l’accroissement du
désordre atomique. On ne peut plus alors distinguer les divers sous-réseaux et les deux phases
sont identiques. Si le désordre s’installe progressivement sur la phase α ' jusqu’à l’obtention
de la phaseα la frontière de séparation des deux phases sera constituée d’une ligne unique
appelée ligne critique. C’est le cas par exemple de la transition α ' ↔ α dans les alliages
(Al,Fe) ou (Fe,Si) . La transition sur la ligne critique se fait alors sans enthalpie latente de
transition ni variation de volume, elle est appelée transition thermodynamique du second
ordre ( A2 ↔ B 2 ). L’étude détaillée de ces transitions critiques sort du cadre de notre étude
du point de vue diagrammatique, leur particularité est l’existence de la ligne critique à la place
du fuseau de lignes conjuguées. La transition β ' ↔ β des laitons (Cu,Zn) est également de ce
type, elle joue d’ailleurs un rôle technologique très important dans la pratique : les tuyaux
- 56 -

d’échangeurs de chaleur sont souvent réalisés en laiton β ' , phase très fragile à l’état ordonné
pour laminer ce matériau il est indispensable de désordonner cette phase, elle devient alors
très malléable à haute température (600 °C). À l’inverse dans les laitons dits de décolletage,
qui peuvent s’usiner à très grande vitesse de coupe, c’est précisément la fragilité de cette
même phase ordonnée β ' qui permet la rupture rapide du copeau sous l’outil de coupe.
Pour d’autres structures ordonnées, la transition ordre ↔ désordre s’accompagne
d’une enthalpie latente de transformation elle est alors du premier ordre thermodynamique,
comme la transitions de phases, telle la fusion.. Par exemple, les congruences signalées dans
les alliages AuCu et AuCu3 sont des transformations ordre ↔ désordre du premier ordre. Les
deux congruences observées sur la même composition CoFe sont de nature très différentes : à
730°C il s’agit d’une transition ordre ↔ désordre congruente du premier ordre sur un réseau
α non compact tandis qu’à 990 °C il s’agit de la transition c.c. ↔ c.f.c. [12]
III-6- Modélisation de la réaction ordre-désordre

Souvent on trouve des phases différentes, comme la phase désordonnée (A2) et la
phase ordonnée (B2) présentées sur la figure-III-7, mais ils ont la même structure cristalline.
L’enthalpie libre de Gibbs de la phase ordonnée peut être modéliser en tenant compte la
relation entre la phase ordonnée et la phase désordonnée. En particulier, lorsque la réaction
ordre-desordre du second ordre se produit entre les deux phases, il est nécessaire d’avoir une
description pour l’enthalpie de Gibbs des deux phases afin de décrire la transition graduelle
d'une phase à une autre. L’expression de l’enthalpie de Gibbs pour une phase ordonnée est
définit comme suit
Gmsl = G mdes (xi ) + ∆Gmord ( yis ) (III-54)

Le premier terme représente la contribution de la phase désordonnée et le second
terme représente la contribution d’ordre. La phase est désordonnée, c-à-d que lorsque les
fractions du réseau sont égales au fraction molaire, l’enthalpie de Gibbs doit être égale à celle
de la phase désordonnée .Ceci est assuré en définissant l’enthalpie molaire comme suit :
∆Gmord ( yis ) = Gmsl ( y is ) − G msl ( yis = xi ) (III-55)
L’enthalpie libre de la phase désordonnée est représentée par le modèle de la solution

substitutionnelle comme nous l’avons déjà vu dans le chapitre précédent. [32]
- 57 -

a b
Figure-III-7 : Représentation schématique des structures A2 et B2.
(a) Tous les sites sont équivalents.
(b) Les sommets de la maille sont majoritairement occupés par les atomes A, le
centre par les atomes B. [32]
III-7-Modélisation de l’ordre à courte distance

Le formalisme général des sous-réseaux décrie les phase, qui présentent un ordre à
longue distance, par différents sous-réseaux. Cependant, sur chaque sous-réseau, une solution
régulière est encore assumée. Ceci signifie que l’ordre à courte distance (OCD) n'est pas pris
en considération. Particulièrement pour les phases de structures cfc, l'ordre à courte distance
peut avoir une contribution importante à toute l'énergie de Gibbs du système. Par conséquent,
le formalisme de sous-résaux ne peut pas décrire d’une manière satisfaisante les deux limites
de phase et les propriétés thermodynamiques. Sundman et al. [33] ont rencontré ce problème
pour le système Au-Cu. Ils ont constaté que dans l'approximation quasi-chimique ∆GmSRO pour
de petites valeurs de paramètres d'ordre à court distance peut s’écrit comme suit:
∆GmSRO = − yi1 y 1j y i2 y 2j
(∆ G)
ref 2
(III-56)
zRT
Pour un système binaire qui est modélisé avec deux sous-résau, Ζ représente le
nombre de coordinance et ∆ref G égal à :
∆ref G = ref Gi: j + ref G j:i − ref Gi::i − ref G j: j (III-57)
L’équation (III-57) indique que l’enthalpie libre d’une solution présentant un ordre à
courte distance (OCD) dépend de la composition comme le cas des paramètres d’interaction
[32].
- 58 -

III-8-Application du modèle aux différentes phases
Dans cette partie on va présenter la modélisation de quelques phases importante pour

notre travail et aussi leurs classement suivant le modèle utilisé pour les décrire (Tableau-III-1)
[34-52].
III-8-1- Modélisation des phases de structure B2
La structure cubique centrée B2 de type CsC1 peut être considérée comme

l’interpénétration de deux sous-réseaux cubiques simples. Elle peut également être décrite
comme une structure ordonnée de la structure cubique centrée A2. La figure-III-7 représente
schématiquement les structures A2 et B2.
Les phases B2 que nous étudierons présentent des défauts ponctuels de type triple
défaut.
-présence d’atomes B sur le site des atomes A (B en antistructure) ;
-présence d’atomes A sur le site des atomes B (A en antistructure) ;
-présence de lacunes sur le site des atomes A ;
-présence de lacunes sur le site des atomes B.
La présence de défauts interstitiels est exclue par la faible différence de rayon
atomique des éléments conduisant à la formation de ces phases.
Nous les décrirons donc grâce au modèle en sous-réseaux selon le schéma suivant:.
Les phases B2 présentant des défauts d’antistructure peuvent être modélisées par (A,B) (A,B)
et pour le défauts triple par le modèle(A,B) (B,Va). Afin de pouvoir combiner les descriptions
de systèmes présentant des types de défauts différents nous introduirons également une
description symétrique équivalente à celle introduite par Krachler [30]:
(A y 'A
, B y ' ,Va y '
B Va
) (A y 'A'
, B y '' ,Va y ''
B Va
)
La composition d’une phase ainsi décrite est liée aux fractions de sites par les
relations suivantes:
y 'A + y 'A' y B' + y B''

B2
= B2
=
(1 − yva' ) + (1 − yva'' ) (1 − yva' ) + (1 − yva'' )
x A x B (III-58)
- 59 -

La figure -III-8 présente l’ensemble des répartitions de site décrit par ce modèle. Elle
représente en particulier les lieux de composition constante x A = 0.3,0.4,0.5,0.6 et 0 .7
Les lieux de composition constante sont définis par des équations de droite du type:
y 'A + x ϕA yVa
''
− 2 x ϕA = 0 (III-59)
Figure -III-8 : Ensemble des répartitions de site décrit par le modèle

(A,B)a (B,Va)b. [30]
L’enthalpie libre molaire, G, d’une phase décrite par le modèle Krachler est définie
par l’équation suivante [53]:
ref
G B 2 + id G B 2 + ex G B 2
B2
Gma = (III-60)
1
2
(
1 − y va' ) + (1 − y va'' )
1
2
avec
G B 2 = ∑∑ yi' y 'j' ° GiB: j2
ref (III-61)
i j
( ) ( )
i j
= RT  y i' Ln ( y i' ) + yi'' Ln ( yi" )

id B2 1 1  (III-62)
G
2 2 
ex
G B2 = ∑∑∑ y i y j y k Li , j:k + ∑∑∑ y i" y 'j y k' LBk:2i , j
' ' " B2
(III-63)
i j >i k i j >i k
+ ∑∑∑ yi' y 'j y k" yl" LBi , 2j:k ,l

i j >i k
- 60 -

III-8-2- Modélisation de la phase σ et χ

La plus part des composées intermétalliques dans les systèmes binaires peuvent
présenter une solubilité considérable et ils peuvent avoir plus de deux sous-réseaux
cristallographique, c’est le cas des composé σ , µ , χ et les phases de laves. Le choix du
modèle pour décrire ces phases est considéré selon la structure interne de ces composés .
Dans le cas de σ (Fig.III-9) qui contient cinq sous-réseaux cristallographiques distincts qui
sont occupés par des atomes suivant le rapport 2:4:8:8:8. Le nombre total d’atomes par cellule
unitaire est alors 30 atomes. Le mélange des composants sur tous les sous-réseaux produit
un nombre très élevé de variables, Ce qui rend un calcul Calphad pour un système multi-
composant presque impossible. Même la Modélisation d’une phase binaire sera extrêmement
difficile et quelques simplifications sont nécessaires. Ceci peut être réalisé pour la phase σ en
considérant deux types de composants de taille différente. L’occupation des sous-réseaux
peut être simplifier ; pour la première possibilité, σ est considérée comme suit (A,
B)16(A,B)10(B)4 où le premier et le troisième sous-réseau ont été intégrés entre eux , même
chose pour le quatrième et cinquième sous- réseau, l’occupation de sites par les atomes A,B
est considéré égale [2]. On peut utilisé la deuxième possibilité (A,B)16(A,B)10(B)4 où le
mélange s’effectue seulement entre le quatrième et le cinquième sous-réseau. Cette dernière
occupation survient lorsque la solubilité de B dans le 2emme sous-réseau est minimale.
Les modèles utilisés pour décrire ce genre de phase sont, soit le modèle
(A,B)16(A,B)10(B)4 ou le modèle (A,B)16(A)10(B)4 en prenant l’exemple du système Co-Mo ou
la phase σ est décrite par (Co,Mo)16(Co)10(Mo)4. L’enthalpie libre pour cette phase est
donnée par la forme suivante :
Gmσ = ° GCo
σ ° σ
[
:Co:Mo + G Mo:Co:Mo + RT 16( y Co Lny Co + y Mo Lny Mo )
1 1 1 1
]
3  ν 
(y − y 1Mo ) 
(III-
+ yCo
1
y 1Mo yCo
2
y Mo  ∑ LνCo ,Mo:Co:Mo 1
Co
 ν 
64)
2
Si y Co = y Mo
3
l’enthalpie libre est celle du modèle de Redlich-Kister . Ce mélange est
simplifié à deux composants occupant un seul sous-réseau [2].
- 61 -

(a)
(b)
Figure-III-9 : a- Structure cristalline de la phase σ à la composition A2BL’atome A
est en gris foncé et l’atome B est en gris clair.
b- Projection des atomes sur le plan x-y [54].
III-8-3- Modélisation de la phase AB2 (phases de laves)

Les composés intermétalliques de type AB2 sont des phases de Laves qui cristallisent
dans une structure cubique ou hexagonale selon le rapport des rayons atomiques et des
concentrations en électrons [55]. Les phases de Laves forment une très grande famille
structurale et de nombreuses substitutions par des éléments réduisant ou augmentant le
volume de maille sont possibles. Le site A est occupé en général par un élément du groupe
IVA (Ti, Zr ou une terre rare) et le site B est occupé par un élément de transition de la
première rangée (i.e. du vanadium au cuivre).
Les composes de stoechiometrie AB2 cristallisent dans l’une des trois structures type :
MgCu2 (structure cubique ou C15), MgZn2 (hexagonale ou C14) ou MgNi2 (structure
hexagonale ou C36) [56,57]. Toutes ces phases sont caractérisées par des polyèdres de
coordination identiques pour l’atome A appelés polyèdres de Friauf. Les trois structures
- 62 -

peuvent être décrites comme un empilement de ces polyèdres. On va décrire les trois
structures types de phases de Laves.
Ø Structure type MgCu2 ou C15

Les atomes B s’arrangent en tétraèdres, les tétraèdres étant joints entre eux par les
sommets. Les atomes A occupent alors les espaces libres entre les tétraèdres formant entre eux
un réseau de type diamant. Le réseau est cubique a faces centrée et il y a 8 motifs AB2 par
maille. La représentation de cette structure est présentée sur la figure-III-10.
Figure-III-10 : structure de type MgCu2 ou C15 [58]
Ø Structure type MgZn2 ou C14

Ce réseau est hexagonal. Les atomes B s’arrangent toujours en tétraèdres mais ces
tétraèdres sont alternativement joints par les bases et par les sommets. Il y a cette fois 4 motifs
par maille (Fig.III-11).
Figure-III-11: Structure de type MgZn2 ou C14 [58].
- 63 -

Ø Structure type MgNi2 ou C36

Cette structure est une combinaison des deux précédentes. Le réseau est hexagonal
avec une maille doublée selon c par rapport a la structure C14. On retrouve les chaînes des
tétraèdres qui existent dans la structure C14 mais ces chaînes ne sont pas infinies et sont
décalés les unes par rapport aux autres (Fig.III-12) [58].
Figure-III-12 : structure de type MgNi2 ou C36 [58].
Généralement les modèles utilisés pour décrire ce genre de phases sont représentés
par trois types de formule:
(A,B)2(A,B), (A,B) 2B, (A,B)6 (A,B)4 A2
Les trois termes de l’enthalpie libre molaires sont donc définie comme suit :
Pour la phaseϕ décrite par (A,B)2(A,B):
ref
G ϕ = ∑∑ yi' y "j °Giϕ: j (III-65)
i j
2 1 
id
G ϕ = RT  ∑ yi' Lnyi' + ∑ y i" Lnyi"  (III-66)
3 i 3 i 
(III-67)
G = y A y B ∑ yi L A, B:i + y A y B ∑ y i Li: A, B
ex ϕ ' ' " ϕ " " ' ϕ
i i
avec
y 'A et y "A respectivement les fractions de site du premier et du second seous-réseau occupés
par l’élément A dans la phase ϕ
- 64 -

Pour la phase C15 décrite par (A,B)2A :

ref
G C15 = ∑ yi ° GiC:B15 (III-68)
i
2
id
G C 15 = RT ∑ y i Lnyi
3 i (III-69)
ex
G C15 = y A y B LCA15
, B:B (III-70)
avec
y A la fraction du premier sous-réseau occupé par l’élément A dans la phase C15
Pour la phase C14 décrite par (A,B)6 (A,B)4 A2:
ref
G 14 = ∑∑ yi' y "j °Gi14: j: A
i j (III-71)
6 4 
id
G 14 = RT  ∑ yi' Lnyi' + ∑ yi" Lnyi" 
12 i 12 i  (III-72)
G = y A y B ∑ yi L A, B:i: A + y A y B ∑ yi Li: A, B: A
ex C14 ' ' " C 14 " " ' C 14
i i (III-73)
avec
y 'A et y "A respectivement les fractions de site du premier et du second sous-réseau
occupés par l’élément A dans la phase C14 [30].
III-8-4 -modélisation de la phase CFC ( L12, A1, L10)

La structure de cette phase peut être décrite par une mise en ordre sur les différents
sites de la structure cubique à faces centrées. Les figures.III-13 a et b présentent une maille
élémentaire pour les structures Al et L12. Dans la structure Al désordonnée, tous les sites sont
équivalents. Dans la structure L12, la répartition des éléments sur les sites aux sommets de la
maille est différente de celle au centre des faces.
(a) (b) (c)
(a) A1 (b)L12 (c)L10

Figure-III-13 : Représentation schématique des structures des phases A1, L10 et L12
- 65 -

Toutefois, la structure c.f.c. peut également s’ordonner en L10. On obtient alors des
plans alternés où la répartition des éléments est différente. La figure-III-13-c présente la
maille élémentaire correspondant à cette structure.
Ansara et al. [59] ont été traité la phase L12 dans le cas du système binaire AlNi par
le modèle en sous-réseaux où les deux éléments peuvent se trouver sur deux sous-réseaux:
(Al, Ni)3 (Al, Ni) .
On peut également décrire la phase L10 avec le modèle à deux sous-réseaux mais
avec un rapport entre le nombre de sites attribué à chaque sous-réseau égal à 1.
(A,B,...) (A,B,...) ou (A,B,...)2 (A,B,...)2[30].
Les phases de structure L12 sont généralement modélisées par deux et quatre sous-
réseaux. L’avantage d’utiliser deux sous-réseaux pour que la modélisation être plus simple.
Alors que en réalité , on utilise le modèle à quatre sous-réseau qui a été extrapolé au
système plus compliqués (exemple L21)d’une part et d’autre part, Certains systèmes font
apparaître les trois phases L12 A1 et L10 (Au-Cu par exemple) [2]. Il est donc intéressant de
décrire ces trois phases avec un même ensemble de paramètres thermodynamiques. Ceci peut
être réalisé en utilisant un modèle à quatre sous-réseaux auxquels on affecte un même nombre
de sites. (A,B,...) (A,B,...) (A,B,...) (A,B,...).
Expression de l’enthalpie libre de la phase L12 décrite par deux sous-réseaux
Pour une structure ordonnée notée (ord) de la phase désordonnée (des) qui décrite par
(A,B,.. )a (A,B,..)b où n éléments identiques occupent les deux sous-réseaux . Nous avons déjà
exprimé l’enthalpie libre d’une phase ainsi décrite par le modèle en sous-réseaux lors de
l’étude de la modélisation des composés intermétallique non-stoechiométrique. Rappelons
simplement :
Gmord = ref G ord + id G ord + ex G ord (III-74)
avec
G ord = ∑∑ yi' y "j °Giord
ref (III-75)
:j
i j
 a 
y i' Ln ( y i' ) + y i" Ln ( yi" )
b
id
G ord = RT  ∑ ∑ (III-76)
a+b i a+b i 
ex
G ord = ∑ ∑ yi' y 'j ∑ y k" Lord
i , j:k + ∑∑ y i y j ∑ y k Lk :i , j
" " ' ord (III-77)
i j >i K i j >i K
+ ∑∑∑∑ y i ' y 'j y k" y l" Lord

i , j : k ,l
i j >i K l >k
- 66 -

Pour que la fonction Gmord décrive les deux phases (ord) et (des), elle devra présenter
un minimum lorsque la phase (ord) correspond à la phase des c’est-à-dire lorsque y ' = y " pour
certaines compositions [30].

Expression de l’enthalpie libre de la phase L12 décrite par quatre sous-réseaux
Les phases ordonnées et désordonnées sont maintenant décrites de la manière
suivante:
(A,B,...) (A,B,...) (A,B,...) (A,B,...).
L’expression de l’enthalpie libre molaire fait intervenir de nombreux paramètres.
Toutefois, afin de pouvoir décrire la phase désordonnée à l’aide de ce modèle, chacun des
quatre sous-réseaux est considéré identique. Il en résulte que les enthalpies de formation des
composés de stoechiométrie identique doivent être égales. Ainsi:
° ° ° ord °
:B:C :D = G A:B:D:C = G A:C :B:D = G A:C :D:B = .......
G Aord ord ord
l’enthalpie libre est donnée par l’expression suivante:

Gmord = ref G ord + id G ord + ex G ord
avec
ref
G ord = ∑ ∑∑∑ yi1 y 2j y k3 yl4 °Giord
: j :k ;l (III-78)
i j k l
RT ∑∑ yis Ln ( yis )
4
1
id
G ord =
4 s= i
(III-79)
ex
G ord = +∑∑ ∑∑∑∑ y ' y i
'
j y k" y l" Lord
i , j : k ,l (III-80)
i j >i K l > k
4′′
+ ∑∑ ∑∑∑∑∑ y i
s
y sj y kr y lt y mu Lord
i , j:k :l :m
s =1 i j >i k l m
Expression de l’enthalpie libre de la phase L12 présentant une transformation

ordre-desordre
La phase c.f.c. est présente dans de nombreux systèmes sans présenter un ordre. Dans
ces systèmes, il serait inutile d’introduire une description comportant deux ou quatre sous-
réseaux pour une simple solution de substitution. Toutefois, il est souhaitable d’avoir une
description identique d’une même structure dans différents systèmes. Dans ce but, l’enthalpie
libre des phases Al et L12 comporte deux termes: équation (III-55) [30].
Il en résulte l’expression suivante pour l’enthalpie libre molaire totale de la phase
ϕ [32] :
( ) ( ) ϕ
Gmϕ = G A1 xiϕ + G L12 y is − G L12 yis = xiϕ ( ϕ
) (III-81)
- 67 -

Tableau-III-1 : Classement des phases selon le modèle de sous-réseaux [30,34-52].
La phase AB
Phase XB % Structure Sous-réseaux réference

cristalline
TiB 0.49-0.51 B27 (Ti)1 (B, Ti)1 [34]
RuTi 45- 52 B2 (Ru, Ti)0.5(Ru, Ti)0.5 [35]

49-55 B2
InPt 48 -51 - (In,Pt)1(In,Pt)1 [36] [37]
γ PdZr 46-51 B2 (Pd,Zr)0.5(Pd,Zr)0.5(Va)3 [44]

49-50
α TiPt 46-54 B19 (Pt, Ti)1(Pt, Ti)1 [36]
β TiPt 46-54 B2 (Pt, Ti)1(Pt, Ti)1 [36]
RuZr 48-52 B2 (Ru, Zr)0.5(Ru, Zr)0.5 [52]
LiZn 50-54 B32 (Li,Zn)05(Li,Zn)05 (Va)3. [36] [38]

AlPd 44-56 B2 (Al,Pd) 1 (Al, Pd)1 [40]
(Al, Pd) (Va, Pd)
TiAl 48-69.5 L10 (Al ,Ti)1 (Al,Ti )1 [36] [40]
AlIr 48-52 B2 (Al,Ir,Va)1\2(Al,Ir,Va)1\2
γNiSb 43-52 B81 SB1\3(Ni,Va)1\3(Va,Ni)1\3 [36] [41]
β CoSb 48-51 B81 SB1\3(Co,Va)1\3(Va,Co)1\3 [42]
AlFe 23.3-55 B2 (Fe,Al,Va)0.5(Fe,Al,Va)0.5 [36] [43]

AlTi 48-69.5 L10 (Al ,Ti)0.5 (Al,Ti )0.5 [30]
AlCo 46-78 B2 (Al,Co,Va) (Al,Co,Va) [30]
(Al,Co) (Co,Va) [30]
β AlNi 42-69.2 B2 (Al,Ni,Va) (Al,Ni,Va) [30]
(Al,Ni) (Co,Ni) [30]
RuZr 49-50 B2 (Ru, Zri)0.5(Ru, Zr)0.5(Va) 3 [52]
HfRu - B2 (Ru, Hf)0.5(Ru, Hf)0.5(Va) 3 [52]
MoPt 41-50 B19 (Mo, Pt)0.5(Mo, Pt)0.5 [45]
NiTi 44-51 B2 (Ni,Ti)(Ni,Va) [30]
(Ti,Ni,Va) (Ti,Ni,Va)
La phase AB2, A2B
TiB2 0.67-0.69 C32 (B, Ti)1 (B, Ti)2 [34]

AlCr2 65.5-71.4 C11b (Al, Cr) 1 (Al, Cr) 2 [36] [46]
AlPd2 65-73 C23 (Al,Pd) 1 (Al, Pd)2 [36] [39]
PdZr2 66-66.7 - PdZr(Pd, Zr) [44]
η, TiAl2 65-68 - (Al,Ti)2 (Al,Ti)1 [36] [40]
- 68 -

Tableau III-1 (Suite)

Pd2Zr 33.3-37 - PdZr(Pd, Zr) [44]
Fe2Hf - - (Fe)0.6667 (Hf,Fe)0.3333 [30]
FeHf2 65-67 - (Fe,He) 0.3333 (Hf ) 0.6667 [36]
Ru2Zr 32-34 C15 (Ru, Zr)0.6667(Ru, Zr)0.3333 [52]
Laves_C1 25-34 C15 (Co, Zr) 2 (Co, Zr)1 [47]
5( )
CoZr2 26-33 - (Co, Zr) 1 (Zr) 2 [47]
Ta2Al 39-49 D8b Ta0.533Al0.333Ta0.133 [30]
Cr2Ta 33-34 C15 (Cr,Ta)2Ta [30]
31-38 (Cr,Ta)2(Cr,Ta)
C14 (Cr,Ta)6(Cr,Ta)4Cr2 [30]
(Cr,Ta)2(Cr,Ta)
MoPt2 65-68 - (Mo,Pt) 0.3333 (Mo,Pt ) 0.6667 [45]
α Cr2Zr 31-36 C15 (Cr,Zr)2(Cr,Zr)1 [30]
β Cr2Zr 30-40 C14 (Cr,Zr)6Cr,Zr)4Cr2 [30]
γ Cr2Zr 31.5-38 C36 (Cr,Zr)6Cr,Zr)4Cr2 [30]
La phase AB3, A3B

CoZr3 ( ) 74.2-75.3 E1a (Co, Zr) 1 (Co, Zr) 1 (Zr) 2 [47]
InPt3 71-76 L12 (In,Pt)025(In,Pt)075(Va)1 [36] [37]

(In,Pt)025(In,Pt)075
Pd3Sb 25-27 B32 (Pd, Sb) 3Sb [48]
Cr3Pt 17-23 A15 (Cr, Pt) 3 (Cr, Pt) [49]
βNi3Ge 22.5-25 L12 (Ni)0.75(Ge, Ni)0.25 [36]

Ti3Pt 22-29 A15 (Pt, Ti)3(Pt, Ti)1 [36]
TiPt3(γ) 75-81 L12 (Ti, Pt)Pt3 [36]
Pd3Zr 17-27 Hcp-Do24 (Pd,Zr)0.75(Pd,Zr)0.25(V a)0.5 [44]
Cr3Ge 20-25.1 A15 (Cr,Ge)0.75(Cr,Ge)0.25 [36] [50]
Ti3Al 22-39 Do19 (Al,Ti)3 (Al,Ti)1 [40]

AlIr3 26.5-28 L12 (Al,Ir,)1/4(Al,Ir,)1/4(Al,Ir,)1/4(Al,Ir,)1/4 [30]
βNi3Sb 23-25 Do3 (Sb)1/4(Ni,V a)1/2(Ni,V a)1/4 [36] [41]
δNi3Sb 23.7-25.1 (Ni)3/4(Ni,Sb)1/4 [41]
AlPt3 67-77 L12 (Al,Pt)0 ,25(Al,Pt)0 ,25(Al,Pt)0 ,25(Al,Pt)0 ,25 [36]
FeAl3 74.5-76.6 Do3 (Al)0.6275 (Fe)0.235(Al ,Va)0.1375 [36][43]
ηNi3Ti 20-25 Do21 (Ni,Ti)3(Ni,Ti)1 [30]
AlNi3 73-76 L12 (Al,Ni)3(Al,Ni)1 [30]
Al 3(Al,Ni)1
Mo 3Pt 38-51 Do19 (Mo, Pt) 3Pt [45]
Mo3Si 24-25 A15 Mo0.75(Mo,Si)0.25 [45]
- 69 -


Divers phases de types AmBn
La phase A2B3 et A3B2
γ Ni3Ge2 41-43 B81 (Ge,Ni)0.6 (Ge)0.4 [36]

βLi2Zn3 50-67 - (Li,Zn)2(Li,Zn)3 [38]
Al3Pd2 38-41.5 D513 (Al,Pd) 3 (Al, Pd)2 [36] [39]
αLi2 Zn3 58-60 - Li2(Li,Zn)3 [36] [38]
β In2Pt3 58-61 B81 (In,Pt)2(In ,Pt)3 [36] [37]
Al3Ni2 36.8-41.5 D513 Al3(Al ,Ni)2(Al ,Ni) [30]
Al3(Al ,Ni)2(Al ,Va)
La phase A5B3
αCr5Ge3 37.41-38.2 D8m (Cr,Ge)0.625(Cr,Ge)0,375 [36] [50]
βCr5Ge3 37.41-38.2 D8m (Cr,Ge)0.625(Cr,Ge)0,375 [36] [50]
εNi5Ge3 35-44 B81 (Ge,Ni)5(Ge,Ni)3 [36]

Pd5Sb3 34.3-38.5 B81 (Pd,Sb)5(Pd,Sb) 3 [48]
Mo5Si3 37-40 D8m (Mo,Si)0,625 (Mo,Si) 0,375 [45]
La phase A9B13
In9Pt13 55-59 - In9(In,Pt)13 [36] [37]

La phase A19B12
Ni19Ge12 38-41 - (Ge,Ni)19(Ge,Ni)12 [36]

La phase A2B5
θNi5Sb2 25.2-28.8 - (Ni)5/7(Ni,Sb)2/7 [36] [41]

β Li2 Zn5 71-92 - (Li,Zn)2Zn5 [36] [38]
αLi2 Zn5 70-73 - (Li,Zn)2Zn5 [36] [38]
Al2Pd5 70.7-71.7 - (Al)2(Al, Pd)5 [36] [39]

La phase A5B6
In5Pt6 51-54.5 - (In,Pt)5(In,Pt)6 [37]
La phase AB4
β LiZn4 50-54 B32 (Li,Zn)02(Li,Zn)08 (Va)0.5. [36] [38]
α LiZn4 77-85 - (Li,Zn)1(Li,Zn)4 [38]
Al4Cr 18.5-20 - Al4(AL,Cr)1 [46]
- 70 -


La phase σ et χ
σ 55-70 D8b Re10(Mo,Re)4(Mo,Re)16 [36] [51]

σ (Co,Mo)16(Co)10(Mo)4 [2]
(Mo,Re)20 (Mo,Re)10 [51]
χ 75.6-78 A15 (Mo,Re)10(Mo,Re)24Re24 [36] [51]
Les phases de la structure L12 et L10
TiAl 48-69.5 L10 (Al, Ti)1 (Al,Ti)1 [40]

InPt3 71-76 L12 (In,Pt)025(In,Pt)075(Va)1 (In,Pt)025(In,Pt)075 [36] [37]
AlIr 48-52 L10 (Al,Ir,)1/4(Al,Ir,)1/4(Al,Ir,)1/4(Al,Ir,)1/4 [36] [30]
Al3Ir - L12 (Al,Ir,)1/4(Al,Ir,)1/4(Al,Ir,)1/4(Al,Ir,)1/4 [30]
AlIr3 - L12 (Al,Ir,)1/4(Al,Ir,)1/4(Al,Ir,)1/4(Al,Ir,)1/4 [30]
AlPt3 67-77 L12 (Al,Pt)0 ,25(Al,Pt)0 ,25(Al,Pt)0 ,25(Al,Pt)0 ,25 [36]
III-9 - conclusion
L'énergie libre de Gibbs d’une phase modélisé par le formalisme CEF, en tenant en
compte l'ordre à longue et courte distance, peut être écrite comme suit :
tot
Gm = Gmdis (xi ) + ∆G mord ( y is ) (III-82)
Gmdis (xi ) : S’écrit comme dans le cas d’une solution non-régulière en utilisant
l'équation de Redlich-Kister.
∆Gmord : S’écrit en utilisant le modèle de sous-réseau.
∆Gmord ( y is )= sl Gm ( yis )− sl G m ( yis = xi ) (III-83)
et
sl
Gm = ∑ PI 0 (Y ).G I0 + RT ∑ N s ∑ y i Ln( y i ) + ∑ PI 1 (Y ).LI 1 + ∑ PI 2 (Y ).LI 2 (III-84)
s i
Où
∑ P (Y ).L
I1
I1
: représente l’interaction des deux éléments présentant dans le même sous-
réseau lorsque les autres sous-réseaux présentant un seul élément.

Et ∑ P (Y ).L
I2
I2
représente les termes d’interactions réciproques en tenant
compte l’ordre à courte distance.

Généralement la structure considérée pour les phases de type AB est représentée par B2 .De
nombreux composés intermétalliques de composition équiatomique présentent cette structure
(FeAl , RuTi, InPt, αTiPt, RuZr, AlPd,…). Leur domaine de stabilité est généralement étendu
- 71 -

en composition. Les déviations à la composition stoechiométrique d’un de ces composés de

structure de B2 sont possibles grâce à la formation de quatre types de défauts ponctuels.
Les phases B2 qu’on a étudiées présentent des défauts ponctuels de type triple
défaut, car il nécessite la présence de trois défauts à la stoechiométrie : un atome en
antistricture et deux lacunes sur l’autre sous-réseau.
Dans notre étude, on a remarqué que la phase équiatomique AB de structure B2 a été
décrite par deux formes de sous-réseaux.
1. la première forme, lorsque la phase ne présente pas de lacune, on peut la
schématiser par la formule suivante : (A,B)m(A,B) m. La structure B19 et L10 ont
été décrites par la même formule.
2. La deuxième forme, lorsque la phase présente de lacunes, on peut la décrire par
trois formes :
v soit en utilisant deux sous-réseaux de la forme (A,B)(A,Va), lorsque les lacune
existent sur le site de B. Ou bien de la forme (A,B)(Va,B), lorsque les lacune
présentent sur le site A.
v soit en utilisant trois sous-réseaux de la forme (A,B)0.5(A,B)0.5 Va3 . La structure
B32 a été décrite aussi par la même formule.
v Soit en utilisant deux sous-résaux de schéma (A,B,Va)(A,B,Va), lorsque les
lacune existent sur les sites A et B(défaut triple).
On a remarqué que le modèle utilisé pour décrire les phases de laves de la
forme AB2 et A2B, notamment les phases C14, C15 et C36, est presque le même selon le
schéma suivant (A,B)1(A,B)2 dans tous les cas qu’on a vue dans le tableau. Toute fois, on peut
aussi modéliser les phases C14 et C36 par (A,B)6(A,B)4A2 et pour la phase C15 par (A,B)2B,
lorsque l’un de constituant ne présente pas assez de solubilité dans l’autre. Notre étude pour
ce type de phases présente des cas exceptionnels comme les phases AlPd2, Ta2Al, PdZr2,
Pd2Zr et Ta2Al qui ont été modélisées différemment aux autres composés tableau (III-1).
Généralement les modèles utilisés pour décrire la phases de type A3B et AB3
sont respectivement de schéma (A,B) 3(A,B) et (A,B) A3. Il y a encore des exceptions dans ce
cas pour les phases : CoZr2, InPt 3, Pd3Zr, AlIr3, Ni3Sb, AlPt3 et FeAl3 tableau (III-1). On
s’intéresse en particulier au modèle (A,B)1/4(A,B,)1/4(A,B,)1/4(A,B)1/4 qui a été discutés en
détaille dans le paragraphe( III-8-4).
Les phases de types L10 et L12 ont été modélisées par deux ou quatre sous-
réseaux pour simplifier la modélisation d’une part et d’autre part pour faire une base de
- 72 -

données cohérente. Il y a une exception pour le composé InPt3 (L12) qui a été modélisé par
trois sous-résau grâce à la présence de lacune (défaut cristallin).
On conclut que le choix d’un modèle pour modéliser une phase intermédiaire
est commandé par deux considérations impératives :
La géométrie cristallographique (solubilité de A dans B) des phases joue un rôle
principal pour choisir un modèle adéquat.
La présence de défauts cristalline influent aussi sur notre choix du modèle.
D’après notre recherche dans la littérature on a trouvé que dans la majorité de
diagrammes de phases qui présentent plusieurs phases différentes, on doit choisir le même
modèle de sous-réseau pour pouvoir construire une base de donnée cohérente.
- 73 -

Chapitre IV
ETUDE
DES SYSTEMES
BINAIRES
(Co-Sb) et (Cd-Sb)
- 73 -

IV-1- Introduction
Pour optimiser les diagrammes de phases on utilise souvent des logiciels pour nous
faciliter la tache. Parmi les logiciels existants pour cet objectif le logiciel, BATNABIN,
PANDAT, TERMOCALC…etc.
Dans ce travail, on a utilisé le programme de BATNABIN pour calculer les
diagrammes de phases binaires (Co-Sb) et (Cd-Sb), en considérant que toutes les phases
intermédiaires sont stoechiométriques. Ce programme permet l’optimisation des grandeurs
thermodynamiques des différentes phases coexistantes dans un diagramme d’équilibre. La
méthode d’optimisation sera étudiée dans le paragraphe IV -3.
Ensuite pour étudier les phases intermédiaires existantes dans les deux systèmes on a
utilisé le programme PANDAT. On a construit pour chaque système un fichier de données
(TDB) à partir des grandeurs d’excès de la phase liquide et les grandeurs de formation des
composés définis déjà calculer par BATNABIN en plus des enthalpies libres d’excès des
phases non stœchiométriques calculées par d’autre auteurs. Ce fichier nous a permet le calcul
du diagramme de phases.
IV-2-Présentation du programme BATNABIN
Le programme BATNABIN permet le calcul des fonctions thermodynamiques des

phases existantes dans un système binaire et la restitution du diagramme de phases en utilisant
différents modèles thermodynamiques.
Après l’introduction des informations expérimentales relatives au système, celles-ci
seront traduites en équations linéaires par rapport aux coefficients ajustables. Pour obtenir des
résultats plus précis on multiplie chaque équation par un facteur poids qui dépend de notre
certitude de la mesure expérimentale lui correspondant. Ce système surdéterminé est résolu
selon la méthode des moindres carrés.
Les coefficients obtenus permettent le calcul des différentes fonctions
thermodynamiques (enthalpie, entropie et enthalpie libre) des phases existantes dans le
système, ainsi que la restitution du diagramme de phases.
Si on constate une incohérence entre les données expérimentales et les résultats de
calcul, le calcul est refait en modifiant les ordres de développement ou/et le poids des données
expérimentales, ce qui permet d’aboutir à une meilleure harmonisation entre les grandeurs
mesurées et celles calculées [6].
- 74 -

IV -3- Méthode de calcul
La détermination des coefficients ai et bi et les enthalpies libres de formation des

composés définis repose sur l'harmonisation d'un ensemble d'informations expérimentales
distinctes concernant les équilibres entre phases d'un système binaire. La condition qui régisse
l'équilibre entre deux phases est basée sur l'égalité des potentiels chimiques de chaque
élément dans les deux phases (condition d'équilibre de Gibbs) [6]. Cette condition se traduit
en deux équations dans le cas de deux phases non Stoechiométriques (équation IV-1) ou en
une seule équation dans le cas d'une phase non stoechiométrique et une phase
stoechiométrique (équation IV-2) [60].
(
 Gα (xα , T ) − Gβ x β , T) ∂Gα (xα , T )
 =
 xα − x β ∂xα

 ∂Gα (xα , T ) = ∂G β (x β , T )
(IV-1)
 ∂xα ∂x β

Gα (xα , T ) − Gc (T ) ∂Gα (xα , T )

= (IV-2)
xα − x c ∂xα
IV-4-Présentation du programme PANDAT

En 1996, le professeur Y. Austin Chang a établit le CompuTherm LLC en se basant
sur le développement d'une nouvelle génération du logiciel de calcul de diagramme de phase
Jusqu aujourd’hui, PANDAT est identifié comme progiciel intelligent (calcul automatique),
fiable, robuste et facile à utiliser et largement répandu dans le milieu universitaire et
l'industrie dans le domaine des sciences de matériaux et de technologie [61, 62].
Ces dernières années, le logiciel PANDAT est devenu un environnement
informatique intégré pour le calcul de diagramme de phase et la simulation des propriétés du
matériau des systèmes à plusieurs éléments basés sur l'approche CALPHAD (CALculation of
PHAse Diagrams) avec une bonne fiabilité et intelligence [63].
PANDAT, est un logiciel de calcul des diagrammes de phases à plusieurs
composants. En donnant un ensemble de paramètres thermodynamiques de toutes les phases
coexistant dans un système ainsi que l’ensemble de contraintes d'utilisateur, PANDAT calcul
automatiquement le diagramme de phases stables sans exiger la connaissance antérieure du
diagramme de phases. Il combine des moteurs puissants de calcul comme :
Ø PanEngine (calcul thermodynamique)
Ø PanOptimizer (optimisation de propriété)
- 75 -

Ø PanPrecipitation (simulation de précipitation)

Ø (PanGUI) Fenêtre windows avec une interface graphique facile à utiliser.
L'architecture du logiciel de PANDAT est schématiquement montrée sur la
figure-IV- 1[64,65].
Figure -IV-1. Les modules du logiciel PANDAT [66]
IV-5-Etude du système Cobalt-Antimoine:
IV-5-1- Les données expérimentales :

Le système Co-Sb est l’un des systèmes binaires de Co les plus importants grâce à ses
propriétés électroniques remarquables.
Le diagramme de phases du système (Co-Sb) comporte trois composés
intermétalliques β CoSb, γ CoSb 2 , δ CoSb3 (Fig IV-2). Les deux réactions eutectique
Liq ↔ (αCo ) + (βCoSb ) et liq ↔ (δCoSb3 ) + ( Sb) sont respectivement à T ≈ 1113°C et
T ≈ 623°C et les deux décompositions péritectique Liq + (βCoSb ) ↔ (γCoSb2 ) et

liq + γCoSb ↔ δCoSb3 sont respectivement à 931°C et 876°C. La phase intermédiaire
β présente une fusion congruente à T=1202°C.
- 76 -

Le liquidus à été déterminé par plusieurs auteurs [67-72]. Lewkonja [67], Lossew
[68] et Geller [69] ont mesurés les températures du liquidus qui sont inférieures à celles de
Feschotte [70] et Hanninger et al. [71] du coté riche en Sb. Ishida et al. [73] ont déterminés la
solubilité de Sb dans ( αCo ) qui est inférieur à celle des anciens données [72,74]. Il n’ya pas
de donnée disponible pour la solubilité de Co dans Sb.
Le diagramme de phases expérimental utilisé pour cette optimisation (Figure IV-2)
est tiré de l’ouvrage de Massalsski [36] ainsi que les données cristallographiques figurant
dans le tableau (IV-1). Les points extraits de ce diagramme utilisés dans l'optimisation sont
présentés dans le tableau (IV-2)
Oelsen [75], Predel and Vogelbein [76] ont mesurés l’enthalpie de formation de
mélange dans la gamme de 0.12 - 0.86 at% Sb à T=873K et de 0.02 - 0.66 at% à T=1060K.
Pour les données thermodynamiques des composés définis, il existe des valeurs
d’enthalpies de formation (tableau IV-5) et celles d’enthalpies libre de formation (tableau IV-
6) déterminées par différents chercheurs [75-79].
Figure IV-2 : Diagramme de phases expérimental du système (Co-Sb) [36]
- 77 -

Tableau IV-1 : Les données cristallographiques du système Co-Sb [80-82].

Phase Composition Pearson Space Strukturbericht Prototype
Sb (at.%) symbol group designation
( ) 0–2.5 cF4 Fm 3 m A1 Cu
( ) 0–2.5 hP2 P63/mmc A3 Mg
( ) 48–51 hP4 P63/mmc B81 NiAs
( ) 64–67 oP6 Pnnm C18 FeS2
( δCoSb3 ) ≈ 75 cI32 Im 3 Do 2 CoAs 3
(Sb) 100 hR2 R3m A7 αAs
IV-5-2-Résultats et discussion
Dans une première étape, on a considéré que la phase βCoSb est stœchiométrique. La
résolution d’un système d’équation surdéterminé construit à partir des données
diagrammatiques et thermodynamique injectées dans le programme BATNABIN pour un
développement de l’enthalpie d’excès de la phase liquide jusqu'à l’ordre 2 nous a permis le
calcul des coefficients ai et bi de la phase liquide et les grandeurs de formation des composés
définis. Les grandeurs thermodynamiques optimisées sont présentées dans le tableau IV-3. Le
diagramme de phases calculé à partir de ces résultats est présenté sur la figure IV-3.
Ensuite dans la deuxième étape on a utilisé les grandeurs calculées de la phase liquide
et les deux composés CoSb2 et CoSb3 ainsi que l’enthalpie libre d’excès de la phase
intermédiaire CoSb calcluée par Zhang et al. [83] pour construire un fichier de données
(TDB). Ce fichier est utilisé par PANDAT pour tracé le diagramme de phases présenté sur la
figure IV-4.
La phase CoSb présente une gamme d’homogénéité considérable au niveau des deux
cotés de la composition équiatomique. La déviation à la stoechiométrique du coté riche de Sb
est réalisée par la présence de sites vacants à la place des atomes de Co, alors que pour le coté
riche en Co elle est effectuée par la présence des atomes de Co en postions interstitielles (Fig.
IV-5).
Les trois sous-réseaux Sb1/3(Co,Va) 1/3 (Va,Co) 1/3 ont été utilisé pour décrire cette
phase [83] :
- 78 -

Tableau IV-2: Données expérimentales relatives au diagramme (Co-Sb)-BATNABIN

T (K) Phase X1(Sb) Phase X2(Sb) W
1573 Liquide 0.012490 Co 0.000000 1
1702 Liquide 0.042990 Co 0.000000 1
1661 Liquide 0.071770 Co 0.000000 1
1613 Liquide 0.102300 Co 0.000000 1
1558 Liquide 0.141200 Co 0.000000 1
1513 Liquide 0.173400 Co 0.000000 1
1465 Liquide 0.212300 Co 0.000000 1
1427 Liquide 0.249500 Co 0.000000 1
1386 Liquide 0.280000 Co 0.000000 1
1386 Liquide 0.280000 ( ) 0.500000 1
1431 Liquide 0.323700 ( ) 0.500000 1
1458 Liquide 0.367400 ( ) 0.500000 1
1482 Liquide 0.412900 ( ) 0.500000 1
1486 Liquide 0.453400 ( ) 0.500000 1
1493 Liquide 0.480300 ( ) 0.500000 1
1486 Liquide 0.522500 ( ) 0.500000 1
1472 Liquide 0.558000 ( ) 0.500000 1
1455 Liquide 0.596800 ( ) 0.500000 1
1424 Liquide 0.639000 ( ) 0.500000 1
1386 Liquide 0.686400 ( ) 0.500000 1
1355 Liquide 0.720200 ( ) 0.500000 1
1325 Liquide 0.750600 ( ) 0.500000 1
1294 Liquide 0.786100 ( ) 0.500000 1
1294 Liquide 0.818300 ( ) 0.500000 1
1256 Liquide 0.818300 ( ) 0.500000 1
1218 Liquide 0.843700 ( ) 0.500000 1
1218 Liquide 0.843700 ( ) 0.660000 1
1170 Liquide 0.884300 ( ) 0.660000 1
1191 Liquide 0.864000 ( ) 0.660000 1
1143 Liquide 0.897800 ( ) 0.660000 1
1143 Liquide 0.897800 ( δCoSb3 ) 0.750000 1
1098 Liquide 0.925000 ( δCoSb3 ) 0.750000 1
1057 Liquide 0.943600 ( δCoSb3 ) 0.750000 1
1002 Liquide 0.962300 ( δCoSb3 ) 0.750000 1
961.1 Liquide 0.975900 ( δCoSb3 ) 0.750000 1
926.9 Liquide 0.984500 ( δCoSb3 ) 0.750000 1
902.9 Liquide 0.991300 ( δCoSb3 ) 0.750000 1
- 79 -

Figure-III-3 : Diagramme de phases calculé du système (Co-Sb) avec des composés définis
Figure-III-4 : Diagramme de phases calculé du système (Co-Sb) présentant la phase

intermédiaire.
- 80 -

Figure IV-5: Illustration de sites de sous-réseau Co de la phase Co1+xSb [83].
Tableau IV- 3 : paramètres thermodynamiques optimisés

phase Paramètres
L0 =2462-32.090T.
Liquide L1=91735-59.401T.
L2= -1407+18.173T.
CoSb ∆G for = −31373+ 4.0191T

°
G Sb:Co:Va = -15370-3.95T
°
GSb:Co:Co = -12044.4+3.2T
°
G Sb:Va:Va = 2010.6+14T
°
G Sb:Va:Co = -7103.3+50T
CoSb
[83] L°Sb:Co ,Va:Co = 1064.8-15.0661T
L°Sb:Co:Co ,Va = 1064.8-15.0661T
L°Sb:Co:Va ,Co = -7522.5+4T
L°Sb:Va:Va,Co = -7522.5+4T
CoSb2 ∆G for = −25545 + 3.727T
CoSb3 ∆G for = −21841+ 3.777T
- 81 -

Sur le tableau (IV-4) on a données les coordonnés des équilibres invariants calculés
dans les deux cas et comparés aux résultats expérimentaux des différents auteurs.
Le composés définis CoSb3 se décompose par une transformation péritectique est en
accord avec le diagramme de phases expérimental. Alors que la transformation péritectique
du composé CoSb2 présente un décalage important plus de 14% de Sb pour la phase liquide
(Fig-IV-4)
L’écart entre les températures calculés et expérimentales de transformations est
acceptable si on tien en compte l’incertitude sur les résultats expérimentaux à haute
température.
On peut dire que les diagrammes de phases calculés (Fig-IV-3) et (Fig-IV-4) sont
globalement cohérent avec le diagramme expérimental avec une légère distorsion de la phase
CoSb entre xSb=46% et xSb=50%
Dans les tableaux (IV-5) et (IV-6) on a comparé les valeurs des enthalpies et des
enthalpies libres de formation des composés définis calculées et obtenues expérimentalement
par d’autres auteurs [75-79].
Tableau IV-4 : Les coordonnées des équilibres invariants calculés comparés aux résultats
expérimentaux.
xsb% T(K) Type de référence
Réaction
transformation
23.6 65 49.6 1391 eutectique [67]
23.3 - - 1366 eutectique [68]
25.2 - - 1375 eutectique [69]
29 22 44 1385 ± 20 eutectique [70]
28 - 43.5 1391 eutectique [71]
- 48 - 1370 eutectique [72]
L ↔ αCo + βCoSb - - 42 1386 eutectique [84]
25.5 25 48.5 1386 eutectique [85]
25.37 38 46.2 1391 eutectique [83]
27 100 50 1381 eutectique Calcul 1

31 100 46 1324 eutectique Calcul 2
- 1466 congruente [67]

L ↔ βCoSb 50 1488 congruente [70]
47 ± 0.5 1493 ± 5 congruente [71]
47 1475 congruente [84]
49.08 49.08 1493 congruente [83]
- 82 -

50 1486 congruente Calcul 1

46-50 1382 congruente Calcul 2
86 52.5 66.7 1209 péritectique [86]

83.1 - 66.7 1170 ± 6 péritectique [67]
81.3 - 66.7 1179 ± 3 péritectique [68]
83 - 66.7 1167 péritectique [69]
89 51.6 66.7 1202 ± 2 péritectique [70]
- 52.5 66.7 1209 péritectique [71]
L + βCoSb ↔ γCoSb2
83 - 66.7 1192 péritectique [84]
- - 66.7 1204 péritectique [87]
89.26 52.53 66.7 1209 péritectique [83]
82 50 67 1235 péritectique Calcul 1

90 66.7 75 1146 ± 5 péritectique [70]
- 66.7 75 1147 péritectique [71]
90 66.7 75 1132 péritectique [84]
90 66.7 75 1149 péritectique [85]
L + γCoSb2 ↔ δCoSb3 - 66.7 75 1149 péritectique [87]
95.31 66.7 75 1147 péritectique [83]

97.1 75 - 889 ± 3 eutectique [67]
100 75 - 898 ± 5 eutectique [68]
97 75 - 891 eutectique [69]
- 75 - 891 ± 2 eutectique [70]
- 75 - 902 eutectique [71]
- 75 - 896 eutectique [84]
L ↔ δCoSb3 + Sb 100 75 - 896 eutectique [85]
99.98 75 - 903.7 eutectique [83]
98 75 - 899 eutectique Calcul 1

98.7 75 - 899 eutectique Calcul 2
Suite tableauIV-4
- 83 -

Tableau IV-5: Enthalpie de formation des composés définis

T(K) référence
Composé ∆H f (kj )
mol
-26.16 [76]
-20.55 873 [75]
CoSb -16.11 850 [77, 78]
-31.373 Notre calcul
-16.11 873 [75]
CoSb2 -15.55 850 [77, 78]
-15 873 [75]
-13.88 850 [77, 78]

CoSb3
Tableau IV-6 : Enthalpie libre de formation des composés définis

Composé T (K) référence
∆G f (kj )
mol
-15.80 [79]
873
CoSb
-27.86 Notre calcul
-16.80 [79]
873
CoSb2
-14.4 [79]
CoSb3 873
-18.54 Notre calcul
- 84 -

IV-6-Etude du système Cadmium-Antimoine

IV-6-2- Les données Expérimentales
Le diagramme de phases du système (Cd-Sb) tiré de l’ouvrage de Massalski présente

une seule phase intermédiaire avec deux réactions eutectiques (Fig- III-6):
Cd + CdSb ↔ L à T=290°C et Sb + CdSb ↔ L à T=445°C et une fusion congruente
à T=456°C. Les données cristallographiques sont présentées dans le tableau (III-8) [36].
La structure cristallographique de la phase intermédiaire CdSb est représentée sur la
figure-IV-7 [14].
Figure -IV-6 : Diagramme de phases expérimental du système (Cd-Sb) [36].
Figure IV-7: La structure cristallographique de la phase CdSb [14].
- 85 -

Tableau III-8 : Les données cristallographiques du système Cd-Sb [36]

Phase Composition Sb Pearson Space Strukturbericht Prototype
(at.%) symbol group designation
Cd 0 hP2 P6 3 /mmc A3 Mg
CdSb 50–51 Op16 Pbca Be CdSb
Sb 100 hR2 R3m A7 αAs
Les points extraits de ce diagramme utilisés dans l'optimisation figurent dans le

tableau (IV-9).
On a également utilisé comme données l’enthalpie et l’entropie de formation de la
phase CdSb déterminé par [88] :
ΔHfor(CdSb) = -7500 J/mole
ΔSfor(CdSb) = -2.3 J/mole. K
Ainsi que des enthalpies et des entropies de mélange de la phase liquide mesurées par
Geffken et al. à 773K [89].
Tableau IV-9 : Données expérimentales relatives au diagramme (Cd-Sb).

T (K) Phase X1(Sb) Phase X2(Sb) W
594.108000 Liquide 0.000000 Cd 1.000000 1
583.000000 Liquide 0.033000 Cd 1.000000 1
573.000000 Liquide 0.058000 Cd 1.000000 1
563.000000 Liquide 0.075000 Cd 1.000000 1
718.000000 Liquide 0.570000 Sb 1.000000 1
773.000000 Liquide 0.683300 Sb 1.000000 1
823.000000 Liquide 0.800000 Sb 1.000000 1
873.000000 Liquide 0.916000 Sb 1.000000 1
903.755000 Liquide 0.990000 Sb 1.000000 1
563.000000 Liquide 0.075000 CdSb 0.500000 1
573.000000 Liquide 0.090000 CdSb 0.500000 1
583.000000 Liquide 0.108000 CdSb 0.500000 1
594.108000 Liquide 0.141600 CdSb 0.500000 1
623.000000 Liquide 0.200000 CdSb 0.500000 1
673.000000 Liquide 0.308000 CdSb 0.500000 1
723.000000 Liquide 0.533 CdSb 0.500000 1
725.000000 Liquide 0.5 CdSb 0.500000 1
713.000000 Liquide 0.558000 CdSb 0.500000 1
718.000000 Liquide 0.570000 CdSb 0.500000 1
- 86 -

IV-3-2- Résultats et discussion

Les fonctions thermodynamiques de la phase liquide et le composé défini optimisé
après le développement de l’enthalpie libre d’excès de la phase liquide jusqu'à l’ordre 2 sont
représentées dans le tableau (IV-10). Le diagramme de phase restitué est présenté sur
la figure IV-8.
Ensuite on a construit un fichier TDB à partir de l’enthalpie libre d’excès de la phase
liquide et l’enthalpie libre de la phase intermédiaire CdSb calculé par Liu et Tedenac. [90]
(Tableau IV-10).
Le diagramme de phases calculé à partir de ce fichier et en utilisant le programme
PANDAT est représenté sur la figure IV-9.
Le modèle des sous-réseaux proposé par Liu et Tedenac [90] est de la forme
Sb(Cd,Va). Ce choix du modèle est complètement compatible avec les résultats qu’on a
présentées dans le troisième chapitre pour un composé équiatomique tableau (III-1).
Les différentes températures des transformations invariantes expérimentales et
calculées de ce système sont donnés dans le tableau (IV-11).
On conclut à partir de ces résultats que le processus d’optimisation (dualité
BATNABIN et PANDAT) utilisé pour calculer le diagramme de phases (Cd-Sb) est très
efficace. On a abouti à une bonne cohérence entre le diagramme expérimental et celui calculé
ainsi que les températures des diverses transformations invariantes qui sont globalement
identiques avec celle de diagramme de phases expérimental.
Les grandeurs thermodynamiques de formations de la phase intermédiaire et de
mélange de la phase liquide, calculées sont en bon concordance avec les résultats
expérimentaux tableau (IV-12, 13).
Tableau IV-10 : paramètres thermodynamiques optimisés.
phase Paramètres
L0 =-9311-3.247T.
Liquide
L1=-1065-3.621T.
L2= 1247611.772T.
CdSb ∆G for = −9964 + 5.563T

°
G Sb:Cd = -24952.46+121.02062T
(CdSb) °
GSb:Co:Co = +10000
[90]
L°Sb:Cd ,Va = +1.70161893+14.1708673
- 87 -

Figure-IV-8 : Diagramme de phases calculé du système (Cd-Sb) avec un composé

définis
Figure-IV-9: Diagramme de phases calculé du système (Cd-Sb) avec une phase

intermédiaire.
- 88 -

Tableau IV-11 : les coordonnés des équilibres invariants

Type de
Réaction Xsb% T (K) Référence
transformation
7 100 50 563 eutectique [36]
7.4 100 50 561.5 eutectique Calcul 1
L ↔ Cd + CdSb
L ↔ CdSb 50 729.1 congruente Calcul 1
50 728.7 congruente Calcul 2
57 51 100 718 eutectique [36]
L ↔ Sb + CdSb 5.89 100 50 721.8 eutectique Calcul 1
Tableau IV-12 : Grandeurs intégrales de formation de composé CdSb.

Le composé Référence
∆H f (kj ) ∆S f (kj .)
mol mol.K
-7.5 -2.3 [88]
CdSb -9.964 -5.55 Notre calcul
Tableau IV-13 : Grandeurs intégrales de mélange de la phase liquide à 773K

x H (kj ) S (kj ) Référence
mol mol.K
-1.81 4 [89]
0.1
-1.96 3.92 Notre calcul
-0.6 5.85 [89]
0.2
-1.45 6.4 [89]
0.3
-2.2 6.5 [89]
0.4
-2.2 6.5 [89]
0.5
-1.6 6.5 [89]
0.6
-1.11 5.9 [89]
0.7
-0.7 4.9 [89]
0.8
-0.3 2.9 [89]
0.9
- 89 -

CONCLUSION
GENERALE
- 89 -

CONCLUSION GNENERALE
La modélisation des phases non-stoechiométriques, qui ont été considérées depuis

longtemps comme stoechiométriques pour simplifier leurs traitement est récemment
effectuée par la méthode des sous-réseaux de Hillert et Staffanson. Cette méthode permet de
décrire une phase à l’aide de plusieurs sous-réseaux (CEF). Ces derniers pouvant admettre
différentes espèces chimiques (ions, molécules, atomes, lacunes).
On a introduit tout d’abord ce modèle dans le cas général en prenant soin de définir
tous les termes impliqués puis on a étudié quelques cas particuliers que le modèle permet de
décrire dans les systèmes métalliques. L’étude de cas simples, lorsque la substitution d’un
seul ou de deux sous-réseaux est envisagée, nous permettra d’introduire des notations plus
explicites que dans le cas général et de représenter les grandeurs thermodynamiques
introduites lors de la présentation du modèle. Le calcule de l’enthalpie libre de Gibbs pour le
modèle de sous-réseau est illustrée dans ce travail pour différents cas et elle peut être divisée
en trois contributions, c’est le cas similaire à celle employé pour les solutions de substitution
mais avec un calcul qui est effectué pour une mole de sites de la phaseϕ au lieu d’une mole
d’atome. En raison il y a une relation entre les deux enthalpies libres qui a été définit dans
notre travail. Mais si la description de la phase ϕ ne fait pas intervenir de lacunes, les deux
enthalpies libres sont égales et la composition des éléments de la phase égale à celle de
fraction de site.
Dans ce travail et après une étude détaillée du modèle des sous-réseaux on a réalisé
un recensement complet sur les différentes phases étendues existantes dans la littérature et
leur modèle de sous-réseau adéquat pour les décrire. On conclut que le choix d’un modèle
pour modéliser une telle phase est commandé par deux considérations impératives :
La géométrie cristallographique (solubilité de A dans B) de phases joue un rôle
principal pour choisir un modèle adéquat.
La présence de défauts cristalline influent aussi sur notre choix du modèle.
D’après notre recherche dans la littérature on a trouvé que la majorité de diagrammes
de phases présentent plusieurs phases différentes, dans ce cas on doit choisir le même modèle
de sous-réseau pour pouvoir construire une base de donnée cohérente.
- 90 -

Enfin est afin d’appliquer ce modèle à des système binaires, on a optimisé les deux
systèmes Co-Sb et Cd-Sb . On premier on a considéré que toutes les phases intermédiaires
sont stoechiométrique pour simplifier les calculs par le programme BATNABIN. Ensuite on
a introduit les enthalpies libres d’excès des phases intermédiaires et les autres enthalpies
libres déjà calculer dans le programme PANDAT pour tracer les diagrammes de phases. Les
résultats obtenus sont en bon accord avec les données expérimentales.
- 91 -

REFERENCE
BIBLIOGRAPHIQUE
91

Bibliographie
[1] M. Hillert «Phase equilibria, phase diagrams and phase transformations, », Cambridge
University Press, (1998).
[2] N. Saunders and A. P. Miodownik, CALPHAD «Calculation of Phase Diagrams»,
University of Surrey, Guildford, UK, (1998).
[3] I. Ansara, C. Chatillon, H. L. Lukas, T. Nishizawa, H. Ohtani, K. Ishida, M. Hillert, B.
Sundman, B. B. Argent, A. Watson, T. G. Chart and T. Anderson, Calphad, 14, 203,
(1994).
[4] J.H. Hildebrand, J. Ame. Chem. Soc, 51, 66-80, (1929).
[5] B. Sundman, B. Jansson and J. O. Anderssson, Calphad, 9, 153, (1985).
[6] Y. Djaballah, thèse de Doctorat, Université de Batna, (2005).
[7] J. D. Eshelby , Solid state physics, 3, 7-79, (1956).
[8] J. Hertz, «Thermodynamique de l’Equilibre Chimique dans la Matière », Edition du
laboratoire de thermodynamique métallurgique, Nancyl (France), (1995).
[9]. J. Benard, A. Michel, J. Philbert et J. Talbot « Métallurgie générale », Massonet Cie,
(1986).
[10] C. Henri, D. Novion et D. Lesueur. « L’état métallique et propriétés atomique-structure
atomique », Technique de l’ingénieur, (1993).
[11]Y. Adda, J. M. Dupouy, J. Philbert, Y. Quere. «Alliage-défauts cristallins. tomIII»,
laboratoire d’optique électronique du C.N.R.S, Toulouse, (1977).
[12] C. H. Novion, D. Lesueur, « Létat métalliques propriétés atomiques», Technique de
l'ingénieur, M 35, (1979).
[13] D. Schwarzenbach et G. Chapuis, « Cristallographie » Presses polytechniques romandes,
(2006).
[14] www.nrL.navy.mil/Lattice/pearson, «Crystal Lattice Structures».
[15] J. Hertz, «Diagramme d’Equilibre, Alliages Binaires », Technique de l’ingénieur, M 70,
(1999).
[16] A. Prince, «Alloy phase equilibria», Elsevier, Amsterdam, (1966) .
[17] http://wwwens.uqac.ca/chimie/Thermochimie/Chap_htm/CHAP_11.html
[18] L. Guillet, «Diagramme de Phases en Métallurgie », Edition Masson & Cie, (1964).
92

[19] J. Hertz. «Application de la thermodynamique à la détermination des diagrammes de

phases », Ecole d’été du C.N.R.S - Aussois (France) ,(1990).
[20] I. Ansara. «Diagramme de phases et stoechiométriex»., Masson et Cie, (1972)
[21] I. Diop, thèse de Doctorat, Université Nancy I, (2008).
[22] J. Manenc, «Thermodynamique Structurale des Alliages », Presse Universitaire de
France, (1972).
[23] P. Haasen, « Physical Metallurgy », Cambridge University Press, (1978).
[24] G. Maurin, «cours thermochimie», (2005).
[25] Y. Djaballah, Mémoire de Magister, Université de Batna, (1999).
[26] D. Jacques, «thermodynamique», ENS de Cachan, (2007)
[27]Y. Austinchany, S. Chen, F. Zhang, X. Yan, F.Xie, R. Scmid-Fetzer and W. A. Oates,
Prog. Mater. Sci., 49, 313-345, (2004).
[28] R. Fourmentin, Thèse de doctorat, Université de Lille, (2004).
[29] O. Redlich and A. Kister, Md. Eng. Chem., 40, 345-348 (1948).
[30] N. Dupin, thèse de doctorat, Institut national polytechnique de Grenoble, (1992)
[31] G. Xiaoting, Master memory, Case Western reserve university, (2008).
[32] D. Keyzer, these de doctorat, Katholieke Universiteit Leuven, (2008).
[33] B. Sundman, S. G. Fries, and W. A. Oates. Calphad, 22(3):335- 354, (1998).
[34] M. Xiaoyan, L. Changrong, D. Zhenmin and Z.Weijing, Journal of Alloys and
Compounds 370, 149–158, (2004).
[35] G. Yongliang, G. Cuiping, L. Changrong, C. Shaohua and D. Zhenmin , Journal of
Alloys and Compounds 479, 148–151, (2009).
[36] T. B. Massalski, «Binary Alloys phase diagrams» ASM International, (1990).
[37] Y. Liang, G. Cuiping, L. Changrong and D. Zhenmin, Journal of Alloys and Compounds
475, 220–226, (2009).
[38] Y. Liang, D. Zhenmin, G. Cuiping and L. Changrong, Journal of Alloys and Compounds
455, 236-242, (2008).
[39] M. Lia, C. Lia, F. Wang, Z. Weijing, Journal of Alloys and Compounds, 14, 39-46,
(2006).
[40] V. T. Witusiewicz, A. A. Bondar, U. Hecht, S. Rexa, T. Ya. Velikanova, Journal of
Alloys and Compounds, 465, 64–77, (2008).
[41] H. Okamoto, Journal of Phase Equilibria and Diffusion 30, 3, (2009).
[42] Y. Zhang, L. Changrong, D. Zhenmin and T. Geng, Calphad, 11,001, (2007).
[43] M. H. G. Jacobs, R. Schmid-Fetzer, Calphad, 33, 170-178, (2009).
93

[44] D. Arias, j. Abriata, L. Gribaudo, Journal of Nuclear Materials 229, 24-28, (1996).
[45] T. Benlahrache, thèse de doctorat, Université Nancy 1, (2008).
[46] Y. Liang, G. Cuiping, L, Changrong and D. Zhenmin, Journal of Alloys and Compounds
460, 314–319, (2008).
[47] A. Durga, K. C. Han Kumar, Calphad 34, 200-205, (2010).
[48] W. Han, L. Mei, Journal of Alloys and Cornpounds, 507, 521-525, (2010).
[49] C. Zhang, J. Zhu , A. Bengtson, D. Morgan, F. Zhang, Y. Yang and Y.A. Chang, Acta
Materialia 56, 5796–5803, (2008).
[50] L. Mei , G. Cuiping, L. Changrong , D. Zhenmin, Journal of Alloys and Compounds 481,
283–290, (2009).
[51] M. Pingli, K. Hana andY. Xin, Journal of Alloys and Compounds 464, 190-196, (2008).
[52] N. David, T. Benlaharche, J.M. Fiorani, M. Vilasi, Intermetallics 15, 1632-1637, (2007).
[53] R. Krachler, H. Ipser, et K. L. Komarek, J. Phys. Chem. Solids, 50, (11), 1127-1135
(1989).
[54] J. M. Joubert, Progress in Materials Science 53, 528- 583, (2008).
[55] G. Renaudin, L. Guenee and K. Yvon, J. Alloys Compd. 350, 145, (2003).
[56] L. Schlapbach and A. Zuttel, Nature 414, 353, (2001).
[57] K. Yano, I. Umehara, T. Miyazawa, Y. Adachi and K. Sato, Physica B 367, 81, (2005).
[58] C. Iuliana, thèse de doctorat, Université de Bordeaux et de Bucarest, (2008).
[59] I. Ansara, B. Sundman and P. Willemin, Acta metall., 36, (4), 977-982, (1988).
[60] N. Djenane, mémoire de magistère, Université de Batna, (2009).
[61] S.L. Chen, S. Daniel, F. Zhang, Y.A. Chang, X.Y. Yan, F.Y. Xie, R. Schmid-Fetzer, and
W.A. Oates, Calphad 26(2), 175-188, (2002).
[62] S.L. Chen, F. Zhang, S. Daniel, F.Y. Xie, X.Y. Yan, Y.A. Chang, R. Schmid-Fetzer, and
W.A. Oates, JOM 55(12), 48-5, (2003).
[63] W. Cao, S. L. Chen, F. Zhang, K. Wu, Y. Yang, Y. A. Chang, R. Schmid-Fetzer and W.
A. Oates, Calphad, 33(2), 328-342, (2009).
[64]S. Daniel, F. Zhang, Y.A. Chang, X. Y. Yan, F. Y. Xie, R. Schmid-Fetzer, W.A. Oatesd ,
Calphad 26, 2 175-188, (2002).
[65] S. L. Chen, F. Zhang, S. Daniel, F. Y. Xie, X. Y. Yan, Y. A. Chang, R. Schmid-Fetzer,
and W.A. Oates, Springer, (2003).
[66] http://computherm.com/main.html#contact.
[67] K. Lewkonja, Z. Anorg. Chem. 59, 305- 312, (1908).
[68] K. Lossew, Zhur. Russ, Fiz.-Khim. Obshchestva. 43, 375 -392, (1911).
94

[69] W. Geller, Arch. Eisenhflttenwes. 13, 263 -266, (1939).

[70] P. Feschotte, D. Lorin, J. Less-Common Met. 155, 255- 269, (1989).
[71] G. Hannninger, H. Ipser, P. Terzieff, K.L. Komarek, J. Less-Common Met. 166, 103 -
114, (1990).
[72] U. Haschimoto, J. Japan. Inst. Met. 1, 177 -190, (1937).
[73] K. Ishida, M. Hasebe, N. Ohnishi, T. Nishizawa, J. Less-Conimon Met. 114, 361- 373,
(1985).
[74] W. Kôster, E. wager, Z. Metallkd. 29, 230 -232, (1937).
[75] F. Kôrber, W. Oelsen, Mitt. Kaiser-Wilhelm-Inst. Eisenforsch. Dûsseldorf 19, 209- 219,
(1937).
[76] B. Predel, W. Vogelbein, Thermochim. Acta. 24, 155 -165, (1978).
[77] V. A. Geiderikh, Y. I. Gerasimov, Doklady Akademii Nauk SSSR 140, 391- 393, (1961).
[78] Y. I. Gerasimov, A. V. Nikol’skaya, «Solid State Phys», 113 -121, (1965).
[79] R. Hultgren, P.D. Desai, D.T. Hawkins, M. Gleiser, K.K. Kelley, «Selected Values of die
Thermodynamic Properties of Binary Alloys», ASM, (1973).
[80] T. Nishizawa, K. Ishida, Bull. Alloy Phase Diagram, 4, 387 -390, (1983).
[81] W. B. Pearson, « Handbook of Lattice Spacings and Structures of Metals and Alloys»,
vol. 2, Pergamon, New York, (1967).
[82] H. W. King, Bull. Alloy Phase Diagrams 2, 401-402, (1981).
[83] Y. Zhang, C. Li, Z. Du, T. Geng, Calphad, 32, 1. 56-63, (2008).
[84] L. D. Dudkin, N. K. Abrikosov, Zh. Neorgan. Khim. 1, 2096-2105, (1956).
[85] K. Ishida, T. Nishizawa, Bull. Alloy Phase Diagrams 11, 243- 248, (1990).
[86] H. Okamoto, J. Phase Equilibria. 12, 244- 245, (1991).
[87] A. Kjekshus, T. Rakke, ACTA Chem. Scand. A 31, 517- 529, (1977).
[88] L. Zabdyr, Arch. Hutnictwa 24, 95, (1979).
[89] R. Geffken, K.L. Komarek, E. Miller, Trans. AIME 239, 1151, (1967).
[90]Y.Liu and J.C. Tedenac, Calphad, 33, 684-694, (2009).
95

Annexe
- 95 -

A1- Fichier TDB du système Co-Sb
$ -------------------------------------------------------------------------
$
ELEMENT /- ELECTRON_GAS 0.0000E+00 0.0000E+00 0.0000E+00!
ELEMENT VA VACUUM 0.0000E+00 0.0000E+00 0.0000E+00!
ELEMENT CO HCP_A3 5.8933E+01 0.0000E+00 0.0000E+00!
ELEMENT SB RHOMBOHEDRAL_A7 1.2175E+02 5.8702E+03 4.5522E+01!
$
$ -------------------------------------------------------------------------
$
FUNCTION GHSERCO 298.15 +310.241+133.36601*T-25.0861*T*LN(T)
-.002654739*T**2-1.7348E-07*T**3+72527*T**(-1); 1768 Y
-17197.666+253.28374*T-40.5*T*LN(T)+9.3488E+30*T**(-9); 6000 N !
FUNCTION GHSERSB 298.15 -9242.858+156.154689*T-30.5130752*T*LN(T)
+.007748768*T**2-3.003415E-06*T**3+100625*T**(-1); 903.78 Y
-11738.83+169.485872*T-31.38*T*LN(T)+1.616849E+27*T**(-9); 2000 N !
FUNCTION LCO 298.15 +15395.278+124.434078*T
-25.0861*T*LN(T)-.002654739*T**2-1.7348E-07*T**3
+72527*T**(-1)-2.19801E-21*T**7; 1768 Y
-846.61+243.599944*T-40.5*T*LN(T); 6000 N !
FUNCTION LSB 298.15 +10579.47+134.231525*T
-30.5130752*T*LN(T)+.007748768*T**2-3.003415E-06*T**3
+100625*T**(-1)-1.74847E-20*T**7; 903.78 Y
+8175.359+147.455986*T-31.38*T*LN(T); 2000 N !
FUNCTION CO_FCC 298.15 +737.832+132.750762*T
-25.0861*T*LN(T)-.002654739*T**2-1.7348E-07*T**3
+72527*T**(-1); 1768 Y
-16770.075+252.668487*T-40.5*T*LN(T)+9.3488E+30*T**(-9); 6000 N !
FUNCTION SB_FCC 298.15 +10631.142+142.454689*T
-30.5130752*T*LN(T)+.007748768*T**2-3.003415E-06*T**3
+100625*T**(-1); 903.78 Y
+8135.17+155.785872*T-31.38*T*LN(T)+1.616849E+27*T**(-9); 2000 N !
FUNCTION SB_HCP 298.15 +10631.142+143.154689*T
-30.5130752*T*LN(T)+.007748768*T**2-3.003415E-06*T**3
+100625*T**(-1); 903.78 Y
+8135.17+156.485872*T-31.38*T*LN(T)+1.616849E+27*T**(-9); 2000 N !
FUNCTION UN_ASS 298.15 0; 300 N !
TYPE_DEFINITION % SEQ *!
DEFINE_SYSTEM_DEFAULT ELEMENT 2 !
DEFAULT_COMMAND DEF_SYS_ELEMENT VA /- !
$
$ --------LIQUID PHASE-----------------------------------------------------------------
$
PHASE LIQUID % 1 1.0 !
CONSTITUENT LIQUID :CO,SB : !
PARAMETER G(LIQUID,CO;0) 298.15 +LCO; 2000 N !
PARAMETER G(LIQUID,SB;0) 298.15 +LSB; 2000 N !
- 96 -

PARAMETER G(LIQUID,CO,SB;0) 298.15 2462-32.090*T; 3000 N !

PARAMETER G(LIQUID,CO,SB;1) 298.15 91735-59.401*T; 3000 N !
PARAMETER G(LIQUID,CO,SB;2) 298.15 -1407+18.173*T; 3000 N !
$ --------CO PHASE -------------------------------------------------
$
PHASE CO_A3 ,, 1 1!
Constituent CO_A3 :CO :!
Parameter G(CO_A3,CO:VA;0) 2.98150E+02 +GHSERCO; 3.00000E+03 N !
$
$ --------SB PHASE -------------------------------------------------
$
PHASE RHOMBOHEDRAL_A7 % 1 1.0 !
CONSTITUENT RHOMBOHEDRAL_A7 :SB : !
PARAMETER G(RHOMBOHEDRAL_A7,SB;0) 298.15 +GHSERSB#; 2000 N !
$
$ Compound phases ----------------------------------------------------------
$
$ COSB sub-lattice model0--------------------------------------------------
PHASE COSB % 3 .3333 0.3333 0.3333 !
CONSTITUENT COSB : SB : CO, VA : VA,CO : !
PARAMETER G(COSB,SB:CO:VA;0) 298.15 -15370-3.95*T
+0.4920984*T*LN(T)+2.260266E-03*T**2-1.0735435E-6*T**3
-20074*T**(-1)+1.5E-21*T**7
+0.333333*GHSERCO#+0.3333*GHSERSB#; 6000 N !
PARAMETER G(COSB,SB:CO:CO;0) 298.15 -12044.4+3.2*T
PARAMETER G(COSB,SB:VA:VA;0) 298.15 2010.6+14*T
+0.3333*GHSERSB; 6000 N !
PARAMETER G(COSB,SB:VA:CO;0) 298.15 -7103.3+50*T
PARAMETER G(COSB,SB:CO,VA:CO;0) 298.15 1064.8-15.0661*T; 6000 N !
PARAMETER G(COSB,SB:CO,VA:VA;0) 298.15 1064.8-15.0661*T; 6000 N !
PARAMETER G(COSB,SB:CO:CO,VA;0) 298.15 -7522.5+4*T; 6000 N !
PARAMETER G(COSB,SB:VA:CO,VA;0) 298.15 -7522.5+4*T; 6000 N !
$ COSB2 Stoechiometrique compound--------------------------------------------------

PHASE COSB2 % 2 .3333 0.6667 !
CONSTITUENT COSB2 : CO : SB : !
PARAMETER G(COSB2,CO:SB;0) 298.15 -
25545+3.727*T+0.333333*GHSERCO#+0.6667*GHSERSB#; 6000 N !
$ COSB2 Stoechiometrique compound--------------------------------------------------
PHASE COSB3 % 2 .25 0.75 !
CONSTITUENT COSB3 : CO : SB : !
PARAMETER G(COSB3,CO:SB;0) 298.15 -
21841+3.777*T+0.25*GHSERCO#+0.75*GHSERSB#; 6000 N !$
- 97 -

A2- Fichier TDB du système CdSb

$ -------------------------------------------------------------------------
$
ELEMENT /- ELECTRON_GAS 0.0000E+00 0.0000E+00 0.0000E+00!
ELEMENT VA VACUUM 0.0000E+00 0.0000E+00 0.0000E+00!
ELEMENT CD HCP_A3 1.1241E+02 6.2509E+03 5.1798E+01!
ELEMENT SB RHOMBOHEDRAL_A7 1.2175E+02 5.8702E+03 4.5522E+01!
$
$ -------------------------------------------------------------------------
$
FUNCTION GHSERCD 298.15 -7083.469+99.506198*T
-22.0442408*T*LN(T)-6.273908E-3*T**2-6966*T**(-1);
594.22 Y -20064.971+256.812233*T-45.1611543*T*LN(T)
+8.832011E-3*T**2-0.899604E-6*T**3+1241290*T**(-1);
1500.00 Y -9027.489+148.20548*T-29.7064*T*LN(T);
1600.00 N !
FUNCTION GHSERSB 298.15 -9242.858+156.154689*T-30.5130752*T*LN(T)

+.007748768*T**2-3.003415E-06*T**3+100625*T**(-1); 903.78 Y
-11738.83+169.485872*T-31.38*T*LN(T)+1.616849E+27*T**(-9); 2000 N !
FUNCTION LCD 298.15 -955.025+89.209282*T

-22.0442408*T*LN(T)-6.273908E-3*T**2-6966*T**(-1);
400.00 Y 21716.884-371.046869*T+53.1313898*T*LN(T)
-115.159917E-3*T**2+28.899781E-6*T**3-1271815*T**(-1);
594.22 Y -3252.303+138.251107*T-29.7064*T*LN(T);
1600.00 N !
FUNCTION LSB 298.15 +10579.47+134.231525*T
-30.5130752*T*LN(T)+.007748768*T**2-3.003415E-06*T**3
+100625*T**(-1)-1.74847E-20*T**7; 903.78 Y
+8175.359+147.455986*T-31.38*T*LN(T); 2000 N !
FUNCTION UN_ASS 298.15 0; 300 N !
TYPE_DEFINITION % SEQ *!
DEFINE_SYSTEM_DEFAULT ELEMENT 2 !
DEFAULT_COMMAND DEF_SYS_ELEMENT VA /- !
$
$ --------LIQUID PHASE-----------------------------------------------------------------
$
PHASE LIQUID % 1 1.0 !
CONSTITUENT LIQUID :CD,SB : !
PARAMETER G(LIQUID,CD;0) 298.15 +LCD; 2000 N !
PARAMETER G(LIQUID,SB;0) 298.15 +LSB; 2000 N !
Parameter G(LIQUID,CD,SB;0) 298.15 -9311-3.247*T; 3.20000E+03 N !
Parameter G(LIQUID,CD,SB;1) 298.15 -1065-3.621*T; 3.20000E+03 N !
Parameter G(LIQUID,CD,SB;2) 298.15 12476-11.772*T; 3.20000E+03 N !
$
$
- 98 -

$ --------CO PHASE (Stoechiométrique) -------------------------------------------------

$
PHASE CD_A3 ,, 1 1!
Constituent CD_A3 :CD :!
Parameter G(CD_A3,CD:VA;0) 2.98150E+02 +GHSERCD; 3.00000E+03 N !
$
$ --------SB PHASE (Stoechiométrique) -------------------------------------------------
$
PHASE RHOMBOHEDRAL_A7 % 1 1.0 !
CONSTITUENT RHOMBOHEDRAL_A7 :SB : !
PARAMETER G(RHOMBOHEDRAL_A7,SB;0) 298.15 +GHSERSB#; 2000 N !
$
$ CDSB sub-lattice model--------------------------------------------------
PHASE CDSB % 2 1 1 !
CONSTITUENT CDSB :SB : CD,VA : !
param G(CDSB,SB:VA;0) 2.98150E+02 +GHSERSB +10000; 6.00000E+03 N 09Liu !
param G(CDSB,SB:CD;0) 2.98150E+02 -24952.46+121.02062*T-
15.588758*T*LN(T)+GHSERCD +GHSERSB ; 1.50000E+03 N !
param G(CDSB,SB:CD,VA;0) 2.98150E+02 1.70161893E+0+1.41708673E+01*T;
6.00000E+03 N !
- 99 -

Alliage Homogene-Héterog-Solution Solide PDF

Transféré par

Droits d'auteur :

Formats disponibles

Alliage Homogene-Héterog-Solution Solide PDF

Transféré par

Informations du document

Titre original

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

Alliage Homogene-Héterog-Solution Solide PDF

Transféré par

Droits d'auteur :

Formats disponibles

REPUBLIQUE ALGERIEN DEMOCRATIQUE ET POPULAIRE

Ministère de l’enseignement supérieur et de la recherche scientifique

Présenté au département des sciences de la matière

Pour obtenir le diplôme de

Etude d’une variété de solutions solides binaires

Président : A. BALGACEM-BOUZIDA Pr U. Batna

PDF created with pdfFactory Pro trial version www.pdffactory.com

Nikola Tesla (1856-1943)

PDF created with pdfFactory Pro trial version www.pdffactory.com

Tous d’abord je remercie le bon DIEU le tout puissant.

Je tiens à remercie très vivement Monsieur Djaballah Yassine maitre de

Je remercie également Monsieur Belgacem Bouzida Aissa, professeur à

Je n’oublie surtout pas monsieur Belbacha El-Djemai, professeur à l’université

Ma profonde reconnaissance va ensuite envers, Madame Bouharkat Fouzia,

Mes remerciements vont également à monsieur Salah Benalia, maître de

Je tiens cependant à remercier tout particulièrement Amel Naceur professeur

PDF created with pdfFactory Pro trial version www.pdffactory.com

ALLIAGE ET DIAGRAMME DE PHASES

PDF created with pdfFactory Pro trial version www.pdffactory.com

LES BASES THERMODYNAMIQUES DES DIAGRAMMES DE PHASES

PDF created with pdfFactory Pro trial version www.pdffactory.com

LE MODELE DES SOUS-RE SEAUX

PDF created with pdfFactory Pro trial version www.pdffactory.com

ETUDE DES SYSTEME BINAIRE PAR LA METHODE DES SOUS-RESEAUX

IV-1- introduction ………………………………………………………………………………74

PDF created with pdfFactory Pro trial version www.pdffactory.com

PDF created with pdfFactory Pro trial version www.pdffactory.com

PDF created with pdfFactory Pro trial version www.pdffactory.com

PDF created with pdfFactory Pro trial version www.pdffactory.com

I-1-1- Alliage homogène

I-1-2- Alliage hétérogène

PDF created with pdfFactory Pro trial version www.pdffactory.com

Figure- I-2 : Alliage hétérogène

I-2-1- Homogénéité de la composition :

PDF created with pdfFactory Pro trial version www.pdffactory.com

D’autre part, la thermodynamique phénoménologique est toujours globalisante par

Figure - I-3 : Principe de l’analyse par microsonde électronique de Castaing.

I-2-2- Structure d’une phase

I-2-3- Aspect densité ou masse volumique d’une phase

PDF created with pdfFactory Pro trial version www.pdffactory.com

I-3-Les solutions solides

I-3-1-Les solutions solides primaires

PDF created with pdfFactory Pro trial version www.pdffactory.com

I-3-2-Les solutions solides d'insertion (ou interstitielles)

I-3-3-Les Solutions Solides de Substitution

PDF created with pdfFactory Pro trial version www.pdffactory.com

Solution désordonnées Solution parfaitement ordonnées Rassemblement

Figure- I-4 : Type de solutions solides de substitution.

I-3-4-Les Solutions solides intermédiaires

PDF created with pdfFactory Pro trial version www.pdffactory.com

I-3-4-1- Solutions ordonnées et désordonnées

• La structure L 10 (type AuCu)

PDF created with pdfFactory Pro trial version www.pdffactory.com

cube, ce qui correspond dans le cas de Au et Cu à la composition stœchiométrique AuCu3. La

I-3-4-2-Composés intermétalliques ordonnés (Composés définis)

(hexagonal), MgCu2(cubique), CuAl2. ou MgNi2 (hexagonal) [10].

I-4- Structures et dénominations des phases métalliques

PDF created with pdfFactory Pro trial version www.pdffactory.com

dans l’une des structures A1, A2, ou A3.