REMERCIEMENT

En achevant ce travail, nous tenons à remercier et à exprimer notre gratitude à

Mr. Moez Trabelsi pour sa collaboration pour sa patience, son effort
d’encadrement et ses conseils tout au long de l’étude.
 I. Introduction

Les silos sont des réservoirs spéciaux destinés à stocker en vrac des
matériaux divers comme les grains, ciments charbon, minéraux, huile etc...
Dans ce bureau d’étude nous allons faire concrétiser nos connaissance vue en
cours à fin d’améliorer nos expérience dans les calculs et les méthodes utilisées
dans un dimensionnement d’un silo.
Ces réservoirs sont de formes diverses en carré, rectangulaire, circulaire
ou polygonaux. Ils peuvent être de très grande hauteur pour des dimensions
limitées en plan et d’une grande surface en plan pour les hauteurs réduites .Ils
peuvent être fermés ou ouverts au sommet .Leurs fonds peuvent être plans ou en
formes spéciales (mamelle ; trémies). Ils peuvent se trouver au niveau du sol,
enterrés ou surélevés. La couverture peut être munie d’accessoires divers
comme aération, dispositifs mécaniques …
II. Données du sujet :

 Hypothèses de calcul :

 Fissuration très préjudiciable / préjudiciable.

Densité du liquide : =1000daN/m3.
Paroi encastrée dans le radier et articulé la coupole

 Données géométriques :

R1=12m / f=3m / e=23m / g=0.65m / e’=0.28m

H=17m/ a=0.4m / d=0.14m/ b=0.4m

  Les charges :

Densité du liquide : =950m

-terre végétale sur la coupole : =1650 daN/m3

-terre de protection de la paroi : =1700daN/m3 et =29

-charge d’exploitation (d’entretien) : q=1KN/m2

-enduit intérieur (paroi et radier) et isolation extérieur : enduit=2cm

=20KN/m3

-béton =25KN/m3

-enrobage= 3cm

Matériaux :

-acier :

FeE 240 pour les aciers transversaux (cadres) .

FeE 400 pour les aciers longitudinaux.

-contrainte admissible du sol : 3.1bar

 I. Calcul de la coupole :

 Poids propre :

25×0,14=3,5daN/m2

 Terre végétale :

1650×0,65=1072,5daN/m2

 Etanchéité :
2×1KN/m2=200daN/m2

PTOTALE=1376daN/m2
2 2 2 2
r +f 12 , 07 +3 , 07
 R 2= = =25 , 26 m
2f 2 ×3 , 07
−1 R1 12 , 07
 sin R 2 = 25 , 26 =28 , 54=¿ φ=28 , 07
−R 2 × P −25 , 26 ×1376
 N φ= 1+cos φ = 1+cos 28 , 54 =−19235 , 59 daN

[ ]
❑
1
 Nθ=a× p
1+cos φ
−cos φ

[ ]
❑
1
¿ 25 , 26 ×1376 −cos 28 , 54
1+cos 28 , 54

¿−12031 ,04 daN /m2

les quatres vérifications :

 −|N φ|
 σ bφ= 100 × e =1,374 MPa< 15 MPa
'

−|Nθ|
 σ bφ= 100 × e =13 , 74 MPa<15 MPa
'

 σ ' bφ=σcr
 σ ' bθ=σcr Amin=0,25 a 0,35%

 Eij=11000 3 × √3 fc 28

¿ 11000 × √ 25=32164 , 2 MPa

3 3

Ee 32164 ,2 ×0 , 14
Fr= = =101 , 52 MPa
 a × √ 3(1+ϑ ❑ ) 25 , 26 × √ 3(1+( 1 )❑2)
2

3 ET 4 vérifié

Amin=100×0,14×0,25=3,5cm2/m

Choix = 10T8 As= 6.28

d=0,14m>0,09m  double nappe

calcul de longueur de ruissellement :

∅ fe
ls= ×
4 τs
2
τs=0 , 6 × ωs × ft 28

ωs=1 , 5

fc28=25 ; ft28=2,1Mpa .
2
τs=0 , 6 ×1 , 5 × 2 ,1=2,835

0,8 400
ls= × =28 ,22 cm<a
4 2,835

a=40cm vérifier

II. Calcul de la ceinture supèrieur :

Dimensionnement de la ceinture :

Calcul de la ceinture
4 4
p(r −f )
H '=
4 f r2

1376 ( 12074 −3 , 074 )

¿ 2
4 ×3 , 07 ×12 , 07

¿ 16256 , 4 daN /m
'
T =R × H

¿ 12 , 07 ×16256 , 4

¿ 1962 ,1 KN

Section d’acier
T
A=
σs

σs=?

Fissuration très préjudiciable

{
fe
σs >min 2
Max { 0 , 4 fe ; 90 √ n × ftj

min
{ 200
Max { 160; 165 }

σs=165 Mpa

1962 , 1
A=
165
2
¿ 0,012 m

A=113cm2 ==> 25HA25  129,75cm2

 Les armatures transversales

Lr=Ls +C
avec C=3+0.8+1.25=5 , 05 cm

Ls=28 , 2cm

Lr=28 , 2+5 , 05=33 ,3 cm

At=4 ∅ 8=2, 01 cm 2

St =40 cm

τs=2,835 Mpa

III. Dimensionnement de la paroi vertical

Cas du Réservoir plein:

 Méthode hangar Soares
1,3165 1,3165
B= = =0,78
√ R e √12 , 23∗√0 , 23
H=17,14m

R=12,23m

K=0,0025

K’=0,000502

K’’=0,834M0

M0=k×H3× δ avec δ =0,95 Kdan/m3

M0=0,0025×17,143×0,95=11,95mkdan/m

M’= -k’×H3× δ =2,4mkdan/m

Nφ max = k’’× δ × R HS×H

=0,834×0,95×12,23×17,14

=166,08Kdan

 Méthode lebelle:

H=17m P0= ρ gh=0,95×17=16,15Kdan/m2

R=12m
1,3165 1,3165
B= = =0 , 82
√ R e √12∗√ 0 ,23
−P0 1
M 0= 2
(1− )
2B Bh

−16 ,15 1
M 0= (1− )
2 ×0 , 82
2
0 ,82 ×17

= -11,25kdan/m
−P0 1
T 0= (2− )
2B Bh

−16 ,15 1
T 0= (2− )
2× 0 , 82 0 ,82 ×17

= 18,99 kdan/m

Mmax=2,4mdan/m

Réservoir vide
 Méthode lebelle

M (0 )=
−P 0
2× B
2 ( (
−1
−1
B× H))
avec

1−sin φ 1−sin 29
K= =
1+sin φ 1+sin 29

¿ 0 , 35

P0=k ×γ × H

¿ 0 , 35 ×1 ,7 × 17=10,115 kdan /m

M (0 )=
−10,115
2× 0 , 82
2
−1
( ( −1
0 , 82 ×17 ))
=7 mkdan/m

−P0 1
T 0= (2− )
2B Bh

T 0=
−10,115
2× 0 , 82
2− ( 1
0 , 82 ×17 )
=11, 9 kdan /m

 Hangar soare
M0=k×H3× δ avec δ =0, 35×1,7=0,595 Kdan/m3

M’= -k’×H3× δ

=0,00502×0,595×17,143=1,5mdan/m

Tableaux recapitulative

reservoir plein reservoir plein

Lebelle hongarSoar Lebelle hongarSoar

M0 KdaN/m -14,64 11,95 -7 7,5

T0 KdaN/m 18,99 ------- 11,9 -------

MAX KdaN/m 2,4 2,4 1,56 1,5

N KdaN/m 163,7 166,08 ------- -------

IV. Calcul ferraillage de la paroi verticale :

 Calcul en travée :

Poids propre du vertical : 25×0,23×1=5,75kn/m2

Poids propre de l’enduit : 20×0,02×1=0,4kn/m2

Poidspropre de l’étanchéité : 2×1=2kn/m2

∑ ¿ 8 ,15 kn /m2
N =8,15×17=183,55 kn/m2

Poids de la ceinture

Poids propre de la ceinture : 25×0,4×0,4 = 4kn/ml

Poids propre de l’étanchéité : 1×(-40-23)+(40-14)=0,9kn/ml

∑ ¿ 4 , 9 kn /ml
Poids du coupole : δ=26 , 5 ° donc φ=166 , 08 ×0 , 14=23 , 25 kdan
'
N p =23 , 25× sin 29=11 ,27 kdan/ml
 '
N p =112 , 7 kn/ml

∑ ¿138 , 55+ 4 , 9+112 , 7=255 , 45 kn /ml

ELS : G + Q =255,45 + 17 = 272,45 kn/ml
M S 2 , 4 10
eS= = =0,088 mkdan/m
N S 272 , 45

' e 0 , 23
e S= = =0,038
6 6

Es ∈ pas { -0,038 ; 0,038 }  donc section partiellement comprimé.

Très préjudiciable ; comprimé.

1 2 α1
M rb = b d σ bc α 1 (1− )
2 3

9 Fc 28
α 1= =0,5824
σ s + 9 fc 28

MIN ¿

σ st =0 ,8 × 201,633=161,306 Mpa

En propose que ∅ 14 
1539
d=0 ,23− −0 , 03=0,199 m
2

e A =e S +V A

h 0 ,23
V A =d− =0,199− =0,084 m=8 , 4 cm
2 2

e A =0,088+0,084=0,172 m

M SA =N S × e A =272 , 45 ×0,172=46 , 86 kNm

σ bc=0 ,6 × 25=15 Mpa

1
M rb = ×1 ×0,199 ×15 × 0,5824 1−
2 (
0,5824
3 )
M rb=0 , 14 Mpa

M rb=140 Kpa
 M SA 46 , 86
U ser= = =78 ,88 KNm
b d σ bc 1 ×0,199 2 ×15

140
U rb= 2
=235 , 7 KNm
1 ×0,199 ×15

U ser <U rbPas d’armature comprimée.

30 M Ser −3
30 × 46 , 86 10
U= 2
= 2
=0 , 22
b d σ st 1 0,199 161,306
3 2
∝1−3 ∝1−3 u ∝1 +3 u=0
3 2
∝1−3 ∝1−(3 × 0 , 22∝1)+30,1985=0

∝10,3922

2
b d ∝1 1 ×0,199 × 0,3922
2
2 2
A st = = =0,001678 m =16 ,78 cm
30(1−∝1 ) 30(1−0,3922)

choix de l’armature 10HA16 (section en travée)

 Calcul de l’appui :
M S 11, 95 10
eS= = =0,4386
N S 272 , 45

d=0 , 9× 0 , 23=0,207 m

d '=0 ,1 × h=0,017 m

h
V A =d− =0,207−0,085=0,122 m
2

e A =e S +V A =0,4386+0,122=0,5606 m

M SA =0,5066 ×272 , , 45=138,023KNm

1
M rb= ×1 ×0,207 2 × 15× 0,5824 1−
2
0,5824
3 (
=0,1508 Mpa )
¿ 150 , 8 kMP

M SA 138,023
U sa = 2
= 2
=214,743 KNm
b d σ bc 1× 0,207 × 15

151
U rb= 2
=214,743 KNm
1 ×0,207 ×15
U Sa<U rb

30 M Ser −3
30 × 46 , 86 10
U= 2
= 2
=599 , 07 Kpa=0,599Mpa
b d σ st 1 ×0,207 161,306
3 2
∝1−3 ∝1−(3 × 0,599 ∝1 )+3 ×0,599=0

∝10,566

b d ∝21 1 ×0,207 × 0,5062

A st = = =0,003576 m2=35 , 76 cm2
30(1−∝1 ) 30(1−0,506)

choix de l’armature 18HA16=36,18(section en travée)

 En fait l’augmentation + 30cm

E + 30 cm

D=0,477
V A =0,212

e A =0,651

M SA =177 , 09

M rb=800 KN

U SA=51 , 88

U rb=234 , 4

U Sa<U rb

 U=144,75
3 2
∝1−3 ∝1−(3 × 0,144 ∝1)+3 ×0,144=0

∝10,33

2
b d ∝1 2
A st = =25 ,8 cm
30(1−∝1 )

choix de l’armature (10HA14+10HA10)=23,24 cm2 (section en travée)

 En fait l’augmentation + 40cm

E + 40 cm

D=0,567
V A =0,252

e A =0 , 69

M SA =188

M rb=234 ,65

U SA=38 , 98

U rb=234 , 65

U Sa<U rb

 U=108,76
3 2
∝1−3 ∝1−(3 × 0,1087 ∝1)+3 ×0,1087=0
3 2
∝1−3 ∝1−0,326 ∝1+ 0,326=0

∝10,29

b d ∝21
A st = =22, 4 cm2
30(1−∝1 )

choix de l’armature (12HA16)=24,12 cm2

 Vérification a ELU d’armature longitudinale en travée

Lf =0 ,7 × 17=11, 9 m

(
e A =max 2 cm ;
L
250 ) (
=max 2 cm;
1700
250 )
=6 , 8 cm

−R 2× Pu −5 , 26 ×1620
Nφ u= = =−20,925166 ,08 KdaN /m
1+ cos φ 1+cos 26 ,5

N u=1 , 35 ( 25 ×0 , 42 ) +25 ( 1 ×0 , 23 ×17 )+ 1, 35 ×2,825+ ¿

Mu
e A= =0,172
Nu

e r=0,208
Lf 20 e
<max ⁡(15 ; )
h h

Lf 11, 9
= =51 , 74
h 0 , 23

(
max 15 ;
20 × 0,208
0 ,23 )
=18 , 06

e maj =0,302

Mua=( e maj + Va ) Nu

h 0 ,22
Va=d− =0 , 9 ×0,022− =0,088
2 2

Mua=66 , 39

M ua
μua= 2
=0,209
b d ×fbu

μ BC =0 , 8
h
d ( h
1−0 , 4 =0,562
d )
μua < μBC

∝=1 , 25 ¿)

∝=1 , 25 ¿)=0,235

Z=d ¿)=0,173

Ast =10,72cm2
Nu 204,026 2
A2= A sr − =9 , 28− =4 , 52 cm
fe 347,826
σe

ft 28 e 0−0,455 d
Amin=0 , 23 bd ×
fe e 0−0,185 d

ft 28=0 , 6+0 , 06 ×25=2 , 1 MP

d=0 , 9× 0 , 23=0,207

b=0 , 23

e 0=0,152
−5 2
Amin=2 , 9210 10 m
A2  Amin 4HA12

A2 =4,52m2 verifier

 Vérification de calcule :

En travée
c=d−e A=0,199−0,172=0,027

90 A 1 90 A 2
p=−3 c −( c−d ) ×
2 '
+(d −c)×
b b

A1=0 ; non pas d’acier comprimé

A2=0,001678 mm2

2 90 A 2
p=−3 c +(d−c)×
b

2 90× 0,001678
p=−3× 0,027 − ( 0,199−0,0027 ) × =0,0238
1

2 90 A2
q=−2 c −( 0,199−0,0027 ) ×
b

2 9 × 00,001678
q=−2 ×0,027 −( d−c ) × =−0,0045
1
3
2 4p
∆=q +
27
3
2 4 × 0,238
¿−0,0045 + =0,002>0
27
3
yc + pyc+ q=0
3
yc +0,0238−0,0045=0

yc=0,056 ; y ser= y c + c

y ser =0,056+ 0,027=0,083

3
1 × y ser 2
I SRH = +15 A 2 (d− y ser)
3
3
1 × 0,023 2
I SRH = +15 × 0,001678(0,199−0,083)
3
I SRH =0,00053 m4

−3
10,27245 10 ×0,056
K= =0,0029
0.0053

σ bc=k × y ser =0,0029 ×0,083=0,00024 <15 Mpa

σ st =15 × k (d − y ser )

σ st =15 × 0,0029 ( 0,199−0,083 )=0,005 Mpa<<161,3Mpa

Donc vérifier

Sur appuis e+40cm :

D=0,567
c=d−e A=0,567−0 ,69=−0,123

2 90 A 2
p=−3 c + ( d−c ) ×
b

2 90 ×0,003576
p=−3×−0,123 −( 0,567−0,123 ) × =0 , 17
1

3 2 90 A 2
q=−2 c −( d −c ) ×
b

3 2 9 ×0,003576
q=−2 ×0,027 −( 0,567−0,123 ) × =−0,157
1
3
2 4p
∆=q + =0,025>0
27
3
yc +0 , 17 c−0,157=0

yc=0 ,32 ; y ser = y c + c

y ser =0 , 32+0,123=0,197
3
1 × y ser 2
I SRH = +15 A 2 (d− y ser)
3
3
1 × 0,197 2
I SRH = +15 × 0,003576(0,567−0,197)
3

I SRH =0,0855m4
 272,452 ×0 , 32
K= =1,018
0.0855

σ bc=k × y ser =1,018 ×0,197=0 , 2 Mpa <15 Mpa

σ st =15 × k (d − y ser )

σ st =15 ×1,018 ( 0,567−0,197 )=5 , 64 Mpa<<161,3Mpa

Donc vérifier

V. Calcul ferraillage des cerces :

Nφ max =166 , 08 KdaN

Nφ max 2
A s= =100 cm
σs
soit le choix32HA 20

VI. Dimensionnement du radier :

 Dimensionnement de la ceinture de radier:
b=e ' =0 , 28

e 0 , 23
R=R 1+ =12+ =12,115 m
2 2

Fr=T ° × R=18 , 99 ×12,115=230,064 KdaN

Ns 230,064 2 2
A s= = =0,0143 m =143 cm soit le choix30HA 25
σ st 161, 3 10 2

Notre section n’est pas suffisante acause du ce nombre de notre choix

d’armature donc en va elarjie une coté

 Calcul du radier

Poids du liquide
2 2
π R h δ=π 12 × 17 ×950=7306 , 09 KdaN /m
Béton de propriété

12,232× 0 ,1 × π ×25=1174, 74 KdaN /m

Poids radier :

12,232× 0 ,28 × π ×25=328 , 93 KdaN /m

Poids du paroi :

π R 2 h−( π ( R+e )2 h ) 25=[π 122 × 0 , 23−( π ( 12+0 , 23 )2 17 ) ]25

¿ 7440 , 8KN/m=744,08Kdan/m

Poids enduit paroi:

2 π R1 h xγee=2 π 12 X 17 x0 , 02 ×2=51 , 27 KdaN

Poids du coupole:
S ×d × γ ba

Il s’agit d’une demi-sphère creuse de rayon R2

R2=25,26m φ=26 , 5 °

φ
S=2 π R
2
2 ∫ sinφ dφ
0

S=2 π R2 [ 1−cos φ ]
2

2 2
S=2 π 25 , 26 0,108=432, 89 m

Poids du coupole :

432,98 x 0,14 x 2,5 =151,543 KdaN /m

Charge exploitation:
Qe=432 , 89× 0 , 1=43,298 KdaN /m

Poids du terre végétale du coupole :

432 , 89 ×0 , 65 ×1 , 65=464 , 37 KdaN /m
Poids du ceinture
a × b ( 2 π × R 1) × 2, 5=0 , 4 × 0 , 4 ( 2 π × 12 ) × 2, 5=30 ,16 KdaN /m

F=∑ Poids=9197,077 KdaN /m

F 9197,077 KdaN
σ= = =20 , 33 =0,2033 Mpa=2 , 03 ¯¿ 3 , 1 ¯¿
S π ×12
2
m

OK vérifier

 Détermination du schéma de ferraillage :

Q = p’ composante verticale du Nφ

Poids du paroi

Poids ceinture

Pour QG = P’G + Poids du paroi + Poids ceinture

Pour QQ =P’Q ; q=1KN /m2

Avec P = G + Q
2
G=0 ,7 × 1650+0 , 1+ 2, 5 ×0 , 14=1,605 kdaN /m

Q =0,1 KdaN /m2

P =1,708

' R2 × p ×sin φ 25 ,26 ×1,705 × sin 26 , 5

PG = = =13,291,605 kdaN /m
1+sin φ 1+sin 26 , 5

' R2 × p ×sin φ 25 ,26 × 0 ,1 ×sin 26 ,5

PQ = = =0 ,78 kdaN /m
1+sin φ 1+sin 26 , 5

Poids de la paroi :
e × h × γ BA =0 , 23 ×17 × 2, 5=9,775 kdaN /m

Ceinture :
a × b × γ BA =0 , 4 × 0 , 4 × 2 ,5=0 , 4 kdaN /m

D’où

QG=13,29 + 9,775+ 0,4 =23,465 kdaN /m

QQ=0,78kdaN /m

QTOT =23,456+0,78= 24,245 kdaN /m

Et d’où la répartition :
Q ×2 π × R 1 2 Q 2 ×24,245 2
q Q= 2
= = =4 , 04 kdaN /m
2π ×R 1
R1 12

2 G 2 ×13 , 29 2
q G= = =2,215 kdaN /m
R1 12
2
Q=( qG + qQ ) × R 1=( 4 , 04+2,215 ) ×12=75 , 06 kdaN /m

Q par un terrain excellent :

6Q
σ A= L=l=2 R 1=∅ R 1
¿
Q R1 75 , 06 ×12
Max= + M 0= +11, 95=112 , 03 KdaNm /m=1, 12 MNm/m
9 9

B=1m et h=0,28m
M ser =M max=1 , 12 MNm

M rb=0,622 MNm< M ser

A sc ≠ 0

 Pour cela on choisit de faire augmenter l'épaisseur de la Radier

on prend e'=50 cm.

1 2 α1
M rb= ∝1 b 0 d σ bc α 1 (1− )
2 3

M rb =0 , 5 ×0 , 58 ×1 ×0 , 45 1−
2
( 3 )
0 ,58
× 15

¿ 0,707 MNm < M ser =1, 12 MNm

A sc ≠ 0

' d ' 0 ,05

δ= = =0 ,11
d 0 , 45
 '
∝1−δ
σ SC =15 σ bc ; σ bc =0 , 6 fc28=15
∝1

0 ,58−0 , 11
σ SC =15 ×15 =182,327 Mpa
0 , 58

' M ser −M rb1 , 12−0,707

2
0,0056 m 2
A = ser = = =56 , 6 c m /m
( d −d ) σ SC ( 0 , 45−0 , 05 ) 182,327
'
m

Choix20HA20 (62, 8c m2)

Z b=d 1− ( ∝1
3 )
=0 , 45 1− (
0 , 58
3
=0,363 m )
M rb ' σ SC 0,707 −4 182,327
A st = +A = +56 , 6 10
Zb × σS σ S 0,363 ×161 , 3 161 ,3
2 2
A st =0,018 m =18 cm

Choix 10HA16 (20,1c m2 /nappe)

ANNEXE