Abaqus
Abaqus
Abaqus
z
T
k
y
T
k
x
T
k
t
T
c
d
dz
d
dy
d
dx
d
d d
Les conditions de frontire peuvent tre indiqus comme suit :
d f d d
d
dz
q m g T T h m g
z
T
k ) ( ) ( )] ( 1 [
1
+
'
p l a q u e t t e l a d e d o m a i n e d u h o r s
p l a q u e t t e l a d e d o m a i n e l e d a n s
m g
0
1
) (
0
z
T
k
d
dz
) (
3 f d d y
d
dy x
d
dx
T T h n
y
T
k n
x
T
k
y
d
dy x
d
dx
n
y
T
k n
x
T
k
q
d
est le flux de la chaleur entrant dans le disque au niveau de la zone de contact, et h
d1
et
h
d3
sont, respectivement, les coefficients de convection hors de la zone de contact et de la
surface du disque.
L'analyse thermique du systme de freinage exige une dtermination prcise de la
chaleur totale de friction produite et ainsi que la distribution de cette nergie entre le
disque et les garnitures. Lors bref dun freinage durgence, toute la chaleur produite
l'interface est gale la chaleur absorbe par le disque et les garnitures.
Ainsi,
q
d
(r,t)=p(r,t)v(r,t)=p(t)(t)r ,
est le coefficient de frottement et est le facteur de dispersion de chaleur reprsentant
la quantit du flux thermique de friction entrant dans le disque, qui est donn comme suit :
d d d p p p
k c k c
/ 1
1
+
,
_
,
_
m
t
t
p t p
exp 1 ) (
0
,
est un paramtre qui caractrise la dure de l'augmentation de la pression de frein de
zro la valeur nominale p
0
.
Lorsquun vhicule freine, une partie de la temprature de friction schappe dans
l'air grce la convection et au rayonnement. Par consquent, la dtermination des
coefficients de convection est trs importante. Il est, cependant, trs difficile de les calculer
9
exactement, car ils sont affects par la forme du systme de freinage, de la vitesse de
dplacement du vhicule et donc, par consquent, de la circulation de l'air.
II.2.b Flux dair en mcanique des fluides
Lcoulement dans le sillage des vhicules automobiles peut-tre spar en deux
structures. La premire issue dune ligne de sparation bidimensionnelle gnre un
coulement de base avec une pression statique constante et une vitesse de transport
voisine de zro. La seconde et issue dune sparation tridimensionnelle de lcoulement et
engendre des zones enroulements tourbillonnaires.
Ces deux effets produisent la trane de
forme laquelle sajoute une trane de pression
rsultant de la distribution de pression sur la
gomtrie du vhicule et une trane de
frottement due au cisaillement du fluide dans la
couche limite. Lensemble constitue la trane
arodynamique globale associe au vhicule.
En arodynamique automobile, les
coulements fluides sont de nature turbulente et
dans les conditions courantes dutilisation, le
dplacement du vhicule seffectue un nombre
de Mach trs infrieur 0.3.
Figure 4 : Flux dcoulement dair sur une moto
De plus, les actions mcaniques de gravit que le fluide exerce sur le vhicule sont
faibles compars aux actions mcaniques de viscosit, de turbulence ou de pression.
Pour ces raisons le fluide est considr incompressible et non pesant.
Pour calculer la trane arodynamique
sous laspect tourbillonnaire on va
sintresser une formulation intgrale
construite sur un bilan de quantit de
mouvement partir des pressions
moyennes au sens de la turbulence.
Formulation thorique :
+
,
_
,
_
s
i i
s
y
s
x
x
d P P
V
d
V
V
d
V
V V
F ] ) (
1 [
2
0
0 0
2
0
0
1
er
terme : Contribution due la vitesse longitudinale
2
me
terme : Contribution due la rotation du fluide
3
me
terme : Contribution due la perte de pression darrt
Coefficient local de frottement Coefficient de trane :
10
Le coefficient local de frottement est dfini partir de la contrainte paritale locale
norme par la pression dynamique soit :
( )
( )
2
2
1
E
P
f
U
x
x C
U S
D
C
D
i
i
x
u
,
_
j
i
j i j
i
j
i
x
u
x x
p
p x
u
u
t
u
,
_
i i i
i p p
x
T
x x
T
u c
t
T
c
La premire quation est lquation de continuit pour un fluide incompressible, la
seconde reprsente lquation de Navier-Stockes pour un fluide Newtonien et enfin
lquation de lnergie pour un petit nombre de Mach et dEckert.
Dans notre cas, lcoulement dair (inflow) est
laminaire jusquau moment o il rencontre le disque.
Une partie de ce flux devient turbulent.
La turbulence est un mode naturel dcoulement
dun fluide visqueux o des mcanismes internes
dchange dnergie assurent la cration et le
maintien de toute une hirarchie de mouvements
chaotiques rpartis continment sur une large gamme
dchelle macroscopiques.
Figure 5 :
II.2.d Mthode de rsolution transient
Pour rsoudre le problme, la technique du transient est employe. Pendant la
simulation de freinage, les diffrentes valeurs du flux thermique sont assignes aux
lments dpendants de la zone de contact chaque pas de temps dans le modle.
11
La localisation de la rgion de contact, chaque pas de temps, doit tre dtermine
par la position relative du disque et de la plaquette.
De cette faon, la distribution de source de chaleur peut tre transfre au modle,
et le problme de conduction peut tre rsolu comme problme tridimensionnel de
transfert thermique avec une source de chaleur immobile.
Au pas de temps suivant, la source de chaleur se dplace aux lments voisins
selon la direction relative de glissement, alors que le problme de conduction thermique
est, plusieurs reprises, rsolu en utilisant un incrment de temps plus petit que les pas
de temps mentionnes incluant la distribution initiale de la temprature correspondant
exactement au rsultat qui a t obtenu au pas de temps prcdent. Le pas de temps est
rgi par la variation de la vitesse de rotation de disque.
Employer la mthode du transient en lments finis permet de simuler le transfert
thermique dans un disque tridimensionnel avec une vitesse variable et un flux transient
non-axisymmetric de la chaleur.
Afin de modliser lensemble de ces phnomnes (mcanique des fluides et
thermique), on utilise le logiciel Star-CD.
Le prochain chapitre dtaillera donc les principales tapes de cration dun modle
pour ce logiciel et de rsolution du problme.
12
III DETERMINATION THEORIQUE DE LA TEMPERATURE
ATTEINTE AU FREINAGE
Pour le problme prsent nous devons savoir la raison pour laquelle la moto
Mtisse est confronte une usure prcoce et in fine une casse du disque de frein. Tout
dabord nous allons donc nous attacher dterminer de faon thorique la temprature
atteinte lors dun gros freinage et vrifier ainsi la cause de la casse. Le but est aussi de
reporter cette temprature sur le modle que nous devons laborer sur Star-CD.
Ces calculs seront effectus avec 2 types de disque avec des matriaux diffrents.
Dans un premier temps nous valuerons la temprature avec de lacier inoxydable,
matire utilise par la Mtisse, et dans un second temps avec du carbone.
Caractristiques thermiques et tempratures limites dutilisation de lacier
inoxydable et du carbone :
On souhaite prdire lvolution de la distribution de temprature T(x,t) de ce milieu
de conductivit thermique uniforme et isotrope. Lquation de bilan thermique dun
domaine limit par des plans est de la forme :
0
1
2
2
t
T
a x
T
Conditions limites :
( ) ( )
i i
T t T T x T , 0 ,
Les donnes cinmatiques du freinage de la Mtisse permettent den apprcier la
violence : la vitesse de la moto (masse de 200 kg , pilote compris) passe
approximativement de 290 km/h 100 km/h en 120m. La dure de freinage sexprime
donc par
v
L
f
,
avec
2
f i
v v
v
+
,
o v
i
= 80,55 m/s et v
f
= 27,77 m/s .
Soit une dure de freinage
f
de 2,2 s .
La dclration est dfinie par
f
f i
v v
.
Soit une dclration suppose constante de 23,9 m/s
-2
qui quivaut 2,4 g .
13
Acier Inoxydable Carbone pyrolytique
K
A
= 20 W.m
-1
.K
-1
A
= 7900 kg.m
-3
C
A
= 560 J.kg
-1
.K
-1
T
lim,A
= 1670 K
K
C
= 65 W.m
-1
.K
-1
C
= 2200 kg.m
-3
C
C
= 1600 J.kg
-1
.K
-1
T
lim,C
= 2270 K
On suppose que lnergie cintique de la moto est totalement convertie en nergie
thermique par frottement des plaquettes (2 par disques) sur les freins, nous avons donc :
( )
f f f f i c
S v v M E 4
2
1
2 2
,
avec S
f
la surface des plaquette en contact avec le disque soit 0.0035 m
2
.
Lnergie cintique de la moto en dclration est donc denviron 571 712 J.
Nous pouvons alors dterminer le flux thermique produite chaque interface plaquette-
disque, soit
f
= 18,56 MW.m
-2
.
Pour un flux impos
0
nous avons :
( )
t
b
T t T
i
0
2 , 0 + ,
avec b leffusivit thermique du matriau, soit pour chaque matriau 1 et 2 :
t
b
T
t
b
T T
i i f
2
2
1
1
2 2 + + .
La densit de flux thermique
f
se rpartit entre la plaquette et le disque selon deux
composantes respectivement notes
1
et
2
, la temprature finale sera alors donne
par :
f f
i f
b b
T T
2 1
2
+
+
.
Calcul de leffusivit thermique de lacier inoxydable et du carbone, respectivement b1 et
b2 :
2
1
1 2 4
1
. . . 10 . 513 , 1 S K m W C K b
A A A
,
2
1
1 2
2
. . . 9406 S K m W C K b
C C C
.
Avec une temprature initiale T
i
de 60 C soit 333 K nous pouvons dsormais dfinir la
temprature T
f
:
C K T
f
1326 1599
.
Cette temprature finale est donc proche de la temprature critique de lacier
inoxydable (1670K) il en rsulte donc un craquage thermique qui se vrifie sur le disque.
Ce rsultat nest pas tonnant au vue des problmes rencontrs sur la Mtisse, lusure est
alors prcoce et le bris de disque est plus probable sur des circuits ayant de fortes
sollicitations en freinage. Cependant lors de ce calcul thorique nous navons pas pris en
compte le refroidissement par air du disque ventil.
Cependant cette temprature est mettre en relation avec un systme disque-
plaquette en matire carbone-carbone, la temprature finale atteinte est alors de
C K T
f
1086 1359
,
alors que sa temprature critique est bien plus haute, 2270 K.
14
IV RESOLUTION DUN PROBLEME DE THERMIQUE DANS
UN ECOULEMENT
Lutilisation de Star-CD est complexe. Celle-ci se fait par utilisation de plusieurs
modules spcifiques.
Nayant pas accs au modeleur de Star-CD, on doit donc crer la gomtrie par un
autre logiciel, et utiliser un format denregistrement compatible.
IV.1 Cration de la gomtrie
La premire tape consiste
crer la bote qui contient le
systme tudier. Celle-ci doit
tre assez grande pour ne pas
gnrer deffets de bords.
Dans une seule et mme pice on
reprsente le disque et la plaquette dans
la configuration relle (plaquette en
contact avec le disque). Cependant on
cre ces lments dans deux corps de
pices diffrents, afin que Star-CD les
diffrencie.
Attention : avant lenregistrement (au
format STL), il ne faut pas laisser le
rendu filaire du modle mais le
remettre en raliste sinon il y aura
une erreur pendant le maillage.
On utilise donc Catia V5R8 afin de crer la zone ou le fluide va scouler.
On applique enfin lopration boolenne retirer et on obtient donc une gomtrie,
enregistrer au format stl.
15
On utilise donc dans lensemble des calculs, afin de mailler notre gomtrie, le
Mesh Wizard de Pro-Am.
IV.2 Maillage de la gomtrie
IV.2.a Importation de la gomtrie
Dans cette tape, on va simplement chercher par Browse le fichier stl cr par Catia,
afin de le stocker en mmoire.
Ltape suivante consiste importer cette gomtrie dans Star-CD. Voici les
paramtres utiliss par le Mesh Wizard de Star-CD :
Le fichier .stl est importable
par Star-CD dans la
premire tape du Mesh
Wizard.
IV.2.b Prparation de la surface
Dans cette tape, Star-CD va prparer la surface afin de faciliter le maillage.
On va donc choisir des tailles maximales et minimales de surfaces que Star-CD va utiliser
pour mailler la surface, ou plutt pour la prparer au maillage.
Dans cette tape, on supprime donc de la gomtrie tous les dtails infrieur la taille
minimale choisie.
16
La deuxime tape consiste
prparer le maillage. En effet, un
premier maillage est appliqu la
gomtrie. Celui-ci servira de base de
travail pour les oprations suivantes.
Lintervalle de longueurs dlments
est ici dfinie. Dans notre cas les
lments ont des longueurs
comprises entre 200 mm (lments
proches de la bote) et 0.5 mm
(lments proches du disque).
IV.2.c Gnration de la couche limite
On va choisir ici lpaisseur de la couche limite, qui sera maill indpendamment
par extrusion du maillage de surface sur une maille, et qui est donc enlever du modle
courant pour les deux tapes suivantes.
Dans la troisime tape, le Mesh Wizard
gnre la couche limite du modle.
17
IV.2.d Gnration de la finesse du maillage
On choisit ici des valeurs maximales et minimales entre lesquelles le mailleur
automatique va grer la finesse du maillage.
On obtient donc la fin de cette tape une cartographie pour lensemble de la
structure de la finesse du maillage qui sera utilise dans ltape suivante.
Cette tape sert dfinir la finesse
du maillage. Ces valeurs dpendent des
longueurs dlments dfinies dans la
deuxime tape. De plus, le Mesh
Wizard gnre ici une cartographie du
maillage utile pour les prochaines
tapes.
IV.2.e Cration du maillage
Dans cette tape, partir de la cartographie de finesse du maillage, le mailleur va
mailler le volume par des mailles de type ttradrique (tetrahedral) ou cubique (trimmed)
au choix de lutilisateur.
Le maillage dfinitif peut maintenant tre
gnr, grce la cartographie obtenue. Deux
types de maillages sont disponibles : Trimmed
(prismes rguliers) et Tetrahedral. Dans notre
cas, un maillage de forme ttradrique est plus
appropri la gomtrie utilise.
18
IV.2.f Extrusion de la couche limite
Ici, le mailleur va reconstituer la couche limite par extrusion du maillage de surface,
afin dobtenir des mailles dpaisseurs gales lpaisseur de la couche limite.
Dans cette tape, le Mesh
Wizard va extruder les volumes entre
les couches limites. Cela a pour but
de diffrencier la bote et la pice, afin
dappliquer un fluide dans les zones
extrudes.
IV.2.g Statistiques
On vrifie ici lensemble des paramtres gnrs par le mailleur automatique, cest
dire le nombre de mailles de fluide, le nombre de mailles gnres par lextrusion lors de
ltape prcdente, ainsi que le nombre de mailles rats .
IV.2.h Vrification du maillage
Le mailleur automatique ntant pas parfait, il est ncessaire de vrifier le maillage
obtenu par le Star-Guide.
On utilise donc le check grid , avec comme option de vrification : severe
check . On obtient alors un listing des vrifications dans la fentre texte de Pro-Am.
IV.2.i Correction du maillage
Maintenant que le maillage est vrifi, il convient de supprimer les mailles
dfectueuses. Pour cela, on va lister celles-ci par le menu lists , et cells . Une
fentre saffiche, listant lensemble des mailles dfectueuses. On les supprime donc, et on
fait de mme pour les couples.
Mais cette opration peut crer de nouveaux problmes dans le maillage. On doit
donc itrer les oprations 9 et 10 jusqu lobtention dune fentre vide dans les deux
listing.
19
IV.3 Modle physique
Dans cette tape, on va renseigner lensemble des caractristiques physique du
fluide, du solide ainsi que des parois.
Le maillage tant effectu, il faut dfinir tous les paramtres des diffrents modles afin de
lancer lanalyse.
IV.3.a Etat stationnaire
Dans un premier temps il faut valider les
modles choisis, donc une analyse en
stationnaire est suffisante : Steady State.
Il ne faut pas oublier dactiver le calcul du transfert de
chaleur : Heat Transfer On dans Thermal Options.
On a ainsi accs au menu Solids qui regroupe les
paramtres des matriaux non fluides.
Vrifier que la gravit sexerce dans la bonne direction
et le bon sens : ici suivant laxe z.
20
IV.3.b Initialisation de la turbulence
Le type de dcoulement doit ensuite tre dtermin. Ici
on choisit un coulement turbulent de type k- en utilisant le
nombre de Reynolds. Ce type dcoulement est utilis
frquemment pour de nombreux problmes. Ceci nous permet
dobserver distinctement une turbulence autour du disque.
Cependant, en ajoutant les paramtres de la thermique nous
avons connus des problmes en lanant lanalyse, donc nous
avons utilis une turbulence constante avec une viscosit de
0.005, bien moins capricieuse pour le calcul.
Vrifier tout dabord que le calcul de la temprature de
lair est activ.
Ensuite initialiser des paramtres de la turbulence : les
deux paramtres ont t rgls 0.5.
21
IV.3.c Dfinition du modle matriau
Tout dabord il est ncessaire dassigner les
matriaux solides utiles au modle. Ainsi pour le disque,
nous avons utilis un acier se trouvant dans la
bibliothque.
Pour simplifier ltude, nous avons utilis ce mme
acier pour la plaquette de frein.
Nous avons donc deux matriaux en tout : lacier pour
les solides et lair pour le fluide.
Pour appliquer ces matriaux aux cellules solides et fluides, il
faut dfinir les types de cellules existantes, ainsi que les matriaux
associs dans Cell Table Editor .
Par dfaut, toutes les cellules sont fluides.
Cependant le contour de la bote et les pices au
milieu ne sont pas classs dans la mme table que
les cellules intrieures de la bote, ce qui explique
la diffrence de couleur.
22
Pour appliquer le
matriau aux pices, il faut
utiliser le bouton Modify
Type dans le Cell Tool, en
ayant pralablement rendu
actif le type de cellules voulu.
On slectionne donc les
cellules correspondantes.
Avec le mode de
slection select zone , des
cellules appartenant la bote
sont parfois slectionnes, il
faut donc changer leur type et
les remettre fluides (ne surtout
pas le supprimer).
Aprs cette opration il est trs probable que certaines cellules soient communes
au disque et la plaquette, et parfois certaines cellules existent en double. Afin de corriger
ce problme, il faut se rendre dans Check Grid dans le Star-Guide puis dans le menu
droulant Check option slectionner severe checks , en vrifiant que put failed
cells into new cell set soit sur yes . Ensuite, dans le menu droulant Lists , puis
Cells apparaissent seulement les cellules critiques (sil y en a), elles doivent tre
supprimes. Ritrer le procd jusqu ce quil ny ait plus aucune cellule dfaillante dans
la liste.
23
IV.3.d Dfinition du modle Conditions Limites
La premire tape consiste slectionner les faces dentre et de sortie du flux de
la bote. Ces options se trouvent dans le menu Create Boundary dans le Star Guide.
La face dentre du flux (Inlet).
La face de sortie du flux (outlet).
Des conditions aux limites sont galement ncessaires sur les deux pices. Nous
avons utilis les conditions Wall . Cela est obligatoire car on peut rgler un certain
nombre de paramtres sur ces C.L. , comme une vitesse de rotation sur le disque, ou
une source de chaleur ainsi que les paramtres de radiation thermique sur la plaquette.
24
IV.3.e Chargements
Lchauffement d au frottement de la plaquette sur le disque est trs compliqu
simuler, cest pourquoi une temprature initiale de 600 K est applique directement sur la
plaquette. Nous avons donc un problme de conduction thermique entre deux solides.
Dans le menu Define boundary regions il est possible de rgler les diffrents
paramtres propres chaque condition.
Dfinition des paramtres :
Inlet : on applique une vitesse U= -58 m.s
-1
pour le flux.
Wall (disque) : on applique une vitesse de rotation
=2000 tr.min
-1
dans le repre n4. On fixe une
temprature initiale de 293 K ainsi quune rsistance de
20 mK/W.
Pour crer un nouveau repre, cliquer sur le bouton
Csys , puis garder les mmes paramtres que pour
le repre 1 en ajoutant une rotation Rzx de 90, puis
cliquer sur New (global) .
Les vitesse de rotation du disque a t pose, pour un
diamtre de la roue de 70 cm, cela donne une vitesse
de 208 km.h
-1
pour le vhicule (ce qui semble raliste),
do une vitesse du flux de 58 m.s
-1
.
Wall (plaquette) : Il faut changer les paramtres
thermiques de cette condition. Dans Wall heat, on choisit Fixed avec pour valeur 600 K
Pour les deux solides les valeurs de lmissivit et de la rflexion thermique sont
respectivement de 0.8 et 0.2.
Radiations :
Ce paramtre doit tre activ afin de visualiser
linfluence de la temprature des pices sur le flux dair.
Cette partie est une des plus dlicate car plusieurs
mthodes de rsolutions sont proposes, il faut donc choisir
la plus approprie notre modle.
25
Dans notre cas, nous avons utilis la mthode
Discrete Transfer Internal VF Calc dans Thermal
Options . Nous avons rgl Beams per patch sur 1600.
Avec cette mthode, Star-CD nous impose dappliquer
des patches sur les conditions aux limites dans longlet
appropri qui se trouve dans Create Boundaries .
Attention, il ne faut pas oublier dactiver le calcul des
radiations pour le solide dans le Cell Table Editor.
(Radiation sur On ).
Avant de lancer une analyse (parfois trs longue), on doit bien sassurer que le
modle ne contient pas derreurs. Cest pourquoi il faut vrifier tout le modle.
Dans Check Model Setup puis Check
Everything , cliquer sur All . Une srie de tests
va contrler les cellules, les conditions limites
On peut suivre en temps rel cette srie de
tests dans la fentre pro-STAR Output . Cela
peut prendre un peu de temps. Le test se termine
par la vrification des rgions, et doit en cas de
succs afficher Region Checking Completed
NO ERROR .
26
IV.3.f Analyse
Autres rglages :
Pour obtenir un meilleur rsultat au niveau du
calcul de la temprature, rgler loption solution
algorithm sur PISO au lieu de SIMPLE . cet
algorithme de calcul fait intervenir plusieurs
correcteurs afin daugmenter la prcision du rsultat.
Afin de visualiser la chaleur, il est ncessaire de
cocher les cases post et print pour ces
paramtres dans le menu Analysis Output .
27
Analyse Stationnaire :
On choisit le nombre ditrations que le
calculateur va effectuer. Plus le nombre est lev, et
plus le rsultat sera prcis, mais galement plus long
calculer.
On dtermine ensuite le facteur dchelle voulu
(0.001 pour les millimtres), et en slectionnant le mode
Double Precision on accrue la fiabilit du rsultat.
On lance lanalyse et on peut suivre la progression et les
messages derreurs dans le STAR-Launch .
28
Pour les premires analyse nous avons prfr visualiser lcoulement de lair
autour du disque. En effet il serait inutile de raliser une analyse complte si lallure du flux
ne semble pas raliste ou si celle-ci provoque des erreurs lors du calcul. Il faut ainsi
procder tape par tape pour connatre les sources derreurs.
Dans la partie Post-
Processing nous avons donc affich
la vlocit totale du flux,
reprsente par des vecteurs. On
distingue nettement le champ de
turbulence cr par la rotation du
disque.
Pour accentuer
cette interprtation, nous
avons ajout un ruban sur
la face inlet. Les rsultats
confirment les
observations effectues
prcdemment. En effet,
le flux contourne le disque
et possde plusieurs
tourbillons.
29
(En ajoutant le chargement thermique dcrit en IV.3.e, le rsultat obtenu nest pas
trs satisfaisant car lchauffement du disque est local, au niveau du contact avec la
plaquette.) JE SAIS PAS SI CEST LE BON MOMENT POUR PARLER DE CA
Donc ici, nous avons appliqu le chargement thermique directement sur le disque
afin dobserver les consquences sur la temprature du flux.
On remarque que le
disque ne subit pas
beaucoup linfluence du flux
dair car il reste
uniformment la mme
temprature. Cependant il
est intressant de voir que
le flux, lui, subit bien la
radiation thermique du
disque, et affiche ainsi une
trane qui a le mme sens
que lcoulement.
30
Analyse transient :
Aprs avoir slectionn le mode Transient dans
Select Analysis Features nous devons dfinir le
domaine temporel tudi.
31
32
Conclusions :
IV.4 Lancement du calcul
Aprs vrification du modle et des conditions aux limites, on lance enfin le calcul dans
cette fentre, avec la bonne unit.
IV.5 Post-traitement
IV.5.a Chargement du fichier de rsultat
Les rsultats sont stocks dans un fichier spcial, il convient donc de le charger en
mmoire.
IV.5.b Chargement des donnes voulues
Dans cette tape, on charge pour laffichage en vecteur et en couleur les donnes voulues
du fichier rsultat.
33
IV.5.c Affichage cran des donnes
A partir des donnes charges ltape prcdente, on les affiche lcran, avec la
visualisation voulue.
V ANALYSE DES RESULTATS
34
VI ETUDE THERMOMECANIQUE DU CONTACT ENTRE UNE
PLAQUETTE EN CARBONE ET UN DISQUE EN ACIER INOX
Cette tude est ralise a laide du logiciel ABAQUS. Cest un logiciel de calcul
lments finies dvelopp par ABAQUS, Inc.
Il se compose de trois produits : ABAQUS/Standard, ABAQUS/Explicit et
ABAQUS/CAE. ABAQUS/Standard est un solveur gnraliste qui recourt un schma
traditionnel d'intgration implicite. Le solveur ABAQUS/Explicit emploie un schma
d'intgration explicite pour rsoudre des problmes dynamiques ou quasi-statiques non-
linaires. ABAQUS/CAE constitue une interface intgre de visualisation et de
modlisation pour les dits solveurs. Chacun de ces produits est complt par des modules
additionnels et/ou optionnels, spcifiques certaines applications.
VI.1 Introduction
Lorsque deux corps entrent en contact avec frottement, il y a dissipation dnergie
et donc de la chaleur produite au niveau du contact, ce qui entrane une dilatation pouvant
accrotre le champ de pression, . Ce phnomne saccrot si les contraintes
tangentielles ainsi que les vitesses relatives de glissement entre les deux corps sont
importantes.
Ces effets thermomcaniques se traduisent la plus part du temps par la formation
de zones localises trs forts gradients thermiques, les points chauds. Il y a apparition
de dformations thermiques et des concentrations de contraintes pouvant gnrer des
fissures, des vibrations.
VI.2 Etude thermomcanique
Nous avons modlis le disque et la plaquette en caractrisant pour chaque pice
les proprits mcaniques et thermiques du matriau, prsentes dans le tableau.
Ensuite, nous avons appliqu les conditions aux limites du modle :
une pression applique sur la face suprieure de la paquette 1000 Pa.
Les nuds de la surface suprieur de la plaquette sont bloqus suivant x et y.
Les nuds du diamtre intrieur du disque sont bloqus suivant z et ont une
vitesse angulaire impose =157rad/min suivant z.
Sans oublier la cration des interactions du contact (frottement, ) entre les deux
pices.
35
Figure 6 : Schma du contact avec frottement entre les deux corps.
Disque plaquette
Matriau Acier inoxydable Carbone
Module de Young 200 000 Mpa 1 000 000 Mpa
Coefficient de poisson 0.3 0.3
Masse volumique 7900 Kg.m
-3
2200 Kg.m
-3
Conductivit 20 W.m
-1
.K
-1
65 W.m
-1
.K
-1
Capacit thermique 560 J.Kg
-1
.K
-1
1600 J.Kg
-1
.K
-1
Coefficient de frottement 0.9
Figure 7 : Tableau des caractristiques mcaniques et thermique.
VI.3 Rsultats
Le contact avec frottement plaquette/disque (lors de la phase du freinage), gnre
de la chaleur qui se diffuse sur la plaquette et le disque. La figure 4 reprsente la
Temprature produite par le frottement.
Pour pouvoir raliser une tude thermomcanique, le type danalyse est une analyse
Dynamique, Temp-disp, Explicite. Le maillage cre est un maillage coupl temprature-
dplacement figure 3. En effet Il est ncessaire de connatre la temprature des nuds
afin de pouvoir prendre en compte la dilatation du maillage mcanique.
36
Figure 8 : couplage lment mcanique-thermique
37
dcollement
glissement
adhrence
Figure 9 : Temprature mise par la surface de contact plaquette/disque.
Le rsultats obtenus, nous montrent la localisation de lapparition de la chaleur ainsi
que sa diffusion au cours du temps lors de la phase de freinage.
Nous remarquons aussi trois zones diffrentes sur la surface de contact. Car lors
dun contact avec frottement, il apparat une surface de glissement, de dcollement et
dadhrence.
38
Conclusion :
Ce sujet dtude a tout dabord servi nous familiariser avec les logiciels STAR-CD et
Abaqus. En suivant les tutoriaux divers, nous avons remarqu que ce sujet inclut
lutilisation de mthodologies prsentes dans la plupart des cas danalyses simples .
En effet, nous avons du associer lcoulement dair, des paramtres de thermique
ainsi que lutilisation de paramtres que lon retrouve dans dautre cas, comme par
exemple la rotation du disque qui a t ralise comme dans les cas dcoulement autour
dun cylindre.
blablabla
.
..
..
..
.
.
.
.
.
Cette tude peut tre bien mieux optimise, en utilisant un maillage spar pour
chaque pice. Il sagit de la mthode utilise pour les cas danalyse de ventilateurs. Les
pices mailles peuvent tre assembles, et en utilisant la C.L. Attach STAR-CD
diffrencie mieux les lments solides en mouvements des solides fixes et du fluide.
Cependant cette mthode requiert des connaissances spcifiques du logiciel et est trs
difficile mettre en place, voici pourquoi nous sommes rests sur notre modle.
39
VII Bibliographie
Transient temperature field analysis of a brake in a non-axisymmetric three-dimensional
model C. H. Gao and X. Z. Lin, Fuzhuo University China, 2002.
StarCD, Pro_Star tutorial, chap. X
StarCD, Method & User Guide, V 3.10
Solutions of 3D Navier-Stockes benchmark problems with adaptive finite elements, M.
Braack and T. Richter, To appear in Computers & Fluids, 2005
Numerical anjd experimental analysis of a pegs-wing ventilated disk brake rotor, wiht pads
and cylinders, P. F. Gotowicki, V. Nigrelli, G. Virzi Mariotti, 2005
Quasi-3D flow simulation of a real disc cavity with conjugate heat transfer, Glen C.
Snedden, 2005
Exploration du code Star-Cd, P. Claudon
A brake disc's temperature estimation module, S. Artus, 2004
Aspect nergtique du contact : couplage thermomcanique, M. Bettaieb, thse
Modlisation du comportement mcanique dengrenages, couplages entre engrenement
et lments de structure flexibles , 2006
Simulation du comportement dun frein disque, LENSIETA, Partenaire du secteur
automobile
Internet
Site dtudiant trs utile sur ltude dune olienne avec les dmarche suivre et les
conseils
http://www.enseeih t.fr/hmf/travaux/CD0203/travaux/optsee/bei/1/conseil.html
Site ddi la moto et particulirement sur la technique adopte par la Mtisse avec son
mono-bras avant :
http://technomoto.blogspot.com/2006/04/metiss.html
Gnralits sur les freins :
http://www.techniques-
ingenieur.fr/dossier/technologie_du_freinage_freins_a_disque/B5574
Application sur StarCD relative aux freins :
http://www.cd-adapco.com/applications/automotive.html#
40
41