Propagation D'une Onde de Choc en Présence D'une Barrière de Protection - Sébastien Eveillard
Propagation D'une Onde de Choc en Présence D'une Barrière de Protection - Sébastien Eveillard
Propagation D'une Onde de Choc en Présence D'une Barrière de Protection - Sébastien Eveillard
barrière de protection
Sébastien Eveillard
Laboratoire PRISME
RAPPORTEURS :
Pr. Abdellah Hadjadj Professeur des Universités à l’INSA de Rouen
Pr. Richard Saurel Professeur des Universités à Polytech Marseille
_________________________________________________________________
JURY :
Pr. Azeddine Kourta Professeur des Universités à Polytech Orléans, Président du jury
Pr. Abdellah Hadjadj Professeur des Universités à l’INSA de Rouen
Pr. Richard Saurel Professeur des Universités à Polytech Marseille
Mr. Hervé Jourdren Ingénieur – chercheur au CEA – DAM
Pr. Isabelle Sochet Professeur des Universités à l’ENSI de Bourges
Dr. Jean-Yves Vinçont Ingénieur – chercheur au CEA – DAM, Encadrant
TABLE DES MATIERES
TABLE DES FIGURES ............................................................................................................................................................ I
REMERCIEMENTS............................................................................................................................................................ VII
NOMENCLATURE ............................................................................................................................................................... IX
INTRODUCTION .................................................................................................................................................................... 1
ANNEXES............................................................................................................................................................................. 221
i
Figure 46 μ Evolution de l’effet protecteur en fonction de la HOB (Borgers (2010)), W = 230 kg de TNT,
H = 4,3 m .............................................................................................................................................................. 63
Figure 47 μ Evolution de l’effet protecteur en fonction de la distance d entre le centre de la charge et le mur
(Borgers 2010), HOB = 1,58 m, W = 230 kg de TNT, H = 4,3 m ......................................................................... 64
Figure 48 μ Représentation schématique d’une barrière de protection pour une charge placée au niveau du
sol .......................................................................................................................................................................... 67
Figure 49 : Diagramme de la démarche des travaux de thèse .............................................................................. 69
Figure 50 : Représentation schématique du banc d’essais des expériences à petites échelles ............................. 71
Figure 51 μ Source d’amorçage de la charge explosive........................................................................................ 71
Figure 52 : Représentation schématique de la distance en vue directe de la charge pour des exemples de
positions de capteurs ............................................................................................................................................. 72
Figure 53 : Photo du centre de calcul TERA 100 ................................................................................................. 74
Figure 54 μ Utilisation de l’AMR, suivi du front de choc et de la phase positive de l’onde de choc .................... 77
Figure 55 μ Utilisation de l’AMR dans les zones d’intérêt .................................................................................... 77
Figure 56 : Variation de la taille des mailles en fonction du temps ..................................................................... 78
Figure 57 μ Dimensions du merlon de l’étude préliminaire du projet ANR BARPPRO ....................................... 79
Figure 58 : Présentation des domaines étudiés .................................................................................................... 80
Figure 59 : Profils de surpression pour les capteurs placés en amont du merlon ................................................ 81
Figure 60 : Profils de surpression au niveau du capteur placé sur la face avant du merlon ............................... 83
Figure 61 : Evolution du maximum de surpression réfléchie le long de la face avant du merlon ........................ 84
Figure 62 : Profils de surpression au niveau du capteur sur le sommet du merlon ............................................. 85
Figure 63 : Evolution du maximum de surpression au sommet du merlon (extraction à 5 mm de la paroi) ........ 86
Figure 64 : Profils de surpression au niveau du capteur sur la face arrière ........................................................ 86
Figure 65 : Evolution du maximum de surpression sur la face arrière (extraction à 5 mm de la paroi) ............. 87
Figure 66 : Profils de surpression au niveau des capteurs placés après le merlon .............................................. 88
Figure 67 : Evolution du maximum de surpression en aval de la barrière de protection ..................................... 89
Figure 68 μ Evolution de l’écart entre les différents domaines ............................................................................. 90
Figure 69 : Profil de surpression de l'étude préliminaire du projet ANR BARPPRO (Eveillard (2010)) ............ 94
Figure 70 : Représentation du capteur de référence............................................................................................. 94
Figure 71 : Profils de surpression pour le capteur de référence (CR1) à 0,1 m du centre de la charge
(soit 0,61 m/MJ1/3) ................................................................................................................................................. 95
Figure 72 : Présentation des obstacles étudiés ..................................................................................................... 96
Figure 73 : Présentations schématiques des configurations étudiées ................................................................... 97
Figure 74 : Représentations schématiques des domaines étudiés en géométrie 2D-axisymétrique ...................... 98
Figure 75 : Profils de surpression au niveau du capteur de référence (CR1) pour les configurations 2 et 3 ...... 99
Figure 76 : Réflexion régulière sur la face avant de l'obstacle (configuration 2, t = 0,28 ms) ............................ 99
Figure 77 μ Evolution de la surpression réfléchie sur la face avant de l’obstacle parallélépipédique ............... 101
Figure 78 μ Réflexion sur une paroi verticale avec et sans influence de l’angle d’incidence (β) ....................... 102
Figure 79 : Profils de surpression pour le capteur F0 situé à mi-hauteur de la face avant de l’obstacle
(configuration 3) ................................................................................................................................................. 103
Figure 80 : Détente de l'onde incidente au sommet du parallélépipède pour chaque configuration
(t = 0,35 ms) ........................................................................................................................................................ 103
Figure 81 : Evolution du maximum de surpression au sommet du parallélépipède ........................................... 104
Figure 82 : Profils de surpression pour le capteur T0 situé sur le sommet de la structure ................................. 105
Figure 83 μ Représentation schématique de la détente et de la réflexion sur le sol après l’obstacle de la
configuration 3 .................................................................................................................................................... 105
Figure 84 μ Représentation schématique des distances à parcourir autour de l’obstacle pour la configuration 3
............................................................................................................................................................................ 106
Figure 85 : Représentation schématique des distances à parcourir autour de l’obstacle pour la configuration 2
............................................................................................................................................................................ 106
Figure 86 μ Représentation schématique du contournement de l’onde de choc pour la configuration 2, vue du
dessus .................................................................................................................................................................. 107
Figure 87 : Profils de surpression pour le capteur de référence et évolution de la surpression le long de la face
arrière du parallélépipède pour la configuration 2 ............................................................................................ 108
Figure 88 : Profils de surpression pour le capteur R0 et évolution de la surpression le long de la face arrière du
parallélépipède pour la configuration 3 ............................................................................................................. 109
Figure 89 : Profils de surpression pour le capteur (CA4) .................................................................................. 110
Figure 90 : Profils de surpression pour le capteur de référence (CR1), configuration 4 de Trélat (2006) ........ 111
Figure 91 : Représentation schématique des deux géométries étudiées pour la configuration 4 ....................... 112
Figure 92 : Profils de surpression pour les capteurs placés sur le cylindre à = 0° ........................................ 113
iii
Figure 137 : Evolution du maximum de surpression réfléchie sur la face avant du merlon 2B ......................... 153
Figure 138 : Visualisation sur les champs de pression lors du phénomène de réflexion sur la face avant de la
barrière de protection (configuration 2B_C6_d10) ............................................................................................ 154
Figure 139 : Visualisation sur les champs de pression du phénomène de détente sur le sommet de la barrière de
protection (configurations 1B_C6_d14 et 2B_C6_d14)...................................................................................... 155
Figure 140 : Evolution du maximum de surpression au sommet du merlon 1B .................................................. 156
Figure 141 μ Evolution du maximum de surpression sur le sommet du merlon 2B placé à 0,07 m d’une charge
gazeuse de rayon (R1) (2B_C6_d7) ..................................................................................................................... 157
Figure 142 : Evolution du maximum de surpression au sommet du merlon 2B .................................................. 158
Figure 143 : Visualisation sur les champs de pression du phénomène de détente sur la face aval de la barrière
de protection (configurations 1B_C6_d14 et 2B_C6_d14) ................................................................................. 159
Figure 144 : Evolution du maximum de surpression sur la face aval du merlon 1B .......................................... 160
Figure 145 μ Evolution du maximum de surpression sur la face aval du merlon 2B placé à 0,07 m d’une charge
gazeuse de rayon (R1) (2B_C6_d7) ..................................................................................................................... 161
Figure 146 : Evolution du maximum de surpression sur la face aval du merlon 2B .......................................... 162
Figure 147 μ Représentation schématique des distances parcourues par l’onde pour les merlons de la série 2 163
Figure 148 μ Evolution du maximum de surpression en aval des merlons 1B et 2B placés à 0,07 et 0,14 m d’une
charge de gaz de rayon (R1) ................................................................................................................................ 164
Figure 149 μ Evolution du facteur d’atténuation pour les merlons de la série 2 et le merlon 1A de la série 1
placés à 0,07 m d’une charge de gaz de rayon (R 1) ............................................................................................ 165
Figure 150 μ Evolution du facteur d’atténuation pour le merlon 1B avec plusieurs valeurs de la distance d
(0,070 et 0,140 m d’une charge de gaz de rayon (R 1))........................................................................................ 166
Figure 151 μ Evolution du facteur d’atténuation pour le merlon 2B avec plusieurs valeurs de la distance d
(0,070 ; 0,085 ; 0,100 ; 0,140 m d’une charge de gaz de rayon (R 1))................................................................. 167
Figure 152 : Synthèse des observations issues de l’étude de différentes géométries de barrière de protection . 170
Figure 153 μ Influence du régime de réflexion sur le calcul de l’angle de déviation .......................................... 174
Figure 154 : Représentation schématique de la configuration étudiée pour le phénomène de détente .............. 175
Figure 155 : Abaque des coefficients de détente à partir de résultats de simulation numérique pour une structure
d’angle d’inclinaison de λ0° ............................................................................................................................... 176
Figure 156 : Définition des secteurs géométriques autour de la barrière .......................................................... 176
Figure 157 μ Représentation des paramètres géométriques pour l’onde de contournement .............................. 178
Figure 158 μ Points d’intersection sur les arêtes de la face avant pour l’onde issue du merlon ........................ 179
Figure 159 μ Détente de l’onde de choc pour le secteur 1 et sa réflexion sur le sol ........................................... 179
Figure 160 μ Détentes de l’onde de choc pour les secteurs 2, 3 et 4 ................................................................... 179
Figure 161 : Définition du secteur 3-bis ............................................................................................................. 180
Figure 162 μ Diagramme de fonction de la méthode d’estimation rapide .......................................................... 181
Figure 163 : Extrapolation linéaire du maximum de surpression par un profil simplifié .................................. 182
Figure 164 : Profils de surpression pour un capteur de référence du merlon « court » 2B ............................... 184
Figure 165 : Dimensions du merlon court de la série 1 (merlon 1A) ................................................................. 185
Figure 166 : Positions des différents capteurs placés autour de la barrière de protection du merlon 1A « court »
............................................................................................................................................................................ 185
Figure 167 : Visualisations des champs et des profils de surpression pour les capteurs S 1 et S2 du
merlon 1A « court » ............................................................................................................................................ 186
Figure 168 : Visualisations des champs et des profils de surpression pour les capteurs S3 et S4 du
merlon 1A « court » ............................................................................................................................................ 187
Figure 169 : Visualisation des champs et des profils de surpression pour le capteur S 5 du merlon 1A « court »
............................................................................................................................................................................ 188
Figure 170 : Visualisations des champs et des profils de surpression pour les capteurs S 6 et S7 du
merlon 1A « court » ............................................................................................................................................ 189
Figure 171 : Dimensions du merlon court de la série 2 (merlon 2B) ................................................................. 191
Figure 172 : Positions des différents capteurs placés autour de la barrière de protection du merlon 2B « court »
............................................................................................................................................................................ 191
Figure 173 : Visualisations des champs et des profils de surpression pour les capteurs S 1 et S2 du
merlon 2B « court » ............................................................................................................................................ 192
Figure 174 : Visualisations des champs et des profils de surpression pour les capteurs S 3 et S4 du
merlon 2B « court » ............................................................................................................................................ 193
Figure 175 : Visualisations des champs et des profils de surpression pour les capteurs S5 et S6 du
merlon 2B « court » ............................................................................................................................................ 194
Figure 176 μ Visualisations des zones d’effets pour le merlon 2B à partir des résultats de simulation numérique
et de la méthode d’estimation rapide, vue aérienne plane au-dessus du sol ....................................................... 196
v
REMERCIEMENTS
Je tiens à remercier Madame la Professeure Isabelle Sochet, directrice de thèse, qui
m’a donné l’opportunité de réaliser cette thèse au sein du Laboratoire Pluridisciplinaire de
Recherche en Ingénierie des Systèmes, εécanique, Energétique sur le site de l’ENSI de
Bourges. Je lui suis reconnaissant pour le temps qu’elle m’a accordé dans la direction de ses
travaux et pour m’avoir encouragé tout le long de ces trois années de thèse. Malgré la
distance, j’ai apprécié de travailler avec elle et son équipe, me donnant ainsi l’occasion de
continuer à parfaire ma formation.
J’exprime également mes plus vifs remerciements à εonsieur Jean-Yves Vinçont et à
Monsieur Nicolas Lardjane, ingénieurs de recherche au CEA – DAM, pour avoir accepté
d’encadrer cette thèse. εon travail de thèse doit beaucoup à leur encadrement avisé et
exigeant. Leur vision globale des problèmes de détonique et de mécanique des fluides m’a été
profitable. C’est auprès de mes encadrants que j’ai également pu obtenir l’aide nécessaire
pour mener à bien cette thèse. Je leur suis reconnaissant de leurs qualités pédagogiques,
scientifiques et humaines. Une thèse n’est pas uniquement une « histoire » scientifique, elle a
également une dimension humaine. J’ai énormément appris à leurs côtés. Je tiens donc à leur
exprimer toute ma gratitude.
J’adresse également mes remerciements aux Professeurs Richard Saurel et Abdellah
Hadjadj pour avoir accepté le rôle délicat de rapporteur et le temps qu’ils ont consacré à
examiner ce travail de thèse. Je remercie également les membres du jury et son président.
Je tiens également à remercier Monsieur Pierre-Franck Piserchia, responsable de
laboratoire au DASE, pour l’accueil dans son laboratoire, sa disponibilité et sa contribution
dans l’organisation de la thèse. J’espère ne pas oublier ses nombreux conseils dans le domaine
des relations humaines et du monde de l’entreprise qui me seront profitables à l’avenir.
J’adresse également mes remerciements à εonsieur Christophe εatignon, responsable
de laboratoire à DCSA, pour son cours particulier au début de ma formation, ainsi que ses
conseils et données qui ont contribué à la réalisation de ces travaux de thèse.
Je souhaite également remercier εonsieur Joël Allain, directeur de l’ENSI de
Bourges, εonsieur Kourta Azeddine, directeur du laboratoire PRISεE, et l’ensemble des
partenaires qui ont participé au projet ANR BARPPRO : Air Liquide, le CEA – DEN,
la DGA TT, l’INERIS, l’IRSN et le PNRI.
J’adresse également mes remerciements à εessieurs εichael Bertin, εartin
Henneton, Olaf Gainville, Yves Grillon, Ludovic Aubry, Christophe Millet, et plus
généralement à toute l’équipe de mon laboratoire du DASE grâce à qui cette thèse s’est
déroulée dans des conditions idéales. J’ai apprécié de travailler au quotidien avec cette équipe.
Je remercie également toute l’équipe du laboratoire PRISεE dont εessieurs Xavier
Rocourt, Pierre-Emmanuel Sauvan, Bastien Pellegrineli et Samer Awamat.
Enfin, je tiens à remercier mes amis et ma famille pour m’avoir soutenu tout au long
de ces trois années de thèse. Je remercie mes parents, Jean-Claude et Elise, et ma sœur,
Marie-Laure, pour leur soutien et leurs conseils. Je remercie également mes amis, Grégoire,
Cécile et Marie-Odile, merci de votre soutien. Je conclurais ces remerciements en les
adressant à Ségolène pour sa patience pendant la rédaction de ce mémoire et son soutien au
quotidien.
viii
« Quand il se présente à la culture scientifique, l'esprit n'est jamais jeune. Il est même
très vieux, car il a l'âge de ses préjugés. Accéder à la science, c'est, spirituellement rajeunir,
c'est accepter une mutation brusque qui doit contredire un passé. […]
Avant tout, il faut savoir poser des problèmes. Et quoi qu'on dise, dans la vie
scientifique, les problèmes ne se posent pas d'eux-mêmes. C'est précisément ce sens du
problème qui donne la marque du véritable esprit scientifique. Pour un esprit scientifique,
toute connaissance est une réponse à une question. S'il n'y a pas eu de question, il ne peut y
avoir connaissance scientifique. Rien ne va de soi. Rien n'est donné. Tout est construit. »
« LA FORMATION de l’esprit scientifique », Gaston Bachelard (1934)
ix
Nomenclature
Caractéristiques mécaniques de l’onde de choc
∆P+ : maximum de surpression positive
∆P- : maximum de surpression négative
t0 μ temps d’arrivée de l’onde
dt+ : durée de phase positive
dt- : durée de phase négative
I+ : impulsion positive
I- : impulsion négative
Caractéristiques de la charge explosive
E : énergie libérée
M : masse de la charge
R ou Rcharge : rayon de la charge explosive
Q : énergie chimique libérée lors de la réaction
BO : balance en oxygène
n : nombre de mol
Mm : masse molaire
Vmol : volume molaire
μ énergie d’amorçage
Eamorçage
: coefficient empirique du modèle de Kamlet-Jacobs
Grandeurs de l’onde de souffle ou du milieu
P : pression hydrostatique
V : volume massique
ρ : masse volumique
u : vitesse matérielle
T : température
c : célérité du son dans le milieu
e : énergie interne massique
D μ célérité de l’onde
h : enthalpie
: coefficient polytropique
cp : chaleur spécifique à pression constante
cv : chaleur spécifique à volume constant
M : nombre de Mach
Z : impédance du milieu
μ angle d’incidence
μ distance parcourue par l’onde (en m/εJ1/3)
z μ distance parcourue par l’onde (en m/kg1/3)
: angle de déviation
hm : hauteur du pied de Mach
Ω μ angle d’inclinaison d’une paroi
μ angle d’observation
Cr : coefficient de réflexion
Cd : coefficient de détente
R : constante molaire des gaz parfaits
x
INTRODUCTION
Sur les sites industriels, publics ou privés, la sécurité des biens et des personnes vis-à-
vis des risques associés à une explosion, qu’elle soit d’origine accidentelle ou malveillante,
est une préoccupation majeure de nos sociétés modernes. Cette crainte est visible à travers les
évolutions réglementaires, où le législateur redéfinit ainsi régulièrement les obligations et le
cadre légal en termes de sécurité applicables aux industriels. Ces mesures visent à réduire
l’occurrence d’apparition ou la gravité d’un accident, mais « le risque zéro n’existe pas ». δes
drames d’AZF (en β001 à Toulouse, France), de la dynamiterie Nitrochimie (en β00γ à Billy-
Berclau, France), d’un dépôt de feux d’artifice à Kolding (en 2004, Danemark) ou de l’usine
d’engrais de West (en β01γ au Texas, États-Unis) sont des exemples qui rappellent cette triste
réalité. Pour limiter l’apparition de nouveaux sinistres, les industriels disposent cependant
d’un arsenal de mesures de sécurité permettant de réduire l’un des paramètres inhérent à un
accident : la prévention ou la protection.
δors de l’accident de Kolding, le site industriel était situé à 3 km du centre-ville. δ’explosion
des piles d’explosif situées sur le site (environ 300 tonnes, source ARIA) a causé le décès
d’un pompier et a blessé 199 personnes. Le bilan des dommages s’élève à 100 millions
d’euros. Le site ne disposait pas de barrières de protection.
La présence de barrières de protection, entre les piles d’explosif, aurait pu empêcher la
propagation de la détonation d’une cellule à l’autre, et donc l’effet domino. En plus de
compartimenter le site, les barrières de protection auraient également permis de limiter les
dommages sur le site lui-même. A la suite de cet accident, le Danemark a modifié sa
réglementation et a adopté les standards européens (notamment dans le classement des
explosifs).
Dans le domaine de la protection face aux explosions, il existe différents types de
barrière de protection physique (murs, remblais et talus) qui peuvent jouer un rôle important
contre ces effets : protection des installations et des personnes contre les éclats, les effets
thermiques et les effets de souffle. Si l’effet protecteur est garanti pour les deux premiers
items, il n’en est pas de même pour le dernier. En effet, l’interaction d’une onde de choc avec
une structure est difficile à prévoir, elle dépend de nombreux paramètres et fait l’objet de cette
étude.
Lors d’une explosion, il existe deux régimes possibles : la détonation ou la
déflagration. Cette étude ne traite que l’aspect détonation qui correspond au régime
d’explosion le plus dangereux et le plus « violent » au niveau du maximum de surpression.
Lors d’une explosion en régime de détonation, l’onde de choc générée a une signature
particulière qui réside dans un pic brutal de surpression. δ’évolution de cette onde dépend de
la nature de la charge explosive qui lui donne naissance. En effet, une onde de choc issue de
l’explosion d’une charge de gaz n’évolue pas dans le temps de la même manière qu’une onde
obtenue à partir d’une charge chimique condensée (par exemple du TNT, Trinitrotoluène).
δ’évolution en champ libre (sans obstacle) des caractéristiques mécaniques de cette onde de
souffle (maximum de surpression, impulsion,…) est bien connue et très documentée. Des
abaques d’évolution ou des formules empiriques sont disponibles dans la littérature pour
respectivement une charge de TNT (Kinney (1985), TM5-1300 (2008)) ou une charge de gaz
(Lannoy (1984),…).
Cependant, il n’en est pas de même pour l’ensemble des phénoménologies présentes lors de
l’interaction d’une onde de choc avec une barrière de protection. En effet, lors de sa
propagation au voisinage d’un merlon ou d’un mur de protection, les paramètres de l’onde
vont évoluer en fonction des phénomènes physiques présents : réflexions, détentes,
diffractions, refocalisation,… L’onde subit alors une modification de ses caractéristiques
2
P(t)
∆P+
∆P
I+
I-
∆P-
Temps
∆P
dt+ dt-
t0 t0 + dt+ t0 + dt+ + dt-
auraient été acquises pour la configuration d’intérêt (échelle 1/1) à partir de celles obtenues à
échelle réduite (échelle 1/k, k facteur de dimensionnement) et inversement. La surpression
(ΔP+) est identique entre les deux échelles pour une même position en distance réduite (z ou
). Les autres paramètres en temps et en espace de l’onde de choc varient selon le facteur de
réduction d’échelle, noté « k » et appelé facteur de Hopkinson, figure 2.
m
k 1/1
1/ 3
m1/ k
pour un équivalent massique, par exemple en kg de TNT (4)
E
k 1/1
1/ 3
E1/ k
pour un équivalent énergétique (5)
Ce facteur peut être déterminé à partir d’un équivalent massique ou énergétique, relations (4)
et (5)
k t0
t0 Temps
dt+
dt
+
dt-
dt
-
t0
d’arrivée
k dt+
Durée de phase
dt+
positive
k x t0 Impulsion k I+ I+
k3 m
Masse de la
m
charge
Energie k3 E E
x dt+
k kdt
+
dt-
k xkdt
- Distance z (ou ) z (ou )
z 2
808 1
P 4.5
z z z
(6)
1 1 1
P0 2 2 2
z 10
980 1
dt 0,54
z 3 z 6 z
(7)
1 1 1
W 1/ 3 2
z
0,067 1
4
0,23
I
z 4
(8)
2
z 1
1/ 3
1,55
z 10
980 1
dt
0,54
W 1/ 3 z 3 z 6 z
(9)
1 1 1
2
I
ln 1 / 3 1,190 0,860 ln 0,190 ln pour 0,βλ < < ζ,β7 m/εJ1/3
E
2
(11)
Extrêmement
graves (blessures
Conséquences sur Effets indirects
mortelles dans Très graves Graves Significatifs
l’homme par bris de verre
plus de 50% des
cas)
Saut de pression
ΔP+ ≥ 0,ζγ bar ΔP+ ≥ 0,β0 bar ΔP+ ≥ 0,1ζ bar ΔP+ ≥ 0,0η bar ΔP+ ≥ 0,0β bar
(en bar)
δe classement des zones d’effets peut être comparé aux valeurs de seuils d’effets de
surpression référencés dans la littérature. Cette comparaison est réalisée pour les deux natures
de charge explosive étudiées dans cette étude : charge chimique condensée (TNT, Kinney
(1985), figure 3) et gazeuse (propane-oxygène à la stœchiométrie, Trélat (2006), figure 4)
Figure 3 : Evolution du maximum de surpression et des seuils d’effets de la surpression pour une charge de TNT
en fonction de la distance réduite (z)
10
Figure 4 : Evolution du maximum de surpression et des seuils d’effets de la surpression pour une charge de gaz
(propane-oxygène) en fonction de la distance réduite ( )
l’énergie effectivement dégagée entre l’état amont et l’état aval, notée Q0 :
Q0 Q eex P0 ,V0ex eg P0 ,V0 g (16)
P P0 0 D V0 V
2
(17)
A chaque célérité D correspond un état fin de détonation. Le lieu de ces points se trouve sur la
courbe de « Crussard » dans le diagramme (P, V). Cette courbe définit les états d’équilibre
possibles de la combustion de l’explosif. Cette courbe peut être déduite des équations de
conservation. Une fois éliminés les termes de vitesse, elle constitue alors l’équation de
Hugoniot (relation entre l’enthalpie (h), la pression (P) et la masse volumique (ρ)), en tenant
compte de l’énergie dégagée par la détonation, équation (18).
h h0 P P0 V V1 Q
1
(18)
2
12
Figure 7 : Représentation des positions possibles de la droite de Rayleigh-εichelson par rapport à l’adiabatique
de Crussard (Matignon (2000))
d log P
CJ (coefficient polytropique au point CJ,
d log S
CJ) (20)
c
2
c
2
1 0
CJ 1 DCJ
DCJ
uCJ (23)
Moles de diazote : nN 2
t
2
La combustion d’une mole d’explosif génère ng moles de gaz brûlés par kilogramme
d’explosif :
y z t
ng 10 3
12 x y 16 z 14 t
4 2 2 (27)
z y
18.10 3 44.10 3 28.10 3
y t
Mmg
2 2 4 2
y z y t
(Mmg en kg.mol-1)
(28)
2 2 4 2
δ’énergie dégagée lors de la réaction d’un kilogramme d’explosif est calculée à partir de sa
composition et des énergies de formation :
capacité propulsive d’un explosif :
E g 150 0,112 0
2
(Eg en J/kg) (33)
15
Remarque :
Une comparaison entre les résultats obtenus par un code de thermochimie du CEA – DAM et
le modèle de Kamlet-Jacobs pour une charge gazeuse a été réalisée. Le principe du modèle de
Kamlet-Jacobs pour un explosif condensé peut être appliqué à une charge gazeuse avec
prudence. Les coefficients empiriques des relations simples de ce modèle ne sont pas adaptés
à une charge gazeuse (exemple : différence d’ordre de grandeur sur la masse volumique (ρ0)
et sur le coefficient polytropique ( )). Le modèle de Kamlet-Jacobs est donc difficilement
applicable à la problématique d’explosion gazeuse (coefficient non adapté pour cette
famille d’explosif).
Le modèle de détonation utilisé lors de cette étude se base sur l’approche de Chapman-
Jouguet pour les deux familles d’explosif utilisées (charges de TNT et de gaz (propane-
oxygène à la stœchiométrie)).
t
C-
M (x,t) C+
DCJ t
x
1 x DCJ t
d
u c sur laquelle : du c 0
dx
C+ : (35)
dt
d
u c sur laquelle : du c 0
dx
C- : (36)
dt
δors de la propagation d’une onde de choc, il se forme en aval du front de choc un réseau
d’ondes composé de C+ et de C-, dites « simples », le long desquelles la variation de la
pression et de la vitesse est nulle (dP+ = dP- = 0 et du+ = du- = 0). Ce train d’ondes permet de
décrire l’écoulement en aval du front de détonation.
Si on considère un écoulement isentropique d’un gaz parfait (avec coefficient
polytropique) :
P V P0 V0
(37)
c0 P0 V0 (38)
1
c c0
0
2
(39)
u c et u c cste
1
dx 2
C+ : (40)
dt
u c et u c cste
1
dx 2
C- : (41)
dt
Dans le cas des produits de détonation (voir figure 8), chaque point M, représentant une
particule des gaz brûlés, peut être décrit à partir des deux droites caractéristiques (C+ et C-).
Les droites « simples » C+ sont centrées sur le point d’amorçage en (x = t = 0) et donc
possèdent la même variation de vitesse, de masse volumique et de célérité du son. Les droites
« simples » C- croisent également le front d’onde de détonation. δes courbes caractéristiques
peuvent alors s’écrire :
u c et u c cste
1
x 2
C+ : (42)
t
PCJ 0 CJ
2
1
D
(44)
uCJ CJ
1
D
(45)
DCJ
cCJ
1
(46)
CJ 0
1
(47)
A partir des relations au point de Chapman-Jouguet pour un gaz parfait, les courbes
caractéristiques C+ et C- peuvent s’écrire :
u c et u c cste
1
x 2
C+ : (48)
t
u c et u c CJ cste
1 1
dx 2 D
C- : (49)
dt
Pour chaque point ε situé au croisement d’une C+ et d’une C- :
uc
x
(50)
t
u c CJ
1 1
2 D
(51)
1
On obtient ainsi un système d’équations à 2 inconnues, u et C en fonction de x et t.
1
c cCJ avec
x
DCJ t
(52)
u uCJ 2 1 (53)
c 1 1 1 1
2 2
0 CJ
c0
(54)
1
2 2
c 1 1 1
P P0 PCJ
c0
(55)
δa trajectoire d’un point ε(x, t) situé dans les produits de détonation est obtenue en intégrant
l’équation différentielle :
u uCJ 2 1
dx x
DCJ t
(56)
dt
18
1
t t 1
x 1
1 t0 t0
x0
(57)
Avec un point limite situé à l’interface entre l’air et les produits de détonation :
t0
x0 R charge
(58)
DCJ DCJ
A partir des relations (57) et (58), il est possible d’estimer l’évolution du déplacement d’une
particule située dans les produits de détonation en tenant compte du réseau d’ondes de détente
(ondes de Taylor) pendant la détonation de la charge.
t C+
M (x,t)
C-
Front de choc
U0 = -Γ × t0 t0
x0
x
Une fois que l’onde sort de la charge explosive, le réseau d’ondes de détente est
généré (fin de la mise en mouvement du piston pour les ondes de détente), figure 9. δ’onde
n’est plus sonique dans la charge (D = u + c, n’est plus valide). En aval de l’onde de choc,
l’écoulement est subsonique (D < u + c). δe réseau d’ondes de détente tend à rattraper le front
de choc et à l’atténuer avec une accélération constante (Γ(t0)) (principe du piston de Taylor).
(t ) uCJ 2
du x
DCJ t 2
(59)
dt
A l’interface entre l’air et les produits de détonation, la mise en mouvement du réseau d’ondes
de détente est déterminée :
t 0 uCJ 2 CJ
D
(60)
Rcharge
Dans le milieu industriel, il est possible de rencontrer deux types de charge explosive :
une charge explosive dite « gazeuse » et une charge explosive dite « chimique condensée ».
Les recommandations et les réglementations dans le domaine de la sécurité sont le plus
souvent émises pour une charge de TNT (charge chimique condensée de référence). Lors de
l’utilisation d’un autre explosif, il est possible d’utiliser des équivalents pour se ramener à
l’explosif de référence (TNT). Dans la littérature, plusieurs coefficients d’équivalence sont
disponibles pour la surpression et l’impulsion positives. Ils sont obtenus par la méthode dite
« des dommages » et ils sont valables sur un intervalle donné de distances réduites. Ces
coefficients (Céquivalent) permettent d’obtenir la masse équivalente de TNT (mTNT) qui génère le
19
même niveau de surpression ou d’impulsion positive pour une masse d’un explosif autre que
le TNT (mexplosif), équation (61).
mTNT mexplosif Céquivalent (61)
Figure 10 : Représentation d’une « onde de choc pure », repère lié à l’onde de choc
dh c p T sachant que
cP
(65)
cv
u1 1 u 2 2
c2 c2
1 2 1
1 1
(67)
2
Il est possible de calculer les variations des grandeurs thermodynamiques à travers la
discontinuité en fonction du nombre de Mach en amont de l’onde de choc (M1). Le nombre de
Mach (Mk) est défini comme le rapport de la vitesse matérielle (uk) sur la célérité du son (ck)
dans le milieu, équation (68).
Mk
uk
(68)
ck
Les expressions des variations de vitesse, de masse volumique et de pression peuvent être
déduites de ce système d’équations, relations (θλ) à (71).
u 2 1 M 12 2
1 M 12
(69)
u1
21
2
1 M 12
1 1 M 12 2
(70)
P2 2 M 12 1
1
(71)
P1
Inversement, à partir de l’équation (71), il est possible de connaître le nombre de Mach de
l’onde de choc en fonction de la variation de pression entre l’état amont et l’état aval,
équation (72).
M1 1 2 1
2
1 P
(72)
P1
Cette dernière équation permet de déterminer le nombre de Mach incident (M1 = Mi) de
l’onde de choc à partir du maximum de surpression incident (∆P1+) et de la pression initiale
(P1), avec P2 = P1 + ∆P1+ et = 1,ζ (gaz parfait), équation (73).
1 P P1
M1 6 1 1
7
(73)
P1
y Mi
Onde incidente
Onde réfléchie
Mi μ nombre de εach de l’onde
incidente
Charge explosive (W) x μ angle d’incidence (angle
entre l’onde incidente et la paroi)
Milieu 1 : Z1 = ρ1 x c1 Milieu 2 : Z2 = ρ2 x c2
Interface
Onde réfléchie Onde incidente Onde transmise
b) cas général
Milieu 1 : Z1 = ρ1 x c1 Milieu 2 : Z2 = ∞
Onde incidente
réfléchissante
parfaitement
Surface
Onde réfléchie
Le diagramme de Ben-Dor (1991) présente les différents régimes de réflexion qui ont
été étudiés dans la littérature, figure 12. δ’appartenance à l’un de ces régimes de réflexion
dépend de l’écoulement incident (nombre de Mach incident (Mi),…) et des caractéristiques
géométriques de la paroi (angle d’incidence ( ),…).
Il est possible de distinguer deux familles de réflexion : la réflexion régulière (écoulement
stationnaire) et irrégulière (écoulement pseudo-stationnaire ou instationnaire). La limite entre
ces régimes est déterminée à partir d’un critère de stabilité1 qui dépend du nombre de Mach
incident (Mi) et de l’angle d’incidence ( ).
1
Le calcul de ce critère de stabilité sera développé par la suite (§ 1.3.1.4)
23
Fonction du critère de
Réflexion régulière stabilité Réflexion irrégulière
Dans le cas d’une réflexion régulière, il existe deux régimes de réflexion : normale ou
oblique. Ces régimes se différencient en fonction de l’angle d’incidence ( ). Si l’onde
incidente est normale à la paroi ( = 0), l’onde subit une réflexion normale. A l’inverse, si
l’onde incidente est inclinée par rapport à la paroi ( ≠ 0), l’onde subit une réflexion oblique.
Dans le cas d’une réflexion irrégulière, il existe plusieurs régimes de réflexion. Les réflexions
dites de « chocs faibles » de Von Neuman, Guderley et Vasilev n’ont pas été observées dans
le cadre de cette étude, aussi ces réflexions ne sont pas développées dans ce mémoire. Lors de
l’apparition d’une réflexion de εach, une troisième onde se forme et se propage
perpendiculairement à la paroi, elle est appelée pied de Mach. Suivant la trajectoire et
l’évolution de ce pied de Mach, il est possible de distinguer trois régimes de réflexion de
Mach : réflexion de Mach directe, stationnaire ou inverse. δors du passage d’une réflexion
régulière à une réflexion de εach, un régime temporaire s’établit, il est appelé « Transition
vers une réflexion de Mach ».
Dans le cadre de cette étude, les régimes de réflexion régulière (réflexion normale et
réflexion oblique) et les régimes de réflexion irrégulière (réflexion de Mach) sont analysés.
Les paragraphes ci-dessous présentent ces différents régimes de réflexion.
24
Etat 2 : M2
1
uT un2
u1
u2 u
Φ= - un1 uT
1
Etat 1 : M1
δ’écoulement est considéré supersonique entre deux états (état 1 et état 2). Il est
possible d’établir les relations entre ces états séparés par une onde de choc oblique, figure 13.
La théorie des chocs obliques implique de considérer l’égalité des vitesses tangentielles
(vitesse uT). A partir des composantes normales (indice n) et tangentielles (indice T) des
vitesses de l’état 1 et de l’état 2, il est possible d’établir les relations géométriques suivantes
(équations (74) et (75)).
tan1
un 2
et tan1 u n1 (74)
uT uT
P2 2 M 12 sin 1 1
1
2
(77)
P1
Pour observer un choc réfléchi, il est nécessaire que le nombre de Mach normal à la paroi soit
supérieur à 1 (écoulement supersonique, M1 > 1). Ainsi, il est possible de déterminer l’angle
d’incidence ( ) minimal et maximal à partir de la relation (75), équations (78) et (79).
Angle minimal d’incidence, min arcsin 1
M1
(78)
l’égalité des vitesses tangentielles (uT), équation (80).
2 1 M12 sin 2 1 tan1
1 arctan
(80)
2 2
1 M 1 sin 1
25
Choc réfléchi : Mr
Choc incident : M1
Etat 2 : c2, P2, T2, u2
Etat 1 : c1, P1, T1, u1 Etat 2 : c2, P2, T2, u2
Etat 3 : cr, Pr, Tr, ur
Pr 3 1 M12 3
partir du nombre de Mach incident (M1), équation (81).
CrN
1 M12 2
(81)
P2
Pour un nombre de Mach (M1) infini, le coefficient de réflexion tend vers une valeur
asymptotique (pour un gaz parfait ( = 1,ζ), si ε1 → ∞, alors CrN → 8).
1.3.1.4 Réflexion oblique (RO)
Si l’angle d’incidence ( 1) est supérieur à 0°, l’onde de choc peut subir une réflexion
oblique. Ce phénomène physique est présenté, figure 15.
1
2
Il est possible de distinguer trois états μ l’état initial (en amont de l’onde incidente, noté 1),
l’état « choqué » (en aval de l’onde incidente, noté β) et l’état réfléchi (en aval de l’onde
réfléchie, noté 3). A partir des relations de Rankine-Hugoniot et des angles de déviation ( 1,
2), il est possible de connaître les variations des grandeurs thermodynamiques entre les trois
états.
26
Comme présenté dans la partie précédente sur l’angle de déviation ( 1), cet angle dépend du
nombre de Mach (M1) et de l’angle d’incidence ( min ≤ 1 ≤ max). Dans le cas d’une réflexion
oblique, il est important d’étudier l’évolution de l’angle de déviation au niveau de l’onde
réfléchie ( 2). Pour un nombre de Mach (M2) donné, il est possible de tracer les valeurs de
l’angle de déviation ( 2), figure 16.
2
Choc fort
’’
’ Choc faible
max 2
Figure 16 : Evolution de l’angle d’incidence de l’onde réfléchie en fonction du nombre de Mach et de l’angle de
déviation de l’état β
Pour un nombre de Mach incident (M2) donné, il existe un angle maximum noté max qui
limite les valeurs possibles de l’angle de déviation. Si 1 > max, alors il n’existe pas de
solution pour les chocs obliques ( 2 < 0). Dans ce cas de figure, le choc est incurvé ou se
détache : passage vers une réflexion de Mach (cas de l’étude, > max, critère de stabilité). Cet
angle d’incidence limite ( 2 lim) est donné par la relation suivante en fonction du nombre de
Une formule empirique établie par Kinney (1985) permet d’estimer rapidement l’angle
d’incidence limite ( 1 lim) entre ces deux régimes de réflexion (réflexion oblique et réflexion
de Mach) à partir du nombre de Mach de l’onde incidente (M1), équation (83). Pour un
nombre de Mach infini, l’angle d’incidence limite tend vers un angle proche de 39,9°. δ’angle
minimum est de 39,23° pour un nombre de Mach de 2,48. Cet angle limite entre une réflexion
régulière et une réflexion de Mach est également appelé « angle de transition ».
1 lim 39
1,75
M 1 1
(83)
Si 1 < max, il existe alors deux solutions ( 2’ et 2’’) en fonction de l’angle de déviation ( 2)
et du nombre de Mach (M2). δa solution ’’ est dite solution de choc « fort » (M2 < 1).
Inversement dans le cas de 2’, la solution est dite de choc « faible » (M2 > 1, pour 1 proche
de max). En pratique, on considèrera la solution de choc faible ( 2’). La solution de choc
« fort » ( 2’’) n’est envisageable seulement si un écoulement supersonique (par exemple,
27
l’action d’un piston) est imposé en aval du front de choc incident, ce qui ne sera jamais le cas
dans la problématique de cette étude.
Si 2 lim > 2 > 0 (solution dans les chocs obliques), le nombre de εach de l’onde réfléchie
(Mr) peut être calculé en fonction de l’angle d’incidence 2 et du nombre de Mach de l’état β
M r M 2 sin2
(M2) à partir des composantes des vecteurs vitesses, équation (84).
(84)
δe coefficient de réflexion dans le cas d’une réflexion oblique (Cro) peut être obtenu à partir
des relations de Rankine-Hugoniot et des angles de déviation entre les trois états, équation
Pr 2 M r 1 2 M1 1 1
(85) (Kinney (1985)).
CrO
1 2 M12 1 12
2 2 2
(85)
P2
Etat 3 : cr, Pr, Tr, ur, Mr 2 1 Etat 1 : c1, P1, T1, u1, M1
« Point Triple »
Détachement de 2
l’onde réfléchie M1
de la paroi 1
3 Mp
χ Etat du pied : cp, Pp, Tp, up, Mp « Pied de Mach »
En fonction de la trajectoire suivie par ce point, appelée « ligne de glissement », il est possible
de distinguer trois phénomènes de réflexion de Mach.
Si la trajectoire du point triple tend à s’éloigner de la paroi (dχ > 0), la hauteur du pied de
Mach augmente. Ce régime de réflexion est appelé « réflexion de Mach directe ».
Si la trajectoire du point triple est parallèle à la paroi (dχ = 0), la hauteur du pied de Mach
est inchangée. Ce régime de réflexion est appelé « réflexion de Mach stationnaire ».
Si le point triple tend à se rapprocher de la paroi (dχ < 0), la hauteur du pied de Mach
diminue. Ce régime de réflexion est appelé « réflexion de Mach inverse ». Le régime de
réflexion devient une réflexion oblique lorsque le point triple atteint la paroi
(Ben-Dor (1991)).
δors d’une réflexion de εach directe, le point d’interaction entre l’onde incidente,
l’onde réfléchie et le pied de Mach est détaché de la paroi. Le glissement de ce point est
représenté par l’angle χ. δe point triple devient le point d’interaction entre l’onde incidente et
28
l’onde réfléchie. δe nombre de εach au niveau du pied de Mach (Mp) est estimé au niveau du
point triple par égalité des vitesses tangentielles à la paroi, équation (86).
Mp
sin 1
M1
(86)
A partir des relations de Rankine-Hugoniot entre l’état 1 et l’état en aval du pied de εach
(indice p), il est possible de déduire la valeur de la pression réfléchie sur la paroi à partir du
nombre de Mach du pied de Mach (Mp), équation (87).
2 M p 1
CRM
1
2
Pp
(87)
P2
λM M
Charge (W)
L
O Ω
ξT ξ ξM F
Figure 19 : Représentation schématique de l’onde incidente et de l’onde réfléchie pour une surface plane
Dans le cas d’une onde sphérique issue de la détonation d’une charge explosive (point O,
centre de la charge), l’onde est supposée plane au point d’impact sur la surface (point I,
[OI] ┴ (onde incidente)). On note F le premier point d’impact sur la paroi. Lors de
l’interaction de l’onde incidente avec la paroi, le régime de réflexion peut évoluer depuis une
réflexion régulière (réflexion normale ou oblique, notée RR) vers une réflexion de Mach
(notée RM, h : hauteur du pied de Mach). Ce passage éventuel est ici matérialisé par le
point M.
Le régime de réflexion le long de la paroi peut alors être déterminé à partir de
plusieurs paramètres géométriques et des formules analytiques présentées précédemment,
figure 20.
Charge (W) L
O Ω
ξ F V
(89)
2
La distance entre le centre de la charge et le point d’impact ( ) permet de déterminer le
maximum de surpression incident au point d’impact à l’aide des formules empiriques
présentées précédemment (partie 1.1.3), relations (90) et (91). Le nombre de Mach incident
est ensuite déduit des relations de Rankine-Hugoniot (équation (73)).
Soit sin cos 1
L L
2 2
d d
(91)
d
δe calcul de l’angle d’incidence ( ) et du nombre de εach incident (Mi = M1) permet de
déterminer le régime de réflexion et d’estimer le maximum de surpression réfléchi à partir des
relations analytiques dans le cas d’une réflexion régulière (réflexion normale ou oblique). Les
abaques de réflexion permettent également d’estimer le maximum de surpression réfléchi
dans le cas d’une réflexion régulière ou irrégulière (réflexion de Mach simple). δ’analyse de
l’évolution de l’angle d’incidence ( ) en fonction des caractéristiques géométriques de la
paroi permet d’évaluer l’impact de ces dimensions sur le phénomène de réflexion.
La figure 21 présente l’évolution de l’angle d’incidence ( ) en fonction du rapport L/d
pour plusieurs inclinaisons (Ω). A partir de ce réseau de courbes, il est possible d’analyser
l’influence de l’inclinaison de la paroi sur l’évolution de l’angle d’incidence ( ) pour un
rapport L/d identique. Si l’angle d’inclinaison (Ω) est élevé (proche de λ0°), l’angle
d’incidence évolue rapidement au début de la paroi (jusqu’à δ/d = 4). A partir de L/d = 6,
l’évolution de l’angle d’incidence ( ) se stabilise dans un intervalle situé entre 80 et λ0°. A
l’inverse pour un angle d’inclinaison faible (proche de 10°), l’angle d’incidence évolue
lentement vers la valeur asymptotique de 90° (le réseau de courbes tend vers cette valeur quel
que soit Ω quand δ/d → ∞).
Dans le cas d’une réflexion de εach, le calcul des paramètres géométriques selon une
approche purement analytique est réalisé au niveau du point triple (point T), figure 22.
βT
T
d
h
I
T
βT
H
λM M
Charge (W)
ξT L
O Ω
F
Les relations (92) à (94) permettent de déterminer l’angle d’observation ( T), l’angle
d’incidence ( T) et la distance entre le centre de la charge et le point triple ( T).
h
sin arctan L2 h 2
L
T arctan
cos arctan L2 h 2 d
(92)
h
L
T T (93)
2
h h
T sin arctan L h 2 cos arctan L2 h 2 d
2 2
L L
(94)
Le calcul de ces paramètres dépend de la hauteur du pied de Mach (h). La hauteur (h) peut
être estimée à partir de formules empiriques issues de la littérature (TM5-1300 (2008) et
Kinney (1985)). Celles-ci estiment la hauteur d’un pied de Mach (h) lors de la réflexion de
l’onde incidente avec le sol pour une charge hémisphérique placée au-dessus du sol
(HOB ≠ 0 m, HOB : Height Of Burst), figure 23.
dC tan HOB
2 HOB2
arctan
RR : réflexion régulière
RM : réflexion de Mach Onde incidente
HOB
O Charge (W) RR Onde réfléchie
Pied de Mach
HOB
h
RM
dc
d
Figure 23 : Représentation schématique de la réflexion de l’onde sur le sol
32
La hauteur du pied de Mach peut être estimée en fonction de la hauteur de la charge explosive
(HOB), de la distance parcourue au niveau du sol par l’onde de choc (d) et de la distance
critique (dc) (distance d’apparition d’un pied de Mach sur le sol, transition entre une réflexion
RR et RM). La formule empirique (52) est issue du Kinney (1985).
d
h 0,07 HOB 1 (Kinney (1985))
2
dc
(95)
δ’abaque issu du Tεη-1γ00 (β008) présente l’évolution de la hauteur du pied de Mach pour
plusieurs valeurs de la HOB, figure 24.
Figure 24 : Evolution de la hauteur du pied de Mach (h) à partir de l’abaque du Tεη-1300 (2008)
0,9 dc 3,6
HOB
2
HOB dc2
W 1/ 3
2
0,9 d 3,6 d 0,9 d
c c
h W 1/ 3 100
5,98.105 1 / 3 3,8.103 1 / 3 0,766
(96)
W 1/ 3 HOB HOB
W W
Les formules empiriques présentées ci-dessus peuvent être adaptées aux cas étudiés, afin
d’estimer la hauteur du pied de Mach (h) en fonction des paramètres de la configuration
étudiée (dimensions de la paroi : d = HOB, dc = Lc et d = L). La figure 25 présente une
comparaison des résultats obtenus à partir de l’abaque du Tεη-1300 (2008), de la formule
empirique du Kinney (1985) et de la formule du NRDC (2001) pour quatre distances d
(d = 0,40 ; 0,99 ; 1,58 et 2,78 m/kg1/3). δ’angle d’incidence de transition (passage entre une
réflexion RR et RM) est estimé à 43,8° pour une paroi inclinée à 90°.
Figure 25 : Evolution de la hauteur du pied de Mach en fonction de la distance parcourue sur la paroi, de la
hauteur de l’explosion et pour une charge de 1 kg de TNT
Les résultats obtenus suggèrent que la première formule empirique (Kinney (1985))
sous-estime la hauteur du pied de εach en champ proche (d ≤ 0,99 m/kg1/3). A l’inverse en
champ lointain, la hauteur du pied de Mach est surestimée par rapport à l’abaque du
TM5-1300 (2008) (d ≥ 2,78 m/kg1/3). Pour une distance d proche de 1,58 m/kg1/3, les
estimations obtenues par la première formule empirique se situent dans un intervalle de - 50%
à + 13% par rapport à l’abaque du TM5-1300 (2008). La seconde formule empirique
(NRDC (2001)) permet une estimation satisfaisante en champ proche (d < 1 m/kg1/3), mais
surestime la hauteur du pied de Mach en champ lointain (d >> 1,58 m/kg1/3).
Dans le cas d’une réflexion de εach, l’approche analytique nécessite d’estimer
l’évolution de la hauteur du pied de εach (h). δ’utilisation des abaques disponibles dans la
34
littérature implique que cette approche se base sur la formule empirique issue du Kinney
(1985) et de la formule du NRDC (2001) pour des configurations situées en dehors du réseau
de courbes du TM5-1300 (2008). Suivant la position de la charge (d < 1 m/kg1/3 et
d > 2,78 m/kg1/3), l’approche analytique présentée peut conduire à une surestimation ou une
sous-estimation du maximum de surpression réfléchi : si d est faible, h est sous-estimée, est
sous-estimé, le maximum de surpression réfléchi est surestimé et inversement si d est grand (h
surestimé, surestimé, maximum de surpression réfléchi sous-estimé). Par conséquent les
abaques de réflexion sont à privilégier dans ce cas de figure. Il vaut mieux en effet majorer les
effets en sécurité pyrotechnique.
Les formulations de Kinney (1985), du TM5-1300 (2008) et du NRDC (2001)
fournissent des évolutions de la hauteur du pied de Mach qui peuvent être différentes en
fonction de la position du point d’observation par rapport au centre de la charge. Cette thèse
étudiant la propagation d’une onde de choc issue de la détonation de deux natures de charge
explosive (gaz et TNT) et ne pouvant définitivement conclure sur la formule la plus
pertinente, on retiendra par la suite celle qui est la plus simple car ne reposant que sur des
paramètres géométriques, c'est-à-dire la formule de Kinney (1985).
1.3.1.8 Estimation de la surpression réfléchie
Le maximum de surpression réfléchi le long de la paroi plane est déterminé en
fonction des paramètres de la configuration (W, d, H et Ω) à partir de l’angle d’incidence ( )
et de la distance réduite ( ). La distance ( ) permet d’estimer le nombre de εach incident à
partir du maximum de surpression incidente (∆Pi+) issu de formules empiriques (TNT :
Kinney (1985) ; propane-oxygène à la stœchiométrie : Trélat (2006)). Lors de la propagation
de l’onde de choc, le régime de réflexion évolue le long de la paroi : RR (Réflexion
Régulière) ou RM (Réflexion de Mach). Le point de transition entre ces régimes de réflexion
peut être déterminé à partir de la relation simplifiée de Kinney (1985). Dans le cas d’une
approche analytique, le calcul des variables ( et ) varie en fonction du régime de réflexion,
figure 26.
Figure 26 : Diagramme des réflexions sur une paroi plane pour une onde incidente sphérique divergente
Dans le cas d’une réflexion de εach, le point à considérer est le point triple (point T) et non
le point d’impact sur la paroi (point I). δe point d’impact sur la surface (point I) est utilisé
dans le cas d’une réflexion régulière. Les abaques de réflexion (TM5-1300 (2008)) tiennent
compte de l’effet du pied de εach, les valeurs des variables et sont donc obtenues à partir
des relations géométriques pour le point d’impact I (figure 20).
35
M>1
min
δ’angle d’inclinaison étant petit ( <<), il est possible de réaliser un développement de Taylor
tan
limité à l’ordre 1 de la relation (80), équations (97) à (100).
2 1 M 12 sin 2
tan
1 M 2 sin 2
(97)
1
tan tan
cos 2
développement limité à l’ordre 1 ( ≈ 0) (98)
tan
2 1 M 12 sin 2 tan
cos 2 1 M 12 sin 2
(99)
M 12 sin 2 1
1 M 2 tan (100)
1
2
δ’angle d’incidence ( ) est proche de l’angle d’incidence minimal ( min), la relation (100)
peut se simplifier à partir de cet angle limite, équations de (101) à (103).
sin min
1
(101)
M1
M 12 sin 2 1
1 M 12
1
(103)
2 M 12
A partir des relations développées ci-dessous, il est possible de réécrire les équations de
variation des grandeurs thermodynamiques (P et u) présentées dans le cas d’un choc oblique
(relations (76) et (77)) pour un choc faible, équations (104) à (106).
P2 M 12
1
M 12 1
(104)
P1
1 g (M 12 ) 1 avec g M 2 1
u 2 1 g M 12 1 2
u1 M 12
M 12 1
1
(105)
2
1 g M 12
1
u1
g M 12
(106)
1
1
u2
Dans le cas d’un choc faible, l’angle ( ) étant petit, un développement limité à l’ordre 1 peut
être utilisé.
1 1
1 g M 12 1 g M 12 1 g M 12 1 ( 2 )
1 1
u1
(107)
u2
M 12
1
u1
M 12 1
(108)
u2
La variation de la norme de la vitesse (ω = || u ||) peut être déterminée pour un choc faible (
petit) à partir de la variation de la vitesse (u) et des relations géométriques sur les
composantes des vitesses développées lors de l’étude de l’angle de déviation, relation (109).
2
1
1 M 12 1
(109)
Les variations des grandeurs thermodynamiques présentées ci-dessus pour un choc faible
peuvent permettre de décrire un écoulement lors d’un phénomène de compression en
considérant une succession de segments faiblement inclinés les uns par rapport aux autres. Le
passage entre deux segments successifs correspond à un passage d’un état uniforme à un autre
avec de faibles variations sur les grandeurs thermodynamiques. Le long de ces segments, la
variation des grandeurs est nulle. Ces lignes séparées par un angle de déviation très faible
(d → 0, choc infinitésimal) sont appelées lignes de Mach. Les relations (110) et (111)
permettent de calculer la variation des grandeurs thermodynamiques au passage d’une ligne
de Mach par rapport au milieu amont du segment courant.
dP M 2 d
M 2 1
(110)
P
d d
M 2 1
(111)
37
Lignes de
M1 Mach
M2
Dans le cas d’une détente, l’angle d’inclinaison correspond à l’angle de déviation ( ). Cet
angle est compté positivement (à l’inverse d’une compression), équations (112) et (113).
d d
M 2 1
(112)
d
Soit d M 2 1
(113)
h0 h
u2
(114)
2
h c p T et M
u
R T
pour un gaz parfait (R = 8,3413 J.mol-1.K-1) (115)
1
1 M2
T0
(116)
T 2
δa variation de la norme de la vitesse peut s’exprimer en fonction du nombre de εach.
cM (117)
d
dc dM
(118)
c M
38
R T 1
1
1 M2
c2
R T0
(119)
2
c0 2
1 M dM
1 M 2 2
dc
(120)
c
d
M 1 M 2 M
dc dM 2 dM
(121)
2
c
M 2 1
d
1 M 2 M
dM (122)
1
2
Cette dernière relation permet d’obtenir la fonction de Prandtl-Meyer qui permet de calculer
le nombre de Mach en aval d’une détente à partir de l’angle de déviation ( ) et du nombre de
Mach en amont de la détente pour un choc droit, équations (123) et (124).
f M
M 2 1
1 M 2 M
dM (123)
1
2
f M
1
1
1
arctan
M 2 1 arctan M 2 1
1
(124)
f M 2 f M1
peut ensuite être déduit par dichotomie de la fonction de Prandtl-Meyer.
(125)
δors d’un phénomène de détente pour un écoulement supersonique isentropique d’un gaz
parfait, il est possible d’exprimer la variation de température et de pression, relations
1 M
(126) et (127).
1
2
1 M 2
T2 1
2 (126)
1
T1
2
2
1 M 2 1
1 1
Cd 2
1
P
P1 2
2 (127)
1 M2
2
Lors d’une détente, la vitesse de l’onde de choc est accélérée (ε1 < M2), la température et la
pression diminue (T1 > T2 et P1 > P2), le fluide se détend en aval.
39
Figure 29 : Evolution du coefficient d’atténuation en fonction du nombre de Mach (M1) en amont de la détente
M1 Détente
Mohanty
2
M2
Le phénomène de détente d’une onde de choc a été étudié par Mohanty et Ripley
(2004) (cité par Miller (2004)). A partir de données expérimentales obtenues lors de l’étude
40
de l’interaction d’une onde de choc avec une structure parallélépipédique, Mohanty (2004) a
établi un coefficient de détente en fonction du nombre de diffraction et de détente subi par
l’onde de choc autour d’un obstacle parallélépipédique (ξ = 90°) (Miller (2004)), relation
(128) et figure 30. Ce coefficient de détente (Cd) permet une estimation rapide de la
surpression en aval du phénomène de détente. Il ne prend néanmoins pas en compte l’angle de
déviation ( ), ni le nombre de εach incident (M1). Ce coefficient est donc établi pour une
configuration particulière.
W
Paramètres :
d
sol (en m)
d e : distance entre le centre de la charge et la face
e
avant (en m)
μ épaisseur au sommet de l’obstacle (en m)
H
α1
: hauteur du merlon (en m)
μ angle d’inclinaison de la face avant
α2 μ angle d’inclinaison de la face arrière
Charge (W)
HOB α1 α2
δes merlons sont des obstacles de protection placés sur le parcours de l’onde de choc
visant à atténuer ses effets (surpression, impulsion…). δes travaux de Allain (1λλζ) et de
Borgers (2010) proposent différentes observations et constatations pour ce type de structure
sans épaisseur au sommet (e = 0 m) et possédant une ou deux pentes inclinées à 45°
(figure 33).
y
Zone en vue directe Zone protégée
de la charge
Charge (0, 0) H
HOB = 0 m
α1
x
« Merlon » n°3 face avant à 45°, Borgers (2010)
Figure 32 : Représentation schématique des configurations étudiées par Borgers (2010) et Allain (1994)
A partir de la figure 32, les paramètres suivants permettent de décrire une configuration de
barrière de type proche du merlon :
d représente la distance entre le centre de la charge et le pied avant de l’obstacle,
H est la hauteur de la barrière de protection,
α désigne l’angle d’inclinaison d’une surface de l’obstacle par rapport au sol,
HOB identifie la hauteur du centre de la charge par rapport au sol.
La zone dite « directe » délimite la zone en vue de la charge, à l’inverse la zone dite
« protégée » identifie une zone où la barrière de protection fait obstacle à l’arrivée de l’onde.
Cette dernière est comprise dans la zone d’intérêt de cette étude.
δ’étude réalisée par Allain (1λλζ) se divise en deux étapes. δa première est constituée
d’une étude expérimentale du merlon n°1 à moyenne échelle (avec α1 = α2 = 45° et
H = 1,5 m) pour des charges de TNT (W = 8 et 37 kg, charge sphérique équivalente en champ
libre) et de composition B2 (η0 kg). δa distance entre la charge et l’obstacle (d) varie de 0,7η
à 17 m/kg1/3. Dans une seconde étape, les données expérimentales sont comparées à des
2
Composition B : explosif militaire composé de 59,5 % de RDX (Cyclonite) et de 40,5 % de TNT (paraffine
comme liant). Il possède une vitesse de détonation de 8050 m/s.
42
Les travaux de Allain (1994) présentent donc une étude paramétrique pour un merlon donné
(H = 1,5 m, e = 0 m et α1 = α2 = 45°). Les variations paramétriques sont effectuées sur la
distance entre la charge et l’obstacle (d) et de la masse de TNT (W).
Etude des écarts entre les résultats d’expériences pour une charge sphérique de 37 kg de TNT
et une référence bibliographique
Les premières configurations sont des barrières positionnées à une distance (d) de
13,3, 18,4, 33,32 et 56,65 m par rapport au centre de la charge. A partir des données en champ
libre de Allain (1994), il est possible de calculer les écarts entre les résultats expérimentaux et
la référence bibliographique sélectionnée, ici Kinney (1985). Les comparaisons sont
présentées sur la figure 33. A partir de ces graphiques, un intervalle maximal des écarts entre
la référence bibliographique (Kinney (1985)) et les données expérimentales de Allain (1994) a
pu être établi pour l’ensemble des essais réalisés, il est de - 10 à + 18% par rapport à
43
Kinney (1985), tableau 4. Les mesures de Allain (1994) restent proches des valeurs attendues
par Kinney (1985).
Numéro de configuration Numéro de figure Ecart maximal par rapport à Kinney (1985)
1 33-a - 10% + 2,7%
2 33-b - 2,5% + 5,2%
3 33-c - 10% + 18%
4 Pas d’information disponible
Tableau 4 : Intervalles des écarts entre la référence bibliographique et les résultats expérimentaux pour une
charge sphérique de 37 kg de TNT en champ libre
Sur la figure 33-c, l’évolution de l’écart entre les résultats du code PATRIC et la
référence bibliographique est présentée. δ’auteur ne précise pas les conditions physico-
numériques de calculs, ni les schémas numériques utilisés dans le code. Les résultats
numériques en champ libre permettent de constater une sous-estimation jusqu’à - 65% par
rapport à Kinney (1985). Les simulations numériques effectuées en présence d’une barrière de
protection seront donc sujettes à des écarts du même ordre de grandeur.
En dépit de ces écarts, les calculs réalisés avec le code PATRIC pour les configurations
expérimentées peuvent être une source de comparaison pour les évolutions de la surpression,
mais ils ne restituent pas l’expérience.
a)
b)
c)
Figure 33 : Analyse des écarts expérimentaux / bibliographiques / simulations avec le code PATRIC (Allain
(1994)) pour une charge sphérique de 37 kg de TNT
44
a)
b)
Figure 34 : Analyse des écarts expérimentaux / bibliographiques pour une charge sphérique de 8 kg de TNT
configurations avec une charge de γ7 kg de TNT et avec une hauteur d’obstacle de 1,η m (soit
0,45 m/kg1/3). Pour expliquer ce résultat, il est possible d’établir un historique des
phénomènes physiques apparus lors de l’interaction de l’onde de choc avec la structure.
a)
b)
c)
d)
Figure 35 : Evolution de la surpression après le passage de l’obstacle pour une charge de γ7 kg de TNT 3
3
Comme précédemment la simulation PATRIC de Allain (1994), configuration n°4, ne restitue pas l’expérience.
46
δe premier phénomène est une réflexion de l’onde de choc sur la face avant de
l’obstacle. A partir des formules empiriques, géométriques et analytiques associées à l’abaque
du TM5-1300 (2008), le régime de réflexion et la surpression réfléchie peuvent être
déterminés et évalués. Le tableau 6 synthétise une analyse réalisée sur trois points de mesure
placés sur la barrière de protection pour une charge sphérique de 37 kg de TNT.
Coordonnées des Surpression Surpression Nombre de
Régime de
points de mesure incidente Nombre Angle réfléchie Mach au
Numéro de Réflexion
(en m) (en bar) de Mach incident (en bar) niveau du
configuration (Kinney
(Kinney incident (en degré) (TM5-1300 pied de
X Y (1985))
(1985)) (2008)) Mach
13,55 0,25 0,448 1,180 45° RR 1,11 /
1 14,05 0,75 0,415 1,175 48° RM 1,11 1,40
14,55 1,25 0,385 1,153 50° RM 1,04 1,38
18,65 0,25 0,236 1,097 45° RR 0,53 /
2 19,15 0,75 0,226 1,093 45° RR 0,51 /
19,65 1,25 0,216 1,089 49° ≈Rε 0,51 1,20
33,57 0,25 0,098 1,041 45° RR 0,21 /
3 34,07 0,75 0,096 1,040 45° RR 0,20 /
34,57 1,25 0,094 1,039 47° RR 0,20 /
56,9 0,25 0,052 1,022 45° RR 0,104 /
4 57,4 0,75 0,051 1,022 45° RR 0,102 /
57,9 1,25 0,051 1,022 45° RR 0,102 /
RR : Réflexion Régulière ; RM : Réflexion de Mach
Tableau 6 : Réflexion sur la face avant de l’obstacle pour une charge sphérique de 37 kg de TNT
Il est possible d’émettre l’hypothèse que la diffraction et la détente sur la face arrière,
du pied Mach, provenant de la réflexion sur la face avant, ou de l’onde incidente conduisent à
conserver le nombre de Mach du front de choc (faible atténuation de la surpression) et à
imposer un angle d’incidence proche de ζη° pouvant approcher les conditions nécessaires à la
formation d’un nouveau pied de Mach sur la fin de la face arrière du merlon pour certaines
configurations. Le nombre de Mach du front de choc doit être supérieur ou proche de 1,3 pour
permettre la formation de ce pied de Mach. Cette éventuelle détente s’accompagne donc
d’une réflexion sur le sol, et peut conduire à la création d’un nouveau pied de Mach plus fort
sur la surpression et l’impulsion, figure 36. A ce stade, cette description n’est qu’une
hypothèse de travail car elle est basée sur l’approche de la détente de Prandtl-Meyer pour un
simple front de choc (§ 1.3.2.2).
47
A partir des essais de Allain (1994), il est possible de constater que le problème n’est
pas linéaire. En effet, l’observation des graphiques de la figure 35 permet de constater que
l’effet « amplificateur » du merlon n’évolue pas linéairement avec la distance entre la charge
et la structure.
Pour les configurations 1 et β, l’augmentation de la surpression maximale se situe à 2 H avec
un écart qui a priori décroît avec l’augmentation de la distance entre la charge et l’obstacle
(d). Cependant, pour les configurations 3 et 4, la surpression croît à l’arrière de l’obstacle
jusqu’à ζ H avec l’augmentation de la distance d. De plus, une singularité est visible sur la
configuration 3 avec une augmentation de la surpression par rapport au champ libre de
manière constante, sans montée brutale, commençant à la valeur du champ libre à 1 H. A
partir de cette distance, et sur une certaine plage de valeurs de la distance d, la détente derrière
le merlon peut être moins marquée du fait d’une instabilité des régimes de réflexion (RR et
RM) sur la face avant, voir tableau 6. δ’auteur ne précise pas la longueur de la barrière, un
effet de contournement ne peut donc pas être totalement exclu à partir d’une certaine
distance d.
Etude paramétrique pour une charge sphérique de 8 kg de TNT
Pour ces essais, Allain (1994) fait varier non seulement la distance par rapport à la
charge, mais il ajoute également un second merlon sur la configuration 6, voir figure 37.
δ’objectif est de vérifier dans quelle mesure une deuxième barricade pourrait augmenter
l’effet protecteur du merlon. Cependant, au regard des résultats précédents, il est possible que
ces configurations puissent produire un effet inverse et finalement accentuer les
caractéristiques de l’onde de choc même après ce deuxième obstacle.
y 6,5 m
2,8 m
1,5 m 1,5 m
Charge (0, 0)
8 kg de TNT
1H
δ’ensemble des résultats pour ces configurations est présenté sur la figure 38. En
observant les données expérimentales obtenues, une augmentation de la surpression est
observée, comme lors de l’analyse précédente. Pour la configuration η, l’effet
d’ « augmentation » de la surpression est visible, mais il tend à se réduire et permet de
48
retrouver les données du champ libre vers 8 H. Cette observation peut s’expliquer avec les
mêmes hypothèses que pour la charge sphérique de 37 kg de TNT, voir tableau 7.
Coordonnées Surpression Nombre Surpression Nombre
Régime de
des points de incidente de Mach Angle réfléchie de Mach
Numéro de Réflexion
mesure (en bar) incident incident (en bar) au niveau
configuration (Kinney
(en m) (Kinney (Kinney (en degré) (TM5-1300 du pied
(1985))
X Y (1985)) (1985)) (2008)) de Mach
1,75 0,25 12,93 3,48 54 RM 47,16 6,43
5 2,25 0,75 6,85 2,62 63 RM 12,90 3,47
2,75 1,25 3,95 2,10 70 RM 5,58 2,40
3,05 0,25 3,83 2,07 50 RM 9,46 3,02
6 3,55 0,75 2,57 1,79 57 RM 4,54 2,21
4,05 1,25 1,79 1,59 62 RM 2,95 1,88
Tableau 7 : Réflexion sur la face avant de l’obstacle pour une charge sphérique de 8 kg de TNT
a) merlon unique
b) merlon double
Figure 38 : Evolution de la surpression après le passage de l’obstacle pour une charge 8 kg de TNT
Cette hypothèse et l’ensemble de celles évoquées lors de cette analyse feront l’objet
d’une comparaison et d’une validation par des calculs CFD avec le code HERA. Le tableau 8
synthétise les effets constatés pour chaque configuration étudiée par Allain (1994). Les
configurations étudiées par l’auteur avec le logiciel PATRIC ne sont pas abordées. En effet,
les simulations numériques menées par l’auteur ne restituent pas assez bien les expériences
pour apporter des informations supplémentaires à la problématique étudiée.
Distance au sol par rapport au pied de la face avant (en m) Distance au sol par rapport au pied de la face avant (en m)
Figure 39 : Evolution des caractéristiques de l’onde de choc en présence d’un obstacle (Borgers β010)
H = 4,3 m ; W = 230 kg ; HOB = 1,58 m ; Merlon n°1 H = 4,3 m ; W = 230 kg ; HOB = 1,58 m ; Merlon n°1
y (n x la hauteur du merlon)
y (n x la hauteur du merlon)
Distance au sol (en n x la hauteur du merlon) Distance au sol (en n x la hauteur du merlon)
H = 4,3 m ; W = 230 kg ; HOB = 1,58 m ; Merlon n°2
H = 4,3 m ; W = 230 kg ; HOB = 1,58 m ; Merlon n°2
y (n x la hauteur du merlon)
y (n x la hauteur du merlon)
Distance au sol (en n x la hauteur du merlon) Distance au sol (en n x la hauteur du merlon)
H = 4,3 m ; W = 230 kg ; HOB = 1,58 m ; Merlon n°3 H = 4,3 m ; W = 230 kg ; HOB = 1,58 m ; Merlon n°3
y (n x la hauteur du merlon)
y (n x la hauteur du merlon)
Distance au sol (en n x la hauteur du merlon) Distance au sol (en n x la hauteur du merlon)
a)surpression positive b)impulsion positive
Figure 40 : Maximum de surpression et d’impulsion positive en aval des merlons de Borgers (2010) (% par
rapport au champ libre)
Il est essentiel de traiter séparément les deux hauteurs de HOB. Une charge possédant
une HOB ≠ 0 m, tout en étant proche du sol, implique généralement la création d’un pied de
Mach par une réflexion au sol. Les deux cas ne possèdent donc pas le même type de réflexion
sur la face avant de l’obstacle. δa méthode d’analyse pour une charge au sol reste donc
identique au cas de Allain (1994). En revanche, pour une charge en hauteur par rapport au sol,
la formule (95) de Kinney (1985) permet d’estimer la hauteur du pied de Mach ainsi formé.
La HOB est égale à 1,58 m pour une masse de TNT de 230 kg. A partir de ces données
d’entrée et en tenant compte de l’hypothèse que la réflexion de εach apparaît sur le sol pour
un angle d’incidence supérieur à 40° (cas pour un nombre de Mach supérieur à 1,6), il est
possible d’estimer la hauteur du pied de εach arrivant sur l’obstacle (voir figure 23),
tableau 9.
Hauteur du pied de Mach,
Numéro du merlon étudié par Distance entre le centre de la charge et
hm (en m)
Borgers (2010) l’obstacle, d (en m)
(Kinney (1985))
1–3 4,3 0,56
2 6,45 1,65
Tableau 9 : Hauteur du pied de Mach suivant les configurations de Borgers (2010)
53
D’après les résultats du tableau 9, la hauteur du pied de Mach (hm) est inférieure à la
hauteur de l’obstacle (hm < H). La formation de ce pied de Mach sur le sol devrait induire un
nombre de Mach incident plus élevé, suite à la réflexion sur le sol, qu’une onde incidente
sphérique (exemple pour une charge au sol). Cette phénoménologie place donc le pied de
εach dans le cas d’un nombre de εach supérieur à γ (minimum du nombre εach observé
pour une onde frappant la paroi avec une charge au sol). Pour les merlons n°1 et γ, l’angle
d’incidence entre l’onde de choc et la structure est supérieur ou égal à 45°. Il est donc difficile
d’estimer avec certitude la nature de la réflexion. Pour la structure possédant une face avant à
90°, le pied de εach subit une réflexion normale à la paroi, seule l’onde incidente provenant
directement de la charge passe par-dessus l’obstacle.
Pour une HOB = 0 m, la bibliographie est plus fournie et donc la compréhension de
l’interaction de l’onde de choc avec la structure est plus simple. Le tableau 10 ci-dessous
résume les résultats, obtenus par les formules analytiques et empiriques, de la réflexion sur la
face avant de l’obstacle suivant la configuration étudiée. δes merlons n°1 et γ étudiés pour
une charge au sol conduisent à la formation d’un premier pied de εach lors de la réflexion
sur la face avant. Concernant le merlon n°2, la réflexion sera régulière avec un angle
d’incidence inférieur à ζ0° (incidence maximale de = γη°).
Coordonnées des Surpression Surpression Nombre de
Régime de
points de mesure incidente Nombre Angle réfléchie Mach au
Numéro de Réflexion
(en m) (en bar) de Mach incident (en bar) niveau du
configuration (Kinney
(Kinney incident (en degré) (TM5-1300 pied de
X Y (1985))
(1985)) (2008)) Mach
4,5 0,2 18,98 4,16 47,5 RM 109,41 9,74
Merlon 1 6,3 2 8,44 2,87 62,6 RM 18,87 4,14
(d = 4,3 m)
8,5 4,2 3,73 2,05 71,3 RM 5,27 2,35
6,45 0,2 8,90 2,94 2 RR 41,17 /
Merlon 2
6,45 2 8,06 2,81 18 RR 38,88 /
(d = 6,45 m)
6,45 4,2 6,00 2,48 33 RR 24,00 /
4,5 0,2 18,98 4,16 47,5 RM 109,41 9,74
Merlon 3
6,3 2 8,44 2,87 62,6 RM 18,87 4,14
(d = 4,3 m)
8,5 4,2 3,73 2,05 71,3 RM 5,27 2,35
Tableau 10 : Réflexion sur la face avant de l’obstacle pour une charge au sol de βγ0 kg de TNT
rapidement dès le pied de la face arrière ce qui peut indiquer la présence d’une réflexion de
type Mach sur le sol.
Dans le cas du merlon n°β, l’onde qui va subir une diffraction et une détente est l’onde
incidente qui passe par-dessus l’obstacle. Cette onde va donc, en accord avec les hypothèses
émises précédemment, se détendre ce qui induit une atténuation du maximum de surpression
de l’onde. Par la suite, l’onde va se réfléchir sur le sol. δ’observation de la figure 39-b pour le
merlon n°β montre que l’impulsion n’augmente rapidement qu’à partir de β H après le
merlon. A cette distance l’angle incident est de 71°, soit vis-à-vis d’une éventuelle formation
d’un pied de Mach, un nombre de Mach faible et proche de 1.
Cette observation rappelle que les phénomènes de diffraction et de détente sont des
phénomènes physiques complexes et non-linéaires. δ’impact de ces phénomènes sur
l’évolution du maximum de surpression dépend de l’intensité de l’onde qui se détend en aval
et de l’angle de déviation de l’écoulement : l’onde incidente possède un nombre de εach plus
faible qu’un pied de Mach issu d’une réflexion de εach.
Concernant le merlon n°γ, la différence avec le merlon n°1 se situe sur l’inclinaison de
sa face arrière, ici à λ0°. δ’onde qui va se détendre est également un pied de Mach dont
l’origine se situe à la réflexion de l’onde sur la face avant de l’obstacle. Cependant, la
trajectoire de la détente n’est pas « guidée »4 pour le merlon n°γ (α = λ0° sur la face arrière) et
la réflexion sur le sol en aval de l’obstacle commence donc dans un premier temps par une
réflexion régulière ( < ζ0°). Puis, dès que l’angle d’incidence est supérieur ou égale à ζ0°, la
réflexion devient une réflexion de type Mach. Ce phénomène est visible sur la figure 39-b
pour le merlon n°γ avec une augmentation de l’impulsion juste avant 1 H après l’obstacle, soit
avant ≈ 45°.
4
Détente « guidée » μ détente d’une onde de choc sur une paroi inclinée, l’obligeant à se détendre suivant une
trajectoire bien définie.
55
HOB = 0 m, - 40%
protecteur en champ proche du
HOB = 0 m, - 45%
protecteur. Pour cette topologie
de merlon l’effet protecteur dure
HOB = 1,58 m, - 60%
8,6 m
2 230
(soit 1,40 m/kg1/3) beaucoup plus longtemps, avec
25% d’atténuation à γ H
encore :
formule (83), ainsi que des caractéristiques mécaniques (εi) et géométriques ( = f(HOB)) de
l’onde incidente (§ 1.3.1.7).
Si cette distance (dc) est inférieure à la distance entre le centre de la charge et le pied
amont du merlon (dc < d), un pied de Mach se propage sur le sol. La hauteur de ce pied
de Mach, au niveau de la face avant de la barrière de protection, peut être estimée à partir
de la formule (95). Le régime de réflexion sur la face avant (réflexion de Mach direct,
inverse ou stationnaire,…) dépend alors de la hauteur (hm) et de l’intensité du pied de
Mach (ΔP+, Mi), ainsi que des dimensions de la paroi (d, α1, H).
Si la distance d’apparition du pied de εach est supérieure à la distance entre le centre de
la charge et la barrière (d < dc), les conditions nécessaires à l’apparition d’un pied de
Mach ne sont pas réunies. δ’onde incidente issue de la charge se réfléchit sur la face
avant du merlon.
Dans le cas d’une charge placée au niveau du sol, l’onde incidente se réfléchit sur la paroi
avant du merlon (HOB = 0 m ou si HOB ≠ 0 m et d < dc). Le régime de réflexion (réflexion
régulière (RR) ou réflexion de Mach (RM)) dépend des caractéristiques mécaniques (Mi) et
géométriques ( = f(d, α1, H)) de l’onde incidente.
δe régime de réflexion au sommet du merlon permet d’identifier la nature et l’intensité de
l’onde qui se propage en aval de la barrière de protection : pied de Mach (RM) ou onde
incidente (RR).
Sur le sommet de la barrière de protection, l’onde incidente (RR) ou le pied de Mach
(RM) subit une ou deux détentes en fonction du dimensionnement de l’épaisseur (e). En effet,
la présence d’une épaisseur non nulle (e ≠ 0 m) sur le sommet de la barrière permet de réaliser
une détente dite « en deux étapes » μ première détente sur l’arête de la face avant et une
seconde sur l’arête de la face arrière. δes recommandations de l’OTAN (GBPP (β00λ))
évaluent une épaisseur minimale de 0,η m (e > 0,η m). Si l’épaisseur est nulle (e ≈ 0 m),
l’onde subit une seule détente sur l’arête entre les faces avant et arrière de la barrière de
protection.
En aval de la barrière de protection, l’onde issue du sommet du merlon se réfléchit sur
le sol. Le régime de réflexion en aval du merlon (RR ou RM) peut être identifié et la
surpression réfléchie estimée à partir des caractéristiques mécaniques (εs, ΔPs) et
géométriques ( = f(α2, H)) de l’onde au niveau du sol.
Il est important de garder à l’esprit que l’interaction d’une onde de choc avec une
structure, atténuation ou augmentation de la surpression ou de l’impulsion, dépend de la
configuration de l’obstacle et de la position de la charge. Ces paramètres condionnent alors le
régime d’écoulement incident qui parvient sur la structure (RR, Rε, détente,…).
57
Création d’un pied de Mach en HOB est la hauteur du centre de la charge par rapport au sol. Réflexion de l’onde incidente
amont de l’obstacle par une Mi représente le nombre de Mach incident. directement issue de la charge sur la
réflexion sur le sol face avant de l’obstacle
ΔP indique la valeur du maximum de surpression incidente.
α est l’angle d’inclinaison de la surface de l’obstacle par
rapport au sol. Suivant εi, ΔP et α
d est la distance entre la charge et l’obstacle.
Réflexion de ce pied de Mach sur la
face avant de l’obstacle
dc est la distance d’apparition du pied de Mach.
Suivant hm, ΔP, α1 H, …
Formation d’un pied Formation d’un pied Formation d’un pied Formation d’un pied Réflexion régulière
de Mach direct de Mach stationnaire de Mach inverse de Mach directe
Détente du front
e>Xm e < X m (cas de e = 0 m)
Phénoménologies au sommet
l’obstacle.
Figure 41 : Diagramme des phénomènes physiques présents lors de l’interaction d’une onde de choc avec une barrière de protection
58
Mur
H
Charge
HOB
Il existe aujourd’hui différents types de compositions pour les murs de protection. Les
premiers murs étaient réalisés avec des sacs de sable. Ils possèdent l’avantage de ne produire
aucun fragment secondaire susceptible de tuer ou de blesser. Les murs en béton ou en acier
sont cependant préférés pour des installations permanentes. En effet, ils résistent mieux aux
épreuves du temps, même s’ils sont plus susceptibles de produire des fragments secondaires
lors de leur endommagement éventuel. Notons que ces dernières années, la recherche sur ces
structures s’oriente vers des barrières composites, alliant le métal et le béton pour un coût
raisonnable et une résistance satisfaisante. Dans ce qui suit, les cas analysés se restreindront à
des murs infiniment rigides, indéformables et donc sans production d’éclats.
asymptotique sur la distance maximale de protection, elle atteint 6 H pour des charges
comprises entre 30 et 100 kg de TNT.
Charge
W = 3 kg
Charge
W = 10 kg
W = 30 kg
W = 100 kg Hauteur de 0,05 m
W = 300 kg W = 1 kg
W = 3 kg
W = 10 kg
W = 30 kg
W = 100 kg
W = 300 kg
Distance au sol par rapport au pied aval (en m) Distance au sol par rapport au pied aval (en m)
Figure 44 μ Evolution de l’effet protecteur suivant la valeur de la masse de la charge W (Borgers (β010)),
HOB = 1,58 m, H = 4,3 m
Pour expliquer ces résultats, une étude sommaire des phénoménologies est présentée
ci-dessous. La méthode et les formules utilisées sont identiques à celles de l’étude
bibliographique réalisée précédemment pour les configurations de type merlon.
δ’étude bibliographique se décompose en deux étapes. δa première consiste à déterminer la
hauteur du pied de Mach (hm) formé par une réflexion sur le sol en amont du mur de
protection. Si la hauteur du mur est supérieure à la hauteur du pied, alors les caractéristiques
de l’onde de choc incidente en haut du mur sont déterminées. δ’épaisseur du mur étant
simplement considérée comme fine par l’auteur (aucune indication de valeur n’est donnée par
l’auteur, e ≈ 0), la diffraction et la détente de cette onde ne se feront qu’après la barrière de
protection.
Pour cette configuration, la hauteur du pied de Mach est négligeable par rapport à la hauteur
du mur. δ’onde subira donc une réflexion régulière et une réflexion de Mach sur la face avant
de l’obstacle. δ’onde qui subira la détente en aval du mur est donc le pied de Mach issu de
61
cette réflexion. Les valeurs sont estimées à partir des formules analytiques et empiriques de
Kinney (1985), voir tableau 13.
Distance entre le centre de la charge et le sommet : 3,38 m
Angle incident de l’onde au sommet de la face avant : 54° (nombre de Mach de transition de 1,12)
Surpression
Surpression incidente Nombre de Mach
Charge étudiée W Nombre de Mach réfléchie au
(en bar) au niveau du pied
(en kg de TNT) incident sommet (en bar)
(Kinney (1985)) de Mach
TM5-1300 (2008)
1 0,63 1,24 1,59 1,54
3 1,42 1,49 3,08 1,91
10 3,61 2,02 7,44 2,72
30 8,29 2,85 26,41 4,86
100 19,31 4,19 78,58 8,27
300 38,83 5,86 195,02 12,97
Tableau 13 : Onde incidente au sommet du mur en fonction de la charge W
Charge
Charge
Hauteur de 0,05 m
H = 0,5 m
H = 1,0 m
H = 1,5 m Hauteur de 0,05 m
H = 2,0 m H = 0,5 m
H = 2,5 m H = 1,0 m
H = 3,0 m H = 1,5 m
H = 2,0 m
H = 3,5 m
H = 2,5 m
H = 4,0 m
H = 4,5 m ect..
H = 5,0 m
Distance au sol par rapport au pied aval (en m) Distance au sol par rapport au pied aval (en m)
Figure 45 : Evolution de l’effet protecteur en fonction de la hauteur du mur (Borgers (β010)), HOB = 1,58 m,
W = 230 kg de TNT
62
Pour comprendre cette observation, une analyse est effectuée selon le même principe
que l’étude précédente (variation de W). La hauteur du pied de Mach sur le sol et une partie
des caractéristiques de l’onde incidente sur le sommet de la barrière de protection sont ainsi
estimées à partir de formules issues de la littérature.
Les résultats sont présentés dans le tableau 14 ci-dessous. Les deux premières réflexions
calculées (H = 0,5 et 1 m) après le mur indiquent la présence d’un pied de Mach en aval du
mur pour ces configurations.
Les valeurs du nombre de Mach au sommet du mur permettent de conclure que même si la
détente derrière le mur permet d’atténuer le maximum de surpression, un pied de Mach se
formera après l’obstacle pour un angle d’incidence égal à 40° compte tenu des « forts »
niveaux de surpression estimés. Ces distances correspondent, sur la figure 45, à une
augmentation brusque de la surpression (x = 1 m). δ’intensité de cette augmentation (due au
phénomène de réflexion) dépend ensuite de la valeur du nombre de Mach au sommet du mur.
Si l’onde doit parcourir une distance plus importante (β,0ζ m pour H = β m et γ,λθ m pour
H = 5 m), le nombre de Mach au sommet sera pour sa part plus faible (M i = 8,23 pour
H = 2 m et Mi = 4,68 pour H = 5 m) ce qui peut retarder la possibilité de formation d’un pied
de Mach après le mur (effet protecteur).
Hauteur du pied de Mach en amont du mur (hm) : 2,8.10-2 m
Si H << HOB, point d’observation placé après le mur de protection, mais si H > HOB alors le point d’observation
devient le sommet du mur.
Surpression Distance
Hauteur
Surpression incidente δ’angle réfléchie après le approximative de
du mur Nombre de
Kinney (1985) (en bar) incident mur création du pied de
étudiée H Mach incident
/ Point d’observation (en degré) (en bar) Mach après le mur
(en m)
(TM5-1300 (2008)) (en m)
32,82 / 1 m après le mur 5,4 62° 88,22 /
0,5
68,74 / sommet du mur 7,54 / / 0,42
12,48 / 3,5 m après le 3,42 /
72° 21,70
1 mur
75,54 / sommet du mur 8,10 / / 0,84
1,5 80,31 / sommet du mur 8,36 1,26
2 77,77 / sommet du mur 8,23 1,68
2,5 68,99 / sommet du mur 7,75 2,10
H > HOB, pas de réflexion derrière
3 58,12 / sommet du mur 7,13 2,52
le mur en vue directe du centre de la
3,5 47,59 / sommet du mur 6,46 2,94
charge
4 38,13 / sommet du mur 5,80 3,36
4,5 30,37 / sommet du mur 5,20 3,78
5 24,42 / sommet du mur 4,68 4,20
Tableau 14 : Etude des configurations en fonction de la hauteur H du mur
Distance au sol par rapport au pied aval (en m) Distance au sol par rapport au pied aval (en m)
Figure 46 : Evolution de l’effet protecteur en fonction de la HOB (Borgers (β010)), W = 230 kg de TNT,
H = 4,3 m
δ’analyse bibliographique est résumée dans le tableau 15. Il présente, suivant la valeur de la
HOB, la hauteur du pied de Mach (hm) et le nombre de εach de l’onde incidente au sommet
du mur.
Variation de la hauteur du pied de Mach avec la variation de la HOB
Hauteur du mur fixée à 4,3 m pour une charge 230 kg de TNT
La distance entre la charge et la barrière de protection (d) est fixée à 2 m
Hauteur de la Hauteur du pied de
Distance Surpression incidente
charge par Mach hm
critique dc (en bar) Nombre de Mach incident
rapport au sol (en m)
(en m) (Kinney (1985))
(en m) (Kinney (1985))
0 / / 17,07 3,95
0,5 0,42 < d 0,500 20,90 4,35
1 0,84 < d 0,133 25,70 4,80
1,5 1,26 < d 0,036 31,96 5,33
2 1,68 < d 0,005 39,83 5,93
2,5 2,09 > d 49,66 6,60
3 2,51 > d Pas de pied de 60,56 7,27
3,5 2,94 > d Mach 71,85 7,91
4 3,36 > d 79,19 8,30
Tableau 15 : Etude des configurations en fonction de la HOB
Pour une HOB comprise entre 2,5 et 4 m, la distance critique de formation d’un pied de Mach
en amont de la barrière est plus grande que la distance entre le centre de la charge et le mur
(dc > d). Pour cette gamme de hauteurs de la charge, il n’y a donc pas formation d’un pied de
Mach sur le sol avant le mur. Et inversement, pour une HOB inférieure à 2,5 m, un pied de
εach se forme et une partie de l’énergie sera réfléchie sur le mur par celui-ci.
δa remarque de Borgers (β010) sur la plage optimale d’effet protecteur (0,η < HOB/H < 1)
peut aussi simplement s’expliquer par la variation de la distance parcourue par l’onde pour
passer au-dessus de l’obstacle. En effet, la distance à parcourir par l’onde avec une
HOB = 0 m est 2,3 fois plus élevée que pour une charge à 4,3 m du sol. Cet écart réduit bien
évidemment les caractéristiques mécaniques de l’onde de choc arrivant au sommet du mur
avec l’effet d’atténuation par la distance parcourue. Ceci est notamment visible pour le
nombre de Mach incident au sommet avec une valeur deux fois plus faible pour une charge au
sol par rapport à une charge à 4,3 m du sol.
Variation de la distance entre le centre de la charge et le mur (H = 4,3 m ; HOB = 1,58 m ;
W = 230 kg de TNT)
δes résultats obtenus indiquent que la distance maximale d’effet est quasiment
identique pour l’ensemble de la gamme des distances d étudiées avec un effet jusqu’à η - 6
64
Distance au sol par rapport au pied aval (en m) Distance au sol par rapport au pied aval (en m)
Figure 47 : Evolution de l’effet protecteur en fonction de la distance d entre le centre de la charge et le mur
(Borgers 2010), HOB = 1,58 m, W = 230 kg de TNT, H = 4,3 m
δ’étude de ces résultats montre deux types de phénoménologies, voir tableau 16. La
première phénoménologie est visible pour une variation de la distance d comprise entre 2 et
7 m. Pour ces distances, la hauteur du pied de Mach formé lors de la réflexion sur le sol reste
inférieure à la hauteur du mur. Cette observation permet de conclure que l’onde qui arrive au
sommet du mur est l’onde incidente directement issue de la charge.
Dans ces variations de d, la singularité à d = 7 m peut potentiellement s’expliquer par la
formation d’un pied de εach avec une hauteur plus grande que les autres. Une plus grande
partie de l’énergie sera réfléchie sur le mur par rapport aux autres configurations (d < 7 m).
La seconde observation concerne les distances d supérieures ou égales à 10 m. Pour ces
configurations, un pied de Mach se forme également en amont du mur. Cependant, la hauteur
du pied de Mach ainsi formé est supérieure à la hauteur du mur, ce qui entraîne que c’est un
pied de Mach et non une onde incidente issue de la charge qui parvient au sommet. Ce régime
d’écoulement ne va générer un effet protecteur que sur de courtes distances avant que le pied
de Mach ne puisse potentiellement se reformer « correctement ».
65
Les observations de Borgers (2010) sont résumées dans le tableau 17. Les travaux de
Borgers (β010) confirment l’intérêt d’étudier les variations des paramètres suivants : d, e, W,
H, α1 et α2.
résolutions de maillage trop grossières pour assurer, avec une précision suffisante, la
restitution des valeurs de crête du maximum de surpression. δa mise en œuvre d’un maillage
suffisamment fin est donc nécessaire afin de se placer dans des conditions de convergence
numérique satisfaisantes. De ce fait, dans les chapitres deux et cinq, un standard de calcul
(gaz et TNT) est mise en place afin d’assurer le calcul des valeurs de crête du maximum de
surpression.
Lors de ces travaux de thèse, la propagation d’une onde de choc issue de la détonation
de deux natures de charge explosive est étudiée : une charge gazeuse (propane-oxygène à la
stœchiométrie) et une charge chimique condensée (TNT, Trinitrotoluène). Ces deux types
d’explosif se différencient, notamment, par leur composition et leur célérité de détonation
(DCJ). Ces deux charges explosives conduisent donc à l’établissement d’ondes de souffles
aériennes différentes.
Tous les problèmes liés à une éventuelle équivalence entre les natures de charge gazeuse et
charge chimique condensée seront évités en simulant de manière différenciée ces deux natures
de sources explosives. Ainsi la détonation d’une charge gazeuse est directement simulée à
partir de ses caractéristiques de détonation suivant le modèle de Chapman-Jouguet et de
l’équation d’état des produits de détonation. Cette équation d’état est obtenue à partir d’un
code de thermochimie du CEA – DAM.
67
Charge (W) H
α1 α2
Figure 48 μ Représentation schématique d’une barrière de protection pour une charge placée au niveau du sol
68
Analyse bibliographique
(chapitre 1)
Etude de configurations
élémentaires (chapitre 3)
Algorithme
simplifié
(chapitre 4)
Etude paramétrique en
géométrie 2D-axisymétrique
δa détonation d’une charge de gaz a une célérité de détonation bien plus faible,
DCJ ≈ β000-3000 m.s-1, que celle d’une charge chimique, DCJ ≈ 7000-9000 m.s-1. Les
pressions engendrées par la détonation sont également bien plus faibles pour une charge
gazeuse, PCJ ≈ γ0-40 bar, que pour une charge chimique, PCJ ≈ 1η0-400 kbar.
δ’énergie libérée par la combustion d’une charge gazeuse peut être estimée à partir des
enthalpies de formation des produits (dioxyde de carbone (CO2) et l’eau (H2O)) et de
l’enthalpie de réaction du propane (C3H8), équation (130) (enthalpies issues du Hanbook of
chemistry and Physics, 78th Edition).
H r* 3 H CO
*
2
4 H H* 2O H C*3H8 5 H O* 2 2,044 MJ (130)
δ’énergie libérée par la réaction d’oxydoréduction du propane peut être ramenée par unité de
volume de la charge (propane-oxygène) aux conditions normales de température et de
pression, équation (131).
H r*
E gaz 15,2 MJ .m 3
n produits Vmol
(131)
Les charges utilisées sont de forme hémisphérique et positionnées sur le sol. La charge de gaz
est confinée dans une bulle de savon (combinaison d’eau savonneuse et de glycérine) dont la
capacité de confinement est suffisante pour assurer la forme hémisphérique et négligeable
devant l’énergie de l’explosion. Pour réaliser la bulle de savon, un cercle métallique est centré
par rapport à la source d’amorçage (centre de la charge explosive). δe gaz sous pression dans
la bouteille de pré-mélange est ensuite injecté dans le mélange aqueux, permettant ainsi de
former la charge gazeuse. Lors des différents essais, deux rayons de charges sont utilisés :
R1 = 0,06 m et R2 = 0,0γ m. δ’énergie libérée par la détonation de ces charges gazeuses est
présentée dans le tableau 18.
Charge de gaz : propane-oxygène à la stœchiométrie (C3H8 + 5 O2)
Rayon de la charge explosive (en m) R1 = 0,06 R2 = 0,03
Energie libérée par la détonation d’une sphère -3
E1 = 13,75.10 E2 = 1,72.10-3
de gaz (en MJ)
Tableau 18 μ Energie libérée par la détonation d’une charge sphérique de gaz pour les rayons (R1) et (R2)
Charge Amplificateur
gazeuse
Amplificateur
Alimentation
haute tension
δ’amorçage de la charge gazeuse est réalisé par le biais d’un fil explosé de cuivre de
diamètre 0,12 mm soudé entre deux électrodes, distantes de 8 mm, reliées à une alimentation
haute-tension. Lors du passage de la décharge électrique, le fil explosé est vaporisé ce qui crée
un point « chaud » dans la charge explosive par la formation d’un plasma métallique. Son
expansion dans le mélange crée alors une onde de choc qui amorce la détonation du mélange
gazeux, figure 51. Il est difficile d’établir précisément l’énergie d’amorçage transmise par
cette onde. Il est néanmoins possible de calculer l’énergie électrique d’amorçage nécessaire à
ce phénomène (Enom, relation (132)).
Energie électrique nominale : Enom Q V 2
1
(132)
2
Dans le cadre des campagnes de tirs, une décharge électrique de 7 kV est utilisée, soit
1λλ,ζγ J déchargé en 1 µs. δ’énergie d’amorçage ainsi estimée est alors très faible en
72
Figure 52 : Représentation schématique de la distance en vue directe de la charge pour des exemples de positions
de capteurs
4
Remarque valable pour la valeur de crête, la décroissance de la phase positive peut néanmoins être perturbée par l’onde
réfléchie sur la face avant de la structure Il est donc uniquement possible de s’assurer de la reproductibilité sur les
caractéristiques mécaniques du maximum de surpression (∆P+) et du temps d’arrivée (t0), mais pas forcément sur la durée de
phase positive (dt+), ni de l’impulsion positive (I+).
73
Les capteurs utilisés sont des capteurs piézoélectriques IPC 113B26 (annexe 2) connectés à
plusieurs amplificateurs (amplificateur 482C05, annexe 3) dont le signal de sortie est relié à
une centrale d’acquisition avec une fréquence de 5 Mhz sur une période de 8 ms. Le
déclenchement de l’enregistrement est donné par l’amorçage de la charge explosive à l’aide
d’un trigger qui relie l’alimentation haute tension et la centrale d’acquisition.
5
Cette méthode, qui correspond à une approche 1D, permet de déterminer la position du front de choc dans la charge
explosive à partir de sa célérité de détonation suivant le modèle de Chapman-Jouguet. Ce front délimite alors les états dits de
gaz frais (inerte) et de gaz brûlés (produits de détonation).
75
1 1 P V0
eP, A1 V0 e r1 A2 V0 e r2
r1 r2
(134)
Les coefficients de la JWL et les caractéristiques de détonation CJ sont obtenus à partir d’un
code de thermochimie du CEA – DAM. Le tableau 19 synthétise les résultats issus du code de
thermochimie pour l’équation d’état d’une charge de TNT.
Charge explosive de TNT (C7H5N3O6)
A1 A2 V0 e0 DCJ PCJ
r1 r2 ω
(Pa) (Pa) (m3/kg) (J/kg) (m/s) (Pa)
3,54283.1011 7,7732.109 4,20463 1,37927 0,225042 6,06061.10-4 5,08515.106 6850,8 2,1315.1010
Tableau 19 : Coefficients de l’équation d’état de type JWδ et caractéristiques de détonation d’une charge
explosive de TNT
Le tableau 20 présente les coefficients de l’équation d’état des produits de détonation dans le
cas d’une charge de gaz (propane-oxygène à la stœchiométrie).
Charge gazeuse (C3H8 + 5 O2)
A1 A2 V0 e0 DCJ PCJ
r1 r2 ω
(Pa) (Pa) (m3/kg) (J/kg) (m/s) (Pa)
3,54283.1011 7,7732.109 4,20463 1,37927 0,225042 0,72372 5,08515.106 2357 34,58.105 1,22
Tableau 20 : Coefficients de l’équation d’état de type JWδ et caractéristiques de détonation d’une charge
gazeuse (propane-oxygène à la stœchiométrie)
76
Il est important de noter que l’équation d’état des gaz non brûlés (en amont du front de
détonation) est différente entre une charge explosive de TNT et une charge gazeuse. Dans le
cas d’une charge de gaz, l’équation des gaz inertes (non brûlés) est l’équation d’état d’un gaz
parfait ( = 1,βλ), alors que pour une charge de TNT, la même équation d’état est utilisée pour
les gaz inertes et pour les produits de détonation.
u
.u u P 0
t
(136)
e
.e P u 0
t
(137)
Ces équations sont complétées par l’équation d’état du fluide étudié. Les produits de
détonation sont modélisés par l’équation d’état de type JWδ et l’air par l’équation d’état d’un
gaz parfait. δ’ensemble de ces équations forme un système de 4 équations à 4 inconnues.
Dans le cas d’une onde de choc aérienne issue de la détonation d’une charge explosive
(chimique ou gazeuse), les conditions des équations d’Euler sont remplies car les forces de
pression sont bien supérieures aux forces dues au comportement visco-élastique.
Phase positive
Maille de Le front de choc
P(t) > P0
niveau 0 occupe moins de 50 %
du volume
Forme « théorique » du
front de choc
δa première zone d’intérêt correspond à la phase positive de l’onde de choc (P(t) ≥ P0,
P0 pression initiale). Lors du passage du front de choc, la maille (de niveau 0) est découpée en
mailles plus petites, figure 54. La phase positive terminée (P(t) < P0), les mailles raffinées sont
regroupées jusqu’à former à nouveau la maille de taille initiale (niveau 0). Cette utilisation de
l’AεR est fonction du gradient de pression et d’une condition sur l’évolution de la phase
positive (P(t) > P0) pilotée par l’utilisateur.
2éme utilisation de l’AεR :
La fonctionnalité AMR du code HERA permet également de piloter la taille de l’arête
des mailles en fonction de leurs positions dans le domaine étudié. Dans cette étude, seule
l’interaction de l’onde de choc avec une structure et l’évolution des caractéristiques
mécaniques de l’onde de souffle en aval du bâtiment sont étudiées, figure 55. Le reste du
domaine n’appartient pas au périmètre de l’étude. Les mailles situées dans cette dernière zone
possèdent une taille d’arête plus importante que les mailles situées dans la zone d’intérêt.
Cette gestion intelligente et totalement automatique au cours du calcul de la taille des mailles
permet de diminuer leur nombre et donc le temps de restitution.
y Lmax
6
Une étude comparative d’un cas de calcul simulé avec, et sans AMR, a été réalisée pour un même niveau de discrétisation
spatiale en géométrie 2D. Les résultats de cette étude permettent de conclure que l’utilisation de l’AMR, dans le cadre de ces
travaux de thèse, n’introduit pas une différence sur l’évolution des caractéristiques mécaniques de l’onde de souffle.
78
grad (P)
L’objectif de cette étude comparative est d’évaluer la marge d’erreur entre les résultats
obtenus par ces deux approches géométriques, géométrie 2D-axisymétrique et géométrie 3D,
à partir d’une configuration expérimentée lors de l’étude préliminaire du projet ANR
BARPPRO (Eveillard (2010)) et donc de pouvoir quantifier les écarts.
Charge
(R1 = 0,06 m)
(0,0) H = 0,20 m
α1 = 45° α2 = 45°
x
L = 0,80 m
x
y
y
z y
Domaine d’étude
y en géométrie Domaine
z d’étud
x
7
Une étude de l’anisotropie sur le code de calcul HERA a été réalisée. Les résultats de cette étude indiquent un écart de
maximum 2-3% pour une approche en géométrie 2D-axisymétrique avec un maillage « fin » et « grossier [a] ». L’écart est
légèrement plus important en géométrie 3D avec une variation au maximum de 7% pour le maillage « grossier [a] ».
81
y
d = 0,14 m
Charge 0,137 m
(0, 0) 0,103 m
α1 = 45° x
Pour le capteur situé à 0,137 m (soit 0,57 m/MJ1/3), le maillage « grossier [a] » ne
permet pas de faire la distinction entre l’onde incidente et l’onde réfléchie (capteur très proche
de la structure). Les deux ondes sont confondues, une seule valeur du maximum de
surpression est alors mise en évidence. Cette valeur est logiquement plus forte que chacune de
celles attendues sur le maximum de surpression de l’onde incidente et de l’onde réfléchie
82
puisque la taille de mailles employée induit la fusion de ces deux ondes. La différence entre la
géométrie 2D-axisymétrique et la géométrie 3D est faible en amont de la barrière avec un
intervalle de - 3% à + 0,6% pour le maximum de surpression (domaine n°5). Les écarts entre
les différents domaines pour le maximum de surpression et l’impulsion positive sont présentés
dans le tableau 21.
Capteurs placés en amont du merlon Capteur à 0,103 m Capteur à 0,137 m
Caractéristiques ∆P+ (en bar) I+ (en bar.ms) ∆P+ (en bar) I+ (en bar.ms)
Données expérimentales 5,18 0,19 3,27 0,21
n°1. 2D-axisymétrique
5,01 0,18 3,41 0,19
maillage « fin »
- 3% - 5% + 3% - 10%
(% par rapport aux données expérimentales)
n°2. 2D-axisymétrique
3,31 0,18 4,80 0,24
maillage « grossier [a] »
- 34% identique + 40% + 26%
(% par rapport 2D « fin »)
Les nouvelles valeurs de références sont les résultats de la géométrie 2D-axisymétrique « grossier [a] »
n°3. Révolution 3D du 2D
3,32 0,19 4,64 0,24
maillage « grossier [a] »
+ 0,3% + 5% - 3% Identique
(% par rapport 2D)
n°4. 3D merlon infini
3,33 0,17 4,65 0,21
maillage « grossier [a] »
+ 0,6% - 6% - 3% - 12,5%
(% par rapport 2D)
n°5. 3D merlon expérimenté
3,33 0,17 4,65 0,21
maillage « grossier [a] »
+ 0,6% - 6% - 3% - 12,5%
(% par rapport 2D)
Tableau 21 : Ecarts sur le maximum de surpression et d’impulsion positive entre les domaines pour les premiers
points de mesures placés en amont du merlon
8
Par la suite, une mousse et plusieurs sacs de sable sont placés dans la structure afin d’éviter la présence de réflexions
multiples internes à la structure creuse qui perturbent la mesure pour les capteurs placés sur les surfaces du merlon. La fixation
des capteurs sur la structure est également renforcée afin d’éviter tout mouvement du capteur sous l’effet de l’impact de l’onde
de choc.
83
y Point de mesure
d = 0,14 m
0,35 m
Charge
(0, 0) x
Figure 60 : Profils de surpression au niveau du capteur placé sur la face avant du merlon
9
Maille mixte : élément de discrétisation qui comprend plus d’un matériau ; nécessite un traitement numérique spécifique par le
code de calcul.
Dans le cadre de ces travaux de thèse, ce traitement est réalisé par une pondération sur la pression (P) et le volume (V) en
fonction de la fraction de présence des matériaux dans la maille mixte.
La structure est considérée comme un matériau indéformable. Ce matériau posséde une masse volumique importante
(1050 kg/m3), soit une impédance qui peut être considérée comme quasi-infinie, et la vitesse (u) est imposée nulle à chaque pas
de temps.
84
y d = 0,14 m
B
Charge
(0, 0) A
x
b) évolution des écarts le long de la face avant du merlon entre les deux approches
géométriques et les deux tailles de mailles utilisées
y
0,17 m
d = 0,14 m
Point de mesure
Charge
(0, 0) x
Pour ce capteur placé au sommet du merlon, il est possible de constater une dispersion de 6%
entre les calculs réalisés pour un maillage « grossier [a] » que la forme du merlon soit torique
ou bien celle réellement expérimentée. La mise en œuvre d’un maillage « fin » en géométrie
2D-axisymétrique permet une convergence numérique satisfaisante avec un écart inférieur de
5% par rapport à la mesure. Pour une même approche, géométrie 2D-axisymétrique, l’écart
entre les deux tailles de mailles utilisées (maillage « fin » et maillage « grossier [a] ») est de
17%. Il est donc préférable de calculer en géométrie 2D-axisymétrique « fin » qu’en
géométrie 3D « grossier [a] » en choisissant un meilleur niveau de convergence numérique
plutôt que la géométrie exacte du merlon. δ’évolution du maximum de surpression tout le
long du sommet du merlon permet de distinguer les différences entre les domaines d’étude,
figure 63 (extraction à 5 mm de la paroi).
86
y d = 0,14 m
A
Charge B
(0, 0)
Dans cette région, le calcul de la géométrie étudiée dans une configuration en géométrie 3D
n’est pas nécessaire. δa configuration en géométrie βD-axisymétrique sur maillage « fin »
apporte plus que la prise en compte de la forme réelle du merlon. Elle permet ainsi d’assurer
un niveau de convergence numérique satisfaisant et conduit à la restitution de la mesure.
y
0,20 m
d = 0,14 m
0,05 m
Point de
Charge mesure
(0, 0) x
Toutes les configurations (géométrie 2D-axisymétrique et géométrie 3D), et les deux tailles
de mailles (« fin » et « grossier [a] ») fournissent des résultats proches de la mesure dans la
zone de détente. δ’évolution du maximum de surpression le long de la face aval du merlon est
présentée sur la figure 65 (extraction à 5 mm de la paroi).
y
d = 0,14 m
A
Charge
(0, 0)
B
x
Les résultats de simulation numérique obtenus pour les différents cas de calcul avec le
maillage « grossier [a] » sont quasi-identiques. δ’approche en géométrie βD-axisymétrique
(merlon torique) et l’approche en géométrie γD (merlon réellement expérimenté) conduisent
donc aux mêmes résultats, aucun effet de compression propre à l’approche en géométrie
2D-axisymétrique n’est visible dans cette zone.
δa mise en œuvre du maillage « grossier [a] » induit une surestimation du maximum de
surpression (maximum 46% par rapport au maillage « fin »). Dans cette zone d’intérêt,
l’utilisation d’un maillage plus convergé numériquement est donc préférable. δ’approche en
géométrie 2D-axisymétrique « fin » permet ainsi d’assurer un niveau de convergence
numérique satisfaisant et donc la restitution de l’expérience.
Analyse de l’écart dans la zone d’intérêt du projet ANR BARPPRO : zone protégée
La figure 66 présente les profils de surpression pour deux points de mesure placés en
aval du merlon et la figure θ7 permet de visualiser l’évolution du maximum de surpression en
aval de la barrière de protection. Le maximum de surpression augmente en aval du merlon
lors de la réflexion de l’onde sur le sol. δ’onde subit dans un premier temps une Réflexion
Régulière (RR) avant de passer en Réflexion de Mach (RM).
88
0,20 m
y d = 0,14 m
0,80 m (4 H)
Charge
(0, 0)
0,40 m (2 H)
L’écart entre tous les cas de calculs avec un maillage « grossier [a] » est faible en aval du
merlon (inférieur à 1% pour les deux capteurs, figure 66). δ’écart entre les données
expérimentales et les résultats de simulation numérique obtenus pour l’approche en géométrie
2D-axisymétrique « fin » est au maximum de 7%, il est donc compris dans l’incertitude des
mesures. Pour une même approche géométrique donnée, 2D-axisymétrique, l’écart est
d’environ 30% entre les maillages « grossier [a] » et « fin ».
Cette région en aval du merlon correspond à la zone d’intérêt du projet ANR BARPPRO.
Dans cette zone, dite « zone protégée », l’approche βD-axisymétrique avec un maillage « fin »
pour la géométrie de merlon étudiée est bien plus satisfaisante qu’une approche en géométrie
3D sur un maillage « grossier [a] » qui respecterait totalement la forme du merlon
expérimenté. Ce cas de calcul, en géométrie 2D-axisymétrique, permet de restituer
l’expérience avec une précision de l’ordre de 10% (intervalle d’incertitude expérimentale).
La conclusion obtenue dans le cas d’une géométrie particulière lde merlon (merlon de
e
l’étude préliminaire du projet ANR BARPPRO, Eveillard (2010)) est étendue (écart de
e
91
maximum 10% entre une approche en géométrie 2D-axisymétrique et une géométrie 3D) à
toutes les configurations expérimentées lors de la constitution des abaques (étude
paramétrique).
δ’étude de l’évolution de l’impulsion positive (I+) ne fait pas partie du périmètre du
projet ANR BARPPRO. Néanmoins, cette caractéristique mécanique de l’onde de souffle est
étudiée lors de la présentation des différents historiques de pression. Le calcul de cette
caractéristique (I+) par l’outil de simulation est moins dépendant de la taille de mailles que le
calcul du maximum de surpression (ΔP+) (écart moyen de 20% par rapport aux données
expérimentales).
10
Convergence numérique : par abus de langage, on appelle une convergence numérique, un résultat qui approche la valeur
numériquement convergée (écart de quelques %, inférieur aux incertitudes (expérimentales ou bibliographiques)).
92
93
Figure 69 : Profil de surpression de l'étude préliminaire du projet ANR BARPPRO (Eveillard (2010))
0,10 m
Charge
(0,0) Capteur de référence (CR1) Structure
étudiée
Figure 71 : Profils de surpression pour le capteur de référence (CR1) à 0,1 m du centre de la charge
(soit 0,61 m/MJ1/3)
18 cm
36 cm
14 cm
40 cm 36 cm
Figure 72 : Présentation des obstacles étudiés
Le tableau 24 présente les tailles de mailles utilisées lors de l’étude des trois configurations de
Trélat (2006). Les simulations numériques en géométrie 2D-axisymétrique « fin » respectent
le standard de calcul défini précédemment (0,13 mm, soit 0,8 mm/MJ1/3 pour les arêtes des
mailles). δa mise en œuvre de ce maillage permet de se placer dans les conditions de
convergence numérique satisfaisante. δ’étude des ondes de contournement autour du cylindre
et du parallélépipède nécessite l’utilisation d’une approche en géométrie 3D. Pour cette
approche géométrique tridimensionnelle, la dimension des arêtes des mailles est fixée à
1,77 mm au passage du front de choc (soit 10,8 mm/MJ1/3). Cette taille de maille (maillage
« grossier [b] ») est plus petite que la taille limite définie précédemment et notée maillage
« grossier [a] ».
Energie de la Type de maillage
Configurations
charge de gaz « fin », arêtes de mailles de 0,13 « grossier [b] », arêtes de mailles de
étudiées
(en MJ) mm (standard, soit 0,80 mm/MJ1/3) 1,77 mm (soit 10,8 mm/MJ1/3)
Géométrie 2D-axisymétrique
1 – capteur CR1 Géométrie 3D « grossier [b] »
« fin »
Géométrie 2D-axisymétrique Géométrie 2D-axisymétrique et
2 – pavé 1
« fin » géométrie 3D « grossier [b] »
4,4.10-3
Géométrie 2D-axisymétrique
3 – pavé 2 Non calculé
« fin »
Géométrie 2D-axisymétrique
4 – cylindre Géométrie 3D « grossier [b] »
« fin »
Tableau 24 : Etude des configurations de Trélat (2006)
97
20 cm
3 cm
Charge 3 cm T3 T0 T7 14 cm
(0, 0, 0)
Lg0
y Ld0
7 cm
Lg1
9 cm
3 cm
x
7 cm
5 cm
CA6
CR1 5 cm
z 9 cm
15 cm
40 cm
a) configuration 2
15 cm 14 cm
T1 T0 T5
Charge
(0, 0, 0) 3 cm
9 cm
R0
y F0
20 cm
7 cm
R1
3 cm
x CR1 6 cm CA4
7 cm
6 cm
z
20 cm
16 cm
18 cm
b) configuration 3
36 cm
T1 T3
15 cm T0 15 cm
F1 R1
3 cm
36 cm
33 cm
33 cm
F2 R2
Charge
18 cm
18 cm
(0, 0, 0)
y
3 cm
F3 R3
3 cm
x 2 cm 2 cm
CR1 CA5
z
12 cm
c) configuration 4
y d2 = 15 cm e2 = 40 cm
H = 14 cm
Capteur de référence (CR1)
x
(0, 0)
a) configuration 2
y d3 = 16 cm e3 = 18 cm H = 14 cm
x
(0, 0)
b) configuration 3
Figure 75 : Profils de surpression au niveau du capteur de référence (CR1) pour les configurations 2 et 3
Maximum de
surpression (bar)
2,130
1,540
0,947
0,354
-0,239
Figure 76 : Réflexion régulière sur la face avant de l'obstacle (configuration 2, t = 0,28 ms)
Mach (passage pour un nombre de Mach proche de 1,4, d’après la formule simplifiée de
Kinney (1985)). Il est donc possible qu’il existe un régime transitoire entre une réflexion
régulière et une réflexion de Mach au sommet. Notons que la configuration 2 est à priori dans
la même situation, l’évolution de la surpression n’indique cependant pas ce changement. Ceci
peut en partie s’expliquer par la complexité, la non-linéarité de cette phénoménologie et la
différence de 7% entre les distances charge-face avant du pavé entre les deux configurations.
Points
d’observation Nombre de Régime de
Numéro de Angle
Mach incident réflexion
configuration (en m) d’incidence ( )
(Trélat (2006)) (Kinney (1985))
X Y
0,150 0,000 0o 1,44 RR
2 0,150 0,126 40o 1,29 RR
0,150 0,140 43o 1,27 ≈ RR
0,160 0,000 0o 1,40 RR
3 0,160 0,134 40o 1,28 RR
0,160 0,140 41o 1,27 ≈ RR
RR : Réflexion Régulière ;
≈ RR : Réflexion Régulière, mais proche d’une Transition vers une Réflexion de Mach (TRM)
Tableau 25 μ Régime de réflexion le long de la face avant de l’obstacle d’après Kinney (1λ85)
Dans ses travaux, Trélat (2006) a établi plusieurs formules empiriques permettant
d’estimer les caractéristiques d’une onde de choc réfléchie en fonction de la distance réduite
entre le centre de la charge et le point d’impact ( ).
La surpression réfléchie peut être évaluée à partir de deux paramètres μ l’angle d’incidence ( )
et le maximum de surpression incidente (ΔP+) (voir § 1.γ.1). Selon la valeur et l’évolution de
ces deux paramètres, il est possible d’estimer le maximum de surpression réfléchie (ΔPr+)
(par exemple, abaque du TM5-1300 (2008)). Les régimes de réflexion et les transitions
possibles depuis une réflexion régulière vers une réflexion de type Mach peuvent être estimés
à partir de la formule simplifiée de Kinney (1985) et de l’équation (7γ) qui détermine le
nombre de Mach incident (avec = 1,ζ).
101
a) configuration 2
b) configuration 3
Si on considère uniquement une surface verticale (α = λ0°) pour une onde incidente sphérique
divergente, alors il est nécessaire de prendre en compte trois paramètres ( , d et ). Pour cette
configuration d’obstacle, l’angle d’incidence ( ) dépend de la distance réduite ( ) et de la
distance entre le centre de la charge et la structure (d), figure 78.
d1
d2
d
Angle incident :
λ2 β2 >> 0o
= cos
Structure 1
d1 > d 2 → 2 > 1
α1 = 90o soit ΔPr1+ ≠ ΔPr2+
β1 = 0o
λ1
Figure 78 : Réflexion sur une paroi verticale avec et sans influence de l’angle d’incidence ( )
La formule proposée par Trélat (2006) qui estime la surpression réfléchie (ΔPr+) se
base sur des données expérimentales. Ces résultats sont obtenus à partir de différentes
configurations pour un angle d’incidence ( ) entre 8,θηo et 47,26o, avec une surface verticale
placée entre 0,10 et 0,β0 m du centre de la charge (d) et des distances réduites ( ) comprises
entre 0,43 et 2,37 m/M1/3, relation (138).
Pr
ln 1,270 1,910 ln( ) 0,290 ln( ) 2
0
(138)
P
Avec P0 = 1 bar, ΔPr+ (en bar) et 0,ζγ ≤ ≤ β,γ7 m/εJ1/3
δ’auteur précise également que cette formule prend en compte les possibles
changements entre les deux principaux régimes de réflexion, régulière ou de Mach.
δes intervalles d’études de Trélat (β00θ) peuvent effectivement conduire à l’observation d’un
changement de régime. Cependant, ce changement n’est pas linéaire et dépend également de
l’angle d’incidence ( ). δ’utilisation d’une formule basée uniquement sur la distance réduite
( ) pour estimer les caractéristiques de l’onde réfléchie (ΔPr+, dtr+, Ir+) pour ces deux régimes
de réflexion, et qui ne tient pas compte de l’angle d’incidence ou du régime de réflexion,
conduit à un modèle qui se révèlera insuffisant.
Les résultats de simulation numérique sont également comparés à la formule (138) qui
estime la surpression réfléchie en fonction de la distance réduite ( ). Cette estimation semble
surestimer la surpression réfléchie (ΔPr+) d’environ β0% par rapport aux résultats de
simulation numérique et du TM5-1300 (2008). Ceci laisse à penser que la formule de Trélat
(2006) pour estimer le maximum de surpression réfléchie (ΔPr+) dans le cas d’une structure
très proche de la charge explosive montre ici ses limites. Elle devrait en effet prendre en
compte l’angle d’incidence ( ) et en l’état n’est donc valable que pour les cas où la surface est
suffisamment éloignée de la charge pour pouvoir faire l’hypothèse d’un rayon de courbure
infini du choc incident.
103
Figure 79 : Profils de surpression pour le capteur F0 situé à mi-hauteur de la face avant de l’obstacle
(configuration 3)
Pour le capteur (F0), les données expérimentales de Trélat (2006) surestiment largement
(> 40%) la surpression réfléchie sur la face avant du parallélépipède par rapport à l’abaque du
TM5-1300 (2008) et aux résultats de la simulation numérique, voir figure 79. Cette
surestimation excède l’intervalle de confiance expérimental établi par l’analyse des mesures
de pression du capteur de référence en champ libre (voir §3.1.1).
Les résultats de simulation numérique sont en accord avec les estimations de la surpression
réfléchie du TM5-1300 (2008) (écart de 3% pour le capteur F0), validant ainsi les conditions
de calculs.
3.1.3.2 Etude de la détente au sommet de l’obstacle
Sur le sommet de l’obstacle, pour les deux positions du pavé, ce n’est pas une onde de
contournement mais l’onde incidente qui, après avoir été réfléchie sur la face avant de
l’obstacle, va subir une diffraction puis une détente le long du sommet, figure 80.
Maximum de Maximum de
surpression (bar) surpression (bar)
1,120 1,320
0,738 0,890
0,357 0,464
-0,025 0,039
-0,407 -0,386
a) configuration 2 b) configuration 3
Figure 80 : Détente de l'onde incidente au sommet du parallélépipède pour chaque configuration (t = 0,35 ms)
voir figure 81. La trajectoire de la détente n’est pas à pente douce comme dans le cas d’un
merlon ayant une épaisseur nulle au sommet (merlons de Allain (1994) et de Borgers (2010),
e = 0 m). En revanche pour cette configuration, il est possible de faire l’hypothèse d’une
détente en deux étapes : détente au sommet de la barrière de protection (e ≠ 0 m), puis une
détente le long de la pente arrière de la structure (cas d’un merlon infini).
a) configuration 2
b) configuration 3
Figure 81 : Evolution du maximum de surpression au sommet du parallélépipède
Il est important de rappeler ici que si la réflexion sur la face avant avait aboutie à la formation
complète d’un pied de Mach, c’est le pied et non l’onde incidente qui aurait subi la détente,
changeant alors son comportement durant le phénomène de détente. Ce cas de figure aurait pu
être envisagé si la hauteur du parallélépipède avait été supérieure (H > 0,14 m), avec un
passage possible vers une réflexion de Mach.
105
a) configuration 2 b) configuration 3
Figure 82 : Profils de surpression pour le capteur T 0 situé sur le sommet de la structure
Cette interprétation est valable uniquement si les ondes de contournement peuvent être
négligées en intensité par rapport à l’onde issue du sommet (∆P+ onde issue du merlon >> ∆P+
onde de contournement). C’est le cas si la distance à parcourir par les ondes de contournement
(C) pour se recombiner en aval est significativement supérieure à la distance parcourue par
l’onde qui passe par le sommet (S) (C >> S). Cette interprétation s’applique à la configuration
3, où C = 0,58 m, soit 3,54 m/MJ1/3, et S = 0,46 m, soit 2,81 m/MJ1/3 (soit C ≈ 1,βθ S),
figure 84.
20 cm
18 cm
18 cm
14 cm
14 cm
Structure
20 cm
C = 58 cm S = 46 cm
C >> S
Figure 84 : Représentation schématique des distances à parcourir autour de l’obstacle pour la configuration γ
9 cm
40 cm Capteur (CA6)
14 cm
14 cm
40 cm
Structure 5 cm
S = 68 cm
9 cm
C = 58 cm S >> C
Figure 85 : Représentation schématique des distances à parcourir autour de l’obstacle pour la configuration β
Une recombinaison des ondes de choc est donc possible après le parallélépipède, permettant
une augmentation de la surpression supérieure à la simple réflexion sur le sol de l’onde issue
du sommet, figure 86. Afin de vérifier cette hypothèse, une simulation en géométrie 3D est
nécessaire pour compléter l’étude de la configuration β. Les conditions de gestion des tailles
des mailles sont alors celles présentées dans le tableau 24. δa simulation s’effectue donc sur
un maillage « grossier [b] » tel que la convergence numérique sera moindre que pour les
107
x
Onde de contournement
Onde incidente
Onde issue de la z
Structure
recombinaison
Onde de contournement
Configuration 2
δa longueur de l’obstacle face à l’explosion est de 0,18 m, les ondes de contournement
vont donc parcourir la même distance que l’onde qui passe par le sommet pour atteindre la
face en aval de l’obstacle. δa distance à parcourir pour atteindre le capteur (CA6) est un peu
plus importante pour l’onde issue du sommet que pour les ondes de contournement (0,69 m
contre 0,63 m). Les résultats de simulation en géométrie 2D-axisymétrique ne permettent
donc pas d’analyser l’évolution du maximum de surpression après l’obstacle car la simulation
ne prend en compte que l’onde issue du sommet. δ’hypothèse d’un obstacle de longueur, face
à l’explosion, infinie n’est plus vérifiée.
Les figures 87-a et 87-b présentent les résultats obtenus pour le capteur de référence (CR1) et
l’évolution du maximum de surpression sur la face aval du pavé. Le capteur de référence
permet d’évaluer une sous-estimation d’environ β0% par rapport à la formule empirique (10)
de Trélat (2006). Cet écart s’explique par la taille des mailles qui est plus importante que le
calcul en géométrie 2D-axisymétrique « fin ». Ce défaut de simulation est purement lié à un
défaut de convergence numérique. Néanmoins, une valeur du maximum de surpression
acceptable vis-à-vis du polynôme champ libre de Trélat (2006) est bien retrouvée. Pour
représenter l’évolution de la surpression maximale sur la face arrière de l’obstacle, l’abscisse
rectiligne (distance au centre de la charge) utilisée habituellement est remplacée par l’abscisse
curviligne.
Le graphique sur la figure 87-b permet de mettre en évidence l’absolue nécessité de
simuler en géométrie 3D les cas de merlons « courts » (faible longueur de la barrière face à la
charge) pour lesquels les effets de contournement ne peuvent pas être négligés. Les ondes de
contournement vont subir les mêmes phénoménologies que l’onde issue du sommet sur les
faces avant et latérales avant de se recombiner toutes les trois sur la face arrière. Après cette
focalisation, l’onde de choc subit une recompression avec une augmentation d’un facteur β
par rapport à une onde seule (simulation 2D). Pour cette configuration, la recombinaison des
ondes après l’obstacle est la phénoménologie prédominante.
En dépit de la mise en œuvre d’un maillage apparemment dégradé, c’est bien le calcul en
géométrie γD qui peut permettre d’assurer la restitution des points de mesure de
Trélat (2006), sur la face aval du pavé.
108
a) capteur de référence (CR1) à 0,10 m du centre de la charge, soit 0,61 m/MJ 1/3
Configuration 3
Dans cette configuration, la plus grande longueur du parallélépipède est face à la
charge, les effets de contournement sont négligeables. La simulation de cette configuration en
géométrie 2D-axisymétrique est suffisante. Elle a donc pu être effectuée avec le maillage
« fin » du standard de calcul.
Les résultats obtenus pour le capteur R0 et l’évolution du maximum de surpression sur
la face aval du pavé sont présentés respectivement sur les figures 88-a et 88-b. Les données
expérimentales de Trélat (2006) pour le capteur (R0) situé au centre de la face arrière (R0)
montrent un écart d’environ + 10% par rapport aux résultats de calculs, voir figure 88-a. La
présence du phénomène de détente conduit à une déformation de la décroissance de
surpression (montée brutale du signal de crête atténuée). D’après ce point de mesure, il est
possible d’observer le passage de l’onde issue du sommet qui subit une détente le long de la
paroi, puis le passage de l’onde réfléchie sur le sol. Ces observations sont également visibles
sur l’évolution du maximum de surpression le long de la face arrière de l’obstacle.
Pour le point de mesures (R1), l’écart entre les données expérimentales et les résultats
de simulation numérique est plus important que pour le capteur précédent (R0), environ
+ 40% par rapport à la simulation numérique. Ceci peut s’expliquer par la présence de l’onde
réfléchie sur le sol en aval de l’obstacle après la détente de l’onde de choc. Cette
recompression serait alors mesurée avec une précision moindre.
En effet, cet écart avait déjà été observé pour le capteur situé sur la face avant (F0, figure 79)
de cette même configuration où la référence bibliographique du TM5-1300 (2008) donnait
raison à la simulation (moins de η% d’écart) aux dépens de la mesure de Trélat (2006). Peut-
être s’agit-il d’un problème de la chaîne de mesures de Trélat (β00θ) lorsqu’un capteur
mesure l’effet de recompression d’une onde de choc.
Le même type de constat que pour le capteur (F0) peut être fait sur le capteur (CA4) au sol en
aval de l’obstacle, figure 8λ. δ’écart observé est ainsi du même ordre de grandeur, + 21%, que
pour le capteur situé sur la face avant de la structure. Un décalage chronométrique de 39%
peut être également constaté.
110
Le second front de choc est issu de la recombinaison des ondes de contournement en aval du
parallélépipède. La présence des ondes de contournement dans cette zone d’intérêt permet de
constater les limites de l’hypothèse d’un pavé de longueur infinie μ nécessité d’une simulation
en géométrie 3D en aval du pavé (capteur CA4).
explosive est proche de la structure. La courbure de l’onde de choc ne peut être alors négligée
comme le suppose Trélat (2006) pour établir la formule (138).
δes capacités restitutives de l’outil de simulation ont donc été confirmées pour une
configuration parallélépipédique lors de la détonation d’une charge gazeuse. Une
configuration cylindrique fait l’objet d’une étude dans la prochaine partie de ce chapitre.
Figure 90 : Profils de surpression pour le capteur de référence (CR1), configuration 4 de Trélat (2006)
sur l’avant du cylindre) de β0%. Cette valeur, majorante en champ proche, est conservée pour
les valeurs du maximum de surpression sur le cylindre et dans son voisinage pour les calculs
en géométrie 3D.
Cette étude pour la configuration 4 se base sur deux cas de calcul : une géométrie
2D-axisymétrique « fin » et une géométrie 3D avec un maillage dégradé, « grossier [b] »,
figure 91. Une étude en géométrie 2D-axisymétrique permet d’étudier le régime de réflexion
et le niveau de surpression le long de l’avant du cylindre ( = 0°). δ’étude de cette
configuration est ensuite complétée par une simulation en géométrie 3D sur un maillage
« grossier [b] » pour tenir compte du phénomène de réflexion autour de la surface cylindrique
(0° ≤ ≤ 90°) ainsi que des effets des ondes de contournement et ceci même si les niveaux de
pression prévus par simulation peuvent être sous-estimés jusqu’à β0% du fait du niveau de
convergence numérique moindre qu’un maillage « fin ».
y d = 12 cm e = 36 cm
H = 36 cm
Charge
(0,0)
x
a) étude en géométrie 2D-axisymétrique (cylindre de forme torique)
d = 12 cm e = 36 cm
= 0°
H = 36 cm
cylindre. Seule la simulation avec un maillage « fin » permet d’y parvenir (calcul en
géométrie 2D-axisymétrique).
Les données expérimentales de Trélat (2006) présentent des écarts qui peuvent être
extrêmement importants avec la référence bibliographique (TM5-1300 (2008)). Cette
différence est supérieure à l’intervalle de confiance, seul le capteur F2 présente une mesure
qui est proche de la valeur estimée par le TM5-1300 (2008). Les deux autres capteurs (F1 et
F3) présentent des écarts très importants. Pour valider de manière discriminante la chaîne de
simulation, la comparaison sera faite essentiellement sur les valeurs fournies par l’abaque de
réflexion du TM5-1300 (2008).
(0,0) (0,0)
x x
Par la suite, l’onde réfléchie qui se forme le long de la génératrice la plus proche de la
charge ( = 0°), qu’elle soit de type réfléchie ou de type εach, devient l’onde incidente de la
surface convexe et influence la nature des ondes de contournement après l’obstacle. δa nature
de l’onde réfléchie à la fin de la surface convexe peut être déterminée à partir des
observations de Ben-Dor (1991), tableau 26.
Charge gazeuse (propane-oxygène à stœchiométrie), E = ζ,ζ.10-3 MJ
Points Surpression Régime de Régime de
Distance Angle Nombre de
d’observation incidente (ΔP+) réflexion à réflexion à
réduite d’incidence Mach incident
(en m) (en bar) = 0° = λ0°
(en m/MJ )1/3
( ) (formule (73))
X Y (formule (10)) (Kinney (1985)) (Ben-Dor (1991))
0,03 o RR → Rε
0,755 14 1,82 1,60 RR
(F3) (figure 95)
0,10 0,956 40o 1,26 1,44 RR RR → Rε
0,12 1,036 45o 1,17 1,41 RM Rε → Rε
0,12
0,18 Rε → Rε
1,310 56o 0,68 1,25 RM
(F2) (figure 96)
0,33
2,141 70 o
0,31 1,12 RM Rε → Rε
(F1)
RR : Réflexion Régulière ; RM : Réflexion de Mach
Tableau 26 : Régimes de réflexion sur l’obstacle cylindrique
(0,0) (0,0)
x x
z z
Réflexion régulière (RR), = 0° (t = 0,15 ms) Réflexion de Mach (RM), > 0° (t = 0,20 ms)
(0,0) (0,0)
x x
z
z
Réflexion régulière (RR), = 0° (t = 0,35 ms) Réflexion de Mach (RM), > 0° (t = 0,40 ms)
(0,0)
x
Trélat (2006) a disposé trois capteurs au sommet de l’obstacle dans l’axe de la charge, noté
(T1), (T0) et (T3). Seul le capteur (T1) situé au sommet de l’obstacle au niveau de la première
détente du pied de Mach est ici étudié. Le maximum de surpression obtenu par simulation est
bon en accord avec les données expérimentales de Trélat (2006) (écart de 16%, dans
l’intervalle de confiance des mesures de Trélat (β00θ)), mais il s’agit peut-être d’un effet
compensatoire. En effet, la convergence numérique n’est pas forcément atteinte ici sur le
maillage « grossier [b] » mis en œuvre (sous-estimation possible jusqu’à 20%, voir l’analyse
sur le capteur de référence (CR1)), figure 98.
Détente et contournement après le cylindre
δe long de la face arrière de l’obstacle pour 90° ≤ ≤ 270°, plusieurs
phénoménologies sont visibles en fonction de la nature de l’onde : l’onde issue du sommet ou
les ondes de contournement, figure 99. Trélat (β00θ) ne fournit pas d’historique de
surpression pour les capteurs placés sur les côtés du cylindre. De ce fait, seule une analyse
qualitative des résultats de simulation est proposée.
δ’onde issue du sommet de l’obstacle est l’onde qui passe par-dessus le cylindre. Le long de
l’arrière du cylindre, elle subit une détente avant de se recombiner avec les ondes de
contournement. Les ondes de contournement correspondent aux ondes qui se propagent
autour du cylindre (surface convexe). Après le passage de l’obstacle, les ondes de
contournement vont se focaliser à l’arrière du cylindre, le long de la génératrice = 180°,
compte tenu de la symétrie de la géométrie.
Capteur (CA5)
(0,0)
Vue dans le plan y= 0 m
z
Passage de l’onde
du sommet et
Recombinaison
des ondes de
contournement
R1
Recombinaison
des ondes de
contournement 0,33 m
, avant 0,26 m
l’arrivée de
l’onde issue du
sommet
Les points de mesure situés en dessous d’une hauteur de 0,26 m le long de la face aval du
cylindre doivent permettre de constater la recombinaison des ondes de contournement et d’en
mesurer les effets. C’est le cas des capteurs (R2) et (R3), figures 101 et 102.
Passage de l’onde
du sommet et
Recombinaison
des ondes de
contournement
Recombinaison
des ondes de
contournement R2
, avant 0,26 m 0,18 m
l’arrivée de
l’onde issue du
sommet
Passage de l’onde
du sommet et
Recombinaison
des ondes de
contournement
Recombinaison
des ondes de
contournement
, avant 0,26 m
l’arrivée de
l’onde issue du R3 0,03 m
sommet
Pour les capteurs (R2) et (R3), il est possible d’observer le passage d’une seule onde de choc
qui est issue de la recombinaison des ondes de contournement et de l’onde provenant du
sommet du cylindre (géométrie tridimensionnelle, points de mesure placés en dessous de
0,26 m). Pour le capteur (R2), les résultats de simulation numérique se situent dans l’intervalle
de confiance des mesures de Trélat (2006) de +/- 20% (écart de 9% pour le maximum de
surpression et de 10% pour le temps d’arrivée). Pour le capteur (R3), l’écart entre les données
expérimentales de Trélat (2006) et les résultats de simulation numérique est très important
(écart de 80% sur le maximum de surpression).
δors de l’étude préliminaire du projet ANR BARPPRO (§ 2.4), une configuration de merlon
creux avait été expérimentée. Après analyse des signaux mesurés, il s’est avéré que certains
capteurs placés sur la structure ont été perturbés par des réflexions multiples internes à la
structure creuse, le merlon jouant ainsi le rôle de caisse de résonance. Ces mesures avaient
donc été retirées et à la place des mesures des valeurs issues de la littérature avaient donc été
utilisées (TM5-1300 (2008), Brode (1968)). δ’allure générale des signaux perturbés par ces
réflexions internes à la structure était similaire au signal acquis par Trélat (2006) pour le
capteur (R3). Les structures expérimentées par Trélat (2006) étant également creuses, il est
possible d’émettre l’hypothèse que ce point de mesure est également entaché du même
problème de mesure : formation de réflexions multiples internes à la structure qui perturbent
l’acquisition de la mesure pour le capteur (R3).
Le capteur (CA5) situé au sol et en aval de la structure subit également le passage de
la recombinaison des ondes de choc issues du contournement de l’obstacle sur la face arrière
du cylindre ( = 180°), figure 103.
120
Passage de l’onde
du sommet et
Recombinaison
des ondes de
contournement
Recombinaison
des ondes de
contournement
Capteur (CA5)
0,02 m
Sens de propagation
Figure 104 : Profils de surpression pour des capteurs placés en champ libre en géométrie 2D-axisymétrique
« fin »
122
E
td
1/ 3
2,081 1 e 0.3361 (141)
E
Avec P0 = 1 bar, ΔP+ (en bar), t0 (en ms), dt+ (en ms), E (en εJ) et 0,8ζ ≤ ≤ 17,η m/MJ1/3
Ces nouvelles formules empiriques peuvent être utilisées jusqu’à 17,5 m/MJ1/3 (limite de la
zone d’étude) avec une incertitude acceptable (+/- 10%), figure 105.
La décroissance de pression peut être estimée à partir du model de Friedlander modifié,
formule (142) (Baker (1973)). Cette relation se différencie du premier modèle (formule (1))
par l’ajout d’un coefficient (C) dont la valeur dépend de la distance réduite ( ).
t t ( ) t t 0 ( )
P( , t ) Pmax ( ) 1 0 exp C ( )
dt ( ) dt ( )
(142)
Figure 105 : Evolution du maximum de surpression (∆P+, en bar), du temps d’arrivée (t0, en ms/MJ1/3), de la
durée de phase positive (td+, en ms/MJ1/3) en fonction de la distance parcourue ( , en m/εJ1/3) pour une charge
gazeuse (propane-oxygène à la stœchiométrie)
124
δa mise en œuvre du maillage « grossier [c] » induit une sous-estimation des résultats
de - 20% sur le capteur A1 placé à 0,83 m/MJ1/3 (au lieu d’environ 10% pour le maillage
« fin », ordre de grandeur des incertitudes expérimentales), figure 106. Cet écart provient d’un
niveau de convergence numérique insuffisant pour le maillage « grossier [c] ».
Figure 106 : Profils de surpression pour des capteurs placés en champ libre en géométrie 3D avec la mise en
œuvre du maillage « grossier [c] »
Les valeurs du maximum de surpression sont extraites sur le maillage. Afin d’éviter le bruit
numérique par la présence de mailles mixtes11 entre l’air et la structure, l’extraction est
réalisée à une distance de 1 mm de la paroi. Notons que les valeurs acquises à moins de
0,4 m/MJ1/3 du centre de la charge explosive peuvent être sujettes à un bruit numérique dû à la
présence de mailles mixtes contenant de l’air et des produits de détonation. Cette zone est
indiquée sur les différents graphiques.
11
Maille mixte : élément de discrétisation qui comprend plus d’un matériau ; nécessite un traitement numérique spécifique par le
code de calcul
125
H 1/k Rayon soit H1/k ≈ 0,1λ m et H1/1 = 2,85 m, charge de rayon (R1)
2 (144)
k
δ’épaisseur au sommet de la barrière de protection peut être estimée à partir de la
relation (145).
emin,1/k
0,5 soit e
min,1/k ≈ 0,0γ m pour k = 15 (145)
k
δa largeur minimale de la barrière de protection au sommet de la pile d’explosif peut être
évaluée à partir de la relation (146).
0,80 m x
1,00 m
1,20 m
1,40 m
1,70 m
1,90 m
z 2,10 m
Figure 107 : Dimensions du merlon 1A pour la série 1
0,80 m x
1,00 m
1,20 m
1,40 m
1,70 m
1,90 m
z 2,10 m
Figure 108 : Dimensions du merlon 2A pour la série 1
Le tableau 28 résume les configurations étudiées à petites échelles pour les merlons de la
série 1 (1A et 2A) avec deux charges gazeuses (R1 = 0,06 m et R2 = 0,03 m).
R (en m)
Numéro du cas d (en m) H (en m) e (en m) L (en m) Remarques
E (en MJ)
Merlon 1A
0,06 Etude du couplage des
1A_C6_d7 0,070 0,190 0,030 0,800
(E1 = 13,75.10-3) différents phénomènes
physiques (réflexion et
0,06
détente d’onde de choc)
Etude de l’effet protecteur
1A_C6_d8,5 0,085 0,190 0,030 0,800
(E1 = 13,75.10-3)
0,06 suivant les recommandations
1A_C6_d10 0,100 0,190 0,030 0,800
(E1 = 13,75.10-3) de l’OTAN
Etude de l’effet protecteur
0,03 pour une charge inférieure
1A_C3_d8,5 0,085 0,190 0,030 0,800
(E2 = 1,72.10-3) aux recommandations de
l’OTAN
Merlon 2A
Etude de l’effet protecteur
0,06 pour une charge supérieure
2A_C6_d8,5 0,085 0,160 0,060 0,800
(E1 = 13,75.10-3) aux recommandations de
l’OTAN
Etude de l’effet protecteur
0,03
2A_C3_d8,5 0,085 0,160 0,060 0,800 suivant les recommandations
(E2 = 1,72.10-3)
de l’OTAN
Tableau 28 : Tableau des configurations étudiées à petites échelles pour la série 1
placées en champ proche de la charge explosive (effet étudié et analysé dans le § 2.4). Cet
effet se traduit par une surestimation du maximum de surpression de l’ordre de 10% pour les
résultats obtenus en géométrie 2D-axisymétrique dans la zone d’intérêt située en aval, dite
« zone protégée » (pour une configuration avec d > 0,58 m/MJ1/3).
La série 1 étudie l’interaction d’une onde de choc avec une barrière de protection
placée en champ proche de la charge gazeuse (0,βλ ≤ d ≤ 0,71 m/MJ1/3). Les configurations,
pour lesquelles la distance d est inférieure à 0,58 m/MJ1/3, sont donc étudiées en géométrie
2D-axisymétrique et en géométrie 3D. Les calculs en géométrie 3D doivent permettre
d’évaluer la surestimation par un effet axisymétrique propre à l’approche en géométrie 2D
(étude pour d = 0,29 m/MJ1/3 (configuration 1A_C6_d7)). Si l’écart entre les deux approches
se situent dans un intervalle de confiance fixé à +/- 20%, les résultats obtenus en géométrie
2D-axisymétrique « fin » pourront être exploités pour la zone étudiée.
3.2.2.2 Réflexion sur la face avant des merlons 1A et 2A
δ’onde incidente sphérique divergente (charge hémisphérique placée au niveau du sol)
se réfléchit sur la face avant de la barrière de protection. Suivant la position et les dimensions
de la barrière (d, H, α1), le régime de réflexion peut être différent : réflexion régulière ou
réflexion de Mach (voir § 1.3.1). δ’analyse bibliographique (formule simplifiée de
Kinney (1985), relation (83)) sur les merlons 1A et 2A et pour une distance d comprise entre
0,29 et 0,70 m/MJ1/3, estime la formation d’un pied de Mach sur la face avant de la barrière.
Cette observation sur les régimes de réflexion est également confirmée lors de l’utilisation de
l’abaque des surpressions réfléchies du TM5-1300 (2008).
δes données expérimentales obtenues pour le merlon 1A placé à 0,07 m d’une charge
gazeuse de rayon (R1) (soit d = 0,29 m/MJ1/3) sont comparées aux résultats de simulation
numérique pour les deux approches géométriques utilisées (géométrie 2D-axisymétrique et
géométrie 3D) avec deux tailles de mailles différentes au passage du front de choc : maillage
« fin » et maillage « grossier [c] » (voir tableau 27).
y B
Gi3
d = 0,07 m
x
A
a) b)
Figure 109 : Evolution du maximum de surpression réfléchie sur la face avant du merlon 1A placé à 0,07 m
d’une charge gazeuse de rayon (R1) (1A_C6_d7, soit 0,29 m/MJ1/3)
y
OI
Interface entre l’air et Air
les produits de Produits de
détonation détonation
OIF OTA
PM
d = 0,07 m OTP
H = 0,19 m
Barrière
ORR
α1 = 45° x
a)Charge (R1) b) t = 0,037 ms
t = 0,0 ms
Phénomène de réflexion sur la face avant du merlon :
OI : onde incidente
PM : pied de Mach
OTP OTA OIF : onde issue de la focalisation des ondes de détente au
OI centre de la charge
Figure 110 : Visualisation sur les champs de pression pour les différents phénomènes de réflexion sur la face
avant de la barrière de protection placée en champ proche de la charge
δ’onde transmise dans les produits de détonation (OTP) est transmise à nouveau à l’air
(OTA). Cette onde rattrape l’onde incidente (OI) (écoulement subsonique en aval de l’onde
incidente (notée OI)) et se recombine avec l’onde réfléchie sur l’interface (ORR). δe couplage
de ces phénomènes physiques dépend des dimensions de la configuration étudiée (d, α1, H) et
de la nature de la charge explosive (vitesse de l’interface, fonction de la célérité de détonation
(DCJ) et du rayon de la charge de gaz (R)). Ce phénomène est étudié plus en détail dans la
série 2. Le digramme sur la figure 111 présente ces premières observations.
Figure 111 : Diagramme des phénomènes présents sur la face avant du merlon
y y
B
Gi3
B Gi3
y d = 0,10 m d = 0,085 m
B
Gi3 x
d = 0,085 m A
x A
x
A
Figure 112 : Evolution du maximum de surpression réfléchie sur la face avant du merlon 1A
132
La phénoménologie présente sur la face avant du merlon 1A est également observée pour le
merlon 2A placé à 0,08η m d’une charge de gaz de rayon (R1), figure 113-a :
effet axisymétrique en géométrie 2D (effet de compression sur la paroi),
premier front de choc confondu avec une onde réfléchie pour les résultats en géométrie
3D, convergence numérique non atteinte (effet compensatoire), mais ces résultats se
situent dans l’intervalle de confiance de +/- 20%,
l’incertitude expérimentale est légèrement plus importante (de +/- 14% à +/- 12%),
augmentation du maximum de surpression réfléchie sur la paroi par la déformation de
l’onde réfléchie avec les produits de détonation (écart entre les données expérimentales-
les résultats de simulation et les estimations du TM5-1300 (2008)).
δe merlon βA situé à 0,08η m d’une charge gazeuse de rayon (R2), soit 0,71 m/MJ1/3, est
étudié en géométrie 2D-axisymétrique « fin ». Les résultats de simulation numérique se
situent dans l’intervalle d’incertitude expérimentale de +/- 12%, figure 113-b.
Merlon 2A placé à 0,085 m de la charge (R1) (2A_C6_d8,5) : calcul en géométrie 3D nécessaire
y B
Gi3
d = 0,085 m
A
x
a) b)
Merlon 2A placé à 0,085 m de la charge (R2) (2A_C3_d8,5) : calcul en géométrie 2D possible
y B
Gi3
d = 0,085 m
x
A
c) d)
Figure 113 : Evolution du maximum de surpression réfléchie sur la face avant du merlon 2A
133
Maximum de
surpression (bar)
Pied de Mach 0,957
Détente 0,665
θ = α1
0,374
t = 0,50 ms -0,083
α1 = 45° -0,208
Figure 114 : Visualisation sur les champs de pression lors du phénomène de détente sur le sommet de la barrière
de protection (configuration 1A_C6_d7)
y Gi4
A B
x
a) b)
Figure 115 : Evolution du maximum de surpression sur le sommet du merlon 1A placé à 0,07 m d’une charge
gazeuse de rayon (R1) (1A_C6_d7)
Pour les géométries de barrière étudiées, un pied de Mach se forme sur la face avant du
merlon. Compte tenu de la distance parcourue par le pied de Mach, il est possible de supposer
que ce front se propage perpendiculairement à la paroi, ainsi l’angle de déviation de la
première détente est égal à l’angle d’inclinaison de la paroi ( = α1).
L’effet axisymétrique en géométrie 2D-axisymétrique est faible pour cette zone
d’intérêt (+ 10% pour le capteur Gi4, figure 115-b). Un calcul en géométrie 3D « fin » aurait
été préférable, mais au vu des incertitudes expérimentales (+/- 10%), l’approche en géométrie
134
2D-axisymétrique « fin » est suffisante et acceptable pour assurer la restitution des données
expérimentales. Cette hypothèse sera appliquée aux autres barrières placées à plus de
0,29 m/MJ1/3 du centre de la charge. Les résultats obtenus en géométrie 2D-axisymétrique
« fin » seront donc comparées aux données expérimentales. Le diagramme sur la figure 116
synthétise ces différents points d’analyse.
Phénomènes de réflexion
Figure 116 : Diagramme des phénomènes physiques présents sur le sommet du merlon
a) Merlon 1A placé à 0,085 m de la charge (R1) b) Merlon 1A placé à 0,10 m de la charge (R1) c) Evolution de la surpression maximale au sommet
(1A_C6_d8,5) (1A_C6_d10) des merlons 1A pour une charge gazeuse (R1)
y Gi4 y Gi4
A B A B
x x
Pour cette zone d’intérêt, la différence entre les deux obstacles (merlons 1A et 2A) se situe sur
l’épaisseur du merlon (épaisseur plus importante pour le merlon 2A, e2A = e1A 2, l’effet
d’atténuation par la distance parcourue est donc plus important), figure 118.
y Gi4
A B
x
a) b)
y Gi4
A B
x
c) d)
Maximum de
surpression (bar)
0,378
θ = α1
Détente
θ = α2 0,210
0,042
Figure 119 : Visualisation sur les champs de pression pour le phénomène de détente sur la face aval de la
barrière de protection (configuration 1A_C6_d7)
y A Gi5
B x
a) b)
Figure 120 : Evolution du maximum de surpression sur la face aval du merlon 1A placé à 0,07 m d’une charge
gazeuse de rayon (R1) (1A_C6_d7)
n’est donc pas nécessaire. Les valeurs du maximum de surpression en champ libre sont plus
importantes que les valeurs sur la face aval du merlon (ΔP+champ libre >> ΔP+barrière). Ces valeurs
en champ libre ne sont donc plus présentées sur les différents graphiques
δ’onde de choc issue du sommet se réfléchit ensuite sur le sol en aval du merlon. La
présence de ce phénomène conduit à une augmentation du maximum de surpression sur le
pied de la face aval. δ’impact de ce phénomène dépend des caractéristiques géométriques de
la barrière de protection (angle d’incidence au niveau du sol ( = f (α2, H)) et de l’impact des
phénomènes précédents (nature de l’onde (réflexion sur la face avant) et de l’atténuation par
la présence des phénomènes de détente). Le diagramme de la figure 121 synthétise les
observations réalisées sur la face aval du merlon.
a) Merlon 1A placé à 0,085 m de la charge (R1) b) Merlon 1A placé à 0,10 m de la charge (R1) c) Merlon 1A placé à 0,085 m de la charge (R1)
(1A_C6_d8,5) (1A_C6_d10) (1A_C3_d8,5)
y A y y
Gi5 A Gi5 A Gi5
B x B x B x
Le couplage des différents phénomènes physiques présents sur la face arrière du merlon 1A
est également observé sur le merlon 2A, figure 123 :
y A Gi5
B x
a) b)
y A Gi5
B x
c) d)
Pour le merlon 2A placé face à une charge gazeuse de rayon (R1), l’historique de pression
pour le capteur Gi5 est parfaitement restitué par l’outil de simulation (configuration
2A_C6_d8,5, figure 123-b).
Pour la seconde configuration du merlon 2A, avec une charge gazeuse de rayon (R2), l’écart
de 15% entre les données expérimentales et les résultats de simulation numérique est
légèrement plus important que l’incertitude expérimentale (incertitude expérimentale mesurée
141
de +/- 12%, configuration 2A_C3_d8,5, figure 123-d). Cet écart supérieur à l’incertitude
expérimentale peut être en partie dû à la qualité de la mesure qu’à la qualité de la simulation
au vu de la faible dynamique signal/bruit de l’ordre de 1γ% pour ce point de mesure.
3.2.2.5 Effet protecteur des merlons 1A et 2A
δa zone en aval de la barrière de protection correspond à la zone d’intérêt du projet
ANR BARPPRO, dite « zone protégée ». Les résultats présentés précédemment ont permis de
confirmer les capacités prédictives de l’outil de simulation numérique avec une précision
maitrisée : restitution de la mesure par la simulation numérique correcte par rapport aux
incertitudes expérimentales et à la qualité de la mesure. Le couplage des différents
phénomènes physiques présents en fonction des caractéristiques géométriques de la barrière a
également pu être analysé, confirmant ainsi les hypothèses de travail évoquées lors de l’étude
de configurations de barrière de protection issues de la littérature (§ 1.3.3.1, figure 41).
t = 1,53 ms -0,032
-0,113
-0,195
Figure 124 : Visualisation sur les champs de pression pour le phénomène de réflexion en aval de la barrière de
protection (configuration 1A_C6_d7)
δ’onde qui passe par le sommet de la barrière de protection se réfléchit sur le sol en
aval du merlon, figure 124. Ce phénomène physique conduit à une augmentation du
maximum de surpression en aval du merlon.
a) b)
« grossier [c] ») que celui entre les géométries. La géométrie 2D-axisymétrique « fin » est
donc suffisante pour étudier cette dernière zone d’intérêt, et des calculs hypothétiques en
géométrie 3D « fin », n’apporteraient pas d’informations supplémentaires.
Dans la zone au sol derrière le merlon, une instrumentation riche de 7 capteurs est disposée
sur le banc d’expérimentation. Ces mesures fournissent une évolution détaillée de la
décroissance du maximum de surpression, GH1 à GH7, figure 125-b. Un bon accord
simulation-mesure est ainsi constaté sur les différents capteurs en tenant compte de
l’incertitude expérimentale liée à l’essai réalisé pour cette configuration (+/- 10%). Une seule
mesure se détache de la simulation numérique, GH4, avec un écart de 18%. Les restitutions
des mesures étant correctes à l’amont du capteur GHζ (γ capteurs) et à son aval (γ capteurs
également) et que dans cette zone le régime d’écoulement ne varie plus, il est possible
d’émettre l’hypothèse que c’est la mesure qui est en défaut pour le capteur GHζ, ou que les
résultats de simulation numérique obtenus en géométrie 2D-axisymétrique « fin » sont
affectés d’une surestimation « locale » induite par l’approche (merlon torique, recombinaison
du premier front de choc avec l’onde issue d’une recombinaison sur l’axe de révolution).
Les longueurs des merlons 1A et βA étant telles qu’ils sont assimilables à des merlons
infinis, les distances parcourues par l’onde pour se réfléchir à leur aval sont très proches
(S1A = 0,527 met S2A = 0,513 m), figure 126. Les merlons 1A et 2A se différencient par leur
hauteur (H1A = 0,19 m et H2A = 0,1θ m) et leur épaisseur au sommet de l’obstacle
(e1A = 0,03 m et e2A = 0,06 m).
H1A = 0,19 m
H2A = 0,16 m
α1 = 45° α2 = 45° α1 = 45° α2 = 45°
x x
Distance parcourue par l’onde issue du sommet pour Distance parcourue par l’onde issue du sommet pour
le merlon 1A (H1A = 0,19 m) : le merlon 2A (H2A = 0,16 m) :
S1A = t1A + e1A + b1A = 0,527 m S2A = t2A + e2A + b2A = 0,513 m
Figure 126 μ Représentation schématique des distances parcourues par l’onde pour les merlons de la série 1
δe facteur d’atténuation permet d’évaluer l’effet protecteur du merlon par rapport à une
configuration champ libre (sans structure), relation (147).
Pr
Ap
Pi
si Ap → 0, effet protecteur maximale (147)
La figure 127 présente l’évolution du facteur d’atténuation pour les deux géométries
analysées (merlons 1A et 2A) pour les deux charges étudiées (R1 ; R2). Il est possible de
distinguer deux zones singulières : zone I et zone II.
Dans la zone I, les résultats de simulation numérique confirment que plus la charge est
importante devant la barrière de protection, plus l’effet protecteur diminue en aval du
merlon (observation identique sur le mur de Borgers (2010), § 1.3.3.2). La barrière de
protection « idéale » est logiquement une barrière avec une hauteur et une épaisseur
importantes afin d’accentuer l’atténuation par la distance parcourue. En effet, l’énergie
relâchée par la détonation d’une charge de gaz de rayon (R1) est plus élevée que pour une
charge de rayon (R2) (E1 = 1,375.10-2 MJ > E2 = 1,719.10-3 MJ) (distance réduite ( )).
Dans la zone II, le dimensionnement du merlon 1A et 2A, suivant les recommandations
de l’OTAN (sans tenir compte des ondes de contournement), permet d’assurer un effet
protecteur jusqu’à quatre fois la hauteur du merlon (zone II, le facteur d’atténuation est
deux fois plus important pour une charge de Rayon (R2) par rapport à celui d’une charge
de rayon (R1) avec le même merlon (1A ou 2A)).
δe dimensionnement d’une barrière de protection dépend donc de la masse de la charge
explosive (W) (abaques en échelle réduite) et de la position de la charge par rapport à la
structure (distance d).
II
Figure 127 μ Evolution du facteur d’atténuation en aval du merlon pour les deux charges gazeuses étudiées
II I
II
Figure 128 μ Evolution du facteur d’atténuation en aval du merlon pour plusieurs positions du merlon 1A
permettent d’analyser l’influence du rapport de forme du merlon pour une même distance
parcourue.
δa variation du rayon de la charge gazeuse permet d’analyser l’influence de la charge
explosive. A merlon donné, plus la charge explosive est importante, moins l’effet
protecteur est important en aval de la barrière (observation zone I, figure 127).
A charge donnée, et pour les deux géométries étudiées, l’effet protecteur obtenu est
quasiment identique (écarts négligeables entre les deux merlons). Les recommandations
de l’OTAN sont donc robustes. Pour une même distance parcourue, les variations
(faibles) sur la forme de la barrière de protection n’induisent pas de différences
significatives sur l’effet protecteur en aval.
δa présence d’une épaisseur au sommet du merlon (e ≠ 0) permet de détendre « en deux
étapes » le pied de Mach issu de la réflexion sur la face avant. Cette recommandation de
l’OTAN est une « bonne pratique » dans le dimensionnement d’une barrière de protection
(e > 0,5 m).
δ’écart entre le merlon 1A et βA est faible en champ proche de la structure (différence de
2%). δ’écart sur le facteur d’atténuation en champ lointain (4 H) est d’environ - 5%
A cette distance (≈ 4 H), pour les merlons placés face à une charge gazeuse de rayon
pour le merlon 1A par rapport au merlon 2A (pour les deux rayons de charge (R1 et R2)).
(R1), l’effet protecteur est très faible. Au-delà de cette distance, les merlons ne protègent
plus suffisamment, l’effet protecteur est donc considéré comme « terminé ».
Figure 129 μ Evolution du facteur d’atténuation pour les deux géométries de merlon étudiées
146
0,80 m x
1,00 m
1,20 m
1,40 m
1,70 m
z 1,90 m
Afin d’évaluer l’impact de l’inclinaison des faces amont et aval sur l’effet protecteur,
le merlon 2B est construit avec un angle d’inclinaison de λ0° pour la face avant et de ζη° pour
la face arrière, figure 131. La comparaison des résultats obtenus pour ces deux barrières de
147
Gi3 Gi4
0,063 m Gi2 Gi5 0,09 m
Charge H = 0,19 m
(0,0) Gi1 Gi6
α1 = 90° α2 = 45° x
0,80 m x
1,00 m
1,20 m
1,40 m
1,70 m
z 1,90 m
Le tableau 29 synthétise les configurations étudiées à petites échelles avec les merlons 1B et
2B pour une charge de gaz (R1 = 0,06 m).
R (en m)
Numéro du cas d (en m) H (en m) e (en m) L (en m) Remarques
E (en MJ)
Merlon 1B (α1 = 45° ; α2 = 90°)
y Maximum de
surpression (bar)
t = 0,12 ms 13,80
x -0,01
y Maximum de
surpression (bar)
10,70
t = 0,14 ms
8,01
Phénomène de réflexion perturbé par
l’onde issue d’une recombinaison sur 5,34
Figure 132 : Visualisation sur les champs de pression pour le merlon 2B en géométrie 2D-axisymétrique
(configuration 2B_C6_d7)
y B
Gi2
d = 0,07 m
x
A
a) b)
Figure 133 : Evolution du maximum de surpression réfléchie sur la face avant du merlon 2B placé à 0,07 m
d’une charge gazeuse de rayon (R1) (2B_C6_d7, soit 0,29 m/MJ1/3)
Simulation en géométrie 3D
y y OI PM
OI
OTA PMR
OR
Interface ORR
Interface
ORR
OTP
(0, 0) (0, 0)
t = 0,05 ms
t = 0,10 ms
a) x b) x
Figure 134 : Visualisation sur les champs de pression pour le merlon βB placé à 0,07 m d’une charge gazeuse
(propane-oxygène à la stœchiométrie, rayon (R1))
familles d’explosif est très différente (DCJ du TNT ≈ β,λ DCJ du gaz). Les valeurs de la
charge chimique condensée (TNT) simulée uniquement. La célérité de détonation de ces deux
L’écart entre la distance parcourue par l’onde et l’interface air/produits de détonation au cours
du temps de propagation de l’onde choc peut être analysé pour ces deux natures de charge
explosive à énergie équivalente. Les simulations numériques sont réalisées en géométrie
2D-axisymétrique avec la mise en œuvre d’un maillage « fin » de 0,8 mm/εJ 1/3. Les
visualisations des champs de pression et de l’interface air/produits de détonation pour les
natures de charge explosive étudiées sont présentées sur la figure 135.
151
Pour un temps de propagation identique à énergie équivalente, l’onde de choc issue d’une
charge chimique condensée possède une célérité plus importante qu’une onde de choc
aérienne émise par la détonation d’une charge gazeuse. Cet écart important provient de la
différence de célérité de détonation entre les deux natures de charge.
La valeur de la célérité de détonation de la charge explosive affecte également la vitesse
de l’interface entre l’air et les produits de détonation. Ainsi la vitesse de l’interface est
plus importante dans le cas d’une charge chimique que dans le cas d’une charge de gaz.
Cette observation est également confirmée par les relations (57) et (58) pour une particule
des produits de détonation située à la limite de la charge explosive.
Des instabilités se forment le long de l’interface lors du mélange d’un milieu à haute
impédance, et d’un milieu à basse impédance (instabilités de type Rayleigh-Taylor). Dans
ce cas de figure, un déplacement moyen de l’interface est calculé pour la charge de TNT.
90° Air
30° 60°
Produits de
détonation
30°
Produits de
détonation
0°
x
0° 1/3
0,48 m/MJ
Charge initiale de TNT x Charge initiale de Gaz
1/3 1/3
R = 0,03 m/MJ1/3 0,36 m/MJ 0,66 m/MJ
R = 0,25 m/MJ1/3
Figure 135 : Visualisation sur les champs de pression et du déplacement de l’interface entre l’air et les produits
de détonation pour une charge de TNT et une charge de gaz (propane-oxygène à la stœchiométrie)
δe calcul du temps d’arrivée de l’onde de choc pour une charge de TNT peut être estimé, en
fonction de la distance réduite (z en m/kg1/3 de TNT), à partir de la vitesse du son dans l’air
(milieu inerte et non choqué, indice 0 (c0 ≈ γγ1,η m/s)) et de la variation du maximum de
surpression (formule (θ)) pour un gaz parfait ( = 1,ζ), relation (1ζλ) (Kinney (1λ85)).
c0 Rc
t0
1 Pi
1 z 1
1
(149)
2 P0
Avec Rc, rayon de la charge de TNT, ∆Pi+, maximum de surpression incident, P0, pression
initial dans le milieu inerte
152
Figure 136 μ Comparaisons des positions du front de choc et de l’interface entre l’air et les produits de détonation
pour une charge de TNT et une charge de gaz (propane-oxygène à la stœchiométrie)
a) Merlon 2B placé à 0,085 m de la charge (R1) b) Merlon 2B placé à 0,10 m de la charge (R1) c) Merlon 1A placé à 0,14 m de la charge (R1)
(2B_C6_d8,5) (2B_C6_d10) (2B_C6_d14)
y y Gi2 B
B
Gi2 d = 0,14 m
x
y d = 0,10 m A
B
x
Gi2 A
d = 0,085 m
x
A
Figure 137 : Evolution du maximum de surpression réfléchie sur la face avant du merlon 2B
154
1,240
0,483
-0,277
α1 = 90° α2 = 45° t = 0,82 ms
Figure 138 : Visualisation sur les champs de pression lors du phénomène de réflexion sur la face avant de la
barrière de protection (configuration 2B_C6_d10)
155
Maximum de
surpression (bar)
0,717
θ = α1
0,482
0,247
0,013
-0,222
α1 = 45° α2 = 90°
t = 0,70 ms
a) Configuration 1B_C6_d14
Maximum de
surpression (bar)
0,717
θ = α1
0,482
0,247
0,013
-0,222
t = 0,40 ms
α1 = 90° α2 = 45°
b) Configuration 2B_C6_d14
Figure 139 : Visualisation sur les champs de pression du phénomène de détente sur le sommet de la barrière de
protection (configurations 1B_C6_d14 et 2B_C6_d14)
Pour le merlon 1B, les données expérimentales sont comparées aux résultats de
simulation numérique en géométrie 2D-axisymétrique « fin ». δ’approche en géométrie
2D-axisymétrique peut induire un faible effet de compression (effet axisymétrique, faible
écart pour cette zone d’intérêt (voir § γ.β.β.γ)), mais cette approche permet la mise en œuvre
d’un maillage « fin » qui assure un niveau de convergence numérique satisfaisant. Cette
approche est donc préférable car elle permet la meilleure restitution possible des mesures,
figure 140.
La phénoménologie présente sur le sommet du merlon 1B est identique aux merlons 1A et
2A. Le niveau de surpression diminue rapidement au début de la paroi du fait de la détente
induite par la rupture de pente au sommet. La géométrie du merlon 1B se différencie des
merlons de la série 1 par une épaisseur au sommet plus importante (e = H) permettant ainsi
d’accentuer l’atténuation par la distance parcourue sur cette surface.
156
y Gi4
A
B
x
a) b)
Merlon 1B placé à 0,14 m de la charge (R1) (1B_C6_d14)
y Gi4
A
B
x
c) d)
y Gi3
A B
x
a) b)
Figure 141 μ Evolution du maximum de surpression sur le sommet du merlon βB placé à 0,07 m d’une charge
gazeuse de rayon (R1) (2B_C6_d7)
Les surpressions maximales obtenues sur le sommet du merlon 2B pour les trois autres
configurations (distance charge – pied amont du merlon, d = 0,085 ; 0,100 ; 0,140 m) sont
présentées sur la figure 142 :
Les résultats obtenus en géométrie 3D ne sont pas convergés numériquement avec le
maillage « grossier [c] » (sous-estimation de 20-βη%). Néanmoins, l’outil de simulation
restitue qualitativement les phénomènes physiques présents et quantitativement les
mesures, environ, dans l’intervalle de confiance de +/- 20%.
Le maximum de surpression est atténué par la présence d’une détente de l’onde de choc
entre la face avant et le sommet de la barrière de protection. δors de la formation d’un
pied de εach sur la face avant du merlon, l’angle de déviation est plus important dans le
cas d’une paroi fortement inclinée ( = α1 = 90°, dimensionnement du merlon 2B).
Sur le sommet de la barrière de protection, le maximum de surpression est également
atténué par la distance parcourue par l’onde sur cette surface.
δors de la série β, les géométries étudiées se différencient par l’inclinaison de la face avant
(α1 = ζη° pour le merlon 1B et α1 = 90° pour le merlon 2B). Cette différence de
dimensionnement modifie la nature et l’intensité de l’onde de choc incidente au sommet du
merlon, ce qui affecte également le phénomène de détente entre la face avant et le sommet du
merlon μ paroi fortement inclinée → maximum de surpression réfléchie plus important sur la
paroi → angle de déviation et nombre de Mach incident important sur le sommet.
158
a) Merlon 2B placé à 0,085 m de la charge (R1) b) Merlon 2B placé à 0,10 m de la charge (R1) c) Merlon 1A placé à 0,14 m de la charge (R1)
(2B_C6_d8,5) (2B_C6_d10) (2B_C6_d14)
Maximum de
surpression (bar)
0,225
0,126
0,027
2
2 -0.072
Onde incidente
α1 = 45° α2 = 90° -0,171
sur le sol
t = 1,20 ms
a) Configuration 1B_C6_d14
Maximum de
surpression (bar)
0,225
0,126
0,027
-0.072
Onde incidente sur le sol -0,171
2
t = 0,90 ms
α1 = 90° α2 = 45° 2
b) Configuration 2B_C6_d14
Figure 143 : Visualisation sur les champs de pression du phénomène de détente sur la face aval de la barrière de
protection (configurations 1B_C6_d14 et 2B_C6_d14)
a) b)
Merlon 1B placé à 0,14 m de la charge (R1) (1B_C6_d14)
c) d)
Après avoir étudié la « forte » détente le long de la face aval du merlon 1B, une
détente plus « douce » est analysée sur le merlon 2B. En effet pour le merlon 2B, la face
arrière est inclinée à 45° ( = ζη°), l’angle de déviation au sommet de la face arrière est donc
plus faible que pour le cas du merlon 1B ( = λ0°). Le long de la face arrière du merlon, la
surpression maximale subit donc une atténuation plus faible par rapport à celle du merlon 1B.
δ’angle d’incidence de l’onde sur le sol au niveau du pied aval du merlon est égale à ζη°, le
régime de réflexion dépend donc de l’évolution du nombre de εach pour une transition
éventuelle vers une réflexion de εach. δe passage de l’onde réfléchie sur la face aval est
moins marqué que dans le cas du merlon 1B. Dans la zone dite « protégée », l’onde issue du
sommet du merlon subit une réflexion régulière oblique au niveau du pied aval de la barrière.
δ’intensité de l’onde réfléchie est plus faible que dans le cas d’une réflexion normale (cas du
merlon 1B).
Pour cette zone d’intérêt, les valeurs du maximum de surpression en champ libre (sans
obstacle) sont plus importantes que les valeurs sur la face aval (ΔP+champ libre >> ΔP+barrière).
161
y A
Gi6
B x
a) b)
Figure 145 μ Evolution du maximum de surpression sur la face aval du merlon βB placé à 0,07 m d’une charge
gazeuse de rayon (R1) (2B_C6_d7)
Les figures 145 et 146 présentent l’évolution du maximum de surpression sur la face
aval du merlon 2B pour plusieurs valeurs de la distance d (0,070 ; 0,085 ; 0,100 et 0,140 m).
Les résultats obtenus en géométrie 3D sous-estiment d’environ 1η-20% la valeur du
maximum de surpression issue des données expérimentales. Cette sous-estimation est
induite par l’utilisation d’un maillage « grossier [c] » (sous-estimation maximale attendue
de l’ordre de β0%).
δ’onde subit une seconde détente lors de son passage sur la face aval du merlon. δa
présence de ce phénomène permet d’atténuer à nouveau le maximum de surpression.
Cette seconde atténuation contribue à l’apparition d’un effet protecteur. δ’angle de
déviation est égal à l’inclinaison de la face aval du merlon ( = α2).
δ’onde se réfléchit sur le sol en aval du merlon. δa présence de ce phénomène conduit à
une recompression de l’onde qui se traduit par une augmentation de la surpression au
niveau du pied aval du merlon. Néanmoins, l’impact de ce phénomène reste limité en
comparaison du merlon 1B par un angle d’incidence plus important que dans le cas du
premier merlon.
δe dimensionnement de la face aval du merlon conditionne la valeur de l’angle de déviation
de la seconde détente à son aval et l’angle d’incidence du phénomène de réflexion sur le sol
dans la zone « protégée ». Ces phénomènes physiques dépendent également de l’évolution du
nombre de εach de l’onde de choc et donc de l’impact des phénomènes qui les précédent :
réflexion sur la face avant de la barrière, détente sur le sommet,…
162
a) Merlon 2B placé à 0,085 m de la charge (R1) b) Merlon 2B placé à 0,10 m de la charge (R1) c) Merlon 1A placé à 0,14 m de la charge( R1)
(2B_C6_d8,5) (2B_C6_d10) (2B_C6_d14)
y y A y
A A
Gi6 Gi6
Gi6 x
x B B x
B
x x
Distance parcourue par l’onde issue du sommet pour Distance parcourue par l’onde issue du sommet pour
le merlon 1B (H1B = 0,19 m) : le merlon 2B (H2B = 0,19 m) :
S1B = t1B + e1B + b1B = 0,649 m S2B = t2B + e2B + b2B = 0,649 m
Figure 147 : Représentation schématique des distances parcourues par l’onde pour les merlons de la série β
Les deux géométries de merlon étudiées se différencient par leur angle d’inclinaison des faces
avant et arrière (merlon 1B μ α1 = ζη° et α2 = 90° ; merlon 2B μ α1 = λ0° et α2 = 45°). Cette
différence affecte donc les valeurs des angles d’incidence (réflexion en amont et en aval du
merlon) et des angles de déviation (face avant-sommet et sommet-face arrière).
Le cas du merlon 1B peut être modélisé en géométrie 2D-axisymétrique avec une face
avant inclinée à ζη° qui limite suffisamment l’effet de confinement autour de l’axe de
révolution. Avec cette approche géométrique, la mise en œuvre d’un maillage « fin », défini
par le standard de calcul, permet d’obtenir des résultats avec un niveau de convergence
numérique suffisant. Les résultats de simulation numérique pour les capteurs au sol, et situés à
l’aval des barrières de protection (GH1 à GH7), se situent alors dans l’intervalle d’incertitude
expérimentale (+/- 10%), figures 148-a et 148-b.
Le merlon 2B ne peut pas être modélisé en géométrie 2D-axisymétrique, les
simulations sont donc réalisées en géométrie 3D. La taille de mailles est dégradée afin de
rendre compatible ce cas de simulation avec les moyens de calcul. Les simulations sont ainsi
réalisées avec le maillage « grossier [c] ». δ’utilisation de ces conditions induit une
sous-estimation des résultats d’environ 20% par une convergence numérique non atteinte. Les
résultats de simulation en géométrie 3D ne permettent pas une restitution complètement
satisfaisante des mesures dans l’intervalle d’incertitude expérimentale (+/- 10%), mais ils sont
qualitativement corrects (décroissance des profils) par rapport aux données expérimentales
(écart proche de l’intervalle de confiance de +/- 20%), figures 148-c et 148-d.
164
a) Merlon 1B placé à 0,07 m de la charge (R1) b) Merlon 1B placé à 0,14 m de la charge (R1)
(1B_C6_d7) (1B_C6_d14)
c) Merlon 2B placé à 0,07 m de la charge (R1) d) Merlon 2B placé à 0,14 m de la charge (R1)
(2B_C6_d7) (2B_C6_d14)
Figure 148 : Evolution du maximum de surpression en aval des merlons 1B et βB placés à 0,07 et 0,1ζ m d’une
charge de gaz de rayon (R1)
Figure 149 μ Evolution du facteur d’atténuation pour les merlons de la série β et le merlon 1A de la série 1 placés
à 0,07 m d’une charge de gaz de rayon (R1)
δa distance parcourue par l’onde le long des merlons 1B et βB est identique, les deux
géométries se différencient par l’angle d’inclinaison de leurs parois (merlon 1B μ α1-1B = 45°
et α2-1B = 90°, merlon 2B μ α1-2B = λ0° et α2-2B = 45°). Cette différence sur le dimensionnement
des merlons affecte les valeurs des angles de déviation pour les deux détentes présentes lors
de l’interaction de l’onde de choc avec la barrière de protection, figure 150.
Le pied de Mach issu du phénomène de réflexion sur la face avant du merlon se détend
sur le sommet de la barrière μ merlon 1B, 1-1B = α1-1B = 45° et merlon 2B,
1-2B = α1-2B = 90°. δ’intensité du pied de Mach est également différente entre les deux
merlons. Le pied de Mach issu du merlon 1B sera plus faible que le pied de Mach issu du
166
Figure 150 : Evolution du facteur d’atténuation pour le merlon 1B avec plusieurs valeurs de la distance d (0,070
et 0,1ζ0 m d’une charge de gaz de rayon (R1))
Les figures 150 et 151 permettent de visualiser l’évolution du facteur d’atténuation pour les
géométries de merlon de la série 2 pour plusieurs valeurs de la distance d (0,07 ; 0,085 ;
0,100 ; 0,140 m du centre d’un charge gazeuse de rayon (R1 = 0,03 m)).
167
Pour le merlon 1B, il est possible de constater comme pour la série 1 : « une barrière de
protection placée en champ proche de la charge permet d’offrir un effet protecteur plus
important au pied aval du merlon » (voir § 3.2.2.5). Ces deux géométries de merlon
possèdent en effet la même inclinaison au niveau de leur face avant (α1 = 45°).
Pour le merlon 2B, le constat est différent. Cette géométrie de merlon se différencie du
merlon 1B par une face avant inclinée à 90° et une face arrière inclinée à 45°. La
variation de la distance d pour cette géométrie permet de constater que plus la barrière est
en champ proche de la charge, plus l’effet protecteur est important sur l’ensemble de la
zone « protégée » (champ proche et champ lointain).
Figure 151 μ Evolution du facteur d’atténuation pour le merlon 2B avec plusieurs valeurs de la distance d (0,070 ;
0,08η ; 0,100 ; 0,1ζ0 m d’une charge de gaz de rayon (R 1))
δ’étude de ces géométries permet de conclure que l’inclinaison de la face avant du merlon
peut influencer également le dimensionnement de la distance d en fonction de l’évolution et
de l’intensité de l’effet protecteur souhaité μ α1 important → d faible ; α1 faible → distance d
dimensionnée suivant la position de la zone à protéger.
168
Dimensionnement idéal : barrière parallélépipédique avec une hauteur et une épaisseur importantes
Figure 152 : Synthèse des observations issues de l’étude de différentes géométries de barrière de protection
171
l’évolution du maximum de surpression d’une onde de choc (∆P+). Les deux prochains
paragraphes décrivent les approches utilisées pour les phénomènes de réflexion et de détente
d’une onde de choc.
y y
RR 1 α1
Point S 2
Point S
RR RM = α1
α1 x α1 x
Charge
Capteurs
h = 1,85 m/kg1/3
b
x
Le domaine total calculé est dimensionné avec une longueur de 8 m. Les conditions de calculs
sont présentées dans le tableau 30. Pour cette étude en champ très proche de la charge
(b ≈ 0,γζη m), la taille des mailles est fixée à 0,5 mm au passage du front de choc jusqu’à
2,4 ms, puis à 1 mm jusqu’à la fin du calcul.
Taille de mailles au Temps de transition Nouvelle taille de mailles
Approche géométrique
début du calcul du maillage jusqu’à la fin du calcul
Géométrie 2D-axisymétrique 0,5 mm 2,4 ms 1 mm
Tableau 30 μ Conditions de calcul de l’étude du phénomène de détente
Figure 155 : Abaque des coefficients de détente à partir de résultats de simulation numérique pour une structure
d’angle d’inclinaison de λ0°
y z
z
Charge (W) x
x
Γ1 Γ
Secteur 4 Γ2
Γ3 (x,z)
Secteur 1 Secteur 3 Γ4
« Champ libre » Secteur 2
Figure 156 : Définition des secteurs géométriques autour de la barrière
177
Le tableau 31 décrit les phénomènes physiques présents dans les différents secteurs étudiés.
Comme évoqué précédemment l’onde de contournement qui se détend sur les coins de la
barrière ajoute deux nouveaux secteurs sur chaque bord du merlon.
Onde de choc « Champ libre » Secteur 1 Secteur 2 Secteur 3 Secteur 4
Réflexion sur la
Pas
face avant + détente Réflexion sur la face avant
Onde issue du d’interaction
sur la face latérale + + détente en deux étapes
merlon directe avec la
réflexion sur le sol + réflexion sur le sol en aval du merlon
barrière
en aval du merlon
Deux détentes de
l’onde de choc
Onde de
Une détente de l’onde sur le coin de la barrière entre les faces
contournement
avant et arrière
du merlon
Tableau 31 : Description des phénomènes physiques présents en fonction des secteurs
Onde de contournement
δ’estimation du maximum de surpression pour l’onde de contournement nécessite de
déterminer le régime de réflexion sur le coin de la face avant (RR ou RM) et les coefficients
de détente (Cdb1 et Cdb2), figure 157. Pour un point au sol situé en (x, z), il faut déterminer le
secteur d’appartenance. Pour les secteurs 1, 2, 3 et 4-bis, l’onde de contournement est soumise
à une unique détente sur le bord A. La diminution du maximum de surpression due à ce
phénomène peut être estimée à partir d’un coefficient de détente (Cdb1). Ce coefficient est
obtenu à partir de l’abaque des coefficients de détente exprimé en fonction de l’angle de
déviation ( b1) et du nombre de Mach en amont du phénomène (Md1). Pour cette première
détente, le nombre de Mach en amont dépend de la surpression réfléchie sur le bord de la face
avant (CrbS).
178
1 1
xd H e
2 )
L/2
tan( ) tan(
d 1
Γ4 x
(x,z)
z
b2
A L
Détente A, Cdb1 Détente B, Cdb2 2 Secteur 4
z
L
b1 Γ1 2 Secteur 4-bis
= Γ1
xd
z d (x,z)
Secteurs 1-2-3
L
A d 2
2
2
(150)
δ’angle d’incidence sur le coin de la barrière est noté d ( d = Γ1). Le coefficient de réflexion,
noté CrbS, est obtenu à partir de l’abaque du Tεη-1300 (2008). δ’angle de déviation est
évalué en fonction du régime de réflexion au point A, formules (151) et (152).
x d pour une réflexion de Mach
b1, MR arctan
z L/2
(151)
Pour un point d’observation placé dans le secteur ζ, l’onde de contournement est soumise à
une seconde détente centrée sur le point B. La surpression incidente au point B est estimée à
partir de la formule (139) en fonction de la distance parcourue par l’onde entre le centre de la
charge et le point B ( B), et des coefficients précédents (Cdb1 et CrbS dans le cas d’une
réflexion de Mach), formule (153).
1
B A e H
tan(1 ) tan( 2 )
1
(153)
δ’angle de déviation ( b2) pour cette seconde détente est obtenu à partir de la formule (154).
L/2 z
b2 arctan
1
(154)
x d e H
tan( ) tan( )
1
1 2
Point S x
x
Point S’ H
Γ
α1 α1
αS’ S αS
S’ z Γ z
Figure 158 : Points d’intersection sur les arêtes de la face avant pour l’onde issue du merlon
x Point de sortie
Point d’observation
sol
(x, z)
h = f(Γ) sol
x - xS
Point S’
sol
S’ 2 sol
z
L
2
z
Secteur 1 k = f(Γ)
Figure 159 μ Détente de l’onde de choc pour le secteur 1 et sa réflexion sur le sol
Point T’ x z
ξ' ξ
sol
ω' x Point T S
x
L ω
2 tan()
H S
H
z sol
z
L
a) Secteur 2 2 b) Secteurs 3 - 4
A titre d’exemple, deux points de sortie sont possibles sur la face avant du merlon,
figure 158. Lorsque la ligne brisée sort de la face avant du merlon, le point de sortie est noté
S’ ou S selon le secteur où se trouve le point d’observation (limite entre le secteur 1 et les
secteurs 2-3-4). Les calculs des différents coefficients de réflexion et de détente nécessitent de
connaître la distance en ligne droite entre le centre de la charge et le point de sortie (S’ ou S).
Les différentes relations géométriques, présentées en annexe 5, sont alors utilisées. Ces
relations permettent notamment de calculer les angles de déviation ( S = f(RR ou RM)), les
angles d’incidence et les différentes distances parcourues. Les coefficients de réflexion (CrSol)
et de détente (Cdw) sont déterminés à partir des abaques (abaque des coefficients de détente et
TM5-1300 (2008)). Cette solution est retenue pour chacun des côtés de la barrière de
180
Avec CrS le coefficient de réflexion sur la face avant, Cdw les coefficients de détente sur
chaque arête de la barrière, Crsol le coefficient de réflexion sur le sol en aval
Le secteur 3 est particulier, en effet il existe une détente supplémentaire possible si le point est
situé dans « l’ombre de la face arrière ». Ce secteur est défini comme le secteur 3-bis,
figure 161.
z
Point S z
Point V
H – h’
h Point T
x
' sol
1 1
Point V l'
xd H e
tan (α1 ) tan (α2 )
x Secteur 3-bis :
Données d’entrée :
caractéristiques géométriques (d, e, H, α1, α2 et L)
nature (gazeuse ou chimique condensée) et masse de la charge (W)
Identification du secteur :
secteur « champ libre », secteur 1, secteur 2, secteur 3, secteur 3-bis, secteur 4 ou secteur 4-bis
si RM, calcul du coefficient de réflexion (Crbs) si RM, calcul du coefficient de réflexion (Crs)
suivant le régime de réflexion (θ = f(RR ou RM)), suivant le régime de réflexion (θ = f(RR ou RM)),
calcul du coefficient de détente (Cdw1)
Attention, la formation d’un pied de Mach sur la
calcul du coefficient de détente (Cdb1)
ondes de contournement)
Estimation des caractéristiques de l’onde de souffle : Estimation des caractéristiques de l’onde de souffle :
temps d’arrivé (t0b) temps d’arrivé (t0w)
durée de phase positive (dt+b) durée de phase positive (dt+w)
maximum de surpression (ΔP+b) maximum de surpression (ΔP+w)
coefficient de décroissance de pression de la coefficient de décroissance de la pression de la
relation de Friedlander modifié (C = f(db)) relation de Friedlander modifié (C = f(dw))
Il est donc possible d’estimer la valeur convergée du maximum de surpression à partir d’un
modèle simplifié au niveau du pic de pression en fonction de résultats de simulation
numérique obtenus pour deux calculs avec des tailles de mailles « raisonnables ». Ce modèle
suppose une décroissance linéaire de la pression sur le sommet du pic et non sur l’ensemble
du profil de pression.
δ’écart entre les calculs quasi-convergés est obtenus à partir de la relation (157) en fonction
du facteur de dimensionnement entre les deux tailles d’arête de mailles (Cm). Ce facteur Cm
est égal à γ dans l’exemple sur la figure 163.
u n 1 Pn 1 Pn Pn Pn 1 un
1 1 (157)
Cm Cm
A partir de la relation (157), il est possible d’établir une suite géométrique pour obtenir la
valeur théorique convergée à partir du facteur Cm et de l’écart entre deux résultats de
simulation numérique, relation (158).
1
1
n
u k P0
n 1
Cm
Pn 1 u1 P0 (158)
k 1
1
1
Cm
Si n tend vers l’infini, la valeur convergée peut être estimée de la formule (159).
Pn 1 C m P1 P 0
1
Cm 1
(159)
Avec Cm le facteur entre les dimensions d’arête de mailles, P0 et P1 les maximums de pression
obtenus par simulation numérique
Cette extrapolation linéaire est réalisée à partir de résultats issus de deux simulations
numériques. Ces simulations nécessitent la mise en œuvre de deux tailles de mailles
différentes. La première simulation est réalisée avec une taille de mailles de 1,35 mm au
passage du front de choc, soit des mailles de 5,6 mm/MJ1/3. Ce premier maillage est noté
« grossier [c] », tableau 32. La seconde taille de mailles est trois fois plus petite que la
première (Cm = γ). δa dimension d’arête des mailles est donc fixée à 0,45 mm au passage du
front de choc, soit des mailles de 1,9 mm/MJ1/3, notée « maillage grossier [d] ».
Taille de maille au passage
Approches Energie de la charge de
Type de maillage du front de choc
géométriques gaz (en MJ)
(en mm / mm/MJ1/3)
Géométrie 3D
Taille intermédiaire 1,35 / 5,6
« grossier [c] »
entre « fin » et
Géométrie 3D
«grossier [a] » 0,45 / 1,9
« grossier [d] » 13,75.10-3 (charge (R1))
Rappel : Géométrie
2D-axisymétrique « fin » 0,20 / 0,80
« fin »
Tableau 32 : Conditions de calculs des « merlons courts »
y 0,085 m
Capteur à 0,60 m
x
δes données expérimentales obtenues pour l’onde incidente sont comparées aux résultats de
simulation numérique. δa mise en œuvre des maillages « grossiers » induit un défaut de
convergence numérique qui se traduit par une sous-estimation du maximum de surpression
d’environ 9 à 10% pour ce point de mesure. δ’extrapolation linéaire réalisée à partir des
résultats de simulation sous-estime de 5%, cet écart se situe alors largement dans l’intervalle
d’incertitude expérimentale de +/- 10% des mesures.
185
d = 0,085 m
Charge H = 0,19 m
(0,0)
α1 = 45° α2 = 45° x
Coupe au sol
L = 0,40 m
Charge
(0,0) x
Les capteurs sont placés autour de la barrière de protection afin d’analyser l’impact des
différentes recombinaisons sur l’effet protecteur du merlon (7 points de mesure, notés de S1 à
S7), figure 166. La visualisation des iso-contours de pression autour de la barrière permet de
rappeler la complexité des phénomènes de recombinaison en aval de la barrière.
Onde réfléchie
Onde issue d’une recombinaison
S1
Capteurs
S5 S2
S6 S4
t = 1,7 ms
S7 S3
Onde issue du merlon x
Figure 166 : Positions des différents capteurs placés autour de la barrière de protection du merlon 1A « court »
Le tableau 33 présente la position des capteurs situés au niveau du sol (y = 0,0 m).
Capteurs S1 S2 S3 S4 S5 S6 S7
x (m) 0,30 0,40 0,60 0,60 0,60 0,60 0,70
z (m) 0,30 0,30 0,30 0,20 0,10 0,00 0,00
Tableau 33 : Positions des capteurs autour du merlon court de la série 1
186
Maximum de
surpression (bar)
0,300
Onde issue du merlon
0,168
0,037
-0,095
x
-0,227 Centre de la charge
(0,0,0)
t = 0,85 ms
Onde réfléchie
Coupe dans l’axe du
a) Capteur S1 capteur S1
0,170
0,060
-0,051
-0,161
Onde réfléchie
t = 1,2 ms
Figure 167 : Visualisations des champs et des profils de surpression pour les capteurs S1 et S2 du
merlon 1A « court »
Pour les capteurs S1 et S2, les profils de pression obtenus par la méthode d’estimation rapide
ne correspondent pas aux données expérimentales et aux résultats de simulation numérique.
La méthode d’estimation rapide estime pour S2 la présence de deux pics de pression. Cette
différence peut être due à l’une des hypothèses de l’approche. En effet, le temps d’arrivée (t0)
de l’onde de choc est estimé à partir de la formule (140) et de la distance parcourue entre le
centre de la charge et le point d’observation. Cette formule est établie à partir de données
expérimentales pour une configuration en champ libre (sans obstacle). δe temps d’arrivée
estimé ne prend donc pas en compte l’impact des différents phénomènes physiques (réflexion
187
et détente) sur cette caractéristique de l’onde de choc. La prise en compte de cette variation
est un des axes d’amélioration de la méthode rapide. Cependant, l’écart sur le maximum de
surpression entre les estimations rapides et les résultats de simulation numérique est
satisfaisant compte-tenu du coût de calcul négligeable de cette approche (intervalle de
+/- 20%).
Onde réfléchie
Maximum de
Onde issue du merlon surpression (bar)
0,220
0,125
0,029
a) Capteur S3
Maximum de
Onde issue du merlon
surpression (bar)
0,150
Onde réfléchie
0,098
0,045
-0,008
-0,060
t = 1,9 ms
Coupe dans l’axe du
capteur S4 x
Première recombinaison
Seconde recombinaison sur l’axe
b) Capteur S4
Figure 168 : Visualisations des champs et des profils de surpression pour les capteurs S3 et S4 du
merlon 1A « court »
188
0,098
0,045
-0,008
-0,060
Coupe dans
l’axe du x t = 1,9 ms
capteur S5
Première recombinaison
Seconde recombinaison sur l’axe
Capteur S5
Figure 169 : Visualisation des champs et des profils de surpression pour le capteur S 5 du merlon 1A « court »
Les figures 168 et 169 présentent les profils de surpression obtenus pour les points de mesure
S3, S4 et S5. La méthode d’estimation rapide inverse l’ordre d’arrivée des ondes de choc pour
le capteur S3. En effet, la méthode rapide estime que le premier pic de surpression correspond
à l’onde de contournement. Cependant, les données expérimentales et les résultats de
simulation permettent d’observer que le premier pic de pression correspond au passage de
l’onde issue du merlon. Cette inversion peut être imputée aux calculs du temps d’arrivée (t0).
Les valeurs du maximum de surpression obtenues par l’approche rapide restent pour les deux
ondes dans un intervalle d’incertitude acceptable de +/- 20%).
Le capteur S4 permet d’observer les différentes recombinaisons d’ondes de choc en
aval du merlon 1A. La position du capteur S4 est particulière. Le capteur est situé à la limite
entre le secteur 4 et le secteur 4-bis. Ces deux secteurs se différencient par la prise en compte
de la seconde détente pour l’onde de contournement. Dans le secteur ζ, l’onde de
contournement n’est soumise qu’à une seule détente (Cdb1). δ’onde de contournement qui
atteint le secteur 4-bis subit une seconde détente (Cdb2). Ce découpage en nombre de secteurs
fixés montre ses limites et peut expliquer certaines petites variations brutales de pression pour
les points de mesure placés sur la limite d’un secteur. Pour le capteur S4, un écart est visible
sur l’onde de contournement (1er pic de pression, - 50% par rapport aux données
expérimentales). Cependant, l’estimation du maximum de surpression reste dans un intervalle
d’incertitude de +/- 20%.
δ’onde issue du merlon se réfléchit en aval du merlon et se focalise avec l’onde de
contournement de chaque côté de la barrière de protection. La somme des pressions au niveau
du point d’observation permet de tenir compte de la présence de ce phénomène physique. δes
ondes issues de cette première recombinaison, de chaque côté du merlon, se focalisent à
nouveau dans l’axe du merlon. Cette seconde recombinaison implique d’ajouter, dans la zone
de recombinaisons, la contribution des deux ondes de contournement. Dans une première
approche, il est possible de considérer le simple doublement du maximum de surpression pour
l’onde de contournement (deux ondes de contournement) dans la zone des recombinaisons
(db > dw) située en aval du merlon (x > distance entre le centre de la charge et le pied aval du
merlon). C’est une approximation qui est faite. δ’augmentation du maximum de surpression
189
de l’onde issue de cette seconde recombinaison, sur l’axe en aval du merlon, n’est pas
linéaire, elle dépend notamment du niveau de pression incident des ondes de contournement.
La reconstruction du profil de surpression pour le capteur S5 prend en compte cette
approche de la seconde recombinaison. δ’écart entre l’estimation rapide et l’extrapolation
linéaire (valeur convergée des résultats de simulation) est de - 27%. Une amélioration
possible de l’approche rapide serait la modélisation de cette seconde recombinaison.
0,088
0,045
-0,003
-0,040
0,083
0,045
0,008
-0,030
t = 2,3 ms
Coupe dans l’axe du
capteur S7
Première recombinaison
b) Capteur S7
Figure 170 : Visualisations des champs et des profils de surpression pour les capteurs S6 et S7 du
merlon 1A « court »
Les profils de surpression pour les capteurs S6 et S7 sont présentés sur la figure 170. δ’onde
issue de la seconde recombinaison rattrape le premier front de choc (onde issue du merlon) et
les ondes issues de la première recombinaison. En champ lointain, ces ondes se recombinent
entre-elles. Ce phénomène est visible sur les profils de surpression issus des mesures ainsi que
190
sur les résultats de simulation numérique. Dans cette zone de recombinaisons, la méthode
d’estimation ne tient pas compte de la variation de la vitesse du son en aval du premier choc.
En champ lointain, les résultats issus de l’estimation rapide peuvent donc conduire à une
surestimation du maximum de surpression. La prise en compte de ce phénomène est
également un axe d’amélioration possible de l’approche.
Le tableau 34 synthétise les résultats obtenus pour le maximum de surpression à partir
des trois approches utilisées lors de cette étude : expérimentale, numérique et estimation
rapide.
Capteurs Maximum de surpression (bar) Ecart entre
Les résultats de Estimation rapide et les
Données Simulation Estimation
simulation et les données données
expérimentales numérique rapide
expérimentales expérimentales
S1 0,309 0,278 0,213 - 10,0% - 31,1%
S2 0,248 0,209 0,197 - 15,7% - 20,5%
S3 0,160 0,137 0,122 - 14,4% - 23,8%
S4 0,155 0,144 0,120 - 7,1% - 22,6%
S5 0,130 0,134 0,097 + 0,3% - 25,4%
S6 0,119 0,104 0,091 - 12,6% - 23,5%
S7 0,106 0,093 0,081 - 12,3% - 23,6%
Tableau 34 : Comparaison entre les données expérimentales, les résultats de simulation numérique et les
estimations pour le merlon 1A « court »
Les estimations obtenues pour les différents points de mesure permettent d’établir pour la
méthode rapide un intervalle de confiance de +/- 30% par rapport aux données
expérimentales.
H = 0,19 m
Charge
α1 = 90° α2 = 45° x
Coupe au sol
L = 0,40 m
Charge x
Pour cette configuration de barrière de protection, plusieurs points de mesure sont placés
autour du merlon. Leurs positions sont présentées sur la figure 172 et données dans
le tableau 35.
S5
S6
x S4 S1
Capteurs S3
S2
t = 1,7 ms Onde issue du merlon
Figure 172 : Positions des différents capteurs placés autour de la barrière de protection du merlon 2B « court »
Capteurs S1 S2 S3 S4 S5 S6
x (m) 0,40 0,60 0,60 0,60 0,60 0,70
z (m) 0,30 0,30 0,20 0,10 0,00 0,00
Tableau 35 : Positions des capteurs autour du merlon court de la série 2
Maximum de
surpression (bar)
0,220
0,160
0,100
x
0,040
-0,020
t = 1,10 ms
z
Coupe dans l’axe du
Onde de contournement capteur S1
a) Capteur S1
Onde réfléchie
0,108
0,065
0,023
Coupe dans
-0,020 l’axe du
capteur S2
t = 1,75 ms
Onde de contournement
b) Capteur S2
Figure 173 : Visualisations des champs et des profils de surpression pour les capteurs S1 et S2 du
merlon 2B « court »
0,108
0,065
t = 1,70 ms
Onde de contournement
a) Capteur S3.
Maximum de
surpression (bar)
0,140
0,100
0,060
0,020
Coupe dans
l’axe du
-0,020
capteur S4
t = 1,75 ms
Onde réfléchie
b) Capteur S4
Figure 174 : Visualisations des champs et des profils de surpression pour les capteurs S3 et S4 du
merlon 2B « court »
premières recombinaisons sur les côtés du merlon (visualisation du capteur S6, figure 175-b).
Cette première version de la méthode rapide ne tient pas compte de ce phénomène physique.
Onde de contournement
Maximum de
surpression (bar)
0,120
Onde issue du merlon
0,053
0,019
-0,015
t = 1,75 ms
Coupe dans
l’axe du
capteur S5
Première recombinaison
a) Capteur S5
Onde de contournement
0,090
0,055
0,020
Coupe
-0,015 dans l’axe
du Seconde recombinaison
t = 1,85 ms capteur S6
Première recombinaison
b) Capteur S6
Figure 175 : Visualisations des champs et des profils de surpression pour les capteurs S5 et S6 du
merlon 2B « court »
Le tableau 36 synthétise les résultats obtenus pour le maximum de surpression selon les trois
méthodes d’évaluation des effets mécaniques de cette étude du merlon « court » 2B.
195
Le niveau de confiance des résultats issus de la méthode d’estimation rapide pour le merlon
2B est identique à celle du merlon 1A : sous-estimation jusqu’à 30% par rapport aux données
expérimentales.
0,085 m
0,085 m
1,10 m
1,10 m
Charge
Charge
Merlon Merlon
x 0,465 m x 0,465 m
1,10 m 1,10 m
Z1 Z2 Z3 Z4 Z5
Figure 176 : Visualisations des zones d’effets pour le merlon βB à partir des résultats de simulation numérique et
de la méthode d’estimation rapide, vue aérienne plane au-dessus du sol
Une grande similarité est observée entre les deux approches pour les zones les plus
dangereuses (Z1, Z2 et Z3). Ces visualisations permettent également de constater une
augmentation du maximum de surpression dans le secteur 3-bis (« zone dans l’ombre de la
face arrière »). En champ lointain de la barrière de protection, il avait déjà été constaté que
l’approche à tendance à surestimer le maximum de surpression. Cette surestimation retarde
l’apparition de la zone d’effets Z5. Cette particularité n’est cependant pas pénalisante dans la
mesure où la problématique est celle de la sécurité pyrotechnique. Il est donc toujours
préférable de majorer les zones d’effets plutôt que de les minorer.
D’autres travaux seraient nécessaires afin d’améliorer la version proposée dans ce
chapitre. δ’amélioration de la méthode de superposition en tenant compte d’une variation de
vitesse du son après le premier choc permettrait d’affiner le calcul du temps d’arrivée du
second pic de pression. δ’introduction des effets de compensation afin d’éviter les
surestimations en champ lointain est également un axe d’amélioration possible de la méthode
d’estimation rapide. Enfin, cette méthode doit être testée sur d’autres configurations afin
d’identifier d’autres limitations et d’autres axes d’améliorations possibles.
En résumé, la méthode rapide présente les avantages suivant :
pas de ressources de calcul nécessaires, un simple ordinateur « de bureau » suffit,
12
Le CEA – DAM dispose de l’intégralité des données (375 abaques). Ces abaques sont disponibles sur demande auprès du
CEA – DAM de Bruyères-Le-Châtel (Essonne (91) - France).
198
y d e
Dans le cadre de ces travaux de thèse, la charge est placée au sol (HOB = 0 m). Les
simulations numériques sont réalisées en géométrie 2D-axisymétrique. Cette approche est
équivalente, en géométrie 3D, à une barrière de protection de forme torique dont la courbure
est fonction de la distance (d) (§ β.β). δa longueur de la barrière de protection n’est donc pas
étudiée lors de l’étude paramétrique (L = ∞). δ’impact de cette dimension sur l’effet
protecteur d’un merlon peut être évalué à partir de simulations numériques en géométrie 3D.
δes intervalles de l’étude paramétrique sont donnés dans le tableau 37. Les intervalles
des expériences à petites échelles sont également représentés (coefficient d’équivalence 2,35
(Trélat (2006)) permet une approximation de la masse équivalente en kg de TNT). Les
barrières de protection étudiées à petites échelles (échelle réduite, explosion gaz) se situent
dans les intervalles de l’étude paramétrique pour plusieurs paramètres (H, e, α1 et α2) à
l’exception de la distance d. En effet, les simulations ne sont pas réalisées en champ proche de
la charge afin de se placer dans la gamme où l’approche en géométrie βD axisymétrique est
licite (l’effet axisymétrique), mais aussi pour éviter une éventuelle déformation de l’onde
réfléchie par la présence des produits de détonation (§ 2.4.2 et § 3.2.3.2). Les charges
présentées sont équivalentes à une charge sphérique en champ libre (soit pour une charge de
1 kg de TNT, une charge hémisphérique de 0,5 kg de TNT placée au niveau du sol).
199
H e d
Légendes 1/3 1/3
α1 et α2
(en m/kg ) (en m/kg ) (en m/kg1/3)
Intervalles de 5,0 5,0 90° 5,4
l’étude 3,6 3,6
paramétrique
45°
Intervalles des
1,5
expériences à 1,4
0,3 1,3
petites échelles 0,4 0,2 20°
0,7
ne sont pas raffinées. Cette zone est en effet en dehors du périmètre de l’étude. δe maximum
de surpression est extrait depuis le pied aval du merlon (distance fixée suivant la géométrie
étudiée) jusqu’à η8 m (soit ∆P+ = 0,014 bar, niveau de surpression inférieur au seuil de la
dernière zone d’effets (Z5 ; 0,020 bar)). La longueur du domaine calculé est de 70 m.
y Lmax = 70 m
Aucun raffinement
Lmax = 70 m
L
40 m
1,602 m 2m
L x
0,002 m
(0, 0) Zone d’extraction, pas de regroupement
Figure 178 μ Représentation schématique d’un domaine de calcul de l’étude paramétrique
Le standard de calcul utilisé pour une charge de 1 kg de TNT est synthétisé dans le
tableau 39.
Taille de mailles au Temps de transition Nouvelle taille de mailles
Approche géométrique
début du calcul du maillage jusqu’à la fin du calcul
Géométrie 2D-axisymétrique
2 mm 1,4 ms 4 mm
« fin »
Tableau 39 : Standard de calcul pour une charge de 1 kg de TNT
Le choix des conditions de calcul dans code de calcul HERA pour une charge de 1 kg
de TNT est établi à partir d’une configuration en champ libre (sans obstacle). Les résultats de
simulation numérique pour le maximum de surpression et le temps d’arrivée sont comparés à
deux références issues de la littérature (Kinney (1985) et de l’abaque du TM5-1300 (2008)),
figure 179.
201
Figure 179 : Evolution du maximum de surpression et du temps d’arrivée de l’onde de choc en champ libre pour une charge de TNT au niveau du sol
202
δ’écart entre les deux références bibliographiques est de - 8,3% à 5 m/kg1/3 et de + 17,6% à
40 m/kg1/3 par rapport à Kinney (1985), tableau 40. Les résultats obtenus par simulation
numérique se situent dans cet intervalle d’incertitude. La variation la plus importante est
observée en champ lointain, pour un maximum de surpression inférieur au seuil de la dernière
zone d’effets (Z5 ; ∆P+ < 0,020 bar). Les conditions physico-numériques de calcul permettent
de se situer dans un niveau de convergence numérique satisfaisant.
Ecart entre Ecart entre les Ecart entre les
Distance Kinney TM5-1300 Simulation
Kinney (1985) résultats de résultats de
réduite (1985) (2008) numérique
et le TM5- simulation et simulation et le
(z, en mkg1/3) (bar) (bar) (bar)
1300 (2008) Kinney (1985) TM5-1300 (2008)
1,45 4,29 3,84 3,84 + 10% - 10% 0%
5,00 0,288 0,302 0,305 - 4,6% + 6,0% + 1,0%
10,00 0,099 0,108 0,107 - 8,3% + 8,0% - 1,0%
25,00 0,034 0,033 0,032 + 3,0% - 6,0% - 3,0%
40,00 0,020 0,017 0,018 + 17,6% - 10% + 6,0%
Tableau 40 : Comparaison des résultats du maximum de surpression entre la simulation numérique et les
différentes références issues de la littérature pour une configuration en champ libre
Figure 180 μ Evolution de la durée de phase positive et de l’impulsion positive en champ libre pour une charge de TNT au niveau du sol
204
5.2 Abaques
Dans le cadre de l’étude paramétrique, les simulations numériques sont réalisées en
géométrie 2D-axisymétrique « fin » pour une charge de 0,5 kg de TNT au sol, soit une charge
de 1 kg de TNT en équivalent champ libre. Les relevés du maximum de surpression (au sol et
à 1,60 m/kg1/3) permettent de construire la base de données. A partir de cette base de données,
375 abaques ont pu être réalisés. Un exemple d’abaque est présenté sur la figure 181.
d,
Paramètres fixés :
H,
α2.
e,
Variables :
α1.
Les abaques sont construits pour trois paramètres fixés, soit 125 abaques : la distance
d, la hauteur H et l’angle d’inclinaison α2. δa variation des paramètres e et α1 est représentée
sur chaque abaque par un réseau de βη courbes. δ’abscisse est en distance curviligne depuis le
205
1
Entrées : W1, d1, H1, e1, α11 et α12
d = 1,4 m 1
e = 0,2 m Sortie : maximum de surpression à 2 m
en aval du merlon
H1 = 1,55 m
1
W = 0,5 kg α12 = 40° 2m
α11 = 30° Point d’observation
Figure 182 μ Exemple d’utilisation des abaques pour une configuration existante
Pour l’exemple 1 présenté sur la figure 182, l’utilisateur souhaite connaître le maximum de
surpression et la zone d’effets associée pour un point d’observation placé à 2 m en aval du
merlon. Dans cette configuration et compte tenu que les abaques décrivent de manière
discrète la variation de 5 paramètres (d, H, e, α1 et α2), l’utilisateur a besoin de β abaques qu’il
faudra interpoler (d = 1,4 m/kg1/3, soit d1 = 1,4 m, H = 1,55 m/kg1/3, soit H1 = 1,55 m, α2 = 30°
et α2 = 45°, pour α12 = 40°), figure 183 et 184.
13
Ce 1er exemple est choisi volontairement très simple sur la valeur de la charge afin de permettre une lecture directe des
abaques sans conversion en mTNT1/3
206
2 m/kg1/3
Figure 183 μ Exemples d’abaque pour une barrière de protection existante (d = 1,4 m/kg1/3, H = 1,55 m/kg1/3, α2 = 30°)
207
2 m/kg1/3
Figure 184 μ Exemples d’abaque pour une barrière de protection existante (d = 1,4 m/kg1/3, H = 1,55 m/kg1/3, α2 = 45°)
208
Les courbes rouges sur les abaques encadrent la configuration d’intérêt (e = 0,β0 m/kg1/3, soit
e1 = 0,20 m et α11 = 30°). Une interpolation linéaire est réalisée entre ces deux courbes pour
un angle d’inclinaison α2 de 40°, figure 185.
Figure 185 : Evolution du maximum de surpression en aval d’une barrière de protection existante (exemple 1)
δ’interpolation linéaire entre les deux courbes permet de calculer le maximum de surpression.
Pour un point d’observation positionné à β m en aval du merlon, le maximum de surpression
est évalué à 0,107 bar, soit la zone d’effets Z4. La présence de la barrière de protection permet
d’atténuer la surpression de γγ% par rapport à une configuration champ libre (sans obstacle).
La version proposée de la méthode rapide permet également d’estimer le maximum de
surpression (∆P+ = 0,07β bar, soit également la zone d’effets Z4). Cette estimation sous-estime
néanmoins le maximum de surpression de 32% (utilisation de la formule empirique de Kinney
(1985) pour le calcul du maximum de surpression incident, formule (6)). Cet écart est proche
de l’intervalle de confiance évalué à +/- 30% pour la méthode d’estimation rapide.
Si la masse de la charge est différente de 0,5 kg de TNT, le calcul des distances réduites sera
différent mais la démarche d’utilisation des abaques sera pour sa part identique (par exemple
si W1 = 13,5 kg de TNT, alors d1 = 1,4 m, soit d = 0,47 m/kg1/3, ramenée à une charge
sphérique).
Entrées :
contrainte d’espace du merlon au
Longueur au sol disponible,
long2 = 6,5 m
sol fixé à 6,5 m, soit 4,5 m/kg1/3,
angle d’inclinaison α2 fixé à 45°,
Zone de stockage
d2 = 2,02 m
H2 = 2,24 m
hauteur (H) limitée à 2,24 m, soit
1,5 m/kg1/3.
Voie Sortie μ géométrie optimale (α1 et e)
d’accès à partir des contraintes.
W2 = 1,5 kg α22 = 45°
Dans l’exemple, une source explosive est située à 11,73 m d’une cuve de stockage de
produits inflammables, soit R2 = 8,13 m/kg1/3. Afin de protéger la cuve des effets mécaniques
d’une éventuelle explosion, l’utilisateur souhaite placer une barrière de protection. La
première contrainte est l’espace disponible au sol pour construire la barrière de protection. La
barrière est placée entre une voie d’accès et une zone de stockage, soit une longueur au sol du
merlon limitée à 6,5 m (long2 = 4,5 m/kg1/3).
Les « bonnes pratiques » sur le dimensionnement d’une barrière de protection recommandent
des angles d’inclinaison des faces avant et arrière de λ0°, afin d’accentuer l’atténuation par la
présence de deux fortes détentes situées en amont et en aval du sommet de la structure. Ce
dimensionnement implique alors de fortes contraintes mécaniques sur la structure : surface
verticale en vue directe de la charge → ondes réfléchies maximales en pression. L’utilisateur
souhaite placer un remblai en terre sur la face aval du merlon (α22 inclinée à 45°) afin d’éviter
la présence de nouveaux fragments issus de la face arrière du merlon. Pour des raisons
financières, la hauteur (H2) est fixée à 2,24 m, soit 1,55 m/kg1/3.
δa contrainte de sécurité impose d’obtenir le plus faible maximum de surpression possible en
aval du merlon. Les contraintes financières et d’espaces sont ici les variables d’entrée des
abaques (H2 = 2,24 m et long2 = 6,5 m). δ’utilisateur sélectionne les abaques qui
correspondent à la configuration d’intérêt avec une abscisse dite en « distance absolue »
(distance entre le centre de la charge et le point d’observation), figure 187. Les abaques
sélectionnés permettront d’identifier la géométrie optimale du merlon pour un maximum de
surpression le plus faible possible au niveau de la cuve et de déterminer de surcroît la zone
d’effets associée.
210
Configurations possibles
Configuration optimale
Géométrie optimale
Pied aval
à 5,9 m/kg1/3 Zone étudiée à 8,13 m/kg1/3
Les contraintes sur le dimensionnement impliquent que le pied aval du merlon soit placé à
8,52 m (soit 5,9 m/kg1/3). Les configurations possibles sont indiquées par un disque vert sur
l’abaque. δe point d’observation est placé à 11,73 m, soit 8,13 m/kg1/3 (position de la cuve de
stockage de produits inflammables). Le dimensionnement étant libre devant le merlon, la
barrière sélectionnée correspond au maximum de surpression le plus faible du réseau de
courbes des configurations dites « possibles ». La géométrie optimale en fonction des
contraintes imposées est donc un merlon placé à 2,02 m (d2, soit 1,4 m/kg1/3) de la charge
explosive, avec une épaisseur au sommet de 3,75 m (e2, soit 2,6 m/kg1/3), une hauteur de
2,24 m (H2, soit 1,55 m/kg1/3), une face avant inclinée à 90° (α21) et une face arrière inclinée à
45° (α22), figure 188.
2,02 m
Voie
d’accès α21 = 90° α22 = 45°
W2 = 1,5 kg
Figure 188 : Géométrie optimale obtenue à partir des abaques en fonction des contraintes
Notons que ce dimensionnement est en accord avec les « bonnes pratiques » présentées
précédemment (angle d’inclinaison de λ0° et une épaisseur la plus importante possible). Le
maximum de surpression incident au sol au niveau de la cuve (à 11,73 m, soit 8,13 m/kg1/3)
est de 0,0ηζ bar, soit la zone d’effets Z4. La méthode rapide estime un maximum de
surpression de 0,059 bar, soit une valeur assez proche. Cette géométrie implique néanmoins
des efforts mécaniques importants sur la face avant du merlon. Le maximum de surpression
réfléchie est en effet estimé à 21,1 bar au début de la face avant du merlon d’après l’abaque de
réflexion du TM5-1300 (2008).
La longueur au sol du cette géométrie de merlon est de 5,99 m (soit 4,15 m/kg1/3).
Cette distance est inférieure à la longueur limite de 6,5 m (long2, soit 4,5 m/kg1/3). En accord
avec les recommandations sur le dimensionnement d’une barrière de protection, l’épaisseur
peut être augmentée afin d’augmenter l’atténuation par la distance parcourue sur le sommet,
soit une épaisseur de 4,26 m (e2bis, soit 2,95 m/kg1/3). Pour cette nouvelle configuration, le
maximum de surpression estimé à l’aide des abaques est de 0,0ζη bar, soit une diminution
moyenne de 16% par rapport à la première proposition, figure 189. Avec cette nouvelle
géométrie de merlon en ce point, la zone d’effets passe en Z5. La méthode rapide estime un
maximum de surpression de 0,052 bar, soit une surestimation de 15% par rapport aux
abaques.
2,02 m
Voie
d’accès α21 = 90° α22 = 45°
W2 = 1,5 kg
Le second paramètre qui peut être modifié est la distance d2. La géométrie optimale de merlon
issue des abaques (1er proposition de merlon) est donc placée à 2,53 m (d2bis, soit
1,75 m/kg1/3). δ’interpolation réalisée sur les différents abaques permet d’estimer un
maximum de surpression de 0,049 bar, soit une diminution de 10% par rapport à la première
proposition de géométrie (estimation de 0,043 bar par la méthode rapide), figure 190.
2,53 m
Voie
α21 = 90° α 2 = 45°
2
d’accès
W2 = 1,5 kg
Si cette barrière de protection ne répond pas aux attentes, il faut revoir les contraintes
(augmentation de la hauteur par exemple) et réitérer l’opération autant de fois que nécessaires
grâce aux abaques établis.
utilisateur pourrait être réalisé, cet outil regrouperait les abaques 14 et l’interpolation des
données ayant permis de les construire ainsi que la méthode rapide. Le CEA – DAM travaille
actuellement à sa réalisation.
Enfin les efforts mécaniques subis par le merlon lors de l’interaction de l’onde de souffle
peuvent également être calculés avec une approche similaire à la méthode d’estimation rapide.
14
Le CEA – DAM dispose de l’intégralité des données (375 abaques). Ces abaques sont disponibles sur demande auprès du
CEA – DAM de Bruyères-Le-Châtel (Essonne (91) - France).
214
215
CONCLUSION
Les travaux présentés dans ce mémoire de thèse concernent la propagation d’une onde
de souffle aérienne générée par la détonation d’une charge gazeuse (propane-oxygène à la
stœchiométrie) et d’une charge chimique condensée (TNT, Trinitrotoluène), en présence
d’une barrière de protection. Ces deux familles d’explosif se différencient par une
composition chimique et une célérité de détonation très différentes. Ces deux natures de
source conduisent donc à l’établissement d’ondes de souffle aériennes qualitativement
semblables et quantitativement différentes. De plus, à énergie équivalente, ces deux natures de
charge explosive dégagent des volumes de produits de détonation différents, et ont des
vitesses d’interface air/produits de détonation (vitesse fonction de la célérité de détonation)
également très différentes.
La présence des produits de détonation en champ proche de la structure et de l’onde de choc
peut donc conduire à une déformation du phénomène de réflexion sur la paroi. Ainsi en
fonction de la nature de la charge explosive, de la position et des dimensions de la structure,
le phénomène de réflexion peut aboutir à la formation d’un double pied de Mach. Ce second
pied de Mach provient de l’interaction de la première onde réfléchie de la paroi sur l’interface
entre l’air et les produits de détonation. δ’apparition de ce nouveau pied dépend notamment
de la nature de la charge explosive.
Dans le cadre de ces travaux de thèse, la détonation de la charge explosive est explicitement
simulée avec l’équation d’état propre à l’explosif étudié (formulation JWδ), en prenant
spécifiquement en compte la nature de la charge : gazeuse ou chimique. δes équations d’état
sont obtenues à partir d’un code de thermochimie du CEA – DAM.
Ce mémoire de thèse présente l’analyse de plusieurs configurations expérimentales
destinées à étudier la propagation d’onde de souffle aérienne issue de la détonation de charge
gazeuse (propane-oxygène à la stœchiométrie) en présence d’obstacle de formes
parallélépipédiques, cylindriques et trapézoïdales (merlons). Les données expérimentales
acquises pour ces configurations servent d’une part à la compréhension des phénomènes
physiques rencontrés et d’autre part à la validation de l’outil de simulation constitué du
solveur HERA et la plateforme cible TERA 100.
Les simulations réalisées s’effectuent selon deux approches géométriques μ l’une en géométrie
2D-axisymétrique et l’autre en géométrie γD. δa simulation d’une géométrie expérimentée
tridimensionnelle à l’aide d’une simulation en géométrie βD-axisymétrique peut conduire, en
fonction des dimensions du merlon (longueur), à un effet de compression (effet
axisymétrique) propre à cette approche 2D (le merlon devient torique par construction
axisymétrique autour de l’axe de révolution). La phénoménologie issue de la simulation peut
alors ne pas respecter la réalité de l’expérience. D’un autre côté, l’utilisation d’une approche
en géométrie 3D impose une dégradation du niveau de discrétisation (maillage « grossier »)
afin d’assurer la compatibilité du cas de simulation avec les moyens de calcul à disposition.
δ’utilisation d’un maillage « grossier » induit un niveau de convergence numérique
insuffisant qui se traduit par une sous-estimation des résultats sur la valeur crête du maximum
de surpression.
Néanmoins dans le cadre de cette étude, il a été montré qu’une approche en géométrie
2D-axisymétrique était valable pour simuler avec une précision satisfaisante (utilisation d’un
maillage « fin » pour atteindre un niveau de convergence numérique suffisant) le cas d’une
barrière en géométrie 3D (impact négligeable des ondes de contournement). Le choix de
simuler en géométrie 2D-axisymétrique permet de réduire le nombre de mailles et le temps de
calcul par rapport à une étude en géométrie 3D pour des conditions numériques par ailleurs
égales tout en s’assurant des niveaux de discrétisation suffisamment fins pour atteindre un
216
bon niveau de convergence numérique. Il s’agit donc de trouver le meilleur compromis entre
les deux approches en fonction de la géométrie étudiée. De ce fait, il a été montré que la
réalisation des abaques, par une étude paramétrique, pouvait s’effectuer avec l’approche en
géométrie 2D-axisymétrique.
La comparaison des données expérimentales et des résultats de simulation numérique
permet de valider les capacités restitutives et donc prédictives de l’outil de simulation pour la
propagation d’une onde de souffle aérienne issue de la détonation d’une charge explosive
(code de calcul HERA, conditions physico-numériques, équation d’état des produits de
détonation,…). δa restitution des données expérimentales est correcte avec la mise en œuvre
d’un maillage « fin » qui permet de se placer à un niveau de convergence numérique
satisfaisant (écart de l’ordre des intervalles d’incertitude expérimentale de +/- 10% à +/- 14%
suivant la configuration). Les résultats de simulation numérique obtenus en géométrie 3D,
pour les configurations qui rendent indispensables cette approche, sous-estiment globalement
le maximum de surpression (effet du maillage « grossier »). La restitution des expériences
concernées n’est donc pas parfaite, mais elle est toujours qualitativement satisfaisante et
quantitativement correcte avec un écart d’environ β0%.
Les résultats de simulation numérique et les données expérimentales acquises pour les
différentes géométries de barrière de protection étudiées ont permis d’analyser l’influence des
différents phénomènes physiques sur l’apparition d’un effet protecteur en aval de la barrière
de protection en fonction de ses dimensions. A partir de cette analyse quelques
recommandations sur le dimensionnement d’une barrière de protection ont pu être établies.
Ces recommandations, ou « bonnes pratiques », complètent ainsi celles existantes : taille
minimale de l’épaisseur au sommet du merlon, influence de l’inclinaison des parois, impact
de la distance entre le centre de la charge et la barrière en fonction de la position de la zone à
protéger,…
δ’étude du couplage des différents phénomènes physiques présents lors de
l’interaction de l’onde de choc avec une barrière de protection a contribué à la mise en place
d’une méthode rapide de pré-dimensionnement d’une barrière de protection, également
appelée « méthode géométrique d’estimation rapide ». Cet outil estime l’évolution du
maximum de surpression et la décroissance de pression autour d’une barrière de protection en
fonction de ses paramètres géométriques pour une géométrie 3D ou une géométrie 2D. Le
temps de restitution de cette méthode d’estimation rapide ne prend que quelques secondes (où
un code CFD nécessitera de nombreuses heures sur supercalculateur) avec une marge d’erreur
de +/- 30% par rapport aux données expérimentales. Néanmoins, la version proposée dans ce
mémoire nécessite encore quelques améliorations : impacts des phénomènes physiques sur le
temps d’arrivée de l’onde, prise en compte de la variation de la vitesse du son en aval du
premier choc,… Cette méthode est également limitée à des configurations de barrière de
protection avec un faible niveau de pression dans la zone des recombinaisons d’ondes de choc
en aval du merlon. Cependant, cette méthode rapide est un outil d’estimation satisfaisant au
regard du temps de calcul (de l’ordre de quelques secondes).
Cet outil d’estimation rapide et les « bonnes pratiques » issues de l’étude de
différentes géométries de barrière de protection sont complétés par la réalisation d’abaques
d’aide au dimensionnement. Ces abaques sont construits à partir de résultats fournis par l’outil
de simulation en géométrie 2D-axisymétrique « fin » pour la détonation d’une charge
sphérique équivalente à 1 kg de TNT (distance réduite en m/kg1/3, équivalente à une charge
hémisphérique de 0,5 kg de TNT placée au niveau du sol).
δes capacités prédictives de l’outil de simulation ont été préalablement validées par l’étude de
configurations expérimentées pour la détonation d’une charge gazeuse (propane-oxygène à la
217
BIBLIOGRAPHIE
[1] L. Allain, « Barricade influence on blast wave propagation », SNPE / Division Défense
Espace, 1994 ;
[2] ARIA, descriptions d’accidents disponible sur le site de l’ARIA du ministère chargé du
développement durable, www.aria.developpement-durable.gouv.fr/ ;
[3] Arrêté du 20/04/2007, « fixant les règles relatives à l'évaluation des risques et à la
prévention des accidents dans les établissements pyrotechniques » ;
[4] P. Ballereau, P. Carribault, F. Duboc, D. Dureau, C. Enaux, H. Jourdren, M. Pérache,
« Méthode de raffinement adaptatif de maillage et modèles avancés de programmation
pour le calcul haute performance », Chocs, 41, p. 81-87, 2012 ;
[5] W. Baker, « Explosions in air », University of Texas Press, Austin, 1973 ;
[6] G. Ben-Dor, « Shock wave reflection phenomena », Springer 2nd edition, 1991 ;
[7] J. Borgers, « Blast walls reviewed », 21th Military Aspects of Blast and Shock, October
4 – 8th 2010 Jerusalem (Israel), 2010 ;
[8] H.L. Brode, « Annual Review of Nuclear Science », vol 18, 153-202, 1968 ;
[9] Circulaire DPPR/SEI2/IH-07-0111 du 20/04/2007, « relative à l’application de l’arrêté
fixant les règles relatives à l’évaluation des risques et à la prévention des accidents dans
les établissements pyrotechniques » ;
[10] F. Duboc, H. Jourdren, C. Matignon, « Approche hydrodynamique en détonique :
méthodes de simulation avancées », Chocs, 38, p. 64-71, 2010 ;
[11] S. Eveillard, « Etude de l’influence d’une barrière de protection de type merlon face à
une explosion », étude préliminaire du projet ANR BARPPRO, CEA – DAM, 2010 ;
[12] S. Eveillard, « Rapport Bibliographique du CEA – DAM, tâche 1 Projet BARPPRO »,
livrable 1 du projet ANR BARPPRO (ANR-10-CSOSG-008), 2011 ;
[13] S. Eveillard, N. Lardjane, J.-Y. Vinçont, I. Sochet, « Towards a fast-running method for
blast wave mitigation by a prismatic blast wall », C.R. Mécanique 341, p. 625–635,
2013 ;
[14] W. Fickett, W. C. Davis, « Detonation », University of CaliforniaPress, Berkeley, Los
Angeles, London, 1979 ;
[15] F.G. Friedlander, « The diffraction of sound pulses », Diffraction by a semi-infinite
plate, Roy. Soc. Lond Q, 186, 322 – 324, 1946 ;
[16] GBPP, Guide SFEPA n°9, « Guide des Bonnes Pratiques en Pyrotechnie », 2009 ;
[17] GFC, « Groupement Français de la Combustion », Dictionnaire des Mots de la
Combustion, www.gfcombustion.asso.fr ;
[18] Green-Book TNO, « Method for the determination of possible damage to people and
objects resulting from releases of hazardous materials», CPR 16E, CIP data of the Royal
Library, First edition 1992 ;
[19] A. Hadjadj, O. Sadot, Shock and blast waves mitigation, Shock Waves 23, pp. 1–4,
2013 ;
[20] Handbook of Chemistry and Physics, 78th Edition, CRC Press, Boca Raton, New-York,
1997-1998 ;
219
[21] F.G. Hulton, « The development of protective walls », 21th Military Aspects of Blast
and Shock, October 4 – 8th 2010 Jerusalem (Israel), 2010 ;
[22] D. Hyde, « Microcomputer Programs CONWEP and FUNPRO, Applications of TM 5-
855-1, Fundamentals of Protective Design for Conventional Weapons (User's Guide) »,
Final rept., 1987 ;
[23] H. Jourdren, « A Hydrodynamic AMR Platform for Multi-Physics Simulations »,
Adaptive Mesh refinement – Theory and Applications, Lecture Notes in Computational
Science and Engineering, volume 41, part III, pp 283-294, 2005 ;
[24] M. J. Kamlet, C. Dickinson, « Evaluation of the Simplified Calculational Method for
Chapman-Jouguet Detonation Pressures on the Basis of Available Experimental
Information », Chemistry of Detonations. III, J. Chem. Phys., 48, 43, 1968 ;
[25] L. Kazandjian, J.F. Danel, « A discussion of the Kamlet-Jacobs Formula for the
Detonation », Propellants, Explosives, Pyrotechnics journal, volume 31, pp 20-24,
2006 ;
[26] G.F. Kinney and J. Kenneth Graham, « Explosive shocks in the air », 2nd edition,
Springer-Verlag, Berlin, 1985 ;
[27] A. Lannoy, « Analyse des explosions air-hydrocarbure en milieu libre : Etudes
déterministe et probabiliste du scénario d’accident », Prévision des effets de
suppression, Bulletin Direct, Etudes et Recherches EDF, A4, 1984 ;
[28] C. Matignon, « Etude de la détonation de deux mélanges stœchiométriques
(CH4/H2/O2/N2 et CH4/C2H6/O2/N2). Influence de la proportion relative des deux
fuels et de la température initiale élevée », Thèse Université de Poitiers, France, 18
Décembre 2000 ;
[29] Ch. Matignon, J.Y. Vinçont, S. Eveillard, « Numerical predictions of blast waves
caused by accidental or intentional detonations of gaseous and condensed explosives in
3D complex geometries », Proc. of the 5th Int. Symp. On Tunnel Safety and Security,
New-York, USA, pp.767-768, 2012 ;
[30] P. Miller, « Towards the modeling of blast loads on structures », thesis for degree of
Master of Applied Science, Department of Civil Engineering, University of Toronto,
2004 ;
[31] B. Mohanty, University of Toronto, Engineering Geosciences Meeting, Toronto, ON,
22 April 2004 ;
[32] C.E. Needham, « Blast waves », Springer-Verlag, New-York, 2010 ;
[33] I. Sochet, « Blast effects of external explosions », 9th international Symposium on
Hazards, Prevention and Mitigation of Industrial Explosions, September 5 – 10th 2010
Yokohama (Japon), 2010 ;
[34] TERA100, http://www-hpc.cea.fr/fr/complexe/tera.htm ;
[35] TM5-1300, U.S. Department of the Army, « Structures to Resist the Effects of
Accidental Explosions », Army TM 5-1300, Navy NAVFAC P-397, AFR 88-22,
Washington, D.C, Departments of the Army, Navy and Air Force, 2008 ;
[36] S. Trélat, « Impact de fortes explosions sur les bâtiments représentatifs d’une
installation industrielle », Thèse Université d’Orléans, France, 7 Décembre 2006 ;
[37] Rapport NRDC, « The US nuclear war plan – a time for change », 2001 ;
220
ANNEXES
ANNEXE 1 : PRESENTATION D’UNE GEOMETRIE DE MERLON (GBPP (2009)) ............................................... 222
ANNEXE 5 : REFLEXION D’UNE ONDE DE CHOC SUR UNE SURFACE CONVEXE ......................................... 226
H2 Merlon
H1 Pile α
Figure A1 μ Représentation schématique des hauteurs d’un merlon placé devant une pile de matière explosive
Merlon
Pile
Figure A2 : Représentation schématique des différentes épaisseurs pour un merlon placé devant une pile de
matière explosive
Longueur L
La longueur du merlon L ne comprend pas les talutages situés aux extrémités du merlon et ne
comprend aucune discontinuité.
223
Les capteurs utilisés sont des capteurs piézoélectriques IPC 113B26, figure A4. Les
caractéristiques techniques des capteurs fournies par le constructeur sont présentées
ci-dessous, figure A5.
δors de l’étude préliminaire du projet ANR BARPPRO (Eveillard (2010)), des capteurs
Kistler 603B ont été utilisés. Ces capteurs ont semblé être un peu plus précis que les capteurs
PCB IPC 113B26 pour l’acquisition des profils de pression dans une zone de forte détente
d’onde de choc (par exemple l’arête aval d’un merlon).
224
Les capteurs piézoélectriques sont connectés à des amplificateurs (figure A6) puis leur
sortie est connectée à la centrale d’acquisition. δ’amplification est réglée directement par la
centrale d’acquisition en fonction de la sensibilité de chaque capteur. Les caractéristiques
techniques des amplificateurs fournies par le constructeur sont présentées sur la figure A7.
De plus pour certains capteurs, la faible dynamique observée se traduit par un bruit important
sur la mesure. Ce bruit induit, malgré les précautions prises, une imprécision dans la lecture
du maximum de surpression.
δa même procédure est reconduite pour les différents historiques de pression d’intérêt.
226
d+R
W : masse de la charge (en kg de TNT de gaz) ;
Paramètres :
La configuration étudiée est une surface convexe placée à une distance d du centre de
la charge (point O), figure A9. δ’obstacle possède un rayon de courbure (R) et l’onde
sphérique divergente est supposée plane au point d’impact (I). δors de l’interaction de l’onde
de choc avec la paroi, l’onde subit une réflexion, figure A10. Ce phénomène physique dépend
227
Maximum de
surpression (bar)
1,80
1,29
0,771
0,257
-0,257
Figure A10 : Visualisation de la propagation d’une onde de choc sur une surface convexe (R = 0,18 m, soit
1,10 m/MJ1/3) pour un obstacle placé à d = 0,12 m (soit 0,73 m/MJ1/3) du centre d’une charge de gaz
(propane-oxygène à la stœchiométrie, E = 4,4.10-3 MJ) (t = 0,28 ms, soit 1,71 ms/MJ1/3)
Possibilités : 1 2 3 4
« Régime à la fin de la Simple réflexion Simple réflexion Simple réflexion Simple réflexion
face avant » à = λ0° de Mach de Mach de Mach de Mach
Figure A11 : Diagramme des réflexions possibles le long d'une paroi convexe
δa réflexion d’une onde de choc sur une surface cylindrique a également été étudiée
par Trélat (β00θ). δ’auteur étudie pour une charge gazeuse (propane-oxygène à la
stœchiométrie) l’interaction d’une onde de choc avec un obstacle cylindrique. A partir de ses
données expérimentales, Trélat (2006) établit des formules empiriques pour estimer
l’évolution du maximum de surpression et ainsi que l’impulsion positive le long de la surface
cylindrique, tableau A1.
228
Caractéristiques de
Formulations de Trélat (2006)
l’onde de choc
P
ΔP+
(en bar)
exp( a b ln( ))
P0
I
exp( A B ln( ))
I+
(en bar.ms) E1/ 3
A, B, a, b : coefficients établis suivant la zone étudiée
P0 : pression atmosphérique (bar)
E : énergie de la détonation de la charge explosive (MJ)
Tableau A1 : Formules empiriques de Trélat (2006)
Cette formulation ne tient cependant pas compte de l’ensemble des paramètres nécessaire à la
description de sa configuration cylindrique : d, , R,… De plus, le phénomène physique de
réflexion est non-linéaire et la formulation étant empirique, il est impossible d’étendre ses
formulations à l’ensemble des surfaces convexes.
Y
d = 0,73 m/MJ1/3 Zone d’étude
Zone d’étude
Onde incidente
X
Centre de la charge (0, 0)
Figure A12 : Représentation schématique de la configuration choisie pour étudier un cas de réflexion sur une
surface convexe
Rayon de la charge de
Distance entre le centre
gaz
Rayon de l’obstacle (R) de la charge et l’obstacle
(en m)
(en m) (d)
(propane-oxygène à
(en m)
stœchiométrie)
0,18 0,041 (E = 4,4.10-3 MJ) 0,12
Approche en géométrie 2D-axisymétrique avec la mise en œuvre d’un maillage
« fin » issu du standard de calcul
Tableau A2 : Surface convexe étudiée par simulation numérique
229
d I*
I R
α R
Charge (W)
α *
w' w
O ξ
* ξ
E F A B C
ξ ξ
X1 X2
Figure A13 : Représentation de la configuration étudiée
La figure A13 présente les différents paramètres permettant de décrire la réflexion sur
une surface convexe : R, , , et d. Deux points d’impact sont ici considérés : I et I*. Le
premier point (I) représente un point d’impact de l’onde de choc sur la surface convexe en vue
directe du centre de la charge (phénomène de réflexion). Le point I* correspond à la fin de la
zone en vue directe de la charge et donc le passage d’un phénomène de réflexion vers un
phénomène de détente du pied de Mach issu de la réflexion sur la surface en vue directe.
A partir des paramètres de la surface d’intérêt (R, , , et d), il est possible d’établir
différentes relations géométriques entre les paramètres.
Angle d’observation ( )
Charge (W) Y
O ξ A C
X1 X2
Figure A14 : Représentation du triangle OCI
Relations géométriques avec les deux triangles (OIA) et (CIA), figure A14 :
X1 X 2 d R (A1)
Résolution du problème :
R 2 X 22 2 X 12 (A5)
R 2 X 22 2 d R X 22
2
(A6)
R 2 2 d R
X2
2 d R
2
(A7)
R 2 2 d R 2
arccos
(A8)
2 R d R
R/d
d (1 R / d ) 1 cos ( ) cos( )
2
2
(A9)
1 R / d
d+R
I*
*
R
Charge (W) *
O
ξ*
C
cos( )
R
*
1 cos 2 cos( * ) cos
2
(A11)
231
α R
Charge (W) Y
w' w
O ξ
E F A B C
Les relations géométriques avec les triangles (ACI), (CBI) et (OBI) permettent de calculer
l’angle d’incidence ( ) en fonction des paramètres géométriques de la configuration,
figure A18 et relation (A12) :
R sin( )
arcsin
(A12)
233
Figure A19 : Evolution de l’angle d’incidence en fonction de l’angle d’observation ( ) et du rapport R/d
Kinney (1985), l’angle de transition ( M) entre une réflexion régulière (RR) et une réflexion
de Mach (RM) est estimé à 42,6° (soit M = 18,6° et M = 1,0 m/MJ1/3), figure A20.
δors de l’interaction d’une onde de choc sur une surface convexe, l’onde subit un
phénomène de réflexion sur la face en vue directe de la charge ( < *). Après ce phénomène
de réflexion, l’onde subit une détente sur le reste de la surface ( > *). Le pied de Mach
formé sur la paroi va se détendre sur la surface non exposée à la charge. Ce phénomène
conduit à une atténuation du maximum de surpression (ΔP+) et de l’impulsion positive (I+) le
long de paroi. δ’angle limite pour la configuration étudiée est de ηγ,1°. A partir de cet angle
limite, l’analyse bibliographique indique qu’il faut envisager un phénomène de détente.
A1.5.2. Résultats et observations
δ’étude de l’interaction de l’onde de choc générée par la détonation d’un mélange
gazeux avec une surface convexe est étudiée par simulation numérique en géométrie
2D-axisymétrique avec un maillage « fin » (standard de calcul).
Capteur de référence
(0, 0) RT
x
Figure A21 : Représentation de la position du capteur de référence, expérience de Trélat (2006)
δ’analyse du capteur de référence permet de comparer les résultats obtenus par simulation
numérique avec la formulation de Trélat (2006), cas d’une configuration en champ libre pour
235
Onde de
choc
a) b)
c)
Figure A23 : Evolution du maximum de surpression réfléchie, profils de surpression et de l’angle d’incidence le long de la paroi
237
δors d’une réflexion régulière (RR) sur une surface convexe, les résultats obtenus pour
le maximum de surpression par simulation numérique se situent dans l’intervalle d’incertitude
bibliographique (+/- 20% par rapport aux estimations du TM5-1300 (2008)). Cet écart peut
être imputé aux choix des conditions physico-numériques de calcul (convergence des
résultats, dimensionnement des arêtes des mailles).
Après la formation du pied de Mach (point (III)), les valeurs du maximum de surpression
réfléchie sont supérieures aux estimations du TM5-1γ00 (β008). Cet écart provient d’un effet
de courbure du front lors de l’interaction de l’onde de choc sur la surface cylindrique que ne
prend pas en compte le TM5-1300 (2008), figure A24 et A25.
b) réflexion de Mach
Onde de choc
a) réflexion régulière
Figure A24 : Représentation schématique de l’interaction d’une onde de choc avec une paroi convexe
≈ λ0o
≈ 80o
≈Cλ0 Maximum de
surpression (bar)
1,34
≈ 80
0,899
0,458
= ζ8o
0,017
μ= ┴
-0,424
= 45o
┴ O
μ
Figure A25 : Visualisation de l’effet de courbure pour la configuration étudiée (t = 0,33 ms)
d
T
h
c
Φ R
Charge (W)
ξc ωc ωc
O Z
F W C
Lors de la propagation d’un pied de εach le long d’une surface, les caractéristiques de
l’onde réfléchie sur la paroi sont obtenues au niveau du point triple (figure A26). Les relations
ci-dessous permettent de calculer l’angle d’incidence ( c) et la distance entre le centre de la
charge et le point triple ( c).
Angle d’incidence ( c)
R h sin
C arcsin
C
(A13)
I’
p
τ
R
τp Φ
M
Charge (W)
O ξp ωp’ ωp
H’ F Z W M
C
Figure A27 : Représentation schématique du point d’impact le long d’une paroi plane équivalente à la surface
convexe
239
La figure A27 propose une équivalence entre une surface convexe et une surface
plane. La méthode considère la paroi plane équivalente comme la tangente au cylindre au
niveau du point M (point de transition entre réflexion régulière et une réflexion de Mach). Les
paragraphes ci-dessous présentent les différentes relations géométriques pour la surface plane
équivalente.
Angle d’incidence ( )
R sin
p M arcsin
cos M p
(A15)
R cos
p R d R d 2
cos M cos M
2
R
(A16)
I’
I’ I
R
M
Charge (W)
C
O ξ ξc
F M
Figure A29 : Diagramme des régimes de réflexion le long d’une surface convexe
Figure A30 : Evolution du maximum de surpression réfléchie par la méthode de la paroi plane équivalente
La méthode proposée ici permet d’estimer rapidement15 la surpression réfléchie à partir d’une
paroi plane, dite « équivalente ». Cette paroi plane correspond à la tangente au point de
transition sur le cylindre (transition entre une réflexion régulière et une réflexion de Mach).
La méthode repose donc sur une hypothèse très forte. Cette approche est une première
ébauche d’une méthode d’estimation qui mériterait d’être poursuivie, afin d’évaluer plus
précisément l’effet de courbure (problème géométrique). En l’état, cette méthode n’est pas
satisfaisante.
En effet, pour des configurations avec un rapport R/d inférieur à 1, l’approche de la paroi
plane équivalente conduit à un écart important dans le calcul de la hauteur du pied de Mach
par une surestimation importante de la surpression réfléchie. Dans ce cas, les distances
parcourues le long de ces deux parois peuvent être différentes pour un même angle
d’observation. Cette méthode n’est donc pas totalement satisfaisante.
Le domaine de validité de cette approximation est limité à des configurations où R/d > 0,75,
ce qui a été vérifié pour des simulations qui ne sont pas présentées ici. Inversement, la
méthode s’est avérée suffisante pour des configurations avec un rapport R/d supérieur à 1.
15
Méthode analytique ne prend que quelques secondes, à l’inverse, la simulation numérique peut prendre plusieurs heures sur
supercalculateur
242
i arctan
L
2 d H 2 e 3 H 4
(A17)
tan1 tan 2
si k ≤ i, k = 1; sinon k = 0
A2.3.2. Onde issue du merlon
δ’estimation de la surpression maximale pour l’onde issue du merlon implique de
connaître le type de réflexion au niveau du point de sortie de la face avant. En fonction de la
position du point d’observation, il est possible de définir deux points de sortie : le point S’
pour le secteur 1 et le point S pour les secteurs 2, 3 et 4, figure A31.
244
Point S x
x
z H
S
x- α1
αS’ xS αS
Poi
sol
nt
Figure A31 μ Représentation des points de sortie S’ et S sur la face avant de la barrière
d tan1
2 tan 2 tan
2 2
L L
S '
cos
(A18)
d tan1
2 tan
L
S ' arctan
(A19)
2 tan cos
L d
d tan1
2 tan
L
S ' S ' arctan
(A20)
L
2
2 sin
Pour le point S, les paramètres géométriques sont obtenus à partir des relations (A21) à (A23).
H
d
tan1
2
S H2
cos
(A21)
S arctan
H 1
H (A22)
tan1 cos
S S arctan
H
d H 1
(A23)
tan1 cos
2
Pour le secteur 1, l’angle de déviation ( w1) pour le calcul du coefficient de détente (Cdw1)
dépend de l’angle d’incidence au niveau du sol ( sol), figure A32. δ’angle d’incidence et le
maximum de surpression incident au sol permettent d’estimer le coefficient de réflexion
245
Γ Point de sortie
Point d’observation
sol (x, z)
sol
Secteur 1 z
k = f(Γ)
δ’angle d’incidence sur le sol ( sol) est calculé à partir de la relation (A24).
d tan1
2 tan
L
sol arctan
2
(A24)
L
x z
2 tan
2 2
2
L
d tan1 z
L
dw1 S ' S ' x
2 tan 2 tan
2 2 2
L L
2
(A25)
Pour le secteur 2, le point de sortie est noté T’. Ce point se situe sur l’arête de la face latérale
du merlon, figure A33-a. Pour les secteurs 3 et 4, le point de sortie T se situe sur l’arête de la
face arrière du merlon, figure A33-b. Les paramètres géométriques au niveau des points T’ ou
T et du point d’observation permettent d’estimer le coefficient de détente (Cdw2) et le
coefficient de réflexion au sol (Crsol).
Point T’ x z
ξ' ξ
ω'
sol
x ω
L
2 tan() H S
H
z sol
z
L
a) Secteur 2 2 b) Secteurs 3 - 4
d
2 tan tan1
L H
'
cos
(A26)
cos
e
(A27)
246
L
2
Secteurs 3 – 4: sol arctan
H
d e
tan1
(A30)
H
2
2
x z cos
2
Les distances parcourues entre le centre de la charge et le point d’observation pour les
secteurs 2, 3 et 4 sont notées respectivement dw2, dw3 et dw4, formules (A31) et (A32).
d e
tan1
2
H
Pour les secteurs 3 et 4 : dw3 dw4 S H 2 x 2 z 2
cos
(A32)
Dans le secteur 3-bis, l’onde issue du merlon est soumise à une troisième détente sur l’arête
de la face arrière du merlon. Le calcul des paramètres géométriques en aval est donc différent
de celui du secteur 3. Le point de sortie V se situe sur l’arête de la face arrière, figure A34.
z
Point S z
Point V
H – h’
h' Point T
x
sol
Point V l'
1 1
Secteur 3-bis :
x xd H e
tan (α1 ) tan (α2 )
L
l ' d e d e tan
2 2
Ces distances permettent de calculer l’angle de déviation de la seconde détente ( w2) ainsi que
la distance ( ) entre le point T et le point V, formules (A35) et (A36).
h'
w2 arctan
l'
(A35)
Le nombre de Mach incident pour la troisième détente (Cdw3) est obtenu à partir des
coefficients établis précédemment (Crs, Cdw1, Cdw2) et de la distance parcourue entre le centre
de la charge et le point V ( V), formule (A37).
V T (A37)
L
2
La distance parcourue entre le centre de la charge et le point d’observation situé dans le
secteur 3-bis est notée dw3-bis, formule (A39).
Pour le secteur 3-bis : dw3bis v H h'2 x 2 y 2 L (A39)
2 2 tan
2 2
L
248
Sébastien Eveillard