SADOK
SADOK
SADOK
THESE
Présentée en vue
d’obtenir le grade de
DOCTEUR
en
Par
Abdelfeteh SADOK
DOCTORAT DE l’UNIVERSITE DE LILLE 1
Titre de la thèse :
Je remercie bien évidemment mes parents d’avoir fait de moi ce que je suis aujourd’hui et
pour leur soutien malgré la distance. Un grand Merci à mes frères Mustapha, Djamel et ma
sœur Imene pour leurs soutiens et leurs encouragements. Je sais que vous êtes fiers de moi et
je ferai en sorte que vous le soyez toujours.
Un merci, aussi distingué que la concernée, s’adresse à ma future épouse Afaf, pour son
soutien indéfectible durant ces années. Merci ma chérie pour ton enthousiasme, ta patience et
tes encouragements.
J’exprime ma reconnaissance envers toutes les personnes qui ont, de près ou de loin,
contribué à la réalisation de ce travail.
RESUME
La valorisation des matériaux alternatifs dans le domaine du génie civil connaît un essor
important ces dernières années, toutefois cette pratique connait certaines limites, notamment
le manque de réglementation claire et surtout le manque d’outils et méthodes de formulation
adaptées et adaptables à l’utilisation des matériaux alternatifs.
Dans ce contexte, le présent travail porte sur la mise en place de méthode de formulation
de matériaux de construction à base de sous-produits industriels. Cette méthode hybride basée
sur des outils de l’intelligence artificielle, combine les Algorithmes Génétique (AG) comme
outils d’optimisation multiobjectifs et la Programmation Génétique (PG) dans ses deux
versions, la PG classique et la MGGP (MultiGène Genetic Programming) comme outils de
modélisation des problèmes complexes par Apprentissage Automatique (Machine Learning).
Des études spécifiques ont été menées sur le fonctionnement de ces outils novateurs, pour
appréhender les avantages et faiblesses de ces outils sur ces cas d’applications dans le Génie
Civil.
The predictions of scientist are increasingly alarming, and the environmental issue has
become a major concern for governments and industry. Effective waste management is part of
the priority actions in order to achieve a green circular economy. This efficient management
imposes first maximum recovery of waste, knowing the large tonnage produced is different
sectors. The field of civil engineering is particularly concerned as it occupies the second place
among the waste sectors producer and largest consumer of raw materials.
The valorization of alternative materials in the field of civil engineering has grown
significantly in recent years. However, this practice knows some limitations, including the
lack of clear regulations and especially the lack of tools and methods suitable for design of
materials including alternatives materials.
In this context, the present work focuses on the development of mix design method of
building materials based on industrial by-products. This hybrid method combines the Genetic
Algorithms (GA) as multi-objective optimization tools and Genetic Programming (GP) in its
two versions, classical GP and MGGP (MultiGene Genetic Programming) as modeling tools
for complex problems by Machine Learning approach.
Specific studies were carried out also or these innovative tools, to demonstrates benefits
and weaknesses of these tools on these applications in civil engineering.
Finally, this method of formulation of building materials based on industrial sub products
proposed in this work was tested on two case studies (design of high performance concrete
and mortars made of alternative materials) and validated by the laboratory tests. The results
are conclusive and promising to generalize the method to other applications of Civil
Engineering.
REMERCIEMENTS ............................................................................................................ 5
RESUME .............................................................................................................................. 6
ABSTRACT ......................................................................................................................... 7
TABLE DES MATIERES ................................................................................................... 8
Liste des symboles et notations .......................................................................................... 12
Liste des figures.................................................................................................................. 15
Liste des tableaux ............................................................................................................... 17
Introduction générale et problématique .............................................................................. 19
CHAPITRE 1:Etat de l’art et contexte ............................................................................... 23
Introduction du chapitre ..................................................................................................... 24
Partie I. La valorisation des matériaux alternatifs dans les bétons ............................... 25
1 Les enjeux du développement durable dans le Génie Civil....................................... 25
2 De l’Eco-matériau à l’économie circulaire ................................................................ 27
3 Les matériaux alternatifs dans le Génie Civil ............................................................ 28
3.1 Les matériaux de démolition .................................................................................. 31
3.1.1 Les catégories de matériaux de démolition .................................................... 32
3.1.2 Paramètres influents sur les caractéristiques des granulats recyclés .............. 33
3.1.3 Avancées scientifiques et techniques en Europe sur les granulats recyclés ... 34
3.2 Boues et sédiments de dragage .............................................................................. 35
3.2.1 Principaux paramètres influents sur les caractéristiques des sédiments ......... 36
3.2.2 Travaux scientifiques et projets relatifs à la valorisation des sédiments ........ 37
3.3 Résidus de procédés thermiques ............................................................................ 37
3.3.1 Les cendres volantes de centrale thermique (CVCT) ..................................... 38
3.3.2 Les mâchefers d’incinération des déchets non dangereux (MIDND) ............ 39
3.3.3 Les sables de Fonderie .................................................................................... 40
3.4 Les particularités des matériaux alternatifs et les contraintes de la valorisation ... 42
3.5 Evaluation environnementale des matériaux alternatifs ........................................ 43
3.6 Principales voies de valorisation des matériaux alternatifs ................................... 43
3.6.1 Dans le matériau béton ................................................................................... 44
3.6.2 Dans le matériau routier ................................................................................. 45
3.6.3 Dans d’autres applications .............................................................................. 46
4 Evolution des méthodes de formulation des bétons .................................................. 46
Partie II. L’intelligence artificielle dans le Génie Civil ................................................. 49
1 Les problèmes d’optimisation ................................................................................... 49
1.1 Vocabulaire et définitions ...................................................................................... 50
1.2 Les Problèmes d’Optimisation Multiobjectifs (POM) ........................................... 51
1.3 La surface de compromis (Frontière de Pareto) ..................................................... 52
2 Les Algorithmes Génétiques ..................................................................................... 53
2.1 Evolution des Algorithmes Génétiques.................................................................. 54
2.2 Avantages et limites des Algorithmes Génétiques ................................................. 55
3 L’apprentissage automatique (Principe et applications) ............................................ 56
3.1 Types d'apprentissage ............................................................................................ 56
3.2 Facteurs de pertinence de l’apprentissage .............................................................. 56
3.3 La Programmation Génétique : du fonctionnement à la modélisation................... 57
3.3.1 Evolution de la Programmation Génétique .................................................... 59
3.3.2 Multi-Gène Genetique Programming (MGGP). ............................................. 59
4 Démarche globale ...................................................................................................... 60
5 Analyse SWOT des Algorithmes Génétiques et la Programmation Génétique ........ 61
Partie III. Contours de la thèse .................................................................................... 62
1 Objectifs et démarches ............................................................................................... 62
2 Les limites de la thèse ................................................................................................ 64
Conclusion du chapitre ....................................................................................................... 65
CHAPITRE 2: Caractérisation des matériaux .................................................................... 67
Introduction du chapitre ..................................................................................................... 68
1 Matériel et mesures .................................................................................................... 69
1.1 La granulométrie .................................................................................................... 69
1.2 Les masses volumiques .......................................................................................... 70
1.3 La compacité .......................................................................................................... 71
1.4 Pertes Au Feu (P.A.F) ............................................................................................ 71
1.5 Surface spécifique Blaine et BET .......................................................................... 71
1.6 Analyse chimique par Fluorescence aux rayons X (FX) ....................................... 72
1.7 Diffraction aux Rayons X (DRX) .......................................................................... 73
2 Caractérisation des matériaux standards.................................................................... 74
2.1 Les granulats .......................................................................................................... 74
2.2 Le ciment (Choix et origine) .................................................................................. 76
2.2.1 Fiche d’identité du ciment .............................................................................. 76
2.2.2 Analyse granulométrique ................................................................................ 77
2.2.3 Composition minéralogique et chimique ........................................................ 77
2.3 Les additions minérales (Les fines) ....................................................................... 79
2.3.1 Les cendres volantes ....................................................................................... 79
2.3.2 Fillers calcaires du boulonnais ....................................................................... 80
2.4 Les adjuvants ......................................................................................................... 81
2.4.1 Superplastifiant ............................................................................................... 81
2.4.2 Agent entraîneur d’air ..................................................................................... 82
3 Caractérisation des matériaux alternatifs ................................................................... 82
3.1 Les Graviers de Béton Recyclé (GBR) .................................................................. 82
3.2 Les sables de fonderie bruts ................................................................................... 85
3.3 Sables de dragage ................................................................................................... 86
3.4 Sédiments de dragage (Marins et fluviaux) ........................................................... 88
3.4.1 Sédiments de Dunkerque ................................................................................ 88
3.4.2 Sédiments de la Deûle .................................................................................... 90
Conclusion du chapitre ....................................................................................................... 91
CHAPITRE 3 : Modélisation de la compacité granulaire .................................................. 93
Introduction du chapitre ..................................................................................................... 94
1 Base de données relative à la compacité des mélanges granulaires binaires ............. 96
2 Etude paramétrique de la compacité granulaire......................................................... 96
2.1 Matériaux de l’étude .............................................................................................. 96
2.2 Modes de compactage ............................................................................................ 97
2.3 Influence du mode de compactage sur la compacité des classes élémentaires .... 101
2.4 Compacité des mélanges granulaires binaires ..................................................... 102
2.5 Effet de l’étendue granulaire sur la compacité des mélanges binaires ................ 108
3 Apprentissage automatique par Programmation Génétique .................................... 110
3.1 Critère d’évaluation des performances du modèle ............................................... 110
3.2 Paramétrage de la PG ........................................................................................... 112
3.3 Modélisation de la compacité granulaire binaire par PG ..................................... 113
3.4 Performances du modèle développé .................................................................... 114
3.5 Modélisation de la compacité granulaire binaire par MGGP .............................. 118
4 L’interface GUImines .............................................................................................. 119
5 Validation de l’approche : L’optimisation du squelette granulaire pour la formulation
de bétons ................................................................................................................................. 120
Conclusion du chapitre ..................................................................................................... 123
CHAPITRE 4 : Apprentissage automatique et optimisation multiobjectifs..................... 125
Introduction du chapitre ................................................................................................... 126
1 Application de la méthode hybride AG&PG dans la formulation de BHP ............. 127
2 Base de données d’apprentissage ............................................................................ 127
3 Modélisation de la résistance mécanique d’un BHP ............................................... 130
3.1 Modélisation de la résistance mécanique d’un BHP par PG ............................... 130
3.2 Modélisation de la résistance mécanique d’un BHP par MGGP ......................... 133
3.3 Performances des modèles développés de la résistance mécanique .................... 134
4 Modélisation de l’affaissement d’un BHP ............................................................... 136
4.1 Modélisation de l’affaissement d’un BHP par PG ............................................... 136
4.2 Modélisation de l’affaissement d’un BHP par MGGP ........................................ 136
4.3 Performances des modèles de l’affaissement développés.................................... 137
5 Formulation d’un BHP par les Algorithmes Génétiques (AG) ............................... 139
5.1 Paramétrage de l’algorithme génétique................................................................ 139
5.2 Validation de l’approche par un cas pratique....................................................... 140
6 Impact de l’écriture des modèles prédictifs sur l’optimisation multiobjectifs ........ 144
Conclusion du chapitre ..................................................................................................... 148
CHAPITRE 5 : Formulation de mortiers à base de matériaux alternatifs ........................ 149
Introduction du chapitre ................................................................................................... 150
1 Plan d’expériences et développement de la base de données du mortier ................ 151
2 Modélisation des propriétés du mortier à base de sable fin de fonderie .................. 156
2.1 Modélisation de la résistance à la compression du mortier par MGGP ............... 156
2.2 Modélisation de la résistance à la flexion du mortier par MGGP ........................ 158
2.3 Modélisation de la maniabilité du mortier par MGGP ........................................ 160
2.4 Modèle de la compacité granulaire ...................................................................... 161
2.5 Modèle du coût .................................................................................................... 162
3 Optimisation multiobjectifs des mélanges de mortier à base de sable de fonderie . 163
3.1 Résultats et discussion ......................................................................................... 164
3.2 Validation au laboratoire des formulations de mortiers ....................................... 165
4 Application de la méthode pour mortiers à base de sable de dragage ..................... 167
4.1 Résultats et discussion ......................................................................................... 167
4.2 Validation au laboratoire des formulations de mortiers ....................................... 168
Conclusion du chapitre ..................................................................................................... 170
Conclusion générale et perspectives................................................................................. 171
Références bibliographiques ............................................................................................ 176
Annexes ............................................................................................................................ 182
Publications dans le même domaine ................................................................................ 185
Liste des symboles et notations
Notations
La valorisation des matériaux alternatifs dans le domaine du génie civil connait un essor
important ces dernières années, toutefois cette pratique connait certaines limites. En effet, en
plus d’un manque de réglementation claire, le manque d’outils adaptés et adaptables à
l’utilisation des matériaux alternatifs, font reculer parfois les industriels, du fait que les
matériaux alternatifs ne présentent pas les mêmes caractéristiques que les matériaux
standards. Aussi, leurs utilisations avec des matériaux standards nécessitent de mettre en place
des méthodes de formulation prenant en compte leurs spécificités.
Cependant, l’efficacité des AG dépend fortement des fonctions d’évaluation des individus
ou des solutions proposées, qui sont des modèles utilisés en guise de fonctions-objectifs
(Park, 2013). Ainsi le choix des fonctions-objectifs sera d’une importance capitale. Il est, au
préalable, primordial de développer ou d’optimiser les modèles prédictifs existants avec des
outils adaptés. Nous avons pour cela utilisé l’Apprentissage Automatique via la
Programmation Génétique « PG » comme algorithme d’apprentissage performant, pour
développer des modèles prédictifs adaptés à l’utilisation des matériaux alternatifs.
La combinaison des AG et la PG, nous a permis de proposer une méthode hybride qui
répond au problème de formulation des matériaux alternatifs pour les applications du Génie
Civil. Cette méthode a été appliquée et validée pour deux cas d’étude : le premier cas, un
béton à hautes performances à base de cendre volantes et le second cas, un mortier à base de
sable de fonderie et sable de dragage.
Des études spécifiques ont été menées également, sur ces outils novateurs, pour
appréhender la répétabilité des résultats de l’apprentissage automatique sur les mêmes
données d’apprentissage et l’impact de la forme de l’écriture mathématique des modèles
développés sur les résultats de l’optimisation par AG. Nous nous sommes intéressés
également à la robustesse des modèles développés par la PG ou la MGGP.
Introduction du chapitre
Dans ce chapitre, nous allons faire un état de l’art sur la production des déchets, la
règlementation actuelle qui régit la gestion de ces derniers, ainsi que les avancées
scientifiques autour de la valorisation des déchets dans le domaine du Génie Civil, ensuite
nous définissons les concepts de base de l’intelligence artificielle en rapport avec les outils
utilisés dans cette thèse, dans le but de proposer une méthode novatrice de formulation de
béton à base de matériaux alternatifs.
24
Chapitre 1 : Etat de l’art et contexte
Dans le contexte actuel, la capacité des générations futures à subvenir à leurs besoins se
trouve incontestablement menacée. Des situations critiques dont les indicateurs les plus
parlants sont le réchauffement climatique identifié, aujourd’hui, comme une conséquence du
rejet en quantité importante de gaz à effet de serre par les activités humaines et la raréfaction
des ressources naturelles, dont certaines sont indispensables au fonctionnement de nos
sociétés. Devant ce schéma, le problème du développement durable ne se pose donc plus sous
l'angle des droits des générations futurs, mais des obligations et des devoirs des générations
actuelles(Wikipédia, 2016).
Dans le cadre des engagements au niveau mondial et européen, les différents pays agissent
prioritairement sur les secteurs économiques de grands impacts sur les consommations
énergétiques ou l’émission de GES, à savoir le transport et l’industrie et surtout le secteur du
bâtiment. La Directive Européenne sur la Performance Energétique des Bâtiments (DPEB) de
2002, constitue le cadre règlementaire au niveau européen qui fixe les exigences concernant le
secteur du bâtiment pour atteindre les engagements pris dans le cadre du protocole de Kyoto.
25
Chapitre 1 : Etat de l’art et contexte
En effet, même si actuellement, l’exploitation des bâtiments prend la plus grande part des
consommations énergétiques, l’évolution des constructions de plus en plus isolées va tendre à
augmenter davantage la part des matériaux de construction dans les consommations
énergétiques et les émissions de CO2. La Figure 1 présente l’évolution du poids relatif des
matériaux de construction et du chauffage dans les consommations énergétiques des
bâtiments. Il passe de 15% à 75% dans des bâtiments consommant respectivement 200 et 15
kWh/m2.an.
26
Chapitre 1 : Etat de l’art et contexte
Le nouveau concept d’une économie qui fonctionne en boucle nous a imposé, en premier
lieu, de nous passer de la notion de déchet et la notion de l’Eco matériau est désormais
privilégiée.
A ce jour, il n'y a pas encore de définition officielle de l'éco matériau. En tenant compte
de la définition d’un matériau alternatif, on peut le définir comme « un matériau renouvelé,
issu d’un usage antérieur, il satisfait naturellement ou après traitement aux cahiers des
charges des matériaux naturels standards, il est de faible impact environnemental et
potentiellement intégrable dans la boucle d’économie circulaire».
27
Chapitre 1 : Etat de l’art et contexte
L’économie circulaire vise à changer la vision du monde par rapport à l’économie dite
linéaire, en limitant le gaspillage des ressources et l’impact environnemental, tout en
augmentant l’efficacité des produits à tous les stades de vie. On peut scinder l’économie
circulaire en deux volets, le premier est la Production et l’offre de biens et de services, qui
porte principalement sur une conception intelligente et fonctionnelle, quant au deuxième volet
il vise la Consommation, la demande et le comportement qui porte sur une consommation
responsable, collaborative et privilégie une consommation en boucle fermée. La Figure 2
récapitule les composantes principales de l’économie circulaire.
La loi de Transition énergétique pour la croissance verte adoptée en Aout 2015 (LOI n°
2015-992 du 17 août 2015 relative à la transition énergétique pour la croissance verte,
2015), se donne comme objectif de « lutter contre les gaspillages et promouvoir l’économie
circulaire ». Dans ce cadre, la politique nationale de prévention et de gestion des déchets
apparaît comme un levier essentiel de cette transition vers l’économie circulaire.
28
Chapitre 1 : Etat de l’art et contexte
29
Chapitre 1 : Etat de l’art et contexte
La production de déchets intervient aussi bien dans les étapes de fabrication de biens
(déchets de fabrication) que dans les étapes de consommation intermédiaires et finales
(déchets post-consommation) d’où leurs différences de nature. La Figure 5 présente la
composition des déchets du BTP produits en France.
Déchets
dangereux
2%
Déchets
industriels
banals
26%
Déchets inertes
72%
Figure 5 : Les types des déchets produits par l’activité de la construction en France (FFB,
2014).
30
Chapitre 1 : Etat de l’art et contexte
Bien que le recyclage des chantiers routiers s’avère très efficace en France, la réutilisation
des déchets de démolition du secteur du bâtiment reste encore limitée. En effet, les déchets
des ouvrages d’art sont facilement recyclables, puisque ils sont souvent constitués uniquement
de béton, alors que dans un bâtiment, le béton cohabite avec d’autres matériaux très divers
(Plâtre, bois, plastiques, métaux, verre …). Le recyclage du béton n’est alors possible qu’à la
suite d’un tri des matériaux, souvent onéreux, on parle ici de déconstruction.
Les utilisations les plus courantes, actuellement, des bétons recyclés sont dans les
remblaiements et les sous-couches de fondation en technique routière. Or, la réutilisation de
ces bétons sur place serait judicieuse d’un point de vue technique et économique, puisque cela
permettra à la fois une économie des ressources naturelles, une limitation de la mise en
décharge, mais surtout une mise en œuvre rapide et une diminution des transports.
31
Chapitre 1 : Etat de l’art et contexte
cette utilisation est limitée dans les proportions qui dépendent étroitement de la qualité des
granulats recyclés et de la classe de résistance souhaitée pour le béton (Laneyrie, 2014).
D’après les chiffres de l’UNPG (Union National des Producteurs de Granulats), en 2012
l’utilisation annuelle de granulats en France a atteint 360 millions de tonnes. Cependant les
importantes sources d’approvisionnement en granulats naturels en France, en comparaison
avec les autres pays européens, représentent un frein pour le développement de la filière si on
néglige le coût environnemental, ce qui est malheureusement le cas actuellement. Le Tableau
1 récapitule les prix des granulats standards et recyclés dans les pays européens. Le rapport ½
exprime le rapport de prix entre agrégat recyclés et naturels.
Pays Prix moyen des agrégats Prix moyen des Rapport (1)/(2)
recyclés agrégats naturels
(1), €/t (2), €//t
France 4,4-9,0 3,6-6,5 1,2-1,4
Belgique 2,4-7,2 2,4-10,8 0,65-1
Allemagne 4,8-8,4 5,6-8,8 0,85-0,95
Pays-Bas 8,7 11,4 0,75
Royaume-Uni 3,3-8,4 4,4-15,2 0,55-0,75
Espagne 3,0 6,0-14,4 0,2-0,5
Tableau 1 : Prix comparés des granulats recyclés et de granulats naturels dans différents
pays de l'Union Européenne (UEPG, 2005).
32
Chapitre 1 : Etat de l’art et contexte
Les catégories (1) et (2) couvrent plus de 90% des approvisionnements. Il existe
cependant quelques étapes primordiales à suivre pour obtenir un produit de recyclage
répondant aux caractéristiques souhaitées à savoir :
Réduction éventuelle des plus gros éléments à l’aide d’un brise roche hydraulique,
Concassage primaire à l’aide d’un concasseur à percussion ou à mâchoires, suivi
d’un déferraillage électromagnétique,
Criblage, visant à éliminer les matériaux indésirables,
Si nécessaire, concassage secondaire sur la fraction supérieure du concassage
primaire.
Quant à la qualité des granulats recyclés, elle se trouve fortement influencée par la
présence de contaminants des sites de leurs provenances. En effet, on trouve parmi les
granulats recyclés, du bois, du plastique, du verre, du bitume, du gypse, des matières
organiques de différentes natures ou des contaminations sous forme de chlorures ou de
33
Chapitre 1 : Etat de l’art et contexte
Les engagements européens sur la question environnementale ont permis la mise en place
de projets expérimentaux de réutilisation du béton recyclé dans plusieurs pays, dont fait partie
le projet français RECYBETON, dont le but est d’aboutir à des normes visant à cadrer
l’utilisation des bétons recyclés, afin de permettre enfin en France leur bonne utilisation dans
le bâtiment (IREX, 2011).
Le projet national sur le recyclage complet des bétons (RECYBETON) est un projet de
recherche et développement collaboratif soutenu par le Ministère de l’Ecologie et du
Développement Durable, d’une durée de 5 ans. Il regroupe des industriels et universitaires
spécialistes du génie civil, et vise à changer la tendance au niveau national en favorisant la
réutilisation de l’intégralité des produits issus des bétons déconstruits. Le projet comporte
d’autres travaux, nous citons :
34
Chapitre 1 : Etat de l’art et contexte
Les boues de curage sont définies comme étant les matières (exceptées les matières
exogènes) enlevées du lit et des berges des cours et plans d’eau ou de leurs ouvrages annexes
du fait des travaux de dragage ou de curage. Les sédiments de dragage marins préoccupent
particulièrement les gestionnaires en raison des quantités énormes extraites chaque année
(Mekerta et al., 2016). En effet, leur valorisation permet de réduire le coût des opérations de
dragage et peuvent combler le manque de matière première (Bourabah, 2012), (Sadok,
2013a).
Les sédiments sont des dépôts qui proviennent de l’altération des roches préexistantes.
Ces dépôts sont ensuite transportés par les cours d’eau et/ou par le vent, pour finalement se
déposer dans des bassins de sédimentation et former des roches sédimentaires. Ces sédiments
appelés aussi vases sont la désignation d’une large famille de matériaux fins, argileux, plus ou
moins organiques, pouvant atteindre des teneurs en eau importantes jusqu’ à 300 % (Dubois,
2006). En général, en plus de l’eau, la constitution des sédiments est la suivante (Khezami,
2014):
L'enquête dragage réalisée par le CETMEF (Centre d’Etudes Maritimes et Fluviales) fait
état d'un volume de 50 millions de m3 de sédiments marins dragués en 2009 en France,
départements d'Outre-Mer compris (Le Guyader and Colin, 2012). Environ 81% de ces
sédiments sont dragués par les 7 Grands Ports Maritimes, qui sont le Havre, Bordeaux,
Rouen, Nantes, Marseille, Dunkerque et La Rochelle (Tableau 3).
35
Chapitre 1 : Etat de l’art et contexte
consommation de ressources naturelles (Sadok, 2013a; Sadok and Mekerta, 2012), ceci
concorde avec le contexte réglementaire actuel. En effet, dans un article publié en mars 2016
concernant le rejet des boues de dragage polluées, la loi pour l’économie bleue, proposée le
24 mars interdit le rejet en mer de boues de dragage polluées, à partir du 1er janvier 2020
(Boughriet, n.d.).
A ces différents phénomènes vont venir s’ajouter les interactions liant les sédiments aux
polluants. Différentes études menées (Dia et al., 2014; Zri et al., 2009) sur le sujet montrent
que les polluants se fixent sur la fraction fine également, à cause des propriétés
électrostatiques des argiles, et de manière moindre, sur les limons. Les polluants peuvent être
de nature inorganique comme les métaux, certains sels (sulfates, phosphates, nitrates…etc) et
les cyanures, ou ils peuvent être de nature organique comme les hydrocarbures aromatiques
polycycliques (HAP) issus de la biosynthèse, les polychlorobiphényles (PCB) ou le
tributylétain (TBT).
Ces polluants se lient étroitement avec les particules fines minérales et les molécules
organiques, ce qui entraîne, dans le cadre d’un processus de valorisation, l’utilisation des
traitements complexes souvent coûteux, pour détruire les polluants afin de satisfaire aux
seuils-limites de valorisation, car selon le type et les caractéristiques des particules fines, le
potentiel de fixation des métaux sera différent.
36
Chapitre 1 : Etat de l’art et contexte
C’est dans ce contexte, que plusieurs études et projets de recherche ont été initiés ces
dernières années. (Sadok, 2013b). Une charte de préfiguration, datée de 2009, manifeste une
forte volonté commune des signataires, de définir et organiser les conditions de lancement
d’un projet «Sédimatériaux» pour aboutir à terme à l’établissement d’une charte
d’engagement définitive impliquant l’ensemble des acteurs concernés. Ce projet est concrétisé
par la charte Sédimateriaux, présentée à Lille au 5e salon Environord en 2009. L’Ecole
des Mines de Douai a créé en avril 2014 une chaire industrielle sur la valorisation des
sédiments, en partenariat avec des acteurs industriels, portuaires et bureaux d'ingénierie.
L'objectif de cette chaire, baptisée ECOSED (ECOnomie circulaire des SEDiments), est de
créer une dynamique scientifique, technologique et partenariale autour de la gestion des
sédiments portuaires et fluviaux en vue de les recycler en technique routière ou en produits de
béton (Mines Douai, n.d.).
Les Résidus de Procédés Thermiques (RPT) sont les cendres et autres résidus
d’incinérations issus des installations individuelles ou industrielles (incinérateurs, centrales
thermiques). Ce sont des produits en partie chimiquement réactifs que l'on cherche à valoriser
pour éviter leur coûteuse mise en décharge.
Les R.P.T. présentent des gisements renouvelables considérables et leur intégration dans
la filière des travaux publics comme substituts des granulats naturels, est une solution de
valorisation intéressante, compte tenu du tonnage très important et de leur composition
élémentaire, souvent proche des analogues naturels. L’appellation de Résidus de Procédés
Thermiques renvoie à toute une famille de déchets d’incinération, on distingue principalement
les cendres volantes de centrale thermique (CVCT), les mâchefers d’incinération, les sables
de fonderie et les laitiers (d’aciérie et de haut fourneaux).
37
Chapitre 1 : Etat de l’art et contexte
Les cendres volantes, ne sont plus considérées comme déchets, et sont désormais un
produit normalisé. La norme NF EN 197-1 distingue deux types de cendres volantes
(siliceuses ; calciques) :
Les cendres volantes siliceuses sont les cendres les plus produites en France
(environ 70% de la production totale) sont essentiellement constituées de dioxyde
de silice et d’oxyde d’aluminium. Elles présentent des propriétés pouzzolaniques.
Les cendres volantes calciques sont essentiellement constituées d’oxyde de
calcium de dioxyde de silice et d’oxyde d’aluminium, elles présentent des
propriétés pouzzolaniques et/ou hydrauliques. Les cendres volantes contiennent en
outre des imbrûlés, provenant de l’incomplète combustion du charbon dont la
proportion n’excède pas 5%.
La norme NF P 11-300 prend uniquement en compte les cendres volantes siliceuses. Elle
les classe dans la gamme des sous-produits industriels (catégorie F2). Les cendres volantes
sont utilisées pour de nombreuses applications, le plus souvent mélangées à d'autres
constituants. L'addition de cendres volantes au ciment Portland améliore sa maniabilité. Elles
38
Chapitre 1 : Etat de l’art et contexte
peuvent rentrer dans la composition de ciments (NF P 15-301), de liants hydrauliques routiers
(NF P 15-108), de bétons (NF EN 450-1) ainsi que dans de nombreux produits brevetés avec
des utilisations diverses.
Figure 7 : Aspect du mâchefer métallurgique de type clinker à la sortie d'un four (Wikipedia).
Avant d’être valorisés, les MIDND peuvent passer par des opérations de préparation
(extraction des métaux ferreux, extraction des métaux non-ferreux et une période de
maturation). Le Tableau 4 récapitule, pour un tonnage de MIDND considéré (49 sites
connaissant les tonnages à chaque étape de préparation des mâchefers), les quantités de
métaux ferreux et non ferreux récupérées. Les proportions de métaux extraits restent stables
entre 2010 et 2011.
39
Chapitre 1 : Etat de l’art et contexte
Les sables de fonderie peuvent différer par le sable d'origine, et par les liants et adjuvants
qui leur assurent la stabilité mécanique nécessaire pour supporter la coulée.
La valorisation des sables de fonderie date des années 1980 en substitution des sables
naturels notamment en technique routière (Devaux and Vecoven, 1984). Les travaux récents
40
Chapitre 1 : Etat de l’art et contexte
portent sur la valorisation des sables fortement phénolés, et l’amélioration de la gestion de ces
déchets par la mise en place des centres spécialisées dans la préparation de sables de fonderie
afin de fournir des stocks pérennes et homogènes pour utilisateurs du domaine du Génie civil.
Cependant, il y a un manque de réglementation claire relatif à la valorisation des sables de
fonderie à l’exception de la circulaire du ministère de l'Environnement, datant de 1991.
41
Chapitre 1 : Etat de l’art et contexte
Le caractère pollué des déchets et la présence des métaux lourds parfois, sont les
principaux handicaps pour une réutilisation intégrale des déchets. Dans le contexte socio-
économique actuel, qui préconise de réemployer au maximum lesdits déchets dans une
démarche d’économie circulaire et donc de limiter les emprunts extérieurs, pour atteindre un
objectif ambitieux « zéro emprunt, zéro déchet ». L’identification et la caractérisation des
matériaux alternatifs est une étape primordiale pour pouvoir les réutiliser efficacement, mais
elle demeure insuffisante pour une exploitation optimale. En effet, en plus de leurs
connaissances il faut s’adapter à leurs particularités et leurs interactions avec les matériaux
standards.
En effet, le caractère perturbateur d’un élément chimique contenu dans les déchets résulte
de la combinaison de multiples paramètres. On peut identifier les propriétés intrinsèques des
matériaux alternatifs (concentration des éléments perturbateurs, forme chimique…), les
caractéristiques du liant (teneur en clinker, chaux, etc.) mais aussi les facteurs
environnementaux ou les conditions de cure (température, humidité). Ces observations
conduisent aussi à s’interroger sur la notion même de seuil de perturbation, seuil à partir
duquel un élément peut entraîner un échec du traitement. L’aptitude d’un matériau au
traitement ne peut donc pas être définie de manière absolue mais au regard de la combinaison
de certains de ces paramètres.
42
Chapitre 1 : Etat de l’art et contexte
En abordant le problème comme tel, nous proposons dans ce travail des outils
informatiques bien adaptés à l’exploitation de quantité importantes de données et la
modélisation des phénomène complexes, avec un avantage majeur qu’est l’exploitation rapide
et continue des données. Ces outils sont des algorithmes d’apprentissage automatiques et
d’optimisation.
L’évaluation environnementale des matériaux alternatifs ne fait pas l’objet de notre étude
particulièrement, cependant il est important de comprendre ces concepts fondamentaux liés à
l’impact environnemental des déchets dans le cadre de leurs valorisations comme matériaux
alternatifs. En effet, un déchet stocké ou valorisé peut avoir un impact sur l'environnement par
le relargage de ses polluants, c’est pourquoi il est primordial d’étudier l’impact des différents
modes de stockage et de valorisation des déchets sur l'environnement, avant de pouvoir les
stocker ou les valoriser.
Les voies de valorisation des matériaux alternatifs sont conditionnées par leurs natures et
leurs propriétés physico-chimiques et mécaniques. La plupart des matériaux du génie civil
43
Chapitre 1 : Etat de l’art et contexte
(Béton, routes, brique) sont souvent une combinaison de plusieurs matériaux de base. Il est
ainsi possible de formuler des matériaux multi déchets si les conditions techniques le
permettent.
Toutefois, la plupart de ces travaux sont restreints à la substitution partielle d’un des
matériaux standards par un matériau alternatif, jusqu’à un dosage optimal de façon à
minimiser les pertes des propriétés mécaniques et de durabilité. A défaut de méthodes de
formulation plus adaptées, cette démarche donne des résultats relativement satisfaisants mais
empêche les praticiens à prendre ces déchets comme des matériaux à part entière en raison
d’un manque de méthode de formulation adaptées, car les méthodes classiques
n’envisageaient pas cette pratique et demeurent inadaptées.
44
Chapitre 1 : Etat de l’art et contexte
C’est dans ce contexte, qu’il devient important de compléter les travaux des scientifiques
dans leurs démarches, par le développement d’outils et de méthodes de formulation de béton,
adaptés à l’utilisation des matériaux alternatifs sachant que ces données de la littérature
peuvent représenter un capital important à exploiter pour proposer de nouvelles méthodes.
Ceci permettra d’optimiser les quantités de déchets valorisés et écarter le frein lié aux
méthodes d’utiliser les dits déchets efficacement. Cette démarche représente un levier notable
pour la valorisation des déchets dans les applications de béton. Dans la section suivante, nous
abordons en détail l’évolution des méthodes de formulation et l’apport de la méthode
proposée dans cette thèse.
Les critères liés aux coûts énergétiques sont utilisés conjointement avec les possibilités
techniques, pour établir une classification globale des déchets pour une utilisation en
technique routière (Tableau 5), cette classification a bien entendu évolué depuis cette étude
menées par le OCDE (OCDE, 1977) comme les boues de dragage qui sont actuellement plutôt
d’un usage courant en raison des avancées scientifiques et l’évolution de la réglementation
comme expliqué précédemment.
45
Chapitre 1 : Etat de l’art et contexte
Différents travaux menés portant sur la valorisation des matériaux alternatifs en technique
routière, nous citons dans ce qui suit les principaux déchets valorisés dans cette application
routière : Les sédiments de dragage marins et fluviaux (Bourabah, 2012; Dubois, 2006;
Khezami, 2014; Raouf, 2013), les mâchefers d’incinération d’ordures ménagères (Becquart,
2007; Dabo, 2008), les laitiers d’aciéries (Waligora, 2008), les granulats de démolition
(Bauchard, 1984; Mathias, 2005), les sables de fonderie et sous-produits (Devaux and
Vecoven, 1984; François, 1996; OCDE, 1977).
46
Chapitre 1 : Etat de l’art et contexte
eau fixée (Sedran, 1999). C’est pourquoi plusieurs méthodes se sont développées se basant sur
l’optimisation de la compacité granulaire via des courbes de références à l’instar de la
méthode Bolomey, de Caquot ou de Faury (Sadok et al., 2016b), s’inscrivant à ce titre dans le
registre des méthodes mono-objectif qui visent à optimiser les mélanges en se souciant
uniquement de la compacité granulaire comme paramètre de grande influence. Cette approche
a été adoptée dans d’autres méthodes de formulation de béton classiques mono-objectif
également, en optimisant la résistance mécanique à la compression, nous citons à ce titre, la
méthode de Féret. D’autres méthodes se sont développées (Figure 9), qui permettent
d’aborder les méthodes de formulation comme un problème bi-objectif, comme la méthode de
Dreux Gorisse, où on optimise le dosage en ciment en fonction de la résistance mécanique
souhaitée et le dosage en eau en fonction de l’ouvrabilité visée, mais cette optimisation n’est
pas simultanée, ce qui nous empêche d’optimiser les dosages des composants avec une
maitrise totale de toutes les propriétés recherchées dans le béton.
Toutefois, les nouvelles exigences des cahiers de charges modernes nécessitent la prise en
considération et la maitrise simultanée d’une multitude de paramètres souvent contradictoires,
allant des propriétés mécaniques classiques aux caractéristiques environnementales et
économiques, faisant désormais des méthodes de formulation des bétons, des Problèmes
d’Optimisation Multiobjectifs (POM). L’optimisation joue un rôle important dans la
recherche opérationnelle, Plusieurs méthodes de résolution de problèmes sont issues de
l'intelligence artificielle dont les Algorithmes Génétique qui sont présentés dans la patrie 2.
47
Chapitre 1 : Etat de l’art et contexte
En plus, la valorisation des matériaux alternatifs dans les bétons, a dévoilé davantage les
limites des méthodes de formulations classiques, en effet, leurs modèles prédictifs sont pour la
plupart développées sur la base de l’expérience acquise sur les matériaux standards produits
en carrière ou traités sur sites. L’intégration de matériaux alternatifs [composé de déchets
générés] rend souvent ces méthodes de formulation classiques inadaptées en raison des
critères de base considérés, liés en général à favoriser le paramètre de compacité et négligeant
par conséquent, les spécificités des matériaux alternatifs souvent de grande influence sur les
propriétés du béton formulé. C’est dans ce contexte que la nécessité de modéliser les
propriétés des bétons à base de matériaux alternatifs s’impose également. La diversité des
déchets et l’état de connaissance actuel ne permettent pas de prendre les spécificités des
déchets dans les modèles physiques prédictifs classiques, cependant l’apprentissage
automatique permettra sur la base d’expérience de proposer des modèles qui prennent
indirectement les spécificités des matériaux alternatifs. Ces outils d’apprentissage
automatique font partie de la grande famille de l’intelligence artificielle également.
48
Chapitre 1 : Etat de l’art et contexte
Où x ∈ Rn, g ( x ) ∈ Rm et h ( x ) ∈ Rp.
Équation 1
Ici, les vecteurs g ( x ) et h ( x ) représentent respectivement m contraintes d’inégalité et
p contraintes d’égalité.
49
Chapitre 1 : Etat de l’art et contexte
Fonction objectif : C’est le nom donné à la fonction f . C’est cette fonction que l’algorithme
d’optimisation va devoir ―optimiser‖ (trouver un optimum).
Variables de décision : Elles sont regroupées dans le vecteur x . C’est en faisant varier ce
vecteur que l’on recherche un optimum de la fonction f .
Minimum global : Un ―point‖ x * est un minimum global de la fonction f si on a :
f ( x *) < v ( x ) quel que soit x tel que
x * ≠ x . Cette définition correspond au point M3 de
Figure 10.
Minimum local fort : Un ―point‖ x * est un minimum local fort de la fonction f si on a :
f ( x *) < f ( x ) Quel que soit x ∈ V ( x *) et
x *≠ x , où V ( x *) définit un ―voisinage‖ de
x *. Cette définition correspond aux points M2 et M4 de la Figure 10
Minimum local faible : Un ―point‖ x * est un minimum local faible de la fonction f si on a :
f ( x *) ≤ f ( x ) quel que soit x ∈ V ( x *) et
x *≠ x , où V ( x *) définit un ―voisinage‖ de
x *. Cette définition correspond au point M1 de la Figure 10
50
Chapitre 1 : Etat de l’art et contexte
Les Problèmes d'Optimisation Multiobjectifs (POM) sont très fréquents dans le monde
réel. Dans un tel problème, les objectifs à optimiser sont généralement conflictuels
(contradictoire : lorsque l’augmentation d’un objectif entraîne une diminution de l’autre
objectif et vice versa), ce qui signifie qu'il n'y a pas une seule solution de ce problème mais
plutôt un ensemble de solutions. Un POM est résolu lorsque toutes ses solutions Pareto
optimales sont trouvées. Cependant, il est impossible de trouver l'ensemble total de solutions
Pareto-optimales, parce que le nombre de solutions, non-dominées, augmente très rapidement
avec l'augmentation des objectives (Farina et al., 2004). En pratique, l'utilisateur n'a besoin
que d'un nombre limité de solutions bien distribuées le long de la frontière optimale de Pareto
(Voir section ci-dessous), ce qui rend la tâche d'un POM relativement plus facile. Le problème
d'optimisation multiobjectif s’écrit de la manière suivante :
Minimiser f 1 ( x ), f 2 ( x )
Avec g (x)≤0
Et h ( x ) = 0
Où x ∈ Rn, f 1 ( x ), f 2 ( x ) ∈ Rk, g ( x ) ∈ Rm et h ( x ) ∈ Rp
Équation 2
Comme on peut le voir ici, on n’a plus un seul objectif à atteindre, mais deux fonctions-
objectifs. Le but que l’on se fixe dans la résolution d’un POM est de minimiser /maximiser
―au mieux‖ les différents objectifs.
« L’intérêt d’effectuer une optimisation multiobjectifs par rapport à une optimisation
monobjectif est la possibilité de visualiser l’amélioration relative d’un objectif par rapport à
un autre sur un problème test simple. Comme pour le cas des bétons, en fonction des
variables (dosage des composants par exemple) voir l’évolution de la résistance mécanique à
la compression, mais au même temps suivre l’évolution de l’ouvrabilité, le coût, la résistance
à la flexion …etc. »(Sadok et al., 2014).
51
Chapitre 1 : Etat de l’art et contexte
Les méthodes qui permettent de résoudre un POM, sont classées en cinq groupes:
les méthodes scalaires,
les méthodes interactives,
les méthodes floues,
les méthodes exploitant une métaheuristique,
les méthodes d’aide à la décision.
Les AG utilisés dans ce travail sont des métaheuristiques. Les métaheuristiques sont,
généralement, des algorithmes stochastiques itératifs qui permettent via des opérateurs
génétiques d’obtenir la frontière de Pareto. La section 2 de ce chapitre aborde les Algorithme
Génétiques en détail.
52
Chapitre 1 : Etat de l’art et contexte
Les Algorithme Génétique (AG) demeurent l’un des outils les plus adaptés au traitement
des problèmes d’optimisation multiobjectifs. Le nombre important d’articles qui ont été
publiés sur ce sujet en témoigne. De plus, ce domaine est très dynamique et ne cesse de se
développer.
Les AG ont été développés suivant le principe de l’évolution naturelle des espèces : au fil
des générations, les gènes conservés par population les gènes les mieux adaptés au milieu
dans lequel elles vivent.
Une première ―population‖ (constituée d’un certain nombre d’ « individus », c’est à dire
un certain nombre de « jeux de variables ») est générée, souvent aléatoirement mais il est
possible d’initialiser le problème manuellement, les critères d’évaluation correspondant à
chaque individu sont calculés, puis la population subit des transformations (mutation,
croisement, etc.) afin d’arriver à une nouvelle génération (c’est à dire une population
constituée de nouveaux individus), qui est à son tour évaluée, jusqu'à N générations.
L’avantage majeur des AG est que l’ensemble des solutions non-dominées, la frontière de
Pareto, est obtenue en une seule fois. La Figure 12 schématise le fonctionnement des AG.
53
Chapitre 1 : Etat de l’art et contexte
Différentes variantes des AG ont été développées, les principales différences étant les
critères de sélection, pour les individus conservés d’une génération à une autre, ou qui sont
touchés par les procédures de croisement et mutation (Al-Shihri, 2010). En effet, les
paramètres qui caractérisent un AG sont :
La taille de la population,
Le mode de génération de la population initiale (généralement de façon aléatoire),
Le choix du taux de survie, c’est à dire le nombre d’individus survivants d’une
génération à la suivante,
Le choix des individus survivants (on utilise fréquemment la roulette biaisée de
Goldberg, ou la méthode de tournoi),
Le taux et la procédure de croisement (généralement, deux parents choisis
aléatoirement dans la population courante génèrent deux enfants),
Le taux et la procédure de mutation,
Le taux d’élitisme (un certain nombre de meilleurs individus de la population parents
est intégré d’office dans la population enfants),
Le nombre de clones tolérés pour un individu,
L’arrêt de la recherche : soit sur un nombre maximal de générations, soit lorsque les
caractéristiques statistiques (moyenne et écart-type de la population) n’évoluent plus.
L’une des premières méthodes utilisant les algorithmes évolutifs est VEGA (Vector
Evaluated Genetic Algorithm) (Schaffer, 1985). Par la suite, différentes méthodes ont été
développées selon ce principe, dont NPGA (Niched Pareto Genetic Algorithm)(Horn et al.,
1994) qui introduit le concept de niche pour préserver la diversité des individus, puis MOGA
(Multi-Objective Genetic Algorithm) (Fonseca et al., 1993), basée sur l’utilisation du concept
de dominance au sens de Pareto, suivies de la méthode NSGA (Non-Dominated Sorting
Genetic Algorithm) (Srinivas and Deb, 1994) qui reprend les principes de la méthode MOGA
avec une différence majeure : les individus sont triés par rang de non-domination, puis la
méthode SPEA (Strengh Pareto Evolutionary Algorithm)(Zitzler and Thiele, 1999), qui
contrairement à MOGA attribue une force selon la domination, et finalement la méthode
NSGA II (Non-Dominated Sorting Genetic Algorithm II), utilisée dans cette étude, et qui est
actuellement une des plus utilisées : dans cette version de l’algorithme Deb et al. (Deb et al.,
2002, 2000) ont corrigé les défauts des méthodes précédentes. L’optimisation multiobjectifs a
54
Chapitre 1 : Etat de l’art et contexte
été effectuée dans ce travail par les AG implémentés dans « matlabR2010.b » via la toolbox
« gamultiobj ».
Dans l’optimisation d’une fonction, les AG, ne nécessitent que la connaissance des
intervalles de variation des variables prises dans la fonction à optimiser. Ils sont donc très
adaptés dans le cas des fonctions à plusieurs optima locaux, contrairement aux méthodes
déterministes qui s’avèrent défaillantes dans ce cas de figure. Sachant que ce type de
problème complexe est fréquent dans le Génie Civil, cette complexité réside par exemple dans
la non linéarité des fonctions ou dans les problèmes avec des interactions entre les variables.
L’avantage des AG réside également dans leur capacité à travailler dans des espaces différents
(continu, discret ou produit cartésien d’espace continu et discret).
Cependant, l’inconvénient des AG reste bel et bien le temps de calcul qui reste élevé et
dans certain cas prohibitif mais le progrès réalisé en calcul parallèle palie relativement à cela.
Reste donc que le paramétrage optimale des AG, en effet dans cette étape l’utilisateur doit
fixer un grand nombre de paramètres (taille de la population, probabilité d’application des
opérateurs de croisement et mutation, les critères d’arrêt) et dont la mise au point de ces
différents paramètres ne peut se faire qu’empiriquement, car il n y a pas de méthode exacte
pour le faire, ceci nécessite donc le plus souvent de procéder à de nombreux tests avant
d’obtenir un résultat optimal, ce qui sera au détriment du coût de calcul encore une fois.
Dans l’optimisation multiobjectifs, les AG se basent sur la fonction fitness qui regroupe
les différentes fonctions-objectifs, ces fonctions sont (dans le cas d’optimisation des mélanges
des bétons) les modèles prédictifs des propriétés fixées par le cahier des charges. Dans
l’absence de ces modèles, l’optimisation par les AG ne peut s’effectuer, mais comme la
formulation des bétons est un problème ancien dans le domaine du génie civil, différents
modèles prédictifs existent développés sur la base des matériaux standards. Pour conclure, la
limite des AG réside aussi dans sa dépendance des modèles prédictifs, ou du moins, la qualité
des résultats des AG dépend de la qualité des modèles prédictifs utilisés, d’où la nécessité
dans notre cas d’étude de béton à base de matériaux alternatifs, de développer des modèles
prédictifs des propriétés des bétons adaptés à l’utilisation des matériaux alternatifs. Ce travail
de modélisation est assez complexes, c’est pourquoi on a fait recours à l’Apprentissage
Automatique (Machin Learning).
55
Chapitre 1 : Etat de l’art et contexte
Nous définissons les deux principaux types d’apprentissage, supervisé et non supervisé,
dans le premier type « supervisé » l’apprentissage se fait après traitement et classification des
données par le superviseur, contrairement à l’apprentissage non supervisé (clustering) ou les
données sont traitées automatiquement par un logiciel avec l'objectif de permettre une
extraction de connaissance organisée à partir de ces données brutes.
D’autre types d’apprentissage existent comme le semi-supervisé, partiellement supervisé,
l'apprentissage par renforcement et l'apprentissage par transfert.
56
Chapitre 1 : Etat de l’art et contexte
plus longue est l'analyse, d’où l’intérêt des plans d’expériences pour élaborer les bases
de données dans le cadre d’une étude paramétrique.
Nombre et qualité des attributs décrivant ces exemples : Les attributs numériques
sont plus faciles à traiter que les attributs catégoriels (couleur, beauté, utilité).
Pourcentage de données renseignées et manquantes : Certains algorithmes
d’apprentissage sont plus sensibles aux données manquant que d’autres.
« Bruit » : Le nombre et emplacement des valeurs douteuses (erreurs potentielles,
valeurs aberrantes) parmi l’ensemble des données.
La qualité de l'apprentissage dépend également de la complexité du modèle (spécifique ou
généraliste), plus il est généralise, plus on a besoin de données et moins l’apprentissage est
pertinent.
La Programmation Génétique (PG) fait partie des algorithmes évolutionnaires (Figure 13).
Elle est basée sur le même principe d'évolution que les AG, mais on l’applique à des
programmes.
57
Chapitre 1 : Etat de l’art et contexte
Comme la plupart des méthodes évolutionnaires, la PG commence son évolution par une
population dite « initiale ». Cette dernière sera améliorée génération après génération en lui
appliquant des opérateurs génétiques (Figure 15)
58
Chapitre 1 : Etat de l’art et contexte
59
Chapitre 1 : Etat de l’art et contexte
Selon les conclusions de Gandomi et Alavi (Gandomi and Alavi, 2011a) , les solutions
proposées par la MGGP sont capables de simuler efficacement les comportements non
linéaire des systèmes étudiés et les résultats de MGGP se trouvent plus précis que ceux de la
PG standard et des modèles proposés par les réseaux de neurones artificiels.
4 Démarche globale
La méthode proposée dans cette thèse est une méthode hybride des deux algorithmes
présentés dans les sections précédentes, les AG comme algorithmes d’optimisation
multiobjectifs combinés avec la PG qui permet de fournir par apprentissage automatique les
modèles prédictifs qui représentent les fonctions d’évaluation pour les AG. La Figure 16
récapitule la méthodologie globale.
60
Chapitre 1 : Etat de l’art et contexte
L’analyse SWOT (des initiales des mots anglais Strenghts, Weaknesses, Opportunities,
Threats) ou analyse des forces, faiblesses, opportunités et menaces est dans notre cas d’étude,
un outil qui permet de cerner l’intérêt d’utiliser ces algorithmes. Il constitue ainsi une aide à la
prise de décision afin d’orienter la stratégie adoptée dans notre démarche.
Le plus souvent, l’analyse SWOT se présente sous la forme d’un tableau. Elle doit être
très factuelle et concise : l’objectif n’étant pas de refaire une liste des opportunités-menaces,
forces-faiblesses mais de mettre en évidence les éléments les plus pertinents et les plus
influents sur notre choix d’utiliser les Algorithmes évolutionnaires (PG Tableau 7 et AG
Tableau 8).
FORCES FAIBLESSES
L’apprentissage aboutit à des modèles Dépendent de la qualité de la base de données
mathématiques pratiques. d’apprentissage.
Pas besoin d’avoir une forme initiale de la
solution.
OPPORTUNITES MENACES
Robustesse face aux données manquantes. Les réseaux de Neurones Artificielles.
Tableau 7 : Analyse SWOT de la Programmation Génétique.
FORCES FAIBLESSES
Optimiser simultanément des objectifs Résultats dépendant des modèles prédictifs
contradictoires en un temps raisonnable utilisés
OPPORTUNITES MENACES
Très adaptés lorsque les espaces de recherches Le recuit simulé
sont importants La recherche taboue
Tableau 8: Analyse SWOT de la Algorithmes Génétique.
61
Chapitre 1 : Etat de l’art et contexte
1 Objectifs et démarches
62
Chapitre 1 : Etat de l’art et contexte
Ainsi le Chapitre 1 : représente un état de l’art sur la valorisation des déchets dans le
domaine du Génie Civil, en situant l’étude dans le contexte générale lié à l’évolution de la
réglementation nationale et internationale sur la gestion des déchets, ensuite ce même chapitre
est scindé en trois parties bibliographiques :
Partie 1 : Recense les principaux déchets valorisés dans le domaine du Génie Civil, avec
des explications sur leurs origines, les quantités produites ainsi que les projets nationaux et
travaux scientifiques menés sur le potentiel de valorisation de chaque déchet, ensuite cette
partie est complétées sur les principales voies de valorisation ainsi que les méthodes de
formulation utilisées notamment pour le cas béton qui fait particulièrement l’objet de cette
étude.
Partie 2 : En continuité aux méthodes de formulation des bétons, cette partie introduit la
nouvelle méthode proposée dans cette thèse. La compréhension de cette méthode requiert une
étude bibliographique sur les deux outils principaux utilisés dedans. A cet effet, les concepts
de base sur les AG comme outils d’optimisation, et la PG comme outil d’apprentissage
automatique, sont présentés, pour finir avec un récapitulatif de la démarche globale. Cette
partie est bouclée par une analyse S.W.O.T sur les algorithmes utilisés.
Partie 3 : Présente les contours et la structuration de la thèse.
Chapitre 2 : Comme chapitre intermédiaire entre la partie bibliographique et la partie
expérimentale et modélisation, ce chapitre porte sur la caractérisation des différents matériaux
standards et alternatifs utilisés dans cette thèse. Les résultats de la caractérisation, sont
comparés aux matériaux similaires évoqués dans la première partie de premier chapitre, et les
résultats de caractérisation sont utilisés comme données nécessaires dans les formulations des
bétons mais également dans le développement des modèles prédictifs, comme paramètres
d’entrée dans certain cas.
Chapitre 3 : Porte sur la modélisation de la compacité granulaire par Apprentissage
Automatique, avec l’intégration du modèle développé dans une interface graphique, ce
chapitre est scindé en deux parties :
Partie 1 : Création d’une base de données expérimentale parallèlement à une étude
paramétrique de la compacité granulaire afin de déterminer les paramètres pertinents sur la
compacité granulaire des matériaux standards, cette partie permet d’évaluer les performances
de la PG dans d’apprentissage automatique.
Partie 2 : Compléter la base de donnés par des données relatives aux matériaux alternatifs
(Granulats recyclés, Sédiments de dragage et Sable de fonderie) de différentes natures, dans le
63
Chapitre 1 : Etat de l’art et contexte
but d’optimiser et généraliser le modèle prédictif de la compacité granulaire sur les différents
types de granulats (standards et alternatifs).
Chapitre 4 : Dans ce chapitre, on combine la PG et Algorithmes Génétique pour proposer
une méthodologie de formulation globale, en commençant par mettre en pratique les
connaissances développées sur la PG, dans la modélisation de la résistance mécanique et
l’ouvrabilité d’un Béton Hautes Performances, sur une base donnée de la littérature
scientifique. Ensuite une optimisation multiobjectif est réalisée par Algorithme Génétique afin
de satisfaire un cahier des charges typique d’un BHP. La méthodologie globale est ensuite
validée expérimentalement au laboratoire par des formulations de BHP.
Dans cette partie, nous menons également une étude spécifique sur la répétabilité des AG
et la PG, ensuite nous nous intéressons à l’écriture mathématique des modèles prédictifs
développés par la PG, et leurs impacts sur l’optimisation multiobjectifs.
Chapitre 5 : Avec la même démarche suivie au chapitre précédent, le chapitre 5 applique
la méthodologie PG & AG pour un cas d’étude des mortiers à base de matériaux alternatifs,
avec en plus, la base de données d’apprentissage réalisée par nos soins, elle comporte 96
formulations de mortiers réalisées au laboratoire selon un plan d’expériences optimisé. Il est à
noter également que contrairement au cas des BHP abordés au chapitre 4 qui est un cas
d’optimisation bi-objectif, le cas de mortier est un problème d’optimisation à cinq objectifs.
64
Chapitre 1 : Etat de l’art et contexte
Conclusion du chapitre
65
CHAPITRE 2: Caractérisation des
matériaux
Chapitre 2 : Caractérisation des matériaux
Introduction du chapitre
Nous avons retenu du chapitre 1 que les matériaux standards et alternatifs sont d’origines
diverses et variées. Ceci leur procure des propriétés très différentes, c’est pourquoi il est
primordial de bien les connaitre via une étude de caractérisation.
L’objectif de ce chapitre est l’établissement des fiches d’identité des différents matériaux
utilisés dans cette étude. Ces essais de caractérisation présentés dans ce chapitre, sont bien
évidemment orientés directement vers l’objectif de l’étude qui porte sur la modélisation des
propriétés les plus recherchées dans les mortiers et les bétons (Compacité du squelette
granulaire, résistance à la compression, résistance à la flexion et maniabilité) par
apprentissage automatique.
Ces fiches d’identité nous permettent de déterminer certaines caractéristiques physico-
chimiques, géotechniques ou minéralogiques de chaque matériau en fonction du besoin. Ces
différents paramètres mesurés peuvent être pris comme des inputs dans les modèles prédictifs
développés, par conséquent ils seront utilisés comme des paramètres d’entrées à prendre en
compte dans une démarche d’optimisation des mélanges.
Enfin, le choix des essais permettra de classer les matériaux, afin de s’assurer de la
pertinence d’un tel matériau dans les formulations des mortiers selon ses propriétés, en se
référant aux normes en vigueur.
Rappelons que les matériaux de l’étude sont des matériaux standards : ciment, fillers,
cendres volantes, graviers et sables de carrières mais aussi des matériaux alternatifs :
sédiments de dragage, sables de fonderie et granulats de béton recyclé.
Tous les essais de caractérisation ont été réalisés au sein du Laboratoire Génie Civil et
Environnemental de Mines Douai.
68
Chapitre 2 : Caractérisation des matériaux
1 Matériel et mesures
Dans cette première section, nous présenterons d’abord les principes de base des
principaux essais, les appareils et les normes utilisées. Dans la seconde section, les résultats
obtenus pour chaque matériau seront discutés.
1.1 La granulométrie
69
Chapitre 2 : Caractérisation des matériaux
Chaque matériau est composé de grains solides, de vides (pores intra granulaires et vide
inter granulaires) et d’eau éventuellement. Chaque phase représente un volume partiel comme
montré sur la Figure 18
70
Chapitre 2 : Caractérisation des matériaux
1.3 La compacité
La compacité mesurée à sec est le rapport du volume solide (Vs) d’un mélange granulaire
sur son volume apparent (Vt), ou en termes de masse volumique elle peut être exprimée
comme décrit dans l’Équation 3 :
Équation 3
Lorsqu’on étudie les matériaux à granulométrie fines (poudres), Les mesures des surfaces
spécifiques sont très importantes. En effet, la surface développée par une poudre conditionne
71
Chapitre 2 : Caractérisation des matériaux
la demande en eau nécessaire pour mouiller l'ensemble des particules et donne une
représentation de la surface offerte aux adjuvants chimiques (Sadok et al., 2016c).
Le principe de l’essai de Blaine est de faire passer une quantité d’air connue au travers
d’un lit de poudre (ciment ou autre matériau fin) et de mesurer le temps de passage de l’air,
dans le but de déterminer la finesse du broyage du matériau. En effet, le temps (t)
d’infiltration de l’air est proportionnel à la surface massique du matériau. Cet essai est décrit
dans la norme NF EN 196-6 en utilisant l’appareil présenté sur la Figure 21 (A gauche).
La méthode BET (Brunauer, Emmett et Teller) est une méthode qui permet la
détermination de la surface BET par adsorption d'azote. L’appareil BET utilisé dans notre
étude est présenté sur la Figure 21 (A droite). Ce matériel est de type micromeritics-
ASAP2010.
72
Chapitre 2 : Caractérisation des matériaux
Les analyses minéralogiques consistent à déterminer les phases minérales d’un matériau
par diffraction aux rayons X. Cette méthode consiste à utiliser la propriété de
diffraction d’un faisceau monochromatique de rayons X par les plans réticulaires. Le
matériel utilisé est de type D8 fabriqué par Siemens, présenté sur la Figure 23.
73
Chapitre 2 : Caractérisation des matériaux
Les matériaux standards caractérisés dans cette partie, englobent les différents granulats
utilisés dans l’étude de la compacité granulaire, ainsi que les matériaux utilisés dans les
formulations des BHP et des mortiers. Le ciment utilisé comme liant hydraulique, les
additions et les adjuvants sont également caractérisés.
La caractérisation des matériaux, permet de bien identifier les caractéristiques des
matériaux utilisés, mais aussi de prendre en compte les plus pertinents comme paramètres
d’entrée (Input) dans les modèles prédictifs.
Les granulats utilisés dans cette étude (Figure 24), sont des granulats calcaires
concassés du Gaurain, commercialisés en tant que classes granulaires : sable Gn 0/4 (A),
gravillons Gn 4/8 (C), Gn 8/14 (B) et Gn4/16. On a utilisé également des granulats calcaires
concassés Bn 14/20 (D) qui proviennent de la carrière du Boulonnais.
Les échantillons ont été préparés suivant les prescriptions de la NF EN 932-1 (méthodes
d'échantillonnage) et les analyses granulométriques faites conformément à la norme NF
EN 933-1 (Analyses granulométriques). Les courbes granulométriques des granulats sont
74
Chapitre 2 : Caractérisation des matériaux
présentées sur la Figure 25 suivi d’un tableau récapitulatif des propriétés granulométriques de
chaque classe (Tableau 9)
100
90
80
70
Gn 0/4
Passant%
60
Gn 4/8
50
Gn 8/14
40
Gn 4/16
30
Bn 14/20
20
10
0
0,01 0,1 1 10
Diamètre (mm)
Dmoy*: Diamètre moyen calculé par la formule de Mellor (Jeanneau, 1967) (Sadok et al.,
2016b) présentée en l’Équation 6
( )( ))
∗ √
Équation 6
75
Chapitre 2 : Caractérisation des matériaux
Pour chaque classe granulaire, des mesures des masses volumiques apparentes, réelles et
absolues ont été effectuées ainsi que les mesures des propriétés mécaniques, les résultats sont
regroupés dans le Tableau 10.
ρ (t/m3) ρr (t/m3) ρs (t/m3) LA MDE
NF EN NF EN NF P EN EN
1097-3 1097-6 94-054 1097-2:1998 1097-1
Gaurain Gn 0/4 1,61 2,66 2,69 18 16
Gaurain Gn 4/8 1,54 2,66 2,68 18 16
Gaurain Gn 8/14 1,56 2,66 2,69 18 16
Gaurain Gn 4/16 1,67 2,67 2,69 18 16
Boulonnais Bn 14/20 1,57 2,69 2,72 22 13
Tableau 10 : Masse volumiques des granulats d’étude.
Les résultats de caractérisation obtenus sur les granulats utilisés révèlent des matériaux de
carrière courants utilisés d’une manière standard pour confectionner des bétons de différentes
classes de résistance.
Le ciment utilisé dans cette thèse est un liant hydraulique, nous avons opté pour le ciment
CEM1 52.5 N MF, de la cimenterie Holcim d’Obourg. C’est un ciment Portland dont le
constituant principal est le clinker Portland (K). La teneur en clinker est supérieure à 95 %.
Mesure Valeurs
Rc2 (MPa) 34
Rc28(MPa) 62
2
Blaine (m /gr) 0,367
2
BET (m /gr) 0,91
3
ρs (t/m ) 3,21
Compacité (NF EN 196-3) 0,526
Tableau 11 : Fiche d’identité du ciment d’étude 52,5 N.
76
Chapitre 2 : Caractérisation des matériaux
D’après les résultats présentés dans le Tableau 11, nous remarquons que le ciment utilisé à
une masse volumique de 3,21 t/m3et une compacité égale à 0,526. Ces résultats sont courants
pour des ciments portland de même classe.
60
50
40
30
Ciment
20
52.5 N
10 63 um
0
0,1 1 10 100 1000
Diamètre (µm)
77
Chapitre 2 : Caractérisation des matériaux
Des essais de FX ont été réalisés également pour détecter les principaux éléments
chimiques qui constituent ce ciment. Les résultats sont récapitulés dans le Tableau 12, les
constituants principaux sont l’oxygène, le calcium, le silicium, le fer et l’aluminium.
Élément Ciment 52.5 N
O 38,9
Na 0,5
Mg 0,7
Al 2,7
Si 7,6
P 0,2
S 1,8
Cl 0,1
K 0,6
Ca 43,5
Ti 0,2
V Traces
Cr Traces
Mn Traces
Fe 2,8
Cu Traces
Zn 0,2
Sr 0,1
Zr Traces
Ba Traces
Total 99,8
Tableau 12 : Composition chimique du ciment obtenue par fluorescence X (% massique).
78
Chapitre 2 : Caractérisation des matériaux
En raison de leur finesse comme de leur réactivité, les additions minérales sont
couramment utilisées dans les bétons, car ils apportent une amélioration très significative des
propriétés des bétons et des produits dérivés. Il était donc important d’envisager leur
utilisation dans les méthodes de formulation des bétons développés dans ce travail. En effet,
nous souhaitons développer tant que possible des modèles fiables autant dans le cas
d’utilisation des gros granulats (gravier, sable) que dans le cas d’utilisation des matériaux fins
(additions minérales).
Cendres volantes
100
90
80
70
Passant %
60
50
40
30
20 Cendres
Volantes
10
63 um
0
0,1 1 10 100 1000
Diamètre (µm)
79
Chapitre 2 : Caractérisation des matériaux
L’analyse granulométrique par granulomètre Laser montre que les cendres volantes ont
une granulométrie bien étalée avec un diamètre maximum ne dépassant pas 1mm (Figure 28).
Notons que la surface spécifique Blaine des cendres utilisées égale à 0,293 m²/gr, cette valeur
est proche de celle enregistrée pour le ciment 0,367 m²/gr. Nous remarquons également une
P.A.F très faible de 0,51%, ceci est assez courant pour les matériaux similaires issus des
procédés thermiques.
Fillers calcaires
100
90
80
70
Passant %
60
50
40
30
20
filler
10
63 um
0
0,1 1 10 100 1000
Diamètre (µm)
Fillers Bn dmin Dmax (mm) D10 (mm) D30 (mm) D50 (mm) D60 (mm) Cu Cc
Dmoy*
(mm) (mm)
0,0005 0,23 0,0008 0,004 0,014 0,021 26,25 0,95 0,14
Tableau 15 : Granulométrie des fillers calcaires.
80
Chapitre 2 : Caractérisation des matériaux
Nous constatons à travers ces résultats, que les fillers utilisés sont de granulométrie très
fine, d’un Dmax ≤1mm avec au moins 85% de passant à 0,063mm. Ceci est en concordance
avec les résultats de la surface spécifique Blaine égale à 0,396 plus importante que celle du
ciment et des cendres volantes respectivement de 0,367 et 0,293.
IL est à noter que la finesse et la morphologie des fillers influent sur leur demande en eau,
sur leurs propriétés rhéologiques et sur leur insertion dans le squelette granulaire du ciment et
par conséquent sur l'existence et l'intensité de l'effet filler.
Nous avons utilisé deux adjuvants dans ce travail plus précisément dans la formulation des
Bétons à hautes performances. Ces adjuvants sont un superplastifiant et un entraineur d’air.
2.4.1 Superplastifiant
Nous avons utilisé CHRYSO® Fluid Optima 224, c’est un superplastifiant haut réducteur
d’eau de nouvelle génération, à base de polycarboxylate modifié. Selon le fabriquant, ce
superplastifiant est particulièrement recommandé et adapté à la formulation de bétons auto-
plaçant homogènes ayant une capacité de remplissage élevée. Il permet de réaliser des bétons
avec un long maintien d’ouvrabilité sans retard de prise et il est compatible avec la majorité
des ciments.
Le Tableau 17 récapitule les propriétés du superplastifiant CHRYSO® Fluid Optima 224
Nature Liquide
Couleur brun
Densité (20° C) 1,05 ± 0,02
pH 4,70 ± 1,00
Teneur en ions Cl¯ ≤0,10 %
Na2O équivalent ≤1,00 %
Extrait sec 21,60 % ± 1,10 %
Tableau 17 : Caractéristiques du superplastifiant utilisé.
81
Chapitre 2 : Caractérisation des matériaux
Nature Liquide
Couleur Jaune claire
Point de congélation 0 °C
Densité (20° C) 1,005 ± 0,005
pH 7,00 ± 1,00
Teneur en ions Cl¯ ≤ 0,10 %
Teneur en Na2O équivalent ≤ 0,50 %
Extrait sec 3,6 % ± 0,3 %
Tableau 18 : Caractéristiques de l’agent entraineur d’air.
Dans le cadre de la présente étude qui porte sur la mise en place de méthode de
formulation adaptée aux matériaux standards et alternatifs, il est primordial de connaitre ces
derniers au même titre que les matériaux standards afin de bien cerner leurs particularités et
de confirmer la possibilité de les utiliser suivant la réglementation en vigueur. Les matériaux
alternatifs utilisés dans cette étude sont les Graviers de Béton Recyclé (GBR), les sables de
fonderie, les sédiments et sables de dragage du Port de Dunkerque (sédiments marins) et les
sédiments de la Deûle (sédiments fluviaux).
Les GBR utilisés sont des Graviers de classe 4/16 fournis par Colas de la plateforme
d’Amiens. Sur la Figure 30 présentant un zoom sur un GBR, on peut voir que ce granulat est
un matériau composite constitué de granulats naturels concassés partiellement et d’une pâte
de ciment hydraté concassée, enrobant les granulats naturels. Le taux de pâte de ciment
présent dans ces GBR, est fonction du béton parent et de la taille des granulats (Zhao, 2014).
82
Chapitre 2 : Caractérisation des matériaux
100
90
80
70
Passant%
60
50
40
30
20
10
GBR 4/16
0
0,1 1 10
Diamètre (mm)
83
Chapitre 2 : Caractérisation des matériaux
Les résultats de l’absorption et les masses volumiques sont récapitulés dans le Tableau 20.
On constate d’après le Tableau 20 que les GBR présentent des propriétés physiques et
mécaniques sensiblement inférieures à celles des granulats naturels. En effet, ceci est relevé à
la comparaison de la masse volumique absolue et la compacité des GBR respectivement de
2,51 t/m3 et 0,585 conte celle obtenues pour les granulats naturels Gn4/16 respectivement de
2,67 t/m3 et 0,620. Ces caractéristiques faibles sont dues à la présence importante des pores
intra granulaire dans les GBR. Quant aux propriétés mécaniques, ceci est traduit par la
résistance à la fragmentation par l’essai Los Angeles et la résistance à l’usure par l’essai
Micro Deval respectivement de 28,7% et 24,6%. Ces résultats sont dus à la pâte durcie du
ciment (PDC) rattachée aux granulats recyclés plus friable, qui affecte les propriétés du GBR
(Berredjem et al., 2015).
Pour une utilisation dans le béton, la valeur de la résistance mécanique, évaluée à partir de
l’essai Los Angeles, pour les GBR reste au-dessous des limites indiquées par les normes les
plus restrictives, c’est-à-dire 40 et 50
L’absorption d’eau des GBR avoisinant les 6% retient particulièrement notre attention en
raison des valeurs habituelles pour les matériaux standards qui varient seulement entre 1 et
3%, ceci est due à la présence de l’ancien mortier dans les GBR qui est poreux, ceci concorde
avec les travaux de nombreux auteurs (Padmini et al., 2009, 2009; Rao et al., 2007; Zhao,
2014) qui estiment l’absorption des GBR entre 3 et 12%.
84
Chapitre 2 : Caractérisation des matériaux
Les sables de fonderie utilisés dans ce travail proviennent l’écopôle Recydem, filiale de
SITA à Lourches (59). La Figure 32 et Tableau 21 récapitulent leurs granulométries
60
50
40
30
20
10 Sab.fond brut
0
0,01 0,1 1 10
Diamètre mm
Figure 32 : Courbe granulométrique des sables de fonderie.
Nous remarquons, d’après ces résultats que les propriétés des sables de fonderie sont assez
similaires des sables standards, la masse volumique absolue (ρs) des sables de fonderie
mesurée est de 2,67(t/m3) contre 2,69 (t/m3) pour le Gn0/4. Cependant nous remarquons que
l’absorption des sables de fonderie est relativement élevée, estimée à 4,5% selon la norme EN
1097-6, ce qui est au dessus des valeurs habituelles enregistrées pour des sables de carrière.
85
Chapitre 2 : Caractérisation des matériaux
Selon l’arrêté du 16 juillet 1991 relatif à l'élimination des sables de fonderie, les sables de
fonderie utilisés dans ce travail sont en dessous du seuil maximal en teneur en phénol
autorisé, fixé à 5mg/kg pour la fabrication de produits à base de liants hydrauliques, et à
1mg/kg pour une utilisation dans les remblais.
60
50
40 Inf à 0,2 mm
30 0.063 mm
20
10
0
0,0001 0,001 0,01 0,1
Diamètre mm
Les sables de dragage utilisés proviennent de Dunkerque, ils contiennent des résidus de
coquillage. L’analyse granulométrique est récapitulée dans Figure 34 et le Tableau 24
86
Chapitre 2 : Caractérisation des matériaux
Sable de dragage
100
90
80
70
Passant %
60
50
40
30
20
10
Sable de dragage
0
0,001 0,01 0,1 1 10
Diamètre (mm)
Nous remarquons à travers ces résultats que le sable de dragage est un sable fin avec plus
de 80% de passant à 1mm. Selon le coefficient d’uniformité Cu le sable de dragage utilisé
présente une granulométrie étalée.
Nous remarquons à travers le Tableau 25, que le sable de dragage est très propre d’après
la valeur d’équivalent de sable égale à 92 ce qui concorde avec l’analyse granulométrique. De
plus la valeur de P.A.F reflète l’absence matière organique dans ce matériau.
87
Chapitre 2 : Caractérisation des matériaux
Afin de couvrir les différentes tranches granulaires, allant des graviers recyclés en passant
par les sables de fonderie et de dragage, des matériaux de granulométrie fine, qui sont les
sédiments de dragage, ont été utilisés. Deux sortes de sédiments sont utilisés, les sédiments de
dragage marins (Sédiments de Dunkerque) et sédiments fluviaux (Sédiments de la Deûle).Une
série d’analyses a été menée pour caractériser les deux types de sédiments.
En raison de la plasticité des sédiments, ces derniers se présentent comme des mottes,
pour les utiliser secs, on a procédé à un séchage ensuite on les a émiettés pour en obtenir une
poudre fine de sédiments. L’opération a été réalisée par le concasseur de laboratoire à
mâchoire présenté à la Figure 35
88
Chapitre 2 : Caractérisation des matériaux
Sédiments de Dunkerque
100
90
80
70
Passant %
60
50
40
30
20
10 Sed. Dunkerque
63 um
0
0,1 1 10 100 1000 10000
Diamètre (µm)
Les résultats obtenus confirment la nature fines des sédiments de dragage de Dunkerque,
avec un Dmax inférieur à 1mm et plus de 60% de passant inférieur à 0,063 mm.
Les résultats obtenus sur les sédiments de Dunkerque présentés sans le Tableau 27,
permettent de voir principalement un forte P.A.F estimée à 10,2 % reflétant un matériau
89
Chapitre 2 : Caractérisation des matériaux
moyennement organique, ceci concorde avec les travaux de (Dia et al., 2014; Khezami, 2014;
Raouf, 2013) sur le caractère pollué des sédiments marins. Les résultats des limites
d’Atterberg révèlent un matériau plastique avec un indice de plasticité Ip=15,8.
Sédiments de la Deûle
100
90
80
70
60
Passant %
50
40
30
20
10 Sed.Deule
63 um
0
0,1 1 10 100 1000 10000
Diamètre (µm)
90
Chapitre 2 : Caractérisation des matériaux
Les résultats obtenus sur les sédiments de la Deûle présentés dans le Tableau 29, révèlent
un matériau non plastique avec un IP=2, de granulométrie plutôt sableuse avec environ 20%
de passant à 0,063 mm et faiblement organique avec une P.A.F égale à 7,8% moins
importante que celle enregistrée pour les sédiments de Dunkerque égale à 10,2%. C’est
résultats sont en concordance avec les travaux antérieurs sur les sédiments fluviaux (Sadok,
2013c; Thanh, 2009)
Conclusion du chapitre
Pour les matériaux standards, nous avons opté pour des matériaux locaux assez communs
et qui sont utilisés d’une manière standard pour confectionner des bétons de différentes
classes de résistance. Pour les matériaux alternatifs, nous avons opté pour des matériaux
d’origines et de caractéristiques très différentes avec un double objectif : être capable de
remplacer différentes classes granulométriques dans les formulations des bétons et mortiers
envisagés, mais également d’utiliser des matériaux alternatifs susceptibles d’être
représentatifs d’autres déchets ou sous-produits industriels abondants.
Il est aussi utile de rappeler l’intérêt d’intégrer dans cette étude de valorisation des
matériaux finis de types sédiments de dragage ou de curage, des sables de fonderie et des
granulats de bétons recyclés qui intègrent des spécificités présentes dans d’autres coproduits
ou sous-produits industriels, comme la présence des matières organiques, l’absorption des
granulats, la non-conformité de la distribution granulométrique ou la forme des matériaux.
Pour les matériaux alternatifs adoptés dans ce travail, il apparait que selon les
caractéristiques de base mesurées, ces derniers peuvent répondre aux exigences d’une
utilisation pour la formulation des bétons ou des mortiers.
91
92
CHAPITRE 3 : Modélisation de la
compacité granulaire
Chapitre 3 : Modélisation de la compacité granulaire
Introduction du chapitre
Avec l’avènement de bétons avec des performances supérieures à celles des bétons
anciens et des usages très techniques, les méthodes basées sur l’utilisation des courbes
maitresses ont montré rapidement leurs limites (Miraoui, 2010). Dans ce contexte, plusieurs
travaux de recherche ont été engagés pour proposer des méthodes ou des modèles de
prédiction plus généralistes.
Ainsi au cours des trois dernières décennies, plusieurs modèles ont été développés pour
prédire la compacité d’un mélange granulaire. Ces modèles souvent complexes nécessitent la
détermination de paramètres parfois difficilement accessibles par l’expérience. Ce qui réduit
considérablement l’intérêt des modèles dans des applications pratiques. De plus, les modèles
développés à ce jour ont été fait sur la base d’études menées sur des matériaux standards et
classiques.
Le modèle ainsi développé peut être comparé au modèle proposé par De Larrard (De
Larrard, 2000) à la base du logiciel de formulation RENE-LCPC. Comme souligné dans
l’introduction générale de ce manuscrit, l’approche adoptée est basée sur la PG.
Pour répondre à cet objectif et au regard des informations disponibles, il était important
dans un premier temps de développer une base de données qui intègre aussi bien des mesures
de compacité sur des matériaux standards que sur des matériaux alternatifs.
Pour s’assurer du caractère général des modèles que nous souhaitons développer,
plusieurs paramètres sont variés dans les mesures effectuées. Parmi ces paramètres, on peut
noter la variation dans une large gamme des rapports des diamètres moyens des mélanges
étudiés, la variation des modes de serrage imposés aux mélanges, la nature des matériaux
(fine, sable, gravillon et gravier) et l’allure des courbes granulométriques.
Dans ce qui suit, après une description rapide de la base de données mise en place et la
description des modes de serrages adoptées. Une analyse des résultats obtenus en termes de
compacité a permis de dégager les paramètres pertinents de la modélisation engagée.
94
Chapitre 3 : Modélisation de la compacité granulaire
Dans un dernier temps, nous avons successivement proposé une modélisation et une
comparaison des modèles issus de la PG ou la MGGP avec les prédictions du logiciel RENE-
LCPC.
Il est important de rappeler que les modèles déduits par la PG ou la MGGP dans ce
chapitre sont adoptés dans le chapitre 5 pour enrichir la méthode de formulation multiobjectif
adoptée dans le cadre de cette étude.
95
Chapitre 3 : Modélisation de la compacité granulaire
Pour assurer un caractère général à cette base de données dans cette étude, plusieurs
paramètres (contrôlables) sont modifiés. L’objectif de ce travail est d’identifier les paramètres
de forte influence sur la compacité d’un mélange granulaire à savoir la nature des matériaux
étudiés et les modes de compactage.
Dans cette étude paramétrique, nous avons utilisé différents matériaux standards et
alternatifs pour la création de la base de données expérimentale (Tableau 30). Le choix des
matériaux a été fait sur la base de leurs caractéristiques granulométriques de façon à étudier le
comportement des matériaux fins jusqu’aux matériaux grossiers. Avec la même logique, le
choix des matériaux alternatifs a été fait. La caractérisation de ces matériaux est donnée en
détail au chapitre 2.
Matériaux
Matériaux standards Matériaux alternatifs
96
Chapitre 3 : Modélisation de la compacité granulaire
Afin d’étudier l’influence du mode de serrage sur la compacité des mélanges binaires
étudiés, quatre modes de compactage ont été adoptés à savoir : versement simple, piquetage,
vibration avec une pression de 10 kPa et table à choc avec pression de 10 kPa.
Il est à noter que ces différents modes imposent des énergies de serrage différentes aux
mélanges. Dans ce cadre, un indice de serrage (K) a été proposé par De Larrard (De Larrard,
2000), par calibrage entre des données théoriques du Modèle d’Empilement Compressible
(MEC) qu’il propose, et les résultats expérimentaux. La valeur de l’indice de serrage K a été
ajustée de façon à obtenir la meilleure coïncidence possible entre prédictions du modèle et
résultats expérimentaux. Les valeurs de K correspondant aux différents modes de serrage
proposés sont récapitulées dans le Tableau 31.
Dans ce travail, nous retenons, les quatre modes de mise en place pour la mesure de la
compacité des classes élémentaires comme celle des mélanges binaires, afin de confirmer les
similitudes observées pour les deux modes : Vibration+10 kPa et Table à choc+10 kPa et les
différences observées en comparaison du versement simple et du piquetage.
Méthode de versement simple (NF EN 1097-3): C’est une méthode qui consiste, à
remplir un récipient de volume V connu, en laissant s'écouler le matériau sous l'effet de
la gravité (Figure 38), sans aucune autre sollicitation. Après arasement et pesée de
l’éprouvette, la masse volumique apparente est calculée. La compacité est calculée selon
l’Équation 3 du chapitre 2.
97
Chapitre 3 : Modélisation de la compacité granulaire
Méthode de piquage ASTM C29-78: cette méthode est inspirée de la méthode ASTM
C29-78 comme montré à la Figure 39. Pour les gravillons, la tige utilisée est en acier
mesurant 600 mm de long avec un diamètre de 16 mm. Pour les sables, la longueur utile de la
tige est de 130 mm avec un diamètre de 6 mm. Les cylindres sont remplis en trois couches
successives, chacune piquée vingt-cinq fois avec la tige adéquate. De même à la fin de la
procédure, la masse volumique de l’éprouvette est calculée puis la compacité est déduite.
98
Chapitre 3 : Modélisation de la compacité granulaire
équipé d’un disque interne de même diamètre qui permet d’appliquer une charge de
10kPa uniformément répartie sur toute la surface des granulats. L’ensemble est solidarisé à
une table vibrante. Après homogénéisation de l’échantillon (Figure 40) ce dernier est versé
dans le cylindre. La vibration est lancée après l’application de la pression (10 kPa) (El
Barrak, 2005). L’échantillon est soumis à une vibration de 150 Hz pendant deux minutes.
A la fin de cette étape, la différence de hauteur provoquée par le compactage par rapport à
la hauteur initiale (Figure 41) est mesurée. La compacité est ainsi déduite.
99
Chapitre 3 : Modélisation de la compacité granulaire
Table à choc avec pression de 10 kPa : Les moules utilisés et les dispositifs de serrage
sont les mêmes que le précédent essai. L’ensemble est solidarisé à une table à choc (Figure
42). Pour appliquer les chocs, il suffit de tourner la manivelle en bas de la table. Les cylindres
sont remplis en trois couches successives, et chaque couche subit 20 chocs. En fin d’essai une
pression totale de 10 kPa sur l’échantillon est appliquée et fait subir 20 chocs à
l’échantillon (De Larrard et al., 2003).On mesure de même le volume de l’échantillon puis on
calcule sa masse volumique apparente et enfin la compacité.
100
Chapitre 3 : Modélisation de la compacité granulaire
Pour ce qui est de l’influence du mode de serrage sur la compacité des classes
élémentaires, on remarque d’après la Figure 43, que la compacité maximale est plus grande
pour les fines, dès qu’elles sont soumises à une sollicitation mécanique. Pour ce qui est de
l’efficacité du compactage, il est clair qu’indépendamment des dimensions des granulats, la
table à choc+10 kPa donne la compacité maximale par rapport aux autres procédés de
compactage testés (versement simple, piquetage, vibration+10kPa).
Nous remarquons également, que contrairement au sable et aux fillers, le piquetage est
plus efficace pour les graviers que la vibration+10kPa. Ceci est dû a la granularité importante
des graviers, causant une friction plus importante entre les grains pendent la mise en place par
vibration, contrairement au compactage par piquetage, où la tige s’enfonce dans le matériau
permettant un réarrangement plus facile des grains. Ceci est en contradiction avec les
conclusions de plusieurs travaux (De Larrard, 2000), sur l’amélioration inconditionnelle de la
compacité des mélanges avec l’augmentation de l’indice de serrage K, logiquement plus
important pour la vibration +10 kPa (K=9) que pour le piquetage (K= 4,5). Nous concluons
donc que certes l’influence du mode de mise en place est confirmée, mais l’efficacité du mode
de compactage est conditionnée, dans certains cas par la granulométrie (ou les
caractéristiques) du matériau également.
0,700
Compacité réelle
0,600
0,500
0,400
0,300
Filler Bn Gn 0/4 Gn 4/8 Gn 8/14 Bn 14/20
Matériaux
Versement simple Piquetage Vibration+10KPa Choc+10KPa
Pour ce qui est de l’influence du mode de mise place sur la compacité des matériaux
alternatifs, nous observons des comportements similaires vis-à-vis de l’énergie de serrage, que
101
Chapitre 3 : Modélisation de la compacité granulaire
les matériaux standards comme exposé sur la Figure 44, où en remarque une augmentation de
la compacité des matériaux en augmentant l’énergie de serrage et en remarque également que
la table à choc induit une compacité presque similaire que la « vibration10+kPa » notamment
pour les sables de fonderie et les granulats des bétons recyclés. Ces premiers résultats
permettent de conclure que le mode de serrage est un paramètre de forte influence sur la
compacité des différentes classes élémentaires, de matériaux standards ou alternatifs.
Il est à noter également que le mode de serrage a encore plus d’influence sur les fractions
fines que grossières. Enfin la comparaison entre les résultats obtenus sur les matériaux
alternatifs et standards est comparable pour les sables et les graviers toutefois très différente
pour les fines.
0,700
Compacité réelle
0,600
0,500
0,400
0,300
Sed-Dunk Sab-Fond GBR
Matériaux
Versement simple Piquetage Vibration+10KPa Choc+10KPa
Nous avons effectué une étude systématique sur des mélanges binaires obtenus par une
combinaison de deux classes élémentaires, préalablement préparées. Pour la réalisation des
mélanges, nous avons procédé à une homogénéisation manuelle. Trois mélanges binaires de
matériaux standards ont été préparés. Pour la réalisation de ces mélanges nous avons choisi
des proportions Y1 (du composé le plus fin du mélange) allant de 10% à 100% avec un
incrément successif de 10%. Les résultats sont présentés dans les Figure 45 à Figure 47.
102
Chapitre 3 : Modélisation de la compacité granulaire
0,700
0,600
0,500
0,400
0 20 40 60 80 100
Pourcentage des fines Y1
Versement Piquetage Vibration+10 kPa Choc+10KPa
Figure 45 : Variation de la compacité en fonction des proportions des grains fins (Y1)
pour des mélanges binaires (Filler Bn-Gn4/8) et du mode de compactage.
Il apparait à partir de la Figure 45 que l’ajout des fillers dans le mélange (Filler Bn-Gn4/8),
améliore sensiblement la compacité jusqu'à un taux avoisinant les 40%, au-delà de ce seuil,
l’augmentation du pourcentage du composé le plus fin (Y1), diminue la compacité. Ceci est
valable pour les différents modes de compactage adoptés, ce résultat s’explique par l’insertion
des grains de la classe de petite dimension (Fillers dans ce cas-là) dans le vide inter-granulaire
de gravier du Gaurain (Gn4/8) augmentant ainsi la compacité du mélange, jusqu’à
comblement de tous les vides à environ 40% de Fillers. Au-delà de ce pourcentage optimal, la
baisse de la compacité proportionnellement au pourcentage Y1, est expliqué par le
phénomène de desserrement des gros grains par les petits grains.
Il est à noter également que la forme des courbes en forme de cloche, s’accentue avec
l’augmentation de l’énergie de serrage, ceci concorde avec les conclusions de Bouterfas
(Bouterfas, 2012) sur l’évolution de la compacité des mélanges binaires en fonction des
modes de serrage.
Les mêmes remarques ont été soulevées sur le mélange binaire Gn0/4-4/8, présenté sur la
Figure 46 avec différents mode de compactage.
103
Chapitre 3 : Modélisation de la compacité granulaire
0,900
Mélange binaire Gn 0/4-4/8
d1/d2=2,45
Compacité 0,800
0,700
0,600
0,500
0,400
0 20 40 60 80 100
Pourcentage des fines Y1
Versement Piquetage Vibration+10 kPa Choc+10KPa
Pour le mélange binaire (Gn 0/4-Gn 4/8), les courbes obtenues se présentent sous forme
de cloche également et les résultats montrent bien que l’optimum de la compacité se trouve à
un taux d’introduction des grains fins avoisinant les 50%. Pour la majorité des mélanges,
indépendamment du mode de compactage, la variation de la compacité est plus rapide du côté
des gros grains dominants (petits pourcentages de grains fins) que du côté des petits grains
dominants. Cette variation s’atténue quand les classes granulaires sont plus proches.
Dans la Figure 47, les résultats du mélange de deux classes élémentaires relativement
proches en granulométrie sont présentés.
104
Chapitre 3 : Modélisation de la compacité granulaire
0,900
Mélange binaire Gn 8/14-Bn 14-20
d1/d2=1,52
Compacité 0,800
0,700
0,600
0,500
0,400
0 20 40 60 80 100
Pourcentage des fines Y1
Versement Piquetage Vibration+10 kPa Choc+10KPa
Figure 47 : Variation de la compacité en fonction des proportions des grains fins pour
des mélanges binaires (Gn8/14-Bn14/20) et du mode de compactage.
Pour ce mélange de classes granulaires de plus grande dimension par rapport aux
précédents mélanges (Gn 8/14-Gn14/20) on remarque bien l’avantage qu’offre le compactage
par piquetage en comparaison avec la vibration+10 kPa. La compacité maximale est encore
une fois atteinte par compactage à la table à choc. Contrairement aux précédents mélanges, on
constate que la proportion de chaque classe élémentaire n’induit pas un changement sur la
compacité. Cependant le mode de compactage reste un facteur important. Ceci nous renvoie,
par conséquent, à étudier l’effet de l’étendue granulaire d2/d1 sur la compacité, qui peut être
un paramètre pertinent.
La même démarche effectuée pour les granulats standards a été entreprise avec des
matériaux alternatifs, où on a réalisé trois mélanges granulaires binaires (Sed.Dunk-GBR,
Sed.Dunk-Sab.Fond et Sab.Fond-GBR) à différentes fractions et compacités avec différents
mode. Nous n’avons pas réalisé le compactage par table à choc+10 kPa sur les matériaux
alternatifs puisque ce mode de compactage s’est avéré d’une énergie assez similaire à la table
de vibration+10 kPa. Les résultats obtenus pour les trois mélanges sont représentés sur les
Figure 48, Figure 49 et la Figure 50.
105
Chapitre 3 : Modélisation de la compacité granulaire
0,900
Mélange binaire Sed.Dunk-GBR
d2/d1=35,86
Compacité 0,800
0,700
0,600
0,500
0,400
0 20 40 60 80 100
Pourcentage des fines Y1
Versement simple Piquetage Vibration+10KPa
Figure 48 : Variation de la compacité en fonction des proportions des grains fins pour
des mélanges binaires (Sed.Dunk-GBR) et du mode de compactage.
0,900
Mélange binaire Sed.Dunk-Sab.Fond
d2/d1=6,46
0,800
Compacité
0,700
0,600
0,500
0,400
0 20 40 60 80 100
Pourcentage des fines Y1
Versement simple Piquetage Vibration+10 kPa
Figure 49 : Variation de la compacité en fonction des proportions des grains fins pour
des mélanges binaires (Sed.Dunk- Sab.Fond) et du mode de compactage.
106
Chapitre 3 : Modélisation de la compacité granulaire
0,900
Mélange binaire Sab.Fond- GBR
d2/d1=5,54
Compacité 0,800
0,700
0,600
0,500
0,400
0 20 40 60 80 100
Pourcentage des fines Y1
Versement simple Piquetage Vibration+10KPa
Figure 50 : Variation de la compacité en fonction des proportions des grains fins pour
des mélanges binaires (Sab.Fond-GBR) et du mode de compactage.
A l’observation des différentes courbes obtenues pour les mélanges granulaires binaires
des matériaux alternatifs, nous tirons les mêmes conclusions que celles des matériaux
standards, à savoir :
107
Chapitre 3 : Modélisation de la compacité granulaire
Afin d’étudier l’effet de l’étendue granulaire sur la compacité des mélanges binaires, les
rapports des diamètres moyens de chaque classe sont rappelés récapitulées dans le Tableau 32
et les résultats des compacités pour un mode de compactage sont représentés sur la Figure 51.
Nous observerons à partir de la Figure 51, l’effet de l’étendue granulaire sur la compacité
des différents mélanges, pour le compactage par vibration+10kPa.
Nous observons à partir de la Figure 51 que pour des mélanges avec un rapport des
diamètres moyens importants, un effet important des proportions des composants fin (Y1) sur
l’augmentation des capacités des mélanges (pour les matériaux standards et alternatifs). Pour
des mélanges impliquant des matériaux grossiers de granulométrie comparable, peu d’effet
sur l’évolution des compacités n’était à signaler. Pour des mélanges impliquant des matériaux
dont les rapports des diamètres moyens sont modérés, il apparait que les formes des courbes
108
Chapitre 3 : Modélisation de la compacité granulaire
granulométriques ou même les classes granulométriques peuvent avoir une influence sur
l’évolution des compacités (Hanini, 2012)
Le bilan de cette étude pour évaluer les paramètres influents sur la compacité des
mélanges binaires permet de confirmer l’impact important de l’énergie de serrage, le rapport
des diamètres moyens des constituants des mélanges mais également la forme des courbes
granulométriques des constituants.
Dans la suite de l’étude, ces paramètres seront intégrés comme des paramètres pertinents
pour l’identification du modèle prédictif par la programmation génétique.
109
Chapitre 3 : Modélisation de la compacité granulaire
Considérons :
– C i: Compacité mesurée de la classe i,
– C'i: Compacité prédite de la classe i par le modèle,
– : Moyenne de compacité mesurée de la classe i,
– : Moyenne de compacité prédite de la classe i par le modèle.
– N: Le nombre de tests réalisés
∑ ( )
√
Équation 7
110
Chapitre 3 : Modélisation de la compacité granulaire
∑ ( )( )
√∑ ( ) ∑ ( )
Équation 8
∑ | |
Équation 9
Comme il n y a pas un seul critère sur quoi se baser pour évaluer la pertinence d’un
modèle, la combinaison de ses trois critères énumérés ci-dessus, permet d’évaluer d’une
manière satisfaisante la capacité de prédiction du modèle développé.
∑ | |
Équation 10
111
Chapitre 3 : Modélisation de la compacité granulaire
3.2 Paramétrage de la PG
Comme tout algorithme évolutionnaire, un ensemble de paramètres doit être ajusté pour
initialiser, faire et arrêter le processus de recherche de solution optimale. La PG n’échappe
pas à ce type de paramétrage.
De plus, ces algorithmes nécessitent de connaitre les paramètres pertinents (appelés dans
ces méthodes les terminaux) et les fonctions en plus de la base de données qui sert de terrain
d’entrainement pour affiner et proposer un modèle mathématique.
Dans le Tableau 33 sont reportés les principaux paramètres fixés dans la toolbox GPLAB:
Paramètres Valeurs
112
Chapitre 3 : Modélisation de la compacité granulaire
L’application de la PG a été entreprise sur les résultats issus de l’étude des matériaux
standards. Ce choix est fait pour comparer sur la même base de données, les prédictions d’un
modèle issu de la PG et la prédiction du logiciel RENE LCPC où le modèle d’empilement
compressible est intégré.
Il est rappelé que les outils utilisés dans ce travail ont la particularité de proposer, pour un
même environnement, d’apprentissage, différentes solutions. Dans le cadre de ce travail
plusieurs simulations ont été réalisées successivement. Pour les différentes simulations, les
critères définis ci-dessus quant à l’évaluation des modèles ont été utilisés pour sélectionner le
meilleur modèle.
Nous présentons dans la Figure 52, le meilleur modèle trouvé sous forme d’arbre
syntaxique, décrivant la compacité résiduelle moyenne d’un mélange en fonction des
paramètres intrinsèques des constituants.
Le meilleur individu retenu présente une Fitness égale à 2,8, un coefficient de corrélation
R de 0,92, un RMSE de 0,03 et un MAE de 0,03. Il aura été insensé de trouver une fitness
nulle, vu que la base de données d’apprentissage a été réalisée expérimentalement, ceci
implique des erreurs de mesures dû à l’opérateur et à la répétabilité de l’essai de compacité
qui reste tout de même satisfaisante avec des écarts types très faibles.
113
Chapitre 3 : Modélisation de la compacité granulaire
[ ( )] ∗ [ ( )]
Équation 11
Avec : ( )( )
Les coefficients A, B et C sont des constantes qui dépendent des propriétés intrinsèques
des composants à savoir : (d2/d1) l’étendue granulaire du mélange et le (Cc1/Cc2) rapport des
coefficients de courbure des deux classes élémentaires.
114
Chapitre 3 : Modélisation de la compacité granulaire
Les deux figures ci-dessous présentent les comparaisons des prédictions du modèle PG et
du logiciel RENE LCPC par rapport aux données expérimentales. Nous remarquons en
premier que les prédictions du logiciel RENE LCPC sont relativement satisfaisantes sur les
matériaux standards, en se basant sur l’erreur absolue moyenne de 0,05 ainsi qu’un RMSE et
un R respectivement égale à 0,07 et 0,86.
115
Chapitre 3 : Modélisation de la compacité granulaire
théoriquement le modèle idéal. Le modèle issu de la PG donne un MAE de 0,03 ainsi qu’un
RMSE et un R respectivement égale à 0,06 et 0,92 qui expriment une liaison fortement
linéaire entre les résultats théorique du modèle et les résultats expérimentaux.
En plus des performances du modèle, nous soulignons le fait que le modèle issu de la PG
présente le grand avantage d’être développé en fonction des paramètres intrinsèques des
composants, facilement accessibles à l’expérience à savoir l’allure des courbes exprimée par
le coefficient de courbure Cc et l’étendue granulaire exprimé par le diamètre moyen,
contrairement au modèle d’empilement compressible dans le logiciel RENE LCPC qui exige
plus de paramètres dont certains sont difficilement accessibles à l’expérience (Hanini, 2012).
Afin d’appréhender et de comparer les performances des modèles développés sur des
matériaux standards quant à la prédiction des compacités mesurées sur des mélanges de
matériaux alternatifs, les mêmes exercices qui précédent on été entrepris. Dans la Figure 55
nous comparons les prédictions du modèle issu de la PG Équation 11) avec les prédictions du
logiciel RENE LCPC sur les matériaux alternatifs sont présentés.
On note que les comparaisons présentées sur la Figure 55, démontre une certaine
dispersion des résultats pour les deux modèles, en comparaison avec leurs prédictions
116
Chapitre 3 : Modélisation de la compacité granulaire
Une analyse plus fine des résultats pointe parfois quelques lacunes (Figure 56) pour la
prédiction des compacités des mélanges qui intègrent des matériaux fins comme les sédiments
de dragage.
0,7
Compacité
0,6
0,5
0,4
0,3
0 10 20 30 40 50 60 70 80 90 100
0,7
0,6
0,5
0,4
PG Expérimental RENE
0,3
0 10 20 30 40 50 60 70 80 90 100
Pourcentage des fines Y1
117
Chapitre 3 : Modélisation de la compacité granulaire
Le modèle obtenu à l’issu de cette optimisation, est présenté dans l’Équation 12:
( ∗ ) ( ∗ ) ( ∗( )) ( ∗ ( ))
∗
( ∗ ( )) ( ∗ )
Équation 12
Les performances du modèle, issu de la MGGP développé par apprentissage sur matériaux
standards et alternatifs, sont confrontées aux résultats de prédiction du modèle précédent issu
de la PG, les résultats sont présentés sur la Figure 57 :
Sur la base des coefficients de corrélation R ainsi que l’erreur absolue moyenne MAE, les
performances du présent modèle sont assez comparables à celui issu de la PG.
118
Chapitre 3 : Modélisation de la compacité granulaire
4 L’interface GUImines
L’interface GUImines (Guide User Interface des mines) présentée dans la Figure 58, a été
créée afin de faciliter l’utilisation du modèle issu de la PG développé dans ce travail.
A défaut d’un modèle plus général, adapté aux mélanges de plus de deux constituants,
nous proposons d’utiliser cet outil pour une optimisation par étapes dans les cas des mélanges
de trois constituants et plus. C’est-à-dire optimiser le dosage des deux grandes classes
granulaires du mélange, et continuer le processus d’optimisation avec la classe des fins,
comme expliqué sur la Figure 59:
119
Chapitre 3 : Modélisation de la compacité granulaire
Nous avons procéder dans cette section à la formulation de quatre bétons ordinaires dont
le cahier des charges est : Rc28 égal à 40MPa et l’affaissement égal à 10 cm. Trois
formulations (F1, F2 et F3) sont à base de matériaux standards et une formulation à base de
matériaux alternatifs (F4).
Il est à noter que le mode de serrage utilisé est le piquage, de coefficient « K=4,5 ». Ce
mode de serrage est le même utilisé dans la formulation de nos bétons.
120
Chapitre 3 : Modélisation de la compacité granulaire
Les matériaux utilisés sont caractérisés dans le chapitre 2 et les résultats des formulations
sont récapitulés dans le Tableau 35.
Les présents résultats nous permettent de constater que les formulations F1 et F2 à base de
matériaux standards (Gn4/16 et Gn0/4), sont assez similaires si ce n’est une légère
amélioration de la résistance mécanique pour la formulation F2, dû à l’amélioration de la
compacité granulaire par le modèle issu de la PG.
L’ajout du filler calcaire comme addition minérale dans la formulation F3 a été réalisée
par le GUImines. L’amélioration de la compacité a engendré une nette augmentation de la
résistance mécanique est constatée qui atteint 52 MPa. Une perte de l’affaissement est
constatée aussi due à la demande en eau des particules fines des fillers. Ceci concorde avec les
différents travaux qui ont démontré que les fillers procurent au béton une meilleure
pompabilité et résistance à la ségrégation et contribuent à l’augmentation de la résistance en
compression (El Hilali et al., 2006). Nous confirmons donc la capacité du modèle PG
implémenté dans GUImines à formuler des mélanges ternaires contentant même des
matériaux fins.
121
Chapitre 3 : Modélisation de la compacité granulaire
d’absorption des GBR estimée à 5,3%, qui n’a pas pu être absorbée entièrement pendent le
malaxage et a comblé la demande en eau des fillers calcaires.
122
Chapitre 3 : Modélisation de la compacité granulaire
Conclusion du chapitre
Pour modéliser la compacité granulaire, nous avons utilisé dans le présent travail la PG
dans ces deux versions (PG et MGGP). Ces outils d’apprentissage automatique, sur la base
d’expériences, permettent de modéliser des problèmes complexes. Nous avons en premier lieu
déterminé les Inputs et les Outputs du système sur la base d’une étude paramétrique. Après
paramétrage de l’algorithme, l’apprentissage est lancé et nous proposons un modèle prédictif
de la compacité d’un mélange granulaire à concentration donnée de ce mélange, à partir de
grandeurs physiques des composants aisément accessibles à l'expérience.
Le modèle de compacité développé par la PG présente une forte capacité de prédiction sur
les matériaux standards avec un coefficient de corrélation R=0,92, un RMSE de 0,06. Dans la
validation du modèle, sur les matériaux alternatifs, de différents nature et dosages, le modèle
de la PG perd légèrement capacité de prédiction en raison de l’apprentissage effectué
uniquement sur les mélanges standards.
Le modèle développé par MGGP (la version améliorée de la PG) par apprentissage sur
matériaux standards et alternatifs, améliore légèrement les performances du modèle issu de la
PG. Il faut toutefois signaler l’avantage de la simplicité des modèles proposés du fait qu’ils
intègrent comme paramètres d’entrées uniquement la granularité (Cc, dmoy), la compacité
propre et le dosage des composants. Ces modèles serviront de fonction-objectif dans
l’optimisation du squelette granulaire dans l’approche d’optimisation multiobjectif des
mélanges de béton proposée dans le chapitre suivant.
Le modèle issu de la PG a été implémenté dans une interface graphique Matlab, pour une
utilisation plus simple pour l’utilisateur. Ceci nous a permis de remplacer la courbe de
référence dans la méthode de formulation de Dreux Gorisse. Des formulations de béton ont
été testées en laboratoire.
123
CHAPITRE 4 : Apprentissage
automatique et optimisation
multiobjectifs
Chapitre 4 : Apprentissage automatique et optimisation multiobjectifs
Introduction du chapitre
Dans la seconde partie présentée dans le chapitre 5, nous utiliserons la même démarche
(méthode hybride PG & AG) sur le cas de la formulation d’un mortier à base de matériaux
alternatifs. Pour ce second cas, nous développerons notre propre base de données
d’apprentissage, et nous introduirons dans le processus d’optimisation d’autres objectifs
comme la compacité granulaire, le coût de formulation en plus de la résistance mécanique et
la maniabilité.
126
Chapitre 4 : Apprentissage automatique et optimisation multiobjectifs
Pour le présent cas d’étude, qui porte sur la formulation des BHP, une base de données de
la littérature scientifique est utilisée (Lim et al., 2004). Cette base servira d’environnement
d’apprentissage dans la modélisation de la résistance mécanique et l’affaissement de ce BHP
par la PG. Ensuite, une optimisation Bi-objectifs sera entamée pour trouver les meilleures
formulations de béton en matière de résistance mécanique et d’affaissement par la PG.
La modélisation par la PG (PG classique ou MGGP), sur les mêmes données, est
susceptible de donner des modèles de formes différentes, ce qui est normal, vu que la PG ou
la MGGP sont des algorithmes stochastiques. La question qui s’impose par conséquent, est de
savoir si ces modèles prédictifs de la même propriété (de formes différentes mais de
performance similaires), donneront des résultats similaires dans la détermination des solutions
optimales par les AG.
Pour répondre à cette question, nous allons donc développer plusieurs modèles pour
prédire la résistance mécanique et plusieurs modèles pour prédire l’affaissement également
sur la même base de données d’apprentissage. Ensuite, nous combinerons les différents
modèles développés dans l’optimisation bi-objectifs et nous comparons les différentes
solutions optimales (frontières de Pareto) obtenues, afin d’évaluer l’impact des formes des
modèles sur le processus d’optimisation.
127
Chapitre 4 : Apprentissage automatique et optimisation multiobjectifs
Outputs Inputs
Les caractéristiques des matériaux utilisés dans les différentes formulations sont
présentées dans le Tableau 37.
Afin de valider la démarche globale par des essais au laboratoire, nous avons adopté des
matériaux de caractéristiques comparables, comme montré sur le Tableau 37.
128
Chapitre 4 : Apprentissage automatique et optimisation multiobjectifs
310 367
SiO2 (%) 21,3 18,2
Al2O3(%) 4,7 5,7
Fe2O3(%) 3,1 4,1
CaO (%) 63,1 62,9
MgO(%) 3 0,9
Dmax 19 20
granulats
Gros
129
Chapitre 4 : Apprentissage automatique et optimisation multiobjectifs
Les modèles issus de la PG sont retenus en fonction de leurs fitness (Maruyama et al.,
n.d.; Sadok, 2016; Zentar et al., 2016) (Figure 60). Le modèle de résistance mécanique à la
compression retenu est représenté en forme d’arbre syntaxique (Figure 61).
130
Chapitre 4 : Apprentissage automatique et optimisation multiobjectifs
131
Chapitre 4 : Apprentissage automatique et optimisation multiobjectifs
(X1) est le rapport eau/liant (%), (X2) le dosage en eau (kg/m3), (X3) le rapport massique des agrégats fins (%), (X4)
cendres volantes (%), (X5) Agent entraîneur d’air (kg/m3) et le (X6) le dosage en super Plastifiant (kg/m3).
Figure 61: Modèle prédictif de la résistance mécanique obtenu par PG présenté en arbre
syntaxique.
Le décodage de l’arbre syntaxique se fait par une méthode bien précise, où globalement il
se fait de la gauche vers la droite et du haut vers le bas (Bouguettaya, 2012) Cette opération
fastidieuse sources importantes d’erreurs a été automatisé dans la MGGP (Sadok et al.,
2016d). Pour l’interprétation de l’arbre présenté dans la Figure 61 en équation mathématique
(Équation 13), nous obtenons un modèle de la résistance mécanique du béton à 28 jours en
fonction du dosage des composants :
( ) ( ) ∗ ∗ (( ) )
( ) ∗ (( ) )
( ) ∗( )
∗ ( )
( )∗ ∗
Équation 13
E : Dosage en eau (kg/m3), E/L : Rapport eau/liant (%), S/G: Rapport massique des agrégats fins (%),
CV: Cendres volantes (%), AE : Agent entraîneur d’air (kg/m3), SP : Super Plastifiant (kg/m3).
132
Chapitre 4 : Apprentissage automatique et optimisation multiobjectifs
Avec la même logique, nous utilisons dans cette partie une variante de la PG, qui est la
MultiGene Genetic Programming (MGGP), dans la modélisation de la résistance mécanique à
la compression d’un BHP. Les mêmes données d’entrée et un paramétrage adéquat à la
MGGP comme décrit sur le Tableau 39 sont utilisés:
Paramètres Valeur utilisée
Taille de la Population 500
Nombre de génération 1000
Nombre maximum de gène dans chaque individu 6
Tournament 20
Elitisme 0.5%
Fonctions utilisées +, −, ×, /, ✓, exp, ln, sin, cos, log, tanh
Tableau 39 : Paramétrage principal de la MGGP.
( ) ( ∗ ) ( ∗ ) ( ∗ ) ( ∗( ))
( ∗( )∗ ∗ ) ( ∗( )∗( )∗( ))
Équation 14
E : Dosage en eau (kg/m3), E/L : Rapport eau/liant (%), S/G: Rapport massique des agrégats fins (%),
CV: Cendres volantes (%), AE : Agent entraîneur d’air (kg/m3), SP : Super Plastifiant (kg/m3).
Avec le même paramétrage décrit dans le Tableau 39, nous avons relancé l’apprentissage
une seconde fois, et nous avons obtenu le modèle Rc-MGGP2 présenté dans l’Équation 15
( ) ( ∗ ) ( ∗ ) ( ∗ ∗ )
∗
( )
( )
∗ ∗( )
( ∗ ∗ ( ))
( )
Équation 15
E : Dosage en eau (kg/m3), E/L : Rapport eau/liant (%), S/G: Rapport massique des agrégats fins (%),
CV: Cendres volantes (%), AE : Agent entraîneur d’air (kg/m3), SP : Super Plastifiant (kg/m3).
133
Chapitre 4 : Apprentissage automatique et optimisation multiobjectifs
Comme évoqué au chapitre 3, les modèles prédictifs développés ont besoin d’être évalués
dans leurs prédictions par rapport aux résultats expérimentaux obtenus comme présentés sur
les Figure 62 à Figure 64. Nous utilisons les même critères définis auparavant pour évaluer les
modèles à savoir ; le coefficient de corrélation (R), l’Erreur Absolue Moyenne (MAE) et la
Racine de l’Erreur Quadratique Moyenne (RMSE).
Performances du modèle
prédictif de la résistance mécanique (Rc28-PG1)
80
R=0,90
MAE=3,38
Prédictions de Rc28 (MPa)
RMSE=4,15
70
60
50
40
30
30 40 50 60 70 80
Valeurs expérimentales de Rc28 (MPa)
Performances du modèle
prédictif de la résistance mécanique (Rc28-MGGP1)
80
R=0,97
MAE=1,81
Prédictions de Rc28 (MPa)
RMSE=2,20
70
60
50
40
30
30 40 50 60 70 80
Valeurs expérimentales de Rc28 (MPa)
134
Chapitre 4 : Apprentissage automatique et optimisation multiobjectifs
Performances du modèle
prédictif de la résistance mécanique (Rc28-MGGP2)
80 R=0,98
MAE=1,32
Prédictions de Rc28 (MPa)
RMSE=1,74
70
60
50
40
30
30 40 50 60 70 80
Valeurs expérimentales de Rc28 (MPa)
Il ressort des résultats montrés (Figure 62 à Figure 64), que la prédiction des modèles
obtenus par les deux approches est très satisfaisante. En effet, le modèle Rc-PG1 présente une
corrélation avec les résultats expérimentaux de R=0,90, une erreur moyenne MAE=3,38 MPa
et un RMSE de 4,15 MPa, ce qui est d’emblée satisfaisant pour un essai similaire. Sur le plan
expérimental, la norme en vigueur relative à la résistance en compression des bétons, admet
une erreur type à 4MPa.
Le modèle Rc-MGGP1, plus simple dans la forme par rapport au modèle Rc-PG1,
présente des résultats meilleurs avec un coefficient de corrélation, une erreur moyenne et un
RMSE respectifs de 0,97, 1,81 MPa et 2.20 MPa par rapport aux résultats expérimentaux. Le
modèle Rc-MGGP2 quant à lui donne de meilleurs résultats que les deux premiers, cependant
il est de forme très compliqué. Il est à noter qu’il est pratiquement impossible d’obtenir des
résultats meilleurs que celui-là, vu que la base de données est générée expérimentalement et
les erreurs obtenues sont dues à l’opérateur et à la répétabilité de l’essai de compression.
135
Chapitre 4 : Apprentissage automatique et optimisation multiobjectifs
Les présents résultats et constats concordent avec les travaux de Gandomi et Alavi
(Gandomi and Alavi, 2011a) sur la meilleure qualité des modèles obtenus par la MGGP en
comparaissant avec la PG classique ou les Réseaux de neurones Artificiels.
Dans la même logique, nous avons procédé à la modélisation de l’affaissement des BHP,
par apprentissage automatique sur les mêmes données que celles adoptées pour la résistance à
la compression. On déduit de façon similaire, un modèle prédictif par PG et un second modèle
prédictif par MGGP.
Après paramétrage de la PG comme décrit sur le Tableau 40, nous avons obtenu un
modèle prédictif de l’affaissement d’un BHP décrit dans l’Équation 16
( ) ∗[ ] ∗ [( ) ∗ ∗ ( )] ( ) ∗
Équation 16
E : Dosage en eau (kg/m3), E/L : Rapport eau/liant (%), S/G: Rapport massique des agrégats fins (%),
CV: Cendres volantes (%), AE : Agent entraîneur d’air (kg/m3), SP : Super Plastifiant (kg/m3).
136
Chapitre 4 : Apprentissage automatique et optimisation multiobjectifs
( ) ( ∗ ) ( ∗ ) ( ∗ ) ( ∗ ) ( ∗ ) ( ∗ | |) (
∗ ( )) ( ) ( ∗ ) ( ∗ ∗ ) ( ∗ ∗ ) ( ∗ ∗ )
( )
( ∗ )
( ∗ ∗ ) ( ) ( ∗ ) ( ∗ ) ( ∗ )
Équation 17
Avec E : Dosage en eau (kg/m3), E/L : Rapport eau/liant (%), S/G: Rapport massique des agrégats fins
(%), CV: Cendres volantes (%), AE : Agent entraîneur d’air (kg/m3), SP : Super Plastifiant (kg/m3), (α, β, γ)
sont des constantes dont α= 0,554, β= 0,0003, γ= 0,0002.
Sur la Figure 65 et la Figure 66, la capacité de prédiction des deux modèles développés,
par PG et MGGP, appelés respectivement Aff-PG1 et Aff-MGGP1 sont comparés. Les
mêmes critères statistiques d’évaluation sont utilisés également pour cette évaluation à
savoir ; le coefficient de corrélation (R), l’Erreur Absolue Moyenne (MAE) et la Racine de
l’Erreur Quadratique Moyenne (RMSE).
Performances du modèle
prédictif de l'affaisement (Aff-PG1)
280
Prédictions de l'affaissement (mm)
R=0,5
MAE=20
RMSE=27,5
240
200
160
120
120 160 200 240 280
Valeurs expérimentaes mesurées Aff (mm)
137
Chapitre 4 : Apprentissage automatique et optimisation multiobjectifs
Performances du modèle
prédictif de l'affaisement (Aff-MGGP1)
280
Prédictions de l'affaissement (mm) R=0,55
MAE=22
RMSE=25,5
240
200
160
120
120 160 200 240 280
Valeur expérimental mesurées (mm)
Au vu des résultats obtenus, nous remarquons que les deux modèles Aff-PG1 et Aff-
MGGP1 présentent des dispersions importantes par rapport aux résultats expérimentaux en
comparaison avec les modèles de la résistance mécanique développés. Ce constat est confirmé
par des coefficients de corrélation (R) relativement bas de l’ordre de 0,5 à 0,55. Cependant, en
se basant sur l’aspect pratique des modèles, on remarque que les deux modèles Aff-PG1 et
Aff-MGGP1 présentent des erreurs absolues moyennes MAE respectivement équivalent à 20
mm et 22 mm, ce qui à l’échelle réelle est très acceptable, sachant que l’essai de
l’affaissement au cône d’Abraham sur un même béton donne des erreurs type similaires.
Dans la suite, les modèles développées seront utilisés comme fonctions-objectifs, dans le
processus de l’optimisation des dosages par AG pour formuler des bétons selon un cahier des
charges bien défini.
138
Chapitre 4 : Apprentissage automatique et optimisation multiobjectifs
Dans cette première partie, pour évaluer la performance des individus dans l’optimisation
par AG, nous avons utilisé les modèles développés par la PG (Rc-PG1 et Aff -PG1).
Pour limiter la recherche à un domaine, les intervalles de variation des inputs ont été
également définis (Tableau 43), Ceci permet aussi de limiter certains ingrédients à des
niveaux raisonnables d’utilisation. Il est à noté que le rapport S/G est le rapport massique du
sable sur gravier et le dosage de cendre volante (Cv) exprimé en pourcentage du dosage en
ciment.
139
Chapitre 4 : Apprentissage automatique et optimisation multiobjectifs
D’après le Tableau 43, si chacun des 6 inputs prend 10 valeurs possibles, à la fin on
obtient 10n cas soit 1 Million de formulations possibles. Ceci est une explosion combinatoire
ce qui confirme l’intérêt d’utiliser les AG pour ce type de problème.
Nous insistons également sur les intervalles de variation des inputs retenus pour ce cas. En
effet l’augmentation des intervalles va forcément agrandir l’espace de recherche de
l’algorithme, par conséquent les résultats de sortie (Résistance et affaissement) peuvent
changer. Inversement ceci peut être un avantage pour l’utilisateur qui peut réduire par
exemple l’intervalle du dosage en ciment s’il veut obtenir un BHP entré de gamme en matière
de résistance à moindre coût. Ces aspects de l’optimisation multiobjectifs par AG, sont traités
en détail dans la dernière section du présent chapitre.
Comme mentionné dans ce travail, la démarche mise en œuvre permet d’identifier une
large gamme de formulation répondant au cahier des charges. La notion d’optimum de
Pareto est utilisée dans ce type de problèmes d’optimisation multi objectifs. Les AG font
converger systématiquement les solutions vers la frontière de Pareto comme on peut le voir
sur la Figure 67. Il est important de noter que dans ce processus d’optimisation nous avons
maximisé simultanément la résistance à la compression et l’affaissement.
140
Chapitre 4 : Apprentissage automatique et optimisation multiobjectifs
Frontière de Pareto
170
169 Formulations
Affaissement (mm) optimales
168
167
166
165
164
163
162
161
160
35 40 45 50 55 60
Rc28 (MPa)
Sur le Tableau 45 quelques compositions des formulations optimales obtenues par les AG
sont récapitulées:
Inputs Outputs
N° E/L E S/G CV AE SP Rc28 Affaissement
(%) (kg/m3) (%) (%) (kg/m3) (kg/m3) (MPa) (mm)
1 36,4 141,4 50,0 13,4 0,04 3,33 53,9 161,9
2 36,2 141,7 49,2 13,4 0,04 3,15 52,1 163,0
3 35,0 141,1 47,8 12,6 0,05 3,01 50,2 164,3
4 35,7 141,8 47,2 13,2 0,05 2,96 49,2 164,9
5 35,3 141,7 46,3 13,0 0,05 2,86 47,7 165,8
6 34,3 141,4 44,7 12,4 0,05 2,69 45,7 166,8
7 33,1 140,7 43,5 11,6 0,06 2,56 43,4 167,7
8 33,4 141,2 42,8 11,9 0,06 2,50 42,1 168,4
9 33,3 141,1 42,6 11,9 0,06 2,48 41,6 168,6
Tableau 45: Quelques solutions optimales de la frontière de Pareto.
On peut voir que le dosage des différents ingrédients est conforme à des dosages
rencontrés pour des BHP. Nous retenons également que les formulations obtenues satisfont
aux critères exigés en matière de résistance mécanique et d’affaissement.
141
Chapitre 4 : Apprentissage automatique et optimisation multiobjectifs
Les matériaux utilisés ainsi que leurs proportions massiques sont présentés sur le Tableau
46. Les matériaux utilisés sont un ciment Portland type I 52,5, une cendres volante (calcique),
un superplastifiant (CHRYSO Fluid Premia 224), des gros granulats concassés (4/20mm), des
granulats concassés (Sable 0/4) et un entraineur d’air (Chryso Air G 100).
Quantité (kg/m3)
Matériau Formulation Formulation Formulation
F1 F4 F9
Eau 141,4 141,8 141,1
Ciment Portland type I 52,5 364,1 372,6 400,3
Cendres volantes 48,6 49,1 47,5
Gros granulats 426,8 435,4 443,9
Granulats fins 217,2 205,7 189,3
Superplastifiant 3,33 2,96 2,48
Entraineur d’air 0,04 0,05 0,06
Tableau 46: Dosage des composants des formulations testées.
Après confection, les éprouvettes sont conservées dans l’eau à température de 20 °C ± 2°C
pendant 28 jours. Les résultats de résistance à la compression sur éprouvettes cylindriques et
affaissements sont reportés dans le Tableau 47.
142
Chapitre 4 : Apprentissage automatique et optimisation multiobjectifs
D’après les résultats du Tableau 47, on remarque une très bonne cohérence entre les
résultats expérimentaux et les prédictions, pour l’affaissement et pour la résistance
mécanique. Ce résultat reflète clairement la bonne prédiction des modèles obtenus par PG.
L’utilisation des matériaux de nature assez similaires de ceux utilisés dans l’élaboration de la
base de données a certainement contribué à la qualité de ces prédictions.
143
Chapitre 4 : Apprentissage automatique et optimisation multiobjectifs
Cependant, devant la variété des formes des modèles prédictifs obtenus par la PG et la
MGGP, il est légitime de se poser la question, sur l’influence de la forme des modèles
prédictifs, sur les résultats de l’optimisation bi objectifs, autrement dit sur la frontière de
Pareto. Cette question est d’autant plus pertinente que deux apprentissages, sur les mêmes
données et le même paramétrage, donneront des modèles de formes différentes et parfois avec
des performances différentes.
Pour répondre à cette question, nous avons procédé à plusieurs optimisations bi-objectif,
en fixant les mêmes intervalles de variation des inputs et en combinant plusieurs modèles,
comme présenté sur le Tableau 48. L’objectif de cette étude étant de regarder l’impact que
cela peut avoir sur l’évolution de la frontière de Pareto.
Propriétés
Model Model
Frontière Pareto Résistance mécanique Affaissement
(Résultats) à 28 Jours (MPa) (mm)
P1 Rc-MGGP1 Aff-PG1
P1s Rc-MGGP1 Aff-PG1
P2 Rc-MGGP2 Aff-PG1
P3 Rc-MGGP1 Aff-MGGP1
P4 Rc-REG Aff-PG1
Tableau 48 : Combinaisons des différents modèles.
L’objectif des choix du Tableau 48 permettra d’évaluer l’impact des modèles développés
par des approches différentes (PG et MGGP) et l’impact des modèles développés
successivement avec la même approche (MGGP1 et MGGP2). Dans la dernière ligne du
tableau, nous intégrons la comparaison avec un modèle proposé par l’auteur de la base de
données. Ce modèle a été développé par régression (Lim et al., 2004).
144
Chapitre 4 : Apprentissage automatique et optimisation multiobjectifs
Dans la partie optimisation, le paramétrage des AG et les intervalles des variations des
inputs sont identiques, respectivement présentés au Tableau 49 et Tableau 50.
Les résultats obtenus en termes de frontière de Pareto sont présentés dans la Figure 68:
Figure 68 : Les différents résultats d’optimisation obtenus par les différents modèles.
145
Chapitre 4 : Apprentissage automatique et optimisation multiobjectifs
les mêmes modèles (Résultats P1 et P1s). Il est important de signaler que ce résultat est
obtenu pour un paramétrage identique de l’AG. Nous remarquons sur la même figure que la
frontière n’est certes pas la même mais la différence entre les deux frontière P1 et P1s est
imperceptible voire inexistante. Ceci confirme la très bonne répétabilité des résultats de
l’optimisation des AG.
La seconde remarque porte sur la similitude des résultats P1 et P2 où les modèles utilisés
sont développés par le même outil (La MGGP). Toutefois, la forme des modèles est
différente. Ceci permet de conclure que la différence d’écriture des modèles n’a pas
d’influence sur le résultat de l’optimisation, à condition que les modèles soient de
performances similaires. Il est important de signaler à ce stade que la performance des
modèles est jugée vis-à-vis des critères mentionnés dans la première partie de ce chapitre à
savoir la fitness, le coefficient de corrélation, la moyenne des erreurs absolues et l’erreur
quadratique moyenne. Ceci a été davantage confirmé par la courbe P4, obtenu par un modèle
de résistance mécanique obtenu par régression multiple (Lim et al., 2004) où on remarque que
les résultats sont aussi proches de la frontière P1 et P2.
Inputs Outputs
N° E/L E S/G CV AE SP Rc28 Affaissement
(%) (kg/m3) (%) (%) (kg/m3) (kg/m3) (MPa) (mm)
1 30,89 160,04 52,35 0,54 0,08 7,07 77,27 198,93
2 30,90 161,36 52,35 0,70 0,08 7,40 76,09 200,92
3 30,89 161,59 52,36 0,72 0,08 8,14 75,82 201,78
4 30,90 164,16 52,36 0,79 0,08 8,02 73,65 204,72
5 30,82 166,82 52,36 0,65 0,08 7,84 71,59 207,37
6 30,92 169,28 52,35 0,72 0,08 8,03 69,36 210,44
7 30,90 171,11 52,36 0,68 0,08 8,08 67,88 212,55
8 30,92 173,96 52,38 0,78 0,08 7,97 65,44 215,65
9 30,92 176,41 52,37 0,69 0,08 7,97 63,43 218,25
10 30,93 178,79 52,40 0,71 0,08 8,11 61,49 221,08
Tableau 51 : Quelques solutions optimales de la frontière de Pareto (P3).
Les matériaux utilisés ainsi que leurs proportions massiques sont présentés sur le Tableau
52. Il est à rappeler que les caractéristiques de ces matériaux sont indiquées au (Tableau 37)
146
Chapitre 4 : Apprentissage automatique et optimisation multiobjectifs
Quantité (kg/m3)
Matériau Formulation Formulation
F1 F10
Eau 160,0 178,8
Ciment Portland type I 52,5 516,7 576,0
Cendres volantes 2,8 4,1
Gros granulats 405,5 392,8
Granulats fins 212,3 205,8
Superplastifiant 7,07 8,11
Entraineur d’air 0,08 0,08
Tableau 52: Dosage des composants des formulations testées.
Après confection, les éprouvettes sont conservées dans l’eau à température de 20 °C ± 2°C
pendant 28 jours. Les résultats de résistance à la compression sur éprouvettes cylindriques et
affaissements sont reportés dans le Tableau 53.
Nous constatons à travers les résultats du Tableau 53, que les résultats expérimentaux
concordent avec les résultats de prédictions notamment pour la résistance mécanique où les
deux formulations testées donnent des résultats très proches de la prédiction du modèle.
Cependant, le modèle prédictif de l’affaissement assure une bonne prédiction une fois sur
deux, ce résultats est du aux performances modestes du modèle de l’affaissement reflétés par
un R et un MAE respectivement de 0,5 et 22, contrairement au modèle de la résistance
mécanique d’un R et MAE respectivement de 0,97 et 1,91.
147
Chapitre 4 : Apprentissage automatique et optimisation multiobjectifs
Conclusion du chapitre
Dans ce chapitre, pour un cas pratique (la formulation d’un BHP) la démarche globale de
formulation a été mise en œuvre et discutée puis validée à différents niveaux.
En termes d’identification des modèles prédictifs comme pour le chapitre qui traite la
compacité granulaire, les performances des algorithmes de type PG ou MGGP ont été
comparées. Il apparait comme pour le chapitre 3, l’avantage d’une méthode sur l’autre n’est
pas flagrant pour les cas étudiés.
En terme d’optimisation des formulations par les AG afin d’identifier une gamme de
béton qui répond à un cahier des charges spécifique, la démarche a permis d’identifier une
gamme de béton qui répond aux critères imposés et dont les formulations sont conformes aux
formulations des BHP.
Dans le chapitre suivant, nous étudierons l’extension de la méthode hybride (PG &AG)
proposée dans ce travail sur un cas d’optimisation à cinq objectifs. Le cas d’étude est un
mortier à base de sable de fonderie, dont les objectifs à optimiser dans le cahier des charges,
sont la résistance mécanique à la compression et à la flexion, la maniabilité, la compacité
granulaire et le coût.
148
CHAPITRE 5 : Formulation de mortiers
à base de matériaux alternatifs
Chapitre 5 : Formulation des mortiers à base de matériaux alternatifs
Introduction du chapitre
Dans l’optique de capitaliser tous les concepts et les résultats obtenus dans les parties
précédentes, le présent chapitre a pour objectif d’appliquer cette nouvelle approche de
formulation des mélanges sur le cas de mortiers à base de matériaux alternatifs. Ce choix a été
fait, d’une part pour développer une base de données représentative avec les matériaux
disponibles. D’autre part faire autant que possible des validations à différentes étapes. Nous
avons également veillé dans ce travail à donner le maximum d’informations sur les
ingrédients et les méthodes pour permettre une ré-exploitation de cette base de données dans
d’autres travaux. Ce travail permettra aussi d’éclaircir davantage la démarche globale dans
son intégralité allant de la base de données expérimentale qui représente la matière première
pour la PG dans la modélisation des propriétés des mortiers, passant par la formulation
multicritères des mortiers par AG pour enfin finir avec la validation au laboratoire des
résultats obtenus. Il est à noter que la base de données a été développée par l’utilisation d’un
plan d’expériences optimisé.
L’étude des différents aspects de l’optimisation multiobjectifs abordée dans le chapitre 4,
notamment l’étude de l’influence des formes des modèles prédictifs sur les résultats de
l’optimisation et la répétabilité des AG, nous a permis d’envisager également d’optimiser les
mélanges de mortier selon 5 objectifs, contrairement au cas des BHP abordé jusqu’ici (Bi-
objectif). En effet, l’optimisation multiobjectifs de ces mortiers prend comme critère à
optimiser la compacité du squelette granulaire, la résistance mécanique à la compression, la
résistance mécanique à la flexion, la maniabilité et le coût à titre indicatif. Ceci permet de
mettre en exergue l’intérêt de la méthode proposée, dans la prise en considération de plusieurs
facteurs simultanément dans le processus de formulation.
Le présent cas d’étude est un mortier à base de sable de fonderie. Ce matériau alternatif a
fait partie des matériaux d’apprentissage du modèle de la compacité développé dans le
chapitre 3, ce choix été fait dans l’optique d’assurer la fiabilité du modèle comme critère
d’optimisation dans ce cas d’étude. Il est important de souligner qu’on a utilisé la fraction
inférieure à 0,2 mm du sable de fonderie, ceci permettra d’adapter cette base de données sur
d’autres matériaux alternatifs comme les sédiments de dragage fluviaux et marins qui ont des
Dmax similaires en granulométrie.
150
Chapitre 5 : Formulation des mortiers à base de matériaux alternatifs
Pour étudier l’impact des dosages des composants (Facteurs) sur une propriété quelconque
du mélange (Réponse), il est plus logique au premier abord de réaliser toutes les combinaisons
possibles. Ceci est nommé dans la discipline des plans d’expériences, un plan factoriel
complet. Ce type de plan permettra sans doute, de réaliser des bases de données complètes
permettant d’extraire toutes les informations nécessaires du phénomène étudié et par
conséquent la possibilité de proposer un modèle prédictif prenant en compte l’importance de
chaque facteur et de leurs interactions.
Cependant, il n’est pas toujours évident de réaliser un plan d’expérience factoriel complet.
En effet, avec l’augmentation rapide du coût des expériences, il est essentiel que le
phénomène soit bien étudié avec le moins d’expériences. Ceci est une des raisons importantes
pour opter à un plan d’expérience optimisé.
Différentes méthodes de plans d’expériences ont été développées et améliorées. Le plus
connu est le plan de Taguchi qui porte le nom de son inventeur.
Pour notre cas d’étude qui porte sur la formulation de mortier à base de coproduits. Nous
avons réalisé un plan d’expérience optimisé développé par le Logiciel XLSTAT. Ce dernier
intègre la méthode des matrices factorielles fractionnaires. Il est à noter que ce logiciel a
permis, par rapport à un plan complet de deviser le nombre d’expériences par 2.
Détermination des facteurs
Avant d’aborder le plan d’expériences, il était nécessaire d’identifier une formulation de
base qui nous permet par la suite d’étudier les variations nécessaires pour cerner le champ des
formulations sélectionnées.
Le mortier formulé est composé de ciment, d’eau, de sable et de matériau alternatif. Le
matériau alternatif de substitution est un sable de fonderie inférieur à 0,2 mm.
Afin d’étudier l’influence des composants sur les propriétés du mortier, nous avons retenu
les plages des variations suivantes (Tableau 54) :
151
Chapitre 5 : Formulation des mortiers à base de matériaux alternatifs
Sur la base du tableau, nous pouvons évaluer le nombre d’expériences pour un plan
factoriel complet égal à 56 expériences.
Plan d’expériences optimisé
Sur la Figure 69 est présentée les expériences d’un plan factoriel complet, sur l’espace de
recherche défini (C, E/C et F/S). Si nous répétons chaque expérience trois fois, le plan
d’expériences atteint 168 expériences pour un plan factoriel complet.
Figure 69: Répartition des expériences dans un plan d’expériences complet pour notre
cas d’étude.
152
Chapitre 5 : Formulation des mortiers à base de matériaux alternatifs
La traduction du plan d’expérience optimisé est reproduite dans le Tableau 55. Le plan est
développé d’expériences témoins au nombre de 8. Ceci traduit un plan factoriel complet pour
le cas d’un mortier sans matériaux alternatifs, et 24 expériences correspondent à un plan
d’expériences demi complet pour l’utilisation d’un matériau alternatif.
Pour valider les résultats et s’assurer de la répétabilité de ces derniers, chaque expérience
est répétée trois fois. Le plan d’expérience permet de générer ainsi une base de données totale
de 96 enregistrements.
153
Chapitre 5 : Formulation des mortiers à base de matériaux alternatifs
154
Chapitre 5 : Formulation des mortiers à base de matériaux alternatifs
C (gr) E/C (%) F/S (%) S (gr) Rc28 (MPa) Rf28 (MPa) Maniabilité (s)
450 40 0 1350 40,43 7,96 55
450 60 0 1350 43,07 9,07 3
450 70 0 1350 29,67 8,98 1
400 50 0 1350 55,03 9,2 10
400 60 0 1350 42,34 8,35 2
400 70 0 1350 38,99 7 1
450 40 5 1350 14,1 2,61 52
450 50 10 1350 49,3 10,15 21
400 60 15 1350 41,93 8,76 11
450 70 20 1350 35,71 8,71 2
400 40 30 1350 1,87 1,52 122
450 60 40 1350 30,44 7,34 126
Tableau 56 : Extrait de la base de données réalisée.
A ce stade du travail, nous n’avons pas fait varier le type de ciment ni l’utilisation
d’adjuvant pour les raisons évoquées ci-dessus à savoir l’explosion du nombre d’expériences
et le temps nécessaire pour réaliser ce travail.
155
Chapitre 5 : Formulation des mortiers à base de matériaux alternatifs
Avec une démarche similaire à celle développée dans le chapitre 4, nous procédons à la
modélisation des propriétés du mortier à base de sable fin de fonderie par apprentissage
automatique, en utilisant la MGGP. L’apprentissage s’effectue sur la base de données
d’apprentissage réalisée ci-dessus. La modélisation intègre la résistance mécanique à la
compression à 28 jours, la résistance mécanique à la flexion à 28 jours et la maniabilité.
Dans cette section, nous présentons le paramétrage de la MGGP appliqué pour chaque cas
de modélisation ainsi que les modèles obtenus pour chaque propriété ainsi que les
performances de chaque modèle. Il est à rappeler que la performance « fitness » de chaque
modèle représente la somme des erreurs absolues, ceci explique pourquoi la fitness de
terminaison de recherche dans le paramétrage de départ est fonction de la propriété modélisée.
Notons que dans cette partie, nous ne reviendrons pas sur l’importance du paramétrage et la
pertinence des différents modèles générés.
156
Chapitre 5 : Formulation des mortiers à base de matériaux alternatifs
∗ ∗( ) ∗ ( ) ∗ ( )∗ ( )
∗( )∗( )
( ) ( )
Équation 18
Où : C : Dosage en ciment (gr), E/C : Rapport eau/Ciment (%), F/S: Rapport massique du sable de fonderie
sur sable (%), S: Dosage en sable (gr)
Performances du modèle
Le modèle développé par la MGGP est évalué en se basant sur le coefficient de corrélation
(R), l’Erreur Absolue Moyenne (MAE) et la Racine de l’Erreur Quadratique Moyenne
(RMSE) comme présenté sur le Tableau 58 et sur la Figure 71 :
Critère Valeur
R 0,86
MAE 5,77
RMSE 9,06
Tableau 58 : Performances du modèle de la résistance mécanique à la compression.
50
40
30
20
10
0
0 10 20 30 40 50 60 70
Mesure expérimentale
Nous remarquons sur cette Figure 71, la concentration des résultats autour de la courbe
optimale (X=Y), reflétant ainsi la grande concordance entre les résultats de prédiction avec les
résultats expérimentaux. Ceci est confirmé par les critères d’évaluation statistique notamment
157
Chapitre 5 : Formulation des mortiers à base de matériaux alternatifs
le coefficient de corrélation égal à 0,86 et par l’erreur absolue moyenne égale à 5,7 MPa. Ce
modèle sera utilisé comme fonction objectif dans la formulation des mortiers par AG.
∗( ) ∗( ) ∗ ∗( ) ∗ ∗ ∗
∗( ) ∗( )
Équation 19
Où : C : Dosage en ciment (gr), E/C : Rapport eau/Ciment (%), F/S: Rapport massique du sable de fonderie
sur sable (%), S: Dosage en sable (gr).
Sur les mêmes critères utilisés précédemment, le modèle développé par la MGGP est
évalué en se basant sur le coefficient de corrélation (R), l’Erreur Absolue Moyenne (MAE) et
158
Chapitre 5 : Formulation des mortiers à base de matériaux alternatifs
12
Prédictions
0
0 3 6 9 12 15
Mesure expérimentale
D’après les résultats obtenus, le modèle développé est performant dans la prédiction de la
résistance mécanique à la flexion du mortier à base de sable de fonderie, au même titre que le
modèle de la résistance à la compression, avec un coefficient de corrélation similaire
équivalent à 0,86, quant à l’erreur absolue moyenne, elle n’excède pas les 2MPa. Ce modèle
représente la seconde fonction-objectif dans la démarche d’optimisation multi critère par les
AG.
159
Chapitre 5 : Formulation des mortiers à base de matériaux alternatifs
Le modèle obtenu de la maniabilité du mortier à base de sable de fonderie est présenté par
l’Équation 20 :
∗
( ∗ ) ( ∗ ( )) ( ∗ ∗ ) ( ∗ ∗ )
Équation 20
Où : C : Dosage en ciment (gr), E/C : Rapport eau/Ciment (%), F/S: Rapport massique du sable de fonderie
sur sable (%), S: Dosage en sable (gr).
160
Chapitre 5 : Formulation des mortiers à base de matériaux alternatifs
Le modèle développé par la MGGP est évalué en se basant sur les mêmes critères
d’évaluation utilisés précédemment. Le Tableau 62 et la Figure 73 récapitulent les résultats
obtenus :
Critères Valeur
R 0,94
MAE 13,9
RMSE 22,4
Tableau 62 : Performances du modèle de la maniabilité.
250
Prédictions
200
150
100
50
0
0 50 100 150 200 250 300
Mesure expérimentale
Se basant sur la forte corrélation entre les valeurs prédites et valeurs expérimentales
traduite par un coefficient de corrélation R=0,94 et la valeur moyenne des erreurs absolues, le
présent modèle de la maniabilité du mortier développé par MGGP est retenu comme fonction
d’évaluation dans le processus d’optimisation multiobjectifs.
161
Chapitre 5 : Formulation des mortiers à base de matériaux alternatifs
( ∗ ) ( ∗ ) ( ∗( )) ( ∗ (( ) ))
∗
( ∗ ( ( ) )) ( ∗ )
Équation 12
Équation 21
P1, P2, P3 et P4: Sont les prix en centime d’euro/gr, respectifs du ciment, eau, sable fin de fonderie et sable.
C, E, F et S: Sont les dosages respectifs du ciment, eau, sable fin de fonderie et sable.
162
Chapitre 5 : Formulation des mortiers à base de matériaux alternatifs
Le développement des modèles prédictifs des différentes propriétés des mortiers est
l’étape la plus sensible, puisque les résultats de l’optimisation par AG dépendent étroitement
de la qualité des fonctions d’évaluation utilisées. En effet, les AG vont procéder à la création
des premières formulations. Chaque formulation qui représente un individu sera évaluée via
les fonctions d’évaluation, ensuite les individus vont être croisés de génération en génération
pour trouver de nouvelles formulations plus intéressantes en matière de propriétés mécaniques
et de coût. L’ensemble de solutions optimales représentent le résultat final des AG. La qualité
de prédiction des modèles développés permet d’envisager l’optimisation par les AG
Comme développé précédemment, il convient dans la seconde étape de paramétrer les AG
pour trouver rapidement les meilleures formulations, ceci implique la définition de l’espace de
recherche de l’algorithme, de définir le nombre de génération et surtout le sens de
l’optimisation de chaque fonction (Maximiser ou Minimiser). Le Tableau 63 récapitule le
paramétrage retenu pour la formulation d’un mortier à base de sable de fonderie.
163
Chapitre 5 : Formulation des mortiers à base de matériaux alternatifs
Dans le choix de l’espace de recherche résumé dans le Tableau 65, nous avons maintenu
le même espace d’apprentissage qui représente l’intervalle de confiance des modèles
développés.
C (gr) E/C(%) F/S (%) Sable (gr)
Min 400,00 40 0,000 1350,00
Max 450,00 70 40,000 1350,00
Tableau 65 : Intervalles des variations des inputs (Dosages potentiel des constituants).
Les prédictions de l’optimisation effectuée par les AG sont présentées dans le Tableau 66.
Inputs Outputs
N C (gr) E/C(%) F/S(%) Sable (gr) Rc28(MPa) Rf28(MPa) Man(s) Comp Prix (Centime)
1 434 69,9 4 1350 58,4 8,4 1 0,594 6,48
2 445 59,9 0 1350 50,2 9,6 2 0,551 7,00
3 442 57,4 1 1350 38,8 9,5 3 0,574 7,10
4 427 68,6 6 1350 55,6 8,5 1 0,600 6,47
5 419 59,1 7 1350 46,1 8,8 5 0,604 6,79
6 441 53,6 0 1350 53,7 9,2 3 0,556 7,32
7 430 69,9 6 1350 57,0 8,3 1 0,599 6,45
8 409 59,3 6 1350 46,3 8,8 5 0,602 6,69
9 421 63,1 2 1350 43,9 9,2 2 0,578 6,63
10 442 55,3 0 1350 49,6 9,4 3 0,563 7,22
11 441 56,1 1 1350 42,3 9,4 3 0,569 7,17
12 404 70,0 4 1350 56,2 8,0 1 0,591 6,22
13 434 58,3 0 1350 48,8 9,5 2 0,547 6,98
14 408 60,0 4 1350 48,5 9,0 3 0,593 6,64
15 408 60,4 6 1350 47,8 8,9 4 0,600 6,63
16 419 59,1 6 1350 47,3 8,9 4 0,601 6,79
17 430 63,2 1 1350 44,0 9,3 2 0,571 6,71
18 405 67,1 2 1350 49,6 8,6 1 0,582 6,33
19 414 57,8 3 1350 45,1 9,1 4 0,589 6,81
20 404 69,9 4 1350 56,3 8,0 1 0,594 6,23
21 423 58,6 7 1350 46,2 8,8 5 0,604 6,85
Tableau 66 : Formulations optimales du mortier à base de sable de fonderie obtenues par
AG.
Les résultats présentés, dans le Tableau 66 et la Figure 74, représentent l’ensemble des
solutions optimales (solution non dominées) obtenu par les AG. Chaque formulation, qui
inclut les quatre inputs, est évaluée selon les cinq critères du cahier des charges (résistance à
la compression, résistance à la flexion, maniabilité, compacité et le coût).
164
Chapitre 5 : Formulation des mortiers à base de matériaux alternatifs
Il apparait d’après ces résultats que l’ensemble des solutions optimales satisfont au cahier
des charges. Nous pouvons remarquer également qu’en matière de résistance mécanique à la
compression, la différence entre les formulations 1 et la formulation 21 est de 12 MPa, ceci
est intéressant d’un point de vue compromis, car l’utilisateur aura une large gamme de
mortier.
Nous remarquons également que les résistances à la flexion demeurent similaires pour
l’ensemble des formulations. Mais nous attirons surtout l’attention sur le taux d’incorporation
des sables de fonderie qui est plafonnée à 8 % dans toutes les solutions optimales, ceci est du
principalement au dosage fixe du sable standard. En effet, le volume des inclusions granulaire
est un paramètre de forte influence sur la résistance mécanique (Yammine et al., 2008), au-
delà d’un volume optimal, la résistance mécanique chute considérablement, c’est pourquoi
l’ajout supplémentaire du sable de fonderie au-delà de 8% n’a pas été possible. A cet effet, il
sera intéressant d’exploiter un intervalle de variation du sable inférieur à 1350 gr où
d’explorer des dosages plus importants en ciment afin d’offrir la possibilité d’augmenter le
volume de la pâte dans les formulations.
165
Chapitre 5 : Formulation des mortiers à base de matériaux alternatifs
8%. Les résultats de la résistance à la compression et à la flexion ainsi que la maniabilité, sont
récapitulés dans le Tableau 67.
Dosage Man Man Rc28 Rc28 Rf28 Rf28
(s) (s) (MPa) (MPa) (MPa) (MPa)
Formulation C E Fines Sable
PRED EXP PRED EXP PRED EXP
N° (gr) (gr) (%) (gr)
10 419 59 7 1350 5 5 46,1 42,5 8,8 9,4
20 442 55 0 1350 3 2 49,6 45,0 9,3 11,7
5 404 70 4 1350 1 3 56,2 47,5 8,0 9,4
Coefficient de corrélation R 0,65 0,98 0,79
Moyenne des erreurs absolue MAE 1 5,6 1,46
Tableau 67 : Résultats obtenus pour les formulations testées au laboratoire.
Au vue des résultats obtenus pour les trois formulations testées, nous confirmons la
grande concordance entre les résultats expérimentaux et les résultats de prédiction pour les
trois propriétés : la maniabilité ainsi que les résistances mécaniques à la compression et à la
flexion. Ceci reflète la grande précision des modèles prédictifs développés par apprentissage
automatique, donnant des moyennes des erreurs absolues équivalentes à 1 seconde pour la
maniabilité, de 5,6 MPa pour la résistance mécanique à la compression et 1,46 MPa pour la
résistance mécanique à la flexion. Ces erreurs restent acceptables pour le cas d’un mortier.
La méthode de formulation proposée dans ce travail est très fiable dans ce cas d’étude
d’optimisation à cinq objectifs. La qualité de ces résultats est due principalement aux
intervalles de variation des inputs qui correspondaient au domaine d’apprentissage. C’est
pourquoi il est intéressant de prospecter la fiabilité des modèles en extrapolation.
166
Chapitre 5 : Formulation des mortiers à base de matériaux alternatifs
Les modèles prédictifs développés sur les mortiers à base de sable de fonderie se sont
révélé des modèles de grande capacité de prédictions en interpolation sur le même type de
mortier. Il est néanmoins plus judicieux de tester la robustesse des modèles pour des mortiers
à base de matériaux alternatifs autres que les sables de fonderie, c’est pourquoi dans cette
partie, nous utilisons les même modèles prédictifs, dans la formulation de mortiers à base de
sable de dragage en plus des matériaux standards à savoir : le ciment 52,5 N, de sable
normalisé et d’eau.
Il est à noter que le paramétrage des AG et les intervalles de variation des inputs, sont
identiques au premier cas d’étude de mortier à base de sable de fonderie comme décrit
respectivement dans le Tableau 63 et Tableau 65.
167
Chapitre 5 : Formulation des mortiers à base de matériaux alternatifs
Au vu de ces résultats, nous remarquons que les intervalles de variation des inputs varient
selon les plages définies au préalable à l’exception du dosage en fines qui est plafonné à 5%.
Ceci nous amène à tirer la même conclusion que le précédent cas d’étude de mortier à base
de sable de fonderie, sur le caractère compromettant du dosage fixe en sable standards égale à
1350 gr qui empêche une incorporation plus importante de matériaux alternatif.
Il est cependant important de valider au laboratoire quelques formulations afin d’étudier la
capacité de prédiction des modèles développés précédemment, sur ce nouveau matériau
alternatif qui n’a pas été utilisé dans la base de données d’apprentissage.
Trois formulations ont été choisies, se basant sur le taux d’incorporation du matériau
alternatif. La « Formulation 2 » avec le maximum de dosage en sable de dragage proche de
5% et la « Formulation 10 » sans sable de dragage incorporé et une dernière formulation
« Formulation 16 » de résistance mécanique à la compression la plus faible. Les résultats de
la résistance à la compression et à la flexion ainsi que la maniabilité, sont récapitulés au
Tableau 69.
Dosage Man Man Rc28 Rc28 Rf28 Rf28
(s) (s) (MPa) (MPa) (MPa) (MPa)
Formulation C E Fines Sable
PRED EXP PRED EXP PRED EXP
N° (gr) (gr) (gr) (gr)
2 435 69,9 5 1350 0 1 58,3 48,8 8,4 7,0
10 450 53,2 0 1350 3 5 53,7 45,0 9,2 9,4
16 441 58,5 1 1350 2 4 44,6 43,1 9,5 9,4
Coefficient de corrélation R 0,99 0,93 0,96
Moyenne des erreurs absolue MAE 2 6,6 0,56
Tableau 69 : Résultats obtenus pour les formulations testées au laboratoire.
168
Chapitre 5 : Formulation des mortiers à base de matériaux alternatifs
Les résultats obtenus pour ce dernier cas d’étude sont très concluants. En effet, les erreurs
absolues moyennes obtenues sont très acceptables au même titre que le mortier à base de
sable de fonderie, avec une erreur type de 2 secondes pour la maniabilité et 6,6 et 0,56 MPa
respectivement pour la résistance mécanique à la compression et à la flexion.
Ces résultats obtenus confirment davantage la capacité de prédiction des modèles en
interpolation, car les intervalles de variation dans ce cas d’étude correspondent bien à l’espace
d’apprentissage des modèles.
Au vue de ces résultats, nous concluons sur la possibilité d’établir une classification
physico-chimique des matériaux alternatifs dont les comportements sont régis les mêmes
modèles prédictifs développés sur des matériaux similaires.
169
Chapitre 5 : Formulation des mortiers à base de matériaux alternatifs
Conclusion du chapitre
Dans ce chapitre, nous avons appliqué la méthode de formulation des mortiers à base de
matériaux alternatifs, développée durant cette thèse. A travers ce cas d’étude d’un mortier à
base de sable de fonderie, nous avons procédé à la modélisation des propriétés de ce dernier
par apprentissage automatique via la MGGP, ensuite nous avons utilisé ces modèles pour
l’optimisation multiobjectifs via les AG.
Nous avons utilisé également le modèle prédictif de la compacité granulaire adapté aux
matériaux alternatifs développé précédemment, comme un des critères d’optimisation. Ceci
permet d’une part de prendre un paramètre important dans le processus d’optimisation qu’est
la compacité granulaire et de prendre indirectement les spécificités des déchets, pris en
compte dans le modèle.
Les résultats obtenus sont très satisfaisants, en se référents à la bonne concordance entre
les résultats prédits et les résultats obtenus expérimentalement, ceci implique la résistance
mécanique à la compression, à la flexion ou la maniabilité des mortiers à base de sable de
fonderie. Ceci nous permettra de conclure que la présente méthode de formulation développée
est adaptée à l’utilisation des matériaux alternatifs.
L’intérêt de cette méthode se confirme d’avantage dans la possibilité de proposer
plusieurs solutions optimales en optimisant simultanément plusieurs objectifs conflictuels.
Dans le présent cas d’étude, nous sommes arrivés jusqu’à cinq propriétés à optimiser en
incluant le coût des formulations.
Nous avons testé les modèles développés sur les mortiers à base de sable de fonderie, sur
des mortiers à base de sable de dragage et les résultats sont de même qualité que le cas
d’apprentissage sous réserve d’appliquer l’optimisation dans le même espace d’apprentissage.
Ceci permet de réduire les cas d’apprentissage pour différents matériaux alternatifs de
propriétés similaires.
170
Conclusion générale et perspectives
Concrètement, l’outil auquel nous aspirions devait être capable de tester numériquement
plusieurs formulations à base de matériaux alternatifs et proposer les meilleures formulations
qui satisfont simultanément à plusieurs critères fixés par le cahier de charge. Ceci parait très
complexe à l’état de connaissances actuel mais l’exploitation des nouveaux algorithmes issus
de l’intelligence artificielle nous a permis de concrétiser cette démarche dans une méthode de
formulation hybride présentée dans ce travail.
Le travail de recherche présenté dans ce manuscrit de thèse, est organisé en cinq chapitres
qui sont organisé comme suit :
Dans le premier chapitre, nous avons rappelé le contexte environnemental mondial actuel
et l’évolution de la règlementation concernant la gestion des déchets, avec la mise en exergue
de l’importance d’encourager la valorisation des déchets dans une démarche d’écologie
171
Conclusion générale et perspectives
circulaire. Nous avons par la suite rappelé leurs natures et les travaux scientifiques relatifs à la
valorisation de chaque déchet dans les différentes applications dans le Génie Civil, en
particulier le béton qui fait l’objet de notre étude.
Le second chapitre quant à lui, est un chapitre de caractérisation des matériaux de l’étude
avec des essais de caractérisation étaient orientés directement vers l’objectif de l’étude.
Dans le troisième chapitre, nous avons procéder à l’application de la Programmation
Génétique dans ces deux versions (PG et MGGP) dans la modélisation par apprentissage
automatique de la compacité granulaire. Ce paramètre a été choisi en premier, en raison de
son importance dans le processus d’optimisation des mélanges des bétons et son influence sur
la plupart des autres propriétés mécaniques.
Nous avons en premier, déterminé les Inputs et les Outputs du système sur la base d’une
étude paramétrique. Après paramétrage de l’algorithme, l’apprentissage est lancé et nous
proposons un modèle prédictif de la compacité d’un mélange granulaire à concentration
donnée de ce mélange, à partir de grandeurs physique des composants aisément accessibles à
l'expérience.
172
Conclusion générale et perspectives
l’affaissement d’un BHP par PG, dont les résultats de validation expérimentale confirment
bien la bonne capacité de prédiction des deux modèles. Ensuite l’optimisation multiobjectifs
par AG, a permis d’obtenir un ensemble de solutions optimales présentées par la frontière de
Pareto, démontrant le grand intérêt de démarche global par le test d’un nombre très importants
de formulation numériquement, qui d’un point de vu pratique, présente un grain financier et
de temps considérables.
Dans la dernière partie de ce chapitre 4, une étude spécifique a été réalisé qui permet
d’étudier la répétabilité des résultats afin de cerner les forces et les faiblesses des outils
utilisés (PG PG & MGGP et AG AG) ainsi que leurs combinaisons. Nous avons conclu que
les résultats obtenus par la PG et la MGGP donnent des modèles de formes différentes à
chaque modélisation. Mais la forme des modèles n’a pas d’impact dans l’optimisation
multiobjectifs, seules les performances des modèles comptent.
Nous avons confirmé dans cette étude que les AG donnent un ensemble de solutions
optimales non dominées, et présentent une très bonne répétabilité, sous réserve de maintenir
un paramétrage identique à chaque processus d’optimisation.
Dans le dernier chapitre (Chapitre5), nous avons capitalisé, les concepts développés sur
l’apprentissage automatique et l’optimisation multiobjectifs, dans un cas pratique qui porte
sur les mortiers à base de matériaux alternatifs. Pour se faire, on a réalisé un plan
d’expériences optimisé de 96 formulations de mortier à base de sable de fonderie. Ce plan
d’expériences a permis de réaliser une base de données qui a servis d’environnement
d’apprentissage pour la MGGP. L’apprentissage a permis d’obtenir des modèles de résistance
mécanique à la compression et à la flexion ainsi que la maniabilité, dont les performances ont
été validées.
En plus du modèle du coût, nous avons intégré le modèle prédictif de la compacité
granulaire adapté aux matériaux alternatifs développé précédemment, comme un des critères
d’optimisation. Ceci a permis d’une part de prendre un paramètre important dans le processus
d’optimisation qu’est la compacité granulaire et de prendre indirectement les spécificités des
déchets, pris en compte dans le modèle.
Les résultats obtenus sont très satisfaisants, en se référant à la bonne concordance entre les
résultats prédits et les résultats obtenus expérimentalement des modèles développés par
MGGP.
173
Conclusion générale et perspectives
Nous avons testé les mêmes modèles développés sur les mortiers à base de sable de
fonderie, sur des mortiers à base de sable de dragage et les résultats sont de qualité similaire
que le cas d’apprentissage, sous réserve d’appliquer l’optimisation dans le même espace
d’apprentissage. Ce dernier cas d’étude a permis de conclure sur la possibilité d’établir une
classification des matériaux alternatifs dont les comportements sont régis par les mêmes
modèles prédictifs développés sur des matériaux similaires.
Les résultats obtenus constituent une base solide pour envisager l'utilisation et le
développement de cette méthode de formulation adaptés aux sous-produits industriels.
Pour continuer d’avancer sur cette voie, les points suivants seront recommandés pour de
futures recherches:
Confronter les résultats de l’optimisation multiobjectifs obtenus par deux types de
fonctions objectifs (modèles physiques et modèles obtenus par apprentissage
automatique)
Capitaliser les données de la littérature scientifique, relatives aux formulations de bétons
et mortiers à base de matériaux alternatifs dans des bases de données et la gestion de ces
dernières par des outils spécifiques de classification et traitement des données.
L’amélioration de la qualité des modèles prédictifs développés par apprentissage
automatique notamment en extrapolation, par l’intégration des paramètres de nature des
constituants, dans les modèles prédictifs. Ceci nécessite des études paramétriques
spécifiques à chaque propriété modélisée, en utilisant des plans d’expériences adéquats.
Etablir une classification des matériaux alternatifs dont l’impact sur les propriétés des
bétons et mortiers est régi par les mêmes modèles de comportement.
La validation du modèle de la compacité granulaire développé dans ce travail, sur
d’autres matériaux alternatifs et son intégration dans une interface graphique simple
d’utilisation.
Appliquer la méthode proposée dans ce travail en technique routière.
L’objectif ultime de ce travail est la concrétisation de cette méthode hybride par la mise en
place d’une interface adaptée pour chaque application du génie civil.
Cette interface connectée à une base de données permettra le perfectionnement des modèles
prédictifs par apprentissage automatique à temps réel sur des bases de données centralisées.
Au second niveau, un logiciel doté de multitudes de modèles prédictifs (propriétés
mécaniques, environnementales et de coûts) offre à l’utilisateur le choix des propriétés à
174
Conclusion générale et perspectives
175
Références bibliographiques
Références bibliographiques
ADEME, 2015. Déchets : Chiffres-clés (No. Réf. 8500), ADEME Éditions 2015.
Al-Shihri, M.A., 2010. Structural Optimization Using a Novel Genetic Algorithm for Rapid
Convergence. International Journal of Civil & Structural Engineering 1, 123–138.
Amar, M., Benzerzour, M., Abriak, N.E., Maherzi, W., 2016. Étude de valorisation des
sédiments de dragage. Presented at the 34èmes Rencontres de l’AUGC du 25 au 27
Mai 2016, Université de Liège.
A.M.O.R.C.E, 2012. Etat des lieux de la gestion des mâchefers en France (Série Technique
No. DT 50).
Bauchard, M., 1984. Utilisation en technique routière de granulats provenant du concassage
de béton de démolition (Bull, liaison Labo P. et Ch. - 134 - No. Réf. 286).
Becquart, F., 2007. 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 (Thèse de doctorat). Université de Lille 1.
Berredjem, L., Arabi, N., Molez, L., 2015. Contribution à l’étude des indicateurs de durabilité
des bétons à base des granulats recyclés, in: 33ÈMES Rencontres Universitaires de
Génie Civil du27 Au 29 Mai 2015. Bayonne-France.
Boughriet, R., n.d. Loi économie bleue : le rejet de boues de dragage polluées est interdit à
compter de 2020 [www document]. Actu-Environnement. URL http://www.actu-
environnement.com/ae/news/loi-economie-bleue-rejet-boues-dragage-interdit-2020-
26496.php4 (accessed 4.1.16).
Bouguettaya, A.W., 2012. La programmation automatique: outils de la programmation
génétique et certaines de ses applications (Thèse de doctorat). Université de Tlemcen.
Bourabah, M.A., 2012. Comportement mécanique des sols fins, Application à la valorisation
des sédiments de barrages en technique routière (Thèse de doctorat). Université de
Tlemcen, Algérie.
Bouterfas, M., 2012. Optimisation de la compacite du squelette granulaire (Thèse de
doctorat). Université de Tlemcen.
Brezocnik, M., Kovacic, M., 2003. Integrated Genetic Programming and Genetic Algorithm
Approach to Predict Surface Roughness. Materials and Manufacturing Processes 18,
475–491. doi:10.1081/AMP-120022023
Brundtland, G.H., 1987. Notre Avenir à Tous, rapport de la commission mondiale sur
l’Environnement et le Développement. Les Editions du Fleuve, Paris (traduction
française de Our Common Future.
Collette, Y., Siarry, P., 2002. Optimisation multiobjectif, Éditions Eyrolles. ed. Paris.
Dabo, D., 2008. Impact environnemental des mâchefers d’incinération d’ordures ménagères
(MIOM) valorisés en technique routière: caractérisation, expérimentations multi-
échelles et modélisation hydro-géochimique (Thèse de doctorat). Université de Paris
ENMP.
De Larrard, F., 2000. Structures granulaires et formulation des bétons, Etudes et recherches
des Ponts et Chaussées. Laboratoire Central des Ponts et Chaussées.
De Larrard, F., Ledee, V., Sedran, T., Brochu, F., Ducassu, J.-B., 2003. Nouvel essai de
mesure de compacité des fractions granulaires à la table à chocs. Bulletin des
laboratoires des ponts et chaussées 246–47, 101–115.
Deb, K., Agrawal, S., Pratap, A., Meyarivan, T., 2000. A fast elitist non-dominated sorting
genetic algorithm for multi-objective optimization: NSGA-II. Lecture notes in
computer science 1917, 849–858.
176
Références bibliographiques
Deb, K., Pratap, A., Agarwal, S., Meyarivan, T., 2002. A fast and elitist multiobjective
genetic algorithm: NSGA-II. Evolutionary Computation, IEEE Transactions on 6,
182–197.
Debieb, F., Kenai, S., 2007. Performances et durabilité du béton à base de granulats recyclés
contaminés par les chlorures et les sulfates (Thèse de doctorat). Université des
sciences et de la technologie Houari Boumedienne.
Devaux, P., Vecoven, J., 1984. Les dechets de sables de fonderie en technique routiere.
Bulletin of the International Association of Engineering Geology-Bulletin de
l’Association Internationale de Géologie de l’Ingénieur 30, 375–378.
Dia, M., Ramaroson, J., Nzihou, A., Zentar, R., Abriak, N.E., Depelsenaire, G., Germeau, A.,
2014. Effect of Chemical and Thermal Treatment on the Geotechnical Properties of
Dredged Sediment. Procedia Engineering, SYMPHOS 2013 - 2nd International
Symposium on Innovation and Technology in the Phosphate Industry 83, 159–169.
doi:10.1016/j.proeng.2014.09.034
Dubois, V., 2006. Etude du comportement physico-mécanique et caractérisation
environnementale des sédiments marins: valorisation en technique routière (Thèse de
doctorat). Université d’Artois.
El Barrak, M., 2005. Contribution à l’étude de l’aptitude à l’écoulement des bétons
autoplaçants à l’état frais (Thèse de doctorat). Université de Toulouse 3.
El Hilali, A., Ghorbel, E., Gonnon, P., 2006. Influence des fillers sur l’ouvrabilité des bétons
autoplaçants. Actes 24èmes Rencontres Universitaires de Génie 155, 1–2.
Farina, M., Deb, K., Amato, P., 2004. Dynamic multiobjective optimization problems: test
cases, approximations, and applications. Evolutionary Computation, IEEE
Transactions on 8, 425–442.
FFB, 2014. Déchets de chantier-Les réponses aux questions que vous vous posez, Fédération
Francaise du Batiment.
Folino, G., Pizzuti, C., Spezzano, G., 2000. Genetic Programming and Simulated Annealing:
A Hybrid Method to Evolve Decision Trees, in: Poli, R., Banzhaf, W., Langdon,
W.B., Miller, J., Nordin, P., Fogarty, T.C. (Eds.), Genetic Programming, Lecture
Notes in Computer Science. Springer Berlin Heidelberg, pp. 294–303.
Fonseca, C.M., Fleming, P.J., others, 1993. Genetic Algorithms for Multiobjective
Optimization: FormulationDiscussion and Generalization., in: ICGA. Citeseer, pp.
416–423.
François, T., 1996. Réutilisation des sables de fonderie dans le bâtiment et les travaux publics
(No. ADEME-94-02-0037).
Gandomi, A.H., Alavi, A.H., 2011a. Multi-stage genetic programming: A new strategy to
nonlinear system modeling. Information Sciences 181, 5227–5239.
doi:10.1016/j.ins.2011.07.026
Gandomi, A.H., Alavi, A.H., 2011b. A new multi-gene genetic programming approach to
nonlinear system modeling. Part I: materials and structural engineering problems.
Neural Comput & Applic 21, 171–187. doi:10.1007/s00521-011-0734-z
Gandomi, A.H., Sajedi, S., Kiani, B., Huang, Q., 2016. Genetic programming for
experimental big data mining: A case study on concrete creep formulation.
Automation in Construction. doi:10.1016/j.autcon.2016.06.010
Gérard Bernier, 2004. Formulation des bétons. Ed. Techniques Ingénieur.
Hanini, W., 2012. Etude de l’incidence de la granularité sur la compacité des mélanges
granulaires (Thèse de doctorat). Université de Tlemcen.
Horn, J., Nafpliotis, N., Goldberg, D.E., 1994. A niched Pareto genetic algorithm for
multiobjective optimization, in: , Proceedings of the First IEEE Conference on
Evolutionary Computation, 1994. IEEE World Congress on Computational
177
Références bibliographiques
178
Références bibliographiques
179
Références bibliographiques
Sadok, A., Mekerta, B., 2012. Valorisation des sédiments dans un Centre de Stockage des
déchets. Presented at the Séminaire de la Technologie Spécialité Génie Civil -STSGC
du 16 Mai 2012, Centre universitaire de Témouchent, Algérie.
Sadok, A., Zentar, R., Abriak, N.-E., 2016a. Determination of water demand of dredged
sediment in concrete formulations. Presented at the 5th International Symposium on
Sediment Management 10 au 13 juillet 2016, Montreal- Canada.
Sadok, A., Zentar, R., Abriak, N.-E., 2016b. Genetic programming for granular compactness
modelling. European Journal of Environmental and Civil Engineering 1–13.
Sadok, A., Zentar, R., Abriak, N.-E., 2016c. L’état Sature Surface Humide (SSH) des
matériaux fins ou a dominance fine dans les formulations de béton. Revue Paralia,
XIVèmes Journées Nationales Génie Côtier – Génie Civil Vol 14.
Sadok, A., Zentar, R., Abriak, N.E., 2016d. Mix design of high-performance concrete by
Applying the Genetic Algorithm and Genetic Programming. Submitted to Engineering
Applications of Artificial Intelligence.
Sadok, A., Zentar, R., Abriak, N.E., 2014. Design method of civil engineering materials with
evolutionary algorithm. Presented at the 4 th International Symposium on Sediment
management CD 17-19 september 2014, Ferrara, Italy, p. 128.
Schaffer, J.D., 1985. Some experiments in machine learning using vector evaluated genetic
algorithms (Thèse de doctorat). Université de Vanderbilt , Nashville, TN (USA).
Sebaibi, N., 2010. Valorisation des composites thermodurcissables issus du recyclage dans
une matrice cimentaire: application aux bétons a ultra-haute performance (Thèse de
doctorat). Université d’Artois.
Sedran, T., 1999. Rhéologie et rhéométrie des bétons. Application aux bétons autonivelants
(Thèse de doctorat). l’Ecole Nationale des Ponts et Chaussées.
Sedran, T., de Larrard, F., 1994. RENE-LCPC–Un logiciel pour optimiser la granularité des
matériaux de Génie Civil. Bulletin des Laboratoires des Ponts et Chaussées 194, 87–
93.
Semcha, A., 2006. Valorisation des sédiments de dragage : Applications dans le BTP, cas du
barrage de Fergoug (Thèse de doctorat). Université de Reims.
Simonin, F., 2000. Comportement thermomécanique de bétons réfractaires alumineux
contenant du spinelle de magnésium (Thèse de doctorat). INSA, Lyon.
Smolczyk, H., 1978. L’emploi du ciment de laitier de haut fourneau dans le béton arme et le
béton précontraint (Partie 3). REV METALL Vol 5.
Srinivas, N., Deb, K., 1994. Muiltiobjective optimization using nondominated sorting in
genetic algorithms. Evolutionary computation 2, 221–248.
Thanh, T.N., 2009. Valorisation de sédiments marins et fluviaux en technique routière (Thèse
de doctorat). Université de Sherbrooke.
Waligora, J., 2008. Obtention de granulats stables de laitier d’aciérie de conversion LD pour
une utilisation en technique routière : étude appliquée au site d’ArcelorMittal
Dunkerque (Thèse de doctorat). Université de Lille 1.
Wikipédia, 2016. Développement durable. Wikipédia.
Wirquin, E., Hadjieva-Zaharieva, R., Buyle-Bodin, F., 2000. Utilisation de l’absorption d’eau
des bétons comme critères de leur durabilité—Application aux bétons de granulats
recyclés. Materials and structures 33, 403–408.
Yammine, J., Chaouche, M., Guerinet, M., Moranville, M., Roussel, N., 2008. From ordinary
rhelogy concrete to self compacting concrete: A transition between frictional and
hydrodynamic interactions. Cement and Concrete Research 38, 890–896.
doi:10.1016/j.cemconres.2008.03.011
180
Références bibliographiques
Zentar, R., Sadok, A., Abriak, N.-E., 2016. Packing density of fine and sandy materials,
Experimental and Modeling approaches. Presented at the 5th International Symposium
on Sediment Management 10 au 13 juillet 2016, Montreal- Canada.
Zhao, Z., 2014. Valorisation des sables de béton recyclé pour la fabrication de mortiers
(Thèse de doctorat). Université de Lille 1.
Zitzler, E., Thiele, L., 1999. Multiobjective evolutionary algorithms: a comparative case study
and the strength Pareto approach. evolutionary computation, IEEE transactions on 3,
257–271.
Zri, A., 2010. Mise en place d’une nouvelle approche de formulation d’une matrice cimentaire
à base de sable de dragage: application aux bétons de sables et granulats (Thèse de
doctorat). Université de Lille 1.
Zri, A., Abriak, N.-E., Zentar, R., 2009. Etude comparative du comportement physico-
mécanique d’un sédiment pollué et non pollué. Déchets Sciences et Techniques 13–19.
181
Annexes
Annexes 1 : Base de données des Bétons à Hautes Performances utilisée.
182
Annexes 2 : Base de données développée des mortiers.
C (gr) E/C (%) F/S (%) S (gr) Rc28 (MPa) Rf28 (MPa) Maniabilité (S)
450 40 0 1350 40,43 7,96 55
450 40 0 1350 36,94 7,59 55
450 40 0 1350 42,18 8,35 55
450 50 0 1350 48,64 9,79 8
450 50 0 1350 48,61 8,4 8
450 50 0 1350 46,76 8,66 8
450 60 0 1350 43,07 9,07 3
450 60 0 1350 41,38 9,99 3
450 60 0 1350 43,36 9,2 3
450 70 0 1350 29,93 7,05 1
450 70 0 1350 28,26 8,04 1
450 70 0 1350 29,67 8,98 1
400 40 0 1350 11,84 2,31 32
400 40 0 1350 10,22 3,03 32
400 40 0 1350 11,4 2,74 32
400 50 0 1350 56,35 8,9 10
400 50 0 1350 55,03 9,2 10
400 50 0 1350 54,79 8,47 10
400 60 0 1350 42,34 7,29 2
400 60 0 1350 42,17 7,47 2
400 60 0 1350 42,34 8,35 2
400 70 0 1350 38,15 6,77 1
400 70 0 1350 39,43 7 1
400 70 0 1350 38,99 7 1
450 40 5 1350 13,62 2,13 52
450 40 5 1350 14,18 2,93 52
450 40 5 1350 14,1 2,61 52
450 50 5 1350 53,26 8,96 7
450 50 5 1350 51,33 5,87 7
450 50 5 1350 50,03 6,83 7
400 60 5 1350 49,6 9 3
400 60 5 1350 48,06 9 3
400 60 5 1350 46,62 9 3
400 70 5 1350 42,34 8 1
400 70 5 1350 43,66 8 1
400 70 5 1350 41,3 8 1
450 40 10 1350 8,08 2,38 57
450 40 10 1350 10,99 2 57
450 40 10 1350 11,51 2,38 57
450 50 10 1350 49,89 10,5 21
450 50 10 1350 49,3 10,15 21
450 50 10 1350 49,3 10,84 21
400 60 10 1350 50,3 7,83 11
400 60 10 1350 50,46 7,07 11
400 60 10 1350 50,23 6,44 11
400 70 10 1350 48,15 7,37 3
400 70 10 1350 47,4 6,68 3
400 70 10 1350 47,7 6,68 3
450 40 15 1350 7,3 0,84 85
450 40 15 1350 6,46 1,21 85
450 40 15 1350 6,4 1,21 85
450 50 15 1350 51,32 7,3 43
450 50 15 1350 50,74 8,03 43
450 50 15 1350 49,94 8,03 43
400 60 15 1350 42,81 8,75 11
400 60 15 1350 40,39 9,01 11
400 60 15 1350 41,93 8,76 11
400 60 10 1350 45,61 9,17 11
400 60 10 1350 45,4 10,06 11
400 60 10 1350 46,49 8,79 11
400 70 15 1350 33,37 8,83 3
400 70 15 1350 34,99 8,23 3
400 70 15 1350 37,46 8,77 3
400 40 20 1350 5,61 2,32 97
400 40 20 1350 5,06 2,48 97
400 40 20 1350 4,12 2,02 97
400 50 20 1350 15,19 5,57 58
400 50 20 1350 18,8 5,29 58
400 50 20 1350 15,78 4,44 58
450 60 20 1350 39,81 8,7 14
450 60 20 1350 37,5 9,53 14
183
450 60 20 1350 39,81 9,42 14
450 70 20 1350 35,71 9,57 2
450 70 20 1350 34,73 8,58 2
450 70 20 1350 35,71 8,71 2
400 40 30 1350 3,03 1,52 122
400 40 30 1350 1,87 1,52 122
400 40 30 1350 2,4 1,52 122
400 50 30 1350 7,83 1,84 103
400 50 30 1350 5,76 1,84 103
400 50 30 1350 8,41 1,84 103
450 60 30 1350 34,42 9,04 71
450 60 30 1350 35,63 7,83 71
450 60 30 1350 35,4 9,16 71
450 70 30 1350 22,48 5,72 71
450 70 30 1350 19,15 5,86 13
450 70 30 1350 20,88 5,72 13
400 40 40 1350 2,11 1 300
400 40 40 1350 1,52 1 300
400 40 40 1350 2,11 1 300
400 50 40 1350 2,44 1 187
400 50 40 1350 2,44 1 187
400 50 40 1350 2,44 1 187
450 60 40 1350 30,5 7,84 126
450 60 40 1350 30,19 8,29 126
450 60 40 1350 30,44 7,34 126
450 57,4 32 1350 41,83 7,07 5
450 57,4 32 1350 42,21 8,24 5
450 57,4 32 1350 39,95 6,97 5
184
Publications dans le même domaine
Sadok, A., Zentar, R., and Abriak, N.-E. (2016). Genetic programming for granular
compactness modelling. European Journal of Environmental and Civil Engineering 1–13.
(Impact Factor: 0,636)
Sadok, A., Zentar, R., and Abriak, N.-E. (2016). L’état Sature Surface Humide (SSH) des
matériaux fins ou a dominance fine dans les formulations de béton. Revue Paralia 14.
Sadok, A., Zentar, R., and Abriak, N.-E. (2016). Design of high-performance concrete
mixture using the Genetic Algorithm and Genetic Programming. Submitted to Advances in
Engineering Software. (Impact Factor: 1,673)
Sadok, A., Zentar, R., Abriak, N.-E., 2016. Determination of water demand of dredged
sediment in concrete formulations. Presented at the 5th International Symposium on Sediment
Management 10 au 13 juillet 2016, Montreal- Canada.
Zentar, R., Sadok, A., Abriak, N.-E., 2016. Packing density of fine and sandy materials,
Experimental and Modeling approaches. Presented at the 5th International Symposium on
Sediment Management 10 au 13 juillet 2016, Montreal- Canada.
Mekerta, B., Semcha, A., Sadok, A., Guessas, H., Rahmani, M.C., 2016. Moyens de lutte
contre la prediction de l’envasement des retenues de barrages d’Ouizert, Bouhanifia et
Fergoug. Presented at the 5th International Symposium on Sediment Management du 10 au 13
juillet 2016, Montreal- Canada.
Sadok, A., 2016. Méthode de formulation d’un BHP par des Algorithmes Génétique et
Programmation Génétique. Presented at the Rencontres Universitaires de Génie Civil 26 au 29
Mai 2016, Liège-Belgique.
Sadok, A., Zentar, R., Abriak, N.E., 2015. The optimization of the concrete granular
compacity by machine learning approaches. Presented at the 9th International SedNet
conference, Solving societal challenges: working with sediments 23-26 September 2015,
Krakow, Poland.
Sadok, A., Zentar, R., Abriak, N.E., 2015. Modélisation de la compacité granulaire par
une approche d’Apprentissage Automatique. Presented at the 33ÈMES Rencontres
Universitaires de Génie Civil du 27 au 29 mai 2015, Bayonne-France.
185
Benhamou, Y., Ameur, M., Mekerta, B., Kazi Aoual, F., Sadok, A., 2015. Pathologie des
sols dans la région de Saida. Presented at the 3ème Colloque International, Les sols non
saturés :du grain à l’ouvrage UNSAT Batna du 16 et 17 novembre 2015, Batna, Algérie.
Sadok, A., Zentar, R., Abriak, N.E., 2014. Design method of civil engineering materials
with evolutionary algorithm. Presented at the 4 th International Symposium on Sediment
management CD 17-19 september 2014, Ferrara, Italy, p. 128.
Rahmani, M.C., Mekerta, B., Allal, A., Sadok, A., 2015. Formulation d’une barrière
ouvragée d’un centre de stockage des déchets à partir de matériaux naturels. Presented at the
2ieme édition de la Conférence Internationale Francophone NoMad du 5 et 6 Novembre 2015,
Douai, France.
Sadok, A., Zentar, R., Abriak, N.-E., 2016. Formulation des matériaux alternatifs dans les
bétons par des approches d’Intelligence Artificielle. Presented at the Séminaire International
sur : Les Matériaux locaux et la construction Durable, Adrar, Algérie.
Sadok, A., Mekerta, B., Allal, A., Rahmani, M.C., 2014. Formulation of a compacted clay
liner for solid waste landfill from a material of sedimentary origin. Presented at the 4 th
International Symposium on Sediment management CD du 17-19 september 2014, Ferrara,
Italy, p. 128.
Sadok, A., Mekerta, B., 2013. Matériaux écologique pour le confinement des déchets.
Presented at the International Congress on Materials & Structural Stability-CMSS du 27 au 30
Novembre 2013, Rabat, Morroco.
Sadok, A., Mekerta, B., 2012. Valorisation des sédiments de dragage en matériaux de
construction. Presented at the 1 er Congrès International du Génie Civil et d’Hydraulique–
LGCH du 10et 11 Décembre 2012, Guelma, Algérie.
Sadok, A., Mekerta, B., 2012. Valorisation des sédiments dans un Centre de Stockage des
déchets. Presented at the Séminaire de la Technologie Spécialité Génie Civil -STSGC du 16
Mai 2012, Centre universitaire de Témouchent, Algérie.
186
RESUME
La question environnementale est désormais une préoccupation majeure pour les gouvernements
et les industriels. La gestion efficace des déchets s’inscrit dans les actions prioritaires dans une optique
d’économie circulaire. Cette gestion efficace préconise en premier lieu la valorisation maximale des
déchets. Le domaine du génie civil est particulièrement concerné.
La valorisation des matériaux alternatifs dans le domaine du génie civil connaît un essor important
ces dernières années, toutefois cette pratique connait certaines limites, notamment le manque de
réglementation claire et surtout le manque d’outils et méthodes de formulation adaptées et adaptables à
l’utilisation des matériaux alternatifs.
Dans ce contexte, le présent travail porte sur la mise en place de méthode de formulation de
matériaux de construction à base de sous-produits industriels. Cette méthode hybride basée sur des
outils de l’intelligence artificielle, combine les Algorithmes Génétique (AG) comme outils
optimisation multiobjectifs et la Programmation Génétique (PG) dans ses deux versions, la PG
classique et la MGGP (MultiGène Genetic Programming) comme outils de modélisation des
problèmes complexes par Apprentissage Automatique.
La méthode de formulation des matériaux de construction à base de sous produits industriels,
proposée dans ce travail a été testée sur deux cas d’étude (la formulation de bétons à hautes
performances et la formulation de mortiers à base de matériaux alternatifs ou multi-déchets) et validée
par des testes au laboratoire. Les résultats obtenus sont concluants et très prometteurs.
Mots Clés : Valorisation, matériaux alternatifs, formulation, béton, Algorithmes Génétiques,
Programmation Génétique, Intelligence artificielle, Apprentissage automatique
ABSTRACT
The environmental issue has become a major concern for governments and industry. Effective
waste management is part of the priority actions in order to achieve a green circular economy. This
efficient management imposes first maximum recovery of waste, knowing the large tonnage produced
in different sectors. The field of civil engineering is particularly concerned.
The valorization of alternative materials in the field of civil engineering has grown significantly in
recent years. However, this practice knows some limitations, including the lack of clear regulations
and especially the lack of tools and methods suitable for design of materials including alternatives
materials.
In this context, the present work focuses on the development of mix design method of building
materials based on industrial by-products. This hybrid method combines the Genetic Algorithms (GA)
as multi-objective optimization tools and Genetic Programming (GP) with his two versions, classical
GP and MGGP (MultiGene Genetic Programming) as modeling tools for complex problems by
Machine Learning approach.
Specific studies were carried out also or these innovative tools, to demontrates benefits and
weaknesses of these tools on these applications in civil engineering.
Finally, this method of formulation of building materials based on industrial sub products
proposed in this work was tested on two case studies (design of high performance concrete and
mortars made of alternative materials) and validated by the laboratory tests. The results are conclusive
and promising to generalize the method to other applications of Civil Engineering.
Keywords: Valorization, alternative materials, mixes design, concrete, Genetic Algorithm,
Genetic Programming, artificial intelligence, machine learning.