Chapitre Un

émondetal15
Résistance des structures à l’impacte d’un projectile 1

Table des matières
Liste des symboles 7
Dédicace 8
Remerciement 9
Résumé 10
abstrat 11
Introduction générale 12
I Dossier technique 14
1 Étude bibliographique 15
1.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
1.2 Projectile . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
1.2.1 Présentation du projectile . . . . . . . . . . . . . . . . . . . . . . . 15
1.2.2 Choix et caractéristique du projectile . . . . . . . . . . . . . . . . . 16
1.3 Classification des impacts . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
1.4 Impacte du projectile sur une structure en béton . . . . . . . . . . . . . . . 19
1.5 Modèles de calcul en dynamique rapide . . . . . . . . . . . . . . . . . . . . 20
1.5.1 Méthode des éléments discrets . . . . . . . . . . . . . . . . . . . . 20
1.5.2 Méthode des éléments finis . . . . . . . . . . . . . . . . . . . . . . . 21
1.6 Modèles de comportement des structures sous Impact . . . . . . . . . . . . 22
1.6.1 Modèle d’endommagement . . . . . . . . . . . . . . . . . . . . . . . 22
1.6.2 Modèles de plasticité . . . . . . . . . . . . . . . . . . . . . . . . . . 25
1.6.3 Modèles couplés élasto-endommageable-plastique . . . . . . . . . . 27
1.7 Force du projectile . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28
1.8 Calcul de la résistance des structures (Béton)aux impacts de projectiles . 31
1.8.1 Généralité . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 31
1.8.2 Corrélations Empiriques . . . . . . . . . . . . . . . . . . . . . . . . 31
1.8.3 Formule génériques (Béton/Acier) pour impact localisé . . . . . . . 32
1.8.3.1 Impact à vitesse supersonique . . . . . . . . . . . . . . . . 32
2
TABLE DES MATIÈRES
1.8.3.2 Impact à vitesse subsonique . . . . . . . . . . . . . . . . . 32

1.9 Conclusion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 34
2 Simulation Numérique 35
2.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 35
2.2 Mise en œuvre des simulations . . . . . . . . . . . . . . . . . . . . . . . . 35
2.2.1 Présentation de l’outil de simulation . . . . . . . . . . . . . . . . . 36
2.2.2 Principaux paramètre de calcul d’impact . . . . . . . . . . . . . . . 36
2.2.2.1 généralité . . . . . . . . . . . . . . . . . . . . . . . . . . . 36
2.2.2.2 Caractéristique du projectile . . . . . . . . . . . . . . . . . 36
2.2.2.3 Caractéristiques de l’objet impacté (paroi) . . . . . . . . . 37
2.2.2.4 Définition des paramètres pour la loi [8] . . . . . . . . . . 38
2.2.2.4.1 Détermination de la déformation équivalente ε̄ . . 38
2.2.2.4.2 Paramètre d’évolution de l’endommagement en
traction et en compression . . . . . . . . . . . . . 38
2.3 Analyse des résultats . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 39
2.3.1 Les effets du projectile sur la paroi . . . . . . . . . . . . . . . . . . 39
2.3.1.1 MAILLAGE . . . . . . . . . . . . . . . . . . . . . . . . . 40
2.3.1.2 Chargement appliqué . . . . . . . . . . . . . . . . . . . . . 40
2.3.1.3 Fissures . . . . . . . . . . . . . . . . . . . . . . . . . . . . 40
2.3.2 Conclusion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 41
II DOSSIER PÉDAGOGIQUE 42
3 Rapport du stage pédagogique 43
3.1 INTRODUCTION . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 43
3.2 Présentation de l’école d’accueil . . . . . . . . . . . . . . . . . . . . . . . . 44
3.2.1 Historique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 44
3.2.2 Situation géographique . . . . . . . . . . . . . . . . . . . . . . . . . 44
3.2.3 Organisation administrative (organigramme) . . . . . . . . . . . . . 44
3.2.3.1 L’organe suprême de gestion scolaire . . . . . . . . . . . . 44
3.2.3.2 Association des parents d’élèves . . . . . . . . . . . . . . . 45
3.2.3.3 La direction . . . . . . . . . . . . . . . . . . . . . . . . . . 45
3.2.3.4 Le conseil des maitres . . . . . . . . . . . . . . . . . . . . 45
3.2.4 Organisation Pédagogique . . . . . . . . . . . . . . . . . . . . . . . 45
3.2.4.1 Groupe pédagogique . . . . . . . . . . . . . . . . . . . . . 45
3.2.4.2 Pédagogie et Didactique appliquées . . . . . . . . . . . . . 45
3.2.4.3 Les filières de formation . . . . . . . . . . . . . . . . . . . 45
3.2.4.4 Les emplois du temps . . . . . . . . . . . . . . . . . . . . 46
3.2.4.5 Les effectifs . . . . . . . . . . . . . . . . . . . . . . . . . . 46
3.2.4.6 Le fonctionnement des classes . . . . . . . . . . . . . . . . 46
3.3 Activités menées . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 46
3.3.1 Activités pédagogiques . . . . . . . . . . . . . . . . . . . . . . . . . 46
3.3.2 Activités parascolaires . . . . . . . . . . . . . . . . . . . . . . . . . 46
3.4 Difficultés rencontrées . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 47
3.4.1 Matérielles . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 47
3.4.2 Pédagogiques . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 47
3.4.3 Socioculturelles . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 47

TABLE DES MATIÈRES
3.5 Solutions apportées . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 47

3.5.1 Initiatives personnelles et apport au d´développement de l’école . . 47
3.5.2 Autres . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 47
3.6 Les défis et perspectives . . . . . . . . . . . . . . . . . . . . . . . . . . . . 48
3.7 Leçons tirées du stage . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 48
3.7.1 Les intérêts . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 48
3.7.2 Les motivations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 48
3.7.3 Points de satisfactions . . . . . . . . . . . . . . . . . . . . . . . . . 48
3.8 Les recommandations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 48
3.8.1 Points de convergence entre le stage et la formation à l’ENETP . . 48
3.8.2 Points de divergences entre le stage et la formation à l’ENETP . . . 48
3.8.3 Besoins en formation pour rendre mon travail plus efficace . . . . . 49
3.8.4 Autres remarques et suggestions . . . . . . . . . . . . . . . . . . . . 49
3.9 La leçon modèle . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 50
3.9.1 Partie administrative . . . . . . . . . . . . . . . . . . . . . . . . . 50
3.9.2 Document élèves . . . . . . . . . . . . . . . . . . . . . . . . . . . . 52
3.9.3 Evaluation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 53
3.9.4 Réponses . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 53
3.10 Conclusion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 54
Conclusion générale 55

Liste des tableaux
1.1 Paramètre de perforation selon la nature du matériau . . . . . . . . . . . . 33

1.2 Épaisseurs de perforation en fonction des caractéristiques des projectiles . . 33
1.3 Paramètre de perforation selon la nature du matériau . . . . . . . . . . . . 34
2.1 classe du béton . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 37

2.2 Valeurs des paramètres de la loi [8] . . . . . . . . . . . . . . . . . . . . . . 39
2.3 Définition du maillage . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 40
3.1 Récapitulatif des effectifs d’élèves par filières . . . . . . . . . . . . . . . . . 46
5
Table des figures
1.1 Modèle mécanique simple d’un impact à travers un système de masse-

ressort(Eibl,1987) (CEB,1988) . . . . . . . . . . . . . . . . . . . . . . . . . 17
1.2 Impact mou (CEB,1987) . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18
1.3 impact dur (Eibl,1987) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18
1.4 classification des impacts (Kœchlin et Potapov, 2009) . . . . . . . . . . . . 19
1.5 Echelle de modélisation (D’addetta,et al.,2001) . . . . . . . . . . . . . . . . 21
1.6 (Comportement uniaxial du béton lors d’un cycle chargement – déchargement
a)( modèle d’endommagement, b) modèle de plasticité avec l’écrouissage
négatif . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27
1.7 (Comportement du béton sous a)( la traction uniaxiale et b) la compression
uniaxiale (Voyiadjis, et al., 2008) . . . . . . . . . . . . . . . . . . . . . . . 27
1.8 Comparaison entre résultats expérimentaux et de simulation pour un mor-
tier sous compression (Burlion, et al., 1998) . . . . . . . . . . . . . . . . . 28
2.1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 38
2.2 Zone de localisation des déformations . . . . . . . . . . . . . . . . . . . . 39
2.3 Mazars : Comportement en traction-compression. . . . . . . . . . . . . . . 39
2.4 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 40
2.5 chargement surfacique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 40
2.6 Propagation des fissures . . . . . . . . . . . . . . . . . . . . . . . . . . . . 41
2.7 Amplitude du déplacement . . . . . . . . . . . . . . . . . . . . . . . . . . 41
3.1 Organigramme de l’établissement. . . . . . . . . . . . . . . . . . . . . . . 44
6
Liste des symboles
g acceleration
hjj oppou
7
Dédicace
Je dédie ce modeste travail :
A mes chers parents qui m’ont guidé durant les moments les plus difficiles de ce long
chemin. Vos prières et vos bénédictions m’ont été un grand secours pour mener à bien
mes études. Aucun mot, aucune dédicace ne saurait être assez éloquente pour exprimer
ce que vous méritez pour tous les sacrifices que vous n’avez cessé de me donner depuis
ma naissance, durant mon enfance et même à l’âge adulte. Je vous dédie ce travail en
témoignage de mon profond amour.
8
Remerciement
Tout d’abord gloire au bon Dieu qui m’a donné la force et le courage d’accomplir ce
modeste travail.
Je tiens à remercier DR. Mamadou Toungara professeur de génie civil à l’école nor-
mal d’enseignement technique et professionnel (l’ENETP) et directeur de ce mémoire,
pour l’aide qu’il a fournie et les connaissances qu’il a su me transmettre. Je le remercie
également pour sa disponibilité et la qualité de ses conseils.
Je souhaite également remercier les membres du jury pour leur présence, pour leur
lecture attentive de mon mémoire ainsi que pour les remarques qu’ils m’adresseront lors
de cette soutenance afin d’améliore mon travail.
Je remercie également toute l’équipe pédagogique du Centre de formation en admi-

nistrations et aux métiers (CFAM), pour tous les conseils qu’elle m’a donnée lors de mon
stage pédagogique.
J’exprime ma gratitude à l’ensemble des enseignants du département génie civil de

l’ENETP qui, par leur parole, leurs écrits, leurs conseils et leurs critiques ont contribué à
la réalisation de ce travail.
Je souhaite également à remercier le directeur des études de l’ENETP, Dr. Yacouba

Daou.
Mes sincères remerciements à tous mes camarades de classe de l’ENETP, pour leur
soutien inconditionnel et leurs encouragements qui ont été d’une grande aide.
Je souhaite également remercier ma femme Gnia Traoré et toutes les membres de sa

famille.
Enfin, je remercie profondément tous mes parents, pour leur soutien constant et leurs
encouragements.
9
Résumé
Le travail présenté dans ce mémoire est constitué d’un dossier technique et d’un dos-
sier pédagogique.
Le dossier technique vise à mieux comprendre le comportement dynamique de structures
en béton armé lors d’un impact de projectile. Pour ce faire, une simulation numérique, par
élément finis, a été réalisé sur une paroi en béton armé avec deux couches d’acier sous l’im-
pact d’un projectile à une force de 42150Kn. Les résultats obtenus ont permis de prédire
le comportement dynamique de la structure à étudier. la simulation est modéliser par la
loi de comportement de mazars. De ce fait une analyse des résultats est effectué pour nous
montrer les dommages causés par le projectile sur la paroi impacté,des déformations et
l’apparition des fissures dans la zone d’impact.
Le dossier pédagogique est constitué d’un rapport de stage pédagogique avec une leçon
modèle. Le stage a été effectué au Centre de Formation en administration et aux métiers
(CFAM) au titre de l’année scolaire 2018-2019.La leçon modèle a porté sur l’assainisse-
ment, dispensée dans une classe de 3ème Année Dessin Bâtiment dans le même établissement.
L’objectif de ce stage était de lier la théorie à la pratique.
10
abstrat
The work presented in this thesis consists of a technical file and a pedagogical file.
The technical file aims to better understand the dynamic behavior of reinforced concrete
structures during a projectile impact. To do this, a numerical simulation, by finite element,
was carried out on a reinforced concrete wall with two layers of steel under the impact of
a projectile at a force of 42150Kn. The results obtained allowed to predict the dynamic
behavior of the structure to be studied. The simulation is modeled by the mazars behavior
law. Thus an analysis of the results is made to show us the damage caused by the projectile
on the impacted wall, deformations and the appearance of cracks in the impact zone. The
pedagogical file is constituted by a pedagogical internship report with a model lesson.
The training course was carried out in the Center of Training in Administration and in
the trades (CFAM) for the school year 2018-2019.
The model lesson was on sanitation, taught in a class of 3rd Year Design Building in the
same institution.
The objective of this course was to link theory to practice.
11
Introduction générale
Le présent document est subdivisé en deux grandes parties :

Une première partie consacrée au dossier technique portant sur la résistance des
structures à l’impact d’un projectile.
Une seconde partie présentant un rapport de stage pédagogique avec une leçon
modèle.
Le béton est le matériau le plus utilisé dans les ouvrages de génie civil (bâtiments,
ponts, barrages, centrales nucléaires,. . . ), ils présentent des propriétés spécifiques liés no-
tamment à l’hétérogénéité de ce matériau, à sa faible résistance en traction (comportement
fragile en traction), à sa faible déformation à la rupture et enfin à sa ductilité en com-
pression.
Fondamentalement, le béton est un matériau hétérogène, mais il est considéré comme un
matériau homogène pour faciliter le calcul. Par conséquent, pour soutenir la bonne qualité
du béton, la production du béton nécessite un bon procédé de fabrication. L’effondrement
du mode de rupture en raison de fissures est un fait important qui se produit dans les
éléments avec des matériaux à base de ciment. La complexité des problèmes qui se posent
en raison de la présence de fissures qui se développent et se propagent aléatoirement, suivis
par des fissures ultérieures. Ce motif de déformation est important en tant que marqueur
de sécurité d’une structure.
Il serait intéressant de passer en revue quelques résultats tirés de l’abondante littérature
qui traite les principaux aspects du comportement mécanique de ce matériau (Mazars et
Bazant, 1988 ;...).
L’une des caractéristiques du béton est un comportement différent en traction et en com-
pression : il est beaucoup plus résistant en compression. C’est la raison pour laquelle
il est souvent considéré comme ayant une résistance nulle en traction lors des calculs
réglementaires de structures. Ce comportement, en traction et en compression, est for-
tement influencé par les mécanismes internes d’altérations qui conduisent à une perte
de raideur. Il est certain que tenir compte de ces mécanismes permettrait de décrire
l’évolution progressive des dégradations dans ce matériau.
Une autre caractéristique du comportement du béton est son caractère unilatéral : lors du
passage d’une sollicitation endommageable en traction à une sollicitation en compression,
macroscopiquement on observe une restauration de la raideur du matériau.
L’effet d’échelle est une autre propriété du béton qui est essentielle en génie civil. Celui-
ci se manifeste par une dépendance de la contrainte nominale maximale, couramment
évaluée à partir de la charge maximale supportée par une structure, vis-à-vis de l’échelle
des structures.
Les observations expérimentales montrent que cette contrainte diminue avec la taille de
la structure (Bazant et Ozbolt, 1990 ; Mazars, Pijaudier-Cabot et Saouridis, 1991 ; Biolzi
et al, 1994).
Ces quelques observations témoignent de la complexité du comportement mécanique du
12
TABLE DES FIGURES
béton dont nous nous apprêtions à étudier une structure à l’impact d’un projectile.
L’objectif de cette étude est de faire une simulation numérique sur une structure porteuse
d’un site industriel suite à une collision de projectile pour voir comment les fissures vont
se propager.
Le plan suivant est proposé pour ce manuscrit :
La première partie comprend deux (2) chapitres :

— Le chapitre 1 présente une étude bibliographique sur le comportement mécanique
du béton sous impact et les différents modèles de comportement capables de re-
produire ces comportements du béton lors de l’impact.
— Enfin le chapitre 2 présente une modélisation numérique capable de simuler le
comportement du béton sous chargement dynamique avec le logiciel COMSOL.
La seconde partie comprend :

— Un rapport de stage pédagogique ;
— Une leçon modèle

Première partie
Dossier technique
14
Chapitre 1
Étude bibliographique
1.1 Introduction
Les structures en béton armé sont les ouvrages les plus fréquents pour la construction
d’ouvrages en Génie Civil. Cependant la connaissance de ses comportements en dyna-
mique rapide reste restreinte. Le béton est un matériau issue d’un milieu hétérogène, son
comportement est complexe à analyser. Plusieurs auteurs ont proposé des lois de compor-
tement avec un degré de complexité plus ou moins élevé pour prédire son comportement
en dynamique rapide.
L’objectif de ce chapitre est de présenter principalement les travaux de recherche sur le
comportement des structures sous l’impact et les modèles numériques capables de repro-
duire ces phénomènes.
1.2 Projectile
1.2.1 Présentation du projectile
Bien avant l’invention des armes à feu, les romains savaient produire des balles de
plomb percées de deux trous (qui sifflait en effrayant l’ennemi). Mais à partir de la fin
du XIXe siècle et pendant le début du XXe siècle, d’énormes progrès ont été faits sur le
matériel d’artillerie aussi bien en ce qui concerne les projectiles, la portée, la précision ou la
rapidité des tirs. Un projectile est un engin explosif conçu pour détruire ou endommager les
structures ou bien pour blesser, tuer ou limiter d’une autre manière les activités humaines.
Il peut provenir d’une arme de jet ou d’une arme à feu, ou bien être lancé à main nue
comme une pierre, un pavé, un marteau ou une lance.
Il existe plusieurs types de projectile à savoir :
- Lance-roquette
- Missile
- Bombe
- Grenade
- Cartouche
- Obus
- Balistique
- Cas accidentel (percutions de voiture, Avion)
- Fragments après l’explosion
15
CHAPITRE 1. ÉTUDE BIBLIOGRAPHIQUE
Ces projectiles sont capables de nuire à la résistance de la cible visée (structure de

bâtiment) .
1.2.2 Choix et caractéristique du projectile

Généralement la règle utilisée de nos jours veut que les :
- engins possédant un guidage soient nommés missiles, quel que soit le système de
propulsion ;
- engins propulsés par un autre moyen qu’une fusée à poudre soient nommés missiles,
à part quelques prototypes datant des environs de la Seconde Guerre mondiale tous
ces engins ont un système de guidage ;
- fusées à poudre sans guidage soient nommées roquettes ;
- projectiles (guidés ou non) se déplaçant sous la surface de l’eau soient nommés
torpilles.
Cependant, il y a des exceptions, tels que les projectiles des lance-roquettes multiples,
qui sont actuellement le plus souvent autoguidés tout en conservant le nom de roquettes,
ou des prototypes datant d’une période où les systèmes électroniques étaient bien plus
coûteux, fragiles et volumineux qu’actuellement. Une telle utilisation de ce terme est
exceptionnelle et, en général, due à un contexte historique particulier.
D’après cette règle le choix portera sur un engin explosif autopropulsé et guidé.
Par définition Un missile est une arme propulsée et guidée conçue pour emporter une
charge militaire. C’est un engin qui lance une ou plusieurs armes en leur donnant une
trajectoire essentiellement balistique. (c’est-à-dire la trajectoire est influencée uniquement
par la gravité et la vitesse acquise par une force d’accélération initiale).
1.3 Classification des impacts

Un projectile est un corps lancé ou projeté pour atteindre une cible, il peut pro-
venir d’une arme de jet (flèche) ou d’une arme à feu (balle/calibre inférieur à 20 mm,
obus/calibre supérieur à 20 mm. . . ), ou bien être lancé à main nue comme une pierre, un
pavé, une dague, un marteau ou une lance.
Il existe deux grands types d’armes de tire qu’on utilise généralement dans les champs de
bataille qui sont :
— armes qui font des tires directes dont la trajectoire est presque rectiligne ;
— armes qui font des tires indirectes dont la trajectoire suit une courbe.
L’importances des impacts sont fonction de la distance que le projectile parcourt pour
atteindre la cible, plus on s’approche de la cible le choc devient grave .
Lors de l’impact entre un projectile et une structure en béton, la structure cible est sou-
mise à différents trajets de chargement, complexes, difficiles à contrôler et dépendant du
type d’impact. La complexité des effets d’un tel impact amène à considérer deux cas li-
mites, impact mou et impact dur, pour l’étude des différents phénomènes physiques liés
à chaque type d’impact. Une fois le type d’impact déterminé, les modes de rupture cor-
respondant de la structure peuvent être prédits pour le dimensionnement de la structure.
La classification des types d’impacts n’est pas facile. Il n’y a pas beaucoup de travaux de
recherche concernant ce problème.
[1] et [2] ont donné une définition qualitative des impacts mous et durs. L’impact étudié
est la collision de deux corps, l’un avec une vitesse initiale vient frapper un autre étant
au repos. L’objet soumis à l’impact est généralement une structure qui doit être conçue

contre les chocs. Ce problème peut être réduit à une collision de deux masses m1 et m2.
Un ressort de contact avec une rigidité k1 est placé entre les deux masses pour simu-
ler la force contact. Un autre ressort avec une raideur k2 représente la relation entre la
déformation et la force résistante de la structure.
Figure 1.1 – Modèle mécanique simple d’un impact à travers un système de masse-
ressort(Eibl,1987) (CEB,1988)
Selon cet exemple de calcul, les équations différentielles représentant les mouvements
des 2 masses-ressorts sont écrites de manière suivante :
m1 ẍ1 (t) + K1 [x1 (t) − x2 (t)] = 0
m1 ẍ2 (t) + K1 [x1 (t) − x2 (t)] + K2 x2 = 0 (1.1)

Si X1 X2 , autrement dit, si la déformation du projectile est beaucoup plus importante
que celle de la structure soumise à l’impact, et en imposant :
F (t) = k1 x1 (1.2)
Les équations [1-1] sont réécrites en [1-3] :
m1 ẍ1 (t) + K1 x1 (t) = 0(a)
m2 ẍ2 (t) + K2 x2 (t) = F (t)(b) (1.3)

Le problème est maintenant découplé, il est possible de résoudre la première équation
[1-3a] indépendamment, déduire ensuite la force F (t) et finalement déterminer la déformation
de la structure à partir de l’équation [1-3b].
Dans ce cas la structure résistante reste indéformée après l’impact. Par conséquent,

Figure 1.2 – Impact mou (CEB,1987)
l’énergie cinétique de l’impacteur déformable est complètement transférée en énergie de

déformation de l’impacteur. Ce cas est appelé impact mou.
L’autre cas limite correspond à X1 X2 , c’est-à-dire à un projectile rigide, ce cas est
appelé l’impact dur par (Eibl,1987). L’énergie cinétique du projectile est complètement
ou partiellement absorbée par la déformation de la structure soumise à l’impact.
Figure 1.3 – impact dur (Eibl,1987)
D’après (Eibl,1987) et (CEB,1988) un autre auteur (Brossard,1997) trouve une autre

méthode simple pour expliquer le phénomène d’impact. Selon Brossard, en considérant
un système de deux masses ponctuelles qui se choquent, s’il n’y a pas de rebond les deux
masses se collent ensemble et deviennent une seule masse ponctuelle après le choc, ce cas
est appelé impact mou. Au contraire, s’il y a un rebond parfait sans dissipation d’énergie,
ce cas est appelé impact élastique, parfois improprement impact dur. Cette définition est
simple et peut également être appliquée au cas de l’impact d’avion sur une structure en
béton.
Sur un exemple de l’impact à l’échelle réelle réalisé par (Sugano, et al., 1993) ont eu à
comparer les classifications d’impacts définies par Brossard et Eibl, avec lequel il n’y a pas
eu de rebond. Cet exemple correspond à un impact mou selon Brossard, mais comme le
déplacement de la structure cible n’est pas négligeable dans ce test, il est considéré comme
un impact dur selon Eibl ! Il est clair pour ce cas que le déplacement de la structure cible

n’est pas un élément clé pour bien distinguer les types d’impacts.
Pour pallier cet inconvénient,( Kœchlin et Potapov, 2009) ont proposé une autre définition
pour distinguer impact mou et impact dur,prenant en compte à la fois les caractéristiques
du matériau et la vitesse du projectile, valable dans le cas où le projectile s’écrase. Selon
cette méthode, elle peut être utilisée pour classifier le type d’impact entre deux corps.
L’idée est de comparer le seuil de rupture de la cible σc avec celui du projectile σp et
de vérifier si la cible va résister à l’impact. Si la cible résiste et le projectile s’écrase, on
considère que c’est un choc mou, et si le projectile pénètre la cible, c’est un choc dur. Lors
de l’impact, la cible est soumise à une contrainte qui a deux composantes, l’une venant
du matériau, l’autre de la vitesse du projectile :
σ = σp + ρp ν02 (1.4)
Avec σp : la contrainte de rupture du projectile,

ρp : la masse volumique du projectile,
ν0 : la vitesse du projectile.
Notons également que σc est la contrainte de rupture du matériau qui constitue la cible,
la limite entre choc dur et choc mou s’établit lorsque :
σc = σp + ρp ν02 (1.5)
C’est-à-dire :
σp ρp ν02
+ =1 (1.6)
σc σc
Selon le critère défini par la formule [1-6], si le terme gauche est inférieur à 1, le cas
est qualifié d’impact mou, et d’impact dur s’il est supérieur à 1.
Figure 1.4 – classification des impacts (Kœchlin et Potapov, 2009)
1.4 Impacte du projectile sur une structure en béton

En cas d’explosion d’une des installations d’un site industriel, les installations voisines
peuvent être impactée par les fragments générés et projetés aux alentours dues à l’explo-
sion tels que des canalisations, des réservoirs atmosphériques ou sous pression ou encore
des structures porteuses.
Les explosifs (projectiles) peuvent les pénétrer partiellement ou les perforer entièrement.
Ce faisant, ces derniers peuvent, à leur tour, constituer des menaces pour les habitants à

un enchaı̂nement de processus catastrophiques successifs appelé effets dominos. L’acciden-

tologie montre que ce phénomène peut être à l’origine de la rupture d’un grand nombre
d’installations et causer de nombreuses victimes. Le risque de ruine de tels équipements
ou l’endommagement d’éléments de structures suite à un impacte de projectile n’est pas
seulement fonction de la vitesse d’impact du projectile, mais également de sa masse, de
ses caractéristiques matériaux comparativement à celles de la cible. Les effets d’endom-
magement suite à l’impact d’un projectile sur une cible en béton peuvent ainsi se traduire
suivant 7 grands phénomènes :
Pénétration : perçage d’un tunnel dans la cible par le projectile, la longueur
de ce tunnel est nommée profondeur de pénétration ;
Fissuration en cône : formation d’un cône lié au cisaillement dans la cible ;
Caractérisation : éjection de petits morceaux de matériau en face avant de la
cible ;
Fissuration : des fissures radiales apparaissent sur l’une ou les deux faces de la
cible et se propagent dans l’épaisseur de la cible ;
Écaillage : éjection des morceaux de matériau en face arrière ;
Perforation : passage complet du projectile le à travers la cible avec ou sans
vitesse résiduelle ;
Réponse globale de la structure : la vibration globale d’une structure peut
parfois induire des endommagements en cisaillement sur toute la cible.
1.5 Modèles de calcul en dynamique rapide

Plusieurs modèles de calculs en dynamique rapide existent afin de résoudre des problèmes
d’impact ou d’explosion. Les différentes méthodes présentés sont :
— méthodes des éléments de frontière

— méthodes sans maillage
— méthodes des éléments finis
— méthodes des éléments discrets
— méthodes basées sur la partition de l’unité
Aujourd’hui, deux types d’approches sont utilisées :
les approches dites continues et celles dites discrètes . La méthode des éléments fi-
nis correspond au premier type d’approches et, largement répandue de nombreux modèles
de comportement des matériaux ont été développés pour ce type d’approches. Fondée sur
l’hypothèse d’un milieu continu elle peut avoir des difficultés à traiter les fortes disconti-
nuités pouvant apparaı̂tre lors d’un impact (fragmentation).
Les méthodes discrètes, non fondées sur l’hypothèse d’un champ de déplacement
continu sont de plus en plus utilisées pour les problèmes de dynamique rapide.
1.5.1 Méthode des éléments discrets

Les approches discrètes, comme la méthode des éléments discrets (Cundall et Strack,
1979) reposent sur une description du milieu en un réseau de points mobiles munis de
masses associés et des géométries simples, par exemple des sphères. Les lois d’interactions
entre les points dépendent du milieu modélisé (contact, loi cohésive, . . . ).
Elles sont particulièrement adaptées aux milieux granulaires mais peuvent également être

utilisées pour traiter des milieux continus, l’apparition des discontinuités (fissure) dans
ces milieux est aisée à mettre en œuvre car il suffit de supprimer l’interaction entre les
points concernés. A l’échelle macroscopique, une structure peut être supposée homogène
et continue, mais son comportement à la rupture résulte de phénomènes discontinus (
[Figure 1-5]). Dans de tels cas, les méthodes discrètes sont plus adéquates que les ap-
proches continues classiques pour modéliser et reproduire les phénomènes. Les approches
Figure 1.5 – Echelle de modélisation (D’addetta,et al.,2001)
aux éléments discrets ont initialement été introduites par (Cundall et Strack, 1979), pour
les matériaux granulaires. Les développements qui ont suivi ont eu pour but d’adapter ce
type de méthodes aux matériaux cohésifs et en particulier les géomatériaux (D’addetta,et
al.,2001). Ces méthodes sont difficiles à mettre en œuvre pour traiter des structures de
grandes tailles. En effet, pour prédire finement le comportement local, ces méthodes
doivent utiliser beaucoup d’éléments discrets ce qui conduit, à l’échelle de l’ouvrage, à
un nombre très important d’éléments et à des temps importants de calculs.
1.5.2 Méthode des éléments finis

La méthode des éléments finis est largement utilisée pour les calculs d’ouvrages. Elle
est fondée sur la mécanique des milieux continus. En considérant que le comportement du
matériau est homogène, cette approche permet d’étudier la réponse globale de la structure
et la réponse locale du matériau. Pour modéliser le comportement à la rupture du béton,
deux types de modèles ont été développés :
Approche continue de l’endommagement : les phénomènes non-linéaires sont traités
à travers des variables d’endommagement, qui varient de 0 (pour le matériau sain)
à 1 (pour le matériau complètement endommagé) (Mazars, 1984).
Théorie de la plasticité : le comportement non-linéaire est pris en compte par
le biais de déformations plastiques qui respectent des lois d’écoulement et des
critères en contrainte. Cependant, lorsque la propagation des fissures joue un rôle
prépondérant, les lois de comportement ne peuvent pas reproduire correctement
les réponses des structures et les endommagements des matériaux.
Pour résoudre ce problème de discontinuité, différentes solutions ont été proposées
(Wœstyn, et al., 2006) :
pour suivre l’évolution d’une fissure, (Ortiz et Pandolfi, 1999) introduit une dis-
continuité sur l’interface reliant deux éléments finis couplée à un algorithme de

maillage adaptatif. Cette méthode introduit une dépendance au maillage utilisé.

une autre technique consiste à introduire la discontinuité en déplacement direc-
tement à l’intérieur des éléments finis (X-FEM (Extended Finite Element Me-
thod)) (Daux et al., 2000). La propagation d’une fissure en quasi-statique est bien
modélisée par cette méthode mais elle permet difficilement de traiter la multi-
fissuration.
le processus d’érosion est une autre méthode souvent utilisée pour le traitement
des problèmes dynamiques d’impact. Un élément fini de la structure sera supprimé
lorsqu’un critère est atteint.
Différents critères d’érosion peuvent être utilisés, souvent le seuil de déformation plastique
de Von-Mises (Agardh et Laine, 1999) est retenu. Certains modèles élasto-endommageables
utilisent un seuil de déformation principale en traction, i.e. le critère du modèle de (Mazars,
1984). Deux méthodes d’intégration temporelle sont utilisées dans les codes aux éléments
finis, la méthode d’intégration implicite (de type Newmark) et la méthode d’intégration
explicite (en général des différences centrées d’ordre 2). Si l’algorithme implicite est sou-
vent utilisé pour modéliser les problèmes quasi-statiques et ceux de dynamique lente ; les
problèmes de dynamique rapide sont généralement traités avec une méthode explicite qui
nécessite d’utiliser un très petit pas de temps d’intégration et qui est plus efficace en
termes de temps de calculs pour les chargements impulsifs (courts).
1.6 Modèles de comportement des structures sous

Impact
Dans ces dernières années plusieurs programmes de recherche ont été réalisés pour le
but d’étudier et valider le comportement des structures en béton armé ou non, leur objec-
tif et de classer par ordre de pertinence les modèles dans le cadre de calcul des structures
soumis à des chargements dynamiques. Les modèles de comportement des structures sous
impact dépendent du matériau qu’on utilise pour réaliser l’ouvrage. Pour ce présent tra-
vail, plusieurs modèles de comportement du béton sont présentés. Avec l’hétérogénéité du
béton, le comportement est complexe et fortement dépendant des conditions de charge-
ment. Il n’existe pas un modèle universel avec un nombre limité de paramètres capable de
reproduire tous les phénomènes physiques, les modèles existants sont donc généralement
spécifiques aux conditions de chargement. Les modèles les plus connus sont les modèles
d’endommagement et les modèles élasto-plastiques, Chaque modèle présente des avan-
tages ainsi des inconvénients pour la simulation des phénomènes complexes apparaissant
lors d’un impact. Le couplage de ces deux classes de modèles permet de couvrir un large
spectre de phénomènes physiques.
1.6.1 Modèle d’endommagement

L’endommagement des structures est dû à la création et la propagation de micro-
fissures sous l’action du chargement. L’endommagement se manifeste par la diminution des
caractéristiques élastiques et le phénomène d’adoucissement des contraintes. Les modèles
numériques basés sur la mécanique de l’endommagement continu sont développés pour
simuler ces phénomènes. Un point commun à tous ces modèles est l’introduction d’une ou
plusieurs variables locales, notées souvent D, décrivant l’évolution de l’endommagement à
partir de la variation relative de raideur du matériau. Les modèles existants peuvent être

distingués en fonction de l’ordre tensoriel des variables d’endommagement utilisées. En

général, un modèle d’endommagement est construit à l’aide des équations constitutives
suivantes :
potentiel d’énergie libre :
1
ρψ = (1 − D)ij Eijkl kl (1.7)
2
Lois d’état correspondant :
∂(ρψ )
σij = = (1 − D)Eijkl kl (1.8)
∂ij
∂(ρψ ) 1
Y =− = ij Eijkl kl (1.9)
∂D 2
Où ij sont les contraintes associées aux déformations ij, et Y le taux de restitution de
l’énergie associé à l’endommagement D. Le taux de dissipation d’énergie φ̇ doit satisfaire
le second principe de la thermodynamique (inégalité de Clausius Duhem) :
φ̇ = σij ij − ρψ̇ = Y Ḋ ≥ 0 (1.10)
Comme Y ≥ 0 ; ( D ≥ 0 ) : Cela signifie que l’endommagement est toujours crois-

sant. Une fonction critère est maintenant nécessaire pour définir la loi d’évolution de
l’endommagement. [8] la définit à partir de la notion de déformation équivalente :
f = − k(D) (1.11)
p
Avec ε̄ = Σε1 (εj )2 la déformation équivalente calculée à partir des déformations prin-
cipales positives εi (< εi >= εi siεi > 0, sinon < εi >= 0siεi 0).
K(D) est la valeur maximale atteinte au cours de l’histoire de chargement, égale initiale-
ment à D0 le seuil d’endommagement (obtenu à la fin de la phase élastique du matériau).
Une fois que ce seuil initial est atteint, l’endommagement est activé, selon le modèle,
l’endommagement est développé seulement en présence de déformations positives dues
aux mécanismes de traction. Comme le comportement du béton est dissymétrique entre
la traction et la compression, les lois d’évolution sont également définies séparément mais
sous une même forme :
εD0 (1 − Ai ) Ai
Dj ε̄ = 1 − − Bi (1.12)
ε̄ e (ε − εD0 )
Avec Di = Dt pour la partie d’endommagement en traction et Di =Dc pour la partie

d’endommagement en compression. Les paramètres Ai et Bi sont obtenus à partir des
tests uni axiaux en compression (Ac , Bc ) et en traction (At , Bt ), εD0 est le seuil d’initiation
de l’endommagement.
L’endommagement total est la combinaison de ces deux endommagements avec la prise
en compte d’un coefficient de distribution α de chaque partie (traction et compression),
calculé en fonction de l’état de déformation du matériau :
D = αtβ Dt + (1 − αt )β Dc (1.13)
Où αt proportionnel aux déformations positives, correspond à la traction,β est un pa-

ramètre désignant l’influence de cisaillement sur l’endommagement final. Un des défauts

du modèle de Mazars est qu’utiliser une seule variable scalaire D pour décrire l’état d’en-
dommagement du béton, même si cette variable comporte une part due à la traction et
une autre part due à la compression. Le caractère unilatéral de la structure (béton) se
traduit notamment par une restauration de la raideur due à la refermeture des fissures
lors du passage de l’état de traction à l’état de compression, selon le trajet de charge-
ment cela nécessite d’augmenter le nombre de variables où de mieux prendre en compte les
différents états possibles. En 1990 par Merabet cet aspect est bien modélisé dans le modèle
incrémenta proposé, ou le modèle d’endommagement de( Frémond et Nedjar, 1995).
Pour enrichir le modèle de Mazars, deux variables représentant des endommagements en
compression Dc ,comp et des endommagements en traction Dt ,tend sont utilisées (Mazars,
et al., 1990) et ( La Borderie, 1991).
Au départ de même potentiel d’énergie libre, mais le modèle d’endommagement proposé
par (Cervera et al., 1996) utilise un critère en contraintes pour définir l’évolution des
variables d’endommagement.
L’équation constitutive du modèle est exprimée en fonction de la contrainte effective σ̄
(sans endommagement) :
σ̄ = D0 : ε
Où D0 est un tenseur élastique linéaire d’ordre 4, ε est un tenseur des déformations. Dans
ce modèle, deux parties de compression et traction sont dissociées à l’aide des équations :
(Pour la partie en traction)
O
σ̄ + = <σ̄> = Σ31 <σ̄i >Pi Pi (1.14)
(Pour la partie en compression)

O
σ̄ = >σ̄< = Σ31 >σ̄i <Pi Pi (1.15)
Où σi représente la ième contrainte principale du tenseur σ̄ , Pi N représente le vecteur

ème
unité associé à la direction i principale respective et le symbole désigne le produit
tensoriel. Les symboles <.> sont les crochets de Macaulay (partie positive) et les symboles
>.< sont définies tels que<X>> + >X< =X. Deux fonctions critères différentes pour la
compression et la traction sont définies dans l’espace des contraintes :
(Pour la partie en traction)
g + (τ̄ + , r+ ) = τ̄ + − r+ ≤ 0 (1.16)
(Pour la partie en compression)
g − (τ̄ − , r− ) = τ̄ − − r− ≤ 0 (1.17)
p
Où(τ̄ + = σ̄ + : D0−1 : σ̄ +
représente
q√la norme de la contrainte effective équivalente en traction, et
−
(τ̄ − = 3K(σ̄ −oct + τ̄oct )
représente la norme de la contrainte effective équivalente en compression.
− −
σ̄oct et τ̄oct sont respectivement des contraintes normale et de cisaillement octahédrale. K
est une propriété du matériau qui dépend du rapport entre les résistances à la compression

uniaxiale et biaxiale du béton, de sorte que

√
K= 2(β − 1)/(2β − 1) (1.18)
(Les valeurs typiques pour le béton sont β = 1.16, K=0.171).
Les variables r+ et r− sont des seuils d’endommagement, introduits pour contrôler la taille
d’expansion des surfaces endommagées. Les premières valeurs des ri sont les résistances en
traction uni axiale r0+ et en compression uni axiale r0− . Les endommagements di augmentent
une fois que les contraintes équivalentes ont atteint les seuils ri .
Les lois d’évolutions de ces variables d’endommagement sont définies telles que :
Partie en traction +
+(1− r + )
+ r0+ A r0
d = 1 − +e (1.19)
r
Partie en compression
−
−(1− r − )
r− r0
d− = 1 − 0− (1 − A− ) − AeB (1.20)
r
Oùd+ etd− sont respectivement les variables d’endommagements en compression et en

traction (0 ≤ di ≥ 1). Les paramètres A+ ,A− , B − sont déterminés à partir des tests
uniaxiaux en compression et en traction.
L’expression de la contrainte finale est maintenant réécrite avec l’introduction des variables
d’endommagement :
σ = (1 − d+ )σ̄ + + (1 − d− )σ̄ − (1.21)
Les quelques modèles d’endommagement isotropes cités ci-dessus donnent la philosophie

générale de leurs concepts. L’hypothèse d’endommagement isotrope pour les structures
béton armée ou non semble valable seulement sous certaines conditions de chargement
(i.e. chargement uniaxial).
En effet, la définition du caractère isotrope du matériau est valable si le comportement
du matériau ne dépend pas des directions de chargement, tandis que les microfissures et
macrofissures apparaissant dans le béton dépendent des directions de chargement ; elles
sont perpendiculaires aux extensionsprincipales. Certains auteurs ont donc introduit le
caractère anisotrope de l’endommagement dans des modèles de comportement du béton
(Dragon & Mroz, 1979), (Chaboche, 1981), (Ortiz, 1985), (Chaboche, et al., 1995), (Ga-
tuingt & Pijaudier-Cabot, 2002), (Chambart, 2009).
Pour cela, une variable d’endommagement tensorielle est utilisée.
1.6.2 Modèles de plasticité

Les modèles élastoplastiques ont été initialement développés pour modéliser les com-
portements des métaux (Hill, 1950), (Bai & Wierzbicki, 2008) (Svoboda, et al., 2010), et
aussi utilisés pour modéliser le comportement du béton pourtant fragile en utilisant des
lois incrémentales.
Avec la partition de la déformation du matériau en deux déformations élastique et plas-
tique, Un modèle de plasticité est fondé par la formule suivante :
dε = dεe + dεp (1.22)

où εe est le tenseur des déformation élastiques,

εp est le tenseur des déformation plastiques,
ε le tenseur des déformations totales
on obtient l’incrément de contrainte à l’aide de la loi de Hooke généralisée :
dσ = Cdεe = C(dε − dεp ) (1.23)
où C est un tenseur de rigidité élastique d’ordre 4 contenant les constantes d’élasticité du
matériau.
L’incrément de déformation plastique est obtenu par l’intermédiaire d’une fonction seuil
(ou critère) de plasticité ainsi qu’une loi d’écrouissage (isotrope, cinématique, ou mixte).Les
différents modèles de plasticité se distinguent par ces deux aspects.
Fonction seuil de plasticité

Pour la structure (béton), le comportement en compression est très différent de celui
en traction. Donc, pour simuler ce type de comportement, plusieurs fonctions critères
peuvent être utilisées dans un même modèle (Voyiadjis, et al., 2008). Dans la littérature,
les différentes fonctions critères proposées se différencient par le nombre de paramètres
qu’ils utilisent. Le critère le plus simple, à un paramètre lié au second invariant du tenseur
des contraintes J2, est par exemple celui de Von-Mises ou celui de Tresca, plus adapté
au comportement des métaux que du béton. Les critères à deux paramètres permettent
de mieux représenter le comportement des géomatériaux par leur dépendance à J2 et à
la pression hydrostatique, utiles pour reproduire le comportement triaxial (Drucker &
Prager, 1952). Mais ce deuxième type de critère ne permet pas encore de bien distinguer
le comportement en compression triaxiale avec le comportement en extension lorsque la
structure est soumise à un chargement multidirectionnel. Par conséquent, des critères à
trois paramètres ont été développés. Outre la dépendance à la pression hydrostatique, un
paramètre prenant en compte le trajet de chargement est nécessaire. Dans le cas du critère
proposé par (William & Warnke, 1974), le troisième invariant du tenseur de contrainte J3
est pris en compte par l’intermédiaire de l’angle de Lode . Ce type de critère a également
été proposé par certains autres auteurs pour le béton (Etse, 1992), (Menétrey & Willam,
1995), (Thabet & Haldane, 2000), (Raveendra, et al., 2006). Des critères à 4 paramètres
(Hsieh, et al., 1982), voire 5 existent aussi (Willam & and Warnke, 1975), dans le but
d’obtenir des fonctions critère à forme parabolique à très fortes pressions hydrostatiques.
Loi d’écrouissage σ Une fois que le matériau atteint son seuil limite d’élasticité, des
déformations plastiques irréversibles apparaissent. Il existe deux modèles principaux ca-
ractérisant le régime plastique : avec et sans écrouissage. - Sans écrouissage, la contrainte
atteint un niveau et reste constante quel que soit la déformation. - Avec écrouissage,
la contrainte évolue avec la déformation plastique (écrouissages isotrope ou cinématique
(Ragueneau & Gatuingt, 2003)). Le modèle sans écrouissage est simple et plutôt facile à
utiliser pour l’acier. Le modèle avec l’écrouissage est nécessaire dans le cas du béton. Il
permet de modéliser à l’échelle macroscopique l’augmentation ou la perte de résistance du
béton associé à l’ouverture ou la fermeture des microfissures lors du chargement. La perte
de la raideur du matériau n’est pas prise en compte dans ce type de modèle contrairement
aux modèles d’endommagement.

Figure 1.6 – (Comportement uniaxial du béton lors d’un cycle chargement –

déchargement a)( modèle d’endommagement, b) modèle de plasticité avec l’écrouissage
négatif
1.6.3 Modèles couplés élasto-endommageable-plastique

Le modèle de couplage a pour avantage de faire la description de l’évolution des
contraintes en prenant en compte à la fois les déformations plastiques irréversibles et
les pertes de raideurs du béton sous chargement axial. Les simulations de tests uni-
axiaux de traction et de compression avec un modèle couplé proposé par (Voyiadjis, et
al., 2008) illustrent bien cette idée ([Figure 1-7]). Le modèle utilise à la fois le tenseur de
déformation plastique et le tenseur d’endommagement lié aux mécanismes des microfis-
sures pour décrire l’état déformé du matériau. Les fonctions critères utilisées sont d’une
part celle du modèle d’endommagement et d’autre part celle du modèle de plasticité. Le
modèle couplé développé par (Jason, et al., 2006) donne aussi des bons résultats en se
basant sur les mêmes principes. Sous un chargement caractérisé par une forte compression
Figure 1.7 – (Comportement du béton sous a)( la traction uniaxiale et b) la compression

uniaxiale (Voyiadjis, et al., 2008)
hydrostatique, la raideur du béton augmente en raison de la fermeture de la porosité. Le

modèle couplé proposé par (Burlion, et al., 1998) prend en compte ce phénomène.
C’est un couplage entre un modèle d’endommagement de type Mazars et un modèle
élastoplastique de type Gurson. Dans ce modèle, le tenseur de déformation plastique
symbole décrivant la fermeture de la porosité et la variable scalaire d’endommagement
D sont utilisés pour décrire l’état de déformation du matériau.Comme (Voyiadjis, et al.,
2008), ces auteurs utilisent deux fonctions critères pour décrire la rupture.

Selon le modèle, l’endommagement de compression est lié au premier invariant de la

déformation plastique irréversible (déformation volumique plastique). A l’inverse, la déformation
plastique n’est pas activée en tension. Par conséquent, sous un chargement de traction,
seul l’endommagement modélisant l’apparition des micro-fissures est activé.
Figure 1.8 – Comparaison entre résultats expérimentaux et de simulation pour un mor-

tier sous compression (Burlion, et al., 1998)
Même si tous ces modèles couplés semblent disposer des ingrédients permettant de
décrire les phénomènes irréversibles sous chargements complexes, une difficulté demeure
à l’intersection des surfaces seuil, une solution pouvant ne pas être régulière (Hansen &
Willam, 2001). Ce dernier problème a été étudié par (Thabet & Haldane, 2000).
1.7 Force du projectile

a) le cas des chocs mou :
dans le cas des chocs mous, l’Eurocode NF EN 1991-1-7 de février 2007 fournit quant à
lui un critère de résistance à l’impact lorsque la condition suivante est remplie :
1
F0 Y0 ≥ Mf V02 (1.24)
2
avec cette condition on peut tirer la formule de la force d’un projectile dans le cas d’un
impact mou.
1Mf V02
F0 ≥ ; (1.25)
2Y0
Mf :masse de l’objet impactant ;
F0 : est la résistance plastique de la structure, c’est à dire la valeur limite de la force
statique d’impact ;
Y0 : la capacité de déformation de la structure, c’est à dire le déplacement du point
d’impact que la structure peut subir.
La vitesse maximum que la cible peut endurer doit donc vérifier :
s
2F0 Y0
V0 ≤ (1.26)
Mf
Dans ce cas l’énergie cinétique est totalement transformée en déformation plastique.

b) le cas des chocs durs :

Dans le cas des chocs durs, l’Eurocode NF EN 1991-1-7 de février 2007 propose les
équations suivantes : p
F = V0 kMf (1.27)
avec
k : rigidité équivalente de l’objet
Mf : masse du projectile
V0 : vitesse de l’objet au moment du choc
Avec : r
Mf
F ∆t=Mf V0 et ∆t=
k
et
EA
k= et Mf = ρc EA
L
L étant la longueur de l’objet impactant, ρc sa masse volumique, A l’aire de la section
transversale et E le module d’élasticité. ∆t représente la durée de l’impact.
En introduisant un impacte dur localisé basé sur la force de pénétration, la demarche
adoptée consiste en la résolution de l’equation de mouvement d’un projectile rigide alors
gouvernée par la seconde loi de newton :
dV
Mf = FR (1.28)
dt
avec les conditions initales suivantes :
V = V0
X = 0 si t = 0
dxp
oùV = xp et v sont respectivement l’épaisseur de pénétration et V la vtesse du pro-
dt
jectile.
Mf :masse du projectile
FR :la force de resistance à la pénétration encore appelée force due à l’interaction projec-
tile/matériau cible durant l’impact.
Elle est généralement donnée sous la forme d’un polynome du second ordre en fonction
de la vitesse du projectile :
FR = FR (V ) = α + βV + γV 2 (1.29)
Elle peut aussi s’exprimer comme une fonction polynominale à deux variables :
x
p
FR = FR (xp ; V ) = g f (V ) (1.30)
d
Le mouvement dynamique du projectile rigide pendant le phénomène d’impact et de
pénétration peut alors être déterminé analytiquement permettant d’estimer l’épaisseur de
pénétration du projectile. Ces modèles peuvent également inclure la forme de la tête du
projectile.
Poncelet ou encore Wen ont proposé les expressions suivantes :
FR = A(a + bV 2 )(P oncelet) (1.31)

p
FR = A(αfc + β ρfc V )(W en) (1.32)
Dans ces formules,A représente la section du projectile, les coefficients (a, b, α, β) sont des

constantes à déterminer liées d’une part à la géométrie du projectile et, d’autre part, aux
propriétés mécaniques de la cible. Ces constantes sont, la plupart du temps, déterminées
expérimentalement.
Un des modèles analytiques force/pénétration les plus usités pour des cibles en béton est
basé sur les modèles de (Forrestal 1994) amélioré par la suite par (Li et Chen 2003) pour
des projectiles de formes diverses. Ce modèle force-pénétration est basé sur la théorie
d’expansion des cavités et décompose le phénomène de pénétration en deux temps.
Dans un premier temps, il considère une phase de pénétration de type cratère où FR
s’exprime comme :
xp πd N ∗ ρ V + Sfc
FR = C.xp pour <k avec C = c 0
d 4k (1 + πkd3 N ∗ ρ (1.33)
3Mf c
Dans un second temps, il considère une zone de tunnel où :

2
d xp
FR = π (S.fc + N ∗ .ρc .V0 ) pour ≥ k (1.34)
2 d
xp :profondeur de pénétration
S :facteur de pénétrabilité de la cible S = 72.0fc−0.5
d : diamètre du projectile
Mf :masse du projectile
fc :résistance à la compression du béton
N ∗ :facteur de forme de missile: 0 ≤ N ∗ ≤ 1
les facteurs de forme du missile :

1 1
1. N ∗ = − Projectile à tête ogivale,
3Ψ 24Ψ 2
R
avec Ψ = R : Rayon de l’ogive ; 0 ≤ N ∗ ≤ 0.5
d
1
2. N ∗ = projectile à tête conique,
1 + 4Ψ 2
h
avec Ψ = h : longueur du nez du projectile ; 0 ≤ N ∗ ≤ 1
d
1
3. N ∗ =1- 2 projectile à tête sphérique,
8Ψ
Rs
avec Ψ = Rs :Rayon de tête sphérique ; 0.5 ≤ N ∗ ≤ 1
d
4. Projectile à tête plate N ∗ =1 ; Projectile à tête hémisphérique N ∗ =0.5
l xp xp
k=0.707+( ) avec l :la longueur du projectile pour < 5, k = 2; pour ≥ 5
d d d
La résolution de l’équation de mouvement conduit alors à une épaisseur donnée par :
xp
Pour ≥ 0.5
d
s
kπ

xp 1 + 4N 4k
= I
I
d 1+ N π

xp
pour <k
d
!
xp 2 1 + NI
= N ln 1 + kπ
+k
d pi 1 + 4N
xp
Pour ≥k
d
xp
Pour ≤ 0.5
d !1.395
4k
xp π
I
= 1.628 I
d 1+ N
Avec
Mf V02

1 1 Mf
I= etN = ∗
S d3 f c N ρ c d3
L’ensemble de ces formules a été validé par un large panel d’essais expérimentaux pour
des projectiles de formes diverses et différentes vitesses d’impact.
1.8 Calcul de la résistance des structures (Béton)aux

impacts de projectiles
1.8.1 Généralité
De nombreuses recherches ont été entreprises pour développer des formules destinées à
décrire simplement les effets. La plupart se concentrent principalement sur l’établissement
de formules déterminant les performances en pénétration, en perforation ou encore en
écaillage. Ces notions sont en effet des facteurs clés en termes de dimensionnement de
protection par exemple. Les lois dites de prédiction de ces différents phénomènes se
répartissent en trois catégories : les corrélations empiriques, les méthodes analytiques,
les méthodes numériques.
1.8.2 Corrélations Empiriques

De nombreux modèles empiriques ont été développés par traduire les effets locaux d’un
impact de projectiles non déformables (durs) sur une cible que celle-ci soit en béton, en
béton armé ou en acier. Ces modèles permettent d’évaluer les phénomènes d’endomma-
gement de pénétration, de perforation ou encore d’écaillage (dans le cas d’une cible en
béton) à travers de différentes grandeurs géométriques :
La profondeur de pénétration Xp :
distance de pénétration du projectile dans la cible ;
L’épaisseur limite de perforation ou épaisseur de juste perforation ep :
épaisseur minimale nécessaire à la cible pour éviter la perforation ; Cette quantité est
parfois caractérisée par :
— La vitesse critique de perforation ou la vitesse de juste perforation Vp qui
exprime la vitesse minimale nécessaire permettant de perforer la cible ; Plusieurs
définitions existent. Backman et Goldsmith (1978) la définissent comme étant la
moyenne entre la vitesse la plus haute engendrant une pénétration partielle et la
plus basse vitesse engendrant une pénétration complète (perforation).

— ou encore l’énergie cinétique critique de perforation du projectile nécessaire afin de

perforer la cible
L’épaisseur limite d’écaillage es :
épaisseur minimale de la cible pour éviter la formation de l’écaillage sur la face arrière.
1.8.3 Formule génériques (Béton/Acier) pour impact localisé

1.8.3.1 Impact à vitesse supersonique
Un savant américain appelé Thor propose des équations empiriques issues d’un très
grand nombre d’essais de perforation de cibles métalliques et non métalliques par des
projectiles en acier. Ces équations permettent notamment d’estimer la vitesse de juste
perforation, la vitesse résiduelle du projectile en cas de perforation.
VP = 0, 3 × 10c1 (6, 1 × 104 H0 A)c2 (1, 55 × 104 Mf )c3 (sec(αi ))c4 (1.35)
5
Vr = V0 (1 − 10c (6, 1 × 104 H0 A)c6 (1, 55 × 104 Mf )c1 (sec(αi ))c8 (3, 3V0 )c9−1 ) (1.36)
L’épaisseur de juste perforation ep (en m) d’une cible en acier par un projectile en acier
massif est quant à elle donnée par TNO :
m1,063 × V01,103
ep = k(cos(αi )1,420 ) ; aveck = 1, 091 ∗ 10−7 kg −1,063 ∗ m1,897 ∗ S 1,103 (1.37)
A
V0 = vitesse d’impact (m/s)
Vr = vitesse résiduelle (m/s)
Vp = vitesse limite de perforation (m/s)
H0 = épaisseur de l’écran cible (m)
m = masse du projectile (kg)
A = surface d’impact moyenne ou surface exposée du projectile lors de l’impact (m2 ).
Lorsque la surface exposée du projectile lors de l’impact n’est pas connu, le TNO propose
de retenir pour A la valeur suivante correspondant à la surface d’un fragment standard
résultant de la fragmentation d’une bombe.
2/3
π m
A= (1.38)
4 0, 645ρacier
Avec ρacier :la masse volumique de l’acier
αi = angle de pénétration (par rapport à la normal de l’écran) (degré)
1
sec(αi ) = sécante de l’angle αi =
cos(αi )
C1 à C2 = coefficients fonction de la nature de la cible donnée.
Ces expressions sont valables pour des impacts à vitesse supersonique. Leur domaine de
validité est donné dans le tableau ci dessous ; il dépend de la nature du matériau consti-
tuant la cible, de son épaisseur, de la masse et la vitesse du projectile ainsi que de l’angle
d’impact du projectile sur la cible.
1.8.3.2 Impact à vitesse subsonique

Le High Pressure Safety code (1985, B.G. COX et G. SAVILLE) (Micado 2002)
propose une formule empirique permettant de calculer un ordre de grandeur des épaisseurs

de pénétration dans l’acier allié et le béton armé pour des gammes de vitesse inférieure à
1000 m/s.
Pour des fragments de masse supérieure à 1 kg et dont le rapport longueur/diamètre est
supérieur à 1, l’épaisseur de pénétration xp (en m) est donnée par :
CM
xp = log(1 + 5.10−5 V02 ) (1.39)
A
avec M : masse du projectile (en kg)
V0 : vitesse du projectile (en m/s)
A : section apparente du projectile lors de l’impact (en m2 )
C : paramètre donné dans le tableau suivant :
Table 1.1 – Paramètre de perforation selon la nature du matériau

Nature du matériau C
Pierre calcaire 7×10−4
Béton non armé- charge de rupture à la compression15MN/m² 10×10−4
Béton armé- charge de rupture à la compression22MN/m² 6×10−4
Béton armé -charge de rupture à la compression44MN/m² 3.5×10−4
Mur de pierres maçonnées 14×10−4
Briquetage 25×10−4
Acier doux 0.5×10−4
Acier allié 0.3×10−4
[Micado 2002]relie les différents types de projectiles, satisfaisant aux conditions précédentes,à
des épaisseurs de perforation.
Table 1.2 – Épaisseurs de perforation en fonction des caractéristiques des projectiles

Dénomination M (kg) V (m/s) A (m2 ) Ratio K acier en (mm) béton en (mm)
Vanne 35 35 0.06 20.5×10−3 0.45 5.3
−3
Vanne 35 100 0.06 58.3×10 3.1 36
−3
Fragment de réacteur 50 150 0.02 375×10 25 290
−3
couvercle de clapet 40 230 0.02 460×10 34 390
couvercle de cylindre 50 150 0.02 675×10−3 44 500
De la lecture du tableau on peut retenir qu’il existe en facteur multiplicatif de l’ordre

de 10 entre le béton et l’acier, en ce qui concerne l’épaisseur du matériau pour laquelle
on observerait une perforation.
Pour de petits fragments (m < 1kg etL/D ∼ = 1),la profondeur de pénétration P (en
m)peut être calculée par la formule :
xp = KM n1 V0n2 (1.40)
avec M = masse de fragment (en kg)

V0 =vitesse du fragment (en m/s)
K,n1 et n2 : constantes données ci-dessous

Table 1.3 – Paramètre de perforation selon la nature du matériau

Matériaux K n1 n2
2 −6
Béton à 35MN/m 18×10 0.4 1.5
−6
Mur de briques 23×10 0.4 1.5
Acier doux 6×10−6 0.33 1
Ces relations ont été établies dans le cadre du dimensionnement de protections d’équipements
et sont corrélées à diverses approches parmi lesquelles figure celle de Petry. Contrairement
aux équations de Thor, elles sont généralement utilisées pour des vitesses d’impact sub-
sonique. Elles sont de fait plus adaptées aux problématiques industrielles.
1.9 Conclusion
En génie civil les structures en béton armé sont les ouvrages les plus fréquents, ce-
pendant la connaissance de son comportement en dynamique rapide reste restreint. Pour
l’étude de son comportement il faut d’abord déterminer l’impact de projectile que la struc-
ture est soumise, vue la complexité des effets d’un tel impact, de ce type de sollicitation
amène à considérer deux cas limites, impact mou et impact dur. Koechlin et Potapov ont
donné une formule pour déterminer l’impact
σp ρp ν02
+ =1
σc σc
Selon leur critère si le terme gauche est inférieure à 1, le cas est qualifier d’impact mou et
impact dur s’il est supérieure à 1. Une fois le type d’impact est déterminé les modes de
rupture correspondant de la structure peuvent être prédits pour le dimensionnement de
la structure. Il existe plusieurs modèles de calculs en dynamique rapide dans ce présent
document les méthodes sont la méthode des éléments finis et la méthode des éléments dis-
crets. Pour calcul en dynamique des structures on étudie le comportement de la structure
sous impact. Avec l’hétérogénéité du béton le comportement est complexe et fortement
dépendant des conditions de chargement, il n’existe pas un modèle universel qui peut re-
produire tous les phénomènes physiques. Les modèles sont : le modèle d’endommagement
de Mazars, le modèle de plasticité et le modèle couplé élasto-endommageable-plastique.
Un point clé révélé par cette étude bibliographique concerne la nécessité de connaı̂tre
les modelés de calcul sous impact ainsi que son comportement dynamique. Différents
modèles de comportement du béton ont été proposés.Il semble qu’aucun modèle ne soit
capable de reproduire tous les phénomènes évoqués.
Les caractéristiques des projectiles sont considérées comme étant des variables aléatoires.
Les auteurs proposent un modèle permettant de d´écrire le mouvement du projectile dans
les deux cas. La résolution des diverses équations, renseigne sur la possibilité d’impact ou
non entre projectile et cible. Enfin, les formules ont été proposés, afin d’estimer la profon-
deur de pénétration de la cible impactée, puis parvenir a‘ l’estimation du risque d’accident.
Cette m´méthodologie permet de dégager une vue globale de l’effet Domino. Les résultats
obtenus soulignent l’importance de prédire cet effet dans le domaine industriel.

Chapitre 2
Simulation Numérique
2.1 Introduction
Un site industriel abrite diverses installations. Si une installation vient à exploser, les
fragments générés et projetés aux alentours peuvent impacter des installations voisines
tels que des canalisations, des réservoirs atmosphériques ou sous pression ou encore des
structures porteuses, ces installations peuvent aussi être percuter par une voiture. Ils
peuvent les pénétrer partiellement, les perforer entièrement ou écailler l’interface.
L’objectif de ce chapitre est de faire une simulation numérique sur une paroi porteuse
d’un site industriel suite à une collision de projectile.
Pour ce faire on procède comme suite :

Mise en œuvre des simulations : Dans cette partie nous allons présenter les
caractéristiques du projectile et la paroi, Déterminer les principaux paramètres de
calcul du projectile (Vitesse et force)
Analyse des résultats :Dans cette seconde partie nous allons donner les effets
du projectile sur la paroi (Force d’interaction et limite d’écaillage sur la paroi).
2.2 Mise en œuvre des simulations

La mise en œuvre est une séquence de calculs effectués sur un ordinateur et reprodui-
sant un phénomène physique. Elle aboutit à la description du résultat de ce phénomène,
comme s’il s’était réellement déroulé. Cette représentation peut être des données d’image
ou même un film vidéo.
Un simulateur peut réagir à des modifications de paramètres et modifier ses résultats en
conséquence. Un simulateur d’impact par exemple, modifie la trajectoire calculée du pro-
jectile en fonction des commandes transmises par l’utilisateur. Une simulation numérique
peut représenter des phénomènes physiques complexes dont la description repose sur un
modèle mathématique comportant des équations aux dérivées partielles.
Bien qu’une pratique courante consiste à représenter une simulation numérique de colli-
sion un projectile contre une paroi.
35
CHAPITRE 2. SIMULATION NUMÉRIQUE
2.2.1 Présentation de l’outil de simulation

Un outil de simulation est un logiciel de simulation numérique basé sur la méthode des
éléments finis. Ce logiciel permet de simuler de nombreuses physiques et applications en
ingénierie, et tout particulièrement les phénomènes couplés ou simulation multi-physiques.
Les logiciels de simulation numérique sont utilisés pour :
Prévoir l’état final d’un système connaissant son état initial (problème direct) ;
Déterminer les paramètres d’un système connaissant un ou plusieurs couples (état
initial - état final) (problème inverse) ;
Préparer des opérateurs à des conditions plus ou moins rares dans leur interaction
avec un système complexe (simulation d’entraı̂nement).
De nombreux logiciels de simulation commerciaux existent sur le marché. Dans le
présent travail, le choix s’est porté sur COMSOL Multiphycs, car il permet de simuler la
réponse mécanique des éléments de structure et de déterminer la déformation de éléments
soumis. Les lois qui gouvernent ces phénomènes sont celles de la thermique et de la
mécanique.
L’objectif principal de la simulation est de modéliser l’impact d’un projectile sur une
paroi porteuse et voir comment les fissures se propagent dans cette paroi.
2.2.2 Principaux paramètre de calcul d’impact

2.2.2.1 généralité
Le modèle 2D et 3D sont a priori plus proche de la réalité et sont opérationnel. Tou-
tefois, les calculs 3D nécessitent beaucoup de temps. Ce modèle permet également de
comprendre l’origine et la nature du comportement non linéaire du béton.
La modélisation est la conception d’un modèle selon son objectif et les moyens utilisés.
Nous avons choisi, dans le cadre de ce travail, le modèle élastique-endommageable de
Mazars (Mazars, 1984). Les étapes de la simulation consistent à faire :
Modéliser la paroi ;
Définir les paramètre de calcul ;
Appliquer la force d’interaction sur la paroi ;
Interpréter les contraintes de déformation sur la paroi
2.2.2.2 Caractéristique du projectile

a. On va réaliser la simulation d’impact de projectile sur une paroi dans le but d’obtenir
la mode de propagation de fissure dans le béton sous un chargement sévère (chargement
dynamique). Le projectile utilisé pour le comportement est le projectile cylindrique (masse
100 kg),la vitesse d’impact est de 8,43 m/s, le diamètre du projectile et on choisira la durée
d’impacte (temps) selon[9] T= 20ms.
b. D’après la seconde loi de Newton nous allons établir une relation entre la quantité
de mouvement P et la force F pour déterminer la force d’interaction du projectile qui
réagira sur la paroi comme une charge surfacique lors de la simulation.
FORMULE :
∆m × V
F = (2.1)
∆T
d’où
∆P
F = (2.2)
∆T

avec P = m × V
P :la quantité de mouvement
m :la masse du projectile
∆T :la durée d’impact
V :la vitesse du projectile
F :la force d’interaction sur la paroi impacté
2.2.2.3 Caractéristiques de l’objet impacté (paroi)
L’impact d’un projectile sur une structure en béton armé peut générer simultanément
des dommages locaux et globaux. Un modèle d’éléments finis 2D où 3D éléments solides
en béton armé est le plus approprié pour reproduire la fissure et la pointe de déformation
complexe généré lors de l’impact sur la cible. De ce fait la cible est un matériau constituer
de béton plus de l’acier.
De sa capabilité de résistance à la compression, nous pouvons tirer et classer les bétons :
Table 2.1 – classe du béton

classe résistance caractéristique (28jours) MPa
Béton ordinaire 16-40
Béton Haute Performances 45-60
Béton Très Hautes Performances 65-100
Béton Ultra Hautes Performances ≥100
Dans le tableau si dessus La résistance en compression du béton est désignée par des
classes de résistances liées à la résistance caractéristique mesurée sur cylindre ou sur cube,
conformément à l’ EN 206-1 . Pour ce présent mémoire la simulation se basera pour un
béton de haute performance C45
- Caractéristique géométrique de la paroi :
La hauteur de la paroi rectangulaire de béton armé composant la cible est de 3 m.
L’épaisseur de la paroi est de 40 cm et la largeur est de 150 cm. L’armature est com-
posée de deux couches d’acier (une de chaque côté de la paroi de béton) avec des barres
de 16 mm de diamètre. D’autres barres de 10 mm de diamètre relient chaque nœud de la
maille d’armature des couches face à face.
- Caractéristique du béton :
La résistance caractéristique à la compression à 28 jours fc28 = 45M P a
La résistance caractéristique à la traction à 28 jours ft28 = 3.3M P a
Le coefficient partiel de sécurité pour le béton γb = 1.5
- Caractéristique de l’acier :
La limite d’élasticité garantie fe = 400MPa
Le coefficient partiel de sécurité pour l’acier γs = 1.15
Remarque : Le comportement mécanique du béton présente des caractéristiques va-
riables de déformation sous diverses conditions et niveaux de charge. Une grande part de la
connaissance des mécanismes de son comportement résulte d’observations expérimentales.
L’accumulation des connaissances empiriques et 1’expérience acquise sont des éléments
très important à prendre en considération lors de la modélisation numérique du compor-
tement du béton.

A : Une portion de la structure B : Squelette de la structure
Figure 2.1 –
2.2.2.4 Définition des paramètres pour la loi [8]

Les paramètres de la loi[8] sont déterminés à l’aide des caractéristiques du matériaux.
En connaissant la valeur de résistance en compression on peut définir le module de Young
par la formule suivante : p
E = 11000 3 fc28 (2.3)
Le coefficient de poisson est le ratio négatif de la souche de transverse à la souche d’axial,
c’est à dire l’effet de Poisson est la circonstance quand un matériau comprimé à une di-
rection et il a la tendance à se développer dans les deux autres directions perpendiculaires
ou parallèles à la direction d’écoulement de chargement. La valeur de ce coefficient est
entre 0.10 et 0.20. Pour le présent travail nous allons choisir ν = 0.2.
β étant un facteur de paramètre désignant l’influence de cisaillement sa valeur est de 1.16
2.2.2.4.1 Détermination de la déformation équivalente ε̄
La formule de la déformation d’équivalente est donnée dans le chapitre 1 [1.6.1]
2.2.2.4.2 Paramètre d’évolution de l’endommagement en traction et en com-

pression
où Ac , Bc , At , et Bt sont des paramètres d’identification du matériau et peuvent être

déterminés à partir des formules.
en traction :
E
Bt = etAt = 0.70 (2.4)
ft28
en compression :
√
fc28 ν 2
1 − ε0 E
Bc = √ etAc = √ √ (2.5)
εc ν 2 εc ν 2%exp(Bc (ε0 − εc ν 2)) − ε0
ft28 √
Avec ε0 = ; ε0c = 0.620E − 3 3 fc28
E
Le tableau récapitulatif des paramètres de la loi [8] :

Table 2.2 – Valeurs des paramètres de la loi [8]

fc28 [MPa] ft28 [MPa] E[MPa] ε0c ν At Bt Ac Bc F [N]
−3
45 3.3 39126 2.21×10 0.2 0.7 11856 1.351 1603.2 42150
2.3 Analyse des résultats

2.3.1 Les effets du projectile sur la paroi
Les résultats numériques obtenus sont présentés sur les Figures 2.2. La localisation de
la déformation est bien illustrée à la Figure 2.2, où une fissure critique se développe le
long de l’axe x de la paroi.
En ce qui concerne la compression, on peut voir sur la Figure [2.2 a] l’apparition d’un
endommagement uniforme jusqu’à’à l’apparition de fissures. On voit nettement le compor-
tement fragile et on peut conclure que le modèle permet de reproduire le comportement
en compression. La valeur de la déformation en compression varie entre −2.23 à 1.98.
Lorsque le béton est endommagé sous une sollicitation extérieure puis déchargé, la fissure
qui apparait sur la structure ne se referment pas totalement. La rugosité de fissure formée
durant le chargement empêchent le retour à l’état initial.cet phénomène permet de nous
montrer la plasticité du matériaux.
a) Déformation équivalente b) Déformation plastique
Figure 2.2 – Zone de localisation des déformations
Figure 2.3 – Mazars : Comportement en traction-compression.
La Figure 2.3 donne l’allure de la loi de Mazars et pour un essai de traction et de

compression, dans se cas la déformation de compression correspondante à la contrainte
maximale (valeur absolue) est fixée aux alentours de 2.21%.

2.3.1.1 MAILLAGE
le maillage est défini dans le tableau suivant :
Table 2.3 – Définition du maillage

Description Valeur
Qualité minimale des éléments 0.007492
Qualité moyenne des éléments 0.253
Prisme 112
Hexaèdre 1190
Triangle 24
Quadrangle 687
Elément d’arêtes 760
Elément ponctuel 136
Figure 2.4 –
2.3.1.2 Chargement appliqué
Figure 2.5 – chargement surfacique
2.3.1.3 Fissures
Lorsque le béton est doté d’une résistance normale est comprimé,les fissures d’adhérence
entre la matrice de mortier et la granulat se propageront autour des granulats.

Les motifs et la cinématique de l’apparition des fissures sont similaires à ceux issus de
l’endommagement. Dans la zone de localisation à l’interface de la paroi, la valeur de l’ou-
verture des fissures de la zone impacté est 5.73 ×10−7 m et à l’arrière face de la zone
impacté la valeur de l’ouverture de fissure est 3×10−7 m .
Figure 2.6 – Propagation des fissures
Lorsque le temps progresse, la fissure s’ouvre mais aucune nouvelle fissure n’apparaı̂t.
Au moment du rebond on remarque une forte dégradation locale du béton sous l’impact.
Avec la vitesse d’impact de projectile utilisé on sent qu’il a un déplacement plus ou moins
grand.
Figure 2.7 – Amplitude du déplacement
En observant la surface d’interaction on peut dire que le projectile a endommagé la

structure. on voit clairement l’apparition de l’écaillage dans la zone impactée.
2.3.2 Conclusion
Ce dernier chapitre conclut donc ce travail par la validation de notre modèle numérique
par la simulation d’impact de projectile contre une paroi en béton armé avec le modèle
de Mazars. la simulation est fait avec un béton de haut performance, dont la classe de
résistance à la compression est 45MPa et la résistance à la traction est 3.3MPa.

Deuxième partie
DOSSIER PÉDAGOGIQUE
42
Chapitre 3
Rapport du stage pédagogique
3.1 INTRODUCTION
L’Ecole Normale d’Enseignement Technique et Professionnel (ENETP) est un établissement
d’enseignement supérieur public situé à la cité universitaire de Kabala (Mali), créé par
ordonnance N° 0032/P-RM du 10 août 2010 à Bamako et ratifié par la loi N° 042 du
15 juillet 2011. L’école vise à donner plus de visibilité à la formation des professeurs
d’enseignements technique et professionnel et à pallier au manque d’enseignants dans ce
domaine, elle regroupe quatre (4) filières de formations qui sont :
— Génie Civil (GC) ;
— Génie Industrie (GI) ;
— Bureautique Communication (BUCO) ;
— Comptabilité Gestion (COGE).
L’École Normale d’Enseignement Technique et Professionnel (ENETP) a pour mission à
savoir :
La formation initiale des professeurs d’enseignement technique et professionnel ;
La formation et le perfectionnement des personnels d’encadrement de l’enseigne-
ment technique ;
La formation poste universitaire en sciences de l’éducation ;
L’élaboration et la production des supports pédagogiques ;
La formation professionnelle continue de formateurs d’organismes de formation
publics et privés ;
La promotion et le développement de la recherche pédagogique et technologique ;
La recherche scientifique.
En fin de cycle l’école prévoit un stage pédagogique pratique pour les élèves professeurs
dans des établissements de l’enseignement technique et professionnel en leur permettant
d’appliquer les théories apprises en classe sur le terrain.
Le stage pédagogique est une période pendant laquelle les élèves professeurs de l’Ecole
Normale d’Enseignement Technique et Professionnel de Bamako sont mis à la disposition
des professionnels de l’enseignement technique et professionnel afin que ceux-ci les ac-
compagnent à la construction des compétences individuelles et collectives. C’est donc une
occasion pour nous, de rencontrer de nombreuses situations apprises en classe, de déceler
les problèmes qui peuvent survenir sur le terrain et d’essayer de résoudre ces problèmes
par l’expérience théorique acquis au cours de notre formation.
C’est dans ce but que j’ai effectué mon stage pédagogique au Centre de Formation en
Administration et aux Métiers (CFAM) sise à Sébéninikoro Derrière le Commissariat de
43
CHAPITRE 3. RAPPORT DU STAGE PÉDAGOGIQUE
police du 9e Arrondissement et près de l’Ecole Fondamentale publique de Kalabambougou

extension en Commune IV-Bamako.
3.2 Présentation de l’école d’accueil

3.2.1 Historique
Le Centre de Formation en Administration et aux Métiers (CFAM) est un établissement
d’enseignement technique et professionnel. Il est créé par l’ordonnance N° 2011 – 2082/MEALN-
SG du 1er juin 2011 et ouvert par l’ordonnance N°2011 – 3528/MEALN-SG du 1er Sep-
tembre 2011.
3.2.2 Situation géographique

Le Centre de Formation en Administration et aux Métiers (CFAM) est situé dans la
commune IV du district de Bamako à Sébénikoro Derrière le Commissariat de police du 9e
Arrondissement et près de l’Ecole Fondamentale publique de Kalabambougou extension.
L’école est bâtie sur une superficie de 400m2 dont elle est constitué de 6 classes et 2 salles
de pratiques plus un bloc qui contienne l’administration. C’est un centre de formation
d’enseignement technique et professionnel, il est placé sous l’autorité de l’académie d’en-
seignement de Bamako rive gauche.
Sa devise est : Discipline – Travail – Solidarité.
3.2.3 Organisation administrative (organigramme)
Figure 3.1 – Organigramme de l’établissement.
3.2.3.1 L’organe suprême de gestion scolaire

Dans cet établissement d’enseignement technique et professionnel, le premier respon-
sable est le Directeur général. Son adjoint est le Directeur des études, en cas d’absence
ou d’empêchement il remplace le directeur général, Sa mission est de se charger à l’orga-
nisation et du fonctionnement des cours. Le chef des travaux et le surveillant général sont
respectivement chargés de l’organisation et du fonctionnement des ateliers et du respect

de la discipline, placée sous l’autorité du directeur des études. L’Econome, placé sous
l’autorité du Directeur général, il est chargé de la gestion financière et matérielle.
3.2.3.2 Association des parents d’élèves

L’association des parents d’élèves doit promouvoir les intérêts de tous les élèves de
l’école en coopération avec le conseil d’administration, le directeur, les enseignants et les
élèves de l’école. Elle doit donc dans ce but conseiller le directeur ou le Conseil d’admi-
nistration sur tout problème concernant l’école.
3.2.3.3 La direction
Le Centre de formation en administration et aux métiers est administré par :
— Un directeur général, qui est le chef de l’établissement ;
— Un directeur des études, qui est chargé des études et de la discipline ;
— Un chef des travaux, responsable des ateliers ;
— Un surveillant général, chargé de la discipline ;
— Un économe responsable de la gestion financière et matérielle au sein de l’établissement.
3.2.3.4 Le conseil des maitres

Il réunit l’administration (Directeur Général, Directeur des études, Surveillant, Secrétaire,
chef des travaux et économe) et les maı̂tres affectés à l’école. Il est présidé par le Directeur
Général, car c’est lui qui assure la police des débats.
3.2.4 Organisation Pédagogique

3.2.4.1 Groupe pédagogique
Un groupe pédagogique est composé d’élèves qui suivent un programme d’études com-
mun. Les élèves sont regroupés en groupes pédagogiques sur la base du plus grand nombre
de cours communs. Un groupe pédagogique désigne la structure pédagogique dans laquelle
chaque élève est inscrit.
Vue les circonstances dans laquelle l’école se trouve des fois (grèves et autres problèmes),
nous les professeurs nous nous regroupons pour discuter de la meilleure stratégie à adapter
pour rehausser le niveau de nos élèves, enfin de compte pour avoir des bons résultats à la
fin de l’année scolaire.
3.2.4.2 Pédagogie et Didactique appliquées

Au CFAM il existe deux (2) approches pédagogiques qui sont : l’approche dite classique
et l’approche par compétence (APC).
3.2.4.3 Les filières de formation

Le centre de formation en administration et aux métiers qui est un lycée d’Enseigne-
ment Technique et Professionnel, offre une bonne qualité de formation en des filières qui
sont :
— Les filières industrielles (Le Dessin Bâtiment, l’électromécanique, l’électricité) ;
— Les filières tertiaires (Technique comptable, Aide-Comptable, Secrétariat de Direc-
tion.

3.2.4.4 Les emplois du temps

A chaque ouverture de classe, le directeur des études met les emplois du temps à la
disposition des enseignants et enseignantes qui interviennent au CFAM.
3.2.4.5 Les effectifs

Le centre de formation en administration et aux métiers compte deux cent quarante
élèves à toutes les filières confondues.
Table 3.1 – Récapitulatif des effectifs d’élèves par filières

CAP Industrie CAP Tertiaire BT Industrie BT Tertiaire
Garçon 40 65 22 40
Filles 4 38 6 25
totaux 44 103 28 65
3.2.4.6 Le fonctionnement des classes

Généralement au Centre de formation en administration et aux métiers, les cours se
déroulent du lundi au Samedi de 07h00 à 17h00 avec une recréation entre 9h et 10h.
Toutefois, les heures diffèrent selon la durée de certaines matières. De plus en plus l’aug-
mentation des effectifs rend la gestion des salles difficile.
3.3 Activités menées

3.3.1 Activités pédagogiques
Pendant mon stage, l’objectif visé était d’appliquer les connaissances théoriques ap-
prises à l’ENETP sur le terrain pour compléter ma formation. Après avoir accepté ma
demande de stage pédagogique on m’a confié à un professeur qui va m’encadrer Mr.
Sissoko Boubacar professeur de Dessin Bâtiment. Mon stage pédagogique au CFAM s’est
déroulé en différentes étapes. En une première étape, j’ai d’abord commencé à assister
mon encadreur Mr. Sissoko Boubacar dans la classe de 2ème année Dessin Bâtiment, en-
suite petit à petit, j’ai commencé à donner des cours dans la même classe.
En appliquant les judicieux conseils de mon encadreur, il me donna plus de responsabilité
au fur de l’année (corriger les feuilles des devoirs, surveillés les devoirs. . . ).
3.3.2 Activités parascolaires

Les activités parascolaires sont essentielles et occupent une place très importante dans
la vie quotidienne des élèves. Vues les occupations des élèves (les cours à l’école, les devoirs
à domicile, les travaux à la maison, etc.), il est important pour les élèves d’avoir des temps
de repos, mais aussi des moments de divertissement, c’est dans ce cadre que chaque année
le CFAM organise un tournoi de football entre les différentes filières et une nuit de soirée
traditionnelle pour la fin d’année.
Ainsi pour une meilleure cohésion et des moments de divertissement, les responsables
des différentes filières se sont concertés pour constitués des équipes de foot masculines et
féminines pour le tournois d’interclasses.

3.4 Difficultés rencontrées

Lors de mon stage pédagogique Les difficultés rencontrées sont nombreuses tant sur le
plan matériel que sur le plan pédagogique.
3.4.1 Matérielles
— Insuffisance de matériels didactiques ;
— Insuffisance de livre dans la bibliothèque ;
— Insuffisance de salle de classe.
3.4.2 Pédagogiques
Les connaissances théoriques et les méthodes d’enseignement que j’ai appris à l’ENETP
m’ont permis de dispenser les cours qu’on m’a confié sans avoir le stress devant les élèves.
Les difficultés rencontrées sont nombreuses, parmi lesquelles on peut citer :
— Les élèves ne font pas de révision avant de venir en classe ;
— La non maitrise du français par les élèves ;
— Les élèves ne font pas de travail individuel à la maison ;
— Les élèves ne s’exercent pas à la maison ;
— L’indiscipline des élèves ;
— La non implication de certains élèves à certaines disciplines ;
— A tout moment il peut effectuer un transfert.
3.4.3 Socioculturelles
Les difficultés rencontrées sont les suivantes :
— Le retard des élèves ;
— Le non suivi des élèves à la maison ;
— Le manque de respect envers les autres.
3.5 Solutions apportées

3.5.1 Initiatives personnelles et apport au d´développement de
l’école
Malgré le mauvais comportement et le faible niveau de certains élèves, j’impliquais
tous les élèves à mon cours. J’ai eu des difficultés pour les faires raisonner afin de trouver
de bonnes solutions pour améliorer le niveau. Avec le respect mutuel (professeur-élèves),
j’établissais un climat de confiance, de dialogue favorable à l’apprentissage de mes élèves.
Même après mes heures de cours j’étais toujours ouvert avec mes élèves à des explications
concernant les cours entrant dans leur formation.
3.5.2 Autres
Lors de mon stage pédagogique, j’ai insinué à mon encadreur Mr. Sissoko si toute fois
après les cours, il pouvait faire une liste de présence de la classe pour la surveillance. Lors

du choix des sujets, j’aidais le professeur Sissoko en lui donnant mes avis et en l’aidant à
faire comprendre les élèves avec les dessins à l’aide du logiciel AutoCad.
3.6 Les défis et perspectives

Durant cette période de stage, mon but était de vérifier ce que j’ai appris à l’Ecole
Normale Enseignement Technique et Professionnel, mon savoir-faire et mon savoir être
aux élèves et qui a été à la hauteur de souhait.
3.7 Leçons tirées du stage

3.7.1 Les intérêts
Ce stage a été très bénéfique pour moi, car il m’a permis :
De lier la théorie à la pratique ;
D’améliorer la qualité de mon travail pour la préparation de mes leçons ;
De m’adapter aux réalités du terrain et ainsi avoir des expériences supplémentaires.
3.7.2 Les motivations

Enseigner est une fonction noble, qui est une de mes raisons de motivation car le
développement d’un pays est en grande partie lié à son éducation. L’amour que je porte
à ce métier est immense vue mes oncles enseignants et mes collègues.
3.7.3 Points de satisfactions

Lors de ce stage, j’ai eu à faire des remarques sur certains points qui m’ont satisfait,
parmi lesquels nous avons :
Le respect mutuel entre les élèves, les enseignants et les personnels de l’adminis-
tration ;
L’égalités est bien respecté entre garçons et filles dans les responsabilités et les
tâches qui leur sont confiées ;
L’échange d’idée entre les enseignants de l’établissement et les stagiaires.
3.8 Les recommandations

3.8.1 Points de convergence entre le stage et la formation à
l’ENETP
Lors de mon stage, j’ai remarqué que la réalité du terrain s’adapte en grande partie
avec la formation reçue à l’ENETP.
3.8.2 Points de divergences entre le stage et la formation à

l’ENETP
Pendant mon stage, j’ai remarqué qu’il y a une nuance entre la formation et le stage à
l’ENETP, c’est à dire la théorie et la pratique était parfois très difficile et parfois quelques

réalités du terrain étaient différentes de celles apprise en classe. La pratique est plus
formative que la théorie, il faut que les autorités de l’état aient le courage et la volonté de
revoir le rapport entre la théorie et la pratique. La théorie est facile mais la pratique est
difficile. Dans les écoles techniques et professionnelles la préparation de la leçon est aussi
une dimension très importante.
3.8.3 Besoins en formation pour rendre mon travail plus efficace

A travers les difficultés rencontrées lors de mon stage, pour améliorer mon travail il
faut :
Doter les écoles techniques et professionnelles de matériels didactiques ;
De la documentation des bibliothèques scolaires ;
De la formation continue.
3.8.4 Autres remarques et suggestions

Il faut :
— L’augmentation de salles de classe compte tenu du manque de salles de classe ;
— L’approvisionnement de l’établissement en matériels didactiques et techniques suf-
fisantes ;
— Accentuer le suivi des élèves en difficultés ;
— Favoriser l’apprentissage des nouvelles technologies en rapport avec la formation
des élèves ;
— L’installation d’une cantine scolaire ;
— Renforcer le cadre de concertation pour faciliter l’intégration des nouveaux au sein
du corps professoral de l’établissement.

3.9 La leçon modèle

3.9.1 Partie administrative
Etablissement : Centre de formation en administrations et aux métiers (CFAM)
Nom de l’enseignant : DOUMBIA
Prénoms : N’faly
Spécialité : DESSIN BÂTIMENT
Année Scolaire : 2018 - 2019
Date : 15 / 08 / 2019
Matière : LA TECHNOLOGIE GéNéRALE (TG)
Leçon n° : 4
Durée prévue : 1H45 MN Classe ou niveau : 3 ème Année Dessin Bâtiment
TITRE : L’ASSAINISSEMENT
OBJECTIF PEDAGOGIQUE : A la fin de la leçon l’élève doit être capable d’identifier
les différentes types d’assainissement
DOMAINE TAXONOMIQUE : Cognitif
NIVEAU TAXONOMIQUE : Connaissance
PRE-REQUIS : Les ouvrages d’assainissement (les puisards, les fosses septiques, les
latrines, les caniveaux etc. . . )
SUPPORTS : Tableau, les craies, la règle
DOCUMENTS ELEVES : Prise de note par les élèves
BIBLIOGRAPHIES :
internet ; V. SAVARY, TECHNOLOGIE DE MAÇONNERIE, quatrième Edition, 290
pages.
PREAMBULE MOTIVANT
Questionnaire :
1. Les eaux de toilettes sont évacuées vert où ?
2. Après la pluie ou la lessive avec quoi l’eau est jetée dans le fleuve (milieu naturel) ?

DEROULEMENT DE LA SEQUENCE
Etapes Plan de la leçon

Micro-objectif 1 : I. Définition(l’assainissement)
— a. Eau pluviale ou Eau de ruissellement
— b. Eaux usées
— c. Eaux usées industrielles
Micro-objectif 2 : II. But de l’assainissement

Micro-objectif 3 : III. Les différents types d’assainissement
— a. L’assainissement individuel
— b. L’assainissement collectif
Micro-objectif 1 : Production élève :

Définir L’assainisse- — l’assainissement c’est de transporter l’eau vert le rejet.
ment — l’assainissement est une opération d’évacuation des
eaux usés dans le fleuve.
Durée : 15 mn
Point clés : Stratégie :
Après avoir utilisé la toilette ou de faire les tâches ménagères
— Opérations dans vos familles vous verrez que l’eau et les déchets sont di-
— Traitement rigés dans un lieu qui s’appelle la fosse septique et si la fosse est
— Milieu remplie il y a une partie qui filtre l’eau pour que les déchets ne
récepteur passent pas pour aller dans le puisard et ce système est appelé
système de drainage.
Transition :
— 1 Pourquoi on fait tous ceci chez nous ?
— 2 Est ce qu’on peut laisser ces eaux stagner devant nos portes ?


But de l’assainisse- — De rendre l’environnement propre
ment — De maintenir la santé
— D’éviter l’inondation dans nos quartiers
Durée : 5 mn
C’est en faisant assainir nos familles et le quartier qu’on évite
— Nuisance plein de sinistre comme : l’inondation, les maladies infectieuses
— Ecologique etc. . .
— Milieu naturel
Transition :
Après avoir définis l’assainissement il Y a deux types d’assainissements
Les différents types — L’assainissement dans l’ensemble et l’assainissement
d’assainissement dans la famille
Durée : 15 mn
je fais un exemple sur une famille et sur un ensemble de famille.
- L’assainissement in-
dividuel
- L’assainissement col-
lectif
3.9.2 Document élèves

L’assainissement
I.Définition :
L’assainissement est l’ensemble des opérations permettant la collecte, l’évacuation des
eaux (pluviale, usées) ; ainsi que leur traitement et rejet dans le milieu récepteur.
a. Eau pluviale ou Eau de ruissellement :
Elle provient de la pluie, des eaux de lavage des voix publics, des eaux de drainages.
b. Eaux usées :
Elles proviennent des toilettes (eaux vannes), dans les cuisines (eaux ménagères).
c. Eaux usées industrielles :
Elles proviennent des diverses usines de fabrication et de transformation.
II. But de l’assainissement :
L’assainissement a pour but de :
Sauvegarder l’équilibre écologique du milieu naturel en particulier les eaux souter-
raines et les eaux de surface ;
D’éliminer la nuisance ;
Préserver la santé de la population.
III. Les différents types de l’assainissement :

a. L’assainissement individuel :
L’assainissement individuel consiste à éliminer les eaux usées et les excrétas à travers les
ouvrages d’assainissement qui sont le regard ou tampon de visite, la fosse septique, le
filtre bactérien, le puisard ou puit perdu, le plateau absorbant, les tranchés filtrant et la
latrine VIP.
b. L’assainissement collectif :
L’assainissement collectif consiste à évacuer et drainer les eaux usées et les excrétas à
travers un réseau appelé réseau de d’égout qui les acheminent vers la station d’épuration
avant le rejet dans le milieu naturel.
Schéma
1 Agglomération
2 Regard
3 Branchement
4 Réseau d’égout
5 Station d’épuration
6 Milieu naturel
3.9.3 Evaluation
1 DEFINIR l’assainissement ?
2 Citer les différents types de l’assainissement et expliquer les ?
3.9.4 Réponses
1 L’assainissement est l’ensemble des opérations permettant la collecte, l’évacuation
des eaux (pluviale, usées) ; ainsi que leur traitement et rejet dans le milieu récepteur.
2 Les différents types d’assainissement sont :
- L’assainissement individuel
- L’assainissement collectif
L’assainissement individuel :
L’assainissement individuel consiste à éliminer les eaux usées et les excrétas à travers les
ouvrages d’assainissement qui sont le regard ou tampon de visite, la fosse septique, le
filtre bactérien, le puisard ou puit perdu, le plateau absorbant, les tranchés filtrant et la
latrine VIP.
L’assainissement collectif :
L’assainissement collectif consiste à évacuer et drainer les eaux usées et les excrétas à

travers un réseau appelé réseau de d’égout qui les acheminent vers la station d’épuration
avant le rejet dans le milieu naturel.
3.10 Conclusion
Ce stage pédagogique a été un grand apport pour moi, car il m’a permis de mettre en
pratique les théories que j’ai apprises à l’ENETP, de mieux comprendre et d’acquérir des
nouvelles connaissances sur le terrain.
Durant ce stage j’ai pu mieux comprendre et améliorer plusieurs notions à savoir :
La procédure de préparation des cours ;
L’esprit d’équipe ;
La méthode de recherche ;
Le sens de responsabilité.
Je reconnais que ce stage m’a formé psychologiquement pour faire face aux différentes
tâches qui m’attendent dans ma vie professionnelle (futur). Force est de savoir que J’ai pu
progresser en termes de la responsabilité, de la ponctualité, de l’assiduité, de la discrétion
et de la courtoisie. Pour mettre fin à cette partie j’attirer l’attention des autorités admi-
nistratives et politiques du pays, qu’ils œuvrent en étroite collaboration afin de mettre
en application les résolutions d’accord avec les enseignants. Si rien n’est fait dans un
proche avenir dans l’amélioration des conditions de travail des enseignants, aussi d’avoir
une bonne politique en matière de système éducatif de l’enseignement fondamental, toute
une génération serait sacrifiée. Cela va affecter davantage l’administration publique en ce
qui concerne les ressources humaines qualifiées.

Conclusion générale
Ce document est composé d’un dossier technique et d’un pédagogique.
Nous avons préparé dans un premier temps le dossier technique qui portait sur la
résistance des structures à l’impact d’un projectile, en évaluant le comportement mécanique
du matériau avec le modèle de Mazars.
Dans la première partie du dossier technique nous avons effectué une étude bibliogra-
phique sur le projectile, sur la classification des impacts , sur les modèles de calcule en
dynamique rapide, sur les modèles de comportement des structures sous impact, sur la
force du projectile et sur le calcul de la résistance des structures aux impacts de projectile.
Dans la seconde partie du dossier technique nous avons effectué une simulation numérique
avec le logiciel COMSOL en utilisant le modèle de Mazars.
cette partie du dossier technique est constitué de deux parties dont nous avons la mise en
œuvre des simulations et l’analyse des résultats.
— La mise en œuvre consiste à déterminer les caractéristiques du projectile, la paroi

et les principaux paramètres de calcule du projectile.
— L’analyse des résultats consiste à donner les effets du projectile sur la paroi.
le projectile est un projectile cylindrique de masse 100kg avec une vitesse de 8.43m/s et
d’après certains auteurs on a pris le temps T = 20ms
l’ouvrage est un élément en béton armé dont la classe du béton choisit nous permet de
définir la résistance en compression et traction.
Dans l’analyse des résultats il y a la courbe de compression et de traction, les différentes
graphiques de déformations et les fissures et le maillage de calcul de la structure.
Le dossier pédagogique présente un rapport du stage pédagogique effectué au Centre

de formation en administration et aux métiers et une leçon modèle sur l’assainissement.
Le stage pédagogique m’a permis d’approfondir mes connaissances sur les préparations
de leçons, d’appliquer les théories apprises à l’École Normale d’Enseignement Technique
et Professionnel.
Ce stage a été pour moi une période d’apprentissage, d’écoute, de partage, de patience,
de coopération et de recherche.
55
Bibliographie
[1] Eibl J. – Soft and hard impact, Proc. FIP Congress, The Concrete Society, Concrete
for hazard protection, Edinburgh, Scotland, sept. 1987, pp. 175186.
[2] CEB-88, Concrete structures under impact and impulsive loading, 1988
[3] Brossard, J., 1997. Mécanique générale - Dynamique : théorie classique du choc, s.l. :
Techniques de l’ingénieur, Dossier A1668.
[4] Sugano, T. et al., 1993. Full-scale aircraft impact test for evaluation of impact force,
s.l. : Nuclear Engineering and Design, vol. 140, n°3, p. 373-385. œchlinPotapovœ-
chlinPotapov
[5] Kœchlin, P. et Potapov, S., 2009. Classification of soft and hard impacts-Application
to aircraft crash, s.l. : Nuclear Engineering and Design 239.
[6] Cundall , P. et Strack , O., 1979. A discrete numerical model for gra nular assemblies,
s.l. : Geotechnique, 29 :47–65.
[7] D’addetta, G., Kun, F., Ramm, E. et Herrmann, H., 2001. From solids to granulates-
discrete element simulations of fracture and fragmentation processes in geomaterials,
s.l. : Lecture Note in Physics, 568 :231–258.
[8] Mazars, J., 1984. Application de la mécanique de l’endommagement au comportement

non linéaire et à la rupture du béton de structure, s.l. : PhD thesis, Doctorat d’état
de l’Université Paris VI, Paris, France, 1984
[9] Pontiroli, C., Rouquand, A. et Mazars, J., 2010. Predicting concrete behavior from
quasi-static loading to hypervelocity impact. An overview of the PRM Model., s.l. :
European Journal of Environmental and Civil Engineering, 2010, vol. 14/6-7, pp.
703727 – doi :10.3166/ejece.14.703-727. œstynalœstynal
[10] Wœstyn, S.et al, 2006. Analyse de la rupture dynamique du béton par un modèle
discret., s.l. : Revue Européenne de Génie Civil, 10(10) :1281– 1308. 51.
[11] Ortiz, M. et Pandolfi, A., 1999. Finite-deformation irreversible cohesive elements

for three-dimensional crack-propagation analysis., s.l. : Int. J. Num. Meth. Eng.,
44 :1267– 1282..
[12] Daux, C. et al., 2000. Arbitrary branched and intersecting cracks with the extended
nite element method, s.l. : International Journal for Numerical Methods in Enginee-
ring, 48 :1741–1760.
56
BIBLIOGRAPHIE
[13] Agardh, L. et Laine, L., 1999. 3D FE-simulation of high-velocity fragment perforation

of reinforced concrete slabs, s.l. : International Journal of Impact Engineering 22,
911-922.
[14] Merabet, O., 1990. Modelisation des structures planes en beton arme sous
chargements monotone et cylique, s.l. : Thèse de doctorat. INSA de Lyon.
émondetalémondetal
[15] Frémond, M.et al, 1995. Damage in concrete : the unilateral phenomenon, s.l. : Nu-
clear Engineering and Design 156, 323-335.
[16] Mazars, J., Berthaud, Y. et Ramtani, S., 1990. The unilateral damage behaviour of
damaged concrete, s.l. : Engineering Fracture Mechanics, 35(4) :629–635. 44.
[17] La Borderie, C., 1991. Phénomènes unilatéraux dans un matériau endommageable :

modélisation et application a l’analyse de structure en béton, s.l. : thèse de doctorat
de l’université PARIS 6, 1991..
[18] Cervera, M., Olivier, J. et Manzoli, O., 1996. A rate-dependent isotropic damage
model for the seismic analysis of concrete dams, s.l. : Earthquake engineering and
structural dynamics, vol 25, 987-1010.
[19] Voyiadjis, et Al-un modèle couplé proposé en 2008
[20] Burlion, N. 1998. Conpaction des Bétons : élélemnts de modélisation et ca-

ractérisation expérimentale. PhD thesis, E.N.S. Cachan.
[21] Dragon, A., Mroz, Z. (1979), “A continuum model for plastic-brittle behaviour of
rock and concrete”, Int. J. Engng. Sci., 17, 121-137.
[22] Jason, L.,Huerta, A.,Pijaudier-Cabot, G. etGhavamian, S. (2006). An elastic plastic

damage formulation for concrete : Application to elementary tests and comparison
with an isotropic damage model. Computer Methods in Applied Mechanics and En-
gineering, 195(52) :7077 – 7092. 13
[23] K. WILLAM et E. P. WARNKE : Constitutive model for the triaxial behavior of

concrete. In IABSE Seminar on Concrete Stractures subjected to triaxial sresses,
volume III-1, 1974.
[24] A. L. GURSON : Continuum theory of ductile rupture by void nucleation and growth :
Part i yield criteria and ow rules for porous ductile media. Engineering Materials and
Technology, 99 :2–15, 1977.
[25] La norme NF EN 1991-1-7
[26] Forrestal MJ, Altman BS, Cargile JD, Hanchak SJ. An empirical equation for pe-
netration depth of ogive-nose projectiles into concrete targets. Int J Impact Eng
1994 ;15(4) :395–405.
[27] Li QM, Chen XW. Dimensionless formulae for penetration depth of concrete target
impacted by a nondeformable projectile. Int J Impact Eng 2003 ;28(1) :93–116.

BIBLIOGRAPHIE
[28] The resistance of various metallic materials to perforation by steel fragments ; empiri-
cal relationships for fragment residual velocity and weight, Technical Report Project
THOR TR No.47, 1961 and No.51, 1963, Ballistic Research Laboratory, Aberdeen
Proving Ground, Washington state, USA.
[29] Methods for the determination of possible damage to people and objects resulting
from the releases of hazardous materials, CPR 16 E, (Green Book), Committee for
the prevention of disasters caused by dangerous substances, the Hague : Directorate-
General of Labour of the Ministry of Social Affairs and Employement, 1992 - 2003
[30] Cox, B.G and Saville, G. (Eds). High Pressure Safety Code. Hign Pressure Associa-
tion, 1975. Reprinted 1985
[31] Méthode pour l’identification et la caractérisation des effets Dominos, mise en place
de la directive SEVESO II, 2002.

Chapitre Un

Transféré par

Droits d'auteur :

Formats disponibles

Chapitre Un

Transféré par

Informations du document

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

Chapitre Un

Transféré par

Droits d'auteur :

Formats disponibles

émondetal15

Résistance des structures à l’impacte d’un projectile 1

Liste des symboles 7

1.8.3.2 Impact à vitesse subsonique . . . . . . . . . . . . . . . . . 32

Résistance des structures à l’impacte d’un projectile 3

3.5 Solutions apportées . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 47

Résistance des structures à l’impacte d’un projectile 4

1.1 Paramètre de perforation selon la nature du matériau . . . . . . . . . . . . 33

2.1 classe du béton . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 37

3.1 Récapitulatif des effectifs d’élèves par filières . . . . . . . . . . . . . . . . . 46

1.1 Modèle mécanique simple d’un impact à travers un système de masse-

3.1 Organigramme de l’établissement. . . . . . . . . . . . . . . . . . . . . . . 44

Je dédie ce modeste travail :

Je remercie également toute l’équipe pédagogique du Centre de formation en admi-

J’exprime ma gratitude à l’ensemble des enseignants du département génie civil de

Je souhaite également à remercier le directeur des études de l’ENETP, Dr. Yacouba

Je souhaite également remercier ma femme Gnia Traoré et toutes les membres de sa

Le présent document est subdivisé en deux grandes parties :

La première partie comprend deux (2) chapitres :

La seconde partie comprend :

Résistance des structures à l’impacte d’un projectile 13

Ces projectiles sont capables de nuire à la résistance de la cible visée (structure de

1.2.2 Choix et caractéristique du projectile

1.3 Classification des impacts

Résistance des structures à l’impacte d’un projectile 16

m1 ẍ1 (t) + K1 [x1 (t) − x2 (t)] = 0

m1 ẍ2 (t) + K1 [x1 (t) − x2 (t)] + K2 x2 = 0 (1.1)

Les équations [1-1] sont réécrites en [1-3] :

m1 ẍ1 (t) + K1 x1 (t) = 0(a)

m2 ẍ2 (t) + K2 x2 (t) = F (t)(b) (1.3)

Résistance des structures à l’impacte d’un projectile 17

Figure 1.2 – Impact mou (CEB,1987)

l’énergie cinétique de l’impacteur déformable est complètement transférée en énergie de

Figure 1.3 – impact dur (Eibl,1987)

D’après (Eibl,1987) et (CEB,1988) un autre auteur (Brossard,1997) trouve une autre

Résistance des structures à l’impacte d’un projectile 18

Avec σp : la contrainte de rupture du projectile,

Figure 1.4 – classification des impacts (Kœchlin et Potapov, 2009)

1.4 Impacte du projectile sur une structure en béton

Résistance des structures à l’impacte d’un projectile 19

un enchaı̂nement de processus catastrophiques successifs appelé effets dominos. L’acciden-

1.5 Modèles de calcul en dynamique rapide

— méthodes des éléments de frontière

1.5.1 Méthode des éléments discrets

Résistance des structures à l’impacte d’un projectile 20

Figure 1.5 – Echelle de modélisation (D’addetta,et al.,2001)

1.5.2 Méthode des éléments finis

Résistance des structures à l’impacte d’un projectile 21

maillage adaptatif. Cette méthode introduit une dépendance au maillage utilisé.

1.6 Modèles de comportement des structures sous

1.6.1 Modèle d’endommagement

Résistance des structures à l’impacte d’un projectile 22

distingués en fonction de l’ordre tensoriel des variables d’endommagement utilisées. En

φ̇ = σij ij − ρψ̇ = Y Ḋ ≥ 0 (1.10)

Comme Y ≥ 0 ; ( D ≥ 0 ) : Cela signifie que l’endommagement est toujours crois-

Avec Di = Dt pour la partie d’endommagement en traction et Di =Dc pour la partie

Où αt proportionnel aux déformations positives, correspond à la traction,β est un pa-

Résistance des structures à l’impacte d’un projectile 23

(Pour la partie en compression)

Où σi représente la ième contrainte principale du tenseur σ̄ , Pi N représente le vecteur

φ̇ = σij ij − ρψ̇ = Y Ḋ ≥ 0 (1.10)