WOW !! MUCH LOVE ! SO WORLD PEACE !
Fond bitcoin pour l'amélioration du site: 1memzGeKS7CB3ECNkzSn2qHwxU6NZoJ8o
  Dogecoin (tips/pourboires): DCLoo9Dd4qECqpMLurdgGnaoqbftj16Nvp


Home | Publier un mémoire | Une page au hasard

 > 

Développement d'un logiciel d'expertise technique d'installation frigorifiques de chambre froides

( Télécharger le fichier original )
par Raoul Ouambo Tobou
ENSAI de Ngaoundere - Cameroun - Ingenieur 2003
  

précédent sommaire suivant

Extinction Rebellion

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,...)

précédent sommaire suivant






Extinction Rebellion





Changeons ce systeme injuste, Soyez votre propre syndic





"Les esprits médiocres condamnent d'ordinaire tout ce qui passe leur portée"   François de la Rochefoucauld