Energie Renouvelable Et Applications Francophone Issue

ISSN: 1737-9334
ICID
International Centre for Innovation & Development
-PET- Vol. 52
Proceedings of Engineering & Technology

-PET-
Energie Renouvelable et Applications
Francophone issue
Editeur:
Dr. Ahmed Rhif
International Centre for Innovation & Development

–ICID–
Editeur en chef Fathi Mohammed (ALG)
Dr. Ahmed Rhif (Tunisie) Sallam Mebrouk (ALG)

Ahmed.rhif@gmail.com
Youcef Soufi (ALG)
Dean of International Centre for Hassane Mahmoudi (MAR)
Innovation & Development (ICID)
Toufik Madani Layadi (ALG)
Comité de redaction
Samira Chader Kerdjou (ALG)
Elmati Khoumri (MAR) Tounzi Abdelmounaïm (FRA)
Driss Youssfi (MAR) Zohra Ameur (ALG)
Fabien Delaleux (FRA)
Georges Descombes (FRA)
Hana Boudjedir (ALG)
Fawaz Massouh (FRA)
Nachida K. Merzouk (ALG)
Houria Siguerdidjane
Khenfer Nabil (ALG)
Mohamed Elamine Djeghlal (ALG)
Manar Slimane (MAR)
Abdelfettah Barhdadi (ALG)
Mohamed Gherbi (ALG)
Mohamed Benbouzid (FRA)
Sara Zatir (ALG)
Philippe Guibert (FRA)
Mustapha Hatti (ALG)
Saoussen Hammami (TUN)
Rahmani Lazhar (ALG)
Tahar Bahi (ALG)

Sommaire
Analyse comparative des perturbations EM conduites BF engendrées par les ponts redresseurs Page 1
triphasés commandés et non commandés.
Saliha Rezini, Zin-eddine Azzouz.
Commande d’un système éolien basé sur une MADA. Page 5

Moulay Fatima, Habbatti Assia, Hamdaoui Habib.
Le biogaz source d’energie renouvelables. Cas des déchets ménagers « biodéchet ». Page 11
Baichata Malika, Tamali Mohamed.
Emissions BF des chargeurs commandés par régulateur à Hystérésis embarqués dans les Page 16
véhicules électriques hybrides rechargeables (V.E.H.R).
L. Bessoltane, Z. Azzouz, D. Toumi.
Stratégie de contrôle du raccord au réseau du générateur photovoltaïque. Page 20

Habbati Bellia Assia, Moulay Benine Fatima.
Analyse de l’écoulement du flux de puissance dans un système éolien en présence d’un creux Page 24
de tension.
Salima Kadri, Mahmoud Habbab, Cherif Benoudjafer.
Écoulement hydromagnetique non lineaire de convection forcee d’un fluide biomagnetique instable. Page 29
Merabet Ayeche Chahra, Kezzar Mohamed, Sari Mohamed Rafik.
Le bâtiment passif pour assurer l'efficience énergétique. Page 35

Mokhnache nardjesse, Alkama djamel.
Le CuInSe2 en couches minces est une structure de base pour une application photovoltaïque. Page 44
Ghalmi Leila, Bensmaine Souhila.
Les intentions bioclimatiques des architectes de la Reconstruction tunisienne d’après-guerre: cas Page 44
des logements collectifs d’El Menzah 1.
Salma Gharbi.
La reproduction de l'habitation troglodyte vernaculaire aux temps 'modernes'. Cas d'un Page 53
projet d'habitat du 'plan de Constantine' a Biskra (sud-est Algérien).
Asma Achraf Zendagui, Azeddine Belakehal.
Comparaison entre les modules d’élasticité d’un béton autoplaçant exposé à hautes température par Page 58
deux méthodes de fréquence de résonance.
Benoudjafer Imane, Benoudjafer Ibtissam.
Amélioration de la performance énergétique d’un bâtiment résidentiel (cas de la ville de Bechar). Page 63
Ibtissame Benoudjafer, Imane Benoudjafer.
Etude comparative entre les différentes méthodes extrapolations verticales de la vitesse du vent Page 70
de la loi à coefficient de puissance avec le langage PYTHON.
Hadjij Merouane, Benfekir Abderrahim, Zoubida Lounis.
Caractérisation physico-chimique et mécanique des argiles de la région d’Adrar et de Timimoun Page 76

en vue d’une utilisation dans les centres de stockages des déchets.
Belaidi Khedidja, Mekerta Belkacem, Allal Mohamed Amine.
Pour une approche environnementale du traitement des eaux usées de la ville de Tamanrasset –Sud Page 80
Algérien- dans le cadre du développement durable.
Bouaoud Mohamed, Bebouche Manel, Hadj Ali Zakaria.
3ème Congrès International sur les Energies Renouvelables et le Développement Durable (ERDD-2019)
Proceedings of Engineering and Technology PET
Analyse comparative des perturbations EM conduites

BF engendrées par les ponts redresseurs triphasés
commandés et non commandés
Saliha Rezini #, Zin-eddine Azzouz *
#
Département de Génie Electrique, Faculté des Sciences et de la Technologie, Université Abdelhamid Ibn Badis
Site 1, Route Belahcel 27000, Mostaganem, Algerie
rezini_s@yahoo.fr
*
Département d’Automatique, Faculté de Génie Electrique Université des Sciences et de la technologie-Mohamed Boudiaf
Laboratoire de Développement des Entraînements Electriques (LDEE)
El Mnaouar, B.P. 1505, Bir El Djir, 31000, Oran, Algérie
zinazzouz@yahoo.fr
Résumé— La pollution électromagnétique (EM) est l'un des doivent être conformes aux normes CEM. Afin d’obtenir cette
problèmes majeurs intervenant dans la dégradation de la qualité conformité, une méthode classique, consiste à réaliser un
de l'énergie électrique dans les réseaux électriques. Cette prototype, puis à mesurer ses émissions. On lui ajoutera enfin
pollution est en partie créée par les dispositifs d’électronique de un filtre au niveau de son alimentation afin de limiter son
puissance notamment les convertisseurs statiques assurant
niveau d’émission en dessous du gabarit imposé par la norme
certaines fonctions dans le réseau électrique. Ces dispositifs sont
considérés comme des charges non linéaires émettant des [8],[9].
perturbations électromagnétique conduites et rayonnées très
contraignantes du point de vue compatibilité électromagnétique Ce processus est long et coûteux, de plus, les filtres CEM
(CEM). sont lourds, volumineux et chers [10]. Cette méthode est loin
Les travaux de recherche présentés dans cet article ont pour d’être optimale car elle ne prend en compte la contrainte CEM
objectif l’étude et la caractérisation de la signature CEM des qu’à la fin du développement du dispositif. Pour s’affranchir
perturbations électromagnétique conduits basses fréquences de ces inconvénients, il est intéressant d’intégrer la contrainte
engendrées par les ponts redresseurs triphasés non commandés et CEM dès la phase de conception en usant d’un outil de
commandés. A noter enfin que tous les résultats de simulation,
modélisation des perturbations EM. Ce dernier servirait à
obtenus à l’aide du logiciel Pspice, ont été validés par des gabarits
de normes spécifiques aux dispositifs étudiés. tester différentes configurations de dispositifs d’électronique
de puissance et de choisir la moins polluante du point de vue
Mots-clés — Compatibilité électromagnétique, electronique de CEM, sans pour autant arriver au stade de réalisation du
puissance, pont redresseur, perturbation EM conduite, basses dispositif [10].
fréquences.
Ainsi, ce travail est consacré à l’étude des perturbations EM
I. INTRODUCTION conduites basses fréquences engendrées par les ponts
Le domaine du génie électrique ne cesse de gagner du redresseurs et transmises au réseau électrique. Pour mettre en
terrain dans les domaines industriels. Cette tendance n’a fait évidence ces perturbations, des simulations dans le domaine
que se renforcer depuis quelques années grâce aux progrès de temporel et fréquentiel ont été effectuées.
l’électronique de puissance qui ont été considérables dans la
maîtrise de la conversion et l’amélioration de la capacité de II. NORMALISATION CEM EN ELECTRONIQUE DE PUISSANCE
transit et de l’énergie électrique distribuée. Cette énergie qui Depuis plusieurs années, les fabricants d'appareils
est produite et distribuée sous forme de tensions sinusoïdales, électriques et ou électroniques sont soumis à des normes de
permet de fournir la puissance électrique nécessaire aux CEM dépendant de leurs produits et leurs marchés d'affaires.
différentes charges connectées au réseau électrique [1]–[4]. Ces normes ont pour but de diminuer le risque de causer des
Lorsque la forme d'onde de la tension n'est plus sinusoïdale, problèmes d'interférences EM. Il existe à ce jour un certain
des perturbations EM peuvent affecter le bon fonctionnement nombre de normes auxquelles doit se soumettre tout produit
de nombreux éléments raccordés au réseau et mettre ainsi en devant être connecté au réseau de distribution [11], [12].
jeu des problèmes de CEM entre ces derniers [5]–[7].
A titre d’exemple, les normes CEI classe A 1000-3-2 et CEI
L'utilisation croissante des dispositifs d'électronique de classe A 1000-2-2 présentent respectivement les limites des
puissance dans les réseaux électriques contribue à la courants injectés dans le réseau public basse tension et les
détérioration de la tension et du courant délivrés par ces niveaux de compatibilité pour les tensions harmoniques sur ce
réseaux. Aussi, les dispositifs d’électronique de puissance même réseau [13]–[15].
Copyright - 2019 Page 1

ISSN 1737-9334
Pour le contrôle du taux de distorsion harmonique (THD)

présent sur les réseaux électriques, on doit utiliser la norme Interrupteurs de puissance
(Diodes, thyristors) Source de perturbations
IEEE Std 519-2014 [16], [17]. Elle impose également les Réseau Charge
limites recommandées de la pollution harmonique engendrée électrique
par les usagers et de la distorsion harmonique totale sur le
réseau électrique.
VRes RRes LRes
III. IDENTIFICATION ET CARACTERISATION DES PERTURBATION 
EM GENEREES PAR LES PONRS REDRESSEURS TRIPHASES
La notion de perturbations basses fréquences comprend les
 Rch Cch
influences EM dans la gamme des fréquences comprises dans 

l’intervalle de 0Hz à 9kHz [18] – [20]. Ces perturbations sont
causées par des utilisateurs possédant des appareils ayant des
caractéristiques non-linéaires connectés au réseau
d’alimentation.
Les ponts redresseurs monophasés présentent un intérêt très

limité parce que les formes d’onde de la tension et du courant Fig. 1 Structure du pont de diode triphasé débitant sur une charge-Modèle
délivrés par ces derniers sont loin de la forme souhaitée basses frequencies
(forme continue). En revanche, les ponts redresseurs triphasés
Les éléments parasites ne sont pas pris en considération et
permettent de redresser les courants et les tensions à une forme
les paramètres de simulation considérés sont consignés sur le
très proche de la forme continue. Les signaux redressés ont
alors une ondulation de fréquence six fois supérieure à celle tableau 1.
des tensions d'entrée et un taux d’ondulation très petit. TABLEAU I
PARAMETRES DE SIMULATION UTILISES
Dans ce paragraphe, nous allons examiner l’influence des

perturbations EM créées par des ponts redresseurs triphasés CCh
VRes RRes LRes RCh
non commandés et commandés. Cette étude passe par la mise (V) (Ω) (µH) (Ω) Pont non Pont
en évidence et la caractérisation des perturbations basses commandé commandé
fréquences générées par ces ponts redresseurs. Ainsi, nous
220 0.4 500 50 900 µF 900 nF
étudions dans un premier temps, les perturbations produites
par un pont redresseur triphasé à diodes. Nous étudions ensuite,
les perturbations générées par un pont redresseur triphasé Résultats et analyse
commandé (pont tout thyristors et pont mixte). Les diodes de Les figures 2 et 3 représentent respectivement la tension à
redressement utilisées sont de type BYT12P-600 et les l’entrée du pont redresseur et le courant de ligne. Les analyses
thyristors sont de type MCR3935-10. spectrales de Fourier de ces deux grandeurs sont présentées
respectivement dans les figures 4 et 5.
Les thyristors fonctionnent avec un retard à l’amorçage de
30° par rapport à la conduction naturelle des diodes (amorçage L’analyse de ces résultats montre clairement l’effet
naturel) dans les configurations de ponts considérées. A noter perturbateur engendré par le pont redresseur triphasé (dans les
que le signal de commande des thyristors est de forme trois configurations étudiées à savoir : le pont redresseur non
trapézoïdale. commandé, le pont redresseur commandé tout thyristors et
enfin le pont redresseur commandé mixte), et sa charge. Cet
Par ailleurs, afin de mettre en évidence le principe de effet est traduit par une déformation des formes d’ondes de
création des perturbations basse fréquence, nous avons tension et du courant de ligne. De plus, on remarque que les
considéré le modèle de base illustré dans la figure ci dessous. formes d’ondes de tension et de courant de ligne obtenues,
Aussi, pour chaque pont redresseur nous avons tracé les correspondant au pont redresseur commandé mixte, présentent
variations temporelles et fréquentielles (analyse spectrale de une dissymétrie entre l’onde positive et l’onde négative.
Fourier) de la tension à l’entrée du pont et du courant de ligne.
Cette déformation se traduit, sur le plan fréquentiel, par
l’apparition de plusieurs harmoniques. Nous constatons aussi
que les harmoniques multiples de trois sont nuls pour les trois
configurations de pont étudiées.
Copyright - 2019
ISSN 1737-9334 Page 2
Pont redresseur triphasé non commandé

Pont redresseur triphasé commandé mixte
Dans les tableau II et III, nous avons reporté, pour chaque
Pont redresseur triphasé commandé tout thyristor configuration, les valeurs du THD, du TDD et du fondamental
800
800V
obtenues. Les valeurs du THD et du TDD sont comparées et
validées par les données de la norme IEEE 519-2014.
400
400V
TABLEAU II
TABLEAU COMPARATIF DES VALEURS DU THD,TDD ET DE L’AMPLITUDE DU
Tension (V)
FONDAMENTAL RELATIF AU PONT REDRESSEUR NON COMMANDE

0
0V
Fondamental (V) 530.5

Tension à Simulation 3.45
-400
-400V
l’entrée du THD
pont (%) Norme (IEEE std 8
-800
519-2014)
-800V
20
20ms 22
22ms 24
24ms 26
26ms 28
28ms 30
30ms 32
32ms 34
34ms 36
36ms 38
38ms 40ms40
V(R5:2,R6:2) V(V6:-,V7:-) V(D7:K,X116:A) Fondamental (A) 11.61
Temps
Time(ms)
Courant de Simulation 86.92

Fig. 2 Variations temporelles de la tension composée à l’entrée du pont ligne TDD
redresseur triphasé (%) Norme (IEEE std 12
519-2014)
80
80A
TABLEAU III
40A
40
TABLEAU COMPARATIF DES VALEURS DU THD,TDD ET DE L’AMPLITUDE DU
FONDAMENTAL RELATIF AU PONT REDRESSEUR COMMANDE
Courant (A)
0
0A Tout
thyristors Mixte
-40
-40A Fondamental (V) 497.8 511.3
Tension à THD Simulation 5 6.42
l’entrée du (%)
-80
-80A
20
20ms 22
22ms 24
24ms 26
26ms 28
28ms 30
30ms 32
32ms 34
34ms 36
36ms 38
38ms 40ms40 pont Norme (IEEE std 8 8
-I(V10) I(R2) I(R9)
Temps
Time (ms)
519-2014)
Fondamental (A) 55.98 37.75
Fig. 3 Variations temporelles du courant de ligne du pont redresseur Courant de TDD Simulation 29.92 55.36
ligne (%)
Norme (IEEE 8 8
600
600V std 519-2014)
A la lumière de ces, il est facile de voir l’effet perturbateur

Spectre de la tension (V)
400
400V
engendré par le pont redresseur triphasé, sur le courant de
ligne et la tension à l’entrée du pont. Ainsi, l’analyse du THD
montre que la valeur de ce paramètre dans les configurations
200
200V
de ponts est inferieur à la valeur limite du THD spécifiée par la
norme IEEE 519-2014.
0
0V
0
0Hz 0.2
0.2KHz
V(R5:2,R6:2)
0.4
0.4KHz
V(V6:-,V7:-)
0.6
0.6KHz 0.8
0.8KHz
V(D7:K,X116:A)
1
1.0KHz 1.2
1.2KHz 1.4
1.4KHz 1.6
1.6KHz 1.8
1.8KHz 2.0KHz 2 De plus, nous constatons que l’amplitude du fondamental,
Fréquence (kHz)
Frequency
correspondant à la tension à l’entrée du pont redresseur, est en
dessous de l’amplitude de la tension fournie par le réseau
Fig. 4 Spectre fréquentiel de la tension composée à l’entrée du pont
redresseur triphasé
électrique. Par ailleurs, l’analyse du TDD dans les trois cas de
configurations de pont étudiées, montre que ce dernier dépasse
60
60A largement les limites recommandées par la norme IEEE 519-
2014.
Enfin nous remarquons, que le pont redresseur triphasé non

Spectre du courant (A)
40
40A
commandé est moins polluant que le pont redresseur triphasé

commandé (tout thyristors, mixte). En effet, les valeurs
20
20A
obtenues du THD et du TDD, correspondant au pont
redresseur non commandé, sont inferieures à celles du pont
redresseur commandé (tout thyristors, mixte) ce qui justifie les
0
0A
valeurs trouvées du fondamental.
0
0Hz 0.2
0.2KHz 0.4
0.4KHz 0.6
0.6KHz 0.8
0.8KHz 1
1.0KHz 1.2
1.2KHz 1.4
1.4KHz 1.6
1.6KHz 1.8
1.8KHz 2.0KHz 2
-I(V10) I(R2) I(R9) Fréquence (kHz)
Frequency
Fig. 5 Spectre fréquentiel du courant de ligne
Copyright - 2019
ISSN 1737-9334 Page 3
IV. CONCLUSIONS [14] G. Eduful, A. Jackson, J. E Cole , “Harmonic emission limits and
selecting PCC location based on the type of distribution system, ”
Le travail présenté dans cet article est une contribution à Proceedings of the world congress on engineering, Vol I, pp 301-306,
l’identification et la caractérisation des perturbations EM July 2-4, London, 2014
engendrées par les dispositifs d’électronique de puissance. [15] Copyright International Electrotechnical, Commission International,
Standard IEC 61000-2-2, 2002.
L’apport de ce travail concerne l’étude des agressions [16] IEEE Std 519-2014, IEEE recommended practices and requirements
électromagnétiques conduites basses fréquences générées par for harmonic control in electrical power systems.
les ponts redresseurs triphasés commandés et non commandés [17] P. Thakur, “Nonlinear load distribution and drives configuration
transmises au réseau électrique. selection for harmonic mitigation as per IEEE 519-2014
standard, ” International conference on power electronics, drives and
energy systems (PEDES), 18-19 December, India, 2018.
Les résultats obtenus après simulation, nous ont permis de [18] J. L. Cocquerelle , CEM et électronique de puissance, Editions Technip,
mettre en évidence les sources et les chemins de propagation Paris, 1999.
des perturbations engendrées par ces convertisseurs. Ces [19] R.T.Kennedy and I. Daut, “Power quality and electromagnetic
compatibility-The simple or not so simple ubiquitous power supply
résultats ont étés comparés et validés par des gabarits de input stage, ” 2nd International conference on power and energy
normes CEM relatifs aux dispositifs étudiés. Nous avons (PECon 08), 1-3 December, Malaysia, 2008.
montré dans cette étude que les ponts redresseurs constituent, [20] C. Beiu, N. Golovanov, R. Porumb and C. Toader, “Power quality in
en basses fréquences, une source non négligeable de relation to the requirements of electromagnetic compatibility, ” 5th
International symposium on electrical and electronics engineering
perturbations dont l’origine est due au fonctionnement des (ISEEE), 20-22 October, Romania, 2017.
diodes et des thyristors.
REFERENCES
[1] E. Félice and P. Révilla, Qualité des réseaux électriques et efficacité
énergétique, Dunod, L'Usine Nouvelle, 2009.
[2] A. M. Abu-Jalala, T. Cox , C. Gerada ,M. Rashed, T. Hamiti and N.
Brown, “ Power Quality Improvement of Synchronous Generators
Using an Active Power Filter, ”: IEEE Transactions on Industry
Applications, pp 4080 - 4090, Vol.54, Issue.5, 19 April, 2018.
[3] C. Sankaran, Power quality, International standard book, Library of
Congress card, Printed in the United States of America, 2002.
[4] L. Angela, N. P. Gabriel and M. D. Corina, “The influence of home
nonlinear electric equipment operating modes on power quality,”
Wseas transactions on systems, pp 357 – 367, E-ISSN: 2224-2678, Vol.
13, 2014.
[5] R. Calvas, Les perturbations électriques en BT, Cahier technique n° 141,
Groupe Schneider, 2001.
[6] M. Deobarro, “ Etude de l'immunité des circuits intégrés face aux
agressions électromagnétiques : proposition d'une méthode de
prédiction des couplages des perturbations en mode conduit, ” Thèse
de doctorat, Université de Toulouse, 2011.
S. Rezini, Z. Azzouz et B. Kouadri, “ Compatibilité électromagnétique
des convertisseurs statiques : action par conduction, ” 5ème Conférence
internationale sur électrotechnique, USTO-MB, Oran, 10-11 Novembre,
2009.
[7] P. M. Nicolae, I.D. Nicolae and M.Ş. Nicolae, “Electromagnetic
compatibility and power quality problems at low frequency for loads
from an urban transportation system and in a power substation, ” IEEE
International EMC Symposium, pp 242 – 247, Colorado Denver, 5-9
Aug, 2013.
[8] F. Costa, C. Gautier, E. Labouré and B. Revol, La compatibilité
électromagnétique en électronique de puissance. Principes et cas
d’études, Hermes-science, Lavoisier, Paris, 2013.
[9] S. Rezini, “Contribution à l’étude et à la caractérisation de
perturbations EM engendrées par des dispositifs d’électronique de
puissance affectant un réseau électrique, ”Thèse de doctorat,
Université des sciences et de la technologie-MB, 2016.
[10] C. Marlier, “Modélisation des perturbations électromagnétiques dans
les convertisseurs statiques pour des applications aéronautiques, ”
Thèse de doctorat, Université de Lille I, 2013.
[11] M. M. Abdusalam, “Structures et stratégies de commande des filtres
actifs parallèle et hybride avec validations expérimentales, ” Thèse de
doctorat, Université Henri Poincaré, Nancy-I, 2008.
[12] A. M. Abu-Jalala, T. Cox , C. Gerada ,M. Rashed, T. Hamiti and N.
Brown, “ Power Quality Improvement of Synchronous Generators
Using an Active Power Filter, ” IEEE Transactions on Industry
Applications, pp 4080 – 4090, Vol.54, Issue.5, 19 April, 2018.
[13] International standard, IEC 61000-3-2, Third edition, 2005.
Copyright - 2019
Page 4
ISSN 1737-9334
Commande d’un système éolien basé sur une MADA

Moulay Fatima#1, Habbatti Assia #1, Hamdaoui Habib #1
#1 Faculty of techology, IRECOM Laboratory
University of Djillali Liabes,Sidi Bel Abbes, Algeria.
1
fatimamoulay66@yahoo.fr
2
bellia_abdeljalil@yahoo.fr
3
hamdaoui_h@outlouk.fr
Résumé - Le présent travail est consacré à une commande Cet article présente une méthode de contrôle d’une
vectorielle de l’énergie active et réactive d’une machine Machine Asynchrone à Double Alimentation, il est organisé
asynchrone à double alimentation(MADA) intégrée à un système comme suit : une introduction puis un modèle mathématique
éolien. Tout d'abord, le modèle mathématique de l'ensemble du de l’ensemble, ensuite, la commande vectorielle et en dernier
système étudié est développé dans les coordonnées de Parc d-q puis
des résultats de simulations sous l’environnement
le transfert de puissance entre le stator et le réseau est réalisé par
action sur les signaux rotoriques via un convertisseur Matlab/Simulink et une conclusion.
bidirectionnel. Le contrôle indépendant de la puissance active et
II. MODELE MATHEMATIQUE DU SYSTEME EOLIEN
réactive est assuré par des régulateurs classiques (PI)). L’étude est
faite par simulation à base du logiciel Matlab/Simulink. Le système étudié est représenté par la figure suivante :
Mots clés: Eolienne - Machine Asynchrone à Double

Alimentation – Commande Vectorielle.
I. INTRODUCTION
Dans le secteur de l’énergie électrique, le monde entier
cherche à réduire considérablement sa dépendance aux
énergies de source fossile, se caractérisant à la fois par une
forte émission de gaz à effet de serre et par des prix instables.
Les opérateurs du secteur de l’énergie électrique s’acharnent
à diversifier leur bouquet énergétique plus spécialement par
une énergie propre et renouvelable de source géothermique, Fig. 1 Système de conversion d’énergie éolienne à base de la MADA
biomasse, solaire ou éolienne [1]. A. Modéle du vent:
L’un des axes actuels de recherche est la génération
La modélisation de l’éolienne exige la modélisation du
d’énergie électrique à l’aide de machine asynchrone à double
vent, du comportement aérodynamique des pales, du
alimentation, utilisant comme moyen d’entraînement
générateur électrique, du convertisseur de puissance et du
l’énergie éolienne. Intégrée dans un système éolien, la
système de commande. Le vent est la variable d’entrée du
MADA permet de fonctionner sur une large plage de vitesse
système éolien. La vitesse du vent est généralement
de vent, et d’en tirer le maximum de puissance possible pour
représentée par une fonction scalaire qui évolue dans le
chaque vitesse de vent. Son circuit statorique est connecté
temps, V = f(t). La vitesse moyenne du vent est donnée par
directement au réseau électrique.
l’expression suivant [3]:
Un second circuit placé au rotor est également relié au
réseau, mais par l’intermédiaire de convertisseurs de V_moy=8+0.2sin(0.1047t)+2sin(0.2665t)+sin(1.2930t)+0.
puissance [2]. Etant donné que la puissance rotorique 2sin(3.6645t) (1)
transitée est moindre, le coût des convertisseurs s’en trouve
réduit en comparaison avec une éolienne à vitesse variable B. Aérodynamique :
alimentée au stator par des convertisseurs de puissance. Durant la recherche bibliographique, plusieurs modèles du
C’est la raison principale pour laquelle on trouve cette coefficient de puissance ont été trouvés [4].
machine pour la production en forte puissance. Une seconde
raison est la possibilité de régler la tension au point de A titre d’exemple, on propose le modèle suivant :
connexion de cette machine.
Copyright - 2019
Page 5
ISSN 1737-9334
C_p (λ,β)=(0.3-0.00167*β) sin⁡[π(λ+0.1)/(10-0.3β)]- Avec :

0.00184.(λ-3).β (2)
Ω2 : vitesse de rotation après multiplicateur.
Il est clair que le coefficient de puissance dans la relation
(2) est en fonction de λ et de β C. Modélisation et application de la MADA dans le
système éolien :
La machine asynchrone à double alimentation se compose
principalement de deux parties. Le stator triphasé identique à
celui des machines asynchrones classique, et le rotor tourne à
l’intérieur de la cavité de la machine et est séparé du stator
par un entrefer. En principe les circuits électriques du stator
sont constitués de trois enroulements identiques couplés en
étoile (ou en triangle) à la seule différence est que celui du
rotor est relié aux bagues sur lesquelles glissent des balais.
Cette machine peut fonctionner comme générateur ou moteur.
Fig. 2 Evolution du coefficient de puissance de l'éolienne pour différente
Le stator de la MADA connecté directement au réseau et le
angle de calage β (°). rotor et connecté à un onduleur [5].
Sur la figure (2), plusieurs courbes sont distinguées mais Pour l’application dans un système éolien, le mode de
nous nous intéressons à celle qui possède le plus haut fonctionnement en génératrice est intéressant. En effet si la
sommet. Cette courbe est caractérisée par le point optimal plage de variation de vitesse ne dépasse pas (±) 30% en de ça
(λ_opt= 4.9,C_(p_max)= 0.3,β=0) qui est le point ou au-delà de la vitesse synchronisme, la machine est capable
correspondant au maximum du coefficient de puissance (C p- de débiter une puissance allant de 0.7 à 0.3 fois la puissance
max) et donc au maximum de la puissance mécanique nominale.
récupérée.
Si les grandeurs électriques triphasées sont équilibrées,
On sait que la puissance produite par un vent traversant alors les grandeurs homopolaires sont égales à zéro. Dès lors,
une surface S dépend du cube de la vitesse du vent V et de la les transformées de Park induisent des vecteurs comportant
densité de l'air p. Cette puissance est donnée par : seulement les deux composantes directe d et quadrature q
1 Pour les équations électriques :

P    S V 3 (3)
2 Vsd=Rs .I sd+(d  sd )/dt- sq  ω0
Où : S  R 2
Vsd=Rs .I sq+(d sq )/dt- sq  ω0
R : Rayon de la turbine (longueur de pale).
La puissance mécanique disponible sur le rotor P_m est

Vsd=Rr .I rd  (drd )/dt-rq ω0  ωr   (6)
généralement exprimée en utilisant le coefficient de puissance
et est donnée par : Vrq=Rr .I rq  (drq )/dt-rd ω0  ωr  
Pm 
1
  S  C _ P ,  V 3 (4) 0  d 0 dt
2
Ω1 : vitesse de rotation avant multiplicateur. 0 La vitesse angulaire arbitraire.
R : Rayon de l’aérogénérateur 35.25𝑚 Pour les équations de flux :
ρ : Densité de l’air 1.225 𝐾𝑔. 𝑚−1

 ds  Ls I ds+MI dr
Compte tenu du rapport du multiplicateur de vitesse, la  ds  Ls I ds  MI dr (7)

puissance mécanique 𝑃𝑚 disponible sur l'arbre du générateur
électrique s'exprime par :  ds  Lr I dr  MI ds
1
Pm  C p  2  R  / K  Vv  p  R 2  V dr  Lr I qr  MI qs
3
(5)
2
Copyright - 2019
Page 6
ISSN 1737-9334
Ls  I s  ms L’inductance cyclique du stator ; E. Modèle complet du système de conversion éolien

Le modèle complet de la chaine de conversion éolienne
Lr  I r  mr L’inductance cyclique du rotor. englobe en outre la turbine et son arbre, la MADA, un
redresseur MLI, un bus continu, un onduleur MLI, un filtre et
La forme générale du couple électromagnétique d’une un transformateur qui assure le couplage au réseau. Ce
machine asynchrone triphasée modélisée dans le repère de système est assez compliqué pour la génération des lois de
Park est donnée par la relation suivante : commande qui le gère, néanmoins il suffit de partager en
parties le système et ainsi suivre le flux de transfert de
Cam  P M sr iqs  idr  ids  iqr  (8) puissance de la turbine jusqu’au réseau. Pour simplifier le
système global un Modèle équivalent de la partie électrique
Puissances statoriques : dans le repère de Park sera construit.
Dans un repère biphasé quelconque, les puissances active III. SIMULATION DE LA CHAINE GLOBALE
et réactive statoriques d’une machine asynchrone s’écrivent Le schéma global d’une chaîne de conversion d’énergie
comme suit : éolienne connectée au réseau électrique est décrit par la figure
(3).
Ps   ds mr   qsiqs
Le système éolien a pour but de convertir l’énergie
Qs   dsiqs   qsids (9) cinétique du vent en énergie mécanique disponible sur un
arbre de transmission, puis en énergie électrique par
D. Principe de fonctionnement d’un redresseur MLI. l’intermédiaire d’une génératrice [8].
La structure de la cascade est celle représentée sur la figure
1. Les deux convertisseurs interposés entre le rotor de la
MADA et le réseau sont de type MLI à deux niveaux,
bidirectionnels en puissance. Les signaux de commande sont
déterminés en appliquant la commande vectorielle à la
MADA. Le convertisseur côté réseau aura le même modèle
que l’autre, sauf que les signaux de commande seront
déterminés en faisant un asservissement de la tension du bus
continu [6].
Dans notre cas, le convertisseur connecté au rotor Fig.3 : structure globale de la chaine de conversion
fonctionnera en redresseur, et celui connecté au réseau
fonctionnera en onduleur. Ce qui nous permet de récupérer la A. Le vent
puissance électrique disponible au niveau du bobinage Nous avons simulé le modèle de profil de vent sous forme
rotorique à travers les contacts glissants et de la réinjecter déterministe par une somme de plusieurs harmoniques. La
dans le réseau. Ce type de convertisseur peut opérer en figure 5 représente Le schéma bloc du vent étudié.
redresseur ou en onduleur. Quand le courant i_res est positif
(Opération redresseur), le condensateur C est déchargé, et le
signal d’erreur demande au bloc de commande plus d’énergie
à partir du réseau, le bloc de commande prend l’énergie
d’alimentation en produisant des signaux appropriés à
l’amorçage des transistors. De cette façon l’écoulement de
courant du côté alternatif vers le côté continu, ainsi, la tension
de condensateur est récupéré. Inversement, quand i_res
devient négatif (Opération Onduleur), le condensateur C est
surchargé, et le signal d’erreur demande à la commande la
décharge du condensateur et renvoyé l’énergie vers le réseau.
Fig. 4 Schéma bloc du vent
L’avantage de la commande MLI nous permet d’avoir une
bonne qualité de signal (formes sinusoïdales), ramenant le On rassemblant les blocs sous simulink, on a obtenu après
contenu harmonique vers des fréquences élevées et par la simulation, la figue 6 de la vitesse du vent en fonction du
suite la facilité de filtrage [7]. temps, qui est un signal aléatoire, obtenu après la sommation
de plusieurs signaux déterministes, tell les sinusoïdes des

ISSN 1737-9334
fréquences variables. Plus la vitesse du vent est plus

importante plus la tension issue de la génératrice est grande.
Fig. 8 Schéma bloc du modèle de la turbine éolienne (Matlab – Simulink).
Fig. 5 Profil du vent appliqué (m⁄s)
B. Pitch control
Pour faire apparaite le role du controle de l’angle de
calage,nous avons appliqué à la turbine un profil de vent varié
Fig. 9 Variation du Cp en fonction du temps
Fig. 6 Schémas bloc du pitch control

Fig. 10 Variation du couple produit en fonction du temps
Fig. 7 la variation de l’angle de calage avec le temps

Fig. 11 Variation de la puissance produite en fonction du temps.
C. Turbine
La turbine est entraînée par une MADA, dont le schéma
bloc correspondant à cette modélisation de la turbine est
représenté dans la figure 8.
Fig. 12 Variation de la vitesse du rotor.

ISSN 1737-9334
La figure (16) présente le schéma bloc de cette machine

présentée sous MATLAB-SIMULUNK .
Le système de commande assure l’optimisation de la
puissance extraite en maintenant le coefficient de la turbine à
sa valeur maximale [9].
D. Association de la MADA avec un onduleur à MLI
Fig. 156 Schéma bloc de la MADA en génératrice dans le référentiel (d-q)
Fig. 13 Génération des Signaux de commande MLI de l'onduleur
Le signal de sortie vaut 1 si la modulante est plus grande

que la porteuse 0, sinon ; le signal de sortie change donc
d’état à chaque intersection de la modulante et de la porteuse
Fig. 167 Courants statoriques
Fig 14 :représente les deux signaux à comparer (modulante et porteuse)
Fig. 178 Zoom sur les courants statoriques
Fig. 145 Tension Va associé à un onduleur MLI
Les courants circulants dans la charge sont parfaitement

sinusoïdaux (figure14) ainsi que les tensions composées
(figure15).Cela montre que ce modèle d’onduleur est
parfaitement adéquat pour notre application son contrôle
repose sur le signal de référence et de la porteuse.
E. Machine asynchrone à double alimentation

Après avoir présenté la modélisation de la machine
asynchrone à double alimentation dans un repère diphasé Fig. 189 Puissances réactives statoriques de la MADA
tournant (d-q), on fera maintenant la simulation d’une MADA

fonctionnant en génératrice [10].

ISSN 1737-9334
Les paramètres de la turbine éolienne utilisée :

Nombre de pale: 𝑁𝑝 = 3
Diamètre d’une pale ∶ 𝑅𝑇 = 3𝑚 Inertie: 𝐽 = 315 𝐾𝑔. 𝑚2
Les paramètres de la turbine éolienne utilisée : Nombre de
pale: 𝑁𝑝 = 3
Diamètre d’une pale ∶ 𝑅𝑇 = 3𝑚 Inertie: 𝐽 = 315 𝐾𝑔. 𝑚2
Paramètre : Résistance du stator :𝑠 = 1.2 𝛺
Résistance du rotor : 𝑅𝑠 = 1.8 𝛺
Inductance stator ∶ 𝐿𝑠 = 0.1554 𝐻
Inductance du rotor :𝐿𝑟 = 0.1568 𝐻
Inductance mutuelle ∶ 𝑀 = 0.15 𝐻
Fig. 20 : Puissances actives statoriques de la MADA Constantes mécaniques : Moment Inertie : 𝐽 = 0.2 𝐾𝑔. 𝑚2
Coefficient de frottement : 𝑓 = 0.001 𝑁. 𝑚. 𝑆/𝑟𝑎𝑑
REFERENCES
[1] Rapport, ‘Le Baromètre Eolien; Systèmes Solaires, Journal des
Energies Renouvelables’, N°183, Union Européen, Février 2008.
[2] F. Poitiers, ‘Etude et Commande de Génératrices Asynchrones pour

l’Utilisation de l’Energie Eolienne’, Thèse de Doctorat, Ecole
Polytechnique de l’Université de Nantes, Décembre 2003.
[3] G. Dreyfus, J. Martinez, M. Samuelides, M.B. Gordon, F. Badran et S.

Thiria, ‘Réseaux de Neurones: Méthodologie et Applications’,
Editions Eyrolles, 2002.
Fig. 191 : variations de la vitesse de rotation
[4] A. Boyette, Contrôle-Commande d’un Générateur Asynchrone à
La figure (17) montre que le courant de la génératrice Double Alimentation avec Système de Stockage pour la Production
asynchrone à double alimentation en phase avec la tension. Eolienne’, Thèse de Doctorat, Université Henri Poincaré, Nancy I,
L’amplitude ainsi que la fréquence s’adaptent malgré la Décembre 2006.
variation de la vitesse.
[5] G. Tapia, A. Tapia and J.X. Ostolaza, ‘Proportional-Integral
Regulator-Based Approach to Wind Farm Reactive Power
On observe une dynamique qui réagit rapidement pour la Management for Secondary Voltage Control’, IEEE Transactions on
puissance active et réactive du stator.la puissance active du Energy Conversion, Vol. 22, N°2, pp. 488 - 498, June 2007.
côté statorique est négative après l’instant t=0.7s ce qui
[6] A. Tapia, G. Tapia, J.X. Ostolaza and J.R. Saenz, ‘Modeling and
signifie que le réseau dans ce cas est un récepteur de l’énergie
Control of a Wind Turbine Driven Doubly Fed Induction Generator’,
fournie par la MADA, la puissance réactive est négligeable, IEEE Transactions on Energy Conversion, Vol. 18, N°2, pp. 194 -
la vitesse de référence et celle mesurée sont tellement 204, June 2003.
confondues qu’il impossible de les différencier et donc il y’a
[7] Y. Djeriri, ‘Commande Vectorielle d’une MADA Intégrée à un
un suivi parfait de la référence. Système Eolien’, Mémoire de Magister, Université de Sidi Bel-Abbes,
Décembre 2009.
IV. CONCLUSION
[8] D. Yathley, ‘Back Propagation Neural Networks for Non-Linear Self
Cette étude a montré la possibilité de contrôle de la Tuning Adaptive Control’, IEEE, Control Systems Magazine, Vol. 10,
machine par les courants rotoriques, afin de déterminer la N°3, pp 44 - 48, January 1990.
référence de la composante directe et en quadrature du
courant rotorique à partir des puissances active et réactive. [9] A. Meroufel*, Y. Djeriri†, A. Massoum et A.
Les résultats montrent que le courant de sortie en phase avec Hammoumi Commande vectorielle par les réseaux de neurones
artificiels de l’énergie d’une MADA intégrée à un système éolien
la tension d’alimentation, en décrivant la bonne performance
Revue des Energies Renouvelables Vol. 13 N°4 (2010) 669 –
du système en boucle fermée en termes de suivi de trajectoire 682
Paramètres du système [10] BOUDJEMA, Z , TALEB, R , DJERIRI, Y , YAHDOU, A . "A novel
direct torque control using second order continuous sliding mode of a
Puissance nominale : 𝑃𝑛 = 4 𝑘𝑊 doubly fed induction generator for a wind energy conversion system".
Tension nominale : 𝑣/𝑈 = 220/380 𝑉 − 50 𝐻𝑧 Turkish Journal of Electrical Engineering and Computer Science 25
La vitesse nominale : Ω𝑛 = 1440 𝑡𝑟/𝑚𝑖𝑛. (2017): 965-975
Nombre de pairs de pôles : 𝑃 = 2

ISSN 1737-9334
LE BIOGAZ SOURCE D’ENERGIE

RENOUVELABLES
Cas des déchets ménagers « biodéchet »
BAICHATA MALIKA Pr TAMALI MOHAMED
Departement de genie civil et hydraulique Department genie electriques
Université de TAHRI MOHAMED
BECHAR Université de TAHRI MOHAMED
rmalika.2007@gmail.com BECHAR
Resume- Cette dégradation biologique implique une

La conversion biologique de la biomasse en méthane a reçu une microflore spécialisée et diversifiée qui exige des
attention croissante au cours des dernières années, la biomasse conditions bio-physico-chimiques spécifiques et adaptées.
végétale ou biodéchet a été explorée pour leur potentiel de Comme la totalité des réactions biologiques, les réactions
digestion anaérobie en méthane et leur productivité ; la impliquées dans la dégradation anaérobie sont réalisées en
caractérisation de ce déchet solide a montré que cette biomasse est présence d’eau, c'est-à-dire en milieu aqueux [4].
riche en matière organique 64.4% avec une ratio optimal de
carbone/azote de 11.4,L’objectif de cette étude est de contribuer à La production de biogaz en Algérie Historiquement la
la diminution de la quantité de déchets solides transporter à la production du biogaz en Algérie remonte à l’ère coloniale. En
décharge, par une digestion anaérobie humide discontinue type " effet, les colons industriel sont utilisé le biogaz issu de la
batch" ,la quantité ainsi que la qualité de biogaz produit est fermentation des déchets de bovins, conçu dans les années 50
estimé a l’échelle de laboratoire. par l’Institut National Agronomique (INA) afin de faire
fonctionner leurs véhicules agricoles.
Actuellement et selon le Ministère de l’Aménagement de
Mots clés biomasse,digestion anaérobie ,humide: biodéchet , Territoire et d’Environnement (MATE) l’Algérie compte
biogaz.
installer plus de 50 centres d’enfouissement technique (CET) à
travers l’Algérie d’ici à 2025. Ceci va permettre le
développement du service public de gestion des déchets
I. INTRODUCTION ménagers dotés par des technologies de la valorisation des
déchets et plus particulièrement la récupération de la fraction
Lorsqu’il est question d’énergies renouvelables, on biodégradable pour la production du biogaz. En ce qui
évoque le plus souvent les éoliennes, les capteurs solaires concerne la production actuelle du biogaz en Algérie aucune
ou encore l’énergie hydraulique. Cependant, il existe estimation n’a été faite à titre principal .Néanmoins, certains
d’autres solutions comme par exemple la production ont été faits à titre d’expérimentation par le Centre de
d’énergie à partir de la biomasse :le bois, les biocarburants Développement des Énergies renouvelables (CDER) afin de
ou encore le biogaz. [1] mettre en évidence la faisabilité technique de ces projets en
Algérie [5].
La population des pays en développement est en pleine
expansion, entrainant une augmentation des besoins Etapes de digestion anaérobie La technique utilisée dans
alimentaires et de la production des déchets et effluents. les systèmes anaérobie s’appelle la digestion anaérobie ou
Le recyclage de déchets participe à la lutte contre la méthanisation. C’est un processus complexe impliquant
pauvreté et à l’assainissement environnemental. plusieurs réactions biologiques. On le décrit généralement
par quatre étapes principales, impliquant quatre consortiums
La méthanisation est un processus de digestion microbiens spécifiques interagissant les uns avec les autres.
anaérobie permettant d’atteindre généralement un double La digestion anaérobie se déroule selon une succession de
objectif de valorisation énergétique par récupération de quatre étapes distinctes, De cette façon, les molécules
méthane (CH4) et de stabilisation des déchets organiques complexes sont métabolisées jusqu’à leur minéralisation,
en vue d’une valorisation matière par sa restitution c’est-à-dire leur transformation en méthane et en CO2. Ces
partielle au sol [3] étapes sont : l’hydrolyse, l’acidogénèse, l’acétogénèse et la
méthanogénèse [6].
Copyright - 2019
Page 11
ISSN 1737-9334
La figure présente l’ensemble des voies métaboliques, avec L’humidité L’eau est considérée comme un facteur
leurs groupes de micro-organismes respectifs, en trois étapes déterminant dans le processus biochimique de dégradation
principales : (1) l’hydrolyse et acidogenèse ; (2) l’acétogenèse anaérobie; les taux d’humidité optimale oscillent selon les
et (3) la méthanogenèse. études, entre 40 et 80 % selon les auteurs d’après [10].
La ratio C/N : Le rendement en méthane d’un substrat

dépend de sa composition chimique .Le rapport C/N (C
non ligneux) du substrat fermenté détermine la
proposition de CH4 dans le biogaz : de l’ordre de 65% de
CH4 pour un rapport C/N =8 ; de l’ordre de 52 %de CH4
pour un rapport C/N =50 [11].
Les prétraitements Les traitements mécaniques (cribles,

tamis, trommel, broyage) permettent de séparer la matière
organique des résidus inertes : verre, plastiques, métaux,
minéraux ainsi que d’autres impuretés. Ces traitements
permettent tout d’abord d’assurer le bon fonctionnement
des réacteurs. De plus, ils ont pour effet de concentrer la
matière organique et aussi de réduire la taille des
particules, améliorant ainsi l’efficacité du procédé de
méthanisation [12-15].
l’inoculum:L’ajout d’un inoculum d’ensemencement lors

du démarrage des tests de biométhanisation favorise un
Fig. 1 Principales étapes de la digestion anaérobie et leurs groupes de état d’équilibre de l’ensemble du processus de digestion
microorganismes respectifs. anaérobie.
II. Matériels et Méthode

Les paramètres d’influence de la digestion anaérobie Estimation des résidus végétaux de restaurant 8 mai 1945 de
Anaérobioses : En l’absence d’oxygène (O2), des bactéries université de BECHAR :
dégradent partiellement la matière organique (MO) [7,8]. L’estimation de bio déchet au niveau de restaurant de la
résidence 8mai 1945 pour un mois de suivi a été révéler que :
Les variations de température L’intensité de l’activité des
La matière végétale (légumes et fruits) journalière entrante est
micro-organismes, dont dépend la production de méthane
estimé de 150Kg à 180 kg, pour 3000-3100 étudiants
est en fonction de la température du milieu.
bénéficiant du repas, cette quantité de la matière végétale
Négligeable à 0°C, elle devient exploitable dés 15°C pour produit un résidu estimé de 40 Kg à 50 Kg par jour
augmenter rapidement a partir de 20°C avec un maximum II. prelevement de biodéchet
à 37-40°C pour les fermentations mésophiles (25-40°C) et Le substrat d’alimentation de digesteur type batch utilisé est
a 55°C pour les fermentations thermophiles. la fraction fermentescible des ordures ménagères : déchets
organiques biodégradables, ou biodéchets (déchets de cuisine
Le pH et TAC : Il est indispensable de maintenir le pH et restauration) ; elles sont récupérés sur le lieu de prés
dans la gamme de la neutralité (6.7 – 7.4), l’optimum étant collecte sélectif des restaurations de université de Béchar.
autour de 7.0 – 7.2 unités pH. Afin de maintenir le II. 2 Prétraitement mécanique de substrat (bio déchet)
réacteur au pH optimal, celui-ci est régulé par l’ajout de Le bio-déchet récupérés est séché a l’air pendant
saude ou de bicarbonate de sodium. Alcalinité TAC : une nuit puis suivi d’un traitement mécanique pour
l’alcalinité est la capacité d’un milieu à absorber un apport réduire la granulométrie (la taille des particules
de protons. Généralement, on considère qu’il faut une organiques de 2cm a 5cm) au moyen d’un broyeur
alcalinité en bicarbonate de 1000 a 3000 mg /l exprimée électrique puis homogénéiser (figure1) enfin le biodéchet
en CaCO3. est humidifié, cette humidité est acquise grâce à l’ajout de
l’eau. Le substrat est inoculé à l’aide d’une quantité des
AGV : Le pH est essentiellement lié à la présence bouses de vache ; au milieu de la fermentation le pH bas
d’acides gras volatils. Lors du bon fonctionnement du est corrigé par l’adition de la chaux éteinte local afin que
digesteur, la variation du pH est tamponnée par la la fermentation méthanique ne soit pas inhiber.
présence des bicarbonates produits par les méthanogènes
[9].
Copyright - 2019
Page 12
ISSN 1737-9334
Détermination de carbone total (COT)
La détermination du Carbone Organique Total

(COT) nous permet d’évaluer la quantité de matière
organique présente dans des échantillons solides. La
détermination de carbone total est basé sur méthode
Walkley-Black modifiée par dosage du carbone
organique dans les sols agricoles de code MA. 1010 –
WB 1.0
bio déchet avant le broyage bio déchet après le broyage
Le contrôle du pH a été effectué les premiers jours
de la fermentation, la chute du pH a été équilibrée en
figure1 Prétraitement de bio déchet
chaulant le digesteur par une solution de chaux éteinte
Fig.3 digesteur de capacité de 120L commercial local.
Détermination de La matière sèche (MS) La température : la température maintenu à 37ºC dans le

domaine mésophile.
La matière sèche (MS) ou matières volatiles sèches
(MVS) : Il s’agit principalement de la détermination de la Description prototype de digesteur au niveau de
matière sèche (MS) par la méthode empirique est laboratoire
généralement basée sur le principe de quantification
gravimétrique (séchage a 105°c pendant 24heures) Norme L’expérience a été réalisé dans deux essais ; dans chaque
NF XP X30-408 essais un erlenmeyer est remplie de substrat homogénéisé de
tel façon que le (2/3) est remplié et le 1/3 de volume reste
Determination de La masse volumique apparente (ρ) vide, on assure l’anaérobiose par un bouchon de silicone muni
La masse volumique apparente (ρ) correspond à la masse d’un trou dans le centre. Le premier erlenmyer est raccorder
de l’unité de volume à l’état sec. La mesure de ‘ρ’ Elle est par un tuyau vers une chambre à air pour faire un test
déterminée en appliquant la formule suivante ρ=m/V d’inflammabilité et l’analyse qualitatif de biogaz produit, le
deuxième erlenmyer est raccordé par un tuyau vers une
Détermination du pH, de la conductivité La mesure du pH et
éprouvette graduéé remplie d’eau acidifier pour suivre la
la conductivité est réalisée selon la norme internationale NF
EN 27888. production journalière de biogaz cumulé fig 2
Le pH est détermines par la méthode potentiométrique d’une On a procédé une agitation magnétique de deux heure
électrode combinée (Consort 861) chaque jours afin d’homogénéiser le substrat.
La Conductivité électrique (λ) est la mesure de la
concentration des ions solubles afin d’apprécier la salinité du
substrat . Elle est mesurée par le conductimètre et elle est
exprimée en (mS/cm).
Détermination de l’azote total (azote de kjeldahl) (NTK)

Dosage de l'Azote Organique par la Méthode
KJELDAHL consiste à transformer l’azote organique en azote
minéral (minéralisation), puis à déplacer l'ammoniac du sel
d'ammonium obtenu (distillation) puis le neutraliser par une
solution acide de titre connu (dosage).
Détermination de la matière organique(MO)
Figure 9: Description prototype au niveau da laboratoire
La détermination de la matière organique (la matière
volatile (MV)) permet la quantification de la matière III. RÉSULTATS ET DISCUSSION
organique dans les déchets ; une méthode gravimétrique
également basée sur la perte de masse d’un échantillon sec
(échantillon issu de la détermination de la MS) après
calcination dans un four à moufle à 550°C pendant 4
heures Norme NF U 44-160.

ISSN 1737-9334
TABLEAU 1
CARACTERISATION PHYSICO-CHIMIQUE DE BIO DECHET (RESIDUS VEGETALES)
Parameters Résultats Unité

MS 21 %
Matière Minéral (MM) 35.6 %

MO 64.4 %
C (COT) 37.4 %/MS
N (NTK) 3.33 %/MS
Ρ 370 Kg/m3/MS
Λ 13500 µs/cm/MS
pH 6.49 /MS
C/N 11.4 Graphe1 : cinétique de la Production de biogaz cumulé
D’après les résultats de caractérisation physico-

chimique de substrat (biodéchet) est riche en matière A partir de graphe on peut distinguer trois phases :
organique et les sels minéraux ( MO 64.4%) avec une
phase 1 : c’est la phase latence :
teneur en l’azote de 33,2 mg/g, avec un faible taux de
l’azote total. avec une ratio optimal de - Les bactéries sont adapté au nouveau milieu et se
carbone/azote de11.4. le test d’inflammabilité est positif développe ; le temps est de 24 h pour une production
figure croissante en fonction du temps
phase 2 :c’est la phase de croissance exponentielle :
- Bactéries mésophiles ont une activité métabolique

optimale. La production de biogaz accumulé est régulière
avec un taux de 90 ml /h
Phase 3 :c’est la phase déclin :
Arrêt de croissance, l’épuisement de substrat la production

est on décroissance jusques l’arrêt total.
On estimé Vtotal=6153ml durant six jours, donc le temps de
séjour est de six jour avec un taux de production moyen de
Figure 22: test d’inflammabilité biogaz est 150ml /MS
l’analyse de biogaz dans la chambre a air est effectue a l’aide
d’un analyseur (Geotech GEM 5000) ; l’analyse de biogaz
cumulé a dévoilé que le biogaz contient 40 % de CH4 et 59% IV CONCLUSION
de CO2 donc nous recommandons des filtres d’adsorption pour La biomasse est une des plus abondantes sources de
une réduire la quantité de CO2 dégagé. carbones renouvelables et aura un rôle prépondérant à jouer
dans la production d’énergie dans le futur.
Le cinétique de la production de la quantité de biogaz cumulé La technologie du biogaz repose sur un concept de
est suivi par la lecture directe de eau verse de l’éprouvette recyclage des déchets lors duquel la masse organique est
graduée équivalente a la quantité de biogaz dégagé en fonction convertie par fermentation méthanique en énergie (méthane :
du temps, donc d’après ses résultats on obtient le graphe CH4) simple, propre et prête à satisfaire un grand nombre
suivant : de besoins énergétiques de manière totalement décentralisée
[21].
La biométhanisation pour le traitement et la valorisation énerg
étique des déchets
Le potentiel de production de biogaz a partir de biomasse

local et avec un digesteur familial peut être une source de

ISSN 1737-9334
l’énergie durable de la société mais il reste que la fermentation durable, de l’Environnement et de la Lutte contre les changements
climatiques du Québec, 2016, 12 p.
anaérobie dépend de plusieurs paramètres tel que la source de
substrat, la température, l’agitation…. etc
REFERENCES
[1] Bouhadiba B,Mezouri F,KehilaY,Matejka G(2010) For an integrated
Management for Municipal solid wast in Algeria.Systemic and
Methodological Approaches.International Review of chemical
Engineering- rapid Communication 2 (3):426-429.
[2] Brahim Djemaci . Recyclage des déchets à travers un système de
consigne : Cas des bouteilles en plastique en Algérie ,Université de
Rouen, Centre d’analyse et de recherche en économie (CARE), France
[3] Joseph o., Rouez m., Métivier-Pignon h., Bayard r., Emmanuel e.,
Gourdon r. Adsorption of heavy metals on to sugar cane bagasse:
Improvement of adsorption capacities due to anaerobic degradation of
the biosorbent. Environmental Technology, 2009, vol. 30, n° 13, pp.
1371-1379.
[4] Moletta,R. et Cansell,F.,(2003).Méthanisation des déchets
organiques :Etude bibliographique,RECORD :193.
[5] Laskri N, Hamdaoui O, Nedjah N.,Traitement et valorization des
déchets par procédé de digestion anaérobie :production du biogaz. Revue
des Energie renouvelables CER’07 Oujda (2007) 23-26.
[6] Godon,J.J.(2008).Aspects biochimiques et microbiologiques de la
méthanisation. In Moletta,R.,La méthanisation (chp.3,p.61-85) .Paris, Les
Editions Eugen Ulmer,120 p.
[7] Chynowetha D.P., Owensa J.M., Legrand R. (2002). Renewable methane
from anaerobic digestion of biomass. Renewable Energy 22, p. 1-8.
[8] Karagiannidis,A., Perkoulidis,G. (2009).A multi-criteria ranking of
different technologie for the anaerobic digestion for energy recovery of
the organic fraction of municipal solid wastes.Bioesource Technology
100 (2009) 2355-2360.
[9] Bitton,G., (1994).wastewater microbiology,Wiley-Liss Editor ,New York.
[10] [Aguilar-Juarez,O. Analyse et modélisation des réaction biologiques
aérobies au cours de la phase d’exploitation d’un casier d’un centre
d’enfouissement technique.These de doctorat spécialité Génie des
Procédés,INSA de Toulouse,2000 ,233 p.
[11] Bertran de La Farge,Le biogaz : Procédés de fermentation méthanique,
Masson, Toulouse,1995, 80 pages ,
[12] Angelidaki I., Ahring B. K. Methods for increasing the biogas
potential from the recalcitrant organic matter contained in manure.
Water Science and Technology, 2000, vol. 41, pp. 189-194.
[13] Forster-Carneiro, T., Perez, M. et Romero, L.I., (2008), Anaerobic
digestion of municipal solid wastes: Dry thermophilic performance.
Bioresource Technology 99: 8181-8184.
[14] Izumi, K., Okishio, Y.K., Nagao, N., Niwa, C., Yamamoto, S. et Toda,
T., (2010),Effects of particle size on anaerobic digestion of food
waste. International Biodeterioration and Biodegradation 64: 601-608.
[15] Pommier, S., Manas Llamas, A. et Lefebvre, X., (2010), Analysis of the
outcome of shredding pretreatment on the anaerobic biodegradability
of paper and cardboard materials. Bioresource Technology 101: 463-
468.
[16] Benkadour O.,Ghezal K. (2007).Contribution à l’étude hydrogéologique
de la nappe du Turonien de Ouakda-Bechar (SW algérien),CU
Bechar,80 p.
[17] Chawakitchareon P., (1990).Etude descriptives de décharge de déchets
et recherche de critères permettant d’estimer le potentiel méthanique des
ordures ménagères et/ou d’orienter leur traitement. Thèse de DEA de
l’INSA,Toulouse,Université de Chulalongkom,Bankok,Thailande.
[18] Degrémont , ‘Mémento Technique de l’Eau’, 10ème Ed, Tome 1.
Lavoisier, p.535- 536, 2005
[19] Makhoukh M , Sbaa M ,Berrahou A , Van.Clooster M .Contribution a
l’etude physico-chimique des eaux superficielles de l’oued Moulouya
(Maroc Oriental) Larhyss journal ,2011 ; 09 pp.149-169
[20] Centre d’expertise en analyse environnementale du Quebec.
Détermination de la demande chimique en oxygène : méthode de reflux
en système fermé suivi d'un dosage par colorimétrie avec le bichromate
de potassium, MA. 315 – DCO 1.1, Rév. 4, Ministère du Développement
Copyright - 2019
ISSN 1737-9334 Page 15
Emissions BF des chargeurs commandés par

régulateur à Hystérésis embarqués dans les
véhicules électriques hybrides rechargeables
(V.E.H.R)
L. BESSOLTANE#, Z. AZZOUZ*, D. TOUMI#
#
Faculté de Sciences Appliquées, Université Ibn Khaldoun, Tiaret.
bess_oltan@yahoo.fr,toumi_dj@yahoo.fr
*
Laboratoire de développement des entrainements électriques -Université des Sciences et de la Technologie d’Oran «MB»
zinazzouz@yahoo.fr
Résumé— Nous abordons dans ce papier l’étude CEM spécialistes dans le domaine du génie électrique notamment
(compatibilité électromagnétique) des chargeurs de batteries des l’électrotechnique.
véhicules électriques hybrides rechargeables (VEHR). Aussi, D’autre part, pour pallier aux problèmes de pollution
nous intéressons aux composantes basses fréquences (BF) des électromagnétique (EM) liés à l’utilisation des équipements
courants harmoniques absorbés par ces dispositifs. Nous connectés au réseau de distribution, des organismes de
abordons également dans cet article les structures des chargeurs
normalisation internationaux ont été institués afin de préserver
(ainsi que leurs stratégies de commande) les plus adoptées dans
les VEH. Une attention particulière (du point de vue CEM) est la qualité d’énergie électrique. Ainsi, la CEI et l’IEEE
réservée aux chargeurs à MLI triphasés à deux niveaux et NPC à (« Institute of Electrical and Electronics Engineers ») sont les
trois niveaux, utilisant la technique de contrôle par hystérésis des deux principaux organismes de normalisation internationaux
courants prélevés sur le réseau avec un réglage classique (PI). dans le domaine de l’électrotechnique. Ils s’occupent de
L’analyse CEM que nous proposons est basée sur la comparaison l’activité de normalisation CEM au niveau mondial [6] [7].
des amplitudes des courants harmoniques des structures étudiées Les chargeurs des batteries des VEHR, connectés
par rapport aux valeurs autorisées par les normes CEM en directement au réseau électrique, absorbent un courant
vigueur. A l’issue de cette comparaison, nous préconisons des supérieur à 16 [1] [4], ce qui nous ramène dans ce cas à
solutions en faveur d’une meilleure CEM.
la norme CEI 61000-3-12. En effet, cette dernière est
MT ayant un courant absorbé supérieur à 16

Mots clés— CEM, VEHR, chargeur à deux niveaux, chargeur à applicable aux appareils connectés aux réseaux publics BT et
ou égal à 75 A
trois niveaux, commande MLI, régulateurs à Hystérésis. et inferieur
par phase [6]. Cette norme décrit les limites
I. INTRODUCTION d’émissions autorisées des courants harmoniques, par le taux
Les technologies de VEH se divisent en trois catégories, de distorsion harmonique (THD) et le rapport de distorsion
dépendant du niveau d’assistance du moteur électrique. On harmonique pondérée partielle (PWHD: «Partial Weighted
parle alors de chaînes de traction Micro-hybride, Mild-hybride Harmonic Distortion ») [8] [7] [2]:
et Full-hybride.
= ℎ
Le véhicule électrique hybride rechargeable (VEHR) est une
variante de la technologie Full-hybride [1]. La nouveauté de
cette technologie est de permettre la recharge de la batterie sur ! "
: Amplitude de l’harmonique au rang ℎ.

une source électrique via des chargeurs. Ces derniers doivent (1)
répondre aux normes CEM qui limitent les interférences
électromagnétiques (I.E.M), l’injection d’harmoniques et de : Courant de court-circuit
courants continus dans les réseaux électriques [2] [3]. D’autre Le PWHD est utilisé pour fournir une marge permettant à
part, ils doivent s’assurer un courant absorbé avec le minimum certaines harmoniques d’ordre élevé de dépasser les limites
d’impact sur la qualité d’énergie et un facteur de puissance spécifiées, tout en garantissant que l’effet des harmoniques
élevé afin de maximiser la puissance réelle disponible [4]. d’ordre élevé ne dépasse pas les limites acceptables:
Le tableau I représente les limites d’émissions de courants
II. CONTRAINTES NORMATIVES ET CEM harmoniques des équipements BT triphasés équilibrés
Au niveau international, des travaux de normalisation du spécifiés par la norme CEI 61000-3-12 classe A [2].
moins doivent être inférieures à 16/ %. Les harmoniques

VEH ont été entrepris depuis de nombreuses au sein des Les valeurs relatives d’harmoniques paires d’ordre 12 ou
comités de la Commission Electrotechnique Internationale
(CEI) [1] qui englobent les constructeurs automobile et les paires de plus d’ordre 12 sont prises en compte dans les

ISSN 1737-9334
valeurs de THD et PWHD comme pour les harmoniques • la structure de redresseur suivi d’un hacheur Boost
impaires [2]. (Correction du facteur de puissance ‘PFC’).
TABLEAU I • la structure de redresseurs à injection de courant.
LIMITES D’EMISSIONS DE COURANTS HARMONIQUES • le redresseur de VIENNA.
SPECIFIES PAR LA NORME CEI 61000-3-12 • le redresseur à Modulation de Largeur d’Impulsions
PQRS d’harmonique TU ⁄TV %

Minimum Courant acceptable individuel Taux de courant maximal (MLI).
Les topologies des chargeurs les plus répandues dans les
?X ?Y ? ?
d’harmonique %
"
‘Short circuit
ratio’ THD PWHD applications à prélèvement sinusoïdal sont les redresseurs de
10. tension à MLI. Leur réversibilité en puissance les rend très
33 7.2 3.1 2 13 22
7 populaires surtout dans les VEHR car ils peuvent être utilisés
66 14 9 5 3 16 25 en redresseur dans le mode recharge de batterie et en onduleur
120 19 12 7 4 22 28 dans le mode de traction [9] [11] [12].
Plusieurs logiques de commande sont envisageables pour cette
250 31 20 12 7 37 38
≥ 350
topologie. Nous nous sommes intéressés dans ce travail à la
40 25 15 10 48 46 commande par hystérésis d’un chargeur à MLI triphasé à deux
niveaux et NPC à trois niveaux.
III. CHARGEURS EMBARQUES DANS LES VEHR
Les chargeurs de batterie peuvent être classés selon leurs IV.1. CHARGEUR A DEUX NIVEAUX COMMANDE PAR
emplacement par rapport au véhicule électrique (VE) : UN REGULATEUR A HYSTERESIS
embarqués (installés à l'intérieur du VE) ou non embarqués La structure globale d’un chargeur triphasé, à MLI deux
(installés à l'extérieur du VE) [5] [9] [4]. niveaux à contrôle par hystérésis des courants prélevés sur le
Il faut savoir, que depuis les années 90 les chercheurs se réseau par un pont redresseur de tension, est présentée sur la
sont orientés vers des chargeurs embarqués qui utilisent Fig. 3 [11] [13] [14].
)
l’onduleur de traction et les enroulements du moteur, ce qui
= - 1
~
permet de s’acquitter des inductances de filtrage
+#
#
()
supplémentaires en entrée afin de réduire le nombre de
~ +$
composants du système. Cependant cette structure présente $
R
~
C
+
des contraintes supplémentaires en termes de CEM [1] [10].
*#$
Nous avons fait le choix dans ce travail, de présenter
uniquement des architectures des chargeurs triphasés
bidirectionnels embarqués. La Fig. 1 montre un schéma M& PLL
simple d’une architecture parallèle d’un VEHR en intégrant le
O#
Triphasé
?#∗
moteur de traction dans le chargeur de batterie [9] [10].
? ∗
#D sin%M&' O$
?$∗
-
? ∗ #D sin%M& − 120' O
+
?∗
Embrayage -
? ∗ #D sin%M& − 240'
+
Réservoir Moteur -
Thermique +
COUPLEUR
de carburant
? ∗ #D √2
?)∗ Régulateur - + ()
∗
×
3. GHII =
Différentiel
de ()
Transmission
∗ ∗
) BC
Convertisseur Moteur
Batterie AC/DC et DC/AC
Electrique
bidirectionnel
Embrayage Fig. 2 Structure globale d’un chargeur triphasé à 2 niveaux commandé
Les trois courants #$ %&' sont contrôlés à l’aide de trois

par régulateur à Hystérésis
Fig. 1 Architecture d’un VEHR intégrant le moteur de traction dans le comparateurs à hystérésis à bande fixe. L’ensemble des trois
chargeur de batterie
sorties de ces comparateurs détermine les ordres de
On distingue deux types de mode de charge à savoir la commande des interrupteurs constituant le pont.
charge lente pour le cas d’un réseau monophasé et la charge Cette structure comprend deux boucles de régulation en
rapide pour un réseau triphasé. Dans le cas d’une alimentation cascade: la boucle externe qui contrôle la tension '() ' et la
boucle interne qui contrôle les courants triphasés ' #$ '.
Les tensions *#$ %&' permettent d’écrire la relation entre les
par un réseau monophasé, l’idée est d’utiliser deux bras de
courants #$ %&' et les tensions de phase +#$ %&' [11]:

l'onduleur de traction côté réseau en redresseur. Le troisième
bras de l'onduleur fonctionne en hacheur série, l’enroulement
+# %&' = *# %&' − - / − 1 # %&' %2'
du moteur est utilisé en tant qu'inductance de filtrage [1] [5]. ). %0'
)0
1
# %&' = 5* %&' − +# %&'6 %3'
- 3+1 #
IV. TOPOLOGIE ET COMMANDE D’UN CHARGEUR
TRIPHASE A PRELEVEMENT SINUSOÏDAL
Plusieurs structures des convertisseurs AC/DC à Le courant à la sortie du pont est exprimé par [5]:
=8 + %4'
prélèvement d’un courant sinusoïdal ont été développées et )9:;
présentées dans la littérature [11] [12] [13] à savoir : ) )0

ISSN 1737-9334
() = %() − ' %5'

D’où : TABLEAU II
PARAMETRES DE SIMULATION
C
85 Ω
Paramètres et Symboles Valeurs
2000gh
courant #$ et de la tension () respectivement en utilisant
Résistance de la charge (1 )
Les figures 3-a et 3-b montrent le principe de contrôle du
0.3 Ω
Capacité du bus continu (8)
Résistance du filtre de raccordement (1I )
*# Inductance du filtre de raccordement (-I ) 14 i
les équations (3) et (5) [15] [12].
50 k
Régulateur de +#
∗ ∗
1
220 =
# #
Fréquence de la tenion alternative (j)
-I . O + 1I
Valeur efficace de la tenion alternative (GHII )
0.8672
a)
+- #
+-
169.375
#
Paramètres du régulateurs IP (lm ), (ln )
Tension de référence (()∗ ) 620 =
∗
#
∗
Régulateur de ) ) 1 () Limites de la bande hystérésis (±∆?) 0.03
+- ( +- 8. 3
1
tension () , le courant prélevé par le chargeur et son spectre

b)
Les figures 6, 7 et 8 illustrent respectivement la
a) du courant #$ b) de la tension ()
Fig. 3 Structure et principe du contrôle de contrôle
fréquentiel.
La commande par hystérésis consiste à contrôler le courant
fixe à deux niveaux, de largeur [? enveloppant sa référence.

prélevé par le redresseur pour le maintenir dans une bande,
Chaque violation de cette bande donne un ordre de

∗
commutation aux interrupteurs. Son principe (Fig. 4) consiste
à établir la différence (\#$ ) entre le courant de référence #$
et le courant absorbé par le chargeur #$ _ H . Cette erreur est
Fig. 6 Tension du bus continu (structure deux niveaux)
ensuite comparée à une bande d’hystérésis afin de fixer les

ordres de commande des interrupteurs [7] [14].
O
∗
#$
vwxR
#$ _ H
− ∆I + ∆I Fig. 7 Courant absorbé par le chargeur (structure deux niveaux)
Fig. 4 Principe de contrôle du courant par hystérésis
pour un chargeur à 2 niveaux
Si ^∗ − ^_ H > +∆? ⇒ O^ = 0 , (k les phases a, b ou c)

L’algorithme de cette stratégie est le suivant [14]:
Si ^∗ − _ H < −∆? ⇒ O^ = 1
Si −∆? < ^∗ − ^_ H < +∆? et d% ^∗ − ^_ H '/d& < 0
⇒ O^ = 1
Si −∆? < ^∗ − ^_ H < +∆? et d% ^∗ − ^_ H '/d& > 0
⇒ O^ = 0
Le rôle de la boucle de régulation de la tension () est de
Fig. 8 Spectre d’harmoniques du courant absorbé par le
chargeur (structure à deux niveaux)
Un régulateur de type ? est souvent utilisé pour la

maintenir cette tension à une valeur de référence constante. IV.2. CHARGEUR MLI NPC A 3 NIVEAUX, COMMANDE
régulation de tension ( ) . Le principe de cette boucle de

PAR UN REGULATEUR A HYSTERESIS
Le chargeur à MLI NPC triphasé à trois niveaux est la
régulation est illustré à la Fig. 5 [11] [12] [13] [15]. structure de chargeur la plus répandue dans le marché des
()∗ ∗
√2 ? #D
∗
? #D √2 1 ()
VEH. Le principe de ces structures est illustré à la Fig. 9 [1]
Régulateur ) ) )
de ( 3GHII 3GHII 8. 3
q#
+- +- [16].
()
1
()
q#r 2
Boucle de
= - 1
~
courant
~ i ()
st&&*u *
~
Fig. 5 Schéma bloc de la boucle de régulation de la tension du
bus continu
A noter que pour mettre en évidence cette technique dans q# ()
2
un contexte CEM, nous avons utilisé MATLAB/Simulink. Les q#r
paramètres de simulation sont consignés dans le tableau II.
Fig. 9 Chargeur MLI NPC triphasé à 3 niveaux

ISSN 1737-9334
Par ailleurs, dans la littérature plusieurs stratégies de fréquentielles du courant absorbé par le chargeur de batterie
commandes du redresseur MLI NPC à 3 niveaux sont montre que le courant absorbé est proche du sinus et que son
développées. Aussi, dans ce travail nous avons adopté la THD est égal à 0.83% pour la structure de 2 niveaux et à
2.65% pour la structure NPC à 3 niveaux ce qui répond
technique de contrôle du courant prélevé par un régulateur à
parfaitement aux exigences de la norme internationale CEI
hystérésis. Cette commande consiste à contrôler le courant 61000-3-12.
largeur [? (Fig. 10) [14].

pour le maintenir dans une bande comporte trois niveaux de
V. CONCLUSION
#_ H
q#
Dans cette étude, nous avons présenté deux topologies des
#_ HI -
chargeurs embarqués dans les VEH utilisant la technique de
q$
+
∆?
$_ H contrôle par hystérésis du courant absorbé. Cette technique
q$
Etats des
$_ HI Logique des Interrupteurs
- présente de nombreux avantages liés à la simplicité de son
q
+
∆?
commutations
_ H
implémentation, au temps de réponse réduit en régime
q
_ HI - dynamique, et à une précision et une stabilité satisfaisantes.
q#
+
∆? D’après les résultats obtenus, nous préconisons l’utilisation

Fig. 10 Principe de contrôle du courant par hystérésis des techniques de commande avancée, telle que la commande
pour un chargeur NPC à 3 niveaux
par hystérésis, pour les chargeurs des VEHR afin de respecter
Le principe de contrôle du courant illustré à la figure 10 les normes CEM basse fréquence spécifiques à ces dispositifs.
nous permet de créer les séquences des vecteurs de commande
suivantes [14]: [1100], [0110] et [0011]. REFERENCES
respectivement les variations temporelles de la tension () , et

[1] M. Marzouk “Développement de chargeurs intégrés pour véhicules
Dans les figures 11, 12 et 13, nous présentons hybrides plug-in”, Thèse Doctorat, université Grenoble, Octobre 2015.
[2] Journal officiel de l’Union Européenne, “Prescriptions uniformes
du courant absorbé par le chargeur ainsi que son spectre relatives à l’homologation des véhicules en ce qui concerne la
compatibilité électromagnétique”, 2012.
d’harmoniques. [3] A. Hanna Nohra “Commande de Filtres Actifs Parallèles sur un Réseau
Fortement Perturbé” Thèse Doctorat, Université Saint -Joseph,
Toulouse, Janvier 2017.
[4] Y. Murat, T. Philip “Review of Battery Charger Topologies, Charging
Power Levels, and Infrastructure for Plug-In Electric and Hybrid
Vehicles”, IEEE Trans. on power electronics, vol. 28, n°5, 2012.
[5] S. LACROIX “Modélisation et commande d’une chaîne de conversion
pour véhicule électrique intégrant la fonction de charge des batteries”,
Thèse Doctorat, Université Paris-Sud, Mai 2013.
[6] EMC limits for harmonic currents of equipment connected with input
Fig. 11 Tension du bus continu (Structure NPC 3 niveaux) current >16A and <75A per phase, IEC61000-3-12, 2004.
[7] N. MESBAHI “Contribution à l’Etude des Performances des
Onduleurs Multiniveaux sur les Réseaux de Distribution”, Thèse
Doctorat, Université Badji Bokhtar, Annaba, 2014.
[8] M. Messaoudi “Stratégie de modulation pour la réduction des
perturbations conduites générées par un variateur de vitesse à
redresseur MLI”, Thèse Doctorat, Ecole Doctorale Sciences pour
l’Ingénieur, Université Lille 1, Décembre 2015.
[9] S. Haghbin “An Isolated Integrated Charger for Electric or Plug-in
Fig. 12 Courant absorbé par le redresseur (Structure NPC 3 niveaux) Hybrid Vehicles”, University of Technology Goteborg, Sweden, 2011.
[10] S. Haghbin “Integrated Motor Drives and Battery Chargers for Electric
or Plug-in Hybrid Electric Vehicles”, Thesis for the degree of doctor of
philosophy, University of Technology Gothenburg, Sweden, 2013.
[11] A. Bouafia, “Techniques de commande prédictive et floue pour les
systèmes d’électronique de Puissance: application aux redresseurs a
MLI”, Thèse Doctorat, université Ferhat Abbas – Sétif, Algérie, 2010.
[12] N. Belhaouchet “Fonctionnement à Fréquence de Commutation
Constante des Convertisseurs de Puissance en Utilisant des Techniques
de Commande Avancées Application : Amélioration de la Qualité de
l’Energie”, Thèse Doctorat, Université Ferhat Abbas-Setif, 2011.
[13] Y. Touzani “Commande aléatoire du convertisseur Buck-Boost
triphasé à prélèvement sinusoïdal du courant”, Thèse Doctorat,
Fig. 13 Spectre d’harmoniques du courant absorbé par le chargeur
Université Blaise Pascal, Marseille, Juillet 2006.
(Structure NPC 3 niveaux)
[14] A. Dwivedi, A.N. Tiwari “Analysis of three-phase PWM rectifiers
D’après les résultats obtenus, nous constatons en premier using hysteresis current control techniques: a survey”, International
Journal of Power Electronics, Vol. 8, N°. 4, January 2017.
lieu que la tension à la sortie du chargeur atteint la tension de [15] A. Chaoui “Filtrage actif triphasé pour charges non linéaires”, Thèse
référence désirée (tension nécessaire pour charger la batterie Doctorat, Université Ferhat Abbas-Setif, Algérie, Octobre, 2010.
de VEH) pour les deux structures étudiées, ce qui montre la [16] Y. Loubna “Contribution à l'étude, la modélisation et la commande des
continu () . D’autre part, l’analyse des allures temporelles et

fiabilité de la boucle de régulation de la tension de bus redresseurs triphasés non-polluants: application au convertisseur trois-
niveaux à point neutre calé”, Thèse Doctorat, école de technologie
supérieure, université du Québec, 2004.
Copyright - 2019
ISSN 1737-9334 Page 19
Stratégie de contrôle du raccord au réseau du

générateur photovoltaïque
Habbati Bellia Assia#1, Moulay Benine Fatima*2
#
ENERGARID Laboratory, Tahri Mohamed University
Bechar, 08000, Algeria
bellia_abdeljalil@yahoo.fr
*
SIRECOM Laboratory, Djilali Liabes University
Sidi Bel Abbes 22000, Algeria
Moulay Fatima66@yahoo.fr
Abstract— La GD avec ses différents types ne cesse de s’intégrer La figure 1 est une représentation Simulink du système PV
dans le réseau électrique à travers le monde selon la connecté au réseau. Pour synchroniser la sortie de l’onduleur
particularité géographique et la disponibilité des ressources de avec le réseau, la tension et le courant au point de connexion
chaque pays. Les gestionnaires des réseaux électriques des pays de la charge sont prélevés afin de pouvoir contrôler
pionniers dans l’intégration de la GD ont établis des normes et
l’onduleur.
des recommandations pour raccorder la GD au réseau. Dans
notre article, nous allons nous intéresser particulièrement à Le sous-système de l’onduleur contient aussi un
l’énergie photovoltaïque (PV). transformateur élévateur de tension car nous n’avons utilisé
Malgré les avantages qu’apporte l’énergie photovoltaïque aux qu’un seul module PWX 500 (49W)
systèmes électriques et aux consommateurs, la fluctuation de la .
puissance de sortie des systèmes PV en raison de la variation
rapide de l'éclairement solaire cause quelques problèmes
opérationnels. [1], [2]. Cet article, présente la stratégie générale
de contrôle du raccord au réseau du générateur PV.
Keywords— photovoltaique, MPPT, Include at least 5 keywords

or phrases
I. INTRODUCTION
Le générateur PV est le dispositif responsable de convertir
la lumière en électricité. Cette électricité est sous forme
continue et variable. Un dispositif de régulation et
d’adaptation au réseau électrique est donc indispensable à Fig. 1 Modèle Simulink du système PV connecté au réseau
l’intégration de cette énergie dans le système électrique. En La boucle de régulation de la tension dont la structure est
général, l’électronique de puissance est très présente dans présentée dans la figure 2 est basée essentiellement sur la
cette interface. L’onduleur se sert de la tension du réseau pour mesure de la tension et du courant du réseau et sur
synchroniser les caractéristiques du générateurs PV. Un filtre l’utilisation de la boucle à verrouillage de phase PLL.
est utilisé pour remédier à la présence de perturbations dans
la tension du courant. La tension triphasée à la sortie du régulateur PI devient la
La littérature offre aujourd’hui une multitude tension de référence pour le générateur d’impulsions à MLI
d’architectures concernant l’assemblage : panneaux- (modulation de largeur d’impulsions) qui commande
convertisseurs DC/DC-onduleurs-filtres. Seules quelques- l’onduleur.
unes sont utilisées dans la pratique en raison du prix, de la
faisabilité technologique ou surtout de la sécurité.
L’intégration du convertisseur DC/DC suscite elle aussi un
intérêt non moindre. De nombreux chercheurs ont publié des
travaux de simulations et pratiques sur bancs d’essais et sur
sites concernant les différentes topologies. Dans cet travail,
On présente les éléments du système électro-énergétique
connectant le générateur PV au réseau électrique.
II. SYSTÈME ELECTROENRGETIQUE
Copyright - 2019
Page 20
ISSN 1737-9334
Mesures triphasées
Vabc et Iabc Les méthodes passives ont pour stratégie de détecter
l'amplitude de la tension ou le changement de fréquence qui
PLL se produisent lorsque le réseau est déconnecté, grâce au
circuit OVP/UVP ( Over Voltage protection/ Under Voltage
Com Protection) et le circuit l'OFP/UFP( Over Frequency
parat Abc/dq
protection/ Under Frequency Protection). L'équation (1)
eur
montre que si la
Vd,ref , V puissance réelle produite par le système PV et la puissance
réelle demandée par la charge ne coïncident pas lorsque
Générateur l'interrupteur est ouvert, la tension RMS au point de couplage
d’impulsions à commun devra, augmenter ou diminuer jusqu'à ce
MLI
Contrôle PI
que . les circuits OVP ou UVP sont alors activé.
U
De même, l’équation (2) indique que si la production de
Onduleur
puissance réactive du système PV et la demande réactive de
la charge ne coïncident pas lorsque le réseau est déconnecté,
la fréquence au point de couplage commun devra changer
Fig. 1 Structure de la boucle de régulation de la tension
jusqu'à ce que . les circuits OFP ou UFP son
alors activés.
Si le réseau est déconnecté par accident ou

intentionnellement, il est nécessaire de débrancher le
générateur PV. Afin de simuler la déconnexion du réseau, (1)
nous avons utilisé un sectionneur triphasé. Le générateur PV,
l’onduleur et la charge se trouvent ainsi en situation d’îlotage,
(2)
il faut alors, ouvrir l’organe de sectionnement qui est
représenté dans la figure 1 par un sectionneur monophasé La déconnexion du réseau a été simulé par l’ouverture du
commandé. Cette commande a été réalisée à l’aide de sectionneur à l’instant 0.05s. Les courants qui alimentaient la
l’application d’une méthode anti îlotage passive basée sur charge à partir du réseau s’annulent aussi à 0.05s. Les figures
l’analyse des paramètres électriques (tension et fréquence) au 4 et 5 montrent respectivement la disparition de la tension du
point de couplage commun. La figure 3 représente la structure réseau et l’annulation de ses courants
de ce sous-système.
Tension triphasée du réseau
Mesure de la tension au PCC 0.8
0.6
0.4
0.2
Détermination de
Tension (pu)
la fréquence 0
-0.2
-0.4
Limites de la Limites de la
tension fréquence -0.6
-0.8
0 0.02 0.04 0.06 0.08 0.1
Temps (s)
& Fig.4.Tension triphasée du réseau
Commande du sectionneur
Fig. 3 Structure de la méthode à seuil de tension

et à seuil de fréquence
Copyright - 2019
Page 21
ISSN 1737-9334
Courants du réseau Courant du générateur PV

8 4
3
6
2
4
Curant (A)
2
courant (A)
0
-1
-2 -2
-4 -3
0 0.02 0.04 0.06 0.08 0.1
temps (s)
-6
0 0.02 0.04 0.06 0.08 0.1
temps (s) Fig.8. Courant du générateur PV
Fig.6. Courant triphasé du réseau
courants-phase-une
15
La sortie de l’onduleur n’étant plus synchrone avec le 10

Icharge
Iréseau
réseau, elle commence à se synchroniser avec la charge qui Ipv
ne reçoit désormais le courant que du générateur PV via 5

l’onduleur. Cette situation va durer dans notre simulation plus
de 0.4s pour que le sectionneur isole le générateur PV.
Courants (A)
Pendant ce temps, nous remarquons une parfaite

0
synchronisation entre les courants de la charge et ceux

délivrés par le générateur PV dans les trois phases comme le -5
montre les figures 8, 9, 10.

L’objectif de cette simulation est d’arriver à isoler le -10
générateur PV quand il se trouve dans une situation d’îlotage

avec la charge. La méthode passive utilisée s’est avérée -15
0 0.02 0.04 0.06 0.08 0.1
efficace mais son temps de réponse est relativement long

Temps (s)
mais reste acceptable car les normes concernant les limites

d’opérations de la tension et de la fréquence sont très Fig.9Les courants dans la phase 1
différentes d’un pays à autre.
Courants-phase-deux
15
Tension PV après le transformateur Iréseau

0.8 Icharge
10
Ipv
0.6
5
0.4
Courants (A)
0.2 0
Tension (pu)
0
-5
-0.2
-10
-0.4
-0.6 -15
0 0.02 0.04 0.06 0.08 0.1
Temps (s)
-0.8
0 0.02 0.04 0.06 0.08 0.1
Temps (s) Fig.10. Les courants dans la phase 2
Fig.7. Tension PV après le transformateur

ISSN 1737-9334
Courants-phase-trois et celle du courant. Toute les deux ont suivi le changement de

15
l‘éclairement.
Iréseau
Icharge
Ipv
Nous avons associé par la suite un étage d’adaptation avec
10
MPPT pour stabiliser la tension délivrée par le générateur PV.
Une boucle de retour a été ajoutée pou réguler la tension en
commandant l’onduleur qui était à son tour associé à un
5
transformateur élévateur de tension. La méthode passive

Courants (A)
0
UOF /UOV a été utilisée pour déconnecter le générateur en
cas de perte du réseau.
-5
Conclusion
-10
Nous avons associé par la suite un étage d’adaptation avec

-15
0 0.02 0.04
Temps (s)
0.06 0.08 0.1 MPPT pour stabiliser la tension délivrée par le générateur PV.
Une boucle de retour a été ajoutée pou réguler la tension en
Fig.11. Les courants dans la phase 3 commandant l’onduleur qui était à son tour associé à un
transformateur élévateur de tension. La méthode passive
Afin de déterminer l’efficacité de ce modèle, il est utilisé UOF /UOV a été utilisée pour déconnecter le générateur en
pour alimenter une charge résistive. Nous avons changé cas de perte du réseau.
l’éclairement à 5s et nous avons tracé la courbe de la tension

ISSN 1737-9334
Analyse de l’écoulement du flux de puissance

dans un système éolien en présence
d’un creux de tension
Salima KADRI1, Mahmoud HABBAB1, Cherif BENOUDJAFER1
1
Département de génie électrique, Université Tarhi Mohammed, Béchar, Algérie
kaddrisalima@yahoo.fr
habbabmahmoud@gmail.com
benoudjafer@gmail.com
Abstract— L’objectif de ce travail est l’étude et la simulation La charge non linéaire est constituée d’un redresseur triphasée
d’un UPQC à commande classique connecté à un éolien. Etant en pont de diode alimentant une charge (R, L).
donner que l'éolienne participe en qualité d'une production
(source de tension) qui alimente avec le réseau les différents A. Modèle de la Turbine
charges connectées, il serait intéressant de voir le comportement La puissance cinétique du vent est exprimée par la
de l'ensemble du réseau avec un filtrage actif série –parallèle relation[1] :
UPQC pour une perturbation en tension (creux de tension) et 1
voir l’impact sur le comportement de l'ensemble du système.
𝑃𝑣 = . 𝜌. 𝑆. 𝑣 3 (1)
2
.
Keywords— Energie éolienne - UPQC – Qualité d’énergie - La puissance aérodynamique apparaissant au niveau du
Réseau rotor de la turbine s’écrit alors :
I. INTRODUCTION 1
𝑃𝑎𝑒𝑟 = 𝐶𝑝 . 𝑃𝑣 = . 𝐶𝑝 (𝜆, 𝛽). 𝜌. 𝑆. 𝑣 3 (2)
L'électricité est une forme d'énergie, qui doit être fournie 2

Où La masse volumique de l’air kg/m3, R la longueur
sous forme d'un système équilibré de tensions sinusoïdales
triphasées, à l'ensemble des divers utilisateurs. Par contre d’une pale en m, v la vitesse du vent m/sec et Cp est le
l’augmentation des problèmes de perturbations sur les réseaux coefficient de puissance qui caractérise le rendement
électriques, affecte de plus en plus cette qualité d'énergie. aérodynamique de la turbine. Il dépend des caractéristiques de
Parmi les perturbations les plus prépondérantes, les creux de la turbine (les dimensions de la pale, le ratio de la vitesse λ et
tension, les harmoniques, qui sont dus à la prolifération des l’angle d’orientation de la pale β). Le ratio de la vitesse est
charges non linéaires dans les équipements industriels et défini par :
Ω .𝑅
domestiques. 𝜆= 𝑡 (3)
𝑣
Diverses solutions de dépollution harmonique des réseaux
électriques ont été déjà proposées dans la littérature. Celles qui
Ou Ω𝑡 : La vitesse de rotation de la turbine,
répondent le mieux aux contraintes industrielles d’aujourd’hui Durant ce travail, ce coefficient sera modélisé par
sont le compensateur actif shunt, série et combinaison shunt - l’expression analytique suivante :[2]
𝑐5
série actifs (aussi appelé UPQC). L’UPQC présente la solution 1
𝐶𝑃 (𝜆, 𝛽) = 𝑐1 . (𝑐2 . − 𝑐3 . 𝛽 − 𝑐4 ) . 𝑒 − 𝐴 + 𝑐6 . 𝜆 (4)
universelle de compensation de tous les types de perturbations 𝐴
1 1 0.035
de courant et de tension sur un réseau électrique. Il associe les = − (5)
𝐴 𝜆+0.08.𝛽 1+𝛽 3
deux structures série et parallèle, le filtre série permet de Avec: c1= 0.5176, c2= 116, c3= 0.4, c4= 5, c5= 21 et c6=
protéger les charges sensibles et de compenser les 0.0068.
perturbations en tensions telles que les harmoniques, creux de La figure (2) représente la variation de ce coefficient en
tension ou les surtensions alors que le filtre parallèle compense fonction du ratio de vitesse λ pour différents angles
les perturbations en courants telles que les harmoniques et le d'orientation des pales β. On remarque que si, pour un angle
déséquilibres. constant 𝛽𝑖, on pouvait maintenir le coefficient de vitesse λ
constant et égal à 𝜆𝑜𝑝𝑡 à chaque instant, la puissance captée par
II. MODELE DE SYSTEME EOLIENNE l’éolienne serait maximale.
on obtient pour notre cas un coefficient de puissance
La structure du système de production d’énergie éolienne maximum de 0.48 pour un ratio de vitesse λ optimal qui vaut
est composée d’une aérogénérateur éolien qui est représenté 8.12 correspondant à un angle d’inclinaison des pales 𝛽 = 00 .
par un générateur asynchrone à cage. Le réseau d’alimentation
est un réseau triphasé représenté par une source programmable.

ISSN 1737-9334
0.5
A. Régulation du courant injecté par le FAP
beta=0 L´onduleur de tension connecté en parallèle au réseau
coefficent de puissance (Cp)
0.4 beta=5 électrique à travers un filtre passif. Le rôle du ce filtre inductif

beta=10
de sortie est de transformer l´onduleur de tension en une
beta=15
0.3
beta=20
source contrôlée de courant vis-à-vis du réseau.[6]
Le schéma de contrôle de la boucle de régulation du courant
0.2 par un régulateur PI peut être représenté comme suit dans la
figure (4).
0.1
0
0 2 4 6 8 10 12 14 16 18
ratio de vitesse spécifique (lamda)
Fig. 1 Coefficient de puissance Cp
Le couple aérodynamique apparaissant au niveau de la

turbine est donc en fonction de cette puissance:[3] Fig.3 : Schéma de régulation du courant injecté par le filtre actif parallèle
𝑃 𝜋
𝐶𝑎𝑒𝑟 = 𝑎𝑒𝑟 = 𝜌. 𝑅3 . 𝑣 2 . 𝐶𝑝 (𝜆, 𝛽) (6)
Ω𝑡 2.𝜆
Le multiplicateur adapte la vitesse de rotation de la turbine à B. Régulation de la tension injectée par le FAS
la vitesse de rotation de la génératrice, il est modélisé Le circuit de connexion entre le réseau électrique et le filtre
mathématiquement par les deux équations suivantes :[4] actif série est généralement à travers un filtre passif du
𝐶 deuxième ordre.[6]
𝐶𝑚 = 𝑎𝑒𝑟 (7)
𝐺 Le schéma global de la boucle de régulation de tension du
Ω𝑚é𝑐 = 𝐺. Ω𝑡 (8)
filtre actif série est donné par la Figure (5).
L’arbre du générateur est modélisé par l’équation
suivante:
𝐽𝑡 𝑑Ω
( + 𝐽𝑚 ) 𝑚𝑒𝑐 + 𝑓𝑣 . Ω𝑚𝑒𝑐 = 𝐶𝑚 − 𝐶𝑒𝑚 (9)
𝐺2 𝑑𝑡
Ainsi, nous pouvons établir le modèle de la turbine dont le
schéma bloc est donné par la figure (3).
Fig. 4 : Schéma de la régulation de la tension du F.A.S
C. Régulation de la tension du bus continu

La tension moyenne aux bornes du condensateur doit
être maintenue à une valeur constante. La cause de variation
Fig. 2 Modélisation de la partie mécanique de l’éolienne
de cette tension est l’échange de puissance active avec le
B. Modèle de la GAS réseau.[7]
Les équations des tensions statoriques et rotoriques de la En tenant compte du régulateur (PI), la boucle de régulation
GAS dans le repère de Park sont définies par :[5] de tension continue peut être représentée par la figure (6).
𝑑
𝑣𝑑𝑠 = 𝑅𝑠 𝑖𝑑𝑠 + ∅𝑑𝑠 − 𝜔𝑠 ∅𝑞𝑠
𝑑𝑡
𝑑
𝑣𝑞𝑠 = 𝑅𝑠 𝑖𝑞𝑠 + ∅𝑞𝑠 + 𝜔𝑠 ∅𝑑𝑠
𝑑𝑡
𝑑 (10)
𝑣𝑑𝑟 = 𝑅𝑟 𝑖𝑑𝑟 + ∅𝑑𝑟 − 𝜔𝑟 ∅𝑞𝑟 = 0
𝑑𝑡 Fig. 5: Boucle de régulation de la tension continue
𝑑
{ 𝑣𝑞𝑟 = 𝑅𝑟 𝑖𝑞𝑟 + 𝑑𝑡
∅𝑞𝑟 + 𝜔𝑟 ∅𝑑𝑟 = 0
III. COMMANDE DE L’UPQC IV. RESULTATS DE SIMULATION

L'objectif de la commande des onduleurs de l'UPQC est de En utilisant le logiciel MATLAB Simulink et le Toolbox
générer les instants d'ouverture et de fermeture des SimPower System afin de voir le comportement de la
interrupteurs qui les composent de façon à obtenir une tension génératrice vis-à-vis au perturbation engendrée dans un réseau
(ou courant) de sortie, côté alternatif, qui suit le plus de distribution, dont les paramètres sont donnés dans le tableau
fidèlement possible la consigne (tension ou courant de 1
référence) en minimisant ses harmoniques non désirés
(harmoniques de commutation) et ce, en optimisant l'amplitude
pour une tension du bus continue donnée.

ISSN 1737-9334
TABLE I
500
PARAMETRE DU SYSTEME EOLIONNE
Vs[V]
Source triphasé (a) 0
Vs, f 220V, 50Hz

-500
Rs, Ls 0.1Ω, 0.2 mH 0 0.1 0.2 0.3 0.4 0.5 0.6
t[s]
500
UPQC
VL[V]
Vdc 735.8V (b) 0
Cdc 17.4mF -500

0 0.1 0.2 0.3 0.4 0.5 0.6
Charge non linéaire t[s]
P_active 8000W 100
Vinj[V]
P_reactive 500Var (c) 0
GAS -100
0 0.1 0.2 0.3 0.4 0.5 0.6
Sn 100KVA t[s]
Vn 380V
f 50HZ Fig. 6 Tension de: a) Source, b) charge, c) injecté
RS 0.016p.u
LS 0.06p.u
100
Rr 0.015p.u
Is[A]
(a) 0
Lr 0.06p.u -100
0 0.1 0.2 0.3 0.4 0.5 0.6
Lm 3.5p.u t[s]
50
P 2
IL[A]
(b) 0
-50
0 0.1 0.2 0.3 0.4 0.5 0.6
Les Figure (6-7-9-11-12) illustre le comportement de t[s]
50
l’UPQC lors de la compensation d’un creux de tension triphasé
Iinj[V]
(c) 0
d’une source sinusoïdale équilibrée non polluée. Le filtre actif -50
0 0.1 0.2 0.3 0.4 0.5 0.6
parallèle assure parfaitement sa fonction et compense bien les t[s]
Iéolien[A]
100
courants harmoniques, Il en résulte un courant source (d) 0
sinusoïdal de THD est égal à 3,44% qui est en phase avec la -100
0 0.1 0.2 0.3 0.4 0.5 0.6
tension source. t[s]
Entre les instants t1=0.25 sec et t2=0.35 sec, un creux de

tension d’une profondeur de 30% est introduit volontairement. Fig. 7 Courant de: a) Source, b) charge, c) injecté, d) éolien
Le filtre série réagit en fournissant la puissance demandée à la
charge et cela par l’injection des tensions de compensation qui 5
Amplitude (% de Fondamentale)
sont égales à la différence entre la tension de référence de la

4
charge et la tension du réseau. Il en résulte une tension
sinusoïdale est égale à la tension désirée aux bornes de la 3
charge même pendant le creux de tension. L’augmentation du
2
courant de source est due à la recharge du bus continu lors de
la compensation du creux de tension. 1
La tension du bus continu Udc suit parfaitement sa
0
référence, d’une façon très rapide avec une légère chute de 0 500 1000 1500 2000 2500
Fréquence [Hz]
tension lors de la compensation du creux. Ismax= 38.11A , THD= 3.44%
Le comportement de la génératrice vis-à-vis un creux de
tension est illustré par l’allure des figures (7.d, 10, 11.b,
12.b).en effet lors d’un creux de tension, la turbine va
Amplitude [% de Fondamentale]
15
accélérer, en plus une démagnétisation de la génératrice. On
remarque aussi à la fin du défaut, lorsque la tension redevient
10
normale, la génératrice peut se re-magnétiser qui peut
provoquer un appel de courant important.
Le comportement de la puissance active et réactive du 5
système est illustré dans la figure (9-10). On note que
l'éolienne participe en qualité d'une production (source de 0
tension) qui alimente avec le réseau la charge connectée. 0 500 1000 1500 2000 2500
Fréquence[Hz]
ILmax = 27.65A , THD= 19.63%
Fig. 8 Spectre des harmoniques du courant source et de charge

ISSN 1737-9334
740 10000
735
Q-source[Var]
5000
Vdc[V]
730
0
725
720 -5000
0 0.1 0.2 0.3 0.4 0.5 0.6 0 0.1 0.2 0.3 0.4 0.5 0.6
t[s] t[s]
Fig. 9 Tension de bus continue Vdc (a)
8000
1.025
6000
1.02
Q-éolien[Var]
4000
1.015 2000
V[p.u]
1.01 0
1.005 -2000
-4000
1 0 0.1 0.2 0.3 0.4 0.5 0.6
0 0.1 0.2 0.3 0.4 0.5 0.6 t[s]
t[s] (b)
Fig. 10 Vitesse mécanique 150
8000
Q-charge[Var]
6000 100
P-source[W]
4000
50
2000
0
0
-2000
0 0.1 0.2 0.3 0.4 0.5 0.6
0 0.1 0.2 0.3 0.4 0.5 0.6 t[s]
t[s] (c)
(a)
2000
Fig. 12 Puissance réactive de : a) Source, b) Eolien, c) Charge
0
V. CONCLUSIONS
P-éolien[W]
l’UPQC permet effectivement de compenser à la fois la

-2000
perturbation harmonique produite par la charge non-linéaire
afin d’avoir un courant source moins polluant et de fournir une
-4000
tension à la charge quasi sinusoïdale de malgré la perturbation
coté source. Ceci est dans le but de protéger et rendre la
-6000
0 0.1 0.2 0.3 0.4 0.5 0.6 tension de charge parfaitement sinusoïdale.
t[s]
(b) ACKNOWLEDGMENT
6000
Ce travail a été effectué au Laboratoire Smart Grids
5000
Energies renouvelables et l’Université Tahri Mohamed Bechar
4000 08000 ALGERIE, dans le cadre du projet scientifique:
P-charge[W]
3000 «Energies renouvelables»; 2015/2018.

2000 REFERENCES
1000 [1] S. El Aimani, “ Modélisation de différentes technologies d’éoliennes
intégrées dans un réseau de moyenne tension”, Thèse de doctorat de
0 l’école centrale de Lille, 2004.
0 0.1 0.2 0.3 0.4 0.5 0.6
[2] N. Laverdure, “Sur l'intégration des générateurs éoliens dans les
t[s]
réseaux faibles ou insulaires”, Thèses de Doctorat, ENS de Cachan,
(c)
2005.
Fig. 11 Puissance active de : a) Source, b) Eolien, c) Charge

ISSN 1737-9334
[3] R. Rouabhi, “Contrôle des puissances générées par un système éolien

à vitesse variable basé sur une machine asynchrone double alimentée ”,
Thèse de Doctorat, l’Université de Batna, 2016.
[4] A. Gaillard, “ Système éolien basé sur une MADA : contribution à
l’étude de la qualité de l’énergie électrique et de la continuité de
service ”, Thèse de Doctorat, l'Université Henri Poincaré, Nancy-I,
2010.
[5] S. Bachir “Contribution Au Diagnostic De La Machine Asynchrone Par
Estimation Paramétrique ”, Thèse de Doctorat, Université de Poitiers
2002.
[6] G.S. Eskander, « Apport de l’UPQC à l’amélioration de la stabilité
transitoire des réseaux électriques », Université Henri Poincaré, thèse
de doctorat, Nancy-I. 2003.
[7] R. Tounsi, «Développement d'un contrôle commande pour un
compensateur série de creux de tension. Validation par simulation du
fonctionnement avec des charges industrielles », Thèse de doctorat de
l’INPT 11 Janvier 1999.
[8]
[9] M. Boulouiha, “Contribution à l’étude et la commande d’une éolienne
pour la production de l’énergie électrique ”, Thèse de Doctorat, 2013

ISSN 1737-9334
3ème Congrès International sur les Energies Renouvelables et le Développement

Durable
30-31 juillet 2019, Tunisie-Monastir
ÉCOULEMENT HYDROMAGNETIQUE NON

LINEAIRE DE CONVECTION FORCEE D’UN FLUIDE
BIOMAGNETIQUE INSTABLE
MERABET AYECHE Chahra*,1, KEZZAR Mohamed1, SARI Mohamed Rafik 2
1
Département de Mécanique, Université 20 Août 1955 de Skikda BP 26, 21000 Algérie.
2
Laboratoire de Mécanique Industrielle, Université Badji Mokhtar d’Annaba BP 12, 23000 Algérie.
* Email : chahram21@yahoo.fr
𝑆 Coefficient de la viscosité de vortex

𝑡 Temps
Résumé— L'objectif de ce travail est d’analyser les 𝑇 Température de la couche limite
caractéristiques du champ dynamique d'un 𝑇𝑓 Température intérieure
écoulement bidimensionnel instable hydromagnétique 𝑇𝑤 Température sur la surface
de la couche limite d'un fluide biomagnétique au- 𝑇∞ Température du fluide ambiant
dessus d'une cale, utilise le modèle du micro polaire 𝑢 Vitesse suivant x
liquide avec une surface convective tenue compte 𝑈 Vitesse du jet libre
l'effet induit de champ magnétique. Le fluide 𝑈0 Vitesse caractéristique
𝑈∗ Vitesse du jet libre sans dimension
biomagnétique le plus connu qui présente de telles
𝑣 Vitesse suivant y
propriétés magnétique est le sang qui a un nombre
𝑋 Longueur caractéristique
de Prandtl=21. En présentant des variables non
𝑥, 𝑦 Coordonnée de la surface
dimensionnelles appropriées le problème est décrit Symboles grecs
par un système couplé et non linéaire des équations 𝜌 Masse volumique
ordinaires avec ses conditions aux limites qui est 𝛽 Angle de cale
résolu numériquement en appliquant la méthode de 𝛿 Balance de longueur
Range-Kita d’ordre 4 (RK4) associée par la méthode 𝜇 Viscosité dynamique
de Tir. 𝜇𝑒 Perméabilité magnétique
𝜐 Viscosité cinématique
Les résultats obtenus sont présentés en forme des Δ Paramètre de viscosité de vortex
graphiques de vitesse adimensionnelle comprenant 𝜉 Paramètre de micro-inertie
l'épaisseur de la couche limite pour différentes valeurs 𝜔 Micro-rotation
des paramètres entrant dans le problème. 𝜎 Conductivité électrique
𝜓 Fonction de courant
Mots clé — Fluide biomagnétique, couche limite, 𝜅 Conductivité thermique.
champs magnétique, surface convective, RK4, la 𝜂 Paramètre de similarité.
méthode de Tir.
NOMENCLATURE 1. INTRODUCTION
𝐶𝑓 Coefficient de frottement superficiel La dynamique des fluides Biomagnétique (BFD) est

𝑐𝑝 Chaleur spécifique à pression constant un domaine relativement nouveau dans la
𝑓 Fonction de courant sans dimensions dynamique des fluides et pendant les dernières
ℎ Paramètre induit de champ magnétique décennies des travaux de recherches étendues ont
ℎ𝑓 Coefficient du transfert de chaleur
été effectué sur la dynamique des fluides
𝐻∗ Intensité induit du champ magnétique
𝐻 𝐻 ∗ sans dimension
biologiques en présence du champ magnétique.
𝐻0 Intensité du champ magnétique appliqué L'influence du champ magnétique sur l'écoulement
𝐻𝑤 Champs magnétique induit des biofluids est intensivement étudiée en raison de
𝑗 Micro-inertie par unité de masse diverses applications dans la technologie biologique
𝐾 Paramètre d’instabilité et la médecine.
𝑀 Paramètre du champ magnétique
𝑚 Exposant de vitesse Le biofluid est existe dans la créature vivante dont
𝑁 Paramètre de micro-rotation sans l'écoulement est influencé par la présence d'un
dimensions
𝑛 Paramètre de micro-rotation champ magnétique nommé un fluide
𝑃 Pression biomagnétique. Le biofluide le plus célèbre est le

ISSN 1737-9334

Durable
sang qui tient les propriétés d'un fluide magnétique d'hydrodynamique ferro (FHD), qui s'occupe sans
en raison de l'existence des oxydes de fer dans les courant électrique induit et considère que
molécules d’hémoglobine, qui est présente à une l'écoulement est affecté par la magnétisation du
concentration uniquement élevée dans les globules fluide dans le champ magnétique (Haik et al 1996,
rouges mûres (Pauling et al 1936, Higashi et al Haik et al 1999 et Haik et al 2002).
1993). Le sang possède la propriété du matériel
magnétique et dans quelques circonstances, peut L'objectif du travail actuel est l'étude numérique du
être considéré comme fluide diamagnétique, transfert de chaleur des écoulements
paramagnétique ou ferromagnétique (Haik et al hydromagnétique bidimensionnel instable de la
1996). couche limite des fluides biomagnétique au-dessus
d'une cale. On utilise le modèle du fluide
La première formulation de la dynamique des micropolaire avec une surface convective tenue
fluides biomagnétiques (BFD), pour la recherche compte l'effet induit de champ magnétique. Les
sur l'écoulement d'un biofluid sous l'influence d'un équations de gouvernement couplées non linéaires
champ magnétique appliqué (écoulement de fluide sont résolues numériquement pour les solutions
biomagnétique), a été développée par (Haik et al semblables locales appliquant la méthode de
1996). Selon cette formulation, le fluide Range-Kita d’ordre 4 (RK4) associée par la
biomagnétique est considéré comme réellement méthode de Tir. Des résultats graphiques pour la
posséder les propriétés magnétiques du sang. Le vitesse adimensionnelle sont présentés pour une
sang dans cette formulation de BFD est considéré gamme des valeurs des paramètres caractéristique.
comme fluide magnétique homogène, newtonien, La discussion d’accompagnement fournit des
électriquement non conducteur. Clairement, le interprétations physiques des résultats.
modèle mathématique de (Haik et al 1996), est
valide pour le flux de sang laminaire dans de grands
v
navires où le comportement newtonien est une u
bonne approximation (Zendehbudi et al 1999, g
Anderson et al 2000]. D'ailleurs, le modèle de y
(BFD) est basé réellement sur celui
x 2
2. FORMULATION DU PROBLEME 2

Horizontal line
Equations de gouvernement :
Nous considérons un écoulement laminaire

bidimensionnel instable d'un fluide biomagnétique
Fig. 1 Configuration d'écoulement
visqueux de masse volumique constante ρ, de
viscosité dynamique µ et de température 𝑇∞ se les forces de corps et l'effet de Joule, peut être
déplaçant au-dessus d'une surface électrique non décrit par les équations de conservation suivantes
conductrice d'une cale, dirigée le long de l'axe des (Rahman et Sattar 2006, Chaudhary et Sharma
abscisses positif, avec une vitesse potentielle U (x , 2006, Rahman 2011).
t) conduit par le gradient de pression du fluide non
visqueux, en présence d'un champ magnétique Equation de continuité
transversal uniforme appliqué, comme illustré dans 𝜕𝑢 𝜕𝑣
+ =0 (1)
la fig.1. L’angle total de la cale est donné par 𝛽𝜋. 𝜕𝑥 𝜕𝑦
On suppose que l'intrados de la cale est chauffé par

Equation de quantité de mouvement
la convection d'un fluide chaud à la température
𝑇𝑓 qui fournit un coefficient de transfert de 𝜕𝑢
+𝑢
𝜕𝑢
+𝑣
𝜕𝑢
=−
1 𝜕𝑃
+ (𝜐 + )
𝑆 𝜕2 𝑢
+
𝑆 𝜕𝜔
+
𝜕𝑡 𝜕𝑥 𝜕𝑦 𝜌 𝜕𝑥 𝜌 𝜕𝑦 2 𝜌 𝜕𝑦
chaleur ℎ𝑓 . Un flux magnétique uniforme 𝐻0 est 𝜇𝑒 𝐻0 𝜕𝐻 ∗
(2)
appliqué en parallèle à l'axe des ordonnées de sorte 𝜌 𝜕𝑦
qu'il sature le fluide biomagnétique.

Équation de moment
Dans le cadre des prétentions au-dessus- 𝜕𝜔 𝜕𝜔 𝜕𝜔 𝜐𝑠 𝜕2 𝜔 𝑆 𝜕𝑢
remarquables, l'écoulement convectif d'un fluide +𝑢 +𝑣 = − (2𝜔 + ) (3)
𝜕𝑡 𝜕𝑥 𝜕𝑦 𝜌𝑗 𝜕𝑦 2 𝜌𝑗 𝜕𝑦
biomagnétique incompressible instable négligeant
Copyright - 2019
ISSN 1737-9334 Page 30

Durable
Équation d'induction magnétique 𝐵𝑖 = √(ℎ𝑓 ⁄𝜅)(𝑥𝜐⁄𝑈0 𝑈∗ )
𝜕𝐻 ∗ 𝜕𝐻 ∗ 𝜕𝐻 ∗ 1 𝜕2 𝐻 ∗ ∗ 𝜕𝑢 𝜕𝑢
+𝑢 +𝑣 = +𝐻 + 𝐻0
𝜕𝑡 𝜕𝑥 𝜕𝑦 𝜇𝑒 𝜎 𝜕𝑦 2 𝜕𝑥 𝜕𝑦  Paramètre induit de champ magnétique
(4)
ℎ = √(𝐻𝑤 ⁄𝑈0 𝑈∗ )(𝜇𝑒 ⁄𝜌)
Conditions aux limites
 Nombre de Reynolds local
Sur la surface de la cale à (𝑦 = 0):
𝜕𝑇 𝜕𝑢 𝑅𝑒 = 𝑈0 𝑥⁄𝜐
𝑢 = 0, 𝑣 = 0, −𝜅 = ℎ𝑓 (𝑇𝑓 − 𝑇), 𝜔 = −𝑛 , 𝐻∗ =
𝜕𝑦 𝜕𝑦
𝐻𝑤 . (5.a)  Nombre de Prandtl magnétique
Assortiment avec le jet libre à (𝑦 → ∞): 𝑃𝑚 = 𝜇𝑒 𝜎 𝜐
𝑢 = 𝑈(𝑥, 𝑡), 𝜔 = 0, 𝑇 = 𝑇∞ , 𝐻 ∗ = 0 (5.b)  Paramètre de viscosité de vortex
Pour obtenir les solutions des équations Δ = 𝑆 ⁄𝜇

gouvernement ; 1, 2, 3 et 4 ainsi que les conditions  Paramètre local de champ magnétique
aux limites Eqs 5.a et 5.b, nous présentons les
transformations de similitude suivante : 𝑀 = √(𝜇𝑒 𝑥⁄𝜇𝑈0 𝑈∗ )
𝑚+1 𝑥 𝑚−1 2𝜐2 𝑥 𝑚+1  Paramètre de micro-inertie

𝜂 = 𝑦√ ,𝜓 = √ 𝑓(𝜂) (6)
2 𝛿 𝑚+1 (𝑚+1)𝛿 𝑚+1
𝜉 = 𝑗𝑈0 𝑈∗ /𝜐𝑥 Avec : 𝑈∗ = 𝑈⁄𝑈0
𝑚+1 𝜐2 𝑥 3𝑚−1 (𝑚+1)𝜌 𝜐𝑥 𝑚
𝜔=√ 𝑁(𝜂), 𝐻 ∗ = √ 𝐻(𝜂) Le système d’équations ; 9, 10 et 11 est semblable
2 𝛿 3𝑚+3 2𝜇𝑒 𝛿 𝑚+1
𝑈0𝑚−1 𝑑𝛿
(7) au terme [ 𝛿𝑚 ] où t apparaît explicitement.
𝜐𝑚 𝑑𝑡
Où 𝜓 est la fonction de courant et f, N et H sont Ainsi la condition de similitude exige que ce doive
respectivement ; la fonction de courant être une constante. Par conséquent, suivant le
adimensionnelle, la microrotation, et le champ travail de Sattar (2011) et des références là-dedans,
magnétique induit. on peut essayer une classe des solutions d'Eqs. 9, 10
et 11 en assumant cela
Depuis : 𝑢 = 𝜓𝑦 et 𝑣 = −𝜓𝑥 , On trouve que :
[ 𝛿𝑚 ] = 𝐾(𝑐𝑜𝑛𝑠𝑡𝑎𝑛) (12)
𝜐𝑚 𝑑𝑡
′ 𝑚+1 𝜐2 𝑥 𝑚−1 𝑚−1
𝑢 = 𝑈(𝑥, 𝑡)𝑓 et 𝑣 = −√ (𝑓 + 𝜂𝑓 ′ )
2 𝛿 𝑚+1 𝑚+1
L’intégrale de l’équation 12, on donne :
(8)
1
Maintenant la substitution des équations ; 6, 7, et 8 𝛿 = [𝐾(𝑚 + 1)𝑈01−𝑚 𝜐 𝑚 𝑡]𝑚+1

dans le système des équations ; 2,3 et 4, nous Après les arguments de Sattar (2011) et d'employer
obtenons les équations différentielles non linéaires Eqs. 9, 10, 11 et 12 devient des équations ordinaires
suivantes : comme suit :
[ 𝛿 𝑚 ] [𝜂𝑓 ′′ + 2𝑓 ′ − 2] + 𝑅𝑒 𝑚−1 [(1 + Δ)𝑓 ′′′ + 𝐾[𝜂𝑓 ′′ + 2𝑓 ′ − 2] + 𝑅𝑒 𝑚−1 [(1 + Δ)𝑓 ′′′ + 𝑓𝑓 ′′ +
𝜐𝑚 𝑑𝑡
𝑓𝑓 ′′ + 𝛽(1 − 𝑓 ′2 ) + Δ𝑁 ′ + 𝑀𝐻 ′ ] = 𝛽(1 − 𝑓 ′2 ) + Δ𝑁 ′ + 𝑀𝐻 ′ ] = 0 (13)
0 (9)
𝐾[3𝑁 + 𝜂𝑁 ′ ] + 𝑅𝑒 𝑚−1 [(1 + Δ⁄2)𝜉𝑁 ′′ − (2 −
𝑈0𝑚−1 𝑑𝛿 𝛽)Δ(2𝑁 + 𝑓 ′′ ) + 𝜉{𝑓𝑁 ′ + (1 − 2𝛽)𝑁𝑓 ′ }] = 0
[ 𝛿 𝑚 ] [3𝑁 + 𝜂𝑁 ′ ] + 𝑅𝑒 𝑚−1 [(1 + Δ⁄2)𝜉𝑁 ′′ − (14)
𝜐𝑚 𝑑𝑡
(2 − 𝛽)Δ(2𝑁 + 𝑓 ′′ ) + 𝜉{𝑓𝑁 ′ + (1 − 2𝛽)𝑁𝑓 ′ }] = 𝐾[2𝐻 + 𝜂𝐻 ′ ] + 𝑅𝑒 𝑚−1 [𝐻 ′′ + 𝑃𝑚𝑓𝐻 ′ + (𝛽 −
0 (10)
1)𝑃𝑚 𝜂𝐻𝑓 ′′ + (2 − 𝛽)𝑃𝑚𝑀𝑓 ′′ ] = 0 (15)
[ 𝑚 𝛿 𝑚 ] [2𝐻 + 𝜂𝐻 ′ ] + 𝑅𝑒 𝑚−1 [𝐻 ′′ + 𝑃𝑚𝑓𝐻 ′ + Ainsi en employant Eqs. 7 et 8, les conditions aux
𝜐 𝑑𝑡
(𝛽 − 1)𝑃𝑚 𝜂𝐻𝑓 ′′ + (2 − 𝛽)𝑃𝑚𝑀𝑓 ′′ ] = 0 (11) limites Eqs. (6.a), et (6.b) deviennent
respectivement
Avec :
𝑓 = 𝑓 ′ = 0, 𝑁 = −𝑛𝑓 ′′ , 𝐻 = √2 − 𝛽 ℎ à 𝜂 = 0 (16.a)
 Nombre de Biot
Copyright - 2019
ISSN 1737-9334 Page 31

Durable
𝑓 ′ = 1, 𝑁 = 0, 𝐻 = 0 à 𝜂 → ∞ (16.b)
3. RESULTATS
L'ensemble des équations ; 13, 14 et 15 sont

fortement non linéaires et couplé et donc le système
ne peut pas être résolu analytiquement. Alors que
cette dernières et les conditions aux limites
correspondants aux équations; 16.a et 16.b, sont
résolus numériquement utilise la méthode de
Range-Kita d’ordre 4 (RK4) associée par la
méthode de Tir.
Fig. 4 Profil de vitesse pour différents valeurs de Re et
Les figures 2, 3, 4, 5, 6 et 7, nous avons montré la β>1.
vitesse adimensionnelle du fluide biomagnétique
dans la couche limite pour différents paramètres. La Les figures 3 et 4 nous avons présenté la vitesse
figure 2 montre l'effet du paramètre d'instabilité k adimensionnelle du fluide biomagnétique pour
sur la vitesse adimensionnelle du fluide différentes valeurs du nombre de Reynolds 0≤β<1
biomagnétique. Pour cette figure (2) nous notons et 1< β<2. Pour la figure 3, on note que l'effet du
que la vitesse du fluide biomagnétique diminue nombre de Reynolds est très semblable au
dans toute la couche limite avec l'augmentation du paramètre d'instabilité K quand β<1. Dans les
paramètre d'instabilité. équations ; 13, 14 et 15 nous notons que Re est
apparaît comme puissance de m-1. Alors que pour
0≤m<1 et 0≤β<1, le nombre de Reynolds est agit
comme le paramètre d'instabilité. De fig.4 nous
constatons que la vitesse du fluide biomagnétique
augmente avec l'augmentation du nombre de
Reynolds dans un certain domaine 𝜂 < 𝜂 critique.
Mais lointain de la surface du corps (𝜂 > 𝜂 critique)
que cette situation décompose et la conséquence est
que le profil de vitesse diminue. Les effets d'un plus
grand nombre de Reynolds sur les profils de vitesse
diminuent.
Fig. 2 Profil de vitesse pour différents valeurs de K.
Fig. 5 Profil de vitesse pour différents valeurs de β.
L'effet des changements du paramètre d'angle de

cale 𝛽 sur la vitesse sans dimensions dans la
couche limite est montré dans la fig.5. La
caractéristique générale de ces solutions est
Fig. 3 Profil de vitesse pour différents valeurs de Re et l'augmentation de la vitesse du liquide
β<1. biomagnétique avec l'augmentation du 𝛽 .
Les résultats prouvent également que les profils de

vitesse sont devenus plus raides pour les valeurs

ISSN 1737-9334

Durable
plus grandes du paramètre d'angle de cale 𝛽. Ce pour déterminer la vitesse du liquide. Les
dernier est une mesure du gradient de pression, et influences du paramètre d'instabilité, du nombre de
ainsi une valeur positive de 𝛽 indique un gradient Reynolds, du paramètre de champ magnétique, du
de pression négative. paramètre d'angle de cale, et du paramètre induit de
champ magnétique sur le champ d'écoulement sont
examinées. Les résultats sont présentés sous forme
de graphiques. De la présente recherche les
conclusions principales suivantes peuvent être
tirées :
• L’instabilité est commande de manière

significative les caractéristiques
d'écoulement et de transfert de chaleur du
fluide biomagnétique.
• la vitesse d'écoulement diminue avec
l'augmentation du paramètre d'instabilité
tandis qu'elle augmente avec
l'augmentation de l'angle de la cale.
Fig. 6 Profil de vitesse pour différents valeurs de M.
• Le champ magnétique induit a un effet
significatif sur l'écoulement de sang.
Ainsi, dans la modélisation du fluide
biomagnétique, l'effet du champ
magnétique induit dans le domaine
d'écoulement devrait être tenu compte.
REFERENCES
[1] Anderson, H. I., Halden, R. et Glomsaker, T. (2000).
Journal of Biomechanics. 33, 1257.
[2] Chaudhary, R.C. et Sharma, B.K. (2006). “Combined heat

and mass transfer by laminar mixed convection flow from a
Fig. 7 Profil de vitesse pour différents valeurs de h. vertical surface with induced magnetic field,”. J Appl Phys,
99:034901–034910.
Pour quelques valeurs représentatives du paramètre
[3] Haik, Y., Chen, J.C. and Pai, V.M. (1996).
M de champ magnétique. La fig. 6 nous voyons que
“Development of bio-magnetic fluid dynamics,”
les profils de vitesse diminuent avec l'augmentation Proceedings of the IX International Symposium on
du paramètre de champ magnétique M, qui indique Transport Properties in Thermal Fluids
Engineering, Singapore, Pacific Center of
que le champ magnétique agit comme une corde qui Thermal Fluid Engineering, S.H. Winoto, Y.T.
retarde le mouvement du liquide. Chew, N.E. Wijeysundera, eds, Hawaii, U.S.A.,
June 25–28, , 121–126.
Dans la fig. 7, l'effet du paramètre induit du champ
[4] Haik, Y., Chen, C. J. et Chatterjee, J. (2002). Journal of
magnétique h sur le champ de vitesse est dépeint Visualization, 5(2), 187.
pour les solutions instables. On le constate que la
vitesse du fluide biomagnétique diminue de [5] Haik, Y., Pai, V. et Chen, C. J. (1999). Biomagnetic Fluid
Dynamics, In: Fluid Dynamics at Interfaces, W. Shyy and R.
manière significative avec l'augmentation du Narayanan (Eds), Cambridge University Press.
paramètre induit du champ magnétique.
[6] Higashi, .T., Yamagishi, A., Takeuchi, T., Kawaguchi, N.,
Sagawa, S., Onishi, S. et Date, M. (1993). Journal of Blood, 82
(4), 1328.
4. CONCLUSION [7] Pauling, L. et Coryell, C. D. (1936). “The magnetic
Properties and Structure of Hemoglobin, Oxyhemoglobin and
L'écoulement hydromagnétique non linéaire de Carbonmonoxy Hemoglobin,” In Proceedings of the National
convection forcée du fluide biomagnétique instable Academy of Science, USA, 22, 210-216.
au-dessus d'une cale avec la condition de la surface [8] Rahman, M.M. (2011). “Heat transfer in biomagnetic fluid
convective en présence du champ magnétique over a wedge with convective surface boundary condition in
the presence of induced magnetic field,” Int J Energy Technol
induit a été analysé théoriquement. L'analyse de ,3(17):1–8.
similitude et la méthode numérique ont été utilisées
Copyright - 2019
Page 33
ISSN 1737-9334

Durable
[9] Rahman, M.M. et Sattar, M.A. (2006).
“Magnetohydrodynamic convective flow of a micropolar fluid
past a continuously moving vertical porous plate in the presence
of heat generation/absorption,” ASME J Heat Tran ,128(2):142–
152.
[10] Sattar, M.A. (2011). “A local similarity transformation for

the unsteady two-dimensional hydrody- namic boundary layer
equations of a flow past a wedge,” Int J Appl Math Mech,
7(1):15–28.
[11] Zendehbudi, G. R. et Moayeri, M. S. (1999). Journal of

Biomechanics. 32, 959.

ISSN 1737-9334
Le bâtiment passif pour assurer l'efficience

énergétique
Mokhnache nardjesse#1, Alkama djamel*2
#
departement d'architecture, loboratoire LGCH université 08 mai 1945
Algerie, Guelma
mokhnache.nardjess@gmail.com
dj.alkama@gmail.com
Resumé— Avec la demande excessive en énergie fossile une donc de modifier le climat a fin de répondre aux besoins des
quantité colossal de gaz a effet de serre est émise dans habitants.
l'atmosphère provoquant d'importants écart de température qui Dans le but d'assurer le confort thermique l'habitant a
engendre de grands besoins de régulation de l'environnement recours a des comportement énergivore pour réguler son
thermique par le billet des énergies fossiles, Le secteur résidentiel
environnement thermique, en puisant dans les énergies fossile
est souvent qualifié d'énergivore suite au mode de vie adopté a fin
d'assurer le confort thermique. D'ou le besoin de penser des a sa disposition, mais toute consommation énergétiques
solutions pour assurer le confort thermique et l'efficience engendre une facture écologique qui se traduit par
énergétique des habitations . Le but de cette recherche est de l'accoisement de l'empreinte écologique, c'est une estimation
démontrer que le bâtiment passif garanti une basse de la surface terrestre nécessaire a l'homme pour subvenir à ses
consommation énergétique tout au long du cycle de vie de besoins : cet outil est une mesure de la pression qu'exerce
l'habitation, Il s'agit aussi de souligner que l'énergie renouvelable l'Homme sur la nature. D'après le Global Footprint Network,
est l'alternative idéal au énergies fossile communément utilisées notre mode de vie nécessite actuellement 1,75 planète pour
par les habitants a fin de diminuer l'empreinte écologique que notre consommation des ressources soit compensée au
engendrée par le mode de vie adopté, il est aussi question de
niveau global. Chaque année, ce réseau calcule le "jour du
mettre l'accent sur le rôle que joue l'efficience énergétique dans le
remboursement de la dette écologique, sur un échantillon dépassement", le point à partir duquel l'humanité a consommé
d'habitation situé dans le climat semi-aride de la ville enclavé de les ressources que la Terre peut compenser en une année. En
Guelma bénéficiant d'un grand potentiel bioclimatique. La 2019, le jour du dépassement était le 29 juillet. Un habitant
corrélation entre l'étude de Givoni ,ainsi qu'une approche par moyen de la Terre n'aurait besoin que de cinq terrains de foot
mesure in situ a l'intérieur de l'habitation : température, pour y vivre, Le mode de vie adopté par l'être humain cause
humidité relative, vitesse de l'air, révèlent des température des répercutions irréversibles sur l'environnent a cause des gaz
intérieur au-delà des plages de confort et les habitant souffre a effet de serre qui cause la hausse de température du globe
situation stress due a l'inadéquation des matériaux de terrestre ce qui accroit la situation d'inconfort.
construction et spécialement le soucis de l'inertie thermique qui
peut être la solution pour améliorer le confort thermique a II. LE CONFORT THERMIQUE
l'intérieur de l'habitation. Les résultats de cette recherche
tendent a démontré qu'il faut produire un bâtiments passif en Le confort thermique est une notion a vaste définition qui
exploitant potentiel bioclimatique du site pour améliorer a la fois dépond de facteurs objectifs, et subjectifs, c'est un état de
le confort thermique et l'efficience énergétique, tout en adoptant satisfaction ressenti par l'être humain vis a vis de son
une approche participative avec les habitants, car le facteur environnement thermique, Selon Milon Le confort thermique
anthropique détermine la consommation énergétique. se définit comme suit: "on est en état de confort thermique
Mots clés — bâtiments passif ; efficience énergétique ; enveloppe lorsque l'on ne perçoit ni sensation de froid ni sensation de
performante ; habitation ; confort thermique chaud». (Marité Milon, 2004)[2]. En effet, Le confort
thermique se définit comme un état d'équilibre et de
I. INTRODUCTION satisfaction général exprimée à l’égard de l’ambiance
La finalité de toute architecture est d'assurer le confort de thermique du milieu environnant.
l'être humain, a fin d'assouvir ses besoins vis a vis de Pour qu’une personne se sente confortable, trois conditions
l'environnement thermique il a recours a des technologies tel doivent être réunies : Le corps doit maintenir une température
que la climatisation au sein son habitation pour réguler son interne stable, La production de sueur ne doit pas être trop
environnement thermique. Elle met l’homme, en premier lieu, abondante et la température moyenne de la peau doit être
à l’abri de diverses agressions vis-à-vis desquelles il éprouve confortable, Aucune partie du corps ne doit être trop chaude ni
le besoin de se protéger,Le rôle de l'habitation est donc trop froide (inconfort local). (Commission de la Santé et de la
d'assurer a l'être humain la protection contre les aléas du climat Sécurité du Travail du Québec, 2002)[2], la notion de confort
pour lui assurer le confort thermique, Hensen (1991) signale est donc la corrélation entre trois paramètres bien définit: une
que le confort thermique est un état dans lequel il n'y a pas température interne stable, un sudation peu ou pas trop
d'impulsions qui conduit à corriger l’environnement de abondante, ainsi que l'absence d'inconfort local.
l’occupant par son comportement .Le but de la construction est
Copyright - 2019
ISSN 1737-9334 Page 35
Pour être en état de confort thermique, une personne doit individus, et également pour un même individu selon
porter une quantité raisonnable de vêtements Il peut être atteint différentes périodes.[5]
seulement lorsque la température, l'humidité et le mouvement
de l'air se situent à l'intérieur des limites de ce qu'on appelle III. L' EFFICIENCE ENERGETIQUE
la « zone de confort » (C.C.H.S.T - Centre Canadien Dans sa recherche de confort l'habitant adopte des
d'Hygiène et de Sécurité au Travail, 2002).[3] comportements énergivore au quotidien pour réguler son
Le confort thermique définit des plages de températures, de environnement thermique par le billet de la climatisation en
vitesse d 'écoulement d'air et des niveaux d'humidité dans les période estivale puisant dans énergies fossile polluantes par les
quelles les habitants ne ressentent pas le besoin de chauffage émissions de Co2 résultant de la consommation énergétique
ou de rafraichissement, et donc ressente le confort thermique des habitations. En domaine des logements, on converse sur
qui dépond des échanges de chaleur entre le corps humain et son isolation ou on peut comparer la consommation d'énergie
son environnement. Selon Givoni Cet "équilibre" n'est entre différents bâtiments, à la base d'une unité de (kW/m²/an).
possible que grâce à l'action constante de mécanismes de Notons bien qu’un bâtiment ne consomme beaucoup d’énergie
régulation. (Givoni, 1978), Givoni en 1978 a établi un c’est un bâtiment bien isolé.[6]. Selon hocine tabbouch et al de
diagramme psychrométrique ou il détermine une zone de nombreux bâtiments n'ont pas été conçus, conçus et construits
confort caractérisé par une température comprise entre 20 °C en tenant compte de la nature et du contexte spécifiques de
et 27°C, et une humidité relative comprise entre 20% et 28% l'environnement où ils se trouvent [7]. L'inadéquation de la
lors de la quelle l'indice d'insatisfaction est relativement nul, construction a son environnent engendre la hausse de
hors de la zone de confort Givoni propose des dispositifs consommation énergétique, c'est pourquoi il est primordiale
architecturaux pour améliorer le confort intérieur selon chaque d'assurer l'efficience énergétique. C'est le rapport entre
cas. Fig. 1 l’énergie utile produite par un système et l’énergie totale
consommée pour le faire bien fonctionner. Cette définition
Fig . 1 diagramme psychrométrique Givoni 1978 peut avoir un élargissement comme elle devient l'ensemble des
technologies et pratiques pour diminuer la consommation
énergétique jusqu'à la mise d'un nivaux de performance
équivalent dont l'objectif est de faire mieux avec moins.
Cette notion est souvent interprétée dans un sens plus large
pour désigner les technologies et pratiques permettant de
diminuer la consommation d’énergie tout en maintenant un
niveau de performance finale équivalent[9] .selon l'agence
international de l'énergie (IEA) l'efficience énergétique est
considérée comme un moyen pour contrôler et diminuer la
consommation énergétique.
l'efficience énergétique est donc l'utilisation de procédés, et de
technologies a fin de répondre au besoins, tout en économisant
le plus d'énergie possible, en effet il s'agit de faire mieux avec
moins d'énergie. De cette manière on assure des lieux de vie
La sensation de confort thermique est procurée par confortable, et on diminue la consommation énergétique et par
l'évacuation de la chaleur du corps. Les mouvement d'aire conséquent diminuer l'empreinte écologique. Il s'agit de
augmente les pertes de chaleur par convection et facilite fournir les mêmes prestations en consommant moins d'énergie
l'évaporation de l'humidité a la surface de la peau. [4] donc le ou même de faire mieux avec moins.
confort thermique est relié au échanges thermiques entre le Si on veut améliorer l'efficience énergétique il faudra adopter
corps humain et son environnement Dans une ambiance la démarche bioclimatique qui rétabli l'architecture a son
modérée, une personne perd sa chaleur selon 5 modes de rapport a son environnement sans le quel le dynamisme entre
transmission . le métabolisme, l'habillement, La température l'homme et l'environnement serai perturbé , il s'agit d'assurer
ambiante de l’air, L'humidité relative de l'air (HR), la vitesse l'équilibre entre la construction et son environnement en
d'écoulement d'air qui ne doit pas dépasser 0.2 m/s. plaçant l'homme au centre de cette boucle.il faut aussi agir sur
Le notion de confort thermique ne dépond pas que de facteurs l'enveloppe extérieur des bâtiments qui est le fil conducteur
physiques objectifs ( température, humidité relative, vitesse vers la consommation énergétique car en diminuant la surface
d'écoulement d'air), mais aussi de facteur subjectif relié a la déprédative on diminue la consommation énergétique.
psychologie des habitants suivant les mœurs, la culture, et le C'est pourquoi l'amélioration de l'efficience énergétique dans
métabolisme, en fait il n'ya pas de conditions types pour un le secteur du bâtiment est un axe prioritaire a développer a fin
homme type. A l’inverse des réponses physiologiques qui de préserver l'environnement, de répondre au besoins, et
peuvent être mesurées de manière objective, la détermination garantir la capacité des générations futures a répondre au leurs.
des réponses subjectives des sens dépend de la propre Pour se faire il est impératif de considéré les deux volet
évaluation de la personne soumise à un environnement donné. habitant et technologie comme étant un phénomène
Cette évaluation n’est pas unique, mais varie avec les d'interaction entre deux.

ISSN 1737-9334
Les technologies de l'énergie et le comportement des Il est aussi important de noter que le bâtiment passif est adapté
occupants ont été traités comme des acteurs distincts dans les a son environnement, les matériaux utilisés, l'orientation,
domaines de confort, de l'ingénierie énergétique et des l'exploitation des énergies renouvelable, pour assurer un
domaines sociaux les efforts récents tentent de les relier, en environnement thermique confortable.
considérant la consommation énergétique comme résultat de Le bâtiment passif a lui seules ne peu assurer l'efficience
l'activité humaine dans un contexte d'interrelations entre énergétique si on ne sensibilise pas les habitants a adopté une
utilisateurs, technologies, compétences, contextes sociaux[10] mode de vie et des pratiques quotidiennes résiliente résilientes,
par le chois des appareils électroménager de catégorie A et B
IV. LE BATIMENT PASSIF sur l'étiquette écologique, et a les utiliser de façon rationnel
le bâtiment passif instaure un mode de vie relativement sans tomber dans l'excès.
résiliant car il diminue la consommation énergétiques en
utilisant l'inertie thermique et l'adéquation des matériaux de V. METHODOLOGIE DE LA RECHERCHE
construction, ainsi que l'orientation et la forme de l'habitation, La corrélation entre l'étude de Givoni qui détermine de larges
en plus d'une panoplie de procédés passifs . Utilisé depuis des zone de sous chauffe et de surchauffe, et l'approche par mesure
siècles dans certains points de la planète Les stratégies in situ a l'intérieur de l'habitation : température, humidité
passives apportent des réponses à la fois économique et Loïc relative, vitesse de l'air, révèlent des température intérieur au-
adapté au milieu où la construction doit avoir le jour Elles delà des plages de confort ce qui cause des situations de stress
prennent en compte la température l'orientation et le vent Il est aux habitants.
aussi possible de profiter de l'orographie local pour assurer la
régulation thermique des intérieurs où faire usage de matériaux A. Critères de choix de la ville Guelma
offrant une masse thermique élevée ou de haute qualité Guelma ville du nord est algérien caractérisée par le climat de
d'isolation.[11]. l’arrière littoral montagne (Zone B) climat semi-humide.
Une maison passive est un bâtiment avec un climat intérieur Déterminé par des hivers plus froids et plus longs et des étés
agréable en hiver comme en été sans installation de chauffage chauds et moins humides que ceux du littoral, bénéficiant d'un
ou de refroidissement conventionnelle »[12] . Cela est rendu grand potentiel durable par sa couverture forestière et agricole,
possible par la réduction drastique des besoins en énergie de ainsi que par le bombardement solaire qui couvre 243.3
chauffage, réalisée essentiellement par des mesures hectare par an aussi par son potentiel hydraulique avec 264.96
architecturales et constructives. millions m³ d'eau mobilisable avec de multiple source d'eau
Le batiment passif rpond a des exigences technique pour chaude, c'est aussi une ville enclavé entouré de quatre
pouvoir etre qualifié de tel. montagne ce qui définit les conditions de vie et le microclimat
 L'isolation thermique renforcée, fenêtres de grande de la ville.
qualité Malgré tout ce potentiel les habitant souffre situation stress,
 La suppression des ponts thermiques comme le démontre l'étude de Givoni selon medejelakh.D il
 L'étanchéité à l'air excellente y'a une zone de sous chauffe très froide pour les mois de
 La ventilation double flux (avec récupération de janvier et février (T°entre 4.5 et 5°C) ainsi qu'une saison de
chaleur) sous-chauffe froide de décembre à février, et de mars (T° entre
 La captation optimale mais passive de l'énergie 5 et 15°C).
solaire et des calories du sol
 La limitation des consommations des appareils
ménagers
Le bâtiment passif ne réalise pas d'économie sue le niveau de
confort mais au contraire l'augmente tout en diminuant la
consommation énergétique, par son enveloppe performante,
selon Selon la norme 90.1-2016 Les nouvelles exigences
d'enveloppe technique comprennent:
Exigences obligatoires pour la vérification de l'enveloppe,
avec un intérêt particulier aux infiltration d'air réduite et des
exigences accrues fuite d'air vers les portes d'enroulement
aériennes.
 Exigences prescriptives plus strictes pour la Fig . 2 situation géographique de
construction métallique les toits et les murs, la Guelma ( source: google )
fenestration et les portes opaques.
 Amélioration de la clarté des définitions des murs B. critère de choix du cas d'étude
extérieurs, de l'orientation du bâtiment, et la clarté
autour de la valeur R. La villa Mokhnache est choisi pour évaluer le confort
thermique des habitants selon différentes pièces, orientation, et
 Nouvelles exigences basées sur l'ajout de la zone
climatique 0.[13]
Copyright - 2019
Page 37
ISSN 1737-9334
fonction, aussi le fait que ce soit mon domicile garanti le Tableau.4 Plan 1er étage
déroulement de la compagne de mesure sans soucis.
Tableau. 1fiche technique cas d'étude
Villa Villa R+2

Mokhnache
Surface 480 m2
totale
Surface 224 m2
bâtie
Coefficient 0.46
d'occupation
au sol
La compagne de mesure s'est déroulée sur 03 jours successives

lors de la vague de canicule, Les pièces sujets a l'investigation
on était choisi selon l'étage, l'orientations et ,l'occupation.
Les plans de la villa dans les tableaux suivants. Tableau.5 Plan 2em étage
Tableau. 2 point de mesures
Tableau.3 Plan RDC

A. Matériel utilisé
Le relevait des températures et de l'humidité sont faite par
un Thermo-Hygromètres portables étanches HANNA 9564.
Fig. 5
La vitesse de l'air était mesuré par un anémomètre a
hélice BA16 Trotec. Fig.6
Fig . 3 Thermo- Fig . 6 anémomètre a

Hygromètre HANNA hélice BA16 Trotec
9565

ISSN 1737-9334
Fig .7 description fonctionnel HANNA 9565
Mesures T°/C et H/% séjour
Fig . 8 Description fonctionnelle anémomètre a

hélice BA16 Trotec
VI. RESULTATS
Les résultats de la compagne de mesure sont traduit dans les
graphs suivants.
La vitesse d'écoulement de l'air lors de toute la duré de la
compagne de mesure est resté constante inferieur a 0.76 m/s
car l'anémomètre ne détecte pas le mouvement d'air inferieur a
0.76 m/s.
Mesures T°/C et H/% cuisine
Copyright - 2019
Page 39
ISSN 1737-9334
Mesures T°/C et H/% chambre 2
Mesures T°/C et H/% chambre parentale

ISSN 1737-9334
Mesures T°/C et H/% salon

ISSN 1737-9334
tandis que le maximum est a34.4°C a 18h le 07/07/2019 au

séjour 2em étage.
Mesures T°/C et H/% séjour
Généralement les températures les plus fraiches sont noté au
RDC tout d'abord car le sous sol participe au rafraichissement
de l'air, mais aussi car l'air chaud est plus lourd que l'aire frai,
c'est pourquoi le RDC est l'étage le plus frai, et les
température les plus élevées sont noté au 2em étage, puis dans
la chambre 1 au 1er étage en raison de l'orientation plein ouest
mais aussi a cause du plafond exposé au rayonnement solaire
car au dessus il s'agit d'un terrasse accessible.
Pour améliorer les conditions thermiques les habitants ont
recours aux climatiseurs dont l'utilisation varie selon l'activité
anthropique. Tableau 6
ETAGE PIECE HORRAIRE REMARQUE
USAGE
Cuisine / Le taux d'insatisfaits
RDC Séjour 16h-19h est a 100%,car tout
rare les habitant atteste
Chambre 00h-12h une situation
1ER 1/2 15h-17h d'inconfort et estime
ETAGE Quotidien qu'il est impossible
Chambre 23h-07h d'utiliser les espaces
parentale Quotidien de vie sans les
Salon 13h-19h climatiseurs pour
22h-10h réguler la température
2EM Séjour 22h-09h 25°C au salon, 24°C
ETAGE rare a chambre 1,et 22°C
chambre1 / a la chambre 2.
VIII. RECOMMANDATION
La solution pour assurer l'efficience énergétique est le
bâtiment passif.
Selon diagramme psychrométrique Givoni 1978 les
corrections a prévoir sont les suivantes:
 masse thermique élevée avec ventilation nocturne
 masse thermique élevé
 ventilation naturel et mécanique
combiné a d'autre stratégies de refroidissement passif on

peut diminuer l'usage des climatiseurs et voir même s'en passé
en adoptant les stratégies suivantes
 emplois de matériaux a forte inertie thermique
pour assurer le déphasage de températures.
 une enveloppe performante
VII. DISCUSSION  une bonne isolation thermique
Les résultats de la compagne de mesure révèlent que la  étanchéité a l'air :réduire les infiltration
vitesse d'écoulement d'air est inferieur a 0.76 m/s et comprise d'air peut économiser 5-40% des besoins
dans la zone de confort 0.2m/s, les taux d'humidité intérieur de refroidissement.
sont comprissent dans les limites de la zone de confort définit  une exposition optimum au soleil
par Givoni 1986,(T° 20°C/27°C°) (H 20M%-80%) le taux  obstruction du vitrage
minimum est de 23.65% a 12h le 07/07/2019 a la chambre  utilisation du double vitrage au lieu du vitrage
parentale, et le taux maximum est de 44.6% a 22h le simple.
05/07/2019 en cuisine, quand a la température c'est toujours  protection solaires fixe et mobile
au-delà de la température de la zone de confort avec un  parement clair pour les murs extérieur, et la toiture
minimum e 27.10°C a 10h au séjour RDC le 06/07/2019,

ISSN 1737-9334
 diminuer le coefficient de forme K pour diminuer REFERENCES

la surface déprédative de l'enveloppe. [1] (Marité Milon, 2004)
 stratégie bioclimatique du froid. [2] (Commission de la Santé et de la Sécurité du Travail du Québec, 2002)
[3] (C.C.H.S.T - Centre Canadien d'Hygiène et de Sécurité au Travail,
 Eviter: enveloppe performante 2002
 Protéger: protections solaires [4] LIEBARD.allin et DE HERDE.andré. traité d'architecture et
 Refroidir: point d'eau, végétation d'urbanisme bioclimatique _ Paris : le moniteur, 2000 p 83 b
 Minimiser les apports internes [5] M’Sellem* , D.alkama. « Le confort thermique entre perception et
évaluation par les techniques d’analyse bioclimatique - Cas des lieux de
 Dissiper: ventilation naturel travail dans les milieux arides à climat chaud et sec ». Revue des
Energies Renouvelables Vol. 12 N°3 (2009) 471 – 488
IX. CONCLUSION [6] BENHARRA Houda: Impact de l'orientation sur la consommation
La situation de stresse climatique ressentie par les habitants énergétique dans le bâtiment. -Cas des zones arides et chaudes-.
Université Mohamed Khider - Biskra pp (12)
est du au défaut de construction, le mauvais choix des [7] Hocine Tebbouche a , Ammar Bouchair b,Said Grimes c . « Towards an
matériaux de construction, l'orientation, les pont thermique, environmental approach for the sustainability ofbuildings in Algeria.
ainsi qu'au manque de geste écologique de la part des habitants ».Energy Procedia 119 (2017) 98-110 pp 05
tant dans le choix des appareils électro ménager, des dispositif [8] BENHARRA Houda: Impact de l'orientation sur la consommation
énergétique dans le bâtiment. -Cas des zones arides et chaudes-.
pour l'éclairage et la climatisation sans se soucier de l'étiquette Université Mohamed Khider - Biskra pp (13)
écologique, que quand dans leurs emplois, c'est pourquoi [9] BENHARRA Houda: Impact de l'orientation sur la consommation
avant de penser a apporter des corrections d'ordre technique il énergétique dans le bâtiment. -Cas des zones arides et chaudes-.
est primordiale d'agir sur le facteur anthropique. Université Mohamed Khider - Biskra pp (13)
[10] Marco A. Ortiz∗ , Stanley R. Kurvers, Philomena M. Bluyssen. « A
L'insatisfaction des habitants vis a vis de leurs review of comfort, health, and energy use: Understanding daily energy
environnement thermique se traduit par une consommation use and wellbeing for the development of a new approach to study
énergétique importantes qui nécessite d'être raisonner, On peut comfort ». Energy and Buildings 152 (2017) 323–335 pp 05
assurer l'efficience énergétique par le billet du bâtiment passif [11] FRECHMANN, Kolon et al. Architecture et énergie un enjeu pour
l'avenir _ Paris : place des victoires pp10
et son enveloppe performante car c'est la partie exposé au aléas [12] Site internet de la Plate-forme Maison Passive asbl
de la nature, elle joue le rôle de régulateur déterminant (www.maisonpassive.be)
l'ambiance intérieur de l'habitation en effet elle crée un micro [13] Ingy El-Darwish *, Mohamed Gomaa. « Retrofitting strategy for
climat avec ses propres paramètres de température, humidité, building envelopes to achieve energy efficiency ». Alexandria
Engineering Journal (2017).pp 02
ventilation, éclairement, ensoleillement... ect qu'elle détermine
par sa forme, ses composants, et son orientation, donc c'est
grace au batiment passif adapté au climat semi-aride de la ville
de Guelma qu'on peut y assurer l'efficience énergétique.

ISSN 1737-9334
Le CuInSe2 en couches minces est une structure de

base pour une application photovoltaïque.
GHALMI Leila 1, BENSMAINE Souhila 2
Département de Physique, Université Abou Bekr Belkaid -Tlemcen, Algérie.
1,2
Unité de Recherche Matériaux et Énergies Renouvelables (URMER). Université de Tlemcen (UABT). BP n°119. Tlemcen.
13000. Algérie.
1 Ghalmileila2017@gmail.com
2 s_bensmaine@yahoo.fr
Résumé — Le CIS (CuInSe2) est l’un des matériaux semi- par exemple, la mis au point en 2007 d’une unité de production
conducteurs les plus prometteurs pour la production des piles avoisinant les 100 MW annuelle [2]. La question qui se pose,
solaires en couches minces, efficaces et moins coûteuse. Des qu’est ce qui doit être fait pour que la technologie CIS atteint
instituts de recherches ont annoncé leurs plans pour les chaines une grande échelle de production d’électricité photovoltaïque?
de production de celles-ci. Mais la technologie CIS requière
encore davantage d’attention concernant l’industrialisation et la II. LE MATÉRIAU CIS
stabilité à long terme. Dans ce travail nous nous intéressons à
l’étude des couches minces du composé ternaire CuInSe2, qui ont
révélé récemment beaucoup d’intérêt dans la communauté des Le diséléniure de Cuivre et d'Indium (CIS) est un matériau
sciences des matériaux et leurs utilisations pour la fabrication des composé de type I-III-VI, de structure chalcopyrite très
cellules photovoltaïques. En effet, dans notre étude nous avons prometteur, puisque le rendement théorique de
donné les principales caractéristiques des matériaux les plus l’hétérojonction (n) CdS-(p) CuInSe2 se situe autour de 25 %
utilisés en couche minces en guise de comparaison avec le [3]. Les cellules à base de composés chalcopyrites quaternaires
CuInSe2, en suite nous nous sommes intéressés aux types de
du type Cu (Ga, In) (Se, S)2 ont récemment atteint des
cellules qui ont été élaboré à base de ce semi-conducteur et les
améliorations qui sont apportées à celles-ci en illustrant leurs rendements de 17 % [4]. Les principales améliorations sont
performances. Par ailleurs, nous avons pris le soin de présenter venues de l’élargissement de la bande interdite du CuInSe2 (Eg
les méthodes et techniques d’élaboration de ce dernier. = 1.02 eV) par l’utilisation d’alliages du type CuGaSe2 et
Finalement, nous présentons les résultats de simulation d’une CuInS2.
cellule typique à base de CIS en utilisant le logiciel de calcul
SCAP-1D tout en variant les deux paramètres suivants :
❖ La variation de la température de 300°C à 350°C
❖ La variation de l’épaisseur de la couche CIS de 0,5µm à
3,0µm
Mots clés — CuInSe2, couches minces, SCAPS-1D, Application
photovoltaïque.
I. INTRODUCTION
L’un des principaux obstacles rencontrés pour que le
photovoltaïque devienne plus populaire est le cout du watt
produit [1]. Il n’est pas encore concurrentiel à celui de
l’électricité produite par les méthodes conventionnelles.
Durant les dernières décennies, la course des différents
groupes de recherches et de sociétés vers la diminution du prix Fig 1 : Schéma d’une coupe transversale d’une hétérojonction CIS /CdS[3].
du watt d’électricité produite par le photovoltaïque, conduit au
développement de nouveaux semi-conducteurs considérés III. TYPES DES CELLULES SOLAIRES A BASE DE CIS
actuellement comme matériaux promotteurs pour les Toutes les cellules sont créées à partir de matériaux semi-
applications photovoltaïques. Parmi ces candidats, les couches conducteurs. La plupart du temps, c’est le silicium (Si) qui est
minces de di-séléniure de Cuivre et d’Indium (CuInSe2). Par utilisé mais on rencontre aussi le sulfure de cadmium (CdS), le
rapport au silicium, les cellules solaires à base de ce semi- tellure de cadmium (CdTe), des alliages de cuivre indium et
conducteur présenteraient une efficacité élevée, une excellente sélénium (CIS), …
stabilité, une résistance aux irradiations et une durée de vie plus
longue. Cependant, après plus de trente années de recherches
et de développement les cellules à base de CIS viennent tout
récemment de passer au stade industriel et commercial, avec

ISSN 1737-9334
A. Cellule solaire à homojonction: MOCVD d’une couche conductrice et transparente de ZnO

Dans le domaine des couches minces de CIS, on constate (1.5 µm) dopée en général à l’aluminium.
qu’il y’a très peu de travaux de recherches pour la réalisation Ainsi c’est l’utilisation de la fenêtre optique ZnO/CdS mince
d’homojonction basées sur ces matériaux. Selon S.Kohiki et al qui, en améliorant la réponse de la photopile dans le bleu
[5], la conductivité des couches minces Cu-In-Se, dont la (figure 3) a permis d’augmenter considérablement le
structure est du type chalcopyrite, change du type n vers le type rendement [9]. C’est donc cette technologie qui a été adopté
p par implantation d’ions des éléments de la colonne V (N, P, par la plupart des équipes travaillant sur les photopiles à base
Sb ou Bi) dans la mesure ou ces ions substituent les atomes de de CIS.
Se dans la couche. A titre d’exemple, une diode à
homojonction p-n a été fabriquée par implantation ionique
d’ions d’azote à 450°C dans une couche CIS (type n), pendant
une minute [6]. Un rendement de 0.35 % sous éclairement AM
1.0 (100 mW/cm2), sur une surface de 4 mm2, a été obtenu sur
cette photodiode. On note également les paramètres suivants :
Vco=0.098 V, Jcc = 9.11 mA/cm2 et un facteur de forme
FF=39 %.
B. Cellules solaire à heterojunctions
Etant donné que les matériaux en général sélectionnés
comme couches minces pour les cellules solaires sont des
semi-conducteurs à gap direct et souvent d’un seul type (n ou
p), les hétérostructures sont généralement utilisées pour
réaliser une barrière de potentiel [7]. Fig 3 : Comparaison des spectres de rendement quantiques des dispositifs
à base de couches minces CuInSe2 avec :
1) La photopile primitive (dite Boeing) : 1) Une fenêtre optique standard (Zn,Cd)S.
Ceux sont les travaux effectués chez Boeing qui ont définit 2) Une fenêtre multicouches ZnO/CdS.
la structure des hétérojonctions à base de CIS (figure 2) [8]. Au
cours des premières recherches l’absorbeur CIS de type p IV. TECHNIQUES D’ELABORATION DU CIS
étaient formés de deux couches évaporées sur un contact Les méthodes utilisées pour le dépôt des couches minces
arrière de Molybdène (le Molybdène étant utilisé à cause de sa peuvent être divisées en deux groupes basés sur la nature du
non réactivité avec le CIS) le rapport Cu/In dans le film de CIS processus physique ou chimique du dépôt. Les méthodes
étant le paramètre clé pour les performances de la cellule, la physiques incluent le dépôt à vapeur physique dite "PVD" (Eng.
composition de la couche composite finale était un léger excès Physical Vapor Deposition), l'ablation laser, épitaxie par jet
de Cu par rapport à In près du contact au Mo et l’inverse c.-à- moléculaire "MBE", et la pulvérisation Cathodique
d. un excès d’In par rapport au Cu sur le dessus (donc près de "Sputtering" [10]. Les méthodes chimiques incluent, les
la jonction). méthodes à dépôt en phase gazeuse et les techniques à solution
(Figure 4). Les méthodes en phase gazeuse sont le dépôt à
vapeur chimique (Chemical Vapor Deposition CVD) [11] et
l'épitaxie à couche atomique (Atomic Layer Epitaxy ALE),
tandis que les méthodes de spray pyrolys, sol-gel, spin coating
et dip-coating emploient des solutions comme précurseurs [12].
Fig 2 : Coupe transversale d’une cellule type CuInSe2/(Cd,Zn)S. a)structure

conventionnelles avec des couches actives et fenêtres optique en bicouches.
b) fenêtre optique ZnO/CdS [7]
2) La cellule améliorée : C’est le groupe Arco Solar

(Brevet Choudary et al en 1986) qui eut l’idée de remplacer le
Fig 4 : Présentation des principaux procédés de dépôt de couches minces
CdZnS par deux couches : une couche très fine (50 à100 Å) de [13].
CdS déposée sur le CIS actif par dépôt en bain chimique (CBD :
Chemical Bath Déposition), suivie du dépôt par sputtering ou

ISSN 1737-9334
V. SIMULATION NUMERIQUE DE LA CELLULE CIS PAR Le courant de court-circuit peut être calculé à partir du
SCAPS-1D. rendement quantique.
∞
Nous nous introduisons le concept de simulation numérique
des semi-conducteurs, particulièrement l’application sur les 𝑰𝒄𝒄 = ∫ 𝑸𝑬(𝝀) . 𝝋(𝝀)𝒅𝑬
𝟎
cellules solaires du type CuInSe2 en utilisant le logiciel de
calcul SCAPS-1D, développé spécialement pour la simulation VI. SIMULATION, RESULTANTS ET DISCUSSION
des cellules à base de CIS.
Les figures 6.a et 6.b montrent les résultats de la simulation
A. La cellule utilisée, (ZnO/CdS/CuInSe2) : des caractéristiques I-V sur la cellule ZnO/CdS/CIS à
On se propose de simuler les propriétés essentielles d’une l’obscurité et sous éclairement respectivement, avec :
cellule à base de CuInSe2 ayant une structure composée d’un
oxyde transparent conducteur (OTC) de type n-ZnO, une
couche tampon n-CdS et d’une couche absorbante de type p-
CuInSe2 (Figure 5).
n-CdS (0,02µm)
Fig 6.a : Caractéristique I-V à Fig 6.b : Caractéristique I-V

l’obscurité sous éclairement.
n-ZnO (0,3µm) p-CuInSe2 (3µm) Une tension à circuit ouvert Vco=0.6617 V, un courant de
court-circuit Jcc= 32.89 mA/cm2, un facteur de forme FF=81.30
et un rendement de conversion η=17.69%.
Fig 5 : La cellule ZnO/CdS/ CuInSe2
Le rendement quantique de la cellule est donné sur la figure
B. Les grandeurs technologiques d’une cellule solaire :
7 montre un pic entre les valeurs 90 -100 % puis une
La caractéristique I=f(V) noté aussi « I-V » est une fonction
diminution due aux pertes par recombinaison.
qui décrit le comportement de la photopile, et à partir de
laquelle plusieurs paramètres relatifs à la cellule peuvent être
calculés.
1) Le courant de court-circuit Icc : correspond à
l’intensité du courant de court-circuit, c.à.d. V=0
2) La tension à circuit-ouvert Vco : correspond à la
tension en circuit ouvert c.à.d. I=0
3) La puissance maximale Pm : C’est la puissance
maximale de sortie correspondant au point pour lequel le
produit de là de la tension et de l’intensité du courant soit
maximal
Pm=Vm. Im
4) Le facteur de forme : C’est le rapport entre la
puissance maximale délivrée par la cellule et le produit entre Fig 7 : Rendement quantique QE de la cellule ZnO/CdS/CuInSe2.
Icc et Vco correspondant à la puissance maximale idéale A. Influence de la température de fonctionnement sur les
𝑽 𝒎 𝑰𝒎 𝑷𝒎 paramètres de la cellule :
𝑭𝑭 = =
𝑽𝒄𝒐 𝑰𝒄𝒄 𝑽𝒄𝒐 𝑰𝒄𝒄 La température joue un rôle très important sur les
5) Le rondement η : mesure le taux de conversion performances de la cellule solaire. La valeur qui
énergétique correspond à un fonctionnement optimal de la cellule est
𝑷𝒎 𝑽𝒄𝒐𝑰𝒄𝒄 T=300 K. les résultats de la simulation montrent que le
𝜼= =
𝑷𝒊𝒏 𝑷𝒊𝒏 paramètre le plus affecté par la variation de la température
6) Le rendement quantique QE : C’est le nombre de est la tension de circuit ouvert Vco, son influence sur cette
pairs électron-trou photo-générés par le nombre de photons dernière est montrée sur la figure 8.
incidents sur la cellule. Il est mesuré en fonction de la
longueur d’onde.

ISSN 1737-9334
Tableau 1: Simulations I-V pour différentes valeurs de la

température.
T (K) Voc (V) Jcc(mA/cm2) FF η (%)
300 0.6617 32.8906 81.30 17.69
310 0.6415 32.9087 80.83 16.98
320 0.6217 32.9266 79.88 16.35
330 0.6013 32.9446 79.05 15.66
340 0.5813 32.9229 78.28 15.00
350 0.5609 32.9816 77.42 14.32
Les valeurs obtenues pour la tension à circuit ouvert

et le courant de court-circuit sont portés sur le tableau 1.
Fig 9 : Rendement quantique QE pour différentes épaisseurs de la couche

CuInSe2. L’épaisseur varie de 3.0μm (rouge) jusqu’à 0.5μm (jaune).
Tableau 2 : Paramètres I-V pour différentes épaisseurs de la

couche CuInSe2.
Epaisseur de la
Jcc
couche CuInSe2 Vco (V) FF η (%)
(mA/cm2)
(μm)
3.0 0.6617 32.8906 81.30 17.69
2.5 0.6595 32.6243 81.18 17.47
2.0 0.6554 32.1956 80.92 17.07
1.5 0.6488 31.4742 80.76 16.49
1.0 0.6405 30.5196 80.11 15.36
0.5 0.6090 27.6702 76.98 12.97
Fig 8 : Les simulations I-V pour différentes valeurs de la température. T
varies de 300K (rouge), à 350K (jaune).
Pour les petites épaisseurs, le phénomène de générations de
En effet, lorsque la température augmente la bande de gap pairs électrons-trous se passe près des interfaces avec le CIS et
devient plus étroite et ceci accélère le phénomène de le contact arrière (forte densité de défauts donc de centres de
recombinaison des pairs électrons-trous entre la bande de recombinaisons) ce qui réduit le nombre de pairs électrons-
conduction et la bande de valence. trous générées et donc du rendement quantique QE.
La forte réduction de Vco conduit à une diminution du
facteur de forme FF ainsi que du rendement de conversion η de En revanche, pour des couches plus épaisses la génération
la cellule solaire. se passe loin des interfaces c’est pourquoi on observe une
B. Influence de l’épaisseur de la couche absorbante CuInSe2 : augmentation du rendement. Une épaisseur de 3.0 μm parait
L’effet de la variation de l’épaisseur de la couche être adéquate pour avoir un rendement optimal (17.69 %).
absorbante CIS et sur les grandeurs I-V de la cellule est donné
VII. CONCLUSIONS
par le tableau 2.
Dans le but de pouvoir améliorer les performances
Les deux valeurs Vco et Jcc diminuent lorsque l’épaisseur de des dispositifs photovoltaïques, ceci nécessite parfois la
la couche absorbante diminue car celle-ci absorbe de moins en présentation de modèles théoriques pour la simulation
moins de photons ayant les longues longueurs d’onde, ce qui numérique. Dans ce travail nous avons introduit le concept
influence négativement sur le taux de génération de pairs de simulation numérique sur la cellule solaire à base de CIS,
électrons-trous. en prenant le soin de présenter le logiciel de calcul utilisé
(SCAPS-1D).
L’effet de la réduction de l’épaisseur de la couche
absorbante CIS sur le rendement quantique observé dans En effet, les propriétés intrinsèques de la cellule sont
l’intervalle des longueurs d’onde λ= 450 à 1100 nm et est très importantes et influencent largement les performances
montré sur la figure 9. de celle-ci telle que la température de fonctionnement, les
épaisseurs des couches composant la cellule. Ceci est
illustré dans les caractéristiques I-V et les courbes obtenues
pour le rendement quantique, avec un rendement de
conversion optimal égal à 17.69% pour une couche
absorbante de 3μm.

ISSN 1737-9334
Finalement, les rendements atteints par les cellules CIS et [4] S. Beghdadi. Etude des performances d’une cellule solaire de type Cu
les performances dont il a fait preuve font de lui l’un des (In,Ga)Se2 (CIGS). Université Abou Bekr Belkaid. Tlemcen. (2016).
[5] S. Kohiki, M.Nishitani, T. Negami, K. Nishikora and T.Kerao. appl
matériaux qui contribueront fortement au potentiel énergétique phys letters. (1991).
solaire du futur. [6] M. Nishitani et al. J. app phys. (1993).
[7] W. H. Bloss and H. W. Schock proc 8th.E.C photovoltaic solar energy
REMERCIEMENTS conf Florence. (1988).
[8] R. A. Mickelsen and W. S. Chen, Proc.16th IEEE photov. SPEC. Conf.
Nos remerciements vont à toute l’équipe de la (1982).
division 1 Matériaux Photovoltaïques (MPV) au sein de notre [9] K. Mitchell, C. Eberspatcher, G. Ermer and D. Pier.proc 20th IEEE
Photovoltaic solar energy conf New York. (1988).
Unité de Recherche Matériaux et Energies Renouvelables [10] S. Bensmaine. Dépôt de couche minces d’oxyde de zinc avec axe c
(URMER). Tlemcen, Algérie. incliné par la pulvérisation magnétron et la réalisation des dispositifs à
ondes acoustique de surface SAW. Doctorat. (2008).
REFERENCES [11] Y. Natsume, H. Sakata, T. Hirayama and H.J. Yanagita, Phys. Status
Solidi. (1995).
[1] I. Dehmous, Etude des couches minces de CuInSe2 pour une application
[12] H. Belkacemi, Dépôt et Caractérisations des Couches minces de ZnO
photovoltaïque. Université Abou Bekr Belkaid, Tlemcen. (2012).
Avec Différentes Sources de Zinc, Université Mohamed kheider-Biskra.
[2] D. Iincot et J. F. Guillemdes, Institut de recherche et développement sur
(2011).
l'énergie photovoltaïque.
[13] A.R. Elshabini-Riad and F.D. Barlow, Thin Film Technology
[3] Ingrid Repins ; Miguel A. Contreras , Brian Egaas , Clay DeHart ,John
Handbook, III series McGraw-Hill. (1998).
Scharf ,Craig L. Perkins, Bobby To, Rommel Noufi. 19·9%‐efficient
ZnO/CdS/CuInGaSe2 solar cell with 81·2% fill factor. (2008).

ISSN 1737-9334
Les intentions bioclimatiques des architectes de la

Reconstruction tunisienne d’après-guerre : cas des
logements collectifs d’El Menzah 1
Salma Gharbi#1
#
Université de Carthage, ERA, Ecole Nationale d’Architecture et d’Urbanisme de Tunis, Tunisie
elgharbisalma@hotmail.com
Résumé— Nous nous intéressons dans cet article à un legs

architectural et urbain de la période de la Reconstruction A. La période hygiéniste
tunisienne d’après-guerre et spécifiquement au quartier d’El La période hygiéniste remonte à la fin du 19ème siècle, où
Menzah 1 à Tunis. La particularité de ce quartier est l’approche un courant de pensée est né, suite aux conditions de vie
environnementaliste de ses concepteurs, se souciant de chaotiques de la population et les découvertes médicales,
l’intégration des logements dans leur environnement naturel et
comme celle de Louis Pasteur. La tâche de l’architecte fut
microclimatique. Il est question ainsi d’étudier les intentions
bioclimatiques de l’équipe de Bernard Zehrfuss, architecte en importante et devait trouver sa juste place entre le savoir
tête de l’équipe reconstructrice et de questionner leur pertinence médical et le pouvoir politique [1]. L’architecte était
de termes d’ambiances lumineuse et aéraulique. redevable de répondre aux exigences des hygiénistes d’une
part, et de respecter le règlement régissant la construction
d’une autre part. La conception architecturale devait être
Keywords— intégration, conditions microclimatiques, intentions
régie par une réglementation administrative qui fut établie
bioclimatiques, ambiance lumineuse, ambiance aéraulique
suite à un état des lieux de l’insalubrité des logements. Ainsi,
I. INTRODUCTION une fois le nouveau règlement énoncé, on pouvait voir
l’application des théories architecturales et urbaines de la
Peut-on parler d’intentions bioclimatiques liées à la
manière de faire de l’équipe de la période de la période hygiéniste sur le bâti et cela s’est manifesté en trois
Reconstruction tunisiennes d’après-guerre ? ce grandes théories : i. Théories des expositions et des
questionnement est basé sur les choix conceptuels entrepris orientations ii. Théories des gabarits des rues iii. Théories des
par l’équipe de Bernard Zehrfuss architecte en chef de la formes et dispositifs des bâtiments, appliquées principalement
Reconstruction. En effet, même si de formation académique par les deux architectes ; Auguste Rey à travers son axe
héliothermique et le plafond lumineux, et Henri Sauvage à
classique, ces architectes ont été longtemps bercé par les
idéaux modernistes du Mouvement Moderne prônant travers ses immeubles à gradins.
plusieurs recommandations liées à l’environnement naturel et B. La période moderne
la prise en compte de son cadre physique. Ainsi nous allons
L’habitat collectif à l’époque moderne rencontre une toute
entreprendre dans cet article de cerner les intentions
autre tournure, celle convoitant « l’idéal solaire de
environnementalistes de ces concepteurs qu’ils soient
l’hygiénisme, conjugué à l’évolution de l’architecture
relatives à la dimension urbaine du quartier d’El Menzah 1,
moderne » [2]. L’intérêt porté à la lumière naturelle n’est plus
œuvre phare de l’équipe, ou à sa dimension architecturale
le même, suite aux découvertes médicales de la vaccination
relative aux immeubles collectifs dressés en barre.
contre certaines maladies. D’autre part, l’utilisation des
II. CADRE HISTORIQUE DE LA PENSEE nouveaux matériaux et les nouvelles techniques de
ENVIRONNEMENTALISTE construction, développées suite à la révolution industrielle, a
permis la conception d’immeubles avec une surface
Pour mettre les choses dans leur contexte, nous devons
importante de pans de verre. Permettant certes l’ouverture sur
remonter aux prémices de la pensée bioclimatique et
l’extérieur, cette surface va aussi créer une surchauffe à
précisément à sa relation avec l’évolution de l’habitat. Ainsi
l’intérieur des logements. C’est ainsi qu’une nouvelle
nous développons dans ce qui suit l’intérêt des concepteurs
problématique se présente aux architectes : celle de la
aux données de l’environnement physique relatives à la
maitrise du flux solaire. Le principal architecte à s’être
lumière et à la ventilation naturelle et ceci dans la période
intéressé à cette question, reste sans doute Le Corbusier [3]. Il
hygiéniste, moderne et contemporaine. La période moderne
essaya de trouver des solutions architecturales et inventa son
dans laquelle s’inscrit l’œuvre de la Reconstruction est une
brise-soleil vers la fin des années 1920. Plus tard, il met aussi
charnière dans l’évolution de l’habitat collectif depuis les
en œuvre un autre dispositif de protection solaire, la loggia,
recommandations hygiénistes de la fin du 19 ème siècle et
comme un espace tampon entre l’intérieur et l’extérieur
début du 20ème siècle jusqu’aux nouvelles tendances
pouvant temporiser les rayons solaires. Il conjugua ces deux
écologistes du 21ème siècle (fig. 1).
dispositifs dans le projet de l’unité d’habitation de Marseille.
Copyright - 2019
Page 49
ISSN 1737-9334
En plus de l’important intérêt accordé à l’ensoleillement, les portaient l’intérêt plus à la ventilation naturelle des logements,
concepteurs du logement social de la période moderne se à leur ensoleillement et spécifiquement à la protection solaire
focalisaient aussi sur une autre composante physique qui est des espaces de vie [6].
l’aération. L’aération ou la ventilation naturelle [4], assurera
le renouvellement de l’air à l’intérieur des pièces qui A. La dimension bioclimatique dans l’espace urbain du
permettra l’extraction de l’air pollué, malodorant et vicié. Un quartier d’El Menzah
air purifié, préservera ainsi un climat intérieur sain, doté Conçue par l’équipe de Zehrfuss, la forme urbaine est
d’une température et d’une humidité adéquates à composée de plusieurs entités qui permettent d’intégrer ce
l’épanouissement de l’habitant. quartier dans la sphère de la pensée moderniste. Le système
viaire, la circulation piétonnière, les espaces verts et le cadre
C. La période contemporaine bâti interagissent avec un intérêt majeur porté à la question de
Le langage contemporain adopté aujourd’hui au niveau des l’ensoleillement et de l’aération, créant un espace urbain
immeubles collectifs s’oppose au langage rationnel, répétitif structuré et aéré. Le quartier présente une densité faible ou
et parfois systématique de celui produit jusqu’aux années moyenne, où les espaces vides enveloppent les espaces bâtis,
soixante-dix. La mise en œuvre de logiciels de rendu et de avec une forte concentration de l’espace libre au centre même
conception en trois dimensions a facilité aux architectes, du quartier. La composante végétale est prédominante dans le
l’accès à un nouveau langage usant d’une volumétrie plus quartier. Elle est présente dans les espaces communs, et
complexe, d’une gamme de couleurs et de textures plus matérialisée dans un îlot dédié à cette fonction. La végétation
recherchée et d’un style formel plus original répondant au domine aussi les limites de l’espace dédié aux logements
développement des techniques de la fin de ce siècle. Outre cet collectifs et ponctue plusieurs rues du quartier. L’habitat
intérêt à la forme, une autre démarche environnementale est individuel est aussi délimité sur le pourtour par un espace
née suite à la naissance de la mouvance écologique et de la extérieur végétal parfois très dense, dissimulant tout le bâti.
politique du développement durable qui tend à trouver un L’ombrage par la végétation peut influer sur la qualité de vie,
juste équilibre entre les données environnementales, sociales la création d’ambiance urbaine et les situations de confort
et économiques. Ainsi, plusieurs architectes conçoivent thermique et visuel des usagers [7]. Une autre caractéristique
aujourd’hui des logements collectifs et même sociaux qui du quartier réside dans le caractère multifonctionnel du lieu
prônent un respect au site, en minimisant les déperditions faisant en sorte d’offrir à part la composante résidentielle,
thermiques et en utilisant des matériaux écologiques. Cette d’autres fonctions comme les équipements scolaires (école
démarche s’intègre dans la création des éco-quartiers qui primaire et collège), les espaces commerciaux regroupés au
offrent un cadre de vie participant au « vivre ensemble » et à centre, et les espaces de jeux pour les enfants. Le tout donne
« la mixité sociale » [5] mais aussi qui répondent aux besoins naissance à un quartier autonome, où l’habitant minimise les
environnementaux en promouvant une bonne gestion des trajets pour les courses quotidiennes à la manière des
ressources et en participant à l’économie de la ville. quartiers écologiques d’aujourd’hui, où la dimension sociale
est favorisée afin de créer une mixité et un partage entre les
voisins.
B. La dimension bioclimatique dans les logements du
quartier d’El Menzah
La tranche du quartier étudiée est constituée de 9
immeubles en barre. Tous les immeubles répondent à la
configuration d’une double exposition, avec une façade, où
s’ouvrent sur des loggias, les espaces de vie orientés Sud-Est
Fig. 1 : Intérêt face à l'ensoleillement et à la ventilation suivant l'époque, et une autre façade sur laquelle ouvrent les espaces de
source : Auteure
services donnant dans le cas du Balkiss sur une coursive
III. LES INTENTIONS BIOCLIMATIQUES DES ARCHITECTES orientée Nord-Ouest. Les façades des immeubles du quartier
DE LA RECONSTRUCTION TUNISIENNE
d’El Menzah 1, ne présentent pas d’ornementation spécifique.
Elles répondent au caractère sobre et épuré de l’architecture
Nous avons entrepris avant d’identifier les intentions moderniste. Cependant, nous remarquons que le traitement de
conceptuelles des architectes un travail de collecte la façade sud-est est différent de celui nord-ouest, où ont été
archivistique basé sur des extraits de revues internationales utilisés des dispositifs architecturaux énoncés par l’équipe de
des années cinquante où l’équipe de Bernard Zehrfuss Zehrfuss comme des dispositifs de protection solaire, se
exprimait leur manière de concevoir les espaces urbains et résumant à des loggias brise-soleil, des lamelles verticales,
architecturaux. Nous avons analysé le contenu et pu dégager des loggias avec claustras et des avancés de maçonnerie.
un champ sémantique relatif à trois dimensions : celle D’autres intentions conceptuelles permettent d’apporter une
lumineuse, aéraulique et visuelle. Ainsi pour la composante plus-value aux espaces en termes de ventilation,
urbaine trois volets intéressaient les architectes, à savoir les d’ensoleillement et d’apport thermique à savoir la double
voiries, les espaces verts et la disposition des immeubles
selon l’orientation. Pour la composante architecturale, ils

ISSN 1737-9334
exposition et la configuration traversante des logements et la

construction en mur de moellon de 50 cm de largeur.
IV. PERTINENCE DES CHOIX CONCEPTUELS DES

ARCHITECTES DE LA RECONSTRUCTION
TUNISIENNE
Nous présentons dans ce qui suit les résultats de l’étude
quantitative (simulations numériques d’une typologie de
logement à l’aide du logiciel Solène) et qualitative (analyse
des dispositifs architecturaux et enquête auprès des habitants).
Rappelons que notre travail vise à évaluer la pertinence des
choix bioclimatiques entrepris dans la conception du quartier
et des logements de la première tranche d’El Menzah 1.
Fig. 2: Niveau d'éclairement le 21 juin à midi du logement de l’immeuble
Balkiss dans son état initial/modifié, source : Auteure
A. Structure de l’espace urbain, orientation et disposition
des immeubles .
La morphologie en barres des immeubles conditionne
l’ambiance dans le quartier. Qu’elle soit lumineuse ou
aéraulique, l’ambiance est directement liée à la structure
urbaine du quartier. La forme en barres possède un double
rôle ambiantal sur les espaces publics. D’une part, elle
masque les rayons solaires et permette d’ombrager un grand
tronçon des voies ; d’autre part, elle engendre des
phénomènes aérauliques pouvant parfois être inconfortables
en saison hivernale. La composante végétale, Qu’elle longe
les rues, ponctue certaines zones ou se densifie dans d’autres,
contribue à la caractérisation ambiantale du quartier en termes
de dimension visuelle, lumineuse, thermique. Cette présence Fig. 3: Ecart au niveau de l'éclairement dû à la suppression de la loggia,
fait du quartier d’El Menzah un « îlot de fraicheur » [8] dans source: Auteure
la ville de Tunis.
B. La façade épaisse L’équipe de Zehrfuss, soucieuse des conditions climatiques
La façade épaisse ou poreuse est une spécificité de du site, a su transposer des doctrines urbanistiques et
l’architecture de la Reconstruction tunisienne. En effet, un architecturales interagissant avec l’environnement physique
interstice est créé entre la limite extérieure du bâtiment et la du quartier, des immeubles et des logements. Toutefois, cette
limite intérieure du logement se déclinant en un dispositif interaction est aujourd’hui déréglée par les modifications
appelé loggia. Ce dispositif constitué aussi de lamelles apportées par les habitants sur ce cadre bâti. En ce qui
horizontales constitue un dispositif de protection solaire concerne la dimension aéraulique, nous pouvons constater
filtrant les rayonnements solaires. qu’en été, la façade de la loggia est directement exposée aux
vents chauds de l’Est-Sud-Est. Un courant d’air dégagerait
1) Rôle du dispositif loggia brise-soleil : Nous avons l’air vicié de l’autre côté. La loggia atténuerait l’effet de ces
simulé, à l’aide du logiciel de simulation solaire Solène, les vents chauds. En hiver, les vents dominants sont plutôt du
durées d’ensoleillement et les niveaux d’éclairement dans les côté de la façade postérieure Nord-Ouest. L’air neuf passe
logements (fig. 2). Ces simulations prennent appui sur deux donc à travers les ouvertures de la cuisine ou de la salle de
supports : les logements comme conçu par les architectes, et bain et s’évacuerait du côté de la loggia. Un courant d’air est
les logements comme modifiés par les habitants. La ressenti par les usagers, donnée confirmée par nos enquêtes
confrontation des deux états, nous a permis de dégager un auprès des habitants.
écart (fig. 3) dans la qualité lumineuse et hygrothermique qui 2) trous d’aération : Un autre dispositif d’aération
en est conséquente. En effet la suppression de la loggia, ayant naturelle a été utilisé par les concepteurs. Ce dispositif est
pour rôle de masquer les rayons solaires parfois pénibles en matérialisé par de petits creux (fig. 4) dans la paroi extérieure
été, dans les logements engendre des zones où les durées du logement du modèle de l’immeuble Virgile donnant sur la
d’ensoleillement et les niveaux d’éclairements sont assez loggia orientée sud-est, ainsi que sur la façade opposée Nord-
élevés. Ceci peut engendrer un inconfort lumineux et Ouest. Ce procédé, souvent utilisé dans l’architecture
hygrothermique pour les habitants vernaculaire, permet l’apport d’un effet-cheminée dans les
pièces et faisant circuler l’air de bas en haut.

ISSN 1737-9334
REFERENCES
[1] Mory, P. (2001). Architecture et hygiénisme à Paris au début du XXe
siècle : l'architecte entre savoir médical et pouvoir politique. Dans P.
Bourdelais, Les hygiénistes ; enjeux, modèles et pratiques (p. 540).
Paris : Belin.
[2] Harzallah, A. (2007, janvier 11). Émergence et évolution des
préconisations solaires dans les théories architecturales et urbaines en
France, de la seconde moitié du XIXe siècle à la deuxième guerre
mondiale. Nantes.
Fig. 4: Trous d'aération permettant la ventilation naturelle de l'espace de vie
[3] Corbusier, L. (1955). Œuvres Complètes 1946-1952 (Vol. 5). Zürich :
dans le logement de l’immeuble Virgile: source: Auteure
GirsbergerR.
C. L’effet traversant des logements [4] Bonneaud, F. (1998). Ventilation naturelle des bâtiments collectifs
dans les départements d'outre-mer : influences de la forme et de la
La double exposition des appartements sur deux façades composante de l'habitat. Nantes-Grenoble : Ecole d'Architecture de
opposées sud-est nord-ouest (fig. 5), permet un apport solaire Nantes, Ecole d'Architecture de Grenoble.
continu à l’intérieur des espaces tout au long de la journée et [5] Athamena, K. (2012). Modélisation et simulation des microclimats
une très bonne ventilation naturelle permettant le dégagement urbains : étude de l’impact de la morphologie urbaine sur le confort
dans les espaces extérieurs. Cas des éco-quartiers. Nantes : Thèse de
de l’air vicié. Les habitants attestent que leur logement est Doctorat, spécialité architecture, Ecole Centrale de Nantes.
« bien conçu » en termes d’exposition et d’ouvertures sur [6] Gharbi S. (2017) ; Identification et caractérisation des ambiances d’un
l’extérieur. La question d’humidité est complétement écartée quartier d’habitat collectif de la période de la Reconstruction
dans leurs propos ce qui affirme encore plus le bon tunisienne de l’après-guerre : cas du quartier d’El Menzah 1 (ex
Crémieux-ville) à Tunis » (thèse de doctorat non publiée), Ecole
ensoleillement et la bonne aération des logements. Nationale d’Architecture et d’urbanisme de Tunis, Tunisie
[7] Gharbi, S., Belakehal, A. (2017). Les ambiances du quartier d’El
Menzah : impacts de la morphologie urbaine et de la végétation sur le
vécu des habitants, XXXème colloque de l’Association Internationale
de Climatologie (A.I.C), « Climat, ville et environnement », Sfax,
Tunisie.
[8] Carrega, P., Charfi, S., & Dahech, S. (2010). Apport de l’imagerie
spatiale dans l’étude multiscalaire de l’ilot de chaleur urbain à Tunis.
23ème Colloque de l’Association Internationale de Climatologie.
Rennes.
Fig. 5: La double exposition sud-est, nord-ouest des logements du quartier

d'El Menzah, source: Auteure
V. CONCLUSIONS
Les architectes de la Reconstruction tunisienne d’après-
guerre étaient soucieux de concevoir, dans le quartier d’El
Menzah 1 à Tunis, un cadre urbain et des logements en
concordance avec les données de l’environnement physique
et spécifiquement celles relatives aux stimuli de la lumière et
du vent. Nous avons ainsi pu dégager plusieurs intentions
conceptuelles et dispositifs engendrant des climats intérieurs
et extérieurs en adéquation avec le soucis environnemental et
usager des habitants comme : la forme de l’espace urbain et
la disposition des immeubles, l’effet masque et régulateur de
la végétation, l’effet isolant des murs extérieurs en moellon
de 50cm, les cages d’escalier ouvertes et aérées, le traitement
en claustra, l’effet masque de la loggia et autres dispositifs
de protection solaire, l’effet traversant des appartements et
les fentes d’aération hautes. Ces intentions ont été évaluées
qualitativement et quantitativement et les résultats de cette
appréciation, attestent d’une pertinence dans le choix
conceptuel des architectes faisant du quartier d’El Menzah,
œuvre de la période moderne de l’histoire de la Tunisie, une
référence en termes d’habitabilité environnementale qui
pourrait bien sure être sujet de réhabilitation énergétique afin
de le rehausser aux standards écologiques internationaux de
notre époque.

ISSN 1737-9334
La reproduction de l'habitation troglodyte

vernaculaire aux temps 'modernes'. Cas d'un
projet d'habitat du 'plan de Constantine' a
Biskra (sud-est Algérien).
Asma Achraf ZENDAGUI 1 et Azeddine BELAKEHAL2
Université de Biskra, Département d’architecture, Laboratoire LACOMOFA, Algérie
1e-mail: zendaguiasma2014@gmail.com
Abstract— Dans ce travail, il y aura lieu de reconstituer et examiner les ambiances générées à
l’intérieur d’une habitation semi troglodytique, proposée pour confronter les conditions sévères du
climat chaud et sec. L’étude architecturale du projet comportant cette proposition date de la fin des
années cinquante. Elle fait partie d’une grande opération de logement dans l’histoire de l’Algérie
appelée ‘plan de Constantine. Le projet s’inscrit dans le registre de l’architecture du mouvement
moderne. l’intérêt d’étudier cette facette environnementale physique des ambiances de ce projet
revient au fait que cette opération de logement par souci économique interdit toute solution de
conditionnement d’air par voie mécanique, entrainant une dépense notable d’énergie. L’ambiance y
est donc tout simplement naturelle.
Keywords— Habitat troglodytique, Architecture passive, Architecture moderne, Climat chaud aride et sec, Ciel
Clair Ensoleillé, Environnements thermiques et lumineux.
I. INTRODUCTION
Ces dernières décennies, la relation entre l’héritage II. L'HABITAT TROGLODYTIQUE
architectural et l’approche environnementale est L’un des premiers abris humains sur terre. C’est
mise en exergue de jour en jour. La mise en valeur une excavation dans la terre ou bien dans la roche
des architectures anciennes commence à prendre qui possède plusieurs formes, et ça selon les
une nouvelle dimension dans le domaine de caractéristiques géologiques du site. On trouve
recherche scientifique ; celle de source inédite de l’habitat excavé en puits ou en grotte (architecture
solutions passives. Dans notre cas d’étude traditionnelle Méditerranéenne, Euromed heritage).
l’habitation troglodytique vernaculaire des anciens L’habitat troglodytique figure parmi les
temps servait de référence aux architectes de la constructions spécifiques aux régions désertiques,
modernité, pour arriver à une architecture de simple elles sont nombreuses au Sahara citant les
confort avec le moindre cout. A notre ère, elle exemples du sud de l’Algérie (Tassili au Hoggar),
présente une référence à la fois économique, aussi ceux du sud de la Tunisie (Matmata…) .les
environnementale et durable. Ce qui va être illustré conditions qu’elles offrent sont peut-être la raison
dans ce texte au moyen d’une étude ambiantale. de leur utilisation dans cette catégorie de zones.
En premier, on commence par une Cette typologie d’habitat a suscité l’intérêt des
contextualisation historique dans laquelle on architectes des temps modernes, tel une solution
présentera un aperçu sur l’habitat troglodytique, et architecturale économique en matière de
les données historiques de genèse du projet objet consommation énergétique. Cet intérêt a été
d’étude, Qui sera objet d’une reconstitution soulevé lors d’une étude monographique d’un
virtuelle.qui nous donne à la fin un objet projet d’habitat inscrit dans le répertoire du
expérimental. Après avoir définir le Protocol mouvement moderne. Ou le maitre d’œuvre s’est
expérimental, la variante reconstituée va être mis à la recherche du type d’habitat le plus
soumise à une simulation des environnements approprié aux contraintes climatiques et
physiques générés. économiques imposés.

ISSN 1737-9334
III. CONTEXTE HISTORIQUE DU PROJET IV. MODELISATION 3D

Pour effectuer la modélisation, il nous a fallu le
passage par trois étapes principales.
Le projet objet d’étude, est un projet d’habitat de i) l’assemblage du corpus documentaire. Ii) la
simple confort « évolutif », qui date de la fin des numérisation des documents. Iii) la modélisation
années cinquante. S’insère dans le cadre 3D.
d'opérations connues par le ‘plan de Constantine’
(1958). Ce plan prévoyait environ 210 000
logements à travers le pays. Le projet a été destiné à
la ville de Biskra, située au sud-est de l’Algérie,
connue par son climat aride et chaud. L’étude été
soumise à des contraintes économiques spéciales tel
que l’interdiction de toutes sortes de climatisation
mécanique. La mission de l’étude architecturale a
été attribuée à l’architecte George Henri
Pingusson, une figure représentative du
mouvement moderne en France. Il a proposé le
recours aux solutions des temps des anciens, qui se
basent sur les principes de physique élémentaire
utilisés pour la réfrigération économique, indiquant Figure 2 Plan RDC
qu’elle est la meilleure réponse aux conditions
climatiques sévères des milieux désertiques. Après
son examen des bâtisses préexistantes et du
contexte géo climatique. Il a suggéré une variante
de logement qui rappel l’habitat troglodytique, il l’a
nommé aussi habitat « géothermique » ou
« souterrain ». Une reproduction moderniste d’un
type d’habitat très ancien. L’étude du projet a
atteint la phase d’exécution, mais sa réalisation fût
abandonnée en raison de l'indépendance de
l'Algérie. Pour cela, nous avons mené, en une
première étape, une reconstitution virtuelle de
l’habitation troglodytique proposée dans ce projet,
et ensuite, soumis cette habitation à une évaluation
des environnements physiques y générés par le
Figure 3 Plan du sous sol
moyen du logiciel ECOTECT.
Figure 1 photo de maquette du projet 540lgts Figure 4 Coupe

ISSN 1737-9334
de comprendre comment s’interagissent les

composantes de cette expérimentation ambiantale.
a. Contexte géo climatique
La ville de Biskra est située à 470 KM au Sud- Est

d'Alger, sur les lignes 34,48° de latitude nord, et
une longitude de 5,73°Est, et à 87 m d’altitude, elle
se situe dans la zone D qualifiée comme étant
semi-aride à climat chaud et sec.
b. Le macroclimat lumineux
Figure 5 Axonométrie montrant l'ancrage de l'habitation Le macroclimat lumineux de la ville de Biskra

dans le sol, source Auteur
ressemble en beaucoup d’aspects à celui des
régions désertiques vue sa situation au nord du
Grand Sahara. Un ciel clair régnant Presque toute
l’année et dont la luminance atteint 100.000 lux
(Satel-light.com).
La portion des jours nuageux est d’environ 6.66%
et les jours ensoleillés constituent une portion
d’environ 73%.
c. Outil de simulation thermique
Pour la simulation de l’environnement thermique

on a utilisé Autodesk® Ecotect™ Analysis 2011
qui est un logiciel de simulation et de conception
au même temps. Développé par Autodesk, Inc, il
combine un modeleur 3D avec des analyses
solaires, thermiques, acoustiques et de coût.
Figure 6 Aspect volumétrique de l'habitation semi
troglodytique, source Auteur
d. Outil de simulation lumineuse
Le programme Radiance (2.0 BETA) est un

ensemble de programmes libres développé par Greg
Ward Larsaon (1985 –1997) au laboratoire national
de Laurent Berkeley pour l'analyse de haute
V. SIMULATION DE précision et la visualisation de l'éclairage dans la
L’ENVIRONNEMENT THERMIQUE conception. C’est le programme qui donne des
ET LUMINEUX valeurs simulées les plus proches de celles réelles
spécifiques aux régions à ciel clair ensoleillé.
La conformation architecturale de l’habitation étant
définie , on procède à l’expérimentation de cette V. RESULTATS DE LA SIMULATION ET
conformation dans son contexte naturel, par le LEURS INTERPRETATIONS
moyen d’une simulation informatique à l’intérieur
des logements. L’objectif est d’explorer les types a. la Simulation de l’environnement
d’ambiances générées à l’intérieur des logements. thermique
Soulignons qu’on va se limiter à simuler deux
environnements physiques, ceux lumineux et La construction simplifiée du logement dans le
thermique. Ces deux dimensions sensorielles sont logiciel Ecotect est la première étape de la
choisies pour être étudiées suite aux intentions de simulation. Pour cela, il fallait définir sa situation,
G.H.Pingusson pour un logement saharien et sa son orientation puis les matériaux de construction
croyance en l’existence d’une proposition idéale prévus (tableau 1). Aussi il est à souligner que la
dans ce contexte. maison est ancrée dans le sol à une profondeur de
En premier, on doit définir un protocole 2.5m.
expérimental dans lequel on expose toutes les Dans la simulation, la terre assiette de la maison
conditions génératrices des deux environnements est saisie comme une zone qui a les caractéristiques
physiques. Ce protocole est structuré sur la base
d’un schéma illustratif (figure 7). Qui nous permet

ISSN 1737-9334
thermiques du sol dans le site de l’étude, un sol Ou bien entre 26° à 28°c selon le CSTB (1958). Il
gypseux en surface, comprend de gros galets s’avère donc, qu’une climatisation mécanique est
d’oued à 4 - 5 mètres de profondeur. nécessaire pour le confort des habitants.
Introduction des matériaux de construction La simulation de l’environnement lumineux

naturel
La construction du modèle avec Ecotect, s’effectue
avec le principe de zones thermiques. Chaque zone L’environnement lumineux est le deuxième
représente un espace architectural bien déterminé environnement physique ciblé par la simulation
avec ses cloisons et ouvertures chaque élément dans ce travail L’opération de simulation, sera
constitutif d’une zone possède des caractéristiques exécutée par le logiciel Radiance. Les étapes de la
thermo physiques qui ont été bien définies. La simulation, cerne en premier l’importation du
simulation des températures cernera la journée de model 3D de l’habitation convoitées. Ensuite, on
pic de chaleur d’été 22 juillet, les résultats sont fixera les points de vision de l’observateur, ce qui
montrés en forme de graphs. permet de constituer une série de scènes. Ces
derniers, seront transférés à Radiance, après un
Températures dans le séjour ensemble de réglages on lance le rendu. Les
Dans le séjour, les températures journalières résultats de la simulation seront présentés sous
varient entre 38c° et 41c°. La température forme d’images comportant les quantités de
minimale est enregistrée de 6:00 h à 10:00 h du lumière.
matin, la maximale à 23:00 h.
On remarque une stabilité des températures d’une Lumière naturelle dans le séjour
moyenne de 39.1c° malgré le changement de la Dans le séjour situé au rez de chaussées, la
température à l’extérieur. Un écart important entre séquence simulée visée depuis le canapé en
T° intérieur et T° de l’extérieur atteint 9.2c° à temps de repos à 16h, et en position assise
14:00h. dirigée vers la table à manger. La source
lumineuse englobe la porte et deux ouvertures
Températures dans la chambre ajourées donnant sur la loggia. Cette séquence
Les températures journalières dans la chambre est rendue en images contenants les luminances
01 située au sous sol, oscillent entre 38c° et 40c°. « human sensitivity ».
La température minimale est enregistrée de 4:00 h à
6:00 h du matin. Elle s’amplifie jusqu’à une valeur
maximale de 40c°, enregistrée à midi et reste
presque la même toute l’après midi jusqu’à la nuit.
On remarque une moyenne de température
journalière de 39.1c° stable malgré le changement
des températures à l’extérieur. Un écart important
dans les heures de pic de chaleur de l’après midi il
atteint 8.7c°
Interprétation et discussion des résultats
On remarque que les températures d’air ambiant

dans la maison, demeurent relativement stables
Figure 7Séquence 16:00h Avril séjour / luminances
pendant toute la journée. La température « human sensitivity »
minimale est enregistrée durant les premières
heures du matin, avec une valeur entre 37 et 38 C°.
La plus élevée est de 42 c°, et ce durant les
premières heures nocturnes. Une moyenne de Lumière naturelle dans La chambre
température journalière est enregistrée de 39.1 C°. Se situe au sous-sol. Elle s’ouvre sur le patio à ciel
ouvert, par la porte et une fenêtre ajourée.
De plus, on note un intervalle important entre la La séquence simulée montre la perception de la
température à l’intérieur de l’habitation et celle à mère au moment du réveil le matin. Elle est assise
l’extérieur qui atteint 9.5 C° pendant les heurs les en face des ouvertures de la chambre, source de
plus chaudes de la journée. Ce qui prouve que la la lumière naturelle. L’image illustrée montre les
maison tienne sa température interne tout au long luminances aux quelles l’homme est sensible,
de la journée. Cependant, la température de confort « human sensitivity ».
n’est pas atteinte. En effet la marge de confort
chez les habitants de la ville de Biskra se situe
entre 18° et 30°c selon Khoukhi et Fezzoui (2012).
Page 56
Copyright - 2019
ISSN 1737-9334
trouvé un écart important entre la température

int/ext qui atteint 9c°, c'est-à-dire en été l’intérieur
de l’habitation est plus frais que l’extérieur. Facteur
qui peut aider à minimiser la consommation
énergétique pour le conditionnement.
Bibliographie
Archives,Apc Biskra. (1960). 540 Logements,Ophlm

Constantine.
Figure 8 Séquence 07:00h Avril chambre / luminances
« human sensitivity » Belakehal, A. (2013). De La Notion D’ambiance.
Courrier Du Savoir . [3]
Interprétation et discussion des résultats
Çelik, Z. (1997). Urban Forms And Colonial
Dans le séjour et la chambre les luminances varient Confrontations Algiers Under French Rule. Berkeley ·
entre 50 et 950 cd/m². Le rapport entre ces deux Los Angeles · Oxford: University Of California Press.
valeurs est conforme aux valeurs recommandées [7]
qui sont de 1/20 à 1/40 (0.05 à 0.025) pour le
Cstb. (1958, Juillet). Comment Construire Au Sahara?
macro-champ visuel. Confort:Traditions,Matériaux,Technologies (N°6) [3]
L’ancrage de la maison dans le sol, est une sorte de
protection des rayons solaires directs. En plus, le Iceb. (2014). Les Guides Bio-Tech L'éclairage Naturel.
fait de surhausser les ouvertures avec des France. [6]
décrochements, a engendré la création d’un Khoukhi, M. F. (2012). Thermal Comfort Design Of
environnement lumineux discret à l’intérieur des Traditional Houses In Hot Dry Region Of Algeria.
espaces. La quantité lumineuse est rationnelle, International Journal Of Energy And Environmental
voire peu suffisante parfois, l’éclairage doit être Engineering [10]
renforcé avec la lumière artificielle
Meliouh, F. (1998). Pratiques Domestiques Feminines
Les séquences simulées de l’environnement Dans Le Logement Collectif:Espaces Et Confort. Biskra:
lumineux à l’intérieur de l’habitation sous un ciel Centre Universitaire M.Khider,Institut D'architecture.
clair ensoleillé, révèlent une multitude [10]
d’atmosphères lumineuses qui différent d’un
espace à l’autre dans la maison. Généralement Picard, A. (S.D.). Architecture Et Urbanisme En Algérie
elles sont caractérisées par un aspect lumineux 1830-1962. Figures De L'orientalisme En Architecture ,
doux et parfois sombre. Cette qualification est Pp. 128-129. [3]
confirmée par la simulation des quantités de Ragon, M. (1986). Histoire Mondiale De L'architecture
Et De L'urbanisme Modernes. Casterman. [1]
lumière naturelle. Ces dernières ont été parfois
Revue Des Mondes Musulmans Et De La Méditéranée.
conformes aux recommandations prescrites et (1994). Figures De L'orientalisme En Architecture.
satisfait le déroulement de l’activité humaine, et (N°73-74). [3]
parfois la quantité de lumière s’avère insuffisante Texier, S. (2006). Geaoge- Henri Pingusson
pour l’accomplissement des taches, donc Architecte,1894-1978- La Poétique Pour Doctrine. Paris:
l’éclairage artificiel sera nécessaire Verdier [1]
Http://Www.Meda-
Corpus.Net/Frn/Portails/Pdf/F1/Es_T01.Pdf[7]
VI. CONCLUSION
Les résultats ont montré que l’environnement

physique lumineux est caractérisé par un aspect
lumineux doux et parfois sombre. Cette
qualification est confirmée par la simulation des
quantités de lumière naturelle. Ces dernières ont
été parfois conformes aux recommandations
prescrites et satisfait le déroulement de l’activité
humaine, et parfois la quantité de lumière s’avère
insuffisante pour l’accomplissement des taches,
donc l’éclairage artificiel sera nécessaire.
L’environnement physique thermique est, quant à
lui, ne montrait pas des températures de confort
sans conditionnement mécanique. Certes, on a

ISSN 1737-9334
Comparaison entre les modules d’élasticité d’un

béton autoplaçant exposé à hautes température par
deux méthodes de fréquence de résonance
BENOUDJAFER imane 1, BENOUDJAFER ibtissam2
1
Laboratoire de Mécanique et structures(LMS), Université de Tahri Mohammed Béchar(UTMB)
imenousse@yahoo.fr
1
Faculté de Technologie, Université de Tahri Mohammed Béchar(UTMB)
Abstract— Cet article a pour but principal de présenter une les conditions fonctionnels qu’il rencontra dans le temps [16].
comparaison entre le module d’élasticité déduit par deux La connaissance de l'état des ouvrages et leur évaluation est
méthodes de fréquence de résonance. Le matériau testé est un donc un enjeu majeur pour décider des interventions de
béton autoplaçant soumis aux différents stades de températures réfection et pour en évaluer le rapport coût / utilité. C’est
élevées jusqu’à 550°C. Cette comparaison repose principalement
pourquoi, une gamme d’essais in situ, appelés essais en place,
sur l’utilisation de variétés procédées expérimentaux, conduisant
à l’extraction des paramètres modaux à partir de l’analyse de a été développée [1]. Ces essais sont traditionnellement
fréquence de résonance. Les résultats expérimentaux obtenus appelés essais non destructifs, étant entendu que quelques
montrent que les modules d’élasticités, déterminés par les deux désordres mineurs peuvent être apportés à la structure, sans
méthodes expérimentales sont proches. toutefois modifier ni sa performance ni son apparence [11].
Une importante caractéristique des essais non destructifs est
Keywords— béton autoplaçant, fréquence de résonance, qu’ils peuvent être refaits au même endroit ou presque, ce qui
température, module d’elasticité, permet de suivre les changements des propriétés du béton
dans le temps [4]. L’utilisation d’essais non destructifs
I. INTRODUCTION entraine une plus grande sécurité et une meilleure
Les travaux menés sur maniabilité des bétons ont permis la planification de la construction, de sorte que l’on peut
mise au point des autoplaçants (BAP) [1]. Cette nouvelle progresser plus rapidement et plus économiquement [5].
catégorie de bétons a été développée afin d’obtenir un
matériau facile à mettre en œuvre, sans faire appel à la C’est dans ce contexte de nécessités de compréhensions,
vibration, présentant de nombreux avantages sur les chantiers d’une part, des données sur le comportement du béton
et en usines de préfabrication [2]. Ces bétons sont formulés autoplaçant exposé à la température pour la prévision de la
différemment des bétons ordinaires : ils contiennent au moins sûreté des bâtiments [2] et des constructions en réponse à
un adjuvant chimique et un ajout minéral en proportions bien certains accidents ou à états particuliers durant son
précises pour satisfaire les exigences de maniabilité et de fonctionnement, et d’autre part, de l’évolution des propriétés
stabilité [15]. intrinsèques du matériau béton et de celle de l’évolution
générale des structures , qu’émergent et s’affirment les
En effet, les essais de laboratoires sont habituellement techniques non destructives de contrôles [3]. Ce qui rend
effectués à une température contrôlée, normalement constante, possible l’étude des évolutions d’une ou plusieurs
généralement choisie entre 20 et 25°C, de sorte que la plupart caractéristiques dans le temps [10]. Le contrôle non destructif
des informations de base sur les propriétés des bétons durcis offre également la possibilité d’évaluer la qualité d’un béton
sont fondées sur la connaissance de comportement du béton à dans un ouvrage existant et suivre la santé du béton dans le
ces températures [11]. Cependant, en pratique, les ouvrages temps [21]. Cette étude a donc pour but principal de
en béton sont toujours susceptibles d’être soumis à une action caractériser un béton autoplaçant formulé à base des
thermique : pour des raisons fonctionnelles, climatiques, matériaux locaux ayant subi des températures élevées jusqu’à
accidentelles et en fin particulières (démolition partielles ou 550°C, à partir de deux méthodes d’analyse mécanique de
locale par procédé thermique) [14]. fréquence de résonance.
II. PROCEDURE EXPERIMENTATLE

Par conséquent, le béton demeure en pratique à différentes
températures, ce qui implique une réflexion sur l'opportunité A. Matériaux formulation et mise en ouvre du béton
de l'état des ouvrages portés aux différents gradients Dans le cadre de ce travail expérimental, le choix des
thermiques et sur la notion de risque lié à leur maintien dans matériaux (agrégats et filler calcaire) a été dicté par leurs

ISSN 1737-9334
disponibilités de la région d’Ain Sefra. Ces matériaux  L'essai d'étalement (BAP) :

proviennent de la carrière de Djebel Ben Hendjir qui se Cet essai permet de démouler un cône normalisé dite « DIN »
trouve sur l’axe routier Ain Sefra et Sfissifa. Nous avons (ou cône d'Abrams), de béton frais testé afin de mesurer le
utilisés deux classes granulaires de gravier (3/8,8/15) de diamètre de la galette formulée. Il faut rappeler, qu’un béton
natures concassés, un seul type de sable (0/3), les fillers est considéré comme autoplaçant, quand la galette mesurée
calcaires qui représentes l’addition calcaire, aussi un ciment forme un diamètre compris entre 60 et 75 centimètres. Le
Portland résistant aux sulfates de classe 42,5 « CPJ-CEM II/A temps d'écoulement intermédiaires a été mesuré permettant
42,5 », et un superplastifiant « SIKAVISCOCRETE 3040 », à d’atteindre un diamètre de 50 cm, noté t50 [AFGC, 00].
haut performance du troisième génération exempté de
chlorures prêt à l’emploi, Il est conçu à base de  L'essai de la boite en « L » :
polycarboxylates modifiés qui se présent sous la forme d’un L’essai de la boite en « L » a été effectué pour tester la
liquide brun. mobilité du béton confectionné et de vérifier que sa mise en
TABLE I
place. pour ce faire, plus de 13 litres de béton ont été mis
caractéristiques physiques et chimiques des fillers calcaires dans la partie verticale puis laisser pour reposer pendant 1
minute. La trappe a été ensuite levée et le béton a été écoulé
Carbonates Insolubles Sulfates Chlorures dans la partie horizontale de la boite à travers le ferraillage. A
86.20 13,70 nul 0.10 la fin de l’écoulement, les différents hauteurs H1, H2 ont été
Bleu de méthylène Densité Densité mesurés pour exprimer en terme de taux de remplissage
Fines<63µm donné par la fraction de (H2/H1).
« VB » absolue apparente
0.41 2.58 1.62 82
 L'essai de stabilité au tamis :
La méthode de formulation utilisée pour concevoir le BAP cet essai consiste à mesurer la proportion π d'éléments fins de
testé dans cette étude expérimentale est une méthode béton (laitance) passant au travers d'un tamis d'ouverture 5
empirique basée sur l'optimisation du squelette granulaire , mm afin de qualifier le béton vis-à-vis la ségrégation et de
permettant , d’une part, d’identifier les différents constituants, déduire si le béton testé possède une stabilité satisfaisante ou
et d’autre part, d’alimenter une base de données efficace non. Les critères d'acceptabilité d'une formulation de BAP
accédant à l’étude des mélanges de diverses combinaisons des sont divisés en trois classes[AFGC, 00]:
constituants afin d’obtenir un béton qui sera satisfaisant à un  0 % <π (laitance) < 15 % : stabilité satisfaisante ;
cahier des charges donné.  15 % <π (laitance) < 30 % : stabilité critique ;
TABLEAU 2  π (laitance) > 30 % : stabilité très mauvaise
Composition de béton formulé (ségrégation systématique, béton inutilisable).
Composition
Composant Les résultats des essais de caractérisation réalisés sur le
(Kg/m3)
Gravier 8/15 582.50 béton confectionné sont récapitulées dans le Tableau 3, avec
Gravier 3/8 311.10 toutes les mesures expérimentales nécessaires pour réaliser
Sable 892.42
les différents essais de caractérisation. À noter que les
Filler calcaire 128
Ciment 350 valeurs reportées dans le tableau constituent une moyenne.
Viscocrete 3040 4.10
L’eau efficace 157.5
Tableau 3
Masse volumique
2425.62
du béton frais Caractérisation du béton étudié.
La caractérisation à l'état frais des bétons a été limitée aux
essais recommandés par l'AFGC , à savoir l’étalement au Dmoy (cm) 62
Etalement
cône, l’écoulement à la boîte en « L » et la stabilité au tamis, T50(s) 2,10
Figure 1.
H2/H1 0,83
L-box
T40(s) 3,60
Stabilité au tamis π (%) 6,32
Il faut mentionner, qu’avant chaque essai, les extrémités des

éprouvettes cylindriques ont été rectifiées par une rectifieuse
des éprouvettes, selon la norme [EN 12390-2].
Figure 1 : Essais pour la caractérisation du béton testé à l’état frais

ISSN 1737-9334
III. TECHNIQUES EXPERIMENTALES utilisés. À noter que nous avons utilisé un agent de couplage
présenté par un gel de transmission ultrasonique afin
A. Analyse par ondes ultrasonores d’analyser l’échantillon testé.
Les méthodes ultrasoniques sont des techniques efficaces
permettant de contrôle non destructif des caractéristiques
mécaniques des matériaux utilisés [6]. En effet, ces outils B. Analyse par Pulse
sont largement adaptés à une très vaste gamme de matériaux. Ce type d’analyse s’appuie sur l’extraction des fréquences
Pour les bétons, qui sont caractérisés par des milieux propres obtenues expérimentalement par l’analyse vibratoire.
hétérogènes à l’échelle microscopique, l’analyse ultrasonore La méthode de mesure consiste à mettre un échantillon de
permet d’autoriser l’évaluation de ses caractéristiques même type et les mesures dynamiques sont réalisées après
mécaniques avec une échelle de la longueur d’onde. excitation de l’échantillon par marteau d’impact et la réponse
est captée par un capteur type accéléromètre. Le modèle
Ce type d’analyse repose sur un principe physique très analytique de Bernoulli-Euler est exploité pour déduire la
simple qui réside sur le déplacement des ondes traversé dans valeur du module d’élasticité à partir de la valeur
un milieu donné [12]. Ainsi, il est intéressant de noter que les expérimentale de la fréquence propre. Ces valeurs des
modifications des paramètres mesurables associés au passage modules d’élasticité sont utilisées pour calculer
d’une onde ultrasonore sont [13] : le temps de propagation, le numériquement les fréquences propres qui seront confrontées
contenu fréquentiel, l’atténuation et la diffusion de l’onde [5]. aux valeurs expérimentales.
Ces différents paramètres peuvent être directement corrélés
aux changements des propriétés physiques et mécaniques Le principe du banc expérimental utilisé pour réaliser des
comme la raideur, le module élastique, etc... [8,23]. Le essais dynamiques est simple. Il comprend une chaîne
principe de base du contrôle ultrasonore consiste à envoyer d’acquisitions permettant d’enregistrer les signaux
une onde élastique par l’intermédiaire d’un transducteur dans d’excitations et de réponses de l’échantillon étudié. Ces
le matériau à contrôler [7,28]. signaux seront par la suite traités et analysés par une carte
d’acquisition et un logiciel de traitement. Comme le
Dans le cadre de notre étude expérimentale, nous avons présentent les Figures V.2 et V.3, ce dispositif est muni d’un
utilisé un banc d’essais avec un contact direct comme montre excitateur (marteau de choc) et de capteurs de mesure.
la figure 2. Ce dispositif expérimental se décompose d’un
générateur d’impulsions, d’un émetteur et d’un récepteur
piézoélectrique permettant de transformer un signal électrique IV. RESULTATS EXPERIMENTAUX ET ANALYSES
en onde ultrasonore. Les modules d’élasticité déduits par les deux méthodes
(ultrasonore(A) et Pulse(B)) de fréquence de résonance, sont
consignés dans le Tableau 1, pour les variétés des poutres en
béton testées portées en températures de (25°C), (105°C),
(175°C), (250°C), (350°C), (450°C), et (550°C). Les valeurs
des modules d’élasticité et des modules d’élasticité relatifs,
définis comme les rapports des modules d’élasticités sur le
module d’élasticité déterminé à T=25°C, sont regroupés sur
les Figure 3 et 4 respectivement .
Tableau 4
Résultats expérimentaux du module d’élasticité obtenus par les deux
méthodes utilisés
Figure 2 : Schéma descriptif de transmission
en contact direct [13]
Module d’élasticité
Température
Module d’élasticité
relatif
Dans ce travail expérimental, nous avons utilisé un (GPa)
(°C)
excitateur qui se caractérise par une fréquence fondamentale

inférieure à 54 kHz. Ce type des ondes émises permet de se Méthode Méthode Méthode Méthode
propager dans un milieu de couplage avec lequel le matériau “A” “B” “A” “B”
25 34.75 32.28 01.00 01.00
est en contact [24]. En revanche, le récepteur piézoélectrique 105 28.52 27.14 00.82 00.84
permet de recevoir le signal ultrasonore transmis au travers 175 18.81 15.36 00.54 00.47
du matériau testé [4,27]. Le signal formé a été mesuré à l’aide 250 12.51 09.02 00.36 00.27
d’une chaîne d’acquisition du signal de sortie reliant par un 350 03.89 02.50 00.11 00.07
ordinateur. 450 01.75 01.45 00.05 00.04
550 00.42 00.38 00.01 00.01
Le banc à contact utilisé permet d’engendrer des ondes

longitudinales et/ou transversales selon le choix des capteurs

ISSN 1737-9334
porosité qui contribue à l'endommagement du matériau

[19,26].
Dans la suite le module diminue atteignant d’environ 60%

à 99°C. Après une exposition à une température de 350°C, les
enchantions de béton présentent des rigidités très faibles. Cet
affaiblissement peut s’expliquer par l'augmentation du
volume poreux dans le béton et aussi à la fissuration de la
zone interface pâte-granulats [11, 18, 20].
L’analyse comparative des résultats expérimentaux de

module d’élasticité nous a permis de faire les observations
suivantes :
Quelque soit la température de mesure, nous pouvons

Figure 3 : Comparaison (deux méthodes)de l’évolution de module élastique observer que le module d’élasticité obtenu par l’analyse
avec la température pour le béton testé vibratoire, est nettement inférieur que celui déterminé par la
fréquence de résonance. L’origine responsable de cette
augmentation n’est pas assez compréhensible, mais
probablement elle est liée à la gamme de fréquence différente
pour chacune de méthode utilisée, les modalités d’essais
distincts et la réponse de l’échantillon testé avec la méthode
utilisée.
L’erreur enregistrée entre les deux modules d’élasticité,

durant la variation de la température (de 25°C à 550°C), est
acceptable. Cependant, pour la température à 250°C, l’erreur
est de l’ordre de 15% ce qui n’est pas admissible. Ceci est
vraisemblablement dû à l’existence d’erreur de mesure
expérimentale, de modèle analytique utilisé, s’appuyant sur la
simplification des hypothèses qui ne correspondent pas à la
réalité du matériau béton.
Ces résultats expérimentaux montrent que les modules

Figure 4 : Comparaison (deux méthodes) de l’évolution de module élastique d’élasticités, déterminés par les deux méthodes
relatif avec la température pour le béton testé expérimentales sont proches.
A partir de ces Figures, nous pouvons observer que
l’évolution du module d’élasticité varie avec la marge de
température, par laquelle le BAP est exposé. Cette variation V. CONCLUSION
comprend principalement deux phases :
Les résultats obtenus montrent que les modules d’élasticités,
déterminés par les deux méthodes expérimentales sont
Dans la phase (I), entre 25°C et 105°C, nous notons une proches.
diminution des modules d’élasticité d’environs 1-18%. Cette
diminution peut avoir comme origine le départ de l’eau du A partir de cette étude expérimentale, nous pouvons constater
matériau [22,25]. Puis, dans la phase (II) entre 105°C et clairement que les résultats montrent une très grande
550°C, nous observons une diminution progressive du concordance entre les valeurs de modules d’élasticité déduites
module d’élasticité au cours de l’échauffement du matériau, par les deux méthodes. Ils sont presque toujours dans
ce qui traduit par une décroissance de la pente des courbes de l’intervalle de confiance fournis pour les essais
l’écart avec la température. expérimentaux. Cependant, nous observons quelques cas,
caractérisant par une erreur importante (14.65% et 15.02%).
A 175°C, une réduction du module d’environ 50% est Ceci peut avoir diverses origines :
enregistrée. Cette diminution peut être expliquée par
l'apparition de micro fissures sur les poutres testées et bien Causes d’erreur dues aux conditions expérimentales : même
entendu par la déshydratation entrainant une élévation de la si nous avons tenté de travailler de façon rigoureuse, lorsque

ISSN 1737-9334
nous avons fait les essais expérimentaux, certaines causes [11] NEVILLE A.M, «Propriétés des bétons », Centre de Recherche
Interuniversitaire sur le Béton, Eyrolles, 2000, 806 p.
d’erreur ont rendu les données obtenues par notre expérience
moins précises. Les données fréquentielles ont une origine [12] PIRANDA J., « Cours des Dynamiques des structure », ENSMM,
expérimentale, ce qui implique l’existence d’erreur de mesure. Besançon, 1993.
En plus, durant la réalisation des essais par analyse vibratoire [13] PIRANDA J., « Analyse modale expérimentale », Mesures et Contrôle,
Techniques de l’ingénieur, 2001
à l’aide de l’instrumentation PULSE, nous avons procédé
qu’un seul essai pour chaque mesure. [14] SCHNEIDER U. (1988)., « Concrete at high temperatures : A general
review. Fire safety Journal », 1988, vol 13, p 55-68.
Causes d’erreur dues aux effets de structure : les effets de [15] REYNOUARD J-M., PIJAUDIER-CABOT G., «Comportement
mécanique du béton », Mécanique et Ingénierie des Matériaux, Sciences
structure étudiée liés au gradient thermique sont limités, par
germes, 2005, 383 p.
le choix de la vitesse de variation de température, de la taille
et forme de la poutre testés. Cependant, ces effets ne sont pas [16] Rheometric Scientific, “ARES instrument manua”l, 1997.
complètement éliminés en sachant que les observations sont [17] RILEM , «Recommandations RILEM CPC8, Modulus of elasticity of
concrete in compression», Materials and Structures, 1972, 6(30), 507-512.
réalisées en régime transitoire. Les propriétés sous
[18] RILEM TC 50 FMC RECOMMANDATION., « Determination of the
température (les effets différent après exposition à la fracture energy of mortar and concrete by means of three-point bend tests on
température, différente de l’ambiante, selon les traitements notched beams ». Materials and Structures, 1985, Vol. 18, n° 106 p 285-290
thermiques). [19] RILEM TC 129 MHT RECOMMENDATION ., « Test methods for
mechanical properties of concrete at high temperatures. Part 3 : Compressive
strength for service and accident conditions ». Materials and Structures,
Causes d’erreur dues au modèle analytique : le modèle 1995, Vol. 28, p 410-414
analytique est construit à partir des hypothèses [20] RILEM TC 129 MHT RECOMMENDATION., « Test methods for
simplificatrices et correspondent à une idéalisation à la réalité mechanical properties of concrete at high temperatures. Part 1: Introduction.
Saclay (France): CEA », March 1997, 22p.,Draft n°3,Report to comittee TC
physique de matériau béton.
129 MH.
[21] RILEM TC 129 MHT RECOMMENDATION., « Test methods for
Causes d’erreur dues aux niveaux de l’environnement réel : il mechanical properties of concrete at high temperatures. Part 6: Thermal
y a de nombreuses sources d’incertitude dans la prédiction strain for service and accident conditions ». London: Imperial College, May
1997, 14p.,.Draft n°11, Report to comittee TC 129 MHT.
des environnements des systèmes vibratoires, la principale
[22] RILEM TC 129 MHT RECOMMENDATION ., « Test methods for
étant la variation caractéristique des niveaux vibratoires d’un mechanical properties of concrete at high temperatures. Part 7: Transient
point à un autre de la structure étudiée. creep for service and accident conditions ». London: Imperial College, March
1997, 13p,.Draft n°9, Report to comittee TC 129 MHT.
[23] ROCARD Y., «Dynamique générale des vibrations », Masson, 1971,
Références bibliographiques 4eme édition, 460 pages.
[24] ROSEA, M., « Vibrations des systèmes mécaniques : méthodes
[1] Acker, P., J. Torrenti, et F. Ulm , «Comportement du béton au jeune âge»,
analytiques et applications», MASSON éditeur, ISBN 2-225-80205-X, 1984,
Traité Mécanique et ingénierie des matériaux, Hermès, 2004.
488 P.
[2] ACKER, P., TORRENTI, J., ULM, F « Comportement du béton au jeune
âge», LAVOISIER, ISBN 2-7462-0980-2, 2005.
[25] SCHNEIDER U., «Behavior of concrete at high temperatures », Paris
[3] ALARCON-RUIZ, L., GALLE, C., MASSIEU, E., « Analyse de
RILEM, 1982, Report to committee n°44-PHT, 72 P.
l’évolution de la porosité des pâtes de ciment à haute température »,
[SCH 88] SCHNEIDER U., «Concrete at high temperatures », Fire safety
Matériaux 2002, Tours
Journal, 1988, Vol 13, P 55-68.
[4] ALVANDI A., « Contribution à l’utilisation pratique de l’évaluation
dynamique pour la detection d’endommagement dans les ponts ». Thèse de
[26] SULLIVAN, P.J.E., SHARSHAR, R. «The performance of concrete at
doctorat, Ecole nationale des ponts et chaussées Paris, 2003.
elevated temperatures: as measured by the reduction in compressive strength
[5] ANDRE P., «Vibrations aléatoires et analyses spectrale», Presse
» Fire Technology, August 1992, pp 240-250.
polytechniques et universitaires romandes, Lausanne, 1990, 343 p.
[27] THOMAS G., «Dynamique des structures : analyse modale numérique»,
Presse polytechniques et universitaires romandes, Lausanne, 1997, 570
[6] BERTHAUD Y. « Mesure de l’endommagement du béton par une
pages.
méthode ultrasonore ». Thèse de Génie Civil :Université de Paris 6, 1988.
175p.
[28] THYAGARAJAN S. K., SCHULZ M.J., PAI P.F., «Detecting structural
damage using frequency response functions», Journal of Sound and
[7] COMBES, B., « Vibrations des structures pour l’ingénieur et le
Vibration, 210:162–170, 1997.
technicien : Théorie et applications», Ellipses Edition Marketing S.A, ISBN
978-2-7298-2079-1, 2009, 251 P.
[8] GARNIER V., CHAIX J.F., CORNELOUP G., « Caractérisation non

destructive du béton par analyse d’ondes ultrasonores », Bulletin des
laboratoires des Ponts et Chaussées, 239:105–118, 2002.
[9] GMUR, T., « Dynamique des structures : Analyse modale numérique »,
Presses polytechniques et universitaires romandes, ISBN 2-88074-333-8 ,
1997.
[10] NECHAD H., « Evaluation de l’endommagement et de la rupture de

matériaux hétérogènes par ultrasons et émission acoustique : estimation de la
durée de vie restante », Thèse de doctorat, Institut National des Sciences
Appliquées de Lyon, 2004.

ISSN 1737-9334
Amélioration de la performance énergétique d’un

bâtiment résidentiel (cas de la ville de Bechar)
Ibtissame BENOUDJAFER#1, Imane BENOUDJAFER*2,
#1
Department of architecture, University of BECHAR
BECHAR, ALGERIA
tissamo@yahoo.fr
#1
Department of architecture, University of BECHAR
BECHAR, ALGERIA
imenousse@yahoo.fr

ISSN 1737-9334
résumé- La maîtrise de l’énergie constitue un véritable enjeu professionnels du bâtiment à savoir: DTR C 3-2 et DTR C
pour les problèmes environnementaux, économiques et 3-4 qui visent la limitation de la consommation énergétique
sociaux et le secteur résidentiel offre des possibilités relative au chauffage et la climatisation des locaux et DTR
importantes en vue de réduire la consommation énergétique. C-3.31 qui fournit les principes généraux qu’il y a lieu
Nous intéressons plus particulièrement au contexte saharien
d’adopter lors de la conception des installations de
où la consommation d’énergie est presque énergivore où
aucune réglementation thermique ou énergétique n’est ventilation naturelle.
appliquée pour ce secteur. Pour ce faire, nous proposerons de Le premier texte est largement inspiré des règles
faire des enquêtes sur le terrain et des questionnaires réalisés françaises dites règles Th.
au cours de cette étude, nous permettent de mettre en évidence Le second est tout simplement une application de la
les défauts de la construction à Bechar et les besoins des méthode de Carrier pour l’Algérie (Larbi Youcef, 2001).
occupants. Ces enquêtes nous amènent à comprendre la
situation réelle et les problèmes sur terrain. En second lieu, Cependant, l’absence d’appliquer la réglementation
après un diagnostic des cas existants, nous proposons une série thermique aux logements construits, issue une situation
de simulation des systèmes passifs pour un type d’habitation,
d’inconfort des habitations et une augmentation des prix
en s’appuyant pour cette dernière sur les spécificités de
fonctionnement et d’usage du bâtiment et de ses équipements des factures d’électricité et des matériaux de construction.
de chauffage, d’eau chaude sanitaire et électroménagers. Elle Dans ce contexte, et plus largement le domaine de
est formée d’une phase d’analyse sur site et de mesures et construction n’est soumis à aucune obligation en termes de
d’une phase de calcul des consommations et du potentiel confort thermique et de consommation énergétique.
d’économie d’énergie. Nous montrons dans cette étude A l’heure actuelle, Le secteur des bâtiments résidentiels
l’intérêt et les difficultés de la prise en compte la spécificité de en Algérie, a cherché tout développement d’une utilisation
la région à travers la littérature et des simulations rationnelle de l’énergie vis-à-vis des ressources
numériques. énergétiques (combustibles, électricité, etc.).
Mots clés : secteur résidentiel, la maitrise d’énergie, les
Ce secteur constitue un gisement potentiel important
performances énergétiques, La consommation
d’économie d’énergie, notamment par la rationalisation de
d’énergie.
l’utilisation de l’énergie finale, et plus particulièrement de
I. INTRODUCTION l’électricité.
C’est dans ce contexte que s’insère notre étude.
Maintenir à un niveau raisonnable reste un défi majeur
Nous nous appuyons dans ce travail de recherche
dans le domaine de construction, un des secteurs clés pour
d’améliorer la performance énergétique du secteur
réformer le modèle énergétique.
résidentiel
Le gisement d’économie d’énergie dans ce secteur est
mieux l’exploiter, permet de diminuer la facture des
ménagers, de gagner le confort et réduire les émissions en
Pour ce type de bâtiment, les simulations thermiques
gaz à effet de serre.
permettent d’une part d’identifier les voies d’amélioration
et d’autre part de justifier l’opération pour la recherche de
Le secteur résidentiel est particulièrement sensible à la
financement.
problématique de la réduction des consommations
Dans le milieu de la recherche, cette démarche de
énergétiques.
réhabilitation passe actuellement par des techniques
Il est apparu essentiel d’améliorer la performance
d’optimisation afin d’aider le concepteur dans la recherche
énergétique dont le secteur résidentiel.
et le choix de solutions [1].
La consommation dans le secteur de la construction est
comprise entre 40% et 55% de la demande totale dans la
Afin de réduire le nombre de simulations à réaliser, nous
ville du sud en Algérie [16].
avons choisi d’étudier les systèmes passifs pour l’enveloppe.
Dans cette étude, on s’emploiera à améliorer la
Les résultats obtenus seront ensuite rassemblés dans le
performance énergétique pour ce secteur, par des scénarios
but d’évaluer les différentes combinaisons d’actions
d’amélioration qui seront élaborés sur la base de
envisageables. La réalisation de ce dernier point nécessite
programme de la performance énergétique cohérent et
de développer une méthode permettant de calculer la
adaptés aux régions à climat chaud et sec, en générale et de
consommation d’électricité des appareils de climatisation
la ville de Béchar, en particulier.
directement à partir des besoins de refroidissement des
bâtiments.
Il convient de souligner que l’Algérie a pris conscience
de valoriser les ressources pour améliorer le cadre
Dans la pratique, des simulations annuelles au pas de
énergétique des constructions.
temps horaire sont réalisées afin d’obtenir une évaluation
Citons la mise en application de la loi 99.09 relative à la
globale du bâtiment. L'indicateur classique d'estimation des
maîtrise de l’énergie dans le secteur du bâtiment, et un
performances énergétiques est le besoin d’énergie en
décret exécutif n°2000-90 qui a pour objectif l’introduction
kWh/m².an. [1]
de l’efficacité énergétique dans les bâtiments neufs à usage
d’habitation et autres et dans les parties de constructions
Parmi ces systèmes étudiés, l’isolation peut se placer de
réalisées comme extension des bâtiments existants.
diverses manières dans un mur : à l’extérieur, en sandwich
Afin d’amorcer une réglementation, on souligne ici trois
ou à l’intérieur.
documents techniques réglementaires à l’usage des

ISSN 1737-9334
Pour ces trois cas, la valeur des transmittances des thermique, qui ne fait pas de ce matériau, un matériau de
parois restent les mêmes et par conséquent les déperditions choix sur le plan du confort thermique et qui par
à travers l’enveloppe également. conséquent, augmente la consommation électrique pour
Cependant, la position de l’isolant modifie l’inertie de la permettre à l’usager d’obtenir des ambiances de qualité.
paroi : seule la masse située à l’intérieur de l’isolation Cette situation provoquée par la pression du besoin
contribue à une inertie utile pour le volume intérieur. [2] d’habiter et la crise du logement, a révélé de nouvelles
contraintes sur le plan du confort (isolation absente ou
Avant l’apparition ces matériaux isolants, la construction insuffisante, nécessité d’une climatisation onéreuse et
traditionnelle avait privilégié deux types de systèmes inesthétique, etc.).[14]
constructifs : la construction massive en pierre ou en brique
(à forte inertie) et la construction légère en bois (à faible
inertie). II. METHODOLOGIE
Ces pourquoi nous avons choisie l’inertie comme un Le choix d’un outil de simulation thermique dynamique
deuxième système. permet d’évaluer les besoins énergétiques et le niveau de
confort thermique des bâtiments de manière plus détaillée
Les proprieties thermiques des matériaux de construction que les méthodes basées sur des bilans mensuels.
ont une influence considérable sur le comportement Le logiciel TRNSYS(Laboratoire ENERGARID,
énergétique du bâtiment que ce soit en termes d’inertie Université TAHRI Mohamed, Bechar)est un outil de
thermique ou de propriété isolante. simulation thermique dynamique de bâtimentdéveloppé
En effet, les propriétés des matériaux et la façon dont ils plusieurs années. Le logiciel calcule le comportement
sont agencés définissent les caractéristiques qu’aura dedifférentes zones thermiques d’un bâtiment en régime
l’enveloppe du bâtiment. dynamique et il permet également desimuler le
Malgré l’importante avancée que connaît aujourd’hui le comportement de l’enveloppe. Il calcule également les
secteur du bâtiment, il n’existe pas une réelle méthodologie consommationspour le chauffage et laclimatisation. Avec
qui permette de mettre en œuvre une sélection convenable unenouvelle bibliothèque TESS (matériaux, éléments,
de ces matériaux. [2] systèmes, équipements, etc.) permet d’avoir la possibilité
de modéliser tous les éléments souhaités dans l’enveloppe
Une étude a été effectuée pour ces régions dans le but de et son système. Celle-ci comporte souvent :
maintenir un confort par des conceptions architecturales Une description géométrique du bâtiment par saisie de
passives des bâtiments résidentiels. Citons l’étude d’une plans, où l’on spécifie les dimensions et orientations des
maison traditionnelle avec ouverture dans les climats pièces, parois, vitrages, ouvrants.
chauds et secs, réalisée par FIZIOUI N. et BENYAMINE les données de localisation correspondantes (latitude,
M, en 2016. Les résultats sont présentés sous forme longitude et altitude).
d'heures chaudes et froides en un an, maximum et minimum. Une description de l’enveloppe du bâtiment, avec
Les températures de confort intérieures et la fréquence de indication de la composition des parois, des vitrages, des
ventilation efficace par personne comme indicateur de portes... et des propriétés physiques (thermiques et optiques)
performance de qualité de l'air dans chaque pièce. [16] des matériaux qui les constituent.
Une description des équipements (chauffage, ventilation,
Plusieurs études ont été réalisées dans le but d’augmenter eau chaude sanitaire, éclairage...) à travers les paramètres
la performance énergétique des bâtiments résidentiels. qui décrivent leur fonctionnement.
Citons l’étude de rénovation d’une maison à Marseille, dans Une description des sollicitations externes (données
le sud de la France, réalisée par C. FLORY CELINI en météorologiques) et internes (scénarios/données de
2008. ventilation, occupation...). [7]
Cette étude montre une réduction des besoins d’un L’objectif est donc l’amélioration du comportement
facteur. Les systèmes passifs adoptés sont : l’isolation thermique des bâtiments par proposition des solutions
transparente, des fenêtres performantes (double vitrage passives. Dans notre contexte d’étude, ces solutions
basse émissivité), un puits canadien, des matériaux à d’amélioration seront évaluées en termes de potentiel
changement de phase (MCP), des stores opaques et d’économie d’énergie et aussi de performance énergétique.
l’isolation de la toiture.[3] Ceci permettraitd’établir des recommandations sur les
meilleures solutions à privilégier en vue de diminuer les
impacts liés au confort d’été.
Une étude a été réalisée par BENOUDJAFER I.,
présente une étude comparative relative à l’efficacité III. SIMULATIONSET RESULTATS
énergétique de deux appartements situés à Béchar, montre
que la majorité des constructions des habitations fait appel à A. Description de l’enveloppe de référence
une nouvelle conception architecturale, prenant en compte
des matériaux nouveaux inadaptés au climat aride des zones Une description complète de l’enveloppe nécessite une
sahariennes [14]. Il s’agit principalement du matériau béton connaissance des éléments entrant dans la construction pour
(structure poteaux-poutres et blocs de parpaings) connu pouvoir établir un modèle du bâtiment. Ces informations
pour ses multiples avantages (plasticité, facilité de mise en sont résumées en trois catégories et illustrées dans le
œuvre, résistance mécanique importante acquise à court tableau 1 :
terme, etc.). Son inconvénient est lié à sa faible résistance
Copyright - 2019
ISSN 1737-9334 Page 65
types et composition des parois opaques (épaisseurs,

matériaux) avec leurs caractéristiques thermo physiques ;
caractéristiques des menuiseries (vitrages et cadres) ;
mode constructif, pour essayer de déterminer les ponts
thermiques linéaires et l’infiltration d’air.
B. Descriptions de différents systèmes étudiés
Pour maintenir un confort pournotre habitation, il

estindispensable de rechercher des solutions techniques
adoptables et cohérentes avec notreville en matière de
conditions d’application (techniques disponibles et savoir-
faire locaux).L'objectif est de réduirela consommation
énergétique et d'améliorer le confort dans les logements
résidentiels. Parmi ces systèmes passifs, nous pouvons citer
les suivantes :
TABLEAU 3.DESCRIPTION DES SYSTÈMES PASSIFS ÉTUDIÉS MATÉRIAUX,
Inertie :l’emploid’une double paroi, avec une lame d’air ÉPAISSEUR ET COEFFICIENT U-VALUE
permetde déphaser et d’amortir lesapports de chaleur
externe ; cette massevolumique peut jouer un rôle
d’accumulateur au cœur du bâtiment. Epaisse U-
Compositi
Description ur en Value(W
on
(m) /m2K)
Fenêtres performantes : en utilisant un double vitrage. Enduit extérieur 0.02
C’est un élément essentiel, avec un rôle important aux Isolation Enduit de finition 0.02
niveaux énergétique et visuel. des parois Enduit en polyuréthane 0.02
opaques Plaque en polystyrène 0.29
Isolation des parois opaques et toiture : en utilisant des 0.08
extérieure expansé
plaques en polystyrène expansive à l’extérieure. Ce système s Brique rouge 0.20
d’amélioration apporte des avantages pour l’hiver comme Enduit en plâtre 0.02
pourl’été. L’isolation extérieure en hiver peutdiminuer carrelage 0.02
l’échange thermique et les déperditions de chaleur avec mortier 0.02
0.77
l’extérieur. Enrevanche, en été, l’isolation joue un rôle dans sable et gravier 0.05
la diminution des apports de chaleurde l’extérieur par les Isolation de
la toiture béton lourd 0.04
parois. «L’isolation est sans doute le facteur le plus Polystyrène expansé 0.04
important puisqu’elle conditionne votre consommation de Enduit de finition 0.02
chauffage en hiver, voir la climatisation en été » [12] Bloc en hourdis 0.16
mortier 0.01
Peinture blanche 0.03
Isolation du carrelage 0.02
plancher mortier 0.02
bas sable et gravier 0.05
béton lourd 0.25 0.27
mortier 0.01
Inertie Enduit Extérieur 0.015
thermique Brique rouge creuse 0.20
Lame d’air 0.02
Brique rouge creuse 0.10 0.69
Enduit en plâtre 0.02
U-value g-Value T.sol Rf- sol Facte
Fenêtre ur
performant solair
e e
2.95 0.77 0.72 0.12 0.45
C. Résultats des simulations des différents systèmes passifs
III.3.1. Indice d’été :« IPS »

La simulation est effectuée, pendant la semaine la plus
froide (du 7 au 10janvier2016) et la semaine la plus chaude
(du 18 juillet au 24 Juillet 2016).
ISSN 1737-9334
Chambre T ° Réf
performantes
T ° isolation
T ° fenêtres
chaude
T ° inertie
thermique
T° toiture
opaque
TABLEAU 4 : VALEURS DE L’IPS CORRESPONDANT À CHAQUE SYSTÈME
parois
D’AMÉLIORATION (ZONE NORD, SEMAINE 33)
T°max 37.54 31.50 32.78 32.88 30.23
Δ T°Max 6.04 4.76 4.66 7.31
T°Moyenne 32.50 31.85 30.63 31.74 31.57 T°max : température intérieure maximale (C°)
Δ Δ T° max : écart entre température maximale référence et
- 0.65 1.87 0.76 0.93
T°moyenne de système
Tconfort 26 26 26 26 26 T moyenne : moyenne de la température intérieure d’une
semaine
IPS - 0.10 0.28 0.11 0.14 ΔT° moyenne : écart entre moyenne de température
intérieure référence et de système
IPS : indice de confort d’été IPS=(T°ref-T° système)/
Pour le système d’isolation des parois, la température
(T°réf-26)
intérieure est constante toute la journée de 10°C.L’isolation
extérieure augmente la température intérieure dans la
III.3.2. Indices de performances des systèmes « Iéco » et
journée de 0,3 °C, et dans la nuit de 1,3 °C en moyenne, ce
« IPR »
qui traduit la conservation de la chaleur à l’intérieur.
Les différentes consommations en énergie pour produire
Pour l’été, nous pouvons constater une diminution de
le chaud et le froid dans notre site d’étude montrent la
température de 0,6 °C en moyenne pour la protection
nécessité de chercher la rentabilité et la pertinence de
solaire et 0,8 °C pour la ventilation nocturne dans la journée
chacun des systèmes. Un nouveau critère va donc nous
et de 1,6 °C pendant la nuit. Le système d’inertie permet
aider à juger chaque système en termes de rentabilité en
aussi une diminution de température maximales plus faibles
valeur de consommation (chauffage et climatisation).
(0,8 C), et des températures minimales plus élevées (0,7 °C).
De la même manière, les fenêtres performantes permettent
La consommation en chauffage et en climatisation est
de réaliser une diminution de la température maximale de
calculée par TRNSYS en kWh. Nous pouvons donc
l’ordre de 0,3 °C, et une augmentation des températures
exprimer l’indice Iéco dans la relation suivant :
minimales de 0,1 °C.
L’observation de tableau 6 montre que la meilleure
Le tableau 4 présente les valeurs de confort d’été pour
performance est celle du système de ventilation naturelle,
chaque système étudié.
puis de l’isolation parois opaques, toiture isolée, l’inertie
thermique, enfin, la protection horizontale.
Nous constatons une performance négative pour le
système des apports internes, avec un effet marginal des
fenêtres performantes.
L’observation de tableau 6 montre que les systèmes

isolation, inertie, et ventilation nocturne sont économes en
termes de consommation énergétique. Par contre, le
système des fenêtres performantes sont moins économe. Ce
résultat permet de valider les choix en systèmes
d’amélioration vis-à-vis de l’économie en consommations
énergétiques.
TABLEAU 6 : VALEURS DE L’INDICE IÉCO SELON CHAQUE SYSTÈME

D’AMÉLIORATION
Copyright - 2019
ISSN 1737-9334 Page 67
B B B toiture B B
Référen isolation fenêtres inertie
Cep-chauffage : consommation d’énergie en chauffage ce parois perform thermi
(KWh/m2an) opaque antes que
Cep-froide : consommation d’énergie en climatisation Cep- 277.7 101.44 185.01 223.00 112,95
(KWh/m2an) chauffage 1
Iéco : Iéco = (Cepchauffage-Réf – Cepchauffage-Sys)+ Cep- 51.73 47.25 37.25 40.33 48,67
(Cep froide-Réf – Cep froide-Sys) froide
IPR % : Indice de réduction TOTAL 329.4 144.69 222.26 263.33 161,62
4
Iéco - 184.75 107.18 66.11 167.82
L’analyse de ces valeurs est synthétisée de la manière
suivante : Réduction 56,08 32,53 20,067 50,94
%
- le système de fenêtre performante : ce système se
IPR % - +54,86 +32,53 +20,06 +50,9
traduit une hausse en besoin de Chauffage en hiver. 4
- L’inertie thermique : C’est un facteur qui agit Jugement - Retenue Retenue iliminé Reten
principalement par le déphasage et il apporte des bénéfices e ue
en été comme en hiver. Plus l’inertie est élevée, plus les contexte spécifique. Nous pouvons d’apprécier les
améliorations sont donc effectives. meilleures solutions de point de vue énergétique,
-Le système d’isolation des parois verticale et économique et de confort.
horizontale : permet une performance importante en confort
d’été, produit par l’isolation, empêchant la pénétration de la
chaleur à l’intérieur du bâtiment.
Cette amélioration effectuée par application des systèmes

passifs disponibles et de mettre en place des valeurs limites
de ces indices que nous pourrions atteindre dans notre

ISSN 1737-9334
REFERENCES BIBLIOGRAPHIQUES
IV. CONCLUSION
[1] RABOUILLE M., thèse de doctorat, Recherche de la performance en
Dans cette étude, nous avons proposé des solutions simulation thermique dynamique : application à la réhabilitation des
bâtiments. Génie civil. Université de Grenoble, 2014. Français
passives innovantes et disponibles sur le marché local, afin [2] CHAHWANE L., thèse de doctorat Valorisation de l’inertie thermique
d’améliorer la performance énergétique des bâtiments pour la performance ´énergétique des bâtiments. Architecture,
résidentiels de Bechar. aménagement de l’espace. Université de Grenoble, 2011. Français
[3] FLORY-CELINI C., thèse de doctorat en modélisation et
positionnement de solutions bioclimatiques dans le bâtiment
Pour ce faire, nous avons effectué une simulation dans le résidentiel existant, université Lyon 1, 2008, p. 31.
but d’améliorer des indices de performance tel que l’indice de [4] DURAN. C.S., Architecture and energy. Un enjeu pour l’avenir pp 10,
confort d’été IPS et deux indices de la consommation 2011, paris, Edition place des victoires383p
d’énergie IPR et l’Eco. [5] ADRA N., Thèse de doctorat, 2001 Proposition d’une procédure de
certification énergétique des logements et application au contexte
libanais
Le résultat de cette simulation a montré également la [6] IBRAHIM M., Improving the buildings envelopes energy performance
synthèse suivante : using aerogel-based insulating mineral rendering, Mohamad Ibrahim
L’isolation extérieure des parois verticales ou horizontales, 2015
[7] BONTEMPS S. Validation expérimentale de modèles : application
présente une amélioration positive en confort d’été et une aux bâtiments basse consommation. Construction durable. Université
consommation d’énergie très importante. Les pertes de Bordeaux, 2015. Français.
thermiques sont réduites en hiver, et en été, l’isolation permet [8] MOTHANNA S.,thèse de doctorat, Développement d'un label
de diminuer les apports de chaleur par les parois. énergétique destine aux bâtiments résidentiels de la région Est-
Méditerranée (Syrie et Liban). Architecture, aménagement de l'espace.
Le système d’inertie thermique qui agit principalement en Conservatoire national des arts et métiers - CNAM, 2014. Français. .
été, marque un déphasage dans notre climat chaud, apporte [9] SDO, Données annuelles d’Electricité et du Gaz.
des bénéfices en été comme en hiver. Plus l’inertie est élevée, [10] Règlementation thermique des bâtiments d’habitation, Règles de
plus les améliorations dans les trois indices sont donc calcul des déperditions calorifiques Fascicule 1(D.T.R C3-2) 1972
[11] EUGENE D.V., Isolation thermique des constructions en Algérie,
effectives. Alger 1986 p225
Par contre, les fenêtres performantes adonc un impact [12] BRIGITTU V., Choisir une énergie renouvelable adaptée à sa maison,
négatif sur le besoin en rafraîchissement. [13] , EYRELLES, 2006 pp5 p93
[14] DJELLOUL A., Simulation du comportement énergétique des
bâtiments résidentiels au sud algérien, Courrier du Savoir – N°17,
Nous pouvons constater que la meilleure solution et la plus Décembre 2013, pp.113-119
économe en consommation énergétique est réalisée par le [15] BENOUDJAFER I., Etude comparative relative à l’efficacité
système d’isolation, avec une réduction de56, 08 %, et 54% énergétique de deux appartements situés à Béchar, Algérie, , Revue
pour l’inertie thermique. Par contre, le système des fenêtres des Energies Renouvelables Vol. 15 N°1 (2012) 1 - 12
[16] RIVALIN L., thèse de doctorat, Vers une démarche de garantie des
performantes ne produit pas de changement important au consommations énergétiques dans les bâtiments neufs : méthodes
niveau de la consommation énergétique globale. d'évaluation des incertitudes associée à la simulation thermique
dynamiquedans le processus de conception et de réalisation. PSL
Il nous semble dès lors nécessaire de revenir aux systèmes ResearchUniversity, 2016, Français.
[17] GRIGNON-M. L., thèse de doctorat Développement d’une
passifs innovants, disponibles et efficaces énergétiquement méthodologie d’analyse coût-bénéfice en vued’évaluer le potentiel de
afin d’améliorer la performance énergétique des bâtiments réduction des impacts environnementaux liés auconfort d’été : cas des
résidentiels existants ou neufs. climatiseurs individuels fixes en France métropolitaine.
[18] FEZIOUI N., et all : The traditional house with horizontal opening: a
trend towards zero energy house in the hot, dry climates, Science
[19]

ISSN 1737-9334
Etude comparative entre les différentes méthodes

extrapolations verticales de la vitesse du vent de la
loi à coefficient de puissance avec le langage
PYTHON
Hadjij Merouane#1, BenfekirAbderrahim*2, ZoubidaLounis#3
LABORATOIRE D'INGENIERIE EN SECURITE INDUSTRIELLE ET DEVELOPPEMENT DURABLE "LISIDD"
UNIVERSITE D'ORAN-2 MOHAMED BEN AHMED
BP 1524 ORAN, Algérie
1
hadjij.marwan@gmail.com
3
lounis_amira @yahoo.fr
LABORATOIRE D'INGENIERIE EN SECURITE INDUSTRIELLE ET DEVELOPPEMENT DURABLE "LISIDD"
UNIVERSITE D'ORAN-2 MOHAMED BEN AHMED
BP 1524 ORAN, Algérie
2
benfkirrahim@yahoo.fr
Abstract—This paper concerns the study of danger about the prévention. En effet, une phase de programmation est
installations of the sites of renewable energies, more exactly wind nécessaire afin de calculer les risques et garantir la sécurité du
farms after we had shown the dangers which could alter a wind produit, cela en se basant sur les équations des deux
turbine and minimize the most commons accidents in the chercheurs Betz et Mikhaiel qui ont établi la relation directe
worldwide, we noticed that the high speed of the wind always
entre la puissance produite, la vitesse du vent, l'aire balayée
returns as main or secondary cause in the incident.
It is in this context that a program with the PYTHON software par le rotor et l'altitude du mât de l'éolienne : tout cela pris en
was established to warn(prevent) the accidents for such wind considération lors du choix du produit. L’ensemble de ces
turbine in such region, in the same time it will help us to compare éléments permet d’identifier la meilleure éolienne pour un site
between different methods so this allows us to view a feasibility précis, une base des données météorologiques étant
and to choice the adequate territory or an adequate wind turbine primordiale pour effectuer les calculs.
in order to preserve the durability of the device by getting it a
performance continues in complete safety . II. LES ACCIDENTS LES PLUS REDOUTES POUR UNE EOLIENNE
L’éolienne, en tant que source d'énergie propre et
Keywords— Wind turbine, Program, Wind, Danger, Risk. renouvelable dépendant principalement de la vitesse du vent,
peut être exposée à des dangers multiples, causés
I. INTRODUCTION essentiellement par le phénomène de survitesse de vent qui
L'énergie éolienne a connu une avancée considérable ces pousse la rotation des pales au bout de leurs limites de
dernières années de telle sorte que ce domaine est investi en conception (Vitesse de déclenchement)
force par plusieurs pays importants tels la Chine, l'Inde, les
USA, le Brésil et les pays européens [1]. • Effondrement d'éolienne
• chute de pales
Cependant, une étude de risque est indispensable afin de • projection de debris
déterminer le site adéquat pour la mise en place des parcs
éoliens sur un territoire donné. Toutefois, le souci premier de III. LA FORMULE D'ALBERT BETZ
tout client ou investisseur dans l'énergie renouvelable réside Albert Betz a réussi à définir la puissance que peut extraire
dans le choix de l'éolienne la plus sûre et la plus performante. l'éolienne de l'énergie cinétique du vent et en a conclu qu’une
partie seulement de la puissance du vent est utilisée par la
L'étude présentée dans ce travail expose les dangers que
peut receler l'installation d’un parc éolienne en cas d'accident turbine d'où le coefficient de Betz de sorte que :
et permet aux clients de bien choisir le type d'éolienne adéquat
pour la région d'installation du projet. Le contenu de cette
étude est en relation avec l’importance des dangers de = . [Watt]
l'installation et des conséquences prévisibles en cas de sinistre.
Le but recherché dans notre étude est la détermination et la
S: Surface balayée par le rotor (m2)

ISSN 1737-9334
ρ: La masse volumique de l'air (1.25 kg/m3) IV. METHODE D’EXTRAPOLATION

V: Vitesse du vent (m/s)
A. Mikhail et Justus [3]

Coefficient de Betz

= 0.5926 59.26% La mesure de la vitesse du vent est effectuée généralement à
des hauteurs manométriques égales à 10 mètres du sol.
la surface du rotor égale à = Toutefois, il serai plus intéressant de pouvoir les produire à
donc on peut écrire l'équation de Betz en fonction du cube de des altitudes dignes d’intérêt telles que les hauteurs des
la vitesse du vent et du carré du diamètre éoliennes. Plusieurs auteurs ont proposé des formules
empiriques permettant l’extrapolation verticale de la vitesse du
vent. La plus connue est la loi de la puissance de Mikhaiel et
16
= Justus. Ayant montré que pour des hauteurs inférieures à 100
27 2 4 m, l’effet de rugosité du sol ne peut être négligé lors de
l’extrapolation de la vitesse du vent, Mikhaiel et Justus. [7] ont
développé une seconde formule empirique, nommée la loi de
dans le cas de ρ=1.25 kg/m3 puissance modifiée.
l'équation de Betz s'écrit Soit une vitesse V1 extrapolée d’une altitude Z1 vers une
altitude Z2, suivant la formule suivante :
= 0.37. .
α1
Z 
A. La puissance électrique récupérée V2 = V1  2 
Outre le coefficient de Betz qui fait diminuer la puissance  Z1 
extraite en raison des frictions des pales et toute résistance
mécanique, appelé coefficient de puissance[2] L’exposant α1 est sous la forme :
= "#. $
 
d'où: 1  0.0881   V1 
α1 = −  ln
=
1
. Z  Z1   6 
$
2 ln 1 − 0.0881*ln
Z0  10 
la puissance électrique de sortie vers le réseau ou au avec:

consommateur est diminuée par le rendement de la boite de Z = exp  ln ( Z1 ) + ln ( Z 2 )  / 2
transmission ƞt et par le rendement du générateur électrique ƞg
sachant que le rendement combiné pour la transmission et le
générateur est égal à 70% de la limite de Betz pour une hélice Z1: 10 mètres
à trois pâles ce qui justifie notre choix Z0: est la rugosité du sol
Alors l'équation pourrait être développée de cette façon : V1:Vitesse du vent à altitude de 10 mètres
V2: Vitesse du vent à altitude de Z2
1 Cela doit être pris en considération du fait que le calcul de

% = "#. ƞ' . ƞ% . . puissance est très important comme la nécessité de connaitre
2
plus exactement la vitesse du vent à laquelle les pales des
éoliennes sont exposées
Alors
B. Modèle de Spera et Richards
()* = % = 0.203. .
D’après l’étude de Spera et Richards [4], le coefficient peut
être calculé en fonction de la rugosité du terrain et de la vitesse
de référence comme suit :
0.2
Z 
α= 0  1 − 0, 55.log (V1 ) 
 Z1 

ISSN 1737-9334
Ce modèle a été dérivé d’un certain nombre d’observations V. CONCEPTION DU PROGRAMME DE SECURITE
faites à travers plusieurs endroits des Etats Unis. De plus il fut
utilisé au centre de recherche Lewis de la Nasa [5] pour Pour concevoir notre programme dans le but d'améliorer la
concevoir des éoliennes à grande échelle. sécurité des éoliennes nous devons prendre en considération
C. Modèle de Mikhail [8] les paramètres de ces dernières tels que le diamètre du rotor, la
longueur du mât, sa vitesse de déclenchement "vitesse cut-
out" : ces données figurent dans la fiche de spécification de
En 1985 Mikhail a repris le modèle de modifié et a proposé toute éolienne fabriquée.
un modèle qui prend en considération des conditions de Il est également indispensable de disposer d’une base de
stabilité moyenne. Le modèle est basé sur le modèle de la loi données météorologiques de la région concernée par
de puissance et le modèle de similitude où le coefficient l'installation (Vitesse maximale journalière pendant 5 ans) et
modifié est donnée par de connaître le coefficient de rugosité du sol de la région visé
α m = α m + b.ln (V1 )
2
avec: A. Objectif du programme
Le programme consiste à étudier chaque éolienne proposée

  par le client dans les circonstances météorologiques de la
1  0.088  région précise en procédant à l’élimination de celles qui ne
αm = +  peuvent supporter la vitesse du vent, les types d’éoliennes
 Zg    Z1   retenues seront classées selon leurs performances ou
ln    1 − 0.088.ln    économiquement en fonction des choix du client
 Z0    10  
Cela va aussi nous aider a comparé entre les quatre
L’expression du coefficient modifié devient ainsi : méthodes d'extrapolation cités ci dessus qui sera la meilleure
pour une bonne prédiction
 
1  0.088 − 0.088.ln (V1 ) 
αm = +  B. Programmation avec le langage PYTHON
2
 Zg    Z1  
ln    1 − 0.088.ln    PYTHON est un langage de programmation à base de C++
 Z0    10   très utilisé à travers le monde en raison de sa compatibilité
avec la majeure partie du web. Cette caractéristique a conduit
au choix de ce langage pour sa faculté à obtenir les
D. Loi de puissance dite à coefficient variable informations du web grâce à des sous programmes spécifiques,
Pour tenir compte des différentes valeurs de la rugosité du la nécessité de disposer de données réelles dans notre
terrain, Knidiri et Laaouina [6] lors de l’établissement de programme étant primordiale afin de mettre en épreuve sa
l’atlas éolien du Maroc en 1986, ont proposé d’écrire crédibilité prévisionnelle.
l’expression sous la forme suivante :
Après la conception du programme avec le PYTHON, une
x − 0.088.ln (V1 ) simulation du programme doit être faite dans le but d'évaluer
α= son efficacité à prédire les accidents. Pour cela nous avons
Z 
1 − 0.088.ln  1  collecté 37 accidents différents enregistrés dans le monde dans
 10  les dernières années et recueilli les données relatives au lieu, la
date et la spécification de l'éolienne, tout en alternant entre les
4 modèles d’extrapolation.
Où le coefficient x varie en fonction de la classe de rugosité
comme suit :
Considérant que le choix de la méthode d'extrapolation
vertical de la vitesse du vent est très important pour le
• x= 0.25 pour 0.005>=Z0>0m
programme, nous allons procéder à une comparaison entre les
• x=0.31 pour 0.05>=Z0>0.005m
quatre modèles cités précédemment ainsi que mettre en
• x=0.37 pour 0.5>=Z0>0.05m
évidence les limites de chaque modèle au cours de cette
• x=0.48 pour 4>=Z0>0.5m
simulation

ISSN 1737-9334
VI. RESULTATS ET INTERPRETA

ATIONS La figure 1 englobe les rrésultats de la simulation du
programme pour différents m modèles d’extrapolations, la
A. Pour le modèle Spera et Richards
ds méthode de Spera& Richards a le meilleur taux de prédiction
suivie par le modèle de Kindiri&
&Laaouinia
La simulation du programme avec le modèle Spera et
VII. ETUDE ET ANALYSE DEES MODELES D'EXTRAPOLATION
Richards a donné comme résultats 26 accideents prédits sur 37
ce qui correspond à un taux de 70.27% .
Pour mieux comprendre et annalyser les quatre modèles tous
La limite concernant ce modèle c'est quee la vitesse du vent basés sur la méthode de la loi dee puissance et diffèrents dans le
n'excède pas les 65.8 m/s (237km/h) assez large pour couvrir calculde l'exposant "alpha", noous avons établi une étude de
les 37 accidents de sorte qu'en aucun cas le programme n’a profil "alpha" en fixant certainss paramètres et variant d'autres
donné des valeurs illogiques telles que la vittesse à une hauteur de manière à exposer les limites de chaque modèle.
élevée est inférieure à celle de la vitesse rééférence ce qui est
impossible physiquement "la vitesse du ventt augmente avec la A. Etude de profil d'aalpha en fonction de la vitesse
hauteur" [9]. réféérence (météo)
B. Pour le modèle Mikhail et Justuss
Lors de la simulation, ce modèle a montréé que son intervalle

de travail est trop étroit : ses limitesétant de
13.6<V1<44.9m/sec
Le programme a donné comme résulltats 24 accidents
prédits sur 37, soit un taux de réussite de 644.86%
C. Pour le modèle Mikhail
Plus précis mais il intervient dans un inntervalle étroit de

l'ordre de 6.2<V1<20.4m/sec.
Ne couvrant pas la majorité des accidentss, le résultat donne
une prédiction de 45.94 % du cas réel
D. Pour le modèle de la loi puissancep dite à Fig. 2 Variation d'alpha en fonctioon de la vitesse météo(3 modèles)
coefficient variable
Le modèle de Kindiri&Laaouiniadonnne d’assez bons Ce graphe montre la variatioon d'alpha des trois modèles en
résultats lors de la simulation du programmee de sorte quela fonction de la vitesse météo (m m/sec) et met en évidence les
prédiction avec ce modèle frôle les 67.56% % comme taux de limites de cohérence de chaque m modèle.
réussite Il est constaté que le modèle de Spera& Richards couvre un
assez large intervalle de bon foonctionnement par rapport aux
deux autres modèles "Mikhail"" et "Mikhail& Justus", avant
d'atteindre sa limite, sachant que
q les paramètres hauteur de
l'éolienne et rugosité ont été fixxés respectivement à He=80m ,
Z0=0.05m
Pour le modèle de Kindiri& &Laaouinia qui dépend de la

variable x selon l’intervalle de la rugosité, un graphe lui est
dédié
Fig. 1 Histogramme des taux de prédiction pour différent

d modèles

ISSN 1737-9334
B. Etude de profile d'alp

lpha en fonction de la rugosité
Fig. 3 Variation d'alpha de Kindiri&Laaouinia en foncttion de la vitesse météo
la figure(3) comprend quatre courbes chacunne représentant un

certain intervalle de rugosité : ainsi on peut dire
d que l'équation Fig. 4 Variation d'alpha een fonction de la rugosité
de l’alpha établie par Kinidiri et Laaouinia enn 1986suit l'une (3 moddèles)
des quatre courbes selon la valeur de la rugossité: Le graphe (Figure 4) monntre l'influence du paramètre
rugosité sur l'exposant "alpha" : ce dernier joue un rôle très
• x= 0.25 pour 0.005>=Z0>0m important dans l'extrapolation verticale
v
• x=0.31 pour 0.05>=Z0>0.005m
• x=0.37 pour 0.5>=Z0>0.05m On constate que le modèle "S Spera & Richards" est presque
• x=0.48 pour 4>=Z0>0.5m indépendant de ce paramètre. Par contre les deux modèles
"Mikhaiel" et "Mikhaiel & Juustus" sont influencés par la
Interprétation rugosité de manière similaire sau
auf pour les valeurs très faibles
On distingue dans le graphe (figure3) que
qu pour la courbe de ce paramètre
x=0.25 la limite est de 17m/sec , pour x=00.31 la courbe est
limitée à 44m/sec, celle de x=0.37 n'excèdde pas les 67m/sec
enfin la dernière peut aller jusqu'a 240 m/ssec mais c'est pour VIII. CONCLUSION
O
une rugosité très élevée ( la hauteur de l'éoolienne est fixée à
80m). Ce programme nous a permiss de choisir l'éolienne la mieux
indiquée pour une région donnéée afin d'avoir la certitude que
Pour le graphe (figure 2) nous constatonss que le modèle de celle-ci ne courra aucun danger, opérera en continu et pourra
Spera& Richards couvre une assez largge étendue avant produire de l'énergie tout au longg de sa durée de vie.
d'atteindre la limite logique d'ordre 65.8mm/sec suivi par le
modèle de Mikhaiel& Justus avec une plagee médiane limitée à Ainsi on déduit que le modèèle "Spera & Richards" à une
45m/sec. En dernier, le modèle de Mikhaieel qui n'excède pas assez bonne prédiction dans le programme
p
20m/sec
Le programme sera amélioréé (Python) afin qu'il puisse être
à la portée du client, facile à utiliser et affichant plus
clairement les résultats ce qui le rendra accessible pour tout
investisseur
IX. PERSPPECTIVE
Nous visons à améliorer le prrogramme de sorte à augmenter
son taux de prédiction entre auttre nous nous engageons aussi
pour modéliser une méthode d''extrapolation qui possède une
bande de travail assez large
Page 74
Copyright - 2019
ISSN 1737-9334
REFERENCES
[1] The World Wind Energy Association, «2014 Half-year
Report»,WWEA Head Office Charles-de-Gaulle- Str. 5 53113
Bonn,2014 Germany.
[2] Albert BETZ "Wind energie und

ihreAusnutzugdurchWindmülen," "Wind Energy and it's
extraction through Windmills, "Germany,1926.
[3] Mikhail et al., ‘Comparaison of Height Extrapolation Models

and Sensitivity Analysis’, Wind Engineering, 1981,Vol. 5, N°2.
[4] Spera DA, Richards TR. Modified power law equations for
vertical wind profiles. In: Conference and workshop on wind
energy characteristics and wind energy siting, Portland, OR,
USA; June 1979.
[5] Spera DA. Wind turbine technology: fundamental concepts of

wind turbine engineering. New York: ASME Press; 1994.
[6] Knidiri F. and Laaouina A. (1986) L'énergie éolienne au Maroc,

CDER.
[7] Mikhail, A. S., and Justus, C. G.,"Comparison of Height

Extrapolation Models and Sensitivity Analysis," Wind
Engineering, Vol. 5,No. 2, 1981
[8] Mikhail, A. S., Height Extrapolation of Wind Data, Journal of

Solar Energy Engineering, 107(1985)10-14
[9] N. KasbadjiMerzouk, Evaluation du gisement énergétique

éolien. Contribution à la détermination du profil vertical de la
vitesse du vent en Algérie, Thèse de doctorat, Université Abou
BekrBelkaid Tlemcen, Unité de Recherche Matériaux et
Energies Renouvelables URMER Tlemcen, 2006, pp.13

ISSN 1737-9334
Caractérisation physico-chimique et mécanique des

argiles de la région d’Adrar et de Timimoun en vue
d’une utilisation dans les centres de stockages des
déchets
BELAIDI Khedidja #1, MEKERTA Belkacem *2, ALLAL Mohamed Amine #3
#1
Enseignant Chercheur, Département science de la technologie, Université Ahmed Draia Route N6-01000, Adrar, Algérie
#3
Professeur Chercheur, Département de génie civil, Université Aboubekr Belkaid Faculté de Technologie,B.P N°230 Pole
Chetouane-13000,Tlemcen, Algérie
khadidja.univ@yahoo.fr
aminallal@yahoo.fr
*2
Professeur Chercheur, , Université Ahmed Draia Route N6-01000, Adrar, Algérie
mekertab@yahoo.fr
Résumé— Nos travaux de recherche ont pour but La présente étude fixe son objectif principal dans la
d’étudier expérimentalement les argiles compactées connaissance des différentes caractéristiques physiques,
saturées utilisées dans la construction de la barrière chimiques et mécaniques des argiles. De ces matériaux
passive du Centre de Stockage des Déchets (CSD). Le particuliers, on a réalisé des essais œdométriques qui ont
choix du matériau utilisé pour la construction de la permis de trouver les paramètres de gonflement et le
barrière passive est une argile qui provient de deux coefficient de perméabilité K, qui est un paramètre
gisements différents, à savoir celui d’Adrar et de fondamental au niveau du choix des CSD dans
Timimoun. Dans la première partie, on s’intéresse à la l’enfouissement des différents déchets (inertes ou dangereux).
connaissance des différentes caractéristiques physiques,
II. QU’EST CE QU’UN CSD
chimiques et mécaniques de ces argiles. La seconde étape
de cette étude, nous a conduit à traiter les essais Un centre de stockage des déchets comprend une surface
œdométriques qui permettront de trouver les paramètres de plusieurs dizaines d’hectares. Il est composé de casiers
de gonflement (taux et pression de gonflement). On a indépendants sur le plan hydraulique dans lesquelles sont
poursuivi nos travaux de recherche dans la détermination confinés les déchets [1] (voir Figure 1).
du coefficient de perméabilité par les deux méthodes de
Casagrande et Taylor.
Mots-clés — Argiles d’Adrar, centre de stockage des déchets,

barrière passive, paramètres de gonflement, perméabilité.
I. INTRODUCTION
La croissance de la population dans les zones urbaines et
les centres ruraux conduisent à la naissance de volumes de
déchets de plus en plus importants. Le confinement des
déchets dans les centres de stockage des déchets (C.S.D) est
la méthode la plus couramment utilisée actuellement. Afin de
préserver les eaux souterraines, les décharges sont construites
sur une couche d’argile naturelle ou reconstituée de faible
perméabilité, connue sous le nom de barrière de sécurité
passive. Pour les déchets ménagères et assimilés, la plupart
des réglementations en vigueur exigent une perméabilité Fig. 1 Schéma type d’un CSD en Algérie [1]
maximum de 1x10-9 m/s.

ISSN 1737-9334
La barrière de sécurité passive est destinée à garantir D’après Skempton, l’activité d’une argile (Ac) est
l’étanchéité des casiers en cas de défaillance de la barrière de définit selon l’équation (1) :
sécurité active sous laquelle elle se trouve (voir Figure 2). La Ac= Ip (%) / % des éléments < 2 µm (1)
valeur de la perméabilité tourne autour de 10 -9 m/s pour Avec Ip : l’indice de plasticité
1mètre d’épaisseur de la barrière passive [2]. Le pourcentage de l’argile considéré est obtenu à partir
des analyses granulométriques données dans le tableau 1.
Suivant les résultats obtenus, on peut conclure que l’argile
d’Adrar à une activité normale, par contre celui de Timimoun
est peu actif.
TABLEAU I
RECAPITULATIF DES RESULTATS DES ESSAIS PHYSICO-
CHIMIQUES ET MECANIQUES SUR LES DEUX TYPES D’ARGILE
Propriétés Argile Argile

des Type d’argile Adrar Timimoun
argiles
Pourcentage des elements < 46 36
2 µm (%)
Analyses Limite de liquidité (%) 79 36
Fig. 2 Principe de l’étanchéité du fond de CSD selon la réglementation Physiques Limite de plasticité (%) 36,74 11,4
française, [2] Indice de plasticité (%) 42,26 24,6
Limite de retrait (%) 12,9 11,8
III. UTILISATION DES ARGILES LOCAUX DE LA Le coefficient d’activité Ac 0,92 0,68
WILAYA D’ADRAR VBS 8,2 5,33
Analyses SST (m²/g) 171,62 111,55
Chimiques CaCO3 (%) 3,59 12,08
Dans nos travaux de recherche, nous allons caractériser
deux types d’argile se trouvant à Adrar dans un rayon de dmax (kN/m3) 17,7 19,9
200Km. Le choix de ces gisements d’argiles est justifié par la Analyses Teneur en eau optimale
proximité des CSD déjà construits (voir Figure 3), à savoir : Mécanique wopt (%) 20,86 10,53
Le premier est proche du CSD déjà réalisé, il est utilisé Indice des vides e0 0,637 0,372
pour les déchets ménagères (Timimoun), et le deuxième Pression de gonflement
gisement sert à alimenter l’usine de briqueterie, il se trouve à Pg (bars) 12,50 7,64
Adrar Nord. Gonflement libre G (%) 39,05 20,10
L’analyse des résultats chimiques nous a conduit à

constater et confirmer que suivant la valeur de bleu, on est en
de sols argileux. D’après les résultats obtenus lors des essais
de la teneur en carbonate de calcium CaCO3 [3], nous
remarquons que l’argile d’Adrar est un sol non calcaire et
l’argile deTimimoun appartient à la classe des sols faiblement
calcaire.
Les résultats du Proctor modifié montrent que l’argile de
Timimoun a une valeur de dmax plus élevée que l’argile
d’Adrar, contrairement à la teneur en eau optimal. L’argile
d’Adrar demande une quantité plus élevée d’eau à cause du
pourcentage des particules fines (% éléments < 2 µm = 46%).
Localisation du site d’Adrar Localisation du site de Timimoun
Fig. 3 Localisation des gisements étudiés IV. ESSAIS OEDOMETRIQUES
A. Caractérisation des sols de la région d’Adrar A. Mesure du potentiel de gonflement

Le tableau 1 résume l’ensemble des résultats des Les argiles gonflantes ont des caractéristiques physiques
essais physico-chimiques et mécaniques réalisés sur les typiques permettant de connaitre leurs potentiels de
argiles étudiés. D’après l’abaque de plasticité de Casagrande gonflement [4]. Le but de cette analyse est de déterminer les
et selon la classification L.C.P.C des sols fins, on a constaté paramètres de gonflement (gonflement libre et pression de
que l’argile d’Adrar appartient à la catégorie des sols très gonflement) des deux argiles sur des échantillons obtenus par
plastiques, par contre l’argile de Timimoun est de consistance compactage statique aux caractéristiques de l’optimum
peu plastique. Proctor. Les essais ont été réalisés sur un oedomètre classique,
l’échantillon est placé dans une cellule cylindrique entre deux
pierres poreuses. Puis suite à l’inondation, on laisse le

ISSN 1737-9334
gonflement se réaliser verticalement sous la pression du méthode de Taylor : Cv= (0,848 x H²)/(4xt90 ) (3)
piston durant plusieurs jours jusqu’à la stabilisation [5].
L’argile d’Adrar a demandé 9 jours de stabilisation, par méthode Casagrande : Cv = (0,197x H²) / t50 (4)
contre l’argile de Timimoun a été stabilisé au 6 ème jour. La Les résultats des valeurs des coefficients de consolidation
variation relative du volume de l’échantillon est notée G (%). des deux argiles sont regroupés dans le tableau 2 ci-dessous.
D’après le tableau 1, la pression de gonflement de l’argile Les valeurs des coefficients de consolidation Cv obtenues par
d’Adrar à une valeur importante, elle est de l’ordre de 12.5 la méthode de Casagrande sont supérieures à celles obtenus
bars, par contre l’argile de Timimoun présente une pression par la méthode de Taylor. On remarque aussi que l’argile de
de gonflement moyennement faible. Timimoun donne une valeur de Cv plus élevée comparée à
La pression de gonflement trouvée est celle qui ramène le l’argile d’Adrar, cela explique que l’argile d’Adrar est une
sol à son volume initial après son gonflement. La valeur de la argile proche des sols très fins.
charge finale lorsque l’échantillon est stabilisé est la pression
TABLEAU II
de gonflement [4]. En ce qui concerne le gonflement libre, la RESULTATS DES COEFFICIENTS DE CONSOLIDATION ET DE
figure 4 montre l’évolution du gonflement libre en fonction PERMEABILITE SUR LES DEUX TYPES D’ARGILES
du temps pour les deux sols. On constate qu’au début de
Type de Type d’argile Argile Argile
l’inondation, les sols gonflent d’une manière similaire en
méthode Adrar Timimoun
fonction du temps, mais par la suite l’argile d’Adrar présente Méthode coefficient de
un taux de gonflement plus important de l’ordre 39.05%. de Taylor consolidation Cv(cm2/s) 0,085.10-3 0,444.10-3
perméabilité k (m/s) 2,03.10-10 5,27.10-10
Méthode coefficient de
de consolidation Cv(cm2/s) 0,171.10-3 0,501. 10-3
Casagrande perméabilité k (m/s) 4,10.10-10 5,94.10-9
D. Variation de la conductivité hydraulique K

On remarque que la conductivité hydraulique de l’argile
d’Adrar déterminée par la méthode de Casagrande est
inférieure à celle de l’argile de Timimoun. Par contre, les
valeurs de la perméabilité des deux argiles obtenues par la
méthode de Taylor sont presque proches. Comme première
proposition dans notre cas, l’argile à prendre en considération
pour la suite de nos travaux de recherche est celle d’Adrar.
Fig. 4 Evolution du gonflement libre en fonction du temps des deux argiles Cependant, on peut poursuivre nos travaux de recherche par
la diminution du gonflement par ajout de la chaux à différents
B. Mesure indirecte de la perméabilité pourcentages et de mesurer par la suite l’évolution de la
L’essai œdométrique permet d’évaluer le tassement et le conductivité hydraulique.
gonflement du matériau au cours du temps, il permet
également de classer les sols vis-à-vis de leur compressibilité. V. CONCLUSIONS ET PERSPECTIVES
On peut estimer la valeur de la perméabilité ksat en Dans cette étude, nous avons réalisé un travail
fonction des paramètres de compressibilité et de expérimental permettant la caractérisation physico-chimique
consolidation selon l’équation suivante (2) : et mécanique des matériaux argileux nommés « argiles
d’Adrar et de Timimoun ». Cette étude a pour objectif
ksat= Cv.γw / E’ (2)
principal l’évaluation de la possibilité d’exploitation d’un
avec : Cv : coefficient de consolidation matériau argileux comme barrière de sécurité passive pour les
E’ : module oedométrique C.S.D, en vérifiant l’argilosité et la perméabilité.
γw : poids volumique de l’eau
Le coefficient de consolidation Cv se calcule en utilisant Le premier critère qui est la granulométrie montre que
les méthodes graphiques de Casagrande et Taylor. On l’argile d’Adrar est la plus favourable, puisqu’elle a une
remarque que pour les mêmes données, ces méthodes teneur en argile de 46%, ce qui assure une bonne
graphiques donnent des valeurs différentes de Cv [6]. perméabilité. L’analyse des résultats chimiques nous a
confirmé que suivant la valeur de bleu, nos sols sont fins et
C. Mesure du coefficient de consolidation Cv que l’argile d’Adrar est un sol très argileux, par contre
Le coefficient de consolidation Cv est un paramètre très l’argile de Timimoun est un sol argileux. D’après les
important pour l’estimation de la conductivité hydraulique, il résultats obtenus lors des essais de la teneur en carbonate de
est fonction de la hauteur de l’échantillon H en cours d’essai calcium CaCO3, nous constatons que l’argile d’Adrar est un
et du temps de consolidation t90 (méthode de Taylor) ou t50 sol non calcaire et l’argile de Timimoun appartient à la classe
(méthode de Casagrande), selon les équations (3) et (4) : des sols faiblement calcaire.

ISSN 1737-9334
La capacité de gonflement joue un rôle important dans coefficient de perméabilité, après traitement aux mélanges à
l’étanchéité de l’argile. Suivant aussi les résultats des limites la chaux 2, 4 et 6%, dans le but du stockage des déchets
d’Atterberg, on peut conclure que l’argile d’Adrar est très dangereux et toxiques au niveau de la wilya d’Adrar
gonflante (référence). (Direction de l’environnement de la wilaya d’Adrar).
Pour le gonflement libre, on constate qu’au début de REFERENCES

l’inondation de l’argile, les sols gonflent d’une manière [1] MC. Rahmani, Durabilité des fondations d’un centre de stockage de
similaire. On remarque que l’argile d’Adrar présente un taux déchets vis-à-vis des sollicitations de service. Mémoire de magister,
Ecole Nationale Polytechnique d’Oran, 111 pages, 2014.
de gonflement important estimé à 39.05%. La pression de [2] CARTAUD F., 2004. Modélisation des écoulements dans les
gonflement de l’argile d’Adrar est élevée, elle est de l’ordre interfaces des barrières d'étanchéité composites d'installations de
de 12.5 bars, par contre l’argile de Timimoun présente une stockage de déchets. Thèse Doctorat de l'Ecole Nationale Supérieure
pression de gonflement moyennement faible (7,64 bars). des Mines de Paris. France.
[3] (2018)Web-1:http://www.lano.asso.fr/web/potentiel_dhydrogene.html
La conductivité hydraulique a été calculée indirectement Calcaire total et calcaire actif. LANO - Laboratoire Agronomique
par les essais œdomètriques par les methodes de Casagrande NOrmand - pH (Potentiel hydrogène).
et Taylor. Les résultats montrent que l’argile d’Adrar est [4] MK. Gueddouda, I. Goual, M. Lamara. et MS. Goual, Comportement
mécanique d’un sol expansif stabilisé par ajout de sable de dune.
relativement imperméable par rapport à l’argile de Timimoun. Revue des Sciences et Sciences de l’Ingénieur, ISSN2170-073. RSSI,
Les deux sols argileux pourront être utilisés comme barriere Vol.02, No.01, p16-26, Décembre 2012.
passive pour le stockage des déchets ménagers et spéciaux en [5] Z. Amraoui, Etude de l’infiltration dans les sols fins non saturé. Thèse
moyennant quelques precautions supplémentaries liées à la de Doctorat de l’INPL , 325 p, 1996 .
[6] M K. Gueddouda, Comportement hydromécanique des sols compactés:
disposition des différentes couches d’un CSD [6]. Application à la conception d’une barrière ouvragée "sable de dune-
Comme perspectives des travaux de recherche, il reste à bentonite". Thèse de Doctorat, Tlemcen,116 pages, 24 janvier 2011.
confirmer le choix de l’argile par rapport aux valeurs du

ISSN 1737-9334
Pour une approche environnementale du traitement

des eaux usées de la ville de Tamanrasset -Sud
Algérien- dans le cadre du développement durable.
Bouaoud Mohamed#1, Bebouche Manel#2, Hadj Ali Zakaria#3
Département Techniques Urbaines et Environnement
Institut de Gestion des Techniques Urbaines
Université Salah Boubnider Constantine 3
Ville universitaire Ali Mendjeli Constantine, Algérie
Mohamed.bouaoud@univ-constantine3.dz
manelgtu1992@gmail.com
zakirroni7@gmail.com
Résumé — Suite à l’accroissement de l’utilisation quotidienne de

l’eau par les citoyens de nos villes et la diminution des ressources I. INTRODUCTION
en ce liquide précieux due au réchauffement climatique, et ça La ville de Tamanrasset, chef-lieu de la plus grande wilaya de
raréfaction dans les zones désertiques et arides, l’homme c’est
l’Algérie par sa superficie qui est plus grande que la France
retrouver dans l’obligation de trouver une alternative à ce
problème soit par les grands transferts (La grande rivière
métropolitaine, est une ville du grand Sahara où le problème
artificielle en Libye, le grand transfert In-Saleh /Tamanrasset en de l’eau était le souci de toutes les générations qui se sont
Algérie) ou bien la désalinisation de l’eau de mer (proche et succédées sur cette ville ; Les autorités publics Algériennes
moyen orient), mais ses deux solutions ne font qu’empirer la ont pu mettre fin à ce problème par un grand transfert de l’eau
situation environnementale par l’assèchement des nappes potable sur une distance de plus de 740Km pour satisfaire au
souterraines qui sont pour le cas du grand Sahara non besoin des citoyens de cette région.
renouvelable et par l’augmentation du degrés de salinité du Et pour une bonne gouvernance de cette eau et garantir une
littorale des mers en rejetant les saumures résidu de la gestion durable de cette ressource hydrique, l’installation de
désalinisation
station d’épuration est impérative pour le traitement et la
Tous cela nous mène à penser à d’autres alternatives pour réutilisation des eaux usées et pour ne pas se retrouver dans
préserver cette ressource pour les générations futures d’où la une catastrophe écologique.
réutilisation des eaux usées à nouveau par leur épuration.
II. DIAGNOSTIC DE LA SITUATION ENVIRONNEMENTALE
Malgré le raccordement de la majorité des quartiers de la ville
de Tamanrasset au réseau d’assainissement et la mise en marche L’Algérie est parmi les pays qui ont lancés un grand
de la station d’épuration en 2011, on trouve que le rejet d’eau programme de traitement des eaux usées (figure 1) avant de
épurée cause beaucoup de problèmes sanitaires et écologiques, les rejetés dans les oueds, rivières et milieu naturel ensuite
cela se traduit par la prolifération de maladies dans les zones où réutilisés ces eaux après la refermeture du cycle de l’eau.
passe le coure d’eau du rejet. On peut dire alors que le processus La ville de Tamanrasset est l’une des villes Algériennes
d’épuration est insuffisant.
caractérisée par sa diversité environnementale et la richesse
Pour améliorer le traitement des eaux usées dans la ville de
Tamanrasset on a opté dans ce travail pour faire un diagnostic
de son écosystème ; cette richesse est menacée par la non
de la situation ainsi que proposer quelques solutions et projets maitrise du domaine hydraulique, et plus précisément
d’avenir pour atteindre et assurer une bonne épuration des eaux l’assainissement des eaux usées et pluviales.
usées de cette ville dans le cadre du développement durable, Tout en étant dans une zone désertique et complètement
toute en améliorant la qualité environnementale et sans aride, chaque ressource hydrique doit être bien exploiter et
influencer sur la santé publique. préserver pour les générations futures.
Mots–clés — Traitement des eaux, Développement durable,

Santé, Ecologie, Tamanrasset.

ISSN 1737-9334
III. EFFET DE LA STATION D’EPURATION DE LA VILLE DE

TAMANRASSET SUR SON ENVIRONNEMENT
La station d’épuration « STEP » de la ville de
Tamanrasset a été mise en service le 1er Novembre 2011, où
la direction d’hydraulique de la wilaya a chargé l’office
nationale d’assainissement de sa gestion et exploitation, elle
se situe à la cité Tahagouine au Sud-Ouest de la ville et à
2Km de l’agglomération, elle a été implantée dans la
direction d’écoulement de Oued Tamanrasset pour faciliter
l’écoulement gravitaire des eaux avant et après traitement.
D’une superficie de 30 hectares et d’une capacité de
126000 m3, la STEP traite un volume journalier de 15500 m3
uniquement de nos jours.
Fig.1 : Volume d’eau épurée dans les stations d’épurations gérées par l’ONA
Malgré l’efficacité du traitement des eaux usées dans la
STEP de Tamanrasset, approuvée par le laboratoire d’analyse
Etant l’un des piliers majeurs du développement durable, la
qui, après chaque cycle de traitement fait un échantillonnage
préservation des ressources pour les générations futures nous
avant de déverser les eaux épurées dans l’oued Tamanrasset,
mène à penser au recyclage des eaux vannes et pluviales pour
cela suivant bien sûr les normes et lois exécutives émise par
les réutiliser principalement dans le secteur agricole.
le ministère d’hydraulique et ressources en eau ou même
Malgré les efforts émis par l’office nationale
décréter dans le journal officiel de la république Algérienne.
d’assainissement « ONA » on trouve que cette ville souffre
Mais au fil du temps et avec la dégradation des
toujours des problèmes écologiques et sanitaires, surtout ceux
équipements mal entretenus de la STEP, nous avons constaté
en relation directe avec l’eau, soit en milieu urbain à cause
qu’il y a plusieurs défaillances au niveau de cette structure :
des fosses septiques mal ou non entretenues et le mal
dimensionnement du réseau d’assainissement ainsi que sa  Au niveau du prétraitement on trouve que les mailles du
vétusté par endroit (figure 2), soit dans les villages avoisinant dégrilleur laissent passer quelques déchets ainsi que
le cours d’eau où passe le rejet de la station d’épuration de la l’absence totale du déshuileur et par conséquence le
ville (figure 3). passage des graisses et huiles dans les bassins d’aération
(figure 4).
Fig.2 : Photos montrant les nuisances causées par le réseau d’assainissement

vétuste et mal dimensionné en milieu urbain Fig.4 : Photos montrant les huiles flottant dans les bassins d’aération
 Au niveau des bassins de décantation parfois on ne donne

pas assez de temps aux eaux usées pour qu’elles puissent
se décanter convenablement.
 Au niveau des bassins d’aération on trouve le problème de
l’entretien des aérateurs et la difficulté de réparation en
cas de pannes.
Et pour appuyer notre hypothèse sur le mal

fonctionnement de la STEP, nous avons mené une enquête
sur le terrain et plus précisément sur le village qui suit
directement la STEP dans la direction de l’écoulement des
eaux rejetées.
Notre enquête a visé les habitants du village de
Fig.3 : Photo montrant le rejet de la station d’épuration de la ville de « Tagrambait » situé à 15 Km de la STEP, elle a pour but de
Tamanrasset connaitre le degré de répercussion sanitaire,

ISSN 1737-9334
environnementale et même économique sur les villageois, des

eaux usées ayant pour origine le rejet de la STEP.
 Commençant par les nuisances causées par les mauvaises
odeurs ou on trouve que 60% des habitants (figure 5)
disent que ce problème vient des eaux polluées qui passe
par l’oued (figure 6).
Fig.8 : Photo montrant la pollution des eaux des puits du village de

« Tagrambait ».
 La pollution a aussi causé un problème économique pour

les agriculteurs de la région après 2011, cela est dû à la
Fig.5 : figure montrant la répartition des habitants sur la question de la mauvaise réputation que leurs produits ont eu à cause de
prolifération des mauvaises odeurs ces eaux usées traitées, surtout que la majeure partie des
eaux est utilisé pour l’irrigation (figure 9).
Fig.6 : Photos montrant la pollution des eaux de l’Oued malgré la mise en

service de la STEP
 L’infiltration des eaux de l’oued dans les puits du village,

car la totalité des habitants souffrent de la pollution de
Fig.9 : figure montrant la répartition des utilisations des eaux de puits
leurs puits (figure 7 & 8) dû au rejet de la STEP, surtout
que les habitants se sont implantés beaucoup plus avant  Du point de vue effets sur la santé publique on trouve
2011 la date de mise en service de la STEP et disent qu’ils principalement les maladies de la peau, ensuite viennent
n’avaient pas ce problème auparavant, surtout qu’ils les maladies respiratoires (Figure 10 & 11).
utilisent largement les eaux de puits pour l’irrigation de
leurs parcelles agricoles.
Fig.10 : Photos montrant la souffrance des habitants du village des maladies

Fig.7 : figure montrant la répartition des habitants sur la question dermatologique.
d’infiltration des eaux polluées aux puits de la région

ISSN 1737-9334
 Ce lac sera conçu de façon qu’il devient un réservoir

d’énergie hydraulique, vu que notre terrain est accidenté
et qu’on peut avoir une inertie hydrique pour produire de
l’électricité pour tous les équipements du site.
 L’installation de poteux pour l’éclairage publique
utilisant l’énergie photovoltaïque.
 Espaces verts irrigués avec les eaux traitées ainsi que des
aires de jeux pour les enfants.
 Un café-restaurant pour le bien être des visiteurs ainsi
qu’un parking
REFERENCES
Fig.11 : figure montrant les effets de la pollution sur la santé des habitants du
village de « Tagrambait ». [1] Z. Hadj Ali, M. Bebouche, “ ‫ دراسة حالة مدينة‬،‫أثر محطة تصفية المياه المستعملة‬
-‫قرية تاقرمبايت‬-‫تمنراست‬,” Mémoire de Master, Institut de Gestion des
Techniques Urbaines, Université Salah Boubnider, Constantine3, Juin
IV. CONCLUSIONS ET PROPOSITIONS POUR L’AMELIORATION 2018.
DU CADRE DE VIE [2] I. Labed, “‫”معالجة مياه الصرف الصحي لمنطقة تقرت بواسطة نباتات منقية محلية‬
Pour l’amélioration du cadre de vie des habitants de la ville Thèse de Doctorat, Faculté des mathématiques et sciences de la
matière, Université Kasdi Merbah de Ouargla, 2015.
de Tamanrasset, nous avons proposé une aire de repos et de [3] A.B.Bisinella de Faria “Development of an integrated approach for
traitement naturel des eaux épurées rejetées par la STEP, tout wastewater treatment plant optimization based on dynamic modelling
en respectant la nature de la région et le milieu saharien dans and environmental assessment ” Thèse de Doctorat, INSA Toulouse,
lequel elle se trouve cela avec l’utilisation des énergies Oct 2016.
[4] F. Dutertre, “Réutilisation des eaux usées traitées pour l’irrigation
renouvelables et les concepts du développement durable. d’un golf selon l’arrêté du 2 Aout 2010: cas de la station d’épuration
de Dax,” Polytech Montpelier 14éme séminaire de l’eau, , p. 215,
Cet espace se trouve près du quartier de Tahagouine au nord- 2015.
ouest de la station d’épuration directement implanté après le [5] A. Nacer, “‫ ”حكم التطهير و استعمال المياه العادمة في الفقه اإلسالمي‬Mémoire de
rejet, d’une surface de 5.8ha il est composé d’ : Magister, Faculté des études supérieures, Université nationale
Ennadjah de Nablous, Palestine, 2012.
 Un lac artificiel d’1.5ha où on a instauré des plantes [6] H. Laribi, “‫تلوث المياه و دور محطات التطهير في التقليل من التلوث المائي‬,”
aquatiques spéciales pour le traitement des eaux usées, Mémoire de Master, Faculté des sciences sociales, Université
Abdelhamid Ibn Badis, Mostaganem, Juin 2011.
d’où la mise en place d’un écosystème qui va limiter la
pollution de l’eau ainsi que la propagation des mauvaises
odeurs et la prolifération des maladies microbiennes et
bactériologiques.
Copyright - 2019
Page 83
ISSN 1737-9334
Biographie de l’éditeur:
Dr. Ahmed Rhif (Tunisia) is a Researcher & Engineer (PhD,

Eng). He has more than 11 years of experience on Scientific
Research, Teaching and industrial projects. He is actually the
Dean of the International Centre for Innovation &
Development (ICID). Ahmed Rhif has worked as a Technical
Responsible Chief in LEONI (International Leader of Wiring
Fibers Companies) and has occuped also the task of Project
Manager and Method Engineer in both SMSI (electronic
development industry) and CABLITEC (Engineering
automobile company). Then he was a Lecturer at both the Private University of
Sousah (UPS) and the High Institute of Applied Sciences and Technologies of
Sousah (ISSATso) and now he is working as Lecturer in the High Institute of
Applied Sciences and Technologies of Al Qayrawan (ISSATk). His research interests
include Modelling, Control Systems and Engineering as well as the implantation of
the international standard of quality (ISO-TS 16949). Dr. Ahmed Rhif is also a
Trainer of Trainers certified from the Canadian Centre of Training in Montreal and
Coach certified NLP from the American Board of Neuro-linguistic programming.

Energie Renouvelable Et Applications Francophone Issue

Transféré par

Droits d'auteur :

Formats disponibles

Energie Renouvelable Et Applications Francophone Issue

Transféré par

Informations du document

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

Energie Renouvelable Et Applications Francophone Issue

Transféré par

Droits d'auteur :

Formats disponibles

ISSN: 1737-9334

Proceedings of Engineering & Technology

Energie Renouvelable et Applications

International Centre for Innovation & Development

Dr. Ahmed Rhif (Tunisie) Sallam Mebrouk (ALG)

Georges Descombes (FRA)

Hana Boudjedir (ALG)

Fawaz Massouh (FRA)

Nachida K. Merzouk (ALG)

Khenfer Nabil (ALG)

Mohamed Elamine Djeghlal (ALG)

Manar Slimane (MAR)

Abdelfettah Barhdadi (ALG)

Mohamed Gherbi (ALG)

Mohamed Benbouzid (FRA)

Sara Zatir (ALG)

Philippe Guibert (FRA)

Mustapha Hatti (ALG)

Saoussen Hammami (TUN)

Rahmani Lazhar (ALG)

Tahar Bahi (ALG)

Commande d’un système éolien basé sur une MADA. Page 5

Stratégie de contrôle du raccord au réseau du générateur photovoltaïque. Page 20

Le bâtiment passif pour assurer l'efficience énergétique. Page 35

Caractérisation physico-chimique et mécanique des argiles de la région d’Adrar et de Timimoun Page 76

Analyse comparative des perturbations EM conduites

Copyright - 2019 Page 1

Pour le contrôle du taux de distorsion harmonique (THD)

influences EM dans la gamme des fréquences comprises dans 

Les ponts redresseurs monophasés présentent un intérêt très

Dans ce paragraphe, nous allons examiner l’influence des

Pont redresseur triphasé non commandé

FONDAMENTAL RELATIF AU PONT REDRESSEUR NON COMMANDE

Fondamental (V) 530.5

Courant de Simulation 86.92

A la lumière de ces, il est facile de voir l’effet perturbateur

Enfin nous remarquons, que le pont redresseur triphasé non

commandé est moins polluant que le pont redresseur triphasé

Fig. 5 Spectre fréquentiel du courant de ligne

Commande d’un système éolien basé sur une MADA

Mots clés: Eolienne - Machine Asynchrone à Double

C_p (λ,β)=(0.3-0.00167*β) sin⁡[π(λ+0.1)/(10-0.3β)]- Avec :

1 Pour les équations électriques :

La puissance mécanique disponible sur le rotor P_m est

R : Rayon de l’aérogénérateur 35.25𝑚 Pour les équations de flux :

ρ : Densité de l’air 1.225 𝐾𝑔. 𝑚−1

Compte tenu du rapport du multiplicateur de vitesse, la  ds  Ls I ds  MI dr (7)

Ls  I s  ms L’inductance cyclique du stator ; E. Modèle complet du système de conversion éolien

Copyright - 2019 Page 7

fréquences variables. Plus la vitesse du vent est plus

Fig. 8 Schéma bloc du modèle de la turbine éolienne (Matlab – Simulink).

Fig. 5 Profil du vent appliqué (m⁄s)

Fig. 9 Variation du Cp en fonction du temps

Fig. 6 Schémas bloc du pitch control

Fig. 7 la variation de l’angle de calage avec le temps

Fig. 12 Variation de la vitesse du rotor.

Copyright - 2019 Page 8

La figure (16) présente le schéma bloc de cette machine

D. Association de la MADA avec un onduleur à MLI

Fig. 156 Schéma bloc de la MADA en génératrice dans le référentiel (d-q)