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 :

1. 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

2. 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

 3. 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 exterieur : 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 quatre 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

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 verifier

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τ

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,08 Kdan

 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

III. Calcul ferraillage de la paroi verticale :

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

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

Poid proprede l’etanchité : 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

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

∑ ¿ 4 , 9 kn /ml
Poids du coupole : δ=26 , 89 ° 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

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

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
 b d ∝21 1 ×0,199 × 0,39222
A st = = =0,001678 m2=16 ,78 cm2
30(1−∝1 ) 30(1−0,3922)

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

Calcul de l’appuis
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 3 (
0,5824
)
=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

b d ∝21
A st = =25 ,8 cm2
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 (10HA14+10HA10)=23,24 cm2 (section en travée)

 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

Verification 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

IV. Calcul ferraillage des cerces :

 Nφ max =166 , 08 KdaN
Nφ max 2
A s= =100 cm
σs
V. Dimensient du 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 choix
σ st 161, 3 102

 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 , 89 °
φ
S=2 π R2∫ sinφ dφ
2

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