RUGC17 - Article2 - Libasse SOW
RUGC17 - Article2 - Libasse SOW
RUGC17 - Article2 - Libasse SOW
RÉSUMÉ. Partant d’une granulométrie de particules de MIOM 0/25 qui a fait l’objet d’une étude expérimentale dans
[BEC 07], nous avons mis en place un modèle numérique 3D dans des volumes élémentaires représentatifs. Ce modèle est
basé sur une stratégie de modélisation dite « multi-échelles hiérarchique ». Les agrégats traités au bitume ont été pris en
compte explicitement dans deux volumes élémentaires représentatifs : sub-mesoscopique (0/6) et mesoscopique (6/25). Les
tests de relaxation de contraintes effectués aux deux échelles ont permis d’obtenir les propriétés thermo-rhéologiques de la
grave de MIOM traités à 3 % de bitume. Ces propriétés sont les poids de la série de Prony à une température de référence de
0 °C ainsi que les constantes de la loi WLF pour décrire le comportement thermo-rhéologiquement simple. Les séries de
Prony sont obtenues par l’intermédiaire des modules élastiques normalisés. Les constantes de loi WLF obtenues sont
C1=20 °C-1 et C2=130 °C. Une validation numérique a été faite en comparant les réponses du matériau hétérogène de grave-
bitume et celles d’un matériau homogène à qui nous avons donné les propriétés tirées du matériau hétérogène.
ABSTRACT. Starting with a particle size of particles MSWI bottom ash 0/25 which was the subject of an experimental study
in the thesis of Frédéric Becquart [BEC 07], a 3D numerical model has been developed. This model is based on a
hierarchical multi-scale modelling strategy. In this study, the aggregates are embedded in bitumen. Two representative
elementary volumes are considered: sub-mesoscopic (0/6) and mesoscopic 6/25). The stress relaxation tests carried out at
the two scales enable one to obtain the thermo-rheological properties of the micro-aggregate treated with 3% bitumen. These
properties are the Prony series at a reference temperature of 0 °C and the WLF law constants. Prony series are obtained by
means of standardized elastic moduli. The obtained WLF law constants are C1=20 °C-1 and C2=130 °C. The numerical
simulations have been validated by comparing the responses of heterogeneous and corresponding homogeneous materials.
MOTS-CLÉS : modélisation multi-échelles hiérarchique, grave routière, MIOM, bitume, viscoélastique, série de Prony, loi
WLF, modules élastiques normalisés
KEY WORDS: hierarchical multi-scale modelling, road aggregates, MSWI bottom ash, bitumen, viscoelastic, Prony series,
WLF law, standardized elastic moduli
1. Introduction
En techniques routières, les différentes couches doivent être constituées de matériaux capables de reprendre et
de repartir les charges provenant des véhicules roulants. Les assises de chaussée jouent un rôle structurel
important du fait qu’elles réduisent les contraintes exercées sur le sol de fondation. Les matériaux constitutifs de
ces couches doivent présenter des caractéristiques mécaniques leur permettant de remplir pleinement leurs
fonctions dans la structure de chaussée. Face à une expansion démographique importante et une consommation
sans cesse croissante des populations, les autorisés publiques, notamment les villes, se trouvent confrontées à des
problèmes d’urbanisation aigus les poussant à revoir leurs politiques de gestion des déchets. Dans un cadre de
recyclage et de développement durable, les Mâchefers d’Incinération d’Ordures Ménagères (MIOM) peuvent
représenter un gisement de substitution de granulats classiques stratégiques pour faire face à l’impossibilité de
mise en décharge continue des déchets et aux difficultés d’ouverture de nouvelles carrières de matériaux routiers.
Plusieurs auteurs se sont intéressés à la modélisation des matériaux bitumineux : [LYT 90], [BUT 99],
[MET 99], [SHA 02], [KIM 04], [SHU 08], [BAN 08], [NGU 08], [DAI 11] et [TEH 12]. Dans sa thèse,
[NGU 08] a fait une étude numérique de la fissuration d’un milieu viscoélastique avec une analyse de l’essai de
rupture sur un bitume viscoélastique pur de classe 50/70. Il a déterminé par un travail expérimental les propriétés
thermo-rhéologiques (paramètres de la série de Prony, les constantes de la loi WLF, la courbe maitresse, entre
autres). [TEH 12] a utilisé les propriétés du bitume de [NGU 08] et a mis en place une stratégie de modélisation
numérique basée sur la méthode des éléments finis pour déterminer le module dynamique des enrobés poreux
(matrice/inclusions). Le matériau considéré est modélisé comme un milieu bi-phasique 2D composé d'inclusions
granulaires avec des propriétés élastiques linéaires et une matrice poreuse bitumineuse contenant un pourcentage
donné de porosité (6 à 5%) présentant un comportement viscoélastique linéaire pour des valeurs de déformation
faibles (Ɛmax ≤ 10-4). Dans ce présent article, nous développons une modélisation 3D de graves routières à base
de MIOM traitées au bitume. Cette modélisation prend en compte explicitement les différentes phases du
matériau granulaire traité (liant, squelette granulaire, porosité) dans des volumes élémentaires représentatifs. Le
liant bitumineux est celui de [NGU 08] et de [TEH 12]. Les MIOM en question (figure 1) sont ceux étudiés
35èmes Rencontres de l’AUGC, ECN/UN, Nantes, 22 au 24 mai 2017 2
Comme tout matériau présentant un squelette granulaire, les MIOM traités peuvent voir leur structure
granulométrique se diviser en différentes échelles. Nous développons une modélisation multi-échelles
hiérarchique consistant à représenter l’hétérogénéité de la grave routière traitée au bitume dans des volumes
élémentaires représentatifs des différentes échelles d’études. L’objectif est de mettre en place une méthodologie
de détermination par simulations numériques des paramètres thermo-rhéologiques simples de la grave routière à
base de MIOM traités à 3% de bitume.
Le comportement mécanique des phases viscoélastiques à chaque échelle est pris en compte en développant
le module de cisaillement et le module volumique sous la forme de séries de Prony :
N
p
GR (t) = G0 (1 − ∑ g̅ k (1 − e−t⁄τk )) [1]
k=1
N
p
K R (t) = K 0 (1 − ∑ k̅ k (1 − e−t⁄τk )) [2]
k=1
Ces expressions montrent que le comportement du matériau peut être phénoménologiquement reproduit par
des associations en parallèle de ressorts et d’amortissements en séries. Les résultats de [NGU 08] sont utilisés :
13 ressorts et amortisseurs sont nécessaires pour bien reproduire le comportement viscoélastique du bitume. En
fonction de la température, le bitume présente un comportement thermo-rhéologiquement simple : il suffit ainsi
de connaître son comportement à une température de référence T ref pour le déduire à n’importe quelle
température. Ce comportement est pris en compte par l’intermédiaire du principe d’équivalence temps-
température. Ce principe permet de préciser la dépendance en température des grandeurs viscoélastiques ; il
stipule que le passage de la température de référence Tref à une température T revient à multiplier l’échelle du
temps par un facteur constant (dit facteur de décalage) aT. La loi WLF est utilisée pour exprimer aT :
−C1(T−Tref )
log aT = C2 +T−Tref
[3]
Nous utilisons les résultats de [NGU 08] à la température de référence Tref = 0 °C. Identifier le comportement
thermo-rhéologiquement simple d’un matériau revient à identifier C1 et C2. La stratégie de simulation adoptée est
celle qui suit :
35èmes Rencontres de l’AUGC, ECN/UN, Nantes, 22 au 24 mai 2017 3
Tableau 4. Ecarts relatifs moyens pour les différentes températures d’essai entre réponse VER homogène de
MGB et VER hétérogène de MGB
Tableau 5. Ecarts relatifs moyens pour les différentes températures d’essai entre réponse VER homogène de
grave-bitume et VER hétérogène de grave-bitume
6. Conclusions et perspectives
En partant d’une granulométrie de MIOM, nous avons mis en place un modèle numérique de modélisation
multi-échelles hiérarchique représentant la grave de MIOM traitée à 3% de bitume. La formulation sur laquelle
nous sommes partis correspond à la courbe granulométrique moyenne du fuseau retenu pour les graves routières.
Avec ce modèle numérique 3D, une stratégie de changement d’échelles avec une méthode de détermination
des paramètres thermo-rhéologiques a pu être testée et validée. Les propriétés thermo-rhéologiques de la micro-
grave-bitume ainsi que celles de la grave de MIOM traitée à 3 % de bitume ont ainsi été déterminées. Nous
sommes parvenus à déterminer : les séries de Prony par l’intermédiaire de modules élastiques normalisés et les
constantes de la loi WLF par le biais de la relaxation de contraintes en fonction du temps.
Ce modèle numérique 3D ouvre des perspectives d’application pour étudier les problèmes courants en
techniques routières tels que l’orniérage, la fissuration, ainsi tous les autres types de dégradations de chaussée.
REFERENCES BIBLIOGRAPHIQUES :
[APP 04] Rapport APPIA. Valorisation des MIOM dans les infrastructures routières. Ecole des mines de Doua,
Laboratoire régional des ponts et chaussées de Lille. Avril 2004.
[BAN 08] R. Bandyopadhyaya, A. Das, S. Basu, Constr. Build. Mater. 22 (2008) 1051– 1058.
[BEC 07] Frédéric Becquart. Première approche du comportement mécanique d’un milieu granulaire issu d’un mâchefer
d’incinération d’ordures ménagères : valorisation en technique routière. PhD thesis. Université des sciences et technologies
de Lille. Mars 2007.
[BUT 99] W.G. Buttlar, D. Bouzkurt, G. Al-khateeb, A.S. Waldhoff, Trans. Res. Rec. 1681 (1999) 157–169.
[DAI 11] Q. Dai, Constr. Build. Mater. 25 (2011) 1102–1114.
[KIM 08] Y. Kim, D.N. Little, J. Mater. Civil Eng. 16 (2) (2004).
[LYT 90] R. Lytton. Strategy High. Resolution Progress, Washington, DC, 1990.
[MET 99] G. Li, Y. Li, J.B. Metcalf, S.S. Pang, J. Mater. Civil Eng. 11 (3) (1999) 236–241.
[NGU 08] Hoai Nam Nguyen. Etude numérique de la fissuration d’un milieu viscoélastique : Analyse de l’essai de rupture
sur bitume. PhD thesis. Ecole Nationale des Ponts et Chaussée. Décembre 2008.
[SHA 02] N. Shashidhar, A. Shenoy, Mech. Mater. 34 (2002) 657–669.
[SHU 08] X. Shu, B. Huang, J. Mater. Civil. Eng. 20 (8) (2008) 530–538.
[TEH 12] F. Fakhari Tehrani, J. Absi, F. Allou, Ch. Petit. Heterogeneous numerical modeling of asphalt concrete through
use of a biphasic approach: Porous matrix/inclusions. Computational Materials Science. Juin 2012.