AVERTISSEMENT

Ce document est le fruit d'un long travail approuvé par le jury de

soutenance et mis à disposition de l'ensemble de la
communauté universitaire élargie.

Il est soumis à la propriété intellectuelle de l'auteur. Ceci

implique une obligation de citation et de référencement lors de
l’utilisation de ce document.

D'autre part, toute contrefaçon, plagiat, reproduction illicite

encourt une poursuite pénale.

Contact : [email protected]

LIENS

Code de la Propriété Intellectuelle. articles L 122. 4

Code de la Propriété Intellectuelle. articles L 335.2- L 335.10
http://www.cfcopies.com/V2/leg/leg_droi.php
http://www.culture.gouv.fr/culture/infos-pratiques/droits/protection.htm
 INSTITUT SUPERIEUR DE GENIE MECANIOUE
ET PRODUCTIOUE
'r.f{ESE
Prêsentéeù

BIBLIOTHEUUTWIVERSITAIRE
_ -vErz

U N IVE RS ITE D E METZ

N'inv. ^3%tG?s
PAR
Cote
s/Su*/rt
ZoubirBELHOUCHET

'.,,
Loc hancrt^T;c-
P'our I'obtention du grade de :
DOCTEUR dE I'UNIVERSITÉ dE METZ

SPECIATITE . Mécanique
OPTION : GénieCivil

COMPORTEMENTET STABILITEDES COQUESCYLINDRIQUtrS

SOUSCHARGESDE TYPE "SILO"
.. ETUDE ANALYTIQUEET NUMERIQUE''

soutenuele zr Septembrer9g8 devani la commissioncomposéede:

J-F. Jullien Professeurà l'I.N.s.A de Lyon. (URGC). Rapporteur

G. Dhatt Professeurà I'I.N.S.Ade Rouen. Rapporteur
M. Potier-Ferry Professeurà l'Universitéde Metz. (LPMM). Exarninateur
J-C. Roth Professeurà I'Universitécle Metz. (LPivIivI). Directeur clethèse
A. Khelil Maître de conférences UHP Nancyl - IUT NB. Co-Directeurde thèse
J. M. Rotter Professeurà I'Universitéde Edinburgh. Examinateur
F. Mohri Maître de conférencesUHP Nancyl - IUT NB. Examinateur
Laboratoire de physique et rnécani<1ue
des matériaux U.R.A. C.N.R.S L2L5.
Université de Metz, ile du Saulcy,57045 ce<IexO1
 Le trauail d,e rechercheprésenté d,ansce mémoire a été réalisé a'u Laboratoire d,e
Physique et Mécaniquedes Matériaux (L.P.M.M) d,eI'Uniuersitéde Metz sous la d,irection
du Professeur Jean-ClaudeRoth.

Tout d,'abord,je tiens à eaprimer toute ma gratitude au Professeur J-C. Rotlt pour
m,'aaoir accue'illiau sein d,eson équipeet surtout pour m'auoir initié au traaa'ild,ereclterche.

J'adresse également rnesrernerciementsau Maître de ConférenceAbdelouahabKhetit

d,eI'[Iniaersité de Henri Poincaré, Nancy 1 pour l'intérêt qu'il a porté au présent traaail
et la précieuseo,ideà l'éIaborationde cettethèse.

Par ailleu,rs,je tiens aiuement à rernercier les ProfesseursJean-François Julien de

I'I.N.S.A de Lyon et Gouri Dh,att de I'I.N.S.A d,eRouen pour aaoir acceptéIa charge d,e
rapporter ce traaail d,ethèse. Ma grand,ereconnaissace0,uprofesseur M. Potier-Ferrey
d'aaoir acceptéd'être président de ce jury.

Je remercie égalernentJ. Michael Rotter Professeur de I'Uniaersi,téde Ed,inburghet

Faudil Mohri Maître de Conférencede l'Uniuersité de Henri Poincaré, Nancy 1 d'auoir
bien aoulu examiner cette thèse et je suis h,eureuxd'auoir I'honneur de les cornpter parrni
les rnernbres d,ujurg.

Toute rnon arnitié aa à mes arnis et collèguesdu L.P.M.M. pour leur soutien permanant.

Enfin, je tiens à eaprirner rna plus profonde gratitud,e à ma farnille qui rn'a toujours
soutenu et encouragé.
 Le traaail de la présenteth,èsese résume en ileuæparties. La première concerneI'étuile
analgti,queet numéri,queilu comportement ile la coquecyli,nilriquesousmi,sea d,escharges
sgmétri,ques combinées (pression hori,zontale et frottement en paroi,s) type si,lo Une
nouaelle rnéthode analyti,quedécri,uantle comportement élasti,quelinéaire d,ela coque est
d,éaeloppé à parti,r iles équationsil'équi,li,bre.
Nous aaons ai,nsiobtenu \esfonctions d,ed,is-
tributi,on des solli,citati,onset d,esiléplacementsiles parois de Ia coque. La modéIisation
ile la structure par élé,mentfini en uti,Iisant le logiciel Abuqus nolrs a permis iJe ualid,er
les enpressi,ons obtenues. Les comparaisonsont montré une bonne concordanceentre les
résultats des si,mulationsnumériqueset ceuaobtenusanalytiquementpar notre approche.
L'analgse ilu cornportementde Ia coquesous chargesnon uniformes complète cette
partie. L'effet ile ces charges,ilues à la aidange ercentrée des silos se trail,ui,tpar iles
moments de fl,eri,ontrès i,mportants.
La deuni,èmeparti,e est consacréeà I'analyseilu aoi,lernentiles coquescyli,ndriquesim-
parfai,tes. Une étude pararnétrique a été effectué,eà I'aiile du programme Abaqus pour
iléterminer l'i,nfl,uencerelatiue ile chaqueparamètre ( épaisseurile Ia coque, ri,gid,i,té d,es
raiili,sseurs... ) affectant la rési,stanceile Ia coque. Sur la base d,esrésultats obtenus,
une rnéthod,esemi,-analytiqueile calcul ile Ia résistanceultime au uoilem,entde la coquea
été établi,epour une coque rai,ilie et, isotrope. Les ili,fférentesconpara'i,sonsétablies entre
Ies résultats, iles forrnules de la littérature, d,eI'Eurocode 3 et, il,e notre approche, aaec
les résultats numéri,quesmontrent Ia ualid,itéile la méthodeproposée. Cependant, cette
contribution presnenteune améIiorati,ondesformules actuellesdonneéspar I'8C3.
Mots clés:
Stabilité des coques cyli,nilriques, éléments fi,nis, formulations d,es coques, charges type
usilo', imperfections géométri,ques,calcul non-linéai,re,effet de mentbrane, effet ile
fleæion.

Thi,s work can be surnmarized,i,n two parts. The first one ileal with a theoretical anil
numeri,calstuili,eof the behaai,orof cyli,ndricalshell submissi,ae to combi,nedsymmetrical
loaili,ng(hori,zontalpressureand,wall friction) silo type. A new analgti,calmethoil,d,éscri,bed
the li,near elasti,cbehauior of sh,ellis ileaeloped, by the equili,briumequations. In th.is wag
we haue get the fonction of stress d,i,stri,bution anil displacementof shell slei,n.Structure
moileling by fi,ni,teelementusi,ngAbaqussoftwareue ûre permi,tto aaliil,atethe erpress'ions
obtai,neil. The comparaisonhaue shows a gooil concorilancebetweenresults of nurneri,cal
si,mi,Iuti,onanil those obtaineil with our analyticalapproch.
The analysis of shell behaaiorsubjecteilto unsymmettri,callaoding cornpletethesepart.
The effect of theseload ilues te eccentricd,ischarge of the silo it self trunslateil by a much
i,mportanted fl,eæi,on mo ments.
The secondpart concernsthe warping analysi,sof imperfect cgli,nilricalsshells. Abaqus
progranLwas used again i,n stuily of th,einfluence of the parameters ( shell thi,knesse,sti,ffner
stiffness ...) affecting shell strength. About basisresults obtaineil, a neu) semi-analytical
rnethoil is ileaelopeilfor calculati,onof the ultimate stress resistance of the sti,ffened anil
i,sotropi,cshells. The di,fferentscomparaisonsmade betweenthe results, of literature for-
mulas, of ECS and, of our appraoch, with the nurnerical results, show a aali,dity of the
suggestmethoil. Eoweaer, this contri,butionpresent a i,mproui,ng for present formulas giuen
by EC?.
Key words:
Cylinilri,cal shells stabili,ty, finite element rnethoil, shells formulati,on, si,Io loailing type,
gemetri,cimperfecti,on,nonlinear corryutation, membraneeffec't,fl,eni,oneffect,.
 SOMMAIRE

INTRODUCTIONGENERALE.. ..........1

CHAPITR^EI: GENERALITESSUR LES SILOSCYLINDRIQUESET

ANALYSE THEORIQUE DE LA STABILITE

1.1. LES SILOS CYLINDRTQUES

( C ON C E P T rON , C H ARGES ET CAUSES DES STNTSTRES) .. ..... 4
1 . 1 . 1 I. n t r o d u c t i o .n. ......4
I.1.2. Flormeet dispositif dessiloscylindriques . . . . .5
.
géométriques"
I.I.2.L Lesdéfautsdeformes"imperfections ........6
résiduelles
L.1.2.2Lescontraintes liéesàla construction ......6
sur les paroisdessiloscylindriques...
1.1.3.Chargesâgissants .. ... ..6
e F P 2 2 - 6 3 0. . . .
1 h a r g esse l o n l a n o û nN
1 . 1 . 3 .C ........7
1,.L.3.2Chargesselonl'Eurocode1.... .....9

L.7.4.Chargesclimatiqueet sismique.... ...11

1.1.5.Causesdessinistres. .. . ".72
 L . 2 .THEORIE DE LA STABILITE DES COQUES CYLINDRIQUES
ET EQUATIONS GENERALES . ... 13

L.2.I.Introduction.. ""'13

7.2.2. Equations d'équilibre et de stabilité des coquescylindriques - . . ... .. 13

I.2.2.7 Equations de la stabilité non-linéaire ' ' ' ' 74

7.2.2.2Equationsde Ia stabilitétinéaire ...'. " " ' 19

L.2.3. Charges critiques de flambage des coques cylindriques. ' ' " ' " '22

1.2.3.7Flambement d'un cylindre sous pressionaxiale ' ' ' ' ' ' '23

1.2.3.2Flambement d'un cylindre simplementmaintenu

à ses deux extrémités et comprimé axialement ' ' " ' ' ' 23

initiales .. . .
7.2.4. Prise en compte des imperfectionsgéométriques . . -.. ...24

1.2.5. La méthode énergétiqueet le critère de stabilité " " ' " '25

I.2.6. Traitement du problème non-linéaire " " " '27

7 . 2 . 7 .C o n c l u s i o n """ " 28

DES COQUESCYLINDRIQUESSOUS
CHAPITRE II: CO]VIPORTE]VIENT
CHARGESDE TYPE ''SILO'' EN ELASTICITE LINEAIRE
2.1. Introduction 29

2.2. Lpproche analytique du calcul des coques cylindriques

sous charges axisymétriques de type "silo" 30

2.2.I. Formulation du Problème 30

2.2.2.Résolutionanalytique 32

2.2.3. Résultats de I'analyse analytique . ' 34

2.3. Approche numérique du calcul des coques cylindriques

sous charges de type "silo" 36

2.3.7. Cas de charges axisymétriques . 37

2.3.2. comparaisondesrésultatsanalytiqueet numérique.. ....38

2.3.3. Interprétation desrésultatsanalytiqueet numérique " " ' 41

2.4. Lnùyse numériquedu comportementdescoquescylindriques
sous charges dissymétriques 41
 2.4.1. Analyse 1.

" v i d a n g e e x c e n t r ésee l o n l e sr é s u l t a i sd e l a b a s e d e C h a r t r e s" . . . . . . . . 4 2

2 . 4 . 1 . LR é s u l t a t sd e l ' a n a l y s e1 .....44

2.4.t.2lnterprétation des résultatsde l'analyse1 .. ....58

2.4.2. Analyse 2.

" v i d a n g e e x c e n t r ésee l o n l a n o r m eN F P 2 2 - 6 3 0" . . ........60

2 . 4 . 2 . 1 R é s u l t a tdse l ' a n a l y s e2 .....61

2.4.2.2lnterprétation
d e sr é s u l t a t sd e I ' a n a l y s e ?. . ....65

2.4.3. Analyse 3.

" Remplissagecentré selon les résultats de la base de Chartres .. . 66

2 . 4 . 3 . LR é s u l t a t sd e I ' a n a l y s e3 .....66

2.4.4. Analyse 4.
t ' R e m p l i s s a geex c e n t r és e l o nl ' E u r o c o d e1 ...... .68
" ..
2 . 4 . 4 . 1 R é s u l t a tdse I ' a n a l y s e
4 .....69

d e sr é s u i t a t sd e l ' a n a l y s e 3e t 4 . . .
2.4.4.2lnterprétation ....70

2.5.Conclusion .,. ...73

CHAPITRE III: VOILEIVIENT DES COQUES CYLINDRIQUES SOUS

CHARGES DE TYPE SILO ET ETABLISEIVIENT D'UNE

FORMULE ANALYTIQUE

3.1. Introduction . ..74

3.2. Flambage des coques cylindriques sous compressionaxiaJe . . . . ..75

3.2.7. Analyse analytique de la stabilité des coquescylindriques

s e l o nT i m o s h e n k o
et Gere .......76

3.2.2. Analyse numérique de la stabilité des coquescylindriques

sous compressionaxiale (calcul linéaire) . ... .79

3.2.3. Comparaison des résultats numérique et analytique

deTimoshenko. .....81

3.3. Analyse numérique de la stabilité des coquescylindriques

sous compressiona-xiale(calcul non-linéaire) .. .....83

 3.3.1. Prise en compte desimperfectionsgéométriques.. ........83
3.3.2. Etapes de créationde la surfaceimparfaite .. ... . . ... ... . 84
3.3.3. Résultats numériques de I'analyse non-linéaire . ... . 85
3.4. Analyse de Ia stabilité des coquescylindriques raidies et non-raidies
sous chargesde type "Silo" (calcul non-linéaire) . . . ... .. .87
3.4.1. Rappel de la solutionanalytiqueseloni'Eurocode3 .... . . ... ... 87
3.4.2. Etude paramétrique par voie numérique (coque raidie et non-raidie) . .90
3.4.2.LChargesappliquées ....91
3.4.2.2Choix du maillageet des conditionsaux limites . .. . . . . .... .92

3.4.3. Résultats de I'analyseanalytique de I'Eurocode 3 . ... . . .. 93

3.4.4. Résultats numériques de l'étude paramétrique . . . .. 98

3.4.4.L Effets des différents paramètres . . .109

3.5. Etablissement d'une formule semi-analytique pour le calcul

des contraintescritiques ...... . 109

3.6. Comparaison des résultats et validation de I'approcheanalytique ....I27

C O N C L U S I O NG E N E R A L E . . ,,,,.I27

BIBLIOGRAPHIE
ANNEXES ......134

Annexe 1. Solution analytique de l'équation différentielle du quatrième ordre... 134

Ar'nexe 2. Résultats numériques du deuxième cas de vida.ngeexcentrée .. L37

Arrexe 3. Présentation de la.cellule "Cl' en vraie grandeure de Chartres . . . . . . .146

Annexe 4. Exemple de vériûcation de la résistance d'une coque cylindrique
raidie verticalement, soumise à des chargesde type silo . . . . ...L47
 II\TRODUCTIOI\ GENERALE

L'analyse des structures étudie I'organisation, le fonctionnement, la stabilité et la

résistance des constructions. Plus particulièrement, elle'examine I'agencement d.esdivers

composants, étudie le jeu des forces, l'équilibre et la déformation des constructions et

détermine ies formes et dimensions à donner aux éléments structuraux pour assurer

la rigidité et la stabilité d'ensemble, compte tenu d'un degré de sécurité convenable.

L'amélioration de la qualité des produits sidérurgiques,plats et longs conduit à l'utilisation

d'éléments particulièrement légers. Les structures à parois minces de type coque employées,

en génie civil, dans l'aérospatiale et dans I'aéronautique, sont en général intensivemeni rai-

dies [5][a6][57][62]... Ce qui a permis aux chercheursde simplifier le problème traité en

supposant que les raidisseurs sont tartinés, c'est -à dire que leurs propriétés sont réparties

sur la totalité de la surfacel et en menant alors des analysessur une coque orthotrope

équivalente. Les résultats obtenus par de telles analysessont en accord avec les résultats

expérimentaux. L'interêt de I'utilisation de ces coques raidies s'est étendu ensuite vers

d'autres types de constructions, telles que les constructions marines et certaines structures

du génie civil, dont les silos, auxquellesnous nous intérressonsprincipalement.

Les silos sont des ouvragesde stokagede matière, présentantun domaine d'application

très vaste. Le calcul de ces structures, est lié au mode d'action des charges agissant sur

les parois cylindriques. De nombreuseséiudes ont été consacréesà la détermination des

actions provoquées par ces matières contre la paroi cylindrique. Une bibliographie im-

portante sur ces études est disponible dans les articles et les livres (voir par exempie

[10][32][45][50]). Par contre peu de recherchesont été réalisées,concernant Ie comporte-

ment et I'instabilité de ces constructions vis-à-vis de la nature réelle des charges sollici-

tant la structure. Ces phénomènesdoivent être contrôlés pour assurer la frabilité de ces
Intoduction générale

constructions. Les quelques études théoriques et expérimentales effectuées sur le com-

portement des coques cylindriques raidies et non raidies ne traitent que des cas de charges

simples. Compte tenu de la complexité des équations d'équilibre, et pour simplifier les

calculs, certains auteurs ne considèrent pas les actions simultanées de la pression interne

et du frottement latéral de la matière ensilée contre la paroi.

Un autre problème particulier concernele comportement de la structure durant la vidange

excentréese rajoute. Cette dernière présentebeaucoup d'avantages,dans I'exploitation des

silos par rapport à la vidange centrée,mais, elle donne naissanceà des pressions quelcon-

ques, non-uniformes, sur toute ia hauteur. Les seulesrésuliats expérimentaux disponibles

sont les mesures expérimentales de la base de Chartres [3a]. L'étude du comportement

et de la stabilité de la structure dans ce cas de charges non uniformes changent, et les

équations seront plus compliquéesà résoudre.

Les objectifs de ce travail de recherche se résument en deux parties. La première

consiste à présenter une approche analytique et numérique concernantle comportement des

coques cylindriques parfaites, en élasticité linéaire, sous chargesde type "silo" (actions de

frottement et pression horizontale). La deuxième partie, est consacréà l'étude du voilement

des coques raidies et non raidies, tout en tenant compte des défauts géométriques.

Nous avons établi, dans le premier chapitre, une revue bibliographique concernant,

et d'autre
d'nne part, les charges de type "silo" donnéespar les règlementsactuels 124)1451,

part les formulations mathématiques sur la stabilité linéaire et non linéaire des coques

cylindriques.

I,e deuxième chapitre est consacréaux comportement élastique linéaire, des coques

cylindriques parfaites et non-raidies, sous I'effet d'actions de type "silo". Nous présentons

une approche analytique qui permet dtanalyser ce comportemeni, sous lteffet d'actions

axisymétriques. Une analyse numérique (en utilisant Abaqus) complète cette étude afin

d'aborder l'effet de la non-uniformité des pressions(effet membranaire et de flexion) sur Ie

comportement.
Intoduction générale

Dans le troisième chapitre nous présentonsune étude du voilement clescoques cylin-

driques raidies et non-raidies sous chargesa.xisymétriques.Cette étude est réalisée en trois

étapes présentant successivement:

1- le voiiement des coques cylindriques parfaites sous compressionaxiale;

2- une étude paramétrique par voie numérique (calcut non-linéaire), permettant

d''analyser i'influence et i'importance des caractéristiques géométriques d.e la structure

et des actions sur la contrainte criiique du voilement.

3- l'établissement d'une formule semi-analytique originale pour le caicul de cette con-

trainte, concernarrt les coquescylindriques imparfaites raid.iesverticalement, tout en tenant

compte des charges réelles existantes en pratique.

Plusieurs comparaisonsmontrent la validité de cette formule ainsi que son efficacité

dans Ie dimensionnement des coquescylindriques sous chargesuniformes d.etype ,,silo".

 ChapitreI

GE}IERALITESSURLES SILOSCYLI}TDRIQUES
ET AI\IALYSETHEORIQUEDE LA STABILITE

r..1. LES SrLOS CYLTNDRTQUES

(CONCEPTION,CHARGES ET CAUSESDES SINISTRES)

L.1.1. Introduction

On désigne par "SILO", tout accumulateur de matière, constitué d'une paroi à

génératrices verticales, terminé inférieurement par des parois inclinées ou par un fond
plat horizontal fig(1. 1).

Les silos en général sont installés en batterie avec r:n certain nombre de cellules de sections
variées, juxtaposées de grande hauteur, terminées à leur base pa,r des trémies de vidange
de forme conique, et fermées à leur partie supérieure par une calotte ou par un plancher
sur lequel sont installés les appareils de remplissage des cellules.

Les cellules de section circulaire sont quelque fois liées entre elles de façon à créer d'autres
cellules dites t' intercalaires " ou " as de carreautt pour une meilleure utilisation de ltespace
disponible.

Les dimensions courantes d'une cellule cylindrique, métallique sont de I'ordre de 4 à 8

mètres de diamètre environ, et de hauteur pouvant atteindre 20 à 30 mètres.
Cha.pi.treI Généralitéssur les siloscyli,ndriqueset analysethéorirlttede la stabilité.

Calotte

Raidisseur
annulaire Cellule

Raidisseur
Ceinture vertical

Trémie

Fig (1.1)

L.7.2. Forlne et dispositif des silos cylindriques [31][42][51][56]

De nombreusesdispositionsd'installation existent,parmis lesquelleson distingue:

- le cas de cellulesreposant sur un certain nombre de poieaux, habituellement de 4 à
- ' est donc renforcé
6. Le périmètre cle base clu cylinclre ( liaison virole trémie poteau )
obligatoirement par une ceinture.
- le cas de cellules poséesdirectement sur le sol par I'interrnédiaire d'une infrastructure
annulaire qui assurele répartition de la charge sur le sol'

conique
Les silos soni égalementmunis d'un orifice de vid.ange,excentréau lieu d'une trémie
cenirée pour les raisons suivantes:
- la vid.angeexcerrtréeest plus rapide que la vidange cenbrée(vérifiée expérirnen'talement

[51]).
Ch,apitreI Généralitéssur les silos cylindriqueset analysethéorique d,ela stabili,té.

- Pour certaines positions du matériel de manutention, la vidange excentrée s'avère plus

convenable.

Enfin, les parois des silos sont généralement réalisées :

- soit avec une tôle d'acier ondulée, raidie ( ou non ) par des montants verticaux.
- soit avec une tôle d'acier lisse, raidie ( ou non ) p* des raidisseurs annulaires et des
éléments verticaux.

7.L.2.L Les défauts de formes 66imperfections géométriquest' [47]

A la construction, la paroi d'un silo n'est jamais parfaite. I1 y a toujours des défauts,
tant dans la direction radiale que dans Ia rectitude des génératrices. Ces irrégulariiés
initiales dans la forme de la coque cylindrique ont une certaine influence sur la stabilité;
elles diminuent sa résistance au flambement.

1,.L.2.2Les contraintes résiduelles liées à la construction

Dans une coque raidie où les raidisseurssont fi.xéssur la tôle par soudure, il se crée à
I'entoure de chaque soudure, dès son refroidissement,un champ de contraintes résiduelles
dans la tôle et dans les raidisseurs.

1.1.3. Charges agissants sur les parois des silos cylindriques

La figure (1.2) présente les forces qui entrent en jeu dans une cellule de silo. Ces
forces comprennent:

* Une pression horizontale P"(z) sur les parois verticales;

+ Une pression verticale Pr(") agissant sur les parois de la cellule;

* Une pression tangentielle verticale

/(z) due aux frottements sur les parois.
La détermination des pressionsexercéespar les matières ensiléessur les parois d'un silo est
complexe. Ceci est lié aux, caractéristiques physiques, aux valeurs des angles de frottement
interne des matières ensilées, du type de parois utilisées et des dimensions de la cellule.
Parallèlement, de nombreuses études théoriques et expérimentales ont été menées pour
définir les sollicitations sur les parois et diverses hypothèses ont été avancéesafin de simpli-
 et analysethéorique de la stabzlité.
Chapitre I Généralitéssur les silos cylind,rr,ques

fier le problème. Parmi toutes cesétudes se trouvent les travaux de Janssen-Koenen(1895)

les études présentéespar
[32], Reimbert ( depuis 1943) [49][50][52],Caquot (1957) [15] et
ie groupe de normalisation P22-630 [45].
Nous n'avons pas l'intention d'entrer dans les détails du développement des équations qui
donnent les charges selon les différents auteurs, car, cela uous éloigne de notre domaine de
travail concernant le calcul de structure.
En revanche, dans la modélisation d.ela coque, nous utilisons les charges données par la
norrne NF P22-630 [45] et I'Eurocode 1 [2a].

plan moyen
deremplissage

z
référencepourle
calculde la cellule
[ ï].,

| -,;1
corpsde la cellule

| +llf+r{l

Géometrie
K/ hessions

@ Fis (1.2)

1.1.3.1 Charges selon la norme NF P22-630 [45]

(i) Pression verticale:

rh
P ,(" ) :1Zs$( z) ; to--- (1-1)
^l-r
Chapi,treI Générali,téssur les silos cylindriques et analgsethéorique de la stabilité.

(1-2)

où:
Poids volumique
6 Angle de frottement de la matière ensilée contre la paroi;
? Angle de frottement interne de la matière ensilée;
p cceffi.cientde frottement contre la paroi F: tan$: f (z)lp^Q);
C Ccefrcient de rugosité de la paroi;
À Rapport des chargeshorizontales verticales À : p"(z)lp";
Th : Rayon hydraulique r1r: AIL;
A Aire de Ia section droite de la cellule;
L Périmètre de la section droite de Ia cellule.
ttr" et
"v" représentent respectivement l'état de remplissageet l'état de vidange.

ô(,) :1 - exp(-:) (1 - B)
\ zo,/

" z" est la côte de profondeur, ayant pour origine le plan moyen de remplissage.
L'angle B est déduit de la formule suivante:

sin2(6,,,)t cos(6,,,)1fsinz(g,,,)- sin2(6,.,,)

co s2 B- 11-4\
sin(p.,r)

Où les signes à adopter sont :

le f pour le remplissageet le - pour la vidange àvêcpa: 0.85g,.
(ii) Pression horizontale:

P"(t): ÀP"(z) (1-5)

(iii) Action de frottement:

f (r): pP"(z) (1-6)

Chapitre I Généralitéssur les silos cylinilriques et analysethéorique de la stabilité. I

(iiii)'Somme des actions de frottement

1z
F(z):lf(z)dz:.y.rhlz-zs.QQ)l ( 1- 7 )
Jo

Les critères de classifi.cationdes modes de vidange sont déterminésen fonction du rapport

(cellule pleine) :
I

, ho-ht
f:T

Ori :
h, : est la hauteur de la zone "morte"'
h1 : est Ia hauteur de Ia trémie;
h" : est Ia hauteur de la partie cylindrique.
k < 3 on adopte les formules de la vidange en masse.

* a f a 3 on adopte les formules de la vidange mixte.

k> 3 on adopte les formules de la vidange en cb.eminée.
Dans le cas de vidange on utilise les formules de base (paragraph 1.1.3.1) avec les
paramètres go et 6, et avec Ie paramètre ) modifié selon :

^: \M; Vidangeen masse

. - (2À'"Àr)
À Vidangemixte
;ffi; A r )
\Aa t

.l:
"- ,-(1À'À1).: vidangeen cheminée
(3À,*À,")'

où:
À. et À, sont calculésà partir de laformute (1-2) avecles paramètresrespectifs(gr,6,) et
(?",6")

L.L.3.2 Charges selon I'Eurocode L lza)

Selon I'Eurocodel, et dans le cas de remplissage centré les formules donnant les
pressions dans les silos n'ont pas changé, à I'exception des valeurs de À et du ccefficientde
frottement p,. Le tableau ci dessousrésume certainesvaleurscaractéristiquesÀ, pl adoptées
par I'Eurocode 1.
Chapitre I Généralitéssur les silos cglind,riqueset analysethéorique d,eIa stabilité. 10

Matière d'ensilage (r)kN/m3 À) Itt Co

Bté 8.5 0.55 0.35 1.35
Ciment 16.0 0.5 0.4 L.40
Clinker de ciment 18.0 0.45 0.45 1.40
Farine 7.0 0.40 0.3 1.45
Mal's d.c 0.50 0.3 1.40
Sable 16 0.45 0.4 1.40
Sucre 9.5 0.50 0.45 1.40

tableau (1-1) Paramètrescaractéristiquesde certainesmatières ensilées

Par contre en vidange centrée les pressionschangent et deviennent :

P..,(t): Cn.P"(z) (1-8)

f " ( r ) : C t . fQ ) (1-e)

Cn el C1 sont des ccefficientsd'augmentation des chargesdo''r'és par:

Cn:Co; Cl:I.1

Les valeurs de Co sont donnéesdans le tableau (1-1).

La figure (1.3) représentela variations des actions, dite type "silo" en fonction de la hauteur
selon les norines pour une cellule cylindrique de quatorze mètres de hauteur, quatre mètres
de diamètre et quatre millimètres d'épaisseur(la cellule C).
Ch,apitreI Généralitéssur les silos cylind,rigueset analysethéori,quede la stabilité. 11

+ prsi6 horiænrale (Ræplissa6c) + græion hciænnle (Roplisge)

-o- liDû@dl(K@Ptste) -O- Éo6@n! (Rofiplbûego)
-{F- grsim horiæ[Éle -+l- pruion horiæntale (vidæge)
(vidaoge æ mæe)
{}- &oaoot (vidæge)
<l-- Ëoaoæt (vidate a Easse)

l2
. Ùfarika æiléo (8Il) .lt"tU" *ife" (BLÉ)
I Cbargæ selon la N F W4n . Chrgæ selm lBreodo I

r0 t0

8 I

o
6 6

4 1

0 -0.1812l.3m24â!2
051015û25

Pressions [kNrho2] ke.ssions [kNrTnf]

Fig (1.3) Pressionh,orizontaleet action d,efrottement d,ansIa cellule "C"

selon la norme NF P22-630 et l'Eurocod,el(oidangeet' remplissagecentré)

L.L.4. Charges climatique et sismique

On disiingue deux grands groupes:

- Le premier est celui des silos isolés ayant une seule cellule cylindrique' potrr lesquels
on peut obtenir facilement la répartition des actions du vent en appliquant les règles N.V
65 [43];

- Le deuxième est celui d.essilos à deux ou plusieurs cellulesjuxtaposées, dans ce cas la

direction du vent et I'espacemententre les cylindres jouent un rôle important. 11n'existe
pas, à l'heure actuelle de règles pouvarlt donner directement la répartition des actions du
vent sur ce type de silo.

La stabilité de ia coque doit être vérifiée dans les deux cas suivants:

- Le cas où la charge du vent agit en même temps et dans le même sens que la charge
due à la vidange excentrée;

- Le cas où la charge d,u vent agit seule sur le cylindre vide. En effet, la charge du vent
 Chapitre I Générq'Iitéssur les silos cylind,rigueset analyse théorique de la. stabitité. 12

peut provoquer un claquage généraiiséde ia paroi [42]. Les raidisseurs annulaires doivent
pouvoir I'empêcher.

On distingue également :

- La pression exercéepar une cellule intercal.airesi de telles cellules existent;

- La sollicitation d'origine thermique, au cas où le p'roduit ensilé a une température

différente de la température ambiante (Clinker);

- La sollicitation d'origine sismique, si les silos sont construits sur des sites exposésà ce
phénomène.

1.1.5. Causes des sinistres

Dans les documents 118][42], des exemples de sinistres intervenus sur des silos
métalliques, ont été présentés. Les causessont diverses,on distingue deu,xgroupes princi-
parl\:

- Les désordres dus à une mauvaise détermination des pressionsdes matières d'ensilage
sur les parois, ou à une sous-estimation des charges du vent;

- Les désordresdus à un mauvais dimensionnementou à des erreurs de conception. Faute

d'une méthode ou de formule appropriée pour le calcul de la stabilité des silos, on emploie
des formules basées sur des hypothèses différentes de ce qui se passe réellement dans les
silos.

En effet les formules développéesdans la litterature pour les cylindres non raidis ([20][22])
et les cylindres raidis [29] considèrent une charge axiale constante et ne prennenr pas en
compte directement la variation des actions (pressionshorizontaleset action de frottement)
en fonction de la profondeur. Le cas de la distribution non uniforme des pressions en
vidange excentrée n'est pas traité entièrement
Chapi,treI Généralitéssur les silos cylindriqueset analysethéorique de la stabilité. 13

1.2. THEORIE DE LA STABILITE DES COQUES CYLINDRIQUES

ET EQUATIONS GENERALES

L.z.L. fntroduction

Lorsqu'un chargement appliqué a pour effet de vaincre " ia rigidité de forme " d'une
structure, celle-ci ne résiste plus que par sa rigidité de fléxion. C'est ie cadre de la théorie
des coques flexibles ou encore des coques souples. En général, on observe de grands
déplacements qui peuvent apparaitre brusquement. C'est le voilement (on dit parfoit
flambage). Le changement de la configuration de la structure peut être graduel, dans ce
cas le point de flambage réel est plutôt arbitraire. Cependant, dans la majorité des cas
(définis par différents types de structures ou différents cas de charges), ce point est net et
facile à identifier.
Après flambage, la résistancede Ia structure peut d.iminuer. Ce phénomène d.épenddu cas
de charge appliqué, du type et de la géométrie de la structure, des niveaux de contraintes
atteints, etc...
Pour les parois minces, il est devenu impératif de prendre en compte ce phénomène dans
le dimensionnement et la vérification des structures. Un cadre simplifié d'étude de la
stabilité des structures à été proposé par Euler pour des colonnes. C'est I'extension aux
coques minces de cette même démarche que nous abordons dans ce chapitre.

L.2.2. Equations d'équilibre et de stabilité des coques cylindriques

Une coque est un solide qui occupe dans ltespaceun volume compris entre deux sur-
faces telles que leur distance soit petite devant les autres dimensions. Cette caractéristique
géométrique permet des simplifications dans les modèles tridimentiorrnels. On peut les
obtenir de différentes manières (hypothèses a priori ou développements asymptotiques).
Mais dans chaque cas , nous sommes conduits aux relations cinématiques de Kirchhoff-
Love [41]. Elles expriment que la normale matérielle à la surface moyenne de la coque ne
subit pas de dilatation et reste, en tant que normale matérielle, orthogonale à la surface
moyenne déformée. En effet, ces conditions impliquent que chaque tranche tranwersale
de matière, se comporte commme un solide rigide, son mouvement étant décrit par un
torseur. On peut ainsi décrire toute la cinématique du milieu tridimentionnel à partir de
grandeurs définies sur la surface moyenne de la coque. En I'occurence, le mouvement se
déduit d'un champ de translations et de rotations. Evidemment, ces hypothèses ne sont
pas toujours compatibles avec les lois de comportement du milieu et on ne devra retenir
Chapitre I Généralitéssur les silos cylind,riqueset analgseth,éorigued,eIo,stabilité. L4

(parmi ces relations de comportement), que celles qui sont en accord avec les équations
d'équilibre. Dans la pratique, le modèle de Kirchhoff-Love [37][41]est relativement délicat
à mettre en æuvre. On préfère souvent faire appel à des versions simplifi.ées qui diffèrent
entre elles sur la façon dont la courbure est prise en compte (modète de Koiter [Sf]' ae
Novozhilov et Donnel [2f], etc...). Dans tous les cas, l'équilibre de la coque est régi par
deux phénomènes:

- l'effet de flexion

- l'effet de membrane

Le premier se manifeste dans un mouvement de rotation de la surface moyenne, tandis que

Ie secondfait intervenir les déformations propres de cette surface.

Les équations des coquescylindriques que nous présentonsdans ce chapitre sont baséessur
la théorie et les simplifications de'Donnell [21]. Elles sont relativement peu compliquées;
elles donnent cependant des résultats très satisfaisantssurtout lorsqu'elles sont employées
dans les problèmes de flambement des panneaux cylindriques relativement plats ou des
cylindres peu profonds (quasi- shallow shells). Pour cette raison, on trouve les équations
de Donnell à la base de la plupart des analysesde stabilité publiées.

Les équations non linéaires d'équiiibre de la coque cylindrique peuvent être obtenues soit :

- directement par la sommation des forces et des moments agissant sur un élément de la
coque cylindrique dans une configuration déformée.

- en déveioppant les équations pour ies coques de forme générale et en déduisant les
équations d'équilibre pour une coque cylindrique comme un cas particulier.

L.2.2.L Equations de la stabilité non-linéaire [6][13][14][36]

La figure(1-4) représenteun élément de coque cylindrique en configuration déformée'

Soit AP" f incrément de la force axiale par unité de surface et P, Ia force radiale par
unité de surface. Ecrivons les équations d'équilibre statique des forces et des moments
agissants sur cet élément. En écrivant successivementque la somme des composantes des
forces dans chacune des directions x, y et z est nulle, et en utilisant les hypothèses de
Donnell (cytindre peut profond, le matériau est homogèneet isotrope), nous aboutissons
aux équations suiva,ntes:
Chapitre I Généralitéssur les silos cylind,riqueset analysethéorique de Ia stabilité.

Suivant z on a:

N","*Nu",r*APr:g (1- 10)

Q u , u : M r , o o- M r v , " v ' (1- 11)

,1
I

qÉAQv/a/dY

N"rrâN'y'ôydy
-311fr616na*ravl
av
NrrâNyô dv

..4*aQ,n"a, N dz
'ùv[àz+
Nz"tôN/dz dz
AftrAwD') d' '\"

Z1
I
I

Fis (1.1)
Chapi.treI Générali'téssur les silos cylind,riqueset analysethéorique ile la stabilité. 16

Suivant y on a:

N o , o* N " s , r : 0 ( 1 - 12)

Q " , " : M " , " " - M y r , n ,, (1- 13)

Suivant x on a:

Q " , , * Q o , o* i * * N r t n , "*, N r , " u , z* N o w , r*s N r , r w , y

-Pr* N r y u ) , " g* N y " , " I ! , y * N y r u , y " { Nnrta,"a LP"w,r:0 (1 - 14)

avec :
ô( )_,,.
-=( a( )_,
a" ),,i fr:{ ),t
En introduisant les équations(1-10) et (1-12) dans l'équation (1-14) on oboutit à l'équation
suivante:

Q " , " * Q n , o* i " , * N " n ) , " r * N y w , y s * N r y w , " s* N s " w , ' s r : P , (1 - 15)

Par ailleurs, les intensités des forces et des moments sont liés aux contraintes internes par
les relations suivantes :

N" : : c(e"*ues)i Ny--

* î) d,a (r .î)dft : c(es+ ue,)
Ilr",(t I:ro,

!,{,!: Ny": r,o(r+ l)a"

Iu ,
Chapitre I Généralitéssur les silos cylindriqueset analyseth,éoriquede Ia stabilité. 17

1i
11[ : - M o r :
"e
,"n(t + i)"a,
{|- z,

- 1ï / f.\ rt /-.0\,
Q,: - Qy: l_rr"ult+=)dx
J| - i . r , n (\t + : ) rd/ r ; J-2

c21 Tzst...etc. étant les valeurs des contraintesà un poirtt quelconquedans l'épaisseurde
la coque ( à une distance x de ia surface moyenne).
Parmis certains hypothèses de Donnell, le rapport i p"nt être négiigé devant l'unité. Et
dans ce cas l'équation (1-15) s'écrit sous la forme suivante :

/t \
M","r-2Mrs,rs* Mn,rrl N 2 u ) , z z a 2 N r n w , " n * N o( ; * . , 0 0 ) : r, (1 - 16)

Les relations entre les d.éformationse; de la coque et les efforts N, par unité d.elongueur
ainsi que les relations entre les courbures K; de la coque et les moments IVI; par unité de
longueur sont données par les relations suivantes:

N": C(t, + r€y) ; M":D(K,*uKu)

Ns : C(e,1- ue") j Mc:D(Ku*uK,) (l - rrj
l a r È \

1-rt 1-rt
It
1 Vz ! : U--7ry i 1vI"!: Oï*"n
z

Où C et D sont réspectivement la rigidité membranaire et la rigidité de flexion données

par les expression suivantes:

C_
(L - ,')
n:ffi-,,1
Les déformations e; et les courbures K; en un point de la coque sont défi.niesen fonction
des déplacements u, v et w à partir des relations de Donnell [21] :

1
e":'Il,z + :r(r,")2 ; K": -11),zz

w1 (1- 18)
Ey : u,s - 1 L, ; Kv: -u,ua
)(*,ù',
'Yrg:U,U-U,ztUrzUrg
I l{ : -u,29
"v
Chapitre I Généralités sur les si,Ioscylinilriques et analyse th,éoriquede lo, stabilité. 18

L'introduction des expressionsdesmoments équation (1-17) et les expressionsde courbures

équation (1-18) dans l'équation (1-16) donne l'équation suivante:

- D ( w , " " " " ] - 2 w , z z s ty r , r n r r ) | N 2 t o , z! z2 N 2 y u , z u * Nn(l - .,rr) (1 - 19)

En résumé, on a donc les trois équations d'équilibre suivântes :

Nr,"*No",r*AP":g
N o , ro N " s , " : o
( 1 - 20)
DVaw- a 2N,rw,,o
* Nu(l - ,,no) *Pr:0
l*,t,""
Les termes Nr, N, ei N", qui figurent dans les équations d'équilibre ci-dessus sont des
fonctions inconnues de z et d.ey. Ils représententà la fois les forces membranaires variables
qui sont dues à I'extension de la surface moyenne, produite par la déflexion et les forces
dues aux charges appliquées.
Nous avons ainsi obtenu trois équations pour quatre inconnues. Il nous faut encore une
quatrième équation qui fait intervenir la déformation du système,autrement dit, l'équation
de compatibilité de déformations.
L'introduction des relations (1-17) et (1-18) dans les équationsd'équiiibre (1-20) conduit
aux équations d'équilibre développéessuivantes :

!-u .I+u ( 7-u z\, (t+u \ AP"

,t,zz*Tu,ur*-Tr,"r*rts,r(.,r' t, .,n, - .,ru :0
; )1.,n \, )+i
(1 - 21a)
!-u l*u / |-u z\ ([lu
o , u u i -- f u , " r * * t,r, \ :U
Tu,"sIr,s (r,r, 7.,""-; )*w,, \. , )
(1 - 21b)
. | / u\ / 1\l
D Y a u- -
* ( u , r * u , , ( ! u ) w , , u ) (+r , r , I u , , , z z . ; ) l
" lu,,(.,",*uu,ss*i)

(u,,-î) . \ * p.:0
(+ .ry) @?,+,.?,)
@?"*,uzy)- (1- 21c)
Ainsi, nous avons trois équations couplées avec trois inconnues u, v et w. Elles sont
suffisantes pour analyser le comportement au flambage et au post-flambage de la coque
cylindrique.
Chapitre I Généralitéssur les silos cyli,ndriqueset analgsethéori,quede Ia stabilité. 19

Dans le cas où la compressionæcialeappliquée est constantele long du cylindre "L,P, :

0". il est possible dans ce cas d'introduire une fonction d'Airy F et d'obtenir ainsi deu-x
équations au lieu de trois avec deux inconnues F et w :

1F,""(l- ,ur)l * a : o
,^. - *lr.r7)),zz-2l,zy'tD,,y
f (r - 22)

I "(r:,,-'t.,zzu,!!'-t,,",):YaF
La fonction d'Airy est supposéesatisfaire les relations suivante:

N,: (tl2)F,ur; Ne : (tl2)F,,"; t{"u : -(tlz)F,,s

Les équations (1-22) sont appelées"les équations de Karman-Donnell".

Remarques:
- Contrairement à ce que I'on observe dans l'étude des plaques et des poutres, même
la branche fondamentale de pré-É,ambageest non-linéaire, pour une coque cylindrique-
Cela vient du fait que, même si le cylindre est parfait, les déformations de pré-flambage
entrainent des rotations des éléments structuraux;
- La solution d.es équations différentielles non - linéaires de l'équilibre donne les
branches d'équilibre primaire et secondaire. Les branches d'équilibre montrent directe-
ment I'existence du point de bifurcation et la valeur de la charge correspondante;
- On ne peut cependant d.éterminer le point de bifurcation que par la résolution des
équations différentielles linéaires que nous alions établir dans le paragraphe suivant.

L.2.2.2 Equations de la stabilité linéaire

La charge correspondante au point de bifurcation pour la coque cylindrique soumise

partir
à une compression a-xiale et une pression latérale conservative est déterminée à
des équations qui sont présentéespar la suite. Ces équations sont obtenues à partir des
équations non - linéaires (1-20) par l'application de Ia méthode dite de l'équilibre adjacent.
Cette méthode est basée sur le fait qu'au point critique, il existe pour la même charge,
d.esconfigurations d'équilibre infiniment voisines à l'état fondarnental. On donne aussi aux
variables de déplacement un petit incrémént :

u:uo+ut
u:uo*ut (1 - 23)
w:wo*wt
Chapitre I Généralitéssur les silos cylind,riqueset analysethéorique de Ia stabilité. 20

(uo, ,o et urs) représentent les composarrtesdu champ cle déplacement correspond.antà la

configuration fondamentale d'équilibre et (u, u et w) une configuration adjacente d'équilibre
pour la même valeur de la charge appliquée;
(ur, ,, et u1) sont les composantes du vecteur représentant les petits incréments. Les
contrainter (If", l/, /f"r) des équations (1-20) prennent la forme suivante :
"t

N,-1r,"+Aff,
rtry-lrrr+aN, ( 1 - 1 -4-))l
\-
Nr!: N r o r o+ A N r s

ff"', Nro et Nroro correspondent aux déplacementsu9, ue et ?r.rsI

Altrr, A.f[, et 6Nzu correspondent aux incréments ur, er et tl1.
En désignant par N"r, Ny, et.f[,,r, les termes de A.lf,, AN, et 6Nzu qui sont linéaires
€n u1, u1 et w1 et en considérantles équations(1-17), les coniraintesgénéraliséesN", Ny
et Nrn peuvent se mettre sous les'formes suivantes :

N, : c ([u,"+ - Y+
ulu,u
]{r,,)') t f;fr,r>'l)
Itrv: c (lu,o- + vlu,"+ }rr,,l'l) ( 1 - 25)
i |{.,r)') *
N,u : c|{rg * t ) , 2* w , , w , y )

En introduisant les équations(1-23) et (7-21dans les équationsnonlinéaires (1-20), tous

les termes €o us, uot 1!0 et P disparaissentparce eu€ us, 1)0,'tlo décrivent la configuration
d'équilibre. Et si en plus, nous négligeonsles termes d'ordre quadratiques et ceu,xd'ordre
supérieurs €rr u1, uLt lDL du fait qu'ils sont très petits nous arrivons aux équations de
stabilité suivantes :

N"r,, * Ng"r,, :0
N a r , s* Nrur,r:0
(
DYawl - (N"owt,r, * N"rwo,r") l2(N"yow1,ry I Nryrwo,"y) | Nnoutl,yy*
I
(

*.(l*-0,r,)
):o-
(1 26)
Chapitre I Généralitéssur les silos cvlinil,riqueset analvsethéorique d'eIa stabilité. 2I

Les équations de stabilité (1-26) permettent de déierminer la charge critique pour un com-
portement linéaire (fig.t" 1.4b) ou non-linéaire (figure 1.4a) avant flambage. Cependant, la
présence des rotations de pré - fla,mbagedans ces équationsde stabi[té introduit une com-
plication substantielle dans les calculs. Heureusement, dans un grand nombre de problèmes
d'instabilités (flambage des coquesminces parfaites), f influence de ces rotations est petite
et négligeable.

(a) (b)

Fig (t./)
leséquations(1-26)deviennent:
de pré-flambage
Lorsqu'onnégligelesrotationsr.!s,2,u.rs,,

N"'' * Ns"'Y- 0
f
/rr r,y I N"sr," : Q
| (t - 27)
,o*-, - (",,r, ,", *2N,uour,,s* r[g,our,y,
* -0
| l"r,)

Les équations (1-27) sont les équations de stabilité dans iesquellesla non-linéarité de pré-
flambage est négligée (voir figure 1.4b).
Les équations (t-27) sous forme de déplacements en utilisant les reiations cinématiques
sont données par : .

ur ,""*
7-u
+
L+u - u : ^u
Tur ,ou Tut,r y iu7,z
ut,yy*
L-u
* -\ Tut,"y
+u - 1 : ^u
n- ut;z iuL,z
(1 t
-Y \ - , ^'
D Y a w l- + ,ur,r) (N"our, ," + 2N rsotrr,rs* Nnowl,nn) : 0
7 (r',0
(1 - 28)
Chapitre I Généralitéssur les silos cylindriqueset analgsethéorique de la stabilité. 22

Ces équations peuvent être découpléesd'aprés les hypothèsesde Donnell. On obtient après
a.rrarrgementles trois équations découpléesde stabilité de Donnell suivantes :

u1-o -Ul
-r'lDl,zzz- r : V
iat,OO
2+u
'ftt,""y :l
* Vnut (1 - 29)
itt,ron:
EJrr,r"""
DV8w1- V4 (Nrowt,r,| 2 N r r o w l , r y* N y o w r , 3 ,*y ) : 0
I

En posant N"so :0, on peut utiliser ces équationspour calculer la charge critique d'un
cylindre soumis à une compressionaxiale et une pression latérale.
En posant N"so : 0 et .lfro : 0, on a le cas d'un cylindre soumis à une compression axiale.
Puisque f influence desrotations d'avant-flambageest négligéedans ces dernièreséquations,
les cæffi.cients(Nro, trg,o et Nroo)'sont gouvernéspar les équations linéaires d'équilibre
suivantes :

( N"o,' * Nu"o* AP' : g

ffro,y * rf"so : o (1 - 30)
{
t1
DV*tor - * P,: g
I ;*u"

Le système d'équations de stabilité possèdeune solution non triviale seulement pour cer-
taines vaJeurs de chargesappliquées. Ces chargesreprésententles valeurs caractéristiques
du système, et à chacunecorrespond un point de bifurcation.. La charge critique correspond
aJors à la plus petite de ces valeurs.

L.2.3. Charges critiques de flambage des coques cylindriques [19]

Ce paragraphe, résume la formulation des charges critiques "Q"" d" flambage des
coques cylindriques. elle est valable uniquement pour les deux modèles suivants :
Modèle 1. Epaisseur très faible devant les rayons de courbure, (les effets de rotation autour
de la normal à la surface moyenne de la coque sont négligés).
Modèle 2. Epaisseur très très faible devant les rayons de courbure. Les déplacements
tangents de la surface moyenne étant eux-mêmes négligeablesdevant les rayons de courbure
de la surface moyenne, comparés au gradient de la flèche. Ce modèle est adapté atr.\ coques
dites peu profondes, c'est-à-dire ressemblant à une plaque. Plusieurs auteurs ont proposé
ce type de modèles, parmi lesquelsNovozhilov et Donnell.
Chapitre I Généralitéssur les silos cylindriqueset analysethéorique de Ia stabilité. 23

f.Z.g.f Flambement dtun cylindre sous pression axiale

En utilisant les notations de la figure (1.6) nous obtenons : (da,nsle cadre des modèles
tet2)
(i) Cas des culindresencasfuét
lzbteê

,n ( aru"\2
-e"1#^('.+) J o , \ ( d f 1 ) '/

rw ( 1- 3 1 )

Fis (1.6)
où:
^ 72(7- u2\
a_:È: ( 1- 3 2 )

En choisissant u, - 1 - cos(z'(l lzil) on déduit:

-e"1ffi('.#) 12
T
( 1 - 33)

(ii) Cas des cyli'ndresencastrés-encastrés

Avec u" - 1 - cos(2r(|lH), la charge critique sera :

-e"1ffi('.#)+ (1-34)
1.Z,B,Z Flambement dtun cylindre simplement maintenu à ses deux extrémités
et comprimé axialement
Chapi,treI Généralitéssur les silos cylindriqueset analysethéorique ile la stabilité. 24

t u plrrs petite charge critique dans ce cas est donnée par I'expression suivante :

E*
tffi
( 1 - 35)

Q" est un effort résultant exercé sur ies couronnes aux extrémités du cylindre dont
l'épaisseur est t. La contrainte critique associéeest donc:

ET
0c: (1- 36)

!.2.4. Prise en compte des imperfections géométriques initiales[4][30][47].

Il existe un désaccordimportant entre la valeur de la charge critique classique "P"1"

et la charge critique expérimentale "P"r', figure (1.7). Or, ni la prise en compte des
déformations de pré-flambage, ni I'influence des conditions aux limites ne peuvent seules
donner une interprétation satisfaisante de i'écart entre I'expérienceet la théorie classique.
En faite, ce sont les imperfections géométriquesinitiales de la coque qui sont principalement
à l'origine de cet écart. Il est nécessairede les prendre en compte [6][36].

Fis (1.7)
Les équations de Karma4:fonnell (L.22) ne sont applicables qu'aux cylindres parfaits.
Pour tenir compte des imperfections intiaJes dans ces équations, quelques modifications
sont à faire :
- le déplacement latéral est composé d'une flèche réelle "ur" produite par la charge ap-
pliquée, et une distorsion intia.le n111
t1<1'correspondantaux imperfections géométriques
initiales.
 d, Io:tùj!ilé'
Ch,api,treI Généralités sur les silos cylind,riqueset analgseth'éoriqu" . ë

- ies premiers termes qui représententles forces de cisaillement dans la première équation
ne subissent aucune
(1.22) dépend.ant seulement de la courbure produite par la flexion,
par pat "r'D*w * " dans
influence d.ecette distorsion. Par contre, on d.oit remplacer "u)"
les autres termes de l'équation"
les équations différentielles gou-
Les équations d.e Karman-Donnell ainsi modifrées seront
initialement imparfaites, et
vernantes pour l,étude de la stabilité des coques cylindriques
se mettent sous la forme suivante :

Y a w - L , l r,n o {r,,, * 7n*,,r -) 2,F,,,( ' ,r o + - *,r v) l

-*l',,"(l
*', ,9! *.*,no)] * * 0 (1-37)

n
r .. * w*ro)-
* w*,"v)2- ('"" + ror*",)(uau
I ''(w,",
* .*,"")fv 4 F
T
l@,r,

1.2.5. La méthode énergétique et le critère de stabilité [19][40][64]

les paragraphes précédents'

Les équations d'équilibre et de stabilité développéesdans
système'
peuvent aussi s'obtenir en analysant l'énergie potentielle du
fI est la somme de son énergie de
Pour une coque cylind.rique, l'énergie potentieile totale
extérieures appliquées'
déformation u et de l'énergie potentielle o des charges

II:[/+O (1-38)

donnéepar la relation :
L'énergied.edéformationu d'un milieu continu et isotropeest

( 1- 3 e )
u:, @,e; +orer-ro)e" +râ'Y,sar*tv,+r;"Y;)dædvdz
| | |
conformément aux hy-
En négligeant les termes d,ecisaillement et de contrainte normale,
les relations constitutives
pothèses ad,optéesd.ansla théorie des coquesl et en introduisant
(1-17) dans l'équation précédante on obtient:

+|t;,')a,aua" (1 - 40)
u:*
| | I (,r'*,;'+2ve"en
Chapitre I Généra|itéssur les silos cylinilriqueset analgseth.éoriquede la stabilité. 26

où:
( e'-e'*xK"
I
, ,i-er*æKu
)
I lâ : ^l'Y | 2xK 'n

e; : Déformation d'un point de la surface moyenne;

e; et 7, : représentent respectivement la déformation normaLeet de cisaillement d'un
point se trouvant à la distance "æ" de la surface moyenne de la coque.
En introduisant ces relations dans I'équation (1-40) et en faisant I'intégration sur "o" qui
varie de -t à t, on trouve i'expression suivante :

U:U*lUo r1- 41\

(r,n: **r?") a,a"

i Il (r, * e?+z,e,e"
u, : + + K?+ 2uK,K,+2(1- u)K!")dædz
| | @?
U* el,U6 Représentent respectivementl'énergie de déformation membranaire et de flexion.
L'énergie potentieiie des cha,rgesappliquéesest l'opposéedu travail fournis par ces charges
lorsque la structure se déforme :

n: - + Pw)d.ad,z (1- - D.\

ll Uu
'-
\ J

Où / et P sont respectivement les composantes axiale et latérale agissant sur la surface

d'un élément coque.
u, ur sont les déplacements d'un point de la surface moyennede la coque.
La substitution des expressionsde I'énergie de déformation [/ et de l'énergie potentielle Q
dans (1-38) donne une expressionde la forme suivante :

n-- [[FsdÆd,z
JJ

En utilisant le principe de I'énergiestationnaire[6II = 0] (l'integrant Fs doit satisfaireles

équationsd'Euler de calcul desvariations[19][16][40]).En introduisantles relations (l-22)
nous retrouvonsles mêmeséquationsd'équilibre(1-25).
 et ano'Iuseth'éoriuuede la stabilité' 27
chapitre I Générali,téssur les si,losculznd,rtgues

incrément virtuel de déplacement

La variation seconde de fI s'obtient , en ajoutant un
de (us' ?0' tl,o)'
(ur,ur,u.r1) sur la configuration d'équiiibre fondamental
Donc:

: + 2e"1esr
* e2n, +|4n,) o'ooo"
înn Z I I I (,t,
.i (N,or?,"*Nuou?,v dtda
*2N"row1|'uL,s)
ll
dxds (1- 43)
*2(1- v)(wl',,u)
| 2uwv"'w7'!a
. "U ('?,', * w?,vc
II
dans la relation précédante'
Les rotations d'avant-flambage uss et ?o6,sont négligées
celui de I'équation (1-43) conduisent aux
Les équations d,Euler avec coûIme intégrant,
le principe de l'équilibre adjacent'
mêmes équations (i-32) déjà trouvées en appiiquant

L.2.6. TYaitement du problème non-linéaire

suivantes :
L,analyse non-Iinéaire est nécessairedans les situations
En effet, la théorie linéarisée de la
1- Pour déterminer la branche d'équilibre secondaire.
du point de bifurcation' Mais elle
stabilité (paragraph L.2.4.) est limitée à la détermination
branche secondaire' or, la connaissance
ne donne aucun renseignementssur la forme de Ia
le comportement de ia structure'
d.ecette branche est importante pour mieux comprendre
sera stabie ( avec une résistance
c,est-à-dire de savoir si la structure aprés bifurcation
croissante ) ou instable'
d'avant-flarnbage et des imperfections
2- Pour tenir compte des déformations non-linéaires
exactitude'
initiales, et de déterminer la charge critique avec
des matériaux'
3- Pour tenir compte du comportement non-linéaire
(L-27) pour Ie cylindre parfait
La solution analytique exacte d.eséquations non-linéaires
pas. oo est généralement anrené à résoudre
ou (1-37) pour Ie cylind.re imparfait n'existe
leproblèmepardesméthodesapproximativesetnumériques.
de ces problèmes sont les suivantes
Les méthodes généralement utilisées d.ansla résoiution

-lesméthodesdedéveloppementenserie(RitzouGalerkin);

- la méthode des différencesfinies;

Chapitre I Généralités sur les silos cyli,ndriqueset analyse tltéorique de Ia stabilité. 28

- la méthode des intégrations numériquesl

- la méthode des éléments ffnis.

Toutes ces méthodes d'analyse numérique réduisent le problème non-linéaire à un problème

de résolution d'un système d'équations algébrique linéaire pour chaque cas de chargement.
Le système d'équations est alors résolu par I'une des méthodesincrémentales ou itératives.

L.2.7 Conclusion

Dans la première partie de ce chapitre on a présentéles chargesde type "silo" agissant

sur les parois cylindiques donnéespar les norrnes,NF P22-630et I'EC1. Nous avons montré
ainsi que les pressions alcialeset latérales, s'exerçant sur un silo cylindrique sont répa.rties
de façon non uniforme le long du cylindre et le long du périmètre.

Dans la deuxième partie on a pr-ésêntéles principales équations de la stabitité linéaire et

non-linéaire. Ce sont les équations de la stabilité de Donnellf2l], eui se résument par:

- les équations de stabilité dans lesquellesles termes de rotations d'avant-flambage sont

retenusl équations ( 1-26);
- les équations de stabilité dans lesquellesde tels termes omis sont données sous trois
formes équivalentes;équations (1-27), (1-28) et (1-29).

Le système des équations algébriques obtenu a été résolu numériquement. Cependant

aucune des études publiées n'a abouti à une solution analytique explicite à usage pratique
pour le dimensionnement du cylindre raidie. Toute fois on peut citer l'étude de Shang,
Marulic et Sturm[26] qui reste valable uniquement pour des cylindres raidis verticalement,
soumis à une compression axiale constante.

La suite de ce travail concerne l'étude anaiytique et numérique du comportement des

coques cylindriques non raidies d'une part et l'analyse de la résistance au voilement des
coques raidies et non raidies dtautre part.
 Chapitre II

DES COQUESCYLil\DRIQUES
COMPORTEME}TT
SOUSCHARGESDE TYPE SILO
E}TELASTICITELI}IEAIRE

2.1. Introduction

Cette partie a pour objectif d'étudier le comportement de ia structure des silos

6{lalliques cylindriques sous I'effet de la matière ensilée. On s'intéresse plus parii-
culièrement à la distribution des contraintes, des moments, des déformations et des
déplacements le long d.esparois d'une cellule à section droite circulaire. Dans la première
partie, nous proposons une solution originale analytique de ce problème en élasticité
linéaire, dans le cas de charges axisymétriques correspondant aux pressions de la matière
ensilée sur les parois du silo en pha.sede remplissage ou de vidange centrée. Dans Ie cas
de ]a vid"..ge excentrée, les pousséesdes grains sur les parois du silo sont dissymétriques.
L'influence de cette distribution non-uniforme, sur le comportement de la structure est
abordée numériquement (en utilisant Abaqus) dans la deuxième partie.
Chapitre II Comportementd,escoquesc\lind,riquessous chargesd,etype silo en E.L. 30

2.2. y'rpproche analytique du calcul des coques cylindriques

66silott
sous charges axisymétriques de type

L'analyse des coques cylindriques minces en élasticité linéaire soumises aux cas de
charges de type silo est complexe. Pour résoudre analytiquement le problème, plusieurs
hypothèses sont à formuler sur la théorie des coqueset sur les charges à considérer.

Les hypothèses ad.optéesdans ia formulation sont les hypothèsesde Kirchhoff-Love. Les

relations entre cléformations €;i de la coque et les efforts Ni par unité de longueur d'une
part et entre les courbures Ki et les moments Mi par unité de surface d'autre part sont
donnéespar les équations (1-17) chapitrel paragraph (I.2.2.1).

Les charges à considérer dans cette analyse sont axisymétriques et correspondent aux
pressionsde Ia matière ensiléesur les parois des silos à la fin du remplissage. Ces pressions
sont calculéesselon la norme NF P22-630 et l'Eurocode 1.

2.2.L, Formulation du Problème

:
Considérons figure (2.1), un silo métallique cylindrique à fond piat de diamètre "(D
2rrr, de hauteur "H" et d'épaisseur"t". On I'assimile à un cylindre à base circulaire;
soumis à des charges de tYPe silo.

Ces charges sont les contraintes de frottement f (") et les pressionshorizontales P"(z) de
"
la matière ensilée agissant sur ia paroi. Elies sont donnéesselonla norrne française P22-630
ou I'Eurocode 1 par les expressionssuivantes :

P.(r): C s {1 - e xpl -C r(H - t)l} ; TQ) : C, {t - expl-C"(H - ")l}

Ori C1, Cz et Ce sont des constantesdéfinies par :

Àp,. 'frn
.t
C1 :yh; Cz: t u3:-
Th I,r

Les valeurs de 7, À et p sont donnéespar la norme NF P22-630 [45]

Chapitre II ComTtortementd,escoqaescglindriquæ so'uschargesde type silo en E.L. 31

virtuellesda.nscecasde charges(chargesa>cisymétriques)
Le principe d.espuissa.nces s'ecrit:

(rHfH
lvu,, I fQ)u,(z)d,z-l N,(z)u,(z)d'z
toro to (2-1)
I r, rt 12 1H
I( vu,, I e.ç"1u,dz--l Noldz--, luI"u,,,"(z)dz
-'
Jo Jo
- r t z J Io

Les équations variationnelles (2-1) conduisent alors après intégration par partie, aux
équations d'équilibre suivantes :

ltr : f (z)
I "''
It12
-- P'(z)
-;", - -r2N1",," (2-2)
1
I IvI"(H): M,,"(H): 0

Mzy+AMzylùù

NI"+àNI"D"ù

Mp€Myy'ôdy

Fi,g(2.1) ModéIisation d,ela coquecylinilri'que

 II
Cha,p,itre Comportementiles coquesc\lind.riquessous clr'argesde type silo en E.L. 32

Par ailleurs les relations de comportement sont :

N":3(u,,"-,+)
Ns:3é+,u",") (2-3)

NI"--{-*r,,"" ,

Les composantes du champ de déplacementsU, et [J" ne dépendent que de la variable z

définissant la côte par rapport à la base du cylindre'
En introduisant les expressions(2-3) dans (2-2), on obtient :

=*lr,u).'+]
(# ++onu,) (2-4)

L,expression (2.4) représente l'équation différentielle gouvernante du comportemeni des

coques cylindriques sous charges de type silo.
Or) a et D sont des constantes définies par :

n:ffi flexionnelle);*^ -
(rigidité 3# (2-5)

2.2.2. Résolution analYtique

L'équation (2-4) peut se mettre sous la forme suivante:

(# : ee)
++onu,) (2-6)

où:
QQ):
ll'"u'."+l
La solution générale de I'expression (2-6), est de la forme suivante (voir annexe 1):

u,(,):"*p(#ll;.<ft)- B.'"(Ë)] - D"r"fffi))r:r:"rr,

*"r-A)
+exv(ffi11t
 Chapitre II Comporternent iles coquescylind,riquæ souschargesde type silo en E.L. 83

vo\ ( z
/ \: @
4aa

Les constantesA, B, C, et D peuvent alors être déterminéesà l'aide des conditionsaux

limites en z : 0 (encastrement)et z:H (bord libre). Da.nsnotre cas, on supposepour
simplifier les calculs que :

nuQ
ot

Ce qui signifie que est trés Bræd, alors exp tend vers I'infini. Les termes ayant
frH ftH
une ercponentiellepositive dans Ur(z) sont alors nuls de façon à satisfaire la condition du
bord libre.

A et B sont données par les expressionssuivantes:

- 1]
e: ffif,", * uc1)[exp(-c2H)
{ (2-8)
-A-ffit,c,*uc1)lrrrt"*oec,H)]
I B-

Ce qui conduit aux expressionsdu moment méridien et de la contrainte orthoradiale en

introduisant I'expression (2-7) dans (2-3) :

M'("):
'\ ) ,, !!.= ,, - B*Uær)*"(,)]e -s)
(e,s"fft)
(L-u2)(Draa) lo,,u.*n(#)
Avec :

s(z): -lnt, * uC1)exp[C2(z-

fr)]

ffy : -* (?) * u.N,(z) (2 - 10)

 Chapitre II Comportement d,escoquescglindriquæ sous ch'argesd,etype silo en E.L. 34

2.2.g'. Résultats de l'analyse analytique

cylin-
Les résultats analytiques exploitent une application numérique aux coques
NF P22-630
d.riques soumises à des charges de type silo, qui sont définis par la norme
et aux actions de
et l,Eurocod.e 1. Ces cha.rgescorrespond,entaux pressionshorizontales
pour un siio de quatorze
frottement. EIIes sont représentéesdans ies tableaux (2-1) et (2-2)
de quatre millimètres (la
mètres de hauteur, de huit mètres de diarnètre et une épaisseur
cellule "C")

Hauteurlm P,(z)-NF /(z)-NF P"(z)-EC.L f(z)-EC.1

0.000 22.L52 11.543 25.770 L4.812
L.250 2 0 .9 1 8 11.068 24.638 14.314
2 . 18 8 1 9 .8 9 2 10.365 23.663 i3.867
3.r25 18.773 9.7818 22.563 13.346
4.063 L7.549 9.7447 2r.323 L2.739
5 . 00 0 16.214 8.4484 19.925 12.031
6.250 14.239 7.4L96 17.780 10.900
7.187 L2.599 6.5651 15.932 9.8873
8 . 12 5 1 0 .8 0 7 5.6309 13.845 d. i u+D

9.063 8.8484 4.6106 11.492 7.3250

1 0 .0 0 6.7122 3.4975 8.8422 5.7L79
11.250 3 .5 5 3 0 1.8514 4.7740 3.1510
L 2 .5 0 0.0000 0.0000 0.0000 0.0000

m2l
tableau (P-1) pre.ssionhorizontale et actionsd'e frottementlkN I
de Remplissagecentré"
ile la matière ensilée(BIé) ilans la cellule"c' selon les normes "cas
Ilauteur[m P"(z)-NF /(z)-NF P"(z)-EC.t f (z)-EC.r
0.000 2 6 .1 0 9 L2.494 34.785 16.293
r.250 24.4L6 rr.779 33.258 75.745
2 . 18 8 23.377 11.187 31.941 r5.254
3.r25 22.033 10.543 30.457 14.681
4.063 20.568 9.8424 28.783 14.013
5.000 1 8 .9 7 6 9.0806 26.896 13.234
L6.632 7.9590 24.00r 11.991
6.250
7.187 1,4.694 7.0315 21.505 10.876
8.r25 12.583 6.0213 18.689 9.5749
9.063 1 0 .2 8 6 4.922L 15.513 8.0575
10.00 7.7894 3.7275 11.936 6.2897
11.250 4 .1 1 3 5 1.9684 6.4443 3.4661
12.50 0.0000 0.0000 0.0000 0.0000

tableau(2-2) pressionhorizontaleet actionsd'efrotternentlkNlrnzl

uid'angecentrée"
d,eIa matièreensilée(BIé) d,ansla cellule"c' selonles normes"cas de
Chapitre II Comportementdes coquæcylinilriquessous ch,argesde type silo en E.L. 35

Les figures (2.2), (2.3) et (2.4) représentent respectivement la variation des fonctions (
U,("), M,(r) et N, ) le long des génératricesdu silo.

'oI
12l
'l
'oI
'I
Èr
()
ff!
k
c)

.É
'I
4l
I

z^ll

,f
I
- -70L
t23456 24681012
Dépl acementradial (rnm) Déplacement radial (mm)

Fi,g(2.2) Déplacernentradial

--- Analpiçe (remplissage)

+ Ànrlrrfimre /widanoc cenltÉe
t4

t2
*cÏl ;elon l a l
l0

8 I
k
a
c) 6
C)
6
=
€
H

4 4

a
2

0 0

_rL- -2
-4
20 -10 0 30 ,10 50
10 20 -200 20 40 60 80
Moment méridien (N.m) Moment méridien (N.m)

Fig (2.3) Moment membro,naireméridien

Chapitre II d'etvpe silo en E.L.
Comporternentdes coquescylind,riquessous ch,arges 36

t4 l4

1t
t2

l0 10

eE 6t
k

.96 do
(€
H
4 4

t
2

0 0

-z -) t_
0102030Æ i0 010203040506070
Contrainte orthoradiale(MPa) '' Containte orthoradiale(MPa)

Fi g ( 2./) Contraint e rnembranaire orth oradiale

2.3. Approche numérique du calcul des coques cylindriques

(6silott
sous charges de tYPe

Dans cette section, et dans un premier temps, on utilise I'approche numérique pour
valider ]e modèle analytique proposé. Ensuite, dans le pragraphe (2.4.) en examine
numériquement I'effet des charges dissymétriques sur le comportement des parois cylin-
driques.

Cinq cas de charges non uniformes, correspondants au-xpressionshorizontales et aux ac-

tions de frottement dans la cellule "C" (annexe 3), seront analysés :

- les deux premiers cas de charges(*), représentent Ia vidange excentrée selon les mesures
expérimentales de la base de Chartres;
- le troisième-représentela vidange excentrée selon la norme NF P22-630;
- le quatrième représente e remplissage centré selon les mesuresexpérimentales de la base
expérimentale de Chartres.
- Ie dernier cas représente le remplissage centré selon I'Eurocode 7 1241,avec la prise en
comptedu "PATCH LOAD";

(*) L" deuxièmecas serareprésentédans I'annexe2-

Chapitre II Cornportementdes coquescylindriquessous chargesde type silo en E.L. 37

Tous nos calculs numériques sont effectués à I'aide du logiciel Abaqus (version 5.5 et
5.6) sous la licence académique de Hibbit, Karsson et Sorensen,Tns.,en considérant le
comportement élastique linéaire de la coque.

2.3.L. Cas de charges axisymétriques

Notre étude porte sur un cylindre à fond plat ( sans trémie ) de 14 métres de hau-
teur, de 8 mètres de diamètre. Le rapport hauteur/diamètre est d'envion I.7 assure un
comportement avec effet de silo. On considèregénéralementque cet effet peut être pris en
compte à partir d'un rapport égai à (0.8) (Norme P22-630). Le silo est construit en tôles
lisses d'une épaisseur de 3.5mm.
Les caractéristiquesm{snniques des tôles sont le module d'Young E :2.te5MPa et le
ccefficientde Poisson u :0.3.
Compte tenu de la symétrie des actions et de la géometrie du silo considéré, nous avons
considéré que le quart de Ia structure. Le maillage est effectué aux moyens d'éléments
coques à huit nceuds(*), voir figure ci-dessous.

Mai,llage ilu quart d,ela coque cglinilrique (éléments coquesS8R5)

(*) Les éléments coques à huit næud (SSR5) donnent un meilleur résultat par rapport aux
éléments à quatres nceuds (S4R5).
Ch,apitreII Cornporternentdes coquescylinilriquæ sous chargesd'etype silo en E.L. 38

2.3.2. Comparaison des résultats analytique et numérique

Le programme Abaqus noun a fournis les déformations, les déplacements et les rota-
tions sur tous les nceuds, les contraintes et les moments de flexion au milieu des éléments.
Les principaux résultats comparés aux résultats analytiques, sont représentés sur les figures
(2.5) et (2.6).

(i\ Cas d,eRemplissage

v
h
lr
a)
J
ct

t2345
Déplacementradial (mm) Déplacementradial (mm)

Fig (2.5a) Euolutiondu déplacementrnd,ialle long d,esgénératrices

E tr
H

g
É

-10010203040 -20 0 20 4a 60
Moment méridien (N.m) Moment méridien (N.m)

Fig (2.5b) Euolution ilu ntoment membranaireméridi,enIe long des génératrices

Chapitre II Comportementd,escoquescglinilriquessous ch,argesile type silo en E.L. 39

l4

12

l0

É6

h
fi

c) 96
<É

0 5 10152025 30 01020304050
Contrainte orthoradiale (IMPa) ConEâinte orthoradiale(MPa)

Fig (2.5c) Eaolution de la contrainte mernbranaire orth,oradialeIe long des génératrices

(ii\Cas d,eaidanoe centrée

|<

g I
=
G' Gt

t234567-024681012
Déplacement radial (rnm) Déplacementradial (mm)

rad,ialle long des génératrices

Fig (2.6a) Eaolution iJu d,éplacement
Ch,apitreII Comporternentd,escoquescylind,riquessous ch'argesd'etype silo en E.L. 40

-r- Abaqus
-+-Analytique
\f - -
r.j r selonilEcl

8
L
(l)
o
d

-'40
-20 -10 0 l0 20 30 ,10 -20 o 20 ,{o 60 80
Moment méridien (N.m) Moment méridien (N.m)

rnéri,dienIe long des génératri'ces

Fi,g(2.6b) Eaotution d,urrLomentn'ternbranaire

t4 I4

t2 t2

l0 l0

8 E

ti F
(l) c)
6

4 4

0010203040 0010203040506070
Contr'âinteorthoradiale(MPa) Contrainte orthoradiale (MPa)

Fig (2.6c) Eaolution de Ia contraintemernbranaireorthorad'ialele long des gé'nératrices

Ch'apitreII Cornportementd,escoquescglindriquessous ch,argesile type silo en E.L. 4l

2.3.3. Interprétation des résultats analytique et numérique

De cette étude comparative on remarque que :

- la coque se comporte essentiellement coûrme une membrane dans notre cas à partir
des figures (2.3), (2.5b) et (2.6b);
- il nty a donc qutune flexion concentréeau voisinagede I'encastrementl
- la contrainte membranaire orthoradiale est maximale,à une cote légèrement supérieure
à celle de l'encastrement. Elle permet le dimensionnement de l'épaisseur de la tôle à Ia
traction, figures (2.4), (2.5c) et (2.6c);
- l'amplitude du déplacement radial de I'ordre de six à onze millimètres est situé
égalementdans Ia partie bassedu silo, figures (2.2), (2.5a) et (2.6a).

Les figures (2.5) et (2.6), montrent une bonne concordanceentre les résultats analytiques
et numériques en élasticité linéaire dans le cas de chargesaxisymétriques. L'erreur entre
les résultats analytiques et numériques est inférieure à 10%. Cette observation montre la
validité de l'approche analytique et prouve que, les chargesont bien été décrites dans notre
étude numérique. Par ailleurs les hypothèsesadoptéesdans la partie théorique sont bien
justifiées.

Nous pouvons donc poursuivre notre analysenumérique en considérant des cas de charges
quelconques, à savoir des pressions horizontales non-uniformes décrivant la vidange ex-
centrée.

2.4. Analyse numériques du comportement des coques cylindriques

sous charges dissymétriques

La distribution des pressions dans les silos peut varier considérablement dans les
deux plans vertical et horizontal. Les réglementsactuels tels que la NF P22-630 ou l'EC1
donnent peu d'informations sur la manière de considérerla variation des pressions dans le
dimensionnement. Après illustration des types de pressionsnon-uniformes qu'on rencontre
dans les silos, on vérifiera si les perturbations liées à Ia non-uniformités des pressions
horizontales et des actions de frottement sont suffisamment importantes pour la prise en
compte dans le dimensionnement.
Une attention particulière doit être portée sur le fait que la distribution des pressionssur les
parois d'un silo cylindrique est rarement uniforme autour du périmètre pour un niveau de
charges donné. Il y a plusieurs causesqui provoquent cette non-uniformité parmi lesquelles
on peut citer :
Chapitre II Comportement des coquescylindriquæ sous chargesde type silo en E.L. 42

) Ia ai,dangeeacentréei
- I'ensoleillement;
- Ia aiilange séquentielle;
- Ia ui,ilangeet Ie remplissagecontinu;
- Ie aentl
- le séisrne;
- Ie remPli'ssageexcentré.

2.4.L. ANALYSE 1 (vidange excentréeselon les résultats de ia base de Chartres) [3a]

de
Notre première analyse porte sur un cylindre à fond plat de 14 mètres de hauteur,
g mètres de dia,mètre et d'épaisseur variable de 2 à 3.5 mitlimètres. Le silo est considéré
encastré à la base et libre à son bord supérieur'
considérons
La structure a perdue sa symétrie suite à la non-uniformité des pressions'nous
donc la structure entière dans ce cas voir figure ci-dessous.
Autrement dit :
- Le long de Ia circonférence' on dispose32 éléments;
- Le long des génératrices et sur une hauteur de 2.5m, nous avons adopté 10 éléments soit
(soit 0'714m
une hauteur d.e0.2bm par élément, et 14 éléments sur le reste de la hauteur
environ par élément).
: 2'1e5MPa et Le
Les caractéristiques nlsnniques d.ela tôle sont ie module d'Young fi
cceffi.cientde Poisson u :0.3.

Maillage ile Io' coquecyli'ndrique entière(éléménts S8R5)

Chapitre II Cornportement d'escoquescylindriques sous chargesde type silo en E.L. 48

Les charges adoptées dans cette analyse sont les résultats expérimentaux de la base de
Chartres. L'ouverture de la zone de répartition des chargesest supposéede 45 degrés. Ces
résultats sont représentéssur la figure (2.8).

I*l
l------Prurænr I

I Gésémtrice 2

. Matière @ilée'BLE-
Ft _l

/a
s
É6 Qe

l,æ*l
1..-.--.nuæ"r I
-
. Géaéatica I

I lr{atiùe æiléa'BLE-

20l'ogt60 Ë30.}t10

Pr6iGtkNb'zl Prryim [ktl/mz]

-h@io! hdiatslo

.------Fdlmd

'Gédræie 3
' MstièF @ilé.

o
â6
É6

û?s30
PrsioBr.ûhl Præim [kN/mz]

;FQ^*"*

Fig (2.8) pression horizontale et action de frottement selon

Ia basede Chartres pour les quatresgénératrices
 clt'argesile type silo en E'L'
ch,api,treII Comportement iles coquescylind'riquet::ut

2.4.L.L Résultats de ltanalYse 1

?s0

200

2 uio
J

E roo
E
sl 50

-50
5t0
Circonférence [ml

Fig (2.9a) Variation d,el'effort tad'i'al

2
t

a
t,
g

.9
o
.9

Circonférence[n]

Fig (2.9b) Variati'ond,el'effort circonférenciel

300

z@

z
J l@

0
€

-100

-200

Circonférence[n]

Fig (2.9c) Variation de l'effort atial

Chapitre II Cornportementiles coquescylinilriquæ sous chargesde type silo en E.L. 45

o
9^
qtr
42

o
o

Circonférence[m]

Fig (2.10a) Variation d,umoment nembranaire méridien

100

O^ 50
a4
d9 0
EE
:p
qE
FF -50
È'Ë

-1m

Circonférence[m]

Fig (2.10b) Variation ilu moment membranairecirconférenciel

80

O^
4
.5É
€-
F4 20
qo
OE
EO
:l:l
0
6t
E'6
Èe -m

{)
5r0r5 z0
Circonférence[m]

Fis (2.10c)Variation d,umoment membranaire de ci,saillernent

 en 46
Ch.apitreII ComPortement des coquescylind,riquessous ch'argesd'etype si,Io E.L.

o
.b
€^
EÆ
5>
É v
Fto
ÉÉ
o.:

FA

Circonférence[m]
e
Fig (2.1la) Variation ile Ia contraint e mernbranairernéri,d,ienn

q^
.:s
E>
È v

6.9
o!)
É )!)

EE

Circonférence [n]

Fig (2.1lb) variation ile la contrai,nternembranaireci'rconférencielle

EO

o
.bî q
E>
2v

oo zo
o9
.:t
5q
Eo
CE

'at
l0 15
Circonférence [m]

Fig (2.11c) Variatzon d,eIa contraintemernbranaired'ecisaillement

Chapitre II Comportement des coquescylind,riquæ sous ch'argesd'etype si,lo en E.L- 47

0.æoE

0.u16
o

0.o04

o 0.æ@,

{.qxz
o

4.0004

{.0006
o 5
chconférence [m]

Fig (2.12a) variati,on de la d,éfonnation nrembranaireméridienne

o
o
I
o
È
p

o
o

{.0005

c:

Ë
{.ær
o 5 tu
ch"oofé'"*u 1-1

Fig (Z.1Lb) Variation ile la déformation membranairecirconférencielle

o 5 'u
cir"ooféo*" 1t1

Fig (2.12c) Variation ile la d,éforrnationmembranaired'eci,saillernent

Cltapitre II Cornportementdes coquescylindriquessous chargesd,etype silo en E.L.

--e-niva0375m ---r-nivau32l4n -o-dvau6.07lm ---+-niva735?m ---+-nivau ll.43m

100

o
a
i
o
ë i
iÂ
Ê a
\
/, l=ù
F
E=l tA
eÉ 0
w
V
V V
:it E-
o
q
\lt

p
-50

V/ V
-l0o I
0 5 10 É
cit"oofé*o"" 1-1

Fig (2.13a) Variati,onilu déplacernentcirconférenciel

+uivau03?5m *aieu32l4m +nivau6.0?lm +nivau7857m --r-dvaull.43o

1æ

50

14
\-n A. ,t,\
v\7
F
A. Jà
o

o
0 Y
\ Y& v
o

-50

-100
5I015?9
Circonférence [m]

Fig (2.19b) Variation du iléplacementradial

o
o
a

5l015j
Circonférence [m]

Fig (2.13c) Variation ilu déplacementaûal

Chapitre II Comportementiles coquescylind,riquessous chargesde type silo en E.L. 49

d)

Ezo
a
6o
2 -zo

50

Circonférence [m]

Fig (2.14a) comparaisondes moments au ni,ueaude 1.25m

2æ

160

80
z
0

€0

10 15
Circonférence [m]

Fig (2.Ub) comparaisondes nronrentsau niaeau d,e9.211m

300

2æ
qlæ
'z

Ëo
É -t*
-2û

-3æ

0510152025
Circonférence [m]

Fig (2.14c) comparaison iles rnoments o,uniùel,u de 7.857m

 d,etype silo en E.L.
chapitre II Comportement iles coquescylind,riquessous charges

80

60

,10
€
g20
cra
Ë'
d

E"o
4

€0 25
0
Circonférence
[m]

Fig (2.15a)comparaisond'escontraintesau ni,aeaude 1.25m

150

100

50
Ë
o
0

-50

-100
25
0
Circonférence [m]

Fi,s(2.15b)comparaison d'escontraintes au ni'aeaud,e9.21/1m

2æ

150

^ l0o

a
.Fv

R -50

-læ

-150

Circonférence [m]

au niaeau de 7'857m
d,escontra,intes
Fig (2.15c) comparo,ison
Ch'apitreII Comportementd,escoquescylind,riquæsous chargesd,etype silo en E.L. b1

offi

û(pl

5 oo@
d

EO

{.@2

.oml

.omd

Circonférence [m]

Fig (2.16a) comparaisoniles d,éformationsau niueau de 1.25m

om

o0@4

ç ow,

x^ 4cf,t2
.oCln4

.{rmt

Circonf&ence [m]

Fig (2.16b) comparaisoniles déformationsau niaeau d,e3.21/m

o(xB

o06

! oooor
Eô
_o

-a.m

.{lm4

{.@6

Circonférence
[m]

Fig (2.16c) comparaison iles déforrnations au niaeau d,e7.857rn

 chargesd'etYPe si,loen E.L. 52
Ch.apitreII CornPortementdes coquescglinilriques sous

L2

t0

k
c) 6

0 -250 -2@ -150 -100 -50

Moment membranaire méridien (N'm)

méridien le long d'esgénératrices

Fig (2.17a) Euolution d,umoment rnembranaire

t2

l0

H
() 6
)

0
-50 0
Moment membranairecirconférenciel (N'm)

Fis (2.17b)Eaolutioni l u m o m e n t m e r n b r a n a i ' r e c i r c o n f é r e n c i e l l e l o n g d e s g é n é r a t r i c e s

 Ch,apitreII Cornportementdes coquescglindriques sous chargesd,etype silo en E.L. 53

-s-G. | +G.2 -.-G.3 -+G.4

L2

l0

|r
c) 6

0
0 50 100
Contrainte membranaire méridienne (MPa)

Fig (2.15a) Eaolution de Ia contrainte mernbranaireméridi,ennele long des génératrices

t2

t0

f{

o 6

Fig (2.18b) Euolution d,ela contrainte membranaire circonférencielleIe long des génératrices
Ch,apitreII Cornporternent iles coquescylind,riquæ sousch'arges
de tgpe silo en E-L- 54

lz

l0

^ E

r{
o 6
6l

0
04Bt'z
Déplacementcirconférenciel (mm)

Ie long d,esgénératrices
Fig (2.19a) Euolution d,ud,éplacementcirconférenciel

t2

10

k
a
o 6
d
-
H

génératrices
Fig (2.1gbi) Eaolutiond,uiléplacenent rad,iallelong iles
 Chapitre II Com,portementd,escoquescylindriquæ solrschargesd,etype si,Ioen E,L.

12

l0

c) 6
=

0
20 40
Moment (N.m)

Fi,g (2.20a) Cornytaraisond,esmoments le long de Ia génératri,ce1

-Momentcirconférenciel -......Ms6sntméridien

t2

l0

H
a
() 6
d

I
i
0
-100

Moment [N.m]

Fig (2.20b) Comparaison iles rnom,entsle long d'eIa génératrice I

I
Chapitre II Comportement d,escoquescylindriquet to* chargesde tgpe silo en E'L' 56

t2

10

^ 8

g 6
.€
4

0
0 100
Contrainte QIPa)

1
Fig (2.21a) ComParaisondes contraintesIe long de Ia génératrice

12

10

^ 8

.9 6
d

0
Conrrainte(MPa)

d'escontraintesIe long de la génératrice3

Fig (2.21b) Cornpara'i,son
Chapitre II Comportement des coquescylind,riquæsous chargesde type silo en E.L. 57

de Ia coquecylind,rique(niaeau 0.375m)
Configuration d,éform,ée

Configuration déformée ile la coquecylind,rique(niaeau 1.25rn)

Configurati,on déformée d,eIa coquecylindrique (niueau 9.211m)

(niueau5.178m)
ile la coquecylind,rique
Configurationiléformé,e
Chapitre II Comportement d,escoquescylindriquessous chargesiJe type silo en E.L. 58

Configuration déform,éede la coquecylindrique (niueau 6.071m')

configuration déformée de Ia coquecylindrique (niueau 7.857rn)

Configuration défonnée ile Ia coquecglindrique (niuea'u 9.375m)

configuration déforméed,eIa coquecylindrique(niueau11.13m)

2.4.l.2Interprétation des résultats de l'analyse 1

Nous constatonsà partir descourbesdesfigures(2.13b)et (2.19b)que I'amplitude

du d.éplacementradial est plus importante dans la pariie supérieuredu cylindre. Son
Chapitre II Cornportementd,escoquescElind,riquusous chargesd,etype silo en E.L. bg

*plitrrd"maximale est de l'ord.re de g cm qui représentepresqueIe centième de la hauteur.

Ce déplacement est responsabledu voiiement des parois.

La même remarque est observée sur les moments de flexions ( moment membranaire
méridien ), d'après les figures (2.10a), (2.I4), (2.I7a) et (2.20). Le maximum est de I'ordre
de 350 N.m se situe au-delà du milieu du cylindre.

I1 en résulte des représentations précédentes,que la flexion n'est plus localisée seulement

en partie basse comme I'indique le résultat de I'analyse avec des charges axisymétriques
du paragraphe (2.3.).

Le calcul du rapport des sollicitations membranaires méridiennes (de flexion) / sollicita-

tions membranaires circonférencielies à différent niveau du cylindre suivant les génératrices
est résumé dans le tableau (2.3)

Nous constatons que les efforts membranaires circonférencielssont plus élevés par rapport
aux efforts membranaires de flexion sur toute la hauteur du silo. Le rapport de ces efforts
est de 0.821 au niveau 3.274m;

Nous constatons également que les moments et ies contraintes membranaires cir-
conférenciellessont prédominantes jusqu'à la côte 3.274 mètres, mais, au-delà les moments
et les contraintes méridiennes (de flexion) induites par les chargessont plus importantsl

Les sollicitations membranaires méridiennes (moments et contraintes) sont pratiquement

trois fois plus grand que les sollicitations membranairescirconférencieilesà la côte de 3.214
mètres.

I1 est donc essentiel de tenir compte de cette flexion, par la prise en compte de la non-
uniformité des pressions horizontales dans ie dimensionnement des silos.

Niweau [m] 0.375 3.2r4 6.071 7.857 11.43

Effort méridien ftN] 35 115 s7 L07 220
Effort circonférenciel ftN] 750 140 s70 400 390
Rapport des effort 0.0466 0.821 o.262 o.2675 0.564
Mornent méridien [N.:n] 13 250 s80 385 L70
Moment circonférenciel [N.m] t4 75 L20 135 60
Rapport des moments 0.928 3.33 3.166 2.85 2.833
Contrainte méridienne [MPr] 13 150 2LO 160 r47
Contrainte circonférencielle [MPa] 125 45 95 95 80
Rapport des contraintes 0.104 3.33 2.21 1.684 1.8375

tableau (2-9) Rapport des sollicitations, Mérid,iennes/ Circonférencielles

le long iles génératrices
 tgpe silo en E.L. 60
Chapitre II Comportement d,escoquescylind'riquessous clrargesd,e

2.4.2. Analyse 2 "vid.angeexcentréeselon la norme NF P22-630" [45]

actions calculées
Ceite analyse porte sur le même silo que l'analyse 1, on considèreles
selon la norme française P22-630 en phase de vidange excentrée.
Formulation des clLarqes
plus ou moins
Lors d.e la vida,nge excentrée, une pression horizontale non-uniforme,
cette pression est donnée
irrégulière, s'exerce sur Ie pourtour de la paroi de la cellule,
par I'expression suivante :
Pn.r,"(t) : C (2- 1i)
"Po,r(z)
de la cellule' 11est
C, dépend de l'excentricité de I'orifice de vidange ainsi que du rayon
donné par la relation suivante

+ 0.a( ;1)
C e : r nax( O.g (2 - L2)

où:

(- 9
r
e : est ltexcentrement de I'orifice de vidange

r : est le rayon de la cellule

Pn,o("): est Ia pressionhorizontale à la vidange'

des actions de frot-
Les autres efforts (pression verticale, action de frottement, et sollune
(voir chapitre I paragraph
tement) restent inchangés par rapport à la vidange centrée
1.1.3.1).
la vidange
La répartition de la pression horizontale selon'la norme française conceroant
excentréeest représentéesur la figure (''")'rr*
0sPb,
(qum)

''xïd;'

Fis (2.22) Répartition en uidangeexcentrée

selonla NF P22'630
 Chapitre II Comporternent des coquescylindriquæ souschargesile type silo en E.L.
61

2.4.2.L Resultats de l,analyse 2

g
to
e
o
e 0 ol
o
=
o
.t0 I
o -æ " oI l
.3{t .30
|
4
*[
to l5 t5
. Cimnfémcc [m] Circonférence[m]

Fig (2.23) Variati,on d,esiléplacernentsle long d,ela circonférence

É
ts
2
â2 €E
69
od
tro tr9
:o.9: : )!l

:9
È';

t0 lJ l0 t5
Circonférence [m]
Circonféræce [m]

Fig (2.2/) Variation d,esmoments le long ile la circonférence

P^
€Ê
Ëë
4s
EE
!E
.58
ËE
c:'6

l0 15
Circonfércocc[m] Circonfémæ [rn]

Fig (2.25) Variation d,escontraintesIe long ile la circonférence

Chapitre II Cornporternent iles coquescglind,riquessous chargesd,etype silo en E.L.

l2
l2 I lc qutrt
t les trois quads
r (qtrycrûre de l0o)

l0

8
E

d6
dc

, a

510152025æ 481216202424

Déplacementradial [mn] Déplacementradial [mm]

Fig (2.26) Eaolution ilu déplacementrad,ialIe long d,esgénératrices

1a
1' t le qurt * 16 trois qutrts

. (owstne de 10P)

10
l0

E
8

66
d6

!
4

0
'10 -20 0 20 40
40 -2o 0 æ 60
Moment membrranaire Moment membranaire
méridien [N.m] méridien[N.ml

génératrices
Fig (2.27) Eaolution iLu moment rnembranaireméridi'enle long des
Chapitre II Comportement des coquescylintlri,quessous chargesde tApe silo en E.L. OB

'le qust

r (owerare de l0o) I l€s kois quans

l0

I
8 t

de g6
F

0t-
-lJo -t@ -50 0 50 læ 150 26 -120 {0 4 0 iO 80 l2O 160 200
Moment membranaire Moment membranaire
circonférenciel [N.m] circonférenciel [N.m]

Fi'g (2.28) Eaolution du moment rnembranaire circonférencielIe long d,esgé,né,ratrices

* le quart * les t$is quarts

* (ouverture de 10")

8 E

=
H
g8
H
A6
=
€

0E
{O€0-{O-æ0æ40 -80 {o 40 -20 0 20
Contrainte membranaire Contrainte membranaire
méridie,nne[MPa] méridienne [MPa]

Fig (2.29) Euolution ile Ia contrai,nterLembranaire mérid,i,ennele long iles génératrices

 de type silo en E.L.
Chapitre II ComPortementdes coquescglindriquæ sous ch,arges

* le çart * lestrois quarts

l2
r,(ouverturede 10o)

10

t
8

o
66

0 -2002040@801@
'-20 4 60 t0 100 r20
o 20
Contrainte membranaire Conirainte membranaire
circonférencielle [MPa] circonférencielle [MPa]
le long d,esgénératrices
Fig (2.30) Euolution d,eIa contrainte membranai'recirconférencielle

4{'

æ
20

l0

I 0
o
o
.10

.20

.æ
4{'

Chconfércnce [m]
niaeau 0'375rn
Fig (2.31) Compataison ilesiléplacements

25î)

200

I tso

5 too

50

0
510r52025
Chconférence [m]

Fis (2.32) Comparaisond,esrnornentsni'aeau0'375m

 Chapitre II Comportement d,escoquescylind,riquæsous ch,arges
de type silo en E.L. 6b

Chconférence [m]
Fig (2.33) Comparaisondes contraintesniaeau 0.375m

2.4.2.2Interprétation des résultats de lranalyse 2

Nous constatons à partir des courbesdes figures (2.23) (2.26) et (2.31) que I'amplitude
du déplacement rad.iai est plus importante dans la partie inférieure du cylindre. Son
amplitude maximale est de I'ordre de 3.8 cm. Ce déplacementest responsabledu voilement
des parois.

La même remarque est observéesur les moments de flexions ( moment méridien ), d'après
les figures (2.24) (2.27) et (2.32). Le ma-ximum est de l'ordre de 75 N.m se situe au niveau
de la base du cylindre.

On remarque également que le rapport sollicitations membranaires,circonférencielles/ sol-

licitations méridiennes au niveau de la base est de deux à quatre fois plus grand, ce qui
prouve que la flexion n'est pas importante, et qu'il faut tenir compte seulement de l'effet
membra.naire de la coque cylindrique pour Ie dimension.rement.

11en résulte des représentations précédentes,que la flexion est toujours localisée en partie
basse courme I'indique le résultat de I'analyse avec des charges axisymétriques paragraphe
(2.3.1). Sauf que les amplitudes des sollicitations ainsi que des déplacements sont plus
élevées.

En comparant les résultats (sollicitations) de I'analyse I et I'analyse II, on trouve que

les charges proposées par la NF P22-630 ignorent des su4>ressions, qui se developpent
au delà de la mi-hauteur de la structure. Ces surpressions donnent naissance à des con-
traintes membranaires de flexion plus importantes que les contraintes membranaires cir-
conférencielles.
 ch'argesde type silo en E'L' 66
chapi,tre II Comportement d,escoquescglind,riquessous

chartres" [34]
2.4.3. Analyse 3 "remplissage centré selon les résultats de la base de

les charges représentant Ie

Dans cette section, on applique sur la même structure
la base de chartres (figure 2'34)'
remplissage centré selon les mesures expérimentales de

,,C,,
cellule
Fig (2.3|) Pressionsen phased,e remplissagecentré dans la
selon Ia basede Chartres

2.4.3.L Résultats de ltanalYse 3

96

'I
I
I

o.mr o.dE 0.@3 0.m4 o.ffi "I

0 o, c4 0.6 0.8 I

Déplacement radial [mm]

D€placene,nt circonféreaciel [mm]

Ie long des génératrices

Fig (2.95) Eaotution d'esd,éplacements
 Chapitre II Comporternentiles coquescylindriquæ sous chargesd,etype silo en E.L. 67

9c

-10 0 10 20
Moment membranaire
méridien[N.m]

Fig (2.36) Eaolution des mornentsle long d,esgénératrices

k
o
a6

-10 0 lO ZO æ ,10 50 4 .30 -20 -10 0 l0 20

Contrainte membnnaire Contrainte membranaire

circonférencielle [MPa] méridierrneMPal

Fig (2.37) Eaolution iles contraintesle long desgénératrices

I
Chapitre JI Comportement ilescoquescylindriquessous ch'argesile tupe silo en E.L. 68

2.4.4. Analyse 4 "remplissage excentré selon l'Eurocode I" l24l

Dans ce paxagraphe on examine I'effet des chargesadoptées par I'Eurocodel sur Ie

comportement d.e la structure "silo" avec Ia prise en compte de la surpression appellée
"Patch Load".
Forrnulation des ch'arges
par les équations
La pression horizontale ainsi que les actions de frottemeni sont données
(i-S) et (i-9), chapitre I.
Le patch load P, est donné Par:

P p :0 '2 0 P n t; 0:114e;fd. ( 2 - 13 )

Ori :
e; et d," représentent respectivement I'excentricité de l'orifi'ce et le diamètre de la cellule.
Cette surpression Po se développe.âone certaine hauteur ,S(voir figure ci-dessous)donnée
Dal :

S :0.2d'"

* Matière ensilée (BLÉ)

' P-o(z) (nnPlissago uæntr6)

avec la prise eB comPte du
'Pâtsh I-oad"

k
a
=o

0r
0 510ut2025
PressionstkN/m'?l
Side elevation

nC" selon l'EC1'

Pressionsen phased,e remplissageeûcentré dansla cellule
 Chapitre II Comporternent d,escoquescylinilriques souschargesile type silo en E.L. 69

2.4.4.LRésultats de ltanalyse 4

Ê
o 6)

€6

, rûf 1 tt? rJ tù7 z tù7 zJ to.?3 10?3J ltr4 to.t 0,t 0.2 0.3 0.4 05 0.6 0.7

déplacementcirconférenciel [m] déplacementradial [mm]

Fig (2.38) Eaolution d,esdéplacementsIe long des génératrices

4-404812 -æ0æ,10
moment membranaire moment membranaire
circonférenciel [N.m] méridien [N.m]

Fig (2.99) Eaolution iles moments Ie long ilesgénératrices

I
Chapitre II Cornportementdes coquescylind,riquettry chargesde type silo en E'L' 70

,6
rl
d

.24 .16 -t 0 E 16
r0 20 30
contsainte membranaire contrainte membranaire
circonf&encielle [MPa] méridienne [MPa]

Fig (2./t0) Euolution des contraintesle long d,esgénératrices

2.4,4.2Interprétation des résultats de ltanalyse 3 et 4

radial est plus

Les figures (2.35) et (2.38) montrent que I'amplitude du déplacement
importante dans Ia pariie inférieure du cylindre'
méridien ), d'après
La même remarque est observée sur les moments d.eflexions ( moment
du cylindre'
Ies figures (2.36) et (2.39). Le maximum est au niveau de ia base
membranaires sont
on remarque également (figures (2.37) et (2.40)) que les contraintes
par sa forme'
plus importantes que les contraintes d.eflexion. La coque résiste uniquement
La flexion est toujours localisée en partie basse'

cylindri'quesousch.argesaxissymétriques
Fi,g(2.11) configurationiléforméeil,eIa coque
ou remplissagecentréselonl'Eurocod,e
(ùi,d,ange l paragraphe[2.3'1' (i) et (ii)])
Chapi'treII Comportemen't,iles coquescylinilriquæ sousch,argesile type silo en E.L. 7!

Fig (2./2) Configuration iléforrnée d,eIa coquecylind,riquesous chargesnon-uniforrnes

(aidange excentreéselon les résultats d,ela basede Chartresparagraphet2.4.Il)

Fig (2.49) Conf'guration déformée d,eIa coquecylind,rique sous ch,argesnon-uniformes

(aid,angeeacentréeselon ta NF P22-690 paragraphe[2.4.2])
 tgpe silo en E'L' 72
Chapitre II Cornportementd,escoquescylinilriquessous chargesile

chargesatisymétriques
Fig (2.40 configuration iLéforméeile la coquecylind'riquesous
(remplissage centré selon les résultats ile Ia based,echartres paragraphe[2'4'9])

ch'argesnon-u'niformes
Fig (2.15) configuration d,éforrnéede la coquecylind'riquesous
d'u load"'paragraph'"[2'/''4'])
(rentplissage eacentré selon I',Eurocod,el auecIa prise en contpte "patch
Chapitre II Comportementd,escoquescylindriquæ sous clrargesde type silo en E.L. Tg

2.5 Conclusion

Au vu des résultats des a,nalysesnumériques présentés dans les paragraphes (2.3)

ef (2.4), nous pouvons conclure que la non-uniformité des pressions horizontales est un
phénomène préjudiciable pour la stabilité des silos cylindriques. La conséquencesur Ia
distribution des sollicitations dans les parois est nette. L'évolution du moment de flexion
suivant la génératrice est totalement quelconque, et le maximum se situe au-delà de la
mi-hauteur du silo. Contrairement à l'étude faite précédemmentdans ie cas de pressions
horizontales uniformes autour de I'a>revertical, paragraphe (2.3.1.), la flexion n'est plus
localisée à la base encastréeet devient prépondérante.

En observa,ntaussi l'évolution du déplacementradial, on constate une répartition qui n'est

plus semblable à celle que nous avons observé dans le cas des pressions uniformes autour
de l'axe vertical, I'amplitude maximale se situe en partie supérieure du silo.

Les résultats numériques obtenus montrent que cette non-uniformité entraîne des con-
traintes membranaires de flexion beaucoup plus importantes que les contraintes rn€rn:
branaires circonférenciellesdans certainessectionsdu cylindre. Dans ce cas, le moment de
flexion est un facteur important da.nsle dimensionnement des silos dans le cas des charges
dissymétriques (non-uniformes) le long du périmètre.
 Chapitre III

VOILEMEI\TTDES COQUESCYLINDRIQUES
SOUSCHARGESDE TYPE SILO ET
ETABLISEVIEI\TD'U}[E FORVIULEA}TALYTIQUE

3.1. Introduction

La coque cylindrique sans ou avec raidisseurs est souvent choisie dans la conception
des structures de nombreux composants industriels: des éIéments de construction nucléaire,
des ouwages de génie civil, tels que les silos, etc...

Le dimensionnement de cette coque est basé sur sa résistanceau phénomène de voilement,

qui peut causer la ruine de celle-ci bien avant I'atteinte de la limite élastique du matériau
la constituant.
Les solutions proposées dans la littérature pour la stabilité des coques concernant des
conditions simples de ùargement (compression axiale constante par exempie), sont bien
éloignées du cas de silo. Notre étude consiste donc a a,r:.éliorerl'état actuel des connais-
sances du voilement des silos cylindriques imparfaits, raidis ou non longitudinalement, à
fonds plats. Nous proposons une approche analytique pour le calcul des contraintes ultimes
de résistance de la structure du silo. Cette approche est baséesur une étude paramétrique
par voie numérique (Abaqus).
Chapitre III uoilernent iles coquæ cylinilriques sous chargesile type si'Io' ,D

ce
Comme nous ne disposons pa"sde résultats dans la littérature qui tiennent compte de
type de charges (charges de type silo), nous avons choisi de faire d'abord un calcul simple
par
sur une coque soumise à une compression axiale uniforme, deja traitée analytiquement
Timoshenko et Gere.

3.2. Flarnbage des coques cylindriques sous compression axiale

Dans cette étude, on se refère à I'exemple déjà traité par Timoshenko et Gere [63].
soumise
I1 s'agit d'une coque cylindrique parfaite, longue, mince et simplement appuyée,
pour que le
à une compression axiaJe uniforme (fis*" 3.1). Elle est suffisamment mince
est
voilement puisse se produire bien avant la limite élastique. Le matériau de Ia coque
homogène et isotroPe.

Cette étude présente un double objectif permettant :

la stabilité
1. de montrer la technique utilisée par Abaqus pour résoudreles problèmes de
étapes :
des structures métailiques. En général l'étude de la stabilité se fait en deux
- analyse aux valeurs propres pour estimer les chargescritiques et les modes critiques
associés;
- analyse non-linéaire en utilisant la méihode incrémentale de Riks ou celle de
imperfections
Newton-Raphson qui permettent d'étudier Ia sensibilité de Ia structure aux
géométriques.

2. de valider Ies résuitats analytiques en les comparant avec des résultats numériques'

la meilleure stratégie
Cet exemple (exemple de Timoshenko et Gere) est étudié pour définir
pour aborder Ie flambage des coques cylindriques sous charges de type silo'

en cornpressionaxiale
Fi,g(3.1) Coque cyli,nd,ique
Chapitre III aoilement d,escoqaescylinilriquæ sous chargesde type silo. 76

3.2.L. Analyse analytique de la stabilité des coques cylindriques

selon Timoshenko et Gere

En utilisant les mêmes notations de la figure (2.1) chapitre II. Les équations
d'équilibre de fla.rnfage sous compression a><ialeconstante sont données par les expres-
sions suivantes :

'***:o
îNy r -_! \0z'-! - + A I / I , ' -_l 0 - ,_ n
$ fj ;u: 1
-+,, ? N " n -* T
,w (3-1)
A0 A" A,
-r*i \ " ô , . , n , a z N I ,+ 0 2 M s , * ô , M n _ 0 2 M " s _ : 0r
U7*iYa*-6"2 ô200' æ- affi
Ces équations peuvent être écrites sous formes de déplacementsu, v et w tout en utilisant
Ies notations suivantes :
*
O':
T'F
N,(L - u2)
ô-
ET

t+u
At:
2

L-u
Az:

#.,o,#-,,y+',ffi:o
+"+,*+'# -'# . "l# - # * 2rAz#l-'r# : o
o,fur*
-ra,ô 2 u s*' ô u + G
0 u - w - | ô3u , ô3u
-,)æ "ô4w * tw1
+r"
aF a" ; " læ*r(2 al æl
-zo, o!^'
-'*'
ar'af _,
io
Soit :
TNTT
^:T
-
 d'etype silo'
Chapitre III uoilement d,escoquescylind,riquæsous chotrges 77

La solution générale du système d'équation (3-2) est de la forme suivante :

( " : A s i n ( n a ) c o sÀ
, - Bcos(na)sinÀ
{
[ , - Csin(na)sinÀ

:
En substituant cette solution dans (3-2) on obtient le systèmesuivant

e(À2+Azn2)*BnAù'*Cu\:o
t . ô _ :0
\ A n A l À + B ( A 2 ^ 2+ n 2 + 2 a A \ 2 + a n 2 - s ' Ô )+ C l n * a n ( n z+ À 2 ) ]
I a u À * B n { t + o 1 r 2+ ( 2 - " ) ^ ' ] } + c l t - ^ ' Ô * a ( \ 2+ n 2 ) 2 :f 0
de :
En posant le d.éterminantdu systèmed'équationségal àzéroon obtieni I'expression {

'Â:!-N'(l-u2)
'sEt

suivantes:
Où E et ^9sont donnéspar les expressions

R: (1 - r')^u +*{(r'+)')n - ( 2 + u ) ( 3 - u ) \ a n z + 2 À 4 ( 1- )'2na(7

+u)
t' ")-
+ (^2n2(3+ u) + nu - zn\)
)

,s- À2{t"'* À')+ 1-trn2)ft+ a(nz+À')'l+ i$' + (L- ")^'l

+ Aznz)ln'
Lr^'
,'r'I
-TI
sous la forme:
Dans le cas d.escoques cylindriques longues et minces { est donnée

( t - u 2 ) \ 4 + a { ( n 4 - n z 1 z+ ^ z l n a- ( n + u ) n 4 + ( 3 i u ) n 2 l (3-3)
ô: ^zn2(n2* L)
est donnée par
La contrainte critique méridienne sous compression axiale uniforme;
I'expression suivante :

ÀrIn 6
l( r-u2)\a+o{(nn-n
2t2 t I n+v)na+(3+ (3-4)
)T ")"'l\
ocr --
w @+ L) ëa
 Chapitre III aoilement iles coquescylindriquessous chargesde type silo. 78

L'analyse analytique de Timoshenko et Gere montre l'existance de "n" ond.esuniformes

autour du périmètre et "m" demi-ondes le long des génératrices. Ils ont d.éterminé les
contraintes critiques en fonction de la combinaison des ondes circonférencielles et longitu-
dinales. Le mode de flambement qui fournit la plus petite contrainte critique sera considéré
comme le premier mode de fla,rnbementpour la coque. Les résultats de leur analyse qui
montre l'évaluation de cette contrainte en fonction de "n" èt "m" sont représentéssur les
figr:res(3.2a) et (3.2b).

1.2 to6
-€- - n=2
---61---n=3
I lod -* -n=4
d
--o-n=5
ë 8 l0' *-n=6
q)
----rn=7
cr
_+_ n=g
k 6 rd ---r- E
o n=9
q) ---{-n=10
cl
k
rd

() 2 rd -- '
- - - .. - - 6 r ' -
----N---
-&---'-
----= -9-^
Ë-- v-

À=mn r/H
Fig (3.2a) Euolution d,ela contrainte critique aæialeen fonction d,e),

€
ë
()
=
()
h

a)

n (nomb,red ondescirconférencielles)

Fig (3.2b) Eaoluti,on iheIa contrainte critique axiale en fonction de n

A travers ces résultats, on constate que la plus petite contrainte critique correspond
 aoilernent iles coquescglindriques sous chargesd'etype silo' 79
chapitre III

,,n :2" et t'7y1:1", ce qui signifie, deux ondes autour du périmètre et

mod.e défini par
une d.emi-onde le long de la génératrice du cylindre'

Le mode 1", correspond.au flambement du cylindre coÛlmele flambement

,rr7-- 0" et "rn:
,les contraintes critiques
Eulerien d'un tube cylindrique. Pour des faibles valeurs de "771"'
des fortes valeurs de "?2" et "rn" '
évoluent trés rapid.ement a,vec"n". Par contre, pour
petite contrainte, et ne varient pas
elles ne sont pas trés grandes par rapport à la plus
sensiblement d'un mode à I'autre'
cylindriques
3.2.2. Analyse numérique de la stabilité des coques
sous compression axiale (calcul linéaire)

par le code d'éléments finis

L,analyse aux valeurs propres de flambement est exécutée
(Abaqus) selon le schéma suivant :
de charge-
- stockage de Ia matrice de régidité à l'étape correspondantà la phase initiale
ment sur la structurel
- application d.'une petite perturbation de charge active'

active est calculée, le calcul

une fois que Ia matrice de contraintes initiaies due à la charge
facteur de charge pour lequel
des valeurs propres est alors exécuté pour déterrriner le
l'instabilité est atteinte.
pendant le calcul de la matrice des contraintes initiales, Ies bords verticaux conserventles
cette contrainte est symétrique'
conditions de symétrie, car Ie pré flambement que provoque
Enfin' les conditions
Ceci s'effectue par I'intermédiaire de I'instruction "BOUNDARY"'
associéesaux charges d.ites "VMS" sont appliquées sur ces bords par le biais du mot
il
clé ,,BouNDARy,LoAD CASE:2". si cette d.euxièmecondition n'est pas explicite,
considère la même condition dans les deux cas'
cylindrique
(i) choix du maillage et des conditions aux limites de la coque
Pour ce problème
Le ccoix du maiilage est trés important dans ce type d'analyse'
et "rn" demi-
particulier, Ie mode de flambement présente "n" ondes circonférencielles
valeurs de "n'" et "rn" pour
ondes longitudinales. Le problème est d.onc de déterminer les
lesquelles on a la contrainte critique 6inirnale'
nceuds, ies dimensions des
Le maillage est effectué au moyen d'éIéments coque à huit
précisément :
éléments peuvent être appréciéessur la figure (3'3)' plus

- suivant la direction circonférencielle,on subdivise le segment en huit;

Chapi,treIII uoilernent iles coquescylindliquæ solrs chargesile type silo. 80

- suivant la génératice,24 rangéesont été adoptées.

Fig (9.3) Maihage d,eIa coquecyli,nilriquepour l'analyse de Ia stabilité

Compte tenu de la symétrie du problème nous avons considéré que le quart du cylindre
dat's Ia direction circonférencielle(f,) et la moitié du cylindre suivant la génératrice. Avec
cette approche, il est nécessairede mener plusieurs analyses,en utilisant différentes valeurs
de "@" et les conditions de symétrie ou d'anti-symétrie dans le plan moyen suivant que le
nombre d'ondes longitudinales est pair ou impaire.

Dans ce cas précis, nous avons fait variet "ot))entre (i) (#), à "n :2"
"t "orr"spondant
et "n:10".

Bord (1)

Ce bord qui sert du point d'application de la charge est simplement appuyé.

Bord (2)

Les conditions de symétrie pour les chargescorrespondent:

- aux déplacements radia.l et axial libres;

- au déplacement circonférenciel nul;

- aux rotations nulles autour de ce bord.

Les conditions d'anti-symétrie sont utilisées pour le flambement, ce qui se traduit par :
Chapitre III uoi,lementdes coquescylinilriquæ solrs chargesde type silo. 81

- des déplacements circonférenciel et axial nuls;

- des rotations non nulles.

Bord (3)

Les conditions appliquéessur ce bord sont des conditions de symétrie ou d'anti-symétrie

selon Ia parité du nombre de demi-ondeslongitudinales.

Bord (a)
Les conditions de symétrie sont appliquéessur ce bord pour les cha,rgesei le flambement.

3.2.3. Comparaison des résultats numérique et analytique de Timoshenko

Dans cette section on établit une comparaison entre les résultats déja obtenus par
Timoshenko et Gere et les résultats numériques, afin de valider le calcul au voilement des
coques sous compression axiale constante sur toute la hauteur. On utilise pour cela trois
coques d'une hauteur de 12.5 mèlres, d'un rayon de 2 à 4 mètres et d'une épaisseur qui
varie de 4 à 8 millimètres.

La figure (3.4) montre les contraintes critiques obtenues numériquement et les courbes
analytiques de Timoshenko et Gere. Les résultats obtenus avec les éléments coques à huit
nceudsvalident les résultats analytiques, I'erreur est de huit à neuf pour cent, ce qui prouve
que le maillage de la coque cylindrique et les conditions d'appuis ont été bien définies ainsi
que les charges (compression axiale constante). Ceci nous permet de poursuivre notre
étude à savoir le voilement des coques cylindriques imparfaites, raidies et non-raidies, sous
charges de type silo.
.......4ùqqB

E. a
B

o
o

F t Rayodc la oço 2oèm

rX . Epaiscu 6 millioètos

' llaeû l2J eètct

tEs2"lc5 MPa

' Næhr doda longiodbalcs

dtE
Alpha

pour r :2m
Fig (3./a) contraintescritiquesrnérid'iennes
Chapitre III uoilernent d,es coquescyli,nilriquæ solls chargesde type silo.
82

-Toebclto ...--.. Abe@

It4

176

l6t Æ
Ë t:
ti
E o
I ld, )
o '5
o
o a

É I Rayo dc la coqæ 3 dm

t4 '&airsæ6-;lli*ùË Â
'llôlæq l2J dtsË
136
t E=21o5 MPs

r NoEbE doldcr loEgiEdi@l€J 1

IB

d4 d6 dl2 dl4 in' dt8 tæ

Alpha

Fig (3./b) contrt,intes critiques mérid,iennespour r : Brn

_ liEdh6k^ Ab.qL,

,tI -Timobcoko .....-.Ab.qu

l!
tm

ll2 "l
"f
I Rsyoqdo la coqÛc4 mèæs

'llamu 12J oèæs

'&21o5
E lol MPa

o 96
o
"f
É
q
o
EE ' Rayo d. la coqû.4 Eèfts
o
c

5
"f
:iI
t Epais6
n millidras o

'!ls@q 12J drçg

. Eg2lc5 MPa
t Noke dodc loqiodiaahs I
a

tl4 da dro rA2 dt4 ''16 d7t tZt

"I
l(m
Alpb.
rlt

Fig (3.,4c) contraintes critiques mériiliennes pour r : nt

Fig (3.0 Comparaisons nurnérique et analytique d,escontraintes critiques nérid,iennes

Chapitre III aoi,Iementiles coquescyli,ndriquessous chargesil,etype silo. 83

axiale
Fig (5.5) Configuration d,éforrnéeile la coquecylind'ri,queparfai,tesous conLpresion

3.3. Analyse numérique de la stabilité des coques cylindriques

sous compression axiale (calcul non-iinéaire)

3.3.1. Prise en compte des imperfections géométriques

des
Au moment d.u dimensionnement d'une structure, la forme et l'amplitude réelles
structure
imperfections ne peuvent être connues. Elles seront connuesseulement lorsque la
forme et
elle-même est réalisée. Pour faire le calcul, il faut cependant considérer une
est
une valeur pout ces imperfections. La forme des imperfections que nous considérons
identique à celle du premier mode de flambement'
entre une tige
Selon l,Eurocod.e 3 [23], les a,mplitudesWs,, de I'imperfection sont mesurées
quelconque entre
droite (fig.3.6a) et un gabarit circulaire (fig.3.6c), disposé en un point
quelconque cercle
cordons de soudure respectivement d'une quelconquegénératrice et d'un
95%
directeur' La longeur de la tige et du gabarit est /' : 4tFrt' Mais sans dépasser
'est une tige de
de la d.istanceentre soud.urecirculaire ou entre soudure sur méridien. C
longueur Ir:25t que I'on doit utiliser pour les soudures circulaires (fig.3.6b).
sur le cercle
L'a^mplitude Wst des imperfections maximales tolérées sur la génératrice et
directeur est définie comme suit:

wor,-0.01/, soitlTor: fiot e fitt

à0.02/,
Chapitre III aoilernent des coquescylindriquæ sous chargesde type silo. 84

ï
T[**
TW-
"LW M
ifL
tr-
(a)
r
+ (b) (c)

Fig (3.6) Lirnitation des imperfectionsd,escoquescylindriques.

3.3.2. Etapes de création de la surface imparfaite [1]

En utilisant les notations de ia fi.gure(3.7), la création de Ia surface imparfaite se fait

selon les étapes suivantes :

.N
H
lc

t-
LIJ
c)
f
ttl

Fis (3.7)

a- On lit les valeurs des déplacementsunitaires du mode de Euler pour tous Ac; et Ay;;

b- On lit les coordonnéesde tous les points du cylindre initial (parfait) n; et gi;

c- En fonction de ot et !;, on calcul I'angle d pour chaquepoint avec : tan?; - l;

d- En fonction de Azl et Âgr; on calcul I'angle,6 pottt chaquepoint avec : t*fu: #;

Chapi,treIII aoilement des coquescylinilriques3ou8chargesd,etype silo. 85

e - O n c a l c u ll a r é s u l t a n t ed . e sd é p l a c e m e n t:sD ; - 1 f f i ;

-
f- On calcul le déplacement radial Dr; en projetant D sur le rayon : Dri : D;cos(0; 0;);

g- On cherche le point "j" où il y a ie plus grand déplacement radial Dr;; et on calcul "1"
telquerq-ffi;

Wstc représente la valeur maximale des imperfections tolérées.

i- On multiplie tous les déplacements du mode de Euler à tous les points par "4", on
obtient ainsi les valeurs des imperfections pour tous les points de Ia surface de la coque
cylindrique;

j- On ajoute les valeurs des imperfections obtenuesaux coordonnéesdes points de la surface

parfaite, on trouve ainsi les nouvelles coordonnées des points de la surface imparfaite.

3.3.3 Résultats numériques de ltanalyse non-linéaire

Le modè|e est traité par la méthode incrémentale, en utilisant la méthode de Riks.

Le maillage ainsi que les conditions aux limites utilisées sont celles du calcul précédant
paragraphe (3.2.2).

Une imperfection de la forme du premier mode obtenue au niveau du calcul linéaire est
considéréecomme étant Ia plus critique.

Les amplitudes maximales des imperfections géométriques utilisées sont celle définies par
I'Eurocode 3 [23].
*STATIC,RIKS"'
L'analyse non-linéaire est exécutéeen utilisant la procédure "

La figure (3.8) présentele déplacementaxial du næud en fonction du paramètre de la charge

appliquée ,,À". Elle montre que les imperfections influent directement sur la diminution
de 1a charge critique de la coque. Nous constatons que pour une valeur d'imperfection de
I'ordre de 0.008 la charge critique est réduite d'environ 15%.

La variation de la contrainte critique méridienne en fonction du rapport rayon/épaisseur'

pour des coques parfaites et imparfaites est représentéeégalementsur la figure (3.8). Nous
constatons que les contraintes critiques de la coquecylindrique imparfaite sont plus faibles.
Chapitre III uoilem,entd,escoquescglindriquessous chargesd,etgpe si,Io. 86

d . RÂyd do Ia €qe 2 ûè@

À ("
..6

s
ll
9u

o5lo15ûér354 &m19049fr
Dépl8eæûtuiÂl lml û
Fig (3.8a) Contrainte et chargecritique pour r :2m

l5
806
/i
/t s
3
lt:
s {i
:
.9

€oe
I
.2

0tlo15æÉm3t4 o€n59@697@ 19

DéPlaæænrdiâl lffiI

Fig (3.8b) Contrainte et ch,argecri'tiquepour r :3m

-- --- crq|F pâ.ÉÀib

t Rryù do Lc.q& 4 dtÉ

îæ
=
& E
o
'=@

€
I

;-.... --."
-.-... .q6

'Rsyod.hcoqr4[ùF

- Coqoo
.--...Coqt

tfotSûËp m fi 7@ 0 ru r@
Déplaæntuiallml rll

Fig (3.8c) Contrainte et chargecritique pour r : 4nt,

Fi,g(3,5) Infl,uence iles imperfecti,onssur la charge etla contrainte ile fl,ambage
 de type si,lo.
Chapitre III aoilement des coquescylindriquæ sous ch.arges d,

3.4. Analyse de la stabilité des coques cylindriques raidies et non-raidies

sous charges de type 36silo" (calcul non-linéaire)

Dans la première partie de ce paragraphe on présente l'étude analytique du voilement

des coques cylindriques imparfaites sous charges de type "silo". La solution' pour le calcul
des contraintes critiques méridiennes dans les silos est donnéepar I'EC3.
Dans la d.euxième partie on présente notre approche a.nalytique, pour le caJcul et la
vérification de la résistance des silos cylindriques imparfaits raidis longitudinalements,
à partir d'une étude paramétrique numérique bien approfondie.

3.4.1 Rappel de la solution analytique selon ltEurocode 3 [23]

(i) Cas des coques cvlindriques imparfaites non-raidies

Contrainte ultime m éri ilienne

o zrk '
c z r d ,: ^ , r n s: 7.7 (3 - 5)
"Yn3

Ot) 7-r est un cæffcient de sécurité.

L'équation (3-5) représente la formule gouvernante du voilement des coques cylindriques
imparfaites sous charges de type silo. Cette solution est adoptée par l'Eurocode 3 pour la
vérification au voilement des silos cylindriques.

ozrk: K"f y,k (3-6)

K":1 si À' < Ào

K":1
- - 0.6i|' ; I'l si Ào< À"< Àp
L À o* À o l
K,:3 si À'àÀp
^z

^,: ,l#i Ào:o'2; \p:'ffi

fu,k : est la limite élastique du matériau;

 Chapitre III aoilement des coquescylinilriques souschargesile type silo. 88

o"r. est la contrainte critique méridienne sous compression axiale constantel donnée par
Itexpression suivante :
Et

ozrc :0.605.8.1 pour / : 0.3

r
o est le cæfficient de voilement sous chargesde type silci (coque sous pression interne et
compression axiale variable le long des génératrices);donné par I'expression :

( 1- o o ) | - u " , = l
ot: oto+ (3-7)
Læ +PJ

(3-8)
t-o rr"

Ori :
P"(r) est la valeur minimalelocalede la pressioninterne.
@sest le cæfficient de voilement lié à la compressionaxiale donné par I'expression:

0.62
ao= (3-e)
L+ L.sû
(ry)''nn
tl: est un cæffi.cient qui tient compte de la non-uniformité des pressions donné par :

,b:#; J:0.2b

t; , 2 : \(f1- ô r ) - L ; t b t : 0 . 4
b
h:0.5V;;

(ffi) représentele rapport des actions de frottement à deux niveaux différents et voisinsl
tlt6 esl un cæffi.cient de réduction effective d'imperfection liée à la flexion globale de la
structure.
Wsk esl' l'arnplitude maximale des imperfections géométriques tolérées; donnée par :

+11
Wo* - ; Q : est la qualité du matériau
â{ î
Ch,apitreIII uoilement des coquescylind.riquessous chargesde tUPesilo. 89

(ii)Cas des coques cvlindriques raidies imparfaites

La contrainte critique méridienne de résistance seion les normes est donnée par la relation
suiva,nte :

@o"''P
o zrd.: (3 - 10)
7m3

Le cæfficient de voilement a est de la forme suivante :

si ffr0.,
IIa-0.65
Io:oo si ff.o.ou etfcoo
La contrairrLe o"r,, est donnée par I'expression suivante :

AtzAzs-Atsâzz À-= L AtzAts-Arr4'zs

-Tæn23 4-
= =-
- ! _--
ffi;Tn$-r (3 - 11)
acr,p:minimum de

Où ies quantitées An, An, Ats, Azz, Azs el A33 sont défi'niespar :

Att: (T)' + G*,(î)' i Azz: (i)' + G,e(Y3'

", ",

Ass:,, (T)n + D,a(tÏ)' (î)' + D,(î)- * #

Atz:(E,e* G,ù(vi) (î) i Azs:'+H i Art:? (T -rcç(T)'

Eç:3r?)*+ :3
i Eço | Ee:**
'.'1 Gt lb" I
G,o : Gvo:; * tl
Lï

Dç:ffi+*+*"#; Da:dl4
uEts , Gf f ô"*,.] * 9C" vç:-b
r)ço:@+ u i
6 L;*tl
Chapitre III aoilernent iles coquescylinilriques sous chargesile type silo. 90

où:
e", f et r4."sont les caractéristiquesgéométriquesde la section transversale d.u raidisseur;
"
6 est la largeur du panneau cylindrique;
ô" est la largeur effi.cacedu panneau cylindrique, comprise entre 6 et 1.9t f E
vil
t" est l'épaisseur fictive donnée pa,r la relation :

'
ts:t*+ (3 - 12)

"n" et "rn" sont respectivement le nombre d'ondes circonférencielleset longitudinales.

La formule simplifiée de la contrainte critique méridienne ocr,pt que peut supporter une
coque cylindrique impanaite sous les actions de type silo est donnée par I'expression suiv-
ante :

12
o"r,p: 8 1 " " 1 aoth6"zl;
E t
ffi*
Pour z :0.3

dcr,p: (g- 13)

ffi*o.6o5asû18!
Isef est la rigidité efficace (coque-raidisseur), calculée en tenant compte d'une largeur
participante "b"" de la coque.

thr: L dans le cas où le raidisseur est situé à I'intérieur de la coque.

,rr
{st: - dans le cas où le raidisseur est situé à I'exterieur de la coque.
L+#

3.4.2. Etude paramétrique par voie numérique (coquesraidies et non raidies)

Cette étude paramétrique pennet de tester I'influence de plusieurs facteurs inter-

vena,ntdans la résistance de la co{ue.

Les paramètres suivants sont analysés:

(1) Epaisseur de Ia coque cylinilrique, "t'

(2) Section des raiilisseu,rs,nAr'

Chapitre III aoilement des coquescylindriquæ sous chargesde type si,Io. 91

(3) Rigidité efi,cace, ul,"1"

(î"
(/1) Rayon d,eIa coquecylinilri'que,

(5) Largeur du panneau cylind,rique, "b"

(6) Nombre d,eraiiliæeu,r, 'h""

(7) La pressionhorizontale "Pn(z)"

(S) Le rapport, '1lt'

(9) Le rapytort, '%ltÇt"

Les dimensions des cylindres étudiés sont définis dans Ie tableau (3-1) suivant:

C aractéristiques Epaisseur[mm Rayon [m] Eauteur fm

Valeur minimale 4 4 r4
Valeur rrùoyenne 6 6 r4

Valeur maaimale I 6 14

cylindresà étudier
tableau(3-l) Dimension des d.i,fférents

Les panneaux cylindriques considéréssont larges, comme ceux étudiés par Syngellakis [65]
et par Watker [66]. La largeur est choisie de façon à casserI'ordre de la répartition des
ondes circonférencielles du cylindre raidi.
Pour cette étude les raid.isseursadoptés sont des profi.lésà section en T avec des sections
trarrsversales variant de 24 à 48 arf . Notons que pour une étude paramétrique, il suffi't de
faire varier la valeur concernée sans que les valeurs choisies correspondent nécessairement
à rp profilé laminé.

3.4.2.L Charges appliquées

Contrairement à toutes les études publiées précédemmentsur les cyiindres raidis et

non raid.is, nos cylindres sont soumis au type de charges réelies agissant sur les parois du
silo.
La détermination de la répartition des charges est faite suivant l'Eurocode 1. Rappelons
que les charges d'ensillages (pression horizontale et somme des actions de frottement)
Chapitre III aoilement des coquescylindriquæ sous chargesde type silo. 92

dépendent :

(1) du iliamètre de Ia ceIIuIe;

(2) de la profond,eu,r"2";
(3) des propriétés de la matière ensilée;
(il d" la ragosité desparois;
(5) de I'état de ûdange ou ile rernplissage.

3.4.2.2 Choix du maillage et des conditions aux limites

(i) Flambage linéaire

Après l'étude de la convergence,nous avons constaté que le maillage effectué au moyen

d'éléments coques à huit nceuds donne des résultats satisfaisants. Grâce à la symétrie du
problème, un quart de la structure est considérédans la présente étude. La discrétisation
est réalisée à l'aide d'un maillage exigeant 2904 degré de liberté.

Les conditions aux limites considéréesau niveau des bords de la struciure (fig 3.3) sont:

le bord inférieur (1) du panneau est encastré, on a donc pour chaque point de ce bord :

U":Uy:Ur-0r:0y:0"-0

Les bords verticaux (2) et (4) sont soumis à des condiiions de symétrie ou d'anti-symétrie.
Au niveau du bord (3) nous avons bloqué tous les déplacementset les rotations sauf Ie
déplacement suivant I'a-xedu cylindre, pour simuler la présenced'un raidisseur de renfort
à ce niveau.
(ii) Flambage non-linéaire

Le maillage est réalisé à partir du maillage précédent, en tenant compte des imper-
fections géométriques tolérées donnéespar I'Eurocode 3. Le calcui pas- à pas est effectué
avec pilotage paf, le déplacement vertical du næud situé au milieu du bord supérieur du
p?.nneau.
A chaque phase d'accroissementde déplacement,la position d'équilibre de la structure est
recherchée pax un processuriitératif suivant la méthode incrémentale de Riks.
 Chapitre III aoilement d,escoquescyli,nilriquessous chargesile t'ypesilo. 93

3.4.3. Résultats de ltanalyse analytique de ltEurocode 3

Dans le paragraphe précédant (3.4.1.) nous avonsiprésenté la formule gouvernante

du voilement d,escoques cylindriques imparfaites sous charges de type silo, données par
I'Egrocode 3. Cette formule d,épendet dans les deux cas (coque raidie et non raidie ) de
plusieurs paramètres :

- ca,ractéristiques géométriques de la coque cylindrique et des raidisseursl

- amplitudes d'imperfections géométriquestolérées;

- les charges de type silo (pression horizontale et somme des actions de frottement).

Da,ns cette section nous présentons les résultas analytiques liés à la solution donnée par
I'Eurocode 3. L'évolution du ccefficient de voilement et de la contrainte ultime méridienne
sont représentés sur les figures (3.9) à (3.16) Les charges considéréescorrespondant
aux états d.e remplissage ou de vidange centrée selon I'Eurocode 1. Sont données par les
expressions suivantes:

P.(r): Cs[1- exP{-c2@ -')}l; rQ):IoC t l t - e x p { - C 2 ( H- z ) } l d z

(i) Cas des coques imparfaites non-raidies

- - -* - - Coqw de 4 mètr6 de
-E- Coquede4mè:tresde ---+. Coquede6mèt6 de
-r- Coquede6mètcsde
--o- Coque de 8 mèt6 de
+- Coquede 8 mètr€sde

€
-ô

Ë o.+
o
o
o
E

'B o.o
lEl
q
r\

Epaisseur-t
"t" [mm]
Ëparsseur lmml ilt

(rft)
Fig (S.g) Eaolution ilu cæfficient ile aoilementen fonction d,eI'épaisseur(t) et d,urapport
Chapi,treIII aoilement des coquescglindriquæ sous cltargesd,etype silo. 94

o o
o o
€ Ë
I

3.F a^o
.ÊÀ .2À læ
o a
o
d

80

60

3æ4æ5æ@
ilt

Fig (3,10) Eaolution de Ia contrainte critique en fonction

d,el'épai,sseur(t) et du rapport (rlt)

(ii) Cas des coques raidies imparfaites

Fig (3.11) Maillage ile Ia coquecyli,nilriqueraidie

Cltapitre III aoilement des coquescylindriquessous chargesde type silo. 95

Ëtr
tE>
b-
o

l!

44551566J775
Epaissew"t" [run] Epaissar"t" [mm]

pour ns:4
Fig (g.12) Euolution d,eIa contrainte critique en fonction d'eI'épaisseur(t)

âE
Èà
o
cl
G

36
u/./rt

d,urapport(blrFî) pour n" = |

Fi,g(3.13) Eaoluti,ond,eIa contraintelritique en foncti,on
Chapitre III uoilement des coquescylindriquæ sous cltargesde type silo. 96

o
E
s

zÆ
.Y çT L?
,--2
b- .=À
q
o
o

4J55.566J7
Epaisseur"t"[m]

Fig (3.14) Eaolution d,eIa contrainte critique en fonction d,el'épaisseur(t) pour n" : tg

o o
o o

g I

gF
.=à
b-
o o

d €
R o

9
u/Vrt

Fig (3.15) Eaolutiond,ela contraintecritiqueen foncti,on

du rapport
h oou, ns :16
Chapitre III aoilement des coquescglind,riques8ou3 chargesile type silo. 97

lgt
o

o t60
rt
t,|{}

iF
.:e r,f)
.=À
e
o

124

Seaion du raidisseur [cm2l

(a) Contrainte critiqued'eaoilementpour ns:4,, t:4mrn'

o 2û

o
180
'9
160
--F
.:F 140
.ËÀ
b-
r20

€
rY

t2x2
Seaiondu raidisseur[co']

- L2, t: 4mm
(b) Contrainte critique d'eaoilementpour ns

220
o
I N
o
€ lEo
ro

160
Ë.E'
.2À t40
q

si r- ar.if**r' t",rft

-
(c) Cont'raint'ecritique de aoilementpour ns 16, t:4mm

Fig (9.16) Eaolution d'eIa contrainte critique en fonction

il'e Ia section iles rai'd'i'sseursA,
Chapitre III uoilement d,escoquescylindriquessous chargesd.etype si,lo. 98

3.4.4 Résultats numériques de ltétude paramétrique

L'évolution de la contrainte critique méridienneest représentéeen fonction :

- de l'épaisseurt figures(3.17)(3.19)et (8.20);
- du rappofi rft, figures(3.18)(3.21)et (8.22);
- du rappo* bltfrt figures(3.23)et
@.2\;
- de la section desraidisseursA, figure (3.25); ,
- du nombre de raidisseursn" figure (3.26);
- de la rigidité efficace/".1 figures(3.27)et (3.2S).
(i) Cas des coquesimparfaitesnon-raidies

l,o
o
q
€rm
ro

o- l@
3_jq

.2À
E- 80
o

8æ

40

41 5 s'5 7 7's t
"t
h"i.."l,, to,-l
Fig (3.17) Eaolution d,ela contrainte critique en fonction de I'épaisseur(t)

l4t)

o
E
u0

# l(It

::>
o to
q
6
@
rY
1{t

3æ ,100 5dr @ 7qt 8æ * l(xl)

,rt

Fig (3.18) Euolution d,ela contraintecritiqueen fonctiondu rapport (rlt)

Ch,apitreIII aoilernent d,escoquescylindriquessous chargesile type silo. 99

(ii) Cas des coques raidies imparfaites

160

o
E l4o
o
!

Eo.o
ær rm
.2à

9eo
G

F60

40
5.5 6 6.5 7 7.5 t

Epaisseur't"[rnrr]

(a) Contrainte critiquepour fls: 4, A, :24crn2

160 -.+-- Coque de4nètræ

o + Coquede6 nè'tres ,...::::t'*"'-"
8uo -- o-- - Coquede8 mètrEs

fuo
O-
zÆ ...---ï. "'r-/
Ë5 '. à
o
o
.E st

360
"'?"'
.......--i-""
44-555.56ô577.58
Epaisseur"t" [mm]

(b) Contrainte critiquepour ns - 12, A" :24c'm2

160 *-- CoquedeI Bet* A" ai"rettil

o +
â Coquede6 ne"o C"ai"re"{
o .- o-- - Coquede8 Eèù€3d€diarnètr{
€ t.o
ro

o: rlo
.Y çT
.2À
b-
o
tæ
i:::r-'a:'l
r""' d
;
Eso ........--i:.::.::..-..
-" ..--.-"i i
rX
60
.---.{r . " -
--ô-----
44.55J.566J17.58
Epaisseu't" [rorn]

(c) Contrainte critiquepour Ds :16, A" :24cm2

Fig (3.19) Euolution ile la contrainte critique en fonction d,el'épaisseur (t)

Chapitre III aoilement des coquescylinilriques sous ch,argesile tgpe silo. 100

o l6{)

l,O
ro

9?
%-
.YF
b- l(I,
o

tlt
n
60

5J.56ôJ77.5E

Epaisseur't" [mm]

(a) Contrainte critiquepour ns:4, A, : 48crn2

l&t

r@
o
€
B l,(,

O-
=Æ um
.=À
o
lq)

80

60
5J66.577.5r
Epaisseu't" [nun]

(b) Contrainte critiquepour ns : L2, Ar: 48cm2

180
o
o t60
.o

l,to

âE
.eÀ ne
o
o
tût

gt

Epaisseur"t" [mm]

(c) Contrainte critique I,our ns - 16, Ar: 48crn2

Fi.g (3.20) Euolution ile Ia contrainte critique en fonction ile I'épaisseur(t)

 ile type silo' 10 1
chapitre III aoilentent d,es coquescglinilriques sousclLarges

l60
o

o lO
€
tm
9?
.Yq
.2à t(xt

o
to

6{l

,|o
3æ ,rtx, 5æ 6æ 7@ t
,/,

(a) Contrai,nte critique Pour rt" : 4, A, :24ctn2

q
!
ro

gF
.2À

o læ

80

/m 50O 600 7æ 8m 900 1000

300
rll

(b) Contrainte cri,ti'quepour tus: 12, A" :24crn2

8tæ
o
E
.Ë l4o

âd no
::>
a 100

læ
û
5æ 600 7æ 8æ 900 lm
3æ 4æ

(c) Contrainte crilique pour rùs: 16, A'--24crn2

d'u rapport ('lt)
Fi,g (3.21) Eaolution d,ela contrainte critique en fonction
Chapitre III aoilement d,escoguescylindriquæ sorrsch,argesd,etype silo. 702

o
160
o

\a L&

zÆ t20
::>
o 100

80

60

300 4æ 500 6m 7æ 800 *ro loæ

(a) Contrainte critique pour rrs:4, Ar: 48crn2

.9 160
E
ro
É
g= ^r t *n
ÉÀ
.:z ,^^
r4v
bH
o
'5 rm
uBo

æ
300 4æ 500 00 7æ 800 900 l0æ
tft

(b) Contrainte critique pour fl" : 12, A, : 48qn2

180
o

o
! 160
'9
O- lto
r_s
::> t20
a

€ læ

80

300 rm 500 600 7æ 8m 900 lm

rh
(c) Contrainte critique?our n" :16, Ar:48crn2

Fig (3.22) Eaolution ile Ia contrainte critique en fonction d,u rapport (r/t)
chapitre III uoilement des coguescglindriquæ sous chargesile type
silo. 10 3

o
E
\a

âE lm
.=À
b-
a 80
6

o @

&
28n36&4brlrt48

(a) Contrainte critiquepour ns : 4t A, :24crn2

o 160
o
l,l0

âEt20
.Ëe
o læ

F
&!

9loll121314 15 16
U/{.t

(b) Contrainte cri,tiquePour rls :12, A, :24crn2

o
E
140
Ê

.2À
EEr20
b-
o læ
6

80

__24crrt'2
(c) Contrainte criti,que[,our fls : L6l A"
rapport (ul'Ft)
Fig (9.23) Euolution d.eIa contrainte cri,tiqueen fonctionilu
Chapitre III uoilement d,escoqaescglind,ri,ques
sous chargesd,etype silo. 104

o
160
o
!
ro l/t0
q-
=.* no
:Ê>
b-
o 100
d

80

60

28 12 36 40 u
ttf
(a) Contrainte criti,quepour rts : 4, A" : 48crn2

.E
E
160

--À1

.â 'à tan

o
'E rm
Q&)

60
e 10 ll 12 13 14 16
Tru"
(b) Contrainte critique pour fls : !2, A" : 48crn2

180
o

€roo
'9
l- ræ
ed
.g>
b- r20
o
'É
loo
u80

78erorrttÏ

(c) Contrainte critique pour ??,: 16, Ar:48crn2

F;g (3.21) Eaolution de Ia contrainte uitique en fonction ilu rapport (Ul'frt)

Ch,apitreIII aoilem,entiles coquescylindriquessous chargesile type silo. 105

O-
ed
E>
b-
o

t2 3rt ,o 4 4E
Seaion ùr raidissetr [cm1
(a) Contrainte critique pour r - 2m, t : Amm

r80
o

o 160
\o
t40
9?
.:>
.c! 120
o
a
lq)

2at24044
Section du raidisseur [cm1
(b) Contraintecritiquepour r - 3m, t:4ntm

t60
q

o l,ttt

ro

O-
um
2.S
.E> læ
o
a

€ g,

rY

Section du raidiss*" [.ott]

(c) Contrainte critique pour r :4rn, t -- 4mm

Fig (3.25) Eaolution d,ela contrainte critique en foncti'on

de Ia section des raid,isseurs"4""
Chapitre III loilement d,escoquescglind,riquessous ch.argesd,etype silo. 106

t,o
o
E rro
o

:E DO

9? tto

Ë€ 't
o
.f 90

,q st

70

Nosrbre de raidisseurs "ns'

(a) Contrainte critique pour A, = 24, t :6mm

150
o
l4{,
o

s tæ

E_Ft20
.2À lt0
a
l(It
d

90
rX

E'

168r0u2
Nombre de raidisseurs 'ns"

(b) Contrainte critiquepour Ar:36,

t60
o

E
F'.
'E l.rc
r!)

g':i tæ

=> l2o
q

IUO

5
€ ldt

Us
gt

(c)Contra'rrrr,
,:;i:'#:;5nr, ,': u**
Fig (3.26) Eaolution d,eIa contrai,ntecritiqueen fonction
ilu nombreile raidisseu,r'hr'
Chapitre III aoilement des coquescylind,riquessous chargesd'etype silo' 107

9160

'!4o
E
!f r:o
::>
!m
.E
Ëæ
o

Rigidiré effætive'IseI' [n'l

Fig (3.27) Euolution d,eIa contrainte Ûi'tiqli'een fonction

-- 36crn2
d,eIa rigid'ité effectiued,esraid'isseurs"Ir"1n pour As

#
.s>

E
Â

Z4 lor Zt lor ,2 lû, 3.6 lor r tO" {r lor

Rigidiré efrætive "Is€f [mI

Fig (3.25a) Eaolution de Ia contrainte critique en fonction

d,ela rigid'ité effectiue d'esraiilisseurc "1""1" pour A, : 48crn2

-lj.Pdz)
- lJ.Pn(?)

g "-"'P{z)

â É

(,'
I
/!
I

/i
f
t.,
Di.rÈF dô !r æq!. 4 dùE

Epairæc 6 nu'lliaràs
t Dirdùæ

'Epaiæq
d! L coqE

6 milltuùæ
6 ûÈ@

'o uô-r 8 tr 16

D€plâc66t ariÂl [m] DQlamatuid[m]

horizontale "Pn(z)"
Fis (3.28b)Eaolutiond,uparamètre d,echarge "^" en fonctionIa pression
Chapitre III uoilenent iles coquescylind,ri,quæ
sous chargesde type silo. 108

-Ë

æ
o

o
o

a
.Ë
o
FO

o
p
/i, _ Ae48
x

o
EP

-...---4e36

'DiaEèæd.lôc@6

tl0l5 5!0
Déplammtaxial [m] Déplacamt axial [-n]

Fig (9.29) Eaolution d,upararnètre d,ech,arge"À' en fonction

du dé,placementaxial pour les sections:A" :24crn2, A, = 36cm2 et A, : 48crn2

1.6

/i
T: 1.4
/,"'
€
z x
712
-5 r s i.
o
E IJ æ
o ,t/
$r i'/
//, €
o
o
o
o
/
I I €or
I
/
I
'2,-"' È-
æ16
e Eor
Æ (,' rl2 Ê

t, t'! t-*r6 -!Fl6

a B!4
1.......;i; 0.4 ----.--æ-12
. _ "' - - - ' -
OJ t qEr
l_._
,' ft
o2
' Diadto cb le coqæ 4 ' Disdtro ' DiaoèÈq do L coquo E
do ls qæ 6

lt0t5 tl0 246t

déplaccment axial [nm] Déplacement arial fnml Déplacemeotaxial [mm]

Fig (3.30) Eaolutiond,uparamètre

de charge"À" en fonction
du d,éplacement
arial pour un numbreile raidisseu,rs
ns:4t Ds: 12 et ns : 16
Chapi,treIII aoilement des coquescylindriquæ sous chargesile type silo. 109

3.4.4.L Effets des différents paramètres

a-Effet de I'épaisseur de la Paroi

L'épaisseur d.e la tôle a un effet primordial sur la contrainte critique. L'amplitude
des ond.esde flarnbement est trés grand,e,lorsque l'épaisseur est petite. Les figures (3.17),
(g.1g) et (3.20) représentent l'évolution de la contrainte critique méridienne en fonction
d,ece paramètre. On remarque son importancel la contrainte ultime est presque Ie double
Iorsque on double l'épaisseur de la tôle.
b-Effet du ra]ton de la coque
Ltinfluence du rayon est nette sur Ia contrainte et sur la charge que peut supporter
la
la coque cylindrique figures (3.17) à (3.30). Quand on double le diamètre de la coque,
contrainte ultime diminue presque de moitié.

plus
A valeur égale de I'inertie efficace"1""f", le raidisseur qui a une aire de section
(3'29)'
grande assureune résistanceplus grande au flambement du cylindre figures (3.25) et
L'influence de "A"" ne dépassepas les 40%.
d-
Les figures (3.26) et (3.30) montrent que la contrainte ainsi que la charge de ruine des
différents cylindres étud,iéscroît avec le nombre de raidisseurs. Le cylindre le plus résistant
est celui qui correspond à ]a coque cylindrique avec l-6 raidisseurs.
t'P-"
e-Effet de la pression interne
a
La pression interne a un effet stabilisateur pour la coque en compression. Elle
tendance à corriger les défauts liés aux imperfections géométriques(figure 3'28b).

B.b. Etablissernent d,une formule semi-analytique pour le calcul des contnaintes

critiques

Le développement de cette formule est basé sur I'effet et I'importance de chaque

paramètre étudié précéd.emment. Pour l'établir plusieurs hypothèses essentielles sont
adoptées :
Hypothèses :
(1) La coque cylindrique ainsi que les raidisseurs sont des éléments minces.
(2) Les raidisseurs sont considéréscomme des piècesrigides (c.à.d. qui sont indéformables
dans le plan de la section transversale).
Chapitre III uoilement d,escoquescylinilriquæ sous ch,argesd,etype silo. 110

(3) La déformation par distorsion de la surface moyenne (section transversale du raidis-

seur) est nulle (c.à.d. que I'a,ngleformé par les lignes de coordonnéesz: constante, reste
droit aprés déformation).

( ) Les hypothèses de Donnel [21] sont toujours applicables.

(5) On suppose que la pression horizontale P.(z) ainsi qu'une partie des actions de
frottement cumulées F(z) sont reprises par la coque cylir.drique. Le reste de la contrainte
de frottement est repris par les raidisseurs.

(6) La force latérale tra.nsféréede la coque cylindrique au raidisseur se fait par un cisaille-
ment caractérisé par un angle de distorsioî "'y^".

(7) Les contraintes dans les raidisseurs varient linéairement suivant la longueur.

(8) La loi du comportement est élastique et linéaire.

(9) Pas d.erisque de fla-bem"rrt'u.r niveau des raidisseurs.

L'hypothèse (6) découle de I'analyse du comportement d'une coque raidie verticalement,

soumise à une compressionaxiale qui varie linéairement sur toute la hauteur (figure 3.31a).
Cette analyse montre bien que le tra.nsfert des efforts, à la jonction âme-raidisseur-coque
se fait par cisaillement direct.

o
a
È

-9 -8 .n .lt -lo .t 0 I l0 -tæ .l@ {t o { .æ qa .224 :88 ./22 -Zzt -Zp .219 -zrE
Cotnirodridi* BNÉfl Cærnhæmbmeùo [tt{lb1 Cominc è cislfcd [klfrhrl

Fig (3.31a) Comportementile Ia coquerai,iliesous cornpressionatiale

Chapitre III uoilement d,escoquescylind,riquessous chargesile type silo. 111

ariale
Fig (5.3lb) Configuration défonné d,eIa coqueraidi,esous con'Lpression
coque
Considérons une coque cylindrique (figure 3.32) soumiseà une compressionaxiale. la
est supposéeencastrée à la base, et raidie verticalement'
Soient :
- H"la longueur du raidisseurl
- A"la section du raidisseurl
- 2H" la hauteur de la coque;
- b Ia largeur du panneau cylindrique.

Le chargement considéré va entrainer une contrainte verticale constante de compression

,ro",,, qui provoquera un raccourcissement "6"oqu"" de la coque cylindrique et t'6"tt du

raidisseur.
moyenne
En se basant sur les hypothèses citées précédemment, on peut estimer la valeur
de la force reprise par le raidisseur, calculéesuivant les étapes suivantes :

(i) valeur moyenne de la force reprise par le raidisseur

traduit
Le transfert des efforts par cisaillementdirect, à la jonction âme-raidisseur-coquese
par l'équation suivante :
^F* ( 3 - 1 4)
7:Jplt:T

où:
,rA, est la section qui participe au cisaillement, elle est donnéepar I'expression suivante :

A : 2.H,t
Chapitre III aoilement d,escoquescylind.riquessous cltargesile type silo. 7L2

L'effort repris par le raidisseur sera donc :

F*:2G1*H"t (3 - 15)

Avec :

(o E) - (F*H,l EA")
"H"l (3 - 16)

F*:o,A"llffi.rl ( 3 - 17)

(ii) Calcul de la contrainte ultime de compression dans la coque cylindrique

En accord avec le modèle (figure 3.32),la déformation axiale diminue dans la direction
verticale du bord supérieur du raidisseurjusqu'à la base. L'augmentation totale moyenne
de la déformation, est calculéeà mis-hauteur du raidisseurdonnéepar I'expression suivante

Ae:7-(y-#",) ( 3 - 18)

( â2,, \
représente la déformation dûe à la rotation du raidisseur vis-à-vis de
\#.")
I'excentrement e, des chargesappliquées.

ôw 6" 02u
: : a"
a, H,; d*
En utilisant I'hypothèse(7); que la variation des contraintesest linéaire sur toute la hauteur,
la déformation relative maximale (coque-raidisseur)s'écrit :

Ae^o,:T-,(#, -t"."") ( 3- 1 e )

Or:
6" Fro
^ F*.e?
a2e": =+
H": EA*; ffi

_+(* _
,,€^o.:rli *)]
Chapitre III aoilernentdes coquescylinilriquessous chargesd,etype silo. 113

Ao*o,:rlo"- F*(* - ( 3 - 20)

*)]
L'introduction de l'équation (3-17) dans (3-20) nous permet d'écrire la relation suivante:

- 2o,['
Lo,n.. - ({'- (3 - 21)
#\, {ffi-'})]

tOz Raidisseurvertical

!ll F*,
IIDi+'(c'que,
A\
t
ô qcoque;

Fig (3.32) ModèIe d,ecalcul d,ela distribution des contraintes

Finalement, I'expression de la contrainte critique d'une coque cylindrique, parfaite raidie

verticalement soumise à une compression axiale, variable linéairement le long de ia hauteur
s'écrit:

oz* ,'o*o,:oztr,,lr-(tt
6z*t,c.R.p: r})] rr- rrr
#\ {ffi*
Chapitre III uoilement des coquescylinilriquessous cltargesd,etype silo. L74

En relplaçarrt "or" par son expression (3-4), donnée par Timoshenko drns le cas d,'une
coque de hauteut 2H, et pour un nombre d'ondes circonférenciellessupérieur à un "(n >
L)", on obtient I'expression suivante:

ozril,c.Rp:Minimumd'e{-,*[t-'({t,#}t{#,ur.t})].

-L .ffi)\
t-(t-,') (w)'
n 2 ( n 2+ I )
(3 - 23)

Remarque: Le minimum de "o rrd,c.R.p)' correspond à n -- 2 et rn : 1.

Une solution semblable a été developpéepar J.C.Shang [26], pour le calcul de cette con-
trainte critique représentéepar I'expressionsuivante:

Formulation de la solution selon Shang [26]

K1K2r2E
ccr:
L2(7- r')(+)'

Kt:1.3 + 0.24!
rt

K -z : 1 - 1 . 7 5 f ( 5 0 9 - 7 ) l ' o , r , 7 < 5 0 0 ; K z : 1 p o u rz > 5 0 0

L krt I

arl :
I, : I0.9I"lbt3 si u :0.3
{btslLz(I - "r)}
La figure (3.33) représente la comparaison des contrajntes critiques, dans le cas d'une
coque cylindrique parfaite raidie verticalement. Les raidisseurs sont des poutrelles à section
transversale en T. L'ensemble (coque-raidisseur)est soumis à une coinpression axiale qui
varie linéairement sur toute la hautèur.
La comparaison entre nos résultats et ceux donnés par Shang, montre que notre solution est
plus proche des résultats numériques avec un écart de 13%. Ce qui justifie nos hypothèses
et permet ainsi la poursuite de cette étude.
Chapitre III uoilement des coquescglindri'quæsous chargesde type silo' 115

lffilo I'Rayondelacoque2mèE€s
l60
f -o-- euaris I ts."tionautaaiseurAs=zo2 ./
S l5o l--t-;;&qr" |'Notut"a"nidisorns=4,/.2
-
//
;a
6
'Ë
r,{,
a'
Bo
o
o
.S r2o

5 rro
lqt

EPaisseur[rnm]

(a) Euoluti'on ile Ia contrainte cri,tiqueen fonction ile I'épaisseur

+ Raym rlela coque2 mères

t@ ' Secriondu raidisseur
As = 24ct'
t Nomb,re !3 = 4
de raidisseur
t50
\
o lro

130
o
o
lm
6

o ll0

lû)

550 @ 650 ?æ 750

b/t

(b) Eaolution d,ela contrai'ntecritique en fonction du raPPort bft

t60

E rro \
î too
â ' RÂyû dela coque2 BèÈ€s
.F ræ r Secrimtlu rafttisseur z
o As = 24d
a t Nonhe dÊraidisseurns= 4
ew
I'
.$ rro
læ

(c) Euotution ile Ia contrai,nteffitique en fonctiondu rapport bltfrî

Analytique)
Fig (3.33) ComParaisond,es contraintes critiques (Shang, Abaquset
Chapi'treIII uoiletnent des coquescylindriques sous chargesile type silo. 116

Pour tenir compte des imperfections géométriquesinitiales, ainsi de la non-uniformité de

la somme des actions de frottemenl "F(z)" agissantssur les parois de la coque cylindrique,
on propose une épaisseur fictive (compte tenu de son influence sur I'augmentation ou la
diminution de la contrainte de flambage). Elle est donnée par I'expression suivante:

Tf":r*W (3 - 24)

où:
Wok: est I'amplitude des imperfections géométriquestoléréesdonnée par I'Eurocode3.

{ : est un cceffi.cientqui tient compte de la non-uniformité de l'effort axial (somme des

actions de frottement dans la cellule), calculé graphiquement, en se basant sur l'évolution
de la contrainte critique méridien+e (coque raidie parfaite et imparfaite) en fonction de la
racine du rapport (rlt), figure (3.34).

À partir de I'analyse graphique, figure (3.34), I'expressionde r/ est de la forme suivante :

7p 1 -0.1055
t' b : 7 . 9 3 4 1 - l P o u r 2 5 0-<t -1 < 1000 ( 3 - 25)
\t/

It0

* Equation(3-23)
170
* (Analytique IT=t+(W0k/b)l)
Æ
= 160
o
t5{,
\.V *Coçeimparfaite
k 140 * Chargement ayial variable
0
c) \..\ . iiiiofitîqù"rrE<wôrrutr
I
l30
t! \.
tRayon&lacoque2mètr€E'\..-
1m
o * Sectiondu raidisselr As = 24cm2 '..*_.-
t Nonbre de raidisseurns = 4
110 '.\--..\_
lq)
215 22
r/(rn)

pour r :2m
(a) Euoluti,onile la contrainte criti,queméri,d,ienne
Ch,apitreIII aoilernent des coquescylindriquessous chargesd,etype silo. I17

læ

* Equation (3-23)
tæ * (Analytique [T=t+(w0k/b)I)
ag

lt0
o
2

lm \â * Coqueimpadaite
o
o \ \r l*9.,1e*d.:**-t:
* ( Numériçe ç1=t+(W0klb)l )
\-
!Y 90 \.
'\
t Rayonde la coque3 mères \
.. '
t Sectiondu raidisseu As = 2zkm2 \\
U t0 ' E\
* Nontre ae raiaisseuros = 4 -

1a a1

./(rlt)

(b) Euotution de la contrainte critique méridiennepour r :3m

96
* Equation(3-23)
(l
* (Airalytique f[=t+(W0k/b)l)
88

o
(t t0
-...>{ i ff:ËH'"fl*'iu?ryq,e
o
\.- * (NuËériçetT=t+(w0Vb)l
)
o
' \\ .
fit
g tRayodelacoque4mères
-É-
t Sectiondu nidisseur As = 2,lcm- \ -.-
I Nonbrc de raidisssrr ns = 4 É- - - . . . . - - - . -
56

ara

.ittltl

(c) Euolution iles contrai'ntescritiques mérid'iennepour r : 4m

Fig (3.3/) Eaolution d'ela contrainte critique mérid'i'enne

en foncti,on de
Chapitre III aoilement des coquescylindriquæ sous ch.argesile type silo. 118

En iotroduisant I'équation (3-24)dans (3-23) tout en remplaçant '1" par "Ty"" on abouti
à I'expression de la contrainte critique méridienne de résistance de la coque imparfaite sous
compression a>cialevariable sur toute la hauteur.

n (t+W)
o zriI,c.R.I : Minimum de
,1ffi-r1
(r-,') (æ)'
L
Ir
(Tt
")2 (," - 1
1212
( 3 - 26)
n 2 ( n2+ t)
(r'+ 1)( æ fntf

)'

Le mode de fl.ambage de la coque cylindrique dans notre cas dépend de Ia pression interne.
11est donc essentiel de tenir compte de son influence, pour aborder le complément de la
contrainte critique méridienne (o r,ra,c.n.t).
Afin de trouver le lien entre les contraintes sous compression axiale variable sur toute
la hauteur et les contraintes sous chargescombinées(compressionplus pression interne).
Nous avons opté pour l'évaluation numérique de cette contrainte pour le cas de charges
combinées en fonction de la résistance de la coque cylindrique soumise à une pression
interne.
Le rapport des contraintes correspondant aux deux analysesnumérique et analytique (fig-
ures (3.36) et (3.37)) nous a permis de formuler le complément de (o,,4,ç.n.r) équation
(3-26) comme suite:

[aP^: - 1 . 0 8 3 .' " (;)-0'067 (3 - 27)

Fig (3.35) Configuration iléforrnée ile Ia coquecylinilrique raiilie sous chargescombinées

Chapitre III uoilement des coquescylindriquessouEchargesde type silo. 119

170 'Raya da laoquc2mèlræ Raym dela coque2 mètre

\ rNobrcdonidiwrs4 lE0 **t"ro.;6;*roJ
. Soti* ao onbscq AF24 cE: Seaim du nidisw AF48@2
r60 \
læ

t50 Equatior(l2o Equation(12Q

\\
6
'.. À160 \
àto \
o
'Ë ræ
\V o.\
9
6
.2
ro
o
'Ê '1
Coqæutsôcgæ e t4o
'1
\ Coqucrucharsæ
o \ cqæs
o \ \ MDEë .g \ \
Ë rzo
F
'.v
o.\

E ,ro
,q ll0 \t\ \ tm
tqt
\\ ll0
'ta
90

tt2 t2 vLa ll.2 a lLa

h.(r^) Pn*(r^)

136
vm I Rayo de la coquc 3 mètse t Ravo de la cæue 3 mèEe
tNcmbrcde nidissgns4 r Ndmbre de nidÎssu nF4
'seaiodunidis A-Aæ2 t t Sætion du midisw As4 o2
tt2

120
? 104 Equrio (126)
-
ù fruatio(3-26)

o96

'E \ 'É 101

8s Coqrcruchcgæ Coque rcu chrge
o o . c@DEees
. c@DE€
'fx
F80

oTz Uts

at
64

16 l7 18 l9
u tl 1314 15 16 17 lE 19
pn*(r^) pn*(r^)

t Raym de la coque 4 mètrs ' Rayon de la coqre 4 mètræ

rNmbcdemidiwns4 ' Nmbre de midissur as=4
I Særiæ <hrnidisq . Særim 6umidiw As=48 q2
As24 oz

Fqutio(}'2O €.F Equari@(!26)

F 1 \
= \
=70
o
o
a =
.Ë 80
\ o \ Coque su charge
o Coqwrughugæ o
g60 r @bE€6
r CûÉbE6
6
F
Ero
rY
50

:
z2 21 26 16 lE 20 2. 24 26
16 lE 20
Pn*(rô Pn*(rlt)

Fig (3.36) Eaoluti,on iles contraintes cri,tiquesméridiennes

en fonction d,u proiluit P" - (i) pour n": {
Chapitre III aoi,Iernentiles coquescylindriquæ sous chargesd,etype silo. L20

t80 r Rayo de Ia ouo 2 mètræ ' Rayon dc la qæ 2 mètsæ

I NdnbrÊ de nidrw Fl6 190 r Ndntre dc nidisæ æ16
170 i Sestim du nidiw As24 q2 I Sætio du nidiscur As4&m2
t80

160 Eqqatim(3-2O Equtioa(12O

,I! 170 \
c.
150
o 160
I
d
'É
140
a q \ Crquruchargæ E rso \ \ CoquesôEcbtrtæ
!)
\ \
æEDB6 g \\ gI9
d
tæ
Et.
n tlo
-

ll0
'o
t20

læ
ll0
728 lt3 t2 tLg 72 8 r3 9.6 ll.2 t2 na
Rt*(rlt) Pn*(r/t)

t Rayo de la coque 3 mètres l4 i Rayo de la coquo 3 uètre

t2a ' Nmbre de aidisw æ16 I Nmbre de midiw æ16
. SætiodunidisurAs24o2 I36 r Sætion du nidiw AsJ8 m2
120
t28
ô Equati<ia (3-2O
^ è Equatim (126)
>ll2
É l2o
; o
È lo4
.F rrz
\
o Coqæsosôagæ Coque rcw chage
9x . CmDm€ o
'e to4
, c@oEees
F

tes o
(J96

80
8E
û
80
16 l7 l8 19 14 15 16 t? t8 19

h*(r^) pn*(r/0

'Rayæ do h coquc 4 mètre t Raym de la coqæ 4 mètræ

rNmbro dc nidisqæ16 'Nmbrc de nidisw asl6
rSærio ùrakliwAs24d r Sætio du nidiw f,5d8 cn2
90

t Equatio (3-2i) Ftæ Equtio (126)

Ë =
o
9so
'5 90
o \ Coqusorschagæ o \ Coque su charges
a r cmDEa o . @mæ
'e
F7o
it
Èto
rY

60

16182022?,26 rrt lE 20 z2 24 26
Er*(rlt) pn*(r^)

Fig (3.3'l) Eaoluti,on iles contraintes criti,quesm,éridiennes

en fonction d,uprod,uitPn * (t) pou, ns : L6
 tou" 'ht T2I
Clua,pitreIII aoilernent d'es
"oqu"t "Vlind'iqutt
que peut supporter une coque
L'expressiea d{ffnitive de la contrainte ultime méridienne
à des charges de type silo est :
cyiindrique, imparfaite, raidie verticalement soumise

rIr
-1 \2
tnTf r
ffi
(3 -28

se simplifie et se met sous la forme

Pour les coques imparfaites non-raid.iescette contrainte
suivante:

I
t (t - r")
t
Im;rr1
\2H" J
:,2

H 12
(n - 1 ) I
TI
t-
,, (n, +1) (r'+ ; tn711
L )
ffi I

-1.088.P,(;)-"") (3 - 2e)

8.6. Comparaison des résultats et validation de ltapproche analytique

une comparaison des résultats

Nous avons présenté sur les figures (3.38) à (3.42)
le choix des paramètres est
numériques et réglementaires. Au sein d.ecette compa^raison
(3'a'2')' Mis à part que nous-nous
le même que celui de l'étud.e para^rrétriqueparagraphe
de quatre mètres dg diemètre'
sommes limités d.ansla présentation à la coque cylindrique
par rapport aux coques de six et huit
Elle présente une grande résistance au voilement
mètres de diamètre.
entre nos résultats et les résultats
ces comparaisons montrent une bonne concordance
surestime certaines valeurs de Ia
numériques (Abaqus). Par contre la formule de I'EC3
pour des épaisseursde quatre à sept mil-
contrainte critique. Cette différence s'accentue
ad'optée par la norme n'enveloppe
Iimètres, par conséquent le choix de l'épaisseur fictive
la pratique'
pas tous les cas possibie d'épaisseursrencontréesdans
rédaction de la norme EC3' partie sIlos'
cette solution sera ainsi proposéeau comité de
Chapitre III aailement des coquescyli,nilriquæ sous ch,argesde type silo. L22

+Erocode3
c)
--- Abaqns
-..o-.- Eqruiion(129)
É
c)
140
I Coque de 4 mètrcs de dianètre

cr-

.Ëâ
e lzo
c)
c)

44555J66J7758

Epaisseur
[mm]

(a) Eaoluti,ond,eIa contrainte critique en fonction de l'épaisseurt

rls +Eurocode3
<' - Abaqus
c) . - -o- -. EquaTon (3-29)

6 140
t Cooue de 4 mèrre de
(u-
Éo;
.f; E tzo
o
q)
\.-r.-,
roo *,*..*-
E
-.*.*,_.*Ë.È:
t0

250 300 350 400 .t50 5fr1

rlt

(b) Euolution ile Ia contrainte critique en fonction ilu rapport rf t

F;g (3.38) Comparaisonnurnériqueet réglementairedes contraintes

critiques m ériili,ennes( co quenon-raidie)
Ch.apitreIII aoi,Iementd,escoquescylindriques sousch'argesd'etgpe silo' 123

o
.€E
a

â-d'
:=à
()
o

È
(J

Epaisseur[mm]

(a) Euolution de la contrai,ntecritique en fonction de I'épaisseurt

Êro +Eurocode3
<. - Abaqus
E l5o - -.o-'- Equation(328)
P l/o

âE''
'Ë à rto
o
c 110
'E rm

5eo
rlt

(b) Eaotution d,ela contrainte critique en fonction du rapport rf t

o
\-e +Abaqus
a.- f,quâf,qq
E 150 -\\ '.. -. -o-' ' Eurocode3
a
l,O

âE'
= > lm
130

b-
o u0
. Coquede4 mèbesde diaûèt(e
tÉ l0o
t Secdondu raidisseurAs=2'lcmz
90 t Nombrede raidisseurrs=4

262A4.n
b/ftt

(c) Eaolution ile la contrainte critique en fonction ilu rapport UltFrt

Fig (9.39) Comparaisonnurnériqueet réglernentaireiles contraintes

:
critiquesmériiliennes (coque raid'ie)pour A, :24crn2, tus 4
Chapitre III aoilernent iles coquescylindriquæ sous chargesde type silo. 124

o
180
+Euroc<lde3
o t70 < - . Abaq8
-o- -- Equatio (3-2E)
1 @ --
!
a
t50
Ë.Fl,O
Êà
130
o
120 r Coquede4 mèE€sde diamètre
t Sectiondu raidisseurAs=48cm2
Ê, 110
t Noobre èraidissurr=4
t00
4s t t'rpui.r.L 7 75 8
rro,i'

(a) Euolution ile Ia contrainte critique en fonction d.e I'épaisseur t

+Eurocode3
r7o <- - Abaqus
€ l@ . - -.- -- Equation(128)
\)
âs
150
'\l----
l,l0
ÊÀ -\' ----*-t

o 130 "'S.- -t-*'a..

\-.
.--.--------..
o
r20 ':aa'\
. Coque de 4 mètres de disEèbe
.-
110 rsectiooduraidisseurAs=48cnf, """'\.- -':\"....-,...-
.Nombrederaidbseurns4
t0o ,

rft

(b) Euolution de Ia contrainte criti,queen fonction d,urapport r f t

180
o

o
€
t1i

lo
+.-
::.\ l=ffiv:
l.--"-.-Equation(3-28)
a
\) .,\)
150
âEl,o
ÊÀ " \ . . . ...\.....-
.\
o 1æ ^_.\-
o
tû
t Coque& 4 EèûEs& dianète ' *'\".
6
.!l '_'\-.--.-
ll0 . Sectiondu raldbseurAs=4Ecnf
' Nombrede raidbseurns=4 Êi :..--
t00

4
b/vil

(c) Eaotution d,ela contrainte critique en fonction ilu rapport bltfrt

Fig (3.40) Comparaison numérique et réglementaire iles contraintes

criti,quesrnéri,diennes(coque raid,ie)pour As - 48crn2, fls : 4
Chapitre III aoilernent des coquescglind'riquæsous chargesile type silo. t25

lE0
o
+Er.nocode3
170 | -4..'
o ---4buq,o | -1.r'-/
'El l@
a

âE'l &
150
1..-.-"
....2'
.âÀ 1æ
.d ."
o _.t/
q)
120 -j-'
'
-- - '.': : t Coquede4 mèEesde diamètre
n0 t SectionduraidisseurAs=2rkm2
-..-'::'y
t Nombrederaidisseurm=16
o l@ -..::?
4555566J1758
Epaisseur[mm]

(a) Eaolution de la contrainte cri'tiqueen fonction d,el'épaisseurt

o
+ Eurocode3
o <-. Abaqus
'tr - - -o- -. Equation(3-28)
a
trr

âE'
' Ë E'
€ rzo
E rro
E i Sectiondu raidiss.to 65=?Acm'
6 1oo t Nombre & nidisseu ns;16
250
tlt

(b) Euolution d,ela contrainte criti"queen fonction du rapport rf t

o +Eurocode3
170 <.. Abaqus
o
E :'.)* -. -c- -. Equation(3-2E)
o
t@ '*.
r50

âE
:2à
l,l0 *..
o
1æ _..,\...
o
lm '
r Coquede4 Eètres ds dianètre ' - \:
B 110
r Sectiondu raidisseurAs=24cnf
o t00 t Nombre de raidisseurns=16
o
71
b/vil

(c) Euolution ile Ia contra'inte critique en fonction du rapport bltÇi

Fig (s.41) comparaison nurnélique et réglertento,iredes contraintes

critiques méridiennes(coque raidie) pour As:24cm2 , ns : 16
Chapitre III uoilement d,escoquescylindriques so,tr.s
chargesde type silo. L26

o
q)

ro

âE
.Ëë
o
o

U
4 7 7s 8
rpars"or lmmjr

(a) Euolution il,e Ia contraintecriti,queen fonction de l'épaisseurt

(b) Euolution d,ela contrainte critique en fonction du rapport r f t

O * Coqæ de 4 mètes de dianèùe

o . Sectionôr raidisseurAs4Ecm
tE0 r Nmbre deraldisseurns=16.
I r70
t60
âE
. Ë ë 150
o 140
o
tæ
120
()o u0
7 73, t EJ
bFvn

(c) Eaolution d,eIa contrainte critique en fonction ilu rapport b/t6t

Fig (3.42) Comparaison numéri,queet réglementaire des contraintes

critiques mériiliennes (coqueraidi,e)pour As:48cm2, n" : 16
Conclusiongénérale 727

CONCLUSION
GENERALE

L'objectif de ce travail est d'une part, I'analyse de I'effet de I'uniformité et de la

non-uniformité des pressions horizontales dues à la matière ensilée sur le comportement

d'une coque cylindrique parfaite non-raidie en élasticité linéaire, et d'autre part, l'étude

de la résistance au voilement des,coquescylindriques raidies verticaLement, sous charges

axisymétriques (pressions exercéespar la matière ensilée dans le cas d'un remplissage ou

d'une vidange centrée).

L'étude analytique menée dans la première partie montre la distribution des sol[citations

dans les parois de la coque pour le cas de charge axisymétrique. Les résultats analytiques

de notre approche montrent une bonne concordance avec les résultats numériques obtenus

en utilisant Abaqus.

Cette analyse de comportement est complétée par l'étude de l'effet des charges non uni-

formes sur la distribution des coatraintes. Nous avons ainsi modélisé la structure dans

Abaqus en considérant d'une part les charges préconiséespar I'Eurocode L et d'autre part

les pressions données par la base de Chartres. Les résultats obteuues montrent une flexion

importante des parois, localisée à un niveau donné du silo. Par conséquent le comporte-

ment membranaire évoqué dans le cas de chargesaxisymétriques n'est pas représentatif de

la situation réelle du silo.

Dans la deuxième partie concernant le voilement, nous avons choisi les actions préconisées
Conclusion générale L28

une étude paramétrique par voie numérique sur Abaqus. Les résultats de cette étude,

sont utilisées pour établir une formule semi-analytique de calcul de la contrainte ultime

résistante au voilement des coques raidies et non raidies longitudinalement. Elle fait in-

tervenir toutes les caractéristiques géométriques qui participent au niveau de la rigidité

prend
d.e I'ensemble (coque-raidisseur), parmis lesquelles figure une épaisseur fictive qui
actions
compte, l,amplitud.es des imperfections géométriquesd'un" part, et la variation d'es

de frottement (frottement de la matière contre la paroi) sur toute la hauteur de la coque

d'autre part. La variation de la pression interne, suivant les formules généralesdes normes

est bien évidemment intégrée dans nos calculs.

La formule analytique que nous avons ainsi établi, peut être utilisée dans le calcul à l'état
et
limite d.e la stabilité d.essilos cylindriques raidis verticalement, encastrés à Ia base

chargés symétriquement. Cette formule a I'avantage de tenir compte de la nature réelle

des charges agissant sur les parois de la structure'

et
En perspectives, nous envisageonsun grand développement de cette formule,
silos avec
I,extension de son champ d'application aux silos reposant sur des poteaux, des
par exemple) et
trémies, des silos soumis à des chargesdissymétriques (vidange excentrée

des silos raid,is dans les deux sensavec la prise en compte des contraintes résiduelles.
 Bibliograpltie L29

BIBLIORAPHIE

[1] anoruE F.: SiIo à parois cylind,riques raidis, prem,ièreétape aers I'étude de
I'instabilité . Mémoire non d.iffusé,à consulteraju S.E.S,.C.E.B.T.P. Saint-rémy
(Mars r98r) p.134.
[2] e.c.t. coMITEE BrJ:Projet denormesaméricaines: Recommend,ed, praticefor designe
and constructionof concretebins,silos and,bunkersfor storing granular material.
A.C.I. Journal,10 (1975).
[3] aluontH B.o.: Infl,uenceof edgecond,itions on the stability of axially cornpressed
cylindricalshell. AIAA Journal.vol.4(Jarrvier1965)p.134.
[4] ennocz J., BABCocKL.w:,,Theeffectof generalimperfections on th,ebucklingof
cylindricalsh,ell.Trans.of ASME (Mars rg6g)p.28-38.
SCHLER.:On the bucklingof stifued imperfectcylindricalsh,ells.AJ;
[5] ennocz a.nd,
A.A.JournaJ,11. (tgZ6).
[6] nennn .1.: Stractures coques" Équation généralesde stabilité ". E.N.S.A.E.
Toulouse. (r98o).
[7] ertcuT G.E.: Defectsin Accepted, Methoilesof Estimating Design Loadingsfor silos.
Proceed.ings,Institution of Civil Engineers, part 1 vol. bf 88 P.1015-1036.
[8] nltcHT G.E.: Temperature - induced load,ingon silo walls. Structural Engineering
review, vol. 4- 1, (r9gz), p.61-71.
[9] nlrcHT G.E.: Responseto loading of tree sirnilar grains silos. Bulk solid handing vol.
5- 2,,(tggz) p. 241-2a3.
[10] nucnT G.E.andcoENERf M.: Effectsof aaring hori,zontaly,ressul:ein cylind,rical silo.
Powde handing and processingvol. 5- 4, (November lgg3).
[11] nr,rcnT G.E.: Design implication of meo,sured,pressuresand stro,insin silo. ASCE
Journal of Structural Engineering, vol. 118- 10, (tggr), p.2720-21742-
[12] nouooN A., KHELIL A., DoREMUSJ.p.,RorH J.c.z Mesure ile déformation d'un silo
métallique Ttendantles ilifférentes ph,asesil'ensilage(Comparaisons entre les mesures
opti,queset les rnesures par jauges d,edéforrnati,ons).Construction métallique.vol. 3
(tggo) p.25-31.
[13] unrtrsn coNSTRUcTIoNAL srEEL woRK AssocIATIoN: lâe th,eory and pratical ile-
sign of bunlcers.LTD -32 (1968).
 Bibliographie r.30

plate and'shelld.McGraw-Hill (1975).

[14] anusg D.o. er ALMoRjrnzBuclclingof bars,
fVe CongrésInternationalde Mécaniquedes
[1b] cequor A.: La pressionilansles silos.
sols-Londres(1957)P.191-195
tlteorg.Prentice-Hall,Inc. Englewood
[16] cnatES A.: Pricipalesof structuratstability
Cliffs, New JersevGgZù.
A. : Calculliné'aireet non li'néaireil,es
[12] cn.lveNT c., coMBEScuREA. er IIoFFMANN
coques. Ra,pportD.E.M.T./S.M.T.S./7g/22,F.E. 5423cours I.P.S.I(Janvier1979)'
Sécuritas et de la SOCOTEC- Path'ologie d'e
[18] oeussy R.: Leschroniquesd,uBureau
la construction métallique.AnnalesI.T.B.T.P.,378,(Dec 1979).
[19] oesrrryNDoR d,escoquesminces.MASSON.Paris Milan Barcelone
p.: Mod,étisation
Mexico, 1ggo.
pourle d,i,mensi'onnernentetla constructi'ondessilos'
[20] orN fi.622/4;Directiaegénérales
('gzs).
poNNELL L.H: A new theorgfor the bucklingof thi'n cylinilresunder aaial compress'i'on
[21]
anil bend,ing.Transaction, Am. Society Mech. Eng, Vol.56. (tgs+)
of imperfecti,onon bucling of thin cylind'resand
[22] DoXNELL L.H., wAN c.C.: Effect
Transaction,A.S.M.E, vol. 72, (tgso).
colutnns und,er aæial cornpressi,on.
sh'ells. Comité Européen de Normalisation'
[23] nUnoCODE 3: Design of steel
SILOS)'
DRAFT prENV(ree3-4-1
CE}r'/TC250/SC3.
action on structures. Pa'rt 4: Action in silos and
[24] eunocoDE r: Basis of d,esignand,
tanks. CEN, ENV (r99r-a).
- ^tilos cylinilriques ouaertset ferrnés soumis à Ia
[2b] essuNGER AHMED er SCHRoEDeR:
Trad'
charge stationnaire ilu aent. Traduii de DER $1a.hlf6.s,12lQg7t-), P'361-368,
v.T.T.72
/L.
of steelbins for Storage of bulk Soli'd's'
[26] cavloRD E.H, GAYLoRDC. N.: Design
prentice-Hall, Inc., InglewoodCliffs, New Jersey07632,(tg8+).
coques.Traduction,Eyrolles(tgg+)'
[27] nf,ÛCGE\À/.: Statiqueet ilgnani,qued,es
in silo'
auf d'ieDruclcaerhô'ltnisse
[2g] nervreEE., KAMINSKIM.: Der Einfl,ubexzentrischer
Boutechnik61, (tg8+), P.73-81.
calculpour les chargesd'euoilement et les
[2g] rfnneERT K.H.: Proposi,tionile basesile
ch,arges limite ilescoques. Der Stahlbau5/(1966),Trad. C'T'I'C'M' (Juillet Lg77)'
anil numericalstuily of tlt'e effect ini'-
[80] nutcHINSoN J.w., BUDIANSKvB.; Analgtical
of struc-
t;al inperfections on tlte inelasti,cbuclclingof a cruciform colwnn Buckling
tures,, edited by Budiansky,springer-verlag.New York. (tgz6).
a, (Avrii 1973),p'156'
[3f] lvrvreENcAn.: .jilo en acier.Acier-Statrl-steel,
 Bi.bliograph,ie 131

[32] leNssEN H. A.: Versuche û,ber Getreidèdruck in Silozellen. Zeitschrift Verein

Deutscher Ingenieur. vol. 39 (Août r8gf) p.1045-1049.
[33] lnNc ET AL: Buclcling of longitud,inally stiffened cglind,ers. Journal, Struct-Div.,
Proc. of 4..S.C.E., (Oct rg6+).
[34] xnnlll, A.: Etud,e du charnp de uitesse et de contraintes dans les silo rnétalliques:
Exploitation et développement des résultats de la base expérimentale d'essais de
silos, en vraie grandeur de Chartres. Thèse. I.N.P.p. Univ. de Nancy (tg8g).
[35] xotreR w.T.: Buckling and,post-bucklingbehaaiourof cylindrical panel under axial
compression. NLR, Rep. S 476, Rep. Trans. Nat. Aero. Res. Inst., vol. 20. (Ma
1956).
[36] xoIrnR w.T.: Elastic stability and post-buclclingbehaa'ior.Proc. Sy*p. non-linea.r
'Wisconsin
probl., University of Press, Madison. (tg63)
[37] xoItnR w. T.: General equationsof elastic stability for thin shells. Proc.Symp. on
the theory of shells to honor Loyd Hamilton Donnell, Univ. Houston, (tg6Z) p.187-
227
(en holland,ais).
[38] xoItnR w. T.: On the stabilityof elasticequi,Iibrum Universityof
Delf thesis. Paris, Amsterda* (rg+S) (Traduction en anglais par NASA Dept. TTF-
10-(re67)).
[39] xottoR w. T., HUTCHINSoN theory.App. Mech. Rev., vol. 23,
J.rûv.:Postbuckling
( t g Z o ),P .1 3 5 3 -1 3 6 6 .
[40] r,eNcHAARH.L.:Energymethod,e in applied,
mechanics. Wiley, New York. (1962).
[41] t ovn A.E.H.: A treatfueon the rnathernatical th.eoryof elasticity. 4th ed. Dover,
New York. (rgrZ).
[42] uonnAu G., soKol L.: Conceptionet calculdes silosen acier. Rapport n7.003-2,
C.T.I.C.M.(Avril leze).
et
[43] N.v. os: Règlesiléfinissantles effet de la neige et ilu aent sur les constr-uctions,
annez,es.Pa,ris.
Réviséesen (tSZZ)
[++] NeooAF osKouEI A.R.: Mod,éIisation par élérnentsfinis des actions des matériaux
granulaires s'trr les parois d,e silos. Thèse, LPMM, ISGMP, University de Metz,
(Juillet 199b).
[45] Nonun ExeERJMENTALES FRANçAISES: Constrtrctonntétallique,SiIo en acier, Cal-
cul iles action ilans les cellulesP22-630,nouvelleédition, (tggr).
[46] onv H., HoFFMANNH.: Stabilityof corrr,gatedshells.ECCS colloquium on stability
of plate and sheli structures,Gheni University,6-8 (April 1987).
ilescoques.Rapport interne,
[47] eotrcR-FERRYM.: Infl,uenceilesd,éfautssurle fl,arnbage
LPMM, ISGMP, Metz, (1986),47 P.
 L32
Bibliographie

projetz Nouuelle aersion ile la Norrne Allernand,eDIN 1055 -teil 6 sur les ch'arges
[4g]
il'ensilage.Trad. C.T.I.C.M, (Avril rgSg)'
les silos'
Manif,festations d,es surpressi,onsd,e ui'd'anged'ans
[49] nUVfBER:f M., A.:
au rnouen de tubes rnétali'ques' Acier-Stahl-
réparation pratique d,e silos acci,d,entés
Steel. 9. (1963).
nouaellessur les efiort's exercés par les matières ensilées
[50] nerrvrBERl M.: Recherche
r$a3).
sur les parois iles silos. Annales I.T.B.T.P. (Mai
théori,eet pratique. Ed. Eyrolles' Paris (tgzt)'
[51] nnwrBERT M., A.: si?os,
Edition
traité théorique et pratique d'u calcul iles si'los'
[52] REIMBERT M., A.: Silos,
Eyrolles. Paris (tgZZ).
load'
ilesign of unstifieneil cylind,rical shells subject to axial
[53] nomER J.M.: on the
of Plate and Sheil Structures'
anil internal pressure.Proc. ECCS Coll. on Siabitity
(Aprih987)' P.539-548.
pccwnntz s.: SheII Stability Hand'book' Crown House' Linton
[54] sevruplsoN L.Â.,
Road, Barking, EssexIG11 8JU, England, (tgg')'
on buck-
serrnuELSON L.Â,.: Buckti,ng of shellswh,i,ttt'local reinforcements' Int' coll'
[55]
Lyon, France' I7-I9 Septem-
ling of shells structures, on land, in the sea and' the air,
ber, (rg9r), Elsevier, London, P'401-408'
méthodesd'ecalcul
si,Iosen béton- Comparaisond'esd'iuerses
[b6] SCffVfOLJ.: CaIcuI d,es
d'ouvrage particulier en
aaec application nurnérique. c.A.c.T., séminaire : Etude
béton - (Novembre r98o).
ilestabilizing
: Buclcling of irnperfect sti,ffeneilcyli'nd'ersund'er
[57] sunrNMAN, SIMITSES
12' (Dec 1977)'
loails i'ncluiling torsi'on' A'J, A'A, Journal, vol' L6''
ensilées'
soror l.: Méthod,e ile calcul pour l'éaaluation d'esacti'onsdes matières
[$]
Revue d.econstruction -{la.llique 2, (rg8g).
stringer reiforced'
sensiti,aityof axially corn'pressed'
[bg] STOeIIENSW.B.: Imperfections
Journal, vol' 9' (Sept 1971)'
cylinilrical panels un6er internal pressure.A.I.A'A,
P.1713-1719.
sh'el(s'Int'
A.C. : Elastic buckling of stiffened',cylindri'cal
[60] svNGELLAKTSS., WALKER
(tgZ8)'
Sy*p. On intergrity of offshore structures Glasgow'
stiffened
WALKER A.C. : Elastic local buckling of longitudinally
[61] svNCELLAKTS S.,
and components' Cardiff'
cylind.ers. Cof. On stability probiem in Eng. Structures
(tgz8).
 Bibliographie 133

[62] runxDRUp p. : Buckling of unstiffened a,nd,ring sti,ffenedcylindrical sltells und,er

axial compression. Int. Journal, Solid and Structures, vol. 9, (tgZg). P.671-873.

[63] truosHENKo s.p., GERE J.M.: Theory of elastic stability. I-nternational student
Édition. Mac Graw Hill (1959)

[64] veln n.: La mécaniqued,esmilieux continus et le calcul d,esstructure. Collection

D.E.R., E.D.F. Eyrolles, Paris vol. 31 (tgSg).

[65] welxER 4.c., sRIDHARANs.: ,4.zalysisof behaaiourof azially cornpressed, stringer-

sti,ffeneilcglinilrical she\Is.Proc. Inst. Civ. E g., part.2, 69, (June rg8o).

[66] welxER 4.c., sRIDHARANS.: Buc&Iing of cornpressed,longitudinally stiffened cylin-

drical shells. BOSS 79,2nd. Inetrnational Conferenceof behaviour of offshor struc-
tures, Imperial College, London. (28 to 31 August LgTù.

[67] wnllER, sINGEn: Further eaperimentalstudie on buckling of integrally ri,ng stiffened

cylindrical shells under axial compression.Experimental Mech. (July ry7ù.

[63] vervreKl N.: Elastic stability of circular cylind,rical sh,ells. North-Holland series in
Appiied Mathematics and Mecha.nics.Amsterdam-New York-Oxford. (i98a)
AIïIÏEXES
Annexes 734

Annexe 1

Solution analytique de l'équation différentielle du quatrième ordre

[équation (2-6) Chapitre II]

On va chercher la solution de I'expression(2-6) on utilisant les étapes suivantes :

L'équation (2-6) est de la forme suivante:

(#++onv,):
*l'".."+]
où:
P " ( r ) - C s { 1- e x p l - C z ( H - " ) ) } ; fQ):Ct{1 -expl-C2(H -
")l}
On pose:

1 | , t : U l ? 1 1 , " : ' L t r 2 i l u z , z - - U g ; ' L t r B , z : ' 1 t r 4 i ' t l 4 , z : O t L ' l . l+L Q Q )

QQ):
*l'"ur-'+]
On peut écrire cette dernièreexpressionsousla forme suivante :

0 100\/"r\

*l::l :[ 3 3âîl[;:)*a("){3)
^/"t\ / /0\

\u4/ \-aa o o o/ \";/ \r/

On pose:

,:Gl)
Anneaes 13 5

Donc:

AU :
AU + QQ)B
A"
où:

A_
(j,itil; B_
fi)
On résoud d'abord. le système homogène # - AU qttt a coûlme solution :

U(z) - Us ù'exP(Az)

( Ut(r): ( I / o ) r e x PÀ 1 2
: (%), exp),22
-^ ) Ut(")
: (%)r exp À32
) Ut(r)
I Unç21: (Uo)+exp),az

Où À1, À2, )e etÀa sont les valeurspropres du vecteur A'

On va faire varier Us :

du : duo Az : AU + QQ)B
+ UoAexp
I ffeænAz
Aprés integration de cette expression Us ce met sous la forme suivante:

uo("): B t c
Q|).exp(-A.z)itz
I
D'ori la solution :

f
u"(") : B [email protected](-A(/ - z))dIf c.exp(A.z)
J

( U " r(" ): c.e x' \t.z

U ,zQ): c'e xPÀz' z
=. I c.e æp\z- z
\ u " r(" ):
- )'a'z
lv"n(r)
: I; Q@. exp(-Àa(I z))dI I c'exp
Anneaes 136

Àprés integration par partie d.ela dernière équation, en tenant compte des conditions
anx limites en z:0 (encastrement) et en z:H (bord libre) c.à.d.:

u,(o):w:o
La solution générale de I'expression (2-6) s'écrit sous la forme suivante :

: .*p(#,lo*uft)- B"t"r#)]
u,(,) +# lr,u,
. +]
Les constantes A et B sont donnéespar les expressionssuivantes:

r1
A: * vc1)lexP(-c2H)- 1l
)l ' l 4'11o'+Qct

I r: -A- #1,t,*uc)1t,fi"*ec,H)]
Anneaes 137

Annexe 2

Résultats numériques du deuxième cas de vidange excentrée

paragraPhe 2.4 ChaPitre II

Dans cette annexe nous présentons les résultats numériques donnés par "Abaqus".
Ces résultats représentent les sollicitations liées à des charges dissymétriques. Ces charges
sont
sont adoptées selon les résultats expérimentaux de la base de Chartres [34], elles
représentées sur la figure ci-dessous.
- È!rsi6

."."'PFbtil

. G{oédriæ 2
. ldatièæ silés -BLE'

c
€6

..-----Proæog

t G€ûéralriæ I

. Matièrc æiléc ?LE-

ror30€J0 ZDÉS

Pesios Br.[Én1 PeiætrNbI

-h8i6 -Pr@id h
-.--.-- Froe@d
...----fusæd
'G#EÈi@3
t G&émtriæ4
' MstièF æiléo
t Marike æilée TLE"

-l
t

o
6
66
e6

loÉû?l
tom?o{o5o60æ

kessio [kNrTn'l Pssimt&NÉ41

pressi,onhorizontaleet acti,ond,efrottementselonles résultatsexpéri'mentn,uî

oC"
de la baseile Chartresilans la cellule
Annexes 138

2fi

2æ
2
ld

3 150
E

l0o
€
tIJ

o 5 tocir"oof&*". 20
1.1tt

Fig (2.1) Variation d,eI'effort rnembranairemériil,ien

E00

z
J
rl00
o
a
o
!!) 0
p
o

t -400

400

5 lo
cin"oof&.o." 1r1tt
Fis (2.2) Vari ati on il e I' eff ort m embr anai,re ci'rconf érenci el

no
150

2
J
roo

€so
,ô0

-50

-t 00

-150

o 5 to*o**."r.rtt 20

F;g (2.9) Vari,ati,onile I'effort mentbranaire de cisaillement

Annexes 139

no

P .|00
Ëî
Ëà
0
Ë.8
qË
3E -100

-200

05101520
Circonf&ence [m]

Fig (2.il Variation d,urnomentrnernbranaire

méridien

160

120

O^ 80
€a
EZ
40
69
EE
ËP
ao
6e -40
>'6
40

-120

0510152025
Circonférence [m]

circonférenciel
Fig (2.5) Variation du rùonlent Tùernbr0'n0,i're

30

2A
O^

9Z
Èv t0
Eg
E9 0
Etg
E'6 -10
È€
-æ

-30

Circonfâence [m]

Fis (2.6) Vo,rio,tiond,umoment membro'nairede cisaillement

Annetes 140

Ltt0

100

qd 50
t! ê.
E>
8o
0
oo
.ËF
Eg
ÉE -50

-100

-150
5101520
Circonférence [m]

Fis (2.7) Variati,onde Ia contrainte rnembranaireméridienne

80

60

.9a
€G
Q

5:
Hv 20
oe
ÉI 0
Eg
.EE
Eq o6 -20
o.b
C)o
4

{0

40

Ckconférence [n]

Fig (2.5) Variation d,eIa contrainte membranaire circonférencielle

60

Q
o
- = 6

ËÊi

=à
5v m
E8
o9
'Ë.9 0
50

5€ ao

40

Cilconf&ence [m]
Fig (2.9) Variation d,ela contrainte membranaire de cisaillernent
Anneaes 141

0.001

âs
qÉ
o.ooos
qq

É.:
.9.6
ËE o
ù
>l
{.0005

o.001
to 15
Circonférence
[m]
Fig (2.10) variation de la d,éfonnationmernbranaireméridienne

0.0015

0.001
!)

F-e
€E E 0.0005
qÉO
qo
eE
ËÉ
tre
€
x {.0005

4.00.|

4.d)15 10
o' 5
Circonf&ence [m]

Fi,g(2.11) variation ile Ia d,éfornationrnembranairecirconférencielle

0.ocp8

0.0006
o

E o.oæ+
Er
E
E o.oo.lz
*'â
È3 Ë!
o
È
â {.ooo2

{.0004

4.m06
0 5 Iu
Circonf&ence [m]

Fig (2.12) Variation d,eta d,éformationrnembranaired'ecisaillement

Annexes r42

60

æ
o
0 20
o
:{l

I 0
o

o
o
-æ
I

B -4

60

Circonférence [n]

Fi,g(2.13) Variation du déplacenentcirconférenciel

o
a

Circonférence [n]

Fig (2.1il Variation d,ud,éplacement

rad,ial
-+niva0Sn #niva3m -#uivauTm --Fuivæulljm

A\
€
x

o
o
a
6
I

-l w\\
A,
qV
1-\
Hgt À#N# NIA li \

Y
n"

8' \y\ J l
\.r
V
-2

:l \
tl \l
\}8

Circonf&eoce [n]

Fig (2.15) Variation d,uiléplacementaaial

Annetes L43

150

1q)

z
0
o

-50

-1æ

u*
o'
Circonférence [m]
Fig (2.16) Comparaisondes rnornents0,uniuel,u ile 0.5m

2&

160

880
z
d

fco
-160

-2N

0 51015?o 25
Circonférence [m]
Fis (2.17)Comparaisoniles mornentsau niueau de 7rn

24

160

q80
z
Ë0

j-æ
-ldr

-240

051015âJ25
Circonférence [m]
Fig (2.15) Comparai,soniles mornentsau ni,ueaude 11.5m
Anneaes 744

&)
@

4
d

20
Ë,
o
Ë0
€
6 -20

{0

€0
0510152025
Circonf&encelml
Fig (2.19) Cornparaisoniles contraintesau niaeau d,e0.5m

150

læ

50

-50

-lm

-lto
à 25
Circonférence
[m]
Fig (2.20) Comparaison d,escontraintes au niueau de 7m

læ

â50
=
î
EO

Q -50

-læ

Chconf&ence[m]
Fig (2.21) Cornpariasoniles contrai,ntesau ni,aeauile 11.5m
Anneaes 745

0.ær2

0.ææ

0.0æ4

{.0æ4

{.0æ8

4.0012

Chconférence
[m]
Fig (2.22) Cornparaisondes déform,o,tions0,uniueau ile 0.5m

0.001

I o.m5

ù
)1
H
0

4.0m5

4.æl
0L

Ckconf&ence[m]
Fig (2.29) Cornparaisoniles d,éfonnations0,uniaeau de 7m

0.æl

0.0005

4.0005

4.æ1
oL
Chconférence [m]
au niaeau d,e11.5m
Fig (2.21) Comparaisoniles d,éfonno,tions
Annexes L46
Annexe 3 Présentation de la cellule 66Ct'en vraie grandeure de Chartres.

REPERAGE
DESGÉNÉRATzuCES

10t65

8580
8580

62gl 6Tn

Jl48 51,18
_-- _--
4W M

28@ 28æ

l7l6 l7t6

r
crixÉnhrnrcr crlÉnhrrucez

ld66

8:t80

6Érl

5148
--
M

aû
t7t6
tn

GEI{ERATRICE
3 GÉNÉRATRICE4
Anneaes L47

Annexe 4

Exemple de véritcation de la résistance dtune coque cylindrique

66silot'
raidie verticalement, soumise à des charges de type

Les caractéristiques géométriques ainsi que les conditions d'appuis de la coque cylin-
drique considérées dans cet exemtile sont les suivantes :

Haateur de Ia coque;2Hs= 14 m

Diamètrede Ia cellule; D = I m

Epaissewde la coqae; t = 4 mm

Nornbre de raidisseur ; ns = 16

Sectiondzs raidisseur.i,'A.t = 13.02c#

Rigidité desraidissews,'13= 835cm4

Conditions d'appuis :
Base encastrée-bord supérieur libre.
Cha,rses considérées :
Les charges (pression horizontale et action de frottement de la matière ensilée) agissant
sgr la paroi d.ela coque cylindrique sont donnéespar la norme NF P22-630 et I'Eurocode
L. Elle sont regroupées da,nsle tableau 41.

(*) Les raidisseursconsidéréssont desdemi-poutrellesHEB-100.

( * ) Le toit et d.tautreséquipementannexesne sont pas pris en compte dans cet exemple.
Annetes 148

.hlauteurlm P,(z)-Nl f'(z)-NF P"(z)-EC.L F(z)-EC.1

0.000 22.t52 M.477 25.770 48.924
1.250 20.918 36.585 24.638 40.243
2 .1 8 8 19.892 31.098 23.663 34.208
3.125 L8.773 26.006 22.563 28.606
4.063 17.549 2L.3L2 21.323 23.443
5.000 16.2t4 77.042 19.925 18.746

tableau (4-1) pression'h,orizontale.(fkNlm2l).etsomrnedes actions d,efrottement (kNlml)

ile Ia mati,ère ensi,lée(BIé) selon les nortnes oCas d,eRemplissagecentré".

Conditions de résistance

,r: 1 ozrd.,c.R.r
P
Où " o"r4,c.R.I" est donnéepar P.*Or.r.ron suivante:

( ,æ\ -n n^z A- (t + 4+)

Czrd',C.R.I:Minimumd'e{_r.oas.r".(;)-,.,",+ry*
t \t/ ,\/3(t-"")

F" lMPal b[m Wor lmm ,b 71" lmm r e ff . Ç . R . rl M P a

tr.t2 L .5 7 7.9r 0.94 /.d 48.267
7.78 L .5 7 7.9L 0.94 7.8 48.265
6.50 r.57 7.9L 0.94 7.8 48.257
5.33 I.D { 7.91 0.94 7.8 48.253
4.26 L .5 7 7.9t 0.94 7.8 48.248
3.01 t.5 7 7.9L 0.94 7.8 48.24r

tableau(l-2)

Le tableau (42) montre que les contraintessont largementvérifi.ées.