Chapitre3 : CHOIX DES EQUIPEMENTS

III. CHOIX DES EQUIPEMENTS :

1. Introduction

Le choix d’équipements consiste à choisir un groupe se rapprochant le plus

possible de la puissance frigorifique effective de l’enceinte à refroidir, compte
tenu des conditions d’utilisation du matériel.

Suivant la gamme offerte par le constructeur, il faut prendre de préférence un

groupe ayant une puissance légèrement au-dessus de celle désirée plutôt qu’au-
dessous .Il développera ainsi cette puissance en un temps de fonctionnement
légèrement réduit.

Les principaux matériaux sont : le compresseur, le condenseur, le détendeur,

l’évaporateur.
Les appareils annexes sont : le filtre des hydrateur, voyant liquide, robinet
solénoïde,

2. Méthode de choix de l’évaporateur :

Pour choisir un évaporateur d’une chambre froide, il doit connaître :

 La puissance frigorifique [kcal/h] ou [watt].(11,5 kW)
 L’écart de température entre la chambre froid et la température
d’évaporation.
 Nature de fluide frigorigène à utiliser.(R404a)
 La puissance la plus proche du catalogue de constructeur à celle qu’on
a calculé.
- Détermination de ∆Tévap :
L’humidité relative pour la conservation de pomme de terre comprise entre
85% à 90% (Annexe : 2.9)

Hr = 85%⇒∆Tévap= 5°C,Hr = 90%⇒∆Tévap= 7°C (Annexe 3.1)

1
Avec : Ta : température de la chambre froide = 7°C

∆Tévap =5c° ⇒ Te=Ta-∆T =7-5=2c°

∆Tévap =7c° ⇒ Te=Ta-∆T =7-7=0c°

Avec la puissance frigorifique de l’évaporateur calcule précédemment et le

∆Tévapcomprise entre 5c°et 7c°, nous allons retirer les valeurs du catalogue

constructeur Friga Bohn pour choisir l’évaporateur qui se rapproche le mieux de

notre puissance frigorifique.(Annexe : 3.2)

Pour ∆Tévap =5c° on Q = 8kW

Pour le R22
Et pour ∆Tévap =7c° on Q = 11kW

Le facteur de correction du fluide R404a =1,06.

∆Tévap =5c° on Q = 8 kW . 1,060 = 8.48kW (Annexe 3. 3)

∆Tévap =7c° on Q = 11 kW.1,060= 11.66kW (Annexe 3.3)

On aura la référence de l’évaporateur : SK11E

3- Méthode du choix de groupe de condensation :

Pour la détermination d’un groupe de condensation qui convient notre

installation il faut d’abord connaitre : (Annexe3.4)

 Puissance frigorifique [kcal/h] ou [watt].

 La température de l’air de condenseur (°c).
 Nature de fluide frigorigène.
 La température de saturation a l’aspiration (T sa) qui est la différence entre
la température d’évaporation et l’abaissement de la température comprise
généralement entre 1et3c°.
 3.1- Choix du groupe compresseur – condenseur :

Le choix du groupe dépend essentiellement de la puissance frigorifique ; de la

température de saturation à l’aspiration (Tsa) ainsi que de la température de l’air
au condenseur.

 Calcule de la température de saturation à l’aspiration(Tsa) :

∆T (aspirations) = Te –Tas : arbitrairement ont considéré que la différence de

température entre l’évaporateur et la température saturation à l’aspiration est de
2c°mais cette valeur sera corrigé par la suite :

Donc : 2c°= Te- Tsa

*quand Te= cela implique que :

Tsa= te-2 Tsa= 2-2 Tsa=0c°

*quand Te = 0 cela implique que :

Tsa=te -2 Tsa = 0-2 Tsa=-2c°

La température de l’air de condenseur est d’environ 33c° à DJENDEL après

interpolation on relève de catalogue BITZER les valeurs suivantes :

 Interpolation Température (0) la zone de température 33c° :

(Annexe3.4)

Te(0) = y1+(𝑦2−𝑦1) ×(x-x1)

(𝑥2−𝑥1)

Te =12,94+ (11.95−12.94) ×(33-27)

(0) (32−27)

Te (0) = 12.94 -1,188

Te(0)=11.76kw
  Interpolation Température (-5) la zone de température 33c° :

Te(-5) = y1+(𝑦2−𝑦1) ×(x-x1)

(𝑥2−𝑥1)

Te(-5) =11,08+ (10.23−11.08) ×(33-27)

(32−27)

Te (-5) =11,08-1.02

Te(-5) =10.06kw

 Interpolation Température (-2) la zone de température 33c° :

Te (-2)= y1+(y2−y1) ×(x-x1)

(X2−X1)

Te =11,76+(10.06−11.76) ×(-2-0)
(-2)
(−5−0)

Te (-2)=11,76-0,68

Te (-2)= 11.08kw

Tsaest donc compris entre :

0c° et -2c°

On aura alors du catalogue BITZER le LH84/2CC-4.2Y

Tsa= 0c° on a une puissance Ø=11,76kw

Tsa = -2c°on a une puissance Ø=11.08kw

3.2- Caractéristiques dugroupe de condensation :

 Puissance. Frigorifique : 11.76KW

 Mode Réfraction et air conditionné :LH84/2CC-4.2Y
 Fluide frigorigène : R404A
 Courant max compresseur : 16.4/9.4 A
  courant max du ventilateur condenseur : 3.08 A
 Débit d’air condenseur : 4577m3/h
 Charge maximal de fluide : 12.5 kg
 Point (stand) : 134kg

4- Equilibrage évaporateur et groupe de condensation :

On trace deux segments caractérisant l’évaporateur et le groupe de

condensation à l’aide de deux différents points de puissance frigorifique en
fonction de ∆T et tsa, l’intersection entre les deux segments donnera le point de
fonctionnement E.

A partir du point E on aura : la puissance réelle d’équilibrage ФEQ, la

température d’évaporation te, température saturante à l’aspiration tsa et la
différence de température ∆T.
 4.1 Graphe d’équilibrage :

 Puissance équilibrage : 11.1 kw

 ∆T équilibrage : 6.8c°

 Température évaporation équilibrage : 0.8c°

 Température saturation a l’aspiration : -1.8c°

5-Choix du détendeur :

Pour sélectionner le détendeur on doit connaitre les critères suivant :

 La puissance de l’installation au point d’équilibrage.

 La température d’évaporation au point d’équilibrage.

 La nature du fluide frigorigène.

 Il fout connaitre également la chute de pression du régime au

détendeur ∆Pd.

5.1- Perte de charge au détendeur :

∆Pd =Pk -∆PL -Pe

∆Pd : chaut de pression au détendeur.

Pk : pression de condensation correspondant (Tk)

Pe : pression d’évaporation correspondant (ta)

∆PL : perte de charge totale dans la conduite liquide.

5.2- perte de charge totale dans la conduite liquide :

C’est la somme des pertes des charges dynamiques singulières les pertes les
pertes de charges dynamiques linéaires ;la perte de charge statique.

∆PL =∆Pd +∆Pst

Avec : ∆Pd =∆Pd sin + ∆Pd Lin

∆Pd Lin : pertes de charge dynamiques linéaires.

∆Pd sin : pertes de charge dynamique singulière

 - Perte de charge dynamique singulière :
On calcul la perte de charge dynamique singulière en même temps qu’on
sélectionne l’organe correspondant ; et on détermine la perte de charge
dynamique singulière de chaque organe.

Maintenant pour faire le choix de détendeur ; il nous faut connaitre la chute de

pression du régime au détendeur qu’on appellera ∆pd.

∆Pd = p entrée – p sortie

P entrée = pk -∆pl

∆pl : perte de charge totale sur la conduite de liquide HP.

P sortie = PE (la pression d’évaporation)

∆Pd = p entrée – p sortie = pk -∆pl –pe.

A partir de la température d’évaporation équilibrage ; on déduit la pression

d’évaporation des tables de tensions de vapeurs du fluide frigorigène r404a

 Interpolation : (AnnexE3.5)

Pe= y1+(y2−y1).(x-x1)
(X2−X1)

Pe= 6,04+(6,44−6,04) ×(0.8-2)

(0−2)

Pe=6.04+0.24

Pe=6,28bar

Te=0 ,8c° après interpolation pe = 6 ,28 bar

  Calcule Pk : (la pression de condensation)

TK=T air cd +∆𝑡tableau

∆T: approximatif au condenseur d’un groupe à air .

Avec Tsa=-1;8c °−20c° < Tsa < 0𝑐° (Annexe3.6)

 Interpolation : (Annexe3.6)

∆T= y1+(y2−y1) ×(x-x1)

(X2−X1)

∆T= 16+(13−16) ×(33-25)

35−25

∆T=16- 2.4 =13,6c°

Tk=33 + 13,6 Tk= 46,6c°

 Interpolation :(AnnexE3.5)

PK= y1+(y2−y1) ×(x-x1)

(X2−X1)

PK= 20,95+(21,95−20,95) ×(46,6-46)

48−46

PK= 20,95+0,3

PK= 21,25bar

ΔPL (Calcule De La Perte De charge Dans La Conduite Liquide) :

ΔPL =ΔP dynamique; singulière + ΔP dynamique; linéaire +ΔP statique

 ΔP dynamique singulière:

Ce sont respectivement :

1 - le filtre Déshydrater

2 - voyant liquide

3– robinet solénoïde
 1- Filtre Déshydrater :
On se base sur deux points pour le choix du filtre déshydrater:

 La puissance nominal (cap n) proche à la puissance de l’installation

 Le diamètre de la conduite liquide tiré du catalogue du groupe.

ΔP nominale = 14kpa ( à partir du catalogue constructeur) (Annexe3.7)

cap (n) : capacité nominale filtre : 20kw (Annexe3.7)

Dimension de la cuivre :( 1/2”.1/2”) (Annexe3.7)

ΔPf =ΔPn. (Øeq / cap (n))²

ΔPf =14.(11,1/20)²

ΔPf = 4,22kpa

On choix le filtre déshydrater suivant : DCL163

2-Voyant Liquide :
Pour le choisir on doit connaitre:

 Le diamètre de conduit liquide.

 Perte de charge nominale singulière pour le voyant est de 10kpa.

ΔPv= 10kpa (sans correction) (Annexe3.8)

SGN10 (1/2” ,1/2”)

3-Robinet solénoïde (électrovanne) :

Pour le choisir on doit connaitre:
 Le diamètre de la conduit liquide.
 La puissance nominaleØn qui est proche de la puissance d’équilibrage
Øeq
 ΔPRSL =15kpa (à partir du catalogue constructeur).(Annexe3.9)
Diamètre de cuivre : (1/2 ”,1/2”)
ΔPRSL =ΔPn. (Øeq / cap (n))²
cap (n) : capacité nominale robinet solénoïde(Annexe3.9)
ΔPRSL = 15.(11 ,1 / 18,9)²
ΔPRSL =5,16kpa

Le choix de l’électrovanne avec bobine se fait sur un catalogue DANFOSS :

ETS50B

Donc : ΔP dynamique singulière = ΔPf + ΔPv + ΔPRSL

ΔP dynamique singulière = 4;22+ 10+5 ,16=19,38 kpa

 ΔP dynamique linéaire:

On cas de présence d’un échangeur

 L’échangeur de chaleur :

la perte de charge à l’échangeur sera exprimée en longueur équivalente en métre

donc elle est ajoutée aux pertes de charge linéaire :

 Coté liquide :

Perte de charge assimilée à une longueur L=3m (du catalogue du constructeur)

 Puissance frigorifique d’équilibrage Øeq.

 Øn : la puissance nominale de l’échangeur de chaleur.

ΔP(ECH) = ΔPNL.(Øeq / cap (n))²

On approxime la conduite de liquide 8m environ

Øeq=11;1kw la conversion en kcal /h

C’est (11;1*3;6)/4.18

Øeq=11,1kw = 9,55kcal/h

 à partir de l’abaque des pertes de charges dans canalisation cuivrées ; on

peut déduire les pertes de charges dynamique linéaire.

Nous allons calcule la perte de charges dans l’échangeur (liquide – vapeur) cette
dernière est exprimée en longueur équivalente avec que rajouterons à la
longueur réelle de la conduite en fonction du diamètre liquide (3/8) et de la
puissance à l’équilibrage.
Øeq=9,55kcal/h

Et température d’évaporateur Te = 0.8c° on choisit l’échangeur (Annexe 3.10)

HE4.0

 la capacité nominale de l’échangeur (liquide –vapeur) est d’environ

(20.3,6)/4,18=17,22kcal/h

 Calcule de la perte de charge dans l’échangeur

ΔPNL : perte de charge nominale exprimée en métré coté liquide 3m

cap (n) : capacité nominale de l’échangeur (kcal /h)

ΔP(ECH) = 3 *(9,55/17,22)²=0 ;90m

 La longueur totale équivalente :

LTEqv =LLiq + ΔPECH(L)= 8+0;90

LTEqv = 8 ,90m

La perte de charge dynamique linéaire est ex primée par l’équation suivante :

ΔPD LIN=ΔPL 30.(L/30) .(100/1 ,019)

L : longueur du conduite liquide + longueur équivalente due à la perte de

charges de l’échangeur (coté liquide) + longueur équivalente due à la perte de
charge des accidents des parcours (en métré)

ΔPL 30 : perte de charges dynamique linéaires dans la canalisation cuivrée pour

le fluide frigorigène de l’installation pour 30m ;qui est déterminée d’après
l’abaque (Annexe3.11)

ΔPD LIN : en ( kpa)

ΔPD LIN = 0,64 .(8 ,90/30).(100/1,019)

ΔPD LIN = 18,21kpa

ΔP statique:

Pour r404a on approxime les pertes de charges statique Dans notre cas ; une
dénivellation de 3m environ cela veut dire que

L’évaporateur détendeur sont situés à un niveau supérieur par rapport au départ

de liquide.

ΔP statique=𝑝. g .h

: masse volumique de R404A

g: 10 kg/s2

h =3m

 Interpolation :(Annexe3.12)

𝑝= y1+(y2−y1) ×(x-x1)
(X2−X1)

𝑝=926,95+(920,21−926,95) ×(46;6-46)
(47−46)

𝑝 =926,95-4,044

𝑝=922,91kg/m3

ΔP statique = 922,91.10.3

ΔP statique =27687, 3psc = 27,68kpsc /m

Maintenant on peut calculer les pertes de charges dans la conduite liquide :

5.3-Perte de charges dans la conduite liquide :

ΔP cl= ΔP dynamique singulière + ΔP dynamique linière +ΔPstat

ΔP cl= 19,38+ 18,21+ 27,68

ΔP cl = 65,27kpa ΔP cl = 65,27/100 ΔP cl = 0,652bar

Nous pouvons maintenant calculer la chute de pression au niveau de détendeur

puisque nous avons calculé la pression à l’entrée c’est –à-dire

(PK-ΔPl) et la pression de la sortie (PE)

∆Pd =Pk -∆Pcl -Pe

∆Pd= 21,25 – 0,652– 6,28

∆Pd = 14,31 bar

Le choix de détendeur se fera en fonction de :

Øeq=9,55kcal/h

Te = 0,8c°

∆Pd = 14,31bar

∆Pd 14 ,31 de la zone de température(10c°) (Annexe 3.13)

∆Pd 14,31= 15;8 kw

∆Pd 14 ;23 de la zone de température(0c°) (Annexe 3.13 )

∆Pd 14,31= 14,9 kw

 Interpolation de Te = 0.8c° (Annexe 3.13)

Te0.8c°
= y1+(y2−y1) ×(x-x1)
(X2−X1)

Te0.8c° = 15.8+(14.9−15.8) ×(0.8-10)

(0−10)

Te0.8c° = 15.8-0,828=14.972

Te0.8c° = 14.972kw *3;6/4;18 Te0.8c° = 12,89kcal/h

Choix de détendeur type de vanne TES 5 - 3.7 (Annexe 3.13)

-Le surdimensionnement ou dimensionnement du détendeur :

(12.89 − 9.55.100)
9.55 = 34,97%
 Pour calcule la perte de charge dans la conduite d’aspiration on doit
connaitre :

 Le diamétre de la conduite d’aspiration.

 La température d’évaporation au point de d’équilibrage.

  La puissance frigorifique d’équilibrage ( Øeq)

La perte de charge dans la conduite liquide d’aspiration est exprimée par

l’équation suivant :

ΔPc asp=ΔPL 30. (L/30). (100/1.019)

L : longueur de la conduite d’aspiration + perte de charge de l’échangeur (coté

vapeur) .

ΔPL 30 : perte de charge dans la canalisation cuivrée pour le R404a par 30m ; qui
est déterminée d’après l’abaque.

On utilise l’abaque de perte de charges pour le fluide R404a on a :

Ø aspiration = 1 3/8 ” (Annexe3.14)

Longueur estimée 8m (la même que pour la conduite liquide)

Te = 0,8c° Øeq = 9,55 kcal / h

Avant de calcule la perte de charge au niveau de l’aspiration ; il faut d’abord

calculer la perte de charge au niveau de l’échangeur coté vapeur

∆Pnv = 1;5m (de l’abaque) ; c’est la perte de charge nominale exprimé en (m)
coté vapeur.

cap(n) = (20,3.6)/4,18=17,22kcal/h (Annexe 3.10)

∆Pech= ∆Pnv*( Øeq / cap(n))²

∆Pech = 1,5 * (9,55/ 17,22)²

∆Pech= 0,45m

La langueur totale équivalente :

L teqv = lvap+ ∆Pech

L teqv = 8 + 0,45

L teqv = 8,45 m
5.4-Calcule de a perte de charge au niveau de l’aspiration :

On relève de l’abaque 0.17 pour 30m pour la conduite d’aspiration (Annex3.11)

∆Pasp= 0.17. (8,45 /30).(100/ 1,019)

∆Pasp = 4,66kpa

A partir de cette perte de charge ; on peut corrige la déférence entre Te et Tsa

On avait pris une valeur arbitraire de 2c°, pour le R404a

Te de 0,8c° et ∆Pasp = 4,12 kpa pour le R404a

On relève de l’abaque la valeur correspondante : 0,3c° (Annexe 3.15)

 Calcule de la nouvelle valeur de Tsa

Te – Tsa =0,3c°

Tsa= Te –0,3c°

Tsa=0,8-0;3

Tsa=0,5c°

 Remarque :

Puisque la position de l’évaporateur est supérieure à celle du groupe de

condensation, nous n’avons pas besoin d’installer un siphon dans la ligne
d’aspiration, car le retour d’huile est assuré par gravité.
5.5- Calcul de la masse du fluide frigorigène :

On calcule la masse du fluide frigorigène d’une installation en la supposant en

fonctionnement et en régime établi
- Caractéristiques du fluide frigorigène R404A :
Nature : Hydrofluorocarbure
Mélange : Zéothrope
Composition : 52% R143 / 44% R125 / 4% R134a
Glissement 0,9 °C à 0,3 bar
Couleur Incolore
Odeur Ethérée légèrement
Température + 72 °C
critique :
Température - 45,8 °C (pression atmosphérique)
ébullition :
Pression critique : 37,4 Bar

- La charge en fluide frigorigène :

La bonne charge en fluide frigorigène est très importante pour le bon
fonctionnement d’une installation frigorifique.
S’il y a un manque de charge, l’alimentation en liquide à la sortie du
détendeur va être insuffisante et cela provoquera une charge excessive à
l’aspiration du compresseur, il va s’ensuivre une baisse de la puissance
frigorifique, des risque de mauvais retour d’huile, une élévation de la
température de refoulement excessive.
Par contre, si on a une charge trop importante, le liquide va prendre un
volume excessif dans le condenseur et provoque une élévation de la pression de
condensation; cela conduit à un mauvais refroidissement, une baisse de la
puissance frigorifique, une augmentation de la puissance absorbée par le
compresseur et une température de refoulement anormale pouvant endommager
le compresseur.
 En fonctionnement, les organes de l’installation frigorifique contiennent un
certain pourcentage de fluide frigorigène (R 404a).
1-25 du volume intérieur de l’évaporateur.
2-25 du volume intérieur de condenseur.
3-10 a15 du volume intérieur du réservoir.
4-Canalisations.

1- Pour évaporateur :

On calcule la masse du fluide pour évaporateur la fourmille suivant :

M=V*F*𝑝

M : masse de fluide en kg

V : volume intérieur en dm3(Annexe3.16)

F : facteur de remplissage 25 pour l’évaporateur à détente séche

𝑝 :la masse volumique de fluide en kg / dm3 (Annexe3.12)

 interpolation :

𝑝= y1+(𝑦2−𝑦1) ×(x-x1)
(𝑥2−𝑥1)

𝑝 = 1150+(1146,2−1150) ×(0.8-0)
(1−0)

𝑝 = 1150-3,04

𝑝 = 1146;96 kg / 𝑚3 =1;146kg / 𝑑𝑚3

Mev = 9,84*0,25*1,146

Mev=2,81kg

2- Pour condenseur :

V : volume intérieur en dm3se 12;5dm3(Annexe3.17)

F : facteur de remplissage 25 pour condenseur à détente séche

𝑝 : Masse volumique de liquide saturé à 25c° (hiver) (Annexe3.12)

𝑝=1044 kg/m3𝑝= 1,044 kg/dm3

Mcd = 12.5*0,25*1,044

Mcd=3.26kg

3- Réservoir haute pression :

V : Volume intérieur du réservoir dm3 se 10 dm3 (Annexe3.18)

𝑝 : Masse volumique du liquide saturé à 15

𝑝=1044 kg/m3𝑝= 1,044 kg/dm3 (Annexe3.12)

Mr= 10.0,10.1,044

Mr=1.044kg

4- Canalisation :

M= L .p. R

M : masse de fluide en kg

L : longueur de la conduite en(m)

: masse du fluide par longueur en dm3

R : remplie de liquide et vapeur 100

4.1- La conduite liquide :

θliq : 1/2 ” = 0,093dm3(Annexe3.19)

Lliq : longueur de la conduite liquide : 8m

: Densité de la masse du haute pression :1;044 kg/dm3 (Annexe3.12) R

: remplie de liquide 100

Mcon-liq= 8.0,093.1.1,044

Mcon-liq= 0,77kg

4.2-La conduite d’aspiration :

Lliq : longueur de la conduite liquide : 8m

θasp : 1 3/8 ” = 0;81dm3

:Densité de la masse du basse pression : 0.81dm3

R : remplie de vapeur 100

 Interpolation :

𝑝 = y1+(y2−y1) ×(x-x1)
(X2−X1)

𝑝 = 30,457+(31,446−30,457) ×(0.8-0)
(1−0)

𝑝 = 30,457+0,7912

𝑝= 31,2482kg / m3 =0,031kg / dm3

Mcon-asp = 8.0,81.0,031.1

Mcon-asp =0,20 kg

4.3-La conduite de refoulement :

elle est négligeable ;il s’agit d’un groupe avec longueur de tube de refoulement
inférieure à 1m le plus souvent.

 calcule de la masse totale :

Mtotale=Mev+Mcd+MR+Mcon-liq+ Mcon-asp

Mtotale=2,81 +3.26+1.044 +0,77 +0,20

Mtotale=8.08kg

 Vérification :

Il faut s’assurer que le réservoir (Hp) peut contenir la totalité de la masse à

injecter 90du recevoir (Hp) représente :

Mmax= 10.0,9.1 ,044

Mmax=9,39kg

Ce qui est parfait car : MTotal<Mmax.