Formalisation Du Savoir Et Des Outils Dans Le Domaine Des Risques Accidentels

Formalisation du savoir et des outils dans le
domaine des risques accidentels

(DRA-35)
-15
Les clatements de rservoirs
Phnomnologie et modlisation des effets
Ministre de lcologie et du Dveloppement
Durable
Direction des Risques Accidentels

Octobre 2004
ECLATEMENT DE RESERVOIR
INERIS-DRA-2004-46055
Les clatements de rservoirs,

Phnomnologie et modlisation des effets
-15
INERIS
Ce document comporte 72 pages (hors couverture).
Direction des Risques Accidentels
Octobre 2004
1/72
Le prsent document forme un ensemble indissociable. Il ne peut tre utilis que de manire intgrale
Auteurs
Qualit
Emargement
L. HEUDIER
Ingnieur lUnit ExplosionDispersion de la Direction des

Risques Accidentels
Sign
Relecteurs
Qualit
Emargement
Ch. PROUST
Ingnieur lUnit ExplosionDispersion de la Direction des

Risques Accidentels, Animateur
scientifique
Sign
J-C. COUILLET
Responsable de lUnit ExplosionDispersion de la Direction des

Risques Accidentels
Sign
F. MERCIER
Ingnieur lUnit Rsistance des

Structures de la Direction des
Risques Accidentels
Sign
B. DEBRAY
Responsable du Programme
Formalisation du savoir et des
outils dans le domaine des risques
majeurs
Sign
C. JOLY
Responsable du suivi du volet 1 du

programme (ralisation dun recueil
global formalisant lexpertise de
lINERIS)
Sign
Vrificateur final
Qualit
Emargement
O. SALVI
Dlgu Scientifique de la Direction

des Risques Accidentels
Sign
Approbateur
Qualit
Emargement
B. FAUCHER
Directeur des Risques Accidentels
Sign
2/72
PREAMBULE
Le prsent document a t tabli :

-
au vu des donnes scientifiques et techniques disponibles ayant fait lobjet dune

publication reconnue ou dun consensus entre experts,
au vu du cadre lgal, rglementaire ou normatif applicable.
Il sagit de donnes et informations en vigueur la date de ldition du document, le 5 octobre

2004.
Le prsent document comprend des propositions ou recommandations. Il na en aucun cas
pour objectif de se substituer au pouvoir de dcision du ou des gestionnaire(s) du risque ou
dtre partie prenante.
3/72
TABLE DES MATIERES
OBJECTIF ET DOMAINE DAPPLICATION
1.1
Contexte gnral .................................................................................................................. 6
1.2
Objectifs ............................................................................................................................... 6
1.3
Rappel succinct sur le phnomne....................................................................................... 7
1.4
Plan retenu........................................................................................................................... 7
RETOUR DEXPRIENCE
2.1
Incendie sous un rservoir ................................................................................................... 8
2.2
Mise sous pression accidentelle.......................................................................................... 11
2.3
Affaiblissement mcanique ................................................................................................ 12
2.4
Explosion interne ............................................................................................................... 12
2.5
Synthse.............................................................................................................................. 13
DESCRIPTION DU PHNOMNE
3.1
Prsentation gnrale du phnomne ................................................................................ 14
3.2
Production et propagation de londe de pression .............................................................. 15
3.3
Mcanismes de rupture...................................................................................................... 16
3.4
La formation des fragments............................................................................................... 20
3.5
Lacclration et la projection des fragments................................................................... 22
3.6
Les effets sur lenvironnement........................................................................................... 24
QUELQUES MTHODES DE PRDICTION
4.1
Prdire la surpression de ruine ......................................................................................... 28
4.2
Modlisation des consquences : les mthodes globales .................................................... 32
4.3
Modlisation des consquences : les mthodes numriques .............................................. 59
4.4
Les limites des modles existants ....................................................................................... 64
MESURES DE PRVENTION ET DE PROTECTION
14
28
65
4/72
5.1
Rduction du risque la source......................................................................................... 65
5.2
Mesures de prvention....................................................................................................... 65
5.3
Mesures de protection........................................................................................................ 66
CONCLUSIONS
67
GLOSSAIRE
69
BIBLIOGRAPHIE
70
LISTE DES ANNEXES
72
5/72
1 OBJECTIF ET DOMAINE DAPPLICATION

1.1 CONTEXTE GENERAL
Depuis lanne 2000, le Ministre en charge de lEnvironnement (anciennement Ministre de
lAmnagement du Territoire et de lEnvironnement devenu Ministre de lcologie et du
Dveloppement Durable) finance un programme dtudes et de recherches, intitul depuis 2003
Formalisation du savoir et des outils dans le domaine des risques accidentels (DRA-35).
Lobjet du premier volet de ce programme est de raliser un recueil global formalisant
lexpertise de lINERIS dans le domaine des risques accidentels. Ce recueil volutif sera
constitu de diffrents rapports consacrs aux thmes suivants :
les phnomnes physiques impliqus en situation accidentelle (incendie, explosion,
BLEVE),
lanalyse et la matrise des risques,
les aspects mthodologiques pour la ralisation de prestations rglementaires (tude
de dangers, analyse critique..).
Chacun de ces documents reoit un identifiant propre du type -X afin de faciliter le suivi
des diffrentes versions ventuelles du document.
In fine, ces documents dcrivant les mthodes pour lvaluation et la prvention des risques
accidentels, constitueront un recueil des mthodes de travail de lINERIS dans le domaine des
risques accidentels.
1.2 OBJECTIFS
Lobjet du prsent document est de prsenter la dmarche adopte par lINERIS pour ltude
du phnomne dclatement de capacit de stockage de gaz ou de liquide d une
augmentation de la pression dans le rservoir ou une dgradation de ses proprits
mcaniques.
Il sinscrit dans une dmarche de valorisation du savoir-faire de lINERIS auprs des pouvoirs
publics, des industriels et du public.
Lobjectif de ce document est de :
rappeler les phnomnes de mise en pression susceptibles dengendrer lclatement dun

rservoir,
faire le point sur les consquences de lclatement dun rservoir par augmentation de sa
pression interne,
prsenter une synthse de quelques mthodes disponibles pour dcrire les effets engendrs
par lclatement dun rservoir, en comparant des simulations des rsultats dessais et au
retour dexprience daccidents passs,
exposer les limites de ces mthodes.

6/72
1.3 RAPPEL SUCCINCT SUR LE PHENOMENE

Lclatement dun rservoir peut tre la consquence de phnomnes affectant les
caractristiques de lenveloppe (la fatigue de lenveloppe, une corrosion excessive), ou de
phnomnes mettant en cause le contenu du rservoir (une explosion interne ou encore une
augmentation plus lente et accidentelle de la pression sous leffet dun chauffement, dun surremplissage...).
Nous ne dcrirons pas en dtail, dans ce rapport, le phnomne lorigine de laugmentation
de pression. Il convient toutefois de souligner que ce phnomne peut avoir une incidence sur
le mode de rupture du rservoir qui par ailleurs fait lobjet dtudes spcifiques lINERIS.
Les consquences de lclatement dun rservoir traites dans ce rapport sont dune part
lmission dune onde de pression et dautre part la projection de fragments.
Londe de pression rsulte de la dtente brutale du gaz contenu dans le rservoir, ou de la
vapeur si le rservoir contient un liquide surchauff.
1.4 PLAN RETENU

Aprs ce premier chapitre introductif, la suite de ce rapport comporte huit autres chapitres :
-
le deuxime chapitre est ddi au retour dexprience sur des clatements accidentels
ou non,
le troisime chapitre est ddi la description du phnomne,
le quatrime chapitre expose quelques mthodes utilises pour la modlisation de ce

phnomne et leurs limites,
le cinquime chapitre prsente succinctement les mesures de prvention et protection

envisageables,
le sixime chapitre expose les principales conclusions de cette tude.
Dans les trois derniers chapitres sont regroups un glossaire des principales abrviations
utilises, une bibliographie et une liste dannexes.
7/72
2 RETOUR DEXPERIENCE
partir de donnes issues de la base ARIA du BARPI (priode 1943-2003),
(www.aria.environnement.gouv.fr ) et de la littrature (Holden, 1988), un recensement des
accidents en lien avec le thme de ce document a t men et dcompos en quatre grandes
catgories :
les clatements conscutifs un incendie sous un rservoir contenant un liquide

vaporisable,
les clatements conscutifs une mise sous pression accidentelle,
les clatements conscutifs un affaiblissement mcanique,
les clatements induits par une explosion interne.
Le parti pris ci-aprs est de dcrire quelques accidents et essais pour lesquels le niveau
dinformation permettrait une exploitation intressante et de complter cette revue
accidentologique partielle en listant les accidents dcrits en annexe 1.
Toutefois, il est difficile dexploiter les informations issues de linterrogation de la base car,
gnralement, elles ne sont pas exhaustives. Par exemple, la nature du produit contenu nest
pas toujours prcise, de mme que la taille du rservoir. Les distances parcourues par les
fragments sont parfois fournies, mais la masse des dbris nest gnralement pas indique.
2.1 INCENDIE SOUS UN RESERVOIR

Une premire situation pouvant conduire lclatement dun rservoir correspond au cas o
celui-ci est pris dans un incendie.
Un essai chelle relle sur un wagon-citerne de 120 m3 de propane a t ralis White
Sands (Nouveau Mexique, USA) le 28/07/1973 (Holden, 1988). Cet essai informe sur le
processus de rupture et fournit des renseignements chiffrs sur les projections de fragments.
Les informations sur les effets de pression sont inexistantes comme pour les accidents
dailleurs.
Le wagon-citerne, spcialement conu pour lessai, a t construit de faon identique ceux
utiliss pour le transport dammoniac, de GPL, de chlorure de vinyle mise part la prsence
dquipements de mesures. Il a t soumis un incendie ; la monte en temprature a entran
une augmentation de la pression interne du rservoir.
Le flux de chaleur moyen reu sur la surface de la citerne a t valu exprimentalement
105 kW/m. La soupape de scurit sest souleve aprs 2,2 min lorsque la pression de 19 bars
a t atteinte. Elle a effectu trois cycles ouverture-fermeture avant de rester ouverte.
La pression interne avait atteint 24 bars au bout de 24,5 min lorsque le rservoir a clat. La
pression et la temprature internes au moment de lclatement taient infrieures celles
8/72
dclenchant les dispositifs de scurit ; ceci peut sexpliquer par laugmentation de temprature
du matriau qui entrane une diminution de son seuil de rupture.
127 fragments provenant du wagon-citerne et des quipements ncessaires aux diffrentes
mesures ont t retrouvs, leurs positions ont t rpertories. Parmi ces fragments, 63
provenaient du wagon et une dizaine provenait de la citerne, les distances et directions de ces
derniers ont t rpertories. La plus grande distance parcourue tait de 407 m.
tant donn que le wagon tait plac dans une fosse de 8 m de profondeur, les fragments qui
sont retombs lextrieur de cette fosse ont peut-tre parcouru des distances plus faibles que
celles que lon aurait observes si lessai avait eu lieu au niveau du sol.
Les thermocouples indiquaient un gradient de temprature considrable le long du toit de la
citerne, au niveau de la rupture. Au point damorce de la fissure, sur lune des extrmits de la
citerne, une temprature de 640C a t releve alors que sur lautre extrmit, elle slevait
seulement 450C. Le point le plus chaud correspond au point probable damorce de la
fracture. Cette observation est cohrente avec les connaissances de la rsistance des matriaux
qui prdisent pour les aciers une diminution drastique de la rsistance ultime partir de 300C.
Le gradient de temprature sur la citerne au moment de la rupture est assez important et
confirme que la rpartition de temprature lorsquun rservoir est pris dans un incendie peut
tre trs htrogne et que ce sont les points les plus chauds quil faut prendre en
considration.
Des accidents correspondant ce phnomne et dont les principales donnes sont reprises dans
le Tableau 1, sont dcrits dans lannexe 1.
Date
Lieu
Description du
rservoir
Nature du
contenu
Effets de
lclatement
Avril 1970
Raffinerie
(Canada ou
USA)
Epurateur de
gaz caustique
Projection des
deux extrmits
du rservoir,
lune plus de
17m
9/72
Aot 1972
Raffinerie
(Canada ou
USA)
Rservoir de
3,5 m de
diamtre et 10 m
de long
Rservoir de
1 m de diamtre
et 5 m de long
clatement 10
15 minutes aprs
le premier
clatement ;
projection de la
plus grosse
portion du
rservoir
08/04/1979
Floride
Wagon-citerne
Ammoniac
clatement de la
citerne en deux
sections
projetes
200 m et 75 m
19/11/1984
Mexique
Sphres
(2400 m3 et
1600 m3 )
Butane
Fragments
projets 600m
Cylindres (de 36
270 m3 )
Butane
Projection de 12
cigares-rockets ;
lun (provenant
dun cylindre de
36 m3 ) a t
retrouv
1200 m
22/07/1990
Core du Sud
Deux rservoirs
Butane
20/08/1991
France,
Charente
Bouteilles de
propane et
butane, camionciterne
Propane et
butane
26/06/1994
Venezuela
Deux rservoirs
(propane et
propylne)
10/72
10/03/1996
Italie
Citerne routire
de 12 m3
Propane
Projection dun
fragment
500 m
Tableau 1 Principales donnes recenses pour des accidents dclatements de rservoirs pris dans des
incendies
2.2 MISE SOUS PRESSION ACCIDENTELLE

La mise sous pression accidentelle dun rservoir peut tre due par exemple la perte de
rgulation de dbit ou une injection accidentelle de produit dans le rservoir.
Le 14/04/1983, en Indonsie, une colonne dchangeur de chaleur de GNL de 47 m et de 5 m
de diamtre maximal a t purge avec un hydrocarbure gazeux chaud et sec afin de la dgivrer
et de la dshydrater. La pression maximale de service de la colonne tait de 2 bars et la
soupape tait tare 4 bars. La source du gaz de purge tait plus de 35 bars. Le
dysfonctionnement dune soupape a entran la pressurisation de la colonne. Une pression de
5 bars a t enregistre ; 20 minutes plus tard, la colonne sest rompue. Le rservoir sest
fragment en trois morceaux principaux, lun a heurt une construction 50 m.
Il convient de remarquer que la pression laquelle la rupture a t observe parat compatible
avec les rgles classiques de conception qui fixent la pression de service environ un tiers de la
pression de rupture du rservoir.
Des accidents correspondant ce phnomne et dont les principales donnes sont reprises dans
le Tableau 2, sont dcrits dans lannexe 1.
Date
Lieu
Description du
rservoir
Nature du
contenu
Effets de
lclatement
Janvier 1966
Raffinerie
(Canada ou
USA)
Rservoir
Projection dun
quipement
28 m
Septembre 1970
Raffinerie
(Canada ou
USA)
Racteur
Azote (essai de Projection dune

mise en pression fraction de la
avec de lazote tte du racteur
27 bars)
81 m
Tableau 2 Donnes principales recenses pour des clatements conscutifs une mise en pression
accidentelle
11/72
2.3 AFFAIBLISSEMENT MECANIQUE

Un affaiblissement mcanique de lenveloppe dun rservoir sous pression peut en provoquer la
rupture. La recherche bibliographique relative au retour dexprience, ralise dans le cadre de
ce document, na pas permis didentifier daccidents de ce type, en dehors de ceux dus une
lvation de temprature sous leffet dun incendie, traits au paragraphe 2.1. Pour autant,
cette catgorie dclatement sera tudie par la suite.
2.4 EXPLOSION INTERNE

Lexplosion conduit souvent la rupture de la capacit dans laquelle elle se produit. Une
explosion est une libration brutale dnergie ; diffrents phnomnes peuvent tre lorigine
dune explosion interne : une combustion, un emballement de raction chimique, un
changement de phase brutal
Un exemple dclatement dune capacit de stockage suite une explosion interne est fourni
par laccident survenu le 20/08/1997 Blaye. Une violente explosion de poussires a largement
dtruit une installation de stockage de crales.
Ce silo, de type vertical, se composait de trois ranges de cellules (44 au total) de section
circulaire (6,2 m de diamtre et 33 m de hauteur) en bton arm (15 cm dpaisseur). Une dalle
en bton arm les recouvrait et servait de plancher une galerie arienne. Les produits
stocker taient introduits par des orifices de 60 cm sur 60 cm pratiqus dans ce plancher. La
galerie de manutention abritait des convoyeurs bandes. Elle communiquait chaque
extrmit avec des tours en bton arm de 50 m de hauteur.
On estime que les fts des cellules taient capables de supporter une pression interne de lordre
de 1 bar. Pourtant, lexplosion a compltement dtruit les cellules de la partie centrale, les
tours, la galerie de manutention et les cellules attenantes lune des deux tours. De gros
morceaux (dimension mtrique ou plus) ont t retrouvs jusqu une distance de lordre de
50 m et des morceaux plus petits (de lordre de 10 cm) jusqu 150 m. Des vitrages ont t
briss jusqu une distance de 500 m du lieu de lexplosion.
Les simulations de cet accident avec le code EFFEX (Proust, 2000) suggrent que la
surpression dans les cellules a t de lordre de 3 4 bars, soit nettement au-del de la pression
standard de ruine des fts. Ceci sexplique essentiellement par le dveloppement extrmement
brutal de lexplosion, dans un laps de temps bien infrieur celui de la rupture et de la
dcharge de la pression par les brches. Cependant, il ne faut pas exclure non plus un effet
direct de la sollicitation mcanique sur le processus de rupture qui pourrait se traduire par une
lvation ou une rduction de la pression de ruine par rapport la situation statique, comme
nous le verrons aux paragraphes 3.3.3 et 4.1.3.
Dautres accidents correspondant ce phnomne dont les principales donnes sont reprises
dans le Tableau 3 ont t recenss et sont dcrits en annexe 1.
12/72
Date
Lieu
Description du
rservoir
1981
Royaume Uni
Schoir
Novembre 1962
Raffinerie
(Canada ou
USA)
Accumulateur
Nature du
contenu
Effets de
lclatement
Hydrosulfate de Propulsion de la
soude, prsence
plus grosse
deau dans le
moiti du
fond
rservoir
travers lusine
Prsence deau
dans le fond
Projection de la
partie principale
93 m,
projection dun
morceau de 3 m
sur 10 m
220 m
Tableau 3 Donnes principales recenses pour des clatements conscutifs une explosion interne
2.5 SYNTHESE
Le retour dexprience effectu sur les accidents de la base ARIA permet de tirer plusieurs
enseignements sur le phnomne dclatement. Parmi les principaux, nous pouvons retenir
linfluence de la temprature et du mode de sollicitation (vitesse daugmentation de la pression)
sur le seuil de rupture dune enceinte.
Ces influences seront dtailles dans le chapitre suivant qui est consacr la description du
phnomne dclatement de rservoir.
Concernant les effets sur lenvironnement, les informations sur les effets de pression ne sont
pas fournies dans les descriptions des accidents ; en revanche, les distances de projection de
fragments sont parfois indiques et la distance maximale atteinte dans les accidents mentionns
prcdemment est de 1200m (rservoir qui a t projet par effet fuse ).
13/72
3 DESCRIPTION DU PHENOMENE
3.1 PRESENTATION GENERALE DU PHENOMENE
Lors de laugmentation de la pression dans un rservoir, les parois les plus fragiles se rompent
lorsque la contrainte ultime (de ruine) est atteinte. La rupture du confinement permet la
libration de la pression contenue qui se traduit par la propagation externe dune onde de
pression arienne (dplacement dune surpression dans lair). Pour utiliser une analogie simple,
le rservoir peut tre assimil un haut-parleur. Le son dlivr par ce haut-parleur est
directement proportionnel lnergie quon lui soumet. Limpression ressentie par un
observateur sera fonction, outre de cette nergie, de la distance et de la taille du haut-parleur.
On comprend ainsi que lamplitude de londe arienne produite par un clatement pneumatique
doit tre fonction de la distance, de la taille du rservoir et de la surpression interne maximale.
Nous verrons que cette amplitude dpend en fait du paramtre (PruptV)1/3 , o Prupt est la
pression de rupture du rservoir et V son volume de gaz.
Pour mmoire, la propagation dune onde de pression arienne dans lenvironnement se traduit
par des variations plus ou moins brusques et de plus ou moins grande amplitude de la pression
en tout point de lespace. En un point donn, ces variations de pression sont notamment
caractrises par :
une phase de surpression, de dure t+ et damplitude maximale note P+,
une phase de dpression, de dure t- et damplitude maximale note P-.

P
P+
Temps
P t+
t-
Figure 1 Profil dune onde de pression
La rupture saccompagne dun clatement du rservoir et de la formation ventuelle de

plusieurs fragments qui peuvent tre propulss sous leffet du diffrentiel de pression entre
lintrieur et lextrieur.
Les fragments suivent une trajectoire dtermine, outre par les conditions initiales, par les
effets de freinage arodynamique, de portance ventuelle et de pesanteur. Notons que
limpulsion est communique au fragment par les jets de fluide issus des brches et que, par
14/72
consquent, leffort de pousse doit se rduire mesure que chaque fragment sloigne de son
point de dpart et dvie de laxe du jet qui le propulse sous leffet de la gravit terrestre.
3.2 PRODUCTION ET PROPAGATION DE LONDE DE PRESSION

Pour valuer la production dune onde de pression, on suppose que les parois disparaissent
instantanment au moment de la rupture. On observe alors une mise en vitesse brutale de
latmosphre environnante opre par une onde de choc qui se propage dans latmosphre.
Dans ces premiers instants, ce phnomne peut tre dcrit qualitativement au moyen de la
thorie des tubes choc (Wright, 1961). Un tube choc est un long rservoir cylindrique
spar en deux compartiments par une cloison perpendiculaire laxe du rservoir. Lun des
compartiments est rempli dun gaz haute pression, lautre dun gaz basse pression.
Chambre haute pression
percuteur
Chambre basse pression
Membrane sparatrice
Figure 2 Tube choc
La Figure 2 reprsente le dispositif. En simplifiant lextrme, le fonctionnement peut tre

dcrit de la faon suivante :
Au moment o le percuteur crve la membrane sparatrice, une onde de choc1

damplitude intermdiaire entre la pression haute et la pression basse part de la
membrane et se propage travers la chambre basse pression avec une clrit au moins
gale celle du son dans la chambre basse pression.
Une onde de dtente se propage travers la chambre haute pression, se rflchit sur le
fond et poursuit ensuite londe de choc. Au bout de trois trajets de londe de dtente, la
Une onde de choc est une onde de pression arienne pour laquelle la monte en pression est trs brutale :
P
P+
Temps
15/72
pression a diminu au point que le fond de la chambre se retrouve la pression

atmosphrique (Baum, 2001).
Les caractristiques des ondes peuvent tre calcules en appliquant la conservation de la masse
et de la quantit de mouvement. La pression de londe de choc est de quelques diximes de la
pression du rservoir au moment de la rupture.
Lors dun accident rel, au moment de lclatement, les parois (qui constituent lenveloppe du
rservoir) ne disparaissent pas instantanment, elles sacclrent progressivement et
ralentissent donc le fluide initialement contenu dans le rservoir du fait de leur inertie. La
dpressurisation est donc plus lente, moins brutale, et londe de choc moins forte. lheure
actuelle, on ne sait pas facilement tenir compte de ce phnomne, on considre donc souvent
des parois sans inertie pour estimer lamplitude des ondes (par excs donc) au moyen par
exemple de la thorie des tubes choc.
Les ondes de dtente finissent par rattraper londe de choc pour former londe de pression
arienne classique dont lamplitude, en vertu des lois de la mcanique des fluides, doit
sattnuer en fonction de la distance au centre lexplosion.
3.3 MECANISMES DE RUPTURE

3.3.1 MODES DE RUPTURE
lchelle macroscopique, mesure que la pression interne augmente, le matriau subit une
dformation lastique rversible tant que les contraintes restent en-de dun seuil critique (la
limite dlasticit) puis, au-del, une dformation plastique permanente peut se mettre en place
avant lapparition de la ruine.
lchelle atomique, les dformations lastiques et plastiques respectent la cohsion de la
matire sachant que le matriau est parsem naturellement de micro-fissures (de lordre du
micron). La rupture survient lorsque cette cohsion est dtruite. Elle sopre par
agrandissement des discontinuits au sein de la matire ; les micro-fissures deviennent des
macro-fissures en atteignant des dimensions de lordre du millimtre, et enfin des fissures
lchelle des structures mcaniques (Mercier, 2003).
La rupture fragile
Si la rupture se produit alors que les dformations sont lastiques (rversibles), le
comportement du matriau et la rupture sont dits "fragiles". Le verre, la fonte grise, les aciers
bruts de trempe, les cramiques, le bton, la plupart des polymres thermodurcissables et
certains alliages haute rsistance, sont des matriaux qui ont un comportement fragile.
La rupture fragile correspond une situation o la propagation des fissures est difficile. Cest
le cas par exemple des alliages mtalliques dans lesquels ont t augmentes fortement les
caractristiques de rsistance. Ainsi, trs gnralement, la fragilit augmente avec la limite
16/72
dlasticit. La rupture fragile se produit lorsque les liaisons inter-atomiques se rompent sans
dformation plastique globale (Mercier, 2003).
Comme le phnomne de rupture est relativement local, un nombre important de petits
fragments peut tre form suite une multiplication des bifurcations de fissures. Il peut tre
admis, par simplification, que la rupture fragile peut tre quasi-instantane.
La rupture ductile
Si la rupture se produit aprs une dformation plastique permanente, le comportement du
matriau et la rupture sont dits "ductiles". La majorit des mtaux, des alliages et certains
polymres thermoplastiques prsentent ce type de comportement.
La rupture ductile rsulte de lagrandissement des micro-fissures et la propagation de la fissure
rsultante.
La rupture ductile, gnratrice de fissures importantes, forme peu de fragments, mais
gnralement de grande taille (Mercier, 2003). La propagation du phnomne est limite par la
vitesse du son dans le matriau.
Transition fragile-ductile
Plusieurs facteurs influencent la fragilit ou la ductilit des mtaux :
La temprature : une temprature assez basse, on observe que la dformation

plastique devient impossible avant rupture du matriau. Ainsi, il existe une temprature
dite de transition au-dessous de laquelle le matriau devient fragile. Au fur et mesure
de laugmentation de temprature, la transition entre la rupture fragile et la rupture
ductile se fait ; il existe un domaine de temprature o la rupture est mixte.
La vitesse de dformation : en gnral, la ductilit diminue si la vitesse de dformation

augmente.
Leffet dentaille : la prsence dun dfaut ou dune entaille entrane une concentration
de contraintes qui peut modifier le rgime de rupture.
La rupture lors de la transition fragile-ductile gnre en gnral peu de fragments de taille

importante.
3.3.2 REGLES DE CONCEPTION DES EQUIPEMENTS MECANIQUES

Le dcret du 13 dcembre 1999 constitue la transposition en droit franais de la directive
europenne de 1997 (97-23-CE) relative aux quipements sous Pression. Le champ
dapplication de ce dcret concerne la conception, la fabrication, lvaluation de la conformit
et la mise sur le march des quipements sous Pression. Larrt du 15 mars 2000 concerne lui
le suivi en service de ces quipements, qui doit tre assur par lexploitant.
17/72
En matire dclatement, les contraintes rglementaires sont essentiellement fonction du

produit PS*V o PS est la pression maximale admissible et V le volume du rservoir considr.
Ce choix parat assez pertinent dans la mesure o ce paramtre reprsente lnergie interne
stocke sous forme de pression.
La rglementation impose pour la construction des appareils pression de gaz que le mtal
utilis soit exempt de fragilits dans le domaine de tempratures et de pressions auquel seront
normalement soumis les quipements. On dfinit la pression de calcul, fixe par le
constructeur, qui est au minimum gale la pression extrme dutilisation (pression maximale
de service).
Il est prvu par la rglementation que la contrainte maximale dans le mtal la pression de
calcul doit tre plus petite que les deux valeurs suivantes :
R
, o R est la rsistance ultime ( la rupture) conventionnelle la traction du
3
mtal,
Re T
0 , 002
, o ReT
0 , 002
est la limite dlasticit conventionnelle, ou contrainte 0,2%

1,6
de dformation du mtal.
En pratique, les rservoirs mcano-souds sont souvent conus pour pouvoir supporter, en
rgime statique, une pression au moins gale trois fois leur pression de calcul.
3.3.3 INCIDENCE DU MODE DE SOLLICITATION

Incidence de la fatigue :
Sous l'action de charges rptes ou alternes, la rupture d'un matriau peut se produire sous
des contraintes infrieures la rsistance du matriau en statique : c'est la rupture par fatigue.
Ce mcanisme de rupture proviendrait dune augmentation de la fissuration lors des
chargements rpts.
Lessai de fatigue le plus simple consiste soumettre des prouvettes des cycles defforts
priodiques, damplitude maximale et frquence constantes, et noter le nombre de cycles N
au bout duquel la rupture se produit. Ce nombre N est report, en gnral sur une chelle
logarithmique, en fonction de la contrainte maximale des cycles. chaque prouvette
correspond un point du plan (, N) et, partir dun lot dprouvettes soumises des
contraintes maximales diffrentes, on obtient une courbe qui a lallure de celle reprsente sur
la Figure 3, connue sous les noms de courbe de Whler et courbe S-N (Stress-Number of
cycles). Dans de nombreux cas, il est possible de tracer une branche asymptotique la courbe
de Whler, lasymptote tant la limite dendurance ou limite de fatigue D (Rabbe et al., 2000).
18/72
Re
D
105
106
107
108
Figure 3 Evolution de la contrainte en fonction du nombre de cycles N (Re est la limite dlasticit et D
la limite de fatigue)
Des expriences menes sur des prouvettes lisses (Rabbe et al., 2000) ont permis dtablir une
relation entre la limite de rsistance la fatigue D et la rsistance la traction. Il en ressort
que la limite de rsistance la fatigue peut tre comprise, grossirement, entre un tiers et deux
tiers de la rsistance la traction2.
Si lon extrapole cette relation au cas des rservoirs soumis des cycles de charge, la pression
de rupture dun rservoir dont lenveloppe serait fatigue peut tre trois fois plus petite que sa
pression de rupture en rgime statique.
Incidence de la temprature :
On rappelle que les caractristiques mcaniques des matriaux dpendent de la temprature. La
figure ci-dessous trace lvolution en fonction de la temprature de la contrainte limite
dlasticit, de la contrainte limite de rupture, du module dlasticit et du coefficient de
dilatation linaire pour les aciers. Ces courbes ont t obtenues partir dessais de traction sur
des prouvettes.
Pour les essais de traction sur des prouvettes, la rsistance la traction peut galement tre appele
contrainte la rupture ou charge de rupture.
19/72
Figure 4 Caractristiques mcaniques de lacier en fonction de la temprature, sur prouvette de traction

(Fruitet, 1979)
La contrainte de rupture des aciers par exemple diminue au-del dune temprature de
quelques centaines de degrs Celsius, si bien quun rservoir pris dans un incendie voit sa
pression interne augmenter et sa contrainte de rupture diminuer. Linfluence de la nature de
lacier tant faible, il est possible de considrer en premire approximation, que les courbes de
la Figure 4 sont valables pour tous les aciers de construction, tant entendu que les ordonnes
sont lues en valeur relative par rapport aux caractristiques de base temprature ordinaire
(Fruitet, 1979).
Incidence de la vitesse daugmentation de pression :

Baker (1983) par exemple, fait remarquer que la rsistance dun quipement dpend de la
vitesse de la sollicitation mcanique. Ce phnomne est pris en compte en ajoutant dans les
matrices de dformation de lquipement considr les termes dinertie. Globalement, le
systme se rsume dans ses configurations les plus simples, un modle masse-ressort o la
raideur reprsente llasticit de la structure et la masse son inertie. Selon le type de
sollicitation, lacclration de la masse peut sajouter (se retrancher) leffort externe, et
augmenter (diminuer) les dformations, si bien que la rupture apparat pour des pressions plus
basses (plus leves) quen rgime statique (Proust, 1997).
3.4 LA FORMATION DES FRAGMENTS

Les diffrentes techniques de dimensionnement des structures (ex : EUROCODES) permettent
destimer la surpression de ruine dune paroi lorsque la surpression est applique de manire
20/72
lente, typiquement lorsque la vitesse de chargement est plus petite que 0,01 bar/s et lorsque la
temprature du matriau est homogne. En rgime de chargement dynamique, tant que la
dure caractristique du chargement reste comparable (ou plus grande) en ordre de grandeur
celle de la dformation de la structure considre, on peut admettre que les fragments
susceptibles dtre forms devraient avoir des dimensions identiques celles des fragments
forms lors de la rupture sous chargement statique. Dans la situation inverse, les rgimes de
vibration sont susceptibles de devenir plus complexes et la taille des fragments ne peut plus
tre prdite simplement (Proust, 2000).
La Figure 5 illustre un exemple extrme de dformation dun rservoir cylindrique contenant
un gaz sous pression initie par une fissure axiale. Il convient de souligner que chacune des
gomtries de cette figure pourrait tre une gomtrie finale. Ces diffrentes gomtries
reprsentent donc les gomtries des missiles potentiels.
La fissure initiale se propage le long de laxe vers chaque extrmit du rservoir o la fracture
bifurque. Ensuite, des fractures circonfrentielles tournent dans des directions opposes autour
de chacune des extrmits. Un grand volet se forme lorsque la fracture progresse autour des
extrmits. Les faces des volets, qui lorigine taient les surfaces extrieures du rservoir,
convergent lune vers lautre. Les fractures circonfrentielles se propagent et se rencontrent
ventuellement entranant le dtachement dune ou deux extrmit(s) du rservoir.
Figure 5 Les dformations successives du rservoir (Baum, 2001)
21/72
Lorsque se superposent des gradients thermiques induits par un feu, la rpartition des
contraintes peut tre htrogne et il parat plus difficile de prdire le nombre et la forme des
fragments qui seront produits lors de lclatement. Le calcul en statique homogne ne peut
fournir quun ordre de grandeur.
3.5 LACCELERATION ET LA PROJECTION DES FRAGMENTS

3.5.1 REPRESENTATION PHYSIQUE
La vitesse et la distance parcourue par les fragments sont principalement fonction de leur
vitesse initiale, de leur masse, et de leur gomtrie.
On peut montrer partir des thories de la mcanique que la vitesse initiale du fragment
dpend essentiellement de limpulsion I ( I p(t )dt o p(t) est la surpression dans le
rservoir) communique lors de la dpressurisation du rservoir aprs clatement. p(t) est une
fonction dcroissante si bien que limpulsion est une grandeur finie. Plus elle est grande, plus la
vitesse initiale est importante, toutes choses gales par ailleurs.
Aprs la phase de pousse, le fragment est soumis aux seules forces de trane et de pesanteur,
plus rarement des effets de portance. En effet, la portance pourrait intervenir si les profils
extrados (surface suprieure du fragment) et intrados (surface infrieure du fragment) taient
dissymtriques.
La trajectoire peut tre dduite des quations de la balistique (Figure 6), connaissant
notamment langle de projection, la hauteur de dpart H0 et la vitesse initiale Vi.
a) Trajectoire du fragment
22/72
b) Bilan des forces ( F1 force de portance, Fd force de trane. P force de pesanteur, V vitesse du
fragment)
Figure 6 Trajectoire du fragment et bilan des forces extrieures exerces sur le fragment
Bien videmment, les fragments ne peuvent pas avoir une vitesse initiale plus grande que la
vitesse matrielle des gaz qui les propulsent. Ce paramtre peut tre extrait de la thorie des
tubes choc. Il correspond la vitesse limite dun fragment de masse nulle.
3.5.2 LORSQUE LE RESERVOIR CONTIENT DES GAZ SOUS PRESSION

Baum a tudi la formation et la mise en vitesse des fragments rsultant dune rupture axiale
dun rservoir cylindrique de gaz sous pression (Baum, 2001). Il a dvelopp un modle simple
permettant de prdire la vitesse atteinte par un rservoir non retenu lorsque sa rupture a t
initie par une fissure axiale. La vitesse prdite par ce modle a ensuite t compare aux
rsultats dune tude exprimentale. La description des essais raliss par Baum ainsi que la
comparaison avec les rsultats thoriques sont dtaills en annexe 2.
Lorsque le gaz haute pression initialement contenu dans le rservoir est expuls par la brche
en cours de dveloppement, le rservoir subit une pousse gale la raction de lexpulsion du
gaz.
Baum a observ que le processus dexpulsion de gaz est entirement termin lorsque la
dformation du rservoir a atteint la forme de la Figure 5b. Ainsi, on peut considrer que
lacclration est termine alors que le rservoir est encore en un seul morceau.
23/72
3.5.3 RESERVOIR CONTENANT UN LIQUIDE

Baum (1999) a tudi la rupture circonfrentielle dun cylindre contenant un liquide surchauff.
Les morceaux du rservoir subissent une acclration semblable celle dune fuse lorsque le
fluide est expuls.
La dpressurisation soudaine induite par la rupture du rservoir initie la vaporisation brutale du
liquide. Les bulles de vapeur entranent laugmentation rapide du volume du fluide. Ce
processus de "flash" maintient la pression lintrieur du rservoir une valeur proche de la
pression de saturation du liquide.
Compte tenu de ltape supplmentaire de vaporisation, lexpansion et lexpulsion du mlange
vapeur-liquide sont plus longues et une plus grande impulsion est communique aux fragments
par rapport la situation prcdente (3.5.2).
3.6 LES EFFETS SUR LENVIRONNEMENT

Lors dun clatement de rservoir, dun point de vue physique, il est cohrent de considrer
deux sortes deffets sur lenvironnement extrieur : les effets de surpression dus londe de
pression, et les effets dimpact des fragments projets.
3.6.1 LES EFFETS DE SURPRESSION

La pression est une force par unit de surface susceptible dinduire des efforts de flexion ou de
cisaillement dans les structures, ventuellement de compression pour le corps humain. Une
onde de pression peut galement propulser des projectiles.
Des seuils deffets sur lhomme sont proposs dans la littrature et, lheure actuelle, les seuils
utiliss dans les tudes rglementaires sont les suivants :
le seuil 14 kPa qui exprime l'apparition d'effets ltaux,
le seuil 5 kPa qui exprime l'apparition des effets irrversibles.
Les mmes seuils sont mentionns dans larrt du 26 septembre 1980, scurit
pyrotechnique (JORF 1988, n1196) qui fait encore autorit aujourd'hui dans son domaine.
Les donnes qui ont servi de base l'tablissement de cette rglementation sont issues de
nombreux essais d'origine militaire, sur des munitions notamment. L'exprience que nous en
avons est que les seuils et valeurs proposs dans cette rglementation en termes de dgts aux
biens sont trs cohrents et proches de l'observation exprimentale (10 20% de marge).
La rglementation pyrotechnique prsente des tableaux qui permettent d'estimer les effets une
distance donne en fonction de la charge d'explosif dtonant. Les abaques des explosifs
24/72
permettent de relier cette charge un niveau de surpression avec une trs bonne prcision
(Baker, Lannoy) et on obtient la lecture suivante :
Niveau de surpression
P > 40 kPa
40 kPa> P > 17 kPa
17 kPa> P > 7 kPa
7 kPa> P > 5 kPa
5 kPa> P
Dommages aux
Blessures mortelles dans
Blessures graves,
Blessures
Possibilits de blessures Trs faibles possibilits
personnes
plus de la moiti des cas ventuellement mortelles
de blessures lgres
Dgts aux biens
Dgts trs graves
Dgts importants
Dgts moyens et lgers
Dgts lgers
Dgts trs lgers
Tableau 4 Dommages aux personnes et dgts aux biens selon les niveaux de surpression
On constate en premier lieu, que le seuil d'apparition des effets irrversibles apparat assez
nettement 5 kPa tandis que le seuil d'apparition des effets ltaux peut tre situ vers 17 kPa.
En second lieu, on observe une association assez troite entre les dgts aux biens et les
dommages aux personnes, ce qui parat finalement assez logique puisque la probabilit est
grande pour que les personnes concernes par une explosion soient au moment de l'accident
l'intrieur d'un btiment.
D'ailleurs, si lon s'intresse aux effets directs d'une onde de surpression sur l'tre humain, on
trouve typiquement que le seuil d'apparition de blessures (rversibles) aux tympans est de plus
de 15 kPa tandis que le seuil d'apparition de la ltalit (hmorragies pulmonaires) serait situ
vers 100 kPa (Green Book, TNO).
On doit donc considrer que les seuils fournis dans la rglementation pyrotechnique sont des
effets indirects correspondant des dgts aux structures environnantes.
On peut interprter lapproche pyrotechnique de la manire suivante. partir de 5 kPa, des
lments lgers et non porteurs des structures peuvent tre dtruits ou endommags,
susceptibles d'occasionner des blessures peu graves. partir de 17 kPa, les infrastructures des
locaux sont fortement endommages et un effondrement partiel est possible ( dgts
importants ) sur les occupants, ce qui engendre le risque ltal.
D'un point de vue pratique, ces deux seuils sont employs plutt pour valuer les effets sur
l'homme et il ne semble pas exister d'approche bien tablie pour ce qui concerne les atteintes
aux structures du point de vue de la conservation de leur fonctionnalit. On peut cependant
penser que le seuil d'apparition des dgts graves aux btiments (17 kPa) de la rglementation
pyrotechnique pourrait galement tre considr comme un seuil pour des dgts graves aux
quipements les rendant impropres l'usage voire dangereux (rupture de rservoirs
d'hydrocarbures et destruction de fondations partir de 20 kPa d'aprs Clancey).
Bien que les donnes qui ont permis d'tablir ces seuils paraissent assez bien tablies et
raisonnables, il est utile de rappeler qu'elles sont associes la dtonation d'explosifs
condenss. Il n'est pas vident que ces seuils soient conservs pour la situation plus classique
des dflagrations et peu d'informations sur ce point sont disponibles puisque les analyses
d'accidents mettant en scne des dflagrations sont gnralement faites sur la base d'un
quivalent TNT .
25/72
On peut considrer que le risque de blessures est susceptible de se matrialiser lorsque les
individus sont frapps par des fragments de vitres, de bois, des objets de dcoration lgers. Ce
risque ne devrait pas tre fatal tant que les structures plus lourdes comme les murs porteurs ne
sont pas atteintes. Le niveau de surpression correspondant pour la dtonation dun explosif
condens est de lordre de 5 kPa (Clancey).
Ds que le risque deffondrement apparat, le risque ltal est prsent, par effet dcrasement ou
de chocs de fragments massifs. On peut estimer que les dommages aux structures peuvent alors
devenir suffisamment importants pour que les moyens de production industriels soient
lourdement affects, notamment les moyens de contrle. On peut situer le dbut de ce risque
vers 14 kPa (Clancey, murs de briques, bton non renforc) pour la dtonation dun explosif
condens.
Dommages aux
personnes
Apparition du risque ltal Apparition du risque de

par crasement
blessure par fragments
lgers
Dgts aux biens
Apparition de dgts Destruction des lments
majeurs aux moyens de lgers et non porteur des
production et aux
btiments
btiments par
effondrement
Niveau de surpression
14 kPa
5 kPa
Tableau 5 Dommages aux personnes et dgts aux biens pour les niveaux de surpression de 50 et
140 mbar
3.6.2 LES EFFETS DIMPACT DES FRAGMENTS

Lors de lclatement dun rservoir, une distinction est faite entre :
les missiles primaires qui proviennent de la source elle-mme, il sagit des fragments
du rservoir ; ils sont dtermins par leur nombre, leur masse, leur forme, leur vitesse et
leur trajectoire,
et les missiles secondaires : il sagit dobjets dans lenvironnement du rservoir qui

sont soulevs par londe de pression de lclatement. Leffet de ces missiles est moindre
car ils possdent beaucoup moins de quantit de mouvement. Cependant, ils doivent
tre pris en compte en cas dabsence de fragments primaires.
Limpact dun projectile est susceptible de produire deux types deffets sur lhomme ou sur des
structures :
-
le choc qui peut tre la cause de fractures sur le corps humain, et entraner la
dformation des structures, ventuellement jusqu leur effondrement et leur
destruction ;
la pntration, que ce soit dans le corps humain ou dans un quipement cible.
Le choc dun projectile anim dune quantit de mouvement suffisante peut entraner la
pntration dans la cible et conduire la ruine des quipements.
26/72
La rglementation pyrotechnique (arrt du 26 septembre 1980 fixant les rgles de

dtermination des distances d'isolement relatives aux installations pyrotechniques) fixe deux
seuils lis aux effets de pntration de fragments de petite taille (moins de 1 kg), exprims en
nergie cintique du projectile :
20 J la limite Z2 / Z33 (blessures graves pouvant tre mortelles) est le seuil de

ltalit ;
8 J la limite Z3 / Z4 (blessures) est le seuil des blessures significatives.
Pour des projectiles massifs (masse suprieure 1 kg), il ne semble pas pertinent de raisonner
en terme dnergie cintique, mais plutt en terme de vitesse dimpact. Lexprience de
lINERIS montre que la vitesse dimpact dun projectile au sol est gale sa vitesse de chute
libre, et est souvent de plusieurs dizaines de m/s. Les donnes disponibles (Baker, 1983 ; TNO,
1989) suggrent de retenir, pour caractriser les effets dun projectile de masse suprieure 1
kg sur lhomme, un seuil de 4 m/s pour lapparition des blessures irrversibles (fracture du
crne).
Pour ce qui est des quipements, il nexiste pas, lheure actuelle, de critre simple qui
permette de caractriser leffet dun projectile sur une structure.
La trs grande varit des quipements cibles rend ncessaire une tude au cas par cas. On
conoit en effet aisment quune maison dhabitation, par exemple, possde une rsistance
mcanique trs diffrente dun quipement industriel conu pour rsister de fortes pressions.
Classification rglementaire des zones de dangers pour une installation pyrotechnique :
Dsignation de
la zone
Z3
Z4
Z5
Blessures
Blessures graves
Dommages
pouvant tre
prvisibles aux mortelles dans
plus de 50% des
mortelles
personnes
cas
Blessures
Possibilits de
blessures
Trs faibles
possibilits de
blessures lgres
Dgts trs
graves
Dgts moyens
et lgers
Dgts lgers
Dgts trs
lgers
Dgts
prvisibles aux
biens
Z1
Z2
Dgts
importants
27/72
4 QUELQUES METHODES DE PREDICTION

Il existe encore peu de mthodes pour prdire les effets conscutifs lclatement dun
rservoir. Nous distinguons les mthodes globales, de type phnomnologique, et les mthodes
plus dtailles qui font appel au support numrique.
4.1 PREDIRE LA SURPRESSION DE RUINE

4.1.1 GRANDEUR DE BASE
Une premire estimation de la surpression de ruine peut tre obtenue sur la base des rgles
habituelles de construction (3.2.1). Ainsi, il est prconis dans louvrage Methods for the
calculation of physical effects resulting from releases of hazardous materials (dit par the
Committee for the Prevention of Disasters, Second edition, 1992) et le guide de lUCSIP
(Union des Chambres Syndicales de lIndustrie du Ptrole) pour llaboration de ltude des
dangers, de retenir un rapport de 2,5 entre la pression de calcul (pression de service effective
maximale admissible) et la surpression de ruine.
Des valeurs plus diffrencies sont proposes dans le CODAP (CODe de construction des
Appareils Pression), qui dpendent de la nature du matriau utilis :
2,4 pour les aciers non allis et les aciers allis non austnitiques4,
3 pour les aciers austnitiques,
4 pour le cuivre et les alliages de cuivre,
2,4 pour laluminium, le nickel et les alliages de nickel,
3,5 pour le titane.
Dans la pratique, il est possible que la surpression de ruine scarte parfois notablement de
cette gamme.
4.1.2 EFFET THERMIQUE

Laccidentologie rvle que la temprature a une influence notable sur la pression de rupture.
On considre par exemple un rservoir cylindrique de 10 m3 dont la pression de calcul est
1,7 bar effectifs entirement rempli de gaz. Le diamtre de ce cylindre est de 1,60 m,
lpaisseur de son enveloppe est de 5 mm en accord avec les rgles du CODAP. La pression de
Les aciers austnitiques sont des aciers inoxydables dont la structure cristalline (agencement des atomes
constituant le matriau) est une structure cubique face centre.
28/72
rupture du rservoir la temprature ordinaire de 25C a t considre pour ce calcul gale

6,1 bar absolus (3*la pression de calcul du rservoir).
On imagine que le cylindre est pris dans un incendie, et que, de ce fait la pression du gaz
lintrieur de lenceinte augmente de faon linaire avec la temprature, selon la loi des gaz
parfaits (Figure 7).
La pression de rupture PR du rservoir est relie la contrainte de rupture R par la relation
classique :
PR
R e
R
o R est le rayon du rservoir et e lpaisseur de son enveloppe. partir de cette relation et

des donnes de la Figure 4, qui relie la contrainte limite de rupture la temprature, il est
possible destimer lvolution de la pression de rupture en fonction de la temprature (Figure
7).
9
Pression du gaz
Pression de rupture
7
Pression [bar]
6
5
4
3
2
1
0
0
100
200
300
400
500
600
700
800
900
Temprature [C]
Figure 7 Evolution de la pression du gaz et de la pression de rupture en fonction de la temprature
On trouve alors que les courbes caractristiques de la rupture et de la pression interne se

croisent vers 550C pour une surpression voisine de la moiti de la surpression nominale
dclatement.
On peut ainsi admettre quun rservoir de gaz pris dans un incendie peut clater la moiti de
sa pression maximale de rupture pour une temprature voisine de 500C. De plus, la Figure 4
(Fruitet, 1979) prsente au paragraphe 3.3.3 montre que pour les aciers, la contrainte de
rupture 500C est proche de la moiti de la contrainte de rupture temprature ordinaire.
Ce rsultat, assez sommairement obtenu dans le cas de lexemple prsent ci-dessus, est
conforme au retour dexprience qui montre quun rservoir pris dans un incendie peut clater
29/72
une pression infrieure sa pression de rupture temprature ambiante, comme cela a t

observ par exemple lors de lessai du wagon-citerne pris dans un incendie (paragraphe 2.1).
4.1.3 EFFET DYNAMIQUE

Il est en thorie possible de dterminer la rsistance dune enceinte soumise une surpression
interne quasi-statique par des calculs classiques de rsistance des structures. Mais on rappelle
que la structure est susceptible de se comporter de manire diffrente lorsque la pression nest
plus applique de manire quasi-statique (typiquement 0,01 bar/s) mais de manire dynamique
(vitesses de chargement de lordre de 1 bar/s ou plus). Le comportement de la structure
dpend du rapport entre sa frquence propre et celle de la monte en pression due
lexplosion. Pour de nombreuses situations, on peut considrer que seul le mode fondamental
de la structure, cest--dire sa frquence caractristique de vibration, est affect (Proust,
2000).
Pour ces situations plus simples, il peut suffire de rduire la structure un systme masseressort comme voqu dans le paragraphe 3.3. Un paramtre caractristique permet de
dterminer assez simplement linfluence dun chargement dynamique sur la rupture de la
capacit. Ce paramtre, DLF (Dynamic Load Factor), reprsente le rapport entre les
chargements dynamique et statique lorigine de la rupture du matriau, cest--dire le rapport
entre les contraintes la rupture pour un chargement dynamique et pour un chargement
statique ; il sagit dun paramtre adimensionnel.
La figure suivante reprsente les variations de la rsistance dune membrane encastre de
dimension caractristique 0,5 m, expose aux effets dexplosion de violences varies
(membrane installe sur une chambre de 10 m3 lintrieur de laquelle ont t amorces des
explosions de poussires).
Figure 8 Dynamic Load Factor pour une membrane fine soumise des explosions de poussires de
violence varies (Proust, 2000)
30/72
On constate alors que leffet dynamique de lexplosion sur une structure peut se traduire soit
par un accroissement soit par une rduction de la rsistance de lquipement considr par
rapport un chargement statique damplitude quivalente (Proust, 2000). En pratique, il est
possible de considrer, de faon forfaitaire, que la pression de rupture en dynamique peut tre
comprise entre la moiti et le double de la pression nominale de rupture en statique.
Ce comportement dynamique augmente le risque sous ses deux aspects. En premier lieu, la
probabilit dclatement augmente puisquil devient possible ds que la pression est proche de
la moiti de la pression de ruine nominale (en statique) alors que, hors ce type de sollicitation
dynamique, le danger dclatement deviendrait possible au voisinage de la pression de ruine
nominale. En second lieu, les effets de lclatement sont susceptibles dtre plus importants car
lclatement peut se produire alors que la pression atteint deux fois la pression de ruine
nominale.
31/72
4.2 MODELISATION DES CONSEQUENCES : LES METHODES GLOBALES

Les mthodes permettant de dcrire les consquences dun clatement de rservoir sont
essentiellement fondes sur les principes de la thermodynamique.
4.2.1 LA METHODE TAC-TNT
Cette mthode a t conue lINERIS au dbut des annes 1990 (Proust, 1991). Une
description de cette mthode a t publie en 1996 (Duplantier, 1996).
On ralise une approximation des ondes produites lors de lclatement en supposant
principalement que :
les gaz sont parfaits et idaux,
toute lnergie de pression sert produire des ondes,
prs du rservoir la thorie des Tubes Choc (TAC) sapplique,
plus loin, les ondes ressemblent celles induites par une libration instantane
dnergie, comme pour un explosif, et peuvent tre reprsentes par les abaques
bien connus de lquivalent TNT.
Brode (Brode, 1959) propose un moyen pour faire le lien entre la zone o la thorie TAC
sapplique (champ proche) et la zone o la dcroissance de type TNT convient (champ
lointain). Daprs Brode, le paramtre qui permet de diffrencier le champ proche du champ
lointain est la masse de gaz mr contenu dans lenceinte avant lclatement. Le champ proche est
dfini par le volume hmisphrique dair V0 , de masse m0 et de masse volumique 0 entourant
lenceinte tel que : m0 10mr
Le rayon de lhmisphre (compt partir du centre du rservoir) est donc gal :
R0 1,7 m r
0
1/ 3
(4.2.1.1)
On peut montrer que ce paramtre peut tre raccord aux allers-retours des ondes de dtente
dans la zone Haute Pression ncessaires pour extraire lnergie.
Le champ proche
En champ proche, la pression du front donde de choc ps0 , appele aussi pression de contact,
peut tre dtermine par la formule des "tubes choc". Cette relation permet de calculer la
pression du front donde de choc qui se propage dans un tube rempli dun mlange gazeux.
Elle peut tre applique au cas des ondes de choc ariennes si lcoulement est
monodimensionnel (onde de choc sphrique ou hmisphrique). Cette relation est la suivante :
32/72
2 1
1 1 a0 a Ps 0
p1
1
Ps 0 1 1
p0
2 0 2 0 0 1Ps 0
1/ 2
1 1
(4.2.1.2)
o :
p1 est la pression de rupture de lenceinte (Pa),
p0 est la pression ambiante (Pa),
Ps 0 est la surpression arienne adimensionnelle aprs lclatement : Ps 0
0 est le rapport des chaleurs spcifiques de lair ambiant (0=1,4),
1 est le rapport des chaleurs spcifiques du gaz comprim,
a0 est la clrit du son dans lair ambiant (340 m/s),
a1 est la clrit du son dans le gaz comprim (m/s).
p s 0 p0
,
p0
Pour un gaz parfait :

a1
a0
TM
1 1 0
0T0 M 1
(4.2.1.3)
o :
T0 est la temprature absolue de lair ambiant (K),
T1 est la temprature absolue du gaz comprim (K),
M0 est la masse molaire de lair ambiant (29,0 g/mol),
M1 est la masse molaire du gaz contenu dans le rservoir (g/mol),
0 et 1 sont les rapports des chaleurs spcifiques.
Lquation (4.2.1.2) est rsolue par itrations successives. Des valeurs calcules de Ps 0
sont fournies sur les diagrammes des
Figure 9 et Figure 10. Cette pression Ps 0 est infrieure la pression des gaz comprims P1 (ou
pression de rupture). Cette relation ne peut pas tre utilise pour calculer les effets de pression
dans le champ lointain.
33/72
Ps0/P0
100
a1/a0=1
a1/a0=2
a1/a0=3
a1/a0=4
10
1
1
10
100
P1/P0
Figure 9 Rapport ps0 /p0 en fonction de p1 / p0 et de (a1 /a0 ) pour un gaz di ou tri atomique(1 =1,4)
Ps0/P0
100
a1/a0=1
a1/a0=2
a1/a0=3
a1/a0=4
10
1
1
10
100
P1/P0
Figure 10 Rapport ps0 /p0 en fonction de p1 / p0 et de (a1 /a0 ) pour un gaz monoatomique (1 =1,66)
34/72
Le champ lointain
Dans le champ lointain, cest--dire pour une distance dobservation R suprieure R0, les
caractristiques de londe de souffle ne dpendent plus que de lnergie de la source. Les
abaques du TM 5-1300 (TM 5-1300, 1990) peuvent tre utiliss pour estimer la pression
maximale de londe de choc en un lieu donn. Ces abaques ont t tablis pour des explosifs
condenss poss au sol tels que le TNT. Comme le montre la Figure 11, lvolution de la
surpression maximale p (p=ps-p0 , o ps est la pression absolue engendre par lexplosion
1/ 3
la distance R) est donne en fonction de la distance rduite R mTNT
(o R est la distance
dobservation et mTNT la masse quivalente de TNT du phnomne).
35/72
10
1
10
100
dp en bar
0.1
0.01
lambda : distance rduite = distance au centre de l'explosion / (masse de TNT

dtonante)1/3
Figure 11 Surpression engendre en champ libre par lexplosion dune charge de TNT au niveau du sol
en fonction de la distance au centre de lexplosion
36/72
La masse dquivalent TNT correspond lnergie de pression disponible au moment de

lclatement de la capacit. Lapplication du premier principe de la thermodynamique londe
qui se dplace permet de montrer que lnergie vhicule dans londe arienne correspond
lnergie dite de Brode 5 (Proust, 1991).
E av
p1 p0 V1
1 1
(4.2.1.4)
o :
p0 la pression de lambiante (Pa),
V1 le volume du ciel gazeux (m3),
et 1 le rapport des chaleurs spcifiques du gaz contenu dans le rservoir.
Cette nergie reprsente laugmentation de lnergie interne du rservoir produite par

laccroissement de la pression. Cet accroissement peut tre obtenu soit par une augmentation
de la temprature des gaz (combustion), soit par linjection de gaz supplmentaire.
Il ny a pas de diffrence de principe entre ce phnomne et lmission donde arienne suite

la dtente brutale des gaz de combustion dun explosif. Il est donc lgitime de sapproprier les
abaques relatifs la propagation des ondes issues de la dtonation dexplosifs sachant que
lnergie de Brode de lexplosif est trs proche de son nergie de combustion. On dfinit
alors un quivalent nergtique classique pour faire le lien avec les abaques :
mTNT Eav
ETNT
(4.2.1.5)
avec mTNT la masse quivalente de TNT et ETNT lnergie spcifique de combustion du TNT
(ETNT=4690 kJ/kg).
La formule de la dtente adiabatique rversible (ou dtente isentropique) est parfois prconise :
1
p 0 p1 p1 1
EA
1 1 p 0 p 0
(4.2.1.6)
Elle reprsente lnergie transmise lenvironnement par une dtente lente et isentropique dun gaz. Elle nest
en gnral pas adapte pour le cas dune onde de choc car ce phnomne est le plus souvent irrversible et
dynamique. Dun point de vue pratique, gnralement on remarque que lnergie libre par une dtente
adiabatique rversible est infrieure celle de Brode. Choisir la valeur la plus leve pour lnergie est une
position conservative.
37/72
tude dun exemple dapplication de la mthode :
On considre lclatement dun rservoir de 1 m3 dair la pression de 40 bars et la

temprature ambiante. On dtermine dabord la pression de contact proximit du
rservoir (
Figure 9, en considrant la courbe a1 /a0 =1) qui est de 4,8 bar absolus. Ensuite, la distance
frontire entre champ proche et champ lointain est dtermine. La relation (4.2.1.1) indique
qu partir dune distance de 5,8 m, labaque TM5-1300 ci-dessus peut tre utilis ; il est alors
ncessaire de dterminer la masse dquivalent TNT, ce qui implique la dtermination de
lnergie de Brode disponible dans la capacit au moment de la rupture. La masse quivalente
de TNT de la source dexplosion considre est de 2 kg.
On a reprsent sur la Figure 12 lvolution de la surpression en fonction de la distance. En
champ proche, la surpression est considre constante et gale la pression de contact ; ds
que la frontire entre champ proche et champ lointain est franchie, nous appliquons la loi de
dcroissance des ondes de choc sphriques contenue dans labaque TNT.
surpression [bar]
10
1
pression de choc
TNT
0.1
0.01
0.1
10
100
distance au centre du rservoir [m]
Figure 12 Evolution de la surpression engendre par un clatement laide de la mthode TAC-TNT
Lexprience montre cependant que la courbe de dcroissance des surpressions de labaque

TM5 en fonction de lloignement est majorante, notamment pour des surpressions infrieures
10-20 mbar.
tant donn que la dcroissance de la surpression nest pas connue en champ proche, nous
considrons, par prudence, que la surpression est constante et gale la pression de choc dans
tout le champ proche.
38/72
4.2.2 LA METHODE PROJEX

LINERIS a dvelopp la mthode TAC-TNT en y apportant des modifications, elle sappelle
aujourdhui PROJEX. Labaque TM5 a t remplac par la courbe Indice 10 de labaque de la
mthode Multi-Energy (V.d.Berg, 1985) (Figure 13) qui fournit des rsultats plus
vraisemblables notamment pour les trs faibles surpressions ariennes. Les effets de pression
engendrs par la mise latmosphre brutale du contenu dun rservoir lors de son clatement
sont ainsi assimils ceux engendrs lors dune dtonation (la courbe indice 10 correspond aux
effets dune dtonation).
P/Pat
Figure 13 Abaque relatif la mthode Multi-Energie donnant des surpressions engendres par des
dflagrations vitesse de flamme constante de volumes explosibles hmisphriques poss au sol
(Mouilleau, Lechaudel, 1999)
On y a adjoint galement un module de balistique pour dterminer la trajectoire des fragments.

Fondamentalement, le module Missile du logiciel Effex de lINERIS (Proust, 2000) a t
intgr en y ajoutant les effets de la compressibilit de lair pour les grandes vitesses de
projection.
39/72
Ce module permet de calculer :
la trajectoire du fragment dduite du bilan des forces extrieures appliques au

fragment,
les caractristiques du point dimpact au sol telles que : la distance, la vitesse et langle
de la trajectoire du fragment par rapport au sol.
Les mouvements de rotation du corps ne sont pas considrs dans la modlisation du

phnomne.
Les quations reprsentatives type utilises sont tout fait classiques :
Mf
d V f (t )
dt
1
gz C x S f 0 U V f (U V f )
2
avec :
U et Vf , les vecteurs vitesse de lcoulement issu de la brche et du fragment dans un

repre fixe,
z, le vecteur unit de direction verticale dirig vers le haut,
Mf et Sf, les masse et surface du fragment (surface projete dans le plan perpendiculaire
lcoulement),
Cx, le coefficient de trane du fragment, typiquement de lordre de 1 pour une plaque

mince (Duplantier, 1996),
0, la masse volumique de latmosphre.
La comparaison des calculs numriques et de quelques rsultats exprimentaux a montr que

cette mthode de calcul permet destimer les caractristiques du point dimpact avec une
prcision satisfaisante.
Outre les caractristiques du fragment, une donne essentielle est la vitesse initiale. Dans
PROJEX, elle est dduite du calcul de limpulsion de pression communique par le jet de gaz
issu de la brche. Limpulsion est une fonction de la surpression maximale et du temps de
dcharge du rservoir. Dans le cas dun gaz sous pression, ce temps rsulte de la propagation
dondes soniques. Dans le cas dun liquide qui se vaporise, ce temps rsulte de la propagation
dune onde de vaporisation .
40/72
4.2.3 LA METHODE DE BAKER

La mthode propose par Baker est dtaille dans son ouvrage "Explosion hazards and
evaluation" (1983) ; elle permet de dcrire les clatements de rservoirs contenant
exclusivement des gaz sous pression.
Comme pour les mthodes TAC-TNT et PROJEX, il est fait une distinction entre le champ
proche et le champ lointain :
pour le champ lointain, Baker fonde sa mthode sur la similitude avec les effets de la
dtonation dun explosif condens,
en champ proche, la mthode est fonde sur les rsultats de calculs numriques
dclatements de rservoirs contenant des gaz parfaits.
4.2.3.1 LES ONDES DE PRESSION

Le champ proche
Baker aboutit, dans le champ proche, un abaque comme celui reprsent sur la Figure 14.
Cest un graphe double entre qui ncessite le calcul dune part de la distance rduite en
abscisse :
p
R rt 0
Eex
1/ 3
(4.2.3.1.1)
avec Eex lnergie de londe de pression, et dautre part lestimation de la pression de choc Ps
sur la base de la thorie des tubes choc. On se repre ainsi la fois en abscisse et en
ordonne, ce qui permet de slectionner une courbe de dcroissance de la pression dans le
champ proche. Pour slectionner cette courbe, on calcule le paramtre :
3V g
r0
2
1/ 3
(4.2.3.1.2)
qui est la dimension quivalente du rservoir de volume Vg et qui permet daccder labscisse
de dpart :
P
R0 r0 0
Eex
1/ 3
(4.2.3.1.3)
On dtermine galement la pression de choc au tout dbut de la dtente Ps0 au moyen des
graphes additionnels (Figure 15 et Figure 16). Le couple ( R0 , Ps 0 ) permet ainsi didentifier la
courbe de dcroissance de la surpression en champ proche.
41/72
Le champ proche est dfini approximativement pour R 2 .

Eex est fonde, comme pour les mthodes INERIS, sur lnergie de Brode. Cependant, les
abaques de Baker (Figure 14 et Figure 17) sont tablis pour des gomtries totalement
sphriques. Dans la ralit, les rservoirs clatent sur le sol, dans un demi-espace, et pour
pouvoir utiliser ces graphes, il faut multiplier lnergie de Brode par 2.
Figure 14 Evolution de la surpression adimensionnelle en fonction de la distance adimensionnelle

(Duplantier, 1996)
42/72
Figure 15 Rapport (a1 /a0 )=1 T1 M0 /0 T0 M1 en fonction de p1 /p0 et de Ps0 (pression de choc au tout dbut
de la dtente) pour un gaz di ou tri-atomique (1 =1,4) (Duplantier, 1996)
Figure 16 Rapport (a1 /a0 ) en fonction de p1 /p0 et de Ps0 pour 1 =1,66 (Duplantier, 1996)
43/72
Le champ lointain
Dans le champ lointain, cest--dire pour R 2 , labaque de la dcroissance de surpression
pour un explosif condens est utilis (Figure 17).
Ps
R
Figure 17 Evolution de la surpression en fonction de la distance pour une explosion de pentolite (TNO,
1997)
Baker a propos des raffinements de sa mthode pour tenir compte plus prcisment de la
gomtrie du rservoir ou de sa position par rapport au sol.
4.2.3.2 LA PROJECTION DES FRAGMENTS

Baker propose une mthode permettant de connatre les trajectoires des fragments en fonction
de la pression du rservoir. Les paramtres dentre sont le nombre de fragments et leur
coefficient de trane (Cx).
44/72
Dans un premier temps, la vitesse initiale des fragments est estime, elle dpend du rservoir et
de la manire dont il se dpressurise. Ensuite, la trajectoire, qui ne dpend que de cette vitesse
initiale, est dtermine.
4.2.3.2.1 Estimation de la vitesse initiale du fragment

Baker a fait une synthse des rsultats exprimentaux et des mthodes numriques antrieures
1983. Cette analyse concerne des essais raliss dans les conditions suivantes :
les enceintes sont soit sphriques soit cylindriques,
lpaisseur de lenceinte est constante,
le rapport longueur sur diamtre (L/D) des enceintes cylindriques est gal 10,
lenceinte est remplie dun gaz parfait,
lclatement de lenceinte est provoqu par linjection sous haute pression dun gaz
(hydrogne, air, argon, hlium, dioxyde de carbone),
les fragments obtenus (2, 10, 100) sont de taille identique. Les diffrentes configurations
dclatement sont reprsentes sur la Figure 18.
45/72
a) Eclatement d'une sphre en 2 fragments gaux
b) Eclatement d'une sphre en plusieurs fragments

gaux
c) Eclatement d'un cylindre en 2 fragments gaux
d) Eclatement d'un cylindre en plusieurs

fragments gaux
e) Eclatement d'un cylindre en 2 fragments ingaux

Figure 18 Fractionnement des sphres et des cylindres en plusieurs fragments (Duplantier, 1996)
Les rsultats de cette analyse ont permis Baker (1983) dlaborer un diagramme reproduit
sur la Figure 19. Ce diagramme permet de dterminer la vitesse initiale adimensionnelle dun
fragment connaissant la pression adimensionnelle de lenceinte P juste avant lclatement
(valable pour P 0, 2 ) dfinie par la relation suivante :
P
p1 p0 V
ma12
o :
p0, la pression ambiante (Pa),
a1, la clrit du son dans le gaz comprim (m/s),
m, la masse totale de lenceinte.

46/72
La vitesse du son a1 dans le gaz comprim est gale :

a12
TR
o :
M
R, la constante des gaz parfaits (8,32 J/Kmol),
T, la temprature absolue du gaz juste avant lclatement (K),
M, la masse molaire du gaz (kg/mol),
, le rapport des chaleurs spcifiques.
Vif
K a1
P
Figure 19 Vitesse du fragment en fonction de la pression et du nombre de fragments (Duplantier, 1996)
47/72
Ces courbes, qui reprsentent le rapport entre la vitesse initiale du fragment Vif et la vitesse du
son dans le gaz a1 en fonction de la surpression adimensionne P , sont traces en fonction du
nombre de fragments attendu. Il est donc ncessaire de connatre le nombre de fragments pour
les exploiter. Lorsque les fragments sont de forme et de masse identiques, le coefficient K est
gal 1. La Figure 20 permet de dterminer ce coefficient K dans le cas de lclatement dun
cylindre (avec Lv/dv=10) en deux fragments ingaux.
Fraction massique du
rapport la masse totale
fragment
par
Figure 20 Coefficient dajustement K pour lclatement dun cylindre (avec Lv/dv=10) en deux fragments
de taille ingale (Duplantier, 1996)
4.2.3.2.2 La distance de projection des fragments

La distance de projection rsulte de simulations numriques des quations classiques de la
balistique. Baker (1983) a reprsent les rsultats sur la Figure 21 pour une altitude de dpart
nulle. Sur ce diagramme, lvolution de la distance adimensionnelle R est reprsente en
fonction de la vitesse adimensionnelle Vi pour diffrents rapports portance sur trane
CLAL/CDAD.
48/72
Figure 21 Distance maximale R en fonction de la vitesse initiale et du rapport des coefficients de trane
et de portance du fragment (Duplantier, 1996)
a C D AD R
Mf
a C D AD vi2
Vi
Mfg
o :
AD est la surface du fragment dans le plan perpendiculaire la trajectoire (m),
AL, la surface du fragment dans le plan parallle la trajectoire (m),
CD, le coefficient de trane,
CL, le coefficient de portance,
g, lacclration de la pesanteur (m/s),
Mf, la masse du fragment (kg),

49/72
R, la distance entre le point djection et le point de retombe au sol (m),
a, la masse volumique de lair ambiant (kg/m3 )
vi, la vitesse initiale du fragment.
Le coefficient de trane pour diffrentes formes de projectiles est donn dans des tables
(Tableau 6, daprs Baker, 1983). En gnral, le coefficient de portance est suppos nul.
50/72
Tableau 6 Coefficient de trane pour diffrentes formes de projectiles (Duplantier, 1996)
4.2.4 LA METHODE DE LUFIP

La mthode prconise par lUnion Franaise des Industries Ptrolires a t initialement
propose par "The Committee for the Prevention of Disasters" dans son ouvrage "Methodes
for the calculation of physical effects resulting from releases of hazardous materials" (1992),
appel communment "TNO Yellow Book". Elle est dveloppe dans le "Guide
Mthodologique UFIP pour la ralisation des tudes de Dangers en raffineries, stockages et
51/72
dpts de produits liquides et liqufis" (UFIP, 2002). Elle est fonde sur une approche
purement nergtique du phnomne.
4.2.4.1 BILAN DENERGIE

Le bilan nergtique du phnomne scrit :
Ei E r E pd E fr
et E r E f E sh
o :
Er est lnergie rsiduelle disponible pour lmission des fragments et de londe de choc,
Ei est lnergie interne du gaz (calcule laide de la relation de Brode),
Epd est lnergie de dformation plastique du rservoir,
Efr est lnergie de rupture du rservoir,
Ef est lnergie disponible pour la propulsion des fragments,
Esh est lnergie dissipe dans londe de choc.
Lnergie absorbe par la rupture de lenceinte (Epd +Efr ) est fonction de la fraction de
dformation du matriau constituant lenceinte la rupture.
Une distribution de lnergie rsiduelle Er est rpartie forfaitairement entre Ef et Esh de la faon
suivante :
E f FE r
E sh 1 F E r
avec :
F=0,2 pour une rupture fragile,
F=0,6 pour une rupture ductile.
52/72
4.2.4.2 LES ONDES DE PRESSION

La mthode de calcul des ondes de pression rsultant de lclatement dune capacit sous
pression prconise par lUFIP utilise le concept dquivalent TNT ainsi que labaque TM-5.
Pour les clatements de rservoirs atmosphriques, la mthode prconise par lUFIP est
dutiliser la courbe 6 de la Figure 22. Cette mthode est proche de la mthode de Baker sauf
que les abaques ont t tablis apparemment pour un demi-espace et non lespace complet, et
comportent la fois le champ proche et le champ lointain (Figure 22). Le paramtre de
labscisse est la distance rduite dfinie par :
R
E sh Pa 1 / 3
Par rapport la mthode de Baker, on constate que les effets de pression sont fonction du
rapport des vitesses du son mme dans le champ lointain, ce qui soulve des interrogations.
53/72
Figure 22 Courbes du modle du TNO pour des pressions dexplosion allant de 5 bars 500 bars
considrant le paramtre Cex /Co , rapport des vitesses du son la pression dclatement et la pression
atmosphrique (daprs UFIP, 2002)
54/72

Lnergie Ef, fraction de lnergie disponible, est suppose tre dissipe sous la forme
dnergie cintique. Par consquent, la vitesse initiale du fragment est gale :
2E f
u
mr
f
0
1/ 2
o mr est la masse du rservoir.
La trajectoire du fragment est ensuite calcule partir des quations classiques de la balistique
prenant en compte les forces de gravit et les forces de dynamique des fluides que sont la
trane et de portance. Une solution simplifie semi-empirique est propose.
55/72
4.2.5 CONCLUSIONS SUR LES METHODES GLOBALES

Il existe dautres mthodes plus confidentielles qui traitent partiellement tel ou tel aspect de
lclatement ; cest le cas en particulier de ce que propose Baum (Baum, 2001 : Annexe 2).
LINERIS propose dans ce paragraphe une comparaison des mthodes prsentes ci-avant en
statuant sparment sur les effets de pression et de projection. Il est important de rappeler que
pour toutes, la pression de rupture de lenceinte est une donne dentre. La dtermination de
ce paramtre est un sujet de recherche en soi.
Les points de comparaison concernent dune part le respect de la physique, dautre part les
approximations choisies ainsi que les limitations associes, et enfin lergonomie dont dpend le
risque derreur.
4.2.5.1 LA PREDICTION DES EFFETS DE PRESSION

Les mthodes TAC-TNT, PROJEX et Baker respectent la physique de lmission et de la
propagation des ondes de pression dans lenvironnement avec une phnomnologie
diffrencie entre le champ proche et le champ lointain. Il existe un questionnement pour la
mthode UFIP pour laquelle lINERIS ne sexplique pas pourquoi les courbes dattnuation
dans le champ lointain restent fonction du rapport des vitesses du son. La mthode UFIP
rpartit, de faon forfaitaire, lnergie disponible pour la pressurisation de lenceinte jusqu sa
rupture (nergie de Brode) entre londe de pression arienne dune part et la projection des
fragments dautre part. Nous ne percevons pas la raison physique pour cela sachant que la
projection des fragments nest pas une question nergtique, ou en tout cas ne peut pas tre
relie lnergie de Brode contrairement aux effets de pression. La prise en compte de ce
rendement peut conduire une nette sous-estimation de lnergie prise en compte pour
calculer le champ de pression.
Nous proposons une comparaison des rsultats sur lexemple suivant :
Rservoir de propane qui clate sous leffet dune explosion interne
Masse molaire du propane : 44,096 g/mol
Contenu : 10 m3 de propane liquide + 90 m3 dun mlange de gaz et dair
Rapport des chaleurs spcifiques des produits brls : 1,314
Temprature des produits brls (temprature de rupture) : 2000C
Pression de rupture : 42,5 bars
Temprature de lair ambiant : 20C
Rapport des chaleurs spcifiques de lair : 1,40
Masse volumique de lair : 1,225 kg/m3
Pression atmosphrique : 1,013 bar

56/72
Dans le tableau suivant sont prsents les rsultats pour les distances aux seuils de 50 mbar et
140 mbar.
Seuil de
Distance calcule Distance calcule
Distance
Distance
surpression
laide de la
laide de la
calcule laide calcule laide
[mbar]
mthode TAC- mthode PROJEX de la mthode
de la mthode
TNT [m]
[m]
de Baker [m]
UFIP [m]
50
140
120
120
90
140
70
50
60
40
Tableau 7 Comparaison des distances calcules avec les diffrentes mthodes aux seuils de 50 et 140 mbar
On constate un bon accord entre les mthodes TAC-TNT, PROJEX et Baker, en logique avec
la proximit des descriptions physiques qui sont faites dans chacune des mthodes. La mthode
TAC-TNT fournit des rsultats lgrement en excs. En effet, labaque TM5 sur-estime
gnralement les niveaux de pression une distance donne, de ce fait, les deux autres
mthodes sont certainement plus prcises. En revanche les rsultats de la mthode UFIP sont
en-dessous en raison de la prise en compte du rendement F.
Ces estimations supposent toutes que les parois de confinement disparaissent instantanment,
cela permet dutiliser lanalogie avec la dtonation des explosifs. En ralit, les brches ne se
forment pas instantanment, lenvironnement est sollicit moins brutalement que dans le cas de
lexplosion dun explosif. Ainsi, la vitesse de sollicitation du milieu est en ralit plus lente, et
leffet de pression devrait par consquent tre plus faible. De ce point de vue, la mthode
PROJEX est certainement plus facile faire voluer car elle fait rfrence la mthode
MULTI-ENERGIE qui prvoit intrinsquement une relation entre les effets de pression
distance et la vitesse de libration de lnergie ; des dveloppements sont donc possibles.
Les mthodes graphiques sont manuelles et prsentent une marge dinterprtation parfois trop
grande, ainsi, par exemple, le choix de la courbe dvolution de pression en fonction de la
distance de la Figure 22 peut tre difficile pour des gammes de pression intermdiaires et des
rapports de vitesse du son compris entre 1 et 5, ce qui reprsente la majorit des situations.

Dun point de vue physique, ces mthodes tiennent a priori bien compte des phnomnes. Une
comparaison sur la base dessais (Baum, 2001 : annexe 2) est propose sur la Figure 23
partir des conditions fournies dans le tableau ci-aprs.
Les essais de Baum ont consist tudier la projection dun rservoir cylindrique contenant un
gaz sous pression et prsentant une fissure axiale initiale.
57/72
Diamtre [m]
Longueur [m]
Masse du
missile [kg]
Pression de
rupture [bar]
Rservoir 1
0.102
0.305
9.54
52.7
Rservoir 2
0.102
0.305
6.32
61.7
Rservoir 3
0.102
0.305
6.32
23.8
Rservoir 4
0.102
0.305
9.22
27.2
Rservoir 5
0.102
0.305
9.1
72.8
Rservoir 6
0.305
1.09
62.35
61
Rservoir 7
0.305
1.22
82.5
76.9
Tableau 8 Caractristiques des diffrents rservoirs utiliss dans les essais de Baum
Pour les cas considrs ici, le missile principal peut tre assimil au rservoir entier, la dtente
du gaz se fait radialement et la temprature dclatement t prise gale 20C.
Vitesse initiale du rservoir [m/s]
120
100
80
essais (Baum)
PROJEX
Baker
UFIP
60
40
20
0
Rservoir Rservoir Rservoir Rservoir Rservoir Rservoir Rservoir
1
2
3
4
5
6
7
Figure 23 Comparaison des vitesses djection dun rservoir entier calcules avec les diffrentes
mthodes globales et celles mesures au cours des essais de Baum
A priori la mthode de Baker paratrait relativement raliste parce que fonde sur des essais,
on saperoit pourtant que les rsultats sloignent parfois significativement de lexprience
mme si les estimations restent conservatives. Les rsultats de la mthode UFIP sont encore
plus conservatifs mais nous mettons en doute le parti pris de raisonner en terme dnergie,
58/72
alors que la projection des fragments est une affaire dimpulsion (Proust, 2000). Finalement, la
mthode PROJEX parat la plus proche.
4.3 MODELISATION DES CONSEQUENCES : LES METHODES NUMERIQUES

Les mthodes numriques permettent de dcrire avec prcision la gomtrie du problme
tudi par lintermdiaire dun maillage, et de prendre en compte la prsence dobstacles. Il est
lgitime de penser que les codes de mcanique des fluides pourraient tre utiliss avec une
certaine confiance pour traiter les questions de propagation donde dans lenvironnement. Il
nexiste pas notre connaissance de codes commerciaux qui intgreraient galement leffet de
rupture et de projection de fragments. Nous ne traitons dans ces paragraphes que de laspect
onde arienne sachant que les aspects de balistique sont assez correctement traits par les
mthodes globales. Deux codes de mcanique des fluides (CFD) ont t mis en uvre sur une
situation qui a t traite exprimentalement. Les codes sont FLUENT et AUTOREAGAS.
4.3.1 AUTOREAGAS
Ce code est propos pour rsoudre des problmes dexplosion et de propagation donde de
pression sur la base du code REAGAS dvelopp par TNO. La partie onde de pression de
AUTOREAGAS propose une rsolution simplifie en volumes finis des quations dEuler
(viscosit nglige) en employant la mthode classique dite FCT (Flux-Corrected Transport)
qui introduit une contre-diffusion artificielle pour limiter les instabilits numriques de la
mthode de rsolution. On note que pour simplifier la rsolution les volutions dans chaque
direction de lespace sont rendues relativement indpendantes les unes des autres en ngligeant
certains termes croiss.
4.3.2 FLUENT
FLUENT est un code de mcanique des fluides vocation industrielle et, comme tel, est
particulirement adapt pour les coulements complexes et turbulents. Il rsout les quations
de Navier-Stockes et propose un maillage adaptatif pour mieux dcrire les angles et les
discontinuits. La rsolution, en volumes finis, est directe sans avoir recours une correction
de la diffusion numrique. Les quations sont rsolues dans leur globalit sans omettre
notamment les termes croiss.
4.3.3 PERFORMANCES DES LOGICIELS

Les logiciels prsents ci-dessus ont t mis en uvre afin de calculer le champ de pression
gnr par lclatement dun rservoir. Les rsultats numriques ainsi obtenus ont ensuite t
confronts aux rsultats de la simulation exprimentale correspondant la mme gomtrie et
la mme chelle, prsente dans le paragraphe suivant.
59/72
4.3.3.1 LA SIMULATION EXPERIMENTALE

La situation test est celle de lclatement dune capacit (rservoir n0) dans un environnement
encombr. La gomtrie de lessai est dcrite sur la Figure 24 ; celle-ci correspond la
reproduction une chelle rduite, dune installation industrielle existante. Le rservoir n0
perd ses deux fonds mais conserve sa virole. Londe de pression qui en rsulte se propage vers
un ensemble de rservoirs paralllpipdiques poss sur le sol. Les pressions ont t mesures
sur les faces latrales et sur le dessus de ces rservoirs (Figure 25).
Rservoir n0
Rservoir n1
Rservoir n2
Rservoir n3
Figure 24 Schma descriptif de la gomtrie considre (dimensions relles mesures en dizaines de

mtres)
Les capteurs sont positionns sur chacun des rservoirs n1, 2 et 3 comme sur le schma
suivant :
Capteur 3
Capteur 2
Capteur 1
Capteur 4
Capteur 5
Figure 25 Emplacement des capteurs sur chacun des rservoirs n1, 2 et 3
Il existe des situations qui ne laissent que peu de marge lincertitude ou tout du moins pour
lesquelles la marge derreur doit tre connue et contrle. Il est alors intressant de songer
mettre profit les lois de similitude des ondes de choc bien tablies (Baker, 1983) pour raliser
des expriences sur maquettes, lchelle rduite.
Ces lois indiquent que les niveaux de surpression incidente Pi ne sont fonction que du rapport
E/X3 o X est la distance au centre de l'explosion. Ds lors, si "D" reprsente la dimension
caractristique de l'installation relle (D par exemple peut-tre la largeur d'un four soit 15-16
60/72
m) et "d", la dimension correspondante de la maquette, alors les niveaux de surpression

mesurs en tout point de la maquette seront semblables ceux atteints sur l'installation relle si
l'nergie de pression Ed dlivre par la source d'explosion, est telle que :
Ed = Eb . (d/D)3
Ainsi, pour la situation relle de la Figure 24, nous avons utilis cette loi pour reprsenter, sur
la maquette, chaque rservoir par un bac en matire plastique, large de 75 cm et haut de 35 cm.
Les trois bacs sont aligns sur un sol plan et disposs 70 cm les uns des autres. La source
dexplosion est dispose dans un cylindre de diamtre 100 mm, de longueur 300 mm, simulant
le rservoir n0, suspendu environ 70 cm laplomb du bac simulant le rservoir n1 (Figure
26). La source de pression dans le cylindre est simule par la dtonation dune petite charge
dexplosif, juste suffisante pour produire instantanment dans ce cylindre la pression requise.
Plaque
simulant
le toit
Cylindre
simulant le
rservoir n0
Bac1 = rservoir 1
Bac2 = rservoir 2
Bac3 = rservoir 3
Figure 26 Installation exprimentale
Les niveaux de surpression atteints sur le toit et les faces latrales de chacun des trois bacs ont
t mesurs. Linterprtation dans le cadre de lapplication industrielle est tout fait simple car
les valeurs mesures sont celles qui ont d se produire dans la ralit. Seule limpulsion des
ondes doit tre calcule pour tenir compte de lnergie rellement disponible. Le niveau de
prcision est estim de lordre de 10%. Les rsultats sont parfaitement cohrents et seront
utiliss dans le paragraphe suivant vis--vis de la comparaison avec les codes CFD.
61/72
4.3.3.2 COMPARAISON AVEC LES RESULTATS NUMERIQUES

Une comparaison des rsultats exprimentaux et des calculs ( la mme chelle que celle de la
simulation exprimentale) est propose sur les Figure 27, Figure 28 et Figure 29.
Surpresssion max calcule/observe (mbar)
180
160
140
120
100
Fluent
Reagas
Essais
80
60
40
20
0
Capteur 1
Capteur 2
Capteur 3
Capteur 4
Capteur 5
Figure 27 Comparaison des niveaux de surpression maximale calculs et observs sur le bac n1
62/72
120
100
80
60
Fluent
Reagas
Essais
40
20
0
Capteur 1
Capteur 2
Capteur 3
Capteur 4
Capteur 5
35
30
25
20
Reagas
Essais
15
10
5
0
Capteur 1
Capteur 2
Capteur 3
Capteur 4
Capteur 5
On constate dune manire gnrale que les tendances sont assez bien respectes par les
logiciels, notamment pour ce qui concerne la prise en compte des rflexions, mais que les
valeurs sont sous-estimes, parfois nettement. Pour AUTOREAGAS, on note une assez
grande dpendance des rsultats la finesse du maillage. Pour FLUENT, une optimisation du
63/72
maillage a permis de limiter cet effet mais laccord nest pas meilleur. Nous pensons que lcart
par rapport lexprience tient assez largement la difficult pour ces codes de capturer
les ondes de pression et les discontinuits. En particulier, londe de choc, infiniment mince, est
contenue au minimum dans une cellule. Plus elle est grande, moindre sera le niveau de
pression. lheure actuelle, les codes commerciaux ne disposent gnralement pas de sousmaillages (traitement interne une cellule de la position des interfaces par exemple) et
semblent donc tous possder ce dfaut difficilement contournable.
4.4 LES LIMITES DES MODELES EXISTANTS

Trois grandes catgories de modles plus ou moins complets ont donc t identifies :
Les modles phnomnologiques (TAC-TNT, PROJEX, Baker, UFIP), qui dcrivent

directement avec plus ou moins dapproximation les phnomnes physiques, ont comme
principal avantage dtre en principe assez faciles daccs. Dans le domaine spcifique de
lclatement de rservoir, on peut estimer quils couvrent bien la physique. Ils ne tiennent
cependant pas compte de la rupture progressive du confinement, ni des dtails
gomtriques qui peuvent renforcer ou au contraire attnuer le champ de pression.
Les modles numriques de type CFD (Computational Fluid Dynamics) sont capables de
dcrire plus prcisment la gomtrie. Nous avons test des codes du commerce, lun ddi
aux explosions (AUTOREAGAS), et lautre prvu pour les questions de mcanique des
fluides (FLUENT). Une comparaison avec lexprience est propose. Les codes semblent
prdire qualitativement les tendances et notamment lattnuation en fonction de la prsence
dobstacles. Cependant, il subsiste un dsaccord trs significatif sur les niveaux de pression
que lon peut attribuer clairement lincapacit de ce type de code capturer les sauts
de pression. La projection des fragments nest pas prise en compte.
La simulation exprimentale met profit les lois de similitude bien tablies des ondes de
choc pour autoriser la ralisation dexpriences lchelle rduite. Un exemple est
propos, qui dmontre la faisabilit du principe. On peut penser que le recours la
maquette serait nanmoins limit des situations pas trop complexes. lheure actuelle, la
question de la projection des fragments nest pas prise en compte mais cela devrait tre
possible.
64/72
5 MESURES DE PREVENTION ET DE PROTECTION

Ce rapport ne traite pas spcifiquement de lanalyse des risques, cest--dire des causes
dclatement. Cependant, le retour dexprience fournit des indications qui permettent de
rappeler quelques mesures permettant de se prmunir contre les scnarios les plus dangereux.
5.1 REDUCTION DU RISQUE A LA SOURCE

On rappelle que les effets de pression sont fonction du paramtre PruptV ; la pression de
1/ 3
rupture Prupt de lenceinte et le volume gazeux V sont les principaux facteurs influenant la
quantification du phnomne et ils ont tous les deux la mme importance. Ainsi, une
diminution de la pression de rupture par deux, volume constant, dans le but de rduire le
potentiel de danger, nentranerait quune rduction modeste des distances deffets (20%). Un
raisonnement analogue peut tre fait propos du volume.
5.2 MESURES DE PREVENTION

Les causes principales de lclatement dun rservoir sont, daprs le retour dexprience :
lincendie, le sur-remplissage, lexplosion interne, la fatigue ou la corrosion.
En dehors des mesures propres limiter le risque dincendie, qui peut tre totalement
indpendant du rservoir, les calculs que nous avons proposs suggrent que la rupture dun
rservoir pris dans un incendie devrait se produire lorsque la pression interne atteint environ la
moiti de la pression de rupture standard. Il est par exemple possible denvisager le calcul
dune soupape tare une pression infrieure ce seuil comme mesure complmentaire.
Pour les situations de sur-remplissage, la mise en place de soupapes ou de siphons pourrait
galement tre prconise.
Pour ce qui concerne les causes mcaniques (fatigue, corrosion), il nest pas vident que les
classiques preuves hydrauliques constituent une rponse adapte. La mise en pression la
pression dpreuve pourrait faire apparatre ou agrandir des fissures sur des quipements
vieillissants et induire ultrieurement sans que rien ne le laisse prsager, une rupture. Dans ces
situations, on pourrait comme alternative proposer une auscultation par chographie, par
mesure dpaisseur et par contrle de ltat de la matire afin de dterminer dans quelles
mesures lquipement a vieilli.
propos de lexplosion interne, en dehors des mesures propres viter ce scnario, il est
possible denvisager la pose dvents de dcharge ou linstallation dun systme dextinction
automatique de faon ce que la pression maximale de lexplosion ne dpasse pas, a priori, la
moiti de la pression de rupture standard.
65/72
5.3 MESURES DE PROTECTION

Il sagit de se protger contre les effets des ondes de pression et des fragments.
La mise en place de merlons permet de diminuer, trs localement, la surpression dun facteur 2
3 (Rglementation pyrotechnique). On pourrait imaginer ventuellement fragiliser des zones
prfrentielles sur l'quipement de manire laminer suffisamment la fuite la brche et
bnficier deffets directionnels, mais il sagit dun sujet dtudes futures.
propos des fragments, on pourrait envisager des zones de rupture prfrentielle permettant
la pose dun systme de retenue. Alternativement, linstallation dun mur de clture pourrait
limiter les angles de projection des incidences quasi-verticales de faon rduire les distances
de projection dun facteur 2 par exemple. Enfin, la protection dinstallations sensibles est
possible grce des crans pare-clats dont le dimensionnement est fonction du nombre de
fragments incidents, de leur vitesse et de leur masse.
66/72
6 CONCLUSIONS
Le retour dexprience montre que lclatement dun rservoir peut tre la consquence de la
fatigue de lenveloppe, dune corrosion excessive, dune explosion interne, ou encore dune
augmentation plus lente et accidentelle de la pression sous leffet dun chauffement, dun surremplissage
Les effets sur lenvironnement de lclatement dun rservoir sont dune part lmission dune
onde de pression, et dautre part la projection des fragments. Londe de pression rsulte de la
dtente brutale du gaz contenu dans le rservoir, ou de la vapeur si le rservoir contient un
liquide surchauff.
Parmi les mthodes de prdiction des effets conscutifs lclatement dun rservoir, on peut
distinguer les mthodes globales des mthodes plus dtailles qui font appel au support
numrique.
LINERIS a dvelopp, au cours des dernires annes, la mthode PROJEX qui permet de
dterminer les consquences dun clatement de rservoir :
La pression de rupture du rservoir est une donne et le mode de rupture est

suppos connu et instantan.
Dans le champ proche, la pression est gale la pression de choc dtermine

classiquement partir de la thorie des tubes choc.
Dans le champ lointain, la surpression dcrot selon la courbe indice 10 de la

mthode multi-nergie.
La projection des fragments est dtermine laide des quations classiques de la

balistique dans lesquelles la vitesse initiale est dduite des conditions de dcharge
du rservoir. Les effets de compressibilit de latmosphre sont pris en compte.
Dans bien des situations cette mthode fournit des estimations trs conservatives. En
particulier, la vitesse de fissuration lors de la rupture, la masse des fragments et la gomtrie du
rservoir ont une incidence sur la formation de londe et les effets directionnels. Ces aspects
devraient faire lobjet de recherches spcifiques. Dautres points devraient galement tre
abords comme la prsence dun milieu biphasique au lieu dun gaz pur. Enfin, il nest pas
vident, et le retour dexprience aprs accident le confirme, que la loi de dcroissance
suppose (multi-nergie indice 10) soit strictement observe dans la ralit surtout grande
distance de lclatement et faible pression en raison par exemple des effets dhtrognit de
latmosphre.
Parmi les mthodes globales prsentes dans ce rapport, les mthodes TAC-TNT, PROJEX, et
la mthode de Baker respectent la physique de lmission et de la propagation des ondes de
pression avec une phnomnologie diffrencie entre le champ proche du rservoir et le champ
lointain. Cette frontire entre le champ proche et le champ lointain dpend de la quantit de
gaz prsente dans le rservoir au moment de lclatement. En champ proche, la surpression
engendre par lclatement est dtermine laide de la thorie des tubes choc. En champ
67/72
lointain, la surpression est suppose suivre, pour ces trois mthodes, des abaques dexplosion
violente. Par ailleurs, des incertitudes demeurent sur certaines hypothses de calcul du modle
TNO utilis par la mthode UFIP.
Une comparaison avec des essais a montr que les estimations de vitesse initiale des fragments
de la mthode de Baker et de la mthode UFIP taient conservatives alors que celles de la
mthode PROJEX taient assez proches des essais.
Enfin, ce type de mthode globale ne permet pas une rponse fine concernant linteraction des
ondes avec les structures environnantes. En principe, les mthodes numriques seraient plus
adaptes mais les codes du commerce ne semblent pas fournir une rponse approprie. Peuttre que lessai sur maquette serait une solution promouvoir pour ces situations.
La mise en uvre de logiciels commerciaux (AUTOREAGAS et FLUENT) a permis de
constater que les tendances dvolution de la surpression taient assez bien respectes, mais
que ces valeurs taient sous-estimes par rapport aux rsultats exprimentaux. La projection
des fragments nest pas prise en compte dans les codes de mcanique des fluides qui ont t
considrs.
Enfin, moyennant la mise en uvre des lois de similitude des ondes de choc, la ralisation
dexpriences sur maquettes peut tre une solution pour la prdiction de leur propagation.
68/72
7 GLOSSAIRE
BARPI
Bureau d'Analyse des Risques et des Pollutions Industrielles
BLEVE
Boiling Liquid Expanding Vapor Explosion
CFD
Computational Fluid Dynamics
GPL
Gaz de Ptrole Liqufi
PPRT
Plan de Prvention des Risques Technologiques
TNO
The Netherlands Organisation of applied Scientific Research
TNT
Trinitrotolune
69/72
8 BIBLIOGRAPHIE
1-
BAKER, COX, WESTIN, KULESZ, STREHLOW, 1983, Explosion hazards and

evaluation, dition Elsevier.
2-
BAUM, 1999, Failure of a horizontal pressure vessel containing a high temperature

liquid : the velocity of end-cap and rocket missiles, Elsevier, Journal of Loss
Prevention in the Process Industries 12, pp.137-145.
3-
BAUM, 2001, The velocity of large missiles resulting from axial rupture of gas
pressurised cylindrical vessels, Elsevier, Journal of Loss Prevention in the Process
Industries 14, pp. 199-203.
4-
BRODE, 1959, Blast wave from a spherical charge, The physics of fluid, volume 2.
5-
CLANCEY, 1972, 6th Int. Meeting of Forensic Sciences, Edinburgh.
6-
DUPLANTIER, 1996, Les phnomnes dexplosion rsultant de la combustion de gaz

de vapeurs et de brouillards dans des appareils clos, Rapport scientifique et technique
de lINERIS, rf. RST 04.
7-
FRUITET, 1979, Scurit contre lincendie dans la construction mtallique,

Techniques de lIngnieur.
8-
HOLDEN, 1988, Assessment of missile hazards : Review of incident experience

relevant to major hazard plant, United Kingdom Atomic Energy Authority.
9-
LEPRETTE, 2002, Le BLEVE, Phnomnologie et modlisation des effets thermiques,

www.ineris.fr, rf. 25427.
10-
MASSON, 1998, Explosion dun silo de crales, Rapport de synthse INERIS,

www.ineris.fr, rf. 21FP30.
11-
MERCIER, 2003, Etude de linfluence du couplage fluide/structure sur la

fragmentation de capacits mtalliques sous pression, BCRD FRACA Rapport
intermdiaire 1re anne, rf. 45013.
12-
MERCIER, PERRETTE, DAUDONNET, CHAINEAUX, 2003, Physics of steam

explosion, 4th ISFEH (International Seminar on Fire and Explosion Hazards) Londonderry - 8/12 Septembre 2003.
13-
MOUILLEAU, LECHAUDEL, 1999, Guide des mthodes dvaluation des effets

dune explosion de gaz lair libre, www.ineris.fr, rf. 20433.
14-
PROUST, 1991, Distances disolement autour des stockages de liquides inflammables

- Explosion de la phase gazeuse, Document technique interne, Ref. 782996, INERIS.
15-
PROUST, 1997, Le mcanisme des explosions de poussires, Communication

prsente au sminaire INERIS-EUROFORUM de mars 1997 Paris.
70/72
16-
PROUST, 2000, Prvoir les effets des explosions de poussires sur lenvironnement,
EFFEX, un outil de simulation, www.ineris.fr, rf. 22751.
17-
RABBE, LIEURADE, GALTIER, 2000, Essais de fatigue, Techniques de lIngnieur.
18-
TM 5-1300, 1990, Structures to resist the effects of accidental explosions,

Departments of the Army, the Navy, and the Air Force, TM 5-1300, NAVFAC, P-397,
AFR 88-22, Washington, DC, Novembre 1990.
19-
T.N.O., 1989, Methods for the Determination of Possible Damage to People and
Objects resulting from Hazardous Materials (Green Book), Report of the committee
for the prevention of Disasters.
20-
T.N.O., 1997, Methods for the Calculation of the Physical Effects of the Escape of
Dangerous Material (Yellow Book), Report of the committee for the prevention of
Disasters.
21-
UFIP, 2002, Guide Mthodologique UFIP pour la ralisation des tudes de Dangers
en raffineries, stockages et dpts de produits liquides et liqufis, Union Franaise
des Industries Ptrolires.
22-
WRIGHT J. K., 1961, Shock Tubes, John Wiley and sons, Londres.
71/72
9 LISTE DES ANNEXES
Repre
Dsignation prcise
Nb pages
Annexe 1
Exemples daccidents ayant donn lieu des clatements de

capacits de stockage
7 pages
Annexe 2
Exprience de Baum
7 pages
72/72
Annexe 1 :
Exemples daccidents ayant donn lieu des
clatements de capacits de stockage
A1/A2
EXEMPLES DACCIDENTS AYANT DONNE

ECLATEMENTS DE CAPACITES DE STOCKAGE
LIEU
DES
1. ECLATEMENTS DE RESERVOIRS PRIS DANS DES INCENDIES

1.1 ESSAI DE RUPTURE DUN RESERVOIR PRIS DANS UN FEU
Le 28/07/73, White Sands Missil Range, New Mexico, USA, un essai de rupture dun
wagon citerne de 120 m3 de propane a t ralis.
Le wagon-citerne, spcialement conu pour lessai, a t construit identique ceux utiliss
pour le transport dammoniac, de GPL, de chlorure de vinyle mise part la prsence
dquipements de mesures. Il a t soumis un incendie ; la monte en temprature a
entran une augmentation de la pression interne du rservoir.
Un flux de chaleur moyen sur la surface mouille de la citerne de 105 kW/m a t mesur.
La soupape de scurit sest souleve aprs 2,2 min lorsque la pression de 19 bars a t
atteinte. Elle a effectu trois cycles avant de rester ouverte.
La pression interne avait atteint 24 bars au bout de 24,5 min lorsque le rservoir a clat.
Les pression et temprature internes atteintes au moment de lclatement taient infrieures
celles dclenchant les dispositifs de scurit.
127 fragments provenant du wagon-citerne et des quipements ncessaires aux diffrentes
mesures ont t retrouvs, leurs positions ont t rpertories. Parmi ces fragments, 63
proviennent du wagon-citerne en lui-mme et une dizaine (9, 10 ou 11) proviennent de la
citerne en elle-mme, les distances et directions de ces derniers ont t rpertories. La plus
grande distance parcourue par un fragment est de 407 m.
tant donn que le wagon tait plac dans une fosse de 8 m de profondeur, les fragments qui
sont retombs lextrieur de cette fosse ont parcouru des distances au sol plus faibles que
celles que lon aurait observes si lessai avait eu lieu au niveau du sol. Toutefois, il est
difficile denvisager des corrections permettant de saffranchir de la prsence de la fosse
tant donn que la distance parcourue par un fragment dpend de deux paramtres (sa
vitesse et son angle djection) et que les renseignements sur les hauteurs relatives des
sources et des points datterrissage ne sont pas disponibles pour la plupart des accidents.
Les thermocouples indiquent un gradient de temprature considrable le long du toit de la
citerne, au niveau de la rupture. Au point damorce de la fissure, sur lune des extrmits de
la citerne, une temprature de 640C a t releve alors que sur lautre extrmit, elle
slevait seulement 450C. Le point chaud correspond au point probable damorce de la
fracture.
Les temprature et pression de rupture ainsi que le temps coul correspondent bien aux
caractristiques de contrainte de rupture du matriau.
1/7
Exploitation :
La citerne a explos une pression infrieure la pression de dclenchement des dispositifs
de scurit. Ceci sexplique par la monte en temprature du matriau qui entrane une
diminution du niveau de contrainte de rupture accept par le rservoir.
Le nombre de fragments provenant de la citerne proprement dite est assez faible (9, 10 ou
11).
Le gradient de temprature sur la citerne au moment de la rupture est assez important ; ceci
nous montre que la rpartition de temprature lorsquun rservoir est pris dans un incendie
peut tre trs htrogne.
1.2 ECLATEMENT DUN EPURATEUR DE GAZ CAUSTIQUE DANS UNE RAFFINERIE (CANADA
OU USA, 04/70)
Lclatement du rservoir est d la chaleur dun feu engendr par linflammation dune
fuite par une chaudire. Le rservoir a clat suite laffaiblissement de lenveloppe au
niveau du ciel gazeux. Les deux extrmits du rservoir se sont envoles, lune a t projete
plus de 17 m avant de percuter le mur dune maison.
Remarque :
Cet accident illustre le fait que lorsquun rservoir contenant du liquide est pris dans un feu,
son enveloppe saffaiblit plus vite au niveau du ciel gazeux. En effet, cette partie de
lenveloppe schauffe plus vite car les changes thermiques seffectuent moins bien avec la
phase gaz.
1.3 ECLATEMENT DE RESERVOIRS PRIS DANS UN FEU DANS UNE RAFFINERIE (CANADA OU
USA, 08/72)
Une fuite de butane sest enflamme au niveau des chaudires situes environ 30 m. Le feu
a ensuite envelopp une tour de fractionnement (= colonne de distillation ?) et deux
rservoirs sous pression situs proximit.
Un bidon de 3,5 m de diamtre et 10 m de long a clat ; puis, 10 15 minutes plus tard, le
second bidon de 1 m sur 5 m a galement clat. La plus grosse portion de ce rservoir a t
jecte ; elle a percut une canalisation quelle a rompue entranant ainsi un rejet de fuel
alimentant le feu. Une tour verticale, affaiblie par la chaleur est tombe en travers des deux
principales canalisations.
2/7
1.4 ECLATEMENT DUN WAGON-CITERNE DAMMONIAC (CRESTVIEW FLORIDE 8/04/79)

Le draillement dun train a engendr un feu suite au dversement de diffrents matriaux
inflammables. Le 59me wagon, contenant de lammoniac, a t entaill de part et dautre de
la citerne et a clat 3 minutes aprs le draillement. La cause attribue la propagation de
la fissure est la dformation subie suite laugmentation de pression interne. La citerne sest
clate en deux sections qui ont t projetes jusqu 200 m et 75 m de part et dautre du
rail.
En considrant la position et lorientation du wagon les plus probables avant la rupture, le
fragment qui a t projet 200 m a probablement t projet selon une direction plus ou
moins axiale. Lautre fragment a t projet sous un angle de 104 par rapport cette
direction.
Le 56me wagon, qui contenait aussi de lammoniac, sest renvers en dplaant son toit et
endommageant ses soupapes au point de bloquer lvacuation. Une extrmit de la citerne a
t prise dans un feu dactone et de mthanol. Au bout de 20 minutes, le 56me wagon sest
rompu en expulsant son contenu et en projetant un fragment environ 200 m et deux autres
fragments moins loin.
1.5 PRODUCTION ET DISTRIBUTION DE COMBUSTIBLES GAZEUX (MEXIQUE, SAN JUAN

IXHUATEPEC, 19/11/84)
Le site de stockage tait constitu de 4 sphres de G.P.L. (mlange 80 % butane, 20 %
propane) d'un volume unitaire de 1600 m3 , de deux sphres de 2400 m3 de volume unitaire
et de 48 rservoirs cylindriques horizontaux de diverses capacits. Au moment de l'accident,
approximativement 11000 12000 m3 de G.P.L. taient stocks sur le site.
Lors d'une phase de remplissage d'un rservoir, une canalisation de 8 pouces (200 mm de
diamtre) sous 24 bar se rompt. Lorsque la hauteur du nuage atteint une hauteur visible
d'environ deux mtres, il s'enflamme, 5 10 minutes aprs le dbut de la fuite, sur une
torchre situe quelques 120 150 mtres du lieu du rejet. Le nuage inflammable ayant
vraisemblablement pntr dans des maisons, son inflammation entrane leur destruction.
Quelques minutes aprs l'inflammation du nuage, deux des plus petites sphres donnent
matire des BLEVE, engendrant la formation d'une boule de feu (d'un diamtre valu,
sans certitude, 350 ou 400 mtres) ainsi que l'jection d'un ou deux rservoirs cylindriques.
Les effets thermiques et l'jection de fragments entranent, par effet domino, d'autres
BLEVE.
En dfinitive, les quatre petites sphres furent dtruites. Les sphres plus grosses restrent
intactes, bien que leur supportage ft bris. Seulement 4 des 48 cylindres demeurrent dans
leur position initiale. Dans un rayon de 300 mtres, toutes les personnes furent tues ou
blesses. Cet accident causa au total le dcs de plus de 500 personnes. Il y eut environ 7000
blesss, et 39000 personnes vacues. Les secours mobilisrent de l'ordre de 4000
personnes. Des fragments de sphres furent retrouvs plus de 600 mtres et 12 cigares3/7
rockets furent projets des distances pouvant atteindre 1200 mtres.
Remarque :
Les cigares-rockets sont constitus par des parties des cylindres horizontaux ; ces fragments
se sont envols en embarquant du liquide qui, en se vaporisant, entrane la formation dune
force de propulsion. Par consquent, ces fragments sont jects plus loin que ceux qui sont
uniquement propulss par lclatement.
1.6 RAFFINAGE DE PETROLE (COREE DU SUD, ULSAN, 22/07/90)

Dans une entreprise ptrolire, une fuite sur un systme d'tanchit d'un stockage de butane
s'enflamme, vraisemblablement cause d'une dcharge d'lectricit statique. Sous l'effet de
l'incendie, le rservoir explose. Plus de 10 000 riverains sont vacus et la circulation sur les
voies alentours est coupe. Un rservoir voisin est touch et explose.
1.7 RUPTURE DUNE CITERNE MOBILE DANS UNE CHAUDRONNERIE (CHARENTE,

CHATEAUBERNARD, 20/08/91)
Une explosion et un incendie se produisent dans une chaudronnerie industrielle. Un camionciterne s'enflamme, le stock de bouteilles de propane et de butane explose et creuse un large
cratre dans un hangar. Le feu se propage une cartonnerie voisine. Les pompiers luttent
2 h pour teindre l'incendie ; 8 000 m de cartons sur bobine sont dtruits. Les murs de
parpaings s'effondrent et la charpente mtallique se dforme sous l'effet de la chaleur. Les
experts valuent la pression interne de la citerne 15 bars lors de son ouverture (la pression
de service tait de 16,4 bars et la pression dpreuve de 24,6 bars) ; un chauffement de sa
paroi a sans doute rduit la rsistance la rupture du mtal. Le mlange air-propane (800
kg) s'est enflamm et aurait form une boule de feu de 37 57 m ; ce phnomne aurait
perdur lgrement moins de 5 s. Les dommages matriels internes et externes sont
respectivement valus 4 et 82 MF.
Remarque :
Cet accident illustre le fait quune augmentation de la temprature entrane une diminution
de la rsistance la rupture du mtal constituant lenveloppe de la citerne.
1.8 INDUSTRIE CHIMIQUE DE BASE (VENEZUELA, EL TABLAZO, 26/06/94)

Une fuite d'un hydrocarbure lourd survient sur une conduite de transfert pendant le
chargement d'un camion-citerne partir d'un rservoir. Le produit entre en contact avec un
4/7
lment chaud d'une station de pompage et s'enflamme. Le feu s'tend alors 2 rservoirs
(propane et propylne) qui explosent. Les rservoirs auraient t en phase de premire rpreuve quinquennale, ou en prparation. L'incendie est teint en 1 heure par le service de
scurit avec l'aide des pompiers de la ville proche d'Altagracia ; 8 sauveteurs sont
lgrement blesss lors de l'explosion.
1.9 COMMERCE DE GROS DE PRODUITS INTERMEDIAIRES NON AGRICOLES (ITALIE,

TREVISE, 10/03/96)
Dans un tablissement de stockage et d'embouteillage de GPL, une fuite de propane survient
durant le dpotage d'une citerne routire de 35 m3 vers un rservoir fixe. L'alarme est
donne. Alors que les pompiers se prparent disperser le nuage avec les lances eau, une
explosion dtruit les bureaux de l'tablissement. La partie suprieure de la citerne s'ouvre et
une petite boule de feu se forme. La chaleur provoque l'explosion d'une citerne routire de
12 m3 et dtruit du matriel des pompiers. Un fragment de la citerne endommage le toit d'un
btiment et atterrit 500 m. Trois employs et 10 pompiers sont blesss, l'un des employs
dcde 4 jours plus tard ; 250 personnes du voisinage sont vacues durant 24h.
2. ECLATEMENTS
DE RESERVOIRS SUITE A UNE MISE SOUS PRESSION
ACCIDENTELLE
2.1
RUPTURE DUNE COLONNE DECHANGEUR DE CHALEUR DE GNL (INDONESIE

14/04/83)
Une colonne dchangeur de chaleur de GNL de 47 m et de 5 m de diamtre maximal a t

purge avec un hydrocarbure gazeux chaud et sec afin de la dgivrer et dshydrater. La
pression maximale de service de la colonne tait de 2 bars et la soupape tait tare 4 bars.
La source du gaz de purge tait plus de 35 bars. Le dysfonctionnement dune soupape a
entran la pressurisation de la colonne. Une pression de 5 bars a t enregistre ; 20 minutes
plus tard, la colonne sest rompue. Le rservoir sest fragment en trois morceaux
principaux et en plusieurs morceaux plus petits.
2.2 ECLATEMENT DUN RESERVOIR ACCIDENTELLEMENT MIS EN PRESSION DANS UNE

RAFFINERIE (USA OU CANADA, 01/66)
Un dchargement accidentel de vapeur haute pression dans un rservoir a entran une
surpression dans le rservoir. Un cleanout port ( ?) a t projet 28 m. Les dgts sur
les quipements situs proximit furent mineurs, bien que le rservoir ait t dplac de sa
5/7
position initiale et quune canalisation ait t dtruite.

2.3 ECLATEMENT DUN REACTEUR DANS UNE RAFFINERIE (CANADA OU USA, 09/70)
Lclatement du racteur est d une fracture fragile initie sur des fissures existantes lors
dun essai de mise en pression avec de lazote 27 bars. Une fraction de la tte du rservoir
a t projete 81 m.
3. ECLATEMENT DE RESERVOIR SUITE A UNE EXPLOSION INTERNE

3.1 ECLATEMENT DUN RESERVOIR PAR CHANGEMENT RAPIDE DE PHASE DANS UNE
RAFFINERIE (USA OU CANADA, 11/62)
De leau tait prsente dans le fond dun accumulateur probablement en raison dune panne
de lalimentation au dmarrage. La prsence des deux phases a entran une vaporisation et
une surpression soudaines. La partie principale de laccumulateur a t projete 96 m et un
morceau de 3 m sur 10 m a t projet 220 m.
3.2 ECLATEMENT DUN SECHOIR SUITE A UN EMBALLEMENT DE REACTION (ROYAUME

UNI, 1981)
Au Royaume Uni, en 1981, un schoir rotatif a clat en deux suite une explosion de
vapeur. Le plus gros morceau a t propuls travers lusine et a atterri sur le parking.
Lexplosion de vapeur a t attribue une raction en chane impliquant de leau
anormalement prsente dans le fond du schoir. En effet, la pompe vide utilise
habituellement pour vacuer leau de rservoir faible temprature na pas fonctionn ce
jour-l. Le schoir contenait de lhydrosulfate de soude (Na2S2 O4 ), qui est un sel qui ragit
lentement, mais de faon exothermique avec leau dans les conditions ambiantes. Du fait de
laugmentation de temprature, la raction sest emballe et leau rsiduelle sest
brutalement transforme en vapeur. (Mercier et al., 2003)
4. EXPLOITATION
Il est difficile dexploiter les informations concernant les accidents passs car, gnralement,
elles ne sont pas exhaustives. En effet, trs peu dclatements pneumatiques dus une simple
monte en pression du rservoir ont t recenss, et ceux qui ont t recenss ne fournissent
pas toujours toutes les informations ncessaires la comprhension du phnomne.
Par exemple, les distances parcourues par les fragments sont souvent donnes, mais la masse
6/7
des fragments nest gnralement pas indique.

De plus, les accidents recenss ne prcisent pas toujours la nature de lclatement.
Notamment dans le cas dun rservoir pris dans un incendie, il est souvent difficile de savoir
sil sagit dun clatement pneumatique ou dun BLEVE.
Dans laccident dcrit ci-dessous, il semble que divers types dclatement ont eu lieu :
Draillement dun train (USA, Mississippi, Laurel, 25/01/69) :
Un train de 51 wagons de propane a draill dans la ville de Laurel. Deux wagons se sont
percs, l'un au niveau de la tte, lautre au niveau de lenveloppe de la citerne ; un feu sest
immdiatement dclar. Douze autres wagons se sont peut-tre rompus de faon violente, en
clatant. Le premier sest rompu entre 6 et 20 minutes aprs le draillement, le dernier entre
45 et 70 minutes aprs le draillement. Parmi ces douze wagons, 8 ou 9 ont subi des
BLEVE, un a t rompu par une extrmit provenant dune autre citerne, et deux ont subi
un impact mcanique lors du draillement qui a entran un feu et donc laugmentation de la
pression interne jusqu la rupture.
Des rapports de laccident mentionnent la prsence de feu dans des btiments situs plus
de 460 m. Ces feux ont t engendrs par des fragments de citernes en feu propulss par la
dtente du propane quils contenaient encore.
7/7
Annexe 2 :
Exprience de Baum
A2/A2
EXPERIENCE DE BAUM
Baum a dvelopp un modle simple (Baum, 2001) permettant de prdire la vitesse du
rservoir cylindrique contenant un gaz sous pression lorsque celui-ci subissait une rupture
axiale. Le phnomne considr dans ce modle est celui dcrit au paragraphe 3.4.
Baum sest donc efforc prdire la vitesse du rservoir lorsque celui-ci tait projet en
raison de la propulsion du gaz sous pression initialement contenu. La dmarche adopte par
Baum est trs intressante car il a compar les rsultats dune tude exprimentale aux
rsultats prdits par son modle.
1. LETUDE EXPERIMENTALE
Le dispositif exprimental est dcrit sur la Figure 30. Le rservoir utilis pour lexprience a
un dfaut axial usin sur la surface extrieure de lenveloppe au milieu de sa longueur.
Alimentation en gaz
Rservoir
Dfaut axial
Figure 30 Dispositif exprimental (Baum, 2001)
Pour chaque exprience, le rservoir est plac horizontalement sur un chssis de manire
ce que le dfaut soit situ sur la partie la plus basse du rservoir. Lapprovisionnement en
gaz se fait par un tuyau flexible et lger afin que lon puisse considrer le rservoir comme
tant non retenu. Le rservoir est pressuris avec de lair.
Le dveloppement de la brche au point le plus bas du rservoir entrane le rservoir tre
pouss verticalement vers le haut.
Les expriences ont t ralises sur deux rservoirs en acier :
un rservoir de 0,102 m de diamtre et 1,6 mm dpaisseur,
un rservoir de 0,305 m de diamtre et 6,4 mm dpaisseur.
Sur les 24 expriences ralises, le rservoir est rest en un seul morceau 18 fois. Trois
1/7
1
expriences ont donn lieu au dtachement dune extrmit, le missile principal est alors
constitu par la robe du rservoir et lextrmit reste attache. Trois expriences ont donn
lieu au dtachement des deux extrmits, le missile principal est alors la robe du rservoir.
Les vitesses observes pour les extrmits dtaches taient gnralement semblables celle
du missile principal. Les donnes des expriences sont rassembls dans le Tableau 9 qui
fournit galement, pour chaque cas, la pression de rupture et la vitesse du missile principal.
Une exprience ralise par Herzerg et al est galement mentionne dans ce tableau.
2/7
Diamtre du Longueur
rservoir
du
[m]
rservoir
[m]
Masse du
rservoir
[kg]
Pression de
rupture
[bar]
Type du missile
principal
Vitesse du
missile
principal
[m/s]
0,305
9,54
52,7
Rservoir entier
2,9
0,305
9,6
119,6
Robe+une extrmit
10
0,305
6,32
61,7
Rservoir entier
8,75
0,305
6,32
23,8
Rservoir entier
5,95
0,305
9,1
45,1
Rservoir entier
6,2
0,305
9,19
35,5
Rservoir entier
4,76
0,305
6,29
35,5
Rservoir entier
7,44
0,305
6,33
20,3
Rservoir entier
5,24
0,305
9,22
27,2
Rservoir entier
5,24
0,305
8,99
32,0
Rservoir entier
5,58
0,305
6,26
37,6
Rservoir entier
9,33
0,305
6,29
38,9
Rservoir entier
9,7
0,305
8,93
38,2
Rservoir entier
6,65
0,305
9,22
48,6
Rservoir entier
7,35
0,305
6,29
41,4
Rservoir entier
8,84
0,305
8,97
55,2
Rservoir entier
7,35
0,305
6,21
62,1
Robe+une extrmit
10,6
0,305
9,10
72,8
Rservoir entier
8,23
1,02
57,6
76,9
Extrmit
39
1,09
62,35
61
Rservoir entier
43,9
1,22
76
121,7
Extrmit
61
1,22
82,5
76,9
Rservoir entier
23,2
1,22
82,5
83,8
Extrmit
54,6
62,5
1,22
82,5
87,2
Robe+une extrmit
36,9
72,5
0,5 (Herzerg 1,33

et al., 1981)
156
25
Rservoir entier
16
0,102
0,305
Tableau 9 Donnes et rsultats des expriences
3/7
Masse du
missile
[kg]
5,66
4,21
47,5
57
2. LE MODELE POUR PREDIRE LA VITESSE DU RESERVOIR

Pour prdire la vitesse du rservoir, il considre limpulsion I applique sur la face interne du
ct oppos la brche. La vitesse du rservoir V est ensuite obtenue en galisant
limpulsion I la quantit de mouvement finale du rservoir :
V
I
o M est la masse du rservoir
M
Lestimation approximative de limpulsion est obtenue en considrant que le rservoir est un

paralllpipde base carre de largeur deux fois le rayon du rservoir, et de longueur la
longueur E du rservoir (Figure 31).
Disparition instantane de
ce ct du rservoir
Figure 31 Reprsentation schmatique du rservoir dfectueux
De plus, le dveloppement de la brche est simplifi par lhypothse de la suppression

instantane dun ct de ce rservoir . Cette disparition gnre une onde de dpression
plane dans le gaz sous pression ; cette onde se propage travers le rservoir en gnrant
un coulement travers la brche.
Cette onde se rflchit sur la face oppose de la brche et retourne ensuite dans le plan de
brche. Ce processus est rpt jusqu ce quune seconde onde de dpression soit initie
cest--dire jusqu ce que la premire onde sorte par la brche.
Daprs Baum, la pression sur la face oppose la brche demeure la valeur initiale P0
2R
jusqu ce que londe de dpression arrive, cest--dire jusquau temps
o R est le
a0
rayon du rservoir et a0 la vitesse du son dans le gaz haute pression.
La pression ensuite diminue pendant le processus de rflexion. La dpressurisation a
2R
6R
essentiellement lieu dans lintervalle de temps compris entre
et
, comme le montre
a0
a0
la Figure 32. Par consquent, la plus grande partie de limpulsion a t applique sur le ct
ferm du rservoir lorsque le temps coul partir de lapparition instantane de la brche
4/7
atteint environ
6R
.
a0
P
P0
Figure 32 Evolution de la pression sur la paroi oppose la brche du rservoir
Limpulsion est donne par :

I PAdt o P est la pression agissant sur le cot ferm du rservoir et A laire du ct
ferm du rservoir, soit 2RE.
Si lon approche la diminution de pression par une chute de pression linaire sur la priode
2R 6R
(
,
), on obtient :
a 0 a0
I
8P0 R 2 E
a0
et
8P0 R 2 E
Ma0
Lapproche ci-dessus est valide pour un vnement rel seulement si le taux de

dveloppement de la brche est tel que le ct du rservoir est compltement ouvert au bout
6R
du temps
aprs linitiation.
a0
Si lon considre un dfaut initial central, cela implique que les fractures axiales se propagent
E
vers les extrmits du rservoir, ce qui prend un temps de lordre de
avec Vp la vitesse
2V p
de fracture axiale. Par consquent, le modle peut tre considr comme valide
Ea0
si :
1.
12 RVp
5/7
Pour de lacier, en supposant des fractures ductiles, la vitesse de fracture axiale est de
lordre de 200 m/s et la vitesse des fractures circonfrentielles est typiquement deux fois plus
petite que la vitesse de fracture axiale.
Ainsi, la sparation de lextrmit par des fractures se propageant dans les deux directions
E
R
autour de lextrmit aura lieu au bout dun temps
compt partir de
2V p V p
2
linitiation de la fracture.
6R
, les fragments spars devraient avoir
a0
des vitesses semblables celle prdite pour le rservoir entier. Les expriences ici dcrites
appartiennent cette catgorie.
Par consquent, si cette priode est suprieure
3. COMPARAISON
ENTRE
MODELE
THEORIQUE
ET
RESULTATS
EXPERIMENTAUX
La Figure 33 compare les vitesses de missiles observes avec la vitesse du rservoir entier
prdite par le modle thorique. Les rsultats de lexprience de Herzerg et al. sont
galement mentionns dans cette tude de Baum.
Vitesse
mesure
(m/s)
8P0 R 2 E
Ma 0
Vitesse prdite (m/s)

Figure 33 Comparaison entre les valeurs des vitesses prdites et des vitesses exprimentales
6/7
Bien que les donnes soient quelque peu parpilles, il apparat que les vitesses mesures
pour le rservoir entier et pour les gros fragments ne sont pas trs loignes.
De plus, on constate que les vitesses observes augmentent approximativement linairement
8P R 2 E
avec les valeurs prdites par le modle thorique et que lquation V 0
fournit une
Ma0
limite suprieure raliste par rapport aux vitesses mesures.
En pratique, la brche souvre sur une priode de temps finie et le taux de dpressurisation
est rduit. Dans ces circonstances, un instant donn, la pression agissant sur la paroi
oppose la brche est plus importante et laire sur laquelle la pression instable agit est
diminue.
Le fait que les vitesses des missiles observes soient infrieures celles prdites par
8P0 R 2 E
lquation V
montre que lexistence dun temps douverture fini entrane une
Ma0
rduction de limpulsion applique sur le missile par rapport limpulsion suppose dans le
modle thorique.
7/7

Formalisation Du Savoir Et Des Outils Dans Le Domaine Des Risques Accidentels

Transféré par

Droits d'auteur :

Formats disponibles

Formalisation Du Savoir Et Des Outils Dans Le Domaine Des Risques Accidentels

Transféré par

Informations du document

Description originale:

Titre original

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

Formalisation Du Savoir Et Des Outils Dans Le Domaine Des Risques Accidentels

Transféré par

Droits d'auteur :

Formats disponibles

Formalisation du savoir et des outils dans le

domaine des risques accidentels

Direction des Risques Accidentels

Les clatements de rservoirs,

Ce document comporte 72 pages (hors couverture).

Direction des Risques Accidentels

Ingnieur lUnit ExplosionDispersion de la Direction des

Ingnieur lUnit ExplosionDispersion de la Direction des

Responsable de lUnit ExplosionDispersion de la Direction des

Ingnieur lUnit Rsistance des

Responsable du suivi du volet 1 du

Dlgu Scientifique de la Direction

Directeur des Risques Accidentels

Le prsent document a t tabli :

au vu des donnes scientifiques et techniques disponibles ayant fait lobjet dune

au vu du cadre lgal, rglementaire ou normatif applicable.

Il sagit de donnes et informations en vigueur la date de ldition du document, le 5 octobre

TABLE DES MATIERES

OBJECTIF ET DOMAINE DAPPLICATION

Contexte gnral .................................................................................................................. 6

Rappel succinct sur le phnomne....................................................................................... 7

Incendie sous un rservoir ................................................................................................... 8

Mise sous pression accidentelle.......................................................................................... 11

Affaiblissement mcanique ................................................................................................ 12

Explosion interne ............................................................................................................... 12

Prsentation gnrale du phnomne ................................................................................ 14

Production et propagation de londe de pression .............................................................. 15

La formation des fragments............................................................................................... 20

Lacclration et la projection des fragments................................................................... 22

Les effets sur lenvironnement........................................................................................... 24

QUELQUES MTHODES DE PRDICTION

Prdire la surpression de ruine ......................................................................................... 28

Modlisation des consquences : les mthodes globales .................................................... 32

Modlisation des consquences : les mthodes numriques .............................................. 59

Les limites des modles existants ....................................................................................... 64

MESURES DE PRVENTION ET DE PROTECTION

Rduction du risque la source......................................................................................... 65

LISTE DES ANNEXES

1 OBJECTIF ET DOMAINE DAPPLICATION

rappeler les phnomnes de mise en pression susceptibles dengendrer lclatement dun

exposer les limites de ces mthodes.

1.3 RAPPEL SUCCINCT SUR LE PHENOMENE

1.4 PLAN RETENU

le troisime chapitre est ddi la description du phnomne,

le quatrime chapitre expose quelques mthodes utilises pour la modlisation de ce

le cinquime chapitre prsente succinctement les mesures de prvention et protection

le sixime chapitre expose les principales conclusions de cette tude.

les clatements conscutifs un incendie sous un rservoir contenant un liquide

les clatements conscutifs une mise sous pression accidentelle,

les clatements conscutifs un affaiblissement mcanique,

les clatements induits par une explosion interne.

2.1 INCENDIE SOUS UN RESERVOIR

2.2 MISE SOUS PRESSION ACCIDENTELLE

Azote (essai de Projection dune

2.3 AFFAIBLISSEMENT MECANIQUE

2.4 EXPLOSION INTERNE

une phase de surpression, de dure t+ et damplitude maximale note P+,

une phase de dpression, de dure t- et damplitude maximale note P-.