Elaboration

CHAPITRE I : DESCRIPTION DU PROCEDE
D’ELABORATION SOUS MICRO-ONDES
-9-
Chapitre I : Description du procédé d’élaboration sous micro-ondes.
1. INTRODUCTION
La maîtrise de l’élaboration de matériaux composites à matrice polymère

fondée sur l’hystérésis diélectrique nécessite de connaître préalablement les
matériaux composites. Cela exige une parfaite maîtrise de l’application de l’énergie
micro-ondes. A la maîtrise de la propagation des ondes électromagnétiques s’ajoute
l’aspect cinétique chimique du procédé (le matériau utilisé est à matrice
thermodurcissable), qui a une influence significative sur le chauffage du matériau.
De plus, le procédé s’accompagne de l’application d’une contrainte de

pression, qui a pour rôle de mettre en forme la pièce en composite. Par conséquent,
une meilleure compréhension de la mise en forme des matériaux composites à
matrice polymère thermodurcissable est important pour l’optimisation du procédé.
Ainsi après un rappel sur les renforts et matrices des matériaux composites
ainsi que leurs propriétés, nous verrons la propagation des ondes
électromagnétiques dans le guide d’onde, le dispositif expérimental, suivi d’une
modélisation cinétique et enfin de la description du transfert de matière lors de la
mise en forme du matériau composite.
2. LES MATERIAUX COMPOSITES
Il existe deux types de matériaux composites :
• les composites grande diffusion (GD)

• les composites hautes performance (HP)
Les composites GD, représentent 95 % des matériaux composites utilisés. En

général, ce sont des plastiques armés ou des plastiques renforcés dont le taux
volumique de renfort avoisine 30 %. Les principaux constituants de base sont les
résines polyesters (95 % des résines thermodurcissables) avec des fibres de verre
(plus de 99 % des renforts utilisés).
- 11 -
Les composites HP sont principalement utilisés dans l’aéronautique. Les

renforts utilisés sont plutôt des fibres longues dont le taux de renfort est supérieur
à 50 %. Contrairement aux composites GD, leurs propriétés mécaniques, résistance
mécanique et rigidité, sont largement supérieurs à celles des métaux.
Les principaux constituants des matériaux composites sont résumés dans
l’organigramme figure 1-1.
La maîtrise de l’élaboration de matériaux composites nécessite la

connaissance des propriétés individuelles des constituants de base des matériaux.
Charges
Craie
Silice
Bille verre
Matrice/liant Renforts/fibres
Phénoliques Celluloses
MATERIAUX
Polyesters COMPOSITES Verre E, R…
A MATRICE Carbone/graphite
Epoxydes ORGANIQUE
Additifs
Catalyseur
Accélérateur
Figure 1 - 1 : Constituants des matériaux composites.
- 12 -
2.1. Les renforts
2.1.1. Généralités
Le renfort a pour rôle d’assurer la tenue mécanique du composite et est

souvent de nature filamentaire : fibres organiques ou inorganiques.
Les fibres sont constituées par plusieurs centaines ou milliers de filaments de

diamètres variant de 5 à 15 µm, ce qui permet de traiter les mèches sur des
machines textiles.
Les différents types de renforts sont cités dans l’organigramme ci-dessous.
Renfort
Organiques Minéraux
Artificiels Végétaux Céramiques Métalliques
Polyesters Aramides Bois Coton Verre Carbone Bore

Papier
Figure 1 - 2 : Les différents types de renfort de base.
- 13 -
Le tableau suivant présente les principales caractéristiques mécaniques de

quelques renforts [I-1] :
Fibre Densité Charge de Charge de Allongement Module Diamètre

(kg/dm3) rupture en rupture en à la rupture d’élasticité du filament
traction compression (%) (MPa) élémentaire
(MPa) (MPa) (µm)
Verre E 2.54 3400 1200 4.8 73000 3 - 30
Verre R 2.48 4400 1300 5.4 86000 3 - 30
Aramide 1.45 3100 500 1 130000 12
haut
module
Carbone 1.8 2200 1300 1.8 400000 7
haut
module
Bore 2.63 3500 3500 0.8 400000 100 - 200
Tableau 1 - 1 : Caractéristiques de quelques renforts.
2.1.2. Les fibres de verre
Le verre possède une très bonne résistance à la rupture dépassant même

celle de certains métaux. Cependant, le verre est un matériau très fragile du fait de
sa forte sensibilité aux microfissures. Cette fragilité diminue lorsque le verre est
sous forme de fibre puisque les microfissures sont moins nombreuses, ce qui lui
confère de bonnes performances.
Il existe différents types de verre : le verre E, C, D, R ou S, qui sont obtenus

suivant une composition différente du mélange de préparation, par conséquent,
leurs propriétés diffèrent également.
Le verre E est employé pour la fibre de verre constituant le renfort du

matériau composite faisant l’objet de ce travail.
- 14 -
2.2. Les matrices
La matrice permet de lier les fibres du renfort fibreux entre elles, ainsi que de
répartir les efforts (résistance à la compression ou à la flexion). La matrice est
facilement déformable et assure la protection chimique des fibres[I-2]. Généralement,
c’est un polymère ou une résine organique.
On peut retrouver dans l’organigramme ci-dessous, les différentes familles de
matrice :
Matrices
Organiques Minérales
Thermodurcissables Thermoplastiques Elastomères

Céramiques Métalliques
Borures Carbures Nitrures
Figure 1 - 3 : Différentes familles de matrices.
Les résines les plus employées dans les matériaux composites sont les
résines thermodurcissables et les résines thermoplastiques.
Les résines thermodurcissables sont des polymères, qui, après un traitement
thermique ou physico-chimique (catalyseur, durcisseur), se transforment en des
produits essentiellement infusibles et insolubles. Ces résines ont donc la
particularité de ne pouvoir être mises en forme qu’une seule fois.
Les résines thermoplastiques, en revanche, peuvent être alternativement ramollies
par chauffage et durcies par refroidissement dans un intervalle de température
spécifique du polymère étudié. De plus ces résines présentent l’aptitude à l’état
liquide de se mouler facilement par plasticité.
La matrice constituant le matériau composite de cette étude est une matrice
thermodurcissable.
- 15 -
Une comparaison entre quelques caractéristiques des résines

thermodurcissables et celles des résines thermoplastiques est présentée dans le
tableau 1-2 :
Critères essentiels Matrices thermoplastiques Matrices thermodurcissables

Etat de base Solide (prêt à l’emploi : Liquide visqueux à
polymérisé) polymériser
Stockage matière de base illimité Temps réduit (précautions à
prendre)
Mouillabilité des renforts difficile aisée
Moulage Chauffage (fusion + Chauffage continu
refroidissement de fixation)
Cycle plus long (refroidissement) court
Critères spécifiques
Tenue au choc assez bonne limitée
Tenue thermique réduite meilleure
Chute et déchets recyclables perdus
Tableau 1 - 2 : Caractéristiques des résines thermodurcissables et des résines

thermoplastiques.
Il existe d’autres catégories de résines : les résines thermostables et les

polymères thermoplastiques, dont l’usage reste spécifique et dont la fabrication
dans les matériaux composite reste limitée.
Les résines thermostables sont caractérisées par des caractéristiques
mécaniques stables sous des pressions et des températures élevées (> 200°C)
appliquées de façon continue. Cette propriété se mesure en déterminant la
température que peut supporter la résine durant 2000 h sans perdre la moitié de
ses caractéristiques mécaniques.
Les élastomères sont caractérisés par une forte élasticité avec un module
d’Young très faible.
2.3. Les matériaux composites structuraux
Les matériaux composites structuraux sont classés en trois catégories :
- 16 -
• les monocouches
• les stratifiés
• les sandwichs
Les monocouches représentent l’élément de base de la structure composite.

La forme du renfort caractérise les différents types de monocouches. Le renfort peut
être à fibres longues : ce sont les plis unidirectionnels (noté parfois UD), courtes ou
à fibres tissées.
Les stratifiés sont constitués d’un empilement de monocouches ayant

chacun une orientation propre par rapport à un référentiel commun aux couches et
désigné comme le référentiel du stratifié.
Figure 1 - 4 : Composite stratifié.
Le choix de l’empilement et plus particulièrement des orientations permettra

d’obtenir des propriétés mécaniques spécifiques.
Les stratifiés peuvent être de trois types :
1) Equilibrés : le stratifié contient autant de couches orientées suivant la
direction + θ que de couches orientées suivant la direction - θ
2) symétriques : les couches du stratifié sont disposées symétriquement
par rapport à un plan moyen
3) orthogonaux : le stratifié comporte autant de couches à 0° que de
couches à 90°
Les sandwichs sont des matériaux possédant deux peaux de grande rigidité
et de faible épaisseur et renfermant un cœur de forte épaisseur et de faible
- 17 -
résistance. L’ensemble forme une structure d’une grande légèreté. Le matériau

sandwich possède une grande légèreté en flexion et c’est un excellent isolant
thermique.
peaux
cœur
Figure 1 - 5 : Composite sandwich.
Les matériaux composites utilisés dans le cadre de ce travail sont

uniquement des monocouches à fibres de verre longues unidirectionnelles associées
à une matrice résine époxyde.
3. ELABORATION DES MATERIAUX COMPOSITES A MATRICE

POLYMERE
Les procédés de mise en œuvre de matériaux composites nécessitent trois

étapes indispensables :
1) une étape d’association : l’imprégnation du renfort par la résine

2) une étape rhéologique de mise en forme de la pièce
3) une étape thermique de solidification du matériau : soit par
refroidissement pour les matrices thermoplastiques, soit par
polycondensation ou réticulation pour les matrices
thermodurcissables (à température croissante)
Il existe différents procédés mais le plus utilisé est le procédé par moulage.
Le moulage peut s’effectuer dans un moule à surface non fermée tel que le cas du
moule ouvert (technique par contact, projection, ou pultrusion) ou à moule fermé
(par exemple, injection de résine thermodurcissable associée à des fibres longues
suivi d’une compression). Les principales opérations d’une mise en forme par
moulage sont rappelées (figure 1-6).
- 18 -
PROCEDE
D’ELABORATION
Résine Renfort
Imprégnation
(mélange)
Mise en place du mélange sur

l’outillage
Compactage
Polymérisation ou
réticulation ou
solidification par
refroidissement
Démoulage
Finition
Figure 1 - 6 : Principales étapes d’un procédé par moulage.
- 19 -
La constitution du moule va dépendre du choix de la technique employée. Le

choix de la technique à utiliser repose sur plusieurs paramètres dont voici quelques
exemples :
! la forme et la taille des pièces à élaborer, car certains procédés se limitent

à l’élaboration de pièces à géométrie simple uniquement.
! la nature de la matrice du matériau composite, qui doit être adaptée à la

technique employée. Par exemple la résine thermodurcissable nécessite
un outillage spécifique puisqu’elle implique des températures et une
pression précises.
! le nombre de pièces à élaborer, qui détermine le coût de l’outillage. Pour

une production en grande série (plusieurs milliers de pièces)
généralement un moule métallique pouvant supporter de très nombreux
cycles est employé. Pour des dizaines de pièces, un moule peu coûteux
(en plâtre par exemple) peut convenir.
! la qualité mécanique des pièces : les pièces à fort taux de fibre ne

s’obtiennent qu’avec certains procédés.
3.1. Procédés classiques et procédés micro-ondes
Il existe une multitude de procédés classiques donc nous citerons

uniquement des procédés le plus couramment employés.
3.1.1. Procédés classiques
Le problème qui nous concerne est la mise en place de la matrice, qui

conduit globalement à une mise en forme de la pièce. Cette mise en forme est
tributaire de la distribution spatiale particulière de la température résultant des
sources micro-ondes. Pour cette raison, nous avons classé les procédés en deux
catégories suivant des critères de température et de contrainte de pression
mécanique appliquée lors de la mise en œuvre.
- 20 -
Contrairement aux procédés à faible pression, les procédés à forte pression

sont très mécanisés et sont utilisés pour les grandes séries.
De plus dans la première catégorie, la résine est mélangée au renfort pendant le
moulage tandis que les matériaux utilisés pour les procédés de la deuxième
catégorie sont déjà sous forme de préimprégnés.
1) Les procédés à faible pression à froid :
! Moulage au contact : le principe consiste à imprégner manuellement les

renforts disposés dans un moule. Les résines employées sont rarement
des résines époxy mais plutôt des résines polyesters. Cette technique est
peu onéreuse et des pièces de formes quelconques peuvent être réalisées
mais à cadence très faible.
! Moulage par projection : technologie similaire au moulage par contact

mais les fibres coupées sont projetées au pistolet cette fois. La limite du
procédé réside dans la qualité moyenne des pièces obtenues.
2) Les procédés à haute pression à chaud :
! Injection de thermodurcissables BMC (Bulk Molding Compound ou

préimprégné en vrac). Procédé discontinu haute pression (de l’ordre de
100 bars). Les principales étapes sont : l’alimentation et dosage du
compound, injection-pression, maintien et réticulation, puis éjection. Les
avantages d’un tel procédé sont : la réalisation de grande série, le faible
coût de la matière, le peu de finition. Les limites sont : la longueur et le
taux des renforts et les propriétés mécaniques du composite obtenu.
! Compression de thermodurcissables SMC (Sheet Molding Compound ou

préimprégnés en feuilles). Le principe consiste à déposer des feuilles de
préimprégnés dans un contre-moule chauffé, de comprimer le matériau
avec un moule chauffé, réticulation puis éjection de la pièce. Les
avantages sont : le coût de la matière, les propriétés thermiques et
mécaniques. Les limites sont : l’aspect de la pièce obtenue, le
dimensionnement des presses, et la finition.
- 21 -
! Pultrusion : est utilisée pour réaliser des matériaux composites à hautes

performances industrielles. Les principales étapes sont : le tirage, la mise
en forme et polymérisation de la matrice d’imprégnation. Les avantages
sont : la production en continu, la possibilité de réaliser des pièces à
section complexe, et d’avoir un taux de renfort élevé. Les limites sont la
lenteur du procédé, et l’obtention, dans la très grande majorité des cas,
de profilés droits à section constante.
! Enroulement filamentaire (ou bobinage) : technologie concernant les

matériaux composites hautes performances. Le principe consiste à
enrouler sous tension, des fibres continues préalablement imprégnées
d’un liant, sur un mandrin tournant autour de son axe. Les avantages
sont la disposition optimale des renforts dans le sens des efforts
appliqués, les très bonnes propriétés mécaniques du fait d’un taux de
fibres élevé (60 à 75 % en masse) et la possibilité de réaliser des pièces de
grandes dimensions avec des surfaces internes lisses. Les limites sont
l’obtention de pièces à formes uniquement convexes, les investissements
importants.
Nous n’avons cité que quelques procédés jugés importants, mais bien sûr il en
existe d’autres tels que : le drapage en autoclave, le thermoformage sous vide ou
encore l’estampage…
3.1.2. Procédés d’élaboration à chauffage par hystérésis

diélectrique
Les avantages (développés ultérieurement) offerts par le procédé micro-ondes,

tels que la rapidité du chauffage en volume, a permis le développement de cette
technologie pour les traitements des matériaux diélectriques, souvent mauvais
conducteurs de la chaleur mais très bons capteurs de l’énergie électromagnétique.
Cette dernière caractéristique est primordiale pour l’efficacité du procédé. C’est
pourquoi la vulcanisation du caoutchouc est le secteur industriel qui a le plus
bénéficié de la technique micro-ondes.
- 22 -
En effet, le caoutchouc est un matériau macromoléculaire comportant des

charges de carbone, ce qui lui confère une très fort coefficient d’absorption
diélectrique, le rendant très sensible à l’énergie électromagnétique.
Si des articles ou brevets [I-3, I-4, I-5, I-6, I-7, I-8, I-9, I-10, I-11, I-12] traitent des procédés
de transformation de matériaux macromoléculaires sous micro-ondes, peu
prennent en compte le problème de fluage de la matrice polymère à travers le
renfort fibreux, puisque l’aspect thermique est toujours privilégié. L’objectif pour
nous est d’étudier le procédé d’élaboration par hystérésis diélectrique de matériaux
composites à matrice polymère thermodurcissable renforcés en fibres de verre, tant
sur le plan de la thermique micro-ondes que sur l’aspect écoulement de la matrice
au cours de la mise en forme de l’objet.
C’est pourquoi notre travail est une contribution à la compréhension des
phénomènes thermiques et rhéologiques mis en jeu ainsi que les phénomènes liés à
l’écoulement de la matière, ces deux types de phénomènes physiques étant couplés.
3.2. Couplages mis en œuvre au cours du procédé fondé sur

l’hystérésis diélectrique
L’élaboration de matériaux composites verre/époxy avec chauffage par

hystérésis diélectrique engendre un couplage entre le champ électromagnétique, la
réaction physico-chimique, et la température et la viscosité du matériau traité
(figure 1-7).
3.2.1. Couplages entre le champ électromagnétique, la

température, et les réactions physico-chimiques au sein du
matériau
La température joue un rôle important au cours du procédé. De plus elle

intervient dans le terme source micro-ondes de l’équation de la chaleur par
l’intermédiaire des caractéristiques diélectriques, ainsi que par l’amplitude du
champ électrique qui se trouve être solution de l’équation de propagation des ondes
(dans laquelle ces caractéristiques interviennent).
- 23 -
Les distributions spatiale et temporelle de la température sont solutions de

l’équation de la chaleur composée d’un terme source micro-ondes et d’un terme
source chimique liée à l’exothermicité de la réaction époxy/amine. La température,
agissant sur les caractéristiques diélectriques du milieu, intervient sur le terme
source micro-ondes dépendant de l’amplitude du champ électrique, solution de
l’équation des ondes dans laquelle interviennent les caractéristiques diélectriques
du milieu. De plus les paramètres thermophysiques tels que la conductivité
thermique k, la chaleur spécifique massique Cp, présents dans l’équation de la
chaleur, dépendent de la température. Il en est de même pour les autres
paramètres physiques intervenant ou non dans l’équation physico-chimique comme
la conversion chimique x, ou la viscosité η, traduisant la réponse du milieu à
l’augmentation de la température, telle que la réaction chimique ou le changement
d’état.
- 24 -
Terme source Terme source

électromagnétique physicochimique
Equation de
l’énergie
k, ρ, Cp
Température Conditions
aux limites
thermiques
ε′, ε″, µ x, η
Equations de Maxwell Equations physico-chimiques
r r
E et H dans le matériau Viscosité, conversion
Conditions aux limites Conditions aux limites

électromagnétiques physico-chimiques
Figure 1 - 7 : Couplages lors de la mise en œuvre.
3.2.2. Couplages entre le matériau et son environnement

thermique et physique
La température au sein du produit chauffé est obtenue par la résolution de

l’équation de la chaleur à laquelle il faut ajouter les conditions aux limites. Ces
conditions aux limites peuvent être des conditions d’échanges thermiques entre la
- 25 -
surface du matériau et son environnement (refroidissement en surface par

exemple). Ainsi la température au sein du matériau résulte des transferts
thermiques surfaciques liés aux conditions thermiques de l’environnement selon les
lois de la diffusion de la chaleur.
La diffusion de la chaleur est un phénomène lent, contrairement à la

propagation électromagnétique. De plus, quelques microsecondes suffisent à
l’établissement de l’équilibre électrique de la relaxation diélectrique contrairement à
l’équilibre thermique (quelques minutes à quelques heures suivant le matériau
considéré et ses dimensions). Tout ceci nous amène à penser que malgré le
couplage entre les phénomènes thermiques et électromagnétiques, ils ne présentent
pas du tout la même évolution dans le temps.
C’est la raison pour laquelle, comme nous verrons par la suite, l’équation de la
chaleur ainsi que l’équation de la cinétique chimique (représentant le phénomène
physico-chimique) seront résolus en régime transitoire tandis que l’équation des
ondes sera résolue en régime permanent.
La détermination du champ de température au sein du matériau élaboré

sous micro-ondes nécessite la résolution de plusieurs équations :
! équations des ondes : équations de Maxwell

! équation de la cinétique chimique de la matrice thermodurcissable choisie
! équation de la chaleur
En effet, la résolution des équations de Maxwell pour des grandeurs physiques

(caractéristiques diélectriques et profil thermique) du milieu fixées initialement
permet d’accéder à la distribution du champ électrique au sein de l’applicateur. La
connaissance de ce champ permet, par la suite, de calculer le terme source micro-
ondes présent dans l’équation de la conduction. Puis s’ajoute le calcul du terme
source d’origine chimique, qui s’obtient en résolvant l’équation de la cinétique
chimique.
La résolution de l’équation de conduction de la chaleur va ainsi conduire à une
nouvelle cartographie de la température. Cette nouvelle distribution de la
température modifie les paramètres physiques initiaux du matériau chauffé,
dépendant de la température. Ainsi un nouveau calcul du champ électrique,
dépendant des grandeurs physiques du milieu, s’impose. Nous sommes donc en
- 26 -
présence d’une résolution en boucle imposée par la nécessité d’actualiser à chaque

itération les grandeurs physiques thermodépendantes, dont est fonction le champ
électrique.
4. CHAUFFAGE PAR HYSTERESIS DIELECTRIQUE
Dans le spectre des ondes électromagnétiques, les micro-ondes, appelées

aussi hyperfréquences, se situent à la frontière entre les ondes électriques et les
rayons optiques. Plus précisément, elles se trouvent entre les ondes hertziennes
utilisées en radiocommunication et le spectre infrarouge. Leurs longueurs d’ondes
varient du millimètre au mètre et les fréquences de 300 MHz à 300 GHz.
Pour éviter toute interférence avec les systèmes de télécommunication, des

fréquences spécifiques de 915 et 2 450 MHz ont été attribuées par la législation
mondiale aux applications micro-ondes. La fréquence fixe utilisée tout au long de ce
travail sera égale à 2 450 MHz soit 2.45 GHz.
4.1. Propagation électromagnétique
La propagation électromagnétique étant régie par les équations de Maxwell,

nous ferons donc un bref rappel de ces équations. Nous verrons ensuite les
différences entre la propagation en espace libre et la propagation guidée ainsi que
les conditions d’existence des modes de propagation.
4.1.1. Equations de Maxwell
Quatre quantités vectorielles décrivent les phénomènes électromagnétiques :
r
1) le champ électrique E (volts par mètre)
r
2) l’induction électrique D (coulombs par mètre carré)
r
3) le champ magnétique H (ampères par mètre)
r
4) l’induction magnétique B (webers par mètre carré)
- 27 -
Ces vecteurs sont dépendants et, pour un milieu diélectrique à pertes, de

1
permittivité complexe ε% (pour le vide ε0 = F /m ) et de perméabilité
36 π ⋅ 109
magnétique complexe µ% (pour le vide µ0 = 4 π ⋅ 10−7 H/m ), ils sont liés par les
relations suivantes :
r r
• D = ε% E (1.1)
r r
• B = µ% H (1.2)
Les équations de Maxwell, dont les champs électromagnétiques sont

solutions, s’écrivent sous forme temporelle en tenant compte des équations (1.1) et
(1.2) :
r r
r r ∂B ∂H
rot E = − = −µ% (1.3)
∂t ∂t
r r
r r r ∂D r ∂E
rot H = J + = J + ε% (1.4)
∂t ∂t
r
div D = ρe (1.5)
r
div B = 0 (1.6)
r
où J est le vecteur densité de courant.
ρe est la densité de charges électriques libres C/m3 ( )
Les équations fondamentales (1.3) et (1.4) résultent respectivement de la

relation de Maxwell-Faraday et de la relation de Maxwell-Ampère. Les deux
équations en divergence (1.5) et (1.6) représentent respectivement le théorème de
Gauss et la caractéristique de l’induction magnétique (soumise à l’impossibilité de
l’existence de charges magnétiques libres).
- 28 -
4.1.2. Équations de la propagation libre
Nous nous intéressons uniquement aux diélectriques non chargés

électriquement ( ρe = 0 ). Les équations de Maxwell, sous forme temporelle et en
l’absence de densité de courant, sont les suivantes :
uuuuuu uv
uvv ∂H
rot E = −µ% (1.7)
∂t
uuuuuu uv
uvv ∂E
rot H = ε% (1.8)
∂t
uv uv
div D = 0 soit div E = 0 si ε% est indépendant de (x, y, z) (1.9)
uv uv
div B = 0 soit div H = 0 si µ% est indépendant de (x, y, z) (1.10)
Ces équations représentent le cas le plus général de la propagation.
Les expressions des champs électrique et magnétique en régime harmonique

(sinusoïdal) sont :
uv uv
E ( x, y, z,t ) = E ( x, y, z ) e jωt (1.11)
uv uv
H ( x, y,z,t ) = H ( x, y, z ) e jωt (1.12)
Les équations de Maxwell en régime harmonique, dans les mêmes conditions que
précédemment, s’écrivent :
uuuuuu
uvv uv
rot E = − jωµ% H (1.13)
uuuuuu
uvv uv
rot H = jωε% E (1.14)
uv
div E = 0 (1.15)
uv
div H = 0 (1.16)
Dans le cas d’une propagation dans le vide, les valeurs de permittivité et de

perméabilité sont celles du vide : µ = µ0 et ε = ε0
- 29 -
Les équations de propagation s’obtiennent en prenant le rotationnel des

équations (1.3) et (1.4) et en tenant compte des relations suivantes :
uuuuuuuuuuu
uvvu uuuuuuuuuuuuuuu
uuuuuu uvuvu uuuuvv uuuuvv uv
( )
rot rot E = grad div E − ∆E = −∆E car div E = 0 d’après l’équation (1.15)
uuuuuuuuuuu
uvuv uuuuuuuuuuuuuuu
uuuuuu uvuvu uuuv
uv uuuv
uv uv
( )
rot rot H = grad div H − ∆H = −∆H car div H = 0 d’après l’équation (1.16)
Ainsi pour un régime temporel nous obtenons deux équations de propagation

uv uv
en E et H :
uv
uv ∂2 E
∆E − µ% ε% 2 = 0 (1.17)
∂t
uv
uv ∂2 H
∆H − µ ε 2 = 0
% % (1.18)
∂t
et pour un régime harmonique ces équations prennent les formes suivantes :
uuuuvv uv
% 2E = 0
∆E − µ% εω (1.19)
uuuv
uv uv
% 2H = 0
∆H − µ% εω (1.20)
où ∆ est le laplacien vecteur.
Sachant que la constante de propagation complexe k% est définie telle que

k% = ω ε% µ% , les précédentes équations de propagation s’écrivent alors :
uv uv
∆E + k% 2 E = 0 (1.21)
uv uv
∆H + k% 2 H = 0 (1.22)
Pour un milieu sans pertes diélectriques et magnétiques on a :
1
v= = λf (1.23)
ε0 ε r µ 0µ r
ω 2πf 2π
k= = = (1.24)
v v λ
- 30 -
avec λ la longueur d’onde, v la vitesse de propagation de l’onde, f la fréquence, ε r
et µ r respectivement la permittivité et la perméabilité relatives du milieu.
4.1.3. La propagation guidée
En propagation guidée, les ondes se propagent dans des structures

métalliques guidantes appelées guide d’ondes. En propagation libre, l’onde se
propage dans l’espace libre non limité, ce qui entraîne une dispersion de l’énergie
donc une impossibilité de mise en œuvre de traitements thermiques de matériaux
dans ce cas. L’utilisation d’un guide d’onde lors de l’élaboration du matériau est
donc d’usage car, contrairement à la propagation libre, en propagation guidée il n’y
a pas de dispersion d’énergie sur les surfaces d’ondes, ainsi le bilan d’énergie est
supérieur à celui obtenu en propagation libre. Cette augmentation du rendement
énergétique provient du fait que l’énergie est confinée à l’intérieur des limites des
« structures guidantes » en question, plus couramment appelées « applicateurs
d’énergie » ou simplement « applicateurs ».
Un « applicateur » se compose d’une enveloppe métallique qui entoure

généralement un remplissage diélectrique et constitue une structure dans laquelle
les ondes électromagnétiques se propagent. Par simplification mathématique, ne
seront prises en compte que les propagations suivant la direction Oz, orthogonale à
la section transversale de l’applicateur, que nous supposerons de forme
rectangulaire, comme les guides d’onde utilisés dans notre cas et qui sont de
section rectangulaire.
Les conditions aux limites sont essentielles et définissent les conditions

uv uv
d’existence des champs E et H sur les parois métalliques du guide. Par définition
les conditions aux limites à la surface d’un métal, supposé conducteur parfait,
sont :
! l’existence d’un champ électrique normal ou nul aux parois

! l’existence d’un champ magnétique tangentiel aux parois de l’applicateur
- 31 -
4.1.3.1. Constante de propagation guidée
Les équations (1.21) et (1.22) donnent les équations d’ondes. Pour que les
solutions harmoniques restent valides, le coefficient k% 2 doit rester positif, sinon les
solutions ne seraient plus sinusoïdales mais décroissantes exponentiellement.
Dans une structure à propagation guidée les composantes des champs électrique et
magnétique se décomposent en deux parties longitudinale et transversale suivant :
uv uv uv
E = E long + Etrans (1.25)
uv uv uv
H = Hlong + Htrans (1.26)
L’équation différentielle de la propagation se décompose donc en deux équations

différentielles dont l’une correspond à la composante longitudinale et l’autre à la
composante transversale.
Lors d’une propagation dans un applicateur rectangulaire, l’équation d’onde

donnée équation (1.21) devient dans un système de coordonnées cartésiennes :
∂ 2E Z ∂ 2E Z ∂ 2E Z
+ + + k 2E z = 0 (1.27)
∂x 2
∂y 2
∂z 2
Une résolution analytique à l’aide de la méthode de séparation des variables[I-13],

permet de déterminer les solutions des équations de propagation. D’après cette
méthode, les solutions de l’équation (1.21) sont de la forme :
Ez = A ( x ) B ( y ) C ( z ) (1.28)
en substituant l’équation (1.28) dans l’équation (1.27), puis en divisant par E z on

obtient :
1 ∂ 2 A 1 ∂ 2B 1 ∂ 2C
+ + + k2 = 0 (1.29)
A ∂x 2 B ∂y 2 C ∂z 2
- 32 -
Sachant que les trois premiers termes de l’équation ne possèdent qu’une seule
coordonnée et que chaque coordonnée est indépendante des deux autres alors on
peut poser les relations suivantes :
1 ∂2 A
= −k 2x (1.30)
A ∂x 2
1 ∂ 2B
= −k 2y (1.31)
B ∂y 2
1 ∂ 2C
= −k 2z (1.32)
C ∂z 2
où k x , k y , et k z sont des constantes indépendantes des coordonnées x, y, z du
point considéré.
Nous obtenons par la suite les trois équations à variables séparées :
∂2 A
+ k 2x A = 0 (1.33)
∂x 2
∂ 2B
+ k 2y B = 0 (1.34)
∂y 2
∂ 2C
+ k 2zC = 0 (1.35)
∂z 2
dont les constantes doivent satisfaire l’équation de séparation :
k 2x + k 2y + k 2z = k 2 (1.36)
En posant k 2c = k 2x + k 2y , traduisant la propagation dans le plan transversal, on
retrouve la définition de la constante de propagation guidée :
k 2 = k 2c + k 2g (1.37)
avec k c la constante de coupure, et k g la constante de propagation guidée (elle-
même égale à k z ), telle que pour k = k c la constante de propagation longitudinale
k g est nulle.
- 33 -
La propagation guidée est donc la superposition d’une propagation

longitudinale de constante k g suivant l’axe z, et d’une propagation transversale de
constante k c dans le plan (x, 0, y).
4.1.3.2. Conditions de propagation guidée
Les solutions de l’équation de propagation dépendent des conditions aux

limites imposées dans les différentes directions du système. Il existe ainsi différents
types de solutions correspondant à des situations physiques différentes. Ces
différentes solutions dépendent essentiellement de la nature de la constante de
propagation k% , qu’elle soit réelle, imaginaire ou complexe. L’argument de ces
solutions est complexe jk% , et se décompose en deux parties de la manière
suivante[I-14] :
! la partie réelle α correspond à l’atténuation de l’amplitude de l’onde au

cours de la propagation au sein du milieu
! la partie imaginaire β marque le caractère propagatif de la solution
Nous en déduisons alors la constante de propagation longitudinale :
γ = α + jβ (1.38)
Par définition, l’onde se propage suivant Oz avec une constante k g , donc
d’après l’équation (1.37) son expression est :
k 2g = k 2 − k 2c (1.39)
Suivant les valeurs de kc , différents types de comportement sont

envisageables[I-15].
Les trois cas de figure suivants peuvent se présenter :
1) k > k c , alors k g = k 2 − k 2c est un nombre réel.
- 34 -
Les solutions sont purement propagatives, c’est-à-dire que les ondes se

propagent dans le guide avec une constante k g sans subir d’atténuation.
2) k = k c , alors k g = 0
La solution est à la coupure, c’est-à-dire que le régime de propagation est

nommé régime de coupure. De ce fait, il n’y a aucune propagation ou
atténuation suivant Oz. Les champs vibrant en phase suivant cette direction,
la propagation se fait uniquement dans le plan transversal du guide.
3) k < k c , alors k g = ± j k 2 − k c2 est un nombre imaginaire pur.
La solution est une fonction exponentielle décroissante en exp ( −αz ) , c’est-à-
dire que les ondes s’atténuent exponentiellement en fonction de la distance z

parcourue dans l’applicateur. Les ondes sont appelées des ondes
évanescentes. Il s’agit du phénomène à l’origine de la mauvaise réception des
ondes de basse fréquence (radio ou téléphoniques par exemple) dans les
tunnels.
En somme, les différents types de propagation des ondes que l’on peut
rencontrer sont les ondes purement propagatives ou propagatives avec atténuation
due aux pertes diélectriques, ou stationnaires ou évanescentes. Ces types de
propagation dépendent d’une fréquence de coupure notée fc définie par :
kc v
fc = v = (1.40)
2π λ c
avec v la vitesse de l’onde en espace libre et λ c la longueur d’onde de coupure.
Une onde électromagnétique se propage dans un guide d’onde seulement si

k 2g > 0 , donc si k g est réel, ce qui revient à dire si sa fréquence d’émission, dans
notre cas : la fréquence du rayonnement micro-ondes, est supérieure à la fréquence

de coupure fc . En raisonnant sur les longueurs d’ondes, cela revient à dire qu’il y a
propagation si la longueur d’onde d’émission est inférieure à la longueur d’onde de
coupure.
2π
La longueur d’onde est définie par : λ g = .
kg
- 35 -
Elle peut aussi être déterminée par l’équation (1.37), dont l’équation équivalente en
longueur d’onde est l’équation ci-dessous :
1 1 1
= 2 + 2 (1.41)
λ0 λc λg
2
1 1 1
d’où l’expression finale : = 2 − 2 (1.42)
λ g λ0 λc
2
dans la situation la plus simple.
De même, on remarque que la constante de coupure élevée au carré est

positive, donc d’après la relation suivante :
k 2c = k 2 − k 2g (1.43)
Pour que ce soit le cas, la constante de propagation en espace libre doit être
supérieure à celle dans un applicateur, ce qui équivaut au fait que la longueur
d’onde guidée dans un applicateur est toujours supérieure à la longueur d’onde en
espace libre.
Comme la fréquence de coupure, la longueur d’onde de propagation dépend

des dimensions de la section droite du guide, du matériau diélectrique remplissant
1
le guide (c’est-à-dire de la vitesse des ondes puisque v = ), enfin de la fréquence
εµ
de l’onde mise en œuvre.
4.1.3.3. Modes de propagation
Lorsqu’une onde se propage au sein d’un guide d’ondes, les multiples

réflexions d’ondes sur les parois internes entraînent des distributions particulières
des champs électromagnétiques. Ces distributions des champs résultant de la
nécessité de satisfaire les conditions aux limites sur les parois du guide sont
appelées modes de propagation.
- 36 -
La propagation en espace libre s’effectue suivant le mode transverse

électromagnétique noté TEM, ce qui signifie avec des vecteurs champs électriques et
magnétiques uniquement transverses.
La propagation d’une onde dans un guide, en revanche, peut s’effectuer de
trois façons :
! seul le champ électrique est transverse, ces modes sont appelés
transverses électriques notés TE, définis tels que E z = 0
! seul le champ magnétique est transverse, ces modes sont appelés
transverses magnétiques notés TM, définis tels que Hz = 0
! ni le champ électrique ni le champ magnétique ne sont transverses, ces
modes sont appelés modes hybrides notés HEM
Suivant le système de coordonnées associé à la géométrie du milieu, ainsi

que les conditions aux limites imposées par la structure guidante, les équations de
uv uv
propagation en E et H peuvent être résolues. Les solutions obtenues sont des
solutions discrètes, dépendant chacune des couples de valeurs entières m et n qui
les caractérisent.
Les distributions des champs électromagnétiques seront alors décrites selon les
modes TE m,n , TMm,n ou HMm,n . Les indices m et n indiquent respectivement le
uv uv
nombre de demi-sinusoïdes que forment les champs E et H , le long du grand et du
petit côté (m pour x et n pour y). Ces modes sont caractérisés par une fréquence de
coupure (fc )m,n définie équation (1.40).
Les modes de propagation permettent de véhiculer dans l’applicateur la

puissance délivrée par le générateur. L’existence ou non d’un mode de propagation
dans un guide va dépendre des conditions de propagation suivantes :
! si f0 > ( fc )m,n ou λ0 < ( λc )m,n : propagation des modes TE m,n et TMm,n
! si f0 < ( fc )m,n ou λ0 > ( λc )m,n : pas de propagation des modes TE m,n et
TMm,n
avec f0 et λ 0 respectivement la fréquence et la longueur d’onde de travail, c’est-à-

dire du rayonnement émis par le générateur micro-ondes.
- 37 -
a z
y
b
Figure 1 - 8 : Applicateur à section rectangulaire.
Pour un applicateur rectangulaire (figure 1-8), le mode TE0,1 est appelé mode
fondamental car il possède la plus basse fréquence de coupure. (figure 1-9).

Pour un applicateur vide ou chargé par un seul matériau diélectrique de
permittivité ε , les composantes électriques pour le mode TE0,1 se réduisent à la
composante E x : fonction sinusoïdale de y.
La fréquence de coupure pour un mode TE m,n s’écrit :
v n2 m 2
(fc )TE = + (1.44)
m,n
2 a 2 b2
donc pour le mode TE0,1 cette expression prend la forme suivante :
v
(fc )TE = (1.45)
0,1
2a
car m = 0 et n = 1
Ex
a
y
b
Figure 1 - 9 : Distribution du champ électrique pour le mode TE0,1.
- 38 -
4.2. Le dispositif micro-ondes et sa conception
4.2.1. Equipements micro-ondes
Trois éléments principaux constituent le dispositif micro-ondes :
! le générateur équipé d’un tube magnétron qui transforme l’énergie

électrique en énergie micro-ondes
! les éléments d’adaptation
! l’applicateur, où le matériau à traiter est placé
Ces éléments sont représentés figure 1-10.
Adaptation
Impedances
d’impédance
adaptation y 1 2 3
circulateur
circulator z
Production
Microwaves
micro-ondes
uv
(EE . 5
production in 4 4
en TE
mode TE0,1
mode 1 2
01
#diélectrique silicone/verre %adaptations diélectriques en biseaux

$diélectrique silicone/mica &absorption totale de l’onde
'matériau composite verre/époxy
Figure 1 - 10 : Représentation conceptuelle du dispositif.
4.2.1.1. Générateur
Le constituant principal du générateur micro-ondes est le magnétron. Le rôle

du magnétron est de produire une tension électrique alternative à la fréquence de
2 450 MHz servant à alimenter une antenne. Le magnétron, en général de structure
cylindrique, est un tube constitué comme une diode classique. Il comporte une
cathode centrale, qui émet des électrons, et une anode concentrique externe percée
de cavités résonnantes.
- 39 -
Un champ électrique de direction radiale est appliqué entre l’anode et la

cathode. Un champ magnétique uniforme de direction parallèle à l’axe du tube est
créé, par des aimants ou bobines. La cathode émet des électrons qui sont accélérés
par le champ électrique suivant des trajectoires radiales vers l’anode. L’association
des deux champs électrique et magnétique confère aux électrons un mouvement de
giration. Ce mouvement provoque la résonance électrique des cavités, et comme les
cavités de l’anode sont couplées, cette résonance est transmise à une antenne
placée dans l’une des cavités, et reliée au milieu extérieur, ce qui permet d’émettre
la puissance électromagnétique.
La fréquence et la puissance émises par l’antenne du magnétron dépendent

de paramètres dimensionnels, de paramètres électriques et de paramètres
magnétiques. Ainsi un générateur micro-ondes est caractérisé par la puissance
émise par l’antenne du magnétron et la fréquence de l’onde émise ainsi que la
puissance qu’elle transporte suivant le régime de propagation dans le volume
entourant l’antenne.
4.2.1.2. Eléments d’adaptation et de mesure
4.2.1.2.1. Guide d’ondes rectangulaire
Le rôle du guide d’ondes est d’assurer le transport des ondes. Le guide utilisé
est un guide standard WR340 dont les dimensions de la section transversale sont
les suivantes : 43.18 mm pour le petit côté noté a, et 86.36 mm pour le grand côté
noté b, section telle que seul le mode fondamental TE01 se propage à la fréquence de
2.45 GHz. Pour un tel guide on a λ g = 170 mm environ.
4.2.1.2.2. L’isolateur
L’isolateur est l’ensemble formé par les deux éléments suivants : le

circulateur et la charge à eau.
- 40 -
Le circulateur est placé entre la source micro-ondes et le circuit d’utilisation,

et sert à protéger le magnétron des réflexions d’ondes, tout en ne perturbant pas
l’onde incidente émise par le magnétron en direction du matériau à traiter. Il est
constitué d’un noyau de ferrite à aimantation permanente.
Le circulateur, dont la structure est en forme d’un té, sert de sens giratoire
pour l’onde, c’est-à-dire que l’onde incidente située dans une des trois voies se
retrouve dans une seule des deux autres voies restantes.
La charge à eau est utilisée dans le but d’assurer un régime d’ondes

progressives au sein de l’applicateur micro-ondes. Son rôle est d’empêcher une
réflexion de l’onde, donc l’établissement d’un régime d’onde stationnaire, par une
absorption totale d’une éventuelle énergie résiduelle en bout d’applicateur. C’est
pourquoi souvent les applicateurs se terminent par une charge à eau.
Généralement, l’association du circulateur à une seconde charge absorbante,
telle qu’une charge à eau, est utilisée dans le cas où le matériau présente une trop
forte réflexion. De cette façon, l’onde réfléchie vers le générateur est déviée par la
ferrite puis absorbée totalement par la charge à eau.
4.2.1.2.3. Le coupleur directif
Le coupleur directif permet la mesure précise de la puissance véhiculée par

l’onde dans un circuit micro-ondes. La méthode consiste à dévier une faible partie
de la puissance électromagnétique par prélèvement de celle-ci par une antenne
reliée à un bolomètre, lui-même relié à un wattmètre. Par exemple pour mesurer la
puissance véhiculée par l’onde arrivant en 1 et sortant en 2, une très faible partie
de cette puissance (1/1 000 000 précisément) sera prélevée par l’antenne puis
mesurée à l’aide du wattmètre.
4.2.1.3. L’applicateur
L’applicateur est l’élément essentiel du dispositif car il permet le transfert de

l’énergie électromagnétique provenant du générateur au sein du matériau à traiter.
Sa conception dépend non seulement de la nature, de la forme, et des dimensions
- 41 -
des produits à traiter mais aussi de la fréquence de travail ainsi que de la

puissance mise en œuvre.
Il existe deux types d’applicateurs : les applicateurs monomodes et les

applicateurs multimodes.
Un applicateur monomode est défini par un seul mode de propagation transportant
la puissance. Un applicateur multimode, comme son nom l’indique, est caractérisé
par un nombre de modes de propagation important. L’apparition de plusieurs
modes de propagation provient de l’augmentation des dimensions de l’applicateur.
De ce fait, la propagation des ondes devient complexe et variée. Généralement, les
enceintes multimodes sont employées pour les grands volumes, dans le but
d’obtenir une certaine distribution des points à champ fort.
La propagation des ondes peut s’y faire suivant différents régimes :
! Régime d’ondes stationnaires : ce régime est défini par des nœuds (points
à champ nul) ou ventres (zones à champ fort) d’énergie séparés d’une
longueur constante et égale à λ g /4 . Le matériau à élaborer est placé
dans un ventre d’énergie où le champ électrique est maximal. Ce mode de

traitement est tout à fait adapté aux objets de faible volume.
! Régime d’ondes progressives : ce régime est caractérisé par une charge

terminant le circuit micro-ondes empêchant toute réflexion d’ondes. Dans
ce cas, la puissance véhiculée par l’onde le long de l’applicateur contenant
le produit à traiter décroît exponentiellement du fait de l’absorption de
l’énergie par le produit.
4.2.2. Conception du dispositif
Le but recherché est l’obtention d’une pièce en composite chauffée de façon

uniforme et en un temps réduit comparé au procédé conventionnel. Nous verrons
comment et par quels moyens expérimentaux il est possible d’atteindre cette
uniformité finale du champ de température.
- 42 -
4.2.2.1. Sélection du mode de propagation fondamental
Le mode fondamental TE 0,1 d’un applicateur parallélépipédique donne les
volumes à distribution spatiale du champ électrique les plus significatifs et les plus
uniformes. En effet, si nous distinguons les trois directions de l’espace, dans la
direction 0x correspondant à l’indice 0 du mode, le vecteur champ électrique,
possède une amplitude constante. Dans la direction y, l’amplitude du champ
électrique suit une loi en sinusoïdale à un seul maximum[I-16]. Enfin, dans la
dernière direction z, l’amplitude du champ électrique suit généralement une
atténuation exponentielle du fait de l’absorption d’énergie dans la direction de
propagation comme nous le verrons en détails ultérieurement. Pour maintenir ce
mode fondamental au sein de l’applicateur et éviter que les modes d’ordre supérieur
ne soient excités, le moule est constitué de matériaux diélectriques spécifiques dont
les propriétés sont spécifiées par la suite.
4.2.2.2. Chargement diélectrique de l’applicateur
L’applicateur est rempli d’un moule dont le diélectrique central est un

composite silicone/mica (car de permittivité voisine du matériau à traiter, ayant
une excellente tenue à la chaleur, ainsi qu’une grande résistance mécanique), et le
diélectrique extérieur un composite silicone/verre. Pour un meilleur traitement, on
place le matériau au milieu de l’applicateur où le champ électrique est maximum.
d I
c II
III e
y IV
c
d V
z
Figure 1 - 11 : Remplissage du moule diélectrique.
- 43 -
Les couches d’épaisseur « d » sont des matériaux diélectriques silicone/verre,

les couches d’épaisseur « c » sont des diélectriques silicone/mica et enfin le
matériau central d’épaisseur « e » est le matériau composite. De ce fait, ce dernier
est placé au centre du moule où le champ électrique est maximum.
4.2.2.3. Adaptation des interfaces
Les ondes électromagnétiques existent dans différents modes suivant

l’appareil utilisé. Comme nous l’avons indiqué précédemment, nous distinguons
deux régimes de propagation d’ondes, le régime d’ondes stationnaires et le régime
d’ondes progressives. Le premier régime est caractérisé par des ventres et des
nœuds correspondant à un champ électrique maximum et nul respectivement, par
conséquent un champ non uniforme. Pour avoir un traitement uniforme au sein du
matériau, il faut éviter le régime d’ondes stationnaires, c’est-à-dire éviter les
multiples réflexions d’ondes au sein du matériau composite et du moule,
conduisant à une superposition de train d’ondes.
Cette réflexion d’ondes est d’autant plus importante que la valeur de la

permittivité du milieu est élevée. De plus la réflexion dépend de l’angle d’incidence
des ondes, ce qui se traduit par la nécessité de construire des structures
électromagnétiques spécifiques.
Outifa [I-23] a étudié ce phénomène de réflexion d’ondes aussi bien à
l’intérieur qu’à l’extérieur du moule, dans le but de réduire la puissance réfléchie, à
l’aide d’un choix convenable de la géométrie des deux transitions diélectriques
adaptées en biseaux (saillants ou rentrants) à l’entrée et à la sortie du guide
d’ondes (figure 1-10). Les résultats ont révélé, d’une part, que les biseaux rentrants
donnent de meilleurs résultats théoriques par rapport aux biseaux saillants.
D’autre part, ces mêmes résultats montrent que le biseau saillant placé à l’entrée
favorise la concentration du champ électrique au milieu de l’applicateur où se
trouve le matériau composite. Cependant des biseaux non identiques aux interfaces
auraient pour conséquence une forte distorsion du champ électrique dans
l’applicateur. Par conséquent, les biseaux sortants ont été retenus aussi bien pour
l’interface d’entrée que pour l’interface de sortie, car ils donnent de bons résultats
expérimentaux et permettent l’obtention d’un taux de réflexion inférieur à 5 %,
entraînant une transmission quasi-totale de l’onde incidente.
- 44 -
Ainsi l’adaptation diélectrique en biseaux évite non seulement l’apparition

d’ondes stationnaires au sein de l’ensemble moule et objet à traiter, mais permet
aussi de maintenir un mode appelé « onde progressive » dont la direction de
propagation est unique.
4.2.2.4. Compensation du phénomène d’atténuation
Le régime d’ondes progressives a pour effet d’améliorer l’uniformité de la

distribution spatiale du champ électrique, c’est-à-dire d’obtenir une meilleure
distribution spatiale des sources de chaleur électromagnétiques qu’en régime
d’ondes stationnaires. Seulement, l’absorption d’énergie dans la direction de
propagation produit une atténuation du champ électrique, et par conséquent
produit une décroissance des valeurs des sources de chaleur dans cette direction.
Pour obtenir une uniformité spatiale des sources de chaleur dans la direction de
propagation, il existe deux solutions :
! l’emploi de sources micro-ondes alternées

! la compensation de l’atténuation par le chargement diélectrique
4.2.2.4.1. Utilisation de sources micro-ondes alternées
La première solution utilise deux sources micro-ondes placées aux deux

extrémités de l’ensemble matériau composite et moule diélectrique. Les deux
directions de propagation de ces deux sources sont de même axe mais de directions
opposées.
Si les deux sources fonctionnaient simultanément, un mode d’onde
stationnaire apparaîtrait résultant de la sommation des deux trains d’ondes. Par
contre, si les deux émissions micro-ondes sont alternées, elles ont pour effet la
sommation des effets thermiques des deux contributions d’énergie. Le système à
deux sources micro-ondes alternées est utilisé dans notre dispositif micro-ondes.
- 45 -
4.2.2.4.2. Compensation par un chargement diélectrique
En réalisant une faible atténuation de l’onde dans le voisinage de l’interface

d’entrée, puis une atténuation progressivement croissante, il est possible de
prélever une part significative de plus en plus importante à la puissance transmise
pour obtenir une valeur absolue constante de la puissance absorbée.
Cette seconde solution utilise les caractéristiques de propagation

électromagnétiques du système moule et objet à traiter, plus précisément les
propriétés particulières du chargement diélectrique, pour compenser l’effet
d’atténuation des ondes.
En effet, l’atténuation de l’onde dans le direction de propagation est une fonction de
l’ensemble des permittivités des matériaux diélectriques constitutifs du moule.
4.3. Distribution du champ électrique à l’intérieur de

l’applicateur
Comme le suggère la formule de la puissance dissipée, la connaissance de la

distribution spatiale du champ électrique est une étape fondamentale pour une
maîtrise thermique du procédé de chauffage par hystérésis diélectrique.
Dans une approche procédé, Outifa[I-17] et al. ont développé des dispositifs micro-
ondes dans lesquels la répartition du champ électrique et de la puissance micro-
ondes sont bien connues et maîtrisées. De plus, ils ont résolu les équations de
Maxwell pour obtenir la distribution spatiale complète du champ électrique et de la
puissance micro-ondes, au sein d’un applicateur avec un chargement symétrique
comprenant deux diélectriques. Le second diélectrique est positionné au centre de
l’applicateur entre les deux couches symétriques du premier diélectrique, similaire
au chargement utilisé dans ce travail.
Bien que notre système soit symétrique, l’objectif ici est d’étudier de manière plus
générale un modèle analytique qui permette d’obtenir la distribution spatiale du
champ électrique d’une onde progressive se propageant en mode TE0,n à travers un
applicateur chargé symétriquement ou non et ce sans faire de simplification sur la

géométrie des matériaux diélectriques.
- 46 -
L’avantage d’un modèle analytique par rapport à un modèle numérique est

que le modèle numérique ne prend pas en compte les caractéristiques du mode de
propagation. Cet aspect est pourtant primordial pour une installation
expérimentale. Une résolution mathématique des équations de Maxwell a été
réalisée dans le but d’obtenir la distribution transverse du champ électrique, c’est-
à-dire suivant la direction 0y.
4.3.1. Mise en équation
En milieu diélectrique parfait, c’est-à-dire en l’absence de charges libres (et

de courant de conduction), les équations de propagation se réduisent aux équations
1-7 et 1-8.
De plus le mode de propagation est le mode fondamental TE0,1, ce qui

implique que E x , E z et Hx sont nuls.
Ainsi les composantes scalaires de ces équations se réduisent à :
∂E x ∂Hy
= −µ% (1.46)
∂z ∂t
∂E x ∂Hz
− = −µ% (1.47)
∂y ∂t
∂Hz ∂Hy ∂E x
− = ε% (1.48)
∂y ∂z ∂t
En intégrant le fait que le champ électrique est dépendant de t et z, suivant
l’expression suivante : e(
jωt −γz )
les trois précédentes équations deviennent :
−γ% E x = − jωµ% Hy (1.49)
∂E x
− = − jωµ% Hz (1.50)
∂y
∂Hz
+ γ%H y = jωε% E x (1.51)
∂y
- 47 -
En utilisant les équations (1-49) et (1-51), on obtient l’équation différentielle

vérifiée par E x :
∂ 2E x
∂y 2
+ (ω ε% µ% + γ )
2 2
Ex = 0 (1.52)
avec γ% = α + jβ
Les précédentes équations permettent également de déterminer les

composantes Hy et Hz du champ magnétique. La perméabilité magnétique de
chaque matériau est constante et égale à µ 0 , alors que la permittivité diélectrique
n’est pas la même pour chacun des matériaux. Nous nous intéressons à la
distribution du champ électrique au sein de l’applicateur, ce qui nous amène à
résoudre l’équation (1.52), pour connaître la composante du champ électrique E x .
Plus précisément, nous souhaitons obtenir le champ électrique transversal E x ( y ) .
Le facteur d’absorption ε′′ a été négligé tout au long du calcul de la

répartition du champ électrique transverse. Par conséquent, lors de la résolution
des équations, la permittivité ε% sera remplacée uniquement par la partie réelle ε′
de chacun des matériaux diélectriques de remplissage. Cette simplification peut
être opérée lorsque le facteur d’absorption des matériaux diélectriques employés
n’excède pas 0.5 pour des permittivités de l’ordre de 5, ce qui est bien notre cas
suivant le tableau suivant :
Silicone/Verre Silicone/Mica Matériau Composite

I et V II et IV III
c = 13.5 mm d = 24 mm e = 5 mm
ε′′ 0.03 0.01 0.5
ε′ 4.5 5.5 4
Tableau 1 - 3 : Caractéristiques diélectriques des matériaux diélectriques de

remplissage de l’applicateur.
- 48 -
4.3.2. Résolution des équations de propagation dans l’applicateur

chargé
Suite aux précédentes hypothèses, puis en considérant les quatre couches de

matériaux diélectriques, ainsi que le matériau composite situé au centre du moule
(figure 1-7 ), et en prenant en compte un champ électrique de la forme :
E x ( y, z,t ) = Eo ⋅ E x ( y ) ⋅ e −αz ⋅ e (
j ωt −βz )
(1.53)
nous sommes amenés à résoudre les équations différentielles suivantes, vérifiées

par E x ( y ) , qui est une fonction relative, dans les cinq régions de l’applicateur
chargé[I-16, I-17] :
d2EIx ( y )
0<y < d :
dy 2
( )
+ ε1′ k 02 − β2 EIX ( y ) = 0 (1.54)
d2EIIx ( y )
d < y < d+c :
dy 2 ( )
+ ε′2 k 02 − β2 EIIx ( y ) = 0 (1.55)
x (y)
d2EIII
d+c < y < d+c+e :
dy 2 ( )
X (y) = 0
+ ε′3 k 02 − β2 EIII (1.56)
x (y)
d2EIV
d + c + e < y < d + 2c + e :
dy 2 ( )
x (y ) = 0
+ ε′2 k 02 − β2 EIV (1.57)
d2EVx ( y )
d + 2c + e < y < a :
dy 2
( )
+ ε1′ k 20 − β2 E Vx ( y ) = 0 (1.58)
Introduisons les grandeurs p, q, r comme l’ont effectué plusieurs auteurs

auparavant[I-18, I-19, I-20, I-21, I-22], ces grandeurs sont telles que :
p ε′ 1
= ε1′ ⋅ k 02 − β2 = 2π 12 − 2 (1.59)
d λ0 λ g
- 49 -
q ε′ 1
= ε′2 ⋅ k 02 − β2 = 2π 22 − 2 (1.60)
c λ0 λ g
r ε′ 1
= ε′3 ⋅ k 02 − β2 = 2π 32 − 2 (1.61)
e2 λ0 λ g
2π 2π
car k 0 = et β = .
λ0 λg
p q 2r
, et sont appelées constantes de propagation transversales
d c e
normalisées dans les régions (I et V), (II et IV) et (III) respectivement. Suivant les
valeurs de permittivité de chaque diélectrique, ces grandeurs peuvent être réelles ou
imaginaires. L’introduction de ces constantes permet la résolution plus facile des
équations (1.54) à (1.58), qui conduit aux expressions suivantes du champ
électrique dans les cinq régions :
p 
EIx ( y ) = sin  y + θ1  (1.62)
d 
q  q 
EIIx ( y ) = A sin  y  + B cos  y  (1.63)
c  c 
2r  2r 
x ( y ) = C sin 
EIII y  + D cos  y (1.64)
 e   e 
q  q 
x ( y ) = E sin 
EIV y  + F cos  y  (1.65)
c  c 
p 
E Vx ( y ) = G sin  y + θ2  (1.66)
 d 
A, B, C, D, E, F, G ainsi que θ1 et θ2 sont des constantes d’intégration qui

seront déterminées ultérieurement grâce aux conditions aux limites suivantes :
• Champ électrique nul aux parois métalliques du guide d’onde

Ex ( y = 0) = 0 (1.67)
- 50 -
Ex (y = a ) = 0 (1.68)
• Continuité des champs électrique E x ( y ) et magnétique Hz ( y ) aux
interfaces de séparation de deux diélectriques, liée d’une part à la

direction du champ électrique parallèle à l’interface et à la constance de la
perméabilité µ r = 1 dans chaque diélectrique
EIx ( y = d ) = EIIx ( y = d ) (1.69)
EIIx ( y = d + c ) = EIII
x (y = d + c) (1.70)
x ( y = d + c + e) = Ex ( y = d + c + e)
EIII IV
(1.71)
x ( y = d + 2c + e ) = E x ( y = d + 2c + e )
EIV V
(1.72)
HIz ( y = d ) = HIIz ( y = d ) (1.73)
HIIz ( y = d + c ) = HIII
z (y = d + c) (1.74)
z ( y = d + c + e ) = Hz ( y = d + c + e )
HIII IV
(1.75)
z ( y = d + 2c + e ) = Hz ( y = d + 2c + e )
HIV V
(1.76)
Les équations (1.67) et (1.68) permettent de déterminer les valeurs des

constantes θ1 et θ2 . Ainsi on en déduit d’après ces équations que :
θ1 = 0 (1.77)
p
θ2 = − a (1.78)
d
En intégrant les équations (1.77) et (1.78) dans les équations (1.62) et (1.66) on
obtient les expressions finales de E x ( y ) :
p 
EIx ( y ) = sin  y  (1.79)
d 
q  q 
EIIx ( y ) = A sin  y  + B cos  y  (1.80)
c  c 
- 51 -
2r  2r 
x ( y ) = C sin 
EIII y  + D cos  y (1.81)
 e   e 
q  q 
x ( y ) = E sin 
EIV y  + F cos  y  (1.82)
c  c 
p 
E Vx ( y ) = G sin  ( y − a )  (1.83)
d 
Les équations (1.69) à (1.72) traduisent la continuité du champ électrique,

tandis que les équations (1.73) à (1.76) traduisent celle du champ magnétique aux
interfaces diélectriques. Sachant que le champ magnétique Hz est proportionnel à
dE x
, on obtient le système d’équations suivants :
dy
EIx ( y = d ) = EIIx ( y = d ) (1.84)
dEIx dEII
l’équation (1.73) devient : (y = d) = x (y = d) (1.85)
dy dy
EIIx ( y = d + c ) = EIII
x (y = d + c) (1.86)
dEIIx dEIII
l’équation (1.74) devient : (y = d + c) = x (y = d + c) (1.87)
dy dy
x ( y = d + c + e) = Ex ( y = d + c + e)
EIII IV
(1.88)
dEIII dEIV
l’équation (1.75) devient : x
(y = d + c + e) = x (y = d + c + e) (1.89)
dy dy
x ( y = d + 2c + e ) = E x ( y = d + 2c + e )
EIV V
(1.90)
dEIV dE V
l’équation (1.76) devient : x
( y = d + 2c + e ) = x ( y = d + 2c + e ) (1.91)
dy dy
Compte tenu des expressions du champ électrique données par les équations
(1.79) à (1.83), les équations (1.84) à (1.91) deviennent :
q ⋅d  q ⋅d 
sin (p ) = A sin   + B cos   (1.92)
 c   c 
- 52 -
p q q ⋅d   q ⋅ d 
cos ( p ) =  A cos   − B s in   (1.93)
d c  c   c 
q  q   2⋅r 2⋅r
A sin  (d + c )  + B cos  (d + c )  = C sin  (d + c )  + D cos  (d + c )  (1.94)
c  c   e   e 
q q  q  2 ⋅ r   2⋅r 2⋅r 
 A cos  (d + c )  − B sin  (d + c )   =  C cos  (d + c )  − D sin  (d + c ) 
c c  c  e   e   e 
(1.95)
2⋅r 2⋅r q q
C sin  (d + c + e )  + D cos  (d + c + e )  = E sin  (d + c + e )  + F cos  (d + c + e ) 
 e   e  c  c 
(1.96)
2⋅r   2⋅r 2⋅r  q 

(d + c + e )  − D sin  (d + c + e )   = E cos  (d + c + e )  − F sin  (d + c + e )  
q q
 C cos 
e   e   e  c  c  c 
(1.97)
q  q 
E sin  (d + 2c + e )  + F cos  (d + 2c + e )  = −G sin (p ) (1.98)
c  c 
q q  q  p
 E cos  (d + 2c + e )  − F sin  (d + 2c + e )   = G cos ( p ) (1.99)
c c  c  d
Les équations (1.92) et (1.93) forment un système de deux équations à deux

inconnues. La résolution de ce système permet d’obtenir les constantes
d’intégration A et B :
q ⋅d  p⋅c q ⋅d
A = sin ( p ) ⋅ sin  + cos (p ) ⋅ cos   (1.100)
 c  q ⋅d  c 
q ⋅d  p⋅c q ⋅d 
B = sin ( p ) ⋅ cos  − cos ( p ) ⋅ sin   (1.101)
 c  q ⋅d  c 
Connaissant A et B, les équations (1.94) et (1.95) constituent à leur tour un

système de deux équations à deux inconnues en C et D, dont la résolution donne le
résultat suivant :
- 53 -
 p⋅c   2⋅r 
C= cos ( p ) sin (q ) + sin ( p ) cos (q )  ⋅ sin  
q ⋅d   e (d + c ) 
(1.102)
e p q   2⋅r 
+  cos ( p ) cos (q ) − sin (p ) sin (q )  ⋅ cos  
2⋅r d c   e (d + c ) 
 p⋅c   2⋅r 
D= cos (p ) sin (q ) + sin (p ) cos (q )  ⋅ cos  
q ⋅d   e (d + c ) 
(1.103)
e p q   2⋅r 
−  cos ( p ) cos (q ) − sin (p ) sin (q )  ⋅ sin  
2⋅r  d c   e (d + c ) 
De même la connaissance des constantes C et D permett à présent de

résoudre le système de deux équations en deux inconnues en E et F que forment les
équations (1.96) et (1.97). En résolvant ce système on obtient :
 p ⋅ c 
E =  cos ( p ) cos (q ) − sin ( p ) sin (q )  ⋅ cos (2r )
 q ⋅ d 
2⋅c ⋅r  p ⋅c    q ⋅ (d + c + e ) 
−  sin ( p ) cos (q ) + cos (p ) sin (q )  ⋅ sin (2r ) ⋅ cos  
q ⋅e  q ⋅d    c 
(1.104)
 p ⋅ c 
+  cos ( p ) sin (q ) + sin ( p ) cos (q )  ⋅ cos (2r )
 q ⋅ d 
e p q    q ⋅ (d + c + e ) 
+  cos ( p ) cos (q ) − sin ( p ) sin (q )  ⋅ sin (2r ) ⋅ sin  
2⋅r  d c    c 
 p ⋅ c 
F =  − cos ( p ) cos (q ) + sin ( p ) sin (q )  ⋅ cos (2r )
 q ⋅ d 
2⋅c ⋅r  p⋅c    q ⋅ (d + c + e ) 
+  sin (p ) cos (q ) + cos (p ) sin (q )  ⋅ sin (2r ) ⋅ sin  
q ⋅e  q ⋅d    c 
(1.105)
 p ⋅ c 
+  cos ( p ) sin (q ) + sin (p ) cos (q )  ⋅ cos (2r )
 q ⋅ d 
e p q    q ⋅ (d + c + e ) 
+  cos ( p ) cos (q ) − sin ( p ) sin (q )  ⋅ sin (2r ) ⋅ cos  
2⋅r  d c    c 
- 54 -
L’équation (1.99) permet de déterminer la constante G puisque :
  q ⋅ (d + 2c + e )   q ⋅ (d + 2c + e )  
E cos   − F sin  
q ⋅ d   c   c  
G= ⋅ (1.106)
p⋅c cos ( p )
De même, l’équation (1.98) fournit une autre expression de la constante G :
  q ⋅ (d + 2c + e )   q ⋅ (d + 2c + e ) 
E sin   + F cos  

  c   c 
G= (1.107)
− sin (p )
En égalant les deux expressions de la constante d’intégration G obtenues par

les équations (1.106) et (1.107), on obtient l’équation suivante :
p   q ⋅ (d + 2c + e)   q ⋅ (d + 2c + e)   q   q ⋅ (d + 2c + e)   q ⋅ (d + 2c + e)  
cotan(p) ⋅ Esin  + F cos    = Fsin  − E cos   
d  c c  
     c   c   c 
(1.108)
Cette équation dépendant uniquement de λ g (puisque les constantes E et F
sont fonction de p, q et r, qui dépendent eux-mêmes de λ g ) est une équation
transcendante. La connaissance des paramètres ε1′ , ε′2 , ε′3 permet de déterminer la
valeur unique de la longueur d’onde du mode de propagation TE0,1 dans
l’applicateur λ g (inconnue de l’équation transcendante).
4.3.3. Simulation du champ électrique sur MapleV
Toutes les constantes d’intégration ayant été déterminées, on en déduit

l’expression complète du champ électrique dans chaque matériau diélectrique :
0 < y < d : E x ( y, z,t ) = E0 ⋅ EIx ( y ) ⋅ e ( ) ⋅ e jωt

− α+ jβ ⋅z
(1.109)
d < y < d+c : E x ( y, z,t ) = E0 ⋅ EIIx ( y ) ⋅ e −(α+ jβ)⋅z ⋅ e jωt (1.110)
- 55 -
d+c < y < d+c+e : E x ( y, z,t ) = E 0 ⋅ EIII

x (y ) ⋅ e
− (α+ jβ )⋅ z
⋅ e jωt (1.111)
d+c+e <y < d+2c+e : E x ( y, z,t ) = E0 ⋅ EIV

x (y ) ⋅ e
− ( α+ jβ )⋅z
⋅ e jωt (1.112)
d+2c+e < y < a : E x ( y, z,t ) = E0 ⋅ E Vx ( y ) ⋅ e−(α+ jβ )⋅z ⋅ e jωt (1.113)
λ g P0 λ g µ 1
avec E0 = ⋅ ⋅ 0 ⋅
p1 b λ 0 ε0 p1 + p2 + p3
d 2
où p1 = ∫ EIx ( y ) dy
0
d +c 2
p2 = ∫ EIIx ( y ) dy
d
a 2 2
p3 = ∫ x ( y ) dy
EIII
d +c
Les expressions de EIx ( y ) , EIIx ( y ) , EIII

x ( y ) sont données par les équations
(1.79) à (1.81).
Sachant que E0 dépend de λ g , pour connaître la distribution du champ
électrique au sein de l’applicateur, il nous faut résoudre préalablement l’équation

transcendante (1-108) en λ g .
4.3.3.1. Résolution de l’équation transcendante
La résolution de la totalité des équations, que ce soit pour déterminer les

constantes d’intégration ou pour résoudre l’équation transcendante, a été effectuée
par un calcul formel sur le logiciel mathématique MapleV.
L’équation transcendante a été résolue pour les conditions expérimentales

résumées dans le tableau 1-3 précédent. Ces conditions fournissent une solution
unique de l’équation transcendante, correspondant à la longueur d’onde du mode
fondamental TE0,1 . La détermination numérique de cette longueur d’onde de
- 56 -
propagation s’effectue par le choix d’intervalles de plus en plus étroits autour de la

valeur recherchée.
Pour des conditions expérimentales différentes, il est possible d’obtenir

plusieurs solutions de l’équation transcendante. Ces solutions correspondent à la
longueur d’onde λ g pour tous les modes TE0,n qui peuvent exister pour ces
conditions. Dans ce cas, un tri parmi les valeurs de λ g solutions de l’équation
transcendante doit être effectué. En effet, seules les solutions physiquement

acceptables sont retenues et les solutions ne possédant pas de sens physique sont
directement éliminées, comme des longueurs d’ondes négatives par exemple. Ainsi,
la plus forte valeur de λ g correspond au mode TE0,1 tandis que la plus faible valeur
de λ g correspond au mode TE 0,n d’ordre le plus élevé qui puisse exister.
Une fois λ g connue, on a accès aux grandeurs de propagation transversales
p, q et r, puis on peut atteindre l’expression du champ électrique, comme cela est

décrit sur l’organigramme figure 1-12 :
- 57 -
Equations de Maxwell
Equations de Propagation
Application des conditions de

Résolution des équations uv uv
uv uv continuité de E et H aux
différentielles en E et H
interfaces des diélectriques
Ecriture des solutions

mathématiques
Equation transcendante
en λ g
Résolution analytique sur

Maple
Composantes du champs
uv
électrique E
Figure 1 - 12 : Organigramme de la méthode de calcul d’un applicateur chargé par des

diélectriques.
- 58 -
4.3.3.2. Distribution du champ électrique pour un applicateur

chargé par cinq matériaux diélectriques
La répartition du champ électrique dans l’applicateur chargé par cinq

matériaux diélectriques est représentée figure 1-13.
12000
10000
Champ électrique (V/m)
8000
6000
4000
2000
0
0.00 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08
y (m)
Figure 1 - 13 : Distribution du champ électrique au sein de l’applicateur.
On observe que les deux maximums du champ électrique sont bien situés
dans les matériaux diélectriques de permittivité la plus élevée. Ceci est à corréler
avec la réduction la plus intense de la vitesse de propagation dans les matériaux à
permittivité la plus élevée. Ce résultat est bien en accord avec celui obtenu par
Outifa[I-24] pour un chargement identique de l’applicateur mais par une méthode de
résolution des équations uniquement adaptée aux chargements symétriques
d’applicateur.
- 59 -
4.4. Interaction de l’onde électromagnétique sur la matière
L’échauffement d’un produit grâce au rayonnement micro-ondes est

provoqué par la dissipation sous forme de chaleur d’une partie de l’énergie
contenue dans cette onde électromagnétique.
4.4.1. Puissance électromagnétique dissipée en chaleur
Les matériaux que l’on traite à l’aide des micro-ondes sont des diélectriques
constitués de dipôles électriques. Le champ électrique joue le rôle d’une contrainte
électrique qui tend à organiser ces dipôles au sein de la matière ; ce phénomène est
la polarisation diélectrique de la matière. Une fois cette contrainte interrompue, les
dipôles retournent vers leur état d’équilibre avec une certaine constante de temps ;
ce phénomène est la relaxation diélectrique. Enfin, le retard de la réponse du
système, considéré comme un ensemble de charges électriques, soumis au champ
électrique variable au cours du temps, est appelé hystérésis diélectrique ; ce
troisième phénomène est la base de la création des sources internes de chaleur.
La caractéristique diélectrique d’absorption ε′′ évalue la densité volumique de

puissance électrique convertie en chaleur au sein du matériau et permet de relier
uv
cette densité Pdiss à l’amplitude du champ électrique E au sein du matériau.
La puissance dissipée par unité de volume est alors donnée par la formule
suivante :
ω 2 2
Pdiss = ⋅ ε′′ ⋅ E = π ⋅ f0 ⋅ ε0 ⋅ ε′′r ⋅ E
2
Cette formule montre que la puissance dissipée, qui s’exprime en W/m3, est
fonction du produit de la fréquence de travail, de la caractéristique d’absorption
diélectrique, de la permittivité du vide et du carré du module du champ électrique.
La caractéristique d’absorption n’est donc pas le seul paramètre à gouverner
la dissipation d’énergie. Ceci explique, par exemple, que l’on puisse décongeler par
micro-ondes l’eau à l’état de glace, bien que le facteur d’absorption à la fréquence
micro-ondes soit beaucoup plus faible qu’à basse fréquence. Cet exemple souligne
- 60 -
l’importance du calcul du champ électrique E x ( y ) , que nous effectuerons
ultérieurement.
Ainsi, lorsque les transferts de chaleur entre le cœur de l’objet à chauffer et

son environnement sont négligeables, la vitesse d’élévation locale de la température
dT
est reliée à la densité volumique de puissance électrique convertie en chaleur
dt
par l’expression :
dT P
= diss
dt ρ ⋅ Cp
avec ρ la masse volumique et Cp la chaleur spécifique massique à pression

constante.
4.4.2. Phénomène d’atténuation de l’onde
La caractéristique d’absorption ε′′ joue aussi un très grand rôle dans la

distribution spatiale des sources internes de chaleur à l’intérieur du matériau. Du
fait de la conversion d’une partie de la puissance véhiculée par les ondes
électromagnétiques en quantité de chaleur créée par unité de temps, l’amplitude
uv uv
des champs E et H est décroissante en fonction de la coordonnée d’espace suivant
la direction de propagation, Oz, dans notre cas.
uv uv
Par exemple, si le champ électrique est de la forme E ( z ) = E0 ⋅ exp ( −αz ) , alors
la puissance véhiculée par l’onde s’écrit P ( z ) = Pi ⋅ exp ( −2αz ) , donc la densité de
puissance dissipée prendra la forme Pdiss = p0 ⋅ exp ( −2αz ) .
Le coefficient d’atténuation α associé à cette atténuation est une fonction de la

caractéristique d’absorption ε′′ et s’écrit :
α=
π
⋅
1
⋅
(ε1′′ ⋅ p1 + ε′′2 ⋅ (p2 − p3 ) + ε′′3 ⋅ p3 )
λ0 λ0 λg p1 + p2
- 61 -
où ε1′′ , ε′′2 , et ε′′3 sont respectivement les caractéristiques diélectriques d’absorption

du diélectrique silicone/verre, du diélectrique silicone/mica et du matériau
composite.
p1 , p 2 , p3 sont identiques à ceux définis au paragraphe précédent.
Les détails concernant ce calcul sont donnés dans la thèse de Outifa[I-23].
Rappelons l’expression de α dans le cas d’une propagation libre[I-24] :
ω ε′′ π ε′′
α= ⋅ = ⋅
2c ε′ λ 0 ε′
5. ASPECT CINETIQUE
Les matériaux utilisés dans ce travail sont des matériaux composites à

matrice thermodurcissable du type résine époxyde associée à des fibres de verre.
Par conséquent, il est nécessaire de connaître l’évolution de la loi cinétique au cours
de la transformation chimique de ces résines.
5.1. Transformation chimique
Lors de l’élaboration des matériaux composites à thermodurcissables résine

époxyde/fibres de verre, on assiste à une réaction de réticulation exothermique [I-25]
dans le moule. Pendant le chauffage du matériau, deux réactions d’addition se
déroulent lorsque la stœchiométrie entre l’époxyde et l’amine est respectée. La
première réaction se produit entre un hydrogène de l’amine et un groupe époxyde,
puis la seconde réaction entre l’autre hydrogène de l’amine et un autre groupe
époxyde. A l’issu de cette réaction, on obtient la formation d’un réseau
tridimensionnel à base d’époxyde.
- 62 -
Figure 1 - 14 : Les réactions d’addition entre une fonction époxyde et une diamine.
Le prépolymère utilisé dans cette étude est le diglycidyléther du bisphénol A,

connu sous le nom de DGEBA, dont la masse molaire est de 340 g/mol pour n = 0.
Figure 1 - 15 : formule chimique de la résine DGEBA.
Le durcisseur employé est le dicyandiamide, connu sous la dénomination DDA.

Ce durcisseur est solide à la température ambiante et sa température de fusion est
de 207°C.
Figure 1 - 16 : Formule chimique du DDA.
Le catalyseur présent dans le mélange est le benzyldiméthylamine, appelé

également BDMA. Ce catalyseur est liquide à température ambiante et sa
température d’ébullition est de 180°C. Il s’agit d’un catalyseur plutôt classique qui
est souvent associé au durcisseur DDA.
- 63 -
5.2. Modélisation cinétique
Les réactions de réticulation des mélanges à base de résine époxyde et

d’amine ont fait l’objet de nombreux travaux conduisant au développement de
plusieurs modèles cinétiques[I-26, I-27, I-28, I-29, I-30, I-31, I-32, I-33, I-34, I-35, I-36].
Dans les travaux rapportés, le terme source chimique est égal à :
dH dx
= ρ ⋅ ∆HT ⋅ (1.114)
dt dt
dont la valeur d’enthalpie totale de la réaction est d’environ ∆HT = 450 J/g, et la
masse volumique ρ = 1 200 kg/m3 pour le système qui nous concerne.
Sachant que le terme source chimique est proportionnel à la vitesse de

dx dx
réaction , alors la connaissance de la loi cinétique (t) au cours de la
dt dt
transformation du matériau est nécessaire pour une bonne maîtrise de l’aspect
thermique du procédé.
5.2.1. Modèle cinétique
Il existe plusieurs modèles cinétiques concernant les réactions de réticulation

de formulations à base de résines thermodurcissables et plus particulièrement de
résines époxydes/amine. Le modèle le plus souvent employé pour les mélanges
époxyde/amine est le modèle cinétique proposé par Kamal et Sourour [I-37], qui
prend en compte l’ensemble des paramètres de la réaction même lorsque
l’hypothèse d’équiréactivité n’est pas vérifiée. L’équation différentielle de Kamal et
Sourour est rappelée ci-dessous :
dx
= (K 2 + K1 ⋅ x m )(1 − x)n (1.115)
dt
- 64 -
où les constantes K1 et K 2 sont les constantes de vitesse de réaction et dépendent

de la température selon des lois d’Arrhenius faisant intervenir les énergies
d’activation E1 et E2 :
Ei
K i = k i exp(− ) avec i = [1,2] (1.116)
R⋅T
avec k 1 = 1.6 ⋅ 106 min−1 , E1 = 68.3 kJ/mol , k 2 = 255 min−1 et E2 = 38.7 kJ/mol .
5.2.2. Hypothèses sur le modèle cinétique et les cas traités :
Le calcul formel sur le logiciel MAPLEV s’est décomposé en deux étapes.

Dans un premier temps, le modèle cinétique a été simplifié pour déterminer
l’influence des paramètres tout en réduisant les temps de calcul. Ensuite, les
résultats obtenus par le modèle simplifié doivent satisfaire les conditions issues de
résultats expérimentaux permettant d’obtenir une bonne interprétation des effets de
la cinétique.
Le modèle proposé par Kamal et Sourour est un modèle général de cinétique

de réaction qui prend en compte l’hypothèse de non équiréactivité à l’aide des
exposants m et n. Etant donné que l’objectif de ce travail n’est pas d’étudier les
mécanismes réactionnels du système DGEBA/DDA, mais uniquement de connaître
l’évolution de la conversion au cours de la réticulation du matériau, nous
supposerons l’hypothèse d’équiréactivité des réactions satisfaite pour simplifier le
modèle. Par conséquent, en prenant les coefficients m et n égaux à 1, le modèle
cinétique devient :
dx
= (K 2 + K1 ⋅ x)(1 − x) (1.117)
dt
Pour résoudre cette équation quatre conditions dénommées C1 à C4 ont été

étudiées :
• C1 : x = x(t, Τ) x est une sigmoïde en temps et en

température
- 65 -
• C2 : x (t, Τ = 120°C) = 0 le début de la transformation a lieu vers

120°C
• C3 : x (t, Τ = 200°C) = 1 la transformation est supposée complète
vers 200°C
dx
• C4 : (t, Τ = 150°C) la vitesse est maximale pour une
dΤ
température de 150°C
Ces conditions ont été proposées pour décrire les principales observations
expérimentales. Les différents paramètres du modèle ont été déterminés de manière
à satisfaire les conditions précitées.
Différents cas traités :
La simulation s’est faite par étapes de difficulté croissante, en considérant

une progressivité des modèles à partir de modèles simples tel que le cas 1 cité ci-
dessous, pour aboutir à un modèle plus élaboré, s’approchant davantage de notre
modèle, qui est le cas .
dx
Cas 1 : = K1 ⋅ x(t)
dt
dx
Cas 2 : = K1 ⋅ x(t) + K 2
dt
dx
Cas 3 : = K1 ⋅ (1 − x(t))
dt
dx
Cas 4 : = K1 ⋅ x(t) ⋅ (1 − x(t))
dt
dx
Cas 5 : = (K1 ⋅ x(t) + K 2 ) ⋅ (1 − x(t))
dt
La résolution des équations des cas 1 et 2 donne des expressions de la

conversion en fonction du temps du type :
x(t, T) = A ⋅ exp (B ) (1.118)
avec A, une constante et B = f(t, T,K i ) .
- 66 -
Or, ce type d’équation admet une pente nulle à l’origine, en conséquence elle
satisfait bien la condition C2, puis elle est quasi-infinie à partir d’une certaine
valeur du fait de la prépondérance de la partie exponentielle de la fonction.
En conséquence, les cas 1 et 2 ne conviennent pas à la condition C3, imposée à la
conversion.
La solution de l’équation du cas 3 est de la forme :
x(t, T) = 1 − C ⋅ exp (D)
avec C, une constante et D = f(t, T,K i ) .
On remarque que contrairement aux cas 1 et 2, la solution du cas 3 est une

fonction convergente. Bien que la conversion converge vers la valeur 1 à l’infini
(condition C3 vérifiée), la pente à l’origine est non nulle. Le cas 3 ne satisfait donc
pas à la condition C2 imposée sur la conversion.
On peut noter que le cas 4 est une réunion des cas 1 et 3, et la conversion
solution de l’équation du quatrième cas est du type :
1
x(t, T) = (1.119)
(1 + E ⋅ exp(F))
avec E, une constante et F = f(t, T,K i ) .
Cette expression admet une pente nulle à l’origine (ce qui signifie que la
condition C2 est satisfaite), puis converge vers la valeur 1 à l’infini (la condition C3
est également vérifiée) et enfin elle admet un maximum sur la dérivée première. Les
valeurs des constantes de Ki ont été ajustées par rapport à ce maximum pour
correspondre aux conditions C4, d’où les valeurs suivantes :
k1 = 3.6 ⋅ 105 min−1

E1 = 68.3 kJ/mol
- 67 -
Ce résultat présente par contre un défaut, qui est une fin de réaction atteinte
très rapidement d’après l’équation de conversion.
Le cas 5 est une réunion des cas 2 et 3. Toutefois le cas 4 traité

précédemment suffit à décrire le phénomène de transformation chimique de notre
formulation à base de résine DGEBA et de durcisseur amine.
En conclusion, le cas 4 est conservé pour être intégré par la suite, au cours
du chapitre 3, dans les programmes ANSYS afin de simuler les effets de la
réticulation pendant l’étape de chauffage par hystérésis diélectrique. La figure ci-
dessous présente les résultats obtenus sur MAPLE à partir du modèle cinétique
simplifié de Kamal et Sourour pour ce cas.
1.0
0.8
Conversion chimique
0.6
0.4
0.2
0.0
0 20 40 60 80 100 120 140 160 180 200 220 240 260
Température (°C)
Figure 1 - 17 : Evolution en fonction de la température de la conversion chimique

issue de la modélisation.
6. DESCRIPTION DE L’ECOULEMENT
L’emploi de sources internes micro-ondes permet la réduction du temps

d’élaboration de pièces en composite thermodurcissable. La réduction du temps de
- 68 -
mise en œuvre et la forte variation de viscosité des matériaux thermodurcissables

nécessitent d’intégrer dans la gestion du procédé de fabrication de ces matériaux
des conditions permettant la maîtrise du fluage de la matrice polymère durant la
phase d’élaboration, en optimisant la distribution de la contrainte de pression
appliquée à l’objet à traiter. Les bases de cette optimisation seront évoquées en
dernière partie du chapitre III.
L’objectif est de savoir si la réduction du temps d’élaboration permet un transfert de
matière satisfaisant. Dans cette partie, nous nous intéresserons essentiellement à
l’aspect théorique de l’écoulement de la matière et à la distribution de la pression
au sein du matériau.
Dans un premier temps, nous présenterons les expériences de mesure du fluage de

la matrice organique réalisée par une mise en température de manière isotherme de
matériaux composites, suivie d’une compression de ces derniers à l’aide d’une
presse hydraulique.
Puis nous présenterons les calculs théoriques de trois cas de figure d’écoulement
(figure 1-18 et 1.20) dans un moule de forme longiligne ouvert à ses deux petites
extrémités sous un dispositif permettant d’appliquer une contrainte de pression :
! écoulement d’un matériau visqueux sur une longueur indéfinie,

! écoulement du même matériau visqueux, avec évacuation de la matière
Li
au delà des cotes z = ±
2
! écoulement d’un matériau composite constitué d’une matrice visqueuse et
de fibres longues, de longueur L i , avec évacuation de la matière au delà
Li
des cotes z = ±
2
le dernier cas étant le plus proche de la réalité.
- 69 -
p (z ) patm − patm
e (t )
l
L z=0 L
− +
2 2
Figure 1 - 18 : Configuration du matériau visqueux.
Ces calculs ont pour but de décrire avec des hypothèses simples des phénomènes
physiques reliés à l’écoulement de matériaux fluides composés de fibres ou non.
6.1. Observation expérimentale du fluage de la matrice

organique
Nous avons reconstitué l’écoulement de la matrice du matériau composite au

cours du procédé, par une mise en température uniforme du composite, en étuve,
suivie d’une compression de ce dernier à l’aide d’une presse hydraulique pour
extraire le mélange époxy/amine non encore réticulé. Cette mise en température
uniforme a été réalisée pour différentes températures allant de la température
ambiante jusqu’à la température de gélification.
6.1.1. Protocole expérimental
Le matériau composite obtenu par empilement de feuilles de préimprégnés

est chauffé en étuve pendant une heure et demie à différentes températures
isothermes allant de 60°C jusqu’à 120°C pour assurer une très bonne uniformité
- 70 -
thermique de l’éprouvette de dimension 20 mm * 86 mm * 300 mm. Il est ensuite

soumis à une pression de 10 bars à l’aide d’une presse hydraulique.
L’expérience consistant ainsi à extraire la matrice du matériau composite aux

différentes températures sera appelée opération « d’essorage ». De même que la
matrice obtenue par cette expérience sera nommée par la suite « résine extraite par
l’opération d’essorage ».
6.1.2. Résultat expérimental
L’opération d’essorage nous permet non seulement d’évaluer la quantité de

résine s’étant écoulée durant la mise sous contrainte de pression du matériau
composite, afin d’identifier les propriétés rhéologiques, mais aussi d’obtenir des
extraits de matrice de compositions éventuellement différentes sur lesquels les
autres propriétés physico-chimiques de la phase organique seront étudiées (voir
chapitre II). La quantité de résine ayant flué est représentée ci-dessous.
7
Fraction résine extraite (%)
0
40 50 60 70 80 90 100 110 120 130
Tessorage (°C)
Figure 1 - 19 : Variation du pourcentage de la fraction de résine extraite en fonction

de la température.
- 71 -
On constate un début de fluage de la résine à partir de 60°C et l’on arrive à extraire

jusqu’à 7 % (le maximum) en masse de résine à 110°C sur une durée totale de
300 s, durée fournissant une quantité de matière satisfaisante.
6.2. Modèles théoriques et écoulement de la matière
Pour estimer le comportement de l’écoulement de résine dans un matériau

composite, la distribution de la pression au sein du matériau doit être connue.
Pour comprendre le mécanisme d’écoulement des monomères (de masse
molaire de 340 g/mol pour la résine DGEBA n = 0) du matériau composite lors de
la phase de mise en forme, nous avons effectué des calculs formels représentant
différentes situations d’écoulement d’un matériau soumis à une contrainte de
pression de direction perpendiculaire à l’écoulement. La démarche que nous avons
employée est progressive, et repose sur une analyse de trois approches différentes
du phénomène d’écoulement dont la troisième permettra une confrontation au cas
qui nous concerne.
La première approche, même simplifiée, consiste à suivre l’écoulement d’un
matériau visqueux, de longueur variable au cours du temps, ce qui sous-entend
que la matière soumise à déplacement s’accumule à l’extrémité de l’objet.
La seconde approche reprend le cas précédent avec cette fois la possibilité
d’une évacuation de la matière aux deux extrémités de l’objet de longueur
constante.
Enfin, dans une troisième approche, nous nous sommes intéressés à un
matériau de longueur constante, constitué de fibres longues unidirectionnelles à
travers lesquelles s’écoule la matière fluide. De même que pour le cas précédent, on
suppose qu’il peut y avoir évacuation de la matière aux deux extrémités de l’objet.
On peut remarquer que cette dernière approche est effectivement la plus voisine du
problème d’écoulement qui nous concerne. En effet, il faut rappeler que nous
sommes dans une configuration de moule ouvert avec la possibilité d’écoulement
aux deux extrémités du matériau composite, constitué de longues fibres de verre
unidirectionnelles. Les différents modèles représentant chaque approche ont été
développés sous condition thermique isotherme, ce qui n’est pas le cas au cours du
procédé d’élaboration, bien évidemment. Il sera donc intéressant, par la suite, de
voir les conséquences d’une telle condition sur le fluage de la matrice organique.
- 72 -
Les objectifs dans cette partie visent à :
! connaître l’évolution du front d’écoulement de la matière L (t ) pour la
première approche
et
dz
! connaître la distribution de vitesse en fonction du temps et de z
dt
! connaître la distribution de pression p ( z ) , en fonction du temps
! déterminer la caractéristique d’écoulement de la matrice
pour les trois approches.
6.2.1. Résolution du problème
Soit un volume de matériau fluide de longueur L, de largeur l constante et

d’épaisseur e ( t ) que l’on suppose dépendante du temps, du fait de l’application de
uv
la force F dans la direction parallèle à la direction de cette épaisseur (figure 1-18).
6.2.1.1. Hypothèses des modèles
Les conditions et hypothèses posées dans les différents modèles sont les suivantes :
h-1) la longueur initiale de la pièce est L i
h-2) l’épaisseur initiale de la pièce est ei

h-3) le volume total du matériau est supposé constant et égal à :
VT = L i ⋅ l ⋅ e i = cons tan te indépendante du temps
h-4) la force totale F appliquée sur la partie supérieure du matériau est

elle-même constante et indépendante du temps
h-5) le matériau reste horizontal une fois cette force appliquée
h-6) la force de pesanteur est négligée
h-7) l’écoulement de la matière est unidirectionnel suivant la direction z
h-8) la résistance à l’écoulement due à l’air est négligeable
h-9) l’écoulement n’est pas modifié par la proximité des parois
- 73 -
h-10) la matrice est incompressible. Toutefois, elle contient des bulles d’air
provenant de l’opération de l’empilement des feuilles de préimprégnés
dans le moule. Ainsi, lors de l’application de la contrainte de pression de
10 bars au cours du procédé, les bulles peuvent rester présentes dans le
matériau mais leur volume sera divisé par 10. Cela nous amène à
supposer que les bulles d’air se déplacent à la même vitesse que la résine
6.2.1.2. Mise en équation
Les équations utilisées dans chaque approche sont :
! l’équilibre des forces

! l’équation de conservation de la matière
! l’équation liant la vitesse de déplacement de la matrice au gradient de
pression hydrostatique en z est une équation phénoménologique du type
de l’équation de Darcy[I-38, I-39, I-40) :
dz dp
= −K ( T ) (1.120)
dt dz
1
avec K ( T ) ∝ où η est la viscosité, et p est la pression de la résine.
η (T )
Quelle que soit l’écoulement considéré, l’équilibre des forces s’écrit de la

façon suivante :
+L 2 +L 2
F = l⋅∫ p ( z ) ⋅ dz = 2 ⋅ l ⋅ ∫ p ( z ) ⋅ dz (1.121)
−L 2 0
L
et en z = ± on a la condition suivante : p ( z ) = patm (C1)
2
où patm est la pression atmosphérique
par symétrie du montage, avec L = L ( t ) dans la première approche et L = L i pour
les deux autres.
- 74 -
6.2.2. 1ère approche : écoulement de matériau fluide de longueur

variable au cours du temps sans évacuation de matière à
ses extrémités
Dans cette approche, on suppose que la matière s’écoule entraînant une

variation de la longueur du matériau au cours du temps. Ainsi, le front
d’écoulement est représenté par L = L ( t ) .
L’équation de conservation de la matière, qui traduit un volume total

constant de la matière visqueuse, s’écrit alors :
ei ⋅ L i ⋅ l = cons tan te = e ( t ) ⋅ L ( t ).l (1.122)
soit e ( t ) ⋅ L ( t ) = cons tan te
d 1 d 1 d
ainsi e (t ) ⋅ L (t ) = 0 ou ⋅ L (t ) = − ⋅ e (t ) (1.123)
dt L ( t ) dt e (t ) dt
Entre les cotes – z et + z, la conservation du volume s’écrit :
2 ⋅ e ( t ) ⋅ z ⋅ l = cons tan te
d e ( t ) ⋅ z 
Ainsi on a =0 (1.124)
dt
dz
Si la tranche de fluide de cote z se déplace à la vitesse de , la relation (1.23)
dt
devient :
1 dz 1 de (t )
⋅ =− ⋅ (1.125)
z dt e ( t ) dt
- 75 -
D’après l’hypothèse disant que le matériau reste horizontal lors de la compression,

on en déduit que le second membre de l’équation (1.124) est une quantité
indépendante de la variable z. Donc en z = L(t) on trouve la relation suivante :
1 dz 1 dL (t )
⋅ = ⋅ (1.126)
z dt L (t ) dt
6.2.2.1. Détermination de L(t)
L’équation (1.120) donnant la relation entre la vitesse de déplacement de la

matière en z et le gradient de pression hydrostatique en z s’écrit aussi :
dp 1 dz
=− ⋅ (1.127)
dz K ( T ) dt
En substituant l’équation (1.126) dans l’équation (1.127) on obtient :
dp 1 1 de ( t )
= ⋅ ⋅ ⋅z (1.128)
dz K ( t ) e (t ) dt
et en intégrant cette dernière équation :
L (t ) 2 1 1 de ( t ) L (t ) 2
∫z
dp =
K (t ) e (t ) dt ∫z
⋅ ⋅ ⋅ z ⋅ dz (1.129)
1 de ( t )
car nous avons vu que la quantité ⋅ est indépendante de z.
e (t ) dt
En prenant en compte la condition (C1), l’équation (1.129) devient :
1 de ( t )  L (t ) 
2
1
p ( z ) = patm − ⋅ ⋅ ⋅ − z2  (1.130)
2 ⋅ K ( T ) e (t ) dt  4 
L (t )  L2 (t ) 
z étant inférieur à , alors  − z 2  > 0 . Le matériau subissant une
2  4 
compression, l’épaisseur est une fonction décroissante du temps :
- 76 -
 1 de ( t )   1 1 de ( t ) 
 ⋅  < 0 , et par conséquent  − ⋅ ⋅ >0
 e (t ) dt   2 ⋅ K ( t ) e ( t ) dt 
 1 de ( t )  L ( t )  
2
1
d’où : − ⋅ ⋅ ⋅ − z2  > 0 .
 2 ⋅ K ( T ) e (t ) dt  4  
Compte tenu de l’équation (1.130), et de l’équation d’équilibre des forces (1.121) :
 1 de ( t )  L ( t ) 
2
F L (t ) 2 L (t ) 2 1
= ∫ p ( z ) ⋅ dz = ∫  patm − ⋅ ⋅  − z 2   ⋅ dz
2⋅l 0 0 
 2 ⋅ K ( T ) e ( t ) dt  4 

et en intégrant de 0 à L(t)/2, on obtient le résultat suivant :
1 de ( t ) L ( t )
3
F 2 1
= patm ⋅ L (t ) − ⋅ ⋅ ⋅ ⋅ (1.131).
l 3 K ( T ) e ( t ) dt 8
En substituant l’équation (1.123) dans l’équation (1.131) on obtient l’équation

différentielle en L(t) suivante :
1 dL (t ) 2 F
⋅ ⋅ L ( t ) + patm ⋅ L (t ) =
12 ⋅ K ( T ) dt l
Cette équation différentielle s’écrit également :
1 dL (t ) F
⋅ ⋅ L (t ) = − patm
12 ⋅ K ( T ) dt l ⋅ L (t )
F
On remarque que le terme est la pression équivalente appliquée par la force
l ⋅ L (t )
F sur l’ensemble de la surface de contact avec le matériau fluide l ⋅ L (t ) .
En supposant la pression atmosphérique négligeable devant la pression moyenne,

F
c’est-à-dire si patm << , on obtient l’équation différentielle simplifiée suivante :
l ⋅ L (t )
- 77 -
1 dL (t ) 2 F
⋅ ⋅ L (t ) = (1.132).
12 ⋅ K ( T ) dt l
12 ⋅ K ( T ) ⋅ F
En écrivant cette dernière équation sous la forme : dL ( t ) ⋅ L2 ( t ) = ⋅ dt
l
et en intégrant cette nouvelle expression de L i à L (t ) l’équation s’écrit :
L (t ) L (t ) 12 ⋅ K ( T ) ⋅ F
∫ dL (t ) ⋅ L2 (t ) = ∫ ⋅ dt
Li Li l
et le résultat est tel que :
36 ⋅ K ( T ) ⋅ F
L (t ) = 3 L3i + ⋅t (1.133)
2⋅l
18 ⋅ K ( T ) F
ou même L (t ) = L i ⋅ 3 1 + ⋅ ⋅t .
L2i l ⋅ Li
dL ( t )
Connaissant l’expression de L(t), à défaut de calculer , nous pouvons à
dt
partir de l’équation (1.132) en déduire :
1 dL ( t ) 12 ⋅ K ( T ) ⋅ F 1
⋅ = ⋅ 3
L ( t ) dt l L (t )
et en remplaçant L(t) dans le membre de gauche par son expression définie par
l’équation (1.133) :
1 dL ( t ) 1
⋅ = (1.134)
L ( t ) dt l⋅L 3
3
+ ⋅t
i
12 ⋅ K ( T ) ⋅ F 2
- 78 -
6.2.2.2. Détermination de p(z)
Pour connaître la distribution de la pression dans un tel cas, revenons à

l’équation (1.126) issue de la conservation de la matière. En substituant l’équation
(1.134) dans cette équation :
1 dz 1
⋅ = (1.135)
z dt l⋅L 3
3
i
+ ⋅t
12 ⋅ K ( T ) ⋅ F 2
cette équation révèle qu’à chaque instant t donné, une relation de proportionnalité
dz
(fonction du temps) existe entre z et .
dt
Le gradient de pression étant lié à la vitesse en z par l’équation (1.120) alors on a :
dp 1 dz z 2⋅ z
=− ⋅ =− =− (1.136)
dz K ( T ) dt  l ⋅ Li
3
3   l ⋅ L3i 
K ( T ) ⋅  + ⋅ t  3 ⋅ K ( T ) ⋅  + t 
 12 ⋅ K ( T ) ⋅ F 2   18 ⋅ K ( T ) ⋅ F 
et en intégrant l’équation (1.136) comme suit :
patm 1 L (t )
∫ ( ) dp = −
p z  l ⋅ L3i 
⋅∫
z
z
3 ⋅ K ( T ) ⋅  + t 
 18 ⋅ K ( T ) ⋅ F 
1
patm − p ( z ) = − ⋅ L2 ( t ) − z 2  (1.137)
 l⋅L 3

3 ⋅ K ( T ) ⋅  i
+ t 
 18 ⋅ K ( T ) ⋅ F 
on obtient à l’aide de l’équation (1.133) :
 18 ⋅ K ( T ) ⋅ F  l ⋅ L3i 
23

1  
p ( z ) = patm + ⋅  ⋅ + t   − z 2  (1.138)
 
 l ⋅ Li   l  18 ⋅ K ( T ) ⋅ F
3
3 ⋅ K ( T ) ⋅  + t      
 18 ⋅ K ( )
T ⋅ F 
- 79 -
D’après l’expression obtenue ci-dessus, on en déduit que p ( z ) est une fonction
parabolique en z, dont le maximum se trouve en z = 0.
Observons la pression en z = 0 et au temps t = 0 :
23
 18 ⋅ K ( T ) ⋅ F 
23
 l ⋅ L3i 
  ⋅  
p ( z = 0,t = 0 ) = patm +
 4⋅l   18 ⋅ K ( T ) ⋅ F  = patm +
L2i 18 ⋅ F
⋅
3 F
= patm + ⋅
 l ⋅ L3i  4 3 ⋅ l ⋅ Li
3
2 l ⋅ Li
3 ⋅ K ( T ) ⋅  
 18 ⋅ K ( T ) ⋅ F 
F
En posant que la pression équivalente à la force F appliquée s’écrit : péq = , on
l ⋅ Li
obtient l’expression finale :
3
p ( z = 0,t = 0 ) = patm + ⋅ péq
2
6.2.2.3. Détermination de K(T) à partir de L(t)
De l’équation (1.133), on peut extraire K(T) comme suit :
18 ⋅ K ( T ) ⋅ F
L3 (t ) − L3i = ⋅t
l
L3 (t ) − L3i  ⋅ l
K (T) = (1.139)
18 ⋅ F ⋅ t
Sachant que le volume total de matière extraite aux deux extrémités du matériau
au bout d’un temps t s’écrit :
Vext = L (t ) − L i  ⋅ l ⋅ e (t ) (1.140)
on en déduit alors la masse de résine extraite :
- 80 -
mext = ρ ⋅ L ( t ) − L i  ⋅ l ⋅ e ( t ) (1.141)
ei ⋅ L i
avec e ( t ) = (conservation du volume),
L (t )
l
mext = ρ ⋅ L ( t ) − L i  ⋅ ⋅ ei ⋅ L i (1.142).
L (t )
Connaissant la masse initiale mi de matériau : mi = ρ ⋅ ei ⋅ l ⋅ L i

on peut exprimer le rapport rm de la matière extraite :
mext L
rm = =1− i (1.143).
mi L (t )
1
L’équation (1.143) est équivalente à : L ( t ) = ⋅ Li (1.144).
1 − rm
Par conséquent en substituant l’équation (1.144) dans l’équation (1.139) on trouve

l’expression finale de K(T) :
l  1  3 
K (T) =   − 1 ⋅ L i
3
18 ⋅ F ⋅ t  1 − rm  

Plus le rapport rm augmente, c’est-à-dire plus on extrait de matière, plus le

1
terme va augmenter, par conséquent K(T) augmente également. Or K(T) est
1 − rm
inversement proportionnel à la viscosité. Cela signifie que plus K(T) augmente, plus
la viscosité diminue d’où l’augmentation de la quantité de matière extraite. Ce
résultat est bien cohérent avec les résultats obtenus par l’expérience.
6.2.3. 2ème approche : écoulement de matériau fluide de longueur

constante avec évacuation de matière à ses extrémités
On suppose dans ce cas que la longueur du matériau L = L i est constante et
Li
est indépendante du temps, car la matière extraite du milieu au delà de z = ± est
2
évacuée du système.
- 81 -
L’équation d’équilibre des forces (1.121) devient alors :
+ Li 2
F = 2⋅l⋅ ∫ p ( z ) ⋅ dz (1.145).
0
Li
Ecrivons la conservation de la masse en z = :
2
dm
soit la perte de masse en totalité aux deux extrémités :
dt
dm  dz 
= 2 ⋅ ρ ⋅ e (t ) ⋅ l ⋅   (1.146)
dt  dt  z = L i
2
Li
En remplaçant z = dans l’équation (1.125) on trouve :
2
 dz  L i 1 de ( t )
  L =− ⋅ ⋅ (1.147)
 dt  z = i 2 e ( t ) dt
2
et en substituant l’équation (1.147) dans l’équation (1.146) puis en intégrant

dm
l’expression de de 0 à t on obtient la masse extraite ∆m :
dt
dm
t t de ( t )
∆m = ∫ 0 dt
⋅ dt = −ρ ⋅ l ⋅ L i ⋅ ∫
0 dt
(1.148)
6.2.3.1. Détermination de la variation de l’épaisseur e(t)
∆m
L’objectif est de déterminer le rapport de masse rm = avec la masse
mi
de (t )
initiale définie par mi = ρ ⋅ ei ⋅ L i ⋅ l . Il est donc nécessaire de déterminer . Pour
dt
Li
cela nous avons intégré l’équation (1.128) de z à et le résultat est le suivant :
2
- 82 -
1 de ( t )  L i (t ) 
2
1
p ( z ) = patm − ⋅ ⋅ ⋅ − z2  (1.149)
2 ⋅ K ( T ) e (t ) dt  4 
En injectant cette expression dans l’équation (1.145) et en calculant l’intégrale de 0

Li
à on trouve :
2
F 1 1 de ( t ) 2
− patm = − ⋅ ⋅ ⋅ Li (1.150)
l ⋅ Li 12 ⋅ K ( T ) e ( t ) dt
de ( t ) 12 ⋅ K ( T ) ⋅ ∆p
Ainsi =− ⋅ dt (1.151)
dt L2i
F
avec ∆p = − patm = péq − patm (1.152)
l ⋅ Li
correspondant à la différence entre la pression équivalente à la force appliquée sur

la surface ( l ⋅ L i ) et la pression atmosphérique.
L’expression de e ( t ) s’obtient en intégrant l’équation (1.151) de 0 à t :
 12 ⋅ K ( T ) ⋅ ∆p 
e (t ) = e i ⋅ exp  − ⋅t (1.153)
 L2i 
On en déduit que plus le paramètre K ( T ) est grand, c’est-à-dire plus la viscosité est
faible, plus l’épaisseur du matériau décroît, et cela de manière exponentielle. En

effet, la compression d’une résine fluide, perpendiculairement à la direction de
l’écoulement sous la condition h-5), entraîne une évacuation de la résine comprimée
perpendiculairement au sens de l’écoulement, vers les extrémités de l’objet,
entraînant une diminution de l’épaisseur de celui-ci.
De même que plus le temps augmente, plus l’épaisseur du matériau décroît de
 m2 
manière exponentielle. On vérifie que le produit K ( T ) ⋅ ∆p a la dimension   , la
 s 
dimension du facteur de temps t dans l’équation (1.153) est vérifiée. Il faut noter
- 83 -
12 ⋅ K ( T ) ⋅ ∆p 1 
que le terme 2
est de dimension   , ce qui fait intervenir une
Li s 
L2i
constante de temps caractéristique de l’écoulement : τe = fonction du
12 ⋅ K ( T ) ⋅ ∆p
carré de la longueur du matériau, exactement comme la caractéristique temporelle

du transfert de chaleur. Par conséquent, de même que pour le transfert de chaleur,
lors du transfert de matière, une certaine constante de temps dépendant de la
longueur au carré du matériau, doit être prise en compte.
En dérivant l’équation (1.153) et en divisant par e ( t ) on a :
1 de ( t ) 12 ⋅ K ( T ) ⋅ ∆p
⋅ =− (1.154)
e (t ) dt L2i
de (t )
Connaissant l’expression de , nous pouvons calculer la masse de résine
dt
extraite ∆m , en insérant l’équation (1.153) dans l’équation (1.148) puis en intégrant
sur le temps :
  12 ⋅ K ( T ) ⋅ ∆p  
∆m = mi ⋅ 1 − exp  − ⋅ t  (1.155).
  L2i  
En prenant le logarithme népérien de cette équation on en déduit :
L2i  ∆m 
K (T ) = − ⋅ ln 1 −  (1.156)
12 ⋅ ∆p ⋅ t  mi 
∆m
Sachant que rm = , on en déduit la relation entre K(T) et rm :
mi
L2i
K (T) = − ⋅ ln (1 − rm ) (1.157)
12 ⋅ ∆p ⋅ t
- 84 -
Pour connaître la pression p(z) on utilise l’équation (1.149), dans laquelle on insère
l’équation (1.154) ainsi :
6 ⋅ ∆p  L2i 
p ( z ) = patm + 2
⋅  − z2  (1.158)
Li 4 
qui donne pour le cas particulier z = 0 :
3
p ( z = 0 ) = patm + ⋅ ∆p (1.159)
2
Li
La valeur moyenne de p ( z ) , soit p ( z ) de z = 0 à z = , s’écrit :
2
2 Li 2 2 Li 2  6 ⋅ ∆p  L2i 
p (z ) = ⋅ ∫ p ( z ) ⋅ dz = ⋅ ∫ patm + 2
⋅  − z 2   ⋅ dz (1.160)
Li 0 Li 0
 Li 4 
soit encore : p ( z ) = patm + ∆p (1.161).
Or ∆p = péq − patm d’après l’équation (1.152), ainsi on en déduit que :
p ( z ) = péq (1.162)
La pression moyenne au sein du fluide est bien identique à la pression équivalente

à la force appliquée sur la surface ( l ⋅ L i ) .
- 85 -
6.2.4. 3ème approche : écoulement de matériau fluide composé de

fibres de longueur constante avec évacuation de matière à
ses extrémités
Dans cette troisième approche, du fait de la considération de la présence des

fibres de verre, des hypothèses supplémentaires à celles énoncées de h-1) à h-10)
doivent être posées :
h-11) soit fi le taux de matrice initial du matériau de longueur L = Li
h-12) soit Vi = 1 − fi le taux de fibre initial
h-13) la fraction volumique de matrice f ( t ) est variable selon le niveau de
charge appliquée donc selon l’épaisseur e ( t ) également
dz
h-14) est la vitesse de déplacement de la matrice
dt
h-15) la vitesse de déplacement des fibres est nulle
h-16) le volume de fibres est constant quel que soit z
e (t )
l
L z=0 L
− +
2 2
Figure 1 - 20 : Configuration du matériau composite.
Pour une section de cote z l’hypothèse h-16) se traduit par :
e (t ) ⋅ l ⋅ 1 − f (t ) = cons tan te (1.163)
- 86 -
En dérivant cette dernière équation, on trouve la relation suivante :
1 de ( t ) 1 df ( t )
⋅ = ⋅ (1.164)
e (t ) dt 1 − f ( t ) dt
Sachant que, pour un volume constant de matrice compris entre les cotes
z = 0 et z, on a, à cause de la conservation de la matière :
e ( t ) ⋅ z ⋅ l ⋅ f ( t ) = cons tan te (1.165)
et si l’on dérive l’équation (1.165) par rapport au temps :
df (t ) dz de (t )
e (t ) ⋅ z ⋅ + e (t ) ⋅ f (t ) ⋅ + z ⋅ f (t ) ⋅ =0 (1.166).
dt dt dt
Introduisons maintenant l’équation (1.164) sous la forme suivante :
de ( t ) de ( t ) df (t )
f (t ) ⋅ = − e (t ) ⋅ (1.167)
dt dt dt
dans l’équation (1.166), cela donne :
dz de ( t )
e (t ) ⋅ f (t ) ⋅ +z⋅ = 0 soit,
dt dt
1 dz 1 de ( t )
⋅ =− ⋅ (1.168).
z dt e ( t ) ⋅ f (t ) dt
En tenant compte de l’équation (1.164) on aboutit à :
1 dz 1 df ( t )
⋅ =− ⋅ (1.169)
z dt f (t ) ⋅ 1 − f (t ) dt
Cette dernière équation va maintenant nous permettre de s’accéder à p(z).
- 87 -
A partir de l’équation (1.120) et en utilisant l’équation (1.169), on trouve :
dp 1 1 df (t )
= ⋅ ⋅ ⋅z (1.170).
dz K ( T ) f ( t ) ⋅ (1 − f ( t )) dt
1 1 df (t )
En posant G (t ) = ⋅ ⋅ (1.171)
K ( T ) f ( t ) ⋅ (1 − f ( t )) dt
on remarque que cette quantité est indépendante de z puisque la fraction

volumique de matrice f(t) est la même quelle que soit z. Ainsi en intégrant l’équation
Li
(1.170) entre z et , tout en tenant compte de la relation (1.171), on obtient de
2
même :
1  L2 
p ( z ) = patm − ⋅ G (t ) ⋅  i − z 2  (1.172)
2 4 
6.2.4.2. Détermination de la variation de la fraction

volumique de matrice f(t)
Pour déterminer f ( t ) , on suit le développement strictement identique à celui
effectué à partir de l’équation (1.149) en tenant compte de l’équation (1.121) :
F 2 1 L2i 1 
= ∫0L i  patm − ⋅ G ( t ) ⋅ + ⋅ G (t ) ⋅ z 2  ⋅ dz (1.173)
2⋅l  2 4 2 
et l’on aboutit à :
F 1 2
= patm − ⋅ L i ⋅ G ( t ) soit,
l ⋅ Li 12
12 ⋅ ∆p
G (t ) = − (1.174)
L2i
- 88 -
compte tenu de la relation (1.152). D’après l’équation (1.174), on en déduit que la

quantité G ( t ) est indépendante du temps.
En remplaçant G(t) par son expression donnée (1.171) :

1 df ( t ) 12 ⋅ K ( T ) ⋅ ∆p
⋅ =− (1.175)
f ( t ) ⋅ 1 − f (t ) dt L2i
et en écrivant que :
1 1 1 1 1
= + = −
f ( t ) ⋅ (1 − f ( t )) f ( t ) (1 − f (t )) f ( t ) (f ( t ) − 1)
on obtient l’équation :
df ( t ) df ( t ) 12 ⋅ K ( T ) ⋅ ∆p
− =− ⋅ dt (1.176)
f (t ) f (t ) − 1 L2i
L2i
et on remarque que le terme est égal à τe , τe étant la constante de
12 ⋅ K ( T ) ⋅ ∆p
temps caractéristique de l’écoulement introduite précédemment.
En intégrant l’équation (1.176) de 0 à t,
 f ( t )  fi − 1 t
ln  ⋅ =−
 f (t ) − 1 fi τe
puis en prenant l’exponentielle de chaque membre de cette égalité on trouve :
f (t ) fi  t 
⋅= ⋅ exp  −  (1.177)
f (t ) − 1 fi − 1  τe 
ou encore :
f (t ) − 1 fi − 1  t 
⋅= ⋅ exp  
f (t ) fi  τe 
d’où l’expression finale de f(t) :
- 89 -
1
f (t ) = (1.178)
 1 − fi   t 
1+   ⋅ exp  
 fi   τe 
fi
On constate que, si t = 0, alors f ( t ) = = fi et donc l’hypothèse h-10)
fi + (1 − fi )
est bien vérifiée.

De même, si t → ∞ , alors d’après l’équation (1.178) f ( t ) → 0 . En effet, pour un
temps long, la matrice totalement fluide va s’écouler sous l’effet de la contrainte de

pression appliquée, pour être ensuite évacuée aux deux extrémités de l’objet, d’où la
diminution de la fraction volumique de matrice dans l’objet. Ainsi on peut conclure
que ce résultat est bien en accord avec l’expérience.
6.2.4.3. Détermination de la vitesse d’écoulement
dz
Avant de déterminer la vitesse d’écoulement , nous allons dans un
dt
1 dz
premier temps calculer ⋅ . Pour cela, on reprend l’équation (1.169) dans
z dt
laquelle on injecte l’équation (1.175) :
1 dz 12 ⋅ ∆p ⋅ K ( T )
⋅ = (1.179)
z dt L2i
ainsi on en déduit l’expression de la vitesse d’écoulement :
dz 12 ⋅ ∆p ⋅ K ( T ) z
= 2
⋅z = (1.180)
dt Li τe
On peut constater que la vitesse d’écoulement de la matrice ne dépend pas

dz
du temps, mais uniquement de z. De plus, est une fonction linéaire croissante
dt
Li
de z, telle que sa valeur soit nulle en z = 0 et maximale pour z = ± . On peut donc
2
- 90 -
en conclure que la vitesse d’écoulement de la matière est nulle au centre de

l’éprouvette et est maximum aux extrémités de l’objet.
6.2.4.4. Détermination de la quantité de matière évacuée
dm L
La détermination de la quantité de matière évacuée en z = ± i passe
dt 2
par la détermination de e ( t ) . En effet, l’équation (1.165) fournit le volume total de
Li L
matrice qui s’est écoulée en z = − et z = + i . Donc la quantité totale de matrice
2 2
Li
évacuée est égale au double de la quantité qui s’est évacuée en z = + :
2
dm  dz 
= 2 ⋅ ρ ⋅ e (t ) ⋅   ⋅ l ⋅ f (t ) (1.181).
dt  dt  z =+ L i 2
En égalant les équations (1.168) et (1.179) :
1 de ( t ) 12 ⋅ ∆p ⋅ K ( T ) 1
− ⋅ = = (1.182)
e (t ) ⋅ f (t ) dt L2i τe
et compte tenu de l’expression suivante de f ( t ) :
fi
f (t ) = (1.183)
 t 
fi + (1 − fi ) ⋅ exp  
 τe 
issue de l’équation (1.178), l’équation (1.182) devient :
1  t 
− ⋅ exp  − 
de ( t ) τe  τe 
= ⋅ e (t ) (1.184).
dt  t  1 − fi
exp  −  +
 τe  fi
- 91 -
En intégrant cette équation de 0 à t d’une part, puis en prenant l’exponentielle du

résultat obtenu, on trouve enfin :
   t  
e ( t ) = ei ⋅ 1 − fi 1 − exp  −    (1.185).
 
   τe   
On remarque que :
si t = 0, alors e ( t ) = ei , l’hypothèse h-2) est bien vérifiée
si t → ∞ , alors e (t ) → e i ⋅ (1 − fi )
dz
Connaissant e ( t ) , nous allons pouvoir calculer le terme e (t ) ⋅ f (t ) ⋅ grâce à
dt
l’équation (1.168) :
dz de ( t )
e (t ) ⋅ f (t ) ⋅ =− ⋅z (1.186).
dt dt
de ( t ) 1  t 
Or d’après l’équation (1.185) : = − ⋅ ei ⋅ fi ⋅ exp  −  (1.187)
dt τe  τe 
d’où en remplaçant dans l’équation (1.186) :
dz 1  t 
e (t ) ⋅ f (t ) ⋅ = ⋅ ei ⋅ fi ⋅ exp  −  ⋅ z (1.188).
dt τe  τe 
Li
En substituant l’équation (1.188) dans l’équation (1.181), on obtient en z = :
2
dm t L  t  mi  t 
= 2 ⋅ l ⋅ ρ ⋅ ⋅ ei ⋅ fi ⋅ i ⋅ exp  − = ⋅ exp  −  (1.189).
dt τe 2  τe  τe  τe 
En intégrant l’équation (1.189) de 0 à t, on obtient la masse totale extraite ∆m ( t )
au bout du temps t :
- 92 -
  t 
∆m ( t ) = mi ⋅ 1 − exp  −  (1.190).
  τe  
∆m
Comme indiqué précédemment : rm = (avec mi la masse initiale du matériau),
mi
on déduit de l’équation (1.190) :
 t 
rm = 1 − exp  −  (1.191).
 τe 
6.2.4.5. Détermination de la caractéristique d’écoulement K(t)
L’équation (1.191) s’écrit également de la façon suivante :
 t 
exp  −  = 1 − rm (1.192)
 τe 
En prenant le logarithme népérien de l’équation (1.192) on obtient :
t
− = ln (1 − rm ) (1.193)
τe
L2i
et en remplaçant τe par l’expression qui suit : = τe , on trouve K ( T ) :
12 ⋅ K ( T ) ⋅ ∆p
L2i
K (T) = − ⋅ ln (1 − rm ) (1.194)
12 ⋅ ∆p ⋅ t
Pour avoir une idée de l’ordre de grandeur de la caractéristique d’écoulement

K ( T ) , nous avons tracé la fraction de résine extraite rm , issue de l’expérience
d’essorage décrit précédemment, pour différentes températures isothermes allant de

60°C à 120°C (figure (1-19)). Les valeurs de K ( T ) sont calculées à partir de
l’équation (1.194) pour les conditions suivantes :
- 93 -
! L i = 0.3 m
! t = 300 s
! ∆p ≈ 9 ⋅ 105 Pa
! fraction massique de matrice égale à 29 % du matériau composite.
4E-12
3.5E-12
3E-12
K (T) (m²/Pa.s)
2.5E-12
2E-12
1.5E-12
1E-12
5E-13
0
60 70 80 90 100 110 120
Température (°C)
Figure 1 - 21 : Variation de la caractéristique d’écoulement K (T) pour chaque température

d’essorage.
Rappelons que la signification physique de K(T) doit être limitée aux

températures inférieures ou égales à 100°C.
7. CONCLUSION
Ce chapitre traite de l’analyse du procédé d’élaboration de matériaux

composites, fondé sur l’hystérésis diélectrique, en trois parties. Avant de nous
intéresser au précédé d’élaboration proprement dit, nous avons décrit le matériau
sujet de l’élaboration. Dans notre cas il s’agit d’un composite à matrice
thermodurcissable (DGEBA/DDA) associée à des fibres de verre unidirectionnelles,
c’est pourquoi une étude de ces deux constituants de base, résine et renfort, a été
menée.
- 94 -
Dans un premier temps, nous nous sommes donc intéressés aux matériaux
composites ainsi qu’aux différents procédés de mise en œuvre de ces matériaux et à
leurs caractéristiques, ce qui nous a amenés à classer ces procédés suivant deux
classes, d’une part les procédés dits « classiques » et de l’autre les procédés fondés
sur l’hystérésis diélectrique. La différence entre ces deux procédés provient du
mode de chauffage.
La description d’un tel procédé repose sur trois aspects : l’aspect
électromagnétique, l’aspect cinétique chimique et enfin l’aspect écoulement de la
matière.
Concernant l’aspect électromagnétique, une comparaison de la propagation
électromagnétique en espace libre ou dans un guide a été rappelée. Puis nous avons
effectué une analyse de la distribution spatiale du champ électrique, à l’origine des
sources internes de chaleur, au sein de l’applicateur chargé par des matériaux
diélectriques et comportant en son centre le matériau composite à traiter. Cette
analyse est fondée sur la résolution des équations de Maxwell par des calculs
formels, à l’aide du logiciel MapleV, pour un applicateur à remplissage aussi bien
symétrique que non symétrique. Les résultats obtenus pour une propagation en
mode TE0,1 sont satisfaisants puisqu’ils montrent bien une variation suivant une
demi-sinusoïde de l’amplitude du champ électrique selon l’axe y, et une amplitude

constante du champ électrique suivant x.
La matrice du matériau composite étant à base de résine époxyde et de

durcisseur aminé, entraînant une forte réaction exothermique au cours du procédé,
une étude de l’aspect cinétique chimique de cette réaction s’est imposée. Cette
analyse nous a permis de connaître l’évolution de la conversion chimique au cours
de la transformation du matériau, cette dernière ayant sa contribution sur la
distribution des sources internes de chaleur.
La dernière partie de ce chapitre concerne la description de l’écoulement de

la matière fluide au cours de sa mise en forme. L’objectif est de mettre en évidence
les grandeurs caractéristiques de l’écoulement. Pour décrire cet écoulement, nous
avons donc créé des situations d’écoulement en isotherme en deux étapes.
D’une part, sur un plan expérimental, nous avons reproduit l’écoulement de
la matrice du matériau composite avec l’expérience appelée « opération d’essorage ».
Ces résultats montrent que la fraction de résine extraite pour différentes
- 95 -
températures isothermes peut atteindre 7 % de la masse totale du matériau

composite.
D’autre part, sur un plan théorique, nous avons étudié des écoulements
particuliers, allant d’écoulement comportant des hypothèses très simples, jusqu’à
des écoulements avec des hypothèses plus complexes. La prise en considération des
configurations d’écoulement à difficulté croissante nous a permis non seulement de
connaître la distribution parabolique de la pression au sein du matériau qui
s’écoule, mais également de décrire l’écoulement réel.
Enfin la confrontation des résultats expérimentaux de la quantité de matière
extraite avec les résultats des calculs formels nous a permis de connaître la
constante caractéristique de l’écoulement K(T) pour différentes isothermes.
Cependant la température n’étant pas constante au cours du procédé d’élaboration,
il nous restera donc à voir la variation de K(T) avec la température. A cet effet, une
discussion portant sur cette grandeur caractéristique de l’écoulement sera menée à
la fin du chapitre III, une fois connue la distribution spatiale de la température au
sein du matériau.
- 96 -

Elaboration

Transféré par

Droits d'auteur :

Formats disponibles

Elaboration

Transféré par

Informations du document

Description originale:

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

Elaboration

Transféré par

Droits d'auteur :

Formats disponibles

CHAPITRE I : DESCRIPTION DU PROCEDE

D’ELABORATION SOUS MICRO-ONDES

La maîtrise de l’élaboration de matériaux composites à matrice polymère

De plus, le procédé s’accompagne de l’application d’une contrainte de

2. LES MATERIAUX COMPOSITES

Il existe deux types de matériaux composites :

• les composites grande diffusion (GD)

Les composites GD, représentent 95 % des matériaux composites utilisés. En

Les composites HP sont principalement utilisés dans l’aéronautique. Les

La maîtrise de l’élaboration de matériaux composites nécessite la

Figure 1 - 1 : Constituants des matériaux composites.

2.1. Les renforts

Le renfort a pour rôle d’assurer la tenue mécanique du composite et est

Les fibres sont constituées par plusieurs centaines ou milliers de filaments de

Les différents types de renforts sont cités dans l’organigramme ci-dessous.

Artificiels Végétaux Céramiques Métalliques

Polyesters Aramides Bois Coton Verre Carbone Bore

Figure 1 - 2 : Les différents types de renfort de base.

Le tableau suivant présente les principales caractéristiques mécaniques de

Fibre Densité Charge de Charge de Allongement Module Diamètre

Tableau 1 - 1 : Caractéristiques de quelques renforts.

2.1.2. Les fibres de verre

Le verre possède une très bonne résistance à la rupture dépassant même

Il existe différents types de verre : le verre E, C, D, R ou S, qui sont obtenus

Le verre E est employé pour la fibre de verre constituant le renfort du

2.2. Les matrices

Thermodurcissables Thermoplastiques Elastomères

Borures Carbures Nitrures

Figure 1 - 3 : Différentes familles de matrices.

Une comparaison entre quelques caractéristiques des résines

Critères essentiels Matrices thermoplastiques Matrices thermodurcissables

Tableau 1 - 2 : Caractéristiques des résines thermodurcissables et des résines

Il existe d’autres catégories de résines : les résines thermostables et les

2.3. Les matériaux composites structuraux

Les matériaux composites structuraux sont classés en trois catégories :

Les monocouches représentent l’élément de base de la structure composite.

Les stratifiés sont constitués d’un empilement de monocouches ayant

Figure 1 - 4 : Composite stratifié.

Le choix de l’empilement et plus particulièrement des orientations permettra

résistance. L’ensemble forme une structure d’une grande légèreté. Le matériau

Figure 1 - 5 : Composite sandwich.

Les matériaux composites utilisés dans le cadre de ce travail sont

3. ELABORATION DES MATERIAUX COMPOSITES A MATRICE

Les procédés de mise en œuvre de matériaux composites nécessitent trois

1) une étape d’association : l’imprégnation du renfort par la résine

Mise en place du mélange sur

Figure 1 - 6 : Principales étapes d’un procédé par moulage.

La constitution du moule va dépendre du choix de la technique employée. Le

! la forme et la taille des pièces à élaborer, car certains procédés se limitent

! la nature de la matrice du matériau composite, qui doit être adaptée à la

! le nombre de pièces à élaborer, qui détermine le coût de l’outillage. Pour

! la qualité mécanique des pièces : les pièces à fort taux de fibre ne

3.1. Procédés classiques et procédés micro-ondes

Il existe une multitude de procédés classiques donc nous citerons

3.1.1. Procédés classiques

Le problème qui nous concerne est la mise en place de la matrice, qui

Contrairement aux procédés à faible pression, les procédés à forte pression

1) Les procédés à faible pression à froid :