2018
DIMENSIONNEMENT D'UN
RESERVOIR CIRCULAIRE
Majdi Ben Amor
Marouane Abid
SOMMAIRE
1 - INTRODUCTION
2 - DIMENSIONNEMENT DU COUPOLE
3 - DIMENSIONNEMENT DU CEINTURE
4 - DIMENSIONNEMENT DU PAROI
5 - DIMENSIONNEMENT DU RADIER
6 - ANNEXE
1 - INTRODUCTION
Les réservoirs sont des ouvrages réalisés en béton armé ou béton précontraint, destiné au stockage des liquides.
On distingue alors plusieurs types de réservoirs, classés suivant plusieurs critères dont, la position par rapport au
sol, la forme de la cuve, le mode de fermeture, la complexité de la construction, l’usage et même la nature du
liquide conservé.
Pour notre mini-projet, nous sommes face un réservoir couvert, reposant sur un radier et avec un sol de remblais
aux alentours de ses parois. On peut donc dès lors le considéré comme étant enterré. Les données de notre projet
sont les suivantes :
Données : R1=14 f=3.1m
e= 24 cm a=40cm
g= 70 cm d=14cm
e’=30 cm b=40cm
h=17 m
Hypothèse de calcul :
Fissuration très préjudiciable.
Densité du liquide : ᵧ =1150 daN/m3
Paroi encastrée liée avec la coupole.
* Charge:
Terre végétal sur la coupole : ᵧ
s=1650 daN/m3
ᵧ
Terre de protection de la paroi : { p=1900 daN/m3 ; =32°} ᵩ
Charge d’exploitation (d’entretien) : {enduit=2cm ;ρenduit=20 kN/m3} ;Etanchéité intérieure et extérieure
(paroi et coupole) : 1 kN/m3
Béton : ρbeton=25 kN/m3
Matériaux Béton : fc28=25 MPa Enrobage : c =3 cm
Acier :
FeE240 pour les aciers transversaux (cadres).
FeE400 pour les aciers longitudinaux. Contrainte admissible du sol : σsol=3.3 bars
2 - DIMENSIONNEMENT DU COUPOLE
Les coupoles sont des coques destinées à couvrir un espace de ferme la plupart du temps circulaires, mais
aussi polygonal ou de toute autre forme.
𝑑 0.07
𝑓+ 3.1+
2 2
tan 𝜃 = 𝑑 = 0.07 =0.2253
𝑅1+ 14+
2 2
=> 𝜃 =12.696 ᵒ
or 𝝋= 2* 𝜃= 25.4 ᵒ
Evaluation des charges sur la coupole :
Poids propre de la coupole Pc beton * d Pc 25x0.14 3.5kN / m2
Etanchéité intérieure et extérieure de la coupole : Et= 2 kN/m2
Poids de la Terre végétale sur la coupole Pt v .g Pt 16.5*0.7 11.55kN / m2
Charge de l’entretien : Q= 1kN/m2
ELU
Pu 1.35.(Pc Et Pt) 1.5.Q Pu 1.35.(3.5 2 11.55) 1.5x1 24.517kN / m2
ELS
Ps (Pc Et Pt) Q => Ps (3.5 2 11.55) 1 18.05kN / m2
Les calculs seront menés à l’Etat limite services (ELS).
Déterminons les efforts N suivant le méridien et N suivant la médiane
− 𝑅2∗𝑃𝑠 𝑅1 14
N1+cos 𝜑 avec R2 = sin 𝜑
=
sin 25.4
= 32.8 m
− 32.8∗18.5
N 1+cos 24.5 N
1 1
N 𝜽= 𝑅2 ∗ 𝑃𝑠 ∗ (1+cos 𝝋 − cos 𝝋) =32.8*18.5*(1+cos 24.5 − cos 24.5)
= -22933.5 daN
Vérification a l’état limite du béton en compression :
Vérification au flambement de la coupole:
Amin 0.35%(b*e) = 0.35 *100 *0.14 = 4.9 cm2/m
=> CHOIX 10 HA8
Longueur de scellement dans la ceinture
𝐿𝑠
= 35.3 => Ls=35.3*0.01=0.353m => 35cm
∅
3 - DIMENSIONNEMENT DU CEINTURE
st= 161.304 N/m2
Ainsi A=2.93*14.2*103/161.304=0.0025793m2= 257.9 cm2
=> CHOIX 50 HA25
4 HA20
Longueur de recouvrement
Lr=0.6*Ls=0.6*(35.3*0.025) =
Longueur développé et longueur réelle
Longueur
Longueur Longueur
développé(m)
projeté(m) réelle
14.54 14,55 14.9
14.25 14.35 14.7
13.84 13.88 14.23
12.81 12.88 13.23
11.16 11.27 11.62
7.85 8.03 8.38
4 - DIMENSIONNEMENT DU PAROI
Réservoir plein:
a/ Méthode de Lebelle:
A l'aide d'un fiche Excel qu'on a déjà programmé pour avoir un tableau
décrivant la variation de moment tout au long de paroi en introduisant les
formules suivantes :
Apres beaucoup d'itérations et ajustement de l'intervalle des x on a pu déterminé
le moment fléchissant maximale et son abscisse
Mmax=3.982 kdaN.m en x=2.137 m
et M0= -17.5476kdaN.m et T0=Tmax= 26.218 kdaN
=15.04kdaN
b/ Méthode de Hangan Soare:
hHS=h+e'/2 =17+(0.3/2) = 17015 m
RHS=R1+e = 14+0.24 = 14.24 m
M k.hHS) k est déterminé sur abaque n°1 et par interpolation a
partir de β*hHS= 0.712*14.24 = 12.21
donc on obtient M= 0.00307*1.15*(17.15)3= 17.808 kdaN.m
X0=[Link] = 0.0603*17.15 = 1.034
ensuite abaque ( 2) nous donne X'0=X0+(π/β)= 1.034+( π/0.712)=5.446m
Vérifions que X0(max)=1.2*( RHS*e)^0.5 =2.2184m > X0 OK
X1=[Link]= 0.127*17.15= 2.178 m
X1(max)= 1.8*( RHS*e)^0.5= 3.327 m > X1 OK
M' -k'.hHS)=-0.0006685*1.15*(17.15)3 = -3.8778 kdaN.m
X2=[Link] = 0.213*17.15=3.652 m
Nmax k ''. .R.h =0.818*1.15*14.24*17.15=229.73 kdaN.m
Réservoir vide:
a/ Méthode de Lebelle:
en utilisant le même fichier Excel on a déterminé Mmax et sa position
Mmax=1.671 kdaN.m
X= 2.16 m
b/ Méthode de Hangan Soare:
on a ka remb =0.254*1.9=0.4826
Moment à l’encastrement inferieur M0 = k* remb *h
=0.00307*0.4826*(17.15)3=7.473 kdaN.m
M' -k'.hHS)
= -0.0006685*0.4826*(17.15)3= -1.627 kdaN.m
TABLEAU RECAPITULATIF
RESERVOIR PLEIN RESERVOIR VIDE
Lebelle HS Lebelle HS
M(x=0) [Link]/m -17,547 17.808 -7.364 7.473
Mmax([Link]/m) 3.982 -3.8778 0,914 0,9866
max (daN/m) 212.381 229.73 ------ ------
kda 26.218 ----- 11 -----
4 - DIMENSIONNEMENT DU PAROI
Dimensionnement sur travée:
Pc = 25*0.7*0.55 = 9.625 KN/ml
V' = 31.89*10*sin(25.4) = 136.678 KN/ml
PP paroi = (25*0.24*1*17)+(20*17*0.02) = 108.8 KN/ml
Etanchéité =1*1*17*2=34KN/ml
Ns=289.205 KN/ml = 0.289 MN.m/ml
a/ Sur travée(ELS):
d=0.9*h = 0.9*24 = 21.6 cm = 0.216 m
eA=eS+Va avec Va = d-(h/2)=21.6-(24/2)=9.6 cm
et eS=Ms/Ns = 39.82/289.205 = 0.137 m
eA=0.137+0.096 = 0.233m
MSA=NS* eA= 289.205*0.233 = 67.384 KN.m
=>
= 161.304 MPa
avec => = 0.582
Mrb= 0.5*1*(0.216)2*15*0.582*(1-(0.582/3))
= 144.39 KN.m > MSA
=(30*67.384*10-3)/(1*0.2162*161.304)
= 0.2686
=> α1=0.424
= (1*0.216*0.4242 )/(30*(1-0.424))
2
= 22.47 cm
b/Vérification sur travée(ELU):
Lf = 0.7*L0 = 0.7*17 = 11.9 m
ea= max( 2cm ; L/250)
= max(2cm ; 1700/250) =6.8 cm
Pu 1.35.(Pc Et Pt) 1.5.Q Pu 1.35.(3.5 2 11.55) 1.5x1
24.5175kN / m2 = 2451.75 daN/m2
− 𝑅2∗𝑃𝑢 𝑅1 14
N avec R2 = = = 32.8 m
1+cos 𝜑 sin 𝜑 sin 25.4
− 32.8∗2451 .75
N daN
1+cos 24.5
Nu= 1.35*[(25*0.7*0.55)+(25*1*0.24*17+20*17*0.02)+(1*1*17*2)]
+[49.97958*sin(25.4)*10]
= 411.891 KN
Mu= Mmax*1.5
= 3.982*1.5 = 5.973 kdaN.m
e0=Mu/Nu
= 59.73/411.891=0.145 m
e1=e0+ea
=0.145+0.068=0.213 m
Lf/h < max(15 ; 20e1/h)
avec Lf/h =11.9/0.24= 49.58
et max(15 ; 20*0.213/0.24)=17.75 => e2 n'existe pas
emoy= e1=0.213 m
Va = d-(h/2)=21.6-(24/2)=9.6 cm
Mua=( emoy + Va )* Nu
= (0.213+0.096)*411.891=127.27 KN.m
fbu=(0.85*fc28)/1.5=14.17MPa
=127*10-3/(1*0.2162*14.17)=0.1925
=0.8*(h/d)*(1-0.4*(h/d)
=0.8*(0.24/0.216)*(1-0.4*(0.24/0.216)
= 0.494
< =>Section partiellement comprimé
< =0.39 Pas d'acier comprimé
α=1.25*(1.(1-2 )0.5)=0.2697
z=d*(1-0.4* α)= 0.216*(1-0.4*0.2697)=0.192
Ast= Mser/(z*fsu)=127.27*10-3/(0.192*(400/1.15))
=19.08 cm²
A2= Ast-( Nu/fsu)=19.05.10-4-(411.891.10-3/(400/1.15))
= 7.208 cm²
c/ Condition de non fragilité C.N.F:
e0-0.455*d
* =11.16 cm²
e0-0.185*d
=>A2=max(22.47 ; 7.208 ; 1.16) = 22.47 cm2
Cerces:
As=max{Nu/(fe/γs) ; Nser/σst }
avec Nu1.5* Ns=344.595 daN
=max{344.595.10-2/(400/1.15) ; 229.7.10-2/161.304}
= max{0.0099071 ; 0.014242}
Ainsi As=142.42 cm² CHOIX 46 HA20
Dimensionnent sur appui:
d/Sur appui ELS:
M0 = 17.808 kdaN.m
eS=M0/Ns = 178.08/289.205 = 0.615 m
d=0.9*h = 0.9*24 = 21.6 cm = 0.216 m
eA=eS+Va avec Va = d-(h/2)=21.6-(24/2)=9.6 cm
eA=0.615+0.096 = 0.711 m
MSA=NS* eA= 289.205*0.711 = 205.624 KN.m
=>
= 161.304 MPa
avec => = 0.582
Mrb= 0.5*1*(0.216)2*15*0.582*(1-(0.582/3))
= 144.39 KN.m < MSA ====> Asc 0
avec d'=0.1*h= 0.024 m
=180.254 MPa
= = 17.83 cm²
=76.42 cm²
A2=Ast-(Ns/σst)
= 76.42.10-4-(289.205.10-3/161.304)
= 58.49 cm²
e/ Gousset:
pour éliminer l'acier comprimé càd avoir MSA< Mrb on a programmé
une fiche Excel et après un calcul itératif on a déterminé un nouveau
épaisseur au niveau bas du paroi ce que se traduit par un gousset
on alors e= 0.4 m ( 0.24m initialement )
Ns=289.205+((0.32*1*0.34)*25)/2
= 293.341 KN
es= 178.08/293.341 = 0.607 m
d=0.9*h = 0.9*0.4 = 0.36 m
Va=0.36-(0.4/2) = 0.16 m
ea= es+ Va= 0.767 m
MSA= ea *Ns = 0767*293.341 = 225 KN.m
Mrb=0.5*15*0.582*0.36²*(1-(0.582/3)) = 0.456 MN.m
Mrb > MSA ====> Pas d'acier comprimé
=(30*0.225)/(1*0.362*161.304)
= 0.323
=> α1=0.455
= (1*0.36*0.4552 )/(30*(1-0.455))
2
= 45.58 cm ====> CHOIX 10 HA25
f/ Vérification ELU:
Lf = 0.7*L0 = 0.7*17 = 11.9 m
ea= max( 2cm ; L/250)
= max(2cm ; 1700/250) =6.8 cm
Pgousset =(0.32*1.034*2500)/2=413.6 dan/m
Pu 1.35.(Pc Et Pt+ Pgousset) 1.5.Q Pu=2865.35 daN/m²
− 𝑹𝟐∗𝑷𝒖 𝑹𝟏 14
N avec R2 = = = 32.8 m
𝟏+𝐜𝐨𝐬 𝝋 𝐬𝐢𝐧 𝝋 sin 25.4
− 32.8∗2865 .35
N daN
1+cos 24.5
Nu= 1.35*[(25*0.7*0.55)+(25*1*0.24*17)+(20*17*0.02)+(1*1*17*2)]
+[58.411*sin(25.4)*10]
= 456.32 KN
Mu= Mmax*1.5
= 17.808*1.5 = 267.12 kdaN.m
e0=Mu/Nu
= 267.12/456.32 = 0.585 m
e1=e0+ea
=0.585+0.068=0.653 m
Lf/h < max(15 ; 20e1/h)
avec Lf/h =11.9/0.4=29.75
et max(15 ; 20*0.213/0.4)=32.65 => e2 n'existe pas
emoy= e1=0.653 m
Va = d-(h/2)=(0.9*40)-(40/2)=16 cm
Mua=( emoy + Va )* Nu
= (0.653+0.016)* 456.32 =371 KN.m
fbu=(0.85*fc28)/1.5=14.17MPa
=371*10-3/(1*0.362*14.17)=0.202
=0.8*(h/d)*(1-0.4*(h/d)
=0.8*(0.4/0.36)*(1-0.4*(0.4/0.36)
= 0.494
< =>Section partiellement comprimé
< =0.39 Pas d'acier comprimé
α=1.25*(1.(1-2 )0.5)=0.285
z=d*(1-0.4* α)= 0.36*(1-0.4*0.285)=0.32
Ast= Mser/(z*fsu)=371*10-3/(32*(400/1.15))
=33.33 cm²
A2= Ast-( Nu/fsu)=33.33.10-4-(456.32.10-3/(400/1.15))
= 20.21 cm²
c/ Condition de non fragilité C.N.F:
e0-0.455*d
* =3.53 cm²
e0-0.185*d
=>A2=max(45.58 ; 20.21 ; 3.53) = 45.58 cm2
===> CHOIX 10 HA25
5 - DIMENTIONNEMENT DU RADIER
a/ Ceinture inferieure:
Fr= T0*R1= 26.218*10*14.12 = 3701.98 KN
Ast= Fr/σst=3701.98.10-3/161.304 = 229.5 cm²
===> CHOIX 48 HA25
Soit une section de ( 1*0.3)
Perimetre =R1*2π = 87.96m
Poids propre de la ceinture inf : 0.3*[Link] = 659.7 KN
Surface : π* R1² = 615.75 m²
Poids propre du radier : 615.75*25.0.3 = 4618.125 KN
Poids propre du liquide:17*615.75*1.15*10 = 120379.125 KN
P. Béton de propreté : 0.05*615.75*22 = 677.325 KN
P de coupole+etanch+paroi+ceinture = 87.96*293.341 =
25802.27 KN
P. Béton de gousset: 4.13*87.96=363.8KN
Ptotal=la somme = 152136.545 KN
σ= Ptotal/S= 247.075 KN/m² = 0.247 MPa
σsol= 3.3 bars = 0.33 MPA > σ ==> verifié
σA=6*Ns/R22 = 6*293.341/32.8²=1.636 Kdan/m²
Mmax=0.11*R2*Ns = 1077.62 KN.m
Moment sur appui : 0.3* Mmax= 323.286 KN.m
Moment sur travée: 0.95*1077.62 = 1023.74 KN.m
=1.023 MN.m
Ferraillage sur travée
d=0.9*0.3 = 0.27 m
Mrb=0.5*15*0.582*0.27²*(1-(0.582/3)) = 0.256 MN.m
====> Asc 0 : on va augmenter le hauteur de radier
pour éliminer l'acier comprimé
Soit h=0.6 m ( 0.3 initialement )
Mrb=0.5*15*0.582*0.54²*(1-(0.582/3)) = 1.026 MN.m
Mrb> Mtra :
=(30*1.023)/(1*0.542*161.304)
= 0.652
=> α1=0.33
= (1*0.54*0.332 )/(30*(1-0.33))
2
= 29.25 cm ====> CHOIX 10 HA16
Ferraillage sur travée
=(30*0.323)/(1*0.542*161.304)
= 0.206
=> α1=0.12
= (1*0.54*0.122 )/(30*(1-0.12))
2
= 21.91 cm ====> CHOIX 12 HA16
6 - ANNEXE
CEINTURE SUP
CEINTURE INF
COUPE SUR LA COUPOLE
coupe vertical de réservoir
