10 juillet 1905 /2019

Mini projet : Silos et Réservoirs

Encadré par : Moez Trabelsi

Elaboré par : Firas Jaoua

Ahmed Mechim
 REMERCIEMENTS
Nous exprimons nos vifs remerciements à Mr MOEZ TRABELSI pour
tout l’aide qu’il nous a fourni, nous sommes lui reconnaissants pour
ses conseils et sa disponibilité.
Il est d'une importance primordiale pour l'ingénieur en génie civil
d’avoir des connaissances. Ce projet nous a donné la possibilité de
compléter notre bagage par des acquis d'ordre techniques et
pratiques qui constituent les premiers pas dans le métier d’ingénieur.
C'est ainsi que nous avons pu constater que la conception est une
étape de travail importante pour l'ingénieur. Le travail réside dans le
choix judicieux des démarches à suivre durant cette étude, et ses
difficultés de réalisation.

P a g e 1 | 28
 INTRODUCTION
Dans ce projet on va s’intéresser à étudier des concepts
importants pour la formation de l’ingénieur.
Pour cela durant la période de ce bureau d’étude on va
s’intéresser à :

Le dimensionnement d’un réservoir circulaire.

L’évaluation des charges sur la couverture.
Le dimensionnement de la coupole.
Le dimensionnement de la ceinture.
Le dimensionnement de la paroi.

P a g e 2 | 28
 Données du projet

Figure Donnée du projet

Données :

Fissuration très préjudiciable

R1=12m f=2.9m

e=0.22m a=0.35m

g=0.7m d=0.14m

e’=0.27m b=0.4m

h=15m

Densité du liquide : ɤ=950daN/m3

Charges : ɤv=1600daN/m3
ɤp= 1600daN/m3 φ=32 °

q=1KN/m²
Étanchéité int et ext =1KN/m² ρ béton =25KN/m3
Matériaux : fc28 =25MPa σsol=3.2 bars
P a g e 3 | 28
P a g e 4 | 28
Coupole et couverture
I- Calcul du rayon de la coupole :

2 2
R1−f 12.07 ²−2.97 ²
R1=12.07 mcos φ= 2 2
= =0.88⇒ φ❑=27.64 °
R +f1
12.07 ²+2.97 ²

r0 12.07
sin φ= =¿> R 2= =26.23 m
r2 0.46

II- Sollicitations :

 Charge permanente :

- Poids propre de la coupole P= 1×1 × 0.14 ×25=3.5 KN

- Etanchéité int et ext 1+1 =2 kn/m ²

- Terre végétale sur la coupole g × γ V =0.7 ×1600=1120daN /m² =11.2KN/m²

P a g e 5 | 28
 Total ¿ 1670 daN /m ²

 Charge d’exploitation :

La coupole est soumise à une charge d’exploitation de 100daN/m².

⇒ La charge par m² à ELS est égal à :

ps =1670+100=1770 daN /m² .

III- Dimensionnement de la ceinture :

Pour dimensionner la ceinture on va utiliser la théorie des coques mince appliqué à une
calotte sphérique de rayon R2=26.23m.

On a :

−a ⋅ p
L’effort normal suivant les méridiens : N ϕ = , avec a=R2
1+ cos φ

L’effort normal suivant les parallèles : N φ =a⋅ p ( 1+cos

1
φ
−cos φ )

a . p . cos φ
L’effort par unité de pourtour : H '=
1+cos φ

L’effort de traction dans la ceinture : T =H ' ⋅ R1

T
Acier de traction : A s=
σa

P a g e 6 | 28
 A.N :

−26.23 ×1770
N φ= =−246.182 daN /m.
1+ cos 27.64

N θ=26.23× 1770 ( 1+ cos127.64 −cos 27.64)=−165.105 daN /m.

26.23 ×1770 cos 27.64
H '= =21808.8 daN /mpourtour .
1+ cos 27.64

T =21808.8 ×12.07=263232.216 daN =263.232 Tonne.

2.6323 −4
A s= =159.56 ⋅1 0 m ²=159 cm² .
164.97

Soit 30HA25+4HA20.

P a g e 7 | 28
 Figure 2 : Ferraillage de la ceinture de la coupole

Vérification :
C.N.F
a ×b × f t 28 0.35 ×0.4 ×2.1
: A s ≥ Amin Amin= 400
=
400
→ 159 cm² ≥ 7.35 cm ²
❑

Dans la zone du recouvrement :

Longueur des barres d’ancrage droit

ϕ fe
L= × τ s=0.6 × ψ 2 × f t 28=2.83 MPa 0.8 × 400 =28.2 cm
4 τs 4 2.83

C = 0 → Lr = Ls
Soit les armatures transversales de diamètre : Ø 8
At (m ×п × ø × τ s )
St
= fet
At × fet 1× 400
St = (m ×п × ø × τ s ) = (11× π × 0.008× 2.835)

= 5.18 cm
→ Soit = 5.2 cm
Dans la zone courante :
St ≤ a → St ≤0.4m
Soit →St = 35 cm

IV- Dimensionnement de la couverture :

1- Calcul de ferraillage :

P a g e 8 | 28
 |Nρ| |Nρ| 246.182
σ b φ= = = =1.76 MPa ≤ σ bc=0.6 × fc 28=15 mp
be 1e 0.14
|Nθ| |Nθ| 165.105
σ bθ= = = =1.18 Mpa ≤ σ bc
be 1e 0.14

 Pas d’armature Asc

σ E 32164.19 ×0.14
cr= = =100.521MPa
a √ 3 (1−ʋ ) 26.23 √❑ 3(1− 1 ² )
()
2

6
σcr > σ b φ
σcr > σ bθ

Acier :

Après en prend 0.3% de la section pour calculer la section minimale de l’acier

Suivants les méridiens
A sc mini mal =1 × 0.14 ×0.3 %=4.2 cm² Soit 10HA8 A=5.024cm²/m

Suivants les parallèles

P a g e 9 | 28
A sc mini mal =1 × 0.14 ×0.3 %=4.2 cm² Soit 10HA8 A=5.024cm²/m

P a g e 10 | 28
 ETUDE DE LA PAROI
On suppose que les parois ne sont pas parfaitement encastrées sur le radier,
c’est le cas auquel se rapporte la résolution par la méthode de Hangan-Soare et
la méthode Méthode Lebelle qui suppose aussi que le radier est parfaitement
rigide.
I- Dimensionnement de la paroi :

Méthode de Lebelle:
a- Réservoir plein :
H = H L = 15m
e 0.22
R = R1 + = 12 +
L
= 12.11 m
2 2

β=
√4 3 = √4 3 = 0.806
√R × e √12.11× 0.22
P0 = H × ρ = 15 × 0.95 = 14.25 kdaN

P0 1 14.25 1
M (x=0) =-
0
× (1− β × H ) = - 2× 0.806² × (1− 0.806 × 15 ) = -10.06
2× β ²
kdaN.m/ml
P0 1 14.25 1
T =
0
× (2− β × H ) = 2× 0.806 × (2− 0.806 × 15 ) = 16.94 kdaN
2× β

P a g e 11 | 28
 Méthode de Hangan-soare
'
e 0.27
Hs = H + = 15 + 2 = 15.135 m
2
R HS = R1 + e = 12 + 0.11 = 12.11 m

β=
√4 3 = √4 3 = 0.8
√R × e √12.11× 0.22
β × H HS = 0.8 × 15.135 = 12.2

e 22
' =
=0.814
e 27

M0 = K × δ × H3

 Interpolation d’après abaque 1

K = 0.0031
M 0 = 0.0031 × 1 × 0.95 × 15.1353 = 10.21 kdaN.m/ml

X 0 = K0 × H

 Après interpolation: K 0 = 0.061 (Abaque 2)

X 0 = 0.061 × 15.135 = 0.92 m

 Vérification :
X 0 max = 1.2 √ R ×e = 1.2 √ 12.11×0.22 = 1.95 m ¿ X 0 = 0.92 m → OK

 Après interpolation : K 1 = 0.126 (Abaque 3)

X 1 = K 1× H = 0.126× 15.135 = 1.90 m

 Après interpolation : k ' = 0.00070 (Abaque 4)

M’= -k’ × δ × H 3 = -0.00070 × 1 × 1.05 × 15.1353 = -2.30 KdaN.m/ml
 Après interpolation : k 2 = 0.24 (Abaque5)
X 2 = 0.24 × 15.135 =3.63 m

 Après interpolation : k ' ' = 0.83 (Abaque6)

P a g e 12 | 28
Nφ
max
= k’’× δ × R × H = 0.83 × 1 × 0.95 × 12.11 × 15.135 = 144.52 KdaN/m

b- Cas d’un reservoir vide:

Méthode de Lebelle:
K0× γ R × H 1
M0 = × [ - (1- ) × θ ( βx ¿+ζ ( βx ¿
2× β ² β×H

γ R = 1600 daN/m3

ρ R = 32°

1−sin ρ R 1−sin 32
K0 = = = 0.267 = 0.25
1+ cos ρ R 1+ cos 32

T=
2β [
P0 −βx
e (cos β x−sin β x )+(1−
1 −βx
βh
)e sin β x T 0= ]
0.25 ×1.6 ×15
2× 0.806
×
(2-0.916) = 7.13

M=
P0
2β
2 [ −(1−
1 −βx
βh
− βx
)e cos β x +e sin β x ]
0.25 ×1.6 ×15
M0 =- 2 × 0.806² × (1-0.916) = 4.23 t.m/ml

Méthode de Hangan-soare:

M = K × δ × H = K ×δ × H × γ R =0.0031 × 15.135 ×1.6× 0.25= 4.3 kdaN.m/ml

0 3 3 3

M’ = -k’ × ( K 0 × γ R ) × H 3 = -0.00071 × 0.25 × 1.6 × 15.1353 = -0.98 KdaN.m

 Poids propre de la paroi :

P a g e 13 | 28
P.P = ρb é ton × 1 ml × H ×e = 25 × 1 × 14 × 0.22 = 77 KN/m

Poids ceinture = b × a × ρbéton × 1 ml = 0.46 × 0.41 × 25 × 1ml = 4.71 KN/m

Enduit voile = ρbéton × 1 ml × e enduit × H = 25 × 1 × 0.02 × 14 = 5.6

V’ = Nφ × sin 27.08 = |-21.154| × sin 27.08 = 9.63 kdaN/m²

Etanchéité = 2 × 1 ml × H = 2 × 1 × 14 = 28 KN/m²
Ns = P.P + poids ceinture + enduit voile + V’ + étanchéité
= 77 + 4.71 + 5.6 + 96.3 + 28 = 211.6 KN = 21.16 kdaN

Réservoir plein Réservoir vide

Lebelle HS Lebelle HS
Mo[MKdaN/m] -10,06 10,21 4,23 4,3
To[KdaN/m] 16.94 7,13
M'[MKdaN/m] 2,28 2,30 0,96 -0,97
Nφmax[KdaN/m] 140.47 144,52
Tableau récapitulatif

Dimensionement du paroi
ferraillage Cerce (traction simple)
T= 144.52 KdaN

N φ max 144.52
A=
❑
=
161.3
= 89.6 cm2

==> Choix d'armature 29HA20 =91.06cm2

P a g e 14 | 28
 Ferraillages en Appuis (flexion composé):

Calcul à l’ELS (fissuration très préjudiciable) :

Ns= poids de la ceinture + poids de la paroi+ V'

Poids d'étanchéité + poids de l’enduit
 Poids propre de la paroi :

=0.22*15*25*1=8.25 KdaN/m

 Poids propre enduit :

=0.02*20*15=0.6KdaN/m
 Poids propre d'etancheité :

=1*2*15=3kdaN
 Poids propre de la ceinture :

= 0.45*0.35*25=0.39 KdaN

 V' :poids propre de coupole :

=N𝛟*sin𝛟 = 246.182*sin27.64=11.42KdaN
===> Ns=8.25+0.6+3+0.39+11.42=23.66 KdaN/m
On a : Ns=23.66 KdaN/m

Mmax=10.21 KdaN/m (en appui)

Me 10.21
eser = Ns = 23.66 = 0.43m

h
Va =d− avec h=0.22 et d=0.9*0.22=0.198m
2

h
Va=d−
2
P a g e 15 | 28
 0.22
Va=0.198 – =0.088 m ==> e ser=0.43> v a=0.088
2

S.P.T= centre de pression est en dehors des armatures

Msa= Ns*ea

h
ea= es+(d- )
2

==>ea=0.43+0.088=0.518m

Msa=23.66*0.518=12.26kdaN.m =122.6Kn.m

bc=0.6 × fc 28=0.6 ×25=15 Mpa

Mrb=0.5∗b∗d 2∗bc∗α 1∗ 1− (
α1
3 )
Mrb=0.5∗1∗0.1982∗15∗0.582∗ 1− ( 0.582
3 )=13.79 KdaN . m

===> M sa< M rb ( A sc =0)

30 × MSA 30 ×0.1226
U= = =0.581
b ×d ² × σ st 1× 0.198² ×161.3

α 1 -3α 1²-3Uα 1+ ¿3U=0

3

α 1 -3α 1²-3×0.581α 1+ ¿3×0.581=0

3

α =0.560

b ×d × α 1 ² −2
1×(19.8 ×10 )×0.560²
Ast= = =47.04×10−4
30×(1−α 1) 30×(1−0.560)

P a g e 16 | 28
 NS 236.6 ×10
−3
A2 = Ast- = 47.04×10−4 - = 32.37 cm²
σ st 161.3

Ferraillages en travée (flexion composé):

Calcul à l’ELS (fissuration très préjudiciable) :

On a : N=23.66KdaN/m

Mser=2.3 mKdaN/m

M ser 2.3
eser= = =0.092 m
Nser 23.66

h
Va =d− = (0.9*0.22)-(0.22/2)=0.088 < eser
2
Donc centre de pression est en dehors des armatures
Msa=Ns ×e a
ea¿ 0.088+ 0.092=0.18m
 Msa=23.66∗0.18=0.0426 MPAS

 Mrb=0.5 × b ×d 2 ×bc ×α 1 × 1−( α1

3 ) α 1=0.582 et bc =15 Mpa

donc Mrb=0.5 × 1× 0.198 ² ×15 ×0.582 × 1− ( 3 )

0.582
=0.1379 MPa

====>Msa< Mrb (Asc=0)

30 × Msa 30 × 0.0472
U= 2
= 2
=0.202
b × d × 1× 0.198 × 161.3
α 13−3 α 1²−3 U α 1+3U = 0 α = 0.37
2 2
b ×d × α 1 1× 0.198 ×0.37
Ast= = =¿ 14.34cm²
30 × (1−α 1 ) 30× ( 1−0.37 )
P a g e 17 | 28
==> Choix d'armature 5HA20 = 15.70cm2

Vérification sur travée (ELU):

f = 0.7*L0 = 0.7*15 = 10.5 m

ea= max( 2cm ; L/250) = max( 2cm ; 1500/250) =6cm
Pu =1.35*(Pc + Et + Pt) +1.5*Q ==>Pu = 1.35*(3.5+ 2 +11.2) +1.5*1
Pu = 24.045kN / m2 = 2404.5 daN/m²
R2 × P −26.23 ×2404.5
Nφ = Nφ = - = 1+ cos 27.64 = -33443.29 daN = -33.44KdaN
1+ cos φ

Nu= 1.35* (8.25+0.6+3+0.39+11.42)=31.941kdan =319.41 kn

Mu= Mmax*1.35 = 2.3*1.35 = 3.105 kdaN.m =31.05kn.m

e0=Mu/Nu = 34.5/319.41=0.097 m
e1=e0+ea =0.097+0.06=0.157 m
Lf/h < max(15 ; 20*e1/h) ==> avec Lf/h =11.2/0.20= 56

et max(15 ; 20*0.157/0.22)=14.27
===> e2 n'existe pas

emoy= e1=0.157 m
Va = d-(h/2)=19.8-(22/2)=8.8 cm
Mua=( emoy + Va )* Nu = (0.157+0.088)*319.41=78.25 KN.m

fbu=(0.85*fc28)/1.5=14.17MPa

P a g e 18 | 28
 Mu A 78.25 ×10
−3
U UA = = = 0.14
b ×d ² × f bu 1× 0.198² ×14.17

U BC =0.8*(h/d)*(1-0.4*(h/d))

U BC=0.8*(0.22/0.198)*(1-0.4*(0.22/0.198))=0.494

U UA < U BC ===>Section partiellement comprimé

Pas d'acier comprimé

α=1.25∗¿=1.06
z=d*(1-0.4* α)= 0.198*(1-0.4*0.106)=0.114
Mser −3
Ast= Z∗fsu = 78.25∗10 = 19.73 cm2
¿¿
−3
319.41∗10
A2=19.73*10−4-( )=10.54 cm2
400/1.15

Ferraillage du Gousset:
Calcul à l’ELS (fissuration très préjudiciable) :
Pour éliminer l'acier comprimé et pour avoir MSA< Mrb on augmente l'épaisseur au niveau
bas du paroi ce que se traduit par un gousset

e= 60 cm (initialement e=22cm)
25∗0.6∗0.92
Ns= +23.66 = 242.29 Kn
2

10.21
es= = 0.421 m avec d=0.9*h = 0.9*0.6 = 0.54 m
24.22
Va=0.54-(0.6/2) = 0.24 m
ea= es+ Va= 0.421+0.24=0.66m
MSA= ea *Ns = 0.660*242,29 = 159.91 KN.m
Mrb=0.5*15*0.542*15*0.582*(1-(0,582/3))= 55.17KdaN.m

P a g e 19 | 28
 Mrb > MSA ====> Pas d'acier comprimé
30 × Msa 30∗0.2429
U= 2 = 2 = 0.13
b×d × 1∗0.54 ∗161.3
Avec α 13−3 α 1²−3 Uα 1+3U =0 α =0.31

2
b×d×α 1 1× 0.54 × 0.31²
A st= = =21.22 cm²
(
30 × 1−α 1 ) 30 × ( 1−0.31 )

Choix d'armature 7HA20 =21.98cm2

Dimensionnement Du Radier
Ceinture inferieure:
Fr= T0*R1= 16.94*10*12.11= 2051.14 KdaN
Fr
Ast= σst = 2051.14.10−3 /161.3 = 127.18 cm²

Choix d'armature 26HA25

Calcul de la surface du radier
F
σ sol =
S

F = poids liquide + Poids radier+ Poids paroi + Poids

ceinture+ Poids béton de propreté

2
3.14∗22
Poids propre du liquide: *15*9.5= 64465.5 KN
4
P a g e 20 | 28
Poids propre du radier : 3.14∗¿ ¿ *0.27*25= 3166.62KN
2
3.14∗24.44
P. Béton de propreté : *0.07*2.2 = 72.25 KN
4

Périmètre =R1*2π = 76.1m

poids propre du paroi = Ns*76.1 =242.9*76.1=18484.7 KN

Ftotale = 86189.1 KN

F
86189.1∗10 ∗4 =0.183 MPA <0.32 Mpa ok
−3
σ sol =
S 2
= 3.14∗24.44

σA = 8¿Q
Mmax = 0.083*Q*l
Q= Ns

σA=8*Ns/R2 = 8*242.9/24.44=7.95 Kdan/m²

Mmax=0.083*R2*Ns =0.083*242.95*24.44= 492.72 KN.m

Moment sur appui : 0.3*Mmax=0.3*492.72= 147.81KN.m=0.14781MN

Moment sur travée: 0.95*Mmax=0.95*492.72=468.084 KN.m=0.468MN

Permillage sur travée

Vérification ELS(fissuration très préjudiciable) :

d=0.9*0.27 = 0.243 m
Mrb=0.5*15*0.582*0.243²*(1-(0.582/3)) = 0.207MN.m
Mrb <Msa doncAsc ≠0 : Redimensionnement du radier

Soit h=0.6m ( 0.27 initialement )

Mrb=0.5*15*0.582*0.54²*(1-(0.582/3)) = 1.026 MN.m

Mrb> Msa

P a g e 21 | 28
 30 × Msa 30∗0.468
U= 2 = 2 = 0.2985
b×d × 1∗0.54 ∗161.3

avec α 13−3 α 1²−3 Uα 1+3U =0 α =¿ 0.44

2
b ×d × α 1 1× 0.54 × 0.44²
Ast= = =¿ 62.22cm²
30 × (1−α 1 ) 30× ( 1−0.44 )

On choisi 13 HA 25 (63.83 cm²)

Vérification ELU(fissuration très préjudiciable) :

−3
468∗10
= 2 = 0.201
(1∗0.54 ∗14.17)

α=1.25∗¿)= 0.283
z=d*(1-0.4* α)= 0.54*(1-0.4*0.283)=0.359

Mser −3
Ast= Z∗fsu = 468∗10 = 37.48 cm2 OK
¿¿

Ferraillage sur appuis

Calcul à l’ELS (fissuration très préjudiciable) :

30 × Msa 30∗0.14781
U= 2 = 2 = 0.167
b×d × 1∗0.54 ∗161.3
avec α 13−3 α 1²−3 Uα 1+3U =0 α =¿ 0.35

2
b ×d × α 1 1× 0.54 × 0.35²
Ast= = =¿ 25.44 cm²
30 × (1−α 1 ) 30 × ( 1−0.35 )

P a g e 22 | 28
Vérification ELU(fissuration très préjudiciable) :
−3
147.81∗10
= 2 = 0.063
(1∗0.54 ∗14.17)

α=1.25∗¿=0.082
z=d*(1-0.4* α)= 0.54*(1-0.4*0.082)=0.391
Mser −3
Ast= Z∗fsu = 147.81∗10 = 10.86 cm2 OK
¿¿

P a g e 23 | 28
Annexe

P a g e 24 | 28
P a g e 25 | 28
P a g e 26 | 28
 CONCLUSION

A la fin de ce mini-projet nous constatons qu’on a

appris beaucoup de choses au niveau pratique de
calcul des réservoirs.
Cependant, nous avons rencontré certaines difficultés
qu’on a pu dépasser grâce à l’aide et le soutien de notre
encadreur.
En fait, il faut reconnaître que ces difficultés nous ont
permis d’apprendre beaucoup de techniques concernant
le calcul, le dimensionnement et le ferraillage des
différents éléments d’un réservoir.

P a g e 27 | 28