Annexes
Annexe 1 : Implémentation des
opérateurs génétiques 55
Annexe 2 : Résultats du bilan thermique
de la chambre expertisée 57
Annexe 3 : Résultats de diagnostic d'un
groupe frigorifique du compartiment principal 58
Annexe 4 : Fonctions exportées par la DLL
des fluides frigorigènes 60
Annexe 5 : Quelques règles fournies avec
la CCF implémentée 61
Annexe 1 : Implémentation des opérateurs
génétiques
Ce code a été testé sur le compilateur
Visual C++ de Microsoft.
// (c) août 2003 OUAMBO TOBOU Raoul. #include
<stdlib.h>
//Longueur en bits du type d'entier utilisé dans le codage
const size_t lng = sizeof(unsigned short)*8;
//Entier max du codage sur 16 bits const unsigned short gmax =
0xFFFF;
//Génère un nombre aléatoire compris entre 0
et 1
inline double Aleatoire(void)
{
return (double) rand()/RAND_MAX;
}
//Fonction de décodage
inline double Decode(unsigned short g, double Xmin, double DX)
{
return Xmin + DX*g;
}
// Sélection par tournoi avec une probabilité de
sélection p
// entre deux individus de fitness respectifs Fit1 et Fit2
inline short SelTournoi(double Fit1, double Fit2, double p)
{
short best;
best = (Fit1 < Fit2) ? 1 : 2;
return (Aleatoire() <= p) ? best : 3 - best;
}
//Croisement en 2 points de *parent1 et *parent2 et leurs
probabilités de mutation inline void Croise2(unsigned short *parent1,
unsigned short *parent2,
unsigned short *fils1, unsigned short *fils2,
unsigned short *proba01, unsigned short *proba02,
unsigned short *proba11, unsigned short *proba12)
{
unsigned short mask1, mask2;
size_t pos1, pos2;
//Postion de la 2ème cission du croisement pos2 = 2 +
(size_t) (Aleatoire() * (lng - 2)); if(pos2 == lng) pos2 = lng - 1;
//Postion de la 1ère cission du croisement
pos1 = 1 + (size_t) (Aleatoire() * (pos2 - 1)); if(pos1 == pos2)
pos1 = pos2 - 1;
// Construction des masques de croisement
mask2 = 0xFFFF;
mask2 <<= (lng - pos2 + pos1);
mask2 >>= (lng - pos2);
mask1 = ~mask2;
// Croisements des individus
*fils1 = (mask1 & *parent1) | (mask2 &
*parent2); *fils2 = (mask2 & *parent1) | (mask1 & *parent2);
// Croisements des probabilités de mutation *proba11 =
(mask1 & *proba01) | (mask2 & *proba02); *proba12 = (mask2 &
*proba01) | (mask1 & *proba02);
}
//Mutation d'un individu et de sa probabilité de
mutation,
//avec une probabilité uniforme p
inline void Mute(unsigned short *ind1, unsigned short *ind2,
unsigned short *pro1, unsigned short *pro2, double p)
{
unsigned short masque, buf;
//Construction du masque de mutation masque = 0;
for(size_t i=0; i < lng; i++)
{
buf = (Aleatoire() <= p) ? 1 : 0; buf <<= i;
masque |= buf;
};
//Mutation bit par bit par l'opération XOR
*ind2 = *ind1 ^ masque; *pro2 = *pro1 ^ masque;
}
//Extrait l'individu N°num décodé
(varialbes d'état) //du tableau codé de la population
inline void Extrait(
unsigned short *popu, //Population
double *var, //Individu (un nvar-uplet de variables)
double *min, //Borne min des varialbles
double *pas, //Pas de discrétisation des intervalles de
variables
int num, //Numéro de l'individu dans la population
int nvar //Nombre de variables dans un individu
)
{
int j, prod;
prod = nvar * num;
for(j=0; j < nvar; j++)
{
};
}
var[j] = Decode(popu[prod + j], min[j], pas[j]);
Annexe 2 : Résultats du bilan
thermique de la chambre expertisée
RAPPORT DU CALCUL DE BILAN THERMIQUE
Date d'édition du rapport : 29/12/2003 Nom du projeteur :
Ouambo
Date de début du projet : 18/10/2003 Client
concerné : Congelcam
Chambre froide : Bonabéri - Ancienne
Pays d'installation de la chambre : Cameroun Ville d'installation
de la chambre : Douala Profil climatique utilisé : Février
Charges du Compartiment 1
Volume du compartiment : 1920 m3
Température désirée dans le compartiment :
-15 °C Humidité relative désirée dans le
compartiment : 85 % Charges à travers les parois verticales : 2801,067 W
Charges à travers les ouvertures : 8853,183 W
Charges à travers le plafond : 1962,562 W
Charges à travers le plancher : 2420,563 W
Charges par renouvellement d'air : 0 W
Charges dues aux denrées entrantes : 4698,958 W Charges
dues à la respiration des denrées : 0 W Charges dues à
l'éclairage : 360 W
Charges dues aux appareillages : 750 W
Charges dues au personnel : 540 W
Charges dues aux moteurs d'évaporateurs : 2656,25 W
Charges dues au dégivrage : 3487,5 W
Puissance frigorifique intermédiaire : 28530,082 W Temps
de fonctionnement des installations : 18 heures
Puissance frigorifique à fournir par les
évaporateurs : 38040,11 W
Charges du Compartiment 2
Volume du compartiment : 8000 m3
Température désirée dans le compartiment :
-20 °C Humidité relative désirée dans le
compartiment : 90 % Charges à travers les parois verticales : 6134,105 W
Charges à travers les ouvertures : 2065,521 W
Charges à travers le plafond : 8998,358 W
Charges à travers le plancher : 11766,624 W
Charges par renouvellement d'air : 0 W
Charges dues aux denrées entrantes : 28193,75 W Charges
dues à la respiration des denrées : 0 W Charges dues à
l'éclairage : 1440 W
Charges dues aux appareillages : 1312,5 W
Charges dues au personnel : 1950 W
Charges dues aux moteurs d'évaporateurs : 10625 W Charges
dues au dégivrage : 13950 W
Puissance frigorifique intermédiaire : 86435,858 W Temps
de fonctionnement des installations : 18 heures
Puissance frigorifique à fournir par les
évaporateurs : 115247,81 W
Annexe 3 : Résultats de diagnostic
d'un groupe frigorifique du compartiment principal
RESULTATS DE LA SIMULATION INTERNE DU GROUPE 1
Date d'édition du rapport : 28/12/2003 Nom du projeteur :
Ouambo
Date de début du projet : 18/10/2003 Client
concerné : Congelcam
Chambre froide : Bonabéri - Ancienne
Pays d'installation de la chambre : Cameroun Ville d'installation
de la chambre : Douala Profil climatique utilisé : Février
PARAMETRES D'ENTREE DE L'ALGORITHME GENETIQUE Probabilité
de sélection : 0,9
Probabilité minimale de mutation : 0,001
Probabilité maximale de mutation : 0,005 Taille de la population : 20
Nombre de génaration simulées : 400
PARAMETRES DE SORTIE DE L'ALGORITHME GENETIQUE Paramètre
d'homogénisation de la population finale : 35 Fitness minimal :
0,92486
Durée des calculs : 31,483 minutes
PARAMETRES DU GROUPE
Fluide frigorigène : R404A
Numéro Evaporateur = 2
Diamètre tuyauterie aspiration (m) = 0,066675
Diamètre tuyauterie de vapeurs surchauffées (m) =
0,022225
Diamètre tuyauterie de liquide (m) =
0,022225 Puissance maximale de compression (W) = 30000
MESURES REALISEES
Température entrée air évaporateur
(°C) (obligatoire) = -17 Température sortie air
évaporateur (°C) (obligatoire) = -18,5
Température entrée air condenseur (°C)
(obligatoire) = 26 Température sortie air condenseur
(°C) (obligatoire) = 30 Pression aspiration compresseur (bars)
(obligatoire) = 2,5 Pression refoulement compresseur (bars) (obligatoire) = 18
Température aspiration compresseur (°C) (obligatoire) =
2 Débit volumique air évaporateur (m3/s) (obligatoire) = 7,5
Débit volumique air condenseur (m3/s) (obligatoire) = 10 Humidité
relative compartiment (%) (obligatoire) = 85 Température refoulement
compresseur (°C) = 70
MESURES OBTENUES PAR CALCUL
Température entrée air évaporateur
(°C) (obligatoire) = -17 Température sortie air
évaporateur (°C) (obligatoire) = -18,49951
Température entrée air condenseur (°C)
(obligatoire) = 26 Température sortie air condenseur
(°C) (obligatoire) = 26,90332 Pression aspiration compresseur
(bars) (obligatoire) = 2,68275 Pression refoulement compresseur (bars)
(obligatoire) = 17,97653 Température aspiration compresseur
(°C) (obligatoire) = 1,58683 Débit volumique air
évaporateur (m3/s) (obligatoire) = 7,5
Débit volumique air condenseur (m3/s) (obligatoire) = 10
Humidité relative compartiment (%) (obligatoire) = 85 Pression
entrée évaporateur (bars) = 3,07802
Pression sortie évaporateur (bars) = 3,01253
Pression entrée condenseur (bars) = 15,92337 Pression
sortie condenseur (bars) = 16,06997 Température entrée
évaporateur (°C) = -20,13993 Température sortie
évaporateur (°C) = -20,13993 Température
entrée condenseur (°C) = 62,11033 Température
sortie condenseur (°C) = 21,82979 Température
refoulement compresseur (°C) = 70,08289
VARIABLES D'ETAT OBTENUES PAR CALCUL Rendement isentropique
compresseur = 0,98854 Efficacité évaporateur = 0,46396
Surchauffe enthalpique évaporateur = -0,011 Surchauffe
aspiration (°C) = 21,72676 Efficacité condenseur =
0,97547
Sous-refroidissement enthalpique condenseur = 0,15341
Débit massique frigorigène (kg/s) = 0,24307 Température
évaporation (°C) = -20,13993
Température condensation (°C) =
34,57298
Pertes de charges ligne d'aspiration (bar) = 0,32978 Pertes de
charges ligne liquide (bar) = 1,79072
Pertes de charges refoulement - entrée condenseur (bar) =
2,05316 Désurchauffe refoulement - entrée condenseur
(°C) = 7,97255 Sous-refroidissement ligne liquide
(°C) = 29,0358
VARIABLES INTERNES OBTENUS PAR CALCUL
Vitesse entrée tuyauterie aspiration (m/s) = 4,46199
Vitesse sortie tuyauterie aspiration (m/s) = 5,70555
Vitesse entrée tuyauterie vapeurs surchauffées
(m/s) = 8,10915 Vitesse sortie tuyauterie vapeurs surchauffées (m/s) =
9,03205 Vitesse entrée tuyauterie liquide (m/s) = 0,59127
Vitesse sortie tuyauterie liquide (m/s) = 0,53188
Production frigorifique (W) = 39451,10032
Chute de pression entrée évaporateur (bar) =
11,05463 Surchauffe évaporateur (°C) = 0
Sous-refroidissement condenseur (°C) =
12,74319
Volume spécifique entrée compresseur (m3/kg) =
0,08196 Volume spécifique refoulement compresseur (m3/kg) = 0,01294
Volume spécifique sortie condenseur (m3/kg) = 0,00094
Annexe 4 : Fonctions exportées par la
DLL des fluides frigorigènes
Nous présentons ci-après un extrait de l'aide
accompagnant les versions 2.0 des bibliothèques de fonctions fournies
par la société SOLVAY FLUOR UND DERIVATE.
The functions of the following properties are implemented:
1. bubble and dew pressure: p' = f(T) und p» = f(T)
2. bubble and dew temperature: T' = f(p) und T» = f(p)
3. specific volume of the liquid phase: v' = f(T)
4. pressure, vapour phase: p = f(T, v)
5. specific volume vapour phase: v = f(p, T)
6. temperature vapour phase: T = f(p, v), T = f(p, h) und T =
f(p, s)
7. specific enthalpy liquid phase: h' = f(T)
8. specific enthalpy vapour phase: h = f(T, v) und h = f(T,
p)
9. specific entropy liquid phase: s' = f(T)
10. specific entropy vapour phase: s = f(T, v) und s = f(T,
p)
11. specific heat capacity liquid phase: c' = f(T)
12. specific heat capacity vapour phase (p=const.): cp = f(T, v)
und cp = f(T, p)
13. specific heat capacity vapour phase (v=const.): cv = f(T, v)
und cv = f(T, p)
14. adiabatic exponent, vapour phase: K = f(T, v) und K = f(T,
p)
15. velocity of sound, vapour phase: w = f(T, v) und w = f(T,
p)
16. thermal conductivity, liquid phase: X ` = f(T)
17. thermal conductivity, vapour phase: X = f(T, p)
18. dynamic viscosity, liquid phase: i ` = f(T)
19. dynamic viscosity, vapour phase: i = f(T, p)
20. surface tension: a = f(T)
21. characteristic data (molar mass, critical temperature,
critical pressure,...)
|
|