CHAPITRE 4 ETUDE MECANIQUE ET HYDRODYNAMIQUE

4.1 Calcul mécanique de vérification

4.1.1 Introduction

Cette partie du chapitre est consacré aux calculs de vérification de la colonne débutaniseur :

- La résistance de la virole aux contraintes dues à la pression de chaque section de la

colonne.
- La résistance des fonds aux contraintes dues à la pression.
- La résistance du fond inférieur lors de l’essai hydrostatique.
- Le calcul des brides.

Ces calculs se font suivant les conditions de service et les dimensions données par le
constructeur, présentées dans les tableaux suivants :

Tableau 4.1. Matériaux de constructions :

Pièces Nom du matériau La limite d’élasticité Limite de la résistance

(MPa) (MPa)
Virole A 516 GR 70 260 485
Fonds A 516 GR 70 260 485
Boulons JIS G 3454 205 450

Tableau 4.2.Paramètres de fonctionnement :

Fluide Hydrocarbure
Pression de service 16.5 bar
Température de servie 280 °c
Pression d’épreuve hydrostatique 28.8 bar

34
CHAPITRE 4 ETUDE MECANIQUE ET HYDRODYNAMIQUE

Tableau 4.3. Données géométrique de la colonne

Le diamètre intérieur de la virole 3 000 mm

L’épaisseur de la virole 24 mm
L’épaisseur du fond inferieur 28 mm
L’épaisseur du fond supérieur 28 mm
La hauteur du fond inferieur 450 mm
La hauteur du fond supérieur 800 mm
La hauteur de la jupe 6500 mm
Le diamètre intérieur de la jupe 3030 mm
L’épaisseur de la jupe 12 mm
Le diamètre extérieur de la jupe 3054 mm
Le diamètre extérieur de l’anneau de base 3330 mm
Le diamètre intérieur de l’anneau de base 2770 mm
La largeur de l’anneau de base 250 mm
L’épaisseur de l’anneau de base 560 mm
La hauteur de la colonne + la jupe 35750 mm
Le poids de la colonne à vide 72500 kg
Le poids de la colonne en service 147000 kg
Le poids de la colonne lors de l’essai hydraulique 317000 kg
La surépaisseur de la corrosion 3,175 mm
Le coefficient de soudure de la virole 1
Coefficient de soudure des fonds 1
Le calorifugeage 50 mm

35
CHAPITRE 4 ETUDE MECANIQUE ET HYDRODYNAMIQUE

Figure 4.1 : Représentation de la colonne

36
CHAPITRE 4 ETUDE MECANIQUE ET HYDRODYNAMIQUE

4.1.2 Résistance de la virole dans chaque section de la colonne

𝑒
Pour le calcul des appareils à paroi mince (𝐷𝑟 < 0,1), on utilise l’équation de LAPLACE
𝑖

[11] :

𝑃 𝜎𝑍 𝜎𝜏 (4.1)
= +
𝑒𝑟 𝜌𝑍 𝜌𝜏
Avec :

P : pression interne

𝑒𝑟 : Épaisseur de la résistance.

𝜎𝑍 : Contrainte axiale parallèle à Z.

𝜎𝜏 : Contrainte tangentielle.

𝜌𝑍 : Rayon de courbure de la direction de 𝜎𝑍 : 𝜌𝑍 = 𝑅𝑚

𝜌𝜏 : Rayon de courbure de la direction de 𝜎𝜏 : 𝜎𝜏 = ∞

Comme dans le cas d’une configuration cylindrique 𝜎𝜏 = ∞ alors le seconds terme de

l’équation de LAPLACE s’annule ce qui nous permet d’écrire :

Dans les différentes sections de la colonne, la condition suivante dois être vérifiée

𝑃𝐶 × 𝑅𝑚 (4.1.a)
𝜎𝑍 = ≤ [𝜎] × 𝜑
𝑒𝑟
𝜑 : Coefficient de diminution de la résistance du au soudage (donnée)

[𝜎] : Caractéristique du matériau, fonction de la limite d’élasticité (𝜎𝑒 ) et de rupture (𝜎𝑟 )

𝜎𝑒𝑡𝑐 𝜎𝑟𝑡𝑐 (4.1.b)

[𝜎] = 𝑚𝑖𝑛 { ; }
𝑛𝑒 𝑛𝑟

𝑛𝑒 , 𝑛𝑟 : Coefficients de sécurité : 𝑛𝑒 = 1,5 et 𝑛𝑟 = 3 (norme AFNOR)

Les contraintes :

 𝜎𝑒𝑡𝑐 = 𝐾 × 𝜎𝑒20°𝑐
 Si 20°𝐶 ≤ 𝑡𝑐 ≤ 350°𝐶 alors 𝜎𝑟𝑡𝑐 = 𝜎𝑒20°𝐶
 Si 𝑡𝑐 ≥ 350°𝐶 alors 𝜎𝑟𝑡𝑐 = 𝐾 × 𝜎𝑒20°𝐶

37
CHAPITRE 4 ETUDE MECANIQUE ET HYDRODYNAMIQUE

Avec :

𝑡𝑐 : Température de calcul

K : coefficient tenant compte de l’influence de la température sur les caractéristiques

mécaniques du matériau, déterminé graphiquement

Figure 4.2 : coefficient K en fonction de la température

a. Section (a) partie inférieur :

𝑒𝑟 = 𝑒𝑣 − (𝑒1 + 𝑒2 + 𝑒3 ) (4.1.c)

Avec :

𝑒𝑟 : Épaisseur de la résistance de la virole.

𝑒𝑣 : Épaisseur réelle de la virole, donnée

𝑒1 : Surépaisseur de la corrosion, donnée

𝑒2 : Surépaisseur de tolérance négative, donnée

𝑒3 : Surépaisseur tenant compte des charges supplémentaires (vent, séisme…) pour les
équipements verticaux, choisi à partir du tableau

38
CHAPITRE 4 ETUDE MECANIQUE ET HYDRODYNAMIQUE

Tableau 4.4. Choix de 𝑒3 en fonction de la hauteur

Hauteur (m) 𝑒3 (mm)

< 20 1÷2
20 ÷ 40 3÷4
> 40 4÷6

Comme la hauteur de la section (a) est H= 8330 mm alors 𝑒3 = 2

En remplaçant ces épaisseurs dans la relation (4.1.c), on obtient :

𝑒𝑟 = 19,025 𝑚𝑚

Le rayon moyen est donné par la relation suivante :

𝐷𝑖 +𝑒𝑟𝑣 (4.1.d)
𝑅𝑚 = = 1509,01𝑚𝑚
2
La pression de calcul (Pc)

𝑃𝑐 = 𝑚𝑎𝑥{𝑃𝑠 + 1𝑏𝑎𝑟; 𝑃𝑠 + 10%𝑃𝑠 } = 1,60 𝑀𝑃a (4.1.e)

Avec : Ps : Pression de service

La température de calcul (Tc)

𝑇𝑐 = 𝑇𝑠 + 15°𝑐 = 224°𝑐 (4.1.f)

Avec : 𝑇𝑠 : Température de service

La contrainte admissible

[𝜎] = 138,6 𝑀𝑃𝑎

Tel que: K= 0,8

Le résultat :

𝜎𝑍 = 127,38 𝑀𝑃𝑎 ≤ [𝜎] × 𝜑 = 138,6 MPa

Alors la résistance de cette partie est vérifiée.

39
CHAPITRE 4 ETUDE MECANIQUE ET HYDRODYNAMIQUE

a. Section (b) partie centre de la colonne :

De la même manière se quand on fait dans la section inferieur (a):

On obtient

 𝑒𝑟 = 19,025 𝑚𝑚
 𝑅𝑚 = 1509,01𝑚𝑚
 𝑃𝑐 = 1,595 𝑀𝑃a
 𝑇𝑐 = 167°𝑐
 [𝜎] = 147,33 𝑀𝑃𝑎

Le résultat :

𝜎𝑍 = 126,51 𝑀𝑃𝑎 ≤ [𝜎] × 𝜑 = 147,33 MPa

Alors la résistance de cette partie est vérifiée

a. Section (c) partie supérieur de la colonne :

De la même manière

 𝑒𝑟 = 18,025 𝑚𝑚
 𝑅𝑚 = 1509,01 𝑚𝑚
 𝑃𝑐 = 1,54 𝑀𝑃a
 𝑇𝑐 = 98°𝑐
 [𝜎] = 161,6 𝑀𝑃𝑎

Le résultat :

𝜎𝑍 = 130,12 𝑀𝑃𝑎 ≤ [𝜎] × 𝜑 = 161,6 MPa

Alors la résistance de cette parie est vérifiée

40
CHAPITRE 4 ETUDE MECANIQUE ET HYDRODYNAMIQUE

4.1.3 Vérification de la résistance des fonds de la colonne

La colonne dé-butaniseur (C102) à des fonds de forme elliptique

Figure 4.3. Représentation d’un fond elliptique

La condition de la résistance des fonds est [11,14]:

𝑌 × 𝑃𝑐 × 𝑅𝑚 (4.2)
𝜎𝜏 = ≤ [𝜎] × 𝜑
𝑒𝑟
Le coefficient de diminution de la résistance du au soudage pour les fonds 𝜑 = 1
𝐷
𝑌 : Le facteur de forme de fond dépend du rapport( 𝐻𝑖 ), sa valeur est donnée par la formule

suivante :
1 𝐷𝑖 2 (4.2.a)
𝑌= [2 + ( ) ]
6 2𝐻𝐹
𝐻
Si 𝑌 = 1 → = 0,25
𝐷𝑖

4.1.3.1 La résistance du fond inférieur

[𝜎] et 𝑃𝑐 : ont les mêmes valeurs de la section (a).

 La Hauteur 𝐻𝐹 𝑖𝑛𝑓 = 800 mm

 𝑌 =0,9

𝑒𝑟 𝐹 𝑖𝑛𝑓 = 𝑒𝑓 𝐹 𝑖𝑛𝑓 − (𝑒1 + 𝑒2 + 𝑒3 ) =18,42 mm (4.2.b)

Avec :

𝑒𝑟 𝐹 𝑖𝑛𝑓 : Épaisseur de la résistance du fond inférieur.

𝑒𝐹 𝑖𝑛𝑓 : Épaisseur du fond inférieur (𝑒𝐹 𝑖𝑛𝑓 = 28 𝑚𝑚)

 𝑅𝑚 =1509,01 mm

41
CHAPITRE 4 ETUDE MECANIQUE ET HYDRODYNAMIQUE

Le résultat :
𝜎𝜏 = 118,38 𝑀𝑃𝑎 ≤ [𝜎] × 𝜑 = 138,66 𝑀𝑃𝑎
On remarque bien que la résistance du fond inférieur est vérifiée.

4.1.3.2 La résistance du fond supérieur

[𝜎] et 𝑃𝑐 : ont le même valeurs de la section (b)

 La hauteur 𝐻𝐹 𝑠𝑢𝑝 = 800 mm

 𝑌 =0,9

𝑒𝑟 𝐹 𝑠𝑢𝑝 = 𝑒𝐹 𝑠𝑢𝑝 − (𝑒1 + 𝑒2 + 𝑒3 ) =17,42 mm (4.2.c)

Avec :
𝑒𝑟 𝐹 𝑠𝑢𝑝 : Épaisseur de la résistance du fond supérieur.

𝑒𝐹 𝑠𝑢𝑝 : Épaisseur du fond supérieur (𝑒𝐹 𝑠𝑢𝑝 = 28 𝑚𝑚)

 𝑅𝑚 =1509,01 mm

Le résultat :

𝜎𝜏 = 121,14 𝑀𝑃𝑎 ≤ [𝜎] × 𝜑 = 156 𝑀𝑃𝑎

On remarque bien que la résistance du fond supérieur est vérifiée.

4.1.4 La résistance du fond inférieur lors de l’essai hydrostatique

La condition de résistance du fond lors de l’essai hydrostatique est la suivante [9] :

𝑃𝑓 𝑖𝑛𝑓 × 𝐷𝑚 (4.3)
𝜎𝑓𝑜𝑛𝑑 = ≤ 0,9 × 𝜎𝑒𝑡
2 × 𝑒𝑓 𝐹 𝑖𝑛𝑓 × 𝜑
Avec :

𝐷𝑚 : Le diamètre moyen du fond.

𝐷𝑚 = 𝐷𝑖𝑛𝑓 + 𝑒𝑟 = 3023,02 𝑚𝑚 (4.3.a)

𝜎𝑒𝑡 : Limite d’élasticité à la température du service (caractéristique mécanique du matériau du

fond inférieur).

𝑃𝑓 𝑖𝑛𝑓 : La pression du fond inférieur.

𝑃𝑓 𝑖𝑛𝑓 = 𝑃𝑠𝑡 + 𝑃𝑒𝑠𝑠𝑎𝑖 = 3,11 𝑀𝑃𝑎 (4.3.b)

42
CHAPITRE 4 ETUDE MECANIQUE ET HYDRODYNAMIQUE

Avec :

𝑃𝑒𝑠𝑠𝑎𝑖 : Pression d’essai hydrostatique, (pour cette colonne𝑃𝑒𝑠𝑠𝑎𝑖 = 2,83 𝑀𝑃𝑎)

𝑃𝑠𝑡 : Pression de l’eau contenue dans la colonne.

Pst = ρgH = 0,3 𝑀𝑃𝑎 (4.3.c)

Avec :

𝐻 = 𝐻𝑣𝑖𝑟𝑜𝑙𝑒 + 𝐻𝑓𝑜𝑛𝑑 𝑖𝑛𝑓 + 𝐻𝑓𝑜𝑛𝑑 𝑠𝑢𝑝 =29650 mm (4.3.d)

D’où :

0,9 × 𝜎𝑒𝑡 =234 MPa (4.3.e)

𝜎𝑒𝑡 =260 MPa

Le résultat :

𝜎𝑓𝑜𝑛𝑑 = 204,16 ≤ 0,9 × 𝜎𝑒𝑡 = 234 𝑀𝑃𝑎

On remarque que le fond inférieur résiste à l’effort de pression lors de l’essai hydrostatique.

4.1.5 Calcule des brides :

Les brides doivent assurer la résistance l’étanchéité nécessaire contre la charge de

fonctionnement, les types des brides utilisées dans cette colonne sont «Slip-on pour les trous
d’homme et Welding-Neck pour les autre raccords ». On prend comme exemple de calcule la
bride de trou d’homme type« Slip-on» [12].

Figure 4.4 : représentation de la bride Slip-On

43
CHAPITRE 4 ETUDE MECANIQUE ET HYDRODYNAMIQUE

Le calcul se fait en deux parties :

 Calcul des boulons ou goujons : boulons pour P ≤ 60bars et goujons P ≥ 60bars.

4.1.5.1 Vérification de l’étanchéité en fonction du pas (t) [11]:

𝜋 × 𝐷𝑏 (4.4)
𝑡=
𝑛𝑏
Avec :

𝑡: Le pas

𝐷𝑏 : Le diamètre de passage de la bride 𝐷𝑏 = 812 mm

𝑛𝑏 : Le nombre de boulons 𝑛𝑏 = 24

Condition d’étanchéité :

𝑡 ≤ (4 ÷ 5) × 𝑑𝑏 𝑠𝑖 𝑃 < 25 𝑏𝑎𝑟 (4.4.a)

{
𝑡 ≤ (2.5 ÷ 3) × 𝑑𝑏 𝑠𝑖 𝑃 ≥ 25 𝑏𝑎𝑟

𝑑𝑏 : Diamètre du boulon 𝑑𝑏 = 41.1𝑚𝑚

𝑃𝑐 : Pression de calcul 𝑃𝑐 = 15,75 bar

On a : 𝑡 = 106.23 mm < 4 × 41,1 = 164.4 mm donc t < (4÷5)× 𝑑𝑏

et on a 𝑃𝑐 = 15,75 bar < 25 𝑏𝑎𝑟

On conclue que la condition d’étanchéité est bien vérifie

4.1.5.2 Condition de résistance des boulons à la traction :

𝑞𝑏 (4.5)
𝜎𝑡𝑟 = ≤ [𝜎]
𝑆𝑏
𝑆𝑏 : La surface du boulon

𝜋𝑑𝑏2 (4.5.a)
𝑆𝑏 = = 1326,02 mm2
4

44
CHAPITRE 4 ETUDE MECANIQUE ET HYDRODYNAMIQUE

 [𝜎] : contrainte admissible à la traction à la température au fond de la colonne tc= 224°c

𝜎𝑒224°𝐶 𝜎𝑟224°𝐶 (4.5.b)
[𝜎] = 𝑚𝑖𝑛 ( , )
𝑛𝑒 𝑛𝑟
D’après le graph (Figure 4.2) on détermine k :
𝜎𝑒224°𝐶 = 𝒌 × 𝜎𝑒20°𝐶 = 0,8 × 205 = 164 𝑀𝑃𝑎 (4.5.c)

𝜎𝑟224°𝐶 = 𝜎𝑟20°𝐶 = 450 𝑀𝑃𝑎 (4.5.d)

Tel que (𝑛𝑒 = 1.5; 𝑛𝑟 = 3)

[𝜎] = 109,33 𝑀𝑃𝑎

𝑞𝑏 : La charge sur un boulon :

𝑄𝑏 (4.5.e)
𝑞𝑏 =
𝑛𝑏

 Calcul de la charge agissant sur les boulons (𝑄𝑏 ) :

𝑄𝑏 = max(𝑄𝑏1 , 𝑄𝑏2 ) (4.5.f)

𝑄𝑏1 : Charge lors d’un fonctionnement normal (force de pression+ force de serrage).

𝑄𝑏2 : Charge due uniquement à de la force de serrage.

a. Calcule de 𝑸𝒃𝟏 :

2
𝜋 × 𝐷𝑚 (4.5.g)
𝑄𝑏1 = (𝛼 × 𝑃 × ) + (𝜋 × 𝐷𝑚 × 𝑏0 × 𝑃 × 𝑚)
4

𝛼 : Coefficient de rigidité de la liaison de la bride il dépend de type et de matériau de joint ;

pour les joints métallique on a : 𝛼= 1.45.

𝑏0 : Largeur efficace de joint pour notre cas :

𝑏0 = 1,2√𝑏 (4.5.h)

45
CHAPITRE 4 ETUDE MECANIQUE ET HYDRODYNAMIQUE

D’où :

𝑏 = 𝐷𝐴𝐵 − 𝐷𝑖 = 812,8 − 616 = 196,8 𝑚𝑚 (4.5.i)

𝑏0 = 16,83 𝑚𝑚

𝑚 : Coefficient d’étanchéité il dépend du matériau utilisé et de type de joint, pour notre cas :
𝑚=3

𝑃: Pression de calcul 𝑃 = 15,75 bar

𝐷𝑚 : Diamètre moyen.

𝐷𝐴𝐵 − 𝐷𝑖 (4.5.j)
𝐷𝑚 = 𝐷𝑖 + = 654,1 𝑚𝑚
2

𝑛𝑏 : Nombre des boulons 𝑛𝑏 = 24

𝑄𝑏1 = 0,930. 106 𝑁

b. Calcul de 𝑸𝒃𝟐 :

𝑸𝒃𝟐 = 𝝅 × 𝑫𝒎 × 𝒒𝒔 × 𝒌 × 𝒃𝟎 (4.5.k)
𝑞𝑠 : La charge spécifique pour le joint de notre cas 𝑞𝑠 = 40 N/mm²

𝒌 : Coefficient obtenu a partir du graphe de figure (4.2) 𝑘 = 0.8

𝑄𝑏2 = 1,106 𝑁
Donc

𝑄𝑏 = max( 𝑄𝑏1 , 𝑄𝑏2 ) = 1,106 𝑁 (4.5.l)

La charge agissant sur un boulon :

𝑄𝑏
𝑞𝑏 = = 46088,93 𝑁
𝑛𝑏
Alors

𝜎𝑡𝑟 = 34,75 𝑀𝑃𝑎 < [𝜎] = 109,33 𝑀𝑃𝑎

On conclue que la résistance des boulons à la traction est vérifiée

46
CHAPITRE 4 ETUDE MECANIQUE ET HYDRODYNAMIQUE

4.2 Calcul hydrodynamique

4.2.1 Introduction

Cette partie est consacrée au calcul reliant la conception des éléments de contact (les
plateaux) dans la colonne en fonctionnement hydrodynamique optimal (vitesse maximale et
perte de charge minimale).

Pour notre cas, nous avons choisis le 8ème plateau pour faire le calcul

Les données de départ nécessaire pour le calcul sont représentées dans les tableaux
suivants :

Tableau 4.5. Dimensions du plateau

Dimensions du plateau Unités Valeurs

Diamètre mm 3000
Type de plateau - A clapet (V-1)
Nombre de passe - 2
Longueur du déversoir mm 2000
Largeur du déversoir mm 2200
Hauteur du déversoir mm 710
Distance entre deux plateaux mm 610
Nombre de clapets - 350
Tableau 4.6. Paramètres de fonctionnement

Vapeur sur le plateau

Température de vapeur C° 83
Débit massique kg/h 154800
Densité kg/m3 36,7
Débit volumique m3/h 4217,98
Liquide sur le plateau
Température de liquide C° 78
Débit massique kg/h 176760
Densité kg/m3 494
Débit volumique m3/h 357,81
Viscosité Pc 0,16

47
CHAPITRE 4 ETUDE MECANIQUE ET HYDRODYNAMIQUE

Figure 4.5 : Schématisation de l’équilibre hydrodynamique d’un plateau

4.2.2 Facteur capacité :

Ce facteur dépend des caractéristiques physico-chimiques des phases liquides et vapeur et de

la géométrie de la colonne, pour le déterminer on utilise la corrélation de FAIR en sachant
l’espace entre deux plateaux et en utilisant la formule suivante [9,11]:

𝐿 𝜌𝑣
√ = 0,6 (4.6)
𝐺 𝜌𝑙

L, G : débits massiques du liquide et vapeur respectivement

ρl , ρv : Masses volumiques du liquide et vapeur respectivement

48
CHAPITRE 4 ETUDE MECANIQUE ET HYDRODYNAMIQUE

Figure 4.6 : Diagramme de corrélation de FAIR

Alors : C = 0,073

4.2.3 Vitesse de vapeur optimale

Calculée selon la relation suivante [4,10] :

𝜌𝐿 − 𝜌𝑉
𝑈𝑠𝑡 = 𝑐√ = 0,27 m/s
𝜌𝑉 (4.7)

La vitesse maximale est obtenue à 80% de l’engorgement en tenant compte du facteur de

moussage 0.9

𝑈𝑚𝑎𝑥 = 0.8 × 0.9 × 𝑈𝑠𝑡 = 0,19 m/s (4.8)

4.2.4 Facteur perte de charge

Il est donné par la formule suivante [10] :

𝐹𝑎 = 𝑈𝑚𝑎𝑥 × (𝜌𝑣 )0.5 = 1,15 kg/𝑚2 𝑠 (4.9)

4.2.5 Efficacité

Elle repère la qualité du contacte liquide vapeur qui conditionne le fonctionnement, calculé
par la formule suivante [9,11] :

𝐸 = 0,17 − 0,616. log(𝜇) = 0,66 (4.10)

Avec : 𝜇 : Viscosité cinématique

49
CHAPITRE 4 ETUDE MECANIQUE ET HYDRODYNAMIQUE

4.2.6 Trajet du liquide sur le plateau (TLP)

Il amène à la détermination du nombre d’orifices sur le plateau

𝐷𝑝 𝐻𝑖 (4.11)
𝑇𝐿𝑃 = (1 + ∑ ) = 1700 𝑚𝑚
𝑁𝑝 𝐷𝑝

Avec : 𝐷𝑝 Diamètre du plateau

𝑁𝑝 Nombre de passe

𝐻𝑖 Largeur de chaque goulotte

𝐷𝑝 − 𝐿𝑑 (4.11.a)
𝐻𝑖 = = 400 𝑚𝑚
2
𝐿𝑑 Largeur du déversoir

4.2.7 L’espace libre d’échange

𝜋 × 𝐷2 (4.12)
𝐴𝑛 = 0,85 × = 6 𝑚2
4
Avec : D : Diamètre de la colonne

4.2.8 Temps de séjour dans les décentes

A sa chute, le liquide passe par le déversoir dont il reste un temps avant d’entrer au plateau
inférieur, ce temps de séjour est calculé comme suit :

𝑆𝑑 × 𝑇 (4.13)
𝑇𝑠 = = 8,72 𝑠
𝑄𝑙
Avec : 𝑆𝑑 La section du déversoir

𝑆𝑑 = 𝑙𝑑 × ℎ𝑑 = 1,42 𝑚2 (4.13.a)

𝑙𝑑 Longueur du déversoir ( 𝑙𝑑 = 2 𝑚)

ℎ𝑑 Hauteur du déversoir (ℎ𝑑 = 710 𝑚𝑚)

𝑇 La distance entre deux plateaux (𝑇 = 610 𝑚𝑚)

𝑄𝑙 La charge du liquide sur le plateau

50
CHAPITRE 4 ETUDE MECANIQUE ET HYDRODYNAMIQUE

4.2.9 Calcul des pertes de charge

4.2.9.1 Hauteur du liquide au-dessus du déversoir

Définie par l’équation hydraulique :

𝑄𝐿 2/3 (4.14)
𝐻𝑜𝑑 = 0,15 ( ) = 1.9 𝑐𝑚
𝐿𝑑
Avec : 𝑄𝐿 Débit volumique du liquide

𝐿𝑑 Largeur du déversoir

4.2.9.2 Hauteur du liquide à la sortie du déversoir

Cette hauteur est définie par :

ℎ𝑔𝑙 (4.15)
𝐻𝑙 = 𝐻𝑜𝑑 + ℎ + = 15,9𝑐𝑚
2
Avec : h Hauteur de barrage déversoir

ℎ𝑔𝑙 Gradient liquide qui est introduit seulement à l’entrée du liquide sur le plateau,
pour notre cas, nous considérons la hauteur à la sortie du plateau donc ℎ𝑔𝑙 = 0

4.2.9.3 Hauteur du liquide clair

C’est le produit de la hauteur du liquide à la sortie du déversoir par le facteur 𝛽

𝐻𝑐𝑙 = 𝐻𝑙 × 𝛽 = 8,6 𝑐𝑚 (4.16)

Avec : 𝛽 : Facteur d’aération

𝑎 + 𝑏 × 𝑈𝑔 × √𝜌𝑣 (4.16.a)
𝛽= = 0.54
𝑎 + 𝑐 × 𝑈𝑔 × √𝜌𝑣

Pour les plateaux à clapets : 𝑎 = 1,063 ; 𝑏 = 2,002 ; 𝑐 = 3,94

𝑈𝑔 = 0,6 En fonctionnement normal

4.2.9.4 Perte de charge au niveau des orifices

La perte de charge des plateaux à clapets est définie par la formule suivante :

𝜌𝑣 (4.17)
𝐻0 = (𝐾1 + 𝐾2 ) × × 𝑈0 2 = 1.35𝑐𝑚
𝜌𝐿

51
CHAPITRE 4 ETUDE MECANIQUE ET HYDRODYNAMIQUE

Avec :

𝐾1 = 0

𝐾2 : Coefficient dépendent du type de clapet et de son épaisseur

Remarque : Type de clapet V-1 et e= 6 mm donc 𝐾2 =16

𝑈0 Vitesse de la vapeur dans l’orifice

𝑄𝑣 (4.17.a)
𝑈0 = = 0,1𝑚/𝑠
𝑁𝑐 𝑆0
𝑄𝑣 : Débit volumique de la vapeur sur le plateau [m3/s]

𝑁𝑐 : Nombre de clapet sur le plateau : 𝑁𝑐 =350

𝑆0 : Section de l’orifice :

𝑆0 = 𝜋. 𝐷𝑜 . ℎ𝑐 = 3 ,14 × 0.1 × 0.1 = 0.0314𝑚2 (4.17.b)

𝐷𝑜 : Diamètre de l’orifice

ℎ𝑐 : Hauteur du clapet

4.2.9.5 Perte de charge dans la section de descente

Représente la perte de charge du liquide sous la chicane de descente à l’entrée du plateau

inférieur, elle se calcule à partir de l’équation suivante :

𝑄𝑙 2 (4.18)
𝐻𝑠 = 16,5 × ( ) = 2,82𝑐𝑚
𝐴𝑠
𝐴𝑠 : Section de passage offerte au liquide entre la base de la tôle de descente et le plateau
inférieur

𝐴𝑠 = 𝑙𝑑 × 𝑒 = 0,240 𝑚2 (4.18.a)

𝑒 : L’espace entre le déversoir et le plateau inférieur (𝑒 = 12 )𝑐𝑚

4.2.9.6 Perte de charge par plateau

La perte de charge par plateau est donnée par la formule suivante [ref] :

𝐻𝑝 = 𝐻0 + 𝐻𝑐𝑙 = 9,94 𝑐𝑚 (4.19)

52
CHAPITRE 4 ETUDE MECANIQUE ET HYDRODYNAMIQUE

4.2.9.7 Bilan pression autour du plateau (𝑯𝒅 )

𝐻𝑑 est la hauteur du liquide clair dans la goulotte de descente, ce bilan est calculé comme
suit :

𝐻𝑑 = 𝐻𝑠 + 𝐻𝑙 + 𝐻𝑝 = 28,66 𝑐𝑚 (4.20)
Pour un fonctionnement normal, la valeur de 𝐻𝑑 ne doit pas dépasser 50% de
l’espacement entre deux plateaux (𝑇) c’est-à-dire [7,11] :

𝐻𝑑 < 0,5 × 𝑇 (4.20.a)

On a :

0,5 × 𝑇 = 30,5 𝑐𝑚
Alors :

28,66 𝑐𝑚 < 30,5 𝑐𝑚 → 𝐻𝑑 < 0,5 × 𝑇

 Conclusion :

Comme les conditions concernant les pertes de charge sont toutes vérifiées, le système
fonctionne normalement.

53
CHAPITRE 4 ETUDE MECANIQUE ET HYDRODYNAMIQUE

4.3 Simulation par solidworks

Solidworks est un logiciel de conception assistée par ordinateur, appartenant à la société

Dassault Systèmes. Il utilise le principe de conception paramétrique et génère trois types de
fichiers qui sont liés: la pièce, l'assemblage, et la mise en plan

4.3.1 Modélisation de la colonne

Nous avons utilisé l’outil solidworks pour modéliser la conception de notre colonne étudie,
selon les étapes suivantes :

 Esquisser les plans initiaux de la colonne (les dimensions de la virole, les fonds et la
jupe)
 Faire des fonctions pour maintenir la colonne en forme 3D (bossage, enlèvement de la
matière, révolution)

Figure 4.7 : Quelques fonctions utilisées dans la modélisation

 Assembler les plateaux et les brides

54
CHAPITRE 4 ETUDE MECANIQUE ET HYDRODYNAMIQUE

Figure 4.8 : la conception de la colonne après l’assemblage (3D)

4.3.2 Vérifications des quelques éléments constitutifs de la colonne

Pour commencer l’étude de vérification il faut préciser le matériau (A516 GR 70), choisir
l’étude statique et mailler la colonne (triangulaire 500mm)

Figure 4.9 : maillage de la colonne

55
CHAPITRE 4 ETUDE MECANIQUE ET HYDRODYNAMIQUE

a. La résistance de la virole

Figure 4.10 : les contraintes de Von Mises sur la virole

On applique la pression P=1,60 MPa sur les parois intérieurs de la colonne on obtient la
contrainte maximale (122,30 MPa) est inférieur à la contrainte admissible (138,6 MPa) donc
la résistance de la virole est vérifiée.

c. La résistance des fonds de la colonne

 Fond supérieur à P = 1,55 MPa

Figure 4.11 : les contraintes de Von Mises au fond superieur

56
CHAPITRE 4 ETUDE MECANIQUE ET HYDRODYNAMIQUE

On remarque que la contrainte maximale (119,40 MPa) est inférieur à la contrainte admissible
(156 MPa) donc la résistance du fond supérieur est vérifiée.

 Fond inferieur lors de l’essai hydrostatique P = 3,11 MPa

Figure 4.12 : les contraintes de Von Mises sur la virole

On remarque que la contrainte maximale (202,88 MPa) est inférieur à la contrainte admissible
(234 MPa) donc la résistance du fond inferieur lors de l’assai hydrostatique est bien vérifiée.

4.4 Comparaison entre les résultats dans le cas calculé et le cas simulé

Tableau 4.7 Comparaison des deux résultats obtenus

contrainte Cas calculé Cas simulé Ecart

maximale (MPa) (SolidWorks) (Mpa) [%]
La virole 130,12 122,30 6
Fond supérieur 121,14 119,40 1,4
Fond inférieur lors de 204,16 202,88 0,6
l’essai hydrostatique

Vue que l’écart entre ces deux méthodes ne dépasse pas 6% alors on peut dire que les
valeurs sont bien proches alors ces deux méthodes sont vérifiées pour ce calcul.

57